مدل‌سازی عددی مخزن زمین‌گرمایی در محیط غیراشباع (مطالعۀ موردی: میدان زمین‌گرمایی شمال‌غرب سبلان)

در این پژوهش، ابتدا یک مدل مفهومی طراحی‌شده از 10 چاه عمیق اکتشافی توسعه داده شده ‌و بر این اساس، مدل عددی سه‌بعدی از سیستم زمین‏گرمایی شمال‏ غرب سبلان ارائه شده است. به همین منظور، برای شبیه‏سازی مدل از مدول eos3 (معادلۀ حالت آب- هوا) کد شبیه‏ساز tough2 استفاده شده است. مدل عددی مخزن با یک منشور مستطیلی با km 11.5 طول و km 8 عرض (km2 92) و km 5.11 عمق بیان شده و با 21 لایۀ افقی با محدودۀ ضخامت 100 تا 1000 متر از ماکزیمم ارتفاع 4110 تا 1000- متر از سطح دریا گسترش داده شده است. تعداد 22 تیپ‏سنگ که توزیع آن‏ها بر اساس توزیع واحدهای سنگی و ساختارهای زمین‏شناسی در 10 چاه اکتشافی عمیق بنا نهاده شده است، به مدل اختصاص داده شد. نفوذپذیری این ساختارهای زمین از 17-10×1 تا 13-10×9 مترمربع متغیر است. مدل‏سازی برای حالت طبیعی مخزن انجام شد و نتایج اعتبارسنجی، تطبیق خوبی بین داده‏های اندازه‏گیری‌شده درون‌چاهی و مدل‏سازی انجام‌شده برای دما و فشار را نشان می‌دهد. در بهترین حالت تطبیق، شبیه‏سازی یک زون جریان دما بالای رو به بالا در بخش جنوب‏ شرقی محدوده (زیر سایت‏های d و e)، فراهم کرد. این جریان از طریق زون‏های نفوذپذیر، گسلی و شکسته انتقال می‌یابد و درنهایت، از طریق سیماهای سطحی (چشمه‏های آب‏گرم) در بخش شمال ‏غربی ناحیه تخلیه می‏شود. نتایج مطالعه مدل مفهومی مناسبی از مخزن زمین‏گرمایی سبلان را ارائه می‏دهد که می‏تواند در توسعه‏های آتی مخزن استفاده شود.

numerical modeling of geothermal reservoir in unsaturated environment (case study: geothermal field northwest of sabalan)

in this research, firstly, a conceptual model designed of ten deep exploratoration wells has been developed, and based on this, a three-dimensional numerical model of the north-west sablan geothermal system has been presented. for this purpose, the eos3 module (water-air equation) of the tough2 simulator code was used to simulate the model. the numerical model of the reservoir is expressed by a rectangular prism with a length of 11.5 km, a width of 8 km (92 km2) and a depth of 5.11 km, and it was expanded with 21 horizontal layers with a thickness range of 100 to 1000 meters from the maximum altitude of 4110 to -1000 masl. the number of 22 rock-types, whose distribution is based on the distribution of rock units and geological structures in 10 deep exploratory wells, was assigned to the model. the permeability of the rock types varies from 1*10e-17 to 9*10e-13 square meters. modeling was done for the natural state of the reservoir and the validation results showed a good match between the data measured in the well and modeled for temperature and pressure. in the best fit, the simulation provided a zone of high temperature upward flow in the southeastern part of the range (sites d and e). this flow is transferred through permeable, faulted and fractured zones and is finally discharged through surface features (hot water springs) in the northwestern part of the area.