最小二乘反演迭代算法在压制海上变深度缆采集数据虚反射中的应用

海水面的虚反射(鬼波)引起海上拖缆采集数据陷波,导致地震记录频带变窄,而近年发展的变深度缆采集技术,具有多样的陷波特征,通过专门的去虚反射处理方法可获得宽频数据.本文基于已有研究成果,将最小二乘反演迭代压制虚反射算法应用于某海上变深度缆宽频处理.基于频率波数域镜像记录生成方法获得镜像炮集记录,并采用最小二乘解从变深度缆原始和镜像炮集记录中提取上行波.针对镜像炮集记录生成受初始速度模型精度的影响,使得某深度缆接收的上行波和下行波之间的实际延迟时间存在误差,采用最小二乘反演迭代算法最优化计算下行波与上行波之间的平均延迟时间和上行波记录,并采用时空数据窗口滑动克服延迟时间随炮检距和目的层深度变化问题.合成数据及某海上实际变深度缆数据处理测试结果表明,该方法能较好地压制变深度缆由海水面产生的虚反射,能达到拓宽地震记录频带目的.     

频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

海上倾斜缆采集技术具有多样的陷波特征,通过去鬼波处理可获得宽频数据.针对海水面波浪起伏及缆深误差引起的鬼波延迟时间估计误差以及崎岖海底和目的层深度变化使得鬼波和一次反射波的振幅差异系数随偏移距的变化而难以给定一个固定值的问题,本文推导出频率慢度域中鬼波滤波算子以及自适应迭代反演求解上行波算法,该鬼波滤波算子与不同水平慢度对应的鬼波和一次反射波的振幅差异系数以及鬼波延迟时间有关.并基于计算出的理论下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代最优计算该振幅差异系数和鬼波延迟时间.合成的及某海上采集的倾斜缆数据去鬼波处理结果表明,频率慢度域自适应迭代反演算法能较好地去除海上变深度缆鬼波,能达到拓宽地震记录频带目的.     

快速3D抛物Radon变换地震数据保幅重建

为解决3D AVO地震数据快速保幅重建问题,在传统3D抛物Radon变换的基础上提出一种3D快速高阶抛物Radon变换方法.该方法将传统抛物Radon变换与正交多项式相结合,通过正交多项式系数描述地震数据AVO信息,确保重建后的地震数据具有良好保幅效果.同时,该方法引入新变量λ_x=q_xf和λ_y=q_yf,通过对q_xf和q_yf的整体采样,消除了3D高阶抛物Radon变换算子对频率的依赖,使变换算子的求逆过程仅需计算一次,大大节省计算时间.理论模型和实际地震资料的处理结果表明,该方法重建效率高,保幅效果良好.     

基于波动方程的上下缆地震数据鬼波压制方法研究

本文发展基于波动方程的上下缆鬼波压制方法,推导了上下缆地震波场频率波数域波动方程延拓合并公式.基于Fourier变换的波场解析延拓确保上下缆资料振幅相位的一致性,消除了长拖缆远偏移距信号的计算误差,同时具有较高的计算效率;上下缆地震波场的波动方程法合并有效解偶鬼波干涉,实现综合利用上下缆地震数据压制鬼波.理论模型数据和实际采集地震数据的测试表明了方法的有效性.     

Evolution of deghosting process for single‐sensor streamer data from 2D to 3D

In marine acquisition, reflections of sound energy from the water–air interface result in ghosts in the seismic data, both in the source side and the receiver side. Ghosts limit the bandwidth of the useful signal and blur the final image. The process to separate the ghost and primary signals, called the deghosting process, can fill the ghost notch, broaden the frequency band, and help achieve high‐resolution images. Low‐signal‐to‐noise ratio near the notch frequencies and 3D effects are two challenges that the deghosting process has to face. In this paper, starting from an introduction to the deghosting process, we present and compare two strategies to solve the latter. The first is an adaptive mechanism that adjusts the deghosting operator to compensate for 3D effects or errors in source/receiver depth measurement. This method does not include explicitly the crossline slowness component and is not affected by the sparse sampling in the same direction. The second method is an inversion‐type approach that does include the crossline slowness component in the algorithm and handles the 3D effects explicitly. Both synthetic and field data examples in wide azimuth acquisition settings are shown to compare the two strategies. Both methods provide satisfactory results.     

基于保幅拉东变换的多次波衰减

为在去除多次波时有效保护地震一次反射波数据的AVO现象,给后续反演、解释提供准确的地震数据,本文提出了一种基于保幅拉东变换的多次波衰减方法,该方法是对常规抛物拉东变换的修改,把常规的稀疏拉东变换在拉东域分成两部分：一部分用于模拟零偏移距处的反射波能量,增加的另一部分用于模拟反射波振幅的AVO特性.该方法不仅考虑了反射波同相轴的形状,还考虑了反射波同相轴振幅幅度的变化,从而可把反射波信息进行有效转换,进而有利于多次波的消除,更好地恢复有效波的能量.在把地震数据由时间域转换到拉东域时,本文采用了IRLS算法实现保幅拉东算子的反演.模型数据和实际地震道集的试算分析表明,与常规拉东变换相比,保幅拉东变换在去除多次波的同时可有效保护一次反射波的AVO现象.     

Fast hyperbolic deconvolutive Radon transform using generalized Fourier slice theorem

The hyperbolic Radon transform has a long history of applications in seismic data processing because of its ability to focus/sparsify the data in the transform domain. Recently, deconvolutive Radon transform has also been proposed with an improved time resolution which provides improved processing results. The basis functions of the (deconvolutive) Radon transform, however, are time-variant, making the classical Fourier based algorithms ineffective to carry out the required computations. A direct implementation of the associated summations in the time–space domain is also computationally expensive, thus limiting the application of the transform on large data sets. In this paper, we present a new method for fast computation of the hyperbolic (deconvolutive) Radon transform. The method is based on the recently proposed generalized Fourier slice theorem which establishes an analytic expression between the Fourier transforms associated with the data and Radon plane. This allows very fast computations of the forward and inverse transforms simply using fast Fourier transform and interpolation procedures. These canonical transforms are used within an efficient iterative method for sparse solution of (deconvolutive) Radon transform. Numerical examples from synthetic and field seismic data confirm high performance of the proposed fast algorithm for filling in the large gaps in seismic data, separating primaries from multiple reflections, and performing high-quality stretch-free stacking.     

Three-parameter Radon transform in layered transversely isotropic media

Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.     

基于编码解码思想的鬼波预测与压制方法

高分辨率的宽带地震勘探技术是最近几年海上油气地震勘探的热点问题,鬼波压制是其中的核心议题.石油工业界,除了提出各种压制鬼波的采集方式外,资料处理过程中压制鬼波的方法也在不断地推陈出新.源、检鬼波的存在使得反射子波的有效频带变窄,成像分辨率降低,且干扰后续的自由表面多次波压制、FWI等.不同于常规的鬼波压制模型,本文基于编码与解码理论框架,用编码建立起鬼波预测模型,在Bayes反演框架下建立起解码估计一次波、从而压制鬼波的方法.基于此鬼波预测模型及相应反演理论的鬼波压制方法对鬼波的压制更为彻底.理论模型数据与实际资料测试结果验证了本文提出的理论框架和方法技术的有效性与优越性.     

Inversion-based deblending using migration operators

In this paper, we compare the denoising- and inversion-based deblending methods using Stolt migration operators. We use Stolt operator as a kernel to efficiently compute apex-shifted hyperbolic Radon transform. Sparsity promoting transforms, such as Radon transform, can focus seismic data into a sparse model to separate signals, remove noise or interpolate missing traces. Therefore, Radon transforms are a suitable tool for either the denoising- or the inversion-based deblending methods. The denoising-based deblending treats blending interferences as random noise by sorting the data into new gathers, such as common receiver gather. In these gathers, blending interferences exhibit random structures due to the randomization of the source firing times. Alternatively, the inversion-based deblending treats blending interferences as a signal, and the transform models this signal by incorporating the blending operator to formulate an inversion problem. We compare both methods using a robust inversion algorithm with sparse regularization. Results of synthetic and field data examples show that the inversion-based deblending can produce more accurate signal separation for highly blended data.     

混合域高分辨率抛物Radon变换及在衰减多次波中的应用

高分辨率Radon变换存在计算效率和分辨率不能兼得的困境.时间域算法可以获得很高的分辨率,但计算效率非常低;频率域算法具有良好计算效率,但分辨率不理想.为此发展了混合域高分辨率抛物Radon变换,即对频率域抛物Radon变换引入时变的稀疏权.本文给出了一种新的混合域高分辨率抛物Radon变换实现方法,并将该算法应用于叠前数据衰减多次波.文中给出了Radon变换和衰减多次波的流程.理论和实际数据算例表明本文方法既有较高的分辨率又有很高的计算效率.     

基于改进线性Radon变换的井孔地震上下行波场分离方法

反射波场分离是井孔地震资料处理中极其重要的一个环节,波场分离的质量直接影响成像结果的精度.不管是VSP还是井间地震资料,其反射波时距曲线都近似直线型,根据这一特征,本文提出一种改进的线性Radon变换方法来进行井孔资料的反射波上下行波场分离.该方法基于频率域线性Radon变换,通过引入一个新的变量λ来消除变换算子对频率的依赖性,避免了求取每一频率分量对应的不同变换算子,显著降低了计算成本;文中在求解该方法对应的最小二乘问题时,引入了发展较为成熟的高分辨率Radon变换技术来进一步提高波场分离的精度.采用本文方法进行井孔地震资料的上下行波场分离可以在保证分离精度的前提下有效地提高计算效率.根据上下行波在λ-f域内分布的特殊性,设计简单的滤波算子就可实现上下行波场的分离.最后通过合成数据试算以及实际资料处理（VSP数据和井间地震数据）验证了该方法的可行性和有效性.     

The simulation of far-field wavelets using frequency-domain air-gun array near-field wavelets

Air-gun arrays are used in marine-seismic exploration. Far-field wavelets in subsurface media represent the stacking of single air-gun ideal wavelets. We derived single air-gun ideal wavelets using near-field wavelets recorded from near-field geophones and then synthesized them into far-field wavelets. This is critical for processing wavelets in marine- seismic exploration. For this purpose, several algorithms are currently used to decompose and synthesize wavelets in the time domain. If the traveltime of single air-gun wavelets is not an integral multiple of the sampling interval, the complex and error-prone resampling of the seismic signals using the time-domain method is necessary. Based on the relation between the frequency-domain phase and the time-domain time delay, we propose a method that first transforms the real near-field wavelet to the frequency domain via Fourier transforms; then, it decomposes it and composes the wavelet spectrum in the frequency domain, and then back transforms it to the time domain. Thus, the resampling problem is avoided and single air-gun wavelets and far-field wavelets can be reliably derived. The effect of ghost reflections is also considered, while decomposing the wavelet and removing the ghost reflections. Modeling and real data processing were used to demonstrate the feasibility of the proposed method.     

广义Radon变换与叠前地震数据处理

在本文中,首先讨论了与几种地震层析成像对应的Radon变换公式,并导出了叠前地震记录的数学模型.在分析叠前地震记录与广义Radon变换的关系的基础上,讨论了速度分析、滤波、动校正、叠加等地震数据处理的数学物理意义.为展示广义Radon变换在地震数据处理中的应用,给出了用于滤波和消除多次波的方法及算例.     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

Deghosting towed streamer data in τ/p domain based on rough sea surface reflectivity

Currently, the deghosting of towed streamer seismic data assumes a flat sea level and a sea-surface reflection coefficient of ?1; this decreases the precision of deghosting. A new method that considers the rough sea surface is proposed to suppress ghost reflections. The proposed deghosting method obtains the rough sea surface reflection coefficient using Gaussian statistics, and calculates the optimized deghosting operator in the τ/p domain. The proposed method is closer to the actual sea conditions, offers an improved deghosting operator, removes the ghost reflections from marine towed seismic data, widens the bandwidth and restores the low-frequency information, and finally improves the signal-tonoise ratio and resolution of the seismic data.     

海上地震资料子波提取鬼波压制技术

海洋地震勘探由于受到海面强反射界面的影响,记录数据中存在鬼波,致使数据频谱中产生陷波效应,限制数据的频带宽度,影响数据的分辨率,干扰地层反射的识别.近年来发展的鬼波压制方法,能够较好的压制地震数据中的鬼波,然而现有鬼波压制技术容易产生延续相位,且计算效率较低.为进一步提高鬼波压制效果和工作效率,本文对海上地震资料特性进行分析,研究认为,海上地震资料中的海底反射信号真实记录了海底有效反射、震源鬼波和检波点鬼波的形态,可以利用提取的海底反射子波对数据进行确定性子波反褶积处理,能够有效的压制震源鬼波和检波点鬼波.模型海上实际数据测试结果表明,该方法具有鬼波压制效果好、计算效率高的优势,处理之后的数据频带得到有效拓宽,低频和高频信息均得到加强.     

Radon变换在非规则观测系统中波场分离的应用

Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了Radon变换在规则及其非规则观测系统中波场分离中应用,讨论了Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型的试算,阐明了Radon变换在缺失道、道距变化情况下依然可以高精度的分离波场,并对线性、抛物Radon变换进行了比较,对隧道超前预报反射波资料进行了处理,给出了隧道超前预报处理流程中提取反射波的算例.     

高分辨率Radon变换方法及其在地震信号处理中的应用

Radon变换方法在地震资料处理中广泛采用,在地震同相轴识别和估计方面具有良好效果．无论是倾斜叠加,还是广义Radon变换方法,一般采用最小二乘反演方法实现．目前,在提高反演算法的效率和分辨率方面仍值得研究．本文从倾斜叠加的定义出发,阐明Radon变换分辨率问题的来源和解决办法．采用最小二乘反演方法研究高分辨率抛物线Radon变换和双曲Radon变换时,给出稀疏约束预条件共轭梯度法求解的高分辨率Radon变换的实现方法,同阻尼最小二乘方法相比,分辨率和精度明显提高,文中给出了模型算例．根据有效波和多次波NMO后剩余时差不同,采用高分辨率抛物线和双曲Radon变换可以压制多次波,分别给出了方法原理,最后给出应用实例．研究表明,稀疏约束预条件共轭梯度法可以有效实现高分辨率Radon变换;数值算例表明,算法计算效率和精度较高,可以更好地实现多次波压制．     

多项式Radon变换

Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了次数为“2”的多项式Radon变换.讨论了多项式正反Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型试算，对多项式、线性、抛物线三种Radon变换进行了比较.对实际资料进行了多项式Radon变换处理，给出了消除地震反射记录中线性干扰的算例.