Next Article in Journal
Data-Driven Model for Real-Time Estimation of NOx in a Heavy-Duty Diesel Engine
Previous Article in Journal
The Contribution of Lean Management—Industry 4.0 Technologies to Improving Energy Efficiency
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Influx Impact on Oil Production in Hamzeh Oil Reservoir in Northeastern Jordan: Case Study

by
Mehaysen Al-Mahasneh
1,
Hussam Elddin Al-Khasawneh
1,*,
Kamel Al-Zboon
2,
Marwan Al-Mahasneh
1 and
Ali Aljarrah
1
1
Chemical Engineering Department, Al-Huson University College, Al-Balqa’ Applied University, Al-Salt P.O. Box 19117, Jordan
2
Water and Environmental Engineering Department, Al-Huson University College, Al-Balqa’ Applied University, Al-Salt P.O. Box 19117, Jordan
*
Author to whom correspondence should be addressed.
Energies 2023, 16(5), 2126; https://doi.org/10.3390/en16052126
Submission received: 22 December 2022 / Revised: 27 January 2023 / Accepted: 9 February 2023 / Published: 22 February 2023

Abstract

:
This paper was conducted to delimit the water influx in the Hamzeh oil reservoir, located in northeastern Jordan approximately 150 km east of Amman. Petroleum reservoirs are frequently encompassed by water aquifers that back up the reservoir pressure through water inflow. When the pressure declines in a petroleum reservoir, the water aquifer responds by providing an influx of water. Gradually, the damage is reduced and then eliminated, and more oil is produced from the reservoir. The material balance equation (MBE) is used as the fundamental method for this study, predicting reservoir performance for a period of 11 years. The results for this study prove that the reservoir has a water drive mechanism and that the original oil in place (OOIP) was 24,958,290 m3. The projected oil recovery factor ranges from 10.9 to 25 percent for the Hummar and Shueib formations, respectively, depending on the areal efficiency assumed in the calculations. The water influx for the 11-year period was predicted by an MBE, an unsteady-state model, and the results of the performance reservoir.

1. Introduction

Petroleum reservoirs are often surrounded by water aquifers that support the reservoir pressure through water influx. In response to a pressure drop at a petroleum reservoir, a water aquifer reacts to offset or retard the pressure decline by providing a water influx or encroachment [1]. This work is important because it defines the relationships between the reservoir pressure, the quantities of oil and gas produced, the oil and gas content of the reservoir, and the characteristics of the reservoir fluids. This study uses an equation that, with sufficient data, makes it possible to calculate the oil content of the reservoir and to quantitatively predict the influence of the gas–oil ratio on the decline of reservoir pressure. Active oil and reservoir energy are the essential parameters for the determination of the rate of water encroachment into a field and its total amount [2,3]. This approach to water influx calculations offers a useful and flexible method for forecasting and analyzing the performance of water drive reservoirs [4].
The data, which provides the skin factors, indicates that the wells have suffered severe wellbore damage according to the assessment of oil well producing by matrix acidizing. The acid treatment was performed on the Hummar and Shueib formations, which improved the productivity of their wells [5].
The aims of this study are as follows:
  • To limit the drive mechanisms and to promote a conformity within suitable water influx models and unknown parameter calculations;
  • To calculate the water influx using an MBE; and
  • To predict the oil production for the next 11 years.
This study was conducted at the Hamzeh oil reservoir to delimit its water influx, to predict its performance under active water drive, and to match its production and pressure decline history to the aquifer model. Water influx models simulate and predict aquifer performance. Most importantly, they predict the cumulative water influx history. When a water influx model is successfully integrated with a reservoir simulator, the net result is a model that effectively simulates the performance of water drive reservoirs [6]. Two geometries are considered: radial- and linear-flow systems. The radial-flow system assumes that the reservoir is a right cylinder and that the aquifer surrounds the reservoir.
An integrated analytical model was developed to simulate the combination of the water drive mechanism and the oil reservoir’s performance. It couples the general material balance equation with equations for the water influx, the water-invaded pore volume, and the gas-invaded pore volume [7]. The method that was chosen for optimizing and accelerating history matching is a material balance analysis. Based on simple measurements such as average pressure, daily production rate, and pressure, volume, and temperature (PVT) tests, a material balance model can be constructed for a three-phase system of a hydrocarbon reservoir [8]. A similar study was conducted to determine the water influx in the Volve oil field, which is situated in the South Viking Graben in the North Sea in block 15/9, approximately 200 km west of Stavanger at the southern end of the Norwegian sector [9].

Description of the Hamzeh Oil Reservoir

The Azraq Basin is located in northeastern Jordan and covers an area of approximately 200 km2. Hydrocarbons were discovered in the Hamzeh field, and the total oil production from this field is approximately 1.49 × 106 bbl = 230,405 m3. The main reservoir consists of the Hummar and Shueib dolomite and limestone formations as show in Figure 1. lithostratigraphic cross section of the Hamzeh field [10]. and contains six wells. Three of these wells produce oil (Well A, Well B, and Well C), and the other three wells need to be worked over.
The properties of the studied wells, such as the rock formation and fluid properties, were extracted and tested. The data are shown in Table 1.

2. Methodology

It has been confirmed that the main driving energy responsible for oil displacement in the Hamzeh reservoir is water influx or a water drive mechanism. Water influx occurs when the reservoir pressure drops due to a long interval of fluid production; therefore, the water influx rate in the general material balance equation must amount to the volumetric withdrawal rate. The OOIP is calculated using a volumetric method and a material balance equation. These methods depend on the PVT data and the analysis of the conventional core samples.
Amerada used other tools to measure the built-up pressure and quartz pressure gauges to predict future reservoir performance [5].
Next, the reserve estimates are performed. These estimates can be classified according to the pre- and post-production stages (static and dynamic methods). Static methods indicate the analogy and volumetric calculations that are used before the reservoir production process, while dynamic methods include the performance techniques that are applied after production has started and typically require production data and pressure analysis of the wells.
A.
Volumetric method
The volumetric method is used to estimate reserves in the early periods of production, before pressure is obtained and the production data is applied to the performance techniques. The hydrocarbon estimation is acquired from the reservoir rock volume and then the final hydrocarbon total can be predicted.
B.
Material Balance Equation (MBE)
The MBE has long been considered a basic method for predicting and analyzing reservoir performance. The MBE can be used to calculate hydrocarbon content (N) and water influx (We) to forecast future reservoir performance, to evaluate reservoir pressure, and to predict the ultimate hydrocarbon recovery under many kinds of drive mechanisms. The volumes of formation waters and crude oils are affected by temperature, pressure, and the quantity of gas in solution. Compressibility of the formation waters, pore volume compressibility, volume factors, predictions of oil displacement by water influx, and decline curve calculations should all be entered into the material balance calculations. These calculations were based on the oil produced by liquid expansion above the bubble point and on the oil produced by rock and fluid compressibility.

3. Calculation to Determine Oil in Place with the Material Balance Equation during Water Influx

3.1. Drive Mechanism

At pressures above the saturation pressure, the material balance equation can be written as follows:
Np   Bo We + Wp   Bw = N   Bio   Coe   Δ P
where the effective compressibility of the oil, Coe, is defined as follows:
Coe = Soi   Co + Sw + Cw + Cf Soii  
and
Co = dBo dP   Boi
Cw = dBw dp   Bwi
With the tank oil in place, one method for estimating N OOIP is to apply Equation (1) at several different successive reservoir pressures. Each calculation set for oil production assumes that the water influx (We) is equal to zero. Afterward, the calculated values of the tank oil in place (Ncalc, oil production [Np]) are plotted.
If there has been water encroachment, the Ncalc vs. Np relationship would slope upward to the right and indicate that the driving mechanisms are the water drive and solution gas drive. If there has been no water encroachment, the relationship would be horizontal, indicating a volumetric-type reservoir with a solution gas drive.
The extrapolation of either curve to the Ncalc axis would give the true N as shown below.
N = NP   BO We + WP   BW Boi   COe   Δ P  
where Coe, the coefficient oil compressibility, is defined as
Coe = Soi   Co + Swi   Cw + Cf Soi
and
Soi = ( 1 Swi )
Co = Bo Boi   Δ P
where NP = cumulative oil produce (m3), N = OOIP = original of oil in place (m3), We = cumulative water influx (reservoir condition), Wp = cumulative water production (m3), BW = water formation volume factor (m3/m3), BoI = initial oil formation volume factor (m3/m3), BO = oil formation volume factor (m3), ΔP = cumulative pressure reduction from initial pressure (bar), Co, Cw, and Cf = compressibility of oil, water, and formation, respectively, Soi = initial oil saturation percentage, Swi = initial water saturation percentage, and Ncalc = calculated values of tank oil in place (m3 × 106). The results are presented in tables for the Hummar formation using Hummar PVT data, the Hummar formation using Shueib PVT data, and the Shueib formation.
The relationship between the calculated values of tank oil in place (Ncalc) and the cumulative oil production (Np) for both the Hummar and Shueib formations is shown on Table 2 and Table 3 and in Figure 2 and Figure 3. These figures explain that the calculated OOIP increases along with the cumulative oil production for both formations. These graphs slope upward, indicating that the drive mechanisms are water drive and solution gas drive.

3.2. Pore Volume Compressibility

The results show that the Hummar and Shueib formations are dual-porosity systems. The secondary porosity in these formations is believed to be attributable to fractures. The total pore volume compressibility of such systems can be expressed as follows:
C rt = F f × C f +   F fm × C m
where Crt = total pore volume compressibility (Pa−1), Ff = the fraction of total pore volume consisting of fractures, Cf = fracture compressibility (Pa−1), Ffm = the fraction of total pore volume that consists of matrix, and Cm = matrix pore volume compressibility (Pa−1).
The values of Cm were approximated by the correlation previously established by Teeuw [11], and the values of Cf were approximated with the following correlation established by Jones [12]. Cf is calculated as follows:
C f = 1 NOB   l   n ( NOB P H )
where NOB = net overburden pressure and PH = fracture healing pressure. The NOB was assumed to be 340 bars. This data came from the evaluation of the special core analyses. For PH, a value of 2721 bars was assumed, as suggested by Jones [12]. The fractured volume (Ff) value is based on the well test interpretation.
Table 4 shows the calculated pore volume compressibility. For an average porosity of 12 percent, the rock compressibility is 0.0310365 Pa−1. A compressibility of 0.02869 Pa−1 means that the volume of one million barrels of reservoir fluid will increase by 0.02869 m3 for every one pascal decrease in pressure.

3.3. Calculation of the Natural Water Influx

The pore volume compressibility is predicted by a simple model and is used to predict the water influx into a reservoir. The compressibility is defined mathematically as follows:
Δ V = C   V   Δ P
When the production of oil is conducted, a drop in the pressure of the reservoir results in the water in the aquifer expanding and entering the inside of the reservoir.
Ni = ( ra 2 re   2 )   h   Ø   5.615
Ne   = ( CW +   Cf )   Wi   ( Pi   P )
where We = cumulative water influx (m3), CW = aquifer water compressibility (Pa−1), Cf = aquifer rock compressibility (Pa−1), Wi = initial volume of water in the aquifer (m3), Pi = initial reservoir pressure (bar), P = current reservoir pressure (bar), ra = radius of the aquifer (m), re = radius of the reservoir (m), h = thickness of the aquifer (m), and Ø = percentage of porosity.
Since water influx occurs when the reservoir pressure drops due to a long period of fluid production, the water influx rate in the general material balance equation must be equal to the volumetric withdrawal rate. This is caused by the aquifer water influx expansion and influx into the reservoir, so the water influx average, or dWe/dt, = [average of active oil volumetric voidage] + [average of free gas volumetric voidage] + [average of water volumetric voidage] [13].
When written in terms of single-phase oil volume factors [14], the formula is as follows:
dWe dt = Bo dNp dt + ( R Rs ) dNp dt Bg + dWp dt Bw
where dNp/dt = daily oil rate (m3/d), dWe/dt = rate of water influx (m3/d), dNp/dt = the daily free gas rate (m3/d), R = current gas–oil ratio (m3/d), Rs = solution gas–oil ratio subtracted from the net daily or current gas–oil ratio (m3/d), and Bg = gas volume factor of the gas voidage term.
In the case of the steady-state aquifer, the rate of water flow dWe/dt is directly proportional to the decline in reservoir pressure (Pi–P). This assumes that the pressure is maintained at the initial Pi value somewhere in the aquifer and that the flow to the reservoir is, according to Darcy’s law, proportional to the pressure differential assuming that the water viscosity, average permeability, and aquifer geometry remain constant. This is calculated as follows:
dWe dt =   K   ( Pi   P )
We =   K   0 1 ( Pi P ) dt
using the technique of Havlena and Odeh [15]. The following material balance equation is for the Hamzeh undersaturated reservoir, which is always producing above the bubble point:
Np   Bo = N ( B 0 Boi ) + We
This will be used as the following equation of a straight line:
If = NEo + We ,   then   We   =   F     NEo
or
F E 0 = N + We E 0  
  F = Np   Bo ( Bo Boi )
Eo = ( Bo Boi )
All the parameters included in the equation are taken from the well test and the core and PVT analyses for the Hummar and Shueib formations. The results of the OOIP calculations are presented in Table 5 and Table 6 for the Shueib and Hummar formations, respectively. The results from the material balance equation calculations match the water influx calculations extremely well.
The plot for F/Eo versus We/Eo should be linear with the intercept F/Eo = N (original oil in place) when We/Eo = 0 and the function has a unit slope. Table 4 and Table 5 show the reserve estimations of water influx Havlena and Odeh [15] for the Hummar and Shueib formations, respectively. Figure 4 and Figure 5 show that the reservoir mechanism is determined by plotting F/Eo versus water influx (We/Eo). From these figures, we can conclude that the reservoir mechanism is water drive because the shape of the curve is a straight line. The material balance equation calculations match the water influx calculations extremely well.
We = 0 t Δ P   dt = Cs Σ Δ P Δ t
where Cs is the water influx constant (m3/kg/m2) and ΔP is the pressure drop (bar).
The water influx is calculated using the unsteady-state water influx equation as illustrated by Van Everdingen, A.F. and W. Hurst [6].
We = U   Dp   WD ( tD )
where WD = the dimensionless water influx, tD = dimensionless time, and
U = 1.119 Ø h Ct rw2 × Θ/360
which is the aquifer constant for radial geometry (m3/kg/m2). Θ = the angle subtended by the reservoir circumferences, h = thickness of the pay zone (m), rw = well radius (m), and Ct = oil compressibility, volume per volume per pound per square inch (Pa−1).
A comparative study of oil recovery factor determination for edge and bottom water drive mechanisms using water influx models shows that aquifer volume and aquifer permeability have a linear relationship with both bottom and edge water drives [16]. Bottom water drive is more efficient than edge water drive, and therefore, bottom water drive reservoirs tend to give higher oil recovery compared to edge water drive reservoirs [16]. The results of this study show that the mechanism of oil production is water, that the relationship is linear, and that this is consistent with other research.
Water influx modeling is highly needed to predict the future performance of the reservoir. Huge uncertainties exist when modeling water influx due to several reasons. First, we rarely drill wells into an aquifer to gain the necessary information about the aquifer reservoir properties. Second, properties are frequently inferred from what has been observed in the reservoir, and finally, geometry and areal continuity of the aquifer itself is a big question [17]. Natural influx of water in oil reservoirs surrounded by water-bearing rocks (aquifers) plays a very important role in increasing oil recovery. Such reservoirs are referred to as water drive reservoirs. The pressure drop in a reservoir caused by production of hydrocarbons prompts the influx of water from the surrounding aquifer to offset it. This pushes more hydrocarbons towards the surface, thus enhancing the oil recovery [18].
The dimensionless water influx (WD), a function of dimensionless time (tD) for a range of the aquifer-to-reservoir radius, is found from the charts done by Van Everdingen, A. and W. Hurst [6]. Successive dimensionless time (tD) is calculated from the following diffusivity equation:
tD = 0.00634 Kx Tx/Ø μCo rw2
where Kx = permeability (millidarcies), T = time (days), Ø = porosity fraction, μ = viscosity centipoise (mPa.s), Co = oil compressibility, volume per pound per square inch (Pa−1), and rw2 = reservoir radius (m).
The water influx is then found using the following equation:
We = UΣΔP × Q(t)
where the Q(t) values are found from the dimensionless time tables. The results are presented in Table 7 and Table 8.
The prediction for the oil and water that is produced as a function of pressure is shown in Figure 6 and Figure 7. The reservoir pressure declines rapidly and continuously while the amounts of produced oil and water increase. This is because the amount of water that encroached into the reservoir is not enough to replace the produced oil and water.

3.4. The Mechanism of Oil Displacement by Water

After the prediction of the future reservoir performance under the active water drive mechanism, it is worthwhile to describe the mechanism of oil displacement by water. The fluid oil and water saturation distribution dip in the normal direction as the flood moves through the reservoir, so at each increment of oil production (Np) corresponding to a pressure step at different time intervals, the water and oil saturation will change at each point in the reservoir. The relative permeability to oil and water, which are themselves a function of the water saturation, will be distributed at the same point in the reservoir.
By depicting this distribution, we construct the mathematical model that describes the movement of water toward the highest point of the reservoir and its influence on the mechanism of oil and water production during the operating life of the reservoir. The parameters for this enable the depiction of the process of oil displacement by water.
In addition to water saturation distribution, the following parameter is calculated: the fractional water flow or the water cut at each time step. The fractional flow equation can qualitatively show the ratio of the flux rate of water to the total flux rate of fluids at definite periods and spaces in a linear system of water injection. The equation demonstrates the relation between the volume of water flux in every zone of the reservoir and the amount of water saturation in a homogenous two-phase isothermal system.
Initially, only oil was produced from the wells in the part of the reservoir where only the connate (immobile) water saturation prevailed. However, in the later stage of production, when the water saturation (Sw) becomes higher than the connate water saturation (Swc), oil and water are produced from the same well. Therefore, in the Hamzeh reservoir, the wells will produce oil only at 11.8% and 11.6% in the Hummar and Shueib formations, respectively. They will then produce oil with water cut percentages corresponding to a higher water saturation (Sw) until the residual oil saturation (Sor) is achieved, and only water will be produced. Thus, the fraction of water flow or water cut corresponding to any water saturation is calculated as follows:
Fw = Qw/Qo + Qw = 1/1+ Ko/Kwx µw/µo
where Ko and Kw are the permeabilities to oil and water at different saturations—as measured in the laboratory—of the core, µo and µw are oil and water viscosities, and Qo and Qw are the oil and water volume factors, which are obtained from the PVT analysis. At each cumulative phase of oil production (Np) and at every corresponding pressure step, the oil saturation that prevails in this part of the reservoir becomes as follows (the water saturation at this step increases as much as the oil saturation decreases):
So = ( N Np )   Bo ( 1 Swi ) NBoi
Thus, the water saturation is calculated as follows:
Sw = 1 So
Ossai, P.G. et al. [19] show the relationship between the relative permeability ratio and water saturations. This was input as data for the prediction of the water flooding scheme performance.
If the water saturation distribution at any point and the relevant water fraction flow or water cut (fw) are known, the water–oil ratio (WOR) can be calculated as follows:
WOR = fw 1 +   fw
A water cut (wco) is the ratio of water production to the production of hydrocarbon fluids. A water cut in exploration wells can be calculated using the following equations:
wco = WOR 1 + WOR
where wco = water cut for oil well and WOR = water–oil ratio, and KRwo = permeability ratio of water to oil fraction.
KRwo = Krw Kro
where Kro = relative permeability for oil and Krw = relative permeability for water.
WOR = fw 1 fw
wc = WOR WOR + 1 = fw fw + 1 fw =   fw
The calculated values for both of the Hummar and Shueib formations are given in Table 9 and Table 10. “Sw” and “fw” are the water saturation and fractional flow at the producing end of the system, and “Pv” is the water invasion at the portion fraction of pore volume.
Table 11 shows the calculations of fractional water flow, fractional flow slope, and water influx.
A.
Water breakthrough time
Water breaks through to the producing wells when a certain value of water saturation is achieved. This is called breakthrough water saturation, and it is obtained using the technique used by S.E. Buckley and M.C. Leverette [20] as well as Welge, H.J [21]. This technique can estimate when the water breakthrough can occur. The time of water breakthrough is a function of the water influx (Wi), itself a function of pore volume during the water flood. The water influx is defined as follows:
Wi = ( 1 / dfw / dsw )   swbt
t = ( Wi × cum .   rock   vol . × porosity ) ( q × Bo × 5.615 × 365 )
where t is the time in years, Wi is the water influx, and Bo is the formation volume factor.
The time of breakthrough of water saturation (swbt) depends on the flow rate q/d and the withdrawal rate from the reservoir, so the time of breakthrough is calculated in the range of 5.50 to 16.50 years and 5.40 to 11 years for the Hummar and Shueib formations, respectively.
B.
Waterfront distance from the water–oil connection or aquifer
Since the oil is dislodged by the water, the water–oil saturation starts to increase from the point of contact in the reservoir. The water saturation at any point in the reservoir is proportional with the net rate of water gain at that point.
At the initial condition before production starts, the water influx (We) equals zero and the time (t) equals zero. Therefore, at any given time after the beginning of production, as the reservoir pressure declines, water flows from the aquifer and moves toward the reservoir through the water–oil contact in the linear system. The rate of the water movement is directly proportional to the cross section of the water breakthrough area and the porosity of the reservoir; therefore, the position of the constant water saturation plane or the distance of the waterfront from the aquifer can be estimated for the Hummar and Shueib reservoirs.
The following equation uses the Buckley–Leverett dfw/dsw:
X = 5.615 × Qt Ø A × dfw dsw ×   sw
This is the Welge equation that is applicable up to the flood front (swf):
X = 5.615 × Qt Ø A × dfw dsw ×   swf
These calculations are for different flow rates at different time steps. We tested the oil in place, the driving mechanism, the water influx, reservoir performance, and oil displacement.
The Hummar and Shueib formations are not homogenous, and in order to predict the flow rates of wells and their effects on pressure distribution, we must consider the structural contour map. The map uses an aerial grid system for both reservoirs and is extremely important for the description of the average petrophysical properties. The hydrocarbons of the both the Hummar and Shueib reservoirs are simulated using an aerial grid system with 200 grids, with an area representing 62,500 sq. ft. or 15.444 acres for each.

3.5. Prediction Pressure and Flow Rate Distribution

At any point in the reservoir, the distribution of pressure with distances from the bore hole in the formation is independent of formation permeability, but the velocity of this distribution is directly dependent on the formation permeability between wells. If one of the wells is producing at rate Q1, it will induce a pressure drop at the other wells located in the same reservoir or drainage area.
Using the Darcy’s law rate equation, a pressure drop in the second well can be calculated with the following:
Pe   Pd = Q 1 ×   Uo   ×   Bo   × Ln   re   rd   7.08 ×   k   ×   h

4. Results

The results obtained during the analysis of the data from Table 1, Table 2 and Table 3 are presented in this text to make a suitable model of reservoir water influx. To understand the effects of water influx on oil production, the following steps must be implemented: delimit the drive mechanism, make a model of the water influx, and predict oil production.
A.
Delimitation of reservoir drive mechanism
Table 2 and Table 3 show the reserve estimations for water influx for the Hummar and Shueib formations, respectively. One method for estimating N = original oil in place is to apply Equation (1) at a number of different successive reservoir pressures. In each calculation, We = 0 and the calculated values of the tank oil in place were plotted, and Ncalc = calculated values of tank oil in place versus Np = cumulative oil produced. If there has been water encroachment, the Ncalc vs. Np relationship would slope upward to the right, indicating water drive and solution gas drive.
Figure 2 and Figure 3 show that the calculated oil in place increases with cumulative production. Since these reservoirs are undersaturated, this indicates water encroachment and that the drive mechanisms for both reservoirs, Shueib and Hummar, are a combination of water drive and expansion drive.
B.
Model, parameter determination, and water influx modeling as shown in Table 4 and Table 5 show the reserve estimation of water influx, Havlen and D.S. Odeh model [15] for Hummar and Shueib formations, respectively. During the interpretation of the data, it appears that the appropriate water influx model for this case is the Pot aquifer model. Therefore, a plot of the term (F/Eo) against (We/Eo) is shown in Figure 4 and Figure 5. This plot produced a straight line. The value of the true N may be substituted in Equation (1), and the equation is solved for the water influx “(We)” at identical reservoir pressures. This would make water influx a function of Np, pressure, or time. Figure 4 and Figure 5 are graphs of the plots for F/Eo versus We/Eo, which determines what the reservoir mechanism is. From these figures, we can tell that the reservoir mechanism is water drive. This is due to the shape of the curve, which is a straight line. The material balance equation calculations match the water influx calculations extremely well. The reservoir energy is from fluid expansion and also accounts for the presence of water influx as source of energy.
C.
Prediction of reservoir behavior
The water influx is calculated using Equation (22), which is for the unsteady-state water influx. The calculation points of F/Eo vs. We/Eo increase as demonstrated. This indicates that the reservoir has been motivated by water influx. The water influx that was calculated with the material balance equation was first calculated utilizing both models of the Schilthuis method ((22) and (23)) and the unsteady-state equation. Results of the water influx are presented in Table 6 and Table 7. Performance predictions were made using the unsteady-state water influx model (23). Figure 8 and Figure 9 show these results, and the points almost make a straight line. Therefore, the water influx that results from each method (MBE and the unsteady-state model) corresponds very well with the reserve estimation that uses water influx calculations along with the D. Havlena and A.S Odeh techniques [15]. for the Hummar and Shueib formation oil reservoir.
D.
Predict oil production
The time of breakthrough of water saturation or the time of achieving such (swbt) depends on the flow rate q/d and the withdrawal rate from the reservoir, so the time of breakthrough is calculated in the range from 5.50 to 16.50 years and 5.40 to 11 years for the Hummar and Shueib formations, respectively.

5. Discussion

Table 1 show the characteristic reservoir; the reservoir is naturally fractured, and log–log plots and type-curve matching indicate the double porosity behavior of the reservoir (transient interporosity flow). Fractures are mostly vertical as obtained by core samples analysis. Initial reservoir pressures in the area are estimated to be 268.36 bar at 9629.73 ft drilled depth in Hummar formation, bottom hole pressure of 264.41 bar 9491.93 ft drilled depth, formation thickness of 16.4 ft, bubble point pressure of 51 bar, oil viscosity 3 cp, formation volume factor 1.075 vol/vol, permeability 132 mD, water saturation 17%, and average flow rate 34.57 m3/day. In the Shueib formation, the fractured permeability around the wellbore is very high compared with Hummar formation. The fractures are horizontal, fractured permeability Kf = 2083 mD, average Permeability Kav = 290 mD, average permeability of all reservoirs about 285 mD, and formation volume factor is 1.141 vol/vol. Initial reservoir pressures in the area are estimated to be 262.58 bar at 9485.37 ft drilled depth in Shueib formation, and bottom hole pressure 244.89 bar at 9508.33 ft drilled depth. Formation thickness is 26.25 ft, bubble point pressure is 68 bar, and production rate 51.62 m3 /day. The gas oil ratio (GOR) in the Shueib formation is relatively higher, 2.148 sm3/m3 compared to 0.402 sm3/m3 in the Hummar formation. There is an indication of water oil contact observed by electrical logging; in addition to that, in most Hamzeh wells the Shueib formation produced water. This study confirmed the existence of water and improved the existing understanding of drive mechanisms and fluid expansion in Hamzeh oil field.
The reservoir power was supplied mainly by rock and fluid expansion assisted by a certain volume of water influx. However, the water influx greatly influences reservoir pressure, thus affecting predictions for production. Oloro and Ukrakpor [22] worked on the determination of water influx in the Niger Delta, and they proposed that, to determine the total water influx into a reservoir at any given time, it is necessary to determine the water influx as a result of each successive pressure drop that has been imposed on the reservoir and aquifer. In calculating accumulative water influx, the total water influx must be calculated from the start. This is required for the different times during which the various pressure declines were effective, which is in line with the approach used in this study.
The water influx was predicted using the material balance equation and the Van Everdingen–Hurst unsteady-state equation, and the results for the cumulative oil production and water influx prediction were greater for the model than for the material balance equation, which does not bear all points obtained from the results in the study [1]. At the same time, a prediction for water influx was made by using the material balance equation and the Van Everdingen–Hurst unsteady-state equation. The results for the accumulative oil production and water influx prediction for the MBE and the unsteady-state model have an average absolute error of 1.214%. The material balance model is 9611.90 m3 and 9287.61 m3 with the unsteady-state model. The MB equation has greater priority than the equation model but not at all points with the results obtained in this study.
A Brownfield development plan was investigated by mining the data pack from the Equinor Volve field repository, and models and simulations were used to establish the quantity of hydrocarbon that is still recoverable therefrom [23]. In this study, models and simulations must be used to establish the quantity of hydrocarbon in the Hamzeh oil field.
A.
The mechanism of oil displacement by water
Utilizing the Buckley–Leverette frontal advance theory, we evaluated the water saturation distribution and residual oil saturation during the immiscible displacement process as well as the oil displacement efficiency. The displacement sweeps at each WOR were calculated from the relative permeability data shown in Table 9 and Figure 6 and Figure 7 for the Shueib and Hummar reservoirs, respectively.
B.
Fractional water flow or the water cut at each time step
By using Equations (32)–(34), the water cut value can be calculated. The resulting fractional flow values from these calculations are summarized in Table 9. In Figure 10, a tangent to the fractional flow curve from irreducible water saturation (Swir) is shown, and the point of tangency represents the water saturation at the front (Swf) and the value of the fractional flow for producing the water cut.
Table 9 shows the resulting relative permeabilities and the fractional flow data of oil. Figure 10 presents the variation in relative permeabilities and the portion of water flux that has the change in water saturation. The fractional water flow (fw) was calculated with Equation (26). The relative permeability diagram has curves of oil and water present in its usual curves, and the endpoint rate of relative water permeability was a fraction of 0.272. At this point, only water is injected into the reservoir.
These results indicate that the fractional flow amount curve was a 0.54 fraction of the oil manifests. The uniform method shows that a major quantity of movable oil in the reservoir is displaced by water injection. In addition, in Figure 10, the Sw versus fw relationship shows a water saturation cutoff. Table 10 presents the resulting relative permeabilities, the fractional flow, and the water cut values of oil.
The diagram of fractional water flux or the water cut versus water saturation was then plotted. Hence, the cutoff value for the water saturation for the whole oil field was realized. The diagram for fractional flow versus water saturation is shown in the following figure. In Figure 11, the Sw versus WC relationship also shows a water saturation cutoff of 40%. Utilizing fractional flow and water cut information together, we were able to notice that there is only small difference between the two values. This evidence shows that fractional flow is equivalent to the water cut; Fw = 0.751 and WC = 0.735. The resulting fractional flow value from the calculation above is summarized in Table 9.
The oil and water relative permeability curves under different displacement pressures were compared, which is shown in Figure 12 and Figure 13. Its ability is visible when the experiment’s displacement pressure is incremental, the range of the two-phase zone raises, and the residual oil saturation lowers. The relative permeabilities of oil and water rise, and the curve shifts to the upper right region. As the displacement pressure increases, the range of the pore sizes that permit the fluid to flow increases, so the remaining oil saturation will decrease and the oil recovery factor will be enhanced.
Next, the diagram of fractional flow or water cut versus water saturation was plotted. Hence, the cutoff value of the water saturation for the field and the fractional flow value equivalent to the water cut value can be obtained, as previously explained. Utilizing both fractional flow and water cut information, it can be observed that there is only a small difference between the two values. This shows that the fractional flow is equal to the water cut.
According to both the cumulative and average of production oil and water, daily oil and water production is calculated for both the Hummar and Shueib formations. Using the different WOR values of 0.79 m3, 3.97 m3, and 15.89 m3, water saturation (Sw) is obtained from the fractional flow curves for each corresponding fractional water flow (fw), or in other words, for each WOR.
Table 9 and Figure 13 show the changes in relative permeabilities and fractional water flow along with the change in water saturation. The fractional water flow (fw) was calculated with Equation (26). Relative permeability curves of oil and water present ordinary curves, and the endpoint rate of relative water permeability was the 0.320 fraction.
At this point, only water was pumped into the reservoir. The results indicate that the fractional flow value curve was 0.673 fraction for oil, manifesting a regular shape, and the application of this curve to the water flooding method shows that a large amount of movable oil in the reservoir is displaced by water injection. In addition, in Figure 14, the Sw versus fw relationship shows a water saturation cutoff of 40%.
In Figure 15, the Sw versus WC relationship also shows a water saturation cutoff of 40%. Utilizing fractional flow and water cut information together, it is possible to notice that there is only one small variance between the two values. This indicates that fractional flow is equal to the water cut.
Figure 10 presents the evidence for the fractional flow curve in the displacement of oil by water: the curve has an elongated fluid (Swi < Sw < 1-Sor). The stage of saturation at the tangent point of a straight line drawn from the irreducible saturation on the fractional flow curve is utilized to define the saturation rate at the waterfront, which conforms to the Buckley–Leverett theory.
A fractional flow curve (Fw versus Sw) is utilized to characterize the immiscible fluid displacement operation. Developing an illustrative fractional flow curve for a specific reservoir perfectly may be challenging after a liquid and special core information analysis has been performed. Then, the technique that measures fractional flow curve from the historic production information must be used.
The curve becomes the primary one for a similar model that permits the estimations of oil rate production and reserves for present or proposed wells. The relative permeability curves can be formed mainly due to the resultant fractional flow curve. A comparison of relative permeability curves can be made from the special core analysis information, which may raise confidence.
By extrapolating this tangent on the fractional flow curve to 1, Figure 10 and Figure 11 will yield the value of the average water saturation at breakthrough. The values that were obtained are as follows: fw = 0.56 for the Hummar formation, fw = 0.674 for the Shueib formation, and SWf = 0.21 and 0.26 for Hummar and Shueib formations, respectively.
Figure 16 shows the extrapolation of the WOR versus Np plot and changes in its decline of the ability to indicate increased oil recovery. Therefore, an inspection of the WOR versus Np plot is helpful in defining the increased recovery due to infill drilling or operational changes, as shown in Figure 16.
The variable slope of the curve indicates increased reserves after infill drilling in additional successful recovery efforts (including completions or treatments to restrain water as well as infill drilling).
As shown in Figure 17, the reservoir pressure declines rapidly and continuously with the amount of production time. Due to this, the water that enters the reservoir is not enough to replace the produced oil.
A.
Water breakthrough time
Table 12 and Figure 18 and Figure 19 show that the oil production rate was varied, beginning from a principle state of 15.89 to 23 m3/day. It moved to 84 m3/day, and then to a range of 31.79 m3/day to 55.64 m3/day. This resulted in six simulation runs for the Hummar and Shueib formations, respectively.
This relation is visible in the diagram that shows performance of and change in oil production rate relative to breakthrough time. This increase in oil rate is inversely proportional to the breakthrough time. This is indicative of the actuality that increasing the production rate will result in a single breakthrough time and coning taking place.
However, the cumulative oil recovery increases with the influx rate, and the volume of oil produced before the synchronous production of the reservoir fluids is affected by the oil column height below and above perforations at breakthrough time. Thus, as the production rate rises, a smaller amount of “water-free” oil and “gas-free” oil will be produced before the ultimate concurrent coning into the well occurs in spite of the oil recovery increase along with the influx rate. This further backs up the actuality that rising the oil production rate will raise the rate of concurrent coning into the well.
The impact of oil influx rate on the slope of the WOR curve is minimal. This indicates that the differences in influx rates influence only the times at which water coning starts, or when water breaks into the well, and changes in the oil influx rate have a minimal effect on gas–oil ratio (GOR) and WOR. The time of breakthrough or the time of achieving the swbt depends on the flow rate q and the withdrawal rate from the reservoir, so the time of breakthrough is calculated in a range from 5.50 to 16.50 years and 5.40 to 11 years for the Hummar and Shueib formations, respectively.
Where time is expressed in years, “We” is water influx, and the average water influx range is between 18.4% and 21.4%. The breakthrough water saturation was 21% to 26% for the Hummar and Shueib formations, respectively. The oil displacement efficiency (ED) for the Hummar formation was 0.85 to 0.94% and was 0.90 to 0.99% for the Shueib formation. The time of water breakthrough, assuming cumulative production rate, with q = 160 to 240 m3/d, was calculated in a range between 5.65 and 16.50 years for the Hummar formation and from 5.40 to 11 years for the Shueib formation. The prediction of oil recovery performance under the energy of water influx is the main portion of the actual reservoir concerning oil displacement by water and recovery factors.

6. Conclusions

The problem that arises is how to know whether the quantity of water influx in the reservoir influenced future performance. With the master method of the material balance equation being applied to a group of available data (pressure, volume, and temperature data, production and reservoir data). Drive mechanisms are determined by the analysis of historical production data, primarily reservoir pressure data and fluid production ratios [24]. A solution to this problem was found using a methodology that was defined according to the main objectives, which consisted of the following:
  • Identify the drive mechanism to discover the presence and size of the aquifer so as to know which energy mainly drives the fluids.
  • Elect a suitable water influx model, water influx modeling method, and oil production prediction. The results showed that the drive mechanism (fluid and rock expansion) for Hamzeh oil reservoir has been confirmed. This study confirms that the reservoir energy comes from fluid expansion and also accounts for the presence of water influx as a source of energy. Similarly, the results of the oil production prediction seem to be confirmed by the total increase in production. The oil production rate was varied, beginning from a principle state of 15.89 to 23 m3/day. It moved to 84 m3/day and then to a range of 31.79 m3/day to 55.64 m3/day, resulting in six simulation runs for the Hummar and Shueib formations, respectively, which increased the recovery factor accompanied by a slight decrease in pressure.
  • The water breakthrough time was calculated to be in a range between 5.65 and 16.50 years for the Hummar formation and 5.40 to 11 years for the Shueib formation.
  • The fractional flow curve (fw) will yield the value of average water saturation at breakthrough. These values obtained were fw = 0.56 for the Hummar formation, fw = 0.674 for the Shueib formation, and SWf = 0.21 and 0.26 for the Hummar and Shueib formations, respectively. A fractional flow curve (Fw versus Sw) is utilized to characterize the immiscible fluid displacement operation.
  • The recovery factor will be increased to 25% in Shueib formation, relatively higher than the 23% in Hummar formation. The results in this study were confirmed with other researchers in the same field. That natural influx of water in oil reservoirs surrounded by water-bearing rocks (aquifers) play a very important role in increasing oil recovery. This study examines the Alwyn North Field and presents a comparative analysis for producing crude oil using water drive, gas drive, and natural drive or depletion drive; profitability index for water injection was 1.73 and was 1.36 for gas injection. These economic tools help us to arrive at the conclusion that the best drive mechanism for the production of the field was water drive [25].

Author Contributions

Conceptualization, M.A.-M. (Mehaysen Al-Mahasneh) and H.E.A.-K.; methodology, K.A.-Z. and M.A.-M. (Marwan Al-Mahasneh); validation, M.A.-M. (Mehaysen Al-Mahasneh), M.A.-M. (Marwan Al-Mahasneh) and A.A.; formal analysis, M.A.-M. (Mehaysen Al-Mahasneh), H.E.A.-K., M.A.-M. (Marwan Al-Mahasneh) and A.A.; investigation, H.E.A.-K. and K.A.-Z.; resources, M.A.-M. (Mehaysen Al-Mahasneh); data curation, M.A.-M. (Mehaysen Al-Mahasneh), M.A.-M. (Marwan Al-Mahasneh) and A.A.; writing—original draft preparation, M.A.-M. (Mehaysen Al-Mahasneh); writing—review and editing, M.A.-M. (Mehaysen Al-Mahasneh), M.A.-M. (Marwan Al-Mahasneh) and H.E.A.-K.; visualization, K.A.-Z. and M.A.-M. (Mehaysen Al-Mahasneh); supervision, M.A.-M. (Mehaysen Al-Mahasneh), M.A.-M. (Marwan Al-Mahasneh), H.E.A.-K. and K.A.-Z.; project administration, M.A.-M. (Mehaysen Al-Mahasneh). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Np Cumulative oil produced, Sm3
PReservoir pressure, bar
Pi initial reservoir pressure, bar
∆p pressure drop, bar
reradius of the reservoir, m
TAbsolute temperature, °C
BWater influx constant, m3/bar
Booil formation volume factor at P, m3/ Sm3
BoiInitial oil formation volume factor, rm3/sm3
BwWater formation volume factor, rm3/sm3
TAbsolute temperature, °C
t Time, days
tDDimensionless time
WeCumulative water influx, sm3
WeDDimensionless water influx
µwViscosity of water in the aquifer, mPa.s
SwWater saturation %
FwFractional flow %
SoOil Saturation %
CfFormation (rock) compressibility, m3/m3
CtTotal compressibility coefficient, Pa−1
CwWater compressibility, pa- 1
Rsis dissolved gas oil ratio, sm3/sm3
Rsigas injection sm3/sm3
w Expansion of formation and water sm3/sm3
Eo Expansion of oil and its originally dissolved gas, sm3/sm3
NOriginal oil in place, sm3
WCWater cut %
reRadius of the reservoir, m
ϕWater injection, sm3
WpCumulative water produced, sm3
θAngle subtended by the reservoir circumference
AArea grid system for both reservoirs, m2
ϕ porosity of the reservoir in fraction
µo Viscosity of oil in the aquifer, mPa.s
SwbtBreakthrough of water saturation
Swf Water saturation at the front %
bPBubble point pressure, bar

References

  1. BinMerdhah, A.B.; Muherei, M.A.; Al-Khudafi, A.M.; Bamumen, S.S. Prediction of Water Influx Effect on Oil Production in Biyad Oil Reservoir. Hadhramaut Univ. J. Nat. Appl. Sci. 2015, 12, 9–21. [Google Scholar]
  2. Schilthuis, R.J. Active Oil and Reservoir Energy. Trans. AIME 1936, 118, 33–52. [Google Scholar] [CrossRef]
  3. Hurst, W. Water Influx Into a Reservoir and Its Application to the Equation of Volumetric Balance. Trans. AIME 1943, 151, 57–72. [Google Scholar] [CrossRef]
  4. Fetkovich, M. A Simplified Approach to Water Influx Calculations-Finite Aquifer Systems. J. Pet. Technol. 1971, 23, 814–828. [Google Scholar] [CrossRef]
  5. Al-Mahasneh, M.; Al Rabadi, S.; Khaswaneh, H. Assessment of oil-producing wells by means of stimulation approach through matrix acidizing: A case study in the Azraq region. J. Pet. Explor. Prod. Technol. 2021, 11, 3479–3491. [Google Scholar] [CrossRef]
  6. Van Everdingen, A.; Hurst, W. The Application of the Laplace Transformation to Flow Problems in Reservoirs. J. Pet. Technol. 1949, 1, 305–324. [Google Scholar] [CrossRef]
  7. Halafawi, M.; Abdel Waly, A. Prediction Modeling for Combination Drive Reservoir Performance. Adv. Pet. Explor. Dev. 2018, 16, 49–71. [Google Scholar] [CrossRef]
  8. Yang, Y.; Yu, J.; Wang, Y.; Ma, C. Optimizing and accelerating history matching progress of numerical res-ervoir simulation by using material balance analysis. MATEC Web Conf. 2021, 336, 01019. [Google Scholar]
  9. Belomo, V.; Nitcheu, M.; Dongmo, E.D.; Chamgoué, A.C.; Kuiatse, G.; Kingni, S.T. Estimation of Water Influx and Oil Reservoir Performance Prediction in the Volve Oil Field. Libyan J. Eng. Sci. Technol. (LJEST) 2022, 2, 69–72. [Google Scholar]
  10. Mahasneh, M.; Banihani, F.; Hussam, K.H. Steam injection in porous media: Case study Wadi-Rajil, Jordan. Jordan. J. Eng. Chem. Ind. (JJECI) 2020, 3, 166–169. [Google Scholar]
  11. Teeuw, D. Prediction of Formation Compaction from Laboratory Compressibility Data (SPH Paper 2973). Soc. Pet. Eng. J. 1975, 11, 263–271. [Google Scholar] [CrossRef]
  12. Jones, F.O.J. A Laboratory Study of the Effects of Confining Pressure on Fracture Flow and Storage Capacity in Carbonate Rocks. J. Pet. Technol. 1975, 27, 21–27. [Google Scholar] [CrossRef]
  13. Ahmed, T. Reservoir Engineering Handbook, 4th ed.; Gulf Professional Publishing: Oxford, UK, 2010. [Google Scholar]
  14. Craft, B.; Hawkins, M.; Terry, R. Applied Petroleum Reservoir Engineering, 2nd ed.; Printice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  15. Havlena, D.; Odeh, A. The Material Balance as an Equation of a Straight Line. J. Pet. Technol. 1963, 15, 896–900. [Google Scholar] [CrossRef]
  16. Nmegbu, C.G.J.; Ebube, O.F.; Edet, E.A. Comparative Study of Oil Recovery Factor Determination for Edge and Bottom Water Drive Mechanism Using Water Influx Models. Eur. J. Eng. Technol. Res. 2021, 6, 1–9. [Google Scholar] [CrossRef]
  17. Yousef, M. Water Influx from the Aquifer: Getting into Basics 2019. Available online: https://www.linkedin.com/pulse/water-influx-from-aquifer-getting-basics-mohamad-yousef-alklih. (accessed on 21 December 2022).
  18. Chauhan, M.; Maheshwari, U. Calculation of Water Influx and Its Impact on Oil Recovery in Edge Water Drive Oil Reservoirs (Rapid Fire). In Proceedings of the 2015mAIChE Annual Meeting, Salt Lake City, UT, USA, 8–13 November 2015. [Google Scholar]
  19. Ossai, P.G.; Onwukwe, S.; Izuwa, N.C.; Noaemeka, N.C. Waterflooding in Niger Delta Oil Reservoir Systems: Concept and Analysis. J. Adv. Res. Pet. Technol. Manag. 2018, 4, 36–41. [Google Scholar]
  20. Buckley, S.; Leverett, M. Mechanism of Fluid Displacement in Sands. Trans. AIME 1942, 146, 107–116. [Google Scholar] [CrossRef]
  21. Welge, H.J. A Simplified Method for Computing Oil Recovery by Gas or Water Drive. J. Pet. Technol. 1952, 4, 91–98. [Google Scholar] [CrossRef]
  22. Oloro, J.; Ukrakpor, E.F. Determination of water influx in reservoir in Niger delta. Cont. J. Eng. Sci. 2011, 6, 37–44. [Google Scholar]
  23. Mcallister, A.E.; Lateef, A. Planning and implementing the redevelopment of volve field, Norwegian north sea using knowledge discovery in databases data mining and numerical simulation techniques. J. Eng. Sci. Technol. 2022, 17, 2567–2593. [Google Scholar]
  24. Rhodes, E.G.; Whitson, C.H.; Golan, M.; Sustakoski, R.J.; Morton-Thompson, D.; Sills, S.R.; Sarem, A.M.S.; Breit, V.S.; Weber, K.; Laudeman, S.K. Reservoir Engineering Methods; American Association of Petroleum Geologists: Tulsa, OK, USA, 1992. [Google Scholar] [CrossRef]
  25. Amadi, A.H.; Raphael, O.; Ebube, O.F. Comparative Analysis of Water Drive, Gas Drive and Natural Drive Mechanism for Oil Production using Alwyn North Field as a Case Study. Eur. J. Eng. Technol. Res. 2020, 5, 479–484. [Google Scholar] [CrossRef]
Figure 1. The lithostratigraphic cross section of the Hamzeh field [10].
Figure 1. The lithostratigraphic cross section of the Hamzeh field [10].
Energies 16 02126 g001
Figure 2. Calculated oil in place vs. cumulative oil production in the Shueib formation.
Figure 2. Calculated oil in place vs. cumulative oil production in the Shueib formation.
Energies 16 02126 g002
Figure 3. Calculated oil in place vs. cumulative oil production in the Hummar formation.
Figure 3. Calculated oil in place vs. cumulative oil production in the Hummar formation.
Energies 16 02126 g003
Figure 4. Reserve estimation from water influx calculations for the Hummar formation.
Figure 4. Reserve estimation from water influx calculations for the Hummar formation.
Energies 16 02126 g004
Figure 5. Shueib formation reserve estimations using water influx calculations.
Figure 5. Shueib formation reserve estimations using water influx calculations.
Energies 16 02126 g005
Figure 6. Pressure versus oil produced.
Figure 6. Pressure versus oil produced.
Energies 16 02126 g006
Figure 7. Pressure versus water produced.
Figure 7. Pressure versus water produced.
Energies 16 02126 g007
Figure 8. Unsteady-state model versus MBE for Hummar formation.
Figure 8. Unsteady-state model versus MBE for Hummar formation.
Energies 16 02126 g008
Figure 9. Unsteady-state model versus MBE for Shueib formation.
Figure 9. Unsteady-state model versus MBE for Shueib formation.
Energies 16 02126 g009
Figure 10. Fractional flow curve for the displacement of oil by water.
Figure 10. Fractional flow curve for the displacement of oil by water.
Energies 16 02126 g010
Figure 11. Fractional flow or water cut diagram curve with water saturation for the displacement from oil by water.
Figure 11. Fractional flow or water cut diagram curve with water saturation for the displacement from oil by water.
Energies 16 02126 g011
Figure 12. Average water–oil relative permeability versus water saturation for Shueib formation.
Figure 12. Average water–oil relative permeability versus water saturation for Shueib formation.
Energies 16 02126 g012
Figure 13. Average water–oil relative permeability versus water saturation for Hummar formation.
Figure 13. Average water–oil relative permeability versus water saturation for Hummar formation.
Energies 16 02126 g013
Figure 14. Fractional flow curve with water saturation in the displacement of oil by water.
Figure 14. Fractional flow curve with water saturation in the displacement of oil by water.
Energies 16 02126 g014
Figure 15. Fractional flow or water cut with water saturation in the displacement of oil by water.
Figure 15. Fractional flow or water cut with water saturation in the displacement of oil by water.
Energies 16 02126 g015
Figure 16. Total production versus water–oil ratio (WOR).
Figure 16. Total production versus water–oil ratio (WOR).
Energies 16 02126 g016
Figure 17. Water–oil ratio (WOR) versus static pressure.
Figure 17. Water–oil ratio (WOR) versus static pressure.
Energies 16 02126 g017
Figure 18. Water breakthrough daily production for the Hummar formation.
Figure 18. Water breakthrough daily production for the Hummar formation.
Energies 16 02126 g018
Figure 19. Water breakthrough daily production for the Shueib formation.
Figure 19. Water breakthrough daily production for the Shueib formation.
Energies 16 02126 g019
Table 1. Rock and fluid characteristics of reservoir.
Table 1. Rock and fluid characteristics of reservoir.
Properties ShueibHummarUnit
Porosity 1511%
Permeability 290132md
Oil gravity31.86°32°API°
Gas gravity0.7410.741sp. gravity
Gas–oil ratio (produced)4.9945350Nm3/m3
Check pressure262.6 260.7 bar
Check temperature90.56 90.55 °C
Water compressibility2.12 × 10−22.12 × 10−2Pa−1
Oil saturation6056%
Water saturation 2017%
Oil density0.86620.8662gr/cm3
Oil viscosity0.00217 0.00217 mPa.s
Oil formation volume factor0.16770.1814m3/m3
Bubble point pressure 6851Bar
Thickness26.516.7ft
Table 2. Calculation of material balance for the Hummar formation using Hummar PVT data, where Swi = 11.8%, Cr = 0.01620 Pa−1, and N = 6.12 × 106 m3..
Table 2. Calculation of material balance for the Hummar formation using Hummar PVT data, where Swi = 11.8%, Cr = 0.01620 Pa−1, and N = 6.12 × 106 m3..
P Bar Δp Bar Np m3 Coe Pa−1Ncalc × 106 m3Np × BoNBoi Coe ΔpWe m3
271.29-------
267.82 3.46 4947.210.0468014.297522,4151,045,260.70 2984.33
267.553.745789.180.0444916.1626113.371,214,199.593799.79
267.343.948955.640.0522820.3739448.201,894,717.826609.04
266.664.6215,898.250.0499932.25916,788.553,363,827.8913,589.00
Table 3. Calculation of material balance for the Shueib formation using Shueib PVT Data, where N = 3.25 × 106 m3 Swi = 11.6%, and Cr = 0.01337 Pa−1.
Table 3. Calculation of material balance for the Shueib formation using Shueib PVT Data, where N = 3.25 × 106 m3 Swi = 11.6%, and Cr = 0.01337 Pa−1.
P Bar Δp Bar Np m3Coe × 10−2 Pa−1Ncalc × 106 m3Np × BoNBoi Coe ΔpWe m3
265.30-------
263.262.04 528.130.044492.72602.59 106,399.29 -
261.224.081901.740.039455.542169.88384,321.07894.05
259.18 6,12 4329.88 0.037118.94 4940.39875,098.953141.97
254.420.88 21,115.42 0.0356825.51 24,092.694,268,167.6521,030.31
Table 4. Calculated pore volume compressibility.
Table 4. Calculated pore volume compressibility.
Formation Porosity % Ff %Crt, V/v/Pa × 10−2
Shueib 12.2500.0137
Shueib12.251.9 0.0239
Hummar 10.4500.0162
Hummar10.751.90.0286
Table 5. Reserve estimation using water influx calculation for Hummar formation.
Table 5. Reserve estimation using water influx calculation for Hummar formation.
Time DaysPavg
Bar
Δ P
Bar
Np m3Bo
m3/m3
F
m3
Eo
m3/m3
F/Eo
m3
tD/tTDUQtWe
m3
We/Eo
0271.290-0.12576---- 121.3---
100267.683.64897.210.125795166.210.00002621.9 × 1080.04974.970121.34.54599.921.7 × 108
117267.343.945805.630.125796124.700.00002862.1 × 1080.04975.815121.355593.231.9 × 108
404266.874.4215,873.670.1257916746.470.00003215.2 × 1080.049720.08121.312.415,545.384.8 × 108
660266.734.5522,822.210.1257924077.360.00003337.2 × 1080.0497-32.80121.317.823,130.956.9 × 108
Table 6. Reserve estimations with water influx calculations using the Havlena–Odeh technique [13] (Shueib formation).
Table 6. Reserve estimations with water influx calculations using the Havlena–Odeh technique [13] (Shueib formation).
Time Days Pavg
Bar
Δ P
Bar
Np
m3
Bo
m3 /m3
F
m3
Eo
m3/m3
F/Eo
m3
tD/tTDU
m3/bar
Qt
m3
WeWe/Eo
0265.300-0.1360-- ---32.58---
8263.601.76528.130.1360602.690.00001195 × 1070.06840.5432.581.063236.83 × 107
33261.903.401901.710.13602170.390.00002508 × 1070.06842.2532.582.515,2687 × 107
76261.224.084329.820.13604941.730.00002981.6 × 1080.06845.2032.584.633,7121 × 108
377260.275.0321,115.080.136124,100.730.00003696 × 1080.068425.732.5815135,5824 × 108
581259.795.5129,996.660.136134,239.200.00004058 × 1080.068439.732.5820.7204,8026 × 108
Table 7. Results of water influx match for Hummar formation.
Table 7. Results of water influx match for Hummar formation.
Time DaysPressure BarMaterial Balance Water Influx m3Calculated Water Influx m3, Steady StateCalculated Water Influx m3,
Unsteady State
134267.89523.04324.95379.00
305265.911325.251395.201513.33
520264.892062.763199.633245.25
540264.422885.163388.813429.65
654263.404538.084588.794551.91
712262.786413.205372.255280.67
Table 8. Results of water influx match for Shueib formation.
Table 8. Results of water influx match for Shueib formation.
Time DaysPressure Bar Material Balance Water Influx m3 Calculated Water Influx m3, Steady StateCalculated Water Influx m3, Unsteady State
25263.12161.52171.22213.65
46261.42582.18574.32687.11
97259.521712.532121.902293.28
112258.502558.142730.022907.58
160257.144590.735002.535095.62
187256.056954.976600.056599.41
237255.109611.909505.259287.61
Table 9. Resulting relative permeabilities, fractional flow data of oil. sw versus fw for Shueib formation, μo = 2.54 cp, μw = 0.3 cp.
Table 9. Resulting relative permeabilities, fractional flow data of oil. sw versus fw for Shueib formation, μo = 2.54 cp, μw = 0.3 cp.
sw%fwwckroKrw
11.60.0000.0001.000.00
20.00.3210.3460.4830.027
25.00.6310.6590.3070.062
30.00.8040.8240.18330.089
35.00,9130.9230.1000.125
40.00.9680.9720.430.1567
45.0 0.9920,9930.0120.193
50.00.9970.9980.0040.230
56.01.000.0.99920.0000.272
Table 10. Resulting relative permeabilities, fractional flow and water cut values of oil. Sw versus fw for Hummar formation, μo = 2.8 cp, μw = 0.3 cp.
Table 10. Resulting relative permeabilities, fractional flow and water cut values of oil. Sw versus fw for Hummar formation, μo = 2.8 cp, μw = 0.3 cp.
swfwwcKroKrw
11.80.0000.4691.0000.000
20.00.45730.6740.5870.053
25.00.66290.8150.4130.087
30.00.80730.8980.2740.123
35.00.89320.9600.1730.155
40.00.95890.9930.0800.200
450.99280.9980.0170..253
521.00000.99940.0000.320
Table 11. Summary calculation of the data of the breakthrough of water saturation.
Table 11. Summary calculation of the data of the breakthrough of water saturation.
FormationSwbtfwdfw/dswWi
Hummar0.21o.5735.4350.184
Shueib0.260.6744.6810.214
Table 12. The results of calculation for breakthrough time.
Table 12. The results of calculation for breakthrough time.
FormationSwCum.rock vol.
m3 × 106
Prod.Rate
m3/d
Water Breakthrough
Time Years for Hummar Formation
Prod.Rate m3/dWater Breakthrough
Time Years for Shueib Formation
Hummar0.11828.1215.8916.5014.8011
23.8411.0021.408.50
Shueib0.11621.2331.798.30029.308.20
39.746.60035.705.65
47.695.50043.655.40
55.644.70052.604.70
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al-Mahasneh, M.; Al-Khasawneh, H.E.; Al-Zboon, K.; Al-Mahasneh, M.; Aljarrah, A. Water Influx Impact on Oil Production in Hamzeh Oil Reservoir in Northeastern Jordan: Case Study. Energies 2023, 16, 2126. https://doi.org/10.3390/en16052126

AMA Style

Al-Mahasneh M, Al-Khasawneh HE, Al-Zboon K, Al-Mahasneh M, Aljarrah A. Water Influx Impact on Oil Production in Hamzeh Oil Reservoir in Northeastern Jordan: Case Study. Energies. 2023; 16(5):2126. https://doi.org/10.3390/en16052126

Chicago/Turabian Style

Al-Mahasneh, Mehaysen, Hussam Elddin Al-Khasawneh, Kamel Al-Zboon, Marwan Al-Mahasneh, and Ali Aljarrah. 2023. "Water Influx Impact on Oil Production in Hamzeh Oil Reservoir in Northeastern Jordan: Case Study" Energies 16, no. 5: 2126. https://doi.org/10.3390/en16052126

APA Style

Al-Mahasneh, M., Al-Khasawneh, H. E., Al-Zboon, K., Al-Mahasneh, M., & Aljarrah, A. (2023). Water Influx Impact on Oil Production in Hamzeh Oil Reservoir in Northeastern Jordan: Case Study. Energies, 16(5), 2126. https://doi.org/10.3390/en16052126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop