Reservoir Simulation of Variable Bubble-Point Problems

1976 ◽  
Vol 16 (01) ◽  
pp. 10-16 ◽  
Author(s):  
L.K. Thomas ◽  
W.B. Lumpkin ◽  
G.M. Reheis

Abstract This paper presents the development of a general beta reservoir simulator that will model conventional (fixed bubble-point) problems as well as problems involving a variable bubble point, such as gas injection projects above the original bubble point and water injection projects resulting in a collapsing gas saturation, Provisions are included for allowing the pressure to cross the bubble point with the same relative ease as in a conventional simulator. Example problems are presented to demonstrate the utility of the model for gas and water injection problems. problems Introduction Many reservoir simulation problems involve treating a variable bubble point throughout the reservoir. For example, when gas is injected into an undersaturated reservoir, gas will go into solution, increasing the bubble point of the oil. As this oil moves away from the injector, the bubble point of surrounding areas also may increase point of surrounding areas also may increase because of mixing. During waterfloods of saturated reservoirs, the gas saturations in regions near the injectors frequently reduce to zero at pressures below the original bubble point. Thus, the bubble point will vary areally throughout the field. point will vary areally throughout the field.Recent publications have discussed certain aspects of the variable bubble-point problem. Most of these papers contain only a brief discussion of this problem. Ridings discusses the need for allowing saturation pressure to vary continuously throughout the reservoir as long as there is free gas associated with the oil. In the model presented by Cook et al., free gas saturation is monitored for saturated systems and the bubble point is set equal to the prevailing reservoir pressure when the gas saturation in a cell disappears. Bubble points for undersaturated cells are allowed to change because of the entrance of free gas and mixing. Spilletta et al. assume that a cell that is saturated or undersaturated at the beginning of a time step will remain so throughout the time step. They then solve their saturation equations for water and gas saturations. The bubble point of any undersaturated cell is adjusted to account for nonzero gas saturation, and the water saturation is modified to conserve oil. Steffensen and Sheffield devote their paper to the reservoir simulation of a collapsing gas saturation during waterflooding. In their model, blocks that have a free gas saturation at the beginning of a time step and have zero or negative gas saturations at the end of a time step are detected, and the bubble points for these cells are set equal to the estimated pressure where Sg reduced to zero. Gas saturation for these blocks is set equal to zero and S is set to 1 - S . The oil saturation in adjacent saturated grid blocks is then adjusted so that oil material balance is maintained. Mixing caused by flow between undersaturated blocks is neglected. This paper presents a comprehensive analysis of modeling variable bubble-point problems.* It treats the specific problems of simulating gas injection above the bubble point as well as waterflooding depleted reservoirs. It differs from previously reported work by accounting for the effect of bubble-point change on computed pressure change during a time step. Also, provisions are included that allow the pressure to cross the bubble point with the same relative ease as in a conventions! simulator In regard to waterflooding, mis paper differs from the work of Steffensen and Sheffield in mat it allows for bubble-point changes caused by mixing. DEVELOPMENT OF FLOW EQUATIONS To simulate the variable bubble-point problem, the expansion of the flow equations above the bubble point must include the effects of pressure and bubble point on fluid properties. Also, special consideration must be given to cells passing through the bubble point if both pressure and material-balance errors are to be eliminated. SPEJ P. 10

2005 ◽  
Vol 53 (6) ◽  
pp. 803-810 ◽  
Author(s):  
José M. Carcione ◽  
Davide Gei ◽  
Giuliana Rossi ◽  
Gianni Madrussani

2021 ◽  
Author(s):  
Chris Boeije ◽  
Pacelli Zitha ◽  
Anne Pluymakers

<p>Geothermal energy, the extraction of hot water from the subsurface (500 m to 5 km deep), is generally considered one of the key technologies to achieve the demands of the energy transition.  One of the main problems during production of geothermal waters is degassing. Many subsurface waters contain substantial amounts of dissolved gasses. As the hot water travels up the production well, the pressure and/or temperature drop will cause dissolved gas to come out of the solution. This causes several problems, such as corrosion of the facilities (due to pH changes and/or degassing-related precipitation) and in some cases even to blocking of the reservoir as the free gas limits the water flow.  To better understand under which conditions free gas nucleates, we need confirmation of theoretical bubble point pressure and temperature, and understand what controls the evolution of the bubble front:  i.e. what are the conditions under which free gas emerges from the solution and at what rate are bubbles created?</p><p>An experimental setup was designed in which the degassing process can be observed visually. The setup consists of a high-pressure visual cell which contains water saturated with dissolved gas at high-pressure. The pressure within the cell can be reduced in a reproducible manner using a back-pressure regulator at the outlet of the system. A high-speed camera paired with a uniform LED light source is used to record the degassing process. The pressure in the cell is monitored using a pressure transducer which is synchronized with the camera. The resulting images are then analysed using a MATLAB routine, which allows for determination of the bubble point pressure and rate of bubble formation.</p><p>The first two sets of experiments at ambient temperatures (~20 <sup>o</sup>C) were carried out using two different gases, N<sub>2</sub> and CO<sub>2</sub>. Initial pressure was 70 and 30 bar for the N<sub>2</sub> and CO<sub>2</sub> experiments respectively. In these first experiments we determined the influence of the initial fluid used to pressurize the system. Using gas as the initial fluid causes a large amount of bubbles, whereas only a single bubble was observed for a system where degassed water is used as the initial fluid. An intermediate system where degassed water is pumped into a system full of air at ambient conditions and is subsequently pressurized yields a number of bubbles in between the two systems described previously. All three methods give reproducible bubble point pressures within 2 bar (i.e. pressure where the first free bubble is formed). There are clear differences in bubble point between N<sub>2</sub> and CO<sub>2</sub>.</p><p>A series of follow-up experiments is planned that will investigate specific properties at more extreme conditions: at higher pressures (up to 500 bar) and temperatures (500 <sup>o</sup>C) and using high-salinity brines (2.5 M).</p>


Author(s):  
Luca Casarsa ◽  
Pietro Giannattasio ◽  
Diego Micheli

A simple and efficient numerical model is presented for the simulation of pulse combustors. It is based on the numerical solution of the quasi-1D unsteady flow equations and on phenomenological sub-models of turbulence and combustion. The gas dynamics equations are solved by using the Flux Difference Splitting (FDS) technique, a finite-volume upwind numerical scheme, and ENO reconstructions to obtain second-order accurate non-oscillatory solutions. The numerical fluxes computed at the cell interfaces are used to transport also the reacting species, their formation energy and the turbulent kinetic energy. The combustion progress in each cell is evaluated explicitly at the end of each time step according to a second-order overall reaction kinetics. In this way, the computations of gas dynamic evolution and heat release are decoupled, which makes the model particularly simple and efficient. A comprehensive set of measurements has been performed on a small Helmholtz type pulse-jet in order to validate the model. Air and fuel consumptions, wall temperatures, pressure cycles in both combustion chamber and tail-pipe, and instantaneous thrust have been recorded in different operating conditions of the device. The comparison between numerical and experimental results turns out to be satisfactory in all the working conditions of the pulse-jet. In particular, accurate predictions are obtained of the device operating frequency and of shape, amplitude and phase of the pressure waves in both combustion chamber and tail-pipe.


2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Zhiwang Yuan ◽  
Zhiping Li ◽  
Li Yang ◽  
Yingchun Zhang

When a conventional waterflooding characteristic curve (WFCC) is used to predict cumulative oil production at a certain stage, the curve depends on the predicted water cut at the predicted cutoff point, but forecasting the water cut is very difficult. For the reservoirs whose pressure is maintained by water injection, based on the water-oil phase seepage theory and the principle of material balance, the equations relating the cumulative oil production and cumulative water injection at the moderately high water cut stage and the ultrahigh water cut stage are derived and termed the Yuan-A and Yuan-B curves, respectively. And then, we theoretically analyze the causes of the prediction errors of cumulative oil production by the Yuan-A curve and give suggestions. In addition, at the ultrahigh water cut stage, the Yuan-B water cut prediction formula is established, which can predict the water cut according to the cumulative water injection and solve the difficult problem of water cut prediction. The application results show Yuan-A and Yuan-B curves are applied to forecast oil production based on cumulative water injection data obtained by the balance of injection and production, avoiding reliance on the water cut forecast and solving the problems of predicting the cumulative oil production of producers or reservoirs that have not yet shown the decline rule. Furthermore, the formulas are simple and convenient, providing certain guiding significance for the prediction of cumulative oil production and water cut for the same reservoir types.


2005 ◽  
Author(s):  
Bernhard Hustedt ◽  
Yuan Qiu ◽  
Dirk Zwarts ◽  
Paul Jacob van den Hoek

1970 ◽  
Vol 10 (03) ◽  
pp. 211-218 ◽  
Author(s):  
J.M. Dumore

Dumore, J.M., Koninklijke/Shell Exploratie En Produktie Laboratorium Rijswijk, The Netherlands Abstract Solution-gas drive is usually described in the literature as a parallel flow of oil and gas, with the gravitational factor being neglected. Under the influence of gravity, however, gas evolving from the oil will migrate upward. The gas then accumulates at the top of the formation where, if not originally present, a gas cap is formed. The gas flows to the producing wells via the gas cap. For an analysis of the drive mechanisms in many oil reservoirs, it is essential to know the magnitude of the gas saturation that develops during solution-gas drive in the oil layer below the gas cap. However, the way in which the gas saturation develops and the parameters on which this development depends are unknown. A study was begun and the results are presented here. Experiments were conducted in packs of glass grains, saturated with a liquid that produced transparence. Gas was injected slowly in the packs via a pore near the bottom. Tests show that packs via a pore near the bottom. Tests show that the upward migration of the gas depends on two conditions related to capillary pressure. At low capillary pressures a conically shaped, gas-saturated region develops, through which the gas is transported upward; whereas at high capillary pressures only one gas channel develops. These pressures only one gas channel develops. These two conditions are called dispersion and nondispersion. Both may occur in oil reservoirs. Under dispersion conditions, the process of solution-gas drive in the oil layer develops in such a way that the entire layer is eventually occupied by disconnected agglomerations of gas bubbles. This results in a-high gas saturation (approximately 20 percent). Under nondispersion conditions, a network of gas channels develops. The lower the pressure-decline rate, the larger the network spacing. pressure-decline rate, the larger the network spacing. Gas saturations of less than 2 percent are often formed. Introduction The solution-gas drive process is usually described in the literature on the basis of the following assumptions:during the process a free-gas saturation develops, such that the gas and oil phases remain in equilibrium;the phases to not segregate by gravity; andrelative permeabilities to gas and oil are unique functions permeabilities to gas and oil are unique functions of the saturations. As a consequence of these assumptions, theoretically, the oil production, hence, the mean-gas saturation developed in the reservoir are related uniquely to the reservoir pressure, and are thus independent of the rate of pressure decline. In reality, gas bubbles are created when a certain critical supersaturation has been reached. As pressure continues to decline, these bubbles grow pressure continues to decline, these bubbles grow and, under the influence of gravity, upward migration of free gas will occur. The gas then accumulates at the top of the formation where, if not originally present, a gas cap is formed. The free gas flows to present, a gas cap is formed. The free gas flows to the producing wells via the gas cap. In order to analyze the drive mechanisms active in many oil reservoirs, we need to know the magnitude of the gas saturation that develops during solution-gas drive in the oil layer below the gas cap. However, the way in which the gas saturation develops and the parameters on which its development depends are unknown. Laboratory experiments have been conducted to investigate this development and to determine whether it is influenced by the rate of pressure decline. UPWARD MIGRATION OF GAS The upward migration of gas was observed in experiments carried out in transparent models. Each lucite model contained a pack of crushed pyrex glass of a narrow sieve fraction. The pack was saturated with a kerosene-Novasol* mixture. SPEJ p. 211


Sign in / Sign up

Export Citation Format

Share Document