盤錦市水環境質量預警預報二維模型的構建與應用

陳 俊

(遼寧省撫順水文局，遼寧 撫順 113015)

隨著社會經濟的快速，人們逐漸意識到水環境安全面臨越來越嚴重的威脅，開展水環境質量預報模型研究已成為管理水質安全的重要環節之一。通過構建市水環境質量預報模型，能預報水域內各斷面的水質變化趨勢，一旦發生重大水污染事件，該模型對事發地可能出現的水質時空變化過程進行快速模擬，評估水污染事件影響的范圍和級別，啟動應急處理預案，從而為水源地保護管理機構提供快速的決策依據，使管理者根據水域的功能和水體自凈能力提出應對水體中污染物濃度超標的相應措施，以避免出現更大的水污染事故，確保當地居民的飲用水安全，防止出現更嚴重的二次污染，確保水體環境的變化在可控范圍之內，從而保障經濟、社會、環境3種效益的有機統一。近些年來，對于河流水質預報預警國內已取得一定研究成果，但是由于河流水質特征不同，其預警等級很難統一，需要結合區域水質實際狀況，對其水環境預警等級進行標準劃分。尤其是北方獨流入海河流的季節性特點，其汛期和非汛期水體中污染物濃度變化較大，因此對于北方獨立入海河流需要結合其河流水質變化特征制定相對應的水環境質量預警等級，通過構建基于二維非恒定流水質模型，分汛期以及非汛期對水中污染物擴散變化進行模擬。當某一斷面污染物超標時，以超標斷面位置為事故發生原點，利用預報模型模擬污染物濃度隨時間變化過程，確定預警開始時間以及預警解除時間，實現污染物濃度變化趨勢的可視化，為平原地區季節性河流水環境水污染事件的應急處理提供數據支撐。

1 水質預警等級劃分

本文依據GB 3838—2002《地表水環境質量標準》中劃分的5類水，結合盤錦市地表水和地下水質實際狀況，將盤錦市水環境質量預警等級確定為安全狀態、初級預警、中級預警、高級預警、嚴重預警5個等級，分別對應水質類別中的Ⅰ—Ⅴ類水，結果見表1。

表1 盤錦地區水環境預警質量等級

2 水質預報模型建立

2.1 模型基本方程

針對平原地區獨流入海季節性河流特點，采用圣維南方程組構建河流水質擴散模型。嚴格意義上講河流水質模型都是三維結構。但在實際應用中，通常根據實際情況將水體的水質計算簡化為二維、一維甚至零維來計算。本文采用河流水質模型基本方程如下：

(1)

式中，c—污染物濃度，mg/L；t—時間，s；ux、uy、uz—不同方向流速，m/s；Dx、Dy、Dz—不同方向擴散系數，m2/s；K—污染物降解速率，d-1。

2.2 水質模型求解

模型求解過程中需要對計算邊界進行確定，在模型邊界不明或者單一條件下，可以對模型進行解析解的直接計算。不過河流實際上具有復雜的邊界，其水質模型的求解需要采用數學方法進行計算，本文主要對2種邊界條件下的求解方法進行分析。

2.2.1無岸邊限制的點源瞬時模型

近似無岸邊限制的主要針對河面較寬的河流，如果污染源任何一點不是連續穩定排放，其濃度主要為時間函數。在x=0及t=0將M克污染物向河流中投放，則表示為點源的瞬時排放。本文無岸邊限制的點源模型主要用于盤錦市水庫擴散模型計算。水體中污染物的濃度為c(x，y，t)，則點源瞬時排放的水質污染模型計算方程為：

(2)

其計算初始條件為：

(3)

其計算邊界條件為：

(4)

(5)

對上述水質模型進行逆變換，二維瞬時點源擴散方程可變換為：

(6)

當a和b均等于0，則原點即為污染排放點，原點假定坐標為(0，0)，其二維點源擴散模型計算方程為：

(7)

當左岸邊為其污染源釋放點，那么點源瞬時排放的二維擴散模型計算方程為：

(8)

2.2.2有岸邊限制的點源瞬時模型

若河流寬度較低，水面計算范圍較窄，污染總量為M克污染物向河流中投放。污染物擴散到岸邊后，岸邊會猶如鏡子一般對污染物有所反射，另一個像源有所呈現，污染程度有所加劇。再次污染會由于岸邊全發射而發生，本文有岸邊限制的河流瞬時模型對盤錦市主要河流進行二維瞬時點源擴散模型的計算：

exp(-Kt)

(9)

當a=0或a=B時，在岸邊排放點源。當a=B/2，在水面中心進行點源排放。模型中排放口為坐標原點，下游為x方向，對岸為y方向，通常在實際計算中只反射1次即可。

式中，x—排放點距離預測點的距離，m；y—排放口距離預測點的橫向距離，m；t—預測時段，s；c—(x，y)位置點污染物的濃度，mg/L；M—污染物單位時間內的排放量，g/s；a—排放口距離河流岸邊的距離，m；H—平均河流水深，m；Dx和Dy—縱向和橫向河流混合(彌散)；ux和uy—橫向和縱向河流流速均值，m/s；B—河流寬度均值，m；c0—污染物在河流濃度本底值，mg/L；n—反射計算次數；π—圓周率。

河流本底疊加濃度并未在方程中進行考慮，若需要考慮則以污染物只降解不擴散原則進行分析，在方程計算結果基礎上將濃度c0exp(-Kt)進行疊加即可。

2.3 模型參數的確定

對模型參數的合理確定對于模型污染模擬精度影響較大，在水質模型中橫向和縱向擴散系數需要在模型建立時進行確定，本文水質模型參數的確定主要通過經驗公式方法進行估算。

2.3.1橫向擴散系數確定

采用經驗公式對天然河流橫向擴散系數Dy進行確定，方程為：

Dy=aHu*

(10)

(11)

式中，a—量綱為1的橫向擴散系數；g—重力加速度，m/s2；I—水力坡降，%；h—水深均值，m。橫向擴散系數a取值具體見表2。

表2 河流不同類型下的橫向擴散系數a參考值

2.3.2縱向離散系數確定

縱向擴散系數要遠大于縱向離散系數，采用表3中經驗方程對河流縱向離散系數進行確定。本文主要采用中國水電科學院推薦的公式進行縱向離散系數的確定。

表3 河流縱向離散系數經驗計算公式

2.3.3斷面設計流速和水深

河流污染物在水體的傳輸需要河流流量、流速以及水面寬度，對于河流有斷面實測數據的，通過流速、水深和斷面流量之間的正向相關性，對計算河段的水力特征進行計算：

u=aQb

(12)

H=αQβ

(13)

式中，Q—河流斷面流量，m3/s；u—河流斷面平均流速，m/s；H—斷面平均水深，m。a、b、α、β—經驗系數，其取值主要通過河段水文數據進行確定。

對方程(12)進行線性轉換，轉化后的方程為：

y=B0+B1x

(14)

因此，y和x的線性回歸方程即可表達為斷面流速均值，和2020年的汛期和非汛期實測水文資料進行相關性分析，得到盤錦地區主要河流流速和水深流速之間的相關方程，具體見表4。

表4 盤錦市河流斷面經驗公式

對于無實測水文資料的河段，河道比降以及糙率可以采用曼寧公式進行平均流速的計算，或者采用內插、類比方法對設計流速和水深進行推求。

3 模型應用

3.1 模型驗證

本文以盤錦地區杜家水文站為驗證計算節點，結合杜家水文站汛期和非汛期2015—2020年實測苯酚采樣監測指標進行模型驗證，選擇在有岸邊限制條件下的瞬時點源二維水質擴散模型進行計算，由方程(9)研究確定苯酚降解系數K=0，河流中心為點源釋放點，污染物揮發酚標準濃度為0.005，污染物濃度C0的本底值0.0029mg/L，汛期和非汛期模型參數見表5。

表5 模型參數設置

表6 模型驗證結果

以苯酚作為模擬污染指標，苯酚主要用于農藥以及化工企業，具有很強的毒性且很難在水體中得到有效降解。從2015—2020年杜家水文站的斷面苯酚監測指標可看出，其各年份苯酚指標汛期和非汛期濃度均低于0.004mg/L，汛期由于進入河流內的污染物有所增多其汛期苯酚濃度總體要高于非汛期。通過建立的點源二維水質擴散模型驗證結果可看出，各年份下模型模擬的苯酚濃度和監測濃度之間的誤差總體低于±30%，滿足水質模擬的規范精度要求，可以用建立的水質模擬進行不同突發水質情景模擬分析。

3.2 突發情景模擬分析

假定向紅旗水庫上游的杜家水文站斷面附近釋放1t的揮發酚污染源，通過建立的瞬時點源二維水質模型，選擇有岸邊限制的邊界條件，結合MATLAB軟件對紅旗水庫上游污染物濃度變化進行模擬，汛期和非汛期不同時刻的污染物濃度擴散模擬結果如圖1—2所示。

圖1 汛期不同時刻揮發酚濃度空間分布

從汛期和非汛期盤錦地區紅旗水庫上游不同時間的揮發酚濃度維空間變化圖可看出，污染團隨著水流向下游漂移，污染物濃度逐步降低但污染范圍有所增加。非汛期污染物濃度降低幅度明顯低于汛期污染物濃度降低幅度。因此，水流流速是污染物移動速率的主要因素，污染物向下游漂移因素的決定性因素是順水方向的河流流速。橫向和縱向擴散系數是污染范圍增加的主要影響因子。水體的自凈能力以及污染物的自身降解是其濃度逐步降低的主因。通過模擬可看出水體中污染物濃度隨著時間的增加逐步減小。當污染物擴散范圍達到下游臨界值后，最高濃度會低于標準值。當擴散范圍低于下游臨界值時，水體污染物濃度會在某時段內超過標準濃度值。對于水質控制點，當該污染物濃度超過標準濃度，對應表1中的預警等級濃度標準，則認為發生突發水質事件，此時進行水質預警。紅旗水庫位于杜家水文站下游2km處，設有應急入盤取水口，本文對不同時段應急取水口的揮發物濃度也進行了模擬，通過模擬分析，當揮發酚濃度超標8.3min后，即發出突發水質預警，在25.3min后濃度回降到標準濃度以下，從發出預警到預警解除時間為17min。非汛期預警開始時間為66.5min，解除預警的時間為106.5min，污染預警持續時間為40min。相比于汛期，非汛期峰值濃度預警時間有所延緩。

圖2 非汛期不同時刻揮發酚濃度空間分布

4 主要結論

(1)在構建水質預報模型時，若需要疊加河流本底濃度時，則以污染物只降解不擴散原則進行分析，在二維點源擴散模型計算方程基礎上將濃度c0exp(-Kt)進行疊加即可。

(2)對于無實測水文資料的河段的水質模擬，可以采用曼寧公式結合河道比降及糙率進行平均流速的計算，或者采用內插、類比方法對設計流速和水深進行推求。

(3)構建好的水質模型，可以為城市水質突發應急事件作為模擬平臺，當預測污染物濃度超過預警等級對應的標準值后，即發出預警，當濃度低于預警等級標準后，即可解除預警，預警時間內污染物濃度時空變化可以為水質突發應急響應措施提供重要的支撐。