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

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

Acoustic reverberation in a logging tool-borehole-saturated porous medium system

The influence of the elastic and hydrodynamic properties of a rock on the time attenuation coefficient of an acoustic wave, which is reflected from the borehole wall, is considered using the methods of the mechanics of saturated porous media. The calculations were performed for a system consisting of an acoustic logging tool of a finite with a finite diameter, a fluid-filled borehole, and a porous permeable rock. The performed simulation showed that in rocks with a low hydrodynamic permeability, the acoustic-reverberation time is determined by the acoustic impedance of the borehole wall (product of the rock density and the longitudinal-wave velocity in it). In the case of rocks with a permeability of about several hundred millidarcy, the time signal’s attenuation coefficient substantially depends on the rock’s permeability.     

Acoustic reflection logging for open and cased holes

In the present work, the waveforms of reflected wave sonic log for open and cased boreholes are calculated. Calculations are performed for a borehole containing an acoustic multipole source (monopole, dipole, or quadrupole). A reflected wave is more efficiently excited at resonant frequencies. These frequencies for all source types are close to the frequencies of oscillations of a fluid column located in an absolutely rigid hollow cylinder. It is shown that the acoustic reverberation is controlled by the acoustic impedance of the rock Z = V_p ρ_s for fixed parameters of the borehole fluid, where V_p is the compressional wave velocity in the rock, and ρ_s is the rock density. This result is correct for all types of acoustic sources (monopole, dipole, or quadrupole). Methods of the waveform processing for determining parameters characterizing the reflected wave are discussed.     

WIND observations of coherent electrostatic waves in the solar wind

The time domain sampler (TDS) experiment on WIND measures electric and magnetic wave forms with a sampling rate which reaches 120 000 points per second. We analyse here observations made in the solar wind near the Lagrange point L1. In the range of frequencies above the proton plasma frequency f_pi and smaller than or of the order of the electron plasma frequency f_pe, TDS observed three kinds of electrostatic (e.s.) waves: coherent wave packets of Langmuir waves with frequencies f ≃ f_pe, coherent wave packets with frequencies in the ion acoustic range f_pi ≥ f ≥ f_pe, and more or less isolated non-sinusoidal spikes lasting less than 1 ms. We confirm that the observed frequency of the low frequency (LF) ion acoustic wave packets is dominated by the Doppler effect: the wavelengths are short, 10 to 50 electron Debye lengths λ_D. The electric field in the isolated electrostatic structures (IES) and in the LF wave packets is more or less aligned with the solar wind magnetic field. Across the IES, which have a spatial width of the order of ≃25_D, there is a small but finite electric potential drop, implying an average electric field generally directed away from the Sun. The IES wave forms, which have not been previously reported in the solar wind, are similar, although with a smaller amplitude, to the weak double layers observed in the auroral regions, and to the electrostatic solitary waves observed in other regions in the magnetosphere. We have also studied the solar wind conditions which favour the occurrence of the three kinds of waves: all these e.s. waves are observed more or less continuously in the whole solar wind (except in the densest regions where a parasite prevents the TDS observations). The type (wave packet or IES) of the observed LF waves is mainly determined by the proton temperature and by the direction of the magnetic field, which themselves depend on the latitude of WIND with respect to the heliospheric current sheet.     

Geometrical spreading factor of head waves in cased and open boreholes

Factors (coefficients) of geometrical spreading of compressional and shear head waves are calculated for an impulse multipole source of elastic oscillations in boreholes. It is shown that the length of the logging tool (i.e., the distance between the source and the nearest receiver) used for sonic measurements and the velocities of elastic waves in the medium both contribute to the factor of geometrical spreading. For a high-velocity formation (the shear wave velocity in the rock is higher than the compressional wave velocity in the fluid that fills the borehole) and a sufficiently long sonic tool with a monopole source, the coefficient of geometrical spreading is approximated by asymptotic formula 1/Z [Roever et al., 1974; Krauklis and Krauklis, 1976], where Z is the length of the tool; i.e., the amplitude of the compressional head wave decreases proportionally to the distance between the source and the receiver. In acoustically soft formations, this approximation is inapplicable even for long tools with length Z > 4 m. Waveforms in cased boreholes have a significant frequency dispersion even in case of good-quality cementing, and the factor of geometrical spreading there depends considerably on the length of the tool and the elastic properties of the rocks.     

THE INFLUENCE OF AN ASYMMETRY IN THE SEQUENCE ROCK/COAL/ROCK ON THE PROPAGATION OF RAYLEIGH SEAM WAVES*

Small offsets in hard coal seams can be detected with the aid of seam (channel) waves. Transmission and reflection of seam waves depend, among other parameters, upon the symmetry properties of the sequence rock/coal/rock. Two typical unsymmetrical sequences are found in European coal deposits: (i) coal seams with roof and floor of differing acoustic impedance and (ii) coal seams interlayered with rock and soil. Two-dimensional analog models with appropriate impedance contrasts are used to study the effect of the unsymmetrical layers upon the propagation of Rayleigh seam waves. Data analysis is based upon amplitude measurements both parallel and perpendicular to the layers and dispersion curves. The effect of unsymmetrical roof (rock 1) and floor (rock 2) was studied with models containing homogeneous coal seams. Leaky mode wave groups with phase velocities (c_R) in the range between the SV-wave velocities (β_{r1, 2}) of the two rock materials, i.e. β_r1≥c_R > β_r2, form a characteristic part of the Rayleigh seam wave signal. Using Knott's energy coefficient calculations it is shown that in that range energy leakage into the surrounding rock by refracted SV-waves is restricted to only one of the two interfaces, namely coal/rock 2. At the other interface, coal/rock 1, all waves are totally reflected. Thus, the high amplitudes of these leaky mode wave groups are explained by “quasi-normal mode” features. The influence of a dirt bed on wave propagation was studied in models where the roof and the floor have the same elastic properties. The maximum thickness of the dirt bed did not exceed 20% of the total seam thickness. The effect of the bed's location within the seam was also investigated. For all recorded normal-mode wave groups either the total seam or the coal layers could be regarded as wave guides. This was shown by the fact that the phases could be associated with the phase velocity dispersion curves calculated for the symmetrical sequence rock/coal/rock. These curves are relevant under the condition that the thickness of the coal layer assumed under the calculation coincides with the thickness of the effective wave guide of the respective wave groups. Wave groups guided in the total seam are not influenced by either the thickness or the position of the dirt bed. On the other hand, for wave groups guided in the coal layers, the quotient of signal amplitudes in the coal layers is influenced by the position of the dirt bed.     

T-phase Observations in Northern California: Acoustic to Seismic Coupling at a Weakly Elastic Boundary

—?Plans for a hydroacoustic network intended to monitor compliance with the CTBT call for the inclusion of five T-phase stations situated at optimal locations for the detection of seismic phases converted from ocean-borne T phases. We examine factors affecting the sensitivity of land-based stations to the seismic T phase. The acoustic to seismic coupling phenomenon is described by upslope propagation of an acoustic ray impinging at a sloping elastic wedge. We examine acoustic to seismic coupling characteristics for two cases; the first in which the shear velocity of the bottom is greater than the compressional velocity of the fluid (i.e., v _p > v _s > v _w ), the second is a weakly elastic solid in which v _s << v _w < v _p . The former is representative of velocities in solid rock, which might be encountered at volcanic islands; the latter is representative of marine sediments. For the case where v _s > v _w , we show that acoustic energy couples primarily to shear wave energy, except at very high slope angles. We show that the weakly elastic solid (i.e., v _s << v _w ) behaves nearly like a fluid bottom, with acoustic energy coupling to both P and S waves even at low slope angles.¶We examine converted T-wave arrivals at northern California seismic stations for two event clusters; one a series of earthquakes near the Hawaiian Islands, the other a series of nuclear tests conducted near the Tuamoto archipelago. Each cluster yielded characteristic arrivals at each station which were consistent from event to event within a cluster, but differed between clusters. The seismic T-phases consisted of both P- and S-wave arrivals, consistent with the conversion of acoustic to seismic energy at a gently sloping sediment-covered seafloor. In general, the amplitudes of the seismic T phases were highest for stations nearest the continental slope, where seafloor slopes are greatest, however noise levels decrease rapidly with increasing distance from the coastline, so that T-wave arrivals were observable at distances reaching several hundred kilometers from the coast. Signal-to-noise levels at the seismic stations are lower over the entire frequency spectrum than at the Pt. Sur hydrophone nearby, and decrease more rapidly with increasing frequency, particularly for stations furthest from the continental slope.     

裸眼井中弹性波传播的非对称模式的数值研究

本文从数值上研究了裸眼井中弹性波传播的非对称模式,给出了合成微地震图和导波（弯曲波）的频散曲线,发现在“硬”地层和“软”地层的井中,导波都是高度频散的,其最大相速度等于地层的横波速度,其截止频率低于对称模式的伪瑞利波的截止频率;在低频（2-3kHz）和长源距（3-4m）的条件下,由非对称的声源（如声偶极子）所产生的微地震图中,初至信号是以横波速度传播的,而以纵波速度传播的信号被抑制。本文的结果对研制横波速度测井仪是有意义的。     

充液井中多极声源的辐射效率

利用反射声场探测井外地质构造的声波远探测方法已逐渐成为石油测井的一门重要应用技术.在远探测测井中声源常采用单极、偶极和四极等.为了考察各种声学辐射器的辐射性能,本文从辐射声场和井孔声场的能流密度出发,提出了一种利用辐射器向井外辐射的能量与沿井筒传播的能量比的大小来评价其辐射效率的方法,考察了不同声源激发的沿井传播的导波能流、地层辐射波能流及声源辐射效率随频率的变化规律.计算结果表明:不同声源的辐射效率和优势激发频段各不相同.单极声源在低频下的辐射效率很低,因为此时声源激发的能量几乎全部被斯通利波带走,相比之下,偶极声源在低频时的辐射效率远大于单极声源.本文的结果说明了偶极声源作为低频远探测声源要优于单极声源.     

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

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

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.     

Nonlinear development of phase-mixed alfvén waves

Abstract

We derive an equation governing the nonlinear propagation of a linearly polarized Alfvén wave in a two-dimensional, anisotropic, slightly compressible, highly magnetized, viscous plasma, where nonlinearities arise from the interaction of the Alfvén wave with fast and slow magnetoacoustic waves. The phase mixing of such a wave has been suggested as a mechanism for heating the outer solar atmosphere (Heyvaerts and Priest, 1983).

We find that cubic wave damping dominates shear linear dissipation whenever the Alfvén wave velocity amplitude δv_y exceeds a few times ten metres per second. In the nonlinear regime, phase-mixed waves are marginally stable, while non-phase-mixed waves of wavenumber k_a are damped over a timescale k_uRe ₀|δ v_y/v_A |^?2, Re ₀ being the Reynolds number corresponding to the Braginskij viscosity coefficient η₀ and v_A the Alfvén speed. Dissipation is most effective where β = (v_s /v_A) ² ≈ 1, v_s being the speed of sound.     

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

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

井中声电转换波场辐射能量定量分析

井中声激发极化效应研究是实现声激发极化测井实用化的技术基础,实质上是声波在双相孔隙介质中向电磁波的转化,其结果就是形成井中声电转换波场,最终表现为井中声场辐射能向电磁场能的转化.通过对双相孔隙介质中声激发极化过程声电转换波场能量的理论分析,研究了井中声电转换波场能量的总体构成,提出了井中低频声激发极化电场电能体密度的概念,并推导出具体表达式.通过采用实际地层参数模型进行井中低频声激发极化电场的模拟计算得到了井中低频声激发极化过程中电能体密度与声波频率与地层参数的相关关系.     

低速地层裸眼井中弹性波的传播

本文在文献[4]的基础上分析了低速地层（地层横波速度低于井内流体速度）裸眼井中由点源激发的弹性波的传播。指出在这种情况下井内传播的波仅有共振纵波、斯通利波和衰减的反射波,共振横波和伪瑞利波将不再出现。由数值计算得出了相速度和群速度频散曲线,给出了势函数的振幅表达式和反射波随井轴的衰减规律。     

软地层横波速度的偶极随钻测井反演

在电缆测井中,可以利用偶极声场中的弯曲波反演软地层的横波速度.然而,在随钻声波测井(LWD)中,钢制钻铤的存在使得井孔结构变得复杂,同时改变了井孔声场,弯曲波也变得难以测定.此外,弯曲波与钻铤波耦合在一起,使得地层横波速度的反演变得困难.本文计算分析了随钻声场的频散曲线和激发曲线,注意到了偶极舒尔特波在较宽频带内速度频散很弱,特别是在本文研究的软地层情况下,偶极舒尔特波速度在3至25kHz的频率范围内几乎为一个常值,并且该值与地层横波速度存在一一对应关系.舒尔特波速度远小于其他模式波速度,与其他模式波在时域上易分离.相对于其他地层参数,舒尔特波对地层横波速度十分敏感.因此,它可以用来反演软地层的横波速度.     

The relationships between the velocities,attenuations and petrophysical properties of reservoir sedimentary rocks1

We have measured the velocities and attenuations of compressional and shear waves in 29 water-saturated samples of sandstones and shales at a confining pressure of 60 MPa and at frequencies of about 0.85 MHz. The measurements were made using a pulse echo method in which the samples (diameter 5 cm, length 1.5 cm to 2.5 cm) were placed between perspex buffer rods inside a high-pressure cell. The velocity of each seismic wave was determined from the traveltime difference of equivalent phase points (corrected for diffraction effects) of the signals reflected from the top and from the base of each sample. Attenuation was determined in a similar way by comparison of the diffraction corrected amplitudes of the signals. The attenuation data are presented as ‘quality factors’: Q_p and Q_s for compressional and shear waves respectively. The results show that Q_s is strongly correlated with V_s, that Q_p is weakly correlated with V_p, and that Q_p is strongly correlated with Q_s. Q_p is strongly dependent on the volume percentage of the assemblage of intra-pore minerals, whether they are clays or carbonates. It is concluded that the attenuation mechanism is due to the local fluid flow arising from the differential dilation of the solid rock frame and the intra-pore mineral assemblage, which is a result of their very different elastic moduli.     

Analysis of a guided wave along a conducting structure in a borehole

A guided wave along a borehole is often observed in borehole radar measurements. These guided waves deform the antenna pattern and can cause artefacts in radar measurements. A water-filled borehole or a conducting logging cable can function as a waveguide for electromagnetic waves under some conditions. We describe the theoretical characteristics of such a guided wave in a borehole and compare them with our experiments. The measured signal discussed was obtained with a directional borehole radar. This radar uses a cylindrical conformal array antenna as receiver and is a model of a conducting structure in a borehole. The induced field around the borehole was compared with the theory. The most fundamental symmetrical and asymmetrical modes were TM₀₁ and HE₁₁, and they were identified in the measured signals using time–frequency distribution analysis and by observation of the azimuthal field distribution of the magnetic field.     

含孔隙、裂隙致密介质中多极子声波的传播特征

非常规油气藏（如致密性地层及蕴藏油气的页岩地层）的重要特征是低孔、低渗,但裂隙或裂缝比较发育.为满足非常规勘探的需求,本文将孔、裂隙介质弹性波传播理论应用于多极子声波测井的井孔声场模拟,重点研究了致密介质中裂隙发育时多极子声波的传播机理以及衰减特征.井孔声场的数值计算结果表明裂隙的存在明显改变了弹性波和井孔模式波的频散、衰减和激发强度,尤其是井壁临界折射纵波的激发谱的峰值随着频率的增加逐渐降低,这与应用经典的Biot理论下的计算结果相反,且裂隙的存在也使得饱含水和饱含气时临界折射纵波激发强度的差异变大.井孔模式波的衰减与地层横波衰减和井壁流体交换有关,井壁开孔边界下致密地层裂隙发育还使得井孔斯通利波和艾里相附近的弯曲波对孔隙流体的敏感性增强,在井壁闭孔边界条件下引起井孔模式波衰减的主要因素是裂隙引起的地层横波衰减造成的,且在截止频率附近弯曲波的衰减与地层的横波衰减一致.数值计算结果为解释非常规油气地层的声学响应特征提供了参考.     

声电效应测井电声比及其与地层渗透率的关系

基于流体饱和孔隙介质中声波-电磁场耦合效应的测井方法具有潜在的应用价值.本文从Pride动电耦合波方程组入手,推导了伴随斯通利波的井孔电场与声压比值(电声比)的频率域表达式.结果表明,在较低频率条件下,电声比幅角的正切值与渗透率呈反比.在此基础上,提出了利用低频的声电效应测井全波反演地层渗透率的方法.针对砂岩地层,从计...