二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

Gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, Northern South China Sea

Based on the comprehensive interpretation and study of the Neogene fracture system and diapiric structure, it can be concluded that the diapiric structures, high-angle fractures and vertical fissure system are the main gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, northern South China Sea. The Neogene fractures widely developed in the study area may be classed into two groups: NW (NNW)-trending and NE (NNE)-trending. The first group was active in the Late Miocene, while the second one was active since the Pliocene. The NE (NNE)-trending fractures were characterized by lower activity strength and larger scale, and cut through the sediment layers deposited since the Pliocene. Within the top sediment layers, the high-angle fracture and vertical fissure system was developed. The diapiric structures display various types such as a turtle-back-like arch, weak piercing, gas chimney, and fracture (or crack, fissure). On the seismic profile, some diapiric structures show the vertical chimney pathway whose top is narrow and the bottom is wide, where some ones extend horizontally into pocket or flower-shaped structures and formed the seismic reflection chaotic zones. Within the overlying sediment layers of the diapiric structure, the tree branch, flower-shaped high-angle fractures and vertical fissures were developed and became the pathway and migration system of the gas-bearing fluid influx. In the study area, the diapiric structures indicate a high temperature/over pressure system ever developed. Closely associated and abundant bright-spots show the methane-bearing fluid influx migrated vertically or horizontally through the diapiric structures, high-angle fractures and vertical fissures. In the place where the temperature and pressure conditions were favor for the formation of gas hydrate, the hydrate reservoir deposition sub-system was developed.

Determination of gold, platinum and palladium in geological samples by inductively coupled plasma-mass spectrometry after separation and preconcentration with C-410 anion exchange resin

A method for the determination of Au, Pt and Pd in geological samples is described. Au, Pt and Pd can be separated and concentrated quantitatively by C-410 anion-exchange resin in the condition of 1.5 mol/L HCl with the adsorption rates of 91.2%, 100.0% and 95.7% respectively. No interference exists from coexisting elements except for Ge(IV), Cr(VI),Ti(IV) in inductively coupled plasma-mass spectrometry. The detection limits are 0.27 mug/L, 0.40 mug/L and 0.19 mug/L for Au, Pt and Pd respectively. The results of these elements in standard geological materials are in agreement with certified values with precision of 19.2% RSD for Au (n = 8), 28.1% RSD for Pt (n=8), and 15.6% RSD for Pd (n=8).

Tectonic geomorphology of the Ryukyu Trench-Arc-Backarc System: geological-geophysical exploration and mapping

OKINAWA TROUGH; BASIN

Recent 2000-year geological records of mud in the inner shelf of the East China Sea and their climatic implications

AXIS(14)C dating and grain-size analysis for Core DD2, located at the north of the Yangtze River-derived mud off the Zhejiang-Fujian coasts in the inner shelf of the East China Sea, provide us a high-resolution grain-size distribution curve varying with depth and time. Data in the upper mud layer of Core DD2 indicate that there are at least 9 abrupt grain-size increasing in recent 2000 years, with each corresponding very well with the low-temperature events in Chinese history, which might result from the periodical strengthening of the East Asian Winter Monsoon (EAWM), including the first-revealed maximum temperature lowering event at around 990 a BP. At the same time, the finer grain size section in Core DD2 agrees well with the Sui-Tang Warming Period (600-1000 a AD) defined previously by Zhu Kezhen, during which the climate had a warm, cold and warm fluctuation, with a dominated cooling period of 750-850 a AD. The Little Ice Age (LIA) can also be identified in the core. It starts around 1450 a AD and was followed by a subsequent cooling events at 1510, 1670 and 1840 a AD. Timing of these cold events revealed here still needs to be further verified owing to some current uncertainty of dating we used in this study.

Simultaneous determination of four phases of copper in geological samples

According to the specific property of respective phases in copper ore, optimal extracting condition for, tour phases of copper with five infusion was examined. After a series of experiments, for example, different weight Of sample, artificial mixed single phases sample and repeated experiment, the extracting ratio constants foil all phases in different infusions were obtained. The mathematical model which is on the basis of the stable extracting constant is established. Based on the purpose-control genetic algorithms, we only need to determine the total copper in the sample in different infusions under the selected condition, then the value of respective phase could be calculated. The-relative standard deviation are (%) free oxidized copper phase:1.5, conjunction oxidized copper phase: 17.6,, secondary copper sulfide phase: 1.9, primary copper sulfide phase : 2.7,total copper: 0.8.