tóm tắt luận án nghiên cứu cơ sở lý thuyết cho việc hiện đại hóa lưới khống chế trắc địa quốc gia ở việt nam bằng hệ thống vệ tinh dẫn đường toàn cầu gnss

28 792 0
tóm tắt luận án nghiên cứu cơ sở lý thuyết cho việc hiện đại hóa lưới khống chế trắc địa quốc gia ở việt nam bằng hệ thống vệ tinh dẫn đường toàn cầu gnss

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT NGHIÊN CỨ Ớ Ngành: Kỹ thuật Trắc địa - Bản đồ Mã số: 62520503 TÓM TẮT LUẬN ÁN TIẾN SỸ KỸ THUẬT – 2014 , Khoa Trắc địa, Trường Đại học Mỏ - Địa chất, Hà Nội Người hướng dẫn khoa học: Trường Đại học Mỏ - Địa chất Phản biện 1: TS Lê Minh Phản biện 2: - Luận án bảo vệ trước Hội đồng đánh giá luận án cấp Trường, họp Trường Đại học Mỏ - Địa chất vào hồi… giờ, ngày… tháng… năm 2014 Có thể tìm hiểu luận án thư viện: Thư viện Quốc Gia – Hà Nội Thư viện Trường Đại học Mỏ - Địa chất MỞ ĐẦU Tính cấp thiết đề tài Từ đời, hệ thống vệ tinh dẫn đường toàn cầu (GNSS) mở kỷ nguyên cho khoa học, cơng nghệ nói chung ngành trắc địa - đồ nói riêng Với ưu cơng nghệ, GNSS làm thay đổi cách quan niệm phương thức giải tốn xác định tọa độ vị trí điểm mặt đất : - Chính xác hóa, tăng dày điểm, hoàn thiện lưới tọa độ truyền thống; - Thiết lập lưới tọa độ quốc gia hoàn toàn tuân theo quan niệm dựa khả công nghệ GNSS Hiện mục đích, người ta thường gọi dạng lưới lưới GNSS CORS Tên gọi mang tính mơ tả cơng nghệ "quan trắc liên tục để kết nối điểm lưới hệ thống mạng trực tuyến" Đối với nước, mật độ điểm phân bố điểm khác Ở Việt Nam thời gian qua, công nghệ GNSS giúp cho việc nâ nhận, đánh dấu giai đoạn khoa học trắc địa đồ nước ta giới, việc tiếp tục xây dựng lưới tọa độ trắc địa túy công nghệ GNSS thay lưới trắc địa truyền thống nhằm thỏa mãn nhu cầu đại tọa độ cần đặt Ứng dụng cơng nghệ GNSS mang tính tồn cầu, thành phần hạ tầng thông tin địa lý xu hướng hình thành giới, phải hòa nhập mà khơng có lựa chọn khác cách có hệ thống vấn đề xây dựng lưới tọa độ trắc địa đại, đa mục tiêu nhằm bảo đảm nhu cầu đặt nay, phù hợp với phát triển thông tin địa lý giới Với tiêu đề “ ” Mục đích nghiên cứu luận án đất, hội nhập thông tin quốc tế,…) dựa công nghệ GNSS Đối tượng nghiên cứu Đối tượng nghiên cứu Nam ứng dụng công nghệ GNSS Phạm vi nghiên cứu Về khoa học: Khả cơng nghệ (tính toàn cầu, liên tục theo thời gian, phương pháp đo, độ xác, độ tin cậy, thành tựu kinh nghiệm quốc tế,…) GNSS Xử ITRF Các giải pháp, phương án xây dựng lưới tọa độ trắc địa nhằm thiết lập sở tọa độ thống trắc địa đồ, giải toán định vị hệ quy chiếu nhiệm vụ nghiên cứu khoa học Trái đất Về không gian: Lưới tọa độ trắc địa quốc gia phạm vi nước có kết nối với lưới tọa độ IGS quốc tế Về thời gian: Kể từ Việt Nam xây dựng lưới tọa độ trắc địa truyền thống Nội dung nghiên cứu : - Đánh giá tổng quan áp dụng công nghệ GNSS vào xây dựng lưới tọa độ trắc địa cấp độ quốc gia - Nghiên cứu khả công nghệ GNSS, phương pháp đo đặc biệt nhóm phương pháp đo động thời gian thực theo khu vực ), phục vụ mục đích quan trắc dịch chuyển vỏ Trái đất phục vụ mục đích xây dựng hệ quy chiếu quốc gia theo quan điểm động thời với số liệu IGS nhằm liên kết t Phương pháp nghiên cứu Luận văn sử dụng phương pháp nghiên cứu sau đây: , phân tích kế thừa: Thu thập : Thu thập số liệu - với lưới IGS, tính chuyển tọa độ, vận tốc hệ tọa độ động điển hình ITRF - Phương pháp so sánh: So sánh ưu điểm chuyên gia - Phương pháp ứng dụng công Ý nghĩa khoa học thực tiễn luận án - Ý nghĩa khoa học: Tiếp cận bước sâu phát triển công nghệ GNSS đưa r - Ý nghĩa thực tiễn: Kết nghiên cứu luận án có đóng góp định cho việc lự lưới tọa độ trắc địa, hệ quy chiếu tọa độ Việt Nam Các luận điểm bảo vệ luận điểm luận án a Các luận điểm bảo vệ - Luận điểm thứ nhất: Lưới tọa độ quốc gia Việt Nam theo quan điểm lư , hạ tầng thông tin, mạng kết nối,…) đáp ứng phương pháp đo đại; bảo đảm cho việc giải nhiệm vụ xác định bảo trì hệ quy chiếu tọa độ; đo đạc chi tiết phục vụ - Luận điểm thứ hai: Với cấu trúc mật độ lưới GNSS CORS, phương pháp thu thập xử lý số liệu phù hợp, lưới toạ độ c nhu cầu công tác trắc địa đồ, kết nối quốc tế thông tin địa lý, thực nghiên cứu khoa học Trái đất Việt Nam theo chuẩn toàn cầu - Luận đ ) b Các điểm luận án khoa học thực tiễn -Đ Kết cấu luận án , trình bày 132 trang, 25 hình vẽ sơ đồ, 34 bảng biểu Chương GNSS giới 1.2 Ứng dụng công nghệ GNSS vào xây dựng lưới trắc địa Việt Nam Từ tranh chung phát triển cơng nghệ GNSS giới q trình áp dụng cơng nghệ GNSS Việt Nam, nhận định phương hướng ứng dụng công nghệ GNSS vào Việt Nam Ngoài ưu điểm thấy rõ ràng công nghệ GNSS thời gian thi công nhanh, chi phí thi cơng thấp, khả cơng nghệ cao, công nghệ GNSS cung cấp điều kiện để tạo nên bước tiến quan trọng làm thay đổi quan niệm dạng thức lưới khống chế tọa độ quốc gia, cụ thể bao gồm: - Dạng thức lưới nhiều cấp với mật độ điểm cần chôn mốc dầy đặc theo cách thức truyền thống thay dạng thức lưới GNSS CORS sở với mật độ điểm cần chôn mốc thưa, tạo nên phương thức xác định số lượng phân bố điểm cuối cần chôn mốc phục vụ đo chi tiết địa hình, địa vật Đây dạng thức lưới khống chế tọa độ quốc gia, chưa đạt tới cách tiếp cận "trắc địa không lưới" đạt khả tiếp cận "lưới trắc địa mật độ thấp" - Lưới trắc địa GNSS động phụ thuộc thời gian cách tiếp cận hệ quy chiếu hệ tọa độ quốc gia, phù hợp với nhu cầu quan trắc Trái đất gắn với hoạt động kiến tạo Việc áp dụng cách tiếp cận vào quốc gia gắn với hệ quy chiếu quốc tế mảng vấn đề cần nghiên cứu chi tiết, cụ thể Chương NGHIÊN CỨU KHẢ 2.1 Những vấn đề chung thay đ ađ t - c, đ GNSS 2.1) GNSS gian SBAS GBAS GPS EGNOS EUPOS GLONASS WAAS SAPOS GALILEO MSAS … COMPASS … … … … [50] - Nguyên tắc l NRTK bắt đầu với tất trạm tham chiếu mạ tinh hệ thống phần mềm cài đặt hệ thống máy chủ NRTK trung tâm Mục đích NRTK để giảm thiểu ảnh hưởng sai số phụ thuộc vào khoảng cách vị trí tính tốn phạm vi giới hạn N - : + + ; ; + + ; : - , 2.1.2 Phương pháp xây dựng lưới GNSS CORS đóng vai lưới tọa độ trắc địa giới khu vực Đông Nam Á Các quốc gia thiết lập lư Nh giảm số lượng mật độ điểm; tự động hoá việc quy chiếu L Bảng 2.1 Số lượng điểm GNSS CORS Việt Nam TT 10 11 K/ GNSS CORS (km) 20 40 50 60 70 80 90 100 150 200 250 D/ (km2) 346 1386 2165 3118 4244 5543 7015 8660 19486 34641 54127 S.lượng trạm GNSS CORS 959 239 153 106 78 60 47 38 17 10 2.3 Xây dựng lưới GNSS CORS hoàn cảnh Việt Nam 2.3.1 Ý tưởng (tier) lưới GNSS CORS quốc gia, cấu trúc, mật độ điể độ xác 2.2 12 CORS Chương NGHIÊN CỨU CHUYỂN DỊCH HIỆN ĐẠI VỎ TRÁI ĐẤT Tại Việt Nam, số lưới quan trắc công nghệ GNSS thực tạo nên sở liệu ngày phong phú Những kết nghiên cứu góp p Bảng 3.1 Tổng hợp số liệu đo TT Ngày DOY 1/12/2010 2/12/2010 3/12/2010 4/12/2010 1/6/2011 2/6/2011 335 336 337 338 152 153 Tên điểm MCRS DNRS PQRS x x x x x x x x x x x x x x x x x x TSRS KUNM NTUS PIMO WUHN x x x x x x x x x x x x x x x x x x x x x x x x x x 13 10 11 12 13 14 3/6/2011 154 x x x x x x x x 4/6/2011 155 x x x x x x x x 5/6/2011 156 x x x x x x x x 1/3/2012 61 x x x x x x x 2/3/2012 62 x x x x x x x 3/3/2012 63 x x x x x x x 4/3/2012 64 x x x x x x x 5/3/2012 65 x x x x x x x : DOY - Day of year, ngày năm 3.4.2 Khai ố liệu liệu hỗ trợ quốc tế Tất số liệu hỗ trợ tệp số liệu vệ tinh, số liệu trục quay Trái đất, mơ hình điện ly tồn cầu, khai thác từ IGS 3.4.3 Phần mềm xử lý phần mềm Bernese 5.0 với quy trình chặt chẽ theo trình đơn 3.4.4 Các bước thực Bảng 3.2 Tọa độ điểm TT X (m) DNRS MCRS PQRS TSRS Bảng 3.3 Z (m) -1915754.4967 -1831172.4798 -1521499.8762 -2354092.2732 5824431.0548 5646418.0114 6094084.7056 5850302.4095 1750989.0477 2325703.8077 1104443.9268 952491.9085 /CORS VB (mm/n) TT Y (m) DNRS MCRS PQRS TSRS VL (mm/n) VH (mm/n) V (mm/n) (0 ) 34.7 30.8 21.1 20.2 -20.5 -7.8 -12.2 -13.5 44.9 -1.8 11.7 -6.0 40.3 31.8 24.4 24.3 120.6 104.2 120.0 123.8 - Như 14 1040 1080 1160 1120 TRUNG Q (MT) (ĐB) MCRS 200 (TB-BTB) DNRS 160 ) CAMPU CHIA ) 120 (NTB-NB) ) PQRS TSRS 80 MALAYSIA 30 mm/năm 100 km MALAYSIA Hình 3.1 DGPS/CORS - Quá trình xử lý tiến hành theo quy trình nghiêm ngặt - Việc tích hợp phương trình chuẩn theo ca đo xác định vận tốc điểm: MCRS (Móng Cái), DNRS (Đà Nẵng), PQRS (Phú Quốc) TSRS (Trường Sa) với giá trị vận tốc cỡ khoảng cm/năm Sai số vận tốc theo chiều thẳng đứng xấp xỉ mm/năm; theo chiều vĩ tuyến kinh tuyến mm/năm So sánh sai số với giá trị vận tốc cỡ cm kết có độ xác cao Chuyển dịch đại vỏ Trái đất tương đươ 15 - Đông Nam, phù hợp với kết khác công bố hướng độ lớn - Bên cạnh giá trị vận tốc, kết đưa tọa độ điểm hệ quy chiếu quốc tế ITRF05 liên quan tới thời điểm 2010 - 03 - 04 00:00:00 Phần mềm cho phép xác định tọa độ liên quan tới thời điểm lựa chọn - Các véc tơ vận tốc nhận qua việc xử lý thực Chương XÂY DỰNG HỆ QUY CHIẾU TỌA ĐỘ QUỐC GIA THEO QUAN ĐIỂM HỆ QUY CHIẾU ĐỘNG 4.1 Quá trình hình thành hệ quy chiếu trắc địa động giới quốc tế (ITRF) với tham số gồm tham số gốc địa tâm, tham số góc định hướng trục tọa độ tham số tỷ lệ chiều dài Mỗi năm, tham số xác định cụ thể với giá trị thực tế tương ứng với năm Do tham số trạng thái động (trái đất tồn trạng thái động) nên người ta bổ sung thêm tham số tốc độ thay đổi tham số này, đưa tổng số lượng tham số thành 14 4.2.1 Khái quát hệ quy chiếu tọa độ quốc gia VN - 2000 Hệ quy chiếu hệ toạ độ quốc gia VN - 2000 h nên 16 - ,n 4.1 - BTB 10 - 8.5 17 25.5 34 42.5 51 59.5 68 76.5 85 - NB 15.9 31.8 47.7 63.6 79.5 95.4 111.3 127.2 143.1 159 : (mm) dmax - 16 32 48 64 80 96 112 128 144 160 - , tọa độ biến dạng vỏ Trái đất 16 32 48 64 80 96 112 128 144 160 (NTB - NB) bị ảnh hưởng đáng kể kiến tạo ,d 17 100 mm 4.1) c max năm S , 4.4.1 Hệ thống công thức chuyển đổi (VB , VL , VH ) ITRF(1) (t - t0) 14 ITRF(t0)(1) ITRF(t)(2) (VX , VY , VZ ) ITRF(1) X(t0)(1) X(t)(1) X(t)(1) 4.1 sang ITRF(2) X(t)(2) a đ t ITRF(1) iđ mt i đ m t0 18 (t - t0) 14 (VB , VL , VH ) ITRF(1) ITRF(t0)(1) ITRF(t)(2) (VX , VY , VZ ) ITRF( ) X(t0)(1) X(t)(1) (VX , VY , VZ ) ITRF( ) (VB , VL , VH ) ITRF( ) 4.2 c t ITRF(1) sang ITRF(2) 4.4 Bảng 4.2 18/7/2012 TT Tên điểm X (m) Y (m) Z (m) XUYO -1606525.3107 5700807.4794 2358910.3445 OAN0 -1566459.7943 5711793.1162 2359391.0194 DOI0 -1554911.1281 5722287.9981 2341248.9526 SOC1 -1621267.4901 5719456.3583 2303158.3259 TAM2 -1601083.2914 5719707.1284 2318672.4007 HUN1 -1571100.2652 5731225.1349 2308717.3752 NTHO -1555548.8832 5730756.5177 2320416.1883 SON1 -1557982.9310 5741852.5798 2291130.9723 LAP1 -1541103.3835 5738426.8455 2311048.3542 10 BAVI -1578207.2652 5740266.8121 2281398.4277 19 11 HOA1 -1561203.3469 5754853.9430 2257478.9710 12 NAM0 -1480436.1083 5741958.4656 2342822.2803 13 NGA1 -1353072.7825 5749668.0049 2402462.2895 14 HAM1 -1355321.3541 5762376.1094 2367436.5756 15 DON1 -1334851.1103 5758491.9047 2387972.1199 16 LEM1 -1335854.0006 5773067.2098 2353215.1336 17 MON1 -1464209.0062 5767238.2661 2291232.1061 18 NOI1 -1457371.9564 5771322.4579 2285116.1053 19 LOT1 -1445752.0185 5771190.8986 2292523.2300 20 QTA2 -1432563.6411 5770244.6106 2303209.1292 21 NAD2 -1458337.9705 5777358.5267 2270039.5502 22 DIEB -1336842.3829 5787988.4739 2315702.2299 23 DOSN -1724757.3380 5714523.9160 2239792.0393 24 NT01 -1726969.5917 5714864.9568 2237081.3908 25 NT03 -1844373.6094 5997105.5965 1142317.0439 26 NT04 -1575936.5479 6075089.2127 1132070.0742 27 QT01 -1339440.8883 5788398.0139 2313170.2549 28 QT03 -1916791.4419 5822974.9368 1754668.6848 29 VUNT -1849617.0276 5995299.9136 1143372.7201 30 CAMP -1772774.5625 5687232.9039 2271331.8300 31 NONN -1921866.9957 5823665.8416 1747139.7398 32 DNRS -1915754.6080 5824431.1345 1750989.0334 33 MCRS -1831172.5471 5646417.9871 2325703.7915 34 PQRS -1521499.9300 6094084.7200 1104443.9070 35 TSRS -2354092.3105 5850302.3783 952491.8786 20 Bảng 4.3 Vận tốc chuyển dịch tuyệt đối điểm GNSS ITRF08 TT Tên điểm VL (mm/n) VB (mm/n) VH (mm/n) V (mm/n) P.vị (0 ) XUYO -10.8 35.1 0.6 36.7 107.1 OAN0 -10.2 33.6 -13.9 35.1 106.9 DOI0 -10.1 33.8 -10.4 35.3 106.6 SOC1 -10.2 33.0 -1.7 34.6 107.2 TAM2 -10.5 32.5 -0.5 34.2 107.9 HUN1 -10.0 33.2 -5.1 34.7 106.8 NTHO -10.5 33.4 -10.6 35.0 107.4 SON1 -10.4 32.9 -3.6 34.5 107.5 LAP1 -10.9 33.0 -4.7 34.8 108.3 10 BAVI -9.4 32.2 -1.2 33.6 106.3 11 HOA1 -9.6 34.0 -2.7 35.3 105.8 12 NAM0 -10.3 35.2 3.3 36.7 106.3 13 NGA1 -8.0 39.2 12.2 40.0 101.5 14 HAM1 -8.8 33.0 15.4 34.2 104.9 15 DON1 -10.4 35.3 15.3 36.8 106.4 16 LEM1 -9.4 34.9 12.0 36.2 105.1 17 MON1 -11.7 32.5 -1.0 34.6 109.8 18 NOI1 -10.4 33.1 -4.3 34.7 107.4 19 LOT1 -11.9 33.6 -6.1 35.7 109.5 20 QTA2 -10.8 34.0 -9.5 35.7 107.6 21 NAD2 -10.4 32.5 -10.1 34.1 107.7 22 DIEB -9.6 27.3 4.0 28.9 109.4 23 DOSN -8.0 27.6 16.7 28.8 106.2 24 NT01 -11.7 36.3 7.6 38.1 107.9 25 NT03 -11.7 27.1 18.8 29.5 113.4 21 26 NT04 -9.8 15.9 -14.5 18.6 121.6 27 QT01 -10.3 30.2 -19.9 31.9 108.8 28 QT03 -15.0 26.9 -2.6 30.8 119.1 29 VUNT -11.1 22.8 1.0 25.4 116.0 30 BDE1 -14.1 22.6 0.5 26.6 122.0 31 DTE1 -11.6 24.1 3.1 26.8 115.7 32 BLO1 -14.1 24.1 2.6 27.9 120.3 33 ALO1 -15.4 29.5 19.5 33.2 117.6 34 LNI1 -13.5 26.1 3.9 29.4 117.3 35 THN1 -13.2 25.5 2.7 28.7 117.4 36 BLV1 -14.7 27.7 -0.3 31.4 118.0 37 LANG -12.5 38.3 -22.6 40.3 108.1 38 DOHO -9.1 24.9 -15.1 26.5 110.1 39 HUES -19.8 30.0 10.2 35.9 123.4 40 STT1 -10.3 21.9 -7.5 24.2 115.2 41 CDA1 -7.4 20.7 -9.6 22.0 109.7 42 HOCM -11.5 22.0 -3.0 24.8 117.6 43 CAMP -4.4 47.7 28.1 47.9 95.3 44 NONN -0.2 44.5 24.6 44.5 90.3 45 DNRS -20.5 35.0 45.0 40.6 120.4 46 MCRS -7.8 31.0 -1.8 32.0 104.1 47 PQRS -12.1 21.4 11.8 24.6 119.5 48 TSRS -13.5 20.4 -5.9 24.5 123.5 * phạm vi lãnh thổ Việt Nam 22 4.3 Phương pháp chuyển đổ áp dụng tính chuyển tọa độ khung quy chiếu quốc gia mang tính động (ví dụ gắn với ITRF) khai hệ tọa độ động tương lai nước ta 1040 1080 1160 1120 200 ) 160 CAMPU CHIA 120 ) 80 30 mm/năm 100 km 4.1 MALAYSIA 23 : v - - 24 K : , phương pháp đo thiết lập phát triển không ngừng Hệ thống tọa độ khống chế dạng lưới điểm GNSS CORS xu tất yếu Bởi vậy, hệ thống hỗ trợ (mạng trạm đo liên tục, trung tâm xử lý, sở hạ tầng thông tin,…) cần thiết lập sớm tốt Với sở hạ tầng này, nhiệm vụ định vị giải công nghệ đại, đáp ứng yêu cầu độ xác, mang lại hiệu kinh tế cao, thúc đẩy q trình hịa nhập tồn diện với giới (2) Tăng cường việc nghiên cứu ứng dụng đa hệ GNSS: Một đặc điểm quan trọng lưới GNSS CORS hoạt động ổn định, liên tục, tin cậy với mức độ tự động hóa tin học hóa cao Về ứng dụng hệ GPS, hệ GLONASS trở lại chế độ hoạt động toàn bộ, hệ GALILEO phóng vệ tinh đầu tiên, phát thu tín hiệu dạng demo Bởi vậy, cần xác định cách tiếp cận có bước chuẩn bị (máy thu, phần mềm, quyền sử dụng,…) để ứng dụng đa hệ, nâng cao hệ số an toàn, tin cậy cho tồn hệ thống, chủ động đối phó với tình khơng thuận lợi DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ (2009), m Global Mapper" - , tr 94-101 (2009), "Nghiên cứu áp dụng phương pháp bình sai lặp đ " 1, tr 37-41 (2010), trạm tham chiếu quan trắc liên tục CORS", 14 (100), tr 52-53 (2011), "Các hệ thống định vị vệ tinh tương lai máy thu GNSS" 15 (125), tr 18-20 (2011), "Sử dụng phần mềm để chuyển đổi liệu dạng Raster hai hệ tọa độ vng góc khơng gian" 15 (125), tr 27-29 (2011), "Khung hệ quy chiếu Trái đất", Tạp c 9, tr 43-47, 53 (2012), " 21 (155), tr 30-32 (2012), "Nghiên cứu áp dụng phương pháp bình sai nhiều giai đoạn việc xử lý số liệu lưới khống chế trắc địa mặt độ cao" ,T (2012), "Tính chuyển vận tốc chuyển dịch tuyệt đối khung quy chiếu Trái đất quốc tế (ITRF)" - Địa chất (2012), " (2013), khung quy chiếu Trái đất quốc tế" - 41, tr 52-56 (2013), " 2013, tr 118-123 (2013), " ăm 2013, tr 223-228 , (DGPS)" ) ... hệ thống vệ tinh dẫn đường toàn cầu (GNSS) mở kỷ nguyên cho khoa học, cơng nghệ nói chung ngành trắc địa - đồ nói riêng Với ưu công nghệ, GNSS làm thay đổi cách quan niệm phương thức giải toán... không gian: Lưới tọa độ trắc địa quốc gia phạm vi nước có kết nối với lưới tọa độ IGS quốc tế Về thời gian: Kể từ Việt Nam xây dựng lưới tọa độ trắc địa truyền thống Nội dung nghiên cứu : - Đánh... công nghệ GNSS giúp cho việc nâ nhận, đánh dấu giai đoạn khoa học trắc địa đồ nước ta giới, việc tiếp tục xây dựng lưới tọa độ trắc địa túy công nghệ GNSS thay lưới trắc địa truyền thống nhằm

Ngày đăng: 03/10/2014, 11:02

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan