共查询到16条相似文献，搜索用时 171 毫秒
1.
本文使用改进后的最短路径算法（MSPM）结合分区多步计算技术实现了三维复杂层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应走时和射线路径的跟踪计算.其原理是将三维复杂层状模型按速度界面分成若干个独立的计算区域,采用分步计算技术进行多次波的跟踪计算.基于多次波是通过速度界面简单的入射、透射、反射及转换波按一定规律及原理的不同组合,因此可实施分区多步计算技术.数值模拟实例及误差分析表明分区多步计算技术具有单步最短路径算法中的诸多优点:算法简单、数值计算稳健、计算精度高、速度快及全球解等,因此是解决多次波跟踪计算行之有效的方法. 相似文献
2.
快速行进法(FMM)是一种求解程函方程数值解计算网格节点走时,然后向后处理进行射线追踪的方法.为了求取任意起伏界面下高精度多震相的地震走时与相应的射线路径,本文采用任意起伏地表条件下的的三维不等距上行差分公式结合分区多步计算技术实现了三维复杂层状起伏介质中多震相(透射、反射、转换波)地震走时的计算,利用上行有限差分公式逐次进行射线路径的追踪,并且通过与较为成熟的不规则最短路径法(ISPM)对比,验证了本算法的计算精度和有效性.数值模拟实例和对比结果表明该算法具有较高的计算精度,数值计算稳健,能灵活处理含任意三维起伏界面模型中多震相地震走时及相应射线路径的追踪问题. 相似文献
3.
在改进后的最短路径算法（MSPM）中引入分区多步计算技术实现了二维层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应的走时和射线路径的跟踪计算.其原理是将二维复杂层状模型按速度界面分成若干个独立的计算区域,采用分步计算技术进行多次波的跟踪计算.基于多次波是通过速度界面的简单入射、透射、反射及转换波按一定规律的不同组合,因此可实施分区多步计算技术.通过某一上、下层界面的透射（或透射转换）波实际上是由上层得到的下行波加上由该界面透射的下行波组成,若为转换波则使用不同的速度模型;而经过某一界面的反射（或反射转换）波实际上是由某层内计算得到的下行（或上行）波再加上由该界面反射的上行（或下行）波组成.这样即可得到分区独立计算,并通过速度界面分步连接达到跟踪多次波的目的.计算结果表明MSPM算法下的分区多步计算技术具有单步SPM算法中的诸多优点,即:算法简单、数值计算稳键、计算精度高、速度快及全球解等,因而是解决多次波跟踪计算行之有效的方法. 相似文献
4.
5.
在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度. 相似文献
6.
《地球物理学进展》2017,(3)
为了解决复杂速度模型中的走时正、反演问题,例如:含不规则起伏地表、不规则地下波阻抗界面、以及不规则速度异常体的复杂地学模型,本文采用三角网格单元模型参数化下的分区多步改进型最短路径算法,实现了多震相地震射线的追踪计算,结合共轭梯度法求解带约束的阻尼最小二乘反演问题,实现了多震相走时联合同时反演成像的方法技术.当界面起伏较大时出现散射,从而造成散射点所在区域射线密度过密,导致该区域内速度和界面的过度更新.为了克服上述问题,我们在同时反演中引入了射线密度的概念,从而有效地解决了上述过度更新问题.数值模拟实验表明:采用三角网格单元进行模型参数化,可保证在复杂模型中的正演计算具有较高的计算精度;同时反演中可以准确地刻画不规则异常体和不规则反射界面.因此,本文提出的走时成像方法技术具有较广的实用价值. 相似文献
7.
8.
球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系下分区多步不规则最短路径算法推广至球坐标系中,进行区域或者全球尺度的多震相射线追踪.然后将其与适合多参数同时反演的子空间算法相结合,形成一种球坐标系下联合多震相走时三参数同时反演的方法技术.与双参数(速度和反射界面或速度和震源位置)同时反演的数值模拟对比分析显示:三参数与双参数的同时反演结果大体接近,并且它们对到时数据中可容许的随机噪声不太敏感.结果说明本文中的同时反演成像为一种提高成像分辨率,同时反演速度、震源位置和反射界面的有效方法. 相似文献
9.
本文针对三维复杂介质的两点射线追踪的适用性和有效性,提出了网格模型和层状模型相结合的建模方式,文中利用规则网格化的层状模型对复杂地质单元进行数值描述,一方面避免了用纯粹复杂地质网格体描述时需要的大量内存,另一方面结合运用规则化的层状模型,必有利于提高计算效率,在计算单元自动生成的射线追踪方法基础上,本文结合波前法提出一种快速、有效的两点射线追踪方法,利用反射波波前在地表形成的射线方位角和倾角坐标体系,能快速、准确地收敛到所有从源点到接收点的射线路径,亦可有效地解决射线路径的多值解难题,从而可适用于复杂地质模型的特点. 相似文献
10.
对于微地震正演模拟,本文以射线追踪的原理为基础,对两点间的射线追踪问题进行了研究,应用二分算法、改进二分算法和微变网格算法对水平层状匀速模型、弯曲层状匀速模型和复杂地质模型进行射线追踪,使得计算效率和适用范围都得到了很大的改善.文中对每种算法误差范围和计算效率进行了对比验证,对于不同的地质模型,选用合适的算法才能在计算速度和精度上得到双重保证,最后正演模拟了多波三分量记录.在模型建立上引入了超薄层概念,并在前人模拟的直达波、透射波、反射波基础上拟了折射波,使正演模拟的多波场信息更丰富.文中的应用实例及模型结果表明:与二分法相比,改进二分法能够对弯曲界面进行射线追踪,并能保证结果的精度.弯曲层状模型中,改进二分法与微变网格法相比计算速度有显著提高,能够应用到资料的反演中. 相似文献
11.
3-D simultaneous inversion for velocity and reflector geometry using multiple classes of arrivals in a spherical coordinate frame 总被引：1，自引：0，他引：1
Guo-jiao Huang Chao-ying Bai Xing-wang Li Stewart Greenhalgh 《Journal of Seismology》2014,18(1):123-135
The two key requirements in conducting 3-D simultaneous traveltime tomography on real data at the regional and global scale with multiple classes of arrival time information are (1) it needs an efficient and accurate arrival tracking algorithm for multiply transmitted, reflected (or refracted) and converted waves in a 3-D variable velocity model with embedded velocity discontinuities (or subsurface interfaces), and (2) a subdimensional inversion solver is required which can easily search for different types of model parameters to balance the trade-off between the different types of model parameter updated in the simultaneous inversion process. For these purposes, we first extend a popular grid/cell-based wavefront expanding ray tracing algorithm (the multistage irregular shortest-path ray tracing method), which previously worked only in Cartesian coordinate at the local scale, to spherical coordinates appropriate to the regional or global scale. We then incorporated a fashionable inversion solver (the subspace method) to formulate a simultaneous inversion algorithm, in which the multiple classes of arrivals (including direct and reflected arrivals from different velocity discontinuities) can be used to simultaneously update both the velocity fields and the reflector geometries. Numerical tests indicate that the new inversion method is both applicable and flexible in terms of computational efficiency and solution accuracy, and is not sensitive to a modest level of noise in the traveltime data. It offers several potential benefits over existing schemes for real data seismic imaging. 相似文献
12.
Seismic wavefront evolution of multiply reflected, transmitted, and converted phases in 2D/3D triangular cell model 总被引：1，自引：0，他引：1
Conventionally grid-cell-based schemes for simulating seismic wavefront propagation, such as the finite difference eikonal equation solver or the shortest-path method, usually adopt regular grids or cells in model parameterization to obtaining (but not exclusively) first arrivals only. However, later arrivals, which often result from the velocity interfaces or discontinuities, can be prevalent and significant (sometimes of large amplitude), making them potentially important additional information to use in practical applications. To better approximate the data acquisition geometry and the irregular interfaces, we exploit a triangular shortest-path method (TSPM; that is to use triangular cells in model parameterization) to simulate seismic wavefront evolution, comprising any kind of transmissions, reflections (or refractions), mode conversions, and combinations thereof, in 2D/3D heterogeneous media. A practical procedure, known as the multistage scheme, was incorporated with the TSPM to propagate seismic wavefronts from one interface (or subsurface in 3D) to the next. By treating each separate layer that the wavefront enters as an independent computational domain, one can simulate wavefront transmission and mode conversion by reinitializing it in the adjacent layer and wavefront reflection (and/or conversion) by reinitializing it in the incident layer. To further improve the computational accuracy, a second level of forward star scheme, previously defined in the grid model, is introduced into the triangular cell model. Several examples (including the Marmousi model) are used to demonstrate the viability and versatility of the multistage TSPM in heterogeneous media, even in the presence of high-velocity contrasts involving interfaces of relatively high curvature. With the introduction of the second level of forward star scheme, the total numbers of nodes are reduced sufficiently, and hereafter the computer memory is less required. Most important is that the computing accuracy with the second-level forward star scheme can be largely improved over those with the first level of forward star scheme applied in the multistage TSPM scheme. 相似文献
13.
三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题. 相似文献
14.
Hui Sun Jian-Guo Sun Zhang-Qing Sun Fu-Xing Han Zhi-Qiang Liu Ming-Chen Liu Zheng-Hui Gao Xiu-Lin Shi 《应用地球物理》2017,14(1):56-63
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models. 相似文献
15.
针对处理起伏地表(或含地下不规则波阻抗界面)条件下发展起来的地震波场数值模拟算法的模拟结果与解析解(大多数情形下无法得到)无法进行对比,且其有效性和正确性难以验证的情况,本文提出了一种可以相互验证波场数值模拟结果与射线追踪数值模拟结果的正确性和有效性的佐证方法,验证了参考射线追踪法.其中,波场数值模拟中采用曲线网格DRP/opt MacCormack有限差分法,射线追踪模拟则采用分区多步三角网格最短路径算法.通过系统对比上述两种方法得到的波场快照、单炮地震记录,以及合成理论地震图的结果显示,本方法相互作证了两种方法所得结果的正确性和有效性.双层和三层起伏层状模型的对比分析结果表明,这种方法不但可以加深理解地震波在复杂介质中的传播规律,同时射线法的引入为清晰识别和标定地震波场数值模拟中各种不同震相提供了一种便捷的途径. 相似文献
16.
传统有限差分方法在处理起伏地表时存在一些困难,而坐标变换法可将起伏地表映射为水平地表以克服此缺点.但同时,地下构造被变换得更加复杂,导致了波传播和成像的不准确.本文提出了一种分层的坐标变换方法,并应用到了弹性波逆时偏移中,此方法既可以克服起伏地表的影响,又可以不破坏地下构造.波场正向延拓、逆时延拓和分离是在辅助坐标系下完成的,而成像是在笛卡尔坐标系下完成的.通过对简单起伏模型和中原起伏模型的试算证明了本文提出方法的准确性.同时,对两种极端起伏地层高程不准确的情况进行测试可以看出:分层坐标变换逆时偏移方法的成像效果远好于传统坐标变换方法. 相似文献