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

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.

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.

胶东地区三类金矿成矿系统地质地球化学及动力学

The East Shandong gold province is located on the southeastern margin of the North China Craton and features uplift in the north and depression in the south. The uplift area is made up of the Archaean Jiaodong Group, the Proterozoic Jingshan Group and Yanshannian granites. Most gold deposits in the uplift area are spatially associated with the Yanshannian granites. Two types of gold mineralization occur in the region: the quartz-vein type hosted in the Linglong granite suite, and the shear zone type hosted by either the Linglong granite or Guojialing granitoid suites. The mineralization ages are 113～126 Ma. The southern part of East Shandong contains the Mesozoic Jiaolai basin, which formed during regional extension. The basin is bounded by the Wulian-Rongcheng fault in the southeast and the Tanlu fault in the west. The Pengjiakuang, Fayunkuang and Dazhuangzi gold deposit occurs on the northeastern margin of the basin. The mineralization ages of these deposits are 110～128 Ma. This paper focuses on a low-angle detachment fault developed between the Proterozoic Jingshan Group metamorphic complex and the northeastern margin of the basin. Our field work shows that the distribution of the Pengjiakuang gold deposit was controlled by the detachment fault. Moreover, the Fayunkuang, Guocheng and Liaoshang gold deposits also occurr in the periphery of the basin, and their features are similar to Pengjiakuang gold deposit. The study of geological geochemistry of the gold deposits has shown: ①three-type gold deposit was situated in the Jiaodong area, including altered rock type (Jiaojia type), quartz vein type (Linglong type) and breccia type (Pengjiakuang type); the ore-forming materials and fluid for Pengjiakuang type gold deposit shows multiple source; ②the ore materials of Jiaojia and Linglong type deposits are mainly from deep source. The author has studied geological-geochemical dynamics of three types deposits in Jiaodong area. The study of tectonic dynamics shows that ore-forming structure differential stress values of Pengjiakuang gold deposit is 100 * 10~6～130 * 10~6 Pa, and that of Jiaojia gold deposit is 100 * 10~5～194 * 10~6 Pa. Dynamics of hydrothermal ore-forming fluid has also been studied in this paper. Author applies Bernoulli equation to dynamic model of hydrothermal fluid motion in brittle fracture and cracks (quartz vein type gold mineralization), and applies Darcy law to dynamic model of hydro thermal fluid motion in porous medium (altered rock type gold mineralization). Author does daring try in order to study quantitativly transport mechanism of hydrothermal ore-forming fluid in this paper. The study of fluid inclusions and crystal dynamics shows that reaction system of hydrothermal ore-forming includes three types, as follows: ore-forming reaction, controlling reaction and buffer controlling reaction. They depend on each other, controlling each other, which form a organic system. Further research shown that formation of ore shoots was controlled by coincidence processes of tectonic dynamic condition and thermodynamic evolution. This paper has summaried reginoal metallogenic laws and seted up metallogenic(dynamics) models for Jiaodong gold ore belt.

Thermal effects on Rayleigh-Marangoni-Benard instability in a system of superposed fluid and porous layers

Thermal effects of the heat transfer at free surface (represented by Biot number) on the Rayleigh-Marangoni-Benard instability in a system of liquid-porous layers with top free surface are investigated numerically. The results indicate that this thermal effect can evidently lead to the mode transition of convection, which is overlooked in previous works. (C) 2010 Elsevier Ltd. All rights reserved.

Characters of Surface Deformation and Surface Wave in Thermal Capillary Convection

In the field of fluid mechanics, free surface phenomena is one of the most important physical processes. In the present research work, the surface deformation and surface wave caused by temperature difference of sidewalls in a rectangular cavity have been investigated. The horizontal cross-section of the container is 52 mmx42 mm, and there is a silicon oil layer of height 3.5 mm in the experimental cavity. Temperature difference between the two side walls of the cavity is increased gradually, and the flow on the liquid layer will develop from stable convection to un-stable convection. An optical diagnostic system consisting of a modified Michelson interferometer and image processor has been developed for study of the surface deformation and surface wave of thermal capillary convection. The Fourier transformation method is used to interferometer fringe analysis. The quantitative results of surface deformation and surface wave have been calculated from a serial of the interference fringe patterns.The characters of surface deformation and surface wave have been obtained. They are related with temperature gradient and surface tension. Surface deformation is fluctuant with time, which shows the character of surface wave. The cycle period of the wave is 4.8 s, and the amplitudes are from 0 to 0.55 mu m. The phase of the wave near the cool side of the cavity is opposite and correlative to that near the hot side. The present experiment proves that the surface wave of thermal capillary convection exists on liquid free surface, and it is wrapped in surface deformation.

Thermocapillary Effect in a Thermal Rheological Jet

The process of die swell in polymer jets is an important feature within polymer processing and can be explained through a study of its rheological effects. The existence of a thermocapillary effect, driven by the gradient of its surface tension, should be considered when examining a thermal jet that has a non-uniform temperature distribution on its free surface, as in various polymer processings. Both the rheological effect and thermocapillary effect on die swell can be studied numerically through a finite element method as used on a two-dimensional and unsteady model, in which a Coleman-Noll second-order fluid model is employed. The results show that the expanding angle depends on both the rheological property of the fluid and the pressure at the vessel exit. Although both the thermocapillary and the rheological effects contribute to the cross-section expansion of the fluid jet, the latter is more important in determining the expansion.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Experimental Study on Surface Deformation and Surface Wave of Fluid Convection with Free Surface

An optical diagnostic system consisting of Michelson interferometer with image processor has been developed for study of the kinetics of thermal capillary convection and buoyancy convection. This optical interferometer has been used to observe and measure surface deformation and surface wave of capillary convection and buoyancy convection in a rectangular cavity with different temperature’s sidewalls. Fourier transformation is used to image processing. The quantitative results of surface deformation and surface wave have been calculated from the interference fringe pattern. With the increasing of temperature gradient, the liquid surface slant gradually. It’s deformation has been calculated, which is related directly with temperature gradient. This is one of the characters introducing convection. Another interesting phenomenon is the inclining direction, which is different when the liquid layer is thin or thick. When the liquid layer is thin, convection is mainly controlled by thermocapillary effect. However, When the liquid layer is thick, convection is mainly controlled by buoyancy effect. Surface deformation in the present experiment are more and more declining in this process. The present experiment proved that surface deformation appears before the appearance of surface wave on fluid convection, it is related with temperature gradient, and the height of liquid layer, and lies on capillary convection and buoyancy convection. The present experiment also demonstrates that the amplitude of surface wave of thermocapillary-buoyancy convection is much smaller than surface deformation, the wave is covered by deformation.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Monte Carlo Simulation Of Thermal Fluctuations Below The Onset Of Rayleigh-Benard Convection

The density fluctuations below the onset of convection in the Rayleigh-Benard problem are studied with the direct simulation Monte Carlo method. The particle simulation results clearly show the connection between the static correlation functions of fluctuations below the critical Rayleigh number and the flow patterns above the onset of convection for small Knudsen number flows (Kn=0.01 and Kn=0.005). Furthermore, the physical nature for no convection in the Rayleigh-Benard problem under large Knudsen number conditions (Kn>0.028) is explained based on the dynamics of fluctuations.