Temporal change in seismic wave attenuation using highly stable vibration sources

We developed a method to detect attenuation changes during seismic wave propagation excited by precisely controlled artificial seismic sources, namely Accurately Controlled Routinely Operated Signal System (ACROSS), and applied it to monitor the temporal changes for in situ data collected by previous studies. Our method, together with the use of the ACROSS sources, is less susceptible to noise level changes, from which conventional methods such as envelope calculation suffer. The method utilizes the noise level that is independently estimated in the frequency domain and eliminates the influence of the noise from the observed signal. For performance testing, we applied this method to a dataset that was obtained in an experiment at Awaji Island, Japan, from 2000 to 2001. We detected a change in amplitude caused by rainfall, variation in atmospheric temperature, and coseismic ground motions. Among them, coseismic changes are of particular interest because there are limited studies on coseismic attenuation change, in contrast to many studies on coseismic velocity decrease. At the 2000 Western Tottori earthquake (MW = 6.6, epicenter distance of 165 km), a sudden decrease in amplitude of up to 5% was observed. The coseismic amplitude reduction and its anisotropic characteristics, which showed a larger reduction in the direction of the major axis of velocity decrease, were consistent with the opening of fluid-filled cracks, as proposed by previous studies. The $$\Delta {Q}^{-1}$$ Δ Q - 1 corresponding to the amplitude change gives similar values to those reported in previous studies using natural earthquakes. Graphical Abstract

Seismic Wave Field Model Excited by a Finite Long Cylindrical Charge

The amplitude-frequency characteristic of a seismic wave excited by explosion sources directly affects the accuracy of seismic exploration. To reveal the effect law related to a cylindrical charge, the research proposes a seismic wavefield model excited by a long cylindrical charge. According to the characteristics of the blasting cavity generated by a finite length cylindrical charge, the seismic wavefield characteristics of a cylindrical charge excitation is obtained by superposing the seismic wavefield excited by a series of spherical charges. Numerical simulation results show that the calculation error of the blasting cavity characteristics of the theoretical model is within 10%. The comparison with field experimental results shows that the error of the model is within 9.4%. The velocity field of the excited seismic wave is almost the same as that of the spherical charge when the explosion distance to the cylindrical charge with finite length is 16-21 times longer than the charge length, but the frequency of the seismic wave is 30% higher than for a spherical charge. Moreover, the explosive velocity has a certain influence on the amplitude-frequency characteristic of the seismic wave excited by the cylindrical charge. The established theoretical model can accurately describe the amplitude-frequency characteristics of the seismic wavefield excited by a cylindrical charge with finite length.

Estimation of Seismic Wave Attenuation from 3D Seismic Data: A Case Study of OBC Data Acquired in an Offshore Oilfield

Previous studies performed in Abu Dhabi oilfields, United Arab Emirates, revealed the direct link of seismic wave attenuation to petrophysical properties of rocks. However, all those studies were based on zero offset VSP data, which limits the attenuation estimation at one location only. This is due to the difficulty of estimating attenuation from 3D seismic data, especially in carbonate rocks. To overcome this difficulty, we developed a workflow based on the centroid frequency shift method and Gabor transform which is optimized by using VSP data. The workflow was applied on 3D Ocean Bottom Cable seismic data. Distinct attenuation anomalies were observed in highly heterogeneous and saturated zones, such as the reservoirs and aquifers. Scattering shows significant contribution in attenuation anomalies, which is unusual in sandstones. This is due to the complex texture and heterogeneous nature of carbonate rocks. Furthermore, attenuation mechanisms such as frictional relative movement between fluids and solid grains, are most likely other important causes of attenuation anomalies. The slight lateral variation of attenuation reflects the lateral homogeneous stratigraphy of the oilfield. The results demonstrate the potential of seismic wave attenuation for delineating heterogeneous zones with high fluid content, which can substantially help for enhancing oil recovery.

Dynamic Response and Failure Process of a Counter-Bedding Rock Slope under Strong Earthquake Conditions

There are massive landslides and potential landslides along the Three Rivers Basin in the Qinghai–Tibet Plateau, which pose a serious threat to the Sichuan–Tibet Railway. A normal shaking table model test was conducted to study the dynamic characteristics and dynamic response of a symmetrical counter-bedding rock slope based on the Zongrong Village landslide. The influences of the dynamic parameters, seismic wave type, and a weak intercalated layer on the slope’s dynamic response were considered. The results showed symmetry between the growth trend of the acceleration amplification factor and other research results. When the input wave amplitude was constant, the acceleration amplification factor increased at first and then decreased as the frequency increased. When the input frequency was near the slope’s natural frequency, the acceleration amplification factor increased at first and then decreased with an increase in the input amplitude and reached the maximum value at 0.3 g. The acceleration amplification factor increased linearly with height in the vertical direction inside the slope but increased slowly at first and then sharply along the slope surface, reaching the maximum value at the slope’s top and exhibiting an obvious “elevation effect”. When sinusoidal waves, Wolong waves, and Maoxian waves with the same amplitude were input, the slope’s amplification effect on the bedrock wave was more obvious. The weak intercalated layer showed the phenomenon of “thin layer amplification” and “thick layer attenuation” in response to the input seismic wave. The slope’s failure process can be roughly divided into three stages: (1) the formation of tensile cracks at the top and shear cracks at the toe; (2) the extension of cracks and the sliding of the slope-surface block; (3) the formation of the main sliding surface.

Effects of Three-Directional Seismic Wave on Dynamic Response and Failure Behavior of High-Steep Rock Slide

The landslide triggered by earthquakes can cause severe infrastructure losses or even fatalities. The high-steep rock slide is the most common type of landslide in the earthquake area. In an earthquake, the ground moves randomly in all directions, two horizontal directions (East-West (EW) direction, North-South (NS) direction) and one vertical direction (Up-Down (UD) direction). Even though extensive studies have been carried out on the earthquake-triggered landslide, the effects of each single seismic wave and the three-directional seismic waves are not considered. This study aims to evaluate the effects of different types of the seismic waves on the dynamic response and failure behavior of the high-steep rock slide. To investigate the effects of each single seismic wave and three-directional seismic wave, this study presents a numerical model with four types of seismic waves, e.g., East-West (EW) direction, North-South (NS) direction, Up-Down (UD) direction, and three-directional wave (EW_NS_UD). The numerical results revealed that the types of the seismic waves have significantly different effects on the dynamic process, failure behavior, run-out distance, velocity, and deposition of the high-steep rock slide.

Linearized Frequency-Dependent Reflection Coefficient and Attenuated Anisotropic Characteristics of Q-VTI Model

Seismic wave exhibits the characteristics of anisotropy and attenuation while propagating through the fluid-bearing fractured or layered reservoirs, such as fractured carbonate and shale bearing oil or gas. We derive a linearized reflection coefficient that simultaneously considers the effects of anisotropy and attenuation caused by fractures and fluids. Focusing on the low attenuated transversely isotropic medium with a vertical symmetry axis (Q-VTI) medium, we first express the complex stiffness tensors based on the perturbation theory and the linear constant Q model at an arbitrary reference frequency, and then we derive the linearized approximate reflection coefficient of P to P wave. It decouples the P- and S-wave inverse quality factors, and Thomsen-style attenuation-anisotropic parameters from complex P- and S-wave velocity and complex Thomsen anisotropic parameters. By evaluating the reflection coefficients around the solution point of the interface of two models, we analyze the characteristics of reflection coefficient vary with the incident angle and frequency and the effects of different Thomsen anisotropic parameters and attenuation factors. Moreover, we realize the simultaneous inversion of all parameters in the equation using an actual well log as a model. We conclude that the derived reflection coefficient may provide a theoretical tool for the seismic wave forward modeling, and again it can be implemented to predict the reservoir properties of fractures and fluids based on diverse inversion methods of seismic data.

Seismic Wave Field Anomaly Identification of Ultra-Deep Heterogeneous Fractured-Vuggy Reservoirs: A Case Study in Tarim Basin, China

Ultra-deep (7500–9000 m) Ordovician tight limestone heterogeneous fractured-vuggy reservoir is an important target of FuMan Oilfield in Tarim Basin. The strike-slip fault controlled reservoir is related to formation fracture and dissolution caused by geological stress. The seismic wave-field anomaly characteristics with different energy and irregular waveform are displayed in the seismic profile. Accurate identification of fractured-vuggy reservoirs wrapped in tight limestone is the direct scheme to improve production efficiency. Therefore, a new combination method flow of seismic wave-field anomaly recognition is proposed. In this process, the seismic data must be preprocessed initially, and on this basis, robust formation dip scanning is carried out. Secondly, the dip data is applied to the transverse smoothing filter to obtain the formation background data. Eventually, the seismic wave-field anomaly data is the residual between background data and original seismic data. This method has been applied in blocks with different structural characteristics and can effectively improve the resolution of strike-slip fault controlled reservoirs. Based on the results, the drilling success rate is increased to more than 95%, and the high-yield rate of oil tests is increased to 75% in 2021. Multiple applications indicate that the method is robust and can be popularized.

Microzonation of Probolinggo Fault Using Microtremor as an Effort for Earthquake Disaster Mitigation

Earthquake can be caused by several things, one of which is due to an active fault. To mitigate earthquake disasters that can be caused by the Probolinggo Fault, measurement of the microtremor geophysical method is carried out to analyze seismic vulnerability. In this study, the microtremor measurements were carried out on 30 measurement points. The data obtained from measurements are then processed using EasyHVSR software to obtain natural frequency (f0) and natural amplification (A0) values. This value is then used to create a microtremor microzonation map, which is seismic vulnerability index, based on natural frequency and amplification . it founded that low natural frequency mostly founded on eastern of measurement area, caused by thick sedimentary from Lamongan volcanic. High amplification mostly founded from middle to western area, and high seismic vulnerability index founded on western of measurement area, include Maron and Krucil Sub-district. It means seismic wave can very destructive on those area.