井孔模式破频散曲线计算.docVIP

中国石油大学（华东）毕业设计（论文） 井孔模式波频散曲线计算 学生姓名：黑 创 学 号 专业班级：信息与计算科学06-3班 指导教师：申淑谦 苏远大 2010年6月24日 摘 要 本文根据井孔中波的传播理论，建立弹性波动方程，理论推导并求出地层的弹性波动解；利用井孔中波的传播理论，模拟了井孔中各种模式波的频散效应；利用摄动积分对频散曲线进行校正，并可根据弯曲波的频散曲线的交叉判断地应力的存在；建立一个纵向无限延伸并充满液体的井孔模型，利用牛顿迭代算法和有限差分编写程序，计算出了频散曲线。 关键词：弹性波动方程；频散曲线；弯曲波；牛顿迭代；有限差分 ABSTRACT In this paper，according to the theory of the spread of waves of well, we establish elastic wave equation, theoretically analyze and obtain solution of formation of elastic wave ; using of transmission theory of waves of well, we simulate models wave in the borehole frequency dispersion;we use perturbation integral to regulate the dispersion curve, and according to flexural wave dispersion curves crossover check the existence of stress;we establish a vertical infinite fluid-filled borehole and model, using Newton iterative algorithm and the finite difference programming to calculate dispersion curve. Keywords: elastic wave equation; dispersion curve; flexural wave; perturbation integral; finite difference 目 录 第1章 前 言 1 第2章 理想无限介质井孔中的频散方程 4 2.1 基本理论 4 2.1.1 波动方程及其势函数 4 2.1.2 单极子源的波动方程及其势函数 5 2.1.3 单极子源和多极子源的波动方程及其势函数 6 2.1.4 位移和应力 8 2.1.5 边界条件及频散方程 9 2.2 数值计算 12 2.2.1 偶极子波频散方程 12 2.2.2 软、硬地层的频散曲线 13 第3章 地应力对井孔模式波频散曲线的影响分析 16 3.1 基本原理 16 3.2 数值计算与结果分析 17 3.3 结论 20 第4章 基于声弹性理论研究不均衡水平地应力引起的地层横波速度分布规律 21 4.1 基本原理 21 4.2 数值计算与结果分析 21 4.3 结论 24 第5章 结 论 25 致 谢 26 参考文献 27 附 录 32 第1章 前 言 自从本世纪五十年代，声波测井技术应用于地球物理勘探以来[1~2]，声波测井的理论、方法和仪器都得到了很快的发展[3~4]，尤其是七十年代后期和八十年代前期，对点声源在井孔中激励的声场分布，在理论上进行了详尽的分析论证[5~13]。利用数字计算机技术发展成为数字全波列声波测井方法[14~16]，从而获得大量的有关井壁岩体的新参数，使声波测井成为井孔物探中最重要的三大方法之一，是探测岩体弹性特性最为有效的手段。 对全波列声波测井理论计算结果表明，点声源在硬地层（指井壁岩体的横波声速大于井液声速的高速地层）井孔中将激励滑行纵波、滑行横波、伪瑞利波和斯通利波；而在软地层（指井壁岩体的横波声速小井液声速的低地层）则仅能激发滑行纵波和斯通利波。这样，在软地层井孔中或硬地层井孔中的软地层带，无法依靠井壁的折射形成滑行横波，所以在这些层段很难精确测定横波声速，而恰恰是这些层段对工程上具有重要意义。另一方面，全波列测井提取横波参数的数据处理方法十分繁杂，尽管采用了长源距的井下声系，纵波和横波在时域上或多或少地发生重叠，为了克服以上不足，人们寻求了新的声波模式作为检测声波的声波测井[17~21]，即多极子源的横波测井。 1967年美国地球物理学家White J. E. [17]教授首先提出了利用偶极子声源直接激励