共查询到20条相似文献，搜索用时 15 毫秒
1.
为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低. 相似文献
2.
提出一种对旅行时进行抛物线插值的地震射线追踪方法（简称PTI方法）,它比基于旅行时线性插值方法（简称LTI方法）计算结果更准确.PTI和LTI方法都是基于2D网格单元模型,用于计算地震波的旅行时和射线路径.首先介绍了相关方法的一些基本概念.旅行时和射线路径都是在网格边界上进行计算的,因此,射线路径在同一恒速网格内是直线.其计算过程有两步.第一步,计算旅行时,第二步追踪射线路径.然后给出了LTI算法的基本公式.因为在炮点网格内可能存在折射波,文章也相应导出了其公式.最后详细推导了PTI 方法的公式.通过模型试算对比说明,用PTI方法较LTI算法更精确、更有效,PTI方法是一种很有发展前途的地震射线追踪算法. 相似文献
3.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
4.
层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行走时层析速度反演是当前比较成熟有效且广泛应用的技术.本文从高斯束偏移成像条件出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域走时扰动与速度扰动的线性关系,建立了成像域走时层析方程及其显式表达的层析核函数.该核函数的本质是有限频层析核函数,利用该核函数替换常规射线层析核函数可以明显提高层析反演精度.该核函数的计算关键是背景波场格林函数的计算,本文利用高斯束传播算子计算格林函数进而得到走时层析核函数,实现方式灵活高效且计算精度较高.基于高斯束传播算子的偏移成像与层析成像相结合进行深度域建模迭代,体现了速度建模与偏移成像一体化的思想.数值计算及实际数据应用证明了基于高斯束传播算子的成像域走时层析方法的有效性. 相似文献
5.
最短路径射线追踪方法及其改进 总被引：35，自引：9，他引：35
综述了用网络最短路径算法求解地震射线追踪问题的原理、方法技术以及存在问题和改进措施。特别介绍了作者在最短路径算法基础上，提出的动态网络最短路径地震射线追踪方法。该方法先采集从炮点到整个模型所有节点上的初至旅行时，其中，在一个单元内，对相邻每对已计算出最小旅行时的节点进行线性插值，并利用Fermat原理计算未知节点的最小旅行时；然后，利用同样的方法，从接收点开始，反向追踪炮点到接收点的射线路径、该方法能适于各种复杂的非均匀介质，极大地提高了射线追踪的精度。 相似文献
6.
近地表速度结构通常是利用射线走时层析或菲涅尔体走时层析等反演方法得到的,但它们的目标函数仍利用射线走时残差构建,导致反演精度不高.为此,本文提出了基于散射积分算法的初至波相位走时层析成像方法.该方法的核心是:(1)提出了依赖于频率的相位走时概念;(2)利用依赖于频率的相位走时信息,而非单一的无限频率射线走时;(3)发展了一种改进的相位展开方法,即通过监测相位不连续性和2π周期判定来消除相位折叠现象;(4)考虑了地震波传播的有限频特征,即基于波动理论而非传统的射线路径或有限空间的菲涅尔体构建核函数.通过利用Overthrust模型的数值实验及与传统射线走时层析和菲涅尔体走时层析的对比表明:本文提出的方法是一种有效的初至波走时反演方法.同时,基于Overthrust模型的数值试验还证明了下列结论,即通过挖掘更多的走时信息的确可以获得更高的反演精度和分辨率. 相似文献
7.
针对传统射线层析存在的种种局限性,菲涅尔带走时层析成像摒弃了传统的数学射线,考虑到地震信号具有一定的频带宽度,中央射线附近的介质对地震波的传播产生不同程度的影响。本文提出了多频段组合菲涅尔带走时层析成像方法。该方法以频率域波动方程Born和Rytov近似为基础,推导出建立在带限地震波理论基础上的波动方程 Rytov 近似走时敏感核函数,实现第一菲涅尔带约束下的波动方程走时层析反演方法。同时由于多个频段的引入,充分利用低频段和高频段的特有优势,从而兼顾菲涅尔带层析的计算效率与分辨率。模型试算结果证明了本方法的有效性和稳定性。 相似文献
8.
本文基于VTI介质拟声波方程,利用散射积分原理,在Born近似下导出了速度与各向异性参数的敏感核函数,同时结合作者前期研究提出的矩阵分解算法实现了一种新的VTI介质多参数全波形反演方法.矩阵分解算法通过对核函数-向量乘进行具有明确物理含义的向量-标量乘分解累加运算实现目标函数一阶方向或二阶方向的直接求取,从而避免了庞大核函数矩阵与Hessian矩阵的存储,该方法同时可以大大降低常规全波形反演在计算二阶方向时的庞大计算量.为了克服不同参数对波场影响程度的不同,本文利用作者前期在VTI介质射线走时层析成像研究中提出的分步反演策略实现了多参数联合全波形反演.理论模型实验表明,本文提出的基于Born敏感核函数的各向异性矩阵分解全波形反演方法可以获得较好的多参数反演结果. 相似文献
9.
《应用地球物理》2019,(3)
变偏移距VSP反射波射线路径是不对称、覆盖次数也不均匀,因此通常使用偏移速度分析、层析反演建立准确的2D、3D速度场,尤其角度域共成像道集最适合偏移速度分析,通常拾取其剩余曲率更新速度,建立准确的速度场。为了建立高精度速度模型,本文推导了VSP角度域共成像道集剩余深度与旅行时残差的关系,用最小二乘求解更新速度,建立了一个四参数层析反演方程,包括地层倾角、入射角、剩余深度、灵敏度矩阵等四个参数。在角度域中将反射波分解为上行、下行两个透射波计算其旅行时。采用"逐步-累加"方式进行复杂地表叠前深度偏移ωx。通过叠前深度偏移得到偏移距域共成像道集,拾取剩余深度;在叠加成像剖面上拾取层位获取地层倾角;利用Runge-Kutta射线追踪计算地层反射点到检波点的射线路径,确定入射角,再计算地层反射点到复杂地表的射线路径;将偏移距域剩余深度映射为角度域,建立层析方程并求解。经模型和塔里木盆地"双复杂"地区Walkaway VSP实测资料验证表明,推导的层析反演方程能建立较准确的速度模型,偏移算法能适应复杂地表。 相似文献
10.
初至波走时层析是获取近地表速度结构的一种常用方法.随着采集技术的不断发展,可使用的数据量迅速增多,传统的基于射线追踪和解方程组的地震走时层析成像方法面临着内存占用大、方程求解不稳定等问题.为了解决这些问题,本文基于前人在波形反演研究中提出的一种改进的散射积分算法,提出了一种预条件最速下降法初至波走时层析.该方法无需存储核函数矩阵与Hessian矩阵即可方便地实现目标函数梯度的计算与预条件,且该方法计算效率高、求解稳定、易于并行.数值实验结果表明,该方法可以获得与传统方法精度相当的反演结果,但所占用的内存大幅减小. 相似文献
11.
12.
随着基于"真地表"叠前深度偏移技术的实施,偏移速度需要由浅表层到深层开始建立,而浅表层的速度一般是利用初至层析技术反演得到,但初至层析技术具有多解性和等效性,这种特性对于静校正量的计算是合理的,但对于叠前深度偏移的速度而言则是不合适的,它直接影响深度偏移射线路径的走向及深层速度迭代的精度.为了提高浅表层速度应用的合理性,本文开展了详细的研究,通过研究发现:参与初至层析的炮检距范围是影响初至反演近地表速度多解性的主要原因,为了选取合适的能够用于叠前深度偏移的近地表速度模型,根据数值模型分析的结果,形成了基于优化炮检距的大炮初至层析近地表速度的选取及应用技术.通过应用该技术,选取的速度与真实速度比较接近,达到了提高近地表速度反演精度和应用合理性的目的,为叠前深度偏移近地表速度模型的合理选取提供了技术依据. 相似文献
13.
根据地震波传播的有限频理论,对于某个特定震相的观测信息,不仅射线路径上的点对该信息具有影响,射线领域上的其他点对接收信息也具有影响,这种影响可以用核函数来表达.本文基于波动方程的Born近似与Rytov近似,给出了非均匀介质情况下初至地震波振幅与走时菲涅尔体层析成像单频、带限层析核函数的计算方法.通过对均匀介质情况下初至地震波菲涅尔体层析成像核函数解析表达式的理论模型实验与分析,给出了不同维度振幅、走时单频菲涅尔体的空间分布范围,进而给出了带限菲涅尔体边界的确定方法.将本文的走时菲涅尔体层析成像理论应用于表层速度结构反演中,理论模型试验与实际资料处理结果表明,初至波菲涅尔体地震层析成像方法比传统的初至波射线层析成像理论具有更高的反演精度. 相似文献
14.
15.
16.
伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点,但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反演结果,而二阶方向的实现又比较困难且费时.在伴随状态法的基础上,将走时差替换为定值,再次进行反演,便可以得到类似于射线密度的矩阵,用该矩阵的逆可以方便地进行预条件.基于该方法,本文提出了一种简单易行的预条件伴随状态法初至波走时层析的实现方法.理论模型和实际资料处理结果都表明,该方法既保留了伴随状态法初至波走时层析的优点,又可以克服一阶方向的局限,获得良好的反演效果. 相似文献
17.
《地球物理学进展》2017,(1)
海底地层速度结构是识别海洋天然气水合物储层的直接依据,本文应用地震反射走时层析成像建立了海底地层速度模型.采用不规则网格对模型进行离散化,使速度单元与反射界面单元完全耦合;利用基于不规则单元波前扩展和走时插值的射线追踪方法,精确确定反射射线路径和反射波走时;在反演中同时使用先验约束、平滑约束、归一化和正则化技术,提高了层析反演的稳定性和结果的可靠性.对南海北部神狐海域SH2井附近的二维地震测线资料,利用多域人机交互法拾取了反射走时,用反射走时层析成像方法获得了SH2井附近含水合物目标区的速度结构.该反演结果与测井声波速度和钻探结果一致,其中的高速带对应水合物储层,表明反射走时层析成像能够有效地得到海洋天然气水合物储层的速度结构,为海洋天然气水合物储层识别提供依据. 相似文献
18.
《地球物理学进展》1990,(3)
本文将线性层析与Herglotz-Wiechert反演公式进行了比较.二维和三维速度结构的层析反演使用了沿射线的线积分,而且可以以Radon交换的形式写出.对球对称介质,Radon变换可简化为Abel交换.因而对于直的射线路径,走时的Abel交换正是一维球对称介质的层析算法.Herglotz-Wiechert,公式应用地震走时资料反演一维地球结构,且通过Abel交换用精确的射线轨进推导.这是有历史意义的,因为它意味若特定的类似于层析的一种方法从本世纪初就应用在地反学中了(Her目心,1907,,八“h口t,1910).用数值例子比较了He叼otz一Wiechert算法和沿宜射线的线性层析.由Her川侧tZ刃犷沮忱h.时算法于特定条件下(无低速层)可应用于非直射线路径,它与层析的联系在评价推广到布曲射线中的层析结果的非唯一性方面非常有用. 相似文献
19.
《地球物理学进展》2016,(4)
使用Zelt和Barton的方法,通过一个计算效率高的有限差分求解eikonal方程,正演计算走时和射线路径.使用最小二乘QR分解法,求解稀疏线性系统方程组.使用正则化层析反演,结合用户给定的最小的、最平坦和最平滑的扰动限制,每一个加权因子随深度变化.结合数据残差和模型粗糙度的最小化,为数据残差提供一个最平滑的近似模型.该反演方法为非线性反演,需要一个初始模型,在每一次迭代时,需要计算新的射线路径.使用二维初至走时数据,对某油田二维井间地震实际资料进行无限频率初至走时层析反演.将反演所得到的速度与井的测井速度曲线相比较,二者吻合程度较高,表明该反演方法所得速度的分辨率比较高.证实了二维无限频率初至走时层析反演可以为全波形反演提供一个分辨率较高的长波长速度模型,从而为全波形反演井间地震实际资料提供了一个比较可靠的初始速度模型. 相似文献
20.
地震初至波中包含着丰富的近地表速度结构信息,如何分阶段、分尺度地利用这些信息进行近地表速度建模是地震勘探中的一个关键问题.在速度反演的不同阶段,综合利用初至波中的不同信息(如走时、包络和波形等)进行联合反演,可以有效地降低反演对初始模型的依赖程度,提高近地表速度模型的反演精度.为此,本文提出了一种统一基于波动方程正演引擎的初至波多信息联合反演方法,该方法同时匹配观测和模拟的初至波走时、包络和波形信息.在不同的反演阶段选择不同的权重因子调节不同信息的权重,这样不仅降低了反演对初始速度模型的依赖,而且自然地实现了多尺度反演.在每轮反演迭代中,一次正演模拟的波场同时应用于初至波走时、包络和波形匹配,无需额外的射线追踪.同时,联合反演方法在一定程度上缓解了串联反演中目标函数漂移问题,提高了近地表速度建模的精度. 相似文献