580 likes | 600 Views
International Conference: Scaling Up and Modeling for Transport and Flow in Porous Media Dubrovnik, Croatia, 13-16 October 2008. Multicomponent two-phase flow in porous media: Macro - kinetics of oscillatory regims. Mojdeh Rassoulzadeh – LEMTA Irina Panfilova – LEMTA/Schlumberger
E N D
International Conference: Scaling Up and Modeling for Transport and Flow in Porous Media Dubrovnik, Croatia, 13-16 October 2008 Multicomponent two-phase flow in porous media: Macro-kinetics of oscillatory regims Mojdeh Rassoulzadeh – LEMTA Irina Panfilova – LEMTA/Schlumberger Michel Panfilov – LEMTA Laboratoire d'Énergétique et de Mécanique Théorique et Appliquée (LEMTA – CNRS UMR 7563)
FLUIDE Subsurface waste storage Components : Gaz Water
FLUIDE Oil and natural gas Components : Gaz Oil
FLUIDE Oil and CO2 Components : Gaz Oil
PT-Diagram Gas Liquid Gas + Liquid
Classic systems Initial state Gas Liquid Gas + Liquid
Initial state Gas Liquid Gas + Liquid Retrograde systems
MAIN PROBLEM OF MULTICOMPONENT FLOW Non-equilibrium behaviour
MAIN PROBLEM OF MULTICOMPONENT FLOW Non-equilibrium behaviour • Oscillatory regimes
MAIN PROBLEM OF MULTICOMPONENT FLOW Non-equilibrium behaviour • Oscillatory regimes 2. Over-saturated zones
PROBLEM 1 : Oscillatory regimes
RETROGRADE GAS-OIL RESERVOIRS Liquid Gas Gas + Liquid
RETROGRADE GAS-OIL RESERVOIRS Theory flow rate composition
RETROGRADE GAS-OIL RESERVOIRS Field data flow rate composition
HYPOTHESES ON THE MECHANISM OF OSCILLATIONS Ganglion character of flow (V. E. Gorbunov, 1990) Each fluid becomes mobile only when it reaches its representative elementary volume (REV) Thermodynamic instability (V. Mitlin, 1990) Stability analysis of the compositional flow model shows that the system becomes instable when is the total mixture density, P is the pressure
OUR THEORY double phase transition: condensation coagulation of liquid internal evaporation internal gas evacuation
OUR THEORY condensation coagulation of liquid P leads to evaporation
OUR THEORY Phase diagram for the initial fluid Phase diagram for the secondary liquid aggregates P condensation liquid coagulation internal evaporation
Double phase transition Initial state Liquid Gas Gas + Liquid
Double phase transition Initial state Liquid Gas Gas + Liquid
Double phase transition Initial state Liquid Gas Gas + Liquid
Capillary condensation Double phase transition Initial state
Capillary condensation Double phase transition Initial state
Double phase transition Liquid coagulation
Double phase transition Liquid coagulation Liquid aggregate
Double phase transition Transition to the second phase diagram
Double phase transition Internal evaporation (boiling) Transition to the second phase diagram
Double phase transition Gas Evacuation
TOTAL COMPOSITION OF THE SYSTEM: 4 PHASES Classic phases 1 2 3 4
MODEL of DOUBLE PHASE TRANSITION Capillary condensation Minimisation of free Gibbs energy Coagulation Smoluchowski + effective media Evaporation Kinetics of Frenkel-Zeldovich Evacuation Gravity segregation + volume exceed mechanism
CAPILLARY CONDENSAION Pore-scale modeling Correlated capillary network Liquid aggregates 1 and dispersed condensate 2, 3
Results of modeling the liquid COAGULATION Dynamics of the averaged size of liquid aggregates
COAGULATON: Effective medium approach Comparison of the effective medium theory and the network simulations Mean vale of particle for power law probability of coagulation kinetic of coagulation
SECONDARY EVAPORATION (BOILING) Evaporation has 2 stages: A : formation and growth of germs of bubbles (Frenkel, Zeldovich) B : coagulation of bubbles • is the mass concentration of the aggregate Is the mass concentration of the boiling gas
EVACUATION: gravity segregation + volume exceed mechanism Internal exchange: formation of gas bubbles leads to the reduction of the liquid mass External exchange: geometrical “volume exceed” gravity-induced uplift of bubbles General kinetic for the external exchange
aggregation evaporation evacuation by volume exceed evacuation by segregation evaporation Volterra generalized model = mass of liquid aggregates = mass of interior gas
Phase portrait Rapid gas evacuation:
Phase portrait Rapid gas evacuation: CENTER
Stable Oscillations (case of rapid gas evacuation)
Phase portrait Slow gas evacuation:
Phase portrait Slow gas evacuation: FOCUS
Attenuating Oscillations System oscillation: 1 sigma theta 0,8 sigma+theta 0,6 0,4 0,2 t 0 0 50 100 150 200 (case of slow gas evacuation)
FOUR-PHASE MODEL: Numerical tests Volterra kinetics production well Total liquid Saturation FLOW Radial coordinate
PSEUDO THREE-PHASE MODEL - Mobile liquid is neglecting - Two-component system (light & heavy components)
LIQUID SATURATION CLASSIC MODEL Flow direction
LIQUID SATURATION CLASSIC MODEL
LIQUID SATURATION MODEL with DOUBLE PHASE TRANSITION
The macroscale oscillations – whether this is possible ?