The design of an optimum spacing between oil wells entails both reservoir characterization and economics considerations. High hydrocarbon recovery requires short distances between wells. However, higher well density leads to a greater development cost. Accordingly, determination of an optimum well spacing is primordial in the development of oil fields. As a matter of fact, the identification of optimum well spacing for heterogeneous sandstone reservoirs undergoing waterflooding requires extensive analytical and numerical studies. The intent of this work is therefore to develop type curves as a quick tool in estimating ultimate recovery and reduce excessive reservoir simulation cost in analog reservoirs. These type curves utilize reservoir heterogeneity and well spacing in the estimating of oil recovery. In this work, we investigated numerically the effects of heterogeneity and well spacing on ultimate recovery using Eclipse black oil simulation and PEEP economic software 2015 and 2009 versions, respectively. The study involved a 50-ft thick Middle Eastern reservoir with porosity variability ranging from 0.2 to 0.9. Corresponding average matrix permeabilities of 1, 10 and 100 md were considered. Type curves relating well spacing and heterogeneity to ultimate oil recovery were developed. Type curves and net present value calculations indicated that there is exists an ultimate well spacing for each of the considered matrix permeabilities.
Proper distribution of wells in a reservoir undergoing waterflooding leads to both low cost and high sweep efficiency, which yields higher recovery. For an injector-producer pair, the sweep efficiency, controlled by heterogeneity, plays a major role in determining the optimum well placement.
Well spacing is defined as the acreage of the productive area divided by the total number of wells in the same area presented in terms of acres/well. Well spacing has been the main concern for the industry for many years. Wu, Laughlin, and Lardon [
It has also been proven that the lower the spacing between the producing wells, the higher the recovery (Bobar [
In the early days, when the concept of proper adjustment of wells for an optimum reservoir development was not yet well defined, the wells were completed with wider spacing. That was because ultimate recovery was considered to be independent of well spacing (Kern [
However, other authors have indicated that closer well spacing leads to an additional oil recovery ranging from 2% to up to 14% depending upon reservoir quality, drive mechanism and production practices. For water flooding, up to 70% of the oil-originally-in-place (OIIP) will be produced when reducing the well spacing to optimized distances (Sloan [
So, it is fair to say that maximizing oil recovery is a complicated and controversial process and does not depend on well spacing alone and that many other factors have to be considered. Many authors have concluded that increase in recovery depends on reservoir characteristics such as reservoir quality, i.e. heterogeneity, permeability, porosity, connectivity, drive mechanism, relative permeability, injectivity, productivity, and production practices (Johen, Al-Qabandi, and Anderson [
A few authors indicated that fluid properties, sweep efficiency and net pay are major factors that could affect oil recovery to some extent. They added that project economics has to also be tied up to rock and fluid properties (Tokunaga and Hise [
It was also concluded that recovery efficiency correlates well with permeability and to lesser degree with connate water saturation and also shows good fit with the size of the reservoir (Kern [
Suarez and Pichon [
So, most authors agree that oil recovery depends on reservoir heterogeneity and well spacing. The following methodology is an illustration of heterogeneity characterization and well spacing definition.
In order to study the effects of heterogeneity on oil recovery, permeability was generated from porosity measurements (Figures 1-3). Neutron porosity logs were ran in many wells and porosities were recorded at different depths. The porosity data bank was converted to permeabilities using the following formula:
k = a 1 0 f b (1)
where “a” and “b” are correlation coefficients; 0.005 and 25, respectively, k is permeability and f is porosity; in md and fractions, respectively. Three average
permeability ranges (classes) of 1, 10 and 100 md have been identified (Figures 1-3).
To characterize heterogeneity, Dykstra-Parsons [
V = k ¯ − k σ k ¯ (2)
Permeability, k, was plotted on a log-probability paper (see
To discretize heterogeneity and for dynamic modelling purposes, permeability maps have been generated for the 3 different permeability classes (1, 10, and 100 md) and exported as input to the developed Eclipse data files.
To study the effect of well spacing on oil recovery, data files using 15 × 40 × 50 grid block system were developed. The injector and producer wells were arranged according to a line-drive process (see
The distance between the injector and producing lines was kept at 820 ft (250 m). Two more wells were added longitudinally each time by infill drilling, reducing well spacing between like wells (injector to injector and producer to producer). At the end, 28 wells were input into the model for a total number of 168 simulation runs.
To study the effects of heterogeneity and well spacing on oil recovery, a sensitivity analysis involving different permeability variability coefficients V (0.0, 0.2, 0.5, 0.7 and 0.9) for average matrix permeabilities of 1, 10 and 100 md at
different well spacing was done.
Oil recovery following waterflooding was simulated for 30 years for all possible combinations (Figures 8-10). Increased variability led to lower recoveries
Well Spacing | Permeability Variability | ||||
---|---|---|---|---|---|
(acres/well) | |||||
185 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
93 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
62 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
46 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
37 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
31 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
26 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
23 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
21 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
19 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
17 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
15 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
14 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
13 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
and this is consistent with Babadagli findings [
The same can be observed for a matrix permeability in the neighborhood of 10 mD. Recovery numbers are, however, higher due to a higher productivity and a stable waterflood displacement front. For this case, recovery ranged from a high of 66.6% at a spacing of 31 cares/well and a 0 variability to a low of 29.9% at a spacing of 185 acres/well and a variability of 0.9.
The same cannot be said for a matrix permeability in the vicinity of 100 mD. Recovery observed a peak of 63.9% at a spacing of 31 acres/well and a 0 variability and a low of 32.6% for a larger spacing of 185 acres/well and a variability of 0.2. This can be triggered by fingering in the high permeability streaks and an unstable water displacement front.
Results were verified with the following series of well logs (see
The viewgraphs (1) to (9) of a post-waterflood resistivity log ran on one of the wells indicate that water displacement (purple) is not uniform due to the high degree of heterogeneity (V = 0.7). The graphs also show that the oil bank (light blue) is bypassed in layers with higher permeability and as a result, recovery is low.
Besides well spacing and the degree of reservoir heterogeneity, oil recovery optimization involves the study of economics. Net present value (NPV) was utilized in the feasibility study of each of the tested scenarios. NPV is defined as the difference between the present value of a future net income and the present value of the total capital expenditure. This parameter was obtained using PEEP.
Figures 12-14 show the resulting relationship between NPV, permeability variability and well spacing. In the figures, only variabilities of 0.2, 0.5, 0.7 and 0.9 were considered. The figures indicate that there exists an optimum NPV for each permeability variability and well spacing.
It is important to say that spacing decision should be based on NPV calculations. It was found that for a matrix permeability in the range of 1 mD, highest NPVs are achieved at a spacing of 37 acres/well for permeability variabilities of 0.2, 0.5, 0.7, and 0.9 (
For matrix permeability in the vicinity of 100 mD,
Well Spacing | Net Present Value | ||||
---|---|---|---|---|---|
(acres/well) | |||||
185 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
93 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
62 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
46 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
37 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
31 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
26 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
23 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
21 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
19 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
17 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
15 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
14 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
13 | V = 0.0 | V = 0.2 | V = 0.5 | V = 0.7 | V = 0.9 |
Type curves have been developed for sandstone reservoirs undergoing waterflooding. The curves apply for analog reservoirs with a porosity range of 0.2 to
0.9 and corresponding permeability classes of 1, 10, and 100 md. A field-scale waterflooding performance model for a Middle Eastern heterogeneous reservoir was used as a base case. Reservoir simulation results indicated that recovery efficiency decreases as permeability variability increases. Based on NPV calculations, it was also concluded that there exists an optimum well spacing for all tested permeability variabilities. Well spacings of 37, 31, and 17 acres/well for permeability matrix values of 1, 10 and 100 mD maximized the field NPV.
Trabelsi, R., Boukadi, F., Lee, J., Boukadi, B., Seibi, A. and Trabelsi, H. (2017) Type Curves Relating Well Spacing and Heterogeneity to Oil Recovery in a Water Flooded Reservoir―A Case Study. Natural Resources, 8, 632-645. https://doi.org/10.4236/nr.2017.810040