河流成因储层沉积微相与剩余油分布及形成机理研究

Based on the study of fluvial sandstone reservoir in upper of Guantao group in Gudao and Gudong oilfields, this paper first introduces A.D.Miall's(1996a) architectural-element analysis method that was summarized from ground outcrop scale into the reservoir formation research of the study area, more subtly divides sedimentary microfacies and establishes sedimentary model of research area.on this base, this paper summarizes the laws of residual oil distribution of fluvial formation and the control effect of sedimentary microfacies to residual oil distribution, and reveals residual oil formation mechanism. These results have been applied to residual oil production, and the economic effect is good. This paper will be useful for residual oil research and production and enhancement of oil recovery in similar reservoir. The major conclusions of this paper are as follows. 1. Using the architectural-element analysis method to the core data, a interfacial division scheme of the first to the dixth scale is established for the studied fluvial formation. 2.Seven architectural-elements are divided in upper of Guantao group of study area. The sandstone group 5～1+2 of Neogene upper Gutao group belongs to high sinuous fine grain meandering river, and the sandstone group 6 is sandy braided river. 3. Inter layer, the residual oil saturation of "non-main layer" is higher than "main layer", but the residual recoverable reserve of former is larger. Therefore, "main layer" is the main body of residual oil distribution. The upper and middle part of inner layer has lower permeability and strong seeping resistance. Addition to gravity effect in process of driving, its driving efficiency is low; residual oil saturation is high. Because of controlling of inside non-permeable interlayer or sedimentary construction, the residual oil saturation of non-driving or lower driving efficiency position also is high. On plane, the position of high residual oil saturation mostly is at element LV, CS, CH (FF), FF etc, Which has lower porosity and permeability, as well as lens sand-body and sand-body edge that is not controlled by well-net, non-perfect area of injection and production, lower press difference resort area of inter-well diffiuent-line and shelter from fault, local high position of small structure. 4.Microscopic residual oil mainly includes the non-moved oil in the structure of fine pore network, oil in fine pore and path, oil segment in pore and path vertical to flow direction, oil spot or oil film in big pore, residual oil in non-connective pore. 5.The most essential and internal controlling factor of fluvial formation residual oil distribution is sedimentary microfacies. Status of injection and production is the exterior controlling factor of residual oil distribution. 6. The controlling effect of formation sedimentary microfacies to residual oil distribution indicates inter-layer vertical sedimentary facies change in scale of injection and production layer-series, planar sedimentary face change and inner-layer vertical sedimentary rhythm and interbed in single layer to residual oil distribution. 7. It is difficult to clear up the inter-layer difference in scale of injection and production layer-series. The using status of minor layer is not good and its residual oil saturation is high relatively. It is obvious that inter-layer vertical sedimentary facies changes control inter-layer residual oil distribution at the same or similar conditions of injection and production. For fluvial formation, this vertical sedimentary facies change mainly is positive gyration. Namely, from down to top, channel sediment (element CHL, LA) changes into over-bank sediment (element LV, CR, CS). 8. In water-injection developing process of transverse connecting fluvial sandstone oil formation, injection water always comes into channel nearby, and breaks through along channel and orientation of high pressure gradient, does not expand into side of channel until pressure gradient of channel orientation changes into low. It brings about that water-driving status of over-bank sedimentary element formation (LV, CR, CS) is not good, residual oil saturation is high. In non-connective abandoned channel element (CH) formation with channel, because this reverse is difficult to control by injection and production well-series, its using status is not good, even terribly not good, residual oil is enrichment. 9. The rhythm and sedimentary structure, sedimentary facies change in single sand body brings about vertical changes of formation character, growth character of inner layer interbed. These are important factor of controlling and affecting vertical water spread volume and inner layer residual oil forming and distribution in single sand body. Positive rhythm, is the principal part of fluvial sandstone inner layer sedimentary rhythm. Namely, from down to upside, rock grain granularity changes from coarse to fine, seeping ability changes from strong to feebleness. It brings about that water-driving status of inner layer upside is not good, residual oil saturation is high. Inner layer interbed has different degree affecting and controlling effect to seeping of oil and water. Its affecting degree lies on interbed thickness, extending scale, position, and jeted segment of production or injection well. The effect of interbed at upside of oil formation to oil and water seeping is less; the effect of interbed at middle of oil formation to oil and water seeping is more. 10. Indoor experiment and research indicate that wettability, permeability step, vertical permeability, position of Kmax and ratio of oil viscousity and water viscousity all have great effect on the water-driving recovery ratio. 11. Microscopic residual oil distribution is affected and controlled by formation pore network structure, pressure field distribution, and oil characteristic. 12.The residual oil forming mechanism: the over-bank sedimentary element and upper part of a positive rhythm sandstone have fine pore and throat network, permeability is low, displacement pressure of pore and throat is high. The water-driving power usually falls short of displacement pressure that brings about injection water does not spread into these pore and throat network, thereby immovable oil area, namely residual oil, is formed. At underside of channel sedimentary element and positive rhythm sandstone, porosity and permeability is relatively high, connecting degree of pore and throat is high, displacement pressure of pore and throat is low. Thereby injection water is easy to enter into pore and throat, driving oil in them. Because the pore space is irregular, the surface of pore wall is coarse and non-flat. That the oil locate on concave hole of pore wall and the dead angle of pore, and the oil attaches on surface of pore wall by surface tension, are difficult to be peeled off, becoming water-driving residual oil (remaining oil). On the other hand, Because flowing section lessens, flowing resistance increase, action of capillary fore, or seeping speed decreases at process of transfer at pass narrow throat path in the course carried by driving water. The "oil drop", "oil bead", or "oil segment" peeled off by driving water is difficult to carry and to drive out by water at less pressure difference. Thereby they are enclosed in pore to form discontinuous residual oil. 13.This results described above have been applied in nine develop blocks of Gudao and Gudong oilfield. Its applying effect is marked through local injection production adjustment, deploying replacement well, repair hole, replacement envelop, block off water and profile control etc. Relative method and technology can be applied to other oil production area of Shengli oilfield, and obtain better economic and societal effect.

准噶尔盆地中拐东北斜坡油气成藏条件研究——西部压性盆地含油气系统理论的研究与实践

This paper builds the model of oil accumulation and achieves the prediction of exploration goal. It uses multiple subject means, the ways of synthetic research and the viewpoint of analyzing genesis, with the academic guidance of sedimentology, structural geology, petroleum geology and geochemistry, the basis of strata sequence frame and structural frame, the frame of "four history" - the burying history, the structural history, the filling history and the evolving history of oil, the masterstroke of hydrocarbon's generation, migration and accumulation, the aim of revealing the genetic relation between mature source rock and oil reservoir in space and time. Some achievements and viewpoints in this study are following. 1. It is proposed that the structural evolution in this area had many periods, and the structural movement of the Xiazijie group telophase formed the structural pattern for the first time. 2. The character of strata sequence in this area is divided by the character of episodic cycle firstly. The study of dividing the facies of single well and the facies of well tie is based on the data of single well. The character of sedimentary facies is con-structed initially. 3. It is believed that Jiamuhe group is the main source rock, which can supply considerable oil and gas resources for the first time. Some criterions of source rock such as the type ,the abundance in Jiamuhe group are analysed. Using the thermal history of source rock, we drawn a conclusion that the original type of source rock in Jiamuhe group is II_1-III, and the abundance achived the level of good source rock, and this set of source rock had contributed to this area. 4. The reservoir strata in this area are assessed and analysed with the reservoir evaluation. There are multi-type reservoirs, such as volcanic lava facies, sedimentary clast facies, continental belch facies. The physical property in reservoir strata is characterized by low porosity and low permeability. The study of diagenetic stage show that the diageneses in Jiamuhe group is A-Bsubage, and the reservoir room is mainly secondary corroded hollow and cleft. 5. The synthetic research on oil system in Jiamuhe group is made for the first time. The type of petroleum system is divided , and we consider that the petroleum system of Jiamuhe group is at the reliable rank. There are two critical time in oil accumulation through studying the critical time of oil accumulation : the early generation of hydrocarbon is oil, and the later is gas. 6. The mechanism of accumulation is analysed. We consider that the accu-mulation of oil in this area has many periods, and the early generated hydrocarbon is expeled by the later , and formed the character of zonal distribution in planar. 7. A bran-new model of oil and gas is proposed. Beneficial enrichment area of oil and gas is analyzed, which can be divided into three sections: Section I can be divided into two sections: I_1 and I_2. The lower subgroup of Jiamuhe is covered by the triassic layer of I_1 section. Fault zone and near the foot wall of fault are charactered with thick phase belt. Then the cover capability in this area is relatively poor, oil can migrate into triassic layer by vertical or lateral migration , and forms I_1 Kelamayi triassic oil pool consequently. The lower subgroup of Jiamuhe is covered by the triassic layer of I_2 section ,which is charactered with thin phase belt. Then the cover capability in this area is relatively good, and forms I_1 Kelamayi triassic oil pool consequently. Section II can be divided into two sections: II_1-I_(I~2). The cover of Jiamuhe group in section II_1 is the low resistivity segment in Wuerhe group, which has thin lithology and poor porosity and permeability. Oil and gas in Jiamuhe group can be covered to form beneficial accumulation area. There are some wells in this area, such as Ke 007 well, 561 well. The thick phase belt layer of Wuerhe high resistivity segment in section II_2 has unconformable relation with Jiamuhe group. The cover ability of the high resistivity segment is poor, petroleum in Jiamuhe can migrate into Wuerhe layer vertically. This area is the beneficial area for accumulating petroleum in Wuerhe layer. there are some wells in this area, such as Ke 75 well, Ke 76 well, Ke 77 well, Ke 78 well, Ke 79 well. Section III can also be divided into two sections: III_1 and III_2. Wuerhe group in section III_1 has unconformable relation with Jiamuhe group. There is thick lithology and poor cover in Wuerhe group, but the strata sequence evolution character of upper subgroup in Jiamuhe group has determined that it has lateral and vertical cover ability. thus, this area is petroleum abundant belt of jiamuhe group, which has the trap. Section III_2 is an area controled by wedgeout of Fengcheng group, Fengcheng group in this area has quite thick lithology so that It has beneficial resevoir phase belt. It can accumulate oil in itself or accept some oil in Jiamuhe group. Jiamuhe group has some oil accumulation condition in this area. Thus, section III_2 is jiamuhe-Fengcheng multiple petroleum accumulation belt, such as Ke 80 well. 8. The goal of exploration is suggested: Depositional trap or combination trap is the important aspect in later exploration. Both types of traps are the goal of the next drilling: Fault block trap in the east of 576 well and the NO. 2 fault block trap in the north of Ke 102 well It is suggested that we should study the law of oil and gas in Jiamuhe group and enhance the study of combination in forming reservoir and trap scale. We do some lithology forecast and reservoir diatropic forecast in order to know the area of oil and gas.

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

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.

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

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.

接收函数方法与中国及邻区上地幔间断面的研究

The receiver function method applied in researching the discontinuities in upper mantle was systematically studied in this paper. Using the theoretical receiver functions, the characteristics of P410S and P660S phases were analyzed, and the influencing factors for detection of these phases were discussed. The stability of receiver function was studied, and a new computational method of receiver function, RFSSMS (Receiver Function of Stack and Smooth of Multi seismic-records at a Single station), was put forward. We built initial reference velocity model for the media beneath each of 18 seismic stations respectively; then estimated the buried depths of 410-km and 660-km discontinuities(simply marked as '410' and '660') under the stations by using the arrive time differences of P410S and P660S with P. We developed a new receiver function inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion), to obtain the whole crust and upper mantle velocity structure and the depths of discontinuities beneath a station. The major works and results could be summarized as follows: (1) By analysis of the theoretical receiver functions with different velocity models and different ray parameters, we obtain the knowledge: The amplitudes of P410S and P660S phases are decreasing with the increasing of epicentral distance A , and the arrival time differences of these phases with P are shorter as A is longer. The multiple refracted and/or reflected waves yielded on Moho and the discontinuities in the crust interfere the identification of P410S. If existing LVZ under the lithosphere, some multiple waves caused by LVZ will interfere the identification of P410S. The multiple waves produced by discontinuity lied near 120km depth will mix with P410s phase in some range of epicentral distance; and the multiple waves concerned with the discontinuity lied near 210km depth will interfere the identification of P660S. The epicentral distance for P4i0s identification is limited, the upper limit is 80° . The identification of P660S is not restricted by the epicenter distance obviously. The identification of P410S and P6gos in the theoretical receiver functions is interfered weakly from the seismic wave attenuation caused by the media absorption if the Q value in a reasonable range. (2) The stability of receiver function was studied by using synthetic seismograms with different kind of noise. The results show that on the condition of high signal-noise-ratio of seismic records, the high frequency background noise and the low frequency microseism noise do not influence the calculating result of receiver function. But the media "scattering noise" influence the stability of receiver function. When the scattering effect reach some level, the identification of P4iOs and P66os is difficult in single receiver function which is yielded from only one seismic record. We provided a new method to calculate receiver function, that is, with a group of earthquake records, stacking the R and Z components respectively in the frequency domain, and weighted smooth the stacked Z component, then compute the complex spectrum ratio of R to Z. This method can improve the stability of receiver function and protrude the P4i0s and P66os in the receiver function curves. (3) 263 receiver functions were provided from 1364 three component broadband seismograms recorded at 18 stations in China and adjacent areas for the tele-earthquakes. The observed arrival time differences of P410S and P660S with P were obtained in these receiver functions. The initial velocity model for every station was built according to the prior research results. The buried depths of '410' and '660' under a station were acquired by the way of adjusting the depths of these two discontinuities in the initial velocity model until the theoretical arrival time differences of P410S and P660S with P well conformed to the observed. The results show an obvious lateral heterogeneity of buried depths of ' 410' and (660' . The depth of '410' is shallower beneath BJI, XAN, LZH and ENH, but deeper under QIZ and CHTO, and the average is 403km . The average depth of '660' is 663km, deeper under MDJ and MAJO, but shallower under QIZ and HYB. (4) For inversing the whole crust and upper mantle velocity structure, a new inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion) has- been developed here. The media beneath a station is divided into segments, then the velocity structure is inversed from receiver function from surface to deep successively. Using PGARFI, the multi reflection / refraction phases of shallower discontinuities are isolated from the first order refraction transform phase of deep discontinuity. The genetic algorithm with floating-point coding was used hi the inversion of every segment, and arithmetical crossover and non-uniform mutation technologies were employed in the genetic optimization. 10 independent inversions are completed for every segment, and 50 most excellent velocity models are selected according to the priority of fitness from all models produced in the inversion process. The final velocity structure of every segment is obtained from the weighted average of these 50 models. Before inversion, a wide range of velocity variation with depth and depth range of the main discontinuities are given according to priori knowledge. PGARFI was verified with numerical test and applied in the inversion of the velocity structure beneath HIA station down to 700km depth.

中强地震震源的破裂过程

With the large developments of the seismic sources theory, computing technologies and survey instruments, we can model and rebuild the rupture process of earthquakes more realistically. On which earthquake sources' properties and tectonic activities law are realized more clearly. The researches in this domain have been done in this paper as follows. Based on the generalized ray method, expressions for displacement on the surface of a half-space due to an arbitrary oriented shear and tensile dislocation are also obtained. Kinematically, fault-normal motion is equivalent to tensile faulting. There is some evidence that such motion occurs in many earthquakes. The expressions for static displacements on the surface of a layered half-space due to static point moment tensor source are given in terms of the generalized reflection and transmission coefficient matrix method. The validity and precision of the new method is illustrated by comparing the consistency of our results with the analytical solution given by Okada's code employing same point source and homogenous half-space model. The computed vertical ground displacement using the moment tensor solution of the Lanchang_Gengma earthquake displays considerable difference with that of a double couple component .The effect of a soft layer at the top of the homogenous half-space on a shallow normal-faulting earthquake is also analyzed. Our results show that more seismic information would be obtained utilizing seismic moment tensor source and layered half-space model. The rupture process of 1999 Chi-Chi, Taiwan, earthquake investigated by using co-seismic surface displacement GPS observations and far field P-wave records. In according to the tectonic analysis and distributions of aftershock, we introduce a three-segment bending fault planes into our model. Both elastic half-space models and layered-earth models to invert the distribution of co-seismic slip along the Chi-Chi earthquake rupture. The results indicate that the shear slip model can not fit horizontal and vertical co-seismic displacements together, unless we add the fault-normal motion (tensile component) in inversions. And then, the Chi Chi earthquake rupture process was obtained by inversion using the seismograms and GPS observations. Fault normal motions determined by inversion, concentrate on the shallow northern bending fault from Fengyuan to Shuangji where the surface earthquake ruptures reveal more complexity and the developed flexural slip folding structures than the other portions of the rupture zone For understanding the perturbation of surface displacements caused by near-surface complex structures, We have taken a numeric test to synthesize and inverse the surface displacements for a pop-up structure that is composed of a main thrust and a back thrust. Our result indicates that the pop-up structure, the typical shallow complex rupture that occurred in the northern bending fault zone form Fengyuan to Shuangji, can be modeled better by a thrust fault added negative tensile component than by a simple thrust fault. We interpret the negative tensile distributions, that concentrate on the shallow northern bending fault from Fengyuan to Shuangji, as a the synthetic effect including the complexities of property and geometry of rupture. The earthquake rupture process also reveal the more spatial and temporal complexities form Fenyuan to SHuangji. According to the three-components teleseismic records, the S-wave velocity structure beneath the 59 teleseismic stations of Taiwan obtained by using the transform function method and the SA techniques. The integrated results, the 3D crustal structure of Taiwan reveal that the thickest part of crustal local in the western Central Range. This conclusion is consistent with the result form the Bouguer gravity anomaly. The orogenic evolution of Taiwan is young period, and the developing foot of Central Range dose not in static balancing. The crustal of Taiwan stays in the course of dynamic equilibrium. The rupture process of 2003)2,24,Jiashi, Xinjiang earthquake was estimated by the finite fault model using far field broadband P wave records of CDSN and IRIS. The results indicate that the earthquake focal is north dip trust fault including some left-lateral strike slip. The focal mechanism of this earthquake is different form that of earthquakes occurred in 1997 and 1998, but similar to that of 1996, Artux, Xinjiang earthquake. We interpreted that the earthquake caused trust fault due to the Tarim basin pushing northward and orogeny of Tianshan mountain. In the end, give a brief of future research subject: Building the Real Time Distribute System for rupture process of Large Earthquakes Based on Internet.

安多地区早中侏罗世岩浆岩岩石地球化学和同位素年代学研究及其对构造演化的制约

Anduo area is located in the Central Tibet, the middle segment of the Bangonghu-Nujiang suture. Anduo Block is the northern part of Lhasa terrane. The relationships among the different geological bodies were determined during the 1: 250000 regional geological surveying. Petrography, petrologic geochemistry, isotopic geochemistry and geochronology of igneous rocks from the suture and granitoids from Anduo Block were analyzed systematically as a whole for the first time. Then, their tectonic setting and history are discussed.Anduo ophiolitic melange consists of metamorphic peridotites, cumulates, plagiogranites, sheeted dykes swarm, pillow lava and radiolarian cherts. The concentration of Cr and Ni in the metamorphic peridotites is very high, with Mg# about 0.94 ~ 0.97, higher 87Sr/86Sr and Pb isotopic ratios, and lower 143Nd/i44Nd ratio. LREE is enriched relative to HREE and positive Eu anomaly is very clear. The REE distribution curve is U shape. Nb and Ta anomalies from cumulate gabbro and sheeted dyke swarm are not clear, while that are slightly negative from pillow lava. Plagiogranite belongs to strong calc-alkaline series with high Si, middle Al, low Fe, Mg and low K contents. Eu anomaly (~ 1.23) from plagiogranites is slightly positive. The character of all components of ophiolite is similar to that of the MORB, while to some extent the ophiolite was influenced by crustal material. Anduo ophiolite formed in a mature back-arc basin. Additionally, intermediate acidity volcanic rocks within Anduo phiolite melange are island arc calc-alkline rocks related to ocean subduction.The early-middle Jurassic plutonic rocks are tonalite, granodiorite bearing-phenocryst, magaporphyritic hornblende monzogranite, magaporphyritic monzogranite, monzogranite bearing-phenocryst and syenogranite in turn. They belong to calc-alkaline series which developed from middle K to high K series temporally. REE distribution curves of all plutonic rocks are similar and parallel to each other. SREE and negative Eu anomaly values decrease. In the multi-element spider diagram, the curves of different plutons are similar to each other, but troughs of Nb, Sr, P and Ti from young plutons become more evident. This suggests that thereare some closely petrogenetic affinities among plutonic rocks which make up amagma plutonism cycle of the early-middle Jurassic. Magma source is mainly crustal,but abundant mafic microgranular enclaves within granitoids indicate that crastalmagma should be mixed with mantle-derived magma and the mantle-derived magmadecreased subsequently. Tonalite has features of I-type granite, magaporphyriticmonzogranite is transition type, and monzogranite bearing-phenocryst is S-typegranite. The characteristic of granitoids from Anduo Block suggest that the formingtectonic setting is active continental margin.Reliable zircon U-Pb SHRIMP ages are obtained in the study area firstly. Plagiogranite from the Anduo ophiolite of the Bangonghu-Nujiang suture is 175.1 Ma, and granitoids from Anduo Block is 172.6-185.4 Ma. Additionally, plagioclase from the plagiogranite dates a 40Ar/39Ar age of 144 Ma, while biotite and hornblend from granitoids of Anduo Block give a 163-165 Ma.Similar cooling ages of plagiogranite from the Anduo ophiolitic melange and granitoids from Anduo Block and the spatial distribution of the ophiolitic rocks between Anduo, Naqu, and Shainzha area suggest that bilateral subduction of the Bangonghu-Nujiang oceanic basin took place in the early-middle Jurassic. During this subduction, Anduo ophiolitic rocks were related to north subduction of the Bangonghu-Nujiang oceanic basin and Anduo back-arc basin spreading, while granitoids from Anduo Block were related to south subduction.

大港油田孔南地区孔一段储层流动单元研究

The main reservoir type in the south of Dagang Oilfield is alluvial reservoir. In this paper, the reservoir structure model and the distribution of connected body and flow barrier were built on base of the study of high-resolution sequential stratigraphic skeleton and fine sedimentary microfacies on level of single sandbody. Utilizing the static and dynamic data synthetically and carrying out the comparision of the classification method for reservoir flow unit in different reservoir, the criterion, which can be used to classify the flow unit in first section of Kongdian formation of Kongnan area, was defined. The qualitative method of well-to-well correlation and the quantitative method of conditional simulation using multiple data are adopted to disclose the oil and water moving regulation in different flow unit and the distribution rule of remaining oil by physical simulation measure. A set of flow unit study method was formed that is suit for the Dagang Oilfield on account of the remaining oil production according to the flow unit. Several outstanding progresses was obtained in the following aspects:It is considered that the reservoir structure of Zao V iow oil group- Zao Vup4 layerand are jigsaw-puzzled reservoir, while ZaoVup3-ZaoVupi layers are labyrinth reservoir,which are studied on base of high-resolution sequential stratigraphic skeleton on the levelof single sandbody in first section of Kongdian formation of Kongnan area and accordingto the study of fine sedimentary microfacies and fault sealeing.When classifying the flow unit, only permeability is the basic parameter using thestatic and dynamic data and, and also different parameters should be chose or deleted, suchas porosity, effective thickness, fluid viscosity and so on, because of the weak or stronginterlayer heterogeneous and the difference of interlayer crude oil character.The method of building predicting-model of flow unit was proposed. This methodis according to the theories of reservoir sedimentology and high-resolution sequencestratigraphic and adopts the quantitative method of well-to well correlation and the quantitative method of stochastic simulation using integrateddense well data. Finally the 3-D predicting-model of flow unit and the interlay er distribution model in flow unit were built which are for alluvial fan and fan delta fades in first section of Kongdian formation of Kongnan area, and nine genetic model of flow unit of alluvial environment that spread in the space were proposed.(4) Difference of reservoir microscopic pore configuration in various flow units and difference of flow capability and oil displacement effect were demonstrated through the physical experiments such as nuclear magnetic resonance (NMR), constant rate mercury penetration, flow simulation and so on. The distribution of remaining oil in this area was predicted combining the dynamic data and numerical modeling based on the flow unit. Remaining oil production measure was brought up by the clue of flow unit during the medium and late course of the oilfield development.

单（多）井抽灌对浅部地温场的影响研究

The technique of energy extraction using groundwater source heat pumps, as a sustainable way of low-grade thermal energy utilization, has widely been used since mid-1990's. Based on the basic theories of groundwater flow and heat transfer and by employing two analytic models, the relationship of the thermal breakthrough time for a production well with the effect factors involved is analyzed and the impact of heat transfer by means of conduction and convection, under different groundwater velocity conditions, on geo-temperature field is discussed.A mathematical model, coupling the equations for groundwater flow with those for heat transfer, was developed. The impact of energy mining using a single well system of supplying and returning water on geo-temperature field under different hydrogeological conditions, well structures, withdraw-and-reinjection rates, and natural groundwater flow velocities was quantitatively simulated using the finite difference simulator HST3D. Theoretical analyses of the simulated results were also made. The simulated results of the single well system indicate that neither the permeability nor the porosity of a homogeneous aquifer has significant effect on the temperature of the production segment provided that the production and injection capability of each well in the aquifers involved can meet the designed value. If there exists a lower permeable interlayer, compared with the main aquifer, between the production and injection segments, the temperature changes of the production segment will decrease. The thicker the interlayer and the lower the interlayer permeability, the longer the thermal breakthrough time of the production segment and the smaller the temperature changes of the production segment. According to the above modeling, it can also be found that with the increase of the aquifer thickness, the distance between the production and injection screens, and/or the regional groundwater flow velocity, and/or the decrease of the production-and-reinjection rate, the temperature changes of the production segment decline. For an aquifer of a constant thickness, continuously increase the screen lengths of production and injection segments may lead to the decrease of the distance between the production and injection screens, and the temperature changes of the production segment will increase, consequently.According to the simulation results of the single well system, the parameters, that can cause significant influence on heat transfer as well as geo-temperature field, were chosen for doublet system simulation. It is indicated that the temperature changes of the pumping well will decrease as the aquifer thickness, the distance between the well pair and/or the screen lengths of the doublet increase. In the case of a low permeable interlayer embedding in the main aquifer, if the screens of the pumping and the injection wells are installed respectively below and above the interlayer, the temperature changes of the pumping well will be smaller than that without the interlay. The lower the permeability of the interlayer, the smaller the temperature changes. The simulation results also indicate that the lower the pumping-and-reinjection rate, the greater the temperature changes of the pumping well. It can also be found that if the producer and the injector are chosen reasonably, the temperature changes of the pumping well will decline as the regional groundwater flow velocity increases. Compared with the case that the groundwater flow direction is perpendicular to the well pair, if the regional flow is directed from the pumping well to the injection well, the temperature changes of the pumping well is relatively smaller.Based on the above simulation study, a case history was conducted using the data from an operating system in Beijing. By means of the conceptual model and the mathematical model, a 3-D simulation model was developed and the hydrogeological parameters and the thermal properties were calibrated. The calibrated model was used to predict the evolution of the geo-temperature field for the next five years. The simulation results indicate that the calibrated model can represent the hydrogeological conditions and the nature of the aquifers. It can also be found that the temperature fronts in high permeable aquifers move very fast and the radiuses of temperature influence are large. Comparatively, the temperature changes in clay layers are smaller and there is an obvious lag of the temperature changes. According to the current energy mining load, the temperature of the pumping wells will increase by 0.7°C at the end of the next five years. The above case study may provide reliable base for the scientific management of the operating system studied.

甘肃金川II号岩体岩浆演化及铜镍硫化物成矿过程探讨

金川矿床是超大型岩浆 Ni-Cu-PGE 硫化物矿床，它赋存有世界第三的镍储量（5.45×106 t，Ni 的平均品位为 1.2%）。该矿床位于华北板块西侧阿拉善地块西南边缘龙首山隆起中。本论文以金川矿区的 II 号岩体为研究对象，在系统整理前人工作的基础上，通过深入细致的野外地质考察和系统采样，选取有代表性的样品，进行详细的显微镜鉴定，并运用元素地球化学等手段，对 II 号岩体的源区特征、岩浆演化过程及其赋存的⑴、⑵号矿体的成矿机制等主要矿床学问题进行了系统成因研究，并尝试提出了二次硫化物熔离成矿模式。本论文主要取得以下几点结论性成果： 1）岩浆起源：起源于上地幔尖晶石－石榴石二辉橄榄岩过渡地带的下部。通过上地幔源岩的柱状部分熔融模式（≥25%）或三角部分熔融模式（10～15%），形成富PGE的S不饱和原始岩浆。 2）地壳同化混染：地壳物质的同化混染主要是由早期下地壳物质的同化混染和晚期上地壳物质的同化混染组成的多阶段同化混染，并且其成分为富钙质围岩。 3）结晶分异演化：金川II号岩体作为橄榄石堆积相，其橄榄石的结晶主要发生在上地壳深部岩浆房，而橄榄石与熔浆的反应及辉石和斜长石的结晶则主要发生在浅部岩浆房，既现今岩体所在位置。 4）S 的饱和机制：岩浆通道内下地壳物质的同化混染，导致岩浆硫饱和，并于岩浆通道中发生了第一次硫化物熔离；上地壳深部岩浆房中地壳物质的同化混染和橄榄石的分离结晶，促使硫饱和的母岩浆发生了第二次硫化物熔离。 5）深部岩浆房的演化：上地壳深部岩浆房中橄榄石分离结晶与硫化物熔离之后，在重力分异作用下，形成了自上而下的硅酸盐岩浆、橄榄石“晶粥”、含硫化物熔体的橄榄石“晶粥”以及矿浆的四层分层格局。同时，底部部分硫化物熔体发生了Pt合金（Pt-Fe）的脱离。 6）侵入期次：上地壳深部岩浆房中，母岩浆分异演化早期，其顶部含有尚未完全成长橄榄石颗粒的“晶粥”在构造应力挤压下首先上侵、成岩形成金川II号岩体顶部、呈“上悬体”的中细粒含辉橄榄岩；分层格局形成之后，橄榄石“晶 VI 粥”、含硫化物熔体的橄榄石“晶粥”以及矿浆依次上侵成岩、成矿，分别形成金川 II 号岩体上部的二辉橄榄岩相、⑴、⑵号矿体的主体以及⑵号矿体底部的块状矿石。

汉语句子和语篇韵律层级边界的认知加工及其ERP效应

The aim of the present study was to explore whether the CPS (Closure Positive Shift) which reflected prosodic processing will be elicited when listeners perceived different hierarchical prosodic boundaries in Chinese sentence and discourse (Quatrain). In addiction, the similarity and difference in amplitude, onset latency and scalp distribution between these CPS were investigated. The nature of the CPS and its relationship to acoustic parameters was also explored systematically. The main results and conclusions of the present study were: (1) Phonological phrase boundaries and intonational phrase boundaries in Chinese sentences both elicited the CPS; however, phonological word boundaries can't evoke it. The CPS induced by phonological phrase boundaries was earlier than the one related to intonational phrase boundaries in onset latency, and the amplitude was also somewhat lower. When the pauses in the vicinity of these two boundaries were removed, the onset latency difference disappeared while amplitude in the new conditions was also lower. This indicates that whenever listeners segment sentence into phrases, the CPS will be elicited. Besides, pause was not the decisive factor to elicit the CPS, but can modify its onset latency and amplitude effectively. (2) The different hierarchical prosodic boundaries in seven character quatrain including phonological phrase boundaries, intonational phrase boundaries and sentence pair boundaries elicited the CPS respectively. Furthermore, just like in the sentence level, onset latency of the CPS induced by the prosodic boundaries in the discourse was also influenced by the length of pause: the shorter the pause was, the earlier the onset latency. For the comparison between the CPS evoked by the same and different hierarchical prosodic boundaries, its amplitude was influenced by the extent to which prosodic representations were activated. Thus, the condition of the CPS elicitation was extended to the prosodic bounaries in discourse, and further indicated that it was influenced by acoustic parameters. (3) No matter what task the participants completed, just like word detection or rythem matching task, the CPS will be evoked. However, its amplitude was larger in the anterior region, when listeners completed the word detection task which needed more attention and higher load of working memory. The present result indicated that the elicitation of the CPS was not influenced by the task the participants completed, but different task influence its scalp distribution. (4) The final syllable of the sentence and quatrain can't elicit the CPS, but a P300-like positive component. Although the scalp distribution was similar to the CPS, it was much higher in amplitude. The present result suggested that only the prosodic boundaries reflecting not only the closure of the former prosodic unit but also integrating the later one will elicit the CPS.