复杂应力环境下断续节理岩体真三轴实验研究及其工程应用

The great deal of joints and faults , existing in the rock mass , are the leading cause of discontinuous rock mass. Structural planes not only destroy the integrality of rock mass, but also lead nonlinearity、heterogeneity、anisotropy and failure mode on mechanical properties of rock mass. Therefore the selection of strength and deformation parameters was very difficult. In practical rock mass engineering, equivalent parameters of rock mass were selected by the method of expert experience and engineering analogy. Based on the fine description of discontinuous joints in the type Ⅳ and Ⅴ rock mass and geological survey datum in situ, models was obtained by generalizing the structure of rock mass by the method of statistical analysis. Model intensity and deformation test were carried out on the true triaxial apparatus. Intermediate principle stress effect, anisotropy and dimension effect of discontinuous rock mass were considered in the model test. 3-D correction to Hoek-Brown empirical criterion was done by analysed the test datum. Detailed works were listed as follows: (1) The factors influenced intensity and deformation of discontinuous joints rock mass were the value of 、continuity, density and included angle of joints and anisotropy of joint plane. True triaxial intensity and deformation tests were carried out by considering above factors. The influence rule was obtained and corresponding relation formulary was established； (2) Based on the true triaxial tests under different stress path and load modes, we obtain intensity and deformation rule of rock mass； (3) Based on a great deal of true triaxial tests and other test datum, correction to the Hoek-Brown empirical criterion was done in the chapter 4. The intermediate principle stress was considered in the corrected formulary. It indicated that the formulary was applicable under a certain condition. In addition, the yield plane form of corrected Hoek-Brown empirical criterion under principle stress space was described in the paper. And the question of corner of yield plane was discussed； (4) Based on the single discontinuity theory, the three-dimensional intensity formulary of discontinuous joint rock mass was established. Correction to the intensity formulary was done considering intermediate principle stress effect. We may obtain the conclusion that the intensity of the discontinuous joint rock mass was influenced on compositive factors. They were 、 、continuity、internal frictional angle and cohesiveness of joint plane and rock； (5) The results of the true triaxial model test was applied into parameters evaluation of dam foundation rock mass of JinPing hydropower station. For there were abundant ophicalcite in the dam foundation, the interval of intensity and formation parameters influenced on continuity were determined based on test datum. (6) Especial mould for prismatic jointing model was designed. True triaxial intensity and deformation tests by Basalt with prismatic jointing were carried out. The influence of intermediate principle stress, stress path, anisotropy effect and dimensional effect to intensity and deformation was discussed in the chapter 6. The work of (3)、(4)、(6) was significative supplement and innovation to current test and theory.

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

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.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

青藏高原东部电性结构特征及其构造意义

Eastern Himalayan Syntaxis (EHS) and its surroundings (eastern margin of Tibet) is one of the most complicated tectonic areas in the world. As the exhaust opening of the balanced materials of the Tibetan Plateau during the collision of Indan and Eurasian plates, the deep structure beneath EHS surrounding region is referred to as the key to the study of the dynamics of the plateau. EHS3D project, sponsored by NSFC, has been proposed to explore the deep electric features of the area. During the first stage of EHS3D(2006-2008), MT+LMT measurements have been conducted along two lines from Chayu to Qingshuihe (EHS3D-3) and Chayu to Ruoergai (EHS3D-2). This paper will discuss the MT models of EHS3D-3 line. By the data procrssing, including distortion analysis, Robust estimation and strike decomposition, rotated apparent resitivities and phases have been obtained for each station. Then conventional 2-D inversion algorithms (NLCG and RRI) were employed to produce 2-D models. The final preferred 2-D model suggests that the upper crust consists of resistive blocks while in mid-lower crust there are two extensive conductive bodies beneath Lhasa block and Qiangtang terrain respectively. Jinshajiang suture is a gradient belt and Bangong-Nujiang suture appear a conductive belt dipping to the north. . We concluded that the formation of the two conductive bodies attributed to the partial melt and fluids in the lower crust. The regional electric strike derived from decomposition analysis indicates that the crust and upper mantle move in different manners. The upper crust moves like slips of rigid blocks along major slip faults while the lower crust creeps as a flow in the conductive channels.

地壳形变的弹性块体位错模型

The velocity field is important to investigate the motion and strain parameters of the block. It is also important to investigate the deformation of the fault, for example, the accumulation of strain and stress, at the boundary of the block. The dislocation model is a classic method to simulate the velocity field. In dislocation model, the aseismic crustal deformation is regarded as the sum of the rigid block motion and the effect of the locked fault. We modify the dislocation model in two aspects. Firstly, the block motion is assumed to be the sum of rotation and linear strain rather than the rigid motion. Secondly, the elastic layered-earth model rather than the homogenous half-space model is applied to calculate the effect of the locked part. The 1990~1995 annually Global Position System (GPS) velocity data of the Taiwan area are used in our dislocation model. The misfit of our modified model is smaller than that of the origin model clearly. Our simulation shows, in eastern Coastal Range, the velocity decreases northward rapidly from Chimei Fault, which may result from the high crustal compressive rate of about 30 mm•a-1 at Chimei Fault. The lock of fault in southern part is stronger than that in northern part generally. In western Taiwan, the most strongly locked faults appear in the southern Coastal Plain where many disaster earthquakes occur frequently. The calculated strain and rotation rates consist with previous results in most areas. The strain rate field reveals the nearly NW-SE compression in most parts of Taiwan with a fan-shaped distribution. The rotation rate field reveals anticlockwise rotation in eastern and southern Taiwan while clockwise rotation in western and northern Taiwan, generally.

山东莒南地区韧性剪切带的变形特征与形成演化

Ju Nan of Shandong province is located at southwest of Sulu UHP (ultrahigh-pressure) metamorphic terrane. It is composed of gneiss, paragneiss, eclogites, ultramafic rocks, marble and quartzite. A large ductile shear zone extends east-west has been found at the Zhubian, The south of Junan county. The Zhubian ductile shear zone is composed of high srain rock and mylonites. The mylonites fall into 3 types: Initial gneiss mylonite, mylonite and altramylonit.obvious lineation of penetration,foliation,S-Cfabrics,porphyroclasts,folds,irregularundulatory,extinction,subgrain boundary, dynamic recrystallization microstructure, core-mantle structure and are common in the ductile shear zone. Based on field work and microstructural analyse, a conclution is arrived: The ductile shear zone is an approximately SE trending faults. The Zhubian ductile shear zone formed at Ep ―Hb facies conditions which could be proved by deformaed and metamorphosed mineral aggregates, Deformation behavior, Ternary-feldspar geothermometry and so on. Zircon MC―ICP―MS U-Pb analysis is performed on the mylonite and have an average age ―835.9±13.9Ma, it’s the primary rocks formed age. The Zhubian ductile shear zone maybe formed at 224-242Ma.

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

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.

复杂地质条件下高精度地震勘探技术研究与应用-以泌阳凹陷为例

In order to carry out high-precision three-dimensional "integration" for the characteristics of the secondary seismic exploration for Biyang Depression, in the implementation process, through a combination of scientific research and production, summed up high-precision seismic acquisition, processing and interpretation technologies suitable for the eastern part of the old liberated areas, achieved the following results: 1. high-precision complex three-dimensional seismic exploration technology series suitable for shallow depression Biyang block group. To highlight the shallow seismic signal, apply goal-based observing system design, trail from the small panel to receive and protect the shallow treatment of a range of technologies; to explain the use of three-dimensional visualization and coherent combination of full-body three-dimensional fine interpretation identification of the 50-100 m below the unconformity surface and its formation of about 10 meters of the distribution of small faults and improve the small block and stratigraphic unconformity traps recognition. 2. high-precision series of three-dimensional seismic exploration technology suitable for deep depression Biyang low signal to noise ratio of information. Binding model using forward and lighting technology, wide-angle observation system covering the design, multiple suppression and raise the energy of deep seismic reflection processing and interpretation of detailed, comprehensive reservoir description, such as research and technology, identified a number of different types of traps. 3. high-precision seismic exploration technology series for the southern Biyang Depression high steep three-dimensional structure. The use of new technology of seismic wave scattering theory and high-precision velocity model based on pre-stack time migration and depth migration imaging of seismic data and other high-precision processing technology, in order to identify the southern steep slope of the local structure prediction and analysis of sandstone bedrock surface patterns provide a wealth of information.

北山及邻区古生代-早中生代增生与碰撞大地构造格局

The Beishan orogenic collage locates at the triple-joint among Xinjiang, Gansu, and Inner Mongolia Provinces, at which the Siberian, Tarim and North China plates join together. It also occupies the central segment of the southern Central Asian Orogenic Belt (CAOB). The main study area in the present suty focused on the southwest part of the Beishan Mountain, which can be subdivided into four units southernward, the Mazhongshan continental block, Huaniushan Arc, Liuyuan suture zone and Shibanshan-Daqishan Arc. 1. The Huaniushan Arc was formed by northernward dipping subduction from the Orcovician to Permian, in which volcanic rocks ranging from basic to acidic with island arc affinity were widely developed. The granitiod intrusions become smaller and younger southward, whichs indicates a southward rollback of slab. The granitiod intrusions are mainly composed of I type granites, and their geochemical compositions suggest that they have affinities of island arc settings. In the early Paleozoic（440Ma-390Ma）. The Shibanshan-Daqishan Arc, however, were produced in the southernward dipping subduction system from Carboniferous to Permian. Volcanic rocks from basic to acidic rocks are typical calcic-alkaline rocks. The granitiod intrusions become smaller and younger northernward, indicating subdution with a northernward rollback. The granitiod intrusions mainly consist of I-type granites, of which geochemical data support they belong to island arc granite. 2. Two series of adakite intrusions and eruptive rocks have been discovered in the southern margin of the Huaniushan Island Arc. The older series formed during Silurian (441.7±2.5Ma) are gneiss granitoid. These adakite granites intruded the early Paleozoic Liuyuan accretionary complex, and have the same age as most of the granite intrusions in the Huanniushan Arc. Their geochemical compostions demonstrate that they were derived from partial melting of the subudcted oceanic slab. These characteristics indicate a young oceanic crust subduction in the early Paleozoic. The late stage adakites with compositons of dacites associate with Nb-enriched basalts, and island arc basalts and dacites. Their geochemistries demonstrate that the adakites are the products of subducted slab melts, whereas the Nb-enriched basalt is products of the mantle wedge which have metasomatized by adakite melts. Such a association indicates the existences of a young ocean slab subduction. 3. The Liuyuan suture zone is composed of late Paleozoic ophiolites and two series of accretionary complexes with age of early Paleozoic. The early Paleozoic accretionary complex extensively intruded by early Palozioc granites is composed of metamorphic clastics, marble, flysch, various metamorphic igneous rocks (ultramafic, mafic and dacite), and eclogite blocks, which are connected by faults. The original compositions of the rocks in this complex are highly varied, including MORB, E-MORB, arc rocks. Geochronological study indicates that they were formed during the Silurian (420.9±2.5Ma and 421.1±4.3Ma). Large-scale granitiods intruded in the accretionary complex suggest a fast growth effect at the south margin of the Huaniushan arc. During late Paleozoic, island arc were developed on this accretionary complex. The late Paleozoic ophiolite has an age of early Permian (285.7±2.2Ma), in which the rock assemblage includes ultra-mafic, gabbros, gabbros veins, massive basalts, pillow basalt, basaltic clastic breccias, and thin layer tuff, with chert on the top.These igneous rocks have both arc and MORB affinities, indicating their belonging to SSZ type ophiolite. Therefore, oceanic basins area were still existed in the Liuyuan area in the early Permian. 4. The mafic-ultramafic complexes are distributed along major faults, and composed of zoned cumulate rocks, in which peridotites are surrounded by pyroxenite, hornblendites, gabbros norite and diorite outward. They have island-arc affinities and are consistent with typical Alaska-type mafic-ultramafic complexes. The geochronological results indicate that they were formed in the early Permian. 5. The Liuyuan A-type granite were formed under post-collisional settings during the late Triassic (230.9±2.5Ma), indicating the persistence of orogenic process till the late Triassic in the study area. Geochronological results suggested that A-type granites become younger southward from the Wulungu A-type granite belt to Liuyuan A-type granite belt, which is in good agreement with the accretionary direction of the CAOB in this area, which indicate that the Liuyuan suture is the final sture of the Paleo-Asin Ocean. 6. Structural geological evidence demonstrate the W-E spreading of main tectonic terrenes. These terrenes had mainly underwent through S-N direction contraction and NE strike-faulting. The study area had experienced a S-N direction compression after the Permian, indicating a collisional event after the Permian. Based on the evidene from sedimentary geology, paleontology, and geomagnetism, our studies indicate that the orogenic process can be subdivided into five stages: (1) the pre-orogenic stage occurred before the Ordovicain; (2) the subduction orogenic stage occurred from the Orcovician to the Permian; (3) the collisional orogenic stage occurred from the late Permian to the late Triassic; (4) the post-collision stage occurred after the Triassic. The Liuyuan areas have a long and complex tectonic evolutional history, and the Liuyuan suture zone is one of the most important sutures. It is the finally suture zone of the paleo-Asian ocean in the Beishan area.

东天山后碰撞镁铁-超镁铁岩镍铜硫化物矿床PGE地球化学示踪与富集规律研究

Eastern Tianshan area, a Paleozoic complex trench-arc-basin system, experienced multi-period sudbuction and collision in geological history. A large number of Early Permian mafic-ultramafic intrusions emplaced along deep faults in post-collision extension tectonic stage and hosted a series of magmatic Cu-Ni sulfide deposits. This paper sets newly-discovered Tianyu magmatic Cu-Ni sulfide deposit related to small intrusion as an example. On basis of the study of ore-forming mechanism of Cu-Ni deposit, we compared PGE content and discussed enrichment mechanism and regularity of mafic-ultramafic rocks and ores in Jueluotage tectonic belt and Central Tianshan Massif. PGE and Cu, Ni, S contents correlate with each other. PGE is mainly controlled by S content.Samilar PGE distribution patterns of mafic-ultramafic show that complex originated from the same parental magma; Pd/Ir and Ni/Cu ratios indicate that high-Mg basaltic magma and deep sulfide segregation happened during magma evolution process. PGE and Cu-Ni ores are enriched in liquid sulfide and only individual samples completely control by monosulfide solid solution. Comparison of all control factors, early segregation of sulfide and quality of supply of magma may be the key factors leading to the Eastern Tianshan magmatic copper-nickel sulfide deposits don’t reach PGE grade, but we do not rule out the possibility of occurrence continuous mafic-ultramafic rocks and PGE-rich layer in deep.

准噶尔盆地南缘台13井区侏罗系油藏表征与建模

Junggar Basin has a large amount of recoverable reserves, However, due to the unfavorable factors, such as bad seismic data quality, complex structure with many faults and less wells, the exploration of oil and gas is still relatively limited, so advanced theoretical guidance and effective technical supports are desirable. Based on the theories of sedimentology, as well as comprehensive studies of outcrops, seismic data, drilling data and setting of this area, the paper establishes the isochronous correlation framework, and analyzes the sedimentary facies types and provenance direction, and obtains the profile and plain maps of the sedimentary facies combined with the logging constrained inversion. Then the paper analyzes the reservoir controlling factors, reservoir lithology attribute, 4-property relationship and sensibility based on the sedimentary facies research, and sets up a 3D geological model using facies controlled modeling. Finally, the paper optimizes some target areas with the conclusions of reservoir, structure and reservoir formation.Firstly, the paper establishs the isochronous correlation framework by the seismic data, drilling data and setting of this area. The sedimentary facies in Tai13 well block are braided river and meandering river according to the analysis of the lithology attribute, logging facies and sedimentary structure attribute of outcrop. The concept of “wetland” is put forward for the first time. The provenance direction of Badaowan and Qigu formation is obtained by the geology setting, sedimentary setting and paleocurrent direction. The paper obtains the profile and plain maps of the sedimentary facies from the sand value of the wells and the sand thickness maps from the logging constrained inversion. Then, this paper takes characteristics and control factors of the Jurassic reservoirs analysis on thin section observation, scanning transmission electron microscope observation and find out the petrology characteristics of reservoir, space types of reservoir and lithofacies division. In this area, primary pores dominate in the reservoir pores, which believed that sedimentation played the most important roles of the reservoir quality and diagenesis is the minor factor influencing secondary porosity. Using stochastic modeling technique，the paper builds quantitative 3-D reservoir Parameter. Finally, combined the study of structure and reservoir formation, the reservoir distribution regularity is concluded: (a) structures control the reservoir formation and accumulation. (b) Locating in the favorable sedimentary facies belt. And the area which meets these conditions mentioned above is a good destination for exploration.

松辽盆地南部大情字井地区CO2埋存体地质格架及储层建模

Increasing attentions have been paid to the subsurface geological storage for CO2 in view of the huge storage capacity of subsurface reservoirs. The basic requirement for subsurface CO2 storage is that the CO2 should be sequestrated as supercritical fluids (physical trapping), which may also interact with ambient reservoir rocks and formation waters, forming new minerals (chemical trapping). In order to the effective, durable and safe storage for CO2, enough storage space and stable sealing caprock with strong sealing capacity are necessitated, in an appropriate geological framework. Up till now, hydrocarbon reservoirs are to the most valid and appropriate CO2 storage container, which is well proven as the favorable compartment with huge storage capacity and sealing condition. The thesis focuses on two principal issues related to the storage and sealing capacity of storage compartment for the Qingshankou and Yaojia formations in the Daqingzijing block, Southern Songliao Basin, which was selected as the pilot well site for CO2-EOR storage. In the operation area, three facies, including deltaic plain, deltaic front and subdeep-deep lake facies associations, are recognized, in which 11 subfacies such as subaqueous distributary channel, river- mouth bar, interdistributary bay, sheet sandbody, crevasse splay and overflooding plain are further identified. These subfacies are the basic genetic units in the reservoir and sealing rocks. These facies further comprise the retrogradational and progradational depositional cycles, which were formed base- level rise and fall, respectively. During the regressive or lake lowstand stage, various sands including some turbidites and fans occurred mostly at the bottom of the hinged slope. During the progradation stage, these sands became smaller in size and episodically stepped backwards upon the slope, with greatly expanded and deeped lake. However, most of Cretaceous strata in the study area, localized in the basin centre under this stage, are mainly composed of grey or grizzly siltstones and grey or dark grey mudstones intercalated with minor fine sandstones and purple mudstones. On the base of borehole and core data, these siltstones are widespread, thin from 10 to 50 m thick, good grain sorting, and have relative mature sedimentary structures with graded bedding and cross- lamination or crossbeds such as ripples, which reflect strong hydrodynamic causes. Due to late diagenesis, pores are not widespread in the reservoirs, especially the first member of Qingshankou formation. There are two types of pores: primary pore and secondary cores. The primary pores include intergranular pores and micropores, and the secondary pores include emposieus and fracture pores. Throat channels related to pores is also small and the radius of throat in the first, second and third member of Qingshankou formation is only 0.757 μm, 0.802 μm and 0.631 μm respectively. In addition, based on analyzing the probability plot according to frequency of occurrence of porosity and permeability, they appear single- peaked distribution, which reflects strong hetero- geneity. All these facts indicate that the conditions of physical property of reservoirs are not better. One reason may be provided to interpret this question is that physical property of reservoirs in the study area is strong controlled by the depositional microfacies. From the statistics, the average porosity and permeability of microfacies such as subaqueous distributary channel, channel mouth bar, turbidites, is more than 9 percent and 1md respectively. On the contrary, the average porosity and permeability of microfacies including sand sheet, flagstone and crevasse splay are less than 9 percent and 0.2md respectively. Basically, different hydrodynamic environment under different microfacies can decide different physical property. According to the reservoir models of the first member of Qingshankou formation in the No. well Hei47 block, the character of sedimentary according to the facies models is accord to regional disposition evolution. Meantime, the parameter models of physical property of reservoir indicate that low porosity and low permeability reservoirs widespread widely in the study area, but the sand reservoirs located in the channels are better than other places and they are the main sand reservoirs. The distribution and sealing ability of fault- fractures and caprock are the key aspects to evaluate the stable conditions of compartments to store CO2 in the study area. Based on the core observation, the fractures widespread in the study area, especially around the wells, and most of them are located in the first and second member of Qingshankou formation, almost very few in the third member of Qingshankou formation and Yaojia formation instead. In addition, analyzing the sealing ability of eleven faults in the three-dimensional area in the study area demonstrates that most of faults have strong sealing ability, especially in the No. well Hei56 and Qing90-27. To some extent, the sealing ability of faults in the No. well Hei49, Qing4-6 and Qing84-29 are worse than others. Besides, the deposition environment of most of formations in the study area belongs to moderately deep and deep lake facies, which undoubtedly take advantage to caprocks composed of mudstones widespread and large scale under this deposition environment. In the study area, these mudstones distribute widely in the third member of Qingshankou formation, Yaojia and Nenjiang formation. The effective thickness of mudstone is nearly ~550m on an average with few or simple faults and fractures. In addition, there are many reservoir beds with widely- developed insulated interbeds consist of mudstones or silty mudstone, which can be the valid barrier to CO2 upper movement or leakage through diffusion, dispersion and convection. Above all, the closed thick mud caprock with underdeveloped fractures and reservoir beds can be taken regard as the favorable caprocks to provide stable conditions to avoid CO2 leakage.

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

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.

纵波方位属性裂缝检测技术综合研究

Fractured oil and gas reservoir is an important type of oil and gas reservoir, which is taking a growing part of current oil and gas production in the whole world. Thus these technologies targeted at exploration of fractured oil and gas reservoirs are drawing vast attentions. It is difficult to accurately predict the fracture development orientation and intensity in oil and gas exploration. Focused on this problem, this paper systematically conducted series study of seismic data processing and P-wave attributes fracture detection based on the structure of ZX buried mountain, and obtained good results. This paper firstly stimulated the propagation of P-wave in weak anisotropic media caused by vertical aligned cracks, and analyzed the rule of P-wave attributes’ variation associated with observed azimuth, such as travel-time, amplitude and AVO gradient and so on, and quantitatively described the sensitive degree of these attributes to anisotropy of fracture medium. In order to further study the sensitive degree of these attributes to anisotropy of fractures, meanwhile, this paper stimulated P-wave propagation through different types and different intensity anisotropic medium respectively and summarized the rule of these attributes’ variation associated with observed azimuth in different anisotropic medium. The results of these studies provided reliable references for predicting orientation, extensity and size of actual complicated cracked medium by P-wave azimuth attributes responses. In the paper, amounts of seismic data processing methods are used to keep and recover all kinds of attributes applied for fracture detection, which guarantee the high accurate of these attributes, thus then improve the accurate of fracture detection. During seismic data processing, the paper adopted the three dimensional F-Kx-Ky field cone filter technique to attenuate ground roll waves and multiple waves, then enhances the S/N ratio of pre-stack seismic data; comprehensively applying geometrical spread compensation, surface consistent amplitude compensation, residual amplitude compensation to recover amplitude; common azimuth processing method effectively preserves the azimuthal characteristics of P-wave attributes; the technique of bend ray adaptive aperture pre-stack time migration insures to obtain the best image in each azimuth. Application of these processing methods guaranteed these attributes’ accuracy, and then improved the accuracy of fracture detection. After comparing and analyzing a variety of attributes, relative wave impedance (relative amplitude) attribute is selected to inverse the orientation of fracture medium; attenuation gradient and corresponding frequency of 85% energy are selected to inverse the intensity of fracture medium; then obtained the fracture distribution characteristics of lower Paleozoic and Precambrian in ZX ancient buried mountains. The results are good accord with the characteristics of faults system and well information in this area.

藏北可可西里碰撞后钾质火山作用和岩石成因

Post-collisional, potassic igneous rocks are widely distributed in the Hoh Xil area of the northern Tibetan Plateau. Based on the field work, petrography, mineral chemistry, K-Ar geochronology, element and Sr-Nd-Pb isotope geochemistry, this thesis systematically studied the spatial and temporal distribution of the volcanic rocks, chemical characteristics, formation mechanism and partial melting mechanism of the magma source region, geodynamic setting of magmatism, as well as crustal assimilation and fractional crystallization (AFC). The results show that: 1. The Miocene (7.77-17.82 Ma) volcanic products dominantly are trachandesite and trachy, and subordinate rhyolites, associated with stike-slip faults and thrust faults, formed morphology of small lava platforms and cinder cones. 2. Phenocrysts in the lavas are augite, andesine, sanidine, calcic amphibole and subordinate orthopyroxene, biotite and Ti-Fe oxides, displaying typical quench texture. Equilibrium temperatures and pressures of clinopyroxene phenocrysts indicate the magma chamber is located in upper-middle crust. 3. Rhyolites are the products of crustal melting and fractionation of shoshonitic magmas. The source region of intermediate magmas is enriched continental lithospheric mantle, which contains residual minerals such as phlogopite, rutile and spinel, and enriched by subducted sediments during earlier multi-episodes of subduction. 4. Upwelling of asthenosphere provides heat for source region melting, and faults provide channels for magma eruption. 5. Northward underthrusting of Indian continental lithosphere and southward of backstop of Asian continental lithosphere resulted in upwelling of hot asthenosphere. Geochemical characteristics of the potassic magmatism in North Tibet are dominantly controlled by source region composition, partial melting, and crustal assimilation and fractional crystallization (AFC).