基于平面波模型重访地震背景噪声互相关及空间自相关(SPAC)

自Aki(1957)提出微震的空间自相关(SPatial AutoCorelation,SPAC)技术以来,SPAC技术一直独立发展,并在工程地震领域获得了广泛应用.近20年来,地震干涉(Seisimic Interferometry,SI)在多个领域引起人们的关注,该技术的核心思想是连续地震背景噪声的互相关函数(Noise Crosscorrelation Function,NCF)可以重建系统的格林函数(Green's Function,GF),对该技术的回溯性研究建立了SPAC和NCF的关系:它们是对同一物理现象的不同描述,SPAC在频率域中描述随机平稳噪声的空间相干,NCF在时间域中描述扩散场的互相关.理论上SAPC和NCF技术要求背景噪声源均匀分布,这样的噪声场可以用平面波叠加来模拟.本文基于平面波模型重访地震背景噪声的互相关和空间自相关技术,从单色平面波的互相关表示出发,对地震背景噪声互相关及空间自相关技术进行评述,试图使这些概念更易于理解.与之前众多研究地震干涉技术的理论相比,本文尤其关注以下几点:(1)基于简单的平面波模型,给出不同维度下,源或台站对方位均匀分布时,平面波互相关对入射波的方位平均和台阵对的方位平均结果,并对格林函数GF和时域互相关函数NCF的关系进行总结.(2)给出声源和(或)交叉台站方位分布不均匀时的互相关表示,指出这种非均匀性对方位的依赖关系,与弱各向异性介质中面波速度的方位依赖关系类似,因此,非均匀源的影响在反演时可能会映射到面波方位各向异性结果中.(3)互相关运算中,哪一个台站是虚拟源.NCF包含因果性和非因果性两部分,NCF的非对称性通常用于研究噪声源的方位分布,但由于源和接收的互易关系,及对互相关运算的不同定义和不同的傅里叶变换习惯,哪一个台站是虚拟源在目前的文献中并不明确.(4)方位平均和时间平均的关系.在SPAC处理中,需要对不同方位分布的台站对进行方位平均,本文从理论上说明,单个平面波入射时,交叉台站互相关系数对台站对的方位平均,等价于单个台站对互相关系数对入射波的时间平均.(5)几种特定分布非均匀噪声源的SPAC表示.包括单独的因果性噪声源和非因果性噪声源给出的互相关函数表示,及由此带来的相移问题.(6)利用SPAC、NCF和面波GF之间的关系,给出交叉分量的空间自相关系数表示.(7)衰减介质的空间相干表示.虽然利用地震干涉技术研究介质衰减在理论上仍然存在一些争议,但人们正试图研究从连续背景噪声记录中提取介质衰减的可能性.本文基于平面波模型,给出了不同坐标选择下,衰减介质的空间相干表示,这种表达的不同,指示了由地震干涉技术提取介质衰减的困难.与众多研究地震干涉的理论相比,比如稳相近似理论、互易定理、时间反转声学等,本文主要考虑均匀介质,不涉及非均匀介质的散射,从最简单的平面波模型,理解背景噪声重建系统格林函数这一地震干涉的核心思想和相应的基本概念.     

Use of the Fourier-Laplace transform and of diagrammatical methods to interpret pumping tests in heterogeneous reservoirs

Advances in computer power and in reservoir characterization allow simulation of pressure transients in complex reservoirs generated stochastically. Generally, interpretation of these transient gives useful information about the reservoir hydraulic properties: a major goal is to interpret these transients in the stochastic context. First we ensemble average the pressure over all the random permeability field realizations to derive an equation which drives the ensemble averaged pressure. We use the Fourier transform in space and the Laplace transform in time, in conjuction with a perturbation series expansion in successive powers of the permeability fluctuations to obtain an explicit solution. The Nth order term of this series involves the hydrodynamic interaction between N permeability heterogeneities and after averaging we obtain an expansion containing correlation functions of permeability fluctuations of increasing order.Next, Feynman graphs are introduced allowing a more attractive graphical interpretation of the perturbation series. Then series summation techniques are employed to reduce the graph number to be summed at each order of the fluctuation expansion. This in turn gives useful physical insights on the homogenization processes involved. In particular, it is shown that the sum of the so-called ‘one-particle irreducible graphs’ gives the kernel of a linear integro-differential equation obeyed by the ensemble average pressure. All the information about the heterogeneity structure is contained in this renormalized kernel, which is a limited range function.This equation on its own is the starting point of useful asymptotic results and approximations. In particular it is shown that interpretation of pumping tests yields the steady-state equivalent permeability after a sufficiently long time for an infinite reservoir, as expected.     

轨道附近地面振动模型中的饱和地层动力格林函数

列车引起场地振动的建模需要能够表达地层的动力格林函数.本文兼顾饱和土的流固两相耦合性、场地土的分层性和波动的三维传播性,构建了半解析的场地动力格林函数.首先,基于Biot方程,在傅里叶变换域求解固体骨架和流体的位移和应力.然后采用传递矩阵方法建立地表位移和应力间的关系,得到格林函数矩阵.进而讨论矩阵的一些固有特征,提出改善竖向位移计算效率的措施.最后利用推导的格林函数计算了几个典型算例.数值结果与文献中其他方法得到的结果十分接近,与场地振动的现场观测试验基本符合.软土场地振动的计算结果高于饱和砂土场地,高速列车场地振动强度高于低速列车.当车速接近场地瑞利波速,模拟结果中显示出马赫锥.数值结果还显示,即使车速略低于瑞利波速,马赫锥也可能出现.本文推导的格林函数将有助于深入理解列车等移动激励作用下层状饱和土场地的振动特征.     

利用广义叠加反演方法估算长岛震群应力降

地震应力降是表征震源特性的一个重要参数,通常由震源谱计算得到.但如何从大量地震的观测谱中校正路径和台站效应,准确地分离出震源谱并计算应力降,仍是一个难题.本文采用广义叠加反演方法通过分阶段叠加从观测谱中迭代分离出震源项、台站项和路径项,并利用分震级的叠加震源项与理论谱的拟合得到经验格林函数,再利用由经验格林函数校正后的真实震源谱估算出应力降.该方法不需要知道仪器响应、场地响应和传播路径上的衰减系数,并且不用事先假设震源具有自相似的恒应力降特征,计算效率很高,适用于反演大量地震的震源参数.我们将广义叠加反演方法应用于长岛震群,利用S波谱估算了1431个ML≥1.0地震的应力降.结果表明,广义叠加反演方法获得的震源参数和用传统方法得到的结果一致;长岛震群地震的应力降较小(中值是0.3 MPa),变化范围很大(0.009～3.04MPa之间),属于低应力降事件集;同时应力降中值随地震矩增大而增大,偏离了自相似理论.应力降空间变化不均匀,高应力降事件主要集中在震群的端部(NW端),到最近断层距离的依赖性不强,在9、10km处的应力降中值略高于其他深度.根据应力降的变化特征和区域构造条件,我们推测地下流体在长岛震群的活动中起到一定的作用.     

Stratigraphic filtering and source penetration depth

Seismic exploration underneath highly heterogeneous layers such as basalt flows is possible by lowering the principal source frequency. Unfortunately this also reduces resolution. Wave‐localization theory is a multiple scattering theory that can be used to study stratigraphic filtering in chaotic lithologies. It predicts the apparent attenuation due to scattering of a plane wave traversing a layer with high velocity fluctuations. It can therefore predict the optimum principle source frequency in the trade‐off between loss of resolution and increased penetration depth. We show how this can be done with the help of a few statistical parameters derived from a well‐log analysis; namely, the average background velocity, the expected standard deviation in the velocity fluctuations, the typical scale length of the heterogeneities and the thickness of the basalt layer. In the likely situation that no local well logs exist, a multitude of scenarios can easily be examined at low cost.     

格林函数的奇异性处理

格林函数在解决电磁场和直流电场问题的积分方程法和边界单元法等方法中有着广泛的应用.在应用格林函数时会遇到格林函数奇异性的问题,即其在源区域的内部,当场点和源点重合时,会碰到对格林函数为无穷大,积分为奇异积分.对于不同问题中遇到的不同形式的格林函数奇异性问题的各种处理方法进行了分析和评述.对于标量格林函数奇异性问题处理方法有：挖去法、级数展开法、绕开法和解析法等.对于并矢格林函数奇异性问题处理方法有：分量处理法、源并矢法和拟源并矢法等处理方法.通过对实际工作中所遇到的新方法新问题中的不同形式的格林函数的研究,提供了奇异性问题的处理的方法和途径.     

A Comparison of Strategies for Seismic Interferometry

The extraction of the response from field fluctuations excited by random sources has received considerable attention in a variety of different fields. We present three methods for the extraction of the systems response that are based on cross-correlation, deconvolution, and the solution of an integral equation, respectively. For systems that are invariant for time-reversal the correlation method requires random sources on a bounding surface only, but when time-reversal invariance is broken, for example by attenuation, a volume distribution of sources is needed. For this reason the correlation method is not useful for diffusive or strongly attenuating systems. We provide examples of the three methods and compare their merits and drawbacks. We show that the extracted field may satisfy different boundary conditions than does the physical field. This can be used, for example, to suppress surface-related multiples in exploration seismology, to study the coupling of buildings to the subsurface, and to remove the airwave in controlled source electromagnetics (CSEM).     

Boundary integral equations for dynamic rupture propagation on vertical complex fault system in half-space:Theory

The boundary integral equation method(BIEM)is now widely used in numerical studies on earthquake rupture dynamics,and is proved to be a powerful tool to deal with problems on complex fault system.However,since this method heavily lies on the specific forms of Green's function and only the Green's function in full-space has a closed analytic expression,it is usually limited to a full-space medium.In this study,as a first step to extend this method to an arbitrary complex fault system in half-space,the boundary integral equations(BIEs)for dynamic strike-slip on vertical complex fault system in half-space are derived based on exact Green's function for isotropic and homogeneous half-space.Effect of the geometry of the complex fault system are dealt with carefully.Final BIEs is composed of two parts:contribution from full-space,which has been thoroughly investigated by Aochi and his co-workers by using the Green's function for full-space,and that from free surface,which is studied in detail in this study.     

Interstation surface wave attenuation by autoregressive deconvolution

A new technique for calculating interstation Green's functions and attenuation coefficients for seismic surface waves is presented.The interstation Green's function is evaluated from the autocorrelation functions of the seismograms, which are obtained from a maximum entropy process.Since a data-invariant time window is not used, the evaluated Green's functions gives reliable information on both the amplitude and the phase spectra of the system.This new technique is compared with other methods by applying them to both synthetic and real data from a path in the Canadian shield.     

Consistency of the spatial autocorrelation method with seismic interferometry and its consequence

We have cross‐checked the conventional theory of the spatial autocorrelation method and the consequence of seismic interferometry: the retrieval of the elastodynamic Green's function. Their mutual consistency is almost complete. The basic formulas of the conventional spatial autocorrelation theory can be derived by an alternative approach based on the retrieval of the elastodynamic Green's function. The only discrepancy is found with the average of the complex coherence function over azimuth in a wavefield dependent on azimuth. It is hypothesized, in discussion, that this discrepancy is due to the way of representing the wavefield in the background theory of seismic interferometry that can produce only wavefields moderately dependent on azimuth and that the mentioned consequence of seismic interferometry can also only make sense in a wavefield moderately dependent on azimuth. Our field experiment with a wavefield dependent on azimuth showed that the consequence of seismic interferometry in the logical framework of the conventional spatial autocorrelation theory is appropriate under such degrees of approximation as the measure proposed in this study, i.e., the deviation of the total dispersion curves is between about 10 and 16 per cent at the maximum from those averaged over azimuth. The acceptance of the retrieval of Green's function gives a proper physical meaning to the complex coherence function: the real part of the elastodynamic Green's function normalized by its zero‐offset version. This makes it possible to take a deterministic approach rather than the statistical one on which the conventional spatial autocorrelation method is based and gives fruitful new aspects and perspectives. For example, the formula for the multi‐mode case is given and the possibility of exploration of two or three dimensional velocity structures is suggested.     

Vector-scattering of elastic waves by directional structural space gradients

A nonstochastic and noniterative theory of vector scattering in inhomogeneous media is presented. The elastodynamic vector wave-equation for 3D inhomogeneous media is solved for a weak heterogeneity at the high-frequency region. It is shown that there exists a forward scattered field which decays slowly along the source-receiver path. Its rate of attenuation depends on the azimuth of the path relative to the direction of the inhomogeneity, but is independent of frequency. The Green's tensor for the above regime is derived in closed form and leads to the quantification of fields of dipolar sources in weak inhomogeneous media. The inhomogeneity at the source creates a source-induced scattering (in addition to path-scattering) having a radiation-pattern that bears the signature of the source. The availability of the analytic Green's tensor, in conjunction with the Huygens-Kirchhoff-Helmholtz formalism, opens new ways to calculate the scattered fields due to various structural inhomogeneities applicable to exploration and earthquake seismology. The theoretical results of this study point to the conclusion that the scalar wave approximation may not always be valid for the propagation of seismic waves in the earth's lithosphere.     

Decoupled elastic prestack depth migration

This paper presents a new decoupled form of the formula for common-shot or common-receiver amplitude-preserving elastic prestack depth migration (PreSDM), which can be used for estimating angle-dependent elastic reflection coefficients in laterally inhomogeneous anisotropic media. The multi-shot or multi-receiver extension of this formula is suitable for automated prestack amplitude-versus-angle (AVA) elastic inversion of ocean-bottom cable (OBC), walkaway VSP (WVSP) or standard towed-cable data at any subsurface location. The essence of the theory is a systematic application of the stationary-phase principle and high-frequency approximations to the basic elastic Green's theorem. This leads to nonheuristic explicit wave mode decoupling and scalarization of vector PreSDM. Used in combination, ray-trace and finite-difference (FD) eikonal solvers create a useful tool to calculate accurate Green's function travel time and amplitude maps. Examples of synthetic OBC data and applications to field WVSP data show that the new imaging technique can produce a clear multi-mode elastic image.     

基于delta波包叠加的时间域深度偏移

delta波包可由高斯波束经傅里叶逆变换得到,是高斯波束在时空域的表达.它最早出现在合成理论地震图的研究中,本文将其应用于偏移领域.通过delta波包叠加表达时间域格林函数,可将高斯波束偏移由频率域转换到时间域,再结合Rayleigh积分和激励时间成像条件,本文给出了基于delta波包叠加的深度偏移算法.该偏移算法可在时间域直接计算,但因包含褶积运算,成像时将耗费大量的计算时间.针对这一问题,本文提出了将褶积简化为乘积的近似公式.近似后的偏移算法,不仅保留了高斯波束偏移的优点,而且计算效率得到显著提升.文中通过两个数值算例验证了上述结论.     

利用互相关函数求地震波衰减

与传统利用单事件波形求介质衰减的方法不同,我们提出一种以一系列事件为源,以源与远处记录波形的互相关函数为基础计算介质品质因子的方法.依据野外现场实验的数据,我们分别以单次和多次激发为源,讨论了事件波形和互相关函数波形的振幅变化及相应的介质品质因子,结果发现,与事件波形一样,互相关函数波形也可以用于计算介质的衰减系数.计算结果表明：（1）实验场地浅层介质在原场条件下的品质因子Q值在40左右;（2）互相关函数的最大幅度与事件波形的振幅随距离具有相同的衰减规律.     

Weakly nonlinear analysis of the α effect

We address mathematical issues raised by the so-called α?effect of dynamo theory, which is a dynamo mechanism arising in conducting flows with small scale fluctuations. Analytical results on the α?effect concern the linear induction equation, and are usually claimed to hold for the whole magnetohydrodynamics (MHD) system, as long as the amplitude of the perturbations is small. We discuss the justification of that claim, in the case of periodic fluctuations of the fields. We show a nonlinear instability result on the MHD system, that predicts dynamo action for a large class of high frequency periodic flows, up to the fully nonlinear regime.     

Constant Q and a Fractal, Stratified Earth

— Frequency-dependent measurements of the quality factor Q typically show a constant behaviour for low frequencies and a positive power-law dependence for higher frequencies. In particular, the constant Q pattern is usually explained using intrinsic attenuation models due to anelasticity with either a single or multiple superposed relaxation mechanisms — each with a particular resonance peak.¶However, in this study, I show using wave localisation theory that a constant Q may also be due to apparent attenuation due to scattering losses. Namely, this phenomenon occurs if the earth displays fractal characteristics. Moreover, if fractal characteristics exist over a limited range of scales only, even an absorption band can be created—in accordance with observations. This indicates that it may be very difficult to distinguish between intrinsic and scattering attenuation on the basis of frequency-dependent measurements of the quality factor only.     

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

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

Stationary phase approximation in the ambient noise method revisited

The method of extracting Green's function between stations from cross correlation has proven to be effective theoretically and experimentally.It has been widely applied to surface wave tomography of the crust and upmost mantle.However,there are still controversies about why this method works.Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space,and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green's function.His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green's function.However,his derivation ignores influence from the two stationary points at infinities,therefore it may fail when attenuation is absent.In order to obtain accurate noise-correlation function due to scatters over the whole space,we compute the total contribution with numerical integration in polar coordinates.Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function,in both amplitude and phase,when the frequency is low with reasonable quality factor Q.Our results argue that the distance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method,and only the station pairs whose distances are above several (5) wavelengths can be used.