تقنيات معالجة وتفسير بيانات MT
بنداهولوان
تُعدّ طريقة المغناطيسية الأرضية (MT) تقنية جيوفيزيائية سلبية تستخدم التغيرات الطبيعية في المجالين الكهربائي والمغناطيسي للأرض لرسم خريطة توزيع المقاومة الكهربائية تحت السطحية. ترتبط المقاومة الكهربائية ارتباطًا وثيقًا بنوع الصخور، ومحتوى السوائل، ودرجة الحرارة، والتركيبات الجيولوجية كالفوالق ومناطق التحول. ولذلك، تُستخدم طريقة المغناطيسية الأرضية على نطاق واسع في استكشاف الطاقة الحرارية الأرضية، والتمعدن، والأحواض الرسوبية (الهيدروكربونات)، ودراسات التكتونيات والقشرة الأرضية.
مع ذلك، نادرًا ما تكون بيانات MT الخام قابلة للاستخدام مباشرةً. فهي تتطلب سلسلة من خطوات المعالجة لتحسين نسبة الإشارة إلى الضوضاء، تليها عملية تفسير تستند إلى نموذج مقاومة كهربائية متوافق جيولوجيًا. تتناول هذه المقالة تقنيات معالجة بيانات MT وتفسيرها، بدءًا من جمع البيانات ومراقبة جودتها، مرورًا بمعالجة المعاوقة، وصولًا إلى عملية الانعكاس والتفسير الجيولوجي.
-
قاعدة البيانات MT: المجال الكهربائي، المجال المغناطيسي، وموتر المعاوقة
في مسوحات MT، تسجل المجسات مركبات المجال الكهربائي (Ex، Ey) ومركبات المجال المغناطيسي (Hx، Hy، وأحيانًا Hz) كدالة للزمن. ومن العلاقة بين المجالين الكهربائي والمغناطيسي في مجال التردد، يُستخلص موتر المعاوقة.
\[
\begin{bmatrix} E_x \\ E_y \end{bmatrix}
=
\begin{bmatrix} Z_{xx} & Z_{xy} \\ Z_{yx} & Z_{yy} \end{bmatrix}
\begin{bmatrix} H_x \\ H_y \end{bmatrix}
\]
ثم تُستخلص من عناصر الموتر هذه معايير تفسيرية رئيسية مثل المقاومة الظاهرية (ρa) والطور (φ). عمومًا، تُظهر البيانات الجيدة اتجاهات سلسة ومتسقة في ρa و φ عبر الترددات، ولها هوامش خطأ واقعية.
-
المرحلة الأولى: مراقبة الجودة والمعالجة المسبقة
تبدأ عملية الترجمة الآلية الجيدة من مراقبة الجودة لأن بيانات الحقلتشمل هذه المرحلة ما يلي:
1. فحص السلاسل الزمنية: البحث عن الارتفاعات المفاجئة، والانحراف، وتشبع المستشعر، أو الاضطرابات الدورية.
2. الضوضاء الثقافية: مصادر الضوضاء الشائعة هي خطوط الطاقة بتردد 50/60 هرتز، والقطارات، والأنشطة الصناعية، والأسوار الكهربائية، والاتصالات السلكية واللاسلكية.
3. حالة الأقطاب الكهربائية والوصلات الأرضية: ستؤدي مقاومة التلامس العالية إلى تدهور جودة Ex/Ey.
4. اتجاه وموضع المستشعر: يمكن أن تتسبب أخطاء سمت المستشعر في تشوهات في التفسير، خاصة في دراسات البنية الاتجاهية.
في هذه المرحلة، يتم عادةً قطع أجزاء البيانات التالفة، ويتم إجراء تصحيح الإزاحة، ويتم إجراء مزامنة الوقت في حالة استخدام محطة مرجعية بعيدة.
-
المرحلة الثانية: تحويل المجال الترددي وتقدير الطيف
لأن تقنية MT تحلل العلاقة بين المجال الكهربائي (E) والمجال المغناطيسي (H) في مجال التردد، يتم تحويل السلسلة الزمنية باستخدام تقنيات طيفية مثل تحويل فورييه. ثم تُقسّم البيانات إلى نوافذ (قطاعات) لضمان الاستقرار الإحصائي. ومن التقنيات الشائعة الاستخدام ما يلي:
– استخدام تقنية التظليل والتضييق (مثل هانينج) لتقليل التسرب الطيفي.
– حساب المتوسط بين النوافذ للحصول على تقدير طيفي أكثر دقة.
– أطياف القدرة المتقاطعة لتحديد العلاقة بين E و H.
الهدف هو الحصول على تقديرات دقيقة للممانعة على نطاق تردد واسع، من الترددات العالية (الضحلة) إلى الترددات المنخفضة (العميقة).
-
المرحلة 3: المعالجة القوية والرجوع عن بعد
من التحديات التي تواجه معالجة البيانات الآلية وجود ضوضاء مترابطة في كثير من الأحيان ضمن قناة معينة. لذا، طُوّرت تقنيات معالجة قوية للحد من تأثير القيم الشاذة وقطاعات البيانات غير المتناسقة. وتستخدم هذه الطرق القوية ترجيحًا تكراريًا لكبح تأثير البيانات غير المتناسقة.
تُعدّ تقنية القياس المرجعي عن بُعد (RR) تقنيةً بالغة الأهمية. تقوم هذه التقنية على قياس المجال المغناطيسي في محطة بعيدة عن مصدر الضوضاء المحلي. ومن خلال ربط قيمتي E وH في المحطة الرئيسية بقيمة H في محطة القياس المرجعي، يُمكن تقليل تأثير الضوضاء المحلية على H. وتُعدّ تقنية القياس المرجعي عن بُعد فعّالة بشكل خاص في المناطق الصاخبة، مثل المناطق السكنية أو البنية التحتية.
الناتج من هذه المرحلة عادةً ما يكون كالتالي:
– منحنيات ρa و φ لـ Zxy و Zyx
– شريط الخطأ (الانحراف المعياري)
– التماسك بين الإشارات الأولية
-
المرحلة الرابعة: إزالة الضوضاء ومعالجة التشوه (التحول الثابت)
بالإضافة إلى التشويش، تتأثر تقنية قياس المقاومة المغناطيسية غالبًا بالانزياح الساكن، وهو انزياح رأسي في منحنى المقاومة الظاهرية نتيجةً لعدم تجانس سطحي (مثل طبقات الطين الرقيقة، أو الحصى، أو ظروف القطب الكهربائي). لا يُغير الانزياح الساكن الطور بشكل ملحوظ، ولكنه يُزيح قيمة المقاومة الظاهرية (ρa) لأعلى أو لأسفل بشكل مضاعف.
يمكن إجراء معالجة الإزاحة الثابتة من خلال:
1. تصحيح قائم على بيانات TDEM/CSAMT كتحكم في المقاومة الضحلة.
2. الانعكاس المشترك لنموذج MT-TDEM لجعل النموذج الضحل أكثر تماسكًا.
3. نهج عكسي قوي للتحول، على سبيل المثال السماح بمعاملات التحول لكل محطة.
بالإضافة إلى ذلك، يتم أيضًا إجراء عملية الكشف عن القيم الشاذة للتردد: عادةً ما يتم تجاهل نقاط معينة على المنحنى تنحرف بشكل حاد أو يتم إعطاؤها خطأ أكبر.
-
المرحلة الخامسة: تحليل الأبعاد والضربة
قبل إجراء عملية الانعكاس، من المهم تحديد ما إذا كان التركيب تحت السطحي أحادي البعد أو ثنائي البعد أو ثلاثي الأبعاد. يؤثر هذا على اختيار طريقة الانعكاس وتفسير النتائج. تشمل التحليلات الشائعة ما يلي:
– معلمات الانحراف (مثل انحراف باهر) لتقييم المستوى ثلاثي الأبعاد.
– موتر الطور لرؤية الاتجاه السائد للهيكل دون التأثر بالانزياح الثابت.
– تحليل اتجاه الضربة لتحديد اتجاه الضربة السائد للهيكل ثنائي الأبعاد.
إذا أظهرت البيانات خصائص ثنائية الأبعاد قوية، يُجرى عادةً تدوير موتر المعاوقة باتجاه الانحدار لتوضيح المكونات الرئيسية (TE/TM). أما إذا كانت الخصائص ثلاثية الأبعاد قوية، فيُفضّل إجراء عملية عكس ثلاثية الأبعاد.
-
الخطوة 6: عكس بيانات MT (أحادي البعد، ثنائي البعد، ثلاثي الأبعاد)
يُجرى التفسير الكمي لبيانات MT عادةً من خلال عملية الانعكاس، والتي تتضمن إيجاد نموذج المقاومة الكهربائية الذي يُطابق استجابة البيانات على أفضل وجه. يُعد انعكاس MT عملية غير خطية وغير محددة جيدًا، مما يتطلب تنظيمًا لمنع النموذج من أن يصبح غير مستقر. أنواع الانعكاس:
1. الانعكاس أحادي البعد: مناسب للمناطق الطبقية (مثل الرواسب الأفقية). سريع، ولكنه محدود.
2. الانعكاس ثنائي الأبعاد: مناسب للهياكل الممتد مثل الصدوع أو الخنادق أو الأنظمة الحرارية الأرضية الاتجاهية.
3. الانعكاس ثلاثي الأبعاد: الأكثر واقعية للجيولوجيا المعقدة، ولكنه يتطلب بيانات كثيفة وحسابات كبيرة ومراقبة جودة صارمة.
تتضمن دالة الهدف العكسية عادةً عدم تطابق البيانات وخشونة النموذج:
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
حيث λ هو معامل التنظيم. يُعد اختيار قيمة λ أمرًا بالغ الأهمية: فالقيمة الصغيرة جدًا تجعل النموذج خشنًا جدًا (مُفرط التخصيص)، والقيمة الكبيرة جدًا تجعل النموذج سلسًا جدًا ويتجاهل السمات الجيولوجية المهمة.
-
الخطوة 7: تقييم نتائج الانعكاس والحساسية
بعد عملية الانعكاس، يجب تقييم النموذج، ليس فقط بصريًا. يشمل هذا التقييم ما يلي:
– عدم التطابق (RMS): ما إذا كان يتطابق مع الهدف (على سبيل المثال RMS ~ 1-2 اعتمادًا على تعريف الخطأ).
– مقارنة منحنيات البيانات مقابل استجابة النموذج في كل محطة.
– اختبار الدقة: على سبيل المثال اختبار رقعة الشطرنج أو خريطة الحساسية.
– تأثير النماذج السابقة: جرب عدة نماذج أولية لمعرفة مدى استقرار النتائج.
إذا تغير النموذج بشكل جذري عند تغيير معلمة صغيرة، فهذا يعني أن التفسير يجب أن يكون أكثر دقة وقد يتطلب بيانات إضافية.
-
التفسير الجيولوجي: ربط المقاومة الكهربائية بالأنظمة تحت السطحية
لا تُعدّ المقاومة الكهربائية نوعاً مباشراً من أنواع الصخور، بل هي استجابة فيزيائية مركبة. ومع ذلك، غالباً ما تُستخدم بعض الأنماط العامة كإرشادات:
– مقاومة منخفضة (موصلة): الطين، والتغير الحراري المائي (غطاء الطين)، ومنطقة السائل الملحي، والجرافيت، أو معادن الكبريتيد.
– المقاومة العالية (المقاومة): الصخور النارية الضخمة، والصخور الجافة، والمناطق المتسيليكية، أو القاعدة البلورية.
في مجال استكشاف الطاقة الحرارية الأرضية، على سبيل المثال، غالباً ما تُظهر النماذج الكلاسيكية ما يلي:
1. غطاء من الطين الموصل في الأعلى
2. يكون الخزان أكثر مقاومة في الأسفل
3. منطقة التدفق الصاعد المتحكم بها هيكليًا (صدع)
4. مصادر الحرارة التي تظهر أحيانًا على شكل شذوذات مقاومة أو معقدة اعتمادًا على الليثولوجيا ودرجة الحرارة.
يرتبط التفسير الجيد دائمًا ببيانات أخرى: الجيولوجيا السطحية، والمظاهر الحرارية، والكيمياء الجيولوجية، والجاذبية، والبيانات الزلزالية، أو بيانات الآبار.
-
غطاء
تتداخل تقنيات معالجة بيانات MT وتفسيرها: بدءًا من مراقبة جودة السلاسل الزمنية، وتقدير الطيف، ومعالجة البيانات المرجعية القوية/البعيدة، وتصحيح الانزياح الثابت، وتحليل الأبعاد، وصولًا إلى الانعكاس ثنائي/ثلاثي الأبعاد وتقييم الدقة. ولا يعتمد نجاح MT على برنامج الانعكاس فحسب، بل أيضًا على جودة القياسات، وفهم التشويش، ودمجه مع السياق الجيولوجي.
بفضل سير العمل المنضبط والتفسير القائم على البيانات المتعددة، تعد تقنية MT أداة قوية لرسم خرائط هياكل المقاومة تحت السطحية والمساعدة في اتخاذ القرارات في الدراسات الاستكشافية والجيولوجية المعقدة.
