温度、压力耦合背景下多相流体运移的数学模型研究

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.

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

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.

Self-consistent kinetic modeling of electron distribution function in a low-pressure inductively coupled plasma

Based upon the spatially inhomogeneous Boltzmann equation in two-term approximation coupled with electromagnetic and fluid model analysis for the recently developed inductively coupled plasma sources, a self-consistent electron kinetic model is developed. The electron distribution function, spatial distributions of the electron density and ionization rate are calculated and discussed.

Pressure behaviour of self-assembled In0.55Al0.45As/Al0.5Ga0.5As quantum dots with multi-modal distribution in size

The pressure behaviour of In0.55Al0.45As/Al0.5Ga0.5As self-assembled quantum dots (QDs) has been studied at 15 K in the pressure range of 0-1.3 GPa. The atomic force microscopy image shows that the QDs have a multi-modal distribution in size. Three emission peaks were observed in the photoluminescence (PL) spectra, corresponding to the different QD families. The measured pressure coefficients are 82, 93 and 98 meV GPa(-1) for QDs with average lateral size of 26, 52 and 62 nm, respectively. The pressure coefficient of small QDs is about 17% smaller than that of bulk In0.55Al0.45As An envelope-function calculation was used to analyse the effect of pressure-induced change of barrier height, effective mass and dot size on the pressure coefficients of QDs. The Gamma-X state mixing was also included in the evaluation of the reduction of the pressure coefficients. The results indicate that both the pressure-induced increase of effective mass and Gamma-X mixing respond to the decrease of pressure coefficients, and the Gamma-X mixing is more important for small dots. The calculated Gamma-X interaction potentials are 15 and 10 meV for QDs with lateral size of 26 and 52 nm, respectively. A type-II alignment for the X conduction band is suggested according to the pressure dependence of the PL intensities. The valence-band offset was then estimated as 0.15 +/- 0.02.

Flow Characteristics of Liquid with Pressure-Dependent Viscosities in Microtubes

It is obvious that the pressure gradient alone, the axial direction in a pipe flow keeps constant according to the Haoen-Poiseuille equation. However, recent experiments indicated that the distribution of the pressure seemed no longer linear for liquid flows in microtubes driven by high pressure (1-30MPa). Based on H-P equation with slip boundary condition and Bridgman's relation of viscosity vs. static pressure, the nonlinear distribution of pressure along the axial direction is analyzed in this paper. The revised standard Poiseuille number with the effect of pressure-dependent viscosity taken into account agrees well with the experimental results. Therefore, the dependence of the viscosity on the pressure is one of the dominating, factors under high driven pressure, and is represented by an important property coefficient et of the liquid.

Generation of long, laminar plasma jets at atmospheric pressure and effects of flow turbulence

Long, laminar plasma jets at atmospheric pressure of pure argon and a mixture of argon and nitrogen with jet length up to 45 fi,Hes its diameter could be generated with a DC are torch by! restricting the movement of arc root in the torch channel. Effects of torch structure, gas feeding, and characteristics of power supply on the length of plasma jets were experimentally examined. Plasma jets of considerable length and excellent stability could be obtained by regulating the generating parameters, including are channel geometry gas flow I ate, and feeding methods, etc. Influence of flow turbulence at the torch,nozzle exit on the temperature distribution of plasma jets was numerically simulated. The analysis indicated that laminar flow plasma with very low initial turbulent kinetic energy will produce a long jet, with low axial temperature gradient. This kind of long laminar plasma jet could greatly improve the controllability for materials processing, compared with a short turbulent are let.

Pressure and thermal effects on elastic properties of single crystal alpha-Al2O3 from Monte-Carlo simulation

A rectangular structural unit cell of a-Al2O3 is generated from its hexagonal one. For the rectangular structural crystal with a simple interatomic potential [Matsui, Mineral Mag. 58A, 571 (1994)], the relations of lattice constants to homogeneous pressure and temperature are calculated by using Monte-Carlo method at temperature 298K and 0 GPa, respectively. Both numerical results agree with experimental ones fairly well. By comparing pair distribution function, the crystal structure of a-Al2O3 has no phase transition in the range of systematic parameters. Based on the potential model, pressure dependence of isothermal bulk moduli is predicted. Under variation of general strains, which include of external and internal strains, elastic constants of a-Al2O3 in the different homogeneous load are determined. Along with increase of pressure, axial elastic constants increase appreciably, but nonaxial elastic constants are slowly changed.

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.

Self-consistent kinetic modeling of low-pressure inductively coupled radio-frequency discharges

An efficient method for solving the spatially inhomogeneous Boltzmann equation in a two-term approximation for low-pressure inductively coupled plasmas has been developed. The electron distribution function (EDF), a function of total electron energy and two spatial coordinates, is found self-consistently with the static space-charge potential which is computed from a 2D fluid model, and the rf electric field profile which is calculated from the Maxwell equations. The EDF and the spatial distributions of the electron density, potential, temperature, ionization rate, and the inductive electric field are calculated and discussed. (C) 1996 American Institute of Physics.

The influence of Saffman lift force on nanoparticle concentration distribution

The lift force on a spherical nanoparticle near a wall in micro/nanofluidics has not received sufficient attention so far. In this letter the concentration of 200 nm particles is measured at 0.25–2.0 m to a wall in a microchannel with pressure-driven de-ionized water flow pressure gradient 0–2000 kPa/m . The measured data show the influence of the lift force on the nanoparticle concentration distribution. By introducing the Saffman lift force into the Nernst–Planck equation near a wall, we find that the lift force is dominant at the range of 2

FBG pressure sensor based on the double shell cylinder with temperature compensation

A novel fiber Bragg grating (FBG) pressure sensor based on the double shell cylinder with temperature compensation is presented. in the sensing scheme, a sensing FBG is affixed in the tangential direction on the outer surface of the inner cylinder, and another FBG is affixed in the axial direction to compensate the temperature fluctuation. Based on the theory of elasticity, the theoretical analysis of the strain distribution of the sensing shell is presented. Experiments are carried out to test the performance of the sensor. A pressure sensitivity of 0.0937 nm/MPa has been achieved. The experimental results also demonstrate that the two FBGs have the same temperature sensitivity, which can be utilized to compensate the temperature induced wavelength shift during the pressure measurement. (C) 2008 Elsevier Ltd. All rights reserved.

Influence of processing temperature on microstructure and microhardness of copper subjected to high-pressure torsion

Cu samples were subjected to high-pressure torsion (HPT) with up to 6 turns at room temperature (RT) and liquid nitrogen temperature (LNT), respectively. The effects of temperature on grain refinement and microhardness variation were investigated. For the samples after HPT processing at RT, the grain size reduced from 43 mu m to 265 nm, and the Vickers microhardness increased from HV52 to HV140. However, for the samples after HPT processing at LNT, the value of microhardness reached its maximum of HV150 near the center of the sample and it decreased to HV80 at the periphery region. Microstructure observations revealed that HPT straining at LNT induced lamellar structures with thickness less than 100 nm appearing near the central region of the sample, but further deformation induced an inhomogeneous distribution of grain sizes, with submicrometer-sized grains embedded inside micrometer-sized grains. The submicrometer-sized grains with high dislocation density indicated their nonequilibrium nature. On the contrary, the micrometer-sized grains were nearly free of dislocation, without obvious deformation trace remaining in them. These images demonstrated that the appearance of micrometer-sized grains is the result of abnormal grain growth of the deformed fine grains.

Instabilities in a non-transferred direct current plasma torch operated at reduced pressure

Using an oscilloscope, a high-speed video camera and a double-electrostatic probe system, the periodicity and amplitude of the fluctuations in arc voltage, jet luminance and ion saturation current of a plasma jet were monitored to investigate various sources of instabilities and their effects in a non-transferred dc plasma torch operated at reduced pressure. The results show that besides a 300 Hz main fluctuation inherited from the power supply, arc voltage fluctuation of 3–4 kHz with an amplitude less than 5% of the mean voltage was mainly affected by the total gas flow rate. The arc voltage fluctuation can affect the energy distribution of the plasma jet which is detectable by electrostatic probes and a high-speed video camera. The steadiness of energy transfer is also affected by the laminar or turbulent flow state of the plasma.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone). V: Estimation of temperature distribution during cure process

A numerical method to estimate temperature distribution during the cure of epoxy-terminated poly(phenylene ether ketone) (E-PEK)-based composite is suggested. The effect of the temperature distribution on the selection of cure cycle is evaluated using a suggested alternation criterion. The effect of varying heating rate and thickness on the temperature distribution, viscosity distribution and distribution of the extent of cure reaction are discussed based on the combination of the here-established temperature distribution model and the previously established curing kinetics model and chemorheological model. It is found that, for a thin composite (<=10mm) and low heating rate (<=2.5K/min), the effect of temperature distribution on cure cycle and on the processing window for pressure application can be neglected. Low heating rate is of benefit to reduce the temperature gradient. The processing window for pressure application becomes narrower with increasing thicknesses of composite sheets. The validity of the temperature distribution model and the modified processing window is evaluated through the characterization of mechanical and physical properties of E-PEK-based composite fabricated according to different temperature distribution conditions.

Spatial distribution of ciliates, copepod nauplii and eggs, Engraulis japonicus post-larvae and microzooplankton herbivorous activity in the Yellow Sea, China

The abundance of anchovy Engraulis japonicus larvae, >20 mum ciliates, copepod eggs and nauplii, and microzooplankton herbivorous activity were studied in the Yellow Sea in June 2000. Anchovy juveniles and larvae were found in only 6 of the 19 stations sampled. The ciliate communities were dominated by 2 species: Laboea strobila and Strombidium compressum. In the surface waters, the abundance of L. strobila ranged between 0 and 560 ind. l(-1). S. compressum only appeared at Stns 15 to 18 (20 to 3300 ind. l(-1)). L. strobila was found mainly in the top 20 m. The abundance of L. strobila was less than 50 ind, l(-1) in waters deeper than 25 m. S, compressum showed subsurface abundance peaks at the salinity abnormality. Tintinnids occurred occasionally with abundance lower than 100 ind. l(-1), The total ciliate abundance fell in the range of 40 to 3420 ind. l(-1). The ciliate biomass in the surface water and the water column ranged between 0,15 and 6.76 mug C l(-1) and 0.4 and 134.4 mg C m(-2), respectively, In the surface waters, the abundance of copepod eggs and nauplii ranged from 0,3 to 3.1 and 1,1 to 15.6 ind, l(-1), respectively. The average abundance of copepod eggs and nauplii in 4 depth (0, 5, 10 and 20 m) fell in the range of 0.2 to 2.8 and 1.0 to 29.4 ind. l(-1), respectively. As a food item of the E. japonicus post-larvae, the abundance of copepod nauplii and eggs appeared to be low. The abundance peaks of ciliate and E, japonicus post-larvae coincided. Although not found in the gut of E, japonicus post-larvae, aloricate ciliates might be ingested by first-feeding anchovy larvae, preventing initial starvation and prolonging the time to irreversible starvation. On the basis of dilution experiments with positive microzooplankton grazing rates, microzooplankton grazed at rates of 0 to 0.61 d(-1). Grazing pressure of microzooplankton on chlorophyll a standing stock (P-i) and potential chlorophyll a primary production (P-p) were 17 to 46% and 35 to 109% d(-1), respectively.