线性拉东域预测反褶积在海洋多次波去除中的应用

海洋多次波是海洋地震资料处理中最难去除的噪声,目前时空域预测反褶积、SRME、高精度拉冬切除等方法均能有效衰减部分多次波,但对于强反射界面(如硬海底,海底以下高速层等)产生的多次波,单独采用以上方法仍会产生能量较强的残余多次波,针对该问题,根据多次波与一次波周期性差异在线性拉东域比时空域更为明显且有利于压制长周期多次波的特点,可以将线性拉东域预测反褶积技术作为SRME、高精度拉东切除技术的补充,取长补短,联合应用消除深水海底多次波.本文对线性拉东域预测反褶积的基本原理进行了阐述,讨论了如何选取参数,并对其应用范围进行了分析.将线性拉东域预测反褶积应用于南海某凹陷海洋资料的多次波衰减处理中,有效的消除了残余海底多次波,验证了该方法的实用性.     

南海南部深部结构的复杂构造地震成像

南海南部海底崎岖、地下构造复杂、二维地震资料中多次波非常发育,在地震数据处理中很难获得令人满意的成像效果,为此本文采用了预测反褶积后F-K域滤波及叠前时间偏移后Radon变换去除多次波的组合方案,有效地衰减了多次波,突出了有效信号能量;采用速度分析点加密法获取了更为准确的速度,实现了陡倾角构造与海山的准确成像;采用百分...     

地表一致性反褶积在地震勘探中的应用及效果

地震资料的反褶积处理是通过改造地震激发子波,进而消除地震激发子波在传播过程中所受的虚反射、层间多次反射和大地滤波等影响的一种地震勘探资料处理方法.反褶积的方法很多,如:脉冲反褶积、预测反褶积、地表一致性反褶积.它们之间主要区别之一在于对地震子波的假设和估计地震子波的方法.所以在处理过程中应根据不同的区域资料特征采取不同的反褶积方法.本文以河南省某煤预查区地震勘探为实例,着重总结和比较地表一致性反褶积技术在地震资料处理中的应用效果.应用研究表明,适当选择时窗和自相关步长进行自相关分析,地表一致性反褶积能够展宽频谱,压缩地震子波,并能校正地震信号的相位谱,输出零相位子波,可以较大程度地提高地震资料的分辨率,提高勘探能力.     

南黄海中古生代残留沉积地层的地震成像

文中展示了南黄海中古生代沉积地层的成像结果,显示南黄海隆起区存在较稳定的中古生代沉积地层残留区.南黄海中古生代残留沉积地层的成像问题一直是地震资料处理的难点,主要原因是以海底鸣震为主的层间多次波和与中古生界顶界面相关的长程多次波非常发育,致使有效反射波无法准确识别.本文利用平面波域层间多次波的可预测性,对地震数据平面波分解之后进行预测反褶积,较好地消除了层间多次波.利用速度差异分离长程多次波,在有效波速度无法确定的条件下,先从低速的、可以确认的多次波开始进行分离,逐步确认并分离多次波,同时逐步进行有效波的确认和成像速度的尝试,最终得到合适的成像速度场和多次波被充分分离的数据.成像结果对南黄海残留盆地研究有一定借鉴意义.     

多次波压制技术在南海北部陆缘深水区的应用

南海北部洋陆过渡带水深变化剧烈,海底地形崎岖,构造复杂多变,地震勘探中多次波十分发育,常规处理难以取得很好的成像效果.由于方法的局限性,任何一种多次波压制方法都很难彻底地压制全部多次波.本文通过对国家基金委共享航次在南海北部陆缘采集的二维反射地震资料的分析,总结了不同多次波压制技术的优缺点和适用范围,结合滤波类方法和波动理论类方法,综合应用了表面相关多次波衰减法(SRME,Surface Related Multiple Elimination)、τ-p域预测反褶积、高分辨率抛物线Radon变换以及叠后压制残余多次波等方法去除多次波,有效突出了一次反射波,证明了方法的有效性和实用性,同时也形成了针对南海深水资料多次波压制的一种有效的处理流程和方法组合.     

抛物线Radon变换方法在压制地震多次波中的应用

多次波的存在会降低地震资料的信噪比,影响地震资料处理效果和后续的地质解释精度,压制多次波干扰是地震资料处理中的一个重要环节。利用有限差分方法正演模拟几个典型地质模型的地震波场,之后进行动校正,再通过抛物线Radon变换方法把t-x域一次反射波和多次波变换到τ-q域,由于存在速度的差异,Radon域记录中的一次反射波和多次波的能量互相分开,在Radon域对多次波的能量切除,保留反射波的有效信息,达到压制多次波的效果,提高资料的信噪比。通过对实际地震数据进行抛物线Radon变换结果也进一步验证了该方法在压制多次波中的应用效果。      

高分辨率数据处理技术在近海工程地震勘探中的应用

在近海工程多道浅地层地震勘探中,高分辨率地震数据处理技术是一个重要的组成部分,也是提高资料分辨率和信噪比的有效途径.由于为工程服务的近海地震勘探对分辨率的要求很高,因此地层速度的拾取和衰减多次反射干扰波成为影响分辨率的主要因素.本文针对我国某海域工程地震资料,运用了高分辨率速度分析技术和K—L变换多次波压制方法,极大地提高了资料的信噪比和分辨率,有效地解决了该地区的工程地质问题.     

海上气田浅层气下覆复杂构造地区Q值补偿及逆时偏移成像——以巴西桑托斯盆地应用实例

针对海上资料特点及研究区上覆碳酸盐岩地层成像差、多次波发育、受浅层气影响的下层构造成像模糊,振幅异常及能量损失、地震反射信号弱、信噪比低的难点,以叠前综合去噪、预测反褶积、CMP道集高精度Radon变换衰减层间多次波为核心前期处理方法,结合基于黏弹性介质地震波衰减模型的Q值补偿技术,利用逆时偏移对复杂断层进行了精确成像,处理后的地震剖面较好地反映出了砂岩、碳酸盐岩储层的形态和构造特征,目的层断块、断点清楚,证明设计流程及方法可以较好地解决该地区目的层成像难的问题,并可应用于相似浅层气地区成像.     

逆散射级数法预测层间多次波的算法改进及其策略

波动方程预测层间多次波有两类算法.一是由Berkhout和Verschuur（1998）提出的基于CFP（共聚焦点）延拓的算法,另一是由Weglein（1997）等人提出的基于逆散射级数法的算法.ISS（逆散射级数）算法具有不依赖速度模型的优点,但是其计算成本很高.本文改进了Weglein提出的1-D预测公式,提高了计算效率.从理论方面,改进算法的计算速度提高了大约12倍.实际结果表明其计算速度更快,提高近80倍左右,主要原因在于后者的空间复杂度小.此外,由于时间-空间域比伪深度-波数域的层间多次波预测有更多的优点,本文推导了1.5-D时间-空间域ISS层间多次波预测算法,该算法有计算速度快、预测噪音小、适应观测系统能力强和不依赖地震子波等优点.在实际应用中,计算成本高是ISS预测层间多次波算法的主要障碍.因此,本文还讨论了预测层间多次波的策略.结合本文提出的改进算法和应用策略,能够在实际地震资料处理中较好地应用ISS预测层间多次波.     

平面波域反数据处理压制多次波方法研究

在地震勘探领域,尤其是海洋地震勘探中,多次波一直是影响地震处理与解释的主要因素之一.本文基于x-t域的反馈模型及多次波衰减的理论,详尽推导了x-t域反数据处理的方法;基于平面波域多次波的产生机制,借鉴x-t域反数据处理的方法,推导出在平面波域进行反数据处理的原理.文中给出一个有限差分的模拟数据进行测试,处理结果表明,本方法可以有效地衰减表层相关多次波,提高地震数据分析的精确性,在保护了一次波能量的同时,可以更加有效快捷地去除多次波.     

SLOPE DATA WATER-BOTTOM MULTIPLE ATTENUATION1

Methods for predicting and attenuating water-bottom multiples by wavefield extrapolation have been discussed by several investigators. Because these prediction methods operate on shot records, boundary conditions must be specified for every shot record. The approach presented operates in the common-offset plane; a model of expected water-bottom multiples is generated from the observed surface wavefield using a finite-difference wave-equation migration algorithm with an offset term. An accurate water-depth profile is required, but there is no restriction on the shape of the water bottom other than a dip limit of approximately 18–20°. In generating a multiple model, the water-bottom primary and each water-bottom multiple reflection of the observed surface wavefield are extrapolated to a higher order. Thus, the extrapolated water-bottom primary of the model is lined up with a water-bottom multiple in the data and each multiple in the model is lined up with a higher-order (or later) multiple in the data. Prestack multiple attenuation is achieved, for one offset at a time, by first adapting the model of expected multiples to the observed data and then subtracting the predicted multiple energy. An error-constrained adaptation algorithm is proposed in order to control instabilities. No assumptions are made about primary reflections and no subwater-bottom velocities are required. Computational efficiency of modelling and adaptation can be improved by applying this method only to near and intermediate offsets as the stacking process usually provides sufficient multiple attenuation at far offsets. A field data example demonstrates the potential of the proposed method for improving the primary-to-multiple ratio in prestack and post-stack data.     

Migration of multiples from the South China Sea

The South China Sea where water depth is up to 5000 m is the most promising oil and gas exploration area in China in the future.The seismic data acquired in the South China Sea contain various types of multiples that need to be removed before imaging can be developed.However,compared with the conventional reflection migration,multiples carry more information of the underground structure that helps provide better subsurface imaging.This paper presents a method to modify the conventional reverse time migration so that multiple reflections can migrate to their correct locations in the subsurface.This approach replaces the numerical impulsive source with the recorded data including primaries and multiples on the surface,and replaces the recorded primary reflection data with multiples.In the reverse time migration process,multiples recorded on the surface are extrapolated backward in time to each depth level,while primaries and multiples recorded on the surface are extrapolated forward in time to the same depth levels.By matching the difference between the primary and multiple images using an objective function,this algorithm improves the primary resultant image.Synthetic tests on Sigsbee2 B show that the proposed method can obtain a greater range and better underground illumination.Images of deep water in the South China Sea are obtained using multiples and their matching with primaries.They demonstrate that multiples can make up for the reflection illumination and the migration of multiples is an important research direction in the future.     

利用多次波提高地震数据分辨率方法

目前对多次波的有效利用仅围绕多次波成像技术展开,通过成像多次波试图获取更丰富的地下构造信息.不同于该思路,本文另辟蹊径,从利用多次波提高地震数据分辨率角度出发,对多次波的有效利用进行了深入挖掘.首先基于聚焦变换思想在聚焦域内实现多次波的降阶,通过理论推导得出聚焦域内多次波表现为原始数据的多维子波反褶积这一重要结论,从理论上证明了本文方法提高地震数据分辨率的可行性;然后采用引入整形正则化的非稳态回归自适应匹配滤波方法将聚焦域内由多次波构建的高分辨率数据分离出来,实现原始数据的高分辨率转换.与常规反褶积模型不同,该方法基于波动理论推导得出,可以适用于任意复杂情况;每一道输出结果中所有炮记录都参与了运算,从空间上加以约束,在提高纵向分辨率的同时可以改善数据的横向分辨率.最后通过模型试算和实际资料处理对本文方法的有效性、适应性和实用性进行了验证.     

一阶多次波聚焦变换成像

将多次波转换成反射波并按传统反射波偏移算法成像,是多次波成像的一种方法.聚焦变换能准确的将多次波转换为纵向分辨率更高的新波场记录,其中一阶多次波转换为反射波.本文对聚焦变换提出了两点改进:1)提出局部聚焦变换,以减小存储量和计算量,增强该方法对检波点随炮点移动的采集数据的适应性;2)引入加权矩阵,理论上证明原始记录的炮点比检波点稀疏时,共检波点道集域的局部聚焦变换可以将多次波准确转换成炮点与检波点有相同采样频率的新波场记录.本文在第一个数值实验中对比了对包含反射波与多次波的原始记录做局部聚焦变换和直接对预测的多次波做局部聚焦变换两种方案,验证了第二种方案转换得到的波场记录信噪比更高且避免了第一个方案中切聚焦点这项比较繁杂的工作.第二个数值实验表明:在炮点采样较为稀疏时,该方法能有效的将一阶多次波转换成反射波;转换的反射波能提供更丰富的波场信息,成像结果更均衡、在局部有更高的信噪比,以及较高的纵向分辨率.     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.     

RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA

Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection.     

A NEW METHOD OF CANCELLING MULTIPLE REFLECTIONS: “SOUSTON”*

Cancelling of multiple reflections by CDP stacking is generally incomplete. When the order of coverage is low, an improvement may be obtained by the use of special filters (e.g. Schneider et al: Geophysics 1965; D'Hoeraene: Geophysics 1966). But when the order of coverage is high, the efficiency of these filters is less obvious. On the other hand, the higher the order of coverage, the higher the efficiency in the solution presented. Schematically, our method includes three steps: Firstly: Searching for the multiple reflections. For this purpose Move Out corrections corresponding to the multiples are applied to each individual record, then CDP stacking is carried out. Secondly: Cleaning up of the individual records. “Decorrections” are firstly applied, then multiples are subtracted. For this, subtraction is weighed versus the energy of the multiples, that is to say versus the correlation between the original record and the pure multiple reflections. Thirdly: Normal processing with the cleaned records which have been obtained. The different steps of this method are illustrated with the help of theoretical examples. Its efficiency improvement in regard to normal CDP stacking is then demonstrated by means of true examples.     

Multiple suppression by 2D filtering in the parabolic τ–p domain: a wave-equation-based method1

The filter for wave-equation-based water-layer multiple suppression, developed by the authors in the x-t, the linear τ-p, and the f-k domains, is extended to the parabolic τ-2 domain. The multiple reject areas are determined automatically by comparing the energy on traces of the multiple model (which are generated by a wave-extrapolation method from the original data) and the original input data (multiples + primaries) in τ-p space. The advantage of applying the data-adaptive 2D demultiple filter in the parabolic τ-p domain is that the waves are well separated in this domain. The numerical examples demonstrate the effectiveness of such a dereverberation procedure. Filtering of multiples in the parabolic τ-p domain works on both the far-offset and the near-offset traces, while the filtering of multiples in the f-k domain is effective only for the far-offset traces. Tests on a synthetic common-shot-point (CSP) gather show that the demultiple filter is relatively immune to slight errors in the water velocity and water depth which cause arrival time errors of the multiples in the multiple model traces of less than the time dimension (about one quarter of the wavelet length) of the energy summation window of the filter. The multiples in the predicted multiple model traces do not have to be exact replicas of the multiples in the input data, in both a wavelet-shape and traveltime sense. The demultiple filter also works reasonably well for input data contaminated by up to 25% of random noise. A shallow water CSP seismic gather, acquired on the North West Shelf of Australia, demonstrates the effectiveness of the technique on real data.     

衰减干涉假象的海洋一次波与多次波同时成像方法

海水与空气间的强波阻抗差使得海洋地震资料普遍发育自由表面相关多次波,如何利用好多次波所携带的有效信息已成为提高海洋地震资料成像品质的新突破点.基于面炮偏移的一次波与多次波同时成像方法能够避免多次波预测精度的影响,但是,正向传播的震源子波与反向延拓的自由表面相关多次波所产生的干涉假象严重制约了该技术的应用,本文提出了一种基于单程波偏移算子,可在成像域压制干涉假象的一次波与多次波同时成像方法.其中包含了三个步骤:第一,传统单程波偏移成像方法中的震源子波替换为一次波、多次波与震源子波,初始上行延拓波场为一次波与多次波,基于单程波算子的波场延拓与互相关成像条件的应用得到包含干涉假象的一次波与多次波同时成像;第二,以子波为震源,自由表面相关多次波为记录,按照传统单程波偏移成像方法得到干涉假象;第三,基于最小二乘匹配滤波算法,将第一步的成像结果与第二步的干涉假象进行匹配相减,得到干涉假象衰减后的一次波与多次波同时成像,避开了由于实际资料子波无法准确提取而造成一次波与多次波对成像能量级的不一致性.Sigsbee2B模型测试验证了本方法的有效性,并在我国某探区深海实际资料处理中得到了成功应用,深层基底得到了清晰刻画,并且照明均衡度明显改善.     

基于波动方程表面多次波预测与 自适应相减方法研究

多次波预测与自适应相减是基于波动方程表面多次波压制的两个重要环节.文中利用具有并行计算优势的GPU加速表面多次波预测,使得预测效率大为提高.在自适应相减算法中,文中将预测的多次波道、预测多次波道的Hilbert变换道、预测多次波道的高频重建道、以及它们的平移道用作自适应相减中的多次波模型道.Hilbert变换道用以补偿预测多次波的相位信息,高频重建道用以改善预测多次波的高频信息,补偿频带能量差异.文中在预测和相减过程中均采用迭代算法,迭代预测,可较好地获得多次波的运动学特性,迭代相减,可较好地获得多次波的动力学特性,迭代预测与相减使预测的多次波与地震数据中实际的多次波更好地匹配.将该方法应用于理论模拟的SMAART模型和实际海洋数据中,测试结果表明,该方法预测多次波效率较高,在保持有效波振幅条件下可有效地压制地震数据中的表面多次波.