0273 ứng dụng viễn thám đánh giá nguy cơ hạn hán khu vực huyện bắc bình tỉnh bình thuận

Trịnh Lê Hùng tgk TẠP CHÍ KHOA HỌC ĐHSP TPHCM _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ ỨNG DỤNG VIỄN THÁM ĐÁNH GIÁ NGUY CƠ HẠN HÁN KHU VỰC HUYỆN BẮC BÌNH, TỈNH BÌNH THUẬN TRỊNH LÊ HÙNG*, ĐÀO KHÁNH HỒI* TĨM TẮT Cơng nghệ viễn thám với ưu điểm bật so với phương pháp nghiên cứu truyền thống ứng dụng hiệu nghiên cứu, giám sát ứng phó với tượng hạn hán Bài báo trình bày kết đánh giá nguy hạn hán khu vực huyện Bắc Bình (tỉnh Bình Thuận) từ tư liệu ảnh vệ tinh đa phổ LANDSAT sử dụng số khô hạn nhiệt độ thực vật (TVDI) Kết nhận sử dụng thành lập đồ nguy khô hạn giảm thiểu thiệt hại hạn hán gây Từ khóa: hạn hán, viễn thám, nhiệt độ bề mặt, độ ẩm đất, số khác biệt thực vật, ảnh đa phổ LANDSAT ABSTRACT Drought risk evaluation using remote sensing: a case study in Bac Binh district, Binh Thuan province Compared to traditional methods, remote sensing technology with advantages has been used effectively in studying, monitoring and reacting to drought This article presents results of drought risk evaluation from LANDSAT multispectral images in BAC BINH district (BINH THUAN province) using temperature vegetation dryness index (TVDI) The results obtained in this study can be used to create the drought risk map and to minimize the damage caused by drought Keywords: drought, remote sensing, land surface temperature, soil moisture, normalized difference vegetation index, LANDSAT multispectral image Mở đầu Hạn hán tượng tự nhiên gây ảnh hưởng nghiêm trọng đến môi trường sống hoạt động sản xuất người dân Hạn hán đánh giá thiên tai gây thiệt hại nặng nề đứng thứ ba sau lũ, bão có xu hướng xảy gay gắt, khó kiểm sốt tác động biến đổi khí hậu Ở Việt Nam, hạn hán xảy hầu khắp nước với mức độ thời gian khác nhau, đặc biệt nghiêm trọng khu vực miền Trung Tây nguyên, gây thiệt hại to lớn kinh tế – xã hội, đặc biệt nguồn nước sản xuất nông nghiệp Hạn hán thường xảy diện rộng, việc quan trắc nghiên cứu phương pháp truyền thống gặp nhiều khó khăn, thực tế đặt trạm quan trắc với mật độ dày đặc chi phí lớn Dữ liệu viễn thám cung cấp thông tin * TS, Học viện Kĩ thuật Quân sự, Hà Nội; Email: trinhlehung125@gmail.com bề mặt Trái Đất kênh phổ khác độ phủ trùm rộng sử dụng hiệu quan trắc giám sát hạn hán Đã có nhiều nghiên cứu giới sử dụng tư liệu viễn thám hồng ngoại nhiệt xác định nhiệt độ độ ẩm đất nhằm đánh giá mức độ khô hạn bề mặt [6, 9, 11, 14] Ở Việt Nam, số nghiên cứu sử dụng ảnh nhiệt MODIS, NOAA/AVHRR xác định độ ẩm đất dựa mối quan hệ nhiệt độ bề mặt loại hình lớp phủ [4] Tuy nhiên, độ phân giải không gian ảnh MODIS, NOAA/AVHRR thấp khơng thích hợp cho nghiên cứu chi tiết Bài báo trình bày kết đánh giá nguy hạn hán khu vực huyện Bắc Bình (Bình Thuận) sử dung tư liệu ảnh hồng ngoại nhiệt LANDSAT Ảnh hồng ngoại nhiệt LANDSAT với độ phân giải không gian trung bình (60 – 120m) cung cấp thơng tin rõ ràng thay đổi độ ẩm bề mặt so với ảnh MODIS, NOAA/AVHRR, sử dụng hiệu nghiên cứu giám sát tượng hạn hán Phương pháp nghiên cứu Các nghiên cứu giới cho thấy, nhiệt độ bề mặt lớp phủ thực vật yếu tố quan trọng cung cấp thông tin độ ẩm bề mặt đất [9, 12, 13, 14] Nhiệt độ tăng lên nhanh trường hợp bề mặt khô hạn trồng bị thiếu nước [12, 13] Để đánh giá mức độ khô hạn bề mặt, nghiên cứu sử dụng số khô hạn nhiệt độ thực vật (temperature vegetation dryness index – TVDI), Saldholt I đưa năm 2002 sở lượng hóa mối quan hệ nhiệt độ bề mặt lớp phủ [13] Chỉ số khô hạn nhiệt độ thực vật TVDI xác định theo công thức sau: 𝑇𝑉𝐷𝐼 = 𝑇�−𝑇�𝑚i� 𝑇�𝑚𝑎𝑥−𝑇�𝑚i (1) � đó: 𝑇�– nhiệt độ bề mặt, 𝑇�𝑚i�, 𝑇�𝑚𝑎𝑥 tương ứng nhiệt độ bề mặt cực tiểu cực đại tam giác không gian nhiệt độ/NDVI 𝑇�𝑚i�, 𝑇�𝑚𝑎𝑥 xác định phương pháp hồi quy tuyến tính giá trị nhiệt độ cực đại khoảng giá trị NDVI Giá trị số TVDI cao tương ứng với nguy khô hạn tăng Tại cạnh khô, số TVDI có giá trị 1, cạnh ướt giá trị TVDI [13] (hình 1) Hình Tam giác khơng gian “nhiệt độ/NDVI” TẠP CHÍ KHOA HỌC ĐHSP TPHCM Trịnh Lê Hùng tgk _ _ _ _ _ _ _ _ _ _ _ _ _ _ N ộ i d u n g n g h i ê n c ứ u v k ế t q u ả t h ự c n g h i ệ m K hu vực thực nghiệ m lựa chọn TẠP CHÍ KHOA HỌC ĐHSP TPHCM Số 5(70) năm 2015 nghiên cứu huyện Bắc số MI – Meteorology index) [5] Trong b) Bình (Bình Thuận) Bắc năm gần đây, biến động Bình hai khó lường khí hậu huyện tỉnh Bình Thuận tác động tiêu cực từ hoạt động có nguy suy thối đất người làm cho tình trạng hạn hán hoang mạc hóa hàng tỉnh Nam Trung Bộ, Tây Nguyên nói đầu khu vực duyên hải chung Bình Thuận nói riêng diễn Nam Trung Bộ, ngày nghiêm trọng Hạn hán diễn nguyên nhân thường xuyên hơn, vào mùa hạn hán gây Với đặc khơ mùa mưa, ảnh điểm khí hậu địa hình tự hưởng lớn đến hoạt động sản xuất nhiên làm cho Bắc Bình đời sống sinh hoạt người dân Do khô nóng quanh năm, hình vậy, ứng dụng tư liệu viễn thám phục vụ thành nên chế độ khí hậu cơng tác giám sát ứng phó với hạn hán bán khơ hạn trở thành vấn đề có ý nghĩa thực tiễn vùng khô Tư liệu viễn thám sử dụng hạn nước Theo số nghiên cứu ảnh vệ tinh quang học liệu thống kê trạm Bàu Landsat chụp vào 05 – 01 – 2002, 13 – 01 Trắng, giai đoạn – 2005, 12 – 02 – 2010 (ảnh Landsat 1960 – 2010, lượng mưa ETM+) 15 trung bình năm – 02 – 2014 (ảnh Landsat OLI) (hình Bắc Bình từ 700 – 1000 1) Ảnh chụp vào mùa khô, thời mm, tập trung chủ gian cao điểm hạn hán khu vực yếu vào mùa mưa (từ tháng miền Trung, Tây Nguyên Với độ phân đến tháng 10), giải khơng gian trung bình (30 m vào mùa khô (từ tháng 11 đến tháng năm sau), lượng mưa [3, 5] Cùng với Tuy Phong, Bắc Bình huyện dự báo có cấp độ khơ hạn nghiêm trọng Bình Thuận (cấp theo kênh đa phổ, 60 – 120m kênh hồng ngoại nhiệt, 15m kênh toàn sắc), đặc biệt cung cấp hoàn toàn miễn phí với chu kì cập nhật 16 ngày, ảnh Landsat nguồn tư liệu quý giá phục vụ nghiên cứu tài nguyên thiên nhiên giám sát môi trường [17] c) d) Hình Ảnh Landsat khu vực Bắc Bình (Bình Thuận) ngày 05 – 01 – 2002 (a), 13 – 01 – 2005 (b), 12 – 02 – 2010 (c) 15 – 02 – 2014 (d) tổ hợp màu RGB = NIR:RED:GREEN Để tính nhiệt độ bề mặt, giá trị số nguyên ảnh Landsat chuyển đổi giá trị thực xạ điện từ (W𝑚−2μ𝑚−1) Đối với tư liệu ảnh Landsat ETM+, việc chuyển đổi thực theo công thức sau [10]: Lλ = L max− L (DN − DN min)+ L DN max− DN min (2) đó, Lλ - giá trị xạ phổ (spectral radiance);Lmax, Lmin – giá trị xạ phổ ứng với DNmax DNmin kênh hồng ngoại nhiệt, cung cấp file metadata liệu ảnh LANDSAT; DNmax, DNmin – giá trị số lớn (=255) nhỏ (=1) [10] Với ảnh LANDSAT OLI, giá trị xạ xác định sau [10]: Lλ = M L Qcal + AL (3) Trong đó, M AL - hệ số chuyển đổi, cung cấp file metadata ảnh L , LANDSAT 8; Qcal – giá trị số kênh ảnh [10] Giá trị xạ phổ tính bước dùng để tính nhiệt độ độ sáng (brightness temperature) theo công thức [10]: K (4) TB = K ln(1 + ) Lλ Trong đó, hệ số �1, �2 cung cấp file metadata ảnh LANDSAT (bảng 1) [10]  _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ B ả n g G i t r ị K , K đ ố i v i d ữ l i ệ u ả n h h n g ngoại nhiệt LANDSAT K ê n h 1 V ệ K t i n h ( W / m L A N D S A T L A N D S A T L A N D S A T cô ) ng th ức sa u [1 , 2, 911 ]: s r µ m ) 7 , , LST T λ − tr on g , λ gi trị bư ớc só ng tr un g tâ m kê nh hồ ng ng oạ i nh iệt ; ρ h 6 , Nhiệt độ độ sáng hiệu chỉnh sở độ phát xạ bề mặt để xác định nhiệt độ bề mặt (land surface temperature) theo J se c c – vậ hằn tố g c Stef án Bolt zma nn ( 8.10 ợ c m/s x h ) c ε – độ đ hằ phá ị t ng n xạ h số bề mặ Pl t t an (su rfa k ce t 6, em issi 62 vit l y) i ệ 10 Đ u ộ − n – 10 σ K σ an * J h p h t x b ề m ặ t c ó sá t ng h (2, ể 99 đ v i ễ n t h m d ự a t r ê n k ế t q u ả p h â n l o i c c l o i h ì n h l p p h ủ h o ặ c c h ỉ s ố t h ự c v ậ t N D V I ; t r o n g đ ó , p h n g p h p d ự a t r ê n c h ỉ số thự c vật N D VI có nhi ều ưu điể m xá c địn h độ ph át xạ chi tiết đế n từn g pix el Để xá c địn h độ ph át xạ bề mặ t, ε = ε v Pv + ε s (1 − Pv ) (6) Pv đó: �𝑣,�� độ tỉ lệ phát xạ đặc trưng thực cho đất thực vật vật nhất, tro Pv có giá trị đối ng với đất trống m khu vực ột pi xe l phủ kín Pv xác định thực theo công thức vật [15] [15]: Giá trị  NDVI − NDVI  v P =   max −  NDVI NDVI  đó, số thực vật NDVI (normalized difference vegetation index) tỉ số hiệu số giá trị phản xạ phổ bề mặt kênh cận hồng ngoại (NIR) kênh đỏ (RED) tổng chúng Chỉ số NDVI ảnh LANDSAT xác định sau [2, 4, 9- 11]: NDVI = ρ NIR −ρ RED +ρ ρ NIR RED Để xác định độ phát xạ sở số NDVI yêu cầu phải biết trước độ phát xạ đất thực vật Nhiều nghiên cứu nước xác định độ phát xạ đất thực vật lấy số liệu phát xạ có sẵn thơng qua kết đo đạc thực nghiệm mẫu đại diện [2, 4] Điều dẫn đến sai số khu vực khác có bề mặt với đặc trưng vật lí khác Để khắc phục vấn đề trên, nghiên cứu này, tác giả tiến nghiệm với 120 vùng mẫu, có 60 mẫu vùng thực vật 60 mẫu vùng đất trống khu vực Bắc Bình (Bình Thuận) liệu ảnh số NDVI chụp vào 05 – 01 – 2002, 13 – 01 – 2005, 12 – 02 – 2010 15 – 02 – 2014 (hình 3) a) b) c) d) Hình Chỉ số thực vật NDVI khu vực Bắc Bình (Bình Thuận) ảnh Landsat ngày 05 – 01 – 2002( a), 13 – 01 – 2005 (b), 12 – 02 – 2010 (c) 15 – 02 – 2014 (d) Trịnh Lê Hùng tgk TẠP CHÍ KHOA HỌC ĐHSP TPHCM _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ K ết thực nghiệ m nhận cho thấy, giá trị NDV I cho đất trống đất phủ kín thực vật ảnh LAN DSA T khu vực nghiê n cứu tươn g ứng 0,127 0,515 Sử dụng phươ ng pháp 11 TẠP CHÍ KHOA HỌC ĐHSP TPHCM Van De Griend [16], độ phát xạ đất thực vật nhận tương ứng 0,912 0,978 theo công thức: ε = 1,0094 + 0,047ln(N DVI) 12 Số 5(70) năm 2015 a) b) c) d) H ình Kết xác định nhiệt độ bề mặt khu vực huyện Bắc Bình (Bình Thuận) ảnh Landsa t ngày 05 – 01 – 2002 (a), 13 – 01 – 2005 (b), 12 – 02 – 2010 (c) 15 – 02 TẠP CHÍ KHOA HỌC ĐHSP TPHCM – 2014 (d) Trịnh Lê Hùng tgk 13 _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ K ết qu ả xá c địn h nhi ệt độ bề mặ t kh u vự c hu yệ n Bắ c Bì nh (Bì nh Th uậ n) đư ợc trì nh bà y trê n hìn h 4(a - d) Phân tích kết nhận cho thấy, vùng có nhiệt độ bề mặt cao phân bố chủ yếu khu vực khơng có thực vật che phủ Sự chênh lệch nhiệt độ khu vực đất rừng khu vực khơng có thực vật lớn, tương ứng với 24,20C, 26,50C, 23,70C, 26,60C ảnh ngày 05–01–2002, 13– 01–2005, 12–02–2010 15–02–2014 Để tính “cạnh khô” 𝑇�𝑚𝑎𝑥 không gian “nhiệt độ/NDVI”, nghiên cứu chia số thực vật NDVI thành 15 khoảng xác định giá trị nhiệt độ cực đại khoảng Cạnh ướt 𝑇�𝑚i� lấy giá trị nhiệt độ bề mặt thấp xác định bước Từ kết hồi quy tuyến tính giá trị nhiệt độ bề mặt cực đại khoảng giá trị NDVI, “cạnh khô” 𝑇�𝑚𝑎𝑥 không gian Ts/NDVI cho ảnh năm 2002, 2005, 2010 2014 xác định sau: 𝑇�(𝑚𝑎𝑥)(2002) = −21,64𝑁𝐷𝑉𝐼 + 313,7 (��) 𝑇�(𝑚𝑎𝑥)(2005) = −24,94𝑁𝐷𝑉𝐼 + 314,8 (��) 𝑇�(𝑚𝑎𝑥)(2010) = −16,39𝑁𝐷𝑉𝐼 + 317,2 (��) 𝑇 �( 𝑚 𝑎 𝑥) ( ) = − , � � � � � � 𝐼 + , ( � � � ) B ả n g P h â n c ấ p m ứ c đ ộ k h ô h n đ ố v ỉ s ố T V D I n h ậ n g i c h t r ị s ố T V D I t r o n g k h TVDo ả n g 0,2 t 0,4 B ả n đ p h â n v ù n g m ứ c đ ộ k h ô h n t n g 0,6 đ ố i 0,8 đ ế n C h k h u v ự c h u y ệ n B ắ c B ì n h ( B ì n h T h u ậ n ) t ỉ l ệ : 0 Số 5(70) năm 2015 TẠP CHÍ KHOA HỌC ĐHSP TPHCM _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ vực dự báo khơng có nguy khơ hạn khơ hạn nhẹ tập trung chủ yếu phía Bắc nơi có thảm rừng che phủ Mặc dù vậy, diện tích khu vực giai đoạn 2002 – 2014 có xu hướng giảm tương đối rõ rệt với suy giảm diện tích rừng, tương ứng năm 2002, 2005, 2010 2014 17,39%, 12,72%, 11,91% 10,74% tổng diện tích huyện Trong đó, diện tích khu vực có nguy khơ hạn nặng nặng tăng mạnh năm 2010 2014 so với năm 2002, 2005 tập trung chủ yếu khu vực đất sản xuất nông nghiệp, đất cát Trong năm 2002 2005, diện tích khu vực có nguy khơ hạn nặng (tương ứng với giá trị số TVDI lớn 0,8) chiếm tỉ lệ thấp (0,16% 0,06%) Tuy nhiên, năm 2010 2014, diện tích tăng lên cách nhanh chóng, ước tính chiếm 12,54% 12,57% Tổng diện tích khu vực có nguy khơ hạn nặng nặng năm 2002, 2005 chiếm 35,67% 38,96% diện tích toàn huyện so với 55,27% 68,92% năm 2010, 2014 (bảng 3) Như vậy, nhìn chung, tượng hạn hán Bắc Bình (Bình Thuận) có xu hướng diễn ngày mạnh, ảnh hưởng nghiêm trọng đến môi trường sống hoạt động sản xuất người dân Kết nhận nghiên cứu cho thấy vai trò quan trọng lớp phủ thực vật giảm nguy hạn hán Ở vùng phủ xanh rừng trồng (thuộc xã Phan Hiệp, Hịa Thắng, Hồng Thái), mức độ khơ hạn từ cao năm 2005 giảm xuống mức trung bình năm 2010 2014 (hình 5) Bảng Phân vùng mức độ khô hạn tương đối khu vực Bắc Bình (Bình Thuận) Diện tích (%) 16 Năm Khơng khơ hạn Khơ hạn nhẹ Khơ hạn trung bình Khơ hạn nặng Khô hạn nặng 2002 0,31 17,08 46,94 35,51 0,16 2005 0,33 12,39 48,32 38,90 0,06 2010 0,05 11,86 32,82 42,73 12,54 2014 0,09 10,75 30,24 46,35 12,57 _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ a) b) c) d) Hình Bản đồ phân cấp mức độ khơ hạn tỉ lệ 1:100 000 khu vực Bắc Bình (Bình Thuận) năm 2002 (a), 2005 (b), 2010 (c) 2014 (d) Kết luận Ảnh hưởng biến đổi khí hậu hoạt động người năm gần làm tăng nguy hạn hán khu vực huyện Bắc Bình (Bình Thuận) Phần lớn diện tích huyện Bắc Bình dự báo có mức độ khơ hạn từ trung bình đến khơ hạn nặng nặng, diện tích khu vực có nguy khơ hạn nặng nặng tăng nhanh năm 2010, 2014 so với năm trước Những khu vực có nguy khơ hạn cao tập trung chủ yếu vùng sản xuất nơng nghiệp vùng khơng có thực vật che phủ Tư liệu ảnh vệ tinh quang học Landsat với ưu điểm độ phân giải khơng gian trung bình, tích hợp kênh hồng ngoại nhiệt đặc biệt cung cấp hồn tồn miễn phí với chu kì cập nhật 16 ngày nguồn tư liệu phong phú quý giá nghiên cứu giám sát tượng hạn hán Kết nhận nghiên cứu sử dụng thành lập đồ nguy hạn hán tỉ lệ 1:100.000, góp phần ứng phó giảm thiểu ảnh hưởng hạn hán đến môi trường sống hoạt động sản xuất người dân TÀI LIỆU THAM KHẢO Trịnh Lê Hùng (2014), “Nghiên cứu phân bố nhiệt độ bề mặt liệu ảnh vệ tinh đa phổ LANDSAT”, Tạp chí Các khoa học Trái Đất, 36(1), 82 – 89 Trịnh Lê Hùng (2014), “Ứng dụng liệu viễn thám hồng ngoại nhiệt LANDSAT nghiên cứu độ ẩm đất sở số khô hạn nhiệt độ thực vật”, Tạp chí Các khoa học Trái Đất, 36(3), 262 – 270 Phan Thị Kim Vân (2006), “Nghiên cứu nguồn nước Bắc Bình sở phân tích yếu tố khí tượng”, Tạp chí Các khoa học Trái Đất, 28(4), 422 – 433 Trần Thị Vân, Hoàng Thái Lan, Lê Văn Trung (2009), “Phương pháp viễn thám nhiệt nghiên cứu phân bố nhiệt độ bề mặt thị”, Tạp chí Các khoa học Trái Đất, 31(2), tr.168 – 177 Phạm Quang Vinh, Phạm Thị Thanh Hương (2012), “Đánh giá hạn nông nghiệp tỉnh Bình Thuận theo kịch biến đổi khí hậu”, Tạp chí Các khoa học Trái Đất, 34(4), 513 – 523 Arnon Karnieli, Nurit Agam, Rachel T Pinker, Martha Anderson, Marc L Imhoff, Garik G Gutman, Natalya Panov, Alexander Goldberg (2009), “Use of NDVI and land surface temperature for drought assessment: merits and limitations”, Journal of Climate, 23, 618 – 633 Chavez P.S (1996), “Image-based atmospheric corrections–revisited and improved”, Photogrammetric Engineering and Remote Sensing, 62(9): pp.1025-1036 Chavez P.S (1988), “An improved dark-object subtraction technique for atmospheric scattering correction of multispectral data”, Remote Sensing of Environment, 24: pp.459-479 Lu Yuan, Tao Heping, Wu Hua (2007), “Dynamic drought monitoring in Guangxi using revised temperature vegetation dryness index”, Wuhan University journal of Natural sciences, 12(4), pp 663 – 668 10 National Aeronautics and Space Administration (NASA), LANDSAT Science data user’s Handbook, 270 pp Trịnh Lê Hùng tgk TẠP CHÍ KHOA HỌC ĐHSP TPHCM _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 11 _ Ru lin da C M., Bij ker W., Ste in A (20 10) , “I ma ge mi nin g for dro ug ht mo nit ori ng in Ea ste rn Afr ica usi ng Me teo sat SE RV IRI dat a”, Int 20 ... tích khu vực có nguy khơ hạn nặng nặng tăng mạnh năm 2010 2014 so với năm 2002, 2005 tập trung chủ yếu khu vực đất sản xuất nông nghiệp, đất cát Trong năm 2002 2005, diện tích khu vực có nguy. .. hạn từ trung bình đến khơ hạn nặng nặng, diện tích khu vực có nguy khơ hạn nặng nặng tăng nhanh năm 2010, 2014 so với năm trước Những khu vực có nguy khô hạn cao tập trung chủ yếu vùng sản xuất... Bình Thuận tác động tiêu cực từ hoạt động có nguy suy thối đất người làm cho tình trạng hạn hán hoang mạc hóa hàng tỉnh Nam Trung Bộ, Tây Nguy? ?n nói đầu khu vực dun hải chung Bình Thuận nói riêng