Luận án tiến sĩ nghiên cứu áp dụng hệ phương pháp phân tích, xử lý hiện đại xác định cấu trúc móng trước kainozoi

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN _ Nguyễn Kim Dũng NGHIÊN CỨU ÁP DỤNG HỆ PHƯƠNG PHÁP PHÂN TÍCH, XỬ LÝ HIỆN ĐẠI XÁC ĐỊNH CẤU TRÚC MÓNG TRƯỚC KAINOZOI TRÊN MỘT SỐ BỂ TRẦM TÍCH THUỘC THỀM LỤC ĐỊA VIỆT NAM THEO TÀI LIỆU TRỌNG LỰC LUẬN ÁN TIẾN SĨ VẬT LÝ Hà Nội - 2018 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN _ Nguyễn Kim Dũng NGHIÊN CỨU ÁP DỤNG HỆ PHƯƠNG PHÁP PHÂN TÍCH, XỬ LÝ HIỆN ĐẠI XÁC ĐỊNH CẤU TRÚC MĨNG TRƯỚC KAINOZOI TRÊN MỘT SỐ BỂ TRẦM TÍCH THUỘC THỀM LỤC ĐỊA VIỆT NAM THEO TÀI LIỆU TRỌNG LỰC Chuyên ngành: Vật lý địa cầu Mã số: 62440111 LUẬN ÁN TIẾN SĨ VẬT LÝ CHỦ TỊCH HỘI ĐỒNG NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Đỗ Đức Thanh TS Hồng Văn Vượng GS.TS Bùi Cơng Quế Hà Nội - 2018 LỜI CAM ĐOAN Tôi xin cam đoan luận án cơng trình nghiên cứu riêng tơi Các kết nêu luận án trung thực chưa công bố cơng trình khác Tác giả luận án Nguyễn Kim Dũng LỜI CẢM ƠN Luận án “ Nghiên cứu áp dụng hệ phương pháp phân tích, xử lý đại xác định cấu trúc móng trước Kainozoi số bể trầm tích thuộc thềm lục địa Việt Nam theo tài liệu Trọng lực ” hoàn thành sở đào tạo Trường Đại học Khoa học Tự nhiên-Đại học Quốc gia Hà Nội, hướng dẫn PGS.TS Đỗ Đức Thanh TS Hoàng Văn Vượng Tác giả bày tỏ lòng biết ơn sâu sắc đến hai Thầy hướng dẫn tận tình bảo hướng dẫn để nghiên cứu sinh hoàn thành luận án Tác giả xin chân thành cảm ơn ban lãnh đạo Viện Địa chất Địa vật lý biển tạo điều kiện cần thiết cho nghiên cứu sinh học tập nghiên cứu Tác giả gửi lời cảm ơn sâu sắc đến đồng nghiệp Viện có ý kiến đóng góp chuyên môn kinh nghiệm cho luận án nghiên cứu sinh Tác giả xin cảm ơn hỗ trợ Thầy Cô Khoa Vật Lý, Bộ mơn Vật lý Địa cầu, cán Phịng Sau đại học giúp nghiên cứu sinh hoàn thiện hồ sơ thủ tục bảo vệ luận án Tác giả xin chân thành cảm ơn đến Thầy, Cô, nhà khoa học công tác viện nghiên cứu, trường Đại học, hội khoa học kỹ thuật địa vật lý Việt Nam, liên đoàn Vật lý địa chất, có đóng góp ý kiến vê chuyên mơn để luận án hồn thiện tốt Nghiên cứu sinh xin dành yêu thương lòng biết ơn tới gia đình người thân ln niềm động viên, khích lệ mạnh mẽ giúp nghiên cứu sinh thực hoàn thành luận án Tác giả luận án Nguyễn Kim Dũng MỤC LỤC Trang Trang phụ bìa Lời cam đoan Lời cảm ơn Mục lục i Danh mục ký hiệu chữ viết tắt iv Danh mục bảng biểu v Danh mục hình vẽ, đồ thị vi MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU VỀ ĐỊA CHẤTĐỊA VẬT LÝ TRÊN THỀM LỤC ĐỊA VIỆT NAM VÀ KẾ CẬN 1.1 Hoạt động điều tra nghiên cứu địa chất địa vật lý biển 1.1.1 Hoạt động điều tra nghiên cứu trước 1975 1.1.2 Hoạt động điều tra nghiên cứu từ 1975 đến 2000 1.1.3 Hoạt động điều tra nghiên cứu từ năm 2000 đến 14 1.2 Tổng quan phương pháp trọng lực nghiên cứu cấu trúc địa chất sâu 22 1.2.1 Tình hình nghiên cứu nước ngồi 22 1.2.2 Tình hình nghiên cứu nước 25 1.3 Kết luận chương 28 i CHƯƠNG 2: NGHIÊN CỨU HỆ PHƯƠNG PHÁP PHÂN TÍCH, XỬ LÝ HIỆN ĐẠI XÁC ĐỊNH CẤU TRÚC BÊN TRONG MÓNG TRƯỚC KAINOZOI 2.1 Phương pháp giải toán ngược 3D xác định phân bố mật độ móng trước Kainozoi 29 2.1.1 Cơ sở lý thuyết 29 2.1.2 Xây dựng chương trình tính tốn 31 2.1.3 Mơ hình số kết tính tốn 33 2.1.4 Nhận xét 36 2.2 Phương pháp xác định biên nguồn gây dị thường giá trị riêng 37 2.2.1 Cơ sở lý thuyết 37 2.2.2 Xây dựng chương trình tính tốn 39 2.2.3 Mơ hình số kết tính tốn 40 2.2.4 Nhận xét 41 2.3 Phương pháp xác định kết hợp vị trí độ sâu đến biên nguồn nhờ đạo hàm tín hiệu giải tích theo hướng giải chập Euler 42 2.3.1 Cơ sở lý thuyết 42 2.3.2 Xây dựng chương trình tính tốn 44 2.3.3 Mơ hình số kết tính tốn 46 2.3.4 Nhận xét 51 2.3.5.Thử nghiệm áp dụng phương pháp xác định vị trí ước tính độ sâu đứt gãy móng trước Kainozoi vùng trũng Sông Hồng 51 2.4 Kết luận chương 56 ii CHƯƠNG 3: KẾT QUẢ ÁP DỤNG HỆ PHƯƠNG PHÁP PHÂN TÍCH, XỬ LÝ HIỆN ĐẠI XÁC ĐỊNH CẤU TRÚC BÊN TRONG MĨNG TRƯỚC KAINOZOI CÁC BỂ TRẦM TÍCH TRÊN THỀM LỤC ĐỊA VIỆT NAM 3.1 Cơ sở số liệu 58 3.1.1 Nguồn số liệu trọng lực 58 3.1.2 Nguồn số liệu trọng lực vệ tinh 61 3.1.3 Nguồn số liệu địa chấn nguồn khác 63 3.2 Xác định cấu trúc bên móng trước Kainozoi bể trầm tích Sơng Hồng.67 3.2.1 Phân bố mật độ đất đá bên móng trước Kainozoi 68 3.2.2 Cấu trúc khối móng trước Kainozoi 72 3.2.3 Hệ thống đứt gãy móng trước Kainozoi 76 3.2.4 Hệ phương pháp xác định cấu trúc bên móng trước Kainozoi 79 3.2.5 Nhận xét chung 82 3.3 Xác định cấu trúc bên móng trước Kainozoi phần thềm lục địa Đông Nam Việt Nam 83 3.3.1 Phân bố mật độ bên móng trước Kainozoi 83 3.3.2 Cấu trúc khối móng trước Kainozoi 89 3.3.3 Hệ thống đứt gãy móng trước Kainozoi 93 3.3.4 Hệ phương pháp xác định cấu trúc bên móng trước Kainozoi 97 3.3.5 Nhận xét chung 100 3.4 Kết luận chương 101 KẾT LUẬN VÀ KIẾN NGHỊ 103 Kết luận 103 Kiến Nghị 104 Danh mục cơng trình khoa học tác giả liên quan đến luận án: 105 Tài liệu tham khảo 106 iii DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT STT Ký hiệu Giải nghĩa 2D Bài toán hai chiều 3D Bài toán ba chiều Vxz Gradient ngang Vzz Gradient thẳng đứng KHCN Khoa học công nghệ ED HGA CGGT BAS 10 GGT Euler Deconvolution (Giải chập Euler) Horizontal gradient amplitude (cường độ gradient ngang) the curvature gravity gradient tensor (độ cong tensor gradient trọng lực) Basement (móng) Gravity Gradient Tensor (tensor gradient trọng lực) iv Ghi DANH MỤC CÁC BẢNG BIỂU STT Tên Bảng Trang Bảng 2.1 Các tham số mơ hình vị trí nguồn 40 Bảng 2.2 Các tham số mơ hình vật thể 46 Bảng 2.3 Các tham số mơ hình vật thể 46 v Ghi DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ STT Tên Hình Hình 1.1 Các bồn trũng trầm tích Đệ tam thuộc Biển Đơng Việt Nam Hình 1.2 Sơ đồ cấu trúc móng Kainozoi bồn dầu khí thềm lục địa Việt Nam Trang 17 18 Hình 1.3: Bản đồ cập nhật cấu trúc móng trước Đệ Tam Biển thềm lục địa Việt Nam (cập nhật đến 19 năm 2016) Hình 2.1 Sơ đồ khối xác định phân bố mật độ móng trước Kainozoi Hình 2.2 Mơ hình ranh giới phân chia thành phần trường tương ứng Hình 2.3 Kết xác định phân bố mật độ đá móng Hình 2.4 Sai lệch dị thường dư với dị thường tính lần lặp cuối tốc độ hội tụ phương pháp 32 34 35 36 Hình 2.5 Sơ đồ khối xác định biên nguồn 39 Hình 2.6 Kết tính tốn mơ hình 41 10 11 12 13 Hình 2.7 Sơ đồ khối xác định vị trí biên ước tính độ sâu biên nguồn Hình 2.8 Kết xác định vị trí biên mơ hình vật thể Hình 2.9 Kết ước tính độ sâu tới nguồn mơ hình vật thể Hình 2.10 Kết xác định vị trí biên mơ hình vật thể vi 45 49 49 50 Ghi TÀI LIỆU THAM KHẢO Tiếng Việt [1] Trần Ngọc Cảnh (2011), Lịch sử ngành dầu khí Việt Nam (đến năm 2010) Tập 1, Nhà xuất trị quốc gia-sự thật, Hà Nội [2] Văn Đức Chương (1991), ” Cấu trúc móng trước Kainozoi vùng trũng Hà Nội”, Tạp chí địa chất, Số 202-203 (1-4)/1991 [3] Nguyễn Thị Dậu, Phan Mỹ Linh, Phan Văn Thắng (2014), “Đánh giá trình vận động thành tạo miocen khu vực phụ trung đông bắc bể Nam Côn Sơn phục vụ lập liệu đầu vào cho mơ hình địa hóa đá mẹ”, Tạp chí dầu khí, số 1, tr.33-43 [4] Vũ Ngọc Diệp, Hồng Ngọc Đang, Trần Mạnh Cường (2011), “Quá trình phát triển thối hóa đá cacbonat tuổi Miocen đới nâng Tri Tơn phần Nam bể trầm tích Sơng Hồng”, Tạp chí dầu khí, số 07, tr.19-27 [5] Vũ Ngọc Diệp nnk (2014), “ Đặc trưng địa chất thành tạo carbonate tuổi Miocen, phần nam bể trầm tích Sông Hồng mối liên quan tới hệ thống dầu khí”, Tạp chí dầu khí, số 01, tr.24-32 [6] Trần Tuấn Dũng nnk (2006), “Kiến tạo đứt gãy theo minh giải tài liệu trọng lực vùng thềm lục địa Đơng Nam Việt Nam”, Tạp chí Khoa học Cơng nghệ biển, T.2, tr.124-133 [7] Trần Tuấn Dũng (2009),“ Ứng dụng phương pháp minh giải tài liệu trọng lực nghiên cứu cấu trúc hệ thống đứt gãy vùng biển Việt Nam lân cận”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập X, tr71-80, NXB khoa học tự nhiên công nghệ, Hà Nội [8] Trần Tuấn Dũng (2010),“ Hệ đứt gãy tây Biển Đơng móng trước Kainozoi khu vực thềm lục địa miền trung theo minh giải tài liệu trọng lực”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập XI, tr.65-72, NXB khoa học tự nhiên công nghệ, Hà Nội [9] Tran Tuan Dung (2012),”Pre-Cenozoic basement structure in the Truong Sa archipelago and sea deep basins”, Tạp chí dầu khí, số 10, tr.17-23 106 [10] Trần Tuấn Dũng, Bùi Công Quế nnk (2013), “Cenozoic basement structure of the south china sea and adjacent areas by modeling and interpreting gravity data”, Russian journal of pacific geology, Vol 7, No.4, pp.227-236 [11] Trần Tuấn Dũng (2014), “Sử dụng phương pháp minh giải kết hợp tài liệu trọng lực, từ địa chấn để xác định phân bố phun trào Bazan núi lửa thềm lục địa Nam trung lân cận”, Tạp chí Dầu khí (8), tr 14-21 [12] Phan Trung Điền nnk (1995), “ Sự phân bố đánh giá triển vọng dầu khí thành tạo địa chất trước Kainozoi thềm lục địa Việt Nam Báo cáo tổng kết đề tài KT-01-17 [13] Phan Trung Dien (1997), “ Pre-Cenozoic basin analysis and petroleum systems on the continental shelf of Vietnam”, Journal of geology, series B, No9-10, pp1-23 http://www.idm.gov.vn/nguon_luc/Xuat_ban/1997/B91.htm [14] Phùng Thị Thu Hằng, Ngơ Gia Thắng, Cao Đình Triều (2014), “Phân khối cấu trúc Tân kiến tạo Kiến tạo đại khu vực Sơng Cả - Rào Nậy”, Tạp chí Địa chất, Loạt A, số 341-345 (3-8/2014) Hà Nội, trang 214-225 [15] Nguyễn Hiệp (2005), Địa chất tài nguyên dầu khí Việt Nam, Tổng cơng ty dầu khí Việt Nam, Hà Nội [16] Nguyễn Thu Huyền, Trịnh Xuân Cường nnk (2016), “Cập nhật đồ cấu trúc móng trước đệ tam sở minh giải tài liệu địa chấn 2D-PVN12”, tạp chí dầu khí, số 11, tr16-20 [17] Phạm Nam Hưng, Lê Văn Dũng (2011), “Cấu trúc địa chất sâu khu vực Hà Nội lân cận sở phân tích tài liệu trọng lực”, Tạp chí Khoa học Trái đất, T.33(2), tr.185 - 190 [18] Lê Chi Mai nnk (2013), “Đặc điểm cấu trúc địa chất bể Nam Côn Sơn sở kết minh giải tài liệu địa chấn”, Tạp chí dầu khí, số 9, tr.16-25 [19] Lê Huy Minh, Lưu Việt Hùng, Cao Đình Triều (2001),”Một vài phương pháp đại phân tích tài liệu từ hàng khơng áp dụng cho vùng Tuần Giáo”, Tạp chí khoa học trái đất T.22(3),tr.207-216, NXB khoa học kỹ thuật, Hà Nội 107 [20] Lê Huy Minh, Lưu Việt Hùng, Cao Đình Triều (2002),“Sử dụng trường véc tơ gradient ngang cực đại việc minh giải tài liệu từ trọng lực Việt nam”, Tạp chí khoa học trái đất, T.24(1), tr.67-80, NXB khoa học kỹ thuật, Hà Nội [21] Lê Huy Minh, Võ Thanh Sơn, nnk (2005),”Phương pháp giải tích tín hiệu giải tích ba chiều ứng dụng minh giải đồ dị thường từ hàng không vùng tuần giáo”, Tuyển tập báo cáo hội nghị KHKT ĐVLVN lần thứ 4, NXB Khoa học kỹ thuật, Hà Nội [22].Trần Nghi (2010), Báo cáo tổng hợp kết khoa học công nghệ đề tài: Nghiên cứu địa tầng phân tập (sequence stratigraphy) bể trầm tích sơng hồng, cửu long, nam sơn nhằm đánh giá tiềm khống sản, Chương trình KHCN cấp nhà nước KC.09-20/06-10 [23] Phùng Văn Phách nnk (2005) “ Bản đồ kiến tạo khu vực Đông Nam Á” [24] Phùng Văn Phách (2009), “ Bàn cấu trúc kiến tạo khu vực biển Nam Trung Bộ”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập X, tr.22-33 [25] Lại Hợp Phòng, Đinh Văn Toàn, Trần Anh Vũ (2011), “Xác định hướng nghiêng số đứt gãy lãnh thổ miền Bắc Việt Nam phương pháp tính cực đại gradient ngang trọng lực”, tạp chí khoa học trái đất, số 33(3ĐB), tr.554-560 [26] Bùi Cơng Quế, Hồng Văn Vượng (1996),“ Nghiên cứu đặc điểm phân bố mật độ đá móng trước Kainozoi khu vực thềm lục địa Đông Nam theo phương pháp mơ hình hóa cấu trúc khối vỏ Trái Đất”, Các cơng trình nghiên cứu địa chất địa vật lý biển,tập II,tr.277-285,Nxb Khoa học kỹ thuật,Hà Nội [27] Bùi Cơng Quế, Hồng Văn Vượng (1998), “Đặc điểm phân bố mật độ bình đồ cấu trúc móng trước Kainozoi khu vực quần đảo Trường Sa”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập III, tr.13-24 [28] Bùi Công Quế (2000), “Nghiên cứu cấu trúc địa chất tiềm khoáng sản vùng biển dải ven biển Việt Nam”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập VI, tr.11-30, NxbKHTNCN Hà Nội 108 [29] Bùi Công Quế, Trần Tuấn Dũng (2008), “Nghiên cứu đặc điểm cấu trúc sâu địa động lực vùng biển Việt Nam kế cận”, Tạp chí khoa học Trái đất, T30, N.4, Tr.481-490 [30] Bùi Công Quế, Trần Tuấn Dũng, Nguyễn Hồng Phương (2016), Nghiên cứu đặc điểm trường địa vật lý cấu trúc sâu vùng biển Việt Nam, Nxb Khoa học tự nhiên công nghệ, Hà Nội [31].Võ Thanh Sơn, Lê Huy Minh,Lưu Việt Hùng (2005),”Xác định vị tri theo phương ngang chiều sâu đến mặt ranh giới mật độ vùng Châu thổ Sông Hồng phương pháp đạo hàm thẳng đứng giải chập Euler tài liệu dị thường trọng lực”, Tạp chí địa chất loạt A, số 287, tr.39-52 [32] Võ Thanh Sơn, Lê Huy Minh nnk (2007),”Xác định vị tri độ sâu ranh giới tương phản từ tính phương pháp tín hiệu giải tích 3D sử dụng đạo hàm bậc cao”, Tuyển tập báo cáo hội nghị khoa học kỹ thuật địa vật lý Việt Nam lần thứ 5, tr.136-147, NXB Khoa học kỹ thuật, Hà Nội [33] Mai Thanh Tân (2003), Chương trình điều tra nghiên cứu biển cấp nhà nước KHCN-06 (1996-2000),Biển Đông III, Địa chất – địa vật lý, Nxb Đại học quốc gia Hà Nội, Hà Nội [34] Nguyễn Thế Tiệp (2006), ATLAS Các điều kiện tự nhiên môi trường vùng biển Việt Nam kế cận, đề tài KC.09-24 thuộc Chương trình Nghiên cứu biển cấp Quốc Gia giai đoan 2000 – 2005 [35] Nguyễn Thế Tiệp (2010), Báo cáo tổng hợp kết khoa học công nghệ đề tài: Nghiên cứu cấu trúc địa chất vùng biển nước sâu (trên 200m nước) nam Việt Nam làm sở khoa học để tìm kiếm tài ngun khống sản liên quan, Chương trình KHCN cấp nhà nước KC.09-18/06-10 [36] Hồng Đình Tiến, Hồng Thi Xn Hương (2012), “Bể Nam Cơn Sơn góc độ địa động lực”, Tạp chí dầu khí, số 8, tr.15-23 [37] Nguyễn Trọng Tín (2010), Báo cáo tổng hợp kết khoa học công nghệ đề tài : Nghiên cứu cấu trúc địa chất đánh giá tiềm dầu khí khu vực 109 Trường Sa Tư – Vũng mây, Chương trình KHCN cấp nhà nước KC.09-25/06-10 [38] Đinh Văn Toàn (1998), “Mối liên quan đứt đoạn trọng lực phân bố đứt gãy kiến tạo lãnh thổ Bắc Việt Nam”, Báo cáo Hội nghị Khoa học nước ASEAN, Hà Nội - tháng 11/1998 [39] Đỗ Đức Thanh (1997), Nghiên cứu áp dụng số phương pháp phân tích xử lý tài liệu dị thường trọng lực thềm lục địa Việt Nam, Luận án Phó tiến sĩ Tốn – Lý, Hà Nội [40] Đỗ Đức Thanh nnk (2013) “ Nghiên cứu áp dụng phương pháp giải toán ngược ba chiều từ trọng lực xác định phân bố mật độ đá móng độ sâu móng từ nhằm phục vụ cơng tác tìm kiếm thăm dị dầu khí thềm lục địa Việt nam”, Đề tài QG 11-04, Đại học QGHN [41] Đỗ Đức Thanh (2003),“Giải toán ngược trọng lực 2D miền tần số xác định bề dày trầm tích Kainozoi trũng Sơng Hồng”, Tạp chí khoa học công nghệ, Trung tâm KHTN CNQG, T.41(4), tr.62-66 [42] Đỗ Đức Thanh (2004), “Giải toán ngược trọng lực 3D miền tần số xác định độ sâu móng trước Kainozoi bể trầm tích Nam Cơn Sơn”, Tạp chí khoa học Trái Đất, T.26(1), tr.86-94 [43] Đỗ Đức Thanh (2004),“ Modelling of the North – Northwest part of Song Hong sedimentary basin from gravity anomalies with exponential density contrast”, Advances in Natural Sciences, Vol.5(3), pp.313-325 [44] Đỗ Đức Thanh (2006), Các phương pháp phân tích, xử lý tài liệu từ trọng lực, Nhà xuất Đại học Quốc gia, Hà Nội [45] Đỗ Đức Thanh, Nguyễn Đình Chiến (2007), “Mơ hình hóa việc giải tốn ngược ba chiều xác định thơng số vật thể độ sâu móng từ”,Tạp chí khoa học công nghệ, tập 45(4), tr.119-124 [46] Đỗ Đức Thanh, Lại Thùy Ngân (2007), “Xây dựng lọc ngược xác định ranh giới phân chia mật độ theo tài liệu dị thường trọng lực”,Tạp chí khoa học công nghệ, tập 45(5), tr.107-115 110 [47] Đỗ Đức Thanh, Giang Kiên Trung (2008),“ Xác định phân bố mật độ đá móng sở kết hợp phương pháp bóc lớp dị thường với việc giải tốn ngược trọng lực”, Tạp chí khoa học cơng nghệ, tập 46(5), tr.121-128 [48] Đỗ Đức Thanh (2011), “Thử nghiệm việc mở rộng toán ngược xác định độ sâu đáy bể trầm tích để xác định độ sâu móng có cấu trúc dạng vịm”, Tạp chí Khoa học Cơng nghệ, T.49 (1), tr.127-134 [49] Cao Đình Triều, Đinh Văn Tồn (1999), “Mơ hình cấu trúc vỏ Trái đất lãnh thổ Việt Nam kế cận sở phân tích tài liệu trọng lực”, Tuyển tập báo cáo khoa học Hội nghị công nghệ biển toàn quốc lần thứ IV, 12-13 tháng 11 năm 1998, Hà nội, tr.854-863 [50] Cao Đình Triều, Phạm Huy Long (2002), Kiến tạo đứt gãy Việt Nam, Nhà xuất Khoa học Kỹ thuật, Hà Nội [51] Cao Đình Triều, Đặng Thanh Hải, Mai Xuân Bách, Ngô Gia Thắng (2003), “Các đới đứt gãy hoạt động phần phía Bắc lãnh thổ Việt Nam”, Tạp chí Địa chất, loạt A, Số 279 (11-12), Hà Nội, trang 8-19 [52] Cao Đình Triều nnk (2005), Báo cáo tổng kết nhiệm vụ: Thành lập đồ đứt gãy biển đông vùng phụ cận tỷ lệ 1:1.000.000 [53] Cao Đình Triều, Phạm Nam Hưng (2008), “ Nâng cao hiệu phương pháp trọng lực nghiên cứu cấu trúc móng trước Kainozoi Việt Nam”, Báo cáo khoa học Hội nghị Dầu khí Việt nam 30 năm phát triển Hội nhập, Tp Hà Nội, trang 557 - 566 [54] Cao Đình Triều nnk (2009), “Một số kết bước đầu nghiên cứu cấu trúc móng trước Kainozoi khu vực Tp Hồ Chí Minh kế cận sở tài liệu trọng lực”, Tạp chí khoa học trái đất, T.31(4), Tr 335-345 [55] Cao Đình Triều, Nguyễn Đức Vinh (2012), “Phân đoạn đứt gãy đánh giá động đất cực đại Việt Nam”, Tạp chí Địa chất, Số 331-332,Hà Nội, tr.59-68 [56] Phan Trọng Trịnh (2012), Kiến tạo trẻ địa động lực đại vùng biển Việt Nam kế cận, Nhà xuất khoa học tự nhiên công nghệ, Hà Nội 111 [57] Nguyễn Như Trung, Nguyễn Thị Thu Hương (2011), “Cấu trúc vỏ Trái đất khu vực Biển Đông Việt Nam”, Tuyển tập báo cáo hội nghị khoa học công nghệ Biển toàn quốc lần thứ V, Quyển 3, Địa lý, địa chất địa vật lý biển, tr43-58, Nxb khoa học tự nhiên công nghệ, Hà Nội [58] Nguyễn Như Trung, Nguyễn Thị Thu Hương (2013), “Xác định nhanh địa hình mặt móng bồn trầm tích theo phương pháp giải ngược trực tiếp 3D số liệu trọng lực”, Tc Các khoa học trái đất, T.26(3), tr258-268 [59] Nguyễn Như Trung, Nguyễn Thị Thu Hương (2013), “Topography of the crust-mantle boundary beneath the East Sea from 3D gravity inversed interpretation”, Acta geophysica, Vol.61(2), pp357-384 [60] Nguyễn Hữu Tuyên, Cao Đình Triều, Phùng Thi Thu Hằng (2012), “Phân khối cấu trúc địa động lực đại vùng Tuần Giáo kế cận”, Tạp chí Địa chất, Số 331-332, Hà Nội, trang 145-154 [61] Hoàng Văn Vượng (1999), “Nghiên cứu yếu tố cấu trúc móng trước Kainozoi thềm lục địa Miền Trung theo tài liệu trọng lực”, Các công trình nghiên cứu địa chất vật lý biển, Tập V, tr143-153 [62] Hoàng Văn Vượng (2000), “Nghiên cứu đặc điểm cấu trúc móng trước Kainozoi mối liên quan đến triển vọng dầu khí thềm lục địa Việt Nam sở xử lý tổng hợp tài liệu trọng lực- từ tài liệu địa vật lý địa chất kiến tạo có khu vực”, Đề tài KHCN cấp Viện Hàn lâm KHCNVN [63] Hoàng Văn Vượng (1999), “ Về khả ứng dụng phương pháp Gradient toàn phần chiều nghiên cứu cấu trúc tầng trầm tích Kainozoi móng trước Kainozoi thềm lục địa Việt Nam”, cơng trình nghiên cứu địa chất địa vật lý Biển, tập V, tr.293.298 [64] Hoàng Văn Vượng (2000), “Nghiên cứu cấu trúc khối-lớp móng trước Kainozoi đới khâu Sông Hồng khu vực vịnh Bắc Bộ theo tài liệu địa vật lý”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập VI, tr.132-142, Nxb khoa học tự nhiên công nghệ Hà Nội 112 [65] Hoàng Văn Vượng (2005), “Biểu ranh giới mật độ theo tín hiệu cực trị GH”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập VI, tr.59-63, NXB khoa học tự nhiên cơng nghệ, Hà Nội [66] Hồng Văn Vượng (2009), “Nghiên cứu cấu trúc sâu bề dày tầng trầm tích khu vực trũng sâu Biển Đơng kế cận theo tài liệu trọng lực”, Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập X, tr.63-70, NXB khoa học tự nhiên công nghệ Hà Nội [67] Hoàng Văn Vượng (2011), “Hiệu phương pháp phân tích tương quan nghiên cứu bề dày trầm tích khu vực trũng sâu Biển Đơng kế cận” Các cơng trình nghiên cứu địa chất địa vật lý Biển, tập XII, tr70-79, NXB khoa học tự nhiên công nghệ, Hà Nội Tiếng Anh [68] Aghajani.H, Ali Moradzadeh and Hualin Zeng (2009), “ Normalized full gradient of gravity anomaly method and its application to the mobrun sulfide body, Canada”, World applied sciences journal, Vol.6(3), pp.393-400 [69] Bagherbandi M (2012), “Moholso: a matlab program to determine crustal thickness by an isostatic and a global gravitational model”, Computers & Geosciences 44 (2012), pp.177-183 [70].Barnes G., and Lumley J (2011),” Proceessing gravity gradient data”, Geophysics, Vol.76 (2), pp.133-147 [71] Barnes G., Barraud J (2012),”Imaging geologic surfaces by inverting gravity gradient data with depth horizons”, Geophysics, Vol 77, No.1 [72].Beiki M (2010), “Analytic signals of gravity gradient tensor and their application to estimate source location”, Geophysics, Vol.75(6), pp.159-174 [73].Beiki M (2011), New techniques for estimation of source parameters, Uppsala Universitet [74].Beiki M., Pedersen L.B (2010), “Eigenvector analysis of gravity gradient tensor to locate geologic bodies”, Geophysics,Vol.75(6), pp.137–149 113 [75].Berezkin.W.M.(1967), “Application of the full vertical gravity gradient to determination to sources causing gravity anomalies (in Russian)”, Expl Geopys.,18,pp.69-76 [76].Bhaskara Rao, D., Prakash, M.I., and Ramesh Babu, N.(1990), “3 and 2 D modelling of gravity anomalies with variable density contrast”, Geophys Prosp, Vol.38, pp 411-422 [77] Bhattacharyya, BK , and Navolio, M.E (1975),“Digital convolution for computing gravity and magnetic anomalies due to arbitrary bodies” Geophysics, Vol.40, pp.981-992 [78] Blakely R.J., and R.W.Simpson (1986), “Approximating edges of source bodies from magnetic or gravity anomalies”, Geophysics,Vol.51(7), pp.1494– 1498, doi: 10.1190/1.1442197 [79] Blakely.R.J.,(1996), Potential theory in gravity and magnetic application, Cambridge University Press [80] Braitenberg C,Susann Wienecke,and Yong Wang (2006),”Basement structures from satellite-derived gravity field: South China Sea ridge”, Journal of geophysical research, Vol 111, B05407, doi:10.1029/2005JB003938 [81].Chai.Y and Hinze, W.J.(1988), “Gravity inversion of interface above which the density contrast varies exponentially with depth”, Geophysics, Vol.53, pp 837-845 [82] Chenot, D and Debeglia, N (1990),“ Three dimensional gravity or magnetic constrained depth inversion with lateral and vertical variation of contrast”, Geophysics,Vol.55, pp.327-335 [83].Chappell A.R, and N J Kusznir (2008), ”Three-dimensional gravity inversion for Moho depth at rifted continental margins incorporating a lithosphere thermal gravity anomaly correction”, Geophys J Int [84] Chakravarthi V., Raghuram H.M., Singh S.B (2002), “3-D forward gravity modeling of basement interfaces above which the density contrast varies continuously with depth”, Computers geosciences, Vol.28, pp.53-57 114 [85] Chakravarthi V., Sundararajan N (2004), “Automatic 3-D gravity modeling of sedimentary basins with density contrast varying parabolically with depth”, Computers geosciences, Vol.30, pp.601-607 [86].Chakravarthi V, Sundararajan N (2008), “TODGINV–A code for optimization of gravity anomalies due to anticlinal and synclinal structures with parabolic density contrast”, Comphuter & Geosciences, Vol 34, pp 955966 [87] Chasseriau P, Chouteau M, (2003), ”3D gravity inversion using a model of parameter covariance”, Journal of Applied Geophysics 52, PP59-74 [88] Cooper.G.R.J., (2008), “Euler Deconvolution with Improved Accuracy and Multiple Different Structural Indices”, Journal of China University of Geosciences, Vol 19(1), p.72–76 [89] Cordell, L and Henderson, R.G (1968),“ Iterative three-dimensional solution of gravity anomaly data using a digital computer” Geophysics, Vol 33, 1968, pp.596- 601 [90] Cordell L, (1979), “Gravi Gravimetric expression of graben faulting in Santa Fe country and the Espanola Basin, New Mexico”, New Mexico Geological Society 30th Annual Fall Field Conference Guidebook, pp.59-64 [91].Cordell L., Grauch V J S (1985), “Mapping basement magnetization zones from aero-magneticdata in the San Juan basin, New Mexico in Hinze W J., Ed The utility of regional gravity and magnetic anomaly maps”, Sot.Explor.Geophys.,tr.181-197 [92].Debeglia N., and J Corpel (1997), “Automatic 3-D interpretation of potential field data using analytic signal derivatives”, Geophysics, Vol.62, pp.87–96, doi:10.1190/1.1444149 [93] Ebrahimzadeh Ardestani, V (2004), “Detection of near – surface anomalies through 2-D normalized full gradient of gravity data”, J.Eath&Space Physics Vol 30, No.2, 2004, P.1-6 115 [94] Epuh.E.E,.et al,(2011),“Basement Depth Estimation of the Gongola Basin Using Second Vertical Derivative Data as Input Anomaly Profile”,European Journal of Scientific Research ISSN1450-216X,Vol.61,No.1(2011),p.172-183 [95] Eshagh M, Romeshkani.M.,(2011), “Generation of vertical–horizontal and horizontal–horizontal gravity gradients using stochastically modified integral estimators”, Advances in Space Research 48 (2011), p.1341–1358 [96] Eshagh M, Abdollahzadeh M., (2012),”Software for generating gravity gradients using a geopotential model based on an irregular semivectorization algorithm”, Computers & Geosciences, 39(2012) , p.152-160 [97] Gosmez-Ortiz D,et al (2005), ”3Dinver.m : a matlab program to invert the gravity anomaly over a 3D horizontal density interface by ParkerOldenburg’s algorithm”,Computers & Geosciences 31, p.513-520 [98].Gottlieb R.G, 1993, Fast gravity, gravity partials, normalized gravity, gravity gradient torque and magnetic field: Derivation, code and data, NASA contractor report 188243 [99] Guo L, Meng X and Lei Shi, (2011), “3D correlation imaging of the vertical gradient of gravity data”, IOP Publishing Journal of Geophysic and Engineering.8 (2011),p.6-12 [100].Guo.C.C.,Xiong Sh.Q.et al., (2014), “Improved Euler method for the interpretation of potential data based on the ratio of the vertical first derivative to analytic signal”, Applied Geophysics,Volume 11(3), p.331–339 [101].Granser, H (1987),“ Nonlinear inversion of gravity using the Schmidt – Lichtenstein approach” Geophysics,Vol.52, pp.88-93 [102] Hidalgo-Gato M.C and Valéria C F Barbosa (2015), ”Edge detection of potential-field sources using scale-space monogenic signal: Fundamental principles”, Geophysics, Vol.80(5), P.J27–J36 [103] Karsli H., Bayrak Y.,(2010), “Application of the normalized total gradient (NTG) method to calculate envelop of seismic reflection signals”, Journal of Applied Geophysics 71, PP.90-97 116 [104] Ma G., Huang D., Liu C (2013), “Application of balanced edge detection filters to estimate the location parameters of the causative sources using potential field data”, Journal of applied geophysics, Vol.99, pp.18-23 [105] Martins.C.M, Valeria C F Barbosa, and João B C SilvaSimultaneous, (2010), “3D depth-to-basement and density-contrast estimates using gravity data and depth control at few points”, Geophysics, Vol75(3), P.I21–I28 [106] Mickus K.L., Juan Homero Hinojosa (2001),”The complete gravity gradient tensor derived from the vertical component of gravity: a Fourier transform technique”, Journal of Applied Geophysics,Vol.46,pp.159-174 [107] Mallick.K (2012),“Bouguer Gravity Regional and Residual Separation: Application to Geology and Environment”, Printed in India Springer [108] Murthy,I.V.R., and Rao, P (1989),“ Two subprograms to calculate gravity anomalies of bodies of finite and infinite strike length with the density contrast differing with depth” Computers & Geosciences.Vol.15(1989), No.8, pp.1265-1277 [109].Nabighian, M N., (1984), “Toward a three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms Fundamental relations”, Geophysics, Vol49, p.780–786 [110] Nagendra R., et al, (1994), “Fortran program based on granser’s algorithm for inverting a gravity field resulting from a density interface”, Computer & Geoscience,Vol.22(3), pp.219-225 [111].Nagendra R.,Prasad.P.V.S, Bhimasankaram.V.L.S.(1996),“ Forward andinverse computer modeling of a gravity field resulting from a density interface using Parker-Oldenberg method”, Computers & Geosciences, Vol22(3),p.227-237 [112] Li Y, Oldenburg D.W, (1998), “3D inversion of gravity data”, Geophysics, Vol.63(1), p.109-119 [113] Li Y, Oldenburg D.W, (1998), “Separation of regional and residual magnetic field data”, Geophysics, Vol.63(2), p.432-439 117 [114] Litinsky, V.A.,1989,”Concept of effective density : key to gravity depth determination for sedimentary basins”, Geophysics, Vol 54, PP.1474-1482 [115] Liu.Yi-Wen, 2012.Hilbert transform and Applications.InTech [116].Oldelbourg, D.W (1974),“ The inversion and interpretation of gravity anomalies”,Geophysics, Vol.39,pp.526-536 [117].Oruc B., Keskinsezer A (2008), ”Structural setting of the Northeastern Biga Peninsula (Turkey) from tilt derivatives of gravity gradient tensors and magnitude of horizontalgravity components”, Pure Appl Geophys, Vol.165, pp.1913–1927 [118].Oruc B., Keskinsezer A.(2008), “Detection of causative bodies by normalized full gradient of aeromagnetic anomalies from east Marmara region, NW Turkey”, Journal of Applied Geophysics 65, pp 39-49 [119].Oruc B, Selim H.H (2011),“Interpretation of magnetic data in the Sinop area of Mid Black Sea, Turkey,using tilt derivative, Euler deconvolution, and discrete wavelet transform”,Journal of Applied Geophysics74,p.194–204 [120] Oruc B, (2012), “Source Location and Depth Estimation Using Normalized Full Gradient of Magnetic Anomalies”, Yerbilimleri, Vol.33 (2), p.177-192 [121].Oruc B.,SertcelikI.,Kafadar O., Selim H.H., (2013), “Structural interpretation of the Erzurum Basin, eastern Turkey, using curvature gravity gradient tensor and gravity inversion of basement relief “, J Appl Geophys, Vol 88, pp.105–113 [122].Pedersen L B., and T M Rasmussen (1990), “The gradient tensor of potential field anomalies: Some implications on data collection and data processing of Geophysics,Vol.55,pp.1558–1566,doi: maps”, 10.1190/1.1442807 [123] Prutkin.I, Vajda P, Tenzer.R, Bielik.M (2011), “3D inversion of gravity data by separation of sources and the method of local correction: Kolarovo gravity high case study”, Journal of Applied Geophysics 75(2011),p.472-478 118 [124] Rajablou R, et al (2012), “Estimation of crustal thickness by combination of two geophysical methods:A case study”,Scientific Research and Essays, Vol 7(24), p.2150-2159, 28 [125].Rao P.Rama., K.V Swamy, I.V.Radhakrishna Murthy (1999),”Inversion of gravity anomalies of three-dimensional density interfaces “, Computer&geosciences, Volume 25(8), PP 887–896 [126] Rapolla A., Cella F., (2002), “Improved techniques in data analysis and interpretation of potential fields: examples of application in volcanic and seismically active areas”, Annals of geophysic, vol.45(6) [127] Ramusen, R and Pedersen, L.B (1979),“ End correction in potential field modelling “ Geophysical prospecting Vol.27, pp.749-760 [128].Reci.H, G.N.Tsokas, C.Papazachos, S.Bushati (2011), ”Conversion of bouguer gravity data to depth, dip, and density contrast with complex attributes analysis technique, in the area of Greece”, Romanian Reports in Physics, Vol 63(1), P 302–320 [129] Reid,A.B, et al, (2013), “Avoidable Euler errors-the use and abuse of Euler deconvolution applied to potential fields”, Geophysical Prospecting [130] Sheng Z, Xiaohong.M (2015), “Improved normalized full-gradient method and its application to the location of source body”, Journal of Applied Geophysics 113, pp.86-91 [131] Shin.Y.H, Choi K.S, Xu.H, (2006), “Three-dimensional forward and inverse models for gravity fields based on the Fast Fourier Transform”, Computers & Geosciences 32, PP 727-738 [132].Talwani,M., Worzel,J.L., Landisman, M., (1959),“Rapid computation of gravitations for two-dimensional bodies with application to the Meandocino submarine Fracture zone”, J Geophys., Res., Vol.64,1959, pp.49-59 [133] Tedla G.E., et al (2011), “A crustal thickness map of Africa derived from a global gravity field model using Euler deconvolution”, Geophysical Journal International, Geophys J Int.187,1-9 119 [134].Troshkov.G.A and Shalaev.S.V.(1961), “Application of the fourier transform to the solution of the reverse problem of gravity and magnetic surveys”, Priki.Geofiz No.30, pp.162-178 [135] Xu.Y, Tianyao Hao, Zhiwei Li, Qiuliang Duan and Lili Zhang, (2009), “ Regional gravity anomaly separation using wavelet transform and spectrum”, Journal of Geophysics and engineering, PP 279-287 [136].Yao L, (2007), “Forward Modeling of Gravity, Gravity Gradients,and Magnetic Anomalies due to Complex Bodies”, Journal of China University of Geosciences, Vol 18(3), p.280–286 [137] Zhang C., M F Mushayandebvu et al (2000), “Euler deconvolution of gravity tensor gradient data”, Geophysics,Vol.65, pp.512–520 [138] Zeng H, Deshu Xu, and Handong Tan (2007), “A model study for estimating optimum upward-continuation height for gravity separation with application to a Bougher gravity anomaly over amineral deposit, Jilin province, northeast, China”,Geophysics, Vol.72(4), P.145-150 [139] Zhdanov M.S., Robert Ellis, et al (2004), “Three-dimensional regularized focusing inversion of gravity gradient tensor component data”, Geophysics, Vol.69(4), p.925-937 [140] Zhou.W., Xiaojuan Du, Jiyan Li, (2013), “The limitation of curvature gravity gradient tensor for edge detection and a method for overcoming it”, Journal of Applied Geophysics 98, p.237-242 [141].William H.Press - Brian P Flannery - Numerical Recipes, Cambridge University Press, 1990 [142] http://topex.ucsd.edu/cgi-bin/get_data.cgi [143] http://www.ngdc.noaa.gov/mgg/mggd.html 120 ... thông tin cho nghiên cứu chi tiết tác giả luận án lựa chọn đề tài luận án tiến sĩ với tiêu đề: “ Nghiên cứu áp dụng hệ phương pháp phân tích, xử lý đại xác định cấu trúc móng trước Kainozoi số... độ xác việc xác định cấu trúc bên móng trước Kainozoi số bể trầm tích thềm lục địa Việt Nam cần thiết 28 CHƯƠNG NGHIÊN CỨU HỆ PHƯƠNG PHÁP PHÂN TÍCH, XỬ LÝ HIỆN ĐẠI XÁC ĐỊNH CẤU TRÚC BÊN TRONG MÓNG... MÓNG TRƯỚC KAINOZOI Nhằm nâng cao hiệu nghiên cứu cấu trúc móng trước Kainozoi theo tài liệu Trọng lực, tác giả luận án đề xuất hệ phương pháp nghiên cứu gồm phương pháp phân tích, xử lý đại tài