Sensitivity Analysis Of The Dimensionless Parameters In Scaling A Polymer Flooding Reservoir

A set of scaling criteria of a polymer flooding reservoir is derived from the governing equations, which involve gravity and capillary force, compressibility of water, oil, and rock, non-Newtonian behavior of the polymer solution, absorption, dispersion, and diffusion, etc. A numerical approach to quantify the dominance degree of each dimensionless parameter is proposed. With this approach, the sensitivity factor of each dimensionless parameter is evaluated. The results show that in polymer flooding, the order of the sensitivity factor ranges from 10(-5) to 10(0) and the dominant dimensionless parameters are generally the ratio of the oil permeability under the condition of the irreducible water saturation to water permeability under the condition of residual oil saturation, density, and viscosity ratios between water and oil, the reduced initial oleic phase saturation and the shear rate exponent of the polymer solution. It is also revealed that the dominant dimensionless parameters may be different from case to case. The effect of some physical variables, such as oil viscosity, injection rate, and permeability, on the dominance degree of the dimensionless parameters is analyzed and the dominant ones are determined for different cases.

国内非牛顿流体力学进展

非牛顿流体;;本构方程;;流变测量;;流动稳定性;;收缩流动;;聚合物加工;;石油工业;;生物流体力学

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

Distributions of temperature, velocity and viscosity at symmetrical axis of cylindrical plume in earth's mantle with variable viscosity

The parameters at the symmetrical axis of a cylindrical plume characterize the strength of this plume and provide a boundary condition which must be given to investigate the structure of a plume. For Newtonian fluid with a temperature-and pressure-dependence viscosity, an asymptotical solution of hydrodynamic equations at the symmetrical axis of the plume is found in the present paper. The temperature, upward velocity and viscosity at the symmetrical axis have been obtained as functions of depth, The calculated results have been given for two typical sets of Newtonian rheological parameters. The results obtained show that the temperature distribution along the symmetrical axis is nearly independent of the theological parameters. The upward velocity at the symmetrical axis, however, is strongly dependent on the rheological parameters.

Pressure Drop Studies of Gas-Liquid Bubbly Flow in Vertical Upward Pipeline

This paper describes the experimental and theoretical studies of gas-liquid bubbly flow in vertical upward pipeline carried out at Institute of Mechanics, Chinese Academy of Sciences. Bubbly flow in a vertical pipe with a 3 m long and 5 cm inner diameter plexiglass pipe was experimentally investigated, and studies carried out on the relationship between superficial velocities of the liquid and gas phases and pressure gradient is described. The developed drift-flux model applied to gas-liquid bubbly flow is presented, and the results are compared against the experimental data measured by ours in air/water vertical pipes.

Simulation of Severe Slugging Formation and Control in Offshore Riser

In this paper, available elimination techniques are assessed. OLGA2000 software is used to simulate severe slugging formation mechanism in certain offshore riser. The simulation results show that pressure fluctuations of riser base and riser top is very large and severe slugging easily forms. Sensibility analysis shows that the measures and methods which include properly reducing pipe riser diameter, reducing water cut increasing terminal pressure, decreasing the height and inclination of riser and increasing GOR can eliminate or control severe slugging in riser pipe.

Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at Present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory. According to gas-liquid two-phase flow theory and beat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Bubble-to-Slug Transition for Low-Re Two-Phase Flows

A critical review on the mechanism and models on the bubble-to-slug transition of two-phase gas-liquid flows are presented in the present paper. It is shown that the most possible mechanism controlled the bubble-to-slug transition is the bubble coalescence. Focusing on the bubble-to-slug transition for the low-Re two-phase flow, a simple Monte Carlo method is used to simulate the influence of the initial bubble size on the bubble-to-slug transition. Some secondary factors, such as the liquid viscosity, the surface tension, and the relative slip between the two phases, are ignored in the present study. It is found that the locus of the dimensionless rate of collision is a universal curve. Based on this curve, it is determined that the bubble initial size can affect the phase distribution and flow pattern when its dimensionless value is in the range from 0.03 to 0.4. A simple relationship between the critical void fraction and the initial bubble size is proposed, which agrees very well with the experimental data.

Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory According to gas-liquid two-phase flow theory and heat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Numerical Solution for Polymer Transient Flows in a Circle Bounded Composite Formation

A new numerical model for transient flows of polymer solution in a circular bounded composite formation is presented in this paper. Typical curves of the wellbore transient pressure are yielded by FEM. The effects of non-Newtonian power-law index, mobility and boundary distance have been considered. It is found that for the mobility ratio larger than 1, which is favorable for the polymer flooding, the pressure derivative curve in log-log form rises up without any hollow. On the other hand, if the pressure derivative curve has a hollow and then is raised up, we say that the polymer flooding fails. Finally, the new model has been extended to more complicated boundary case.

Microfluidic Rheometry

The development and growth of microfluidics has stimulated interest in the behaviour of complex liquids in micro-scale geometries and provided a rich platform for rheometric investigations of non-Newtonian phenomena at small scales. Microfluidic techniques present the rheologist with new opportunities for material property measurement and this review discusses the use of microfluidic devices to measure bulk rheology in both shear and extensional flows. Capillary, stagnation and contraction flows are presented in this context and developments, limitations and future perspectives are examined. (C) 2008 Elsevier Ltd. All rights reserved.

Flow-induced vibrations of long circular cylinders modeled by coupled nonlinear oscillators

The dynamics of long slender cylinders undergoing vortex-induced vibrations (VIV) is studied in this work. Long slender cylinders such as risers or tension legs are widely used in the field of ocean engineering. When the sea current flows past a cylinder, it will be excited due to vortex shedding. A three-dimensional time domain model is formulated to describe the response of the cylinder, in which the in-line (IL) and cross-flow (CF) deflections are coupled. The wake dynamics, including in-line and cross-flow vibrations, is represented using a pair of non-linear oscillators distributed along the cylinder. The wake oscillators are coupled to the dynamics of the long cylinder with the acceleration coupling term. A non-linear fluid force model is accounted for to reflect the relative motion of cylinder to current. The model is validated against the published data from a tank experiment with the free span riser. The comparisons show that some aspects due to VIV of long flexible cylinders can be reproduced by the proposed model, such as vibrating frequency, dominant mode number, occurrence and transition of the standing or traveling waves. In the case study, the simulations show that the IL curvature is not smaller than CF curvature, which indicates that both IL and CF vibrations are important for the structural fatigue damage.

Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

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

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.

Non-relativistic conformal symmetries in fluid mechanics

The symmetries of a free incompressible fluid span the Galilei group, augmented with independent dilations of space and time. When the fluid is compressible, the symmetry is enlarged to the expanded Schrodinger group, which also involves, in addition, Schrodinger expansions. While incompressible fluid dynamics can be derived as an appropriate non-relativistic limit of a conformally invariant relativistic theory, the recently discussed conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry. This is explained by the subtleties of the non-relativistic limit.