首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 440 毫秒
1.
介绍了在倾斜均匀各向同性层状介质情况下,根据地震反射走时反演二维地层结构和地层速度的地震射线层析成像方法——走时反演。正演模型射线追踪是根据Fermat原理实现的,即求解满足该原理的非线性方程组得到射线与界面的交点,进而计算相应的走时。反演是先假设一初始模型,用最优化方法使射线追踪走时与观察走时的残差极小。最后计算了分辨矩阵和信息密度矩阵,以评价反演结果。对有噪情形也进行了反演。  相似文献   

2.
地震走时层析成像法是一种利用不同地震台站连续接收到的地震事件,利用地震走时断层扫描反演地球内部速度分布的方法,主要分为正演(射线追踪)和反演两个环节。首先正演利用基于程函方程的快速行进法,实现了对入射波走时场和反射波走时场的计算,接着进行了射线追踪与走时灵敏度矩阵的计算,最后反演采用了子空间算法实现了反射波走时反演,反演成像结果可以分辨出异常体的形状和位置。  相似文献   

3.
为解决深海水体速度建模问题,建立了两个利用Munk公式和海底反射走时反演深海水体速度的反演理论,一个是射线走时反演理论,另一个是波动方程走时反演理论。具体地讲,是建立了地震波走时关于Munk参数的Fréchet导数和海森矩阵公式。为得到射线走时的Fréchet导数和海森矩阵,利用了3种不同的方法,即变分法、微扰法和Taylor级数展开法。这3种方法在理论上完全等价,其最终结果也完全相同。为得到波动方程走时反演理论中的Fréchet导数和海森矩阵,利用复合函数求导的链式法则对文献中的有关结果进行变换,以将其从面向速度本身的Fréchet导数和海森矩阵变换为面向Munk参数的Fréchet导数和海森矩阵。此外,还提出了实现利用Munk公式和海底反射走时构建深海水体声速剖面的5个策略,包括利用解析公式计算海底反射走时的由下至上的射线追踪策略,选用最简单的海底边界条件把数值模拟区域限定为水体的策略,用已知的深海声道轴深度和交替反演策略将四参数反演问题转简化为单参数反演问题的策略,由点到面,利用一维声速剖面公式实现深海水体的维速度建模的策略,以及用最简单的算法实现走时反演的最优化计算策略等。与常规走时反演相比,利用Munk公式进行的海底反射走时反演最多只需要反演4个参数,并可利用交替反演将四参数反演化为单参数反演,大大地减少了走时反演的计算量。  相似文献   

4.
作为一种模拟地震波走时的传统途径,射线追踪不仅对地震波理论研究具有重要意义。而且也可直接应用于地震波反演及偏移成像等过程。在一般条件下,射线追踪只能用数值方法来实现,但当介质上具有垂向不均匀性时,射线路径及其走时可以由解析式给出。根据适当的解析公式,可建立反射点和射线路径以及走时之间的关系,并用VC^ 6.0编制了基于图形界面的计算程序。这些工作对于逐步开发地震波的偏移和反偏移程序具有基础性的意义。  相似文献   

5.
井间地震层要成象在环境工程勘探中的应用前景广阔,走时层析成象的射线追踪方法和在型稀疏矩阵的方法的方法是其中的关键。文中讨论了LSQR算法系列,数值模拟证明了它的优越性,最后应用于鹏飞大厦场地,取得了良好的效果。  相似文献   

6.
张建南 《物探与化探》2006,30(4):319-321,326
论述了用数字方法设计出地层弹性波速度模型的原理和方法,阐述了地震射线追踪法的算法并用其对数字模型的反射、折射波时距曲线做了计算,通过实例分析了复杂速度结构与地震剖面的对应关系,可指导解释人员根据实际地质情况任意设计地层模型进行计算,很快得出结果。  相似文献   

7.
井间地震层析成像技术能够提供钻孔之间地质体的高分辨图像.LTI法与SIRT法相结合是目前使用较广的一种层析成像算法,在此基础上考虑到一个较好的初始模型能提高算法的准确性及效率性,以及在SIRT算法当中如何分配残差能使反演的结果更好,提出了改进模型的层析成像算法.此方法基于LTI法进行射线追踪得到每个接收点在初始模型上的理论走时后,就利用这个理论走时与实际拾取走时的差值以及射线追踪所得到的射线路径,对初始模型进行修正.通过数值模拟实验发现,该算法的效果较原始算法有一定的提高.  相似文献   

8.
跨孔雷达走时层析成像主要利用雷达波的走时进行反演,走时提取的正确与否将直接影响到层析成像的效果。数字图像分割法基于凸集投影(POCS)方法,使用能量比彩色图像分割技术准确提取走时。数字图像分割法提取走时首先应用在折射地震波的数据处理中。笔者首次将数字图像分割法提取走时的方法应用到跨孔雷达走时层析成像中,使用迭代线性反演算法重建了雷达波速度场。反演过程中,使用最小二乘QR分解法(LSQR)求解线性方程组,利用弯曲射线追踪技术构建雅可比矩阵,走时的计算值则由多模板快速推进算法(MSFM)得到。为了验证数字图像分割法在走时层析成像中的效果,使用一组合成数据及一组实测数据分别对比了基于数字图像分割法和基于传统能量比法获得的层析成像反演结果。对比结果表明,使用数字图像分割法得到的层析成像结果更为精确,误差更小,能为判断地下雷达波速度场提供更为有力的帮助。  相似文献   

9.
位于土耳其西部伊兹密尔市hltlndag地区的滑坡是其他自然灾害的诱发因素之一,利用电阻率成像法(ERT)和地震折射层析法(SRT)对该滑坡区开展地球物理勘查。在滑坡体上沿南-北和东-西剖面采用了温纳-施兰贝尔电极排列方式做了4个剖面的电阻率层析成像勘查和采用纵波地震检波器沿着与南-北向电阻率剖面一致的一条剖面开展地震折射层析勘查。使用最小二乘法反演技术处理电阻率和地震数据。利用一种不基于射线追踪的方法对地表折射数据进行初至波走时反演。这种方法通过对波走时的函数描述导出雅可比矩阵,而该矩阵由基于单元慢度扰动的有限差分近似法计算得出。利用程函方程求解走时。通过利用这些方法均获得了有关滑坡体内部结构、物理特征和滑动面几何形状的有价值的结论。滑坡体物理的特点是电阻率和地震波速度均较低。电阻率成像结果同时表明,这些低电阻率区与滑坡体中水、粘土含量较高有关。把固结的碎屑岩层假定为滑坡区基岩,该碎屑岩层的物理特点是电阻率相对较高(中等)和地震波速度很高。通过南一北向剖面的电阻率和地震折射数据综合解译,有助于我们确定滑动面的几何形状和滑坡体厚度的变化。沿该剖面存在一个不断变化的滑动面,且在剖面中部和北部(坡脚区)滑坡体的厚度较大。此外,我们认为滑坡体的水含量在滑坡体运动中起着关键作用。  相似文献   

10.
井间地震速度反演的遗传算法   总被引:1,自引:0,他引:1  
井间地震速度反演具有重要的意义,不仅对于准确圈定地下构造的传统地震勘探是重要的,而且速度的变化可以指示非背斜型的地层岩性圈闭.传统的射线追踪方法一般是先给出入射角,然后根据折射定律修改入射角以达到反演的目的.该方法对每条射线均需通过扫描确定入射角,计算量较大,对初始地质模型的要求较高.遗传算法具有全局收敛性,是一种不用梯度信息的优化方法,特别适用于大型的组合优化问题.采用多项式展开表示界面深度和速度,通过弯曲射线追踪算法来计算射线的旅行时间,利用遗传算法进行迭代优化,可以同时反演出地下复杂构造的界面形态和速度变化.该方法与网格化速度反演相比,减少了未知参数,较好地克服了多解性问题,试验结果表明,该方法可以达到较高的反演精度.  相似文献   

11.
利用线性旅行时插值射线追踪对近地表模型进行正演计算,可以快捷、准确地获得初至波走时和射线路径。由于该算法计算的初至波不局限于折射波,因此很好地解决了浅层折射勘探中的低速"隐蔽层"问题;而且,由于该算法是基于网格划分和线性插值,因此它不仅可以追踪任意复杂介质的初至波,而且可以使得追踪的初至波射线路径逼于真实,避免了同类算法直接连接网络节点形成射线路径的缺陷(路径过于弯折,计算走时偏大)。将LTI算法同其他几种算法的追踪结果进行的对比和分析表明,LTI算法在计算初至波走时和射线路径方面较其他算法更为精准、稳定,是一种有效的射线追踪方法。  相似文献   

12.
为了实现适应崎岖海底、大陡坡、海底火山等复杂海底地质条件且灵活、稳定、精度高的射线追踪方法,并基于该方法详细分析复杂海底条件下各种地震波型的运动学特征,综合多种算法的优势,实现了一种快速推进迎风插值射线追踪方法。首先,采用混合网格法剖分复杂海底地质模型;其次,通过融入迎风差分思想的线性插值策略来构建精度高、无条件稳定且灵活的局部走时和射线路径计算公式;然后,综合应用这些公式和多级次快速推进法,灵活计算整个模型各种类型地震波的走时,并基于逆向追踪方法计算射线路径;最后,对该方法进行了精度分析,发现其能够获得相对高的走时和射线路径计算精度,且反射波的计算精度远高于入射波。此外,计算实例分析显示,初至波中富含折射波和陡倾构造的反射波在很大偏移距还能被接收,崎岖海底各种波型传播路径复杂;基于此提出加大采集排列长度和采用直达波走时可提高复杂构造成像质量等复杂海底地震数据采集与处理方面的思考与建议。  相似文献   

13.
对于山地地震勘探而言,地形起伏,不能直接使用Herglotz-Wiechert公式反演近地表速度模型,这里采用不同于前人在炮域处理的思路,而是在CMP域建立了一组基于CMP浮动基准面对起伏地表上的炮检对走时和偏移距进行校正的公式,首先用本文提出的公式去除起伏地表对初至波的影响,然后应用Herglotz-Wiechert公式快速反演地形起伏地区近地表速度模型。该方法采用地震反射波的CMP浮动基准面校正思路,包含三个关键步骤:①产生CMP浮动基准面;②实测时距曲线的校正;③计算射线参数和速度-深度函数的计算。此方法无需射线追踪,无需初始模型,计算速度快,反演得到的模型具有和地震数据相同的水平分辨率,该模型可用于计算基准面静校正量,或用于叠前深度偏移速度场的表层建模,也可以作为回转初至波层析反演或全波形反演的初始模型。  相似文献   

14.
本文提出的地震反射层析成像方法,基于射线原理,应用叠前反射波的走时资料,逐层反演地下介质结构。它是在局部范围内,将地层各界面视为平面,求取界面、速度参数,全剖面的成像只需反复进行局部反演。用该方法处理了三板溪水电站工区的两个地震剖面,处理结果与钻孔资料对比,一致性较好。  相似文献   

15.
等时叠加波动方程叠前正演   总被引:6,自引:3,他引:6  
非零炮检距地震正演模拟在复杂区叠前地震资料处理和解释中有着重要意义。通常使用射线追踪法和全程波动方程法来实现这一模拟。通过提出等时叠加原理和数学检波器表达式,形成了非零炮检距地震正演的新一类方法——单程波动方程法。该方法与射线追踪法相比能保持波的动力学特征,正演结果与实际地震记录相似性好。并且在计算过程中,人工干预的程度低,操作简便。与全程波动方程法相比,该方法计算效率高,且不产生直达波、多次波等干扰,信噪比高,与目的层相对应的一次波突出,波型识别难度不大。因此,单程波动方程法特别适合于复杂油气藏反射地震资料的正、反演研究。  相似文献   

16.
从射线在整条路径上满足Fermat原理出发,采用一阶Taylor展开等办法,推出一个适用于任意层状地质结构下近似求解射线路径的对称正定三对角矩阵方程,通过求解该方程并结合迭代技术,从而实现对整条射线路径的全路径迭代计算。模型计算表明,本方法对于任意复杂结构和随机速度分布情况,其计算速度比传统的打靶法快,精度可根据需要预先确定,为复杂结构的地震偏移和层析成像提供了一种有效的射线追踪方法。   相似文献   

17.
快速精确计算走时在层析成像反演中起重要作用,本文采用有限差分法来计算。首先通过变换,将程函方程化成一守恒形式;然后再差分离散,得到一阶和二阶精度的差分格式,并证明了差分格式条件稳定。模型试算表明:本文算法简明稳定,能快速高精度地计算地震走时,是一种值得注意和采用的走时计算方法。  相似文献   

18.
作者在文中对三维任意界面情况下的两点间射线追踪问题进行了研究,从Fermat原理出发,对射线路径上的旅行时,采用一阶Taylor不完全展开等办法,推出了三维情况下计算射线路径的整体迭代算法。并通过模型计算证明了本方法的稳健性以及计算的高精度和高效率,为三维结构下地震偏移和层析成像提供了一种可行的射线追踪方法。  相似文献   

19.
三维层状地质结构下整体迭代射线追踪法   总被引:2,自引:1,他引:2  
作者在文中对三维任意界面情况下的两点间射线追踪问题进行了研究,从Fermat原理出发,对射线路径上的旅行时,采用一阶Taylor不完全展开等办法,推出了三维情况下计算射线路径的整体迭代算法。并通过模型计算证明了本方法的稳健性以及计算的高精确度和高效率。为三维结构下地震偏移和层析成像提供了一种可行的射线追踪方法。  相似文献   

20.
层析反演静校正   总被引:8,自引:0,他引:8  
表层模型层析反演是一种非线性模型反演技术,它利用地震初至波射线的走时和路径反演介质速度结构,不受地表及近地表结构纵横向变化的约束。根据正演初至时间与实际初至时间的误差,修正速度模型,经反复迭代,最终达到要求的误差精度。求取静校正时采用射线法计算炮点到检波点的旅行时,从而得到基准面校正量。  相似文献   

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

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