海上多次波的联合衰减法

地震勘探尤其是海上地震勘探中存在着各类多次波,由于多次波的存在严重影响了速度分析、叠加、偏移成像等地震资料处理,海上多次波主要有全程多次波和层间多次两大类,为了压制海上不同类型的多次波,本文首先分析波场外推、预测反褶积和拉东变换衰减不同多次波的理论基础,然后联合采用这些方法,分别衰减全程多次波和层间多次波;即：首先对炮记录或者接收点记录进行波场外推,建立海底多次波模型,预测并减去全程多次波,然后利用预测反褶积衰减掉周期性明显的多次波,最后将数据转换到τ-p域,用拉冬变换根据同一时间多次波和有效波在速度等方面的差异,进一步分离层间多次波和剩余的全程多次波,并在该域中切除分离出的多次波,从而实现联合多次波衰减处理.通过对悉尼海区和里海等实际地震资料的处理证明,文提出的联合多次波衰减方法在海洋地震资料的处理中有着广泛的应用,联合衰减多次波处理流程具有快速、简洁、易于实现的特点,经处理后的地震资料有效信号损失小、保福性好的特点,有利于速度分析、叠加、叠前偏移等的后续处理工作.     

深水陡坡带绕射多次波压制方法研究

在深水陡坡带和复杂海底构造区域,即使采用长电缆,由绕射产生的自由表面相关多次波因为传播方向的异变,仍然难以避免在自由表面发生下行反射的反射点位置位于电缆长度范围之外,造成预测多次波的近偏移距信息不足或者负偏移距信息缺失等问题,以致不能准确地预测和压制陡坡带的绕射多次波.本文在常规反馈迭代自由表面多次波压制方法的基础上做出相应改进,将邻炮的炮记录依照炮点、检波点互换原理引入到反馈迭代循环当中,补充绕射多次波预测所需要的近偏移距和负偏移距信息.经过Sigsbee2b模型和实际资料测试,达到较好压制深水陡坡带的自由表面相关绕射多次波的目地.     

改进的多道预测算子压制浅水多次波方法

在浅水情况下,由于观测数据中缺少近偏移距信息,水层多次波的压制面临挑战.利用多道预测算子压制水层多次波是浅水环境下压制多次波的重要方法之一,这种方法先从输入数据中估计出多道预测算子,再将预测算子和输入数据做褶积预测出水层相关多次波.然而,估计的多道预测算子很容易受噪声污染,从而影响多次波模型的精度.所以,我们提出了改进的多道预测算子压制浅水多次波方法.该方法先从数据中估计出多道预测算子,并利用估计的算子构建出精确的水层模型;然后,通过计算算子的走时信息、估计振幅信息、合成新算子三个步骤来修正原始的多道预测算子.修正的算子不仅不受噪声影响,还含有精确的走时信息、可靠的振幅信息;最后,该方法用修正的算子来预测多次波,并结合自适应相减,将预测的多次波从输入数据中去除.通过合成数据和实际资料的验证表明,相比于原始的多道预测算子压制多次波方法,改进的方法能够取得更好的压制效果.     

井间地震物理模拟研究

井间地震资料具有极高的分辨率,但是,其波场十分复杂,特别是有效一次反射波场能量较弱,且往往被管波等强相干干扰淹没,波场的识别与分离难度较大.地震物理模拟技术是认识复杂地震波场的有效手段.我们通过单一地质体模型、复杂地质体模型和真实井间模拟模型等多种逼真地质模型的物理模拟,揭示了井间地震观测下直达波、透射波、反射波、折射波,多次波和导波等各类波型的特征.物理模拟试验表明,在炮检域平面上分析直达波的信噪比与能量分布,可以使处理人员对旅行时的误差有更直观的认识.所以识别初至与拾取旅行时间时,不仅要在共震源点道集(CSG)上分析拾取,而且还要在共炮检距道集(COG)上分析拾取,这样就提高了直达波识别的可靠性与拾取时间的精度.     

COHERENCY WEIGHTING—AN EFFECTIVE APPROACH TO THE SUPPRESSION OF LONG LEG MULTIPLES*

Long leg multiples can be suppressed by a method which provides an alternative to weighted common-depth-point stacking and multichannel stacking filtering. The suppression is achieved by coherency weighting whereby the time-dependent weighting factor decreases as the semblance of the multiple reflections increases. The algorithm of the method is described. Its efficiency is discussed in relation to the input data and results of its application to marine seismic data are presented. For practical application, the stacking velocity of the multiples has to be known. As the process is based on stacking velocities, different types of multiples can be handled, for instance water-bottom multiples or internal multiples. The parameter analysis shows that the degree of multiple suppression can easily be controlled by adapting the parameters of the procedure to the field conditions. During the suppression of multiples, the primaries are saved according to the moveout differences between the two. The non-linear behaviour of the process causes signal suppression and distortion effects, which have to be corrected by AGC normalization and low-pass filtering. Among the various applications available, only the suppression of long leg water-bottom multiples is treated here. The results show that their suppression on the basis of moveout differences is efficient even when standard length streamers are used in regions with water depth of up to 1500 m and more, if the stacking velocity of the primaries is about 10 to 20% higher than that of the multiples. Even if those parts of the primaries which are masked by the multiples are suppressed in the individual common-depth-point gathers by the procedure, the remaining primaries in the AGV stacked section are largely uncovered by the multiple suppression.     

径向道变换压制相干噪声方法研究

径向道变换(Radial Trace Transform)是将地震道集振幅值从偏移距一双程旅行时坐标系变换到视速度-双程旅行时坐标系,通过这种坐标系的变换,使相干噪声与有效信号在视速度和频率方面都有效分离.本文在介绍RT变换基本原理基础上,分析了RT变换中两种常用插值方法及其特点.并利用对模拟地震资料的处理,证明了RT域模拟-相减法较其带通滤波法在相干噪声压制与反射信号保持方面具有明显优势.最后,根据噪声特点,通过选择合理RT滤波参数,对实际地震资料进行处理试验,获得了较好的去噪效果,明显提高了资料信噪比,验证了研究方法的有效性.     

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

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

地震多次波成像

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

MULTI-OFFSET ACOUSTIC INVERSION OF A LATERALLY INVARIANT MEDIUM: APPLICATION TO REAL DATA1

The aim of seismic inversion methods is to obtain quantitative information on the subsurface properties from seismic measurements. However, the potential accuracy of such methods depends strongly on the physical correctness of the mathematical equations used to model the propagation of the seismic waves. In general, the most accurate models involve the full non-linear acoustic or elastic wave equations. Inversion algorithms based on these equations are very CPU intensive. The application of such an algorithm on a real marine CMP gather is demonstrated. The earth model is assumed to be laterally invariant and only acoustic wave phenomena are modelled. A complete acoustic earth model (P-wave velocity and reflectivity as functions of vertical traveltime) is estimated. The inversion algorithm assumes that the seismic waves propagate in 2D. Therefore, an exact method for transforming the real data from 3D to 2D is derived and applied to the data. The time function of the source is estimated from a vertical far-field signature and its applicability is demonstrated by comparing synthetic and real water-bottom reflections. The source scaling factor is chosen such that the false reflection coefficient due to the first water-bottom multiple disappears from the inversion result. In order to speed up the convergence of the algorithm, the following inversion strategy is adopted: an initial smooth velocity model (macromodel) is obtained by applying Dix's equation to the result of a classical velocity analysis, followed by a smoothing operation. The initial reflectivity model is then computed using Gardner's empirical relationship between densities and velocities. In a first inversion step, reflectivity is estimated from small-offset data, keeping the velocity model fixed. In a second step, the initial smooth velocity model, and possibly the reflectivity model, is refined by using larger-offset data. This strategy is very efficient. In the first step, only ten iterations with a quasi-Newton algorithm are necessary in order to obtain an excellent convergence. The data window was 0–2.8 s, the maximum offset was 250 m, and the residual energy after the first inversion step was only 5% of the energy of the observed data. When the earth model estimated in the first inversion step is used to model data at moderate offsets (900 m, time window 0.0–1.1 s), the data fit is very good. In the second step, only a small improvement in the data fit could be obtained, and the convergence was slow. This is probably due to the strong non-linearity of the inversion problem with respect to the velocity model. Nevertheless, the final residual energy for the moderate offsets was only 11%. The estimated model was compared to sonic and density logs obtained from a nearby well. The comparison indicated that the present algorithm can be used to estimate normal incidence reflectivity from real data with good accuracy, provided that absorption phenomena play a minor role in the depth interval considered. If details in the velocity model are required, large offsets and an elastic inversion algorithm should be used.     

基于局部余弦基小波束的观测系统沉降法叠前深度偏移

将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子.本算法中,先对频率域的共点源和共点检波器道集做局部余弦小波束分解,然后分别沿共小波束源和共小波束检波器在深度方向延拓得到下一层波场.每个深度的波场,都等效于把源和检波器放在该层后所能接收到的地震记录,每点的像值由炮点和检波点重合时的零时刻波场值给出.通过二维SEG/EAGE盐丘模型和Marmousi模型的偏移成像结果验证该方法理论推导的正确性.另外,结果显示该方法继承了小波束域波场延拓在速度扰动较大情况下波传播及成像精度高的优点.     

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.     

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

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

利用微震数据对震源下方构造直接成像

现有的微震记录直接成像方法是将微震记录既当作入射记录,也当作散射记录,从而实现偏移成像.但此方法并不能突出透射波所携带的来自震源下方的深层散射波信息.本文在假设已知微震位置与子波的前提下,提出了对微震下方构造进行逆时偏移的成像方法.该方法类似于常规的逆时偏移,只是震源位置在地下.这使得在成像时,地下更深部的入射波场相比震源在地表时会更为精确,因此能够获得更加准确的成像结果.该方法会给成像结果带来一种尾波高频干扰：地下的震源发出的上行波与上方介质作用后,所产生的多级散射波会干扰反传波场.对此,在成像过程中,对入射场和散射场都进行左右行波分离,以压制该噪声.而在子波信息未知,无法重构入射场时,使用了激发时间成像条件,也能够实现同等效果的偏移成像,且不会出现尾波高频干扰.利用数值实验验证了本文方法的有效性.     

Internal multiple attenuation by Kirchhoff extrapolation

Surface‐related multiple elimination is the leading methodology for surface multiple removal. This data‐driven approach can be extended to interbed multiple prediction at the expense of a huge increase of the computational burden. This cost makes model‐driven methods still attractive, especially for the three dimensional case. In this paper we present a methodology that extends Kirchhoff wavefield extrapolation to interbed multiple prediction. In Kirchhoff wavefield extrapolation for surface multiple prediction a single round trip to an interpreted reflector is added to the recorded data. Here we show that interbed multiples generated between two interpreted reflectors can be predicted by applying the Kirchhoff wavefield extrapolation operator twice. In the first extrapolation step Kirchhoff wavefield extrapolation propagates the data backward in time to simulate a round trip to the shallower reflector. In the second extrapolation step Kirchhoff wavefield extrapolation propagates the data forward in time to simulate a round trip to the deeper reflector. In the Kirchhoff extrapolation kernel we use asymptotic Green's functions. The prediction of multiples via Kirchhoff wavefield extrapolation is possibly sped up by computing the required traveltimes via a shifted hyperbola approximation. The effectiveness of the method is demonstrated by results on both synthetic and field data sets.     

基于广义的炮偏移方法实现地表多次波和一次波联合成像

多次波是地下反射层的多次波反射,也蕴含了地下反射界面的信息,因此并不是绝对地只能被当做噪音来处理.为实现对地下构造的准确成像,本文基于广义概念上的炮偏移成像算法,对常规一次波偏移方法从用于向下延拓的上、下行场以及成像条件方面进行了改进,将同时含有表层多次波的炮记录与脉冲震源之和作为下行延拓的震源波场,将同时含有表层多次...     

Utilization of multiple scattering: the next big step forward in seismic imaging

Surface removal and internal multiple removal are explained by recursively separating the primary and multiple responses at each depth level with the aid of wavefield prediction error filtering. This causal removal process is referred to as “data linearization.” The linearized output (primaries only) is suitable for linear migration algorithms. Next, a summary is given on the migration of full wavefields (primaries + multiples) by using the concept of secondary sources in each subsurface gridpoint. These secondary sources are two‐way and contain the gridpoint reflection and the gridpoint transmission properties. In full wavefield migration, a local inversion process replaces the traditional linear imaging conditions. Finally, Marchenko redatuming is explained by iteratively separating the full wavefield response from above a new datum and the full wavefield response from below a new datum. The redatuming output is available for linear migration (Marchenko imaging) or, even better, for full wavefield migration. Linear migration, full wavefield migration, and Marchenko imaging are compared with each other. The principal conclusion of this essay is that multiples should not be removed, but they should be utilized, yielding two major advantages: (i) illumination is enhanced, particularly in the situation of low signal‐to‐noise primaries; and (ii) both the upper side and the lower side of reflectors are imaged. It is also concluded that multiple scattering algorithms are more transparent if they are formulated in a recursive depth manner. In addition to transparency, a recursive depth algorithm has the flexibility to enrich the imaging process by inserting prior geological knowledge or by removing numerical artefacts at each depth level. Finally, it is concluded that nonlinear migration algorithms must have a closed‐loop architecture to allow successful imaging of incomplete seismic data volumes (reality of field data).     

一阶多次波聚焦变换成像

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

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

A LAYER-STRIPPING TECHNIQUE FOR THE SUPPRESSION OF WATER-BOTTOM MULTIPLE REFLECTIONS*

A new method to suppress water-bottom multiples (water-bottom reverberations) uses the fact that in the domain of intercept time and ray parameter (τ–p domain) the water-bottom reverberations are strictly periodical for a horizontal flat sea bottom. Using this property a comb filter can be designed. The window of the filter should be approximately equal to the duration of a source pulse. The algorithm finds the maximum of the periodical energy throughout the τ–p domain and then designs the comb filter which eliminates the water bottom reverberations from each trace in the τ– p domain. This process can be repeated for higher order reverberations. Finally the τ–p domain with attenuated multiples is transformed back to the conventional x -- t space. The method is illustrated on a variety of synthetic data and on a set of real marine CMP data acquired in the North Sea near the Norwegian shore.     

LONG PERIOD SEA-FLOOR MULTIPLES AND THEIR SUPPRESSION*

Multiple sea-floor reflections in deep water often are not effectively suppressed by either CDP stacking nor standard predictive deconvolution methods. These methods fail because the reflection coefficient varies markedly with angle of incidence and also because of the variation of arrival time with offset and because of dip. For a reasonablly flat sea-floor, multiples of various orders and the primary sea-floor reflection which have all been reflected at nearly the same angle lie along a straight line through the origin in time-offset space. This line is called the “radial direction.” The multiples which lie along this line show a systematic relationship because they all experience the same water-bottom reflection effect. In other words, multiples behave in a stationary manner along the radial directions on multi-trace seismic records. A technique of multi-channel predictive deconvolution, called “Radial Multiple Suppression,” utilizes this aspect to design Wiener operators for the prediciton and suppression of water bottom multiples. The effectiveness of the technique is demonstrated by the study of field records, autocorrelations, velocity analyses, and stacked sections before and after Radial Multiple Suppression processing.