भूजल प्रवाह समीकरण

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

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

द्रव्यमान संतुलन

क्षणिक भूजल प्रवाह समीकरण पर पहुंचने के लिए बड़े पैमाने पर संतुलन किया जाता है और इसे डार्सी के नियम के साथ प्रयोग किया जाना चाहिए। यह संतुलन ऊष्मा समीकरण में आने के लिए ऊष्मा हस्तांतरण में प्रयुक्त ऊर्जा संतुलन के अनुरूप होता है। यह मात्र लेखांकन का एक बयान है, कि किसी दिए गए नियंत्रण मात्रा के लिए स्रोतों या सिंक के अतिरिक्त द्रव्यमान को बनाया या नष्ट किया जा सकता है। द्रव्यमान के संरक्षण में कहा गया है कि समय की एक निश्चित वृद्धि (Δt) के लिए सीमाओं के पार बहने वाले द्रव्यमान और आयतन के भीतर के स्रोतों के बीच का अंतर भंडारण में परिवर्तन होता है। जिसे इस रूप में दिखाया जाता है,

${\frac {\Delta M_{stor}}{\Delta t}}={\frac {M_{in}}{\Delta t}}-{\frac {M_{out}}{\Delta t}}-{\frac {M_{gen}}{\Delta t}}$

प्रसार समीकरण (क्षणिक प्रवाह)

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

$S_{s}{\frac {\partial h}{\partial t}}=-\nabla \cdot q-G.$

इसे अन्य क्षेत्रों में प्रसार समीकरण या ऊष्मा समीकरण के रूप में जाना जाता है, यह एक परवलयिक आंशिक अंतर समीकरण (पीडीई) के रूप में होता है। यह गणितीय कथन इंगित करता है कि बायीं ओर समय के साथ हाइड्रोलिक हेड में परिवर्तन फ्लक्स के नकारात्मक विचलन के बराबर होता है और स्रोत शर्तों से इस समीकरण में हेड और फ्लक्स अज्ञात रूप में होते हैं, लेकिन डार्सी का नियम फ्लक्स को हाइड्रोलिक हेड्स से संबंधित होता है, इसलिए इसे फ्लक्स (q) के लिए प्रतिस्थापित करने से होता है

$S_{s}{\frac {\partial h}{\partial t}}=-\nabla \cdot (-K\nabla h)-G.$

अब अगर हाइड्रोलिक चालकता (K) स्थानिक रूप से एकसमान है और टेन्सर के अतिरिक्त आइसोट्रोपिक है, तो इसे स्थानिक व्युत्पन्न से बाहर निकाला जा सकता है, जिससे उन्हें लाप्लासियन में सरल बनाया जा सके, यह समीकरण बनाता है।

$S_{s}{\frac {\partial h}{\partial t}}=K\nabla ^{2}h-G.$

विशिष्ट भंडारण (S_s) द्वारा विभाजित करके, दाहिनी ओर हाइड्रोलिक विसरण (α = K/S_sया समकक्ष, α = T/S) के रूप में होता है। हाइड्रोलिक विसरण उस गति के समानुपाती होती है जिस पर एक परिमित दबाव पल्स प्रणाली के माध्यम से α के बड़े मान संकेतों के तेजी से प्रसार के लिए प्रसारित होता है और इस प्रकार भूजल प्रवाह समीकरण बन जाता है।

${\frac {\partial h}{\partial t}}=\alpha \nabla ^{2}h-G.$

जहां सिंक/स्रोत शब्द G, में अब समान इकाइयों के रूप में हैं, लेकिन उपयुक्त भंडारण अवधि से विभाजित है जैसा कि हाइड्रोलिक विसरण प्रतिस्थापन द्वारा परिभाषित किया गया है।

आयताकार कार्टेसियन निर्देशांक

मॉडफ्लो में प्रयुक्त त्रि-आयामी परिमित अंतर ग्रिड

विशेष रूप से आयताकार ग्रिड परिमित अंतर मॉडल का उपयोग करते है उदाहरण के लिए यूएसजीएस द्वारा बनाए गए मॉडफ्लो कार्टेशियन निर्देशांक का वर्णन करते है। इन निर्देशांकों में सामान्य लाप्लासियन ऑपरेटर विशेष रूप से तीन आयामी प्रवाह के लिए बन जाता है।

${\frac {\partial h}{\partial t}}=\alpha \left[{\frac {\partial ^{2}h}{\partial x^{2}}}+{\frac {\partial ^{2}h}{\partial y^{2}}}+{\frac {\partial ^{2}h}{\partial z^{2}}}\right]-G.$

मॉडफ्लो कोड गवर्निंग ग्राउंडवाटर फ्लो इक्वेशन के एक ओर्थोगोनल 3-डी फॉर्म को अलग करता है और अनुकरण करता है। चूँकि, अगर उपयोगकर्ता ऐसा करना चाहता है तो उसके पास अर्ध-3D मोड में चलने का विकल्प होता है; इस स्थिति में नमूना k और S_s के अतिरिक्त लंबवत औसत T और S से संबंधित होता है। अर्ध-3डी मोड में रिसाव की अवधारणा का उपयोग करके 2डी क्षैतिज परतों के बीच प्रवाह की गणना की जाती है।

परिपत्र बेलनाकार निर्देशांक

एक अन्य उपयोगी समन्वय प्रणाली 3डी बेलनाकार निर्देशांक के रूप में है, सामान्यतः जहां एक पंपिंग कुआं Z अक्ष के समानांतर मूल पर स्थित एक लाइन स्रोत के रूप में होता है, जिससे अभिसरण रेडियल प्रवाह होता है। इन शर्तों के अनुसार उपरोक्त समीकरण r रेडियल दूरी और θ कोण के रूप में बन जाता है।

${\frac {\partial h}{\partial t}}=\alpha \left[{\frac {\partial ^{2}h}{\partial r^{2}}}+{\frac {1}{r}}{\frac {\partial h}{\partial r}}+{\frac {1}{r^{2}}}{\frac {\partial ^{2}h}{\partial \theta ^{2}}}+{\frac {\partial ^{2}h}{\partial z^{2}}}\right]-G.$

अनुमान

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

• एक्वीफर सामग्री असंपीड्य है मैट्रिक्स में कोई बदलाव नहीं है दबाव उर्फ ​​अवतलन में परिवर्तन के कारण होते है
• पानी निरंतर घनत्व असंपीड्य के रूप में होते है
• एक्वीफर पर कोई बाहरी भार जैसे, ओवरबर्डन वायुमंडलीय दबाव स्थिर रूप में होते है
• 1डी रेडियल समस्या के लिए पम्पिंग कुआँ पूरी तरह से एक गैर रिसाव वाले एक्वीफर में प्रवेश के रूप में है
• भूजल धीरे-धीरे बह रहा है और इस प्रकार रेनॉल्ड्स संख्या से कम होता है
• हाइड्रोलिक चालकता (k) एक समदैशिक अदिश भौतिकी के रूप में है

इन बड़ी धारणाओं के अतिरिक्त भूजल प्रवाह समीकरण स्रोतों और सिंक के क्षणिक वितरण के कारण एक्वीफर में प्रमुखों के वितरण का प्रतिनिधित्व करने का अच्छा काम करता है।

लाप्लास समीकरण (स्थिर अवस्था प्रवाह)

अगर एक्विफायर में रिचार्जिंग सीमा की स्थितियां हैं तो एक स्थिर स्थिति तक पहुंचा जा सकता है या इसे कई स्थिति में अनुमान के रूप में उपयोग किया जा सकता है और प्रसार समीकरण लाप्लास समीकरण को सरल करता है।

$0=\alpha \nabla ^{2}h$

यह समीकरण बताता है कि हाइड्रोलिक हेड एक हार्मोनिक फलन है और अन्य क्षेत्रों में इसके कई एनालॉग हैं। लाप्लास समीकरण को प्रद्यौगिकीय का उपयोग करके हल किया जा सकता है, ऊपर बताई गई समान मान्यताओं का उपयोग करते हुए, लेकिन एक स्थिर अवस्था प्रवाह क्षेत्र की अतिरिक्त आवश्यकताओं के रूप में होती है।

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

एक पम्पिंग कुएं में स्थिर-अवस्था का प्रवाह जो वास्तव में कभी नहीं होता है, लेकिन कभी-कभी एक उपयोगी सन्निकटन के रूप में होता है जिसे सामान्यता थिएम समाधान कहा जाता है।

द्वि-आयामी भूजल प्रवाह

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

डुपिट-फोर्चहाइमर धारणा को लागू करके भूजल प्रवाह समीकरण का एक वैकल्पिक सूत्रीकरण प्राप्त किया जा सकता है, जहां यह माना जाता है कि शीर्ष ऊर्ध्वाधर दिशा में भिन्न नहीं होते हैं (अर्थात, $\partial h/\partial z=0$). एक क्षैतिज जल संतुलन क्षेत्र के साथ एक लंबे ऊर्ध्वाधर स्तंभ पर लागू होता है $\delta x\delta y$ एक्वीफर आधार से असंतृप्त सतह तक विस्तार होता है। इस दूरी को संतृप्त मोटाई, b के रूप में जाना जाता है। एक सीमित एक्वीफर में, संतृप्त मोटाई एक्वीफर H की ऊंचाई से निर्धारित होती है और दाब शीर्ष, हर जगह गैर-शून्य होता है। एक असीमित एक्वीफर में, संतृप्त मोटाई को जल सारिणी की सतह और एक्वीफर आधार के बीच ऊर्ध्वाधर दूरी के रूप में परिभाषित किया जाता है। यदि $\partial h/\partial z=0$ और एक्वीफर आधार शून्य आधार पर है, तो असंबद्ध संतृप्त मोटाई शीर्ष के बराबर है, अर्थात, b=h।

हाइड्रोलिक चालकता और प्रवाह के क्षैतिज घटकों दोनों को मानते हुए एक्विफायर की संपूर्ण संतृप्त मोटाई के साथ समान रूप में होती है (अर्थात, $\partial q_{x}/\partial z=0$ और $\partial K/\partial z=0$), हम एकीकृत भूजल निर्वहन,Q_x और Q_y:के संदर्भ में डार्सी के नियम को व्यक्त करते हैं,

$Q_{x}=\int _{0}^{b}q_{x}dz=-Kb{\frac {\partial h}{\partial x}}$

$Q_{y}=\int _{0}^{b}q_{y}dz=-Kb{\frac {\partial h}{\partial y}}$

यदि द्रव्यमान संतुलन अभिव्यक्ति के रूप में सम्मिलित होते है, तो असम्पीडित संतृप्त भूजल प्रवाह के लिए सामान्य 2D के रूप में समीकरण प्राप्त करते हैं

${\frac {\partial nb}{\partial t}}=\nabla \cdot (Kb\nabla h)+N.$

जहाँ n एक्विफायर सरंध्रता के रूप में होते है। स्रोत शब्द, N लंबाई प्रति समय ऊर्ध्वाधर दिशा में पानी के अतिरिक्त जैसे, पुनर्भरण का प्रतिनिधित्व करता है और इस प्रकार संतृप्त मोटाई, विशिष्ट भंडारण और विशिष्ट परिणाम के लिए सही परिभाषाओं को सम्मिलित करते है, और इसे सीमित और अपरिमित स्थितियों के लिए दो अद्वितीय रूप के समीकरणों में बदल सकते हैं

$S{\frac {\partial h}{\partial t}}=\nabla \cdot (Kb\nabla h)+N.$

(सीमित), जहां S=S_sb एक्वीफर भंडारण के रूप में होते है

$S_{y}{\frac {\partial h}{\partial t}}=\nabla \cdot (Kh\nabla h)+N.$

(अपरिबद्ध), जहां S_y एक्विफायर की विशिष्ट रूप में होते है।

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

$\nabla \cdot (K\nabla h^{2})=-2N.$

या, सजातीय जलवाही स्तर के रूप में है,

$\nabla ^{2}h^{2}=-{\frac {2N}{K}}.$

यह फॉर्मूलेशन हमें असीमित प्रवाह के स्थिति में रैखिक पीडीई को हल करने के लिए मानक विधियों को लागू करने की अनुमति देता है। बिना पुनर्भरण वाले विषम एक्वीफर के लिए मिश्रित सीमित/अपरिबद्ध स्थिति के लिए संभावित प्रवाह विधियों को लागू किया जा सकता है।

यह भी देखें

• विश्लेषणात्मक तत्व विधि
  • आंशिक अंतर समीकरणों के लिए संख्यात्मक विधियों के लिए प्रयुक्त एक संख्यात्मक विधि के रूप में होती है
• डुपिट-फोर्चाइमर धारणा
  • ऊर्ध्वाधर प्रवाह के संबंध में भूजल प्रवाह समीकरण का सरलीकरण
• भूजल ऊर्जा संतुलन
  • ऊर्जा संतुलन पर आधारित भूजल प्रवाह समीकरण
• रिचर्ड्स समीकरण

अग्रिम पठन

• H. F. Wang and M.P. Anderson Introduction to Groundwater Modeling: Finite Difference and Finite Element Methods
  • An excellent beginner's read for groundwater modeling. Covers all the basic concepts, with simple examples in FORTRAN 77.
• Freeze, R. Allan; Cherry, John A. (1979). Groundwater. Prentice Hall. ISBN 978-0133653120.

बाहरी संबंध

• USGS groundwater software — free groundwater modeling software like MODFLOW
• Groundwater Hydrology (MIT OpenCourseware)