中祁连花岗岩年代学、地球化学及其构造意义

The Qilian Orogenic Belts had undergone very complicated evolutional histories and play an important role in understanding the tectonic evolutions of old terrains in northwestern China, in which granitiods formed during Proterozoic-early Mesozoic are widely outcropped. Detailed studies of these granitiods can shed some light on the tectonic evolution of this region. In this thesis, we have conducted geochronological and geochemical studies on eight selected granitic plutons to unravel their emplacement ages and petrogenesis. Furthermore, their tectonic implications were also discussed based on these results. In Neo-Proterozoic, our results suggest that two stages of magmatic activities were taken place in Central Qilian Block, GroupⅠ(750-790Ma) and Group Ⅱ(845－ 930Ma). In Neo-Paleozoic, most granitic plutons were emplaced from Ordovician to Devonian, whereas granitiods with Triassic ages have also been discovered in South Qilian Belt. Inherited zircons with old ages of 1.7Ga, 2.1Ga and 2.7Ga have also been obtained in our study. Geochemical studies suggest that the Proterzoic granites were produced under high pressures and low temperatures from metamorphosed protolith rocks with compostions from basic to intermediate. This implies that some hot sources were underplated beneath lithosophere via mantle-derived magmatism. In combination with regional geological data, we propose that the Cental Qilian block was an old arc terrene during Precambrian, and two stage granitoids were formed under a back-arc extensional setting. Granitic rocks emplaced in early Paleozoic belong to strong peraluminous S-type granites, which were derived from metagreywacke having strong relationships with collisional process. Together with previous data, our results indicate that granitoids in Qilian Orogenic Belt formed during early Paleozoic have different petrogenesis and emplaced ages, which reflect that Qilian Orogenic Belt had underwent complicated multi-stage subduction-collusional processes in early Paleozoic. On the other hand, granitic rocks in South Qilian Belt with Triassic ages were formed by subduction of East Kulun during early Paleozoic-Late Mesozoic, which represent another orogenic episode in the northern margin of Tibetan Plateau.

班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆¬流体-成矿作用

Duobuza copper deposit, newly discovered typical gold-rich porphyry copper deposit with superlarge potential, is located in the Tiegelong Mesozoic tectonic -magmatic arc of the southern edge of Qiangtang block and the northern margin of Bangonghu-Nujiang suture. Quartz diorite porphyrite and grandiorite porphyry, occurred in stock, are the main ore-bearing porphyries. As the emplacement of porphyry stock, a wide range of hydrothermal alteration has developed. Within the framework of the ore district, abundant hydrothermal magnetite developed, and the relationship between precipitation of copper and gold and hydrothermal magnetite seems much close. Correspondingly, a series of veinlets and network veinlets occurred in all alteration zones. Therefore, systematic research on such a superlarge high-grade Duobuza gold-rich porphyry copper deposit can fully revealed the metallogenic characteristics of gold-rich porphyry copper deposits in this region, establish metallogenetic model and prospecting criteria, and has important practical significance on the promotion of regional exploration. In addition, this research on it can enrich metallogenic theory of strong oxidation magma-fluid to gold-rich porphyry copper deposit, and will be helpful to understand the metallogenic characteristics in early of subduction of Gangdese arc stages and its entire evolution history of the Qinghai-Tibet Plateau, the temporal and spatial distribution of ore deposits and their geodynamics settings. Northern ore body of Duobuza copper deposit have been controlled with width (north-south) about 100 ~ 400 m, length (east-west) about 1400 m, dip of 200 °, angle of dip 65 °~ 80 °. And controlled resource amount is of 2.7 million tons Cu with grade 0.94% and 13 tons Au with 0.21g/tAu. Overall features of ore body are large scale, higher grade copper, gold-rich. Ore occurred in the body of granodiotite porphyry and quartz diorite porphyrite and its contact zone with wall rock. Through the detailed mapping and field work studies, some typies of alteration are identificated as follows: albitization, biotititation, sericitization, silication, epidotization, chloritization, carbonatization, illitization, kaolinization and so on. The range of alteration is more than 10km2. Wall alteration zone can be divided into potassic alteration, moderate argillization alteration, argillization, illite-hydromuscovite or propylitization from ore-bearing porphyry center outwards, but phyllic alteration has not well developed and only sericite-quartz veins occurred in local area. Moreover, micro-fracture is development in ore district , and correspondingly a series of veinlets are development as follows: biotite vein (EB type), K-feldspar-biotite-chalcopyrite-quartz vein, magnetite-antinolite-K-feldspar vein, quartz-chalcopyrite-magnetite veins (A-type), quartz-magnetite-biotite-K-feldspar vein, chalcopyrite veinlets in potassic alteration zone; (2) chalcopyrite occurring in the center vein–quartz vein (B type), chalcopyrite veinlets, chalcopyrite-gypsum vein in intermediate argillization alteration; (3) chalcopyrite- pyrite-quartz vein, pyrite-quartz vein, chalcopyrite-gypsum veins, quartz-gypsum- molybdenite-chalcopyrite vein in argillization alteration; (4) gypsum veins, quartz-(molybdenite)-chalcopyrite vein, quartz-pyrite vein, gypsum- chalcopyrite vein, potassium feldspar veinlets, Carbonate veins, quartz-magnetite veins in the wall rock. In short, various veins are very abundant within the framework of the ore district. The results of electronic probe microscopy analysis (EMPA) indicate that Albite (Ab 91.5~99.7%) occurred along the rim of plagioclase phenocryst and fracture, and respresents the earliest stages of alteration. K-feldspar (Or 75.1~96.9%) altered plagioclase phenocryst and matrix or formed secondary potassium feldspar veinlets. Secondary biotite occurred mainly in phenocryst, matrix and veinlets, belong to magnesium-rich biotite formed under the conditions of high-oxidation magma- hydrothermal. Chloritization developed in all alteration zones and alterd iron- magnesium minerals such as biotite and hornblende and then formed chlorite veinlets. As the temperature rises, Si in the tetrahedral site of chlorite decreased, and chlorite component evolved from diabantite to ripiolite. The consistent 280℃~360℃ of formation temperature hinted that chlorite formed on the same temperature range in all alteration zones. However, formation temperature range of chlorite from the gypsum-carbonate-chlorite vein was 190℃~220℃, and it may be the product of the latest stage of hydrothermal activity. The closely relationship between biotite and rutile indicate that most of rutiles are precipitated in the process of biotite alteration and recrystallization. In addition, the V2O3 concentration of rutile from ore body in Duobuza gold-rich porphyry copper deposit is >0.4％, indicate that V concentration in rutile has important significance on marking main ore body of porphyry copper deposit. Apatites from Duobuza deposit all are F-rich. And apatite in the wall rock contained low MnO content and relatively high FeO content, which may due to the basaltic composition of the wall rocks. The MnO in apatite from altered porphyry show a strong positive correlation with FeO. In addition, Cl/F ratio of apatite from wall rock was highest, followed by the potassic alteration zone and potassic alteration zone overprinted by moderate argillization alteration was the lowest. SO2 in Apatite are in the scope of 0 to 0.66%, biotite in the apatite has the highest SO2, followed by the potassic alteration zone, potassic alteration zone overprinted by moderate argillization alteration, and the lowest in the surrounding rocks, which may be caused by the decrease of oxygen fugacity of hydrothermal fluid and S exhaust by sulfide precipitation in potassic alteration. Magnetite in the wall rock have higher Cr2O3 and lower Al2O3 features compared with altered porphyry, this may be due to basalt wall rock generally has high Cr content. And magnetites have higher TiO2 content in potassic alteration than moderate argillization alteration overprinted by potassic alteration, argillization and wall rock, suggested that its formation temperature in potassic alteration was the highest among them. The ore minerals mainly are chalcopyrite and bornite, and Au contents of chalcopyrite, bornite, and pyrite are similar with chalcopyrite slightly higher. The Eu* negative anomaly of disseminated chalcopyrite was relatively lower than chalcopyrite in veinlets. Within a drill hole, the Eu* negative anomaly of disseminated chalcopyrite was gradually larger from bottom to top. Magnetite has the same distribution model, with obvious negative Eu* abnormal, and ΣREE in great changes. The gypsum has the highest ΣREE content and the obvious negative anomaly, and biotite obviously has the Eu* abnormal. Based on the petrographic and geochemical characteristics, five series of magmatic rocks can be broadly classified; they are volcanic rocks of the normal island arc, high-Nb basaltic rocks, adakites, altered porphyry and diorite. The Sr, Nd, Hf isotopes and geochemistry of various series of magmatic rock show that they may be the result of mixing between basic magma and various degrees of acid magma coming from lower crust melted by high temperature basic underplating from partial melting of the subduction sediment melt metasomatic mantle wedge. Furthermore S isotope and Pb isotope of the sulfide, ore-bearing porphyries and volcanic rocks indicated ore-forming source is the mantle wedge metasomatied by subduction sediment melt. Oxygen fugacity of magma estimated by Fe2O3/FeO of whole rock and zircon Ce4+/Ce3+ indicated that the oxidation of basalt-andesitic rocks is higher than ore-forming porphyry, and might imply high-oxidation characteristics of underplated basic magma. Its high oxidative mechanism is likely mantle sources metasomatied by subduction sediment magma, including water and Fe3+. And such high oxidation of basaltic magma is conducive to the mantle of sulfides in the effective access to melt. And the An component of dark part within plagioclase phenocryst zoning belong to bytownite (An 74%), and its may be a result of magma composition changes refreshment by basaltic magma injection. SHRIMP zircon U-Pb and LA-ICP-MS zircon U-Pb geochronology study showed that the intrusions and volcanic rocks from Duobuza porphyry copper deposit belong to early Cretaceous magma series (126~105Ma). The magma evolution series are as follows: the earliest diorite and diorite porphyrite → ore-bearing porphyry and barren grandiorite porphyry →basaltic andesite → diorite porphyrite → andesite → basaltic andesite, and magma component shows a evolution trend from intermediate to intermediate-acid to basic. Based on the field evidences, the formation age of high-Nb basalt may be the latest. The Ar-Ar geochronology of altered secondary biotite, K-feldspar and sericite shows that the main mineralization lasting a interval of about 4 Ma, the duration limit of whole magma-hydrothermal evolution of about 6 Ma, and possibly such a long duration limit may result in the formation of Duobuza super-large copper deposit. Moreover, tectonic diagram and trace element geochemistry of volcanic rocks and diorite from Duobuza porphyry copper deposit confirm that it formed in a continental margin arc environment. Zircon U-Pb age of volcanic rocks and porphyry fall in the range of 105~121Ma, and Duobuza porphyry copper deposit locating in the north of the Bangonghu- Nujiang suture zone, suggested that Neo-Tethys ocean still subducted northward at least early Cretaceous, and its closure time should be later than 105 Ma. Three major inclusion types and ten subtypes are distinguished from quartz phenocrysts and various quartz veins. Vapor generally coexisting with brine inclusions, suggest that fluid boiling may be the main ore-forming mechanism. Raman spectrums of fluid inclusions display that the content of vapor and liquid inclusion mainly contain water, and vapor occasionally contain a little CO2. In addition, the component of liquid inclusions mainly include Cl-, SO42-, Na+, K+, a small amount of Ca2+, F-; and Cl- and Na+ show good correlation. Vapor mainly contains water, a small amount of CO2, CH4 and C2H6 and so on. The daughter minerals identified by Laman spectroscopy and SEM include gypsum, chalcopyrite, halite, sylvite, rutile, potassium feldspar, Fe-Mn-chloride and other minerals, and ore-forming fluid belong to a complex hydrothermal system containing H2O-NaCl-KClFeCl2CaCl2. H and O isotopic analysis of quartz phenocryst, vein quartz, magnetite, chlorite and gypsum from all alteration zones show that the ore-forming fluid of Duobuza gold-rich porphyry copper deposit consisted mainly of magmatic water, without addition of meteric water. Duobuza gold-rich porphyry copper deposit formed by the primary magmatic fluid (600-950C), which has high oxidation, ultra-high salinity and metallogenic element-rich, exsolution direct from the magma, and it is representative of the typical orthomagmatic end member of the porphyry continuum. Moreover, the fluid evolution model of Duobuza gold-rich porphyry copper deposit has been established. Furthermore, two key factors for formation of large Au-rich porphyry copper deposit have been summed up, which are ore-forming fluids earlier separated from magma and high oxidation magma-mineralization fluid system.

荣成基底片麻岩原岩岩石成因及其构造意义

Extensive high to ultrahigh pressure metamorphic rocks are outcropped in the the Dabie-Sulu UHP orogenic belt. Disputes still exist about for protolith nature of metamorphic rocks, petrogenesis, tectonic setting, and influence on upper mantle during the Triassic deep subduction. In this study, a combined study of petrology, geochemistry, isotope geochemistry and zircon chronology was accomplished for high-grade gneisses in the basement of the ultrahigh-pressure metamorphic Rongcheng terrane to reveal protolith nature and petrogenesis of the gneisses and to disucss the magmatic succession along the northern margin of the Yangtze block in Neoproterozoic. Gneisses in the Rongcheng terrane are characterized by negative Nb, Ta, P and Ti anomalies, relatively low Sr/Y ratios and relatively high Ba/La, Ba/Nb and Ba/Zr ratios, mostly displaying geochemical affinity to Phanerozoic volcanic arc. Neoproterozoic protolith ages (0.7 ~ 0.8 Ga) and Paleoproterozoic average crustal residence time (1.92 ~ 2.21 Ga) favour a Yangtze affinity. The gneisses mostly display characteristics of enrichment of LREE, flat heavy rare earth elements (REE) patterns, moderately fractionation between LREE and HREE and slight negative or positive Eu anomalies, probably reflecting that melting took place in the middle to low crust (26 ~ 33 km), where amphibole fractionated from the melts and/or inherited from source material as major mineral phases in the source area. Sr-Nd isotopic composition of the gneisses supports this conclusion. According to εNd(t) and εHf(t) values, the gneisses can be divided into three groups. Gneisses of group I have the highest εNd(t) and εHf(t) values, corresponding to the range of -6 ~ -3 and -2.9 ~ 13.4, respectively. This suggests obvious influx of depleted mantle or juvenile crust in the formation of protoliths. Gneisses of group II have medium εNd(t) (-9 ~ -7) and εHf(t) values (-15.8 ~ -1.4), corresponding to relatively high TDM2(Nd) (1.99 ~ 2.31 Ga) and TDM2(Hf) (1.76 ~ 2.67 Ga) , respectively. This suggests these gneisses were formed by partial melting of Paleoproterozoic crust. Gneisses of group III have the lowest εNd(t) (-15 ~ -10) and εHf(t) values (-15.8 ~ -1.4), corresponding to the largest TDM2(Nd) (1.99 ~ 2.31 Ga) and TDM2(Hf) ( 1.76 ~ 2.67 Ga), respectively. This indicates that gneisses of group III were formed by remelting of Archean crustal material and further demonstrates existence of an Archean basement probably of the Yangtze affinity beneath the Rongcheng terrane. Gneisses of three groups have also certain different geochemical characteristics. Contents of REEs and trace elements reduce gradually from group I to group III. Zirconium saturation temperatures also show similar tendency. Compared to gneisses of group II and group III, gneisses of group I display geochemical feature similar to extensional tectonic setting, having relatively little influence by the source area. Therefore, geochemical characteristics for gneisses of group I can indictate that the protoliths of the Rongcheng gneisses formed in an extensional rifting tectonic setting. This conclusion is supported by the results of eclogites and gabbros previously reported in the Dabie-Sulu orogenic belt. Statistical results of the protolith ages of the Rongcheng gneisses show two age peaks around ~728 Ma and ~783 Ma with an about 50 Ma gap. Extensive magatism in abou 750 Ma along the northern margin of the Yangtze block can hardly be observed in the Rongcheng terrane. This phenomenon likely suggests discontinuous Neoproterozoic magmatism along the northern margin of the Yangtze block.

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

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.

松辽盆地南部大情字井地区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.

急倾斜金属矿体开采岩体移动规律与变形趋势预测研究

Rockmass movement due to mining steep metallic ore body is a considerable question in the surface movement and deformation issue caused by underground mining. Research on coal mining induced rockmass movement and its prediction problem have been performed for a long-term, and have achieved great progress at home and abroad. However, the rockmass movement caused by mining steep metal mine is distinctivly different from coal seam mining.. Existing surface movement laws and deformation prediction methods are not applicable to the rockmass movement caused by mining steep metal mine. So far the home and abroad research to this theory is presently at an early stage, and there isn’t mature theory or practical prediction method, which made a great impact on production. In this paper, the research object—Jinchuan nickel mine, which is typical steep metal mine, characterized by complex geological conditions, developed faults, cracked rockmass, high geostress, and prominent engineering stability problems. In addition, backfill mining method is used in the mine, the features of rockmass movement caused by this mining method are also different from other mining methods. In this paper, the laws of rock mass movement, deformation and destroy mechanism, and its prediction were analyzed based on the collection of data, detailed in-sit engineering geology survey, ground movement monitoring by GPS, theoretical analysis and numerical simulation. According to the GPS monitoring of ground surface movement, ground subsidence basin with apparent asymmetry is developing, the influence scope is larger in the upper faulted block than in the lower faulted block, and the center of ground movement is moving along the upper faulted block direction with increasing depth of mining. During the past half and seven years, the largest settlement has amounted to 1287.5mm, and corresponding horizontal displacement has amounted to 664.6mm. On the ground surface, two fissure belts show a fast-growing trend of closure. To sum up, mining steep metal mine with backfill method also exist the same serious problem of rockmass movement hazards. Fault, as a low intensity zone in rockmass, when it located within the region of mining influence, the change of potential energy mainly consumed in fault deformation associated with rockmass structure surface friction, which is the essence of displacement and stress barrier effects characterized by fault rupture zone. when steep fault located in the tensile deformation region incurred by underground excavation, no matter excavation in hangingwall or in footwall of the fault, there will be additional tensile stress on the vertical fault plane and decrease in the shear strength, and always showing characteristics of normal fault slip, which is the main reason of fault escarpment appeared on the ground surface. The No.14 shaft deformation and failure is triggered by fault activation, which showed with sidewall move, rupture, and break down features as the main form of a concentrated expression of fault effects. The size and orientation of principal stress in surrounding rock changed regularly with mining; therefore, roadway deformation and damage at different stages have different characteristics and distribution models. During the process of mining, low-intensity weak structures surface always showed the most obvious reaction, accompany with surface normal stress decrease and shear strength bring down, to some extent, occurred with relative slide and deformation. Meanwhile, the impact of mining is a relatively long process, making the structure surface effect of roadway deformation and damage more prominent than others under the influence of mining. Roadway surrounding rockmass deformation caused by the change of strain energy density field after excavation mainly belongs to elastic deformation, and the correspondented damage mainly belongs to brittle rupture, in this circumstance, surrounding rockmass will not appear large deformation. The large deformation of surrounding rockmass can only be the deformation associated with structure surface friction or the plastic deformation of itself, which mainly caused by the permanent self-weigh volume force，and long-term effect of mining led to the durability of this deformation Good pitting fill effect and supporting effect of backfill, as well as the friction of rockmass structure surface lead to obvious macro-rockmass movement with long-lag characteristics. In addition, the loss of original intensity and new structure surface arisen increased flexibility in rockmass and fill deformation in structure surface, which made the time required for rockmass potential energy translate into deformation work associated with plastic deformation and structure surface friction consumed much, and to a large extent, eliminated the time needed to do those plastic work during repeated mining, all of which are the fundamental reason of rockmass movement aftereffect more significant than before. Mining steep deposits in high tectonic stress area and in gravity stress area have different movement laws and deformation mechanism. The steep deposit, when the vertical size of the mining areas is smaller than the horizontal size of the orebody, no matter mining in gravity stress area or in high tectonic stress area, they have similar features of ground movement with mining horizontal orebody; contrarily, there will appear double settlement centers on the ground surface under the condition of mining in high tectonic stress area, while there will always be a single center under the other condition. Meanwhile the ground movement lever, scale of mining influence area and macro features of ground movement, deformation and fracture are also different from mining in gravity stress area, and the fundamental reason lies in the impact of orientation of the maximum principal stress on rock movement features in in-site rock stress field. When mining thick and steep deposit, the ground surface movement and deformation characteristic curves are significantly different from excavating the horizontal ore bed and thin steep deposit. According to the features of rockmass movement rate, the development process of mining-induced rockmass movement is divided into three stages: raising stage, steadily stage and gradually decay stage. Considering the actual exploitation situation, GPS monitoring results and macro-characteristics of surface movement, the current subsidence pattern of Jinchuan No.2 mine is in the early stage of development. Based on analysis of surface movement rate, surface subsidence rate increase rapidly when mining in double lever at the same time, and reach its peak until the exploitation model ended. When double lever mining translate into single, production decreased, surface subsidence rate suddenly start to reduce and maintain a relatively low value, and the largest subsidence center will slowly move along with the hangingwall ore body direction with increasing depth of mining, at the same time, the scope and extent of subsidence in footwall ore body will begin magnify, and a sub-settlement center will appear on ground surface, accompanied with the development and closure trend of ground fissure, the surrounding rockmass of shaft and roadway will be confronted to more frequent and severe deformation and failure, and which will have a negative impact on the overall stability of No.2 mine mining. On the premise of continuity of rockmass movement, gray system model can be used in ground rockmass movement prediction for good results. Under the condition of backfill mining step by step, the loose effect of compact status of the hard, broken rockmass led to lower energy release rate, although surrounding rockmass has high elastic energy, loose and damage occurred in the horizontal ore body, which made the mining process safety without any large geological hazards. During the period of mining the horizontal ore body to end, in view of its special “residual support role”, there will be no large scale rockmass movement hazards. Since ground surface movement mainly related to the intensity of mining speed and backfill effect, on the premise of constant mining speed, during the period of mining the horizontal ore body to end, the rate of ground surface rockmass movement and deformation won’t have sudden change.

松辽盆地南部岩性油藏成藏模式与分布规律

According to the basic geologic conditions, the paper is directed by the modem oil-gas accumulation theory and petroleum system in which typical oil pools are analyzed and the shape of lithologic trap and geologic factors are pointed out. The process during which oil and gas migrate from source rock to lithologic trap is rebuilt, and the accumulation model of oil pool is set up. With the comprehensive application of seismic geologic and log data and paying attention to the method and technology which is used to distinguish lithologic accumulation. Promising structural-lithofacies zones are got and the distribution rule of various lithologic accumulation is concluded. With making use of the biologic mark compound, different reservoirs are compared. As a result, the oil and gas in HeiDimiao come from Nenjiang Group's source rocks; in SaErTu from QingShenkou Group's and Nenjiang Group's, and in PuTaohua. GaoTaizi and FuYang from QingShankou Group's. According to the development and distribution of effective source rock, oil distribution and the comparison in the south of SongLiao basin, the characteristic of basin structure and reservoir distribution is considered, and then the middle-upper reservoir of SongLiao basin south are divided into two petroleum system and a complex petroleum system. Because of the characteristic of migration and accumulation, two petroleum systems can furtherly be divided into 6-7 sub-petroleum systems,20 sub-petroleum systems in all. As a result of the difference of the migration characteristic, accumulation conditions and the place in the petroleum system, the accumulation degree and accumulation model are different. So three accumulation mechanism and six basic accumulation model of lithologic trap are concluded. The distribution of lithologic pools is highly regular oil and gas around the generation sag distribute on favorable structural-lithofacies zones, the type of lithological pool vary regularly from the core of sandstone block to the upper zone. On the basic of regional structure and sedimentary evolution, main factors which control the form of trap are discovered, and it is the critical factor method which is used to discern the lithologic trap. After lots of exploration, 700km~2 potential trap is distinguished and 18391.86 * 10~4 tons geologic reserves is calculated. Oil-water distribution rule of pinch-out oil pool is put up on plane which is the reservoirs can be divided into four sections. This paper presented the law of distribution of oil and water in updip pinch-out reservoir, that is, hydrocarbon-bearing formation in plane can be divided into four zones: bottom edge water zone, underside oil and water zone, middle pure oil zone and above residual water zone. The site of the first well should be assigned to be middle or above pure oil zone, thus the exploration value of this type of reservoir can be recognized correctly. In accordance with the characteristics of seism and geology of low permeability thin sandstone and mudstone alternation layer, the paper applied a set of reservoir prediction technology, that is: (1)seism multi-parameter model identification; (2) using stratum's absorbing and depleting information to predict reservoir's abnormal hydrocarbon-bearing range. With the analysis of the residual resource potential and the research of two petroleum system and the accumulation model, promising objective zones are predicted scientifically. And main exploration aim is the DaRngZi bore in the west of ChangLin basin, and YingTai-SiFangZi middle-upper assembly in Honggang terrace.

华南地壳岩石弹性和变形性质的高温高压实验研究

The South continent of China lies to southeast of Eurasia block. It is an active area from the view of crust growth and continent spread and is a transition zone between continental crust and oceanic crust. The compressional wave velocities and anisotropies of typical crustal metamorphic rocks were determined at high temperature (up to 1000 ℃) and high pressure(up to 800MPa). The experimental results show that the velocities generally increase with pressure, and is unaffected by temperature up to around 550 ℃. But the velocities of all experimental samples start to drop above a temperature point. For an example, this greatly reduce the speed of wave propagation in amphibolite and serpentinite above 760 ℃ and above 550 ℃ respectively, which may be due to dehydrate of amphibole and serpentine. P-wave anisotropy coefficients of those rocks range from 2% to 10% at 800MPa and 500 ℃. The anisotropies decrease with increasing pressure at room temperature, but hardly change as function of temperature at constant 800MPa or 600MPa pressure. The average velocity of the six crustal rocks is 6.28km/s under the condition of 800MPa and 550 ℃, which is consistent with the result of deep seismic sounding data. Based on this experimental result, we deduce there may exist a lot of felsic granulites and amphibolites at the depth of 15-25km underground. With increasing temperature and pressure, the deformation behavior of the rocks undergoes from localized brittle fracture, semi-brittle deformation (cataclastic flow or semi-brittle faulting, semi-brittle flow) to homogeneous crystal-plastic flow. This transition is associated with mechanical behavior and micro-mechanism. It is very important to understanding earthquake source mechanics, the strength of the lithosphere and the style of deformation. The experiments were conducted at temperature of 600-1000 ℃, confining pressure of 500MPa, and stain rates of 10~(-4)-10~(-6) S~(-1). For fine-grained natural amphibolite, the results of experiments show that brittle faulting is major failure mode at temperature <600 ℃, but crystal-plastic deformation is dominate at temperature >800 ℃, and there is a transition with increasing temperature from sembrittle faulting to cataclastic flow and sembrittle flow at temperature of 670-750 ℃. For medium-grained natural Felsic granulite, the results of experiments show that brittle faulting is major failure mode at temperature <500 ℃, but crystal-plastic deformation is dominate at temperature >700 ℃, and there is a transition with increasing temperature from semibrittle faulting to cataclastic flow and sembrittle flow at temperature of 500-600 ℃.

铂族元素分析方法研究及其在有关地质问题中的应用

Characterization of Platinum Group Elements (PGE) has been applied to earth, space and environmental sciences. However, all these applications are based on a basic prerequisite, i.e. their concentration or ratio in the research objects can be accurately and precisely determined. In fact, development in these related studies is a great challenge to the analytical chemistry of the PGE because their content in the geological sample (non-mineralized) is often extremely low, range from ppt (10~(-12)g/g) to ppt (10~(-9)g/g). Their distribution is highly heterogeneous, usually concentrating in single particle or phase. Therefore, the accurate determination of these elements remains a problem in analytical chemistry and it obstructs the research on geochemistry of PGE. A great effort has been made in scientific community to reliable determining of very low amounts of PGE, which has been focused on to reduce the level of background in used reagents and to solve probable heterogeneity of PGE in samples. Undoubtedly, the fire-assay method is one of the best ways for solving the heterogeneity, as a large amount of sample weight (10-50g) can be hold. This page is mainly aimed at development of the methodology on separation, concentration and determination of the ultra-trace PGE in the rock and peat samples, and then they are applied to study the trace of PGE in ophiolite suite, in Kudi, West Kunlun and Tunguska explosion in 1908. The achievements of the study are summarized as follows: 1. A PGE lab is established in the Laboratory of Lithosphere Tectonic Evolution, IGG, CAS. 2. A modified method of determination of PGE in geological samples using NiS Fire-Assay with inductively coupled plasma-mass spectrometry (ICP-MS) is set up. The technical improvements are made as following: (1) investigating the level of background in used reagents, and finding the contents of Au, Pt and Pd in carbonyl nickel powder are 30, 0.6 and 0.6ng/g, respectively and 0.35, 7.5 and 6.4ng, respectively in other flux, and the contents of Ru, Rh, Os in whole reagents used are very low (below or near the detection limits of ICP-MS); (2) measuring the recoveries of PGE using different collector (Ni+S) and finding 1.5g of carbonyl nickel is effective for recovering the PGE for 15g samples (recoveries are more than 90%), reducing the inherent blank value due to impurities reagents; (3) direct dissolving nickel button in Teflon bomb and using Te-precipitation, so reducing the loss of PGE during preconcentration process and improving the recoveries of PGE (above 60% for Os and 93.6-106.3% for other PGE, using 2g carbonyl nickel); (4) simplifying the procedure of analyzing Osmium; (5)method detection limits are 8.6, 4.8, 43, 2.4, 82pg/g for 15g sample size ofRu, Rh, Pd, Ir, Pt, respectively. 3. An analytical method is set up to determine the content of ultra-trace PGE in peat samples. The method detection limits are 0.06, 0.1, 0.001, 0.001 and 0.002ng/mL for Ru, Rh, Pd, Ir and Pt, respectively. 4. Distinct anomaly of Pd and Os are firstly found in the peat sampling near the Tunguska explosion site, using the analytical method. 5. Applying the method to the study on the origin of Tunguska explosion and making the following conclusions: (1) these excess elements were likely resulted from the Tunguska Cosmic Body (TCB) explosion of 1908. (2) The Tunguska explosive body was composed of materials (solid components) similar to C1 chondrite, and, most probably, a cometary object, which weighed more than 10~7 tons and had a radius of more than 126 m. 6. The analysis method about ultra-trace PGE in rock samples is successfully used in the study on the characteristic of PGE in Kudi ophiolite suite and the following conclusions are made: (1) The difference of the mantle normalization of PGE patterns between dunite, harzburgite and lherzolite in Kudi indicates that they are residual of multi-stage partial melt of the mantle. Their depletion of Ir at a similar degree probably indicates the existence of an upper mantle depleted Ir. (2) With the evolution of the magma produced by the partial melt of the mantle, strong differentiation has been shown between IPGE and PPGE; and the differentiation from pyroxenite to basalt would have been more and more distinct. (3) The magma forming ophiolite in Kudi probably suffered S-saturation process.

东营凹陷辛镇复杂断块油藏形成机制和剩余油分布规律研究

In this paper, the complex faulted-block oil reservoir of Xinzhen area in Dongying depression is systematically studied from basic conditions forming faulted-block oil and gas reservoir integrating geology, seismic, logging and reservoir engineering information and computer; guided by petroleum geology, geomechanics, structural geology and geophysics and other theories. Based on analysis of background condition such as regional strata, structure and petroleum geology, structural research on geometry, kinemaitcs and dynamics, oil-controlling fault research on the seal features, sealing mechanism and sealing pattern, and research on enrichment rules and controlling factors of complex faulted-block oil reservoir are carried out to give out the formation mechanics of oil reservoir of Xinzhen complex faulted-block oil reservoir. As a result, the reservoir formation pattern is established. At the same time, through dissecting the characteristics and hydrocarbon enrichment law of complex faulted-block oil reservoir, and studying its distribution law of remaining oil after entering extra high water-cut period, a set of technologies are formed to predict complex faulted-block oil reservoir and its remaining oil distribution and to enhance oil recovery (EOR). Based on the time relationship between migration of hydrocarbon and trap formation, accumulating period of Xinzhen oil reservoir is determined. The formation of Xinzhen anticlinal trap was prior to the primary migration. This is favorable to formation of Xinzhen anticlinal hydrocarbon reservoir. Meanwhile, because anticline top caving isn't at the sane time as that of moving or faulted-trap forming inner anticline, oil and gas migrated many times and Xinzhen complex faulted-block oil reservoir formed from ES_3~(upper) to EG. Accumulating law and controlling factors of complex faulted-block reservoir are analyzed from many aspects such as regional structure background controlling hydrocarbon accumulating, plastic arch-open structure controlling oil-bearing series and reservoir types, sealing-opening of fault controlling hydrocarbon distribution and structure pattern controlling enriched trap types. Also, we established the structure pattern in Xinzhen a'ea: the arch-open of underlying strata cause expanding fracture. The main block groups developed here are shovel-like normal fault block group in the north area of Xinzhen and its associated graben block group. Block groups dominate the formation and distribution of reservoirs. We studied qualitatively and quantitatively the sealing characteristics, sealing history and sealing mechanism of faults, too. And, the sealing characteristics are evaluated and the distribution pattern of hydrocarbon controlled by faults is researched. Due to movement intensity of big faults, deep falling of downthrown block, high degree of repture and development of fracture, shallow layers close to the downthrown block of secondary faults are unfavorable to hydrocarbon accumulation. This is confirmed by the exploration practice in Xinzhen anticline. In terms of the downthrown blocks of sencondary contemporaneous faults lied in the south and north area of Xinzhen, hydrocarbon is poor close to fracture belt, while it is relatively abundant in tertiary companion faults. Because of long-term movement of faults that control hydrocarbon, fi'om ES3 to EG, six set of oil-bearing series formed. And their opening causes the inhomogeneity in hydrocarbon abundance among each block--in two flanks of anticline reservoirs are abundant while in the axial area, oil and gas are sporadic. There the sealing characteristics control oil-bearing area of oil/gas accumulation and the height of oil reservoir. Longitudinally, oil and gas are enriched in dip-flat areas in mid-plane of faults. It is established that there are four types of accumulating patterns in complex faulted-block oil reservoirs in Xinzhen. The first is accumulating pattern of lithologic oil reservoirs in E~S_3~(mid-lowwer), that is, self-generating-self-reserving-self-covering lithologic trap pattern. The second is drag-anticline accumulating pattern in Xinzhen. The structure traps are drag anticlines formed by the contemporaneous faults of the second basement in the north of Xinzhen, and the multiple source rocks involve Ek_2, Es_4, Es_3 and Es_1 members. The reservoirs are fluvial-delta sandstones of the upper member of Shahejie formation and Guantao formation, covered by regional thick mudstone of the upper member of Guantao formation and MingHuazhen formation. The third is the accumulating pattern of reverse listric fault, the third-degree fault of Xinzhen anticline limb and the reservoirs form reservoir screened by reverse listric faults. The forth is accumulating pattern of crossing faults which form closing or semi-closing faulted-blocks that accumulate hydrocarbon. The technologies of predicting remaining oil in complex faulted-block reservoir during the mid and late development stage is formed. Remaining oil in simple large faulted-blocks enriches in structural high, structural middle, structural low of thick bottom water reservoirs, points near bent edge-fault oftertiary faults and part the fourth ones with big falling displacement, microstructure high place of oil-sandbodies and areas where local well pattern isn't perfect. While that in small complex faulted-blocks enriches near small nose, small high point, angle of small faults, small oil-bearing faulted-blocks without well and areas with non-perfect well pattern. The technologies of enhancing recovery factor in complex faulted-block reservoir during the mid and late development stage is formed as follows: fine reservoir description, drilling adjust wells, designing directional wells, sub-dividing layer series of development, improving flooding pattern, changing water-injection direction and enhancing swept volume, cyclic waterflooding and gas-injection, etc. Here, directional wells include directional deflecting wells, lateral-drilling wells, lateral-drilling horizontal wells and horizontal wells. The results of this paper have been used in exploration and development of Shengli oilfield, and have achieved great social and economic profit, especially in predicting distribution of complex faulted-block reservoir, remaining oil distribution during middle and late stage of development, and in EOR. Applying the achievement of fault-closure research, new hydrocarbon-bearing blocks are discovered in flanks of Dongying central uplift and in complex blocks with proved reserves 15 million tons. With the study of remaining oil distribution law in complex faulted-block reservoirs, recovery factors are increased greatly in Dongxin, Xianhe and Linpan complex faulted-block reservoirs and accumulated oil production increment is 3 million tons.

地震波散射非线性反演问题不动点理论及方法的研究

How to create a new method to solve the problem or reduce the influence of that the result of the seismic waves scattering nonlinear inversion is not uniqueness is a main purpose of this research work in the paper. On the background of research into the seismic inversion, new progress of the nonlinear inversion is introduced at the first chapter in this paper. Especially, the development, basic theories and assumptions on some major theories of seismic inversion are analyzed, discussed and summarized in mathematics and physics. Also, the problems faced by the mathematical basis of investigations of the seismic inversion are discussed, and inverse questions of strongly seismic scattering due to strong heterogeneous media in the Earth interior are analyzed and viewed. What the kernel of paper is that gathers all our attention making a new nonlinear inversion method of seismic scattering. The paper provides a theory and method of how to introduce the fixed-point theory into the nonlinear seismic scattering inversion and how to obtain the solution, and gives the actually method to create a serials of contractive mappings of velocity parameter's in the mapping space of wave. Therefore, the results testify the existence of fixed point of velocity parameter and give the method the find it. Further, the paper proves the conclusion that the value obtained by taking the fixed point of velocity parameter into wave equation is the fixed point of the wave of the contractive mapping. Thence, the fixed point is the global minima since the stabilities quality of the fixed point. Based on the new theory, in the chapter three, many inverse results are obtained in the numerical value test. By analysis the results one could find a basic facts that all the results, which are inversed by the different initial model, are tended to the true value in theoretical true model. In other words, the new method mostly eliminates the non-uniqueness that which is existed in seismic waves scattering nonlinear inversion in degree. But, since the test results are quite finite now, more test is need here to positive our theory. As a new theoretical method, it must be existed many weaken in it. The chapter four points out all the questions which is bother us. We hope more people to join us to solve the problem together.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

港西油田现代油藏管理研究与实践

Through years of practice, reservoir management has already become the basic mode of foreign oil companies to realize the high-efficient development of the oil field. From the view of reservoir development and technological economy, reservoir management regards the study of the reservoir engineering, designs of reservoir projects and the dynamic analysis of the reservoir's performance as a system. In the fields of reservoir description, the establishment of the geological models and development models, the dynamic simulations of reservoir exploitation and the design of the oil engineering, reservoir management emphasizes the cooperation of the geology and the engineering, the combination of the engineering technology and the economic evaluation. In order to provide the means and basis for the reservoir geology study, reservoir evaluation, reserves calculation, numerical simulation, development plan and risk analysis, it adopts the reservoir management activities(team work) to make and implement the optimized oil field development management strategies so that secientific and democratic decision making can be achieved. Under the planned economic system for a long time, the purpose of Chinese reservoir development has been to fulfill the" mandatory" production task. With the deepening of the reform, the management organization of Chinese petroleum enterprises has been gradually going through the transition and reforms to the operational entity and the establishment of the mode of oil companies under the socialist market economy system. This research aims at introducing the advanced reservoir management technique from foreign countries to further improve the reservoir development results and wholly raise the economic benefits of Chinese mature land facieses sandstone reservoirs in the later stage of the water flooding. We are going to set up a set of modern reservoir management modes according to the reservoir features, current situation and existing problems of GangXi oil field of DaGang oil company. Through the study and implementation of the reservoir description and numerical simulation technology effectively, we plan to work out integrated adjustment projects, to study the related technology of oil recovery; to set up the effective confirmable data procedure and data management system of the reservoir management, to establish the coordinated model and workbench related to geology, engineering and economy in order to realize the real time supervision and evaluation on the process of reservoir development. We hope to stipulate modernization management tools for GangXi oil fields to rationally utilize various kinds of existing technological methods and to realize the economic exploitation and achieve the maximum benefits from the reservoir. The project of the modem reservoir management will be carried out on the GangXi oil field of DaGang oil company for this oil field is typical and has integrated foundamental materials and perfect networks. Besides, it is located in the good geographical position enjoying very convenient traffic. Implementing modern reservoir management will raise the recovery ratio, reduce the production cost and improve the working efficiency. Moreover, the popularization of modern reservoir management will improve the comprehensive benefits of DaGang oil company and even the whole Petro China. Through the reserch of this project, the following technical indicators can be reached: Establishing the concept of modern reservoir management. Establishing a set of integrated data information management system adapt to the features of GangXi reservoir. 3. Forming technical research modes of modern reservoir management suitable for mature reservoirs in the later developing stage. 4. Advancing projects of GangXi reservoir which are maxium optimized in engineering technique and economic benefits of oil exploitation. Besides, this set of technology, research principle and method can guide the mature reservoir of DaGang oil field and even the whole PetroChina to develop the further research of reservoir adjustment and improve the reservoir recovery factor and developing level constantly.

利用宽频带台阵资料研究华北地区地壳及上地幔间断面结构

Intense tectonic renovation has occurred in the eastern continent of china since Mesozoic, as evidenced by the high heat flow, widespread magma extrusion and volcanic activities, and development of large sedimentary basins. To explain the cause and mechanism for the tectonic process in this period, some researchers have put forward various models, such as mantle plume, subduction of the Pacific slab, Yangtze Block-North China Block collision, etc. Their seismological evidence, however, is still scarce..During the period from 2000 to 2003, large temporary seismic arrays were established in North China by the Institute of the Geology and Geophysics, Chinese Academy of Sciences. Total 129 portable seismic stations were linearly emplaced across the western and eastern boundaries of the Bohai Bay Basin, and accumulated a large amount of high-quality data. Moreover, abundant data were also collected at the capital digital seismic network established in the ninth five-year period of national economic and social development. These provide an unprecedented opportunity for us to study the deep structure and associated geodynamic mechanism of lithospheric processes in North China using seismological techniques.Seismology is a kind of observation-based science. The development of seismic observations greatly promotes the improvement of seismologic theory and methodology. At the beginning of this thesis, I review the history of seismic observation progress, and present some routine processing techniques used in the array seismology. I also introduce two popular seismic imaging methods (receiver function method and seismic tomography).Receiver function method has been widely used to study the crustal and upper mantle structures, and many relevant research results have been published. In this thesis I elaborate the theory of this method, including the basic concept of receiver functions and the methodology for data pre-processing, stacking and migration. I also address some problems often encountered in practical applications of receiver function imaging.By using the teleseismic data collected at the temporary seismic arrays in North China, in particular, the traveltime information of P-to-S conversion and multiple reverberations of the Moho discontinuity, I obtain the distributions of the crustal thickness and the poisson ratio at the northwest boundary area of the Bohai Bay Basin and discuss the geological implications of the results.Through detailed intestigations on the crustal structural feature around the middle part of the Tanlu fault, considerable disparity in poisson ratios is found in the western and eastern sides of the Tanlu fault. Moreover, an obvious Moho offset is coincidently observed at the same surface location. A reasonable density model for the Tanlu fault area is also derived by simulating the observed gravity variations. Both receiver function study and gravity anomaly modeling suggest that the crustal difference between the western and eastern sides of the Tanlu fault is mainly resulted from their different compositions.With common conversion point imaging of receiver functions, I estimate the depths of the upper and lower boundaries of the mantle transition zone, i.e., the 410 and 660 km discontinuities, beneath most part of the North China continent The thickness of the transition zone (TTZ) in the study area is calculated by subtracting the depth of .410 km discontinuity from that of the 660km discontinuity. The resultant TTZ is 10-15 km larger in the east than in the west of the study area. Phase transitions at the 410 km and the 660 km discontinuities are known to have different Clapeyron slopes. Therefore, the TTZ is sensitive to the temperature changes in the transition zone. Previous studies have shown that the TTZ would be smaller in the mantle plume areas and become larger when the remnants of subducted slabs are present The hypothesis of mantle plume cannot give a reasonable interpretation to the observed TTZ beneath North China, Instead, the receiver function imaging results favor a dynamic model that correlates the thermal structure of the mantle transition zone and associated upper mantle dynamics of North China to the Pacific plate subduction process.