横波速度径向层析成像方法及应用研究

近来三维声波成像技术的一个重要进展是利用偶极子弯曲波反演地层径向横波速度剖面。确定地层径向横波速度变化对评价井壁稳定性、估算地应力和优化油藏生产等有着重要的意义。偶极横波测井激发的弯曲波低频探测较深,高频探测较浅,本文根据偶极弯曲波的频散特征,采用约束反演的方法求取地层径向横波速度变化。在实际测井中,经常会遇到偶极子弯曲波的频散曲线在高频时缺失或者不准确的情况,这会使约束反演方法产生误差。针对这一现象,本文提出利用理论频散数据进行高频约束,并通过数值模拟验证了改进后约束反演方法的正确性。本文给出了一个致密气层的横波速度径向剖面综合应用实例,发现地层径向横波速度变化程度与岩石脆性有很好的对应关系,最终的压裂效果由压前、压后正交偶极横波各向异性的变化得到了印证,压裂后油气产量提高。     

流体-孔隙介质圆柱界面波传播特性

研究流体-孔隙介质圆柱界面波传播特性,分析孔隙介质孔隙率等参数对频散曲线的影响。理论上建立了无限大流体包裹孔隙介质圆柱界面波的模型,利用孔隙介质弹性波动理论,通过数值模拟计算得到流体-孔隙介质圆柱的频散曲线及时域波形,并分析了孔隙介质为开孔和闭孔状态下孔隙介质圆柱半径、孔隙率及渗透率对频散曲线的影响。结果表明,时域上斯通利波可以被明显区分开,孔隙介质圆柱半径的变化改变了圆柱尺度,孔隙率的变化改变了孔隙介质的纵、横波波速,因此对于斯通利波频散曲线的影响较大。而渗透率的变化既不改变圆柱的尺度也不改变孔隙层的纵、横波速度,因此对斯通利波频散曲线影响较小。      

利用简正波频散特征映射方法进行海洋环境参数反演

浅海波导中传播的脉冲信号可以通过频散特征映射方法(warping变换)进行简正波的有效分离,这使得频散曲线的获取更加容易和准确.利用简正波频散特征映射方法分析2010年南海声层析实验中单个阵元接收的声学信号,有效分离出各号简正波,并进一步通过Wigner-Ville分布处理提取出波导中各号简正波的频散曲线.将提取的频散曲线与理论计算到达时差进行匹配,实现了利用单水听器信号的海洋环境参数反演.实验结果表明:将简正波分离与Wigner-Ville分布相结合可有效消除Wigner-Ville分布在估计多分量信号时的交叉干扰,提取的频散特征可很好地用于同时反演海水和海底声速剖面。为评价反演结果的有效性,对反演结果进行了后验概率分析,并将一小时内的海水声速剖面反演结果与温度链实测值作对比,二者符合较好.      

利用气枪声源数据的地声参数反演*

本文通过对东海采集的气枪声源数据进行处理,基于地声参数对不同声场观测量的敏感性差异,利用分步反演的策略进行分析。首先选择距离接收阵七个不同的距离点数据,通过warping变换准确提取第三阶至第八阶简正波的频散曲线;利用海底声速对频散曲线敏感的特性来反演声速;由Hamilton经验公式求得海底密度;通过传播损失拟合获得海底衰减。不同距离点数据反演的海底声速与密度一致性较高。实验提取的频散曲线和和反演参数仿真结果、实验获取的传播损失与反演参数获得传播损失均拟合较好。     

分层介质半空间瑞利波的时频分析

对分层介质半空间瑞利波的频散特性,用一种时频分析方法——重排的平滑伪魏格纳维尔分布(RSPWVD,Reassign- ment of Smooth Pseudo-Wigner Ville Distribution)进行了分析和研究。对均匀半空间和两层介质半空间的理论和实验研究表明,由于层状介质中瑞利波的频散曲线存在多个模式,所获得的群速度频散曲线在不同的频段显示出来的模式是对应表面位移幅度占主导作用的模式。频散曲线的这种模式判定对利用层状半空间的瑞利波反演介质参数是必须预先了解的。     

脱粘缺陷对粘接结构频散特性的影响

采用全局矩阵理论,通过引入内聚强度弱化模型和弹簧模型,计算了脱粘缺陷对粘接结构频散特性的影响。分析了内聚强度缺陷、界面缺陷及混合缺陷等因素对粘接结构频散特性的影响,发现频散曲线的演化规律与材料声学性质密切相关。不同材料和不同组合,其规律有很大不同。针对各向同性粘接结构,可能通过考察特定频率段内频散模态的变化,实现对脱粘缺陷的检测。     

利用井中偶极弯曲波反演渗透率的影响因素及灵敏度分析

本文以各向同性饱和流体孔隙介质充液井孔中声场传播理论为基础,系统考察了影响偶极弯曲波特征的主要因素,分析了快慢地层条件下弯曲波衰减和速度对不同参数的敏感度,对比结果表明,弯曲波的速度基本不受地层渗透率的影响,其衰减对渗透率的敏感性要远高于速度。弯曲波中心频率会随着波的衰减而发生偏移,说明利用弯曲波衰减和频移特征联合反演地层渗透率是可行的。通过对比快慢地层下的弯曲波衰减和速度灵敏度曲线,发现在慢速地层中反演渗透率的结果要好于快速地层。同时,本文也指出在低孔低渗储层中,还可以考虑利用弯曲波衰减识别气层。     

含孔隙分层半空间中瑞利波的传播特性研究

为了研究含孔隙介质分层半空间中瑞利波的传播规律,分析孔隙介质参数对瑞利波频散曲线的影响,本文进行了数值模拟研究。采用传递矩阵算法,计算了含孔隙分层半空间中一定频率范围内瑞利波所有模式的频散及激发强度曲线,并与均匀弹性固体分层半空间情况作了类比分析,在含孔隙覆盖层的两层模型和含低速孔隙夹层的三层模型下,详细研究了孔隙度、渗透率、层厚度等参数对瑞利波各模式的影响,发现孔隙度及层厚度的变化对频散曲线影响较大,而渗透率的变化对频散曲线影响较小。      

孔隙介质包裹的充液管道结构中导波传播特性

研究孔隙介质包裹的充液管道中纵向导波传播特性,分析孔隙介质参数对频散曲线的影响。建立了孔隙介质包裹充液管道的结构模型,利用孔隙介质弹性波动理论,建立对应的频散方程,数值模拟计算得到该模型的频散曲线和时域波形,并分析了孔隙介质参数以及管壁厚度对频散曲线的影响。结果表明孔隙介质的渗透率对于导波频散的影响较小,孔隙度的改变对时域波形的位移幅度影响较大。同时,导波存在衰减,且衰减随孔隙度增大而增大。所得结论为埋地管道无损检测方面提供一定理论参考。      

根据非线性贝叶斯理论的界面波频散曲线反演

通过时频分析法从海底环境噪声数据中提取界面波频散曲线,进而采用非线性贝叶斯反演方法估算海底沉积物厚度、剪切波速度、压缩波速度和密度等参数及其不确定性。参数的最大后验概率(MAP)估计值和边缘概率分布分别通过自适应单纯形模拟退火法和Metropolis-Hastings采样法在各参数先验区间内搜索获得,采用贝叶斯信息准则(BIC)从不同参数化模型中选择最优模型。界面波频散曲线反演结果表明:满足实测数据的最优海底模型结构为3层均匀分布剪切波速度剖面结构,海底深度的反演精度在800m以内,比起压缩波速度和密度,剪切波速度的不确定性更小,对界面波频散曲线更敏感。      

Stoneley wave generation by an external seismic point source in an infinite fluid-filled borehole embedded in a transversely isotropic formation: Basic relationships

The method of integral transformations is used to obtain a long-wave solution to the problem of tube wave excitation by an external point source in an infinite fluid-filled borehole embedded in a transversely isotropic formation. The external field that occurs in the formation gives rise to waves in the borehole fluid. The waves generated in the borehole include the lowest mode of the Stoneley wave (tube wave), which is the borehole eigenmode, and the qP-and qSV-waves. It is shown that the Stoneley wave is determined by the contributions of two poles in the complex plane of horizontal slowness. According to the asymptotic solution, the Stoneley wave can be described by one of three different waveforms depending on the relationship between the elastic parameters of the surrounding anisotropic formation and the borehole fluid. An analysis of the results of calculations shows that the shape and polarity of the Stoneley wave strongly depend on the sign of the nonellipticity parameter of the elastic medium, which offers a possibility of estimating the anisotropy of the borehole environment from observations of the waveform of the Stoneley wave.     

Tube-wave propagation in a fluid-filled borehole generated by a single point force applied to the surrounding formation

Propagation of tube waves in an infinite fluid-filled borehole, generated by a single-force point source placed in the elastic surrounding formation, is analyzed in the long-wave approximation. Integral representations of the precise solution are obtained both for fast and slow formations. An asymptotic analysis of tube-wave propagation in the fluid-filled borehole is performed on the basis of these two integral representations. The complete asymptotic wave field in the borehole fluid for a fast formation consists of P and SV phases and the lowest eigenmode of the Stoneley wave (tube wave). For a slow formation the conical Stoneley wave (Mach wave) is generated. It appears only behind the critical angle defined by the ratio of the S wave velocity in the formation to the low-frequency Stoneley wave velocity and decays weakly with an offset. Asymptotic wave forms are in good agreement with wave forms obtained by straightforward calculations.     

基底上覆分层薄膜超声Scholte波特性与薄膜定征

分层薄膜-基底结构广泛应用于微电子器件等诸多领域,但薄膜材料参数超声测量尤其是横波速度的定征是一个困难的问题。本文对液固界面Scholte界面波的频散特性和脉冲激励的声压响应进行了理论分析。结果表明,液固界面Scholte波频散与分层膜-基底结构的速度分布密切相关。薄膜材料各层的厚度和横波速度对界面波频散特性有显著影响。基于Scholte界面波的频散特性,提出了一种多层膜的多参数反演定征方法。首先针对理论信号进行薄膜参数反演,验证了该方法的可行性。后续对不同类型的多层膜材料样品进行了液固界面波激发与采集实验,实验信号的薄膜参数反演结果进一步验证了该方法的可行性和有效性。     

Inversion of guided-wave dispersion data with application to borehole acoustics

The problem of inferring unknown geometry and material parameters of a waveguide model from noisy samples of the associated modal dispersion curves is considered. In a significant reduction of the complexity of a common inversion methodology, the inner of two nested iterations is eliminated: The approach described does not employ explicit fitting of the data to computed dispersion curves. Instead, the unknown parameters are adjusted to minimize a cost function derived directly from the determinant of the boundary condition system matrix. This results in an efficient inversion scheme that, in the case of noise-free data, yields exact results. Multimode data can be simultaneously processed without extra complications. Furthermore, the inversion scheme can accommodate an arbitrary number of unknown parameters, provided that the data have sufficient sensitivity to these parameters. As an important application, we consider the sonic guidance condition for a fluid-filled borehole in an elastic, homogeneous, and isotropic rock formation for numerical forward and inverse dispersion analysis. We investigate numerically the parametric inversion with errors in the model parameters and the influence of bandwidth and noise, and examine the cases of multifrequency and multimode data, using simulated flexural and Stoneley dispersion data.     

Measurement of cylindrical Rayleigh surface waves using line-focused PVDF transducers and defocusing measurement method

Line-focused PVDF transducers and defocusing measurement method are applied in this work to determine the dispersion curve of the Rayleigh-like surface waves propagating along the circumferential direction of a solid cylinder. Conventional waveform processing method has been modified to cope with the non-linear relationship between phase angle of wave interference and defocusing distance induced by a cylindrically curved surface. A cross correlation method is proposed to accurately extract the cylindrical Rayleigh wave velocity from measured data. Experiments have been carried out on one stainless steel and one glass cylinders. The experimentally obtained dispersion curves are in very good agreement with their theoretical counterparts. Variation of cylindrical Rayleigh wave velocity due to the cylindrical curvature is quantitatively verified using this new method. Other potential applications of this measurement method for cylindrical samples will be addressed.     

国产的多极子阵列声波测井仪(MPAL)在油田开发生产中的应用

本文对国产的多极子阵列声波测井(MPAL)资料进行了分析,验证了仪器有较好的稳定性,测试结果表明:硬地层的纵波、横波、斯通利波波形清楚,幅度变化正常;通过其在岩性识别、裂缝识别、气层识别、岩石机械特性分析及各向异性的应用进行研究,能够较为准确地求取地层的纵波、横波、斯通利波的声波时差及幅度和衰减曲线,计算岩石力学参数,并进行地层各向异性评价。     

Perturbation Analysis on Guided Waves in a Fluid-Filled Borehole Surrounded by a Cubic Crystal Anisotropic Medium

The perturbation method is employed to analyse the guided waves in a borehole surrounded by a cubic crystal medium for the first time. The cubic crystal medium is regarded as a reference unperturbed isotropic state added to the perturbation. The dispersion characteristics of Stoneley wave, pseudo-Rayleigh wave, flexural wave, and screw wave are investigated in detail. It is found that dispersion of the guided waves excited by monopole and dipole sources does not depend on the azimuth of the source, whereas the dispersion of screw wave excited by quadrupole source is significantly related to the azimuth of the source. Screw waves propagated along different azimuth in the borehole can be split. This is different from screw waves in transversely isotropic media （hexagonal crystal）, which have been widely studied.     

Guided Waves in a Multi-Layered Cylindrical Elastic Solid Medium

We investigate the guided waves in a multi-layered cylindrical elastic solid medium. The dispersion function of guided waves is usually complex and the dispersion curves of all modes are not conveniently obtained. Here we present an effective method to obtain the dispersion curves of all modes. First, the dispersion function of the guided waves is transformed into a real function. The dispersion curves are then calculated for all the modes of the guided waves by the bisection method. The modes with the orders n = 0, 1, and 2 are analysed in two- and three-layer media. The existence condition of Stoneley wave is discussed. The modes of the guided waves are also investigated in a two-layer medium, in which the velocity of shear wave in the outer layer is less than that in the inner layer.