A New Insight into Mathematical modeling for sand production prediction due to high-speed fluid flow in loosely Consolidated Sandstones

तकनीकी सारांश: रेत उत्पादन के पूर्वानुमान के लिए गणितीय मॉडलिंग में एक नई अंतर्दृष्टि

समस्या विवरण
खराब रूप से संकुचित तेल और गैस भंडारों (अनकन्फाइंड कंप्रेसिव स्ट्रेंथ < 1000 psi) में रेत का उत्पादन (sand production) पेट्रोलियम उद्योग के लिए एक गंभीर चुनौती पेश करता है। यह वेलबोर अस्थिरता, उपकरणों का क्षरण, उत्पादन में व्यवधान, और महत्वपूर्ण पर्यावरणीय एवं आर्थिक लागतों का कारण बनता है। जबकि सैंड कंट्रोल विधियाँ मौजूद हैं (यांत्रिक, रासायनिक और प्रतिबंधित उत्पादन दर), प्रभावी नियंत्रण सटीक पूर्वानुमान पर अत्यधिक निर्भर करता है। मौजूदा मॉडलिंग दृष्टिकोण आम तौर पर स्ट्रेन-आधारित, इरोजन-आधारित या पार्टिकल-आधारित श्रेणियों में आते हैं। हालाँकि, एक ऐसे व्यापक मॉडल की आवश्यकता है जो उच्च-वेग वाले तरल प्रवाह के तहत कणों के अलग होने (detachment), परिवहन (transport) और पुन: आसंजन (re-adhesion) के भौतिकी को एकीकृत कर सके ताकि आयतनी रेत उत्पादन और परिणामस्वरूप भंडार गुणों में होने वाले परिवर्तनों का पूर्वानुमान लगाया जा सके।

कार्यप्रणाली
यह शोध पत्र एक संयुक्त इरोजन-आधारित और पार्टिकल-आधारित संख्यात्मक मॉडल प्रस्तावित करता है ताकि एक-आयामी (रैखिक) और अक्षीय सममित (रेडियल) प्रवाह ज्यामिति में रेत उत्पादन का पूर्वानुमान लगाया जा सके। कार्यप्रणाली इस प्रकार संरचित है:

1. कण पृथक्करण का भौतिकी (Physics of Particle Detachment): अध्ययन एक पोर वॉल (pore wall) से जुड़े कण पर कार्य करने वाले बलों के मोमेंट बैलेंस (moment balance) के आधार पर एक नवीन सैंडिंग मानदंड विकसित करता है। इसमें ड्रैग, लिफ्ट, उछाल (buoyancy) और DLVO (डेरजागिन-लैंडौ-वेरवे-ओवरबीक) बलों पर विचार किया गया है। विश्लेषण निष्कर्ष निकालता है कि विशिष्ट रेत के आकार के लिए, हाइड्रोडायनामिक बलों की तुलना में DLVO बल नगण्य हैं। एक क्रिटिकल वेलोसिटी ($U_{cr}$) व्युत्पन्न की गई है; रेत का उत्पादन केवल तभी होता है जब तरल का वेग इस दहलीज (threshold) को पार कर जाता है, जिससे टॉर्क बैलेंस बाधित होता है और कण अलग हो जाता है।
2. गणितीय सूत्रीकरण (Mathematical Formulation):
   • रैखिक प्रवाह (1D): फ्लुइडाइज्ड, अटैच्ड और री-अडेयर्ड कणों के मास बैलेंस के लिए गवर्निंग समीकरण स्थापित किए गए हैं। मॉडल तीन प्रमुख ट्यूनिंग पैरामीटर पेश करता है: एक फिल्ट्रेशन गुणांक ($\lambda_1$), एक प्रथम-इरोजन गुणांक ($\lambda_2$) जो सस्पेंडेड पार्टिकल फ्लक्स द्वारा संचालित होता है, और एक द्वितीय-इरोजन गुणांक ($\lambda_3$) जो तरल प्रवाह वेग द्वारा संचालित होता है। पोरोसिटी (porosity) परिवर्तनों को एक संशोधित वर्दालकिस (Vardoulakis) समीकरण के माध्यम से इरोजन से जोड़ा गया है, और पारगम्यता (permeability) को कारमन-कोज़ेनी (Carman-Kozeny) समीकरण का उपयोग करके अपडेट किया गया है।
   • रेडियल प्रवाह (अक्षीय सममित): वेलबोर से दूरी के साथ तरल वेग के परिवर्तन को ध्यान में रखते हुए, उत्पादन कुएं का अनुकरण करने के लिए रैखिक समीकरणों को रेडियल निर्देशांकों के लिए अनुकूलित किया गया है।
3. संख्यात्मक समाधान (Numerical Solution): कण सांद्रता और पोरोसिटी पर समय के साथ गवर्निंग समीकरणों को हल करने के लिए एक एक्सप्लिसिट फाइनाइट डिफरेंस स्कीम का उपयोग किया जाता है। मॉडल इनकंप्रेसिबल फ्लो, समान कण आकार, और तरल एवं सस्पेंडेड कणों के समान वेग को मानता है।
4. कैलिब्रेशन और सत्यापन: मॉडल मापदंडों ($\lambda_1, \lambda_2, \lambda_3$) को MATLAB में लीस्ट-स्क्वायर विधि का उपयोग करके पापामिचोस एट अल. [40] के प्रयोगात्मक डेटा के साथ मिलान करके कैलिब्रेट किया गया था।
5. संवेदनशीलता विश्लेषण (Sensitivity Analysis): संचयी रेत उत्पादन, कण सांद्रता, पोरोसिटी, पारगम्यता और दबाव गिरावट (pressure drop) पर प्रभाव देखने के लिए आयतनी प्रवाह दर और तीन इरोजन/फिल्ट्रेशन गुणांकों को बदलकर मॉडल के व्यवहार का आकलन किया गया।

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

• नवीन सैंडिंग मानदंड: हाइड्रोडायनामिक और इंटर-पार्टिकल बलों के टॉर्क बैलेंस से व्युत्पन्न एक क्रिटिकल वेलोसिटी फॉर्मूला का परिचय, जो कण पृथक्करण के लिए एक भौतिक दहलीज प्रदान करता है।
• एकीकृत संख्यात्मक मॉडल: एक एकीकृत मॉडल का विकास जो कण पृथक्करण, परिवहन और पुन: आसंजन को ध्यान में रखता है, और पोरोसिटी एवं पारगम्यता के गतिशील अपडेट के साथ जुड़ा हुआ है।
• दोहरी-ज्यामिति दृष्टिकोण: प्रयोगशाला डेटा को फील्ड-स्केल वेल भविष्यवाणियों में अनुवादित करने की अनुमति देने के लिए दोनों 1D रैखिक और अक्षीय सममित रेडियल समाधान प्रस्तुत किए गए।
• पैरामीटर पहचान: प्रयोगात्मक डेटा से मेल खाने के लिए मॉडल गुणांकों का सफल ट्यूनिंग, जो प्रयोगात्मक रेत उत्पादन गतिकी को दोहराने की मॉडल की क्षमता को प्रदर्शित करता है।

परिणाम

• मॉडल की सटीकता: संख्यात्मक मॉडल ने प्रयोगात्मक डेटा के साथ उच्च सहमति दिखाई, दो अलग-अलग परीक्षण मामलों के लिए $R^2$ मान क्रमशः 0.9735 और 0.9961 प्राप्त किए।
• प्रवाह दर का प्रभाव: संवेदनशीलता विश्लेषण से पता चला कि उत्पादन दर बढ़ाने से रेत उत्पादन की शुरुआत और दर तेज हो जाती है, लेकिन कुल संचयी रेत उत्पादन निर्माण की क्षमता द्वारा सीमित होता है। एक बार जब क्षरण योग्य सामग्री समाप्त हो जाती है या वेग पोरोसिटी परिवर्तनों के कारण क्रिटिकल स्तर से नीचे गिर जाता है, तो उत्पादन स्थिर हो जाता है।
• निकट-वेलबोर प्रभाव (Near-Wellbore Effects): रेडियल मॉडल में, रेत का क्षरण वेलबोर के पास सबसे गंभीर पाया गया जहाँ तरल वेग उच्चतम होता है। क्षरण क्षेत्र वेलबोर से कुछ मीटर तक ही सीमित था।
• गुणों में परिवर्तन: रेत के क्षरण के कारण निकट-वेलबोर क्षेत्र में पोरोसिटी और पारगम्यता में महत्वपूर्ण वृद्धि देखी गई।
• दबाव गिरावट (Pressure Drop): अध्ययन में पाया गया कि रेत उत्पादन विशिष्ट मॉडल गुणांकों के लिए निकट-वेलबोर क्षेत्र में दबाव गिरावट को लगभग 36% तक कम कर सकता है। यह कमी पारगम्यता में वृद्धि के कारण होती है, जो सस्पेंडेड कणों के कारण होने वाली घनत्व वृद्धि से अधिक प्रभावी है।
• शमन (Mitigation): परिणाम संकेत देते हैं कि तरल वेग को क्रिटिकल डिटैचमेंट थ्रेशोल्ड से नीचे रखने के लिए उत्पादन दरों को कम करके रेत उत्पादन को महत्वपूर्ण रूप से कम किया जा सकता है।

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

यह अध्ययन एक मध्यम दायरे को बनाए रखता है, यह नोट करते हुए कि मॉडल वर्महोल (wormhole) विकास जैसे जटिल भौतिकी को ध्यान में नहीं रखता है। इसके अलावा, लेखक सुझाव देते हैं कि हालांकि वर्तमान एक्सप्लिसिट संख्यात्मक समाधान अच्छी सहमति प्रदान करता है, भविष्य के कार्य में सटीकता को और बेहतर बनाने के लिए इंप्लिसिट या विश्लेषणात्मक समाधान विकसित करने पर ध्यान केंद्रित किया जा सकता है। शोध पत्र निष्कर्ष निकालता है कि उत्पादन दरों को कम करना रेत उत्पादन को कम करने के लिए एक व्यवहार्य रणनीति है, जिसका समर्थन मॉडल द्वारा तरल वेग और कण पृथक्करण के बीच संबंध के प्रदर्शन से होता है।