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

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”.

高温高压下H2O-CH4体系PVTx性质的分子动力学模拟研究

Geological fluids are important components in the earth system. To study thephysical chemistry properties and the evolution of fluid system turns out to be one of the most challenging issues in geosciences. Besides the conventional experimental approaches and theoretical or semi-theoretical modeling, molecular level computer simulation(MLCS) emerges as an alternative tool to quantificationally study the physico-chemical properties of fluid under extreme conditions in order to find out the characteristics and interaction of geological fluids in and around earth. Based on our previous study of the intermolecular potential for pure H2O and thestrict evaluation of the competitive potential models for pure CH4 and the ab initio fitting potential surface across H2O-CH4 molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH4 and the H2O-CH4 mixtures. Comparison of 1941 simulations with experimental PVT data for pure CH4 shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the H2O-CH4 mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase,indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region. After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the H2O-CH4 system covering 673 to 2573 K and 0.01 to 10 GPa. Isochores for compositions < 4 mol% CH4 up to 773 K and 600 MPa are also determined in this thesis.

承德公路隧道围岩分级与塌方机理分析

During tunnel constriction the classification of rock mass is widely used in tunnel design and construction. Moreover it offers the base information about tunnel investment and security. The quick classification of rock mass is very important for not delaying tunnel construction. Nowadays the tunnel engineers usually use initial survey files which are obtained by probe drilling to design a tunnel. It brings the problem that initial surrounding rock classification is usually much different from the real condition during the tunnel construction. Because initial surrounding rock lack credibility, it need us to make real time surrounding rock classification during the tunnel construction, and feed back the result to designers and constructors. Therefore, to find a quick wall rock classification method is very important not only for the time limit for a project but also for not delaying tunnel construction. Not all but many tunnels and underground constructions do suffer form collapse during the period of construction. Although accidental collapse in a large project in civil and geotechnical engineering sometimes appears to be a local event, if it occurred, it can bring about casualties, disrupted，production, construction delay, environmental damage, capital cost etc，therefore, it has been a difficult problem ,both in theory and in practice, establishing how to prevent underground structures form collapse and how to handle such an event in case in occurs. It is important to develop effective solutions and technical measures to prevent and control the collapse. According to the tunnel collapse occurred in Cheng De this paper analyze the main collapse mechanism leading to tunnel collapse and summon up the disposal method when collapse happened. It may be useful for tunnel construction in Cheng De in future. This paper is base on tunnel surrounding rock classification and tunnel support tasks during the tunnel construction in Cheng De area. It aims at solving 4 important problems in tunnel design and construction. 1) The relationship between rock rebound strength and rock single axle compression strength. First we go to the face wall and do rebound test on the tunnel face, then we chose some pieces of rock and do point loading test. Form the tests record we try to find the relationship between rock rebound strength and rock single axle compression strength. 2) The relationship between the value [BQ] and the value Q. First in order to obtain the information of rock character, rock strength, degree of weathering, the structure of rock mass, the joint condition, underground water condition and so on, we go to the tunnel face to do field investigation. And then we use two kinds of rock classification method to make surrounding rock classification. Base on the works above, finally we analyze the relationship between the value [BQ] and the value Q. 3) Sum up the mechanism leading to tunnel collapse and it disposal method in Cheng De area According to the tunnel collapse occurred in Cheng De this paper analyze the main reasons leading to the tunnel collapse and sum up the disposal method when collapse happened. 4) Obtain the properties of steel frame grid by numerical simulation. First we establish the 3D numeral model of steel frame grid by ADINA, and then find the mechanics properties by numerical simulation in ADINA. Second Based on the rock mass geological structure model, we established steel frame grid numeral model which is installed in the tunnel by FLAC3D and simulated the progress of tunnel construction. We hope that the support effect in tunnel can be evaluated from the numerical simulation.

基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

积分偏移的计算几何与并行实现

The Second Round of Oil & Gas Exploration needs more precision imaging method, velocity vs. depth model and geometry description on Complicated Geological Mass. Prestack time migration on inhomogeneous media was the technical basic of velocity analysis, prestack time migration on Rugged surface, angle gather and multi-domain noise suppression. In order to realize this technique, several critical technical problems need to be solved, such as parallel computation, velocity algorithm on ununiform grid and visualization. The key problem is organic combination theories of migration and computational geometry. Based on technical problems of 3-D prestack time migration existing in inhomogeneous media and requirements from nonuniform grid, parallel process and visualization, the thesis was studied systematically on three aspects: Infrastructure of velocity varies laterally Green function traveltime computation on ununiform grid, parallel computational of kirchhoff integral migration and 3D visualization, by combining integral migration theory and Computational Geometry. The results will provide powerful technical support to the implement of prestack time migration and convenient compute infrastructure of wave number domain simulation in inhomogeneous media. The main results were obtained as follows: 1. Symbol of one way wave Lie algebra integral, phase and green function traveltime expressions were analyzed, and simple 2-D expression of Lie algebra integral symbol phase and green function traveltime in time domain were given in inhomogeneous media by using pseudo-differential operators’ exponential map and Lie group algorithm preserving geometry structure. Infrastructure calculation of five parts, including derivative, commutating operator, Lie algebra root tree, exponential map root tree and traveltime coefficients , was brought forward when calculating asymmetry traveltime equation containing lateral differential in 3-D by this method. 2. By studying the infrastructure calculation of asymmetry traveltime in 3-D based on lateral velocity differential and combining computational geometry, a method to build velocity library and interpolate on velocity library using triangulate was obtained, which fit traveltime calculate requirements of parallel time migration and velocity estimate. 3. Combining velocity library triangulate and computational geometry, a structure which was convenient to calculate differential in horizontal, commutating operator and integral in vertical was built. Furthermore, recursive algorithm, for calculating architecture on lie algebra integral and exponential map root tree (Magnus in Math), was build and asymmetry traveltime based on lateral differential algorithm was also realized. 4. Based on graph theory and computational geometry, a minimum cycle method to decompose area into polygon blocks, which can be used as topological representation of migration result was proposed, which provided a practical method to block representation and research to migration interpretation results. 5. Based on MPI library, a process of bringing parallel migration algorithm at arbitrary sequence traces into practical was realized by using asymmetry traveltime based on lateral differential calculation and Kirchhoff integral method. 6. Visualization of geological data and seismic data were studied by the tools of OpenGL and Open Inventor, based on computational geometry theory, and a 3D visualize system on seismic imaging data was designed.

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

高分辨率非线性岩性物性联合反演方法和应用

In the prediction of complex reservoir with high heterogeneities in lithologic and petrophysical properties, because of inexact data (e.g., information-overlapping, information-incomplete, and noise-contaminated) and ambiguous physical relationship, inversion results suffer from non-uniqueness, instability and uncertainty. Thus, the reservoir prediction technologies based on the linear assumptions are unsuited for these complex areas. Based on the limitations of conventional technologies, the thesis conducts a series of researches on various kernel problems such as inversions from band-limited seismic data, inversion resolution, inversion stability, and ambiguous physical relationship. The thesis combines deterministic, statistical and nonlinear theories of geophysics, and integrates geological information, rock physics, well data and seismic data to predict lithologic and petrophysical parameters. The joint inversion technology is suited for the areas with complex depositional environment and complex rock-physical relationship. Combining nonlinear multistage Robinson seismic convolution model with unconventional Caianiello neural network, the thesis implements the unification of the deterministic and statistical inversion. Through Robinson seismic convolution model and nonlinear self-affine transform, the deterministic inversion is implemented by establishing a deterministic relationship between seismic impedance and seismic responses. So, this can ensure inversion reliability. Furthermore, through multistage seismic wavelet (MSW)/seismic inverse wavelet (MSIW) and Caianiello neural network, the statistical inversion is implemented by establishing a statistical relationship between seismic impedance and seismic responses. Thus, this can ensure the anti-noise ability. In this thesis, direct and indirect inversion modes are alternately used to estimate and revise the impedance value. Direct inversion result is used as the initial value of indirect inversion and finally high-resolution impedance profile is achieved by indirect inversion. This largely enhances inversion precision. In the thesis, a nonlinear rock physics convolution model is adopted to establish a relationship between impedance and porosity/clay-content. Through multistage decomposition and bidirectional edge wavelet detection, it can depict more complex rock physical relationship. Moreover, it uses the Caianiello neural network to implement the combination of deterministic inversion, statistical inversion and nonlinear theory. Last, by combined applications of direct inversion based on vertical edge detection wavelet and indirect inversion based on lateral edge detection wavelet, it implements the integrative application of geological information, well data and seismic impedance for estimation of high-resolution petrophysical parameters (porosity/clay-content). These inversion results can be used to reservoir prediction and characterization. Multi-well constrains and separate-frequency inversion modes are adopted in the thesis. The analyses of these sections of lithologic and petrophysical properties show that the low-frequency sections reflect the macro structure of the strata, while the middle/high-frequency sections reflect the detailed structure of the strata. Therefore, the high-resolution sections can be used to recognize the boundary of sand body and to predict the hydrocarbon zones.

东南沿海火山碎屑岩40Ar/39Ar年代学研究

The Mesozoic pyroclastic rocks cover a vast coastal area in southeastern China, which is an important part of volcanic belt around Pacific Ocean. However, the previous dates for these rocks are confused and conflicted with each other, which have limited the further researches. It is difficult to date pyroclastic rocks, for almost all the dating methods, due to the multiple sources during their formation. Single crystal laser probe 40Ar/39Ar method is a powerful means to date pyroclastic rocks with complex sources. By distinguishing the xenocrysts and altered materials, Singe crystal total fusions of CO2 lasing on the sanidine separates could yield rational 40Ar/39Ar results and distinguish their sources. Timing on formations of the Moshishan Group, after avoiding the exotic and altered grains by lasing on the single sanidine separate, was reported in this study. 40Ar/39Ar incremental heating results of Cement part of the pyroclastic rocks show that the age spectrums are too complex to interpret for geological significance because of the alteration and 39Ar recoil. Incremental heating on K-feldspar separate from pyroclastic rocks give reliable data. Combining 40Ar/39Ar incremental heating ages and laser total fusion ages, we suggest the age was 125Ma-120Ma for Dashuang Fm, 120-119 Ma for Gaowu Fm, 119-114Ma for Xishantou Fm, 114-112Ma for Chawan Fm and 112-110Ma for Jiuliping Fm. The most intense period of volcanic activity in eastern Zhejiang Province was at about 120Ma. These new ages are much younger than the previous ones, suggesting that these thick volcanic formations had been formed in very short durations.

地震层析成像与天山地区的壳幔深部结构研究

The Tien Shan is the most prominent intracontinental mountain belt on the earth. The active crustal deformation and earthquake activities provide an excellent place to study the continental geodynamics of intracontinental mountain belt. The studies of deep structures in crust and upper mantle are significantly meaningful for understanding the geological evolution and geodynamics of global intracontinental mountain belts. This dissertation focuses on the deep structures and geodynamics in the crust and upper mantle in the Tien Shan mountain belt. With the arrival time data from permanent and temporal seismic stations located in the western and central Tien Shan, using seismic travel time tomographic method, we inversed the P-wave velocity and Vp/Vs structures in the crust and uppermost mantle, the Pn and Sn velocities and Pn anisotropic structures in the uppermost mantle, and the P-wave velocity structures in the crust and mantle deep to 690km depth beneath the Tien Shan. The tomographic results suggest that the deep structures and geodynamics have significant impacts not only on the deformations and earthquake activities in the crust, but also on the mountain building, collision, and dynamics of the whole Tien Shan mountain belt. With the strongly collision and deformations in the crust, the 3-D P-wave velocity and Vp/Vs ratio structures are highly complex. The Pn and Sn velocities in the uppermost mantle beneath the Tien Shan, specially beneath the central Tien Shan, are significantly lower than the seismic wavespeed beneath geological stable regions. We infer that the hot upper mantle from the small-scale convection could elevate the temperature in the lower crust and uppermost mantle, and partially melt the materials in the lower crust. The observations of low P-wave and S-wave velocities, high Vp/Vs ratios near the Moho and the absences of earthquake activities in the lower crust are consistent with this inference. Based on teleseismic tomography images of the upper mantle beneath the Tien Shan, we infer that the lithosphere beneath the Tarim basin has subducted under the Tien Shan to depths as great as 500 km. The lithosphere beneath the Kazakh shield may have subducted to similar depths in the opposite direction, but the limited resolution of this data set makes this inference less certain. These images support the plate boundary model of converge for the Tien Shan, as the lithospheres to the north and south of the range both appear to behave as plates.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

陕北何家集区延长组长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.

松辽盆地南部腰英台区块低渗透储层整体压裂改造技术研究

The Yaoyingtai Block is located within the northeastern Changling Depression of southern Songliao Basin, where the reservoir sandstones are petrophysically characterized by very low permeability, which results in the low success probability of artificial fracturing, and the low oil yield by water injection in the course of oil production. In order to improve the situations as stated above, this research aims to work out an integral fracturing technology and strategy applicable to the low permeable reservoirs in Yaoyingtai Block. Under the guidance of geological theory, reservoir engineering and technology, the subsurface occurrences of natural and hydraulic fractures in the reservoirs are expected to be delineated, and appropriate fracturing fluids and proppants are to be optimized, based on the data of drilling, well logging, laboratory and field experiments, and geological data. These approaches lay the basis of the integral fracturing technology suitable for the low permeable reservoir in the study area. Based on core sample test, in-situ stress analysis of well logging, and forward and inversion stress field modeling, as well as fluid dynamic analysis, the maximum in-situ stress field is unraveled to be extended nearly along the E-W direction (clustering along N85-135°E) as is demonstrated by the E-W trending tensional fractures. Hydraulic fractures are distributed approximately along the E-W direction as well. Faulting activities could have exerted obvious influences on the distribution of fractures, which were preferentially developed along fault zones. Based on reservoir sensitivity analysis, integrated with studies on rock mechanics, in-situ stress, natural fracture distribution and production in injection-production pilot area, the influences of primary fractures on fracturing operation are analyzed, and a diagnostic technology for primary fractures during depressurization is accordingly developed. An appropriate fracturing fluid (hydroxypropyl guar gum) and a proppant (Yixing ceramsite, with a moderate-density, 0.45-0.9mm in size) applicable to Qingshankou Formation reservoir are worked out through extensive optimization analysis. The fracturing fluid can decrease the damage to the oil reservoir, and the friction in fracturing operation, improving the effect of fracturing operation. Some problems, such as sand-out at early stage and low success rate of fracturing operations, have been effectively solved, through pre-fracturing formation evaluation, “suspension plug” fracturing, real-time monitoring and limited-flow fracturing. Through analysis of fracture-bearing tight reservoir with variable densities and dynamic analysis of influences of well patterns on fracturing by using numerical simulation, a fracturing operation scheme for the Qingshankou Formation reservoir is proposed here as being better to compress the short factures, rather than to compress the long fractures during hydraulic fracturing. It is suggested to adopt the 450m×150m inverted 9-spot well pattern in a diamond shape with wells placed parallel to fractures and a half fracture length of 60-75m.

准噶尔盆地西北缘三叠系层序地层及有利勘探目标选择

Based on outcrop, borehole, seismic and regional geological data, the sequence stratigraphy, sedimentary facies of the Triassic in the western margin of the Zhugaer basin was studied, and favorable exploration target was forecasted. The major achievements include: (1) the Triassic in the western margin of the Zhugaer basin can be divided into 1 second-order sequence and 5 third-order sequences, which are, in ascending order, TSQ1, TSQ2, TSQ3, TSQ4, and TSQ5. TSQ1 is equivolent to Baikouquan formation, TSQ2 is equivolent to lower Kelamayi formation, TSQ3 is equivolent to upper Kelamai formation, TSQ4 is equivolent to lower and middle Baijiantan formation, and TSQ5 is equivolent to upper Baijiantan formation. Each sequence is divided into transgressive and regressive system tracts. Thus the sequence correlation framework is established. (2) The factors controlling development of sequences are analyzed, and it is believed that tectonic is the major controlling factor. Model of sequence development is summarized. (3)Through study on sedimentary facies, 6 types of facies are recognized: alluvial fan, fan delta, braided river, braided delta, delta and lake. Their microfacies are also recognized. In this study, it is proposed that the upper and lower Kelamayi formation（TSQ2、 TSQ3）is deposited by braided river instead of alluvial fan. This conclusion is of important theoretical and practical significance.(4) The sedimentary facies map of each sequence is compiled, and the sedimentary facies developed in each sequence is determined. In TSQ1, the sedimentary facies developed is alluvial fan and fan delta. In TSQ2, the sedimentary facies developed is mainly alluvial fan and fan delta in the north, and braided river and braided delta in the south. In TSQ3, the sedimentary facies developed is mainly braided river and braided delta. In TSQ4, the sedimentary facies developed is mainly braided delta in the north, and meandering delta in the south. In TSQ5, the sedimentary facies developed is mainly braided river and braided delta. (5) In the framework of sequence stratigrahpy, favorable areas for concealed traps are forecasted, and different types of traps are developed in different system tracts. (6) Favorable areas for future exploration are predicted.

油气勘探风险评价与决策技术研究

PetroChina and other national petroleum incorporations need rigorous procedures and practical methods in risk evaluation and exploration decision at home and abroad to safeguard their international exploration practice in exploration licence bidding, finding appropriate ratio of risk sharing with partners, as well as avoiding high risk projects and other key exploration activities. However, due to historical reasons, we are only at the beginning of a full study and methodology development in exploration risk evaluation and decision. No rigorous procedure and practical methods are available in our exercises of international exploration. Completely adopting foreign procedure, methods and tools by our national incorporations are not practical because of the differences of the current economic and management systems in China. The objective of this study is to establish a risk evaluation and decision system with independent intellectual property right in oil and gas exploration so that a smooth transition from our current practice into international norm can take place. The system developed in this dissertation includes the following four components: 1. A set of quantitative criteria for risk evaluation is derived on the basis of an anatomy of the parameters from thirty calibration regions national wide as well as the characteristics and the geological factors controlling oil and gas occurrence in the major petroleum-bearing basins in China, which provides the technical support for the risk quantification in oil and gas exploration. 2. Through analysis of existing methodology, procedure and methods of exploration risk evaluation considering spatial information are proposed. The method, utilizing Mahalanobis Distance (MD) and fuzzy logic for data and information integration, provides probabilistic models on the basis of MD and fuzzy logic classification criteria, thus quantifying the exploration risk using Bayesian theory. A projection of the geological risk into spatial domain provides a probability map of oil and gas occurrence in the area under study. The application of this method to the Nanpu Sag shows that this method not only correctly predicted the oil and gas occurrence in the areas where Beibu and Laoyemiao oil fields are found in the northwest of the onshore area, but also predicted Laopu south, Nanpu south and Hatuo potential areas in the offshore part where exploration maturity was very low. The prediction of the potential areas are subsequently confirmed by 17 exploration wells in the offshore area with 81% success, indicating this method is very effective for exploration risk visualization and reduction. 3. On the basis of “Methods and parameters of economic evaluation for petroleum exploration and development projects in China”, a ”pyramid” method for sensitivity analysis was developed, which meets not only the need for exploration target evaluation and exploration decision at home, but also allows a transition from our current practice to international norm in exploration decision. This provides the foundation for the development of a software product “Exploration economic evaluation and decision system of PetroChina” (EDSys). 4. To solve problem in methodology of exploration decision, effort was made on the method of project portfolio management. A drilling decision method was developed employing the concept of geologically risked net present value. This method overcame the dilemma of handling simultaneously both geological risk and portfolio uncertainty, thus casting light into the application of modern portfolio theory to the evaluation of high risk petroleum exploration projects.