เทคนิคการประมวลผลและการตีความข้อมูล MT

เทคนิคการประมวลผลและการตีความข้อมูล MT

เพนดาฮูหวน
วิธีการแมกนีโตเทลลูริก (MT) เป็นเทคนิคทางธรณีฟิสิกส์แบบพาสซีฟที่ใช้ประโยชน์จากการเปลี่ยนแปลงตามธรรมชาติของสนามไฟฟ้าและสนามแม่เหล็กโลกเพื่อสร้างแผนที่การกระจายตัวของความต้านทานใต้ผิวดิน ความต้านทานมีความสัมพันธ์อย่างใกล้ชิดกับชนิดของหิน ปริมาณของเหลว อุณหภูมิ และโครงสร้างทางธรณีวิทยา เช่น รอยเลื่อนหรือเขตการเปลี่ยนแปลง ดังนั้น MT จึงถูกนำมาใช้กันอย่างแพร่หลายในการสำรวจพลังงานความร้อนใต้พิภพ การเกิดแร่ แอ่งตะกอน (ไฮโดรคาร์บอน) และการศึกษาทางด้านธรณีแปรสัณฐานและเปลือกโลก

อย่างไรก็ตาม ข้อมูล MT ดิบนั้นแทบจะไม่สามารถนำไปใช้ได้โดยตรง จำเป็นต้องผ่านกระบวนการประมวลผลหลายขั้นตอนเพื่อปรับปรุงอัตราส่วนสัญญาณต่อสัญญาณรบกวน จากนั้นจึงตีความโดยใช้แบบจำลองความต้านทานที่สอดคล้องกับทางธรณีวิทยา บทความนี้จะกล่าวถึงเทคนิคการประมวลผลและการตีความข้อมูล MT ตั้งแต่การได้มาซึ่งข้อมูลและการควบคุมคุณภาพ ไปจนถึงการประมวลผลอิมพีแดนซ์ ผ่านการผกผันและการตีความทางธรณีวิทยา

-

ฐานข้อมูล MT: สนาม E, สนาม H และเทนเซอร์อิมพีแดนซ์
ในการสำรวจแบบ 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: การควบคุมคุณภาพ (QC) และการเตรียมการเบื้องต้น
กระบวนการประมวลผล MT ที่ดีเริ่มต้นด้วย QC ซึ่งเริ่มจากข้อมูลภาคสนาม ขั้นตอนนี้ประกอบด้วย:

1. การตรวจสอบอนุกรมเวลา: มองหาค่าผิดปกติ การเปลี่ยนแปลงอย่างฉับพลัน การอิ่มตัวของเซ็นเซอร์ หรือความผิดปกติเป็นระยะๆ
2. เสียงรบกวนจากกิจกรรมทางวัฒนธรรม: แหล่งกำเนิดเสียงรบกวนทั่วไป ได้แก่ สายส่งไฟฟ้าแรงสูง 50/60 เฮิรตซ์ รถไฟ กิจกรรมทางอุตสาหกรรม รั้วไฟฟ้า และการสื่อสารโทรคมนาคม
3. สภาพของอิเล็กโทรดและจุดต่อลงดิน: ความต้านทานการสัมผัสสูงจะทำให้คุณภาพของ Ex/Ey แย่ลง
4. การวางแนวและตำแหน่งของเซ็นเซอร์: ข้อผิดพลาดในการวางแนวของเซ็นเซอร์อาจทำให้เกิดความคลาดเคลื่อนในการตีความ โดยเฉพาะอย่างยิ่งในการศึกษาโครงสร้างแบบมีทิศทาง

ในขั้นตอนนี้ โดยปกติแล้วส่วนข้อมูลที่เสียหายจะถูกตัดออก มีการปรับแก้ค่าชดเชย และมีการซิงโครไนซ์เวลาหากมีการใช้สถานีอ้างอิงระยะไกล

-

ขั้นตอนที่ 2: การแปลงโดเมนความถี่และการประมาณค่าสเปกตรัม
เนื่องจาก MT วิเคราะห์ความสัมพันธ์ระหว่าง E และ H ในโดเมนความถี่ อนุกรมเวลาจึงถูกแปลงโดยใช้เทคนิคสเปกตรัม เช่น การแปลงฟูริเยร์ จากนั้นข้อมูลจะถูกแบ่งออกเป็นช่วง (เซกเมนต์) เพื่อให้มั่นใจถึงเสถียรภาพทางสถิติ เทคนิคที่ใช้กันทั่วไปบางส่วน ได้แก่:

– การกำหนดขอบเขตและการลดทอน (เช่น วิธี Hanning) เพื่อลดการรั่วไหลของสเปกตรัม
– การหาค่าเฉลี่ยระหว่างช่วงเวลาต่างๆ เพื่อให้ได้ค่าประมาณสเปกตรัมที่แม่นยำยิ่งขึ้น
– สเปกตรัมกำลังไขว้เพื่อสร้างความสัมพันธ์ระหว่าง E กับ H

เป้าหมายคือการประมาณค่าอิมพีแดนซ์ได้อย่างแม่นยำในช่วงความถี่กว้าง ตั้งแต่ความถี่สูง (ตื้น) ไปจนถึงความถี่ต่ำ (ลึก)

-

ขั้นตอนที่ 3: การประมวลผลที่แข็งแกร่งและการอ้างอิงระยะไกล
หนึ่งในความท้าทายของการแปลด้วยเครื่องคือสัญญาณรบกวนที่มีความสัมพันธ์กันภายในช่องสัญญาณเดียวกัน ดังนั้นจึงมีการพัฒนาเทคนิคการประมวลผลที่แข็งแกร่งเพื่อลดอิทธิพลของค่าผิดปกติและส่วนข้อมูลที่ไม่ดี วิธีการที่แข็งแกร่งเหล่านี้ใช้การถ่วงน้ำหนักแบบวนซ้ำเพื่อลดผลกระทบของข้อมูลที่ไม่สอดคล้องกัน

เทคนิคที่สำคัญอีกอย่างหนึ่งคือการอ้างอิงระยะไกล (Remote Referencing: RR) แนวคิดคือการวัดสนามแม่เหล็กที่สถานีซึ่งอยู่ห่างไกลจากแหล่งกำเนิดสัญญาณรบกวน โดยการหาความสัมพันธ์ระหว่างค่า E และ H ที่สถานีหลักกับค่า H ที่สถานี RR จะช่วยลดผลกระทบของสัญญาณรบกวนในพื้นที่ต่อค่า H ได้ เทคนิค RR มีประสิทธิภาพเป็นพิเศษในพื้นที่ที่มีสัญญาณรบกวนสูง เช่น บริเวณใกล้พื้นที่อยู่อาศัยหรือโครงสร้างพื้นฐาน

ผลลัพธ์จากขั้นตอนนี้โดยทั่วไปคือ:
– เส้นโค้ง ρa และ φ สำหรับ Zxy และ Zyx
– แถบแสดงความคลาดเคลื่อน (ค่าเบี่ยงเบนมาตรฐาน)
– ความสอดคล้องกันระหว่างสัญญาณหลัก

-

ขั้นตอนที่ 4: การลดเสียงรบกวนและการจัดการความผิดเพี้ยน (Static Shift)
นอกจากสัญญาณรบกวนแล้ว MT มักได้รับผลกระทบจากค่าเบี่ยงเบนคงที่ ซึ่งเป็นการเลื่อนในแนวตั้งของเส้นโค้งความต้านทานปรากฏเนื่องจากความไม่สม่ำเสมอในระดับตื้น (เช่น ชั้นดินเหนียวบางๆ กรวด หรือสภาพของอิเล็กโทรด) ค่าเบี่ยงเบนคงที่ไม่ได้เปลี่ยนแปลงเฟสอย่างมีนัยสำคัญ แต่จะทำให้ค่า ρa เลื่อนขึ้นหรือลงแบบทวีคูณ

การจัดการการเปลี่ยนเกียร์แบบคงที่สามารถทำได้โดยวิธีดังต่อไปนี้:
1. การแก้ไขโดยใช้ข้อมูล TDEM/CSAMT เพื่อควบคุมค่าความต้านทานไฟฟ้าในระดับตื้น
2. การผกผันร่วมของ MT-TDEM เพื่อทำให้แบบจำลองระดับตื้นมีความแม่นยำมากขึ้น
3. แนวทางการผกผันที่แข็งแกร่งสำหรับการปรับเปลี่ยน เช่น การอนุญาตให้กำหนดพารามิเตอร์การปรับเปลี่ยนต่อสถานี

นอกจากนี้ ยังมีการตรวจจับค่าผิดปกติของความถี่ด้วย โดยปกติแล้ว จุดบางจุดบนเส้นโค้งที่เบี่ยงเบนอย่างรวดเร็วจะถูกตัดทิ้งหรือกำหนดค่าความคลาดเคลื่อนที่มากขึ้น

-

ขั้นตอนที่ 5: การวิเคราะห์มิติและการโจมตี
ก่อนทำการวิเคราะห์แบบผกผัน จำเป็นต้องตรวจสอบว่าโครงสร้างใต้ผิวดินเป็นแบบ 1 มิติ 2 มิติ หรือ 3 มิติ ซึ่งจะมีผลต่อการเลือกวิธีการวิเคราะห์แบบผกผันและการตีความ การวิเคราะห์ที่ใช้กันทั่วไป ได้แก่:

– พารามิเตอร์ความเบี่ยงเบน (เช่น ความเบี่ยงเบนของ Bahr) เพื่อประเมินระดับ 3 มิติ
– เทนเซอร์เฟสเพื่อดูทิศทางหลักของโครงสร้างโดยไม่ได้รับผลกระทบจากการเลื่อนคงที่
– การวิเคราะห์แนวการวางตัวของหินเพื่อกำหนดทิศทางการวางตัวของหินที่เด่นชัดของโครงสร้าง 2 มิติ

หากข้อมูลแสดงลักษณะ 2 มิติที่ชัดเจน มักจะทำการหมุนเทนเซอร์อิมพีแดนซ์ตามทิศทางการกระทบเพื่อชี้แจงส่วนประกอบหลัก (TE/TM) หากลักษณะ 3 มิติชัดเจน จะนิยมใช้การผกผัน 3 มิติมากกว่า

-

ขั้นตอนที่ 6: การกลับข้อมูล MT (1 มิติ, 2 มิติ, 3 มิติ)
การตีความเชิงปริมาณของ MT โดยทั่วไปจะทำผ่านกระบวนการผกผัน ซึ่งเกี่ยวข้องกับการค้นหาแบบจำลองความต้านทานที่เหมาะสมที่สุดกับผลตอบสนองของข้อมูล การผกผัน MT เป็นแบบไม่เชิงเส้นและไม่สมบูรณ์ จึงจำเป็นต้องมีการปรับเสถียรภาพเพื่อป้องกันไม่ให้แบบจำลองกลายเป็น "ผิดปกติ" ประเภทของการผกผัน:

1. การวิเคราะห์แบบ 1 มิติ: เหมาะสำหรับพื้นที่ที่มีชั้นตะกอน (เช่น ตะกอนแนวนอน) รวดเร็ว แต่มีข้อจำกัด
2. การวิเคราะห์แบบ 2 มิติ: เหมาะสำหรับโครงสร้างที่มีลักษณะยาว เช่น รอยเลื่อน ร่องลึก หรือระบบความร้อนใต้พิภพแบบมีทิศทาง
3. การผกผันแบบ 3 มิติ: วิธีนี้สมจริงที่สุดสำหรับธรณีวิทยาที่ซับซ้อน แต่ต้องการข้อมูลจำนวนมาก การคำนวณขนาดใหญ่ และการควบคุมคุณภาพอย่างเข้มงวด

โดยทั่วไป ฟังก์ชันวัตถุประสงค์ผกผันจะรวมเอาความคลาดเคลื่อนของข้อมูลและความหยาบของแบบจำลองไว้ด้วย:

\[
\Phi = \Phi_d + \lambda \Phi_m
\]

โดยที่ λ คือพารามิเตอร์การปรับค่าให้เรียบ การเลือกค่า λ นั้นมีความสำคัญมาก: ถ้าค่า λ น้อยเกินไปจะทำให้แบบจำลองหยาบเกินไป (เกิดการโอเวอร์ฟิตติ้ง) ถ้าค่า λ มากเกินไปจะทำให้แบบจำลองเรียบเกินไปและพลาดรายละเอียดทางธรณีวิทยาที่สำคัญ

-

ขั้นตอนที่ 7: การประเมินผลลัพธ์การผกผันและความไว
หลังจากทำการกลับด้านแบบจำลองแล้ว จะต้องมีการประเมินแบบจำลองนั้น ไม่ใช่แค่การประเมินด้วยสายตาเท่านั้น การประเมินนี้รวมถึง:

– ค่าความคลาดเคลื่อน (RMS): ว่าตรงกับเป้าหมายหรือไม่ (เช่น RMS ~ 1–2 ขึ้นอยู่กับนิยามของความคลาดเคลื่อน)
– การเปรียบเทียบเส้นกราฟข้อมูลกับผลตอบสนองของแบบจำลอง ณ แต่ละสถานี
– การทดสอบความละเอียด: ตัวอย่างเช่น การทดสอบแบบตารางหมากรุก หรือแผนที่ความไวต่อการสัมผัส
– อิทธิพลของแบบจำลองก่อนหน้า: ลองใช้แบบจำลองเริ่มต้นหลายแบบเพื่อดูความเสถียรของผลลัพธ์

หากแบบจำลองเปลี่ยนแปลงอย่างมากเมื่อเปลี่ยนค่าพารามิเตอร์เล็กน้อย นั่นหมายความว่าการตีความควรทำอย่างระมัดระวังมากขึ้น และอาจต้องใช้ข้อมูลเพิ่มเติม

-

การตีความทางธรณีวิทยา: ความสัมพันธ์ระหว่างค่าความต้านทานกับระบบใต้ดิน
ค่าความต้านทานไฟฟ้าไม่ใช่ "ชนิดของหิน" โดยตรง แต่เป็นการตอบสนองทางกายภาพแบบผสมผสาน อย่างไรก็ตาม รูปแบบทั่วไปบางอย่างมักใช้เป็นแนวทางได้:

– ความต้านทานต่ำ (นำไฟฟ้า): ดินเหนียว, การเปลี่ยนแปลงทางความร้อนใต้ดิน (ชั้นดินเหนียวปกคลุม), เขตของเหลวเกลือ, กราไฟต์ หรือแร่ซัลไฟด์
– ค่าความต้านทานสูง (ต้านทาน): หินอัคนีขนาดใหญ่ หินแห้ง บริเวณที่มีซิลิกา หรือฐานหินผลึก

ตัวอย่างเช่น ในการสำรวจพลังงานความร้อนใต้พิภพ แบบจำลองแบบดั้งเดิมมักแสดงให้เห็นว่า:
1. ปิดทับด้วยดินเหนียวที่มีคุณสมบัติเป็นตัวนำไฟฟ้า
2. อ่างเก็บน้ำมีความต้านทานมากกว่าด้านล่าง
3. เขตการไหลขึ้นที่ควบคุมโดยโครงสร้าง (รอยเลื่อน)
4. แหล่งความร้อนที่บางครั้งปรากฏเป็นค่าความต้านทานหรือความผิดปกติที่ซับซ้อน ขึ้นอยู่กับลักษณะทางธรณีวิทยาและอุณหภูมิ

การตีความที่ดีนั้นจำเป็นต้องเชื่อมโยงกับข้อมูลอื่นๆ เสมอ เช่น ธรณีวิทยาพื้นผิว ปรากฏการณ์ทางความร้อน ธรณีเคมี แรงโน้มถ่วง แผ่นดินไหว หรือข้อมูลจากบ่อเจาะ

-

ปิด
เทคนิคการประมวลผลและการตีความข้อมูล MT นั้นเชื่อมโยงกัน ตั้งแต่การควบคุมคุณภาพอนุกรมเวลา (QC) การประมาณค่าสเปกตรัม การประมวลผลอ้างอิงที่แข็งแกร่ง/ระยะไกล การแก้ไขการเลื่อนคงที่ การวิเคราะห์มิติ ไปจนถึงการผกผัน 2 มิติ/3 มิติ และการประเมินความละเอียด ความสำเร็จของ MT ไม่ได้ขึ้นอยู่กับซอฟต์แวร์การผกผันเท่านั้น แต่ยังขึ้นอยู่กับคุณภาพของการวัด ความเข้าใจเกี่ยวกับสัญญาณรบกวน และการบูรณาการเข้ากับบริบททางธรณีวิทยาด้วย

ด้วยขั้นตอนการทำงานที่เป็นระบบระเบียบและการตีความข้อมูลหลายแหล่ง MT จึงเป็นเครื่องมือที่มีประสิทธิภาพสำหรับการทำแผนที่โครงสร้างความต้านทานใต้ผิวดินและช่วยในการตัดสินใจในการสำรวจและงานวิจัยทางธรณีวิทยาที่ซับซ้อน