特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

三维非均匀各向异性介质中的纵波走时反演研究

Many observations show that seismic anisotropy is very common in the crust and upper mantle of the Earth. Seismic anisotropy can provide some clue about the changing and transporting process inside the earth. in recent years, abundant earthquake travel time data are accumulated, computers become more powerful, and these make the inversion of earthquake travel time data practical. In this thesis we studied the theory of elastic wave in anisotropic media, some formule for travel time inversion were derived. We present an iterative procedure to determine 21 elastic parameters from qP wave travel times. No a priori assumptions about heterogeneity and anisotropy of the model are made. The procedure is suitable for the case when we know nothing about the symmetry of anisotropy of the media, as well as for the case of earthquake travel time inversion which may contain various symmetry of anisotropy. The procedure is tested with a synthetic multiple-source offset VSP experiment. The results proved that the formulae are correct, and the procedure is practical. The results and the related theory indicate that the anisotropic inversion needs more rays than isotropic case. For a 2-D weak anisotropic (WA) medium, we need at least 5 rays in different directions to retrieve the elastic parameters on one grid point, and for a 3-D WA medium we need at least 15 rays in different directions to retrieve the elastic parameters on one grid point. The results also indicate that the starting background velocity has no influence on the final results, at least for the model we specified. Our results also show that insufficient illumination coverage will slow down the convergence rate, and make the results more sensitive to noise. We apply the procedure to a set of field travel time data. The data is from an artificial seismic observation. This observation is for locating micro-seismic events around a tunnel, its purpose is to find out if the digging process and the stress condition around the tunnel can generate micro-cracks. The size of this area is around 100m. The anisotropy derived from qP travel times is the same as the anisotropy showed by apparent velocities, and is also consistent with the anisotropy derived from S-wave splitting phenomena.

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

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

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.

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

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.

古潜山油藏非均质三维地质建模

The discovery of the highly productive Renqiu buried hill reservoir in Bohai Bay Basin in 1975 started the high tide of finding buried hill reservoirs in China and their research. As the advance of E&P technologies, the study of buried hill reservoir in China had a qualitative leap. The reservoir description and some other aspects of development have reached or approached to the international leading level. However, some core techniques for reservoir study such as structure & faulting system study, formation prediction and connection study and heterogeneous model's construction could not completely carry out the quantitative or accurate reservoir description, e. g. the areal distribution of porosity, permeability and oil saturation. Especially, the modeling for reservoir simulation is still wandering in the stage of simplicity. The inaccurate understanding of geology could not derive 3D heterogeneous geological model that can reveal the actual underground situation thus could not design practical and feasible oilfield development plan. Therefore, the problems of low oil recovery rate, low recovery factor and poor development effectiveness have not been solved. The poor connection of the reservoir determined that waterflooding could not get good development effect and the production had to depend on the reservoir elastic energy, and this will bring big difficulty for development modification and improvement of oil recovery. This study formed a series of techniques for heterogeneous model research that can be used to construct heterogeneous model consistent with the reservoir geology. Thus the development effectiveness, success ratio of drilling and percent of producing reserves can be enhanced. This study can make the development of buried hill reservoir be of high recovery rate and high effect. The achievements of this study are as follows: 1. Evaluated the resources, summarized the geological characteristics and carried out the reservoir classification of the buried hill reservoirs in Shengli petroliferous area; 2. Established the markers for stratigraphical correlation and formed the correlation method for complex buried hill reservoirs; 3. Analyzed the structural features of the buried hill reservoirs, finished the structure interpretation and study of faulting system using synthetic seismograms, horizontal slices and coherent analysis, and clarified structural development history of the buried hill reservoirs in Shengli petroliferous area; 4. Determined the 3 classes and 7 types of pore space and the main pore space type, the logging response characteristics and the FMI logging identified difference between artificial and natural fractures by the comprehensive usage of core analysis, other lab analyses, conventional logging, FMI logging and CMR logging; 5. Determined the factors controlled the growth of the fractures, vugs and cavities, proposed the main formation prediction method for buried hill reservoir and analyzed their technical principium and applicability, and formed the seismic method and process for buried hill reservoir description; 6. Established the reserve calculation method for buried hill reservoirs, i. e. the reserves of fractures and matrix are calculated separately; the recoverable reserves are calculated by decline method and are classified by the SPE criteria; 7. Studied restraining barriers and the sealing of the faults thus clarified the oil-bearing formations of the buried hill reservoirs, and verified the multiple reservoir forming theory; 8. Formed reasonable procedure of buried hill reservoir study; 9. Formed the 3 D modeling technology for buried hill reservoirs; 10. Studied a number of buried hill blocks on the aspects of reservoir description, reservoir engineering and development plan optimization based on the above research and the profit and social effect are remarkable.

丘陵油田三间房组储层三维地质模型和砂体预测

This paper is belonging to Chinese Petrochemical Corporation's science and technology project. Although it is difficult, it has important theoretical and practical value. The study was aimed to reveal inhomogeneity of two kinds of reservoirs of fan-shaped delta and braided river by using new theories, new methods and new technology about 3-D model building and reservoir knowledge repository throughout the world, and to build reservoir knowledge repository and 3-D geological model which would predict the type of sand body forming reason and distribution rule in order to improve exploration result in Qiuling oil fields. Multi-discipline theories such as petroleum structure geology, reservoir geology, petroleum geology, sequence geology, logging geology, geomathematics and so on are used as guide. The information of geology, seism, logging and production test is combined. Outcrop area and overlap area are combined. By making full use of computer, stable structure, reservoir geometric shape, spatial distribution and inhomogeneity of bed of interest are investigated, described and characterized. Petroleum pool 3-D static geological model of reservoir knowledge repository was built. Sand body distribution was predicted. It has guided oil development, lowed the investment and improved development benefits. Several results are achieved as follows: (1) Strata framework of Sanjianfang group in Qiuling oil field has been established. (2) Geometric shape, spatial distribution and evolve rule of two different forming reason's reservoir of fan-shaped delta and braided river of Sanjianfang group in Qiuling oil field are discussed. (3) The two kinds of reservoirs have lower pore and permeability and very strong inhomogeneity. (4) Reservoir knowledge repository of two different forming reasons has been built of Sanjianfang group, which includes 5 geological knowledge sublibrary. (5) 3-D geological model of two kinds of forming reason's reservoirs has been built. (6) That same sequence instruction a simulation and probability field were used to predict sand body of Sanjianfang group was put forward. Coincidence rate is high after production test. It shows this method has great popularity value. (7) A set of theories, methods and technologies of knowledge repository of two kinds of reservoir of braided river and fan-shaped delta and 3-D geological model building were finished. (8) A set of theories, methods and technologies of investigating, describing, characterizing and predicting two kinds of oil pool were developed. It gets noticeable economic benefit after exploration. Theory and method about extrusion basin are developed.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

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

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 foreland basin on the northern margin of the lower reach of the Yangtze river (the lower Yangtze foreland basin) is tectonically situated in the basin-mountain transitional area along the southeastern flank of the Dabie mountains. The early formation and development of the basin is closely related to the open-up of the Mian-Lue paleo-oceanic basin on the southern margin of the Central Orogenic System represented by Qinling-Dabei orogenic belt, while the tectonic evolution of the middle-late stage of the basin is mainly related to development of the Mian-Lue tectonic zone that occurred on the basis of the previous Mian-Lue paleo-suture. The foreland basin of the northern rim of the lower reach of the Yangtze river was formed during the middle-Triassic collision between the Yangtze and North China plates and experienced an evolution of occuirence-development-extinction characterized by marine facies to continental facies and continental margin to intracontinent in terms of tectonic setting.The foreland basin (T2-J2) was developed on the basis of the passive continental marginal basin on the south side of the Mian-Lue paleo-ocean and superimposed by late Jurassic-Tertiary fault basin. The tectonic setting underwent a multiple transformation of rifting-collisional clososing-tensional faulting and depression, which resulted in changes of the property for the basin and the final formation of the superposed compose basin in a fashion of 3-story-building. According to the tectonic position and evolution stages of plate collision happening on the southeastern margin of the Dabie mountains, and tectono-tratigraphic features shown by the foreland basin in its main formational period, the evolution of the foreland basin can be divided into four stages: 1) pre-orogenic passive margin (P2-Ti). As the Mian-Lue ocean commenced subduction in the late-Permian, the approaching of the Yangtze and North China plates to each other led to long-periodical and large-scale marine regression in early Triassic which was 22 Ma earlier than the global one and generated I-type mixed strata of the clastic rocks and carbonate, and I-type carbonate platform. These represent the passive stratigraphy formed before formation of the foreland basin. 2) Foreland basin on continental margin during main orogenic episode (T2.3). The stage includes the sub-stage of marine foreland basin (T2X remain basin), which formed I-type stratigrphy of carbonate tidal flat-lagoon, the sub-stage of marine-continental transition-molasse showing II-type stratigraphy of marine-continental facies lake - continental facies lake. 3) Intracontinental foreland basin during intracontinental orogeny (Ji-2)- It is characterized by continental facies coal-bearing molasses. 4) Tensional fault and depression during post-orogeny (J3-E). It formed tectono-stratigraphy post formation of the foreland basin, marking the end of the foreland evolution. Fold-thrust deformation of the lower Yangtze foreland basin mainly happened in late middle-Jurassic, forming ramp structures along the Yangtze river that display thrusting, with deformation strength weakening toward the river from both the Dabie mountains and the Jiangnan rise. This exhibits as three zones in a pattern of thick-skinned structure involved the basement of the orogenic belt to decollement thin-skinned structure of fold-thrust from north to south: thrust zone of foreland basin on northern rim of the lower reach of the Yangtze river, foreland basin zone and Jiannan compose uplift zone. Due to the superposed tensional deformation on the earlier compressional deformation, the structural geometric stratification has occurred vertically: the upper part exhibits late tensional deformation, the middle portion is characterized by ramp fault -fold deformation on the base of the Silurian decollement and weak deformation in the lower portion consisting of Silurian and Neo-Proterozoic separated by the two decollements. These portions constitutes a three-layered structural assemblage in a 3-D geometric model.From the succession of the lower reach of the Yangtze river and combined with characteristics of hydrocarbon-bearing rocks and oil-gas system, it can be seen that the succession of the continental facies foreland basin overlies the marine facies stratigraphy on the passive continental margin, which formed upper continental facies and lower marine facies hydrocarbon-bearing rock system and oil-gas forming system possessing the basic conditions for oil-gas occurrence. Among the conditions, the key for oil-gas accumulation is development and preservation of the marine hydrocarbon-bearing rocks underlying the foreland basin. The synthetic study that in the lower Yangtze foreland basin (including the Wangjiang-Qianshan basin), the generation-reservoir-cover association with the Permian marine facies hydrocarbon-bearing rocks as the critical portion can be a prospective oil-gas accumulation.Therefore, it should aim at the upper Paleozoic marine hydrocarbon-bearing rock system and oil-gas forming system in oil-gas evaluation and exploration. Also, fining excellent reservoir phase and well-preserved oil-gas accumulation units is extremely important for a breakthrough in oil-gas exploration.

滇西上芒岗卡林型金矿床地球化学特征

岩石化学与微量元素研究表明，在强经泥质灰岩中，SiO2和Au、As、Sb含量增高，CaO、CO2和MgO含量降低；在强碳酸盐化泥质灰岩中，AgO和CO2含量增高，SiO2含量降低；不同蚀变与矿化岩石的稀土含量和配分模式与未蚀变泥质灰岩基本一致。流体包裹体研究揭示，成矿溶液为富Cl^-型和富SO4^2-型，爆裂温度为183－378℃。同位素研究指示，δ^34S=(6.471-21.965)×10^-3，δ^18O=(9.3-11.6)×10^-3，δD=(-70--77)×10^-3,^206Pb=19.551-20.05,^027Pb/^204Pb=15.690-16.240,^208Pb/^204Pb=38.710-40.137；成矿流体具有岩浆水和大气降水相混合的特点，成矿物质具有壳－幔混合特征。

有机质在贵州晴隆大厂锑矿成矿过程中的作用

锑成矿中有机作用的研究国内外进行得极少，本文首先运用（有机）包裹体探讨了大厂锑成矿流体的一系列物理化学性质和主要成矿特征：1．成矿温度与盐度分别为150-180℃ 和10-12 Wt%，2．主成矿阶段水质类型属 F-SO_4Ca－K型，气相成份中甲烷含量较高(C_1/C_2～C_4 > 1)，3. 成矿介质水主要来自大气降水，δD＝-63～-66‰, δ~(18)O = -3.5～-0.6‰)。有机地球化学研究显示岭好组有机碳含量高（平均 Corg = 0.71% > 0.5%)，完全具备生油能力，是一重要的烃源层。有机质反射率（Ro）达1.87-2.93%，演化程序较高，热变沥青多呈充填状，具明显的运移特征。模拟实验证明石油对锑有较强的萃取能力，可作为 Sb 迁移的良好载体；且按萃取实验结果，有机迁移与无机迁移同等重要。对腐殖酸（质）和高硫原油作了较详细的讨论－－腐殖酸（质）在形成矿源层中有重要作用，高硫原油热裂解本身可以产生H_2S，提供部份硫；而辉锑矿的最终形成则缘于大气降水的加入，温度降低和氧化等因素综合作用的结构。总结不同阶段不同有机质的不同作用，提出了有机－无机综合成矿模式。

An industrial scale dehydration process for natural gas involving membranes

An industrial scale dehydration process based on hollow fiber membranes for lowering the dew point of natural gas is described in this paper. A pilot test with the feed flux scale of 12x10(4) Nm(3)/d was carried out. Dew points of -8 degreesC-13 degreesC at a gas transport pressure in the pipeline of 4.6M Pa and methane recovery of more than 98% were attained. The water vapor content of the product gas could be maintained around 0.01 vol% during a continuous run of about 700 hours. The effects of feed flux and operation pressure on methane recovery and water vapor content were also investigated. Additionally, some auxiliary technologies, such as a full-time engine using natural gas as fuel and the utilization of vent gas in the process, are also discussed. A small amount of the vent gas from the system was used as a fuel for an engine to drive vacuum pumps, and the heat expelled from the engine was used to warm up the natural gas feed. The whole system can be operated in a self-sustainable manner from an energy point of view, and has a relatively high efficiency in the utilization of natural gas.

A two-stage process with temperature-shift for enhanced anthocyanin production in strawberry cell suspension cultures

A two-stage process with temperature-shift has been developed to enhance the anthocyanin yield in suspension cultures of strawberry cells. The effect of the temperature-shift interval and the shift-time point was investigated for the optimization of this strategy. In this process, strawberry cells were grown at 30 degrees C (the optimum temperature for cell growth) for a certain period as the first stage, with the temperature shifted to a lower temperature for the second stage. In response to the temperature shift-down, anthocyanin synthesis was stimulated and a higher content could be achieved than that at both boundary temperatures but cell growth was suppressed. When the lower boundary temperature was decreased, cell growth was lowered and a delayed, sustained maximum anthocyanin content was achieved. Anthocyanin synthesis was strongly influenced by the shift-time point but cell growth was not. Consequently, the maximum anthocyanin content of 2.7 mg.g-fresh cell(-1) was obtained on day 9 by a temperature-shift from 30 degrees C, after 3-d culture, to 15 degrees C. The highest anthocyanin yield of 318 mg.L-1 on day 12 was achieved when the temperature was shifted from 30 degrees C, after 5-d culture, to 20 degrees C. For a global optimization of both the yield and productivity, the optimum anthocyanin yield and productivity of 272 mg.L-1 and 30.2 mg.L-1.d(-1) on day 9 were achieved by a two-stage culture with a temperature-shift from 30 degrees C after 3 d to 20 degrees C.

应急疏散语音告警信号设计参数的研究

Voice alarm plays an important role in emergency evacuation of public place, because it can provide information and instruct evacuation. This paper studied the optimization of acoustic and semantic parameters of voice alarms in emergency evacuation, so that alarm design can improve the evacuation performance. Both method of magnitude estimation and scale were implemented to investigate participants' perceived urgency of the alarms with different parameters. The results indicated that, participants evaluated the alarms with faster speech rate, with greater signal to noise ratio (SNR) and under louder noises more urgent. There was an interaction between noise level and content of voice alarm. Signals with speech rate below 4 characters / second were evaluated as non urgent at all. Intelligibility of the voice alarm was investigated by evaluating the key pointed recognition performance. The results showed that, speech rate’s effect was a marginal significance, and 7 characters / second has the highest intelligibility. It might because that the faster the signal spoken, the more attention was paid. Gender of speaker and SNR did not have a significant effect on the signals’ intelligibility. This paper also investigated impact of voice alarms' content on human behavior in emergency evacuation in a 3-D virtual reality environment. In condition of "telling the occupants what had happened and what to do", the number of participants who succeeded in evacuation was the largest. Further study, in which similar numbers of participants evacuate successfully in three conditions, indicated that the reaction time and evacuation time was the shortest in the aforesaid condition. Although one-way ANOVA shows that the difference was not significant, the results still provided some reference to the alarm design. In sum, parameters of voice alarm in emergency evacuation should be chosen to meet needs from both perceived urgency and intelligibility. Contents of the alarms should include "what had happened and what to do", and should vary according to noise levels in different public places.