Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (588)

Search Parameters:
Keywords = the seismic network

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
28 pages, 9654 KiB  
Article
Time Series Foundation Models and Deep Learning Architectures for Earthquake Temporal and Spatial Nowcasting
by Alireza Jafari, Geoffrey Fox, John B. Rundle, Andrea Donnellan and Lisa Grant Ludwig
GeoHazards 2024, 5(4), 1247-1274; https://doi.org/10.3390/geohazards5040059 - 21 Nov 2024
Viewed by 103
Abstract
Advancing the capabilities of earthquake nowcasting, the real-time forecasting of seismic activities, remains crucial for reducing casualties. This multifaceted challenge has recently gained attention within the deep learning domain, facilitated by the availability of extensive earthquake datasets. Despite significant advancements, the existing literature [...] Read more.
Advancing the capabilities of earthquake nowcasting, the real-time forecasting of seismic activities, remains crucial for reducing casualties. This multifaceted challenge has recently gained attention within the deep learning domain, facilitated by the availability of extensive earthquake datasets. Despite significant advancements, the existing literature on earthquake nowcasting lacks comprehensive evaluations of pre-trained foundation models and modern deep learning architectures; each focuses on a different aspect of data, such as spatial relationships, temporal patterns, and multi-scale dependencies. This paper addresses the mentioned gap by analyzing different architectures and introducing two innovative approaches called Multi Foundation Quake and GNNCoder. We formulate earthquake nowcasting as a time series forecasting problem for the next 14 days within 0.1-degree spatial bins in Southern California. Earthquake time series are generated using the logarithm energy released by quakes, spanning 1986 to 2024. Our comprehensive evaluations demonstrate that our introduced models outperform other custom architectures by effectively capturing temporal-spatial relationships inherent in seismic data. The performance of existing foundation models varies significantly based on the pre-training datasets, emphasizing the need for careful dataset selection. However, we introduce a novel method, Multi Foundation Quake, that achieves the best overall performance by combining a bespoke pattern with Foundation model results handled as auxiliary streams. Full article
Show Figures

Figure 1

Figure 1
<p>Illustration of the construction of a nowcast model for California. The nowcast is a 2-parameter filter on the small earthquake seismicity [<a href="#B42-geohazards-05-00059" class="html-bibr">42</a>,<a href="#B43-geohazards-05-00059" class="html-bibr">43</a>]. (<b>a</b>) Seismicity in the Los Angeles region since 1960, M &gt; 3.29. (<b>b</b>) Monthly rate of small earthquakes as cyan vertical bars. The blue curve is the 36-month exponential moving average (EMA). (<b>c</b>) Mean rate of small earthquakes since 1970. (<b>d</b>) Nowcast curve that is the result of applying the optimized EMA and corrections for the time-varying small earthquake rate to the small earthquake seismicity. (<b>e</b>) Optimized receiver operating characteristic (ROC) curve (red line) used in the machine learning algorithm. Skill is the area under the ROC curve and is used in the optimization. Skill trade-off diagram shows the range of models used in the optimization.</p>
Full article ">Figure 2
<p>Image showing the application of the trained QuakeGPT transformer to an independent, scaled nowcast validation curve (green shading), followed by prediction of future values beyond the end of the nowcast curve (magenta shading). In this model, 36 previous values are used to predict the next value. Dots show the predictions and the solid line shows the nowcast curve whose values are to be predicted. Green dots show the predictions of the transformer up to the last 37 values. The 36 blue dots are predictions that were made and then fed back into the transformer to predict the final point (red dot). In this model, 50 members of an ensemble of runs were used to make the predictions. The dots represent the mean predictions. Brown areas represent the 1-sigma standard deviations to the mean values. In this model, 2021 years of simulation data were used to train the model.</p>
Full article ">Figure 3
<p>Distribution of earthquake epicenters in Southern California (32° N to 36° N, −120° to −114°) from USGS data (1986–2024). The scatter plot shows the spatial density of seismic events used to analyze and optimize spatial bins for earthquake nowcasting.</p>
Full article ">Figure 4
<p>The 500 most active and vulnerable spatial bins, marked in blue, selected for analysis out of the total 2400, based on the frequency of earthquakes from 1986 to 2024. This selection focuses on high-risk areas.</p>
Full article ">Figure 5
<p>Six time series from randomly selected spatial bins, highlighting earthquakes of magnitude greater than 5.</p>
Full article ">Figure 6
<p>The final graph structure representing the 500 most active bins, created using an epsilon of 0.15 degrees. Initially forming a multi-component graph, components are linked to ensure full connectivity.</p>
Full article ">Figure 7
<p>Released energy time series plots for six randomly selected spatial bins, comparing model predictions (GNNCoder one-layer, DilatedRNN, TiDE, iTransformer-M4) against actual observed seismic activities. The brown line represents our GNN approach, which shows a closer match with the actual time series, capturing crucial upward slopes that may signal an impending earthquake. The green and red lines occasionally miss these trends, making more errors where even slight changes in seismic activity are critical. The purple line from the iTransformer-M4 model fails to accurately capture the time series values and exhibits excessive fluctuations.</p>
Full article ">Figure 8
<p>This plot illustrates the spatial bins overlaid on the fault lines to assess the extent to which the fault lines are captured by the bins (graph nodes). It highlights the limitations of the current graph, where some critical fault lines fall outside the spatial bins, impacting the performance of deeper GNN models like the GNNCoder 3-layer model.</p>
Full article ">
22 pages, 10038 KiB  
Article
Analytical Fragility Surfaces and Global Sensitivity Analysis of Buried Operating Steel Pipeline Under Seismic Loading
by Gersena Banushi
Appl. Sci. 2024, 14(22), 10735; https://doi.org/10.3390/app142210735 - 20 Nov 2024
Viewed by 199
Abstract
The structural integrity of buried pipelines is threatened by the effects of Permanent Ground Deformation (PGD), resulting from seismic-induced landslides and lateral spreading due to liquefaction, requiring accurate analysis of the system performance. Analytical fragility functions allow us to estimate the likelihood of [...] Read more.
The structural integrity of buried pipelines is threatened by the effects of Permanent Ground Deformation (PGD), resulting from seismic-induced landslides and lateral spreading due to liquefaction, requiring accurate analysis of the system performance. Analytical fragility functions allow us to estimate the likelihood of seismic damage along the pipeline, supporting design engineers and network operators in prioritizing resource allocation for mitigative or remedial measures in spatially distributed lifeline systems. To efficiently and accurately evaluate the seismic fragility of a buried operating steel pipeline under longitudinal PGD, this study develops a new analytical model, accounting for the asymmetric pipeline behavior in tension and compression under varying operational loads. This validated model is further implemented within a fragility function calculation framework based on the Monte Carlo Simulation (MCS), allowing us to efficiently assess the probability of the pipeline exceeding the performance limit states, conditioned to the PGD demand. The evaluated fragility surfaces showed that the probability of the pipeline exceeding the performance criteria increases for larger soil displacements and lengths, as well as cover depths, because of the greater mobilized soil reaction counteracting the pipeline deformation. The performed Global Sensitivity Analysis (GSA) highlighted the influence of the PGD and soil–pipeline interaction parameters, as well as the effect of the service loads on structural performance, requiring proper consideration in pipeline system modeling and design. Overall, the proposed analytical fragility function calculation framework provides a useful methodology for effectively assessing the performance of operating pipelines under longitudinal PGD, quantifying the effect of the uncertain parameters impacting system response. Full article
(This article belongs to the Section Civil Engineering)
Show Figures

Figure 1

Figure 1
<p>Pipeline subjected to longitudinal PGD: (<b>a</b>) 3D view; (<b>b</b>) 2D schematic representation.</p>
Full article ">Figure 2
<p>Pipeline response to longitudinal PGD according to analytical model in [<a href="#B11-applsci-14-10735" class="html-bibr">11</a>], assuming symmetric material behavior for tension and compression: (<b>a</b>) case I; (<b>b</b>) case II.</p>
Full article ">Figure 3
<p>Schematic representation of operating pipeline response subjected to longitudinal PGD: (<b>a</b>) pipeline displacement subjected to longitudinal soil block movement (case II); (<b>b</b>) soil–pipeline system behaving like a pull-out test under tension (region I) and compression (region IV).</p>
Full article ">Figure 4
<p>Schematic representation of the axial constitutive behavior of the steel pipe material, defined within the associated von Mises plasticity with isotropic hardening [<a href="#B30-applsci-14-10735" class="html-bibr">30</a>].</p>
Full article ">Figure 5
<p>The comparison between the numerical, the conventional [<a href="#B8-applsci-14-10735" class="html-bibr">8</a>,<a href="#B11-applsci-14-10735" class="html-bibr">11</a>,<a href="#B13-applsci-14-10735" class="html-bibr">13</a>], and the proposed analytical models, evaluating the pipeline performance under longitudinal PGD (<span class="html-italic">L<sub>b</sub></span> = 300 m) in terms of maximum tensile and compressive pipe strain as a function of the ground displacement <span class="html-italic">δ</span>.</p>
Full article ">Figure 6
<p>The variation of the critical soil block length, <span class="html-italic">L<sub>cr</sub></span> = (<span class="html-italic">F<sub>t,max</sub></span> − <span class="html-italic">F<sub>c,max</sub></span>)/<span class="html-italic">f<sub>s</sub></span>, as a function of the ground displacement <span class="html-italic">δ</span>, with an indication of the critical values (<span class="html-italic">δ<sub>cr,i</sub></span>, <span class="html-italic">L<sub>cr,i</sub></span>) associated with the achievement of the pipeline performance limit states.</p>
Full article ">Figure 7
<p>The peak axial strain magnitude in the pressurized pipeline (<span class="html-italic">P<sub>i</sub></span>/<span class="html-italic">P<sub>max</sub></span> = 0.75, Δ<span class="html-italic">T</span> = 50 °C) as a function of the PGD length <span class="html-italic">L<sub>b</sub></span> and displacement <span class="html-italic">δ</span> for (<b>a</b>) tension and (<b>b</b>) compression. The dashed horizontal curves represent the strain isolines corresponding to the NOL and PIL performance limit states.</p>
Full article ">Figure 8
<p>The peak axial strain magnitude in the unpressurized pipeline (<span class="html-italic">P<sub>i</sub></span>/<span class="html-italic">P<sub>max</sub></span> = 0, Δ<span class="html-italic">T</span> = 0 °C) as a function of the PGD length <span class="html-italic">L<sub>b</sub></span> and displacement <span class="html-italic">δ</span> for (<b>a</b>) tension and (<b>b</b>) compression. The dashed horizontal curves represent the strain isolines corresponding to the NOL and PIL performance limit states.</p>
Full article ">Figure 9
<p>Fragility surface of buried pipeline (<span class="html-italic">H<sub>c</sub></span> = 1.5 m) for (<b>a</b>) Normal Operability Limit (NOL) and (<b>b</b>) Pressure Integrity Limit (PIL).</p>
Full article ">Figure 10
<p>Schematic representation of the performance assessment of the buried pipeline subjected to the PGD demand (<span class="html-italic">δ</span>, <span class="html-italic">L<sub>b</sub></span>), using the deterministic and fragility analysis framework.</p>
Full article ">Figure 11
<p>Fragility surface of buried pipeline for different cover depths and performance limit states: (<b>a</b>) <span class="html-italic">H<sub>c</sub></span> = 1.0 m, NOL; (<b>b</b>) <span class="html-italic">H<sub>c</sub></span> = 1.0 m, PIL and (<b>c</b>) <span class="html-italic">H<sub>c</sub></span> = 2.0 m, NOL; and (<b>d</b>) <span class="html-italic">H<sub>c</sub></span> = 2.0 m, PIL.</p>
Full article ">Figure 12
<p>The comparison of the first-order and total-order sensitivity indices of the system input parameters for the (<b>a</b>) NOL and (<b>b</b>) PIL performance limit states.</p>
Full article ">Figure A1
<p>Response of the pressurized pipeline (<span class="html-italic">P<sub>i</sub></span>/<span class="html-italic">P<sub>max</sub></span> = 0.75, Δ<span class="html-italic">T</span> = 50 °C) to longitudinal PGD with block length <span class="html-italic">L<sub>b</sub></span> = 200 m (case I): (<b>a</b>) pipe axial force; (<b>b</b>) pipe axial stress; (<b>c</b>) soil friction; (<b>d</b>) ground displacement; (<b>e</b>) pipe axial displacement; (<b>f</b>) pipe axial strain vs. distance from tension crack.</p>
Full article ">Figure A2
<p>Response of the unpressurized pipeline (<span class="html-italic">P<sub>i</sub></span>/<span class="html-italic">P<sub>max</sub></span> = 0, Δ<span class="html-italic">T</span> = 0 °C) to longitudinal PGD with block length <span class="html-italic">L<sub>b</sub></span> = 200 m (case I): (<b>a</b>) pipe axial force; (<b>b</b>) pipe axial stress; (<b>c</b>) soil friction; (<b>d</b>) ground displacement; (<b>e</b>) pipe axial displacement; (<b>f</b>) pipe axial strain vs. distance from tension crack.</p>
Full article ">Figure A3
<p>Response of the pressurized pipeline (<span class="html-italic">P<sub>i</sub></span>/<span class="html-italic">P<sub>max</sub></span> = 0.75, Δ<span class="html-italic">T</span> = 50 °C) to longitudinal PGD with block length <span class="html-italic">L<sub>b</sub></span> = 300 m (case II): (<b>a</b>) pipe axial force; (<b>b</b>) pipe axial stress; (<b>c</b>) soil friction; (<b>d</b>) ground displacement; (<b>e</b>) pipe axial displacement; (<b>f</b>) pipe axial strain vs. distance from tension crack.</p>
Full article ">Figure A4
<p>Response of the pressurized pipeline (<span class="html-italic">P<sub>i</sub></span>/<span class="html-italic">P<sub>max</sub></span> = 0, Δ<span class="html-italic">T</span> = 0 °C) to longitudinal PGD with block length <span class="html-italic">L<sub>b</sub></span> = 300 m (case II): (<b>a</b>) pipe axial force; (<b>b</b>) pipe axial stress; (<b>c</b>) soil friction; (<b>d</b>) ground displacement; (<b>e</b>) pipe axial displacement; (<b>f</b>) pipe axial strain vs. distance from tension crack.</p>
Full article ">
12 pages, 9300 KiB  
Article
Field Experiments of Distributed Acoustic Sensing Measurements
by Haiyan Shang, Lin Zhang and Shaoyi Chen
Photonics 2024, 11(11), 1083; https://doi.org/10.3390/photonics11111083 - 18 Nov 2024
Viewed by 288
Abstract
Modern, large bridges and tunnels represent important nodes in transportation arteries and have a significant impact on the development of transportation. The health and safety monitoring of these structures has always been a significant concern and is reliant on various types of sensors. [...] Read more.
Modern, large bridges and tunnels represent important nodes in transportation arteries and have a significant impact on the development of transportation. The health and safety monitoring of these structures has always been a significant concern and is reliant on various types of sensors. Distributed acoustic sensing (DAS) with telecommunication fibers is an emerging technology in the research areas of sensing and communication. DAS provides an effective and low-cost approach for the detection of various resources and seismic activities. In this study, field experiments are elucidated, using DAS for the Hong Kong–Zhuhai–Macao Bridge, and for studying vehicle trajectories, earthquakes, and other activities. The basic signal-processing methods of filtering and normalization are adopted for analyzing the data obtained with DAS. With the proposed DAS technology, the activities on shore, vehicle trajectories on bridges and in tunnels during both day and night, and microseisms within 200 km were successfully detected. Enabled by DAS technology and mass fiber networks, more studies on sensing and communication systems for the monitoring of bridge and tunnel engineering are expected to provide future insights. Full article
(This article belongs to the Section Lasers, Light Sources and Sensors)
Show Figures

Figure 1

Figure 1
<p>The map for the DAS test with the optical fiber cables along the Hong Kong–Zhuhai–Macao Bridge in Guangdong–Hong Kong–Macau Greater Bay Area, China. The inset shows an example of the fiber optic cabling along the bridge corridor.</p>
Full article ">Figure 2
<p>The results of the original wave swarm within 21 min after UTC 2021-08-06 08:10:06, in the sea area near the Hong Kong–Zhuhai–Macao Bridge Port. The horizontal axis in the figure represents the measurement time with unit seconds (s). The vertical axis is along the fiber direction from Zhuhai Port (<b>bottom</b>) to Hong Kong Port (<b>top</b>). The unit of fiber length is meters (m). Gray dashed lines mark the sections of the coastal region, bridge, and ocean within the figure.</p>
Full article ">Figure 3
<p>The results of the strain-rate data in the first section of 2.8 km. (<b>a</b>) The original wave waterfall swarm plot, which marks the 0–500 m section with a red dashed box; (<b>b</b>) the calculated f-k spectrum; (<b>c</b>) the zoomed-in view of the 0–500 m range corresponding to the red dashed box in (<b>a</b>).</p>
Full article ">Figure 4
<p>The examples of the original channel wave-time and the spectrum-frequency plot at fiber distances of 210 m, 270 m, and 280 m. The vertical axis represents the strength. The horizontal axis in wave-time plot and spectrum-frequency plot represents time and frequency, respectively. The insets provide the enlarged details.</p>
Full article ">Figure 5
<p>The results of the signal output after filtering and normalization measured during the day.</p>
Full article ">Figure 6
<p>The examples of the results for the Yangxi earthquake swarm recorded at UTC 2021-08-11 17:16:36. The solid vertical lines mark the time of three microseisms ML = 1.1 at UTC 2021-08-11 17:18:11 (white line), 2.6 at UTC 2021-08-11 17:18:17 (red line), and 1.2 at UTC 2021-08-11 17:19:33 (orange line).</p>
Full article ">Figure 7
<p>The recorded results for the bridge section. (<b>a</b>) A waterfall plot of the data in the bridge section; the rectangular box marks the vibration signal observed after spectral filtering. (<b>b</b>) The examples of the recorded vibration signals.</p>
Full article ">Figure 8
<p>The results of the data for no microseisms and with microseism swarms in the last ocean section. Original plots of (<b>a</b>) no microseisms and (<b>b</b>) with microseisms. The signal output of (<b>c</b>) no microseisms and (<b>d</b>) with microseisms after filtering and normalization. The solid vertical lines mark the time of three microseisms.</p>
Full article ">
14 pages, 13165 KiB  
Article
Detection and Monitoring of Mining-Induced Seismicity Based on Machine Learning and Template Matching: A Case Study from Dongchuan Copper Mine, China
by Tao Wu, Zhikun Liu and Shaopeng Yan
Sensors 2024, 24(22), 7312; https://doi.org/10.3390/s24227312 - 15 Nov 2024
Viewed by 308
Abstract
The detection and monitoring of mining-induced seismicity are essential for understanding the mechanisms behind earthquakes and mitigating seismic hazards. However, traditional underground seismic monitoring networks for mining-induced seismicity are challenging to install and operate, which has limited their widespread application. In recent years, [...] Read more.
The detection and monitoring of mining-induced seismicity are essential for understanding the mechanisms behind earthquakes and mitigating seismic hazards. However, traditional underground seismic monitoring networks for mining-induced seismicity are challenging to install and operate, which has limited their widespread application. In recent years, an alternative approach has emerged: utilizing dense seismic arrays at the surface to monitor mining-induced seismicity. This paper proposes a rapid and efficient data processing scheme for the detection and monitoring of mining-induced seismicity based on the surface dense array. The proposed workflow includes machine learning-based phase picking and P-wave first-motion-polarity picking, followed by rapid phase association, precise earthquake location, and template matching for detecting small earthquakes to enhance the completeness of the earthquake catalog. Additionally, it also provides focal mechanism solutions for larger mining-induced events. We applied this workflow to the continuous waveform data from 90 seismic stations over a period of 27 days around the Dongchuan Copper Mine, Yunnan Province, China. Our results yielded 1536 high-quality earthquake locations and two focal mechanism solutions for larger events. By analyzing the spatiotemporal distribution of these events, we are able to investigate the mechanisms of the induced seismic clusters near the Shijiangjun and Lanniping deposits. Our findings highlight the excellent monitoring capability and application potential of the workflow based on machine learning and template matching compared with conventional techniques. Full article
(This article belongs to the Special Issue Sensors and Sensing Technologies for Seismic Detection and Monitoring)
Show Figures

Figure 1

Figure 1
<p>Workflow diagram showing the detection and monitoring of mining-induced earthquakes.</p>
Full article ">Figure 2
<p>(<b>a</b>) Location of Dongchuan Copper Mine in China. (<b>b</b>) Distribution of deposits in Dongchuan Copper Mine and seismic stations used in this study. Blue dots indicate the epicenter of the regional network catalog from 2009 to 2021. Abbreviations: SKS, Sikeshu; YKS, Yikeshu; YM, Yinmin; LX, Luoxue; SJJ, Shijiangjun; LNP, Lanniping; BXL, Baixila.</p>
Full article ">Figure 3
<p>An example of machine learning-based phase picking. (<b>a</b>) A segment of 30 s waveforms starting from 02:36:30. (<b>b</b>) Probabilities of P-wave phase (blue) and S-wave phase (red). The picking probabilities threshold is set to 0.3 in this study. The event near 20:36:33 can be detected due to its high probability; however, the event within the red rectangle at 02:36:50 cannot be detected.</p>
Full article ">Figure 4
<p>(<b>a</b>) The 1D velocity model used for phase association. (<b>b</b>) Travel time–hypocentral distance curves of 856 associated earthquakes.</p>
Full article ">Figure 5
<p>Earthquake catalog comparison between (<b>a</b>) phase association, (<b>b</b>) absolute location, (<b>c</b>) relative location, and (<b>d</b>) template matching. Yellow dots indicate the Cu deposits. Red dots indicate the epicenter of seismic events. Open triangles indicate the short-period stations. Green triangles indicate the broadband stations. Black solid triangle indicates the reginal station.</p>
Full article ">Figure 6
<p>Magnitude–time plot of seismicity during the entire study period.</p>
Full article ">Figure 7
<p>Comparison of magnitude completeness between regional network catalog and dense array catalog obtained in this study.</p>
Full article ">Figure 8
<p>High-precision earthquake catalog (same as <a href="#sensors-24-07312-f005" class="html-fig">Figure 5</a>d) around the Dongchuan Copper Mines using a dense seismic array, machine learning, and template matching. (<b>a</b>) Map view. (<b>b</b>) West–east cross-section. (<b>c</b>) North–south cross-section. (<b>d</b>) Enlarged view of SJJ cluster. (<b>e</b>) Enlarged view of LNP cluster. Beach balls indicate the focal mechanism. Yellow dots indicate the Cu deposits. Red dots indicate the epicenter of seismic events. Open triangles indicate the seismic stations.</p>
Full article ">Figure 9
<p>(<b>a</b>) The 3D view of the SJJ (red) and LNP (orange) clusters. (<b>b</b>) The projections of the SJJ and LNP clusters on each plane in 3D space.</p>
Full article ">Figure 10
<p>Cumulative number of seismicity and seismicity rate per day for (<b>a</b>) SJJ cluster and (<b>b</b>) LNP cluster, respectively.</p>
Full article ">
19 pages, 10355 KiB  
Article
A Case Study Comparing Methods for Coal Thickness Identification in Complex Geological Conditions
by Tao Ding, Yanhui Wu, Lei Wang, Zhen Nie and Lei Zhang
Appl. Sci. 2024, 14(22), 10381; https://doi.org/10.3390/app142210381 - 12 Nov 2024
Viewed by 377
Abstract
This study compares the effectiveness of different methods for coal thickness identification, aiming to identify the most accurate approach and provide a reference for intelligent coalmine development. Focused on the No. 2 coal seam in a mining area in Shanxi, China, the analysis [...] Read more.
This study compares the effectiveness of different methods for coal thickness identification, aiming to identify the most accurate approach and provide a reference for intelligent coalmine development. Focused on the No. 2 coal seam in a mining area in Shanxi, China, the analysis employs well log-constrained impedance inversion and seismic multi-attribute techniques. The results show that the back propagation (BP) neural network model, as part of the seismic multi-attribute approach, delivers prediction accuracy comparable to the well log-constrained inversion method. Specifically, after applying proper static corrections, a four-layer BP neural network was constructed using four optimized sensitive attributes as the input layer, achieving an error range of 0.11% to 1.36%, compared to 0.03% to 6.59% for the logging-based method. The BP neural network demonstrated strong applicability in complex geological environments. Empirical analysis further validated the BP neural network’s geological reliability and practicality in systematic coal thickness determination. Full article
Show Figures

Figure 1

Figure 1
<p>Typical landforms of the exploration area.</p>
Full article ">Figure 2
<p>Typical single-shot record and seismic time section in this area.</p>
Full article ">Figure 3
<p>Coverage in the exploration area. The color scale represents the number of stacking iterations.</p>
Full article ">Figure 4
<p>Comparison of static correction before and after: (<b>a</b>) before static correction; and (<b>b</b>) after static correction.</p>
Full article ">Figure 5
<p>Workflow diagram for well log-based impedance inversion.</p>
Full article ">Figure 6
<p>Neural network model.</p>
Full article ">Figure 7
<p>Extraction of wavelet and spectrum for the first time.</p>
Full article ">Figure 8
<p>Final wavelet and spectrum.</p>
Full article ">Figure 9
<p>Comparison between synthetic records and seismic profiles: (<b>a</b>) density log; (<b>b</b>) velocity log; (<b>c</b>) seismogram synthesis; and (<b>d</b>) seismic trace.</p>
Full article ">Figure 10
<p>Section of impedance.</p>
Full article ">Figure 11
<p>Cross-plot analysis of effectiveness for attributes of the No. 2 coal seam.</p>
Full article ">Figure 12
<p>Coal thickness prediction results.</p>
Full article ">Figure 13
<p>Comparison of error results.</p>
Full article ">Figure 14
<p>Error replacement diagram for one attribute.</p>
Full article ">Figure 15
<p>Error diagram for replacement of two attribute combinations.</p>
Full article ">Figure 16
<p>Elevation distribution map of the entire area.</p>
Full article ">Figure 17
<p>Statics correction distribution map of the entire area.</p>
Full article ">
18 pages, 18719 KiB  
Article
Seismic Vibration Control and Multi-Objective Optimization of Transmission Tower with Tuned Mass Damper Under Near-Fault Pulse-like Ground Motions
by Ying Lin and Tao Liu
Buildings 2024, 14(11), 3572; https://doi.org/10.3390/buildings14113572 - 10 Nov 2024
Viewed by 457
Abstract
Although the wind load is usually adopted as the governing lateral load in the design of transmission towers, many tall transmission towers may be damaged or even collapse in high seismic intensity areas, especially under near-fault pulse-like ground motions. To study the seismic [...] Read more.
Although the wind load is usually adopted as the governing lateral load in the design of transmission towers, many tall transmission towers may be damaged or even collapse in high seismic intensity areas, especially under near-fault pulse-like ground motions. To study the seismic vibration control effect of a tuned mass damper (TMD) attached to transmission tower, parametric analyses are conducted in SAP2000 through CSI OAPI programming, including TMD parameters such as the mass ratio μ from 0.5% to 10%, the frequency ratio f from 0.7 to 1.2, and the damping ratio ξ from 0.01 to 0.2. Based on the obtained analysis results, artificial neural network (ANN) is trained to predict the vibration reduction ratios of peak responses and the corresponding vibration reduction cost. Finally, the NSGA-III algorithm is adopted to perform the multi-objective optimization of a transmission tower equipped with TMD. Results show that the vibration reduction ratios first increase and then decrease with the increase of frequency ratio, but first increase and then remain stable with the increase of mass ratio and damping ratio. In addition, ANN fitting can accurately predict the nonlinear relationship between TMD parameters and objective functions. Through multi-objective optimization with the NSGA-III algorithm, TMD can simultaneously and significantly reduce different peak responses of transmission towers under near-fault pulse-like ground motions in a cost-effective manner. Full article
(This article belongs to the Section Building Structures)
Show Figures

Figure 1

Figure 1
<p>Research flowchart of this study [<a href="#B44-buildings-14-03572" class="html-bibr">44</a>].</p>
Full article ">Figure 2
<p>Finite element model and the first three vibration modes of transmission tower.</p>
Full article ">Figure 3
<p>Schematic diagram of transmission tower with TMD.</p>
Full article ">Figure 4
<p>Acceleration time histories of different near-fault ground motions.</p>
Full article ">Figure 5
<p>Acceleration response spectra of different near-fault ground motions [<a href="#B44-buildings-14-03572" class="html-bibr">44</a>].</p>
Full article ">Figure 6
<p>Effect of mass ratio on vibration reduction ratios.</p>
Full article ">Figure 7
<p>Effect of frequency ratio on vibration reduction ratios.</p>
Full article ">Figure 8
<p>Effect of damping ratio on vibration reduction ratios.</p>
Full article ">Figure 9
<p>Vibration reduction ratio and cost of transmission tower with TMD: (<b>a</b>) <span class="html-italic">VRR<sub>d</sub></span>; (<b>b</b>) <span class="html-italic">VRR<sub>v</sub></span>; (<b>c</b>) <span class="html-italic">VRR<sub>a</sub></span>; (<b>d</b>) <span class="html-italic">VR<sub>c</sub></span>.</p>
Full article ">Figure 9 Cont.
<p>Vibration reduction ratio and cost of transmission tower with TMD: (<b>a</b>) <span class="html-italic">VRR<sub>d</sub></span>; (<b>b</b>) <span class="html-italic">VRR<sub>v</sub></span>; (<b>c</b>) <span class="html-italic">VRR<sub>a</sub></span>; (<b>d</b>) <span class="html-italic">VR<sub>c</sub></span>.</p>
Full article ">Figure 10
<p>ANN fitting effects of vibration reduction ratios and cost: (<b>a</b>) <span class="html-italic">VRR<sub>d</sub></span>; (<b>b</b>) <span class="html-italic">VRR<sub>v</sub></span>; (<b>c</b>) <span class="html-italic">VRR<sub>a</sub></span>; (<b>d</b>) <span class="html-italic">VR<sub>c</sub></span>.</p>
Full article ">Figure 10 Cont.
<p>ANN fitting effects of vibration reduction ratios and cost: (<b>a</b>) <span class="html-italic">VRR<sub>d</sub></span>; (<b>b</b>) <span class="html-italic">VRR<sub>v</sub></span>; (<b>c</b>) <span class="html-italic">VRR<sub>a</sub></span>; (<b>d</b>) <span class="html-italic">VR<sub>c</sub></span>.</p>
Full article ">Figure 11
<p>Optimal variables obtained by NSGA-III algorithm.</p>
Full article ">Figure 12
<p>Pareto front of <span class="html-italic">VRR<sub>d</sub></span>, <span class="html-italic">VRR<sub>v</sub></span>, <span class="html-italic">VRR<sub>a</sub></span>, and <span class="html-italic">VR<sub>c</sub></span>.</p>
Full article ">Figure 13
<p>Vibration reduction ratios of peak responses under different near-fault ground motions.</p>
Full article ">Figure 14
<p>Energy response of transmission tower: (<b>a</b>) without TMD and (<b>b</b>) with TMD.</p>
Full article ">Figure 15
<p>Time history comparisons of peak response at the top of transmission tower under RSN 6962 ground motion: (<b>a</b>) peak displacement, (<b>b</b>) peak velocity, (<b>c</b>) peak acceleration.</p>
Full article ">
14 pages, 15387 KiB  
Article
Optimization and Numerical Verification of Microseismic Monitoring Sensor Network in Underground Mining: A Case Study
by Chenglu Hou, Xibing Li, Yang Chen, Wei Li, Kaiqu Liu, Longjun Dong and Daoyuan Sun
Mathematics 2024, 12(22), 3500; https://doi.org/10.3390/math12223500 - 9 Nov 2024
Viewed by 434
Abstract
A scientific and reasonable microseismic monitoring sensor network is crucial for the prevention and control of rockmass instability disasters. In this study, three feasible sensor network layout schemes for the microseismic monitoring of Sanshandao Gold Mine were proposed, comprehensively considering factors such as [...] Read more.
A scientific and reasonable microseismic monitoring sensor network is crucial for the prevention and control of rockmass instability disasters. In this study, three feasible sensor network layout schemes for the microseismic monitoring of Sanshandao Gold Mine were proposed, comprehensively considering factors such as orebody orientation, tunnel and stope distributions, blasting excavation areas, construction difficulty, and maintenance costs. To evaluate and validate the monitoring effectiveness of the sensor networks, three layers of seismic sources were randomly generated within the network. Four levels of random errors were added to the calculated arrival time data, and the classical Geiger localization algorithm was used for locating validation. The distribution of localization errors within the monitoring area was analyzed. The results indicate that when the arrival time data are accurate or the error is between 0% and 2%, scheme 3 is considered the most suitable layout; when the error of the arrival time data is between 2% and 10%, scheme 2 is considered the optimal layout. These research results can provide important theoretical and technical guidance for the reasonable design of microseismic monitoring systems in similar mines or projects. Full article
(This article belongs to the Special Issue Numerical Model and Artificial Intelligence in Mining Engineering)
Show Figures

Figure 1

Figure 1
<p>Scheme 1 of sensor monitoring network. (<b>a</b>) Sensor level 1, <span class="html-italic">z</span> = 0 m. (<b>b</b>) Sensor level 2, <span class="html-italic">z</span> = 15 m. (<b>c</b>) Sensor level 3, <span class="html-italic">z</span> = 30 m. (<b>d</b>) Sensor level 4, <span class="html-italic">z</span> = 45 m.</p>
Full article ">Figure 2
<p>Scheme 2 of sensor monitoring network. (<b>a</b>) Sensor level 1, <span class="html-italic">z</span> = 0 m. (<b>b</b>) Sensor level 2, <span class="html-italic">z</span> = 15 m. (<b>c</b>) Sensor level 3, <span class="html-italic">z</span> = 30 m. (<b>d</b>) Sensor level 4, <span class="html-italic">z</span> = 45 m.</p>
Full article ">Figure 3
<p>Scheme 3 of sensor monitoring network. (<b>a</b>) Sensor level 1, <span class="html-italic">z</span> = 0 m. (<b>b</b>) Sensor level 2, <span class="html-italic">z</span> = 15 m. (<b>c</b>) Sensor level 3, <span class="html-italic">z</span> = 30 m. (<b>d</b>) Sensor level 4, <span class="html-italic">z</span> = 45 m.</p>
Full article ">Figure 4
<p>The distribution of locating errors under different schemes when the arrival times are accurate.</p>
Full article ">Figure 5
<p>The distribution of locating errors under different schemes when the arrival time errors range from 0% to 2%.</p>
Full article ">Figure 6
<p>The distribution of locating errors under different schemes when the arrival time errors range from 2% to 5%.</p>
Full article ">Figure 7
<p>The distribution of locating errors under different schemes when the arrival time errors range from 5% to 8%.</p>
Full article ">Figure 8
<p>The distribution of locating errors under different schemes when the arrival time errors range from 8% to 10%.</p>
Full article ">
15 pages, 12447 KiB  
Article
Iterative Separation of Blended Seismic Data in Shot Domain Using Deep Learning
by Liyun Ma, Liguo Han and Pan Zhang
Remote Sens. 2024, 16(22), 4167; https://doi.org/10.3390/rs16224167 - 8 Nov 2024
Viewed by 331
Abstract
Accurate deblending techniques are essential for the successful application of blended seismic acquisition. Deep-learning-based deblending methods typically begin by performing a pseudo-deblending operation on blended data, followed by further processing in either the common-shot domain or a non-common-shot domain. In this study, we [...] Read more.
Accurate deblending techniques are essential for the successful application of blended seismic acquisition. Deep-learning-based deblending methods typically begin by performing a pseudo-deblending operation on blended data, followed by further processing in either the common-shot domain or a non-common-shot domain. In this study, we propose an iterative deblending framework based on deep learning, which directly addresses the blended data in the shot domain, eliminating the need for pseudo-deblending and domain transformation. This framework is built around a unique architecture, termed WNETR, which derives its name from its W-shaped network structure that combines U-Net and Transformer. During testing, the trained WNETR is incorporated into the iterative framework to extract useful signals iteratively. Tests on synthetic data validate the effectiveness of the proposed deblending iterative framework. Full article
Show Figures

Figure 1

Figure 1
<p>Blended acquisition.</p>
Full article ">Figure 2
<p>Velocity model. (<b>a</b>) Velocity models for generating the training set. (<b>b</b>) Velocity model for generating the test set.</p>
Full article ">Figure 3
<p>The generation of blended data. The four colored dots indicate seismic sources classified into four groups.</p>
Full article ">Figure 4
<p>The blended data and their labels. (<b>a</b>) Blended seismic data from four sources. (<b>b</b>) Deblending results as training label 1. (<b>c</b>) Blending noise as training label 2.</p>
Full article ">Figure 5
<p>The blended data and their labels. (<b>a</b>) Blending noise in <a href="#remotesensing-16-04167-f004" class="html-fig">Figure 4</a>c serves as the blended seismic data for the three sources. (<b>b</b>) Deblending results as training label 1. (<b>c</b>) Blending noise as training label 2.</p>
Full article ">Figure 6
<p>The blended data and their labels. (<b>a</b>) Blending noise in <a href="#remotesensing-16-04167-f005" class="html-fig">Figure 5</a>c serves as the blended seismic data for the two sources. (<b>b</b>) Deblending results as training label 1. (<b>c</b>) Blending noise as training label 2.</p>
Full article ">Figure 7
<p>The architecture of WNETR.</p>
Full article ">Figure 8
<p>The iterative framework for deblending. <math display="inline"><semantics> <mrow> <mi>N</mi> </mrow> </semantics></math> represents the maximum iteration number and <math display="inline"><semantics> <mrow> <mi>n</mi> </mrow> </semantics></math> represents the current number of iterations.</p>
Full article ">Figure 9
<p>(<b>a</b>) A test sample with a mixing degree of 2. (<b>b</b>) The predicted deblending results of the first source. (<b>c</b>) The predicted deblending results of the second source.</p>
Full article ">Figure 10
<p>(<b>a</b>) Residuals of the deblending results of the first source and the original unblended data. (<b>b</b>) Residuals of the deblending results of the second source and the original unblended data.</p>
Full article ">Figure 11
<p>The blended data containing four sources.</p>
Full article ">Figure 12
<p>The result of the first iteration. (<b>a</b>) Deblending results. (<b>b</b>) Blending noise.</p>
Full article ">Figure 13
<p>The result of the second iteration. (<b>a</b>) Deblending results. (<b>b</b>) Blending noise.</p>
Full article ">Figure 14
<p>The result of the third iteration. (<b>a</b>) Deblending results. (<b>b</b>) Blending noise.</p>
Full article ">Figure 15
<p>Comparison of deblending results. (<b>a</b>) The useful signal of the first source obtained by our method. (<b>b</b>) The useful signal of the first source obtained by a method without self-supervised learning. (<b>c</b>) The useful signal of the second source obtained by our method. (<b>d</b>) The useful signal of the second source obtained by a method without self-supervised learning.</p>
Full article ">
15 pages, 10808 KiB  
Article
A Strong Noise Reduction Network for Seismic Records
by Tong Shen, Xuan Jiang, Wenzheng Rong, Lei Xu, Xianguo Tuo and Guili Peng
Appl. Sci. 2024, 14(22), 10262; https://doi.org/10.3390/app142210262 - 7 Nov 2024
Viewed by 369
Abstract
Noise reduction is a critical step in seismic data processing. A novel strong noise reduction network is proposed in this study. The network enhances the U-Net architecture with an improved inception module and coordinate attention (CA) mechanism, suppressing noise and enhancing signal clarity. [...] Read more.
Noise reduction is a critical step in seismic data processing. A novel strong noise reduction network is proposed in this study. The network enhances the U-Net architecture with an improved inception module and coordinate attention (CA) mechanism, suppressing noise and enhancing signal clarity. These enhancements improve the network’s capability to distinguish between signal and noise in the time–frequency domain. We trained and tested our model on the STEAD dataset, which eliminated noise across various frequency bands, improved the signal-to-noise ratio (SNR) of seismic records, and reduced the waveform distortion significantly. Comparative analyses against U-Net, DeepDenoiser, and DnRDB models, using signals with SNRs ranging from −14 dB to 0 dB, demonstrated our model’s superior performance. At the same time, we demonstrated that the Inception Conv Block has a significant impact on the denoising ability of the network. Furthermore, validation using the “Di Ting” dataset and real noisy signals confirmed the model’s generalizability. These results show that the proposed model significantly outperforms the comparative methods in terms of the SNR, correlation coefficient (r), and root mean square error (RMSE), delivering higher-quality seismograms. The enhanced phase-picking accuracy underscores the potential of our approach to advance in geophysics applications. Full article
(This article belongs to the Topic Remote Sensing and Geological Disasters)
Show Figures

Figure 1

Figure 1
<p>The detailed structure diagram of the proposed network. (<b>a</b>) The structure of our network. The inputs are the real and imaginary parts of the noise data after the STFT, and the outputs are the masks corresponding to the signal and noise. The different colored squares represent different modules in the network and the numbers represent the size of the feature map after the output of each layer of modules. (<b>b</b>) The inception module used in our network. (<b>c</b>) The CA module used in our network.</p>
Full article ">Figure 2
<p>The left panels in (<b>a</b>,<b>b</b>,<b>d</b>–<b>f</b>) display the clean signal, noisy signal, denoised signal, and recovered noise, respectively, in the time domain. The right panels present the corresponding spectrograms. (<b>c</b>) The mask of the signals and noises depicted in (<b>a</b>,<b>b</b>).</p>
Full article ">Figure 3
<p>Training and validation loss during the training process.</p>
Full article ">Figure 4
<p>Denoising performance of low-frequency noise: (<b>I</b>,<b>III</b>,<b>V</b>) are the original signal, the noise signal, and the denoised signal in the time domain; and (<b>II</b>,<b>IV</b>,<b>VI</b>) are the time–frequency domain data of (<b>I</b>,<b>III</b>,<b>V</b>).</p>
Full article ">Figure 5
<p>Denoising performance of high-frequency noise: (<b>I</b>,<b>III</b>,<b>V</b>) are the original signal, the noise signal, and the denoised signal in the time domain; and (<b>II</b>,<b>IV</b>,<b>VI</b>) are the time–frequency domain data of (<b>I</b>,<b>III</b>,<b>V</b>).</p>
Full article ">Figure 6
<p>Denoising performance of mixed-frequency noise: (<b>I</b>,<b>III</b>,<b>V</b>) are the original signal, the noise signal, and the denoised signal in the time domain; and (<b>II</b>,<b>IV</b>,<b>VI</b>) are the time–frequency domain data of (<b>I</b>,<b>III</b>,<b>V</b>).</p>
Full article ">Figure 7
<p>Denoising performance of spike noise: (<b>I</b>,<b>III</b>,<b>V</b>) are the original signal, the noise signal, and the denoised signal in the time domain; and (<b>II</b>,<b>IV</b>,<b>VI</b>) are the time–frequency domain data of (<b>I</b>,<b>III</b>,<b>V</b>).</p>
Full article ">Figure 8
<p>Comparison of the denoising effects of different networks: (<b>a</b>) improvement of the SNR; (<b>b</b>) correlation coefficient; (<b>c</b>) root mean square error; and (<b>d</b>) arrival time error.</p>
Full article ">Figure 9
<p>Waveform comparison of the noise reduction results of various methods for an SNR of 10 dB: (<b>a</b>) noisy signal; (<b>b</b>) original signal; (<b>c</b>) U-Net denoising result; (<b>d</b>) DnRDB denoising result; (<b>e</b>) DeepDenoiser denoising result; and (<b>f</b>) our model’s denoising result. The red line represents the true first arrival time of the P-wave, whereas the blue line represents the result of picking up the first arrival of the denoised signals from each network using the STA/LTA algorithm [<a href="#B53-applsci-14-10262" class="html-bibr">53</a>].</p>
Full article ">Figure 10
<p>P-wave picking results of various models after denoising real noisy signals.</p>
Full article ">
22 pages, 23011 KiB  
Article
Removing Instrumental Noise in Distributed Acoustic Sensing Data: A Comparison Between Two Deep Learning Approaches
by Xihao Gu, Olivia Collet, Konstantin Tertyshnikov and Roman Pevzner
Remote Sens. 2024, 16(22), 4150; https://doi.org/10.3390/rs16224150 - 7 Nov 2024
Viewed by 503
Abstract
Over the last decade, distributed acoustic sensing (DAS) has received growing attention in the field of seismic acquisition and monitoring due to its potential high spatial sampling rate, low maintenance cost and high resistance to temperature and pressure. Despite its undeniable advantages, DAS [...] Read more.
Over the last decade, distributed acoustic sensing (DAS) has received growing attention in the field of seismic acquisition and monitoring due to its potential high spatial sampling rate, low maintenance cost and high resistance to temperature and pressure. Despite its undeniable advantages, DAS faces some challenges, including a low signal-to-noise ratio, which partly results from the instrument-specific noise generated by DAS interrogators. We present a comparison between two deep learning approaches to address DAS hardware noise and enhance the quality of DAS data. These approaches have the advantage of including real instrumental noise in the neural network training dataset. For the supervised learning (SL) approach, real DAS instrumental noise measured on an acoustically isolated coil is added to synthetic data to generate training pairs of clean/noisy data. For the second method, the Noise2Noise (N2N) approach, the training is performed on noisy/noisy data pairs recorded simultaneously on the downgoing and upgoing parts of a downhole fiber-optic cable. Both approaches allow for the removal of unwanted noise that lies within the same frequency band of the useful signal, a result that cannot be achieved by conventional denoising techniques employing frequency filtering. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Noise recorded on an acoustically isolated coil of fiber; (<b>b</b>) corresponding strain rate spectral density.</p>
Full article ">Figure 2
<p>Flowchart for constructing semi-synthetic DAS data using synthetic data and recorded instrumental noise.</p>
Full article ">Figure 3
<p>Schematics of the borehole optic cable used for the acquisition of the Noise2Noise training dataset. The cable is looped at the bottom of the well; the downgoing and upgoing parts of the cable are depicted in red and orange, respectively.</p>
Full article ">Figure 4
<p>Illustration of the data preparation workflow and neural network training for the N2N approach.</p>
Full article ">Figure 5
<p>Denoising performance of the SL neural network trained on semi-synthetic data. (<b>a</b>) Clean data, (<b>b</b>) input semi-synthetic data, (<b>c</b>) denoised result and (<b>d</b>) removed noise. To quantitatively assess the denoising performance, SNR sections are computed for (<b>e</b>) the input data, (<b>f</b>) the denoised result and (<b>g</b>) the removed noise. SNR values at the top of the plots indicate the average SNR value for the corresponding section.</p>
Full article ">Figure 6
<p>Comparison of strain rate spectral densities calculated for the input data (blue line), clean data (red line) and the data denoised using the SL network trained on semi-synthetic data (orange line). The removed noise (green dotted line) and the noise recorded on an acoustically isolated coil (black dotted line) are also displayed. The corresponding gathers for the data are shown in <a href="#remotesensing-16-04150-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 7
<p>Denoising performance of the N2N-trained neural network on the Otway semi-synthetic dataset: (<b>a</b>) clean data, (<b>b</b>) input gather, (<b>c</b>) denoised result and (<b>d</b>) removed noise. Corresponding SNR attribute sections are shown in panels (<b>e</b>–<b>g</b>). The average SNR values are displayed on top of each section.</p>
Full article ">Figure 8
<p>Strain rate spectral densities of the noisy input (blue line), N2N-denoised (orange line), and clean input (red lines) data for the Otway semi-synthetic record example. The spectral densities of the removed noise (dotted green line) and the recorded noise (dotted black line) that were originally added to the clean synthetic data are also compared.</p>
Full article ">Figure 9
<p>SL neural network application to the DAS-VSP data acquired with a low-power active source: (<b>a</b>) raw DAS shot record, (<b>b</b>) denoised result and (<b>c</b>) removed noise. The corresponding SNR sections are shown in figures (<b>d</b>–<b>f</b>). SNR values at the top of the plots indicate the average SNR value for the corresponding section.</p>
Full article ">Figure 10
<p>Comparison of strain rate spectral densities calculated for the input data (blue line), the data denoised using the SL supervised network (orange line) and the removed noise (dotted green line). The corresponding gathers for the input data, denoised data and removed noise data are shown in <a href="#remotesensing-16-04150-f009" class="html-fig">Figure 9</a>.</p>
Full article ">Figure 11
<p>N2N neural network application to the DAS VSP data example used in training: (<b>a</b>) raw DAS shot record, (<b>b</b>) denoised result and (<b>c</b>) removed noise. The corresponding SNR sections are shown in figures (<b>d</b>–<b>f</b>). SNR values at the top of the plots indicate the average SNR value for the corresponding section.</p>
Full article ">Figure 12
<p>Comparison of strain rate spectral densities calculated for the input data (blue line), the data denoised using the N2N approach (orange line) and the removed noise (green dotted line). The corresponding gathers for the input data, denoised data and removed noise data are shown in <a href="#remotesensing-16-04150-f011" class="html-fig">Figure 11</a>.</p>
Full article ">Figure 13
<p>Denoising performance of the SL neural network on an example of the microseismic DAS data: (<b>a</b>) input data, (<b>b</b>) denoised result and (<b>c</b>) removed noise. To quantitatively assess denoising ability, SNR sections are computed for (<b>d</b>) input data, (<b>e</b>) denoising results and (<b>f</b>) removed noise. The average SNR values are displayed at the top of each section.</p>
Full article ">Figure 14
<p>Denoising performance of the N2N-trained neural network on a microseismic event recorded in Otway CRC-7 well: (<b>a</b>) input, (<b>b</b>) denoised result and (<b>c</b>) removed noise. The black rectangles delineate the magnified sections shown in <a href="#remotesensing-16-04150-f011" class="html-fig">Figure 11</a>. Corresponding SNR attribute sections are shown in panels (<b>d</b>–<b>f</b>). The average SNR values are displayed at the top of each section.</p>
Full article ">Figure 15
<p>Magnified sections of denoising results of the N2N-trained neural network on a microseismic event recorded in Otway CRC-7 well: (<b>a</b>) input, (<b>b</b>) denoised result and (<b>c</b>) removed noise.</p>
Full article ">Figure 16
<p>Comparison of the denoising results obtained for the microseismic event recorded in Otway CRC-7 well using (<b>a</b>,<b>e</b>) the SL approach, (<b>b</b>,<b>f</b>) N2N trained neural network, (<b>c</b>,<b>g</b>) bandpass filtering and (<b>d</b>,<b>h</b>) FCDNet. The upper row shows the denoised sections, while the lower row shows the removed noise sections. The black rectangles in panels (<b>a</b>–<b>d</b>) delineate the magnified sections shown in <a href="#remotesensing-16-04150-f017" class="html-fig">Figure 17</a>.</p>
Full article ">Figure 17
<p>Magnified sections of the denoised results obtained using (<b>a</b>) the SL-trained neural network, (<b>b</b>) the N2N-trained neural network, (<b>c</b>) bandpass filtering and (<b>d</b>) FCDNet for the microseismic event recorded in Otway CRC-7 well.</p>
Full article ">Figure 18
<p>Comparison between (<b>a</b>,<b>d</b>) N2N trained with high-frequency semi-synthetic data, (<b>b</b>,<b>e</b>) the SL approach trained with high-frequency semi-synthetic data (160 Hz Ricker wavelet) and (<b>c</b>,<b>f</b>) the SL approach trained with lower frequency band semi-synthetic data (100 Hz Ricker wavelet). The upper row shows the denoised sections, while the lower row shows the removed noise sections.</p>
Full article ">Figure 19
<p>Comparison between the current N2N-trained network (<b>a</b>,<b>c</b>) and the original DAS-N2N network trained by Lapins et al. [<a href="#B23-remotesensing-16-04150" class="html-bibr">23</a>] (<b>b</b>,<b>d</b>). Figure (<b>a</b>,<b>b</b>) show the denoised sections, while figure (<b>c</b>,<b>d</b>) show the removed noise sections.</p>
Full article ">
12 pages, 2745 KiB  
Article
Single-Shot Time-Lapse Target-Oriented Velocity Inversion Using Machine Learning
by Katerine Rincon, Ramon C. F. Araújo, Moisés M. Galvão, Samuel Xavier-de-Souza, João M. de Araújo, Tiago Barros and Gilberto Corso
Appl. Sci. 2024, 14(21), 10047; https://doi.org/10.3390/app142110047 - 4 Nov 2024
Viewed by 503
Abstract
In this study, we used machine learning (ML) to estimate time-lapse velocity variations in a reservoir region using seismic data. To accomplish this task, we needed an adequate training set that could map seismic data to velocity perturbation. We generated a synthetic seismic [...] Read more.
In this study, we used machine learning (ML) to estimate time-lapse velocity variations in a reservoir region using seismic data. To accomplish this task, we needed an adequate training set that could map seismic data to velocity perturbation. We generated a synthetic seismic database by simulating reservoirs of varying velocities using a 2D velocity model typical of the Brazilian pre-salt ocean bottom node (OBN) acquisition, located in the Santos basin, Brazil. The largest velocity change in the injector well was around 3% of the empirical velocity model, which mimicked a realistic scenario. The acquisition geometry was formed by the geometry of 1 shot and 49 receivers. For each synthetic reservoir, the corresponding seismic data were obtained by estimating a one-shot forward-wave propagation using acoustic approximation. We studied the reservoir illumination to optimize the input data of the ML inversion. We split the set of synthetic reservoirs into two subsets: training (80%) and testing (20%) sets. We point out that the ML inversion was restricted to the reservoir zone, which means that it was inversion-oriented to a target. We obtained a good similarity between true and ML-inverted reservoir anomalies. The similarity diminished for a situation with non-repeatability noise. Full article
(This article belongs to the Section Earth Sciences)
Show Figures

Figure 1

Figure 1
<p>P-wave velocity model and acquisition geometry used to model the synthetic database. The rectangle indicates the reservoir region.</p>
Full article ">Figure 2
<p>Illumination map produced by a single seismic source. This color map should be interpreted as areas where the energy of the seismic wave pass by and is captured by the receptors: (<b>a</b>) traces of sequence numbers from 11 to 20 (inclusive), (<b>b</b>) the 5 smallest offset traces, (<b>c</b>) traces from 21 to 30, and (<b>d</b>) traces from 31 to 40. The yellow rectangle indicates the target reservoir region as in <a href="#applsci-14-10047-f001" class="html-fig">Figure 1</a>. The illumination map reveals that this simple acquisition geometry provides reasonable information to invert the target region.</p>
Full article ">Figure 3
<p>Block diagram of the proposed ML inversion methodology. The input set consists of the seismic time-lapse (TL) difference and the output consists of the inverted reservoir anomaly. The red dashed rectangles indicate the subsets of time-lapse differences considered to calculate input and target data: we extract from seismic differences the time window concentrating most of the reflection energy coming from the reservoir, and from the velocity differences we consider the spatial region of the reservoir.</p>
Full article ">Figure 4
<p>Proposed neural network architecture. The input of the ML is the seismic time-lapse difference within the reservoir time window and the output is the velocity anomaly of the target region.</p>
Full article ">Figure 5
<p>Loss curves computed on the training and validation subsets during the training of the neural network employed for velocity inversion. The loss function is calculated with velocities in the scale of the 4D anomalies (0–100 m/s range).</p>
Full article ">Figure 6
<p>Inversion results of individual samples for the perfect repeatability scenario. The yellow pattern in the panels represent the velocity anomaly, the horizontal and vertical dimensions reproduce the reservoir region indicated in the rectangle of <a href="#applsci-14-10047-f001" class="html-fig">Figure 1</a>. Comparison of true (first row) and inverted (second row) reservoir anomaly for five samples at key points of the SSIM distribution on the test subset: minimum (worst case), 25th percentile, median, 75th percentile, and maximum (best case). The last two rows compare, respectively, the central vertical and central horizontal velocity profiles of the velocity anomalies.</p>
Full article ">Figure 7
<p>Inversion results of individual samples for a scenario with non-repeatability, modeled by randomly moving the receivers in the lateral directions, with maximum perturbations equal to ±0.1 m. The yellow pattern in the panels represent the velocity anomaly; the horizontal and vertical dimensions reproduce the reservoir region indicated by the rectangle in <a href="#applsci-14-10047-f001" class="html-fig">Figure 1</a>. Comparison of true (first row) and predicted (second row) time-lapse velocity anomalies in the target region for specific samples of a test dataset contaminated with geometry non-repeatability noise. The 4D noise was modeled by randomly shifting the receivers in the lateral direction, with maximum perturbations equal to ±0.1 m. The shown samples are located at key percentiles of the SSIM distribution on the referred dataset: minimum (worst case), three quartiles, and maximum (best case). The last two rows compare, respectively, the central vertical and central horizontal velocity profiles of the velocity anomalies.</p>
Full article ">Figure 8
<p>Inversion results of individual samples for a scenario with non-repeatability, modeled by randomly moving the receivers in the lateral directions, with maximum perturbations equal to ±0.5 m. The yellow pattern in the panels represent the velocity anomaly; the horizontal and vertical dimensions reproduce the reservoir region indicated by the rectangle in <a href="#applsci-14-10047-f001" class="html-fig">Figure 1</a>. Comparison of true (first row) and predicted (second row) time-lapse velocity anomalies in the target region for specific samples of a testing dataset contaminated with geometry non-repeatability noise. The illustrated samples are located at key percentiles of the SSIM distribution on the referred dataset: minimum (worst case), three quartiles, and maximum (best case). The last two rows compare, respectively, the central vertical and central horizontal velocity profiles of the velocity anomalies.</p>
Full article ">Figure 9
<p>Spatial distribution of <math display="inline"><semantics> <mrow> <mi>Δ</mi> <mi>v</mi> </mrow> </semantics></math> prediction errors for the test scenarios with (<b>a</b>) perfect repeatability, (<b>b</b>) ±0.1 m geometry 4D noise, and (<b>c</b>) ±0.5 m geometry 4D noise.</p>
Full article ">
16 pages, 3839 KiB  
Article
Hybrid Duplex Medium Access Control Protocol for Tsunami Early Warning Systems in Underwater Networks
by Sung Hyun Park, Ye Je Choi and Tae Ho Im
Electronics 2024, 13(21), 4288; https://doi.org/10.3390/electronics13214288 - 31 Oct 2024
Viewed by 486
Abstract
Tsunamis are devastating natural phenomena that cause extensive damage to both human life and infrastructure. To mitigate such impacts, tsunami early warning systems have been deployed globally. South Korea has also initiated a project to install a tsunami warning system to monitor its [...] Read more.
Tsunamis are devastating natural phenomena that cause extensive damage to both human life and infrastructure. To mitigate such impacts, tsunami early warning systems have been deployed globally. South Korea has also initiated a project to install a tsunami warning system to monitor its surrounding seas. To ensure reliable warning decisions, various types of data must be combined, but efficiently transmitting heterogeneous data poses a challenge due to the unique characteristics of underwater acoustic communication. Therefore, this paper proposes a Hybrid Duplex Medium Access Control (HDMAC) protocol designed for a tsunami warning system, with a specific focus on heterogeneous data transmission. HDMAC efficiently handles both seismic and environmental data by utilizing hybrid duplexing, which combines frequency duplex for seismic data with time duplex for environmental data. The protocol addresses the distinct transmission requirements for each data type by optimizing channel utilization through a group Automatic Repeat request (ARQ) scheme and packet size adjustment. Theoretical analysis predicts that HDMAC can achieve a channel utilization of up to 0.91 in smaller networks and 0.64 in larger networks. HDMAC is validated through simulations, and the simulation results closely match these predictions. The simulation results demonstrate the efficiency of HDMAC in supporting real-time submarine earthquake monitoring systems. Full article
(This article belongs to the Special Issue New Advances in Underwater Communication Systems)
Show Figures

Figure 1

Figure 1
<p>Planned network configuration for HDMAC design. (Explanations related to latency and bandwidth are discussed in <a href="#sec3dot1-electronics-13-04288" class="html-sec">Section 3.1</a>).</p>
Full article ">Figure 2
<p>The DART system. (<b>a</b>) The DART 4th generation system [<a href="#B19-electronics-13-04288" class="html-bibr">19</a>]. It consists of one BPR and one surface buoy. (<b>b</b>) The map indicates where the DART system is deployed [<a href="#B18-electronics-13-04288" class="html-bibr">18</a>]. The yellow diamonds indicates a station with recent data and the red diamond indicates a station with no data in last eight hours. (<b>c</b>) Areas of responsibility of the DART system used for damage prevision [<a href="#B20-electronics-13-04288" class="html-bibr">20</a>]. PTWC stands for Pacific Tsunami Warning Center and NTWC stands for National Tsunami Warning Center.</p>
Full article ">Figure 3
<p>HDMAC channel band allocation plan.</p>
Full article ">Figure 4
<p>HDMAC general process diagram: one frame time (<span class="html-italic">Tf</span>) example. The diagonal lines indicate the propagation delays of each packet.</p>
Full article ">Figure 5
<p>Group ARQ process of HDMAC. The sequence number of the seismic data packet is indicated by i in this figure.</p>
Full article ">Figure 6
<p>A group ARQ process example when the group size (<span class="html-italic">M</span>) is five and when all the transmissions have been successfully transmitted to the buoy. The two SNs are omitted in this figure.</p>
Full article ">Figure 7
<p>HDMAC control packet bit plan: first two bytes for the seismometer and last two bytes for the two SNs.</p>
Full article ">Figure 8
<p>Predicted HDMAC performance when <span class="html-italic">Be</span> is 4000 bits. (<b>a</b>) Theoretical channel utilization of HDMAC as a function of the distance between the buoy and SN1. (<b>b</b>) ARQ group size of HDMAC increases as a function of the distance between the buoy and SN1.</p>
Full article ">Figure 9
<p>The optimized packet size of HDMAC under the project requirement conditions. The optimized packet size decreases as the distance between the buoy and SN1 increases.</p>
Full article ">Figure 10
<p>Predicted HDMAC performance when the optimized packet size is applied: (<b>a</b>) Theoretical channel utilization of HDMAC as a function of the distance between the buoy and SN1. (<b>b</b>) ARQ group size of HDMAC as the distance between the buoy and SN1 increases.</p>
Full article ">Figure 11
<p>The simulation results compared with the theoretical HDMAC performance in <a href="#sec3dot3-electronics-13-04288" class="html-sec">Section 3.3</a>. (<b>a</b>) Channel utilization when the <span class="html-italic">Be</span> is 4000 bits. (<b>b</b>) The group size when <span class="html-italic">Be</span> is 4000 bits. (<b>c</b>) Channel utilization using the optimized packet size. (<b>d</b>) The group size using the optimized packet size.</p>
Full article ">Figure 11 Cont.
<p>The simulation results compared with the theoretical HDMAC performance in <a href="#sec3dot3-electronics-13-04288" class="html-sec">Section 3.3</a>. (<b>a</b>) Channel utilization when the <span class="html-italic">Be</span> is 4000 bits. (<b>b</b>) The group size when <span class="html-italic">Be</span> is 4000 bits. (<b>c</b>) Channel utilization using the optimized packet size. (<b>d</b>) The group size using the optimized packet size.</p>
Full article ">
16 pages, 2116 KiB  
Article
Visibility Graph Investigation of the Shallow Seismicity of Lai Chau Area (Vietnam)
by Luciano Telesca, Anh Tuan Thai, Dinh Trong Cao and Thanh Hai Dang
Entropy 2024, 26(11), 932; https://doi.org/10.3390/e26110932 - 31 Oct 2024
Viewed by 451
Abstract
In this study, the topological properties of the shallow seismicity occurring in the area around the Lai Chau hydropower plant (Vietnam) are investigated by using visibility graph (VG) analysis, a well-known method to convert time series into networks or graphs. The relationship between [...] Read more.
In this study, the topological properties of the shallow seismicity occurring in the area around the Lai Chau hydropower plant (Vietnam) are investigated by using visibility graph (VG) analysis, a well-known method to convert time series into networks or graphs. The relationship between the seismicity and reservoir water level was analyzed using Interlayer Mutual Information (IMI) and the Frobenius norm, both applied to the corresponding VG networks. IMI was used to assess the correlation between the two variables, while the Frobenius norm was employed to estimate the time delay between them. The total seismicity, which resulted in an M0.8 with a b-value of 0.86, is characterized by a kM slope of ≈9.1. Analyzing the variation of the seismological and topological parameters of the seismicity relative to the distance from the center of the Lai Chau reservoir revealed the following features: (1) the b-value fluctuates around a mean value of 1.21 at distances of up to 10–11 km, while, for distances larger than 25–30 km, it tends to the value of 0.86; (2) the maximum IMI between the monthly number of earthquakes and the monthly mean water level occurs at a distance of 9–11 km, showing a distance evolution similar to that of the b-value; (3) at these distances from the center of the reservoir, the time lag between the earthquake monthly counts and the monthly water level mean is 9–10 months; (4) the relationship between the b-value and the kM slope suggests that the kM slope depends on the number of earthquakes within a 22 km radius from the center of the dam. Our study’s findings offer new insights into the complex dynamics of seismicity occurring around reservoirs. Full article
(This article belongs to the Special Issue Time Series Analysis in Earthquake Complex Networks)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of earthquake epicenters during the period of September 2014 to June 2021 (dark circles). Lai Chau reservoir is in light blue. The blue and red building symbols indicate the local seismic stations and national seismic stations, respectively. Geology and identified geologically mapped faults (red) with dip and slip direction as given by [<a href="#B18-entropy-26-00932" class="html-bibr">18</a>]. (Modified from [<a href="#B22-entropy-26-00932" class="html-bibr">22</a>]).</p>
Full article ">Figure 2
<p>(<b>a</b>) Earthquake sequence and the links among the magnitudes defined by the VG. (<b>b</b>) Graph of the links among the nodes defined by the VG.</p>
Full article ">Figure 3
<p>Comparison between two sinusoids, <math display="inline"><semantics> <mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo form="prefix">sin</mo> <mfenced separators="" open="(" close=")"> <mstyle scriptlevel="0" displaystyle="true"> <mfrac> <mrow> <mn>2</mn> <mi>π</mi> <mi>t</mi> </mrow> <mn>10</mn> </mfrac> </mstyle> </mfenced> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <msub> <mi>y</mi> <mi mathvariant="normal">d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo form="prefix">sin</mo> <mfenced separators="" open="(" close=")"> <mstyle scriptlevel="0" displaystyle="true"> <mfrac> <mrow> <mn>2</mn> <mi>π</mi> <mo>(</mo> <mi>t</mi> <mo>−</mo> <mn>3</mn> <mo>)</mo> </mrow> <mn>10</mn> </mfrac> </mstyle> </mfenced> </mrow> </semantics></math> (red), both with period <math display="inline"><semantics> <mrow> <mi>T</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> and the second one delayed by <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Froebenius norm <math display="inline"><semantics> <mrow> <mrow> <mo>∥</mo> </mrow> <msub> <mi>A</mi> <mi>y</mi> </msub> <mo>−</mo> <msub> <mi>A</mi> <msub> <mi>y</mi> <mi>d</mi> </msub> </msub> <msub> <mrow> <mo>∥</mo> </mrow> <mn>2</mn> </msub> </mrow> </semantics></math> between the two sinusoids plotted in <a href="#entropy-26-00932-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 5
<p>Frequency–magnitude distribution of earthquakes occurring at depths of up to 10 km and within a 40 km radius from the center of the dam.</p>
Full article ">Figure 6
<p>Sequence of earthquakes occurring within 40 km of the center of Lai Chau reservoir. The links among the magnitudes are defined by the NVG (<b>a</b>) and HVG (<b>b</b>).</p>
Full article ">Figure 7
<p>Sequences of degree <span class="html-italic">k</span> for the NVG and HVG applied to the seismic dataset shown in <a href="#entropy-26-00932-f006" class="html-fig">Figure 6</a>.</p>
Full article ">Figure 8
<p><span class="html-italic">k-M</span> relationship between the degree and the magnitude for the whole seismic dataset. The slope of the regression line is 9.07.</p>
Full article ">Figure 9
<p>Distribution of the <span class="html-italic">k-M</span> slope for the randomized earthquake sequences. The vertical line indicates the <span class="html-italic">k-M</span> slope of the original sequence.</p>
Full article ">Figure 10
<p>Monthly earthquake counts (red) and mean water level (blue) during the investigation period.</p>
Full article ">Figure 11
<p>Number of earthquakes in the complete seismic catalog (<b>a</b>,<b>b</b>) showing completeness magnitude versus distance from the center of the reservoir.</p>
Full article ">Figure 12
<p>Variation of the <span class="html-italic">b</span>-value with the distance from the center of the dam (blue). The red horizontal line represents the <span class="html-italic">b</span>-value calculated for the tectonic seismicity observed prior to the reservoir impoundment, corresponding to a completeness magnitude of 0.7, as indicated in [<a href="#B34-entropy-26-00932" class="html-bibr">34</a>]. The error on <span class="html-italic">b</span> is indicated by the vertical bars, while the red dotted horizontal lines delimit the error band on <span class="html-italic">b</span> calculated in [<a href="#B34-entropy-26-00932" class="html-bibr">34</a>].</p>
Full article ">Figure 13
<p>IMI between monthly number of earthquakes and monthly mean water level calculated by NVG (red) and HVG (blue).</p>
Full article ">Figure 14
<p>Variation with the time lag <math display="inline"><semantics> <mi>τ</mi> </semantics></math> of the Frobenius norm of the difference between the adjacency matrix of the monthly number of earthquakes and that of the <math display="inline"><semantics> <mi>τ</mi> </semantics></math>-shifted monthly mean of water calculated by using the NVG (<b>a</b>) and the HVG (<b>b</b>).</p>
Full article ">Figure 15
<p>Variation of the <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>−</mo> <mi>M</mi> </mrow> </semantics></math> slope with the distance from the center of the dam.</p>
Full article ">Figure 16
<p>Relationship between the <span class="html-italic">k</span>-<span class="html-italic">M</span> slope and the <span class="html-italic">b</span>-value for the Lai Chau dataset (for distances from the dam center less than, larger than, or equal to 22 km) compared with that of the seismic datasets analyzed in previous studies (Iran [<a href="#B12-entropy-26-00932" class="html-bibr">12</a>], Italy and Taiwan [<a href="#B33-entropy-26-00932" class="html-bibr">33</a>], Mexico [<a href="#B11-entropy-26-00932" class="html-bibr">11</a>], Pannonia [<a href="#B38-entropy-26-00932" class="html-bibr">38</a>], R1 and R2 [<a href="#B39-entropy-26-00932" class="html-bibr">39</a>], and ST2 [<a href="#B37-entropy-26-00932" class="html-bibr">37</a>]).</p>
Full article ">Figure 17
<p>Relationship between the <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>−</mo> <mi>M</mi> </mrow> </semantics></math> slope and the <span class="html-italic">b</span>-value for the seismic datasets extracted at distances from the center varying from 6 to 40 km. The vertical black line separates the values relative to distances less than, greater than, or equal to 22 km from the center of the dam.</p>
Full article ">Figure 18
<p>Relationship between the <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>−</mo> <mi>M</mi> </mrow> </semantics></math> slope and the number of events for the seismic datasets extracted at distances from the center varying from 6 to 40 km. The vertical black line separates the values relative to distances less than, greater than, or equal to 22 km from the center of the dam.</p>
Full article ">Figure 19
<p>Variation of the ratio between <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>−</mo> <mi>M</mi> </mrow> </semantics></math> slope and the <span class="html-italic">b</span>-value with catalog size for Lai Chau dataset (red circles) compared with the relationship proposed in [<a href="#B15-entropy-26-00932" class="html-bibr">15</a>] (blue filled circles).</p>
Full article ">
14 pages, 16241 KiB  
Article
Seismic Random Noise Attenuation Using DARE U-Net
by Tara P. Banjade, Cong Zhou, Hui Chen, Hongxing Li, Juzhi Deng, Feng Zhou and Rajan Adhikari
Remote Sens. 2024, 16(21), 4051; https://doi.org/10.3390/rs16214051 - 30 Oct 2024
Viewed by 546
Abstract
Seismic data processing plays a pivotal role in extracting valuable subsurface information for various geophysical applications. However, seismic records often suffer from inherent random noise, which obscures meaningful geological features and reduces the reliability of interpretations. In recent years, deep learning methodologies have [...] Read more.
Seismic data processing plays a pivotal role in extracting valuable subsurface information for various geophysical applications. However, seismic records often suffer from inherent random noise, which obscures meaningful geological features and reduces the reliability of interpretations. In recent years, deep learning methodologies have shown promising results in performing noise attenuation tasks on seismic data. In this research, we propose modifications to the standard U-Net structure by integrating dense and residual connections, which serve as the foundation of our approach named the dense and residual (DARE U-Net) network. Dense connections enhance the receptive field and ensure that information from different scales is considered during the denoising process. Our model implements local residual connections between layers within the encoder, which allows earlier layers to directly connect with deep layers. This promotes the flow of information, allowing the network to utilize filtered and unfiltered input. The combined network mechanisms preserve the spatial information loss during the contraction process so that the decoder can locate the features more accurately by retaining the high-resolution features, enabling precise location in seismic image denoising. We evaluate this adapted architecture by applying synthetic and real data sets and calculating the peak signal-to-noise ratio (PSNR) and structural similarity index measure (SSIM). The effectiveness of this method is well noted. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>DARE U-Net architecture.</p>
Full article ">Figure 2
<p>Residual connection.</p>
Full article ">Figure 3
<p>Local residual connection within each layer of an encoder.</p>
Full article ">Figure 4
<p>Structure of residual dense block.</p>
Full article ">Figure 5
<p>A sample of the training data. (<b>a</b>) Noise-free data. (<b>b</b>) Noisy data.</p>
Full article ">Figure 6
<p>Test on four sets of seismic data. First to fifth column: noise-free data, noisy data, denoised by wavelet, U-Net, and DARE U-Net.</p>
Full article ">Figure 7
<p>(<b>a</b>) Noise-free data. (<b>b</b>) Noisy data. (<b>c</b>) Denoised by wavelet. (<b>d</b>) Denoised by U-Net. (<b>e</b>) Denoised by DARE U-Net.</p>
Full article ">Figure 8
<p>FK spectrum comparisons. (<b>a</b>) Noise-free data. (<b>b</b>) Noisy data. (<b>c</b>) Denoised by wavelet. (<b>d</b>) Denoised by U-Net. (<b>e</b>) Denoised by DARE U-Net.</p>
Full article ">Figure 8 Cont.
<p>FK spectrum comparisons. (<b>a</b>) Noise-free data. (<b>b</b>) Noisy data. (<b>c</b>) Denoised by wavelet. (<b>d</b>) Denoised by U-Net. (<b>e</b>) Denoised by DARE U-Net.</p>
Full article ">Figure 9
<p>Real data test. (<b>a</b>) Noise-free data. (<b>b</b>) Noisy data. (<b>c</b>) Denoised by wavelet. (<b>d</b>) Denoised by U-Net. (<b>e</b>) Denoised by DARE U-Net.</p>
Full article ">Figure 10
<p>Residual section of denoised real data. (<b>a</b>) Wavelet. (<b>b</b>) U-Net. (<b>c</b>) DARE U-Net.</p>
Full article ">
15 pages, 7689 KiB  
Article
Supervised Deep Learning for Detecting and Locating Passive Seismic Events Recorded with DAS: A Case Study
by Emad Al-Hemyari, Olivia Collet, Konstantin Tertyshnikov and Roman Pevzner
Sensors 2024, 24(21), 6978; https://doi.org/10.3390/s24216978 - 30 Oct 2024
Viewed by 462
Abstract
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 [...] Read more.
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. Full article
Show Figures

Figure 1

Figure 1
<p>A high-level diagram demonstrating the application of supervised deep learning.</p>
Full article ">Figure 2
<p>(<b>a</b>) A diagram of the four wells used to record induced events at the Otway Stage 3 CO<sub>2</sub> injection monitoring site, where green stars mark event locations. (<b>b</b>) Configuration schematics of the fiber-optic cable installations.</p>
Full article ">Figure 3
<p>Upscaled 1D model of (<b>a</b>) P-wave and S-wave velocity and (<b>b</b>) density from vertical CRC-3 well. (<b>c</b>) An extended 2D model over 3000 m of offset with an overlay of source and receiver locations used for modelling, focusing on an area of interest around the induced event locations.</p>
Full article ">Figure 4
<p>Three 1.5 s long strain-rate seismograms with matching geometries of (<b>a</b>) synthetic data for a source located at a depth of 1450 m and an offset of 720 m from CRC-3 well, (<b>b</b>) a strong induced event at an estimated depth of 1470 m and an offset of 720 m from CRC-3 well, and (<b>c</b>) a noise record from CRC-3 well.</p>
Full article ">Figure 5
<p>A schematic demonstrating the application of the deep-learning approach for microseismic event detection, location, and subsurface property estimation using ResNet50 architecture.</p>
Full article ">Figure 6
<p>Training and application data preparation workflow highlighting common steps in blue, steps unique to training data in purple, and steps unique to application data in green.</p>
Full article ">Figure 7
<p>(<b>a</b>) Synthetic events location predictions compared to the ground truth locations used for training. (<b>b</b>) Reconstructed subsurface properties were overlayed on the true models.</p>
Full article ">Figure 8
<p>Predictions of event locations for data from vertical CRC-3 well. Green stars represent the actual locations of induced events relative to the well positioned at zero offset. Stars denoting predictions of induced event locations are in red, and predicted noise locations are in blue. The magenta dots show the actual locations of synthetic data used for the training.</p>
Full article ">Figure 9
<p>The effect of not including time shifts in the training data on (<b>a</b>) event location predictions with time-shift annotations, and (<b>b</b>) the resulting large prediction errors, as compared to the effect of including time shifts in the training data on (<b>c</b>) event location predictions, and (<b>d</b>) the resulting reasonably low prediction errors.</p>
Full article ">Figure 10
<p>The effect of varying synthetics signal-to-noise ratios on detection.</p>
Full article ">Figure 11
<p>Predictions of event and noise locations for data from slightly deviated (<b>a</b>) CRC-4 and (<b>b</b>) CRC-6 wells. Green stars represent the actual locations of induced events relative to the well positioned at zero offset.</p>
Full article ">
Back to TopTop