火成岩储层裂缝识别、预测理论与方法研究－以济阳坳陷火成岩油藏为例

Since 1970s, igneous reservoirs such as Shang741, Bin674 and Luol51 have been found in Jiyang depression, which are enrichment and heavy-producing. Showing good prospect of exploration and development, igneous reservoirs have been the main part of increasing reserves and production in Shengli oilfield. As fracture igneous reservoir being an extraordinary complex concealed reservoir and showing heavy heterogeneity in spatial distribution, the study of recognition, prediction, formation mechanism and the law of distribution of fracture is essential to develop the reservoir. Guided by multiple discipline theory such as sedimentology, geophysics, mineralogy, petroleum geology, structural geology and reservoir engineering, a set of theories and methods of recognition and prediction of fractured igneous rock reservoir are formed in this paper. Rock data, three-dimensional seismic data, log data, borehole log data, testing data and production data are combined in these methods by the means of computer. Based on the research of igneous rock petrography and reservoir formation mechanism, emphasized on the assessment and forecast of igneous rock reservoir, aimed at establishing a nonhomogeneity quantification model of fractured igneous rock reservoir, the creativity on the fracture recognition, prediction and formation mechanism are achieved. The research result is applied to Jiyang depression, suggestion of exploration and development for fractured igneous rock reservoir is supplied and some great achievement and favourable economic effect are achieved. The main achievements are gained as follows: 1. The main facies models of igneous rock reservoir in JiYang depression are summarized. Based on data and techniques of seism, well log and logging,started from the research of single well rock facies, proceeded by seismic and log facies research, from point to line and line to face, the regional igneous facies models are established. And hypabyssal intrusion allgovite facies model, explosion volcaniclastic rock facies model and overfall basaltic rocks facies model are the main facies models of igneous rock reservoir in JiYang depression. 2. Four nonhomogenous reservoir models of igneous reservoirs are established, which is the base of fracture prediction and recognition. According to characteristics of igneous petrology and spatial types of reservoir, igneous reservoirs of Jiyang depression are divided into four categories: fractured irruptive rock reservoir, fracture-pore thermocontact metamorphic rock and irruptive rock compound reservoir, pore volcanic debris cone reservoir and fracture-pore overfall basaltic rock reservoir. The spatial distribution of each model's reservoir has its features. And reservoirs can be divided into primary ones and secondary ones, whose mechanism of formation and laws of distribution are studied in this paper. 3. Eight geologic factors which dominate igneous reservoirs are presented. The eight geologic factors which dominates igneous reservoirs are igneous facies, epigenetic tectonics deformation, fracture motion, intensity of intrusive effect and adjoining-rock characters, thermo-contact metamorphic rock facies, specific volcano-tectonic position, magmatic cyclicity and epigenetic diagenetic evolution. The interaction of the eight factors forms the four types nonhomogenous reservoir models of igneous reservoirs in Jiyang depression. And igneous facies and fracture motion are the most important and primary factors. 4. Identification patterns of seismic, well log and logging facies of igneous rocks are established. Igneous rocks of Jiyang depression show typical reflecting features on seismic profile. Tabular reflection seismic facies, arc reflection seismic facies and hummocky or mushroom reflection seismic facies are the three main facies. Logging response features of basic basalt and diabase are shown as typical "three low and two high", which means low natural gamma value, low interval transit-time, low neutron porosity, high resistivity and high density. Volcaniclastic rocks show "two high and three low"-high neutron porosity, high interval transit-time, low density, low-resistance and low natural gamma value. Thermo-contact metamorphic rocks surrounding to diabase show "four high and two low" on log data, which is high natural gamma value, high self-potential anomaly, high neutron porosity, high interval transit-time and low density and low-resistance. Based on seismic, well log and logging data, spatial shape of Shang 741 igneous rock is described. 5. The methods of fracture prediction and recognition for fractured igneous reservoir are summarized. Adopting FMI image log and nuclear magnetic resonance log to quantitative analysis of fractured igneous reservoir and according to formation mechanism and shape of fracture, various fractures are recognized, such as high-angle fracture, low-angle fracture, vertical fracture, reticulated fracture, induced fracture, infilling fracture and corrosion vug. Shang 741 intrusive rock reservoir can be divided into pore-vug compound type, pore fracture type, micro-pore and micro-fracture type. Physical properties parameters of the reservoir are computed and single-well fracture model and reservoir parameters model are established. 6. Various comprehensive methods of fracture prediction and recognition for fractured igneous reservoir are put forward. Adopting three-element (igneous facies, fracture motion and rock bending) geologic comprehensive reservoir evaluation technique and deep-shallow unconventional laterolog constrained inversion technique, lateral prediction of fractured reservoir such as Shang 741 is taken and nonhomogeneity quantification models of reservoirs are established.

孤南-三合村洼陷含油气系统研究

The development petroleum geology has made people from studying and studying and predicting in statically and respectively the pool-forming conditions of an area such as oil source bed, reservoir, overlying formation, migration, trap and preservation, etc. to regarding these conditions as well as roles of generation, reservation and accumulation as an integrated dynamic evolution development system to do study .Meanwhile apply various simulating means to try to predict from quantitative angle. Undoubtedly, the solution of these questions will accumulate exploration process, cut down exploration cost and obtain remarkable economic and social benefits. This paper which take sedimentology ,structural geology and petroleum geology as guides and take petroleum system theory as nucleus and carry out study thinking of beginning with static factor and integration of point and face as well as regarding dynamic state factor as factor and apply study methods of integration of geology, Lab research and numerical modeling proceed integrated dissect and systematic analysis to GuNan-SanHeCun depression. Also apply methods of integration of sequence stratigraphy, biostratigraphy, petrostratigraphy and seismic data to found the time-contour stratigraphic framework and reveal time-space distribution of depositional system and meantime clarify oil-source bed, reservoir and overlying distribution regular patterns. Also use basin analysis means to study precisely the depositional history, packed sequences and evolution. Meanwhile analyze systematically and totally the fracture sequence and fault quality and fault feature, study the structural form, activity JiCi and time-space juxtaposion as well as roles of fault in migration and accumulation of oil and gas of different rank and different quality fault. Simultaneously, utilize seismic, log, analysis testing data and reservoir geology theory to do systematic study and prediction to GuNan-SanHeCun reservoir, study the reservoir types macroscopic distribution and major controlling factors, reservoir rock, filler and porosity structural features as well as distribution of reservoir physical property in 3D space and do comprehensive study and prediction to major controlling and influential factors of reservoir. Furthermore, develop deepingly organic geochemistry comprehensive study, emphasis on two overlaps of oil source rock (ESI, ES3) organic geochemistry features, including types, maturity and spatial variations of organic matter to predict their source potential .Also apply biological marks to proceed oil-to-source correlation ,thereby establish bases for distribution of petroleum system. This study recover the oil generation history of oil source rocks, evaluate source and hydrocarbon discharge potential ,infer pool-forming stages and point out the accumulation direction as well as discover the forming relations of mature oil-source rock and oil reservoir and develop research to study dynamic features of petroleum system. Meanwhile use systematic view, integrate every feature and role of pool forming and the evolution history and pool-forming history, thereby lead people from static conditions such as oil source bed, reservoir, overlying formation, migration, trap and preservation to dynamically analyzing pool-forming process. Also divide GuNan-SanHeCun depression into two second petroleum system, firstly propose to divide second petroleum system according to fluid tress, structural axis and larger faults of cutting depression, and divide lower part of petroleum system into five secondary systems. Meanwhile establish layer analysis and quantitative prediction model of petroleum model, and do quantitative prediction to secondary petroleum system.

准噶尔盆地侏罗系层序格架和油气成藏条件研究

Jurassic is an important hydrocarbon-bearing formation in Junggar Basin. Analyzing in strata sequence stratigraphy and hydrocarbon formation has both theoretical and practical values. First of all, strata sequence stratigraphy of continental facies is a new development and supplement in the theory of stratigraphic geology. Stratum of continental facies, unlike sea facies, has rich sup-plements, rapid facies changes, and was influenced slightly by sea level changes. The structural background and sedimentary environment of the basin in west China are greatly different from those of the basins in east China. So it is important to build the patterns of strata sequence stratigraphy in west China basins. Secondly, it is also of significance to find out all kinds of traps, for the dominant types are structural ones so far. After 50 years exploration, the stratigraphic or litholigic traps have become the main concern. This desertation is mainly focused on establishing the isochronal strtaum frame for Junggar Basin to show the evolvement characters of the basin sediment system and the regionalstuctrue background. By analyzing the conditions and patterns of the regional oil and gas bearing formations with typical cross-sections, we have established the patterns of sedimentary conditions for different reserviors. By authur's study, several fruitful results have been obtained in the following: Strata sequence frame and evolvement characteristics of Jurassic: By studing strata sequence, Jurassic has been divided into 2 second rank strata sequences and 3 third rank strata sequences based on the interface unconformities. Only 2 fourth rank strata sequences were grouped in BaDaoWan group. Also different seismic facies and sediment units have been recognized with the establishment of the of sediment system model. The oil-gas system characteristics in Jurassic: We conclude that hydrocar bon resources have the best oil potential. Potential of coal, carbonaceous and dark mudstone were reduced in turn. In this thesis we have made the evaluation of three hydrocarbon sources and the distribution oil-gas resource, and studied the potentials of hydrocarbon and evolvement for each kind of micro-component of the two main resource rocks. Prediction of paleo-temperature: In Junggar basin the evolvement of paleo-ground temperature can be divided into three stages. From Carboniferous to early Permian grads of ancient ground temperature was 8-5 ℃/100m, 5-3 ℃/100m from later period of Permian to end Trias, 3-2 ℃/100m from Jurassic to early Tertiary. Patterns of Jurassic hydrocarbon-bearing reserviors: There were two kinds of hydrocarbon source of Permian and Jurassic. They form different hydrocarbon - bearing systems. Six fundamental hydrocarbon - bearing trap modeS have been established. Directions for later exploration: There were two kinds of regional belts in Jurassic, One is structural belt caused by Yanshan and Ximalaya process, and the other was the stratum one caused by paleostructural rises.

弱各向异性介质中传播的地震波及其速度重建

We begin our studies to make the best of information of seismic data and carry out the description of cracks parameters by extracting anisotropic information. The researching contents are: (1) velocity and polarization anomaly of seismic wave (qP and qSV wave) in weak anisotropic media; (2) reflection seismic synthetic record in anisotropic media; (3) multiple scattering induced by cracks; (4) anisotropic structure inversion and velocity reconstruction with VSP (Vertical Seismic Profile) data; (5) multi-parameters analysis of anisotropy in time-domain and depth-domain. Then we obtain results as follows: (1) We achieve approximate relation of qP and qSV wave's velocity and polarization property in weak anisotropic media. At the same time, we calculate anisotropic velocity factors and polarization anomaly of several typical sedimentary rocks. The results show there are different anisotropic velocity factors and polarization anomaly in different rocks. It is one of the primary theoretical foundation which is expected to identify lithology; (2) We calculate reflection seismic synthetic record with theoretical model; (3) We simulate scattering induced by cracks with Boundary Element Method. Numerical studies show that in the presence of cracks; spatial and scale-length distributions are important and cannot be ignored in modeling cracked solids; (4) From traveltimes information of VSP data, we study the velocity parameter inversion of seismic wave under isotropic and anisotropic models, and its result indicate that the inversion imaging under anisotropic model will not destroy the original features of isotropic model, but it will bring on some bigger error if we adopt the method of isotropic model for anisotropic model data. Further more, basing on the study we develop the CDP mapping technology of reflecting structure under isotropic and anisotropic models, and we process real data as a trial of the methods; (5) We study the problem of initial model reconstruction of anisotropic parameters structure represented by Anderson parameter in depth domain for surface data.

山东东营凹陷浊积岩油藏形成机制及预测技术

More and more exploration practises and research results prove that deep-water turbidites and slump turbidites in abrupt slope belt as well as slump turbidites in deposit center are favorable places to find subtle reservoirs. Based on sedimentology, sequence stratigraphy, tectonics and petroleum geology and with the tool of computer, this paper focused a set of questions such as the forming mechanism and m&a of turbidite reservoir in Dongying Sag, and obtained a systematic research result from available information of geology, seismic survey, well logging and testing, etc. Since the structural pattern and tectonic stress are controlling factors on distribution of turbidite reservoirs, deep-water sequence and shallow-water sequence have close relation with various turbidites. The author primarily disclose the regular pattern of the turbidites in the abrupt slope trouth, and predict that the abrupt slope is mainly connected to near-bank slump turbidite fans and deep-water turbidite fans; whereas the depocenter area usually developed with front delta slump fans. According to this, seismic-geology models corresponding to various deposit systems can be set up as a principle for predicting the turbiditea in faulted lacustrine basins. Based on comprehensive study of the turbidite reservoir forming condition and main controlling factors, a differential accumulation model and a multiple migration&accumulation model were also set up for abrupt slope turbidities. Referenced the m

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

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

东营凹陷北部陡坡带砂砾岩体性油气藏形成机理及分布规律研究

The north steep slope zone of Dongying Depression has great potential in oil resource and as the usage of 3-d seismic data in the last decade, the exploration of oil and gas has get into the stage of sandy glavel body lithological oil-gas pool exploration. In this thesis, writer take the north steep slope zone of Dongying Depression as target area and take Sha-III and Sha-[V Menber as purpose stratum, study on sequence stratigraphy, depositional system, reservoir description, emphasesing on analyzing of forming of lithological oil-gas pool, especially the dynamics principle of oil and gas preliminary movement from the source rock to the reservoir form lithological oil-gas pools. The aim of this work is to give some quantitatively explanation for the mechanism of lithological oil-gas pool forming, and set up the theory of pool form with characteristic terrestrial faulted basin. There are main conclusions and views as follow. 1. Applying with principle of sequence stratigrapgy, according to the depositional cycles of Dongying Depression, the sequence stratigraphical partition of Tertiary was finished, stressing on dismembering Sha-III and Sha-IV Menber as 5system tracts. 2. The structure of Dongying Depression especially of the north steep slope zone has accomplished, including the analyzing the structural cortroling to depositional condition of the north steep slope zone of Dongying Depression, discussed relationship between the structure of the north steep slope zone and the pool-forming. 3. The horizontal and vertical exchanges of ancient climates and ancient physiognomy of the all stratum units and studies on characteristic of depositional system distribution have been finished, found that there are five depositional systems in the north steep slope zone of Dongying Depression as fluvial, delta (tan-delta), sub-water fluvial fan lacustrine, gravitive flow, and seven formations of sandy glavel body, and forecasting of all kinds of sandy glavel body has been made. 4. Seismic stratigraphy and log stratigraphy have been made, described and forecasted all kinds of reservoir of objective stratum by means of physical geography method, setup a series means of sandy glavel body description suit to target area. 5. The pool-forming system has been studied, analyzing all the elements in petroleum sub-system of Sha-III and Sha-IV Menber of Dongying Depression with view of source controlling, estimated the petroleum system applying source rock potential index combining with distribution ofreservior. 6.Through studying types of pool, the controlling factors of pool-forming of sandy glavel body were discussed by deposition stages, formation types, structure ect. as a conclusion that the characteristics of pool forming in the north steep slope zone of Dongying Depression are, the controlling factor of the pools is mainly lithology, petrophysics of oil sands vary greatly, with a large heterogeneity, all kind of reservoir with different formation has different pool-forming conditions, and as a result, formed various pools of sandy glavel body along the steep slope with regular combination, distribution and constituted the multiple petroleum accumulative pattern. 7. It's the first time to cauculate and estimate the fluid pressure in source rock of Dongying Depression, set up the stratum fluid pressure in Dongying Depression, and firstly use equivalent charging pressure and reservoir forming index to quantitatively evaluate the pool-forming condition of lithological pool.8. Above all studies, follow up the scent of the exploration combined with practice a lot of explorative targets were found, and got geat economic and social benefit.

叠前深度偏移速度模型的建立方法研究

Seismic exploration is the main tools of exploration for petroleum. as the society needs more petroleum and the level of exploration is going up, the exploration in the area of complex geology construction is the main task in oil industry, so the seismic prestack depth migration appeared, it has good ability for complex construction imaging. Its result depends on the velocity model strongly. So for seismic prestack depth migration has become the main research area. In this thesis the difference in seismic prestack depth migration between our country and the abroad has been analyzed in system. the tomographical method with no layer velocity model, the residual curve velocity analysical method based on velocity model and the deleting method in pre-processing have been developed. In the thesis, the tomographysical method in velocity analysis is been analyzed at first. It characterized with perfection in theory and diffculity in application. This method use the picked first arrivial, compare the difference between the picked first arrival and the calculated arrival in theory velocity model, and then anti-projected the difference along the ray path to get the new velocity model. This method only has the hypothesis of high frequency, no other hypothesis. So it is very effective and has high efficiency. But this method has default still. The picking of first arrival is difficult in the prestack data. The reasons are the ratio of signal to noise is very low and many other event cross each other in prestack data. These phenomenon appear strongly in the complex geology construction area. Based on these a new tomophysical methos in velocity analysis with no layer velocity model is been developed. The aim is to solve the picking problem. It do not need picking the event time contiunely. You can picking in random depending on the reliability. This methos not only need the pick time as the routine tomographysical mehtod, but also the slope of event. In this methos we use the high slope analysis method to improve the precision of picking. In addition we also make research on the residual curve velocity analysis and find that its application is not good and the efficiency is low. The reasons is that the hypothesis is rigid and it is a local optimizing method, it can solve seismic velocity problem in the area with laterical strong velocity variation. A new method is developed to improve the precision of velocity model building . So far the pattern of seismic prestack depth migration is the same as it aborad. Before the work of velocity building the original seismic data must been corrected on a datum plane, and then to make the prestack depth migration work. As we know the successful example is in Mexico bay. It characterized with the simple surface layer construction, the pre-precessing is very simple and its precision is very high. But in our country the main seismic work is in land, the surface layer is very complex, in some area the error of pre-precessing is big, it affect the velocity building. So based on this a new method is developed to delete the per-precessing error and improve the precision of velocity model building. Our main work is, (1) developing a effective tomographical velocity building method with no layer velocity model. (2) a new high resolution slope analysis method is developed. (3) developing a global optimized residual curve velocity buliding method based on velocity model. (4) a effective method of deleting the pre-precessing error is developing. All the method as listed above has been ceritified by the theorical calculation and the actual seismic data.

山西堡子湾金矿角砾岩带控矿规律、成矿时代及隐伏矿体定位预测

The Puziwan gold deposit is located at the northern boundary of the North China platform where there are very favorable ore-forming conditions. The deposit is a medium-sized gold deposit associated with silver and copper polymetallic minerallization which were closely related to cryptoexplosive processes and fluidization and underwent multi-epoch superposition of mineralization. The mineralization consists of cryptoexplosive breccia-type, fluidizing-type, porphyry-type, quartz vein-type, etc. The ore-controlling law of ore-bearing brecciated zone is systematically studied in the paper. The shape, scale and attitude of ore bodies are dominat in the upper, the ore bodies of fluidizing breccia-type are dominant in the upper and the middle part, the porphyry-type ore bodies are dominant in the bottom. The quartz vein-type ore bodies are impenetrated in all the brecciated zone. The metallogenic epoch of Puziwan gold deposit is collated and stipulated. The Rb-Sr isochron age of quartz porphyry (wall rock) is 233 Ma, refering to the Ar-Ar age of gold-bearing quartz, excluded the former conclusion that the so-called metallogenic epoch (245.9Ma) is in Indosinian epoch. The nonage metallogenic age of Buziwan gold deposit shoule be in the Yanshanian epoch (142.5Ma). By applying the sub-specimen sampling technique, the Rb-Sr isochron age of gold-bearing pyrites in late mineralization epoch is dated to be 64 Ma. In conclusion, the metallogenetic epoch of Puziwan gold deposit is in late Yanshanian-early Himalayan epoch. On the above basis, the metallogenic model of the "train type" and new texture model of ore deposit are established. by applying the methods of geology, shallow seismic exploration, remote sensing, partial geochemical extraction and the study of inclusions in the late granite porphyry, the author has made the metallogenic progosis in the depth and the prephery of Puziwan gold deposit and eight prospecting targets are proposed.

Lg尾波衰减成像研究——北京及邻区衰减成像

With the development of the technology of earthquake observation, more and more researchers work at many fields' of seismicity using seismic kinetic property, as the result, the study of attenuation has also made great progress, especially in the mechanism of the attenuation and the physical process. Aki put forward single back scattering theory to explain the forming of the seismic coda wave in 1969. Then, researchers started to develop the study in seismic scattering and attenuation. My thesis is also based on that theory. We assume that the Lg wave is a superposition fo high-mode surface waves, the coda of Lg is caused by scattering. Sato proposed Single Isotropic Scattering model (SIS model) to interpreted the scatter property, and he also formulated the geometrical spreading term. Then Xie (1988) developed the single spectral-ratio (SSR) method to obtain the Lg coda Q and the frequency dependent factor n. Later, he get to lateral images in the area of scatter ellipse. SSR method is explored and used in the study of Lg coda waves of regional earthquakes in my thesis. Choosing the earthquakes records with high ratio of signal-noise ,which were recorded at the stations from 1989 to 1999, we obtain the single trace Lg coda Q and its frequency dependent factor n. The results proved that SIS model is the reasonable model to explain the Lg coda wave, and SSR method also can be used to process Lg coda of regional earthquakes to get to the satisfied Lg coda Q. Based on the Lg coda Q we obtained using the former method, we explore the programs to inverse the regional Lg coda Q independently, and then make use of them to inverse the Lg coda Q of Beijing and adjacent area. The inversion result is satisfied. We conclude that the distribution of Qo (Q in lHz) is marked by the inhomogeneity, which is related to the tectonic structure: The value of Qo in uplift area, for example, Yanshan uplift, Taihang uplift, Luxi Uplift, is higher than the depression area, for example, Jizhong depression, Huanghua depression, and Jiyang depression, and the border between the higher Q area and lower Q area is very clear; Lg coda Q is also related to the velocity structure, higher velocity area is also with higher Q, lower velocity area is with lower Q; and higher heat-flow area is companied with lower Q. All in all, the value of Q reflects the difference of characteristics of lithofaces, porosity, the liquid content between the pores and heat flow. So, the Q value difference between uplift area and depression area reveals the difference of tectonic structure, lithology and physical character of the rock. So, the study of Lg coda Q is help to understand the earthquakes propagation mechanism through the inhomogenous medium, the cause of the coda, attenuation mechanism of the coda. Making use of the lateral images of Q, with velocity images, heat flow results, and other experimental result, we will be promoted to understand the complex structure of the crust, its inhomogenous character, and so on.

地震勘探三维可视化方法研究与实现

This thesis mainly studies the technologies of 3-D seismic visualization and Graphic User Interface of seismic processing software. By studying Computer Graphics and 3-D geological modeling, the author designs and implements the visualization module of seismic data processing software using OpenGL and Motif. Setting seismic visualization flow as the subject, NURBS surface approximation and Delaunay Triangulation as the two different methods, the thesis discusses the key algorithms and technologies of seismic visualization and attempts to apply Octree Space Partitioning and Mip Mapping to enhance system performance. According to the research mentioned above, in view of portability and scalability, the author adopts Object-oriented Analysis and Object-oriented Design, uses standard C++ as programming language, OpenGL as 3-D graphics library and Motif as GUI developing tool to implement the seismic visualization framework on SGI Irix platform. This thesis also studies the solution of fluid equations in porous media. 2-D alternating direction implicit procedure has been turned into 3-D successive over relaxation iteration, which possesses such virtues as faster computing speed, faster convergence rate, better adaptability to heterogeneous media and less memory demanding.

DAS24-3B数据采集系统控制系统软件的研制与开发

Under the auspices of the 'knowledge-Innovation Program' of CAS, Institute of Geology and Geophysics has established the Broadband Seismic Laboratory. A new kind of 24-bit high-resolution seismograph DAS24-3B has been designed and manufactured in an effort of developing China's own technology of seismic array. Since these instruments will primarily be used in field operation, there is a need to optimize the system software of data acquisition system (DAS) to enhance its stability, compatibility and maintenance. The design ideas of the system software of DAS24-3B are partly learned from the advanced DAS 72A-08. In this system there are two exclusive communication programs DNAPI-COM1 and DNAPI-LPT1, which are suitable for all standard industrial computers with ECP parallel port and serial port. By these exclusive parallel and serial communication interface the system software is split into three parts, acquisition program, user's control program and graphical display program, which can function well in separate units and can run correctly in whole. The three parts of DAS24-3B's system software possess different functions and advantages. The function of acquisition program is to control the process of seismic data acquisition. DAS24-3B system reduced its power and harddisk read-write disturbance by using the extended memory attached to its CPU, which functions as enlarging the data buffer of system and lessening the times of harddisk read-write operations. Since GPS receiver of DAS is strongly sensitive to the around environment and has the possibility of signal loss the acquisition program has been designed with the ability to automatically trail the GPS locked time. The function of user's controlling program is to configure the system's work environment, to inform the user's commands to DAS, to trail the status of DAS in real-time. The function of graphical display program is to illustrate data in figures, to convert data file into some common formatted file, to split data file in parts and combine data files into one. Both user's control program and graphical display program are API (Application Programming Interface) in window 95/98 system. Both possess the features of clearness and friendship by use of all kind of window controls, which are composed by menu, toolbar, statusbar, dialogue box, message box, edit box, scrollbar, time control, button and so on. Two programs of systemic exception handles are provided to treat the trouble in field. The DAS24-3B DAS has been designed to be easier to use-better ability, more stable and simpler. It has been tested in field and base station and has been proved more suitable for field operation of seismic array than other native instruments.

青藏高原壳幔形变与层圈耦合数值模拟

Based on multi-principle (such as structures, tectonics and kinematics) exploratory data and related results of continental dynamics in the Tibetan plateau, the author reconstructed the geological-geophysical model of lithospherical structure and tectonic deformation, and the kinetics boundary conditions for the model. Then, the author used the numerical scheme of Fast Lagrangian Analysis of Continua (FLAC), to stimulate the possible process of the stress field and deformational field in the Tibetan plateau and its adjacent area, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. With the above-mentioned results, the author discussed the relationship between crustal movement in shallow layer and the deformational process in interior layers, and its possible dynamic constraints in deep. At the end of the paper, an integrative model has been put forward to explain the outline images of crust-mantle deformation and coupling in the Tibetan Plateau. (1) The characteristics of crust-mantle structure of the Tibetan plateau have been shown to be very complex, and vertical and horizontal difference is significant. The general characteristics of crust-mantle of the Tibetan plateau may be that it's layering in depth direction, and shows blocking from south to north and belting from east to west, mainly according to the results of about 20 seismic sections, such as wide-angle seismic profiles, CMP, seismic tomography and so on. (2) The crust had shortened about 2200km, while the shortening is different for different block from south to north in the Tibetan plateau. It is about 11.5mm/a in Himalayan block, about 9.0mm/a in Lhas-Gangdese block, about 7.0mm/a in Qiangtang block and Songpan-Ganzi-Kekexili block, about 8.0mm/a in Kunlun-Qaidam, and about ll.Omm/a in Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. Which - in demonstrates the shortening rate decreases from south to north, but this rate increases near the north edge of the Tibetan plateau. The crust thickening rate is about 0.4mm/a in the whole Tibetan plateau; and this rate is about 0.5mm/a in Himalayan block, about 0.4mm/a in Lhas-Gangdese block, about 0.3mm/a in Qiangtang block, about 0.2mm/a in Songpan-Ganzi-Kekexili block and about O.lmm/a in Kunlun-Qaidam-Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. This implies that the thickening rate decreases in the blocks of the Tibetan plateau. From south to north, the displacement of eastern boundary in the Tibetan plateau is about 37mm/a in Himalayan block, about 45mm/a in Lhas-Gangdese block, about 47mm/a in Qiangtang block, about 43mm/a in Songpan-Ganzi-Kekexili block, and about 35mm/a in Kunlun-Qaidam-Qilian block, since the collision-matching between the Indian continent and Eurasia continent had happened about 50Ma ago. This implies that the rate of eastward displacement is biggest in the middle of plateau, and decreases to both sides. The transition of S-N compression stress field in Tibetan Plateau, since about 28Ma+ ago, may be caused by two reasons: On one hand, the movement direction of Eurasia continent changed from northward to southward about 28Ma± ago in the northern plateau. On the other hand, the front belt that is located between India continent's and Eurasia continent's convergence-collision, had moved southward to high Himalayan from Indus-Brahmaputra suture almost at the same time in southern plateau. Affected by the stress field, the earlier tectonics rotated clockwise, NE and NW conjugate strike-slip faults developed, and the SN rift formed. This indicated that the EW movement started. The ratio between upper crust and lower crust of different blocks from south to north in the Tibetan plateau during the process of deformation are as following: about 3.5~5:1 in Himalayan block, about 1~5: 3-4 (which is about 1:3o--4 in south and about 4~5:3 in north) in Lhas-Gangdese block, about 1:3~447mm/a in these blocks: Which is located to the north of Banggong-nujiang suture.

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

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

沾化凹陷泥质岩油气藏成藏特征研究

Mudstone reservoir is a subtle reservoir with extremely inhomogeneous, whose formation is greatly related to the existence of fracture. For this kind of reservoir, mudstone is oil source rock, cover rock and reservoir strata, reservoir type is various, attitude of oil layer changes greatly, and the distribution of oil and gas is different from igneous or clastic rock reservoir as well as from carbonate reservoir of self-producing and self-containing of oil and gas. No mature experience has been obtained in the description, exploration and development of the reservoir by far. Taking Zhanhua depression as an example, we studied in this thesis the tectonic evolution, deposit characteristics, diagenesis, hydrocarbon formation, abnormal formation pressure, forming of fissure in mudstone reservoir, etc. on the basis of core analysis, physical simulation, numerical simulation, integrated study of well logging and geophysical data, and systematically analyzed the developing and distributing of mudstone fissure reservoir and set up a geological model for the formation of mudstone fissure reservoir, and predicted possible fractural zone in studied area. Mudstone reservoir mainly distributed on the thrown side of sedimentary fault along the sloping area of the petroleum generatiion depression in Zhanhua depression. Growing fault controlled subsidence and sedimentation. Both the rate of subsidence and thickness of mudstone are great on the thrown side of growing fault, which result in the formation of surpressure in the area. The unlocking of fault which leads to the pressure discharges and the upward conduct of below stratum, also makes for the surpressure in mudstone. In Zhanhua depression, mudstone reservior mainly developed in sub-compacted stratum in the third segment of Shahejie formation, which is the best oil source rock because of its wide spread in distribution, great in thickness, and rich in organic matter, and rock types of which are oil source mudstone and shale of deep water or semi-deep water sediment in lacustrine facies. It revealed from core analysis that the stratum is rich in limestone, and consists of lamina of dark mudstone and that of light grey limestone alternately, such rock assemblage is in favor of high pressure and fracture in the process of hydrocarbon generation. Fracture of mudstone in the third segment of Shahejie formation was divided into structure fracture, hydrocarbon generation fracture and compound fracture and six secondary types of fracture for the fist time according to the cause of their formation in the thesis. Structural fracture is formed by tectonic movement such as fold or fault, which develops mainly near the faults, especially in the protrude area and the edge of faults, such fracture has obvious directivity, and tend to have more width and extension in length and obvious direction, and was developed periodically, discontinuously in time and successively as the result of multi-tectonic movement in studied area. Hydrocarbon generation fracture was formed in the process of hydrocarbon generation, the fracture is numerous in number and extensively in distribution, but the scale of it is always small and belongs to microfracture. The compound fracture is the result of both tectonic movement and hydrocarbon forming process. The combination of above fractures in time and space forms the three dimension reservoir space network of mudstone, which satellites with abnormal pressure zone in plane distribution and relates to sedimentary faces, rock combination, organic content, structural evolution, and high pressure, etc.. In Zhanhua depression, the mudstone of third segment in shahejie formation corresponds with a set of seismic reflection with better continuous. When mudstone containing oil and gas of abnormal high pressure, the seismic waveform would change as a result of absorb of oil and gas to the high-frequency composition of seismic reflection, and decrease of seismic reflection frequency resulted from the breakage of mudstone structure. The author solved the problem of mudstone reservoir predicting to some degree through the use of coherent data analysis in Zhanhua depression. Numerical modeling of basin has been used to simulate the ancient liquid pressure field in Zhanhua depression, to quantitative analysis the main controlling factor (such as uncompaction, tectonic movement, hydrocarbon generation) to surpressure in mudstone. Combined with factual geologic information and references, we analyzed the characteristic of basin evolution and factors influence the pressure field, and employed numerical modeling of liquid pressure evolution in 1-D and 2-D section, modeled and analyzed the forming and evolution of pressure in plane for main position in different periods, and made a conclusion that the main factors for surpressure in studied area are tectonic movement, uncompaction and hydrocarbon generation process. In Zhanhua depression, the valid fracture zone in mudstone was mainly formed in the last stage of Dongying movement, the mudstone in the third segment of Shahejie formation turn into fastigium for oil generation and migration in Guantao stage, and oil and gas were preserved since the end of the stage. Tectonic movement was weak after oil and gas to be preserved, and such made for the preserve of oil and gas. The forming of fractured mudstone reservoir can be divided into four different stages, i.e. deposition of muddy oil source rock, draining off water by compacting to producing hydrocarbon, forming of valid fracture and collecting of oil, forming of fracture reservoir. Combined with other regional geologic information, we predicted four prior mudstone fracture reservoirs, which measured 18km2 in area and 1200 X 104t in geological reserves.