共查询到20条相似文献，搜索用时 171 毫秒
1.
地壳和上地幔三维速度与界面结构层析成像——理论部分 总被引：2，自引：0，他引：2
提出了一种三维模型参数化方法,使得正问题和反问题可以用解析式表达．地壳和上地幔速度与界面结构分别由一组非常光滑的函数描述．在正问题和反问题中,使用了初值射线追踪方法,走时对模型网格节点参数的偏微商在射线追踪的同时进行了解析计算．由于速度和界面函数具有二阶连续偏微商,由射线路径的发散和聚焦在地表引起的影区可以避免．射线追踪之后,可以获得由矩阵和矢量组成的反问题．解决反问题使用了具有阻尼因子的奇异值分解法．一组由几条纵剖面组成的合成数据被用来进行数值模拟计算,结果表明方法是健全的．与二维方法比较,在地壳模型非常复杂的情况下,三维反演方法可以给出准确的速度结构和界面位置． 相似文献
2.
传统地震定位方法利用震源轨迹确定震源位置,但基于均匀或横向均匀介质模型必然导致定位误差。为此对传统方法进行改进,发展适用于三维复杂地壳速度模型的地震定位交切法。利用最小走时树射线追踪技术,以离散方式准确计算三维复杂地壳速度模型中的震源轨迹,将震源定位于震源轨迹交汇的密集点。将该方法应用于云南地区地震重定位,得到较高定位精度。 相似文献
3.
利用中国地震局地球物理勘探中心1998年在伽师地震区记录的天然地震资料，通过地震记录波形资料的衰减来反演伽师地震区地壳介质的三维Ｓ波Q值结构，并初步了解该地震区内地壳深部断层的大致形态．首先，对精选的450多次地震记录的Ｓ波求出其振幅谱，称为观测振幅谱；然后，假设震源符合布伦圆盘位错模式，在对仪器和场地效应校正之后，用非线性阻尼最小二乘法对理论振幅谱和观测振幅谱进行拟合，从而得到观测吸收特征时间；最后，利用三维射线追踪方法，对理论吸收特征时间与观测吸收特征时间进行非线性阻尼最小二乘拟合．在拟合过程中，速度模型保持不变．速度模型是通过三维走时层析成像得到的，只修改Q值模型，当拟合达到给定精度后，便得到最终的Q值模型．结果表明，在10～18ｋｍ 的深度范围内，存在一个北东向的低Q值条带；在12～18ｋｍ 的深度范围内，存在一个北西向的低Q值条带．这与他人所做的北东向、北西向的两个低速条带较接近．不同的是，低Q值条带深度范围比低速带的更大一些． 相似文献
4.
5.
利用地震波走时联合反演算法(改进型最短路径算法)进行三维弯曲地震射线追踪正演, 以及共轭梯度法求解带约束的阻尼最小二乘问题进行反演, 同时更新速度模型和地震震中位置, 结合地方震和区域地震走时资料得到了黄土高原(含汾渭断陷盆地)及邻区地壳三维P波速度结构. 其横向变化结果表明, 研究区地壳内的P波高速异常区与其内的地震活动构造带相一致, 地震多发生在P波高速异常区的边缘或高、 低速异常区的交汇处. 秦岭山区和鄂尔多斯块体东南区为P波低速异常区. 而垂向变化结果则表明研究区存在低速异常区. 相似文献
6.
《地球物理学进展》1990,(2)
多数井间地震走时层析成象法是将震源到接收器间的射线路径近似为直线,但这只有当研究区域速度变化较小时才有效.实际上许多地区都可观测到速度变化达到10—20%以上,引起显著的射线弯曲.虽然在研究井间成象问题时有人已考虑到这一非线性效应.但多数结果并不令人满意.本文提出一种二维射线追踪的迭代反演方法并成功地应用于野外实测数据. 本方法是先进行迭代射线追踪然后来修改速度模型.每次迭代都把在当前模型中经射线追踪计算得到的理论走时与观测走时之差通过线性方程组与未知的速度扰动联系起来,并通过阻尼最小二乘法求出用来修改模型的速度扰动.迭代一直进行到理论走时与观测资料拟合到给定误差,或走时拟合无法再改进时为止.理论模拟表明本方法经几次迭代后就可收敛到正确解. 利用这一方法处理了结晶岩地区的实测井间地震实验数据.地震波的频率范围是1—6.6 kHz,允许的分辨率为几米.得到的速度图象与其它地质和地球物理资料十分吻合.根据求解的速度模型计算出的马斯洛夫(Maslov)理论地震图与实际波形图相当一致,这也有力地验证了本反演方法的有效性. 相似文献
7.
8.
Pg回折波的上地壳三维层析成像 总被引：7，自引：1，他引：7
提出利用人工地震Ｐｇ回折波走时重建区域上地壳三维速度分布的方法。正问题利用经典最小二乘迭代求解垂向速度梯度模型的参数值及相应走时残差。反问题采用泛函空间的最小二乘准则来实现。对速度不需要预先进行参数化处理，可以计算反演区内任意点的速度值。 相似文献
9.
基于图形结构的三维射线追踪方法 总被引：39，自引：16，他引：23
在地震层析成像研究中,为了克服最小走时射线路径追踪方法存在的问题,对该方法计算过程中的关键步骤进行了改进．在节点走时的计算中引入Ｂｒｅｓｅｎｈａｍ画线算法;在最小走时节点查寻中,结合使用快速排序算法与插入排序算法,替代以往方法中多采用的堆排序算法;所采用的节点设置方式,可以引入速度界面,还可以实现反射波射线追踪．模型计算证明,改进的最小走时射线路径方法具有精度高,速度快的特点,所提出的三维空间反射波射线追踪算法简便易行。 相似文献
10.
为了尽快定位地震,在“着未着”定位算法基础上,尝试引入三维地壳模型,构建具有三维空间格点分布的走时表.根据已触发台站的到时和未触发台站的位置信息设计概率分布函数,通过八叉树搜索方法,快速给出震源在三维空间的可能位置.使用波前追踪算法,计算中国几个地区三维地壳模型的走时网格,利用中国地震台网资料,对区域内发生的地震进行定位分析.结果显示,在一定的台网密度条件下,三维实时定位方法能在震后数秒给出震源位置,可满足地震预警要求. 相似文献
11.
ＩｎｄｉｒｅｃｔａｐｐｒｏａｃｈｍｅｔｈｏｄｆｏｒｓｐｅｃｉｆｉｅｄｅｎｄｐｏｉｎｔｓｓｅｉｓｍｉｃｒａｙｔｒａｃｉｎｇｉｎｔｈｒｅｄｉｍｅｎｓｉｏｎａｌｉｎｈｏｍｏｇｅｎｅｏｕｓｍｅｄｉａＣＨＡＯＦＡＮＸＵ（徐朝繁）ＸＩＡＮＫＡＮＧＺＨ... 相似文献
12.
Ｉｎｔｒｏｄｕｃｔｉｏｎ３－Ｄｓｅｉｓｍｉｃｔｏｍｏｇｒａｐｈｙｈａｓｂｅｅｎａｐｐｌｉｅｄｔｏｖａｒｉｏｕｓｇｅｏｐｈｙｓｉｃａｌｐｒｏｂｌｅｍｓ．ＡｋｉａｎｄＬｅｅ（１９７６）ａｎｄＨａｗｌｅｙｅｔａｌ．（１９８１）ｉｎｖｅｒｔｅｄ３－Ｄｍｏｄｅ... 相似文献
13.
TTI介质是石油地震勘探领域最常用的各向异性介质,快速计算TTI介质射线路径和走时信息有重要的研究意义.TTI介质传统运动学射线追踪方法一般基于任意弹性介质射线方程,利用Bond变换或者四阶张量变换来处理复杂的21个弹性参数,因而非常耗时.实际野外对称轴统一的TTI介质模型,一般可以看成VTI介质模型旋转一定角度获得.为此,本文推导了三维VTI介质射线追踪方程,提出先在本构坐标系中进行VTI介质射线追踪,再通过坐标旋转将射线路径旋转至观测坐标系中,获得TTI介质射线路径.数值模型计算表明该方法高效和精确,较传统方法效率提高了近4倍.在强各向异性等特殊情况下,体波波前面都与理论群速度面一致. 相似文献
14.
对复杂山地介质的非均质性以及介质中地震波运动学特征进行深入研究,对于提高复杂山地区域地震勘探的效果有着重要的理论意义和实际价值.为了研究复杂山地非均质性和该介质中地震波的一些运动特性,提出了一种复杂山地随机介质的建模方法和一种新的射线追踪算法.与常规算法相比,复杂山地随机介质的生成方法采用更贴近实际介质特点的梯度介质作为背景介质,并在模型生成过程中加入地形修正步骤;新提出的GMM-ULTI射线追踪算法,充分融合群推进法、迎风思想、走时插值法的优势,采用先计算走时后追踪射线路径的两步策略完成射线追踪.算法分析与计算实例表明:复杂山地随机介质的生成方法能灵活、精细且更贴近实际地刻画复杂山地介质的非均质特点;新射线追踪算法兼顾精度和效率、能无条件稳定且灵活地适应复杂山地随机介质的特点;同时基于对几个模型试算结果的分析也得出了复杂山地随机介质中的地震波的一些传播规律. 相似文献
15.
16.
利用在青藏高原东部及其邻近地区记录到的１万余条近震到时资料,反演该地区的地壳上地幔三维速度结构。采用网格点模型描述三维速度结构,模型维数为２２２２６,网格点间距水平向为１００ｋｍ,垂直向为２０ｋｍ,网格点之间的速度值通过线性插值给出。采用改进了的快速三维射线追踪方法,确定三维非均匀介质中的地震射线路径和理论走时。反演结果显示,青藏高原南部的上地壳中（３０ｋｍ左右的深度）存在一低速区,这和面波反演的结果一致,羌塘块体下地壳有明显的低速异常带,青藏公路沿线的垂直速度剖面显示出岩石层受挤压增厚的构造特征。 相似文献
17.
Gaussian beam is an important complex geometrical optical technology for modeling seismic wave propagation and diffraction in the subsurface with complex geological structure. Current methods for Gaussian beam modeling rely on the dynamic ray tracing and the evanescent wave tracking. However, the dynamic ray tracing method is based on the paraxial ray approximation and the evanescent wave tracking method cannot describe strongly evanescent fields. This leads to inaccuracy of the computed wave fields in the region with a strong inhomogeneous medium. To address this problem, we compute Gaussian beam wave fields using the complex phase by directly solving the complex eikonal equation. In this method, the fast marching method, which is widely used for phase calculation, is combined with Gauss–Newton optimization algorithm to obtain the complex phase at the regular grid points. The main theoretical challenge in combination of this method with Gaussian beam modeling is to address the irregular boundary near the curved central ray. To cope with this challenge, we present the non-uniform finite difference operator and a modified fast marching method. The numerical results confirm the proposed approach. 相似文献
18.
19.
合成三维横向非均匀介质远震体波接收函数的Maslov方法 总被引：5，自引：3，他引：2
三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要．本文发展了基于Ｍａｓｌｏｖ理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪,并采用递归算法组织射线的追踪过程,从而在保证追踪精度的同时提高了射线追踪的效率,并可完全避免传统射线追踪过程中的射线编码问题．正演计算结果表明,此法是一种高效、实用的三维横向非均匀介质接收函数计算方法． 相似文献
20.