Laboratory measurements of guided‐wave propagation within a fluid‐saturated fracture

A fluid‐saturated flat channel between solids, such as a fracture, is known to support guided waves—sometimes called Krauklis waves. At low frequencies, Krauklis waves can have very low velocity and large attenuation and are very dispersive. Because they propagate primarily within the fluid channel formed by a fracture, Krauklis waves can potentially be used for geological fracture characterization in the field. Using an analogue fracture consisting of a pair of flat slender plates with a mediating fluid layer—a trilayer model—we conducted laboratory measurements of the velocity and attenuation of Krauklis waves. Unlike previous experiments using ultrasonic waves, these experiments used frequencies well below 1 kHz, resulting in extremely low velocity and large attenuation of the waves. The mechanical compliance of the fracture was varied by modifying the stiffness of the fluid seal of the physical fracture model, and proppant (fracture‐filling high‐permeability sand) was also introduced into the fracture to examine its impact on wave propagation. A theoretical frequency equation for the trilayer model was derived using the poroelastic linear‐slip interface model, and its solutions were compared to the experimental results.     

黏弹性介质瑞雷波有限差分模拟与特性分析

本文通过数值模拟研究了介质黏弹性对瑞雷波传播的影响.模拟采用结合了交错Adams-Bashforth时间积分法、应力镜像法和多轴完美匹配层的标准交错网格高阶有限差分方案.通过模拟结果和理论结果对比,测试了方法的精度,验证了结果的正确性.在均匀半空间模型中,分别从波场快照、波形曲线及频散能量图三个角度,对黏弹性介质瑞雷波衰减和频散特性进行了详细分析.两层速度递增模型被用于进一步分析瑞雷波在黏弹性层状介质中的特性.结果表明：由于介质的黏弹性,瑞雷波振幅发生衰减,高频成分比低频成分衰减更剧烈,衰减程度随偏移距增大而增强;瑞雷波相速度发生频散,且随频率增大而增大,频散能量的分辨率有所降低;黏弹性波动方程中的参考频率,不会影响瑞雷波振幅衰减和相速度频散的程度,但决定了黏弹性和弹性介质瑞雷波相速度相等的频率位置.本研究有助于人们更好地理解地球介质中瑞雷波的行为,并为瑞雷波勘探的应用和研究提供了科学和有价值的参考.     

微动多阶瑞雷波SPAC系数反演方法及应用研究

为充分利用微动信号中基阶和高阶模式瑞雷波,本文研究了基于多阶瑞雷波SPAC系数直接反演的方法.该方法首先基于地层介质响应计算多阶瑞雷波的能量占比,考虑实际观测台阵有限台站个数对SPAC系数影响,正演计算多阶瑞雷波SPAC系数,再采用快速模拟退火算法对其反演以获得地下介质横波速度结构.在此基础上,本文通过数值模拟验证该方法的可靠性,分别选取三种典型地质模型,基于模式叠加算法合成理论微动信号,采用本文方法计算其理论多阶瑞雷波SPAC系数并反演,给出反演结果与真实模型对比.我们将该方法应用于上海中心城区的地质调查中,通过与钻探结果对比,进一步验证该方法的有效性.本文理论与实际应用研究表明,基于多阶瑞雷波SPAC系数直接反演的微动探测方法有助于提高反演结果的可靠性,尤其对含软硬夹层的复杂地层介质,可提高探测精度.     

Surface waves in media having lateral inhomogeneities

Summary On the basis of Alsop's method (1966), approximated reflection and transmission coefficients are obtained for surface waves obliquely incident at a vertical discontinuity. Rayleigh and Love waves are arranged in such a way that they form a homogeneous eigenfunction system and their mutual conversion into each other is compatible with this theory. The structure of the eigenvalue equation is responsible for the mathematical complications in many surface wave problems. On the basis of a computer routine, first result concerning the angular dependence of the reflection and transmission coefficients including mode conversion are presented. The problem of phase jump is discussed. A comparison of theoretical and experimental results shows the applicability of Alsop's method for 90°-corners, too.     

Rayleigh波频散曲线“交叉”及多模式耦合作用研究

Rayleigh波可以用来反演近地表结构,在工程物探、石油物探、地球内部结构探测中均有重要意义.数值计算得到的含低速层的层状介质对应的Rayleigh波频散曲线会出现看似“交叉”的现象,但是对于这种现象目前还没有进行系统的研究.事实上可以验证,有些看似交叉的频散曲线实际上不相交.改变低速层的厚度和横波速度发现低速层越明显（即低速层速度越低或层厚越厚）频散曲线越不容易相交.凡友华等在2007年提出频散曲线对应着四种基本模式,在频散曲线发生“交叉”现象的区域实际上存在两个以上模式的频散曲线.本文主要研究了存在R模和S₂模的区域内频散曲线的“交叉”现象.首先利用竖直本征振动曲线研究R模和S2模Rayleigh 波的振动特点,发现R模对应的本征振动主要集中在地表,随着深度变化能量快速衰减,S₂模对应的本征振动主要集中在第2层.研究“交叉点”附近频散点对应的本征振动曲线发现这一区域有些Rayleigh波同时具有R模和S₂模的振动特点,对应着一种耦合模式.通过对实例的研究发现,在“交叉点”附近,若两条频散曲线不发生交叉,则每条曲线对应的模式会发生R模和S₂模之间经由耦合模式的转变,本文称这种现象为两种模式发生耦合;若两条频散曲线相交,则同一条频散曲线上的Rayleigh波模式几乎相同,只是在离交点很近的区域会存在一些耦合模式,本文称此时两种模式不发生耦合.本文研究结果主要供Rayleigh波对低速层结构的反演研究参考.     

Phase and group velocities and Q of mantle Love and Rayleigh waves of the first two modes and their azimuthal dependences for the 1963 Kurile Islands earthquake

Phase and group velocities and Q of mantle Love and Rayleigh waves from the 1963 Kurile Islands earthquake (M_w = 8.5) were determined over 37 great circle paths by a time variable filtering technique, in a period range 100–500 s for the fundamental modes and 100–275 s for the first higher modes. The preliminary reference Earth model (PREM) explains reasonably well the average dispersion results for the fundamental Love and Rayleigh waves. There exists a small, but significant inconsistency between the observation and the model for the first higher Love and Rayleigh waves. The Q structure of PREM is inconsistent with the observation for the fundamental Love waves, but explains other observations reasonably well. The dispersion of each mode shows a clear azimuthal dependence from which the four azimuthal windows were established. The phase and group velocity measurements for each window were, in general, shown to be mutually consistent. The azimuthal variations are largest for the first higher Rayleigh waves, indicating strong lateral heterogeneity in the structure of the low velocity zone. The first of the four windows is characterized by the largest fraction of Precambrian shields and the second window by the largest fraction of normal oceans. A comparison of these two windows may give some insight into deep lateral heterogeneity between continents and oceans. The observed phase and group velocities of the first window are systematically higher than those of the second window for the fundamental Love and Rayleigh waves at periods up to 400 s, and for the first higher Love and Rayleigh waves up to 175 s. Their differences are greatest for the first higher Rayleigh waves and least for the fundamental Rayleigh waves. Although the fundamental Rayleigh waves show the least velocity differences, their persistence up to a period of longer than 300 s is in striking contrast with some of the pure path phase velocities derived earlier for continents and oceans. A set of models for continents and oceans. PEM-C and PEM-O are not consistent with our observation. The third azimuthal window is characterized by trench-marginal seas and the fourth window by mountainous areas, typically the Asian high plateaus from northern China to the Middle East through Tibet. A comparison of these two windows gives some information about deep structural differences between subduction zones and continental collision zones, both belonging to plate convergence zones. For the fundamental and the first higher Love waves, the phase and group velocities for the third window are markedly low, whereas those for the fourth window are somewhat comparable to those for the second window. Slow Rayleigh waves are evident for two windows, with the fourth window apparently being the slowest for the fundamental Rayleigh above 200 s and for the first higher Rayleigh. For the fundamental Rayleigh waves, the third window is very slow below 200 s, but becomes progressively fast as the period increases and tends to be the fastest window around 400 s, suggesting a deep seated high velocity anomaly beneath trench-marginal seas. The dispersion characteristics of the fourth window indicate a thick high velocity lid with an extensive low velocity zone beneath it. The shield-like lithosphere, coupled with an extensive low velocity zone, may be a characteristic feature of continental collision zones. The particle motion of the fundamental Love waves was found not to be purely transverse to a great-circle connecting the epicenter to a station. The departure from the purely transverse motion is systematic among different periods, different G arrivals (G₂, G₃,…) and different stations, which may be interpreted as being due to lateral refraction.     

Dispersion Energy Analysis of Rayleigh and Love Waves in the Presence of Low-Velocity Layers in Near-Surface Seismic Surveys

High-frequency surface-wave analysis methods have been effectively and widely used to determine near-surface shear (S) wave velocity. To image the dispersion energy and identify different dispersive modes of surface waves accurately is one of key steps of using surface-wave methods. We analyzed the dispersion energy characteristics of Rayleigh and Love waves in near-surface layered models based on numerical simulations. It has been found that if there is a low-velocity layer (LVL) in the half-space, the dispersion energy of Rayleigh or Love waves is discontinuous and ‘‘jumping’’ appears from the fundamental mode to higher modes on dispersive images. We introduce the guided waves generated in an LVL (LVL-guided waves, a trapped wave mode) to clarify the complexity of the dispersion energy. We confirm the LVL-guided waves by analyzing the snapshots of SH and P–SV wavefield and comparing the dispersive energy with theoretical values of phase velocities. Results demonstrate that LVL-guided waves possess energy on dispersive images, which can interfere with the normal dispersion energy of Rayleigh or Love waves. Each mode of LVL-guided waves having lack of energy at the free surface in some high frequency range causes the discontinuity of dispersive energy on dispersive images, which is because shorter wavelengths (generally with lower phase velocities and higher frequencies) of LVL-guided waves cannot penetrate to the free surface. If the S wave velocity of the LVL is higher than that of the surface layer, the energy of LVL-guided waves only contaminates higher mode energy of surface waves and there is no interlacement with the fundamental mode of surface waves, while if the S wave velocity of the LVL is lower than that of the surface layer, the energy of LVL-guided waves may interlace with the fundamental mode of surface waves. Both of the interlacements with the fundamental mode or higher mode energy may cause misidentification for the dispersion curves of surface waves.     

Rayleigh wave scattering at wedge corners with major wedge angles

Theoretical calculations were performed on reflected and transmitted Rayleigh waves and scattered body waves, in the case where a two-dimensional Rayleigh wave is incident to a wedge-shaped medium having a wedge angle between 250° and 290° and arbitrary value of Poisson ratio. The reflection and transmission coefficients of Rayleigh waves were also experimentally measured in cases of wedges with 190° to 330° wedge angles. The method of theoretical analysis and the techniques of experiment are based on those developed in our preceding research (W-1. W-2 and W-3). Compared with the results where the wedge angle is smaller than 180° (W-1 and W-2), all features show consistent variation with wedge angle.     

Analytical partial derivatives of the phase- and group velocities for Rayleigh waves propagating in a layer on a half-space

The dispersion relation for Rayleigh waves in a layer on a half-space is modified by introducing the quadratic wave number instead of the phase velocity. The implicit function theorem is then used to derive analytical formulae for the group velocity and for the phase- and group velocity partial derivatives with respect to the parameters of the medium. As two examples, the method is applied to the interpretation of dispersion curves for short-period Rayleigh waves observed in a sedimentary basin of the West Carpathians, and in the Moravo-Silesian region, Czech Republic.     

大别-苏鲁及其邻近地区基于背景噪声的瑞雷波和勒夫波相速度层析成像

采用与作者2014年发表的“大别-苏鲁及其邻近地区基于背景噪声的勒夫波群速度成像”文章相同的资料,用频时分析提取5 000余条瑞雷波和4 000余条勒夫波相速度频散曲线,反演得到了8—32 s的瑞雷波和勒夫波相速度分布图像．结果显示,瑞雷波与勒夫波相速度分布具有很好的一致性．8 s的相速度分布与地表构造特征相吻合,造山带与隆起区均表现为高速,盆地因其规模不同而显示不同程度的低速．随着周期的增大,大别 苏鲁的高速带由强变弱,但始终存在．16—24 s的高速可能主要受到中地壳高速的控制,而32 s的高速则可能与上地幔顶部的高速有关．比较大别造山带与苏鲁造山带的平均频散曲线,发现大别造山带和苏鲁造山带的勒夫波频散曲线均高于AK135模型计算的理论频散曲线,而瑞雷波则没有这一现象. 这可能意味着两个地区有比较强烈的径向各向异性．      

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

Approximation to Cutoffs of Higher Modes of Rayleigh Waves for a Layered Earth Model

A cutoff defines the long-period termination of a Rayleigh-wave higher mode and, therefore is a key characteristic of higher mode energy relationship to several material properties of the subsurface. Cutoffs have been used to estimate the shear-wave velocity of an underlying half space of a layered earth model. In this study, we describe a method that replaces the multilayer earth model with a single surface layer overlying the half-space model, accomplished by harmonic averaging of velocities and arithmetic averaging of densities. Using numerical comparisons with theoretical models validates the single-layer approximation. Accuracy of this single-layer approximation is best defined by values of the calculated error in the frequency and phase velocity estimate at a cutoff. Our proposed method is intuitively explained using ray theory. Numerical results indicate that a cutoffs frequency is controlled by the averaged elastic properties within the passing depth of Rayleigh waves and the shear-wave velocity of the underlying half space.     

Crustal structure beneath Liaoning province and the Bohai Sea and its adjacent region in China based on ambient noise tomography

The velocity structure of the crust beneath Liaoning province and the Bohai sea in China was imaged using ambient seismic noise recorded by 73 regional broadband stations. All available three-component time series from the 12-month span between January and December 2013 were cross-correlated to yield empirical Green's functions for Rayleigh and Love waves. Phasevelocity dispersion curves for the Rayleigh waves and the Love waves were measured by applying the frequencytime analysis method. Dispersion measurements of the Rayleigh wave and the Love wave were then utilized to construct 2D phase-velocity maps for the Rayleigh wave at8–35 s periods and the Love wave at 9–32 s periods,respectively. Both Rayleigh and Love phase-velocity maps show significant lateral variations that are correlated well with known geological features and tectonics units in the study region. Next, phase dispersion curves of the Rayleigh wave and the Love wave extracted from each cell of the 2D Rayleigh wave and Love wave phase-velocity maps,respectively, were inverted simultaneously to determine the3 D shear wave velocity structures. The horizontal shear wave velocity images clearly and intuitively exhibit that the earthquake swarms in the Haicheng region and theTangshan region are mainly clustered in the transition zone between the low-and high-velocity zones in the upper crust, coinciding with fault zones, and their distribution is very closely associated with these faults. The vertical shear wave velocity image reveals that the lower crust downward to the uppermost mantle is featured by distinctly high velocities, with even a high-velocity thinner layer existing at the bottom of the lower crust near Moho in central and northern the Bohai sea along the Tanlu fault, and these phenomena could be caused by the intrusion of mantle material, indicating the Tanlu fault could be just as the uprising channel of deep materials.     

Crustal Structure of the Korean Peninsula Using Surface Wave Dispersion and Numerical Modeling

Surface wave dispersion is studied to obtain the 1-D average velocity structure of the crust in the Korean Peninsula by inverting group- and phase-velocities jointly. Group velocities of short-period Rayleigh and Love waves are obtained from cross-correlations of seismic noise. Multiple-filter analysis is used to extract the group velocities at periods between 0.5 and 20 s. Phase velocities of Rayleigh waves in 10- and 50-s periods are obtained by applying the two-station method to teleseismic data. Dispersion curves of all group and phase velocities are jointly inverted for the 1-D average model of the Korean Peninsula. The resultant model from surface wave analysis can be used as an initial model for numerical modeling of observations of North Korean events for a velocity model appropriated to the Korean Peninsula. The iterative process is focused especially on the surface sedimentary layer in the numerical modeling. The final model, modified by numerical modeling from the initial model, indicates that the crust shear wave velocity increases with depth from 2.16 km/s for a 2-km-thick surface sedimentary layer to 3.79 km/s at a Moho depth of 33 km, and the upper mantle has a velocity of 4.70 km/s.     

中国东部海域地壳-上地幔瑞利波速度结构研究

为了进一步了解中国东部沿海及相邻海域的地壳－上地幔结构特征,对该区域的构造演化历史、地震活动及深部构造等方面研究提供一些基础资料,利用３１个数字地震台记录的高质量瑞利波资料,采用一种新的混合路径频散的网格反演方法（Ｏｃｃａｍ方法）,对中国东部海域瑞利波群速度横向不均匀分布进行了初步研究．根据反演得到的１０-１５０ｓ共３６个中心周期的群速度分布特征,以及几个典型地点的剪切波速度结构的深度变化,对研究区域内各构造单元的划分以及它们在速度结构和上地幔低速层埋深等方面的特征进行了讨论。     

CROSS-BOREHOLE OBSERVATION OF MODE CONVERSION FROM BOREHOLE STONELEY WAVES TO CHANNEL WAVES AT A COAL LAYER1

Conversion of borehole Stoneley waves to channel waves was observed in data from a seismic cross-borehole experiment conducted between wellbores penetrating a thin coal layer at 2022 m depth, near Rifle, Colorado. Traveltime moveout observations show that borehole Stoneley waves underwent partial conversion to channel waves at the coal layer. The channel waves were detected directly in an adjacent borehole 35 m away at receiver positions within the coal. Stoneley waves, subsequently produced by partial conversion of channel waves, were also detected at receiver positions located up to 50 m above and below the coal layer in the adjacent borehole. We infer the channel wave to be the first-higher Rayleigh mode by comparing the observed group velocity with theoretically derived dispersion curves. Identifying the conversion between borehole and stratigraphically guided waves is significant because coal penetrated by multiple wells may be detected without placing a transmitter or receiver in the coal itself.     

PARAMETER ESTIMATION AND FAULT DETECTION BY THREE-COMPONENT SEISMIC AND GEOELECTRICAL SURVEYS IN A COAL MINE*

Three-component seismic and geoelectrical in-mine surveys were carried out in Lyukobanya colliery near Miskolc, Hungary to determine the in situ petrophysical parameter distributions and to detect inhomogeneities in the coal seam. The seismic measurements comprise an underground vertical seismic profile, using body waves, and an in-seam seismic amplitude-depth distribution and transmission survey, using channel waves. The geoelectrical measurements are based on the drift- and seam-sounding method. Interval traveltime-, amplitude-, multiple-filter- and polarization analysis methods are applied to the seismic data. They lead to a five-layer model for the strata including the coal seam. The coal seam and two underlying beds act as a seismic waveguide. The layer sequence supports the propagation of both normal and leaky mode channel waves of the Love- and Rayleigh type. A calculation of the total reflected energy for each interface using Knott's energy coefficients shows that the velocity ranges of high reflection energy and of normal and leaky mode wavegroups coincide. The excitation of wavegroups strongly depends on the seismic source. A simultaneous inversion of a geoelectrical drift- and seam-sounding survey prevents misinterpretations of the seismic data by clearly identifying the low-velocity coal seam as a high-resistivity bed. Calculations of dispersion and sounding curves improve the resolution of the slowness and resistivity in each layer. Both diminished amplitudes and distortions in the polarization of transmission seismo-grams and decreasing resistivities in a geoelectrical pseudosection of the coal seam are related to an inhomogeneity. A calculation of synthetic seismograms for Love and Rayleigh channel waves with the finite-difference and the Alekseev-Mikhailenko method agrees well with the field data for the main features, i.e., particular arrivals in the wave train, wavegroups, velocities and symmetries or asymmetries. This in-mine experiment demonstrates that the simultaneous acquisition, processing and interpretation of seismic and geoelectrical data improve the lithological interpretation of petrophysical parameter distributions. Coal seam inhomogeneities can also be detected more reliably by the two independent surveys than by one alone.     

Estimation of the quality factor in shallow soil using surface waves

IntroductionVelocityanditsattenuationinformationiscloselylinkedwiththeoreticalstudiesonthegroundmovementsduringearthquake.Therecentstudy(Malagnini,1996)showedthatthevelocitystructureofshearwaveinshallowsoilabove30mplaysanimpoftantroletoestimatestfonggroundmotionofsite.However,itishardtopreciselymeasurethesoilstructuresanddynamiccharacteristics.First,theloosesoilabsorbstheseismicwaveswithhighfrequencies;Second,theeffectsoffocusinganddispersioncausedbylocallyinhomogeneoussitecannotbeneglectedin…     

Anomalous Rayleigh-Wave Propagation Along Oceanic Trench

A clear later phase of amplitude larger than the direct surface wave packet was observed at stations in Hokkaido, Japan, for several events of the December 1991 off-Urup earthquake swarm in the Kuril Islands region. From its particle motion, this phase is likely to be a fundamental Rayleigh wave packet that arrived with an azimuth largely deviated from each great-circle direction. As its origin, Nakanishi (1992) proposed that the sea-trench topography in this area as deep as 10 km may produce a narrow zone of low velocity for Rayleigh waves of periods around 15 sec. Following this idea, we compute ray paths and estimate how Rayleigh waves would propagate if we assume that lateral velocity variations are caused only by seafloor topography. We confirm that thick sea water in the trench indeed produces the phase velocity of Rayleigh waves to be smaller than in a surrounding area by the degree over 100%. Such a low-velocity zone appears only in a period range from 12 to 20 sec. Although this strong low-velocity zone disturbs the direction of Rayleigh wave propagation from its great circle, the overall ray paths are not so affected as far as an epicentre is outside this low-velocity zone, that is, off the trench axis. In contrast, the majority of rays are severely distorted for an event within the low-velocity zone or, in other words, in the neighborhood of the trench axis. For such an event, a part of wave energy appears to be trapped in this zone and eventually propagates outwards due to the curvature or bend of trench geometry, resulting in very late arriving waves of large amplitude with an incident direction clearly different from great circles. This phenomenon is observed only at a very limited period range around 16 sec. These theoretical results are consistent with the above mentioned observation of Nakanishi (1992).     

Validity of the long-wave approximation in periodically layered media

In seismic modelling, a stack of thin layers is often replaced by an effective equivalent anisotropic homogeneous slab. For waves with finite wavelength, this is an approximation, and the error thus introduced can be quantified by considering the relative error in the phase velocity between the layer stack and the effective medium. For periodic layering, the relative phase-velocity error can be expressed in closed form as a function of wavelength, reflection coefficients and layer thicknesses. By comparing the relative phase-velocity error with laboratory measurements and numerical simulations, we find that the difference in seismic response between a periodic layer stack and an equivalent effective medium depends not only on wavelength, but it also depends significantly on reflection coefficients and the ratio between layer thicknesses. For a 1% relative error in the phase velocity, and if all layers have the same thickness measured in vertical traveltime, we find that the wavelength must be larger than approximately three times the layer period for a reflection coefficient of 0.1, but this increases to 13 times the layer period for a reflection coefficient of 0.9, which is highly unrealistic in a geological setting.