整合多尺度地球物理数据建立高分辨率储层模型的方法及应用

Most of the fields in China are in the middle-late development phase or are mature fields. It becomes more and more difficult to develop the remaining oil/gas. Therefore, it is import to enhance oil/gas recovery in order to maintain the production. Fine scale modeling is a key to improve the recovery. Incorporation of geological, seismic and well log data to 3D earth modeling is essential to build such models. In Ken71 field, well log, cross-well seismic and 3D seismic data are available. A key issue is to build 3D earth model with these multi-scales data for oil field development.In this dissertation, studies on sequential Gaussian-Bayesian simulation have been conducted. Its comparison with cokriging and sequential Gaussian simulation has been performed. The realizations generated by sequential Gaussian-Bayesian simulation have higher vertical resolution than those generated by other methods. Less differences between these realization and true case are observed. With field data, it is proved that incorporating well log, cross-well seismic and 3D seismic into 3D fine scale model is reliable. In addition, the advantages of sequential Gaussian-Bayesian simulation and conditions for input data are demonstrated. In Ken71 field, the impedance difference between sandstone and shale is small. It would be difficult to identify sandstone in the reservoir with traditional impedance inversion. After comparisons of different inversion techniques, stochastic hillclimbing inversion was applied. With this method, shale content inversion is performed using 3D seismic data. Then, the inverted results of shale content and well log data are incorporated into 3D models. This demonstrates a procedure to build fine scale models using multi scale seismic data, especially 3D seismic amplitude volume.The models generated through sequential Gaussian-Bayesian simulation have several advantages including: (1) higher vertical resolution compared with 3D inverted acoustic impedance (AI); (2) consistency of lateral variation as 3D inverted AI; (3) more reliability due to integration cross-well seismic data. It is observed that the precision of the model depends on the 3D inversion.

Application of rough set-based analysis to extract spatial relationship indicator rules: An example of land use in Pearl River Delta

Spatial relations, reflecting the complex association between geographical phenomena and environments, are very important in the solution of geographical issues. Different spatial relations can be expressed by indicators which are useful for the analysis of geographical issues. Urbanization, an important geographical issue, is considered in this paper. The spatial relationship indicators concerning urbanization are expressed with a decision table. Thereafter, the spatial relationship indicator rules are extracted based on the application of rough set theory. The extraction process of spatial relationship indicator rules is illustrated with data from the urban and rural areas of Shenzhen and Hong Kong, located in the Pearl River Delta. Land use vector data of 1995 and 2000 are used. The extracted spatial relationship indicator rules of 1995 are used to identify the urban and rural areas in Zhongshan, Zhuhai and Macao. The identification accuracy is approximately 96.3%. Similar procedures are used to extract the spatial relationship indicator rules of 2000 for the urban and rural areas in Zhongshan, Zhuhai and Macao. An identification accuracy of about 83.6% is obtained.

Energy conversion in shape memory alloy heat engine - Part II: Simulation

In this paper, a theoretical model proposed in Part I (Zhu et al., 2001a) is used to simulate the behavior of a twin crank NiTi SMA spring based heat engine, which has been experimentally studied by Iwanaga et al. (1988). The simulation results are compared favorably with the measurements. It is found that (1) output torque and heat efficiency decrease as rotation speed increase; (2) both output torque and output power increase with the increase of hot water temperature; (3) at high rotation speed, higher water temperature improves the heat efficiency. On the contrary, at low rotation speed, lower water temperature is more efficient; (4) the effects of initial spring length may not be monotonic as reported. According to the simulation, output torque, output power and heat efficiency increase with the decrease of spring length only in the low rotation speed case. At high rotation speed, the result might be on the contrary.

Molecular dynamics simulation of the role of dislocations in microcrack healing

The molecular dynamics method is used to simulate microcrack healing during heating or/and under compressive stress. A centre microcrack in Cu crystal would be sealed under compressive stress or by heating. The role of compressive stress and heating in crack healing was additive. During microcrack healing, dislocation generation and motion occurred. When there were pre-existing dislocations around the microcrack, the critical temperature or compressive stress necessary for microcrack healing would decrease, and, the higher the number of dislocations, the lower the critical temperature or compressive stress. The critical temperature necessary for microcrack healing depended upon the orientation of the crack plane. For example, the critical temperature for the crack along the (001) plane was the lowest, i.e. 770K.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Structural Failure Analysis and Numerical Simulation of Micro-Accelerometers under Impulsive Loading

Micromachined accelerometer is a kind of inertial MEMS devices, which usually operate under intensive impact loading. The reliability of micromachined accelerometers is one of the most important performance indices for their design, manufacture and commer

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.

Temperature measurement of reflected shock wave by using chemical indicator

This report describes a new method for measuring the temperature of the gas behind the reflected shock wave in shock tube, corresponding to the reservoir temperature of a shock tunnel, based on the chemical reaction of small amount of CF4 premixed in the test gas. The final product C2F4 is used as the temperature indicator, which is sampled and detected by a gas chromatography in the experiment. The detected concentration of C2F4 is correlated to the temperature of the reflected shock wave with the initial pressure P-1 and test time tau as parameters in the temperature range 3 300 K < T < 5 600 K, pressure range 5 kPa < P1 <12 kPa and tau similar or equal to 0.4 ms.

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Two-dimensional simulation of an electron cyclotron resonance plasma source with self-consistent power deposition

Using a refined two-dimensional hybrid-model with self-consistent microwave absorption, we have investigated the change of plasma parameters such as plasma density and ionization rate with the operating conditions. The dependence of the ion current density and ion energy and angle distribution function at the substrate surface vs. the radial position, pressure and microwave power were discussed. Results of our simulation can be compared qualitatively with many experimental measurements.

Formation of Fivefold Deformation Twins in Nanocrystalline Face-Centered-Cubic Copper Based on Molecular Dynamics Simulations

Fivefold deformation twins were reported recently to be observed in the experiment of the nanocrystalline face-centered-cubic metals and alloys. However, they were not predicted previously based on the molecular dynamics (MD) simulations and the reason was thought to be a uniaxial tension considered in the simulations. In the present investigation, through introducing pretwins in grain regions, using the MD simulations, the authors predict out the fivefold deformation twins in the grain regions of the nanocrystal grain cell, which undergoes a uniaxial tension. It is shown in their simulation results that series of Shockley partial dislocations emitted from grain boundaries provide sequential twining mechanism, which results in fivefold deformation twins. (c) 2006 American Institute of Physics.

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