模型约束三维折射静校正方法研究

本文详细的介绍了模型约束三维折射静校正技术的基本原理和计算步骤,具体给出了三维抽线,目的层质量监控,三维速度内插,三维底界内插等技术的实现方法,并结合了西部某工区的具体实例加以说明.结果表明,这种由微测井、小折射求取近地表深度-时间速度模型,由大炮记录折射初至旅行时求取各点延迟时,两者联合反演低降速层底界的厚度和高程,最终求取低降速层静校正量的方法,可明显提高低降速带的反演精度,较好地解决表层结构复杂地区的静校正问题.     

川东北复杂山地三维静校正应用及实例分析

川东北复杂山地由于地表高程变化大、风化层厚度不均匀、低降速带速度横向变化大、高速层底界不稳定等因素,得到的地震资料具有严重的静校正问题,因此,如何解决山地静校正问题是后续地震资料处理的关键之一,本文首先分析高程静校正、折射静校正、层析静校正的基本原理和适用条件,结合实际复杂山地三维地震资料的特点,进行试验对比,提出了进行山地静校正的基本思路,即：首先进行高程静校正,这样可以利用高程静校正更容易拾取初至时间,然后利用折射静校正结合微测井等资料建立近地表速度-深度模型,以此速度-深度模型作为层析静校正的初始模型进行迭代处理,最后得到最终的近地表速度-深度模型和静校正值.根据以上处理流程,我们建立了适合于川东北山地三维复杂地表地震资料的静校处理正方法,并在实际生产过程中取得了良好的效果.     

三维复杂地形近地表速度估算及地震层析静校正

在地表一致性模型的基础上提出一种可适用于宽线剖面、弯曲测线、传统的二维和目前广泛使用的三维地震观测.在地形及近地表低降速带地质结构复杂的探区，低降速带厚度及速度估算的精度是静校正处理的关键.本研究根据三维地震观测的初至走时数据，利用最小平方与QR分解相结合的算法，在三维空间重建近地表低降速带速度模型，根据重建速度模型实现了静校正长波长分量与短波长分量的同步计算.分析了复杂的近地表低降速带模型初至波的性质，在观测值的自动拾取以及理论值的计算中充分考虑了可能成为初至波的直达波、折射波和反射波的利用，提高了低降速带速度模型反演的精度.在初至走时观测数据的拾取中，本研究采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差，实现了三维初至拾取的大规模全自动化运算.在射线路径与初至波理论走时的计算中，本研究采用一种计算量与模型复杂程度无关的三维射线追踪方法，该方法以最小走时射线路径保证了与观测数据有同等意义的初至波的射线追踪及理论走时的计算.野外实际资料的处理结果表明了方法的有效性.     

黄土塬区模型约束折射静校正技术应用研究

鄂尔多斯盆地东缘黄土塬区地形起伏剧烈,低降速层厚度和速度横向变化大,表层结构复杂,折射波静校正表层模型建模困难,静校正问题突出.本文通过黄土塬区表层结构规律分析,提出了一种基于初至走时层析求取近地表速度、以等速度界面近似高速顶界面的表层模型约束建模方法;在基准面静校正概念的基础上,对比分析了基岩顶界面、层析速度收敛面等不同折射面解释方案对静校正量、标志层反射时间及平均速度的影响,结合实际钻井资料证实了在黄土塬区以基岩顶界面作为折射面能够获得正确的静校正量、标志层反射时间与平均速度.     

起伏地表静校正浮动基准面选取方法研究

本文分析了静校正中不同浮动基准面的特点及确定方法。通过理论模型试验对平滑地表、平均静校正量与最小静校正误差等浮动基准面的静校正效果进行比较,验证了在最小静校正误差基准面上得到的叠加速度仅取决于低速带底界下伏地层的速度,而与地形起伏、低速带结构无关,得到的叠加剖面具有较好的同相叠加效果。另外为了更加符合实际资料处理情况,本文中采用波动方程模拟数据进行理论模型试验,并由此给出了一种新的基于波形的目标函数计算方法。同时修改了最小静校正误差计算公式,使其适用于起伏底界非均匀速度模型的实际资料静校正处理。最后部分对实际资料的处理进一步证明了该方法的优越性。     

厚风化层覆盖区转换波静校正方法

P-SV转换波处理与传统的P-P波处理有很大的不同,如S波静校正、CCP叠加、P-SV速度分析和偏移等,其中最大的难题就是S波静校正问题.S波速度基本不受潜水面的影响,与纵波静校正没有直接相关性,有时横波静校正量能达到纵波静校正量的十倍,用纵波静校正量乘以比例系数来解决横波静校正问题将导致较大误差.同一接收点X和Z分量存在一定的初至时差,该时差代表了P波和S波在低降速带的走时差,可以利用该时差和近地表纵横波速度比信息去除低降速带对横波的影响,得到准确的静校正量.本文利用多分量初至时差推导了较为精确的横波静校正公式,再结合共检波点叠加求取剩余静校正量的方法,形成了完整的转换波静校正配套方法.利用该方法对苏里格气田二维及三维多波地震资料进行了实际处理,数据处理结果证明了该方法的有效性,该方法尤其适用于其他方法难以奏效的风化层较厚地区的横波静校正量求解,该方法也同时考虑了长波长横波静校正问题.     

复杂山区初至波层析反演静校正

提高静校正精度是取得复杂山区良好地震成像的一个重要条件.而建立在水平折射面假设基础之上折射波静校正方法,无论是假设前提还是实际应用效果,都不适应于地表剧烈起伏,速度纵、横向变化大的复杂区.为此本文提出使用初至波层析反演静校正方法,即利用地震记录中初至旅行时反演出表层速度模型,计算出炮点和检波点的静校正量.通过正演模拟数据和实际资料的验证,很好的解决了复杂地表引起的静校正问题.     

川东南复杂山地灰岩区静校正技术综合应用研究

复杂山地灰岩出露区勘探存在严重的静校正问题,解决静校正问题难度较大,由于静校正问题的存在会降低地震资料的垂向分辨率、出现假构造、影响速度分析,严重影响地震资料成像质量,降低勘探成功率.为解决复杂山地灰岩区勘探存在的静校正问题,本文研究了高程静校正、折射静校正和初至层析静校正,分析对比不同方法的优缺点和应用效果,在此基础上提出并研究融合静校正技术,在本研究区中,该方法是采用折射静校正和初至层析静校正量的优势区域进行融合,充分发挥不同方法的优点和适应性,较好的解决一次静校正;针对剩余静校正问题,提出了基于反射波剩余静校正、成像域射线束剩余静校正和速度建模的循环迭代综合剩余静校正处理技术,分步逐级迭代解决剩余静校正问题,该方法集成了反射波剩余静校正、成像域射线束剩余静校正的优点,明显提高地震剖面的细节成像质量.通过本文研究,建立了一套针对复杂山地灰岩区低信噪比资料的静校正处理技术序列并形成一套处理流程,在川东南复杂山地灰岩区低信噪比数据处理中,验证了该技术序列的可行性,取得比较好的效果.     

三维地震与地面微地震联合校正方法

由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景.     

结合基准面重建的叠前时间偏移方法

提出了一种结合虚拟界面、瑞利积分和相移法的混合的基准面重建方法.通过与叠前时间偏移方法结合,形成了针对起伏地表采集数据的叠前时间偏移方法和新流程.该方法能正确考虑波在近地表传播的实际路径,克服了高速层出露时静校正方法的误差;它也能自己确定虚拟层速度,避免了现行基于波场延拓的基准面重建方法需要准确近地表速度的困难.文中分别用近地表存在明显低速层和近地表有高速层出露这两类模型的理论数据,验证了所发展方法和流程的有效性.     

基于共姿态道集的静校正方法

静校正问题是地震勘探的关键问题,直接影响地震勘探精度和准确性.实际地震采集过程中,当在相同接收点位置上不同时间内插拔布设了不同的检波器时,对于目前基于地表一致性理论假设的基准面静校正和剩余静校正,以及非地表一致性剩余静校正都不具备适用条件.为解决这一问题,本文提出了基于共姿态道集的静校正方法,将相同接收点位置上不同时间布设的检波点所接收的地震数据抽成不同的共姿态道集,在共姿态道集内实施地表一致性静校正;当某接收点位置上具有若干个共姿态道集时,该接收点位置上可能会存在多个检波点静校量;炮点静校问题仍然采用地表一致性静校正方法解决.该方法解决了同一接收点位置上不同共姿态道集之间的非地表一致性静校正问题,同时也解决了全区的检波点和炮点的地表一致性静校正问题,在实际数据应用效果明显.     

CRS-stack-based seismic imaging considering top-surface topography

In this case study we consider the seismic processing of a challenging land data set from the Arabian Peninsula. It suffers from rough top‐surface topography, a strongly varying weathering layer, and complex near‐surface geology. We aim at establishing a new seismic imaging workflow, well‐suited to these specific problems of land data processing. This workflow is based on the common‐reflection‐surface stack for topography, a generalized high‐density velocity analysis and stacking process. It is applied in a non‐interactive manner and provides an entire set of physically interpretable stacking parameters that include and complement the conventional stacking velocity. The implementation introduced combines two different approaches to topography handling to minimize the computational effort: after initial values of the stacking parameters are determined for a smoothly curved floating datum using conventional elevation statics, the final stack and also the related residual static correction are applied to the original prestack data, considering the true source and receiver elevations without the assumption of nearly vertical rays. Finally, we extrapolate all results to a chosen planar reference level using the stacking parameters. This redatuming procedure removes the influence of the rough measurement surface and provides standardized input for interpretation, tomographic velocity model determination, and post‐stack depth migration. The methodology of the residual static correction employed and the details of its application to this data example are discussed in a separate paper in this issue. In view of the complex near‐surface conditions, the imaging workflow that is conducted, i.e. stack – residual static correction – redatuming – tomographic inversion – prestack and post‐stack depth migration, leads to a significant improvement in resolution, signal‐to‐noise ratio and reflector continuity.     

CRS-stack-based residual static correction

In the case of onshore data sets, the acquired reflection events can be strongly impaired due to rough top‐surface topography and inhomogeneities in the uppermost low‐velocity layer, the so‐called weathering layer. Without accounting for these influences, the poor data quality will make data processing very difficult. Usually, the correction for the top‐surface topography is not perfect. The residuals from this correction and the influence of the weathering layers lead to small distortions along the reflection events. We integrated a residual static correction method into our data‐driven common‐reflection‐surface‐stack‐based imaging workflow to further eliminate such distortions. The moveout‐corrected traces and the stacked pilot trace are cross‐correlated to determine a final estimate of the surface‐consistent residual statics in an iterative manner. As the handling of top‐surface topography within the common‐reflection‐surface stack is discussed in a separate paper in this special issue, the corresponding residual static correction will be explained in more detail. For this purpose, the results obtained with a data set from the Arabian Peninsula will be presented.     

THE RESOLUTION OF NARROW LOW-VELOCITY ZONES WITH THE GENERALIZED RECIPROCAL METHOD1

The depth to the surface of a refractor and the seismic velocity within the refractor are very often intimately related. In the shallow environment, increased thicknesses of weathering occur in areas of jointing, shearing or lithological variations, and these zones of deeper weathering can have lower subweathering refractor velocities. This association is important in geotechnical investigations and in the measurement of weathering thicknesses and sub-weathering velocities for statics corrections for reflection seismic surveys. Algorithms, which employ forward and reverse traveltime data and which explicitly accommodate the offset distance through the process known as refraction migration, are necessary if detailed structure on a refractor and rapid lateral variations of the seismic velocity within it are to be resolved. These requirements are satisfied with wavefront construction techniques, Hales’ method and the generalized reciprocal method (GRM). However, these methods employ refraction migration in fundamentally different manners. Most methods compute an offset distance with an often imprecise knowledge of the seismic velocities of the overlying layers. In contrast, the GRM uses a range of offset distances from less than to greater than the optimum value, with the optimum value being selected with a minimum-variance criterion. The approach of the GRM is essential where there are undetected layers and where there are rapid variations in the depth to a refractor and the seismic velocity within it. In the latter situations the offset distance necessary to define the seismic velocities can differ considerably from the value required to define depths. The efficacy of the GRM in resolving structure and seismic velocity is demonstrated with three model studies and two field examples.     

基于偏移成像道集的剩余静校正方法

针对陆上地震资料处理的静校正问题,提出了一种基于偏移成像道集的剩余静校正方法.与传统的由动校正后的CMP道集中拾取剩余时差不同,本文基于偏移成像道集求取剩余时差,避免了复杂情况下同相轴归位不准确导致的剩余时差拾取误差.通过生成随炮点和检波点位置变化的偏移道集,实现了由偏移道集中直接拾取炮、检点的地表一致性剩余时差;该炮、检点偏移道集只在指定的局部时窗生成,并不增加大的计算量.二维和三维实际数据测试表明了该方法的有效性和实用性.     

低信噪比转换波地震资料静校正(英文)

Converted waves have slow velocity and low signal-to-noise ratio. It is also difficult to pick first-breaks and bin the common-conversion-points （CCP）. Some statics methods, which work well for P-wave data, can＇t be effectively used for solving convertedwave statics problems. This has become the main obstacle to breakthroughs in convertedwave data processing. To improve converted-wave static corrections, first, a statics method based on the common-receiver-point （CRP） stack is used for the initial receiver static correction to enhance the coherency of the CRP stack. Second, a stack-power-maximization static correction which improves the continuity of the CCP stack is used for detailed receiver statics. Finally, a non-surface-consistent residual moveout correction of the CCP gathers is used to enhance the stack power of reflection signals from different depths. Converted-wave statics are solved by the joint use of the three correction methods.     

Directional statics in common reflection surface stack for rugged surface topography

Seismic data acquired along rugged topographic surfaces present well‐known problems in seismic imaging. In conventional seismic data processing, datum statics are approximated by the surface consistence assumption, which states that all seismic rays travel vertically in the top layer. Hence, the datum static for each single trace is constant. In case this assumption does not apply, non‐constant statics are required. The common reflection surface (CRS) stack for rugged surface topography provides the capability to deal with this non‐vertical static issue. It handles the surface elevation as a coordinate component and treats the elevation variation in the sense of directional datuming. In this paper I apply the CRS stack method to a synthetic data set that simulates the acquisition along an irregular surface topography. After the CRS stack, by means of the wavefield attributes, a simple algorithm for redatuming the CRS stack section to an arbitrarily chosen planar surface is performed. The redatumed section simulates a stack section whose acquisition surface is the chosen planar surface.     

黄土塬覆盖区的层析反演静校正方法研究

我国北方地区黄土塬覆盖区的静校正问题是地震数据处理中的难点问题之一.黄土塬表层覆盖巨厚黄土,高差起伏较大,地震静校正问题严重;而且黄土塬覆盖区的潜水面普遍较深,常规折射静校正方法无法取得令人满意的处理效果.本文针对黄土塬覆盖区的静校正难题,研究层析反演静校正方法在黄土塬地区的适用性和可靠性.层析反演静校正利用地震波初至走时数据、通过迭代反演的方法构建速度模型,进而依据所得的近地表速度模型对地震数据进行静校正处理.本文的迭代反演采用同步迭代重构算法（SIRT）,并且对同步迭代重构算法进行了改进,使得层析反演的迭代过程趋于稳定.但是,因为黄土塬覆盖区地表高程的横向变化剧烈,相邻检波点的高差及其静校正量有时差异很大,在运用层析静校正求取长波长静校正量的同时,还需采用初至波剩余静校正方法求取短波长静校正量.实例证明,综合应用依据初至波走时数据的层析静校正和剩余静校正方法,同时计算长波长和短波长的静校正量,能够有效地解决黄土塬覆盖区实际地震资料的静校正问题.     

三维波动方程基准面校正方法的应用研究

波动方程基准面校正处理被认为是当地表高程变化剧烈、地表一致性假设又不成立的情况下对于常规高程基准面校正的必要替代. 在二维波动方程基准面校正方面,前人已经作了大量工作,并且在实际应用中获得了良好的效果,应该将其进一步推广到三维. 本文采用三维频率空间域有限差分波场延拓算子以“逐步－累加”的方式实现了三维波动方程基准面校正,并对实际数据进行了处理. 对于西部某三维数据的实际计算结果表明,相对于传统的模型法高程静校正,波动方程基准面校正更合理地实现了基准面校正,有助于提高后续的速度分析精度和成像质量.