The effects of carbon dioxide leakage on fractures, and water quality of potable aquifers during geological sequestration of CO2

Since industrial revolution, the ＂greenhouse effect＂ is one of the most important global environmental issues. Of all the greenhouse gases, CO2 is responsible for about 64% of the enhanced ＂greenhouse effect＂, making it the target for mitigation, so reducing anthropogenic discharge of carbon dioxide attracts more and more attention. Geological sequestration of CO2 in deep saline aquifers is one of the most promising options. But because unknown fractures and faults may exist in the caprock layers which can prevent the leakage of CO2, CO2 will leak upward into upper potable aquifers, and lead to adverse impacts on the shallow potable aquifers. In order to assess the potential effect of CO2 leakage from underground storage reservoirs on fractures and water quality of potable aquifers, this study used the non-isothermal reactive geochemical transport code TOUGHREACT developed by Xu et al to establish a simplified 2-D model of CO2 underground sequestration system, which includes deep saline aquifers, caprock layers, and shallow potable aquifers, and study and analyze the changes of mineral and aqueous components. The simulation results indicated that the minerals of deep saline aquifers and fractures should be mainly composed of aluminosilicate and silicate minerals, which not only enhance the mass of CO2 sequestrated by mineral trapping, but also decrease the porosity and permeability of caprock layers and fractures to prevent and reduce CO2 leakage. The results from deep saline aquifers showed that the mass of carbon dioxide trapped by minerals and solution phases is limited, the rest remained as a supercritical phase, and so once the caprock aquifers have some unknown fractures, the free carbon dioxide phase may leak from CO2 geologic sequestration reservoirs by buoyancy.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Modeling of the classification behavior of the diaphragms used in multi-chamber cement mills

A new approach has been developed for describing the classification that arises at the diaphragm in multi-chamber ball mills. Fourteen industrial scale cement mills have been sampled in this study and after mass balancing, model fits of the equipment were done. Successive mill product predictions for the sampled conditions were achieved, and requirements to build up an optimization capable model are discussed.     

虚拟样机技术在地质工程机械领域的应用研究

基于对国内地质工程机械的研发现状的考虑，根据虚拟样机技术的特点，讨论了在地质工程机械领域实施虚拟样机技术的必要件，并提出了建立地质工程机械虚拟样机的策略。联合PRO／E、ANSYS及ADAMS建立了气动潜孔锤的虚拟样机模型，对气动潜孔锤虚拟样机进行了仿真，阐明了虚拟样机技术在地质工程机械领域的应用前景。     

Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.     

Assessment of mercury mobility and bioavailability by fractionation method in sediments from coastal zone inundated by the 26 December 2004 tsunami in Thailand

The 26 December 2004 tsunami covered significant portion of a coastal zone with a blanket of potentially contaminated sediments. In this report are presented results on mercury concentrations in sediments deposited by the tsunami in a coastal zone of Thailand. Since the total mercury concentrations are insufficient to assess mercury mobility and bioavailability in sediment, its fractionation was applied. Sediments were sampled within 50 days after the event and analyzed by sequential extraction method. The procedure of sequential extraction involved five subsequent stages performed with solutions of chloroform, deionized water, 0.5 M HCl, 0.2 M NaOH, and aqua regia. The mean concentration of total mercury in sediments was 119 ± 50 ng g⁻¹ dry mass (range 66–230). The fractionation revealed that mercury is mainly bound to the least bioavailable sulphides 75 ± 6% (range 62–86), organomercury compounds 14 ± 7% (range 4–26), and humic matter 9 ± 7% (range 1–27). The lowest contributions bring fractions of water-soluble mercury 0.8 ± 1.0% (range 0.1–3.6) and acid soluble mercury 0.9 ± 0.5% (range 0.2–2.1). Although, the total mercury content is similar in a reference sample and in the tsunami sediments, the highly toxic organomercury fraction contribution is higher in the latter. The results were compared with chemical and sedimentological properties of the sediments but no significant correlations were obtained between them.     

Atomistic calculations of structural and elastic properties of serpentine minerals: the case of lizardite

The physical properties of the hydrous phyllosilicate lizardite have been investigated by atomistic simulation using the GULP code based on transferable semi-empirical interatomic potentials. Lizardite behavior was first investigated during structure relaxation at room temperature. The Helmholtz free energy is minimum for an equilibrium structure that is in agreement with experiment. The bulk, shear, and Young modulii for lizardite were calculated along with the Poisson ratio. From the shear and bulk modulii, we also calculated translational and longitudinal acoustic wave velocities that are important quantities for tectonophysics models. As expected, lizardite is stiffer in the a direction parallel to the layers than in the c perpendicular direction; the variation of the unit cell parameters with pressure is in good agreement with experiment. The cohesive energy between two successive layers along c direction was calculated at 0.33 eV (i.e., 0.11 eV per OH bond) in good agreement with recent ab initio calculations. Upon pressure and temperature variations, we evidenced that structural changes are mainly pressure induced; pressure being accommodated by a decrease of the c parameter up to 10 GPa. We also found that the change of slope in the derivative of the c cell parameter with respect to pressure occurring around 2 GPa originates from the bending of the interlayer hydroxyl groups with respect to the layer normal direction.     

Distribution of cadmium, chromium, copper, lead and zinc in marine sediments in Hong Kong waters

Partitioning of heavy metals (Cd, Cr, Cu, Pb, Zn) in marine sediments collected from various sites in Hong Kong waters were determined using sequential extraction method. Sediments from Kellette Bank, located in Victoria Harbour, had higher metal concentrations especially Cu and Zn than most other sites. Slightly over 20% of total Cu and Cr existed as readily available forms in Peng Chau and Kellette Bank. At most sampling sites, over 15% of the Cu existed as the exchangeable form indicating that Cu could be readily released into the aqueous phase from sediments. A significantly higher percentage of Pb and Zn was associated with the three non-residual fractions. Hence, there is a greater environmental concern for remobilization of Pb and Zn compared with Cr. The high amount of residual Cd (>50%) and the relatively lower Cd content indicate that little environmental concern is warranted for the remobilization of Cd. Distribution of metals in sediments collected from different depth at Kellette Bank shows that metal concentrations decreased with profile depth. The levels of Pb and Zn associated with the two readily available fractions increased sharply in the surface sediment. These metals represented the pollutants, which were introduced into the area in the mid-eighties through early nineties as a result of rapid economic and industrial development in the territory. As significant portions of these metals were bound to the readily available phases in the surface sediments, metal remobilization could be a concern. An erratum to this article can be found at     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

考虑滑脱效应和温度场的煤层气渗流数值模拟

由于煤层气的解吸热效应,煤层气的运移过程是一个非等温过程。因此,温度场对煤层气渗流有着重要的影响。前人在研究煤层气渗流规律的研究中并没有同时考虑滑脱效应和温度场因素的影响,而在实际的深部开采中综合考虑滑脱效应和温度场的影响对研究深部煤层气运移规律有着重要的意义。因此本文建立了考虑滑脱效应及温度场的煤层气渗流数学模型,利用有限元数值方法研究了考虑滑脱效应和温度场耦合的煤层气运移规律;研究了考虑滑脱效应和温度场对压力场分布的影响;对考虑滑脱效应及温度场因素的煤层气产量进行了预测。得出了随着温度的升高煤储层的压力在降低,温度的升高对煤层气的产量有着负面的影响这一重要结论。