劉軍政,猶 偉,祁蕓泉,趙鐘楠,黃火建,張新華

(1.四川大學(xué)水力學(xué)及山區(qū)河流開發(fā)保護(hù)國家重點實驗室，四川 成都 610065；2.四川大學(xué)水利水電學(xué)院，四川 成都 610065； 3.水利部水利水電規(guī)劃設(shè)計總院，北京 100011)

當(dāng)前，新實施的《中華人民共和國環(huán)境保護(hù)法》(自2015年1月1 日起施行)和中國共產(chǎn)黨第十九次全國代表大會均對實施生態(tài)補償提出了明確的要求，生態(tài)補償制度已成為推進(jìn)科學(xué)發(fā)展、加快生態(tài)文明制度建設(shè)的重要舉措和生態(tài)環(huán)境保護(hù)管理體制改革創(chuàng)新的重要內(nèi)容。我國從20世紀(jì)90年代末開始關(guān)注流域生態(tài)補償，如今已經(jīng)形成了多種補償類型和補償方式并存的局面[1]。目前，江蘇、福建、浙江、河北、河南、遼寧、山東、山西等省在全省范圍內(nèi)開展了污染賠償型生態(tài)補償試點工作[2-3]。其中，生態(tài)補償標(biāo)準(zhǔn)是當(dāng)前研究的熱點與難點，眾多學(xué)者對此進(jìn)行了研究，如劉桂環(huán)等[4]、冉濤等[5]總結(jié)我國生態(tài)補償?shù)膶嵺`，提出了2種對于污染賠償型的生態(tài)補償標(biāo)準(zhǔn)核算方法，一種是按照水質(zhì)濃度超標(biāo)倍數(shù)扣繳一定數(shù)額的補償金；另一種是按照考核斷面超標(biāo)污染物通量的削減成本確定補償金額。前者對生態(tài)補償?shù)暮怂愫啽?，但誤差大。已有的研究和實踐多數(shù)主張采用污染物通量核算補償量，污染物通量全面考慮了水質(zhì)與總量的關(guān)系，具有較強的科學(xué)性[6-8]。然而，目前該方法中單位污染物削減成本統(tǒng)一采用某一固定值，沒有考慮工業(yè)廢水、城市生活用水等不同廢水處理成本的差異，與實際值可能存在偏差。同時，在污染物通量核算上忽略了年內(nèi)各水期的實際水文、水質(zhì)特點，故也存在一定偏差。因此，筆者在基于Webb等[9]和郝晨林等[10]研究污染物通量核算模型的基礎(chǔ)上，以瀨溪河渝川緩沖區(qū)為研究對象，從污染物削減成本和污染物通量兩個方面進(jìn)行改進(jìn)研究。

1 研究方法

1.1 基于污染物通量超標(biāo)排放量的生態(tài)補償標(biāo)準(zhǔn)核算方法

該方法首先求得時段內(nèi)超過目標(biāo)水質(zhì)的排放總量，再與單位排放量補償金額進(jìn)行乘積求得補償總金額，該方法從宏觀上把握了河流上下游斷面污染物的遷移規(guī)律。補償資金計算公式為

(1)

式中：P單為單因子補償資金；W超為考核斷面實測水質(zhì)超過目標(biāo)水質(zhì)的污染物排放量，若實測水質(zhì)濃度低于目標(biāo)水質(zhì)濃度則為負(fù)，即此時下游向上游反向補償；C(t)實測為t時刻考核斷面實測水質(zhì)濃度；C(t)目標(biāo)為t時刻考核斷面目標(biāo)水質(zhì)濃度；Q(t)為t時刻通過考核斷面流量；S為單位超標(biāo)污染物排放通量的補償金額，可視為單位污染物削減的處理成本。

考慮到實際工作中很難取得河流每個瞬時的流量與污染物濃度值，且我國當(dāng)前多數(shù)水文水質(zhì)監(jiān)測站監(jiān)測數(shù)據(jù)不匹配[11]，多數(shù)水文站可達(dá)到流量每日一測，而大多數(shù)水質(zhì)監(jiān)測站則只能做到每月一測。因水文、水質(zhì)數(shù)據(jù)時間跨度大、不連續(xù)，所以需對污染物通量估算公式進(jìn)行離散。

(2)

式中：C(i)實測為考核斷面計算時期內(nèi)第i次水質(zhì)實測瞬時濃度值；C(i)目標(biāo)為考核斷面計算時期內(nèi)第i次水質(zhì)目標(biāo)濃度值；Q(i)為第i次監(jiān)測代表時間段的流量；n為代表監(jiān)測時段內(nèi)污染物監(jiān)測次數(shù)。

若需考核多因子指標(biāo)，則最終補償金額為將各單因子補償金額進(jìn)行正負(fù)疊加即可。

1.2 污染物通量估算方法

時段污染物通量是污染賠償型生態(tài)補償?shù)闹匾罁?jù)，它的計算可簡化為以下5 種方法[9]。

(3)

(4)

(5)

(6)

(7)

由于實際工作中很難取得河流每個瞬時的流量與污染物濃度值，數(shù)據(jù)相對比較稀缺。因此，如何根據(jù)離散的水量、水質(zhì)數(shù)據(jù)選擇恰當(dāng)?shù)耐抗浪惴椒?，對提高通量估算的可靠性具有十分重要的意義。筆者在富國[12]和郝晨林等[10]對5種通量估算公式適用性研究的基礎(chǔ)上，對研究區(qū)域劃分不同水期，根據(jù)不同的水文、水質(zhì)條件，判定不同水期內(nèi)的污染物來源，從而選擇水期內(nèi)相應(yīng)的通量估算公式，進(jìn)行通量估算，以此作為計算補償金額的依據(jù)(表1)。

表1 5種通量估算方法特征與適用性

1.3 單位污染物超標(biāo)排放量補償金額的確定

本文以單位污染物削減成本計算確定單項污染物的補償金額。王佳偉等[13]、王曉青[14]、李燁楠等[15]的研究表明：不同來源的廢水污染物削減成本不同，工業(yè)廢水顯著高于城市生活污水削減成本。而目前污染物削減成本統(tǒng)一采用定值，沒有考慮工業(yè)廢水、城市生活污水等不同類型廢水處理成本的差異性。因此，筆者將污染物削減成本按照污染物來源分為點源污水(主要為工業(yè)廢水、城市生活污水)和農(nóng)村非點源污水(主要為農(nóng)村生活污水、農(nóng)業(yè)生產(chǎn)污水)削減成本兩部分。根據(jù)不同廢水來源的處理成本，按污染物排放量比重加權(quán)平均得到削減成本。

2 實例研究與結(jié)果分析

2.1 研究區(qū)域概況

瀨溪河為沱江左岸一級支流，發(fā)源于重慶市大足縣，流經(jīng)榮昌縣和四川省瀘縣，最終在瀘縣注入沱江，屬于長江流域沱江水系。河流全長約190 km，平均坡降約0.11%，其中重慶市境內(nèi)河長123 km，四川省境內(nèi)長65.5 km。流域總面積3 257 km2，多年平均徑流量11.67億m3，多年平均流量37 m3/s，枯水期多年平均流量4.14 m3/s。瀨溪河渝川緩沖區(qū)起始斷面為重慶市榮昌縣，終止斷面為四川省瀘縣(位于川渝省界交界處)，總長16.5 km。瀘縣福集鎮(zhèn)有一處水文站，緩沖區(qū)出境斷面設(shè)有水質(zhì)監(jiān)測斷面。瀨溪河徑流年內(nèi)分配不均勻，汛期較長，多出現(xiàn)在6～10月，徑流量占年徑流總量的70%左右；平水期較短，多為4～5月和11月，徑流量占年徑流總量的20%左右；其余月份為枯水期，徑流量很小，只占年徑流總量的10%左右。瀨溪河面臨工業(yè)點源和農(nóng)業(yè)非點源污染較重的情況，水質(zhì)狀況較差，常年達(dá)不到Ⅲ類水質(zhì)標(biāo)準(zhǔn)。流域主要超標(biāo)水質(zhì)項目為高錳酸鹽指數(shù)、氨氮和CODcr，屬于有機(jī)污染型河流。污染來源除受上游來水影響外，還受流域內(nèi)密集工礦區(qū)廢水及農(nóng)村面源污染影響，其中點源污染來源主要為工業(yè)廢水和城市生活污水，非點源污染來源主要為農(nóng)村生活污水、農(nóng)田徑流污染、畜禽養(yǎng)殖污染等。2015年調(diào)查數(shù)據(jù)顯示，瀨溪河流域瀘縣境內(nèi)廢水排放3 061萬t，高錳酸鹽指數(shù)排放量4 870 t，氨氮排放量1 244 t。

2.2 污染物通量估算及其改進(jìn)

2.2.1 污染物通量估算及存在的問題

四川省人民政府辦公廳關(guān)于《在岷江沱江流域試行跨界斷面水質(zhì)超標(biāo)資金扣繳制度的通知》(川辦函[2011]200號)[16]中規(guī)定：岷江、沱江流域，水質(zhì)補償項目選擇高錳酸鹽指數(shù)和氨氮兩項指標(biāo)。故筆者采用這兩項污染物指標(biāo)進(jìn)行上、下游生態(tài)補償金額的計算，所以只核算這兩種污染物指標(biāo)的通量。

瀨溪河渝川緩沖區(qū)重慶出境斷面處有一福集水文監(jiān)測站，并有多年逐日流量監(jiān)測數(shù)據(jù)，緩沖區(qū)下游四川入境斷面設(shè)置有水質(zhì)監(jiān)測斷面，監(jiān)測頻率為每月一次。筆者從多年流量系列中選取2015年監(jiān)測的逐日流量數(shù)據(jù)和每月一次的水質(zhì)數(shù)據(jù)，采用式(3)～(7)分別對渝川緩沖區(qū)進(jìn)行污染物通量估算，具體結(jié)果見表2。

表2 5種污染物通量估算方法核算結(jié)果

由基礎(chǔ)資料可知，該省界緩沖區(qū)僅在下游斷面布設(shè)有水文、水質(zhì)站，下游斷面要求水質(zhì)達(dá)到Ⅲ類標(biāo)準(zhǔn)，根據(jù)我國《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)中的規(guī)定可知，當(dāng)水質(zhì)類別為Ⅲ類時，高錳酸鹽指數(shù)和氨氮最大限值分別為6 mg/L和1 mg/L，2015年年均流量為35.5 m3/s，相應(yīng)的最大容許污染物的年通量如表3所示。

將5種估算方法核算的結(jié)果與下游斷面容許限值進(jìn)行對比(表4)，各污染物超標(biāo)數(shù)據(jù)可作為上下游生態(tài)補償賠償金額的依據(jù)。

表3 下游斷面水質(zhì)要求及容納負(fù)荷通量

表4 5種估算方法污染物超標(biāo)情況

注：“正值”代表上游污染物超標(biāo)，“負(fù)值”代表上游污染物未超標(biāo)，下同。

表4為采用不同污染物通量估算方法核算各污染物超標(biāo)排放量，通過對比可知：式(3)、式(4)、式(6)、式(7)高錳酸鹽指數(shù)核算結(jié)果均不超標(biāo)(即表4中的W1、W2、W4、W5對應(yīng)值)；而式(5)中高錳酸鹽指數(shù)核算結(jié)果超標(biāo)471 t(即表4中W3對應(yīng)值)。不同方法對于同一區(qū)域污染物超標(biāo)核算結(jié)果大不相同，甚至出現(xiàn)矛盾。對于氨氮而言，雖然5種估算方法核算污染物是否超標(biāo)未出現(xiàn)矛盾，但是超標(biāo)總量核算結(jié)果差別較大。所以不同通量估算方法在補償中會出現(xiàn)不同結(jié)果，甚至出現(xiàn)矛盾。

2.2.2 污染物通量估算改進(jìn)

將該流域全年劃分為豐、平、枯3個水期，考慮不同水期實際的水文、水質(zhì)特征，選擇相應(yīng)的通量估算公式。最終得到高錳酸鹽指數(shù)和氨氮的超標(biāo)排放量，具體結(jié)果見表5。

表5 分水期改進(jìn)后污染物超標(biāo)情況

方法改進(jìn)后，將全年劃分為豐水期、平水期、枯水期，根據(jù)不同水期污染物來源不同，采用相應(yīng)的污染物通量估算公式，核算的高錳酸鹽指數(shù)、氨氮年通量分別為4 957 t、1 253 t，與下游斷面2種污染物年通量容許限值6 735 t、1 122 t對比可知，高錳酸鹽指數(shù)未超標(biāo)，氨氮超標(biāo)131 t。根據(jù)2015年調(diào)查資料顯示：高錳酸鹽指數(shù)、氨氮年排放量分別為4 870 t、1 144 t，可以看出污染物超標(biāo)排放量核算結(jié)果比較接近實際調(diào)查情況，因此核算的補償金額結(jié)果也會更加真實可信。

2.3 單位污染物超標(biāo)排放量補償金額的確定及其改進(jìn)

四川省對岷江、沱江流域的主要污染物處理成本規(guī)定為：高錳酸鹽指數(shù)1 300元/t，氨氮17 000元/t，污染物削減成本統(tǒng)一采用定值，沒有考慮不同類型廢水處理成本的差異性。因此，筆者將污染物削減成本按照污染物來源分為點源污水(主要為工業(yè)廢水、城市生活污水)和農(nóng)村非點源污水(主要為農(nóng)村生活污水、農(nóng)業(yè)生產(chǎn)污水)削減成本兩部分。根據(jù)不同廢水來源的處理成本，按污染物排放量比重加權(quán)平均得到削減成本。

工業(yè)廢水處理成本按照重慶市典型企業(yè)污染物削減成本的加權(quán)均值分別為：高錳酸鹽指數(shù) 1 600元/t、氨氮18 000元/t；城市生活污水2種污染物削減成本分別為160元/t和9 150元/t；農(nóng)村生活污水2種污染物削減成本分別為190元/t和10 000元/t；農(nóng)業(yè)生產(chǎn)污水分別為1 300元/t、15 000元/t。不同污染源單位削減成本及排放比例見表6

表6 不同污染源削減成本

按照污染物排放量比例加權(quán)平均最終得到點源和非點源的單位削減成本。其中，點源中高錳酸鹽指數(shù)總成本為1 283元/t，氨氮總成本14 903元/t;非點源中高錳酸鹽指數(shù)總成本為490元/t，氨氮總成本12 350元/t。由計算結(jié)果可知：不同廢水來源的污染物削減成本不同，且對于高錳酸鹽指數(shù)與氨氮而言，點源污染比非點源污染削減成本高。

2.4 補償金額計算結(jié)果

2.4.1 統(tǒng)一削減成本補償金額計算結(jié)果

方法改進(jìn)后估算的污染物通量比較接近實際調(diào)查量，為了比較削減成本改進(jìn)前后的效果，采用污染物通量估算方法改進(jìn)后的污染物超標(biāo)量核算補償金額，具體補償金額計算結(jié)果見表7。

表7 統(tǒng)一削減成本補償金額計算結(jié)果

由表7可看出，采用統(tǒng)一削減成本核算補償金額，該省界緩沖區(qū)上游斷面重慶市因超標(biāo)排放氨氮需賠償下游斷面四川省223萬元，少排放高錳酸鹽指數(shù)又可獲得補償231萬元，所以疊加后重慶市獲得補償8萬元。

2.4.2 分水期改進(jìn)后補償金額計算結(jié)果

將該區(qū)域全年劃分為豐水期、枯水期、平水期，分別計算不同水期內(nèi)的污染物超標(biāo)量和污染物削減成本，通過不同水期的補償金額進(jìn)行疊加，得到全年的補償金額，具體污染物補償結(jié)果見表8。

表8 分水期改進(jìn)后補償金額計算結(jié)果

由表8可看出，該省界緩沖區(qū)上游斷面重慶市因少排放高錳酸鹽指數(shù)可獲得補償124萬元，超標(biāo)排放氨氮又需賠償下游斷面四川省205萬元，所以疊加后需賠償四川省81萬元。而削減成本改進(jìn)前四川省需賠償重慶市8萬元，顯然不符合上游污染下游的實際情況。兩者補償金額結(jié)果出現(xiàn)差異性的主要原因是豐水期、平水期高錳酸鹽指數(shù)主要污染來源為非點源污染，單位污染物削減成本較低，而全年采用統(tǒng)一削減成本，會導(dǎo)致補償金額計算較大。所以根據(jù)不同水期內(nèi)實際水質(zhì)特征，判定不同水期污染物來源，從而確定不同污染源的削減成本，不采用統(tǒng)一固定值，這樣核算補償金額更加合理。

3 結(jié) 論

通過分析研究污染物通量超標(biāo)排放量核算補償量方法，根據(jù)年內(nèi)水文、水質(zhì)變化特征對污染物通量的估算和削減成本的確定進(jìn)行了以下改進(jìn)：

a. 根據(jù)流域的實際水文、水質(zhì)特點劃分豐、平、枯水期，判定不同水期的污染物來源，從而提出了根據(jù)不同污染物通量估算方法的適用條件選擇不同水期內(nèi)相應(yīng)的污染物通量估算方法。

b. 根據(jù)不同廢水來源其污染物削減成本不同，提出了補償中可根據(jù)實際水文條件劃分水期，判定不同水期污染物來源，按照污染物排放量比例加權(quán)平均確定不同污染源削減成本的方法。

通過以上改進(jìn)計算的生態(tài)補償金額比按常規(guī)方法計算的補償金額更加科學(xué)、合理，對上下游的補賠償也會更加公平。因此，按筆者提出的方法實施污染型生態(tài)補償能夠促進(jìn)流域水環(huán)境保護(hù)與水生態(tài)文明城市建設(shè)，保障上下游經(jīng)濟(jì)社會的可持續(xù)發(fā)展。