涇河甘陜緩沖區(qū)納污能力淺析

王海強(qiáng)，馬 康

（咸陽水文水資源勘測中心，陜西 咸陽712000）

1 概況

涇河甘陜緩沖區(qū)內(nèi)有豐富的煤礦資源，隨著煤炭資源開發(fā)利用，區(qū)域內(nèi)用水需求不斷增加，隨之排水加大了甘陜緩沖區(qū)水資源承載力，同時涇河省界的污染主要是涇河上游入境及馬蓮河入境的影響，造成對涇河甘陜緩沖區(qū)水資源水質(zhì)變化加劇，在保護(hù)涇河水環(huán)境前提下，高質(zhì)量發(fā)展，符合現(xiàn)有的發(fā)展趨勢，隨著水污染防治工作的不斷推進(jìn)，涇河水質(zhì)的不斷改善，根據(jù)目前經(jīng)濟(jì)和水生態(tài)的發(fā)展布局需要，明確水功能區(qū)達(dá)標(biāo)情況，結(jié)合景村水文站的長系列水文資料，綜合考慮河流支流匯入、排污口設(shè)置等因素，對涇河甘陜緩沖區(qū)納污能力及限排總量進(jìn)行了水文分析計算，并核定該甘陜緩沖區(qū)納污能力及限排總量。

2 納污水域概況

涇河是黃河的二級支流，是渭河的最大支流，發(fā)源于寧夏回族自治區(qū)涇源縣六盤山東麓，于陜西省高陵縣境內(nèi)注入渭河，全長455 km，流域面積45421 km2，是彬長礦區(qū)的最大過境河流。涇河干、支流河道呈羽狀展布，從整個流域而言，絕大部分支流分布于干流中、上游地區(qū)。涇河流域橫跨黃土塬區(qū)，整個流域大致呈扇形分布，總體地形西北高、東南低。

涇河流域的水功能區(qū)劃涉及寧夏回族自治區(qū)、甘肅省和陜西省。根據(jù)黃河流域水功能區(qū)劃成果，涇河流域一級水功能區(qū)有32個，涉及涇河干流及11條支流；二級水功能區(qū)22個，涉及涇河干流及8條支流。

按照《全國重要江河湖泊水功能區(qū)劃（2011-2030）》《陜西省水功能區(qū)劃》（2004）涇河水功能區(qū)在咸陽劃分為兩段，分為水功能一級區(qū)和水功能二級區(qū)。涇河水功能區(qū)見圖1。

圖1 涇河水功能區(qū)劃圖

3 涇河甘陜緩沖納污能力計算

涇河甘陜緩沖區(qū)納污能力計算步驟：①分析設(shè)計流量，確定區(qū)段設(shè)計流量Q和流速U；②確定來水濃度C0、出水濃度Cs。③分析模型計算參數(shù)k值；④分析計算區(qū)段現(xiàn)狀，選取計算模型并計算納污能力。

根據(jù)《水域納污能力計算規(guī)程》規(guī)定：河流流量基本滿足Q≤150 m3/s，縱向擴(kuò)散項遠(yuǎn)小于平流輸送項，從而忽略縱向擴(kuò)散項，污染源是連續(xù)點源，可僅考慮污染物的降解；排污口入河污染物在橫斷面上均勻混合。涇河甘陜段基本滿足以上幾個條件，因此采用一維水質(zhì)模型進(jìn)行計算。對于功能區(qū)下斷面，其污染物濃度計算公式為：

相應(yīng)的水域納污能力計算公式為：

式中：W為納污能力，g/s；CX為水體質(zhì)量標(biāo)準(zhǔn)，mg/L；C0為污染物背景濃度，mg/L；X為水域計算距離，m；Q為水域流量，m3/s；Qp為排污口排放流量，m3/s；U為水域平均流速，m/s；K為污染物衰減系數(shù)，1/s。

3.1 設(shè)計水文條件

（1）在計算水功能區(qū)下斷面的納污能力時，下斷面流量是計算關(guān)鍵，設(shè)計流量的大小對納污能力的計算結(jié)果影響很大，流量資料系列太短則無法反映水文規(guī)律，資料太長則無法反映人類活動對水資源造成的影響，特別是對枯水期小流量的影響。因此，在計算胡家河村斷面流量時，采用涇河景村水文站實測月平均流量推算胡家河村斷面處的月平均流量。

依據(jù)《水域納污能力計算規(guī)程》（GB/T 25173-2010）要求，有水文長系列資料時，現(xiàn)狀設(shè)計流量的確定選用設(shè)計保證率的最枯月平均流量，采用頻率計算法計算。無水文長系列資料時，可采用近10年系列資料中的最枯月平均流量作為設(shè)計流量。無水文資料時，可采用內(nèi)插法、水量平衡法、類比法等方法推求設(shè)計流量。

胡家河村斷面距景村水文站43.2 km，選取景村水文站作為參證站?，F(xiàn)狀條件下，一般采用最近10年最枯月平均流量或90%保證率最枯月平均流量作為設(shè)計流量。依據(jù)景村水文站近10年枯水月平均流量值，采用比擬法計算胡家河村斷面處的枯水月平均流量。設(shè)計流量計算是基于景村水文站的實測流量資料，對歷年月流量進(jìn)行統(tǒng)計，將其從大到小排列，經(jīng)驗頻率計算公式：

式中：p為設(shè)計保證率；m為樣本秩數(shù)；n為樣本總數(shù)。

將頻率計算結(jié)果點繪成對數(shù)頻率曲線，用P-Ⅲ型曲線與之適配，據(jù)p-III型曲線查90%保證率所得流量即為設(shè)計流量，表1為景村站近11年每月平均流量統(tǒng)計表。

表1 景村站2009年～2019年月平均流量統(tǒng)計表 單位：m3/s

按照10年最枯月平均流量：選用水文年1月～3月、11月～12月為枯水期，計算Q月均為13.7 m3/s。按照90%保證率最枯月平均流量：選用景村站20年～30年最枯月流量進(jìn)行統(tǒng)計，用P-Ⅲ型曲線與之適配，見圖2景村站最枯流量P-Ⅲ型曲線圖。

圖2 景村站最枯流量P-Ⅲ型曲線圖

表2 景村站最枯流量不同頻率流量成果表

通過兩種方法計算，11年最枯月平均流量13.7 m3/s，90%保證率最枯月平均流量13.6 m3/s，兩種方法計算出的流量基本一致。為了反映該段（涇河甘陜緩沖區(qū)）水文年際周期變化和其中長期發(fā)展趨勢，流量資料系列視資料情況取10年枯水期月平均流量值13.7 m3/s。

以景村站斷面流量為依據(jù)，采用水文比擬法，確定涇河彬州市胡家河村斷面流量，見表3。

表3 胡家河村斷面流量成果表

2）斷面設(shè)計流速確定

按下式計算：

式中：V為設(shè)計流速；Q為設(shè)計流量；A為過水?dāng)嗝婷娣e。

對涇河胡家河村河道斷面實測，利用曼寧公式推求該斷面水位～流量和水位～面積關(guān)系，曼寧公式：

式中：Q為流量；A為斷面面積；n為糙率；R為水力半徑；i為水面比降。

依據(jù)曼寧公式推求該斷面水位～流量和水位～面積關(guān)系，計算出胡家河村斷面設(shè)計流速V，見表4。

圖3 胡家河村斷面水位-流量關(guān)系曲線

圖4 胡家河村斷面水位-面積關(guān)系曲線

表4 胡家河村斷面設(shè)計流速成果表

通過計算涇河甘陜緩沖區(qū)胡家河村斷面10年枯水期月平均流量值為11.8 m3/s，流速為0.674 m/s（采用景村水文站資料）。

3.2 控制參數(shù)

（1）控制因子及水質(zhì)標(biāo)準(zhǔn)

參照最嚴(yán)格的水資源管理制度納污控制考核的要求，結(jié)合我省多年來在污染控制方面的實際，確定COD、NH3-N為本次納污能力計算的因子。

表5 COD、NH3-N水質(zhì)標(biāo)準(zhǔn)

（2）綜合衰減系數(shù)K的確定

為簡化計算，在水質(zhì)模型中，將污染物在水環(huán)境中的物理降解、化學(xué)降解和生物降解概化為綜合衰減系數(shù)，確定合理的污染物綜合衰減系數(shù)是模擬污染物在水體中濃度衰減的關(guān)鍵之一。分別采用以下方法計算，并進(jìn)行綜合確定：

1）實測法

為較準(zhǔn)確地掌握河流的綜合衰減系數(shù)K，在景村水文站設(shè)置了監(jiān)測斷面，同步監(jiān)測水文及水質(zhì)等各項要素，并與長系列資料比較，經(jīng)分析，確定了較可靠的干、支流、排污口的水文、水質(zhì)監(jiān)測資料，利用實測資料反推法確定綜合袞減系數(shù)K值，對取得數(shù)據(jù)進(jìn)行合理性和可靠性分析后采用計算公式如下：

式中：K為污染坊綜合降解系數(shù)，1/d；CA為河段上斷面污染物濃度，mg/L；CB為河段下斷面污染物濃度，mg/L；u為河段平均流速，m/s；x為上下斷面間的距離，km。

計算得到?jīng)芎痈闪鞯腃OD的k值為0.18d-1～0.24d-1，NH3-N的k值為0.20d-1～0.30d-1。

2）分析借用法

本次納污能力計算，綜合衰減系數(shù)（K）的確定依據(jù)2015年通過水利部審查的《黃河流域（片）重要江河湖泊水功能區(qū)納污能力核定和分階段限制排污總量控制方案》。分析確定涇河甘陜緩沖區(qū)干流納污能力計算中綜合衰減解系數(shù)（K）的取值為：化學(xué)需氧量的K值為0.20d-1，氨氮的K值為0.28d-1。涇河化學(xué)需氧量的K值為0.21d-1，氨氮的K值為0.29d-1?？紤]涇河上下游水力特性、污染狀況等條件的差異，此值可作為參考值使用。

背景濃度Co和控制斷面濃度CX是計算納污能力的另外兩個參數(shù)，它們反映了水功能區(qū)來水和出水質(zhì)量，其取值主要考慮功能區(qū)的水質(zhì)目標(biāo)。

根據(jù)2015年～2019年對胡家河村斷面以上6.1 km處的政平斷面每月水質(zhì)監(jiān)測，監(jiān)測的結(jié)果，在政平斷面處的COD的監(jiān)測值年平均值都比目標(biāo)值大。NH3-N的監(jiān)測值年平均值都比目標(biāo)值小。按照三條紅線考核的相關(guān)規(guī)定，跨界水功能區(qū)超標(biāo)原因是受上游來水影響的可以對污染物量進(jìn)行核減的原則，本次納污能力計算C0值采用水功能區(qū)水質(zhì)目標(biāo)的上限，分別為COD 20 mg/L，氨氮1.0 mg/L，CS采用論證水功能區(qū)的水質(zhì)目標(biāo)Ⅲ類水的限制，分別為20.0 mg/L及1.0 mg/L。

通過計算涇河甘陜緩沖區(qū)納污能力分別為：COD 1254.7 t/a，NH3-N 89.3 t/a，主要參數(shù)見表6。

表6 涇河甘陜緩沖區(qū)納污能力計算表

4 涇河甘陜緩沖限制排放總量計算

限制排污總量控制計算方法，水功能區(qū)限制排污總量在核定水域納污能力的基礎(chǔ)上，結(jié)合《全國水資源綜合規(guī)劃》和黃河流域綜合規(guī)劃修編成果、區(qū)域經(jīng)濟(jì)技術(shù)水平、河流水資源配置等因素，嚴(yán)格控制入河排污總量，綜合確定重要江河湖泊水功能區(qū)分階段限制排污總量方案。

以2019年為例，分為有、無污染物入河量資料兩種情況，涇河甘陜緩沖區(qū)有污染物排放。依據(jù)《全國水資源綜合規(guī)劃》黃河干流、主要支流及西北主要河流的保護(hù)區(qū)、省界緩沖區(qū)、飲用水源區(qū)及其他重要水功能區(qū)，原則上應(yīng)在2015年達(dá)到水功能區(qū)水質(zhì)目標(biāo)要求，根據(jù)污染布局，以核定的納污能力或小于納污能力綜合確定2015年限制排污總量。

預(yù)計2020年、2030年仍不能實現(xiàn)水功能區(qū)水質(zhì)達(dá)標(biāo)的，按照從嚴(yán)控制、未來有所改善的要求，確定水功能區(qū)2020年、2030年限制排污總量。具體可以考慮如下幾種控制方案：

1）根據(jù)水功能區(qū)污染程度，考慮社會經(jīng)濟(jì)發(fā)展水平、污染治理水平及其可達(dá)性，按一定的入河消減百分比提出階段污染物限排總量。

2）考慮地區(qū)水資源條件、水功能區(qū)現(xiàn)狀水質(zhì)、現(xiàn)狀污染物入河排放量以及污染治理水平等因素，采取75%、50%或者其他保證率設(shè)計條件計算的納污能力作為階段污染物限排總量。

2020年、2030年的水功能區(qū)限制排污總量方案參照以上方法。

依據(jù)2019年對政平斷面的監(jiān)測的結(jié)果（見表7，分析功能區(qū)分階段限制排污總量，采用2019年COD及NH3-N的監(jiān)測值年枯水月最最小值為作為C0值。即，分別為COD 16.4 mg/L，氨氮0.41 mg/L為作為C0值。

表7 2019年政平斷面COD及氨氮實際監(jiān)測值表 單位：mg/L

方案1：采用入河消減的75%提出階段污染物限排總量

依據(jù)水功能區(qū)污染程度，采用政平斷面2019年的平均值COD及氨氮的數(shù)據(jù)，計算涇河甘陜緩沖區(qū)水環(huán)境容量，計算公式：

式中：W為水環(huán)境容量，t/a；K為污染物綜合降解系數(shù)，1/d；C0為河段上斷面污染物濃度，mg/L；Ey為橫向彌散系數(shù)，m2/s；Cs為功能區(qū)水質(zhì)標(biāo)準(zhǔn)濃度，mg/L；u為河段平均流速，m/s；L為河道長度，km；h為河段平均水深，m；86.4為單位換算系數(shù)。

通過計算涇河甘陜緩沖區(qū)水環(huán)境容量分別為：COD 2474 t/a，NH3-N 312 t/a，主要參數(shù)見表8。

通過涇河甘陜緩沖區(qū)水環(huán)境容量計算，得到該段緩沖區(qū)入河消減量見表9。

表8 涇河甘陜緩沖區(qū)水環(huán)境容量計算表

表9 涇河甘陜緩沖區(qū)入河消減量成果表

按照入河削減的75%提出階段污染物限排總量，分別COD 206.3 t/a，NH3-N 26.3 t/a。

方案2：取75%保證率設(shè)計條件計算的納污能力作為階段污染物限排總量。

涇河甘陜緩沖區(qū)胡家河村斷面75%保證率下的設(shè)計流量為13.7 m3/s，流速為0.698 m/s（采用景村水文站資料）。

在政平斷面處的COD的監(jiān)測值年平均值都22.85 mg/L，氨氮0.871 mg/L為作為C0值。

通過計算涇河甘陜緩沖區(qū)限制排放分別為：COD 167.7 t/a，NH3-N 89.3 t/a，主要參數(shù)見表10。

表10 涇河甘陜緩沖區(qū)限制排放計算表

通過方案1和方案2比較，方案1考慮到水環(huán)境容量，考慮因素較多，采用方案1比較合理。依據(jù)2019年限制排放計算，2020年、2030年的水功能區(qū)限制排污總量見表11。

表11 涇河甘陜緩沖區(qū)限制排放成果表

5 結(jié)論

本文通過計算涇河甘陜緩沖區(qū)的納污能力COD、氨氮入河量分別為1254.7 t/a、89.3 t/a；確定2020年～2030年的水功能區(qū)限制排污總量染物COD、氨氮入河量分別為206.3 t/a、26.3 t/a。