东天山红山高硫型浅成低温铜-金矿床氧化带的形成机理和硫酸盐矿物的系列发现及意义

Located in the Paleozoic uplift along the southern margin of Tu-Ha basin in eastern Xinjiang, the newly discovered Hongshan Cu-Au deposit occurs in the superimposed Mesozoic volcanic basin upon the north section of later Paleozoic Dananhu-Tousuquan accretionary arc. Kalatage Cu-Au orebelt is controlled by NWW-trend faults, and includes Hongshan and Meiling Cu-Au deposits. The host rocks of Hongshan ore district are mainly rhyolitic-dacitic ignimbrites, whereas Cu-Au mineralization is closely related to quartz porphyry, rhyolitic porphyry and granitic porphyry. Mineralization styles are dominantly veinlet-disseminated and veinlet, occasionally stockwork. The mineral association is chalcopyrite, pyrite, bornite, chalcocite and sphalerite. The hydrothermal alteration consists of silicfication, sericitization, alunitization, pyrophylitization, illitization, hydromuscovitization, and chloritization. Hongshan Cu-Au deposit, on the edge of the desert, is one of the driest areas in eastrn Tianshan. Moreover， the highest temperature has been up to 60℃, and the average rainfall receives only 34.1mm/y. The light rainfall and rapid evaporation in the vicinity of this deposit have allowed the formation of a great variety of water-soluble sulfates. Oxidization zone of this deposit lies on the upper part of primary sulfide orebodies appearing with a depth of 50-60m, which is dominant in sulfate minerals. 1. Based on the field observation, the volcanic and sub-volcanic rock composition, hydrothermal alteration, ore structure and mineralization characteristics, this paper proposed that the Hongshan Cu-Au deposit belongs to a transitional type from high-sulfide epithermal to porphyry Cu-Au deposit, which corresponds with the typical HS-epithermal deposit such as Zijinshan Au-Cu deposit in Fujian Province, SE-China. 2. The Hongshan copper-gold deposit was controlled by the tectonic, stratum, magma activity and volcanic apparatus, whereas Au mineralization is closely related to quartz porphyry, rhyolitic porphyry and fine grained pyritization in hydrothermal activity, and Cu mineralization is closely related to quartz porphyry and hydrothermal explosive breccia. 3. Oxidation zone of Hongshan Cu-Au deposit lies on the upper part of primary sulfide orebodies deposit. 23 sulfate minerals were identified in this work. The results of samples XRD and chemical analysis were furthermore confirmed through thermal, infrared spectrum and mössbauer spectrum analysis. Among those, nine minerals as Ferricopiapite, Cuprocopiapite, Rhomboclase, Parabutlerite, Krausite, Yavapaiite, Metasideronatrite Kroehnkite and Paracoquimbite were founded in China for the first time. And Paracoquimbite was secondly reported in the world (first case reported at 1938 in Chile). 4. EPMA analysis shows that Al impurity in crystal lattice is important to polytype formation of paracoquimbite and coquimbite besides stack fault. 5. Compared with Meiling Cu-Au deposit in the same Kalatage ore belt from the characteristics of δ34S of barite, lithofacies, hydrothermal alteration and homogeneous temperature, Hongshan Cu-Au deposit belongs to the same metallogenic system of HS-epithermal type as Meiling Cu-Au deposit. But Hongshan Cu-Au deposit has less extensive alteration and shallower denudation. 6. Sulfur isotope analyses show that δ34S values of pyrites vary in the range of +1.86‰~+5.69‰, with an average of 3.70‰, mostly in the range of +1.86‰~+3.20‰, and δ34Scp<δ34Spy. Therefore ore-forming fluid of porphyry comes from mantle and was contaminated by the earth’s crust. Sulfur isotope has reached balance in ore-forming process. 7. Sulfur isotope analyses show that δ34S values of sulfates vary in the range of +2.15‰~+6.73‰, with an average of +3.74‰, mostly equals as δ34S values of primary sulfides in Hongshan Cu-Au deposit. So supergene sulfates inherit sulfur of primary sulfide. δ34S values are mostly same in different sulfates. As well as pyrite and chalcopyrite, volcanic hot spring and associated native sulfur underground also provide water medium and sulfur during the formation process of sulfate. 8. According to the EPMA of sample chalcopyrite and pyrite in Hongshan Cu-Au, the value of Cu/Ni is 0.98-34.72, mostly close to the value of 5, which shows that Hongshan deposit is a typical volcanogenic magmaic hypothermal deposit. Au and Ag, Zn, Te and Bi are positive correlation, Cu and Hg, Se, Sb are positive correlation, indicates Au and Cu don’t locate in the factor of mineralization of same mineralization groups. The reasons of gold concentration in the oxidation zone are: 1). Change of redox potential (Eh) makes gold to deposit from the liquid of mineralization zone; 2). PH is one of the most factors of gold’s deposition; 3). Soluble complex and colloid of gold can be adsorbed easily. 9. The biotite and hornblende K-Ar isotopic ages from the wall rock－quartz diorite, biotite granite and monzonite granite are 231.99±3.45Ma, 237.97±2.36Ma and 296.53±6.69Ma respectively. The ore-bearing rhyolitic breccia lava contains breccia of the biotite granite which indicates the volcanism and related Cu-Au mineralization occurred later than the granite, possibly in Mesozoic. K-Ar ages of granitoids in Sanya, Baishiquan and Hongliugou area and Molybdenite Re-Os age of Baishan Mo deposit all are in Triassic. Besides late Paleozoic magmatism, igneous magmatic event of Mesozoic was widespread in eastern Tianshan. 10. The K-Ar age dating indicates that the K-Ar age of Voltaite occurred below surface 1m is 56.02±3.98Ma, K-Ar age of Ferricopiapite occurred below surface 1.5m is 8.62±1.12Ma, K-Ar age of Yavapaiite occurred below surface 14 m is 4.07±0.39Ma, and K-Ar age of Voltaite occurred below surface 10 m is 14.73±1.73Ma. So the age interval of oxidation zone of Hongshan copper-golden bed is between 60 -3.38Ma. Oxidization occurred at Caenozoic era (from 65Ma), which can be identified through comparing with different deposits oxidation zone in other countries. The coupling between global tectonic event and climatic change event which occur from Caenozoic era has some effect on epigeosphere system, which can act on the surface of bed oxidation zone similarly. It induces that the age mentioned above coincide with collision of India-Asia and multistage uplifting of Qinhai-Tibet Plateau happened subsequently. Bed oxidation zone is the effect and record of collision and uplifting of Tibet Plateau. The strong chemical weathering of surface accumulation to which was leaded by PETM event occurred Paleocene and Eocene is the reason of Voltaite sharply rises. On the contrary, Ferricopiapite formed due to the global cold weather. The predecessor did much research through biota, isotopes, susceptibility, but this paper try to use different sulfate mineral instead of climatic change. So the research of sulfate minerals not only indicates a great deal of oxidized zone feature, but also the intergrowth of sulfate minerals may be used to trace paleoenviroment and paleoclimate of oxidation zone. 11. Analysis of the information of alteration and mineralization features of four bore cores, induced activity polarization well logging and Eh-4 geophysical section, deep mineralization anomaly objects of Hongshan ore districts shows low resistance, middle and high polarization, measurements of Eh-4 consecutive conductance section show the existing of concealed porphyry ore body deeper than 450m, on the top of and around rock body there are low resistance body ranged from 100-300Ω•m, this area may be the ore-bearing part. In a word, Hongshan Cu-Au deposit deposit is a combine of upper HS-style epithermal Au deposit and deeper porphyry mineralization system. It has great potential to find large HS-style epithermal-porphyry Au-Cu deposits. This paper consists of seven chapters and twenty seven sections. The geological character of deposit is basic condition in this work. Constitute of oxidation zone, research of sulfate mineral, relation between oxidation and primary zone, K-Ar ages of potassic sulfate are key parts of thesis. Genesis of ore deposit is the further expansion of this research. Analysis of ore-controlling factors is the penetration above basic. Analysis of potential is application of exploration.

均衡重力异常计算与综合解释——以黄海地区为例

Petroleum and Natural Gas is an important strategic resources. The reserves of Petroleum and Natural Gas can’t meet the need of our country, which also blocks the development of economy and threatens the safety of national. Therefore, it makes a great sense to bring “the second round of oil & gas exploration” into effect and study the exploration of oil and gas of Pre-Cenozoic residual basins in China. The integrated geophysical exploration is the main way to research the Pre-Cenozoic residual basins. Gravity exploration is one of the most important exploration methods, which has played an important role in oil and gas prospecting, such as compartmentalizing geotectonic elements, delineating the distribution range of sedimentary basins, searching oil and gas structure, abstracting oil and gas information, and so on, from its naissance. The isostatic gravity anomalies is significant for exploration, which can help us research deep crustal structure, the equilibrium state of earth, the geologic structure of shallow crust, the basement shape of sedimentary basins and the genetic evolution of sedimentary basins. In the paper, we stress the implication and physical meanings systemically, and discuss the calculation theory. On the basis of previous work, we test different isostatic compensation models and parameters to find out their influences to the result of isostatic gravity anomalies. In addition, we improve the method of isostatic gravity anomalies calculation and give a system of isostatic gravity anomalies calculation which is proved has satisfying effect. From the research above, we find that the results of Platt model and Airy model are consistent, which have similar form and almost the same value. However, by contrast, the Airy model is proved has better adaptability than Platt model. The two main parameters——crust thickness and density difference of crust and mantle, both have influence to the isostatic gravity anomalies, but the latter have more. Finally, we adopt the regional field extending edge method to make the result more of actual geologic condition. On the methods above, we calculate the isostatic gravity anomalies field in Yellow Sea area from the Bouguer gravity anomalies and the water depth and altitude data. And then the isostatic gravity anomalies character is analyzed and the integrated geological-geophysical interpretation is made on the basis of summarizing the previous research result systemically and analyzing other geophysical data and geological information. From the research, we find that the Yellow Sea area belongs to continental type crust equilibrium regions, where the isostatic gravity anomalies field is placid and has less fluctuation values, which implies that the area is in equilibrium state to different extends.

黄海海域前新生代残留盆地综合地质地球物理研究

China’s annual oil import volume has been increasing in recent years, but the oil price in the international market fluctuates and poses a severe threat to China’s economic development and national security. Therefore, it is of great importance to study the gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea. Yellow Sea has widespread and thick Mesozoic and Paleozoic strata that contain multilayer source rock. Hence, Yellow Sea Mesozoic and Paleozoic strata have good conditions of forming Pre-Cenozoic hydrocarbon reservoirs. Pre-Cenozoic Residual Basins are usually buried deep and then transformed many times in its long evolutional history. These characteristics make it difficult to apply a single method in exploring Pre-Cenozoic Residual Basins. On the other hand, it is highly effective to solve key problems of gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea by using integrated geological and geophysical methods which make full use of the advantages of various exploring techniques. Based on the principle of “the region controls the local; the deep restricts the shallow,” this study focuses on Pre-Cenozoic Residual Basins in Yellow Sea to describe the structure frame of its distribution, with gravity, magnetic, seismic, drill-hole and geological data and previous research findings. In addition, the distribution characteristics of Pre-Cenozoic Residual Basins in Yellow Sea are also analyzed. This paper explores the characteristics of error between gravity forward with constant density and gravity forward with variable density through the study on 2-D and 3-D gravity forward in frequency domain. The result shows that there is a linear relationship between error and depth of 2-D geological model but there is a nonlinear relationship between error and depth of 3-D geological model. The error can be removed according to its linear characteristics or statistical nature of nonlinear characteristics. There is also error between gravity inversion with constant density and gravity inversion with variable density due to variable density and edge-effect. Since there are not noticeable rules between the error and the two causes as variable density and edge-effect, this study adopts gravity inversion with variable density and methods to eliminate the edge-effect in basement inversion to improve inversion accuracy. Based on the study on the rock physical properties and strata distribution of Yellow Sea and adjacent regions, this study finds that there is a big density contrast between Cretaceous-Jurassic strata and their substratum. The magnetic basement of south Yellow Sea is regarded as top of Archeozoic-Proterozoic early strata, and there are double magnetic basements in north Yellow Sea. Gravity and magnetic data are used to inverse the gravity basement and magnetic basement of Yellow Sea, with seismic and drill-hole data as constrains. According to data of gravity and magnetic basement distribution, the depth of Cenozoic strata and previous research findings, this paper calculates the thickness of the Mesozoic and Pre-Mesozoic Residual Basins, draws the distribution outline of Pre-Cenozoic Residual Basins in Yellow Sea, and analyzes its macro-distribution characteristics. Gravity inversion is applied on a typical geological profile in Yellow Sea to analyze the characteristics of its fractures and magnetic basements. The characteristics of Pre-Cenozoic Residual Basins distribution outline in Yellow Sea and the fractures and magnetic basements of its typical profile shown by profile inversion provides new geophysical evidence for these structure views such as “the South Yellow Sea and the North Yellow Sea belong to different structural units” and “Sino-Korea and Yangtze blocks combine along Yellow Sea East Fractured Zone in Yellow Sea”.

地球化学场的多重分形分析与分解

Geochemical and Geophysical anomaly play an important role in mineral exploration，their spatial structure character include singularity and self－similar。The singularity of an anomaly reflects the enrichment characters of the geochemical element ,The anomaly separation by multifractal model is useful in mineral anomaly assessment。In recent years, The continuous multifractal mode of the geochemical fields was studied ,it can be separated into the simple continuous multifractal mode and the high concentration multifractal mode, and the C－A、C－D、 S－A、MSDV、W－A method to decompose the anomaly were presented。Those are succeeded in interpretation of Geochemical and Geophysical anomaly。 This study makes a summarization to these method， we present a multifractal method based on wavelet transform to analyze the multifractal fields 。The singularity and spectrum are calculated through tracing the wavelet maximum modulus in different measure，and then the fields can be decomposed by the characters of the singularity。 It is demonstrated to be useful in interpretation of Au anomaly in Gekou－Shicheng region Rushang Shandong Province 。 Based on the multifractal theory , Using the concentration—area（C—A）method ,We study two geochemical fields in Chifeng area , Inner Mongolia。The results show that the geochemical fields have three different multifractal modes。Based on these ,we discuss the enrichment mode of the geochemical elements and their distributions in space and get the anomaly lower limit ,then the geochemical backgrounds、regional anomalies and local anomalies are distinguished

地幔压力下CaSiO3 Perovskite状态方程和弹性常数的第一性原理研究

In the past decade density functional theory (DFT) has made its way from a peripheral position in quantum chemistry to center. Of course the often excellent accuracy of the DFT based methods has provided the primary driving force of this development. This dissertation is devoted to the study of physical and chemical properties of planetary materials by first-principle calculation. The concerned properties include the geometry, elastic constants and anisotropy. In the first chapter, we give a systematic introduction to theoretical background and review its progress. Development of quantum chemistry promotes the establishment of DFT. Theorem of Hohenberg-Kohn is the fundament of DFT and is developed to Kohn-Sham equation, which can be used to perform real calculations. Now, new corrections and extensions, together with developed exchange-correlation, have made DFT more accurate and suitable for larger systems. In the second chapter, we focus on the calculational methods and technical aspects of DFT. Although it is important to develop methods and program, external package are still often used. At the end of this chapter, we briefly some widely used simulation package and the application of DFT. In the third chapter, we begin to focus on properties of real materials by first principles calculation. We study a kind of minerals named Ca perovskite, investigate its possible structure and anisotropy at Earth’s mental condition. By understanding and predicting geo-physically important materials properties at extreme conditions, we can get the most accurate information to interpret seismic data in the context of likely geophysical processes.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

奇异性分析方法及其在储层预测中的应用研究

In research field of oil geophysical prospecting, reservoir prediction is refers to forecasting physical properties of petroleum reservoir by using data of seismic and well logging, it is a research which can guide oil field development. Singularities of seismic and logging data are caused by the heterogeneity of reservoir physical property. It's one of important methods that using singularity characteristics of seismic and logging data to study the reservoir physical property in recently. Among them, realization of reservoir quantitative prediction by analyzing singularity of the data and enhancing transition description of data is difficulty in method research. Based on wavelet transform and the fractal theory, the paper studied the singularity judgment criterion for seismic and logging data, not only analyzed quantitative relation between singularity data and reservoir physical property, but also applied it in practical reservoir prediction. The main achievements are: 1. A new method which provides singular points and their strength information estimation at only one single scale is proposed by Herrmann (1999). Based on that, the dissertation proposed modified algorithm which realized singularity polarity detection. 2. The dissertation introduced onset function to generalize the traditional geologic boundaries variations model which used singularity characteristics to represent the abruptness of the lithologic velocity transition. We show that singularity analysis reveals generic singularity information conducted from velocity or acoustic impedance to seismogram based on the convolution seismic-model theory. Theory and applications indicated that singularity information calculated from seismic data was a natural attribute for delineating stratigraphy boundaries due to its excellent ability in detecting detailed geologic features. We demonstrated that singularity analysis was a powerful tool to delineate stratigraphy boundaries and inverse acoustic impedance and velocity. 3. The geologic significances of logging data singularity information were also presented. According to our analysis, the positions of singularities indicate the sequence stratigraphic boundary, and there is subtle relationship between the singularity strength and sedimentary environment, meanwhile the singularity polarity used to recognize stratigraphic base-level cycle. Based on all those above, a new method which provided sedimentary cycle analysis based on the singularity information of logging data in multiple scales was proposed in this dissertation. This method provided a quantitative tool for judging interface of stratum sequence and achieved good results in the actual application.

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

南襄盆地泌阳凹陷断裂斜坡带构造特征与油气成藏

The faulted slope zone of Biyang depression, a multiple hydrocarbon accumulation zone lying in a rich oil depression of Nanxiang basin, is a structural-sedimentary compounded slope, which is developed in Yanshanian period and has an area of 500 km2. From the ‘bottom up’, the developed strata may be divided into Yuhuangding formation in Neogene, Dacang Fang, Hetao-yuan plus Liaozhuang formations and Fenghuangzheng plus Pingyuan formation in Neogene, while Hetao-yuan formation is the main hydrocarbon-bearing target. Because of transtensional stress fields formed by persistent action of large-scale faulting in the south of the depression, sedimentary differential compaction in different stages, and tectonic inversion in later developing stage of the depression, a series of nose structure zones cut by different strike faults are developed. Therefore, the reservoir migration and accumulation are controlled by the complex faulted-nose structural zone, reservoir types are dominated by faulted-noses, faulted-blocks and fault-lithology, while lithology and stratigraphic unconformable reservoirs are locally developed. In combination with demands of practical production, applying with a new research approach of systematology and a combination with dynamic and static modes, guided by modern petroleum geologic theory, and based on previous data and studies, new techniques, methods of geophysical exploration, various computer simulation and forecasting techniques are applied in the new research of this paper. Starting from the structural features and formation mechanism, the forming mechanism of faulted structure, conditions and controlling factors of hydrocarbon accumulation, as well as various space-time allocation relationships in the process of accumulation are analyzed in the research. Besides that, the hydrocarbon migration, accumulation mechanism and dynamic evolution process are also discussed in the paper. Through the research, the accumulation rule of the faulted slope zone in faulted lake basin, the distribution and enrichment regularity of different reservoir controlling factors are systematically summarized. The summarizations indicate that the faulted slope is a favorable orientational zone, hydrocarbon is accumulated in nose structures and enriched in the main body of nose structures, faulted transformation zone and the ascent direction of laddering faulted blocks, the faults are the controlling factors, hydrocarbon accumulation zones controlled by fault-lithology are distributed along the faulting direcion. In the end, hydrocarbon migration and accumulation models of complex faulted-nose blocks are established. 1) Down cut model—‘flank-sheet’: the hydrocarbon is migrated like ‘sheet’ along a series of faults with parallel distribution and accumulated in the flank of nose structures; 2）Cross cut --‘axis-string’ model: the hydrocarbon cutting across the faults is migrated like ‘string’ and accumulated in the axis of nose structures. In view of different distribution models, reservoir forming combination patterns are divided and hydrocarbon reservoir evaluation exploration is carried out, which achieves good results in application. Key words: faulted slope zone; migration and accumulation model; reservoir controlling mechanism; reservoir-forming combination

陕北何家集区延长组长4＋5－长6储层特征研究

Abstract： Hejiaji area lies on eastern part of Shanbei Slope in Ordos Basin and the primary oil-bearing bed is Chang 4＋5 and Chang 6 of Yanchang Formation. It is indicated that the sedimentary facies and reservoir characteristics restricted the hydrocarbon accumulation regularity by the geological information. Therefore, Applied with outcrop observation，core description, geophysical logging interpretation, thin section determination, Scanning Electron Microscope, reservoir lithology and physical property analysis and other analytic machinery, the sedimentary facies ,micro-characteristic and master control factors on hydrocarbon reservoir of Yanchang Formation in Hejiaji area are studied deeply by means of sedimentology，reservoir geology and petroleum geology and provide a reliably reference for later prospect . Delta facies are identified in Hejiaji area and of which distributary channels in delta plain microfacies controlled the distribution of sand bodies and accumulation of oil and gas.The distribution of sand bodies distributed from northeast to southwest are dominated by sedimentary facies . It was shown that the sandstones are medium to granule arkose，which the mud matrix is r and including，calcite，the content of matrix is lower and that mostly are cements which are mainly quartz and feldspar overgrowths and chlorite films, in the second place are hydromica and ferrocalcite. All the sandstones have entered a period of late diagenetic stage in which the dominant diagenesis types in the area are compaction, cementation and dissolution. Remnant intergranular porosity and feldspar dissolved pore are main pore types which are megalospore and medium pore. Medium-fine throat, fine throat and micro-fine throat are the mainly throat type. Pore texture can be classified as megalospore and fine throat type, medium-pore and micro-fine throat type mainly, and they are main accumulate interspace in research region. The reservoir of Yanchang Formation in Hejiaji area is low- pore and low- permeability in the mass which have strong heterogeneity in bed, interbedded and plane. Studying the parameter of pore and permeability comprehensively and consulting prevenient study results of evaluation of reservoir, the reservoir is classifiedⅡ，Ⅲ and Ⅳ three types in which the Ⅱand Ⅲ can be divided into Ⅱa and Ⅱb, Ⅲa and Ⅲb respectively. Ⅱb and Ⅲa are the main reservoir type in Hejiaji area which are about 72.73％and 80％percent of whole reservoir and effective reservoir respectively.

陕北地区石家湾-李家岔探区三叠系延长组石油成藏条件分析

Shijiawan –Lijiacha area, lying on the northeastern part of the Shanbei Slope of Ordos Basin, was selected as studying area. The previous explorations proved that the 2nd segment and 6th segment of the Yanchang Formation are the most important oil-bearing formations. It is indicated that the sedimentary facies and reservoir characteristics restricted the hydrocarbon accumulation regularity. Therefore, with petrology methodologies, such as outcrop observation, core description, geophysical logging interpretation, thin section determination, scanning electron microscope, as well as rock property analysis, the reservoirs was were systematically studied and characterized. The sedimentary micro-facies, seals, reservoir-seal combines, migration pathways and entrapping modes were taken into account. The author tempted to establish a base for further studies on reservoirs and on petroleum geology, and to provide some reliably geological evidences for later prospect activities. It was found that the sediments in the 2nd and 3rd segments of the Yanchang Formation in Shijiawan –Lijiacha area were deposited in braided rivers, and most sandy-bodies were identified as channel sandbars. The 4+5th and 6th segments were principally deposited in deltaic-plain environment, consisting of corresponding sub-facies such as distributary channels, natural levee, crevasse-splay and marsh. The skeleton sandy-bodies were identified as sandy sediments of distributary channels. The sand grains in reservoir in studied area possess generally low mineralogical maturity and moderate structural maturity, and the form of pores may be classified into intergranular types and dissolved types. Most reservoirs of Yanchang Formation in Shijiawan –Lijiacha area belong to extreme low-porosity low-permeability ones (type III), and the 2nd sediments belongs to low permeability one (type II) and the 6th segment belong to super low-permeability one(type Ⅳ). The reservoirs in the 2nd segment behave more heterogeneous than those in the 6th segment. The statistic analysis results show that, for 6th and 4+5th segments, the high quality reservoir-seal combines may be found everywhere in the studied area except in the northwest and the southwest parts; and for 1st and 2nd segments, in the northeast, central and southwest parts Petroleum migration happened in the duration of the Early Cretaceous period in both lateral and vertical directions. The migration paths were mainly constructed by permeable sandy-bodies. The superimposed channel sandy-bodies consist of the principal part of the system of carriers. the vertical fractures, that may travel through the seals between reservoirs, offered the vertical paths for migrating oil. It may be synthesized that oil coming from south kitchens migrated first laterally in carriers in the 6th segment. When arrived at the studied area, oil will migration laterally or/and vertical within both the sandy-bodies and fractures, in a climbing-stair way. The results demonstrate that the oil was entrapped in traps structure-lithology and/or lithology traps. In some cases, the hydrodynamic force may help to trap oil. Accumulation of oil in the area was mainly controlled by sedimentary facies, seals, structure, and heterogeneity of reservoir in the 2nd, 4+5th and 6th segments. Especially, the oil distributions in both the 2nd and 6th segments were obviously influenced by seals in the 4+5th segment. The existence of seals in 1st segment seems important for accumulation in the 2nd segment.

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

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.

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

叠前地震多角度三参数同步反演方法

Pre-stack seismic inversion has become the emphasis and hotspot owing to the exploration & exploitation of oil field and the development of seismic technology. Pre-stack seismic inversion has the strongpoint of making the most of amplitude versus offset compared with the post-stack method. In this dissertation, the three parameters were discussed from multi-angle reflectance of P-wave data based on Zoeppritz’s and Aki & Richard’s equation, include P-wave velocity, S-wave velocity, and density. The three parameters are inversed synchronously from the pre-stack multi-angle P-wave data, based on rockphysics model and aimed at the least remnant difference between model simulation and practical data. In order to improve the stability of inversion and resolution to thin bed, several techniques were employed, such as the wavelet transform with multi-scale function, adding the Bayesian soft constraint and hard constraints (the horizon, structure and so on) to the inversion process. Being the result, the uncertainty of the resolution is reduced, the reliability and precision are improved, the significance of parameters becomes clearer. Meeting to the fundamental requirement of pre-stack inversion, some research in rockphysics are carried out which covered the simulation and inversion of S-wave velocity, the influence of pore fluids to geophysical parameters, and the slecting and analyzing of sensitive parameters. The difference between elastic wave equation modeling and Zoeppritz equation method is also compared. A series of key techniques of pre-stack seismic inversion and description were developed, such as attributes optimization, fluid factors, etc. All the techniques mentioned above are assembled to form a technique sets and process of synchronous pre-stack seismic inversion method of the three parameters based on rock physics and model simulation. The new method and technology were applied in many areas with various reservoirs, obtained both geological and economic significance, which proved to be valid and rational. This study will promote the pre-stack inversion technology and it’s application in hidden reservoirs exploration, face good prospects for development and application.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.