Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study     

F.A. Radu  N. Suciu  J. HoffmannA. Vogel  O. Kolditz  C.-H. ParkS. Attinger 《Advances in water resources》2011,34(1):47-61

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.  相似文献   

Influence of background heterogeneity on traveltime shifts for compacting reservoirs     

Rodrigo F. Fuck  Ilya Tsvankin  Andrey Bakulin 《Geophysical Prospecting》2011,59(1):78-89

Compaction induced by pore‐pressure decrease inside a reservoir can be monitored by measuring traveltime shifts of reflection events on time‐lapse seismic data. Recently we introduced a perturbation‐based formalism to describe traveltime shifts caused by the 3D stress‐induced velocity field around a compacting reservoir. Application of this method to homogeneous background models showed that the offset variation of traveltime shifts is controlled primarily by the anisotropic velocity perturbations and can provide valuable information about the shear and deviatoric stresses. Here, we model and analyse traveltime shifts for compacting reservoirs whose elastic properties are different from those of the surrounding medium. For such models, the excess stress is influenced primarily by the contrast in the rigidity modulus μ across the reservoir boundaries. Synthetic examples demonstrate that a significant (25% or more) contrast in μ enhances the isotropic velocity perturbations outside the reservoir. Nevertheless, the influence of background heterogeneity is mostly confined to the reservoir and its immediate vicinity and the anisotropic velocity changes are still largely responsible for the offset dependence of traveltime shifts. If the reservoir is stiffer than the host rock, the background heterogeneity reduces anisotropic velocity perturbations inside the reservoir but increases them in the overburden. As a result, in this case, the magnitude of the offset variation of traveltime shifts is generally higher for reflections from interfaces above the reservoir. We also study compaction‐induced stress/strain and traveltime shifts for a stiff reservoir embedded in a softer layered model based on velocity profiles from the Valhall Field in the North Sea. Despite producing discontinuities in strain across medium interfaces, horizontal layering does not substantially alter the overall behaviour of traveltime shifts. The most pronounced offset variation of traveltime shifts is observed for overburden events recorded at common midpoints close to the reservoir edges. On the whole, prestack analysis of traveltime shifts should help better constrain compaction‐induced velocity perturbations in the presence of realistic background heterogeneity.  相似文献   

Chemical fluid transport (CFT): A window into Earth and its development     

Werner Schrön 《Chemie der Erde / Geochemistry》2013

Gas phase transport according to chemical fluid transport (CFT) in Earth's crust as well as in the solar nebula is characterized by very high transport efficiency. Systematic investigations of mobilization, transport and deposition of gaseous MeX (Me = metal, X = F or Cl) compounds by solid gas equilibrium reactions are suitable to explain numerous extensive accumulations of minerals and ores. More than 40 of the considered chemical elements form volatile MeX compounds. Some elements tend to form MeF compounds, whereas others are more likely to form MeCl compounds. Silicon reacts with HF to form SiF₄ and replaces other elements to form MeF compounds at low temperature ranges. Accumulations caused by SiF₄ transport explain the formation of numerous quartz varieties and silicate minerals in Earth's crust. Iron most likely reacts with HCl to form FeCl₂ as well as FeCl₃ and explain the formation of iron or iron compounds. Thermodynamically directed transport from cool to hot areas in connection with cyclic processes increases the transport efficiency of MeX-species. Such species are SiF₄, Al₂F₆, POF₃, Cu₃Cl₃, SnCl₄, BF₃, GeF₄, GeCl₄, Ga₂Cl₆, ZrF₄, NbF₅ and TiF₄. The transport gases SiF₄ and POF₃ often react with environmental compounds forming pneumatolytic and metasomatical mineral accumulations. CFT is the “motor” of pneumatolytic and metasomatical processes.  相似文献   

The coral reef crisis: The critical importance of <350 ppm CO₂     

J.E.N. Veron  O. Hoegh-Guldberg  C.R.C. Sheppard  M.G. Stafford-Smith 《Marine pollution bulletin》2009,58(10):1428-1436

Temperature-induced mass coral bleaching causing mortality on a wide geographic scale started when atmospheric CO₂ levels exceeded ∼320 ppm. When CO₂ levels reached ∼340 ppm, sporadic but highly destructive mass bleaching occurred in most reefs world-wide, often associated with El Niño events. Recovery was dependent on the vulnerability of individual reef areas and on the reef’s previous history and resilience. At today’s level of ∼387 ppm, allowing a lag-time of 10 years for sea temperatures to respond, most reefs world-wide are committed to an irreversible decline. Mass bleaching will in future become annual, departing from the 4 to 7 years return-time of El Niño events. Bleaching will be exacerbated by the effects of degraded water-quality and increased severe weather events. In addition, the progressive onset of ocean acidification will cause reduction of coral growth and retardation of the growth of high magnesium calcite-secreting coralline algae. If CO₂ levels are allowed to reach 450 ppm (due to occur by 2030-2040 at the current rates), reefs will be in rapid and terminal decline world-wide from multiple synergies arising from mass bleaching, ocean acidification, and other environmental impacts. Damage to shallow reef communities will become extensive with consequent reduction of biodiversity followed by extinctions. Reefs will cease to be large-scale nursery grounds for fish and will cease to have most of their current value to humanity. There will be knock-on effects to ecosystems associated with reefs, and to other pelagic and benthic ecosystems. Should CO₂ levels reach 600 ppm reefs will be eroding geological structures with populations of surviving biota restricted to refuges. Domino effects will follow, affecting many other marine ecosystems. This is likely to have been the path of great mass extinctions of the past, adding to the case that anthropogenic CO₂ emissions could trigger the Earth’s sixth mass extinction.  相似文献   

Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures     

Hussein Hoteit  Abbas Firoozabadi 《Advances in water resources》2008

Contrast in capillary pressure of heterogeneous permeable media can have a significant effect on the flow path in two-phase immiscible flow. Very little work has appeared on the subject of capillary heterogeneity despite the fact that in certain cases it may be as important as permeability heterogeneity. The discontinuity in saturation as a result of capillary continuity, and in some cases capillary discontinuity may arise from contrast in capillary pressure functions in heterogeneous permeable media leading to complications in numerical modeling. There are also other challenges for accurate numerical modeling due to distorted unstructured grids because of the grid orientation and numerical dispersion effects. Limited attempts have been made in the literature to assess the accuracy of fluid flow modeling in heterogeneous permeable media with capillarity heterogeneity. The basic mixed finite element (MFE) framework is a superior method for accurate flux calculation in heterogeneous media in comparison to the conventional finite difference and finite volume approaches. However, a deficiency in the MFE from the direct use of fractional flow formulation has been recognized lately in application to flow in permeable media with capillary heterogeneity. In this work, we propose a new consistent formulation in 3D in which the total velocity is expressed in terms of the wetting-phase potential gradient and the capillary potential gradient. In our formulation, the coefficient of the wetting potential gradient is in terms of the total mobility which is smoother than the wetting mobility. We combine the MFE and discontinuous Galerkin (DG) methods to solve the pressure equation and the saturation equation, respectively. Our numerical model is verified with 1D analytical solutions in homogeneous and heterogeneous media. We also present 2D examples to demonstrate the significance of capillary heterogeneity in flow, and a 3D example to demonstrate the negligible effect of distorted meshes on the numerical solution in our proposed algorithm.  相似文献   

Stochastic analysis of bounded unsaturated flow in heterogeneous aquifers: Spectral/perturbation approach     

Ching-Min Chang  Hund-Der Yeh 《Advances in water resources》2009

This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84].  相似文献   

摩擦时间依从的地震活动性细胞自动机模型   总被引：11，自引：1，他引：10       下载免费PDF全文

刘桂萍  傅征祥  刘杰 《地球物理学报》2000,43(2):203-212

设计了一个改进的单断层地震孕育过程细胞自动机（ＣＡ）模型,通过设计外界通过 施加应力与模型间进行的能量交换和模型的细胞之间存在的非线性力学作用,试图理解地震 活动特性的力学机制．与早期细胞自动机模型相比,改进了参数的取值方式,将摩擦时间依 从的理论引进模型,使该细胞自动机模型更接近实际的孕震系统．研究表明参数取值方式对 人工地震序列和各细胞破裂事件的非均匀时间特性有重要影响,较小震级和较大震级范围中 的事件分别遵从明显不同的累积频度一震级关系．  相似文献   

非均匀套管井中的过套管电阻率测井响应   总被引：13，自引：2，他引：11       下载免费PDF全文

高杰  刘福平  包德洲  程希  杨长春 《地球物理学报》2008,51(4):1255-1261

套管井壁的电势分布包含了金属套管和地层模型的信息,套管电阻率的变化对过套管电阻率测井有重要影响,本文考虑了套管电阻率的变化对传输线方程系数的影响,给出了非均匀套管壁电势分布的计算方法．利用套管井壁电势分布实现了非均匀套管测井响应的计算,考察了套管电阻率变化对过套管电阻率测井响应的影响,模拟结果表明,套管的变薄或变厚对测井响应均有一定影响,当套管电阻突然变小（套管变厚）或突然变大（套管变薄）时过套管电阻率测井异常发生在套管电阻突变界面附近．模拟结果还显示电极距对测井结果有明显影响．  相似文献   

Review of Methods for Simulating Land Surface Temperature at the Pixel Scale Based on Ground Measurements over Heterogeneous Surface     

Peng Zhixing  Zhou Ji  Li Mingsong 《地球科学进展》2016,31(5):471-480

Remotely sensed Land Surface Temperature (LST) is a key parameter for studying the global climate changes and the exchanges of water and energy. Acquiring LST accurately is important to diagnose the change of environment on earth. Quantifying the uncertainty of remotely sensed LST is the first step of its application. However, due to the difficulties in obtaining the ground truth of LST at the pixel scale, it is difficult to validate the remotely sensed LST. Here, methods for simulating the LST at the pixel scale based on ground measurements over heterogeneous area were reviewed. From the way to construct the ground scene, these methods were classified into three types, including the Modified Geometric Projection model (MGP), realistic structural three-dimensional model, and other model. The advantages and disadvantages of these models were examined and compared. Finally, some issues in simulating LST at the pixel scale over heterogeneous area needed to be solved and on-going directions in the future were summarized.  相似文献   

Physiological response to high temperature in the Tropical Eastern Pacific coral Pocillopora verrucosa          下载免费PDF全文

Alma P. Rodríguez‐Troncoso  Eugenio Carpizo‐Ituarte  Amílcar L. Cupul‐Magaña 《Marine Ecology》2016,37(5):1168-1175

The ecological and economic values of coral reef communities have encouraged efforts to implement periodic visual survey programs to secure their conservation. To date, visual monitoring‐based approaches have detected bleaching events weeks or months after the initial onset. An evaluation of the stress response of coral colonies, as well as their ability to resist and recover from the stress events, may increase our understanding of the physiological processes underling the stress and/or acclimation responses. Coral bleaching events, caused primarily by abnormally high temperatures, are continuously affecting coral communities worldwide. To evaluate the stress and recovery responses of the reef‐building coral Pocillopora verrucosa, a laboratory study was conducted herein. Coral nubbins were exposed to high temperatures to induce bleaching, and their ability to resist and recover from bleaching was subsequently monitored after returning the corals to ambient temperature. Lipid and chlorophyll concentrations, as well as Symbiodinium spp. density, decreased in samples exposed to nearly 31°C, slowly recovering to near‐control levels upon return to non‐stressful conditions. The present data set allows for an assessment of the vulnerability and ability of this common species to resist sub‐optimal environmental conditions, particularly the thermal stress events that will occur more commonly in their habitats as global seawater temperatures continue to rise.  相似文献