李吉學(xué) 紀(jì)杰善 田在豐

（1山東省濟寧水文水資源勘測局 濟寧 272019 2山東省濟寧市任城區(qū)水務(wù)局 濟寧 272073）

MIKE11水動力模塊在入河排污口設(shè)置計算中的應(yīng)用研究

李吉學(xué)1紀(jì)杰善1田在豐2

（1山東省濟寧水文水資源勘測局 濟寧 272019 2山東省濟寧市任城區(qū)水務(wù)局 濟寧 272073）

本文以王樓礦為研究區(qū)，在全面了解項目區(qū)取用水和污水排放情況的基礎(chǔ)上，采用DHI的MIKE11一維河道、河網(wǎng)綜合模擬軟件建立了排污河段的水流、水質(zhì)數(shù)值模型，分析了入河排污口設(shè)置后污水排放對水功能區(qū)的影響范圍，對水質(zhì)、水生態(tài)、第三者權(quán)益的影響以及入河排污口設(shè)置的合理性，提出了入河排污口設(shè)置的位置和滿足水功能區(qū)水質(zhì)和南水北調(diào)沿線水質(zhì)要求的排放措施，為確保南水北調(diào)工程水質(zhì)安全，保護(hù)水環(huán)境，促進(jìn)水資源的可持續(xù)利用提供了科學(xué)依據(jù)。

水文 入河排污口 設(shè)置 水質(zhì)模型

1 問題的提出

開展排污口設(shè)置及環(huán)境影響的技術(shù)研究是為更好地保護(hù)水資源，改善水環(huán)境，為促進(jìn)水資源可持續(xù)利用提供技術(shù)支撐，也是貫徹落實水利部22號令《入河排污口監(jiān)督管理辦法》和水利部水資源〔2005〕79號文《關(guān)于加強入河排污口監(jiān)督管理工作的通知》精神的重要體現(xiàn)。要按照水功能區(qū)的需求確定保護(hù)目標(biāo)，再核定水域納污能力，提出限制排污總量的意見，管理入河排污口，逐步使水域達(dá)標(biāo)。

本文以王樓礦井為項目區(qū)，利用丹麥DHI公司的水環(huán)境數(shù)值模擬軟件MIKE 11建立了模擬河段的水流水質(zhì)數(shù)值模擬系統(tǒng)，開展排污口設(shè)置及環(huán)境影響的技術(shù)研究。為各級水行政主管部門審批入河排污口設(shè)置方案以及建設(shè)單位合理設(shè)置入河排污口提供科學(xué)依據(jù)。

2 水質(zhì)模型（MIKE 11）建立

水質(zhì)模型提供一種預(yù)測自然過程和人類活動對河流水庫系統(tǒng)中水的物理特性、化學(xué)特性和生物特性影響的手段，廣泛應(yīng)用于評價來自廢水處理的廢水負(fù)荷的影響，或者來自其他各種點源和非點源污染物負(fù)荷的影響。

通過深入分析研究區(qū)域水文氣象、自然地理、河道和水流特性以及水污染的特點，應(yīng)用DHI公司的水環(huán)境數(shù)值模擬軟件MIKE 11來建立模擬河段的一維水流水質(zhì)數(shù)值模擬系統(tǒng)。結(jié)合實測河道大斷面資料及2004年、2005年實測水文、水質(zhì)數(shù)據(jù)。應(yīng)用模型參數(shù)自動優(yōu)選方法率定水動力模型和水質(zhì)模型的參數(shù)值，進(jìn)行模型的率定和驗證工作。

2.1 MIKE11軟件簡介

目前，在世界上處理地表水問題的專業(yè)軟件中，Mike 11是最受歡迎的一個河流模擬系統(tǒng)。它是DHI20多年經(jīng)驗與開發(fā)的結(jié)晶。以求解圣維南方程組作為理論基礎(chǔ)，并帶有對流擴散、水質(zhì)生態(tài)、泥沙傳輸、降雨徑流、洪水預(yù)報、實時操作及潰壩模擬等多種模塊。主要用于河口、河流、灌溉系統(tǒng)和其他內(nèi)陸水域的水文學(xué)、水力學(xué)、水質(zhì)和泥沙傳輸模擬，在防汛洪水預(yù)報、水資源水量水質(zhì)管理、水利工程規(guī)劃設(shè)計論證均得到廣泛應(yīng)用。

主要模型有：水動力學(xué)模型（HD model）、對流擴散及粘性輸沙模型、非粘性沙傳導(dǎo)模型、NAM降雨徑流模型（NAMmodel）、單位線模型、洪水實時預(yù)報模型（FF）、地理信息系統(tǒng)。

2.2 建立模型

（1）水動力模型。水動力模型所用的描述一維非恒定水流運動規(guī)律的基本方程為圣維南方程組，其數(shù)學(xué)表達(dá)式為：

水質(zhì)方程中需要確定的參數(shù)有兩個：擴散系數(shù)D和綜合衰減系數(shù)K。擴散系數(shù)是反應(yīng)河流縱向混合特性的重要參數(shù)，它主要受水流條件、斷面特征及河道形態(tài)等因素的影響。污染物綜合衰減系數(shù)是對污染物在水體中物理、化學(xué)和生化反應(yīng)過程的一個綜合描述，包括了河道自凈、沉降、吸附等復(fù)雜的反應(yīng)過程。

2.3 水流水質(zhì)模擬模型計算單元概化和模型輸入

一維水流、水質(zhì)基本方程加上實際河流的邊界情況，即構(gòu)成了研究河段的數(shù)值模擬系統(tǒng)。對于此次研究的洙趙新河梁山閘以下至南陽湖入口河段，總長度為24km，考慮到在入湖口上游500m處建有劉官屯控制閘，人為調(diào)控著洙趙新河與南陽湖的水流連通關(guān)系，建立的一維水流水質(zhì)數(shù)值模擬模型的模擬范圍為梁山閘下至劉官屯閘位置，總模擬河段長23.5km。利用收集到的26個實測大斷面資料（見圖1、圖2），對模擬河段進(jìn)行水下地形的數(shù)值概化，河道縱向空間計算步長為500m；采用有限差分法，對水流、水質(zhì)基本方程進(jìn)行離散求解。

上邊界：采用梁山閘處實測水質(zhì)數(shù)據(jù)和梁山閘下泄流量。

糙率n也是水流數(shù)學(xué)模型中唯一需要確定的參數(shù)，它主要反映了水流、泥沙、河道特性等多種因素的綜合阻力作用。

（2）水質(zhì)模擬。描述物質(zhì)在水體中輸運的一維非恒定流對流擴散基本方程為：

3 模型參數(shù)率定

3.1 模型參數(shù)率定

使用2005年3月1日～11月30日實測水流資料和水質(zhì)資料進(jìn)行模型參數(shù)的率定工作。

通過模型試算，對實測數(shù)據(jù)和模型中的計算值進(jìn)行對比分析，來率定水質(zhì)模型中的參數(shù)，使兩者擬合最好，（見圖3～圖6）。通過參數(shù)率定，最終確定整個模擬河段內(nèi)糙率取值為n=0.028，擴散系數(shù)D=525m2/s，BOD5綜合衰減系數(shù)K=0.003/h，CODCr綜合衰減系數(shù) K=0.0082/h。

3.2 模型驗證

驗證計算是模型應(yīng)用前必須做的工作，它是驗證模型可用性的重要環(huán)節(jié)。利用率定好的水流、水質(zhì)模型參數(shù)，進(jìn)行2004年水流、水質(zhì)的數(shù)值模擬，并對比分析模型計算值與實測值。結(jié)合模型率定成果和驗證情況，在現(xiàn)有資料基礎(chǔ)上，可以認(rèn)為水質(zhì)模型的參數(shù)取值還是比較合理，可以應(yīng)用于水質(zhì)預(yù)測計算。

4 水功能區(qū)水質(zhì)數(shù)學(xué)模擬分析

4.1 水質(zhì)預(yù)測設(shè)定條件

預(yù)測河段上邊界來水流量受梁山閘調(diào)度控制，采用梁山閘1980～2005年26年實測資料進(jìn)行頻率計算，求得頻率90%的枯水年放水量3535萬m3，并選用最不利典型年

分別按洙趙新河現(xiàn)狀水質(zhì)條件下和規(guī)劃水質(zhì)條件下的水質(zhì)，進(jìn)行預(yù)測分析，現(xiàn)狀水質(zhì)采用2004年和2005年實測水質(zhì)的平均值，規(guī)劃水質(zhì)為地表水Ⅲ類水質(zhì)標(biāo)準(zhǔn)限值。

4.2 水質(zhì)預(yù)測方法

考慮南水北調(diào)工程要求，分為調(diào)水期11月~翌年5月和汛期6~10月兩個水期。預(yù)測河段受梁山閘和劉官屯閘兩個節(jié)制閘運用的影響，為受人控制河道，閘門調(diào)度使用對該河段水量影響較大。

在頻率90%枯水年中11月、12月和3月出現(xiàn)了上游梁山閘關(guān)閉沒有下泄水量的情況，同時考慮到冬季1月、2月存在的河道封冰情況，對調(diào)水期間11月~翌年3月采用完全混合水質(zhì)模型來進(jìn)行水質(zhì)預(yù)測計算，從偏安全角度出發(fā)，不考慮污染物降解作用。

調(diào)水期4月、5月，頻率90%枯水年中梁山閘有少量放水，同時存在農(nóng)業(yè)灌溉用水，模擬河段呈現(xiàn)一定流動狀態(tài)，采用一維非恒定流水質(zhì)模型進(jìn)行水質(zhì)預(yù)測計算。

汛期6～10月劉官屯閘開啟泄洪放水，洙趙新河與南陽湖連通，一維流動特征比較明顯，采用一維非恒定流水質(zhì)模型進(jìn)行水質(zhì)預(yù)測計算。

排污條件：排污口作為源項處理，考慮排污口布設(shè)在距離入湖口 1km、2km、3km、5km、10km、15km、20km等 7種情況；廢污水為連續(xù)均勻排放，排放濃度考慮經(jīng)處理后正常排放濃度和異常未處理污水濃度兩種情況。

4.3 計算工況

根據(jù)南水北調(diào)工程要求，水質(zhì)預(yù)測計算中分為調(diào)水期和汛期兩個水期。按照采用的水質(zhì)預(yù)測計算方法不同，將調(diào)水期分為11月～翌年3月、4～5月兩個時期。在不同計算水期水質(zhì)預(yù)測中，分別考慮河道本底水質(zhì)為現(xiàn)狀水質(zhì)和規(guī)劃水質(zhì)兩種狀態(tài)，排污口布設(shè)在距離入湖口1km、2km、3km、5km、10km、15km、20km位置，排污水平為正常排放和事故排放兩種情況，共有58種計算工況。

4.4 預(yù)測結(jié)果及分析

根據(jù)上述58種計算工況下，王樓礦井入河排污口設(shè)置從不同水期模擬計算結(jié)果分析，排污口設(shè)置后對洙趙新河水質(zhì)沒有明顯的影響。在正常排污條件下，污染物濃度貢獻(xiàn)率BOD5為0.03 mg/L左右，CODcr為0.06 mg/L左右。在上游來水符合功能區(qū)地表水Ⅲ類水質(zhì)目標(biāo)時，在正常排污和事故排污條件下，洙趙新河CODcr和BOD5濃度值均不超過Ⅲ類水標(biāo)準(zhǔn)值，對洙趙新河水功能區(qū)（水域）水質(zhì)、生態(tài)環(huán)境和第三者權(quán)益不會造成影響。

5 結(jié)語

水污染防治是一項復(fù)雜的系統(tǒng)工程，涉及面廣，治理難度大，為確保東線工程水質(zhì)達(dá)到Ⅲ類水標(biāo)準(zhǔn)，必須加強入河排污口的監(jiān)督管理。要通過法律、行政、經(jīng)濟、技術(shù)和工程等綜合措施，實施以節(jié)水為本，清潔生產(chǎn)，治污為先，配套截污導(dǎo)流、污水資源化和流域綜合整治工程，形成“治理、截流、導(dǎo)流、回用、整治”一體化的治污工程體系，實現(xiàn)污水資源化。

MIKE11水動力模塊在入河排污口設(shè)置計算中的應(yīng)用是采用MIKE11一維河道、河網(wǎng)綜合模擬軟件，運用有限差分法對圣維南方程組進(jìn)行數(shù)值求解，對模擬河段進(jìn)行水下地形的數(shù)值概化，對水流、水質(zhì)基本方程進(jìn)行離散求解。解決了南四湖湖西河道由于地形復(fù)雜，長期無法率定參數(shù)建立水流水質(zhì)模型的難題，解決了投資2.635億元、設(shè)計年生產(chǎn)能力45萬t的王樓礦井入河排污口設(shè)置的難題，既滿足了生產(chǎn)急需，又保護(hù)了南四湖的水環(huán)境，同時促進(jìn)了當(dāng)?shù)亟?jīng)濟與社會發(fā)展，為山東省排污口設(shè)置管理提供了科學(xué)的技術(shù)方法。MIKE11水動力模塊為促進(jìn)入河排污口設(shè)置的規(guī)范化管理，優(yōu)化水資源的可持續(xù)利用提供了科學(xué)依據(jù)，經(jīng)濟、社會和環(huán)境效益巨大，具有重大的推廣應(yīng)用價值