MÔ HÌNH CHẤT LƯỢNG NƯỚC - CHƯƠNG 3 pps

Chương 3
THIẾT LẬP VÀ LỰA CHỌN MÔ HÌNH
CHẤT LƯỢNG NƯỚC
3.1. Các bước thực hiện xây dựng mô hình chất lượng nước
3.1.1.Bước đầu
- Mục đích của việc xây dựng mô hình, vấn đề đặt ra để nghiên cứu là vấn
đề gì, tính cấp thiết của vấn đề cần nghiên cứu.
- Vấn đề nghiên cứu: cái gì sẽ chuyển biến, thay đổi và sẽ xảy ra trong dòng
chảy mà chúng ta nghiên cứu.
- Khả năng cho phép chúng ta nghiên cứu. Các yếu tố cần xác định cũng
như khả năng và điều kiện cho phép chúng ta có thể xác định được
- Các vấn đề chính cần xác định: độ nhớt, tỷ trọng, dòng chảy, hệ số chuyển
hoá hay là các quá trình phú dưỡng hoá cần xem xét Các vấn đề cần đi sâu
nghiên cứu, các vấn đề cần giới hạn.
3.1.2. Bước tiếp theo
- Cân nhắc lựa chọn mô hình để nghiên cứu : mô hình vật lý hay mô hình
số (numericl model).
Mô hình vật lý: tốn kém do việc phải xây dựng một mô hình thu nhỏ,
chi phí vận hành và hiệu chỉnh các thông số mô phỏng
điều kiện tự nhiên tốn kém và đòi hỏi các nhóm
chuyên gia có kinh nghiệm lâu năm trong lĩnh vực này.
Mô hình số: cân nhắc lựa chọn giữa mô hình 1D,2D và 3D. Khoảng
thời gian tính toán với các hệ số. Trung bình trong một
chu kỳ thời gian hay tức thời.
Mô hình 1 chiều (1D) : kinh tế thoả mãn yêu cầu. Việc tính toán đơn giản và có
thể tính tay. Giúp chúng ta có thể xác định nhanh và có thể sơ bộ xác định được
các các hệ số trong quá trình khảo sát và đo đạc thực nghiệm để hiệu chỉnh mô
hình.
Các loại mô hình này rất phổ biến do số liệu đầu vào đơn giản hơn rất nhiều.
Thích hợp cho các nghiên cứu ban đầu.
Khi vấn đề nghiên cứu trở nên phức tạp hơn rất nhiều, mô hình 1D không

đáp ứng được yêu cầu của vấn đề đặt ra thì trên cơ sở mô hình 1D có thể dễ dàng
phát triển lên 2D, 3D.
Nếu cảm thấy đủ sử dụng rồi thì có thể dừng lại ở đây. Hay nói một cách
khác mô hình một chiều là mô hình cơ sở của mô hình chất lượng nước. Đặc biệt
đối với các dòng chảy sông, khi không có sự phân vùng trong dòng chảy theo các
mặt cắt ngang. Còn theo chiều thẳng đứng ở đây thì thường có thể bỏ qua do chiều
sâu thường nhỏ hơn rất nhiều so với chiều rộng.
Các mô hình 2,3D: kết quả cho chi tiết hơn, tuy nhiên khi xây dựng thường
rất khó xác định các hệ số Khi sử dụng đòi hỏi rất nhiều các hiểu biết sâu và cụ
thể về vấn đề mà chúng ta nghiên cứu. Chưa chắc cho chúng ta các kết quả nghiên
cứu tốt hơn nhưng chắc chắn cho chúng ta hình ảnh sẽ đẹp hơn và dễ nhìn hơn
3.2. Lựa chọn mô hình chất lượng nước
3.2.1. Lựa chọn mô hình
Việc lựa chọn mô hình chất lượng nước : cần xem xét mô hình nào là thích
hợp cho đối tượng mà chúng ta nghiên cứu dựa trên cơ sở lý thuyết.
Thường vấn đề này đã có người khác lựa chọn vấn đề này giúp chúng ta
rồi. Đó là các mô hình đã sử dụng rất rộng rãi gần hơn 100 năm của lịch sử phát
triển MHCLN.
Tuy nhiên cần lưu ý : trong các mô hình chất lượng nước đã sử dụng rộng rãi thì
các hệ số, các sai số trong điều kiện thực tế của họ đã được hiệu chỉnh trong điều
kiện cụ thể của họ.
Như vậy, chúng có sự thừa nhận trong con mắt của rất nhiều người (những
người thừa nhận mô hình). Nhưng khi áp dụng vào điều kiện thực tế của chúng ta
sẽ xuất hiện các vấn đề:
-Kinh phí đầy đủ để mua toàn bộ chương trình tính toán.
-Các dữ liệu, các hệ số có thể thu thập và xác định trong điều kiện cụ thể
trên đối tượng của bạn hay không. Nếu không có đủ sẽ xảy ra vấn đề là : đối tượng
nghiên cứu của bạn một phần sẽ trong điều kiện của địa phương nghiên cứu và
một phần sẽ nằm tại địa điểm của người xây dựng chương trình này.

Qúa trình lan truyền
chất
Dữ liệu
Chuổi số
li
ệu

ứng dụng s
ơ b
ộ

Kiểm tra,so sánh
với các số liệu
Hi
ệu chỉnh

áp dụng vào tính
toán

Ki
ểm tra

Sửa
ch
ữa

Chuổi số
li
ệu

Chuổi số
liệu
Chuổi số
liệu
Mô
hìn

h
so
sán
Phát
triể
n lý
thuy
Phương
pháp s
ố

Lý
thuy
ết
Mô hình
m
ới

Tổng quát tính
toán

Dự báo chất
lư
ợng n
ư
ớc

Chất lượng nước
th
ực tế

Hình 3.1. Sơ đồ thiết lập mô hình chất lượng
nư
ớc


Vấn đề này, có thể ví dụ như một chiếc ô tô có thể là rất thích hợp, rất tốt trong
điều kiện ôn đới nhưng nó có thể trở thành không có giá trị trong điều kiện nhiệt
đới. Không có sản phẩm thương mại nào mà có thể áp dụng cho mọi đối tượng.
Như vậy khi muốn sử dụng một chương trình có sẵn cần thiết phải có hiểu
rất sâu về các giới hạn, các điều kiện biên, các hệ số của chương trình để có thể cải
tiến và hiệu chỉnh cần thiết.
Điều kiện tiên quyết ở đây là: bạn phải xây dựng và lựa chọn các hệ số
trong điều kiện cụ thể của bạn.
3.2.2. Sự phát triển lý thuyết của vấn đề
Việc đầu tiên cần phát triển lý thuyết của vấn đề của bạn với đối tượng cụ
thể mà cần xem xét và đánh giá. Phát triển lý thuyết các vấn đề: phương trình
chuyển động, phương trình liên tục và phương trình lan truyền chất.
Quan điểm của Albert Einstein: mỗi mô hình cần phải thoả mãn yêu cầu
của thực tiễn và càng đơn giản càng tốt, trừ phi không thể đơn giản hơn.
3.2.3. Kiểm định và hiệu chỉnh
Bằng các chuỗi số liệu thực đo.
Ví dụ: mô hình tính toán sự lan truyền các chất từ các nguồn thải điểm vào dòng
chảy, điều kiện chọn lựa là mùa hè. Các chuỗi số liệu được kiểm định, các hệ số sẽ
thay đổi so với điều kiện các mùa khác sẽ không cho chúng ta giá trị xấp xĩ tốt để
hiệu chỉnh.
Sau khi sử dụng các chuỗi số liệu để hiệu chỉnh các hệ số sẽ thay đổi. Khi
hiệu chỉnh cần lựa chọn một vài hệ số ổn định để hiệu chỉnh các hệ số khác và sau
đó tiến hành hiệu chỉnh các hệ số cố định hoặc hiệu chỉnh tất cả các hệ số đồng
thời.

Các điều kiện cần làm rõ khi hiệu chỉnh:
- Điều kiên biên và tải trọng.
- Điều kiện ban đầu.
Các quá trình đặc trưng cho đối tượng (quá trình mà chúng ta cần mô tả và làm
rõ). Các chất hữu cơ : quá trình vật lý, hoá học hay sinh học trong dòng chảy.
Các hệ số cần kiểm định: động học, hệ số, các hệ số đặc trưng cho quá trình
chuyển hoá các chất trong dòng chảy.
3.3. Các phương pháp số tính toán sự lan truyền các chất ô nhiễm trong
dòng chảy
3.3.1.Các phương pháp số trong nghiên cứu mô hình thủy lực
Trong tự nhiên dòng chảy thường không ổn định, ở đó các yếu tố thủy lực
tại một mặt cắt thay đổi theo thời gian (chuyển động không dừng). Chuyển động
không ổn định được chia thành hai loại: chuyển động không ổn định thay đổi gấp,
và chuyển động không ổn định thay đổi chậm dần.
Phương trình cơ bản của dòng chảy là tập hợp hệ thống hai phương trình
đạo hàm riêng phi tuyến có các hệ số biến đổi dạng Hyperbolic. Phương pháp số
để giải loại phương trình này gặp một số trở ngại trong việc đặt các điều kiện biên.
Việc tích phân hệ phương trình này gặp nhiều khó khăn, chỉ có thể tìm
được nghiệm trong một vài trường hợp đặc biệt như kênh có tiết diẹn hình chử
nhật, đáy nằm ngang và bỏ qua sức cản. Các điều kiện này rất khác xa so với điều
kiện thực tế.
Hiện nay, các phương pháp số được sử dụng để tính gần đúng cho dòng ổn
định thay đổi chậm dần bao gồm : phương pháp giải tích, phương pháp đường đặc
trưng và phương pháp sai phân hữu hạn.
Phương pháp tích phân toán học chặt chẽ được sử dụng để tìm nghiệm giải
tích của hệ phương trình trên. Theo phương pháp này, cần có những giả thiết để
đưa hệ phương trình Saint-Vennnant về dạng đơn giản, rồi tích phân các phương
trình này với các giả thiết là mặt cắt sông hình chử nhật hay lăng trụ, độ dốc bằng
không hoặc bằng hằng số, độ nhám cố định, bỏ qua sức cản
Nhờ các giả thiết trên, các phương trình vi phân đạo hàm riêng phi tuyến được

biến đổi thành phương trình vi phân đạo hàm riêng tuyến tính. Để giải được các
bài toán dạng không ổn định.
Phương pháp sai phân hữu hạn
Phương pháp sai phân là biến đổi từ hệ phương trình vi phân đạo hàm riêng
Saint-Vennant thành hệ phương trình đại số phi tuyến do tính chất phi tuyến của
hệ ơng trình Saint-Venannt. Sau đó thực hiện việc tuyến tính hóa hệ phương trình
đại số phi tuyến và giải hệ phương trình thu nhận được các giá trị về các yếu tố
thủy lực cần biết trên các nút lưới tính toán.
Phương pháp sai phân để giải bài toán thủy động lực một chiều đã được các
tác giả nghiên cứu, như J.A.Cung, C.Lai, V.Bellos, F.Ionescu với các phương
pháp sai phân theo các sơ đồ khác nhau.
Có hai phương pháp sai phân : phương pháp sai phân theo sơ đồ ẩn và
phương pháp sai phân theo sơ đồ hiện. Mỗi một phương pháp có những mặt mạnh
khác nhau.
Phương pháp sai phân sơ đồ hiện
Nếu sau phép sai phân hệ hai phương trình đại số với hai ẩn số và có thể
giải ngay được hai ẩn số đó hoặc giá trị của hàm ẩn tại các nút lưới được xác định
một cách riêng lẽ không cần các phương trình tại các nút lưới khác thì sơ đồ sai
phân được gọi là sơ đồ hiện. Các sơ đồ sai phân : tam giác cân ngược, tam giác
cân thuận, sơ đồ hình thoi trung tâm ưu điểm của sơ đồ hiện là sự đơn giản của
việc tính toán và lập chương trình. Nhưng lạ không cho chúng ta tính toán với các
bước thời gian lớn, vì để sơ đồ sai phân ổn định giá trị bước thời gian bị hạn chế
bởi điều kiện CFL(Courant-Friedric-Levi) dẫn tới việc lựa chọn thời gian tính toán
bé không thích hợp cho các bài toán yêu cầu tính toán với các chu kỳ dài.
Phương pháp sai phân sơ đồ ẩn
Trong trường hợp sau khi sai phân, số hàm ẩn lớn hơn hai thì phải áp dụng
sai phân cho một lớp các điểm rời rạc. Cùng với các điều kiện biên sau khi tuyến
tính hóa ta được hệ phương trình đại số tuyến tính đóng kín. Giải hệ phương trình
đại số này ta có nghiệm ở một loạt các điểm mà chúng ta cần tính toán. Với sơ đồ
ẩn việc tính toán các giá trị hàm ẩn tại một nút nào đó đòi hỏi phái giải đồng thời

một hệ phương trình của tất cả các nút ở cùng một thời điểm, còn bước thời gian
được chọn theo quan điểm của độ chính xác chứ không phải theo tính ổn định. Các
sơ đồ ẩn trở nên rất thuận tiện cho việc tính toán dòng chảy không ổn định trong
mạng kênh sông có lòng dẫn phức tạp, nhược điểm là việc giải hệ phương trình tại
mỗi bước thời gian cần phải tính lặp, có khối lượng tính toán lớn, đòi hỏi bộ nhớ
và thời gian tính toán lớn hơn.
Với sơ đồ ẩn việc tính toán các giá trị hàm ẩn tại một nút nào đó đòi hỏi
phái giải đồng thời một hệ phương trình của tất cả các nút ở cùng một thời điểm,
còn bước thời gian được chọn theo quan điểm của độ chính xác chứ không phải
theo tính ổn định. Các sơ đồ ẩn trở nên rất thuận tiện cho việc tính toán dòng chảy
không ổn định trong mạng kênh sông có lòng dẫn phức tạp, nhược điểm là việc
giải hệ phương trình tại mỗi bước thời gian cần phải tính lặp, có khối lượng tính
toán lớn, đòi hỏi bộ nhớ và thời gian tính toán lớn hơn.
3.3.2. Phương pháp số giải bài toán lan truyền chất
Các phương pháp giải tích hầu như không thể áp dụng để tìm nghiệm tổng
quát của bài toán Tải-khuếch tán trong trường hợp tổng quát. Chỉ được áp dụng
trong trường hợp với giả thiết dòng chảy với chế độ ổn định Cách giải phổ biến
là các phương pháp số : phương pháp sai phân hữu hạn, phương pháp phần tử hữu
hạn. Hai phương pháp này có những ưu và hạn chế khác nhau.
Phương trình tải-khuếch tán là phương trình dạng Parabolic, được giải bằng một
số phương pháp số quen thuộc, như phương pháp sai phân hữu hạn, phương pháp
phần tử hữu hạn, phương pháp đường đặc trưng. Các phương pháp sai phân cho
bài toán loại này : phương pháp Lax-Richardson, phương pháp Upwin, phương
pháp Leonard, phương pháp Duford-Frankel, phương pháp Lax-Wendroff, phương
pháp Crank-Nicolson, sơ đồ Brian-Stone, phương pháp sai phân theo hướng,
phương pháp sai phân trung tâm.
3.4.Các mô hình BOD & D (DO) trong dòng chảy
3.4.1. Phương trình cơ bản
Phương trình cơ bản mô tả sự lan truyền và phân bố vật chất trong dòng
chảy được gọi là phương trình tải-phân tán vật chất trong dòng chảy. Trên cơ sở

định luật bảo toàn khối lượng ta có :

)(sF
x
C
E
zy
C
E
yx
C
E
xx
C
u
x
C
u
x
C
u
t
C
xyxzyx














































(1.1)

Trong đó :
C -Nồng độ các chất trong dòng chảy, mg/l.
u
x
,u
y
,u
z
-Thành phần vận tốc trong dòng chảy, m/s.
E
x
, E
y
, E
z
-Hệ số phân tán rối vật chất trong dòng chảy, m
2
/s.
F(s) -Hàm số mô tả sự thay đổi các chất trong dòng chảy do
các quá trình vật lý, hóa học và sinh học.

Với giả thiết, bỏ qua sự tải do vận tốc dòng chảy, quá trình khuếch tán theo
các phơng oy, oz từ phương trình (1.1) ta có phơng trình mô tả quá trình lan
truyền-khuếch tán các chất hữu cơ tại một điểm trong dòng chảy.
)(sF
x
C
E
xx
C
u
t
C
xx
















(1.2)

Trong các nghiên cứu ban đầu cũng như ứng dụng các MHCLN cho các
dòng chảy sông rộng, cửa sông thường đợc xem xét trong không gian một chiều.
3.4.2. Phương trình cổ điển Streeter-Phelps
Các nghiên cứu về sự lan truyền, chuyển hóa các chất hữu cơ trong dòng
chảy đều dựa trên cơ sở phơng trình cổ điển Streeter-Phelps với các giả thiết sau:
-Tốc độ phân huỷ các chất hữu cơ tuân theo quy luật phản ứng bậc nhất.
-Sự thiếu hụt oxy hoà tan (D) trong dòng chảy do sự phân huỷ các chất hữu
cơ.
-Chế độ lan truyền chất trong dòng chảy ổn định.
Năm 1925, Streeter và Phelps thiết lập phương trình toán học mô tả quá
trình tiêu thụ oxy do sự phân hủy các chất hữu cơ trong dòng chảy.
Phương trình Streeter-Phelps :
u
dL
dx
k L
d
  (1.3a)
DkLk
dx
dD
u
ad
 (1.3b)
Trong điều kiện ổn định, giải hệ phương trình (1.3a,b) bằng phương pháp
tích phân cho kết quả về phân bố nồng độ các chất hữu cơ và độ thiếu hụt oxy hoà
tan (BOD &D) ở vùng hạ lưu nguồn thải.

u
xk

o
d
eLL

 (1.4a)










u
xk
u
xk
da
od
u
xk
o
ada
ee
kk
Lk
eDD (1.4b)
Trong đó :

u -Vận tốc trung bình của dòng chảy, m/s.
D
0
- Độ thiếu hụt oxy, mg/l.Tại điểm có toạ độ x=0 (điểm xáo
trộn nước thải với nước sông).
D - Độ thiếu hụt oxy, mg/l.
L
0
- Nồng độ các chất hữu cơ theo BOD, mg/l. Tại điểm có tọa
độ x=0 (điểm xáo trộn nước thải với nước sông).
L - Nồng độ các chất hữu cơ theo BOD, mg/l.
k
d
- Hằng số tốc độ tiêu thụ oxy do sự phân huỷ hợp chất hữu
cơ, ngày
-1
.
k
a
- Hằng số tốc độ hòa tan oxy qua bề mặt thoáng, ngày
-1
.
Bài tập :
1. tính toán, vẽ đường cong BOD và DO hoặc D trong dòng chảy ổn định.
3.4.3. Các nghiên cứu phát triển trên cơ sở phương trình Streeter-Phelps
Mô hình Streeter-Phelps (1.3a,b &1.4ab) chỉ kể đến các yếu tố tiêu thụ oxy
do quá trình hô hấp của vi khuẩn hiếu khí trong tầng nước và hòa tan oxy qua bề
mặt thoáng khi trong nước có sự thiếu hụt oxy hòa tan. Khi áp dụng cho các đối
tượng khác nhau cho thấy, độ tin cậy của mô hình còn thấp do chưa đề cập đến các
yếu tố ảnh hưởng đến phân bố BOD & DO trong dòng chảy như : quá trình quang

hợp và hô hấp của hệ thực vật nước, hô hấp trong tầng cặn đáy, các quá trình
ammôn hoá, nitơrit, nitơrat cũng như ảnh hưởng của dòng chảy đến sự lan
truyền các chất trong dòng chảy như : quá trình lắng các chất lơ lửng, sự phân tán
các chất do dòng chảy rối, ảnh hưởng của độ mặn và gió
Nhằm mục đích hoàn thiện hơn, nâng cao độ tin cậy của mô hình Streeter-
Phelps đối với các dòng chảy trong tự nhiên, các tác giả O’Connor, Dobbin’s,
Thomat, Camps, Di Toro bằng các nghiên cứu thực nghiệm bổ sung thêm các
quá trình sinh thái-chất lượng nước nhằm mục đích mô tả một cách tổng quát, đầy
đủ các yếu tố có liên quan đến phân bố, thay đổi BOD&DO trên dòng chảy.
Qúa trình hô hấp của lớp bùn đáy
Giả thiết dòng chảy có lớp nước nông với sự tích luỹ nhiều các chất hữu cơ
ở đáy dòng chảy, từ các phương trình (1.3a,b) ta có:
H
S
Dk
x
D
u
t
D
a






(1.5)
Phân bố BOD & D trong điều kiện ổn định là (1.4) và:























u
xk
Hk
S
u
xk
DD
a
a

a
o
exp1exp
(1.6)
Trong đó :
S

- Nhu cầu oxy của lớp cặn đáy, mgO/m
2
.ngày.
H - Chiều sâu lớp nước, m.
Qúa trình oxyhoá các hợp chất chứa nitơ
Bổ sung thêm quá trình tiêu thụ oxy do quá trình oxy hoá các hợp chất hữu
cơ chứa nitơ (quá trình nitrat hoá), từ (1.3a,b) ta có:
NkDk
dx
dD
u
na
 (4.5)
Độ thiếu hụt oxy D được xác định theo:
































u
xk
u
xk
kk
Nk
u
xk

DD
an
na
ona
o
expexpexp
(4.6)
Trong đó :
N - Nồng độ các hợp chất chứa Nitơ theo NBOD, mg/l
k
n
-Hằng số tốc độ tiêu thụ oxy do quá trình nitrat hóa, ngày
-1
.
Qúa trình quang hợp-hô hấp
Các phương trình toán học mô tả sự thiếu hụt oxy do sự phân huỷ các chất
hữu cơ và quá trình quang hợp-hô hấp của hệ thực vật nước.
tba
PRDk
x
D
u
t
D
)( 






(1.9)
Sự thiếu hụt oxy trong điều kiện ổn định :
























u
xk
k
PR

u
xk
DD
a
a
a
o
exp1exp

Trong đó:
P - Sản phẩm sơ cấp của quá trình quang hợp, mg/l.ngày
R - Lượng oxy tiêu hao do quá trình hô hấp, mg/l.ngày
Cho ví dụ tính toán
Qúa trình lắng các chất lơ lửng
Gỉa thiết, trong quá trình lan truyền một phần các chất hữu cơ dạng lơ lững,
phân tán lắng đọng, trong trờng hợp này, sự lan truyền các chất hữu cơ đợc mô tả
theo các phương trình sau:
Lk
dx
dL
u
r
 với k
r
= k
d
+ k
s

(3.7)

DkLk
dx
dD
u
ad
 (3.8)

Tích phân phương trình (3.7) và thế vào (3.8) ta có:








u
xk
LL
r
O
exp (3.9)






























u
xk
u
xk
kk
Lk
u
xk
DD

ar
ra
Oda
o
expexpexp
(3.10)
Trong đó:
k
r
-Hệ số tốc độ chuyển hoá các chất hữu cơ trong dòng chảy,
ngày
-1
.
k
s
-Hệ số tốc độ thay đổi nồng độ các chất hữu cơ do các quá
trình lắng đọng, hấp thụ các chất hữu cơ trong dòng chảy,
ngày
-1
.

Qúa trình phân tán vật chất trong dòng chảy
Các sông rộng và vùng cửa sông với sự ảnh hởng của gió, dòng triều có chế
độ xáo trộn hết sức phức tạp việc mô phỏng cần xem xét và đánh giá đồng thời
ảnh hởng của các quá trình chính lên sự phân tán các chất ô nhiễm trong dòng
chảy. Các quá trình chính bao gồm : sự phân tán rối vật chất do chế độ xáo trộn
phức tạp. Các phương trình cơ bản của quá trình :
Lk
x
x

E
x
L
u
t
L
r









2
2
(1.10)
DkLk
x
D
E
x
D
u
t
D
ar










2
2
(1.11)
Phân bố BOD & D trong điều kiện ổn định :




















2
2
4
11
2
exp
4
1
u
Ek
E
ux
u
Ek
Q
W
L
d
d


 
   























2
2
1
1
1
2
exp
1
1
2
exp
1
m
E

ux
m
m
E
ux
mQkk
Wk
D
da
d


Với
2
1
4
1
u
Ek
m
d
 ;
2
2
4
1
u
Ek
m
a


Trong đó:
Q

- Lưu lượng dòng chảy, m
3
/s.
W - Tải trọng thải của các chất ô nhiễm từ nguồn thải, mg/s.
E - Hệ số phân tán, m
2
/s.
Cho ví dụ tính toán
Phương trình tổng quát
Mô hình tổng quát cho dòng chảy kênh, sông, sông rộng và cửa sông.

L
A
Q
dt
dL
x
L
E
x
L
A
Q
t
L
w










2
2
(1.4)
DO
A
Q
dt
dDO
x
DO
E
x
DO
A
Q
t
DO
w










2
2
(1.5)

Với :

 
Lkk
dt
dL
sd


   
PR
H
S
NkLkDODOk
dt
dDO
ndBHa

Trong đó :
Q - Lưu lợng dòng chảy sông, m

3
/s.
Q
w
- Lưu lợng dòng chảy gia nhập có chứa chất ô nhiễm, m
3
/s.
L - Nồng độ các chất hữu cơ theo BOD, mg/l.
k
r
- Hệ số tốc độ chuyển hóa các chất hữu cơ trên đoạn sông,
ngày
-1
.
k
a
- Hằng số tốc độ phân hủy các chất hữu cơ, ngày
-1
.
k
s
- Hệ số tốc độ lắng các chất hữu cơ trong dòng chảy, ngày
-1
.
k
d
- Hằng số tốc độ phân hủy các chất hữu cơ, ngày
-1
.
DO

BH
- Nồng độ oxy hòa tan ở trạng thái bão hòa, mg/l
DO - Nồng độ oxy hòa tan, mg/l.
k
n
- Hằng số tốc độ phân hủy các hợp chất chứa Nitơ ngày
-1

N - Nồng độ các hợp chất chứa Nitơ theo NBOD, mg/l
P - Sản phẩm sơ cấp của quá trình quang hợp, mg/l.ngày
R - Lợng oxy tiêu hao do quá trình hô hấp, mg/l.ngày
S - Nhu cầu oxy của lớp cặn đáy, mg/m.ngày.
H - Chiều sâu trung bình của dòng chảy, m.
E
x
- Hệ số phân tán dọc theo chiều dòng chảy, m
2
/s
Các phương trình (1.4), (1.5) mô tả một cách tổng quát và phản ánh một
cách đầy đủ nhất các quá trình có liên quan đến sự thay đổi của BOD&DO trong
dòng chảy. Từ các phương trình tổng quát trên, tùy thuộc từng hoàn cảnh cụ thể,
mục đích nghiên cứu mà lựa chọn mô hình thích hợp đối với các đối tợng nghiên
cứu khác nhau. Việc lựa chọn mô hình cho một đối tợng nghiên cứu cụ thể, dựa
trên cơ sở đặc điểm : chế độ thủy lực dòng chảy, các quá trình sinh thái chất lợng
nớc diễn ra trong dòng chảy. Các mô hình chất lợng nớc cho các dòng chảy khác
nhau bao gồm : kênh, mương thoát nớc thải đô thị, dòng chảy suối, sông nhỏ, sông
rộng và cửa sông.