Microsoft Word - 211.docx CHEMICAL ENGINEERING TRANSACTIONS VOL. 61, 2017 A publication of The Italian Association of Chemical Engineering Online at www.aidic.it/cet Guest Editors: Petar S Varbanov, Rongxin Su, Hon Loong Lam, Xia Liu, Jiří J Klemeš Copyright © 2017, AIDIC Servizi S.r.l. ISBN 978-88-95608-51-8; ISSN 2283-9216 Geothermal Power Potential Assessment in North Eastern Morocco Alae-Eddine Barkaouia,*, Yassine Zarhloulea, Petar Sabev Varbanovb, Jiří Jaromír Klemešb a Laboratory of Hydrogeology-Environment, Faculty of Sciences, University Mohamed Ist, Oujda, Morocco. b Sustainable Process Integration Laboratory (SPIL), NETME Centre, Faculty of Mechanical Engineering, Brno University of Technology, Brno, Czech Republic a.barkaoui@ump.ac.ma The volumetric method is one of the most applicable methods of low-temperature geothermal resource assessment. While applying volumetric method, the values of uncertain parameters should be determined. An add-in software program to Microsoft EXCEL, @RISK, is used as a tool to define the uncertainties of the parameters in the volumetric equation. Monte Carlo simulation is used as the probabilistic approach for the assessment of low temperature in north-eastern Morocco. In this region, the utilisation of this resource is limited. Assessment studies using triangular and uniform distribution type functions for each parameter are used to give the mean values of recoverable heat energy of the field. 1. Introduction Comparing to other energy sources, geothermal energy is a clean, renewable, constant and available worldwide. Due to its reduction capacity in greenhouse gas emissions, this energy is already being used for electricity generation and direct utilization (Barkaoui, 2016). For Morocco, the use of geothermal energy, among other renewable resources, is primordial in order to break the current dependence on fossil fuel and electricity imports (Barkaoui, 2013). The Kingdom of Morocco is the only North African country with no natural oil resources and is the largest energy importer in the region with 96 % of its energy needs being sourced externally. Morocco has small quantities of gas and it has large reserves of oil shale. However, in the absence of a proven specific industrial process that can produce oil and gas from this unconventional source, Morocco has turned to implementing a number of strategies that promote renewable energy and energy efficiency. In 2009, the total installed capacity and the electricity generation in Morocco reached the levels of 6,370 MW and 21 TWh, respectively. 4,6 TWh was imported from Spain to recover the power demand which reached 25 TWh. In 2008 Morocco launched the national energy strategy, with renewable energy and energy efficiency plan as the main pillars. The country has one of the most ambitious renewable energy programs in the region. It expects 42 % (equivalent to about 6,000 MW) of its total energy mix to come from renewable sources by 2020. In Morocco, thermal waters are mainly hosted within sedimentary reservoirs, consisting of Liassic limestones with a thickness up to 500 m. The geothermal fluid is characterized by a complex deep circulation and it ascends through complex fault systems. The Liasic reservoir of North-eastern province is considered as the most important geothermal aquifer in the country. This reservoir feeds more than twenty-four thermal manifestations, with temperatures ranging from 26 to 54 °C. Some of these hot aquifers, e.g. Fezouane, near Berkane and Hammam Ben Kachour, at Oujda play an important role in the economy of the area – as shown in the investigation on geothermal potentials in Morocco overall (Zarhloule, 1999) and a study on the surface geothermal potentials (Zarhloule et al., 2001). Geodynamic studies linked the zones showing geothermal gradient and heat flow exceeding 50 °C/km and 100 mW/m 2 respectively, to Neogene - quaternary volcanic and neotectonic activities. However, these thermal phenomena are still not developed and their exploitation limited to drinkable water distribution or to balneotherapy “ancient Hamam”. DOI: 10.3303/CET1761269 Please cite this article as: Barkaoui A.-E., Zarhloule Y., Varbanov P.S., Klemeš J.J., 2017, Geothermal power potential assessment in north eastern morocco, Chemical Engineering Transactions, 61, 1627-1632 DOI:10.3303/CET1761269 1627 The aim of this work is to approach the power potential in the study area. After successful exploration of a geothermal prospect, stakeholders are always eager to have those results. This comes as early as after completion of surface geo-scientific exploration or even after initial exploration drilling. The earlier estimates of power potential give confidence to the project owners to source for more resources to undertake subsequent stages of development. With high uncertainty and scanty data available during initial stages of exploration, stochastic and risk analysis methods are frequently used to estimate the range and probable distribution of stored heat reserves and hence, exploitable energy base of the newly explored geothermal prospect or fields. These methods have been borrowed from the oil industry where they have been used for a long time to estimate probabilistic hydrocarbon-in-place and oil and gas reserves in sedimentary basins. 2. Geothermal potential Geological and hydrogeological data from boreholes show the Liassic carbonates to be the main hydrogeothermal reservoir in the region. This reservoir is highly variable in thickness. The meteoric waters penetrate from the surface through the outcrops of the Liasic limestones in the southern part of the Angad plain, continues flowing downward through the same formation that becomes deeper going to the north. According to Zarhloule (1999), the hot temperature and the artesian rise of most of the thermal springs are due to groundwater circulating at depth within a framework of a recent volcanic area and a system of basement faults, forming horsts and grabens. Winckel (2002) performed a thorough geochemical analysis of the main thermal water sources in Morocco and found that eleven of them release CO2 and are partially of deep origin. These water sources are mainly located on a NE-SW line from Nador to Taza, and from Fes (Moulay Yacoub) to Oulmes south of the Rif frontal thrust, along the so-called Moroccan Hot Line (MHL). Tassi et al. (2006) confirmed that CO2-rich thermal waters with 3He anomalies are likely related to MHL. The contemporary presence of 3He anomalies and minor recent basalt outcrops indicate that CO2 originates from mantle degassing or deep hydrothermal systems in these thermal discharges. The regional pattern highlights heat flux increasing north-eastward, from less than 60 (north Mauritania) to more than 80-90 mW·m-2 in the eastern Rif, north-eastern Morocco, Alboran Sea and north-western Algeria. The largest values of the geothermal gradient are fond in the north-eastern part of Morocco, where they can reach 50 °C·km-1 To understand better the behaviour of the thermal water inside the liasic geothermal reservoir, many water boreholes were logged, especially in the north-eastern part of the country. Among the recorded thermal profiles, Figure 1 shows one interesting example for well 1624/7, located west of Berkane. This hole is characterized by an increase in geothermal gradient at 300 m depth from 29 to 127 °C km-1. At the same depth, the lithology changes from clay to dolomite. At about 470 m depth, the temperature is about 50 ºC. The shape of the thermal profile suggests a conductive thermal regime both in the upper (clay) and in the lower (carbonate) section of the hole. The dolomitic formation continues until the hole bottom (1,042 m depth). By extrapolating the thermal gradient inferred in the lowermost section of the hole, a bottom temperature of about 120 °C is inferred. The lithology change cannot explain the increase of geothermal gradient. As dolomite is expected to have much greater thermal conductivity than clay (see e.g. thermal conductivity data for NW Morocco rocks (Zarhloule et al., 2007) and the recent compilation by Pasquale et al., 2011), one would expect the geothermal gradient to decrease. An explanation for the anomalous pattern of the thermal gradient might be found in the advective heat transfer, which can occur at depth in the carbonate formation. We may argue that heat advection occurring in the main deep thermal aquifer, encountered at 1,042 m depth, can yield the increase of thermal gradient observed in the overlying dolomitic layers. 3. Thermal energy calculation An important factor in a geothermal assessment is an assessment of the volume of the geothermal system in question using the volumetric method. We assume, for simplicity, that the volume is a box with a surface area A in the xy plane and height (thickness) z1-z0 along the z-axis, where z1 and z0 are the lower and upper limit of the geothermal system, respectively. When the volume of the geothermal system has been assessed, the choice has to be made on how to calculate the useable heat that the system contains. For simplicity, it can be assumed that the heat capacity and temperature are homogeneous in the xy plane and are only dependent on depth. The volumetric method is used to estimate the amount of energy in a geothermal resource. This method involves calculating the amount of thermal energy contained in a given volume of rock and water and then estimating how much of this energy maybe recoverable given a reference temperature. The volumetric method uses the volume of the rock, the specific heat and temperature of rock to calculate the energy (Pálmason, 2005). This method is patterned from the work applied by the USGS to the Assessment of 1628 Geothermal Resources of the United States (Muffler and Cataldi, 1978). The equation used in calculating the thermal energy for a liquid dominated reservoir is as follows: QT = Qr + Qw (1) Qr = A.h * (ρr.Cr * (1-Ø) * (Ti-Tf) (2) Qw = A.h * (ρwCw * Ø * (Ti-Tf)) (3) Where: QT = total thermal energy, kJ/kg Qr = heat in rock, kJ/kg Qw = heat in water, kJ/kg A = area of the reservoir, m2 h = average thickness of the reservoir, m Cr = specific heat of rock at reservoir conditions, kJ/kgK Cw = specific heat of water at reservoir conditions, kJ/kgK Ø = porosity Ti = average temperature of the reservoir, °C Tf = final or abandonment temperature, °C ρr = rock density, kg/m3 ρw = water initial density, kg/m3 Figure 1: Thermal profile of the borehole 1624-7 located in the region of Berkane 3.1 Power plant sizing The above calculations only provide for the total thermal energy in place in the reservoir. To size the power that could be supported by the resource, the following equation is further introduced. 1629 = × × × (4) where P = Power potential (MWe); Rf = Recovery factor; Ce = Conversion efficiency; Pf = Plant factor; and t = Time in years (economic life): 3.2 The Monte Carlo Simulation The reserves estimation is done using commercial software that provides for a probabilistic approach of calculating uncertainty in the occurrence of events or unknown variables. The most common commercial software is @Risk which is used in assessing risks in investment, pharmaceuticals, petroleum reserves and mining evaluation. Monte Carlo simulation can also be programmed using an Excel spreadsheet. In this study, to obtain a good representation of the distribution, sampling is done through 1,000 iterations with continuous calculation. 3.3 Parameters used in the model Based on the geological and geophysical data, the geothermal reservoir underlies the entire town and its suburbs (200 km2), but its permeability in the major part is not high enough for reasonable geothermal development. The maximum temperature input into the calculations was the highest recorded temperature 52 °C in surface. Higher temperature can be reached if drilling deeper. By extrapolating those values and using the geothermal gradient of 49 °C/km a temperature of 120 °C can be expected at a depth of 1,042 m. The fluid density and specific heat capacity into the simulator were obtained from steam tables based on the reservoir temperatures. The values for those two parameters are respectively 800 Kg/m3 and 4248 J/Kg°C. The deep geothermal reservoir in North-Eastern Morocco is considered to be composed mostly of sedimentary rocks ; therefore, the value for the heat capacity of the rock, is equal to 920 J/Kg°C,. The porosity φ of the Limestone rocks in the study area is of the order of 10 %. The average density of the rock is set at 2,900 kg/m3 on the basis of the gravity available data. The geothermal recovery factor (Rf) refers to the fraction of the stored heat in the reservoir that could be extracted to the surface. It is dependent on the fraction of the reservoir that is considered permeable and on the efficiency by which heat could be swept from these permeable channels. This factor is assumed to be around 0.25 as a most likely value. The Conversion efficiency takes into account the conversion of the recoverable thermal energy into Electricity, in this study the value of this parameter is around 6 %. The plant factor refers to the plant availability throughout the year taking into consideration the period when the plant is scheduled for maintenance, or whether the plant is operated as a base-load or peaking plant. The good performance of many geothermal plants around the world places the availability factor to be from 90-97 %. In this study, a 95 % load factor is used. The economic life of the project is the period it takes the whole investment to be recovered within its target internal rate of return. This is usually 25-30 y. 4. Results and discussion The estimated generation capacities are only preliminary estimates since they are based on parameters with considerable uncertainties. At the present time, geothermal water is used directly in North Eastern Morocco. A small geothermal power plant of 2.09 MWe might be installed there in the next few years. An estimate of the electric power, which could be produced from the recoverable heat with the cut-off temperature of 120 °C from the geothermal reservoir, has been calculated. This was done for 25 y production time scenario. The results are presented as a relative frequency plot and discreet cumulative probability distribution. Each simulation consists of up to 1,000 iterations. From these random outcomes, miscellaneous statistical information can be found. These include the likeliest outcome, 90 % confidence interval, mean and median of the outcomes, standard deviation and where the 90 % limit for the cumulative probability lies. According to the statistics of the probability distribution, it is seen that the volumetric model predicts that with 90 % confidence the power production is around 0.861 MWe for 25 y. The minimum energy that can be produced is around 0.36 MWe, while the maximum that can be produced is around 4.93 MWe (Figure 2). It should be emphasised that the great range of values resulting from the Monte Carlo calculations simply reflects the uncertainty in the results obtained by the volumetric assessment method. It is primarily caused by uncertainty in the size, temperature and recovery factor for the geothermal reservoir resource. 1630 The probability that the output is greater than or equal to 0.861 MW is 90 % and the probability that the output is greater than or equal to 3.672 MW is 5 % (Figure 3). These results imply that the field could initially support a 0.861 MW power plant for 25 y; possible expansion will be subject to further delineation drilling and availability of field performance data. The risk that the field could not sustain 3.672 MW is equal to or less than 5 %. Figure 2: Relative frequency plot of the volumetric reserves estimation in NorthEastern Morocco Figure 3: Cumulative frequency plot of the volumetric reserves estimation in NorthEastern Morocco 5. Conclusion A Monte Carlo volumetric capacity assessment, based on the available data, has been performed. An estimate for the electric power, which can be produced from the recoverable heat in Northeastern Morocco has been calculated. According to the results, 2.09 MWe can be produced. For more certitude and with 90 % confidence, 0.861 MWe can be produced in North-Eastern Morocco if the recoverable heat is used for 25 y. The wide range of these estimates simply reflects the uncertainty in the size, temperature and recovery factor of the study area. 1631 Acknowledgments The authors would like to acknowledge the financial support by the EC and Croatian Ministry of Science Education and Sports projects “INTERGEORES” (NEWFELPRO Grants Agreements 43), as well as the support of the project "Sustainable Process Integration Laboratory – SPIL", project No. CZ.02.1.01/0.0/0.0/ 15_003/0000456 funded by EU "CZ Operational Programme Research and Development, Education", Priority 1: Strengthening capacity for quality research. References Barkaoui A.E., Boldyryev S., Duic N., Krajacic N., Guzović Z., 2016, Appropriate integration of geothermal energy sources by Pinch approach: Case study of Croatia. Applied Energy, In Press. doi:10.1016/j.apenergy.2016.04.112. Barkaoui A.E., Nemet A., Varbanov P.S., Klemeš J.J, Zarhloule Y., Rimi A., 2013, Integration of Geothermal Energy in the Case of North Eastern Morocco. Chemical Engineering Transactions, 32, 247-252, DOI: 10.3303/CET1332042. Muffler P., Cataldi R., 1978, Methods for regional assessment of geothermal resources, Geothermics 7, 53. Pálmason G., 2005, Geothermal book -Nature and use of resource (in Icelandic: Jarðhitabók - Eðli og nýting auðlndar). Reykjavík: Hið íslenska bókmenntafélag, Iceland. Pasquale V., Verdoya M., Chiozzi P., 2011, Groundwater flow analysis using different geothermal constraints: The case study of Acqui Terme area, northwestern Italy, Journal of Volcanology and Geothermal Research, 99, 38–46. Tassi F., Vaselli O., Moratti G., Piccardi L., Minissale A., Poreda R., Delgado Huertas A., Bendkik A., Chenakeb M., Tedesco D., 2006, Fluid geochemistry versus tectonic setting: the case study of Morocco. Geological Society, Special Publications, 262, 131-145, London, United Kingdom. Zarhloule, Y., 1999, Geothermal potentials of Morocco: integrated treatment by deep temperatures and surface indices. (in French), Doctorat d‘Etat, Fac Sci Oujda, Maroc. Zarhloule, Y., Lahrach, A., Ben Aabidate, L., Bouri, S., Boukdir A., Khattach, D. and Ben Dhia, H., 2001, Surface geothermal prospecting in Morocco: hydrodynamics, geothermal anomalies and surface indices. (in French), Journal African of Earth Sciences, 32, 851-867. Zarhloule Y., Verdoya M., El Mandour A., Chiozzi P., Boughriba M., Lahrach A., 2007, Hydrogeothermal Characters of the Moroccan Atlas, Proceedings, IUGG XXIV General Assembly “Earth, Our Changing Planet”, 2-13 July 2007, Perugia, Italy. Winckel A., 2002, Establishment of a typology of thermal waters by a hydrochemical, isotopic and tectonic approach. Example of Morocco. (in French), Thèse de Doctorat Université Paris Sud, Paris, France. 1632