2.5-D Electrical Resistivity Forward Modelling with Undulating Topography using a Modified Half-Space Analytical Solution

यहाँ शोध पत्र "2.5-D Electrical Resistivity Forward Modelling with Undulating Topography Using a Modified Half-Space Analytical Solution" का विस्तृत तकनीकी सारांश दिया गया है।

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

डायरेक्ट करंट (DC) रेसिस्टिविटी इमेजिंग का व्यापक रूप से उपसतह प्रोफाइलिंग के लिए उपयोग किया जाता है, लेकिन यह अक्सर लहरदार स्थलाकृति (undulating terrain) पर आयोजित की जाती है। विश्वसनीय व्युत्क्रमण (inversion) के लिए सटीक फॉरवर्ड मॉडलिंग आवश्यक है, फिर भी जटिल स्थलाकृति को शामिल करना महत्वपूर्ण संख्यात्मक चुनौतियाँ प्रस्तुत करता है:

• विलक्षणता हटाना (Singularity Removal): मानक संख्यात्मक विधियाँ पॉइंट सोर्स इलेक्ट्रोड पर इलेक्ट्रिक पोटेंशियल की विलक्षणता (singularity) के साथ संघर्ष करती हैं। "प्राइमरी-सेकेंडरी" पोटेंशियल स्प्लिट तकनीक का आमतौर पर उपयोग किया जाता है, जहाँ विलक्षण प्राइमरी पोटेंशियल को विश्लेलेषणात्मक (analytically) रूप से हल किया जाता है, और सुचारू सेकेंडरी पोटेंशियल को संख्यात्मक रूप से हल किया जाता है।
• मौजूदा विधियों की सीमाएँ: वर्तमान दृष्टिकोण (जैसे, Penz et al., 2013) एक समतल सजातीय अर्ध-स्थान (flat homogeneous half-space) के लिए व्युत्पन्न एनालिटिकल प्राइमरी पोटेंशियल पर निर्भर करते हैं। यह धारणा केवल तभी मान्य है जब स्रोत के पास का भूभाग चिकना और समतल हो।
• मुख्य मुद्दा: प्राकृतिक भूभाग में, जहाँ तीखे कोने, उच्च वक्रता (curvature), या ढलान विच्छेदन (slope discontinuities) होते हैं, या यहाँ तक कि जब चिकने भूभाग को लीनियर फिनेट एलिमेंट्स (linear finite elements) का उपयोग करके विविक्त (discretized) किया जाता है, तो फ्लैट-सरफेस धारणा विफल हो जाती है।
  • लीनियर तत्व घुमावदार सतहों को पीसवाइस लीनियर सेगमेंट के रूप में अनुमानित करते हैं, जिससे स्रोत पर कृत्रिम "V-आकार" के कोने बन जाते हैं।
  • फ्लैट-हाफ-स्पेस एनालिटिकल समाधान इन V-आकार की ज्यामिति पर न्यूमैन सीमा स्थिति (Neumann boundary condition - शून्य फ्लक्स) को संतुष्ट नहीं करता है।
  • यह ज्यामितीय बेमेल (geometric mismatch) महत्वपूर्ण मॉडलिंग त्रुटियों (अक्सर >1%) की ओर ले जाता है, जो मेश रिफाइनमेंट या भूभाग के ज्यामितीय स्मूथिंग के बावजूद बनी रहती हैं।

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

लेखक एक संशोधित विलक्षणता हटाने की रणनीति (modified singularity removal strategy) प्रस्तावित करते हैं जो स्रोत की स्थानीय ज्यामिति के अनुरूप एक नया एनालिटिकल प्राइमरी पोटेंशियल व्युत्पन्न करती है।

• शासी समीकरण (Governing Equations): यह अध्ययन DC रेसिस्टिविटी के 2.5-D सूत्रीकरण का उपयोग करता है, जो 3-D पॉइसन समीकरण को 2-D समस्या (x-z प्लेन में) में बदलने के लिए फूरियर कोसाइन ट्रांसफॉर्म लागू करता है।
• प्राइमरी-सेकेंडरी स्प्लिट: कुल पोटेंशियल $u$ को एक प्राइमरी पोटेंशियल ($u_p$) और एक सेकेंडरी पोटेंशियल ($u_s$) में विभाजित किया जाता है।
• संशोधित प्राइमरी पोटेंशियल का व्युत्पत्ति:
  • फ्लैट हाफ-स्पेस (ठोस कोण $S = 2\pi$) मानने के बजाय, लेखक एक V-आकार के वेज (V-shaped wedge) (एक स्फेरिकल ल्यून) के शीर्ष पर स्थित पॉइंट सोर्स के लिए एक नया एनालिटिकल समाधान व्युत्पन्न करते हैं।
  • यह समाधान स्रोत पर स्थानीय स्थलाकृति द्वारा बनाए गए ठोस कोण (solid angle, $S$) को ध्यान में रखता है।
  • नया एनालिटिकल अभिव्यक्ति वाला प्राइमरी पोटेंशियल है:
    $$u_p(r) = \frac{I}{\sigma_0 S r}$$
    2.5-D फूरियर डोमेन में:
    $$\tilde{u}_p(x, k_y, z) = \frac{I}{\sigma_0 S} K_0(k_y r_{xz})$$
    जहाँ $S = 4\pi - 2\beta$ (जहाँ $\beta$ वेज एंगल है), और एक फ्लैट सतह के लिए, $S$ पुनः $2\pi$ हो जाता है।
• कार्यान्वयन (Implementation):
  • संशोधित प्राइमरी पोटेंशियल का उपयोग सेकेंडरी पोटेंशियल के लिए सीमा स्थितियों (boundary conditions) को परिभाषित करने के लिए किया जाता है।
  • इसे लीनियर (Q4) और उच्च-क्रम (Q9, Q16) तत्वों का उपयोग करते हुए एक फिनेट एलिमेंट मेथड (FEM) ढांचे के भीतर लागू किया गया है।
  • ठोस कोण $S$ की गणना स्रोत से जुड़े तत्वों की ज्यामिति के आधार पर स्पष्ट रूप से की जाती है, जो लीनियर और कर्व्ड विविक्तीकरण (discretization) दोनों को समायोजित करती है।

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

1. ज्यामितीय बेमेल की पहचान: यह शोध पत्र कठोरता से प्रदर्शित करता है कि लहरदार स्थलाकृति में मॉडलिंग त्रुटियाँ, फ्लैट-सरफेस एनालिटिकल प्राइमरी फील्ड और वास्तविक ठोस कोण (जो विविक्त स्थलाकृति द्वारा बनाया गया है) के बीच बेमेल के कारण होती हैं (विशेष रूप से लीनियर तत्वों के साथ)।
2. नवीन एनालिटिकल सूत्रीकरण: एक पॉइंट सोर्स के लिए V-आकार के वेज में प्राइमरी पोटेंशियल के लिए एक क्लोज्ड-फॉर्म एनालिटिकल समाधान का व्युत्पत्ति, जो मानक हाफ-स्पेस समाधान का सामान्यीकरण करता है।
3. लीनियर तत्वों के साथ मजबूती: यह विधि कोर्स लीनियर फिनेट एलीमेंट मेश के उपयोग की अनुमति देती है, जिसके लिए कृत्रिम स्थलाकृति स्मूथिंग या स्रोत के पास अत्यधिक मेश रिफाइनमेंट की आवश्यकता नहीं होती है, जो मौजूदा विधियों में कम्प्यूटेशनल रूप से महंगे समाधान हैं।
4. त्रुटि स्वतंत्रता: सैद्धांतिक विश्लेषण से पता चलता है कि पारंपरिक फ्लैट-सरफेस दृष्टिकोण में त्रुटि पूरी तरह से वास्तविक ठोस कोण और $2\pi$ के अनुपात ($|1 - S/2\pi|$) पर निर्भर करती है, जिससे त्रुटि स्रोत से दूरी के बावजूद स्थिर रहती है।

4. परिणाम (Results)

प्रस्तावित विधि को फ्लैट, V-आकार की ट्रेंच (trench), और साइनुसोइडल हिल-वैली (hill-valley) मॉडलों पर संख्यात्मक प्रयोगों के माध्यम से सत्यापित किया गया था:

• फ्लैट सजातीय और स्तरित मॉडल (Flat Homogeneous & Layered Models): इस विधि ने एनालिटिकल समाधानों के विरुद्ध < 0.1% सापेक्ष त्रुटि प्राप्त की, जो कार्यान्वयन की शुद्धता की पुष्टि करती है।
• V-आकार की ट्रेंच मॉडल:
  • पारंपरिक विधि (फ्लैट धारणा): स्रोत पर तीव्र ढलान विच्छेदन के कारण ~16.5% की त्रुटि उत्पन्न हुई।
  • प्रस्तावित विधि: मोटे मेश (1 मीटर तत्व आकार) और लीनियर तत्वों के साथ भी < 0.1% की त्रुटि प्राप्त की।
  • त्रुटि पूरे डोमेन में स्थिर रही, जो इसके सैद्धांतिक व्युत्पत्ति की पुष्टि करती है।
• साइनुसोइडल हिल-वैली मॉडल:
  • पारंपरिक विधि: लीनियर तत्वों के साथ त्रुटि 1% से 2% के बीच थी और ढलान जंक्शन पर स्रोत होने पर, उच्च-क्रम तत्वों के साथ भी, यह 10% तक पहुँच गई।
  • प्रस्तावित विधि: सभी मेश प्रकारों (लीनियर और कर्व्ड) और स्रोत स्थानों (वैली, फ्लैट, और स्लोप जंक्शन) में लगातार ≤ 0.1% की त्रुटि प्राप्त की।
  • विशेष रूप से, प्रस्तावित विधि ने कुछ परिदृश्यों में कर्व्ड हायर-ऑर्डर तत्वों के साथ पारंपरिक विधि को भी पीछे छोड़ दिया।

5. महत्व और निष्कर्ष (Significance and Conclusion)

• कम्प्यूटेशनल दक्षता: स्रोत के पास घने मेश रिफाइनमेंट या मैनुअल टेरेन स्मूथिंग की आवश्यकता को समाप्त करके, यह विधि फॉरवर्ड मॉडलिंग और व्युत्क्रमण (inversion) की कम्प्यूटेशनल लागत को काफी कम कर देती है।
• जटिल स्थलाकृति में सटीकता: यह उन क्षेत्रों में DC रेसिस्टिविटी मॉडलिंग के लिए एक मजबूत समाधान प्रदान करता है जहाँ तीखे टोपोग्राफिक फीचर्स (जैसे चट्टानें, खाइयाँ, या कटक) होते हैं जहाँ पिछली विधियाँ विफल रही थीं।
• व्यावहारिक अनुप्रयोग: उच्च सटीकता बनाए रखते हुए लीनियर तत्वों (जो गणना में तेज़ और बनाने में आसान हैं) का उपयोग करने की क्षमता इस दृष्टिकोण को बड़े पैमाने के 2.5-D व्युत्क्रमण ढांचों के लिए अत्यधिक उपयुक्त बनाती है।
• भावी कार्य: लेखक इस सॉलिड-एंगल-आधारित दृष्टिकोण को पूर्ण 3-D फॉरवर्ड मॉडलिंग तक विस्तारित करने का सुझाव देते हैं।

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