Luận chứng hiệu quả áp dụng phương pháp phân cực kích thích trong tìm kiếm quặng đồng , vàng, chì , kẽm ở việt nam

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT - TRẦN BÌNH TRỌNG LUẬN CHỨNG HIỆU QUẢ ÁP DỤNG PHƯƠNG PHÁP PHÂN CỰC KÍCH THÍCH TRONG TÌM KIẾM QUẶNG ĐỒNG - VÀNG, CHÌ - KẼM Ở VIỆT NAM LUẬN ÁN TIẾN SĨ ĐỊA CHẤT Hà Nội - 2008 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT - TRẦN BÌNH TRỌNG LUẬN CHỨNG HIỆU QUẢ ÁP DỤNG PHƯƠNG PHÁP PHÂN CỰC KÍCH THÍCH TRONG TÌM KIẾM QUẶNG ĐỒNG VÀNG, CHÌ - KẼM Ở VIỆT NAM Chuyên ngành: Các phương pháp Địa vật lý Địa hóa học tìm kiếm thăm dị mỏ khống sản Mã số: 1.06.17 LUẬN ÁN TIẾN SĨ ĐỊA CHẤT NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Nguyễn Trọng Nga TS Nguyễn Ngọc Loan Hà Nội - 2008 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết luận án trung thực chưa công bố công trình Nghiên cứu sinh Trần Bình Trọng MỤC LỤC DANH MỤC CÁC BẢNG .4 DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ .5 MỞ ĐẦU Chương 1- ĐẶC ĐIỂM, ĐỊA VẬT LÝ CÁC LOẠI HÌNH QUẶNG ĐỒNG VÀNG, CHÌ KẼM Ở VIỆT NAM .16 1.1 Đặc điểm địa chất loại hình quặng đồng - vàng, chì - kẽm Việt Nam 16 1.2 Đặc điểm tham số địa vật lý loại quặng đồng - vàng, chì - kẽm Việt Nam 35 Chương 2- NGHIÊN CỨU QUÁ TRÌNH SUY GIẢM HỆ SỐ PHÂN CỰC KÍCH THÍCH THEO THỜI GIAN 41 2.1 Các quy luật suy giảm độ phân cực thực nghiệm 42 2.2 Tham số suy giảm độ phân cực số khoáng vật quặng kim loại phổ biến 45 2.3 Tham số suy giảm độ phân cực số khoáng vật quặng kim loại theo kết đo mẫu tham số phòng .48 2.4 Tham số suy giảm độ phân cực theo kết đo tham số vết lộ quặng đồng (vàng) Tả Phời, Lào Cai 52 2.5 Tham số suy giảm độ phân cực theo kết đo tham số vết lộ quặng chì kẽm, vùng Ba Xứ ,Tuyên Quang 60 2.6 Minh chứng hiệu tham số suy giảm mơ hình 2D .63 Chương 3- NGHIÊN CỨU ÁP DỤNG PHƯƠNG PHÁP PHÂN CỰC KÍCH THÍCH 3D NHẰM NÂNG CAO HIỆU QUẢ TRONG ĐÁNH GIÁ QUẶNG ĐỒNG VÀNG - CHÌ KẼM Ở VIỆT NAM 67 3.1 Cơ sở lý thuyết toán thuận ba chiều .67 3.2 Kỹ thuật đo phân cực kích thích 3D 77 3.3 Cơ sở lý thuyết phương pháp toán ngược chương trình ứng dụng xử lý số liệu đo điện 3D 86 3.4 Áp dụng kết tốn thuận phân cực kích thích 3D để tính tốn cho số mơ hình khống sản chì - kẽm đồng - vàng Việt Nam 99 3.5 Phương pháp xử lý tương quan cặp .125 3.6 Quy trình cơng nghệ đo 3D phân cực kích thích: 126 Chương 4- KẾT QUẢ ÁP DỤNG THỬ NGHIỆM PHƯƠNG PHÁP PHÂN CỰC KÍCH THÍCH 3D TRONG ĐÁNH GIÁ QUẶNG ĐỒNG VÀNG, CHÌ KẼM Ở VIỆT NAM 130 4.1 Cơ sở lựa chọn đối tượng diện tích đo thử nghiệm .130 4.2 Đặc điểm địa chất 131 4.3 Lựa chọn thiết bị đo thử nghiệm 139 4.4 Kết áp dụng thử nghiệm phương pháp phân cực kích thích 3D đánh giá quặng đồng (vàng) vùng Tả Phời, Lào Cai 144 4.5 Kết áp dụng thử nghiệm phương pháp phân cực kích thích 3D đánh giá quặng chì - kẽm vùng Ba Xứ, Tuyên Quang .156 KẾT LUẬN VÀ KIẾN NGHỊ 165 CÁC CƠNG TRÌNH KHOA HỌC CỦA TÁC GIẢ ĐÃ CƠNG BỐ CÓ LIÊN QUAN ĐẾN NỘI DUNG LUẬN ÁN 167 TÀI LIỆU THAM KHẢO 169 PHỤ LỤC 173 DANH MỤC CÁC BẢNG Bảng 1.1 Tổng hợp tham số điện phân cực kích thích loại quặng chì - kẽm đá vây quanh quặng Việt Nam [14] .37 Bảng 1.2 Tổng hợp tham số điện phân cực kích thích loại quặng đồng vàng đá vây quanh quặng Việt Nam [14] 38 Bảng 2.1 Tổng hợp kết đo độ phân cực kích thích theo thời gian số khoáng vật dẫn điện gây hiệu ứng phân cực phổ biến (Colect, 1959) [27] 46 Bảng 2.2 Bảng tính tốn tham số đặc trưng cho quy luật suy giảm số khoáng vật đặc trưng 47 Bảng 2.3 Bảng tính tốn tham số đặc trưng cho quy luật suy giảm số khoáng vật quặng sulfur đa kim Việt Nam 49 Bảng 2.4 Bảng tính tốn tham số đặc trưng cho quy luật suy giảm đá vây quanh quặng Việt Nam 49 Bảng 2.5 Bảng tính tổng hợp hệ số tương quan số khoáng vật quặng sulfur đa kim đá vây quanh Việt Nam .50 Bảng 2.6 Tổng hợp tính tham số suy giảm theo quy luật hàm loga, hàm e mũ độ phân cực vết lộ đá quặng đồng vùng Tả Phời, Lào Cai 56 Bảng 2.7 Tổng hợp tính tham số suy giảm theo quy luật hàm mũ, hàm hyperbol độ phân cực vết lộ đá quặng đồng vùng Tả Phời, Lào Cai 56 Bảng 2.8 Tổng hợp hệ số tương quan quy luật suy giảm vết lộ đá quặng đồng vùng Tả Phời, Lào Cai .59 Bảng 3.1 Tổng hợp kết tính tốn độ phân cực kích thích mơ hình chì - kẽm đá phiến vôi theo thiết bị đo hướng 104 Bảng 3.2 Tổng hợp kết tính tốn điện trở suất mơ hình chì - kẽm đá phiến vơi theo thiết bị đo hướng 105 Bảng 3.3 Tổng hợp nhận định kết tính tốn mơ hình chì - kẽm đá phiến vôi thiết bị 3D phân cực kích thích với mạng lưới .109 Bảng 3.4 Tổng hợp nhận định kết tính tốn mơ hình đồng đá phiến, điorit 119 Bảng 3.5 Tổng hợp nhận định kết tính tốn mơ hình mơ hình vàng thạch anh sulfur đá bazan 123 Bảng 4.1 Bảng thống kê tham số suy giảm độ phân cực 145 Bảng 4.2 Hệ số tương quan hàm lượng đồng tham số suy giảm độ phân cực kích thích 145 Bảng 4.3 Tổng hợp sai số dự báo hàm lượng Cu theo kết đo KCKT 3D 147 Bảng 4.4 Thống kê tham số phân cực quặng chì - kẽm Ba Xứ - Tuyên Quang 157 Bảng 4.5 Hệ số tương quan 157 DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1.1 Bản đồ phân bố khống sản chì - kẽm miền Bắc Việt Nam 20 Hình 1.2 Bản đồ phân bố khống sản đồng miền Bắc Việt Nam 21 Hình 1.3 Bản đồ phân bố khống sản vàng miền Bắc Việt Nam 22 Hình 1.4 Phân bố thân quặng đồng nguồn gốc nhiệt dịch (mỏ Sinh Quyền) [2]31 Hình 1.5 Hình thái thân quặng mỏ quặng chì - kẽm Việt Nam 32 Hình 1.6 Hình thái thân quặng mỏ chì - kẽm nhiệt dịch Kazaxtan [13] .33 Hình 1.7 Hình thái thân quặng mỏ vàng Thèn Sin, Phong Thổ, Lai Châu [3] 34 Hình 2.1 Đường cong đo độ phân cực suy giảm đường cong tính theo phương trình thực nghiệm 43 Hình 2.2 Biểu đồ biểu diễn quy luật suy giảm độ phân cực số khoáng vật quăng kim loại đặc trưng 47 Hình 2.3 Biểu đồ so sánh hệ số tương quan số liệu thực tế tính toán quy luật suy giảm độ phân cực số khoáng vật đặc trưng 48 Hình 2.4 Biểu đồ so sánh hệ số tương quan số liệu thực tế tính tốn lý thuyết hàm số suy giảm độ phân cực kết đo mẫu tham số miền Bắc Việt Nam 51 Hình 2.5 Tham số suy giảm độ phân cực theo thời gian đo vết lộ quặng đá vây quanh mỏ đồng Tả Phời, Lào Cai (Số liệu Nguyễn Văn Sâm -Liên đoàn InterGeo) 54 Hình 2.6 Biểu đồ nhận dạng quặng đồng đá vây quanh mỏ đồng Tả Phời, Lào Cai theo tham số Be điện trở suất đo vết lộ 54 Hình 2.7 Biểu đồ nhận dạng quặng đồng đá vây quanh mỏ đồng Tả Phời, Lào Cai theo tham số Be Ae tính từ số liệu PCKT đo vết lộ 55 Hình 2.8 Biểu đồ nhận dạng quặng đồng đá vây quanh mỏ đồng Tả Phời, Lào Cai theo tham số Bl Al tính từ số liệu PCKT đo vết lộ 55 Hình 2.9 So sánh độ phân cực (Eta), thời gian suy giảm 1/2 biên độ (St) điện trở suất theo loại quặng đá mỏ đồng (vàng) Tả Phời, Lào Cai 57 Hình 2.10 So sánh tham số suy giảm độ phân cực quy luật loga (Al, Bl) loại quặng đá mỏ đồng (vàng) Tả Phời, Lào Cai 57 Hình 2.11 So sánh tham số suy giảm độ phân cực quy luật hàm mũ số e (Ae, Be) loại quặng đá mỏ đồng (vàng) Tả Phời, Lào Cai 57 Hình 2.12 So sánh tham số Tham số suy giảm độ phân cực quy luật hàm mũ (Ax, Bx) loại quặng đá mỏ đồng (vàng) Tả Phời, Lào Cai 58 Hình 2.13 So sánh tham số suy giảm độ phân cực quy luật hàm hyperbol (Ah, Bh) loại quặng đá mỏ đồng (vàng) Tả Phời, Lào Cai 58 Hình 2.13 So sánh hệ số tương quan quy luật suy giảm theo loại quặng đá mỏ đồng (vàng) Tả Phời, Lào Cai 59 Hình 2.14 So sánh độ phân cực (η1) điện trở suất loại quặng đá đo mỏ chì – kẽm vùng Ba Xứ, Tuyên Quang 61 Hình 2.15 So sánh tham số suy giảm độ phân cực theo quy luật loga (Al, Bl) loại quặng đá đo mỏ chì – kẽm vùng Ba Xứ, Tuyên Quang 61 Hình 2.16 So sánh tham số suy giảm độ phân cực theo quy luật hàm mũ số e (Ae, Be) loại quặng đá đo mỏ chì – kẽm vùng Ba Xứ, Tuyên Quang.61 Hình 2.17 So sánh tham số suy giảm độ phân cực theo quy luật hàm mũ (Ax, Bx) loại quặng đá đo mỏ chì – kẽm vùng Ba Xứ, Tuyên Quang 62 Hình 2.18 So sánh tham số suy giảm độ phân cực theo quy luật hàm hyperbol (Ah, Bh) loại quặng đá đo mỏ chì – kẽm vùng Ba Xứ, Tuyên Quang 62 Hình 2.19 So sánh tham số suy giảm ½ biêm độ (St) tổng tích phân (Cha) loại quặng đá đo mỏ chì – kẽm vùng Ba Xứ, Tuyên Quang 62 Hình 2.19 Mặt cắt địa vật lý tuyến theo kết đo sâu phân cực kích thích .64 Hình 2.20 Mặt cắt xử lý định lượng Al, Bl theo quy luật logarit 65 Hình 2.21 Mặt cắt xử lý định lượng Ae, Be theo quy luật hàm e mũ 66 Hình 3.1 Rời rạc mơi trường thành khối phần tử hình chữ nhật 72 Hình 3.2 Vị trí điện cực tương ứng với điện cực dịng việc bố trí thiết bị sử dụng [28] 79 Hình 3.3 Sử dụng phương pháp đo chiếu để khảo sát lưới 10×10 với hệ thống máy đo điện trở suất 50 điện cực [28] .80 Hình 3.4 - Bố trí điện cực phép đo cực số đo (a) - theo chiều thuận, (b) với khoảng cách lớn điện cực (c) theo chiều ngược .81 Hình 3.5 Sơ đồ dịch chuyển hệ thiết bị cực - lưỡng cực đo đạc 3D 82 Hình 3.6 Sơ đồ bố trí hệ thiết bị lưỡng cực nghiên cứu 3D 85 Hình 3.7 Các mơ hình sử dụng giải định lượng 98 Hình 3.8 Mơ hình khơng gian 3D (điện trở suất độ phân cực kích thích) loại hình khống sản chì - kẽm đá phiến vơi 101 Hình 3.9 Kết tính mơ hình khơng gian 3D tham số điện trở suất loại hình chì kẽm dạng ổ, thấu kính đá phiến vôi .102 Hình 3.10 Kết tính mơ hình khơng gian 3D tham số độ phân cực kích thích loại hình chì - kẽm dạng ổ, thấu kính đá phiến vơi 103 Hình 3.11 Kết xử lý định lượng tham số η từ số liệu tính mơ hình lý thuyết quặng chì - kẽm dạng ổ, thấu kính theo hệ thiết bị đo khác 106 Hình 3.12 Mơ hình biểu diễn chiều theo kết xử lý định lượng theo số liệu tốn mơ hình lý thuyết hệ thiết bị cực - lưỡng cực 107 Hình 3.13 Mơ hình khơng gian 3D (điện trở suất độ phân cực kích thích) loại hình khống sản chì - kẽm đứt gãy 111 Hình 3.13a Tổng hợp kết tính tốn độ phân cực kích thích mơ hình chì - kẽm đứt gãy 112 Hình 3.14 Mơ hình khơng gian 3D (điện trở suất độ phân cực kích thích) loại hình khống sản chì - kẽm giả tầng đá vôi 114 Hình 3.14a Tổng hợp kết tính tốn độ phân cực kích thích mơ hình chì - kẽm giả tầng đá vôi 115 Hình 3.15 Mơ hình khơng gian 3D (điện trở suất độ phân cực kích thích) loại hình khống sản đồng (vàng) đá phiến, diorit .117 Hình 3.15a Tổng hợp kết tính tốn độ phân cực kích thích mơ hình khống sản đồng (vàng) đá phiến, diorit 118 Hình 3.16 Mơ hình khơng gian 3D (điện trở suất độ phân cực kích thích) loại hình khống sản vàng - thạch anh - sulfur đới dập vỡ 121 Hình 3.16a Tổng hợp kết tính tốn độ phân cực kích thích mơ hình khoáng sản vàng - thạch anh - sulfur đới dập vỡ 122 Hình 3.17 Sơ đồ quy trình cơng nghệ đo phân cực kích thích 3D 128 Hình 4.1 Bản đồ địa chất khoáng sản đồng vùng Tả phời [1] 133 Hình 4.2 Sơ đồ địa chất khoáng sản đồng vùng Tả phời thể chiều 133 Hình 4.3 Mặt cắt đặc trưng thân quặng đồng (vàng) vùng Tả Phời, Lào Cai [1] 134 Hình 4.4 Sơ đồ bố trí mạng lưới đo thử nghiệm 3D vùng Tả Phời 134 Hình 4.5 Sơ đồ địa chất khống sản tiểu khu Đỉnh Mười, vùng Ba Xứ [15] 137 Hình 4.6 Sơ đồ địa chất khoáng sản chiều diện tích đo thử nghiệm, vùng Đỉnh Mười 137 Hình 4.7 Sơ đồ bố trí mạng lưới đo thử nghiệm 3D vùng Đỉnh Mười 143 Hình 4.8 Mặt cắt 3D theo diện tích với độ sâu khác 148 Hình 4.9 Mặt cắt 3D theo mặt cắt theo trục OY 149 Hình 4.10 Biểu diễn kết mơ hình ba chiều (3D) 150 Hình 4.11 Mặt cắt 3D theo diện tích với độ sâu khác 151 Hình 4.12 Mặt cắt 3D theo trực OY tham số suy giảm quy luật hàm e mũ 152 Hình 4.13 Biểu diễn kết mơ hình ba chiều (3D) tham số Ae 153 Hình 4.14 Biểu diễn kết mơ hình ba chiều (3D) tham số Be 153 Hình 4.15 Thiết đồ lỗ khoan diện tích đo thử nghiệm [1] 154 Hình 4.16 Sơ đồ dự báo hàm lượng đồng không gian 3D theo kết đo điện phân cực kích thích 155 Hình 4.17 Mặt cắt 3D theo kết xử lý định lượng chương trình RES3DINV 160 Hình 4.18 Kết xử lý định lượng số liệu đo 3D vùng Đỉnh Mười, Ba Xứ 161 Hình 4.19 Kết xử lý toán 3D tham số suy giảm độ phân cực theo quy luật hyperbol tiểu khu Đỉnh Mười, Ba Xứ 162 Hình 4.20 Biểu diễn mơ hình ba chiều (3D) theo quy luật suy giảm hyperbol 163 Hình 4.21 Thiết đồ lỗ khoan 1, tiểu khu Đỉnh Mười, Ba Xứ, Tuyên Quang [15] 164 MỞ ĐẦU Tính cấp thiết đề tài Ở Việt Nam khoáng sản đồng - vàng, chì - kẽm nguyên liệu ngành luyện kim màu có vị trí vai trò quan trọng kinh tế quốc dân nhu cầu sử dụng với đà phát triển ngày lớn Tuy nhiên, mặt địa chất diện tích khơng lớn, lại chịu tác động nhiều chu kỳ kiến tạo nên thành tạo địa chất, mỏ quặng sulfur đa kim Việt Nam có thay đổi phức tạp Hiện nay, tổ hợp phương pháp điều tra, đánh giá thăm dị chì kẽm đồng - vàng thiếu vắng phương pháp địa vật lý, đặc biệt thăm dò điện Trong phương pháp thăm dị điện phương pháp có hiệu để phát hiện, theo dõi thân quặng đồng - vàng, chì - kẽm phương pháp phân cực kích Phương pháp phân cực kích thích khơng có nhiều ưu điểm phát loại quặng dẫn điện điện tử có dạng mạng mạch, mà quặng xâm tán, phương pháp điện trở, điện trường thiên nhiên hiệu phát Với trường hợp khống sản đồng - vàng, chì - kẽm liên quan với khống hố sulfur thường có dị thường phương pháp truyền thống độ phân cực biểu kiến (ηk) điện trở suất biểu kiến (ρk) có bề rộng lớn, gây khó khăn cho việc xác định vị trí cơng trình khai đào tỷ lệ cơng trình gặp quặng chưa cao, đặc biệt khả dự báo triển vọng khoáng sản theo dị thường đặc trưng PCKT chưa thực sản xuất Nhằm khắc phục tồn nêu phương pháp phân cực kích thích với chu ý phương pháp có chi phí lớn, nghiên cứu sinh nghiên cứu sử dụng tham số suy giảm độ phân cực kích thích (A), độ phân cực kích thích thời điểm ngắt dòng (B) để nâng cao khả phân biệt dị thường phân cực kích thích quặng đồng - vàng, chì - kẽm với đá vây quanh Với máy đo phân cực kích thích nay, đo nhiều thời gian nên xác định đường cong suy giảm độ phân cực kích thích theo thời gian 164 Hình 4.21 Thiết đồ lỗ khoan 1, tiểu khu Đỉnh Mười, Ba Xứ, Tuyên Quang [15] 165 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Từ kết nghiên cứu luận án đưa kết luận sau: Các tham số tốc độ suy giảm phân cực kích thích (A) có mối quan hệ với hàm lượng khoáng vật quặng chủ yếu thân quặng đồng - vàng, chì - kẽm, độ phân cực kích thích thời điểm ngắt dịng (B) có độ phân giải cao so với tham số ρk, ηk truyền thống Qua sử dụng chúng để phân tích, xử lý tài liệu phân cực kích thích 3D nâng cao hiệu phương pháp phân cực kích thích Qua nghiên cứu xây dựng mơ hình địa chất - địa vật lý tính tốn tốn thuận phương pháp phân cực kích thích 3D cho sáu hệ thiết bị cực - cực, cực lưỡng cực, lưỡng cực - lưỡng cực, lưỡng cực - lưỡng cực phương vị, Wenner Wenner – Schlumberger phù hợp với số mơ hình khống sản chì - kẽm đồng - vàng Việt Nam chứng minh khả áp dụng có hiệu phương pháp phân cực kích thích 3D Đã nghiên cứu, áp dụng thành cơng phương pháp phân cực kích thích 3D đánh giá, thăm dị khống sản đồng - vàng, chì - kẽm Việt Nam Đã xây dựng phương trình tính hàm lượng kim loại đồng theo tham số điện trở suất, tốc độ suy giảm độ phân cực kích thích thời điểm ngắt dịng Kết sử dụng để dự báo cho vùng mỏ đồng Tả Phời - Lào Cai vùng khác có điều kiện tương tự Xây dựng quy trình cơng nghệ đo, xử lý phân tích tài liệu phân cực kích thích 3D đánh giá thăm dị quặng đồng - vàng, chì - kẽm Quy trình có tính hợp lý, thực khả thi điều kiện trang thiết bị có Việt Nam Kiến nghị Phương pháp phân cực kích thích 3D phương pháp có hiệu giới, song chưa áp dụng vào thực tiễn sản xuất Việt Nam Trong trình thử nghiệm phương pháp, luận án sử dụng máy đo với thiết bị cải tiến từ máy đo phân cực kích thích thơng thường, suất lao động chưa cao 166 Vì vậy, để áp dụng phương pháp 3D cần phải có chuyên gia có kinh nghiệm thành thạo việc triển khai phương pháp 2D Mặt khác, thời gian có hạn nên số vùng thử nghiệm hạn chế, chưa mở rộng khoáng sản sulfur đa kim khác Để nhanh chóng đưa phương pháp phân cực kích thích 3D vào thực tiễn sản xuất Việt Nam, nghiên cứu sinh có số kiến nghị sau: Tiếp tục nghiên cứu bổ sung hồn thiện quy trình cơng nghệ đo, xử lý phân tích tài liệu phân cực kích thích 3D đánh giá thăm dị quặng đồng vàng, chì - kẽm Đồng thời mở rộng phạm vi áp dụng phương pháp phân cực kích thích 3D khống sản khác Tổng hợp, xây dựng Atlat phương pháp phân cực kích thích 3D cho khống sản đồng - vàng, chì - kẽm Việt Nam Tiếp tục ứng dụng triển khai kết luận án vào công tác đánh giá thăm dị khống sản đồng - vàng, chì - kẽm sulfur đa kim khác lãnh thổ Việt Nam 167 CÁC CƠNG TRÌNH KHOA HỌC CỦA TÁC GIẢ ĐÃ CƠNG BỐ CĨ LIÊN QUAN ĐẾN NỘI DUNG LUẬN ÁN Trần Bình Trọng (2000), "Khai thác thông tin miền thời gian theo qui luật suy giảm logarit phương pháp phân cực kích thích dịng chiều tìm kiếm khống hố thạch anh - sulfur -vàng", Hội thảo khoa học Địa chất Khoáng sản, Cục Địa chất Khoáng sản Việt Nam, Hà Nội Nguyễn Phương, Trần Bình Trọng, Lương Quang Khang (2002), "Lựa chọn phương pháp đánh giá trạng môi trường phóng xạ mỏ xạ - vùng Tây Bắc - Việt Nam", Báo cáo khoa học HNKH lần thứ 17, Trường đại học Mỏ - Địa chất, Hà Nội Nguyễn Quang Hưng, Trần Bình Trọng, Nguyễn Phương (2002), "Kết ban đầu nghiên cứu môi trường phóng xạ miền trung Việt Nam", Báo cáo khoa học HNKH lần thứ 17, Trường đại học Mỏ - Địa chất, Hà Nội Đào Văn Thịnh, Trần Bình Trọng, n.n.k (2005), "Nội dung điều tra địa chất môi trường tai biến địa chất", Báo cáo khoa học HNKH kỷ niệm 60 năm thành lập ngành địa chất, Cục Địa chất Khoáng sản Việt Nam, Hà Nội Trần Bình Trọng, Vũ Văn Bích (2005), "Một số kết việc áp dụng máy phổ alpha để điều tra trạng mơi trường phóng xạ Việt Nam", Báo cáo khoa học HNKH kỷ niệm 60 năm thành lập ngành địa chất, Cục Địa chất Khoáng sản Việt Nam, Hà Nội Nguyễn Đăng Thành, Trần Bình Trọng (2005), "Một số kết nghiên cứu đặc điểm phân bố quặng chì - kẽm Na Sơn", Báo cáo khoa học HNKH kỷ niệm 60 năm thành lập ngành địa chất, Cục Địa chất Khoáng sản Việt Nam, Hà Nội Nguyễn Văn Nam, Bùi Tất Hợp, Vũ Văn Bích, Trần Bình Trọng (2006), “Nghiên cứu áp dụng phương pháp khí phóng xạ dùng máy phổ alpha RAD7 xác định đứt gãy kiến tạo”, Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất (14), (42006) 168 Trần Bình Trọng , Đồn Thế Hùng (2006), "Đặc trưng suy giảm theo thời gian hệ số phân cực kích thích, quan hệ chúng với quặng hóa đá vây quanh", Báo cáo khoa học HNKH lần thứ 17, Trường đại học Mỏ - Địa chất, Hà Nội Trần Bình Trọng, Trịnh Đình Huấn, Nguyễn Phương (2007), "Điều tra trạng mơi trường phóng xạ tụ khống Đơng Pao, Thèn Sin – Tam Đường (Lai Châu), Mường Hum (Lào Cai), Yên Phú (Yên Bái), Thanh Sơn (Phú Thọ), An Điềm, Ngọc Kinh – Sườn Giữa (Quảng Nam)", Tạp chí Địa chất, Loạt A, (298), 1-2/2007) 10 Trần Bình Trọng, Trịnh Đình Huấn (2007), “Đánh giá ảnh hưởng mơi trường phóng xạ mỏ đất – phóng xạ n Phú – n Bái” Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất (18), 4-2007 11 Trần Bình Trọng (2007), Nghiên cứu ứng dụng phương pháp phân cực kích thích 3D (dịng chiều) điều tra đánh giá khoáng sản kim loại, áp dụng thử nghiệm số vùng Việt Nam, Báo cáo đề tài NCKH, Bộ Tài nguyên Môi trường, Hà Nội 169 TÀI LIỆU THAM KHẢO Bùi Xuân Ánh, (2007), Đánh giá triển vọng quặng đồng khoáng sản khác khu vực Tả Phời, thị xã Cam Đường, tỉnh Lào Cai, Báo cáo địa chất, Hà Nội Tạ Việt Dũng, Đỗ Hải Dũng, Trần Tất Thắng (1995), "Tài nguyên khoáng sản kim loại Cu, Pb- Zn Việt Nam" Báo cáo HNKH Địa chất VN lần thứ III : 31-44, Hà Nội Bùi Tất Hợp, Trần Bình Trọng n.n.k (2000), Báo cáo đánh giá khoáng sản vàng dải Thèn Sin - Phong Thổ - Lai Châu, Lưu trữ địa chất, Hà Nội Dương Đức Kiêm (2002), Nghiên cứu kiến tạo sinh khoáng Bắc Bộ, Lưu trữ địa chất, Hà Nội Nguyễn Ngọc Loan (1996), Đánh giá triển vọng thân quặng sulfur đa kim theo đặc trưng dị thường phân cực kích thích dịng chiều, Luận văn phó tiến sỹ, Viện Địa chất Khống sản, Hà Nội Nguyễn Quang Luật (2003), Điều kiện thành tạo quy luật phân bố quặng hóa đồng, vàng Tây Bắc Việt Nam bối cảnh kiến tạo mảng, Trường Đại học Mỏ-Địa chất Nguyễn Quang Luật (2005), Điều kiện thành tạo quy luật phân bố quặng đa kim chì-kẽm Đơng Bắc Việt Nam, Trường Đại học Mỏ-Địa chất Nguyễn Nghiêm Minh (1997), "Một số nét tiềm khoáng sản Tây Bắc Việt Nam", Tạp chí Địa chất, A/Phụ trương, 4-9, Hà Nội Nguyễn Trọng Nga (1996), "Tổ hợp phương pháp địa vật lý tìm số quặng ẩn điển hình Việt Nam", Báo cáo khoa học HNKH lần thứ 12, trường đại học Mỏ - Địa chất, Hà Nội 10 Nguyễn Trọng Nga (1997), Giáo trình thăm dị điện phân giải cao, Trường đại học Mỏ - Địa chất, Hà Nội 11 Trần Thanh Rỹ (2003), Hoàn thiện công nghệ khai thác thông tin địa vật lý để điều tra phát thân quặng sunfua đa kim vùng Tòng Mụ - Chợ Đồn Bắc Kạn, Luận án Tiến sỹ địa chất, Trường đại học Mỏ - Địa chất, Hà Nội 170 12 Mai Thanh Tân - “Giáo trình Địa Vật Lý đại cương” - Đại học Mỏ - Địa Chất, năm 2004 13 Nguyễn Đăng Thành, Trần Bình Trọng (2005), "Một số kết nghiên cứu đặc điểm phân bố quặng chì - kẽm Na Sơn", Báo cáo khoa học HNKH kỷ niệm 60 năm thành lập ngành địa chất, Cục Địa chất Khống sản Việt Nam, Hà Nội 14 Nguyễn Hữu Trí (1998), Sách tra cứu tính chất vật lý đá số loại quặng chủ yếu Việt Nam, Cục Địa chất Khoáng sản Việt Nam 15 Phùng Quốc Trị (2007), Báo cáo bước địa chất Đánh giá triển vọng quặng chì kẽm vùng Ba Xứ,Yên Sơn, tỉnh Tun Quang, Liên đồn Địa chất Đơng Bắc 16 Trần Văn Trị, n.n.k (2000), Tài nguyên khoáng sản Việt Nam Bộ Công nghiệp, Hà Nội 17 Nguyễn Như Trung (1996), Nghiên cứu áp dụng phương pháp sai phân hữu hạn để xử lý phân tích số liệu đo sâu điện trở suất môi trường bất đồng ngang, Luận án tiến sĩ, Trường đại học Mỏ - Địa chất, Hà Nội 18 Nguyễn Như Trung (2002), Giáo trình cao học Phương pháp thăm dò điện dòng chiều môi trường bất đồng ngang, Trường đại học Mỏ Địa chất, Hà Nội 19 Trần Bình Trọng (2000), "Khai thác thông tin miền thời gian theo qui luật suy giảm logarit phương pháp phân cực kích thích dịng chiều tìm kiếm khống hố thạch anh - sulfur -vàng", Hội thảo khoa học Địa chất Khoáng sản, Cục Địa chất Khoáng sản Việt Nam, Hà Nội 20 Trần Bình Trọng , Đồn Thế Hùng (2006), "Đặc trưng suy giảm theo thời gian độ phân cực kích thích, quan hệ chúng với quặng hóa đá vây quanh", Báo cáo khoa học HNKH lần thứ 17, Trường đại học Mỏ - Địa chất, Hà Nội 171 Tiếng Anh: 21 Denne, Rebecca, (2001), "A new survey design for 3D IPinversion modelling at Copper Hill", presented at Conference of the Australian Society of Exploration Geophysicists, Brisbane, Australia 22 Dey A Morrison H.F (1979), "Resistrivity modeling for arbitrarily shaped two - dimensional structures", Geophysical Prospecting 27, P 106 - 136 23 Dey A Morrison H.F (1979), "Resistrivity modeling for arbitrarily shaped three dimensional structures", Geophysics; Vol 44, 4, P 753-780 24 Loke M.H 2004., 3-D resistivity & IP forward modeling using the finitedifference and finite-element methods., (CDROM) 25 Loke M.H 2004, Tutorial : 2-D and 3-D electrical imaging surveys http://www.geoelectrical.com/coursenotes.zip 26 Loke M.H (2000), A practical guide to 2-D and 3-D surveys, http://www.agiusa.com 27 John Coggon, (1995), Induced polarisation methods, WA School of Mines 28 Kowalczyk, Peter (2000), "Inversion of induced polarization data from Donlin Creek, Alaska", 70th Ann Internat Mtg., Soc Expl Geophys., Expanded Abstracts 29 Li, Y and Oldenburg, D.W (2000), "3D inversion of induced polarization data", Geophysics, Vol65, #6 30 Medeiros W.E.; Silva J.B.C (1996), "Geophysical inversion using approximate equality restraints", Geophysics; 61/6; 1678-1688 USA: Society of Exploration Geophysicists 31 Mohan N.L; Babu L.A (1995), "An analysis of 3-D analytic signal", Geophysics; 60/2; 531-536 USA: Society of Exploration Geophysicists 32 Haber E.; Ascher U.M.; Oldenburg D.W (2004), "Inversion of 3D electromagnetic data in frequency and time domain using an inexact all-atonce approach", Geophysics; 69/5; 1216-1228 USA: Society of Exploration Geophysicists 172 33 Prem V Sharma, (1997), Environmental and engineering geophysics, Cambridge University Prees 34 Smith J.T (1996), "Conservative modeling of 3-D electromagnetic fields, Part I: properties and error analysis", Geophysics; 61/5; 1308-1318 USA: Society of Exploration Geophysicists 35 Telford W M., L Geldart P., Sheriff R.E., (2001), Applied Geophysics Second Edition, Cambridge University Preess 36 Wijns, C., and Kowalczyk, P., (2004), "Human interaction in geophysical inversion", ASEG 17th Conference & Exhibition, Sydney, Extended Abstracts 37 Zhao S.; Yedlin M.J (1996), "Some refinements on the finite-difference method for 3-D DC resistivity modeling", Geophysics; 61/5; 1301-1307 USA: Society of Exploration Geophysicists Tiếng Nga: 38 Комаров В А (1980) Электроразведкa Методом поляризации, Издательство Недра, Ленинград вызванной 173 PHỤ LỤC Chương trình tính tham số suy giảm độ phân cực kích thích * Trần Bình Trọng - 6/2007 * Nghiên cứu suy giảm theo thời gian đo điện phân cực kích thích loại máy DWJ-2 * Mơ hình lựa chọn a.Ln(t+1)+b, b.e^(at), b.(t+1)^a, b/(1+a.t), tích phân (Ch), tốc độ suy giảm (Ch_e1, mV/s) * Số cửa sổ tham gia tính tuỳ ý *set defa to C:\ddll1\luanan clea clos data set talk off @ 2,2 say "Chọn File số liệu : " font ".VntimeH", 12 tf=getfile("DBF") @ 2,2 say "File số liệu sử dụng: " + tf font ".VntimeH", 12 use &tf * Tạo them cột suy giam copy to tg copy stru exte to tg1 use tg1 inde on field_name to t1 seek "AL" ?found() if !found() clea set inde to appe blan repl field_name with "Al" repl field_type with "N" repl field_Len with 10 repl field_dec with appe blan repl field_name with "Bl" repl field_type with "N" repl field_Len with 10 repl field_dec with appe blan repl field_name with "Rxyl" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Tl" repl field_type with "N" repl field_Len with 174 repl field_dec with appe blan repl field_name with "Ae" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Be" repl field_type with "N" repl field_Len with 10 repl field_dec with appe blan repl field_name with "Te" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Rxye" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Ax" repl field_type with "N" repl field_Len with 10 repl field_dec with appe blan repl field_name with "Bx" repl field_type with "N" repl field_Len with 10 repl field_dec with appe blan repl field_name with "Tx" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Rxyx" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Ah" repl field_type with "N" repl field_Len with 10 175 repl field_dec with appe blan repl field_name with "Bh" repl field_type with "N" repl field_Len with 10 repl field_dec with appe blan repl field_name with "Th" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Rxyh" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Ch" repl field_type with "N" repl field_Len with repl field_dec with appe blan repl field_name with "Ch_E1" repl field_type with "N" repl field_Len with repl field_dec with @ 4,20 say "Đã thực xong thêm trường số liệu" font ".Vntime", 12 else @ 4,20 say "Không thực hiên thêm trường số liệu " font ".Vntime", 12 endi use crea &tf from tg1 sele use &tf appe from tg dele file tg1.dbf sele use tg * =================== Bat dau tinh khong chuan hoa ve Eta1 ============ sele use kieuthz disp all 176 kieutz="kieu1 " sbc=10 @ 22,2 say "Kiểu thời gian cài đặt (Mod 64: 1; mẫu tham số: 2; Trích từ đồ thị: 10; mod 62: 3/4/5/6/7/8/9)" font ".Vntime", 12 @ 24,2 say "Lựa chọn kiểu thời gian" font ".Vntime", 11 get kieutz font ".Vntime", 11 read dime tbt1[7,sbc] for i=1 to 10 go i tbt1[1,i]=&kieutz endf use sele sbd=1 sbc=8 scs=sbc-sbd+1 @ 24,2 say "Lựa chọn cửa sổ " font ".Vntime", 12 @ 26,2 say "Từ Eta" font ".Vntime", 11 get sbd font ".Vntime", 11 @ 26,18 say " đến Eta" font ".Vntime", 11 get sbc font ".Vntime",11 read for j=sbd to sbc tbt1[2,j]=log(tbt1[1,j]) && log thời gian endf * Tính tham số trung bình, phương sai thời gian tz11=0 && trung bình thời gian thường tz12=0 && trung bình thời gian log tz21=0 && độ lệch thời gian thường tz22=0 && độ lệch thời gian log for j=sbd to sbc tz11=tz11+tbt1[1,j]/scs tz12=tz12+tbt1[2,j]/scs endf for j=sbd to sbc tz21=tz21+(tbt1[1,j]-tz11)**2 tz22=tz22+(tbt1[2,j]-tz12)**2 endf tz21=sqrt(tz21/scs) tz22=sqrt(tz22/scs) ? "tz11= ",tz11 ? "tz12= ",tz12 ? "tz21= ",tz21 ? "tz22= ",tz22 for j=sbd to sbc tbt1[5,j]=tbt1[1,j]-tz11 && Tz thường tbt1[6,j]=tbt1[2,j]-tz12 && Tz Log endf 177 coun to nn for i=1 to nn go i rxy1=0 && X1,Y1 Hệ số tương quan Eta thường log(tz) [3,6] rxy2=0 && X2,Y2 Hệ số tương quan log(Eta) thường tz [4,5] rxy3=0 && X3,Y3 Hệ số tương quan Eta thường tz*Eta [3,7] rxy4=0 && X4,Y4 Hệ số tương quan log(Eta) log(tz) [4,6] etb1=0 && Y1tb = Eta trung bình thường etb2=0 && Y2tb = Eta trung bình log etb3=0 && Y3tb = t.Eta trung bình for j=sbd to sbc ttr="eta"+alltrim(str(j)) tbt1[3,j]=&ttr && Y1 (hàng 3, 6) - Mơ hình log tbt1[4,j]=log(tbt1[3,j]) && Y2=Log(Y) hàng (4, 5)- Mơ hình E xponential tbt1[7,j]=tbt1[1,j]*tbt1[3,j] && X3 mo hinh hybecbolit etb1=etb1+tbt1[3,j]/scs etb2=etb2+tbt1[4,j]/scs etb3=etb3+tbt1[7,j]/scs endf edl1=0 && Y1dl = Eta độ lệch thường edl2=0 && Y2dl = Eta độ lệch log edl3=0 && Y3dl = Eta*tz độ lệch for j=sbd to sbc edl1=edl1+(tbt1[3,j]-etb1)**2 edl2=edl2+(tbt1[4,j]-etb2)**2 edl3=edl3+(tbt1[7,j]-etb3)**2 endf edl1=sqrt(edl1/scs) edl2=sqrt(edl2/scs) edl3=sqrt(edl3/scs) for j=sbd to sbc rxy1=rxy1+tbt1[6,j]*(tbt1[3,j]-etb1) rxy2=rxy2+tbt1[5,j]*(tbt1[4,j]-etb2) rxy3=rxy3+(tbt1[3,j]-etb1)*(tbt1[7,j]-etb3) rxy4=rxy4+tbt1[6,j]*(tbt1[4,j]-etb2) endf rxy1=rxy1/tz22/edl1/scs rxy2=rxy2/tz21/edl2/scs rxy3=rxy3/edl1/edl3/scs rxy4=rxy4/tz22/edl2/scs tct1=abs(rxy1)/(1-rxy1**2)*sqrt(scs) tct2=abs(rxy2)/(1-rxy2**2)*sqrt(scs) aa1=rxy1*edl1/tz22 bb1=etb1-aa1*tz12 178 aa2=rxy2*edl2/tz21 bb2=etb2-aa2*tz11 aa2=abs(1/aa2) bb2=exp(bb2) aa3=rxy3*edl1/edl3 bb3=etb1-aa3*etb3 aa3=-aa3 aa4=rxy4*edl2/tz22 bb4=etb2-aa4*tz12 bb4=exp(bb4) aa5=0 && he so tong nap (Chargeability) aa6=0 for j=sbd to sbc-1 aa5=aa5+abs(tbt1[3,j]+tbt1[3,j+1])*abs(tbt1[1,j+1]-tbt1[1,j])/2 aa6=aa6+abs(tbt1[3,j]-tbt1[3,j+1])/abs(tbt1[1,j+1]-tbt1[1,j]) endf repl al with aa1 repl bl with bb1 repl tl with tct1 repl ae with aa2 repl be with bb2 repl te with tct2 repl ah with aa3 repl bh with bb3 repl ax with aa4 repl bx with bb4 repl ch with aa5 repl ch_e1 with aa6*1000 repl rxyl with rxy1 repl rxye with rxy2 repl rxyh with rxy3 repl rxyx with rxy4 endf * ============== Ket thuc tinh khong chuan hoa ============= sele use dele file tg.dbf set talk on sele brow fiel eta1,ch,ch_e1,al,bl,rxyl,ae,be,rxye,ax,bx,rxyx,ah,bh,rxyh RETU ... BÌNH TRỌNG LUẬN CHỨNG HIỆU QUẢ ÁP DỤNG PHƯƠNG PHÁP PHÂN CỰC KÍCH THÍCH TRONG TÌM KIẾM QUẶNG ĐỒNG VÀNG, CHÌ - KẼM Ở VIỆT NAM Chuyên ngành: Các phương pháp Địa vật lý Địa hóa học tìm kiếm thăm dị... dựng luận khoa học thực tiễn chứng minh hiệu áp dụng phương pháp phân cực kích thích phân cực kích thích 3D tìm kiếm quặng đồng - vàng, chì - kẽm Việt Nam Đối tượng phạm vi nghiên cứu - Hiệu áp dụng. .. dụng phương pháp phân cực kích thích tìm kiếm quặng đồng - vàng, chì - kẽm Việt Nam - Phạm vi nghiên cứu luận án là: mỏ đồng - vàng, chì - kẽm Việt Nam, tập trung chủ yếu miền Bắc Việt Nam 11