地震波干渉法に基づく地震波散乱特性の時空間変化推定のための自己完結的手法の構築

Tohoku University中原恒課

INVERSE ATTENUATION-FILTERING

When seismic waves propagate through the Earth, they are affected by numerous inelastic effects of the medium. These effects are usually characterized by the concept of the Q-factor and lead to variations of spectra of the signal and shapes of the waveforms, which further affect the results of reflection seismic imaging. Attenuation compensation, also often called the inverse Q filtering is a signal-processing procedure broadly used to compensate both of these effects of attenuation in reflection sections or volumes. The objective of this thesis is to present and investigate a new attenuation-compensation approach that is much more general than the conventional inverse Q filtering

Chapitre 2 • Well seismic surveying

Approaches that are typically applied in deep exploration geophysics, combining different seismic and logging methods, can be technically adapted for certain geotechnical or hydrogeological surveys or some site characterizations in the framework of seismic hazard studies. Currently it is entirely feasible to implement this type of geophysical surveying if the situation requires. After reviewing the current state of knowledge regarding borehole measurements of subsurface shear velocities applied to the geotechnical field, this book illustrates the feasibility of carrying out vertical seismic profiles (VSPs) and logs in this field. This approach also illustrates the value of combining velocity measurements of formations provided by borehole seismic tools (VSP) and acoustic (sonic) tools. An innovative example of the application of borehole seismic and logging methods is then presented in the case study of a relatively near-surface (from 20 to 130 m) karst carbonate aquifer. It shows how a multi-scale description of the reservoir can be carried out by integrating the information provided by different 3D-THR surface seismic methods, full waveform acoustic logging, VSP with hydrophones, borehole optical televiewer and flow measurements. In this book the authors provide readers with guidelines to carry out these operations, in terms of acquisitions as well as processing and interpretation. Thus, users will be able to draw inspiration to continue transferring petroleum techniques and other innovative methods for use in near-surface studies

Geophysical data registration using modified plane-wave destruction filters

I propose a method to efficiently measure local shifts, slopes, and scaling functions between seismic traces using modified plane-wave destruction filters. Plane-wave destruction can efficiently measure shifts of less than a few samples, making this algorithm particularly effective for detecting small shifts. When shifts are large, amplitude-adjusted plane-wave destruction can also be used to refine shift estimates obtained by other methods. Amplitude-adjusted plane-wave destruction separates estimation of local shifts and amplitude weights, allowing the time-shift to be measured more accurately. This algorithm has clear applications to geophysical data registration problems, including time-lapse image registration, multicomponent image registration, automatic gather flattening, automatic seismic-well ties, and image merging. The effectiveness of this algorithm in predicting shifts associated with fluid migration, wave mode conversions, and anisotropy and amplitude gradients associated with amplitude variations with offset or angle is demonstrated by applying the algorithm to a synthetic trace, a time-lapse field data example from the Cranfield CO₂ sequestration project, a multicomponent field data example from West Texas, and the Mobil AVO prestack seismic data. Finding correspondence between different parts of the same dataset falls into the same category of problems as local shift estimation. Computation of structure-oriented amplitude gradients for attribute-assisted interpretation requires the estimation of local slopes by correlating reflections between neighboring seismic traces in an image. One of the major challenges of interpreting seismic images is the delineation of reflection discontinuities that are related to geologic features, such as faults, channels, salt boundaries, and unconformities. Visually prominent reflection features often overshadow these subtle discontinuous features which are critical to understanding the structural and depositional environment of the subsurface. For this reason, precise manual interpretation of these reflection discontinuities in seismic images can be tedious and time-consuming, especially when data quality is poor. Discontinuity enhancement attributes are commonly used to facilitate the interpretation process by enhancing edges in seismic images and providing a quantitative measure of the significance of discontinuous features. These attributes require careful pre-processing to maintain geologic features and suppress acquisition and processing artifacts which may be artificially detected as a geologic edge. The plane-wave Sobel filter cascades plane-wave destruction filters with plane-wave shaping in the transverse direction to compute an enhanced discontinuity attribute. The plane-wave Sobel attribute can be applied directly to a seismic image to efficiently and effectively enhance discontinuous features, or to a coherence image to create a sharper and more detailed image. I demonstrate the effectiveness of this method by applying it to two field data sets from offshore New Zealand and offshore Nova Scotia with several faults and channel features and compare the results to other coherence attributes.Geological Science

Scientific Rationale and Requirements for a Global Seismic Network on Mars

Following a brief overview of the mission concepts for a Mars Global Network Mission as of the time of the workshop, we present the principal scientific objectives to be achieved by a Mars seismic network. We review the lessons for extraterrestrial seismology gained from experience to date on the Moon and on Mars. An important unknown on Mars is the expected rate of seismicity, but theoretical expectations and extrapolation from lunar experience both support the view that seismicity rates, wave propagation characteristics, and signal-to-noise ratios are favorable to the collection of a scientifically rich dataset during the multiyear operation of a global seismic experiment. We discuss how particular types of seismic waves will provide the most useful information to address each of the scientific objectives, and this discussion provides the basis for a strategy for station siting. Finally, we define the necessary technical requirements for the seismic stations

Seismic Attenuation Anisotropy in the Southernmost Part of the Taupo Volcanic Zone, North Island, New Zealand

We investigate the mechanisms of seismic anisotropy and attenuation (1/Q) beneath the southernmost part of the Taupo Volcanic Zone (TVZ) by computing variations in S-wave attenuation factors with the direction of wave polarization. We rotate pairs of horizontal components in steps of 22.5◦ from 0◦ to 67.5◦ and into the radial and transverse directions to search for the optimal separation of the attenuation curves and thereby determine an anisotropy symmetry system. The frequency dependence of Q for the rotated S-waves is estimated by means of the non-parametric generalized inversion technique (GIT) of Castro et al. (1990) using shallow earthquakes (< 40 km depth) recorded by GeoNet within 100 km of Mt. Ruapehu. To analyze the effects on computed attenuation properties of source locations, we divide our dataset into two groups: a “TVZ” group containing earthquakes within the TVZ in a distance range of 5–55 km and a “non-TVZ” group containing earthquakes outside the TVZ in a distance range of 5–50 km. To measure Q, we compute the spectral amplitude decay with distance in terms of empirical functions at 20 separate frequencies in the frequency bands 2–10 Hz and 2– 12 Hz for the TVZ and non-TVZ datasets respectively. We construct homogeneous and two-layer Q models for the TVZ dataset based on characteristic features of the attenuation function, while for outside TVZ we only analyse a homogeneous Q model. The homogeneous Q models obtained for the two datasets indicate that S-waves are more attenuated within the TVZ than outside. The homogeneous Q model for the TVZ dataset reveals that the S-wave is anisotropic at high frequencies ( f > 6 Hz) along N–S/E– W directions with the relation QSE ( f ) = (6.15±1.22) f (1.73±0.12) and QSN ( f ) = (4.14± 1.26) f (2.06±0.14), while the non-TVZ dataset shows a weak frequency dependence of attenuation anisotropy at low frequencies in NE–SW/SE–NW directions giving the power law function QSNE ( f ) = (50.93±1.18) f (0.20±0.10) and QSSE ( f ) = (22.60±1.10) f (0.53±0.06). Here, the uncertainty estimates are 95% confidence intervals. To investigate the variation of attenuation anisotropy with depth within the TVZ, we first calculate Q along propagation paths (< 25 km, which corresponds to a maximum turning point depth of 9 km ) and then using paths of 25–55 km length. Small attenuation anisotropy with low attenuation in the N–S direction for the upper crust of TVZ may be related to heterogenous structure as reported by previous studies. Attenuation anisotropy in the northwest direction yielding lower attenuation inferred for the deeper crust suggests the presence of connected melt aligned with the extension direction of TVZ

Seismic imaging of the Cocos plate subduction zone system in central Mexico

Broadband data from the Meso-America Subduction Experiment (MASE) line in central Mexico were used to image the subducted Cocos plate and the overriding continental lithosphere beneath central Mexico using a generalized radon transform based migration. Our images provide insight into the process of subducting relatively young oceanic lithosphere and its complex geometry beneath continental North America. The converted and reverberated phase image shows complete horizontal tectonic underplating of the Cocos oceanic lithosphere beneath the North American continental lithosphere, with a clear image of a very thin low-velocity oceanic crust (7–8 km) which dips at 15–20 degrees at Acapulco then flattens approximately 300 km from the Middle America Trench. Farther inland the slab then appears to abruptly change from nearly horizontal to a steeply dipping geometry of approximately 75 degrees underneath the Trans-Mexican Volcanic Belt (TMVB). Where the slab bends underneath the TMVB, the migrated image depicts the transition from subducted oceanic Moho to continental Moho at ∼230 km from the coast, neither of which were clearly resolved in previous seismic images. The deeper seismic structure beneath the TMVB shows a prominent negative discontinuity (fast-to-slow) at ∼65–75 km within the upper mantle. This feature, which spans horizontally beneath the arc (∼100 km), may delineate the top of a layer of ponded partial mel

Well seismic surveying and acoustic logging

Approaches that are typically applied in deep exploration geophysics, combining different seismic and logging methods, can be technically adapted for certain geotechnical or hydrogeological surveys or some site characterizations in the framework of seismic hazard studies. Currently it is entirely feasible to implement this type of geophysical surveying if the situation requires. After reviewing the current state of knowledge regarding borehole measurements of subsurface shear velocities applied to the geotechnical field, this book illustrates the feasibility of carrying out vertical seismic profiles (VSPs) and logs in this field. This approach also illustrates the value of combining velocity measurements of formations provided by borehole seismic tools (VSP) and acoustic (sonic) tools. An innovative example of the application of borehole seismic and logging methods is then presented in the case study of a relatively near-surface (from 20 to 130 m) karst carbonate aquifer. It shows how a multi-scale description of the reservoir can be carried out by integrating the information provided by different 3D-THR surface seismic methods, full waveform acoustic logging, VSP with hydrophones, borehole optical televiewer and flow measurements. In this book the authors provide readers with guidelines to carry out these operations, in terms of acquisitions as well as processing and interpretation. Thus, users will be able to draw inspiration to continue transferring petroleum techniques and other innovative methods for use in near-surface studies

Multiazimuth velocity analysis using velocity-independent seismic imaging

textMultiazimuth seismic data contains information about how the Earth’s seismic response changes with azimuthal direction. Directional-dependence of the seismic response can be caused by anisotropy or heterogeneity, associated with subsurface features such as fractures, stresses, or structure. Characterizing azimuthal variations is done through velocity analysis, which provides a link between an acquired data set and its image, as well as between the image and subsurface geology. At the stage which conventional velocity analysis is applied, it is difficult to distinguish the geologic cause of observed azimuthal velocity variations. The inability to distinguish the similar effects of anisotropy and heterogeneity leads to positioning errors in the final image and velocity estimates. Regardless of the cause, azimuthally variable velocities require at least three parameters to characterize, as opposed to the conventional single-parameter isotropic velocity. The semblance scan is the conventional tool for seismic velocity analysis, but it was designed for the isotropic case. For multiple parameters, the semblance scan becomes computationally impractical. In order to help address the xiissues of geologic ambiguity and computational efficiency, I develop three methods for multiazimuth seismic velocity analysis based on “velocity-independent” imaging techniques. I call this approach, velocity analysis by velocity-independent imaging, where I reverse the conventional order of velocity estimation followed by image estimation. All three methods measure time-domain effective-velocity parameters. The first method, 3D azimuthally anisotropic velocity-independent NMO, replaces the explicit measurement of velocity with local slope detection. The second method, time-warping, uses local slope information to predict traveltime surfaces without any moveout assumption beforehand, and then fit them with a multiparameter velocity model. The third method, azimuthal velocity continuation, uses diffraction image focusing as a velocity analysis criterion, thereby performing imaging and velocity analysis simultaneously. The first two methods are superior to the semblance scan in terms of computational efficiency and their ability to handle multi-parameter models. The third method is similar to a single multi-parameter semblance scan in computational cost, but it helps handle the ambiguity between structural heterogeneity and anisotropy, which leads to better positioned images and velocity estimates.Geological Science

Rapid magma ascent beneath La Palma revealed by seismic tomography

Data availability The seismic catalogue of IGN is publicly available at: https:// www. ign. es/ web/ ign/ portal/ sis- catal ogo- terre motos. The seismic catalogue of INVOLCAN is available under request to Dr. Luca D’Auria ([email protected]). The LOTOS code is publicly available at: www. ivan- art. com/ scien ce/ LOTOS. An online version of the code with the La Palma dataset is available in: Koulakov Ivan. (2022). Data and program codes to reproduce the results of seismic tomography for La Palma Island [Data set]. Zenodo. https:// doi. org/ 10. 5281/ zenodo. 65893 67. The digital elevation model used in all figures and historical lava flows of Figs. 1 and 3 were downloaded from the public graphic repository of GrafCan (www. grafc an. es). The 2021 lava flow was downloaded from the European agency Copernicus Emergency Management Service (httts://emergency.copernicus.eu/mapping/list-of-components/ EMSR546). The software used to generate Fig. 1, Figs. S1, S2 and S3 was QGIS 3.22 (https:// www. qgis. org). The software used to generate Figs. 3, 4 and 6, Figs. S4, S5 and S6 is the LOTOS code.Acknowledgements JP and JMI were partially supported by the FEMALE project of the Spanish Government (Grant No. PID2019-106260GB-I00). IK was supported by the Russian Science Foundation (Grant No. 20-17-00075). The INVOLCAN team was supported by the projects VOLRISKMAC II (MAC2/3.5b/328), co-financed by the EC Cooperation Transnational Program MAC 2014-2020, and “Cumbre Vieja Emergencia”, financed by the Spanish Ministry of Science and Innovation. English language editing was performed by Tornillo Scientific, UK.Supplementary Information The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-21818-9.For the first time, we obtained high-resolution images of Earth's interior of the La Palma volcanic eruption that occurred in 2021 derived during the eruptive process. We present evidence of a rapid magmatic rise from the base of the oceanic crust under the island to produce an eruption that was active for 85 days. This eruption is interpreted as a very accelerated and energetic process. We used data from 11,349 earthquakes to perform travel-time seismic tomography. We present high-precision earthquake relocations and 3D distributions of P and S-wave velocities highlighting the geometry of magma sources. We identified three distinct structures: (1) a shallow localised region (< 3 km) of hydrothermal alteration; (2) spatially extensive, consolidated, oceanic crust extending to 10 km depth and; (3) a large sub-crustal magma-filled rock volume intrusion extending from 7 to 25 km depth. Our results suggest that this large magma reservoir feeds the La Palma eruption continuously. Prior to eruption onset, magma ascended from 10 km depth to the surface in less than 7 days. In the upper 3 km, melt migration is along the western contact between consolidated oceanic crust and altered hydrothermal material.FEMALE project of the Spanish Government (Grant No. PID2019-106260GB-I00)Russian Science Foundation (Grant No. 20-17-00075)INVOLCAN team was supported by the projects VOLRISKMAC II (MAC2/3.5b/328)EC Cooperation Transnational Program MAC 2014-2020Spanish Ministry of Science and Innovatio