ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Phan Thành Bắc MÔ PHỎNG QUÁ TRÌNH LAN TRUYỀN VẬT CHẤT Ô NHIỄM DƯỚI TÁC ĐỘNG CỦA CÁC YẾU TỐ ĐỘNG LỰC TẠI VỊNH CAM RANH BẰNG MÔ HÌNH SỐ Chuyên ngành: Hải dương học Mã số: 60.44.97 LUẬN VĂN THẠC SĨ KHOA HỌC Người hướng dẫn khoa học: PGS TS Nguyễn Minh Huấn Hà Nội LỜI CẢM ƠN Tác giả xin dành lời bày tỏ lòng biết ơn sâu sắc tới thầy, cô giáo khoa Khí tượng, Thủy văn - Hải dương học (trường Đại học Khoa học tự nhiên Hà Nội) nhà khoa học viện Hải dương học tận tình giúp đỡ, truyền thụ, trao đổi kiến thức chuyên môn tác giả thời gian qua Luận văn hoàn thành nỗ lực làm việc thân có công lớn thầy giáo PGS.TS Nguyễn Minh Huấn, người không ngừng đôn đốc, động viên truyền thụ kiến thức Tác giả xin gửi lời biết ơn chân thành sâu sắc đến thầy Tác giả xin gửi lời cảm ơn tới ThS Nguyễn Chí Công tất cán nghiên cứu phòng Vật Lý Biển nói riêng, Viện Hải Dương học – nơi tác giả công tác nói chung, giúp đỡ nhiệt tình nguồn số liệu sử dụng Bên cạnh đó, tác giả xin gửi lời cảm ơn đến dự án “Nghiên cứu khả tự làm sạch, đề xuất giải pháp nhằm bảo vệ cải thiện chất lượng môi trường đầm Thủy Triều - vịnh Cam Ranh” PGS.TS Bùi Hồng Long ThS Nguyễn Hữu Huân đồng chủ nhiệm, cho phép sử dụng nguồn số liệu phục vụ cho luận văn Suốt trình học tập nghiên cứu luận văn, tác giả giúp đỡ từ dự án chống biến đổi khí hậu CLIMEEViet, hợp tác nghiên cứu Viện Hải Dương học với phủ Đan Mạch, mà đứng đầu PGS.TS Nguyễn Ngọc Lâm Tác giả xin trân trọng cảm ơn dự án tài trợ mặt kinh phí, thiết bị hỗ trợ nghiên cứu nguồn số liệu tham khảo vô quí giá Qua đây, tác giả xin chân thành cảm ơn giúp đỡ thân tình bạn bè, thân hữu trình học tập nghiên cứu trường Đại học Khoa học tự nhiên MỤC LỤC MỞ ĐẦU CHƯƠNG MÔ HÌNH SỐ TRỊ 1.1 TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU 1.1.1 Tổng quan tình hình nghiên cứu giới 1.1.2 Tổng quan tình hình nghiên cứu nước 1.2 MIKE 21 HD 1.2.1 Cơ sở toán học 1.2.2 Phương pháp số 11 1.3 MÔĐUN ECOLAB 14 1.3.1 Cơ sở lý thuyết 14 1.3.2 Ôxy hòa tan (DO) nhu cầu ôxy sinh hóa (BOD) 16 1.3.3 Các hợp phần Nitơ .19 1.3.4 Hợp phần Photpho .22 CHƯƠNG TỔNG QUAN VÙNG NGHIÊN CỨU 23 2.1 TỔNG QUAN VỀ ĐIỀU KIỆN TỰ NHIÊN 23 2.1.1 Vị trí địa lí 23 2.1.2 Đặc điểm gió .24 2.1.3 Đặc điểm thủy, hải văn .24 2.1.4 Đặc điểm nhiệt - muối 26 2.1.5 Đặc điểm dòng chảy 27 2.1.6 Đặc điểm thủy triều dao động mực nước .28 2.2 ĐẶC ĐIỂM KINH TẾ - XÃ HỘI 28 2.3 HIỆN TRẠNG MÔI TRƯỜNG VỊNH CAM RANH 29 2.3.1 Các nguồn thải 29 2.3.2 Chất lượng nước vịnh Cam Ranh .30 CHƯƠNG ÁP DỤNG MÔ HÌNH VÀ KẾT QUẢ 31 3.1 THIẾT LẬP CÁC THÔNG TIN ĐẦU VÀO CHO MÔ HÌNH 31 3.1.1 Thu thập số liệu 31 3.1.2 Địa hình đáy .32 3.1.3 Thiết lập lưới tính .32 3.1.4 Điều kiện biên điều kiện ban đầu 34 3.2 HIỆU CHỈNH MÔ HÌNH 36 3.3 MỘT SỐ KẾT QUẢ TÍNH TOÁN 39 3.3.1 Kết tính toán cho mùa khô 39 3.3.2 Kết tính toán cho mùa mưa 70 3.4 ĐÁNH GIÁ TÁC ĐỘNG CỦA CÁC YẾU TỐ Ô NHIỄM 104 KẾT LUẬN 125 KIẾN NGHỊ 126 TÀI LIỆU THAM KHẢO 127 MỞ ĐẦU Trong năm gần đây, khu vực đầm thuỷ triều đứng trước nguy ô nhiễm nguồn nước Đầm Thủy Triều nằm vịnh Cam Ranh, thuộc địa bàn huyện Cam Lâm thành phố Cam Ranh, tỉnh Khánh Hòa Nơi phong phú đa dạng số lượng trữ lượng thủy sản Trong tương lai, đầm Thủy Triều mắt xích quan trọng việc phát triển du lịch tỉnh Khánh Hòa vịnh Cam Ranh tỉnh quy hoạch thành trung tâm du lịch biển tầm cỡ quốc gia quốc tế đến năm 2025 Theo nhận định người dân nơi đây, vòng gần chục năm nay, tôm, cá loại nghêu, ốc đầm thường chết hàng loạt, chí “sống dai” loài sá sùng biển (gọi trùn biển) phải chết trắng đầy đầm, môi trường đầm ngày trở nên ngột ngạt, đục ngàu, nước đầm có mùi hôi thối nồng nặc theo thời gian làm cho hệ sinh thái đầm bị biến dạng, nguồn lợi thủy sản không Do vậy, việc đánh khai thác nguồn lợi đầm không hiệu quả, đời sống nhân dân lại khốn khó Một nguyên nhân gây ô nhiễm đầm nhà máy đường Cam Ranh Trong trình vận hành nhà máy, khối nước thải từ nhà máy sau xử lí sẽ đổ đầm qua cống xả thải Các kết từ phân tích mẫu nước vị trí cống xả thải khu vực xung quanh nhà máy ghi nhận vượt ngưỡng thông số môi trường xung quanh khu vực Khi khối nước thải xả đầm, trình thuỷ động lực (dòng chảy, gió, trình xáo trộn,…) làm khuếch tán chất đồng thời mang khối nước thải lên phía bắc hay xuống phía nam theo dòng chảy triều lên triều xuống Vì vậy, trình động lực khu vực đóng vai trò quan trọng việc phân bố, truyền tải, pha loãng, làm vùng đầm thuỷ triều Một cách tiếp cận để nghiên cứu ảnh hưởng khối nước thải từ nhà máy đường sử dụng mô hình tính toán để tính toán mô trình vật lý (dòng chảy) mô hình sinh hoá diễn khu vực đầm có tác động khối nước thải Các kết tính toán từ mô hình kết hợp với số liệu khảo sát mô cách liên tục trình động lực trình truyền tải vật chất mô kịch khác điều kiện động lực khác điều kiện xả thải khác Việc mô kịch ô nhiễm khác giúp nhà quản lý phản ứng linh hoạt hơn, hiệu tốn Từ đưa kế hoạch, chiến lược để quy hoạch, khai thác cách hiệu tài nguyên khu vực đầm việc kiểm soát điều tiết nguồn thải hợp lý Nhận thức mức độ cấp thiết vấn đề môi trường vịnh Cam Ranh, học viên lựa chọn hướng nghiên cứu với đề tài: “Mô trình lan truyền vật chất ô nhiễm tác động yếu tố động lực vịnh Cam Ranh mô hình số” để mô số vật chất có khả ảnh hưởng đến chất lượng môi trường Có nhiều kỹ thuật đánh giá mức độ ô nhiễm nước dựa vào giá trị thông số chọn lọc Các kỹ thuật sử dụng số để thực mức độ ô nhiễm Trong nêu số số công nhận như: Chỉ số ô nhiễm dinh dưỡng (NPI) dựa vào thông số NH 4+, NO3-, NO2-, tổng P, pH, chlorophyll, độ dẫn điện độ đục Chỉ số ô nhiễm hữu (OPI) dựa vào thông số BOD, COD, nhiệt độ DO Với nguồn số liệu có từ số đề tài thực Viện Hải dương học đề tài cấp Cơ sở phòng Vật lý biển, phòng Thủy địa hóa, đề tài cấp Viện Khoa học Công nghệ, Các Dự án hợp tác quốc tế, tác giả sử dụng gói phần mềm MIKE 21 HD, ECO Lab để mô trình lan truyền số vật chất gây ô nhiễm từ nguồn thải khu công nghiệp, nuôi trồng thủy sản khu dân cư mùa: mùa mưa mùa khô Trong khuôn khổ luận văn, mục tiêu học viên tính toán, mô phỏng, đưa tranh trình động lực trình truyền tải vật chất gồm BOD, DO, NO3-, PO4+, NH3+ Một kịch mô lan truyền vật chất ô nhiễm với giả thiết có gia tăng cực đại nồng độ chất gây ô nhiễm từ số liệu thực đo cống xả thải công suất tính thời điểm khảo sát từ nguồn thải khu công nghiệp, nuôi trồng thủy sản khu dân cư để đánh giá mức độ lan truyền ảnh hưởng vật chất tới chất lượng nước bãi tắm khu vực Cam Ranh Các kết nghiên cứu luận văn góp phần bổ sung thêm thông tin khoa học nghiên cứu, đánh giá vai trò tác động từ nguồn thải khu công nghiệp, nuôi trồng thủy sản khu dân cư tác động ngược lại khu vực nuôi trồng thủy sản, du lịch sinh thái bãi tắm CHƯƠNG MÔ HÌNH SỐ TRỊ 1.1 TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU 1.1.1 Tổng quan tình hình nghiên cứu giới Sử dụng mô hình số để tính toán, mô phỏng, đánh giá chất lượng môi trường nước khu vực gần bờ, khu bãi tắm, khu nuôi trồng thủy sản thực phổ biến giới Tùy thuộc vào đối tượng mục đích nghiên cứu, việc áp dụng loại mô hình tính toán khác Có thể liệt kê số mô hình thường áp dụng để đánh giá chất lượng nước giới Mô hình WASP7 (Water Quality Analysis Simulation Program 7) mô hình xây dựng dựa mô hình trước (WASP – xây dựng Di Toro, 1983; Connolly vaf Winfield, 1984; Ambrose, R.B, 1988) Mô hình sử dụng để mô tả dự báo chất lượng nước giúp nhà quản lý đưa định, giải pháp đối phó với tượng ô nhiễm tự nhiên người Mô hình cho phép người sử dụng áp dụng không gian 1D mô tựa 2D 3D cách chia hộp với đa dạng thành phần chất ô nhiễm Mô hình WASP liên kết với mô hình thủy động lực vận chuyển trầm tích để thu trường dòng chảy, nhiệt độ, độ muối thông lượng trầm tích Mô hình WASP sử dụng để mô trình yếm khí vịnh Tampa; Cung ứng Photpho cho hồ Okeechobee; Quá trình yếm khí cửa sông Neuse River; Ô nhiễm vật chất hữu dễ phân hủy cửa sông Delaware, ô nhiễm kim loại nặng sông Deep, bắc Carolina Mô hình AQUATOX mô hình mô hệ sinh thái thủy sinh Mô hình dự báo trình suy tàn nhiều loại chất gây nhiễm môi trường dinh dưỡng, hóa học hữu cơ, ảnh hưởng chúng lên hệ sinh thái, bao gồm loài cá, động vật không xương sống loài thực vật thủy sinh AQUATOX công cụ hữu hiệu cho nhà môi trường học, sinh học, nhà mô hình hóa chất lượng nước cần quan tâm tới việc đánh giá rủi ro suy giảm hệ sinh thái thủy sinh Mô hình QUAL2K (hay Q2K) (River and Stream Water Quality Model) nâng cấp từ mô hình trước QUAL2E (hay Q2E (Brown Barnwell 1987)) Đây mô hình mô chất lượng nước suối sông chiều có tham gia trình xáo trộn rối bên Một đặc điểm linh hoạt mô hình chạy môi trường Visual basic môi trường Excel Mô hình có đặc điểm sau: tính toán phân đoạn sông nhánh sông Mô hình tính toán chu trình Nitơ Thông qua chu trình chuyển hóa nitơ để biểu diễn hợp chất cacbon (loại ôxy hóa nhanh chậm), loại cacbon hữu không sống (các phân tử cacbon, nitơ, phôtpho hợp chất hóa học) Các trình thiếu hụt ôxy gần tới giá trị không trình ôxy hóa, trình khử nitơ bước tương tác Tính toán thông lượng trao đổi ôxy hòa tan dinh dưỡng trầm tích nước DELFT 3D Viện nghiên cứu thuỷ lực Hà Lan cho phép kết hợp mô hình thuỷ lực chiều với mô hình chất lượng nước Ưu điểm mô hình việc kết hợp module tính toán phức tạp để đưa kết tính mô cho nhiều chất nhiều trình tham gia SMS Trung tâm nghiên cứu phát triển kỹ thuật quân đội Mỹ xây dựng cho phép kết hợp mô hình thuỷ lực 1, chiều với mô hình chất lượng nước, module RMA4 mô hình số trị vận chuyển yếu tố chất lượng nước phân bố đồng theo độ sâu Nó tính toán tập trung thành phần bảo toàn không bảo toàn tính toán theo lưới chiều chiều ECOHAM (phiên 2) mô hình số 3D kết hợp module thủy lực với module sinh thái phát triển nhóm nghiên cứu Trường đại học Hamburg (Đức) Mô hình chủ yếu tính toán dựa chu trình hợp phần Nitơ Photpho có tính đến thực vật động vật phù du nước biển ECOSMO (ECOSystem MOdel) mô hình cặp ba chiều thủy động lực – băng biển – sinh địa hóa Mô hình phát triển dựa mô hình thủy động lực HAMSOM (HAMburg shelf Ocean Model) liên kết mô đun động lực nhiệt động lực biển - băng (Schrum Backhaus, 1999) môđun sinh học (Schrum, 2006) Môđun sinh học NPZD dựa trình chuyển đổi mức thứ hai chuỗi thức ăn điểu khiển thông lượng Nitơ, Photpho Silic Điều quan trọng tính toán mô hình thống giới hạn chu trình dinh dưỡng vĩ mô động vật phù du mô hình chuẩn đoán biến đổi cho tương tác phi tuyến hệ sinh thái mức thứ thứ hai chuỗi thức ăn Thêm vào đó, mô hình tính toán biến đổi mảnh vụn ôxy để đánh giá lượng lại trình ôxy hóa Các tính toán sinh khối sơ cấp thứ cấp Mô hình ECOSMO áp dụng cách thành công việc mô tả khu vực có động lực dinh dưỡng yếu khu vực Biển Bắc BASINS EPA nhằm trợ giúp đánh giá kiểm tra hệ thống liệu thông tin môi trường, giúp hệ thống phân tích môi trường phân tích phương án quản lý Một điểm bật BASINS đưa vào cách tiếp cận dựa tảng lưu vực sông, có kết hợp quản lý liệu không gian thông qua hệ thông tin địa lý GIS BASINS dùng cho mục đích sau: Mô điều kiện lưu vực đánh giá trạng chất lượng nước; Mô tác động việc thay đổi sử dụng đất có tính đến cân nước, mô kịch nguồn ô nhiễm điểm diện, xây dựng phát triển cách quản lý lưu vực Các nhóm tham số mô hình bao gồm: Các hợp chất dinh dưỡng Nitơ Photpho, DO, BOD, thuốc trừ sâu, thuốc bảo vệ thực vật, bùn Bộ phần mềm MIKE Viện Thuỷ lực Đan Mạch (DHI) phát triển thương mại hoá Một đặc điểm mạnh MIKE dễ sử dụng với giao diện Windows, kết hợp chặt chẽ với GIS (hệ thống thông tin địa lý) MIKE tích hợp module thuỷ lực (HD) chất lượng nước (ECO Lab), bao gồm: thuỷ lực, truyền tải - khuếch tán chất lượng nước MIKE mô hình với nhiều tính mạnh, khả ứng dụng rộng rãi cho nhiều dạng thuỷ vực khác 1.1.2 Tổng quan tình hình nghiên cứu nước Ở nước ta, năm gần đây, hướng nghiên cứu, xây dựng sử dụng mô hình nghiên cứu thủy động lực – môi trường quan tâm Trong nghiên cứu, điều tra, tính toán ô nhiễm môi trường vũng vịnh khu vực ven biển - khu vực tập trung chủ yếu hoạt động kinh tế CV CR CB N DC CP CT D CD TT SB C DL CL DT CL 0.04 032 0.05 318 0.03 983 0.08 048 0.03 769 0.07 816 0.01 154 0.00 478 Tên điể m Cực tiểu CV CR CB N DC CP CT D CD TT SB C DL CL DT CL 0.09 933 0.08 144 0.08 765 0.06 755 0.09 946 0.10 739 0.08 722 0.07 386 0.13 0.09 821 148 0.1 0.5 0.12 0.09 386 141 0.1 0.5 0.12 0.09 752 055 0.1 0.5 0.10 0.09 889 964 0.1 0.5 0.10 0.07 087 919 0.1 0.5 0.11 0.09 418 999 0.1 0.5 0.18 0.12 032 382 0.1 0.5 0.10 0.07 86 087 0.1 0.5 NO3-N (mg/l) QCV QCV Trun N N Cực g nước nước đại bình NTT bãi S tắm 0.11 0.10 079 607 0.09 0.08 69 942 0.10 0.09 303 442 0.09 0.08 134 138 0.11 0.11 459 043 0.11 0.11 489 135 0.12 0.11 998 409 0.09 0.08 584 523 0.00 644 0.05 286 0.04 636 0.02 939 0.05 475 0.03 182 0.04 681 0.03 943 0.05 69 0.06 919 0.05 353 0.06 939 0.08 536 0.04 802 0.08 387 0.04 66 0.03 612 0.05 629 0.04 937 0.05 036 0.06 674 0.03 923 0.05 929 0.04 179 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 0.1 0.2 Hình 3.66: Biểu đồ biểu diễn giá trị cực đại, cực tiểu trung bình yếu tố tính toán vị trí tuyến điểm, tháng 05/2011, thời kỳ mùa khô, kỳ triều cường 113 114 115 Hình 3.67: Biểu đồ biểu diễn giá trị cực đại, cực tiểu trung bình yếu tố tính toán vị trí tuyến điểm, tháng 05/2011, thời kỳ mùa khô, kỳ triều kiệt 116 117 118 Hình 3.68: Biểu đồ biểu diễn giá trị cực đại, cực tiểu trung bình yếu tố tính toán vị trí tuyến điểm, tháng 10/2011, thời kỳ mùa mưa, kỳ triều cường 119 120 121 Hình 3.69: Biểu đồ biểu diễn giá trị cực đại, cực tiểu trung bình yếu tố tính toán vị trí tuyến điểm, tháng 10/2011, thời kỳ mùa mưa, kỳ triều kiệt 122 123 124 KẾT LUẬN Qua kết thu từ mô hình tính toán, sở phân tích, thống kê, tác giả đưa số kết luận sau: - Việc sử dụng gói công cụ MIKE đáp ứng mục tiêu nội dung đặt luận văn Các kết tính toán từ mô hình mô trình dòng chảy, trình lan truyền chất gây ô nhiễm khu vực vịnh Cam Ranh từ nguồn thải - Dòng chảy khu vực vịnh Cam Ranh chủ yếu dòng triều chịu tác động chế độ gió mùa khu vực Trường phân bố dòng chảy từ cửa vịnh truyền vào chi làm hai hướng từ đông nam lên phía bắc truyền vào đầm thủy triều từ phía đông nam sang tây nam truyền vào vịnh Trong thời gian triều lên triều xuống có xuất xoáy thuận khu vực vịnh - Trong thời kỳ mùa khô, thành phần vật chất gây ô nhiễm khu vực vịnh Cam Ranh không ảnh hưởng tới chất lượng bãi tắm dọc theo bãi biển Cam Ranh, bắt đầu xuất dấu hiệu đạt ngưỡng giới hạn tiêu ô nhiễm từ khu vực cảng Ba Ngòi đến khu vực nhà máy chế biến thực phẩm bán kính 2km Trong thời kỳ mùa mưa, tiêu chất tiệm cận giới hạn cho phép có phần đáng quan ngại mùa khô lan truyền chất gây ô nhiễm vươn tới khu bãi tắm khu vực dân sinh xung quanh khu vực nói phạm vi lớn - Theo tiêu chuẩn chất lượng nước nuôi trồng thủy sản khu vực vịnh Cam Ranh có nhiều khu vực ven bờ bị ô nhiễm NH4, điển khu vực cảng Ba Ngòi, sân bay Cam Ranh cửa đầm Thủy Triều - Do yếu tố chủ quan nhà qui hoạch nhằm phục vụ phát triển thành phố Cam Ranh yếu tố khách quan đặc điểm khí tượng 125 thủy văn, địa hình, cấu tạo địa chất chủ yếu cát nên lượng nước sông truyền vào vịnh thấp, ảnh hưởng từ cửa sông tới chất lượng nước bãi tắm mùa mưa khô không đáng kể Nhưng khả truyền nước từ địa hình xung quanh vào vịnh qua đường nước ngầm lớn hạn chế điều kiện số liệu mô hình nên tác giả chưa tính toán - Tuy số hạn chế, song kết nghiên cứu luận văn góp phần làm sáng tỏ thêm tranh thủy động lực – môi trường vịnh Cam Ranh Góp phần giải vấn đề ô nhiễm môi trường, quy hoạch, sử dụng khai thác hợp lý tài nguyên môi trường KIẾN NGHỊ Mặc dù trạng môi trường vịnh Cam Ranh chưa vượt ngưỡng tiêu cho phép song đáng quan ngại Nhiều yếu tố môi trường bắt đầu tiệm cận ngưỡng giới hạn nội dung nghiên cứu luận văn trình bày Do qui hoạch thành phố tương lai cần thiết phải bổ sung thêm nhà máy xử lý chất thải nước thải để xử lý trước đổ vào vịnh tránh trường hợp phần lớn xả trực tiếp vào vịnh Bên cạnh đó, địa hình chủ yếu cát nên khả thấm nước mùa mưa cao nên cần phải có nghiên cứu chất lượng nguồn nước ngầm đổ vào vịnh nhằm phản ánh xác chất lượng môi trường vịnh Cam Ranh Ngoài ra, giới hạn kinh phí khả chuyên môn nên chưa có số liệu xác toàn diện mặt không gian thời gian nhằm tiếp tục nghiên cứu sâu để tiến tới xây dựng kịch lan truyền vật chất ô nhiễm số vị trí nhạy cảm đưa dự báo, cảnh báo, tư vấn cho ngành nghề liên quan đến môi trường biển 126 TÀI LIỆU THAM KHẢO Tiếng Việt Đoàn Văn Bộ (2001), “Hóa học biển”, NXB Đại học Quốc gia Hà Nội Bùi Hồng Long - Trần Văn Chung (2006), “Tính toán thử nghiệm dòng chảy ba chiều (3-D) cho vùng vịnh Vân Phong”, Tạp chí Khoa học Công nghệ biển, 6(1), tr 12-27 PGS.TS Bùi Hồng Long, ThS Nguyễn Hữu Huân (2011), “Nghiên cứu khả tự làm sạch, đề xuất giải pháp nhằm bảo vệ cải thiện chất lượng môi trường đầm Thủy Triều - vịnh Cam Ranh” GS.TS Mai Trọng Nhuận (2008), “Điều tra đánh giá tài nguyên môi trường vũng vịnh trọng điểm ven bờ phục vụ phát triển kinh tế - xã hội bảo vệ môi trường” Phạm Văn Thơm (2005), “Hiện trạng môi trường vịnh Cam Ranh” Phạm Văn Thơm (2008), “Điều tra đánh giá trạng môi trường vịnh Cam Ranh làm sở phát triển kinh tế” UBND tỉnh Khánh Hòa (2003), “Địa chí Khánh Hòa”, NXB Chính trị quốc gia UBND tỉnh Khánh Hòa (2010), “Quy hoạch phát triển tỉnh Khánh Hòa đến năm 2015 ” Sở khoa học công nghệ tỉnh Khánh Hòa (2004), “Đặc điểm khí hậu – thủy văn tỉnh Khánh Hòa” 10 http://www.khanhhoa.gov.vn Tiếng Anh 11 Mike Flow model (DHI 2007), Hydronamic module: Scientific Documentation 12 Mike Flow model (DHI 2007), ECO Lab module: Scientific Documentation 127