王晨婉 李笑晉 （天津理工大學(xué)環(huán)境科學(xué)與安全工程學(xué)院 天津 300384)

張會(huì) 高丹 （山東省濟(jì)南市排水管理服務(wù)中心 山東濟(jì)南 250100)

雨水徑流污染是城市河道水環(huán)境污染的重要因素。在降雨過程中，雨水及所形成的徑流排入河流，對(duì)地表水乃至地下水造成污染。污染的發(fā)生具有隨機(jī)性、污染物排放不確定性、污染負(fù)荷的時(shí)空差異性以及監(jiān)測(cè)、控制、管理難度大等特點(diǎn)。目前，應(yīng)用數(shù)學(xué)模型研究流體水量水質(zhì)特性成為各種研究機(jī)構(gòu)和管理部門的重要方式，通過建立河道污染數(shù)學(xué)模型，了解流體的紊動(dòng)輸移規(guī)律和污染物在水體中的輸移擴(kuò)散規(guī)律，從而進(jìn)行合理的污染物濃度預(yù)測(cè)。我國(guó)在城市雨水徑流污染模型研究方面己經(jīng)取得了一定有價(jià)值的研究成果。[1，2]

本文旨在通過監(jiān)測(cè)雨水徑流及受納水體的水量水質(zhì)，基于河道模擬軟件MIKE11，模擬河道水量水質(zhì)的變化規(guī)律，探究建立排入雨水情況與受納水體之間的內(nèi)在關(guān)系。在此基礎(chǔ)上，結(jié)合雨水水質(zhì)水量，對(duì)降雨形成的污染物遷移擴(kuò)散進(jìn)行預(yù)測(cè)預(yù)報(bào)。該研究能夠?yàn)樗廴绢A(yù)測(cè)預(yù)警系統(tǒng)提供技術(shù)支撐，可提高城市主管部門應(yīng)對(duì)轄區(qū)內(nèi)水環(huán)境各類風(fēng)險(xiǎn)預(yù)警和水污染防治的科學(xué)決策能力。

1 研究河網(wǎng)概況

研究河網(wǎng)位于濟(jì)南市老城區(qū)內(nèi)，所在河系是一條南北向貫通的河道，上游為羊頭峪東溝，中游為歷山路邊溝，下游為研究河網(wǎng)。上游羊頭峪東溝起源于濟(jì)南千佛山山脈，下游匯入小清河。研究河網(wǎng)現(xiàn)狀河形比較規(guī)整，斷面為矩形河床，全長(zhǎng)為13.2km，河系匯水面積17.78km2。研究河網(wǎng)屬于典型的城市河流，它是濟(jì)南市中心城區(qū)主要排洪河道，近年來受河道侵占、棚蓋及雨污混流的影響，河水流動(dòng)性差，河道污染嚴(yán)重，行洪斷面嚴(yán)重不足，防洪標(biāo)準(zhǔn)僅為20年一遇。[3]

研究河網(wǎng)由一條主干流及左右兩條分支組成，如圖1所示，主干流長(zhǎng)約8400m，左支流長(zhǎng)約2200m，右支流長(zhǎng)約3000m。研究河網(wǎng)兩側(cè)居民區(qū)較多，排污口較多，水流動(dòng)性差，且相對(duì)集中，非常具有代表性。

圖1 研究河段的河網(wǎng)形狀

2 研究河網(wǎng)水力水質(zhì)模型（MIKE11）的建立

MIKE11是由DHI公司研究和開發(fā)的用于水體模擬的軟件，并帶有對(duì)流擴(kuò)散、水質(zhì)生態(tài)、泥沙傳輸、降雨徑流、洪水預(yù)報(bào)、實(shí)時(shí)操作及潰壩模擬等模塊。

2.1 水動(dòng)力模型

2.1.1 HD水動(dòng)力模型簡(jiǎn)介 河網(wǎng)水動(dòng)力模擬的基本目的是提供河道各個(gè)斷面、各個(gè)時(shí)刻的水位和流量等水文要素信息，并模擬泵站和閘門調(diào)度規(guī)則對(duì)河道水文條件的影響，為水質(zhì)模型提供基礎(chǔ)信息。水動(dòng)力學(xué)模型采用的是MIKE11HD模型，[4]基本原理為圣維南方程組。其差分格式采用了六點(diǎn)中心隱式格式，數(shù)值計(jì)算采用傳統(tǒng)的追趕法，即雙掃算法。計(jì)算網(wǎng)格由水位點(diǎn)和流量點(diǎn)交叉組成，[5,6]圣維南方程組為：

式中：Q為流量；A為斷面面積；q為旁側(cè)入流；h為水深；C為謝才阻力系數(shù)；R為水力半徑；α為動(dòng)量系數(shù)，一般取值為1。

2.1.2 水動(dòng)力模型的構(gòu)建 根據(jù)實(shí)測(cè)的水文要素信息進(jìn)行流域描述、河道斷面形狀定義、模型邊界的水文數(shù)據(jù)輸入及填加水工構(gòu)筑物，完成水動(dòng)力模型構(gòu)建。首先依據(jù)現(xiàn)有地形圖進(jìn)行河網(wǎng)概化，其原則是能基本反映天然河網(wǎng)的水力特性，即概化后河網(wǎng)、湖泊的輸水能力和調(diào)蓄能力與實(shí)際河網(wǎng)、湖泊相近或基本一致。[7]其次，定義河網(wǎng)斷面文件，確定河道相關(guān)參數(shù)。本文中進(jìn)行河道形狀定義主要考慮以下因素：河道形狀定義為梯形；參照資料文獻(xiàn)定義河底高程。通過河網(wǎng)斷面文件的輸入，結(jié)合河網(wǎng)文件，完成模擬河網(wǎng)的物理形狀。本文中河道形狀較規(guī)則，坡度變化不大。河流初始條件和河床糙率n，也是水流數(shù)學(xué)模型中需要確定的參數(shù)。河床糙率主要反應(yīng)了水流、泥沙、河道特性等多種因素的綜合阻力作用，[8]依據(jù)相關(guān)參考文獻(xiàn)，河網(wǎng)初始水深設(shè)為0.3m，由于河流底部比較均勻，設(shè)置糙率值（n）為0.025。再次，設(shè)定邊界條件及模擬時(shí)間。結(jié)合研究河網(wǎng)的河道斷面資料、2011年實(shí)測(cè)水文數(shù)據(jù)及相關(guān)研究資料，確定研究河網(wǎng)的水位與流量。設(shè)定模擬時(shí)間為24h，其中排放雨水前模擬時(shí)間9h，排放雨水模擬時(shí)間2h，排放雨水后模擬時(shí)間13h。時(shí)間步長(zhǎng)為1min，即每分鐘記錄一次模擬結(jié)果。到此就完成水動(dòng)力模擬文件的輸入。最后，結(jié)合模擬結(jié)果及實(shí)地監(jiān)測(cè)數(shù)據(jù)，確定并調(diào)整水動(dòng)力模型的參數(shù)值，完成水動(dòng)力模型的構(gòu)建。

2.1.3 水動(dòng)力模擬結(jié)果 水動(dòng)力模擬結(jié)果可以觀察河網(wǎng)水位流量、斷面水位等水力參數(shù)的變化情況。圖2和圖3分別為研究河網(wǎng)的水位和流量動(dòng)態(tài)模擬結(jié)果。

圖2 研究河網(wǎng)水位動(dòng)態(tài)模擬

圖3 研究河網(wǎng)流量動(dòng)態(tài)模擬

圖4為研究河網(wǎng)某一時(shí)刻的縱截面水位圖，藍(lán)綠色部分代表河道此刻的水位情況，紅色間斷線表示模擬區(qū)間水位最大值。研究河網(wǎng)流量較穩(wěn)定，因此水位模擬結(jié)果隨時(shí)間變化不大，藍(lán)綠色部分的最大值與紅色間斷線的距離很近。

圖4 研究河網(wǎng)主干流水位動(dòng)態(tài)模擬

2.2 水質(zhì)模型的構(gòu)建

2.2.1 AD水質(zhì)模型的簡(jiǎn)介 MIKE11AD可對(duì)水體中的可溶性物質(zhì)和懸浮性物質(zhì)對(duì)流擴(kuò)散過程進(jìn)行模擬，根據(jù)HD模塊產(chǎn)生的水動(dòng)力條件，應(yīng)用對(duì)流擴(kuò)散方程進(jìn)行計(jì)算。一維河流水質(zhì)模型的基本方程為：

式中：C為模擬水質(zhì)指標(biāo)濃度；D為擴(kuò)散系數(shù)；Q為流量；A為斷面面積；K為綜合衰減系數(shù)；C2為源匯項(xiàng)濃度；q為旁側(cè)入流。

水質(zhì)方程中需要確定的參數(shù)有兩個(gè)：擴(kuò)散系數(shù)D和綜合衰減系數(shù)K。

2.2.2 水質(zhì)模型建立 擴(kuò)散系數(shù)D是反應(yīng)河流縱向混合特性的重要參數(shù)，它主要受水流條件、斷面特征及河道形態(tài)等因素的影響，這里采用經(jīng)驗(yàn)值5m2/s。

污染物綜合衰減系數(shù)K是對(duì)污染物在水體中物理、化學(xué)和生化反應(yīng)過程的一個(gè)綜合描述，包括了河道自凈、沉降及吸附等復(fù)雜的反應(yīng)過程。本文中分別在黃臺(tái)南路邊溝西口南岸、柳行河中斷四號(hào)橋北河西側(cè)、南全福大街、柳行河北頭一號(hào)橋北河西側(cè)這4個(gè)站點(diǎn)，進(jìn)行表面取水的實(shí)測(cè)數(shù)據(jù)，運(yùn)用DRB200反應(yīng)器、分光光度計(jì)等儀器監(jiān)測(cè)。依據(jù)實(shí)地監(jiān)測(cè)數(shù)據(jù)，基于MIKE11軟件的數(shù)學(xué)原理，確定COD的衰減系數(shù)為0.12/d，NH3-N的衰減系數(shù)為0.08/d。

2.2.3 水質(zhì)模擬結(jié)果 將基于MIKE11模型的模擬結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)值進(jìn)行對(duì)比。從圖5、圖6可以看出：隨著研究河網(wǎng)主河道長(zhǎng)度的增加，COD值和NH3-N值整體是呈緩慢下降的趨勢(shì)，在主河道2400m和4000m左右COD值和NH3-N值驟然上升是由于左右支流匯入的緣故。模擬結(jié)果表明，模擬值與實(shí)測(cè)值基本吻合，可為河網(wǎng)水質(zhì)模型和暴雨模擬計(jì)算提供準(zhǔn)確的水質(zhì)條件。

圖5 研究河網(wǎng)的COD率定結(jié)果

圖6 研究河網(wǎng)的NH3-N率定結(jié)果

3 暴雨模擬與分析

3.1 暴雨模擬

研究河網(wǎng)上游主河道西側(cè)2400m處（膠濟(jì)鐵路南邊）設(shè)有雨水排水口，該排水口建有雨水泵站，雨水排放量較大。本文針對(duì)夏季暴雨模擬研究河網(wǎng)水量變化及污染物的遷移規(guī)律。模擬暴雨由上午9點(diǎn)開始?xì)v時(shí)2小時(shí)，至11時(shí)結(jié)束，基于MIKE 11軟件構(gòu)建的河網(wǎng)模型，預(yù)測(cè)暴雨事件造成的河道污染的危害范圍、時(shí)間與程度。

3.2 模擬結(jié)果與分析

圖7 研究河網(wǎng)主河道COD的變化情況，高線與低線分別代表河道擬時(shí)段內(nèi)COD濃度最大值與最小值，中間線為模擬時(shí)刻的COD濃度情況

圖8 研究河網(wǎng)主河道NH3-N的變化情況，高線與低線分別代表河道擬時(shí)段內(nèi)NH3-N濃度最大值與最小值，中間線為模擬時(shí)刻的NH3-N濃度情況

由圖7、圖8可以看出：在降雨開始初期，雨水及所形成的徑流流經(jīng)城市地面等沖刷、聚集了一系列污染物，研究河網(wǎng)在雨水口處污染物急劇增加并向下游擴(kuò)散達(dá)到縱坐標(biāo)最高線，見圖7（b）及圖8（b）；在11:00降雨停止后，污染物繼續(xù)向中下游擴(kuò)散，研究河網(wǎng)主河道中下游的COD、NH3-N含量持續(xù)增加，見圖7（c）與圖8（c）；隨著時(shí)間的推移和左右支流流量的匯入，由于河流本身的自凈和納污能力，研究河網(wǎng)中下游河道污染物的濃度逐漸下降，在15:00左右研究河網(wǎng)河道污染物濃度達(dá)到穩(wěn)定，之后基本保持穩(wěn)定，見圖7（d）與圖8（d）。模擬結(jié)果顯示，暴雨結(jié)束4個(gè)小時(shí)后，暴雨引起的河網(wǎng)水質(zhì)惡化影響逐漸消除，河網(wǎng)水質(zhì)恢復(fù)到暴雨前的狀態(tài)。

4 結(jié)語

本文根據(jù)研究河網(wǎng)實(shí)測(cè)資料與歷史數(shù)據(jù)，基于MIKE11軟件構(gòu)建了河流水動(dòng)力與主要水質(zhì)指標(biāo)COD、NH3-N遷移規(guī)律模型；并利用實(shí)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行率定和驗(yàn)證，模擬值和實(shí)測(cè)值吻合良好，該模型能夠準(zhǔn)確模擬研究河網(wǎng)水量水質(zhì)情況。

在模型基礎(chǔ)上，進(jìn)一步模擬城市暴雨對(duì)研究河網(wǎng)河道的影響規(guī)律，模型結(jié)果顯示：受降雨影響河道水環(huán)境質(zhì)量惡化，相關(guān)水質(zhì)指標(biāo)劇增，在位于膠濟(jì)鐵路南邊河西側(cè)的雨水口處，NH3-N值由 1.3mg/L增加到最高達(dá) 4.4mg/L，COD值由21mg/L增加到最高達(dá)46.2mg/L，由于前期雨水的污染物含量較大，因此排入河網(wǎng)雨水后，河網(wǎng)污染物的含量迅速增高；后期雨水帶來的污染物較少，此刻河網(wǎng)整體的污染物濃度逐漸降低，暴雨影響持續(xù)4個(gè)小時(shí)之后，在河道自凈作用下水質(zhì)質(zhì)量恢復(fù)正常；模擬結(jié)果與實(shí)測(cè)值相符。該模型能夠用于降雨徑流對(duì)研究河網(wǎng)水量水質(zhì)影響的預(yù)測(cè)。

城市水環(huán)境系統(tǒng)預(yù)測(cè)是城市水環(huán)境管理的重要手段和決策依據(jù)，它是利用先進(jìn)的信息技術(shù)，建立城市河道信息發(fā)布系統(tǒng)和應(yīng)急處置體系，以有效支撐城市水環(huán)境的科學(xué)保護(hù)。針對(duì)突發(fā)的暴雨，預(yù)測(cè)的時(shí)間越短，就越能提高應(yīng)急處置的主動(dòng)性和及時(shí)性，減少災(zāi)害損失。此外，河道水環(huán)境預(yù)測(cè)河道水污染控制系統(tǒng)規(guī)劃、系統(tǒng)分析、截流前期雨水方案制定及河網(wǎng)污染控制管理具有重要意義。■

[1]劉俊勇，張?jiān)疲迾浔?東江三角洲水環(huán)境綜合模型及其應(yīng)用研究[J].人民珠江，2008（6）:4-8.

[2]馬鐵民.松花江水質(zhì)模型簡(jiǎn)介[J].東北水利水電，2007，25（9）：48-51.

[3]譚振東，張忠學(xué).RMA2模型在研究河段水動(dòng)力模擬中的應(yīng)用[J].黑龍江水利科技，2010（2）：115-116.

[4]梁彬銳.MIKE11模型在沙井河片區(qū)防洪排澇工程中的應(yīng)用[J].中國(guó)農(nóng)村水利水電，2008（7）:81-83.

[5]王領(lǐng)元.丹麥MIKE11水動(dòng)力模塊在河網(wǎng)模擬計(jì)算中的應(yīng)用研究[J].中國(guó)水運(yùn)（學(xué)術(shù)版），2007，7（2）：106-107.

[6]顧玨蓉，徐祖信，林衛(wèi)青.蘇州河水系水動(dòng)力模型建立及應(yīng)用[J].上海環(huán)境科學(xué)，2002，21（10）：606-609.

[7]徐祖信，盧士強(qiáng).平原感潮河網(wǎng)水質(zhì)模型研究[J].水動(dòng)力學(xué)研究與進(jìn)展（A輯），2003，l8（2）:182-188.

[8]李吉學(xué)，汪中華.MIKE11在入河排污口設(shè)置研究中的應(yīng)用[A].中國(guó)水利學(xué)會(huì).中國(guó)水利學(xué)會(huì)2010學(xué)術(shù)年會(huì)論文集（上冊(cè)）[C].鄭州：黃河水利出版社，2010.