CHƯƠNG 2 ĐỐI TƯỢNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU
2.2. Phương pháp nghiên cứu
2.2.2. Phương pháp điều tra, khảo sát thực địa
Đã tiến hành khảo sát thực tế tại khu vực nghiên cứu để điều tra vị trí chính xác các nguồn thải, vị trí lấy mẫu quan trắc chất lượng nước, hệ thống thu gom
nước thải và quan sát, lựa chọn các vị trí dự kiến đặt hố ga, trạm bơm, trạm xử lý
nước thải.
Mẫu nước gồm 12 vị trí lấy mẫu nước mặt (tổng số có 14 mẫu) dọc theo tuyến kênh Than và kênh Cầu Trắng với điểm đầu là xã Ninh Hải, điểm cuối thuộc xã Hải Thanh và 4 vị trí lấy mẫu nước thải tại khu vực trung tâm thị trấn Tĩnh Gia. Thời gian lấy mẫu vào tháng 10 năm 2017. Sơ đồ vị trí lấy mẫu được thể hiện trên Hình 2.2.
Hình 2.2. Sơ đồ các vị trí lấy mẫu quan trắc mơi trường nước 2.2.3. Phương pháp lấy mẫu, bảo quản và phân tích mẫu
Các mẫu nước được lấy và bảo quản theo TCVN 5992 - 1995 và TCVN 5993 – 1995, được gửi phân tích tại Công ty TNHH 1TV Kỹ thuật Tài nguyên và
Bảng 2.2. Một số các phương pháp phân tích các chỉ tiêu chất lượng nước Chỉ tiêu Phương pháp Chỉ tiêu Phương pháp pH TCVN 6492:2011 Tổng Fe TCVN 6184:1996 Tổng Mn TCVN 6002: 1995 Độđục TCVN 6177: 1996 Amoni (NH4+) TCVN 6179-1:1996 Tổng Colifrom TCVN 6187-1:2009
Nhu cầu oxy sinh hóa (BOD5) TCVN 6001-1: 2008 Nhu cầu oxy hóa học (COD) TCVN 6491:1999
Tổng rắn lơ lứng (TSS) TCVN 6625: 2000 DO TCVN 6492:2011 P-PO43- SMEWW4500P.D-2012 N-NH4+ SMEWW5210B-2012 Thủy ngân (Hg) TCVN 7877:2008 Chì (Pb) TCVN 6193:1996 Asen (As) TCVN 6626:2000
2.2.4. Phương pháp đánh giá nhanh
Bằng cách quan sát và tính tốn các đại lượng trung bình trong trường hợp cần thiết, phương pháp này cho phép hoàn chỉnh thêm số liệu về môi trường của khu vực trong thời gian nghiên cứu như màu và mùi của nước, các điểm xả thải trên các tuyến của lưu vực kênh Than; hiện trạng hệ thống thu gom: hố ga, cống thu
gom, hướng tuyến đặt cơng trình; hiện trạng vị trí dự kiến xây dựng trạm xửlý nước thải: loại đất, cao trình vị trí, mức độ giải phóng mặt bằng…
2.2.5. Phương pháp phân tích, so sánh và đánh giá tổng hợp
Để đánh giá khu vực nghiên cứu, đã sử dụng các Quy chuẩn kỹ thuật Quốc gia vềmôi trường dưới đây:
+ QCVN 08-MT:2015/BTNMT: Quy chuẩn kỹ thuật Quốc gia về chất lượng
+ QCVN 14:2008/BTNMT: Quy chuẩn kỹ thuật Quốc gia về chất lượng nước thải sinh hoạt;
Để đề xuất biện pháp xây dựng trạm thu gom, xử lý nước thải, sử dụng các tiêu chuẩn và quy phạm dưới đây:
+ TCXD - 33/2006: Cấp nước - Mạng lưới đường ống và cơng trình – Tiêu chuẩn thiết kế; TCVN 7957 - 2008: Thốt nước – Mạng lưới và cơng trình bên ngồi – Tiêu chuẩn thiết kế;
+ Quy chuẩn xây dựng Việt Nam tập II, III ban hành theo văn bản số
439/BXD-CSXD ngày 25/9/1997;
+ QCXDVN 01:2008/BXD: Quy chuẩn Kỹ thuật quốc gia về quy hoạch xây dựng; QCVN 07:2010/BXD: Quy chuẩn kỹ thuật quốc gia về các cơng trình hạ tầng kỹ thuật đơ thị.
Các tài liệu, số liệu và kết quảđược tổng hợp theo các nội dung phù hợp của luận văn để có báo cáo hồn chỉnh.
2.2.6. Phương pháp mơ hình tốn
Với mục tiêu mơ phỏng và tính tốn chếđộ thủy văn, thủy lực và chất lượng
nước cho lưu vực kênh Than, đã lựa chọn áp dụng bộ mơ hình MIKE 11 với các mơ
đun thủy lực (HD), mô đun tải - khuếch tán (AD) và mô đun sinh thái (Ecolab). Các
đầu vào (Input), đầu ra (Output) của mơ hình được thể trong Bảng 2.3.
Bảng 2.3. Các yếu tố đầu vào, đầu ra của mơ hình MIKE 11
Đầu vào (Input) Đầu ra (Output)
Biên trên là quá trình lưu lượng thực đo của các trạm (Q ~ t);
Biên dưới là quá trình mực nước thực đo
của các trạm (H ~ t);
Biên kiểm tra là quá trình lưu lượng hoặc mực nước thực đo của các trạm trong hệ
thống;
Thành phần hạt bùn cát, kết cấu thành phần hạt (dính hay khơng dính);
Thành phần các chất dinh dưỡng (Nitơ,
Phốt pho, dầu mỡ…).
Q trình lưu lượng, mực nước tại vị
trí cần tính tốn;
Hình dạng lịng dẫn;
Các điều kiện ban đầu và điều kiện biên
Các điều kiện ban đầu
Khi sử dụng các mơ hình thủy động lực (HD model), mơ hình lan truyền, khuếch tán (AD model) và mơ hình chất lượng nước (WQ model) là: (i) Mực nước;
(ii) Lưu lượng và (iii) Hàm lượng/nồng độ các thành phần chất. Khi bắt đầu tính tốn cần phải chọn cách xác định điều kiện ban đầu. Đối với mô đun thủy lực (HD), MIKE 11 có 3 dạng điều kiện ban đầu để lựa chọn: dạng Steady-state: Hệ thống kênh, sông ở trạng thải ổn định ln có nước; dạng Prameter file: dựa vào file thông số do người sử dụng đã khai báo trong file *.hd11; dạng Hotstart: dựa vào file kết quả thủy lực đã có sẵn (lấy từ kết quả mô phỏng trước) [20].
Các điều kiện biên: Để tính tốn và mơ phỏng chất lượng nước, các điều kiện
biên được xác định cho các mô đun trong MIKE 11 như sau:
Mơ đun thủy lực (HD): Có 3 dạng điều kiện biên có thể lựa chọn để đưa vào tính tốn trong mơ đun HD:
- Biên mực nước (H) và lưu lượng (Q) không đổi theo thời gian; - Biên mực nước (H) và lưu lượng (Q) biến đổi theo thời gian; - Biên quan hệ H/Q.
Mỗi loại điều kiện biên đưa vào tính tốn dựa trên các phương trình nút khác
nhau. Việc lựa chọn điều kiện biên chô mô đun HD tùy thuộc vào hiện trạng thực tế được mơ phỏng và sự sẵn có của số liệu.
- Biên lưu lượng (Q) được sử dụng tại các vị trí: thượng lưu của sơng; dịng bên ra nhập (lateral inflow); hạlưuvà lưu lượng điều tiết (bơm).
- Biên mực nước (H) được sử dụng tại các vị trí: hạlưu của sơng và cửa sơng. - Biên quan hệ Q/H chỉ áp dụng ở hạlưu của sơng hoặc dịng ra tới hạn từ mơ hình tính tốn.
Mô đun lan truyền, khuếch tán và mơ đun chất lượng nước
Có các dạng điều kiện biên sau đây để lựa chọn sử dụng trong tính tốn vận chuyển và chuyển hóa các thành phần chất:
- Biên nồng độ/ hàm lượng mở dòng vào và dòng ra: Các biên nồng độ mở trong mô đun AD tương ứng với các biên lưu lượng, mực nước hoặc các biên Q-h
trong mô đun thủy động lực HD;
- Biên vận chuyển mở: Các điều kiện biên vận chuyển mở chỉ sử dụng tại những biên khi có dịng vào. Chuyển vận vào khu vực mơ hình được tính bằng cách sử dụng lưu lượng tại biên được tính bằng mơ đun thủy lực HD nhân với nồng độ biên được xác định theo công thức: C = Cbf Q.C = Qbf.Cbf trong đó Cbf - Nồng
độ biên (giá trịđược chỉ ra ởđầu vào mơ hình);
- Biên kín: Điều kiện biên kín là tại những điểm khơng có chất gì được vận chuyển vào và ra khỏi mơ hình tính tốn.
- Điều kiện biên nhiệt độ: Khi sử dụng mơ hình chất lượng nước kết nối với mơ hình lan truyền, khuếch tán thì trong mơ hình chất lượng nước cần thiết phải có biên nhiệt độ (T) [21].
Các thơng số và số liệu tính tốn
Các thơng số, số liệu tính tốn thủy động lực
Thông số thủy động lực cho phép người sử dụng xác định các giá trị cho một số biến trong mơ hình tính tốn thủy động lực (HD).
- Mực nước và lưu lượng ở điều kiện ban đầu: Hai thơng số này có thểđược áp dụng trên tồn bộ các nút tính tốn trong mơ hình hoặc tại một số nút trên các nhánh sông hoặc kênh dẫn.
- Hệ số lực cản đáy (Bed Resistance): Có thể lựa chọn 1 trong 3 dạng sau: (i) Hệ số Manning’n; (ii) Hệ số Manning M (M = 1/n) - hệ số này đã được mặc định; (iii) Hệ số Chezy (C = R1/6/n = M. R1/6).
- Ứng suất di chuyển của gió (wind shear stress): τw = tfac.Cw.ρa.V102. Trong
đó: Cw - Hệ số ma sát gió (3,24.10-6); tfac - Hệ số phụ thuộc địa hình; V10 - Vận tốc gió trên bề mặt (10 m); ρa - Mật độkhơng khí. Khi người sử dụng đưa thơng số dịch chuyển của gió vào tính tốn thì một điều kiện biên là biến thời gian (t) cho trường gió phải được đưa vào mô phỏng. Điều kiện biên trường gió bao gồm các mơ tả về hướng gió và vận tốc gió.
- Thơng số kiểm sốt dịng giảổn định (Quasi steady).
- Hệ số lực cản bãi ngập lũ (Flood Plain Resistance): Hệ số lực cản bãi ngập
được đưa vào bằng cách chỉnh các hệ số lực cản tương đối trong mô tả dữ liệu mặt cắt hoặc có thểđiều chỉnh hệ số lực cản bãi ngập lụt bằng cách điều chỉnh hệ số lực cản (Resistance Factor) trong dữ liệu mặt cắt đã xử lý.
- Bán kính thủy lực hoặc bán kính trở lực:
+ Bán kính trở lực (Resistance radius): √𝑅𝑅 =1𝐴𝐴∫ 𝐻𝐻0𝐵𝐵 3/2𝑑𝑑𝑏𝑏
Trong đó: A - Diện tích mặt cắt ướt (m2); B - Độ rộng mặt nước (m); H - Độ
sâu dòng chảy (water depth) m.
+ Bán kính thủy lực (Hydraulic Radius): 𝑅𝑅ℎ =𝐴𝐴𝑃𝑃 với P - chu vi ướt (m). Cần
xác định nên dùng phương pháp nào để tính bán kính thủy lực dùng tổng diện tích dịng hoặc bán kính thủy lực dùng diện tích dịng hữu hiệu).
- Mặt thủy chuẩn (datum): Việc điều chỉnh mặt thủy chuẩn sẽ được đưa vào
cho tất cả các tọa độ z và tất cả các cao trình [21].
Các thơng số, số liệu tính tốn lan truyền, khuếch tán
- Nồng độ ban đầu của các thành phần cần tính tốn.
- Hệ số phát tán dọc trục: D = f.Vex trong đó: D - hệ số phát tán (m2/s); V - vận tốc dòng chảy (m/s); f - Hệ số; ex - sốmũ khơng có thứ nguyên. Hệ số phát tán tối thiểu và tối đa được xác định sao cho không vượt quá giới hạn mà hệ số khuếch
tán được tính có thay đổi.
- Hệ số phân hủy (K): Nồng độ các thành phần phân hủy được xác định theo công thức:
𝑑𝑑𝑐𝑐
𝑑𝑑𝑡𝑡 =𝐾𝐾.𝐶𝐶
Trong đó: K - hệ số phân hủy (giờ-1); C - nồng độ thành phần phân hủy.
Các thơng số, số liệu tính tốn chất lượng nước
- Hệ số nhiệt độArrhennius (Θ(T-20)) - Nồng độ ôxy bão hòa (Cs):
- Nhiệt độ của nước (0C)
- Hằng số bổ cập ôxy ở 200C (ngày-1), K2
- Hằng số phân hủy chất hữu cơ dạng hòa tan ở 200C (ngày-1), Kd3 - Hằng số phân hủy chất hữu cơ dạng lơ lửng ở 200C (ngày-1), Ks3 - Tốc độ lắng đọng (m/ngày), K5
- Tốc độ pha loãng (resuspension) (g/m2/ngày), S1
- Hằng số phân hủy chất hữu cơ dạng lắng đọng ở 200C (ngày-1), Kb3 - Tốc độ Nitrat hóa ở 200C (ngày-1), K4
- Tốc độ phản Nitrat ở 200C (ngày-1), K6
- Hằng số tốc độ giải phóng photpho hịa tan, K8 - Hệ số hấp thụ, U3
- Năng suất quá trình quang hợp thực tế (g O2/m2/ngày), P - Tốc độ q trình hơ hấp thực tế (g O2/m2/ngày), R
- Hệ số sản lượng photpho (thể hiện lượng photpho giải phóng từ q trình phân hủy), Y3
- Hằng số tốc độ giải phóng photpho từ quá trình phân hủy sinh học, K3 - Tốc độ pha loãng (resuspension) của các hạt hữu cơ, S2
- Tốc độ lắng đọng trầm tích của các hạt hữu cơ, K6
- Tốc độ phân hủy tác động đến biến đổi Coliform hoặc faecal Coliform (ngày-1), Kd - Tốc độ phân hủy ở 200C (ngày-1), khơng tính đến độ mặn và ánh sáng, K
- Hệ số mặn, Θs
- Độ mặn, Sal - Hệ sốánh sáng, ΘI
- Cường độ ánh sáng (Kw/m2), I - Nhiệt độ của nước (0C), T
Sau khi đã có đủ các thơng số, số liệu, các điều kiện ban đầu và điều kiện biên, giải các phương trình của mơ hình sẽ cho kết quả phân bố vận tốc, nồng độ
các cấu tử theo không gian và thời gian [22]. MIKE 11 sử dụng các phương pháp và cơ sở lý thuyết của mơ hình như sau:
Mơ đun thủy lực (HD)
Mô đun thủy lực được xây dựng trên cơ sở hệ phương trình Saint Venant
một chiều cho trường hợp dịng khơng ổn định, gồm hai phương trình: + Phương trình liên tục: q t A x Q = ∂ ∂ + ∂ ∂ + Phương trình động lượng có dạng: α t Q ∂ ∂ + x ∂ ∂ (β A Q2 ) + gA x h ∂ ∂ + g RA 2 C | Q | Q = 0
trong đó: Q - Lưu lượng qua mặt cắt (m3/s); A - Diện tích mặt cắt ướt (m2); t - Thời gian tính tốn (s); α - Hệ sốđộng năng; g - Gia tốc trọng trường g= 9.81 m/s2; R - Bán kính thủy lực; q - Lưu lượng nhập lưu; x - Chiều dài theo dòng chảy (m); β - Hệ số phân bốlưu tốc; C - Hệ số Sê-di.
Hệphương trình này nói chung khơng giải được bằng phương pháp giải tích,
mà được giải bằng phương pháp số với lược đồ sai phân hữu hạn 6 điểm ẩn (Abbott - Inoescu) [21].
Các mô đun chất lượng nước
Để giải quyết vấn đề chất lượng nước có liên quan đến những phản ứng sinh hóa, mơ hình MIKE 11 sử dụng đồng thời hai mô đun là mô đun tải - khuếch tán
(AD) và mơ đun sinh thái (Ecolab) trong tính tốn.
- Mơ đun truyền tải khuếch tán được dùng để mô phỏng vận chuyển một chiều của chất huyền phù hoặc hoà tan (phân huỷ) trong các lịng dẫn hở dựa trên
phương trình để trữ tích luỹ với giả thiết các chất này được hoà tan trộn lẫn. Quá
trình này được biểu diễn qua phương trình sau:
trong đó: A - diện tích mặt cắt ngang; C - nồng độ; C2 - nồng độ nguồn; D - hệ số
khuếch tán; K - hệ số phân huỷ tuyến tính; q - dịng gia nhập; t - thời gian; x - khoảng cách. q C AKC x C AD x x QC t AC 2 + − = ∂ ∂ − ∂ ∂ − ∂ ∂ + ∂ ∂
Phương trình phản ánh cơ chế vận chuyển bình lưu/ đối lưu bởi dịng chảy trung bình và vận chuyển khuếch tán bởi gradient nồng độ [13].
- Mô đun sinh thái (Ecolab) trong mơ hình MIKE 11 giải quyết khía cạnh chất lượng nước trong kênh tại những vùng bị ảnh hưởng bởi các hoạt động dân sinh kinh tế… Mô đun này phải được đi kèm với mô đun tải – khuếch tán, điều này
có nghĩa là mơ đun chất lượng nước giải quyết các quá trình biến đổi sinh học của hợp chất trong sơng cịn mô đun tải – khuếch tán được dùng để mô tả quá trình truyền tải của các hợp chất đó.
Phần tính tốn chất lượng nước với 6 mức độ khác nhau mơ phỏng tính tốn các q trình chuyển hóa các hợp phần. Phần tính tốn sinh thái cho phép tính tới 13 thơng số sinh thái và các quá trình sinh học phức tạp. Các tác động của các yếu tố
vật lý và địa lý như nhiệt độ, bức xạ mặt trời, điều kiện ánh sáng, độ muối và vị trí
địa lý cũng được đưa vào mơ hình [16].
+ Tính tốn DO với các q trình tương tác với ơxy khí quyển trên bề mặt, các q trình hơ hấp và quang hợp của sinh vật, q trình ơxy chuyển amoni thành nitrat, nhu cầu ôxy đáy.
+ Tính tốn BOD5 thông qua các thơng số lựa chọn, có thể tính tốn được BOD5 dạng lơ lửng và dạng hòa tan trong nước, BOD5 trong lớp bùn đáy. Mơ hình
cịn cho phép tính tốn các q trình sinh hóa của BOD5 là phân rã BOD5, các q trình chuyển hóa giữa các hợp phần BOD5.
+ Tính tốn phốt pho mơ hình cho phép tính tốn hai hợp phần riêng biệt là
phốt pho vô cơ và hữu cơ, các quá trình sinh hóa xảy ra như nhận phốtpho từ phân rã BOD5, phốt pho chuyển hóa vào sinh vật, phân hủy phốt pho do vi khuẩn.
+ Tính tốn amoni với quá trình sinh ra do phân hủy BOD5, tiêu hao do chuyển hóa thành nitrat, thực vật và vi khuẩn hấp thụ.
+ Tính tốn nitrat với các q trình sinh ra do chuyển hóa amoni, suy giảm do chuyển hóa thành nitơ tự do.
+ Tính tốn các hợp phần kim loại nặng và các chất rắn lơ lửng trong nước, các hợp phần kim loại nặng được phân ra và tính tốn tồn tại trong dạng kết dính