Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization   总被引：2，自引：0，他引：2  

Satomi Suzuki  Guillaume Caumon  Jef Caers 《Computational Geosciences》2008,12(1):105-119

This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response. The proposed workflow consists of searching for models that match production history from a large set of prior structural model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic reservoir examples are presented to prove the applicability of the proposed method.  相似文献   

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments     

D. Pardo  C. Torres-Verdín  M. Paszynski 《Computational Geosciences》2008,12(1):83-89

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.  相似文献   

Geochemical Eccentricity of Ground Water Allied to Weathering of Basalts from the Deccan Volcanic Province,India: Insinuation on CO<Subscript>2</Subscript> Consumption     

N. J. Pawar  J. B. Pawar  Suyash Kumar  Ashwini Supekar 《Aquatic Geochemistry》2008,14(1):41-71

Analyses of 72 samples from Upper Panjhara basin in the northern part of Deccan Plateau, India, indicate that geochemical incongruity of groundwater is largely a function of mineral composition of the basaltic lithology. Higher proportion of alkaline earth elements to total cations and HCO₃>Cl + SO₄ reflect weathering of primary silicates as chief source of ions. Inputs of Cl, SO₄, and NO₃ are related to rainfall and localized anthropogenic factors. Groundwater from recharge area representing Ca + Mg–HCO₃ type progressively evolves to Ca + Na–HCO₃ and Na–Ca–HCO₃ class along flow direction replicates the role of cation exchange and precipitation processes. While the post-monsoon chemistry is controlled by silicate mineral dissolution + cation exchange reactions, pre-monsoon variability is attributable chiefly to precipitation reactions + anthropogenic factors. Positive correlations between Mg vs HCO₃ and Ca + Mg vs HCO₃ supports selective dissolution of olivine and pyroxene as dominant process in post-monsoon followed by dissolution of plagioclase feldspar and secondary carbonates. The pre-monsoon data however, points toward the dissolution of plagioclase and precipitation of CaCO₃ supported by improved correlation coefficients between Na + Ca vs HCO₃ and negative correlation of Ca vs HCO₃, respectively. It is proposed that the eccentricity in the composition of groundwater from the Panjhara basin is a function of selective dissolution of olivine > pyroxene followed by plagioclase feldspar. The data suggest siallitization (L < R and R _k) as dominant mechanism of chemical weathering of basalts, stimulating monosiallitic (kaolinite) and bisiallitic (montmorillonite) products. The chemical denudation rates for Panjhara basin worked out separately for the ground and surface water component range from 6.98 to 36.65 tons/km²/yr, respectively. The values of the CO₂ consumption rates range between 0.18 × 10⁶ mol//km²/yr (groundwater) and 0.9 × 10⁶ mol/km²/yr (surface water), which indicates that the groundwater forms a considerable fraction of CO₂ consumption, an inference, that is, not taken into contemplation in most of the studies.  相似文献   

Numerical simulation of a stratigraphic model     

Véronique Gervais  Roland Masson 《Computational Geosciences》2008,12(2):163-179

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.  相似文献   

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.  相似文献   

Environmental tracers as indicators of karst conduits in groundwater in South Dakota,USA   总被引：3，自引：1，他引：2  

Andrew J. Long  J. Foster Sawyer  Larry D. Putnam 《Hydrogeology Journal》2008,16(2):263-280

Environmental tracers sampled from the carbonate Madison aquifer on the eastern flank of the Black Hills, South Dakota, USA indicated the approximate locations of four major karst conduits. Contamination issues are a major concern because these conduits are characterized by direct connections to sinking streams, high groundwater velocities, and proximity to public water supplies. Objectives of the study were to estimate approximate conduit locations and assess possible anthropogenic influences associated with conduits. Anomalies of young groundwater based on chlorofluorocarbons (CFCs), tritium, and electrical conductivity (EC) indicated fast moving, focused flow and thus the likely presence of conduits. δ¹⁸O was useful for determining sources of recharge for each conduit, and nitrate was a useful tracer for assessing flow paths for anthropogenic influences. Two of the four conduits terminate at or near a large spring complex. CFC apparent ages ranged from 15 years near conduits to >50 years in other areas. Nitrate-N concentrations >0.4 mg/L in groundwater were associated with each of the four conduits compared with concentrations ranging from <0.1 to 0.4 mg/L in other areas. These higher nitrate-N concentrations probably do not result from sinking streams but rather from other areas of infiltration.  相似文献   

Estimating groundwater evapotranspiration rates using diurnal water-table fluctuations in a semi-arid riparian zone   总被引：2，自引：1，他引：1  

Laura K. Lautz 《Hydrogeology Journal》2008,16(3):483-497

In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ET_G) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ET_G, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ET_G to potential evapotranspiration (PET), K _c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K _c, in semi-arid rangelands. Before vegetation senescence, K _c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K _c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ET_G is primarily from estimates of specific yield (S _y), which is difficult to quantify in heterogeneous sediments. ET_G values may be more reliable because the range of acceptable S _y is smaller than K _c and S _y does not change with vegetation type or soil moisture.  相似文献   

Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)     

Francois Larroque  Wiktor Treichel  Alain Dupuy 《Hydrogeology Journal》2008,16(2):215-233

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.

  相似文献   

Contributions of boron isotopes to understanding the hydrogeochemistry of the coastal detritic aquifer of Castellón Plain,Spain   总被引：1，自引：0，他引：1  

Elena Giménez Forcada  Ignacio Morell Evangelista 《Hydrogeology Journal》2008,16(3):547-557

The Castellón Plain alluvial aquifer, Spain, is intensively exploited to meet the demand for agricultural irrigation and industrial water supply. The geochemistry of its groundwater shows complex salinization in the northern and southern parts of the aquifer, with significant pollution from human origin in the central portion. Boron content and B isotope geochemistry are useful for distinguishing between various sources of pollution and their relative importance in different parts of this aquifer. Boron concentrations in the groundwater vary between 0.01 and 0.85 mg/L. In the more saline groundwaters, found at the northern and southern ends of the study area, the presence of B is linked to inputs from seawater and water with a calcium-magnesium sulphate facies, which feed the aquifer and clearly influence the chemistry of its waters. Evidence of B adsorption processes in some samples is shown by the low B/Cl ratios and the high values of δ¹¹B. In the central portion of the aquifer, the high B/Cl ratios and the strongly negative δ¹¹B are related to pollution of human origin.  相似文献   

Water balance in the playa-lakes of an arid environment,Monegros, NE Spain   总被引：1，自引：0，他引：1  

Carmen Castañeda  Miguel Ángel García-Vera 《Hydrogeology Journal》2008,16(1):87-102

The playa-lakes of the Monegros desert in north-east Spain are saline wetlands in an arid environment, a rare phenomenon in Europe. These extremely valuable habitats are threatened by changes associated with agricultural expansion and incorporation of new irrigated areas. An understanding of the present hydrologic regime will enable changes to be identified, particularly those brought about by flooding and pollution caused by irrigation surplus. This study sets out to show the results of applying a daily water balance in three selected playa-lakes. The balance was in two parts and consisted of: (1) the average balance for all the endorheic basin using the BALAN_11 program, and (2) the water balance in some playa-lakes, applying discharge flows obtained from the previous balance. The resulting volumes of water were converted to water depths and contrasted with reference volumes taken from field and Landsat images. The model was calibrated by applying various hypotheses of function which enabled the results to be adjusted. The proposed balance is an acceptable reproduction of field water measurements during this period, and underlines the consistency of the conceptual model. The methodology used is appropriate for understanding the playa-lakes function and for monitoring them for conservation purposes.  相似文献