首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 687 毫秒
1.
对波达波夫和Pride震电波方程组的对比分析   总被引:2,自引:1,他引:2       下载免费PDF全文
用Biot介质参数说明了波达波夫震电波方程组中弹性动力学 参数的含义,解释了第一类和第二类震电效应的意义,在忽略第一类震电效应条件下将该方 程组与Pride方程组进行比较,说明了二者在描述第二类震电效应方面的异同点. 同时指出 :波达波夫方程组忽略了流体与固体的耦合质量;方程中的黏性耗散项丢掉了一个孔隙度因 子,依据该方程组计算出的弹性波和转换电场的幅度将偏大;边界条件之一存在错误,会影 响对波在界面上的反射透射规律的描述.  相似文献   

2.
The coupling–ray–theory tensor Green function for electromagnetic waves or elastic S waves is frequency dependent, and is usually calculated for many frequencies. This frequency dependence represents no problem in calculating the Green function, but may pose a significant challenge in storing the Green function at the nodes of dense grids, typical for applications such as the Born approximation or non–linear source determination. Storing the Green function at the nodes of dense grids for too many frequencies may be impractical or even unrealistic. We have already proposed the approximation of the coupling–ray–theory tensor Green function, in the vicinity of a given prevailing frequency, by two coupling–ray–theory dyadic Green functions described by their coupling–ray–theory travel times and their coupling–ray–theory amplitudes. The above mentioned prevailing–frequency approximation of the coupling ray theory enables us to interpolate the coupling–ray–theory dyadic Green functions within ray cells, and to calculate them at the nodes of dense grids. For the interpolation within ray cells, we need to separate the pairs of prevailing–frequency coupling–ray–theory dyadic Green functions so that both the first Green function and the second Green function are continuous along rays and within ray cells. We describe the current progress in this field and outline the basic algorithms. The proposed method is equally applicable to both electromagnetic waves and elastic S waves. We demonstrate the preliminary numerical results using the coupling–ray–theory travel times of elastic S waves.  相似文献   

3.
再论地震数据偏移成像   总被引:1,自引:1,他引:0       下载免费PDF全文
利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.  相似文献   

4.
双变参数标量纵波方程正演模拟方法   总被引:1,自引:0,他引:1       下载免费PDF全文
常见弹性波动理论的建立是基于介质均匀这一基本假设,实际介质的非均匀性非常普遍.为研究连续介质中波的传播特征,本文从弹性力学中建立弹性波动方程的三个基本方程出发,考虑连续介质弹性参数的空变特征,建立非均匀介质的弹性波动方程,利用Alkhalifah声学近似思想建立位移表征的纵波波动方程,利用本征值问题求解方法建立标量波频率-波数域传播算子,从而建立描述纵波传播的标量波方程,其中波函数为纵波位移的散度,不同于均匀介质标量波方程的波函数为位移势.随后推导含PML边界波动方程差分格式并建立不同模型数值模拟进行数值试算,与均匀假设标量波方程和变密度方程对比证明本方法的准确性和稳定性.  相似文献   

5.
When there exists anisotropy in underground media, elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media, four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.  相似文献   

6.
A general solution is deduced of the differential equations describing the propagation of elastic waves in a dissipative liquid-filled viscoelastic porous solid. The velocities of three existing waves have been expressed in convenient form using the moduli of the solid phase and by introducing the frequency-dependent equivalent mass densities. The solution is then used to examine some of the phenomena which arise when each of the three-body waves, in turn, are incident on a traction-free plane boundary. Analytic expressions for the reflection coefficients are obtained. Numerical calculations have been made, for a particular model, in case of incidentP I wave. Effect of viscoelasticity and viscosity on the reflection coefficients has also been exhibited.  相似文献   

7.
基于GSLS模型TI介质衰减拟声波方程   总被引:2,自引:2,他引:0       下载免费PDF全文
随着计算机硬件技术的发展以及高分辨率勘探需求的增加,我们希望能够更准确地模拟地下介质,得到更丰富的地层信息.然而,传统的声学假设并不能描述实际地层所存在各向异性和黏滞性,使得成像分辨率较低.为了实现深部储层的高精度成像,本文同时考虑了介质的各向异性和黏滞性,从TI介质弹性波的基本理论出发,结合各向异性GSLS理论,并通过声学近似方法导出基于GSLS模型的各向异性衰减拟声波方程.数值模拟表明该方程既能准确地描述各向异性介质下的准P波运动学规律,又能体现地层的吸收衰减效应;模型逆时偏移结果表明,在实现成像过程中考虑各向异性和黏滞性的影响,能对高陡构造清晰成像,且剖面振幅相对均衡,分辨率较高.  相似文献   

8.
孔隙地层中存在着动电耦合现象,弹性波作用在声阻抗或者电化学性质差异地层会产生辐射电磁波.本文探讨了随钻条件下,利用孔隙地层中的动电耦合效应探测井旁地质体的可行性.为此模拟了井外存在声阻抗差异和电化学差异界面时声波诱导的辐射电场响应特征.求出柱坐标下介质中声场和电场函数表达式,由此计算并比较了均匀介质和弹性差异界面情况下井中接收到的电场.发现随钻条件下,井外的弹性差异界面会产生以光速传播的辐射电磁波.该电磁波反映了井外地质信息,并先于井中直达声波的伴随电磁场到达接收器,因而可以直接观测到.此外,还发现界面辐射电磁波可以用来探测油、水等电化学性质差异大而弹性阻抗差异小的界面;声源频率以及界面两侧地层性质差异程度是辐射电磁波幅度的影响因素.  相似文献   

9.
孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen(LAC)广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.  相似文献   

10.
越来越多的研究表明来自基岩的地震波并不是垂直地面向上传播的。地震波在斜入射与垂直入射时所产生的场地效应有很大不同,由于存在全反射现象,SV波在斜入射时产生的场地效应更为复杂。文章基于均匀弹性半空间地震波传播理论,分别推导得到SV波入射角在小于、大于等于临界角时地震动的计算表达式,通过模型算例研究SV波全反射作用下的地震动特性。研究发现:由SV波产生的地震动主要由入射波和反射波构成,滑行波作用可以忽略;地面运动轨迹具有面波旋转振动特点;随着入射角的增大,地面震动从以水平方向振动为主逐渐过渡到以垂直方向振动为主,两个方向的振动均小于入射波峰值的2倍。  相似文献   

11.
针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题。以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型。通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱和多孔热弹性介质边界上的四种反射波的振幅反射率的理论表达式。在此基础上进行数值计算,分别讨论平面S波的入射频率、入射角和热膨胀系数等参数对四种反射波的振幅反射率的影响情况。结果表明:各反射波的振幅反射率分别随频率和热膨胀系数的增大而增大,同时也受到平面S波入射角变化的影响。该结论对于土动力学的理论研究及其相关的工程勘探具有一定的指导意义。  相似文献   

12.
本文作为基础理论研究, 从热弹性波方程组出发, 求出了热弹性波在层状介质中传播的解析解, 并给出了热弹性波对于平面夹层的反射波, 透射波及其伴随的具有相同传播速度的温度波。此时, 反射系数、 透射系数均为复数表明, 反射波、 透射波的振幅及相位都受介质的物性参数、 夹层厚度及入射波频率的影响。这些与弹性波不同的结果会在地震勘探、 地震工程等有关问题的研究中得到应用。  相似文献   

13.
各向异性介质qP波传播描述I:伪纯模式波动方程   总被引:6,自引:6,他引:0       下载免费PDF全文
程玖兵  康玮  王腾飞 《地球物理学报》2013,56(10):3474-3486
地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征.各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用.为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程.本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程.文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势.  相似文献   

14.
应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法,将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质,给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波,数值结果表明具有很高的计算精度。  相似文献   

15.
The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.  相似文献   

16.
Expansion formulae associated with the interaction of oblique surface gravity waves with a floating flexible plate in the presence of a submerged horizontal flexible structure are derived using Green’s integral theorem in water of finite and infinite water depths. The associated Green’s functions are derived using the fundamental solution associated with the reduced wave equation. The integral forms of the Green’s functions and the velocity potentials are advantageous over the eigenfunction expansion method in situation when the roots of the dispersion relation coalesce. As an application of the expansion formulae, diffraction of oblique waves by a finite floating elastic plate in the presence of a submerged horizontal flexible membrane is investigated in water of finite depth. The accuracy of the numerical computation is demonstrated by analysing the convergence of the complex amplitude of the reflected waves and the energy relation. Effect of the submerged membrane on the diffraction of surface waves is studied by analysing the reflection and transmission coefficients for various parametric values. Further, the derivation of long wave equation under shallow water approximation is derived in a direct manner in the appendix. The concept and methodology can be easily extended to deal with acoustic wave interaction with flexible structures and related problems of mathematical physics and engineering.  相似文献   

17.
With the multi-wave, multi-component seismic wave exploration, one can apply the anisotropy of fracture media to analyze the attributes of the fracture media, including the fracture azimuth. In the meantime, the techniques of full-polarimetric electromagnetic wave, including full-polarimetric borehole radar, can also be used to analyze the attributes of the fracture. However, the analysis precision of both the multi-component elastic wave exploration and full-polarimetric electromagnetic wave exploration is prone to the influence of noise and other factors. So far, some researchers have conducted studies on the joint inversion of electromagnetic waves and seismic waves. This paper develops evaluation techniques of fracture azimuth by electromagnetic wave, elastic wave, and joint analysis of coincident elastic reflection and electromagnetic data. Firstly, based on the shear wave splitting of elastic waves, this paper develops a statistical analysis technique which applies Pearson correlation coefficient to count and analyze the azimuth angle of fracture. Secondly, based on the information of electromagnetic polarization rotated by fracture, this paper develops a statistical analysis method of full-polarimetric electromagnetic waves which applies the maximum amplitude ratio between the co-polarization and cross-polarization to analyze the azimuth angle of fracture. Furthermore, based on the analysis result of the elastic wave and full-polarimetric electromagnetic wave, this paper develops a joint analysis technique which adopts the standard deviation. At last, authors in this study conduct joint detection experiments on the coincident fracture medium by using the ultrasonic and full-polarimetric ground penetrating radar. The experimental result indicates that both single geophysical methods are capable of analyzing the fracture azimuth angle, but the joint analysis is more accurate.  相似文献   

18.
A problem of reflection and transmission of elastic waves at a plane interface between a uniform elastic solid half-space and a porous elastic half-space containing two immiscible fluids is investigated. The theory developed by Lo, Sposito and Majer for porous media containing two immiscible fluids is employed to find out the reflection and transmission coefficients. The incident wave is assumed to propagate through the uniform elastic half-space and two cases are considered. In the first case, a beam of plane longitudinal wave is assumed to be incident and in the second case, a beam of transverse wave is assumed to be incident at the interface. By taking granite as impervious elastic medium and columbia fine sandy loam containing air-water mixture as porous medium, reflection and transmission coefficients are obtained. By neglecting the inertial coupling coefficients, these coefficients are reduced to those obtained by Tomar and Arora using the theory of Tuncay and Corapcioglu. It is found that the inertial coupling parameters significantly affect the phase speeds and the amplitude ratios of the transmitted waves.  相似文献   

19.
从双相横向各向同性介质弹性波波动方程出发 ,利用边界上的 4个连续性条件 ,计算双相横向各向同性介质分界面上弹性波反射和透射系数 .计算表明 ,快纵波在双相横向各向同性介质分界面上 ,要产生反射快纵波、反射转换 qSV波、反射转换慢纵波和透射快纵波、透射转换qSV波、透射转换慢纵波 .反射转换慢纵波振幅和透射转换慢纵波振幅均较小 .频率、耗散和各向异性大小影响着快纵波反射系数的大小 .  相似文献   

20.
Chen's technique of computing synthetic seismograms,which decomposes every vector with a set of basis of orthogonality and completeness before applying the Luco-Apsel-Chen(LAC)generalized reflection and transmission coefficients method,is confirmed to be efficient in dealing with elastic waves in multi-layered media and accurate in any frequency range.In this article,we extend Chen's technique to the computation of coupled seismic and electromagnetic(EM)waves in layered porous media.Expanding the involved mechanical and electromagnetic fields by a set of scalar and vector wave-function basis,we obtain the fundamental equations which are subsequently solved by using a recently developed version of the LAC generalized reflection and transmission coefficients method.Our approach and corresponding program is validated by reciprocity tests.We also show a numerical example of a two-layer model with an explosion source.The P-to-EM conversion waves radiated from the interface may have potential application.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号