Residual Saturation under Pressure-Controlled Drainage

तकनीकी सारांश: दबाव-नियंत्रित जल निकासी के तहत अवशिष्ट संतृप्ति (Residual Saturation)

समस्या विवरण
यह शोध पत्र सरंध्र माध्यमों (porous media) के भौतिकी में एक मौलिक अनसुलझे प्रश्न को संबोधित करता है: अनंत-आकार की सीमा (infinite-size limit) में दबाव-नियंत्रित जल निकासी (pressure-controlled drainage) के तहत अवशिष्ट संतृप्ति का व्यवहार। जबकि मानक आक्रमण पारगमन (invasion percolation - IP) मॉडल भविष्यवाणी करते हैं कि जैसे-जैसे सिस्टम का आकार अनंत की ओर बढ़ता है, आक्रमणकारी चरण की संतृप्ति लुप्त हो जाती है (ब्रेकथ्रू क्लस्टर की फ्रैक्टल प्रकृति के कारण ही रक्षात्मक तरल का एक सीमित अंश फंसा रह जाता है), यह स्पष्ट नहीं है कि क्या एक दबाव-नियंत्रित प्रक्रिया—जो आक्रमणकारी चरण को दिए गए कैपिलरी थ्रेशोल्ड पर इनलेट से जुड़े सभी छिद्रों तक पहुँचने की अनुमति देती है—अनंत सीमा में रक्षात्मक चरण की एक परिमित अवशिष्ट संतृप्ति (finite residual saturation) का परिणाम देती है। यह अंतर क्षेत्र-स्तर की घटनाओं जैसे वॉटरफ्लडिंग (waterflooding) को समझने के लिए महत्वपूर्ण है, जहाँ यह प्रश्न उठता है कि क्या जलाशय के आयामों के बढ़ने पर अवशिष्ट तेल संतृप्ति परिमित रहती है।

कार्यप्रणाली
लेखक दबाव-नियंत्रित जल निकासी को एक पोर-नेटवर्क ग्राफ पर ट्रैपिंग के साथ बॉन्ड परकोलेशन (Bond Percolation with Trapping - BPT) के रूप में प्रतिपादित करते हैं।

• ग्राफ प्रतिनिधित्व: सरंध्र माध्यम को एक अनडिरेक्टेड ग्राफ के रूप में मॉडल किया गया है जहाँ नोड्स पोर क्षेत्रों का प्रतिनिधित्व करते हैं और किनारे (edges) प्रवेश त्रिज्याओं (कैपिलरी थ्रेशोल्ड) के साथ थ्रोट्स (throats) का प्रतिनिधित्व करते हैं।
• प्रक्रिया परिभाषा: मानक IP के विपरीत, जो स्थानीय रूप से सबसे सुलभ पथ के माध्यम से थ्रोट-दर-थ्रोट आगे बढ़ता है, BPT एक वैश्विक दबाव प्रतिबंध के तहत कार्य करता है। जैसे-जैसे कैपिलरी दबाव बढ़ता है, "सक्रिय ग्राफ" (जिसमें वर्तमान दबाव से नीचे के प्रवेश थ्रेशोल्ड वाले थ्रोट्स शामिल हैं) का विस्तार होता है।
• ट्रैपिंग तंत्र: आक्रमण को इनलेट पर सक्रिय ग्राफ के साथ जुड़ाव (connectivity) के रूप में परिभाषित किया गया है। ट्रैपिंग तब होती है जब रक्षात्मक चरण का एक हिस्सा आउटलेट सेट से अपना जुड़ाव खो देता है। एल्गोरिदम सिस्टम को अपडेट करता है: (1) नए थ्रोट्स को सक्रिय करने के लिए दबाव बढ़ाना; (2) आक्रमण क्षेत्र की पहचान करना (इनलेट-कनेक्टेड सक्रिय पोर); (3) रक्षात्मक क्षेत्र की पहचान करना (आउटलेट-कनेक्टेड गैर-आक्रमण वाले पोर); और (4) शेष डिस्कनेक्टेड पोर को "ट्रैप्ड" (trapped) के रूप में वर्गीकृत करना।
• सिमुलेशन: अध्ययन विभिन्न 2D लैटिस (स्क्वायर, ट्राइएंगुलर, वोरोनोई) और 3D लैटिस (डायमंड, विभिन्न समन्वय संख्याओं के साथ सिंपल-क्यूबिक) पर BPT की तुलना मानक आक्रमण पारगमन (IP) से करता है। लॉग-नॉर्मल वितरित थ्रोट रेडियस और परिमित-आकार प्रभावों (finite-size effects) का विश्लेषण करने के लिए एन्सेम्बल एवरेजिंग का उपयोग किया गया है।

प्रमुख योगदान और परिणाम

1. अनंत सीमा में परिमित अवशिष्ट संतृप्ति:
   प्राथमिक निष्कर्ष यह है कि दबाव-नियंत्रित जल निकासी एक स्पेस-फिलिंग आक्रमण अवस्था की ओर ले जाती है जहाँ रक्षात्मक चरण की अवशिष्ट संतृप्ति, जैसे-जैसे सिस्टम का आकार $L \to \infty$ होता है, एक गैर-शून्य, परिमित सीमा ($0 < S_\infty < 1$) की ओर बढ़ती है। यह मानक IP के बिल्कुल विपरीत है, जहाँ ब्रेकथ्रू के समय आक्रमणकारी संतृप्ति लुप्त हो जाती है ($S_{IP} \to 0$) क्योंकि आक्रमणकारी क्लस्टर फ्रैक्टल होता है और रक्षात्मक तरल के बड़े क्षेत्रों को बायपास कर देता है।

2. दो आयामी (2D) में परिमित-आकार स्केलिंग:
   2D में, ऊष्मागतिकी सीमा से अवशिष्ट संतृप्ति का विचलन एक परिमित-आकार स्केलिंग नियम का पालन करता है:
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   जहाँ घातांक (exponent) $\delta \approx 0.25$ है। यह घातांक विभिन्न लैटिस प्रकारों (स्क्वायर, ट्राइएंगुलर, वोरोनोई) में सार्वभौमिक है और सूक्ष्म विवरणों से स्वतंत्र है। धीमा अभिसरण ($\Delta S_r \sim L^{-0.25}$) यह दर्शाता है कि 2D में प्रतिनिधि एलीमेंटरी वॉल्यूम (representative elementary volumes) प्राप्त करने के लिए अत्यंत बड़े सिस्टम आकार की आवश्यकता होती है (उदाहरण के लिए, सिस्टम के आकार को दोगुना करने से परिमित-आकार पूर्वाग्रह केवल $\approx 1/16$ के कारक से कम होता है)।

3. तीन आयामी (3D) में आयामीता और समन्वय संख्या प्रभाव:
   3D में, परिमित-आकार सुधार काफी तेजी से घटते हैं, जिसमें घातांक $\delta \approx 0.75$ है। फलस्वरूप, 3D में अनंत-आकार सीमा का अनुमान लगाने के लिए आवश्यक प्रतिनिधि आयतन 2D की तुलना में बहुत छोटा है। इसके अलावा, जबकि 3D में स्पर्शोन्मुखी (asymptotic) अवशिष्ट संतृप्ति परिमित रहती है, इसका मान लैटिस की समन्वय संख्या (coordination number) पर निर्भर करता है; उच्च कनेक्टिविटी अधिक वैकल्पिक आक्रमण पथ प्रदान करती है, जिससे फंसा हुआ अंश कम हो जाता है, हालांकि यह इसे पूरी तरह से समाप्त नहीं करता है।

4. इनवेशन परकोलेशन (Invasion Percolation) से अंतर:
   यह शोध पत्र स्पष्ट करता है कि BPT और IP के बीच का अंतर केवल वक्र के आकार या बर्स्ट डायनेमिक्स का मामला नहीं है, बल्कि मौलिक रूप से अलग स्पर्शोन्मुखी अवस्थाओं का मामला है। IP ब्रेकथ्रू पर एक फ्रैक्टल आक्रमण क्लस्टर देता है, जबकि BPT एक स्पेस-फिलिंग आक्रमण क्षेत्र देता है जिसमें अनंत सीमा में एक परिमित फंसा हुआ अंश सह-अस्तित्व में रहता है।

महत्व
यह शोध पत्र यह प्रदर्शित करके परकोलेशन सिद्धांत और दबाव-संतृप्ति संबंधों के बीच एक सीधा सैद्धांतिक संबंध स्थापित करता है कि जुड़ाव (connectivity) और ट्रैपिंग दबाव-नियंत्रित जल निकासी के तहत अवशिष्ट संतृप्ति के निर्धारण कारक हैं। यह दिखाकर कि अनंत-सिस्टम सीमा में भी छिद्र स्थान का एक परिमित अंश फंसा रह जाता है, यह कार्य मानक आक्रमण-परकोलेशन चित्र को ब्रेकथ्रू अवस्था से आगे बढ़ाता है। यह मैक्रोस्कोपिक कंस्ट्रिट्यूटिव लॉ (constitutive laws) की व्याख्या करने के लिए एक कठोर ढांचा प्रदान करता है और यह सुझाव देता है कि "प्रतिनिधि एलीमेंटरी वॉल्यूम" की समस्या आयामीता के प्रति अत्यधिक संवेदनशील है, जो 2D में 3D की तुलना में बहुत अधिक प्रतिबंधात्मक है। परिणाम दर्शाते हैं कि दबाव-नियंत्रित परिदृश्यों में, रक्षात्मक तरल का एक परिमित अंश सिस्टम के आकार के बावजूद अप्राप्य रहता है, जो उपसतह (subsurface) द्रव पुनर्वितरण और रिकवरी प्रक्रियाओं के लिए प्रत्यक्ष निहितार्थ रखता है।