层状介质的声波波动方程反演

基于广义反射透射系数矩阵正演方法 ,讨论了层状介质的声波波动方程反问题 .推导出波数频率域中的雅可比矩阵的解析表达式 ,其计算在正演过程中求出 .采用最小二方法可得到层介质参数 .数值结果表明反演方法的正确有效性 .     

海底强地面运动计算

本文将计算综合地震图的广义反射、透射系数矩阵和离散波数方法进行了推广,使之适合于计算海底地面运动的位错点源格林函数。通过对数值例子的分析表明,理论模拟方法可以用于研究海水层和沉积层性质对海底运动的影响,以根据陆地上的资料来定性地获得海底强地面运动的信息。     

计算分层介质中位错点源静态位移场的广义反射、透射系数矩阵和离散波数方法

本文将计算合成地震图的广义反射、透射系数矩阵和离散波数方法应用于静态问题,给出了分层介质中位错点源产生的静态位移场的计算方法。该方法保留了反射、透射矩阵方法的原有优点,通过将数值结果与解析解的对比表明,该方法具有较高的数值精度,可用于研究地震断层活动造成的地表形变。     

计算分层介质中位错点源静态位移场的广义反射、透射系数矩阵和离散波数方法

本文将计算合成地震图的广义反射、透射系数矩阵和离散波数方法应用于静态问题,给出了分层介质中位错点源产生的静态位移场的计算方法。该方法保留了反射、透射矩阵方法的原有优点,通过将数值结果与解析解的对比表明,该方法具有较高的数值精度,可用于研究地震断层活动造成的地表形变。     

计算综合地震图的广义反射、透射系数矩阵和离散波数方法（二）——对不同深度点源的算法

本文讨论了计算近场综合地震图时使用的广义反射、透射系数矩阵和离散波数方法,并对不同深度点源的情况给出了新的算法及相应的计算实例。用本文提供的新算法,与原算法相比,可以减少60％的计算时间。在需要计算大量点源格林函数时,该方法是很有实用价值的。     

广义反射-透射系数算法的无量纲化

水平层状介质中瑞利波频散曲线的计算一直是受关注的问题.陈晓非提出的广义反射-透射系数方法虽然具备很好的精确性和稳定性,但是算法中起决定作用的矩阵E存在有量纲元素,并且不同元素之间的数量级差异很大.本文作者通过把广义反射-透射系数算法无量纲化,得到了简洁的瑞利波频散函数的公式体系.分别计算了改进后算法和原算法中矩阵E的最大模元素与最小模元素的比值ε,发现前者的矩阵E中最大模元素与最小模元素最多只相差1个数量级,而后者的矩阵E中最大模元素与最小模元素之间最大达11个数量级的差异,尽管后者可以通过选取变量的单位减小ε的值,但新算法精度更高,从而完善了利用广义反射-透射系数方法求解瑞利波频散曲线问题的理论和算法.     

充液井孔中的相控线阵声波辐射器在地层中产生的声场

应用有限元法研究了充液井孔中的相控线阵声波辐射器在井外地层中产生的声场,探讨了增强声源向地层的某一特定方向辐射的声波能量的方法. 采用相控线阵声波辐射器并调节其参数可以控制由充液井孔向地层中辐射的声波能量的方向,使声源向地层中的某个特定方向辐射的纵波能量得到明显提高. 在井间地震勘探和反射声波成像测井等领域中应用相控线阵声波辐射器,将有利于增大探测深度、提高探测分辨率和接收信号的信噪比.     

波在多层弹性介质中一些特性的研究

本文根据弹性力学理论计算出波在多层弹性介质上的反射和透射系数,给出计算反射和透射系数的源程序和框图,并对波在多层弹性介质中的特性进行了一些研究。     

基于快速广义反射透射系数方法的面波群速度计算

在利用实际地震资料反演地球内部结构时,群速度往往比相速度更有意义.本文在水平层状均匀介质模型中,用快速广义反射透射系数方法（Pei,2008）计算了基阶瑞雷面波的群速度,并推导了Chen(1993) 4×4 层矩阵E的逆的解析表达式.文中用数值实验的方法合成了两个典型模型的群速度频散曲线,并与复合矩阵方法计算的结果进行了对比,结果表明本文的方法是切实可行的,与广义反射透射系数方法相比明显地提高了计算效率,尤其是在高频、层厚度大和层数多的情况下表现得更为突出.     

用滑移界面理论模拟多裂缝体系的偶极横波远探测声场

偶极横波远探测已成为探测储层裂缝的重要测井技术,裂缝及其充填物质的弹性性质决定了其对远探测声场的响应特征.本文运用滑移界面理论来描述井外地层中裂缝对声场的响应特征,考察了SH横波的反射和透射系数随充填物的变化规律.建立了井外双裂缝远探测声场的几何模型,以双裂缝为例给出了偶极横波远探测声场的快速模拟方法,并将所得结论应用于现场远探测裂缝的识别和解释.本文的研究结果为偶极横波远探测数据的处理和解释提供了理论依据.     

用广义射线和有限差分计算近场理论地震图的混合方法

本文提出一种结合广义射线和有限差分理论而形成的混合方法（以下简称混合方法）．在介质均匀区（包括平行介质分层区）应用广义射线理论，在介质不均匀区应用有限差分理论．在这种结合过程中，把均匀区广义射线的结果作为不均匀区有限差分区波场连续的输入．对于均匀介质模型，混合方法理论地震图和广义射线理论地震图、广义反射透射系数矩阵及离散波数法地震图的对比，取得了满意的结果，在数值上证明了该方法的正确性．计算了１９９０年８月５日Ｈａｋｏｎｅ地震前震在日本Ａｓｈｉｇａｒａ山谷盲测点的响应，理论地震图和实际资料的比较同样取得好的结果．这种方法可用于地震勘探、工程地震、强地面运动的研究．     

利用广义反射-透射系数方法求解含低速层水平层状介质模型中面波频散曲线问题

在计算面波频散问题的广义反射-透射系数方法基础上研究了含有低速层的均匀平层固体地球模型的面波的频散曲线以及简正振型的计算问题. 当不存在低速层时,仅利用陈晓非1993年提出的单一的久期函数对求解这两个问题已经足够. 但是当模型中存在低速层时,该单一久期函数在数值求解较低阶模式很困难. 本文提出用久期函数族代替该单一久期函数、对不同的简正振型模式采用对应的久期函数来求解的方法,从而更加容易有效地求解较低阶的模式,完善了利用广义反射-透射系数方法求解面波频散问题的理论和算法.     

合成地震图的广义反射透射系数矩阵方法

简述了在均匀层状介质中，合成地震图的广义反射透射系数矩阵方法，并通过与有人的计算结果对比，检验了经过修改调试后的合成地震图的计算程序。计算与对比表明，经修改调试后的程序对于各种因素，包括震源类型、震源深度、震中、频率和慢度都有较好的普适性，适用于实际的地震波形研究。     

薄圆管压电换能器瞬态响应及其在井中激发声场

针对声波测井中采用的柱状薄圆管压电换能器,提出一种研究换能器响应的解析方法.该方法考虑了轴向与径向振动之间的耦合作用,导出置于耦合液体中压电换能器径向位移满足的微分方程和频域传递函数,比较分析了不同激励源作用下换能器的位移响应特性.利用波数域实轴积分法,计算了换能器在弹性地层井孔中激发的声场,研究了不同地层中的波形差异.与前人在假定声源强度的模型下只能求出波场的相对幅度不同,本文得出在一定激励电压作用下换能器在井孔中激发声场的幅度值.这为说明声波测井响应信号的强弱及仪器的声系设计提供了依据,对估计动电效应测井信号的强弱有参考意义.     

A note on the equivalence of three major propagator algorithms for computational stability and efficiency

It is shown in this note that the three methods, the orthonormalization method, the minor matrix method and the recursive reflection-transmission matrix method are closely related and solve the numerical instability in the original Thomson-Haskell propagator matrix method equally well. Another stable and efficient method based on the orthonormalization and the Langer block-diagonal decomposition is presented to calculate the response of a horizontal stratified model to a plane, spectral wave. It is a numerically robust Thomson-Haskell matrix method for high frequencies, large layer thicknesses and horizontal slownesses. The technique is applied to calculate reflection-transmission coefficients, body wave receiver functions and Rayleigh wave dispersion.     

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

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

TUNNELLING AND THE GENERALIZED RAY METHOD IN PIECEWISE HOMOGENEOUS MEDIA1

The contributions of the tunnelled constituents to a seismic wave are analysed in two different configurations pertaining to homogeneous acoustic media: a thin high-velocity layer, present in a plane-layered configuration, and a thin layer in media separated by dipping interfaces. The generalized ray method in the far-field is used to determine them. We expand around the relevant ray parameters in order to determine the characteristics of the tunnelling ray and find that the most important feature of this type of ray is a phase (in terms of asymptotic ray theory) which has a real and an imaginary part. Numerical results illustrate this.     

井筒声场研究回顾与展望

声波测井是测井学科的三大方法之一,在储层评价中起到了重要作用．随着陆上油气勘探目标趋向孔隙结构复杂、非均质性强的油气储集层,基于地层各向同性介质的测井理论显露弊端;因此研究复杂地质条件下的井筒声场分布特征、为储层解释和发展测井新方法提供理论依据成为必然．本文对井筒声场研究进行了回顾,并就当前井筒声场研究现状进行了总结,对其发展趋势和应用前景进行了分析。     

Resolution of prestack depth migration

The resolution of a general 3-D common-shot elastic prestack depth migration in a heterogeneous anisotropic medium is studied approximately, using the ray theory. It is demonstrated that the migrated section can approximately be represented by the convolution of the reflectivity function with the corresponding local resolution function. Alternatively, it can also be approximately represented by the convolution of the spatial distribution of the weak-contrast displacement reflection-transmission coefficient with the corresponding local resolution function. The derived explicit approximate equations enable us to predict the migration resolution approximately without doing the whole and expensive migration. The equations are applicable to 3-D elastic migrations in 3-D isotropic or anisotropic, heterogeneous velocity models. Both the reflectivity function and the spatial distribution of the weak-contrast displacement reflection-transmission coefficient approximately determine the linear combination of the perturbations of elastic moduli and density to which the migrated section is sensitive. The imaged linear combination of the perturbations of elastic parameters depends on the selection of the polarizations (wave types) of the incident and back-propagated wavefields and on the directions of propagation. The resolution of the linear combination of the perturbations of elastic moduli and density in the migrated section is determined by the above mentioned local resolution functions. The local resolution functions depend on the aperture and on the imaging function. The imaging function is determined by the source time function and by the form of the imaging functional. The local resolution functions are considerably sensitive to heterogeneities. The local resolution functions in elastic media differ from their acoustic counterparts, especially by the existence of converted scattered waves in elastic media.     

Study on the simulation of acoustic logging measurements in horizontal and deviated wells

The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.