Composite analytical solutions for a soil vapour extraction system

Soil vapour extraction (SVE) is a common remediation technique for cleaning up unsaturated soils contaminated by volatile organic compounds (VOCs). Analytical solutions, which result from simple mathematical models, can allow the fast approximation of the time‐dependent effluent concentration and the gaining of insight into the processes that take place during soil remediation. Deriving the analytical solutions to advection–dispersion equations that simultaneously take into account the mechanical dispersion and molecular diffusion is very difficult because of the variable dependence of governing equations' coefficients. In this study, we first present two simplified analytical solutions that only consider mechanical dispersion or molecular diffusion. The two developed analytical solutions are compared with the numerical solution that simultaneously considers both mechanical dispersion and molecular diffusion to examine the applicability of the two simplified analytical solutions and distinguishes the individual contribution of the mechanical dispersion and molecular diffusion to total VOCs transport in an SVE system. Results show that dispersion plays an important role during SVE decontamination and neither the diffusion‐dominated solution nor the dispersion‐dominated solution can agree well with the numerical solution when both mechanical dispersion and molecular diffusion have significant contributions to the total VOCs transport flux. A composite analytical solution that linearly couples the diffusion‐ and dispersion‐dominated analytical solutions, which is proposed herein to eliminate the discrepancy between the analytical solutions and the numerical solution. Results indicate that the proposed composite analytical solution agrees well with the numerical solution and is an effective tool for quickly and accurately evaluating the time‐dependent effluent concentration for parameters of the different ranges of interest in an SVE remedial system. Copyright © 2006 John Wiley & Sons, Ltd.     

Identifying parameters controlling soil delayed behaviour from laboratory and in situ pressuremeter testing

The aim of this paper is to present a methodology for identifying the soil parameters controlling the delayed behaviour from laboratory and in situ pressuremeter tests by using an elasto‐viscoplastic model (EVP‐MCC) based on Perzyna's overstress theory and on the elasto‐plastic Modified Cam Clay model. The influence of both the model parameters and the soil permeability was studied under the loading condition of pressuremeter tests by coupling the proposed model equations with Biot's consolidation theory. On the basis of the parametric study, a methodology for identifying model parameters and soil permeability by inverse analysis from three levels of constant strain rate pressuremeter tests was then proposed and applied on tests performed on natural Saint‐Herblain clay. The methodology was validated by comparing the optimized values of soil parameters and the values of the same parameters obtained from laboratory test results, and also by using the identified parameters to simulate other tests on the same samples. The analysis of the drainage condition and the strain rate effect during a pressuremeter test demonstrated the coupled influence of consolidation and viscous effects on the test results. The numerical results also showed that the inverse analysis procedure could successfully determine the parameters controlling the time‐dependent soil behaviour. Copyright © 2008 John Wiley & Sons, Ltd.     

Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.     

Accelerating the convergence of coupled geomechanical‐reservoir simulations

The pressure variations during the production of petroleum reservoir induce stress changes in and around the reservoir. Such changes of the stress state can induce marked deformation of geological structures for stress sensitive reservoirs as chalk or unconsolidated sand reservoirs. The compaction of those reservoirs during depletion affects the pressure field and so the reservoir productivity. Therefore, the evaluation of the geomechanical effects requires to solve in a coupling way the geomechanical problem and the reservoir multiphase fluid flow problem. In this paper, we formulate the coupled geomechanical‐reservoir problem as a non‐linear fixed point problem and improve the resolution of the coupling problem by comparing in terms of robustness and convergence different algorithms. We study two accelerated algorithms which are much more robust and faster than the conventional staggered algorithm and we conclude that they should be used for the iterative resolution of coupled reservoir‐geomechanical problem. Copyright © 2006 John Wiley & Sons, Ltd.     

Water source characterization through spatiotemporal patterns of major,minor and trace element stream concentrations in a complex,mesoscale German catchment

The link between spatiotemporal patterns of stream water chemistry and catchment characteristics for the mesoscale Dill catchment (692 km²) in Germany is explored to assess the catchment scale controls on water quality and to characterize water sources. In order to record the spatiotemporal pattern, ‘snapshot sampling’ was applied during low, mean and high flow, including 73 nested sites throughout the catchment. Water samples were analysed for the elements Li, B, Na, Mg, Al, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Rb, Sr, Mo, Ba, Pb and U using inductively‐coupled‐plasma mass spectrometry, and for electric conductivity and pH. Principle component analysis and hierarchical cluster analysis were used to find typical element associations and to group water samples according to their hydrochemical fingerprints. This revealed regional hydrochemical patterns of water quality which were subsequently related to catchment attributes to draw conclusions about the controls on stream chemistry. It was found that various lithologic signals and anthropogenic point source inputs controlled the base flow hydrochemistry. During increased flows, stream waters were diluted causing additional hydrochemical variability in response to heterogeneous precipitation inputs and differences in aquifer storage capacities. The hydrochemical patterns further displayed in‐stream mixing of waters. This implied, that stream waters could be apportioned to the identified water sources throughout the catchment. The basin‐wide hydrochemical variability has the potential to outrange the tracer signatures typically inferred in studies at the hillslope scale and is able to strongly influence the complexity of the catchment output. Both have to be considered for further catchment scale tracer and modelling work. Despite the likelihood of non‐conservative behaviour, the minor and trace elements enhanced the rather qualitative discrimination of the various groundwater types, as the major cations were strongly masked by point source inputs. Copyright © 2007 John Wiley & Sons, Ltd.     

Storm‐event rainfall–runoff modelling approach for ungauged sites in Taiwan

This work develops a top‐down modelling approach for storm‐event rainfall–runoff model calibration at unmeasured sites in Taiwan. Twenty‐six storm events occurring in seven sub‐catchments in the Kao‐Ping River provided the analytical data set. Regional formulas for three important features of a streamflow hydrograph, i.e. time to peak, peak flow, and total runoff volume, were developed via the characteristics of storm event and catchment using multivariate regression analysis. Validation of the regional formulas demonstrates that they reasonably predict the three features of a streamflow hydrograph at ungauged sites. All of the sub‐catchments in the study area were then adopted as ungauged areas, and the three streamflow hydrograph features were calculated by the regional formulas and substituted into the fuzzy multi‐objective function for rainfall–runoff model calibration. Calibration results show that the proposed approach can effectively simulate the streamflow hydrographs at the ungauged sites. The simulated hydrographs more closely resemble observed hydrographs than hydrographs synthesized using the Soil Conservation Service (SCS) dimensionless unit hydrograph method, a conventional method for hydrograph estimation at ungauged sites in Taiwan. Copyright © 2008 John Wiley & Sons, Ltd.     

Lie group variational integrators for the full body problem in orbital mechanics

Equations of motion, referred to as full body models, are developed to describe the dynamics of rigid bodies acting under their mutual gravitational potential. Continuous equations of motion and discrete equations of motion are derived using Hamilton’s principle. These equations are expressed in an inertial frame and in relative coordinates. The discrete equations of motion, referred to as a Lie group variational integrator, provide a geometrically exact and numerically efficient computational method for simulating full body dynamics in orbital mechanics; they are symplectic and momentum preserving, and they exhibit good energy behavior for exponentially long time periods. They are also efficient in only requiring a single evaluation of the gravity forces and moments per time step. The Lie group variational integrator also preserves the group structure without the use of local charts, reprojection, or constraints. Computational results are given for the dynamics of two rigid dumbbell bodies acting under their mutual gravity; these computational results demonstrate the superiority of the Lie group variational integrator compared with integrators that are not symplectic or do not preserve the Lie group structure.     

Coupled surface water–groundwater model and its application in the arid Shiyang River basin,China

The planning and management of water resources in the Shiyang River basin, China require a tool for assessing the impact of groundwater and stream use on water supply reliabilities and improving many environment‐related problems such as soil desertification induced by recent water‐related human activities. A coupled model, integrating rule‐based lumped surface water model and distributed three‐dimensional groundwater flow model, has been established to investigate surface water and groundwater management scenarios that may be designed to restore the deteriorated ecological environment of the downstream portion of the Shiyang River basin. More than 66% of the water level among 24 observation wells have simulation error less than 1·0 m. The overall trend of the temporal changes of simulated and observed surface runoff at the Caiqi gauging station remains almost the same. The calibration was considered satisfactory. Initial frameworks for water allocation, including agricultural water‐saving projects, water diversion within the basin and inter‐basin water transfer, reducing agricultural irrigation area and surface water use instead of groundwater exploitation at the downstream were figured out that would provide a rational use of water resources throughout the whole basin. Sixteen scenarios were modelled to find out the most appropriate management strategies. The results showed that in the two selected management options, the groundwater budget at the Minqin basin was about 1·4 × 10⁸ m³/a and the ecological environment would be improved significantly, but the deficit existed at the Wuwei basin and the number was about 0·8 × 10⁸ m³/a. Water demand for domestic, industry and urban green area would be met in the next 30 years, but the water shortage for meeting the demand of agricultural water use in the Shiyang River basin was about 2·2 × 10⁸ m³/a. It is suggested that more inter‐basin water transfer should be required to obtain sustainable water resource use in the Shiyang River basin. Copyright © 2009 John Wiley & Sons, Ltd.     

柴达木盆地侏罗系油气成藏模式与油气聚集规律

针对柴达木盆地侏罗系油气勘探目标选择的难点,从烃源岩与圈闭时空配置的角度,建立侏罗系油气早期聚集早期成藏、早期聚集晚期成藏、晚期聚集晚期成藏3种模式。盆地油气来源、充注时期、油藏调整及破坏等方面存在规律性,这决定了油气分布的有序性,呈现出“断褶区聚气,山前、潜伏区聚油,凸起顶、凸起斜坡区油气并存”的基本格局和油气富集程度依次变好的序列。据此预测油气富集范围,从中优选冷湖4号东斜坡、乌南—巴依凸起斜坡、怀头他拉—德北斜坡、潜伏Ⅰ号等新靶区。