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

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.

内蒙古白云鄂博超大型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.

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

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 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.