Seismic Field Research Articles

Supervised Deep Learning for Detecting and Locating Passive Seismic Events Recorded with DAS: A Case Study

Exploring shallow mineral resources requires acquiring denser seismic data, for which Distributed Acoustic Sensing is an effective enabler and relevant to mining operations monitoring. Passive seismic can be of interest in characterizing the subsurface; however, dealing with large amounts of data pushes against the limits of existing computational systems and algorithms, especially for continuous monitoring. Hence, more than ever, novel data analysis methods are needed. In this article, we investigate using synthetic seismic data, paired with real noise recordings, as part of a supervised deep-learning neural network methodology to detect and locate induced seismic sources and explore their potential use to reconstruct subsurface properties. Challenges of this methodology were identified and addressed in the context of induced seismicity applications. Data acquisition and modelling were discussed, preparation workflows were implemented, and the method was demonstrated on synthetic data and tested on relevant seismic monitoring field dataset from the Otway CO2 injection site. Conducted tests confirmed the effects of time shifts, signal-to-noise ratios, and geometry mismatches on the performance of trained models. Those promising results showed the method’s applicability and paved the way for potential application to more field data, such as seismic while drilling.

Quasilinear Earthquake Chains in Groups of Seismic Events Baikal Rift System

Identification of a large number of quasi-linear chains of earthquakes in the epicentral seismicity field of the Baikal region and their study showed that among these chains there may be not only chains of “migration” earthquakes, but also chains formed during a random spatio-temporal distribution of earthquakes. In this work, using a statistical analysis of the distribution of distances between epicenters, the possibility of forming chains of “migrations” within groups of seismic events is shown and their distribution is studied. It is noted that chains in groups of earthquakes are distinguished not only by the distribution of distances, but also by the time between them. The formation of chains of earthquakes was established in the areas of the following groups of earthquakes: Busiyingol and 1976 and 1991. in the same area, the South Baikal, Kyakhtinsky, Kichersky earthquakes of the Tompudinsky series, Oldongsinsky and Charudinsky groups. It is shown that these chains were formed during the implementation of these groups. In the identified areas of concentration of chains of grouping seismicity, a connection can be traced between the location and direction of the chains with the strike of fault zones, near-fault cracks and the orientation of nodal planes in earthquake sources.

Seismic Imaging of the Arctic Subsea Permafrost Using a Least-Squares Reverse Time Migration Method

High-resolution seismic imaging allows for the better interpretation of subsurface geological structures. In this study, we employ least-squares reverse time migration (LSRTM) as a seismic imaging method to delineate the subsurface geological structures from the field dataset for understanding the status of Arctic subsea permafrost structures, which is pertinent to global warming issues. The subsea permafrost structures in the Arctic continental shelf, located just below the seafloor at a shallow water depth, have an abnormally high P-wave velocity. These structural conditions create internal multiples and noise in seismic data, making it challenging to perform seismic imaging and construct a seismic P-wave velocity model using conventional methods. LSRTM offers a promising approach by addressing these challenges through linearized inverse problems, aiming to achieve high-resolution, subsurface imaging by optimizing the misfit between the predicted and the observed seismic data. Synthetic experiments, encompassing various subsea permafrost structures and seismic survey configurations, were conducted to investigate the feasibility of LSRTM for imaging the Arctic subsea permafrost from the acquired seismic field dataset, and the possibility of the seismic imaging of the subsea permafrost was confirmed through these synthetic numerical experiments. Furthermore, we applied the LSRTM method to the seismic data acquired in the Canadian Beaufort Sea (CBS) and generated a seismic image depicting the subsea permafrost structures in the Arctic region.

Leveraging Bayesian methods for addressing multi-uncertainty in data-driven seismic liquefaction assessment

When assessing seismic liquefaction potential with data-driven models, addressing the uncertainties of establishing models, interpreting cone penetration tests (CPT) data and decision threshold is crucial for avoiding biased data selection, ameliorating overconfident models, and being flexible to varying practical objectives, especially when the training and testing data are not identically distributed. A workflow characterized by leveraging Bayesian methodology was proposed to address these issues. Employing a Multi-Layer Perceptron (MLP) as the foundational model, this approach was benchmarked against empirical methods and advanced algorithms for its efficacy in simplicity, accuracy, and resistance to overfitting. The analysis revealed that, while MLP models optimized via maximum a posteriori algorithm suffices for straightforward scenarios, Bayesian neural networks showed great potential for preventing overfitting. Additionally, integrating decision thresholds through various evaluative principles offers insights for challenging decisions. Two case studies demonstrate the framework's capacity for nuanced interpretation of in situ data, employing a model committee for a detailed evaluation of liquefaction potential via Monte Carlo simulations and basic statistics. Overall, the proposed step-by-step workflow for analyzing seismic liquefaction incorporates multifold testing and real-world data validation, showing improved robustness against overfitting and greater versatility in addressing practical challenges. This research contributes to the seismic liquefaction assessment field by providing a structured, adaptable methodology for accurate and reliable analysis.

Neotectonic evolution of the Caucasus: recent vertical movements and mechanism of crustal deformation

It is generally recognized that the formation of the fold-and-thrust tectonic structures of mobile belts on continents is associated with crushing and narrowing of the Earth’s crust as a result of collision of lithospheric plates. The deformation of the Caucasian lithosphere in the neotectonic time is generally consistent with these ideas. However, the block differentiation of the Caucasian lithosphere introduces specific features in the directivity of modern vertical and horizontal movements. In this paper, we analyze vertical movements of the Caucasus estimated by means of high-precision leveling over more than a century and consider their spatial correlation with tectonics, seismicity, stress-strain state, and geophysical fields. A clear relationship indicating the deep tectonic nature of the long-term uplifting of the Caucasus crust is revealed. Due to the differentiation of the Arabian plate movement, the territory of the Caucasus is divided into provinces that differ from each other in the pattern of modern movements, in the orientation of faults, and in the stress-strain state. The seismic regime in these provinces also has differences in the number of seismic events and focal mechanisms of the earthquakes. We propose a model of the deformation mechanism of the Greater Caucasus, which takes into account the long-term trend of the Caucasus uplifting in the conditions of general shortening of the Earth’s crust. The results of the analysis are used as a basis for discussion of a probable mechanism of tectonic evolution of the Greater Caucasus in the neotectonic time, which can be used in the assessment of seismic hazard in the North Caucasus.

A new seismic wave input method for canyon sites based on the finite element method-indirect boundary integral equation method

The seismic input is the basis for the seismic safety analysis of dams, but current seismic input methods have not reasonably considered the influence of canyon scattering effects on the dynamic response of dams. In this paper, the total motion field of the canyon is deconstructed into the free field of a uniform half space and the scattering field of the canyon and they are obtained separately. The indirect boundary integral equation method (IBIEM) is used to determine the scattering field of the canyon, which is superimposed with the free field to obtain the seismic ground motion field (total motion field) of the canyon. The accuracy of the total motion field of the canyon is verified through reference solutions. In the finite element analysis of seismic ground motion field in the canyon, the total motion field at the truncation boundary of the canyon foundation is transformed into the equivalent seismic input loads and the interior of the canyon are simulated using the finite element method (FEM). Then, based on the wave input for the free field combined with the viscoelastic artificial boundary, a new wave input method for the total motion field combined with the viscoelastic artificial boundary is established. The seismic wave input method based on the finite element-indirect boundary integral equation method is applied to numerically determine the seismic response of a trapezoidal canyon. The differences in the seismic response of the canyon under free-field input and total motion field input are discussed. The results show that there is a certain deviation between the displacement and waveform obtained from the free-field input and the solutions obtained from the indirect boundary integral equation method. There is a large error for a large angle of incidence, and the errors increase with increasing seismic wave frequency. When the incident angle is 60°, the maximum error in the frequency domain reaches 92.3 %, and the maximum error in the time domain reaches 250 %. The displacement amplitude and waveform of the canyon obtained from the total motion field input showed agreement with the results of the indirect boundary integral equation method. The wave input method established in this paper has high calculation accuracy and reasonably considers the scattering effect of canyons, which provides a basis for more accurate predictions of the seismic response of dams on canyon foundations.

Velocity model-based adapted meshes using optimal transport

In geophysical numerical models using the finite-element method or its variant, the spectral-element method, to solve seismic wave equations, a mesh is employed to discretize the domain. Generating or adapting a mesh to complex geological properties is a challenging task. To tackle this challenge, we develop an r-adaptivity method to generate or adapt a two-dimensional mesh to a seismic velocity field. Our scheme relies on the optimal transport theory to perform vertices relocation, which generates good-shaped meshes and prevents tangled elements. The mesh adaptation can delineate different regions of interest, like sharp interfaces, salt bodies, and discontinuities. The algorithm has a few user-defined parameters that control the mesh density. With typical seismic velocity examples (e.g., Camembert, SEAM Phase, Marmousi-2), mesh adaptation capability is illustrated within meshes with triangular and quadrilateral elements, commonly employed in seismic codes. Besides its potential use in mesh generation, the method developed can be embedded in seismic inversion workflows like multiscale full waveform inversion to adapt the mesh to the field being inverted without incurring the I/O cost of re-meshing and load rebalancing in parallel computations. The method can be extended to three-dimensional meshes.

Post-earthquake rapid seismic demand estimation at unmonitored locations via Bayesian networks

Post-earthquake safety assessment of buildings and infrastructure poses significant challenges, often relying on time-consuming visual inspections. To expedite this process, safety criteria based on a demand-capacity model are utilized. However, rapid assessment frameworks require accurate estimations of intensity measures (IMs) to estimate seismic demand and assess structural health. Unfortunately, post-earthquake IM values are typically only available at monitored locations equipped with sensors or monitoring systems, limiting broader assessments. Simple spatial interpolation methods, while possible, struggle to consider crucial physical factors such as earthquake magnitude, epicentral distance, and soil type, leading to substantial estimation errors, especially in areas with insufficient or non-uniform seismic station coverage. To address these issues, a novel framework, BN-GMPE, combining a Bayesian network (BN) and a ground motion prediction equation (GMPE), is proposed. BN-GMPE enables inference and prediction under uncertainty, incorporating physical parameters in seismic wave propagation. A further novelty introduced in this work regards separating the near and far seismic fields in the updating process to attain a clearer understanding of uncertainty and more accurate IM estimation. In the proposed approach, a GMPE is employed for the estimation, and the bias and standard deviation of the prediction error are updated after any new information is entered into the network. The proposed method is benchmarked against a classic Kriging interpolator technique, considering some recent earthquake shocks in Italy. The proposed BN framework can naturally extend for estimating the probability of failure of various structures in a targeted region, which represents the ultimate aim of this research.

Retrogradation of carbonate platforms on a rifted margin: the late Ediacaran record of the northwestern Yangtze Craton (SW China)

The Late Ediacaran Dengying Formation, located in the Sichuan Basin of the northwestern Yangtze Craton, is of significant interest in oil and gas exploration due to its abundant pores and vugs within microbial mound-shoal complexes. However, there is still uncertainty regarding the spatiotemporal distribution and controlling factors of the platform margin. This study comprehensively analyzes the retrogradation pattern of the Dengying Formation platform margin using seismic data, well logs, field outcrops, and petrological characteristics. Our findings reveal that the Dengying Formation strata surrounding the rift basin at the northwestern of the Yangtze Craton can be divided into three main depositional facies: basin facies, slope facies, and platform margin facies. Additionally, based on the integration of lithological, log, and seismic characteristics, the Dengying Formation is subdivided into four third-order sequences, with five sequence boundaries and three seismic facies identified. Supported by sequence stratigraphy and geophysical data, we have reconstructed the tectono-sedimentary evolution of the multiple platform margins on the eastern side of the Deyang-Anyue rift in the Sichuan Basin during the late Ediacaran. Our findings indicate that the platform underwent two phases of retrogradation. The second-stage platform margin underwent retrogradation towards the interior, spanning a distance between 10 and 80 km, based on the initial configuration established by the first-stage platform margin. The main controls for progradation and retrogradation of carbonate platforms are eustatic sea-level changes and tectonic activity. Eustatic sea-level changes can be divided into constructive and destructive phases. Constructive phases are commonly observed in highstand systems tracts, while destructive phases are often associated with transgressive systems tracts and are related to platform retrogradation processes. However, sea-level changes alone cannot fully control the process of platform retrogradation. The thermal subsidence following mantle plume events likely played a significant role in the retrogradation of the platform in the study area. During this period, tectonic processes controlled the geometry of the platform and the deposition of carbonates in the platform margin-slope-basin environment. Additionally, karst-related mound-shoal complexes developed extensively along the platform margin of the Dengying Formation in the northwestern Yangtze Craton. The Lower Cambrian dark shales represent high-quality hydrocarbon source rocks, while the Dengying Formation exhibits an optimal source-reservoir configuration.

Performance-based multi-hazard engineering (PB-MH-E): The case of steel buildings under earthquake and wind

Focusing on natural hazards as wind and earthquake, the goal of this work is to understand if interference occurs in the design choices for steel buildings of a certain height, thus determining which design requirement and hazard is predominant, between Ultimate Limit State (ULS) requirements under earthquakes, or comfort of the occupants and the Serviceability Limit State (SLS) of the structure under wind. Thus, in order to compare structural responses, both in linear (i.e., SLSs) and non-linear field (i.e., ULS) a simplified analysis procedure that could be also implemented in Standards and in design practice (so-called SAC-FEMA method, originally introduced in the seismic field by Cornell in the early 2000s and more recently extended to the wind) is adapted, leading to a true optimal Performance-Based Multi-Hazard Design (PB-MH-D) of the structure considering the two hazards. The procedure is applied to two case-study steel buildings, 17 and 60 floor heigh; the approach is shown to efficiently lead to a design solution who consistently tackles with the same level of reliability with the two hazards.

Active Tectonic Deformation of the Qilian Shan, Northeastern Tibetan Plateau

Abstract —The Qilian Shan (or Qilian Mountains), located on the northeastern margin of the Tibetan Plateau, is an actively growing orogenic belt resulting from the far-field impact of the India–Eurasia collision. The northward penetration of the Indian Plate is responsible for intense crustal shortening in the Qilian Shan. However, the tectonic deformation pattern in response to the crustal shortening remains unclear. In this study, we present the regional seismicity, fault activity, and GPS crustal movement velocity field to characterize the active tectonic deformation of the Qilian Shan based on historical data over the past two decades. The results suggest that the western Qilian Shan is characterized by distributed north–south crustal shortening, while the eastern Qilian Shan is dominated by blocklike eastward extrusion of crust along major strike-slip faults coupled with clockwise rotation. North–south crustal shortening and east–west lateral extrusion, two deformation modes responding to the India–Eurasia convergence, match the crustal deformation in the Qilian Shan. The tectonic deformation of the western Qilian Shan is largely in agreement with the former, while the eastern Qilian Shan corresponds closely to the latter. Lower crustal flow beneath the central Tibetan Plateau provides the potential driving force to induce the eastward extrusion of crustal material out of the plateau and the growth of some boundary mountain ranges, such as the Qilian Shan.

Seismic and Potential Field Constraints on Upper Crustal Architecture of Inner Bering Shelf, Offshore Southwestern Alaska

Southwestern Alaska encompasses a group of fault-bounded tectonostratigraphic terranes that were accreted to North America during the Mesozoic and Paleogene. To characterize the offshore extension of these terranes and several significant faults identified onshore, we reprocessed three intersecting multichannel deep seismic reflection profiles totaling ∼750 line-km that were shot by the R/V Ewing across part of the inner Bering continental shelf in 1994. Since the uppermost seismic section is often contaminated by high amplitude water layer multiples from the hard and shallow seafloor, the migrated reflection images are supplemented with high-resolution P wave velocity models derived by traveltime tomography of the recorded first-arrivals to depths of up to 2000 m. Additionally, other geophysical datasets such as seismicity, well logs, high resolution satellite-altimetry gravity, air-borne magnetics, ship-board gravity and magnetics, are also incorporated into an integrated regional interpretation. We delineate the offshore extension of the major mapped geological elements, including the Togiak fault (TGF), East Kulukak fault (EKF), Chilchitna fault (CF), Lake Clark fault (LCF), Togiak terrane (TT), Goodnews terrane (GT), Peninsular terrane (PT), Northern Kahiltna flysch (NKF) and Southern Kahiltna flysch (SKF) deposits, and the regional suture zone. The geophysical evidence from this study suggest that the major faults and terrane boundaries of southwestern Alaska, excluding the LCF, not only extend beneath the Bering shelf offshore but also appear to rotate, forming a trend parallel to the inactive Beringian margin. The LCF extends offshore but likely terminates in the northeastern part of the Bristol Bay basin. Additionally, the constraints from seismicity data indicates that while the major faults in southwestern Alaska exihibit some activity onshore, they remain dormant in the offshore region. These new findings will contribute to a better understanding of terrane accretion processes and existing fault models of southwestern Alaska.

VSP wavefields separating method based on parallel local slope scanning

Due to the detectors being distributed in the target medium, vertical seismic profile (VSP) is a seismic observation method that has the advantages of high resolution and a high signal-to-noise ratio. Separating the mixed wavefields into upgoing and downgoing waves can obtain more obvious dynamic and kinematic characteristics of seismic waves, guiding subsequent imaging and interpretation. The traditional method is mainly based on the pickup of first breaks. Flattening the global seismic data by first breaks can enhance the seismic events through techniques like median filter and singular value decomposition (SVD). However, this method relies on high-precision pickup of first breaks and is limited by zero offset. To address this limitation, we introduce an improved median filtering separation method. This method employs separations of the local dip angles into positive and negative through multi-window scanning (MWS). Due to the high accuracy and robustness of this method in 2-D VSP data, we propose to use the local dip angle of the wavefields to median filter the wavefield through positive and negative angles to obtain upgoing and downgoing waves. This method for wavefield separation is optimized by iteratively identifying the directions of the seismic data. However, this optimization comes at the cost of increased computational requirements, especially under high-precision conditions. To help alleviate this problem, we use multi-thread parallel computing technique on a multi-core central processing unit (CPU) to improve computational efficiency. Finally, we validate the proposed method by testing it on synthetic seismic data and field VSP data, respectively. The results show that this wavefield separation method has advantages in terms of accuracy and robustness compared to the median filtering method based on the first break picking.

Development of a novel seismic acquisition system based on fully autonomous ocean bottom nodes

It is now established that ocean bottom node (OBN) surveys provide superior seismic data compared to traditional towed streamers, but deployment of the nodes is time-consuming and expensive. As an alternative, we are developing a novel seismic acquisition system that utilises a fleet of fully autonomous underwater vehicles (AUVs) as OBNs. The objective is to eliminate the reliance on remotely operated vehicles (ROVs) or ropes for deployment, making ocean bottom seismic acquisition significantly more efficient and cost-effective. To succeed, several criteria need to be met, including recording of good geophysical data, system affordability, low power consumption, units’ ability to act autonomously (e.g. steer, land, listen, re-position, surface) and to respond to remote commands (e.g. cease activity, return), a failure recovery protocol, automated deployment and retrieval, and system-wide communications. The ability to record good seismic data was assessed in October 2022, when we performed active-source seismic field trials in the UK. Five AUVs were deployed adjacent to commercial nodes (all units carried identical geophysical payloads) and a series of 2D lines were acquired at varying orientations. Analysis of the data confirmed that our units achieved satisfactory coupling with the seabed. The unit’s in-water capabilities were assessed during multiple field trials in 2022 and 2023, offshore UK and Australia. Tests consisted of simulating seismic acquisition cycle (flights, landings, take-offs and repositioning) where it was verified that missions were reliably executed, including landing within 10 m from pre-planned locations and recording good quality passive seismic data.

Rapid determination of seismic influence field based on mobile communication big data-A case study of the Luding Ms 6.8 earthquake in Sichuan, China.

Smartphone location data provide the most direct field disaster distribution data with low cost and high coverage. The large-scale continuous sampling of mobile device location data provides a new way to estimate the distribution of disasters with high temporal-spatial resolution. On September 5, 2022, a magnitude 6.8 earthquake struck Luding County, Sichuan Province, China. We quantitatively analyzed the Ms 6.8 earthquake from both temporal and geographic dimensions by combining 1,806,100 smartphone location records and 4,856 spatial grid locations collected through communication big data with the smartphone data under 24-hour continuous positioning. In this study, the deviation of multidimensional mobile terminal location data is estimated, and a methodology to estimate the distribution of out-of-service communication base stations in the disaster area by excluding micro error data users is explored. Finally, the mathematical relationship between the seismic intensity and the corresponding out-of-service rate of communication base stations is established, which provides a new technical concept and means for the rapid assessment of post-earthquake disaster distribution.

Spectral analysis of seismic data based on a combined generalized S transform

The S transform is widely used to delineate seismic spectral anomalies associated with hydrocarbons during reservoir characterization. To overcome the relatively poor time-frequency resolution due to the fixed varying pattern of the window used in the S transform, we present a combined generalized S transform (CGST) for the spectral analysis of seismic data. First, the CGST decomposes the signal into a series of intrinsic mode functions (IMFs) using a complete ensemble empirical mode decomposition (CEEMD). An objective function is defined to quantitatively determine the spectral energy concentration criteria for optimal parameter selection of the improved Gaussian windows. Then, the time-frequency spectra for all the IMFs within narrower frequency bands are calculated using the respective Gaussian windows with normalized energy and individual optimal controlling parameters. We applied the CGST to spectral analyses of both synthetic signals and seismic field data. The results indicate that the CGST exhibits good performance with a higher time-frequency resolution. Field examples illustrate that the CGST can be used to characterize the spectral behavior of thin beds and indicate low-frequency shadows associated with hydrocarbon reservoirs.

Strain field reconstruction from helical-winding fiber distributed acoustic sensing and its application in anisotropic elastic reverse time migration

Optical fiber-based distributed acoustic sensing (DAS) technology has been a popular seismic acquisition tool due to its easy deployment, wide bandwidth, and dense sampling. However, the sensitivity of straight optical fiber to only single-axis strain presents challenges in fully characterizing multicomponent seismic wavefields, making it difficult to use these data in elastic reverse time migration (ERTM). The helical-winding fiber receives projecting signals projected onto the fiber from all seismic strain field components and has the potential to reconstruct those strain components for ERTM imaging. Here, we give detailed mathematical principles of helical fiber-based DAS with crucial parameters such as pitch angle, gauge length, and rotating angle. At least six points of DAS responses are required in one or several winding periods to rebuild the strain fields within the seismic wavelength. The projecting matrix of conventional regular helical-winding fiber is singular and ill conditioned, which results in computation challenges for the inverse of the Hessian matrix for strain component reconstruction. To tackle this problem, we develop a nonregular variant pitch-angle winding configuration for helical fiber. Our winding design is validated using the rank and condition number of the projecting matrix, which is proven to be an important tool in the reconstruction of the original seismic strains. The recovered strain components from the DAS response are then used to backward propagate the receiver wavefields in ERTM with an efficient P/S decoupled approach. To summarize, we develop a novel winding design of helical fiber to recover the strain fields and then develop an efficient 3D anisotropic P/S wave-mode decomposition method for generating vector P and S wavefields during their propagation. Both methods are applied to build an anisotropic DAS-ERTM workflow for producing PP and PS images. Two synthetic examples demonstrate the effectiveness of our approach.

An Analytical Review of the Recent Crustal Uplifts, Tectonics, and Seismicity of the Caucasus Region

This paper analyzes and reviews the rapid uplifts of the Earth’s crust in the Caucasus that occurred over the last century. The uplifts were registered by precise repeated state leveling and reflected on officially published maps of vertical movements of the Earth’s crust. This study summarizes information on the region’s vertical movements over more than a century. The present study describes the technology for creating maps of recent vertical movements of the Earth’s crust using precision leveling data. This paper summarizes cases of recording uplifts of the Earth’s surface in other regions of the world in connection with seismic activity. The authors carried out intercomparison of vertical movements with tectonics, seismicity, and geophysical fields, which discovered their apparent mutual correspondence. This indicates the deep tectonic nature of the observed uplifts of the Earth’s crust. Spatial and temporal agreement with the distribution of strong earthquakes showed a natural relationship. It has been shown that strong earthquakes are confined to the boundaries of zones of rapid uplift. They occur predominantly in areas of transition between uplifts and subsidence. The results obtained demonstrate the role of the study and observations of vertical movements of the Caucasus in assessing periods and areas of increased seismic hazard.

3D seismic Fault Detection via Contrastive-Reconstruction Representation Learning

Fault detection is a critical step in structural modeling and characterizing reservoirs. 3D seismic fault labeling is almost impossible to obtain, while networks trained on synthetic fault data have limited generalization to real data. We use self-supervised representation learning to enable networks to learn seismic field data features during pre-training, improving generalization. However, widely popular ViT/SwinViT-based methods cannot capture low-level fault features well. CNN-based have limited capability in transferring to downstream tasks. Tackle this, we designed the Tiny Self-Attention and extensively embedded it into the HRNet. It merges the advantages of convolutional and self-attention, allowing for in-depth learning of seismic representation information during the representation learning phase, thus achieving better inter-class separation in downstream tasks. We crafted two proxy tasks. The contrastive task aims to minimize the projected feature distance of overlapping areas from two distinct views. To tackle the issue of memory overflow and training suspension caused by the high spatial and temporal complexity of 3D high-resolution feature matching, we introduced sparse distance matching and an adaptive feature aggregation. In reconstruction task, we designed a masking strategy tailored to the unique attributes of 3D seismic data. It process named “Fault Detection via Contrast-Reconstruction Representation Learning” (FaultCRL). Experimentally, we compared the latest fault detection methods, and representation learning techniques like MAE and SwinUNETR, showcasing the notable advantage of FaultCRL in fault detection tasks. The appendix provides all qualitative results from mainstream geophysical journals regarding The Netherlands-F3 data in recent years, indicating that our approach has reached the current state-of-the-art level.

Seismo-acoustic anomalous interference pattern in shallow water

In shallow water layered leaky waveguide, the sound energy below the cut-off frequency is coupled to the seismic sound field. Very low frequency (VLF) sound is propagated in both water-borne and bottom-trapped modes. Under these circumstances, the interference characteristics of sound field and its application are different from those of medium and high frequency, which is the subject of this study. In a multi-channel seismic exploration experiment in the South Yellow Sea, an anomalous interference pattern for the VLF band, e.g., the waveguide invariant of interference striations changes from positive to negative with increase of frequency, was observed in the range-frequency domain. Numerical simulations and modal analyses show that the anomalous interference pattern corresponds to the interference between waterborne modes and bottom-trapped modes, because the group velocity for these two types normal modes change differently in the VLF band.