A Two-Dimensional Areal Model for Density Dependent Flow Regime

A two-dimensional (2D) plane model of saltwater intrusion was developed, for the simulation of groundwater level and the average solute concentration in a 2D horizontal plane, together with the estimation of the saltwater depth. The proposed approach is of particular interest when assessing the effect of different regional pumping scenarios on groundwater level and its quality. The corresponding MEL2DSLT code was developed on the basis of the Modified Eulerian–Lagrangian (MEL) method to overcome difficulties arising from hyperbolic behavior of flow and transport equations, due to the advective nature of solute transport and heterogeneity of the soil characteristics (permeabilities and dispersivities). The code was verified against the 2D cross sectional model SUTRA and the three-dimensional (3D) model SWICHA. Simulation was conducted concerning the problem of saltwater intrusion in the Khan Yunis portion of the phreatic coastal aquifer of Gaza Strip. After calibrating the model for the aquifer parameters, we investigated its predictions resulting from various regional pumping scenarios using the actual pumping intensity from the year 1985 and extrapolating on the basis of 3.8% annual population growth. Results show a considerable depletion of groundwater level and intrusion of seawater due to excessive pumping.     

Salt Water Intrusion in a Three-dimensional Groundwater System in The Netherlands: A Numerical Study

Salt water intrusion is investigated in a coastal groundwater system in the northern part of the province Noord-Holland, The Netherlands. Density dependent groundwater flow is modeled in three-dimensions with MOCDENS3D. This computer code is a version of MOC3D (Konikow et al., 1996) that has been adapted to simulate transient density-driven groundwater flow. Results from the model suggests that in this Dutch hydrogeologic system a severe and irreversible salinisation is already occurring. Within a few tens to hundreds of years, the salinity of the shallow aquifer is estimated to increase substantially. This salinisation process is a result of human activities such as the reclamation of the low-lying areas during the past centuries. Without changing the present boundary conditions, seepage into the low-lying areas will decrease slightly because of predicted increases in groundwater salinity. However, the rate in salt load through the Holocene aquitard into the low-lying areas will increase significantly due to an increase in salinity in the shallow aquifer. In addition, a relative sea level rise of 0.5m per century will intensify the salinisation process, causing an enormous increase in salt load in all low-lying areas in this part of The Netherlands.     

Simulation of flow around rigid vegetation stems with a fast method of high accuracy

This paper develops the virtual cylinder model (VCM) to simulate flows through vegetation canopies at low to medium high Reynolds numbers in wetlands. Using a combined Eulerian and Lagrangian approach, the VCM is capable of simulating the flow field around each vegetation stem (modeled as a cylinder) as well as a bulk of numerous stems in a fluid. Numerical results of flows through regularly and randomly arranged cylinders have been compared with direct numerical simulations. Good agreement has been achieved. This model maintains an excellent balance between accuracy and efficiency in modeling interactions between the fluid and vegetation stems. Simulation results demonstrate that this model is promising for investigating the hydrodynamics and vegetation resilience in wetlands that serve as a first line of defense against the sea level rise, storm impacts and coastal erosion.     

The effects of mean sea level rise and strengthened winds on extreme sea levels in the Baltic Sea

Mean sea level rise and climatological wind speed changes occur as part of the ongoing climate change and future projections of both variables are still highly uncertain. Here the Baltic Sea's response in extreme sea levels to perturbations in mean sea level and wind speeds is investigated in a series of simulations with a newly developed storm surge model based on the nucleus for European modeling of the ocean(NEMO)-Nordic. A simple linear model with only two tunable parameters is found to capture the changes in the return levels extremely well. The response to mean sea level rise is linear and nearly spatially uniform, meaning that a mean sea level rise of 1 m increases the return levels by a equal amount everywhere. The response to wind speed perturbations is more complicated and return levels are found to increase more where they are already high. This behaviour is alarming as it suggests that already flooding prone regions like the Gulf of Finland will be disproportionally adversely affected in a future windier climate.     

地下水位上升下黄土斜坡稳定性分析

黄土高原一些地区,由于塬上引水灌溉使得地下水位不断抬升,造成黄土滑坡频繁发生。地下水位变化严重影响着黄土斜坡的稳定性。基于饱和-非饱和渗流理论和延伸的摩尔-库伦破坏准则,结合室内饱和和非饱和试验结果,针对泾阳南塬一典型黄土斜坡,考虑地下水位上升情况下,对其进行了瞬态饱和-非饱和渗流分析;然后将计算得到的瞬态孔隙水压力分布用于斜坡的极限平衡分析。结果表明:地下水位上升对暂态渗流场和斜坡稳定性有明显影响;考虑非饱和渗流和吸力强度的边坡稳定分析方法更加符合实际情况。     

黄河黑山峡大柳树坝坝址岩体渗漏量及灌浆量数值模拟研究

大坝坝址岩体的渗漏研究对于大坝的安全运营至关重要。在进行渗漏治理的措施中,灌浆是常见的一种,但目前确定灌浆的各项参数,如配合比,孔距等,需通过灌浆试验花费大量的人财物力及时间才能确定。为快速查清每一种灌浆方案所灌浆量,本文针对模拟坝址水库渗漏灌浆效果的三维数值模拟开展了计算研究。在对大柳树坝坝址岩体地质情况及渗漏情况分析的基础上,将蓄水前的地下水运移规律模拟出,与实际通过钻孔平硐资料及分析得出的实际地下水运移规律做比较,验证对地下水运移规律的假设; 然后将蓄水后的地下水运移规律模拟出,确定地下水的运移规律,了解地下水运移的方向,验证与根据实际情况初步判断的地下水运移方向是否吻合; 最后通过数值模拟确定在一定灌浆浆液配合比的情况下所需灌浆量与实际的灌浆试验结果做比较。得出结果与实际灌浆试验结果吻合,证明通过数值模拟,建立浆液配合比与灌浆效果之间的关系,在类似情况快速地选取配合比是可行的。     

Modeling Biogeochemical Processes in Leachate-Contaminated Soils: A Review

During subsurface transport, reactive solutes are subject to a variety of hydrological, physical and biochemical processes. The major hydrological and physical processes include advection, diffusion and hydrodynamic dispersion, and key biochemical processes are aqueous complexation, precipitation/dissolution, adsorption/desorption, microbial reactions, and redox transformations. The addition of strongly reduced landfill leachate to an aquifer may lead to the development of different redox environments depending on factors such as the redox capacities and reactivities of the reduced and oxidised compounds in the leachate and the aquifer. The prevailing redox environment is key to understanding the fate of pollutants in the aquifer. The local hydrogeologic conditions such as hydraulic conductivity, ion exchange capacity, and buffering capacity of the soil are also important in assessing the potential for groundwater pollution. Attenuating processes such as bacterial growth and metal precipitation, which alter soil characteristics, must be considered to correctly assess environmental impact. A multicomponent reactive solute transport model coupled to kinetic biodegradation and precipitation/dissolution model, and geochemical equilibrium model can be used to assess the impact of contaminants leaking from landfills on groundwater quality. The fluid flow model can also be coupled to the transport model to simulate the clogging of soils using a relationship between permeability and change in soil porosity. This paper discusses the different biogeochemical processes occurring in leachate-contaminated soils and the modeling of the transport and fate of organic and inorganic contaminants under such conditions.     

Effects of Clay Dispersion on Aquifer Storage and Recovery in Coastal Aquifers

Cyclic injection, storage, and withdrawal of freshwater in brackish aquifers is a form of aquifer storage and recovery (ASR) that can beneficially supplement water supplies in coastal areas. A 1970s field experiment in Norfolk, Virginia, showed that clay dispersion in the unconsolidated sedimentary aquifer occurred because of cation exchange on clay minerals as freshwater displaced brackish formation water. Migration of interstitial clay particles clogged pores, reduced permeability, and decreased recovery efficiency, but a calcium preflush was found to reduce clay dispersion and lead to a higher recovery efficiency. Column experiments were performed in this study to quantify the relations between permeability changes and clay mineralogy, clay content, and initial water salinity. The results of these experiments indicate that dispersion of montmorillonite clay is a primary contributor to formation damage. The reduction in permeability by clay dispersion may be expressed as a linear function of chloride content. Incorporating these simple functions into a radial, cross-sectional, variable-density, ground-water flow and transport model yielded a satisfactory simulation of the Norfolk field test – and represented an improvement over the model that ignored changes in permeability. This type of model offers a useful planning and design tool for ASR operations in coastal clastic aquifer systems.     

Lattice model applied to the fracture of large strain composite

An improved lattice model is developed to simulate the fracture behavior of large strain composite. Based on equivalent relation between the continuum and the lattice model for small deformation, the equivalent relation between large strain continuum and improved lattice model is established by introducing large strain elastic law into the lattice system. The theory can simulate large deformation. The program of large strain lattice model simulates several representative problem of large strain elasticity. The results agree with the theoretical results. Assumed failure criterion is used to describe the fracture process of large strain elasticity and large strain composite. The improved lattice model provides an effective method for fracture simulation of large strain composite.     

海上溢油迁移转化的双层数学模型

海上溢油的双层数学模型将泄漏在海域中的溢油考虑为表层油膜和分布在整个水体中的悬浮 油滴两层组成. 采用Lagrangian追踪法模拟油膜的输移扩散. 在此方法中,油膜可用大量的小油滴表示,对每个油滴都规定随时间而变化的坐标系,油滴 的运动受到风、潮流和周围油的浓度影响. 油的迁移过程包括对流、扩展、湍动扩 散、附着在岸边以及沉降到海底等过程. 转化过程包括挥发、溶解、乳化等. 此外,光化学反应、生物降解能够改变油的特征以及减小油的污染. 该模型可以用于瞬时和连续溢油的情况,不仅可以用于模拟油,也可以模拟其它与油的密度 相近的有害物质的泄漏. 水动力学模型采用美国普林斯顿海洋模式(POM). 该溢油模型已应用于渤海海峡突发性溢油事故的模拟与预报.     

海岸波浪场模型研究进展

从建模原理、波浪在近岸区域传播的众多机制、模型的类别、优势、局限性以及模型在未来的发展趋势等方面,综述了在海岸工程实践中广泛运用的以下两大类海岸波浪场预测模型的最新研究进展:(1)能量平衡模型.它一般用来预测海洋深水波候,已发展到相当完善的阶段,例如,最为著名的WAM3G模型.这种模型在海岸工程中的作用就在于可以模拟施加在波浪上的随时间变化的风场效应.(2)质量、动量守恒模型.它在海岸工程中应用最为普遍,并且内容丰富,数值技巧多样.目前包含了以下代表性的模型:缓坡方程、抛物型方程、非线性浅水方程、高阶Boussinesq型方程、Green-Naghdi理论.      

露天矿边坡饱和渗透系数随机场数值模型研究

为分析强降雨入渗及渗透系数空间变异性对闭坑露天矿边坡渗流场的影响程度,基于非饱和渗流理论和随机场理论,采用非侵入式随机方法,通过FISH语言编写非饱和区单元饱和度、渗透系数与基质吸力的修正函数,建立饱和-非饱和渗流随机场模型,开展强降雨作用的高大陡深岩质边坡渗流特征研究。研究结果表明,修正饱和-非饱和渗流随机场模型能够准确地描述露天矿边坡的降雨入渗过程。且降雨入渗主要影响到露天矿边坡浅层渗流场。随降雨持时变化,在坡面与地下水位线之间形成包围的且逐渐缩小的非饱和区。坡面最早出现暂态饱和区且厚度逐渐增加,但增幅逐渐放缓。同时坡面点孔隙水压力最早达到稳定值零,离坡面越远的点孔隙水压力达到稳定值零的用时会越长。该结论可为闭坑露天矿边坡的地质灾害风险预测提供参考意义。     

A MULTI-LEVEL NUMERICAL MODEL OF COASTAL UPWELLING: A DIAGNOSTIC STUDY

A two-dimensional baroclinic model is described for coastal upwelling in a vertical plane perpendicular to the coast. The model consists of equations of motion, continuity and turbulence energy along with equations for salinity and thermal energy and an equation of state. The role of density gradient in the baroclinic pressure gradient is investigated to understand the dynamics during the upwelling process. To represent the surface and bottom boundaries corresponding to a fixed computational level in the discretized equations, a set of non-dimensional co-ordinates is used. These co-ordinates are then transformed onto logarithmic co-ordinate axes to resolve effectively the boundary layers. The first experiment is carried out with a flat bottom to understand the dynamics of the upwelling and the structural features of the process by diagnostic analysis of the balance between various terms of the momentum equation. Starting from a state of rest, a spatially uniform alongshore wind stress corresponding to the mean monthly wind stress for the month of May is applied and held constant thereafter. The fluid is assumed to be incompressible and stratified, with the initial temperature and salinity having no horizontal variations but a uniform vertical gradient. As the upwelling phenomenon is transient in nature and keeping in mind the additional computational overheads, the response of the model is studied day-wise up to 4 days. In the second experiment the model is applied to study the upwelling off the east coast of India in a plane normal to the coast of Visakhapatnam. The analysis area extends to 100 km offshore with real topography. The results are presented day-wise for 4 days, comparing the balance between various terms in the upwelling region and in the open sea, and the dynamics of the baroclinic coastal jet is explained. © 1997 John Wiley & Sons, Ltd.     

河口混合与泥沙输运

根据振荡边界层理论和波流分解方法,导出了河口往复水流的流速垂向分布廓线,据此建立了河口垂向准二维水流、盐度、泥沙运动模型。对泥沙输运,完整地考虑了其对流、扩散、起动和沉降的动力学过程。模拟结果与实测资料进行了对比。应用该模型研究河口泥沙输运,分析了河口混合对泥沙输运的影响及最大浑浊带的时空变化规律。     

Investigation of Groundwater Contaminant Discharge into Tidally influenced Surface-water Bodies: Theoretical Analysis

Data from an one-dimensional homogeneous sand column, which is utilized to investigate the effect of tides on the concentration of groundwater contaminants discharging to a surface-water body, demonstrate that the tidal fluctuations in water level elevation create concentration oscillations upgradient of the groundwater discharge locations and there is a resulting decrease in average contaminant concentration at the point of groundwater discharge to a surface-water body. The further upgradient an observation point is located, the smaller the amplitude of the tidally induced concentration oscillations. In addition, an excessive upstream migration of concentration oscillations is observed although there is a net downgradient flow. As the classical groundwater flow and transport model could not reproduce this phenomena, a multi-mobility model is proposed with one highly mobile liquid phase, one less mobile liquid phase and a solid phase. Averaging theory is applied in a first step to develop the macroscopic mass conservation equation from its microscale counterpart and then, in a second step, averaging is again used to reduce dimensionality to one-dimensional governing equations defined along the axis of the column. The simulation confirms the existence of an enhanced tidally induced mixing process and the suitability of our mathematical-physical representation of it.     

Finite volume solution to integrated shallow surface–saturated groundwater flow

This paper describes development of an integrated shallow surface and saturated groundwater model (GSHAW5). The surface flow motion is described by the 2‐D shallow water equations and groundwater movement is described by the 2‐D groundwater equations. The numerical solution of these equations is based on the finite volume method where the surface water fluxes are estimated using the Roe shock‐capturing scheme, and the groundwater fluxes are computed by application of Darcy's law. Use of a shock‐capturing scheme ensures ability to simulate steady and unsteady, continuous and discontinuous, subcritical and supercritical surface water flow conditions. Ground and surface water interaction is achieved by the introduction of source‐sink terms into the continuity equations. Two solutions are tightly coupled in a single code. The numerical solutions and coupling algorithms are explained. The model has been applied to 1‐D and 2‐D test scenarios. The results have shown that the model can produce very accurate results and can be used for simulation of situations involving interaction between shallow surface and saturated groundwater flows. Copyright © 2005 John Wiley & Sons, Ltd.     

沥青海砂作为软岩相似材料的实验研究

为了解决岩石工程相似模型实验中传统相似材料难以模拟深部软岩低强度、大变形和强流变的问题,本文选用70#沥青(产量大、成本低、粘弹性、易成型、易软化等)和海砂作为深部软岩的相似模拟材料.采用连续加压法对其进行了单轴压缩实验,对不同配比下试件的密度、抗压强度及初始弹性模量进行分析,得到了该相似材料有关的原始数据和规律.本文初步验证了使用该相似材料模拟深部软岩的可行性,并针对不足之处给出一些合理改进建议,为进一步模型模拟实验提供指导.     

A Decomposition Method for Solving Coupled Multi–Species Reactive Transport Problems

Concerns over the problems associated with mixed waste groundwater contamination have created a need for more complex models that can represent reactive contaminant fate and transport in the subsurface. In the literature, partial differential equations describing solute transport in porous media are solved either for a single reactive species in one, two or three dimensions, or for a limited number of reactive species in one dimension. Those solutions are constrained by many simplifying assumptions. Often, it is desirable to simulate transport in two or three dimensions for a more practical system that might have multiple reactive species. This paper presents a decomposition method to solve the partial differential equations of multi–dimensional, multi–species transport problems that are coupled by linear reactions. A matrix method is suggested as a tool for describing the reaction network. In this way, the level of complexity required to solve the multi–species reactive transport problem is significantly reduced.     

Conditional model for sorption in porous media with fractal properties

In this study, the conditional moment closure approach, which is proven to be very useful for modelling of reactions in turbulent flows, is extended to characterise adsorbing, desorbing or reacting flows in porous media. A complete specification of the porous distance conditioned moment closure model, which is formulated in terms of single-conditioned expectations, is presented. The closure of the model equations is obtained assuming the diffusion approximation for fluxes of the reactive species. The model simulates complex multi-cascade processes of convective and diffusive transport of species between pores in a continuous and consistent manner and is a generalisation of dual (or triple) porosity concept. The model addresses the major difficulty of describing transport, entrapment and sorption processes in porous media with fractal properties, where distant transport occurs in the largest pores or fractures, while the adsorbing or desorbing surface is mainly allocated in small pores. The model is able to simulate various regimes of methane replacement by CO₂ in a coal sample, which makes it useful for optimising the design and parameters of enhanced coal bed methane recovery operations. It is demonstrated that the power-low decrease in downstream methane concentration, which has been observed experimentally, can be accurately reproduced by the model.     

Multiscale Model for Assessing Effect of Bacterial Growth on Intrinsic Permeability of Soil: Model Description

A multiscale model is proposed for assessing the impact of bacterial growth attached to soil particles on intrinsic permeability of soil. The model combines two modules. One of the modules solves the equations describing multiphase reaction–diffusion–transport model at the macroscopic [Representative Elementary Volume (REV)] scale, augmented with an equation for attached biomass growth. This module is used as an upscaling environment. The other module in the coupling set is a cellular automaton (CA) model for bacterial growth attached to soil particles. The coupling is achieved by having the CA module model bacterial growth at the pore level and integrating the result to provide information at the REV level. The details of the two modules, the links between them and an algorithm for simulation using the coupled model are described. The two modules are designed to work alternately till the end of a simulation period. The continuum module (reaction–diffusion–transport augmented with attached biomass growth equation) calculates the change in attached biomass. This is followed by the determination of a change in the states of cells in the CA module. Any change in the states necessitates a recalculation of permeability values to be used for the next time step of the continuum module. A companion paper gives results on using this coupled model to simulate a column experiment.