A Method for Handling Spatially Varying Fluid Properties in a Simulation Model for a Fissured Reservoir

1970 ◽  
Vol 10 (01) ◽  
pp. 25-32
Author(s):  
W.P. Sibley

Abstract The high relief, fractured carbonate reservoirs of the Asmari formation in Iran have extremely thick oil columns, with a large vertical change of reservoir temperature. This large change results in a significant effect on reservoir fluid properties. In contrast, there is a remarkable uniformity of PVT properties on a horizontal plane. Therefore, to properties on a horizontal plane. Therefore, to obtain meaningful results from reservoir engineering studies, PVT properties must be carefully weighted vertically. Intensive study of the known oil recovery mechanisms within these highly fissured systems resulted in a sophisticated reservoir simulation model. The model is programmed to include these recovery processes, which occur essentially in horizontal layers or zones. It also includes a technique for volumetrically weighting the large vertical variation of the PVT data. A description of this weighting process is the primary purpose of this paper. Introduction Oilfield structures found in Iran have been described as very long, high relief, assymmetrical anticlines that contain unusually thick oil columns (see Fig. 1). The oil reservoirs generally are capped with large primary gas accumulations and are often affected by natural water drive. Producing formations include the Asmari, Bangestan, and Khami carbonates, with the Asmari being by far the most common and prolific. Although some Asmari reservoirs have been discovered that contain sandstones interbedded within the limestones and dolomites, most reservoirs in Iran contain the bulk of the oil in compact carbonates that have been contorted and highly fractured during structural deformation. The resulting anticlinal oil accumulations are produced mainly from complex fracture systems, which apparently exist quite uniformly throughout the dense matrix. Study of surface rocks, cores and well producibility show that well developed fissure producibility show that well developed fissure systems are responsible for the excellent fluid communication. Because the fissures contain relatively little oil, maintenance of such prolific rates is dependent upon rate of oil replenishment from the adjacent matrix. The various productive mechanisms determine the rate of replenishment and duration. RESERVOIR AND MODEL CONSIDERATIONS Engineering studies of the fissured Iranian reservoirs have led to a mathematical model that is even more sophisticated than an earlier one. Many complex features have been included in the current model, making it a useful aid for advising management. The model consists of a master digital computer program that encompasses 54 subprograms and runs on any of the IBM 7040, 7090 and 360/65 computers. The known oil recovery mechanisms included in the model are based on consideration of practical oil recovery observations. For example, oil recovery from the expanding gas cap is predicted by using recognized gravity drainage formations, but if necessary, retrograde condensation recovery from the gas cap can also be included. Fluid and rock expansion, solution gas drive, and water displacement are also included, and vary according to the reservoir in question. But the manner in which gas production differs from what is usually observed in a solution gas drive must be known to include a proper simulation of behavior. Experience has shown that the producing gas-oil ratio continues to decline during the entire producing life of a well unless the proximity of the gas cap results in gas coning, at which time the well is shut in. Material balance calculations show that the initial gas-oil separation singe actually occurs in the reservoir, with the fissure system readily enabling the liberated gas to segregate toward the gas cap. SPEJ P. 25

SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 170-179 ◽  
Author(s):  
Songyan Li ◽  
Zhaomin Li

Summary Foamy-oil flow has been successfully demonstrated in laboratory experiments and site applications. On the basis of solution-gas-drive experiments with Orinoco belt heavy oil, the effects of temperature on foamy-oil recovery and gas/oil relative permeability were investigated. Oil-recovery efficiency increases and then decreases with temperature and attains a maximum value of 20.23% at 100°C. The Johnson-Bossler-Naumann (JBN) method has been proposed to interpret relative permeability characteristics from solution-gas-drive experiments with Orinoco belt heavy oil, neglecting the effect of capillary pressure. The gas relative permeability is lower than the oil relative permeability by two to four orders of magnitude. No intersection was identified on the oil and gas relative permeability curves. Because of an increase in temperature, the oil relative permeability changes slightly, and the gas relative permeability increases. Thermal recovery at an intermediate temperature is suitable for foamy oil, whereas a significantly higher temperature can reduce foamy behavior, which appears to counteract the positive effect of viscosity reduction. The main reason for the flow characteristics of foamy oil in porous media is the low gas mobility caused by the oil components and the high viscosity. High resin and asphaltene concentrations and the high viscosity of Orinoco belt heavy oil increase the stability of bubble films and prevent gas breakthrough in the oil phase, which forms a continuous gas, compared with the solution-gas drive of light oil. The increase in the gas relative permeability with temperature is caused by higher interfacial tensions and the bubble-coalescence rate at high temperatures. The experimental results can provide theoretical support for foamy-oil production.


1961 ◽  
Vol 1 (03) ◽  
pp. 142-152 ◽  
Author(s):  
J.S. Levine ◽  
M. Prats

Abstract Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semi-steady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-terminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semi-steady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curve was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative permeability properties of the reservoir are fixed, the producing gas-oil ratio is found to be a function of the fraction of oil recovered.


2021 ◽  
Author(s):  
Hilario Martin Rodriguez ◽  
Yalda Barzin ◽  
Gregory James Walker ◽  
Markus Gruenwalder ◽  
Matias Fernandez-Badessich ◽  
...  

Abstract This study has double objectives: investigation of the main recovery mechanisms affecting the performance of the gas huff-n-puff (GHnP) process in a shale oil reservoir, and application of optimization techniques to modelling of the cyclic gas injection. A dual-permeability reservoir simulation model has been built to reproduce the performance of a single hydraulic fracture. The hydraulic fracture has the average geometry and properties of the well under analysis. A history match workflow has been run to obtain a simulation model fully representative of the studied well. An optimization workflow has been run to maximize the cumulative oil obtained during the GHnP process. The operational variables optimized are: duration of gas injection, soaking, and production, onset time of GHnP, injection gas flow rate, and number of cycles. This optimization workflow is launched twice using two different compositions for the injection gas: rich gas and pure methane. Additionally, the optimum case obtained previously with rich gas is simulated with a higher minimum bottom hole pressure (BHP) for both primary production and GHnP process. Moreover, some properties that could potentially explain the different recovery mechanisms were tracked and analyzed. Three different porosity systems have been considered in the model: fractures, matrix in the stimulated reservoir volume (SRV), and matrix in the non-SRV zone (virgin matrix). Each one with a different pressure profile, and thus with its corresponding recovery mechanisms, identified as below: Vaporization/Condensation (two-phase system) in the fractures.Miscibility (liquid single-phase) in the non-SRV matrix.Miscibility and/or Vaporization/Condensation in the SRV matrix: depending on the injection gas composition and the pressure profile along the SRV the mechanism may be clearly one of them or even both. Results of this simulation study suggest that for the optimized cases, incremental oil recovery is 24% when the gas injected is a rich gas, but it is only 2.4% when the gas injected is pure methane. A higher incremental oil recovery of 49% is obtained, when injecting rich gas and increasing the minimum BHP of the puff cycle above the saturation pressure. Injection of gas results in reduction of oil molecular weight, oil density and oil viscosity in the matrix, i.e., the oil gets lighter. This net decrease is more pronounced in the SRV than in the non-SRV region. The incremental oil recovery observed in the GHnP process is due to the mobilization of heavy components (not present in the injection gas composition) that otherwise would remain inside the reservoir. Due to the main characteristic of the shale reservoirs (nano-Darcy permeability), GHnP is not a displacement process. A key factor in success of the GHnP process is to improve the contact of the injected gas and the reservoir oil to increase the mixing and mass transfer. This study includes a review of different mechanisms, and specifically tracks the evolution of the properties that explain and justify the different identified mechanisms.


1998 ◽  
Vol 1 (05) ◽  
pp. 416-420 ◽  
Author(s):  
G.E. Petrosky ◽  
F. Farshad

This paper (SPE 51395) was revised for publication from paper SPE 26644, first presented at the 1993 SPE Annual Technical Conference and Exhibition, Houston, 3-6 October. Original manuscript received for review 25 October 1993. Revised manuscript received 1 October 1997. Paper peer approved 28 January 1998. Summary New empirical pressure-volume-temperature (PVT) correlations for estimating bubblepoint pressure, solution gas-oil ratio (GOR), bubblepoint oil formation volume factor (FVF), and undersaturated isothermal oil compressibility have been developed as a function of commonly available field data. Results show that these PVT properties can be predicted with average absolute errors ranging from 0.64% for bubblepoint oil FVF to 6.66% for undersaturated isothermal oil compressibility. P. 416


SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1236-1253 ◽  
Author(s):  
Tae Wook Kim ◽  
E.. Vittoratos ◽  
A. R. Kovscek

Summary Recovery processes with a voidage-replacement ratio (VRR) (VRR = injected volume/produced volume) of unity rely solely on viscous forces to displace oil, whereas a VRR of zero relies on solution-gas drive. Activating a solution-gas-drive mechanism in combination with waterflooding with periods of VRR less than unity (VRR < 1) may be optimal for recovery. Laboratory evidence suggests that recovery for VRR < 1 is enhanced by emulsion flow and foamy (i.e., bubbly) crude oil at pressures under bubblepoint for some crude oils. This paper investigates the effect of VRR for two crude oils referred to as A1 (88 cp and 6.2 wt% asphaltene) and A2 (600 cp and 2.5 wt% asphaltene) in a sandpack system (18-in. length and 2-in. diameter). The crude oils are characterized with viscosity, asphaltene fraction, and acid/base numbers. A high-pressure experimental sandpack system (1 darcy and Swi = 0) was used to conduct experiments with VRRs of 1.0, 0.7, and 0 for both oils. During waterflood experiments, we controlled and monitored the rate of fluid injection and production to obtain well-characterized VRR. On the basis of the production ratio of fluids, the gas/oil and /water relative permeabilities were estimated under two-phase-flow conditions. For a VRR of zero, the gas relative permeability of both oils exhibited extremely low values (10−6−10−4) caused by internal gas drive. Waterfloods with VRR < 1 displayed encouraging recovery results. In particular, the final oil recovery with VRR = 0.7 [66.2% original oil in place (OOIP)] is more than 15% greater than that with VRR = 1 (55.6% OOIP) with A1 crude oil. Recovery for A2 with VRR = 0.7 (60.5% OOIP) was identical to the sum of oil recovery for solution-gas drive (19.1% OOIP) plus waterflooding (40.1% OOIP). An in-line viewing cell permitted inspection of produced fluid morphology. For A1 and VRR = 0.7, produced oil was emulsified, and gas was dispersed as bubbles, as expected for a foamy oil. For A2 and VRR < 1, foamy oil was not clearly observed in the viewing cell. In all cases, the water cut of VRR = 1 is clearly greater than that of VRR = 0.7. Finally, three-phase relative permeability was explored on the basis of the experimentally determined two-phase oil/water and liquid/gas relative permeability curves. Well-known algorithms for three-phase relative permeability, however, did not result in good history matches to the experimental data. Numerical simulations matched the experimental recovery vs. production time acceptably after modification of the measured krg and krow relationships. A concave shape for oil relative permeability that is suggestive of emulsified oil in situ was noted for both systems. The degree of agreement with experimental data is sensitive to the details of gas (gas/oil system) and oil (oil/water system) mobility.


1987 ◽  
Vol 109 (4) ◽  
pp. 214-217 ◽  
Author(s):  
D. A. Obomanu ◽  
G. A. Okpobiri

Existing correlations for predicting solution gas oil ratio, Rs, and oil formation volume factor, Bo, gave standard deviations as high as 50 and 12 percent, respectively, for Nigerian crudes. New correlations developed using 503 Pressure-Volume-Temperature (PVT) data points from 100 Nigerian crude oil reservoirs of the Niger Delta Basin are presented. The correlations for Rs and Bo predict values from different reservoirs within 6 and 2 percent standard deviations, respectively, and will apply to crudes of specific gravity range 0.811 to 0.966. These correlations are applicable to other crudes with characteristics similar to those of Nigerian crudes.


1962 ◽  
Vol 14 (06) ◽  
pp. 595-598
Author(s):  
Thomas G. Roberts ◽  
H. Edison Ellis

Sign in / Sign up

Export Citation Format

Share Document