High-resolution simulation and characterization of density-driven flow in CO2 storage in saline aquifers

Simulations are routinely used to study the process of carbon dioxide (CO₂) sequestration in saline aquifers. In this paper, we describe the modeling and simulation of the dissolution–diffusion–convection process based on a total velocity splitting formulation for a variable-density incompressible single-phase model. A second-order accurate sequential algorithm, implemented within a block-structured adaptive mesh refinement (AMR) framework, is used to perform high-resolution studies of the process. We study both the short-term and long-term behaviors of the process. It is found that the onset time of convection follows closely the prediction of linear stability analysis. In addition, the CO₂ flux at the top boundary, which gives the rate at which CO₂ gas dissolves into a negatively buoyant aqueous phase, will reach a stabilized state at the space and time scales we are interested in. This flux is found to be proportional to permeability, and independent of porosity and effective diffusivity, indicative of a convection-dominated flow. A 3D simulation further shows that the added degrees of freedom shorten the onset time and increase the magnitude of the stabilized CO₂ flux by about 25%. Finally, our results are found to be comparable to results obtained from TOUGH2-MP.     

Comment on the effect of anisotropy on the onset of convection in a porous medium

The effect of anisotropy on the onset of convection in a saturated porous medium is discussed. In particular, the case of time-dependent density-driven convection is examined. The applicability of the value of an equivalent Rayleigh number as the criterion for the onset of convection is discussed.     

An extended stability criterion for density-driven flows in homogeneous porous media

Stability of density-driven flows is a challenging problem with current applications in major areas like energy exploration, water pollution, nuclear and oil industries. The mathematical model for such flows is a system of coupled non linear partial differential equations. To study the physical stability of the system, we consider steady-state flow and perturb the solution of the full system of equations (without Boussinesq approximation) and investigate how it evolves in time: if the solution does not grow indefinitely, the system is called stable. The perturbations are treated as being the result of sub-scale interactions between the velocity field and the solute mass. Making use of a two-scale expansion of the solution, we derived extended stability criteria that include the effects of density, viscosity and flow velocity in flow configurations aligned parallel as well as orthogonal to gravity forces. Numerical simulations with the numerical simulator d³f$d^{3}f$ are presented to test the theoretical stability criteria.     

The stability of density-driven flows in saturated heterogeneous porous media

This work concludes the investigations into the stability of haline flows in saturated porous media. In the first part [33] a stability criterion for density-driven flow in a saturated homogeneous medium was derived excluding dispersion. In the second part [34], the effects of dispersion were included. The latter criterion made reasonable predictions of the stability regimes (indicated by the number of fingers present) as a function of density and dispersivity variations. We found out that destabilising variables caused an increase in the number of fingers and vice versa. The investigation is extended here for the effects of the medium heterogeneity. The cell problem derived via homogenization theory [20] is solved and its solution used to evaluate the elements of the macrodispersion tensor as functions of time for flow aligned parallel to gravity. The longitudinal coefficient exhibits asymptotic behaviour for favourable and moderately unfavourable density contrasts while it grows indefinitely for higher density contrasts. The range of densities stabilised by medium heterogeneities can thus be estimated from the behaviour of the coefficient. The d³f software program is used for the numerical simulations. The code uses the cell-centred finite volume and the implicit Euler techniques for the spatial and temporal discretisations respectively.     

Optimal design of pumping tests in leaky aquifers for stream depletion analysis

We analyze the optimal design of a pumping test for estimating hydrogeologic parameters that are subsequently used to predict stream depletion caused by groundwater pumping in a leaky aquifer. A global optimization method is used to identify the test’s optimal duration and the number and locations of observation wells. The objective is to minimize predictive uncertainty (variance) of the estimated stream depletion, which depends on the sensitivities of depletion and drawdown to relevant hydrogeologic parameters. The sensitivities are computed analytically from the solutions of Zlotnik and Tartakovsky [Zlotnik, V.A., Tartakovsky, D.M., 2008. Stream depletion by groundwater pumping in leaky aquifers. ASCE Journal of Hydrologic Engineering 13, 43–50] and the results are presented in a dimensionless form, facilitating their use for planning of pumping test at a variety of sites with similar hydrogeological settings. We show that stream depletion is generally very sensitive to aquitard’s leakage coefficient and stream-bed’s conductance. The optimal number of observation wells is two, their optimal locations are one close to the stream and the other close to the pumping well. We also provide guidelines on the test’s optimal duration and demonstrate that under certain conditions estimation of aquitard’s leakage coefficient and stream-bed’s conductance requires unrealistic test duration and/or signal-to-noise ratio.     

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

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].     

Acoustic emission and volumetric strain induced in coal by the displacement sorption of methane and carbone dioxide

The aim of this study was to assess whether acoustic emission (AE) could carry information on preferential sorption/desorption of CH₄ or CO₂ in coal. AE and expansion/contraction of two nearly identical cylindrical coal samples were continuously monitored during displacement sorption experiments. One sample was subjected to presorption of CH₄, followed by sorption of CH₄/CO₂ mixture. With the other one, presorption of CO₂ preceded sorption of the mixture. The results obtained are the following: first, AE and stain kinetics show that the affinity of the coal tested is higher for CO₂ than for CH₄; second, methane is preferentially desorbed after presorption of CH₄ — sorption of mixture of CH₄ and CO₂; third, during displacement sorption, kinetics of AE and sample swelling/shrinkage bring out the importance of presorption and the sorbate used. It matters whether the coal is first exposed to CH₄ or to CO₂. The present study has demonstrated that injection of CO₂ into the coal previously exposed to CH₄ causes considerable swelling of the coal. On desorption after CH₄/CO₂ exchange sorption, initial shrinkage is followed by swelling of the coal. These results could have implications for the sequestration of CO₂ in coal seams and CH₄ recovery from coalbeds (ECBM). Swelling/shrinkage of the coal matrix should be included in models used to predict coal permeability and gas flow rates. They also show that the AE technique can give more insights into coal matrix-gas interactions.     

地震作用下边坡稳定性分析方法

地震引起的滑坡已成为大地震引起的严重次生灾害现象,如2008年汶川地震引起的滑坡灾害．本文对地震作用下边坡稳定性分析的研究现状及特点进行了分析总结．地震边坡稳定性分析方法总体上可分为定性分析和定量分析两类．本文分别从工程地质与灾害调查分析、 灾害调查与地震烈度及地震动参数的统计分析、 简化模型理论分析、 数值计算分析和模型试验分析等方面,介绍了相应边坡稳定性分析理论和方法的提出及其沿革,分析了不同理论和方法的特点及适应性,并进一步提出了相关研究存在的问题和发展趋势．     

Post‐stack stratigraphic inversion workflow applied to carbon dioxide storage: application to the saline aquifer of Sleipner field‡

In the Norwegian North Sea, the Sleipner field produces gas with a high CO₂ content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO₂) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO₂ plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO₂ over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO₂ plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO₂ injection and the 2006 vintage carried out after a CO₂ injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO₂ in a reservoir. With the post‐stack inversion workflow applied to CO₂ storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.     

Evaluating the impact of soil redistribution on the in situ mineralization of soil organic carbon

Soil erosion, transport and deposition by water drastically affect the distribution of soil organic carbon (SOC) within a landscape. Moreover, soil redistribution may have a large impact on the exchange of carbon (C) between the pedosphere and the atmosphere. One of the large information gaps within this research domain, concerns the fate of SOC after erosion by water. According to different (mainly laboratory) studies, soil redistribution leads to aggregate breakdown, thereby exposing the contained SOC to mineralization. Our study aims to quantify the extent to which such increased mineralization occurs in a real field situation. Carbon dioxide (CO₂)‐efflux was measured in the field after an important erosion event for a continuous period of 112 days. The specific situation on the field ensured that almost none of eroded SOC was exported from the field. Measurements of CO₂‐efflux were done in areas with sediment deposition, as well as in comparable areas without sedimentation. Comparison of these measurements allowed the net effect of soil deposition on CO₂‐efflux to be assessed. Field data were complemented by measurements on incubated, undisturbed soil core samples, in order to disentangle the contribution of environmental factors (moisture, temperature) from any erosional effect on CO₂‐efflux. Results of these measurements on the field showed that CO₂‐efflux was regulated by a complex interplay of different factors (mostly soil porosity, soil moisture and soil temperature). In combination with the incubation measurements, it could be concluded that the processes of erosion and transport indeed led to an increased mineralization of SOC, as a result of aggregate breakdown and exposure of previously encapsulated SOC. This effect was, however, much smaller than observed in previous laboratory studies. Moreover, it was only important in the first weeks, immediately after the erosion event. The calculated net erosional effect on CO₂‐efflux represented a mere 1·6% of total SOC, originally present in the soil. Copyright © 2010 John Wiley & Sons, Ltd.     

60年来大气中二氧化碳浓度数据的趋势方程研究

通过讨论已有的60年来大气中CO₂浓度数据的分布状态,采用趋势分析方法,给出了具体趋势方程形式.与冰芯分析或观测数据对比结果表明,趋势方程曲线与已有数据基本符合,随后初步给出了2010年至2016年间大气中CO₂浓度预测值.     

国产磁通门磁力仪稳定性分析

国产磁通门磁力仪(GM-3磁通门磁力仪)在全国各地磁台站推广应用,两年来,该仪器在武汉地震台工作比较稳定。本文对该仪器的基线值进行了多层线性回归分析。结果表明仪器性能比较稳定。     

地震边坡稳定性分析的拟静力方法适用性探讨

应用极限平衡拟静力方法进行地震边坡稳定性分析时,对粘聚力较小、内摩擦角较大的浅层滑坡而言,随着地震水平设计加速度值的增大,其临界滑动面基本保持不变;然而,对粘聚力较大而摩擦角较小的深层滑坡而言,随着地震水平设计加速度值的增大,临界滑动面越深、越宽,会呈现出发散现象,选取两个典型的均质土坡,设计了由浅及深的多条滑动面,分别采用极限平衡瑞典方法进行了相应安全系数的计算,同时比较了自重、地震力、粘聚力在抗滑力矩以及滑动力矩中所占比重的变动情况。研究发现:对于浅层滑坡,随着滑动面由浅及深的变化,粘聚力产生的抗滑力矩所占百分比下降幅度较大,自重产生的滑动力矩与总体抗滑力矩的比值下降幅度较大;而对于深层滑坡,随着滑动面由浅及深变化,自重产生的抗滑力矩所占百分比增长的幅度较大,地震力产生的滑动力矩与总体抗滑力矩比值增长幅度也较大。     

动力反应分析的显式积分方法及其稳定性

Newm ark-更新精细积分法是动力方程求解的隐式的时域逐步积分法,其稳定性条件非常容易满足。与隐式方法相比较,显式积分方法不需要求解耦联的方程组,可以有效地减少内存占用和机时耗费。因此,根据显式积分方法的特点和优点,基于Newm ark-更新精细积分法的基本思想,提出其显式积分格式。对显式积分方法的精度与稳定性进行了初步的分析,指出该显式积分方法具有极好的稳定性,其精度比隐式积分方法的精度稍低。随着时间步长的增加,其精度优于传统的方法。     

Molecular and carbon isotopic compositions of gas inclusions of deep carbonate rocks in the Tarim Basin

Gaseous components of gas inclusions in deep carbonate rocks (>5700 m) from the Tacan 1 well were analyzed by online mass spectrometry by means of either the stepwise heating technique or vacuum electromagnetism crushing. The carbon isotopic compositions of gases released by vacuum electromagnetism crushing were also measured. Although the molecular compositions of gas inclusions show differences between the two methods, the overall characteristics are that gas inclusions mainly contain CO2, whilst hydrocarbon gases, such as CH4, C2H6 and C3H8, are less abundant. The content of CO is higher in the stepwise heating experiment than that in the method of vacuum electromagnetism crushing, and there are only minor amounts of N2, H2 and O2 in gas inclusions. Methane δ13C values of gas inclusions in Lower Ordovician and Upper Cambrian rocks (from 5713.7 to 6422 m; -52‰-63‰) are similar to those of bacterial methane, but their chemical compositions do not exhibit the dry character in comparison with biogenic gases. These characteristics of deep gas inclusions may be related to the migration fractionation. Some deep natural gases with light carbon isotopic characteristics in the Tazhong Uplift may have a similar origin. The δ13C1 values of gas inclusions in Lower Cambrian rocks (7117-7124 m) are heavier (-39‰), consistent with highly mature natural gases. Carbon isotopic compositions of CO2 in the gas inclusions of deep carbonate rocks are similar (from -4‰ to -13‰) to those of deep natural gases, indicating predominantly an inorganic origin.     

The effect of dispersion on the stability of density-driven flows in saturated homogeneous porous media

This work is a continuation of Musuuza et al. [37] in which a stability criterion for density-driven flow in a saturated homogeneous medium was derived. The criterion predicted the stability of a system as a function of the density and viscosity contrasts, the magnitude of the flow velocity and the concentration gradients for flow aligned orthogonal to gravity. It could not accurately predict stability transition with increasing velocity, a failure we attributed to dispersion effects that were not included. Small-scale dispersion and molecular diffusion, the main stabilising mechanisms in homogeneous media can stabilise flow if the instability wavelengths are below a certain cutoff. The width of the mixing zone is also central in controlling the range of wavelengths that persist into fingers. We propose a method of quantifying the cutoff wavelength and the width of the mixing zone, which are incorporated into the earlier criterion as constituents of the dispersive part. The earlier criterion is reformulated in terms of the Rayleigh number and with the dispersive part added, we attempt to predict the number of fingers formed which is directly related to the physical stability of the system. The inclusion of the Rayleigh number and dispersion into a single stability criterion provides new insight in the way dispersion affects vertical flow systems. Stability numbers computed with the new criterion are in reasonable agreement with numerical simulations for a range of physical variables. The numerical computations are performed with the software package d³f, which uses the cell-centred finite volume and the implicit Euler methods for the spatial and temporal discretisations, respectively. The admission of the density and dispersivities as inputs into the criterion makes it usable in practical problems.     

Hydrologic and geologic control of carbonate water chemistry in the subarctic nahanni karst,canada

There are six classes of water and five geologic environments in the subarctic Nahanni karst. During the summers of 1972 and 1973, 214 water samples were collected from 15 of the 30 hydrogeologic categories. Linear discriminant function analysis, using five measured and two derived chemical variables, indicates that there are statistically significant (0.005 level) differences in water chemistry between similar waters in different geologic environments, between waters in the same geologic environment, and between waters in different hydrogeological categories. Geological environment labels a natural water because it determines the availability of soil CO₂ and of soluble minerals. Measurements indicate that mean soil log PCO₂ is greatest in areas of shale mantled by till (?2.39), and least in areas of sandy fluvioglacial drift (?3.27). Low values on the sandy drift are due to the sparse shrub vegetation, and to the high degree of soil aeration; soils in areas of shale are clay-rich and support a dense boreal forest. Hydrology influences water chemistry because it determines how much CO₂ natural waters pick up from the environment and how much they subsequently lose to the atmosphere, and as a result, whether they dissolve or deposit soluble materials. The similarity between mean calculated log PCO₂ in natural waters (?2.92) and mean measured soil log PCO₂ (?2.80) suggests that natural waters in Nahanni are dose to equilibrium with mean soil CO₂ levels.     

广东新丰江水库表层水体CO2分压及其影响因素

于2012年7月和2013年1月定点采集新丰江水库表层水样,测定水样的理化及生物学参数,计算水体中二氧化碳分压(p(CO_2))大小并分析其时空变化,探讨新丰江水库p(CO_2)的影响因素及其CO_2源/汇机制.结果表明:丰水期p(CO_2)变化范围为16~3545μatm,均值为999μatm,从水库上游到坝前p(CO_2)逐渐升高;枯水期p(CO_2)变化范围为399~1355μatm,均值为756μatm,从水库上游到坝前p(CO_2)呈下降趋势.丰水期p(CO_2)受温度影响较小,与营养盐(NO-3、DSi)浓度呈正相关,与叶绿素a(Chl.a)、溶解有机碳(DOC)浓度呈负相关,与溶解无机碳(DIC)浓度没有明显相关性;枯水期p(CO_2)受温度影响也较小,受碳酸盐体系的影响,与NO-3、DSi、Chl.a、DIC浓度呈正相关,与DOC浓度没有明显相关性.新丰江水库相对于大气来说是一个通量值偏低的CO_2源.     

Long‐term carbon storage and hydrological control of CO2 exchange in tundra ponds in the Hudson Bay Lowland

Carbon dioxide fluxes and water balance were examined in 43 tundra ponds in the northern portion of the Hudson Bay Lowland near Churchill, Manitoba. Most of the ponds were hydrologically disconnected from their catchments during dry periods throughout the post‐melt season. However, episodic reconnection occurred following large precipitation events where depression storage was exceeded. Significant shifts in pond chemistry were observed following precipitation events, with the degree of CO₂ saturation increasing during these periods. Pond CO₂ concentrations rapidly fell to pre‐event levels following events, suggesting that hydrological connectivity can affect the magnitude and direction of CO₂ gas fluxes in tundra ponds. Atmospheric CO₂ invaded ponds with highly organic sediments for most of the summer, suggesting that terrestrially derived inorganic carbon was insufficient to meet the demands of algal net production. In contrast, ponds with highly mineral sediments continued to evade CO₂ during the summer. In a subset of 11 ponds, long‐term rates of carbon accumulation in sediment ranged from 0·6 to 2·2 mol C m^?2 year^?1. Very strong correlations existed between average sediment accumulation rates and pond perimeters and basin areas suggesting that peat may be a major source of sediment carbon. Aeolian transport is also a potentially large source of sediment carbon. Copyright © 2004 John Wiley & Sons, Ltd.