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

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.

Molecular docking and 3D-QSAR study of pyranmycin derivatives against 16S rRNA A site

To understand pharmacophore properties of pyranmycin derivatives and to design novel inhibitors of 16S rRNA A site, comparative molecular field analysis (CoMFA) approach was applied to analyze three-dimensional quantitative structure-activity relationship (3D-QSAR) of 17 compounds. AutoDock 3.0.5 program was employed to locate the orientations and conformations of the inhibitors interacting with 16S rRNA A site. The interaction mode was demonstrated in the aspects of inhibitor conformation, hydrogen bonding and electrostatic interaction. Similar binding conformations of these inhibitors and good correlations between the calculated binding free energies and experimental biological activities suggest that the binding conformations of these inhibitors derived from docking procedure were reasonable. Robust and predictive 3D-QSAR model was obtained by CoMFA with q(2) values of 0.723 and 0.993 for cross-validated and noncross-validated, respectively. The 3D-QSAR model built here will provide clear guidelines for novel inhibitors design based on the Pyranmycin derivatives against 16S rRNA A site. (c) 2005 Elsevier B.V. All rights reserved.

Experimental and field studies of Mesocyclops notius (Copepoda : Cyclopoida) from a subtropical Chinese lake, Lake Donghu

This paper reports experimental and field studies on the cyclopoid Mesocyclops notius from subtropical Lake Donghu close to the Yangtze river. Mesocyclops notius, a dominant crustacean Zooplankter throughout tropical Australia, was previously considered to be endemic to Australia, but recently, Mesocyclops leuckarti in Lake Donghu was re-identified as M. notius. Laboratory culture experiments were conducted to reveal the effect of temperature (15, 20, 25 and 30 degreesC) on the development, growth and reproduction of M. notius. Temperature was inversely related to development times of eggs, nauplii and copepodites, body length and physiological longevity of adults, and brood size. Body length and physiological longevity of females were greater than those of males at the same temperature. No reproduction took place at 15 degreesC. Production and seasonal cycles of M. notius during 1980-1982 were studied at two sampling stations of Lake Donghu. At the mid-lake station, the annual production and production/biomass (P/B) ratio of M. notius varied between 6.0 and 18.0 g dry wt m(-2) year(-1) and 74.6 and 95.5, respectively. Mesocyclops notius reached their highest density peaks in the warm months (July-October), with a maximum density of 1256 individuals l(-1) at a littoral site. No reproduction and recruitment by AL notius took place during the cold months (December March) when the temperature of the lake water was < 15 degreesC. Mesocyclops notius were more abundant at a littoral station than at a pelagic station, possibly due to different food availability. The higher male:female sex ratio of M. notius at the littoral station was most likely caused by size-selective fish predation on larger females.

Classical two-site fidelity and quantum phase transitions in the Ising model and the extended Hubbard model

In this Letter, the classical two-site-ground-state fidelity (CTGF) is exploited to identify quantum phase transitions (QPTs) for the transverse field Ising model (TFIM) and the one-dimensional extended Hubbard model (EHM). Our results show that the CTGF exhibits an abrupt change around the regions of criticality and can be used to identify QPTs in spin and fermionic systems. The method is especially convenient when it is connected with the density-matrix renormalization group (DMRG) algorithm. (C) 2008 Elsevier B.V. All rights reserved.

Simulated morphological landscape of polymer single crystals by phase-field model

The novel phase field model with the "polymer characteristic" was established based on a nonconserved spatiotemporal Ginzburg-Landau equation (TDGL model A). Especially, we relate the diffusion equation with the crystal growth faces of polymer single crystals. Namely, the diffusion equations are discretized according to the diffusion coefficient of every lattice site in various crystal growth faces and the shape of lattice is selected based on the real proportion of the unit cell dimensions.

An investigation of low-field magnetoresistance in the double perovskites Sr2Fe1-xZnxMoO6, x=0, 0.05, 0.15 and 0.25

The electrical, magnetic and transport properties of Zn doped polycrystalline samples of Sr2Fe1-xZnxMoO6 ( x = 0, 0.05, 0.15 and 0.25) with the double perovskite structure have been investigated. The subtle replacement of Fe3+ ions by Zn2+ ions facilitates the formation of a more ordered structure, while further substitution leads to disordered structure because of the presence of a striped phase. Analysis of the x-ray powder diffraction patterns based on Rietveld analysis indicates that the replacement of Fe3+ by Zn2+ ions favours the formation of Mo6+ ions. The spin-glass behaviour can be explained on the basis of the competition between the antiferromagnetic superexchange and the ferromagnetic double-exchange interaction. The low-field magnetoresistance was moderately enhanced at x = 0.05, and its origin was found to be the competition between the decrease of the concentration of the itinerant electrons and the weaker antiferromagnetic superexchange in the antiphase boundaries. An almost linear negative magnetoresistance in moderate field has been observed for x = 0.25. A possible double-exchange mechanism is proposed for elucidating the observations; it also suggests a coexistence of (Fe3+, Mo5+) and (Zn2+, Mo6+) valence pairs.

Infrared heater arrays for warming ecosystem field plots

There is a need for methodology to warm open-field plots in order to study the likely effects of global warming on ecosystems in the future. Herein, we describe the development of arrays of more powerful and efficient infrared heaters with ceramic heating elements. By tilting the heaters at 45 degrees from horizontal and combining six of them in a hexagonal array, good uniformity of warming was achieved across 3-m-diameter plots. Moreover, there do not appear to be obstacles (other than financial) to scaling to larger plots. The efficiency [eta(h) (%); thermal radiation out per electrical energy in] of these heaters was higher than that of the heaters used in most previous infrared heater experiments and can be described by: eta(h) = 10 + 25exp(-0.17 u), where u is wind speed at 2 m height (m s(-1)). Graphs are presented to estimate operating costs from degrees of warming, two types of plant canopy, and site windiness. Four such arrays were deployed over plots of grass at Haibei, Qinghai, China and another at Cheyenne, Wyoming, USA, along with corresponding reference plots with dummy heaters. Proportional integral derivative systems with infrared thermometers to sense canopy temperatures of the heated and reference plots were used to control the heater outputs. Over month-long periods at both sites, about 75% of canopy temperature observations were within 0.5 degrees C of the set-point temperature differences between heated and reference plots. Electrical power consumption per 3-m-diameter plot averaged 58 and 80 kW h day(-1) for Haibei and Cheyenne, respectively. However, the desired temperature differences were set lower at Haibei (1.2 degrees C daytime, 1.7 degrees C night) than Cheyenne (1.5 degrees C daytime, 3.0 degrees C night), and Cheyenne is a windier site. Thus, we conclude that these hexagonal arrays of ceramic infrared heaters can be a successful temperature free-air-controlled enhancement (T-FACE) system for warming ecosystem field plots.