基于快速推进迎风双线性插值法的三维地震波走时计算

三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先, 在常规双线性插值公式推导过程中, 充分利用平面波双线性假设的结论, 获得了二元极小值超越方程的解析解, 进而推导出了准确的局部走时计算公式, 同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次, 引入迎风差分的基本思想, 提出迎风双线性插值的局部走时计算策略, 该计算策略能简化算法、提高效率且保证无条件稳定性;然后, 将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合, 获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后, 通过精度和效率分析检验了新算法的精度、效率和正确性, 并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性.     

三维复杂介质中转换波走时快速计算

复杂介质中转换波走时计算是多波勘探地震学中重要内容之一.本项研究利用惠更斯原理和费玛原理,获得了三维复杂介质中转换波快速计算的改进型最小走时树方法.其中,在保证精度的条件下,为了提高三维转换波走时计算效率,首先对初至波最小走时树基本算法进行了改进.本方法通过将转换波分为上、下行波分别进行射线追踪以实现三维转换波走时的快速计算.模型计算表明,方法的计算速度快,而且稳定性强,对多波地震勘探具有较大的应用价值.     

三维起伏地表条件下地震波走时计算的不等距迎风差分法

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型.     

三维起伏地表条件下地震波走时计算的不等距迎风差分法

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型.     

三维复杂黄土塬区变加密不等距网格初至走时计算与特征分析

三维复杂黄土塬区初至走时计算与特征分析,对于分析复杂黄土塬区地震波运动学特征、观测系统设计、初至拾取及层析反演等均具有重要意义.针对该问题：首先,采用三维分区局部变加密不等距网格剖分复杂黄土塬速度模型;然后,建立基于费马原理和迎风格式的双线性插值局部走时算法,采用修正后的群推进法作为波前扩展方式;最后,以理论模型考察了算法的精度与效率、分析了三维复杂黄土塬区地震初至波的走时特征.结果表明：变加密网格策略能在增加很少计算量的情况下,大幅降低因震源奇异性引起的计算误差问题,同时还能无条件稳定且灵活的适应三维复杂黄土塬区的复杂地震地质条件;和层状均匀模型相比,黄土塬区地震初至波的走时场和时距曲线较为复杂.

     

3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching

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.     

A case history of secondary 3D seismic exploration in the Liaohe basin: An example from the Rongshengpu-Qianjin area

The working area is located in the industrially developed region of Rongshengpu-Qianjin, where a surface water system is developed, surface-layer lithology is complicated, and various kinds of hydrocarbon traps are buried at depth. The seismic data acquired previously couldn’t be interpreted due to the complex surface and geological conditions. Taking secondary 3D seismic from the Rongshengpu-Qianjin area as an example, this paper describes a set of techniques designed to overcome these difficulties and improve the quality of seismic data. The applied techniques included flexible acquisition geometry, low-noise receiver conditions, quantitative quality control, and so on. First Author Tian Qingbo, Engineer, graduated in Applied Geophysics Department from Chengdu Geology College in 1991. He engaged in acquisition and processing of seismic data, research of acquisition method for a long time. Now, he works in the department of Project and Technique in Geophysical Exploration Company of Liaohe Petroleum Exploration Bureau.     

芦山震区地壳三维P波速度精细结构及地震重定位研究

联合芦山地震序列5285个地震的50711条P波初至绝对到时数据及7294691条高质量的相对到时数据,利用双差地震层析成像方法联合反演了芦山震源区高分辨率的三维P波速度精细结构及5115个地震震源参数.反演结果表明,芦山主震震中为30.28°N,103.98°E,震源深度为16.38 km,主震南西段余震扩展长度约23 km,余震前缘倾角较和缓,主震北东段余震扩展长度约12 km,余震前缘呈铲形,倾角较陡.芦山震源区P波三维速度结构表现出明显的横向不均匀性,近地表处的P波速度异常与地形起伏及地质构造密切相关:宝兴杂岩对应明显的高速异常,此异常由地表延伸到地下15 km深度附近,而中新生代岩石表现为低速异常;大兴附近区域亦显示出小范围的大幅度高速异常,宝兴高速异常与大兴高速异常在10 km深度附近相连,进而增加了芦山震源区的高低速异常对比幅度.在芦山主震的南西、北东两段速度结构存在着较大差异,芦山主震在水平向位于宝兴及大兴高速异常所包围的低速异常的前缘.主震南西段余震主要发生在倾向北西的高低速异常转换带上并靠近低速一侧,其下盘为低速异常,上盘为高速异常.而芦山主震北东段的余震主要分布在宝兴高速体与大兴高速体之间,主发震层向北西倾斜,主发震层上方的宝兴高速异常下边界出现一条南东倾向的反冲地震带,两地震带呈"y"型分布.     

基于印模法的地面磁电阻率数据三维非线性共轭梯度反演

本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识：（1）由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;（2）采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;（3）在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;（4）确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中.