共查询到16条相似文献，搜索用时 187 毫秒
1.
本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。 相似文献
2.
本文研究了横向非均匀吸收介质中的Maslov渐近理论地震图方法.对于如何提高计算精度及计算速度作了比较深入的探讨.本方法先计算远场格林函数,由褶积得到一般震源函数的弱吸收理论地震图.对于高Q引起的数值不稳定性,通过引入平滑函数得以解决.由此产生的新的被积函数具有非常好的收敛性,从而提高了计算精度及计算速度.以此理论编写的程序在模拟完全弹性介质以及一维吸收介质时与反射率法作了比较,并取得了令人满意的效果.对于二维横向非均匀吸收介质,计算方法与高斯光束方法作了对比,二者在主要特征上基本一致.数值对比中出现的一些差异也作了必要的讨论及说明. 相似文献
3.
研究了横向非均匀介质中的Ｍａｓｌｏｖ面波渐近理论，在横向弱非均匀介质的假设下，介质的纵向非均匀性反应在局部本征函数中，以局部本征函数近似真本征函数是射线理论的直接推论．由此，三维结构下的面波计算退化为准二维问题．由于本文方法属于慢度法，面波的频散使得在一般情况下得不到与体波ＷＫＢＪ方法相似的褶积结果；在震源函数为高斯波包的假设下，得出了与二维体波Ｍａｓｌｏｖ理论图形式上完全相同的褶积结果．还讨论了吸收介质中的面波波包理论图计算，最后结果与二维体波吸收介质中的结果相似． 相似文献
4.
5.
6.
合成三维横向非均匀介质远震体波接收函数的Maslov方法 总被引：5，自引：3，他引：2
三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要．本文发展了基于Ｍａｓｌｏｖ理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪,并采用递归算法组织射线的追踪过程,从而在保证追踪精度的同时提高了射线追踪的效率,并可完全避免传统射线追踪过程中的射线编码问题．正演计算结果表明,此法是一种高效、实用的三维横向非均匀介质接收函数计算方法． 相似文献
7.
Maslov渐近理论与辛几何算法 总被引：26，自引：9，他引：17
为了克服地震层析成像中的焦散问题,本文系统地研究了Ｍａｓｌｏｖ渐近理论与辛几何算法,同时提出了一种基于辛几何算法计算Ｍａｓｌｏｖ波场的数值计算方法,并就其中射线追踪这一重要环节,利用辛算法和非辛Ｒｕｎｇｅ－Ｋｕｔｔａ方法进行了数值模拟．计算结果表明,这两类算法在精度上相差无几,但辛算法的速度要快;在Ｈａｍｉｌｔｏｎ量保持方面,辛算法具有非辛Ｒｕｎｇｅ－Ｋｕｔｔａ方法无可比拟的优越性． 相似文献
8.
作为近地表横波速度结构成像的主要手段之一,面波多道分析法的正问题研究对现场观测系统设计及后续反演计算具有重要意义.目前面波频散曲线的正演主要分为两类:一是对水平层状介质中面波的本征值问题进行求解,该类方法计算效率高但较难考虑地下介质在横向上的不均匀性;二是基于波动方程的全波场模拟,该类方法在理论上可考虑任意复杂的地质模型但计算成本相对较高.本文基于振幅归一化加权的聚束分析,提出了一种适用于横向非均匀介质模型的多道瑞雷波频散曲线正演方法.首先,基于聚束分析的计算公式推导得到了经振幅归一化加权后输出功率谱中相速度与局部相速度之间的关系,然后通过黄金分割极值搜索算法计算得到了多道瑞雷波数据的理论频散曲线.数值分析结果表明,该算法能够快速地实现横向非均匀介质中多道瑞雷波频散曲线的正演计算,所求取的频散曲线与采用二维弹性波时间域有限差分模拟分析得到的结果误差较小,这在一定程度上说明了该计算方法的可靠性,从而可为面波多道分析法中的观测系统快速优化设计以及横向非均匀介质中频散曲线的反演解释提供理论支撑. 相似文献
9.
采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。 相似文献
10.
除介质各向异性之外,地球内部介质的横向非均匀性也是控制面波速度变化的重要因素.本文基于振型耦合和多重散射的地震波传播理论,数值模拟并分析了在具有均匀介质背景的三维异常体——上地幔横向非均匀介质中传播时,地震面波的振幅与偏振等动力学响应参数;其中分别模拟了不同周期入射、不同角度入射和不同尺度非均匀介质模型等多种情形下面波波场,并对横向非均匀性诱导的面波偏振异常进行了分析.结果表明,相对于长周期面波而言,短周期面波的振幅和偏振方向受横向非均匀性的影响更大,特别是偏振方向对地球结构的非均匀性更为敏感;切向分量存在横向非均性引起的Rayleigh与Love面波耦合现象;异常体边界处表现出强的面波波场响应. 相似文献
11.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7％左右. 相似文献
12.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7％左右. 相似文献
13.
本文提出一种结合广义射线和有限差分理论而形成的混合方法（以下简称混合方法）．在介质均匀区（包括平行介质分层区）应用广义射线理论，在介质不均匀区应用有限差分理论．在这种结合过程中，把均匀区广义射线的结果作为不均匀区有限差分区波场连续的输入．对于均匀介质模型，混合方法理论地震图和广义射线理论地震图、广义反射透射系数矩阵及离散波数法地震图的对比，取得了满意的结果，在数值上证明了该方法的正确性．计算了１９９０年８月５日Ｈａｋｏｎｅ地震前震在日本Ａｓｈｉｇａｒａ山谷盲测点的响应，理论地震图和实际资料的比较同样取得好的结果．这种方法可用于地震勘探、工程地震、强地面运动的研究． 相似文献
14.
The imaginary wavefront of a common midpoint gather of seismic data is analysed in terms of elementary wavefronts. The wavefronts emitted by the sources are split into a sequence of elementary circular wavefronts. Each such wavefront is independently transmitted through the subsurface and received by a group of receivers associated with the sources. The theory is first demonstrated in the case of a one-layer medium and is then extended to the case of a multilayer medium. We show that pairs of imaginary stationary points F and F' occur in the first medium through which all the rays, corresponding to the group of sources and receivers associated with the elementary wavefront pass. The location of the points is determined by 2n equations with 2n unknowns; (2n - 1) of these equations correspond to n layers of the medium, and an additional equation is given by the geometry of the CMP gather. Applying Fermat's principle between the points F and F', we obtain an algorithm for CMP ray tracing. The wavefront sampling problem is considered but not theoretically solved. Results obtained in the general case (a multilayer medium with lateral velocity variation) are presented and compared with those obtained by other modelling techniques. 相似文献
15.
长白山天池火山区上地壳Q值结构 总被引：4，自引：0，他引：4
同天然地震资料相比，人工地震剖面观测资料具有震源球对称且位置已知，观测点位密集、震源和观测点呈直线排列、速度结构已知等突出优点。本文提出了一种利用人工地震振幅资料反演分层介质中Q值的新方法。在利用该方法得到Q值深度函数的基础上，综合各炮点的结果，可研究Q值沿剖面的横向变化。利用上述方法计算了长白山L1剖面上地壳的Q值值结构。发现长白山天池火山区下方Q值明显低于邻近地区，该区域位于天池地区壳内低速异常体的上部，对应于已往的火山喷发通道。 相似文献