Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam.BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ ĐỊA CHẤT NGUYỄN THÀNH LÊ NGHIÊN CỨU XÁC ĐỊNH DỊ THƯỜNG ĐỘ CAO TRÊN CƠ SỞ KẾT HỢP CÁC DỮ LIỆU MẶT ĐẤT VÀ VỆ TINH, ÁP DỤNG CHO KHU VỰC MIỀN TRUNG VIỆT NAM NGÀNH.
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT NGUYỄN THÀNH LÊ NGHIÊN CỨU XÁC ĐỊNH DỊ THƯỜNG ĐỘ CAO TRÊN CƠ SỞ KẾT HỢP CÁC DỮ LIỆU MẶT ĐẤT VÀ VỆ TINH, ÁP DỤNG CHO KHU VỰC MIỀN TRUNG VIỆT NAM NGÀNH: KỸ THUẬT TRẮC ĐỊA - BẢN ĐỒ MÃ SỐ: 9.520503 TÓM TẮT LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI - 2023 Công trình hồn thành tại: Bộ mơn Trắc địa cao cấp, Khoa Trắc địa - Bản đồ Quản lý đất đai, Trường Đại học Mỏ - Địa chất, Hà Nội Người hướng dẫn khoa học: PGS TS Nguyễn Văn Sáng, Trường Đại học Mỏ - Địa chất GS TS Bùi Tiến Diệu, Trường Đại học Southeast, Na Uy Phản biện 1: GS TS Võ Chí Mỹ, Hội Trắc địa - Bản đồ - Viễn thám Việt Nam Phản biện 2: TS Vũ Văn Đồng, Cục Bản đồ - Bộ Tổng Tham mưu Phản biện 3: GS TSKH Hoàng Ngọc Hà, Trường Đại học Mỏ - Địa chất Luận án bảo vệ trước Hội đồng đánh giá luận án cấp Trường, họp Trường Đại học Mỏ - Địa chất, vào hồi ngày tháng năm 2023 Có thể tìm hiểu luận án tại: Thư viện Quốc Gia Việt Nam Thư viện Trường Đại học Mỏ - Địa chất MỞ ĐẦU Tính cấp thiết đề tài Nghiên cứu xác định dị thường độ cao, xây dựng mơ hình Geoid/Quasigeoid nhiệm vụ quan trọng quốc gia, khu vực toàn cầu Hiện nay, nhờ có cơng nghệ GNSS, việc xác định độ cao trắc địa cho độ xác cao cao với độ cao chuẩn, cịn nhiệm vụ khó khăn, việc xác định độ cao phải thực trực tiếp bề mặt tự nhiên Trái Đất Theo lý thuyết, mối liên hệ Ellipsoid bề mặt Quasigeoid giá trị dị thường độ cao Bài toán đặt ra, xác định dị thường độ cao với độ xác cao mật độ đủ dày, hồn tồn xác định độ cao chuẩn điểm dựa mối liên hệ độ cao trắc địa mặt Ellipsoid giá trị dị thường độ cao Như vậy, toán xác định độ cao chuẩn so với mơ hình Quasigeoid trở thành toán xác định dị thường độ cao phạm vi tồn lãnh thổ Sóng Quasigeoid lại chia thành thành phần bước sóng (dài, trung, ngắn cực ngắn) Mỗi loại liệu lại phản ánh thành phần bước sóng khác Vì vậy, để xác định xác dị thường độ cao cần phải kết hợp liệu với Cùng với phát triển khoa học cơng nghệ, có nhiều dự án đo cao vệ tinh gradient trọng lực từ vệ tinh Dữ liệu trọng lực (trên mặt đất, tàu, đo cao vệ tinh dự án gradient trọng lực từ vệ tinh) nguồn liệu đầu vào toán xác định dị thường độ cao, cần nghiên cứu sử dụng nguồn liệu trọng lực kết hợp toán xác định dị thường độ cao Việt Nam Đồng thời, liệu đo GNSS-Thủy chuẩn cho phép ta xác định trực tiếp dị thường độ cao với độ xác cao Vì kết xác định dị thường độ cao từ nguồn liệu khác cần đánh giá làm khớp với liệu dị thường độ cao GNSS-Thủy chuẩn 2 Mục tiêu, đối tượng phạm vi nghiên cứu - Xác lập sở khoa học phương pháp luận xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh nói chung thực nghiệm tính tốn thành cơng cho khu vực miền Trung Việt Nam, nói riêng - Lựa chọn sở toán học bao gồm hệ quy chiếu, hệ tọa độ hệ triều thống cho loại liệu thu thập từ nguồn gốc khác nhau, đảm bảo thuận tiện tính tốn đánh giá kết - Xác định dị thường độ cao cho khu vực miền Trung Việt Nam sở kết hợp liệu mặt đất vệ tinh từ nguồn liệu hạn chế Thành lập đồ dị thường độ cao hệ quy chiếu hệ tọa độ WGS84 quốc tế cho khu vực thực nghiệm Nội dung nghiên cứu Nghiên cứu phương pháp xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh theo kỹ thuật Loại bỏ-Tính tốn- Phục hồi phương pháp Collocation bình phương nhỏ từ liệu dị thường trọng lực Xây dựng quy trình module chương trình máy tính phục vụ xác định dị thường độ cao theo kỹ thuật RCR phương pháp LSC Thực nghiệm đánh giá độ xác kết xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh cho khu vực miền Trung Việt Nam hệ quy chiếu hệ tọa độ WGS84 quốc tế Thành lập đồ dị thường độ cao trọng lực cho khu vực thực nghiệm hệ quy chiếu hệ tọa độ WGS84 quốc tế Phương pháp nghiên cứu Phương pháp thống kê, phân tích, thực nghiệm, so sánh, tốn học, ứng dụng tin học phương pháp chuyên gia Ý nghĩa khoa học thực tiễn luận án Ý nghĩa khoa học: Phát triển, hoàn thiện phương pháp quy trình xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh; Góp phần phát triển lý thuyết nghiên cứu hình dạng, kích thước trọng trường Trái Đất Ý nghĩa thực tiễn: Kết nghiên cứu luận án ứng dụng thực tiễn xác định dị thường độ cao từ liệu kết hợp mặt đất vệ tinh phục vụ việc xây dựng Quasigeoid Việt Nam Dị thường độ cao xác định miền Trung sử dụng để tính độ cao chuẩn từ độ cao đo công nghệ GNSS Các luận điểm bảo vệ - Luận điểm thứ nhất: Kết hợp kỹ thuật “loại bỏ–tính tốn–phục hồi” phương pháp collocation bình phương nhỏ giải pháp phù hợp cho phép xác định dị thường độ cao từ nguồn liệu khác khu vực địa hình đặc trưng miền Trung Việt Nam - Luận điểm thứ hai: Quy trình module chương trình tính dị thường độ cao từ kết hợp liệu khác cho phép nâng cao độ xác, mức độ tự động hóa thuận lợi ứng dụng điều kiện thực tế Việt Nam Các điểm luận án - Luận án nghiên cứu thiết lập sở khoa học xác định dị thường độ cao sở kết hợp đa dạng nguồn liệu mặt đất vệ tinh theo kỹ thuật “Loại bỏ – Tính tốn – Phục hồi” phương pháp Collocation bình phương nhỏ cách tường minh mặt lý thuyết, kết thể đắn thực nghiệm áp dụng với đặc điểm địa hình miền Trung Việt Nam - Xây dựng quy trình module chương trình xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh từ nguồn liệu hạn chế miền Trung Việt Nam Quy trình xây dựng chặt chẽ, sát với quy định kỹ thuật, module lập trình có tính tự động hóa, tạo điều kiện thuận lợi cho việc ứng dụng kết nghiên cứu thực tế công tác trắc địa Việt Nam - Thực nghiệm xác định dị thường độ cao trọng lực cho khu vực miền Trung Việt Nam, sở kết hợp liệu mặt đất vệ tinh, độ lệch chuẩn đạt ±0,078 m Thành lập đồ dị thường độ cao trọng lực cho khu vực thực nghiệm hệ quy chiếu hệ tọa độ WGS84 quốc tế, sở cho việc ứng dụng xác định độ cao chuẩn công nghệ đo cao GNSS Cấu trúc nội dung luận án Luận án gồm ba phần: mở đầu, chương nội dung kết luận CHƯƠNG 1: TỔNG QUAN VỀ NGHIÊN CỨU XÁC ĐỊNH DỊ THƯỜNG ĐỘ CAO Ở VIỆT NAM VÀ TRÊN THẾ GIỚI 1.1 Tổng quan xác định dị thường độ cao giới Xác định dị thường độ cao nghiên cứu xây dựng mơ hình Geoid/Quasigeoid nhiệm vụ quan trọng khơng tồn cầu, khu vực mà đặc biệt quan trọng với quốc gia giới Điển hình cho việc nghiên cứu xác định dị thường độ cao xây dựng mơ hình Geoid/Quasigeoid quốc gia giới luận án tìm hiểu xác định dị thường độ cao xây dựng mơ hình Quasigeoid Slovakia, khu vực Châu Âu, Hoa Kỳ, Canada, Ba Lan, Uganda, Australia, New Zealand, Iran xác định mơ hình Geoid hỗn hợp Đài Loan Dựa sở phân tích, đánh giá nghiên cứu xác định dị thường độ cao, xây dựng mơ hình Geoid/Quasigeoid quốc gia khu vực giới, kết đạt thành tựu sau: • Điểm bật tất nước có liệu trọng lực mặt đất dày đặc chi tiết với độ xác cao, mắt lưới thực với kích thước 1′ × 1′ Dữ liệu trọng lực vệ tinh từ dự án CHAMP, GOCE GRACE nhiều quốc gia sử dụng xây dựng mơ hình Geoid/Quasigeoid Dữ liệu trọng lực hàng không liệu trọng lực tàu sử dụng phục vụ cho việc xây dựng mơ hình Geoid/Quasigeoid • Việc xác định dị thường độ cao không thực trực tiếp mà dựa liệu dị thường trọng lực Xác định dị thường độ cao thực theo phương pháp RCR; cơng thức Stokes bình phương nhỏ cải tiến kết hợp số hiệu chỉnh (LSMSA) phương pháp UNB Trường Đại học kỹ thuật New Brunswick, Canada • Các module phần mềm GRAVSOFT nhiều quốc gia sử dụng xác định dị thường độ cao, xây dựng mơ hình Geoid/Quasigeoid • Về liệu tính tốn gồm: liệu trọng lực chi tiết đo mặt đất, biển thuộc phạm vi lãnh thổ, vùng ngồi lãnh thổ sử dụng liệu từ mơ hình trường trọng lực tồn cầu; mơ hình số độ cao, mơ hình số địa hình; mơ hình trọng trường dự án đo trọng lực vệ tinh; mơ hình trọng trường Trái Đất Việc làm khớp đánh giá kết dựa điểm GNSS-Thủy chuẩn • Hiện nay, việc xác định dị thường độ cao xây dựng mơ hình Geoid/Quasigeid với u cầu độ xác ngày cao thực nhiều quốc gia giới Những nước có mơ hình u cầu mơ hình phải có độ xác cao mơ hình cũ Độ xác mơ hình Quasigeoid đạt 10 cm, với số nước đạt cm 1.2 Tổng quan nghiên cứu xác định dị thường độ cao Việt Nam Ở Việt Nam, việc nghiên cứu xác định dị thường độ cao hay xây dựng mơ hình Quasigeoid cho lãnh thổ nhiều tác giả nghiên cứu Tuy nhiên nội dung nghiên cứu chủ yếu thử nghiệm mơ hình trọng trường Trái Đất có cho phù hợp với lãnh thổ Việt Nam cải tiến mơ hình trọng trường Trái Đất có cách hiệu chỉnh giá trị chỉnh vào hệ số mô hình Các cơng trình nghiên cứu cụ thể việc xây dựng mơ hình hồn thiện từ tổng hợp nguồn liệu cịn hạn chế Một vài cơng trình có ứng dụng kết từ sản phẩm dự án đo Gradien trọng lực vệ tinh Các cơng trình nghiên cứu bật Việt Nam kể đến: Xây dựng mơ hình Geoid địa phương lãnh thổ Việt Nam tác giả Lê Minh; Chính xác hóa dị thường độ cao EGM2008 dựa liệu GNSS-Thủy Chuẩn tác giả Nguyễn Duy Đơ Trên sở phân tích, đánh giá cơng trình nghiên cứu nước cơng trình giới xác định dị thường độ cao, xây dựng mơ hình Quasigeoid, nhận thấy vấn đề tồn Việt Nam: • Các nội dung nghiên cứu chủ yếu thử nghiệm mô hình trọng trường Trái Đất có cho phù hợp với lãnh thổ Việt Nam cải tiến mô hình trọng trường Trái Đất có cách hiệu chỉnh giá trị chỉnh vào hệ số mơ hình Các cơng trình nghiên cứu xác định dị thường độ cao, xây dựng mơ hình Quasigeoid cở sở kết hợp đầy đủ loại liệu mặt đất vệ tinh hạn chế • Phương pháp quy trình xác định dị thường độ cao, xây dựng mơ hình Quasigeoid sở kết hợp liệu mặt đất vệ tinh phù hợp với điều kiện liệu có Việt Nam chưa hồn thiện Chưa có phần mềm riêng xác định dị thường độ cao từ liệu kết hợp 1.3 Những vấn đề cần tiếp tục nghiên cứu Từ vấn đề nghiên cứu tổng quan giới Việt Nam, vấn đề cần tập trung nghiên cứu luận án xác định sau: • Nghiên cứu phương pháp kết hợp liệu, bao gồm dị thường trọng lực mặt đất, dị thường trọng lực vệ tinh, dị thường trọng lực từ dự án đo cao vệ tinh, mơ hình số địa hình, mơ hình trọng trường Trái Đất, GNSS-Thủy chuẩn để xác định dị thường độ cao Nghiên cứu tính số hiệu chỉnh như: hiệu chỉnh ảnh hưởng triều, hiệu chỉnh ảnh hưởng địa hình đến liệu trọng lực trước đưa vào tính dị thường độ từ liệu kết hợp Đề xuất việc thống sở toán học hệ quy chiếu, hệ tọa độ hệ triều cho loại liệu thu thập từ nguồn gốc khác tính tốn nội suy dị thường độ cao Cơ sở phương pháp tính chuyển dị thường độ cao hệ tọa độ VN-2000 hệ tọa độ WGS84 quốc tế Nghiên cứu ứng dụng phương pháp LSC, kỹ thuật RCR tính tốn dị thường độ cao sở kết hợp nhiều loại liệu • Nghiên cứu ứng dụng phương pháp LSC, kỹ thuật RCR tính tốn dị thường độ cao sở kết hợp nhiều loại liệu Xây dựng quy trình tính tốn xác định dị thường độ cao từ liệu kết hợp mặt đất liệu vệ tinh • Thực nghiệm tính tốn đánh giá độ xác kết xác định dị thường độ cao cho khu vực miền Trung Việt Nam sở kết hợp liệu mặt đất vệ tinh Thành lập đồ dị thường độ cao hệ quy chiếu hệ tọa độ WGS84 quốc tế, sở cho việc ứng dụng xác định độ cao chuẩn công nghệ đo cao GNSS cho khu vực thực nghiệm Kết luận chương Trong chương 1, luận án nghiên cứu tìm hiểu tình hình xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh Việt Nam giới Từ vấn đề tìm hiểu, nghiên cứu, luận án xác định độ xác đạt mơ hình Geoid/Quasigeoid quốc gia giới, phương pháp kỹ thuật quốc gia sử dụng xác định dị thường độ cao sở kết hợp liệu Đồng thời, xác định trạng tồn Việt Nam, từ đó, xác định nội dung cần thực nghiên cứu luận án CHƯƠNG 2: NGHIÊN CỨU PHƯƠNG PHÁP XÁC ĐỊNH DỊ THƯỜNG ĐỘ CAO TRÊN CƠ SỞ KẾT HỢP CÁC DỮ LIỆU MẶT ĐẤT VÀ VỆ TINH CHO KHU VỰC MIỀN TRUNG VIỆT NAM 2.1 Dị thường độ cao đại lượng liên quan 2.2 Các liệu sử dụng để xác định dị thường độ cao 2.3 Nguyên tắc xử lý liệu từ nguồn khác 2.4 Loại bỏ thành phần bước sóng dài từ mơ hình trọng trường Trái Đất liệu dị thường trọng lực liệu dị thường độ cao 2.5 Loại bỏ ảnh hưởng địa hình liệu dị thường trọng lực liệu dị thường độ cao theo phương pháp RTM 2.6 Tính phần dư dị thường độ cao từ dị thường phần dư trọng lực theo phương pháp Collocation bình phương nhỏ 2.6.1 Cơ sở phương pháp Collocation bình phương nhỏ Như biết, kỹ thuật RCR, nhiễu T chia thành ba phần: T = TGGM + TRTM + TRes (2.65) đó: TGGM ảnh hưởng mơ hình trường trọng lực toàn cầu (Global Geopotential Model - GGM); TRTM ảnh hưởng địa hình từ mơ hình địa hình cịn dư (RTM) [80] TRes trường trọng lực dư Moritz (1980) [45], coi T hàm không gian, đó, mơ hình tuyến tính trị đo thể dạng: (2.66) yi =AiT X + Li ( T ) + ei đó: yi trị đo thứ i; X vectơ tham số; AiT ma trận chuyển vị hàng thứ i ma trận A liên quan đến X; Li số hạng tự nhiễu T tương ứng với trị đo thứ yi, ei sai số trị đo [55] Theo kỹ thuật RCR, mơ hình tuyến tính cơng thức (2.66) biến đổi: δ yi = AiT X + Li ( Tres ) + ei (2.67) đó: δyi trị đo cịn dư Công thức (2.67) công thức phương pháp LSC [35] Biểu diễn dạng ma trận, ta có: δY = AX + L( Tres ) + E (2.68) 12 2.9 Xây dựng quy trình xác định dị thường độ cao sở kết hợp liệu cho khu vực miền Trung Việt Nam Kết luận chương Kết luận chương Trong nội dung chương 2, luận án nghiên cứu phương pháp xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh Việt Nam Các nội dung trình bày chi tiết, sở khoa học chặt chẽ, công thức tường minh Dựa sở phương pháp, luận án xây dựng quy trình tính tốn với cấu trúc quy trình 13 thực chặt chẽ Đây sở cho việc thực nghiệm xác định dị thường độ cao sở kết hợp liệu mặt đất vệ tinh áp dụng cho khu vực miền Trung Việt Nam kỹ thuật RCR, phương pháp LSC Trong nội dung chương, dựa kiến thức lý thuyết tích lũy q trình nghiên cứu, luận án lập trình module chương trình cải tiến có tính tự động hóa “Xác định phương sai mức”, “Xác định dị thường trọng lực dị thường độ cao từ hệ số hàm điều hịa cầu”, “Xác định ảnh hưởng địa hình theo phương pháp RTM”, “Xác định hiệp phương sai thực nghiệm” “Xác định hàm hiệp phương sai lý thuyết” Việc xây dựng sơ đồ quy trình module chương trình cải tiến giúp cho luận án hiểu sâu, nắm sở lý thuyết ứng dụng sở lý thuyết xác định dị thường độ cao kỹ thuật RCR phương pháp LSC Các chương trình kiểm chứng với chương trình gốc Trong nội dung thực nghiệm tính tốn luận án, ngồi module “GEOCOLˮ phần mềm GRAVSOFT luận án phải sử dụng, module chương trình cịn lại thực chương trình nghiên cứu sinh tự lập trình CHƯƠNG 3: THỰC NGHIỆM XÁC ĐỊNH DỊ THƯỜNG ĐỘ CAO CHO KHU VỰC MIỀN TRUNG VIỆT NAM TRÊN CƠ SỞ KẾT HỢP CÁC DỮ LIỆU MẶT ĐẤT VÀ VỆ TINH 3.1 Khu vực thực nghiệm liệu sử dụng 3.1.1 Đặc điểm khu vực thực nghiệm Giới hạn khu vực thực nghiệm có vĩ độ từ 15° ÷ 20°B, kinh độ t 104 ữ 109, kớch thc vựng tớnh ì 5° Khu vực thực nghiệm có địa hình dạng hình tuyến, trải dọc từ Bắc xuống Nam 3.1.2 Dữ liệu sử dụng Dữ liệu gồm: liệu trọng lực, liệu mơ hình trọng trường Trái Đất, liệu mơ hình số địa hình độ phân giải cao liệu điểm song trùng GNSS-Thủy chuẩn thuộc phạm vi khu vực tính Dữ liệu dị thường trọng lực đất liền dị thường trọng lực khơng khí tự cung cấp từ đề tài nghiên cứu khoa học cấp Bộ, Bộ Giáo dục 14 Đào tạo (mã số: B2021-MDA-06) gồm 9.897 điểm, liệu hệ quy chiếu hệ tọa độ WGS84 quốc tế, hệ triều hệ không phụ thuộc triều Phạm vi Biển Đông liệu dị thường trọng lực khơng khí tự trích suất từ mơ hình trường trọng lực toàn cầu độ phân giải cao DTU17GRA [60], gồm 1.034 điểm Mơ hình trọng trường Trái đất mơ hình EIGEN-6C4 Mơ hình số địa hình mơ hình số kết hợp mơ hình số địa hình thuộc phần đất liền SRTM3arc_v4.1 độ phân giải cao 3ʺ × 3ʺ mơ hình số độ sâu GEBCO 2020 độ phân giải 15ʺ × 15ʺ trích suất theo phạm vi thực nghiệm Dữ liệu điểm song trùng GNSS-Thủy chuẩn gồm 131 điểm song trùng thuộc đề tài nghiên cứu khoa học cấp Bộ Các liệu tính quy chiếu thống hệ tọa độ hệ quy chiếu WGS84; hệ triều hệ không phụ thuộc triều (tide free), liệu trọng lực thuộc liệu dị thường trọng lực khơng khí tự để đảm bảo thuận lợi tính tốn đánh giá độ xác kết sau tính Mặt khác, để đánh giá quy trình xây dựng chặt chẽ xác, luận án chia quy trình xây dựng thành ba phương án tính để thực thực nghiệm Mục đích đánh giá kết phương án kết hợp liệu xác định dị thường độ cao Ba phương án thực nghiệm gồm: * Phương án 1: Sử dụng liệu mơ hình EIGEN-6C4 kết hợp liệu GNSS-Thủy chuẩn; * Phương án 2: Sử dụng liệu dị thường trọng lực kết hợp liệu mơ hình EIGEN-6C4 liệu GNSS-Thủy chuẩn; * Phương án 3: Sử dụng liệu dị thường trọng lực kết hợp liệu mơ hình EIGEN-6C4, liệu địa hình liệu GNSS-Thủy chuẩn; 3.2 Xác định dị thường độ cao sở kết hợp liệu mơ hình EIGEN-6C4 liệu GNSS-Thủy chuẩn (PA1) 15 Mục đích xác định dị thường độ cao từ mơ hình EIGEN-6C4 kết hợp với liệu GNSS-Thủy chuẩn Dị thường độ cao khu vực thực nghiệm tính từ mơ hình EIGEN-6C4 làm khớp với liệu dị thường độ cao GNSS-Thủy chuẩn Kết tính tốn đánh giá cách so sánh với liệu GNSS-Thủy chuẩn khơng tham gia tính tốn 3.2.1 Các bước thực PA1 Việc xác định dị thường độ cao đánh giá độ xác thực theo bước: Tính tốn - Sử dụng mơ hình EIGEN-6C4 tính dị thường độ cao cho lưới (3ʹ × 3ʹ) phạm vi khu vực thực ngiệm 101 điểm GNSS-Thủy chuẩn - Xác định độ lệch hệ thống kết tính liệu GNSS-Thủy chuẩn - Làm khớp GNSS-Thủy chuẩn với dị thường độ cao tính từ EIGEN-6C4 - Xác định dị thường độ cao từ dị thường độ cao EIGEN-6C4 GNSS-Thủy chuẩn làm khớp Đánh giá độ xác kết thực nghiệm dựa điểm GNSS-Thủy chuẩn không tham gia tính 3.2.2 Kết thực nghiệm PA1 Độ xác xác định dị thường độ cao sở kết hợp liệu mơ hình EIGEN-6C4 liệu GNSS-Thủy chuẩn theo PA1 đạt ±0,183 m Vậy, xác định dị thường độ cao sở kết hợp liệu mơ hình EIGEN-6C4 liệu GNSS-Thủy chuẩn (PA1) so với việc đánh giá dị thường độ cao trực tiếp từ mơ hình EIGEN-6C4 với dị 16 thường độ cao điểm GNSS-Thủy chuẩn, độ xác nâng lên, từ ±0,194 m lên ±0,183 m 3.3 Xác định dị thường độ cao sở kết hợp liệu trọng lực, liệu mơ hình EIGEN-6C4 liệu GNSS-Thủy chuẩn (PA2) Mục đích nội dung xác định dị thường độ cao từ liệu dị thường trọng lực kết hợp với liệu mơ hình EIGEN-6C4 liệu dị thường độ cao GNSS-Thủy chuẩn Dị thường trọng lực dị thường độ cao GNSS-Thủy chuẩn loại bỏ thành phần bước sóng dài từ mơ hình EIGEN-6C4 Dữ liệu phần dư dị thường trọng lực kết hợp với phần dư liệu dị thường độ cao GNSS-Thủy chuẩn sử dụng để tính phần dư dị thường độ cao trọng lực theo phương pháp LSC Kết tính tốn đánh giá cách so sánh với liệu dị thường độ cao GNSS-Thủy chuẩn khơng tham gia tính 3.3.1 Các bước thực PA2 Trong nội dung 3.3, việc xác định thực theo bước sau: Loại bỏ thành phần bước sóng dài từ mơ hình EIGEN-6C4 cho liệu dị thường trọng lực, dị thường độ cao GNSS-Thủy chuẩn Tính tốn a Tính dị thường độ cao từ liệu dị thường trọng lực loại bỏ thành phần bước sóng dài: 17 - Xác định tham số số hàm hiệp phương sai lý thuyết cho liệu dị thường trọng lực dư - Dựa tham số tìm được, thực tính dị thường độ cao cho tập liệu điểm song trùng sử dụng để tham gia tính dị thường độ cao cho lưới chi tiết - Xác định độ lệch hệ thống kết tính liệu gốc điểm GNSS-Thủy chuẩn tham gia tính - Thực hiệu chỉnh giá trị độ lệch hệ thống cho điểm lưới dị thường độ cao trọng lực b Tính dị thường độ cao từ liệu dị thường độ cao trọng lực hiệu chỉnh độ lệch hệ thống kết hợp liệu GNSS-Thủy chuẩn: - Xác định tham số hàm hiệp phương sai lý thuyết cho liệu dị thường độ cao trọng lực hiệu chỉnh độ lệch hệ thống - Dựa tham số tìm được, thực tính dị thường độ cao cho tập điểm GNSS-Thủy chuẩn sử dụng để đánh giá tính hiệu chỉnh lại cho lưới dị thường độ cao trọng lực sở kết hợp liệu dị thường độ cao hiệu chỉnh độ lệch hệ thống liệu tập dị thường độ cao GNSS-Thủy chuẩn tham gia tính Phục hồi thành phần bước sóng dài từ mơ hình EIGEN-6C4 cho kết dị thường độ cao Đánh giá độ xác 3.3.2 Kết thực nghiệm PA2 Kết xác định dị thường độ cao sở kết hợp liệu trọng lực, liệu mô hình EIGEN-6C4 liệu GNSS-Thủy chuẩn, độ xác đạt là: ±0,092 m Độ xác PA2 cao so với PA1, từ ±0,183 m lên ±0,096 m Như vậy, thực xác định dị thường độ cao sở kết hợp liệu trọng lực, liệu mơ hình EIGEN-6C4 liệu GNSS- 18 Thủy chuẩn cho kết tốt kết hợp liệu mơ hình EIGEN-6C4 với liệu GNSS-Thủy chuẩn 3.4 Xác định dị thường độ cao sở kết hợp liệu trọng lực, liệu mơ hình EIGEN-6C4, liệu địa hình liệu GNSS-Thủy chuẩn (PA3) Mục đích xác định dị thường độ cao từ liệu dị thường trọng lực kết hợp với liệu mơ hình EIGEN-6C4, liệu địa hình liệu dị thường độ cao GNSS-Thủy chuẩn Dị thường trọng lực dị thường độ cao GNSS-Thủy chuẩn loại bỏ thành phần bước sóng dài từ mơ hình EIGEN-6C4, thành phần bước sóng ngắn, bước sóng trung ảnh hưởng địa hình cịn lại liệu phần dư dị thường trọng lực phần dư dị thường độ cao GNSS-Thủy chuẩn Phần dư liệu dị thường trọng lực kết hợp với phần dư liệu dị thường độ cao GNSS-Thủy chuẩn sử dụng để tính phần dư dị thường độ cao trọng lực theo phương pháp LSC Kết tính toán đánh giá cách so sánh với liệu dị thường độ cao GNSS-Thủy chuẩn không tham gia tính Các bước tính tốn thực theo PA3 quy trình luận án xây dựng 3.4.1 Tính ảnh hưởng địa hình theo phương pháp RTM 3.4.2 Thực nghiệm PA3 3.4.2.1 Loại bỏ