Tổng quan nghiên cứu

Trong xu thế hiện đại hóa công tác trắc địa và xây dựng bản đồ, công nghệ định vị vệ tinh toàn cầu GNSS đã trở thành công cụ chủ đạo, chiếm tới hơn 85% các dự án thành lập mạng lưới khống chế tọa độ. Tuy nhiên, trên thực tế thi công tại các khu vực đô thị có mật độ nhà cao tầng dày đặc, địa hình rừng núi hiểm trở hoặc thung lũng sâu, tín hiệu vệ tinh thường xuyên bị che khuất hoặc suy giảm nghiêm trọng. Để khắc phục triệt để nhược điểm này, giải pháp tối ưu là tích hợp đồng thời trị đo vệ tinh GNSS với các trị đo mặt đất truyền thống từ máy toàn đạc điện tử như góc bằng, góc đứng và chiều dài cạnh.

Vấn đề cốt lõi đặt ra là làm thế nào để xử lý toán học và bình sai chặt chẽ tập hợp dữ liệu không gian đa chiều hỗn hợp này mà không làm biến dạng mạng lưới hình học. Đề tài tập trung giải quyết bài toán bình sai kết hợp trị đo GNSS và trị đo mặt đất trên hai bề mặt toán học cơ bản: mặt phẳng quy chiếu và mặt Ellipsoid quy ước WGS-84. Mục tiêu cụ thể là xây dựng hoàn chỉnh cơ sở lý thuyết chuyển đổi ma trận hiệp phương sai, thiết lập quy trình tính toán tự động và phát triển chương trình máy tính chuyên dụng nhằm tối ưu hóa độ chính xác vị trí điểm.

Nghiên cứu được triển khai thực nghiệm trên dữ liệu mạng lưới địa chính và trắc địa công trình tại Việt Nam trong giai đoạn 2015 - 2016. Kết quả thực nghiệm chứng minh rằng phương pháp bình sai hỗn hợp giúp nâng cao độ chính xác vị trí điểm từ 15% đến 25%, kiểm soát sai số vị trí điểm ở mức dưới 3 mm đến 5 mm, đồng thời rút ngắn khoảng 30% đến 40% thời gian đo đạc ngoài thực địa so với các phương pháp thành lập lưới độc lập cổ điển.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Cơ sở toán học của nghiên cứu dựa trên nguyên lý bình sai gián tiếp theo phương pháp số bình phương nhỏ nhất, với điều kiện tối ưu là tổng bình phương số hiệu chỉnh nhân với trọng số đạt giá trị cực tiểu. Để tích hợp các nguồn dữ liệu khác nhau về bản chất hình học, nghiên cứu áp dụng lý thuyết chuyển đổi tọa độ không gian 3 chiều và mô hình biến đổi vi phân ma trận hiệp phương sai. Các hệ tọa độ cơ bản được liên kết chặt chẽ bao gồm: hệ tọa độ vuông góc không gian địa tâm (X, Y, Z), hệ tọa độ trắc địa (B, L, H), hệ tọa độ địa diện chân trời (N, E, U) và hệ tọa độ phẳng vuông góc (x, y) theo các phép chiếu Gauss-Kruger hoặc UTM với hệ số biến dạng kinh tuyến trục tương ứng là 1,0000 hoặc 0,9996 cho múi chiếu 6 độ và 0,9999 cho múi 3 độ.

Nghiên cứu làm rõ 4 khái niệm và công cụ then chốt trong xử lý số liệu:

  1. Vector đường đáy baseline GNSS gồm 3 thành phần gia số tọa độ không gian địa tâm đi kèm ma trận hiệp phương sai phản ánh mức độ tương quan sai số nội tại.
  2. Ma trận chuyển đổi vi phân tọa độ bậc một được thiết lập thông qua khai triển chuỗi Taylor để biến đổi ma trận phương sai - hiệp phương sai giữa các không gian hình học.
  3. Sai số trung phương trọng số đơn vị dùng làm tiêu chuẩn đánh giá độ tin cậy và chuẩn hóa trọng số giữa hai nhóm trị đo có bản chất vật lý khác biệt.
  4. Elip sai số vị trí điểm và sai số tương hỗ giữa các cặp mốc khống chế, phản ánh trực quan độ chính xác hình học sau bình sai.

Phương pháp nghiên cứu

Nguồn dữ liệu nghiên cứu sử dụng tập hợp số liệu thực nghiệm đo đạc từ lưới địa chính và lưới đường chuyền cấp IV, bao gồm 8 điểm khống chế mặt đất, hàng chục góc đo ngang, cạnh đo bằng máy toàn đạc điện tử có độ chính xác đo góc 2 giây, đo cạnh 2 mm + 2 ppm, kết hợp với các vector baseline đo bằng máy thu GNSS đa tần số. Phương pháp chọn mẫu tập trung vào các đồ hình lưới hỗn hợp điển hình: lưới đường chuyền chèn cạnh GNSS, lưới tam giác đo góc kết hợp baseline và lưới liên kết hai mảng mặt bằng độc lập.

Lý do lựa chọn phương pháp phân tích gián tiếp tuyến tính hóa là khả năng tự động hóa vượt trội trên máy tính, dễ dàng thiết lập ma trận hệ số và xử lý hệ phương trình chuẩn đối xứng kích thước lớn. Quá trình xử lý dữ liệu được tổ chức theo quy trình 2 bước nghiêm ngặt: bước đầu tiên bình sai độc lập lưới GNSS trong không gian địa tâm để ước lượng phương sai đơn vị tiên nghiệm và lọc sai số thô thông qua sai số khép hình; bước thứ hai đưa các trị đo đã chuẩn hóa vào bình sai tổng thể trên mặt phẳng chiếu hoặc mặt Ellipsoid. Toàn bộ chương trình tính toán được lập trình, kiểm thử và hoàn thiện liên tục trong suốt 11 tháng nghiên cứu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán thực nghiệm và phân tích dữ liệu đã đem lại các kết quả nổi bật sau:

Thứ nhất, việc đưa đầy đủ ma trận hiệp phương sai của các vector baseline GNSS vào quá trình bình sai có ý nghĩa quyết định đến tính trung thực của lưới. Khi sử dụng ma trận trọng số đầy đủ thay vì giả định các trị đo độc lập (ma trận đường chéo), sai số trung phương vị trí điểm giảm từ 12,4 mm xuống còn 4,2 mm, tương ứng mức cải thiện độ chính xác hơn 66%.

Thứ hai, so sánh kết quả bình sai hỗn hợp trên mặt phẳng chiếu và trên mặt Ellipsoid cho thấy độ lệch tọa độ phẳng giữa hai phương pháp là rất nhỏ, dao động trong khoảng từ 0,5 mm đến 2,8 mm đối với khu đo có bán kính dưới 15 km. Cả hai mô hình toán đều đảm bảo tính hội tụ chặt chẽ, tuy nhiên phương pháp bình sai trên Ellipsoid thể hiện ưu thế vượt trội khi khu vực đo có sự thay đổi độ cao địa hình lớn.

Thứ ba, nghiên cứu thực nghiệm đánh giá ảnh hưởng của độ cao trung bình khu đo tại hai mức 200 m và 250 m so với mặt Ellipsoid chuẩn. Kết quả chỉ ra rằng nếu không hiệu chỉnh quy đổi chiều dài cạnh từ thực địa về mặt Ellipsoid, sai số biến dạng cạnh có thể vượt quá 3,5 cm trên mỗi 1.000 m đo đạc, làm suy giảm nghiêm trọng độ chính xác của mạng lưới hạng cao.

Thứ tư, mô-đun phần mềm HHNet được xây dựng thành công đã tự động hóa 100% các khâu: chuyển đổi tọa độ B-L-H sang X-Y-Z và N-E-U, biến đổi vi phân ma trận hiệp phương sai, lập phương trình số hiệu chỉnh và giải hệ phương trình chuẩn kích thước lên tới hàng trăm ẩn số chỉ trong vài giây.

Thảo luận kết quả

Nguyên nhân tạo nên sự vượt trội của phương pháp bình sai hỗn hợp là sự bù trừ hoàn hảo giữa hai công nghệ. Trị đo GNSS cung cấp các liên kết vector không phụ thuộc vào điều kiện thông hướng mặt đất với độ chính xác cạnh dài vượt trội, trong khi trị đo toàn đạc điện tử bổ sung liên kết chi tiết tại các điểm bị che khuất tán cây hoặc công trình. Kết quả này hoàn toàn nhất quán và làm sâu sắc thêm các kết luận của nhiều công trình nghiên cứu trước đây trong và ngoài nước, đồng thời khắc phục được hạn chế của các giải pháp cũ vốn chỉ xét đến chiều dài cạnh mà bỏ qua trị đo góc ngang.

Về mặt biểu diễn trực quan, dữ liệu thực nghiệm được minh họa rõ nét qua các bảng so sánh tọa độ trước và sau bình sai, bảng phân tích số gia tọa độ và biểu đồ elip sai số vị trí điểm. Trên đồ hình biểu diễn, các elip sai số sau khi bình sai hỗn hợp có diện tích thu hẹp rõ rệt, hình dạng tròn đều hơn so với elip sai số thu được từ phương pháp đo toàn đạc đơn thuần, khẳng định sự phân bố độ chính xác đồng đều theo mọi hướng trong không gian mạng lưới.

Đề xuất và khuyến nghị

Nhằm ứng dụng hiệu quả các kết quả nghiên cứu vào thực tiễn sản xuất và công tác đào tạo chuyên ngành, 4 khuyến nghị trọng tâm được đề xuất như sau:

  1. Ban hành quy định chuẩn hóa quy trình xử lý dữ liệu trắc địa hỗn hợp cho các đơn vị đo đạc địa chính và khảo sát công trình. Mục tiêu là bắt buộc áp dụng việc chuyển đổi đầy đủ ma trận hiệp phương sai trong xử lý số liệu GNSS, đảm bảo 100% các đồ án lưới cấp II và cấp III đạt sai số khép tương đối tốt hơn 1/25.000. Thời gian triển khai thực hiện trong vòng 12 tháng bởi Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam.

  2. Hoàn thiện và thương mại hóa mô-đun phần mềm tính toán bình sai tự động HHNet, tích hợp giao diện người dùng thân thiện và tương thích với hệ tọa độ quốc gia VN-2000. Mục tiêu giảm 50% thời gian tính toán nội nghiệp cho kỹ sư trắc địa, hoàn thành việc nâng cấp phần mềm trong thời hạn 6 tháng do nhóm tác giả phối hợp cùng các viện nghiên cứu chuyên ngành.

  3. Triển khai áp dụng mô hình bình sai trên mặt Ellipsoid cho các công trình xây dựng có quy mô lớn, cầu vượt nhịp dài, hầm giao thông và nhà cao tầng tại các khu vực đồi núi có độ cao trên 150 m. Chủ đầu tư và các tổng thầu xây dựng cần quy định chỉ tiêu này trong hồ sơ kỹ thuật nhằm kiểm soát sai số biến dạng cạnh dưới mức 2 mm/km.

  4. Cập nhật bài toán bình sai lưới hỗn hợp 3D và 2D vào chương trình giảng dạy đại học và sau đại học tại các cơ sở đào tạo kỹ thuật trắc địa. Mục tiêu giúp 100% sinh viên và học viên cao học làm chủ thuật toán xử lý dữ liệu hiện đại, triển khai định kỳ từ các năm học tiếp theo dưới sự chủ trì của các khoa chuyên ngành.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị thực tiễn và học thuật cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  • Kỹ sư trắc địa công trình và trắc địa địa chính: Tài liệu cung cấp cẩm nang chi tiết về quy trình thiết kế, đo đạc và xử lý số liệu lưới hỗn hợp trong các điều kiện thực địa phức tạp như đô thị nén, khu công nghiệp hoặc vùng đồi núi hiểm trở, giúp nâng cao độ chính xác định vị lên tới từng milimet.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Trắc địa - Bản đồ và Địa tin học: Luận văn là tài liệu tham khảo chuẩn mực về mặt lý thuyết toán xử lý số liệu trắc địa, cung cấp hệ thống công thức vi phân chuyển đổi tọa độ và ma trận hiệp phương sai được chứng minh mạch lạc.
  • Các lập trình viên và doanh nghiệp phát triển phần mềm Trắc địa - GIS: Nguồn tham khảo giá trị về cấu trúc thuật toán, sơ đồ khối xử lý và các đoạn mã nguồn then chốt trong việc giải hệ phương trình chuẩn kích thước lớn và tự động hóa tính toán trắc địa.
  • Chuyên viên quản lý kỹ thuật tại các cơ quan quản lý đất đai và đo đạc bản đồ: Căn cứ khoa học xác đáng để xây dựng quy chuẩn kỹ thuật, định mức kinh tế kỹ thuật và nghiệm thu các dự án đo đạc thành lập bản đồ địa chính có sử dụng công nghệ phối hợp.

Câu hỏi thường gặp

Tại sao cần phải bình sai kết hợp trị đo GNSS với trị đo toàn đạc điện tử? Trong thực tế đo đạc, công nghệ GNSS bị hạn chế tại nơi bị che khuất tầm nhìn vệ tinh, trong khi toàn đạc điện tử bị giới hạn về tầm nhìn thông hướng giữa các mốc xa nhau. Việc bình sai kết hợp giúp tận dụng ưu thế đo xa chính xác của GNSS và tính linh hoạt vùng khuất của máy toàn đạc, nâng độ chính xác vị trí điểm lên khoảng 20% và tiết kiệm 30% chi phí xây dựng mốc trung gian.

Ma trận hiệp phương sai của vector baseline GNSS có vai trò gì trong bình sai? Ma trận hiệp phương sai chứa đựng thông tin về phương sai sai số và mức độ tương quan phụ thuộc giữa 3 thành phần gia số tọa độ không gian. Nếu bỏ qua các giá trị ngoài đường chéo chính, tọa độ và sai số vị trí điểm sau bình sai có thể bị sai lệch từ 5 mm đến hơn 10 mm, làm sai lệch bản chất hình học của mạng lưới.

Bình sai trên mặt phẳng chiếu và trên mặt Ellipsoid khác nhau như thế nào? Bình sai trên mặt phẳng chiếu đòi hỏi phải chuyển đổi toàn bộ vector baseline và trị đo mặt đất về hệ tọa độ phẳng phẳng trước khi tính toán. Ngược lại, bình sai trên Ellipsoid giữ nguyên bản chất mặt cong của Trái đất, giúp loại bỏ triệt để sai số biến dạng phép chiếu đối với các mạng lưới có độ chênh cao trên 200 m.

Tại sao cần chuẩn hóa trọng số trị đo theo quy trình 2 bước? Do trị đo GNSS và trị đo góc cạnh mặt đất có đơn vị đo và bản chất sai số khác nhau hoàn toàn, quy trình 2 bước cho phép xác định chính xác sai số trung phương trọng số đơn vị của lưới GNSS ở bước 1, làm căn cứ hiệu chỉnh và cân bằng trọng số khách quan giữa hai loại trị đo ở bước 2.

Phần mềm HHNet giải quyết bài toán bình sai hỗn hợp theo nguyên lý nào? HHNet được lập trình dựa trên phương pháp bình sai gián tiếp số bình phương nhỏ nhất. Phần mềm tự động lập hệ phương trình số hiệu chỉnh tuyến tính, thành lập ma trận trọng số khối từ ma trận hiệp phương sai và giải hệ phương trình chuẩn đối xứng để tìm số hiệu chỉnh tọa độ với độ chính xác đạt tới 0,1 mm.

Kết luận

  • Luận văn đã giải quyết toàn diện bài toán xử lý toán học cho mạng lưới trắc địa hỗn hợp, tích hợp thành công dữ liệu không gian GNSS 3D và dữ liệu toàn đạc mặt đất 2D.
  • Thiết lập hoàn chỉnh hệ thống công thức giải tích chuyển đổi tọa độ và ma trận hiệp phương sai giữa 4 hệ tọa độ cơ bản phục vụ bình sai chặt chẽ.
  • So sánh định lượng chi tiết giữa hai phương pháp bình sai trên mặt phẳng chiếu và trên mặt Ellipsoid, khẳng định tính chính xác và phạm vi ứng dụng tối ưu của từng mô hình.
  • Xây dựng thành công phần mềm HHNet, hiện thực hóa quá trình tự động hóa xử lý số liệu lưới địa chính và trắc địa công trình với độ tin cậy đạt 100%.
  • Kiểm chứng thực nghiệm thành công trên dữ liệu thực địa, chứng minh hiệu quả giảm thiểu sai số vị trí điểm xuống mức 3 mm đến 5 mm và tối ưu hóa năng suất lao động.

Đóng góp cốt lõi của công trình là đã khép kín khoảng trống giữa lý thuyết bình sai không gian hiện đại và thực tiễn sản xuất đo đạc tại Việt Nam. Kế hoạch nghiên cứu tiếp theo sẽ tập trung mở rộng mô hình bình sai 4D có tính đến dịch chuyển vỏ Trái đất và tích hợp dữ liệu quét laser 3D trong vòng 24 tháng tới. Hãy liên hệ với tác giả hoặc bộ môn chuyên ngành để tiếp cận phần mềm HHNet và bộ dữ liệu thực nghiệm phục vụ công tác nghiên cứu và triển khai sản xuất.