Kirchhoff型真振幅反偏移加权函数分析

Kirchhoff型真振幅反偏移的基本目的是通过等时线叠加得到偏移的输入场,是将深度偏移的结果再反映射到时间域去,而真振幅加权函数是Kirchhoff型真振幅反偏移的核心内容,它没有唯一的数学形式,目前发表在文献当中的加权函数公式是在物理直观下的最优表述.本文在前人研究的基础上推导了真振幅反偏移加权函数公式,对影响加权...     

逆时偏移与Kirchhoff偏移的抽道集对比研究

共成像点道集(CIG)可用于速度建模、偏移质量控制以及地下属性解释等。对于复杂碳酸盐岩储层,叠前AVO反演技术在储层识别、裂缝预测和流体检测方面的作用日益显著,叠前反演的精度很大程度上依赖于共成像点道集的AVO特性是否符合地下介质的真实地震响应。因此,叠前偏移不仅要实现反射点的准确归位,还必须得到振幅保真的共成像点道集,从而为AVO/AVA反演提供保幅的共成像点道集输入数据。本文在理论上对真振幅Kirchhoff叠前深度偏移和逆时偏移(RTM)抽道集技术的保幅性进行分析:真振幅Kirchhoff偏移抽道集通过改变加权函数实现振幅保持,仅适用于横向速度变化不大的情况;基于双程波动方程的逆时偏移采用互相关成像条件抽取共成像点道集,成像精确且无倾角限制,相对保幅性优于前者,适应任意复杂速度模型。模型和实例测试结果表明:在简单构造区域,真振幅Kirchhoff保幅叠前深度偏移与保幅逆时偏移抽取的共成像点道集都能起到良好的保幅效果;但在复杂构造区,基于逆时偏移抽取的共成像点道集,保留了全部波场的路径与振幅信息,其保幅性与成像精度优于Kirchhoff保幅叠前深度偏移。因此,针对探区的复杂程度差别及勘探精度的不同要求,应选择不同的抽道集方法,保证振幅精度的同时,兼顾生产效率。      

面炮成像、控制照明与AVA道集

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点.面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题.本文采用面炮成像技术进行叠前深度偏移,通过对面炮震源下行波场的质量控制和优选射线个数和范围,以达到最佳的成像效果.采用控制照明技术,较大地提高了目标地层的成像精度.与此同时,得到振幅随入射角变化(AVA)道集,有利于叠前振幅解释和储层岩性预测.数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算.     

起伏地表稳相偏移方法及其在准噶尔盆地高密度地震勘探中的应用

稳相偏移方法为基于菲涅耳带的散射叠加偏移方法.叠前时间偏移实践中,稳相偏移对振幅保真、偏移噪声压制、陡倾角构造正确成像以及计算效率提升均有较大帮助.基于浮动基准面假设,提出了起伏地表倾角道集计算方法.菲涅耳带可直观的展示于这一偏移域倾角道集中.采用人机交互拾取菲涅耳带边界后,可在起伏地表叠前时间偏移中结合地震道在成像点的倾角计算,剔除菲涅耳带之外的噪声干扰,从而实现起伏地表稳相偏移技术流程.准噶尔盆地阜东斜坡某试验区的单点高密度数据处理试验表明,起伏地表稳相偏移方法由于提供了空变、时变的偏移孔径,因而可兼顾陡倾角构造成像和偏移噪声的压制.进一步的分析表明,稳相偏移方法可提高单点高密度采集数据高频段信号的信噪比,为结合叠前时间偏移补偿介质对高频信号的吸收效应,提高地震成像的分辨率,进而发挥单点高密度采集数据集的宽频带优势提供重要辅助作用.     

基于旅行时梯度场的Kirchhoff PSDM角度道集生成方法研究（英文）

角度域共成像点道集(ADCIGS)是偏移速度分析和振幅随角度变化分析(AVA)的基础数据。传统Kirchhoff叠前深度偏移(KPSDM)按偏移距组织数据,能方便的输出偏移距域共成像点道集(ODCIGS),其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为:(1)利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;(2)根据旅行时场的梯度方向计算反射张角;(3)在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是:(1)按照炮数据来组织输入数据;(2)旅行时场的输入与单炮覆盖范围相联系以节省内存;(3)多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。     

Kirchhoff积分叠前时间偏移全方位角度道集生成方法与实现

方位角度域共成像点道集能够客观反映地下介质的速度、各向异性参数异常以及振幅随角度变化（AVA）和裂缝信息。传统Kirchhoff PSTM通常输出偏移距域共成像点道集,对于速度分析、各向异性分析、AVA分析、裂缝识别等均存在诸多不便。本文提出了基于走时梯度的Kirchhoff叠前时间偏移全方位角度集输出方法并提出工业上切实可行的实现方案。通过走时场梯度计算波场传播方向矢量,形成能够反映观测系统参数和波场传播情况的全方位角度域共成像点道集。为了在大规模地震数据Kirchhoff积分叠前时间偏移中输出全方位角度道集,本文给出基于输入道方式的偏移实现方法,采用逐条inline线进行线偏移成像,从而大大降低了全方位角度道集输出对计算机内存的压力,显著提高了Kirchhoff积分时间偏移输出全方位角度道集的可行性。三维盐丘模型测试和海上某区块三维实际资料试验证明了本文方法的正确性。      

共炮检距道集波动方程保幅叠前深度偏移方法

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.     

检波点与震源能量比值在偏移成像条件中的作用分析

利用检波点能量与震源能量比值构造衰减系数,研究波场各频率成分在自相关与反褶积两类成像条件中偏移成像的作用.分析检波点能量与震源能量比值分布特征,有利于克服反褶积型成像条件中分母趋于零时计算不稳定现象,也有助于平滑检波点波场与震源波场;检波点能量与震源能量比值超出正常范围容易导致偏移结果中出现照明不均匀现象.本文方法有助于改善反褶积型成像条件的偏移计算,也有助于应用照明技术指导偏移成像结果的振幅校正.     

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

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

逆时偏移与Kirchhoff偏移的对比研究

Kirchhoff偏移和逆时偏移(RTM)是目前最广泛应用的两种叠前深度偏移方法.两种方法的基础均为波动方程,理论上没有倾角限制.抽取的共成像点道集(CIG)可用于速度建模、叠前反演以及地下属性解释等.叠前反演的精度依赖于共成像点道集的AVO特性是否符合地下介质的真实地震响应.因此,叠前偏移不仅要实现反射点的准确归位,得到的共成像点道集必须要振幅保真.本文在理论上阐述两种偏移方法的基本原理,并对抽道集技术的保幅性进行对比.Marmousi和Sigsbee2a模型结果表明:在简单构造区,两者的成像效果相当;在复杂构造区,不论在成像精度还是保幅性上逆时偏移都比Kirchhoff偏移优越,特别是高陡倾角断层、岩下构造等复杂构造成像时.然而,在碳酸盐岩储层实际资料处理时发现,由于速度场准确性的限制,两种偏移方法各有优势.因此,在实际应用时,由于无法明确得到地下介质的确切信息,较难对两种偏移方法给出明确的优劣评价.     

The stationary phase analysis of the Kirchhoff-type demigrated field

In the literature, stationary phase analysis of Kirchhoff-type demigrated fields is carried out mainly under the following two conditions: (1) The considered isochrone and the target reflector are tangential to each other; (2) The spatial duration of the wavelet of the depthmigrated image is short. For the isochrones that are not tangential to the target reflector and for the depth-migrated images that have a large spatial duration, the published stationary phase equation for the demigrated field will become invalid. For performing the stationary phase analysis of the Kirchhoff-type demigrated field under the conditions that the considered isochrone and the target reflector are not tangential to each other and that the spatial duration of the wavelet of the depth-migrated image is not short (the general conditions), I derive the formulas for the factors appearing in the stationary phase formula in two dimensions, from which I find that for different isochrones the horizontal coordinates of the stationary point of the depth difference function are different. Also, the equation for the Kirchhoff-type demigrated field consists of two parts. One is the true-amplitude demigrated signal and the other is the amplitude distortion factor. From these facts the following two conclusions can be drawn: (1) A demigrated signal is composed of many depth-migrated images and one depth-migrated image trace provides only one sample to the demigrated signal; and (2) The amplitude distortion effect is an effect inherent in the Kirchhoff-type demigration and cannot be eliminated during demigration. If this effect should be eliminated, one should do an amplitude correction after demigration.     

Time-related capture zones for radial flow in two dimensional randomly heterogeneous media

We consider the effect of randomly heterogeneous hydraulic conductivity on the spatial location of time-related capture zones (isochrones) for a non-reactive tracer in the steady-state radial flow field due to a pumping well in a confined aquifer. A Monte Carlo (MC) procedure is used in conjunction with FFT-based spectral methods. The log hydraulic conductivity field is assumed to be Gaussian and stationary, with isotropic exponential correlation. Various degrees of domain heterogeneity are considered and stability and accuracy of the MC procedure is examined. The location of an isochrone becomes uncertain due to heterogeneity, and it is strongly influenced by hydraulic conductivity variance. The probability that a particle released at a point in the aquifer is pumped by the well within a given time is identified. We propose a new expression for the probabilistic spatial distribution of isochrones, which is formally similar to the analytical solution for a uniform medium and takes into account the effects of heterogeneity.     

偏移成像中的垂直分辨率

Kirchhoff真振幅偏移研究表明,偏移对地震子波具有拉伸作用,拉伸程度与地震波速度、反射角和地层倾角有关.子波拉伸将改变子波宽度(子波延续长度),对垂直分辨率有重要影响.本文通过对衰减Sinc子波的实际计算,分析了地震波速度、反射角和地层倾角分别与子波宽度的关系.从分析来看,陡倾角和大反射角（大炮检距）造成偏移成像垂直分辨率严重下降.高主频子波对各参数均有较低敏感度.     

How rough sea affects marine seismic data and deghosting procedures

Most seismic processing algorithms generally consider the sea surface as a flat reflector. However, acquisition of marine seismic data often takes place in weather conditions where this approximation is inaccurate. The distortion in the seismic wavelet introduced by the rough sea may influence (for example) deghosting results, as deghosting operators are typically recursive and sensitive to the changes in the seismic signal. In this paper, we study the effect of sea surface roughness on conventional (5–160 Hz) and ultra‐high‐resolution (200–3500 Hz) single‐component towed‐streamer data. To this end, we numerically simulate reflections from a rough sea surface using the Kirchhoff approximation. Our modelling demonstrates that for conventional seismic frequency band sea roughness can distort results of standard one‐dimensional and two‐dimensional deterministic deghosting. To mitigate this effect, we introduce regularisation and optimisation based on the minimum‐energy criterion and show that this improves the processing output significantly. Analysis of ultra‐high‐resolution field data in conjunction with modelling shows that even relatively calm sea state (i.e., 15 cm wave height) introduces significant changes in the seismic signal for ultra‐high‐frequency band. These changes in amplitude and arrival time may degrade the results of deghosting. Using the field dataset, we show how the minimum‐energy optimisation of deghosting parameters improves the processing result.     

Focusing in Prestack Isochrone Migration Using Instantaneous Slowness Information

—Prestack migration finds increasing application in processing crustal seismic data. However, less effort has been made to incorporate slowness information in the imaging process. The combination of slowness information with migration leads to an improved image in the depth domain, especially by reducing migration artefacts and noise. A slowness-driven isochrone migration scheme is introduced for migration of 2-D seismic data. Instantaneous slowness information p(x, t) is extracted from the data using correlation analysis in moving time and space windows. Slowness values resulting from spatial coherent energy (signal) and incoherent background noise are distinguished by the simultaneous evaluation of an instantaneous coherence criterion g(x, t). In slowness-driven isochrone migration this information is used for locally weighting the amplitude A(x, t) smearing on the isochrone surface. In particular, slowness p and coherence criterion g determine position and sharpness of a Gaussian weighting function. The method is demonstrated using two synthetic data examples and is subsequently applied to two deep crustal data sets, one wide-angle (along DEKORP4) and one steep-angle reflection seismic observation (KTB8506). Both data sets were collected in the surroundings of the KTB drill site, Oberpfalz, as part of the German DEKORP project.     

广义付里叶变换及其在磁共振成像中的应用

许多共振成像（ＭＲＩ）的应用,要求人体组织磁化强度的空间分布Ｆｏｕｒｉｅｒ变换的幅度和相位的精确分离,及利用正交双通道采集时间和空间高分辨图象序列。在传统的基于Ｆｏｕｒｉｅｒ变换的成象方法中,图象序列从一帧空间编码独立地重建得到,因此给定空间分辨率的情况下,每帧空间编码的数量限制了图象序列的时间分辨力;又常通过牺牲空间分辨率来提高图象序列的空间分辨录。由于传统Ｆｏｕｒｉｅｒ变换的核函数是线性相位,若变换ＭＲＩ技术必须很好地满足主磁场均匀性和梯度场线性的条件,因此传统Ｆｏｕｒｉｅｒ变换成象方法对临床应用中常出现ＭＲ信号相位失真极为敏感。本文提出：利用主磁场中组织的磁化强度不变,能很好消除图象的相位失真。     

复小波包域的微震信号分析提取方法

提出了一种在复小波包域分析提取微震信号的新方法。此方法先采用复小波包分解同时获得信号的幅值信息和相位信息,再采用一种称为“比值加权”的复小波包重构法,重构出信噪比大大增强的微震信号。实际应用表明这种新方法在提取微震信号,提高信噪比方面效果显著。     

Enhancement of Signal-to-noise Ratio in Natural-source Transient Magnetotelluric Data with Wavelet Transform

—For audio-frequency magnetotelluric surveys where the signals are lightning-stroke transients, the conventional Fourier transform method often fails to produce a high quality impedance tensor. An alternative approach is to use the wavelet transform method which is capable of localizing target information simultaneously in both the temporal and frequency domains. Unlike Fourier analysis that yields an average amplitude and phase, the wavelet transform produces an instantaneous estimate of the amplitude and phase of a signal. In this paper a complex well-localized wavelet, the Morlet wavelet, has been used to transform and analyze audio-frequency magnetotelluric data. With the Morlet wavelet, the magnetotelluric impedance tensor can be computed directly in the wavelet transform domain. The lightning-stroke transients are easily identified on the dilation-translation plane. Choosing those wavelet transform values where the signals are located, a higher signal-to-noise ratio estimation of the impedance tensor can be obtained. ? In a test using real data, the wavelet transform showed a significant improvement in the signal-to-noise ratio over the conventional Fourier transform.     

SOURCE ARRAY SCALING FOR WAVELET DECONVOLUTION*

A seismic source array is normally composed of elements spaced at distances less than a wavelength while the overall dimensions of the array are normally of the order of a wavelength. Consequently, unpredictable interaction effects occur between element and the shape of the far field wavelet, which is azimuth-dependent, can only be determined by measurements in the far field. Since such measurements are very often impossible to make, the shape of the wavelet—particularly its phase spectrum—is unknown. A theoretical design method for overcoming this problem is presented using two scaled arrays. The far field source wavelets from the source arrays have the same azimuth dependence at scaled frequencies, and the far field wavelets along any azimuth are related by a simple scaling law. Two independent seismograms are generated by the two scaled arrays for each pair of source-receiver locations, the source wavelets being related by the scaling law. The technique thus permits the far field waveform of an array to be determined in situations where it is impossible to measure it. Furthermore it permits the array design criteria to be changed: instead of sacrificing useful signal energy for the sake of the phase spectrum, the array may be designed to produce a wavelet with desired amplitude characteristics, without much regard for phase.