ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Dương Thị Hoài Thu XÁC ĐỊNH PHÔNG KHU VỰC TRƯỜNG TRỌNG LỰC TRÊN MỘT SỐ BỂ TRẦM TÍCH KAINOZOI THỀM LỤC ĐỊA VIỆT NAM LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – Năm 2012 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Dương Thị Hoài Thu XÁC ĐỊNH PHÔNG KHU VỰC TRƯỜNG TRỌNG LỰC TRÊN MỘT SỐ BỂ TRẦM TÍCH KAINOZOI THỀM LỤC ĐỊA VIỆT NAM Chuyên ngành: Vật lý địa cầu Mã số: 60 44 15 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC PGS.TS ĐỖ ĐỨC THANH Hà Nội – Năm 2012 LỜI CẢM ƠN Trong thời gian thực luận văn tốt nghiệp em gặp nhiều khó khăn, bỡ ngỡ, nhờ động viên, giúp đỡ nhiều người em hoàn thành luận văn Do thời gian thực luận văn trình độ thân nhiều hạn chế nên thực luận văn chắn không tránh khỏi sai sót, nên mong đóng góp thầy cô người quan tâm Luận văn thực hướng dẫn PGS.TS Đỗ Đức Thanh Thầy tận tình bảo, hướng dẫn, giúp đỡ em trả lời nhiều thắc mắc đề tài suốt thời gian làm luận văn Em xin chân thành cảm ơn thầy! Em xin cảm ơn thầy cô môn Vật lý địa cầu trường ĐH Khoa học Tự nhiên truyền đạt cho em kiến thức chuyên ngành vô quý báu suốt thời gian học môn, tạo tiền đề cho em làm luận văn Lời tiếp theo, xin cảm ơn người bạn, anh chị đồng hành, giúp đỡ em trình tìm tài liệu, trao đổi học thuật để em hoàn thành luận văn cách tốt Một lần xin chân thành cảm ơn giúp đỡ quý báu trên! Hà Nội, ngày 28 tháng 11 năm 2012 Học viên Dương Thị Hoài Thu MỤC LỤC MỞ ĐẦU CHƯƠNG PHƯƠNG PHÁP XÁC ĐỊNH DỊ THƯỜNG TRỌNG LỰC ĐỐI VỚI CÁC VẬT THỂ CÓ DẠNG HÌNH HỌC ĐỀU ĐẶN 10 1.1.Những khái niệm .10 1.2 Các biểu thức tích phân tổng quát đạo hàm trọng lực 12 1.3 Bài toán thuận cho vật thể có dạng hình học 15 1.3.1 Hình cầu điểm vật chất 15 Hình 1.3: Xác định đạo hàm vật thể hình cầu 16 Hình 1.4: Trường trọng lực hình cầu .17 1.3.2 Thanh vật chất nằm ngang, hình trụ tròn nằm ngang .17 1.3.3 Nửa mặt phẳng vật chất nằm ngang 18 Hình 1.6: Xác định đạo hàm nửa mặt phẳng vật chất nằm ngang 19 Hình 1.7: Trường trọng lực nửa mặt phẳng vật chất nằm ngang 20 1.3.4 Hình hộp vuông góc 20 1.3.5 Lăng trụ thẳng đứng 21 1.3.6 Bậc thẳng đứng 22 Hình 1.8 a, 1.8 b: Bậc thẳng đứng 22 1.3.7 Bậc nghiêng .24 CHƯƠNG XÁC ĐỊNH HIỆU ỨNG TRỌNG LỰC THEO PHƯƠNG 28 PHÁP GIẢI TÍCH VÀ PHƯƠNG PHÁP PHỔ .28 2.1 Xác định dị thường trọng lực theo phương pháp giải tích 28 2.1.1 Xác định dị thường trọng lực hình trụ tròn nằm ngang 28 2.1.2 Xác định dị thường trọng lực vật thể 3D có dạng hình lăng trụ đứng 29 Hình 2.2: Mô hình lăng trụ chiều 30 2.1.3 Xác định dị thường trọng lực bể trầm tích .30 2.2 Xác định dị thường trọng lực theo phương pháp phổ 31 2.2.1 Xác định dị thường trọng lực vật thể 3D có dạng hình lăng trụ đứng 31 2.2.2 Xác định dị thường trọng lực bể trầm tích .32 2.2.3 Nâng cao độ xác việc tính dị thường trọng lực miền tần số phương pháp “trượt mẫu” (Shiftampling) 33 CHƯƠNG XÁC ĐỊNH PHÔNG KHU VỰC TRÊN MỘT SỐ BỂ TRẦM .36 TÍCH KAINOZOI THỀM LỤC ĐỊA VIỆT NAM .36 3.1 Xây dựng đồ Bughe khu vực thềm lục địa Việt Nam vùng Biển Đông lân cận 36 3.1.1 Kết tính toán .36 .37 Hình 3.1: Bản đồ dị thường Fye tỉ lệ 1: 200.000 khu vực thềm lục địa 37 Việt Nam vùng Biển Đông kế cận 37 Việt Nam vùng Biển Đông kế cận 38 Hình 3.3: Bản đồ hiệu chỉnh địa hình đáy biển tỉ lệ 1: 200.000 khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận .39 Hình 3.4: Bản đồ dị thường trọng lực Bughe tỉ lệ 1: 200.000 khu vực thềm 40 lục địa Việt Nam vùng Biển Đông kế cận 40 Hình 3.5: Trường dị thường trọng lực Bughe Biển Đông tỉ lệ 1:1.000.000 41 (Nguồn Viện Địa chất Địa Vật lý biển) .41 3.1.2 Nhận xét 42 3.2 Xác định phông khu vực trường trọng lực số bể trầm tích Kainozoi thềm lục địa Việt Nam .44 3.2.1 Phông khu vực phần đông nam thềm lục địa Việt nam 44 3.2.2 Phông khu vực dị thường dư bể trầm tích Cửu long Nam Côn sơn 45 Hình 3.6: Bản đồ dị thường trọng lực Bughe tỉ lệ 1: 200.000 khu vực 46 Đông Nam thềm lục địa Việt Nam 46 Hình 3.7: Bản đồ trường phông tỉ lệ 1: 200.000 khu vực 47 Đông Nam thềm lục địa Việt Nam 47 Hình 3.8: Bản đồ trường phông khu vực Cửu Long tỉ lệ 1:200.000 48 48 Hình 3.9: Bản đồ trường dư khu vực Cửu Long tỉ lệ 1:200.000 48 .49 Hình 3.10: Bản đồ trường phông khu vực Nam Côn Sơn tỉ lệ 1:200.000 49 Hình 3.11: Bản đồ trường dư khu vực Nam Côn Sơn tỉ lệ 1:200.000 49 KẾT LUẬN .49 TÀI LIỆU THAM KHẢO 50 DANH MỤC CÁC HÌNH VẼ MỞ ĐẦU CHƯƠNG PHƯƠNG PHÁP XÁC ĐỊNH DỊ THƯỜNG TRỌNG LỰC ĐỐI VỚI CÁC VẬT THỂ CÓ DẠNG HÌNH HỌC ĐỀU ĐẶN 10 Hình 1.3: Xác định đạo hàm vật thể hình cầu 16 Hình 1.4: Trường trọng lực hình cầu .17 Hình 1.6: Xác định đạo hàm nửa mặt phẳng vật chất nằm ngang 19 Hình 1.7: Trường trọng lực nửa mặt phẳng vật chất nằm ngang 20 Hình 1.8 a, 1.8 b: Bậc thẳng đứng 22 CHƯƠNG XÁC ĐỊNH HIỆU ỨNG TRỌNG LỰC THEO PHƯƠNG 28 PHÁP GIẢI TÍCH VÀ PHƯƠNG PHÁP PHỔ .28 Hình 2.2: Mô hình lăng trụ chiều 30 CHƯƠNG XÁC ĐỊNH PHÔNG KHU VỰC TRÊN MỘT SỐ BỂ TRẦM .36 TÍCH KAINOZOI THỀM LỤC ĐỊA VIỆT NAM .36 .37 Hình 3.1: Bản đồ dị thường Fye tỉ lệ 1: 200.000 khu vực thềm lục địa 37 Việt Nam vùng Biển Đông kế cận 37 Việt Nam vùng Biển Đông kế cận 38 Hình 3.3: Bản đồ hiệu chỉnh địa hình đáy biển tỉ lệ 1: 200.000 khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận .39 Hình 3.4: Bản đồ dị thường trọng lực Bughe tỉ lệ 1: 200.000 khu vực thềm 40 lục địa Việt Nam vùng Biển Đông kế cận 40 Hình 3.5: Trường dị thường trọng lực Bughe Biển Đông tỉ lệ 1:1.000.000 41 (Nguồn Viện Địa chất Địa Vật lý biển) .41 Hình 3.6: Bản đồ dị thường trọng lực Bughe tỉ lệ 1: 200.000 khu vực 46 Đông Nam thềm lục địa Việt Nam 46 Hình 3.7: Bản đồ trường phông tỉ lệ 1: 200.000 khu vực 47 Đông Nam thềm lục địa Việt Nam 47 Hình 3.8: Bản đồ trường phông khu vực Cửu Long tỉ lệ 1:200.000 48 48 Hình 3.9: Bản đồ trường dư khu vực Cửu Long tỉ lệ 1:200.000 48 .49 Hình 3.10: Bản đồ trường phông khu vực Nam Côn Sơn tỉ lệ 1:200.000 49 Hình 3.11: Bản đồ trường dư khu vực Nam Côn Sơn tỉ lệ 1:200.000 49 KẾT LUẬN .49 TÀI LIỆU THAM KHẢO 50 MỞ ĐẦU Địa vật lý thăm dò lĩnh vực khoa học nghiên cứu Trái đất Có nhiều phương pháp địa vật lý sử dụng rộng rãi thăm dò như: phương pháp thăm dò trọng lực, phương pháp thăm dò điện, phương pháp phóng xạ Dựa vào việc nghiên cứu trường trọng lực trường từ trái đất thông qua phép đo tiến hành mặt đất, biển, không giúp việc nghiên cứu trái đất rõ ràng hơn, giúp ta khám phá hình dạng cấu trúc vỏ trái đất, tìm kiếm mỏ khoáng sản phục vụ cho tồn phát triển nhân loại Trong năm gần đây, với phát dầu đá móng bể trầm tích thuộc phần Biển Đông, việc xác định độ sâu tới móng kết tinh, việc nghiên cứu cấu trúc mà trước hết nghiên cứu phân bố dị thường Bughe đặc biệt trở nên quan trọng thu hút quan tâm nhiều nhà Địa Vật lý nước Để góp phần giải vấn đề này, luận văn, tiến hành xây dựng đồ Bughe tỉ lệ 1:200.000 sở việc giải toán thuận theo phương pháp giải tích tính hiệu ứng trọng lực gây lăng trụ thẳng đứng Kết thu sử dụng để xác định phông khu vực từ xác định phần trường dư gây số bể trầm tích thuộc thềm lục địa Việt nam Chương trình tính toán viết ngôn ngữ FORTRAN MATLAB Luận văn chia làm ba chương sau: - Chương 1: Phương pháp xác định dị thường trọng lực vật thể có dạng hình học đặn - Chương 2: Xác định hiệu ứng trọng lực theo phương pháp giải tích phương pháp phổ - Chương 3: Xác định phông khu vực số bể trầm tích Kainozoi thềm lục địa Việt nam CHƯƠNG PHƯƠNG PHÁP XÁC ĐỊNH DỊ THƯỜNG TRỌNG LỰC ĐỐI VỚI CÁC VẬT THỂ CÓ DẠNG HÌNH HỌC ĐỀU ĐẶN 1.1.Những khái niệm Sau tu chỉnh số liệu đo đạc máy trọng lực, ta thành lập đồ đồ thị đạo hàm bậc nhất, bậc hai trọng lực (dị thường trọng lực nói chung) Giải thích địa chất dị thường trọng lực bao gồm phân tích quy luật phân bố mặt đất (hoặc gần mặt đất) mối liên hệ để giải nhiệm vụ khác Giải thích địa chất dị thường trọng lực hình thành sau: Dựa vào số liệu trọng lực đo đạc số liệu địa chất, địa vật lý sẵn có, vào kinh nghiệm giải thích trọng lực vùng tương đương, ta đưa kết luận địa chất vùng cho trước tương ứng với nhiệm vụ địa chất đề Nhiệm vụ giải thích dị thường trọng lực phân loại hai hình thức: Phân tích định tính định lượng Khi giải thích định tính cần xác định: - Các yếu tố địa chất chắn ảnh hưởng lên trường trọng lực trường vật lý khác (nếu phương pháp Địa vật lý khác áp dụng) 10 Hình 3.1: Bản đồ dị thường Fye tỉ lệ 1: 200.000 khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận 37 Hình 3.2: Bản đồ địa hình đáy biển tỉ lệ 1: 200.000 khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận 38 Hình 3.3: Bản đồ hiệu chỉnh địa hình đáy biển tỉ lệ 1: 200.000 khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận 39 Hình 3.4: Bản đồ dị thường trọng lực Bughe tỉ lệ 1: 200.000 khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận 40 Hình 3.5: Trường dị thường trọng lực Bughe Biển Đông tỉ lệ 1:1.000.000 (Nguồn Viện Địa chất Địa Vật lý biển) 41 3.1.2 Nhận xét - So với kết tính trường dị thường trọng lực Bughe Biển Đông chương trình điều tra nghiên cứu biển cấp nhà nước KC-09-02 (hình 3.5) ta thấy đồ trường trọng lực Bughe thu theo phương pháp tính hiệu chỉnh lăng trụ kể (hình 3.4) có kết phù hợp - Trường trọng lực Biển Đông nhìn chung có cấu trúc phức tạp phân dị, có hướng cấu trúc áp đảo Đông Bắc – Tây Nam kinh tuyến, đặc biệt rõ loạt dải bất thường bậc thang với biên độ hàng trăm mgal kéo dài hàng nghìn km theo hướng Bắc Đông Bắc Các dải bất thường bậc thang lớn chạy song song tạo thành nhóm quan sát thấy khu vực sau: * Hướng Đông Bắc gồm có: - Dải chạy qua phía Tây Bắc quần đảo Hoàng Sa lên đến gần Đài Loan - Dải chạy ven bờ biển đảo Kalimantan Palawan - Ba dải chạy dọc theo vùng trung tâm Biển Đông - Một dải chạy men theo phía Tây Bắc quần đảo Trường Sa * Hướng kinh tuyến bật lên dải bất thường bậc thang chạy dọc theo ven biển Trung Bộ Việt Nam dải chạy song song ven biển phía Tây Philippin Các nhóm bất thường bậc thang (gradien) hướng Tây Bắc – Đông Nam (cắt ngang trục Biển Đông) hướng vĩ tuyến thấy thể rõ biên độ không lớn độ dài hạn chế, bị đứt đoạn đan chéo với dải gradien hướng chủ đạo Đông Bắc – Tây Nam Nổi bật nhóm dải gradien hướng vĩ tuyến phần Đông Bắc Biển Đông, dải bất thường bậc thang hướng vĩ tuyến chạy từ phía Nam đảo Hải Nam qua khu vực trung tâm Biển Đông sát vĩ độ 160N Dải bất thường bậc thang nằm phía Nam quần đảo Trường Sa trải theo vĩ 42 độ 60N số dải bất thường hướng vĩ tuyến phần phía Đông Biển Đông cắt vuông góc với dải gradien hướng kinh tuyến ven biển Philippin Các bất thường bậc thang hướng Tây Bắc – Đông Nam biên độ không lớn, thường bị gián đoạn biến dạng cắt qua dải lớn có hướng Đông Bắc sườn trung tâm Biển Đông Những vùng có cấu trúc trường bình ổn tương đối so với vùng rìa trung tâm Biển Đông khu vực thềm lục địa Đông Nam Trung Quốc Phần Đông Bắc từ vĩ độ 21 0N trở lên, toàn vùng vịnh Thái Lan thềm lục địa Đông Nam Việt Nam Bắc Indonexia phía Tây dải Gradien kinh tuyến 1090 – 1100E Ở vùng nói giá trị bất thường Bughe biến đổi giới hạn từ ± 10, ± 20 đến + 50, +60 mgal, số vùng trường có giá trị âm cỡ vài ba chục mgal Bức tranh bất thường biến động từ từ, dải gradien mạnh kéo dài, hướng cấu trúc thể tương đối rõ vùng vịnh Thái Lan, thềm lục địa Đông Nam Việt Nam thềm lục địa phía Bắc Indonexia Ở vùng quần đảo Hoàng Sa trường trọng lực có cấu trúc tương đối phức tạp, phân dị thành bất thường dạng hẹp kéo dài Hướng cấu trúc Hoàng Sa áp đảo hướng vĩ tuyến Ở quần đảo Hoàng Sa giá trị trọng lực Bughe biến đổi từ +30 đến +120 mgal Các vùng rìa phía Tây Bắc, Nam phía Đông có dải bất thường gradien lớn bao bọc Vùng trung tâm Biển Đông kéo dài theo phương Đông Bắc có đặc điểm cấu trúc trường trọng lực khác hẳn đặc trưng, giá trị bất thường Bughe lớn, biến đổi giới hạn từ +200 đến +300, +320 mgal Toàn vùng giới hạn dải bất thường bậc thang biên độ lớn gồm đường đẳng mức chạy sít Bên vùng cấu trúc trường tương đối bình ổn tạo nên số bất thường nhỏ có dạng khối đặn Ở phía Đông Bắc hình thành dải bất thường bậc thang cắt theo hướng vĩ tuyến Tây Bắc – Đông Nam Các vùng phía đông thềm lục địa miền Trung Việt Nam phía Đông Bắc 43 quần đảo Hoàng Sa trường trọng lực có cấu trúc phức tạp tượng đan chéo dải gradien mạnh theo nhiều hướng khác Ở đường đẳng trị trọng lực biến đổi đột ngột, phức tạp tạo nên tranh phân dị rõ hướng cấu trúc chủ đạo Vùng thềm phía Tây Bắc đảo Palawan nơi đan chéo dải bất thường gradien theo hướng khác nhau: Đông Bắc, kinh tuyến vĩ tuyến, tạo nên tranh bất thường phức tạp Giá trị bất thường Bughe vùng biến đổi từ +50 mgal đến +180, +200 mgal Cùng với vùng quần đảo Trường Sa Hoàng Sa vùng bất thường tạo nên đai chuyển tiếp bao quanh vùng trung tâm, có đặc điểm cấu trúc phức tạp giá trị bất thường tương đối lớn, biến đổi giới hạn rộng 3.2 Xác định phông khu vực trường trọng lực số bể trầm tích Kainozoi thềm lục địa Việt Nam 3.2.1 Phông khu vực phần đông nam thềm lục địa Việt nam Để xác định phông khu vực từ tính dị thường dư số bể trầm tích Kainozoi thuộc thềm lục địa Việt Nam, phần này, xem phần trường khu vực phạm vi phần Đông Nam thềm lục địa Việt Nam có dạng phông bậc hai: ( ) F ( x, y ) = a ( x − x0 ) + b ( x − x0 ) ( y − y0 ) + c ( y − y0 ) + d x − x + e ( y − y0 ) + f (3.1) 2 Trong Q(x0, y0) điểm chọn làm gốc tọa độ Thực chất việc xác định phông khu vực bậc theo phương pháp bình phương tối thiểu phải xác định hệ số khai triển a, b, c, d, e, f cho: N ∑ P ∆g ( x, y ) − F ( x, y )  k =1 k k k = (3.2) Trong ∆g k ( x, y ) kí hiệu dị thường trọng lực quan sát điểm thứ k, N số điểm quan sát lấy, Pk gọi hàm trọng số Nó xác 44 γ định sau:  R − d2  Pk =  k ÷  dk −η  (3.3) với R bán kính khu vực tính, d k khoảng cách từ gốc tọa độ tới điểm quan sát thứ k η γ hệ số Ở đây, việc xây dựng phông bậc hai phạm vi phần Đông Nam thềm lục địa Việt Nam thực dựa vào đồ dị thường Bughe khu vực Đông Nam thềm lục địa Việt Nam phạm vi kinh độ từ 105 đến 1100; vĩ độ từ 50 đến 110 (hình 3.6) tách từ tờ đồ Bughe khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận (hình 3.4) Từ 276x331 giá trị trọng lực Bughe khu vực này, giải phương trình (3.2) theo phương pháp bình phương tối thiểu ngôn ngữ lập trình Matlab, ta kết sau hệ số phông bậc hai khu vực: a = 1.2056181e-005; b = 3.4920913e-005; c = -1.0053530e-004; d = 4.2472383e-002; e = 1.5007690e-001; f = -6.7869089e+001 Kết trình bày hình vẽ 3.7 3.2.2 Phông khu vực dị thường dư bể trầm tích Cửu long Nam Côn sơn Sau đó, từ hình vẽ tiến hành tách để lấy riêng trường phông khu vực phạm vi bể trầm tích Cửu Long, có kinh độ từ 106 055' đến 108040'; vĩ độ từ 9015' đến 10030', bể trầm tích Nam Côn Sơn có kinh độ từ 106 030' đến 109020'; vĩ độ từ 6025' đến 9015' từ tính phần trường dư gây bể trầm tích Kết tính toán trình bày hình vẽ: 3.8, 3.9, 3.10 3.11 45 Hình 3.6: Bản đồ dị thường trọng lực Bughe tỉ lệ 1: 200.000 khu vực Đông Nam thềm lục địa Việt Nam 46 Hình 3.7: Bản đồ trường phông tỉ lệ 1: 200.000 khu vực Đông Nam thềm lục địa Việt Nam 47 Hình 3.8: Bản đồ trường phông khu vực Cửu Long tỉ lệ 1:200.000 Hình 3.9: Bản đồ trường dư khu vực Cửu Long tỉ lệ 1:200.000 48 Hình 3.10: Bản đồ trường phông khu vực Nam Côn Sơn tỉ lệ 1:200.000 Hình 3.11: Bản đồ trường dư khu vực Nam Côn Sơn tỉ lệ 1:200.000 KẾT LUẬN 49 Qua việc xây dựng đồ Bughe khu vực thềm lục địa Việt Nam vùng Biển Đông kế cận tỉ lệ 1:200.000 để từ xác định phông khu vực số bể trầm tích phạm vi thềm lục địa Việt Nam, rút số nhận xét sau: - So với đồ Bughe có, Bản đồ Bughe xây dựng dựa theo việc tính hiệu ứng trọng lực lăng trụ có tỉ lệ lớn có độ xác cao - Đối với bể trầm tích Cửu Long, phông khu vực tăng dần theo hướng chủ đạo hướng Tây Bắc - Đông Nam, với gradien tăng dần - Đối với bể trầm tích Nam Côn Sơn, phông khu vực tăng theo hướng Tây Tây Nam - Đông Đông Bắc với gradien tăng dần Như vậy, thấy phạm vi thềm lục địa vùng biển Đông kế cận, phông khu vực bậc hai có hướng cấu trúc phù hợp với hướng cấu trúc chủ đạo địa hình bề mặt Moho khu vực thềm lục địa toàn biển Đông (Bùi Công Quế, Nguyễn Hiệp, 1990) [2] TÀI LIỆU THAM KHẢO 50 Tiếng Việt [1] Tôn Tích Ái (2003), Trọng lực thăm dò trọng lực, NXB Đại học Quốc gia Hà Nội, Hà Nội [2] Bùi Công Quế, Nguyễn Hiệp n.n.k (1990), Đặc điểm trường Địa vật lý thềm lục địa Việt Nam vùng kế cận Báo cáo tổng kết đề tài 48-B 03.02, Chương trình nghiên cứu biển 48-B, Hà Nội [3] Mai Thanh Tân (2003), Biển Đông Địa chất-Địa vật lý biển, NXB Đại học Quốc gia Hà Nội, Hà Nội [4] Đỗ Đức Thanh (2006), Các phương pháp phân tích, xử lý tài liệu từ trọng lực, NXB Đại học Quốc gia Hà Nội, Hà Nội, Hà Nội [5] Hoàng Văn Vượng, Đỗ Chiến Thắng (2003), Về khả minh giải tổng hợp tài liệu trọng lực, từ nghiên cứu móng trước Kainozoi TLĐ Việt Nam, Tạp chí Dầu khí số 3/2003, Hà Nội Tiếng Anh [6] Bhaskara Rao, Geophys.J, Roy.Astr.Soc (1986), D Modelling of sedimentary basins from gravity anomalies with variable density contrast,Vol.84,pp.207212 [7] Bhaskara Rao, D., Prakash, M.I., and Ramesh Babu (1990), and 2 D modelling of gravity anomalies with variable density contrast, Geophys Prosp, Vol.38, pp 411-422 [8] Banerje, B., and Gupta, S.P.D (1977), Gravitational attraction of a regtangular parallelepiped Geophysis, Vol.42, pp.1053-1055 [9] Chai, Y and Hinze, W.J., (1988), Gravity inversion of interface above which the density contrast varies exponentially with depth Geophysics, Vol.53, pp 837845 [10] http://topex.ucsd.edu/cgi-bin/get_data.cgi 51 [...]... CHƯƠNG 3 XÁC ĐỊNH PHÔNG KHU VỰC TRÊN MỘT SỐ BỂ TRẦM TÍCH KAINOZOI THỀM LỤC ĐỊA VIỆT NAM 3.1 Xây dựng bản đồ Bughe khu vực thềm lục địa Việt Nam và vùng Biển Đông lân cận 3.1.1 Kết quả tính toán Dựa vào phương pháp xác định dị thường trọng lực trong miền không gian đã trình bày ở chương 2, trong phần này, chúng tôi tiến hành xây dựng bản đồ Bughe khu vực thềm lục địa Việt Nam và vùng Biển Đông kế cận trong... = Z1 (2.10) 2.1.3 Xác định dị thường trọng lực của bể trầm tích Để xác định dị thường trọng lực, bể trầm tích được chia thành các lăng trụ đặt cách nhau ở phía dưới mỗi điểm quan sát, có độ rộng bằng khoảng cách giữa các điểm quan sát và có chiều cao bằng độ sâu của bể trầm tích tại điểm đó Dị 30 thường trọng lực gây ra bởi bể trầm tích được tính bằng cách lấy tổng dị thường trọng lực của tất cả các... không để ra phương pháp giải tích để xác định các thông số của chúng mà thông thường so sánh dạng đường cong của chúng từ bậc nghiêng với các góc α khác nhau 27 CHƯƠNG 2 XÁC ĐỊNH HIỆU ỨNG TRỌNG LỰC THEO PHƯƠNG PHÁP GIẢI TÍCH VÀ PHƯƠNG PHÁP PHỔ 2.1 Xác định dị thường trọng lực theo phương pháp giải tích 2.1.1 Xác định dị thường trọng lực của hình trụ tròn nằm ngang Các vật thể địa chất dạng này là các cấu... Hinze[9] đã thu được công thức tính phổ dị thường trọng lực của lăng trụ này là: ~ (u , v) = 2πfσ (0)(4ab sin 2ua )(sin 2vb)( dg m,n 1 )x λ m , n +2πs x {1-exp[- (λm,n+2πs) Zm,n]}exp [-2πi(uαm+vβn )] (2.15) 2.2.2 Xác định dị thường trọng lực của bể trầm tích Việc tính phổ dị thường trọng lực của bể trầm tích được thực hiện bằng cách lấy tổng phổ dị thường trọng lực của tất cả các lăng trụ này: g~ (u , v)... của lớp trầm tích bề mặt a1 là hệ số suy giảm ở đây R = x 2 + y 2 + z 2 còn ∆x, ∆y tương ứng là khoảng cách giữa các điểm quan sát theo các trục x và y Cuối cùng, dị thường trọng lực của bể trầm tích được xác định theo công thức sau: g ( x, y ) = M N i =1 j =1 ∑ ∑ dg ( x y) (2.12) Trong đó: M là số lăng trụ được chia theo trục x N là số lăng trụ được chia theo trục y 2.2 Xác định dị thường trọng lực theo... trường địa chất dùng để phân tích dị thường trọng lực - Tác dụng của một trong những yếu tố địa chất gây nên dị thường trội hơn Để đảm bảo yêu cầu này, trong thực tế người ta sử dụng các phương pháp biến đổi trường - Các yếu tố địa chất trong vùng tương đối ổn định có thể sử dụng một hoặc tổ hợp phương pháp phân tích chung - Các số liệu đo đạc có độ chi tiết và chính xác cao 11 - Cơ sở lý thuyết phân tích. .. =2πfσ(0)(4absin2ua)(sin2vbxexp[-2πi(uαm+vβn)] (2.16) Dị thường trọng lực của bể trầm tích được xác định nhờ phép biến đổi Fourier ngược theo công thức sau: 32 +∞ +∞ gp(x,y) = ∫ ∫ g~ (u,v) exp [2πi(ux +vy)]dudv (2.17) −∞ −∞ Trên thực tế sự suy giảm mật độ dư của bể trầm tích theo chiều sâu có xu hướng tiệm cận tới giá trị σd không đổi nào đó Bởi vậy sự suy giảm mật độ dư của bể trầm tích theo chiều sâu có thể được xấp xỉ bằng... độ sâu có thể xấp xỉ bằng một hàm bậc hai: σ(z) = a0 + a1z + a2z2 (2.8) Trong đó z là độ sâu; a 1, a2 là các hệ số suy giảm; a 0 là mật độ dư của lớp trầm tích bề mặt Để xác định dị thường trọng lực, bể trầm tích được chia thành các lăng trụ đặt cách nhau ở phía dưới mỗi điểm quan sát, có độ rộng bằng khoảng cách giữa các điểm quan sát và có chiều cao bằng độ sâu của bể trầm tích tại điểm đó 29 Hình... kể trên thu được tại chính điểm quan sát đó Khu vực nghiên cứu được chia thành mạng lưới ô vuông có kích cỡ ∆x = ∆y = 2 km ; Cả khu vực nghiên cứu được chia thành 990x1045 lăng trụ Hiệu ứng trọng lực tại mỗi điểm được tính bằng cách lấy tổng hiệu ứng trọng lực của tất cả các lăng trụ trên Hiệu ứng trọng lực của mỗi lăng trụ được tính theo công thức (2.12) Kết quả tính toán được trình bày lần lượt trên. ..- Vị trí của yếu tố địa chất hoặc vật quặng - Vùng hoặc khu vực cần phải tiến hành nghiên cứu tỉ mỉ hơn - Điểm hoặc vùng nhỏ tại đó có thể đặt được các lỗ khoan hoặc đào hầm lò - Khả năng và điều kiện để phân tích định lượng Trong trường hợp tổng quát, có bốn yếu tố địa chất chính gây nên dị thường trọng lực: - Cấu tạo các lớp trầm tích - Địa hình mặt nền kết tinh - Cấu tạo bên