Computers & Geosciences最新文献_第3页

Novel empirical curvelet denoising strategy for suppressing mixed noise of microseismic data 抑制微地震数据混合噪声的新经验小曲线去噪策略

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-07 DOI: 10.1016/j.cageo.2024.105751

Liyuan Feng , Binhong Li , Huailiang Li , Jian He

{"title":"Novel empirical curvelet denoising strategy for suppressing mixed noise of microseismic data","authors":"Liyuan Feng , Binhong Li , Huailiang Li , Jian He","doi":"10.1016/j.cageo.2024.105751","DOIUrl":"10.1016/j.cageo.2024.105751","url":null,"abstract":"<div><div>We present a novel denoising strategy based on empirical curvelet transform (ECT) for noisy microseismic data. Our approach can simultaneously suppress high-frequency, low-frequency, and shared-bandwidth noises and preserve detailed information on the noisy microseismic data. Initially, we design a new threshold estimation method by adding a scale factor for ECT threshold denoising. Subsequently, we construct an adaptive parameter model employing the similarity standard deviation for the non-local means (NLM) algorithm. Then, we divide the coefficients obtained from the ECT decomposition into two sets based on the energy spectrum, subjecting each set to improved adaptive thresholding and improved NLM denoising algorithms. Eventually, we reconstruct the denoised signals using the empirical curvelet inverse transform. Our results demonstrate that under a signal-to-noise ratio (SNR) of <span><math><mo>−</mo></math></span>10 dB, the proposed strategy achieves a correlation coefficient of 0.9524, a root mean square error of 0.198, an SNR of 1.36 dB, and reduces the first arrival picking error to 0.00382 s. Furthermore, application on the real microseismic data further confirms that the proposed method can clarify the corresponding first arrival.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105751"},"PeriodicalIF":4.2,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659016","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}

引用次数: 0

New fast imaging techniques for electrical source transient electromagnetic data: Approaches and application 电源瞬态电磁数据的新型快速成像技术：方法与应用

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-07 DOI: 10.1016/j.cageo.2024.105770

Yu-lian Zhu , Wei-ying Chen , Wan-ting Song , Si-xu Han

{"title":"New fast imaging techniques for electrical source transient electromagnetic data: Approaches and application","authors":"Yu-lian Zhu , Wei-ying Chen , Wan-ting Song , Si-xu Han","doi":"10.1016/j.cageo.2024.105770","DOIUrl":"10.1016/j.cageo.2024.105770","url":null,"abstract":"<div><div>The rapid imaging of electrical source transient electromagnetic (TEM) data involves two essential processes: the calculation of apparent resistivity and the conversion of time to depth. Traditionally, the definition of full-time apparent resistivity is defined by considering solely the vertical magnetic field, which is predicated on the monotonic relationship between the resistivity and the electromagnetic field response. Based on the concept of peak time, we have developed distinct methodologies for calculating the apparent resistivity for both the horizontal electric field (<em>e</em><sub>x</sub>) and the vertical induced voltage (<em>v</em><sub>z</sub>), which demonstrated accuracy across the entire time range examined. We also introduced a formula to address discrepancies in apparent resistivity arising from the non-dipole size effect of the source, thereby ensuring that the algorithm can adapt to any transmitting and receiving configuration. Furthermore, we provided straightforward and precise time-depth conversion equations applicable to both <em>e</em><sub><em>x</em></sub> and <em>v</em><sub><em>z</em></sub>, which facilitate the rapid imaging of observational data. Multiple numerical examples were employed to illustrate the effectiveness and robustness of this approach. Finally, we applied this imaging technique to the data processing of actual measured data from a survey area conducted in Ningxia Province, and the imaging results accurately reflected the distribution of the electrical structure of the subsurface strata. The innovative imaging technique presented in this study holds considerable potential for the expedited processing and analysis of ground-based and semi-aerial electrical source transient electromagnetic survey data, which are widely employed in contemporary applications.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105770"},"PeriodicalIF":4.2,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659030","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}

引用次数: 0

Digital rock physics: Calculation of effective elastic properties of heterogeneous materials using graphical processing units (GPUs) 数字岩石物理学：利用图形处理器（GPU）计算异质材料的有效弹性特性

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-06 DOI: 10.1016/j.cageo.2024.105749

Yury Alkhimenkov

{"title":"Digital rock physics: Calculation of effective elastic properties of heterogeneous materials using graphical processing units (GPUs)","authors":"Yury Alkhimenkov","doi":"10.1016/j.cageo.2024.105749","DOIUrl":"10.1016/j.cageo.2024.105749","url":null,"abstract":"<div><div>An application based on graphical processing units (GPUs) applied to 3-D digital images is described for computing the linear anisotropic elastic properties of heterogeneous materials. The application can also retrieve the property contribution tensors of individual inclusions of any shape. The code can be executed on professional GPUs as well as on a basic laptop or personal computer Nvidia GPUs. The application is extremely fast: a calculation of the effective elastic properties of volumes consisting of about 7 million voxel elements (191<sup>3</sup>) takes less than 4 s of computational time using a single A100 GPU; 3 min for 100 million voxel elements (479<sup>3</sup>) using a single A100 GPU; 14 min for 350 million voxel elements (703<sup>3</sup>) using a single A100 GPU. Several comparisons against analytical solutions are provided. In addition, an evaluation of the anisotropic effective elastic properties of a 3-D digital image of a cracked Carrara marble sample is presented. The software can be downloaded from a permanent repository Zenodo, the link with a doi is given in the manuscript.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105749"},"PeriodicalIF":4.2,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659019","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}

引用次数: 0

Backpropagation-based inference for spatial interpolation to estimate the blastability index in an open pit mine 基于反向传播推理的空间插值法估算露天矿的可爆性指数

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-06 DOI: 10.1016/j.cageo.2024.105756

Yakin Hajlaoui , Richard Labib , Jean-François Plante , Michel Gamache

引用次数: 0

Million-core scalable 3D anisotropic reverse time migration on the Sugon exascale supercomputer 在曙光超大规模超级计算机上实现百万核级可扩展三维各向异性反向时间迁移

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-05 DOI: 10.1016/j.cageo.2024.105754

Sihai Wu , Jiubing Cheng , Jianwei Ma , Tengfei Wang , Xueshan Yong , Yang Ji

{"title":"Million-core scalable 3D anisotropic reverse time migration on the Sugon exascale supercomputer","authors":"Sihai Wu , Jiubing Cheng , Jianwei Ma , Tengfei Wang , Xueshan Yong , Yang Ji","doi":"10.1016/j.cageo.2024.105754","DOIUrl":"10.1016/j.cageo.2024.105754","url":null,"abstract":"<div><div>Reverse time migration (RTM) plays a crucial role in high-resolution seismic imaging of the Earth’s interior. However, scaling it across millions of cores in parallel to process large-scale seismic datasets poses significant computational challenges, because the conventional storage solutions are insufficient to deal with the I/O and memory bottlenecks. To address this issue, we present a highly scalable 3D RTM algorithm for vertically transverse isotropic (VTI) media, optimized for the Sugon exascale supercomputer, utilizing over 1,024,000 cores with optimal weak-scaling efficiency. Through cache optimizations tailored for the new deep computing unit (DCU) accelerator architecture, our approach achieves a maximum speedup of 6x compared to conventional methods on a single accelerator. Moreover, based on the lossy compression and boundary-saving techniques, we reduce storage requirements by 266 times, which allows for the effective utilization of million-core computing resources and ensures scalability efficiency when handling large-scale datasets for complex geophysical tasks. Finally, when applied to a industrial dataset, the method demonstrates robust scalability and high efficiency, making it well-suited for large-scale geophysical exploration.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105754"},"PeriodicalIF":4.2,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659028","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}

引用次数: 0

Geothermal modeling in complex geological systems with ComPASS 利用 ComPASS 进行复杂地质系统的地热建模

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-05 DOI: 10.1016/j.cageo.2024.105752

A. Armandine Les Landes , L. Beaude , D. Castanon Quiroz , L. Jeannin , S. Lopez , F. Smai , T. Guillon , R. Masson

{"title":"Geothermal modeling in complex geological systems with ComPASS","authors":"A. Armandine Les Landes , L. Beaude , D. Castanon Quiroz , L. Jeannin , S. Lopez , F. Smai , T. Guillon , R. Masson","doi":"10.1016/j.cageo.2024.105752","DOIUrl":"10.1016/j.cageo.2024.105752","url":null,"abstract":"<div><div>In deep geothermal reservoirs, faults and fractures play a major role, serving as regulators of fluid flow and heat transfer while also providing feed zones for production wells. To accurately model the operation of geothermal fields, it is necessary to explicitly consider objects of varying spatial scales, from the reservoir scale itself, to that of faults and fractures, down to the scale of the injection and production wells.</div><div>Our main objective in developing the ComPASS geothermal flow simulator, was to take into account all of these geometric constraints in a flow and heat transfer numerical model using generic unstructured meshes. In its current state, the code provides a parallel implementation of a spatio-temporal discretization of the non-linear equations driving compositional multi-phase thermal flows in porous fractured media on unstructured meshes. It allows an explicit discretization of faults and fractures as 2D hybrid objects, embedded in a 3D matrix. Similarly, wells are modeled as one dimensional graphs discretized by edges of the 3D mesh which allows arbitrary multi-branch wells. The resulting approach is particularly flexible and robust in terms of modeling.</div><div>Its practical interest is demonstrated by two case studies in high-energy geothermal contexts.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105752"},"PeriodicalIF":4.2,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659029","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}

引用次数: 0

Regularization by double complementary priors for full waveform inversion 全波形反演的双互补先验规范化

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-02 DOI: 10.1016/j.cageo.2024.105753

Hongyu Qi , Zhenwu Fu , Yang Li , Bo Han , Longsuo Li

{"title":"Regularization by double complementary priors for full waveform inversion","authors":"Hongyu Qi , Zhenwu Fu , Yang Li , Bo Han , Longsuo Li","doi":"10.1016/j.cageo.2024.105753","DOIUrl":"10.1016/j.cageo.2024.105753","url":null,"abstract":"<div><div>Full-waveform inversion (FWI) represents an advanced geophysical imaging technique focused on intricately depicting subsurface physical properties by iteratively minimizing the differences between the simulated and observed seismograms. Unfortunately, the conventional FWI utilizing a least-squares loss function suffers from various drawbacks, including the challenge of local minima and the necessity for human intervention in parameter fine-tuning. It is particularly problematic when handling noisy data and inadequate initial models. Recent works have exhibited promising performance in two-dimensional FWI by integrating structural sparse representation to procure adaptive dictionaries. Drawing inspiration from the competitiveness of structural sparse representation, we introduce a paradigm of group sparse residuals that integrates two types of complementary prior information by harnessing both the internal and external subsurface media models. The proposed algorithm is based on an alternate minimization algorithm to guarantee workflow flexibility and efficient optimization capabilities. We experimentally validate our method for two baseline geological models, and a comparison of the results demonstrates that the proposed algorithm faithfully recovers the velocity models and consistently outperforms other traditional or learning-based algorithms. A further benefit from the group sparse coding used in this method is that it reduces the sensitivity to data noise.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105753"},"PeriodicalIF":4.2,"publicationDate":"2024-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659017","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}

引用次数: 0

Latent diffusion models for parameterization of facies-based geomodels and their use in data assimilation 用于基于面的地质模型参数化的潜在扩散模型及其在数据同化中的应用

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-02 DOI: 10.1016/j.cageo.2024.105755

Guido Di Federico, Louis J. Durlofsky

{"title":"Latent diffusion models for parameterization of facies-based geomodels and their use in data assimilation","authors":"Guido Di Federico, Louis J. Durlofsky","doi":"10.1016/j.cageo.2024.105755","DOIUrl":"10.1016/j.cageo.2024.105755","url":null,"abstract":"<div><div>Geological parameterization entails the representation of a geomodel using a small set of latent variables and a mapping from these variables to grid-block properties such as porosity and permeability. Parameterization is useful for data assimilation (history matching), as it maintains geological realism while reducing the number of variables to be determined. Diffusion models are a new class of generative deep-learning procedures that have been shown to outperform previous methods, such as generative adversarial networks, for image generation tasks. Diffusion models are trained to “denoise”, which enables them to generate new geological realizations from input fields characterized by random noise. Latent diffusion models, which are the specific variant considered in this study, provide dimension reduction through use of a low-dimensional latent variable. The model developed in this work includes a variational autoencoder for dimension reduction and a U-net for the denoising process. Our application involves conditional 2D three-facies (channel-levee-mud) systems. The latent diffusion model is shown to provide realizations that are visually consistent with samples from geomodeling software. Quantitative metrics involving spatial and flow-response statistics are evaluated, and general agreement between the diffusion-generated models and reference realizations is observed. Stability tests are performed to assess the smoothness of the parameterization method. The latent diffusion model is then used for ensemble-based data assimilation. Two synthetic “true” models are considered. Significant uncertainty reduction, posterior P<sub>10</sub>–P<sub>90</sub> forecasts that generally bracket observed data, and consistent posterior geomodels, are achieved in both cases.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105755"},"PeriodicalIF":4.2,"publicationDate":"2024-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142659015","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}

引用次数: 0

Rock-type classification: A (critical) machine-learning perspective 岩石类型分类：机器学习的（关键）视角

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-01 DOI: 10.1016/j.cageo.2024.105730

Pedro Ribeiro Mendes Júnior , Soroor Salavati , Oscar Linares , Maiara Moreira Gonçalves , Marcelo Ferreira Zampieri , Vitor Hugo de Sousa Ferreira , Manuel Castro , Rafael de Oliveira Werneck , Renato Moura , Elayne Morais , Ahmed Esmin , Leopoldo Lusquino Filho , Denis José Schiozer , Alexandre Ferreira , Alessandra Davólio , Anderson Rocha

{"title":"Rock-type classification: A (critical) machine-learning perspective","authors":"Pedro Ribeiro Mendes Júnior , Soroor Salavati , Oscar Linares , Maiara Moreira Gonçalves , Marcelo Ferreira Zampieri , Vitor Hugo de Sousa Ferreira , Manuel Castro , Rafael de Oliveira Werneck , Renato Moura , Elayne Morais , Ahmed Esmin , Leopoldo Lusquino Filho , Denis José Schiozer , Alexandre Ferreira , Alessandra Davólio , Anderson Rocha","doi":"10.1016/j.cageo.2024.105730","DOIUrl":"10.1016/j.cageo.2024.105730","url":null,"abstract":"<div><div>We investigate machine-learning techniques for rock-type classification. A throughout literature review (considering the machine-learning technique, number of classes, rock types, and image types) presents a diversity of datasets employed and a wide range of classification results as well as multiple problem formulations. Throughout the discussion of the literature, we highlight some common machine-learning pitfalls and criticize the decisions taken by some authors on the problem formulation. We present an experimental contribution by evaluating the classification of seven types of rocks found in carbonate reservoirs along with state-of-the-art Convolutional Neural Networks (CNNs) architectures available through a well-known open-source library. For this experimentation, we detail the preparation of the dataset of drill core plugs (DCPs), the experimental setup itself, and the obtained results considering the normalized accuracy and the traditional accuracy as metrics. We performed the manual background segmentation of the employed dataset of DCPs; so the results reported are not influenced by the background of the images. We evaluate top-1, top-2, and top-3 performance for the problem. We apply fusion of multiple CNNs for richer classification decisions. We also contribute by presenting the manual classification — human labeling by looking at the image on the computer screen — of the same seven-class dataset, performed by six non-geologist volunteers. Finally, we present a conclusion for the results obtained with our experiments and share valuable advice for researchers applying machine learning to rock classification.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"193 ","pages":"Article 105730"},"PeriodicalIF":4.2,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142552516","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}

引用次数: 0

A hybrid inversion algorithm to obtain the resistivity of the uninvaded zone based on the array induction log 一种基于阵列感应测井的混合反演算法，用于获取未侵蚀区的电阻率

IF 4.2 2区地球科学

Computers & Geosciences Pub Date : 2024-11-01 DOI: 10.1016/j.cageo.2024.105766

Xinmin Ge , Mohmmed Ishag , Haiyan Li , Jundong Liu , Cuixia Qu , Badreldein Mohamed

{"title":"A hybrid inversion algorithm to obtain the resistivity of the uninvaded zone based on the array induction log","authors":"Xinmin Ge , Mohmmed Ishag , Haiyan Li , Jundong Liu , Cuixia Qu , Badreldein Mohamed","doi":"10.1016/j.cageo.2024.105766","DOIUrl":"10.1016/j.cageo.2024.105766","url":null,"abstract":"<div><div>This study investigates the impact of the drilling mud invasion on the borehole-measured resistivity. The primary objective is to retrieve the true resistivity of the formation, which helps in identifying different fluids in the reservoir. To achieve this goal, We proposed a hybrid inversion approach integrating the Levenberg-Marquardt and Markov Chain Monte Carlo algorithms with a five-parameter formation resistivity model. Synthetic and real-world data are utilized to assess the method's robustness and reliability. The simulated result indicated that the method is reliable when the data noise level is less than 5%.</div><div>The method applied to real-world data revealed that the resistivity profile on the water zone showed a slight increase in the inverted resistivity from measured resistivity. Meanwhile, in the oil zone, the calculated resistivity revealed a high deviation from the measured resistivity, indicating the effects of mud invasion. The introduced methods are only applicable when the invasions of mud occur within the range of the logging tool's depth of investigation. Moreover, the method may give no reliable result when the invasion exceeds the tool's investigation depth. It indicates its limitation.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"194 ","pages":"Article 105766"},"PeriodicalIF":4.2,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142593369","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}

引用次数: 0