柴达木盆地北缘西段油气运聚动力学研究

Hydrocarbon migration and accumulation are the important process to form reservoirs in sedimentary basins, and their researches are usually very difficult to be done in petroleum geology. In this paper, the west segment of northern margin of the Qaidam Basin was selected as study area. The concept of fault open coefficient, that combines multi-factors dealing with fault sealing, was applied to estimate semi-quantitatively the sealing characteristics of six faults which were considered controlling the hydrocarbon migration and accumulation. The data from boreholes were investigated to appraise the permeable characteristics of lithology combinations upon and beneath the unconformity surface. The result suggests that the basal conglomerates consist frequently the carriers. The data from boreholes and outcrops were collected to describe the sand carrier system. In order to eliminate the influence of inverse activities of the basin that made the formations be very steep, author adopts the phase method to build the basin models: for the steps before Pliocene the recovered true thickness maps were used to build the basin block; for the steps after Pliocene, the structure maps of today were used to build the basin block. During the modeling process, the results were calibrated by various measured data . the modeled results includes the dynamic evolvement course of trap form phase, vitrinite reflectance mature, the source rock expelled hydrocarbon intensity and fluid potential and petroleum plays. Author integrates the source rock expelled hydrocarbon intensity, fluid potential and carrier system and apply the migration technology based on percolation theory to simulate the oil and gas migration and accumulation course in the main accumulation times. The dominant pathways of oil and gas may show clearly the prospect distribution. Based on the hydrocarbon migration characteristics, the main control factors were synthesized, that including the effective source rock distribution, the match relationship of structural trap forming and hydrocarbon expelling from source rocks, the unconformity of Mesozoic and Cenozoic, the structures and the faults movement at Quaternary Finally, the author figures out the prospect plays in the study area.

北部湾盆地迈陈凹陷地震勘探方法研究

Maichen Depression lie between Leizhou Peninsula and Qiongzhou Strait. Oil and gas have been discovered in Weixinan Depression, Wushi Depression and Fushan Depression, which pertain to a same basin — North Sea Basin along with Maichen Depression.Jiangsu Oil started exploration at 2002. The first well began to drill at November, 2004 after gravity survey, electric method prospecting and 2D seismic exploration had been finished. Generating rock and hydrocarbon shows have been verified by the drilling. Low yield oil stream has been tested. And we started 3D seismic exploration at November, 2005. My thesis topic came from the actual needs of our exploration in the Maichen Depression. In the thesis, I give emphasis to analyse the own seismic geologic conditions of Maichen Depression. By real tests, we choosed the means to overcome or weaken the unfavorably impress owing to the own coditions in Maichen Depression. Finally, we obtained the usable seismic data. 1. Owing to the multiphase eruptive rock during the Quaternary Period, the near surface layers are very inhomogeneous. By simultaneous testing at same point with short refraction, uphole surveys of radial source and of surface source, the most appropriate method had been sorted out. Radial source uphole survey has been regarding the best practicable means in the complex area. Accurate surficial geology was very helpful to choosing of acquirement means and parameters. Basically the appropriate method of seismic acquirement has been built at Maichen area. 2. The seismic primary data has many, very strong and complex noise. By noise characteristic analysis in different domain, many means of denoising had been paralleled individual and joint application researched. As a result, the pre-stack multidomain joint denoise flow was the appropriate method. It can improve the seismic signal-to-noise ratio. 3. The problem of seismic static correction at Maichen Depression is very conspicuous. Many static correction methods had been tested individual and joint researched. The seismic data quality has been improved after choosing the appropriate combination of static correction flows. 4. Although the above-mentioned process are resultful, the seismic profile quality is just passable. Some reflector continuity and fault zone imagery are ambiguity. So it was the useful method to reduce the structural ambiguity during seismic interpretation that built-up geologic model in accord with real geologic character by areal structure study upon backbone seismic profiles. In the same way, traps have been assessed and drill targets have been selected.

非自喷层测试及评价解释技术研究与应用

Now low porosity and low permeability reservoir is one of the main targets of exploration for the onshore oilfields of China. Most of the reservoirs are none flowing because of bad formation percolation condition, poor gas oil ratio , low formation pressure coefficient and other factors. In the recent years, a number of domestic oilfields have carried out some research work and achieved some success on oil testing and production technology in such formation. But by now, there is still no systematic and mature technology, particularly testing technology in none flowing formation is still needed further study. Based on study the key problem of well testing and interpretation technology in none flowing formation, solve the important problems in well testing technology, continuously improve and innovate geological information acquisition technology for none flowing reservoir, accurately acquire boundary information and evaluate reservoir flow characteristics. Its wide application remarkable result has shown. The main results and cognitions obtained from research are as follows: 1. This new technology research results help solve the occurrent problems in well testing process for none flowing formations, such as small investigation radius, poor representative of interpretation results from the poor data, low level application of interpretation results. This new technology helps create favorable conditions for early precise reservoir evaluation and reduction of the risk of exploration. 2. The technological difficulties for none flowing well testing are successfully solved by using none flowing formation combined mechanical tool string .This method has been proved by its applications to be able to improve the efficiency of the testing and the quantity of the acquired test data ,and so as to enhance the application of the interpretation results of the test in development of oil fields. 3. The application of the rotary formation tester, selective test valve, well testing string and their allier tools help to resolve problems such as the operation of opening and shutting-in the well under different well conditions, to broaden the scope of well test technology for none flowing formations. 4. Refined Testing Technique for production Wells has greatly shortened the testing dwration and improved the efficiency and accuracy of operation, enriched test results, and at the same time created conditions for conducting multi-well interference well testing.

内蒙古白云鄂博超大型REE-Nb-Fe矿床脉状矿化时代与成矿流体特征

Bayan Obo REE-Nb-Fe ore deposit is the largest REE deposit in the world. Owing to its unique type and tremendous economic value, this deposit has widely attracted interests from geological researchers and vast amount of scientific data have been accumulated. However, its genesis, especially ore-forming age and REE sources, have been under dispute for a long time. On the basis of previous research works, this paper mainly conducts studies on the Early Paleozoic ore-forming event in the Bayan Obo deposit. The following results and conclusions can be suggested: Sm-Nd isotopic analytical results of bastnaesite, beloeilite, albite and fluorite samples from a coarse-crystalline ore lode present an isochron age of 436±35Ma. Besides, Rb-Sr isotope dating of the coarse-crystalline biotite lode that intruded into banded ores gives an isochron age of 459±39Ma. The two ages verify the exist of Early Paleozoic ore-forming event at Bayan Obo, which characterized by extensive netted mineralization of REE fluorocarbonates, aeschynite and monazite, accompanied by widely fluorite-riebeckite-aegirine-apatite alteration. Sr-Nd isotope composition of vein minerals is located between EMI and ancient lower crust component in the ISr(t)-εNd(t) correlation diagram, indicating that there is a crustal contamination during veined mineralization. A large area late Paleozoic granitoids are distributed in the southeast region of east open pit of the mine. The granitoids intruded directly into the ore-bearing dolomite, and produced intense skarnization. Moreover, at 650-660m of the drill core on 22 line and 1598m level flat in the south of East Open Pit, we firstly found skarnization rocks. Single grain and low background Rb-Sr isochrone dating on phlogopite in skarn gives 309±12Ma. Considering the intruded contacting relationship, the late Paleozoic granitoids, already extended to the under part of REE ore bodies, must be posterior to the latest intense REE mineralization, and is only a destructive tectonic and magmatic activity. Fluid inclusion types of fluorite in the Bayan Obo deposit consist of multiphase daughter mineral-bearing inclusion, two or three phase CO2-bearing inclusion and two phase aqueous inclusion. Petrography, laser Raman analysis and microthermometry study indicate that the fluids involving in REE-Nb-Fe mineralization at Bayan Obo might be mainly of H2O-CO2-NaCl-(F-REE) system. The presence of REE-carbonate as a daughter mineral in fluid inclusions shows that the original ore-forming fluids are rich in REE elements.

塔里木盆地塔北、塔中地区早古生代白云岩形成机理

In Tarim Basin, extensive carbonates of Lower Paleozoic occur, in which thick Cambrian and Lower Ordovician dolostones are widespread and show a potential perspective in hydrocarbon exploration. So they are viewed as an important target for exploration. Tarim Basin is a poly phase composite basin, which underwent multiphase tectonic modification and volcanic activities; these exerted significant influences on the basin-fills and basin fluid evolution, thereby the diagenetic history, particularly on the deep-buried Lower Paleozoic dolostones. Referring to the classification of dolomite texture proposed by Gregg & Sibley (1984) and Sibley & Gregg (1987). In view of crystal size, crystal shape, crystal surface and contact relation, eight genetic textures of dolomite crystals are identified, based on careful petrographic examinatoins. These textures include: 1) micritic dolomite; 2) relict mimetic dolomite; 3)finely crystalline, planar-e(s), floating dolomite; 4)finely crystalline, planar-e(s) dolomite; 5) finely-coarse crystalline, nonplanar-a dolomite; 6)coarse crystalline, nonplanar saddle dolomite; 7) finely-medium crystalline, planar-e(s) dolomite cement; 8) coarse crystalline, nonplanar saddle dolomite cement, in which the former six textures occurs as in matrix, the latter two in the cements. Detailed geochemistry analysis is carried out on the basis of genetic textures of dolomite and related minerals such as quartz and calcite. The result showed that the calcite has the highest average content in Sr, which can be sorted into two groups; micritic dolomite has the highest average content in Sr among all kinds of dolomites; the REE patterns of all kinds of dolomites is similar to those of marine limestone samples. Saddle dolomite cement has δ13C values from -2.44‰ to 1.27‰ PDB, and δ18O values from -13.01‰ to -5.12‰ PDB, which partially overlap with those of matrix dolomite (δ13C values from -2.83‰ to 2.01‰ PDB, δ18O values from -10.63‰ to -0.85‰ PDB). Saddle dolomite cement has 87Sr/86Sr ratios from 0.7086 to 0.7104, which totally overlap with those of matrix dolomite (0.7084 ~ 0.7116). Compared with saddle dolomite derived from other basins all over the world, the saddle dolomites of Tarim Basin have similar δ13C, δ18O and 87Sr/86Sr ratios values with those of matrix dolomite. This scenario reflects the unusual geological setting and special dolomitizing liquid of Tarim Basin. The values of δ18O, δ13C and 87Sr/86Sr ratios of calcite also can be sorted out two groups, which may been resulted from the one stage of extensive uplift of Tarim Basin from Mesozoic to Cenozoic. Fluid inclusion microthermometry data of the diagenetic mineral indicates that matrix dolomite has relatively low homogenization temperatures (Th) of 80～105oC and salinities of 12.3％ (wt% NaCl equivalent); saddle dolomite has highest Th values, which concentrate in 120～160oC and salinities of 13.5～23.7％ (wt% NaCl equivalent); quartz has relatively low Th of 135～155oC and salinities of 17.8～22.5％ (wt% NaCl equivalent); calcite has relatively low Th of 121～159.5oC and salinities of 1.4～17.5％ (wt% NaCl equivalent). These data suggest that the saddle dolomites could have formed in thermal brine fluids. Based on comprehensive petrographical study, detailed geochemistry and fluid inclusion microthermometry analysis on Lower Paleozoic dolomite of Tarim Basin, three types of dolomitisation mechanism are proposed: Penecontemporaneous dolomitisation (Sabkha dolomitisation & Reflux dolomitisation); Burial dolomitisation (shallow-intermediate burial dolomitisation & Deep burial dolomitisation ); Hydrothermal cannibalized dolomitisation. In view of host-specified occurrences of hydrothermal dolomite, the low abundance of saddle dolomite and high geochemical similarities between saddle dolomite and host dolomite, as well as highest Th and high salinities , the hydrothermal dolomite in Tarim Basin is thus unique, which could have been precipitated in modified fluid in the host dolomite through intraformational thermal fluid cannibalization of Mg ions from the host. This scenario is different from the cases that large scale dolomitizing fluid migration took place along the fluid pathways where abundant saddle dolomite precipitated. Detailed observations on 180 petrographic and 60 casting thin sections show original pores in Lower Paleozoic dolomite were almost died out by complicated diagenetic process after a long time geologic evolution. On the other hand, deep-buried dolomite reservoirs is formed by tectonic and hydrothermal reforming on initial dolomites. Therefore, the distribution of structure-controlled hydrothermal dolomite reservoirs is predicted in Tabei and Tazhong Area of Tarim Basin based on the geophysical data.

石油二次运移机理物理模拟研究

Oil and gas migration is very important for theoretical hydrocarbon geology study and exploration practice, but related research is weak. Physical simulation is a main method to study oil migration. Systematic experiments were done to quantitatively describe the migration patterns, path characters and oil saturation by adjusting the possible dynamic factors respectively. The following conclusions were drawn. 1. Darcy velocity and pore throat diameter were calculated according to seepage cross-sectional area and glass beads arrangement. With such normalized Darcy velocity and pore throat diameter, the date from one and two dimensional experiments can be reasonably drawn in two phase diagrams. It is found that the migration pattern can be identified using only one dimensionless number L which is defined as the ration of capillary number and Bond number. 2. Oil saturated in the pores between glass beads was used as calibration and oil saturation in the path was measured by magnetic resonance imaging. The results show that oil saturation in the center of migration path can reach 100%, is higher than oil saturation in the edge of migration path. 3. Percolation backbone during secondary oil migration was identified experimentally using Hele-Shaw cell. The backbone formed mainly because of the spatial variation of the cluster conductivity caused by oil saturation heterogeneity, main resistant force change, and path shrinkage and snap-off. Percolation backbone improves hydrocarbon migration efficiency and is a favorable factor for reservoir forming. 4. In the three dimensional filling models, the thickness of the secondary migration path is mall. It is only 2.5cm even for the piston pattern. Inclination of the model is the main influencing factor of the secondary path width.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

三角洲储层孔隙结构模型与剩余油分布研究——以临南油田沙河街组储层为例

The study of pore structure in reservoir was paid attention to in the early reservoir research and now a systematic research methodology is set up. On the limits of tools and conditions, methodologies and technologies on formation condition and distribution laws of pore structure and the relationship between remaining oil distribution and pore structure are uncertain and some knownage about it is also uncertain. As the development of petroleum industry, the characterization of pore structure and the prediction of remaining oil are the hot spot and difficult point in the research of oil development. The author pays a close attention to this subject and has done much research on it. In a case study in Linnan oilfield Huimin sag Jiyang Depression Bohai Bay basin by using a new method, named varied scale comprehensive modeling of pore structure, the author builds pore structure models for delta reservoir, reveals the remaining oil distribution laws in delta facies, and predicts the distribution of remaining oil in Linnan oilfield. By the application of stratigraphy, sedimentology and structure geology. the author reveals the genetic types of sandbody and its distribution laws, builds the reservoir geological models for delta sandstone reservoir in Shahejie group in Linnan oilfield and points out the geological Factors that control the development of pores and throats. Combining petrology and the reservoir sensitive analysis, the author builds the rock matrix models. It is the first time to state that rocks in different sentimental micro facies have different sensitive .response to fluid pressed into the rocks. Normally. the reservoirs in the delta front have weaker sensitivity to fluid than the reservoirs in delta plain, In same subfacies, the microfacies that have fine grain, such as bank and crevasse splay, have stronger reservoir sensitivity than the microfacies that have coarse grains, such as under-water branched channel and debauch bar. By the application of advanced testing, such as imagine analysis, scan electronic microscope, and morphology method, the author classifies the pore structure and set up the distribution models of pore, throat and pore structure. By the application of advanced theory in well-logging geology, the author finds the relationship between microscope pore structure and macroscopic percolation characteristics, and then builds the well-logging interpretation formulae for calculating pore structure parameters. By using the geostatistics methods, the author reveals the spatial correlative characteristics of pore structure. By application of conditional stochastic simulation methods, the author builds the 3D models of pore structure in delta reservoir. It is the base of predicting remaining oil distribution. By a great deal of experiments and theoretical deduction, The author expounds the laws of percolation flow in different pore structures, and the laws by which the pore structure controls the micro distribution of remaining oil, and then, states the micro mechanism of remaining oil distribution. There are two types of remaining oil. They are by-pass flow caused by micro-fingering and truncation caused by non-piston movement. By new method, the author states the different pore structure has different replacement efficiency, reveals the formation condition and distribution laws of remaining oil. predicts the remaining oil distribution in Linnan oil field, and put forward some idea about how to adjust the oil production. The study yielded good results in the production in Linnan oilfield.

温度、压力耦合背景下多相流体运移的数学模型研究

The Mathematical modeling of multiphase fluid flow is an important aspect of basin simulation, and also is a topic of geological frontier. Based on coupling relation of temperature, pressure and fluid flow, this dissertation discusses the modeling which conform to geological regularities of fluid migration. The modeling that is multi-field and multiphase includes heat transport equation, pressure evolvement equation, solution transport equation and fluid transport equation. The finite element method is effective numerical calculation methods. Author applies it to solve modeling and implements the finite element program, and the modeling is applied to Ying-Qiong Basin. The channels of fluid vertical migration are fault, fracture and other high penetrability area. In this thesis, parallel fracture model and columnar channel model have been discussed, and a characteristic time content and a characteristic space content been obtained to illustrate the influences of stratigraphic and hydrodynamic factors on the process. The elliptoid fracture model is established and its approximately solution in theory is gotten. Three kinds of modeling are applied to analyze the transient variation process of fluid pressure in the connected permeable formations. The elliptoid fracture model is the most similar geology model comparing with the other fracture models so the research on this fracture model can enhance the understanding to fluid pressure. In the non-hydrodynamic condition, because of the difference between water density and nature gas density, nature gas can migrate upon by float force. A one-dimension mathematical model of nature gas migration by float force is established and also applied to analyze the change in the saturation of gas. In the process of gas migration its saturation is non-continuous. Fluid flow is an important factor which influences the distribution of the temperature-field, the change of temperature can influence fluid property (including density, viscidity, and solubility),a nd the temperature field has coupling relations to the fluid pressure field. In this dissertation one-dimension and two-dimension thermal convection modeling is developed and also applied to analyze convective and conductive heat transfer. Author has established one-dimension and two-dimension mathematical modeling in which fluid is a mixture of water and nature gas based on the coupling relation between temperature and pressure, discussed mixture fluid convection heat transfer in different gas saturation, and analyzed overpressure form mechanism. Based on geothermal abnormity and pore pressure distribution in Dongfong 1-1, Yinggehai Basin, South China Sea, one-dimension mathematical modeling of coupling temperature and pressure is established. The modeling simulates the process that fluid migrates from deep to shallow and overpressure forms in shallow. When overpressure is so large that fractures appear and overpressure is released. As deep fluid flow to shallow, the high geothermal then forms in shallow. Based on the geological characteristics in Ya13-1, two-dimension mathematical modeling of coupling temperature and pressure is established. Fluid vertically flows in fault and then laterally migrates in reservoir. The modeling simulates the geothermal abnormity and pore pressure distribution in reservoir.

地质体科学计算可视化研究及其应用

Describing visually space-time properties of geological phenomena consists of one of the most important parts in geology research. Such visual images are of usually helpful for analyzing geological phenomena and for discovering the regulations behind geological phenomena. This report studies mainly three application problems of scientific visualization in geology: (Dvisualizing geological body A new geometric modeling technique with trimmed surface patches has been eveloped to visualize geological body. Constructional surfaces are represented as trimmed surfaces and a constructional solid is represented by the upper and lower surface composed of trimmed surface patches from constructional surfaces. The technique can completely and definitely represent the structure of geological body. It has been applied in visualization for the coal deposit in Huolinhe, the aquifer thermal energy storage in Tianjin and the structure of meteorite impact in Cangshan et al. (2)visualizing geological space field Efficient visualization methods have been discussed. Marching-Cube algorithm used has been improved and is used to extract iso~surface from 3D data set, iso-line from 2D data set and iso-point from ID data set. The improved method has been used to visualize distribution and evolution of the abnormal pressures in Zhungaer Basin. (3)visualizing porous space a novel way was proposed to define distance from any point to a convex set. Thus a convex set skeleton-based implicit surface modeling technique is developed and used to construct a simplified porous space model. A Buoyancy Percolation numerical simulation platform has been developed to simulate the process of migration of oil in the porous media saturated with water.

石油二次运移的物理与数值模拟实验及其应用

Facing the problems met in studies on predominant hydrocarbon migration pathways, experiments and numerical simulating were done in this thesis work to discuss the migration mechanisms. The aim is to analyze quantitatively the pathway pattern in basin scale and to estimate the hydrocarbon loss on the pathway that offer useful information for confirming the potential hydrocarbon accumulation. Based on our understandings on hydrocarbon migration and the fluid dynamic theory, a series of migration experiments were designed to observe the phenomena where kerosene is used as draining phase driven only by buoyancy force that expulses pore water. These experiments allow to study the formation of migration pathways, the distribution of non-wetting oil along these pathways, and the re-utilizing of previously existing pathways marked by residual traces etc. The types of pattern for migration pathways may be characterized by a phase diagram using two dimensionless numbers: the capillary number and the Bond number. The NMR technique is used to measure the average saturation of residual oil within the pathways. Based our experiment works and percolation concept, a numerical simulation model were proposed and realized. This model is therefore called as BP (Buoyancy Percolation) simulator, since buoyancy is taken as the main driving force in hydrocarbon migration. To make sure that BP model is applicable to simulate the process of oil secondary migration, the experimental phenomena are compared with those simulated with BP model by fractal method, and the result is positive. After then, we use BP simulator to simulate the process of migration of oil in the porous media saturated with water at different scale. And the results seem similar to those cited in literatures. In addition, our software is applied in Paris basin to predict the pathway of hydrocarbon migration happened in the Middle Jurassic reservoirs. It is found that the results obtained with our BP model are generally agree with Hindle (1997) and Bekeles'(1999), but our simulated migration pathway pattern and migration direction seem more reasonable than theirs.

大别（超）包产到户铁岩和阜新火山岩铂族元素斌存规律研究

The platinum-group elements (PGE), including Os, Ir, Ru, Rh, Pt and Pd, axe strongly siderophile and chalcophile. On the basis of melting temperature, the PGE may be divided into two groups: the Ir group (IPGE, >2000°C) consisting of Os, Ir and Ru, and the Pd group (PPGE, <20GO°C) consisting of Rh, Pt and Pd. Because of their unique geochemical properties, PGE provide critical information on global-scale differentiation processes, such as core-mantle segregation, late accretionary history, and core-mantle exchange. In addition, they may be used to identify magma source regions and unravel complex petrogenetic processes including partial melting, melt percolation and metasomatism in the mantle, magma mixing and crustal contamination in magma chambers and melt crystallization.Compared with other rocks, (ultra)mafic rocks have lower REE content but higher PGE content, so PGE have advantages in studying the petrogeneses and evolution of them. In this study, we selected (ultra)mafic rocks collected in Dabie Orogen and volcanic rocks from Fuxin Region. Based on the distribution and behaviour of platinum-group elements, combined with other elements, we speculate the magma evolution and source mantle of these (ultra)mafic rocks and volcanic rocks.Many (ultra)mafic rocks are widely distributed in Dabie Region. According to their deformation and metamorphism, we classed them into three types. One is intrusive (ultra)mafic rocks, which are generally undeformed and show no or little sign of metamorphism, such as (ultra)mafic intrusions in Shacun, zhujiapu, Banzhufan, qingshan, Xiaohekou, Jiaoziyan, Renjiawan and Daoshichong. The other one is ultrahigh pressure metamorphic (ultra)mafic rocks, some of them appeared as eelogites, such as complex in Bixiling and adjacent Maowu. Another one is intense deformed and metamorphic, termed as tectonic slice, alpine-type (ultra)mafic rocks. The most representative is Raobazhai and Dahuapin. However, there are many controversies about the formation of those (ultra)mafic rocks. Here, we select typical rocks of the three types. The PGE were determined by inductively coupled plasma mass spectrometry (ICP-MS) ater NiS fire-assay and tellurium co-precipitation.The PGE tracing shows that three components are needed in the source of the cretaceous (uitra)mafic intrusions. They could be old enriched sub-continental lithospheric mantle, lower crust and depleted asthenospheric mantle. The pattern of PGE also shows the primitive magma of these intrusions underwent S saturation. According to palladium, we can conclude that the mantle enrich in PGE. Distribution of PGE in Bixiiing and Maowu (ultra)mafic rocks display they are products of magmas fractional crystallization. The (ultra)mafic rocks in Bixiiing and Maowu are controlled by various magmatic processes and the source mantle is depleted in PGE. Of interest is that the mantle produced UHP (ultra)mafic rocks are PGE-depleted, whereas the mantle of cretaceous (ultra)mafic intrusions are enrich in PGE. This couldindicate that the mantle change from PGE-enriched to PGE-depleted during120-OOMa, which in accord with the time of tectonic system change in the East China. At the same time, (ultra)mafic intrusions in cretaceous took information of deep mantle, which means the processes in deep mantle arose structural movement in the crust The character of PGE in alpine-type (ultra)mafic rocks declared that the rocks had experienced two types of metasomatic processes - hydrous melt derived from slab and silicate melt. In addition, we analyze the platinum-group elements in volcanic rocks on the northern margin of the North China Craton, Fuxin. The volcanic rocks characterized by negative anomalies of platinum. This indicates that platinum alloys, which may host some Pt resided in the mantle. The PGE patterns also show that Jianguo alkali basalts derived from asthenospheric mantle source, but wulahada high-Mg andesites derived from lithospheric mantle.

Biomass?Diversity Responses and Spatial Dependencies in Disturbed Tallgrass Prairies

Javier G. P. Gamarra and Ricard V. Sole (2002). Biomass-diversity responses and spatial dependencies in disturbed tallgrass prairies. Journal of Theoretical Biology, 215 (4) pp.469-480 RAE2008

Safe Compositional Specification of Networking Systems

The Science of Network Service Composition has clearly emerged as one of the grand themes driving many of our research questions in the networking field today [NeXtworking 2003]. This driving force stems from the rise of sophisticated applications and new networking paradigms. By "service composition" we mean that the performance and correctness properties local to the various constituent components of a service can be readily composed into global (end-to-end) properties without re-analyzing any of the constituent components in isolation, or as part of the whole composite service. The set of laws that would govern such composition is what will constitute that new science of composition. The combined heterogeneity and dynamic open nature of network systems makes composition quite challenging, and thus programming network services has been largely inaccessible to the average user. We identify (and outline) a research agenda in which we aim to develop a specification language that is expressive enough to describe different components of a network service, and that will include type hierarchies inspired by type systems in general programming languages that enable the safe composition of software components. We envision this new science of composition to be built upon several theories (e.g., control theory, game theory, network calculus, percolation theory, economics, queuing theory). In essence, different theories may provide different languages by which certain properties of system components can be expressed and composed into larger systems. We then seek to lift these lower-level specifications to a higher level by abstracting away details that are irrelevant for safe composition at the higher level, thus making theories scalable and useful to the average user. In this paper we focus on services built upon an overlay management architecture, and we use control theory and QoS theory as example theories from which we lift up compositional specifications.

JOSTLE: multilevel graph partitioning software: an overview

In this chapter we look at JOSTLE, the multilevel graph-partitioning software package, and highlight some of the key research issues that it addresses. We first outline the core algorithms and place it in the context of the multilevel refinement paradigm. We then look at issues relating to its use as a tool for parallel processing and, in particular, partitioning in parallel. Since its first release in 1995, JOSTLE has been used for many mesh-based parallel scientific computing applications and so we also outline some enhancements such as multiphase mesh-partitioning, heterogeneous mapping and partitioning to optimise subdomain shape