共查询到18条相似文献，搜索用时 312 毫秒
1.
合成三维横向非均匀介质远震体波接收函数的Maslov方法 总被引：5，自引：3，他引：2
三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要．本文发展了基于Ｍａｓｌｏｖ理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪,并采用递归算法组织射线的追踪过程,从而在保证追踪精度的同时提高了射线追踪的效率,并可完全避免传统射线追踪过程中的射线编码问题．正演计算结果表明,此法是一种高效、实用的三维横向非均匀介质接收函数计算方法． 相似文献
2.
3.
本文研究了横向非均匀吸收介质中的Maslov渐近理论地震图方法.对于如何提高计算精度及计算速度作了比较深入的探讨.本方法先计算远场格林函数,由褶积得到一般震源函数的弱吸收理论地震图.对于高Q引起的数值不稳定性,通过引入平滑函数得以解决.由此产生的新的被积函数具有非常好的收敛性,从而提高了计算精度及计算速度.以此理论编写的程序在模拟完全弹性介质以及一维吸收介质时与反射率法作了比较,并取得了令人满意的效果.对于二维横向非均匀吸收介质,计算方法与高斯光束方法作了对比,二者在主要特征上基本一致.数值对比中出现的一些差异也作了必要的讨论及说明. 相似文献
4.
Maslov渐近理论与辛几何算法 总被引：26，自引：9，他引：17
为了克服地震层析成像中的焦散问题,本文系统地研究了Ｍａｓｌｏｖ渐近理论与辛几何算法,同时提出了一种基于辛几何算法计算Ｍａｓｌｏｖ波场的数值计算方法,并就其中射线追踪这一重要环节,利用辛算法和非辛Ｒｕｎｇｅ－Ｋｕｔｔａ方法进行了数值模拟．计算结果表明,这两类算法在精度上相差无几,但辛算法的速度要快;在Ｈａｍｉｌｔｏｎ量保持方面,辛算法具有非辛Ｒｕｎｇｅ－Ｋｕｔｔａ方法无可比拟的优越性． 相似文献
5.
用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算． 相似文献
6.
本文提出了体波单震相测定震源机制的方法。按理论地震图公式，直达Ｐ波和Ｓ波地表垂直分量位移由地震矩，震源幅射和介质响应之积组成。我们将其变换成地震矩，假定断裂参数的单位地震矩单震相理论子波和待守断裂参数与假定断裂参数辐射函数之商的乘积，在各台理论与实际体波单震相最大振幅的最小二乘法中求解出震源参数，利用ＣＤＳＮ数字台网资料，我们分析了青海共和地震的震源机制。结果表明，１９９０年４月２６日Ｍｓ＝６．７ 相似文献
7.
波动方程Maslov射线解的有限差分计算 总被引：2，自引：0，他引：2
借助于守恒型程函方程的有限差分解法，构造出一种地震波场Ｍａｓｌｏｖ渐近解的数值计算方法，并针对不同坐标系进行了讨论。文中通过不同模型分别计算了透射波和反射波的射线路径、走时及相应的合成地震记录，结果表明本算法具有快速、准确的特点。本文算法可用于叠前Ｋｉｒｃｈｈｏｆｆ深度偏移及非均匀介质逆散射问题中Ｇｒｅｅｎ函数的计算，并就算法本身的局限性作了讨论。 相似文献
8.
早期的褶积微分算子都是基于正反傅立叶变换而实现的,其精度比四阶有限差分的精度稍高,本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与16阶有限差分的精度相当,远优于错格伪谱法的精确度.另外,2.5维数值模拟比二维模拟可以更真实地模拟三维介质的臬个剖面的波场,并且2.5维地震波模拟的计算量比三维模拟的计算量及计算耗时要大大减少.本文利用基于Forsyte广义正交多项式褶积微分算子法计算2.5维非均匀介质地震波场,模拟结果表明,该算法的计算速度快,计算精度高,能够直观、高效地反映复杂介质中波场的传播规律,并且2.5维波场数值模拟具有更高的计算效率,是一种非常值得深入研究并广泛应用的方法. 相似文献
9.
除介质各向异性之外,地球内部介质的横向非均匀性也是控制面波速度变化的重要因素.本文基于振型耦合和多重散射的地震波传播理论,数值模拟并分析了在具有均匀介质背景的三维异常体——上地幔横向非均匀介质中传播时,地震面波的振幅与偏振等动力学响应参数;其中分别模拟了不同周期入射、不同角度入射和不同尺度非均匀介质模型等多种情形下面波波场,并对横向非均匀性诱导的面波偏振异常进行了分析.结果表明,相对于长周期面波而言,短周期面波的振幅和偏振方向受横向非均匀性的影响更大,特别是偏振方向对地球结构的非均匀性更为敏感;切向分量存在横向非均性引起的Rayleigh与Love面波耦合现象;异常体边界处表现出强的面波波场响应. 相似文献
10.
11.
本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。 相似文献
12.
本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。 相似文献
13.
The application of Maslov asymptotic theory in a general 3-D mixed subspace of 6-D complex phase space is proposed to obtain the integral superpositions of Gaussian packets and beams. The ray method and the superposition of plane waves (Maslov method of Chapman and Drumond [7]) are special limiting cases of the above mentioned approach. The same high-frequency asymptotic expansion formulae for seismic body waves were derived previously in [8] using the Gaussian beam method. 相似文献
14.
A. B. Druzhinin 《Pure and Applied Geophysics》1996,148(3-4):637-683
A comprehensive approach, based on the general nonlinear ray perturbation theory (Druzhinin, 1991), is proposed for both a fast and accurate uniform asymptotic solution of forward and inverse kinematic problems in anisotropic media. It has been developed to modify the standard ray linearization procedures when they become inconsistent, by providing a predictable truncation error of ray perturbation series. The theoretical background consists in a set of recurrent expressions for the perturbations of all orders for calculating approximately the body wave phase and group velocities, polarization, travel times, ray trajectories, paraxial rays and also the slowness vectors or reflected/transmitted waves in terms of elastic tensor perturbations. We assume that any elastic medium can be used as an unperturbed medium. A total 2-D numerical testing of these expressions has been established within the transverse isotropy to verify the accuracy and convergence of perturbation series when the elastic constants are perturbed. Seismological applications to determine crack-induced anisotropy parameters on VSP travel times for the different wave types in homogeneous and horizontally layered, transversally isotropic and orthorhombic structures are also presented. A number of numerical tests shows that this method is in general stable with respect to the choice of the reference model and the errors in the input data. A proof of uniqueness is provided by an interactive analysis of the sensitivity functions, which are also used for choosing optimum source/receiver locations. Finally, software has been developed for a desktop computer and applied to interpreting specific real VSP observations as well as explaining the results of physical modelling for a 3-D crack model with the estimation of crack parameters. 相似文献
15.
本文提出了在计算机上实现分块三维地壳模型及利用加权最小二乘拟合生成平缓光滑的三维速度函数的方法,给出了适用于分块、块内速度连续变化的三维模型中Cauchy射线追踪的新算法,简介了基于上述方法反射线的基本理论所编制的合成三维理论地震图的程序包RSSGTD.给出的两个盆地状模型的算例表明,所使用的模型生成方法具有模拟复杂地壳结构的能力;与三维样条函数方法比较,最小二乘拟合方法能给出更加适合射线方法合成地震图计算的速度函数,并且内存小、计算速度快;所给出的Cauchy射线追踪算法能够适合块状模型中任何体波射线的追踪. 相似文献
16.
TTI介质是石油地震勘探领域最常用的各向异性介质,快速计算TTI介质射线路径和走时信息有重要的研究意义.TTI介质传统运动学射线追踪方法一般基于任意弹性介质射线方程,利用Bond变换或者四阶张量变换来处理复杂的21个弹性参数,因而非常耗时.实际野外对称轴统一的TTI介质模型,一般可以看成VTI介质模型旋转一定角度获得.为此,本文推导了三维VTI介质射线追踪方程,提出先在本构坐标系中进行VTI介质射线追踪,再通过坐标旋转将射线路径旋转至观测坐标系中,获得TTI介质射线路径.数值模型计算表明该方法高效和精确,较传统方法效率提高了近4倍.在强各向异性等特殊情况下,体波波前面都与理论群速度面一致. 相似文献
17.
M. CAMPILLO 《Geophysical Prospecting》1987,35(3):236-249
An approach that relies on a discrete representation of seismic wavefields allows the computation of synthetic SH-seismograms in a laterally varying medium with plane and curved interfaces in the two dimensional (2-D) case. The diffracting interface is represented by an array of body forces located along the interface at equal spacing. The numerical treatment is limited to the irregular boundary while the propagation in flat layered zones is obtained by the reflection-transmission matrix method. As an example we have studied the case of a dome in a stratified medium. The solutions obtained verify the reciprocity theorem with good accuracy. The computation of vertical profiles and of surface reflection profiles illustrates the effects of diffraction and the importance of lateral propagation in such a structure. 相似文献
18.
本文提出一个新算法,用来高精度计算三维不均匀地球模型中地震位错引起的地表以及空间固定点同震重力变化.具体地说,我们首先把实际三维不均匀地球分解成球对称地球模型和对应的横向不均匀增量,分别进行计算,二者对应的计算结果分别称为球对称解和三维响应.由于球对称解可直接利用球对称地球模型位错理论计算得到,本文的目标是计算三维响应即地球的横向不均匀结构对同震重力变化的影响.然后,我们把三维响应再分为震源的响应和地球横向不均匀结构的响应,它们可分别借助对震源函数的扰动以及对平衡方程式的变分求解.本文推导出六个特殊点源位错引起的地表以及空间固定点同震重力变化计算公式(一个垂直走滑位错,两个相互垂直的倾滑位错,三个开裂位错),对这些公式进行适当组合就可以计算任意位置任意类型位错产生的同震重力变化,对应的计算公式同步给出.接着,依据36阶P波速度模型,我们利用岩石试验经验关系式推导出三维S波速度模型,密度模型,位场模型以及重力模型.最后,利用上述三维模型,本文计算出三种典型类型的点源位错产生的同震重力变化,结果显示三维响应与位错类型,震源深度都有关系,其最大响应占球对称解的0.5%左右,且在所有影响因素中S波速度模型影响最大.数值结果同时表明,三维响应中震源的响应与地球横向不均匀构造的响应处于同一量级.本文给出的地表和空间固定点同震重力变化计算公式可分别高精度解析地表重力和卫星重力观测数据(GRACE、GOCE等),提高大地测量数据理论解析水平. 相似文献