基于數(shù)理統(tǒng)計的水功能區(qū)達標評價方法

王　佳, 魏文杰, 董　飛,姜國珍，張新華

(1. 四川大學水力學及山區(qū)河流開發(fā)保護國家重點實驗室，四川 成都　610065；

2. 中國水利水電科學研究院水環(huán)境研究所，北京　100038； 3. 四川大學后勤集團，四川 成都 　610065)

?

基于數(shù)理統(tǒng)計的水功能區(qū)達標評價方法

王佳1, 魏文杰1, 董飛2,姜國珍3，張新華1

(1. 四川大學水力學及山區(qū)河流開發(fā)保護國家重點實驗室，四川 成都610065；

2. 中國水利水電科學研究院水環(huán)境研究所，北京100038； 3. 四川大學后勤集團，四川 成都 610065)

摘要：在水功能區(qū)達標評價工作中，針對相同的水質(zhì)數(shù)據(jù)，采用年均值法和頻次法這兩種水功能區(qū)達標評價方法評價，其評價結(jié)果可能截然不同。為了能夠科學、合理地選擇評價方法，筆者引入P-Ⅲ型曲線等數(shù)理統(tǒng)計分析方法，從水質(zhì)監(jiān)測數(shù)據(jù)結(jié)構(gòu)特征進行了分析研究，得到了年均值法和頻次法的適用條件和范圍，當水質(zhì)數(shù)據(jù)滿足≤0.5且≤2采用年均值法，反之采用頻次法。此結(jié)論在大凌河流域和沱江流域的應(yīng)用中得到了驗證，同時解決了不同方法評價所產(chǎn)生評價結(jié)果的矛盾性，為水資源管理規(guī)劃工作提供了有力的依據(jù)。

關(guān)鍵詞：水功能區(qū)；達標評價；頻次法；年均值法；P-Ⅲ型曲線

2011年國務(wù)院批復(fù)了《全國重要湖泊水功能區(qū)劃》，其中水功能區(qū)達標評價方法研究是全國重要湖泊水功能區(qū)評估體系的重要組成部分[1-2]。隨著水資源開發(fā)利用的要求越來越嚴格，以及人們環(huán)保意識的增強,水功能區(qū)的監(jiān)督和管理制度也越來越全面深入，其達標評價方法的選取及評價結(jié)果也備受關(guān)注。

水功能區(qū)達標評價包含單個水功能區(qū)達標評價和流域(區(qū)域)水功能區(qū)達標評價。單個水功能區(qū)達標評價又包括單次水功能區(qū)達標評價、單次水功能區(qū)主要超標項目評價、水期或年度水功能區(qū)達標評價、水期或年度水功能區(qū)主要超標項目評價。

單個水功能區(qū)達標評價應(yīng)以其水質(zhì)管理目標對應(yīng)的水質(zhì)標準為評價標準。目前，水質(zhì)評價方法種類繁多，主要有單因子法、綜合指數(shù)法、模糊綜合法、灰色關(guān)聯(lián)度法、主成分分析法等[3-4]。但每種水質(zhì)評價方法所依托的數(shù)學原理不盡相同，各個方法都具有優(yōu)點和不足之處[5-8]。在實際水功能區(qū)達標評價過程中，水質(zhì)評價方法的選取將直接影響水質(zhì)評價結(jié)果，從而影響水功能區(qū)達標評價的結(jié)果。本研究涉及的水質(zhì)評價方法為單因子法。

單次水功能區(qū)達標評價首先需要根據(jù)規(guī)定的水質(zhì)類別管理目標對水質(zhì)類別進行達標評價。只有所有參評水質(zhì)項目均滿足水質(zhì)類別管理目標要求的水功能區(qū)才達標，有任何一項不滿足水質(zhì)類別要求的水功能區(qū)均不達標。水期或年度水功能區(qū)達標評價應(yīng)在各水功能區(qū)單次達標評價成果基礎(chǔ)上進行。關(guān)于水功能區(qū)達標評價方法暫未形成統(tǒng)一的體系，根據(jù)SL395—2007《地表水資源質(zhì)量評價技術(shù)規(guī)程》及全國重要江河湖泊水功能區(qū)劃，單個水功能區(qū)的水期或年度水功能區(qū)達標評價方法，常采用年均值法和頻次法。然而，對于相同的水質(zhì)數(shù)據(jù)，選取年均值法評價和頻次法評價得出的結(jié)果可能會截然不同，因此研究年均值法和頻次法選取的合理性顯得尤為重要。殷世芳[9]曾對水功能區(qū)達標評價及其影響因素進行了分析，指出評價方法的確定帶有一定的主觀性，科學、合理地確定評價方法需要充分考慮各種客觀影響因素；藍忠華[10]闡述了水功能區(qū)達標評價的原則、標準以及頻次法的定義；喬倩倩等[11]指出了達標評價方法的盲點及其給水環(huán)境評價工作帶來的問題，并提出了改進建議。但是，目前國內(nèi)外對水功能區(qū)達標評價方法的合理性分析和適用范圍的研究還比較少。

本研究應(yīng)用P-Ⅲ型曲線對離散的水質(zhì)數(shù)據(jù)進行擬合，并運用數(shù)理統(tǒng)計等原理和方法，推導單個水功能區(qū)水質(zhì)達標評價的年均值法和頻次法的適用條件與范圍，為科學合理地選擇正確的評價方法提供理論和技術(shù)支撐，同時考慮到中國南方和北方河流可能存在較大差異，選擇北方的大凌河流域及南方的沱江流域作為研究對象，分析驗證本文提出的方法與結(jié)果的合理性與可靠性。

1水功能區(qū)達標評價方法比較

1.1　年均值法

年均值法的水功能區(qū)達標評價基于年度(水期)內(nèi)水質(zhì)指標評價達標結(jié)果，其數(shù)學表達式為

(1)

依次求出該水功能區(qū)全部水質(zhì)指標的平均值，根據(jù)各個指標的平均值來判斷水功能區(qū)是否達標。

1.2　頻次法

基于頻次法的水功能區(qū)達標評價是在各水功能區(qū)單次達標評價成果基礎(chǔ)上進行的。在評價年度(水期)內(nèi)，達標率大于(含等于)80%的水功能區(qū)為年度(水期)達標水功能區(qū)[12]。年度(水期)水功能區(qū)達標率計算式為

(2)

式中：D為年度(水期)水功能區(qū)達標率；G為年度(水期)水功能區(qū)達標次數(shù)；N為水期或年度水功能區(qū)評價的總次數(shù)。

1.3　兩種方法的區(qū)別

年均值法表示年度(水期)內(nèi)水質(zhì)指標的平均情況，它可以反映出某年度(水期)內(nèi)各個水質(zhì)指標整體水平的高低。當水質(zhì)指標數(shù)據(jù)整體的離散程度較小時，年均值可以較好地體現(xiàn)出年度(水期)內(nèi)水質(zhì)指標的整體情況。并且年均值法只需對水質(zhì)指標整體的平均值進行一次評價即可，操作簡單且省時。

頻次法需要對每次監(jiān)測的水質(zhì)數(shù)據(jù)進行水質(zhì)評價，能夠比較完整地描述水質(zhì)變化過程。當水質(zhì)指標數(shù)據(jù)的離散程度較大時，年均值法不能很好地體現(xiàn)和代表該水功能區(qū)的水質(zhì)變化情況，然而頻次法運用達標頻次卻能夠較好地描述水功能區(qū)年度(水期)內(nèi)的變化狀況。

兩種評價方法相互獨立，并未形成統(tǒng)一的評價體系，實際使用時，對于相同的水質(zhì)監(jiān)測數(shù)據(jù)，年均值法和頻次法的評價結(jié)果可能不同，在達標評價方法的選擇中，具有一定的主觀性和隨意性。

2水功能區(qū)達標評價方法的適用范圍分析

2.1　達標評價方法選取的決定因素

根據(jù)上述兩種達標評價方法的特點及差異，可依據(jù)水質(zhì)數(shù)據(jù)的離散程度來確定單個水功能區(qū)的達標評價方法。故選用偏差系數(shù)Cv、偏態(tài)系數(shù)Cs作為達標評價方法選取的主要因素。

2.2　水質(zhì)擬合曲線的選擇

實際的水質(zhì)監(jiān)測工作中測得的均為離散的水質(zhì)數(shù)據(jù)，為了方便研究，需要對離散的水質(zhì)數(shù)據(jù)進行擬合，最終通過對擬合函數(shù)曲線的特點分析年均值法和頻次法的適用條件。

P-Ⅲ型曲線常用于我國的水文現(xiàn)象研究與應(yīng)用中，然而河流水質(zhì)主要受這些水文條件的驅(qū)動和影響，并隨不同水文情況發(fā)生改變，水質(zhì)監(jiān)測數(shù)據(jù)服從P-Ⅲ分布已有學者進行了充分論證[13-16]。因此，本文選用P-Ⅲ型曲線對水質(zhì)數(shù)據(jù)進行擬合。

圖1　P-Ⅲ型密度曲線形狀變化

2.3　達標評價方法適用范圍的推導

首先推導年均值法的適用范圍。為此，從概率密度函數(shù)曲線的特性和概率密度函數(shù)本身兩個方面進行分析。

2.3.1從概率密度函數(shù)曲線的特性考慮

P-Ⅲ型密度曲線的形狀主要取決于參數(shù)Cs，當Cs>0時，可分為圖1所示的4種情況[17-18]：

a. 當2

b. 當Cs=2時，為指數(shù)型密度曲線，仍呈“乙”形，但左端交于x=a0，右端趨近于x軸，如圖1(b)；

圖2　Cs對密度曲線的影響

2.3.2從概率密度函數(shù)的特性考慮

年均值法要求均值能夠較好地代表整體數(shù)據(jù)的情況，因此選取概率密度函數(shù)的極值點對應(yīng)的樣本值(正偏時為xf,負偏時為xg)作為研究對象，下面簡述xf和xg的推導過程。

正偏時：

(3)

負偏時：

(4)

式中：Γ(α)為α的伽瑪函數(shù)；x為樣本濃度；α、β、a0分別為P-Ⅲ型分布的形狀、尺度和位置參數(shù)α>0，β>0。

若函數(shù)取極值時，對應(yīng)的樣本值有：

由曲線圖和參數(shù)關(guān)系可知：

(5)

(6)

(7)

可得：

(8)

(9)

下面針對Cv、Cs、λ值的確定進行探討。

2.4　參數(shù)估計

上述方法所得出的適用條件涉及P-Ⅲ型密度曲線中包含的Cv、Cs。由于水質(zhì)情況是連續(xù)不斷變化的，但實際水質(zhì)觀測是定時觀測，因此監(jiān)測工作是用有限的樣本數(shù)據(jù)去估計總體分布線型中的參數(shù)，故稱為參數(shù)估計。

如何合理地估計參數(shù)，將直接影響水質(zhì)達標評價的結(jié)果，參數(shù)估計在達標評價中至關(guān)重要。目前，樣本估計總體參數(shù)的方法很多，由于水功能區(qū)達標評價的觀測樣本容量較小，為了提高Cs的精度，因此選取權(quán)函數(shù)法計算P-Ⅲ型頻率曲線參數(shù)Cs，選取矩法估計計算參數(shù)Cv，詳細計算可見有關(guān)參考文獻。

2.5　分界點λ值的確定方法

3單個水功能區(qū)達標評價方法選取原則

綜合上述推導分析，可以得出流域水功能區(qū)達標評價的年均值法和頻次法的適用范圍及確定原則：

4實例分析

4.1　以北方的大凌河流域為例

將上述研究應(yīng)用于實際水功能區(qū)的水質(zhì)監(jiān)測數(shù)據(jù)時，需要考慮多個水質(zhì)指標影響下的水功能區(qū)達標評價。以遼寧省大凌河宮山嘴水庫飲用水源區(qū)為例，首先利用該水功能區(qū)2007年的8項水質(zhì)指標監(jiān)測結(jié)果(表1)，分析驗證本文提出的水功能區(qū)達標評價方法的選取原則是否合理。評價流程為：首先運用單因子法進行水質(zhì)評價，基于水質(zhì)評價結(jié)果(是否達到水質(zhì)管理目標)進行分析討論水功能區(qū)達標評價情況。

由表1可見，該水功能區(qū)主要的污染物為COD。其中COD在1月、3月、4月、7月的監(jiān)測結(jié)果都超標，達標頻率僅67%，故用頻次法評價該水功能區(qū)不達標；然而，該水功能區(qū)的8項水質(zhì)指標的年

表1　2007年大凌河宮山嘴水庫飲用水源區(qū)水質(zhì)及達標評價

注：*為所測得的水質(zhì)指標超標；

圖3　大凌河值數(shù)量分布函數(shù)

4.2　以南方的沱江流域為例

圖4　沱江流域值數(shù)量分布函數(shù)

表2　2009年沱江金堂段飲用水源區(qū)水質(zhì)及達標評價

注：*為所測得的水質(zhì)指標超標；

5分析討論

5.1　評價方法適用范圍確定的優(yōu)勢

5.2　評價方法選取的合理性討論

5.3　關(guān)于λ值的確定

6結(jié)論與建議

然而λ值的確定需要大量的水質(zhì)數(shù)據(jù)，筆者僅以南北河流的兩個典型流域為研究對象，提出了λ=0.5，鑒于作者掌握的資料有限，該值還有待于在更多流域進行驗證。λ值確定后，即可在每個流域內(nèi)根據(jù)頻次法和年均值法的適用條件，科學合理地選取水功能區(qū)達標評價方法。

參考文獻：

[ 1 ] 水利部. 水功能區(qū)管理辦法[J]. 中國水利, 2003(7):61-62.(The Ministry of Water Resources.Measures for the management of water function areas[J].China Water Resources,2003(7):61-62.(in Chinese))

[ 2 ] 張學峰,尚曉成,張軍獻. 水功能區(qū)管理淺談[J].人民黃河,2004,26(4):21-23.(ZHANG Xuefeng, SHANG Xiaocheng, ZHANG Junxian. Introduction to water management[J]. Yellow River,2004,26(4):21-23.(in Chinese))

[ 3 ] 朱靜平. 幾種水環(huán)境質(zhì)量綜合評價方法的探討[J].西南科技大學學報,2002,17(4):62-67.(ZHU Jingping. Comparison of several comprehensive assessment methods of water environmental quality[J].Journal of Southwest University of Science and Technology,2002,17(4):62-67.(in Chinese))

[ 4 ] 黃興國,劉秀花. 水環(huán)境質(zhì)量評價中幾種方法的對比[J]．地下水，2005，27(2):125 -126．(HUANG Xingguo, LIU Xiuhua. Comparison on several methods of quality evaluation on water environment[J]. Ground Water,2005,27(2):125-126.(in Chinese))

[ 5 ] 劉臣,王淑文. 水環(huán)境質(zhì)量評價三種方法應(yīng)用淺析[J].東北水利水電,2004(6):44-46.(LIU Chen, WANG Shuwen. Superficial analysis of 3 methods for water environment quality evaluation[J]. Water Resources & Hydropower of Northeast China,2004(6):44-46.(in Chinese))

[ 6 ] 毛興華. 常用水質(zhì)評價方法的選擇[J].水科學與工程技術(shù),2006(1):22-23.(MAO Xinghua. Selection on assessment methods of typical water quality[J]. Water Science and Engineering Technology,2006(1):22-23.(in Chinese))

[ 7 ] 賴曉珍,張炎齋,馮英. 關(guān)于水功能區(qū)水質(zhì)評價方法的分析[J]. 調(diào)查與評價, 2009(12):8-9.(LAI Xiaozhen, ZHANG yanzhai, FENG ying. Analysis of water function area water quality evaluation methods[J]. Research and Assessment , 2009(12):8-9. (in Chinese))

[ 8 ] 劉發(fā)根,郭玉銀. 一種水功能區(qū)水質(zhì)達標評價的新方法 [J].人民長江, 2014(18): 28-32.(LIU Fagen,GUO Yuyin. A new water quality up to standard assessment methodology for water function zone[J]. Yangtze River, 2014(18):28-32.(in Chinese))

[ 9 ] 殷世芳. 水功能區(qū)達標率評價及影響因素分析[J]. 人民黃河,2012, 34(5):38-39.( YIN Shifang. Functional areas of water compliance rate evaluation and its influencing factors[J]. Yellow River, 2012, 34(5):38-39.(in Chinese))

[10] 藍忠華. 水功能區(qū)達標評價的方法研究[J]. 中國新科技新產(chǎn)品, 2010(2):1.(LAN Zhonghua. Research on water function area evaluation[J].China New Technologies and Products, 2010(2):1.(in Chinese))

[11] 喬倩倩,許鑫,駱素娜. 水功能區(qū)水質(zhì)達標評價方法分析[J]. 東北水利水電, 2013(8):5-7. (QIAO Qianqian, XU Xin, LUO Suna. Analysis of evaluation method of water quality reaching standard for water function zones[J]. Water Resources & Hydropower of Northeast China, 1988, 1(4):1-7. (in Chinese))

[12] SL395—2007地表水資源質(zhì)量評價技術(shù)規(guī)程[S].

[13] 辛江. 試用概率方法進行大遼河水質(zhì)評價[J]. 東北水利水電,2001(8):44-46.(XIN Jiang. Daliao River water quality assessment probability method[J]. Water Resources & Hydropower of Northeast China,2001(8):44-46. (in Chinese))

[14] 劉秀花,黃興國.基于皮爾遜Ⅲ型曲線的渭河水質(zhì)歷時評價研究[J].灌溉排水學報,2004,23(6):49-51.(LIU Xiuhua, HUANG Xingguo. Evaluation of duration of pollution based on the distribution curve of Pearson Ⅲ in the Shanxi Section of Weihe River [J].Journal of Irrigation and Drainage,2004,23(6):49-51. ( in Chinese))

[15] 喬明葉,周進,任金亮.皮爾遜Ⅲ型曲線在灤河水質(zhì)歷時評價中的應(yīng)用[J].人民黃河,2012(5):61-62. (QIAO Mingye,ZHOU Jin,REN Jinliang. Evaluation of duration of pollution of Luanhe River based on distribution curve of Pearson Ⅲ [J]. Yellow River，2012(5):61-62. (in Chinese))

[16] 姚允龍. 頻率曲線法及皮爾遜Ⅲ型曲線在水質(zhì)評價中的應(yīng)用探討[J]. 水文, 2001, 21(1):42-46. (YAO Yunlong. Frequency curve and Pearson Ⅲ curve applied in water quality evaluation[J]. Journal of China Hydrology, 2001, 21(1):42-46. (in Chinese))

[17] 葉守澤, 詹道江. 工程水文學[M]. 3版.北京： 中國水利水電出版社, 2007.

[18] 葉守澤. 水文水利計算[M]. 北京：中國水利水電出版社, 1992.

[19] 蔡體錄. 負偏頻率曲線的計算[J]. 東海海洋, 1988, 1(4):1-7. (CAI Tilu. Calculation of egative bias frequency curve[J]. Donghai Marine Science, 1988, 1(4):1-7. (in Chinese))

[20] 李興拼,鄭江麗. 淺談P-Ⅲ型負偏頻率曲線計算方法[J]. 廣東水利水電, 2010(9):17-18. (LI Xingpin, ZHENG Jiangli. Study on the calculation method of Pearson Ⅲ frequency curve in negative bias station[J]. Guangdong Water Resources and Hydropower, 2010(9):17-18. (in Chinese))

[21] 劉天厚,施為光. 沱江水環(huán)境容量及水質(zhì)管理規(guī)劃探討[J]. 上海環(huán)境科學, 1987(9):17-18.(LIU Tianhou, SHI Weiguang. Study on environment capacity and its application of Tuojiang River[J]. Shanghai Environmental Science, 1987(9):17-18.( in Chinese))

Research on water function area evaluation based on mathematical statistics method

WANG Jia1, WEI Wenjie1, DONG Fei2, JIANG Guozhen3, ZHANG Xinhua1

(1.StateKeyLaboratoryofHydraulicsandMountainRiverEngineering,Chengdu610065,China;

2.ChinaInstituteofWaterResourcesandHydropowerResearch,WaterEnvironment

ResearchInstitute,Beijing100038,China;

3.LogisticsGroup,SichuanUniversity,Chengdu610065,China)

Abstract：Different assessment results may be obtained through annual average method and frequency method evaluating water function zones with same water quality data, during water function area compliance evaluation work. In order to choose the evaluation method scientifically and reasonably, Mathematical statistics analysis methods including Pearson Ⅲ curve was introduced and the characteristic of water quality monitoring data was analyzed. In the end, the applying conditions and range of annual average method and frequency method are acquired. When the water quality data meeting≤2, annual average method should be applied, otherwise frequency method should be chosen. Two case studies have been conducted in the northern Daling River basin and the southern Tuojiang River basin, respectively, to verify the conclusion. The contradiction between the evaluation results of different methods is solved, providing strong basis for water resource management and planning.

Key words：water function area; compliance evaluation; frequency method; annual average method; Pearson Ⅲ curve

(收稿日期：2015-01-05編輯：徐娟)

中圖分類號：TV212.4

文獻標志碼：A

文章編號：1004-6933(2016)01-0154-07

作者簡介：王佳(1990—)，女，碩士研究生，研究方向為水環(huán)境及水信息。E-mail: wozhixiwang@126.com通信作者：張新華，教授。E-mail: xhzhang@scu.edu.cn

項目基金：中央分成水資源費項目(1261330111042)；國家自然科學基金(51379137)

DOI：10.3880/j.issn.1004-6933.2016.01.027