Numerical study of seismic scattering and waveguide excitation in faulted coal seams

Finite‐difference P‐SV simulations of seismic scattering characteristics of faulted coal‐seam models have been undertaken for near‐surface P‐ and S‐wave sources in an attempt to understand the efficiency of body‐wave to channel‐wave mode conversion and how it depends on the elastic parameters of the structure. The synthetic seismograms clearly show the groups of channel waves generated at the fault: one by the downgoing P‐wave and the other by the downgoing S‐wave. These modes travel horizontally in the seam at velocities less than the S‐wavespeed of the rock. A strong Airy phase is generated for the fundamental mode. The velocity contrast between the coal and the host rock is a more important parameter than the density contrast in controlling the amplitude of the channel waves. The optimal coupling from body‐wave energy to channel‐wave energy occurs at a velocity contrast of 1.5. Strong guided waves are produced by the incident S‐sources for source angles of 75° to 90° (close to the near‐side face of the fault). As the fault throw increases, the amplitude of the channel wave also increases. The presence of a lower‐velocity clay layer within the coal‐seam sequence affects the waveguiding characteristics. The displacement amplitude distribution is shifted more towards the lower‐wavespeed layer. The presence of a ‘washout’ zone or a brecciated zone surrounding the fault also results in greater forward scattering and channel‐wave capture by the coal seam.     

Impact–entrainment relationship in a saltating cloud

The problem of impact–entrainment relationship is one of the central issues in understanding saltation, a primary aeolian transport mode. By using particle dynamic analyser measurement technology the movement of saltating particles at the very near‐surface level (1 mm above the bed) was detected. The impacting and entrained particles in the same impact–entrainment process were identified and the speeds, angle with respect to the horizontal, and energy of the impacting and entrained sand cloud were analysed. It was revealed that both the speed and angle of impacting and entrained particles vary widely. The probability distribution of the speed of impacting and entrained particles in the saltating cloud is best described by a Weibull distribution function. The mean impact speed is generally greater than the mean lift‐off speed except for the 0·1–0·2 mm sand whose entrainment is significantly influenced by air drag. Both the impact and lift‐off angles range from 0° to 180°. The mean lift‐off angles range from 39° to 94° while the mean impact angles range from 40° to 78°, much greater than those previously reported. The greater mean lift‐off and especially the mean impact angles are attributed to mid‐air collisions at the very low height, which are difficult to detect by conventional high‐speed photography and are generally ignored in the existing theoretical simulation models. The proportion of backward‐impacting particles also evidences the mid‐air collisions. The impact energy is generally greater than the entrainment energy except for the 0·1–0·2 mm sand. There exists a reasonably good correlation of the mean speed, angle and energy between the impacting and entrained cloud in the impact–entrainment process. The results presented in this paper deserve to be considered in modelling saltation. Copyright © 2002 John Wiley & Sons, Ltd.     

Preliminary analysis on the source properties and seismogenic structure of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>7.0 Jiuzhaigou earthquake

At GMT time 13:19, August 8, 2017, an M_s7.0 earthquake struck the Jiuzhaigou region in Sichuan Province, China, causing severe damages and casualties. To investigate the source properties, seismogenic structures, and seismic hazards, we systematically analyzed the tectonic environment, crustal velocity structure in the source region, source parameters and rupture process, Coulomb failure stress changes, and 3-D features of the rupture plane of the Jiuzhaigou earthquake. Our results indicate the following: (1) The Jiuzhaigou earthquake occurred on an unmarked fault belonging to the transition zone of the east Kunlun fault system and is located northwest of the Huya fault. (2) Both the mainshock and aftershock rupture zones are located in a region where crustal seismic velocity changes dramatically. Southeast to the source region, shear wave velocity at the middle to lower crust is significantly low, but it rapidly increases northeastward and lies close to the background velocity across the rupture fault. (3) The aftershock zone is narrow and distributes along the northwest-southeast trend, and most aftershocks occur within a depth range of 5–20 km. (4) The focal mechanism of the Jiuzhaigou earthquake indicates a left-lateral strike-slip fault, with strike, dip, and rake angles of 152°, 74° and 8°, respectively. The hypocenter depth measures 20 km, whereas the centroid depth is about 6 km. The co-seismic rupture mainly concentrates at depths of 3–13 km, with a moment magnitude (M_w) of 6.5. (5) The co-seismic rupture also strengthens the Coulomb failure stress at the two ends of the rupture fault and the east segment of the Tazang fault. Aftershocks relocation results together with geological surveys indicate that the causative fault is a near vertical fault with notable spatial variations: dip angle varies within 66°–89° from northwest to southeast and the average dip angle measures ~84°. The results of this work are of fundamental importance for further studies on the source characteristics, tectonic environment, and seismic hazard evaluation of the Jiuzhaigou earthquake.     

Seismic migration and absorbing boundaries with a one-way wave system for heterogeneous media1

A first-order one-way wave system has been created based on characteristic analysis of the acoustic wave system and optimization of the dispersion relation. We demonstrate that this system is equivalent to a third-order scalar partial-differential equation which, for a homogeneous medium, reduces to a form similar to the 45° paraxial wave equation. This system describes accurately waves propagating in a 2D heterogeneous medium at angles up to 75°. The one-way wave system representing downgoing waves is used for a modified reverse time migration method. As a wavefield extrapolator in migration, the downgoing wave system propagates the reflection events backwards to their reflectors without scattering at the discontinuities in the velocity model. Hence, images with amplitudes proportional to reflectivity can be obtained from this migration technique. We present examples of the application of the new migration method to synthetic seismic data where P-P reflections P-SV converted waves are present. Absorbing boundaries, useful in the generation of synthetic seismograms, have been constructed by using the one-way wave system. These boundaries absorb effectively waves impinging over a wide range of angles of incidence.     

非一致地震波动输入下大型钢筋混凝土框架结构-地基体系地震响应规律

为探索非一致地震波动输入对大型钢筋混凝土框架结构地震响应的影响,基于OpenSees软件平台建立二维钢筋混凝土框架结构\|地基动力相互作用有限元模型。将El-Centro地震波按P波波形分别以0°、15°、30°和35°角入射该有限元模型进行计算,对比分析框架柱内力和楼层层间位移的地震响应。研究发现非一致地震波输入方法对于大型钢筋混凝土框架结构建筑动力响应影响明显,随着地震波入射角的增大,钢筋混凝土框架结构底层柱的轴力幅值减小,剪力幅值增大,而弯矩幅值变化较小,楼层层间位移幅值也随之增大。研究结果对于大型钢筋混凝土框架结构抗震设计具有参考意义。     

Re-examination of the earth's free oxcillations excited by the Kamchatka earthquake of November 4, 1952

Benioff's suggestion that the 58-min period sinusoidal oscillation found on a Pasadena strain seismogram after the Kamchatka earthquake of November 4, 1952 may represent the earth's gravest normal mode is re-examined in terms of a slow large-scale post-seismic deformation. The mechanism and the seismic moment of the main shock of the Kamchatka earthquake are determined by using the amplitude and the initial phase of G₂ and R₂ recorded at Pasadena and R₆ recorded at Palisades. By constraining the dip angle and the strike of the fault at 30° (towards NW) and N34°E, respectively, on the basis of the geometry of the Benioff zone, the slip angle is determined as 110° which represents 74% thrust and 26% right-lateral faulting. The direction of the slip angle agrees with the slip direction of the Pacific plate. A seismic moment of 3.5 · 10²⁹ dyn cm is obtained. If a fault area of 650 · 200 km² is assumed, an average dislocation of 5 m is obtained. Spectral analyses of the Pasadena strain records show that the 58-min sinusoidal oscillation in fact consists of a spectral peak near 54 min which is very close to the ₀S₂ mode and other high-frequency peaks which can be correlated to the earth's normal modes. The records from two independent recording galvanometers correlate with each other very well, indicating that the recorded oscillation represents a real strain and not instrumental noise. The phase relation between the NS and EW components is consistent with the strain field associated with ₀S₂ mode. Although these results provide positive evidence for a slow post-seismic deformation, the cause of the abrupt termination of the oscillation and the excitation mechanism remain unresolved.     

地震各向异性流体检测技术在胜利油田垦71井区的应用(英文)

根据Chapman理论模型,在各向异性介质(如HTI介质)中,当入射角在0-45。范围内,慢横波会发生较大的衰减和频散,且对流体粘度敏感,而P波和快横波则比较小。对于沿裂隙法向传播的慢横波,其振幅受流体影响很大。因此,在P波响应对流体不敏感的情况下,可利用慢横波来获得裂隙型油气藏的流体信息。本文分析了胜利油田垦71地区三维三分量地震数据,检测出的慢横波振幅和旅行时异常与该区的测井资料十分吻合。分析结果还发现,与含油区相比,含水区会产生更高的横波分裂。在含水区,慢横波振幅会产生明显变化,而在含油区则几乎没有变化。     

Earthquake strong motions recorded by a large near-source array of digital seismographs

Specially designed arrays of strong-motion seismographs near to the earthquake source are required for seismological and engineering studies of the generation and near-field properties of seismic waves. The first such large digital array, called SMART 1 (with radius 2 km and 37 accelerometers), to record substantial ground motion (up to 0·24g horizontal acceleration) became operational in late 1980 in a highly seismic region of Taiwan. During the first 6 months of operation, SMART 1 recorded nine earthquakes with magnitudes ranging from M_L 3·8 to 6·9. Three were located directly below the array at focal depths of 59 to 76 km. The remaining six had shallow depths and epicentral distances from 7 to 193 km. Digital records from 27 three-component accelerographs were obtained from a magnitude 6·9 (M_L) local earthquake on 29 January 1981. Representative measurements are described of seismic wave coherency and power spectrum as a function of wave number, frequency, azimuth of propagation and wave type. Acceleration waveforms varied significantly across the array for each event. On average, peak acceleration of horizontal components was about three times that of the vertical component. Relative spectral changes from earthquake to earthquake were large.     

Impact of incident angles of P waves on the dynamic responses of long lined tunnels

The impact of the incident angle of earthquake motion on the seismic response of the long lined tunnels is studied. Based on the time‐domain finite element method with the viscous‐spring artificial boundary condition, the earthquake motion of oblique incidence is transformed into the equivalent nodal forces acting on the truncated boundary of finite element model. In the present work, the formulas of equivalent nodal forces for the plane P wave with arbitrary incident angle are deduced and implemented into the commercial software abaqus   1 . The effectiveness of the formulas and its implementation are demonstrated by two numerical examples with the reference solutions. The proposed method is applied to investigate the seismic responses of the long lined tunnels under the obliquely incident P waves. The numerical results indicate that the seismic responses of the long lined tunnels are highly affected by the incident angles of P waves. Copyright © 2016 John Wiley & Sons, Ltd.     

WAVE FIELD EXTRAPOLATION TECHNIQUES FOR INHOMOGENEOUS MEDIA WHICH INCLUDE CRITICAL ANGLE EVENTS. PART I: METHODS USING THE ONE-WAY WAVE EQUATIONS*

Part I of this series starts with a brief review of the fundamental principles underlying wave field extrapolation. Next, the total wave field is split into downgoing and upgoing waves, described by a set of coupled one-way wave equations. In cases of limited propagation angles and weak inhomogeneities these one-way wave equations can be decoupled, describing primary waves only. For large propagation angles (up to and including 90°) an alternative choice of sub-division into downgoing and upgoing waves is presented. It is shown that this approach is well suited for modeling as well as migration and inversion schemes for seismic data which include critical angle events.     

Effect of the angle of seismic incidence on the fragility curves of bridges

This study examines the effect of the angle of seismic incidence θ on the fragility curves of bridges. Although currently, fragility curves of bridges are usually expressed only as a function of intensity measure of ground motion (IM) such as peak ground acceleration, peak ground velocity, or S_a(ω₁), in this study they are expressed as a function of IM with θ as a parameter. Lognormal distribution function is used for this purpose with fragility parameters, median c_m and standard deviation ζ to be estimated for each value of θ chosen from 0 < θ < 360°. A nonlinear 3D finite element dynamic analysis is performed, and key response values are calculated as demand on the bridge under a set of acceleration time histories with different IM values representing the seismic hazard in Los Angeles area. This method is applied to typical straight reinforced concrete bridges located in California. The results are validated with existing empirical damage data from the 1994 Northridge earthquake. Even though the sample bridges are regular and symmetric with respect to the longitudinal axis, the results indicate that the weakest direction is neither longitudinal nor transverse. Therefore, if the angle of seismic incidence is not considered, the damageability of a bridge can be underestimated depending on the incidence angle of seismic wave. Because a regional highway transportation network is composed of hundreds or even thousands of bridges, its vulnerability can also be underestimated. Hence, it is prudent to use fragility curves taking the incident angle of seismic waves into consideration as developed here when the seismic performance of a highway network is to be analyzed. Copyright © 2012 John Wiley & Sons, Ltd.     

伴随等离子体密度下降的磁声波与辐射带电子的波粒相互作用及其散射效应

利用范阿伦卫星的高质量观测数据,我们报道了伴随等离子体密度下降的磁声波现象.通过选取分别发生于2013年7月26日(事件A)和2013年9月19日(事件B)的两个相应事件进行细致分析,我们开展试验粒子模拟计算了磁声波对辐射带电子的散射系数,并求解二维福克-普朗克扩散方程量化了磁声波散射导致的辐射带电子动态变化.结果表明,事件A中的磁声波的散射作用主要发生于投掷角范围为60°～80°、能量范围为20～200keV的辐射带电子,而事件B中的磁声波的散射作用主要发生于投掷角范围为50°～80°、能量范围为20～400keV的辐射带电子;两个事件中的磁声波均能导致辐射带电子的蝴蝶状投掷角分布,但是由于事件B的磁声波幅度更强,形成的电子蝴蝶状分布更明显.     

Application of common reflection angle migration for imaging deformation structures in an inner accretionary wedge,Nankai Trough,Japan

To better image deformation structures within the inner accretionary wedge of the Nankai Trough, Japan, we apply common reflection angle migration to a legacy two-dimensional seismic data set acquired with a 6 km streamer cable. In this region, many seismic surveys have been conducted to study the seismogenic zone related to plate subduction. However, the details of the accreted sediments beneath the Kumano forearc basin are still unclear due to the poor quality of seismic images caused by multiple reflections, highly attenuated signals, and possibly complex geological structures. Generating common image gathers in the subsurface local angle domain rather than the surface offset domain is more advantageous for imaging geological structures that involve complex wave paths and poor illumination. By applying this method, previously unseen structures are revealed in the thick accreted sediments. The newly imaged geometric features of reflectors, such as the folds in the shallow part of the section and the deep reflectors with stepwise discontinuities, imply deformation structures with multiple thrust faults. The reflections within the deep accreted sediments (approximately 5 km) are mainly mapped to far angles (30°–50°) in the common reflection angles, which correspond to the recorded offset distances greater than 4.5 km. This result indicates that the far offset/angle information is critical to image the deformation structures at depth. The new depth image from the common reflection angle migration provides seismic evidence of multiple thrust faults and their relationship with the megathrust fault that is essential for understanding the structure and evolution of the Nankai Trough seismogenic zone.     

Rupture process of the M L=4.1 earthquake in Huailai Basin on July 20, 1995

On July 20, 1995, an earthquake of M _L=4.1 occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, seismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocentres were precisely located. About 2 hours after the occurrence of the main shock, a smaller event of M _L=2.0 took place at 40.323°N, 115.447°E with a focal depth of 5.0 km, which is very close to the main shock. Using the M _L=2.0 earthquake as an empirical Green’s function, a regularization method was applied to retrieve the far-field source-time function (STF) of the main shock. Considering the records of HDSN are the type of velocity, to depress high frequency noise, we removed instrument response from the records of the two events, then integrated them to get displacement seismogram before applying the regularization method. From the 5 field stations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained from each seismic phases are in good agreement, showing that the M _L=4.1 earthquake consisted of two events. STFs from each station demonstrate an obvious “seismic Doppler effect”. Assuming the nodal plane striking 37° and dipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through a trial and error method, the following results were drawn: Both of the events lasted about 0.1 s, the rupture length of the first one is 0.5 km, longer than the second one which is 0.3 km, and the rupture velocity of the first event is 5.0 km/s, larger than that of the second one which is about 3.0 km/s; the second event took place 0.06 s later than the first one; on the fault plane, the first event ruptured in the direction γ=140° measured clockwise from the strike of the fault, while the second event ruptured at γ=80°, the initial point of the second one locates at γ=−100° and 0.52 km from the beginning point of the first one. Using far-field ground displacement spectrum measurement method, the following source parameters about the M _L=4.1 earthquake were also reached: the scalar earthquake moment is 3.3×10¹³ N·m, stress drop 4.6 MPa, rupture radius 0.16 km. Contribution No. 99FE2022, Institute of Geophysics, China Seismological Bureau. This study is supported by the Chinese Joint Seismological Science Foundation (95-07-411).     

Effect of fracture fill on frequency‐dependent anisotropy of fractured porous rocks

In fractured reservoirs, seismic wave velocity and amplitude depend on frequency and incidence angle. Frequency dependence is believed to be principally caused by the wave‐induced flow of pore fluid at the mesoscopic scale. In recent years, two particular phenomena, i.e., patchy saturation and flow between fractures and pores, have been identified as significant mechanisms of wave‐induced flow. However, these two phenomena are studied separately. Recently, a unified model has been proposed for a porous rock with a set of aligned fractures, with pores and fractures filled with two different fluids. Existing models treat waves propagating perpendicular to the fractures. In this paper, we extend the model to all propagation angles by assuming that the flow direction is perpendicular to the layering plane and is independent of the loading direction. We first consider the limiting cases through poroelastic Backus averaging, and then we obtain the five complex and frequency‐dependent stiffness values of the equivalent transversely isotropic medium as a function of the frequency. The numerical results show that, when the bulk modulus of the fracture‐filling fluid is relatively large, the dispersion and attenuation of P‐waves are mainly caused by fractures, and the values decrease as angles increase, almost vanishing when the incidence angle is 90° (propagation parallel to the fracture plane). While the bulk modulus of fluid in fractures is much smaller than that of matrix pores, the attenuation due to the “partial saturation” mechanism makes the fluid flow from pores into fractures, which is almost independent of the incidence angle.     

基于地震活动性资料估计海原断裂倾角

本文以海原断裂带区域活动构造为基础,将海原断裂划分为西、中、东3段.基于1999年12月26日至2010年7月26日间的精定位小震目录,估计了海原断裂带各段的倾角.考虑海原断裂沿走向可能存在南倾与北倾两种情况,将倾角的范围设置为.首先运用网格搜索法确定了平面断层模型,其次以特征深度节点为基础数据运用多项式构建了曲面断层模型.结果表明：当进行平面拟合时,海原断裂西段与中段、东段的倾向不同,西段为南倾,其倾角值为71°,而中段、东段为北倾,其倾角值分别为72°、65°,各段的倾角值均由地表以下8 km地震资料确定.当进行曲面拟合时,在8 km深度以内海原断裂西段、中段、东段的倾角均处于80°左右,即接近陡立.西段的倾角在深度为9 km处出现转换,之后倾角接近陡立;中段的倾角在深度为16 km处出现转换,之后倾角逐渐减小,当深度为18 km时倾角为30°;东段的倾角在深度为11 km处出现转换,倾角为42°,在深度为16 km处出现第二次转换,倾角为55°,之后倾角逐渐减小.结合震源机制解和大地测量观测资料反演拟合的合理性,验证了本文所估计倾角的可靠性.     

计算地震初至波视出射角方法

在理论上证明了使用地震初至波位移、速度或加速度记录资料在时间域或频率域计算视出射角,使视出射角计算不受位移记录及时间域限制。在干扰信号频率段与有效信号频率段存在差异时分别在时间域和频率域进行了视出射角计算,结果表明在频率域中得到的视出射角基本与理论值相符,与方位角无关;而在时间域中获得的视出射角与理论值波动很大。应用实际地震波记录资料计算了不同台站的视出射角,其结果同样表明在频率域中获得的视出射角比时间域视出射角更加稳定和可信。     

Dynamics of the midlatitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at sunrise

The behavior of the F2 layer at sunrise has been studied based on vertical-incidence ionospheric sounding data in Almaty (76°55′E, 43°15′N). Records with small amplitudes of electron density background fluctuations were selected in order to exactly estimate the onsets of a pronounced increase in the electron density at different altitudes. It has been indicated that the electron density growth rate is a function of altitude; in this case, the growth rate at the F2 layer maximum is much lower than such values at fixed altitudes of ～30–55 km below the layer maximum. The solar zenith angle (χ) and the blanketing layer thickness (h ₀) at the beginning of a pronounced increase in the electron density at altitude h are linearly related to the h value, and these quantities vary within ～90° < χ < 100° and 180 km < h ₀ < 260 km, respectively.     

Phase functions of polar mesospheric cloud ice as observed by the CIPS instrument on the AIM satellite

The Cloud Imaging and Particle Size (CIPS) instrument on the Aeronomy of the Ice in the Mesosphere (AIM) spacecraft is a 4-camera nadir pointed imager with a bandpass centered at 265 nm and a field of view of 120°×80°. CIPS observes polar mesospheric clouds (PMCs) against the sunlit Rayleigh-scattered background. At individual polar locations approximately 5 km×5 km in area, CIPS observes the same volume of air seven times over a range of scattering angles from about 35° to 150°. These multi-angle observations allow the identification and extraction of the PMC scattered radiance from the Rayleigh-scattered background. We utilize the fact that the former has a highly asymmetric phase function about 90° scattering angle, while the latter has a phase function that is symmetric. The retrieved PMC phase function can then be interpreted to obtain PMC particle size distributions. We describe a technique for identification of PMCs in the CIPS observations through the separation of the Rayleigh and PMC radiances. PMC phase function results are shown for the first season of CIPS observations. Assuming the particles are oblate spheroids with an axial ratio of 2, and a Gaussian distribution of width 14 nm, we find the phase functions are consistent with mean radii between 50 and 60 nm. These results are similar to those discussed by Hervig et al. [2009. Interpretation of SOFIE PMC measurements: cloud identification and derivation of mass density, particle shape, and particle size. J. Atmos. Sol. Terr. Phys., in review.] in this issue from the Solar Occultation for Ice Experiment (SOFIE) which also flies on the AIM satellite.     

Advanced sensitivity calibration of the Los Angeles strong motion array

Results are presented of recent sensitivity calibration of 76 accelerographs (SMA-1) of the Los Angeles Strong Motion Array. These have pendulum-like transducers and optical recording system. One characteristic of their design is off-axis sensitivity, which is magnified by transducer misalignment. A new calibration procedure was applied, which considers off-axis sensitivity and measures the angles of misalignment (φ and ψ), as well as the incident angle of the light beam onto the film (θ₀). These are required (1) for accurate estimation of sensitivity, and (2) for proper instrument correction of recorded accelerograms which considers also cross-axis sensitivity and misalignment. These effects are important near large acceleration peaks (approaching and exceeding 1g), e.g. like the ones recorded near the source of the 1994 Northridge earthquake (M_L=6·4). This earthquake was recorded by 65 stations of the Los Angeles Strong Motion Array, at epicentral distances from 2 to 85 km. Histograms showing distribution of the misalignment angles, light beam incidence angle θ₀ (for unloaded position) and the transducer sensitivities are presented. These indicate that the misalignment angles are typically 1–1·5°, but may also be 3–4°. Angle θ₀ (usually neglected), is mostly between ±8°, but may reach ±12°. Assuming θ₀=0 leads to systematically smaller values of the measured sensitivity (e.g. by ∼3% for θ₀=8° and ∼4% for θ₀=12°). Comparison of the newly measured sensitivities with those measured prior to installation (in 1979/1980), s^old, shows that, in general, the new values are systematically smaller. The difference is typically within 5 per cent, but in some cases is as large as 10 per cent. Other principal sources of the observed differences and their mechanisms are discussed. Those include long-term changes in the transducers (e.g. change of stiffness, reflected in changes of the natural frequency) and differences in the calibration procedure (e.g. errors associated with manual reading film records with tilt test data, and with transducer and instrument housing misalignment). The presented results may be considered typical of similar strong motion arrays worldwide. © 1998 John Wiley & Sons, Ltd.