三维数字岩心透射超声波模拟与速度精度分析（英文）

基于数字岩心的透射超声波模拟是数字岩石物理的重要组成部分,可用于研究波在孔隙岩石中的传播规律,计算岩石的等效速度等。在三维数字岩心中模拟超声波传播时,沿波传播方向在数字岩心两侧附加均匀的镶边层,设置震源和检波器。震源激发后检波器记录数字岩心的入射波和透射波,拾取两者波峰时间的差值,计算超声波在数字岩心中的传播速度(即数字岩心的速度)。为评估透射超声波模拟的精度,将一个数字岩心的孔隙分别饱和气、油和水,计算它们的速度。当震源频率较高时,速度随频率的降低而升高,这应当是波场散射的结果。当孔隙流体按气→油→水方向变化时,不同频率对应的速度具有相似的变化规律,并与线弹性静力学模拟输出的速度的变化规律相似,但数值有差异。考虑到线弹性静力学模拟已被广泛应用,若将其作为基准,则可认为透射超声波模拟具有较高的相对精度。     

含裂缝VTI介质的弹性波传播特性研究

有利油气储集的层状沉积地层中广泛存在裂缝,属于含裂缝VTI介质.综合考虑层状介质和裂缝介质的特点,基于几何相似原理,设计并制作了含定向排列裂缝的VTI介质物理模型,实验研究了测量频率、裂缝特征尺寸、裂缝密度对弹性波透射的影响.结果表明,波垂直裂缝面传播时,纵波首波速度随测量频率的增加而线性增大,峰值速度基本不变,单纯由裂缝导致的衰减以良好的乘幂规律降低;纵波波速随裂缝密度、裂缝特征尺寸的增大而线性增大;衰减随裂缝密度增大而增大,随裂缝特征尺寸的增大而减小.当接收探头波长与特征尺寸比为1~12时,纵波首波速度变化显著,处于射线理论与等效介质理论描述的过渡区间.纵横波速比对裂缝特征尺寸变化敏感,可作为典型的裂缝监测参数,也为基于弹性波相应的预测岩石孔隙结构提供了依据.     

热弹性P波在有流体夹层介质中的传播问题及其基本特征*

本文作为基础理论研究,讨论了热弹性P波在有流体夹层介质中的传播问题。研究结果表明:热弹性波在流体夹层的传播中,不仅存在反射波、透射波,而且,还伴随有具有相同传播速度的温度波;且其反射系数、透射系数均为复数,并与介质的物性参数及夹层厚度有关。同时指出反射波、透射波的振幅、位相均受介质的物性参数及夹层厚度的影响。此外,还表明夹层中往返震荡的层间波是一系列正传热弹性波和反传热弹性波的叠加。     

热弹性P波在有流体夹层介质中的传播问题及其基本特征

本文作为基础理论研究,讨论了热弹性P波在有流体夹层介质中的传播问题。研究结果表明:热弹性波在流体夹层的传播中,不仅存在反射波、透射波,而且,还伴随有具有相同传播速度的温度波;且其反射系数、透射系数均为复数,并与介质的物性参数及夹层厚度有关。同时指出反射波、透射波的振幅、位相均受介质的物性参数及夹层厚度的影响。此外,还表明夹层中往返震荡的层间波是一系列正传热弹性波和反传热弹性波的叠加。     

饱和多孔微极介质的波动方程及其势函数方程

土是由一定尺寸大小颗粒所构成的多孔介质,具有明显的颗粒特性,当土颗粒间的孔隙被流体（如水或油）充满时则成为饱和土.利用微极理论和Biot波动理论的研究成果,把饱和土中多孔固体骨架部分近似地视为微极介质,孔隙中的流体部分视为质点介质,获得饱和多孔微极介质的弹性波动方程.借鉴Greetsma理论,建立了饱和多孔微极介质弹性本构方程力学参数与相应单相介质弹性参数的相互关系,使饱和多孔微极介质弹性波动方程中的物理参数具有明确的物理意义,易于在试验中确定.运用场论理论把饱和多孔微极介质的波动方程简化为势函数方程,建立了饱和多孔微极介质中五种弹性波的弥散方程,数值分析了五种简谐体波在无限饱和多孔微极介质中的传播特性. 结果表明,P1波、P2波和剪切S1波的波速弥散曲线与经典饱和多孔介质基本相同,当频率小于临界频率ω0时旋转纵波θ波和横波S2波不存在,当频率大于临界频率ω0时,θ波和S2波的传播速度随频率增加而减小.     

P波入射下断层参数对地表地震动的影响

基于应用透射人工边界条件的显式有限元法计算断层破碎带宽度及力学参数变化、地震动入射角变化时二维断层场地模型P波入射下地表地震动场的分布。结果表明:(1)低速度破碎带的存在导致整个场地都有P波转换为SV波的分量,且在断层破碎带的区域出现断层陷波;(2)低速度破碎带的存在使输入场地恒定的能量向破碎带集聚放大,随着破碎带宽度增大或其介质波速降低集聚放大效应增大;(3)场地放大效应是频率相关的,宽度较宽或介质波速较低的断层破碎带对输入地震动中较低的频率成份放大显著;(4)竖向断层破碎带能阻隔斜入射地震P波,随着入射角增加隔震效应更显著。     

实验样品波速测量的尺度效应和环境效应

本文用超声波地震模型实验方法,对引起有机玻璃样品P波速变化的尺度因素和环境因素进行了实验研究。结果表明,当样品厚度小于一个波长时,样品尺度和环境介质对样品P波速都有明显影响。当样品厚度由使用换能器的1/16个波长增加到3个波长时,P波速变化约16％。在低于样品P波速的液体中测量,样品的P波速增加,增加的幅度与液体的波速有关。文中还对纵波和拉伸波的转换现象作了初步探讨。     

昆明场地脉动优势频率、等效剪切波速度和场地类别判定

以云南省昆明地区为例,对28个钻孔分别以20 m、25 m、30 m厚度计算等效剪切波速和卓越频率,同时测定场地脉动优势频率.结果显示:以20 m、25 m、30 m厚度计算的等效剪切波速,其后者一般都大于前者.对多数钻孔,用25 m厚等效剪切波速和卓越频率判定的场地土类别一致;少数钻孔在靠近30 m时二者判定结果一致.经测定,场地脉动优势频率与20 m厚波速卓越频率相近,但却明显高于25 m厚波速卓越频率.脉动优势频率与不同计算厚度的等效剪切波速度相关性基本相同,对同一厚度(深度)脉动优势频率随等效剪切波速度增加而增加.若等效剪切波速度相等,则深度小的脉动优势频率高.由此推出,脉动优势频率主要由地表层20 m厚岩土力学性质决定,而且越靠近表层的岩土力学性质对脉动优势频率的影响越大.本文从弹性力学理论证明了脉动优势频率和剪切波速度的关系式.通过进一步分析证明,用25 m厚等效剪切波速判定场地土类别更可靠,用脉动优势频率判定场地土类别可作为有效的辅助方法.它们将影响对场地类别的判定.     

零偏VSP反演Q值CFS方法及影响因素研究

当地震波在地层中传播时,能量的衰减是传播介质非弹性性质的反映.这种介质所固有的衰减特性通常用品质因子(Quality Factor, Q)来描述.本文基于地震波在传播过程中主频降低的衰减特性,利用质心频率偏移法(Centroid Frequency Shift,简称 CFS)对零偏VSP正演模型进行Q值计算,分析薄层、上下行波场、频带宽度、反演波速、震源位置等因素的影响.模型结果显示:CFS法比频谱比法、振幅衰减法能更加准确识别薄层界面,可准确反演出厚度10m的地层;上行波场的加入降低反演结果准确度,尤其对高Q值层;当有效波频带增宽时,高频部分衰减明显,Q值反演结果接近理论值;反演波速误差对反演准确性影响很大.最后,利用CFS对实际井旁地面地震资料进行Q值反演,反Q滤波后其主频由25Hz提高到35Hz,分辨率得到改善,证明了CFS方法的实用性.     

透射波波场重建和介质层厚反演研究

本文以P波和SV波波场比值,提出层状介质情况下透射波波场的反向重建算法原理,并据此给出在已知介质速度值时,求介质层厚度的反演方法。文中还以深源远震记录波形资料给出反演算例,证实本文的算法是可行的。     

岩石声学参数的实验测量及研究

通过对天然岩芯及人工岩芯的测量，研究２４０ｋＨｚ-１．ＳＭＨｚ频段内的声频散、温度对岩石声速的影响、岩石含水饱和度对岩石声速和声衰减的影响，给出岩石中声衰减与声波频率间关系的最佳拟合，并讨论了裂缝对岩石中声波传播的影响。     

地层波速的电缆和随钻动电测井实验研究

基于孔隙介质动电效应的测量方法在油气勘探领域有着巨大的潜在应用价值.本文针对理论预测的随钻动电测井可削弱钻铤波对地层声波干扰的问题,设计了室内随钻动电测井实验装置和测量系统,并在小尺寸砂岩模型井中开展了随钻动电测井实验研究.为了突出随钻动电测井技术的特点,本文先在砂岩模型井中进行了电缆声波测井和电缆动电测井实验测量,获得了地层声波的传播速度.然后又在同一模型井中进行了随钻声波测井和随钻动电测井实验,记录了这两种情况下的测井全波波形,比较了电缆声波/动电测井和随钻声波/动电测井的差异性,并进一步分析了钻铤波和井中动电转换信号的传播特性.本文还通过随钻动电测井的时域波形提取到被测地层的纵横波速度,从实验角度验证了随钻动电测井技术测量地层波速的可行性.此外,本文实验结果对随钻动电测井的仪器设计及现场测量也有参考意义.     

Numerical analysis of one-dimensional nonlinear acoustic wave

Numerical investigations on one-dimensional nonlinear acoustic wave with third and fourth order nonlinearities are presented using high-order finite-difference (HFD) operators with a simple flux-limiter (SFL) algorithm. As shown by our numerical tests, the HFDSFL method is able to produce more stable, accurate and conservative solutions to the nonlinear acoustic waves than those computed by finite-difference combined with the flux-corrected-transport algorithm. Unlike the linear acoustic waves, the nonlinear acoustic waves have variable phase velocity and waveform both in time-space (t-x) domain and frequency-wavenumber (f-k) domain; of our special interest is the behaviour during the propagation of nonlinear acoustic waves: the waveforms are strongly linked to the type of medium nonlinearities, generation of harmonics, frequency and wavenumber peak shifts. In seismic sense, these characteristics of nonlinear wave will introduce new issues during such seismic processing as Normal Moveout and f-k filter. Moreover, as shown by our numerical experiment for a four-layer model, the nonlinearities of media will introduce extra velocity errors in seismic velocity inversion.     

Acoustic log simulation in a viscoelastic formation: Cole–Cole model

It is known that sedimentary rocks demonstrate velocity dispersion in the acoustic log frequency range. In this paper we have calculated the waveforms of sonic log for a borehole located in a viscoelastic medium. The acoustic field in the borehole has been obtained for acoustic multipole sources. To describe the viscoelastic properties of a rock we used the Cole–Cole model. This model describes the dispersion of acoustic wave velocities and quality factors in a wide frequency range. To solve the acoustic log direct problem we have applied the double integral Fourier transform (RAI method). The results obtained have shown the feasibility of S-wave velocity dispersion estimation from acoustic dipole waveform processing.     

基于简化的Pride理论模拟声电效应测井响应

针对声电效应测井问题，提出了一种全波列数值模拟方法.该方法忽略声电效应测井 时转换电场对声场的影响,并将电场视为似稳场.采用点声源模型，依据Biot理论得出了井外 孔隙介质声场的表达式，运用这些表达式和似稳电场近似方法，导出了声电效应测井时转换 电场的计算公式.在计算出的转换电场波形中，有伴随斯通利波的电场、伴随纵波和横波的 电场、和临界折射电磁波场.在25kHz以下的频率范围内，依据这种方法计算出来的声电转换 波波形与依据完整Pride理论计算的波形一致.     

随钻动电测井中声诱导电场的理论模拟

为防止沿钻铤传播的声信号掩盖来自地层的声信号,目前需在钻铤上刻槽以实现对地层纵、横波速度的测量.但刻槽不仅占用钻铤空间,且降低了钻铤的强度和刚度.本文探讨在不刻槽条件下,利用动电效应获得地层纵横波速度的可行性.为此,模拟了随钻声波测井信号及由岩石动电效应产生的电场信号.先视地层为黏弹性介质,模拟随钻声波测井信号,再利用围压与孔隙流体压强的关系获得了声波在地层中引起的孔隙压强,然后计算由孔隙地层动电效应产生的声诱导电磁场.对多极声源激发的声场波形及其诱导电场的波形进行比较,发现随钻声诱导电场中钻铤波的相对幅度明显小于随钻声场中钻铤波的相对幅度.此外,发现在慢速地层条件下,偶极声源激发的电场中出现了比较明显的地层横波信号.     

多极源随钻声波测井实验分析

本文针对随钻声波测井中钻铤波干扰地层声波测量的问题,设计了小尺寸随钻声波测井探头,在实验室内开展了多极源(单极源、偶极源和四极源)随钻声波测井实验研究.先在水池中对裸露的探头进行了随钻测量,记录到了沿钻铤传播的直达钻铤波,并获得了不同声源激发的钻铤波速度.之后在砂岩和有机玻璃模型中进行了随钻声波测井实验,观测到多极源随钻测井的全波波形,并通过与水池中实验结果的对比,分析了井中钻铤波的传播特性.特别地,在偶极随钻测井实验中不仅记录到了偶极一阶钻铤波,还观测到对弯曲波干扰较大的偶极二阶钻铤波,进而研究了不同声源频率下该波群的响应特性及其对弯曲波测量的影响.此外,基于单极源和偶极源随钻声波测井实验数据,本文发现:随着声源频率的增加,单极和偶极钻铤波的传播特性不同,但它们在测井全波中的相对幅度均降低,进而可从测井全波中较好地提取地层的声波速度.本文实验结果对随钻声波测井仪器设计及测井数据解释具有重要意义.     

套管井贴壁声源激发的非轴对称声场的有限差分模拟及结果分析

套管井贴壁声波测井仪器在固井质量评价中的应用越来越广泛.本文利用三维交错网格高阶有限差分方法首次研究了贴壁声源在套管井中激发的非轴对称声场特征,给出了测井声源在贴壁固体弹性介质中激发时的数值模拟方法,模拟的全波波形与测井仪器在实际井中记录的波形吻合较好,这为利用数值模拟方式开发贴壁声波仪器的潜在应用奠定了基础.通过对套管井不同胶结状况下的数值模拟,得到了贴壁声源在套管井中激发的声场特征和传播规律,模拟结果显示加载声源的贴壁滑板的存在大大压制了套管中弯曲型Lamb波的幅度;套管中传播的拉伸波的幅度相对较强,其泄漏到水泥中的声波传播到水泥和地层界面时发生反射,反射波的幅度和相位特征携带了水泥环第二界面的胶结信息.数值计算结果为利用贴壁声源激发接收的波形实现水泥环第二界面井周不同角度上的胶结评价奠定了理论基础.     

Experimental study on the performance of an azimuthal acoustic receiver sonde for a downhole tool

A new azimuthal acoustic receiver sonde with a body and corresponding circuits was designed for a downhole tool. The 64‐sensor receiver sonde holds eight receiver stations that can be combined into at least 64 three‐sensor receiver subarrays. As a result, the receiver sonde can use different sensor combinations instead of different transducer types to produce multiple modes, including a phased azimuthal reception mode and conventional monopole, dipole, and quadruple modes. Laboratory measurements were conducted to study the performance of the azimuthal acoustic receiver sonde for a downhole tool, and the experimental results indicate that the receiver sonde provides a consistent reception performance. Individual sensors receive similar time‐domain waveforms, and their corresponding frequency bands and sensitivities are consistent within the measurement errors of around 5%. The direction of the reception main lobe is approximately parallel to its exterior normal direction. In addition, a receiver subarray with three sensors receives waveforms that have higher energy and narrower beamwidths. For individual sensors, the angular width of the dominant reception lobe is 191.3^° on average, whereas that of the individual receiver subarrays is approximately 52.1^° on average. The amplitude of the first arrival received by the receiver subarray centred at the primary sensor directly pointing to the source is approximately 2.2 times the average amplitude of the first arrivals received by the other receiver subarrays in the same receiver station. Thus, the maximum amplitude of the waveforms received by the receiver subarrays can be used to determine the direction of the incident waves. This approach represents a promising method for determining the reflector azimuth for acoustic reflection logging and three‐dimensional acoustic logging.     

Modelling Seismic Waves Around Underground Openings in Fractured Rock

—?The potential for large excavation-induced seismic events may be recognised, even if the timing of an event may be inherently unpredictable. In this case, modelling the wave propagation from a potential event could allow the dynamic motions around an excavation to be projected, and for areas of danger to be anticipated. However, the above and other potential applications require accurate models of wave interaction with the openings, as well as with the fractured rock which surrounds such excavations. This paper considers real recorded waveforms and how well these waveforms are modelled by explicit mechanical models of the source, the medium and the excavation. Models of experiments at three different scales of the problem are presented: small and large amplitude waveforms recorded around a deep-level mining tunnel in a synthetic rockburst experiment; waveforms from laboratory experiments of waves through plates of steel representing fractures; waveforms from active pulses in an acoustic emission experiment in a small volume of fractured rock at the surface of an underground excavation. The results show that elastic wave propagation around an excavation was a first approximation for small amplitude waves, but was less successful for modelling large amplitude waves and more fractured rock. Fractures in the models were represented explicitly with displacement discontinuities. Waveforms through known fracture geometries were particularly well-reproduced, and indicate the importance of fracture stiffness, the in situ stress state, and stress-dependence of the fractures in such models. Overall, the models are sufficiently successful at representing recorded behaviour, to be encouraging for the goal of representing accurate wave motions around excavations.