趙建國(guó)　李洪波　劉存歧　陳新勇　李霄宇

摘要：根據(jù)現(xiàn)場(chǎng)采集調(diào)查、實(shí)測(cè)取得河段的水體水質(zhì)和水文數(shù)據(jù)，使用QUAL2K模型模擬永定河懷來(lái)段DO、NH₄⁺-N和COD三種水質(zhì)參數(shù)，并進(jìn)行驗(yàn)證，結(jié)果表明：QUAL2K可以準(zhǔn)確模擬永定河懷來(lái)段水體水質(zhì)，模擬擬合值和實(shí)測(cè)值的相關(guān)性較好，相對(duì)誤差均在20%以內(nèi)；對(duì)永定河懷來(lái)段水質(zhì)參數(shù)模擬精度最高的是DO，其次是COD和NH₄⁺-N；官?gòu)d水庫(kù)上游退地水生態(tài)濕地保護(hù)與修復(fù)綜合治理工程建設(shè)完成后，理論上洋河張家口市八號(hào)橋控制單元DO、NH₄⁺-N能夠達(dá)到Ⅲ類水體標(biāo)準(zhǔn)。

關(guān)鍵詞：QUAL2K模型；永定河懷來(lái)段；水質(zhì) 中圖分類號(hào)：X522 文獻(xiàn)標(biāo)識(shí)碼：A

永定河屬于海河流域，洋河、桑干河在懷來(lái)縣朱官屯匯合后，形成永定河，而后懷來(lái)京西潔源污水處理廠產(chǎn)生的污水經(jīng)東沙河混入永定河，進(jìn)入官?gòu)d水庫(kù)，永定河懷來(lái)段是永定河進(jìn)入官?gòu)d水庫(kù)前最后的河段。河道每年接納大量來(lái)自上游和區(qū)域內(nèi)排放的污染物，導(dǎo)致富營(yíng)養(yǎng)化問(wèn)題嚴(yán)重，區(qū)域河流水質(zhì)下降。根據(jù)環(huán)保部印發(fā)的《國(guó)家環(huán)境保護(hù)標(biāo)準(zhǔn)“十三五”發(fā)展規(guī)劃》，洋河張家口市八號(hào)橋控制單元由Ⅳ類升為Ⅲ類。永定河是官?gòu)d水庫(kù)入庫(kù)最后一個(gè)階段，河流水質(zhì)直接影響官?gòu)d水庫(kù)水環(huán)境質(zhì)量。為改善官?gòu)d水庫(kù)上游永定河水環(huán)境質(zhì)量，在永定河大秦鐵路上游至宋家營(yíng)村南河灘地范圍內(nèi)實(shí)施退地水生態(tài)濕地保護(hù)與修復(fù)綜合治理工程，凈化永定河干流及東沙河水體。

QUAL2K模型是一個(gè)一維的綜合性河流水質(zhì)模型。模型適用于支狀河流，允許沿河有多個(gè)排污口、取水口、支流，也允許人流量有緩慢的變化，以分析人流點(diǎn)、面源負(fù)荷（包括數(shù)量質(zhì)量和位置）對(duì)受納河流水質(zhì)的影響。根據(jù)現(xiàn)場(chǎng)采集調(diào)查、實(shí)測(cè)取得河段的水體水質(zhì)和水文數(shù)據(jù)，模擬永定河懷來(lái)段水質(zhì)。

1永定河懷來(lái)段水質(zhì)優(yōu)化管理模型

1.1QUAL2K模型

QUAL2K模型是在QUAL2E模型的基礎(chǔ)上經(jīng)美國(guó)Tufts大學(xué)Steve Chapra和Greg Pelle-tier開(kāi)發(fā)的。模型基本假定為：河流斷面為梯形，污染物平流和彌散作用只在主流方向發(fā)生，水量和污染物質(zhì)量守恒。模板基本控制方程為包括源匯項(xiàng)的對(duì)流擴(kuò)散反應(yīng)方程?？刂品匠掏ㄟ^(guò)在時(shí)間和空間中隱式向后差分方法求解。模型既可以用作穩(wěn)態(tài)模型，也可以用作時(shí)變的動(dòng)態(tài)模型[3]。

QUAL2K模型基本方程為1維平流一擴(kuò)散物質(zhì)遷移方程，該方程綜合考慮了平流擴(kuò)散、稀釋、水質(zhì)組份自身反應(yīng)以及相互作用、組份的外部源和匯對(duì)組份濃度的影響。對(duì)于任意一種水質(zhì)組分，有：

式中：C為污染物濃度，mg/L；z為河流縱向坐標(biāo)，m；t為時(shí)間，s；Az為河流過(guò)水?dāng)嗝婷娣e，m²；DL為河流縱向離散系數(shù)，m²/s；u為河流斷面平均流速，m/s；S為外部的源和匯，g/s；V為計(jì)算單元的體積，m²。方程右邊的4項(xiàng)分別代表擴(kuò)散、平流、組分反應(yīng)和組分的外部源和匯。

1.2河流的分段結(jié)構(gòu)

應(yīng)用QUAL2K河流水質(zhì)模型模擬永定河懷來(lái)段河流水質(zhì)，首先將河道劃分為為18個(gè)1km的單元，由于同一河段具有相同的水利和水質(zhì)參數(shù)，將研究水系分為3個(gè)河段，其中1～3單元為第一河段，為洋河與桑干河匯合前的河段，5～12為第二河段，為洋河與桑干河匯合后但沒(méi)有和東沙河匯合前的河段，14～18為第三河段，是干流與東沙河匯合后的河段。單元及河段劃分見(jiàn)圖1。

1.3參數(shù)率定

1.3.1水力學(xué)參數(shù) 河流的流量、流速、河寬主要用于計(jì)算河流污染物在水體中的輸送。根據(jù)水文數(shù)據(jù)的實(shí)測(cè)值，采用流量系數(shù)法，模擬流量變化下的河流水文資料：

式中V為平均速率，m/s；H為水深，m；A為水流橫截面積，m²；Q為流量，m³/s；a、b、α、β為經(jīng)驗(yàn)系數(shù).b和β的和必須小于等于1，自然河道曼寧系數(shù)取值0.05。確定永定河懷來(lái)段水力學(xué)參數(shù)值，見(jiàn)表1。

1.3.2水質(zhì)參數(shù) 模型所需的水質(zhì)參數(shù)根據(jù)所選擇模擬的指標(biāo)有DO、NH₄⁺-N和COD，參數(shù)因子中復(fù)氧系數(shù)Ka選取O，Connor-Dobbins公式計(jì)算得到，耗氧系數(shù)K6的確定是利用高錳酸鉀法確定水體耗氧系數(shù)，天然水體的縱向彌散系數(shù)Exi根據(jù)QUAL2K模型提供的經(jīng)驗(yàn)系數(shù)確定，其他參數(shù)根據(jù)模型推薦值，參照實(shí)測(cè)數(shù)據(jù)適當(dāng)調(diào)整，經(jīng)過(guò)反復(fù)調(diào)試后，確定水質(zhì)參數(shù)，見(jiàn)表2。

2模擬結(jié)果分析及驗(yàn)證

QUAL2K模型的率定期選定2016年11月到2017年10月的監(jiān)測(cè)數(shù)據(jù)，對(duì)其進(jìn)行統(tǒng)計(jì)并取平均值，和項(xiàng)目監(jiān)測(cè)得到的水文數(shù)據(jù)，選取DO、NH₄⁺-N和COD作為因子對(duì)永定河懷來(lái)段水體水質(zhì)進(jìn)行模擬，結(jié)果見(jiàn)圖2。模擬結(jié)果與2017年11月水質(zhì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證，確定模型可以準(zhǔn)確模擬水質(zhì)變化情況。度沿河段呈下降趨勢(shì)，且模擬結(jié)果相對(duì)誤差最大為13.94%，模型精度較高；洋河和桑干河匯合后COD值明顯下降，說(shuō)明桑干河水質(zhì)優(yōu)于洋河，干流與東沙河匯合后，COD值顯著增高，說(shuō)明東沙河來(lái)水水質(zhì)惡劣，匯合后降低說(shuō)明，河流對(duì)污染物有一定的自凈作用，模擬結(jié)果的相對(duì)誤差平均值為0.95%；DO質(zhì)量濃度沿河段呈下降趨勢(shì)，但干流與東沙河匯合后，DO明顯上升，說(shuō)明東沙河DO值高于永定河來(lái)段干流，模擬結(jié)果相對(duì)誤差最大為-1.69%。

應(yīng)用QUAL2K模型對(duì)永定河懷來(lái)段DO、NH₄⁺-N和COD模擬結(jié)果相對(duì)誤差均在20%以內(nèi)，說(shuō)明QUAL2K模型能夠很好的對(duì)永定河懷來(lái)段進(jìn)行水質(zhì)模擬。

3優(yōu)化方案

為改善官?gòu)d水庫(kù)上游永定河水環(huán)境質(zhì)量，在永定河大秦鐵路上游至宋家營(yíng)村南河灘地范圍內(nèi)實(shí)施退地水生態(tài)濕地保護(hù)與修復(fù)綜合治理工程。項(xiàng)目區(qū)面積約467hm²，濕地治理面積333hm²，主要建設(shè)內(nèi)容包括：河道疏灘浚工程、濕地兩側(cè)防護(hù)林工程、人工濕地工程、濕地恢復(fù)工程等。項(xiàng)目生態(tài)系統(tǒng)修復(fù)后，出水水質(zhì)COD、BOD、NH₄⁺-N、總磷的處理效率均為33%。濕地的凈化效率帶入模型，水質(zhì)改善實(shí)施效果模擬，COD、NH₄⁺-N的質(zhì)量濃度到八號(hào)橋時(shí)分別為19.96mg/L和0.46mg/L，理論上洋河張家口市八號(hào)橋控制單元能夠達(dá)到Ⅲ類水體標(biāo)準(zhǔn)。

4結(jié)論

QUAL2K可以準(zhǔn)確模擬永定河懷來(lái)段水體水質(zhì)，模擬擬合值和實(shí)測(cè)值的相關(guān)性較好，相對(duì)誤差均在20%以內(nèi)。

對(duì)永定河懷來(lái)段水質(zhì)參數(shù)模擬精度最高的是DO，其次是COD和NH₄⁺-N。

官?gòu)d水庫(kù)上游退地水生態(tài)濕地保護(hù)與修復(fù)綜合治理212程建設(shè)完成后，理論上洋河張家口市八號(hào)橋控制單元DO、NH₄⁺-N能夠達(dá)到Ⅲ類水體標(biāo)準(zhǔn)。

（收稿日期：2017-12-29）