首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 874 毫秒
1.
赖晓玲  张先康 《地震学报》1997,19(5):506-516
研究利用反射波的走时反演三维弯曲界面和介质层速度的计算方法.各层界面利用分段非完全三次多项式描述.正问题采用一种快速的三维射线追踪方法.反演过程采用变阻尼最小二乘法.数值模拟结果表明,解很快收敛到真模型.处理了通过唐山震区的实测资料,重建了该震区莫霍界面三维构造形态,并揭示了该区域构造与地震活动的关系.   相似文献   

2.
Pg回折波的上地壳三维层析成象   总被引:1,自引:0,他引:1  
提出利用人工地震Pg回折波走时重建区域上地壳三维速度分布的方法。正问题利用经典最小二乘迭代求解垂向速度梯度模型的参数值及相应走时残差。反问题采用泛函空间的最小二乘准则来实现。对速度不需要预先进行参数化处理,可以计算反演区内任意点的速度值。数值计算结果表明了该方法的可靠性。处理了滦县地震区人工地震Pg走时资料,获得了该区上地壳顶部的三维速度分布。  相似文献   

3.
三维多层介质重力-地震同步联合反演   总被引:1,自引:0,他引:1       下载免费PDF全文
联合反演是地球物理勘探的重要解释手段,能够提高模型参数的反演精度.本文在归纳和分析重力与地震资料联合反演的研究和应用现状的基础上,利用三维多层介质模型的地震走时和重力正演公式,推导了地震走时和重力异常对界面深度的雅可比矩阵,实现了三维重力-地震同步联合反演界面成像.最后进行了数值理论模型模拟和实例计算,结果表明地震走时和重力同步联合反演很好的重建了三维多层介质界面.  相似文献   

4.
横向均匀介质可由若干水平分层足够好地逼近,本文引入零厚度层来描述成层介质,从而建立了回折波、反射波、首波和直达波统一的时距关系以及地震走时随模型参数变化的公式。本文详细讨论了水平分层介质里地震走时曲线的特征,合成地震走时的试验表明了由地震走时,用阻尼最小二乘法可以成功地反演速度模型参数,所反演的模型含有速度梯度层,也可以含有低速层。  相似文献   

5.
层状各向异性介质转换波克希霍夫叠前时间偏移   总被引:5,自引:2,他引:5       下载免费PDF全文
在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.  相似文献   

6.
模拟退火方法在三维速度模型地震波走时反演中的应用   总被引:5,自引:3,他引:2  
采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。  相似文献   

7.
三维复杂介质中转换波走时快速计算   总被引:22,自引:9,他引:13       下载免费PDF全文
复杂介质中转换波走时计算是多波勘探地震学中重要内容之一.本项研究利用惠更斯原理和费玛原理,获得了三维复杂介质中转换波快速计算的改进型最小走时树方法.其中,在保证精度的条件下,为了提高三维转换波走时计算效率,首先对初至波最小走时树基本算法进行了改进.本方法通过将转换波分为上、下行波分别进行射线追踪以实现三维转换波走时的快速计算.模型计算表明,方法的计算速度快,而且稳定性强,对多波地震勘探具有较大的应用价值.  相似文献   

8.
正交各向异性介质P波走时分析及Thomsen参数反演   总被引:10,自引:3,他引:10  
对于包含有垂向裂缝的横向各向同性地层或含有多组正交裂缝的各向同性地层,正交各向异性介质模型是最简单的与实际地层相符的方位各向异性模型.本文对单层水平反射界面正交各向异性模型采用射线追踪法计算了全方位角变化的P波走时,时距曲线表现出强方位各向异性.采用小生境遗传算法,对三条成一定角度的测线的走时信息进行速度和各向异性参数反演.模型算例表明,此方法可以得到高精度的裂缝方位角、P波垂直速度和较高精度的Thomsen各向异性参数.  相似文献   

9.
三维地质模型中地震波共轭梯度非线性走时反演   总被引:2,自引:0,他引:2  
地震体波走时层析成像是探测地球内部速度结构的重要方法之一。基于三维块状建模以及三角形拼接的界面描述方式,结合快速高效的逐段迭代射线追踪方法,获得三维复杂地质模型中的地震射线路径与走时信息,采用共轭梯度非线性反演算法,进行地震波走时反演。实验结果表明共轭梯度反演算法在三维层状模型中具有较高的有效性。  相似文献   

10.
地震走时数据可以反演近地表的速度,但不能反演出隐蔽层和低速层。航空电磁数据可以反演近地表的高电阻率和低电阻率,但是对垂直方向的分辨率低。联合三维地震走时数据和航空电磁数据反演了近地表的速度和电阻率结构,并采用棋盘模型测试和实际数据应用展示了联合反演出的速度模型优于单独的走时反演出的速度模型。结果表明:该算法可以应用到处理数据量大的勘探地震成像中,能提供优化的速度结构。  相似文献   

11.
Pg回折波的上地壳三维层析成像   总被引:7,自引:1,他引:7  
提出利用人工地震Pg回折波走时重建区域上地壳三维速度分布的方法。正问题利用经典最小二乘迭代求解垂向速度梯度模型的参数值及相应走时残差。反问题采用泛函空间的最小二乘准则来实现。对速度不需要预先进行参数化处理,可以计算反演区内任意点的速度值。  相似文献   

12.
计算最小走时和射线路径的界面网全局方法   总被引:27,自引:16,他引:27       下载免费PDF全文
用慢度分块均匀正方形模型将介质参数化,仅在正方形单元的边界上设置计算结点,这些结点构成界面网.根据Huvsens和Fermat原理,由不断扩张、收缩的波前点扫描代替波前面搜索,在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似,通过比较确定最小走时和相应的次级源位置,记录在以界面网点位置为指针的3个一维数组中.借助这些数组通过向源搜索可计算任意点(包括界面网以外的点)上的全局最小走时和射线路径.这一方法不受介质慢度差异大小限制,占内存少,计算速度较快,适于走时反演和以Maslov射线理论为基础的波场计算.  相似文献   

13.
我们研究了由近震得到的地震走时资料估计地震波速度三维结构的反演方法。这种方法由同时确定速度界面的二维深度分布、三维速度分布、台站校正和震源来表征。速度界面的深度分布和慢度扰动的分布以纬度、经度和深度的幂级数来拟合。这种三维地震结构的表示不仅有助于减少未知参数的个数,也有助于解析地求取界面深度、速度和走时关于未知参数偏导数的值。连接震源和台站之间的射线轨迹,是在具有确切界面深度分布的简化的速度分布条件下由称为试射法的射线追踪方法确定的。数值试验证明,我们的反演方法可成功地估计三维速度结构。  相似文献   

14.
华北克拉通中东部岩石圈减薄研究已经成为近年固体地球科学研究的热点,同时该区域也是中国强震的多发区.在华北及邻近地区已完成了二维人工地震宽角反射/折射探测(简称DSS剖面探测)剖面42条,文章收集了这些二维探测剖面的研究成果,在对二维剖面速度结构和界面结构进行网格化处理的基础上,采用克里金(Kriging)插值方法构建了华北克拉通中东部地壳三维速度结构模型HBCrust1.0.采用走时拟合的办法对模型进行了检验,结果表明,用HBCrust1.0地壳模型计算得到的初至波走时与实际台站观测到的初至波走时之间拟合良好.研究结果表明,上地壳是主要的孕震层,C界面(上、下地壳分界面)附近区域是地壳内部的脆韧转换带.唐山地震震源区下方Moho面深度有明显的变化,其低速结构一直延续到下地壳,推断唐山地震的应力积累与地幔物质运移和变形密切相关.克拉通东部华北盆地中心区域的下地壳平均速度明显要高于周围地区,可能是地幔岩浆长期底侵作用的结果.  相似文献   

15.
快速行进法(FMM)是一种求解程函方程数值解计算网格节点走时,然后向后处理进行射线追踪的方法.为了求取任意起伏界面下高精度多震相的地震走时与相应的射线路径,本文采用任意起伏地表条件下的的三维不等距上行差分公式结合分区多步计算技术实现了三维复杂层状起伏介质中多震相(透射、反射、转换波)地震走时的计算,利用上行有限差分公式逐次进行射线路径的追踪,并且通过与较为成熟的不规则最短路径法(ISPM)对比,验证了本算法的计算精度和有效性.数值模拟实例和对比结果表明该算法具有较高的计算精度,数值计算稳健,能灵活处理含任意三维起伏界面模型中多震相地震走时及相应射线路径的追踪问题.  相似文献   

16.
本文讨论了射线和射束理论地震图对速度模型和弯曲界面变化的敏感性.Nowack和Lutter(1989)以前的成果给出走时和射线振幅相对于光滑变化的速度变化导数,在实施构造的线性化最大似然反演时这些导数是需要的.在射线近似中,局部应用Snell定律、校正波前曲率以及使用局部平面波反射/透射系数来体现光滑界面,利用费马原理可沿着原始射线轨迹直接计算走时的偏导数.对于反射/透射射线的振幅扰动,必须考虑扰动后的两点射线轨迹的射线偏移.本文遵循的方法是利用没有额外射线追踪的扰动理论近似计算扰动后的两点射线,然后沿着这条近似的两点射线直接计算扰动的射线振幅,包括修正过的反射/透射系数和几何扩展.进行了若干次数值实验,利用走时和振幅反演速度和界面形状来检验所得出的偏导数算子、走时和振幅反演结果,还与对大尺度性质不太敏感而对不均匀曲率较敏感的振幅作了比较.  相似文献   

17.
中国地震台网初至P波区域三维走时表的建立   总被引:1,自引:0,他引:1       下载免费PDF全文
基于LLNL-G3Dv3全球P波三维速度模型,应用FMM软件包计算并建立了中国地震台网990个台站的初至P波区域三维走时表.该走时表覆盖了以台站为中心的水平向20°×20°、 垂直向-5.1—80 km (向下为正)的三维空间. 其水平向间隔为0.2°,垂直向间隔为5 km.这样对于任一深度小于80 km的震源,均可以应用此三维走时表计算其到周围10°范围内台站的走时.中国地震台网初至P波区域三维走时表的建立,对于改善区域初至P波走时预测,提高地震定位精度有一定现实意义.   相似文献   

18.
高斯波包反射走时速度反演方法   总被引:1,自引:1,他引:0       下载免费PDF全文
李辉  殷俊锋  王华忠 《地球物理学报》2017,60(10):3916-3933
扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.  相似文献   

19.
天山-准噶尔地区地震层析成像与壳幔结构   总被引:2,自引:0,他引:2       下载免费PDF全文
通过对东西天山两条天然地震剖面的远震P波与近震P波走时资料以及2001-2007年新疆地方台网的近震走时资料分析,获得了天山-准噶尔地区的三维地震层析图像.层析反演使用的台站数140个,地震事件1132个,P波走时数24904条.检测板试验结果指出,东西天山剖面、天山轴部和东准噶尔地区P波速度扰动恢复能力相对高.纵波波...  相似文献   

20.
利用面波频散数据进行各向同性和各向异性横波速度结构成像是研究全球和区域构造变形、获取地壳上地幔和近地表精细结构的一种非常有效的方法.传统的基于频散数据的面波成像通常采用两步法,即首先反演二维相速度或群速度分布图,然后再基于格点的纯路径频散反演格点下方的一维横波速度模型,之后再拼合成三维速度模型.在本文中我们较为详细地综述新发展的基于面波频散走时的三维面波直接成像方法体系,即面波一步法成像.该方法体系包括通过所有路径、不同频率的面波频散走时直接反演三维各向同性横波速度结构的方法(DSurfTomo)、直接反演三维方位各向异性横波速度结构的方法(DAzimSurfTomo),以及直接反演三维径向各向异性横波速度结构的方法(DRadiSurfTomo).新的面波直接成像方法计算不同周期的面波射线路径,从而更好地考虑复杂结构下面波射线路径弯曲对成像精度的影响.随后我们介绍面波直接成像方法的应用实例,包括地壳上地幔及地壳浅部的多尺度各向同性和各向异性横波速度结构的成像,这些成像研究为认识区域构造演化、孕震构造、断裂带浅部结构、成矿构造、城市地下结构等提供了重要的约束.最后我们讨论面波直接成像方法...  相似文献   

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

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