A three-dimensional numerical model to simulate Iranian NW Sabalan geothermal system

Mirmahdi Seyedrahimi-Niaraq, Faramarz Doulati Ardejani, Younes Noorollahi, Soheil Porkhial, Ryuichi Itoi, Saeid Jalilinasrabady

研究成果: ジャーナルへの寄稿記事

2 引用 (Scopus)

抄録

The plan for electricity generation from north-west (NW) Sabalan geothermal field, Iran has been presented since 1994. Now construction of the first pilot plant is running with a capacity of 5 MWe in this site. A three-dimensional (3D) numerical model of fluid flow and heat transfer of the geothermal reservoir was developed using Tough2 simulator code on the basis of the updated conceptual model from the field recent data. The aim of this study was to present an updated 3D model of the geothermal reservoir in this area by taking ten deep wells data into account and validate an existing model in which the data obtained from only three wells were included. The model covers an area of 92 km2 and extends to an approximate depth of 5.11 km from the land surface. The model includes 17 horizontal layers, (AA through QQ). The thickness of the layers range between 100 m and 1000 m. The model extends from a minimum level of -1000 to a maximum level of 4110 masl. The model has 12,784 grid blocks, in total. Each layer has 752 grid blocks with horizontal dimensions of 250 m × 250 m and 500 m × 500 m. Altogether, 21 rock types were used in the model to assign different horizontal permeability values from 1.0 × 10−17 to 9.0 × 10-13 m2 based on the conceptual model. The model was first calibrated by varying the permeability, the total magnitude and location of the recharges and well bottom pressure of discharge zones. It was further validated by the profiles of temperature and pressure distribution at 10 deep exploration wells data. A close agreement was obtained between the model predictions and the measured downhole temperature and pressure data in the exploration wells. In the best-fitted model, the natural state simulation proved the presence of a high temperature upflow zone in the southern part (Below the ground, between Pads D and E). High temperature fluid at a rate of 108 kg/s and with an enthalpy of 1056.43 kJ/kg (244 °C) recharged over an area of 2.25 km2 from 36 grid blocks from the bottom layer in southern part of the field. The overall natural discharge is about 43 kg/s from hot springs which was simulated using the pressure dependent deliverability method. The results show that there is mostly one inflow of geothermal water at the level of -1000 m at southeast below the land surface between the Pad D and Pad E. The fluid flow moves upward by the fracture zones and faults toward the northwest of area. The faults NW5, NW3 and NNW2 play a major role in this system. Finally, the geothermal fluid discharging by the surface manifestations can be seen as hot springs at the land surface. This model can be used to simulate future production scenarios to evaluate the sustainability of the reservoir.

元の言語英語
ページ(範囲)42-61
ページ数20
ジャーナルGeothermics
77
DOI
出版物ステータス出版済み - 1 1 2019

Fingerprint

geothermal system
Numerical models
Hot springs
well
land surface
thermal spring
fluid flow
Flow of fluids
permeability
Geothermal fields
bottom pressure
Fluids
fluid
electricity generation
Pilot plants
fracture zone
enthalpy
Pressure distribution
Temperature
heat transfer

All Science Journal Classification (ASJC) codes

  • Renewable Energy, Sustainability and the Environment
  • Geotechnical Engineering and Engineering Geology
  • Geology

これを引用

A three-dimensional numerical model to simulate Iranian NW Sabalan geothermal system. / Seyedrahimi-Niaraq, Mirmahdi; Doulati Ardejani, Faramarz; Noorollahi, Younes; Porkhial, Soheil; Itoi, Ryuichi; Jalilinasrabady, Saeid.

:: Geothermics, 巻 77, 01.01.2019, p. 42-61.

研究成果: ジャーナルへの寄稿記事

Seyedrahimi-Niaraq, M, Doulati Ardejani, F, Noorollahi, Y, Porkhial, S, Itoi, R & Jalilinasrabady, S 2019, 'A three-dimensional numerical model to simulate Iranian NW Sabalan geothermal system', Geothermics, 巻. 77, pp. 42-61. https://doi.org/10.1016/j.geothermics.2018.08.009
Seyedrahimi-Niaraq, Mirmahdi ; Doulati Ardejani, Faramarz ; Noorollahi, Younes ; Porkhial, Soheil ; Itoi, Ryuichi ; Jalilinasrabady, Saeid. / A three-dimensional numerical model to simulate Iranian NW Sabalan geothermal system. :: Geothermics. 2019 ; 巻 77. pp. 42-61.
@article{085ba615dbd14a75bac443a625e63ab6,
title = "A three-dimensional numerical model to simulate Iranian NW Sabalan geothermal system",
abstract = "The plan for electricity generation from north-west (NW) Sabalan geothermal field, Iran has been presented since 1994. Now construction of the first pilot plant is running with a capacity of 5 MWe in this site. A three-dimensional (3D) numerical model of fluid flow and heat transfer of the geothermal reservoir was developed using Tough2 simulator code on the basis of the updated conceptual model from the field recent data. The aim of this study was to present an updated 3D model of the geothermal reservoir in this area by taking ten deep wells data into account and validate an existing model in which the data obtained from only three wells were included. The model covers an area of 92 km2 and extends to an approximate depth of 5.11 km from the land surface. The model includes 17 horizontal layers, (AA through QQ). The thickness of the layers range between 100 m and 1000 m. The model extends from a minimum level of -1000 to a maximum level of 4110 masl. The model has 12,784 grid blocks, in total. Each layer has 752 grid blocks with horizontal dimensions of 250 m × 250 m and 500 m × 500 m. Altogether, 21 rock types were used in the model to assign different horizontal permeability values from 1.0 × 10−17 to 9.0 × 10-13 m2 based on the conceptual model. The model was first calibrated by varying the permeability, the total magnitude and location of the recharges and well bottom pressure of discharge zones. It was further validated by the profiles of temperature and pressure distribution at 10 deep exploration wells data. A close agreement was obtained between the model predictions and the measured downhole temperature and pressure data in the exploration wells. In the best-fitted model, the natural state simulation proved the presence of a high temperature upflow zone in the southern part (Below the ground, between Pads D and E). High temperature fluid at a rate of 108 kg/s and with an enthalpy of 1056.43 kJ/kg (244 °C) recharged over an area of 2.25 km2 from 36 grid blocks from the bottom layer in southern part of the field. The overall natural discharge is about 43 kg/s from hot springs which was simulated using the pressure dependent deliverability method. The results show that there is mostly one inflow of geothermal water at the level of -1000 m at southeast below the land surface between the Pad D and Pad E. The fluid flow moves upward by the fracture zones and faults toward the northwest of area. The faults NW5, NW3 and NNW2 play a major role in this system. Finally, the geothermal fluid discharging by the surface manifestations can be seen as hot springs at the land surface. This model can be used to simulate future production scenarios to evaluate the sustainability of the reservoir.",
author = "Mirmahdi Seyedrahimi-Niaraq and {Doulati Ardejani}, Faramarz and Younes Noorollahi and Soheil Porkhial and Ryuichi Itoi and Saeid Jalilinasrabady",
year = "2019",
month = "1",
day = "1",
doi = "10.1016/j.geothermics.2018.08.009",
language = "English",
volume = "77",
pages = "42--61",
journal = "Geothermics",
issn = "0375-6505",
publisher = "Elsevier Limited",

}

TY - JOUR

T1 - A three-dimensional numerical model to simulate Iranian NW Sabalan geothermal system

AU - Seyedrahimi-Niaraq, Mirmahdi

AU - Doulati Ardejani, Faramarz

AU - Noorollahi, Younes

AU - Porkhial, Soheil

AU - Itoi, Ryuichi

AU - Jalilinasrabady, Saeid

PY - 2019/1/1

Y1 - 2019/1/1

N2 - The plan for electricity generation from north-west (NW) Sabalan geothermal field, Iran has been presented since 1994. Now construction of the first pilot plant is running with a capacity of 5 MWe in this site. A three-dimensional (3D) numerical model of fluid flow and heat transfer of the geothermal reservoir was developed using Tough2 simulator code on the basis of the updated conceptual model from the field recent data. The aim of this study was to present an updated 3D model of the geothermal reservoir in this area by taking ten deep wells data into account and validate an existing model in which the data obtained from only three wells were included. The model covers an area of 92 km2 and extends to an approximate depth of 5.11 km from the land surface. The model includes 17 horizontal layers, (AA through QQ). The thickness of the layers range between 100 m and 1000 m. The model extends from a minimum level of -1000 to a maximum level of 4110 masl. The model has 12,784 grid blocks, in total. Each layer has 752 grid blocks with horizontal dimensions of 250 m × 250 m and 500 m × 500 m. Altogether, 21 rock types were used in the model to assign different horizontal permeability values from 1.0 × 10−17 to 9.0 × 10-13 m2 based on the conceptual model. The model was first calibrated by varying the permeability, the total magnitude and location of the recharges and well bottom pressure of discharge zones. It was further validated by the profiles of temperature and pressure distribution at 10 deep exploration wells data. A close agreement was obtained between the model predictions and the measured downhole temperature and pressure data in the exploration wells. In the best-fitted model, the natural state simulation proved the presence of a high temperature upflow zone in the southern part (Below the ground, between Pads D and E). High temperature fluid at a rate of 108 kg/s and with an enthalpy of 1056.43 kJ/kg (244 °C) recharged over an area of 2.25 km2 from 36 grid blocks from the bottom layer in southern part of the field. The overall natural discharge is about 43 kg/s from hot springs which was simulated using the pressure dependent deliverability method. The results show that there is mostly one inflow of geothermal water at the level of -1000 m at southeast below the land surface between the Pad D and Pad E. The fluid flow moves upward by the fracture zones and faults toward the northwest of area. The faults NW5, NW3 and NNW2 play a major role in this system. Finally, the geothermal fluid discharging by the surface manifestations can be seen as hot springs at the land surface. This model can be used to simulate future production scenarios to evaluate the sustainability of the reservoir.

AB - The plan for electricity generation from north-west (NW) Sabalan geothermal field, Iran has been presented since 1994. Now construction of the first pilot plant is running with a capacity of 5 MWe in this site. A three-dimensional (3D) numerical model of fluid flow and heat transfer of the geothermal reservoir was developed using Tough2 simulator code on the basis of the updated conceptual model from the field recent data. The aim of this study was to present an updated 3D model of the geothermal reservoir in this area by taking ten deep wells data into account and validate an existing model in which the data obtained from only three wells were included. The model covers an area of 92 km2 and extends to an approximate depth of 5.11 km from the land surface. The model includes 17 horizontal layers, (AA through QQ). The thickness of the layers range between 100 m and 1000 m. The model extends from a minimum level of -1000 to a maximum level of 4110 masl. The model has 12,784 grid blocks, in total. Each layer has 752 grid blocks with horizontal dimensions of 250 m × 250 m and 500 m × 500 m. Altogether, 21 rock types were used in the model to assign different horizontal permeability values from 1.0 × 10−17 to 9.0 × 10-13 m2 based on the conceptual model. The model was first calibrated by varying the permeability, the total magnitude and location of the recharges and well bottom pressure of discharge zones. It was further validated by the profiles of temperature and pressure distribution at 10 deep exploration wells data. A close agreement was obtained between the model predictions and the measured downhole temperature and pressure data in the exploration wells. In the best-fitted model, the natural state simulation proved the presence of a high temperature upflow zone in the southern part (Below the ground, between Pads D and E). High temperature fluid at a rate of 108 kg/s and with an enthalpy of 1056.43 kJ/kg (244 °C) recharged over an area of 2.25 km2 from 36 grid blocks from the bottom layer in southern part of the field. The overall natural discharge is about 43 kg/s from hot springs which was simulated using the pressure dependent deliverability method. The results show that there is mostly one inflow of geothermal water at the level of -1000 m at southeast below the land surface between the Pad D and Pad E. The fluid flow moves upward by the fracture zones and faults toward the northwest of area. The faults NW5, NW3 and NNW2 play a major role in this system. Finally, the geothermal fluid discharging by the surface manifestations can be seen as hot springs at the land surface. This model can be used to simulate future production scenarios to evaluate the sustainability of the reservoir.

UR - http://www.scopus.com/inward/record.url?scp=85052902225&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85052902225&partnerID=8YFLogxK

U2 - 10.1016/j.geothermics.2018.08.009

DO - 10.1016/j.geothermics.2018.08.009

M3 - Article

AN - SCOPUS:85052902225

VL - 77

SP - 42

EP - 61

JO - Geothermics

JF - Geothermics

SN - 0375-6505

ER -