首页 > 最新文献

Computers & Geosciences最新文献

英文 中文
PoreViT: Automated pore typing in carbonate rocks using vision transformers and neighborhood features PoreViT:利用视觉变压器和邻域特征对碳酸盐岩进行自动孔隙分型
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-30 DOI: 10.1016/j.cageo.2025.106071
Yemna Qaiser , Mohammed Ishaq , Mohammed Yaqoob , Mohammed Yusuf Ansari , Isaac Sujay , Talha Khan , Harris Rabbani , Juan Carlos Laya , P.J. Moore , Thomas Daniel Seers
The classification of pores into their intrinsic depo-diagenetic or petrophysical morphotypes is a fundamental practice within carbonate petrography, providing linkage between pore-scale textures and their associated petrophysical signatures and/or paragenetic histories. Typically, pore classification is performed manually in a qualitative/semi-quantitative manner, which is hampered by inefficiency, subjectivity, and a lack of scalability. Though aimed at addressing the limitations of manual pore classification, efforts to automate petrographic pore-typing through artificial intelligence and computer vision techniques are limited by the inability of models to classify pores into genetic classes solely based upon simplistic size and shape features, which have been the focus of the existing literature. To address this nuanced classification problem, we present PoreViT: a Vision Transformer (ViT) model used to classify macropores observed in thin-sections into their respective Lucia classes (interparticle, touching vug, separate vug). The core novelty of PoreViT lies in its Feature Fusion block, which integrates ViT features, enhanced by a Global Token Addition layer, with spatial features extracted from a Convolutional Neural Network (CNN). Critically, our classifier leverages neighborhood information to provide the model with localized pore system topology, recognizing that pore types need to be identified not just by shape but also by their local spatial context. Trained and tested using 4115 labels obtained from 25 high-resolution thin-section scans, PoreViT provides an accurate, automated classification of carbonate macropores, achieving precision and recall values of 0.92 and 0.93 (macro-F1 0.92) corresponding to absolute improvements of +4.0% and +4.0%, and relative gains of +4.54% and +4.5%, respectively, over the best-performing CNN model (DenseNet121). The high throughput pore-textural classification capabilities demonstrated herein offer unprecedented opportunities in the integrated quantitative characterization of carbonates.
将孔隙划分为其内在的沉积成岩或岩石物理形态是碳酸盐岩岩石学的基本实践,它提供了孔隙尺度结构与其相关的岩石物理特征和/或共生历史之间的联系。通常,孔隙分类以定性/半定量的方式手动进行,效率低、主观性强、缺乏可扩展性。虽然旨在解决人工孔隙分类的局限性,但通过人工智能和计算机视觉技术自动化岩石学孔隙分型的努力受到模型无法仅根据简单的大小和形状特征将孔隙划分为遗传类的限制,这已经成为现有文献的重点。为了解决这个细致入微的分类问题,我们提出了PoreViT: Vision Transformer (ViT)模型,用于将薄切片中观察到的大孔隙分类为各自的Lucia类(颗粒间、接触空隙、分离空隙)。PoreViT的核心新颖之处在于其特征融合块,该块集成了由全局令牌添加层增强的ViT特征,以及从卷积神经网络(CNN)中提取的空间特征。关键的是,我们的分类器利用邻域信息为模型提供局部孔隙系统拓扑,认识到孔隙类型不仅需要通过形状来识别,还需要通过其局部空间环境来识别。使用从25个高分辨率薄层扫描中获得的4115个标签进行训练和测试,PoreViT提供了准确的、自动的碳酸盐大孔隙分类,与表现最好的CNN模型(DenseNet121)相比,精度和召回率分别达到0.92和0.93 (macro-F1 0.92),分别提高了+4.0%和+4.0%,相对增益分别为+4.54%和+4.5%。本文所展示的高通量孔隙结构分类能力为碳酸盐的综合定量表征提供了前所未有的机会。
{"title":"PoreViT: Automated pore typing in carbonate rocks using vision transformers and neighborhood features","authors":"Yemna Qaiser ,&nbsp;Mohammed Ishaq ,&nbsp;Mohammed Yaqoob ,&nbsp;Mohammed Yusuf Ansari ,&nbsp;Isaac Sujay ,&nbsp;Talha Khan ,&nbsp;Harris Rabbani ,&nbsp;Juan Carlos Laya ,&nbsp;P.J. Moore ,&nbsp;Thomas Daniel Seers","doi":"10.1016/j.cageo.2025.106071","DOIUrl":"10.1016/j.cageo.2025.106071","url":null,"abstract":"<div><div>The classification of pores into their intrinsic depo-diagenetic or petrophysical morphotypes is a fundamental practice within carbonate petrography, providing linkage between pore-scale textures and their associated petrophysical signatures and/or paragenetic histories. Typically, pore classification is performed manually in a qualitative/semi-quantitative manner, which is hampered by inefficiency, subjectivity, and a lack of scalability. Though aimed at addressing the limitations of manual pore classification, efforts to automate petrographic pore-typing through artificial intelligence and computer vision techniques are limited by the inability of models to classify pores into genetic classes solely based upon simplistic size and shape features, which have been the focus of the existing literature. To address this nuanced classification problem, we present PoreViT: a Vision Transformer (ViT) model used to classify macropores observed in thin-sections into their respective Lucia classes (interparticle, touching vug, separate vug). The core novelty of PoreViT lies in its Feature Fusion block, which integrates ViT features, enhanced by a Global Token Addition layer, with spatial features extracted from a Convolutional Neural Network (CNN). Critically, our classifier leverages neighborhood information to provide the model with localized pore system topology, recognizing that pore types need to be identified not just by shape but also by their local spatial context. Trained and tested using 4115 labels obtained from 25 high-resolution thin-section scans, PoreViT provides an accurate, automated classification of carbonate macropores, achieving precision and recall values of 0.92 and 0.93 (macro-F1 0.92) corresponding to absolute improvements of +4.0% and +4.0%, and relative gains of +4.54% and +4.5%, respectively, over the best-performing CNN model (DenseNet121). The high throughput pore-textural classification capabilities demonstrated herein offer unprecedented opportunities in the integrated quantitative characterization of carbonates.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106071"},"PeriodicalIF":4.4,"publicationDate":"2025-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145467036","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
Dual watermarking algorithm for trajectory data based on vector decomposition 基于矢量分解的轨迹数据双水印算法
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-27 DOI: 10.1016/j.cageo.2025.106072
Heyan Wang , Luanyun Hu , Yuchen Hu , Changqing Zhu , Wang Zhang , Na Ren
Trajectory data security demands simultaneous copyright protection and spatiotemporal integrity authentication, yet existing watermarking methods struggle with irreconcilable trade-offs between robustness, reversibility, and temporally-aware tamper detection. Traditional approaches often introduce irreversible geometric distortions under affine transformations and lack mechanisms to validate timestamp authenticity, leaving critical vulnerabilities in precision-sensitive applications. This paper proposes a dual watermarking framework integrating geometric-invariant reversible embeddings and semi-fragile temporal authentication. The method decomposes trajectory coordinates into affine-invariant α/β coefficients using vector decomposition, enabling high-capacity copyright watermark embedding resilient to geometric attacks. Temporal integrity verification is achieved through Pearson-correlation analysis of timestamp sequences, detecting malicious deletions with complete accuracy while tolerating legitimate temporal shifts. A bidirectional quantization modulation scheme guarantees lossless coordinate recovery, reducing restoration errors to sub-microradian precision by synchronizing geometric invariants with relative distance preservation. Comprehensive evaluations across 68,632 real-world trajectories demonstrate superior performance: the framework achieves normalized correlations above 0.94 against 45° rotation/200 % scaling attacks, 21 % higher robustness than state-of-the-art methods under compression, and full-error detection of timestamp tampering. By unifying geometric invariance, temporal causality verification, and reversible modulation within a single architecture, this work establishes a new paradigm for dual-function security in spatiotemporal data management, with direct applicability to GIS.
轨迹数据安全需要同时保护版权和时空完整性认证,然而现有的水印方法在鲁棒性、可逆性和时间感知篡改检测之间难以调和。传统的方法通常会在仿射变换下引入不可逆的几何扭曲,并且缺乏验证时间戳真实性的机制,这在精度敏感的应用中留下了严重的漏洞。提出了一种结合几何不变可逆嵌入和半脆弱时间认证的双水印框架。该方法利用矢量分解方法将轨迹坐标分解为仿射不变的α/β系数,使高容量版权水印嵌入能够抵御几何攻击。时间完整性验证是通过时间戳序列的皮尔逊相关分析来实现的,在允许合法时间偏移的同时,可以完全准确地检测恶意删除。双向量化调制方案保证了无损的坐标恢复,通过同步几何不变量和相对距离保持将恢复误差降低到亚微弧度精度。对68,632个真实世界轨迹的综合评估显示了卓越的性能:该框架在45°旋转/ 200%缩放攻击下实现了0.94以上的归一化相关性,比最先进的压缩方法高21%的鲁棒性,以及对时间戳篡改的全错误检测。通过在单一架构中统一几何不变性、时间因果验证和可逆调制,本工作建立了时空数据管理双功能安全的新范式,可直接适用于GIS。
{"title":"Dual watermarking algorithm for trajectory data based on vector decomposition","authors":"Heyan Wang ,&nbsp;Luanyun Hu ,&nbsp;Yuchen Hu ,&nbsp;Changqing Zhu ,&nbsp;Wang Zhang ,&nbsp;Na Ren","doi":"10.1016/j.cageo.2025.106072","DOIUrl":"10.1016/j.cageo.2025.106072","url":null,"abstract":"<div><div>Trajectory data security demands simultaneous copyright protection and spatiotemporal integrity authentication, yet existing watermarking methods struggle with irreconcilable trade-offs between robustness, reversibility, and temporally-aware tamper detection. Traditional approaches often introduce irreversible geometric distortions under affine transformations and lack mechanisms to validate timestamp authenticity, leaving critical vulnerabilities in precision-sensitive applications. This paper proposes a dual watermarking framework integrating geometric-invariant reversible embeddings and semi-fragile temporal authentication. The method decomposes trajectory coordinates into affine-invariant α/β coefficients using vector decomposition, enabling high-capacity copyright watermark embedding resilient to geometric attacks. Temporal integrity verification is achieved through Pearson-correlation analysis of timestamp sequences, detecting malicious deletions with complete accuracy while tolerating legitimate temporal shifts. A bidirectional quantization modulation scheme guarantees lossless coordinate recovery, reducing restoration errors to sub-microradian precision by synchronizing geometric invariants with relative distance preservation. Comprehensive evaluations across 68,632 real-world trajectories demonstrate superior performance: the framework achieves normalized correlations above 0.94 against 45° rotation/200 % scaling attacks, 21 % higher robustness than state-of-the-art methods under compression, and full-error detection of timestamp tampering. By unifying geometric invariance, temporal causality verification, and reversible modulation within a single architecture, this work establishes a new paradigm for dual-function security in spatiotemporal data management, with direct applicability to GIS.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106072"},"PeriodicalIF":4.4,"publicationDate":"2025-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145417519","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
RIPPL, a Python-based InSAR stack and tropospheric delay software package RIPPL,基于python的InSAR堆栈和对流层延迟软件包
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-22 DOI: 10.1016/j.cageo.2025.106069
G. Mulder, F.J. van Leijen, P. Lopez-Dekker, R.F. Hanssen
Interferometric Synthetic Aperture Radar (InSAR) has a wide range of applications, including the monitoring of solid-earth and cryospheric geophysical processes and the monitoring of the built environment. The use of InSAR for atmospheric applications is less thoroughly developed. To perform such analyses the atmospheric phase delay of the SAR signal between different overpasses is used, which needs to be disentangled from other phase constituents, such as displacements and topography, which requires stack processing of large data volumes. Typically, initial atmospheric delays are predicted using existing numerical weather prediction (NWP) models, but InSAR processing and NWP model delay estimation software are not well integrated. Here we present a pure Python-based software package that integrates the automatic downloading and processing of InSAR and NWP model data to create time-series of unwrapped InSAR interferograms and InSAR equivalent tropospheric delays from NWP models. By combining the geometry of the InSAR radar signals with different NWP model datasets the tropospheric delays can accurately be derived on a pixel by pixel basis.
干涉合成孔径雷达(InSAR)具有广泛的应用,包括固体地球和冰冻圈地球物理过程的监测以及建筑环境的监测。InSAR在大气应用方面的应用发展得不太彻底。为了进行这种分析,使用了不同立交桥之间SAR信号的大气相位延迟,需要将其与其他相位成分(如位移和地形)分离,这需要对大数据量进行叠加处理。通常,使用现有的数值天气预报(NWP)模式预测初始大气延迟,但InSAR处理和NWP模式延迟估计软件没有很好地集成。本文提出了一个基于python的软件包,该软件包集成了InSAR和NWP模型数据的自动下载和处理,以创建来自NWP模型的解包裹InSAR干涉图的时间序列和InSAR等效对流层延迟。通过将InSAR雷达信号的几何形状与不同的NWP模式数据集相结合,可以精确地逐像元导出对流层延迟。
{"title":"RIPPL, a Python-based InSAR stack and tropospheric delay software package","authors":"G. Mulder,&nbsp;F.J. van Leijen,&nbsp;P. Lopez-Dekker,&nbsp;R.F. Hanssen","doi":"10.1016/j.cageo.2025.106069","DOIUrl":"10.1016/j.cageo.2025.106069","url":null,"abstract":"<div><div>Interferometric Synthetic Aperture Radar (InSAR) has a wide range of applications, including the monitoring of solid-earth and cryospheric geophysical processes and the monitoring of the built environment. The use of InSAR for atmospheric applications is less thoroughly developed. To perform such analyses the atmospheric phase delay of the SAR signal between different overpasses is used, which needs to be disentangled from other phase constituents, such as displacements and topography, which requires stack processing of large data volumes. Typically, initial atmospheric delays are predicted using existing numerical weather prediction (NWP) models, but InSAR processing and NWP model delay estimation software are not well integrated. Here we present a pure Python-based software package that integrates the automatic downloading and processing of InSAR and NWP model data to create time-series of unwrapped InSAR interferograms and InSAR equivalent tropospheric delays from NWP models. By combining the geometry of the InSAR radar signals with different NWP model datasets the tropospheric delays can accurately be derived on a pixel by pixel basis.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106069"},"PeriodicalIF":4.4,"publicationDate":"2025-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145617890","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
openKARST: A novel open-source flow simulator for karst systems openKARST:一个新颖的开源岩溶系统流动模拟器
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-20 DOI: 10.1016/j.cageo.2025.106066
Jannes Kordilla, Marco Dentz, Juan J. Hidalgo
We introduce the open-source Python-based code openKARST for flow in karst conduit networks. Flow and transport in complex karst systems remain a challenging area of hydrogeological research due to the heterogeneous nature of conduit networks. Flow regimes in these systems are highly dynamic, with transitions from free-surface to fully pressurized and laminar to turbulent flow conditions and Reynolds numbers often exceeding one million. These transitions can occur simultaneously within a network, depending on conduit roughness properties and diameter distributions. openKARST solves the transient dynamic wave equation using an iterative scheme and is optimized through an efficient vectorized structure. Transitions from free-surface to pressurized flows in smooth and rough circular conduits are realized via a Preissmann slot approach in combination with an implementation of the Darcy–Weisbach and Manning equations to compute friction losses. To mitigate numerical fluctuations commonly encountered in the Colebrook–White equation, the dynamic switching from laminar to turbulent flows is modeled with a continuous Churchill formulation for the friction factor computation. openKARST supports common boundary conditions encountered in karst systems, as and includes functionalities for network import, export and visualization. The code is verified via comparison against several analytical solutions and validated against a laboratory experiment. Finally, we demonstrate the application of openKARST by simulating a synthetic recharge event in one of the largest explored karst networks, the Ox Bel Ha system in Mexico.
本文介绍了基于python的开放源代码openKARST岩溶管道网络中的流动。由于管道网络的非均质性,复杂岩溶系统的流动和输运仍然是水文地质研究的一个具有挑战性的领域。这些系统中的流动状态是高度动态的,从自由表面到全压,层流到湍流状态的转变,雷诺数通常超过一百万。根据导管的粗糙度和直径分布,这些转变可以同时发生在管网中。openKARST采用迭代格式求解瞬态动力波动方程,并通过有效的矢量化结构进行优化。通过Preissmann槽法结合Darcy-Weisbach和Manning方程计算摩擦损失,实现了光滑和粗糙圆形管道中从自由表面到加压流动的过渡。为了减轻Colebrook-White方程中经常遇到的数值波动,用连续Churchill公式模拟了从层流到湍流的动态转换,以计算摩擦系数。openKARST支持喀斯特系统中常见的边界条件,包括网络导入、导出和可视化功能。该代码通过与几个分析解的比较和与实验室实验的验证来验证。最后,我们通过模拟墨西哥Ox Bel Ha系统中最大的岩溶网络之一的综合补给事件来演示openKARST的应用。
{"title":"openKARST: A novel open-source flow simulator for karst systems","authors":"Jannes Kordilla,&nbsp;Marco Dentz,&nbsp;Juan J. Hidalgo","doi":"10.1016/j.cageo.2025.106066","DOIUrl":"10.1016/j.cageo.2025.106066","url":null,"abstract":"<div><div>We introduce the open-source Python-based code <span>openKARST</span> for flow in karst conduit networks. Flow and transport in complex karst systems remain a challenging area of hydrogeological research due to the heterogeneous nature of conduit networks. Flow regimes in these systems are highly dynamic, with transitions from free-surface to fully pressurized and laminar to turbulent flow conditions and Reynolds numbers often exceeding one million. These transitions can occur simultaneously within a network, depending on conduit roughness properties and diameter distributions. <span>openKARST</span> solves the transient dynamic wave equation using an iterative scheme and is optimized through an efficient vectorized structure. Transitions from free-surface to pressurized flows in smooth and rough circular conduits are realized via a Preissmann slot approach in combination with an implementation of the Darcy–Weisbach and Manning equations to compute friction losses. To mitigate numerical fluctuations commonly encountered in the Colebrook–White equation, the dynamic switching from laminar to turbulent flows is modeled with a continuous Churchill formulation for the friction factor computation. <span>openKARST</span> supports common boundary conditions encountered in karst systems, as and includes functionalities for network import, export and visualization. The code is verified via comparison against several analytical solutions and validated against a laboratory experiment. Finally, we demonstrate the application of <span>openKARST</span> by simulating a synthetic recharge event in one of the largest explored karst networks, the Ox Bel Ha system in Mexico.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106066"},"PeriodicalIF":4.4,"publicationDate":"2025-10-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145364038","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
Enhanced global soil moisture prediction through a sampling-weighted sensitive learning strategy applied to various LSTM-based models 基于lstm模型的采样加权敏感学习策略增强了全球土壤湿度预测
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-17 DOI: 10.1016/j.cageo.2025.106068
Xiaoning Li , Zhichao Zhong , Qingliang Li , Cheng Zhang , Hongwei Zhao , Xiaofeng Li , Jinlong Zhu , Sen Yan
Soil moisture (SM) plays a critical role in land-atmosphere interactions, influencing both water and carbon cycles. Accurate SM predictions are essential for effective disaster response, optimized irrigation practices, and progress in environmental research. Deep learning (DL) models have become increasingly popular for predicting SM. However, many existing approaches overlook the imbalance in observed data—where moderate moisture levels are far more common than extreme dry or wet conditions. This skewed distribution limits the models' ability to accurately capture rare but critical extremes, ultimately reducing their overall effectiveness. To overcome this limitation, we propose a Sampling-Weighted Sensitive Learning Strategy that improves model generalization by assigning greater importance to rare samples during training. We evaluated this approach using three widely used DL architectures: Long Short-Term Memory Network (LSTM), Bidirectional Long Short-Term Memory Network (BiLSTM), and Gated Recurrent Unit (GRU). To ensure consistency across experiments, the same random seed was applied throughout. Our results demonstrate notable improvements in prediction accuracy when applying the proposed strategy. The BiLSTM model, in particular, showed the most significant gains: unbiased Root Mean Square Error (ubRMSE) decreased by 7.38 %, and Bias was reduced by 11.64 %. Its Kling-Gupta Efficiency (KGE) improved by 2.73 %—slightly below the 5.35 % gain observed for the unidirectional LSTM—but regional results were particularly strong. In data-scarce areas, especially North Africa and Western Asia, BiLSTM KGE improvements frequently exceeded 20 %. Models trained with the proposed strategy also produced narrower 95 % confidence intervals during high-variability periods (e.g., summer and dry seasons), indicating greater predictive robustness under challenging environmental. These findings underscore the importance of addressing sample imbalance in training data and demonstrate the effectiveness of our strategy in enhancing DL models for SM prediction.
土壤水分在陆地-大气相互作用中起着至关重要的作用,影响着水和碳的循环。准确的SM预测对于有效的灾害响应、优化灌溉实践和环境研究进展至关重要。深度学习(DL)模型在预测SM方面越来越受欢迎。然而,许多现有的方法忽略了观测数据的不平衡——适度的湿度水平远比极端干燥或潮湿的条件更常见。这种倾斜的分布限制了模型准确捕捉罕见但关键的极端情况的能力,最终降低了模型的总体有效性。为了克服这一限制,我们提出了一种抽样加权敏感学习策略,该策略通过在训练过程中赋予稀有样本更大的重要性来提高模型泛化。我们使用三种广泛使用的深度学习架构来评估这种方法:长短期记忆网络(LSTM)、双向长短期记忆网络(BiLSTM)和门控循环单元(GRU)。为了确保实验的一致性,在整个过程中使用相同的随机种子。我们的结果表明,当应用所提出的策略时,预测精度显着提高。特别是BiLSTM模型,显示出最显著的增益:无偏均方根误差(ubRMSE)降低了7.38%,偏差降低了11.64%。它的克林-古普塔效率(KGE)提高了2.73%,略低于单向lstm观察到的5.35%的增益,但区域结果特别强劲。在数据匮乏的地区,特别是北非和西亚,BiLSTM KGE的改进经常超过20%。采用该策略训练的模型在高变异性时期(如夏季和旱季)也产生了更窄的95%置信区间,表明在具有挑战性的环境下具有更强的预测稳健性。这些发现强调了解决训练数据中样本不平衡的重要性,并证明了我们的策略在增强深度学习模型用于SM预测方面的有效性。
{"title":"Enhanced global soil moisture prediction through a sampling-weighted sensitive learning strategy applied to various LSTM-based models","authors":"Xiaoning Li ,&nbsp;Zhichao Zhong ,&nbsp;Qingliang Li ,&nbsp;Cheng Zhang ,&nbsp;Hongwei Zhao ,&nbsp;Xiaofeng Li ,&nbsp;Jinlong Zhu ,&nbsp;Sen Yan","doi":"10.1016/j.cageo.2025.106068","DOIUrl":"10.1016/j.cageo.2025.106068","url":null,"abstract":"<div><div>Soil moisture (SM) plays a critical role in land-atmosphere interactions, influencing both water and carbon cycles. Accurate SM predictions are essential for effective disaster response, optimized irrigation practices, and progress in environmental research. Deep learning (DL) models have become increasingly popular for predicting SM. However, many existing approaches overlook the imbalance in observed data—where moderate moisture levels are far more common than extreme dry or wet conditions. This skewed distribution limits the models' ability to accurately capture rare but critical extremes, ultimately reducing their overall effectiveness. To overcome this limitation, we propose a Sampling-Weighted Sensitive Learning Strategy that improves model generalization by assigning greater importance to rare samples during training. We evaluated this approach using three widely used DL architectures: Long Short-Term Memory Network (LSTM), Bidirectional Long Short-Term Memory Network (BiLSTM), and Gated Recurrent Unit (GRU). To ensure consistency across experiments, the same random seed was applied throughout. Our results demonstrate notable improvements in prediction accuracy when applying the proposed strategy. The BiLSTM model, in particular, showed the most significant gains: unbiased Root Mean Square Error (<em>ubRMSE</em>) decreased by 7.38 %, and <em>Bias</em> was reduced by 11.64 %. Its Kling-Gupta Efficiency (<em>KGE</em>) improved by 2.73 %—slightly below the 5.35 % gain observed for the unidirectional LSTM—but regional results were particularly strong. In data-scarce areas, especially North Africa and Western Asia, BiLSTM KGE improvements frequently exceeded 20 %. Models trained with the proposed strategy also produced narrower 95 % confidence intervals during high-variability periods (e.g., summer and dry seasons), indicating greater predictive robustness under challenging environmental. These findings underscore the importance of addressing sample imbalance in training data and demonstrate the effectiveness of our strategy in enhancing DL models for SM prediction.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106068"},"PeriodicalIF":4.4,"publicationDate":"2025-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145321635","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
Invertible neural network for real-time inversion and uncertainty quantification of ultra-deep resistivity measurements 超深电阻率测量实时反演与不确定度量化的可逆神经网络
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-17 DOI: 10.1016/j.cageo.2025.106067
George Bittar , Sihong Wu , Yawei Su , Shubin Zeng , Jiajia Sun , Xuqing Wu , Yueqin Huang , Jiefu Chen
Real-time geosteering, formation evaluation, and wellbore placement decisions hinge on the ability to invert electromagnetic (EM) well logging measurements in a fast manner while understanding the associated uncertainties. Conventional deterministic inversion methods, such as the Levenberg–Marquardt algorithm (LMA) and Occam’s inversion, often get trapped in local minima and yield a single optimal solution, neglecting the impact of the non-uniqueness of solutions. Bayesian approaches like Markov Chain Monte Carlo (MCMC) can provide the posterior distribution but are computationally expensive, making them impractical for real-time inversion. In this study, we develop a deep learning-based invertible neural network (INN) that performs rapid approximate Bayesian inversion under a specific likelihood and provides uncertainty quantification (UQ) for ultra-deep resistivity measurements. Synthetic tests demonstrate that the INN recovers the posterior distribution and generates an ensemble of predictions to quantify uncertainty within seconds. We compare its performance with conventional inversion algorithms, including LMA and Occam’s inversion, evaluating accuracy and inference efficiency. The results show that the INN delivers reliable resistivity inversion with uncertainty information at a fraction of the computational cost, highlighting its potential for real-time geosteering and other drilling-related decision-making tasks.
实时地质导向、地层评价和井筒布置决策取决于能否快速反演电磁测井数据,同时了解相关的不确定性。传统的确定性反演方法,如Levenberg-Marquardt算法(LMA)和Occam反演,往往陷入局部最小值,只得到一个最优解,而忽略了解的非唯一性的影响。马尔可夫链蒙特卡罗(MCMC)等贝叶斯方法可以提供后验分布,但计算成本高,无法实现实时反演。在本研究中,我们开发了一种基于深度学习的可逆神经网络(INN),该网络在特定似然下执行快速近似贝叶斯反演,并为超深电阻率测量提供不确定性量化(UQ)。综合测试表明,INN可以在数秒内恢复后验分布并生成预测集合以量化不确定性。我们将其性能与传统的反演算法,包括LMA和Occam的反演进行了比较,评估了精度和推理效率。结果表明,INN能够以很小的计算成本提供可靠的不确定信息电阻率反演,突出了其在实时地质导向和其他钻井相关决策任务方面的潜力。
{"title":"Invertible neural network for real-time inversion and uncertainty quantification of ultra-deep resistivity measurements","authors":"George Bittar ,&nbsp;Sihong Wu ,&nbsp;Yawei Su ,&nbsp;Shubin Zeng ,&nbsp;Jiajia Sun ,&nbsp;Xuqing Wu ,&nbsp;Yueqin Huang ,&nbsp;Jiefu Chen","doi":"10.1016/j.cageo.2025.106067","DOIUrl":"10.1016/j.cageo.2025.106067","url":null,"abstract":"<div><div>Real-time geosteering, formation evaluation, and wellbore placement decisions hinge on the ability to invert electromagnetic (EM) well logging measurements in a fast manner while understanding the associated uncertainties. Conventional deterministic inversion methods, such as the Levenberg–Marquardt algorithm (LMA) and Occam’s inversion, often get trapped in local minima and yield a single optimal solution, neglecting the impact of the non-uniqueness of solutions. Bayesian approaches like Markov Chain Monte Carlo (MCMC) can provide the posterior distribution but are computationally expensive, making them impractical for real-time inversion. In this study, we develop a deep learning-based invertible neural network (INN) that performs rapid approximate Bayesian inversion under a specific likelihood and provides uncertainty quantification (UQ) for ultra-deep resistivity measurements. Synthetic tests demonstrate that the INN recovers the posterior distribution and generates an ensemble of predictions to quantify uncertainty within seconds. We compare its performance with conventional inversion algorithms, including LMA and Occam’s inversion, evaluating accuracy and inference efficiency. The results show that the INN delivers reliable resistivity inversion with uncertainty information at a fraction of the computational cost, highlighting its potential for real-time geosteering and other drilling-related decision-making tasks.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106067"},"PeriodicalIF":4.4,"publicationDate":"2025-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145364037","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 new elastic wave equation for decoupling P-wave and S-waves and its application 纵波与横波解耦的弹性波动方程及其应用
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-14 DOI: 10.1016/j.cageo.2025.106065
Meng Guo , Bingshou He , Qianqian Ci
The imaging of P-wave and S-wave in reverse time migration (RTM) of elastic waves is often achieved by cross-correlating P-waves or S-waves with different propagation directions. This requires us to obtain the Poynting vector or optical flow vector of each imaging point at different times during the wavefield extrapolation process and use it to indicate the direction of wave propagation. But we can only obtain the Poynting vector of the mixed wavefield of P-wave and S- waves, and we cannot obtain the Poynting vector of pure P-wave or pure S-wave when using the existing velocity-stress elastic wave equations for the wavefield extrapolation process. Therefore, the propagation direction obtained is also a mixed wavefield rather than pure P-wave or pure S-wave, and this does not meet the requirements for elastic wave RTM and will cause errors. The existing first-order velocity-dilation-rotation elastic wave equation, although it overcomes the aforementioned issues, cannot accurately describe the law of wave propagation at the wave impedance interface due to the assumption of a homogeneous medium. Especially when the interface of P-wave and S-wave velocities is not consistent, it will lead to errors in the reflection, transmission, and conversion wavefields when using this equation for elastic wavefield extrapolation. In addition, severe energy leakage effects will occur at the interface of S-wave velocity when using this equation, which will lead to inaccurate S-wave imaging. In this paper, we propose a new elastic wave equation for decoupling P-wave and S-waves based on the assumption of an inhomogeneous medium, which not only gives the propagation direction of pure P-wave and pure S-wave, but also completely overcomes the above problems. Using the new equation of the Poynting vector in the elastic wave field to perform cross-correlation imaging, the model calculations show that the imaging results eliminate the noise generated by RTM, demonstrating the accuracy and applicability of the equation.
弹性波逆时偏移(RTM)中的纵波和横波成像通常是通过不同传播方向的纵波或横波相互关联来实现的。这就要求我们在波场外推过程中,获取每个成像点在不同时刻的坡印亭矢量或光流矢量,并用它来指示波的传播方向。但我们只能得到纵波和横波混合波场的Poynting矢量,而用现有的速度-应力弹性波方程进行波场外推时,无法得到纯纵波或纯横波的Poynting矢量。因此,得到的传播方向也是混合波场,而不是纯p波或纯s波,这不符合弹性波RTM的要求,会产生误差。现有的一阶速度-膨胀-旋转弹性波动方程虽然克服了上述问题,但由于假设介质均质,无法准确描述波在波阻抗界面处的传播规律。特别是当纵波和横波速度界面不一致时,用该方程进行弹性波场外推时,会导致反射、透射和转换波场出现误差。此外,使用该方程时,在横波速度界面处会产生严重的能量泄漏效应,导致横波成像不准确。本文基于非均匀介质的假设,提出了一种新的纵波与横波解耦的弹性波动方程,不仅给出了纯纵波和纯横波的传播方向,而且完全克服了上述问题。利用弹性波场中新的Poynting矢量方程进行互相关成像,模型计算表明,成像结果消除了RTM产生的噪声,证明了该方程的准确性和适用性。
{"title":"A new elastic wave equation for decoupling P-wave and S-waves and its application","authors":"Meng Guo ,&nbsp;Bingshou He ,&nbsp;Qianqian Ci","doi":"10.1016/j.cageo.2025.106065","DOIUrl":"10.1016/j.cageo.2025.106065","url":null,"abstract":"<div><div>The imaging of P-wave and S-wave in reverse time migration (RTM) of elastic waves is often achieved by cross-correlating P-waves or S-waves with different propagation directions. This requires us to obtain the Poynting vector or optical flow vector of each imaging point at different times during the wavefield extrapolation process and use it to indicate the direction of wave propagation. But we can only obtain the Poynting vector of the mixed wavefield of P-wave and S- waves, and we cannot obtain the Poynting vector of pure P-wave or pure S-wave when using the existing velocity-stress elastic wave equations for the wavefield extrapolation process. Therefore, the propagation direction obtained is also a mixed wavefield rather than pure P-wave or pure S-wave, and this does not meet the requirements for elastic wave RTM and will cause errors. The existing first-order velocity-dilation-rotation elastic wave equation, although it overcomes the aforementioned issues, cannot accurately describe the law of wave propagation at the wave impedance interface due to the assumption of a homogeneous medium. Especially when the interface of P-wave and S-wave velocities is not consistent, it will lead to errors in the reflection, transmission, and conversion wavefields when using this equation for elastic wavefield extrapolation. In addition, severe energy leakage effects will occur at the interface of S-wave velocity when using this equation, which will lead to inaccurate S-wave imaging. In this paper, we propose a new elastic wave equation for decoupling P-wave and S-waves based on the assumption of an inhomogeneous medium, which not only gives the propagation direction of pure P-wave and pure S-wave, but also completely overcomes the above problems. Using the new equation of the Poynting vector in the elastic wave field to perform cross-correlation imaging, the model calculations show that the imaging results eliminate the noise generated by RTM, demonstrating the accuracy and applicability of the equation.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106065"},"PeriodicalIF":4.4,"publicationDate":"2025-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145364036","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
Parallel finite element forward modeling of 3-D magnetotelluric conductivity and permeability anisotropy with coupled PML boundary conditions 耦合PML边界条件下三维大地电磁导电性和渗透率各向异性的平行有限元正演模拟
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-10 DOI: 10.1016/j.cageo.2025.106064
Shuaiying Qiao , Tiaojie Xiao , Junjun Zhou , Chunye Gong , Bo Yang , Jie Liu , Yun Wang , Qinglin Wang
Magnetotelluric sounding (MT) is a crucial geophysical exploration method, with its response primarily influenced by two physical parameters: conductivity and magnetic permeability. MT forward modeling typically presents as a large-scale, open-domain problem, necessitating boundary truncation and computational acceleration for the simulation area. Compared to traditional boundary conditions, the Perfectly Matched Layer (PML) offers a more efficient and accurate truncation method. However, the current application of the PML is confined to scenarios involving variations in conductivity alone, and is unable to accommodate simultaneous variations in both conductivity and permeability, as well as complex anisotropic models. Therefore, this paper proposes a PML that accounts for both conductivity and permeability parameters, as well as anisotropy, making it suitable for complex anisotropic models. Furthermore, by integrating the Multi-Processing Interface (MPI) to design a multi-level parallel processing scheme, we have achieved parallel vector finite element forward modeling of three-dimensional magnetotelluric conductivity and permeability anisotropy with coupled PML boundary conditions. In comparison with previous results, the PML boundary conditions have been validated to possess the advantages of high efficiency, high precision, and stable performance. Numerical experimental results indicate that, compared with traditional boundary conditions, the PML reduces the degrees of freedom (DOFs) by over 85%, and decreases both computation time and memory usage by more than 90%. Compared with the conventional method with 888,822 DOFs, the proposed method, which integrates the PML and a multi-level parallelization strategy, achieves a speedup of approximately 85.24 for a single frequency using 32 processes and approximately 649.63 for 8 frequencies using 512 processes. The PML boasts a wider range of applicability, better performance, and thus holds broader prospects for application.
大地电磁测深是一种重要的地球物理勘探方法,其响应主要受电导率和磁导率两个物理参数的影响。MT正演建模通常是一个大规模的开放域问题,需要对模拟区域进行边界截断和计算加速。与传统的边界条件相比,完美匹配层(PML)提供了一种更有效、更准确的截断方法。然而,目前PML的应用仅限于仅涉及电导率变化的场景,无法适应电导率和渗透率的同时变化,以及复杂的各向异性模型。因此,本文提出了一种同时考虑电导率和渗透率参数以及各向异性的PML,使其适用于复杂的各向异性模型。此外,通过集成多处理接口(MPI)设计多层次并行处理方案,实现了具有耦合PML边界条件的三维大地电磁导电性和渗透率各向异性的并行矢量有限元正演模拟。通过与已有结果的比较,验证了PML边界条件具有效率高、精度高、性能稳定等优点。数值实验结果表明,与传统边界条件相比,该边界条件的自由度降低了85%以上,计算时间和内存占用均降低了90%以上。与传统方法的888,822 dof相比,该方法集成了PML和多级并行化策略,使用32个进程实现了单个频率约85.24的加速,使用512个进程实现了8个频率约649.63的加速。PML的适用范围更广,性能更好,具有更广阔的应用前景。
{"title":"Parallel finite element forward modeling of 3-D magnetotelluric conductivity and permeability anisotropy with coupled PML boundary conditions","authors":"Shuaiying Qiao ,&nbsp;Tiaojie Xiao ,&nbsp;Junjun Zhou ,&nbsp;Chunye Gong ,&nbsp;Bo Yang ,&nbsp;Jie Liu ,&nbsp;Yun Wang ,&nbsp;Qinglin Wang","doi":"10.1016/j.cageo.2025.106064","DOIUrl":"10.1016/j.cageo.2025.106064","url":null,"abstract":"<div><div>Magnetotelluric sounding (MT) is a crucial geophysical exploration method, with its response primarily influenced by two physical parameters: conductivity and magnetic permeability. MT forward modeling typically presents as a large-scale, open-domain problem, necessitating boundary truncation and computational acceleration for the simulation area. Compared to traditional boundary conditions, the Perfectly Matched Layer (PML) offers a more efficient and accurate truncation method. However, the current application of the PML is confined to scenarios involving variations in conductivity alone, and is unable to accommodate simultaneous variations in both conductivity and permeability, as well as complex anisotropic models. Therefore, this paper proposes a PML that accounts for both conductivity and permeability parameters, as well as anisotropy, making it suitable for complex anisotropic models. Furthermore, by integrating the Multi-Processing Interface (MPI) to design a multi-level parallel processing scheme, we have achieved parallel vector finite element forward modeling of three-dimensional magnetotelluric conductivity and permeability anisotropy with coupled PML boundary conditions. In comparison with previous results, the PML boundary conditions have been validated to possess the advantages of high efficiency, high precision, and stable performance. Numerical experimental results indicate that, compared with traditional boundary conditions, the PML reduces the degrees of freedom (DOFs) by over 85%, and decreases both computation time and memory usage by more than 90%. Compared with the conventional method with 888,822 DOFs, the proposed method, which integrates the PML and a multi-level parallelization strategy, achieves a speedup of approximately 85.24 for a single frequency using 32 processes and approximately 649.63 for 8 frequencies using 512 processes. The PML boasts a wider range of applicability, better performance, and thus holds broader prospects for application.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106064"},"PeriodicalIF":4.4,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145321623","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
Eigenvector decomposition for joint analysis of spatial characteristics in the North Atlantic from 1979 to 2024 1979 - 2024年北大西洋空间特征联合分析的特征向量分解
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-09 DOI: 10.1016/j.cageo.2025.106062
Andrey K. Gorshenin , Anastasiia A. Osipova , Konstantin P. Belyaev
The extension of the use of Itô stochastic differential equations (SDEs) for joint analysis of spatio-temporal characteristics in the North Atlantic region, such as sea surface temperature (SST), the sum of sensible and latent heat fluxes, and surface atmospheric pressure for the period between 1979 and 2024 is introduced. Previously, this model was used only for the fluxes. The joint point estimates for the random coefficients of SDEs as multidimensional matrices (the drift vector and the diffusion matrix) are obtained for the entire considered period. The numerical estimations of these values were carried out using high-performance computing equipment with software implementation in Python language using the reanalysis data from the ERA5 database. Developed methods and tools are used for the statistical analysis of the temporal evolution of the coefficients of the Itô equation, analysis of joint and marginal diffusion matrices, their finite-dimensional Karhunen–Loéve’s decomposition into eigenvalues and eigenvectors, determination of their interrelations, temporal trends, as well as dynamic visualization on geographical maps of the region under study. The spatial structure of the eigenvectors of the diffusion matrix, their time evolution and the relationship to jet streams and large-scale heat waves that determine latitudinal heat transfer in the North Atlantic are shown. It is also demonstrated that there is a positive trend in the interannual variability in drift and diffusion coefficients. This indicates a quantitative and qualitative increase in the air–sea interaction and the relationship between heat fluxes and ocean surface temperature. It also makes it possible to quantify the energy exchange between the ocean and atmosphere on an interannual scale. The way of using quantities from a stochastic model to improve the neural network forecasts is also discussed.
介绍了Itô随机微分方程(SDEs)在1979 - 2024年北大西洋地区海温、感热通量和潜热通量和大气压力等时空特征联合分析中的推广应用。以前,该模型仅用于通量。在整个考虑周期内,得到了SDEs随机系数作为多维矩阵(漂移向量和扩散矩阵)的结合点估计。利用ERA5数据库的再分析数据,利用高性能计算设备和Python语言软件实现对这些值的数值估计。已开发的方法和工具用于统计分析Itô方程系数的时间演变,分析联合和边际扩散矩阵,分析它们的有限维karhunen - lo分解为特征值和特征向量,确定它们的相互关系,时间趋势,以及在研究区域的地理地图上动态可视化。给出了扩散矩阵特征向量的空间结构、时间演化及其与决定北大西洋纬向传热的急流和大尺度热浪的关系。漂移系数和扩散系数的年际变化也呈正趋势。这表明海气相互作用以及热通量与海洋表面温度之间的关系在数量和质量上都有所增加。这也使得在年际尺度上量化海洋和大气之间的能量交换成为可能。本文还讨论了利用随机模型中的量来改进神经网络预测的方法。
{"title":"Eigenvector decomposition for joint analysis of spatial characteristics in the North Atlantic from 1979 to 2024","authors":"Andrey K. Gorshenin ,&nbsp;Anastasiia A. Osipova ,&nbsp;Konstantin P. Belyaev","doi":"10.1016/j.cageo.2025.106062","DOIUrl":"10.1016/j.cageo.2025.106062","url":null,"abstract":"<div><div>The extension of the use of Itô stochastic differential equations (SDEs) for joint analysis of spatio-temporal characteristics in the North Atlantic region, such as sea surface temperature (SST), the sum of sensible and latent heat fluxes, and surface atmospheric pressure for the period between 1979 and 2024 is introduced. Previously, this model was used only for the fluxes. The joint point estimates for the random coefficients of SDEs as multidimensional matrices (the drift vector and the diffusion matrix) are obtained for the entire considered period. The numerical estimations of these values were carried out using high-performance computing equipment with software implementation in Python language using the reanalysis data from the ERA5 database. Developed methods and tools are used for the statistical analysis of the temporal evolution of the coefficients of the Itô equation, analysis of joint and marginal diffusion matrices, their finite-dimensional Karhunen–Loéve’s decomposition into eigenvalues and eigenvectors, determination of their interrelations, temporal trends, as well as dynamic visualization on geographical maps of the region under study. The spatial structure of the eigenvectors of the diffusion matrix, their time evolution and the relationship to jet streams and large-scale heat waves that determine latitudinal heat transfer in the North Atlantic are shown. It is also demonstrated that there is a positive trend in the interannual variability in drift and diffusion coefficients. This indicates a quantitative and qualitative increase in the air–sea interaction and the relationship between heat fluxes and ocean surface temperature. It also makes it possible to quantify the energy exchange between the ocean and atmosphere on an interannual scale. The way of using quantities from a stochastic model to improve the neural network forecasts is also discussed.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106062"},"PeriodicalIF":4.4,"publicationDate":"2025-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145269659","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
TorchTEM3D: PyTorch-Driven forward modeling platform for fast 3D transient electromagnetic modeling and efficient sensitivity matrix calculation TorchTEM3D: pytorch驱动的正演建模平台,用于快速3D瞬变电磁建模和高效灵敏度矩阵计算
IF 4.4 2区 地球科学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-10-04 DOI: 10.1016/j.cageo.2025.106063
Ziteng Li , Hai Li , Keying Li , Ahmed M. Beshr
The three-dimensional (3D) forward modeling of transient electromagnetic (TEM) data is often computationally demanding due to its high complexity and limited hardware acceleration, which also affects the efficiency of sensitivity matrix calculation. In recent years, deep learning frameworks, particularly PyTorch, have been widely used in various fields due to their high flexibility, parallel computing capabilities, and powerful automatic differentiation function. In this paper, we develop a time-domain finite-difference forward modeling platform for 3D TEM, named TorchTEM3D, based on the powerful parallel computing and GPU acceleration capabilities of PyTorch. By fully utilizing the automatic differentiation function of PyTorch, we achieve efficient and fast calculation of sensitivity matrix (the gradient of the electromagnetic response to the geoelectric model). Compared with existing open-source Python computing platforms such as SimPEG and custEM, our method improves computing speed by 15–60 times. Furthermore, high-precision sensitivity matrices can be obtained with a single forward modeling run.
瞬变电磁(TEM)数据的三维正演建模由于其高复杂性和有限的硬件加速,往往需要大量的计算量,这也影响了灵敏度矩阵计算的效率。近年来,深度学习框架,特别是PyTorch,由于其高灵活性、并行计算能力和强大的自动微分功能,在各个领域得到了广泛的应用。本文基于PyTorch强大的并行计算和GPU加速能力,开发了三维瞬变电磁法时域有限差分正演建模平台TorchTEM3D。充分利用PyTorch的自动微分功能,实现了灵敏度矩阵(电磁响应对地电模型的梯度)的高效快速计算。与现有的开源Python计算平台SimPEG和custEM相比,我们的方法将计算速度提高了15-60倍。此外,单次正演模拟可以获得高精度的灵敏度矩阵。
{"title":"TorchTEM3D: PyTorch-Driven forward modeling platform for fast 3D transient electromagnetic modeling and efficient sensitivity matrix calculation","authors":"Ziteng Li ,&nbsp;Hai Li ,&nbsp;Keying Li ,&nbsp;Ahmed M. Beshr","doi":"10.1016/j.cageo.2025.106063","DOIUrl":"10.1016/j.cageo.2025.106063","url":null,"abstract":"<div><div>The three-dimensional (3D) forward modeling of transient electromagnetic (TEM) data is often computationally demanding due to its high complexity and limited hardware acceleration, which also affects the efficiency of sensitivity matrix calculation. In recent years, deep learning frameworks, particularly PyTorch, have been widely used in various fields due to their high flexibility, parallel computing capabilities, and powerful automatic differentiation function. In this paper, we develop a time-domain finite-difference forward modeling platform for 3D TEM, named TorchTEM3D, based on the powerful parallel computing and GPU acceleration capabilities of PyTorch. By fully utilizing the automatic differentiation function of PyTorch, we achieve efficient and fast calculation of sensitivity matrix (the gradient of the electromagnetic response to the geoelectric model). Compared with existing open-source Python computing platforms such as SimPEG and custEM, our method improves computing speed by 15–60 times. Furthermore, high-precision sensitivity matrices can be obtained with a single forward modeling run.</div></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"207 ","pages":"Article 106063"},"PeriodicalIF":4.4,"publicationDate":"2025-10-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145269658","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