射线法在粘—弹性介質中的推广

一、引言在完全弹性介貭中,弹性波的传播問題已經研究得比較深入,尤其用射綫法来解决波的传播強度問題,有了一套完善的方法。实际上,大地岩层近于粘—弹性介貭,地震波在传播过程中要受到內摩擦的作用,振幅随距离逐漸衰減。衰減規律,現在均根据經驗公式来确定。在理論上,虽然可以通过解波动方程的方法解决这問题,但所得結果,即使极簡单的問題也是十分复杂的。在完全弹性介貭中,也存在同样的情形,但如应用射线法,則变得較为方便。     

波在飽水孔隙弹性介貭中的传播

一、引言弹性体中波的传播問題迄今已有較广泛的研究.对飽水孔隙弹性介貭中波的传播問題研究的尚少。但此种問題在地球物理、地震、土木工程和声学等研究中均有重要意义.     

波在飽水孔隙弹性介貭中的传播

一、引言弹性体中波的传播問題迄今已有較广泛的研究.对飽水孔隙弹性介貭中波的传播問題研究的尚少。但此种問題在地球物理、地震、土木工程和声学等研究中均有重要意义.     

弹性波反射法在桩基检测中的理论与应用

根据波动理论探讨了弹性波在砼桩中的传播衰减规律及反射特征，提出了判断桩基质量的有关理论依据，并结合工程实例应用弹性波反射法讨论了时域、频域分析法的基本原理及存在问题．     

瞬态弹性波散射的变换域积分方程法

针对瞬态弹性波散射的问题,从弹性动力学问题的积分表示定理出发,采用Laplace变换的方法,得到了变换域内均质体位移场的积分方程表示;在此基础上推导了适合瞬态弹性波对异质体散射求解的变换域位移场积分方程。     

弹性波反射法在探测软弱层中的应用

相似文献   

有限单元法弹性波偏移

本文在弹性动力学问题的有限单元解的基础上,提出了一种弹性波偏移的新方法——有限单元法弹性波偏移。本方法不仅能使反射纵波得到偏移,而且能同时对反射纵波和反射横波进行偏移。 用本方法已对一种层状模型成功地进行了检验。说明了它对反射纵波和转换横波同时向下延拓和偏移是稳定的,所得到的偏移结果是令人满意的。     

有限单元法弹性波偏移

本文在弹性动力学问题的有限单元解的基础上,提出了一种弹性波偏移的新方法--有限单元法弹性波偏移。本方法不仅能使反射纵波得到偏移,而且能同时对反射纵波和反射横波进行偏移。 用本方法已对一种层状模型成功地进行了检验。说明了它对反射纵波和转换横波同时向下延拓和偏移是稳定的,所得到的偏移结果是令人满意的。     

伪谱法在常Q粘弹介质地震波场模拟中的应用效果

本文采用伪谱法推导了常Q粘弹介质中地震波场的运动方程和波动方程,并进行了地震波场数值正演模拟,模拟结果达到了其它模拟方法同样的效果.但与其它方法相比,采用伪谱法不仅简洁明了,更便于理解;更重要的是它比以往的方法更具有通用性.由此可以说,如果采用伪谱法模拟其它更复杂的粘弹介质模型时,完全可以避免不必要的而且烦杂的数学推导过程;直接将体积松弛模量和剪切松弛模量代入即可得其波动方程.     

粘弹性介质中地震体波Q值的测定方法

本文以常用的几种粘弹性体模拟地球介质的粘弹性,建立了该种介质中的波动方程组;在运用合理的近似处理方法基础上,得到了基本解;提出了校正其它衰减因素的功率谱切比雪夫拟合方法,由此给出了粘性Q值的概念和算法,最后讨论了该方法的可行性。     

AN EXTENSION OF THE BORN INVERSION METHOD TO A DEPTH DEPENDENT REFERENCE PROFILE*

An extension of the multidimensional Born inversion technique for acoustic waves is described. In earlier work, a perturbation in reference sound velocity was determined by assuming that the reference velocity was constant. In this extension, we allow the reference velocity to be a function of the depth variable z. The output of this method is a high-frequency bandlimited reflectivity function of the subsurface. The reflectivity function is an array of bandlimited singular functions scaled by the normal reflection strength. Each singular function is a Dirac delta function of a scalar argument which measures distance normal to a reflecting interface. Thus, the reflectivity function is an indicator map of subsurface reflectors equivalent to the map produced by migration. In addition to the assumption of small perturbation, the method requires that the reflection data reside in the high frequency regime in a well-defined sense. The method is based on the derivation of an integral equation for the perturbation in sound velocity from a known reference velocity. When the reference velocity is constant, the integral equation admits an analytic solution as a multifold integral of the reflection data. Further high frequency asymptotic analysis simplifies this integral considerably and leads to an extremely efficient numerical algorithm for computing the reflectivity function. The development of a computer code to implement this constant-reference-velocity solution is published elsewhere. For a reference velocity c(z) we can no longer invert the integral equation exactly. However, we can write down an asymptotic high-frequency approximation for the kernal of the integral equation and an asymptotic solution for the perturbation. The computer implementation of this result is designed along the same lines as the code for constant background velocity. In tests the total processing time for this algorithm with depth-dependent background velocity is usually considerably less than that required by a standard Kirchhoff migration algorithm. The method is implemented as a migration technique and compared with alternative migration algorithms on the flanks of the salt dome.     

计算横向非均匀介质体波理论地震图的Maslov渐近法

本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。     

计算横向非均匀介质体波理论地震图的Maslov渐近法

本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。     

计算横向非均匀介质体波理论地震图的Maslov渐近法

本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。     

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.     

一种反演不均匀介质大地电磁测深资料的方法

通过反演受近地表不均匀介质影响的大地电磁测深资料的新方法,可获得较准确的一维地电断面参数、不均匀体埋藏深度及其电性特征。本方法计算速度快,单测ρxy,ρyx可独立反演并相互验证结果的准确性     

THE INSENSITIVITY OF REFLECTED SH WAVES TO ANISOTROPY IN AN UNDERLYING LAYERED MEDIUM1

Propagation in the plane of mirror symmetry of a monoclinic medium, with displacement normal to the plane, is the most general circumstance in anisotropic media for which pure shear-wave propagation can occur at all angles. Because the pure shear mode is uncoupled from the other two modes, its slowness surface in the plane is an ellipse. When the mirror symmetry plane is vertical the pure shear waves in this plane are SH waves and the elliptical SH sheet of the slowness surface is, in general, tilted with respect to the vertical axis. Consider a half-space of such a monoclinic medium, called medium M, overlain by a half-space of isotropic medium I with plane SH waves incident on medium M propagating in the vertical symmetry plane of M. Contrary to the appearance of a lack of symmetry about the vertical axis due to the tilt of the SH-wave slowness ellipse, the reflection and transmission coefficients are symmetrical functions of the angle of incidence, and further, there exists an isotropic medium E with uniquely determined density and shear speed which gives exactly the same reflection and transmission coefficients underlying medium J as does monoclinic medium M. This means that the underlying monoclinic medium M can be replaced by isotropic medium E without changing the reflection and transmission coefficients for all values of the angle of incidence. Thus no set of SH seismic experiments performed in the isotropic medium in the symmetry plane of the underlying half-space can reveal anything about the monoclinic anisotropy of that underlying half-space. Moreover, even when the underlying monoclinic half-space is stratified, there exists a stratified isotropic half-space that gives the identical reflection coefficient as the stratified monoclinic half-space for all angles of incidence and all frequencies.     

ON THE ROLE OF PARTIAL RAY EXPANSION IN THE COMPUTATION OF RAY SYNTHETIC SEISMOGRAMS FOR LAYERED STRUCTURES1

This paper presents results of testing an efficient ray generation scheme needed whenever ray synthetic seismograms are to be computed for layered models with more than 10‘ thick’layers. Our ray generation algorithm is based on the concept of kinematically equivalent waves (the kinematic analogs) having identical traveltimes along different ray-paths between the source and the receiver, both located on the surface of the model. These waves, existing in any medium composed of laterally homogeneous parallel layers, interfere at any location along the recording surface, thereby producing a composite wavelet whose amplitude and shape depend directly on the number of kinematic analogs (the multiplicity factor). Hence, explicit knowledge of the multiplicity factor is crucial for any analysis based on the amplitude and shape of individual wavelets, such as wavelet shaping, Q estimation, or linearized wavelet inversion. For unconverted waves, such as those discussed in this paper, the multiplicity factor can be computed analytically using formulae given in the Appendix; for converted waves, the multiplicity factor should be computed numerically, using the algorithm employed for the computation of the seismograms presented in a previous paper by one of the authors.