1 引言

　　水質(zhì)達(dá)標(biāo)評(píng)價(jià)(Water Quality Compliance Assessment)是根據(jù)一定時(shí)段內(nèi)的水質(zhì)監(jiān)測數(shù)據(jù)，評(píng)價(jià)水體狀態(tài)是否符合其功能對(duì)應(yīng)的水質(zhì)標(biāo)準(zhǔn)的決策過程，是流域水污染防治和水質(zhì)管理的重要科學(xué)依據(jù).區(qū)別于水質(zhì)評(píng)價(jià)研究中對(duì)水質(zhì)綜合指標(biāo)、時(shí)間序列趨勢、季節(jié)性特征等的關(guān)注，水質(zhì)達(dá)標(biāo)評(píng)價(jià)的關(guān)注點(diǎn)是水質(zhì)是否達(dá)標(biāo)，并可以由此評(píng)估已有的流域治理措施效果、制定新的治理規(guī)劃(如TMDL).在流域水污染防治資金緊缺的條件下，準(zhǔn)確識(shí)別水體是否達(dá)標(biāo)至關(guān)重要：如果評(píng)價(jià)方法過于保守，會(huì)導(dǎo)致水生態(tài)系統(tǒng)健康受損;若評(píng)價(jià)方法過于嚴(yán)格，會(huì)導(dǎo)致水體過于保護(hù)，水體自凈能力無法得到有效利用及高的治理資金投入.科學(xué)合理的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法是進(jìn)行水質(zhì)達(dá)標(biāo)評(píng)價(jià)的基礎(chǔ).我國地表水環(huán)境質(zhì)量評(píng)價(jià)采用單因子評(píng)價(jià)方法，對(duì)斷面時(shí)空監(jiān)測數(shù)據(jù)求算術(shù)平均值判斷其水質(zhì)類別，根據(jù)斷面水質(zhì)類別比例對(duì)流域水質(zhì)進(jìn)行分級(jí).美國清潔水法303(d)條款規(guī)定對(duì)受損水體實(shí)施TMDL計(jì)劃，EPA導(dǎo)則推薦使用監(jiān)測數(shù)據(jù)的90%分位數(shù)評(píng)價(jià)水體是否受損，即監(jiān)測數(shù)據(jù)的超標(biāo)個(gè)數(shù)不能超過10%.在歐盟，采用每年12個(gè)月監(jiān)測數(shù)據(jù)的算術(shù)平均值作為評(píng)價(jià)依據(jù)，少數(shù)成員國將監(jiān)測數(shù)據(jù)的變異性考慮在內(nèi)，給出統(tǒng)計(jì)置信區(qū)間.

　　水質(zhì)指標(biāo)監(jiān)測值存在不確定性已經(jīng)被廣泛地認(rèn)識(shí)和研究.我國水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法未考慮水質(zhì)指標(biāo)的不確定性，美國和歐盟已經(jīng)將不確定性納入到水質(zhì)達(dá)標(biāo)評(píng)價(jià)中.水質(zhì)指標(biāo)的不確定性包括固有的變異性(如監(jiān)測數(shù)據(jù)的季節(jié)性差異)和可減小的不確定性(如測量誤差)2部分.通過合理選擇監(jiān)測斷面、監(jiān)測指標(biāo)以及質(zhì)量保證和質(zhì)量控制(QA/QC)可以減小不確定性，卻不能完全消除.忽視不確定性可能會(huì)帶來巨大的決策風(fēng)險(xiǎn)，科學(xué)合理的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法必須能夠表征監(jiān)測數(shù)據(jù)的不確定性，將由水質(zhì)指標(biāo)不確定性帶來的決策風(fēng)險(xiǎn)定量化，為決策提供有用信息.

　　在統(tǒng)計(jì)學(xué)視角下，將水質(zhì)指標(biāo)視為隨機(jī)變量，可有效地表征不確定性.水質(zhì)指標(biāo)具有時(shí)空分異性，在同一時(shí)間獲得空間所有點(diǎn)位的監(jiān)測數(shù)據(jù)或者在同一點(diǎn)位進(jìn)行連續(xù)的監(jiān)測是不可能的，只能通過有限的監(jiān)測數(shù)據(jù)對(duì)水質(zhì)進(jìn)行評(píng)價(jià)，即根據(jù)樣本數(shù)據(jù)估計(jì)總體特征.統(tǒng)計(jì)學(xué)假設(shè)檢驗(yàn)方法已經(jīng)廣泛應(yīng)用于環(huán)境系統(tǒng)中，假設(shè)檢驗(yàn)的2類錯(cuò)誤可以定量表征水質(zhì)指標(biāo)的不確定性帶來的決策風(fēng)險(xiǎn).針對(duì)美國EPA提出的最大超標(biāo)率為10%，Smith等和McBride等認(rèn)為10%應(yīng)該指10%的時(shí)間，即總體分布的90%分位數(shù)，并采用二項(xiàng)分布檢驗(yàn)法(Binomial Test)進(jìn)行了水質(zhì)達(dá)標(biāo)評(píng)價(jià).該方法是一種針對(duì)二元數(shù)據(jù)的非參數(shù)方法，通過假設(shè)檢驗(yàn)給出在一定允許超標(biāo)率和樣本容量時(shí)的最大超標(biāo)個(gè)數(shù)，并與實(shí)際監(jiān)測得到的超標(biāo)個(gè)數(shù)對(duì)比，評(píng)價(jià)水質(zhì)是否超標(biāo).盡管研究者指出二項(xiàng)分布檢驗(yàn)法存在信息損失的問題，并提出了參數(shù)檢驗(yàn)方法，然而McBride的研究表明，在大部分情況下，二項(xiàng)分布檢驗(yàn)法由于不受限于原始數(shù)據(jù)的分布而更具有魯棒性.本文以基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法為基礎(chǔ)，分析其決策過程及決策風(fēng)險(xiǎn)的定量表征，并以海河流域?yàn)槔�，�?duì)比不同評(píng)價(jià)方法的差異性，分析其適用性和靈活性.

　　2 研究方法

　　2.1 基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法

　　在二項(xiàng)分布檢驗(yàn)法之前，需將監(jiān)測數(shù)據(jù)轉(zhuǎn)化為0~1變量，方法為：對(duì)于某一特定監(jiān)測數(shù)據(jù)x，將其與水質(zhì)標(biāo)準(zhǔn)θ對(duì)比，對(duì)比后的結(jié)果記為y，如果x>θ，即超過水質(zhì)標(biāo)準(zhǔn)，記y=1;如果x≤θ，記y=0.令p(0θ的概率，則y服從成功概率為p的0~1分布，其分布律為：

　　假設(shè)收集到的N個(gè)監(jiān)測數(shù)據(jù)用于評(píng)價(jià)水質(zhì)是否達(dá)標(biāo)，記隨機(jī)變量X=x1，x2，...，xN，將X中的每個(gè)元素轉(zhuǎn)化為0~1變量，便得到隨機(jī)變量Y=y1，y2，...，yN，則Y服從參數(shù)為(N，p)的二項(xiàng)分布，記為Y~B(N，p).令n為Y成功的總次數(shù)，即為N個(gè)樣品中超標(biāo)總個(gè)數(shù)，其分布律為：

　　二項(xiàng)分布的分布律表征了水質(zhì)變量轉(zhuǎn)化為二元數(shù)據(jù)后超標(biāo)總個(gè)數(shù)不確定性.欲檢驗(yàn)總體分布的90%分位數(shù)是否達(dá)標(biāo)，即超標(biāo)率是否小于10%，可令原假設(shè)為水質(zhì)達(dá)標(biāo)，即H0:p≤0.1，備擇假設(shè)為H1:p>0.1.原假設(shè)為真時(shí)，可得p=0.1時(shí)超標(biāo)個(gè)數(shù)n的分布函數(shù)F;實(shí)際監(jiān)測數(shù)據(jù)中超標(biāo)個(gè)數(shù)nc不應(yīng)太大，在給定的顯著性水平α下，應(yīng)有F(nc)=P(n≤nc)≤1-α;否則拒絕原假設(shè)，判定為水質(zhì)超標(biāo).

　　2.2 棄真錯(cuò)誤和取偽錯(cuò)誤

　　在經(jīng)典統(tǒng)計(jì)學(xué)的假設(shè)檢驗(yàn)過程中，不可避免地要討論2類錯(cuò)誤：棄真錯(cuò)誤和取偽錯(cuò)誤.在水質(zhì)達(dá)標(biāo)評(píng)價(jià)中，這2類錯(cuò)誤是由于水質(zhì)變量的不確定性導(dǎo)致的.棄真錯(cuò)誤為當(dāng)原假設(shè)為真時(shí)，放棄原假設(shè)而犯的錯(cuò)誤(falsely inferring a breach of standard)，棄真錯(cuò)誤的概率記為α;取偽錯(cuò)誤為當(dāng)原假設(shè)為假時(shí)，相信原假設(shè)而犯的錯(cuò)誤(falsely inferring compliance)，棄偽錯(cuò)誤的概率記為β.對(duì)于二項(xiàng)分布，$\alpha =P\left( n>{{n}_{s}} \right)=\sum\limits_{j={{n}_{s}}+1}^{N}{C_{N}^{j}{{P}^{j}}}{{\left( 1-P \right)}^{N-j}}$.計(jì)算β時(shí)，需要給定效應(yīng)值η(備擇假設(shè)和原假設(shè)之間的差距)，二項(xiàng)分布檢驗(yàn)法一般取η=0.15，即p′=p+η=0.25來計(jì)算β，計(jì)算公式為β=G-1(ns)，其中G(n)=$\sum\limits_{k=1}^{N}{C_{N}^{k}{{P}^{'k}}}{{\left( 1-P' \right)}^{N-k}}$，為成功概率為p′的二項(xiàng)分布的分布函數(shù)，G-1(n)為G(n)的反函數(shù).

　　按照上述α和β的計(jì)算公式，可以得到N一定時(shí)，以不同的ns(0≤ns≤N)作為決策準(zhǔn)則(即nc>ns判定為超標(biāo)，nc≤ns判定為達(dá)標(biāo))時(shí)，α和β的變化情況.由圖 1可知N一定時(shí)，對(duì)于特定的ns存在唯一的(α，β)組合;α隨ns減小，β隨ns增加.圖中豎直虛線假設(shè)變量服從對(duì)稱分布，算術(shù)平均值等于50%分位數(shù)，則平均值法等價(jià)于最大允許超標(biāo)率為50%.圖中豎直實(shí)線代表采用10%超標(biāo)樣品數(shù)目(以下稱為“超標(biāo)比例法”)得到的ns，當(dāng)N=12時(shí)，ns=1，α=0.341;當(dāng)N=50時(shí)，ns=5，α=0.384.可見對(duì)于相同樣本容量，平均值法棄偽錯(cuò)誤的概率很大，超標(biāo)比例法棄真錯(cuò)誤的概率較大，而二項(xiàng)分布檢驗(yàn)法則可根據(jù)決策者需要選擇不同的(α，β)組合.相對(duì)于選擇特定的顯著性水平值做決策，(α，β)組合能量化決策準(zhǔn)則的風(fēng)險(xiǎn)，為決策者提供更大的決策空間.α和β隨N增加均有減小的總體趨勢，參見Smith等.

　　圖 1N=12(a)和50(b)時(shí)2類錯(cuò)誤概率隨最大允許超標(biāo)個(gè)數(shù)的變化

　　2.3 決策過程

　　基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)的決策過程包括2個(gè)環(huán)節(jié)：首先確定ns，然后根據(jù)監(jiān)測樣品的超標(biāo)個(gè)數(shù)nc判斷水質(zhì)是否達(dá)標(biāo).一般而言，決策過程應(yīng)該使得α和β最小.在N一定時(shí)，由于α、β隨ns單調(diào)性相反(圖 1)，對(duì)于特定的ns不能同時(shí)使α和β最小.作為決策者，在選擇ns時(shí)：一方面，需要根據(jù)風(fēng)險(xiǎn)偏好和科學(xué)判斷，選擇最大允許超標(biāo)率p;另一方面，需要權(quán)衡各相關(guān)方的利益，如：作為受體人群和生態(tài)系統(tǒng)的及作為治理投資者的政府或企業(yè).在同等條件下，ns隨著p的增大而不嚴(yán)格單調(diào)遞增;在p和N一定時(shí)，根據(jù)決策者的利益權(quán)衡，可以得到不同的ns.

　　下面介紹3種典型的(α，β)組合選擇策略：①保證α在一定的低水平，例如α≤0.05，使得β最小;②保證β在一定的低水平，例如β≤0.05，使得α最小;③控制α≈β，即minα-β.圖 2給出了當(dāng)p=0.1，N為4~100時(shí)，這3種策略得到的最大允許超標(biāo)個(gè)數(shù)ns和最大允許超標(biāo)率η=nsN(散點(diǎn)圖，連續(xù)線為點(diǎn)過于密集的結(jié)果).可見當(dāng)N較小時(shí)，控制α≤0.05得到的ns最大，控制β≤0.05得到的ns最小，且當(dāng)N≤10時(shí)，無論如何選擇ns均不能使β≤0.05;隨著N變大，控制α≤0.05得到的ns最小，控制β≤0.05得到的ns最大;策略③的結(jié)果介于①和②之間.總體上，控制α≤0.05得到的η呈遞減趨勢，由于假設(shè)檢驗(yàn)是p=0.1時(shí)二項(xiàng)分布的95%分位數(shù)，因此η恒大于0.1(圖 2d虛線);控制β≤0.05得到的η呈遞增趨勢，由于假設(shè)檢驗(yàn)是p=0.25時(shí)二項(xiàng)分布的5%分位數(shù)，因此η恒小于0.25(圖 2e虛線).綜上，基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)決策過程可表示為圖 3.決策者可以根據(jù)風(fēng)險(xiǎn)偏好和利益權(quán)衡，選擇p值和(α，β)組合，進(jìn)而確定N一定時(shí)的ns，結(jié)合水質(zhì)監(jiān)測數(shù)據(jù)中的實(shí)際超標(biāo)個(gè)數(shù)nc，對(duì)水體的達(dá)標(biāo)狀況進(jìn)行判斷;通過p和(α，β)可將決策的各種風(fēng)險(xiǎn)定量化.

　　圖 2(Fig. 2)

　　圖 2最大允許超標(biāo)個(gè)數(shù)和最大允許超標(biāo)率隨樣本容量的變化

　　圖 3基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)決策過程

　　3 案例研究(Case study)3.1 研究對(duì)象和方法

　　選擇海河流域的CODMn作為達(dá)標(biāo)評(píng)價(jià)的對(duì)象，數(shù)據(jù)來源為環(huán)保部網(wǎng)站，收集2014年7月22日-2015年7月14日的周監(jiān)測數(shù)據(jù)，共有北京密云古北口(51)、北京門頭溝沿河城(51)、天津三岔口(51)、天津果河橋(51)、河北張家口八號(hào)橋(45)、河北石家莊崗南水庫(51)、山東聊城秤鉤灣(42)等7個(gè)監(jiān)測斷面，括號(hào)中為不同監(jiān)測斷面的有效監(jiān)測樣品數(shù).

　　當(dāng)N=42、45、51時(shí)，采用平均值法、二項(xiàng)分布檢驗(yàn)法的3種策略、超標(biāo)比例法確定α、β、ns(對(duì)于平均值法，求各個(gè)樣本中大于平均值的樣品個(gè)數(shù)，作為平均值法的ns)，對(duì)比監(jiān)測數(shù)據(jù)與地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)中各水質(zhì)類別的濃度限值θi(i=1，2，...，5)，得到超過θi的樣品個(gè)數(shù)nci，水質(zhì)類別判定遵循如下過程：當(dāng)i=1時(shí)，若nc1≤ns則水質(zhì)類別判定為I類;當(dāng)i=5時(shí)，若nc5>ns則判定為劣V類;若nci≤ns

　　3.2 水質(zhì)評(píng)價(jià)結(jié)果

　　各種方法對(duì)應(yīng)的α、β、ns結(jié)果見表 1，各種方法的嚴(yán)格程度為：超標(biāo)比例法>B(2)>B(3)>B(1)>平均值法(“a>b”表示a嚴(yán)格于b，B(i)表示二項(xiàng)分布檢驗(yàn)法的第i種策略);采用平均值法β接近于1，超標(biāo)比例法得到的α很大(>0.4)，二項(xiàng)分布檢驗(yàn)法的3種策略得到的2類錯(cuò)誤之和小于0.3，平均在0.2左右，且均處于相對(duì)較低水平.水質(zhì)評(píng)價(jià)結(jié)果見表 2，在7個(gè)斷面中，石家莊崗南水庫水質(zhì)最好，聊城秤鉤灣水質(zhì)最差;石家莊崗南水庫和門頭溝沿河城2個(gè)斷面采用各種方法的評(píng)價(jià)結(jié)果相同，另外5個(gè)斷面有差異.根據(jù)斷面的水質(zhì)類別得到的流域分級(jí)結(jié)果也有差異.斷面水質(zhì)類別的評(píng)價(jià)結(jié)果，除了與評(píng)價(jià)方法的嚴(yán)格程度有關(guān)，也與水質(zhì)監(jiān)測數(shù)據(jù)有關(guān).采用相對(duì)嚴(yán)格方法評(píng)價(jià)達(dá)標(biāo)的斷面，在不嚴(yán)格方法下自然達(dá)標(biāo);采用相對(duì)不嚴(yán)格方法評(píng)價(jià)達(dá)標(biāo)的斷面，在嚴(yán)格的方法下不一定達(dá)標(biāo)，例如密云古北口采用平均值法評(píng)價(jià)為I類水質(zhì)，而采用其它方法評(píng)價(jià)時(shí)卻為II類水質(zhì).總體上看，評(píng)價(jià)結(jié)果符合嚴(yán)格方法得到的水質(zhì)類別較低的規(guī)律.

　　表 1 3種樣本容量對(duì)應(yīng)的各種評(píng)價(jià)方法的2類錯(cuò)誤概率和最大允許超標(biāo)個(gè)數(shù)

　　表 2 海河流域CODMn達(dá)標(biāo)評(píng)價(jià)結(jié)果

　　綜合表 1和表 2，對(duì)于二項(xiàng)分布檢驗(yàn)法的3種方法，控制β≤0.05有利于保護(hù)受體的利益，對(duì)水體水質(zhì)要求嚴(yán)格，流域分級(jí)為“輕度污染”;控制α≤0.05有利于減少治理投資，對(duì)流域分級(jí)為“良好”;二項(xiàng)檢驗(yàn)的3種策略的評(píng)價(jià)結(jié)果對(duì)比，說明在采用二項(xiàng)分布檢驗(yàn)法進(jìn)行水質(zhì)達(dá)標(biāo)評(píng)價(jià)時(shí)，即使對(duì)于完全相同的水質(zhì)監(jiān)測數(shù)據(jù)，決策時(shí)傾向于不同利益群體，可能得出不同的水質(zhì)達(dá)標(biāo)評(píng)價(jià)結(jié)果.二項(xiàng)分布檢驗(yàn)法為水質(zhì)達(dá)標(biāo)評(píng)價(jià)的決策提供了較為廣闊的決策空間.

　　案例研究中，取最大超標(biāo)率p=0.1，采用二項(xiàng)分布檢驗(yàn)法的3種典型策略進(jìn)行水質(zhì)達(dá)標(biāo)評(píng)價(jià);在實(shí)踐中，決策過程有更大的靈活性：對(duì)于p，決策者可以根據(jù)風(fēng)險(xiǎn)偏好選擇不同值.對(duì)于策略①和②，決策可設(shè)定不同的錯(cuò)誤概率限值，例如令α≤0.1使得β最小，再比如，可同時(shí)限定α，β的最大值，例如令α≤0.2同時(shí)β≤0.1，即可得到多組(α，β)組合;對(duì)于策略③，可根據(jù)需要設(shè)定限制條件，例如出于保護(hù)受體的目的令α<β，或者給定最大可接受的差值δ，令α-β≤δ，亦可得到不同的(α，β)組合.對(duì)得到的多組(α，β)組合，決策者可結(jié)合其它約束條件進(jìn)行優(yōu)選確定最終的ns.

　　4 討論

　　4.1 先驗(yàn)知識(shí)與假設(shè)檢驗(yàn)

　　研究中采用的原假設(shè)是水質(zhì)達(dá)標(biāo)，α代表水質(zhì)達(dá)標(biāo)而判定為不達(dá)標(biāo)的概率，β代表水質(zhì)未達(dá)標(biāo)而判定為達(dá)標(biāo)的概率，二項(xiàng)分布檢驗(yàn)法的策略①實(shí)質(zhì)上是求成功概率為p的二項(xiàng)分布的95%分位數(shù);Gibbons指出如果有很強(qiáng)的先驗(yàn)知識(shí)認(rèn)為水體超標(biāo)，則假設(shè)檢驗(yàn)的原假設(shè)應(yīng)為水質(zhì)未達(dá)標(biāo)，此時(shí)α和β的含義互換，而二項(xiàng)分布檢驗(yàn)法的策略①變?yōu)榍蟪晒Ω怕蕿閜的二項(xiàng)分布的5%分位數(shù).由計(jì)算可得，當(dāng)先驗(yàn)知識(shí)認(rèn)為水質(zhì)超標(biāo)時(shí)，得到的ns較小.如當(dāng)效應(yīng)值為0.15，樣本容量N=50，最大超標(biāo)率p=0.25時(shí)，按策略①，原假設(shè)為水質(zhì)達(dá)標(biāo)時(shí)ns=18，原假設(shè)為水質(zhì)超標(biāo)時(shí)ns=7，則對(duì)于超標(biāo)個(gè)數(shù)在8~18的樣本按照2種原假設(shè)得到的評(píng)價(jià)結(jié)果不同，例如當(dāng)nc=10時(shí)，若原假設(shè)為水質(zhì)達(dá)標(biāo)，則判定為達(dá)標(biāo)，若原假設(shè)為水質(zhì)超標(biāo)，則判定為超標(biāo).可見，水質(zhì)現(xiàn)狀的先驗(yàn)知識(shí)對(duì)評(píng)價(jià)結(jié)果具有顯著影響，在先驗(yàn)知識(shí)充足時(shí)，水質(zhì)達(dá)標(biāo)評(píng)價(jià)具有“寬進(jìn)嚴(yán)出”的規(guī)律，即將先驗(yàn)認(rèn)為達(dá)標(biāo)的水體列為不達(dá)標(biāo)狀態(tài)時(shí)對(duì)超標(biāo)個(gè)數(shù)的要求較為保守，而將先驗(yàn)認(rèn)為未達(dá)標(biāo)水體列為達(dá)標(biāo)狀態(tài)時(shí)對(duì)超標(biāo)個(gè)數(shù)的要求較為嚴(yán)格.β值與效應(yīng)值的選取有關(guān)，本文參照文獻(xiàn)做法取效應(yīng)值為0.15，即備擇假設(shè)的p=0.25.在同等條件下，盡管根據(jù)不同效應(yīng)值可以得到不同的β，但是其隨N和ns的變化規(guī)律一致.假設(shè)檢驗(yàn)的效果隨樣本容量的增加而增大，對(duì)海河流域的研究顯示，采用二項(xiàng)分布檢驗(yàn)法可以控制各個(gè)斷面的總錯(cuò)誤概率平均值在0.2左右，評(píng)價(jià)結(jié)果具有很高的可信度.

　　4.2 最大允許超標(biāo)概率的選擇

　　對(duì)于研究中最大超標(biāo)率設(shè)定的目的，研究者給出了不同的解釋：Smith等、McBride等認(rèn)為水質(zhì)變量不確定性造成的某次監(jiān)測數(shù)據(jù)超標(biāo)是難以避免的，因此對(duì)全部水質(zhì)監(jiān)測數(shù)據(jù)給出一個(gè)可以接受的超標(biāo)率(10%)是合理的;此時(shí)10%可理解為“容忍度”.Gibbons則認(rèn)為對(duì)90%分位數(shù)進(jìn)行假設(shè)檢驗(yàn)是為了保證總體的均值小于標(biāo)準(zhǔn)值;此時(shí)，樣本均值與90%分位數(shù)之間的差異可理解為“安全因子”.Smith等針對(duì)的是一類分位數(shù)標(biāo)準(zhǔn)(percentile-based standards)，而Gibbons針對(duì)的則是平均值標(biāo)準(zhǔn)，平均值代表流域水質(zhì)變量波動(dòng)性的平均水平，與室內(nèi)實(shí)驗(yàn)或野外試驗(yàn)的穩(wěn)定濃度相對(duì)應(yīng).本文建議，在使用二項(xiàng)分布檢驗(yàn)法進(jìn)行水質(zhì)達(dá)標(biāo)評(píng)價(jià)時(shí)，需要首先區(qū)分評(píng)價(jià)變量的標(biāo)準(zhǔn)類型，如果是分位數(shù)標(biāo)準(zhǔn)，則根據(jù)要求選擇較小的p值(例如0.1)是合理的;如果是平均值標(biāo)準(zhǔn)，過大的安全因子會(huì)導(dǎo)致水質(zhì)目標(biāo)難以實(shí)現(xiàn)，需要選擇較大的p值(例如0.25).

　　結(jié)合4.1和4.2小節(jié)，最大允許超標(biāo)概率和假設(shè)檢驗(yàn)的2類錯(cuò)誤率均是由于水質(zhì)監(jiān)測數(shù)據(jù)的不確定性導(dǎo)致的.最大允許超標(biāo)概率體現(xiàn)了由于監(jiān)測數(shù)據(jù)波動(dòng)性，而對(duì)其采用的“容忍度”或者“安全因子”，決策者對(duì)最大允許超標(biāo)概率的選擇體現(xiàn)了風(fēng)險(xiǎn)偏好;假設(shè)檢驗(yàn)的2類錯(cuò)誤體現(xiàn)了由于水質(zhì)變量的不確定性，使用樣本估計(jì)總體特征的可信程度，決策者對(duì)2類錯(cuò)誤的選擇體現(xiàn)了對(duì)不同群體的利益權(quán)衡.通過這3個(gè)參數(shù)，基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法將水質(zhì)監(jiān)測數(shù)據(jù)的不確定性對(duì)決策過程的影響定量化，對(duì)決策的各種風(fēng)險(xiǎn)進(jìn)行了定量表征，使得決策更具科學(xué)性.具體參見污水寶商城資料或http://wlmqsb.com更多相關(guān)技術(shù)文檔。

　　5 結(jié)論

　　1) 水質(zhì)變量存在不確定性，平均值法忽視了水質(zhì)變量的不確定性，會(huì)造成很大的取偽錯(cuò)誤概率，而超標(biāo)比例法會(huì)導(dǎo)致很大的棄真錯(cuò)誤概率;基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法可以將由水質(zhì)變量的不確定性導(dǎo)致的決策風(fēng)險(xiǎn)定量化，為決策者提供科學(xué)依據(jù).

　　2) 本文在基于二項(xiàng)檢驗(yàn)的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法的基礎(chǔ)上，納入決策者的風(fēng)險(xiǎn)偏好和對(duì)利益的權(quán)衡，提出了水質(zhì)達(dá)標(biāo)評(píng)價(jià)的決策框架.二項(xiàng)分布檢驗(yàn)法對(duì)原始數(shù)據(jù)的分布沒有要求，是一種方便快捷的非參數(shù)方法.本文提出的決策框架，可為決策者提供較大的決策空間，決策者可以根據(jù)風(fēng)險(xiǎn)偏好和利益權(quán)衡選擇或限制最大超標(biāo)率、效應(yīng)值、棄真錯(cuò)誤概率和取偽錯(cuò)誤概率組合，確定最大允許超標(biāo)個(gè)數(shù)，進(jìn)行水質(zhì)達(dá)標(biāo)評(píng)價(jià).

　　3) 對(duì)海河流域的水質(zhì)評(píng)價(jià)結(jié)果表明，采用基于二項(xiàng)分布檢驗(yàn)法的水質(zhì)達(dá)標(biāo)評(píng)價(jià)方法，相對(duì)于平均值法和超標(biāo)比例法，其評(píng)價(jià)結(jié)果具有很高的可信度，不同的利益權(quán)衡會(huì)得到不同的評(píng)價(jià)結(jié)果，驗(yàn)證了二項(xiàng)分布檢驗(yàn)法的適用性和靈活性.