声电效应测井的有限差分模拟

本文研究声电效应测井波场的有限差分模拟算法.忽略井外地层中诱导电磁场对孔隙弹性波的影响,将求解动电耦合波方程组的问题解耦,先计算孔隙弹性波,再计算其诱导电磁场.基于轴对称柱坐标系下的速度-应力交错网格,采用时域有限差分计算井孔流体声波和井外地层孔隙弹性波.将电磁场近似看作似稳场,基于轴对称柱坐标系下的5点式有限差分网格,求解不同时刻的电位Poisson方程,计算诱导电场.结果表明:本文算法可准确模拟频率6.0 kHz的声电效应测井全波;在声波测井频率范围内,电导率、动电耦合系数和动态渗透率的低频近似对伴随电磁场的计算影响不大;地层水平界面导致伴随反射斯通利波的电场和显著的界面电磁波,后者对于探测地层界面具有潜在的应用价值.     

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

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

含流体裂缝介质中地震波场数值模拟

油气勘探开发实践证明,裂缝常常是油气藏存储的空间或运移的通道,因此,裂缝各向异性介质中地震波场的研究越来越倍受关注,国内外很多岩石物理学者、地球物理专家等对裂缝信息的描述提出了很多理论认识与方法技术.本文根据Eshelby-Cheng各向异性裂缝介质模型理论,求取各向异性裂缝介质的弹性参数,并建立Eshelby-Cheng各向异性裂缝介质的波动方程,利用时间错格伪谱法对含流体裂缝介质进行数值模拟,模拟结果表明,采用时间错格伪谱法能有效解决各向异性介质的波场传播,利用时间错格有限差分算子替代普通的差分算子来求解时间导数,利用快速傅氏变换求解空间导数,大大提高了正演模拟的计算精度与计算效率.并且与各向同性介质相比,地震波在含流体裂缝各向异性介质中的传播要复杂得多,各向同性介质层中的波是纯的,其横波不会发生分裂,而在各向异性介质层中,横波将发生分裂.     

利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播

烃类储集层是一种复合多相介质，在固体颗粒的空隙中含有气体或液体. 研究弹性波在该类地层中的传播规律对于油气勘探开发，特别对于全波列声波测井有重要意义. 为了提高孔隙弹性介质数值模拟的计算效率，本文采用改进显式交错网格有限差分算法取代常用的空间域四阶和时间域二阶的速度 - 应力有限差分算法，算法的空间域为八阶、时间域为二阶. 虽然计算的时间步长略小于空间域四阶的情形，但高阶有限差分算法可以选择较粗糙的网格，因此补偿了计算的低效；同时高阶交错网格有限差分算法的空间频散性比低阶算法小. 利用该算法计算了一个两层模型的波场，同时还模拟了等效弹性和孔隙弹性模型中波的传播. 结果表明慢波及其影响明显，尽管慢波衰减很快，但被某一界面反射后，转换形成的P波和S波仍以正常的方式传播，且比慢波衰减小.     

裂缝诱导各向异性双孔隙介质波场传播特征

基于裂缝诱导各向异性和双相介质理论,对裂缝诱导的具有水平对称轴的横向各向同性（HTI）双孔隙介质的本构关系进行了研究,与等效连续介质模型相结合,综合考虑裂缝系统和基质孔隙系统的两种孔隙度和两种渗透率参数,得到裂缝诱导HTI双孔隙介质的等效孔隙度和等效渗透率,进而得到介质的运动平衡方程;并进一步推导出介质的一阶速度-应力方程.采用交错网格高阶有限差分法对模型进行了数值模拟,结果揭示了介质中两套系统的存在对其波场传播特征的影响,为进一步研究实际地球介质的波场特征奠定了基础.     

TTI介质的交错网格伪P波正演方法

研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果.     

任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟

给出了任意倾斜各向异性介质中二维三分量一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式及其稳定性条件,并推导出了二维任意倾斜各向异性介质完全匹配吸收层法边界条件公式和相应的交错网格任意偶数阶精度差分格式. 数值模拟结果表明,该方法模拟精度高,计算效率高,边界吸收效果好. 各向异性介质中弹性波波前面形态复杂, 且qP波波速不总是比qS波波速快. qS波波前面和同相轴的三分叉现象普遍, 且其同相轴一般不是双曲线型. 当TI介质倾斜时,3个分量上均能够观测到横波分裂现象, 而且各波形的同相轴变得不对称.      

径向分层TI孔隙介质井孔中激发的模式波的数值研究

横向各向同性介质是地层中普遍存在的一种各向异性介质.本文对径向分层TI孔隙介质包围井孔中激发的斯通利波和弯曲波的传播特性进行了理论计算,发现模式波在低频时更多的是反应原状地层的信息,而随着频率的增加侵入带参数逐渐起控制作用;Biot理论描述的地层衰减比速度更容易受井壁附近地层参数的影响.利用灵敏度曲线定量研究了不同频率下地层各个参数对相速度和衰减系数的贡献大小,主要结果显示模式波的衰减受水平渗透率影响明显,而垂直渗透率的变化对模式波几乎无影响;斯通利波对水平向传播的横波速度比弯曲波的灵敏度高.从单极子和偶极子声源在井孔中激发的全波波形也可发现,声波测井仪器较宽的声源频带和合适的源距设置有利于对不同径向深度上的地层声学参数进行成像.     

各向异性双相介质弹性波场褶积算法数值模拟

从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到.     

基于斯通利波能量衰减的裂缝渗透性评价方法研究

塔里木库车地区气藏目的层主要为致密砂砾储层,裂缝是流体主要的运移通道.该地区多采用油基泥浆钻井,导致电成像效果差,降低裂缝有效性评价精度.斯通利波对井中流体运动的敏感性,为利用斯通利波评价裂缝渗透性提供了可能:阵列声波资料提取的斯通利波能量经归一化后计算衰减系数,然后根据实验得到的衰减系数与裂缝宽度的关系反演出裂缝宽度,同时根据裂缝宽度、裂缝倾角反演得到裂缝渗透率.结果表明,水基泥浆井中利用该方法反演得到的衰减系数、裂缝宽度、裂缝渗透率曲线与电成像解释成果图对应良好,在裂缝发育段曲线呈现明显高值,在裂缝不发育段呈现低值;同时与试油结果吻合,可以用来评价裂缝渗透性.最后将该法推广到油基泥浆井中,结合电成像测井可提高裂缝有效性评价的精度.     

BOREHOLE STONELEY WAVE PROPAGATION ACROSS PERMEABLE STRUCTURES1

This study investigates the propagation of borehole Stoneley waves across permeable structures. By modelling the structure as a zone intersecting the borehole, a simple 1D theory is formulated to treat the interaction of the Stoneley wave with the structure. This is possible because the Stoneley wave is a guided wave, with no geometric spreading as it propagates along the borehole. The interaction occurs because the zone and the surrounding formation possess different Stoneley wavenumbers. Given appropriate representations of the wavenum-ber, the theory can be applied to treat a variety of structures, including a fluid-filled fracture. Of special interest are the cases of permeable porous zones and fracture zones. The results show that, while Stoneley wave reflections are generated, strong Stoneley wave attenuation is produced across a very permeable zone. This result is particularly important in explaining the observed strong Stoneley wave attenuation at major fractures where it has been difficult to explain the attenuation in terms of the single planar fracture theory. In addition, by using a simple and sufficiently accurate theory to model the effects of the permeable zone, a fast and efficient method is developed to characterize the fluid transport properties of a permeable fracture zone.     

Stoneley wave generation by an incident P-wave propagating in the surrounding formation across a horizontal fluid-filled fracture

This study investigates the generation of the low-frequency borehole Stoneley wave (tube wave) by a plane P-wave propagating through the surrounding elastic formation, which is intersected by a fluid-filled fracture. A model is constructed taking into account the dynamic fluid coupling between the borehole interior and the fluid-filled fracture of infinite extent with parallel walls. The basic mechanism of such coupling is due to the contraction of the fracture walls by the incident P-wave, leading to seismic radiation into the fracture. The dynamic fluid flux from the fracture into the borehole interior, and vice versa, is the source of the low-frequency Stoneley wave. An expression for the monopole pressure source, exciting the tube wave, is obtained. The tube-wave equation in the long-wave approximation is derived in the presence of a fluid-filled fracture of infinite extent. Amplitudes and waveforms of Stoneley waves are analysed in the seismic wavelength range for P-wave pulses of various shapes. It is shown that the amplitude and waveform of the Stoneley wave depends significantly on the two dimensionless parameters of the problem: (1) the ratio of the borehole radius to the dominant wavelength of the incident pulse; (2) the ratio of the fracture width to the borehole radius. It is found that the amplitude of the generated Stoneley wave can be of the order of the P-wave amplitude in the borehole fluid. Stoneley waveforms are found to be completely different from those of the incident pulse.     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.     

Simulation of full-waveform log in saturated cracked formations using Hudson's approach

We study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix. We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes). Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.     

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

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

Analytical expressions for the borehole seismic signatures excited by an explosive source buried in an elastic half‐space

Borehole guided waves that are excited by explosive sources outside of the borehole are important for interpreting borehole seismic surveys and for rock property inversion workflows. Borehole seismograms are typically modelled using numerical methods of wave propagation. In order to benchmark such numerical algorithms and partially to interpret the results of modelling, an analytical methodology is presented here to compute synthetic seismograms. The specific setup is a wavefield emanating from a monopole point source embedded within a homogeneous elastic medium that interacts with a fluid‐filled borehole and a free surface. The methodology assumes that the wavelength of the seismic signal is much larger than the borehole radius. In this paper, it is supposed that there is no poroelastic coupling between the formation and the borehole. The total wavefield solution consists of P, PP, and PS body waves; the surface Rayleigh wave; and the low‐frequency guided Stoneley wave (often referred as the tube wave) within the borehole. In its turn, the tube wave consists of the partial responses generated by the incident P‐wave and the reflected PP and PS body waves at the borehole mouth and by the Rayleigh wave, as well as the Stoneley wave eigenmode. The Mach tube wave, which is a conic tube wave, additionally appears in the Mach cone in a slow formation with the tube‐wave velocity greater than the shear one. The conditions of appearance of the Mach wave in a slow formation are formulated. It is shown that the amplitude of the Mach tube wave strongly depends on Poisson's ratio of the slow surrounding formation. The amplitude of the Mach tube wave exponentially decreases when the source depth grows for weakly compressible elastic media with Poisson's ratio close to 0.5 (i.e., saturated clays and saturated clay soils). Asymptotic expressions are also provided to compute the wavefield amplitudes for different combinations of source depth and source‐well offset. These expressions allow an approximate solution of the wavefield to be computed much faster (within several seconds) than directly computing the implicit integrals arising from the analytical formulation.     

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

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

Estimating the permeability of the Nojima Fault Zone by a hydrophone vertical seismic profiling experiment

Abstract A multi-offset hydrophone vertical seismic profiling (VSP) experiment was done in a 747 m deep borehole at Nojima Hirabayashi, Hyogo prefecture, Japan. The borehole was drilled to penetrate the Nojima Fault, which was active in the 1995 Hyogo-ken Nanbu earthquake. The purpose of the hydrophone VSP is to detect subsurface permeable fractures and permeable zones and, in the present case, to estimate the permeability of the Nojima Fault. The analysis was based on a model by which tube waves are generated when incident P-waves compress the permeable fractures (or permeable zones) intersecting the borehole and a fluid in the fracture is injected into the borehole. Permeable fractures (or permeable zones) are detected at the depths of tube wave generation, and fracture permeability is calculated from the amplitude ratio of tube wave to incident P-wave. Several generations of tube waves were detected from the VSP sections. Distinct tube waves were generated at depths of the fault zone that are characterized by altered and deformed granodiorite with a fault gouge, suggesting that permeable fractures and permeable zones exist in the fault zone. Tube wave analysis shows that the permeability of the fault gouge from 624 m to 625 m is estimated to be approximately 2 × 10⁻¹² m².     

利用井中低频偶极横波进行声波远探测的新方法

为了突破目前声波远探测技术存在的局限性,提出了一种新的偶极横波远探测方法,即利用井中偶极子产生的井中弯曲波存在低频截止频率的现象,在声源截止频率以下激发偶极声波.通过对比分析井中偶极声源分别在截止频率上、下激发时,井孔内外产生的辐射声场,明确了截止频率以下井中偶极声源的远场辐射特征和低频截止频率激发偶极横波的优势,结合数值模拟,进一步对其反射声场进行了分析.结果表明,该方法可以避免艾里相的巨大振幅对数据量化产生的"饱和"效应,相比传统的远探测测井方式更具优势,常规源距即可满足专门的远探测测井仪器需求.     

利用随钻正交偶极子声波测井评价地层各向异性的数值研究

在随钻测井条件下,由于钻铤占据了井孔内的大部分空间,充液井孔中沿着井轴方向传播的模式波的特性与电缆测井非常不同.本文建立了随钻正交偶极子测井声学模型,采用三维有限差分方法模拟了偶极子声源在随钻条件下各向异性地层井孔内激发的声场,研究了地层的声学各向异性在随钻正交偶极子声波测井中的响应特征.数值模拟结果表明,在随钻测井条件下,对于井轴同TI地层对称轴垂直的情况,弯曲波分裂现象仍然存在,通过正交偶极子测量方式和合适的反演算法能够准确有效地确定地层的快横波方位角,可以考虑采用同正演理论相结合的反演算法来获得地层的快、慢横波速度及声学各向异性信息;对于井轴同介质对称轴呈一定夹角的TI地层井孔,情况变得非常复杂,不同井斜倾角下弯曲波的速度的变化趋势并非同对应的地层横波速度的变化趋势完全一致,不过在一定的频段内,地层横波速度仍然是弯曲波的最主要控制因素.对于本文研究的模型,当井轴同介质对称轴的夹角大于大于60°时,此时获得的弯曲波的各向异性值基本能够反映对应角度下地层横波速度的各向异性信息.