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

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.

改善三维波动方程叠前深度偏移处理效果的方法研究

3D wave equation prestack depth migration is the effective tool for obtaining the exact imaging result of complex geology structures. It's a part of the 3D seismic data processing. 3D seismic data processing belongs to high dimension signal processing, and there are some difficult problems to do with. They are: How to process high dimension operators? How to improve the focusing? and how to construct the deconvolution operator? The realization of 3D wave equation prestack depth migration, not only realized the leap from poststack to prestack, but also provided the important means to solve the difficult problems in high dimension signal processing. In this thesis, I do a series research especially for the solve of the difficult problems around the 3D wave equation prestack depth migration and using it as a mean. So this thesis service for the realization of 3D wave equation prestack depth migration for one side and improve the migration effect for another side. This thesis expatiates in five departs. Summarizes the main contents as the follows: In the first part, I have completed the projection from 3D data point area to low dimension are using de big matrix transfer and trace rearrangement, and realized the liner processing of high dimension signal. Firstly, I present the mathematics expression of 3D seismic data and the mean according to physics, present the basic ideal of big matrix transfer and describe the realization of five transfer models for example. Secondly, I present the basic ideal and rules for the rearrange and parallel calculate of 3D traces, and give a example. In the conventional DMO focusing method, I recall the history of DM0 process firstly, give the fundamental of DMO process and derive the equation of DMO process and it's impulse response. I also prove the equivalence between DMO and prestack time migration, from the kinematic character of DMO. And derive the relationship between DMO base on wave equation and prestack time migration. Finally, I give the example of DMO process flow and synthetic data of theoretical models. In the wave equation prestak depth migration, I firstly recall the history of migration from time to depth, from poststack to prestack and from 2D to 3D. And conclude the main migration methods, point out their merit and shortcoming. Finally, I obtain the common image point sets using the decomposed migration program code.In the residual moveout, I firstly describe the Viterbi algorithm based on Markov process and compound decision theory and how to solve the shortest path problem using Viterbi algorithm. And based on this ideal, I realized the residual moveout of post 3D wave equation prestack depth migration. Finally, I give the example of residual moveout of real 3D seismic data. In the migration Green function, I firstly give the concept of migration Green function and the 2D Green function migration equation for the approximate of far field. Secondly, I prove the equivalence of wave equation depth extrapolation algorithms. And then I derive the equation of Green function migration. Finally, I present the response and migration result of Green function for point resource, analyze the effect of migration aperture to prestack migration result. This research is benefit for people to realize clearly the effect of migration aperture to migration result, and study on the Green function deconvolution to improve the focusing effect of migration.

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

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.

波动方程辛几何算法三维叠前深度偏移及其并行实现

Oil and scientific groups have been focusing on the 3D wave equation prestack depth migration since it can solve the complex problems of the geologic structure accurately and maintain the wave information, which is propitious to lithology imaging. The symplectic method was brought up by Feng Kang firstly in 1984 and became the hotspot of numerical computation study. It will be widely applied in many scientific field of necessity because of its great virtue in scientific sense. This paper combines the Symplectic method and the 3-D wave equation prestack depth migration to bring up an effectual numerical computation method of wave field extrapolatation technique under the scientific background mentioned above. At the base of deep analysis of computation method and the performance of PC cluster, a seismic prestack depth migration flow considering the virtue of both seismic migration method and Pc cluster has formatted. The software, named 3D Wave Equation Prestack Depth Migration of Symplectic Method, which is based on the flow, has been enrolled in the National Bureau of Copyright (No. 0013767). Dagang and Daqing Oil Field have now put it into use in the field data processing. In this paper, the one way wave equation operator is decompounded into a phase shift operator and a time shift operator and the correct item with high rank Symplectic method when approaching E exponent. After reviewing eliminating alias frequency of operator, computing the maximum angle of migration and the imaging condition, we present the test result of impulse response of the Symplectic method. Taking the imaging results of the SEG/EAGE salt and overthrust models for example and seeing about the imaging ability with complex geologic structure of our software system, the paper has discussed the effect of the selection of imaging parameters and the effectuation on the migration result of the seismic wavelet and compared the 2-D and 3-D prestack depth migration result of the salt mode. We also present the test result of impulse response with the overthrust model. The imaging result of the two international models indicates that the Symplectic method of 3-D prestack depth migration accommodates great transversal velocity variation and complex geologic structure. The huge computing cost is the key obstruction that 3-D prestack depth migration wave equation cannot be adopted by oil industry. After deep analysis of prestack depth migration flow and the character of PC cluster ,the paper put forward :i)parallel algorithms in shot and frequency domain of the common shot gather 3-D wave equation prestack migration; ii)the optimized setting scheme of breakpoint in field data processing; iii)dynamic and static load balance among the nodes of the PC cluster in the 3-D prestack depth migration. It has been proven that computation periods of the 3-D prestack depth migration imaging are greatly shortened given that adopting the computing method mentioned in the paper. In addition,considering the 3-D wave equation prestack depth migration flow in complex medium and examples of the field data processing, the paper put the emphasis on: i)seismic data relative preprocessing, ii) 2.5D prestack depth migration velocity analysis, iii)3D prestack depth migration. The result of field data processing shows satisfied application ability of the flow put forward in the paper.

中国大陆及邻区地震活动背景、环境触发与地震预测研究

This dissertation that includes most of the P. PH.D research work during 2001~2002 covers the large-scale distribution of continental earthquakes in mainland China, the mechanism and statistic features of grouped strong earthquakes related to the tidal triggering, some results in earthquake prediction with correlativity analysis methods, and the flushes from the two strong continental earthquakes in South Asia in 2001. Mainland China is the only continental sub-plate that is compressed by collision boundaries at its two sides, within which earthquakes are dispersive and distributed as seismic belts with different widths. The control capability of the continental block boundaries on the strong earthquakes and seismic hazards is calculated and analyzed in this dissertation. By mapping the distribution of the 31282 ML:3s2,0 earthquakes, I found that the depth of continental earthquakes depend on the tectonic zonings. The events on the boundaries of relatively integrated blocks are deep and those on the new-developed ruptures are shallow. The average depth of earthquakes in the West of China is about 5km deeper than that in the east. The western and southwestern brim of Tarim Basin generated the deepest earthquakes in mainland China. The statistic results from correlation between the grouped M7 earthquakes and the tidal stress show that the strong events were modulated by tidal stress in active periods. Taking Taiwan area as an example, the dependence of moderate events on the moon phase angles (￡>) is analyzed, which shows that the number of the earthquakes in Taiwan when D is 50° ,50° +90° and 50° +180° is more than 2 times of standard deviation over the average frequency at each degree, corresponding to the 4th, 12th and 19th solar day after the new moon. The probability of earthquake attack to the densely populated Taiwan island on the 4th solar day is about 4 times of that on other solar days. On the practice of earthquake prediction, I calculated and analyzed the temporal correlation of the earthquakes in Xinjinag area, Qinghai-Tibet area, west Yunnan area, North China area and those in their adjacent areas, and predicted at the end of 2000 that it would be a special time interval from 2001 to 2003, within which moderate to strong earthquakes would be more active in the west of China. What happened in 2001 partly validated the prediction. Within 10 months, there were 2 great continental earthquakes in south Asia, i.e., the M7.8 event in India on Jan 26 and M8.1 event in China on Nov. 14, 2001, which are the largest earthquake in the past 50 years both for India and China. No records for two great earthquakes in Asia within so short time interval. We should speculate the following aspects from the two incidences: The influence of the fallacious deployment of seismic stations on the fine location and focal mechanism determination of strong earthquakes must be affronted. It is very important to introduce comparative seismology research to seismic hazard analysis and earthquake prediction research. The improvement or changes in real-time prediction of strong earthquakes with precursors is urged. Methods need to be refreshed to protect environment and historical relics in earthquake-prone areas.