共查询到20条相似文献,搜索用时 9 毫秒
1.
三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先, 在常规双线性插值公式推导过程中, 充分利用平面波双线性假设的结论, 获得了二元极小值超越方程的解析解, 进而推导出了准确的局部走时计算公式, 同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次, 引入迎风差分的基本思想, 提出迎风双线性插值的局部走时计算策略, 该计算策略能简化算法、提高效率且保证无条件稳定性;然后, 将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合, 获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后, 通过精度和效率分析检验了新算法的精度、效率和正确性, 并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性. 相似文献
2.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
3.
地震波走时信息在叠前偏移、叠前速度分析、地震层析成像、走时反演及地震定位等中都有重要应用.快速推进法因其理论完善、精确灵活,无条件稳定,近年来已在走时计算领域得到广泛应用.快速扫描法作为求解一阶非线性双曲型偏微分方程的高效方法,已在图像处理、计算机视图、控制论等领域得到有效应用,且在走时计算方面有所应用且展现了广泛的应用前景.本文介绍了两种方法的基本原理且(通过均匀介质模型、局部低速体模型和Marmousi模型)把两种方法做了详细对比.研究结果表明:1)基于逆风差分格式的快速推进法和快速扫描法对纵横向速度变化很大的不均匀介质依然有很好的稳定性和适用性,均可以准确地计算地震波初至走时;2)对于相同的模型和在相同的计算条件下,两种方法的精度相当,但快速扫描法所耗的CPU时间较快速推进法明显减少,效率显著提高. 相似文献
4.
本文使用最小二乘线性迭代反演方法对跨孔雷达直达波初至时数据进行反演,每次迭代过程中,用有限差分法求解走时程函方程,并用高精度快速推进方法(HAFMM)进行波前扩展,通过追踪波前避免了进行射线追踪.为了验证该方案,我们对三组合成数据进行了测试,分析了单位矩阵算子、一阶差分算子和拉普拉斯算子等三种不同模型参数加权算子对模型的约束和平滑效果;讨论了FMM和HAFMM对反演精度的影响;测试了LSQR,GMRES和BICGSTAB等三种矩阵反演算法的反演效果.此外,我们还对一组野外实测数据进行了反演,对比了基于本方案以及基于平直射线追踪和弯曲射线追踪的走时层析成像反演效果.对比分析结果表明,使用拉普拉斯算子和HAFMM进行反演能较好地进行目标体重建,而三种矩阵反演方法对反演效果的影响差别不大;并且通过对波前等时线图的分析可以定性地判断异常体的性质和位置;而在对实测数据目标体的重建上,本方案能达到甚至优于弯曲射线算法的重建效果. 相似文献
5.
Subsurface rocks (e.g. shale) may induce seismic anisotropy, such as transverse isotropy. Traveltime computation is an essential component of depth imaging and tomography in transversely isotropic media. It is natural to compute the traveltime using the wavefront marching method. However, tracking the 3D wavefront is expensive, especially in anisotropic media. Besides, the wavefront marching method usually computes the traveltime using the eikonal equation. However, the anisotropic eikonal equation is highly non‐linear and it is challenging to solve. To address these issues, we present a layer‐by‐layer wavefront marching method to compute the P‐wave traveltime in 3D transversely isotropic media. To simplify the wavefront tracking, it uses the traveltime of the previous depth as the boundary condition to compute that of the next depth based on the wavefront marching. A strategy of traveltime computation is designed to guarantee the causality of wave propagation. To avoid solving the non‐linear eikonal equation, it updates traveltime along the expanding wavefront by Fermat's principle. To compute the traveltime using Fermat's principle, an approximate group velocity with high accuracy in transversely isotropic media is adopted to describe the ray propagation. Numerical examples on 3D vertical transverse isotropy and tilted transverse isotropy models show that the proposed method computes the traveltime with high accuracy. It can find applications in modelling and depth migration. 相似文献
6.
三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题. 相似文献
7.
由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景. 相似文献
8.
为克服矩形网速度间断面及边界描述精度差的缺点,本文提出了一种基于波前最小走时单元的三角网射线追踪算法.该算法用三角单元剖分速度模型,把与波前面相邻、处于振动状态的三角单元称为波前单元,从源点所在单元开始构造由波前单元组成的波前单元域;1)选域中最小走时三角单元为当前单元;2)它先与相邻三角单元振动传递,再进行相邻单元各自内部节点的传递、并将处于振动状态的相邻单元扩展为波前单元;3)然后确认当前单元的振动状态、决定其在波前单元域中的去留,以此三步循环进行波前单元域的扩展及节点的最小走时和次级源搜索,直至波前单元域在剖分区域的消亡;最后利用各节点最小走时和次级源,通过最小走时搜索,拾取从接收点到源点的射线路径.数值模拟表明,该算法的模型剖分灵活性强,对速度间断面的描述准确,追踪结果精度高. 相似文献
9.
三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果. 相似文献
10.
地震体波走时层析成像是探测地球内部速度结构的重要方法之一。基于三维块状建模以及三角形拼接的界面描述方式,结合快速高效的逐段迭代射线追踪方法,获得三维复杂地质模型中的地震射线路径与走时信息,采用共轭梯度非线性反演算法,进行地震波走时反演。实验结果表明共轭梯度反演算法在三维层状模型中具有较高的有效性。 相似文献
11.
波前构建法是一种能够快速计算地震波走时、射线路径及振幅的算法,该算法考虑从整个波场发出的射线族的追踪计算,为几何射线理论在勘探地震学中的应用带来了决定性的变化,使得能够同时计算成千上万务射线以及附在上面的走时和振幅成为可能.本文在收集、整理国内外研究现状的基础上,对波前构建法目前的应用介质模型、求解射线追踪方程组的算法... 相似文献
12.
线性走时插值法(LTI)在走时的计算中,由于射线方向考虑不全,计算得到的节点走时不一定最小,导致追踪的射线路径无法满足最小走时.针对这一问题,本文提出了一种改进的射线追踪算法,通过采用多方向的循环计算,得到所有计算节点的最小走时,使追踪到的射线路径能真正满足最小走时,以确保射线追踪的精度.模拟实验结果表明,在介质速度变化剧烈的结构中,该算法与传统的LTI算法相比,有效地提高了射线追踪的精度. 相似文献
13.
伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点.但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反演结果,而二阶方向的实现又比较困难且费时.在伴随状态法的基础上,将走时差替换为定值,再次进行反演,便可以得到类似于射线密度的矩阵,用该矩阵的逆可以方便地进行预条件.基于该方法,本文提出了一种简单易行的预条件伴随状态法初至波走时层析的实现方法.理论模型和实际资料处理结果都表明,该方法既保留了伴随状态法初至波走时层析的优点,又可以克服一阶方向的局限,获得良好的反演效果. 相似文献
14.
Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is required to sample the interfaces that are undulating in depth direction, such as the Moho and the 660 km discontinuity.Here we propose a new finite difference algorithm to solve the eikonal equation on non-orthogonal grid(irregular grid).To demonstrate its efficiency and accuracy, a test was conducted with a two-layer model. The test result suggests that the similar accuracy of a regular grid with ten times grids could achieve with our new algorithm, but the time cost is only about 0.1 times. A spherical earth model with an undulant660 km discontinuity was constructed to demonstrate the potential application of our new method. In that case, the traveltime curve fluctuation corresponds to topography. Our new algorithm is efficient in solving the first arrival times of waves associated with undulant interfaces. 相似文献
15.
目前,有关伴随状态法初至波走时层析成像方法的文献,基本上都是基于面积分来定义目标函数,由此得到的伴随方程也都依赖于地表的法向量.这样,一方面会因为伴随变量计算的不准确而造成梯度的不合理,另一方面也无法合理地处理井中观测问题.本文从理论或数值试验角度指出了这些问题,并提出了不依赖地表法向量的改进的伴随状态法走时层析成像方法.主要改进包括:(1)采用体积分定义目标函数,避免了传统方法不能较好处理井中观测数据的缺陷,可以适应任意地表或井中观测系统.(2)采用摄动法得到了新的伴随方程,克服了传统方法中伴随场计算需要依赖于地表法向量的缺陷,使得检波点处的走时残差可以正确地反传播至地下,进而得到更加合理的速度修正方向,提高了速度反演的精度. 相似文献
16.
Kirchhoff 3D prestack migration, as part of its execution, usually requires repeated access to a large traveltime table data base. Access to this data base implies either a memory intensive or I/O bounded solution to the storage problem. Proper compression of the traveltime table allows efficient 3D prestack migration without relying on the usually slow access to the computer hard drive. Such compression also allows for faster access to desirable parts of the traveltime table. Compression is applied to the traveltime field for each source location on the surface on a regular grid using 3D Chebyshev polynomial or cosine transforms of the traveltime field represented in the spherical coordinates or the Celerity domain. We obtain practical compression levels up to and exceeding 20 to 1. In fact, because of the smaller size traveltime table, we obtain exceptional traveltime extraction speed during migration that exceeds conventional methods. Additional features of the compression include better interpolation of traveltime tables and more stable estimates of amplitudes from traveltime curvatures. Further compression is achieved using bit encoding, by representing compression parameters values with fewer bits. 相似文献
17.
Accurate calculations of travel times and raypaths of reflection waves are important for reflection travel time tomography. The multistage shortest path method(MSPM) and multistage fast marching method(MFMM) have been widely used in reflection wave raytracing, and both of them are characterized by high efficiency and accuracy. However, the MSPM does not strictly follow Snell's law at the interface because it treats the interface point as a sub-source, resulting in a decrease in accuracy. The MFMM achieves high accuracy by solving the Eikonal equation in local triangular mesh. However, the implementation process is complex. Here we propose a new method which uses linear interpolation to compute the incident travel time of interface points and then using Snell's law to compute the reflection travel time of grid points just above the interface. Our new method is much simpler than the MFMM; furthermore, numerical simulations show that the accuracy of the MFMM and our new method are basically the same, thus the reflection tomography algorithms which use our new method are easier to implement without decreasing accuracy. Besides, our new method can be extended easily to other grid-based raytracing methods. 相似文献
18.
针对波前构建法射线追踪中非规则的波前四边形网格与规则的矩形网格节点之间的相对定位问题,采用矢量叉积判断方法从四个方面较详细的讨论了规则的矩形网格节点与非规则的波前四边形的相对位置关系,并针对四种不同的位置关系采用不同的插值算法进行其属性的转换计算,误差分析及计算实例表明本文所采用的网格定位方法以及属性转换的插值算法从理论到实际应用都是正确的、切实可行的. 相似文献
19.
We present a new method of three-dimensional (3-D) seismic ray tracing, based on an improvement to the linear traveltime interpolation (LTI) ray tracing algorithm. This new technique involves two separate steps. The first involves a forward calculation based on the LTI method and the dynamic successive partitioning scheme, which is applied to calculate traveltimes on cell boundaries and assumes a wavefront that expands from the source to all grid nodes in the computational domain. We locate several dynamic successive partition points on a cell's surface, the traveltimes of which can be calculated by linear interpolation between the vertices of the cell's boundary. The second is a backward step that uses Fermat's principle and the fact that the ray path is always perpendicular to the wavefront and follows the negative traveltime gradient. In this process, the first-arriving ray path can be traced from the receiver to the source along the negative traveltime gradient, which can be calculated by reconstructing the continuous traveltime field with cubic B-spline interpolation. This new 3-D ray tracing method is compared with the LTI method and the shortest path method (SPM) through a number of numerical experiments. These comparisons show obvious improvements to computed traveltimes and ray paths, both in precision and computational efficiency. 相似文献
20.
通过引入模糊均值聚类(FCM)模型约束函数对电阻率与速度进行约束,开展二维直流电阻率法与地震初至波走时成像法联合反演研究.在地下浅层结构勘探中,通常低电阻率的地质体具有低速特征,较高电阻率的地质体表现为较高的地震波速度.直流电阻率法因为低电阻率区域吸引电流而对其敏感,地震走时成像法因为射线集中在高波速区而对高速体敏感,因此,两者联合成像能够大幅度提高反演效果.合成数据反演表明,直流电阻率法和地震初至波走时联合反演对于两类地质体的分辨能力均有提升,能够优势互补.尤其是引入FCM模型约束进行联合反演,根据已知物性进行监督学习,进一步提高了反演质量,改善了成像模型的分辨率. 相似文献
|