Phối hợp nghiên cứu để xây dựng và hoàn thiện hệ thống trọng lực phục vụ công tác nghiên cứu lãnh thổ và khai thác tài nguyên, bảo vệ môi trường ở Việt nam

Phối hợp nghiên cứu để xây dựng và hoàn thiện hệ thống trọng lực phục vụ công tác nghiên cứu lãnh thổ và khai thác tài nguyên, bảo vệ môi trường ở Việt nam.

BTNMT VNCĐC Bộ tài nguyên môi trờng viện nghiên cứu địa Đờng Hoàng Quốc Việt, Cầu Giấy, Hà Nội báo cáo tổng kết khoa học kỹ thuật đề tài hợp tác quốc tế khoa học công nghệ môi trờng: Phối hợp nghiên cứu để xây dùng vµ hoµn thiƯn hƯ thèng träng lùc phơc vơ công tác nghiên cứu lÃnh thổ khai thác tài nguyên, bảo vệ môi trờng Việt nam Chủ nhiệm: TS Lê Minh 6736 12/2/2008 Hà Nội, 12-2007 Bản quyền thuộc VNCĐC Đơn xin chép toàn phần tài liệu phải gửi đến Viện trởng VNCĐC trừ trờng hợp sử dụng với mục đích nghiên cứu BTNMT VNCĐC Bộ tài nguyên môi trờng viện nghiên cứu địa Đờng Hoàng Quốc Việt, Cầu Giấy, Hà Nội báo cáo tổng kết khoa học kỹ thuật đề tài hợp tác quốc tế khoa học công nghệ môi trờng: Phối hợp nghiên cứu để xây dựng hoàn thiện hệ thống trọng lực phục vụ công tác nghiên cứu lÃnh thổ khai thác tài nguyên, bảo vệ môi trờng Việt nam Hà Nội, ngày tháng năm 2007 Chủ nhiệm Đề tài Hà Nội, ngày tháng năm 2007 Cơ quan chủ trì Đề tài Q.Viện trởng Viện Nghiên cứu Địa TS Lê Minh TS Ngun Dịng TiÕn Hµ Néi, 12-2007 danh sách ngời thực STT Họ tên Cơ quan TS Lê Minh Trung tâm Viễn thám KS Nguyễn Tuấn Anh Trung tâm Viễn thám KS Điều Văn Vân Viện Nghiên cứu Địa KS Trần Đình ấu Viện Nghiên cứu Địa KS Lê Thanh Hải Cục Địa chất Khoáng sản KS Phan Xuân Hậu Viện Nghiên cứu Địa TS Trần Đình Lữ Viện Nghiên cứu Địa TS Maiorov MA Viện Nghiên cứu Trắc địa Bản đồ Liên bang Nga TS Kopachevski Xí nghiệp Trắc địa ảnh Mátcơva - Liên bang Nga 10 TS Nhiverov Xí nghiệp Trắc địa ảnh Mátcơva - Liên bang Nga Bài tóm tắt Đề tài Hợp tác Quốc tế Khoa học công nghệ Môi trờng Phối hợp nghiên cứu để xây dựng hoàn thiện hệ thống trọng lực phục vụ công tác nghiên cứu lnh thổ khai thác tài nguyên, bảo vệ môi trờng Việt Nam đề tài hợp tác với Liên hiệp Trắc địa Bản đồ ảnh Liên bang Nga nhằm mục đích xây dựng hoàn thiện hệ thống trọng lực quốc gia nớc ta Mục tiêu đề tài phối hợp hợp tác với quan đo đạc đồ Liên bang Nga nhằm thực dự án Xây dựng hoàn thiện hƯ thèng träng lùc qc gia ë ViƯt Nam” ®· đợc Bộ Tài nguyên Môi trờng phê duyệt đa vào thực từ năm 2003 Ngoài ra, xây dựng sở khoa học cho việc thống hệ thống độ cao khu vực thông qua việc ứng dụng phơng pháp đo cao vệ tinh Đề tài đà giải đợc nội dung sau: - Xây dựng yêu cầu cho việc xây dựng hƯ thèng träng lùc nhµ n−íc bao gåm l−íi träng lực sở lới trọng lực hạng I Đa yêu cầu độ xác lới trọng lực phục vụ cho việc xác định thay đổi trọng lực với đại lợng thay đổi hàng năm từ 0,01 - 0,03 mGal cần phải xác địh độ xác giá trị trọng lực tuyệt đối 0,005 mGal Ngoài ra, xem xét đến độ xác trung bình dị thờng trọng lực theo ô chuẩn x phục vụ xây dựng mô hình Geoid độ xác cao Việt Nam Các yêu cầu tọa độ độ cao cần phải xác định cho mốc trọng lực sở trọng lực hạng I đợc nghiên cứu có khoa học ã ã - Đối với l−íi träng lùc c¬ së: MX,Y≤ 0,5 m, Mh ≤ 2,0 cm Đối với lới trọng lực hạng I: MX,Y 2,0 m, Mh 5,0 cm Trong báo cáo đà đề cập đến việc thu thập đánh giá tổng quan vỊ hƯ thèng träng lùc hiƯn cã ë n−íc ta CHDCND Lào Trên sở đánh giá đà đa yêu cầu kỹ thuật việc xây dựng hệ thống trọng lực nhà nớc bao gồm: ã Yêu cầu kỹ thuật xây dựng lới trọng lực sở; ã Yêu cầu kỹ thuật xây dựng lới trọng lực hạng I Trong phần đà xây dựng quy trình đo xử lý kết ®o träng lùc (®o träng lùc tut ®èi vµ ®o trọng lực tơng đối), đà xây dựng phần mềm tính toán bình sai trọng lực tuyệt đối lới trọng lực (chơng trình bình sai lới trọng lực đà đợc áp dụng cho lới trọng lực điểm tựa) Trên sở nghiên cứu đà đề xuất phơng án xây dựng lới trọng lực nhà nớc Việt Nam bao gồm: ã Lới trọng lực sở: có 12 điểm ã Lới trọng lực hạng I: có 28 điểm Lới trọng lực sở sử dụng máy trọng lùc tut ®èi GBL cđa Nga ®Ĩ ®o L−íi träng lực hạng I sử dụng máy lắc tơng đối AGAT (Nga) để đo Các lới trọng lực nhà nớc có hợp tác với Liên bang Nga để thực Trong thời gian thực Đề tài đà phối hợp với Xí nghiệp Trắc địa ảnh Matscơva đo đợc điểm trọng lực tuyệt đối điểm trọng lực hạng I Đề tài đà xây dựng phơng án thống hệ thống độ cao khu vực thông qua việc nghiên cứu phơng pháp xác định độ chênh h mặt nớc biển trung bình mặt Geoid ë ®iĨm khëi tÝnh ®é cao (ë n−íc ta ë điểm Hòn Dấu độ xác khoảng 0,9 m) Đà xây dựng phơng pháp thống hệ thống độ cao thông qua việc xử lý hỗn hợp liệu GPS - TC độ cao Geoid sở xây dựng lới GPS mốc độ cao quốc gia Đề tài đà đề xuất phơng pháp bình phơng tối thiểu Collocation việc xử lý liệu Các kết đề tài đà đợc Viện Nghiên cứu Địa đa vào áp dụng việc thực dự án Xây dựng hoàn thiện hệ thống trọng lực nhà nớc Mục lục Đặt vấn ®Ò Các yêu cầu cho l−íi träng lùc quèc gia 10 1.1 Các yêu cầu cho lới träng lùc quèc gia 10 1.2 Yêu cầu độ xác lới trọng lực qc gia 11 HiƯn tr¹ng cđa hệ thống trọng lực Việt Nam khả sử dụng vào mục đích trắc địa 15 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 HiƯn tr¹ng vỊ hƯ thèng träng lùc ë ViÖt Nam .15 §iĨm gèc träng lùc .16 Đờng đáy trọng lực .16 L−íi träng lùc c¬ së 17 L−íi träng lùc h¹ng I 17 L−íi träng lùc h¹ng II 18 Lới trọng lực hạng III träng lùc chi tiÕt 18 Các liệu trọng lực biển trọng lực vũ trụ .19 Các liệu trọng lùc ë Lµo 20 X©y dùng L−íi träng lùc qc gia 23 3.1 C¸c yêu cầu giải pháp kỹ thuật xây dựng lới trọng lực hạng cao 24 3.1.1 Yêu cầu kỹ thuật lới trọng lực quốc gia .24 3.1.2 Yêu cầu kỹ thuật lới trọng lực hạng I 25 Phơng án đo trọng lực độ xác cao 27 4.1 Đo trọng lực tuyệt đối 27 4.1.1 4.1.2 4.1.3 4.1.4 Nguyên tắc chung 27 Mét sè m¸y träng lùc tuyệt đối 28 Đo trọng lực tuyệt đối máy Lazer GBL (Nga) 30 Trình tự đo trọng lực máy GBL 33 4.2 §o träng lùc tơng đối 34 4.2.1 Nguyên tắc chung 34 TÝnh to¸n gi¸ trị trọng lực đánh giá độ xác 37 5.1 Tính toán giá trị trọng lực tut ®èi .37 5.2 Tính giá trị trọng lực tơng đối 40 5.2.1 Tính giá trị trọng lực từ kết đo máy trọng lực lắc tơng đối 40 5.2.2 Đánh giá độ xác kết đo máy trọng lực lắc tơng đối .41 5.2.3 Tính giá trị trọng lực từ kết đo máy trọng lực tĩnh .43 5.2.4 Số hiệu chỉnh kết đo vỊ t©m mèc 47 B×nh sai l−íi träng lùc 47 6.1 Nguyên tắc chung 47 Phơng án xây dựng lới trọng lực quốc gia ë ViÖt nam 48 7.1 7.2 7.3 7.4 7.5 ThiÕt kÕ l−íi träng lùc c¬ së 48 ThiÕt kÕ l−íi träng lùc h¹ng I 50 Đo trọng lực hạng I máy “Agat” 51 KiÓm tra kết đo ngoại nghiệp 53 Đo điểm trọng lực vệ tinh .53 KÕt đo trọng lực tuyệt đối 54 8.1 Kh¸i qu¸t chung 54 8.2 Kết đo trọng lực tuyệt đối 55 Phơng án đo thống hệ thống độ cao 56 9.1 Giíi thiƯu chung 56 9.2 C¬ së khoa häc cđa hƯ thèng ®é cao 56 9.3 Phơng pháp thống hệ ®é cao .63 10 Kết luận Kiến nghị 67 10.1 KÕt luËn 67 10.2 KiÕn nghÞ 68 Tài liệu tham khảo 70 Phụ lục Đặt vấn đề Trọng lực ngành khoa học trái đất nghiên cứu xác định trờng trọng lực trái đất xác định hình dạng trái đất Trọng lực ngành khoa học chung trái đất nhng phần tách biệt đợc Trắc địa Trọng lực Trắc địa đóng vai trò to lớn, giải nhiều toán Trắc địa nh xác định độ lệch dây dọi dị thờng độ cao trọng lực để chỉnh lý trị đo góc trị đo cạnh mặt ellipsoid quy chiếu Dị thờng độ cao trọng lực đóng vai trò quan trọng việc xác định mô hình Geoid cục toàn cầu Dị thờng trọng lực phục vụ cho việc xác định số cải trọng lực cho việc thành lập hệ độ cao chuẩn Nhờ có giá trị trọng lực thời kỳ năm 70 80 kỷ trớc đà xác định đợc độ xác độ lệch dây dọi 0,5 dị thờng độ cao đạt độ xác 3,0m đáp ứng đợc việc chỉnh lý hệ thống trắc địa công nghệ cổ truyền Ngày nay, với phát triển khoa học trắc địa phạm vi toàn cầu, nhiều toán trắc địa đợc giải víi sù tham gia cđa nhiỊu n−íc ®· cho kÕt có độ xác khả ứng dụng cao vào nhiều nhiệm vụ nghiên cứu khoa học phát triển kinh tế nớc Hệ thống số liệu trọng lực toàn cầu đà đợc thiết lập với độ xác cao nhiều công nghệ phơng pháp khác Trớc hết, công nghệ trọng lực mặt đất ngày thiết bị đo trọng lực đất liền biển đà có cải tiến đáng kể Với thiết bị đo mặt đất, hệ thống máy trọng lực tuyệt đối hệ đời có kết hợp công nghệ điện tử tin học cho phép xác định độ xác giá trị trọng lực tuyệt đối tới 2àKGal máy trọng lực tuyệt đối FG-5 (Mỹ) 5àKGal với máy trọng lực tuyệt đối GBL (Nga) Độ xác mở ứng dơng hÕt søc quan träng cđa träng lùc cho nghiªn cứu xây dựng mô hình Geoid độ xác cao tới 1cm, xác định thay đổi giá trị träng lùc sù biÕn ®ỉi thÊt th−êng cđa vËt chất nh chuyển động vỏ trái đất sở để phát động đất Sự biến ®ỉi cđa khÝ hËu thêi tiÕt, sù chun ®éng cđa vật chất lòng đại dơng thấy đợc qua việc xác định biến đổi trọng lực Hệ thống trọng lực tơng đối có cải tiến phát triển đáng kể Thế hệ máy trọng lực tơng đối Lacoste Rombert (G) (Mỹ) có biên độ đo tới 7000 mGal độ xác đạt 0,01 mGal, máy Ascania (Tây Đức) có biên độ tới 5000 mGal, độ xác đạt tới 0,01 mGal cho phép thiết lập mạng lới träng lùc ph¹m vi tõng quèc gia, tõng khu vực với độ xác cao với mật ®é ®iĨm t ý Cịng cÇn l−u ý r»ng sù dao động (xê dịch) điểm máy nhỏ, vòng 24h đạt 0,1 - 0,2 mGal Hệ thống máy đo trọng lực biển đà đợc cải tiến đáng kể, trớc máy đo trọng lực biển đặt tàu hiệu ứng EVOST đạt độ xác từ 3-5 mGal ngày máy trọng lực đo đợc biển không đạt độ xác nhỏ mGal, loại máy hÃng Lacoste Rombert (Sea air Gravity) Độ xác đo trọng lực biển đợc nâng lên phần nhiều phụ thuộc vào việc tự động hoá hoàn toàn trình đo xử lý số liệu đo Hệ thống giám sát kiểm soát liệu đợc tự động hoá cao, hiệu ứng EVOST đà đợc xử lý tính toán đạt độ xác tới phần mời mGal Ngoài ra, việc định vị tàu công nghệ GPS có độ xác tới - 2m làm tăng đáng kể độ xác dị thờng trọng lực Hệ thống trọng lực toàn cầu đợc thiết lập không nhờ vào việc đo trọng lực đất liền biển mà nhờ vào công nghệ vũ trụ Nhờ vào kết đo liên tục khoảng cách từ vệ tinh tới mặt nớc biển ngời ta xác định đợc độ cao Geoid biển với độ xác tới decimeter xác định dị thờng trọng lực với độ xác nhỏ 10 mGal Trong giai đoạn nay, mô hình trọng trờng trái đất đà đợc thiết lập nh mô hình EGM-96 (Mỹ) GAO-97 (Nga) xác định đợc nhờ hệ số điều hoà cầu có bậc n = m =360, đồng thời nhờ hệ số điều hoà xác định dị thờng trọng lực toàn cầu với độ xác cho vùng đất liỊn tõ 20 - 30 mGal, cho vïng biĨn tõ 10mGal Với độ xác dị thờng trọng lực cho phép xác định độ cao Geoid với độ xác từ 0,5 - 1,5m phạm vi toàn cầu Từ năm 2004, Mỹ đà đa vào hoạt ®«ng hƯ thèng träng l−c vƯ tinh Grace cho phÐp xác định giá trị trọng lực phạm vi toàn cầu với độ xác cao gấp nhiều lần so với mô hình trọng trờng EGM-96 Độ xác trọng lực đạt tới 5-10 mgal Ngoài ra, ngày với việc sử dụng thiết bị đo trọng lực vũ trụ nh gradiometer xác định giá trị trọng lực với độ xác tới - mGal Việt Nam, từ năm 1998 đà tiến hành nghiên cứu sử dụng liệu trọng lực đo đợc đất liền, biển liệu từ mô hình trọng trờng trái đất phục vụ cho việc xây dựng mô hình Geoid đáp ứng cho nhiệm vụ chỉnh lý hệ thống lới toạ độ Việt Nam Hệ thống trọng lực Việt Nam đợc xây dựng từ năm 70 kỷ trớc, có điểm gốc trọng lực Hà Nội đợc đo nối với hệ thống trọng lực quốc tế thông qua điểm PoCoVo (Nga) với độ xác đo tơng đối 0,04 mGal Hệ thống trọng lực đà đợc truyền tới điểm trọng lực hạng I, hạng II phủ trùm nớc Năm 1998, hệ thống trọng lực Việt Nam đà đợc đại hoá, nhiên độ xác không đợc tăng lên Hiện nay, hệ thống trọng lùc cị ®· cã nhiỊu thay ®ỉi tíi 70% sè mốc trọng lực hạng I hạng II đà bị thực địa Các giá trị trọng lực, toạ độ độ cao cha đợc xác định với độ xác cần thiết Yêu cầu đòi hỏi phải xây dựng hệ thống trọng lực với độ xác cao theo tiêu chuẩn quốc tế đà đợc đặt cho công tác trắc địa Việt Nam để giải toán xây dựng mô hình Geoid độ xác cao cỡ vài centimeter, chí 1cm; vấn đề nghiên cứu chuyển dịch vỏ trái đất, phát sớm nguyên nhân động đất; vấn đề biến đổi vật chất, chuyển động thềm lục địa biến đổi môi trờng nh toán thăm dò dầu khí, khoáng sản, tài nguyên thiên nhiên đòi hỏi phải xây dựng hệ thống trọng lực có độ xác tới phần trăm mGal phần nghìn mGal đáp ứng thoả mÃn đợc nhiệm vụ nghiên cứu lÃnh thổ phát triển kinh tế, phòng tránh thiên tai thời gian năm Hệ thống trọng lực quốc gia nớc ta đợc xây dựng hoàn thành thời gian 2003- 2007 Trong dự án trình bày chuyên mục sau: Yêu cầu cho lới trọng lực quốc gia; Hiện trạng hệ thống trọng lực Việt Nam khả sử dụng trắc địa; Các giải pháp kỹ thuật công nghệ đo trọng lực hạng cao: ã Đo trọng lực tuyệt đối ã Đo trọng lực tơng đối Xây dựng phơng án xử lý toán học số liệu trọng lực hạng cao; Xây dựng phơng án đo trọng lực độ xác cao; Xây dựng quy trình công nghệ đo xử lý kết đọ trọng lực hạng cao; Chơng trình xử lý, tính toán kết đo trọng lực; Xây dựng phơng án đo xử lý thống hệ thống độ cao phụ thuộc vào loại đất có giá trị khoảng 17 àKGal meter h, hTB - Giá trị độ sâu độ sâu trung bình tầng nớc ngầm so với mặt đất Thông tin mặt tầng nớc ngầm đợc nhận từ tổ chức đào giếng khoan gần điểm đo 4.2 Tính giá trị trọng lực tơng đối Tính giá trị trọng lực từ kết đo máy trọng lực lắc tơng đối: Tính toán kết đo bao gồm bớc sau: ã Tính chu kỳ dao động lắc số hiệu chỉnh ã Tính hiệu giá trị trọng lực ã Đánh giá độ xác kết đo Chu kỳ dao động lắc đợc hiƯu chØnh tÝnh theo c«ng thøc sau: S = S ' + ∆S a + ∆S t + ∆S B + S f + S o (4.10) đây: S - Chu kỳ dao động lắc đo đợc xác định công thức S'= QS N f0 QS - Số xung động thạch anh phát khoảng thời gian tơng ứng với N dao động lắc f - Giá trị tần số chuẩn động Các số cải S a , ∆S t , ∆S B , ∆S f , S đợc lần lợt tính nh sau: - S a số cải biên độ để hiệu chỉnh chu kỳ dao động lắc tới biên độ vô bé, đợc tính công thức sau: S a = − S' 1024 ∆Y χ F sin 2 πQa (4.11) K a S ' f Trong đó: Qa - Tổng xung động phát khoảng thời gian tơng ứng với số K a dải sóng (vùng sóng) qua từ lắc đầu khe hở màng ngăn (tấm chắn) đo biên độ y - Khoảng cách đầu khe hở chắn tÝnh b»ng milimeter F - Tiªu cù kÝnh vËt tÝnh milimeter - Số chùm ánh sáng bị phản lại từ gơng cặp lắc - S t : Số cải nhiệt độ - S B : Số cải áp suất không khí tồn (sót lại) máy - S f : Số cải sai lệch tần số máy phát f với tần số chuẩn f0 19 S f = S ' f0 − f (4.12) f0 - ∆S : Sè c¶i chÝnh ¶nh h−ëng søc hút mặt trăng, mặt trời S đợc tính theo c«ng thøc: ∆S = − S' 2g ⎡ 1⎞ 1⎞ ⎤ ⎛ ⎛ 2 ⎢0,0760⎜ cos Z − ⎟ + 0,1647⎜ cos Z e − ⎟ K p ⎥ K y (4.13) ⎝ ⎠ ⎝ đây: Z0, Ze - Khoảng cách thiên đỉnh mặt trời mặt trăng Ky - Hệ số đàn hồi trái đất Kp - Thị sai mặt trăng - g số cải Khoncasalo đợc tính bổ sung vào kết tính toán giá trÞ träng lùc δg = 0,030569(1 − sin ϕ )Ky (4.14) ϕ - VÜ ®é ®iĨm ®o HiƯu giá trị trọng lực g (mGal ) điểm đo so với điểm khởi tính đợc xác định công thøc: ⎛ S2 ⎞ ∆g = g ⎜ 02 − 1⎟ (4.15) ⎜S ⎟ ⎝ ⎠ Hc: ∆g = λ (S − S ) + µ (S − S )2 + ε (S − S )3 (4.16) = đây: 2g0 3g 4g , = 20 , ε = − 30 S0 S0 S0 S, S0 - Giá trị trung bình chu kỳ dao động lắc điểm đo điểm khởi tính, đợc cải tất số hiệu chỉnh g, g0 - Giá trị trọng lực điểm đo điểm khởi tính Bình sai lới trọng lực Nguyên tắc chung Đối với tính toán bình sai trọng lực có trờng hợp sau: Đối với đo trọng lực tuyệt đối, phơng trình sai số có d¹ng: g k − g k = Vk (5.1) Trong đó: g k : giá trị trọng lực cần xác định điểm k g k : giá trị trọng lực đo đợc điểm k Đối với đo trọng lực tơng đối, phơng trình sai số có dạng: gk - gl - ∆gkl = Vkl (5.2) Trong ®ã: gk, gl: giá trị trọng lực cần xác định điểm k điểm l gkl: gia số trọng lực đo đợc điểm k điểm l Vk, Vkl: số hiệu chỉnh đến trị đo 20 Với n trị đo lập phơng trình sai số có dạng: V = AX - L (5.3) Giải phơng trình (6.3) theo phơng pháp bình phơng tối thiểu [pvv] = Ta có phơng trình chuẩn: AT pL = AT PAX Từ có: đây: AT pL = N n AT PA N = AT PA X = n = AT pL víi sai sè trung ph−¬ng sè đơn vị: = V T pv nu Sai số giá trị trọng lực i: i = µ Qll −1 Trong ®ã: Qll = (AT pv ) Giá trị trọng lực cần tính: g K = g K + δg K δg K : Èn số cần tìm Phơng án xây dựng lới trọng lùc qc gia ë ViƯt nam 6.1 ThiÕt kÕ l−íi trọng lực sở Nh đà trình bày vai trò lới trọng lực sở xác định xác giá trị trọng lực tuyệt đối g điểm trọng lực nhà nớc lÃnh thổ Việt Nam liên kết với hệ thống trọng lực quốc tế làm sở để nghiên cứu sù biÕn thiªn cơc bé cịng nh− sù biÕn thiªn khu vực toàn cầu giá trị trọng lực phục vụ cho toán xác định xác kích thớc hình dạng trái đất, cấu trúc vỏ trái đất, biến thiên vật chất lòng trái đất phục vụ cho mục đích quân việc xác định quỹ đạo chuyển động đạn pháo tầm xa quỹ đạo đờng tên lửa hành trình Mặt khác, vào trạng công nghệ khả thực công tác trọng lực nớc ta, việc xây dựng hệ thống điểm trọng lực tuyệt đối có khoảng cách điểm từ 200km - 400km phơng án khả thi đáp ứng đợc yêu cầu đặt trớc mắt lâu dài công tác trọng lực nớc ta Lới trọng lực sở đợc đo máy trọng lực tuyệt đối GBL Liên bang Nga có độ xác 0.005 mGal Tổng số điểm lới dự kiến ®−ỵc thiÕt kÕ bao gåm 12 ®iĨm, ®ã: Các điểm đà đo năm 2006 (4 điểm) Điểm gốc Hà Nội Điểm Điện biên Phủ Điểm Đà Nẵng Điểm thành phố Hồ Chí Minh 21 Các điểm đo (8 điểm ) §iĨm Lµo Cai (Sa Pa) − §iĨm Cao B»ng − §iĨm TP Vinh − §iĨm Nha Trang − §iĨm Bu«n Mê Thuật Điểm Lao Bảo Điểm Côn Đảo Điểm Phú Quốc Xung quanh điểm trọng lực sở đo thêm điểm trọng lực vệ tinh Các điểm trọng lực vệ tinh đợc đo máy GNUKV máy trọng lực có độ xác tơng đơng, độ xác ®iĨm träng lùc vƯ tinh ± 0.02 mGal ViƯc x©y dựng thêm điểm trọng lực tuyệt đối nớc ta tạo đợc lới trọng lực nhà nớc gồm 12 điểm, phân bố nớc có độ xác cao gấp 10 lần lới trọng lùc hiƯn cã ë n−íc ta (xem Phơ lơc 8) Lới trọng lực sở đáp ứng việc nghiên cứu biến thiên trọng lực hàng năm từ 0.020.04 mGal ë mäi khu vùc trªn l·nh thỉ n−íc ta phục vụ trực tiếp cho việc phát sớm tai biến thiên tai nh động đất, núi lửa giúp cho công tác cảnh báo thiên tai giảm thiểu thiệt hại cho đất nớc Trong số 12 điểm trọng lùc tut ®èi cã mét sè ®iĨm träng lùc tut đối đợc đặt khu vực nhạy cảm hay xảy động đất nh Điện Biên Phủ, Lào Cai hay gần TP Hồ Chí Minh Các điểm trọng lực tuyệt đối cần đợc đo lặp thờng xuyên với tần suất - năm/lần Ngoài ra, xung quanh điểm cần xây dựng mạng lới điểm trọng lực vệ tinh có từ - 10 điểm đợc phân bổ xung quanh điểm (khoảng cách không vợt 80 km) Ngoài ra, có điểm trọng lực tuyệt đối đợc xây dựng đảo Côn Đảo Phú Quốc phục vụ cho đo trọng lực biển phục vụ cho việc nghiên cứu biến thiên vỏ trái đất vùng biển nớc ta Điểm trọng lực tuyệt đối Lao Bảo việc phục vụ nghiên cứu biến thiên trọng lực khu vực làm sở để phát triển lới trọng lực quốc gia cho CHDCND Lµo 6.2 ThiÕt kÕ l−íi träng lùc hạng I Lới trọng lực hạng I đợc thiết kế dựa vào điểm lới trọng lực sở, khoảng cách điểm từ 100 km - 150 km, với mật độ trung bình từ 22 10.000km2 - 20.000km2/điểm Tổng số điểm trọng lực hạng I cần xây dựng gồm có 28 điểm, bao gồm điểm sau: STT Tên điểm Ghi Hà Giang xây mốc Lai Châu (Thị xà Cam Đờng) Lạng Sơn Yên Bái Nà Sản Thảo Nguyªn Tiªn Yªn -nt- KiÕn An -nt- Thái Bình (Tiền Hải) -nt- -nt-ntmốc cũ (đo trùng vào mèc träng lùc h¹ng I cị) -ntSư dơng mèc cị 10 Bạch Long Vĩ -nt- 11 Thọ Xuân -nt- 12 Mờng Xén -nt- 13 Đồng Hới -nt- 14 Phú Bình -nt- 15 Đảo Ly Sơn -nt- 16 Kon Tum -nt- 17 Phù Cát -nt- 18 Pleiku -nt- 19 Tuy Hòa -nt- 20 Liên Khơng -nt- 21 Xa Mát -nt- 22 Phan ThiÕt -nt- 23 Phó Q -nt- 24 Vịng Tµu sử dụng mốc cũ 25 Châu Đốc mốc 26 Cần Thơ -nt- 27 Cà Mau -nt- 28 Hà Tiên -nt- 23 Kết đo trọng lực tuyệt đối 7.1 Khái quát chung Theo kết hợp tác Viện Nghiên cứu Địa - Bộ Tài nguyên Môi trờng Việt Nam Xí nghiệp Trắc địa ảnh Matscơva - Liên bang Nga từ ngày 26 tháng 10 năm 2004 đến ngày 31 tháng 12 năm 2005 Phía Nga đà tiến hành đo trọng lực tuyệt đối máy lắc GBL Việt Nam Viện Nghiên cứu Địa đà tiến hành chọn điểm chôn mốc trọng lực tuyệt đối điểm, bao gồm: Điểm trọng lực gốc (Láng Hà Nội) Điểm trọng lực Điện Biên Phủ Điểm trọng lực Đà Nẵng Điểm trọng lực Nha Trang §iĨm träng lùc gèc TP Hå ChÝ Minh (28 Ngun Văn Trỗi) Trong điểm đà sử dụng mốc cũ là: Điểm trọng lực gốc Láng - Hà Nội Điểm trọng lực hạng I Đà Nẵng, điểm đợc đặt nhà làm việc Viện bảo tàng, nhà đà tôn cao nên mốc trọng lực bị thấp, tu sửa đà nâng thêm 15 cm Điểm trọng lực Nha Trang lúc đầu đặt Trờng Đại học Thủy sản nhng gần biển (cách biển 300 m), thời gian đo tháng 10 tháng 11 biển động mốc bị rung, không đo đợc nên phải chọn điểm mốc huyện Diên Khánh (cách biển khoảng 10 - 12 km) Điểm Diên Khánh đợc chọn chôn mốc Điểm trọng lực Điện Biên Phủ đợc chọn chôn mốc khu đất trạm DGPS Điện Biên Phủ Cục Đo đạc Bản đồ quản lý Điểm trọng lực TP Hồ Chí Minh sử dụng mốc điểm trọng lực gốc tầng nhà số 28 phố Nguyễn Văn Trỗi Phía Nga đà tiến hành kiểm nghiệm máy Liên bang Nga cuối tháng 10 năm 2005 đà đa máy trọng lực tuyệt đối GBL sang đo Việt Nam Nguyên tắc hoạt động máy trọng lực tuyệt đối GBL đo khoảng thời gian Ti mà vật rơi tự đợc khoảng cách cho trớc Si Đo khoảng cách Si mà vật rơi tự đợc nhờ máy giao thoa kế laser Kết đo máy trọng lực GBL nhận đợc số đọc khoảng thời gian Ti mà vật rơi tự đợc khoảng cách cho trớc Si Các giá trị Ti Si đợc tính theo thời ®iÓm L = 1, 2, N, sè ®äc N thay đổi từ 150 đến 600 phụ thuộc vào loại máy 24 Các giá trị lực trọng trờng gi đợc tính phơng pháp bình phơng tối thiểu để đạt đợc giá trị trọng lực gần tốt so với giá trị gia tốc trọng trờng điểm đo 7.2 Kết đo trọng lực tuyệt đối Trong tháng đà đo đợc giá trị trọng lực tuyệt đối điểm với độ xác nh sau: Điểm gốc trọng lực Láng - Hà Nội: = 0,0034 mGal Điểm Điện Biên Phủ: = 0,0035 mGal Điểm Đà Nẵng: = 0,0054 mGal Điểm TP Hồ Chí Minh (28 Nguyễn Văn Trỗi): = 0,0052 mGal Điểm Diên Khánh (Nha Trang): = 0,0187 mGal Nh vậy, có điểm trọng lực Láng (Hà Nội), Điện Biên Phủ, Đà Nẵng TP Hồ Chí Minh, độ xác gia tốc trọng trờng đạt àkGal Các điểm đạt tiêu chuẩn điểm trọng lực sở quốc gia Điểm Nha Trang đạt độ xác 19 àkGal, đạt độ xác điểm trọng lực hạng I Phơng án đo thèng nhÊt hƯ thèng ®é cao 8.1 Giíi thiƯu chung quốc gia hệ thống tọa độ xây dựng hệ thống độ cao quốc gia Hệ thống độ cao ba thành phần quan trọng trắc địa Trong báo cáo này, trình bày sở khoa học hệ thống độ cao phơng án thống hệ thống độ cao khác 8.2 Cơ sở khoa häc cđa hƯ thèng ®é cao Nh− ®· biÕt, ®Ĩ xây dựng hệ thống độ cao quốc gia cần xác định mặt khởi tính độ cao (mặt Geoid) Theo định nghĩa, mặt Geoid trọng lực mặt đẳng trọng trờng trái đất qua điểm độ cao (điểm khởi tính độ cao) trùng với mặt nớc biển trung bình Trong thực tế, mặt nớc biển trung bình yếu tố đại dơng nh dòng chảy, nhiệt độ, độ mặn áp suất nớc biển thay đổi mà chênh với mặt Geoid tới 1m Nh vậy, thực tế mặt Geoid trọng lực không trùng với mặt nớc biển trùng bình dùng làm mặt khởi tính độ cao Hai mặt chênh đại lợng h thay đổi theo thời gian Theo Hình vẽ 8.1 ta có: 25 A h C B mặt đất mặt n−íc biĨn trung b×nh D Geoid Ellipsoid Ĩ H×nh 8.1 đây: Đoạn AB = h: độ cao thủy chuẩn Đoạn CD = N: độ cao Geoid Đoạn AD = H: độ cao trắc địa Đoạn BC = h :độ cao chênh lệch Geoid mặt nớc biển trung bình nớc ta, trạm nghiệm triều Hòn Dấu đợc chọn trạm quan trắc để xác định mặt nớc biển trung bình làm mặt khởi tính độ cao cho hƯ thèng ®é cao ë n−íc ta Mèc ®é cao gốc Trạm thủy văn (Hòn Dấu) Geoid h Mặt nớc biển trung bình khu vực Mặt nớc biển tức thời Mặt nớc biển thủy triều mức cạn (trong khoảng thời gian T) Hình 8.2 Theo mô tả Hình 8.2 cho thấy đại lợng h thay đổi độ cao mặt nớc biển trung bình thay đổi Nh vậy, để xây dựng hệ thống độ cao quốc gia cần thiết phải xây dựng điểm gốc ®é cao qc gia (®iĨm khëi tÝnh ®é cao) trªn sở sử dụng số liệu quan trắc trạm thủy văn liên tục hàng chục năm độ cao mực nớc biển để xác định đợc độ cao trung bình mặt nớc biển khu vực quan trắc làm sở khởi tính độ cao cho hệ thống lới độ cao quốc gia Việc nghiên cứu xác 26 định độ chênh h mặt nớc biển trung bình độ mặt Geoid có ý nghĩa quan trọng việc nghiên cứu thống hệ thống độ cao khác việc xác lập mô hình Geoid độ xác cao Để xác định độ cao có phơng pháp: - Phơng pháp đo thủy chuẩn hình học; - Phơng pháp đo cao vệ tinh 8.3 Phơng pháp thống hệ độ cao quốc gia xây dựng cho hệ thống độ cao riêng Mỗi hệ thống độ cao có điểm khởi tính độ cao riêng (mặt nớc biển trung bình) mà độ cao cho (H0 = 0) Điều mặt nớc biển trung bình trùng với mặt Geoid điểm khởi tính độ cao Tuy nhiên, thực tế mặt nớc biển trung bình chênh với mặt Geoid đại lợng h Chính điều dẫn đến khác biệt hệ thống độ cao có điểm khởi tính độ cao khác Do độ chênh mặt nớc biển trung bình hi so với mặt Geoid điểm khởi tính độ cao đại lợng khác Ngày nay, phát triển hệ thống định vị vệ tinh GPS cho phép xác định tọa độ X, Y, Z độ cao trắc địa H Ellipsoid tham chiếu WGS-84 với độ xác cao Theo báo cáo Hội nghị lới trắc địa Châu - Thái Bình Dơng tổ chức Thành phố Hồ Chí Minh độ xác MX,Y,Z 1,0mm mH ≤ - cm Ngoµi ra, tr−êng träng lùc toàn cầu đà đợc xác định ngày xác mô hình trọng trờng EGM-96 cho phép xác định độ cao Geoid phạm vi toàn cầu với độ xác mN 1,0m vùng núi, vùng đồng mN 0,30 - 0,50m, vùng biển độ xác Geoid đạt 0,10 - 0,15m Tuy nhiên, quốc gia với số liệu trọng lực đo chi tiết mặt đất đủ mật độ độ xác sử dụng kết hợp với mô hình trọng trờng toàn cầu nh mô hình trọng lực EGM-96 có độ xác tới 0,10m, chÝ tíi 0,05m (Mü, Ph¸p ) ë n−íc ta, víi kết xây dựng mô hình GraGeoid việc chỉnh lý mô hình Geoid trọng lực điểm GPS - TC cho phép xác định độ xác độ cao Geiod vùng đồng 0,10m vùng núi 0,20 - 0,30m Với kết nh trên, để thèng nhÊt ®é cao khu vùc (hƯ thèng ®é cao nhiều nớc khác nhau) việc đo nối hệ thống độ cao nớc với (có đo trọng lực dọc tuyến độ cao) để thống hệ thống độ cao mặt quy chiếu, xây dựng lới độ cao thống khu vực công nghệ GPS Nội dung toán nh sau: Trên K điểm thủy chuẩn có độ cao hK hệ thống độ cao nớc khu vực tiến hành đo GPS xác định độ cao trắc địa HK Dựa vào mô hình 27 Geoid khu vực xác định đợc độ cao Geoid trọng lực NK Nh vậy, K điểm ta xác lập đợc hệ phơng trình có dạng: T H K hK − N K = a K X + V K (8.5) X: c¸c Èn sè VK: sè hiƯu chØnh (sai số ngẫu nhiên hệ quy chiếu) Thờng dïng biÓu thøc: T a K X = x + x1 cos ϕ K cos λ K + x cos ϕ K sin λ K + x3 sin ϕ K + x sin ϕ K C¸c ẩn: x0, x1, x2, x3, x4 cần tìm Lu ý r»ng biÓu thøc: HK - hK - NK ≠ nguồn sai số sau: - Không đồng hệ quy chiếu nguồn sai số hệ thống khác bao gồm: sai số xác định độ cao Geoid N, sai sè ®o thđy chn h, sai số khác mặt Geoid trọng lực mặt tham chiếu thủy chuẩn; - Sai số lấy giá trị gần để tính giá trị H N (giá trị trọng lực chuẩn gần đúng, độ cao địa hình gần đúng, giá trị gần ®óng cđa mËt ®é vËt chÊt, sù sai lƯch gi÷a mặt địa hình đáy biển mặt nớc biển trung bình trạm nghiệm triều, - Sai số biến động bề mặt tợng địa động học trái đất, v.v Nh vậy, điểm đo GPS ta có biểu thức: 0 H K − hK − N K = L0 K BiĨu thøc (8.5) cã thĨ viÕt d−íi d¹ng: a T X + L0 = V K K K (8.6) Víi ®iỊu kiƯn [Vα V T ] = tìm đợc nghiệm x0, x1, x2, x3, x4 Nh− vËy, ®é cao cđa mét ®iĨm bÊt kú khu vực xác định công thøc: T hK = H K − N K + AK X Trong đó: HK, NK: độ cao trắc địa độ cao Geoid điểm K X = (x , x1 , x , x3 , x ) : hệ số tính đợc từ giải phơng trình Kết tính toán tham sè thèng nhÊt hƯ thèng ®é cao n−íc Đông Dơng: Dữ liệu độ cao: Name NT02 19 HN00 VT02 φ (Vĩ tuyến) 10.38975 10.88831 09.30724 21.04534 10.33952 λ (kinh tuyến) 104.48365 105.78426 105.72037 105.78368 107.07191 H WGS-84 -10.155 -4.331 -1.870 -18.962 1.834 28 H TC 2.138 2.719 1.422 9.010 3.019 N EGM96 -12.021 -8.364 -4.671 -28.232 -1.882 Residuals -0.272 1.314 1.379 0.260 0.697 QN01 HN01 NT01 MN01 HOUA PHUT 42617 PAKSE VIENTIANE M.XAI 13.77498 18.80139 20.66847 21.43570 20.26185 16.61458 16.55053 15.12864 18.02551 20.69312 109.25436 105.77662 106.81424 107.96586 100.43721 106.50061 104.75513 105.78458 102.51559 101.99107 3.265 -20.394 -22.345 -20.696 373.287 258.794 117.896 87.556 192.537 607.636 3.697 4.092 2.541 2.348 409.483 275.730 141.470 104.428 223.560 640.764 -1.071 -23.906 -23.926 -21.873 -36.848 -17.474 -24.561 -18.309 -32.040 -33.898 0.639 -0.580 -0.960 -1.171 0.652 0.538 0.987 1.437 1.017 0.770 Các liệu đợc lấy từ kết đo GPS lới châu - Thái Bình Dơng Việt Nam, Lào Campuchia năm 1999 - 2000 số kết đo GPS độ xác cao việc xây dựng mô hình Geoid đồng sông Cửu Long Theo kết đánh giá thì: Độ xác độ cao trắc địa mH 5-6 cm Độ xác độ cao thuỷ chuẩn đạt từ hạng III trở lên Sơ đồ phân bố điểm: Theo sơ đồ hình 8.5 dới đây: (Kích thớc ô lới sơ đồ phân bố điểm 100 km x 100 km) Hình 8.5 29 - Kết tính toán: Nh vậy, độ cao điểm khu vực Đông Dơng tính c«ng thøc: hi = H i − N i + 47,877550 + 22,733979 cos ϕ i cos λ i − 46,412790 cos ϕ i sin λ i + 36,032440 sin ϕ i − 107,318615 sin ϕ i KÕt luận Kiến nghị 9.1 Kết luận Dự án Phối hợp nghiên cứu để xây dựng hoàn thiện hệ thống trọng lực phục vụ công tác nghiên cứu lnh thổ khai thác tài nguyên, bảo vệ môi trờng Việt Nam đà hoàn thành đợc nội dung sau: Nghiên cứu sở khoa học biến thiên trọng lực yêu cầu độ xác giá trị trọng lực cho việc xác định mô hình Geoid độ xác cao Việt Nam với kết nghiên cứu sau: ã Sự biến thiên trọng lực Việt Nam khoảng g 0,002 mGal/N ã Để xác định độ xác độ cao Geoid với độ xác cm mg= mGal, với độ xác cm mg= mGal cho ô chuẩn kích thớc 5km x 5km Thu thập, đánh giá tỉng quan hƯ thèng träng lùc hiƯn ë ViƯt Nam Lào Việt Nam, mật độ điểm trọng lực chi tiết chung nớc đạt từ 10 - 15 km2/điểm Tuy nhiên, vùng đồng trung du mật độ đạt từ - km2/điểm Vùng núi cao vùng biên giới mật độ đạt từ 20 - 30 km2/điểm, chí lớn Độ xác trung bình ô chuẩn 5km x 5km g 3,4 mGal Đối với Lào mật độ điểm trọng lực đạt 70 - 80 km2/điểm Nghiên cứu sở khoa học yêu cầu kỹ thuật cho việc xây dùng vµ hoµn thiƯn hƯ thèng träng lùc nhµ n−íc: ã Độ xác xác định tọa độ, độ cao điểm trọng lực tuyệt đối trọng lực hạng I: + Đối với trọng lực tuyệt đối: mX,Y 0,5 m mh ≤ 2,0 cm 30 + §èi víi träng lùc h¹ng I: mX,Y ≤ 2,0 m, mh ≤ cm ã Các yêu cầu kỹ thuật đo trọng lực tuyệt đối; ã Các yêu cầu kỹ thuật xây dựng lới trọng lực hạng I; ã Quy trình xử lý, tính toán ®¸nh gi¸ ®é chÝnh x¸c cđa ®o träng lùc tut đối trọng lực hạng I Xây dựng chơng trình tính toán, xử lý bình sai điểm trọng lực tuyệt đối lới trọng lực: ã Chơng trình tính giá trị trọng lực tuyệt đối; ã Chơng trình bình sai lới trọng lực Xây dựng phơng ¸n hoµn thiƯn hƯ thèng träng lùc nhµ n−íc ë Việt Nam: ã Xây dựng lới trọng lực tuyệt đối: 12 điểm ã Xây dựng lới trọng lực hạng I: 28 điểm Đo trọng lực tuyệt đối: Đo ®iĨm träng lùc tut ®èi víi ®é chÝnh x¸c mg 0,005 mGal, gồm điểm: ã Điện Biên Phủ ã Hà Nội ã Đà Nẵng ã TP Hồ Chí Minh Xây dựng phơng án xử lý hỗn hợp liệu mặt đất vệ tinh, đề xuất phơng án thống hệ thống lới độ cao khu vực phơng pháp GPS - TC - Geoid sở xây dựng lới GPS vào mốc độ cao quốc gia Xây dựng phơng án bình phơng tối thiểu Collocation cho việc đánh giá lới GPS - TC - Geoid phơc vơ cho viƯc x©y dùng mô hình Geoid độ xác cao Việt Nam 9.2 Kiến nghị - Hợp tác với Liên bang Nga ®Ĩ ®o tiÕp ®iĨm träng lùc tut ®èi víi ®é chÝnh x¸c 0,005 mGal b»ng c¸c m¸y träng lùc tuyệt đối GBL Liên bang Nga Đồng thời, hợp tác với Liên bang Nga để đo trọng lực hạng I máy AGAT (máy trọng lực tuyệt đối lắc) với độ xác mg 0,02 mGal - áp dụng phơng pháp bình phơng tối thiểu Collocation ®Ĩ ®¸nh gi¸ l−íi GPS - TC - Geoid phơc vụ cho việc xây dựng mô hình Geoid độ xác cao Việt Nam - Nghiên cứu áp dụng phơng pháp đo cao GPS mô hình Geoid độ xác cao để xác định độ chênh cao mặt nớc biển trung bình mặt Geoid điểm khởi tính độ cao (độ chênh mặt đẳng trọng lực (mặt Geoid) mặt 31 địa hình đáy biĨn) phơc vơ cho viƯc chØnh lý hƯ thèng ®é cao quốc gia - Tiếp tục hợp tác nghiên cứu với Liên bang Nga xử lý, tính toán bình sai lới trọng lực hạng I xử lý hỗn hợp liệu GPS - TC mô hình Geoid träng lùc viƯc x©y dùng thèng nhÊt hƯ thống độ cao xây dựng mô hình Geoid độ xác cao - Phối hợp với Liên bang Nga nghiên cứu xây dựng poligon địa động học nhằm nghiên cứu chuyển dịch vỏ trái đất phơng pháp vệ tinh - thủy chuẩn trọng lực 32 Tài liệu tham khảo 10 11 12 13 14 Quy ph¹m xây dựng lới trọng lực nhà nớc năm 1988 Liên Xô Xây dựng sở liệu trờng trọng lực toàn cầu thiết lập mô hình Geoid độ xác cao lÃnh thổ Việt Nam - GS TS KH Đặng Hùng Võ, TS Lê Minh, TS Phạm Hoàng Lân, KS Nguyễn Tuấn Anh Máy trọng lực tuyệt đối GBL 1989 Dự án: Xây dựng hoàn thiện hệ thống trọng lực nhà nớc - Viện Nghiên cứu Địa năm 2003 Height Determination by GPS A practical Experiment in central Sweden - Stag Gioran Martensson th¸ng 5/2001 Advanced physical Geodesy GMoris 1983 Space Geodesy and Model Geodinamika Moscow 1996 Gravimetry New York 1989 - Wolfgung Torge Gravity and Geoid Symposium 1994 - No 113 Reference Frames in Astronomy and Geophysics - Jean Kovalevsky, Ivan I Mueller th¸ng 7/1988 Evalution of Geoid Models and their use in combined GPS/Levelling/Geoid Height network Adjustment - Gfotopvulos, C.Kotsakis th¸ng 6/1989 Precision Geoid determination spherical FFT in and around the Korean peninsula - Hong Sic Yun th¸ng 10/1988 Gravimetry instruments Moscow 1999 Spectral Analysis of gravity field data and errors in view of sub-decimetre Geoid determination in Canada by Ohigbenga Esan th¸ng 5/2000 33 ... Đề tài Hợp tác Quốc tế Khoa học công nghệ Môi trờng Phối hợp nghiên cứu để xây dựng hoàn thiện hệ thống trọng lực phục vụ công tác nghiên cứu lnh thổ khai thác tài nguyên, bảo vệ môi trờng Việt. .. hệ thống trọng lực Việt Nam Lới đợc đo nối với hệ thống trọng lực quèc tÕ vµ khu vùc nh»m thèng nhÊt vµ chuÈn hoá hệ thống trọng lực Lới trọng lực sở quốc gia làm khống chế để phát triển hệ thống. .. lới trọng lực hạng thấp; Các điểm trọng lực vệ tinh dùng để lu giữ giá trị trọng lực độ xác cao điểm trọng lực sở trọng lực hạng I bị dùng để khôi phục lại điểm Để nghiên cứu biến thiên trọng