共查询到19条相似文献，搜索用时 171 毫秒
1.
地震走时层析成像方法是解决复杂近地表模型速度建模问题的重要技术.该方法是一种迭代反演方法,在反演过程中需要反复计算地震射线走时.故而,高效高精度且能适应复杂模型的走时计算方法是地震走时层析成像实用化的关键技术之一.本文引入医学成像领域研究的MSFM(Multi-stencils Fast Marching Methods)用于地震层析反演中的走时计算.该方法在标准FMM(Fast Marching Methods)基础上利用坐标旋转生成新的FMM计算模板,使计算网格点对角方向邻点参与计算,改善了标准FMM存在对角方向误差大的缺陷.本文分析对比了MSFM和标准FMM的计算精度和计算效率;针对地震层析成像技术解决的起伏地表模型建模问题,研究了起伏地表模型地震走时计算的MSFM实现方法;采用炮点邻近区域局部细分网格技术只需增加很少的计算量即可大幅提高计算精度.理论分析和模型试算表明MSFM算法明显改善了FMM的计算精度,同时保持了FMM算法的高效性.文章通过对崎岖地表模型的正演和层析反演试算,验证了基于MSFM的地震走时计算方法对复杂模型有很强的适应能力.研究表明该方法作为地震走时层析反演中高效高精度的正演算法,有很好的应用价值. 相似文献
2.
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力. 相似文献
3.
本文采用Ｈａｍｉｌｔｏｎ理论方法,研究反射地震记录域中地下同一反射点对应的反射走时随炮检距变化的特性,这在叠前地震资料处理中具有重要意义．由于走时函数可以表示成炮点位置、检波点位置、地震波射线出射角和接收角的函数,通过对共反射点走时与炮点位置变化关系的分析,提出了复杂介质中共反射点轨迹可用Ｈａｍｉｌｔｏｎ正则方程描述．在线性变速介质中,走时函数可以解析给出,由此可导出该介质中共反射点轨迹的Ｈａｍｉｌｔｏｎ正则方程．文中结合常速度、横向变速和线性变速介质模型,计算了地下不同反射面上各点相应的共反射点轨迹,阐明Ｈａｍｉｌｔｏｎ方法的一些特点,并首次引入Ｈａｍｉｌｔｏｎ理论方法描述共反射点走时随炮检距的变化特性,开创了该领域研究的一条新途径． 相似文献
4.
在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度. 相似文献
5.
针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性. 相似文献
6.
射线类正演方法以其高效性和灵活性的特点,被广泛应用于地震勘探中.然而,普通射线类正演方法存在焦散区、阴影区和多值走时等问题,计算精度不够理想.为此,本文在传统的各向同性介质声波高斯束正演的基础上,推导了各向异性介质运动学和动力学射线追踪方程,发展了各向异性介质弹性多波射线追踪算法;并将该算法应用到高斯束正演模拟中,实现了一种各向异性介质弹性多波高斯束正演模拟方法.各向异性VTI介质断块和VTI介质复杂构造模型试算的结果表明:本文研究的方法能够对各向异性介质构造进行正演模拟,在保证计算精度的前提下,具有较高的计算效率,模型试算的结果说明了方法的有效性和正确性. 相似文献
7.
8.
《地球物理学进展》2016,(4)
使用Zelt和Barton的方法,通过一个计算效率高的有限差分求解eikonal方程,正演计算走时和射线路径.使用最小二乘QR分解法,求解稀疏线性系统方程组.使用正则化层析反演,结合用户给定的最小的、最平坦和最平滑的扰动限制,每一个加权因子随深度变化.结合数据残差和模型粗糙度的最小化,为数据残差提供一个最平滑的近似模型.该反演方法为非线性反演,需要一个初始模型,在每一次迭代时,需要计算新的射线路径.使用二维初至走时数据,对某油田二维井间地震实际资料进行无限频率初至走时层析反演.将反演所得到的速度与井的测井速度曲线相比较,二者吻合程度较高,表明该反演方法所得速度的分辨率比较高.证实了二维无限频率初至走时层析反演可以为全波形反演提供一个分辨率较高的长波长速度模型,从而为全波形反演井间地震实际资料提供了一个比较可靠的初始速度模型. 相似文献
9.
三维地质模型中地震波共轭梯度非线性走时反演 总被引：2，自引：0，他引：2
地震体波走时层析成像是探测地球内部速度结构的重要方法之一。基于三维块状建模以及三角形拼接的界面描述方式,结合快速高效的逐段迭代射线追踪方法,获得三维复杂地质模型中的地震射线路径与走时信息,采用共轭梯度非线性反演算法,进行地震波走时反演。实验结果表明共轭梯度反演算法在三维层状模型中具有较高的有效性。 相似文献
10.
秦自耕 李子轩 庹先国 李怀良 杜小龙 QINZi-geng LIZi-xuan TUOXian-guo LIHuai-liang DUXiao-long 《地球物理学进展》2015,(2):769-775
地震走时层析建模技术是解决复杂近地表模型速度建模问题的重要技术.该方法是一种迭代反演方法,在地震反演过程中需反复计算地震射线走时.本文采用交互初至拾取方法,通过层析反演出近地表模型,对比分析了走时层析反演表层模型与野外表层模型差别,从而引入了控制点成果约束的层析反演综合建模技术.最终将野外调查成果同走时层析反演建模结合起来,运用非线性走时层析反演技术,通过可视图片、数据、曲线的比较,得出如下结论:在速度纵横向变化大的复杂地区,通过野外资料与室内处理信息交互约束,可以为建立符合实际情况的近地表模型奠定基础,进一步为地震采集资料静校正处理方法的确定和叠加剖面品质提高创造了条件. 相似文献
11.
12.
三维复杂地形近地表速度估算及地震层析静校正 总被引：18，自引：6，他引：18
在地表一致性模型的基础上提出一种可适用于宽线剖面、弯曲测线、传统的二维和目前广泛使用的三维地震观测.在地形及近地表低降速带地质结构复杂的探区，低降速带厚度及速度估算的精度是静校正处理的关键.本研究根据三维地震观测的初至走时数据，利用最小平方与QR分解相结合的算法，在三维空间重建近地表低降速带速度模型，根据重建速度模型实现了静校正长波长分量与短波长分量的同步计算.分析了复杂的近地表低降速带模型初至波的性质，在观测值的自动拾取以及理论值的计算中充分考虑了可能成为初至波的直达波、折射波和反射波的利用，提高了低降速带速度模型反演的精度.在初至走时观测数据的拾取中，本研究采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差，实现了三维初至拾取的大规模全自动化运算.在射线路径与初至波理论走时的计算中，本研究采用一种计算量与模型复杂程度无关的三维射线追踪方法，该方法以最小走时射线路径保证了与观测数据有同等意义的初至波的射线追踪及理论走时的计算.野外实际资料的处理结果表明了方法的有效性. 相似文献
13.
井间地震数据直达波走时层析成像 总被引：1，自引：0，他引：1
本文利用矢量射线追踪正演模拟技术计算地震波直达波传播的路径及走时,进而利用射线走时及路径的内插,发展了弯曲射线迭代反演技术.该方法可用来重建井间地层的速度图像.基于所发展的方法,我们对较为复杂的典型地质模型进行了井间速度重建.结果表明该方法是一项快速、高精度的走时层析成像技术. 相似文献
14.
15.
Field static corrections in general need be applied to all onshore seismic reflection data to eliminate the disturbing effects a weathering layer or near-surface low velocity zone has on the continuity of deep seismic reflections. The traveltimes of waves refracted at the bottom of the low velocity zone (or intermediate refracting interfaces) can often be observed as first breaks on shot records and used to develop a laterally inhomogeneous velocity model for this layer, from which the field static corrections can then be obtained. A simple method is described for computing accurate field statics from first breaks. It is based on a linearization principal for traveltimes and leads to the algorithms that are widely and successfully applied within the framework of seismic tomography. We refine an initial model for the low velocity layer (estimated by a standard traveltime inversion technique) by minimizing the errors between the observed first arrivals on field records and those computed by ray theory through an initial model of the low velocity layer. Thus, one can include more lateral velocity variations within the low velocity layers, which are important to obtain good field static corrections. Traditional first break traveltime inversion methods cannot, in general, provide such refined velocity values. The technique is successfully applied to seismic data from the Amazon Basin. It is based on a simple model for the low velocity layer that consists of an undulating earth surface and one planar horizontal refractor overlain by a laterally changing velocity field. 相似文献
16.
Fei Li Tao Xu Minghui Zhang Zhenbo Wu Chenglong Wu Zhongjie Zhang Jiwen Teng 《地震科学（英文版）》2014,27(2):127-136
Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model. 相似文献
17.
描述地震波衰减特征的品质因子Q对地震数据处理和油藏描述非常重要,在地震勘探领域,Q值一般通过垂直地震剖面(VSP)数据或地面地震数据得到.由于叠前地面地震数据具有复杂的射线路径且存在噪声、调谐干涉效应等影响,从叠前地震数据中准确估计Q值相对困难.本文以地震波射线传播为基础,根据同相轴局部斜率和射线参数的映射关系,将多射线波形频谱同时带入谱比法联合反演估计Q值,提出了基于多射线联合反演的速度无关叠前Q值估计方法.该方法通过局部斜率属性避开了速度对Q值估计的影响,局部斜率携带地震波传播的速度信息,具有相同局部斜率的地震反射波具有相同的传播射线参数.同相轴局部斜率是地震数据域的属性,而速度是模型域的参数,在估计Q值中采用数据域的属性参数可以直接应用于数据的联合反演,而不需要通过速度对其做进一步的转化,从而提高了Q值估计的精度.同时,本方法采用预测映射(predictive mapping)技术将非零炮检距反射信息映射到零炮检距处,从而获得零偏移距走时对应的Q值.模拟和实际算例验证了本文方法的有效性. 相似文献
18.
19.
全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性. 相似文献