Nghiên cứu ứng dụng số liệu trọng lực đo trực tiếp và từ đo cao vệ tinh để hợp nhất kết nối hệ thống đứt gãy trên đất liền và biển khu vực đới bờ đông bắc việt nam

N g U â ttiíé NGHIÊN CỨU ỨNG DỤNG SĨ LIỆU TRỌNG L ự c ĐO TRựC TIÉP VÀ TỪ ĐO CAO VỆ TINH ĐE HỢP NHẤT KÉT NỔI HỆ THỔNG ĐỨT GÃY TRÊN ĐẤT LIỀN VÀ BIỂN KHU V ự c ĐỚI BỜ ĐÔNG BẮC VIỆT NAM TRẦN TUẤN DƯƠNG*1), TRẦN TUẤN DŨNG(1’2), TRẦN TRỌNG LẬP*1), NGUYỄN PHƯƠNG NAM*1), TRẰN HOÀNG TÂM*3), ĐÀO TUẤN LINH*4) (I>Viện Địa chất Địa vật lý biển, Viện Hàn lâm Khoa học Công nghệ Việt Nam (2>H ọc viện Khoa học Công nghệ, Viện Hàn lâm Khoa học Công nghệ Việt Nam (ỉ)Trường Đại học Khoa học Thải Nguyên (4)Trường Đại học Mỏ - Địa chất Tóm tắt: Trên khu vực đới bờ, nói chung, khảo sát đo đạc trọng lực trực tiếp mức độ thưa thớt, chưa đồng đều, nhiều khu vực thực địa hình chuyển tiếp phức tạp Số liệu dị thường trọng lực vệ tinh cỏ độ xác chưa cao so với sổ liệu đo trực tiếp, đặc biệt khu vực ven biển, khu vực đảo, quần đảo khu vực cỏ nhiều bãi ngẩm Do đỏ, dị thường trọng lực vệ tỉnh cần phải có hiệu chỉnh theo sổ liệu đo trực tiếp nhằm nâng cao tính đồng độ xác Trong nghiên cứu này, phương pháp bình phương toi thiểu Collocation áp dụng hiệu chỉnh kết nối sổ liệu trọng lực biển đẩt liền dựa theo sổ liệu trọng lực đo trực tiếp từ đo cao vệ tinh Phương pháp góc nghiêng gradient ngang dị thường trọng lực sử dụng để xác định vị tri đặc điểm cẩu trúc hệ thống đứt gãy khu vực Ket nghiên cứu cho thấy tính thực tiễn, hữu ích tích hợp dị thường trọng lực đo trực tiếp với trọng lực vệ tinh; Nó thể tính hiệu minh giải tài liệu trọng lực xác định, hợp kết hệ thong đứt gãy đất liền biển khu vực đới bờ Đơng Bắc Việt Nam Từ khóa: Trọng lực vệ tinh; trọng lực đo trực tiếp; đới bờ Đông bẳc Việt Nam; phương pháp Collocatỉon Giới thiệu chung Khu vực nghiên cứu bao gồm khu vực đất liền khu vực đới bờ Đông Bắc Việt Nam Khu vực đất liền nằm vùng chuyển tiếp từ địa hình núi cao phía Đơng Bắc xuống vùng đồng thấp phía Tây Nam làm nên riêng biệt cấu trúc địa chất Cho đến việc điều tra, nghiên cứu cấu trúc địa chất cna vùng chuyển tiếp khu vực núi cao đồng bằng, nori mà tương hỗ phức tạp thể Ngày nhận bài: 1/5/2022, ngày chuyển phản biện: 5/5/2022, ngày chấp nhận phàn biện: 9/5/2022, ngày chấp nhận đăng: 28/5/2022 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ s ố 52-6/2022 từ sâu vỏ Trái đất nhiều vấn đề chưa đề cập đển Khu vực đởi bờ Đông Bẳc Việt Nam, trải qua nhiều năm nghiên cửu nhiều nơi chi thực mức độ sơ lược [14] Một điều thuận lợi để khắc phục hạn chế trên, đỏ nguồn số liệu đo cao vệ tinh (cùng với tổ hợp phương pháp địa vật lý) khai thác cách hiệu lấp đầy khoảng trống số liệu mà khảo sát tàu chưa thực Đã có nhiều cơng trình nghiên cửu xác định, nâng cao mức độ xác, độ phân giải dị thường trọng lực vệ tính dị thường trọng lực vệ tính (từ vệ tinh Geosat, Topex/Poseidon, ERS1, Envisat, Jason-l, Ciyosat-2 SARAL/Altika ) dị thường trọng lực đo trực tiếp đất liền, biển để tạo mạng lưới số liệu trọng lực cho đại dương giới [4], [13] Tiêu biểu kể đến Sandwell, D T nnk (1999), Huang, M T., nnk (2006), Zhang, s., (2016) [6], [20] Ở Việt Nam có Bùi Cơng Quế, Trần Tuấn Dũng nnk (2008), Nguyễn Văn Sáng nnk (2012) đẵ thực phép so sảnh tích hợp dị thường trọng lực vệ tinh dị thường đo trực tiếp tàu [3], [10], [15] Các nghiên cứu cho thấy rằng, khu vực gần bờ, trường trọng lực đạt có độ xác khơng cao Đây vấn đề cần phải khắc phục sử dụng số liệu trọng lực khu vực đới bờ Như vậy, nói cịn thiếu cơng trình nghiên cứu đầy đủ hợp kết nối hệ thống đứt gãy chuyển tiếp khu vực đất liền biển Hiện nay, có nhiều phương pháp phát triển để xác định ranh giới nguồn gây H dị thường trọng lục Một số kể đến phương pháp góc nghiêng gradient ngang trọng lực đề xuất Miller Singh, 1994 [8] Sau đó, phương pháp góc nghiêng Wijns nnk (2005) [19], Ferreira nnk (2013) cải tiến thuật toán, nâng cao mức độ chi tiết xác định ranh giới nguồn gây dị thường Nghiên cửu xác định đặc trưng kết nối hệ thống đứt gãy đẩt liền biển góp phần không nhỏ nâng cao hiệu nghiên cứu cấu trúc địa chất số loại hình tai biến liên quan Qua đặc trưng đó, luận giải mối quan hệ mật thiết cấu trúc đứt gãy với số loại hình tai biến địa chất trượt lở đất, xói lở, bồi tụ bờ biển khu vực đới bờ Đông Bắc Việt Nam Nguồn tài liệu sử dụng Trong nghiên cứu này, điểm đo trọng lực trực tiếp sử dụng từ chuyến đo tàu RV Professor Gagarinsky qua khảo sát khoa học hợp tác Việt Nam Liên bang Nga năm 1990, 1992, với nguồn số liệu đo đạc từ dự án điều ừa khảo sát Tổng cục Biển Hải đảo, Tổng cục Địa chất Khoáng sần, Viện Địa chất Địa vật lý biển Ngồi cịn có số lượng tài liệu trọng lực lưu trữ Cục Trắc địa Pháp Các điểm đo trọng lực trực tiếp đất liền 673 điểm biển 404 điểm (trường trọng lực bỉnh thường tính chuyển theo công thức quốc tế hệ Potsdam; nguồn số liệu biểu diễn hệ quy chiếu WGS84 theo độ vĩ độ kinh trắc địa) (hình la) [3], [7], [15], [16], [17] TẠP CHÍ KHOA HQC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Hình a, Minh họa điểm đo trọng lực trực tiếp đẩt liền biển ((!)> Viện Địa chất Khoảng sản; (2)> Cục Trắc địa Pháp; (3)>Tàu Gagarinsky; (4)> Tổng cục Biển Hải đảo; (5)> Viện Địa chất Địa vật lỷ biển)); b, Dị thường trọng lực vệ tinh DỊ thường trọng lực xác định từ đo chênh lệch so với dị thường trọng lực đo trực cao vệ tinh với lưới grid l'x l' (sau tiếp đất liền biển Ở đây, mục đích gọi tắt trọng lực vệ tinh), đưa chinh phương pháp bình phương tối thiểu Sandwell, D T nnk (được cập nhật bổ sung Collocation tích hợp hiệu chỉnh số liệu thời điểm biểu diễn ừọng lực vệ tinh theo số liệu trọng lực đo trực hệ quy chiếu WGS84 theo độ vĩ độ tiếp Phương pháp mơ tả cách kinh trắc địa) (hinh lb) [5], [13] tổng quát sau: Phương pháp nghiên cứu Giả sử khu vực nghiên cứu có k giá trị dị thường trọng lực vệ tinh 3.1 Phương pháp bình phương tối thiểu Collocation A g f , A g f , , A g f m giá trị dị thường Từ số liệu đo cao vệ tình, nhà khoa học tính tốn xác định dị thường trọng lực phạm vi tồn cầu, có vùng biển Việt Nam [13] Tuy nhiên, dị thường trọng lực có độ chỉnh xác chưa cao có K (\ ■) - Ma trận hiệp phương sai dị thường trọng lực (giữa dị thường trọng lực vệ tinh với dị thường trọng lực điểm p K (A g alt,A gP); dị thường trọng lực đo trực tiếp vói dị thường trọng lực điểm p K (A g ăo,A gP); giũa dị thường trọng lực vệ tinh điểm khác K (A g alt,A g ữlt); dị thường trọng lực đo trực tiếp điểm khác ỉ' Trong đó: < ■K(Ag‘*,Ag°l') + CA, K r (Aga*,Ag‘fo) > T Khi đó, theo phương pháp Collocatìon, dị thường trọng lực điểm p tính theo cơng thức [ló], [12]: ẩ 'K(Aga\ A g P) _K{Agđ°,AgP)_ trọng lực đo trực tiếp A gf°, A g ị° , , Ag ị ° K(Ag‘fo,Ạg‘fc) + CA^ _ -1 Ag * _Ag‘fo K (A gđo,A g đo) ; dị thường trọng lực vệ tinh với dị thường trọng lực đo trực tiếp K (A g alt,A g đ0)) KT(-, ■) - Ma trận chuyển vị ma trận hiệp phương sai K(-, •) C aa - Ma trận hiệp phương sai cùa sai số đo (cả đo trực tiếp đo vệ tinh) Các giá trị hiệp phương sai dị thường trọng lực tính theo phương pháp Tscheming cc, Rapp RH, 1974, cơng thức TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐƠ SỐ 52-6/2022 u tính sau [18]: K(Agi,AgJ) = a ỵ d ì (l-ì)2 1=2 ri T j y+1 / d,2 y+1 co _ Ả £ ^ ) P,(cosự)+ Pi (cosự) , l~^N+i - ) ( / + ố) #}.ry v ^ v Trong đó: /](c o s ^ ) - Đa thức Legendre bậc ỉ; iỵ - Khoảng cách cầu điểm i j; n Tị - Khoảng cách đến điểm ì j tỉnh từ gốc tọa độ; R - Bán kính trung bình Trải đất; a - Tham số bổ sung; d i- Phương sai hệ số đến bậc N; b - Hằng số, thường chọn 4; Á - Hằng số có đơn vị (m/s)4; R b - Bán kính hình cầu Bjerhammar Các tham số a, di, N, A, R b xác định cách làm khớp hàm hiệp phương sai lý thuyết với giá trị hiệp phương sai thực nghiệm, chúng xác định theo phần mềm cùa nghiên cứu [11] 3.2 Phương pháp góc nghiêng gradien ngang trọng lực a Xác định góc nghiêng gradỉent ngang trọng lực Phương pháp thực phép chuẩn hóa đạo hàm thẳng đứng gradient ngang trọng lực Góc nghiêng gradient ngang trọng lực (TDR) tỉnh theo công thức: dH[Aff(x,y)] dz TDR {H[A#(;c,y)]} = atan (dn[Ag(x,y)]Ý /dH [A ,gO ,y)]ý dx M dy ) Trong H[Ag(x, y)] gradient ngang trọng lực, tính sau [2]: ịd A g { x , y ) ^ + Ợ A g ( x , y ) j H[Ag(x,y)] dầ.g(x,y) , d&g{x,y) , —^dx— — dy — đao hàm bâc nhât trường dị thường ừọng lực theo phương X, y Ag(x, y) dị thường ừọng lực TDR (H[A^(x,y)]} hàm lượng giác atan nên giá trị biến đổi từ —- đến + - Các điểm cực đại TDR phản ánh ranh giói mật độ đất đá - với góc độ đó, ừong minh giải trọng lực, gọi chúng hệ thống đứt gãy Có thể nói, phương pháp có tính khách quan độ tin cậy xác định hệ thống đứt gãy ranh giới khối cấu trúc nâng hạ [1] b n Xác định điểm cực đại (TDR(H)max) Xác định vị trí giá trị cực đại TDR mạng lưới số liệu tiến hành cách so sánh giá trị TDR(h ) (có tọa độ (ỉ,j)) với điểm xung quanh theo bất đẳng thức [2]: )< J > JD *(Í7„W) TDlAp,^ )< T D ft,j) > T D ấ H ,„ ) TDRựỉ,^ )< T D ấ Ị il ) > TDR(h ,_,m ) )< H M H )> Z « ỉ ( / W , ) TDR TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SÔ 52-6/2022 Nghiên diễn cách chi tiết kết trung gian Sau sử dụng 1077 điểm đo trọng lực trực tiếp để hiệu chỉnh, tích hợp với số liệu dị thường trọng lực vệ tinh, kết so sánh dị thường trọng lực vệ tinh với dị thường họng lực đo trực tiếp trước sau hiệu chỉnh, tích hợp sau: đại TDR Một biến đếm N (0 < N < 4) sử dụng để xác định khả đạt cực đại điểm Cứ bất đẳng thức ứên thỏa mãn thi biến đếm N tăng thêm đơn vị Trong nghiên cứu này, lựa chọn N >2, đó, vị trí giá trị cực đại TD R^H ^) xác định: bd_ ■^max — ~ 2a Trong đó: a, b hệ số phụ thuộc vào TDr {h ì j ); d khoảng cách nút lưới số liệu theo phương tính tốn TDR cực đại điểm Xmax xác định hàm bậc 2: TDR{H)mm = axịw +bxw a+TDR(HiJ) Với bất đẳng thức thỏa mãn, xác định vị trí Xmax tương ứng với giá trị T D R (H )m a x - Kết nghiên cứu 4.1 D ị thường trọng lực hiệu chỉnh hợp Dị thường họng lực trước sau hiệu chỉnh so sánh với Độ lệch chúng sau hiệu chỉnh biểu diễn hình 3a Ở đây, ừong khuôn khổ giới hạn báo tiêu đề báo nêu, khơng sâu vào trình bày, biểu - Độ lệch trung phương dị thường ừọng lực trước hiệu chỉnh với dị thường họng lực đo trực tiểp: aT ầg = ± 8,233 mGal - Độ lệch trung phương dị thường trọng lực sau hiệu chỉnh với dị thường trọng lực đo trực tiếp: aịg = ± 3,373 mGal Có thể thấy sau hiệu chỉnh độ lệch trung phương cải thiện từ ± 8,233 mGal xuống ±3,373 mGal Sai lệch dị thường ữọng lực trước sau áp dụng phương pháp hiệu chỉnh bình phương tối thiểu Collocation từ -4,34 mGal đến +3,40 mGal biểu diễn hình 3a Như vậy, với phương pháp Collocation, dị thường trọng lực vệ tinh hiệu chỉnh, tích hợp với dị thường trọng lực đo trực tiếp Độ chinh xác dị thường trọng lực vệ tinh nâng lên, nâng cao độ xác Kết đạt mạng lưới số liệu trọng lực với kích thước ' X r khu vực đới bờ Đông Bắc Việt Nam lân cận (hình 3b) Dị thường trọng lực sau hiệu chỉnh phương pháp Collocation có giá trị biến đổi khoảng từ -71,5 mGal đến +39,3 mGal; So sánh với nguồn số liệu ban đầu có dải biến đổi tích từ -75,8 mGal đển +40,5 mGal, tương ứng Từ kết so sánh hiệu chỉnh thấy rằng, sau hiệu chỉnh giá trị dị thường trọng lực thay đổi đáng kể xung quanh vị trí điểm đo tàu biển, cách xa giá trị hiệu chỉnh dị thường trọng lực thay đổi Điều nói lên rằng, tốn hiệu chinh dị thường trọng lực, TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐC SỐ 52-6/2022 13 Nghiên cứu mức độ xác nâng cao mật độ điểm đo trực tiếp tàu biển lớn phản bố ưên khu vực nghiên cứu Hình a, Sai lệch dị thườiĩg trọng lực vệ tinh trước sau khỉ hiệu chỉnh, tích họp phương pháp bình phương tối thiểu Collocatỉon; b, DỊ thường trọng lực vệ tinh hiệu chỉnh tích họp (free-air) 4.2 H ệ thông đứt gãy hợp nhât đới phù họp với kêt công trình [9] bờ Đơng Bắc Việt Nam Bằng việc phân tích phân bố khơng gian Trong q trình tính tốn, qua bước trung gian, dị thường trọng lực free-air tính chuyển sang dị thường trọng lực Bughe Trên hình 4a phân bố góc nghiêng građient ngang trọng lực bughe (TDR) cực đại góc nghiêng (TDR(H) max)* Trên phưong diện nghiên cứu đứt gãy theo minh giải tài liệu trọng lực coi phân bố hệ thống đứt gãy, ranh giói mật độ Hệ thống đứt gãy thể vị trí hướng dải cực đại TDR(H)max Mặc dù phân bố cực đại TDR(H)max cịn phức tạp (hình 4a) thể rõ đứt gãy ranh giới mật độ đất đá khu vực Các đứt gãy biểu thị dải cực đại TDR(H)max kéo dài hướng, đặc biệt rõ nét đứt gãy có qui mơ mang tính khu vực Nếu sâu vào nghiên cửu phân bố TDR(H) max cac mức sâu khác (bằng phép giải tích trường lên nửa khơng gian phía lọc trường với bước sóng khác nhau) xác định chi tiết đặc trưng phân bố không gian hệ thống đứt gãy (góc nghiêng, hướng cắm, phương độ sâu phát triển ) Trên phần đất liền, đứt gãy 14 giá trị TDR TDR(H)max với tổ họp phương pháp thích họp, dải điểm cực đại kết nối với tạo nên tranh hệ thống đứt gãy cà vị trí khơng gian quy mơ phát triển Tuy nhiên, phạm vi báo vào mơ tả cách khái qt hình thái cấu trúc số đứt gãy khu vực (hình 4b), là: - Hệ đứt gãy hướng Tây Bắc - Đông Nam thể rõ dải TDR(H)max phát triển liên tục, có giá trị > +0,5 rad, phát triển chù yếu phía Tây Nam khu vực nghiên cửu Trên vùng ven biển Nam Định - Thái Bình - Hải Phịng, đút gãy chuyển hướng chạy song song với đường bờ biển, có phương phát triển Đơng Bắc - Tây Nam Ngoài ra, dọc theo bờ biển Quảng Ninh tồn hệ đứt gãy nhỏ có phương kéo dài song song với đường bờ; đứt bị phân cắt đứt gãy chạy từ đất liền hướng biển - Hệ đứt gãy phương Á Kinh tuyến, chúng phát triển chủ yếu Trung tâm phía Nam khu vực Các đứt gãy có qui mơ nhỏ, phát triển tầng trầm tích nơng TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Nghiên cứu Chúng thường bị cắt bới đứt gãy phương Tây Bắc - Đông Nam Ở khu vực Trung tâm, hầu hết đứt gãy phương Á Kinh tuyến bị không chế đứt gãy có qui mơ lớn (được biểu rõ nét dải TDR(H)max phát triển liên tục, có giá trị > +0,7 rad) phương Tây Bắc - Đông Nam - Phía Bắc khu vực nghiên cứu, vùng núi cao địa hình phức tạp, hệ đứt gãy có phương phát triển đa dạng chạy theo xu địa hình Cụ thề chúng có phương Á Vĩ tuyến, Đơng Bắc - Tây Nam Tây Bắc Đơng Nam Nhìn chung đứt gãy có qui mơ khơng lớn Hình a, Góc nghiêng điểm cực đại góc nghiêng gradỉent ngang dị thưịng trọng lực b, Phân bố hệ thống đứt gãy đất liền biển khu vực Đông Bắc Việt Nam Kết luận Số liệu trọng lực biển đất liền hiệu chỉnh kết nối tăng độ xác Độ lệch trung phương dị thường trọng lực vệ tinh dị thường trọng lực đo trực tiếp cải thiện từ ±8,233 mGal lên đến ± 3,373 mGal; Hệ thống đứt gãy đất liền biển khu vực đới bờ Đông Bắc Việt Nam xác định, họp kết tổ họp phương pháp, mơ hình minh giãi dị thường trọng lực; Kết nghiên cứu cho thấy tính thực tiễn, hữu ích tích họp dị thường trọng lực đo trực tiếp với trọng lực vệ tinh; Nó tính hiệu minh giải tài liệu ưọng lực xác định, họp kết nối hệ thổng đứt gãy đất liền biển khu vực Đông Bắc Việt N am o Lòi cảm on Các tác giả chân thành cảm ơn Viện Địa chất Địa vật lý biển - VAST tạo điều kiện để hoàn thành cơng trình nghiên cửu Tài liệu tham khảo [1] A Eshaghzadeh, A Dehghanpou and R.A Kalantari., 2018 Application of the tilt angle of the balanced total horizontal derivative íilter for the interpretation of potential fi eld data Bollettino di Geoíĩsica Teorica ed Applicata, Vol 59, n.2, pp 161178; https://doi.org/10.4430/bgta0233 [2] Blakely, R, J and Simpson, R.w 1986 Approximating edges of source bodies írommagnetic or gravity anomalies Geophysics, 51, 1494-1498 [3] Bùi Công Quế., 2008 Nghiên cứu xác định đặc trưng cấu trúc địa động lực hệ đứt gãy thềm lục địa Việt Nam Biển Đơng Tuyển tập cơng trình nghiên cứu vật lý địa cầu 2008 Nxb Khoa học tự nhiên Công nghệ Hà Nội Tr 231-245 [4] Emmanuel TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 s Garcia, David T 15 Nghiên cứu Sandwell and Walter H.F Smith Retracking CryoSat-2, Envisat and Jason-l radar altimetry waveforms for improved gravity field recovery Geophysical Joumal International Advance Access published lanuary 2, 2014 doi: 10.1093/gji/ggt469 from double retracked satellite altimetry Joumal o f Geodesy, Volume 84, Issue 3, pp.191-199, DOI 10.1007/s00190-009-03559, Springer [13] Sandwell, D T., Mũller, R Smith, w H F.,Garcia, E., Francis, R., 2014 New global marine gravity model from [5] https://topex.ucsd.edu/ and Jason-l reveals buried tectonic [6] Huang, M T., Zhai, G J., Ouyang, CryoSat-2 Y z., et al., 2006 Recovery of marine gravity structure, Science, Vol 346, no 6205, pp 65íield using integrated data from multi-satellite 67, doi: 10.1126/science 1258213 missions Sci Surv Mapp 31 (6), 37-39 [14] Trần Nghi nnk, 2020 Nghiên c diễn biển dự báo xu biến động địa [7] Liên đoàn Vật lý Địa chất, 2011 Bản đồ trường dị thường trọng lực Việt Nam hệ khu vực ven biển từ cửa Ba Lạt đến cửa (Phần đất liền) tỷ lệ 1:500.000, niên đại 2011 Đáy từ Holocen đến Báo cáo đề tài KC.09.02/16-20 Chương trình nghiên cứu [8] MillerH.G, Singh V., 1994 Potential KH&CN phục vụ quản kỷ biển, hải đảo íield tilt a new concept for location of phát triển kinh tế biển potential fíeld sources Joumal o f Applied [15] Tran Tuan Dung nnk., 2019 Geophysics 32:213-217 Improving accuracy o f altimeter-derived [9] Nguyễn Công Lượng nnk., 1999 marine gravity anomalies in the East Vietnam Bản đồ địa chất khoáng sản Việt nam, tỷ lệ Sea deep-basin and adjacent area Vietnam 1:200.000 Viện Địa chất Khoáng sản Việt Joumal of Marine Science and Technology; Nam Vol 19, No 3B: 43-53 [10] Nguyễn Văn Sáng, 2013 Làm khớp https://doi.org/10.15625/1859dị thường trọng lực tính từ số liệu đo cao vệ 3097/19/3B/14497 tinh với sổ liệu đo trọng lực trực tiếp Tạp chí [16] Trần Tuấn Dũng nnk., 2021 Đ Khoa học kỹ thuật Mỏ-Địa chất, 44/10-2013, tra, đánh giá đặc điểm địa chất-trầm tích Đệ tứ tr.44-48 muộn liên quan đến xói lở bồi tụ khu vực [11] Nielsen, J., Tscheming, C.C., ven biển Thái Bình-Quảng Ninh Đe tài Jansson, T.R.N and Forsberg, R 2012, UQĐTCB.02/21-22 Development and User Testing of a Python [17] Trung tâm Địa chất Khống s Interíace to the GRAVSOFT Gravity Field biển Điều tra đặc điểm địa chất, địa động lực, Programs s Kenyon et al (eds.), Geodesy for địa chất khống sản, địa chất mơi trường Planet Earth, International Association of dự báo tai biến địa chất vùng biển Hải Phòng Geodesy Symposia 136, DOI 10.1007/978-3- Quảng Ninh từ - 30m nước 2007 - 2010 642-20338-1 54, Springer-Verlag Berlin [18] Tscheming, c.c and R H Rap Heidelberg 2012 1974 Closed Covariance Expressions for [12] Ole Baltazar Andersen, Per gravity Anomalies, Geoid Undulations, and Knudsen, Philippa A M Berry, 2010 The Deílections o f the Vertical Implied by DNSC08GRA global marine gravity fíeld 16 TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 Nghiên cứu Anomaly Degree Variance Models Reports of the Department of Geodetic Science No 208, The Ohio State University, Columbus, Ohio [19] Wijns c, Perez c, Kowalczyk p., 2005 Theta map: Edge detection in magnetic data Geophysics 70: 39-43 [20] Zhang, s., Sanđvvell, D.T., Taoyong Jin, Dawei Li, 2016 Retracking o f SARAL/AltiKa radar altimetry waveforms for optimal gravity ííeld recovery Mar Geod http.7/dx.doi.org/10.1080/01490419.2016.126 5032.0 Summary Research and applỉcation of the direct measured gravity and satellite-derived gravity data for consolỉdating and connecting the faults System onland and inshore in the Vỉetnam Northeast coatal zone Tran Tuan Duong, Institute o f Marỉne Geology and Geophvsics, Vỉeừiam Academy o f Science and Technology; Graduate University o f Science and Technology, Vietnam Academy o f Science and Technology Tran Trong Lap, Instỉtute o/Marine Geology and Geophysics, Vỉetnam Academy o f Science and Technology Nguyên Phuong Nam, Institute o f Marine Geoỉogy and Geophysics, Vieừiam Academy o f Science and Technology Tran Hoang Tam, Thai Nguyên Universitv o f Sciences Dao Tuan Linh, Hanoi Universỉty o f Mining and Geology In the Coastal zone, dửect gravity survey is still sparse and uneven; in many areas, it cannot even be conducted due to complicated topography Meanwhile, this area is covered entirely by satellite-derived gravity data However, the satellite-derived gravity anomalies are not highly accurate compared with dừectly measured gravity anomalies, especially in Coastal zones, islands, archipelagos, or areas with many undemeath islands Thereíore, satellite-derived gravity anomalies need to be corrected to improve the accuracy for geological structure research o f these regions, here in Vietnam northeastem Coastal zone and the adjacent areas In this study, the least squares collocation method is applied to coưect and consolidate the gravity data on land and sea based on direct measurement and satellite gravity anomalies The title angle method o f the total horizontal gradient o f the gravity anomaly is used to determine the locations and structural characteristics o f the faults System in the area Research results have shown the practìcality and useM ness o f integrating the satellite-derived gravity anomalies with dừect measured gravity anomalies; It also shows the efficiency in the interpretation o f the gravity data for determining and consolidating the land-sea íaults System in the Vietnam Northeast Coastal zone.O Key words: Satellite-derìved gravity; Direct measured gravity; Vietnam Northeastem Coastal zone; Least square Collocation method TẠP CHÍ KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ SỐ 52-6/2022 17 ... tích họp dị thường trọng lực đo trực tiếp với trọng lực vệ tinh; Nó tính hiệu minh giải tài liệu ưọng lực xác định, họp kết nối hệ thổng đứt gãy đất liền biển khu vực Đông Bắc Việt N am o Lòi cảm... tích hợp với dị thường trọng lực đo trực tiếp Độ chinh xác dị thường trọng lực vệ tinh nâng lên, nâng cao độ xác Kết đạt mạng lưới số liệu trọng lực với kích thước ' X r khu vực đới bờ Đông Bắc Việt. .. chung đứt gãy có qui mơ khơng lớn Hình a, Góc nghiêng điểm cực đại góc nghiêng gradỉent ngang dị thưịng trọng lực b, Phân bố hệ thống đứt gãy đất liền biển khu vực Đông Bắc Việt Nam Kết luận Số liệu