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

(1. 四川大學(xué)水力學(xué)及山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，四川 成都　610065；

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

?

基于數(shù)理統(tǒng)計(jì)的水功能區(qū)達(dá)標(biāo)評價(jià)方法

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

(1. 四川大學(xué)水力學(xué)及山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，四川 成都610065；

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

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

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

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

水功能區(qū)達(dá)標(biāo)評價(jià)包含單個水功能區(qū)達(dá)標(biāo)評價(jià)和流域(區(qū)域)水功能區(qū)達(dá)標(biāo)評價(jià)。單個水功能區(qū)達(dá)標(biāo)評價(jià)又包括單次水功能區(qū)達(dá)標(biāo)評價(jià)、單次水功能區(qū)主要超標(biāo)項(xiàng)目評價(jià)、水期或年度水功能區(qū)達(dá)標(biāo)評價(jià)、水期或年度水功能區(qū)主要超標(biāo)項(xiàng)目評價(jià)。

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

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

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

1水功能區(qū)達(dá)標(biāo)評價(jià)方法比較

1.1　年均值法

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

(1)

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

1.2　頻次法

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

(2)

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

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

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

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

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

2水功能區(qū)達(dá)標(biāo)評價(jià)方法的適用范圍分析

2.1　達(dá)標(biāo)評價(jià)方法選取的決定因素

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

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

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

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

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

2.3　達(dá)標(biāo)評價(jià)方法適用范圍的推導(dǎo)

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

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

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

a. 當(dāng)2

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

圖2　Cs對密度曲線的影響

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

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

正偏時：

(3)

負(fù)偏時：

(4)

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

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

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

(5)

(6)

(7)

可得：

(8)

(9)

下面針對Cv、Cs、λ值的確定進(jìn)行探討。

2.4　參數(shù)估計(jì)

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

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

2.5　分界點(diǎn)λ值的確定方法

3單個水功能區(qū)達(dá)標(biāo)評價(jià)方法選取原則

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

4實(shí)例分析

4.1　以北方的大凌河流域?yàn)槔?/h2>

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

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

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

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

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

4.2　以南方的沱江流域?yàn)槔?/h2>

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

表2　2009年沱江金堂段飲用水源區(qū)水質(zhì)及達(dá)標(biāo)評價(jià)

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

5分析討論

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

5.2　評價(jià)方法選取的合理性討論

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

6結(jié)論與建議

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

參考文獻(xiàn)：

[ 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 ] 張學(xué)峰,尚曉成,張軍獻(xiàn). 水功能區(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ì)量綜合評價(jià)方法的探討[J].西南科技大學(xué)學(xué)報(bào),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ì)量評價(jià)中幾種方法的對比[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ì)量評價(jià)三種方法應(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ì)評價(jià)方法的選擇[J].水科學(xué)與工程技術(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ì)評價(jià)方法的分析[J]. 調(diào)查與評價(jià), 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ì)達(dá)標(biāo)評價(jià)的新方法 [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ū)達(dá)標(biāo)率評價(jià)及影響因素分析[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] 藍(lán)忠華. 水功能區(qū)達(dá)標(biāo)評價(jià)的方法研究[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] 喬倩倩,許鑫,駱?biāo)啬? 水功能區(qū)水質(zhì)達(dá)標(biāo)評價(jià)方法分析[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ì)量評價(jià)技術(shù)規(guī)程[S].

[13] 辛江. 試用概率方法進(jìn)行大遼河水質(zhì)評價(jià)[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ì)歷時評價(jià)研究[J].灌溉排水學(xué)報(bào),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] 喬明葉,周進(jìn),任金亮.皮爾遜Ⅲ型曲線在灤河水質(zhì)歷時評價(jià)中的應(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ì)評價(jià)中的應(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] 葉守澤, 詹道江. 工程水文學(xué)[M]. 3版.北京： 中國水利水電出版社, 2007.

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

[19] 蔡體錄. 負(fù)偏頻率曲線的計(jì)算[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-Ⅲ型負(fù)偏頻率曲線計(jì)算方法[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)境科學(xué), 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

文獻(xiàn)標(biāo)志碼：A

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

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

項(xiàng)目基金：中央分成水資源費(fèi)項(xiàng)目(1261330111042)；國家自然科學(xué)基金(51379137)

DOI：10.3880/j.issn.1004-6933.2016.01.027