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

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.

泥河湾盆地大长梁剖面岩石磁学和环境磁学研究

The fluvio-lacustrine sequence in the Nihewan Basin is an important archive of late Pliocene-Pleistocene climate and environment changes in temperate northern China, which provides excellent sources of early human settlements in high latitude East Asia. The recent years have witnessed a considerable progress in the paleomagnetic dating of its stratigraphy, which has notably increased our understanding of a series of important issues such as the early human occupation in the Old World, the infilling history of the Nihewan Basin, and the chronological sequence of the Nihewan faunas. Up to now, the long-term paleoenvironmental changes directly retrieved from this basin, which might influence the evolution and expansion of early humans in the Nihewan Basin, are still poorly constrained, although several paleoclimatic records have been retrieved from this area. In this study, a combined mineral-magnetic and geochemical investigation was carried out on the fluvio-lacustrine sequence from the Dachangliang section at the eastern margin of the basin in order to reveal its rock magnetic and environmental magnetic characteristics and its implications for early human evolution in East Asia. The major findings and conclusions are listed as the following: First, there is an increased cooling coupled with an intensified aridification recorded in the fluvio-lacustrine sequence of the Dachangliang section. The cooling is related to an up-section decrease in propensity to chemical weathering as inferred from an increase in low-field susceptibility after cycling to 700 °C. Close to 700 °C, reacting chlorite is providing the iron source for newly formed very fine-grained ferrimagnetic minerals which enhances the susceptibility signal. The reactivity of chlorite after annealing at temperatures above 600 °C is documented with X-ray diffraction. Second, degrees of chemical weathering in the Nihewan Basin are further estimated by clay mineralogy (i.e. chlorite and illite contents and chlorite/illite ratio) and a series of major element proxies (i.e. Na2O/Al2O3 versus K2O/Al2O3 diagram, Al2O3-(CaO + Na2O)-K2O ternary diagram (A-CN-K), chemical index of alteration (CIA), (CaO + Na2O + MgO)/TiO2, (CaO + Na2O + MgO + K2O)/(TiO2 + Al2O3), CaO/Al2O3 and CaO/TiO2). The up-section decrease in propensity to chemical weathering suggested by the aforementioned rock mangetic measurement is further confirmed by these geochemical analyses. Combining the chemical weathering records from the Nihewan Basin, Chinese Loess Plateau, South China Sea and eastern China, we find that the consecutive decreasing trend in chemical weathering intensity during the late Cenozoic is ubiquitous across China. This pattern may result from a long-term decreasing East Asian summer monsoon and increasing East Asian winter monsoon, and thus a consecutive increasing of aridification and cooling in Asia during the Quaternary. Furthermore, the chemical weathering intensity increased from South China to North China during the Quaternary, in line with the decreasing East Asian summer monsoon and increasing East Asian winter monsoon and thus the gradually intensified aridification and cooling from South China to North China. Third, a combined mineral-magnetic and geochemical investigation provides evidences that the large-amplitude alterations of concentration of magnetic minerals mainly result from preservation/dissolution cycles of detrital magnetic minerals in alternately oxic and anoxic depositional environments. The preservation/dissolution model implies that the high-magnetic and low-magnetic cycles of this sedimentary sequence represent glacial and interglacial climate cycles, respectively. This contribute significnatly to our understanding of the link between climate and magnetic properties. Finally, the paleoclimatic implications of these rock magnetic and geochemical characteristics significantly increase our understanding of the general setting of early humans in high northern latitude in East Asia. We propose that the cold and dry climate may have contributed significantly to the expansion and adaptation of early humans, rather than bringing hardship, as is often thought. The relationship between magnetic properties and climate possibly provides valuable information on the climatic context of the Paleolithic sites in the basin, especially whether the occupation occurred during an interglacial or glacial period.

中国北方风成堆积的风化作用与环境变迁—U-Th、~(10)Be及元素地球化学的研究

The surface of the Earth is continuously undergoing changes as a result of weathering-erosion, plate tectonics and volcanic processes. Continental weathering-erosion with its complex rock-water interactions is the central process of global biochemical cycling of elements, and affects the long-term ocean atmosphere budget of carbon dioxide both through the consumption of carbonic acid during silicate weathering and through changes in the weathering and burial rates of organic carbon. Rates of the weathering-erosion depend on a variety of factors, in particular rock properties and chemical composition, climate (especially rainfall), structure, and elevation. They are quite variable on a regional scale. Thus, environmental changes in a region could be indicated by the history of weathering-erosion in the region. Recent attention has focused on increased silicate weathering of tectonically uplifted areas in the India-Asia collision zone as a possible cause for falling atmospheric CO_2 levels in the Cenozoic era. The wind blown dust deposits in the Loess Plateau is derived from the arid and semiarid regions in northwestern China, in turn, where the deposits have been derived from the Qinghai-Xizang Plateau and the high mountains around. Therefore, geochemistry of the wind blown loess-paleosol and red clay sequences may provide insight both to paleoenvironmental changes on the Loess Plateau, and to the uplift and weathering-erosion histories of the Qinghai-Xizang Plateau. In this paper, uranium-thorium series nuclides and cosmogenic ~(10)Be have been employed as tracers of weathering intensities and histories of the dust sediments in the Loess Plateau. Major elements, such as Na, Al, Fe etc., are also used to estimate degree of chemical alteration of the dust sediments and to rebuild the history of weathering on the Loess Plateau. First of all, using a low-level HPGe γ-ray detector, we measured U and Th series nuclides in 170 loess and paleosol samples from five sites in the Loess Plateau, going back 2.6 Ma. The results show that ~(238)U activities are disequilibrium with its daughter nuclide ~(230)Th in young loess-paleosol sequence, indicating that weathering was happened both in dust deposition site and in dust source regions. Using concentrations of ~(238)U and ~(232)Th in the samples, we estimated the amounts of ~(238)U leached out of from paleosols due to weathering. Further, based on analyses of ~(230)Th in paleosols deposited in the past ca. 140 ka, we determined when the paleosols weathered in the source regions. We conclude that most of the weathering in the dust-source regions may have occurred during the interglacials before dust deposition.