首页 > 最新文献

Computers & Geosciences最新文献

英文 中文
A study on grille structure modeling algorithm for fault-controlled fractured-cavity reservoirs: A case study of the shunbei no. 5 fault zone 断控缝洞型油藏格栅结构建模算法研究——以顺北油田为例。5断裂带
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-12 DOI: 10.1016/j.cageo.2025.105988
Haojie Shang , Shuyang Chen , Yunfeng He , Lixin Wang , Yanshu Yin , Pengfei Xie
Establishing high-resolution 3D geological models of fault-controlled reservoirs is crucial for optimizing well placement and development plan. Carbonate fault-controlled fracture-cavity reservoirs in the Shunbei area of the Tarim Basin, Northwest China, exhibit complex heterogeneity. These reservoirs typically comprise fracture planes, caves and disordered bodies. Fracture planes are narrow and banded, with caves and disordered bodies distributed around them. Within fracture planes and caves, crush belts and bedrock belts alternate to form grille structures. These pose significant challenges for traditional modeling algorithms to characterize accurately. To address this, we proposed a hierarchical object-based modeling algorithm to reproduce fault-controlled fracture-cavity body's grille structural trends and shapes. Using seismic data from the Shunbei No.5 Fault Zone (with a resolution of 25∗25m) and well logging data (including drilling fluid loss data and resistivity logging data), conduct research on grille structure modeling algorithms. First, the fault-controlled fracture-cavity reservoirs are distinguished by fracture planes, caves, and disordered bodies, and contour models are established via seismic attributes threshold truncation. Second, statistics on the scale of development of crush belts and breccia belts under 100 m of fracture planes and caves in different stress sections by logging data. A regional growth tracking algorithm are applied to identify fracture planes trend lines, which can be classified into single, multi, convergent, and branching forms based on contour characteristics. Third, cumulative probability sampling is used to determine the number and scale of the crush and breccia belts. Grille structure models were constructed at three levels: bedrock, crush, and breccia belts. Results indicate successful identification of trend lines matching the structural contours, establishing accurate grille structure models by employing hierarchical simulation strategy under trend line constraints. The models established by traditional methods exhibit significant randomness, making it difficult to control both the variable developmental trajectories of individual belts and the relative positional relationships among multiple belts. Based on these geological facies models, corresponding physical property models were generated, achieving high accuracy in reserve calculations and numerical simulations with less than 10 % error, thus providing valuable guidance for oil and gas development. In the future, more compatible contour models can be established through methods like multi-attribute fusion and deep learning. By integrating production data, the developmental positions and connectivity of grille belts can be constrained.
建立断控油藏的高分辨率三维地质模型对于优化井位和开发方案至关重要。塔里木盆地顺北地区碳酸盐岩断控缝洞型储层具有复杂的非均质性。这些储层通常由裂缝面、溶洞和无序体组成。断裂面狭窄,呈带状,周围分布有溶洞和无序体。在裂隙面和溶洞内,破碎带和基岩带交替形成格栅结构。这些对传统建模算法的准确表征提出了重大挑战。为了解决这一问题,我们提出了一种基于分层对象的建模算法来再现断层控制缝洞体的格栅结构趋势和形状。利用顺北五断裂带地震资料(分辨率为25 * 25m)和测井资料(包括钻井液漏失和电阻率测井资料),开展格栅结构建模算法研究。首先,通过裂缝面、洞室和无序体对断控缝洞型储层进行区分,并通过地震属性阈值截断建立等高线模型;其次,利用测井资料统计不同应力剖面100 m裂缝面和溶洞下破碎带和角砾岩带的发育规模。采用区域增长跟踪算法识别裂缝面趋势线,根据裂缝面趋势线的轮廓特征将裂缝面趋势线划分为单一、多、收敛和分支形式。第三,采用累积概率抽样法确定破碎带和角砾岩带的数量和规模。在基岩带、破碎带和角砾岩带三个层次建立了格栅结构模型。结果表明,在趋势线约束下,采用分层仿真策略,成功识别出与结构轮廓相匹配的趋势线,建立了精确的格栅结构模型。传统方法建立的模型具有显著的随机性,难以控制单个带的可变发展轨迹和多个带之间的相对位置关系。根据这些地质相模型建立相应的物性模型,实现了储量计算和数值模拟的高精度,误差小于10%,为油气开发提供了有价值的指导。未来,可以通过多属性融合和深度学习等方法建立更加兼容的轮廓模型。通过整合生产数据,可以约束格栅带的发展位置和连通性。
{"title":"A study on grille structure modeling algorithm for fault-controlled fractured-cavity reservoirs: A case study of the shunbei no. 5 fault zone","authors":"Haojie Shang ,&nbsp;Shuyang Chen ,&nbsp;Yunfeng He ,&nbsp;Lixin Wang ,&nbsp;Yanshu Yin ,&nbsp;Pengfei Xie","doi":"10.1016/j.cageo.2025.105988","DOIUrl":"10.1016/j.cageo.2025.105988","url":null,"abstract":"<div><div>Establishing high-resolution 3D geological models of fault-controlled reservoirs is crucial for optimizing well placement and development plan. Carbonate fault-controlled fracture-cavity reservoirs in the Shunbei area of the Tarim Basin, Northwest China, exhibit complex heterogeneity. These reservoirs typically comprise fracture planes, caves and disordered bodies. Fracture planes are narrow and banded, with caves and disordered bodies distributed around them. Within fracture planes and caves, crush belts and bedrock belts alternate to form grille structures. These pose significant challenges for traditional modeling algorithms to characterize accurately. To address this, we proposed a hierarchical object-based modeling algorithm to reproduce fault-controlled fracture-cavity body's grille structural trends and shapes. Using seismic data from the Shunbei No.5 Fault Zone (with a resolution of 25∗25m) and well logging data (including drilling fluid loss data and resistivity logging data), conduct research on grille structure modeling algorithms. First, the fault-controlled fracture-cavity reservoirs are distinguished by fracture planes, caves, and disordered bodies, and contour models are established via seismic attributes threshold truncation. Second, statistics on the scale of development of crush belts and breccia belts under 100 m of fracture planes and caves in different stress sections by logging data. A regional growth tracking algorithm are applied to identify fracture planes trend lines, which can be classified into single, multi, convergent, and branching forms based on contour characteristics. Third, cumulative probability sampling is used to determine the number and scale of the crush and breccia belts. Grille structure models were constructed at three levels: bedrock, crush, and breccia belts. Results indicate successful identification of trend lines matching the structural contours, establishing accurate grille structure models by employing hierarchical simulation strategy under trend line constraints. The models established by traditional methods exhibit significant randomness, making it difficult to control both the variable developmental trajectories of individual belts and the relative positional relationships among multiple belts. Based on these geological facies models, corresponding physical property models were generated, achieving high accuracy in reserve calculations and numerical simulations with less than 10 % error, thus providing valuable guidance for oil and gas development. In the future, more compatible contour models can be established through methods like multi-attribute fusion and deep learning. By integrating production data, the developmental positions and connectivity of grille belts can be constrained.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105988"},"PeriodicalIF":4.2,"publicationDate":"2025-06-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144291288","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
Deep learning assisted denoising of fast polychromatic X-ray micro-CT imaging of multiphase flow in porous media 深度学习辅助多孔介质中多相流快速多色x射线微ct成像去噪
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-11 DOI: 10.1016/j.cageo.2025.105990
E.S. Mathew , S.J. Jackson , D. Wildenschild , P. Mostaghimi , K. Tang , R.T. Armstrong
Understanding the flow of fluids in the subsurface and their interaction with different solid surfaces is crucial for addressing challenging geological applications such as CO2 sequestration, enhanced oil recovery, and environmental remediation of polluted aquifers. Synchrotron-based 3D X-ray micro-computed tomography (micro-CT) has enabled the visualization of dynamic pore-filling events in multiphase flow experiments at sub-second time resolutions. However, the limited accessibility of synchrotron facilities has driven the use of low-flux polychromatic micro-CT systems, which often produce relatively noisy images during fast scans. To overcome this limitation, we propose a deep learning workflow using a cycleGAN network trained on unpaired datasets as no direct pixel-wise correspondence exists between the noisy domain and the high-quality domain. This approach transforms noisy fast polychromatic micro-CT scans into high-quality images, enabling detailed analysis of multiphase flow dynamics. The effectiveness of the denoising process was verified using blind image quality evaluators and Minkowski functionals for the non-wetting phases. The results indicate that the cycleGAN network achieves an average 1 to 6 percentage error difference for 3D morphological analysis parameters and outperforms other filtering methods such as non-local means and the adaptive Weiner filter, demonstrating its potential as a reliable technique for restoring noisy fast scans from polychromatic micro-CT systems.
了解地下流体的流动及其与不同固体表面的相互作用对于解决具有挑战性的地质应用至关重要,例如二氧化碳封存、提高石油采收率和污染含水层的环境修复。基于同步加速器的三维x射线微计算机断层扫描(micro-CT)能够以亚秒级的时间分辨率可视化多相流实验中的动态孔隙填充事件。然而,同步加速器设施的有限可及性促使了低通量多色微型ct系统的使用,这些系统在快速扫描期间经常产生相对嘈杂的图像。为了克服这一限制,我们提出了一种深度学习工作流,使用在未配对数据集上训练的cycleGAN网络,因为噪声域和高质量域之间不存在直接的像素对应关系。该方法可将嘈杂的快速多色微ct扫描图像转换为高质量图像,从而实现多相流动力学的详细分析。使用盲图像质量评估器和非润湿阶段的Minkowski函数验证了去噪过程的有效性。结果表明,cycleGAN网络在3D形态分析参数上实现了平均1到6个百分点的误差差异,并且优于其他滤波方法,如非局部均值和自适应Weiner滤波器,证明了它作为一种可靠的技术,可以从多色微ct系统中恢复有噪声的快速扫描。
{"title":"Deep learning assisted denoising of fast polychromatic X-ray micro-CT imaging of multiphase flow in porous media","authors":"E.S. Mathew ,&nbsp;S.J. Jackson ,&nbsp;D. Wildenschild ,&nbsp;P. Mostaghimi ,&nbsp;K. Tang ,&nbsp;R.T. Armstrong","doi":"10.1016/j.cageo.2025.105990","DOIUrl":"10.1016/j.cageo.2025.105990","url":null,"abstract":"<div><div>Understanding the flow of fluids in the subsurface and their interaction with different solid surfaces is crucial for addressing challenging geological applications such as CO<sub>2</sub> sequestration, enhanced oil recovery, and environmental remediation of polluted aquifers. Synchrotron-based 3D X-ray micro-computed tomography (micro-CT) has enabled the visualization of dynamic pore-filling events in multiphase flow experiments at sub-second time resolutions. However, the limited accessibility of synchrotron facilities has driven the use of low-flux polychromatic micro-CT systems, which often produce relatively noisy images during fast scans. To overcome this limitation, we propose a deep learning workflow using a cycleGAN network trained on unpaired datasets as no direct pixel-wise correspondence exists between the noisy domain and the high-quality domain. This approach transforms noisy fast polychromatic micro-CT scans into high-quality images, enabling detailed analysis of multiphase flow dynamics. The effectiveness of the denoising process was verified using blind image quality evaluators and Minkowski functionals for the non-wetting phases. The results indicate that the cycleGAN network achieves an average 1 to 6 percentage error difference for 3D morphological analysis parameters and outperforms other filtering methods such as non-local means and the adaptive Weiner filter, demonstrating its potential as a reliable technique for restoring noisy fast scans from polychromatic micro-CT systems.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105990"},"PeriodicalIF":4.2,"publicationDate":"2025-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144338653","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
Spatio-temporal deep networks with feature disentangling for advancing earthquake monitoring 基于特征解缠的时空深度网络推进地震监测
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-10 DOI: 10.1016/j.cageo.2025.105974
Fanchun Meng , Tao Ren , Xinyu He , Zhuoran Dong , Xinyue Wang , Hongfeng Chen
Earthquake monitoring is essential for assessing seismic hazards and involves interconnected tasks such as phase picking and location estimation. Existing single-parameter estimation methods suffer from error accumulation caused by task interdependencies and typically rely on empirical values. Multi-parameter estimation methods often depend on data from multiple stations, posing challenges in modeling and revealing the inter-station relationships. To address these challenges, this study proposes a novel neural network, SINE, designed to simultaneously estimate key parameters in earthquake monitoring, including P-phase arrival time, location, and magnitude. SINE develops a multi-task framework that incorporates Graph Neural Networks (GNNs) and Bidirectional Long Short-Term Memory networks (BI-LSTM) to extract spatio-temporal features, effectively mitigating error accumulation across the tasks. Unlike previous GNN-based models, SINE incorporates a feature disentanglement structure to automatically identify multiple potential relationships between seismic stations. Additionally, the CNN-based parsing unit is employed to regress multiple seismic parameters simultaneously. Evaluation on datasets from Southern California and Italy shows that SINE outperforms existing DL models and traditional seismological methods. Furthermore, SINE effectively reduces inter-task dependencies, enhancing robustness in earthquake monitoring.
地震监测是评估地震灾害的关键,它涉及相位选择和位置估计等相互关联的任务。现有的单参数估计方法存在任务相互依赖导致的误差累积,并且通常依赖于经验值。多参数估计方法往往依赖于来自多个站点的数据,这给建模和揭示站间关系带来了挑战。为了应对这些挑战,本研究提出了一种新的神经网络——正弦正弦神经网络,旨在同时估计地震监测中的关键参数,包括p相到达时间、位置和震级。SINE开发了一个多任务框架,该框架结合了图神经网络(gnn)和双向长短期记忆网络(BI-LSTM)来提取时空特征,有效地减少了任务间的错误积累。与之前基于gnn的模型不同,SINE采用了特征解纠缠结构,可以自动识别地震台站之间的多个潜在关系。此外,采用基于cnn的解析单元对多个地震参数同时进行回归。对南加州和意大利数据集的评估表明,SINE优于现有的DL模型和传统的地震学方法。此外,正弦函数有效地降低了任务间依赖性,增强了地震监测的鲁棒性。
{"title":"Spatio-temporal deep networks with feature disentangling for advancing earthquake monitoring","authors":"Fanchun Meng ,&nbsp;Tao Ren ,&nbsp;Xinyu He ,&nbsp;Zhuoran Dong ,&nbsp;Xinyue Wang ,&nbsp;Hongfeng Chen","doi":"10.1016/j.cageo.2025.105974","DOIUrl":"10.1016/j.cageo.2025.105974","url":null,"abstract":"<div><div>Earthquake monitoring is essential for assessing seismic hazards and involves interconnected tasks such as phase picking and location estimation. Existing single-parameter estimation methods suffer from error accumulation caused by task interdependencies and typically rely on empirical values. Multi-parameter estimation methods often depend on data from multiple stations, posing challenges in modeling and revealing the inter-station relationships. To address these challenges, this study proposes a novel neural network, SINE, designed to simultaneously estimate key parameters in earthquake monitoring, including P-phase arrival time, location, and magnitude. SINE develops a multi-task framework that incorporates Graph Neural Networks (GNNs) and Bidirectional Long Short-Term Memory networks (BI-LSTM) to extract spatio-temporal features, effectively mitigating error accumulation across the tasks. Unlike previous GNN-based models, SINE incorporates a feature disentanglement structure to automatically identify multiple potential relationships between seismic stations. Additionally, the CNN-based parsing unit is employed to regress multiple seismic parameters simultaneously. Evaluation on datasets from Southern California and Italy shows that SINE outperforms existing DL models and traditional seismological methods. Furthermore, SINE effectively reduces inter-task dependencies, enhancing robustness in earthquake monitoring.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105974"},"PeriodicalIF":4.2,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144291286","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
Controlled-source electromagnetic signal detection using a hybrid deep learning model: Convolutional and long short-term memory neural networks 使用混合深度学习模型的可控源电磁信号检测:卷积和长短期记忆神经网络
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-07 DOI: 10.1016/j.cageo.2025.105986
Yecheng Liu , Diquan Li , Jin Li , Yanfang Hu , Zijie Liu , Xian Zhang
Controlled-source electromagnetic (CSEM) method using a periodic transmitted signal source suppresses random noise by superimposing and averaging the recorded signal over multiple periods. However, it still faces great challenges in strong interference environments. Here we present a CSEM signal detection method based on a hybrid architecture of convolutional neural network (CNN) and long short-term memory neural network (LSTM). Our method improves the quality of the recorded signal by filtering out all the noise periods and retaining the useful signal periods. The core lies in combining the spatial feature extraction capability of CNN with the time series modeling capability of LSTM, which can deeply excavate the feature differences between noise and CSEM useful signal. Meanwhile, a classification model of signal and noise is constructed using a large-scale training dataset. This hybrid model exhibits superior performance compared to other deep learning models such as CNN or LSTM. Also, we propose a novel signal detection mechanism that not only maintains the periodicity of CSEM signal, but also greatly enhances processing efficiency. The results of synthetic and measured data demonstrate that our method can obtain useful signals from CSEM data containing different noise types and significantly improve the quality of sounding curves. In particular, our method is useful for CSEM signal detection in strong interference.
控制源电磁(CSEM)方法采用周期发射信号源,通过对记录信号在多个周期内的叠加和平均来抑制随机噪声。然而,在强干扰环境下,它仍然面临着很大的挑战。本文提出了一种基于卷积神经网络(CNN)和长短期记忆神经网络(LSTM)混合架构的CSEM信号检测方法。该方法滤除所有噪声周期,保留有用的信号周期,提高了记录信号的质量。其核心在于将CNN的空间特征提取能力与LSTM的时间序列建模能力相结合,能够深度挖掘噪声与CSEM有用信号之间的特征差异。同时,利用大规模训练数据集构建了信号和噪声的分类模型。与CNN或LSTM等其他深度学习模型相比,该混合模型表现出优越的性能。此外,我们还提出了一种新的信号检测机制,既保持了CSEM信号的周期性,又大大提高了处理效率。合成数据和实测数据的结果表明,该方法可以从不同噪声类型的CSEM数据中获得有用的信号,显著提高了测深曲线的质量。该方法特别适用于强干扰条件下的CSEM信号检测。
{"title":"Controlled-source electromagnetic signal detection using a hybrid deep learning model: Convolutional and long short-term memory neural networks","authors":"Yecheng Liu ,&nbsp;Diquan Li ,&nbsp;Jin Li ,&nbsp;Yanfang Hu ,&nbsp;Zijie Liu ,&nbsp;Xian Zhang","doi":"10.1016/j.cageo.2025.105986","DOIUrl":"10.1016/j.cageo.2025.105986","url":null,"abstract":"<div><div>Controlled-source electromagnetic (CSEM) method using a periodic transmitted signal source suppresses random noise by superimposing and averaging the recorded signal over multiple periods. However, it still faces great challenges in strong interference environments. Here we present a CSEM signal detection method based on a hybrid architecture of convolutional neural network (CNN) and long short-term memory neural network (LSTM). Our method improves the quality of the recorded signal by filtering out all the noise periods and retaining the useful signal periods. The core lies in combining the spatial feature extraction capability of CNN with the time series modeling capability of LSTM, which can deeply excavate the feature differences between noise and CSEM useful signal. Meanwhile, a classification model of signal and noise is constructed using a large-scale training dataset. This hybrid model exhibits superior performance compared to other deep learning models such as CNN or LSTM. Also, we propose a novel signal detection mechanism that not only maintains the periodicity of CSEM signal, but also greatly enhances processing efficiency. The results of synthetic and measured data demonstrate that our method can obtain useful signals from CSEM data containing different noise types and significantly improve the quality of sounding curves. In particular, our method is useful for CSEM signal detection in strong interference.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105986"},"PeriodicalIF":4.2,"publicationDate":"2025-06-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144254094","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 novel algorithm and software for efficient global gravimetric forward modeling in the spherical coordinate system 一种在球坐标系下实现全球重力正演模拟的新算法和软件
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-07 DOI: 10.1016/j.cageo.2025.105985
Wenjin Chen , Xiaoyu Tang , Robert Tenzer
We present a novel algorithm and complementary software for the global gravimetric forward modeling that accommodates masses with complex shapes and density distributions defined in the frame of spherical coordinates. Traditional gravimetric forward modeling techniques often face significant challenges when dealing with irregularly shaped bodies and complex density variations, leading to high computational costs and long processing times. To address these practical limitations, we introduce an innovative algorithm that divides the 3-D mass-density body into spherical concentric rings with equal intervals in the radial direction, while discretizing the latitudinal and longitudinal directions into a grid with equal intervals. After discretization, we assume that each spherical volumetric mass ring has a laterally varying density and constant upper and lower bounds. Based on this approach, the gravitational field at any point outside the Earth is evaluated as the sum of the gravitational contributions generated by each concentric ring. This discretization allows applying the Fast Fourier Transform (FFT) technique to drastically improve computational efficiency of the spherical harmonic analysis and synthesis. Numerical results are validated against corresponding solutions obtained using the tesseroid method in the spatial domain. The comparison of results reassures a high accuracy of proposed method, with relative differences between results obtained from both methods less than 1 % and 4 % for the gravitational attraction and gradient respectively, while significantly improving the numerical efficiency. When modelling the gravitational field quantities of very complex structures by means of their geometry and density distribution, such as the Earth's crustal density structure, the numerical efficiency improved several orders of magnitude.
我们提出了一种新的全球重力正演建模算法和补充软件,该算法和软件可以适应在球坐标框架中定义的具有复杂形状和密度分布的质量。传统的重力正演建模技术在处理不规则形体和复杂密度变化时往往面临重大挑战,导致计算成本高,处理时间长。为了解决这些实际限制,我们引入了一种创新的算法,该算法将三维质量密度体在径向上划分为等间隔的球形同心圆,同时将纬度和纵向离散为等间隔的网格。离散化后,我们假设每个球形体积质量环具有横向变化的密度和恒定的上下边界。基于这种方法,地球外任何一点的引力场被评估为每个同心圆产生的引力贡献的总和。这种离散化允许应用快速傅里叶变换(FFT)技术大大提高球谐分析和合成的计算效率。数值结果与曲面法在空间域中得到的相应解进行了验证。结果对比表明,该方法具有较高的精度,两种方法的引力和梯度计算结果的相对差异分别小于1%和4%,同时显著提高了计算效率。在对非常复杂的结构,如地球的地壳密度结构,利用其几何和密度分布来模拟引力场量时,数值效率提高了几个数量级。
{"title":"A novel algorithm and software for efficient global gravimetric forward modeling in the spherical coordinate system","authors":"Wenjin Chen ,&nbsp;Xiaoyu Tang ,&nbsp;Robert Tenzer","doi":"10.1016/j.cageo.2025.105985","DOIUrl":"10.1016/j.cageo.2025.105985","url":null,"abstract":"<div><div>We present a novel algorithm and complementary software for the global gravimetric forward modeling that accommodates masses with complex shapes and density distributions defined in the frame of spherical coordinates. Traditional gravimetric forward modeling techniques often face significant challenges when dealing with irregularly shaped bodies and complex density variations, leading to high computational costs and long processing times. To address these practical limitations, we introduce an innovative algorithm that divides the 3-D mass-density body into spherical concentric rings with equal intervals in the radial direction, while discretizing the latitudinal and longitudinal directions into a grid with equal intervals. After discretization, we assume that each spherical volumetric mass ring has a laterally varying density and constant upper and lower bounds. Based on this approach, the gravitational field at any point outside the Earth is evaluated as the sum of the gravitational contributions generated by each concentric ring. This discretization allows applying the Fast Fourier Transform (FFT) technique to drastically improve computational efficiency of the spherical harmonic analysis and synthesis. Numerical results are validated against corresponding solutions obtained using the tesseroid method in the spatial domain. The comparison of results reassures a high accuracy of proposed method, with relative differences between results obtained from both methods less than 1 % and 4 % for the gravitational attraction and gradient respectively, while significantly improving the numerical efficiency. When modelling the gravitational field quantities of very complex structures by means of their geometry and density distribution, such as the Earth's crustal density structure, the numerical efficiency improved several orders of magnitude.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105985"},"PeriodicalIF":4.2,"publicationDate":"2025-06-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144291287","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
An innovative adaptive mineral segmentation method via augmentation and fusion of single and orthogonal polarized images: AMS-p/xpl 基于单偏振和正交偏振图像增强和融合的自适应矿物分割方法:AMS-p/xpl
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-06 DOI: 10.1016/j.cageo.2025.105987
W. Ma , P. Lin , S. Li , Z.H. Xu
An adaptive mineral segmentation method is proposed to address the challenges of traditional thin section identification, which is typically expert-dependent, highly subjective, and time-consuming. The method is based on the augmentation and fusion of single polarized (PPL) and orthogonal polarized (XPL) images using an enhanced Deeplabv3-based segmentation model. The Depth Separable Convolution (DSC) is introduced to strengthen edge and texture features from PPL images, and a ColorBoost module is designed to enhance color information from XPL images. An adaptive feature fusion mechanism is employed to integrate complementary polarized features and dynamically adjust their contribution weights. The results demonstrate that the proposed model achieved the highest segmentation performance on the test set, with a mean intersection over union (mIoU) of 89.0 % and an accuracy of 96.7 %. Compared to widely used semantic segmentation networks such as FCN, it demonstrates a notable improvement in mIoU, with a maximum gain of 32.9 %. Additionally, through the integration of feature augmentation and an adaptive fusion mechanism, the model outperforms the baseline DeepLabv3 by 5 % in mIoU. The proposed method provides a more efficient and automated solution for thin section mineral identification, reducing reliance on expert knowledge and improving applicability in practical and non-specialist settings.
针对传统薄片识别方法中专家依赖性强、主观性强、耗时长等问题,提出了一种自适应矿物分割方法。该方法基于基于增强的deeplabv3分割模型对单极化(PPL)和正交极化(XPL)图像进行增强和融合。引入深度可分离卷积(DSC)增强PPL图像的边缘和纹理特征,设计ColorBoost模块增强XPL图像的颜色信息。采用自适应特征融合机制对互补极化特征进行融合,并动态调整互补极化特征的权重。结果表明,该模型在测试集上获得了最高的分割性能,平均相交比(mIoU)为89.0%,准确率为96.7%。与FCN等广泛使用的语义分割网络相比,它在mIoU方面有了显著的改进,最大增益为32.9%。此外,通过集成特征增强和自适应融合机制,该模型的mIoU比基线DeepLabv3高出5%。该方法为薄层矿物识别提供了一种更有效和自动化的解决方案,减少了对专家知识的依赖,提高了在实际和非专业环境中的适用性。
{"title":"An innovative adaptive mineral segmentation method via augmentation and fusion of single and orthogonal polarized images: AMS-p/xpl","authors":"W. Ma ,&nbsp;P. Lin ,&nbsp;S. Li ,&nbsp;Z.H. Xu","doi":"10.1016/j.cageo.2025.105987","DOIUrl":"10.1016/j.cageo.2025.105987","url":null,"abstract":"<div><div>An adaptive mineral segmentation method is proposed to address the challenges of traditional thin section identification, which is typically expert-dependent, highly subjective, and time-consuming. The method is based on the augmentation and fusion of single polarized (PPL) and orthogonal polarized (XPL) images using an enhanced Deeplabv3-based segmentation model. The Depth Separable Convolution (DSC) is introduced to strengthen edge and texture features from PPL images, and a ColorBoost module is designed to enhance color information from XPL images. An adaptive feature fusion mechanism is employed to integrate complementary polarized features and dynamically adjust their contribution weights. The results demonstrate that the proposed model achieved the highest segmentation performance on the test set, with a mean intersection over union (<em>mIoU</em>) of 89.0 % and an accuracy of 96.7 %. Compared to widely used semantic segmentation networks such as FCN, it demonstrates a notable improvement in <em>mIoU</em>, with a maximum gain of 32.9 %. Additionally, through the integration of feature augmentation and an adaptive fusion mechanism, the model outperforms the baseline DeepLabv3 by 5 % in <em>mIoU.</em> The proposed method provides a more efficient and automated solution for thin section mineral identification, reducing reliance on expert knowledge and improving applicability in practical and non-specialist settings.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105987"},"PeriodicalIF":4.2,"publicationDate":"2025-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144262984","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 gradient-optimized least-squares reverse time migration based on the safe type-I anderson acceleration 基于安全i型安德森加速度的梯度优化最小二乘逆时偏移
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-06-05 DOI: 10.1016/j.cageo.2025.105984
Yingming Qu , Chongpeng Huang
Least-squares reverse time migration (LSRTM) can generate preferable images for complex media, but faces substantial computational challenges in field data applications, especially in 3D cases. Many optimization algorithms have been proposed to alleviate this problem, such as the limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method, the restarted generalized minimal residual method, and Anderson acceleration (AA). AA is a popular gradient optimization algorithm that has been widely used in many fields due to its ability to greatly accelerate the convergence of fixed-point iterations and considerably reduce the computational cost. According to Broyden's method, AA is divided into type-I AA (AA-I) and type-II AA (AA-II), with most implementations favoring AA-II due to residual oscillation issues observed in AA-I during data residual minimization. To address the residual vibration issue of AA-I and expedite the convergence of LSRTM, we apply a safe AA-I method to LSRTM, incorporating Powell-type regularization, re-start checking, and safe guarding steps. The Powell-type regularization guarantees the non-singularity of AA-I, while the re-start checking preserves its strong linear independence, both contributing to the stability of AA-I. The safe guarding steps examine the data residual reduction and accelerate the convergence. Our analysis reveals that the optimal step length for the safe AA-I method is approximately 5 times or 10 times the initial steepest descent (SD) iteration. We also derive an exponential scaling law for the safe AA-I step length. In addition, the safe AA-I has faster data residual convergence speed, less computational cost, and higher quality images than SD, conjugate gradient (CG), AA-II, and LBFGS. The safe AA-I is approximately twice as efficient as the LBFGS. Field validation through land seismic data processing shows that LSRTM based on the safe AA-I delivers enhanced structural resolution with sharper imaging events and improved stratigraphic continuity relative to LBFGS-based implementations.
最小二乘逆时偏移(LSRTM)可以为复杂介质生成更好的图像,但在现场数据应用中面临着巨大的计算挑战,特别是在3D情况下。为了解决这一问题,人们提出了许多优化算法,如有限记忆Broyden-Fletcher-Goldfarb-Shanno (LBFGS)法、重新启动广义最小残差法和Anderson加速(AA)法。AA是一种流行的梯度优化算法,由于其能够大大加快不动点迭代的收敛速度,大大降低计算成本,在许多领域得到了广泛的应用。根据Broyden的方法,AA分为i型AA (AA- i)和ii型AA (AA- ii),由于AA- i在数据残差最小化过程中观察到残差振荡问题,大多数实现倾向于AA- ii。为了解决AA-I的残余振动问题,加快LSRTM的收敛速度,我们将一种安全的AA-I方法应用于LSRTM,包括鲍威尔式正则化、重新启动检查和安全防护步骤。powell型正则化保证了AA-I的非奇异性,而重新启动检验保持了AA-I的强线性无关性,都有助于AA-I的稳定性。安全防护步骤检查数据残差减少,加快收敛速度。我们的分析表明,安全的AA-I方法的最佳步长大约是初始最陡下降(SD)迭代的5倍或10倍。我们还推导了安全AA-I步长的指数标度定律。此外,与SD、共轭梯度(CG)、AA-II和LBFGS相比,安全的AA-I具有更快的数据残差收敛速度、更少的计算成本和更高的图像质量。安全的AA-I的效率大约是LBFGS的两倍。通过陆地地震数据处理的现场验证表明,与基于lbfgs的方法相比,基于安全AA-I的LSRTM可以提高构造分辨率,成像事件更清晰,地层连续性更好。
{"title":"A gradient-optimized least-squares reverse time migration based on the safe type-I anderson acceleration","authors":"Yingming Qu ,&nbsp;Chongpeng Huang","doi":"10.1016/j.cageo.2025.105984","DOIUrl":"10.1016/j.cageo.2025.105984","url":null,"abstract":"<div><div>Least-squares reverse time migration (LSRTM) can generate preferable images for complex media, but faces substantial computational challenges in field data applications, especially in 3D cases. Many optimization algorithms have been proposed to alleviate this problem, such as the limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method, the restarted generalized minimal residual method, and Anderson acceleration (AA). AA is a popular gradient optimization algorithm that has been widely used in many fields due to its ability to greatly accelerate the convergence of fixed-point iterations and considerably reduce the computational cost. According to Broyden's method, AA is divided into type-I AA (AA-I) and type-II AA (AA-II), with most implementations favoring AA-II due to residual oscillation issues observed in AA-I during data residual minimization. To address the residual vibration issue of AA-I and expedite the convergence of LSRTM, we apply a safe AA-I method to LSRTM, incorporating Powell-type regularization, re-start checking, and safe guarding steps. The Powell-type regularization guarantees the non-singularity of AA-I, while the re-start checking preserves its strong linear independence, both contributing to the stability of AA-I. The safe guarding steps examine the data residual reduction and accelerate the convergence. Our analysis reveals that the optimal step length for the safe AA-I method is approximately 5 times or 10 times the initial steepest descent (SD) iteration. We also derive an exponential scaling law for the safe AA-I step length. In addition, the safe AA-I has faster data residual convergence speed, less computational cost, and higher quality images than SD, conjugate gradient (CG), AA-II, and LBFGS. The safe AA-I is approximately twice as efficient as the LBFGS. Field validation through land seismic data processing shows that LSRTM based on the safe AA-I delivers enhanced structural resolution with sharper imaging events and improved stratigraphic continuity relative to LBFGS-based implementations.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105984"},"PeriodicalIF":4.2,"publicationDate":"2025-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144262986","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
Lithological mapping using Spatially Constrained Bayesian Network (SCB-Net): A deep learning model for generating field-data-constrained predictions with uncertainty evaluation using remote sensing data 使用空间约束贝叶斯网络(SCB-Net)的岩性制图:一种深度学习模型,用于生成现场数据约束预测,并使用遥感数据进行不确定性评估
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-05-30 DOI: 10.1016/j.cageo.2025.105964
Victor Silva dos Santos , Erwan Gloaguen , Shiva Tirdad
Geological maps are an important source of information for the Earth sciences. These maps are created using numerical or conceptual models that use geological observations to extrapolate data. Geostatistical techniques have traditionally been used to generate reliable predictions that take into account the spatial patterns inherent in the data. However, as the number of auxiliary variables increases, these methods become more labor-intensive. Additionally, traditional machine learning methods often struggle to capture spatial context and extract valuable non-linear information from geoscientific datasets. To address these challenges, we developed the Spatially Constrained Bayesian Network (SCB-Net)—an architecture designed to effectively integrate auxiliary variables while generating spatially constrained predictions. SCB-Net employs a late-fusion strategy, processing auxiliary data and ground-truth information through two parallel encoding paths. Additionally, it leverages Monte Carlo dropout as a Bayesian approximation to quantify model uncertainty. The SCB-Net has been tested on two real-world datasets from northern Quebec, Canada, demonstrating its effectiveness in generating field-data-constrained lithological maps while providing uncertainty estimates for unsampled locations. Our method outperformed the Attention U-Net – a widely used model in image segmentation – by at least 4.7% in accuracy across all tested datasets.
地质图是地球科学的重要资料来源。这些地图是使用数值或概念模型创建的,这些模型使用地质观测来推断数据。传统上使用地质统计技术来产生考虑到数据中固有的空间模式的可靠预测。然而,随着辅助变量数量的增加,这些方法变得更加劳动密集型。此外,传统的机器学习方法往往难以捕获空间背景并从地球科学数据集中提取有价值的非线性信息。为了应对这些挑战,我们开发了空间约束贝叶斯网络(SCB-Net),这是一种架构,旨在有效地集成辅助变量,同时生成空间约束预测。SCB-Net采用后期融合策略,通过两条并行编码路径处理辅助数据和真值信息。此外,它利用蒙特卡罗dropout作为贝叶斯近似来量化模型的不确定性。SCB-Net已经在加拿大魁北克北部的两个真实数据集上进行了测试,证明了它在生成现场数据受限的岩性图方面的有效性,同时为未采样的位置提供了不确定性估计。我们的方法在所有测试数据集上的准确率至少比注意力U-Net(一种广泛使用的图像分割模型)高出4.7%。
{"title":"Lithological mapping using Spatially Constrained Bayesian Network (SCB-Net): A deep learning model for generating field-data-constrained predictions with uncertainty evaluation using remote sensing data","authors":"Victor Silva dos Santos ,&nbsp;Erwan Gloaguen ,&nbsp;Shiva Tirdad","doi":"10.1016/j.cageo.2025.105964","DOIUrl":"10.1016/j.cageo.2025.105964","url":null,"abstract":"<div><div>Geological maps are an important source of information for the Earth sciences. These maps are created using numerical or conceptual models that use geological observations to extrapolate data. Geostatistical techniques have traditionally been used to generate reliable predictions that take into account the spatial patterns inherent in the data. However, as the number of auxiliary variables increases, these methods become more labor-intensive. Additionally, traditional machine learning methods often struggle to capture spatial context and extract valuable non-linear information from geoscientific datasets. To address these challenges, we developed the Spatially Constrained Bayesian Network (SCB-Net)—an architecture designed to effectively integrate auxiliary variables while generating spatially constrained predictions. SCB-Net employs a late-fusion strategy, processing auxiliary data and ground-truth information through two parallel encoding paths. Additionally, it leverages Monte Carlo dropout as a Bayesian approximation to quantify model uncertainty. The SCB-Net has been tested on two real-world datasets from northern Quebec, Canada, demonstrating its effectiveness in generating field-data-constrained lithological maps while providing uncertainty estimates for unsampled locations. Our method outperformed the Attention U-Net – a widely used model in image segmentation – by at least 4.7% in accuracy across all tested datasets.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105964"},"PeriodicalIF":4.2,"publicationDate":"2025-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144222039","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
AI-based geological subsurface reconstruction using sparse convolutional autoencoders 基于稀疏卷积自编码器的人工智能地质地下重建
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-05-29 DOI: 10.1016/j.cageo.2025.105981
Rodrigo Uribe-Ventura , Yoan Barriga-Berrios , Jorge Barriga-Gamarra , Patrice Baby , Willem Viveen
Subsurface reconstruction is critical for geological modeling and resource exploration. Conventional spatial interpolation methods are limited by stationarity and spatial isotropy assumptions, while advanced geostatistical techniques require specialized datasets. Deep learning approaches often need large datasets, which is impractical for geoscientific applications. This study presents an AI-based methodology using a sparse convolutional autoencoder for robust subsurface modeling under data constraints and integrating secondary data sources such as Vertical Electrical Sounding (VES) data. A four-stage testing framework was implemented: (1) emulating conventional interpolation for baseline performance; (2) reconstructing subsurface geometries from synthetic data; (3) incorporating geophysical constraints through VES forward modeling; and (4) validating the methodology using a real-world case study from the Huancayo tectonic basin in the Peruvian Andes, using 41 VES measurements across two cross-sections (12 and 14 km long). Results demonstrate that the proposed model effectively emulates kriging interpolation (mean squared error: 1.5 × 103 to 1.2 × 103 with 100–800 training examples) through transfer learning from an inverse-distance, pre-trained model. In subsurface reconstruction, the model outperforms kriging (37.4–61.7 % improvement across 1–15 % sampling densities) through its ability to adapt to non-stationary conditions. When incorporating synthetic VES data, the model effectively reconstructed subsurface geometries with error reduction from 4.1 × 101 to 9.1 × 103 as stations increased from 1 to 40, demonstrating diminishing returns beyond this point. Application to the Huancayo basin case study validated the model's practical applicability by successfully identifying previously unmapped features including the contact between basement and sedimentary infill, folds and faults. The methodology demonstrates the AI's capability to enhance geological understanding in complex tectonic settings, revealing subtle features and refining existing assumptions about subsurface architecture.
地下重建是地质建模和资源勘探的关键。传统的空间插值方法受到平稳性和空间各向同性假设的限制,而先进的地质统计学技术需要专门的数据集。深度学习方法通常需要大型数据集,这对于地球科学应用来说是不切实际的。本研究提出了一种基于人工智能的方法,使用稀疏卷积自编码器在数据约束下进行鲁棒地下建模,并集成了垂直电测深(VES)数据等辅助数据源。实现了一个四阶段的测试框架:(1)模拟常规插值的基线性能;(2)利用合成数据重建地下几何形状;(3)通过地震探测正演模拟纳入地球物理约束;(4)利用秘鲁安第斯山脉万卡约构造盆地的实际案例研究验证了该方法,该研究使用了两个横截面(12和14公里长)的41个VES测量值。结果表明,该模型通过反距离预训练模型的迁移学习,有效地模拟了克里格插值(均方误差为1.5 × 10−3至1.2 × 10−3,训练样例为100-800)。在地下重建中,该模型通过适应非平稳条件的能力优于克里格法(在1 - 15%的采样密度下提高37.4 - 61.7%)。当结合合成的VES数据时,该模型有效地重建了地下几何形状,当站点从1个增加到40个时,误差从4.1 × 10−1减少到9.1 × 10−3,显示出超过该点的收益递减。在Huancayo盆地的应用验证了该模型的实用性,成功地识别了基底与沉积充填体之间的接触、褶皱和断层等以前未映射的特征。该方法展示了人工智能在复杂构造环境中增强地质认识的能力,揭示了微妙的特征,并完善了关于地下结构的现有假设。
{"title":"AI-based geological subsurface reconstruction using sparse convolutional autoencoders","authors":"Rodrigo Uribe-Ventura ,&nbsp;Yoan Barriga-Berrios ,&nbsp;Jorge Barriga-Gamarra ,&nbsp;Patrice Baby ,&nbsp;Willem Viveen","doi":"10.1016/j.cageo.2025.105981","DOIUrl":"10.1016/j.cageo.2025.105981","url":null,"abstract":"<div><div>Subsurface reconstruction is critical for geological modeling and resource exploration. Conventional spatial interpolation methods are limited by stationarity and spatial isotropy assumptions, while advanced geostatistical techniques require specialized datasets. Deep learning approaches often need large datasets, which is impractical for geoscientific applications. This study presents an AI-based methodology using a sparse convolutional autoencoder for robust subsurface modeling under data constraints and integrating secondary data sources such as Vertical Electrical Sounding (VES) data. A four-stage testing framework was implemented: (1) emulating conventional interpolation for baseline performance; (2) reconstructing subsurface geometries from synthetic data; (3) incorporating geophysical constraints through VES forward modeling; and (4) validating the methodology using a real-world case study from the Huancayo tectonic basin in the Peruvian Andes, using 41 VES measurements across two cross-sections (12 and 14 km long). Results demonstrate that the proposed model effectively emulates kriging interpolation (mean squared error: 1.5 <span><math><mrow><mo>×</mo></mrow></math></span> 10<span><math><mrow><mo>−</mo><mn>3</mn></mrow></math></span> to 1.2 <span><math><mrow><mo>×</mo></mrow></math></span> 10<span><math><mrow><mo>−</mo><mn>3</mn></mrow></math></span> with 100–800 training examples) through transfer learning from an inverse-distance, pre-trained model. In subsurface reconstruction, the model outperforms kriging (37.4–61.7 % improvement across 1–15 % sampling densities) through its ability to adapt to non-stationary conditions. When incorporating synthetic VES data, the model effectively reconstructed subsurface geometries with error reduction from 4.1 <span><math><mrow><mo>×</mo></mrow></math></span> 10<span><math><mrow><mo>−</mo><mn>1</mn></mrow></math></span> to 9.1 <span><math><mrow><mo>×</mo></mrow></math></span> 10<span><math><mrow><mo>−</mo><mn>3</mn></mrow></math></span> as stations increased from 1 to 40, demonstrating diminishing returns beyond this point. Application to the Huancayo basin case study validated the model's practical applicability by successfully identifying previously unmapped features including the contact between basement and sedimentary infill, folds and faults. The methodology demonstrates the AI's capability to enhance geological understanding in complex tectonic settings, revealing subtle features and refining existing assumptions about subsurface architecture.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"204 ","pages":"Article 105981"},"PeriodicalIF":4.2,"publicationDate":"2025-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144178633","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
Quadrangular adaptive mesh for elastic wave simulation in smooth anisotropic media 光滑各向异性介质弹性波模拟的四边形自适应网格
IF 4.2 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-05-23 DOI: 10.1016/j.cageo.2025.105946
Marius Rapenne , Paul Cupillard , Guillaume Caumon , Corentin Gouache
Smooth anisotropic media are often met when implementing effective medium theory, full waveform inversion or seismic imaging. However, computational overburden is often a recurring problem when working with high frequencies or when quantifying uncertainties. In this context, adaptive meshes constitute, in principle, an attractive representation to maximize simulation accuracy while minimizing the computational cost. However, such meshes are difficult to create in the context of smooth anisotropic media as the optimal local size of the elements is not clearly defined. In this work, we present a two-step algorithm to efficiently mesh these media for spectral element method (SEM) simulation in the 2D elastic case. Our algorithm yields quadrangular only meshes which adapt the size of the element to the local and directional S-wave velocity. It relies on a quadtree division introduced by Maréchal (2009) to divide the mesh until the size of each element edge is adapted to the local minimum wavelength that will be propagated. Then, a Laplacian smoothing is applied to further optimize the size of the elements, increasing the global time step and makes the SEM simulation faster while keeping a good accuracy and even improving it in some cases. An application of our method on a 2D section of the homogenized Groningen area shows that simulation time can be reduced by a factor up to 7.
在进行有效介质理论、全波形反演或地震成像时,往往会遇到光滑的各向异性介质。然而,在处理高频或量化不确定性时,计算覆盖常常是一个反复出现的问题。在这种情况下,自适应网格在原则上构成了一种有吸引力的表示,以最大化模拟精度,同时最小化计算成本。然而,在光滑的各向异性介质中,由于单元的最佳局部尺寸没有明确定义,这种网格很难创建。在这项工作中,我们提出了一种两步算法来有效地网格化这些介质,以便在二维弹性情况下进行谱元法(SEM)模拟。我们的算法只产生四边形网格,使单元的大小适应局部和定向s波速度。它依靠marsamchal(2009)引入的四叉树划分来划分网格,直到每个元素边缘的大小适应将被传播的局部最小波长。然后,应用拉普拉斯平滑进一步优化单元的大小,增加全局时间步长,使SEM模拟速度更快,同时保持良好的精度,甚至在某些情况下提高精度。将该方法应用于均匀化格罗宁根区域的二维截面上,结果表明,模拟时间最多可减少7倍。
{"title":"Quadrangular adaptive mesh for elastic wave simulation in smooth anisotropic media","authors":"Marius Rapenne ,&nbsp;Paul Cupillard ,&nbsp;Guillaume Caumon ,&nbsp;Corentin Gouache","doi":"10.1016/j.cageo.2025.105946","DOIUrl":"10.1016/j.cageo.2025.105946","url":null,"abstract":"<div><div>Smooth anisotropic media are often met when implementing effective medium theory, full waveform inversion or seismic imaging. However, computational overburden is often a recurring problem when working with high frequencies or when quantifying uncertainties. In this context, adaptive meshes constitute, in principle, an attractive representation to maximize simulation accuracy while minimizing the computational cost. However, such meshes are difficult to create in the context of smooth anisotropic media as the optimal local size of the elements is not clearly defined. In this work, we present a two-step algorithm to efficiently mesh these media for spectral element method (SEM) simulation in the 2D elastic case. Our algorithm yields quadrangular only meshes which adapt the size of the element to the local and directional S-wave velocity. It relies on a quadtree division introduced by Maréchal (2009) to divide the mesh until the size of each element edge is adapted to the local minimum wavelength that will be propagated. Then, a Laplacian smoothing is applied to further optimize the size of the elements, increasing the global time step and makes the SEM simulation faster while keeping a good accuracy and even improving it in some cases. An application of our method on a 2D section of the homogenized Groningen area shows that simulation time can be reduced by a factor up to 7.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"203 ","pages":"Article 105946"},"PeriodicalIF":4.2,"publicationDate":"2025-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144148034","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
期刊
Computers & Geosciences
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1