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

Short columns with molecularly imprinted monolithic stationary phases for rapid separation of diastereomers and enantiomers

Three molecularly imprinted monolithic columns with different length but almost identical column volume had been prepared. It was observed that the separation factors of diastereomers and enantiomers were almost unaffected by column length. However, the short column with dimension of 38 mm x 8 mm W. showed much lower resistance to flow rate so that it could be operated at much higher flow rates. By combining stepwise gradient elution with elevated flow rate, the diastereomers of cinchonine and cinchonidine and the enantiomers of Cbz-DL-Trp and Fmoc-DL-Trp were successfully separated within 3 min on the short column with dimension of 38 mm. x 8 mm i.d.. Based on the above results, a cinchonine imprinted monolithic disk with dimension of 10 mm x 16 mm W. was further developed. The SEM image and the pore size distribution profile showed that large flow-through pores are present on the prepared monolith, which allowed mobile phase to flow through the disk with very low resistance. Chromatographic performances on the monolithic disk were almost unchanged compared with the long columns. A rapid separation of cinchonine and cinchonidine was achieved in 2.5 min at the flow rate of 9.0 ml/min. Furthermore, it was observed that there was almost no effect of the flow rate on the dynamic binding capacity at high flow rates. In addition, the effect of the loading concentration of analytes on the dynamic binding capacity, namely adsorption isotherm, was also investigated. A non-linear adsorption isotherm of cinchonine was observed on the molecularly imprinted monolith with cinchonine as template, which might be a main reason to result in the peak tailing of template molecule. (C) 2004 Elsevier B.V. All rights reserved.

A simple preparation of stationary phase of alkylamide reversed-phase liquid chromatograph for the separation of basic substances

A simple preparation process of alkylamide phase for reversed-phase HPLC (RP-HPLC) is described. The process includes aminopropyltrimethoxysilane firstly reacted with octanoyl chloride, then the intermediate was coupled onto porous silica. The resultant bonded silica has a reproducible ligand surface concentration and homogenous bonded ligand distribution on the porous silica. Characterization of prepared packing was carried out with elemental analysis, solid-state C-13 NMR and Fourier transform infrared (FT-IR). Chromatographic evaluations were carried out by using a mixture of organic compounds including acidic, basic and neutral analytes under methanol/water as binary mobile phase. The results showed that the stationary phase have excellent chromatographic properties and can be efficiently used for the separation of basic compounds.

Intensity Distribution Design for Laser-Induced Thermal Loading Based on Numerical Simulation

To accomplish laser-induced thermal loading simulation tests for pistons,the Gaussian beam was modulated into multi-circular beam with specific intensity distribution.A reverse method was proposed to design the intensity distribution for the laser-induced thermal loading based on finite element(FE) analysis.Firstly,the FE model is improved by alternating parameters of boundary conditions and thermal-physical properties of piston material in a reasonable range,therefore it can simulate the experimental resul...

Effects of the orientation of microgravity on the concentration distribution in crystallization from solution

A quasi-steady state growth and dissolution in a 2-D rectangular enclosure is numerically investigated. This paper is an extension to indicate the effects of the orientation of gravity on the concentration field in crystallization from solution under microgravity, especially on the lateral non-uniformity of concentration distribution at the growth surface. The thermal and solute convection are included in this model.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Concentration distribution in solution crystal growth: effect of moving interface conditions

The linear diffusion-reaction theory with finite interface kinetics is employed to describe the dissolution and the growth processes. The results show that it is imperative to consider the effect of the moving interfaces on the concentration distribution at the growth interface for some cases. For small aspect ratio and small gravity magnitude, the dissolution and the growth interfaces must be treated as the moving boundaries within an angle range of 0 degrees < gamma < 50 degrees in this work. For large aspect ratio or large gravity magnitude, the effect of the moving interfaces on the concentration distribution at the growth interface can be neglected except for gamma < - 50 degrees.

Water Content Distribution in the Surface Layer of Maoping Slope

The water content distribution in the surface layer of Maoping slope has been studied by testing the water content at 31 control sites. The water content profiles at these sites have also been determined. The water content distributions at different segments have been obtained by using the Kriging method of geostatistics. By comparing the water content distributions with the landform of the slope, it was shown that the water content is closely dependent on the landform of the slope. The water content distribution in the surface layer provided a fundamental basis for landslide predication and treatment.

Thermodynamical versus log-Poisson distribution in turbulence

The thermodynamical model of intermittency in fully developed turbulence due to Castaing (B. Castaing, J. Phys. II France 6 (1996) 105) is investigated and compared with the log-Poisson model (Z-S, She, E. Leveque, Phys. Rev. Lett. 72 (1994) 336). It is shown that the thermodynamical model obeys general scaling laws and corresponds to the degenerate class of scale-invariant statistics. We also find that its structure function shapes have physical behaviors similar to the log-Poisson's one. The only difference between them lies in the convergence of the log-Poisson's structure functions and divergence of the thermodynamical one. As far as the comparison with experiments on intermittency is concerned, they are indifferent.

Effects on Dense-phase Zone Temperature of Fluidized Bed with Uneven Air Distribution

由于采用非均匀布风,内旋流流化床的移动区空气量不足,导致燃烧不充分,温度较低。当移动区未流化时,密相区内存在较明显的温度不均匀性。随着移动区流速的提高,温度差迅速减小。当移动区流速超过2.0#mu#m后,密相区温度基本均匀一致。流动区流速对密相区温度均匀有一定的影响,流速越高,温度越均匀。

In situ formation by laser cladding of a TiC composite coating with a gradient distribution

An in situ method was developed to produce an Ni alloy composite coating reinforced by in situ reacted TiC particles with a gradient distribution, using one-step laser cladding with a pre-placed powder mixture on a 5CrMnMo steel substrate. Dispersed and ultra-fine TIC particles were formed in situ in the coating. Most. of the TiC particles, with a marked gradient distribution, were uniformly distributed within interdendritic regions because of the trapping effect of the advancing solid-liquid interface. In addition, the TiC-gamma-Ni interfaces generated in situ were found to be free from any deleterious surface reaction. Finally, the microhardness also showed a gradient variation, with the highest value of 1250 Hv0.2 and the wear properties of the coating were significantly enhanced.

Experimental Investigation On The Holdup Distribution Of Oil-Water Two-Phase Flow In Horizontal Parallel Tubes

An experimental investigation was conducted to study the holdup distribution of oil and water two-phase flow in two parallel tubes with unequal tube diameter. Tests were performed using white oil (of viscosity 52 mPa s and density 860 kg/m(3)) and tap water as liquid phases at room temperature and atmospheric outlet pressure. Measurements were taken of water flow rates from 0.5 to 12.5 m(3)/h and input oil volume fractions from 3 to 94 %. Results showed that there were different flow pattern maps between the run and bypass tubes when oil-water two-phase flow is found in the parallel tubes. At low input fluid flow rates, a large deviation could be found on the average oil holdup between the bypass and the run tubes. However, with increased input oil fraction at constant water flow rate, the holdup at the bypass tube became close to that at the run tube. Furthermore, experimental data showed that there was no significant variation in flow pattern and holdup between the run and main tubes. In order to calculate the holdup in the form of segregated flow, the drift flux model has been used here.

Reconstructing The Vertical Distribution Of The Aeolian Saltation Mass Flux Based On The Probability Distribution Of Lift-Off Velocity

The probability distribution of lift-off velocity of the saltating grains is a bridge to linking microscopic and macroscopic research of aeolian sand transport. The lift-off parameters of saltating grains (i.e., the horizontal and vertical lift-off velocities, resultant lift-off velocity, and lift-off angle) in a wind tunnel are measured by using a Phase Doppler Particle Analyzer (PDPA). The experimental results show that the probability distribution of horizontal lift-off velocity of saltating particles on a bed surface is a normal function, and that of vertical lift-off velocity is an exponential function. The probability distribution of resultant lift-off velocity of saltating grains can be expressed as a log-normal function, and that of lift-off angle complies with an exponential function. A numerical model for the vertical distribution of aeolian mass flux based on the probability distribution of lift-off velocity is established. The simulation gives a sand mass flux distribution which is consistent with the field data of Namikas (Namikas, S.L., 2003. Field measurement and numerical modelling of acolian mass flux distributions on a sandy beach, Sedimentology 50, 303-326). Therefore, these findings are helpful to further understand the probability characteristics of lift-off grains in aeolian sand transport. (c) 2007 Elsevier B.V. All rights reserved.

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.