Giới thiệu dự án

Trong tiến trình hiện đại hóa công tác quản lý đất đai và xây dựng cơ sở dữ liệu địa chính số tại Việt Nam, việc thành lập mạng lưới khống chế mặt bằng đóng vai trò tiên quyết nhằm đảm bảo độ chính xác hình học cho toàn bộ bản đồ địa chính. Theo thống kê của ngành Tài nguyên và Môi trường, việc ứng dụng công nghệ trắc địa truyền thống bằng đường chuyền kinh vĩ hoặc toàn đạc điện tử tại các khu vực địa hình bán sơn địa thường tốn từ 60% đến 70% thời gian cho công tác thông hướng, dọn tuyến và chuyển trạm máy, đồng thời dễ tích lũy sai số đo góc trên các tuyến chuyền kéo dài.

Xã Nhạo Sơn, huyện Sông Lô, tỉnh Vĩnh Phúc có tổng diện tích tự nhiên 366,99 ha (theo số liệu kiểm kê đất đai), địa hình đặc thù gồ ghề, đất đồi thấp xen kẽ ruộng trũng và các khu dân cư manh mún. Việc đo vẽ ranh giới thửa đất phục vụ cấp Giấy chứng nhận quyền sử dụng đất tại đây gặp nhiều trở ngại do tầm nhìn bị che khuất và cơ sở hạ tầng mốc khống chế cũ bị suy giảm mật độ. Đề tài tập trung giải quyết bài toán: Ứng dụng công nghệ định vị vệ tinh GPS (Global Positioning System) tĩnh tương đối để thiết kế, đo đạc và bình sai mạng lưới khống chế đo vẽ mặt bằng đạt chuẩn kỹ thuật quốc gia.

Mục tiêu của dự án:

  1. Khảo sát và đánh giá toàn diện: Phân tích đặc điểm địa hình, kinh tế - xã hội tại 7 thôn thuộc xã Nhạo Sơn và thu thập hệ thống 7 mốc tọa độ gốc (03 điểm GPS quốc gia: 91478, 91516, 91517 và 04 điểm địa chính cấp huyện: SL35, SL36, SL37, SL41).
  2. Thiết kế và chôn mốc mạng lưới đo vẽ: Thành lập mới 62 điểm khống chế mặt bằng (ký hiệu từ KV-01 đến KV-62) dạng mạng tam giác dày đặc phủ kín 366,99 ha.
  3. Thu nhận và xử lý tín hiệu GPS tĩnh: Áp dụng phương pháp đo tương đối tĩnh pha sóng tải L1/L2 với máy thu GPS-Huace X20, xử lý vector cạnh đáy bằng phần mềm Compass và bình sai trên hệ tọa độ quốc gia VN-2000.
  4. Biên tập và nghiệm thu kỹ thuật: Xuất 07 bảng biểu chuẩn hóa theo Thông tư số 25/2014/TT-BTNMT qua phần mềm DPsurvey 2.8, bảo đảm sai số trung phương vị trí điểm $M_p \le 5{,}0\text{ cm}$.

Phạm vi và giới hạn: Dự án triển khai thực nghiệm đo đạc mặt bằng trong thời gian từ ngày 10/08/2015 đến 30/11/2015 trên địa bàn xã Nhạo Sơn. Kết quả độ cao trong nghiên cứu được xác định đồng thời từ dữ liệu xử lý GPS-Leveling kết hợp với mốc độ cao nhà nước hiện có.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi công nghệ GNSS/GPS được phổ cập, công tác thành lập lưới địa chính đo vẽ chủ yếu dựa vào phương pháp đường chuyền đa giác hoặc mạng lưới tam giác đo góc - cạnh truyền thống.

Tiêu chí so sánh Đo đường chuyền toàn đạc Mạng tam giác đo góc - cạnh Đo GPS tĩnh tương đối
Yêu cầu thông hướng Bắt buộc phải thông hướng giữa tất cả các điểm liền kề Bắt buộc thông hướng giữa các đỉnh tam giác Không đòi hỏi thông hướng trực tiếp giữa các trạm máy
Ảnh hưởng thời tiết Phụ thuộc mạnh vào khúc xạ ánh sáng, sương mù, mưa gió Phụ thuộc tầm nhìn xa và điều kiện khí tượng Hoạt động liên tục 24/7, ít chịu tác động thời tiết
Tốc độ thi công Chậm (phải phát quang, dịch chuyển máy liên tục) Trung bình (cần xây dựng tiêu ngắm cao tầng) Rất nhanh (40 - 60 phút/ca đo đồng thời nhiều máy)
Tích lũy sai số Tăng dần theo chiều dài tuyến đường chuyền Phụ thuộc hình dạng hình học tam giác Độc lập theo khoảng cách vector cạnh, sai số phân bố đều
Độ chính xác mặt bằng $1/5.000 - 1/10.000$ $1/10.000 - 1/20.000$ Đạt $> 1/25.000$ (cỡ mm đến cm)

Yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must-have (Bắt buộc): Chuyển đổi dữ liệu sang hệ tọa độ VN-2000 (kinh tuyến trục $105^\circ00'$, múi chiếu $3^\circ$). Sai số trung phương vị trí điểm sau bình sai không vượt quá $5\text{ cm}$ (Lưới cấp 1) và $7\text{ cm}$ (Lưới cấp 2).
  • Should-have (Nên có): Lọc bỏ tự động các vệ tinh có góc cao vệ tinh $< 15^\circ$ hoặc hệ số suy giảm độ chính xác $GDOP > 4$.
  • Could-have (Có thể có): Đồng bộ tự động trút tệp thô .HCN sang định dạng chuẩn quốc tế RINEX để kiểm tra chéo trên các phần mềm mở.
  • Won't-have (Không áp dụng): Đo động thời gian thực RTK tầm xa khi chưa có trạm tham chiếu hoạt động liên tục (CORS) phủ sóng toàn diện.

Thiết kế hệ thống

Kiến trúc xử lý định vị không gian được tổ chức chặt chẽ từ thu tín hiệu vệ tinh đến xuất cơ sở dữ liệu địa chính:

[Chùm vệ tinh NAVSTAR GPS] 
[Phần mềm Compass v3.0]
[Phần mềm DPsurvey v2.8]

Thông số thiết lập hệ tọa độ VN-2000:

  • Ellipsoid quy chiếu: WGS-84 sửa đổi ($a = 6378137\text{ m}$, độ dẹt $\alpha = 1 / 298{,}257223563$).
  • Kinh tuyến trục (Central Meridian): $105^\circ00'$ (áp dụng cho tỉnh Vĩnh Phúc).
  • Múi chiếu: $3^\circ$ (Hệ số biến dạng chiều dài $k_0 = 0{,}9999$).
  • Gốc tọa độ quy ước: $X_0 = 0\text{ m}$, $Y_0 = 500.000\text{ m}$.
  • Bộ thông số tính chuyển 3 tham số (Translation): $\Delta X = +191{,}904\text{ m}$, $\Delta Y = +39{,}302\text{ m}$, $\Delta Z = +111{,}450\text{ m}$.

Methodology

Quy trình triển khai áp dụng phương pháp luận khảo sát trắc địa nghiêm ngặt chuẩn quốc gia gồm 5 giai đoạn:

Implementation và kết quả

Development process

1. Thu nhận tín hiệu ngoại nghiệp:

Tổ chức đo đạc gồm 04 máy thu GPS-Huace X20 hoạt động đồng thời trong 53 ca đo độc lập. Thời gian thu tín hiệu tại mỗi trạm máy đạt từ $40$ đến $60\text{ phút}$. Chiều cao anten được đo bằng thước thép chuyên dụng với độ chính xác đến $1\text{ mm}$ tại hai thời điểm: bắt đầu ca đo và kết thúc ca đo (lấy giá trị trung bình).

2. Thuật toán xử lý sai phân và giải vector cạnh đáy:

Mô hình toán học của phương pháp đo tương đối dựa trên phương trình sai phân bậc hai (Double Difference - DD) loại trừ triệt để sai số đồng hồ vệ tinh và đồng hồ máy thu:

$$\Delta \nabla \Phi_{AB}^{jk}(t) = \left( \Phi_B^k(t) - \Phi_B^j(t) \right) - \left( \Phi_A^k(t) - \Phi_A^j(t) \right) = \frac{1}{\lambda} \Delta \nabla \rho_{AB}^{jk}(t) + \Delta \nabla N_{AB}^{jk} - \Delta \nabla I_{AB}^{jk} + \Delta \nabla T_{AB}^{jk} + \epsilon_{\Delta \nabla \Phi}$$

Trong đó:

  • $\Phi_A^j, \Phi_B^k$: Pha sóng tải quan trắc từ trạm $A, B$ tới vệ tinh $j, k$.
  • $\Delta \nabla \rho_{AB}^{jk}$: Sai phân bậc hai khoảng cách hình học từ vệ tinh đến trạm thu.
  • $\Delta \nabla N_{AB}^{jk}$: Số nguyên đa trị (Ambiguity) sai phân bậc hai.
  • $\Delta \nabla I, \Delta \nabla T$: Ảnh hưởng khúc xạ tầng điện ly và đối lưu (được triệt tiêu trên các cạnh ngắn $< 5\text{ km}$).

Hệ phương trình số gia hiệu chỉnh tọa độ được giải bằng phương pháp Số bình phương nhỏ nhất (Least Squares Adjustment):

$$V = AX - L \quad \text{với điều kiện} \quad V^T P V = \min$$

Nghiệm của vector tọa độ điều chỉnh $X$:

$$X = (A^T P A)^{-1} A^T P L$$

import numpy as np

def least_squares_network_adjustment(A: np.ndarray, P: np.ndarray, L: np.ndarray):
    """
    Tính toán bình sai lưới trắc địa mặt bằng theo nguyên lý số bình phương nhỏ nhất
    :param A: Ma trận hệ số phương trình số gia (m x n)
    :param P: Ma trận trọng số trị đo (m x m)
    :param L: Vector số hạng tự do (m x 1)
    :return: Vector ẩn số tọa độ X, vector số hiệu chỉnh V, sai số trọng số đơn vị sigma_0
    """
    # 1. Thành lập hệ phương trình chuẩn: (A^T * P * A) * X = A^T * P * L
    N = A.T @ P @ A
    U = A.T @ P @ L
    
    # 2. Giải nghiệm ẩn số tọa độ
    X = np.linalg.inv(N) @ U
    
    # 3. Tính vector số hiệu chỉnh trị đo
    V = A @ X - L
    
    # 4. Tính sai số trung phương trọng số đơn vị sigma_0
    m, n = A.shape
    degrees_of_freedom = m - n
    vpv = V.T @ P @ V
    sigma_0 = np.sqrt(vpv / degrees_of_freedom)
    
    # 5. Ma trận hiệp phương sai của các ẩn số
    Q_xx = np.linalg.inv(N)
    sigma_X = sigma_0 * np.sqrt(np.diag(Q_xx))
    
    return {
        "X_coords": X,
        "Residuals_V": V,
        "Sigma_0": sigma_0,
        "Sigma_X": sigma_X
    }

# Minh họa ma trận trị đo kiểm thử cho một tam giác độc lập
if __name__ == "__main__":
    A_mock = np.array([[1.0, 0.0], [0.0, 1.0], [1.0, -1.0]])
    P_mock = np.diag([1.0, 1.0, 0.5])
    L_mock = np.array([[0.012], [-0.008], [0.019]])
    res = least_squares_network_adjustment(A_mock, P_mock, L_mock)
    print(f"Số gia tọa độ bình sai (m): \n{res['X_coords']}")
    print(f"Sai số đơn vị trọng số: {res['Sigma_0']:.4f}")

Testing và validation

Quá trình kiểm định chất lượng số liệu được thực hiện qua 3 cấp kiểm tra độc lập:

  1. Kiểm tra đồ hình vệ tinh và SNR: Loại bỏ các vệ tinh có tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio) $SNR < 30\text{ dB-Hz}$ hoặc bị phản xạ tín hiệu do tán cây, đường dây điện 110kV chạy qua xã.
  2. Kiểm tra sai số khép hình học: Sai số khép tương đối của các tam giác độc lập ($W_s$) đều thỏa mãn tiêu chuẩn: $W_s \le 1 / 15.000$.
  3. Kiểm tra độ chính xác theo quy chuẩn Thông tư 25/2014/TT-BTNMT:
Chỉ tiêu kỹ thuật Quy chuẩn Thông tư 25/2014/TT-BTNMT (Cấp 1) Kết quả thực tế đạt được tại xã Nhạo Sơn Đánh giá
Sai số trung phương vị trí điểm ($M_p$) $\le 5{,}0\text{ cm}$ Lớn nhất: $2{,}14\text{ cm}$ (tại điểm KV31) Vượt chuẩn 57,2%
Sai số tương đối cạnh sau bình sai ($M_s/S$) $\le 1 / 25.000$ $1 / 27.834$ (KV31-KV62) & $1 / 28.245$ (KV55-KV57) Đạt yêu cầu
Sai số khép tam giác lớn nhất $\le 1 / 10.000$ $1 / 18.255$ (Tam giác KV27-KV30-KV62) Đạt yêu cầu
Chiều dài cạnh đo trung bình ($S_{tb}$) $250 - 500\text{ m}$ $342{,}45\text{ m}$ (Cạnh dài nhất: $1123{,}79\text{ m}$) Tối ưu hình học

Kết quả đạt được

  • Hoàn thành thiết lập và xác định tọa độ chính xác cao cho 62 điểm khống chế đo vẽ mới (KV-01 đến KV-62), liên kết đồng bộ với 7 điểm gốc nhà nước.
  • 100% vector cạnh được tính toán thành công dạng Fixed Baseline (lời giải nguyên số đa trị cố định).
  • Xuất thành công hồ sơ kỹ thuật gồm 07 bảng biểu chuẩn mực từ phần mềm DPsurvey: Bảng tọa độ điểm gốc, Bảng kết quả đo đại lượng góc - cạnh, Bảng tính sai số khép hình học, Bảng kết quả tọa độ bình sai VN-2000, Bảng độ cao trắc địa và Sơ đồ lưới khống chế mặt bằng.

Đổi mới và đóng góp

  1. Ứng dụng giải pháp cắt lọc vệ tinh cục bộ (Satellite Deselection): Đối với các trạm máy đặt gần địa vật che khuất (rừng cây bạch đàn, đồi cao tại xóm Đồi), quy trình tiến hành phân tích biểu đồ Residuals và chủ động cắt bỏ các vệ tinh có tín hiệu suy hao, chỉ giữ lại 4-6 vệ tinh có chùm sóng sạch nhất, giúp tỷ lệ hội tụ giải pháp Fixed đạt 100%.
  2. Quy trình kết nối kép giữa phần mềm Compass và DPsurvey 2.8: Khắc phục nhược điểm của phần mềm ngoại quốc (Compass không hỗ trợ xuất báo cáo tiếng Việt theo mẫu Thông tư 25) bằng cách trích xuất tự động tệp cấu trúc .isb chuyển sang DPsurvey, rút ngắn 85% thời gian lập hồ sơ nghiệm thu.
  3. Cắt giảm chi phí và nhân lực thi công:
Phương pháp đo toàn đạc truyền thống:

Công nghệ GPS tĩnh tương đối (Đề tài ứng dụng):

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế:

  • Đo vẽ bản đồ địa chính tỷ lệ 1/500 và 1/1000: 62 điểm mốc KV được phân bố đều với bán kính khống chế $< 300\text{ m}$, cho phép đặt máy toàn đạc điện tử hoặc máy đo GPS RTK đo chi tiết ranh giới từng thửa đất mà không cần phát triển thêm trạm phụ.
  • Phục vụ công tác cấp Giấy chứng nhận quyền sử dụng đất: Cung cấp tọa độ pháp lý chuẩn xác để giải quyết tranh chấp đất đai tại các khu vực giáp ranh đồi rừng phức tạp của xã Nhạo Sơn.

Yêu cầu triển khai và khả năng mở rộng:

  • Trang thiết bị phần cứng: Tối thiểu 03 máy thu GNSS 1 hoặc 2 tần số (Huace X20, Trimble, Leica hoặc South), máy tính cá nhân cấu hình tối thiểu Intel Core i3, 4GB RAM.
  • Hệ thống phần mềm: Hệ điều hành Windows 7/10/11, Compass Geodetic Processing Software, DPsurvey 2.8+, AutoCAD/MicroStation SE/v8i kết hợp Famis.
  • Khả năng mở rộng: Mạng lưới khống chế 69 điểm này sẵn sàng làm mốc cơ sở để thiết lập hệ thống trạm Base RTK phục vụ công tác quy hoạch xây dựng nông thôn mới, cắm mốc quy hoạch giao thông Tỉnh lộ 307B.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật:

  • Hiện tượng đa đường dẫn tín hiệu (Multipath): Tại các vị trí thung lũng sâu hoặc điểm đặt sát vách đá đồi, tín hiệu GPS phản xạ gây kéo dài thời gian giải đa trị.
  • Phụ thuộc vào mốc gốc: Độ chính xác bình sai ràng buộc phụ thuộc hoàn toàn vào chất lượng hình học và tính ổn định của các mốc địa chính gốc huyện Sông Lô (SL35, SL36, SL37, SL41).

Hướng phát triển:

  • Tích hợp đa hệ vệ tinh GNSS: Nâng cấp sử dụng máy thu đa hệ thống vệ tinh (GPS + GLONASS + Galileo + BeiDou) để tăng số lượng vệ tinh quan sát đồng thời lên $> 20$ vệ tinh, giảm thời gian đo xuống còn $15 - 20\text{ phút/ca}$.
  • Chuyển giao công nghệ RTK mạng VN-CORS: Khai thác dữ liệu từ các trạm CORS quốc gia của Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam để đo vẽ trực tiếp tọa độ ranh giới thửa đất theo thời gian thực mà không cần thiết lập mốc đo vẽ trung gian.

Đối tượng hưởng lợi

  • Sinh viên chuyên ngành Quản lý đất đai / Kỹ thuật Trắc địa - Bản đồ: Tiếp cận toàn diện quy trình ngoại nghiệp - nội nghiệp chuẩn tắc, làm tài liệu tham khảo chất lượng cao cho các khóa luận tốt nghiệp.
  • Kỹ sư Trắc địa công trình: Nắm vững quy trình xử lý triệt tiêu sai số vệ tinh yếu trên phần mềm Compass và kỹ thuật bình sai lưới tự do/ràng buộc.
  • Ủy ban nhân dân xã Nhạo Sơn và Phòng TN&MT huyện Sông Lô: Sở hữu hệ thống mốc khống chế mặt bằng bền vững, phục vụ công tác chỉnh lý biến động đất đai và đền bù giải phóng mặt bằng nhanh chóng.

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

1. Cần điều kiện gì về mặt bằng khi chọn vị trí chôn mốc GPS đo vẽ?

Vị trí chọn mốc phải thông thoáng, góc mở lên bầu trời từ độ cao $15^\circ$ trở lên không có vật cản lớn (nhà cao tầng, cây rậm rạp). Điểm mốc phải cách xa trạm phát sóng vô văn, đường dây điện cao thế tối thiểu $50\text{ m}$ để tránh nhiễu điện từ và đặt trên nền đất ổn định, không bị xói lở.

2. Tại sao phải đo chiều cao anten 2 lần vào đầu và cuối ca đo?

Để loại trừ sai số do chân máy bị lún trong quá trình quan trắc hoặc sai số do thao tác đọc số của người đo. Hai lần đo lệch nhau không quá $2\text{ mm}$, giá trị trung bình cộng được nhập trực tiếp vào phần mềm Compass để tính chuyển chính xác tâm pha anten về tâm mốc trắc địa mặt đất.

3. Tại sao cần kết hợp cả bình sai tự do và bình sai ràng buộc?

Bình sai tự do nhằm đánh giá độ chính xác nội tại của các trị đo vector cạnh GPS, phát hiện các sai số thô mà không bị ảnh hưởng bởi sai số của các điểm gốc. Sau khi lưới tự do đạt tiêu chuẩn, tiến hành bình sai ràng buộc vào 7 điểm gốc nhà nước để đưa toàn bộ hệ thống tọa độ về quy chuẩn thống nhất VN-2000.

4. Phần mềm Compass xử lý hiện tượng trượt chu kỳ pha (Cycle Slip) như thế nào?

Compass sử dụng thuật toán tổ hợp tuyến tính không chứa tầng điện ly (Ionosphere-free Linear Combination) và sai phân bậc ba (Triple Difference) theo thời gian để phát hiện điểm nhảy bước sóng, tự động khôi phục số nguyên đa trị hoặc ngắt đoạn dữ liệu để giải riêng từng phân đoạn.

5. Chi phí đầu tư thiết bị và hiệu quả kinh tế so với toàn đạc điện tử?

Mặc dù chi phí đầu tư ban đầu cho máy thu GPS cao hơn máy toàn đạc điện tử khoảng 1,5 - 2 lần, nhưng hiệu suất đo đạc tăng gấp 3 lần, tiết giảm 50% nhân sự tại hiện trường và loại bỏ hoàn toàn chi phí phát quang thông hướng, giúp hoàn vốn đầu tư chỉ sau 2-3 dự án đo vẽ cấp xã.

Kết luận

Đề tài "Ứng dụng công nghệ GPS thành lập lưới khống chế mặt bằng đo vẽ tại xã Nhạo Sơn, huyện Sông Lô, tỉnh Vĩnh Phúc" đã giải quyết triệt để bài toán khó khăn về xây dựng mạng lưới trắc địa tại địa bàn bán sơn địa phức tạp. Với 62 điểm khống chế mới được thành lập đạt độ chính xác vị trí điểm lớn nhất $M_p = 2{,}14\text{ cm}$ (vượt xa giới hạn $5{,}0\text{ cm}$ của Thông tư 25/2014/TT-BTNMT), công trình khẳng định tính ưu việt vượt trội của công nghệ định vị vệ tinh tĩnh tương đối so với các phương pháp trắc địa quang học truyền thống.

Kết quả nghiên cứu không chỉ cung cấp cơ sở dữ liệu tọa độ không gian chính xác phục vụ công tác đo vẽ bản đồ địa chính và cấp Giấy chứng nhận quyền sử dụng đất tại địa phương, mà còn cung cấp một khung quy trình kỹ thuật hoàn chỉnh từ khâu khảo sát, phân ca đo, xử lý vector cạnh trên Compass đến biên tập tự động trên DPsurvey 2.8. Đây là nguồn tài liệu tham khảo có giá trị ứng dụng thực tiễn cao cho các kỹ sư, nhà quản lý và sinh viên ngành Quản lý đất đai và Kỹ thuật Trắc địa.