Đang tải... (xem toàn văn)
Công nghệ radar giao thoa đã được biết đến như một công nghệ để giám sát sự thay đổi các yếu tố trên bề mặt trái đất từ nhiều năm nay. Có khá nhiều phương pháp Radar giao thoa trong đó phương pháp giao thoa Radar tán xạ cố định (PSInSAR) sử dụng một sê ri ảnh để xác định biến ddoognj địa hình khá tốt. Tuy nhiên với khu vực có nhiều thực phủ thì số lượng điểm tán xạ cố định sẽ bị hạn chế, vì vậy trong bài báo này chúng tôi lựa chọn một phương pháp cải tiến từ phương pháp PSInSAR đó là sử dụng một sê ri ảnh đa thời gian với khoảng cách đường đáy ngắn (SBAS). Với chuỗi ảnh đa thời gian ALOS PalSAR 1 từ năm 2007 đến năm 2010 chúng tôi đã phát hiện được nhiều điểm trượt lở đất ở khu vực huyện Bát xát, Lào cai. Các vị trí được phát hiện từ ảnh radar khá phù hợp với các vị trí có nguy cơ trượt lở đất tại cùng khu vực do viện Địa chất và Khoáng sản cung cấp.
XÁC ĐỊNH TRƯỢT LỞ ĐẤT KHU VỰC LÀO CAI VIỆT NAM BẰNG PHƯƠNG PHÁP RADAR GIAO THOA ĐA THỜI GIAN Trần Vân Anh, Trần Hồng Hạnh, Lê Thanh Nghị Khoa Trắc địa Bản đồ Quản lý đất đai – Trường ĐH Mỏ - Địa chất Tóm tắt: Cơng nghệ radar giao thoa biết đến công nghệ để giám sát thay đổi yếu tố bề mặt trái đất từ nhiều năm Có nhiều phương pháp Radar giao thoa phương pháp giao thoa Radar tán xạ cố định (PSInSAR) sử dụng sê ri ảnh để xác định biến ddoognj địa hình tốt Tuy nhiên với khu vực có nhiều thực phủ số lượng điểm tán xạ cố định bị hạn chế, báo lựa chọn phương pháp cải tiến từ phương pháp PSInSAR sử dụng sê ri ảnh đa thời gian với khoảng cách đường đáy ngắn (SBAS) Với chuỗi ảnh đa thời gian ALOS PalSAR -1 từ năm 2007 đến năm 2010 phát nhiều điểm trượt lở đất khu vực huyện Bát xát, Lào cai Các vị trí phát từ ảnh radar phù hợp với vị trí có nguy trượt lở đất khu vực viện Địa chất Khoáng sản cung cấp Từ khóa: Trượt lở đất, SBAS, ALOS-PalSAR, Radar đa thời gian Đặt vấn đề Tai biến địa chất đại nói chung, tai biến trượt lở nói riêng gây nhiều tổn thất to lớn người cho nhiều quốc gia giới đặc biệt vùng miền núi Việt Nam Tai biến vấn đề thời sự, quan tâm đặc biệt nhà quản lý, nghiên cứu nhiều quốc gia giới Việt Nam số quốc gia nằm khu vực thường xuyên bị ảnh hưởng sạt lở đất khí hậu nhiệt đới ẩm gió mùa địa hình đồi núi nhiều Trong điều kiện biến đổi khí hậu tồn cầu diễn gay gắt sạt lở đất ngày trở nên nguy hiểm, khốc liệt Theo nghiên cứu Lê Quốc Hùng(Lê Q.H, nnk, 2015), năm nước ta vào mùa mưa xảy sạt lở đất, gây thiệt hại lớn người tài sản Ở Việt Nam có nhiều nghiên cứu sạt lở đất phải kể đến là: Nghiên cứu tai biến trượt lở điểm dân cư vùng thủy điện Hòa Bình (Bùi Khôi Hùng, 1992), hay dự báo tai biến thiên nhiên tỉnh Hòa Bình (Nguyễn Ngọc Thạch, 2003) Nghiên cứu nguy trượt lở miền núi Bắc Bộ giải pháp phòng tránh (Nguyễn Quốc Thành nnk, 2005) Ngồi nhiều nghiên cứu khác (Lê Thị Nghinh nnk, 2003), (Nguyễn Trọng Yêm nnk, 2004), (Trần Thanh Hà, 2004) Trong nghiên cứu có đề tài Lê Quốc Hùng (2012-2015) với đề án điều tra, phân vùng cảnh báo nguy trượt lở đất đá vùng miền núi nghiên cứu công phu đưa số đồ cảnh báo nguy trượt lở đất Các nghiên cứu chủ yếu xác định trượt lở đất cảnh báo trượt lở số phương pháp địa chất đồ phân bố trượt lở đất thường tạo điều tra trực tiếp thực địa kết hợp với giải đoán ảnh hàng không Từ năm đầu kỷ XXI, nhà khoa học giới sâu nghiên cứu vấn đề trượt - lở đất cho cơng bố nhiều cơng trình nghiên cứu có giá trị lĩnh vực (Lee S et al, 2001, 2002, 2005) Trong đó, nhiều cơng trình sử dụng tư liệu viễn thám vào việc xác định điểm trượt - lở đất, đới phá hủy kiến tạo, trạng lớp phủ thực vật, yếu tố có ảnh hưởng đến q trình trượt - lở đất Công nghệ viễn thám trở nên cơng cụ hữu ích xác định trượt lở đất chúng cung cấp nhìn tổng hợp lặp lại nhiều thời gian khác Đặc biệt, phép đo giao thoa vệ tinh SAR phương pháp đánh giá thay đổi bề mặt trái đất ứng dụng từ khoảng 20 năm (Anh T.V nnk, 2016) Differential SAR Interferometry (DInSAR) phương pháp sử dụng hai ảnh hai thời điểm khác vị trí trước sau có thay đổi địa hình xảy sử dụng để tìm dịch chuyển cách đo độ lẹch pha hai chu kỳ thu ảnh Tuy nhiên phương pháp có nhiều hạn chế khơng loại bỏ số ảnh hưởng: ảnh hưởng khí quyển, ảnh hưởng số đặc tính tán xạ đối tượng bề mặt Cách tiếp cận khác, PSI đề xuất Ferretti 2001 (Ferretti, A, 2001) , dựa việc sử dụng loạt ảnh SAR đa thời gian vị trí để chiết tách số điểm có tán xạ phản hồi cố định từ tìm biến động địa hình Các kỹ thuật PSI PSInSAR (Ferretti, A, 2001) từ dần thay phương pháp DInSAR Để tăng số lượng điểm tán xạ cố định số phương pháp khác phát triển SqueeSAR (Ferretti, A, 2011), Phương pháp STAMPS Stanford (Hopper, A nnk, 2010) cho xác định tán xạ cố định tán xạ phân bố, phương pháp đường đáy ngắn (SBAS) (Berardino nnk, 2002), công nghệ pixel tương quan (PSP, (Herrera, G., nnk , 2011), Kỹ thuật SAR giao thoa đa thời gian cho thấy khả cung cấp thông tin biến dạng mặt đất khu vực rộng có độ xác đạt tới milimet Phương pháp phù hợp với khảo sát khu vực có địa hình thực phủ nhiều độ dốc cao Đặc biệt, với kho lưu trữ ảnh SAR lớn tạo cho PSI khả đo giám sát tượng thay đổi địa hình khứ (Colesanti, C., nnk, 2003) Ngồi có số ảnh SAR cho miễn phí nguồn ảnh hữu ích cho nghiên cứu biến động Trong nghiên cứu sử dụng phương pháp Small Baseline (SBAS) dạng phương pháp SAR giao thoa đa thời gian cho xác định trượt lở đất khu vực Lào Cai Việt Nam với loại ảnh ALOS-PalSAR Lý chọn phương pháp với số lượng ảnh sử dụng khơng q nhiều bên cạnh phương pháp tối ưu hóa với việc sử dụng cặp ảnh có đường đáy ảnh ngắn Dữ liệu phương pháp nghiên cứu 2.1 Khu vực nghiên cứu Khu vực nghiên cứu huyện Bát xát thuộc tỉnh Lào cai (hình 1) Bát xát nằm phía tây bắc tỉnh Lào Cai, phía tây bắc đông bắc giáp tỉnh Vân Nam (Trung Quốc), phía tây giáp huyện Phong Thổ(Lai Châu), phía nam huyện Sa Pa thành phố Lào Cai, phía đơng nam thành phố Lào Cai Tồn địa hình Bát Xát kiến tạo nhiều dải núi cao, bật hai dải núi tạo nên hợp thuỷ: Ngòi Phát, suối Lũng Pơ, Suối Quang Kim Địa hình cao dần, điểm cao có độ cao 2945m, điểm thấp có độ cao 88m Kiến tạo địa hình Bát Xát hình thành hai khu vực Tuy nhiên, hai khu vực (vùng thấp gồm xã thị trấn, vùng cao gồm 16 xã) có chung đặc điểm: Vùng núi cao có độ chia cắt lớn, thung lũng hẹp khe sâu, độ dốc lớn Vùng thấp (ven sông Hồng, bồn địa nhỏ) nơi tập trung dải đồi thấp, thoải địa hình tương đối phẳng Hàng năm từ tháng đến tháng khoảng thời gian mưa nhiều thường xảy trượt lở đất lũ quét cho khu vực miền núi Việt Nam, huyện Bát xát điểm nóng trượt lở đất Theo nghiên cứu Lê Quốc Hùng (Lê Q.H, 2015) xã Phìn Ngan, xã Quang Kim hay quốc lộ 4D nơi thường xuyên bị trượt lở đất Hình 1: Bản đồ hành huyện Bát xát tỉnh Lào cai 2.2 Dữ liệu nghiên cứu Trong nghiên cứu này, sử dụng 13 cảnh ảnh ALOS PALSAR-1 với quỹ đạo lên, độ phân giải không gian 16m Bộ liệu sử dụng hình ảnh với đường quét 478 hàng quét 440 tạo thành tập hợp liệu khoảng thời gian ba năm, từ tháng năm 2007 đến tháng 11 năm 2010 Bảng thông tin tập hợp liệu ALOS PalSAR -1 Các liệu PALSAR-1 xử lý mức 1.1 chế độ phân cực đôi (FBD) xử lý để chọn kiểu phân cực HH Bên cạnh mơ hình số độ cao SRTM với độ phân giải 90 m sử dụng để loại bỏ pha địa hình hiệu chỉnh hệ tọa độ địa lý khu vực Việt Nam Bảng 1: Thông tin tập liệu STT 10 11 12 13 Ngày thu ảnh 2007/08/10 2007/09/25 2007/11/10 2008/05/12 2008/06/27 2008/08/12 2009/06/30 2009/08/15 2009/09/30 2010/07/03 2010/08/18 2010/10/03 2010/11/18 Mode FBD FBD FBD FBD FBD FBD FBD FBD FBD FBD FBD FBD FBD Orbit Ascending Ascending Ascending Ascending Ascending Ascending Ascending Ascending Ascending Ascending Ascending Ascending Ascending 2.3 Phương pháp nghiên cứu Việc sử dụng kỹ thuật giao thoa radar (InSAR) để đo bề mặt trái đất, bao gồm bề mặt địa hình biến dạng địa hình, chứng minh thành công hai thập kỷ qua Phương pháp giao thoa truyền thống cho phép tạo hình ảnh giao thoa pha dịch chuyển hai ba hình ảnh thu thời điểm khác khu vực Phương trình (1) minh họa pha biến dạng bề mặt đất đề xuất (Zebker nnk, 1986) Giả thiết có sẵn pha mơ hình số độ cao DEM, φTopo mô loại bỏ từ pha giao thoa tổng hợp ΔφInt ta thu gọi pha biến động địa hình (DInSAR) ΔφD-Int (Cơng thức 1) (1) Trong φTopo_simu thành phần địa hình mơ phỏng, có chứa thành phần pha phẳng Lưu ý quỹ đạo ảnh hưởng đến thành phần địa hình mơ phỏng, q trình làm phẳng khơng thực cách rõ ràng Cơng thức (1) tóm tắt ngun lý DInSAR, cho phép xác định dịch chuyển địa hình tạo từ hai ảnh SAR phức Cơng thức (1) đại diện cho phương trình xác định biến động địa hình từ DIsSAR đơn giản Phương trình đầy đủ biểu diễn (2): (2) Trong φTopo_res thành phần sai số tồn dư (RTE), φAtm thành phần pha khí thời điểm thu cảnh ảnh φOrb thành phần pha lỗi quỹ đạo cảnh ảnh φNoise pha nhiễu Thuật ngữ cuối cùng, 2k π, k giá trị số nguyên lần pha, kết chất pha bị đóng ΔφD-Int, nghĩa pha DInSAR bị giới hạn dãy (-π, π) Mục đích kỹ thuật DInSAR lấy φDispl từ ΔφD_Int Điều đồng nghĩa với việc tách φDispl từ thành phần pha khác phương trình số (2) Một điều kiện cần thiết để thực việc phân chia để phân tích điểm ảnh có nhiễu φNoise , thường liên quan đến hai tán phản xạ: nơi phản hồi radar vật thể phản chiếu mạnh không đổi theo thời gian gọi tán xạ cố định (Permanent Scatterer, PS ) phản hồi mà đáp ứng không đổi theo thời gian, đối tượng phân tán nhỏ khác (Distribatter Scatterers, DS) Những hạn chế DInSAR bao gồm: (i) tương quan thời gian hình học có ảnh hưởng đến thành phần φNoise (Hanssen, 1996); (ii) việc giải mở pha liên quan đến ước tính giá trị k (Ghiglia Pritt, 1998); (iii) thành phần khí (Zebker cộng sự, 1997) PSI đại diện cho phương pháp cơng nghệ SAR giao thoa, khai thác nhiều hình ảnh SAR thu khu vực, quy trình xử lý phân tích liệu thích hợp để tách riêng φDispl từ thành phần pha khác mô tả công thức (2) Phương pháp đường đáy ngắn (SBAS) ban đầu đề xuất (Berandino nnk.) SBAS sử dụng để phân tích mục tiêu phân tán Kết đầu giống với kết tạo phân tích giao thoa DInSAR, chúng có liên quan đến chuỗi thời gian SAR lớn thay liên quan đến hai (tối đa cảnh ảnh phương pháp DInSAR) Đối với trình xử lý PS, cách tiếp cận SBAS nhạy cảm với số lượng ảnh đầu vào, SBAS khai thác tương quan phân bố khơng gian thay dựa giá trị điểm ảnh; SBAS số giả định thực mối tương quan với điểm biến dạng địa hình Tuy nhiên trường hợp nhiều liệu đầu vào cho kết đầu tốt thành phần khí pha giao thoa ước tính loại bỏ tốt có nhiều thời kỳ ảnh Phương pháp nhằm mục đích giảm thiểu phân cách khoảng thời gian phạm vi tần số Doppler cặp ảnh thu nhận để tăng mối tương quan cặp giao thoa (Hopper cộng sự, 2010) phát triển phương pháp lọc pha tương quan thấp (SDFP) cho khoảng thời gian giãn cách ngắn Do đó, điểm ảnh SDFP bị lập bởi điểm ảnh liên quan xác định để xử lý Kết nghiên cứu thảo luận Bằng cách kết hợp cảnh ảnh thời điểm khác bảng trên, 12 cặp ảnh với đường đáy ngắn tạo Sự phân bố đường đáy ảnh vng góc thời gian biểu diễn Hình Việc xử lý giao thoa thực phần mềm SARscape với pha địa hình loại bỏ cách sử dụng SRTM DEM Tại khu vực Lào Cai, có hàng trăm khu vực không ổn định độ dốc với kích thước khác tốc độ sạt lở đất khác Chúng chọn huyện Bát Xát có đặc trưng địa hình địa chất điển hình: địa hình đa dạng, nhiều sạt lở đất tồn tốc độ khác nhau, phá hủy trực tiếp vào đường Sạt lở đất ảnh hưởng đến khu dân cư xung quanh Dền Sáng, Dền Thàng, Bản Xèo Sau tiến hành xử lý để quan sát thay đổi độ trượt vị trí phạm vi tập hợp ảnh ALOS PalSAR khu vực Bát xát, thấy trượt khu vực tiến triển nhanh, hình mơ tả phát triển trượt lở đất suốt thời gian nghiên cứu Lở đất lớn giai đoạn hiển thị vị trí hình ảnh đánh dấu hình elip hình Trong kết (Hình 3a) sau 46 ngày quan sát, thay đổi nhỏ, mầu xanh lam thể mức độ trượt lở 5cm (bên elip) Trong kết thứ hai (hình 3b) sau lần quan sát, trượt lở đất bắt đầu mạnh lần đầu tiên, tốc độ trượt tối đa khoảng 10 cm năm Sau năm quan sát (hình 4f), sụt lở đất xảy khơng vị trí đánh dấu mà nhiều vị trí khác Sáu số liệu Hình biểu loạt kết cho thấy vụ lở đất khác khu vực nghiên cứu Trượt lớn tổng cộng khoảng 150cm khu vực Lấy trung bình năm kết xử lý, tỷ lệ trượt lở đất trung bình lớn khoảng 50cm năm vị trí Dền Thắng, Dền Sáng, Bản Vược, Mường Vị Bản Xèo Hình 2: Thời gian đường đáy vng góc -50cm 50cm Hình 3: Sạt lở đất giai đoạn 2007-2010 (a): Thời gian 10/8/2007 25/9/2007, (b): 25/9/2007 10/11/2007; (c) : 5/12/2008 27 /6/2008, (d): 8/12/2008 30/6/2009, (e) 15/8/2009 3/7/2010, (f): 3/7/2010 11/11/ 2010 Hình 4: Sạt lở đất trung bình giai đoạn 2007-2010 ảnh Radar, điểm nguy trượt lở viện Địa chất Khoáng sản cung cấp Bản đồ địa điểm nguy trượt lở đất thực Viện Địa chất Khoáng sản vào năm 2015 (Lê Q.H, 2015) Những điểm đánh dấu hình tam giác mầu xanh (hình 4) vị trí điểm có nguy trượt lở chạy từ mơ hình dự báo nguy trượt lở đất khu vực Lào Cai Trong khu vực đánh dấu Dền Thắng, Dền Sáng khu vực sạt lở đất lớn với vị trí trượt phù hợp với điểm đánh dấu Mơ hình dự báo khơng khơng phải liệu đo theo chu kỳ nên xác định khối lượng trượt lở đất mà đánh dấu vị trí có nguy trượt lở cao, việc so sánh khó khăn Tuy nhiên, liệu đánh giá khả sử dụng ảnh vệ tinh việc xác định sạt lở đất Do thời kỳ xác định trượt lở đất ảnh Radar đồ xác định nguy trượt lở đất mơ hình khơng trùng thời gian với việc xác định trượt lở ảnh vệ tinh Radar xác định vị trí mà theo đồ đánh dấu nơi có khả trượt lở cao Việc sử dụng ảnh vệ tinh ALOS-PalSAR-1 với band ảnh L có bước sóng dài, có khả xuyên qua khu vực có thực phủ làm tăng khả xác định vị trí bị trượt lở đất mà bị che khuất tán cây, loại ảnh phù hợp với khu vực miền núi có thực phủ dày đặc Kết luận Việc sử dụng ảnh vệ tinh Radar đa thời gian giúp hiểu rõ mơ hình hóa tiến trình trượt lở đất Ngoài ra, dựa liệu ALSOS PalSAR liệu ảnh Radar khác, chúng tơi tính tốn vận tốc trượt lở đất thời điểm định, biến dạng nhỏ (cm / năm) Phương pháp SBAS khơng đòi hỏi nhiều ảnh phương pháp PSInSAR, thêm vào đó, phương pháp khắc phục số hạn chế PSInSAR, có điểm PS nhiều nơi che phủ thực vật Tuy nhiên, phương pháp có số hạn chế khả xác định độ dịch chuyển địa hình với liệu ALOS phụ thuộc nhiều vào lựa chọn quỹ đạo phù hợp với khu vực theo dõi Đối với khu vực che phủ thảm thực vật dày đặc, khơng có điểm tán xạ cố định (PS) theo quan sát thời điểm khác kích thước điểm PS khơng đủ lớn để phản ánh dịch chuyển hay thay đổi địa hình lúc phương pháp khơng thể thực Để khắc phục việc điểm tán xạ giả nên lắp đặt khu vực nghiên cứu, hoạt động điểm PS tín hiệu tán xạ điểm tốt cho xác định biến động địa hình Lời cảm ơn Nhóm tác giả xin gửi lời cảm ơn nhận tư liệu từ đề tài: "Nghiên cứu sở khoa học đề xuất giải pháp sử dụng ảnh Radar đa thời gian giám sát biến động lớp phủ bề mặt biến động địa hình phục vụ cho dự báo nguy biến động địa hình khu vực miền núi Việt Nam.", 2017-2019, Bộ Tài nguyên Môi trường Việt Nam Tài liệu tham khảo Anh, T.V., Cuong, T.Q., Anh, N.D., Khac, D.V., 2016, Study of Subsidence detection by DinSAR and evaluation of some factors to the outcome VIETNAM JOURNAL OF EARTH SCIENCES 37,344-354 Berardino, P.; Fornaro, G.; Lanari, R.; Sansosti, E., 2002, "A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms", IEEE Transactions on Geoscience and Remote Sensing, vol.40, no.11, pp 2375- 2383 Colesanti, C.; Wasowski, J., 2006, Investigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry Eng Geol., 88, 173–199 Colesanti, C.; Ferretti, A.; Prati, C.; Rocca, F., 2003, Monitoring landslides and tectonic motions with the Permanent Scatterers Technique.Eng Geol., 68, 3–14 Ferretti, A., Prati, C., Rocca, F.: 2001, Permanent scatterers in SAR interferometry IEEE Transactions on geoscience and remote sensing 39,8-20 Ferretti, A.; Fumagalli, A.; Novali, F.; Prati, C.; Rocca, F.; Rucci, A., 2011, A New Algorithm for Processing Interferometric Data-Stacks: SqueeSAR, IEEE Transactions on Geoscience and Remote Sensing, Vol.49, Issue: Ghiglia, D C and Pritt, M D., 1998, Two-Dimensional Phase Unwrapping, Theory, Algorithms, and Software, John Wiley & Sons, Inc., New York, USA, 493p Trần Thanh Hà, Đặng Văn Bào, Nguyễn Hiệu, Ngô Văn Liêm (2004), Ứng dụng phương pháp đánh giá đa tiêu chuẩn GIS nghiên cứu trượt lở đất tỉnh Lào Cai, Tuyển tập công trình khoa học, Hội nghị khoa học Địa lý – Địa chính, Đại học Quốc gia Hà Nội, 11-2004 Herrera, G.; Notti D.; Garcia-Davalillo, J.C.; Mora, O.; Cooksley, G.; Sanchez, M.; Arnaud, A.; Crosetto, M 2011, Analysis with C- and X-band satellite SAR data of the Portalet landslide area Landslides, 8, 195–206 10 Hooper, A., Spaans, K., Bekaert, D., Cuenca, M.C., Arıkan, M., Oyen, A., 2010, StaMPS/MTI manual Delft Institute of Earth Observation and Space Systems Delft University of Technology, Kluyverweg 1,2629 11 Hanssen, R F and Feijt, A., 1996, A first Quantitative Evaluation of Atmospheric Effects on SAR Interferometry, Proceedings of the ‘Fringe 96’ Workshop on ERS SAR Interferometry, Zurich, Switzerland, pp 277-282 12 Bùi Khôi Hùng, 1990, Nghiên cứu tình hình trượt lở điểm dân cư vùng hồ thủy điện Hòa Bình, Cơng ty khảo sát thiết kế điện 1, Hà Nội 13 Kampes, B.M.; Hanssen, R.F., 2004, Ambiguity resolution for permanent scatterer interferometry IEEE Trans.Geosci Remote Sens, 42, 2446–2453 14 Lanari, R.; Casu, F.; Manzo, M.; Zeni, G.; Berardino, P.; Manunta, M.; Pepe, A., 2007, An overview of the smallbaseline subset algorithm: A dinsar technique for surface deformation analysis Pure Appl Geophys,164, 637–661 15 Lee S, Min K, 2001: Statistical analysis of landslide susceptibility at Yongin, Korea Environmental Geology 40: 1095-1113 16 Lee S, Chwae U, Min K, 2002: Landslide susceptibility mapping by correlation between topography and geological structure: the Janghung area, Korea Geomorphology 46: 49-162 17 Lee S, Nguyen Tu Dan, 2005: Probabilistic landslide susceptibility mapping in the Lai Chau province of Vietnam: Focus on the relationship between tectonic fractures and landslides Environ Geol 48 778-787 SprigerVerlag 18 Le Quoc Hung, 2015, Investigation and delineation of landslides in Vietnam mountainous areas, Ministry of Natural Resources and Environment 19 Lanari, R.; Mora, O.; Manunta, M.; Mallorqui, J.J.; Berardino, P.; Sansosti, E., 2004, A small baseline approach for investigating deformation on full resolution differential SAR interferograms.IEEE Trans Geosci Remote Sens., 42, 1377–1386 20 Lê Thị Nghinh, 2003 Nghiên cứu đánh giá tai biến trượt lở đất khu vực tỉnh miền núi phía bắc giải pháp phòng tránh, Đề tài cấp Nhà nước 21 Mora, O.; Mallorqui, J.J.; Broquetas, A., 2006, Linear and nonlinear terrain deformation maps from a reduced set of interferometric SAR images.IEEE Trans Geosci Remote Sens, 41, 2243–2253 22 Nguyen Quoc Thanh, 2005, Initial results of the landslides and sliding disaster in mountainous areas in North Vietnam and recommendations for some preventive measures Scientific conference: Slides and floods - Flood mud, prevention solutions in the Northern Mountains (09/06/2005) Ministry of Science and technology, Scientific and technological program KC-08.Hanoi, p12-30 23 Nguyễn Trọng Yêm, 2006: Báo cáo Nghiên cứu đánh giá trượt - lở, lũ quét - lũ bùn đá số vùng nguy hiểm miền núi Bắc Bộ, kiến nghị giải pháp phòng tránh, giảm nhẹ thiệt hại Đề tài độc lập cấp nhà nước Lưu trữ Viện Địa chất, Viện KH&CN Việt Nam, Hà Nội 24 Nguyễn Ngọc Thạch nnk (2013), Kết thực nhiệm vụ hợp tác quốc tế khoa học công nghệ theo nghị định thư, Tăng cường lực nghiên cứu, đào tạo viễn thám Hệ thôngtin địa lý việc nghiên cứu, quản lý tai biến lũ lụt, lũ quét trượt lở đất, nghiên cứu điển hình Vĩnh Phúc Bắc Kạn, Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội 25 Perissin, D.; Wang, T., 2012, Repeat-pass SAR interferometry with partially coherent targets IEEE Trans Geosci Remote Sens., 50, 271–280 26 Zebker, H.A., Rosen, P.A and Hensley, S., 1997,Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps, Journal of Geophysical Research, 102(B4), pp.7547-7563 27 Zebker, H A and Goldstein, R M., 1986, Topographic mapping from interferometric synthetic aperture radar observations, Journal of Geophysical Research, 91(B5), pp 4993-4999 28 Paolo Berardino, Gianfranco Fornaro, Riccardo Lanari, and Eugenio Sansosti, 2002, A New Algorithm for Surface deformation Monitoring Based on Small Baseline Differential SAR Interferograms, IEEE transactions on Geoscience and Remote sensing, vol 40, no 11 DETERMINATION OF LANDSLIDE OF LAOCAI VIETNAM BY MULTI- TEMPORAL RADAR INTERFEROMETRY Tran Van Anh, Tran Hong Hanh, Le Thanh Nghi Geomatics and Land Administration – Hanoi University of Mining and Geology Abstract: Radar Interferometric technology has been known as a technology for monitoring land deformation for many years There are quite a number of methods of radar interferometry in which permanent scatter interferometric radar (PSInSAR) using a series of images to determine the land displacement very well However, in the case of dense vegetation, the number of permanent scattering PS will be limited, therefore in this paper we have chosen an improvement method from the PSInSAR, which is to use a series of images with short baseline (SBAS) With the ALOS PalSAR-1 multi-temporal images from 2007 to 2010 we have detected many landslides in Bat xat district, Lao Cai Positions detected from the radar image are well suited to landslide risk locations in the same areas provided by the Geological and Mineral Institute Keywords: Landslide, SBAS, ALOS-PalSAR, Multi-temporal Radar Thơng tin tác giả/nhóm tác giả Họ tên: Trần Vân Anh - Học hàm/họcvị/chức danh: PGS.TS - Cơ quan/Đơn vị: Khoa Trắc địa – Bản đồ-Quản lý đất đai, Trường Đại học Mỏ Địa chất - Địa chỉ:18 Phố Viên, Phường Đức Thắng, quận Bắc Từ Liêm - Số điện thoại liên lạc:0986811971 - E-mail:tranvananh@humg.edu.vn Họ tên: Trần Hồng Hạnh - Học hàm/họcvị/chức danh: TS, - Cơ quan/Đơn vị: Khoa Trắc địa – Bản đồ-Quản lý đất đai, Trường Đại học Mỏ Địa chất Lê Thanh Nghị - Học hàm/họcvị/chức danh: ThS - Cơ quan/Đơn vị: Khoa Trắc địa – Bản đồ-Quản lý đất đai, Trường Đại học Mỏ Địa chất ... sử dụng phương pháp Small Baseline (SBAS) dạng phương pháp SAR giao thoa đa thời gian cho xác định trượt lở đất khu vực Lào Cai Việt Nam với loại ảnh ALOS-PalSAR Lý chúng tơi chọn phương pháp với... định trượt lở đất ảnh Radar đồ xác định nguy trượt lở đất mơ hình khơng trùng thời gian với việc xác định trượt lở ảnh vệ tinh Radar xác định vị trí mà theo đồ đánh dấu nơi có khả trượt lở cao Việc... (hình 4) vị trí điểm có nguy trượt lở chạy từ mơ hình dự báo nguy trượt lở đất khu vực Lào Cai Trong khu vực đánh dấu Dền Thắng, Dền Sáng khu vực sạt lở đất lớn với vị trí trượt phù hợp với điểm đánh