Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Production Test Analysis of XYZ-Well at Dieng Geothermal Field Using Horizontal Discharge Lip Pressure Method with Russel James Equation and Hiriart Equation Aris Tristianto Wibowo 1, Muhammad Thasril 2, and Puji Sirait 3 1 Petroleum Engineering, ITB, E-mail: aris.wibowo@outlook.com 2 Halliburton, E-mail: m.thasril@gmail.com 3 Halliburton, E-mail: pujigwa@yahoo.com Keywords: production test, horizontal discharge lip pressure method, Russel James equation, Hiriart equation ABSTRACT Production test is one of very crucial activities in the geothermal field development process because the condition of the well, reservoir characteristics, and production capacity of the well can be determined by the test. Production test data should be processed carefully to get the best possible determination of production capacity approaching the actual condition of the well, so that future development plans of production wells may be optimized and monitored. Errors in production test data processing would decrease the production efficiency of the wells based on the production wells development design that have been made. In the other hand, the production can be analyzed each time based on available baseline condition prior to mitigate the problems. This paper describes the production analysis results based on case studies of XYZ-well production test at Dieng geothermal field that uses the Horizontal Lip Pressure Method with two silencers and two different types of weir boxes simultaneously at the same time of production test. Data processing is accompanied by data selection process, lip pressure correction, and correction of the gas presence. Total rate and total enthalpy of the fluid, which become the objectives of this production test, will be determined by Russel James equation and Hiriart equation. The results from both equations are also compared to know the sensitivity of each equation toward the available production test data. All results will be presented in graphs, so that reservoir characteristics and XYZwell condition can be specified based on those graphs. 1. INTRODUCTION Dieng is a water-dominated geothermal field located in Dieng Plateau, Central Java and developed at an altitude of about 0-2100 meters above sea level with temperatures between 10-20 C. Dieng field is made up of a series of volcanoes with andesite composition (Layman, 2). Geothermal manifestations that appear on the surface are in the form of solfatar, hot sping, mudpool, fumarole, and moffet. (Calibugan, et al., 0). Dieng field reservoir temperatures reach over 225 C and it is categorized as a hightemperature reservoir system. Based on the study results conducted by Layman (2), in general, Dieng geothermal field can be divided into two sectors, Sileri and Sikidang. Reservoir fluid of Sileri sector has moderate salinity, ph neutral, and low gas content characteristics, while fluid of Sikidang sector has high gas content with moderate feed zone depth and high enthalpy. Currently, Dieng geothermal field is managed by PT Geo Dipa Energy Dieng Unit (Persero). The field was developed with seven production wells and four injection wells for Unit 1, resulting in a production of 60 MW 2. THEORY 2.1 Horizontal Lip Pressure Method In this method, reservoir fluid will be ejected horizontally from the well through a lip pipe with certain diameter to the silencer. Measured parameters as input data for this method are wellhead pressure, lip pressure, and lip pipe diameter. Lip pressure is measured at the end of the lip pipe and fluid flow rate from silencer is determined using a weir box. 2.2 Weir Box Equation Weir box is used to determine the water mass flow rate from the silencer. In this case, there are two types of weir box that will be explained, V-Notch Weir Box and Rectangular Weir Box. 2.2.1 V-Notch Sharp Crested Weir Box Here are the equations used to determine water mass flow rate through the V-Notch Sharp Crested Weir Box. ( ) ( ) (1) (2) (3) 1
where Q, C, θ, h, and k are water flow rate (m 3 /hour), discharge coefficient, V-Notch angle (degree), weir height (m), and weir height correction (m), respectively. 2.2.2 Rectangular Sharp Crested Weir Box Figure 1: V-Notch Sharp Crested Weir Box. Slightly different from the equations before, these following equations are used for Sharp Crested Rectangular Weir Box. (4) where Q, C e, b e, h e, K, and g are water flow rate (m 3 /hour), discharge coefficient, space between weir plate (m), weir height (m), correction factor (m), and gravity acceleration (m/s 2 ), respectively. (5) (6) Figure 2. Rectangular Sharp Crested Weir Box. Discharge coefficient (C e ) and correction factor can be determined by following graphs. Figure 3. Graph to determine the value of discharge coefficient (C e ). 2
2.3 Russel James Equation Figure 4. Graph to determine the correction factor (K b ). Russel James equation connecting mass flow rate, enthalpy, lip pipe cross-sectional area, and lip pressure as follows. ( ) (7) (8) where W atm, H, P, A, h f, h g, h fg, and M are water mass flow rate (ton/hour), fluid enthalpy (kj/kg), lip pressure (bara), enthalpy of liquid-phase (kj/kg), enthalpy of vapor-phase (kj/kg), latent heat (kj/kg), and total mass flow rate (ton/hour), respectively. 2.4 Hiriart Equation Gerardo Hiriart creates an equation to determine the steam flow rate from vapor-dominated geothermal wells as follows. ( ) (9) This equation is quite simple, because it only requires pressure and lip pipe diameter as the variables. Hiriart equation negates the fluid enthalpy parameter, so that the accuracy level is lower than Russel James equation. Assumption applied in the Hiriart equation is steam flow to the atmosphere at sonic velocity condition as in perfect gas. However, Equation 6 can only be applied for vapor-dominated wells, so that to determine steam flow rate of the two-phase wells, Equation 6 needs to be corrected into the following equation. ( ) (10) where Q s, P, D, Q w, and Q s are vapor mass flow rate (ton/hour), lip pressure (psia), lip pipe diameter (inch), water mass flow rate at atmospheric pressure (ton/jam), and corrected vapor mass flow rate (ton/hour). 2.5 Lip Pressure Correction Equation Lip pressure obtained from measurement activity with a measuring device still needs to be corrected with a calibration factor based on lip pipe equation made by California Energy Company (CEC). ( ) (11) where P lip, P atm, P lip correction, and Corr. P lip are measured lip pressure (psig), atmospheric pressure (psia), calibration correction factor (inch), and corrected lip pressure (psia) 2.6 Gas Correction Equation Malcolm A. Grant develops an equation, which is a modification of the previous equation, to correct the results of enthalpy values due to the presence of gas in the well. The presence of gas that carried with the production fluid will give an impact on the results of enthalpy values. The effect can be very significant on certain conditions, so that the previous results need to be corrected by the following equation. ( ) ( ) ( ) (12) where ΔH, H, f lip, and H are enthalpy correction (kj/kg), initial enthalpy (kj/kg), gas mass fraction at lip pressure, and corrected enthalpy value (kj/kg). (13) 3
1240 1310 1335 1405 1435 1505 1550 1630 1730 1830 1930 2100 2300 100 300 500 700 900 1100 1300 1500 1700 0 2 0 2 120 300 500 700 900 1100 1300 1425 1500 1700 1900 2100 2300 100 300 500 700 900 1100 1330 1430 1700 1900 0 2 0 2 0 2 WHP (psig) Aris and Thasril 2.7 Two-Phase Fluid Enthalpy Equation Production test by using horizontal lip pressure method can also be conducted with more than one silencer for one well. If that condition exists, the enthalpy of the fluid can t be directly summed from the calculation of each silencer. Values that can be directly combined are the total fluid flow rate and water flow rate. By knowing the combination of both parameters from available silencers, enthalpy of the fluid can be determined from production test by the following equation. where H, x, W atm, M, h f, and h fg are fluid enthalpy (kj/kg), mass fraction of vapor, combined water flow rate (ton/hour), combined total fluid flow rate (ton/hour), enthalpy of liquid phase (kj/kg), and latent heat (kj/kg), respectively. 3. CASE STUDY Production test with horizontal lip pressure method, which performed on XYZ-well at Dieng geothermal field, used two lips and two different weir boxes. Lips and weir boxes specifications used are summarized in Table 1. Table 1: Lips and weir boxes specification used in XYZ-well production test. (14) (15) Static Data: Atmospheric Pressure: 11.5 psia ID of James Tube: 6.00 inch 8.00 inch Height from gauge to James Tube: 125.00 cm 145.00 cm Plip Correction: 49.21 inch 57.09 inch Left Right Rectangular Weir Box Type: V-Notch B = 1.20 m θ = 90 b = 0.34 m p = 0.73 m The production test was performed from March 24 th to 31 st, 1997. Production test profile during the time interval is shown in Figure 5. Production test measurements graph obtained shows inconsistent results in some particular time. Theoretically, the larger the valve is opened the wellhead pressure will be smaller and will be stabilized at a certain pressure during valve is opened at the same conditions. Therefore, production test data obtained need to be selected and validated more in accordance with the theory. Open valve 3.5" (10.7%) Open valve 5.25" (29.1%) WHP coming up Water turn to black Choke throtle valve due to sinker bar can not flow pass the fow tee Open valve 6.25" (40%) Open valve 7" (47.5%) 0 Open back to original position Weir on separator two overflow Open valve 7.25" (50.1%) Date and Time 24-03-1997 25-03-1997 26-03-1997 27-03-1997 28-03-1997 29-03-1997 30-03-1997 31-03-1997 Figure 5: Wellhead pressure profile of XYZ-well production test. With the valid production test data, subsequent data processing can be conducted that begins with the mass flow rate calculation process of each weir box. Mass flow rate of both weir boxes are next combined to obtain total mass flow rate, so that it can be used to calculate total enthalpy of the fluid produced from XYZ-well. Total mass flow rate and total enthalpy of the fluid are determined by two equations, Russel James equation and Hiriart equation, so the results of both equations can be compared. Production potential of XYZ-well and its reservoir conditions can be predicted based on the Wellhead Pressure vs. Total Rate and Wellhead Pressure vs. Enthalpy plot. Analysis of the two plots is explained in Grant, Malcolm A. and Bixley, Paul F.: Geothermal Reservoir Engineering Second Edition, Elsevier Inc., (2011), Oxford. Plots for XYZ-well are shown in Figure 6 and Figure 7. According to Wellhead Pressure vs. Total Rate plot, it can be analyzed that XYZ-well producing two-phase fluid and the result is appropriate with the actual condition. Then, according to Wellhead Pressure vs. Enthalpy plot, it can be analyzed that XYZ-well producing two-phase fluid from the reservoir which has a pretty good permeability. Total mass flow rate and total enthalpy of the fluid obtained from the Russel James equation and the Hiriart Equation have a little difference. It happens because the two equations have different calculation assumptions. Russel James equation has more complex variables, so that the Russel James equation has better accuracy rate, while Hirirat equation is simpler because it only includes lip pressure and lip pipe diameter as the variables. 4
Figure 6: XYZ-well flow test plot before data selection process. Figure 7: XYZ-well flow test plot after data selection process. The data processing results based on the Hiriart equation always deliver greater value compared with the calculation results based on the Russel James equation, both for total flow rate and total enthalpy of fluid. Statistically, the average percentage of difference value resulting from the Russel James equation and the Hiriart equation is 5.34% for the total flow rate value and 3.42% for the total enthalpy value of the fluid. Hiriart equation, which is simple, can be modified in order to generate value that approaching Russel James equation, so that the modified Hiriart equation further can be easily used and produce an accurate value. Modified Hiriart equation is formulated by adding a constant in the earlier Hiriart equation. The constant is determined by calculating the average value ratio of the results obtained from the equations of line representing Russel James equation and Hiriart equation. The equations are interpolation results of the equations of line corresponding to the graphs in Figure 7. The following equations are equations of line to determine the total flow rate of fluid. Russel James: Hiriart: (16) (17) The following equations are equations of line to determine the total enthalpy of fluid. Russel James: Hiriart: 5 (18) (19)
Table 2: Calculation results based on Russel James equation and Hiriart Equation (a) Total Rate; (b) Enthalpy Total Rate (ton/hr) Pwh Difference (psig) Russel Hiriart (%) James 202.05 211.96 4.91 210 201.81 211.73 4.91 220 201.53 211.45 4.92 230 201.22 211.14 4.93 240.87 210.79 4.94 250.48 210.40 4.95 260.06 209.97 4.95 270 199.59 209.50 4.96 280 199.09 208.99 4.97 290 198.56 208.44 4.98 300 197.98 207.86 4.99 310 197.37 207.23 5.00 320 196.72 206.57 5.01 330 196.03 205.87 5.02 340 195.31 205.13 5.03 350 194.54 204.35 5.04 360 193.74 203.53 5.05 370 192.91 202.67 5.06 380 192.03 201.77 5.07 390 191.12.84 5.08 190.17 199.86 5.10 410 189.18 198.85 5.11 420 188.16 197.80 5.12 430 187.10 196.71 5.14 440 186.00 195.58 5.15 450 184.86 194.41 5.16 460 183.69 193.20 5.18 470 182.48 191.95 5.19 480 181.23 190.67 5.21 490 179.94 189.34 5.23 500 178.62 187.98 5.24 510 177.26 186.58 5.26 520 175.86 185.14 5.28 530 174.42 183.66 5.30 540 172.95 182.14 5.32 550 171.44 180.58 5.34 560 169.89 178.99 5.36 570 168.30 177.35 5.38 580 166.68 175.68 5.40 590 165.02 173.97 5.42 163.32 172.21 5.45 610 161.59 170.42 5.47 620 159.81 168.59 5.49 630 158.00 166.73 5.52 640 156.15 164.82 5.55 650 154.27 162.87 5.58 660 152.35 160.89 5.61 670 150.39 158.86 5.64 680 148.39 156.80 5.67 690 146.35 154.70 5.70 700 144.28 152.56 5.74 710 142.17 150.38 5.77 720 140.03 148.16 5.81 730 137.84 145.91 5.85 740 135.62 143.61 5.89 750 133.36 141.27 5.94 760 131.06 138.90 5.98 770 128.73 136.49 6.03 780 126.36 134.04 6.08 790 123.95 131.55 6.13 121.50 129.02 6.19 Average 172.82 181.96 5.34 Conversion Factor 0.950 (a) 6 Enthalpy (kj/kg) Pwh Difference (psig) Russel Hiriart (%) James 1640.10 1688.10 2.93 210 1638.10 1686.24 2.94 220 1636.12 1684.40 2.95 230 1634.17 1682.59 2.96 240 1632.24 1680.81 2.98 250 1630.34 1679.05 2.99 260 1628.46 1677.33 3.00 270 1626.61 1675.63 3.01 280 1624.78 1673.97 3.03 290 1622.97 1672.33 3.04 300 1621.19 1670.72 3.06 310 1619.44 1669.14 3.07 320 1617.70 1667.59 3.08 330 1616.00 1666.06 3.10 340 1614.31 1664.57 3.11 350 1612.65 1663.10 3.13 360 1611.02 1661.67 3.14 370 1609.41 1660.26 3.16 380 1607.82 1658.88 3.18 390 1606.26 1657.53 3.19 1604.72 1656.20 3.21 410 1603.21 1654.91 3.22 420 1601.72 1653.64 3.24 430.26 1652.41 3.26 440 1598.82 1651.20 3.28 450 1597.40 1650.02 3.29 460 1596.01 1648.87 3.31 470 1594.64 1647.75 3.33 480 1593.30 1646.65 3.35 490 1591.98 1645.59 3.37 500 1590.69 1644.55 3.39 510 1589.42 1643.54 3.41 520 1588.17 1642.57 3.42 530 1586.95 1641.62 3.44 540 1585.76 1640.69 3.46 550 1584.59 1639.80 3.48 560 1583.44 1638.94 3.50 570 1582.32 1638.10 3.53 580 1581.22 1637.29 3.55 590 1580.14 1636.51 3.57 1579.09 1635.77 3.59 610 1578.07 1635.04 3.61 620 1577.07 1634.35 3.63 630 1576.09 1633.69 3.65 640 1575.14 1633.05 3.68 650 1574.21 1632.45 3.70 660 1573.30 1631.87 3.72 670 1572.42 1631.32 3.75 680 1571.57 1630.80 3.77 690 1570.74 1630.31 3.79 700 1569.93 1629.84 3.82 710 1569.15 1629.41 3.84 720 1568.39 1629.00 3.86 730 1567.66 1628.63 3.89 740 1566.95 1628.28 3.91 750 1566.27 1627.96 3.94 760 1565.61 1627.67 3.96 770 1564.97 1627.40 3.99 780 1564.36 1627.17 4.02 790 1563.77 1626.96 4.04 1563.21 1626.79 4.07 Average 1594.47 1648.99 3.42 Conversion Factor 0.967 (b)
Table 3: Calculation results based on Russel James equation and Modified Hiriart Equation (a) Total Rate; (b) Enthalpy Total Rate (ton/hr) Pwh Difference (psig) Russel Hiriart (%) James 202.05 201.31 0.36 210 201.81 201.09 0.36 220 201.53.83 0.35 230 201.22.53 0.34 240.87.20 0.33 250.48 199.83 0.33 260.06 199.42 0.32 270 199.59 198.97 0.31 280 199.09 198.49 0.30 290 198.56 197.97 0.29 300 197.98 197.41 0.29 310 197.37 196.82 0.28 320 196.72 196.19 0.27 330 196.03 195.52 0.26 340 195.31 194.82 0.25 350 194.54 194.08 0.24 360 193.74 193.30 0.23 370 192.91 192.49 0.22 380 192.03 191.64 0.21 390 191.12 190.75 0.19 190.17 189.82 0.18 410 189.18 188.86 0.17 420 188.16 187.86 0.16 430 187.10 186.82 0.15 440 186.00 185.75 0.13 450 184.86 184.64 0.12 460 183.69 183.49 0.11 470 182.48 182.31 0.09 480 181.23 181.09 0.08 490 179.94 179.83 0.06 500 178.62 178.54 0.04 510 177.26 177.21 0.03 520 175.86 175.84 0.01 530 174.42 174.43 0.01 540 172.95 172.99 0.02 550 171.44 171.51 0.04 560 169.89 170.00 0.06 570 168.30 168.44 0.08 580 166.68 166.85 0.10 590 165.02 165.23 0.13 163.32 163.56 0.15 610 161.59 161.86 0.17 620 159.81 160.12 0.19 630 158.00 158.35 0.22 640 156.15 156.54 0.25 650 154.27 154.69 0.27 660 152.35 152.80 0.30 670 150.39 150.88 0.33 680 148.39 148.92 0.36 690 146.35 146.93 0.39 700 144.28 144.89 0.42 710 142.17 142.83 0.46 720 140.03 140.72 0.49 730 137.84 138.57 0.53 740 135.62 136.39 0.57 750 133.36 134.18 0.61 760 131.06 131.92 0.66 770 128.73 129.63 0.70 780 126.36 127.30 0.75 790 123.95 124.94 0.80 121.50 122.54 0.85 Average 172.82 172.82 0.278 (a) Enthalpy (kj/kg) Pwh Difference (psig) Russel Hiriart (%) James 1640.10 1632.28 0.48 210 1638.10 1630.48 0.47 220 1636.12 1628.70 0.45 230 1634.17 1626.95 0.44 240 1632.24 1625.23 0.43 250 1630.34 1623.53 0.42 260 1628.46 1621.87 0.41 270 1626.61 1620.23 0.39 280 1624.78 1618.62 0.38 290 1622.97 1617.03 0.37 300 1621.19 1615.48 0.35 310 1619.44 1613.95 0.34 320 1617.70 1612.45 0.32 330 1616.00 1610.97 0.31 340 1614.31 1609.53 0.30 350 1612.65 1608.11 0.28 360 1611.02 1606.72 0.27 370 1609.41 1605.36 0.25 380 1607.82 1604.02 0.24 390 1606.26 1602.72 0.22 1604.72 1601.44 0.20 410 1603.21.19 0.19 420 1601.72 1598.96 0.17 430.26 1597.77 0.16 440 1598.82 1596.60 0.14 450 1597.40 1595.46 0.12 460 1596.01 1594.35 0.10 470 1594.64 1593.26 0.09 480 1593.30 1592.20 0.07 490 1591.98 1591.17 0.05 500 1590.69 1590.17 0.03 510 1589.42 1589.20 0.01 520 1588.17 1588.25 0.00 530 1586.95 1587.33 0.02 540 1585.76 1586.44 0.04 550 1584.59 1585.58 0.06 560 1583.44 1584.74 0.08 570 1582.32 1583.93 0.10 580 1581.22 1583.15 0.12 590 1580.14 1582.40 0.14 1579.09 1581.68 0.16 610 1578.07 1580.98 0.18 620 1577.07 1580.31 0.21 630 1576.09 1579.67 0.23 640 1575.14 1579.05 0.25 650 1574.21 1578.47 0.27 660 1573.30 1577.91 0.29 670 1572.42 1577.38 0.31 680 1571.57 1576.87 0.34 690 1570.74 1576.40 0.36 700 1569.93 1575.95 0.38 710 1569.15 1575.53 0.41 720 1568.39 1575.14 0.43 730 1567.66 1574.77 0.45 740 1566.95 1574.43 0.48 750 1566.27 1574.12 0.50 760 1565.61 1573.84 0.53 770 1564.97 1573.59 0.55 780 1564.36 1573.36 0.58 790 1563.77 1573.16 0.60 1563.21 1572.99 0.63 Average 1594.47 1594.47 0.281 (b) 7
The results obtained by the four equations above for some P wh values are shown in Table 2 (a) and (b). From the tables, it is known that the average value ratio obtained for the total flow rate of fluid is 0.950 with error percentage 5.34%, while for the total enthalpy of fluid is 0.967 with error percentage 3.42%. Based on these results, the constants used to determine the total flow rate and the total enthalpy of fluid by modified Hiriart equation are 0.950 and 0.967 respectively, so that modified Hiriart equation will be obtained as follows. ( ) (20) ( ) (21) Table 3 (a) and (b) show the errors of difference results between Russel James equation and modified Hiriart equation. According to the tables, it is known that the average error percentage of fluid total flow rate and fluid total enthalpy are reduced to 0.278% and 0.281% respectively. The great reduction in the average error shows that the results obtained from the modified Hiriart equation are very close to the results based on Russel James equation, as shown in Figure 8 and Figure 9. (a) Figure 8: Comparison of the total flow rate based on Russel James equation with (a) Hiriart equation; (b) modified Hiriart equation (b) (a) Figure 9: Comparison of the total enthalpy based on Russel James equation with (a) Hiriart equation; (b) modified Hiriart equation (b) 4. CONCLUSION The XYZ-well produces two-phase fluid from a geothermal reservoir that has a pretty good permeability. Total mass flow rate and total enthalpy value of the fluid based on Hiriart equation gives value that is always greater than the calculation results based on the Russel James equation. It is resulting average error percentage 5.34% for the total flow rate value and 3.42% for the total enthalpy value of the fluid. Modified Hiriart equation gives very close results according to Russel James equation with errors 0.278% for the total flow rate and 0.281% for the total enthalpy of the fluid. 8
REFERENCES Borromeo, CMR., and Orizonte, RG.: Output Computation for SNGP Wells Discharging Dry Steam, (1997). Cahyono, Yanuaris Dwi: The Application of Modified Hiriart for Fluids Flow Measurement in Geothermal Wells Flow Test Using Horizontal Lip Pressure Methode, (2013), Jakarta. Grant, Malcolm A. and Bixley, Paul F.: Geothermal Reservoir Engineering Second Edition, Elsevier Inc., (2011), Oxford. Hiriart, Gerardo: Steam Flow Rate Calculation by a Very Simple Equation, GRC6, 269-271, (1982). Saptadji, Nenny Miryani: Teknik Panas Bumi, ITB. 9