Giáo trình Ứng dụng mô hình Mike 11 mô phỏng quá trình lan truyền chất ô nhiễm do nuôi trồng thủy sản trên một số sông lớn tỉnh Quảng Trị- Nguyễn Vũ Anh Tuấn

Nội dung text: Giáo trình Ứng dụng mô hình Mike 11 mô phỏng quá trình lan truyền chất ô nhiễm do nuôi trồng thủy sản trên một số sông lớn tỉnh Quảng Trị- Nguyễn Vũ Anh Tuấn

1. Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất và Môi trường, Tập 32, Số 3S (2016) 250-255 Ứng dụng mô hình Mike 11 mô phỏng quá trình lan truyền chất ô nhiễm do nuôi trồng thủy sản trên một số sông lớn tỉnh Quảng Trị Nguyễn Vũ Anh Tuấn1,*, Nguyễn Quang Hưng2, Nguyễn Thanh Sơn2, Nguyễn Thị Liên2 1Tổng cục Hậu cần - Kỹ thuật, Bộ Công an, 80 Trần Quốc Hoàn, Cầu Giấy, Hà Nội, Việt Nam 2Khoa Khí tượng Thủy văn và Hải dương học, Trường Đại học Khoa học Tự nhiên, ĐHQGHN, 334 Nguyễn Trãi, Thanh Xuân, Hà Nội, Việt Nam Nhận ngày 08 tháng 8 năm 2016 Chỉnh sửa ngày 26 tháng 8 năm 2016; Chấp nhận đăng ngày 16 tháng 12 năm 2016 Tóm tắt: Bài báo trình bày các kết quả mô phỏng chất lượng nước trên các sông chính của tỉnh Quảng Trị có xét đến tác động của xả thải từ các điểm nuôi trồng thủy sản sử dụng mô hình Mike 11. Số liệu mưa và mực nước năm 2005 và 2009 được sử dụng để xây dựng mô hình thủy văn và thủy lực. Các số liệu quan trắc chất lượng nước hơn 30 điểm trên hệ thống sông của các năm 2013 và 2014, số liệu thực đo khảo sát các đơn vị nuôi trồng thủy sản được sử dụng để thiết lập mô hình chất lượng nước, mô phỏng quá trình lan truyền chất ô nhiễm trên sông. Các kết quả mô phỏng đã được so sánh với các số liệu thực đo để khẳng định độ tin cậy của mô hình và phân tích, làm sáng tỏ hiện trạng chất lượng nước cũng như làm cơ sở để xây dựng công cụ quy hoạch nuôi trồng thủy sản trên khía cạnh bảo vệ môi trường. Từ khóa: Chất lượng nước, nuôi trồng thủy sản, Mike 11 WQ, Quảng Trị. 1. Khái quát khu vực nghiên cứu1 sinh kinh tế của người dân phụ thuộc rất lớn vào chế độ thủy văn cũng như chất lượng nước Quảng Trị nằm trong phạm vi từ 16018 đến của các hệ thống sông này đặc biệt là các thành 17010 vĩ độ Bắc; 106032 đến 107024 kinh độ phần kinh tế nông nghiệp và nuôi trồng thủy Đông. Phía Bắc giáp với tỉnh Quảng Bình, phía sản (NTTS). Nam giáp với tỉnh Thừa Thiên Huế, phía Tây Cũng như một số tỉnh miền Trung khác, giáp tỉnh Savanakhet và Salavan - nước nuôi trồng thủy sản ở tỉnh Quảng Trị có xu CHDCND Lào, phía Đông là biển Đông với hướng phát triển mạnh trong khoảng thời gian chiều dài bờ biển là 75km [1]. 10 năm trở lại đây. Tổng diện tích nuôi trồng Quảng Trị là một tỉnh nghèo miền Trung, thủy sản của Quảng Trị năm 2003 là 1730 ha, dân số chủ yếu sống tập trung ở dải ven biển và năm 2006 là 2380,4 ha và đến năm 2013 là hơn dọc theo các hệ thống sông chính là sông Bến 3,1 nghìn ha [2, 4]. Theo số liệu thống kê thu Hải và Thạch Hãn, vì thế các hoạt động dân được của Sở Nông nghiệp và phát triển nông thôn tỉnh Quảng Trị, trong năm 2013 tính riêng ___ diện tích nuôi trồng thủy sản nước ngọt là 2185 * Tác giả liên hệ. ĐT.: 84-983250586 ha và nước mặn, nước lợ là 1022 ha [4]. Email: anhtuan86bca@gmail.com 250
2. N.V.A. Tuấn và nnk. / Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất và Môi trường, Tập 32, Số 3S (2016) 250-255 251 2. Vấn đề ô nhiễm do nuôi trồng thủy sản diễn ra quy mô lớn ở vùng mặn hóa ven biển làm gia tăng xâm nhập mặn ở các vùng ven Bên cạnh các lợi ích kinh tế mà NTTS biển. Nhiều ao hồ nuôi tại các khu đất ở sinh mang lại thì sự phát triển của lĩnh vực này cũng hoạt, hoặc các khu đất ven biển, ven sông, được có những tác động tiêu cực đến môi trường lót đáy, nên nước nuôi không trực tiếp ngấm trong và ngoài ao nuôi. xuống đất. Tuy nhiên, khi xả thải không xử lý Đối với các ao hồ nuôi thủy sản nước ngọt, sẽ gây ảnh hưởng đến nguồn nước ngầm. Đồng hình thức nuôi là bán thâm canh hoặc nuôi thả, thời, quá trình nuôi trồng đòi hỏi bổ sung một nước thải ra chủ yếu trong các đợt thay nước lượng lớn nước sạch, chủ yếu từ các giếng nuôi và khi thu hoạch, được xả trực tiếp vào khoan dẫn đến khả năng làm sụt mực nước môi trường không qua xử lý. Các hộ nói chung ngầm, sụt lở đất [1, 3]. không sử dụng các loại thuốc kháng sinh, tuy Theo Cơ quan bảo vệ môi trường nhiên ô nhiễm vẫn xuất hiện do thức ăn thừa (Environemtal Protection Agency - EPA) của gây nên [1]. Hoa Kỳ - hoạt động nuôi trồng thủy sản phát sinh nhiều yếu tố có thể gây nguy hại đến môi Đối với các khu nuôi nước mặn, lợ thì nước trường thủy sinh nếu xả thải vào môi trường thải được xả thẳng ra môi trường với lưu lượng nước tiếp nhận một lượng đáng kể [6]. và tần suất lớn. Đối với tôm sú và tôm chân trắng, thời gian thay nước có thể lên tới 20% Chất rắn lơ lửng tăng làm tăng độ đục, gây lượng nước nuôi trồng trong 1 ngày, mang theo cản trở sự xâm nhập của ánh mặt trời, qua đó thức ăn thừa và các chất thuốc kháng sinh, ảnh làm giảm hoạt động quang hợp và tiêu thụ ô-xy hưởng lớn đến môi trường xung quanh [1]. của thực vật thủy sinh và phù du, lượng ô-xy hòa tan trong nước sẽ giảm và khiến động vật Thành phần lớp bùn trong các đầm, ao (tôm, cá, ) chết theo. Việc gia tăng chất lơ NTTS chủ yếu là các chất hữu cơ như prôtêin, lửng còn khiến nhiệt độ nước bề mặt tăng cao lipids, axit béo, photpholipids, Sterol, vitamin do hấp thụ nhiều ánh mặt trời, cũng sẽ làm D3, các hoocmon, carbohydrate, chất khoáng và giảm lượng ô-xy hòa tan trong nước. Chất rắn vitamin, vỏ tôm lột xác, Lớp bùn này luôn ở lơ lửng làm xước xát mang và da của động vật trong tình trạng ngập nước, yếm khí, các vi sinh thủy sinh nhạy cảm, tăng nguy cơ nhiễm trùng vật yếm khí phát triển mạnh, phân huỷ các hợp và bệnh tật. Khi lắng xuống đáy, các hạt rắn còn chất trên tạo thành các sản phẩm là có thể làm hỏng trứng cá, gây cản trở cho sự hydrosulphua (H2S), a-mô-ni-ắc (NH3), khí sinh sản của các loài thủy sản. metan (CH4), rất có hại cho thuỷ sinh vật. Khí Ni-tơ cùng với Phốt-pho là hai dưỡng chất a-mô-ni-ắc (NH3) cũng được sinh ra từ quá trình chính trong các khu NTTS được xả thải ra môi phân huỷ yếm khí thức ăn tồn dư gây độc trực tiếp trường. Việc các nguồn nước tiếp nhận dư thừa cho tôm, làm ảnh hưởng đến độ pH của nước và quá nhiều dưỡng chất sẽ dẫn đến hiện tượng kìm hãm thực vật phù du phát triển [5]. phú dưỡng, gây hậu quả như là thực vật phát Việc sử dụng các chất hóa học trong việc triển quá mức, đục nước, DO thấp, cá chết, suy bảo vệ sức khoẻ động vật thuỷ sản đang diễn ra giảm hệ sinh vật, kích thích hoạt động của vi phổ biến. Ngành NTTS công nghiệp sử dụng sinh vật, trong đó có một số chủng loài có hại một lượng lớn các loại thuốc kháng sinh và các đến sức khỏe con người. loại thuốc khác và sự lạm dụng hoặc dùng Chất hữu cơ được thải từ các khu NTTS không đúng cách thì tồn dư của chất hóa học chủ yếu do phân và thức ăn thừa. Nồng độ chất này có thể kéo theo những hậu quả nhất định hữu cơ gia tăng sẽ làm tăng nhanh hiện tượng đến hệ sinh thái thủy sinh và nguồn nước tiếp phú dưỡng và suy giảm ô-xy (vi sinh vật cần nhận xung quanh [3, 5]. tiêu thụ ô-xy để phân hủy chất hữu cơ). Quá trình chuyển dịch sử dụng đất từ các Sự gia tăng về diện tích NTTS nước mặn, lợ diện tích trồng lúa sang nuôi trồng thủy sản dẫn tới sự gia tăng tải lượng các chất hữu cơ và
3. 252 N.V.A. Tuấn và nnk. / Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất và Môi trường, Tập 32, Số 3S (2016) 250-255 cặn lơ lửng. Tuy nhiên do khả năng tự làm sạch Lăng rồi nhập với hệ thống sông Ô Lâu trước hiện thời của hai con sông là lớn nên chưa gây khi đổ ra biển. Sơ đồ thủy lực mô phỏng phần ra hậu quả nghiêm trọng. Tuy nhiên, nồng độ hạ lưu được tính toán của các sông này được chất ô nhiễm hữu cơ và cặn lơ lửng đặc biệt lớn minh họa trên Hình 1. tại các vị trí xả thải của các ao nuôi và trong thời gian các ao nuôi xả thay nước. 3.3. Số liệu 3. Ứng dụng mô hình MIKE 11 mô phỏng Toàn bộ mạng tính toán 1 chiều được thiết quá trình lan truyền chất ô nhiễm và tự lập với 140 mặt cắt, 398 nút tính toán. Để mô làm sạch phỏng thủy lực cho hệ thống sông Bến Hải – Thạch Hãn, số liệu lưu lượng tại 6 biên trên và Để đánh giá tình hình ô nhiễm chất lượng số liệu mực nước tại 3 biên dưới được sử dụng nước trên các hệ thống sông tỉnh Quảng Trị, bao gồm: phương pháp được sử dụng là phương pháp mô Số liệu biên trên: Cầu Sa Lung, Gia Vòng, hình tính thủy động lực kết hợp với lan truyền Cam Tuyền, Đăkrông, Hải Sơn, Phò Trạch. Căn chất và quá trình biến đổi sinh học của các hợp cứ trên điều kiện số liệu quan trắc, trong số các chất trong sông. Nhóm nghiên cứu đã lựa chọn biên trên duy nhất có biên Gia Vòng lấy giá trị bộ mô hình MIKE 11 với các mô đun thủy động thực đo lưu lượng, các biên còn lại sử dụng tài lực (HD), mô đun lan truyền chất (AD) và mô liệu dòng chảy tính toán từ mưa bằng mô hình đun chất lượng nước (WQ). thủy văn (NAM) với bộ thông số đã được hiệu 3.1. Lý thuyết mô hình MIKE 11 [7, 8 ] chỉnh và kiểm định. Số liệu biên dưới: Mực nước triều tại Cửa Mô đun thủy lực (HD): được xây dựng trên Tùng, Cửa Việt, Cửa Lác. Mực nước triều tại cơ sở hệ phương trình Saint Venant bao gồm hai Cửa Việt, Cửa Lác không có trạm đo, hiệu phương trình, để xác định lưu lượng và mực nước chỉnh theo mực nước tại Cửa Việt. tại các mặt cắt trên hệ thống sông, đây chính là cơ Các số liệu chất lượng nước được đưa vào sở để xây dựng hầu hết các mô đun khác. tính toán mô phỏng trong nghiên cứu bao gồm Mô đun lan truyền chất (AD): được dùng để DO, BOD, Tổng N. Nguồn số liệu đo đạc chất mô phỏng vận chuyển một chiều của chất lượng nước trên sông hệ thống sông Bến Hải - huyền phù hoặc phân hủy trong các lòng dẫn hở Thạch Hãn 10/2014 của Sở Tài nguyên và Môi dựa trên phương trình khuếch tán. trường tỉnh Quảng Trị được sử dụng để mô Mô đun chất lượng nước (WQ): giải quyết phỏng nền chất lượng nước. khía cạnh chất lượng nước trong sông. Mô đun này luôn đi kèm với mô đun lan truyền chất, có Nguồn gây ô nhiễm xét đến trong nghiên nghĩa là mô đun chất lượng nước giải quyết các cứu là các cơ sở nuôi trồng thủy sản trên địa quá trình biến đổi sinh học của các hợp chất bàn tỉnh. Dựa trên các đợt khảo sát, vị trí xả trong sông còn mô đun lan truyền chất được thải, lưu lượng xả và nồng độ tại các vị trí xả đã dùng để mô phỏng quá trình tải khuếch tán của được thu thập. Điểm xả thải của một số khu các hợp chất đó. được tập hợp thành các kênh xả đổ vào hệ 3.2. Thiết lập mạng thủy lực thống sông chính. Tải lượng các chất ô nhiễm gia nhập từ các kênh xả đó được đưa vào mô Khu vực nghiên cứu có hai hệ thống sông hình như là các biên nội của hệ thống. chính là Bến Hải và Thạch Hãn. Nối kết giữa hai hệ thống sông này là sông Cánh Hòm. 3.4. Hiệu chỉnh và kiểm nghiệm mô hình thủy lực Ngoài ra, trên địa bàn tỉnh Quảng Trị còn có sông Vĩnh Định, nối từ cống Việt Yên thuộc xã - Số liệu mực nước thực đo dùng để hiệu Triệu An chảy qua các huyện Triệu Phong, Hải chỉnh mô hình: 17/9 - 5/12/2004.
4. N.V.A. Tuấn và nnk. / Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất và Môi trường, Tập 32, Số 3S (2016) 250-255 253 - Số liệu mực nước thực đo dùng để kiểm định mô hình: 1/9 - 30/11/2009 (Hình 2, 3). Kết quả so sánh trên hình 2 -3 cho thấy giá trị tính toán từ mô hình phù hợp với giá trị thực đo cả về giá trị và pha dao động. Đường quá trình mực nước gần như trùng khít. Từ kết quả này cho thấy có thể ứng dụng mô hình cho các tính toán tiếp theo. 3.5. Thiết lập mô hình chất lượng nước Kiểm Định Hình 2. Kết quả hiệu chỉnh và kiểm định mô hình Điều kiện ban đầu: dựa vào các số liệu quan thủy lực tại trạm Đông Hà. trắc môi trường tại các vị trí trên các sông tỉnh Quảng Trị năm 2014 để thiết lập các điều kiện khởi tạo ban đầu của hệ thống và hệ số khuếch tán D (m2/s). Biên lưu lượng nhập bên: gồm 12 biên nhập bên, lưu lượng thải từ các nguồn gây ô nhiễm. Biên nồng độ: nồng độ các chất tại biên trên, biên dưới và biên nhập bên Mô hình đã chạy với thời gian từ 1/10 - 31/10. Hình 9, 10, 11 là các bản đồ mô phỏng nồng độ các chất ô nhiễm Hiệu chỉnh cao nhất trên hệ thống sông thuộc tỉnh Quảng Trị (Hình 4, 5, 6). Kiểm định Hình 3. Kết quả hiệu chỉnh và kiểm định mô hình thủy lực tại trạm Thạch Hãn. Hình 1. Sơ đồ mạng lưới tính toán thủy lực hệ thống sông Bến Hải - Thạch Hãn. Hình 4. Bản đồ phân bố nồng độ BOD trên hệ thống Hiệu chỉnh sông tỉnh Quảng Trị.
5. 254 N.V.A. Tuấn và nnk. / Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất và Môi trường, Tập 32, Số 3S (2016) 250-255 thời gian khoảng 1-2 tiếng thì lượng chất ô nhiễm sẽ trôi xa khoảng 11km trên sông Bến Hải, 5 km trên hệ thống sông Thạch Hãn, 3,5 km trên sông Sa Lung. Thời gian để quá trình tự làm sạch và pha loãng xáo trộn đưa chất lượng nước về các giá trị cân bằng ban đầu là khoảng một đến vài ngày, tùy thuộc vào các điểm xả thải (có các giá trị nồng độ và lưu lượng xả thải khác nhau). Nghiên cứu sử dụng các số liệu đo đạc trong năm 2014, tuy nhiên với sự thay đổi nhanh chóng diện tích nuôi trồng thủy sản, chất Hình 5. Bản đồ phân bố nồng độ DO trên hệ thống sông tỉnh Quảng Trị. lượng nước trong các sông cũng sẽ thay đổi rất nhanh theo chiều hướng xấu. Vì thế, sử dụng mô hình kết hợp với các số liệu đo đạc mới nhất có thể là một công cụ tốt để quy hoạch nuôi trồng thủy sản cho các địa phương trên tỉnh, dựa vào kết quả tính toán tải lượng ô nhiễm tối đa trên sông. Tài liệu tham khảo [1] Nguyễn Quang Hưng, Nguyễn Thanh Sơn, Nguyễn Vũ Anh Tuấn (2015), “Tổng quan các phương pháp xử lý có khả năng áp dụng để xử lý Hình 6. Bản đồ phân bố nồng độ tổng N trên hệ nước thải nuôi trồng thủy sản tại tỉnh Quảng thống sông tỉnh Quảng Trị. Trị”, Tạp chí Khoa học ĐHQGHN: Khoa học Tự 4. Kết luận nhiên và Công nghệ, Tập 31, Số 1S, tr. 39 - 47. [2] Nguyễn Tiền Giang và mk (2007), Đánh giá hiện Các kết quả mô phỏng cho thấy mô hình trạng ô nhiễm nguồn nước do nuôi trồng thủy Mike 11 đã đáp ứng khá tốt việc tính toán mô sản, vấn đề xâm nhập mặn tỉnh Quảng Trị và đề phỏng thủy văn, thủy lực và chất lượng nước xuất các giải pháp góp phần phát triển kinh tế xã hội và bảo vệ môi trường. trên các sông thuộc địa bàn tỉnh Quảng Trị. [3] Nguyễn Tiền Giang, Trần Ngọc Anh, Nguyễn Việc ứng dụng chạy mô phỏng cho một thời Thanh Sơn và mk (2007), Đánh giá hiện trạng ô gian dài liên tục 30 ngày giúp mô hình có tính nhiễm nguồn nước do nuôi trồng thủy sản, vấn đề ổn định cao, mô phỏng được đầy đủ các quá xâm nhập mặn tỉnh Quảng Trị và đề xuất các giải trình lan truyền chất và tự làm sạch trên sông. pháp góp phần phát triển kinh tế xã hội và bảo vệ môi trường, Trường Đại học Quốc gia Hà Nội. Từ các bản đồ trên cho thấy, hầu hết nồng [4] Sở nông nghiệp và phát triển nông thôn tỉnh độ BOD, DO, N đều nằm trong giới hạn A1. Quảng Trị (2013), Báo cáo tình hình nuôi trồng Tuy nhiên, nồng độ BOD, DO tại các vị trí ao thủy sản 2013. Kế hoạch nuôi trồng thủy sản nuôi và khu vực xung quanh các ao nuôi có xu năm 2014, Quảng Trị. hướng tăng dần và vượt mức giới hạn B1. Giá [5] Trịnh Ngọc Tuấn (2005), Nghiên cứu hiện trạng trị nồng độ tổng N cũng tăng dần tại các vị trí khai thác, nuôi trồng thủy sản ở Việt Nam và đề ao nuôi, nhưng đều nằm trong giới hạn A1, A2. xuất phương pháp xử lý nước thải. Kết quả trên mô hình cho thấy, tại các điểm [6] Environmental Protection Agency (2004), “Chapter 6: Water use, wastewater xả, môi khi nước thải NTTS được xả ra trong characterization, and pollutants of concern”.
6. N.V.A. Tuấn và nnk. / Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất và Môi trường, Tập 32, Số 3S (2016) 250-255 255 Technical development document for the final [7] Denmark Hydraulic Institute (DHI), MIKE 11 effluent limitations guidelines and new source User Guide, DHI, 2004 performance standards for the concentrated [8] Denmark Hydraulic Institute (DHI), Reference aquatic animal production point source category Manual, MIKE 11 - A modeling system for (Revised August 2004). rivers and channels, DHI, 2004. Application of MIKE 11 Model to Simulate the Spread of Contaminants caused by Aquaculture Acitivities on Some Major Rivers in Quang Tri Province Nguyen Vu Anh Tuan1, Nguyen Thanh Son2, Nguyen Quang2, Nguyen Thi Lien2 1Department of Logistic and Engineering, Ministry of Public Security, 80 Tran Quoc Hoan Street, Cau Giay, Hanoi, Vietnam 2VNU University of Science, 334 Nguyen, Thanh Xuan, Hanoi, Vietnam Abstract: By means of using MIKE 11 Model with the consideration to the impact of the discharge from aquaculture sites, this paper show the simulation results on water quality of some major rivers in Quang Tri Province. The rainfall and water level data in the year of 2005 and 2009 were used to run hydrology and hydraulic models. In addition, the monitoring data of water quality at more than 30 points as well as the measured data at the aquaculture ponds, which were taken in the year of 2013 and 2014, were used to establish water quality model and simulate the spread of contaminants on rivers. The simulation results were compared with the measured data in order to confirm the reliability of the modelling process. Some analysis were conducted in order for demostrating current state of water quality as well as making fundament to research and create an aquaculture plan in terms of environmental protection in Quang Tri Province. Keywords: Water quality, Aquaculture, MIKE 11, Quang Tri.