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

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.

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

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

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

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

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

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.

断层封闭性及其与油气成藏的关系研究-以沾化凹陷孤岛地区为例

Godao area, located in the east of the Zhanhua depression of Jirang sag in Bohai Bay Basin, is the studying area in my dissertation. It is first time that fault sealing properties and the related relationship with the pool forming are studied in Gudao area. On the base of the analysis of the regional tectonics, the author has studied the tertiary structural evolution of the Gudao area and distinguished the fault's level and put forward the distinguishing principle. The geometrical feature, mechanical characters, developmental mechanism and history of the boundary faults in the tectonic unit of this area are all studied and emphasized especially. The buried history of oil-generating depression (that is Gudao depression) and the history of oil and gas migration simultaneously are discussed, the juxtaposition relationship between boundary fault evolutionary history and oil and gas migrated history are expatiated. To the geological condition of the Gudao area, three level faults sealing properties of this area were discussed in detail. Their characteristics of behavior and the intrinsic relationship between their sealing and oil and gas migrated reservoir are elucidated. The pool-forming models related to fault seal are exposed. The author has studied the lithologies of different order of faults, the relationship of occurrence assemblage analysis, normal stress of fault plane in different depth and shale smear factor faults. Then analysis their role in the various faults sealing and confirms fault sealing marks of three different orders faults and exposes the mechanism of fault sealing. Shale smear zone formed by first order fault in lasting activities is one different type of fault breccia and mainly controlled factor to its entrapment of petroleum. Effective sealing threshold value and fault displacement is ascertained. Mainly controlled factor of second fault sealing is bigger compressive stress loaded on fault plane. According to this, quantitatively evaluated index is given. Shale smear zone is necessary condition for second fault stress entrapment. The juxtaposed relationship between the different lithologies within third order faults is most important controlled factor for its sealing. Based on various order of fault sealing features and mechanism in Gudao area, the author proposed three orders of fault sealing models. Shale smear zone sealing model, normal stress sealing model and lithologies juxtaposed sealing model are suggested to first, second and third order fault respectively. The conclusion of this studying has not only the very important theoretical significance and practical value in Gudo area but also the very important guiding role for other areas of related aspects.

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.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.

砂砾岩储层流动单元四维模型研究

This article is an important part of "95" technological subject of SINOPEC. It has a large number of difficulties and workloads, and has significant theoretical meanings and practical value. The study area is composed of sandstone & conglomerate reservoir of alluvial fan & fan delta, which belong to Sha3 lower member and Sha4 upper member of lower tertiary of Yong'an Town Oilfield in Dongying Depression. The target stataum develops in the hanging wall of the synsedimentary fault in the scarp zone of Dongying Depression. The frequently intense movements result in the variation of sandstone and conglomerate reservoir and the evolution of the time and space of Sha3 lower member and Sha4 upper member in Yong'an Town Oilfield. As a result, it is difficult for the individual reservoir correlation at the root of fan, which bring about a tackle problem for the exploitation of oilfield. In this background, the research of fluid units will be more difficult. In this article, the new concepts, the new methods, and the new techniques of sedimentology, petroleum geology, reservoir geology, physics of crystal surface, dynamic & static state reservoir description and well logging geology are synthetically applied, and the computer technology are made full uses of, and the identifying, dividing and appraising of the two-formation-type sandstone & conglomerate reservoir fluid units of Sha3 lower member and Sha4 upper member systemically analyzed in Yong'an Town Oilfield, Dongying Depression. For the first time, the single-well model, the section model, the plane model, the nuclear magnetism log model, the microcosmic network model, the 4-D geology model and the simulation model of the two-formation-type reservoir fluid units of the of sandstone & conglomerate reservoir of Sha3 lower member and Sha4 upper member are established, and the formative mechanism and distributing & enrichment laws of oil-gas of the two type of sandstone and conglomerate reservoir fluid units are revealed. This article established the optimizing, identifying, classifying and appraising standard of the two-formation-type reservoir fluid units of the of sandstone and conglomerate reservoir of Sha3 lower member and Sha4 upper member, which settles the substantial foundations for static state model of the fluid units, reveals the macroscopic & microcosmic various laws of geometrical static state of the fluid units, and instructs the oil exploitation. This article established static state model of the two-formation-type sandstone and conglomerate reservoir fluid units by using the multi-subject theories, information and techniques, and reveals the geometrical configuration, special distribution and the oil-gas enrichment laws of the sandstone and conglomerate reservoir fluid units. For the first time, we established the nuclear magnetism log model of the two-formation-type sandstone and conglomerate reservoir of Sha3 lower member and Sha4 upper member, which reveals not only the character and distributing laws of the porosity and permeability, bat also the formation and distribution of the movable fluid. It established six type of microcosmic net model of the two-formation-type sandstone and conglomerate reservoir of Sha3 lower member and Sha4 upper member in the working area by using the advanced theories, such as rock thin section, SEM, image analysis, intrusive mercury, mold, rock C.T. measure & test image etc., which reveals the microcosmic characteristic of porosity & throat, filterate mode and microcosmic oil-gas enrichment laws of the sandstone and conglomerate reservoir. For the first time, it sets up the 4-D model and mathematic model of the sandstone and conglomerate reservoir, which reveals the distributing and evolving laws of macroscopic & microcosmic parameters of the two-formation-type sandstone and conglomerate reservoir and oil-gas in 4-D space. At the same time, it also forecasts the oil-gas distribution and instructs the oilfield exploitation. It established reservoir simulation model, which reveals the filterate character and distributing laws of oil-gas in different porosity & throat net models. This article established the assistant theories and techniques for researching, describing, indicating and forecasting the sandstone and conglomerate reservoir fluid units, and develops the theories and techniques of the land faces faulted basin exploitation geology. In instructing oilfield exploitation, it had won the notable economic & social benefits.

濮城油田沙三中6-10高分辨率层序地层学研究

The foundation of reservoir model and residual oil prediction have been the core of reservoir detailed description for improved oil production and enhanced oil recovery. The traditional way of sandstone correlation based on the geometrical similarity of well-logs which emphasizes "based on the cycle and correlating from larger to smaller" has shown its theoretical limits when explaining the correlating and the scale, geometry, continuity, connectivity of sandstones and the law of the reservoir property. It has been an urgent and difficult subject to find new theory and methods to solve the reservoir correlation and property prediction. It's a new way to correlate strata and found framework of reservoir through the process-response analysis in the base-level cycles. And it is also possible to analyze the reservoir property in reservoir framework. Taking the reservoir of zonation 6-10 in S3~2 of Pucheng Oil Field in Henan Province as an example, we founded the detailed reservoir stratigraphic framework through base-level correlation. In the strata frame, sediment distribution and its development are discussed based on sediment volume partitioning and facies differentiation analysis. Reservoir heterogeneities and its relation to base-level are also discussed. The analysis of primary oil distribution shows the base-level controlled oil distribution in reservoir. In this paper, subjects as following are discussed in detail. Based on the analysis of sedimentary structure and sedimentary energy, the facies model was founded. Founding stratigraphy framework through base level analysis In the studying zone, one long term cycle, 6 middle term cycles and 27 short term cycles was identified and correlated. 3 Predicting the property of reservoir for improving oil development The base level controlled the property of sandbody. The short and very short term cycle controlled the pattern of heterogeneities in sandbody, and the middle and long term cycle controlled the area and inter-layer heterogeneities. On the lower location of the middle and long term base level, the sandbody is well developed, with a wide area and large thickness, while on the high location of base level, there is an opposite reservoir character. 4 The studying of reservoir development response and oil distribution making a solid base for development adjustment Primary oil distribution is controlled by base level location. It tells that the sandbody on the high base level location was poor developed for its difficulty to develop. While on the low location of the base level, the sandbody is well developed for its relative easy to develop and dominant role in the development, but high residual oil for its high original oil content.

滚石机理及滚石灾害评价研究

Rockfall is a geological evolution process involving detachment of blocks or boulders from a slope face, then their free falls, bouncing, rolling or sliding, and finally deposition near the toe of the slope. Many facts indicate that the rockfall can cause hazards to peoples, and it can be regarded as a geological hazard. A rockfall event may only involve a boulder or rock, and also several ones. When there are peoples, buildings, or other man-made establishments within the scope of rockfall trajectory, losses will be possibly induced in tenns of human lives or damages to these facilities. Researches into mechanism, kinematics, dynamics, hazard assessment, risk analysis, and mitigation measures of rockfalls are extremely necessary and important. Occurrence of rockfall is controlled by a lot of conditions, mainly including topographical, geomorphic, geological ones and triggering factors. The rockfall especially in mountainous areas, has different origins, and occurs to be frequent, unexpected, uncertain, in groups, periodic and sectional. The characterization and classification of the rockfalls not only increase knowledge about rockfall mechanism, but also can instruct mitigation of the hazards. In addition, stability of potential rockfalls have various sensitivity to different triggering factors and changes of geometrical conditions. Through theoretical analyses, laboratory experiments and field tests, the author presents some back-analysis methods for friction coefficients of sliding and rolling, and restitution coefficients. The used input data can be obtained economically and accurately in the field. Through deep studies on hazard assessment methods and analysis of factors influencing rockfall hazard, this paper presents a new assessment methodology consisting of preliminary assessment and detailed one. From the application in a 430 km long stretch of the Highway, which is located between Paksho and Nyingtri in Tibet, the methodology can be applicable for the rockfall hazard assessment in complex and difficult terrains. In addition, risk analyses along the stretch are conducted by computing the probability of encountering rockfalls and life losses resulting from rockfall impacts. Rockfall hazards may be mitigated by avoiding hazardous areas, clearness of dangerous rocks, reinforcement, obstructing the rockfalls, leading the rockfalls, warning and monitoring for rockfalls, etc. Seen from present remedial level of rockfall hazards, different mitigation measures, economical and effective buffering units, monitoring tecliniques and consciousness of environmental protection for rockfall mitigations should be further developed.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

地震反射界面形态与介质各向异性属性联立成像：地震反射波旅行时与偏振信息综合反演

It has been a difficult problem faced by seismologists for long time that how exactly to reconstruct the earth's geometric structure and distribution of physical attributes according to seismic wave's kinematical and dynamic characteristics, obtained in seismological observation. The jointing imaging of seismic reflector and anisotropy attributes in the earth interior is becoming the research hot spot. The limitation of shoot and observation system makes that the obtained seismic data are too scarce to exactly reconstruct the geological objects. It is popular that utilizing only seismic reflection traveltimes or polarizations information make inversion of the earth's velocity distribution by fixing seismic reflector configuration (vice versa), these will lead to the serious non-uniqueness reconstruction due to short of effective data, the non-uniqueness problem of reconstructing anisotropy attributes will be more serious than in isotropy media. Obviously it is not enough to restrict the media structure only by information of seismic reflection traveltimes or polarizations, which even sometimes will lead to distorted images and misinterpretation of subsurface structure. So we try to rebuild seismic reflection structure (geometry) and media anisotropic structure (physics) in the earth interior by jointing data of seismic wave kinematics and dynamics characteristics, we carry out the new experiment step by step, and the research mainly comprises of two parts: one is the reconstruction of P-wave vertical velocity and anisotropic structure(Thomsen parameter s and 8) in the transversely isotropic media with vertical symmetrical axis(VTI) by fixing geometrical structure, and the other is the simultaneous inversion of the reflector surface conformation and seismic anisotropic structure by jointing seismic reflection traveltimes and polarizations data. Simulated annealing method is used to the first research part, linear inversion based on BG theory and Simulated annealing are applied to the second one. All the research methods are checked by model experiments, then applied to the real data of the wide-angle seismic profile from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. The results are as following The inversion results based on jointing seismic PP-wave or PSV-wavereflection traveltimes and polarizations data are more close to real model than themodels based simply on one of the two data respectively. It is shown that the methodwe present here can effectively reconstruct the anisotropy attributes in the earth'sinterior when seismic reflector structure is fixed.The layer thickness, P-wave vertical velocity and Thomsen anisotropicparameters {s and 8) could be resolved simultaneously by jointing inversion ofseismic reflection traveltimes and polarizations with the linear inversion methodbased on BG theory.The image of the reflector structure, P-wave vertical velocity and theanisotropy parameters in the crust could be obtained from the wide-angle seismicprofile from Tunxi (in Anhui Province), to Wenzhou (in Zhejiang Province). Theresults reveal the difference of the reflector geometrical structure and physicalattributes in the crust between Yangtze block and Cathaysia block, and attempt tounderstand the characteristics of the crustal stress field in the areas.

疏水性有机污染物的非线性吸附过程及滞后现象的研究

该论文测定了疏水性有机污染物在地质吸附剂上的吸附-解吸等温线,选用合适的模型拟合了这些等温线,并且研究了非线性吸附过程、解吸滞后现象与壤/沉积物异质性之间的关系;得到以下几点认识:1、利用红外光谱仪和核磁共振谱仪研究了土壤/沉积物有机质、土壤胡敏酸的组成.2、测定了菲、蔡、1,3,5-三氯苯在土壤/沉积物及处理样品上的吸附等温线.所有地质吸附剂上的吸附等温线表现出明显的非线性.Freundlic上模型能很好的拟合了有机物的吸附等温线,双区位反应模型(DRDM)能很好的拟合有机物在部分土壤/沉积物及处理样品上的吸附等温线,并且DRDM模型清楚的反映出有机物在低浓度和高浓度下不同的吸附特征.3、对相同土壤/沉积物样品,萘的吸附等温线要比菲的吸附等温线的线性要好些,这与它们的物化性质有关.4、对碱提取后残留土壤测定了菲、茶、1,3,5-三氯苯的吸附-解吸等温线.5、菲和萘在土壤/沉积物及处理样品上的吸附-解吸过程存在明显的滞后现象.

Tensor density matrix theory for determining population and alignment of symmetric top molecule using rotationally unresolved fluorescence

General expressions used for transforming raw laser-induced fluorescence (LIF) intensity into the population and alignment parameters of a symmetric top molecule are derived by employing the density matrix approach. The molecular population and alignment are described by molecular state multipoles. The results are presented for a general excitation-detection geometry and then applied to some special geometries. In general cases, the LIF intensity is a complex function of the initial molecular state multipoles, the dynamic factors and the excitation-detection geometrical factors. It contains a population and 14 alignment multipoles. How to extract all initial state multipoles from the rotationally unresolved emission LIF intensity is discussed in detail.