盆地地下水资源开发的数值模拟研究－以山西灵石静升盆地为例

Groundwater basin is important for water supply in northern China. The paper took the Jingsheng Basin in Lingshi County, Shanxi Province as a case to study the basin groundwater system by numerical modeling. The hydrogeological characteristics were analysed basing on the field investigation, and a three-dimensional groundwater flow model was established to describe the groundwater flow system in the Jingsheng groundwater basin. The boundary of the model was determined by using geophysics and GIS data, and the lumped parameter model of runoff was used to depict the transform between the surface water and groundwater, and the groundwater dating data was used to calibrate the model. All these methods were used to improve the model. The Software Visual MODFLOW 2000 was applied to set up the numerical groundwater flow model. The groundwater flow pattern in the average year, the high-water year and the low-water year were simulated respectively by the model. Some new cognition to the groundwater movement in Jingsheng Basin was obtained in the paper. The difficult problems were resolved when using the conventional and theoretical analysis to forecast and appraise the exploitation of the groundwater, and supplies the instructional technology base for the reasonable exploitation and optimization collocation. The numerical model will improve evaluation of the basin groundwater resources.

填土边坡水文作用机理与破坏过程研究—以深圳羊宝地填土滑坡为例

A full understanding of failure mechanism, critical hydrological condition, and process of mobilization and deposition of a landslide is essential for optimal design of stabilization measure and forecasting of landslide hazard. This requires a quantitative study of hydrological response of a slope to rainfall through field monitoring, laboratory test and numerical modelling. At 13:40 on September 18, 2002, a fill slope failed following a period of prolonged rain in Shenzhen, resulting in 5 fatalities and 31 injuries. The failed mass with a volume about 2.5×104m3 traveled about 140m on level ground. Field monitoring, laboratory test, theoretical analysis and numerical modelling were carried out to undestand the hydrological response and failure mechanism of this fill slope. This thesis mainly focuses on the following aspects: (1) The hydrological responses and failure processes of slopes under rainfall infiltration were reviewed. Firstly, the factors influencing on the hydrological responses of slopes were analysed. Secondly, the change of stress state of slope soil and modelling methods of slope failure under rainfall infiltration were reviewed. (2) The characteristics of the Yangbaodi landslide and associated rainfall triggering the failure were presented. The failure was characterized by shallow flowslide, due to an increase of ground water table caused by rainfall infiltration. (3) A fully automated instrumentation was carried out to monitor rainfall, and saturated – unsaturated hydrological response of the fill slope, using a raingauge, piezometers, tensiometers and moisture probes. A conceptual hydrogeological model was presented based on field monitoring and borehole data. Analysis of monitoring data showed that the high pore water pressure in fill slope was caused by upward flow of semiconfined groundwater in the moderately decomposed granite. (4) Laboratory and in-situ testing was performed to study the physical and mechanical properties of fills. Isotropically consolidated undrained compression tests and anisotropically consolidated constant shear stress tests were carried out to understand the failure mechanism of the fill slope. It is indicated that loosely compacted soil is of strain-softening behaviour under undrained conditions, accompanied with a rapid increase in excess pore water pressure. In anisotropically consolidated constant shear stress tests, a very small axial strain was required to induce the failure and the excess pore water pressure increased quickly at failure. This indicated that static liquefaction caused by rise in groundwater table due to rainfall infiltration occurred. (5) The hydraulic conductivity of the highly and moderately decomposed granite was estimated using monitering data of pore water pressure. A saturated – unsaturated flow was modeled to study the hydrological response of the fill slope using rainfall records. It was observed that the lagged failure was due to the geological conditions and the discrepancy of hydraulic conductivity of slope soils. The hydraulic conductivity of moderately decomposed granite is relatively higher than the other materials, resulting in a semiconfied groundwater flow in the moderately decomposed granite, and subsequent upward flow into the upper fill layer. When the ground water table in the fill layer was increased to the critical state, the fill slope failed. (6) Numerical exercises were conducted to replay the failure process of the fill slope, based on field monitoring, laboratory and in-situ testing. It was found that the fill slope was mobilized by a rapid transfer of the concentrated shear stress. The movement of failure mass was characterized by viscosity fluid with a gradual increase in velocity. The failure process, including mobilization and subsequent movement and deposition, was studied using numerical methods.

地下水流动数值模拟的并行计算方法研究

Numerical modeling of groundwater is very important for understanding groundwater flow and solving hydrogeological problem. Today, groundwater studies require massive model cells and high calculation accuracy, which are beyond single-CPU computer’s capabilities. With the development of high performance parallel computing technologies, application of parallel computing method on numerical modeling of groundwater flow becomes necessary and important. Using parallel computing can improve the ability to resolve various hydro-geological and environmental problems. In this study, parallel computing method on two main types of modern parallel computer architecture, shared memory parallel systems and distributed shared memory parallel systems, are discussed. OpenMP and MPI (PETSc) are both used to parallelize the most widely used groundwater simulator, MODFLOW. Two parallel solvers, P-PCG and P-MODFLOW, were developed for MODFLOW. The parallelized MODFLOW was used to simulate regional groundwater flow in Beishan, Gansu Province, which is a potential high-level radioactive waste geological disposal area in China. 1. The OpenMP programming paradigm was used to parallelize the PCG (preconditioned conjugate-gradient method) solver, which is one of the main solver for MODFLOW. The parallel PCG solver, P-PCG, is verified using an 8-processor computer. Both the impact of compilers and different model domain sizes were considered in the numerical experiments. The largest test model has 1000 columns, 1000 rows and 1000 layers. Based on the timing results, execution times using the P-PCG solver are typically about 1.40 to 5.31 times faster than those using the serial one. In addition, the simulation results are the exact same as the original PCG solver, because the majority of serial codes were not changed. It is worth noting that this parallelizing approach reduces cost in terms of software maintenance because only a single source PCG solver code needs to be maintained in the MODFLOW source tree. 2. P-MODFLOW, a domain decomposition–based model implemented in a parallel computing environment is developed, which allows efficient simulation of a regional-scale groundwater flow. The basic approach partitions a large model domain into any number of sub-domains. Parallel processors are used to solve the model equations within each sub-domain. The use of domain decomposition method to achieve the MODFLOW program distributed shared memory parallel computing system will process the application of MODFLOW be extended to the fleet of the most popular systems, so that a large-scale simulation could take full advantage of hundreds or even thousands parallel processors. P-MODFLOW has a good parallel performance, with the maximum speedup of 18.32 (14 processors). Super linear speedups have been achieved in the parallel tests, indicating the efficiency and scalability of the code. Parallel program design, load balancing and full use of the PETSc were considered to achieve a highly efficient parallel program. 3. The characterization of regional ground water flow system is very important for high-level radioactive waste geological disposal. The Beishan area, located in northwestern Gansu Province, China, is selected as a potential site for disposal repository. The area includes about 80000 km2 and has complicated hydrogeological conditions, which greatly increase the computational effort of regional ground water flow models. In order to reduce computing time, parallel computing scheme was applied to regional ground water flow modeling. Models with over 10 million cells were used to simulate how the faults and different recharge conditions impact regional ground water flow pattern. The results of this study provide regional ground water flow information for the site characterization of the potential high-level radioactive waste disposal.

水-岩土反应的力学与环境效应研究

The research on mechanical effects of water-rock and soil interaction on deformation and failure of rocks and soils involves three aspects of mechanics, physics and chemistry. It is the cross between geochemistry and rock mechanics and soil mechanics. To sum up, the mechanical effects of water-rock and soil interaction is related to many complex processes. Research in this respect has been being an important forward field and has broad prospects. In connection with the mechanism of the effects of the chemical action of water-rock on deformation and failure of rocks and soils, the research significance, the present state, the developments in this research domain are summarized. Author prospects the future of this research. The research of the subject should be possessed of important position in studying engineering geology and will lead directly to a new understand on geological hazard and control research. In order to investigation the macroscopic mechanics effects of chemical kinetics of water-rock interaction on the deformation and failure, calcic rock, red sandstone and grey granite reacting chemically with different aqueous solution at atmospheric temperature and atmospheric pressure are uniaxially compressed. The quantitative results concerning the changes of uniaxially compressive strength and elastic modulus under different conditions are obtained. It is found that the mechanical effects of water on rock is closely related to the chemical action of water-rock or the chemical damage in rock, and the intensity of chemical damage is direct ratio to the intensity of chemical action in water-rock system. It is also found that the hydrochemical action on rock is time-dependent through the test. The mechanism of permeation and hydrochemical action resulting in failure of loaded rock mass or propagation of fractures in rocks would be a key question in rock fracture mechanics. In this paper, the fracture mechanical effects of chemical action of water-rock and their time- and chemical environment-dependent behavior in grey granite, green granite, grey sandstone and red sandstone are analyzed by testing K_(IC) and COD of rock under different conditions. It is found that: ①the fracture mechanical effect of chemical action of water-rock is outstanding and time-dependent, and high differences exist in the influence of different aqueous solution, different rocks, different immersion ways and different velocity of cycle flow on the fracture mechanical effects in rock. ②the mechanical effects of water-rock interaction on propagation of fractures is consistent with the mechanical effects on the peak strength of rock. ③the intensity of the mechanical fracture effects increases as the intensity of chemical action of water-rock increases. ④iron and calcium ion bearing mineral or cement in rock are some key ion or chemical composition, and especially iron ion-bearing mineral resulting in chemical action of water-rock to be provided with both positive and negative mechanical effects on rock. Through the above two tests, we suggest that primary factors influencing chemical damage in rock consist of the chemical property of rock and aqueous solution, the structure or homogeneity of rocks, the flow velocity of aqueous solution passing through rock, and cause of formation or evolution of rock. The paper explores the mechanism on the mechanical effects of water-rock interaction on rock by using the theory of chemistry and rock fracture mechanics with chemical damage proposed by author, the modeling method and the energy point of view. In this paper, the concept of absorbed suction between soil grains caused by capillary response is given and expounded, and the relation and basic distinction among this absorbed suction, surface tension and capillary pressure of the soil are analyzed and established. The law of absorbed suction change and the primary factors affecting it are approached. We hold that the structure suction are changeable along with the change of the saturation state in unsaturated soils. In view of this, the concept of intrinsic structure suction and variable structure suction are given and expounded, and this paper points out: What we should study is variable structure suction when studying the effective stress. By IIIy κHH's theory of structure strength of soils, the computer method for variable structure suction is analyzed, the measure method for variable structure suction is discussed, and it reach the conclusions: ①Besides saturation state, variable structure suction is affected by grain composition and packing patter of grains. ②The internal relations are present between structure parameter N in computing structure suction and structure parameter D in computing absorbed suction. We think that some problems exit in available principle of effective stress and shear strength theory for unsaturated soil. Based on the variable structure suction and absorbed suction, the classification of saturation in soil and a principle of narrow sense effective stress are proposed for unsaturated soils. Based on generalized suction, the generalized effective stress formula and a principle of generalized effective stress are proposed for unsaturated soils. The experience parameter χ in Bishop's effective stress formula is defined, and the principal factors influencing effective stress or χ. The primary factor affecting the effective stress in unsaturated soils, and the principle classifying unsaturated soils and its mechanics methods analyzing unsaturated soils are discussed, and this paper points out: The theory on studying unsaturated soil mechanics should adopt the micromechanics method, then raise it to macromechanics and to applying. Researching the mechanical effects of chemical action of water-soil on soil is of great importance to geoenvironmental hazard control. The texture of soil and the fabric of soil mass are set forth. The tests on physical and mechanical property are performed to investigate the mechanism of the positive and negative mechanical effects of different chemical property of aqueous solution. The test results make clear that the plastic limit, liquid limit and plasticity index are changed, and there exists both positive and negative effects on specimens in this test. Based on analyzing the mechanism of the mechanical effects of water-soil interaction on soil, author thinks that hydrochemical actions being provided with mechanical effects on soil comprise three kinds of dissolution, sedimentation or crystallization. The significance of these tests lie in which it is recognized for us that we may improve, adjust and control the quality of soils, and may achieve the goal geological hazard control and prevention.The present and the significance of the research on environmental effects of water-rock and soil interaction. Various living example on geoenvironmental hazard in this field are enumerated. Following above thinking, we have approached such ideals that: ①changing the intensity and distribution of source and sink in groundwater flow system can be used to control the water-rock and soil interaction. ②the chemical action of water-rock and soil can be used to ameliorate the physical and mechanical property of rocks and soils. Lastly, the research thinking and the research methods on mechanical effects and environmental effects of water-rock and soil interaction are put forward and detailed.

单（多）井抽灌对浅部地温场的影响研究

The technique of energy extraction using groundwater source heat pumps, as a sustainable way of low-grade thermal energy utilization, has widely been used since mid-1990's. Based on the basic theories of groundwater flow and heat transfer and by employing two analytic models, the relationship of the thermal breakthrough time for a production well with the effect factors involved is analyzed and the impact of heat transfer by means of conduction and convection, under different groundwater velocity conditions, on geo-temperature field is discussed.A mathematical model, coupling the equations for groundwater flow with those for heat transfer, was developed. The impact of energy mining using a single well system of supplying and returning water on geo-temperature field under different hydrogeological conditions, well structures, withdraw-and-reinjection rates, and natural groundwater flow velocities was quantitatively simulated using the finite difference simulator HST3D. Theoretical analyses of the simulated results were also made. The simulated results of the single well system indicate that neither the permeability nor the porosity of a homogeneous aquifer has significant effect on the temperature of the production segment provided that the production and injection capability of each well in the aquifers involved can meet the designed value. If there exists a lower permeable interlayer, compared with the main aquifer, between the production and injection segments, the temperature changes of the production segment will decrease. The thicker the interlayer and the lower the interlayer permeability, the longer the thermal breakthrough time of the production segment and the smaller the temperature changes of the production segment. According to the above modeling, it can also be found that with the increase of the aquifer thickness, the distance between the production and injection screens, and/or the regional groundwater flow velocity, and/or the decrease of the production-and-reinjection rate, the temperature changes of the production segment decline. For an aquifer of a constant thickness, continuously increase the screen lengths of production and injection segments may lead to the decrease of the distance between the production and injection screens, and the temperature changes of the production segment will increase, consequently.According to the simulation results of the single well system, the parameters, that can cause significant influence on heat transfer as well as geo-temperature field, were chosen for doublet system simulation. It is indicated that the temperature changes of the pumping well will decrease as the aquifer thickness, the distance between the well pair and/or the screen lengths of the doublet increase. In the case of a low permeable interlayer embedding in the main aquifer, if the screens of the pumping and the injection wells are installed respectively below and above the interlayer, the temperature changes of the pumping well will be smaller than that without the interlay. The lower the permeability of the interlayer, the smaller the temperature changes. The simulation results also indicate that the lower the pumping-and-reinjection rate, the greater the temperature changes of the pumping well. It can also be found that if the producer and the injector are chosen reasonably, the temperature changes of the pumping well will decline as the regional groundwater flow velocity increases. Compared with the case that the groundwater flow direction is perpendicular to the well pair, if the regional flow is directed from the pumping well to the injection well, the temperature changes of the pumping well is relatively smaller.Based on the above simulation study, a case history was conducted using the data from an operating system in Beijing. By means of the conceptual model and the mathematical model, a 3-D simulation model was developed and the hydrogeological parameters and the thermal properties were calibrated. The calibrated model was used to predict the evolution of the geo-temperature field for the next five years. The simulation results indicate that the calibrated model can represent the hydrogeological conditions and the nature of the aquifers. It can also be found that the temperature fronts in high permeable aquifers move very fast and the radiuses of temperature influence are large. Comparatively, the temperature changes in clay layers are smaller and there is an obvious lag of the temperature changes. According to the current energy mining load, the temperature of the pumping wells will increase by 0.7°C at the end of the next five years. The above case study may provide reliable base for the scientific management of the operating system studied.

遵义市岩溶水地下水流和污染物运移数值模拟研究

地下水数值模拟技术已成为评估人类活动对地下水质和量的影响、评价地下水资源、预测地下水污染发展趋势等的最主要的方法和手段。喀斯特含水层由于含水介质和地下水流场的非均质性和各向异性，对其进行地下水流的数值模拟一直是水文地质学界的难题。 遵义市地处我国西南喀斯特发育区，为贵州省第二大工业城市，属重度缺水地区，地下水资源的开发利用极大缓解了区内的缺水危机。但长期以来，由于对地下水资源的开发利用缺乏合理的统筹规划和强有力的管理，引发了一些环境地质问题，如地下水降落漏斗、岩溶塌陷、地下水质恶化等。因此选择遵义市进行地下水流和污染物运移数值模拟研究具有理论和实际意义。 通过ArcGIS平台建立了研究区的水文地质信息数据库，对研究区地下水的水位动态以及水化学特征做了简单分析。概要总结和阐述了高桥－河溪坝块段的自然地理、地质概况和水文地质条件，建立了水文地质概念模型；在水文地质概念模型的基础上，利用Groundwater Vistas软件建立了枯水期和丰水期的二维非均质各向异性稳定流模型，三维有限差分地下水流模拟程序MODFLOW用于模拟地下水水流，三维溶质运移模块MT3DMS用于模拟污染物在对流弥散情况下的迁移。根据分析和模拟结果可以得出如下几点结论： 1、基岩裂隙水水位峰值滞后大气降水峰值2～3个月，属渐峰型动态；岩溶水水位、地下河出口和泉流量变化步调与降水强度一致，对降水响应敏感。 2、对NH4+、NO3-、NO2-、SO42-、Mn五种组分含量进行了时空分析，结果表明地下水污染物的含量可能受人为活动输入物质的不均匀性和降雨等各方面因素控制，各组分每年的污染面积不一致，没有明显的规律性；受污染的一般是岩溶水，尤其是在石灰岩溶洞、地下河强烈发育而三废排放量大的居民集中地区面积较大。 3、为了有效地进行地下水资源管理，论文对高桥－河溪坝岩溶含水系统进行了一定的概化，将岩溶含水介质近似作为等价多孔介质（Equivalent Porous Media, EPM）模型来进行研究，采用MODFLOW的六个子程序模拟含水层系统的源汇项：降水子程序包RCH模拟降水入渗量、井流子程序包WEL模拟抽水量、通用水头子程序包GHB模拟侧向补给/排泄量、排水沟渠子程序包DRN模拟地下河出口流量、河流子程序包RIV模拟河流与地下水的交换量和已知水头边界子程序包CHD。从水位观测点和地下水位等势面两者结合来校正模型，结果表明能够达到相应国家标准规定的要求。因此EPM模型是可以适用于我国西南喀斯特地区的地下水流模拟的。 4、通过稳定流模型识别了枯水期和丰水期的渗透系数。在高桥和茅草铺附近渗透系数较高，枯水期介于100～400 m/d，而丰水期在高桥最高可达到3220m/d；其余单元渗透系数低于100 m/d，大多数小于10m/d。总体来说，由于丰水期含水层的饱水度大，渗透系数要高于枯水期。 5、通过地下水均衡计算，确定了各补给项和排泄项的水量。枯水期最重要的补给来源是研究区东北角的侧向补给量，占总补给量的70％，人工开采是最大的排泄项；丰水期最重要的补给源是西部的已知水头边界，占总补给量的49％，东北角的侧向补给量是第二补给源，占39％，地下河出口是最主要的排泄方式，达到排泄总量的74％。 6、对水文地质参数和源汇项敏感度分析的结果表明，不管是枯水期还是丰水期，对研究区水位影响最大的是渗透系数，外部源汇项中则是抽水量对地下水流形态的影响最大。 7、研究区岩溶地下水流速很大，污染物的运移是一个对流占绝对优势的问题，弥散的作用则相对很小。通过在茅草铺地区假设污染源，用MT3DMS程序模拟了地下水污染物在时间和空间上的迁移特征。结果发现：污染羽的形状和扩散方位主要受地下水流场的控制，而污染物的浓度与水量多少相关。

Relationships between precipitation, soil water and groundwater at Chongling catchment with the typical vegetation cover in the Taihang mountainous region, China

Vegetation cover plays an important role in the process of evaporation and infiltration. To explore the relationships between precipitation, soil water and groundwater in Taihang mountainous region, China, precipitation, soil water and water table were observed from 2004 to 2006, and precipitation, soil water and groundwater were sampled in 2004 and 2005 for oxygen-18 and deuterium analysis at Chongling catchment. The soil water was sampled at three sites covered by grass (Carex humilis and Carex lanceolata), acacia and arborvitae respectively. Precipitation is mainly concentrated in rainy seasons and has no significant spatial variance in study area. The stable isotopic compositions are enriched in precipitation and soil water due to the evaporation. The analysis of soil water potential and isotopic profiles shows that evaporation of soil water under arborvitae cover is weaker than under grass and acacia, while soil water evaporation under grass and acacia showed no significant difference. Both delta O-18 profiles and soil water potential dynamics reveal that the soil under acacia allows the most rapid infiltration rate, which may be related to preferential flow. In the process of infiltration after a rainstorm, antecedent water still takes up over 30% of water in the topsoil. The soil water between depths of 0-115 cm under grass has a residence time of about 20 days in the rainy season. Groundwater recharge from precipitation mainly occurs in the rainy season, especially when rainstorms or successive heavy rain events happen.