基于地震图像校准的共成像点道集增强技术

由共成像点道集抽取的共偏移距道集可以当成相同地下成像的多次观测.由于偏移误差的影响,在不同共偏移距道集上,同一采样点存在水平和垂直方向上的错位.本文提出一种基于地震图像校准的共成像点道集增强技术,实现了不同偏移距道集在时间和空间上的逐点匹配对齐.在本文中以2D局部归一化互相关来表征共偏移距道集和叠加道集在相同时空位置的相似度,假设在不同水平和垂直移动量时的互相关满足连续凸函数,利用求导方法估计共偏移距道集在该位置处的非整数校正量,最后采用双线性内插方法估计成像振幅.传统道集拉平技术在垂直方向进行校正量消除,本文方法能有效估计共偏移距道集中的水平和垂直校正量,并在亚像素域估计正确的成像振幅.模型数据和实际数据的处理结果表明,本文方法能有效增强共成像点道集中同相轴的一致性,提高叠加结果的分辨率.     

基于波动方程的广义屏叠前深度偏移

地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具.     

一种倾角时差校正和叠前成象方法

本文从测量射线参数出发进行反向射线追踪,导出倾角时差校正（DMO）的公式。经过DMO后,可以从一组等炮检距剖面得出共分角线点道集。用于对这些道集进行叠加的速度值与界面倾角无关。对经过DMO的资料的等时切片进行叠前成象（PSI）,就可以把分布在圆上的绕射能量沿圆弧加起来,并放在圆弧上对应于最大炮检距的位置。经过这两种处理,再应用标准的速度分析和叠加方法,就可得出偏移后的剖面。这两种处理均与速度无关。最后用物理模型试验说明了DMO和PSI的效果是好的。     

Kirchhoff型偏移反偏移串联地震数据炮检距转换技术

零炮检距数据在海洋地震资料处理中有很多优势和用途,然而,受海洋水平拖缆地震采集作业方式的限制,接收数据的最小炮检距一般在200 m左右,小于该炮检距的数据是无法直接获得的.通常的做法是通过对较大炮检距数据进行动校,通过外推来变相获得零炮检距（包括小炮检距）数据,其外推的精度会受到地震资料信噪比、动校正速度的精度等因素影响,并且保幅性较差.本文通过一种基于Kirchhoff真振幅偏移和反偏移串联的技术,在反偏移过程中改变观测系统,有效实现了大炮检距反射地震数据向零炮检距（包括小炮检距）数据之间的转换,且很好地保持了零炮检距（包括小炮检距）数据的振幅特性.同时,经偏移-反偏移串联处理后,有效压制了地震数据中的随机噪声,地震资料信噪比和成像精度均得到显著提高.     

台湾西南海域地震数据处理及天然气水合物地震属性

对台湾西南海域增生楔部位长排列多道地震数据进行地震成像、速度分析、AVO分析、AVO反演处理,获得了天然气水合物多属性地震特征.在偏移剖面上,BSR与海底近似平行,极性与海底相反,穿越沉积层.AVO分析显示,强BSR振幅部位,BSR振幅随偏移距增大而增大.精细速度分析表明强BSR振幅下方存在纵波低速层.对应于强BSR振幅部位,AVO反演的P波、G波为相对高负值区,位于P、G交会图的第三象限,该部位泊松比变化率为负值,横波反射系数接近于零.以上多属性地震特征均预示着该区域可能存在天然气水合物层,且天然气水合物层下方可能存在游离甲烷气层.     

振幅保真的单程波方程偏移理论

本报告综述了近年来发展起来的真振幅单程波方程偏移理论,对各种基于单程波方程的偏移振幅作出了系统的分析介绍.为了得到正确的偏移振幅,叠后偏移之前必须进行球面扩散校正.但是,作为叠后相位移偏移算法的推广,Dubrulle氏的共炮检距偏移除了零炮检距和平层两种特例外,它不能给出正确的偏移振幅.通过分析,我们发现传统的单程波单炮偏移不是保振幅算法.利用真振幅单程波方程分解和校正地表的初值条件,我们可以将常规偏移方法改造为真振幅算法并且证明在高频渐近的意义下,新方法的偏移振幅等价于Kirchhoff反演的结果.进一步研究发现可以利用单平方根算子和双平方根算子输出真振幅共反射角剖面.我们的分析指出,这时正确的成像条件为乘积pUp*D.与真振幅共炮偏移所需的商型成像条件pU/pD相比,共反射角偏移的计算稳定性大为改善并且在实际地震资料处理中有重要的应用.     

利用改进的共偏移距—共反射面(COCRS)叠加压制地滚波（英文）

我们改进了共偏移距-共反射面(COCRS)法,可用以衰减地滚波,即由于低速、低频/高振幅瑞雷波通常产生的相干噪声。COCRS算子是基于双曲线,因此它可以拟合双曲走时的同相轴,如叠前数据中的反射同相轴。相反,地滚波在共中点(CMP)和共炮点道集中是线性的并可以可以利用COCRS算子鉴别与压制。因此,我们在共偏移距剖面之前共炮道集中搜索反射倾斜和曲率。因为这对反射振幅的危害最小化是最理想的,我们只对在地滚波区多次覆盖的数据进行叠加。在CO剖面前搜索CS道聚集是对常规COCRS叠加的另一个改进。我们使用合成和真实数据集测试了所提出的方法,数据采自伊朗西部地区。我们将本方法压制地滚波的结果与f-k滤波和f-k滤波后常规COCRS叠加压制地滚波的结果进行了比较。结果表明,该方法对真伪滚压制效果优于F-K滤波与传统CRS叠加。然而,计算时间高于其他常规的方法,如f滤波。     

地震勘探中广义弹性阻抗的正反演

常规的地震道反演方法建立在反射P波垂直入射假设 的基础上,而实际地震资料采集时多数是非零炮检距的，反射振幅是共中心点道集叠加的结 果 . 因此，利用常规地震道反演方法就不能得到可靠的波阻抗或其他岩性信息. 本文利用Patr ick Connolly弹性阻抗的思想，通过对Zoeppritz方程的进一步简化，推导出适合常规叠后 资料的、非零炮检距条件下纵波反射系数递推公式，提出了广义弹性阻抗的概念，解决了非 零炮检距条件下，常规叠后地震道正反演的关键问题. 广义弹性阻抗不仅包含波阻抗，还包 含了纵横波速度等岩性信息，具有很好的实用价值. 进行广义弹性阻抗的反演，能较常规地 震道反演获得更多、更可靠的流体、孔隙度、砂泥含量等信息，有助于解释常规地震道反演 和道积分剖面中的假象，降低反演的多解性，提高储层预测的精度.     

地震勘探空间分辨力分析

未偏移地震剖面的空间分辨力研究通常采用第一菲涅尔带的解释方法,它是以几何地震学的射线理论为基础,仅考虑了地震波的运动学特征,不能详细描述地震波的动力学特征.本文从物理地震学的广义绕射叠加观点出发,综合分析地震波的运动学和动力学特征,由绕射叠加效应推导出了零炮检距地震道空间分辨力的计算公式.该公式的计算结果与第一菲涅尔带半径在数值上相当,并通过正演模拟验证了计算公式的合理性.因此利用物理地震学的原理研究地震勘探空间分辨力更符合地震波客观实际的传播特征.     

基于多次波和自适应预测误差滤波器的地震近炮检距数据重建方法

地震数据采集过程中,受野外施工条件的制约,往往很难获得完整的地震波场,尤其近炮检距数据的缺失尤为严重.当前,很多地震数据处理方法的应用都依赖于近炮检距数据,如何对近炮检距的缺失数据进行重建,是一个重要的研究课题.本文通过多次波和一次波的互相关构建虚拟一次波,利用数据本身的波场信息,对缺失的近炮检距数据进行插值重建.由于通过多次波构建的一次波与真实的一次波存在振幅和相位方面的差异,提出通过基于非平稳自回归过程的自适应预测误差滤波器来表征虚拟一次波的能量谱,利用最小二乘反演方法重建近炮检距缺失数据,自适应预测误差滤波器通过求解正则化约束下的数学欠定问题来实现局部自适应特征.通过对Sigsbee2B模型和实际数据的测试结果表明新方法可以合理地重建复杂的近炮检距缺失数据.     

Pre-stack Kirchhoff depth migration local space-shift imaging condition: synthetic and data examples

Extracting accurate common image angle gathers from pre-stack depth migrations is important in the generation of any incremental uplift to the amplitude versus angle attributes and seismic inversions that can lead to significant impacts in exploration and development success. The commonly used Kirchhoff migration outputs surface common offset image gathers that require a transformation to angle gathers for amplitude versus angle analysis. The accuracy of this transformation is one of the factors that determine the robustness of the amplitude versus angle measurements. Here, we investigate the possibility of implementing an extended imaging condition, focusing on the space-lag condition, for generating subsurface reflection angle gathers within a Kirchhoff migration. The objective is to determine if exploiting the spatial local shift imaging condition can provide any increase in angle gather fidelity relative to the common offset image gathers. The same restrictions with a ray-based approach will apply using the extended imaging condition as both the offset and extended imaging condition method use travel times derived from solutions to an Eikonal equation. The aims are to offer an alternative ray-based method to generate subsurface angle gathers and to understand the impact on the amplitude versus angle response. To this end, the implementation of the space-shift imaging condition is discussed and results of three different data sets are presented. A layered three-dimensional model and a complex two-dimensional model are used to assess the space shift image gathers output from such a migration scheme and to evaluate the seismic attributes relative to the traditional surface offset common image gathers. The synthetic results show that the extended imaging condition clearly provides an uplift in the measured amplitude versus angle over the surface offset migration. The noise profile post-migration is also improved for the space-lag migration due to the double summation inside the migration. Finally, we show an example of a space-lag gather from deep marine data and compare the resultant angle gathers with those generated from an offset migration and a time-shift imaging condition Kirchhoff migration. The comparison of the real data with a well log shows that the space-lag result is a better match to the well compared to the time-lag extended imaging condition and the common offset Kirchhoff migration. Overall, the results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers with an incremental improvement in amplitude versus angle fidelity and lower noise but comes at a higher computational cost.     

Investigating a time‐shift extended imaging condition in a Kirchhoff pre‐stack depth migration algorithm

Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study.     

叠前地震数据的平面波深度偏移法

提出了一套基于平面波分解的波动方程叠前地震数据深度偏移方法. 通过对共炮点道集和共偏移距道集地震数据的平面波分解，分别得到适用于单平方根波场外推方程和双平方根波场外推方程的共ps（炮点坐标平面波参数）平面波道集和共ph（偏移距坐标平面波参数）平面波道集. 在对共炮点道集和共偏移距道集地震数据的平面波分解时，不需要进行通常意义下的τ p变换计算. 通过对共ps平面波道集和共ph平面波道集的偏移效果对比，我们认为在速度弱横向变化介质中，两种平面波道集偏移方法的效果相当，但对于速度强横向变化介质，共ps平面波道集偏移方法的效果要优于共ph平面波道集偏移方法. 在计算效率方面，共ps平面波道集偏移方法与共ph平面波道集偏移方法基本相同.     

Enabling affordable omnidirectional subsurface extended image volumes via probing

Image gathers as a function of subsurface offset are an important tool for the inference of rock properties and velocity analysis in areas of complex geology. Traditionally, these gathers are thought of as multidimensional correlations of the source and receiver wavefields. The bottleneck in computing these gathers lies in the fact that one needs to store, compute, and correlate these wavefields for all shots in order to obtain the desired image gathers. Therefore, the image gathers are typically only computed for a limited number of subsurface points and for a limited range of subsurface offsets, which may cause problems in complex geological areas with large geologic dips. We overcome increasing computational and storage costs of extended image volumes by introducing a formulation that avoids explicit storage and removes the customary and expensive loop over shots found in conventional extended imaging. As a result, we end up with a matrix–vector formulation from which different image gathers can be formed and with which amplitude‐versus‐angle and wave‐equation migration velocity analysis can be performed without requiring prior information on the geologic dips. Aside from demonstrating the formation of two‐way extended image gathers for different purposes and at greatly reduced costs, we also present a new approach to conduct automatic wave‐equation‐based migration‐velocity analysis. Instead of focusing in particular offset directions and preselected subsets of subsurface points, our method focuses every subsurface point for all subsurface offset directions using a randomized probing technique. As a consequence, we obtain good velocity models at low cost for complex models without the need to provide information on the geologic dips.     

用于AVO分析的振幅保真平面波叠前时间偏移

To support amplitude variation with offset （AVO） analysis in complex structure areas, we introduce an amplitude-preserving plane-wave prestack time migration approach based on the double-square-root wave equation in media with little lateral velocity variation. In its implementation, a data mapping algorithm is used to obtain offset-plane-wave data sets from the common-midpoint gathers followed by a non-recursive phase-shift solution with amplitude correction to generate common-image gathers in offset-ray-parameter domain and a structural image. Theoretical model tests and a real data example show that our prestack time migration approach is helpful for AVO analysis in complex geological environments.     

Eliminating 3D pre‐stack migration artefacts by 6D filtering

State‐of‐the‐art 3D seismic acquisition geometries have poor sampling along at least one dimension. This results in coherent migration noise that always contaminates pre‐stack migrated data, including high‐fold surveys, if prior‐to‐migration interpolation was not applied. We present a method for effective noise suppression in migrated gathers, competing with data interpolation before pre‐stack migration. The proposed technique is based on a dip decomposition of common‐offset volumes and a semblance‐type measure computation via offset for all constant‐dip gathers. Thus the processing engages six dimensions: offset, inline, crossline, depth, inline dip, and crossline dip. To reduce computational costs, we apply a two‐pass (4D in each pass) noise suppression: inline processing and then crossline processing (or vice versa). Synthetic and real‐data examples verify that the technique preserves signal amplitudes, including amplitude‐versus‐offset dependence, and that faults are not smeared.     

三维偏移距平面波有限差分叠前时间偏移

本文提出了中点-半偏移距域内的三维偏移距平面波(offset plane-wave)方程,并给出了其有限差分解法.偏移距平面波可通过对CMP道集进行平面波分解(倾斜叠加或线性Radon变换)生成,然而这样做会产生严重的噪音干扰.本文提出了局部倾斜叠加方法(local slant-stacking)来消除离散线性Rado...     

Kirchhoff叠前时间偏移角度道集

提出三维Kirchhoff叠前时间偏移角度域共像点道集的改进算法,克服传统角度求取算法局限,可计算相对倾斜地层法线入射角;与Kirchhoff直射线叠前时间偏移求角度算法相比,本文方法考虑射线弯曲效应,包含层速度,角度范围加大,更接近真实入射角;计算走时采取弯曲射线或者适应线性横向变速介质的非对称走时等算法,角度道集在大角度处得到拉平;采用相对保幅的权因子以及覆盖次数校正技术,有利于叠前AVA反演.模型测试结果表明：叠前时间偏移角度道集,相对CMP、CRP所转化角度道集,更准确反应AVA效应;实际三维数据测试表明本文方法可以提供品质优良的角度道集,适用于AVA分析、反演,提高叠前反演分辨率.