Poly(phenylene sulfide-tetraaniline): The soluble conducting polyaniline analogue with well-defined structures

A novel conducting polymer poly(phenylene sulfide-tetraaniline) (PPSTEA), with tetraaniline (TA) and phenylene sulfide (PS) segments in its repeat unit, has been synthesized through an acid-induced polycondensation reaction of 4-methylsulfinylphenyl-capped tetraaniline. The new polymer, which represents the first soluble conducting polyaniline analogue with well-defined structure, has high molecular weight, good solubility in common solvents, and good film-forming properties. Its electrical property is analogous to polyaniline. The conductivity of preliminarily, protonic-doping PPSTEA is up to 10 degrees S/cm. This synthetic strategy appears to be general for developing novel well-defined polyaniline analogue containing much longer fixed conjugation length.

塔克拉玛干沙漠地球化学研究

The Taklamakan Desert, lying in the center Tarim Basin of sourthern Xinjiang, is the largest sand sea in China and well known in the world as its inclemency. For understanding the formation and evolution of the Taklamakan Desert, it is very important to identify the provenance of aeolian sediments in the extensive dune fields, but the opinions from earlier studies are quite different. In this study, we examined the major- and trace-element compositions, mineral compositions and grain-size distributions of some Quaternary aeolian and nonaeolian sediments collected from the Taklamakan Desert, together with the variation of chemical and mineralogical compositions of different grain-size fraction. At the same time, we also studied the chemistries of some natural water samples (river water and groundwater) with the items of TDS, pH, Alkalinity, conductivity and major cation and anion compositions. Our results of analysis show some significant opinions as follow: Most of the frequency-distribution curves of grain size of dune sand samples are simgle peak, but that of the river and lacustrine sediment are most double peak or multi-peak. The grain-size distribution of dunefield sand changed gradually from north to south with the major wind direction in large scale, but there are many differences in regional scale. The major, REE, trace element compositions and mineral compositions are very different among the coarse, fine fraction and bulk samples due to the influence of grain-size. Most of the fine fractions are geochemically homogenous, but the coarse fractions and bulk samples are heterogenous. All the surface and ground waters are limnetic or sub-salty, their chemical compositions are mainly controlled by rock-weathering and crystallization- evaporation processes, and mainly come from the evaporate, while the contributions of the carbonate and silicate are little, excluding the influence on oasis water by carbonate. The mineral compositions of selected samples are stable, mainly composed of the strongly resistant mineral types. The mineral maturity of them are more immature at whole compared with other sandy sediments in the world, and they have experienced less degree of chemical weathering and recycling, being lying in the early stage of continental weathering. Among these sediments, the river sediments are relatively primitive. The sources of these sediments are maybe mainly terrigenous, silicic and subaluminous/ metalunious rocks, such as the granodiorite and its metamorphic rock. The geochemical compositions of dunefield sand are similar with those of the river sediments and dune sands near the river way; There are not only the resemblances but also the differences on geochemistry and granularity between north and south dunesands; The sediments from same section have different age but same trace-element compositions; The sediments from the south edge of Tarim Basin are all somewhat geochemically similar with the palaeo-river-sediment on the south edge of studying area. The REE data support the idea that the south dunesands are a little older than the north dunesands, and the tectonic settings of source area are mainly active continental margin based on the major-element compositions, so they indicate that the sediment of Taklamakan Desert maybe come mainly from the rock-weathering production of north part of the Kunlun Mountains. Compared with the sands of other dune field in north of China, the sands of the Taklamakan Desert are distinct by REE composition, but similar with the Luochuan loess, center China, and the two sandy dusts of Beijing, eastern China.

水利水电岩溶发育深度的定量化理论研究

The distribution in outcrop of carbonate rocks in China is about 3.44×106km2 and accounts for over one third of the total area of the country, mainly in the Yangtze valley and the Yungui tableland ,and there is a great plenty of Water and hydropower resources in these carbonate area. A large number of projects will be built in order to develop the Water and hydropower resources. They are facing amount of complicated problems of karsts, specially the problem of the depths of karsts below water, which is a key problem that hasn’t been solved well theoretically. So, systematically research in quantitative analysis of the problem is necessary, which is very important to solve the inconsistency between current theories and facts, and is helped to foresee the depth of karsts below water and decrease the expense of the engineering. For the problem mentioned above, the thesis makes a detailed research on the length of corrosion of karsts water, based on the kinetics of corrosion, and founds the primary theory about the processes of four typical geologic elements, which are single fracture, single conduit, porous limestone and fractured limestone. And it has done a deep research upon the characterties and kinetics of corrosion of the four typical geologic elements by simulation using the programme edited by the author. The thesis also makes a discussion on the general model and process of the form of caves by using the founded theory. According to the characteristic of Water Resource and Hydropower engineering, the thesis creates three representative flow-dissolution models for three types of common geologic conditions under the induction of the theory kinetics of corrosion of karsts water. It has done quantitative research upon the process of dissolution and brought out primary theory about quantitative analysis of the depth of karsts below water as well. It found that the depth of karsts below water relates to the specific geologic conditions and the time of dissolution by simulation based on typical parameters and deep analysis of the result. That is to say there are karsts caves in any depths in flow areas of groundwater in specific geologic conditions and appropriate time. The thesis also discourses three basic problems that frequently encountered in the reconnaissance of Water Resource and Hydropower engineering, which are karsts base level, essential conditions of karsts, the comparability of caves and terrace, and indicates the importance of time in karsts. Finally, as an application, the karsts of the scheming Qianzhong water resource engineering is analyzed by the primary theory founded of quantitative analysis. The result of the application comes to a conclusion that the theory accords with the facts properly.

地下水位上升造成的黄土远程滑坡形成机理与数值模拟研究

Based on field survey, laboratory testing and numerical modeling, engineering characteristics of undisturbed loess and the mechanism of long-runout loess landslides caused by underground water level rise, as well as the formation conditions and spatial distribution of landslides, are systematically studied and analyzed. Loess landslides at south Plateau of Jingyang County are mainly classified as flowslide, slide and fall. Flowslide is the main type characteristic of high velocity, long runout and multi-stages. The steep relief composed of loose structured loess-old aged soil serials and the rise of groundwater table are the predominant conditions for landslides in the study area. To study loess mechanic poperties and loess landslides mechanisims, isotropically and anisotropically consolidated undrained compression(ICU and ACU) tests and constant-deviator-drained compression (CQD) tests were carried out on undisturbed samples. The results of undrained compression tests performed at the in-situ stress level show that the soils are of consistently strain-softening in the stress-strain relations and cause high excess pore pressure. The steady-state line and the potential region of instability are obtained from ICU and ACU test results. A necessary condition for liquefaction is that the soil state initially lies in or is brought into the potential instability region. In addition, a strong strain-softening model is also formed. CQD tests demonstrate that the mobilized friction angle is far less than the steady-state angle and that the soil experiences undrained contractive failure suddenly at very small strains when its stress path during drained loading tries to cross the potential instability region，thus validates the proposed instability region. Based on the location of the region of potential instability and the stress state of slope soil, a method of static liquefaction analysis is proposed for loess landslides caused by rise in groundwater table. Compared with other liquefaction analysis methods, this method overcomes the limitations inherent in conventional slope stability method and undrained brittleness index method. Triaxial tests composed of constant water content (CW) and wetting tests at constant deviator stress are performed on undisturbed unsaturated samples. The stress-strain relation of CW tests takes on strain-hardening behavior; The results of wetting tests at constant deviator stress designed to study the mechanics of failure of unsaturated loess caused by an increase in the degree of saturation (wetting) shows that a contractive failure occurs in the undisturbed samples. On the basis of the above triaxial test results, the initiation of static liquefaction is presented for long-runout loess landslides caused by rise in groundwater table, that is, the loess slope soil gradually transfer from unsaturated to saturated state under the infiltration of irrigation. A contractive failure occurs in the local region at very small strain by increasing the pore-water pressure at constant deviator stresses under drained conditons. It is the contractive failrue resulting from rise of pore pressure that leads to high excess pore pressure in the neighbour soil which reduces shear resistance of soil. The neighbour soils also fail due to the rapid increase in pore-water pressure. Thus a connected failure surface is developed quickly and a flowslide occurs. Based on the saturated-unsaturated seepage theory, transient seepage is computed using the finite element method on loess slope under groundwater table rise. Pore-water pressure distribution for every time step after irrigation are obtained. The phreatic surface in the slope increases with the groundwater table. Pore-water pressure distribution within 8m above the phreatic surface changes very quickly，but the water content and pore water pressure in the region ranging from 8m above the phreatic surface up to ground surface is almost not affected and the matric suction usually is kept at 100～120 kPa. Based on the results of laboratory tests and seepage flow analysis, the development process of loess landslide is modeled considering groundwater table rise. The shearing plastic zone first occurs at the slope toe where the soil is soaked for long term during rise in groundwater table. As irrigation continues, the shearing plastic zone gradually extends to the interior soils, with the results that the tensile plastic zone occurs at the slope crown. As time goes on, both the shearing plastic zone and tensile plastic zone continue to extend. Then a connected plastic zone is formed and fowslide occurs. In comparision to laboratory test results, the results of numerical simulation quite well verify the presented mechanism of static liquefaction of long-runout loess landslides caused by rise in groundwater table.

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

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.

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

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.

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

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

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

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.