基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

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

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.

巴丹吉林沙漠及其南缘地区的水化学和环境同位素特征与环境演变

The Badain Jaran Desert lies on the Alashan Plateau in western Inner Mongolia. Because of huge dunes, permanent lakes and on the northern fringe of the Asian summer monsoon, the Badain Jaran Desert has been drawing attentions of many experts. And they have made great progress in dune’s geomorphology, botany in desert, paleoclimate change and other study areas. We analyzed environmental isotope and ion chemistry in lakes and groundwater of the desert and southeastern area, and collected some other evidences from 14C dating, fossils and archeology. According to chemical analysis, we discuss the difference spatial character of ion chemistry and environmental isotope in lakes and groundwater of the desert and adjacent. Contrasting with ion chemistry and isotope results in other arid area, we argue origin of groundwater and lakes in the desert area, and get a preliminary understanding of desert lakes’ evolution during Holocene. Some main conclusions were drawn as follows: 1. It has a obvious difference in hydrophysical parameters between lakes and groundwater in the desert and margin. 2. The results of ion analysis show that Na+ and Cl- are dominant in most lakes of the desert. Meanwhile, Na+ 、Cl- and HCO3- are dominant in groundwater of the desert and adjacent, and alsoMg2+、Ca2+、and NO3- have more percentage than in lakes. 3. Owing to different solubilities, the conten of main ions in water varies with the content of TDS. Whereas the content of TDS is over 100 g/L, the content of SO42-、HCO3-、Mg2+and Ca2+ in lakes descend. 4. The result of isotope analyzing indicate the lakes and groundwater in southeast desert have a similar vaporing trend with the groundwater in the southeast margin of the desert. It imply there would have some kind of contact between groundwater in margin and lakes of southeast desert. 5. Contrasting with isotope results of groundwater in other arid area, it show that the groundwater in the desert and Yabulai area should be phreatic water which have a low water table. Therefore, we conclude that the groundwater in southeast part of the desert and southern margin mainly are recharged by precipitation of local abundant rainfall and groundwater of low mountain of southern area. 6. And all of these evidences, which are different from salinity, the content of CO32- and geological data, show that the bigger northern lake group and southeastern lake group in the desert have different groundwater replenishing system because a fold belt lie between of the two group lakes and obstruct them in landform. and HCO3- 7. The 14C dating results of fossil and lacustrine deposits show that there maybe have a wider range of shoreline during early and middle Holocene than today. 8. By the discovery and study of some pieces of pottery and fine stoneware, we preliminary conclude that there maybe have some certain amount of early human activities in the Badain Jaran Desert.

基于盆地地下水模型与山区地表水模型耦合的济源盆地地下水模拟研究

The aim of this study is to increase the precision of groundwater modeling. The way is use the distributed model calculate the mountain basin groundwater lateral discharge and the river runoff. With appropriate technique help, the groundwater model can couple the distributed model results. This paper’s study object is makeing the distributed hydrological model HEC-HMS coupled to the popular groundwater model Visual MODFLOW. The application example is Jiyuan basin which is a typical basin of North China. HEC-HMS can calculate the surface runoff and subsurface runoff at mountain-pass. The subsurface runoff can turn to recharge well straightly. The water level - runoff course and Trial method is used to back analyze the parameters of surface runoff to Visual MODFLOW. So the distributed hydrological model can coupled to the groundwater model. The research proved that base on couple the distributed surface water model the groundwater model’s results are notability improved. The example is Jiyuan basin where use the distributed model coupled to the groundwater model. On the base of the coupled model applied to Jiyuan basin groundwater modeling. The paper estimates the groundwater change in the study area. Then, by use the water resources integrated planning results, the article calculate the basin groundwater can be development and utilization quantity and potential.

基于流域—库区—坝区模型“三重边界嵌套”方法的河南济源蟒河口水库渗漏数值模拟研究

Constructing reservoir in mountain area, in order to avoid dangerous reservoir and ensure long-term stable reservoir storage, study on reservoir leakage has been one of research hotspots. How to accurately predict the capacity of leakage and the effect of anti-seepage after reservoir impoundment under different anti-seepage measures is one of the most difficult problems of study on reservoir leakage. The paper took Manghekou Reservoir as a case study, which is located in Jiyuan City, Henan Province, China. Based on the system analysis of the geological and hydrogeological characters in the research area, the hydrogeological conceptual model was established. Then, according to the triplex boundary nested grid model method, using three-dimensional numerical simulation technique, the 3-D groundwater seepage flow numerical model was established in the study area. Based on the model, the change of seepage flow field, the capacity of leakage and the effect of anti-seepage are predicted after Manghekou reservoir impoundment under different condition. The reasonable and reliable anti-seepage measure is selected according to the model. All these provides scientific basis for the design and assessment of anti-seepage measure for the reservoir in mountain area.

海底金属矿床充水条件及矿坑突水的断层效应研究

The Xinli mine area of Sanshandao mine is adjacent to the Bohai Sea and its main exploitable ore deposit occurs in the undersea rock mass. The mine is the biggest undersea gold mine of China after production. The mine area faces a latent danger of water bursting, even sudden seawater inrush. There is no mature experience in undersea mining in China so far. The vein ore deposit is located in the lower wall of a fault; its possible groundwater sources mainly include bittern, Quaternary pore water and modern seawater. To ensure the safety of undersea mining, to survey the flooding conditions of the ore deposit using proper measures and study the potential seawater inrush pattern are the key technical problems. With the Xinli mine area as a case study, the engineering geological conditions of the Xinli mine area are surveyed in situ, the regional structural pattern and rock mass framework characteristics are found out, the distribution of the structural planes are modeled by a Monte Carlo method and the connectivity coefficients of rock mass structural planes are calculated. The regional hydro-geological conditions are analyzed and the in-situ hydro-geological investigation and sampling are performed in detail, the hydrochemistry and isotopes testing and groundwater dynamic monitoring are conducted, the recharge, runoff, discharge conditions are specified and the sources of flooding are distinguished. Some indices are selected from the testing results to calculate the proportion of each source in some water discharge points and in the whole water discharge of the Xinli mine area. The temporal and spatial variations of each water source of the whole ore deposit flooding are analyzed. According to the special project conditions in the Xinli mine area, the permeability coefficient tensors of the rock mass in Xinli mine area are calculated based on a fracture geometry measurement method, in terms of the connectivity and a few hydraulic testing results, a modified synthetic permeability coefficient are calculated. The hydro-geological conceptual and mathematical model are established，the water yield of mine is predicted using Visual Modflow code. The spreading law of surrounding rock mass deformation and secondary stress are studied by numerical analysis; the intrinsic mechanism of the faults slip caused by the excavation of ore deposit is analyzed. The results show that the development of surrounding rock mass deformation and secondary stress of vein ore deposit in the lower wall of a fault, is different from that in a thick-big ore deposit. The secondary stress caused by the excavation of vein ore deposit in the lower wall of a fault, is mainly distributed in the upper wall of the fault, one surface subsidence center will occur. The influences of fault on the rock mass movement, secondary stress and hydro-geological structures are analyzed; the secondary stress is blocked by the fault and the tensile stress concentration occurs in the rock mass near the fault, the original water blocking structure is destructed and the permeable structure is reconstructed, the primary structural planes begin to expand and newborn fissures occur, so the permeability of the original permeable structure is greatly enhanced, so the water bursting will probably occur. Based on this knowledge, the possible water inrush pattern and position of the Xinli mine area are predicted. Some computer programs are developed using object-oriented design method under the development platform Visual Studio.Net. These programs include a Monte Carlo simulation procedure, a joint diagrammatizing procedure, a structural planes connectivity coefficient calculating procedure, a permeability tensor calculating procedure, a water chemical formula edit and water source fixture conditions calculating procedure. A new computer mapping algorithm of joint iso-density diagram is raised. Based on the powerful spatial data management and icon functions of Geographic Information System, the pit water discharge dynamic monitoring data management information systems are established with ArcView.

巴丹吉林沙漠沙丘、植被和湖泊演化

The formation and evolution of magedunes and lakes in Badain Jaran desert is focused by scholars for a long time. In this paper, detailed investigation was made with the Remote Sensing, Geographic Information Systems, climatological data of research region to study the evolution of magedunes, vegetation and lakes. The main conclusions are presented as follows: (1) The magedunes are mainly located in south and northeast of the desert. The density of magedunes decreases from south to north, east to west, and greatest in southeast of the desert. Most of magedunes strike NE30°~60°, and which are NE-SW and NNE-SSW directions. The distribution of magedunes’ strikes indicates that the paleowind of the desert is northwest and west wind. (2) The vegetation cover is low in the Badain Jaran Desert, but there are some well-vegetated areas in the desert, and the value of NDVI is 0.6107. In southeast of the desert, the vegetation degradated seriously from 1987 to 2000. (3) At about 2000, the total area of permanent lakes was 2231.64 hm2 (±5.9%) in southeast of Badain Jaran Desert. From 1960 to about 2000, the lakes changed as followings: in the early 30 years, the area of lakes increased, but after 1990 it decreased. When lakes were large, the toral area of lakes is 20648 hm2 in Badain Jaran Desert, much larger than that of present lakes. (4) From 1961 to 2001, the mean annual evaporation of lakes was 1039.8mm (±10%) in Badain Jaran Desert, and the difference between years was large. The monthly evaporation of lakes was mainly affected by rainfall, temperature. And it was largest in June, 140.9mm(±10%). (5) Based on water balance calculation of present lakes and old lakes in southeast of Badain Jaran Desert, the total precipitation is more than the total evaporation in this area. So it is probable that the precipitation in the lake area supplies water to the lakes, and also to groundwater in west and north part of the desert.

陕西关中盆地地热水循环的同位素研究

Geothermal resource is rich in Guanzhong Basin, but as to its cycle characteristic, there has been lack of systematic study so far. Blind exploitations lead to water-temperature reducing, the decrease of spring flow rate and so on. Based on groundwater system and hydrogeological and hydrological geochemical theory, this paper studied the recycling type of geothermal water and analyzed the resources of dissolved inorganic carbon (DIC) and sulfate. The origin of the internal geothermal water is ice and snow in Qinling Mountain and Liupan Mountain above 1400m. The precipitation and surface water entered the deep part of the basin along piedmont faults, heated and water-expansion increased. The karst groundwater comes from meteoric water of the bare carbonate rock area in the North Mountains. Geothermal-water DIC mainly came from the dissolution of carbonate rock in the deep part of Guanzhong Basin, sulfate of Xi’an depression and Lishan salient came from the dissolution of continental evaporate , and sulfate of Gushi depression and Xianli salient came from co-dissolution of continental and marine evaporate. The above results supply science basis for reasonable exploitation and sustainable utilization of the geothermal water in Guanzhong Basin.

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

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.

滑坡型泥石流转化机理研究—以四川省达县青宁乡滑坡为例

Slide-debris flow is debris flow which is transformed from landslide consecutively in a short time, it comprises of two phases: First, Landslide starts to slide; Second, landslide changes to debris flow. Slide-debris flow which brings great property and life loss happens frequently at home and abroad. In order to forecast the happening possibility and scope of slide-debris flow, transfromation mechanism of Slide-debris flow must be studied. Research on transformation mechanism of slide-debris flow is intersectant science of landslide kinetics and debris flow starting theory, It is a fringe problem as well as front problem of geological hazard. This paper takes Qingning slide-debris flow in Da County, Sichuan Province for example and has studied the mechanism of its instability and transfromation into debris flow through indoor test (including usual soil test and ring shear test) and digital modeling method.The research gets the following conclusions. Qingning Landslide took place mainly because of confined water head arising from rainfall infiltration. Before Landslide occurring, it rained continuously for 22 days, accumulated precipitation arrived at 521.6mm.Investigation shows that strata of Qingning Landslide contains quaternary loose accumulation, slip soil and highly weathered bedrock, which is a good condition for formation of confined water in the slope. Further more, groundwater seepage in the slope body and corresponding slope safety factor before landslide occurring have been computed through finite element method. The result shows that because of infiltration of rainfall, confined water head in the slope arose sharply, accordingly, the safety factor of the slope declined quickly. The result also shows that force put on the slide body by the rock mass detached from Dazhaiyan mountain was the direct factor for landslide occurring. Qingning slide-debris transformation mode has been summarized, the process the landslide changed into debris flow is divided into three phases in the prospective of macroscopic geological condition: landslide occurring, transformation and debris flow. Landslide occurring phase is from slope’ local creeping slide to Landslide occurring; transformation phase contains slide body sliding on the slide bed after slide occurring and sliding on the slope after shearing opening; debris flow phase is that slide body breaks up completely and flows downward into the ditches. The transformation mechanism of Qingning slide-debris flow has been studied through indoor ring shear test of slip soil. The result shows that transformation mechanism contains two points: first, during slide body sliding on the slide bed and slope after shearing opening, shearing shrinkage, grain crushing and grain layering brought about declining of its volume and produced excess pore water pressure, and because producing velocity of excess pore water pressure is much greater than its dissipating velocity, shear strength of slide body decreased sharply because of accumulated pore water pressure. Second, grains crushing and grains layering during slide body sliding brought about thick liquefied layer at the bottom of the slidebody, liquefied layer contained high water content and its shear strength was very low, its thickness increased as the sliding displacement increasing. Liquefied layer makes slide body sliding fast and easily break down to debris flow. Excess pore water pressure and liquefied layer made shear strength of slidebody became very low, furthermore, water in the pit of slope joining in the slidebody was also a facter that made slidebody accelerate the transformation. Influence of slide body thickness and fine grains content to transformation of slide-debris flow has been studied through ring shear test. The result reaches two conclusions. First, thickness of slide body affects transformation of slide-debris flow by two ways, porewater pressure and effect of “soft base” increases as thickness of slide body increasing.so the thicker slide body is ,the easier transformation is. Second, actual dissipating velocity of porewater pressure should be considered when studying the influence of fine grains content to tranformation of slide-debris flow. There should be a critical content of fine grains which makes the difference of producing and dissipating velocity of water pore pressre greatest, this value is the best for slide-debris transformation. The whole process of slide-debris flow transformation is reproduced through discrete element method. Transformation mechanism of slide-debris flow is studied through monitoring various parameters including pore water pressure, grain crushing and grain layering in the slide body during the transformation. The result confirms and supplements the transformation mechanism of slide-debris flow got from ring shear test well.

关中盆地东部浅层地下水水化学与CFC年龄研究

Hydrochemistry, isotope and CFCs were used to determine groundwater transport in the eastern part of the Guanzhong Basin. In this paper, we systematically collected water samples and measured major ions, 2H-18O and CFCs in surface water and shallow groundwater. Groundwater in this region can be divided into three categories based on total dissolved solids (TDS): fresh water with TDS < 1g/L, brackish water with TDS between 1～3g/L, and saline water with TDS > 3g/L. Saline water is mainly located in the north of the Wei River, and saline groundwater is not in the south. Tributaries in the south of the Wei River and underlain groundwater had similar 2H-18O values, indicating a close hydraulic connection between them. Tributaries in the north of the Wei River characterized certain extent of evaporation, and 2H-18O values deviated to a differing extent between surface water and groundwater, indicating that surface water in the north bank of the Wei River has little hydraulic connection with underlain groundwater. The CFCs age of groundwater from the piedmont recharge area was young, and became older toward the Wei River valley. Vertically, the CFCs age of groundwater increased with well depth. The shallow groundwater is mainly composed of young water with ages < 60 years and old water with ages > 60 years. Young water is in a larger proportion. The NO3-N concentration positively correlates with the CFC-12 concentration in the groundwater, which indicates that young water is subjected to be contaminated. Keyword: Guanzhong Basin , shallow groundwater, isotope, CFC

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

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.