International Journal of Electrical and Computer Engineering (IJECE) Vol. No. February 2019, pp. ISSN: 2088-8708. DOI: 10. 11591/ijece. Combined heat and power-optimal power flow based on thermodynamic model with associated petroleum and wet gas utilization constraints Priambudi Pujihatma1. Sasongko Pramono Hadi2. Sarjiya3. Tri Agung Rohmat4 1,2,3Department of Electrical Engineering and Information Technology. Universitas Gadjah Mada. Indonesia 4Department of Mechanical & Industrial Engineering. Universitas Gadjah Mada. Indonesia Article Info ABSTRACT Article history: Oil fields produce associated petroleum and wet gas, which can be mixed with commercial natural gas as fuel. Associated petroleum and wet gas are a low cost, low quality fuel, whereas commercial natural gas is the opposite. Two parameters are affected by this mixture: the fuel cost and the power Ae steam output of gas turbine Ae heat recovery steam generators. This research develops a Unit Commitment and Optimal Power Flow model based on Mixed Integer Nonlinear Programming to optimize combined heat and power cost by considering the optimal mixture between associated petroleum - wet gas and commercial natural gas. A thermodynamic model is used to represent the performance of gas turbineAeheat recovery steam generators when subjected to different fuel mixtures. The results show that the proposed model can optimize cost by determining the most efficient power Ae steam dispatch and optimal fuel mixture. Furthermore, the optimization model can analyse the trade-off between power system losses, steam demand and associated - wet gas utilization. Received Feb 1, 2017 Revised Jul 31, 2018 Accepted Aug 13, 2018 Keywords: Associated petroleum gas Combined heat and power Mixed integer nonlinear Optimal power flow Unit commitment Copyright A 2019 Institute of Advanced Engineering and Science. All rights reserved. Corresponding Author: Sasongko Pramono Hadi. Department of Electrical Engineering and Information Technology. Universitas Gadjah Mada. Jl. Grafika No. Yogyakarta. Indonesia. Email: sasongko@ugm. INTRODUCTION Associated petroleum gas (APG) is natural gas that is produced as an oil field by-product. APG can be dissolved in oil or accumulated as a gas cap inside the oil reservoir . In addition to APG, oil fields produce wet gas from small size/marginal gas reserves. Due to the high gas processing cost and low gas quantity, wet gas and APG are used as fuel in oil field operations . , . In a conventional operation. APG used to be flared. Due to environmental constraints, oil companies must find a way to monetize APG and wet gas, with gas to wire (GTW) electricity generation as one promising alternative . APG quality is much lower compared to commercial natural gas. It has low heating content, high inert, high hydrogen sulfide and unstable volume and composition. The direct use of APG and wet gas as power generation fuel is limited by its low methane content and unstable gas composition. Wet gas contains higher hydrocarbons which can cause thermal regime failure in gas turbines . The low heat content of APG and wet gas may affect the gas turbineAos performance. A gas turbine requires a higher mass flow of APG and wet gas compared to natural gas to produce the same power output . Thus. APG and wet gas can be combined with commercial natural gas to minimize the effect on turbine parts and maintain reliable Several studies investigate the use of low quality gas for power generation. Panjaitan . conducted an experiment of running a generator with biogas fuel. The result showed that the generator output power has Journal homepage: http://iaescore. com/journals/index. php/IJECE Int J Elec & Comp Eng ISSN: 2088-8708 lower frequency when compared to high quality fuel. In terms of associated and wet gas utilization. Anosike . conducted a techno-economic study to evaluate APG utilization for power plants. Watanabe . proposed a gas turbine combined cycle system with high-voltage direct current (HVDC) transmission to monetize associated and stranded gas reserves. Arutyunov . proposed the conversion of APG to syngas for better performance. Gorbachev . 5 and 6 MW gas turbine generators fueled by APG. Guryanov . conducted a thermo-physical analysis of a combustion chamber and engine operation and made recommendations on conditions and restrictions on gas turbines running with APG. Vanadzina . proposed power generation as an effective solution for APG utilization in the reformed electricity market. Rajovic . conducted a life cycle assessment of APG utilization for combined heat power plants and heat boilers in an oil field. Although extensive research has been done, to the best of our knowledge, there has not been any study related to power & steam dispatch optimization that considers multiple fuel streams: APG, wet gas and commercial natural gas. A mature oil field utilizes an enhanced recovery method where the produced steam is injected into the ground to reduce oil viscosity and ease the oil flow to the surface. Combined heat and power (CHP) system produces steam by utilizing hot exhaust gas from the gas turbine. The thermal efficiency of this system is very high, up to 80 % compared to 36-50% of open cycle power generation . This scheme can be implemented by connecting heat recovery steam generators (HRSG) to a gas turbine exhaust. Some HRSGs are equipped with additional duct burners (DB) to increase steam production. In addition to the CHP system, steam is also produced from gas-fired boilers. In previous studies. CHP optimization is based on empirical objective functions and constraints. Roijers and Amerongen . model the objective function, fuel vs. power and steam, as a second order The thermodynamic process constraint is modeled by Guo . using feasible operating regions (FOR) which comprise several linear equations. Heuristic and deterministic methods have been used to solve CHP optimization. Some examples of heuristic algorithms to address the CHP problem are as follows: ant colony . , evolutionary programming . , . , particle swarm optimization . , . , self-adaptive genetic algorithm . , bee colony . , opposition-based group search . , differential evolution . and gravitational search algorithm . , . Some deterministic methods that have been implemented for CHP optimization are linear programming . Lagrangian relaxation . , and BenderAos decomposition . Shi . and Abdullah . included power system losses as constraints of CHP optimization. Kim . introduced a thermodynamic model in a combined heat and power economic dispatch (CHPED) in which objective function is to maximize revenue by exporting electricity in a deregulated market. These studies use a constant and predefined fuel price with a single fuel source. The effect of fuel mixture has not been discussed as yet and it is assumed that all gas turbines and HRSGs have a single fuel source of a constant This research focuses on combined heat and power (CHP) optimization with multiple fuel sources: commercial natural gas. APG and wet gas in a mature oil field operation. For convenience, the term Aufield gasAy will be used to represent APG and wet gas. The term Aunatural gasAy will be used to represent commercial natural gas. The contribution of this research is a Unit Commitment and Optimal Power Flow analysis with the following novelties: The use of thermodynamic model which considers multiple streams of fuel input: natural and field gas. For processes that operate over a wide range of operating condition, empirical models are unsuitable to represent the process dynamics . The empirical model on fuel vs. power vs. heat produced can no longer give a close approximation, since a field and natural gas mixture may have a large variance of gas compositions and heating content. Previous unit commitment problems have dealt with load demand, renewable generation, unit availability and energy prices as uncertainties . This optimization model will deal with new uncertainty factor, which is the fuel mixture between high-quality high-price natural gas with lowquality low-price field gas. The fuel cost will depend on optimum fuel mixture. A model that is able to analyze the tradeoffs between steam production, power system losses and field gas utilization. Steam production depends on power dispatch, which will affect the power system voltage profile and losses. Power dispatch depends on fuel quality, which is determined by mixture between natural and field gas. The case studies use a simplified CHP system in one of South East AsiaAos oil field. Combined heat and power-optimal power flow based on thermodynamicA (Priambudi Pujihatm. A ISSN: 2088-8708 RESEARCH METHOD System overview The system consists of 21 gas turbines, 5 HRSGs without duct burners and 3 HRSGs with duct The gas turbines supply an island power system with 450 MW average load, 13 substations, three voltage levels: 230, 115, 13. 8 kV, and 280 km length of transmission lines. Table 1 shows gas turbine capabilities, base load heat rate and steam production capabilities. Figure 1 shows the single line diagram. Table 1. Gas Turbine Generator Capability. Steam Production Capability and Heat Rate Unit NO. UNITS TYPE DG 1-3 Gas Turbine Type 1 AG1-5 Gas Turbine Type 2 Gas Turbine Type 2 MG1-5 Gas Turbine Type 3 MG6-8 Gas Turbine Type 2 MG9-11 Gas Turbine Type 4 Gas Boiler NG = Natural Gas FG = Field Gas ACTIVE POWER (MW) Reactive Power (MVAR) Base load Heat rate (GJ/MWH) Steam Prod. 3 cwep. 15,900 1,590 Fuel NG FG NG FG NG FG NG FG Figure 1. Simplified single line diagram There are four gas turbine types in the system: Types 1-4, with different capabilities and base load heat rate. The base load heat rate indicates gas turbine efficiency. The most efficient gas turbine is DG1-3 (Gas Turbine Type . , while the least efficient is MG1-5 (Gas Turbine Type . The gas turbines are divided into six groups. each group comprises identical gas turbines. For example. MG 1-5 consists of five identical Type 3 gas turbines. Gas turbines AG 1-5 and DG 1-3 have the capability to produce steam. DG 1-3 have additional duct burners installed to produce more steam. Ten gas-fired boilers (GB) are also utilized to produce steam. Gas turbines AG1-5. MG1-5. MG6-8. MG9-11 have two fuel sources: natural gas and field gas. Other gas turbines and gas boilers only have natural gas as single fuel source. Each field gas stream has an average production volume, shown in Table 2. Total field gas production volume is 37,357 cubic meters per hour at standard condition . oC temperature, 1 atm pressur. Table 2 shows significant difference between natural and field gas in terms of heating value. Table 2. Natural and Field Gas Data towards AG1-5. MG1-5. MG6-11 Parameter Unit Density . Volume . Low Heating Value . kg/m3 m3/hour kJ/kg Natural Gas 44,821 Int J Elec & Comp Eng. Vol. No. February 2019 : 42 - 54 Field Gas AG1-5 15,612 32,394 Field Gas MG1-5 4,130 30,932 Field Gas MG6-11 17,615 20,466 Int J Elec & Comp Eng ISSN: 2088-8708 Gas and thermodynamic analysis Firstly, the optimization model analyzes the composition of each gas stream . atural gas and field The net heating value per mass unit is calculated as follows: yayaycO , yayaycO , . yayaycO . Subscript NG and FG correspond to natural gas and field gas, respectively. Mass flow for natural and field gas is the optimizationAos control variable. Then, the LHV value becomes the input for GT-HRSG thermodynamic equations. GT-HRSG thermodynamic equations are based on the model proposed by Kim . The model calculates the final output of the gas turbines and HRSG in terms of power and steam. Gas turbine model The gas turbine is modeled using the following formulas. Air mass flow entering compressor . ycO = Air volumetric flow . ycO =ycO . Compressor outlet temperature (K): ycN =ycN 1 Turbine inlet temperature (K): yc, =yc yc The gas composition and mass flow play an important role in determining LHVm,mix value. Turbine outlet temperature (K) ycN =ycN 1Oe 1Oe ycN =ycN . ycu = ycEycI ycu = ycEycI Compressor and turbine efficiency is assumed constant throughout the operating region. Gas turbine power output (MW): ycO yc , ya ycN Oe ycN Oe ycO ya . cN Oe ycN ) AE1000 Fuel consumption can be converted to volume . tandard meter cubic per secon. as given: ycO =yc AEyuU ycO =yc . AEyuU HRSG and duct burner model HRSG and the duct burner are modelled with these formulas: HRSG inlet temperature (K): ycN =ycN yu =yuC ya . yu =yuC yuC yayaycO , . for HRSGs which are not equipped with duct burners, the value of wf,DB=0. Combined heat and power-optimal power flow based on thermodynamicA (Priambudi Pujihatm. A HRSG steam production . ISSN: 2088-8708 HRSG steam production can be converted to cold water equivalent meter cubic per second by using cold water density value at 1000 kg/m3: AE1000 = ycO Gas-fired boiler model The gas-fired boiler is modelled as follows: Gas boiler steam production . ycO yu =yuC . yayaycO . Gas boiler steam production converted to cold water equivalent meter cubic per second: = ycO AE1000 Optimization formulation In general, optimization is divided into two stages: Unit Commitment: to determine which gas turbine and gas-fired boiler units should be on/off based on predicted load and steam demand. The gas turbine. HRSG and gas boiler on / off status becomes input for the next optimization stage (Optimal Power Flo. Optimal Power Flow: to calculate the most optimal output of the gas turbines. HRSG and gas-fired boilers while still meeting power system constraints, such as generator active and reactive power limits, transformer taps and bus voltages. Optimal power flow analysis is needed to understand the tradeoff between power system losses, voltage and power-steam dispatch. Objective functions The objective function is the minimum cost for power and steam generation, while still meeting power system constraints. The total fuel consumption for gas turbines, duct burners and gas-fired boilers is formulated as follows: Total fuel for gas turbine ($): ya yc , [ yc , yayaycO , , ycyyc , yayaycO , , ycyyc where NGT is the total number of gas turbines. Total fuel for duct burners ($): ya =Oc =Oc yc yc, , yayaycO , ycyyc where NDB is the total number of duct burners. Total fuel for gas-fired boilers ($): ya =Oc yc yc, , yayaycO , ycyyc where NGB is the total number of gas boilers. The objective function can be formulated as: min ya ya ya . The sGT,i, sDB,i and sGB,i variables are binary numbers, which indicate the running status for gas turbines, duct burners and gas boilers. Thus, the optimization is a mixed integer non-linear programming After securing the status of the gas turbines, duct burners and gas boilers, optimization is refined using the Combined Heat and Power Optimal Power Flow model. The objective function is still the same Int J Elec & Comp Eng. Vol. No. February 2019 : 42 - 54 Int J Elec & Comp Eng ISSN: 2088-8708 with equation . , only the binary variables have been attained from 1st stage optimization. Natural and field gas prices are assumed to be 5$/GJ and 1$/GJ, respectively. Constraints The constraints for 1st stage optimization are: Spinning reserve OOc OuycO OycE OycE OycO OycO. OycO. Turbine inlet temperature ycN, ycO Steam production capability limit Field gas availability limit =ycE where NHRSG is the total number of HRSG. Real power capability limit Steam demand and supply balance Oc ycIycI Real power balance Oc OycN. OycN. Fuel mass flow limit yc, . Oyc . Oyc . The above . constraints are still maintained in the 2nd stage of optimization, with the following additional constraints: Power Flow Balance OeOc OO ycE OeOc OO ycE = ycO ya ycO ycO ya ycaycuyc yuE Oe yuE yaA ycycnycu yuE Oe yuE = OeycO yaA ycO ycO ya ycycnycu yuE Oe yuE Oe yaA ycaycuyc yuE Oe yuE OeycE OeycE Reactive power capability limit OycE OycE Combined heat and power-optimal power flow based on thermodynamicA (Priambudi Pujihatm. A Voltage levels ycO ISSN: 2088-8708 OycO OycO Limit of transformer tap changer ycN . OycN . OycN . For unit commitment analysis, the model uses the BONMIN solver . which is provided in OPTI . , a MATLAB-based optimization add-in. For the Combined Heat and Power Optimal Power Flow, the model uses three solvers: Interior Point (IPOPT) . provided in OPTI. Active Set (MATLAB Optimization Toolbo. and Sequential Quadratic Programming, provided in TOMLAB. These solvers are used for comparison purposes. Case studies are simulated based on steam and electrical load. Steam demand for high, medium and low categories are 2,319. 1,656 and 662 m3 cold water equivalent per hour, respectively. Electrical load for high and low categories are 450 and 350 MW. Thus, the simulation covers six scenarios, as follows: Scenario 1: high steam demand and high electrical load Scenario 2: medium steam demand, high electrical load Scenario 3: low steam demand, high electrical load Scenario 4: high steam demand and low electrical load Scenario 5: medium steam demand, low electrical load Scenario 6: low steam demand, low electrical load RESULTS AND ANALYSIS First, we discuss the unit commitment result shown in Table 3, which calculates the number of gas turbines running. DG-Type 1 gas turbines have the lowest heat rate. thus they are the most efficient. Furthermore, these units can also deliver steam. Hence, it is recommended that all three DG Type 1 gas turbines run in every scenario. Table 3. Number of Gas Turbine Running Scenario DG - Gas Turbine Type 1 AG Ae Gas Turbine Type 2 MG- Gas Turbine Type 3 MG- Gas Turbine Type 2 MG Ae Gas Turbine Type 4 Total Unit High SteamAe High Load Med. SteamAe High Load Low SteamAe High Load High SteamAe Low Load Med SteamAe Low Load Low SteamAe Low Load AG Type 2 gas turbines are always prioritized to run due to their steam production capabilities. AG Type 2 is only less prioritized when the steam demand is low. MG Type 4 gas turbines are the second most efficient. During high load scenarios, it is recommended that all three run. In low load scenarios. MG Type 4 priority becomes lower/similar with AG Type 2 units, which have higher heat rate but are capable of producing steam. Only one out of five MG-Type 3 units that is recommended to run in most scenarios. Even though MG Type 3 is the least efficient unit, it still has higher prioritization compared to the other gas turbines which have lower heat rates. This prioritization is due to the availability of dedicated field gas for MG-Type 3, as shown in Table 2. Further simulation shows that the power dispatch for MG Type 3 is equal to its available field gas supply. Figure 2 shows the real power dispatch. DG Type 1 units are running on base load in all scenarios due to their low heat rate. AG Type 2 can produce steam and their power dispatch depends on the steam MG units, which do not have a steam production capability, only provide an additional power supply to meet electrical load demand. All MG unitsAo power dispatch increases significantly in low steam demand scenarios. Table 4 shows that there is a significant change in power system losses in low steam demand scenario. As shown in Appendix Tables 10 and Table 11, the majority of the load is in F. M substations, which are located closer to the MG units. When steam demand is low, there is no need to generate much power from Int J Elec & Comp Eng. Vol. No. February 2019 : 42 - 54 Int J Elec & Comp Eng ISSN: 2088-8708 AG. Thus, the optimization model would rather prioritize MG units, which are located close to the load center, to significantly reduce power system losses. DG-Type1 Real Power (MW) AG-Type2 MG-Type3 MG-Type2 MG-Type4 CG-Type2 High Steam High Load Med Steam High Load Low Steam High Load High Steam Low Load Med Steam Low Load Low Steam Low Load Figure 2. Real power dispatch (MW) Table 4. Power System Losses (MW) Scenario High Steam Ae High Load Med. Steam Ae High Load Low Steam Ae High Load High Steam Ae Low Load Med Steam Ae Low Load Low Steam Ae Low Load Losses Table 5 and Table 6 show the optimal fuel composition in terms of volumetric flow and mixture ratio, respectively. As shown in Table 2, the available field gas production from AG units is 15,612 standard m3 per hour. Table 5 shows that AG units utilize all available field gas in all scenarios. Steam production from AG units depends on their power dispatch. Natural gas acts as a Auback upAy fuel to provide additional power output to meet steam demand. During low steam demand, natural gas consumption decreases Table 5. Fuel Volumetric Flow (Standard m3 per Hou. AG Type 2 High Steam Ae High Load 30,813 15,612 Med Steam Ae High Load 30,813 15,612 Low Steam Ae High Load 1,250 15,612 High Steam Ae Low Load 6,810 15,612 Med Steam Ae Low Load 1,044 15,612 Low Steam Ae Low Load 1,044 15,612 NG= Natural Gas FG = Field Gas Scenario MG Type 3 4,130 4,130 2,831 4,130 4,130 4,130 MG Type 2&4 7,671 17,615 7,671 17,615 30,181 16,493 13,458 15,265 15,265 Table 6. Optimal Mixture of Consumed Natural Gas: Field Gas by Volume Scenario AG Type2 MG Type 3 MG Type 2&4 High Steam Ae High Load 66%:34% 0%:100% 30%:70% Med. Steam Ae High Load 66%:34% 0%:100% 30%:70% Low Steam Ae High Load 7%:93% 41%:59% 65%:35% High Steam Ae Low Load 30%:70% 0%:0% 0%:100% Med Steam Ae Low Load 6%:94% 0%:100% 0%:100% Low Steam Ae Low Load 6%:94% 0%:100% 0%:100% Combined heat and power-optimal power flow based on thermodynamicA (Priambudi Pujihatm. A ISSN: 2088-8708 MG Type 3 units are the most inefficient unit based on heat rate value. In most scenarios, power dispatch on MG-Type 3 is recommended at 5. 7 MW. The available field gas for MG Type 3 based on Table 2 is 4,130 standard m3 per hour. Table 5 and Table 6 show that in most scenarios. MG-Type 3 consumes only the available field gas production, without consuming additional natural gas. Thus, this unit is running merely to consume field gas in the majority of scenarios. This result shows that apart from heat rate, low cost field gas availability also affects power dispatch optimization. Inefficient units can be dispatched at the exact amount of equivalent available field gas energy. A different condition is found in the low steam-high load scenario. The optimization model recommends MG Type 3 consumes natural gas, with the purpose of boosting MG Type 3 output. As stated in the previous section, boosting the power output from MG units, including Type 3, will have a significant impact on loss reduction. In the high steam-low load scenario, it is not recommended that MG Type 3 to run. In this scenario, units that are capable to produce steam are prioritized. With high heat rate values. Type 3 is the unit given the lowest priority. MG Type 2 and Type 4 share the same field gas supply, as shown in Table 2. In the first two scenarios . igh steamAehigh load and med steam-high loa. , field gas becomes the priority fuel. However, in low steam-high load scenarios, the optimization model chooses natural gas instead of using all the available field gas. Since natural gas has a better LHV value, selecting natural gas boosts the output of MG Type 4 and Type 2. Since MG units are closer to the load center, boosting the power output of MG Type 4 and 2 reduces power system losses significantly from 9. 47 to 5. 78 MW, as shown in Table 4. In all low loads scenarios. MG Type 2 and 4 only consume field gas. However, the consumed volume is less than the maximum available volume stated in Table 2. During low load scenarios, units that are capable of producing steam still have significant power dispatch to boost steam production. The remaining load that is supplied from MG Type 2 and 4 is not enough to consume all available field gas. Table 7 summarizes the percentage of total utilized field gas compared to available field gas production. Of the six scenarios, it is only in the high steam-high load and medium steamAehigh load scenarios that all the available field gas production is utilized as fuel. Table 7. Percentage of Field Gas Utilization by Volume in Each Scenario High Steam Ae High Load 37,357 37,357 Scenario Utilized Field Gas m3/hr @ STD Available Field Gas m3/hr @ STD Field Gas % Utilization Med. Steam Ae High Load 37,357 37,357 Low Steam Ae High Load 36,235 37,357 High Steam Ae Low Load 29,070 37,357 Med Steam Ae Low Load 35,006 37,357 Low Steam Ae Low Load 35,006 37,357 Figure 3 shows that steam demand has a significant impact on power system voltage quality. To boost steam production, power dispatch on DG and AG units must be increased. These units are located remotely from the load center. Thus, during high and medium steam demand, power system losses tend to be higher compared to the low steam demand scenario. This causes the voltage profile in high and medium steam demand scenarios to be lower compared to low demand scenario. High Steam - High Load Low Steam - High Load Med Steam - Low Load Voltage Profile . Med Steam - High Load High Steam - Low Load Low Steam - Low Load Substation Figure 3. Voltage profile . er uni. Int J Elec & Comp Eng. Vol. No. February 2019 : 42 - 54 Int J Elec & Comp Eng ISSN: 2088-8708 Steam production from HRSG. HRSG-DB and GB is shown in Table 8. It is recommended that GBs run only in high steam demand scenarios. Table 9 shows the objective function in each scenario. All three solvers produce the same optimal cost in all scenarios. Table 8. Steam Output . 3/hou. Scenario High Steam Ae High Load 1,829 Med. Steam Ae High Load 1,346 Low Steam Ae High Load High Steam Ae Low Load 1,829 Med Steam Ae Low Load 1,553 Low Steam Ae Low Load Med Steam Ae Low Load 25,410 25,410 25,410 Low Steam Ae Low Load 17,050 17,050 17,050 Table 9. Optimal Cost ($/hou. Scenario OPTI Matlab Tomlab High Steam Ae High Load 37,075 37,075 37,075 Med. Steam Ae High Load 29,781 29,781 29,781 Low Steam Ae High Load 22,935 22,935 22,935 High Steam Ae Low Load 33,092 33,092 33,092 CONCLUSION This paper presented an optimization model based on mixed integer nonlinear programming which can analyze the tradeoff between cost, steam production, losses and field gas utilization. The utilization of multiple fuel streams is a unique problem, the fuel mixture must be optimized to achieve the lowest fuel cost. The unit commitment result is not solely dependent on heat rate. rather, it also depends on field gas A less efficient unit can be prioritized to run if there is a sufficient amount of field gas supply. There is a strong relationship between power dispatch, steam dispatch, field gas utilization and power system losses. With high steam demand, high losses are accepted as a consequence of maximizing steam production from gas turbines with HRSG. Steam demand also has a significant impact on power system voltage quality. Since the majority of the load is far from the steam-producing gas turbines, the voltage profile in high steam demand scenarios is lower compared to low steam demand scenarios. The simulation is able to find the most optimal fuel mixture. The optimal mixture does not always mean that all available low-cost field gas must be consumed. This demonstrates that there is a need to limit field gas production to just the right amount to avoid over production and ensure optimal cost. The optimization model is also able to give the most optimal steam dispatch to HRSG. HRSG-DB and GB. The simulation shows that GB is the least prioritized steam producer. The case studies in this work only focus on single objective optimization. For future work, this can be expanded to multi-objective optimization, with cost vs. loss, cost vs. emission or cost vs. system voltage stability as the objective functions. Regarding gas, this research only focuses on gas volume. In future work, this can be expanded to handle constraints in the gas distribution network. REFERENCES