Computational Geosciences最新文献_第6页

A generalized time-domain velocity-stress seismic wave equation for composite viscoelastic media with a topographic relief and an irregular seabed 具有地形起伏和不规则海床的复合粘弹性介质的广义时域速度应力地震波方程

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2024-02-01 DOI: 10.1007/s10596-024-10273-2

Chao Jin, Bing Zhou, Mohamed Kamel Riahi, Mohamed Jamal Zemerly

{"title":"A generalized time-domain velocity-stress seismic wave equation for composite viscoelastic media with a topographic relief and an irregular seabed","authors":"Chao Jin, Bing Zhou, Mohamed Kamel Riahi, Mohamed Jamal Zemerly","doi":"10.1007/s10596-024-10273-2","DOIUrl":"https://doi.org/10.1007/s10596-024-10273-2","url":null,"abstract":"Accurate seismic wave modeling of viscoelastic anisotropic medium is a fundamental tool for seismic data processing, interpretation and full waveform inversion. Also, free water surface, topographic relief and irregular seabed are often encountered in practical seismic surveys. Thus, basing on the General Maxwell Body, we proposed a generalized matrix form of the velocity-stress seismic wave equation, which becomes valid for composite viscoelastic anisotropic media and satisfies the boundary conditions in presence of topographic free surfaces and irregular fluid–solid interfaces. We theoretically show that the viscoelastic effect of a medium may be considered as the intrinsic body sources accumulated in wavefield history and computed by a recursive convolution formula accurately and efficiently. We also demonstrated that such a generalized viscoelastic wave equation may be solved with the curvilinear MacCormack finite difference method and validated the accuracy and feasibility of the proposed method. The modeling results in homogeneous and heterogeneous media match well with the analytical solutions and the references yielded by the spectral element solutions.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"27 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139665337","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A benchmark study on reactive two-phase flow in porous media: Part I - model description 多孔介质中反应两相流的基准研究：第一部分--模型描述

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2024-01-29 DOI: 10.1007/s10596-024-10268-z

引用次数: 0

Ensemble-based history matching of the Edvard Grieg field using 4D seismic data 利用四维地震数据对埃德瓦德-格里格油气田进行基于集合的历史匹配

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2024-01-27 DOI: 10.1007/s10596-024-10275-0

Rolf J. Lorentzen, Tuhin Bhakta, Kristian Fossum, Jon André Haugen, Espen Oen Lie, Abel Onana Ndingwan, Knut Richard Straith

{"title":"Ensemble-based history matching of the Edvard Grieg field using 4D seismic data","authors":"Rolf J. Lorentzen, Tuhin Bhakta, Kristian Fossum, Jon André Haugen, Espen Oen Lie, Abel Onana Ndingwan, Knut Richard Straith","doi":"10.1007/s10596-024-10275-0","DOIUrl":"https://doi.org/10.1007/s10596-024-10275-0","url":null,"abstract":"The Edvard Grieg field is a highly complex and heterogeneous reservoir with an extensive fault structure and a mixture of sandstone, conglomerate, and shale. In this paper, we present a complete workflow for history matching the Edvard Grieg field using an ensemble smoother for Bayesian inference. An important aspect of the workflow is a methodology to check that the prior assumptions are suitable for assimilating the data, and procedures to verify that the posterior results are plausible and credible. We thoroughly describe several tools and visualization techniques for these purposes. Using these methods we show how to identify important parameters of the model. Furthermore, we utilize new compression methods for better handling large datasets. Simulating fluid flow and seismic response for reservoirs of this size and complexity requires high numerical resolution and accurate seismic models. We present a novel dual-model concept for a better representation of seismic data and attributes, that deploy different models for the underground depending on simulated properties. Results from history matching show that we can improve data matches for both production data and different seismic attributes. Updated parameters give new insight into the reservoir dynamics, and are calibrated to better represent water movement and pressure.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"16 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139583729","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Extraction of weak geochemical anomalies based on multiple-point statistics and local singularity analysis 基于多点统计和局部奇异性分析的弱地球化学异常提取方法

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2024-01-27 DOI: 10.1007/s10596-024-10272-3

Wenyao Fan, Gang Liu, Qiyu Chen, Laijun Lu, Zhesi Cui, Boxin Zuo, Xuechao Wu

{"title":"Extraction of weak geochemical anomalies based on multiple-point statistics and local singularity analysis","authors":"Wenyao Fan, Gang Liu, Qiyu Chen, Laijun Lu, Zhesi Cui, Boxin Zuo, Xuechao Wu","doi":"10.1007/s10596-024-10272-3","DOIUrl":"https://doi.org/10.1007/s10596-024-10272-3","url":null,"abstract":"Traditional interpolations might cause smoothing effect on geochemical anomaly detection due to the moving weighted average properties. Since Multiple-Point Statistics (MPS) is a kind of stochastic simulation based on regional variables statistical patterns in a certain space, it can reduce the smoothing effect and quantify the element distribution uncertainties effectively. However, due to the insufficient Training Images (TIs) in geochemical exploration fields, simulation processes cannot be directly applied on the original data. Meanwhile, element spatial distribution patterns cannot be finely characterized under single scale, with uncertainty exists during the attribute information prediction in some regions. In addition, due to the stochastic properties, it is difficult to identify geochemical anomalous information accurately based on various simulation results. Therefore, a hybrid framework combined MPS and Local Singularity Analysis (LSA) are mainly introduced in this paper. Firstly, rasterization algorithms are used to construct geochemical TI to ensure the MPS simulation processes. Then, two-step simulation, including large-scale and small-scale simulation, is applied to finely represent the geochemical element distribution patterns. Based on various simulation results, LSA and information fusion are finally introduced to construct the probability map of geochemical anomalies. The stream sediment geochemical data was mainly used in this paper to verify the feasibility of proposed methods. Results show that comparing with the Kriging-based ones, smoothing effect of different geochemical anomalous fields is significantly reduced, which shows a closer spatial correlation with the known deposits according to the ROC curve analysis. Based on the anomaly identification results, some mineralization indices can be preliminarily determined to offer some theoretical supports for further mineral exploration.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"154 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139583480","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Numerical simulation on staggered grids of three-dimensional brinkman-forchheimer flow and heat transfer in porous media 在交错网格上对多孔介质中的三维布林克曼-福克海默流动和传热进行数值模拟

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2024-01-05 DOI: 10.1007/s10596-023-10266-7

Wei Liu, Yingxue Song, Yanping Chen, Gexian Fan, Pengshan Wang, Kai Li

{"title":"Numerical simulation on staggered grids of three-dimensional brinkman-forchheimer flow and heat transfer in porous media","authors":"Wei Liu, Yingxue Song, Yanping Chen, Gexian Fan, Pengshan Wang, Kai Li","doi":"10.1007/s10596-023-10266-7","DOIUrl":"https://doi.org/10.1007/s10596-023-10266-7","url":null,"abstract":"In this paper, three-dimensional numerical algorithm is constructed to simulate the behavior of the Brinkman-Forchheimer flow and thermal fields. Numerical results of velocity, pressure and temperature are obtained by applying the efficient modified two-grid marker and cell (MAC) algorithm on staggered grids with the second-order backward difference formula (BDF2) time approximation. The modified-upwind idea is introduced to convective heat transfer equations for improving accuracy without any numerical oscillation. The second-order convergence rate can be achieved for pressure, velocity and temperature of considered three-dimensional model. Some numerical experiments are presented to illustrate the efficiency of algorithm. The numerical example with analytical solution is used to validate the effectiveness and accuracy of the algorithm by comparing with the results of traditional MAC algorithm. A time-dependent test is proposed to show a detailed sensitivity analysis to indicate the influence of parameters including the (varepsilon ), Forchheimer number, Brinkman number and thermal diffusivity on the physical properties of Brinkman-Forchheimer flow and heat transfer in porous media.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"4 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2024-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139376192","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

XFVM modelling of fracture aperture induced by shear and tensile opening. XFVM 模拟由剪切和拉伸开口引起的断裂孔径。

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2024-01-01 Epub Date: 2023-07-25 DOI: 10.1007/s10596-023-10214-5

Giulia Conti, Stephan Matthäi, Patrick Jenny

{"title":"XFVM modelling of fracture aperture induced by shear and tensile opening.","authors":"Giulia Conti, Stephan Matthäi, Patrick Jenny","doi":"10.1007/s10596-023-10214-5","DOIUrl":"10.1007/s10596-023-10214-5","url":null,"abstract":"In reservoir simulation, it is important to understand the mechanical behaviour of fractured rocks and the effect of shear and tensile displacements of fractures on their aperture. Tensile opening directly enhances the fracture aperture, whereas shear of a preexisting rough-walled fracture creates aperture changes dependent on the local stress state. Since fracture dilatation increases reservoir permeability, both processes must be included in a realistic and consistent manner into the mechanical reservoir simulation model. Here, we use the extended finite volume method (XFVM) to conduct flow and geomechanics simulations. In XFVM, fractures are embedded in a poroelastic matrix and are modelled with discontinuous basis functions. On each fracture segment the tractions and compressive forces are calculated, and one extra degree of freedom is added for both the shear and tensile displacement. In this particular XFVM implementation we assume that linear elasticity and steady state fluid pressure adequately constrain the effective stress. In this paper, shear dilation is not calculated a posteriori, but it enters the equations such that aperture changes directly affect the stress state. This is accomplished by adding shear dilation to the displacement gradients and therefore ascertains a consistent representation in the stress-strain relations and force balances. We illustrate and discuss the influence of this extra term in two simple test cases and in a realistic layer-restricted two-dimensional fracture network subjected to plausible in situ stress and pore pressure conditions.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"1 1","pages":"227-239"},"PeriodicalIF":2.5,"publicationDate":"2024-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11080645/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"52164191","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Power law slip boundary condition for Navier-Stokes equations: Discontinuous Galerkin schemes 纳维-斯托克斯方程的幂律滑移边界条件非连续 Galerkin 方案

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2023-12-29 DOI: 10.1007/s10596-023-10265-8

J. K. Djoko, V. S. Konlack, T. Sayah

引用次数: 0

On spatially correlated observations in importance sampling methods for subsidence estimation 沉降估算重要采样方法中的空间相关观测

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2023-12-06 DOI: 10.1007/s10596-023-10264-9

Samantha S. R. Kim, Femke C. Vossepoel

{"title":"On spatially correlated observations in importance sampling methods for subsidence estimation","authors":"Samantha S. R. Kim, Femke C. Vossepoel","doi":"10.1007/s10596-023-10264-9","DOIUrl":"https://doi.org/10.1007/s10596-023-10264-9","url":null,"abstract":"The particle filter is a data assimilation method based on importance sampling for state and parameter estimation. We apply a particle filter in two different quasi-static experiments with models of subsidence caused by a compacting reservoir. The first model considers uncorrelated model state variables and observations, with observed subsidence resulting from a single source of strain. In the second model, subsidence is a summation of subsidence contributions from multiple sources which causes spatial dependencies and correlations in the observed subsidence field. Assimilating these correlated subsidence fields may trigger weight collapse. With synthetic tests, we show in a model of subsidence with 50 independent state variables and spatially correlated subsidence a minimum of (varvec{10^{13}}) particles are required to have information in the posterior distribution identical to that in a model with 50 independent and spatially uncorrelated observations. Spatial correlations cause an information loss which can be quantified with mutual information. We illustrate how a stronger spatial correlation results in lower information content in the posterior and we empirically derive the required ensemble size for the importance sampling to remain effective. We furthermore illustrate how this loss of information is reflected in the log likelihood, and how this depends on the number of model state variables. Based on these empirical results, we propose criteria to evaluate the required ensemble size in data assimilation of spatially correlated observation fields.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"38 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138520456","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Parameterizing the fluid forces on limpet shells in unidirectional flow 单向流动中帽贝壳上流体力的参数化

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2023-11-25 DOI: 10.1007/s10596-023-10263-w

Carley Walker, Julian Simeonov, Ian Adams

引用次数: 0

Comparative calibration of 1D+2D and 3D hydrogeological watershed models 1D+2D与3D水文地质流域模型对比定标

IF 2.5 3区地球科学

Computational Geosciences Pub Date : 2023-11-23 DOI: 10.1007/s10596-023-10261-y

Gillien Latour, Pierre Horgue, François Renard, Romain Guibert, Gérald Debenest

{"title":"Comparative calibration of 1D+2D and 3D hydrogeological watershed models","authors":"Gillien Latour, Pierre Horgue, François Renard, Romain Guibert, Gérald Debenest","doi":"10.1007/s10596-023-10261-y","DOIUrl":"https://doi.org/10.1007/s10596-023-10261-y","url":null,"abstract":"In this work, we study the calibration of the parameters of a hydrogeological watershed model by comparing a 1D+2D approach that combines unsaturated 1D columns and a saturated 2D model, with a full 3D approach. In a first step, a heterogeneous permeability field is estimated by an inversion procedure for each model (2D saturated and 3D unsaturated). The fields obtained are similar but the calculation time is obviously much higher in the case of the 3D model: the 2D model seems therefore sufficient and more efficient to evaluate permeability fields using piezometric measurements in the case of vertically homogeneous aquifers. The second step focuses on the calibration of the hydraulic parameters by adjusting the hydraulic heights either derived from a 1D+2D reference model at several fictitious points distributed over the entire domain, or measured in a dozen real piezometers. Both approaches provide a good fit to the piezometric measurements, but the parameter values differ significantly: the van Genuchten alpha coefficient is unrealistic in the 1D+2D approach, reflecting a poorer consideration of the modeling unsaturated zone, while the porosity value is higher in the 3D approach, which can probably be remedied by developing a more suitable cost function.","PeriodicalId":10662,"journal":{"name":"Computational Geosciences","volume":"11 1","pages":""},"PeriodicalIF":2.5,"publicationDate":"2023-11-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138542803","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0