地震波在垂向不均匀介质中传播问题的某些探讨

本文计算了含有高速夹层介质中首波的理论地震图。通过分析得到,当高速夹层薄到一定程度时,就会产生干涉型首波,从而从一个侧面证明了射线理论的局限性。通过对地震波反射—折射系数能量守恒关系的分析,探讨了反射—折射系数大于1的可能性。最后,介绍了一种计算垂向不均匀介质中拉梅问题理论地震图的数值方法——有限差分法。     

用波函数展开法求解界面圆孔的SH波散射问题

研究了ＳＨ波绕界面圆孔散射的波函数展开理论。由入射,反射和透射波组成的自由波场与孔的散射场叠加成总波场。将定义于两个半平面的散射场按一定方式延拓到全平面。通过Ｈａｎｋｅｌ－Ｆｏｕｒｉｅｒ型积分来满足界面与孔缘处的边界条件,得到了确定特定系数的封闭代数方程组。给出了具体算例,计算了孔缘处的动应力系数分布。     

横向各向同性地层(VTI)井孔声弹效应对弯曲波的影响

本文发展了建立在地层参考状态为各向同性介质假定下的现行井孔声弹性理论,就井外为横向各向同性面与井轴垂直的、具有9个独立三阶弹性模量的横向各向同性介质(VTI井况),水平面内受双轴应力作用下给出了一个简洁的与井内压力、应力差、应力和以及多极源偏振方位角有关的井孔弯曲波声弹公式,并且导出了平面纵、横波速度的声弹公式.数值考察了弯曲波速度之改变量的灵敏系数随频率的变化、受井外水平双轴应力作用时两种偏振的偶极弯曲波频散曲线以及对应不同方位径向偏振的横波速度.研究结果表明弯曲波声弹公式与5个二阶弹性系数以及7个独立的三阶弹性模量有关;而且由内压引起的井孔弯曲波声弹性公式中的三阶弹性模量仅与6个独立的弹性模量有关.横向各向同性介质井孔弯曲波速度的交叉现象仍是判断地应力存在的标志;一个重要的认识是受双轴应力作用的弯曲波速度变化在低频区主要与c₁₄₄和c₁₅₅两个三阶弹性模量有关,而且此认识与径向偏振的平面横波一致.在缺乏足够实验条件情况下,对VTI情况,以c₁₄₄,c₁₅₅和c₁₂₃三个独立的量进行测量,然后可暂不考虑三阶弹性模量的各向异性,建立简化的应力反演公式.反之,如果已知地层的地应力信息,由简化的声弹公式可以反演三阶弹性模量c₁₄₄,c₁₅₅和c₁₂₃.     

Determining maximum shear stress in confined substrate from elastic wave reflection coefficient

The relationship between the maximum shear stress in a substrate solid and the elastic wave reflection coefficient from the interface between the substrate solid and an overlying solid half-space is investigated. Both substrate and overlying solid media are assumed to be initially isotropic and stress-free. Then as the substrate is subjected to horizontal confined stresses it becomes anisotropic. It is shown that longitudinal and shear wave reflection coefficients are related to the degree of stress induced anisotropy in the substrate medium. From this relation the confined stress level and the maximum shear stress generated on the vertical planes of the substrate are estimated. Authors in their previous investigation computed plane wave reflection coefficient in a biaxially compressed solid substrate immersed in a fluid. This paper reports for the first time how the maximum shear stress in a biaxially compressed substrate medium can be measured from the plane wave reflection coefficients when the overlying medium is also a solid half-space.     

汶川地震断裂带东北端浅部结构的人工地震探测

结合汶川地震断裂带动态监测,利用快速响应探测系统,开展了断层带浅部结构人工地震探测.针对地震断裂带动态监测条件下的复杂波场和低信噪比的情况,在f-k波场分离的基础上,分别利用了折射波共中心点成像、面波速度反演、反射波叠加成像方法,进行了浅层断层和构造成像处理,并对处理结果进行了综合解释,给出了断裂带浅部断层分布和速度特征.为汶川地震龙门山断裂带东北端动态监测提供了基础结构信息,所发展的断裂带快速响应探测技术对于地震应急动态监测具有重要意义.     

Reflection and refraction of waves at the interface of an isotropic medium over a highly anisotropic medium

In this paper, we have considered the reflection and refraction of a plane wave at an interface between two half-spaces. The lower half-spaces is composed of highly anisotropic triclinic crystalline material and the upper half-space is homogeneous and isotropic. It has been assumed that due to incidence of a plane quasi-P (qP) wave, three types of waves, namely, quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated in the lower half space whereas P and S waves will be generated in the upper half space. The phase velocities of all the quasi waves have been calculated. It has been assumed that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Some specific relations have been established between directions of motion and propagation, respectively. The expressions for reflection coefficients of qP, qSV, qSH and refracted coefficients of P and SV waves are obtained. Results of reflection and refraction coefficients are presented.     

声波在两种多孔介质界面上的反射和透射

本文导出了声波在两种多孔介质界面上反射、透射的一般计算公式.作为例子,数值计算了P₁波入射于界面时,P₁、P₂和S波的反射、透射系数与声波频率、入射角等量之间的关系.结果表明,各种模式波的反射、透射系数与入射角、多孔介质性质有关,在Biot特征频率附近与频率有关,并用界面两侧的法向能流相等验证了结果的正确性.若把多孔介质当作均匀固体处理,将会得到显著不同的结果.     

Reflection and transmission of plane SH-waves in an initially stressed inhomogeneous anisotropic magnetoelastic medium

This study investigates the reflection and transmission of plane SH-waves in two semi-infinite anisotropic magnetoelastic media. The lower half-space is considered as initially stressed and inhomogeneous. The density of lower half-space is taken exponentially varying with depth. The solutions for half-spaces are obtained analytically. The expressions for reflection and transmission coefficient are obtained in the closed form subject to continuity conditions at the interfaces of anisotropic magnetoelastic half-spaces and the Snell’s law. It is found that these coefficients depend on the initial stress, inhomogeneity parameter, the magnetoelastic coupling parameter, and the angle at which wave crosses the magnetic field of the half-spaces. Numerical computations are performed for these coefficients for a specific model of two different anisotropic magnetoelastic half-spaces. The numerical results are illustrated by the graph of reflection and transmission coefficient versus the angle of incidence. In general, as the initial stress increases the reflection and transmission coefficient increases, the affect is more prominent for more than 10 GPa. Inhomogeneity in the density of the material also increases the reflection and transmission coefficient. The anisotropic magnetoelastic parameter and the angle at which the wave crosses the magnetic field for both the half-spaces have a quite significant effect on the reflection and transmission coefficient.     

考虑河谷场地效应的大跨CFST拱桥地震响应分析

本文对某特大跨上承式铁路钢管混凝土拱桥建立“V”型河谷场地有限元动力计算模型,结合时空解耦的时域动力有限元方法,研究“V”型河谷场地效应对其地震响应的影响,并考虑地震入射角、行波速度及地震动峰值等因素对场地效应的影响。结果表明:海拔高且坡度均匀的地形使地震波在反射过程中能量增大,海拔较低且坡度下陡上缓的地形使地震波在反射过程中能量减少,“V”型河谷场地效应使大部分钢管混凝土拱圈截面内力测量值放大作用明显;“V”型河谷场地地形效应与地震波的入射角、行波速度及加速度峰值等因素均有关联,且CFST拱肋面外受弯抗震性能受行波速度的影响较大。因此,钢管混凝土拱桥抗震设计中应注意场地地形效应及行波速度的影响。     

Analytical solutions for dynamic pressures of coupling fluid-solid-porous medium due to P wave incidence

Wave reflection and refraction in layered media is a topic closely related to seismology, acoustics, geophysics and earthquake engineering. Analytical solutions for wave reflection and refraction coefficients in multi-layered media subjected to P wave incidence from the elastic half-space are derived in terms of displacement potentials. The system is composed of ideal fluid, porous medium, and underlying elastic solid. By numerical examples, the effects of porous medium and the incident wave angle on the dynamic pressures of ideal fluid are analyzed. The results show that the existence of the porous medium, especially in the partially saturated case, may significantly affect the dynamic pressures of the overlying fluid.     

2.5D modelling, inversion and angle migration in anisotropic elastic media

2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D .     

界面和速度反射联合成像──理论与方法

依据波逆行原理重新推导了由地震波反射走时资料反演界面深度的关系式──走时对界面偏导数关系，它不仅适用于反射波，也适用于透射波．还给出了任意多个复杂界面情况下走时对界面偏导数关系的离散形式．另外，正交算子投影法被推广用来解决速度与界面的同时成像问题．数值模拟的结果表明，本文的方法是有效的．     

The seismic velocity distribution in the vicinity of a mine tunnel at Thabazimbi, South Africa

Analysis of the refracted arrivals on a seismic reflection profile recorded along the wall of a tunnel at an iron mine near Thabazimbi, South Africa, shows variations in P-wave velocity in dolomite away from the de-stressed zone that vary between 4.4 and 7.2 km/s, though values greater than 5.8 km/s predominate along most of the profile. The seismic velocities at the tunnel wall, however, vary between 4.2 and 5.2 km/s. Time–depth terms are in the range from 0.1 to 0.9 ms, and yield thicknesses of the zone disturbed by the tunnel excavations of between 2 and 9 m. The very low seismic velocities away from the tunnel wall in two regions are associated with alcoves or ‘cubbies’ involving offsets in the wall of up to 10 m. The large variations in seismic velocity resolved over distances less than 15 m with signals of wavelength around 6–9 m are attributed to variations in the sizes and concentrations of fracture systems and cracks, and in the degree of groundwater saturation of the fracture systems. The results suggest that seismic velocity variations from reflection surveys may also assist modelling studies of the stress regime in deep mines, particularly if both P and S wave velocity variations can be determined. The seismic velocity variations inferred also show that application of refraction static corrections in the processing of ‘in-mine’ seismic reflection profiles is as important as in surface surveys, because of the higher frequencies of the seismic energy recorded in the deep mine environment.     

速度光滑对反射地震信号的影响分析

光滑处理使得单界面成为非均匀薄层,界面反射转变为层反射.为了探讨光滑处理的影响,以平面波作为入射波场,首先利用过渡层反射系数推导了反射信号的理论公式,然后就非均匀薄层下反射系数的计算问题,给出了具体的实现算法,并通过与经典Epstein过渡层反射系数解析结果的对比说明了算法的精度.最后在单界面及其被光滑后界面的对比分析中,得出了几点重要结论：随着光滑次数的增加,反射信号的到时基本保持不变,而反射信号的主频与能量呈减少趋势,其中信号能量在低光滑次数的衰减速率明显大于高光滑次数.     

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.     

球面波的反射P波AVO分析

本文对球面波反射P波反射系数计算公式进行了推导,根据推导公式计算出第I类AVO和第Ⅲ类AVO介质中反射系数及相位随入射角的变化,并与基于平面波的Zoeppritz方程计算的结果进行对比.结果显示,对于第I类AVO介质,球面波反射系数在临界角附近较基于平面波的Zoeppritz方程计算结果更精确;同时受界面深度的影响较大,随深度的增加,球面波AVO趋势接近平面波AVO.最后,通过数值模拟技术,对模型数据和实际数据进行模拟,对平面波AVO道集和球面波AVO道集进行了验证分析.浅层的第I类AVO现象,在近临界角和超临界角处折射引起反射系数与相位的变化较大,对于第Ⅲ类AVO现象,由于不存在临界角问题,球面波模拟结果与基于平面波的Zoeppritz计算结果差别较小.上述计算分析,可为实际资料的大偏移距道集的AVO分析提供理论基础.     

Introduction to seismic travel time methods in anisotropic media

Summary Section 1 (and 11) develops the concepts of the front velocity, the front gradient, the travel time in space and on seismometric profiles, the profile velocity and the profile gradient in connection with the propagation of the fronts of elastic waves in solid isotropic and anisotropic media. The sectional velocity and the sectional gradient are defined in terms of the motion of the curve of intersection of a front with a fixed surface. Section 2 (and 12) relates the coefficients of elasticity of the medium, the front types, and their respective rays. In section 12, the theory of fronts of arbitrary shape and of the corresponding rays for any anisotropic, homogeneous or inhomogeneous solid medium is summarized. In section 3 (and 13), the law of reflection and refraction of fronts on surfaces of discontinuity of arbitrary shape is presented. Sections 4 to 6 (and 14 to 16) treat some elementary applications of seismic travel time methods to homogeneous, uniaxially anisotropic media (=transverse isotropy) in greater detail. In section 4 (and 14), the travel time of a direct front generated by a point source is considered and it is shown how the coefficients of elasticity of the medium can be found based on travel time measurements. The seismic prospection of a plane reflector and of a reflecting boundary of arbitrary shape and position are discussed in section 5 (and 15). In section 6 (and 16), the seismic refraction method is used to locate a plane boundary between a homogeneous, uniaxially anisotropic and a homogeneous isotropic medium, where the boundary is perpendicular or at an arbitrary angle to the direction of anisotropy.     

Sensitivity analysis and application of time‐lapse full‐waveform inversion: synthetic testing and field data example from the North Sea,Norway

Time‐lapse refraction can provide complementary seismic solutions for monitoring subtle subsurface changes that are challenging for conventional P‐wave reflection methods. The utilization of refraction time lapse has lagged behind in the past partly due to the lack of robust techniques that allow extracting easy‐to‐interpret reservoir information. However, with the recent emergence of the full‐waveform inversion technique as a more standard tool, we find it to be a promising platform for incorporating head waves and diving waves into the time‐lapse framework. Here we investigate the sensitivity of 2D acoustic, time‐domain, full‐waveform inversion for monitoring a shallow, weak velocity change (?30 m/s, or ?1.6%). The sensitivity tests are designed to address questions related to the feasibility and accuracy of full‐waveform inversion results for monitoring the field case of an underground gas blowout that occurred in the North Sea. The blowout caused the gas to migrate both vertically and horizontally into several shallow sand layers. Some of the shallow gas anomalies were not clearly detected by conventional 4D reflection methods (i.e., time shifts and amplitude difference) due to low 4D signal‐to‐noise ratio and weak velocity change. On the other hand, full‐waveform inversion sensitivity analysis showed that it is possible to detect the weak velocity change with the non‐optimal seismic input. Detectability was qualitative with variable degrees of accuracy depending on different inversion parameters. We inverted, the real 2D seismic data from the North Sea with a greater emphasis on refracted and diving waves’ energy (i.e., most of the reflected energy was removed for the shallow zone of interest after removing traces with offset less than 300 m). The full‐waveform inversion results provided more superior detectability compared with the conventional 4D stacked reflection difference method for a weak shallow gas anomaly (320 m deep).     

Reflection and refraction of plane harmonic waves at an interface between elastic solid and porous solid saturated by viscous liquid

A general solution of Biot's field equations governing small motions of a porous solid saturated by viscous liquid is employed to study the reflection and refraction at the interface between an elastic solid and a liquid-saturated porous solid. The incident wave is assumed to be plane and homogeneous, propagating through the isotropic elastic solid. The poroelastic solid is considered to be a dissipative one. Amplitude and energy ratios are computed numerically for a particular model. With first-order corrections for the porosity of solid and viscosity of liquid, the limiting cases of low and high frequencies are computed.     

THE GENERALIZED SPECTRAL FUNCTION OF A PLANE LAYERED MEDIUM*

The spectral function of a plane layered medium, which represents the net downgoing energy in the first layer due to a normally incident impulsive plane wave, plays an important role in the solution of the one-dimensional inverse problem in reflection seismology. Hitherto the extension to non-normal incidence was known only for a medium with free surface. By giving the extension for arbitrary surface reflection coefficients, this paper fills a gap.