三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

The toxic effects of microcystin-LR on embryo-larval and juvenile development of loach, Misguruns mizolepis Gunthe

Microcystin-LR, a specific and potent hepatotoxin, was tested for its effects oil loach embryo-larval and juvenile development, The results of this study showed that loach embryos were more sensitive when exposed to microcystin-LR at a later than at an earlier stage of development, Juveniles were far less sensitive to MC-LR than were embryos and larvae. Mortality and developmental abnormality were proven to be dose-dependent and to be stage-specific sensitive. Among the abnormal changes noted were: pericardial edema and tubular heart, bradycardia, homeostasis, poor yolk resumption. small head, curved body and tail, and abnormal hatching, Liver and heart were the main targets of microcystin-LR toxicity. Ultrastructural analysis documented a complex set of sublethal effects of microcystin-LR on loach hepatocytes, chiefly including morphological alteration in nuclear and RER of loach liver cells. fit addition, microcystin-LR was lethal to loach juvenile in the subacute (7 days) exposure (LC50) = 593.3 mug/l). (C) 2002 Elsevier Science Ltd. All rights reserved.

Thermally induced evolution of phase transformations in gas hydrate sediment

Thermally induced evolution of phase transformations is a basic physical-chemical process in the dissociation of gas hydrate in sediment (GHS). Heat transfer leads to the weakening of the bed soil and the simultaneous establishment of a time varying stress field accompanied by seepage of fluids and deformation of the soil. As a consequence, ground failure could occur causing engineering damage or/and environmental disaster. This paper presents a simplified analysis of the thermal process by assuming that thermal conduction can be decoupled from the flow and deformation process. It is further assumed that phase transformations take place instantaneously. Analytical and numerical results are given for several examples of simplified geometry. Experiments using Tetra-hydro-furan hydrate sediments were carried out in our laboratory to check the theory. By comparison, the theoretical, numerical and experimental results on the evolution of dissociation fronts and temperature in the sediment are found to be in good agreement.

Comparison of comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry and gas chromatography-mass spectrometry for the analysis of tobacco essential oils

A sample of tobacco essential oil was analyzed using gas chromatography-mass spectrometry (GUMS)and comprehensive two-dimensional gas chromatography coupled to a time-of-flight mass spectrometry (GC x GC/TOFMS), respectively. In the GUMS analysis, serially coupled columns were used. By comparing the GUMS results with GC x GC/TOFMS result,,, many more components in the essential oil could be found within the two-dimensional separation space of GC x GC. The quantitative determination of components in the essential oil was performed by GC x GC with flame ionization detection (FID), using a method of multiple internal standards calibration, (c) 2005 Elsevier B.V. All rights reserved.

A Research on New Method of Transportation Technics for Heavy Oil and Practice in Oil Field

介绍一种可应用于高粘度稠油管输的新工艺。即用自行研制的蒸汽引射器采用无界引射方式,将蒸汽直接注入到输油管道中,利用蒸汽释放的热量提高稠油温度降低粘度,从而达到降低稠油输送压降的目的,它比间接加热输送工艺所用的蒸汽量或耗煤量大大减少。方法在辽河油田输油管线上进行了工业现场试验,取得了很好的效果。

Experimental Investigation On The Slip Between Oil And Water In Horizontal Pipes

This work is devoted to study of the slip phenomenon between phases in water-oil two-phase flow in horizontal pipes. The emphasis is placed on the effects of input fluids flow rates, pipe diameter and viscosities of oil phase on the slip. Experiments were conducted to measure the holdup in two horizontal pipes with 0.05 m diameter and 0.025 m diameter, respectively, using two different viscosities of white oil and tap water as liquid phases. Results showed that the ratios of in situ oil to water velocity at the pipe of small diameter are higher than those at the pipe of big diameter when having same input flow rates. At low input water flow rate, there is a large deviation on the holdup between two flow systems with different oil viscosities and the deviation becomes gradually smaller with further increased input water flow rate. (C) 2008 Elsevier Inc. All rights reserved.

Preparation, characterization and gas sensing properties of high soluble metal (II) phthalocyanine thin films by spin-coating method

Three novel metal (II) phthalocyanine complexes were synthesized by cyclic tetramerisation reaction of a dicyano benzene component and different metal ions (Pd2+, Co2+, Zn2+). The structure of complexes was confirmed by elemental analysis, mass and IR spectrum. The excellent solubility of the complexes in benzene enabled us to obtain films by a spin-coating method. The films were characterized by IR, electronic spectral and AFM. The gas sensing properties to NO2 of the metal (II) phthalocyanine complex films were studied. In addition, the effects of different metal ions and the gas sensing temperature on the sensing properties were studied. (C) 2005 Elsevier B.V. All rights reserved.

Preparation, characterization and gas sensing properties of lead tetra-(tert-butyl)-5,10,1 5,20-tetraazaporphyrin spin-coating films

A spin-coated film of lead tetra-(tert-butyl)-5,10,15,20-tetraazaporphyrin complex (PbTAP(t-Bu)(4)) was obtained and characterized by IR spectra, absorption spectra and atomic force microscopy. The response and recovery characteristics of the film to NH3, NO2 and C2H5OH vapor were investigated at room temperature. In addition, the reversibility and stability of the film to NH3 were also studied. The results indicate that the PbTAP(t-Bu)(4) derivative can be exploited as an NH3 sensor at room temperature. (c) 2007 Elsevier B.V. All rights reserved.

Preparation, characterization and gas sensing properties of nickel octa-iso-pentyloxynaphthalocyanine spin-coating films

Spin-coated films of nickel 1,6,10,15,19,24,28,33-octa-iso-pentyloxy-2,3-naphthalocyanine complex were obtained and characterized by UV-vis absorption spectroscopy. A linear relationship between the absorbance and solution concentration was observed. Low concentration solutions could afford smooth and homogeneous film surfaces as indicated by atomic force microscopy. The film structure was studied by small angle X-ray diffraction. The films were used for NO2 sensing experiments. The results indicate that the elevation of sensing temperature can shorten the response time and increase recovery ratio and response magnitude of the sensing films. High NO2 concentration can also shorten response time. (C) 2008 Elsevier B.V. All rights reserved.