CA2634757C - Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive - Google Patents

Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive Download PDF

Info

Publication number
CA2634757C
CA2634757C CA2634757A CA2634757A CA2634757C CA 2634757 C CA2634757 C CA 2634757C CA 2634757 A CA2634757 A CA 2634757A CA 2634757 A CA2634757 A CA 2634757A CA 2634757 C CA2634757 C CA 2634757C
Authority
CA
Canada
Prior art keywords
gas
correlation
reservoir
baseline
relative permeability
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CA2634757A
Other languages
French (fr)
Other versions
CA2634757A1 (en
Inventor
Mridul Kumar
Frederic Gadelle
Akshay Sahni
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chevron USA Inc
Original Assignee
Chevron USA Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chevron USA Inc filed Critical Chevron USA Inc
Publication of CA2634757A1 publication Critical patent/CA2634757A1/en
Application granted granted Critical
Publication of CA2634757C publication Critical patent/CA2634757C/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/10Locating fluid leaks, intrusions or movements

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

A method a system and a program storage device for predicting a property of a fluid from a subterranean reservoir containing heavy oil entrained with gas The methoo includes developing a baseline correlation between gas relative permeability k(rg) and gas saturation S(g), and capillary number N(c) calculated from a plurality of cells and at least one of critical gas saturation S(gc) and gas relative permeability k(rgro) based upon plurality of depletion rates The baseline correlation is then adjusted to comport with at least one of S(gc) and k(rgro) to produce a plurality of corresponding adjusted baseline correlations where gas relative permeabilities k(rg) are selected A reservoir simulation is then run utilizing the selected relative permeabilities k(rg) to predict a property of at least one fluid

Description

METHOD, SYSTEM AND PROGRAM STORAGE
DEVICE FOR RESERVOIR SIMULATION
UTILIZING HEAVY OIL SOLUTION GAS DRIVE
TECHNICAL FIELD
The present invention relates generally to methods and systems for reservoir simulation predicting the flow of fluids in an underground reservoir, and more particularly, to enhancing reservoir performance forecasting by accounting for fluid flow effects due to heavy oil solution gas drive.
BACKGROUND OF THE INVENTION
Reservoir simulation is used to predict the flow of fluids in an underground reservoir.
The fluid flow may include oil, gas and water. Such reservoir forecasting is important in reservoir management and estimating the potential recovery from a reservoir.
Reservoir simulation is well known. throughout the oil industry and in the scientific literature. A good primer on the principles behind reservoir simulation is K.
Aziz and A. Settari, Petroleum Reservoir Simulation, Elsevier Applied Science Publishers, London (1979). Another description of how reservoir simulation is generally performed is described in U.S. Patent No. 6,052,520 to Watts III et al.
The following are general steps taken in a conventional reservoir simulation.
First, a reservoir is selected for which the rock and fluid properties are to be modeled and simulated. The reservoir=is modeled and discretized into a plurality of cells.
Nonlinear governing equations are constructed for each cell, generally in the form of finite difference equations, which are representative of properties of rocks and fluids in the reservoir. Examples of rock properties include porosity, capillary pressure, and relative permeability for each phase of fluid (oil, water, gas.) Examples of fluid properties include oil viscosity, oil formation factor (80), and pressure, temperature, and saturation in each of the cells. Nonlinear terms in these equations are linearized to - - =

arrive at a set of linear' equations for each timestep of the simulation.
These linear equations can then be solved to estimate solutions for unknowns such as pressure and saturation in the cells. From these values of pressure and saturation other properties can be estimated including the overall production of oil, gas and water from the reservoir in a timestep. The aforementioned steps are repeated over many such timesteps to simulate fluid flow over time in the reservoir.
One of the key properties needed in reservoir simulation is the permeability of a rock to flow. Absolute permeability K is a measure of a rock's ability to transmit flow and can vary greatly throughout a reservoir and surrounding formations. When gas, oil and water move through porous rock, they do not move at equal velocities.
Rather, the fluids compete with one another. Relative permeability, kr, is the ratio of the effective permeability, Ice, when more than one fluid is present, to the absolute permeability K.
Effective permeability ice is the measured permeability of a porous medium to one fluid when another is present. The relationship between relative permeability kr and saturation S depends on the reservoir rock and fluid and may vary between formations. Also, the relative permeability kr depends on the relative proportion of the fluids present, i.e. fluid saturations.
FIG. 1 illustrates a typical relative permeability kg versus saturation Sg curve for gas.
Gas cannot flow at any appreciable rate until gas saturation reaches a minimum threshold value. Looking to FIG. 1, this threshold value is referred to as critical gas saturation 4, and begins at a value of approximately 0.03 or about 3%
saturation. At the other end of the curve is an endpoint relative permeability kr% rõ which is the gas relative permeability value krg at which movement of residual oil remaining in the rock is minimal. As reservoir rock will always contain a minimal amount of residual oil, gas saturation cannot reach 100%. The total percentage of saturation must add up to 100%. In this case, there is a maximum 76% gas saturation Sr and 24 %
residual oil saturation Sors . As seen in FIG. 1, the maximum relative permeability, k occurs at a saturation of approximately 0.76 with kr= 0.40. These values of S.
and
-2-.

j4,0 shall be referred to as endpoint baseline values for gas saturation S and relative permeability krg=
Ideally, relative permeability curves are developed through laboratory experiments on core samples taken from reservoirs for which reservoir simulation is to be performed.
For example, displacement tests may be used to develop the relative permeability krg vs. saturation Sg curves. Such tests are well known. Particularly well known displacement test procedures are described in E.F..Johnson, D.P. Bossier, and V.O. Naumann, Calculations of Relative Permeability from Displacement Experiments, Trans. Am. Inst. Mining Engineers, Volume 216, 1959, pp. 370-378 and S.C. Jones and W.O. Roszelle, Graphical Techniques for Determining Relative Permeability from Displacement Experiments, Journal of Petroleum Engineering, Volume 30, pp. 807-817 (1978). These displacements experiments are usually conducted at slow depletion rates as it is commonly accepted that permeability curves are generally independent of how fast gas flows through reservoir rock.
Alternatively, if core samples are not available, the relative permeability kg versus saturation Sg curves can be theoreticially created. For example, the curves may be developed from comparable analogue reservoirs.
Once relative permeability krg versus saturation Sg curves have been obtained, then the relative permeabilities kg to be used in a reservoir simulation can simply be obtained from these curves assuming saturations Sg in the cells of the reservoir model are known. The saturations S8 are generally known either from initial conditions established at the beginning of a simulation, from the last timestep in the simulation = or else from calculations within an iteration in a timestep.
The production of heavy oil is initially driven primarily by oil pressure.
Heavy oil may be considered to include oil having an API gravity 20 or less.
Significant quantities of gas are often entrained within the heavy oil while under high reservoir = pressures. After sufficient production of heavy oil from a reservoir, the pressure in
- 3 -
4 PCT/US2006/049095 portions of the reservoir may drop below the bubble point pressure. At this pressure, gas readily comes out of solution from the heavy oil. Once sufficient gas has been released from the oil, the gas is believed to form a continuous phase and gas can flow through the reservoir and the rate of production of gas is significantly enhanced. As indicated above, the saturation Se at which there is an initiation of gas flow is referred to as the critical gas saturation or S. . FIG. 11 shows a graph of cumulative gas produced from a core sample versus time in minutes. The breakpoint in the curve shown there represents See .
Tests have shown that the amount of oil recovery from a heavy oil reservoir is dependent upon the rate of depletion of the reservoir. Often higher rates of depletion will lead to an overall enhanced oil recovery. As the mechanisms of heavy oil solution gas drive are not well understood, reservoir simulators typically utilize static gas relative permeability kra versus saturation Sg curves, such as the one seen in FIG. 1, which are independent of fluid flow or depletion rates. Once these curves are developed for respective types of rock which are to be modeled, the curves will remain the same (i.e., endpoints of S c and kg ,.0 remain fixed) throughout the reservoir simulation regardless of the rate of flow through the reservoir cells. Such assumptions that permeability-curves are static for general reservoir simulation of hydrocarbon bearing subterranean formations containing non-heavy oil are generally satisfactory.
However, in the case of heavy oil, non-equilibrium solution gas drive ("Foamy Oil") is a significant production mechanism affecting critical gas saturation Se, and oil recovery. Currently, understanding of heavy oil solution gas drive at all scales (pore, core and field) is limited. Conventional reservoir simulators fail to accurately account for this solution gas drive in forecasting fluid flow in a reservoir. This is a significant shortcoming often resulting in forecasts which underestimate heavy oil production.
The present invention overcomes this shortcoming by accounting for the effects of heavy oil solution gas drive.

SUMMARY OF THE INVENTION
A method of predicting a property of at least one fluid in a subterranean reservoir containing heavy oil entrained with gas is disclosed. For example, the property might include the overall production of fluids from the reservoir, i.e., oil, gas and water. The prediction is made using a reservoir simulator which uses a reservoir model having a plurality of cells representative of the reservoir. For at least some of the cells and for at least some of the iterations of the reservoir simulation, gas relative permeability krg is dependent upon the local fluid velocities vo in the cells.
In a preferred embodiment of this method, a baseline correlation is developed for gas relative permeability I% versus gas saturation Sg, typically based on displacements tests performed at slow depletion rates. Next, a capillary number Arca dependent correlation is developed between at least one of, and most preferably, both of critical gas saturation So and capillary number Nca and endpoint of gas relative permeability Icgro and capillary number ka. Non-limiting examples of how this correlation may be expressed include, by way of example and not limitation, using a mathematical equation which describes a curve or by creating a corresponding look-up table.
These experimentally derived capillary number Ncadependent correlations can then be used, in conjunction with reservoir simulation, to capture the effects that heavy oil solution gas drive and depletion rates have on the production of heavy oil and gas entrained therein. Capillary numbers Ne are calculated for a plurality of cells in the reservoir model representative of the subterranean reservoir for which fluid properties are to be simulated. Sgc and/or krg,.õ values are selected from the capillary number dependent correlations based upon the capillary numbers Nc calculated for the cells.
Adjusted baseline correlations are then developed. For example, the original endpoints of the baseline curve, i.e. Sg , and k, .o,õ are replaced with the new capillary number dependent So and kgro values and the curve therebetween adjusted, such as
- 5 -by linear scaling. FIG. 2 suggests that an adjusted baseline curve can be developed by changing the original endpoint values S. and kr gõ, to other values of Sgõ, and kgõ
which are based, in part, upon the velocity of oil vo flow through the cells.
corresponding adjusted baseline correlations. These relative permeabilities kg are then used in a reservoir simulation to predict a property of at least one fluid in a subterranean reservoir containing heavy oil entrained with gas. This predicted property may be the production of oil, water or gas. Preferably, once saturation 5', in adjusted baseline correlation for that cell is fixed for the remaining simulation time-steps. This fixing of the adjusted baseline correlation once gas begins to flow assists in maintaining stability during the solution of the system of equations modeling the reservoir.
One or both of the capillary dependent correlations of S ge or k can be used in adjusting the baseline correlation to come up with an adjusted baseline correlation.
These adjusted baseline correlations, through the use of the capillary numbers N0, capture the effects that the depletion rate/fluid velocity flow and viscosity have on necessary, it is possible to theoretically predict what such capillary number dependent correlations should be.
Relative permeabilities krg can be selected which are dependent upon capillary numbers N0 calculated at the beginning of a time step in a reservoir simulation.
Alternatively, the capillary numbers N0 can be repeatedly calculated throughout iterations in a timestep to provide constant updating of relative permeability curves
- 6 -permeability curves of a cell is preferably stopped once the saturation Sg in a cell remains at or above the critical gas saturation Sg, during simulation.
It is an object of an aspect of the present invention to enhance reservoir performance forecasting by better accounting for fluid- flow effects due to heavy oil solution gas drive than in conventional reservoir simulators thereby improving the predictive capability of reservoir simulations involving heavy oil flow in subterranean formations which can lead to improved reservoir management strategies.
It is another object of an aspect to experimentally determine values for critical gas saturations 4, and/or for endpoint of gas relative permeability kgõ for a core sample at a number of different depletion rates and correlate these values against capillary numbers N,,, to create capillary number dependent correlations. These capillary number N ca dependent correlations can be used in conjunction with a reservoir model, and calculated capillary numbers A/c. calculated during a reservoir simulation, to more accurately estimate relative permeabilities kg to be used in the reservoir simulation of heavy oil.
According to another aspect, there is provided a method of predicting a property of at least one fluid in a subterranean reservoir containing heavy oil entrained with gas, the method comprising:
(a) utilizing a baseline correlation of gas relative permeability kg versus gas saturation S, in a reservoir simulator;
- 7 -(b) utilizing a capillary number dependent correlation between at least one of critical gas saturations Sõ and capillary numbers Nõ and endpoint gas relative permeabilities kgõ and capillary numbers Nar in the reservoir simulator;
(c) calculating capillary numbers N, for a plurality of cells in a reservoir model representative of the subterranean reservoir;
(d) adjusting the baseline correlation to comport with at least one of S, and lc,õ selected from the capillary number dependent correlation of step (b) using the capillary number N, calculated in step (c) to produce a plurality of corresponding adjusted baseline correlations;
(e) selecting relative permeabilities krg for the plurality of cells from the corresponding adjusted baseline correlations of step (d); and (1) running a reservoir simulation utilizing the selected relative permeabilities Icrg of step (e) to predict a property of at least one fluid in the subterranean reservoir containing heavy oil entrained with gas.
According to a further aspect, there is provided a method for simulating the flow of heavy oil in a subten-anean reservoir, the method comprising:
creating a reservoir model representative of a subterranean reservoir for which fluid flow is to be simulated, the reservoir model including a plurality of reservoir cells;
- 7a -determining velocity-dependent relative permeabilities for the reservoir cells that account for a velocity of the fluid flow through the reservoir cells, the velocity-dependent relative permeabilities for the reservoir cells being determined using a gas relative permeability versus gas saturation correlation constructed by scaling a baseline gas relative permeability versus gas saturation correlation to at least one of a critical gas saturation and an endpoint gas relative permeability; and running a reservoir simulation utilizing the velocity-dependent relative permeabilities to simulate the flow of heavy oil in the subterranean reservoir.
According to another aspect, there is provided a system for simulating the flow of heavy oil in a subterranean reservoir, the system comprising:
a program storage device carrying computer instructions to perform a method of reservoir simulation; and a reservoir simulator for executing the computer instructions of the program storage device to perform the method of reservoir simulation, the method of reservoir simulation comprising the steps of (a) utilizing a baseline correlation of gas relative permeability icrg versus gas saturation 5', ;
(b) utilizing a capillary number dependent correlation between at least one of critical gas saturations Sg and capillary numbers Nõ and endpoint gas relative permeabilities and capillary numbers Nõ;
- 7b -(c) calculating capillary numbers /V, for a plurality of cells in a reservoir model representative of the subterranean reservoir containing heavy oil;
(d) adjusting the baseline correlation to comport with at least one of Sg, and lc, selected from the capillary number dependent correlation of step (b) using the capillary number Nc calculated in step (c) to produce a plurality of corresponding adjusted baseline correlations;
(e) selecting relative permeabilities krg for the plurality of cells from the corresponding adjusted baseline correlations of step (d); and (f) simulating the flow of the heavy oil in the subterranean reservoir utilizing the selected relative permeabilities krg of step (e).
According to a further aspect, there is provided a method for predicting a property of at least one fluid in a subterranean reservoir containing heavy oil entrained with gas, the method comprising:
(a) inputting a baseline correlation of gas relative permeability krg versus gas saturation Sg that defines a baseline critical gas saturation Sg, and a baseline endpoint gas relative permeability krg, into a reservoir simulator;
(b) calculating at least one of an updated critical gas saturation Sge and an updated endpoint gas relative permeability k,,.0 based upon a plurality of depletion rates;
- 7c -(c) adjusting the baseline correlation to comport with at least one of the updated critical gas saturation Sg, and the updated endpoint gas relative permeability krgro;
(d) selecting relative permeabilities kg for the plurality of cells from the adjusted baseline correlation of step (c); and (e) running a reservoir simulation utilizing the selected relative permeabilities krg of step (d) to predict a property of at least one fluid in the subterranean reservoir containing heavy oil entrained with gas.
BRIEF DESCRIPTION OF THE DRAWINGS
These and other objects of aspects, features and advantages of the present invention will become better understood with regard to the following description, pending claims and accompanying drawings where:
FIG. 1 shows a conventional gas relative pen-neability ki g versus saturation Sg curve;
FIG. 2 depicts adjusting the conventional curve of FIG. 1 by modifying the endpoints of S'g, and k , gõ to coincide with values of Sg, and k, g,, selected from capillary number dependent correlations of Sg, versus Na, and k Kw versus Nea;
- 7d -FIG. 3 shows a flowchart of steps taken in a preferred embodiment of the present invention for carrying out reservoir simulation which utilizes gas relative permeabilities krg which are dependent upon local velocities võ of fluid flow in cells;
.. FIG. 4 shows a schematic drawing of an experimental setup used to determine gas saturation Sg from core and sandpack samples;
FIG. 5 depicts a graph of average sandpack pressure and pressure differential versus time across a sandpack sample in a fast depletion experiment;
FIG. 6 illustrates a graph of cumulative oil and gas produced in the fast depletion experiment of FIG. 5;
FIG. 7 shows a graph of average sandpack sample pressure and effluent density for a .. slow depletion experiment;
FIG. 8 depicts a graph of average sandpack sample pressure and cumulative oil produced for a slow depletion experiment;
.. FIG. 9 shows a graph of oil recovery as a function of average pore pressure for sandpack experiments at depletion rates of 0.3 and 0.03 cc/min.;
FIG. 10 is a graph of oil recovery as a function of average pore pressure for core experiments at depletion rates of 0.082, 0.08, and 0.002 cc/min, respectively;
and FIG. 11 is a graph of cumulative gas produced (measured) and cumulative solution gas produced (calculated) vs. time.
- 8 -DETAILED DESCRIPTION OF THE INVENTION
I. Introduction The present invention accounts for the effects of heavy oil solution gas drive, and more particularly, for the effects that the rates of fluid depletion have on heavy oil production. Velocity or depletion rate dependent relative permeability values krg are utilized in a heavy oil reservoir simulation to provide for more accurate reservoir simulation forecasts than are achieved with conventional reservoir simulation.
In a preferred embodiment, capillary numbers N, which are dependent on oil velocities va , are calculated for reservoir cells. These capillary numbers N.
are used to adjust baseline relative permeability correlations to account for the velocity or depletion rate effects on relative permeability krg. In this preferred embodiment, capillary number N dependent critical gas saturations Sgc and/or endpoint relative permeabilities kwõ correlations are first developed, preferably based on laboratory experiments. Then values of Sgc and/or kw, , corresponding to the capillary number N c calculated for a cell, are used to adjust the baseline relative permeability correlation for that cell. Relative permeability krg values are then selected from these capillary number adjusted baseline relative permeability correlations based upon the saturations Sg in the cells.
FIG. 3 provides an exemplary flowchart of steps which may be used to implement the heavy oil solution gas drive reservoir simulation of the present invention. In step 100, a baseline correlation is created between krg and Sg. Correlations are then developed between Sg, and N and/or kw, and Ncõ in step 110. For a number of cells in a reservoir model, capillary numbers N c are calculated in step 120. For each of these cells, adjusted baseline correlations between krg and Sg are established in step 130 which are dependent upon N c and the correlations developed in step 110. Gas relative permeabilities krg are then selected in step 140 for each of the cells from the adjusted baseline correlations between kg and Sg using saturation Sg values from the
- 9 -cells. These capillary number dependent permeabilities krg are then used in step 150 in a reservoir simulation to predict properties of fluid flow in the reservoir model.
A description of an exemplary test method for establishing correlations between Sgc and Nca and between k rgre and Na will be described. Then, modifications will be described which are made to a conventional reservoir simulator to incorporate the depletion rate/capillary number dependent Sg, and/or krg,õ correlations for selecting relative perrneabilities krg when conducting a reservoir simulation.
II. Establishing Correlations A. Baseline Gas Relative Permeability Ictg vs. Saturation Sg Correlations Correlations between gas relative permeability krg and saturation Sg are established so that relative permeability values kg can be utilized by a reservoir simulator based upon known saturations values Sg in cells of a reservoir model. Ideally, these correlations are experimentally developed from core samples from the reservoir for which the reservoir simulation is to be performed. Alternatively, representative sand packs and/or synthetic oil may also be used to develop the correlations. The preferred methods to establish these baseline correlations are the methods of Johnson, Bossier, and Naumann or else the method Jones and Roszell, which were cited above in the background section and are well known to those skilled in establishing permeability curves. Alternatively, there are many other well known schemes for establishing gas relative permeability krg versus saturation Sg curves for reservoir rocks and fluids.
Typically, krg is going to be a function of Sg. For practical reasons, one often normalizes the gas saturation used in the krg correlation. One such normalization is described in Eqn.12. Such normalization allows the simulator to readily evaluate krg for changing end-points (e.g., Sg,, and S.) If core samples are not available, then the correlations between relative permeability kg and saturation Sg can be theoretically estimated. As a non-limiting example, an
- 10-analogous formation maybe used to initially establish baseline curves. Non-limiting examples of correlations may take several forms such as curves, mathematical expressions, 'look-up tables, etc.
FIG. 1 is an exemplary baseline curve or correlation of gas relative permeability krg versus saturation Sg. A baseline value for Sg , is shown at about 0.03 or 3%.
Above this value, it is expected that gas will begin to flow freely rather than being primarily trapped within the the porous medium. The maximum gas saturation Sg is about 76% with there being a 24% saturation of residual oil saturation S. It is assumed there is very little presence of water for this example. At the maximum gas saturation Sg= 76%, the maximum gas relative permeability kr1,, is approximately 0.4%.
B. Correlations Between Sgõ vs. Nca and kõ,, vs. Nca Laboratory experiments were conducted at various depletion rates to establish Sg, vs. N0 and krg,õ vs. Nca correlations. Sg, is obtained in a method to be described below. N ca is calculated using Eqn. (8) below. From the experiments and history matching using reservoir simulations on core or sandpack samples, values of Sge, k,r.õ, and Nct, for each depletion rate were obtained. Then correlations between Sg, and ,Aca and between krõ,.,, and A fcc, were obtained by curve fitting the Sg, , krgra and Na data. History matching of production data on the core samples may be used to enhance the accuracy of the correlations.
1. Live Oil Preparation =
Live oil was prepared by combining unfiltered dead oil and methane. The water content of the oil was negligible. PVT (Pressure, Volume, and Temperature) data: Gas-Oil-Ratio (Rd), Oil Formation Volume Factor (BO and Gas Formation Volume Factors (Bg), were determined through a combination of experiments
- 11 -(constant composition expansion, flash, density measurement) and tuning of equation of states. Live oil viscosity was measured in a capillary viscometer (ID ------0.05 in) at reservoir temperature. Table 1 lists relevant properties of the live oil at 178 F.
Bubble Point Pressure (Psia) 1350 Solution GOR (cc/cc) 20 Bo at Bubble Point Pressure 1.0918 Live Oil Viscosity (cp) 240 Dead Oil Viscosity (cp) 1300 Properties Of Crude Oil 2. Depletion Experiments Depletion experiments were conducted at constant depletion rates in either a horizontal 80-cm long sandpack or in a 29-cm horizontal composite core (4 plugs).
The sand used in the sandpack experiments was clean Ottawa sand ranging in size from 75 to 125 .tn. The sand was packed in a specially made Viton sleeve equipped with pressure ports. The sandpack and composite core porosities were measured with a helium porosimeter. Sandpack and composite core properties are listed in Table 2:
Composite Sandpack core Temperature, F 178 178 Length, cm 80 29 Diameter, cm 5.04 5.04 Overburden Pressure, psia 2050 2050 Porosity 0.33 0.27 Pore Volume, cm3 560 162 Live Oil Permeability, md 2000 1850 Range of Depletion Rates, cmi/min 0.002 to 0.3 0.0003 to 0.03 Sandpack And Composite Core Properties
- 12 -The depletion rate was controlled using one or two ISCO pumps operating in a refill mode. FIG. 4 shows a schematic of the experimental set-up. During the depletion, the pressure (inlet, outlet, and at several points along the core), the production of oil and gas, and the density of the effluent was monitored. The coreholder was placed in a Siemens Somatom HiQ CT scanner to monitor spatial and temporal gas saturation.
FIG. 4 shows Live oil (410), Inlet side (412), Float Head Teflon Spacer (reference for CT Scanner)(414), Overburden Pressure (416), Pressure Transducers (418), Outlet Side (420), Fixed Head (422), and Gas Collection (cylinder under vacuum) (424).
3. Procedure The dry sandpack was initially CT (Computer Tomography) scanned at reservoir conditions (i.e., under overburden stress and at temperature). The core was then flushed with CO2, evacuated and saturated with kerosene at a back pressure of-psia. The sandpack (or composite core) permeability was measured with kerosene at several flow rates. The kerosene-saturated sandpack was also CT-scanned. The sandpack porosity was calculated using the wet and dry CT-scans and CT number of air and kerosene. Live oil was then slowly injected into the core to displace the kerosene. Permeability of the sandpack was also measured with live oil at several flow rates. The live-oil injection rate was then reduced so that the differential pressure across the core was less than 2 psi.
The live-oil saturated sandpack was CT-scanned to record initial conditions.
Depletion was started at a pressure of ¨1500-1700 psia (about 150-350 psi above the bubble point pressure). The inlet valve was closed and the downstream Isco pump A
- 13-was operated at a constant withdrawal rate. After a given depletion time, the pumps were switched and Isco pump B withdrew fluids while Isco pump A delivered oil and gas into the collection system. The pump cycle was repeated until the outlet pressure decreased to about 200 psia. Pressures, temperatures and fluid accumulation in the collection system were continuously recorded using conventional data acquisition software. The density of the produced fluid was continuously measured using an in-line density meter. The sandpack was also periodically scanned to determine directly gas saturation, Sg , as a function of time and position.
- 13a -4. CT-Scanning A Siemens Somatom HiQ CT scanner was used to monitor spatial and temporal gas saturation. This third generation CT-scanner has 768 stationary detectors and a rotating X-ray source. Scans were conducted at 133 kV and the scan time was 2.7 seconds. The voxel size was approximately 0.625 mm3 for a scan thickness of mm and the uncertainty in saturation measurement was +/- 1.5 saturation units.

Scan thicknesses of 10 mm and/or 5 mm were acquired.
=
10 5. Results During the course of experiments pressure information along the core and at the closed core inlet and open core outlet, the amount of oil and gas produced, the effluent density and gas saturation (via the CT-scanner) were acquired. The typical responses observed during an experiment are shown in FIGS. 5 and 6. FIG. 5 shows . the average sandpack pressure and pressure differential across the sandpack during a fast depletion experiment. FIG. 6 illustrates the cumulative oil and gas produced during a fast depletion experiment.
While not wishing to be held to a particular theory, it is believed that at an early time, production is through oil and formation expansion only (there is no free gas in the system) and the pressure falls rapidly. At the (apparent) bubble point pressure, gas bubbles start to nucleate. As the pressure decreases below the bubble point pressure, gas bubbles slowly grow in size and oil production is dominated by gas expansion. As can be seen from FIG. 5, the rate of pressure decrease was significantly reduced. Oil was the only moving phase and the gas collected was by liberation of dissolved gas in the collection system. At the critical gas saturation Sg, , gas bubbles are connected throughout the sandpack and gas starts to flow freely. Note that there is a significant increase in gas production while the oil production tapered off (see the sharp break in the cumulative gas production plot at ¨ 270 minutes).
For the slower depletion rate experiments in the sandpacks and for the core experiments, the effluent density was also measured. FIGS. 7 and 8 show typical
- 14 -responses which were observed with this instrument. FIG. 7 illustrates the average sandpack pressure and effluent density for a slow depletion experiment. FIG. 8 depicts the average sandpack pressure, cumulative oil produced (collected in the separator and inferred based on the effluent density) for a slow depletion experiment.
6. Rate Effect The main effect observed during the depletion experiments was that oil recovery is highly sensitive to the depletion rate. This phenomenon was observed with both large sandpack experiments (FIG. 9) and small core experiments (FIG. 10).
FIG. 9 illustrates oil recovery as a function of average pore pressure (sandpack experiments ¨ rates = 0.3 and 0.03 cc/min). FIG. 10 shows oil recovery as a function of average pore pressure. (Core experiments ¨ rates = 0.082, 0.08, and 0.002 cc/min.) In addition to the rate effect, note that the overall oil recovery observed in these experiments is quite large (up to ¨ 30% 00IP). Such high recovery and this dependency on depletion rates can not be readily explained by traditional physics.
Moreover, this phenomenon is not modeled properly with current commercial simulators.
7. Data Analysis ¨ Sg And So Determination The critical gas saturation So is the saturation at which the cumulative gas produced starts to increase significantly. FIG. 11 shows the cumulative gas produced (measured) and cumulative solution gas produced (calculated) vs. time. The critical gas saturation So can also be determined based on the effluent density.
With the set-up described in FIG. 4, there are several ways to determine the gas saturation:
(1) direct in situ measurement with the CT-scanner;
- 15-(2) material balance using the amount of fluids collected in the collection system;
and (3) material balance using the density of the effluent stream.
Methods 2 and 3 require the use of PVT data (namely formation volume factor and density as a function of pressure).
Material balance:
Sx = 1 (1) =
(N ¨ N p ) X Bo = ________________________________________________________________________ (2) N x Boi x(1¨ c f (Pi ¨P)) where N is the oil in place (stb) at the beginning of the experiment and at pressure P
Np is the cumulative oil produced (stb) at pressure P (Np is measured with the collection system), Bo and B01 are the oil formation volume factors at P and P

respectively and cf is the rock or sandpack compressibility (I/psi).
Above the bubble point, oil is produced through oil and formation expansion only.
That is Bo.
N p = (co + c x (Pi ¨ x N
(3) Bo where the oil compressibility is given by B0 ¨B01 C0¨ (4) Boi x (Pi ¨ P) With c, known, the sandpack and composite core compressibility are calculated using Eqn. (3).
- 16-As noted above, Np is measured through the collection system. Alternatively, the amount of oil produced can be based on the effluent density, p,ff :
N = depletion _rate x P = ¨ Pg 1 x At x + Np-1 (5) Po ¨ Pg Both porosity and gas saturation can be calculated using the CT-scanner.
Porosity is given by CTsaturated _ core CTdry _core (6) CTliq ¨ CT gas where CTsaturated core is the CT number for the sandpack saturated with kerosene (at initial pressure), and CTd,y_core is the CT number of the sandpack saturated with gas.
CTN, and CTga, are the CT numbers for kerosene and air, respectively.
Similarly, the gas saturation is obtained with the following equation:
CT p C -Tsaturated _core S R = ____________________________________________________________________ (7) Crdry _core CTsaturated _core where CT,,, is the CT number measured during the depletion (at pressure P), CTsalurated core is the CT number for the sandpack saturated with live oil (at initial pressure), and CTdry core is the CT number of the sandpack saturated with air and at initial pressure.
8. Data Analysis ¨ Capillary Number Calculation For each experiment, the average capillary number (N,,,) was calculated using the pressure differential recorded during the depletion. The capillary number can be
- 17 -calculated in several ways. In this preferred embodiment, the following formula was used:
K x AP
Ca _______________________________________________________________________ (8) o- x L
where K is the permeability of the core or sandpack, cr is the gas-oil surface tension (estimated to be 80 dynicm for the oil used in the experiment), L is the sandpack length, and AP is the pressure differential observed before the gas is becomes mobile.
9. Data Analysis ¨ Sgc and crgro as a function of A r ce, Based on the above analysis, Sg, is plotted as a function ofN,., for all the available experiments. The data is then curve fit, preferably, using an exponential function (Eqn. (9)) to interpolate/extrapolate the missing data. The coefficient "a"
and exponent "b" values are specific to each oil/rock system.
By way example, and not limitation, the preferred mathematical correlations between Sg, and k rgro as functions of N,,, are as follows:
S = a = log(i \ r.õ)+ b (9) and krgro = C = (.1\ T )d (10) Kragro and kr ,,,õ are "conventional" critical gas saturation and end-point of gas relative permeability values, respectively, as described above in the background and as shown in FIG. 1.
- 18 -Reservoir simulations conducted on core samples at various depletion rates are used to determine the values for k rg,õ . For each simulation run, the critical gas saturation Sg, is known, so this endpoint on a gas relative permeability krg versus saturation Sg is known. Various estimates are made for the other endpoint of the curve Icrgõ
. A trial and error method is then used to determine which estimated value of k rgro matches the experimental production output from the core sample at a particular depletion rate.
This history matching of experimental production results with simulated runs is used to determine kg, at a number of depletion rates, which correspond to N õ
values.
These values of krg,õ versus Nõ are then curve fit to arrive at a capillary number dependent correlation. Most preferably, this correlation is in the form of Eqn. (10) with values of "c" and "d" being determined.
III. Reservoir Simulation Utilizing Heavy Oil Solution Gas Drive Functional forms of Sg, and k rgõ vs. Na , obtained from experimental data, are implemented in this exemplary embodiment, preferably, using a modified implicit algorithm in a reservoir simulator. By way example, and not limitation, the preferred forms for Sg, and krg, are input as functions of Nõ using Eqns. (9) and (10) from above. The parameters a, b, c and d are user's input to the reservoir simulator. Note in function of c, d, and N . In the preferred embodiment of this invention, the following are default values: a = 104; b = 1.0; c = 104 and d = 1. Ideally, the calculated S and k p.m values are limited to user's specified maximums and minimums, respectively. For example, maximum Sg, = 0.1 and minimum value of krgm = le may be used. Since IV, is directional, Sg, and k rgõ are calculated for each cell face and thus are directional too.
To reduce oscillation and convergence problems, a modified implicit algorithm of the preferred embodiment is implemented to calculate Sg, and kg, . When the gas-phase
- 19 -is not mobile, i.e., saturation Sg 5_ Sgõ Sge and krgro are calculated, for example, using Eqns. (9) and (10), respectively. When the gas-phase is flowing, Sgõ and k,8,0 become invariant ¨ neither increase nor decrease. Their values are calculated using the capillary number N, at the beginning of the time-step when the gas-phase becomes mobile and fixed for all remaining time-steps.
A. Calculation of Cell Capillary Numbers N, In this preferred exemplary embodiment, a modified expression for capillary number N, is preferably incorporated into the reservoir simulator using the following expression:
k k = V(Põ ¨ põgD) ' 0 _ = (11) crag aog where crõg is oil-gas interfacial tension, K is rock permeability, (1)0 is oil-phase potential, Po is the change in pressure across a face of a cell, po = density of oil, g= gravitational constant, and D = change in depth from a datum.
This modified definition of /V, leaves out oil relative permeability in the equation.
Since N, is ideally computed implicitly, this greatly simplifies the calculation of derivatives for gas relative permeability (k,.9) as a function of primary variables during Jacobian generation. Also, the potential gradient in the N calculation is directional and is based on the gradient across the face of the two adjacent grid blocks. For each Newton iteration, a capillary number N, is calculated for each grid-block face. In a 3-D model, there will be six directional N, for each grid block.
Each N, corresponds to one of the six values at the cell faces. The use of directional AT, results in a Jacobian that can be easily solved by conventional linear equation solvers. For wells, in this preferred embodiment, an averaged AT, from all grid-block faces is calculated.
-20 -B. Adjusting Baseline Relative Permeability Correlations Each cell is assigned a particular rock type or fades. Each of these rock types or facies corresponds to particular baseline gas relative permeability krg vs.
saturation Sg curve, such as the one shown in FIG. 1. These respective baseline curves are adjusted for each respective cell. This is accomplished for each cell by replacing the original values of Se!, and kr with capillary number dependent values of Sgc and krgõ calculated using Eqns. (9) and (10) and the particular capillary number Nc calculated for each respective cell. The baseline curve connecting these endpoints is preferably adjusted by scaling. Scaling of the relative permeability could be done using several methods. Equation (12) shows such a method:
Krg =F( _______________________________________ ) (12) 1¨ ¨ S 0 rg Eqn. (12) simply states that krg is a function of Sg , Sg, and Sorg . (For gas saturation greater than Sorg the oil phase is immobile ¨ i.e., Kro = 0). The function F
could be (but is not limited to) a simple power law:
S ¨S
Krg K rgro ( 8 __ 8c 2 ) (13) 1¨Sg, rg In the conventional treatment of gas relative permeabilities, Sge in Eqn. 12 or 13 is equal to S; . However, with this formulation, Sgc in Eqns. 12 and 13 is now a function of the capillary number.
Additionally, if the endpoint of gas relative permeability krg, is decreased by 10%
relative to the original kr g,õ of the baseline curve, then all gas relative permeability values on the correlation or curve will be decreased by 10%. Those skilled in the art will appreciate that many other ways of adjusting the baseline curve to reflect changes
- 21 -in the updated values of endpoints Sg, and/or k ,..gro can be used and are within the scope of this invention as well.
C. Selecting Gas Relative Permeabilities kte For Incorporation Into The Reservoir Simulator Saturation values Sg may come from initial conditions when the reservoir simulation is first started, from the previous time step, or else from values calculated during iterations within a time step. The saturation Sg of each reservoir cell is then examined and the corresponding relative permeability k,. is selected from the adjusted baseline correlation. As described above, if Sg Sg, , then the correlation from the previously calculated curve is used to determine kg.
=
D. Running Reservoir Simulation Using Selected Gas Rrelative Permeabilities krg Finite difference equations are solved to determine unknowns, such as pressure P or saturation 5g . These finite difference equations rely upon the latest updated relative permeabilities kr, including the capillary number dependent gas permeabilites krg for =
the reservoir cells. Such finite difference equations are well known those skilled in the art of reservoir simulation. Examples of well known solution methods for such equations include: (1) Fully Explicit; (2) Implicit Pressure, Explicit Saturation (IMPES); (3) fully Implicit; (4) Sequential Implicit (SEQ), Adaptive Implicit (AIM);
and Cascade. In the preferred embodiment, a fully Implicit method is used to solve these equations.
If the solutions to a state variable, i.e. pressure or saturation, are within a satisfactory tolerance range during an iteration, then final fluid properties will be established for a timestep. Volumes of production of gas, water and oil during the tirnestep can be established from these fluid properties, as is conventionally done with reservoir simulators. The reservoir simulator may then run over many more timesteps until a predetermined length of time is met. The cumulative production over these tirnesteps provides an estimation of the production from the subterranean formation.
-22 -The present invention also include a system for carrying out the above reservoir simulation using relative permeabilities krg that are dependent upon depletion rate/fluid velocity and viscosities of crude oil. Further, the present invention also includes a program storage device which carries instructions for carrying out this reservoir simulation using fluid velocity dependent relative permeabilities.
While in the foregoing specification this invention has been described in relation to certain preferred embodiments thereof, and many details have been set forth for purpose of illustration, it will be apparent to those skilled in the art that the invention is susceptible to alteration and that certain other details described herein can vary considerably without departing from the basic principles of the invention.
-23 -=
Nomenclature a = coefficient for calculating Sg, ;
= exponent for calculating Sgc ;
Boi = oil formation volume factor at Pi;
Bo = oil formation volume factor at P;
Bg = gas formation volume factor at P;
= coefficient for calculating krgro;
Cf = rock or core sample compressibility (1/psi);
co = oil sample compressibility (1/psi);
Crthy_core = CT number of a sample saturated with gas;
CTsaturated core= CT number of a sample saturated with kerosene (at initial pressure);
CTp = CT number measured during depletion at pressure P;
CTiiq = CT number for kerosene;
CTgos = CT number for air;
= exponent for calculating k rgõ;
= change in depth from a datum;
= gravitational constant;
ke = effective permeability;
kr = relative permeability, dimensionless;
krg gas relative permeability, dimensionless;
krg, = gas relative permeability with minimum residual oil; dimensionless;
kr gro = endpoint gas relative pemeability with minimum residual oil, dimensionless;
kro= oil relative permeability, dimensionless;
K = rock permeability;
= slope of the solution-gas curve, psi-I
N = capillary number calculated for a particular cell of a reservoir model;
N0 = capillary number;
AP = change in pressure (psi);
L = length of test chamber (inches);
= oil in place (stb) at initial conditions;
= cumulative oil produced (stb) at pressure P (cm3);
(1) 0 = oil-phase potential, Pi = pressure at time i, psi;
d P, the change in pressure across a face, peff = effective density;
= density of gas;
Pg = density of oil;
R, - Gas-to-oil ratio;
S = saturation, dimensionless;
Sg= gas saturation, dimensionless;
Sge = critical gas saturation, dimensionless;
= oil saturation, dimensionless;
- 24 -Sg ,, = endpoint critical gas saturation, dimensionless;
Sam = residual oil saturation to gas for a particular rock region, dimensionless;
stb = stock tank barrel;
= interfacial tension;
aog = oil-gas interfacial tension; and vo = velocity of oil.
- 25 -

Claims (20)

WHAT IS CLAIMED IS:
1. A method of predicting a property of at least one fluid in a subterranean reservoir containing heavy oil entrained with gas, the method comprising:
(a) utilizing a baseline correlation of gas relative permeability k rg versus gas saturation S g in a reservoir simulator;
(b) utilizing a capillary number dependent correlation between at least one of critical gas saturations S gc and capillary numbers N ca and endpoint gas relative permeabilities k rgro and capillary numbers N ca in the reservoir simulator;
(c) calculating capillary numbers N c for a plurality of cells in a reservoir model representative of the subterranean reservoir;
(d) adjusting the baseline correlation to comport with at least one of S gc and k rgro selected from the capillary number dependent correlation of step (b) using the capillary number N c calculated in step (c) to produce a plurality of corresponding adjusted baseline correlations;
(e) selecting relative permeabilities k rg for the plurality of cells from the corresponding adjusted baseline correlations of step (d); and (f) running a reservoir simulation utilizing the selected relative permeabilities of step (e) to predict a property of at least one fluid in the subterranean reservoir containing heavy oil entrained with gas.
2. The method of claim 1 wherein:
the step of utilizing the capillary number dependent correlation includes utilizing a correlation between S gc and N ca.
3. The method of claim 1 wherein:
the step of utilizing the capillary number dependent correlation includes utilizing a correlation between k rgro and N ca.
4. The method of claim 1 wherein:
the step of utilizing the capillary number dependent correlation includes utilzing a correlation between S gc and N ca and utilizing a correlation between k rgro and N ca.
5. The method of claim 4 wherein:
the step of adjusting the baseline correlations comports with values of S gc selected from the correlation between S gc and N ca and with values of k rgro selected from the correlation between k rgro and N ca.
6. The method of claim 1 wherein:
the step of utilizing the capillary number dependent correlation includes conducting depletion experiments on core samples from the subterranean reservoir.
7. The method of claim 1 wherein:
the step of utilizing the capillary number dependent correlation includes predicting the correlation without conducting depletion experiments on a core sample from the subterranean reservoir.
8. The method of claim 1 wherein:
the step of developing the capillary number dependent correlation includes conducting depletion experiments on sandpack samples.
9. The method of claim 1 wherein:
the capillary number N ca which is calculated for the plurality of cells remains fixed in a time step iteration.
10. The method of claim 1 wherein:
the capillary numbers which are calculated for the plurality of cells are updated during iterations of a time step conducted in the reservoir simulation.
11. The method of claim 1 wherein:
the capillary number dependent correlation is a look-up table.
12. The method of claim 1 wherein:
the capillary number dependent correlation is a mathematical function.
13. The method of claim 1 wherein:
the reservoir simulator uses a fully implicit method to solve equations.
14. A method for simulating the flow of heavy oil in a subterranean reservoir, the method comprising:
creating a reservoir model representative of a subterranean reservoir for which fluid flow is to be simulated, the reservoir model including a plurality of reservoir cells;
determining velocity-dependent relative permeabilities for the reservoir cells that account for a velocity of the fluid flow through the reservoir cells, the velocity-dependent relative permeabilities for the reservoir cells being determined using a gas relative permeability versus gas saturation correlation constructed by scaling a baseline gas relative permeability versus gas saturation correlation to at least one of a critical gas saturation and an endpoint gas relative permeability; and running a reservoir simulation utilizing the velocity-dependent relative permeabilities to simulate the flow of heavy oil in the subterranean reservoir.
15. A system for simulating the flow of heavy oil in a subterranean reservoir, the system comprising:
a program storage device carrying computer instructions to perform a method of reservoir simulation; and a reservoir simulator for executing the computer instructions of the program storage device to perform the method of reservoir simulation, the method of reservoir simulation comprising the steps of:
(a) utilizing a baseline correlation of gas relative permeability k rg versus gas saturation S g;
(b) utilizing a capillary number dependent correlation between at least one of critical gas saturations S gc and capillary numbers N ca and endpoint gas relative permeabilities k rgro and capillary numbers N ca;
(c) calculating capillary numbers N c, for a plurality of cells in a reservoir model representative of the subterranean reservoir containing heavy oil;
(d) adjusting the baseline correlation to comport with at least one of S gc and k rgro selected from the capillary number dependent correlation of step (b) using the capillary number N c calculated in step (c) to produce a plurality of corresponding adjusted baseline correlations;
(e) selecting relative permeabilities k rg for the plurality of cells from the corresponding adjusted baseline correlations of step (d); and (f) simulating the flow of the heavy oil in the subterranean reservoir utilizing the selected relative permeabilities k rg of step (e).
16. A method for predicting a property of at least one fluid in a subterranean reservoir containing heavy oil entrained with gas, the method comprising:
(a) inputting a baseline correlation of gas relative permeability k rg versus gas saturation S g that defines a baseline critical gas saturation S gc and a baseline endpoint gas relative permeability k rgro into a reservoir simulator;
(b) calculating at least one of an updated critical gas saturation S gc and an updated endpoint gas relative permeability k rgro based upon a plurality of depletion rates;
(c) adjusting the baseline correlation to comport with at least one of the updated critical gas saturation S gc and the updated endpoint gas relative permeability k rgro;
(d) selecting relative permeabilities k rg for the plurality of cells from the adjusted baseline correlation of step (c); and (e) running a reservoir simulation utilizing the selected relative permeabilities of step (d) to predict a property of at least one fluid in the subterranean reservoir containing heavy oil entrained with gas.
17. The method of claim 14 wherein:
the baseline gas relative permeability versus gas saturation correlation defines a baseline critical gas saturation; and scaling the baseline gas relative permeability versus gas saturation correlation to at least one of the critical gas saturation and the endpoint gas relative permeability comprises adjusting all gas relative permeability values on the baseline gas relative permeability versus gas saturation correlation by a ratio of the critical gas saturation to the baseline critical gas saturation.
18. The method of claim 14 wherein:
the baseline gas relative permeability versus gas saturation correlation defines a baseline endpoint gas relative permeability; and scaling the baseline gas relative permeability versus gas saturation correlation to at least one of the critical gas saturation and the endpoint gas relative permeability comprises adjusting all gas relative permeability values on the baseline gas relative permeability versus gas saturation correlation by a ratio of the endpoint gas relative permeability to the baseline endpoint gas relative permeability.
19. The method of claim 14 wherein the critical gas saturation is calculated using a capillary number dependent correlation.
20. The method of claim 14 wherein the endpoint gas relative permeability is calculated using a capillary number dependent correlation.
CA2634757A 2005-12-22 2006-12-20 Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive Active CA2634757C (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US75390005P 2005-12-22 2005-12-22
US60/753,900 2005-12-22
PCT/US2006/049095 WO2007076044A2 (en) 2005-12-22 2006-12-20 Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive

Publications (2)

Publication Number Publication Date
CA2634757A1 CA2634757A1 (en) 2007-07-05
CA2634757C true CA2634757C (en) 2014-12-16

Family

ID=38218675

Family Applications (1)

Application Number Title Priority Date Filing Date
CA2634757A Active CA2634757C (en) 2005-12-22 2006-12-20 Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive

Country Status (4)

Country Link
CN (1) CN101366041B (en)
BR (1) BRPI0620170A2 (en)
CA (1) CA2634757C (en)
WO (1) WO2007076044A2 (en)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2009238481B2 (en) 2008-04-22 2014-01-30 Exxonmobil Upstream Research Company Functional-based knowledge analysis in a 2D and 3D visual environment
US8892407B2 (en) 2008-10-01 2014-11-18 Exxonmobil Upstream Research Company Robust well trajectory planning
WO2010053618A1 (en) 2008-11-06 2010-05-14 Exxonmobil Upstream Research Company System and method for planning a drilling operation
EP2457187B1 (en) 2009-09-02 2020-06-03 Landmark Graphics Corporation System and method of hydrocarbon formation modeling
WO2011096964A1 (en) 2010-02-03 2011-08-11 Exxonmobil Upstream Research Company Method for using dynamic target region for well path/drill center optimization
US8452580B2 (en) * 2010-02-26 2013-05-28 Chevron U.S.A. Inc. Method and system for using multiple-point statistics simulation to model reservoir property trends
WO2011112221A1 (en) 2010-03-12 2011-09-15 Exxonmobil Upstream Research Company Dynamic grouping of domain objects via smart groups
WO2013169429A1 (en) 2012-05-08 2013-11-14 Exxonmobile Upstream Research Company Canvas control for 3d data volume processing
CN103590807B (en) * 2012-08-13 2016-06-08 中国石油天然气集团公司 A kind of method determining heavy oil thermal recovery steam chamber state
WO2014200685A2 (en) 2013-06-10 2014-12-18 Exxonmobil Upstream Research Company Interactively planning a well site
CN104481523B (en) * 2014-11-11 2018-07-13 中国石油天然气股份有限公司 Simulation experiment system and method for thick oil dissolved gas drive development
CN104481504B (en) * 2014-11-11 2017-09-01 中国石油天然气股份有限公司 Simulation experiment system and method for thick oil dissolved gas drive development
FR3034894B1 (en) * 2015-04-09 2018-08-10 Services Petroliers Schlumberger OIL STORAGE SATURATION TANK AND PERMEABILITY MODELING
FR3036820B1 (en) 2015-06-01 2021-12-31 Services Petroliers Schlumberger MODELING OF OIL FIELD RESERVOIR SATURATION AND PERMEABILITY
FR3038408B1 (en) 2015-06-30 2017-08-25 Services Petroliers Schlumrberger MODELING THE SATURATION AND PERMEABILITY OF PETROLEUM FIELD RESERVOIR
CN105443080A (en) * 2015-12-04 2016-03-30 中国石油天然气股份有限公司 Simulation device for failure development of bottom-edge water gas-containing oil reservoir
CN107101773B (en) * 2017-04-14 2019-07-09 中国石油天然气股份有限公司 Method and device for testing bubble point pressure value and quasi-bubble point pressure value
CN107944599B (en) * 2017-10-31 2020-10-09 中国石油天然气股份有限公司 Method for predicting yield of oil-gas horizontal well
CN109281662B (en) * 2018-08-16 2022-03-29 中国石油天然气股份有限公司 Method for determining relative compressibility index of reservoir at different section positions of horizontal well
US11492895B2 (en) 2018-11-13 2022-11-08 Saudi Arabian Oil Company Relative permeability ratio from wellbore drilling data
CN110671084B (en) * 2019-09-04 2021-07-06 中国石油化工股份有限公司 Method for simulating flooding numerical value of common thick oil water-soluble viscosity reducer
CN110863806B (en) * 2019-11-28 2021-07-23 西安石油大学 Carbon dioxide flooding gas front dynamic change prediction method
WO2024144813A1 (en) * 2022-12-30 2024-07-04 Landmark Graphics Corporation Adjusting fluid injection into a wellbore based on relative permeability of a formation

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6052520A (en) * 1998-02-10 2000-04-18 Exxon Production Research Company Process for predicting behavior of a subterranean formation
US6227296B1 (en) * 1998-11-03 2001-05-08 Exxonmobil Upstream Research Company Method to reduce water saturation in near-well region
FR2837572B1 (en) * 2002-03-20 2004-05-28 Inst Francais Du Petrole METHOD FOR MODELING HYDROCARBON PRODUCTION FROM A SUBTERRANEAN DEPOSITION SUBJECT TO DEPLETION
US6810332B2 (en) * 2003-01-31 2004-10-26 Chevron U.S.A. Inc. Method for computing complexity, confidence and technical maturity indices for reservoir evaluations
CN1302386C (en) * 2003-10-17 2007-02-28 大庆油田有限责任公司 Three-dimensinal composite drive computer simulating method for low-concentration surface active agent and phase state joint

Also Published As

Publication number Publication date
WO2007076044A2 (en) 2007-07-05
BRPI0620170A2 (en) 2011-11-01
CN101366041B (en) 2011-10-12
CN101366041A (en) 2009-02-11
CA2634757A1 (en) 2007-07-05
WO2007076044A3 (en) 2007-12-27

Similar Documents

Publication Publication Date Title
CA2634757C (en) Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive
US9020793B2 (en) Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive
US7289942B2 (en) Performance prediction method for hydrocarbon recovery processes
Wu et al. Flow and displacement of Bingham non-Newtonian fluids in porous media
Wu et al. Analysis of multiphase non-Darcy flow in porous media
CN114372352B (en) Method for predicting peak regulation capacity of gas storage of complex fault block oil reservoir through seepage-temperature double-field coupling numerical simulation
Behmanesh et al. Treatment of rate-transient analysis during boundary-dominated flow
Hussain et al. A semi-analytical model for two phase immiscible flow in porous media honouring capillary pressure
Baker Streamline technology: reservoir history matching and forecasting= its success, limitations, and future
Kumar et al. Effect of depletion rate on gas mobility and solution gas drive in heavy oil
Hossain et al. Numerical approach for the prediction of formation and hydraulic fracture properties considering elliptical flow regime in tight gas reservoirs
Fattah et al. New Inflow Performance Relationship for solution-gas drive oil reservoirs
Ito The introduction of the microchanneling phenomenon to cyclic steam stimulation and its application to the numerical simulator (sand deformation concept)
Wu Non‐Darcy displacement of immiscible fluids in porous media
He et al. CO2 storage capacity estimation by considering CO2 Dissolution: A case study in a depleted gas Reservoir, China
Sui et al. Determining multilayer formation properties from transient temperature and pressure measurements in gas wells with commingled zones
Beraldo et al. Compressible streamline-based simulation with changes in oil composition
Takbiri-Borujeni et al. Upscaling the steam-assisted-gravity-drainage model for heterogeneous reservoirs
Yang et al. A Practical Approach to History Matching Premature Water Breakthrough in Waterflood Reservoir Simulation: Method and Case Studies in South Belridge Diatomite Waterflood
Ali Developments in measurement and interpretation techniques in coreflood tests to determine relative permeabilities
Jones Jr et al. A predictive model for water and polymer flooding
Al-Mossawy A simple approach to identify the proper relative permeability model
Worawutthichanyakul Unstable immiscible displacement study in oil-wet rocks
Camacho-V Constant-pressure production in solution-gas-drive reservoirs: transient flow
Azim et al. 3D numerical and experimental study on upscaling two phase relative permeability curve of naturally fractured reservoirs

Legal Events

Date Code Title Description
EEER Examination request