Bối cảnh và vấn đề nghiên cứu
Trong công tác trắc địa hiện đại, việc xác định vị trí không gian của các điểm trên bề mặt Trái Đất (tọa độ mặt bằng $X, Y$ và độ cao $H$) đã có bước chuyển biến căn bản nhờ sự phát triển của công nghệ vệ tinh dẫn đường toàn cầu (GNSS). Tuy nhiên, độ chính xác của kết quả định vị GNSS phụ thuộc chặt chẽ vào nhiều yếu tố sai số nguồn, bao gồm ảnh hưởng của tầng điện ly, tầng đối lưu, mô hình trọng trường Trái Đất và đặc biệt là sai số quỹ đạo vệ tinh.
Đề tài tốt nghiệp "Khảo sát kết quả xử lý số liệu GNSS sử dụng lịch vệ tinh quảng bá và lịch vệ tinh chính xác" do sinh viên Phạm Thanh Bình thực hiện tại Trường Đại học Mỏ - Địa chất (năm 2022) tập trung phân tích vai trò của các dạng thông tin quỹ đạo vệ tinh trong quy trình xử lý số liệu trắc địa vệ tinh độ chính xác cao.
Mục tiêu và nhiệm vụ nghiên cứu của đồ án bao gồm:
- Tổng hợp và hệ thống hóa quy trình kỹ thuật thành lập và xử lý số liệu mạng lưới trắc địa bằng công nghệ GNSS từ khâu thực địa đến văn phòng.
- Trình bày cơ sở toán học của các thuật toán xử lý cạnh (baseline), kiểm tra sai số khép hình và bình sai gián tiếp mạng lưới GNSS không gian 3D.
- Nghiên cứu bản chất vật lý, cơ sở thời gian và các mô hình biểu diễn quỹ đạo vệ tinh nhân tạo theo định luật Kepler.
- Phân tích sự khác biệt về độ chính xác và cấu trúc tham số giữa lịch vệ tinh dự báo (Almanac), lịch vệ tinh quảng bá (Broadcast ephemeris) và lịch vệ tinh chính xác (Precise ephemeris do IGS cung cấp).
Đối tượng nghiên cứu là số liệu đo pha sóng tải và mã giả khoảng cách GNSS, cùng các tập tin lịch vệ tinh tương ứng. Phạm vi nghiên cứu tập trung vào quy trình xử lý số liệu trên phần mềm chuyên dụng Trimble Business Center (TBC) phiên bản 2.0 và các quy chuẩn kỹ thuật quốc gia hiện hành về đo đạc lưới trắc địa.
Cơ sở lý thuyết và phương pháp
Khung lý thuyết của đồ án dựa trên các nguyên lý cơ bản của trắc địa vệ tinh, thiên văn học và giải tích sai số:
- Cơ học thiên thể: Mô hình chuyển động Kepler 2 vật thể trong trường trọng lực trung tâm, xác định qua 6 tham số quỹ đạo cơ bản $(\Omega, i, \omega, a, e, T_0)$.
- Lý thuyết hệ thời gian: Khung quan hệ giữa thời gian quốc tế (UT1), thời gian phối hợp quốc tế (UTC), thời gian nguyên tử quốc tế (IAT) và giờ hệ thống GPS (GPST).
- Lý thuyết bình sai trắc địa: Nguyên lý số bình phương nhỏ nhất Gauss-Markov ($V^T P V = \min$) áp dụng cho bài toán bình sai gián tiếp lưới không gian 3 chiều có trọng số tương quan.
- Lý thuyết hệ quy chiếu và tính chuyển tọa độ: Mô hình tính chuyển 7 tham số (3 tham số tịnh tiến gốc $dX, dY, dZ$; 3 góc xoay $\epsilon_x, \epsilon_y, \epsilon_z$; 1 hệ số tỷ lệ $m$) giữa hệ tọa độ không gian địa tâm toàn cầu WGS-84 và hệ tọa độ trắc địa địa phương.
Phương pháp nghiên cứu được tác giả sử dụng là phương pháp phân tích lý thuyết kết hợp hệ thống hóa quy trình công nghệ và thực nghiệm xử lý dữ liệu:
- Phương pháp phân tích - tổng hợp: Phân tích cấu trúc các loại tệp tin dữ liệu vệ tinh tiêu chuẩn (RINEX, DAT, YUMA), cấu trúc thông số dẫn đường và quy chuẩn kỹ thuật lưới hạng II, hạng III quốc gia.
- Phương pháp mô hình hóa toán học: Thiết lập hệ phương trình số hiệu chỉnh sai phân pha bậc hai, bậc ba trong giải cạnh và hệ phương trình chuẩn bình sai mạng lưới.
- Nguồn dữ liệu: Tài liệu quy chuẩn kỹ thuật trắc địa Việt Nam, dữ liệu quỹ đạo của Cục Đo đạc bờ biển Hoa Kỳ, Trung tâm Dịch vụ GNSS Quốc tế (IGS), và các cấu hình mô phỏng xử lý trên phần mềm Trimble Business Center 2.0.
Nội dung chính theo từng chương
Chương 1: Quy trình xử lý số liệu lưới GNSS
Chương 1 trình bày một cách tuần tự các bước công nghệ từ thiết kế, thu nhận tín hiệu ngoài thực địa đến các thuật toán tính toán trong phòng nội nghiệp.
Quy trình ngoài thực địa bắt đầu từ công tác thiết kế mạng lưới, chọn điểm, chôn mốc, lập lịch đo và thiết kế ca đo. Khi đo tĩnh, các yêu cầu kỹ thuật nghiêm ngặt gồm: định tâm, cân bằng máy, đo độ cao ăng-ten chính xác đến 1 mm (phân biệt độ cao nghiêng Uncorrected Vertical và độ cao thẳng đứng True Vertical), ghi chép sổ đo và bảo đảm khoảng thời gian đo đồng thời giữa các trạm máy. Tác giả nhấn mạnh nguyên tắc thời gian chung ($\Delta t_{1-2} = T_3 - T_2$) và số lượng vệ tinh quan sát đồng thời giữa các trạm để bảo đảm chất lượng giải cạnh baseline.
Công tác xử lý nội nghiệp được chia thành các giai đoạn cụ thể:
- Trút số liệu: Sử dụng các phần mềm chuyên dụng (như Data Transfer, GPload của Trimble hoặc Topcon Link) để đưa dữ liệu trị đo pha ($L_1, L_2$), khoảng cách giả ($C_1, P_1, P_2$), Doppler ($D_1, D_2$) và lịch vệ tinh từ máy thu vào máy tính qua cổng USB hoặc cáp COM.
- Xử lý vectơ cạnh (Baseline processing): Giải cạnh qua hai bước (tham số mặc định và chế độ can thiệp như lọc bỏ vệ tinh xấu, tăng góc ngưỡng trên $15^\circ$, cắt bỏ thời gian nhiễu). Đánh giá chất lượng cạnh qua phương sai chuẩn (Variance - giá trị tối ưu từ 1 đến 8 hoặc 9 đối với máy 1 tần số), tỷ số Ratio (phải đạt từ 2 đến 3 trở lên để đạt lời giải cố định Fixed, tránh dùng lời giải động Float).
- Kiểm tra sai số khép hình: Tính toán sai số khép các thành phần $f_X, f_Y, f_Z$, sai số khép toàn phần $f_{(X,Y,Z)} = \sqrt{f_X^2 + f_Y^2 + f_Z^2}$ và sai số khép tương đối $1/T = f_{(X,Y,Z)} / \sum D$. Văn bản tổng hợp các hạn sai theo quy chuẩn quốc gia cho lưới hạng II và hạng III:
| Tổng chiều dài vòng khép ($\sum D$) |
Lưới Hạng II |
Lưới Hạng III |
| $< 5\text{ km}$ |
— |
$f_s < 5\text{ cm}$ |
| $5\text{ km} - 10\text{ km}$ |
— |
$1/100.000$ |
| $10\text{ km} - 25\text{ km}$ |
$1/300.000$ |
$1/150.000$ |
| $25\text{ km} - 50\text{ km}$ |
$1/500.000$ |
$1/300.000$ |
| $> 50\text{ km}$ |
$1/1.000.000$ |
$1/500.000$ |
- Thuật toán bình sai gián tiếp lưới 3D: Thiết lập hệ phương trình số hiệu chỉnh cho từng cạnh đo $i-j$:
$$V = A \Delta X + L$$
Với ma trận trọng số $P = M^{-1}$ lập từ các ma trận hiệp phương sai $3 \times 3$ của từng baseline. Bình sai được thực hiện theo nguyên lý $V^T P V = \min$ qua 2 giai đoạn: bình sai tự do (1 điểm gốc tối thiểu) để kiểm tra mô hình trọng số thông qua hệ số phương sai đơn vị Reference Factor và kiểm định Chi-Square; sau đó bình sai phụ thuộc với toàn bộ điểm gốc liên kết mô hình Geoid (*.ggf).
- Tính chuyển tọa độ: Trình bày mô hình tính đổi từ tọa độ địa tâm $(X,Y,Z)$ sang tọa độ trắc địa $(B,L,H)$, tính chuyển sang tọa độ phẳng $(x,y)$ theo phép chiếu hình trụ ngang Transverse Mercator (phép chiếu UTM với hệ số biến dạng $m_0=0{,}9996$ cho múi $6^\circ$, $m_0=0{,}9999$ cho múi $3^\circ$; phép chiếu Gauss-Krüger với $m_0=1$) và xác định 7 tham số chuyển đổi tọa độ giữa hai hệ quy chiếu bằng phương pháp số bình phương nhỏ nhất với số điểm song trùng $n \ge 3$.
Chương 2: Khái quát và phân tích lịch vệ tinh trong công nghệ GNSS
Chương 2 đi sâu vào bản chất khoa học của lịch vệ tinh, các hệ thời gian tham chiếu và sự tác động của sai số quỹ đạo đến độ chính xác vị trí điểm.
Lịch vệ tinh là tập hợp số liệu mô tả vị trí của vệ tinh trong không gian theo thời gian. Đồ án phân loại 3 dạng lịch vệ tinh chính với mức độ chính xác tương ứng:
| Loại lịch vệ tinh |
Độ chính xác sai số |
Đặc điểm và nguồn gốc |
| Niên giám (Almanac) |
Vài km |
Dùng để định vị gần đúng, tìm vệ tinh; cập nhật tối thiểu 6 ngày một lần. |
| Quảng bá (có SA) |
$2 - 50\text{ m}$ |
Khi Bộ Quốc phòng Mỹ kích hoạt chế độ làm giảm độ chính xác chủ định (SA). |
| Quảng bá (không có SA) |
$1 - 5\text{ m}$ |
Phát trực tiếp từ vệ tinh theo chu kỳ 1 giờ, sử dụng trong khoảng 4 giờ. |
| Chính xác (Precise) |
$0{,}05 - 0{,}2\text{ m}$ |
Do IGS cung cấp sau vài ngày đến 2 tuần; nghiệm cuối cùng đạt $\pm 5\text{ cm}$. |
Tác giả phân tích chi tiết các hệ đo thời gian đóng vai trò quyết định trong việc đồng bộ tín hiệu vệ tinh:
- Thời gian thiên văn và thời gian Trái Đất: Thời gian quốc tế UT1 dựa trên sự quay thực tế của Trái Đất. Giờ UTC duy trì bằng đồng hồ nguyên tử nhưng được chèn thêm giây nhuận khi $|dUT1| > 0{,}9\text{ s}$.
- Thời gian GPS (GPST): Không có giây nhuận, giữ độ lệch cố định 19 giây so với thời gian nguyên tử quốc tế ($IAT = GPST + 19\text{s}$) và trùng với UTC tại thời điểm chuẩn 00:00:00 ngày 06/01/1980. Đến tháng 06/2000, sau khi chèn 32 giây nhảy ($n=32$), GPST chạy trước UTC đúng 13 giây.
- Ngày Julian và tuần GPS: Trình bày thuật toán chuyển đổi giữa ngày dân sự $(Y, M, D, UT)$ sang Ngày Julian (JD), Ngày Julian sửa đổi (MJD) và số tuần GPS ($\text{WEEK} = \text{INT}[(\text{JD} - 2444244{,}5)/7]$), lưu ý hiện tượng quay vòng số tuần GPS sau chu kỳ 1024 tuần (lần đầu diễn ra vào đêm 21-22/08/1999).
Về chuyển động quỹ đạo, đồ án trích xuất phương trình vi phân chuyển động 2 vật thể:
$$\ddot{r} + \frac{G(m_1 + m_2)}{r^3} r = 0$$
với tham số hấp dẫn chuẩn WGS-84 là $\mu = GM_E = 3986004{,}418 \times 10^8 \text{ m}^3\text{s}^{-2}$.
Để cung cấp lịch chính xác, mạng lưới của Trung tâm Dịch vụ GNSS Quốc tế (IGS) vận hành hơn 220 trạm thu 2 tần số toàn cầu gắn với các hệ thống đo giao thoa vô tuyến liên lục địa (VLBI) và đo khoảng cách laser vệ tinh (SLR).
Đồ án đưa ra công thức đánh giá tương quan giữa sai số lịch vệ tinh và sai số độ dài baseline:
$$\Delta x \approx \frac{l}{d} \Delta X$$
(trong đó $l$ là chiều dài cạnh đo, $d \approx 20.000\text{ km}$ là khoảng cách từ vệ tinh đến bề mặt Trái Đất, $\Delta X$ là sai số vị trí vệ tinh).
Bảng phân tích ảnh hưởng sai số quỹ đạo trong văn bản cho thấy: khi sai số lịch vệ tinh là $2{,}5\text{ m}$, sai số gây ra cho cạnh đo dài $1000\text{ km}$ lên tới $100\text{ mm}$ ($1 \times 10^{-7}$). Khi sử dụng lịch vệ tinh chính xác của IGS (sai số cỡ vài cm), sai số cạnh được triệt tiêu xuống mức một vài milimét, điều này khẳng định tính bắt buộc của việc tích hợp lịch IGS đối với mạng lưới trắc địa hạng cao và mạng lưới quan trắc địa động lực.
Cuối cùng, chương 2 phân tích chi tiết cấu trúc định dạng của tệp lịch dự báo YUMA (gồm 11 tham số chính như PRN, Health, Eccentricity, Time of Applicability, $\sqrt{a}$, $\Omega$, $\omega$, $M_0$, $a_0, a_1$) và cấu trúc 16 tham số của thông điệp lịch quảng bá vệ tinh (bổ sung các hệ số hiệu chỉnh điều hòa bậc hai $C_{uc}, C_{us}, C_{rc}, C_{rs}, C_{ic}, C_{is}$ do trường hấp dẫn Trái Đất không cầu và áp suất bức xạ Mặt Trời gây ra).
Kết quả và đóng góp
Đồ án đã tổng hợp và chuẩn hóa các kết quả nghiên cứu lý thuyết và thực hành xử lý số liệu GNSS:
- Hệ thống hóa quy trình công nghệ chuẩn: Thiết lập sơ đồ chi tiết các bước xử lý số liệu từ thực địa đến báo cáo kỹ thuật trên phần mềm Trimble Business Center 2.0, chỉ rõ tiêu chí đánh giá cạnh thông qua hệ số Ratio ($\ge 2-3$) và Variance ($1-8$).
- Lượng hóa ảnh hưởng của nguồn số liệu quỹ đạo: Đưa ra bảng số liệu định lượng về độ chính xác của các loại lịch vệ tinh, chứng minh bằng toán học công thức tỷ lệ $\Delta x \approx (l/d)\Delta X$. Với các cạnh đo dài trên $50\text{ km}$ hoặc mạng lưới hạng II, việc dùng lịch quảng bá thông thường (sai số $1-5\text{ m}$) không đáp ứng được yêu cầu sai số khép tương đối $1/1.000.000$, do đó bắt buộc phải sử dụng lịch vệ tinh chính xác IGS với độ chính xác $\pm 5\text{ cm}$.
- Đóng góp học thuật: Đồ án là tài liệu tổng hợp đầy đủ các công thức tính chuyển toán học giữa hệ WGS-84, hệ tọa độ trắc địa elipsoid và hệ tọa độ phẳng UTM/Gauss-Krüger, kết hợp mô hình xác định 7 tham số Bursa-Wolf và hiệu chỉnh Geoid địa phương, giúp người xử lý số liệu hiểu rõ bản chất giải thuật phía sau các giao diện phần mềm thương mại.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế: Nội dung văn bản tập trung chủ yếu vào cơ sở lý thuyết, mô hình toán học và quy trình vận hành phần mềm TBC 2.0; phần thực nghiệm so sánh đối chiếu số liệu đo đạc thực tế của một mạng lưới cụ thể giữa hai phương án dùng lịch quảng bá và lịch chính xác chưa được trình bày trọn vẹn trong phần văn bản này.
- Hướng nghiên cứu tiếp: Mở rộng nghiên cứu thực nghiệm đo đạc một lưới GNSS cụ thể ngoài thực địa có chiều dài cạnh đa dạng (từ cạnh ngắn dưới $5\text{ km}$ đến cạnh dài hàng trăm km), tiến hành xử lý song song với hai loại tệp lịch (Broadcast ephemeris *.nav và IGS Precise ephemeris *.sp3), từ đó lập bảng so sánh sai số thực tế tọa độ không gian $(X,Y,Z)$ và độ cao thủy chuẩn sau bình sai.
Giá trị tham khảo
Đồ án là tài liệu kỹ thuật hữu ích cho:
- Sinh viên ngành Kỹ thuật Trắc địa - Bản đồ, Địa tin học (Geoinformatics), Công trình: Dùng làm tài liệu tham khảo chi tiết về quy trình đo ngắm GNSS tĩnh, phương pháp trút số liệu, các bước bình sai lưới và cách vận hành các module tính toán trong phần mềm Trimble Business Center.
- Kỹ sư trắc địa thực địa và nội nghiệp: Tham khảo các tiêu chuẩn kiểm tra sai số khép hình cho lưới hạng II, hạng III; quy tắc can thiệp giải cạnh (loại bỏ vệ tinh góc cao thấp, khắc phục trượt chu kỳ Cycle Slip); và nguyên lý lựa chọn tệp Geoid phù hợp để chuyển đổi độ cao trắc địa sang độ cao thủy chuẩn.
Câu hỏi thường gặp
1. Tỷ số Ratio và Phương sai chuẩn (Variance) trong phần mềm xử lý cạnh GNSS có ý nghĩa gì?
Phương sai chuẩn (Variance) là sai số trung phương trọng số đơn vị của lời giải cạnh (không có đơn vị), giá trị lý tưởng bằng 1; đối với máy 1 tần số có thể chấp nhận từ 1 đến 8 hoặc 9. Tỷ số Ratio phản ánh tỷ lệ phương sai giữa lời giải tốt nhất và lời giải kế cận; tỷ số này phải đạt tối thiểu từ 2 đến 3 trở lên để đạt được lời giải cố định số nguyên đa trị pha (Fixed solution) có đủ độ tin cậy.
2. Lịch vệ tinh chính xác (Precise ephemeris) của IGS có độ chính xác bao nhiêu và mất bao lâu để có được?
Lịch vệ tinh chính xác IGS đạt độ chính xác từ $0{,}05\text{ m}$ đến $0{,}2\text{ m}$, trong đó lời giải cuối cùng (Final orbit) đạt độ chính xác xấp xỉ $\pm 5\text{ cm}$. Dữ liệu này được các trung tâm phân tích tính toán lại và cung cấp sau thời điểm quan sát từ 2 ngày đến 2 tuần.
3. Tại sao khi chiều dài cạnh đo GNSS càng lớn thì việc sử dụng lịch vệ tinh chính xác càng trở nên bắt buộc?
Theo công thức $\Delta x \approx (l/d)\Delta X$, sai số thành phần cạnh đo $\Delta x$ tỷ lệ thuận với chiều dài cạnh $l$. Với khoảng cách vệ tinh $d \approx 20.000\text{ km}$, sai số lịch quảng bá $2{,}5\text{ m}$ sẽ gây sai số tới $100\text{ mm}$ trên cạnh dài $1000\text{ km}$. Do đó, để đạt độ chính xác cỡ vài milimét trên các cạnh đo dài, bắt buộc phải sử dụng lịch vệ tinh chính xác để triệt tiêu sai số quỹ đạo $\Delta X$.
4. Giờ hệ thống GPS (GPST) khác biệt như thế nào so với Giờ phối hợp quốc tế (UTC)?
Giờ GPST là thang thời gian nguyên tử liên tục, không chèn giây nhuận, giữ độ lệch cố định 19 giây so với Thời gian nguyên tử quốc tế ($IAT = GPST + 19\text{s}$). GPST trùng khớp với UTC vào ngày 06/01/1980. Do UTC liên tục chèn các giây nhuận theo chu kỳ quay của Trái Đất, tính đến tháng 06/2000 (với 32 giây nhảy $n=32$), GPST chạy trước UTC đúng 13 giây.
Kết luận
Đồ án tốt nghiệp của sinh viên Phạm Thanh Bình đã trình bày một cách khoa học, chi tiết toàn bộ cơ sở lý thuyết toán học và quy trình công nghệ xử lý số liệu lưới trắc địa GNSS. Nghiên cứu đã làm rõ cấu trúc dữ liệu, cơ chế đồng bộ thời gian và giải thuật quỹ đạo Kepler của các loại lịch vệ tinh. Công trình khẳng định tầm quan trọng của việc kết hợp phần mềm Trimble Business Center với lịch vệ tinh chính xác IGS nhằm bảo đảm độ chính xác tối ưu cho các mạng lưới trắc địa quốc gia và quan trắc biến dạng công trình.
