海上多次波的联合衰减法

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

STABILIZATION OF NORMAL-INCIDENCE SEISMOGRAM INVERSION REMOVING THE NOISE-INDUCED BIAS*

A main problem in computing reflection coefficients from seismograms is the instability of the inversion procedure due to noise. This problem is attacked for two well-known inversion schemes for normal-incidence reflection seismograms. The crustal model consists of a stack of elastic, laterally homogeneous layers between two elastic half-spaces. The first method, which directly computes the reflection coefficients from the seismogram is called “Dynamic Deconvolution”. The second method, here called “Inversion Filtering”, is a two-stage procedure. The first stage is the construction of a causal filter by factorization of the spectral function via Levinson-recursion. Filtering the seismogram is the second stage. The filtered seismogram is a good approximation for the reflection coefficients sequence (unless the coefficients are too large). In the non-linear terms of dynamic deconvolution and Levinson-recursion the noise could play havoc with the computation. In order to stabilize the algorithms, the bias of these terms is estimated and removed. Additionally incorporated is a statistical test for the reflection coefficients in dynamic deconvolution and the partial correlation coefficients in Levinson-recursion, which are set to zero if they are not significantly different from noise. The result of stabilization is demonstrated on synthetic seismograms. For unit spike source pulse and white noise, dynamic deconvolution outperforms inversion filtering due to its exact nature and lesser computational burden. On the other hand, especially in the more realistic bandlimited case, inversion filtering has the great advantage that the second stage acts linearly on the seismogram, which allows the calculation of the effect of the inversion procedure on the wavelet shape and the noise spectrum.     

来自海底高速层径向波的理论地震图研究

本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波（简称径向波）,这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义.     

一种改进的基于非高斯性最大化的预测反褶积算法

The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.     

NON-LINEAR ESTIMATION OF REFLECTION COEFFICIENTS FROM SEISMIC DATA*

For years, reflection coefficients have been the main aim of traditional deconvolution methods for their significant informational content. A method to estimate seismic reflection coefficients has been derived by searching for their amplitude and their time positions without any other limitating assumption. The input data have to satisfy certain quality constraints like amplitude and almost zero phase noise—ghosts, reverberations, long period multiples, and diffracted waves should be rejected by traditional processing. The proposed algorithm minimizes a functional of the difference between the spectra of trace and reflectivity in the frequency domain. The estimation of reflection coefficients together with the consistent “wavelet’ is reached iteratively with a multidimensional Newton-Raphson technique. The residual error trace shows the behavior of the process. Several advantages are then obtainable from these reflection coefficients, like conversion to interval velocities with an optimum calibration either to the well logs or to the velocity analysis curves. The procedure can be applied for detailed stratigraphic interpretations or to improve the resolution of a conventional velocity analysis.     

Combining two seismic experiments to attenuate free-surface multiples in OBC data*

The current inverse scattering solution used for multiple attenuation of marine seismic reflection data assumes that sources and receivers are located in the water. To adapt this solution to the ocean-bottom cable (OBC) experiment where receivers are located on the sea-floor, we have proposed combining the conventional marine surface seismic reflection data (streamer data) with OBC data. The streamer data add to the OBC data some of the wave paths needed for multiple attenuation. This combination has allowed us to develop a multiple attenuation method for OBC data which does not require any knowledge of the subsurface and which takes into account all free-surface multiples, including receiver ghosts. A non-linear synthetic data example consisting of pressure and particle velocity fields is used to illustrate the procedure.     

一阶多次波聚焦变换成像

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

MULTIPLE REFLECTIONS AND HEAD WAVES IN THE GULF OF SUEZ*

Inferior reflection quality in the Gulf of Suez at the target depth interval is attributable in part to surficial multiple reflections. An excellent example of the latter is observed on a typical seismic line in the northern portion of the Gulf. An increase in prominence of the multiple reflections appears associated with decreasing depth to a dipping highvelocity layer. Inversion of a second-order polynomial time-distance function, fitted to the observed refraction onset time-distance values, gives the velocity-depth function for sediments between the water bottom and a high-velocity layer. Velocities thus determined increase non-linearly with depth from a value near water velocity at the water bottom. Depths to the high-velocity layer are obtained from the associated head-wave linear time-distance function and by ray tracing in the overlying sediments. As the high-velocity layer approaches the water bottom from sub-water depths exceeding 0.6 km to a depth of 56 m, intensity of the multiple reflections increases to the extent of completely dominating individual records to a time of at least 3 s. The estimated plane-wave normal-incident reflection coefficient at the top of the high-velocity layer increases with decreasing depth to this layer, approaching 0.5 at the shallowest depth. This strong reflection coefficient further substantiates the existence of multiple reflections between the high-velocity layer and water layer. However, existence of water-layer multiples cannot be ruled out. The estimated water-bottom reflection coefficient is approximately 0.3, a substantial value. Multiple reflections of considerably less intensity are apparent where the high-velocity layer is deepest, and it is likely that such are waterlayer multiple reflections. Unfortunately, water-layer multiple reflections and multiple reflections between the water surface and high-velocity layer cannot be separated by their coincidence with time-distance (normal moveout) curves, the configuration of each visibly matching the curves equally well.     

地震资料叠前去噪技术的现状与未来

地震资料叠前去噪是勘探地震资料处理的关键问题之一，但这个问题长期以来一直没有得到很好的解决，是提高地震资料分辨率的一个主要障碍，其中，如何有效地消除产生与地层间的多次反射波又是地震资料叠前去噪的核心问题。为此，人们进行了长期不懈的努力，以更好地消除多次波。本文着重概述现今实生常中常用的几种比较有效的消除多次波方法，包括拉冬变换，F－K方法和聚束滤波方法；同时也概述了预测反积消除水层混响和用K－L变换方法消除随机噪音提高信噪比等方法；最后，还讨论了近年发展的地震资料叠前去噪新方法及其发展趋势。     

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

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

基于广义高斯分布的地震盲反褶积方法研究

Bussgang算法是针对褶积盲源分离问题提出的,本文将其用于地震盲反褶积处理.由于广义高斯概率密度函数具有逼近任意概率密度函数的能力,从反射系数序列的统计特征出发,引入广义高斯分布来体现反射系数序列超高斯分布特征.依据反射系数序列的统计特征和Bussgang算法原理,建立以Kullback-Leibler距离为非高斯性度量的目标函数,并导出算法中涉及到的无记忆非线性函数,最终实现了地震盲反褶积.模型试算和实际资料处理结果表明,该方法能较好地适应非最小相位系统,能够同时实现地震子波和反射系数估计,有效地提高地震资料分辨率.     

ADAPTIVE PREDICTIVE DECONVOLUTION OF SEISMIC DATA*

The Wiener prediction filter has been an effective tool for accomplishing dereverberation when the input data are stationary. For non-stationary data, however, the performance of the Wiener filter is often unsatisfactory. This is not surprising since it is derived under the stationarity assumption. Dereverberation of nonstationary seismic data is here accomplished with a difference equation model having time-varying coefficients. These time-varying coefficients are in turn expanded in terms of orthogonal functions. The kernels of these orthogonal functions are then determined according to the adaptive algorithm of Nagumo and Noda. It is demonstrated that the present adaptive predictive deconvolution method, which combines the time-varying difference equation model with the adaptive method of Nagumo and Noda, is a powerful tool for removing both the long- and short-period reverberations. Several examples using both synthetic and field data illustrate the application of adaptive predictive deconvolution. The results of applying the Wiener prediction filter and the adaptive predictive deconvolution on nonstationary data indicate that the adaptive method is much more effective in removing multiples. Furthermore, the criteria for selecting various input parameters are discussed. It has been found that the output trace from the adaptive predictive deconvolution is rather sensitive to some input parameters, and that the prediction distance is by far the most influential parameter.     

Optimal Estimation of Two-Dimensional Contaminant Transport

The state-space estimation technique presented herein provides a method for obtaining optimal estimates of concentrations for two-dimensional plumes in ground water. The concentration of a plume was defined as the state variable. The technique uses the Kalman filter and involves combining two independent estimates of plume concentrations. One estimate is called the process modeling data and the other is called the “measurement” data. The process modeling data is obtained from a numerical model. The “measurement” data is obtained from field measurements; however, for illustration in this paper it was generated by a different numerical model than the one used to obtain the process modeling data. The state-space technique produces a distribution of contaminant concentrations that is more accurate than either of the distributions generated by the process modeling or the “measurement” data. An example is presented to show that the technique produces significant improvements in the prediction of plume concentration distributions.     

OBC地震物理模型实验中的广角反射现象(英文)

广角地震采集作为获取野外高质量地震数据的一种重要方法,愈来愈受到人们的关注。因此,对广角反射的深入研究具有重要实际意义。本文通过层状物理模型海底电缆(0BC)地震勘探实验研究,比较详细地讨论了广角反射现象。一些实验结果并不支持理论数值模拟结论。主要实验结论有:1)广角反射波的振幅能量较非广角反射波强(约为1倍多),但和理论计算的15倍关系差异甚大,过临界角时反射能量缓慢增大,而不是理论计算的那样急剧增大;2)大炮检距时,反射同相轴仍保持近似双曲线趋势;3)广角反射波的主频较非广角反射波低(降低20-30%),且偏移距增大有所降低,而非广角反射在小偏移距段主频变化不明显;4)在临界角前后反射波波形没有发生突变,极性无变化;5)反射波组特征在临界角前后也未发生变化。6)水中直达波、多次波和水底折射波对海底层广角反射有一定的影响。     

3D targeted multiple attenuation

Short-period multiple reflections pose a particular problem in the North Sea where predictive deconvolution is often only partially successful. The targeted multiple attenuation (TMA) algorithm comprises computation of the covariance matrix of preflattened prestack or post-stack seismic data, the determination of the dominating eigenvectors of the covariance matrix, and subtraction of the related eigenimages followed by reverse flattening. The main assumption made is that the flattened multiple reflections may be represented by the first eigenimage(s) which implies that the spatial amplitude variations of primaries and associated multiples are similar. This assumption usually limits the method to short-period multiple reflections. TMA is applicable post-stack or prestack to common-offset gathers. It is computationally fast, robust towards random noise, irregular geometry and spatial aliasing, and it preserves the amplitudes of primaries provided they are not parallel to the targeted multiples. Application of TMA to 3D wavefields is preferable because this allows a better discrimination between primaries and multiples. Real data examples show that the danger of partially removing primary energy can be reduced by improving the raw multiple model that is based on eigenimages, for example by prediction filtering.     

OFFSET CONTINUATION OF SEISMIC SECTIONS*

Geometrical acoustic and wave theory lead to a second-order partial differential equation that links seismic sections with different offsets. In this equation a time-shift term appears that corresponds to normal moveout; a second term, dependent on offset and time only, corrects the moveout of dipping events. The zero-offset stacked section can thus be obtained by continuing the section with maximum offset towards zero, and stacking along the way the other common-offset sections. Without the correction for dip moveout, the spatial resolution of the section is noticeably impaired, thus limiting the advantages that could be obtained with expensive migration procedures. Trade-offs exist between multiplicity of coverage, spatial resolution, and signal-to-noise; in some cases the spatial resolution on the surface can be doubled and the aliasing noise averaged out. Velocity analyses carried out on data continued to zero offset show a better resolution and improved discrimination against multiples. For instance, sea-floor multiples always appear at water velocity, so that their removal is simplified. This offset continuation can be carried out either in the time-space domain or in the time-wave number domain. The methods are applied both to synthetic and real data.     

Marine seismic wavefield measurement to remove sea-surface multiples

We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated. Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain.     

大容量气枪震源陆地反射地震探测——以长江中下游铜陵地区为例

为了实验大容量气枪震源陆地水体流动激发反射地震探测效果,在长江中下游安徽省铜陵段,采用气枪船长江航道流动激发、沿江岸布设反射地震仪器接收的非纵弯线工作方式,得到了反映测线经过地区地壳深部结构和构造特征的反射地震数据。原始资料信噪比较低,但部分资料不同部位仍可辨认出来自地壳及莫霍面反射波组。就传播距离而言,地震波传播的水平距离最大可达21km,垂直深度可达30km以上。在数据处理中,根据原始资料特点,针对性采用了非纵弯线面元定义、三维层析静校正、叠前多域去噪及组合反褶积技术,最终得到的叠加时间剖面上具有丰富的壳内反射波组。结果显示,测线经过地区的地壳结构为双层结构,总厚度为30.0～36.0km。上地壳呈现隆坳相间的反射特征,下地壳存在多组叠层状弧型反射波组,莫霍面反射特征清晰,由2～3个反射同相轴组成,呈现SW端向NE段抬升的形态。剖面经过地区存在一个切穿下地壳和莫霍面的深部断裂,应该是长江深断裂的反映。研究结果充分说明,大容量气枪震源可应用于陆地流动水体地壳精细结构的深地震反射探测。     

A framework for sourcing of evaporation between saturated and unsaturated zone in bare soil condition

ABSTRACT

Sourcing subsurface evaporation (E_ss) into groundwater (E_g) and unsaturated zone (E_u) components has received little scientific attention so far, despite its importance in water management and agriculture. We propose a novel sourcing framework, with its implementation in dedicated post-processing software called SOURCE (used along with the HYDRUS1D model), to study evaporation sourcing dynamics, define quantitatively “shallow” and “deep” water table conditions and test the applicability of water table fluctuation (WTF) and “bucket” methods for estimation of E_g and E_u separately.

For the “shallow” and “deep” water table we propose E_g?>?0.95E_ss and E_g = 0 criteria, respectively. Assessment of the WTF method allowed sourcing of very small fluxes otherwise neglected by standard hydrological methods. Sourcing with SOURCE software was more accurate than the standard “bucket” method mainly because of greater flexibility in spatio-temporal discretization. This study emphasized the dry condition relevance of groundwater evaporation which should be analysed by applying coupled flow of heat, vapour and liquid water.

Editor D. Koutsoyiannis; Associate editor S. Kanae     

分形脉冲反褶积方法

解地震反演问题的脉冲反褶积方法是基于反射系数白噪和子波为最小相位的假设下提出的．近几年的研究证明反射系数并不都是白噪,而是某种分形噪声,如果用一类分形反褶积方法,则将地震反演问题化为难以求解的非线性方程组．本文用反射系数的分形性质,推导出一个更为简单易解的线性方程组,称为分形脉冲反褶积．数值计算表明,本文的方法是有效的．