共查询到20条相似文献，搜索用时 160 毫秒
1.
快速行进法(FMM)是一种求解程函方程数值解计算网格节点走时,然后向后处理进行射线追踪的方法.为了求取任意起伏界面下高精度多震相的地震走时与相应的射线路径,本文采用任意起伏地表条件下的的三维不等距上行差分公式结合分区多步计算技术实现了三维复杂层状起伏介质中多震相(透射、反射、转换波)地震走时的计算,利用上行有限差分公式逐次进行射线路径的追踪,并且通过与较为成熟的不规则最短路径法(ISPM)对比,验证了本算法的计算精度和有效性.数值模拟实例和对比结果表明该算法具有较高的计算精度,数值计算稳健,能灵活处理含任意三维起伏界面模型中多震相地震走时及相应射线路径的追踪问题. 相似文献
2.
3.
地震直达波走时层析成像可归结为求解一个大型的、稀疏的、常常是病态的线性方程组.求解方程组常用的迭代法,需要一个比较合理的初始猜测解,也即是初始速度模型.初始模型关系到反演的效率甚至成像的正确性.本文在前人研究基础上提出一种生成模型网格节点初始速度方法,假定震源到检波点路径为直线,记录每条射线穿过的单元和统计每个网格单元穿过的射线数目、自动拾取网格节点所在单元的数目等.实例中,由原始数学模型的正演旅行时资料生成节点初始速度模型,效果可以.最后,分别采用均匀模型和本文方法生成的初始模型进行迭代反演,通过比较,证实该自动生成节点初始模型的可行性和可靠性,并对存在的问题进行讨论和解释. 相似文献
4.
人工源宽角反射(折射)地震资料具有偏移距较大、信噪比较低等特点,通常用于地震波走时反演重建地壳速度结构。逆时偏移成像方法作为勘探地震学领域获取地下构造形态的有效手段之一,可以有效弥补走时反演方法的不足。本文针对大偏移距宽角反射(折射)地震实验,利用四边形网格谱元法进行波场模拟,结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度获取模型合成地震记录,后采用逆时偏移成像方法将合成地震记录偏移归位,获取地壳几何结构,验证了逆时偏移成像方法在宽角地震资料处理及结果解释中的适用性,为后期实际地震资料的偏移成像提供了理论依据和支持。 相似文献
5.
射线追踪方法是研究地震波在横向非均匀地壳介质中传播的重要方法.本文推导了理论走时对网格化节点速度的偏导数公式,提出了针对深地震宽角反射/折射剖面数据反演的联合迭代法,并使用该方法对横跨龙门山断裂带中段的一条深地震宽角反射/折射剖面进行了反演和解释.首先,对每一炮的观测数据进行一维反演,在此基础上插值出一个粗略的二维速度模型;然后,使用射线追踪方法计算理论走时,再根据理论走时与观测走时的拟合程度对二维模型进行调整,以获得更加接近实际的二维速度模型;最后,利用联合迭代法对观测走时进行反演,经反复迭代使所有接收点理论走时与实测走时的残差平方和最小,最终获取该剖面的二维地壳速度结构.反演结果表明:测线东段的沉积盖层明显厚于中段褶皱带和西部高原,中部褶皱带部分地区出现基岩裸露;构造转换带两侧的地层分界面近于水平层状分布,其西侧的中、下地壳内各存在一个层间速度间断面;构造转换带内存在薄厚不等的低速层,自西向东有增厚趋势.此外,龙门山断裂带的3条主断裂向下深切结晶基底,这是由于西部松潘—甘孜地块自西向东运动,受到刚性扬子地块的阻挡,沿铲式断裂向上爬升所致;而在断层上盘距地表约15 km深处出现的最大剪应力极值区,正是发生汶川MS8.0地震的震源位置. 相似文献
6.
最短路径射线追踪方法是计算地震波走时的主要方法之一,该方法基于惠更斯原理和费玛原理,具有稳健、适于复杂介质模型的优点.为处理方便,最短路径方法中的介质模型通常以规则网格进行剖分,界面节点(界面与网格的交点)以其邻近的模型单元节点(即边界单元节点)近似表示.界面近似将导致计算误差,对于反射波尤为严重.反射波的走时精度可通过减小网格的尺寸提高,但这样会大大增加计算时间,为高精度和高效率地计算地震反射波走时,我们提出了一种基于规则网格的走时校正技术.地震波传播至或起始于边界单元节点的走时校正为地震波传播至或起始于该边界单元节点所对应的界面节点的走时.数值模型计算结果表明,走时校正方法可使反射波的走时精度提高约1～2个数量级,而其计算时间则和常规算法基本上在相同量级. 相似文献
7.
8.
在地震定位中常常需要求解震源轨迹,但由于复杂介质中的震源轨迹较为复杂,难以给出其解析解,因此震源轨迹的计算通常仅限于简单介质模型.本文基于最小走时树射线追踪技术,提出了一种计算复杂介质中震源轨迹的方法.为回避发震时间问题,以观测到时差作为震源轨迹的约束条件.首先从模型节点中选出少量理论到时差与观测到时差之绝对差,即双重时差较小的点作为震源轨迹的代表点,然后以其中双重时差最小的点为初始点,在双重时差场中利用最小走时树射线追踪方法计算出初始点到其他震源轨迹代表点的射线路径作为震源轨迹.当选的震源轨迹代表点较多时,得到的震源轨迹较为粗略,此时可去掉射线经过次数较少的代表点的射线路径使震源轨迹更为精细.为减少计算量,对最小走时树射线追踪方法的终止条件做了修正.以一个复杂介质模型中的地震为例,计算了包括速度扰动、到时扰动等不同情况下的震源轨迹,结果表明所提出的震源轨迹计算方法切实可行. 相似文献
9.
本文采用Ｈａｍｉｌｔｏｎ理论方法,研究反射地震记录域中地下同一反射点对应的反射走时随炮检距变化的特性,这在叠前地震资料处理中具有重要意义．由于走时函数可以表示成炮点位置、检波点位置、地震波射线出射角和接收角的函数,通过对共反射点走时与炮点位置变化关系的分析,提出了复杂介质中共反射点轨迹可用Ｈａｍｉｌｔｏｎ正则方程描述．在线性变速介质中,走时函数可以解析给出,由此可导出该介质中共反射点轨迹的Ｈａｍｉｌｔｏｎ正则方程．文中结合常速度、横向变速和线性变速介质模型,计算了地下不同反射面上各点相应的共反射点轨迹,阐明Ｈａｍｉｌｔｏｎ方法的一些特点,并首次引入Ｈａｍｉｌｔｏｎ理论方法描述共反射点走时随炮检距的变化特性,开创了该领域研究的一条新途径． 相似文献
10.
在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度. 相似文献
11.
12.
Wide-angle reflections are now routinely recorded in high resolution explosion seismics to study the crustal structure. Use of Dix's hyperbolic approximation to the nonhyperbolic wide-angle reflection travel times causes major errors in the determination of interval velocities and layer thicknesses of a stack of horizontal velocity layers. Here we propose a layer stripping method to directly calculate the interval velocities and layer thicknesses in a vertically heterogeneous earth from the strong and reliable wide-angle reflected events. Synthetic reflection travel times, at wide-angle range, for a given velocity model, contaminated by some random errors, have been used to demonstrate the reliability of the algorithms to determine the interval velocities and thicknesses of various layers. The method has also been tested on two field examples along two deep seismic sounding (DSS) profiles with well identified wide-angle reflection travel times, which illustrates the practical feasibility of the proposed method. 相似文献
13.
The coupling–ray–theory tensor Green function for electromagnetic waves or elastic S waves is frequency dependent, and is usually calculated for many frequencies. This frequency dependence represents no problem in calculating the Green function, but may pose a significant challenge in storing the Green function at the nodes of dense grids, typical for applications such as the Born approximation or non–linear source determination. Storing the Green function at the nodes of dense grids for too many frequencies may be impractical or even unrealistic. We have already proposed the approximation of the coupling–ray–theory tensor Green function, in the vicinity of a given prevailing frequency, by two coupling–ray–theory dyadic Green functions described by their coupling–ray–theory travel times and their coupling–ray–theory amplitudes. The above mentioned prevailing–frequency approximation of the coupling ray theory enables us to interpolate the coupling–ray–theory dyadic Green functions within ray cells, and to calculate them at the nodes of dense grids. For the interpolation within ray cells, we need to separate the pairs of prevailing–frequency coupling–ray–theory dyadic Green functions so that both the first Green function and the second Green function are continuous along rays and within ray cells. We describe the current progress in this field and outline the basic algorithms. The proposed method is equally applicable to both electromagnetic waves and elastic S waves. We demonstrate the preliminary numerical results using the coupling–ray–theory travel times of elastic S waves. 相似文献
14.
V. Meshbey E. Ragoza D. Kosloff U. Egozi T. Wexler 《Pure and Applied Geophysics》2002,159(7-8):1563-1582
— We present a new travel-time calculation method based on Fermat's principle. In the method, travel times are recursively calculated on horizontal planes of increasing depth. For typical configurations of exploration geophysics, the distance between the planes is on the order of 100 m. The travel times on the first plane are calculated by connecting straight ray segments from the source point to the grid points on the plane and integrating the slowness along each segment. The times on the other planes are calculated by finding the minimum of the combination of the times on the plane above plus the additional time along segments connecting grid points on the two planes. The travel-time calculation method is designed for calculating either the first arrival times, or the time of the shortest travel path arrival. The method is extended to handle vertically transverse isotropic (VTI) media by an approach which increases the computing time only slightly. The algorithm is tested against synthetic examples for isotropic and VTI wave propagation. 相似文献
15.
使用震源轨迹确定震源位置不仅稳健而且直观,但当介质复杂时震源轨迹难以给出解析解.基于最小走时树射线追踪技术计算震源轨迹的方法(以轨迹所在的残差场中残差最小的点(初始点)至残差较小的点(震源轨迹代表点)的射线路径表示震源轨迹)适用于复杂速度模型,但尚不能正确计算由多段组成的震源轨迹,同时兼顾计算轨迹的完整性和精细性较为困难,计算参数设置烦琐不适于大批量数据的自动处理.针对该方法存在的问题,本文对其进行了改进:(1)采用一种"削皮"算法选取震源轨迹所经过的模型单元的节点作为轨迹代表点;(2)将残差较小的区域作为震源轨迹计算区域(该区域依轨迹分布自适应地划分为若干个连通区域),从未计算的轨迹代表点中选取残差最小者作为射线路径初始点,利用最小走时树算法依次计算所有连通区域内的震源轨迹;(3)通过去掉较短的不再分叉的射线路径使震源轨迹更为精细.虚拟和真实事件的算例表明,改进方法有效克服了原方法的不足,可便捷地计算复杂速度模型中事件的震源轨迹,计算的轨迹精细且较完整. 相似文献
16.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
17.
最短路径射线追踪算法,用预先设置的网络节点的连线表示地震波传播路径,当网络节点稀疏时,获得的射线路径呈之字形,计算的走时比实际走时系统偏大. 本文在波前扩展和反向确定射线路径的过程中,在每个矩形单元内,通过对某边界上的已知走时节点的走时进行线性插值,并利用Fermat原理即时求出从该边界到达其他边界节点的最小走时及其子震源位置和射线路径,发展了相应的动态网络算法. 从而克服了最短路径射线追踪算法的缺陷,大大提高了最小走时和射线路径的计算精度. 相似文献
18.
三维复杂地形近地表速度估算及地震层析静校正 总被引：18，自引：6，他引：18
在地表一致性模型的基础上提出一种可适用于宽线剖面、弯曲测线、传统的二维和目前广泛使用的三维地震观测.在地形及近地表低降速带地质结构复杂的探区，低降速带厚度及速度估算的精度是静校正处理的关键.本研究根据三维地震观测的初至走时数据，利用最小平方与QR分解相结合的算法，在三维空间重建近地表低降速带速度模型，根据重建速度模型实现了静校正长波长分量与短波长分量的同步计算.分析了复杂的近地表低降速带模型初至波的性质，在观测值的自动拾取以及理论值的计算中充分考虑了可能成为初至波的直达波、折射波和反射波的利用，提高了低降速带速度模型反演的精度.在初至走时观测数据的拾取中，本研究采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差，实现了三维初至拾取的大规模全自动化运算.在射线路径与初至波理论走时的计算中，本研究采用一种计算量与模型复杂程度无关的三维射线追踪方法，该方法以最小走时射线路径保证了与观测数据有同等意义的初至波的射线追踪及理论走时的计算.野外实际资料的处理结果表明了方法的有效性. 相似文献
19.
全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性. 相似文献