Xiaolei Tu, Esteban Jeremy Bowles-Martinez, Adam Schultz
{"title":"Massively parallel modeling of electromagnetic field in conductive media: An MPI-CUDA implementation on Multi-GPU computers","authors":"Xiaolei Tu, Esteban Jeremy Bowles-Martinez, Adam Schultz","doi":"10.1016/j.cageo.2024.105710","DOIUrl":"10.1016/j.cageo.2024.105710","url":null,"abstract":"<div><p>Numerical modeling of electromagnetic (EM) fields in a conductive marine environment is crucial for marine EM data interpretation. During marine controlled-source electromagnetic (MCSEM) surveys, a variety of transmitter locations are used to introduce electric currents. The resulting electric and magnetic fields are then concurrently logged by a network of receivers. The forward simulation of MCSEM data for a subsea structure whose electrical properties vary in all three dimensions is computationally intensive. We demonstrate how such computations may be substantially accelerated by adapting algorithms to operate efficiently on modern GPUs with many core architectures. The algorithm we present features a hybrid MPI-CUDA programming model suitable for multi-GPU computers and consists of three levels of parallelism. We design the optimal kernels for different components to minimize redundant memory accesses. We have tested the algorithm on NVIDIA Kepler architecture and achieved up to 105 × speedup compared with the serial code version. We further showcased the algorithm's performance advantages through its application to a realistic marine model featuring complex geological structures. Our algorithm's significant efficiency increase opens the possibility of 3D MCSEM data interpretation based on probabilistic or machine learning approaches, which require tens of thousands of forward simulations for every survey.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105710"},"PeriodicalIF":4.2,"publicationDate":"2024-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142157755","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Combining deep neural network and spatio-temporal clustering to automatically assess rockburst and seismic hazard – Case study from Marcel coal mine in Upper Silesian Basin, Poland","authors":"Adam Lurka","doi":"10.1016/j.cageo.2024.105709","DOIUrl":"10.1016/j.cageo.2024.105709","url":null,"abstract":"<div><p>Mine induced seismic events are a major safety concern in mining and require careful monitoring and management to reduce their effects. Therefore, an essential step in assessing seismic and rock burst hazards is the analysis of mine seismicity. Recently, deep neural networks have been used to automatically determine seismic wave arrival times, surpassing human performance and allowing their use in seismic data analysis such as seismic event location and seismic energy calculation. In order to properly automate the rockburst and seismic hazard assessment deep neural network phase picker and a spatio-temporal clustering method were utilized. Seismic and rockburst hazards were statistically quantified using two-way contingency tables for two categorical variables: seismic energy level of mine tremors and number of clusters. Correlations between several spatio-temporal clusters and a statistical association between two categorical variables: seismic energy level and cluster number indicate an increase of seismic hazard in the Marcel hard coal mine in Poland. A new automated tool has been elaborated to automatically identify high-stress areas in mines in the form of spatio-temporal clusters.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105709"},"PeriodicalIF":4.2,"publicationDate":"2024-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0098300424001924/pdfft?md5=2bf4c7ce15a9d7979aa62ba8147334ed&pid=1-s2.0-S0098300424001924-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142136832","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"DRRGlobal: Uncovering the weak phases from global seismograms using the damped rank-reduction method","authors":"Wei Chen , Yangkang Chen","doi":"10.1016/j.cageo.2024.105687","DOIUrl":"10.1016/j.cageo.2024.105687","url":null,"abstract":"<div><p>Some target seismic signals in the earthquake data can be very weak compared with interfering phases, and are thus difficult to detect, which further hinders the effective usage of these weak phases for subsequent high-resolution imaging of earth interiors. The strong ambient noise makes this situation even more troublesome since the weak signals can be mostly buried in the noise. Here, we present an open-source package for uncovering the weak phases from global seismograms. We adopt a two-step scheme to reconstruct and denoise array data. The first step is weighted average interpolation which puts the data into irregular grids. The second step adopts the weighted projection-onto-convex sets based on damped rank-reduction to further interpolate and denoise for the binned data. Taking the complexity of the weak signal into consideration, we adopt the automatic strategy to select an appropriate rank in different localized windows. We conduct several synthetic tests to carefully investigate the performance regarding effectiveness, robustness, and efficiency, and compare the algorithm with the frequency–wavenumber-domain projection onto convex sets method that is already used in the global seismology literature. Finally, the proposed framework is validated via a recorded array data set of the 1995 May 5 Philippines earthquake.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105687"},"PeriodicalIF":4.2,"publicationDate":"2024-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142050260","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Forward modeling of single-sided magnetic resonance and evaluation of T2 fitting error based on geometric analytical method","authors":"Ruixin Miao, Yunzhi Wang, Qingyue Wang, Yan Zheng, Xiyu He, Chunpeng Ren, Chuandong Jiang","doi":"10.1016/j.cageo.2024.105705","DOIUrl":"10.1016/j.cageo.2024.105705","url":null,"abstract":"<div><p>Single-sided magnetic resonance (SSMR) offers advantages of portability and noninvasive measurement for water detection, with significant potential applications in groundwater exploration, petroleum well logging, and soil moisture monitoring. However, the inherent highly inhomogeneous static magnetic field and radiofrequency (RF) field in SSMR necessitate the utilization of the Carr–Purcell–Meiboom–Gill (CPMG) sequence measurement scheme. To accelerate forward modeling during pulse excitation, we introduce a Geometric Analysis Method (GAM) and assess <span><math><msub><mrow><mi>T</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span> error using its primary parameters. The GAM involves applying spatial geometric rotations on the magnetization vector, leading to an analytical solution to the Bloch equation that disregards relaxation effects. Compared with the rotation matrix (RM) method, the GAM demonstrates high accuracy and reduces computational time by approximately 20.9%. By analyzing the primary parameters governing the magnetization vector in the analytical formula, we evaluated their impact on the transverse relaxation time (<span><math><msub><mrow><mi>T</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span>) obtained through fitting the SE signal. Ultimately, the forward modeling results of the CPMG sequence within the region of interest (ROI) of a single-sided Halbach magnet array are validated. The <span><math><msub><mrow><mi>T</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span> fitting error increases as the primary parameters deviate from the ideal values, highlighting their significant role in the <span><math><msub><mrow><mi>T</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span> fitting results. This study provides a theoretical foundation for optimizing the design of SSMR magnets and RF coils.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105705"},"PeriodicalIF":4.2,"publicationDate":"2024-08-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142044498","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Fukai Zhang , Zhengli Yan , Chao Liu , Haiyan Zhang , Shan Zhao , Jun Liu , Ziqi Zhao
{"title":"Enhanced taxonomic identification of fusulinid fossils through image–text integration using transformer","authors":"Fukai Zhang , Zhengli Yan , Chao Liu , Haiyan Zhang , Shan Zhao , Jun Liu , Ziqi Zhao","doi":"10.1016/j.cageo.2024.105701","DOIUrl":"10.1016/j.cageo.2024.105701","url":null,"abstract":"<div><p>The accurate taxonomic identification of fusulinid fossils holds significant scientific value in palaeontology, paleoecology, and palaeogeography. However, imbalanced image samples lead to the model preferring to learn features from categories with many samples while ignoring fewer sample categories, greatly reducing the prediction accuracy of fusulinid fossil identification. Moreover, the textual description of fusulinid fossils contains rich feature information. We collected and created an order fusulinid multimodal (OFM) dataset for research. We proposed a transformer-based multimodal integration framework (TMIF) using deep learning for fusulinid fossil identification. Compared to traditional neural networks, the transformer can create global dependencies between features at different locations. TMIF incorporates image and text branches dedicated to extracting features for both modalities, and a pivotal cross-modal integration module that allows visual features to learn textual semantic features sufficiently to obtain a more comprehensive feature representation. Experimental evaluation using the OFM dataset shows that TMIF achieves a prediction accuracy of 81.7%, which is a 2.8% improvement over the only image-based method. Further comparative analyses across multiple networks affirm that the TMIF performs optimally in addressing the taxonomic identification of fusulinid fossils with imbalanced samples.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105701"},"PeriodicalIF":4.2,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142006696","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Semantic segmentation of coastal aerial/satellite images using deep learning techniques: An application to coastline detection","authors":"Pietro Scala, Giorgio Manno, Giuseppe Ciraolo","doi":"10.1016/j.cageo.2024.105704","DOIUrl":"10.1016/j.cageo.2024.105704","url":null,"abstract":"<div><p>A new CNN based approach supported by semantic segmentation, was proposed. This approach is frequently used to carry out regional-scale studies. The core of our method revolves around a CNN model, based on the famous U-Net architecture. Its purpose is to identify different classes of pixels on satellite images and later to automatically detect the coastline. The recently launched Coast Train dataset was used to train the CNN model. Traditional coastline detection was improved (“water/land” segmentation) by means of two new aspects the use of the Sobel-edge loss function and the segmentation of the satellite images into several categories like built-up areas, vegetation and land besides beach/sand and water classes. The approach used ensures a more precise coastline extraction, distinguishing water pixels from all other categories. Our model adeptly identifies features, such as cliff vegetation or coastal roads, that some models might overlook. In this way, coastline localization and its drawing for regional scale study, have minor uncertainties. The performance of the CNN-based method, achieving 85% accuracy and 80% IoU (Intersection over Union) in the segmentation process. The ability of the model to extract the coastline was validated on a Sicilian case study, notably the San Leone beach (Agrigento). The model's results align closely with the ground truth, moreover, its reliability was further confirmed when it was tested on other Sicilian coastal regions.</p><p>Beyond robustness, the model offers a promising avenue for enhanced coastal analysis potentially applicable to coastal planning and management.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105704"},"PeriodicalIF":4.2,"publicationDate":"2024-08-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142012036","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Wei Li , Megha Chakraborty , Claudia Quinteros Cartaya , Jonas Köhler , Johannes Faber , Men-Andrin Meier , Georg Rümpker , Nishtha Srivastava
{"title":"SAIPy: A Python package for single-station earthquake monitoring using deep learning","authors":"Wei Li , Megha Chakraborty , Claudia Quinteros Cartaya , Jonas Köhler , Johannes Faber , Men-Andrin Meier , Georg Rümpker , Nishtha Srivastava","doi":"10.1016/j.cageo.2024.105686","DOIUrl":"10.1016/j.cageo.2024.105686","url":null,"abstract":"<div><p>Seismology has witnessed significant advancements in recent years with the application of deep learning methods to address a broad range of problems. These techniques have demonstrated their remarkable ability to effectively extract statistical properties from extensive datasets, surpassing the capabilities of traditional approaches to an extent. In this study, we present SAIPy, an open-source Python package specifically developed for fast seismic data processing by implementing deep learning. SAIPy offers solutions for multiple seismological tasks, including earthquake signal detection, seismic phase picking, first motion polarity identification and magnitude estimation. We introduce upgraded versions of previously published models such as CREIME_RT capable of identifying earthquakes with an accuracy above 99.8% and a root mean squared error of 0.38 unit in magnitude estimation. These upgraded models outperform state-of-the-art approaches like the Vision Transformer network. SAIPy provides an API that simplifies the integration of these advanced models, including CREIME_RT, DynaPicker_v2, and PolarCAP, along with benchmark datasets. It also, to the best of our knowledge, introduces the first fully automated deep learning based pipeline to process continuous waveforms. The package has the potential to be used for real-time earthquake monitoring to enable timely actions to mitigate the impact of seismic events. Ongoing development efforts aim to further enhance SAIPy’s performance and incorporate additional features that enhance exploration efforts, and it also would be interesting to approach the retraining of the whole package as a multi-task learning problem. A detailed description of all functions is available in a supplementary document.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105686"},"PeriodicalIF":4.2,"publicationDate":"2024-08-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0098300424001699/pdfft?md5=181c42ee7372a7ceb6bfb0f6134f713e&pid=1-s2.0-S0098300424001699-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142040326","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Geo-Hgan: Unsupervised anomaly detection in geochemical data via latent space learning","authors":"Liang Ding , Bainian Chen , Yuelong Zhu , Hai Dong , Guiyang Chan , Pengcheng Zhang","doi":"10.1016/j.cageo.2024.105703","DOIUrl":"10.1016/j.cageo.2024.105703","url":null,"abstract":"<div><p>Reconstructing geochemical data for anomaly detection using Generative Adversarial Networks (GANs) has become a prevalent method in identifying geochemical anomalies. However, injecting random noise into GANs can induce model instability. To mitigate this issue, we propose a novel anomaly detection model, Geo-Hgan, which integrates a dual adversarial network architecture with a Latent Space Adversarial Module (LSAM) to learn the distribution of latent variables from arbitrary data and optimize the sample reconstruction process, thereby alleviating instability during GAN training. Additionally, an encoder guided by the LSAM-pretrained GAN is employed to extract variational features, facilitating rapid and effective sample mapping into the latent space defined by LSAM. Experimental results demonstrate that under unsupervised conditions, Geo-Hgan achieves an Area Under the Curve (AUC) score of 85% across three geochemical datasets, outperforming similar models in accuracy and reconstruction capabilities. To assess its versatility and generalization ability, we extend Geo-Hgan to anomaly detection tasks in computer vision, where it achieves an average AUC score of 98.7% on the MvtecAD dataset, setting a new state-of-the-art performance in the domain. Furthermore, we propose AnomFilter, a method for setting anomaly thresholds based on the clustering hypothesis. AnomFilter identifies high-confidence anomaly samples identified by Geo-Hgan in the source domain and iteratively transfers them to the target domain. These high-confidence anomaly samples, combined with a small number of known positive samples in the target domain, enhance the accuracy of supervised geochemical anomaly detection in the target domain, which achieved an AUC score of 94%. The utilization of anomaly detection models for sample transfer learning offers a novel perspective for future work.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105703"},"PeriodicalIF":4.2,"publicationDate":"2024-08-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142012133","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
An-Bo Li , Shi-Yu Xu , Xian-Yu Liu , Guo-Nian Lü , Xian-Li Xie , Matthew Fox
{"title":"Relative dating of fault activity using the principle of cross-cutting relationships: An automated approach","authors":"An-Bo Li , Shi-Yu Xu , Xian-Yu Liu , Guo-Nian Lü , Xian-Li Xie , Matthew Fox","doi":"10.1016/j.cageo.2024.105702","DOIUrl":"10.1016/j.cageo.2024.105702","url":null,"abstract":"<div><p>Fault dating plays an essential role in understanding deformation histories and modeling the tectonic evolution of orogenic belts. However, direct fault dating methods via different isotope geochronological techniques are expensive, and their use is often limited in many cases, making it essential to develop a fast and low-cost fault relative dating method. Therefore, on the basis of knowledge graphs and knowledge reasoning technology, this study proposes an automatic method to relatively date periods of fault activity using the principle of cross-cutting relationships between faults and strata. The method mainly involves (1) generating the knowledge graph based on a digital geological map; (2) using the knowledge reasoning algorithm to interpret the cross-cutting relationships amongst faults and generating the temporal sequence of fault activity; (3) relative dating the faults based on the cross-cutting relationships between faults and strata; and (4) according to the temporal sequence of fault activity, the relationship between faults can be revealed, and relative dating can be optimized. Results for cases in western Nevada and Qixia Hill of Nanjing illustrate the effectiveness of this method for interpreting the period of fault activity. The accuracy rates of the recognition results in the two cases were 90.24% and 80.77%, respectively, which means that the proposed method has the potential to relatively date fault activity across large areas. The algorithm is an effective supplement to the existing direct method of fault dating. The algorithm can efficiently infer the development sequence and the age of fault activity based on geological maps and geological cross-sections, which is of great significance for understanding regional tectonic history, evaluating earthquake disasters, and modeling tectonic evolution processes.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105702"},"PeriodicalIF":4.2,"publicationDate":"2024-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142002356","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xinpeng Ma , Yunhe Liu , Changchun Yin , Jingru Li , Jun Li , Xiuyan Ren , Shiwen Li
{"title":"Global mantle conductivity imaging using 3-D geomagnetic depth sounding with real earth surface conductivity constraint","authors":"Xinpeng Ma , Yunhe Liu , Changchun Yin , Jingru Li , Jun Li , Xiuyan Ren , Shiwen Li","doi":"10.1016/j.cageo.2024.105697","DOIUrl":"10.1016/j.cageo.2024.105697","url":null,"abstract":"<div><p>The water content in the Earth's interior is of great significance for material circulation and the dynamic evolution of the planet. The water content in mantle minerals significantly affects their conductivities. By measuring the variations in conductivity within the Earth, we can infer the water content in the mantle and study the movement and processes of materials within the Earth. The geomagnetic depth sounding is a widely used method for imaging the mantle conductivity as it has large sounding depth. However, the ocean induction effects can seriously impact geomagnetic data that can't be well corrected using conventional methods. Here, we present a novel three-dimensional inversion method for geomagnetic depth sounding to overcome the ocean induction effects by directly adopt the real earth surface conductivity into the inverse model. In this method, the unstructured tetrahedral grids are used to represent the model in multi-scale and the vector finite-element method is adopted to accurately compute the geomagnetic responses. The synthetic model tests show that the earth surface conductivity has serious effect on the inversion results, but it can be well suppressed by directly modeling it in the inverse model. We further invert the data from 128 geomagnetic stations around the world and obtain a more accurate new model of global mantle conductivity.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"192 ","pages":"Article 105697"},"PeriodicalIF":4.2,"publicationDate":"2024-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141978533","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}