Luận án tiến sĩ kỹ thuật xây dựng thuật toán xử lý số liệu GNSS dạng RINEX - Nguyễn Gia Trọng, Trường Đại học Mỏ - Địa chất

Hệ thống vệ tinh nhân tạo liên tục được nâng cấp với nhiều tín hiệu mới. Sự xuất hiện của các chùm vệ tinh thế hệ mới mở rộng khả năng phủ sóng toàn cầu. Độ tin cậy định vị phụ thuộc trực tiếp vào số lượng vệ tinh quan sát được. Việc kết hợp dữ liệu đa hệ thống giúp tăng độ chính xác vị trí. Tốc độ thu nhận số liệu ngày càng nhanh với mật độ dày đặc. Các ứng dụng đòi hỏi độ chính xác đến hàng milimét trong giám sát biến dạng công trình. Định vị vệ tinh không chỉ phục vụ đo vẽ bản đồ mà còn ứng dụng trong giao thông và dẫn đường. Quá trình này đòi hỏi các giải pháp xử lý số liệu GNSS linh hoạt và hiện đại.

Phần lớn các đơn vị đo đạc trong nước vẫn phụ thuộc vào phần mềm thương mại ngoại nhập. Các gói phần mềm này thường đóng kín thuật toán và khó tùy biến. Chi phí đầu tư ban đầu cùng phí duy trì hàng năm tạo gánh nặng lớn. Việc làm chủ thuật toán xử lý số liệu GNSS giúp chủ động trong mọi tình huống kỹ thuật. Các nhà khoa học có thể linh hoạt hiệu chỉnh mô hình theo đặc thù khí quyển từng khu vực. Nhu cầu xây dựng phần mềm nội địa trở nên cấp thiết hơn bao giờ hết. Giải pháp tự phát triển giúp đảm bảo an toàn thông tin và bảo mật dữ liệu địa không gian. Đây là bước đi chiến lược nhằm nâng cao năng lực công nghệ trắc địa bản đồ quốc gia.

Định dạng RINEX là chuẩn quốc tế để lưu trữ và trao đổi dữ liệu quan trắc vệ tinh. Tệp RINEX bao gồm phần đầu chứa thông tin trạm đo và phần thân chứa các trị đo chi tiết. Phiên bản RINEX 2.11 phổ biến với dữ liệu GPS và GLONASS truyền thống. Trong khi đó, định dạng dữ liệu RINEX 3.0x mở rộng hỗ trợ đa hệ vệ tinh và nhiều tần số mới. Cấu trúc tệp được tổ chức khoa học, phân chia rõ ràng giữa tệp trị đo và tệp đạo hàng. Việc đọc và phân tích cú pháp tệp RINEX là bước khởi đầu quan trọng trong chương trình tính toán. Một bộ giải mã RINEX chuẩn xác giúp chương trình tiếp nhận dữ liệu từ mọi loại máy thu trên thị trường.

Tín hiệu vệ tinh bị bẻ cong và trễ pha khi xuyên qua các tầng khí quyển của Trái Đất. Mô hình sai số tầng điện ly và đối lưu (Troposphere & Ionosphere Delay Modeling) là bắt buộc trong tính toán chính xác. Đối với tầng điện ly, tổ hợp không chứa tầng điện ly (Ionosphere-Free) được sử dụng để khử bậc nhất sai số. Mô hình Klobuchar hoặc bản đồ điện ly toàn cầu GIM giúp hiệu chỉnh cho máy thu đơn tần số. Đối với tầng đối lưu, các mô hình Saastamoinen hoặc Hopfield kết hợp hàm chiếu Niell và VMF1 giúp ước lượng độ trễ chính xác. Việc mô hình hóa đúng giúp giảm thiểu sai số vị trí từ vài mét xuống còn vài xentimét.

Hiện tượng trượt chu kỳ xảy ra khi máy thu mất khóa pha tín hiệu vệ tinh trong khoảnh khắc. Nguyên nhân do chướng ngại vật che khuất hoặc tỷ số tín hiệu trên nhiễu giảm đột ngột. Kỹ thuật phát hiện và khắc phục trượt chu kỳ pha (Cycle Slip Detection) là khâu bắt buộc khi xử lý trị đo pha sóng tải. Phương pháp tổ hợp hình học tuyến tính và tổ hợp Melbourne-Wübbena thường được áp dụng hiệu quả. Thuật toán kiểm tra bước nhảy pha trên từng kênh đo và tự động chèn số nguyên chu kỳ bị mất. Việc phục hồi liên tục pha sóng tải đảm bảo tính liền mạch của chuỗi quan trắc. Điều này giúp tối ưu hóa quá trình giải đa trị nguyên trong các bước tiếp theo.

Định vị điểm chính xác cao PPP dựa trên mô hình toán học tổ hợp không điện ly. Phương pháp loại bỏ hoàn toàn ảnh hưởng bậc nhất của tầng điện ly trên cả hai trị đo mã và pha. Véc tơ tham số cần ước lượng bao gồm tọa độ trạm, sai số đồng hồ máy thu, độ trễ đối lưu thiên đỉnh và số nguyên đa trị pha trôi tự do. Do các tham số được ước lượng đồng thời, thời gian hội tụ ban đầu thường kéo dài từ vài chục phút đến vài giờ. Sau giai đoạn hội tụ, độ chính xác tọa độ đạt mức xentimét hoặc milimét. PPP là giải pháp tối ưu cho công tác đo đạc tại các vùng sâu, vùng xa hoặc trên biển.

Bài toán định vị điểm đơn tiêu chuẩn (SPP) và định vị động yêu cầu thuật toán xử lý dữ liệu theo thời gian thực. Bộ lọc Kalman mở rộng (Extended Kalman Filter) là công cụ toán học tối ưu để giải quyết tính chất phi tuyến của phương trình quan sát. Bộ lọc thực hiện chu trình dự báo trạng thái và cập nhật trạng thái liên tục qua từng kỷ nguyên quan sát. Ma trận hiệp phương sai sai số được cập nhật để phản ánh đúng độ tin cậy của các trị đo. Thuật toán lọc Kalman giúp loại bỏ nhiễu ngẫu nhiên, làm mượt quỹ đạo di chuyển và tự động phát hiện các trị đo thô. Đây là nền tảng quan trọng trong các hệ thống dẫn đường và tự hành hiện đại.

Dữ liệu lịch thiên văn phát sóng trong thông điệp vệ tinh thường chứa sai số tọa độ từ một đến hai mét. Để nâng cao độ chính xác, việc tích hợp quỹ đạo chính xác và hiệu chỉnh đồng hồ vệ tinh (IGS Precise Ephemeris & Clock Products) là điều kiện tiên quyết. Dịch vụ GNSS Quốc tế (IGS) cung cấp các sản phẩm quỹ đạo SP3 và tệp đồng hồ vệ tinh với độ chính xác đến vài xentimét và pico-giây. Chương trình xử lý tự động nội suy tọa độ vệ tinh bằng đa thức Lagrange bậc cao. Việc sử dụng sản phẩm IGS chính xác giúp loại bỏ hoàn toàn sai số quỹ đạo và sai số đồng hồ vệ tinh trong mô hình tính toán PPP.

Phương pháp hiệu trị đo là nền tảng của định vị tương đối GNSS. Hiệu đơn giữa hai máy thu triệt tiêu sai số đồng hồ vệ tinh và giảm ảnh hưởng của sai số quỹ đạo. Hiệu kép giữa hai máy thu và hai vệ tinh triệt tiêu hoàn toàn sai số đồng hồ máy thu. Phương trình hiệu kép chỉ còn chứa tọa độ tương đối, độ trễ khí quyển dư và số nguyên đa trị pha sóng tải. Dạng phương trình này giúp đơn giản hóa mô hình toán học và nâng cao tính ổn định của ma trận hệ số. Trị đo pha sóng tải sau khi lập hiệu kép trở thành đại lượng quan trắc chính xác nhất để xác định véc tơ cạnh đáy.

Giải đa trị nguyên pha sóng tải (Carrier Phase Ambiguity Resolution - LAMBDA) là bước quyết định chất lượng định vị tương đối. Đại lượng đa trị nguyên ban đầu có tính chất thực và tương quan cao giữa các tham số. Phương pháp LAMBDA áp dụng phép biến đổi không gian nguyên để làm giảm tương quan ma trận hiệp phương sai sai số. Quá trình tìm kiếm số nguyên được thu hẹp trong một siêu elipsoid tối ưu. Sau khi tìm được bộ nghiệm nguyên, kiểm định tỷ số Ratio Test được thực hiện để đánh giá độ tin cậy. Khi nghiệm nguyên được cố định thành công (Fix), độ chính xác tọa độ tăng vọt từ hàng đề-xi-mét lên hàng mi-li-mét.

Sau khi xử lý từng véc tơ cạnh đáy riêng lẻ, toàn bộ lưới GNSS được đưa vào khối bình sai tổng thể. Phương pháp bình sai bình phương tối thiểu (Least Squares Adjustment) được áp dụng để tìm lời giải tối ưu không thiên vị. Thuật toán thiết lập hệ phương trình số hiệu chỉnh và tìm cực tiểu của hàm mục tiêu tổng bình phương số gia. Ma trận trọng số phản ánh chính xác độ tin cậy của từng cạnh đo thông qua ma trận hiệp phương sai cạnh. Quá trình bình sai lưới tự do và bình sai lưới ràng buộc cung cấp tọa độ chính xác cao và sai số vị trí điểm. Kết quả đánh giá sai số trung phương vị trí điểm đảm bảo độ tin cậy theo quy chuẩn kỹ thuật.

Kiến trúc phần mềm được phân chia thành ba tầng chức năng rõ rệt. Tầng giao diện tiếp nhận dữ liệu đầu vào và thiết lập các tùy chọn thông số kỹ thuật. Tầng xử lý trung tâm chứa các mô đun giải mã RINEX, mô hình hóa sai số khí quyển, lọc Kalman và bình sai bình phương tối thiểu. Tầng xuất dữ liệu cung cấp báo cáo chi tiết về tọa độ, đồ thị sai số và các chỉ số thống kê độ chính xác. Mã nguồn được tối ưu hóa cho các phép toán đại số tuyến tính phức tạp. Kiến trúc mở giúp dễ dàng nâng cấp thêm các thuật toán mới hoặc mở rộng hỗ trợ các hệ vệ tinh tương lai.

Phần mềm được kiểm chứng qua số liệu thực tế tại các trạm quan sát liên tục CORS như Đan Phượng (DANP), Côn Minh (KUNM) và Manila (PIMO). Trong chế độ PPP, tọa độ điểm đạt độ lệch dưới ba xentimét sau thời gian hội tụ hai giờ. Trong chế độ định vị tương đối, các cạnh đáy thử nghiệm có độ lệch chỉ vài milimét so với phần mềm chuẩn quốc tế Bernese 5.0. Kết quả tính toán khẳng định thuật toán giải đa trị nguyên LAMBDA và mô hình sai số hoạt động rất ổn định. Chương trình hoàn toàn đáp ứng các yêu cầu kỹ thuật trong xây dựng lưới khống chế trắc địa và quan trắc dịch chuyển công trình tại Việt Nam.