Kirchhoff型反偏移场稳相分析(英文)

在文献中,Kirchhoff型反偏移场的稳相分析主要是在下列两个条件下进行的:(1)等时面和目标反射面相切;(2)深度偏移像场信号的长度接近于零。对于与目标反射面不相切的等时面和长度远大于零的深度偏移像场子波,已有的结果将不再成立。为了在等时面和目标反射面不相切和深度偏移像场的子波长度远大于零的条件(一般条件)下对Kirchhoff型反偏移场进行稳相分析,我推导了出现在二维稳相分析公式中的诸因子的计算公式,并从中发现:(1)对于不同的等时面,距离差函数的稳相点具有不同的水平坐标;(2)Kirchhoff型真振幅反偏移的输出场由两部分(真振幅反偏移信号与振幅畸变因子)的乘积组成。由此得到下列两个结论:(1)一个给定的反偏移信号由多个深度偏移信号上的采样点组装而成,反偏移信号上的采样点个数等于对于这种组装有贡献的偏移信号的个数。(2)振幅畸变效应是Kirchhoff型反偏移中的固有效应,靠反偏移本身无法消除。如果一定要消除这种振幅畸变效应,必须对反偏移结果进行振幅校正。     

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子波的实际计算,分析了地震波速度、反射角和地层倾角分别与子波宽度的关系.从分析来看,陡倾角和大反射角（大炮检距）造成偏移成像垂直分辨率严重下降.高主频子波对各参数均有较低敏感度.     

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.     

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

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

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

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

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.     

Phase inversion deconvolution for long and short period multiples attenuationl1

A new approach to deconvolution has been developed to improve the attenuation of multiple energy. This approach to deconvolution is unique in that it not only eliminates the usual assumptions of a minimum phase lag wavelet and a random distribution of impulses, but also overcomes the noise limitation of the homomorphic deconvolution and its inherent instability to phase computation. We attempt to analyse the continuous alteration of the acoustic waveform during the propagation through a linear system. Based on the results of this analysis, the surface-related measurements are described as a convolution of the impulse response of the system with the non-stationary forward wavelet which includes all multiple terms generated within the system. The amplitude spectrum of the forward wavelet is recovered from the amplitude spectrum of the recorded signal, using the difference between the rate of decay of the source wavelet and the duration of the measurement. The phase spectrum of the forward wavelet is estimated using the Hilbert transform and the fact that the mixed phase lag wavelet can be presented as a convolution of the minimum and maximum phase lag wavelets. The multiples are discriminated from primaries by comparison of the phase spectrum of the seismic signal and the inverse of the forward wavelet. Therefore, the technique is called phase inversion deconvolution (PID). This approach requires no velocity information in order to recognize and attenuate multiple energy. Therefore, primary energy is recovered in the near-offset region where the velocity differential between primary and multiple energies is very small.     

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.     

美国中东部地震台阵波场的三分量波形梯度

波形梯度法是一种全新的台阵数据处理技术,该方法利用子台阵中的波形差异可以得到一些基本的地震波传播参数.本文首次将其应用于美国中东部地震台阵面波的三分量研究当中.首先通过垂向分量的面波进行波形梯度分析,得到垂直分量面波的相速度、传播方向、几何扩散和辐射花样.再通过垂向分量得到的传播方向进行坐标旋转,从而得到了径向分量以及切向分量的面波,再将其应用于波形梯度分析,分别得到径向分量以及切向分量面波的相速度、传播方向、几何扩散和辐射花样.利用2012年8月27日发生在中美洲西海岸地震事件的三分量地震数据的面波波形得到的结果显示,在相同周期,研究区域的三个分量面波的相速度分布横向差异显著.切向分量面波的相速度分布特征差异较大,可能是由于径向各向异性造成的.三个分量的传播方向变化都不大,且切向分量的传播方向变化大于垂直分量与径向分量,说明地震波的切线分量在传播过程中受到的影响更大,同时还可以看出传播方向的变化呈现出条带状的特征.几何扩散和辐射花样都是与地震波的振幅项有关的信息,三个分量的几何扩散特征基本一致.但是由于切向分量传播方向变化相对较大,可能导致了切向分量面波的辐射花样有所差异.     

基于F-K偏移和反偏移的地震道插值方法研究

野外采集的地震数据经常存在空道或者坏道的情况,为了满足地震数据处理精度的要求,就必须首先进行插值.本文提出了一种新的地震数据插值方法：F-K偏移和反偏移插值法,该方法是通过F-K偏移和反偏移的串联使用来实现的.与其他的插值方法相比,这种插值方法的优点在于计算速度快,没有经过近似处理,精度高.F-K偏移/反偏移能够实现道插值的原理和Kirchhoff偏移/反偏移插值方法类似,也是由于数据是有限带宽的原因引起的.通过对模型的试算,可以看到F-K偏移和反偏移插值方法有着比较好的插值效果,是一种高效、准确、可行的方法.     

Synthesis of a seismic virtual reflector*

We describe a method to process the seismic data generated by a plurality of sources and registered by an appropriate distribution of receivers, which provides new seismic signals as if in the position of the receivers (or sources) there was an ideal reflector, even if this reflector is not present there. The data provided by this method represent the signals of a virtual reflector. The proposed algorithm performs the convolution and the subsequent sum of the real traces without needing subsurface model information. The approach can be used in combination with seismic interferometry to separate wavefields and process the reflection events. The application is described with synthetic examples, including stationary phase analysis and with real data in which the virtual reflector signal can be appreciated.     

On the aperture effect in 3D Kirchhoff-type migration

It is well known that the migrated image given by a Kirchhoff-type (diffraction-stack) migration with limited aperture is always accompanied by some events which depend on the migration aperture. Although these events may severely affect the quality of migration, they have been studied only in 2D cases. Here, the events due to the migration aperture in 3D situations are investigated using a new method of analysing the reconstructed wavefield. It is found that a finite migration aperture results in a reconstructed wavefield with two components. One comes from the tangent points and curves between the traveltime surfaces of reflected and point-diffracted rays and is independent of the migration aperture, and the other is from the boundary of the migration aperture and depends strongly on the location and size as well as on the shape of the migration aperture. It is this last component that describes the aperture effect in migration. If the migration aperture is not sufficiently large, and if the input for migration is not zero on the boundary of the migration aperture, the boundary component may partially or totally cancel the migration signal. Furthermore, for synthetic data, the aperture effect cannot be eliminated by enlarging the migration aperture because, except for the common-shotpoint data, the aperture effect always exists however large the migration aperture becomes. This leads to the conclusion that the published Kirchhoff-type operators are not the exact inverse operators of the Fresnel–Kirchhoff integral if the input data are synthetic.     

THE DETERMINATION OF THE SEISMIC QUALITY FACTOR Q FROM VSP DATA: A COMPARISON OF DIFFERENT COMPUTATIONAL METHODS1

Ten methods for the computation of attenuation have been investigated, namely: amplitude decay, analytical signal, wavelet modelling, phase modelling, frequency modelling, rise-time, pulse amplitude, matching technique, spectral modelling and spectral ratio. In particular, we have studied the reliability of each of these methods in estimating correct values of Q using three synthetic VSP seismograms for plane P-waves with different noise contents. The investigations proved that no single method is generally superior. Rather, some methods are more suitable than others in specific situations depending on recording, noise or geology. The analytical signal method has been demonstrated to be superior if true amplitude recordings are available. Otherwise spectral modelling or, in the ‘ noise-free’ case the spectral ratio method, is optimal. Finally, two field VSPs in sediments are investigated. Only in the case of the highest quality VSP can significant information be deduced from the computed attenuation.     

Effects of wave-wave and wave-mean flow interactions on the evolution of a baroclinic wave

Abstract

The weakly nonlinear evolution of a free baroclinic wave in the presence of slightly supercritical, vertically sheared zonal flow and a forced stationary wave field that consists of a single zonal scale and an arbitrary number of meridional harmonics is examined within the context of the conventional two-layer model. The presence of the (planetary-scale) stationary wave introduces zonal variations in the supercriticality and is shown to alter the growth rate and asymptotic equilibrium of the (synoptic-scale) baroclinic wave via two distinct mechanisms: The first is due to the direct interaction of the stationary wave with the shorter synoptic wave (wave-wave mechanism), and the second is due to the interaction of the synoptic wave with that portion of the mean field that is corrected by the zonally rectified stationary wave fluxes (wave-mean mechanism). These mechanisms can oppose or augment each other depending on the amplitude and spatial structure of the stationary wave field. If the stationary wave field is confined primarily to the upper (lower) layer and consists of only the gravest cross-stream mode, conditions are favorable (unfavorable) for nonzero equilibrium of the free wave.

In addition to the time dependent heat flux generated by baroclinic growth of the free wave, its interaction with a stationary wave field consisting of two or more meridional harmonics generates time dependent heat fluxes that vary with period of the free wave. However, if the stationary wave field contains several meridional harmonics of sufficiently large amplitude, the free baroclinic wave is destroyed.     

小波变换的过零点特性与地震勘探信号的信噪比和分辨率

利用小波变换研究地震勘探信号小波变换的过零点特性，本文提出了用小波变换的过零点特性和地震勘探信号相邻道的横向相关性提高信号分辨率和信噪比的新方法．该方法包括两个主要步骤：①利用相邻地震道信号具有很好相关性，而噪音相关性差的特点以及小波变换的过零点特性得到有效反射波同相轴随空间坐标的变化信息．②利用奇异值分解和最小二乘（SVD-TLS）方法沿同相轴对振幅进行多项式拟合去噪并增加信号高频提高信号分辨率．     

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

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

Bayesian frequency-domain blind deconvolution of ground-penetrating radar data

Enhancing the resolution and accuracy of surface ground-penetrating radar (GPR) reflection data by inverse filtering to recover a zero-phased band-limited reflectivity image requires a deconvolution technique that takes the mixed-phase character of the embedded wavelet into account. In contrast, standard stochastic deconvolution techniques assume that the wavelet is minimum phase and, hence, often meet with limited success when applied to GPR data. We present a new general-purpose blind deconvolution algorithm for mixed-phase wavelet estimation and deconvolution that (1) uses the parametrization of a mixed-phase wavelet as the convolution of the wavelet's minimum-phase equivalent with a dispersive all-pass filter, (2) includes prior information about the wavelet to be estimated in a Bayesian framework, and (3) relies on the assumption of a sparse reflectivity. Solving the normal equations using the data autocorrelation function provides an inverse filter that optimally removes the minimum-phase equivalent of the wavelet from the data, which leaves traces with a balanced amplitude spectrum but distorted phase. To compensate for the remaining phase errors, we invert in the frequency domain for an all-pass filter thereby taking advantage of the fact that the action of the all-pass filter is exclusively contained in its phase spectrum. A key element of our algorithm and a novelty in blind deconvolution is the inclusion of prior information that allows resolving ambiguities in polarity and timing that cannot be resolved using the sparseness measure alone. We employ a global inversion approach for non-linear optimization to find the all-pass filter phase values for each signal frequency. We tested the robustness and reliability of our algorithm on synthetic data with different wavelets, 1-D reflectivity models of different complexity, varying levels of added noise, and different types of prior information. When applied to realistic synthetic 2-D data and 2-D field data, we obtain images with increased temporal resolution compared to the results of standard processing.