专利名称：一种泄漏模式波反演软地层横波速度的方法
技术领域：
本发明属于应用地球物理声波测井技术，特别是涉及超大井径、疏松 浅层或扩径严重的软地层勘探时，利用频散的泄漏模式波来反演地层横波 速度的方法。
背景技术：
随着石油勘探和开发技术的发展，地震勘探和钻井工程等领域均迫切 要求声波测井提供准确的地层横波速度剖面，如地震弹性波阻抗反演、地
震AVO属性处理、钻井地层压力预测和井壁稳定性评价。目前，井孔声波 测井技术是提供连续地层横波速度剖面的最有效手段。但在地层横波速度 小于井内泥浆声速的软地层，传统声波测井的对称单极子声源不能激发起 地层的横波，而是利用单极激发的斯通利波反演地层横波速度。近年来， 发展了非对称的偶极子声源，通过激发井孔弯曲模式波实现测量任意地层 横波速度的偶极子声波测井技术。实际测井时，单极斯通利波和偶极弯曲 波的激发强度和频散特性均受地层性质和井径等参数影响严重，甚至在大 井径、疏松的浅层或扩径较严重的软地层，偶极声源往往不能激发起弯曲 波，单极声源也不能激发起斯通利波，如图1中430 435米和446 452 米井段所示。因此，在大井径、疏松的浅层或扩径较严重的软地层，从现 有偶极声波测井技术测量的角度和单极斯通利波反演的角度，均无法保证 提供地层的横波速度。
从大量现场多极子(单极和偶极)阵列声波测井资料中发现，大井径、 疏松的浅层或扩径较严重的软地层偶极测井无法直接测量地层的横波速 度。建立井孔多极子声波测井理论模型，数值模拟研究软地层的单极、偶 极声波测井响应特征表明，在理论上扩径严重的软地层单极、偶极均激发 幅度较大的泄漏模式波，而用于测量横波速度的单极斯通利波和偶极弯曲 波幅度很小，甚至完全消失，如图2和图3所示。软地层单极、偶极激发 的泄漏模式波是一种频散模式波，即相速度随频率变化而变化，如图4所示。泄漏模式波的频散特性主要敏感于地层纵波速度、井内泥浆声速和地
层横波速度，如图5所示，其中，该模式波对地层横波速度的敏感程度达 到10%。因此，在单极斯通利波和偶极弯曲波均没有激发起来的软地层情 况，可考虑用相对比较稳定的单极或偶极的泄漏模式波的频散特性反演软 地层的横波速度。
发明内容
本发明的目的在于提供一种泄漏模式波反演软地层横波速度的方法, 以解决利用频散的泄漏模式波对大井径、疏松浅层或扩径较严重的软地层 的横波速度进行反演的技术问题。
本发明的目的是基于内存256M以上的微型计算机的Windows环境， 利用Microsoft ￥0++软件开发工具，通过如下方式实现
首先，将扩径严重的软地层井段一个深度点的由多极子阵列声波测井 仪器单极或偶极工作模式测量的泄漏模式波的波形数据和对应的该深度点 由常规测井测量的地层纵波速度、井径、泥浆声速、泥浆密度参数输入给 计算机的几个存储器。
其次，利用频谱加权相关法、Prony方法、修正的Matrix Pencil法或 同态处理方法的阵列波形频散曲线提取技术，由软地层的单极子阵列波形 或偶极子阵列波形提取出泄漏模式波的相速度频散曲线，并确定单极泄漏 模式波或偶极泄漏模式波的主要能量频段。
再次，根据Kurkjian等人提出的充液井孔声场理论，建立井孔多极子 声波测井模型，模型中的地层纵波速度、密度、井径、泥浆声速和泥浆密 度由传统的测井方法测量已知，改变地层横波速度用迭代法或抛物线法由 理论的频散方程"("H^"^)^解出软地层单极或偶极泄漏模式 波的相速度随频率变化的频散曲线，式中，"2^^为波数，/一频率，^ 一相速度，^一地层横波速度，F—地层纵波速度，^一井内泥浆声速，P 一地层密度，^一泥浆密度，及一井眼半径。
然后，在单极或偶极泄漏模式波主要能量频段内，计算不同i也层横波 速度下理论计算的泄漏模式波相速度与现场阵列波形资料提取的泄漏模式 波相速度的方差之和，得到随地层横波速度变化的方差曲线。利用拉格朗 日法对方差曲线精确插值，确定出方差为零时对应的地层横波速度，即为 该深度点反演得到的地层横波速度。最后，利用上述反演流程对扩径严重的软地层井段进行连续处理，获 得软地层井段的地层横波速度成果图，将成果图显示在显示器上，或打印 出图。
本发明克服了现有多极子阵列声波测井仪器在扩径严重的疏松软地层 井段不能直接测量地层横波速度的缺点，利用单极或偶极阵列波形中频散 的泄漏模式波来反演地层横波速度，这一方法拓宽了井孔声波测井技术在 石油勘探与开发过程中的应用领域，提高了其应用价值。
图1是疏松浅层井段单极和偶极声波测井的波形资料及对比示意图。
图2是软地层理论合成的单极子声波测井的阵列时域波形示意图。图 3是软地层理论合成的偶极子声波测井的阵列时域波形示意图。 图4是软地层单极子声源激发的泄漏模式波的相速度频散曲线图。 图5是软地层单极子泄漏模式波对地层横波速度、纵波速度和井内泥
浆声速的敏感性曲线图。
图6是依据本发明所提供的一种频散的泄漏模式波反演地层横波速度
方法的流程框图。
图7是理论合成的软地层的单极子阵列时域波形示意图。
图8是软地层单极子阵列波形对应的频谱示意图。
图9是利用频谱加权相关法从单极子阵列波形得到的泄漏模式波的频
谱加权相干图和提取出的频散曲线(粗实线所示)。 图10是井孔单极子声波测井理论模型示意图。
图11是依据本发明所述方法从阵列波形中提取的泄漏模式波相速度 频散曲线(小圆圈)和调整地层横波速度理论计算的泄漏模式波相速度频散 曲线(实线)图。
图12是依据本发明所述方法在单极泄漏模式波主要能量频段内，不同 地层横波速度下计算的泄漏模式波频散曲线与提取频散曲线的方差之和， 由两者方差等于零时反演出软地层横波速度的示意图。
图13是依据本发明所述方法对现场浅层井段(地层纵波速度在2500米 /秒左右，地层横波速度在1200左右的)多极子阵列声波测井资料的单极泄 漏模式波反演得出的地层横波速度成果图。
具体实施方式
本发明关于附图的说明如下
图1反映了疏松浅层井段单极和偶极声波测井的波形资料及对比，该 井段为软地层井段，图中第l道给出自然伽马、钻头尺寸和井径曲线，第
2、 3道分别给出偶极和单极波形。从430 435米和446 452米两个井段 的单极和偶极波形可看出由于扩径严重导致单极和偶极声波测井资料只 有一种泄漏模式波，没有激发起单极斯通利波和偶极弯曲波。
图2反映了软地层理论合成的单极子声波测井的阵列时域波形，从图 中可看出，该软地层情况单极声源只激发起泄漏模式波，没有斯通利波。
图3反映了软地层理论合成的偶极子声波测井的阵列时域波形，从图 中可看出，该软地层情况偶极声源只激发起泄漏模式波，没有弯曲波。
图4反映了软地层单极子声源激发的泄漏模式波的相速度频散曲线， 其中横坐标为频率，纵坐标为相速度。
图5反映软地层单极子泄漏模式波对地层横波速度、纵波速度和井内 泥浆声速的敏感性曲线，从图中可看出泄漏模式波的频散特性在低频段主 要敏感于地层纵波速度，高频段主要敏感于井内泥浆声速。同时发现在主 要能量频带该模式波的频散特性对地层横波速度也比较敏感，敏感系数达 到10%左右。
图7反映了理论合成的软地层的单极子阵列时域波形，其中只有一种 泄漏模式波，没有斯通利波。
图8反映了软地层单极子阵列波形对应的频谱，主要能量集中在5千 赫兹 10千赫兹。
下面结合附图和实施例来详细描述本发明。
在扩径严重的软地层，单极声源和偶极声源均激发起幅度很强的频散 的泄漏模式纵波，该模式波相对比较稳定，且非常敏感于地层横波速度。 因此，本发明在井径大于30厘米、地层横波速度小于1400米/秒的软地层， 可采用单极泄漏模式波反演地层横波速度和偶极泄漏模式波反演地层横波 速度两种具体实施方案。
方案一，单极泄漏模式波反演地层横波速度的具体步骤为
(1).将扩径严重的软地层井段一个深度点的由多极子阵列声波测井仪 器单极工作模式测量的泄漏模式波的波形(如图7所示)数据和对应的该 深度点由常规测井测量的地层纵波速度、井径、泥浆声速、泥浆密度参数输入给计算机的几个存储器。
(2).利用快速傅立叶变换对如图7所示的单极子阵列波形做频谱分析， 获得每道波形的频谱，如图8所示。在此基础上，利用下述公式由阵列波 形获得频谱加权相关图，并由相关系数极大值提取出单极泄漏模式波的相 速度频散曲线，如图9所示，<formula>formula see original document page 7</formula>
式中/为频率，"2^/^为波数，Kp为相速度；W为接收器的个数，d是 相邻两个接收器的间距；X"(/)是第"个接收器的波形对应的频谱，"*"表 示取复共轭。
(3) .建立井孔单极子声波测井理论模型，如图10所示，根据井孔内、 外波动方程和边界条件推导出井孔单极子模式波的理论频散方程 Z)(/,A:，^，)^,^，Ap,，i )-0，其中，A^2; /VFp为波数，/一频率，^一相速度， K一地层横波速度，^一地层纵波速度，r,一井内泥浆声速，p—地层密 度，p,—泥浆密度，i 一井眼半径。模型中地层纵波速度、密度、井径、 泥漿声速和泥浆密度由其他测井方法测量已知，改变地层横波速度用牛顿 迭代法由频散方程求解单极泄漏模式波的相速度频散曲线，如图11中实线 所示。
(4) .在单极泄漏模式波主要能量频段内，即5千赫兹 10千赫兹，计 算不同地层横波速度下理论计算的泄漏模式波相速度(图11中的实线)与阵 列波形资料提取的泄漏模式波t目速度(图11中的小圆圈)的方差之和，得到 随地层橫波速度变化的方差曲线，如图12所示。
(5) .通过对方差曲线插值，确定方差为零时对应的地层横波速度，即 为该深度点反演得到的地层横波速度，如图12所示。按本发明所述方法得 到的该地层横波速度为507.32m"，实际的地层横波速度为508m"。
方案二，偶极泄漏模式波反演地层横波速度的具体步骤为 (1).将扩径严重的软地层井段一个深度点的由多极子阵列声波领拼仪 器偶极工作模式测量的泄漏模式波的波形数据和对应的该深度点由常规测 井测量的地层纵波速度、井径、泥浆声速、泥浆密度参数输入给计算机的 几个存储器。(2) .禾拥快速傅立叶变换对偶极子阵列波形做频谱分析，获得每道波 形的频谱。在此基础上，利用频谱加权相关法提取出偶极泄漏模式波的相 速度频散曲线。
(3) .建立井孔偶极子声波测井理论模型，根据井孔内、外波动方禾呈和 边界条件推导出井孔偶极子模式波的理论频散方程。模型中地层纵波速度、 密度、井径、泥浆声速和泥浆密度由其他测井方法测量已知，改变地层横 波速度用牛顿迭代法由频散方程求解偶极泄漏模式波的相速度频散曲线。
(4) .在偶极泄漏模式波主要能量频段内，即2千赫兹 5千赫兹，计 算不同地层横波速度下理论计算的泄漏模式波相速度与现场偶极阵列波形 资料提取的泄漏模式波相速度的方差之和，得到随地层横波速度变化的方 差曲线。
(5) .通过对方差曲线插值，确定方差为零时对应的地层横波速度，即 为该深度点反演得到的地层横波速度。
综上所述，用本发明所述方法处理浅层井段的多极子声波测井资料， 输出连续的地层横波速度剖面，如图13所示。该浅层井段的地层纵波速度 在2500米/秒左右，地层横波速度在1200左右，属于典型的软地层。如图 中第1道钻头尺寸、井径和自然伽马曲线所示，该井使用9.5英寸钻头打 井在430 435米和446 453米井段出现严重扩径现象，井径达到50厘米。 图中第2和3道分别给出多极子阵列声波测井资料的单极和偶极的全波列 波形，从全波列波形可看出，在扩径严重的两个井段单极和偶极均激发幅 度很强的泄漏模式波，而没有激发用于软地层测量横波速度的单极斯通利 波和偶极弯曲波。用传统的慢度时间相关法处理偶极阵列波形资料，得到 相关系数平面和对应的地层横波速度如图中第4道所示。从相关系数平面 图看出，只有在435 445米规则小井径井段偶极弯曲波的相关系数较大， 可由偶极弯曲波提取出可靠的地层横波速度，而在其他两个扩径严重的井 段偶极弯曲波的相关性很差甚至没有，此时无法由偶极资料提取出地层横 波速度。图中第5道给出利用本发明方法利用单极泄漏模式波频散特性反 演的地层横波速度与规则井段偶极提取的横波速度、地层纵波速度对比。 对比结果表明，本发明方法在井况复杂的疏松浅层可由泄漏模式波反演获 得连续的地层横波速度剖面。在435 445米井径规则的软地层井段，偶极 弯曲波波形较好，本发明所述方法的反演结果(粗实线)与偶极提取的横波速度结果(小圆圈)吻合的很好；在430 435米和446 453米扩径严重井 段，没有激发起偶极弯曲波，从偶极资料无法提取地层横波速度，本发明 所述方法由频散的单极泄漏模式波反演出的地层横波速度稳定、可靠，与 地层纵波速度趋势相一致且纵横波时差比合理，可用于石油勘探和开发的 工程应用。
1、 频谱加权相关法参考文献唐晓明，郑传汉著，2004年，定量测 井声学，石油工业出版社，p 52-54。
2、 Prony方法参考文献Lang， S. W. Kurkjian， J. H.等，1987， Estimating slowness dispersion from arrays of sonic logging waveforms, Geophysics, Vol. 52， p 530-544。
3、 Matrix Pencil法参考文献Ekstrom， M. P" 1995，Dispersion estimation from borehole acoustic arrays using a modified matrix pencil algorithm, 29th Asilomar Conference on Signals Systems and Computers, p 449-453 。
4、 同态处理方法参考文献Ellefsen, K. J. Cheng, C. H. Tubman, K. M. ,1989， Estimating phase velocity and attenuation of guided waves in acoustic logging data, Geophysics, VOL, 54,p 1054-1059。
5、 井孔声场建模理论参考文献Kurkjian, J. H. and Chang，S.K.， 1986， Acoustic multipole sources in fuid-filled boreholes, Geophysics.51 (1) 148-163。
权利要求
1、 一种泄漏模式波反演软地层横波速度的方法，它是基于内存256M以上的微型计算机的Windows环境，利用Microsoft VC++软件开发工具，通过如下方式实现首先，将扩径严重的软地层井段一个深度点的由多极子阵列声波测井仪器单极或偶极工作模式测量的泄漏模式波的波形数据和对应的该深度点由常规测井测量的地层纵波速度、井径、泥浆声速、泥浆密度参数输入给计算机的几个存储器；其次，利用频谱加权相关法、Prony方法、修正的Matrix Pencil法或同态处理方法的阵列波形频散曲线提取技术，由软地层的单极子阵列波形或偶极子阵列波形提取出泄漏模式波的相速度频散曲线，并确定单极泄漏模式波或偶极泄漏模式波的主要能量频段；再次，根据Kurkjian等人提出的充液井孔声场理论，建立井孔多极子声波测井模型，模型中的地层纵波速度、密度、井径、泥浆声速和泥浆密度由传统的测井方法测量已知，改变地层横波速度用迭代法或抛物线法由理论的频散方程D(f，k，Vs，Vc，Vf，ρ，ρf，R)＝0解出软地层单极或偶极泄漏模式波的相速度随频率变化的频散曲线，式中，k＝2πf/Vp为波数，f-频率，Vp-相速度，Vs-地层横波速度，Vc-地层纵波速度，Vf-井内泥浆声速，ρ-地层密度，ρf-泥浆密度，R-井眼半径；然后，在单极或偶极泄漏模式波主要能量频段内，计算不同地层横波速度下理论计算的泄漏模式波相速度与现场阵列波形资料提取的泄漏模式波相速度的方差之和，得到随地层横波速度变化的方差曲线；利用拉格朗日法对方差曲线精确插值，确定出方差为零时对应的地层横波速度，即为该深度点反演得到的地层横波速度；最后，利用上述反演流程对扩径严重的软地层井段进行连续处理，获得软地层井段的地层横波速度成果图，将成果图显示在显示器上，或打印出图。
全文摘要
本发明提供了一种泄漏模式波反演软地层横波速度的方法，属于应用地球物理声波测井技术。具体步骤如下输入扩径严重的软地层井段的单极或偶极泄漏模式波波形和对应井段地层的常规测井资料，再由阵列波形中提取出泄漏模式波的相速度频散曲线，并确定泄漏模式波的主要能量频段，计算不同地层横波速度下的理论计算的与现场阵列波形提取的泄漏模式波相速度的方差，确定方差为零时对应的地层横波速度，作连续井段处理并输出反演横波速度成果图。本发明解决了现有多极子阵列声波测井仪器，在超大井径、疏松浅层或扩径严重的软地层井段不能直接测量地层横波速度的问题，拓宽了井孔声波测井技术在石油勘探与开发过程中的应用领域和提高了其应用价值。
文档编号E21B49/00GK101285381SQ20071009038
公开日2008年10月15日 申请日期2007年4月9日 优先权日2007年4月9日
发明者孙建孟, 庄春喜, 汤丽娜, 秦瑞宝, 苏远大, 邓启才, 陈雪莲 申请人:中国石油大学(华东);中海石油研究中心