Unu GTP 2019 09 PDF
Unu GTP 2019 09 PDF
Unu GTP 2019 09 PDF
Jeffrey M. Andal
Department of Energy
Energy Center, Rizal Drive
Bonifacio Global City, Taguig City 1632
PHILIPPINES
andal_jef@yahoo.com.ph
ABSTRACT
The growth of the geothermal industry in the Philippines remains slow despite the
initiatives of the government. Based on the study made by the Philippines’
Department of Energy, one of the major factors affecting the growth is the high cost
of development for low-to-intermediate temperature geothermal resources. To
provide a deeper understanding of this cost, we performed power plant modelling
and thermoeconomics modelling for a prospective geothermal power plant in
Montelago, Philippines. The models in this report were used to determine the
following: the optimal power output of the system, the cost of the power plant
development, the minimum cost of electricity generation, and exergetic evaluation.
The model indicates that the cost of installed power to develop the Montelago
geothermal prospect is 3574 USD/kWe and the minimum cost of electricity
generation and the break-even point is 6.22 USDcents/kWh and 8.26
USDcents/kWh, respectively. This suggests that further development of the prospect
could be economical.
1. INTRODUCTION
Utilization of geothermal energy for power generation plays a very important role in the Philippines.
According to the 2018 Power Statistics Report of the Philippines’ Department of Energy (DOE, 2019a),
geothermal energy provides a total of 10 270 GWh. This is about 10.88% of the total gross generation
of the country. So, for every ten lightbulbs, it may be assumed that one lightbulb is powered by
geothermal energy. However, despite of this important role in meeting the energy demand of the
country, the rate of development of this resource seems to be slowing down. As a response, the
Philippines’ government continues to promote its utilization.
The country has been exploring geothermal energy for more than fifty years. Its development started
on 1967 when the first light bulb using electricity from geothermal, was lit in Tiwi, Albay. Now, the
country is the third (3rd) largest power producer in the world using geothermal energy, after the USA
and Indonesia, with a total installed capacity of 1918.16 MW (DOE, 2019b). Based on the internal
report of the DOE (2019b), since the start of exploitation of geothermal energy for power production in
1977 until 2018, the country has already generated almost 300 000 GWh of electricity using geothermal
83
Andal 84 Report 9
providing the country with a clean, reliable, and affordable source of electricity. Based on the DOE´s
report the geothermal energy development in the Philippines is illustrated in Figure 1.
Figure 1 shows that there have been two leaps of development of geothermal energy utilization. These
leaps of development have mainly been driven by governmental intervention. The first intervention was
when the government addressed the oil-crisis in the mid-1970’s by developing indigenous resources
such as geothermal energy. This led to the development of Tiwi, Makban, Tongonan, and Palinpinon
geothermal power plants. The other spurt of development was in 1990’s where the government enacted
a law, allowing the Build-Operate-Transfer scheme. This led to the development of the Mt. Apo
geothermal field and the expansion of the Tongonan geothermal field (Clemente et al., 2016).
Additionally, it can be seen in Figure 1 that the intervention to further develop the utilization of
geothermal energy does not stop in 1990’s. In 2008, the Philippines’ government enacted the Renewable
Energy Act, which was intended to promote the development, utilization and commercialization of
renewable energy resources in the country. However, despite of commercialization of new power plants
and expansion of existing geothermal fields, the development of utilization of geothermal energy
remains slow.
According to the DOE (2019b) internal report, the following are the barriers that affects the development
of geothermal energy in the Philippines: (1) environmental and socio-cultural concerns in protected
areas and ancestral lands; (2) high cost of development for low-to-intermediate temperature and acidic
geothermal resources; (3) low-level awareness on non-power applications of geothermal energy; and (4)
numerous permits and time-consuming processes to get permits. The identified barriers prompted the
government to issue policies and projects to address it. Some of the policies and projects are shown in
Table 1 below:
The policies and projects listed aim to assist the investors in tackling the barriers in the development of
a geothermal resource. The main barrier is the high cost of development for low-to-intermediate
temperature and acidic geothermal resources. While there are few known acidic resources in the
Philippines, there are about 428 MW known low-to-intermediate temperature geothermal resources in
the country and this number might continue to increase as further studies in several fields are being
conducted (DOE, 2019b). Utilizing this low-to-intermediate temperature geothermal resources would
be beneficial to the country as these resources are clean and indigenous.
Report 9 85 Andal
Policies Projects
Executive Order No. 30 (Approved June 28, 2017): Detailed assessment of selected low-
Creating the Energy Investment Coordinating Council enthalpy geothermal resources in the
(EICC) in order to streamline the regulatory procedures Philippines (2010 – 2015).
affecting energy projects.
Republic Act No. 11032 (Approved May 28, 2018): Comprehensive resource assessment of
An act promoting Ease of Doing Business and efficient Philippine low-enthalpy geothermal areas
delivery of government services. (2015 – 2017).
Republic Act No. 11234 (Approved March 08, 2019): Philippine geothermal resource inventory
An act establishing the Energy Virtual one-stop shop for and assessment (ongoing).
streamlining the permitting process of power generation,
transmission, and distribution projects.
Among the areas with low-to-intermediate temperature geothermal resources, one of the most advanced
fields in the exploration stage is the Montelago geothermal prospect located in Mindoro Island,
Philippines (Figure 2). Thus, this report uses the Montelago geothermal prospect as a case study to
provide a better understanding of the cost of development of this type of resource using power plant and
thermoeconomics modelling. However, this report will not cover a comparison of development using
emerging technologies and this report will not include the barriers that affect the development of a
project.
A resource with a temperature of less than 150°C is considered a low-temperature resource. These
resources can be used for heat pumps or direct use application. Resources with temperatures ranging
from 150 to 200°C are moderate or intermediate temperature resources. These intermediate resources
can be used for electricity generation, but the wells might require pumping as these wells are not capable
of producing large quantities of fluids which are required for large scale electric production (DiPippo,
2016).
According to the DOE´s internal report (2019a), it is estimated that the Philippines has a geothermal
energy potential of about 4407 MW. As shown on Appendix I,, the known estimated potential from the
low-to-intermediate temperature geothermal resources is about 428 MW and some fields requires further
study to determine their potential.
The Montelago geothermal prospect is in the north western part of Mindoro Island, Philippines between
the Naujan lake and Tablas strait. The island of Mindoro is the seventh largest island in the Philippines
with a total area of 10 571 km2. It has a total population of 1 331 473 (2015). The economy on this
island is largely based on agriculture and tourism (Wikepedia.org, 2019).
The prospect has a tropical climate with significant rainfall in Naujan that averages to 2003 mm per year
and a short dry season. The average annual temperature in the prospect is 27.3°C (Climate-Data.org,
2019).
For 2016, Mindoro island has e total peak demand of more than 64 MW based on the combined total
load of the electric cooperatives located in Mindoro (NGCP, 2019). The cost of electricity in Mindoro
remains expensive as some part of the island relies on generators fuelled by diesel. According to Ahmed
(2019), the true cost of diesel in Mindoro is ranging from 9.60 to 38.47 PHP/kWh or 18 to 74
USDcents/kWh, based on the foreign exchange rate of 1 USD to PHP 52.32 dated 22 August 2019 (XE
Corporation, 2019), while the existing subsidized approved generation rate is 11 USDcents/kWh (5.64
PHP/kWh).
Several exploration studies had been carried out since 1979 until present. However, these studies made
have inconsistencies. To get a good estimate of the potential of the geothermal prospect, this section
will describe and adapt information from the recent resource assessment made on the prospect.
However, this is not intended to provide a complete review and or to question other studies made on the
prospect, but instead is a summary of details that are important for the elaboration of the power plant
modelling. Additionally, this report is independent of other previous studies and not meant by any
means to criticize the other studies. This report is purely aimed at the evaluation of the resource using
modelling.
The prospect has low to moderate elevation and is located between the Naujan lake and Tablas strait. It
is bounded by 19 fault lines (Figure 3) and has four eruptive centres, which are Mt. Montelago, Mt.
Pungao, Mt. Buloc, and Mt. Matabang Bundok (Asmin et al., 2016).
Report 9 87 Andal
FIGURE 3: Fault system in the Montelago geothermal prospect (Asmin, et al., 2016)
At the beginning of the exploration, in 1978-1979, eight shallow temperature gradient holes with depths
ranging from 195 to 305 m were drilled. Two of them (NGH-4 and NGH-6) are near the prospect. NGH-
4 and NGH-6 have thermal gradients of 23°C / 100 m and 10°C / 100 m, respectively. The temperature
obtained in NGH-4 and the chemical composition of the Pungao springs encouraged further exploration
of the prospect (SKM, 2011). From 2014 to 2016, four new wells were drilled, that is two slim holes
and two deep wells. Table 2 shows the information about these wells:
The wells drilled in the prospect provided insights into the resource conditions. Among these wells,
slim hole SH-02 was the only well that was able to produce. This slim hole showed an elevated
concentration of NCG (non-condensable gas) but is remains unclear whether the NCG enter the wellbore
at shallower depth (which could be cased-off when drilling a deeper well) or if a gas rich zone is present
at greater depth. Further, the well also encountered significant calcite scaling. The formation of
Andal 88 Report 9
amorphous silica scale is unlikely at temperatures above approximately 110°C, at lower temperatures
silica oversaturation occurs. However, it is believed that the fluid produced during the test is not
representative of the deeper reservoir fluid (GeothermEx, 2017). On the other hand, wells MN-01 and
MN-02 appears to be the out of the reservoir zone while SH-01 confirmed the intermediate temperature
of the reservoir but has too low permeability to produce.
Based on the additional information obtained from the wells, GeothermEx´s (2017) report provided a
conservative estimate of reservoir area of about 0.8 km2 to 2.4 km2. This estimate considers areas that
might be acidic or too low in permeability and nonetheless implies that the reservoir is economically
exploitable. Results from wells MN-01 and MN-02 show marginally economic temperatures for self-
flowing wells that are too hot to pump. Additionally, based on the depths with considerable permeability
and to avoid the excess risk of a long open-hole section, the report estimates the average reservoir
thickness to be between 1300 m and 1700 m. The report also estimated a range of temperatures which
lies between 175°C and 205°C. Other than these, the report also described the following parameters to
estimate the recoverable energy of the resource as seen in Table 3:
With the information above, GeothermEx estimated that the capacity of the resource is around 15 MW
and could potentially be up to 24 MW.
For this paper, the author has adopted most of the estimates made by GeothermEx (2017) while some
information was changed to accommodate the design of a binary geothermal power plant that has a plant
life of 30 years. Additionally, to simplify the calculation, the most likely reservoir temperature was
adopted from Axelsson and Halldórsdóttir (2015) which is 190°C. It was also assumed that the
thermodynamic properties of the geothermal fluid were pure water properties. To estimate the
reinjection temperature the author applies the silica “rule of thumb” which entail that it is only possible
to cool the water by around 100°C without risking scaling (Thórhallsson, 2005). Thus, the author uses
Report 9 89 Andal
90°C as the reinjection temperature in this study. This report also assumes that the conversion efficiency
of the system is based on the typical thermal efficiency of binary plants which is 8 to 12% (DiPippo,
2016). The parameters are summarized in Table 4.
This section aims to discuss the cycle of the system and to determine the optimal capacity of a generating
unit that could be built with the described system. Additionally, in this section we will determine the
cost of developing the prospect. The characters, symbols and abbreviations used in this model are
defined in the nomenclature at the end of this report.
Andal 90 Report 9
Most binary geothermal systems use an Organic Rankine Cycle (ORC). This system is analogous to a
thermal power plant that uses a Rankine Cycle, except that the thermal energy used for the binary
geothermal system came from a geothermal reservoir instead of a boiler in Rankine Cycles (Fujii et al.
2011). There are many variations of this system and which the best arrangement is determined mainly
by the constraints of the resource and the environment.
One of the constraints of the system is the reinjection temperature. According to DiPippo (2016), when
reinjection temperature is limited, heat recuperation can increase power output resulting in
improvements in the cycle thermal efficiency. Since this report considers a reinjection temperature limit
of 90°C, a recuperator was considered in the model.
The prospect is situated between the Naujan lake and the Tablas strait. This does seem to be an ideal
location to use an evaporative heat reinjection system due to presence of water near the prospect.
However, using this type of condenser is not ideal as the climate in the area is humid and the dry period
is very short. Therefore, air-cooled condensers (ACC) are an ideal choice for the cooling system.
In this model, the heat extracted from the fluid (Step s-1) is transferred through the heat exchangers to
preheat (Step 2 - 3) and vaporize (Step 3 - 4) the pressurized secondary fluid or working fluid circulated
in a closed loop, while the geothermal fluid is reinjected into a reinjection well (Step s-3). The vaporized
secondary fluid is expanded in a turbine coupled to a generator to generate electricity (Step 4 - 5). The
expanded vapour then passes through a recuperator (Step 5 - 6) where the unused heat is transferred to
the working fluid that is in between the preheater and the feed pump (Step 1 - 2). The working fluid
vapour is then condensed in the condenser (Step 6 - 7). Then the working fluid is pumped at high
pressure (Step 7 - 1) to be reheated, thus completing the cycle.
Report 9 91 Andal
In the schematic diagram of the model (Figure 4) several major components were identified. This
section will discuss and describe the other calculations made.
As mentioned in Section 2.3 of this report, in the Montelago prospect there are eight shallow gradients,
two slim holes and two deep wells. However, none of these well is a good candidate for power
production regarding the design objectives in this report. Thus, with the parameters from the resource
assessment, Table 5 summarizes the main parameters of the four producing wells that will be used in
this report.
below the well head and the required discharge pressure is 3.5 bar.
With those assumptions made, a line-shaft pump with a motor size of 200 kW must be used to give a
flow rate of 40 kg/s. Additionally, it was assumed that the pump at 300 m depth will be kept the well
running during the desired plant life (Lýdur Skúlason, project manager at Deilir Technical Service, pers.
comm., 2019).
Heat exchangers are devices used to transfer heat from a hot fluid to a colder fluid. In the system shown
in Figure 5, there are four heat exchangers (evaporator, preheater, recuperator and condenser). In the
preheater and evaporator, the heat from the geothermal fluid is used to heat-up and vaporize the working
fluid. In the recuperator, the heat from the exhaust from the turbine is used to heat-up the working fluid
before the preheater inlet. In the condenser, the heat from the working fluid is removed and condensed
by cold air. The heat transfer in the heat exchangers can be calculated based on the energy balance of
hot and cold fluid as listed in Table 6.
The turbine uses the enthalpy from the pressurized working fluid vapour and converts it into mechanical
energy that will be converted into electrical energy in the generator. The power generated in the turbine
is calculated with this formula:
𝑊 𝑚 ℎ ℎ (1)
When the working fluid enters the turbine, it undergoes an isentropic expansion process. However, this
process does not occur without any losses, so the values are corrected using the turbine efficiency to
determine the true value at the turbine outlet. The following equation is used:
ℎ ℎ
𝜂 (2)
ℎ ℎ _
The mechanical power generated by the turbine is then converted into electricity by the generator. The
following equation is used to determine the generated power:
𝑊 𝑊 ∗𝜂 (3)
3.3.4 Pump
The pumps in the system are the main drivers to create flow. In this report, the modelled system has
two pumps, a feed pump and a wellhead pump. The power of these pumps is calculated as follows:
𝑣 , ∗ Δ𝑃 ∗ 100
𝑊 (4)
𝜂
Report 9 93 Andal
3.3.5 Condenser
There are various types of condensers that can be used in ORC power plants. Choosing the suitable
type of condenser is dependent on the weather and the availability of water.
The power requirements of these units were calculated using the following formula:
𝑣 , ∗ Δ𝑃 ∗ 100
𝑊 (5)
𝜂 ∗𝜂
When choosing a working fluid, several factors need to be considered. According to DiPippo (2016),
the performance of a fluid is dependent on the geothermal conditions, the type of binary cycle, and any
other operating or design constraint. Other factors are flammability, toxicity, chemical aggressiveness,
potential hazards to the environment, and cost. Based on the factors mentioned, this report considered
using either isopentane or N-pentane as a working fluid as the critical temperatures of these fluids fit
the resource temperature.
For this report, the author adopts the seven years’ cost estimation made by the DOE (2019c) for the pre-
development stage or exploration stage of a geothermal project in the Philippines as illustrated in Figure
8. This cost range, represented by the shaded area on Figure 8, was based on the records and data
obtained from the internal report from the DOE and varies depending on the size of the prospect area
and location.
Andal 94 Report 9
FIGURE 8: Geothermal exploration cost in the Philippines for the first seven years
Based on the report, the pre-development cost in the Philippines for seven years varies in the range USD
9 877 294 - 28 330 658 (PHP 516 780 000 - 1 482 260 000). The costs included in the report are as
follows:
Year 1 - Preparatory activities and acquisition of permits and clearances;
Year 2 - Preliminary geoscientific studies;
Year 3 - Detailed geoscientific surveys;
Year 4 - Drilling preparations;
Year 5 - Confirmatory drilling;
Year 6 - Delineation/development drilling; and
Year 7 - Project review and planning.
For this report, the author uses the minimum value listed in the DOE´s internal report as the modelled
prospect is a small binary geothermal power plant. Additionally, the computed cost estimate of a
conventional well is used rather than the one provided in the report as the computed value better reflects
on the cost of the project itself. Based on this, the total cost of investment for the pre-development stage,
without the cost of drilling of conventional well, is USD 2 232 000 (PHP 116 780 000).
3.4.2 Drilling
The author uses the total number of wells needed using the estimation made by the International Finance
Corporation (IFC). According to IFC (2013), the average well drilling success rate differs according to
the phase of the project: during the exploration phase an average of 59% of wells are successful; during
the development phase the average is 74%; and during the operation phase it increases to 83%. For this
study, the author determines the total number of wells based on a conservative well drilling success rate
of about 60%. It is also assumed that the unsuccessful or none producing wells could be used for
reinjection. In a real scenario this is not always possible as some wells may not be useful for either
production or reinjection. For the cost of drilling of each well, the author uses the cost estimation
provided in Hance (2005):
𝐶 240,785 210 𝐷 0.019069 𝐷 (6)
To get a better cost estimate of the major components of the power plant, the author explored several
methods that include using empirical formulas for cost estimation from the literature and consulting
Report 9 95 Andal
experienced professional with a wide background in such project. The cost estimation is summarized
in Tables 7 and 8:
Component, k Formula
𝑙𝑜𝑔 𝐶 k k 𝑙𝑜𝑔 𝑊 k 𝑙𝑜𝑔 𝑊
1
Turbine k1 = 2.6259; k2 = 1.4398; k3 = -0.1776
2
k1 = 2.2476; k2 = 1.4965; k3 = -0.1538
3
k1 = 2.7050; k2 = 1.4400; k3 = -0.1770
𝑙𝑜𝑔 𝐶 k k 𝑙𝑜𝑔 𝐴 k 𝑙𝑜𝑔 𝐴
1
k1 = 4.6656; k2 = -0.1557; k3 = 0.1547
Heat exchangers 2
k1 = 4.6656; k2 = -0.1557; k3 = 0.1547
3
k1 = 4.3247; k2 = -0.3030; k3 = 0.1634
𝑙𝑜𝑔 𝐶 k k 𝑙𝑜𝑔 𝑊 k 𝑙𝑜𝑔 𝑊
1
Pumps k 1 = 3.3892; k2 = 0.0536; k3 = 0.1547
2
k1 = 3.3892; k2 = 0.0536; k3 = 0.1547
3
k1 = 3.8700; k2 = 0.3160; k3 = 0.1220
1
Formula from heat exchanger was used
.
Condenser 2
𝐶 12300 ∗
3
Formula from heat exchanger was used
1
No formula available
.
Generator 2
𝐶 1,850,000
,
3
No formula available
1
El-Eman & Dincer (2013); 2 Lemmens (2016); 3 Jing, et al. (2017)
Component, k Formula
Turbine and generator 𝐶 𝑊 ∗ 400
Heat exchangers 𝐶 𝐴 ∗ 300
Pumps 𝐶 𝑊 ∗ 400
Air-cooled condenser 𝐶 𝐴 ∗ 600
For the cost estimation of turbine, generator and pump, the calculations could be easily done based on
the results of the power plant modelling. The area of the heat exchangers is calculated using the
following formula:
𝑄 𝑈 ∗ 𝐴 ∗ 𝐿𝑀𝑇𝐷 (7)
LMTD can be calculated following this formula:
𝑇 , 𝑇 , 𝑇 , 𝑇 ,
𝐿𝑀𝑇𝐷
𝑇 , 𝑇 , (8)
𝑙𝑛
𝑇 , 𝑇 ,
The subscripts “hot” and “cold” refer to the fluids listed in Table 9 and the subscripts “in” and “out”
refers to the inlet and outlet of the heat exchanger.
The heat transfer coefficient varies widely as this value is highly dependent on the materials and the
fluids in the heat exchanger. To simplify this part of the report, the following values were assumed (Dr.
Andal 96 Report 9
Páll Valdimarsson, adjunct professor at Reykjavík University, pers. comm., 2019) (Table 9):
𝒌𝑾
Process Heat transfer coefficient
𝒎𝟐 𝑲
Geothermal heat to isopentane (vapour) 0.8
Geothermal heat to isopentane (boiling) 1.2
Geothermal heat to isopentane (liquid) 0.5
Cooling air to isopentane (vapour) 0.5
Cooling air to isopentane (dew) 1.0
Isopentane (vapour) to isopentane (liquid) 0.3
The result of cost estimates from various source are presented in Appendix II. The results obtained from
different literature sources vary making it harder to get a good cost estimate for the component.
In consultation with the advisors, the author uses the cost estimates presented in Table 8, as this approach
is simple and updating the cost estimate according to the market price is easier compared to the formulas
presented in Table 7 of which some were already outdated.
Cost of O&M is 20% to 25% of the purchased equipment cost (PEC) (El-Eman and Dincer, 2013). For
this report, the author assumes that the annual cost of O&M is equal to the 20% of PEC.
The total cost of investment is the sum of the Fixed Capital Investment (FCI) and other outlays. The
FCI is the total of direct cost and indirect cost. The direct cost (DC) is composed but not limited to the
following: purchased equipment cost, purchased equipment installation, piping, instrumentations and
controls, electrical equipment and materials, land, civil, structural, and architectural work, and service
facilities. Indirect cost, on the other hand, is composed but not limited to the following: engineering and
supervision, construction costs including contractor´s profit, and contingencies. The other outlays costs
are composed of the following: start-up cost, working capital, cost of licensing, research and
development, and allowance for funds used during construction. The cost for the FCI used in this report
is summarized in Table 10.
In each project, it is important to evaluate the profitability of the investment before capital is invested.
However, calculating the profitability of an investment is complex as there are associated risks and
uncertainties. Analysing those risk and uncertainties is not covered in this report. Therefore, to
determine the profitability y of the model, the author adopts the deterministic investment analysis
presented in Bejan et al. (1996). In this analysis, the following was assumed:
1. There is a perfect capital market: the supply of funds is unrestricted;
2. There is complete certainty about investment outcomes;
3. Investment projects are indivisible; and
4. Profitability of one project does not in any way affect the profitability of any other project.
After setting up the conditions for the profitability evaluation, the Net Present Value (NPV) method and
the Internal Rate of Return (IRR) method are used for analysis.
Report 9 97 Andal
The selected prospect is currently one of the most advanced projects in the pre-development stage in the
Philippines. Despite of this, assumptions have to be made for the power plant model:
1. Well characteristics:
a. No. of required producing wells = 4;
b. Reservoir temperature = 190°C;
c. Reservoir pressure = 12.55 bar.
2. The temperature in the system plays a very important role in economics and the performance of
Andal 98 Report 9
the power plant. For this report, the following temperatures were assigned in the model:
a. Tcondensation = 47°C;
b. Tpinch = 4°C;
c. Tvap,sh = 2°C;
d. Tboilingmargin = 2°C;
e. Tcond,air,in = 27°C;
f. Tcond,air,out = 37°C;
g. Trecuperator,out = 52°C.
3. Efficiencies used in this report were assumed to be:
a. Turbine = 82%;
b. Generator = 95%;
c. Pump = 75%;
d. Condenser fan = 65%;
e. Fan motor = 98%.
4. As discussed in Section 3.3.6, the working fluid used in this report is isopentane at high pressure
(13 bars).
5. Due to high air humidity in the prospect area, using an ACC would be more cost effective
compared to evaporative cooling condenser.
6. Change in air pressure in the condenser is assumed to be 0.0017 bar.
7. As mentioned on Section 3.3, each well will be equipped with a well head pump that has a mass
flow rate of 40 kg/s and can deliver up to 3.5 bar.
8. While a high amount of NCG is present in fluids from SH-02, according to GeothermEx (2017),
it is believed that the fluid produced during the test is not representative of the deeper reservoir
fluid.
9. Thus, for the modelling the author assumed that the NCG content of the reservoir is negligible.
10. The geothermal fluid was assumed to be pure water.
11. To avoid silica saturation, the rejection temperature of the geothermal brine should be higher than
90°C.
12. Pressure loss due to friction in the system was disregarded.
13. For the profitability equation, the following values were assumed:
a. Effective discount rate = 10%;
b. Annual operations period = 340 days.
According to Clarke (2014), a binary geothermal power plant has four (4) major constraints that must
be met to ensure effective plant operations. The constraints are as follows:
1. The lower limit of the quality of vapour as it exits the turbine should be greater than 97% to prevent
excessive damage of the turbine blades.
2. The upper limit of evaporator pressure is two (2) MPa or 20 bars to limit the mechanical stresses
on the plant components and to ensure that the fluid remains sub-critical throughout the cycle.
3. The outlet temperature of the geothermal brine should exceed the working fluid temperature at the
inlet of the preheater by at least the pinch point temperature difference.
4. The problem of scaling should be considered.
4. THERMOECONOMIC MODELLING
In this section, the author describes the background of thermoeconomic modelling which is used to
determine the minimum cost of generation and to give an insight on which components could be
optimized to further improve the system. The parameters used in this model are defined in the
nomenclature in the back of this report.
Report 9 99 Andal
Thermoeconomics is the combination of energy analysis and economic principles and is used to get
information that is not available through conventional energy analysis and economic evaluations but
crucial to the design and operation of a cost-effective system (Bejan et al., 1996). It is also used to
balance expenditure or capital cost and exergy cost to estimate the minimum cost of the plant product
(Kotas, 1985). The cost balance expresses that the cost rate of the product, p, is equal to the sum of the
rate of expenditures or capital cost and the cost rate of fuel, f:
𝐶 , 𝐶 , 𝑍 𝑍 (9)
Each component has two cost rates, capital investment, CI and O&M. The cost rate of the capital
investment of a component is the product of the present worth factor for specific components of plant
equipment and the capital recovery factor (CRF) divided by the annual period of operation time 𝜏 of the
plant:
𝑃𝑊 ∗ 𝐶𝑅𝐹
𝑍 (10)
𝜏
The present worth factor is:
𝑆
𝑃𝑊 𝐶 (11)
1 𝑖
To obtain the O&M cost rate of a component, the component’s O&M cost should be expressed first to
its levelized value. This concept of levelization is defined as the use of time-value-of-money arithmetic
to convert a series of varying quantities to a financially equivalent constant quantity or annuity over a
specified time interval. The levelized value of O&M is expressed through the following equation:
𝑘∗ 1 𝑘
𝐶, & 𝐶 & ∗ 𝐶𝐸𝐿𝐹 𝐶 & ∗ ∗ 𝐶𝑅𝐹 (13)
1 𝑘
where k is calculated as:
1 𝑟
𝑘 (14)
1 𝑖
In this report, it is assumed that the nominal rate rn is 4%. The cost rate of O&M for a specific component
is:
𝐶, & ∗ 𝐶𝑅𝐹 𝐶
𝑍 ∗ (15)
𝜏 𝐶
Exergy analysis is an important tool for design and analysis of thermal systems and is used as the basis
of thermoeconomics. Exergy analysis is used to improve the effectiveness of energy resource use, as it
enables the determination of the location, cause, and true magnitude of waste and loss (Bejan et al.,
Andal 100 Report 9
1996). Exergy is defined as the maximum portion of energy that can be converted into work
(Valdimarsson, 2010) while the none converted part is called anergy. This means that the specific
exergy rate for a stream is limited to the environmental state of the process and is expressed as:
𝑒 ℎ ℎ 𝑇 𝑠 𝑠 (16)
and the exergy rate is expressed as:
𝐸 𝑚 ∗𝑒 (17)
In this report, the environmental state or the dead state, which is represented by subscript 0 in Equation
16, is identified as the ambient temperature.
In thermoeconomics, it is important to take note of the hidden cost. This cost is associated with exergy
destruction which can be revealed through thermoeconomic analysis. The exergy destruction is
expressed as:
𝐸 , 𝐸 , 𝐸, 𝐸 , (18)
The Exergy rate of component associated with fuel, 𝐸 , , refers to the exergy coming from the
geothermal system. The Exergy rate of component associated with product, 𝐸 , , refers to the exergy
produced by the component. The Exergy rate of loss of component, 𝐸 , , refers to the exergy stream
that flows from the outside of the component and was not used by other components in the system. The
Exergy rate of destruction of component, 𝐸 , , refers to the exergy destroyed in the stream.
Exergy costing is an approach in thermoeconomics in which it is believed that exergy is the only rational
basis for assigning costs to the interactions that a thermal system experiences with its surroundings and
with the sources of inefficiencies within it (Bejan et al., 1996). It is also an effective tool to evaluate
the cost effectiveness of thermals systems, used to evaluate and enhance the performance from both an
economic and exergetic point of view (Adefila et al., 2015).
Exergy costing is associated with the entering and exiting of streams and their associated rates of exergy
transfer. This involved cost balances of entering and exiting streams plus the appropriate charges due
to capital investments and operating and maintenance expenses on each component separately. This is
expressed in the following equation.
𝐶 , 𝐶 , 𝐶 , 𝐶, 𝑍 (19)
𝑐 , 𝐸 , 𝑐 , 𝐸 , 𝑐 , 𝐸 , 𝑐, 𝐸, 𝑍 (21)
When the component receives power, 𝐶 , would be moved to the right-hand side of the equation. 𝐶 ,
would be transferred to the left side if there is a heat transfer from the component.
When analysing a component, it is important to remember that the cost of exergy is based on the stream
that is entering and exiting a component. We can assume that the exergy cost per unit is known for all
Report 9 101 Andal
entering streams since they are the cost of the stream exiting the previous component and/or the cost of
investment of the component. The remaining unknown variable is the exergy cost that exits the
component. Additionally, when calculating exergy costing, it is important to consider the cost per
exergy unit. Lastly, since the exergy exiting a component contains the cost of stream exiting a previous
component, it is important to define a break point to break the loop. In this report, the author selected
the inlet of the pump (point 7 on Figure 5) as the break point of the system.
The first stream to analyse is in the feed pump. The first stream has one inlet stream, one outlet stream
and one power stream for work done by the pump on the system. The cost balance is calculated using
the following equation:
𝐶 , 𝐶 , 𝐶 , 𝑍 (22)
As mentioned earlier, the pump inlet is the break point. Therefore, the cost at pump inlet equals zero.
Additionally, to determine the cost rate of pump power, an auxiliary equation is used to calculate the
cost per exergy unit for the net power exported from the system while power input in the pump remains
constant.
𝐶 , 0 (23)
𝐶 , 𝐶 ,
(24)
𝑊 𝑊
The second stream of the working fluid is in the recuperator. In this stream, there are two (2) inlet
streams and two outlet streams. The cost balance in the recuperator is as follows:
𝐶 , , 𝐶 , ,
(25)
𝐶 , 𝐶 , 𝑍
In this stream, there are two outlet streams, a high-pressure side (hps) and a low-pressure side (lps),
therefore, an additional equation must be added. According to Bejan et al. (1996), the purpose of a heat
exchanger is to heat the cold stream and exergy is removed from the hot stream. This means that the
cost per exergy unit in the hot stream remains constant.
𝑐 , , 𝑐 , (26)
Before analysing the stream in the preheater and evaporator, the cost of stream of the source should be
calculated since itis added to the cost in these two components. When analysing the cost rate in the
production well, it should be noted that the cost rate in the reservoir is zero since the geothermal fluid
from the reservoir is natural to the environment and assumed to be free of cost. With this, the following
equation for the cost balance in the production well can be formulated:
𝐶 , , 𝑍 , 𝐶 , (27)
After passing the reservoir, the stream goes to the evaporator and the geothermal fluid enters a heat
exchanger. Similarly, to what was described for the second stream of the working fluid, the cost per
exergy unit in this hot stream remains constant:
𝑐 , , 𝑐 , , (28)
The next stream of the geothermal fluid is in the preheater. In the preheater, the cost rate of geothermal
fluid entering the preheater is already identified and the cost rate exiting the preheater is equal to the
cost rate of the reinjection well. The cost rate of the reinjection well is expressed as:
𝐶 , , 𝐶 , 𝑍 , 0 (29)
Andal 102 Report 9
According to the cost balance of the reinjection well, all cost rates are entering the component. This
means that the value obtained at the preheater outlet is negative. This suggests that the cost of the
reinjection well would also be carried out in the preheater in the system.
As the cost rates in the geothermal fluid side are now defined, the cost balance in the preheater in the
working fluid side can be calculated. In the preheater, there are two outlet streams and two (2) inlet
streams:
𝐶 , , 𝐶 , ,
(30)
𝐶 , , 𝐶 , , 𝑍
After the preheater, the next stream is in the evaporator. This evaporator has two outlet streams and two
(2) inlet streams:
𝐶 , , 𝐶 , ,
(31)
𝐶 , , 𝐶 , , 𝑍
After the evaporator, the next stream is in the turbine. In this stream, there is one outlet stream, one (1)
inlet stream and one (1) work done on the turbine by the system. The cost balance in the turbine is
expressed as:
𝐶 , 𝐶 , 𝐶 , , 𝑍 (32)
In this equation, there are two (2) exiting streams, one is the turbine outlet and the work done in the
turbine. In this scenario, an auxiliary relation is required. According to Bejan et al. (1996), since the
purpose of a turbine is to generate power, the exergy rate spent to generate the power and the exiting
exergy rate at the turbine would not change since cost would change only if exergy was added to the
working fluid during the turbine expansion. Therefore:
𝑐 , , 𝑐 , (33)
After the turbine, the stream enters the low-pressure side of the recuperator which was described earlier.
After the recuperator, the stream enters the last component of the system, the air-cooled condenser
(ACC), before the cycle is repeated. In the condenser, there are two outlet streams, two inlet streams
and one work done by the fan in the turbine. The cost balance in this component is expressed as:
𝐶 , , 𝐶 , , 𝐶 , , 𝐶 , , 𝐶 , 𝑍 (34)
In this cost balance, it could be concluded that the 𝐶 , , contains all the cost rate of the system,
except for the cost rate of the work done in the turbine and other cost rates outside the major components.
With this, the minimum cost of product or electricity can be expressed as:
𝐶 𝐶 , , 𝐶 , (35)
and in USD/kWh, that is:
𝐶 ∗ 3600
𝑀𝑖𝑛𝑖𝑚𝑢𝑚 𝑐𝑜𝑠𝑡 𝑜𝑓 𝑔𝑒𝑛𝑒𝑟𝑎𝑡𝑖𝑜𝑛 𝑈𝑆𝐷/𝑘𝑊ℎ (36)
𝑊
As mentioned in section 4.3, there are costs that are associated with the cost of exergy loss and exergy
destruction. The cost of exergy loss can be expressed as:
𝐶 , 𝐶 , 𝐶, 𝑍 (37)
Report 9 103 Andal
In Equation 39a, it is assumed that the product 𝐸 , is fixed and that the unit cost of fuel 𝑐 , is
independent to the exergy destruction. This implies that the cost of exergy destruction is equal to the
product of cost of fuel and exergy destruction:
𝐶 , 𝑐 , 𝐸 , (40a)
In Equation 39b, on the other hand, it is assumed that the fuel 𝐸 , is fixed and that the unit cost of the
product 𝑐 , is independent of the exergy destruction and expressed as:
𝐶 , 𝑐 , 𝐸 , (40b)
For this report, the author uses Equation 40b since the fuel rate from the geothermal system is assumed
to be fixed.
The exergy destruction and exergy loss provide thermodynamic measures of system inefficiencies.
These inefficiencies can be compared to the total exergy rate of the fuel of the system which is given as:
𝐸 ,
𝑦 (41)
𝐸,
𝐸,
𝑦 (42)
𝐸 ,
The relative cost difference 𝑟 expresses the relative increase of cost per exergy unit between the fuel
and product. This variable is used for iterative cost optimization where the objective is to minimize the
relative cost difference instead of minimizing the cost per exergy unit of the component to reveal the
real cost sources. This variable is can be expressed as:
1 𝜖 𝑍 , 𝑍 ,
𝑟 (43)
𝜖 𝑐 , 𝐸 ,
The exergetic efficiency is used to determine how well the exergy was utilized in the system. Exergetic
efficiencies of a binary power plants is usually within the range from 17.2 to 53.9% which is roughly
three to five times higher than the thermal efficiency which ranges from 2.1 to 10.3% (Haraldsson,
2016).
Andal 104 Report 9
The ratio of the non-exergy-related cost to the total cost is the exergo-economic factor (𝑓 ):
𝑍
𝑓 (45)
𝑍 𝑐 , 𝐸 , 𝐸,
A low exergo-economic factor implies that cost saving in the entire system could be achieved by
improving the components’ efficiency through the reduction of exergy destruction even if the capital
investment for this component increases. On the other hand, a high value of this factor suggests that the
capital investment could be lessened even if this would decrease the exergetic efficiency. For heat
exchangers the factor is typically lower than 55%, for compressor and turbines it lies between 35 to
75%, and for pumps it is typically above 70% (Bejan et al., 1996).
In geothermal power plant modelling the power plant design is mainly limited, among other factors, by
the geothermal resource and climate in the area. The parameters described in Section 3.10 and the
constraints set in Section 3.11 served as the foundation of the model. With these considerations, the
author developed a model using the Engineering Equation Solver (EES) software. The result of the
modelling shows the optimal power output of the system and is illustrated the power plant block diagram
in Figure 9.
The wells were considered to deliver 160 kg/s of 190°C hot water to the power plant. The heat from the
geothermal fluid is transferred to the pressurized isopentane through the evaporator and preheater before
it is reinjected. After the heat transfer, the isopentane boils and becomes slightly superheated reaching
temperatures of up to 131.5°C. This superheated isopentane is used to drive the turbine to generate 9390
kW. After passing the turbine, the fluid remains superheated, but its pressure drops to 1.877 bar. This
superheated vapour will then go to the recuperator to transfer heat to the working fluid that is located
between the feed pump and the preheater. The recuperator also helps to decrease the temperature in the
system prior to the ACC. Using ambient air, the ACC will now condense the isopentane at 47°C. Then,
this condensed fluid is pumped back at 13 bar and the cycle is repeated. The cost of exergy C will be
discussed in the thermoeconomic model.
Report 9 105 Andal
The modelled power plant has a parasitic load of 1327 kW in the ACC and of 413.1 kW in the feed
pump which reduces the power output to 7650 kW. With this, the thermal efficiency of the system is
calculated to be 5.92%. To further understand the thermal efficiency a Sankey diagram is used to
illustrate the heat flow in the system (Figure 10).
The Sankey diagram is a representation of the heat flow in the system where the length of the triangle
on the side of the box represents the amount of heat flow in each component. This allows a visualization
of the areas that mainly affects the thermal efficiency of the power plant. One of the significant
observations that can be made in the diagrams is that most of the unused heat is the reinjected heat in
the reinjection well and the ACC. However, optimizing the system to use this waste heat to generate
more electricity would not be easy because of the constraints on the geothermal fluid and the
environment that were mentioned. In the diagram, it can also be seen that the recuperator plays a
significant role through recovering waste heat and feeding it back into the system. Based on this model,
the cost of developing the prospect was also calculated. The calculated cost is summarized in Table 11.
The total cost of investment and the annual O&M cost of the prospect is 33.6 and 1.56 Million USD
(MUSD), respectively. The major cost of development is the drilling of wells followed by the total
PEC. Therefore, the cost of development of the prospect adds up to 3574 USD/kW. This cost analysis
will be further discussed in the thermoeconomic modelling.
After determining the generated net power and the cost of development, the NPV and IRR of the system
were calculated (Figure 11). Based on the profitability evaluation, the project will have a positive NPV
when the tariff is higher than 8.26 USDcents/kWh and a calculated IRR of 10% which is the break-even
point of the system.
Report 9 107 Andal
To create a thermoeconomic
model, the author uses the
information gathered in the
power plant model and the
parameters discussed in
section 4.8. Like the power
plant modelling, the author
uses the EES to generate the
thermoeconomic model. The
result of the modelling for the
cost of components and
exergy analysis is presented in
Table 12.
Zci ZO&M Zk
Ein Eout Ew Ed
Component k (103, (103, (103,
(kW) (kW) (kW) (kW)
USD/s) USD/s) USD/s)
Production well 56.086 0.3461 56.432 0.00 22 699.00 - -
Reinjection well 23.44 0 23.44 4 522.00 0.00 - -
Feed pump 0.5967 0.1786 0.7753 270.30 586.30 413.10 96.99
Recuperator 2.644 0.7916 3.4356
High-pressure side 586.30 1 415.00 - 277.60
Low-pressure side 4 980.00 3 874.00 - -
Preheater 2.626 0.7862 3.4122
GF side2 10 416.00 4 522.00 - -
WF side3 1 415.00 6 473.00 - 835.30
Evaporator 1.643 0.492 2.135
GF side2 22 699.00 10 416.00 - -
3
WF side 6 473.00 16 704.00 - 2 053.00
Turbine 13.56 4.061 17.621 16 704.00 4 980.00 -9 885.004 1 839.00
Condenser 7.148 2.14 9.288
WF side3 3 874.00 270.30 - 519.00
Air side 0.00 4 412.00 1 327.00 -
1 2 3
Cost of O&M at the well head pump; Geothermal fluid side; Working fluid side
4
The negative value indicates the value of exergy going out of the component.
*The parameters are defined in the nomenclature at the back of this report.
Andal 108 Report 9
The Grassman diagram shows the exergy flow in the system. The length of the triangle on the side of a
box represents the amount of exergy flow in each component. The diagram shows that the effectiveness
of the system in converting the exergy is high and reflects the exergetic efficiency of the power plant
which is 43.55%. It could also be observed that there was significant exergy destruction in the
evaporator and the turbine.
The analysis made with the Grassman diagram is important because, as mentioned in previous sections,
it is used to balance between expenditure or capital cost and exergy cost to estimate the minimum cost
of the plant product. With this definition, the author used the exergy cost, presented in Figure 9, to
determine the minimum cost of generation. To illustrate the relationship between the exergy cost and
the minimum cost of generation, the author uses the Value Flow diagram in Figure 13.
The Value Flow diagram shows the cost of exergy in the system. The length of the triangle on the side
of a box represents the cost of exergy in each component. This diagram identifies the major source of
cost in the system. The major contributor is the cost of production and reinjection, followed by the
cooling system and its parasitic load, and the turbine investment. As described earlier on this section,
the production and reinjection wells have the largest total investment and O&M cost rates in the system
followed by the turbine but, when the parasitic load was considered, the cost rate in the cooling system
becomes the second largest contributor.
The Value Flow diagram also shows that the minimum cost of generation is the sum of the cost of exergy
in the stream in the turbine and the stream in the condenser which is equivalent to 132.17x10-3 USD/sec
or 6.22 USDcents/kWh. This value represents the minimum cost of generating electricity of the system
based on the cost balance. This cost only covers the cost of the power generating system, therefore, the
real cost of generated electricity is expected to be higher as per the value determined using NPV analysis,
8.26 USDcents/kWh. Nevertheless, the minimum cost of generating electricity is much lower than the
true cost of diesel in Mindoro that ranges from 18 to 74 USDcents/kWh or the existing subsidized
approved generation rate of 11 USDcents/kWh.
The result for the cost of exergy is positive. Through evaluation of the Grassman diagram and the
calculated thermoeconomics variables, the system could be further optimized. To do this, the
components were arranged in descending order of the sum of cost of exergy destruction, investment and
O&M cost rate (Cd+Z) (Table 13).
From the thermoeconomics viewpoint, the turbine and the preheater have the highest value of Cd+Z and
are, therefore, the most important component of the system. The low value of exergoeconomic f of the
preheater and the evaporator shows that the cost associated with these components are comprised mostly
of the exergy destruction. Thus, cost saving might be achieved by improving the components efficiency
through the reduction of the exergy destruction even if the capital investment for this component will
increase. However, as for the turbine and the feed pump, an improvement of the evaporator would not
be possible as the quality of this component is based on the quality provided by the manufacturer.
Therefore, cost saving can obtain through improvements in the preheater.
6. CONCLUSION
Power plant modelling and thermoeconomics modelling was used to determine the optimal capacity of
the system, cost of development, minimum cost of generation, and exergetic evaluation.
The result of the modelling shows that the Montelago geothermal prospect could be capable of a
generating 9390 kW electricity with a net generating capacity of 7650 kW. The system has thermal
efficiency and exergetic efficiency of 5.92 and 43.55%, respectively. The obtained efficiencies are well
within the range of 2.1 to 10.3% for thermal efficiency and 17.2 to 53.9% for exergetic efficiency.
Andal 110 Report 9
Additionally, it was determined that the total cost of developing the prospect is 3574 USD/kW while
the minimum cost of generation is 6.22 USDcents/kWh. However, this cost only covers the cost of the
power generating system and does not cover the costs on offsite costs, indirect cost and other outlays,
therefore, the real cost of generated electricity is expected to be higher than this value. That is why it is
important to additionally analyse the system using the profitability evaluation. The profitability
evaluation shows that the project starts to have a positive NPV when the tariff is higher than 8.21
USDcents/kWh with IRR of 10% which is the break-even point of the system. This is way lower than
the current price of electricity available in the region which his 18 to 74 USDcents/kWh and still lower
than the existing subsidized approved generation rate which is 11 USDcents/kWh.
Additionally, the exergetic analysis found that cost saving in the entire system might be achieved by
improving the preheater´s exergetic efficiency through the reduction of the exergy destruction even if
the capital investment for this component will increase.
Therefore, it is concluded that, based on the described parameters, further developing the geothermal
prospect is economical. However, it needs to be taken into consideration that as the prospect is further
developed, new data might not be like the described parameters in this report. This means that the design
considerations and constraints used in this report cannot be directly applied. Nevertheless, the study
made on the Montelago geothermal prospect provided a wide overview using power plant and
thermoeconomics modelling. Additionally, the model could be used as a base line for other authors who
have an interest in conducting a similar study.
ACKNOWLEDGEMENTS
I want to express my sincerest gratitude to Lúdvík S. Georgsson, Ingimar Gudni Haraldsson and to the
rest of the staff of the UNU-GTP for giving me the opportunity to be part of this training and for
providing me with necessary support in the entire period of the program. Also, I would like to thank
my advisers, Dr. Páll Valdimarsson and Dr. María Sigrídur Gudjónsdóttir, for their guidance and
supervision during the entire project period, and Lýdur Skúlason for providing assistance and guidance
with my report.
Also, I would like to thank my colleagues at the Philippines Department of Energy and the Mindoro
Geothermal Power Corp. for providing me with all the necessary information and data that I needed to
complete this report.
Lastly, I would also like to thank my family, the UNU Fellows of 2019 and my friends for their support.
This report would not have been possible without these people.
NOMENCLATURE
A = Area (m2).
Ck = Cost of component (USD).
𝐶 = Cost of exergy (USD/s).
cp,k = Cost per exergy unit associated with fuel (USD/kJ).
ck = Cost per exergy unit (USD/kJ).
cp,k = Cost per exergy unit associated with product (USD/kJ).
𝐶 , = Cost of exergy destruction (USD/s).
𝐶, = Cost of exergy stream exiting (USD/s).
Report 9 111 Andal
REFERENCES
Adefila, S.S., Alam, M.M., Fagbenle R.O., and Oyedepo, S.O., 2015: Exergy costing analysis and
performance evaluation of selected gas turbine power plants. World J. Engineering, 12-2, 161-176.
Ahmed, S. J., 2019: The Philippine energy transition: building a robust power market to attract
investment, reduce prices, improve efficiency and reliability. Institute for Energy Economics and
Financial Analysis, report, 14 pp.
Asmin, S.N., Catigtig, D., Pratama, A.B., Regandara, R., Saputra, M.P., and Suryantini, 2016:
Geological, isothermal, and isobaric 3-D model construction in early stage of geothermal exploration.
Proceedings of the 5th ITB International Geothermal Workshop, Bandung, Indonesia,14 pp.
Axelsson, G., and Halldórsdóttir, S., 2015: Revised volumetric assessment of Montelago geothermal
system. ÍSOR – Iceland GeoSurvey, Reykjavík, report, ISOR-15020, 12 pp.
Bejan, A., Moran, M., and Tsatsaronis, G., 1996: Thermal design and optimization. John Wiley & Sons,
Inc., Canada, 533 pp.
Clarke, J., 2014: Optimal design of geothermal power plants. Virginia Commonwealth University,
Richmond, Virginia, 192 pp.
Clemente, J. V., Alcober, H.C., de Guzman, R.C., and Bayrante, L.F., 2016: Country update on
geothermal utilization and barriers affecting its growth – Philippines. Proceedings of the 11th Asian
Geothermal Symposium, Chiangmai, Thailand, 6 pp.
DOE, 2019b: GEMD-NREP 2018-2040 report. Department of Energy, Taguig, Philippines, unpublished
internal report.
DOE, 2019c: Cost of geothermal exploration in the Philippines. Department of Energy, Taguig,
Philippines, unpublished internal report.
DiPippo, R., 2016: Geothermal power generation: developments and innovation. Woodhead Publishing,
Kidlington, UK, 822 pp.
El-Eman, R.S., and Dincer, I., 2013: Exergy and exergoeconomic analyses and optimization of
geothermal organic Rankine cycle. Applied Thermal Engineering, 59-1/2, 435-444.
Fujii, H., Itoi, R., Jalilinasrabady, S., and Valdimarsson, P., 2011: Energy and exergy analysis of Sabalan
binary geothermal power plant. J. Geothermal Research Society of Japan, 33-3, 113-121.
GeothermEx, 2017: Resource evaluation for the Montelago geothermal project, Oriental Mindoro, The
Philippines. GeothermEx, Richmond, CA, report, 89 pp.
Hance, C.N., 2005: Factors affecting costs of geothermal power development. Geothermal Energy
Association for the US Department of Energy, 64 pp.
Report 9 113 Andal
Haraldsson, I.G., 2016: Efficiency in geothermal utilization processes. Presented at “SDG Short Course
I on Sustainability and Environmental Management of Geothermal Resource Utilization and the Role
of Geothermal in Combating Climate Chage”, organized by UNU-GTP and LaGeo, in Santa Tecla, El
Salvador, 66 pp.
IFC, 2013: Success of geothermal wells: a global study. International Finance Corporation, World Bank
Group, Washington, DC, 76 pp.
Jing, D., Liu, P., Tian, H., Shu, G., and Wang, X., 2017: Operational profile based thermal-economic
analysis on an organic Rankine cycle using for harvesting marine engine´s exhaust waste heat. Energy
Conversion and Management, 2017, 107-123 pp.
Kotas, T.J., 1985; The exergy method of thermal plant analysis. Anchor Brendon Ltd, Tiptree, Essex,
Great Britain, 296 pp.
Lemmens, S., 2016: Cost of engineering techniques and their applicability for cost estimation of organic
Rankine cycle systems. Energies, 9-7, 18 pp.
Lukawski, M., 2009: Design and optimization of standardized organic Rankine cycle power plant for
European conditions. The School for Renewable Energy Science, Akureyri, Iceland, 76 pp.
NGCP, 2019: Transmission development plan 2016-2040 final report volume 1: major network
development. NGCP.ph, website: www.ngcp.ph/Attachment-Uploads/TDP%202016-2040%20Final
%20Report%20Volume%201%20Major%20Network%20Development-2019-05-14-16-43-41.pdf
SKM, 2011: Assessment of the Montelago prospect, Philippines. Green Earth Energy Ltd, report, 38 pp.
Thórhallsson, S., 2005: Common problems faced in geothermal generation and how to deal with them.
Presented at “Workshop for Decision Makers on Geothermal Projects and Management”, organized
by UNU-GTP and KenGen, Naivasha, Kenya, 12 pp.
Valdimarsson, P., 2010: Production of electricity from a geothermal source. In: Geothermal energy,
Macedonian Geothermal Association, Skopje, Macedonia, 50-180.
van Leeuwen, W.A., 2016: The Montelago geothermal prospect. In: Geothermal exploration using the
magnetotelluric method. Utrecht University, Dept. of Earth Sciences, the Netherlands, 133-194.
APPENDIX II: Major component cost estimates in the power plant, in USD