多极子阵列声波测井数据的纵横波慢度联合反演方法

多极子阵列声波测井仪器采集的单极和偶极数据受到地层、井孔、仪器测量系统的影响.在处理实际声波测井数据时,必须考虑多极子模式波的频散效应,以及测井仪器在其中的影响.根据仪器等效理论和相位匹配方法,本文提出了一种从多极子阵列声波测井数据中同时获得纵、横波慢度的联合反演方法.这种方法的关键在于利用相同仪器-地层模型计算多极子模式波频散曲线,以此来匹配频域内纵波与横波数据的相位.相对于将泄漏纵波和弯曲波频散效应分开处理的其他方法,该方法不仅可以减少纵横波速度反演的不确定性,而且还避免了从声波数据中提取频散数据的繁琐过程.通过理论分析和现场数据处理证明了本文联合反演方法的准确性和有效性.

     

随钻斯通利波测井反演地层渗透率的理论、方法及应用

针对孔隙渗透地层的随钻声波测井问题,用Biot-Rosenbaum孔隙弹性波测井理论推导了孔隙地层的随钻井孔声场表达式.据此考察了随钻条件下井中斯通利波的波形、相速度频散、衰减以及相速度对渗透率的灵敏度,并与电缆测井中的情况进行了对比.数值模拟结果表明,随钻条件下斯通利波对地层渗透率的灵敏度相对于电缆测井有明显增加,更有利于用来反演地层渗透率.为快捷有效地处理现场测井数据和反演计算,采用简化Biot-Rosenbaum理论和钻铤的等效模型,对随钻斯通利波的频移和时滞进行联合反演.结果表明,随钻斯通利波反演的渗透率与核磁渗透率和岩心覆压测试渗透率符合较好,并且与常规测井曲线所反映的储层性质具有较好的一致性,证明了利用随钻斯通利波评价地层渗透率的有效性.

     

瑞利型槽波品质因子及特性分析(英文)

为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子Q_R反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子Q_R曲线;采用控制变量法,分析了瑞利型槽波品质因子Q_R随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子Q_R曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子Q_R最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子Q_R随煤层横波品质因子Q_(S2)的增大而增大。利用瑞利型槽波品质因子Q_R曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。     

煤层厚度变化时地震槽波理论频散曲线计算方法及频散特征分析

槽波地震勘探利用槽波的频散特性反演煤层的结构特征,故理论频散曲线的计算是一个重要方面.使用水平层状模型假设下的面波频散曲线计算方法能够计算煤层厚度恒定模型地震槽波频散曲线;但当煤层厚度变化时该方法不再适用.基于前人水平层状均匀介质模型的面波理论频散曲线计算方法,对于含煤三层模型,本文发展了煤层厚度变化情况下的地震槽波理论频散曲线计算方法,并使用该方法计算分析了不同厚度函数模型的频散曲线形态特征.研究表明：与稳定厚度煤层相比,煤层厚度变化使得地震槽波群速度成为与频率及传播射线在水平面投影路径相关的二元函数;射线路径上煤层厚度的变化使得频散曲线在群速度方向上压缩,群速度变化范围变小,且使处于最小值位置的埃里相群速度增大;而煤层厚度的线性变化模型频散曲线只与射线首、尾处的煤层厚度有关,与煤层厚度恒定模型相比,曲线形态不发生改变;煤层厚度呈非线性变化时,频散曲线形态上可能发生改变.

     

基于八阶NAD算子的保辛分部Runge-Kutta方法及其波场模拟(英文)

基于声波方程扩充的哈密尔顿系统,本文给出了空间精度为八阶的近似解析离散化(NAD)保辛分部Runge-Kutta方法,简称八阶NSPRK方法。该方法采用八阶精度的近似解析离散算子近似空间高阶偏微分算子,并使用二阶精度的辛分部Runge-Kutta方法进行时间离散。我们从理论和数值计算两个方面研究了八阶NSPRK方法的稳定性条件和数值频散关系,并同四阶NSPRK方法、八阶Lax-Wendroff(LWC)方法和八阶交错网格(SG)方法进行了比较。结果表明八阶NSPRK方法压制数值频散的能力显著优于传统数值计算方法。与四阶NSPRK方法和传统四阶辛格式(SPRK)方法相比,八阶NSPRK方法具有最小的数值误差和最高的计算效率:在达到同样消除数值频散的前提下,八阶NSPRK方法的计算速度约为四阶NSPRK方法的2.5倍、为四阶SPRK方法的3.4倍;八阶NSPRK方法的存储量仅为四阶NSPRK方法的47.17%、为四阶SPRK方法的49.41%。在双层介质、非均匀介质和Marmousi等复杂速度模型中,八阶NSPRK方法模拟得到的波场快照非常清晰,无可见数值频散。这些结果表明,八阶NSPRK方法在粗网格条件下能有效地压制数值频散,从而能够极大地节省计算内存,提高计算速度。总体而言,八阶NSPRK方法是一种在地震探测领域和地震学研究中有着巨大应用潜力的数值计算方法。     

声波测井中的相速度与群速度讨论(英文)

声波测井过程中获取的速度到底是相速度还是群速度,目前仍存在一些争议,本文从理论分析和数值模拟的角度,使用三种模型对这一问题进行了研究。首先,构造一个相速度与群速度可调的稳态声波传播模型——不同声速的两个平面波叠加模型,利用慢度时间相关(STC)方法提取声波波速,数值模拟结果表明,无论相速度较大或是群速度较大,STC方法提取出来的波速都是相速度;其次,通过频散分析和割线积分得到刚性壁圆柱流体模型中的频散曲线与分波波形,使用STC方法得到的速度与相速度的频散曲线吻合较好,而直接读取波至获得的速度与群速度的频散曲线趋势一致;最后,利用频散分析和实轴积分方法,获得偶极子在慢地层中激发的模式及全波波形,得到的结果再次验证了刚性壁圆柱流体模型中的结论。     

背景噪声频散曲线测定及其在华北地区的应用

利用华北地震台阵垂直分量的观测资料,采用滑动绝对平均方法对资料进行预处理,通过互相关方法从背景噪声中提取瑞利面波的格林函数,开发了群速度频散曲线的自动提取工具,测量了位于华北盆地、燕山隆起和太行山隆起的3条频散曲线,反演得到了3个区域的S波速度结构．研究分析表明,滑动绝对平均方法可以有效降低地震和台站附近干扰源的影响．为了得到可靠的层析成像结果,应计算格林函数的信噪比,选择高信噪比的格林函数测量其频散曲线,进行层析成像反演．当信噪比大于7时,一般都能得到稳定可靠的频散曲线．群速度频散曲线的最大可信周期(Tmax)与台站间距有关,华北地区最大可信周期以不超过台站间距的1/12为宜,周期大于Tmax时不同月份测得的频散曲线变化较大．     

等效介质理论中裂缝密度的适用性研究

为了比较裂缝等效介质理论的预测精度,本文利用物理模拟的方法制作了一组9个裂隙密度模型,其中裂隙密度从0增大到12%,裂隙厚度和直径保持不变,分别为0.12 mm和3 mm.使用0.5 MHz的纵横波换能器,采用脉冲透射法得到纵横波的速度和各向异性.然后,将基质的模量、裂缝填充物模量和裂缝参数代入等效介质理论中,计算出纵横波速度和各向异性随着裂隙密度的变化情况.通过纵横波速度的实测值与等效介质理论的定量对比,发现当裂隙密度小于12%时,Hudson理论和线性滑动理论对纵波速度的预测与实测值非常吻合.当裂隙密度小于8%时,Hudson理论和线性滑动理论对快横波速度的预测与实测值非常吻合.当裂隙密度小于9%时,各向异性自洽理论和NIA（非相互作用近似）理论对慢横波速度的预测与实测值非常吻合.通过纵横波速度各向异性实测值与理论计算值定量对比,发现当裂隙密度小于12%时,Hudson理论和线性滑动理论对纵波速度各向异性的预测与实测值更吻合;当裂隙密度小于6%时,各向异性自洽理论和NIA理论对横波速度各向异性的预测与实测值更吻合.

     

我国大陆表面沉积层的研究

我国西部地震常常激发较强的短周期瑞雷和乐夫型面波,具有明显的波散。这些波主要在沉积层内传播。根据频散曲线确定了沉积层的厚度:青藏高原约为5—6公里;塔里木盆地为9—10公里;中部为5公里;北部平均为6—8公里。     

改进的等效半空间法及瑞雷波频散曲线反演

层状模型的理论频散曲线表明,层状模型表面瑞雷波同一频率对应的相速度具有多值特点,这就是所谓的瑞雷波多阶性. 面波的多阶性导致了实测频散曲线的复杂性,当地层中存在软弱夹层时这一问题尤为突出,从而给实测频散曲线的解释带来了困难. 本文从工程实用的角度出发,根据等效半空间理论,提出了一种计算理论频散曲线的新算法——改进的等效半空间法,避开了面波多阶性这一复杂问题. 由此方法计算得到的面波相速度并不同于某个具体阶数的面波,而是对多阶面波的综合反映. 按照改进的等效半空间法编制了拟合反演程序,并在工程中进行了运用,取得较好的效果.     

基于超级混合虚折射干涉法的阵列声波处理方法及其在套管井中的应用

在声学测井勘探时,复杂多变的测井条件（如胶结不好的套管井和随钻测井等）下的强直达干扰信号会导致地层信号的信噪比降低,使得地层有用信号的提取十分困难.本文在超级虚折射法的基础上,利用声源与接收器的互易原理,对原始波场信号进行互相关和褶积重建,并对所有可能的位形产生的信号进行叠加,将该叠加波场与超级虚折射波场进行叠加,重构出超级混合虚折射波场,使不同接收器上的折射波信噪比在超级虚折射的基础上进一步得到提高,显著提高了有用信号的信噪比,同时解决了超级虚折射近场能量不足的问题.套管井测井的理论模拟及实际处理表明,该方法可以很好地压制套管波信号、提高地层信号的信噪比.

     

基于波干涉原理的套管井阵列声波处理方法及其应用

在套管井测井中,沿着钢管产生的拉伸波(即管波)在声波数据中产生很强的干扰,使得地层声速的提取十分困难.基于管波与地层声波的干涉原理,本文提出了一种声波处理方法.管波的传播特征是已知的,且管波和地层波都由同一声源激发,因此,可以用管波来改造测井声波数据.对改造后数据采用希尔伯特变换作进一步处理,得到与地层信息相关的解析信号,从该信号中即可提取出地层声速.通过对理论模拟数据的处理和现场测井数据的应用,验证了新方法的正确性和有效性.     

Study on 3D simulation of wave fields in acoustic reflection image logging

The borehole acoustic reflection imaging logging is a newly developed acoustic logging method that has attracted many interests. These converted and reflected waves for imaging are usually mixed up with borehole guided waves and therefore difficult to be clearly identified. To improve the downhole tool design and develop more sophisticate data processing and interpretation algorithms,studies on precisely numerical modeling of the wave fields in the acoustic reflection imaging logging are neces-sary and critical. This paper developed a parallelized scheme of 3D finite difference (3DFD) with non-uniform staggered grid and PML absorbing boundary to simulate the acoustic wave fields in isotropic and anisotropic formations. Applications of this scheme to the typical cases of isotropic and anisot-ropic formations and comparison with the results from published analytical solutions have demon-strated the validation and efficiency of the scheme. Higher accuracy and lower computation cost (3.5 times faster than the conventional schemes) have been achieved with this scheme for modeling such a complex wave fields of 60 dB dynamic range with higher frequency (10 kHz). This simulating program provides a quantitative analytical means for studying acoustic reflection imaging tool and development of the data processing and interpretation methods.     

低渗砂泥岩地层的声波测井波形特征

用于低渗砂岩油气开发的压裂需要地层的纵横波时差参数;压裂形成的人工裂缝或砂泥岩中的自然裂缝是油气流入井中的主要通道.声波测井是测量地层的纵横波时差、探测连通裂缝的有效方法.本文给出了新设计的偶极子组合声波测井仪器在低渗砂泥岩地层所测量的波形.新仪器利用声波在井内传播的固有频率设计探头,利用探头阻抗随频率的变化曲线设计匹配电路,有效地提高了声发射功率,在井下采用16位A/D转换器提高了原始波形的采样精度.测量的波形比通常的声波测井波形所含的地层信息丰富.单极子长源距声波测井波形中除了通常的纵、横波和Stoneley波外,还发现了起始于纵波首波的反射波以及随深度改变其到达时间的后续波.通过与其它测井曲线综合对比发现:其反射波与井眼扩径有关,使测量的声波时差加大;后续波在显示地层特征的同时,与横波混叠在一起,使横波时差的提取变得困难.另外,在偶极子所测量的低频波形中也发现了明显的反映地层连通裂缝的反射波,该反射波可以用于指示低渗砂岩地层中的连通裂缝.     

偶极子声波测井在双层套管井中激发的声场特征

多层套管井中的套管-水泥以及水泥-地层在径向上的分布较单层套管井界面多且径向分布范围更深.目前还未有针对多层套管井固井质量评价的有效方法,偶极子阵列声波测井具有工作频率低径向探测深的优势,本文通过实轴积分和三维有限差分计算了偶极子声源在不同胶结状况的双层套管井中激发的声场,旨在探讨偶极子阵列声波在双层套管井中评价外层套管和地层之间胶结状况的可行性.计算结果表明,在两层套管之间充填泥浆时,全波波形中可明显观测到地层弯曲波和套管弯曲波,若外层套管和地层之间充填水泥套管弯曲波的慢度明显高于充填泥浆时的慢度值;在外层套管和地层之间胶结的水泥出现缺失,套管弯曲波的频散曲线会处于完全胶结与完全无胶结之间,其变化量与水泥在外层环空的占比近似成线性关系,几乎不受水泥缺失的径向位置影响.由于偶极子声源的指向性,套管弯曲波的响应特征主要受与声源振动指向的方位扇区胶结状况的影响;在两层套管之间充填水泥时,若外层套管和地层之间充填泥浆,仍可观测到套管弯曲波,模拟波列提取的套管弯曲波的频散曲线与解析解计算的频散曲线吻合,两层套管通过水泥胶结在一起,其弯曲振动的传播速度稍有增加,类似单层套管厚度加厚的响应特征.两种模拟方法得到的结果均表明套管弯曲波对外层套管和地层之间的充填介质敏感,验证了偶极子阵列声波测井在双层套管井固井质量评价中应用的可行性.

     

阵列感应(AIT)测井响应的连续反演

构建了基于分辨率匹配信号处理的AIT反演方法.该方法用数值模式匹配方法计算测井响应,并进行真分辨率合成和分辨率匹配处理,再将分辨率匹配后的数据直接用于反演计算.该方法实现了对测井曲线的自动分层和连续反演,而且分辨率匹配曲线具有相同的分辨率和较好的围岩校正,有利于提高反演的可信度.     

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.     

地震资料波阻抗多尺度融合反演

本文依据最佳约束反演策略,应用波动方程多尺度反演方法反演中频段波阻抗.利用神经网络技术实现对中频段波阻抗的标定及高、中、低频波阻抗的融合,避免了因高、中、低频波阻抗频段重叠及中频段波阻抗精度低所造成的反演结果的误差.这种方法具有对初始模型依赖性低、收敛速度快、反演结果稳定、抗噪能力强等特点.