平面波在粘滞性界面上的反射特征研究

当介质中存在按照一定方向排列的裂隙时,会在该区域形成接触程度不好的界面.当有地震波传播时,界面两侧的质点不再满足位移连续的条件,位移差与该点的应力有关.从而引起质点间的滑动,形成粘滞性界面.界面的这种特性可以用切向柔量和法向柔量来表示,柔量值越大,界面的接触程度和质点位移的连续性都越差.文中详细推导了平面波在粘滞性界面上的反射和透射系数公式,通过数值算例证明了粘滞性界面上的反射系数与入射波频率有关的性质.并说明了界面的粘滞性会在一定程度上影响到反射波的AVO特性.且当平面P波入射时,法向柔量对反射P波和转换波的影响均大于切向柔量的影响.     

Seismic modelling for geological fractures

Based on knowledge of a commutative group calculation of the rock stiffness and on some geophysical assumptions, the simplest fractured medium may be regarded as a fracture embedded in an isotropic background medium, and the fracture interface can be simulated as a linear slip interface that satisfies non‐welded contact boundary conditions: the kinematic displacements are discontinuous across the interface, whereas the dynamic stresses are continuous across the interface. The finite‐difference method with boundary conditions explicitly imposed is advantageous for modelling wave propagation in fractured discontinuous media that are described by the elastic equation of motion and non‐welded contact boundary conditions. In this paper, finite‐difference schemes for horizontally, vertically, and orthogonally fractured media are derived when the fracture interfaces are aligned with the boundaries of the finite‐difference grid. The new finite‐difference schemes explicitly have an additional part that is different from the conventional second‐order finite‐difference scheme and that directly describes the contributions of the fracture to the wave equation of motion in the fractured medium. The numerical seismograms presented, to first order, show that the new finite‐difference scheme is accurate and stable and agrees well with the results of previously published finite‐difference schemes (the Coates and Schoenberg method). The results of the new finite‐difference schemes show how the amplitude of the reflection produced by the fracture varies with the fracture compliances. Later, comparisons with the reflection coefficients indicate that the reflection coefficients of the fracture are frequency dependent, whereas the reflection coefficients of the impedance contrast interface are frequency independent. In addition, the numerical seismograms show that the reflections of the fractured medium are equal to the reflections of the background medium plus the reflections of the fracture in the elastic fractured medium.     

平面声波在粗糙界面上的反射特征研究

基于有关粗糙界面的Rayleigh假设,讨论了平面声波按余弦规律快速变化的小尺度粗糙界面上的反射特征.研究表明：这类界面与位于该位置的一个过渡地层的作用相当.该过渡层的厚度为粗糙界面的起伏幅度,速度和密度为上下两层介质相应量的平均值.研究了埋藏很深的微粗糙界面所引起的地震反射（绕射）波的频散特性和走时的构成,即包含零炮检距反射时间、正常时差和界面粗糙时差三部分内容.该粗糙时差与空间坐标和时间坐标无关,与绕射波的阶次有关,绕射波尾随在反射波之后以某一固定的时差出现.且只有当界面的粗糙波长与地震波的波长相当时,才能观测到这类绕射波.该结论为粗糙界面地震反射资料的处理方法提供了理论依据.     

Calculating complex‐valued P‐wave first‐arrival traveltimes in attenuative vertical transversely isotropic media with an irregular surface

The complex‐valued first‐arrival traveltime can be used to describe the properties of both velocity and attenuation as seismic waves propagate in attenuative elastic media. The real part of the complex‐valued traveltime corresponds to phase arrival and the imaginary part is associated with the amplitude decay due to energy absorption. The eikonal equation for attenuative vertical transversely isotropic media discretized with rectangular grids has been proven effective and precise to calculate the complex‐valued traveltime, but less accurate and efficient for irregular models. By using the perturbation method, the complex‐valued eikonal equation can be decomposed into two real‐valued equations, namely the zeroth‐ and first‐order traveltime governing equations. Here, we first present the topography‐dependent zeroth‐ and first‐order governing equations for attenuative VTI media, which are obtained by using the coordinate transformation from the Cartesian coordinates to the curvilinear coordinates. Then, we apply the Lax–Friedrichs sweeping method for solving the topography‐dependent traveltime governing equations in order to approximate the viscosity solutions, namely the real and imaginary parts of the complex‐valued traveltime. Several numerical tests demonstrate that the proposed scheme is efficient and accurate in calculating the complex‐valued P‐wave first‐arrival traveltime in attenuative VTI media with an irregular surface.     

Analytical wavefront curvature correction to plane‐wave reflection coefficients for a weak‐contrast interface

Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients.     

基于层状双孔介质的地震波反射和透射系数频散特性研究

本文研究了纵波垂直入射情况下两种介质分界面处的纵波反射和透射系数的频散特性,分界面上下两侧分别为层状双孔页岩介质和层状双孔砂岩介质.当纵波沿垂直于分界面的方向传播至分界面处时,会在上层双孔介质中产生三类反射纵波,在下层双孔介质中产生三类透射纵波.基于层状双孔介质的特性,给出了分界面处的六个边界条件.根据层状双孔介质的波动方程,利用平面波分析得到了纵波的反射和透射系数.结果表明：当多孔介质中存在流体时,纵波的反射和透射系数与频率相关,即存在频散现象.波致流体流动是造成纵波反射和透射系数频散的主要原因.此外,结果还表明局部流体流动引起地震频带内反射和透射系数的频散,宏观Biot流引起超声频带内反射和透射系数的频散.本文同时对岩石参数对反射和透射系数频散曲线的影响进行了研究.     

Joint PP‐ and PSV‐wave amplitudes versus offset and azimuth inversion for fracture compliances in horizontal transversely isotropic media

We propose a robust approach for the joint inversion of PP‐ and PSV‐wave angle gathers along different azimuths for the elastic properties of the homogeneous isotropic host rock and excess compliances due to the presence of fractures. Motivated by the expression of fluid content indicator in fractured reservoirs and the sensitivity of Lamé impedances to fluid type, we derive PP‐ and PSV‐wave reflection coefficients in terms of Lamé impedances, density, and fracture compliances for an interface separating two horizontal transversely isotropic media. Following a Bayesian framework, we construct an objective function that includes initial models. We employ the iteratively reweighted least‐squares algorithm to solve the inversion problem to estimate unknown parameters (i.e., Lamé impedances, density, and fracture compliances) from PP‐ and PSV‐wave angle gathers along different azimuths. Synthetic tests reveal that the unknown parameters estimated using the joint inversion approach match true values better than those estimated using a PP‐wave amplitude inversion only. A real data test indicates that reasonable results for subsurface fracture detection are obtained from the joint inversion approach.     

具有离散裂缝空间分布的二维固体中地震波传播的有限差分模拟

结合有限差分方法和等效介质理论,模拟了离散分布裂缝介质中地震波的传播. 基于等效介质理论,利用二维有限差分实现封闭裂缝的离散分布;裂缝可以处理成固体岩石中的高度柔性界面,并可以用线性滑动或者位移间断模型进行裂缝的物理描述. 对于含有多组裂隙的破裂固体,其有效柔度可以认为是固体骨架背景柔度和裂缝附加柔度之和. 在一阶近似条件下,固体骨架和裂缝参数可以通过有效各向异性系数联系起来,有效各向异性系数决定了各向异性（裂缝效应）对于地震波传播的影响. 通过与射线理论方法的对比检验,说明本文提出的模拟方法的有效性,并通过几个数值算例说明本方法可有效模拟不同的裂缝分布效应. 结果表明,即使在裂缝密度很小的情况下,具有相同裂缝密度的不同的空间分布可以产生不同的波场特征. 同时,也验证了不同裂缝尺度对波长的不同影响,以及裂缝尺度具有幂率分布（分形）时,尺度对波场的影响. 最后得出结论：在运用建立在等效介质理论基础上的地震各向异性概念来描述裂缝固体的特征时,要倍加小心,等效介质理论中尚未合理处理的裂缝尺度和空间分布对波的传播特征具有重要的影响.     

Seismic reflection and transmission coefficients at an air-water interface of saturated porous soil

Based on the modified Biot＇s theory of two-phase porous media, a study was presented on seismic reflection and transmission coefficients at an air-water interface of saturated porous soil media. The major differences between air-saturated soils and water-saturated soils were theoretically discussed, and the theoretical formulas of reflection and transmission coefficients at an air-water interface were derived. The characteristics of propagation and attenuation of elastic waves in air-saturated soils were given and the relations among the frequency, the angle of incidence and the reflection, transmission coefficients were analyzed by using numerical methods. Numerical results show that the propagation characteristic of the wave in air-saturated soils is great different from that in water-saturated soils. The frequency and the angle of incidence can have great influences on the reflection and transmission coefficients at interface. Some new cognition about the wave propagation is obtained and the study suggests that we may carefully pay attention to the influence of air on the dynamic analysis of seismic wave.     

Reflection and Transmission Coefficients for an Incident Plane Shear Wave at an Interface Separating Two Dissimilar Poroelastic Solids

Using Biot’s poroelasticity theory, we derive expressions for the reflection and transmission coefficients for a plane shear wave incident on an interface separating two different poroelastic solids. The coefficients are formulated as a function of the wave incidence angle, frequency and rock properties. Specific cases calculated include the boundary between water-saturated sand and water-saturated sandstone and the gas–water interface in sand. The results show a very different interface response to that of an incident P wave. Plane SV wave incidence does not significantly excite the Biot slow P wave if the frequency of the wave is below the transition frequency. Above this frequency, an incident plane SV wave can generate a mode-converted slow Biot P wave which is actually a normal propagating wave and not highly attenuating as in the usual (diffusive) case. For an incident SV wave onto a gas–water interface, even at very high frequency, there is no significant Biot second P wave produced. For small incident angles, the gas–water interface is essentially transparent. With increasing angles, there can arise an unusual "definitive angle" in the reflection/transmission coefficient curves which is related to the change of fluid viscosity on both sides of the interface and provides a possible new means for underground fluid assessment.     

Characterization of dipping fractures in a transverselyisotropic background

Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities.     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上，利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右，地壳厚度的减薄量主要体现在下地壳，剖面北段下地壳厚度约为30 km，剖面南段下地壳厚度仅为15 km左右；上地幔顶部局部位置P波速度值偏低，一般为76～78 km/s，反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育，其中莫霍强反射出现在景云桥下方；在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.     

Seismic characterization of vertical fractures described as general linear-slip interfaces

Fluid flow in many hydrocarbon reservoirs is controlled by aligned fractures which make the medium anisotropic on the scale of seismic wavelength. Applying the linear‐slip theory, we investigate seismic signatures of the effective medium produced by a single set of ‘general’ vertical fractures embedded in a purely isotropic host rock. The generality of our fracture model means the allowance for coupling between the normal (to the fracture plane) stress and the tangential jump in displacement (and vice versa). Despite its low (triclinic) symmetry, the medium is described by just nine independent effective parameters and possesses several distinct features which help to identify the physical model and estimate the fracture compliances and background velocities. For example, the polarization vector of the vertically propagating fast shear wave S₁ and the semi‐major axis of the S₁‐wave normal‐moveout (NMO) ellipse from a horizontal reflector always point in the direction of the fracture strike. Moreover, for the S₁‐wave both the vertical velocity and the NMO velocity along the fractures are equal to the shear‐wave velocity in the host rock. Analysis of seismic signatures in the limit of small fracture weaknesses allows us to select the input data needed for unambiguous fracture characterization. The fracture and background parameters can be estimated using the NMO ellipses from horizontal reflectors and vertical velocities of P‐waves and two split S‐waves, combined with a portion of the P‐wave slowness surface reconstructed from multi‐azimuth walkaway vertical seismic profiling (VSP) data. The stability of the parameter‐estimation procedure is verified by performing non‐linear inversion based on the exact equations.     

Reflection and refraction at an imperfectly bonded interface between poroelastic solid and cracked elastic solid

Porous solid is in contact with a cracked elastic solid at a plane interface between them. For the presence of vertically aligned microcracks, the elastic solid behaves transversely isotropic to wave propagation. The coefficients of elastic anisotropy depend on the crack density and crack porosity in the medium. A loose bonding is considered between the two solids so that a limiting case could be the welded contact. At the plane interface, the imperfection in welded bonding is represented by tangential slipping and, hence, results in the dissipation of a part of strain energy. Three types of waves propagate in an isotropic fluid-saturated porous medium, which are considered for incidence at the interface. Incidence of a wave results in three reflected waves and two refracted waves. Partition of incident energy among the reflected and refracted waves is studied for each incidence, varying from normal to grazing directions. Numerical example calculates the energy shares of reflected and refracted waves at the plane interface between water-saturated sandstone and basalt. These energy shares are computed and analyzed for different values of crack parameters as well as loose bonding parameter.     

Research Note: Experimental measurements of Q‐contrast reflections

While seismic reflection amplitudes are generally determined by real acoustical impedance contrasts, there has been recent interest in reflections due to contrasts in seismic‐Q. Herein we compare theoretical and modelled seismic reflection amplitudes for two different cases of material contrasts. In case A, we examine reflections from material interfaces that have a large contrast in real‐valued impedance () with virtually no contrast in seismic‐Q. In case B, we examine reflections from material interfaces that have virtually no contrast in but that have very large seismic‐Q contrasts. The complex‐valued reflection coefficient formula predicts non‐zero seismic reflection amplitudes for both cases. We choose physical materials that typify the physics of both case A and case B. Physical modelling experiments show significantly large reflections for both cases – with the reflections in the two cases being phase shifted with respect to each other, as predicted theoretically. While these modelling experiments show the existence of reflections that are predicted by theory, there are still intriguing questions regarding the size of the Q‐contrast reflections, the existence of large Q‐contrast reflections in reservoir rocks and the possible application of Q‐reflection analysis to viscosity estimation in heavy oilfields.     

用地震反射剖面研究汤阴地堑上地壳结构与断裂特征

汤阴地堑位于太行山脉与华北平原的过渡带,是太行山前重要的地质构造单元。为研究汤阴地堑上地壳结构和断裂特征,利用安阳市与新乡市活断层探测获得的深、浅地震反射剖面,结合研究区已有地震、地质资料,对汤阴地堑浅部结构、汤东断裂特征进行分析研究。结果表明,汤阴地堑北部与南部地壳结构差异明显。地堑北部是由汤东断裂控制的半地堑,地堑内反射震相丰富、反射波层组关系清晰,多组新生代、古生代沉积层强反射不整合地覆盖在自西向东倾伏的结晶基底反射T_g上,且随着深度增加,地层倾角增大,显示出明显的多期掀斜运动特征。地堑南部表现为由汤东、汤西断裂共同控制的断陷型地堑,地堑内新近系底界面反射波TN自西向东倾伏,其下为一些横向呈水平或东倾、延续性较短的反射震相。地堑南部与北部不同的反射震相特征表明,测线控制区域内的汤阴地堑沉积环境与运动特征可能有所差异。汤东断裂为走向NE,倾向NW的铲型正断层,汤西断裂为走向NE、倾向SE的正断层,北部中深层地震反射剖面上未发现该断层的存在。     

An approach to combined rock physics and seismic modelling of fluid substitution effects

The aim of seismic reservoir monitoring is to map the spatial and temporal distributions and contact interfaces of various hydrocarbon fluids and water within a reservoir rock. During the production of hydrocarbons, the fluids produced are generally displaced by an injection fluid. We discuss possible seismic effects which may occur when the pore volume contains two or more fluids. In particular, we investigate the effect of immiscible pore fluids, i.e. when the pore fluids occupy different parts of the pore volume. The modelling of seismic velocities is performed using a differential effective‐medium theory in which the various pore fluids are allowed to occupy the pore space in different ways. The P‐wave velocity is seen to depend strongly on the bulk modulus of the pore fluids in the most compliant (low aspect ratio) pores. Various scenarios of the microscopic fluid distribution across a gas–oil contact (GOC) zone have been designed, and the corresponding seismic properties modelled. Such GOC transition zones generally give diffuse reflection regions instead of the typical distinct GOC interface. Hence, such transition zones generally should be modelled by finite‐difference or finite‐element techniques. We have combined rock physics modelling and seismic modelling to simulate the seismic responses of some gas–oil zones, applying various fluid‐distribution models. The seismic responses may vary both in the reflection time, amplitude and phase characteristics. Our results indicate that when performing a reservoir monitoring experiment, erroneous conclusions about a GOC movement may be drawn if the microscopic fluid‐distribution effects are neglected.     

长白山天池火山区三维地壳结构层析成像

利用长白山天池火山区三维空间深地震测观测系统所采集的反射P波走时资料,采用层析成像技术,重建了该区地壳界面构造形态和速度分析图像。界面成像结果表明：研究区地壳界面总体上显示了由北西向东南加深的趋势;马鞍山—三道白河断裂和富尔河—红旗河断裂是本区两条主要的深部构造,尤其是马鞍山—三道白河断裂,北北东向穿越天池火山口,其两侧的地壳界面存在明显的错断,预示了该处地壳厚度陡变或深大断裂带的存在,速度成像结果显示在10km深度,明显的P波低速异常分布在天池周围;15km深度上它表现为一个近南北向的P波低速异常条带,其延展尺度南北向为80~90km,东西向30~40km;随着深度的增加,P波低速异常分布在天池西侧,其尺度有明显的缩小,分布范围更加集中,而且低速扰动幅度更大。这种P波速度异常的变化图像在一定程度上反映了天池火山口下方壳内岩浆系统的空间分状况。     

THE SCATTERING OF SH-CHANNEL WAVES BY A FAULT IN A COAL SEAM*

The imaging of faults in coal seams by the in-seam seismic method has now become standard practice. In the UK over 300 surveys have been undertaken and the technique is now part of the exploration arsenal of colliery planners. From these users comes the pressure for two major improvements, namely an increase in range and target identification. This paper is directed towards the latter problem. It has long been recognized that the reflected channel waves must contain information on the fault structure that caused the reflection, and model experiments have been undertaken to investigate the reflection process. Only recently, however, have attempts been made to quantify the reflection process. Calculations using both the finite-difference and finite-element techniques have been carried out, and estimates of the reflection coefficient as a function of frequency have been obtained. The object of this paper is to extend these considerations by calculating analytically the scattering matrix of an SH-channel wave after interaction at a fault plane. The scattering matrix is calculated as a function of frequency, hade angle, and fault throw. The method employed is based on the decomposition of the incident SH-channel wave into Fourier components, the calculation of plane wave reflection and transmission coefficients within the constraints of geometrical acoustics, and finally the synthesis of the scattering matrix by application of the Helmholtz-Kirchhoff integral. The calculation throughout is restricted to normal modes.     

地震信号随机噪声压制的双树复小波域双变量方法

有效地压制地震信号中的噪声是地震信号解释和后续处理的重要环节之一.本文建立两种双树复小波域双变量模型对地震信号中的随机噪声进行压制.地震信号经双树复小波变换后,同一方向实部与虚部系数、实部(或虚部)系数与对应的模之间存在较强的相关性.鉴于此,对同一方向实部与虚部小波系数建立双变量模型,从含噪地震信号小波系数中估计原始信号的小波系数,再基于双树复小波逆变换重构得到降噪后的地震信号.进一步对同一方向实部(或虚部)系数与对应的模建立双变量模型,得到地震信号随机噪声压制的第二种双树复小波域双变量方法.最后对合成地震记录和实际地震资料中的随机噪声进行压制的实验结果证实本文两种方法都能够有效地压制地震信号中的随机噪声.