首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 44 毫秒
1.
IndirectapproachmethodforspecifiedendpointsseismicraytracinginthredimensionalinhomogeneousmediaCHAOFANXU(徐朝繁)XIANKANGZH...  相似文献   

2.
二步法射线追踪   总被引:12,自引:2,他引:12  
马争鸣  李衍达 《地球物理学报》1991,34(4):501-508,T001
本文提出的二步法射线追踪,能够在复杂地质模型下有效、准确地确定射线在炮点的出射角,使得射线恰好到达检波点.射线到达地面的位置是射线在炮点的出射角的函数,称为出射函数.二步法射线追踪通过搜索出射函数的极值点来确定其单调变化的区间;然后通过比较检波点位置与出射函数极值的大小,确定与检波点相对应的出射角究竟落在出射函数的哪一个单调变化的区间里;最后,利用出射函数在局部区域里单调变化的性质,迅速、准确地确定与检波点相对应的出射角.出射函数的极值点,对于共炮集中每一个检波点和检波点中每一个不同方向的入射射线都是适用的,因此,二步法射线追踪只需要进行一次一维全局搜索(极值点搜索).另外,本文采用B-样条的线性组合来表示速度函数,通过提高B-样条的阶数,可以改善速度函数的光滑程度,从而既满足射线追踪对于速度函数光滑性的要求,又加强射线在炮点的出射角对于射线的控制作用.  相似文献   

3.
对于微地震正演模拟,本文以射线追踪的原理为基础,对两点间的射线追踪问题进行了研究,应用二分算法、改进二分算法和微变网格算法对水平层状匀速模型、弯曲层状匀速模型和复杂地质模型进行射线追踪,使得计算效率和适用范围都得到了很大的改善.文中对每种算法误差范围和计算效率进行了对比验证,对于不同的地质模型,选用合适的算法才能在计算速度和精度上得到双重保证,最后正演模拟了多波三分量记录.在模型建立上引入了超薄层概念,并在前人模拟的直达波、透射波、反射波基础上拟了折射波,使正演模拟的多波场信息更丰富.文中的应用实例及模型结果表明:与二分法相比,改进二分法能够对弯曲界面进行射线追踪,并能保证结果的精度.弯曲层状模型中,改进二分法与微变网格法相比计算速度有显著提高,能够应用到资料的反演中.  相似文献   

4.
迭代法在三维网络模型最短路径射线追踪中的应用   总被引:2,自引:2,他引:0       下载免费PDF全文
目前有许多基于二维网格模型的射线追踪算法,但应用与三维网格模型的射线追踪算法较少,主要是因为三维网格射线追踪算法在算法实现、编程以及图形显示方面都有一定难度.本文提出了一种适用于三维网格射线追踪的迭代算法,并应用于几种理论模型,同时实现了三维绘图显示,通过研究证明迭代算法能在一定程度上提高三维网格射线追踪的精度和速度.  相似文献   

5.
TTI介质是石油地震勘探领域最常用的各向异性介质,快速计算TTI介质射线路径和走时信息有重要的研究意义.TTI介质传统运动学射线追踪方法一般基于任意弹性介质射线方程,利用Bond变换或者四阶张量变换来处理复杂的21个弹性参数,因而非常耗时.实际野外对称轴统一的TTI介质模型,一般可以看成VTI介质模型旋转一定角度获得.为此,本文推导了三维VTI介质射线追踪方程,提出先在本构坐标系中进行VTI介质射线追踪,再通过坐标旋转将射线路径旋转至观测坐标系中,获得TTI介质射线路径.数值模型计算表明该方法高效和精确,较传统方法效率提高了近4倍.在强各向异性等特殊情况下,体波波前面都与理论群速度面一致.  相似文献   

6.
复杂地质层面的网格化模型及两点射线追踪方法   总被引:3,自引:3,他引:0       下载免费PDF全文
本文针对三维复杂介质的两点射线追踪的适用性和有效性,提出了网格模型和层状模型相结合的建模方式,文中利用规则网格化的层状模型对复杂地质单元进行数值描述,一方面避免了用纯粹复杂地质网格体描述时需要的大量内存,另一方面结合运用规则化的层状模型,必有利于提高计算效率,在计算单元自动生成的射线追踪方法基础上,本文结合波前法提出一种快速、有效的两点射线追踪方法,利用反射波波前在地表形成的射线方位角和倾角坐标体系,能快速、准确地收敛到所有从源点到接收点的射线路径,亦可有效地解决射线路径的多值解难题,从而可适用于复杂地质模型的特点.  相似文献   

7.
复杂山地随机介质GMM-ULTI法射线追踪   总被引:2,自引:2,他引:0       下载免费PDF全文
对复杂山地介质的非均质性以及介质中地震波运动学特征进行深入研究,对于提高复杂山地区域地震勘探的效果有着重要的理论意义和实际价值.为了研究复杂山地非均质性和该介质中地震波的一些运动特性,提出了一种复杂山地随机介质的建模方法和一种新的射线追踪算法.与常规算法相比,复杂山地随机介质的生成方法采用更贴近实际介质特点的梯度介质作为背景介质,并在模型生成过程中加入地形修正步骤;新提出的GMM-ULTI射线追踪算法,充分融合群推进法、迎风思想、走时插值法的优势,采用先计算走时后追踪射线路径的两步策略完成射线追踪.算法分析与计算实例表明:复杂山地随机介质的生成方法能灵活、精细且更贴近实际地刻画复杂山地介质的非均质特点;新射线追踪算法兼顾精度和效率、能无条件稳定且灵活地适应复杂山地随机介质的特点;同时基于对几个模型试算结果的分析也得出了复杂山地随机介质中的地震波的一些传播规律.  相似文献   

8.
横向各向同性介质中地震波走时模拟   总被引:15,自引:0,他引:15       下载免费PDF全文
横向各向同性介质是地球内部广泛分布的一种各向异性介质.针对这种介质,我们对各向同性介质的最小走时树走时模拟方法进行了推广,推广后的方法可适用于非均匀、对称轴任意倾斜的横向各向同性介质模型.为保证计算效率,最小走时树的构建采用了一种子波传播区域随地震波传播动态变化的改进算法.对于弱各向异性介质,我们使用了一种新的地震波群速度近似表示方法,该方法基于用射线角近似表示相角的思想,对3种地震波(qP, qSV和qSH)均有较好的精度.应用本文地震波走时模拟方法对均匀介质、横向非均匀介质模型进行了计算,并将后者结果与弹性波方程有限元方法的模拟结果进行了对比,结果表明两者符合得很好.本文方法可用于横向各向同性介质的深度偏移及地震层析成像的深入研究.  相似文献   

9.
A fast and robust method for two-point ray tracing in one-dimensional layered media is presented. This method is applicable to layered models with constant or linearly varying isotropic layer velocity. For given model properties and source and receiver positions, a ray path can be uniquely determined once its ray parameter (i.e. horizontal slowness) is known. The ray parameter can be obtained by numerically solving the nonlinear offset (i.e. source–receiver horizontal distance) equation using Newton's method, which generally works well at near and mid offsets. However, Newton's method becomes hard to converge at large offsets due to the oversensitivity of offset to ray parameter. Based on the analysis of the characteristic of the offset equation, a modified ray parameter is proposed and used to replace the generic ray parameter in numerical calculation. Numerical experiments show that the iteration process becomes stable and converges rapidly with the modified ray parameter. Moreover, a rational function that asymptotically approximates the shape of the offset equation is introduced for obtaining good initial estimates of the modified ray parameter. Numerical tests show that this method is robust in any situation, and an accurate ray parameter can be obtained within two or three iterations for a wide range of model velocity structure and source–receiver distance. Furthermore, the proposed two-point ray tracing method is easy to implement.  相似文献   

10.
The stable-beam method for forward modelling of seismic data is introduced. The method is applicable to geometries which may be approximated by a series of single-valued (in depth) interfaces separating constant-velocity layers. For models of this restricted type, the results are of similar accuracy to those from waveequation-based methods whilst the run times are similar to, or better than, those for simple ray-tracing approaches. The basis of the method is to approximate interfaces by a series of straight-line segments. This allows very rapid and stable ray tracing through the model. Pseudodiffractions are then added from all of the interface discontinuities formed between adjacent segments. These pseudodiffractions have the effect of correcting for the phase, amplitude and wavefront continuity errors introduced by the interface approximation. Comparison of the stable-beam results to analytical, Kirchhoff, finite-difference and physical model results confirm the accuracy ofthe technique.  相似文献   

11.
I introduce a new explicit form of vertical seismic profile (VSP) traveltime approximation for a 2D model with non‐horizontal boundaries and anisotropic layers. The goal of the new approximation is to dramatically decrease the cost of time calculations by reducing the number of calculated rays in a complex multi‐layered anisotropic model for VSP walkaway data with many sources. This traveltime approximation extends the generalized moveout approximation proposed by Fomel and Stovas. The new equation is designed for borehole seismic geometry where the receivers are placed in a well while the sources are on the surface. For this, the time‐offset function is presented as a sum of odd and even functions. Coefficients in this approximation are determined by calculating the traveltime and its first‐ and second‐order derivatives at five specific rays. Once these coefficients are determined, the traveltimes at other rays are calculated by this approximation. Testing this new approximation on a 2D anisotropic model with dipping boundaries shows its very high accuracy for offsets three times the reflector depths. The new approximation can be used for 2D anisotropic models with tilted symmetry axes for practical VSP geometry calculations. The new explicit approximation eliminates the need of massive ray tracing in a complicated velocity model for multi‐source VSP surveys. This method is designed not for NMO correction but for replacing conventional ray tracing for time calculations.  相似文献   

12.
In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.  相似文献   

13.
分块三维速度模型生成及理论地震图的计算   总被引:2,自引:0,他引:2       下载免费PDF全文
本文提出了在计算机上实现分块三维地壳模型及利用加权最小二乘拟合生成平缓光滑的三维速度函数的方法,给出了适用于分块、块内速度连续变化的三维模型中Cauchy射线追踪的新算法,简介了基于上述方法反射线的基本理论所编制的合成三维理论地震图的程序包RSSGTD.给出的两个盆地状模型的算例表明,所使用的模型生成方法具有模拟复杂地壳结构的能力;与三维样条函数方法比较,最小二乘拟合方法能给出更加适合射线方法合成地震图计算的速度函数,并且内存小、计算速度快;所给出的Cauchy射线追踪算法能够适合块状模型中任何体波射线的追踪.  相似文献   

14.
Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the coupling ray theory for S waves. The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known, are summarized in a form which represents a complete algorithm suitable for coding and numerical applications.  相似文献   

15.
跨孔观测地震数据的速度重建   总被引:2,自引:0,他引:2  
利用矢量射线追踪正演模拟技术 ,计算地震波传播的路径及走时 ,进而利用射线走时及路径的内插 ,发展了弯曲射线迭代反演技术 .该方法可用来重建井间地层的速度图像 .基于所发展的方法 ,我们对两种较为复杂的典型地质模型进行了井间速度重建 .结果表明该方法是一项快速、高精度的跨孔数据速度重建技术  相似文献   

16.
三维复杂介质的块状建模和试射射线追踪   总被引:33,自引:9,他引:24       下载免费PDF全文
为了解决三维复杂介质的射线追踪,本文改变了传统的三维层状地层的建模描述方式,提出了块状结构的建模描述方法,结合三角形面片来描述地质界面,可以构造非常复杂的三维地质模型.为了满足射线追踪的需要,本文对模型界面内的法向量进行光滑处理,光滑后的法向量在界面内是连续变化的.在块状模型的基础上,本文运用三角形的面积坐标,提出了几种试射角度的修正方法:细分三角形法、分割三角形法和子三角形法,计算表明子三角形法最好.文中给出了三维块状模型和射线追踪实例.  相似文献   

17.
复杂介质小波多尺度井间地震层析成像方法研究   总被引:15,自引:1,他引:15       下载免费PDF全文
复杂介质井间地震层析成像是一个很复杂的非线性反 演问题,常规的线性化反演方法无法得到好的解. 采用基于图形的弯曲射线追踪方法, 并将小波多尺度思想引入到井间层析成像,提出了小波多尺度井间地震层析成像方法,很好 地解决了非线性成像的难题,提高了图像的质量和分辨率. 物理模型实验结果表明, 该方法适合于复杂介质成像,并具有良好的实用性和效果.  相似文献   

18.
起伏地表条件下各向异性地震波最短路径射线追踪   总被引:2,自引:1,他引:1       下载免费PDF全文
在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.  相似文献   

19.
In this paper, a global inversion method is developed for seismic moment tensor inversion by using the body wave forms. The algorithm depends on neither the selection of starting model nor the forms of objective function and constraints. When the error function, measure of the difference between the observed and synthetic waveforms, is chosen as the objective function, the best fitting source model is found; when a certain combination in seismic moment tensor elements is selected as the objective function and the values of error function are constrained in a suitable bound, the extreme source models can be produced by minimizing or maximizing this combination. By changing the form of the combination of moment tensor elements, a variety of different source characteristics can be considered. Therefore the extreme solution provides an estimation of the uncertainty in the best fitting source model. The seismic waveform data was used to evaluate the effectiveness of this algorithm. This research was supported by the National Natural Science Foundation of China.  相似文献   

20.
三维非均匀地质模型中的逐段迭代射线追踪   总被引:8,自引:4,他引:4       下载免费PDF全文
地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性.  相似文献   

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

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