Phương pháp thăm dò điện trở suất là một trong những phương pháp địa vật lý được sử dụng phổ biến nhất. Phương pháp này sử dụng các loại hệ điện cực khác nhau tùy vào mục đích và điều kiện thi công thực địa, mỗi loại hệ điện cực có ưu và nhược điểm riêng.
52 Journal of Mining and Earth Sciences Vol 61, Issue (2020) 52 - 60 Inversion of multiple data sets acquired by different array configuration of geoelectrical resistivity method Thong Duy Kieu * Faculty of Oil and Gas, Hanoi University of Mining and Geology, Vietnam ARTICLE INFO Article history: Received 15th Dec 2019 Accepted 6th Jan 2020 Available online 28th Feb 2020 Keywords: Inversion, Resistivity, Integration, Electrode array ABSTRACT The geoelectrical resistivity method is one of the most commonly used geophysical methods This method uses different electrodes configuration, electrode array, depending on the purpose and conditions of the field, each type of array has its advantages and disadvantages Due to the development of data acquisition technology, it is common for geoelectrical instruments enable to record data arising from different electrode arrays with negligible real-time construction However, current software’s only allows to process for each individual electrode array Inverted models of different electrode array can be integrated to build a common earth model However, due to the nature of the geophysical inversion is non-unique solutions, it means that there will be an infinite of models that can be suitable for a measurement in a certain noise level This leads to the same measurement data in an area with different electrode array may produce different geoelectrical models making the dificulty for integration process To solve this problem, we utilise the simultaneous joint inversion algorithm of data sets arising from multiple electrode arrays The test results on synthetic data show that this combination is better than the solution of each individual electrode array The best result is a combination of pole - dipole (PD), dipole - pole (DP) and dipole - dipole (DD) Copyright © 2020 Hanoi University of Mining and Geology All rights reserved _ *Corresponding author E-mail: kieuduythong@humg.edu.vn DOI: 10.46326/JMES.2020.61(1).06 Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất Tập 61, Kỳ (2020) 52 - 60 53 Giải ngược tổ hợp tài liệu hệ điện cực phương pháp điện trở suất 2D Kiều Duy Thông* Khoa Dầu khí, Trường Đại học Mỏ - Địa chất, Việt Nam THƠNG TIN BÀI BÁO Q trình: Nhận 15/12/2019 Chấp nhận 06/01/2020 Đăng online 28/02/2020 Từ khóa: Giải ngược, Điện trở suất, Tổ hợp, Hệ điện cực TÓM TẮT Phương pháp thăm dò điện trở suất phương pháp địa vật lý sử dụng phổ biến Phương pháp sử dụng loại hệ điện cực khác tùy vào mục đích điều kiện thi công thực địa, loại hệ điện cực có ưu nhược điểm riêng Do phát triển công nghệ đo ghi số liệu, máy thăm dò điện trở suất đa kênh ghi đồng thời số liệu nhiều hệ điện cực với thời gian thi công thực địa không đáng kể so với đo ghi hệ điện cực Nhưng phần mềm xử lý tài liệu cho phép thực cho hệ điện cực riêng biệt Nếu muốn tận dụng ưu điểm hệ điện cực, tổng hợp mơ hình giải ngược hệ điện cực Tuy nhiên, toán giải ngược địa vật lý đa nghiệm, nghĩa số liệu đo khu vực với hệ điện cực khác nhau, sau giải ngược cho mơ hình địa điện khác làm cho trình xử lý tổng hợp gặp khó khăn Để giải vấn đề này, chúng tơi sử dụng thuật tốn giải ngược tổ hợp đồng thời số liệu nhiều hệ điện cực Kết thử nghiệm mơ hình cho thấy giải ngược tổ hợp tốt giải ngược hệ điện cực riêng biệt Các kết tốt kết hợp hệ điện cực - lưỡng cực (PD), lưỡng cực - cực (DP) lưỡng cực - lưỡng cực (DD) © 2020 Trường Đại học Mỏ - Địa chất Tất quyền bảo đảm Mở đầu Địa vật lý (ĐVL) cơng cụ hữu ích phục vụ cho mục đích khác Điểm mạnh phương pháp ĐVL khảo sát với nhiều tỉ lệ, không cần tiếp cận trực tiếp đối tượng Tài liệu địa vật lý thường xử lý để xây dựng mô hình phân bố tính chất vật lý liên quan đến đối _ *Tác giả liên hệ E - mail: kieuduythong@humg edu DOI: 10.46326/JMES.2020.61(1).06 tượng mặt đất, trình gọi giải ngược tài liệu địa vật lý, báo để ngắn gọn gọi tắt giải ngược Vấn đề nan giải giải ngược tính đa trị tốn Có nghĩa có vơ số mơ hình phân bố tính chất vật lý mặt đất giải thích cho số liệu đo mặt đất khoảng nhiễu định Do vậy, trình giải ngược địa vật lý cần có thêm thơng tin để hạn chế tính đa trị, hay làm hẹp miền nghiệm toán Một phương pháp hữu hiệu sử dụng rộng rãi giới sử dụng tổ hợp tài liệu phương pháp ĐVL để giải ngược (Gallardo and Meju, 2011; Heincke et 54 Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 al., 2017; Lines et al., 1988; Moorkamp, 2017; Moorkamp et al., 2011; Paasche and Tronicke, 2007; Sun and Li, 2015; Thông, 2018; Vozoff and Jupp, 1975a) Hơn nữa, độ phân giải phương địa vật lý khác nhau, nên việc kết hợp phương pháp trình giải ngược tổ hợp làm tăng độ phân giải kết giải ngược (Vozoff and Jupp, 1975b) Tổng quan giải ngược tổ hợp giới thiệu (Thông, 2018) Trong vấn đề giải ngược hai điểm quan trọng cần phải giải Thứ tìm mối liên hệ tham số vật lý mơ hình, mơ hình tham số vật lý phương pháp thành viên liên hệ với mối liên kết cấu trúc (structural link) thạch học (petrohysical link) (Thông, 2018) Thứ hai cách thức tiến hành giải ngược tổ hợp số liệu, giới gồm hai cách thức chủ yếu (1) giải ngược đồng thời tài liệu địa vật lý (joint inversion); số liệu phương pháp ĐVL khác tiến hành trình (2) Giải ngược tài liệu địa vật lý (co-operative inversion) Giải ngược phương pháp nối tiếp nhau, kết trình giải ngược phương pháp trước đưa vào đầu vào cho trình giải ngược phương pháp sau Mỗi cách thức liên kết mơ hình q trình giải ngược có ưu điểm nhược điểm riêng áp dụng tùy theo tốn cụ thể (Thơng, 2018) Trong nghiên cứu áp dụng giải ngược tổ hợp tài liệu địa vật lý cho toán điện trở suất 2D Đây phương pháp sử dụng phổ biến với hệ điện cực thông dụng bao gồm: lưỡng cự - lưỡng cực (DD), cực - lưỡng cực (PD), lưỡng cực - cực (DP) Wenner Schlumberger (WS) (Hình 1) Mỗi loại hệ điện cực có ưu nhược điểm riêng (Athanasiou et al., 2007) Với phát triển công nghệ đo nghi, máy đo điện đa cực thường có khả đo liên tục hệ điện cực khác với chi phí thời gian không nhiều Do vậy, câu hỏi đặt liệu xử lý giải ngược tổ hợp hệ điện cực với nhằm tận dụng ưu điểm hệ điện cực có đem lại kết tốt kết xử lý truyền thống hệ điện cực hay không Do đặc điểm toán thực phương pháp điện trở suất nên áp dụng giải ngược tổ hợp đồng thời (joint inversion) Đã có nghiên cứu thực vấn đề (Athanasiou et al., 2007), nhiên nghiên cứu tác giả đánh giá dựa mơ hình đơn giản, mơ hình gồm khối đồng khơng tính đến ảnh hưởng bất đồng lớp Trong nghiên cứu sử dụng mơ hình gần với thực tế tính đến bất đồng gần mặt đất Tổ hợp hệ điện cực sử dụng khác so với nghiên cứu nghiên cứu Athanasiou et al., (2007) Giải ngược tài liệu địa vật lý 2.1 Lý thuyết giải ngược Giải ngược tài liệu ĐVL trình xác định phân bố tính chất vật lý khơng gian mặt đất từ tài liệu đo Thông thường giải ngược thực trình giải lặp, cực tiểu hóa hàm sau: ∅ = ∅𝒅𝒅 + 𝛽𝛽∅𝒎𝒎 (1) ∅𝑑𝑑 = ‖𝑾𝑾𝑾𝑾 − 𝑾𝑾𝑾𝑾(𝒎𝒎)‖𝟐𝟐𝟐𝟐 (2) Trong đó: ϕd - cơng thức (2), độ lệch (misfit) hay sai số (error) tài liệu đo tài liệu lý thuyết tính tốn từ mơ hình; ϕm - cơng thức (3), đánh giá tham số mơ hình, độ trơn (smoothness) nghiên cứu Constable et al (1987), bao gồm thêm kỹ thuật hạn chế nghiệm khác công trình Menke (2015); Sun and Li (2014, 2015); Tarantola and Valette (1982); β - (Lagarian regulerization parameter) tham số cân sai số tham số cấu trúc mơ hình (Tikhonov and Arsenin, 1977) 𝟐𝟐 �𝒎𝒎� ∅𝒎𝒎 = �𝝏𝝏 𝟐𝟐 (3) Trong ‖ ‖𝟐𝟐 - tiêu chuẩn Euclidian, d=(d1, d2, …, dM)T m=(m1, m2, …, mN)T số liệu tham số mơ hình, f(m) - tốn tử mơ hình hóa, địa vật lý toán tử thường hàm phi tuyến; W ma trận đường chéo NxN phương sai số liệu � toán tử hiệu (Different operator) σ(4), 𝝏𝝏 𝑾𝑾 = 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑{1⁄𝜎𝜎1 , 1⁄𝜎𝜎2 , … , 1⁄𝜎𝜎𝑁𝑁 }, (4) Để giải toán giải ngược với hàm phi tuyến địa vật lý, thuật toán thường sử dụng giải lặp (Constable et al., 1987) Hàm f(m) sấp xỉ thành chuỗi Taylor sau: Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 𝒇𝒇�𝒎𝒎(𝟏𝟏) + ∆𝒎𝒎� = 𝒇𝒇(𝒎𝒎(𝟏𝟏) ) + 𝑱𝑱(𝟏𝟏) ∆𝒎𝒎 (5) Trong ∆m=m(2)-m(1) véc tơ thể thay đổi tham số mơ hình, gọi tắt véc tơ biến đổi mơ hình, J(1) ma trận đạo hàm bán phần hàm f(m) theo vec tơ tham số mơ hình m(1), Jacobian Từ phương trình (1), (2) (3), (4) (5) ta có: 𝟐𝟐 �𝒎𝒎(𝟐𝟐) �𝟐𝟐 (6) ∅ = �𝑾𝑾𝒅𝒅(𝟏𝟏) − 𝑾𝑾𝑱𝑱(𝟏𝟏) 𝒎𝒎(𝟐𝟐) � + 𝛽𝛽�𝝏𝝏 Trong 𝟐𝟐 𝒅𝒅(𝟏𝟏) = 𝒅𝒅 − 𝒇𝒇�𝒎𝒎(𝟏𝟏) � + 𝑱𝑱(𝟏𝟏) 𝒎𝒎(𝟏𝟏) 𝟐𝟐 (7) Đạo hàm hàm số (6) theo tham số mơ hình m(2): 𝝏𝝏∅ = 2�(𝐉𝐉(𝟏𝟏) 𝐓𝐓 𝐀𝐀𝐉𝐉(𝟏𝟏) + 𝛽𝛽𝐇𝐇�𝐦𝐦(𝟐𝟐) 𝝏𝝏𝐦𝐦(𝟐𝟐) − 2𝐉𝐉(𝟏𝟏) 𝐓𝐓 𝐀𝐀𝐝𝐝(𝟏𝟏) (8) �𝑻𝑻 𝝏𝝏 � A=WTW, cho đạo Trong đó: 𝑯𝑯 = 𝝏𝝏 hàm khơng, ta tìm véc tơ tham số mơ hình thỏa mãn điều kiện phiếm hàm (6) cực tiểu hóa Sau vài phép biến đổi ta tìm véc tơ biến đổi mơ hình sau bước lặp thứ i (9) (9) −𝟏𝟏 ∆𝒎𝒎(𝒊𝒊+𝟏𝟏) = �(𝑱𝑱(𝒊𝒊) 𝑻𝑻 𝑨𝑨𝑱𝑱(𝒊𝒊) + 𝛽𝛽𝑯𝑯� �𝑱𝑱(𝒊𝒊) 𝑻𝑻 𝑨𝑨∆𝒅𝒅(𝒊𝒊) � Trong ∆𝒅𝒅(𝒊𝒊) = 𝒅𝒅 − 𝒇𝒇�𝒎𝒎(𝒊𝒊) � sai số số liệu thực tế đo số liệu lý thuyết tình từ mơ hình với tham số m(i) 2.2 Giải ngược tổ hợp tài liệu điện trở suất từ hệ điện cực Trong phương pháp điện trở suất sử dụng nhiều hệ điện cực khác (Szalai and Szarka, 2008) Mỗi loại hệ điện cực có ưu điểm nhược Wenner - Schlumberger (WS) C1 P1 P2 C2 Dipole - Dipole (DD) C1 C2 P1 P2 55 điểm riêng, tùy thuộc vào mơ hình địa điện (Athanasiou et al., 2007; Dahlin and Zhou, 2004, 2006; Szalai et al., 2011) Trên thực tế phương pháp điện trở suất thường sử dụng hệ điện cực mơ tả Hình Trong trường hợp tổng quát, hệ điện cực WS thường có độ phân giải đứng tốt, ngược lại độ phân giải ngang hệ điện cực DD, PD DP tốt (Dahlin and Zhou, 2004) Chiều sâu nghiên cứu hệ điện cực PD DP tốt hệ điện cực DD hệ điện cực WS (Szalai et al., 2011) Tuy nhiên ảnh hưởng bất đồng địa phương gần mặt đất, thường không liên quan đến đối tường nghiên cứu lại tăng dần từ hệ điện cực WS, đến PD DP cao hệ điện cực DD (Dahlin and Zhou, 2004) Do vậy, giải ngược tổ hợp tài liệu hệ điện cực cho kết tốt giải ngược tài liệu hệ điện cực riêng lẻ, tận dụng ưu điểm hệ điện cực để bù đắp nhược điểm hệ điện cực Do điện trở suất biểu kiến mơ hình địa điện phụ thuộc vào hệ điện cực sử dụng Để tránh việc giá trị điện trở suất biểu kiến hệ điện cực chiếm vai trò định q trình giải ngược tổ hợp Ví dụ trường hợp hệ điện cực DD thường có ảnh hưởng định đến q trình giải ngược tổ hợp (Athanasiou et al., 2007), nên phải thiết lập hệ số cân Trong cơng trình Athanasiou et al (2007), tác giả dựa vào độ nhạy mơ hình, hay ma trận Jacobian để tính hệ số Kết cơng trình cho thấy, khó để đánh giá hiệu sử dụng hệ số cân Chúng tơi cho ngun nhân việc tính hệ số theo độ nhạy mơ hình khơng phù hợp, chất hệ số cân phải liên quan đến điện trở suất biểu kiến Do vậy, nghiên cứu chúng tơi sử dụng cách tính hệ số cân dựa vào điện trở C1 Pole - Dipole (PD) P1 P2 Dipole - Pole (DP) C1 C2 P2 Hình 1: Các hệ điện cực thường dùng: Wenner - Schlumberger (WS), lưỡng cực trục (DD), ba cực thuận (PD), ba cực nghịch (DP) Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 56 suất biểu kiến, kỹ thuật đề xuất Vozoff and Jupp (1975b) Kết thảo luận 3.1 Mơ hình (a) Khoảng cách tuyến (m) Res1 (e) Chiều sâu (m) Chiều sâu (m) (d) Khoảng cách tuyến (m) Res2 Res3 Bảng Giá trị điện trở suất mơ hình (Hình 2) Trong rand hàm ngẫu nhiên Giá trị điện trở suất nhóm ký hiệu Res1, Res2, Res3, Res4 Nhóm Cơng thức tính Giá trị điện trở suất (Ωm) 1.01 𝑅𝑅𝑅𝑅𝑠𝑠1 1.07 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 � −0.5� 1.29 = 10 1.08 45.21 𝑅𝑅𝑅𝑅𝑠𝑠2 42.15 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 � +1� 46.24 = 10 31.25 314.43 𝑅𝑅𝑅𝑅𝑠𝑠3 339.80 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 � +2� 473.00 = 10 326.74 1622.43 𝑅𝑅𝑅𝑅𝑠𝑠4 1371.76 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 � +3.1� 1941.15 = 10 1973.72 Chiều sâu (m) Chiều sâu (m) Trong nghiên cứu chúng tơi thiết lập mơ hình lý thuyết (Hình 2) Giá trị điện trở suất mơ hình chia làm nhóm tính Bảng Có thể thấy nhóm điện trở bao gồm nhóm điện trở suất thấp (Res1), nhóm điện trở suất trung bình (Res2), nhóm điện trở suất cao (Res3) nhóm điện trở suất cao (Res4) Trong mơ hình (Hình 1), giá trị khối mơ hình xếp ngẫu nhiên Lớp gồm giá trị điện trở suất hai nhóm Res2 Res3 Lớp đặc trưng cho lớp phủ, thường có có giá điện trở suất thay đổi mạnh Lớp thứ hai có giá trị điện trở suất nhóm Res2, lớp cuối có giá trị điện trở suất nhóm Res3 Hai khối bất đồng địa phương O1 O2 có giá trị điện trở suất thấp Res1 cao Res4 Số liệu tính từ mơ hình cộng thêm phần nhiễu ngẫu nhiên, 5% giá trị điện trở suất biểu kiến, cho số liệu mô gần với thực tế Res4 O1 (b) Khoảng cách tuyến (m) (c) Khoảng cách tuyến (m) Điện trở suất (Ωm) O2 L1 L2 L3 H1 H2 Hình Mơ hình địa điện (e) bao gồm ba lớp L1, L2 L3 phân chia hai ranh giới H1 H2 (các đường màu đỏ), có hai khối bất đồng địa phương O1 O2.Giá trị điện trở suất mơ hình thể Bảng Giá trị điện trở suất biểu kiến tính từ mơ hình cộng với 5% nhiễu ngẫu nhiên hệ điện cực DD (a), PD (b), DP (c) WS (d) thực phần mềm Res2dmod [23] Khoảng cách hai điện cực liên tiếp 4.0 m Trong hình (a), (b), (c) (d) dải màu thể giá trị logarit số 10 điện trở suất biểu kiến; chấm màu đen thể vị trí điểm lấy số liệu Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 Trong mơ hình việc sử dụng tài liệu đo điện trở suất để xác định ranh giới H1 H2 tương đối dễ ràng, nhiên việc xác định đối tượng O1 O2 khó khăn, O1 Thứ ảnh hưởng bất đồng lớp L1 Thứ hai, đối tượng O1 có điện trở suất khác biệt lớn so với môi trường vây quanh, nằm chiều sâu lớn nằm lớp L2 có điện trở suất thấp Do vậy, điểm quan trọng toán đánh giá xem liệu hệ điện cực phát đối tượng O1 O2 3.2 Kết giải ngược 3.2.1 Kết giải ngược cho hệ điện cực Để kiểm chứng hiệu việc kết hợp hệ điện cực Đầu tiên chạy giải ngược cho hệ điện cực, sử dụng phần mềm sử dụng rộng rãi giới, Res2dinv [23] Kết thể Hình Các kết cho thấy nhìn chung ranh giới H1 H2 phân chia lớp L1, L2 L3 thể mơ hình giải ngược bốn hệ điện cực Đối tượng O2 nhìn thấy mơ hình này, đối tượng O1 tất DD 57 trường hợp không rõ ràng Đối tượng O1 tạo thành đới dị thường điện trở suất thấp kết hệ điện cực DD, PD DP giống dị thường đới đứt gãy, chí O1 khơng thấy thể kết hệ điện cực WS 3.2.2 Kết giải ngược tổ hợp Các kết giải ngược thực chương trình Matlab Chúng tơi thay đổi chương trình từ chương trình gốc cơng bố Akca (2016) cho q trình giải ngược báo Tổ hợp số liệu hệ điện cực trình bày Bảng Trên Hình thể độ lệch (misfit) số liệu tính tốn lý thuyết tài liệu đo thực tế (ở tài liệu tính theo mơ hình cộng với nhiễu ngẫu nhiên) Giá trị độ lệch nhỏ thể số liệu tính lý thuyết số liệu đo gần Độ lệch suy giảm cách ổn định theo bước lặp thể chương trình giải ngược làm việc tốt Nhìn chung khác biệt độ lệch tổ hợp số liệu điện cực không nhiều, đề nằm mức nhiễu 5% Cần lưu ý độ lệch nhỏ không đảm bảo kết giải ngược tốt hơn, mà cần độ lệch tương đương mức nhiễu (Constable et al., 2015) DP Chiều sâu (m) Chiều sâu (m) Khoảng cách tuyến (m) Chiều sâu (m) Khoảng cách tuyến (m) PD Chiều sâu (m) WS Khoảng cách tuyến (m) Khoảng cách tuyến (m) Hình Kết giải ngược tài liệu hệ điện cực khác nhau, sử dụng phần mềm Res2dinv [23] Các đường màu trắng thể vị trí ranh giới H1 H2 hai đối tượng địa phương O1 O2 Bảng Ký hiệu tổ hợp số liệu từ hệ điện cực khác Lưu ý viết tắt hệ điện cực xem Hình TT Ký hiệu Tổ hợp số liệu hệ điện cực TT Ký hiệu Tổ hợp số liệu hệ điện cực PD-DD PD DD PD-DP-DD PD; DP; DD PD-DP PD DP PD-DP-WS PD; DP WS PD-WS PD WS PD-DD-WS PD; DD; WS WS-DD WS DD PD-DP-DD-WS PD; DP; DD; WS Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 58 Các mơ hình giải ngược tổ hợp số liệu thể Hình Có thể rút số nhận xét sau kết Thứ nhất, so với kết giải ngược tài liệu hệ điện cực riêng rẽ kết tổ hợp nhìn chung tốt Lưu ý thang màu hai kết giải ngược Hình đặt giống Ngoài hai ranh giới H1 H2 đối tượng dễ xác định, đối tượng địa phương O1 O2 xác định trình giải ngược tổ hợp điện cực tốt so với sử dụng riêng hệ điện cực Thứ 2, tổ hợp có sử dụng hệ điện cực WS cho kết không tốt hệ điện cực lại Tổ hợp hệ điện cực PD DD PD DP cho kết tốt, kết giải ngược tốt với tổ hợp PD, DP DD Điều cho thấy kết Athanasiou et al (2007) thực tổ hợp hai hệ điện cực WS với DD WS với PD chưa thực hiệu Đây khác biệt đóng góp quan trọng nghiên cứu so với nghiên cứu cơng bố Hình Độ lệch (misfit) tài liệu đo tài liệu tính lý thuyết thay đổi theo bước lặp (Iterations) PD-DD PD-WS Chiều sâu (m) Chiều sâu (m) Khoảng cách tuyến (m) Khoảng cách tuyến (m) Chiều sâu (m) Chiều sâu (m) PD-DP Khoảng cách tuyến (m) Khoảng cách tuyến (m) Chiều sâu (m) Chiều sâu (m) PD-DP-DD Khoảng cách tuyến (m) Khoảng cách tuyến (m) PD-DP-WS Khoảng cách tuyến (m) Chiều sâu (m) Chiều sâu (m) PD-DP-DD-WS DD-WS PD-DD-WS Khoảng cách tuyến (m) Hình Kết giải ngược tổ hợp tài liệu hệ điện cực khác (Bảng 2), sử dụng chương trình Matlab Kết tổ hợp PD-DP-DD PD-DP-DD-WS tốt Các đường màu trắng thể vị trí ranh giới H1 H2 hai đối tượng địa phương O1 O2 Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 Kết luận Kết từ mơ hình lý thuyết cho thấy giải ngược tổ hợp số liệu hệ điện cực cho mơ hình tốt so với mơ hình hệ điện cực riêng lẻ Sử dụng tổ hợp số liệu hệ điện cực đòi hỏi chi phí thực địa tính tốn cao Tuy nhiên không thực vấn đề lớn nhờ phát triển công nghệ đo đạc sức mạnh máy tính Kết tổ hợp hệ điện cực: cực - lưỡng cực (PD), lưỡng cực - cực (DP) lưỡng cực - lưỡng cực (DD) với mô hình dạng cho kết tốt Trong nghiên cứu tiến hành thử nghiệm tài liệu lý thuyết, mô gần với điều kiện thực tế Trong nghiên cứu thử nghiệm tài liệu đo thực tế để đánh giá xác khả áp dụng vào thực tế sản xuất phương pháp Lời cảm ơn Bài báo hoàn thành kết đề tài cấp sở Trường Đại học Mỏ - Địa chất, mã số T18-36 Bản quyền phần mềm Matlab tài trợ công ty RoqSense Chúng chân thành cảm ơn Akca tác giả báo “ELRIS2D: A MATLAB Package for the 2D Inversion of DC Resistivity/IP Data” cung cấp mã nguồn cho chương trình giải ngược tài liệu điện trở hai chiều Chúng trân trọng cảm ơn tiến sỹ Loke, người phát triển phần mềm Res2dmod Res2dinv cung cấp dùng miễn phí với mục đích nghiên cứu Tài liệu tham khảo Akca, I., (2016) ELRIS2D: A MATLAB Package for the 2D Inversion of DC Resistivity/IP Data, Acta Geophysica 64 443 Athanasiou, E N., Tsourlos, P I., Papazachos, C B., and Tsokas, G N., (2007) Combined weighted inversion of electrical resistivity data arising from different array types: Journal of Applied Geophysics 62(2) 124 - 140 Constable, S C., Parker, R L., and Constable, C G., (1987) Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data: Geophysics 52(3) 289 - 300 59 Constable, S., Orange, A., and Key, K., (2015) And the geophysicist replied: “Which model you want?”: Geophysics 80(3) E(197-E212) Dahlin, T., and Zhou, B., (2004) A numerical comparison of 2D resistivity imaging with 10 electrode arrays: Geophysical Prospecting, v 52, no 5, p 379-398.-, (2006), Multiplegradient array measurements for multichannel 2D resistivity imaging: Near Surface Geophysics 4(2) 113 - 123 Gallardo, L A., and Meju, M A., (2011) Structurecoupled multiphysics imaging in geophysical sciences: Reviews of Geophysics 49(1) Heincke, B., Jegen, M., Moorkamp, M., Hobbs, R W., and Chen, J., (2017) An adaptive coupling strategy for joint inversions that use petrophysical information as constraints: Journal of Applied Geophysics 136 279 - 297 Lines, L., Schultz, A., and Treitel, S., (1988) Cooperative inversion of geophysical data: Geophysics 53(1) - 20 Menke, W., (2015) Review of the Generalized Least Squares Method: Surveys in Geophysics 36(1) - 25 Moorkamp, M., (2017) Integrating Electromagnetic Data with Other Geophysical Observations for Enhanced Imaging of the Earth: A Tutorial and Review: Surveys in Geophysics Moorkamp, M., Heincke, B., Jegen, M., Roberts, A W., and Hobbs, R W., (2011) A framework for 3-D joint inversion of MT, gravity and seismic refraction data: Geophysical Journal International 184(1) 477 - 493 Paasche, H., and Tronicke, J., (2007) Cooperative inversion of 2D geophysical data sets: a zonal approach based on fuzzy c-means cluster analysis: Geophysics 72(3) 35 - 39 Sun, J., and Li, Y., (2014) Adaptive Lp inversion for simultaneous recovery of both blocky and smooth features in a geophysical model: Geophysical Journal International 197(2) 882 899; (2015), Multidomain petrophysically constrained inversion and geology differentiation using guided fuzzy c-means clustering: Geophysics 80(4) ID1-ID18 60 Kiều Duy Thơng/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 61 (1), 52 - 60 Szalai, S., and Szarka, L., (2008) On the classification of surface geoelectric arrays: Geophysical Prospecting 56(2) 159 - 175 Szalai, S., Novák, A., and Szarka, L., (2011) Which geoelectric array sees the deepest in a noisy environment? Depth of detectability values of multielectrode systems for various twodimensional models: Physics and Chemistry of the Earth, Parts A/B/C 36(16) 1398 - 1404 Tarantola, A., and Valette, B., (1982) Generalized nonlinear inverse problems solved using the least squares criterion: Reviews of Geophysics 20(2) 19 - 232 Tikhonov, A N., and Arsenin, V I A., (1977) Solutions of ill-posed problems, Winston Thông, K D., (2018) Tổng quan giải ngược tổ hợp tài liệu địa vật lý, Hội nghị toàn quốc khoa học trái đất phát triển bền vững (ERSD) Trường, Đại học Mỏ - Địa chất Vozoff, K., and Jupp, D L B., (1975a) Joint Inversion of Geophysical Data: Geophysical Journal of the Royal Astronomical Society, v 42, no 3, p 977-991.-, (1975b, Joint Inversion of Geophysical Data 42(3) 977 - 991 ... 28/02/2020 Từ khóa: Giải ngược, Điện trở suất, Tổ hợp, Hệ điện cực TÓM TẮT Phương pháp thăm dò điện trở suất phương pháp địa vật lý sử dụng phổ biến Phương pháp sử dụng loại hệ điện cực khác tùy vào... số liệu nhiều hệ điện cực Kết thử nghiệm mơ hình cho thấy giải ngược tổ hợp tốt giải ngược hệ điện cực riêng biệt Các kết tốt kết hợp hệ điện cực - lưỡng cực (PD), lưỡng cực - cực (DP) lưỡng cực. .. hệ điện cực cho kết tốt giải ngược tài liệu hệ điện cực riêng lẻ, tận dụng ưu điểm hệ điện cực để bù đắp nhược điểm hệ điện cực Do điện trở suất biểu kiến mơ hình địa điện phụ thuộc vào hệ điện