环渤海花岗岩体地下油库建设的工程地质力学研究

As a kind of strategic resource，petroleum play an very important role in current social stability, economic development and state safety. Since 1993 China has turned from a net oil exporter into a net oil importer, the figure of imported oil increased from then on. In 2004 China's total energy consumption exceeded Japan’s, and ranked in the second place, just inferior to America. Today China is the world’s third-largest importing nation, accounting for 6% of world imports and 8% of world consumption. Comparing with other strategic petroleum reserve schemes, underground oil storage possess many advantages, such as security, economy, less pollution, save land, suited for strategic reserve and so on, so it is the most ideal form for strategic petroleum reserve. In the background of China Strategic Petroleum Reserve Program started just now, this paper choose Circum-Bo sea region as a study area, and do some system study on the underground oil storage caverns constructed in inter-large granite rock masses in Circum-Bo sea region. On the foundation of a great amount of information come from both home and abroad, firstly this paper analysed the principle, economy, cavern shape, profile dimension, and gain some cognizances and logos, as follows: ①Hard rock mass such as granite is the major rock, in which underground oil storage are constructed; ②Unlined underground oil storage caverns had been wide spread used as a sort of oil storage form abroad, there already exist a suit of skilled experience and technologies to prevent oil product from leaking; ③Compared with surface tanks, underground oil storage cavern possess predominance in economy clearly. In general, it will be more economical when the storage capacity exceed 50000m3. The quality of rock mass is the most important factor for underground storage cost, however such as hydrogeology, storage capacity, the number of storage galleries, the length, storage product, mechanical equipments, geographic location also influent the cost. In designed depth of the underground storage, the rock mass of Jinzhou mainly belong to class Ⅱ, but parts with dykes, clayization alteration, and dense joints are Ⅲ, Ⅳ; ④Now, there are few underground oil storages span more than 25m in both abroad and home. The examples of some ancient underground works and modern underground excavation with wide span surely give us many precious elicitations to construct more great unlined storage caverns, when the rock mass quality is good, cavern shape and construction method also are proper, it is quite possible to construct underground oil storage cavern with span more than 30m . The main axis orientation of Jinzhou underground oil storage cavern is NW direction, the cavern's elevation locate between -53msl and -76msl. The storage's total volume is about 3×106m3, composed of 8 parallel galleries with 950m length, the pillars between them are 45m, and every two of galleries form one unit, which can deposit 75×104m3 for each unit. The product will be stored are Saudi light and Saudi medium crude oil, the main cavern's section is 411.5m2, with 23m height and 19m width. According to the principle and technique of engineering geomechanics, this study supply a sort of system scientific thinking and method for sitting location of underground oil storage in granite region: ① On the foundation of the earth crust stability sub-zone appraise of Circum-Bo sea region, farther research concerning granite distribution, genesis, geological period and fault structure are conducted in stable areas, generally, this paper select Liaoxi, east shore of Liaotung peninsula and Jiaotung peninsula as target areas for underground oil storage regions, where Mesozoic granite is magnitude; ②After roundly comparison in facts of geologic structure, engineering geology, hydrogeology, topography, transportation and so on of three granite distributed areas, at last, selecting Jingzhou granite zone in Liaoxi out as an ideal construction area; ③ Detailed investigation is conducted in the southeast of Baimashi in Jingzhou development district, the final field. Ultrasonic Borehole Television, as a major way to collect original information of borehole rock mass were used, which is very effective to appraise the quality of deep rock mass; ④ According to the field data of tectonic stress, rock mass quality, the spatial distribution of fracture water, some optimum designs in cross section, axial direction and cavern span have been designed for the underground oil storage cavern layout in Jinzhou. To understand the characteristics of swelling alteration rock in Jinzhou granite mass, collected abundant swelling alteration rock engineering examples in granite, which study them in detail, concluded the swelling alteration rock distribute nearly everywhere in China, intruded medium-basic dykes alteration, along discontinuities and mineral hydrothermal alteration with genesis of granite are three main forms clayization alteration rock in granite rock mass. In Jinzhou field, from macro to micro studied the swelling rock which induced by mid-basic dyke intrusion, with weak swelling. In conclusion, this paper conclude the distribution rule and features of expansion alteration rock in filed, and advise some technical suggestions for excavation at swelling alteration rock part. The main features of this paper: ①In the process of site selection, investigation and design, a suit of technique and method of engineering geomechanics metasynthesis were formed, which is significative to guide the large scale underground oil storage cavern sitting location, investigation and design in granite rock mass; ②The detailed discussion on the engineering geology problems in granite mass, such as weathering crust, faults, dykes and clayization alteration rock, are useful for other projects in aspects of site selection, engineering geology evaluation and stability estimation; ③The summary and integration of the genesis, type, countermeasure relate to swelling alteration rock, also is likely to be used for other underground oil storage caverns constructed in swelling alteration granite. In conclusion, this study is meaningful for guiding the large scale underground oil storage for site selection, investigation and design in granite rock mass.

基于GIS的滑坡危险度区划定量研究

Abstract With the Development of the West Regions, the contradiction between economy and geological hazards was once again brought to our face in the Loess Plateau where was ,and is liable to geological hazards for the unique conditions of geology, hydrogeology, geography and meteorology. The goal to realize harmonious development between human and the earth was always there, and landslide hazard zoning provided us an effective way against geological hazards and damage. In the background of the construction of 750KV transformer substation in TianShui, we summarized some theories, methods and development of landslide hazard zoning and discussed the application of information value model in landslide hazard zoning. A called “judgement matrix” like in AHP was introduced to the information value model to solve the key point of landslide hazard zoning — choice of factors and weight of each factor. GIS was applied in the landslide hazard zoning, with its comprehensive function on data management, spatial analysis and mapping. A zonation map of landslide hazard was worked out on MAPGIS aimed to have something of reference and instruction on the construction of the transformer substation. Key words: landslide; hazard zoning; information value model; GIS; judgement matrix

流体地球化学：成因、流动和流体-岩石相互作用及塔里木盆地实例研究

Geofluid in sedimentary basins is related to petroleum generation, migration, accumulation and preservation, and is a topic of geological frontier. By integrating the multi-discipline methods of petroleum geochemistry, sedimentology, hydrogeology, petroleum geology and experimental geochemistry, the thesis has carried out experiments of microcline dissolution in solutions with organic acids, crude oil, brines with high total dissolved solids (TDS), and has dealt with Al distribution between the crude oil and the brines after the experiments. Cases for study includes Central Tarim, Hetianhe Gas Field and Kucha forland basin with data containing fluid chemistry and isotopic compositions, thin sections of sandstones and carbonates, homogenization temperatures and salinities of fluid inclusions, isotopic compositions of bulk rock and autigenic minerals. The aims are to elucidate fluid origin and flow in the three areas, effect of hydrocarbon emplacement on diagenesis, and to show occurrence of microbe-mediated, and thermochemical sulfate reduction in the Tarim Basin. Microcline dissolution experiments show that after 100 hour, part of the dissolved Al distributes in the crude oil, and the Al concentrations in the crude oil rise when organic acids are added. The result can be used to explain that most oilfield waters in the Tarim Basin are characterized by less than 3mg/L Al. Crude oil added to the solutions can enhance microcline dissolution, which is also observed in the case - Silurian sandstones with early crude oil emplacement in the Central Tarim. Al and Si have higher concentrations in the experiments of oxalic acid than of acetic acid under the same pH conditions, suggesting that there exist Al-oxalate and Si-oxalate complexes. Presence of acetate can enhance the activity of Ca and Al, but Al concentrations have not been increased significantly due to formation of small Al-acetate complex during the experiments. Relationships between δD and δ~(18)O in conjunction with chemistry of oilfield waters show that the waters are evaporated connate waters, which subsequently mixed with meteoric water, and were influenced by water-rock interactions such as salt dissolution, dolomitization of calcite, albitization of feldspar. In the Hetianhe Gas Field where salt dissolution took place, δD and δ~(18)O values can be used to trace nicely meteoric water recharge area and flow direction, but TDS can not. Part of the waters have high TDS but very light δD and δ~(18)O. When combined with paleo-topography, or fluid potentials, meteoric water is suggested to flow eastward in the Hetianhe Gas Field, which is the same with the Central Tarim. Whist in the Kuche forland basin, meteoric water may have permeated Cambrian-Ordovician strata. Relationship between ~(87)Sr/~(86)Sr and 1/Sr can be used to indicate migration and mixing of brines from carbonate strata (low ~(87)Sr/~(86)Sr ratio but high Sr content), clastic strata (high ~(87)Sr/~(86)Sr ratio but low Sr content) and crystalline basement (high ~(87)Sr/~(86)Sr ratio and heavy δ~(18)O value). Using this approach, it can be found that ~(87)Sr-depleted brine from Ordovician carbonates have migrated up to and mixed with ~(87)Sr-enriched waters from Silurian and Carboniferous sandstones, and that Silurian brines have mixed with meteoric water. In the Kuche forland basin, brines from the Cambrian and Ordovician carbonates have higher ~(87)Sr/~(86)Sr ratios than those from the overlying sandstones, when combined with chemistry, δ~(15)N and ~3He/~4He ratios of the coexisting natural gases, suggesting that the brines were derived from the basement. There exists some debate on the effect of hydrocarbon emplacement on mineral diagenesis. Case-study from Silurian sandstones in the Central Tarim show that quartz has kept overgrowing secondarily when oil saturation was decreased by meteoric water flushing subsequently to hydrocarbon emplacement. Silicon precipitates on the water-wet quartz surface, leading to decreased Si concentration close to the surface. A Si grads can result in Si diffusion, which supplies Si for quartz overgrowth. Hydrocarbon oxidation-sulfate reduction is an important type of organic-inorganic interaction. Not only can it make secondary alteration of hydrocarbons, but generate H_2S and CO_2 gases which can improve reservoir property. Thermochemical sulfate reduction took place at the temperatures more than 125 ℃ to 140 ℃ in the Cambrian-Ordovician carbonates, the products - H_2S and CO_2 gases migrated up to the Silurian, and precipitated as pyrite and calcite, respectively. The pyrite has an average δ~(34)S value close to those of Ordovician seawater and anhydrite, and calcite has δ~(13)C value as low as -21.5‰. In the Hetianhe Gas Field, sulfate reduction bacteria carried by meteoric water flowing eastward may have preferentially depleted ~(12)C of light hydrocarbon gases, and results in heavier δ~(13)C values of the residual hydrocarbon gases and higher molar CO_2 in the natural gases in the west than in the east. Coexisting pyrite has δ~(34)S values as low as -24.9‰.

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

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.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN

溶蚀岩体渗透水力学特性研究

Seepage control in karstic rock masses is one of the most important problems in domestic hydroelectric engineering and mining engineering as well as traffic engineering. At present permeability assessment and leakage analysis of multi-layer karstic rock masses are mainly qualitative, while seldom quantitative. Quantitative analyses of the permeability coefficient and seepage amount are conducted in this report, which will provide a theoretical basis for the study of seepage law and seepage control treatment of karstic rocks. Based on the field measurements in the horizontal grouting galleries of seepage control curtains on the left bank of the Shuibuya Hydropower Project on the Qingjiang river, a hydraulic model is established in this report, and the computation results will provide a scientific basis for optimization of grouting curtain engineering. Following issues are addressed in the report. (1) Based on the in-situ measurements of fissures and karstic cavities in grouting galleries, the characteristics of karstic rock mass is analyzed, and a stochastic structural model of karstic rock masses is set up, which will provide the basis for calculation of the permeability and leakage amount of karstic rock mass. (2) According to the distribution of the measured joints in the grouting galleries and the stochastic results obtained from the stochastic structural model of karstic rock mass between grouting galleries, a formula for computation of permeability tensor of fracturing system is set up, and a computation program is made with Visual Basic language. The computation results will be helpful for zoning of fissured rock masses and calculation of seepage amount as well as optimization of seepage control curtains. (3) Fractal theory is used to describe quantitatively the roughness of conduit walls of karstic systems and the sinuosity of karstic conduits. It is proposed that the roughness coefficient of kastic caves can be expressed by both fractal dimension Ds and Dr that represent respectively the extension sinuosity of karstic caves and the roughness of the conduit walls. The existing formula for calculating the seepage amount of karstic conduits is revised and programmed. The seepage amount of rock masses in the measured grouting galleries is estimated under the condition that no seepage control measures are taken before reservoir impoundment, and the results will be helpful for design and construction optimization of seepage curtains of the Shuibuya hydropower project. This report is one part of the subject "Karstic hydrogeology and the structural model and seepage hydraulics of karstic rock masses", a sub-program of "Study on seepage hydraulics of multi-layer karstic rock masses and its application in seepage control curtain engineering", which is financially supported by the Hubei Provincial key science and technology programme.

三峡库区滑坡的时空分布、成因及其预报模式

Landslides are widely distributed along the main stream banks of the Three Gorges Reservoir area. Especially with the acceleration of the human economic activities in the recent 30 years, the occurrence of landslide hazards in the local area trends to be more serious. Because of the special geological, topographic and climatic conditions of the Three Gorges areas, many Paleo-landslides are found along the gentle slope terrain of the population relocation sites. Under the natural condition, the Paleo-landslides usually keep stable. The Paleo-landslides might revive while they are influenced under the strong rainfall, water storage and migration engineering disturbance. Therefore, the prediction and prevention of landslide hazards have become the important problem involving with the safety of migration engineering of the Three Gorges Reservoir area.The past research on the landslides of the Three Gorges area is mainly concentrated on the stability analysis of individual landslide, and importance was little attached to the knowledge on the geological environment background of the formation of regional landslides. So, the relationship between distribution and evolution of landslides and globe dynamic processes was very scarce in the past research. With further study, it becomes difficult to explain the reasons for the magnitude and frequency of major geological hazards in terms of single endogenic or exogenic processes. It is possible to resolve the causes of major landslides in the Three Gorges area through the systematic research of regional tectonics and river evolution history.In present paper, based on the view of coupling of earth's endogenic and exogenic processes, the author researches the temporal and spacial distribution and formation evolution of major landslides(Volume^lOOX 104m3) in the Three Gorges Reservoir area through integration of first-hand sources statistics, .geological evolution history, isotope dating and numerical simulation method etc. And considering the main formation factors of landslides (topography, geology and rainfall condition), the author discusses the occurrence probability and prediction model of rainfall induced landslides.The distribution and magnitude of Paleo-landslides in the Three Gorges area is mainly controlled by lithology, geological structure, bank slope shape and geostress field etc. The major Paleo-landslides are concentrated on the periods 2.7-15.0 X 104aB.R, which conrresponds to the warm and wettest Paleoclimate stages. In the same time, the Three Gorges area experiences with the quickest crust uplift phase since 15.0X 104aB.P. It is indicated that the dynamic factor of polyphase major Paleo-landslides is the coupling processes of neotectonic movement and Quaternary climate changes. Based on the numerical simulation results of the formation evolution of Baota landslide, the quick crust uplift makes the deep river incision and the geostress relief causes the rock body of banks flexible. Under the strong rainfall condition, the pore-water pressure resulted from rain penetration and high flood level can have the shear strength of weak structural plane decrease to a great degree. Therefore, the bank slope is easy to slide at the slope bottom where shear stress concentrates. Finally, it forms the composite draught-traction type landslide of dip stratified rocks.The susceptibility idea for the rainfall induced landslide is put forward in this paper and the degree of susceptibility is graded in terms of the topography and geological conditions of landslides. Base on the integration with geological environment factors and rainfall condition, the author gives a new probabilistic prediction model for rainfall induced landslides. As an example from Chongqing City of the Three Gorges area, selecting the 5 factors of topography, lithology combination, slope shape, rock structure and hydrogeology and 21 kinds of status as prediction variables, the susceptibility zonation is carried out by information methods. The prediction criterion of landslides is established by two factors: the maximum 24 hour rainfall and the antecedent effective precipitation of 15 days. The new prediction model is possible to actualize the real-time regional landslide prediction and improve accuracy of landslide forecast.