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

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.

降雨触发滑坡的非线性模型与斜坡演化的混沌效应

This paper analyzes landsliding process by nonlinear theories, especially the influence mechanism of external factors (such as rainfall and groundwater) on slope evolution. The author investigates landslide as a consequence of the catastrophic slide of initially stationary or creeping slope triggered by a small perturbation. A fully catastrophe analysis is done for all possible scenarios when a continuous change is imposed to the control parameters. As the slip surface continues and erosion due to rainfall occurs, control parameters of the slip surface may evolve such that a previously stable slope may become unstable (e.g. catastrophe occurs), when a small perturbation is imposed. Thus the present analysis offers a plausible explanation to why slope failure occurs at a particular rainfall, which is not the largest in the history of the slope. It is found, by analysis on the nonlinear dynamical model of the evolution process of slope built, that the relationship between the action of external environment factors and the response of the slope system is complicatedly nonlinear. When the nonlinear action of slope itself is equivalent to the acting ability of external environment, the chaotic phenomenon appears in the evolution process of slope, and its route leading to chaos is realized with bifurcation of period-doublings. On the basis of displacement time series of the slope, a nonlinear dynamic model is set up by improved Backus generalized linear inversion theory in this paper. Due to the equivalence between autonomous gradient system and catastrophe model, a standard cusp catastrophe model can be obtained through variable substitution. The method is applied to displacement data of Huangci landslide and Wolongsi landslide, to show how slopes evolve before landsliding. There is convincing statistical evidence to believe that the nonlinear dynamic model can make satisfied prediction results. Most important of all, we find that there is a sudden fall of D, which indicates the occurrence of catastrophe (when D=0).

煤层气储层工程的地下水行为研究

The practice of coalbed methane development from home and abroad demonstrated Hydrogeological factor is one of the important geological factors influencing the coalbed methane productivity. The grasp of groundwater behavior feature is the prerequisite to success development of coalbed methane. Through researching the mechanism by which hydrodynamics factors control the storage and transportation of coalfen methane, the ground- water behavior reflecting the feature of coalbed, and mathematics model describing the production process of coalbed methane, this paper devoted to finding the law of groundwater behavior during the course of storage and production and gave hydrogeology theoretical support to the development of coalbed methane. This paper firstly studied hydrodynamic factors influencing the productivity of coalbed methane, based on the analysis of the relative feature of coalbed methane and that of it's reservoir. The productivity of coalbed methane is controlled by reservoir pressure、permeability and recharge conditions. Reservoir pressure, the key factor controlling gas content of coalbed, is ruled by the history of hydrodynamic and current hydrogeological conditions. It indirectly controls the poductivity through influencing the permeability. The permeability of coalbed is the direct factor controlling the productivity. The recharge conditions controls the productivity through influencing initial reservoir pressure and the descend of reservoir pressure during development of coalbed methane. The field of hydrodynamic and the field of hydrochemistry can be used to identified the flow model of groundwater and the coalbed feature can be deducted by the hydraulic gradient、pressure compartment and hydrochemistry. The production of coalbed methane is a complex physical process which including the mutual action between water、solid and gas. This paper studied the mechanism of water-solid action and that of water-gas action, conducted the controlling equation describing the complex process and gave the corresponding mathematics model with its solution by finite-Element method. Finally, this paper analysised the prospective of coalbed methane development of the south section of Hongguo area in Yizikong basin and put emphasis on the analysis of productivity of liangshan and jingzhuping blocks.

工程地质三维建模与可视化技术应用开发

Synthetic Geology Information System (SGIS) is an important constituent part of the theory of Engineering Geomechanics Mate-Synthetic (EGMS), and is the information system more suited for the collection, storage, management, analysis and processing to the information coming from engineering geology,' geological engineering and geotechnical engineering. Its contents involve various works and methods of the investigation, design, and construction in different stages of the geological engineering. Engineering geological and three-dimensional modeling and visualization is the fundamental part of the SGIS, and is a theory, method and technique by which, adopting the computer graphics and image processing techniques, the data derived from engineering geological survey and the calculated results obtained from the geomechanical numerical simulation and analysis are converted to the graphics and images displayed on the computer screen and can be processed interactively. In this paper, the significance and realizing approaches of the three-dimensional modeling and visualization for the complex geological mass in the engineering geology are discussed and the methods of taking advantage of the interpolation and fitting for the scattered and field-surveyed data to simulate the geological layers, such as the topography and earth surface, the groundwater table and the stratum boundary, are researched into. At the mean time, in mind the characteristics of the structure of the basic data for three-dimensional modeling, its visual management can be resolved into the engineering surveyed database management module, plot parameter management module and data output module and the requirement for basic data management can be fulfilled. In the paper, the establishment and development of the three-dimensional geological information system are probed tentatively, and an instance of three-dimensional visual Engineering Distribution Information System (EDIS), theConstruction Management Information System for an airport, in which the functions, such as the real-time browse among the three-dimensional virtual-reality landscapes of the airport construction from start to finish, the information query to the airport facility and the building in the housing district and the recording and playback of the animation sets for the browse and the takeoff and landing of the planes, is developed by applying the component-mode three-dimensional virtual-reality geological information system (GIS) software development kits (SDK), so the three-dimensional visual management platform is provided for the airport construction. Moreover, in the gaper, integrated with the three-dimensional topography visualization and its application in the Sichuan-Tibet Highways, the method of the digital elevation model (DEM) data collection from the topographic maps is described, and the three-dimensional visualization and the roaming about the terrain along the highway are achieved through computer language programming. Understanding to the important role played by the varied and unique topographical condition in the gestation and germination of the highly-dense, frequently-arising and severely-endangered geological hazards can be deepened.

东天山地区金矿床成矿流体地球化学及其成矿机制

The mechanism of gold ore formation in the eastern Tianshan Mountains, Xinjiang Uygur Autonomous Region, that has been dealt with from various aspects, remains unclear. On the basis of investigations of regional geology, ore deposit geology, and microscopic observations of ores and related rocks of the Jinwozi, the 210, and the Mazhuangshan gold ore deposits, this thesis made a systematic research on the microthermometry of gangue quartz-hosted fluid inclusions, gas, liquid ion and rare earth element compositions and hydrogen, oxygen isotope compositions of sulfide- and quartz-hosted fluid inclusions, and sulfur and lead isotope compositions of sulfide ore minerals from the major ores in the three deposits. On the basis of the above synthetic studies, sources of ore-forming fluids and metals, and mechanism of gold ore formation in the region were discussed. Gas compositions of pyrite- and sphalerite-hosted fluid inclusions were first analyzed in this thesis. Compared with gangue quartz-hosted fluid inclusions, the sulfide-hosted ones are richer in gaseous species CO2, CO, and CH4 etc. Both gas and liquid CO2 are commonly observed in fluid inclusions, whereas halite daughter minerals rarely occur. Ore-forming fluids for the three gold ore deposits are characteristically of medium to low temperatures, medium to low salinities, are rich in CO2 and Na+, K+, Cl" ions. Gas covariation diagrams exhibit linear trends that are interpreted as reflecting mixing between the magamtic fluid and meteoric-derived groundwater. Regarding rare earth element compositions, the Jinwozi and the 210 deposits show moderate to strong LREE/HREE fractionations with negative Eu anomalies. However, the Mazhuangshan deposit shows little LREE/HREE fractionation with positive Eu anomalies. Hydrogen and oxygen isotope compositions of pyrite-hosted fluid inclusions that were first analyzed in this thesis indicate the presence of magmatic water. Hydrogen and oxygen isotope compositions of pyrite- and quartz-hosted fluid inclusions suggest mixing between magmatic water and meteoric-derived groundwater. Sulfur and lead isotope compositions of sulfide ore minerals indicate multi-sources for the metallogenetic materials that range from the crust to the mantle. On the basis of the above synthetic studies, genesis of the gold ore deposits in the eastern Tianshan Mountains was approached. From the Middle-Late Hercynian to Early Indosinian, geodynamic regime of the region was changing from the collisional compression to the post-collisional extension. During the period, magmas were derived from the crust and the mantle and carried metallogenetic materials. Magma intrusion in the upper crust released the magmatic fluids, and drove circulation of groundwater. Mixing of magmatic fluid with groundwater, and extraction of metallogenetic materials from the country rocks are the mechanism for the gold ore formation in the eastern Tianshan Mountains.

抽水引起岩溶塌陷的机理及非线性预测研究

Karst collapse is one of the most important engineering geology hazards in Karst district, which seriously endangers the living of humankind and the environment around us, as well as the natural resources. Generally speaking, there exist three processes of overburden karst collapse:the formation of soil cavity, the expansion of soil cavity and the fall of the cavity roof. During these processes, groundwater is always the most active factor and plays a key role. Pumping will bring into the great change of groundwater in flow state, flowrate, frequency of fluctuation as well as hydraulic gradient and will speed the fall. Statistics shows that most of the man-made karst collapse are induced by pumping, so studying the mechanism of Karst collapse induced by pumping will provide theoretical base for the prediction and precaution of collapse. By theoretically studying the initial condition for the forming and expanding of a soil cavity, Spalling step by step the essential mechanism of Karst collapse induced by pumping is put forward. The catastrophe model for the collapse induced by pumping is set up to predict the fall probability of a cavity roof, and the criterion for the collapse is determined. Simultaneously, Karst collapse induced by pumping is predicted with manmade neural network theory. Finally, the appropriate precaution measurements for the collapse induced by pumping are provided. The creative opinions of the paper is following: The initial condition of forming a soil cavity is put forwarded as formula (4-1-5), (4-1-24),(4-1-25) and (4-1-27); which provide theoretical base for foreclosing the formation of a soil cavity and defending collapse. Spaliing step by step as the essential mechanism of Karst collapse induced by pumping is put forward. The spaliing force is defined as formula (4-2-15). The condition for the expanding of a soil cavity is that spaliing force is greater than tensile strength of soil. The stability of a soil cavity is first studied with catastrophe theory. It is concluded that the process of development up to ground collapse of a small cavity is continuous, however, the process of a big cavity is catastrophic. It is feasibility that the Karst collapse be predicted with manmade neural network theory as a new way.

水库滑坡复活机理及其发展趋势预测研究

Various numbers of ancient landslides of various scales are frequently distributed on both banks of reservoirs, especially large reservoirs, both in China and abroad. During inundation and operation of theses reservoirs, some of the landslides are reactivated, which caused losses of people's lives and properties to various extents, some even disasters. Systematic studies are, however, very few on the reservoir-induced reactivation mechanism and development tendency prediction. Based on investigation of reservoir-induced reactivation phenomena of ancient landslides and relevant existing research problems, a systematic study is carried out on the field identification, induced reactivation mechanism, development tendency prediction, risk decision-making and treatment of reservoir-related ancient landslides, through analysis of large numbers of engineering geological investigation results, scientific experimental and research results, in combination with prevention and treatment practices of reservoir-related landslides both in China and abroad, and a series of research results have been obtained. 1. On the basis of study of the distribution features, genesis mechanism of ancient landslides on river banks, a set of scientific methods are summarized on field identification of ancient landslides, and a significant method named "lithologic sequence method" or "indicator layer method", is proposed, which is proved to be very useful. 2. A detail study is made on the reservoir-induced hydraulic effects and material mechanic effects (or softening effects) on the ancient landslide through model and case studies, which concludes that the magnitude and properties of reservoir-induced hydraulic effects are related to the shapes of sliding planes, water content and permeability of landslide materials and variation rate and magnitude of reservoir levels; the magnitude of material mechanic effects are related to the material composition (including mineral composition and grain size), natural water content and saturation state of sliding zones. Also a sensitive analysis is made on the factors that are related to the stabilities of the landslides, which indicate that the stability of a landslide is more sensitive to the groundwater head h_w in the slides and the inner friction angleψof sliding zones than others. 3. The joint inducing mechanism of rainfall and reservoir is also discussed in the paper through model analysis and case study, which proves that reservoir inundation increases firstly the sensitivity of a landslides to rainfall through reduction of its stability or cracking deformation which will increase the rainfall infiltration to the slide body, and then rainfall triggers reactivation or intensifies the reservoir-induced deformation of a landslide. 4. Based on rheologic test results of sliding zones of several reservoir-related ancient landslides, the rheologic characteristics of sliding zones have been discussed in detail and several typical rheologic models have been set up, which well explains the dynamic process of slide deformation. The response types to reservoir inundation and development tendency of reservoir -related ancient landslide are discussed in the paper based on field investigation results. And prediction methods for reservoir-related landslides have been studied based on the Mate-Synthetic principle of quantitative and qualitative analysis, as well as combination of computation and internal mechanism analysis, and a rheologic analytical method is proposed which is proved very useful for prediction of the landslide development tendency. 6. In disaster-prevention and treatment of reservoir-related landslides, risk decision-making has been proved very significant both in engineering and economics. Based on the practices in disaster-prevention and treatment of reservoir-related landslides both in China and abroad, the disaster-prevention risk decision-making for reservoir-related landslides has been proposed in terms of philosophy, methods and procedures, and well put into practice. A summary is also made through case study of the experiences of treatment of reservoir-related landslides both in China and abroad in terms of principle, methods and technical lines. 7 A detail study is made as a case study of the reactivated Maoping ancient landslide on the left bank of the Geheyan Reservoir on Qingjiang river in Hubei province, China, including its field identification features, reservoir-induced reactivation characteristics and mechanism, development tendency prediction and proposed counter measures based on risk analysis.

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

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程序模拟了地下水污染物在时间和空间上的迁移特征。结果发现：污染羽的形状和扩散方位主要受地下水流场的控制，而污染物的浓度与水量多少相关。

Is China's fifth-largest inland lake to dry-up? Incorporated hydrological and satellite-based methods for forecasting Hulun Lake water levels

Hulun Lake, China’s fifth-largest inland lake, experienced severe declines in water level in the period of 2000-2010. This has prompted concerns whether the lake is drying up gradually. A multi-million US dollar engineering project to construct a water channel to transfer part of the river flow from a nearby river to maintain the water level was completed in August 2010. This study aimed to advance the understanding of the key processes controlling the lake water level variation over the last five decades, as well as investigate the impact of the river transfer engineering project on the water level. A water balance model was developed to investigate the lake water level variations over the last five decades, using hydrological and climatic data as well as satellite-based measurements and results from land surface modelling. The investigation reveals that the severe reduction of river discharge (- 364±64 mm/yr, ~70% of the five-decade average) into the lake was the key factor behind the decline of the lake water level between 2000 and 2010. The decline of river discharge was due to the reduction of total runoff from the lake watershed. This was a result of the reduction of soil moisture due to the decrease of precipitation (-49±45 mm/yr) over this period. The water budget calculation suggests that the groundwater component from the surrounding lake area as well as surface run off from the un-gauged area surrounding the lake contributed ~ net 210 Mm3/yr (equivalent to ~ 100 mm/yr) water inflows into the lake. The results also show that the water diversion project did prevent a further water level decline of over 0.5 m by the end of 2012. Overall, the monthly water balance model gave an excellent prediction of the lake water level fluctuation over the last five decades and can be a useful tool to manage lake water resources in the future.

Potential impacts of leakage from deep CO2 geosequestration on overlying freshwater aquifers.

Carbon Capture and Storage may use deep saline aquifers for CO(2) sequestration, but small CO(2) leakage could pose a risk to overlying fresh groundwater. We performed laboratory incubations of CO(2) infiltration under oxidizing conditions for >300 days on samples from four freshwater aquifers to 1) understand how CO(2) leakage affects freshwater quality; 2) develop selection criteria for deep sequestration sites based on inorganic metal contamination caused by CO(2) leaks to shallow aquifers; and 3) identify geochemical signatures for early detection criteria. After exposure to CO(2), water pH declines of 1-2 units were apparent in all aquifer samples. CO(2) caused concentrations of the alkali and alkaline earths and manganese, cobalt, nickel, and iron to increase by more than 2 orders of magnitude. Potentially dangerous uranium and barium increased throughout the entire experiment in some samples. Solid-phase metal mobility, carbonate buffering capacity, and redox state in the shallow overlying aquifers influence the impact of CO(2) leakage and should be considered when selecting deep geosequestration sites. Manganese, iron, calcium, and pH could be used as geochemical markers of a CO(2) leak, as their concentrations increase within 2 weeks of exposure to CO(2).

Climate change, water resources, and the politics of adaptation in the Middle East and North Africa

Through an examination of global climate change models combined with hydrological data on deteriorating water quality in the Middle East and North Africa (MENA), we elucidate the ways in which the MENA countries are vulnerable to climate-induced impacts on water resources. Adaptive governance strategies, however, remain a low priority for political leaderships in the MENA region. To date, most MENA governments have concentrated the bulk of their resources on large-scale supply side projects such as desalination, dam construction, inter-basin water transfers, tapping fossil groundwater aquifers, and importing virtual water. Because managing water demand, improving the efficiency of water use, and promoting conservation will be key ingredients in responding to climate-induced impacts on the water sector, we analyze the political, economic, and institutional drivers that have shaped governance responses. While the scholarly literature emphasizes the importance of social capital to adaptive governance, we find that many political leaders and water experts in the MENA rarely engage societal actors in considering water risks. We conclude that the key capacities for adaptive governance to water scarcity in MENA are underdeveloped. © 2010 Springer Science+Business Media B.V.

Methane contamination of drinking water accompanying gas-well drilling and hydraulic fracturing.

Directional drilling and hydraulic-fracturing technologies are dramatically increasing natural-gas extraction. In aquifers overlying the Marcellus and Utica shale formations of northeastern Pennsylvania and upstate New York, we document systematic evidence for methane contamination of drinking water associated with shale-gas extraction. In active gas-extraction areas (one or more gas wells within 1 km), average and maximum methane concentrations in drinking-water wells increased with proximity to the nearest gas well and were 19.2 and 64 mg CH(4) L(-1) (n = 26), a potential explosion hazard; in contrast, dissolved methane samples in neighboring nonextraction sites (no gas wells within 1 km) within similar geologic formations and hydrogeologic regimes averaged only 1.1 mg L(-1) (P < 0.05; n = 34). Average δ(13)C-CH(4) values of dissolved methane in shallow groundwater were significantly less negative for active than for nonactive sites (-37 ± 7‰ and -54 ± 11‰, respectively; P < 0.0001). These δ(13)C-CH(4) data, coupled with the ratios of methane-to-higher-chain hydrocarbons, and δ(2)H-CH(4) values, are consistent with deeper thermogenic methane sources such as the Marcellus and Utica shales at the active sites and matched gas geochemistry from gas wells nearby. In contrast, lower-concentration samples from shallow groundwater at nonactive sites had isotopic signatures reflecting a more biogenic or mixed biogenic/thermogenic methane source. We found no evidence for contamination of drinking-water samples with deep saline brines or fracturing fluids. We conclude that greater stewardship, data, and-possibly-regulation are needed to ensure the sustainable future of shale-gas extraction and to improve public confidence in its use.

Denitrification and inference of nitrogen sources in the karstic Floridan Aquifer

Aquifer denitrification is among the most poorly constrained fluxes in global and regional nitrogen budgets. The few direct measurements of denitrification in groundwaters provide limited information about its spatial and temporal variability, particularly at the scale of whole aquifers. Uncertainty in estimates of denitrification may also lead to underestimates of its effect on isotopic signatures of inorganic N, and thereby confound the inference of N source from these data. In this study, our objectives are to quantify the magnitude and variability of denitrification in the Upper Floridan Aquifer (UFA) and evaluate its effect on N isotopic signatures at the regional scale. Using dual noble gas tracers (Ne, Ar) to generate physical predictions of N2 gas concentrations for 112 observations from 61 UFA springs, we show that excess (i.e. denitrification-derived) N2 is highly variable in space and inversely correlated with dissolved oxygen (O2). Negative relationships between O2 and δ15N NO3 across a larger dataset of 113 springs, well-constrained isotopic fractionation coefficients, and strong 15N:18O covariation further support inferences of denitrification in this uniquely organic-matter-poor system. Despite relatively low average rates, denitrification accounted for 32 % of estimated aquifer N inputs across all sampled UFA springs. Back-calculations of source δ15N NO3 based on denitrification progression suggest that isotopically-enriched nitrate (NO3-) in many springs of the UFA reflects groundwater denitrification rather than urban- or animal-derived inputs. © Author(s) 2012.

Geochemical and isotopic evidence for seawater contamination of the hydrothermal system of Taal volcano, Luzon, the Philippines

The hydrologic structure of Taal Volcano has favored development of an extensive hydrothermal system whose prominent feature is the acidic Main Crater Lake (pH<3) lying in the center of an active vent complex, which is surrounded by a slightly alkaline caldera lake (Lake Taal). This peculiar situation makes Taal prone to frequent, and sometimes catastrophic, hydrovolcanic eruptions. Fumaroles, hot springs, and lake waters were sampled in 1991, 1992, and 1995 in order to develop a geochemical model for the hydrothermal system. The low-temperature fumarole compositions indicate strong interaction of magmatic vapors with the hydrothermal system under relatively oxidizing conditions. The thermal waters consist of highly, moderately, and weakly mineralized solutions, but none of them corresponds to either water-rock equilibrium or rock dissolution. The concentrated discharges have high Na contents (>3500 mg/kg) and low SO4/Cl ratios (<0.3). The Br/Cl ratio of most samples suggests incorporation of seawater into the hydrothermal system. Water and dissolved sulfate isotopic compositions reveal that the Main Crater Lake and spring discharges are derived from a deep parent fluid (T≃300°C), which is a mixture of seawater, volcanic water, and Lake Taal water. The volcanic end member is probably produced in the magmatic-hydrothermal environment during absorption of high-temperature gases into groundwater. Boiling and mixing of the parent water give rise to the range of chemical and isotopic characteristics observed in the thermal discharges. Incursion of seawater from the coastal region to the central part of the volcano is supported by the low water levels of the lakes and by the fact that Lake Taal was directly connected to the China sea until the sixteenth century. The depth to the seawater-meteoric water interface is calculated to be 80 and 160 m for the Main Crater Lake and Lake Taal, respectively. Additional data are required to infer the hydrologic structure of Taal. Geochemical surveillance of the Main Crater Lake using the SO4/Cl, Na/K, or Mg/Cl ratio cannot be applied straightforwardly due to the presence of seawater in the hydrothermal system.