水質泥沙一維數學模型在黃河下游應用分析

張海鋒 陳莉紅 李蘭濤

（黃河水利委員會水文局，河南 鄭州 450004）

界定一定水質狀況對應的水體容納污染物的數量為水環境容量，則可用水環境容量表征水環境資源的多少［1］。河流是污染物遷移載體，污染物是借助于河流中的水流和泥沙來輸運的。黃河來流流量的大小影響著河流攜帶污染物的能力，也影響著河流的水環境容量及其自凈能力［2］。進入黃河的泥沙，給黃河造成一定污染，同時由于其組成含有相當數量的黏土礦物、無機膠體和一定數量的有機膠體、有機無機復合膠體，對排入黃河的種類繁多的污染物具有顯著的吸附作用，某種程度上表現出凈化水質的效果。

泥沙對河流水質模擬的影響主要體現在以下幾個方面。①泥沙本身作為具有很大不確定性的面污染源的產物，還具備能夠解吸污染物并可能對水體產生二次污染的特性，大大增加水質模擬的難度。②水體含沙量變化的隨機性和泥沙與污染物相互作用的復雜性，不經使河流水質模型的結構變得復雜，而且水質模型的精度也會相應下降。

本文分析多泥沙河流水流含沙量對河流污染物遷移轉化突出影響，將一維水質數學模型中加入泥沙因素，在簡單河段進行多泥沙河流的氨氮遷移轉化水質模型應用，取得了一些效果。

1 理論基礎

污染物在水體中的遷移轉化與眾多因素有關，是個十分復雜的過程。在諸多影響因素中，河道的河床形態、河流水文、水力條件、河道水體溫度、pH 值、水體初始污染濃度、泥沙吸附與解析、泥沙沉淀與再懸浮是重要的影響條件［3］。

河流內由于各種作用（河段內污染物的降解與增生、沉淀與再懸浮、泥沙的吸附與解吸、引入與引出、降雨與滲漏作用等）使水體的某項污染物在單位時間內的變化量稱為源漏項∑Si，則

式中∑Si為河段水流污染物的源漏項，mg/（L.d）；K1為污染物的一級降解系數，d-1；Ks為含沙量變化對污染物濃度引起的增減速率，d-1；ΔS泥沙變化量，mg/L。

因為河流得離散系數Ed一般要比分子擴散系數Em、紊動擴散系數Et大得多，后者與前者相比，常?？梢院雎圆挥嫞瑒t

于是，可得河流一維水質遷移轉化基本方程形式

式中C 為河段中某污染物濃度，mg/L；t 為時間，d；x為河水得流動距離，km；u為河段水流的平均流速，km/d;E為河段水流的縱向離散系數，km2/d；C1為流入反應單元的水流污染物濃度，mg/L。

對于一個均勻多泥沙河段，可得下面的一維均勻河段的水質遷移轉化基本方程：

河流一維水質遷移轉化基本方程形式：

式中C 為河段中某污染物濃度，mg/L；t 為時間，d；x為河水得流動距離，km；u為河段水流的平均流速，km/d;E為河段水流的縱向離散系數，km2/d;C為流入反應單元的水流污染物濃度，mg/L；

對于一般不受潮汐影響的內陸河流，擴散、離散作用相對于移流作用常常很小，此時可直接解得排污引起的下游濃度變化為：

式中C0為上游端水體污染物濃度，mg/L；S0為上游端水體泥沙濃度，mg/L；S1為上游端水體引起的下游端水體泥沙濃度，mg/L。

2 初步應用驗證

黃河高村站到濼口站河段是黃河典型河段，屬地上懸河，河段基本無入流河道，入河排污口少，河道屬性簡單，為模型驗證提供良好的試驗條件。

黃河下游氨氮降解系數及含沙量變化對污染物濃度引起的增減速率Ks［4］值見表1：

表1 豐枯水期黃河下游K1，KS值

為驗證模型的合理性和可靠性，用2006年山東水文水資源局提供的監測數據進行驗證。氨氮實測值與預測值對比見表2。

表2 氨氮實測值與預測值對比表

由表2可知，采用式4模擬的氨氮濃度值與實測值對比，相對誤差在5%以內，模型所用參數是合理的，模型滿足預測要求，模型可作為水質預測的適用工具。

3 結論

3.1 在多泥沙河流中，泥沙因素是不可忽視的因素，因此本文把源漏項的泥沙項單獨列出來作為一個重點考慮項，十分必要。

3.2 上述數據是在斷面水文條件相似、流量比較穩定、河道均勻、實驗室監測背景相同的條件下取得的，具有一致性，所以計算結果比較好。對于水文條件的河道復雜要考慮離散作用，應對系數k1、ks進行擬合檢驗。

3.3 對于上游來水突然增大或減少，由于泥沙的沉淀與再懸浮作用對系數k1影響很大，在有條件的情況下應對系數ks進行驗證。

3.4 實測資料相對匱乏，在條件允許的情況下，應對系數k1、ks進行校核，探索在復雜條件下泥沙對水質遷移轉化的作用。

3.5 黃河下游山東河道相對簡單，并且污染物僅考慮了氨氮。有關方面應給予支持，進行基于污染物和不同河段經濟研究，探討污染物在復雜黑路條件中的演進規律，泥沙與污染物相互影響的機理，為黃河水資源保護，為保護黃河健康生命做出科學技術支持。

［1］吳澤寧.基于生態經濟的區域水質水量統一優化配置研究［D］.南京：河海大學，2004.

［2］雒文生.水環境分析及預測［M］.武漢：武漢水利電力大學出版社，2000.

［3］白玉川，顧元棪，邢煥政，等.水流泥沙水質數學模型理論及應用［M］.天津：天津大學出版社，2005.

［4］胡國華.多泥沙河流水質模型研究［J］.安全與環境學報，2004.08.