Heat Transfer Modeling in Enhanced Geothermal Energy: A Three-Temperature Approach for Solid, Injected, and Residing Fluids

यहाँ शोध पत्र "Heat Transfer Modeling in Enhanced Geothermal Energy: A Three-Temperature Approach for Solid, Injected, and Residing Fluids" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

एन्हांस्ड जियोथर्मल सिस्टम्स (EGS) गहरे, गर्म, फ्रैक्चरयुक्त सरंध्र चट्टानों (porous rock) में ठंडे वर्किंग फ्लूइड्स को इंजेक्ट करने पर निर्भर करते हैं ताकि ऊष्मा (heat) निकाली जा सके। जलाशय के थर्मल विकास (thermal evolution) का सटीक पूर्वानुमान लगाना बिजली उत्पादन और जलाशय के जीवनकाल का अनुमान लगाने के लिए महत्वपूर्ण है।

मौजूदा मॉडलों की सीमाएँ:

• लोकल थर्मल इक्विलिब्रियम (LTE): यह मानता है कि चट्टान और तरल पदार्थ एक ही तापमान साझा करते हैं। यह तब विफल हो जाता है जब प्रवाह दर (flow rates) के सापेक्ष ऊष्मा विनिमय (heat exchange) धीमा होता है।
• मानक लोकल थर्मल नॉन-इक्विलिब्रियम (LTNE): ठोस और तरल तापमान के बीच अंतर करता है, लेकिन तरल को एक एकल, औसत घटक के रूप में मानता है।
• मुख्य मुद्दा: EGS में, ठंडा इंजेक्ट किया गया तरल और गर्म निवासी (resident) तरल एक ही पोर स्पेस (pore space) में सह-अस्तित्व में होते हैं। मानक मॉडल जलाशय के माध्यम से यात्रा करते समय इंजेक्ट किए गए तरल के थर्मल इतिहास को स्पष्ट रूप से हल नहीं कर पाते हैं। वे बल्क औसत (bulk averages) पर निर्भर करते हैं, जिससे थर्मल चैनलिंग (thermal channeling), फिंगरिंग अस्थिरता (fingering instabilities) और उत्पादन कुओं तक पहुँचने से पहले इंजेक्ट किए गए तरल के विशिष्ट हीटिंग पथों जैसी महत्वपूर्ण घटनाओं को समझना कठिन हो जाता है।

2. कार्यप्रणाली (Methodology)

लेखक एक नवीन थ्री-वे लोकल थर्मल नॉन-इक्विलिब्रियम (Three-Way LTNE) फ्रेमवर्क प्रस्तावित करते हैं जो उन्नत संख्यात्मक विविक्तकरण (numerical discretization) तकनीकों के साथ जुड़ा हुआ है।

A. गणितीय मॉडलिंग (Mathematical Modeling)

नियंत्रित प्रणाली डार्सी फ्लो (Darcy flow) को तीन अलग-अलग ऊर्जा समीकरणों और एक सांद्रता (concentration) समीकरण के साथ जोड़ती है:

1. प्रवाह (Flow): दबाव $p$ और वेग $u$ द्वारा नियंत्रित सिंगल-फेज डार्सी फ्लो।
2. सांद्रता ($z$): एक स्केलर वेरिएबल ($0 \le z \le 1$) जो इंजेक्ट किए गए तरल के आयतन अंश (volume fraction) को ट्रैक करता है।
3. तीन-तापमान क्षेत्र (Three-Temperature Fields):
   • $T_i$: इंजेक्ट किए गए तरल का तापमान।
   • $T_r$: निवासी तरल का तापमान।
   • $T_s$: ठोस मैट्रिक्स (solid matrix) का तापमान।
4. कपलिंग तंत्र (Coupling Mechanism):
   • कुल तरल एन्थैल्पी को विभाजित किया जाता है: $C_f T_f = C_i T_i + C_r T_r$, जहाँ $C_i$ और $C_r$, सांद्रता $z$ पर निर्भर करते हैं।
   • ऊष्मा विनिमय तीनों चरणों के बीच होता है:
     • इंजेक्ट किया गया तरल $\leftrightarrow$ ठोस ($\Lambda_{is}$)
     • निवासी तरल $\leftrightarrow$ ठोस ($\Lambda_{rs}$)
     • इंजेक्ट किया गया तरल $\leftrightarrow$ निवासी तरल ($\Lambda_{ir}$)
   • कपलिंग स्ट्रेंथ को सांद्रता $z$ और सरंध्रता (porosity) $\phi$ द्वारा स्केल किया जाता है।

B. संख्यात्मक विविक्तकरण (Numerical Discretization)

यह प्रणाली दृढ़ता से जुड़ी हुई (strongly coupled) और एडवेक्शन-प्रधान (advection-dominated) है, जिसके लिए मजबूत संख्यात्मक योजनाओं की आवश्यकता होती है ताकि गैर-भौतिक दोलन (non-physical oscillations) को रोका जा सके और संरक्षण (conservation) सुनिश्चित किया जा सके।

1. स्थानिक विविक्तकरण (Spatial Discretization): एनरिच्ड गैलरकिन (Enriched Galerkin - EG) विधि

   • यह निरंतर $Q_1$ (लीनियर) तत्वों को विच्छिन्न (discontinuous) $Q_0$ (पिसवाइस कांस्टेंट) एनरिचमेंट के साथ जोड़ता है।
   • लाभ: यह सुनिश्चित करता है कि स्थानीय द्रव्यमान संरक्षण (जो परिवहन समस्याओं के लिए महत्वपूर्ण है) बना रहे, जबकि पूर्ण डिस्कंटीन्यूअस गैलरकिन (DG) विधियों की तुलना में अपेक्षाकृत कम डिग्री ऑफ फ्रीडम (degrees of freedom) बनाए रखता है।
   • इसे दबाव, तापमान और सांद्रता समीकरणों पर लागू किया गया है।

2. स्थिरीकरण (Stabilization): फ्लक्स-करेक्टेड ट्रांसपोर्ट (FCT)

   • इसे विशेष रूप से सांद्रता और तापमान एडवेक्शन पदों पर लागू किया गया है।
   • तंत्र: यह एक प्रेडिक्टर-करैक्टर दृष्टिकोण का उपयोग करता है। एक लो-ऑर्डर मोनोटोन स्कीम एक स्थिर आधार प्रदान करती है, और हाई-ऑर्डर सुधार जोड़े जाते हैं लेकिन डिस्क्रीट मैक्सिमम प्रिंसिपल को लागू करने और स्प्यूरियस ऑसिलेशन (spurious oscillations) को दबाने के लिए एक FCT लिमिटर द्वारा सीमित किए जाते हैं।

3. समय एकीकरण (Time Integration): IMPES-प्रकार का स्प्लिटिंग

   • इम्पलेस प्रेशर, एक्सप्लिसिट कंसंट्रेशन, इम्पलिस हीट (IMPRESS):
     • दबाव/वेग: स्पष्ट रूप से (implicitly) हल किया जाता है।
     • सांद्रता ($z$): एक स्ट्रॉन्ग-स्टेबिलिटी-प्रिजर्विंग रनगे-कुटा (SSP-RK2) स्कीम का उपयोग करके स्पष्ट रूप से (explicitly) हल किया जाता है।
     • तापमान ($T_i, T_r, T_s$): स्टिफ डिफ्यूसिव टर्म्स को संभालने के लिए BDF2 (बैकवर्ड डिफरेंशिएशन फॉर्मूला ऑफ ऑर्डर 2) स्कीम का उपयोग करके स्पष्ट रूप से (implicitly) हल किया जाता है।

3. मुख्य योगदान (Key Contributions)

1. नवीन थ्री-वे LTNE फॉर्मूलेशन: यह सांद्रता वेरिएबल का उपयोग करके एक निरंतर फ्रेमवर्क के भीतर इंजेक्ट किए गए और निवासी तरल तापमान को स्पष्ट रूप से अलग करने वाला पहला मॉडल है। यह इंजेक्ट किए गए तरल के हीटिंग पथों को सीधे ट्रैक करने की अनुमति देता है।
2. ऊष्मागतिक स्थिरता (Thermodynamic Consistency): मॉडल को कठोरता से व्युत्पन्न किया गया है ताकि यह सुनिश्चित हो सके कि अनंत ऊष्मा विनिमय की सीमा में, यह मानक LTE में बदल जाता है, और अनंत इंटरफेशियल क्षेत्र की सीमा में, यह मानक टू-टेम्परेचर LTNE में बदल जाता है।
3. मजबूत संख्यात्मक फ्रेमवर्क: EG, FCT, और SSP-RK2/BDF2 स्प्लिटिंग का एकीकरण एडवेक्शन-प्रधान जियोथर्मल समस्याओं के लिए एक स्थिर, स्थानीय रूप से संरक्षित और दोलन-मुक्त सॉल्वर प्रदान करता है।
4. छिपे हुए थर्मल स्ट्रक्चर का समाधान: यह फ्रेमवर्क उन थर्मल डायनेमिक्स (जैसे सांद्रता फ्रंट्स पर तीव्र तापमान जंप) को प्रकट करता है जिन्हें पारंपरिक सिंगल-टेम्परेचर मॉडल में औसत (average) कर दिया जाता है।

4. संख्यात्मक परिणाम (Numerical Results)

लेखकों ने तीन संख्यात्मक प्रयोगों के माध्यम से अपने मॉडल को सत्यापित किया:

• अभिसरण परीक्षण (Convergence Test): चिकनी समाधानों (smooth solutions) के लिए युग्मित EG-FCT स्कीम के लिए इष्टतम अभिसरण दर ($O(h^2)$ $L^2$ नॉर्म में) प्रदर्शित की गई, जो कार्यान्वयन की सटीकता की पुष्टि करती है।
• LTE बनाम LTNE तुलना:
  • तेज़ ऊष्मा विनिमय: दिखाया गया कि जब इंटरफेशियल हीट ट्रांसफर अधिक होता है, तो थ्री-वे मॉडल LTE समाधान की ओर अभिसरित होता है, जो इसके एसिम्प्टोटिक विश्लेषण को मान्य करता है।
  • धीमी ऊष्मा विनिमय: यह प्रदर्शित किया गया कि धीमी विनिमय (जो EGS में आम है) के तहत, LTE चरणों के बीच तापमान अंतर को पकड़ने में विफल रहता है। LTNE मॉडल तरल और ठोस तापमान के बीच महत्वपूर्ण विचलन की सही भविष्यवाणी करता है, जो एडवेक्शन द्वारा बढ़ाया जाता है।
• विषम फ्रैक्चरयुक्त माध्यम (Heterogeneous Fractured Medium - EGS परिदृश्य):
  • कम पारगम्यता वाली चट्टान में अंतर्निहित 34 उच्च-पारगम्यता वाले फ्रैक्चर चैनलों वाले डोमेन के माध्यम से प्रवाह का अनुकरण किया गया।
  • मुख्य निष्कर्ष: मॉडल ने थर्मल ब्रेकथ्रू (thermal breakthrough) व्यवहार को सफलतापूर्वक हल किया। इसने दिखाया कि फ्रैक्चर चैनलों के साथ इंजेक्ट किया गया तरल तेजी से गर्म होता है जबकि निवासी तरल ठंडा होता है।
  • एन्थैल्पी ट्रैकिंग: मॉडल ने इंजेक्ट किए गए तरल ($T_i$) और मिश्रित तरल तापमान के बीच अंतर किया। इसने खुलासा किया कि जबकि कुछ क्षेत्रों में $T_i$ उच्च हो सकता है, इंजेक्ट किए गए तरल द्वारा ले जाई जाने वाली एन्थैल्पी सांद्रता द्वारा भारित (weighted) होती है, जो ऊर्जा निष्कर्षण की अधिक सटीक तस्वीर प्रदान करती है।
  • तुलना: मानक औसत मॉडल तीव्र तापमान ग्रेडिएंट और इंजेक्ट किए गए तरल के विशिष्ट हीटिंग प्रक्षेपवक्र (trajectories) को पकड़ने में विफल रहे।

5. महत्व और प्रभाव (Significance and Impact)

• बेहतर EGS डिज़ाइन: इंजेक्ट किए गए तरल के थर्मल विकास को स्पष्ट रूप से हल करके, इंजीनियर थर्मल ब्रेकथ्रू समय (जब ठंडा तरल उत्पादन कुएं तक पहुँचता है) का बेहतर पूर्वानुमान लगा सकते हैं और ऊष्मा निष्कर्षण को अधिकतम करने के लिए इंजेक्शन रणनीतियों को अनुकूलित कर सकते हैं।
• भौतिक शुद्धता (Physical Fidelity): यह मॉडल फ्रैक्चरयुक्त मीडिया में औसत तरल तापमान की "ब्लैक बॉक्स" प्रकृति को संबोधित करता है, जो चैनलिंग और स्थानीयकृत हॉट स्पॉट जैसी मल्टीस्केल घटनाओं का अध्ययन करने के लिए एक उपकरण प्रदान करता है।
• कंप्यूटेशनल दक्षता: Enriched Galerkin विधि का उपयोग DG विधियों के सख्त संरक्षण और निरंतर गैलरकिन विधियों की दक्षता के बीच एक संतुलन प्रदान करता है, जिससे बड़े पैमाने पर 3D सिमुलेशन संभव हो जाता है।
• भविष्य का कार्य: लेखक नोट करते हैं कि जबकि वर्तमान प्रसार गुणांक (diffusion coefficients) अत्यधिक एडवेक्शन प्रभुत्व से बचते हैं, भविष्य के कार्य के लिए अत्यधिक संवहन-प्रधान शासन (highly convection-dominated regimes) के लिए अतिरिक्त स्थिरीकरण और पॉइंट सोर्स (वेल) के समावेश की आवश्यकता हो सकती है।

संक्षेप में, यह शोध पत्र एक भौतिक रूप से कठोर तीन-तापमान मॉडल पेश करके जियोथर्मल जलाशय सिमुलेशन में एक महत्वपूर्ण प्रगति प्रस्तुत करता है, जो इंजेक्ट किए गए तरल, निवासी तरल और रॉक मैट्रिक्स के बीच जटिल थर्मल इंटरैक्शन को कैप्चर करता है, जिसे एक स्थिर और संरक्षित संख्यात्मक सॉल्वर द्वारा समर्थित किया गया है।