Điều kiện biên và điều kiện ban đầu

Một phần của tài liệu 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ố (Trang 43 - 46)

Với môđun MIKE 21 HD, tại biên lỏng phía biển là biên cửa lớn của vịnh nằm phía Tây Nam, thông với biển Khánh Hòa. Mô hình sử dụng dao động mực nước biển làm điều kiện biên. Các phần mềm dự báo mực nước hiện nay đang rất được sử dụng như như gói phần mềm Tide42 hoặc gói phần mềm đã được tích hợp sẵn trong MIKE đều cho kết quả tương đối chính xác. Tuy nhiên, để tiện lợi, tác giả sử dụng số liệu mực nước tại các biên lỏng từ kết quả dự báo thủy triều của gói phần mềm TMD (Tide Model Driver). Đây là một gói phần mềm cho phép dự báo thủy triều trên toàn cầu do viện ESR (Earth & Space Research, Mỹ), là viện nghiên cứu phi lợi nhuận về không gian và trái đất xây dựng.

Trường gió đưa vào mô hình là trường gió trung bình tháng của hướng gió thịnh hành nhất trong tháng đó.

Bảng 3.1:Tốc độ gió trung bình tháng và hướng gió thịnh hành trong tháng khu vực Cam Ranh (Số liệu được lấy từ Đài Khí tượng Thủy văn Nam Trung bộ)

Yếu tố I II III IV V VI VII XIII IX X XI XII

Tốc độ gió (m/s) 4.7 4.5 4.8 4.6 4.7 3.6 3.7 3.8 3.9 4.1 4.8 5.1

Hướng gió N NE NE SE SE SW SW SW SE NE N N

Với mô đun ECO Lab, các giá trị biên đưa vào dựa trên các kết quả phân tích mẫu trong các tháng đại diện cho hai mùa. Giá trị biên ngoài được lấy từ các số

39

liệu quan trắc tại các trạm xa bờ, nơi có thể xem không còn ảnh hưởng của các yếu tố môi trường có nguồn gốc từ lục địa xét trên quy mô mùa.

Tại thời đi ểm bắt đầu tính của mô hình , mực nước trên toàn miền tính là giá trị trung bình mực nước trên các biên tại thời điểm bắt đầu tính . Các thành phần vâ ̣n tốc dòng chảy bằng không . Giá trị nồng độ các chất ô nhiễm trên toàn bộ miền tính được nội suy từ số liệu quan trắ c ta ̣i 18 trạm phân bố trên toàn vịnh và các giá trị tại các biên (hình 3.2, hình 3.3), các kết quả này được sử dụng làm điều kiện ban đầu của các yếu tố ô nhiễm . Ngoài ra giá trị đầu vào của các yếu tố ô nhiễm môi trường tại biên là các kết quả của trạm đo liên tục tại cửa vịnh Cam Ranh của chuyến khảo sát trong tháng 5 và tháng 10/2011(bảng 3.2, bảng 3.3)

Bảng 3.2: Trung bình các giá trị tại biên của các thành phần vật chất trong mô hình tính tính toán trong mùa khô

STT Temp (°C) Sal (‰) DO (mg/l) BOD (mg/l) PO4 (mg/l) NH4 (mg/l) NO3 (mg/l)

1 33.9 29.25 5.92 0.63 24.77 864.521 4.32

Bảng 3.3: Trung bình các giá trị tại biên của các thành phần vật chất trong mô hình tính tính toán trong mùa mưa

STT Temp (°C) Sal (‰) DO (mg/l) BOD (mg/l) PO4 (mg/l) NH4 (mg/l) NO3 (mg/l)

1 31.23 27.26 5.98 0.45 32.05 27.22 0.87

Bảng 3.4: Giá trị các biến trạng thái sử dụng trong môđun ECO Lab

Stt Mô tả Kiểu Giá trị Thứ nguyên

1 Vĩ độ Biến số File đi ̣a

hình

(0), (m)

2 Quá trình BOD: Tốc độ phân rã bậc 1 tại 20oC (dạng hòa tan) Hằng số 0.161 (/ngày) 3 Quá trình BOD: Tốc đô ̣ phân rã bâ ̣c 1 tại 20oC (dạng lơ lửng) Hằng số 0.05 (/ngày) 4 Quá trình BOD: Tốc đô ̣ phân rã bâ ̣c 1 tại 20oC (dạng trầm tích) Hằng số 0.05 (/ngày) 5 Quá trình BOD: Hê ̣ số nhiê ̣t của tốc đô ̣ phân rã (dạng hòa tan) Hằng số 1.07 Phi thứ nguyên 6 Quá trình BOD: Hê ̣ số nhiê ̣t của tốc đô ̣ phân rã (dạng lơ lững) Hằng số 1.07 Phi thứ nguyên 7 Quá trình BOD: Hê ̣ số nhiê ̣t của tốc đô ̣ phân rã (dạng trầm tích) Hằng số 1.07 Phi thứ nguyên 8 Quá trình BOD: Nồng đô ̣ ôxy bán bão hòa Hằng số 2 mg/l

9 Quá trình tái lơ lửng: Vâ ̣n tốc dòng chảy tới ha ̣n Hằng số 0.3 m/s 10 Quá trình tái lơ lửng: Tỉ lệ tái lơ lửng của BOD (dạng trầm

tích)

Hằng số 0 (/ngày)

40

12 Quá tình lắng đọng: Tỉ lệ lắng đo ̣ng của BOD (dạng trầm tích) Hằng số 0.2 (/ngày) 13 Quá trình đạm hóa: Tốc đô ̣ phân rã bâ ̣c 1 ở 200C Hằng số 0.05 (/ngày) 14 Quá trình đạm hóa: Hê ̣ số nhiê ̣t của tốc đô ̣ phân rã Hằng số 1.088 Phi thứ nguyên 15 Quá trình đạm hóa: Nhu cầu ôxy cho quá trình đa ̣m hóa Hằng số 4.57 g O2/g NH4-N 16 Quá trình đạm hóa: Nồng đô ̣ ôxy bán bão hòa Hằng số 2 mg/l

17 Quá trình ôxy: Cực đa ̣i sản xuất ôxy buổi trưa, m2

Hằng số 2 (/ngày) 18 Quá trình ôxy: Độ sâu đĩa Secchi Hằng số 0.4 m 19 Quá trình ôxy: Hê ̣ số hiê ̣u chỉnh vào buổi trưa Hằng số 0 giờ 20 Quá trình ôxy: Tốc đô ̣ hô hấp của thực vâ ̣t, m2 Hằng số 0 (/ngày) 21 Quá trình ôxy: Hê ̣ số nhiê ̣t quá trình hô hấp Hằng số 1.08 Phi thứ nguyên 22 Quá trình ôxy: Nồng đô ̣ bán bão hòa trong quá trình hô hấp Hằng số 2 mg/l

23 Quá trình ôxy: Nhu cầu ôxy trầm tích trên m2 Hằng số 0.5 (/ngày) 24 Quá trình ôxy: Hê ̣ số nhiê ̣t của SOD Hằng số 1.07 Phi thứ nguyên 25 Quá trình ôxy: Nồng đô ̣ bán bão hòa của SOD Hằng số 2 mg/l

26 Quá trình ammoniac : Tốc đô ̣ giải phóng ammoni từ quá trình phân rã BOD (dạng hòa tan)

Hằng số 0.3 g NH4-N/g BOD

27 Quá trình ammoniac : Tốc đô ̣ giải phóng ammoni từ quá trình phân rã BOD (dạng lơ lửng)

Hằng số 0.3 g NH4-N/g BOD

28 Quá trình ammoniac : Tốc đô ̣ giải phóng ammoni từ quá trình phân rã BOD (dạng trầm tích)

Hằng số 0.3 g NH4-N/g BOD

29 Quá trình ammoniac: Tổng NH3-N hấp thu ̣ bởi thực vâ ̣t Hằng số 0.066 g N/g DO 30 Quá trình ammoniac: Tổng NH3-N hấp thu ̣ bởi vi khuẩn Hằng số 0.109 g N/g DO 31 Quá trình ammoniac: Nồng đô ̣ bán bão hòa của Nitơ hấp thụ Hằng số 0.05 mgN/l 32 Quá trình Nitrat: Tốc đô ̣ khử nitơ bâ ̣c nhất ở 200

C Hằng số 0.1 (/ngày) 33 Quá trình Nitrat: Hê ̣ số nhiê ̣t đô ̣ của tốc đô ̣ khử nitơ Hằng số 1.16 Phi thứ nguyên 34 Quá trình Phốtpho: Lươ ̣ng phốtpho chứa trong BOD hòa tan Hằng số 0.06 g P/g BOD 35 Quá trình Phốtpho: Lươ ̣ng phốtpho chứa trong BOD lơ lửng Hằng số 0.06 g P/g BOD 36 Quá trình Phốtpho: Lươ ̣ng phốtpho chứa trong BOD trầm tích Hằng số 0.0091 g P/g BOD 37 Quá trình Phốtpho: Lươ ̣ng PO4-P hấp thu ̣ bởi thực vâ ̣t Hằng số 0.06 g P/g DO 38 Quá trình Phốtpho: Lươ ̣ng PO4-P phân hủy bởi vi khuẩn Hằng số 0.015 g P/g DO 39 Quá trình Phốtpho: Nồng đô ̣ bán bảo hòa của phốtpho hấp thụ Hằng số 0.005 mgP/l

41

Một phần của tài liệu 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ố (Trang 43 - 46)

Tải bản đầy đủ (PDF)

(130 trang)