引調(diào)水工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)研究

徐　聰，馮民權(quán)

(西安理工大學(xué) 西北旱區(qū)生態(tài)水利工程國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地，陜西 西安 710048)

?

引調(diào)水工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)研究

徐聰，馮民權(quán)

(西安理工大學(xué) 西北旱區(qū)生態(tài)水利工程國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地，陜西 西安 710048)

【目的】 分析引調(diào)水工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)，為引調(diào)水工程的順利建設(shè)提供科學(xué)依據(jù)?！痉椒ā?利用水質(zhì)綜合指數(shù)評(píng)價(jià)法對(duì)水質(zhì)進(jìn)行評(píng)價(jià)，分別利用水質(zhì)綜合指數(shù)和污染最為嚴(yán)重的指標(biāo)與水量數(shù)據(jù)進(jìn)行Copula水量水質(zhì)聯(lián)合分布函數(shù)的構(gòu)建，并以引沁入汾工程為例對(duì)引調(diào)水工程水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)進(jìn)行研究。【結(jié)果】 引沁入汾工程水源地水質(zhì)狀況良好，總體均能達(dá)到水質(zhì)功能區(qū)Ⅳ類(lèi)水要求；通過(guò)水量與水質(zhì)綜合指數(shù)計(jì)算得到的聯(lián)合風(fēng)險(xiǎn)可以看出：就整個(gè)調(diào)水時(shí)期而言，有利于調(diào)水的概率為99.29%，不利于調(diào)水的概率為2.01%；通過(guò)水量與污染最為嚴(yán)重的水質(zhì)指標(biāo)TN計(jì)算得到的聯(lián)合風(fēng)險(xiǎn)可以看出：就整個(gè)調(diào)水時(shí)期而言，有利于調(diào)水的概率為32.16%，不利于調(diào)水的概率為69.16%?！窘Y(jié)論】 引沁入汾工程水源地適合作為該工程的取水源地，但以后需加強(qiáng)對(duì)TN的控制和治理。

引沁入汾工程；水質(zhì)綜合指數(shù)評(píng)價(jià)法；Copula函數(shù)；水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)

隨著國(guó)民經(jīng)濟(jì)的快速發(fā)展和人民生活質(zhì)量的不斷提高，人們對(duì)水資源的需求和環(huán)境要求越來(lái)越高，一些原有的人類(lèi)生存區(qū)域和空間的水資源已不能滿足本區(qū)域國(guó)民經(jīng)濟(jì)發(fā)展和人類(lèi)生活的需求。由于我國(guó)人均水資源總量不足及水資源時(shí)空分布不均勻，使我國(guó)許多地區(qū)的經(jīng)濟(jì)社會(huì)發(fā)展受制于水資源的供給能力，因而引調(diào)水工程應(yīng)運(yùn)而生。然而，引調(diào)水工程勢(shì)必會(huì)對(duì)水源地的水量水質(zhì)等帶來(lái)一定的風(fēng)險(xiǎn)，對(duì)引調(diào)水工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)的研究可為工程建設(shè)提供科學(xué)依據(jù)。

目前，有關(guān)Copula函數(shù)在引調(diào)水工程方面的研究已有很多，但這些研究多基于徑流豐枯遭遇[1-9]、降水豐枯遭遇[10-13]和洪水遭遇[14-16]層面，少有基于水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)的研究。少數(shù)作者利用Copula函數(shù)僅對(duì)水量水質(zhì)聯(lián)合分布進(jìn)行了研究，但并未對(duì)聯(lián)合風(fēng)險(xiǎn)進(jìn)行具體分析[17]?；诖?，本研究利用Copula函數(shù)從2個(gè)角度對(duì)水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)進(jìn)行研究，在構(gòu)建Copula水量水質(zhì)聯(lián)合分布函數(shù)時(shí)，首先選取多個(gè)水質(zhì)指標(biāo)并引入水質(zhì)綜合評(píng)價(jià)的思想，基于水質(zhì)綜合指數(shù)評(píng)價(jià)法使用簡(jiǎn)便、評(píng)價(jià)結(jié)果直觀及精確度和信息利用率高的優(yōu)點(diǎn)，利用此方法對(duì)水源地水質(zhì)進(jìn)行綜合評(píng)價(jià)，以便較完整地反映水污染程度[18]，并利用水質(zhì)綜合指數(shù)作為水質(zhì)基礎(chǔ)數(shù)據(jù)，進(jìn)行聯(lián)合分布函數(shù)的構(gòu)建，從而進(jìn)行聯(lián)合風(fēng)險(xiǎn)分析；其次選擇水質(zhì)指標(biāo)中污染最為嚴(yán)重的一個(gè)指標(biāo)作為水質(zhì)基礎(chǔ)數(shù)據(jù)，進(jìn)行聯(lián)合分布函數(shù)的構(gòu)建及聯(lián)合風(fēng)險(xiǎn)分析。最后以引沁入汾工程為例進(jìn)行分析，旨在通過(guò)對(duì)取水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)的綜合評(píng)價(jià)，為取水源地的選擇和后期水量水質(zhì)的監(jiān)測(cè)管理提供參考。

1　理論基礎(chǔ)

1.1Copula函數(shù)

Copula函數(shù)是邊緣分布為[0,1]區(qū)間均勻分布的聯(lián)合分布函數(shù)，總體上可以分為橢圓型、阿基米德型和二次型3類(lèi)，其中生成元為1個(gè)參數(shù)的阿基米德型Copula函數(shù)的應(yīng)用最為廣泛，其構(gòu)造形式為：

Cθ(u,v)=φ-1[φ(u)+φ(v)],0≤u,v≤1。

(1)

式中：Cθ(u,v)為Copula函數(shù)，u、v為邊界分布函數(shù)，θ為參數(shù)，φ為Copula函數(shù)的生成元。

由式(1)所示的阿基米德型Copula函數(shù)的定義可知，只要算子的表示形式已知，就可確定一種相應(yīng)的Copula函數(shù)形式。不同算子的選擇會(huì)產(chǎn)生不同類(lèi)別的阿基米德型Copula函數(shù)。幾種常見(jiàn)的二維阿基米德型Copula函數(shù)[19]如表1所示。

表 1　幾種常見(jiàn)的二維阿基米德型Copula函數(shù)Table 1　Several common two-dimensional Archimedean Copula functions

注：Frank Copula函數(shù)τ與θ關(guān)系式中的t(t∈[0,θ])表示變量。

Note：t(t∈[0,θ])represents a variable in the relationshipτandθin Frank Copula function.

表1中的τ是可以描述變量之間非線性相關(guān)關(guān)系的Kendall相關(guān)系數(shù)，可由下式計(jì)算：

(2)

1.2Copula函數(shù)的擬合優(yōu)度檢驗(yàn)

對(duì)一組數(shù)據(jù)系列，有多個(gè)不同Copula函數(shù)可供選擇，其中最優(yōu)的模型需要通過(guò)擬合優(yōu)度檢驗(yàn)進(jìn)行優(yōu)選[7]。Copula函數(shù)的擬合優(yōu)度檢驗(yàn)有多種方法，本研究擬選用以下方法進(jìn)行擬合優(yōu)度檢驗(yàn)。

應(yīng)用雙樣本K-S進(jìn)行檢驗(yàn)，雙樣本K-S檢驗(yàn)統(tǒng)計(jì)量的定義為：

(3)

式中：Ci為聯(lián)合觀測(cè)樣本(xi,yi)的Copula值，ni為聯(lián)合觀測(cè)樣本中滿足x≤xi、y≤yi的聯(lián)合觀測(cè)個(gè)數(shù)，n為樣本個(gè)數(shù)。

為了檢驗(yàn)Copula函數(shù)的擬合優(yōu)度，采用普通最小二乘法(OLS)的最小準(zhǔn)則進(jìn)行評(píng)價(jià)，離差平方和OLS的定義如下：

(4)

式中：Pei為聯(lián)合分布的經(jīng)驗(yàn)頻率，Pi為聯(lián)合分布的理論頻率。

1.3水質(zhì)綜合指數(shù)評(píng)價(jià)法

水質(zhì)綜合指數(shù)評(píng)價(jià)法參考GB 3838-2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》的5級(jí)劃分原則進(jìn)行，定義轉(zhuǎn)化指數(shù)[20]為：

(5)

式中：Si max為標(biāo)準(zhǔn)中i指標(biāo)的最高級(jí)限值；Si min為i指標(biāo)可能的最低值，假定Si min=0。

將所有評(píng)價(jià)項(xiàng)目的各等級(jí)濃度限值做統(tǒng)一的規(guī)范化處理，規(guī)范化公式如下：

(6)

式中：i=1,2,…，n；h=1,2,…，m；sih為第i個(gè)項(xiàng)目的h級(jí)限值；rij為第i個(gè)項(xiàng)目的轉(zhuǎn)化指數(shù)。

計(jì)算各級(jí)評(píng)價(jià)標(biāo)準(zhǔn)的綜合指數(shù)，有：

(7)

確定實(shí)測(cè)樣本的指數(shù)為：

(8)

式中：Xij為i個(gè)指標(biāo)的第j個(gè)實(shí)測(cè)值。

計(jì)算實(shí)測(cè)樣本綜合指數(shù)為：

(9)

2　水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)分析

本研究以引沁入汾工程為背景，對(duì)其水源地水量水質(zhì)的聯(lián)合風(fēng)險(xiǎn)進(jìn)行研究。引沁入汾工程是在安澤縣和川鎮(zhèn)嶺南村興建和川取水輸水工程，開(kāi)鑿分水嶺-草峪嶺隧洞，經(jīng)新建五馬水庫(kù)調(diào)節(jié)，調(diào)水至汾河流域的大型跨流域自流引水工程，不僅可有效緩解汾東地區(qū)嚴(yán)重缺水的供需矛盾，改善當(dāng)?shù)厣鷳B(tài)環(huán)境，還可為汾河下游段實(shí)現(xiàn)清水復(fù)流做出積極貢獻(xiàn)。引沁入汾工程平面示意圖如圖1所示。

圖 1　山西臨汾市引沁入汾工程平面示意圖

2.1基于水質(zhì)綜合指數(shù)的水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)分析

2.1.1基礎(chǔ)數(shù)據(jù)(1)水量數(shù)據(jù)。采用引沁入汾工程水源地飛嶺水文站1980-2010年的年徑流量數(shù)據(jù)，作為構(gòu)建Copula水量水質(zhì)聯(lián)合分布函數(shù)的基礎(chǔ)數(shù)據(jù)。該水文站的實(shí)測(cè)年徑流量時(shí)序如圖2所示。

圖 2　飛嶺水文站實(shí)測(cè)年徑流量時(shí)序圖

(2)水質(zhì)數(shù)據(jù)。采用引沁入汾工程水源地飛嶺水文站1980-2010年水質(zhì)實(shí)測(cè)數(shù)據(jù)中的6個(gè)指標(biāo)(高錳酸鹽指數(shù)及COD、TP、NH3-N、TN與揮發(fā)酚含量)作為水質(zhì)評(píng)價(jià)的基礎(chǔ)數(shù)據(jù)，運(yùn)用水質(zhì)綜合指數(shù)評(píng)價(jià)法對(duì)水質(zhì)進(jìn)行評(píng)價(jià)，得出水質(zhì)綜合指數(shù)值I。然后利用表2的綜合指數(shù)法地表水環(huán)境質(zhì)量等級(jí)劃分標(biāo)準(zhǔn)對(duì)飛嶺站進(jìn)行水質(zhì)綜合指數(shù)評(píng)價(jià)，得到的結(jié)果如表3所示。

表 2　地表水環(huán)境質(zhì)量等級(jí)劃分標(biāo)準(zhǔn)Table 2　Division standard of surface water environmental quality

表 3　基于水質(zhì)綜合指數(shù)評(píng)價(jià)法的沁河飛嶺水文站1980-2010年的水質(zhì)評(píng)價(jià)結(jié)果Table 3　Water quality assessment at Feiling station of Qin River from 1980 to 2010 based on synthetic index method

由于引沁入汾工程的輸水工程壩址位于安澤縣城北的沁河干流上，而沁河為黃河的一級(jí)支流，根據(jù)《山西省地表水域環(huán)境管理區(qū)域方案》，工程所在區(qū)域?yàn)楣I(yè)用水功能，執(zhí)行《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)中Ⅳ類(lèi)水標(biāo)準(zhǔn)。因而，由表2和表3可以看出，引沁入汾工程水源地水質(zhì)狀況良好，總體均能達(dá)到水質(zhì)功能區(qū)Ⅳ類(lèi)水要求。表3中的水質(zhì)綜合指數(shù)I值，即為構(gòu)建聯(lián)合分布函數(shù)所用的水質(zhì)數(shù)據(jù)。

2.1.2邊緣分布Copula函數(shù)方法不受變量邊際分布形式的限制，對(duì)于單變量水文數(shù)據(jù)系列的分析，在我國(guó)常假定水文變量服從皮爾遜Ⅲ型(P-Ⅲ)分布。因此，本部分選用P-Ⅲ型，其分布表達(dá)式為：

(10)

采用矩法對(duì)P-Ⅲ型分布進(jìn)行參數(shù)估計(jì)，得到飛嶺水文站水量和水質(zhì)分布參數(shù)的估計(jì)值如表4所示。

表 4　飛嶺水文站水量和水質(zhì)二維邊緣分布函數(shù)的參數(shù)估計(jì)值Table 4　Parameter estimation of two-dimensional edge distribution function

2.1.3Copula函數(shù)參數(shù)估計(jì)和類(lèi)型選擇計(jì)算所得到的Kendall系數(shù)τ=-0.165 6，呈弱負(fù)相關(guān)關(guān)系。從各Copula函數(shù)對(duì)變量相關(guān)關(guān)系的要求來(lái)看：Frank Copula、Ali-Mikhail-Hap Copula函數(shù)適用于本研究的數(shù)據(jù)序列，可以根據(jù)這2種函數(shù)來(lái)確定θ、D、OLS。各參數(shù)具體計(jì)算結(jié)果如表5所示。

表 5　Copula函數(shù)各參數(shù)計(jì)算結(jié)果Table 5　Parameter calculation of Copula function

由表5可知，F(xiàn)rank Copula函數(shù)的OLS值為0.091 6，Ali-Mikhail-Hap Copula函數(shù)的OLS值為0.101 3，所以，本研究選用Frank Copula函數(shù)作為水量水質(zhì)聯(lián)合分布函數(shù)。

采用以下公式計(jì)算二維聯(lián)合分布的經(jīng)驗(yàn)頻率，即：

H(x，y)=P(X≤xi,Y≤yi)=

(11)

式中：P為滿足X≤xi、Y≤yi的經(jīng)驗(yàn)概率，Nmn為X≤xi且Y≤yi的數(shù)目，N為樣本容量。

圖3為利用Frank Copula函數(shù)模型得出的研究區(qū)域在不同時(shí)間尺度下水量水質(zhì)聯(lián)合分布的理論分布與經(jīng)驗(yàn)點(diǎn)的擬合情況，其復(fù)相關(guān)系數(shù)R2=0.939 9，擬合狀況良好引沁入汾工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)進(jìn)行研究，可將水量水質(zhì)的聯(lián)合分布分為18種情況，其聯(lián)合分布概率的計(jì)算結(jié)果如表6所示。

圖 3　引沁入汾工程水源地水量水質(zhì)聯(lián)合分布理論分布與經(jīng)驗(yàn)點(diǎn)的擬合情況

對(duì)是否有利于調(diào)水，本研究規(guī)定，在水量水質(zhì)聯(lián)合分布的18種情況中，若水源地水質(zhì)達(dá)標(biāo)(即優(yōu)于或等于Ⅳ類(lèi)水)，則有利于調(diào)水；反之，則不利于調(diào)水。由表6可以看出，豐水期、平水期、枯水期水質(zhì)達(dá)標(biāo)的分布概率分別為33.67%，15.46%，50.16%，因水質(zhì)原因，有利于調(diào)水的概率為99.29%。同樣可得，豐、平、枯水期水質(zhì)不達(dá)標(biāo)的概率分別為0.68%，0.31%，1.02%，因水質(zhì)不達(dá)標(biāo)而不利于調(diào)水的概率為2.01%?？傮w而言，引沁入汾工程水源地徑流量無(wú)論處于哪種狀態(tài)，水質(zhì)不達(dá)標(biāo)的概率均較小。在不利于調(diào)水的組合中，處于枯水期且水質(zhì)不達(dá)標(biāo)為最不利調(diào)水情況，其概率為1.02%；在有利于調(diào)水的組合中，處于豐水期且水質(zhì)達(dá)標(biāo)為最有利調(diào)水情況，其概率為33.67%，處于枯水期且水質(zhì)達(dá)標(biāo)為有利調(diào)水最少的情況，其概率為50.16%。

若已知徑流量的豐枯狀態(tài)，可用條件概率的計(jì)算來(lái)預(yù)測(cè)水質(zhì)情況。若已知徑流處于豐水期，則預(yù)測(cè)水質(zhì)超標(biāo)的概率為：

P(Y>62.777 8|X≥1.473 8)=

若徑流處于平水期，則預(yù)測(cè)水質(zhì)超標(biāo)的概率為：

P(Y>62.777 8|1.089 4

若徑流處于枯水期，則預(yù)測(cè)水質(zhì)超標(biāo)的概率為：

P(Y>62.777 8|X≤1.089 4)=

從條件概率的計(jì)算可以看出：無(wú)論徑流量處于哪種狀態(tài)，所預(yù)測(cè)水質(zhì)的超標(biāo)概率均較小。

2.2基于TN的水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)分析

采用引沁入汾工程水源地飛嶺水文站1980-2010年的年徑流量數(shù)據(jù)和污染最為嚴(yán)重的水質(zhì)指標(biāo)TN作為構(gòu)建Copula水量水質(zhì)聯(lián)合分布函數(shù)的基礎(chǔ)數(shù)據(jù)。經(jīng)過(guò)對(duì)比分析，水量采用P-Ⅲ型分布，水質(zhì)指標(biāo)TN采用Gumbel分布。經(jīng)計(jì)算Frank Copula函數(shù)適合作為該聯(lián)合分布的函數(shù)。

根據(jù)所建立的模型對(duì)引沁入汾工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)進(jìn)行研究，可將水量與TN的聯(lián)合分布分為18種情況，具體情況及計(jì)算結(jié)果如表7所示。在表7水量水質(zhì)聯(lián)合分布的18種情況中，若水源地TN達(dá)標(biāo)(即優(yōu)于或等于Ⅳ類(lèi)水中TN的標(biāo)準(zhǔn))，則有利于調(diào)水；反之，則不利于調(diào)水。經(jīng)過(guò)分析可以得出：就整個(gè)調(diào)水時(shí)期而言，有利于調(diào)水的概率為32.16%，不利于調(diào)水的概率為69.16%。在不利于調(diào)水的組合中，處于枯水期且TN不達(dá)標(biāo)為最不利調(diào)水情況，其概率為34.96%；在有利于調(diào)水的組合中，處于豐水期且TN達(dá)標(biāo)為最有利調(diào)水情況，其概率為10.91%，處于枯水期且TN達(dá)標(biāo)為有利調(diào)水最少的情況，其概率為16.23%。

表 7　基于TN的引沁入汾工程水源地水量水質(zhì)風(fēng)險(xiǎn)的聯(lián)合分布概率Table 7　Joint distribution probability of water quantity and quality risk of water sources of water diversion project from Qin River to Fen River based on TN 　%

3　結(jié)　論

(1)引沁入汾工程水源地水質(zhì)狀況良好，總體均能達(dá)到水質(zhì)功能區(qū)Ⅳ類(lèi)水要求。

(2)通過(guò)水量與水質(zhì)綜合指數(shù)計(jì)算得到的聯(lián)合風(fēng)險(xiǎn)可以看出：就整個(gè)調(diào)水時(shí)期而言，有利于調(diào)水的概率為99.29%，不利于調(diào)水的概率為2.01%。在不利于調(diào)水的組合中，處于枯水期且水質(zhì)不達(dá)標(biāo)為最不利調(diào)水情況，其概率為1.02%；在有利于調(diào)水的組合中，處于豐水期且水質(zhì)達(dá)標(biāo)為最有利調(diào)水情況，其概率為33.67%，處于枯水期且水質(zhì)達(dá)標(biāo)為有利調(diào)水量最少的情況，其概率為50.16%。

(3)通過(guò)水量與污染最為嚴(yán)重的水質(zhì)指標(biāo)TN計(jì)算得到的聯(lián)合風(fēng)險(xiǎn)可以看出：就整個(gè)調(diào)水時(shí)期而言，有利于調(diào)水的概率為32.16%，不利于調(diào)水的概率為69.16%。在不利于調(diào)水的組合中，處于枯水期且TN不達(dá)標(biāo)為最不利調(diào)水情況，其概率為34.96%；在有利于調(diào)水的組合中，處于豐水期且TN達(dá)標(biāo)為最有利調(diào)水情況，其概率為10.91%，處于枯水期且TN達(dá)標(biāo)為有利調(diào)水最少的情況，其概率為16.23%。

由此可以看出：引沁入汾工程水源地適合作為該取水工程的取水源地，但以后需加強(qiáng)對(duì)TN的控制和治理，該研究結(jié)果既可為引調(diào)水工程水源地水量水質(zhì)聯(lián)合風(fēng)險(xiǎn)的研究提供借鑒，也可為引調(diào)水工程的建設(shè)提供科學(xué)依據(jù)。

[1]莫淑紅,沈冰,張曉偉,等.基于Copula函數(shù)的河川徑流豐枯遭遇分析 [J].西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,37(6):131-136.

Mo S H,Shen B,Zhang X W,et al.Synchronous-asynchronous encounter probability of rich-poor runoff based on Copula function [J].Journal of Northwest A&F University(Nat Sci Ed),2009,37(6):131-136.

[2]丁志宏,馮平,張永.基于Copula模型的豐枯頻率分析：以南水北調(diào)西線工程調(diào)水區(qū)徑流與黃河上游來(lái)水的豐枯遭遇研究為例 [J].長(zhǎng)江流域資源與環(huán)境,2010,19(7):759-764.

Ding Z H,Feng P,Zhang Y.Rich poor frequency analysis of runoff based on Copula model：a case study on encounter probability of rich-poor runoff between rivers in area of western route of south to north water transfer project and that of the upper reach of the Yellow river [J].Resources and Environment in the Yangtze Basin,2010,19(7):759-764.

[3]張忠波,張雙虎,王浩.基于Copula函數(shù)的三峽工程供水期豐枯遭遇分析 [J].人民長(zhǎng)江,2012,43(3):5-8.

Zhang Z B,Zhang S H,Wang H,Analysis of wetness:dryness encounter in TGP water supply period based on Copula function [J].Yangtze River,2012,43(3):5-8.

[4]謝華,羅強(qiáng),黃介生. 基于三維copula函數(shù)的多水文區(qū)豐枯遭遇分析 [J].水科學(xué)進(jìn)展,2012,23(2):186-193.

Xie H,Luo Q,Huang J S.Synchronous asynchronous encounter analysis of multiple hydrologic regions based on 3D copula function [J].Advances in Water Science,2012,23(2):186-193.

[5]何慶升,朱永忠.應(yīng)用高維Copula函數(shù)的月徑流量聯(lián)合概率分布研究 [J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,27(8):112-120.

He Q S,Zhu Y Z.Studying of annual flow joint probability distribution based on high dimensions Copula function [J].Journal of Chongqing University of Technology(Natural Science Edition),2013,27(8):112-120.

[6]Jeong C,Lee T.Copula-based modeling and stochastic simulation of seasonal intermittent streamflows for arid regions [J].Journal of Hydro-environment Research,2015,9(4):604-613.

[7]馮民權(quán),郝竹林,張園園,等.禹門(mén)口調(diào)水工程水源區(qū)與受水區(qū)徑流豐枯遭遇風(fēng)險(xiǎn)分析研究 [J].自然災(zāi)害學(xué)報(bào),2012,21(6):156-163.

Feng M Q,Hao Z L,Zhang Y Y,et al.Risk analysis of runoff rich poor encounter between water source area and receiving area of Yumenkou water transfer project [J].Journal of Natural Disasters,2012,21(6):156-163.

[8]Yu K X,Xiong L H,Lars G.Derivation of low flow distribution functions using copulas [J].Journal of Hydrology,2014,508:273-288.

[9]牛軍宜,馮平,丁志宏.基于多元Copula函數(shù)的引灤水庫(kù)徑流豐枯補(bǔ)償特性研究 [J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2009,39(6):1095-1100.

Niu J Y,Feng P,Ding Z H.Study of the wetness-dryness compensation characteristics of Luan River diversion reservoirs’ inflow based on multivariate Copula Functions [J].Journal of Jilin University(Earth Science Edition),2009,39(6):1095-1100.

[10]鄭紅星,劉昌明.南水北調(diào)東中兩線不同水文區(qū)降水豐枯遭遇性分析 [J].地理學(xué)報(bào),2000,55(5):523-532.

Zheng H X,Liu C M.Analysis on a synchronism-synchronism of regional precipitation in planned South-to-North water transfer areas [J].Acta Geographica Sinica,2000,55(5):523-532.

[11]熊其玲,何小聰,康玲.基于Copula函數(shù)的南水北調(diào)中線降水豐枯遭遇分析 [J].水電能源科學(xué),2009,27(6):9-11，112.

Xiong Q L,He X C,Kang L.Synchronous asynchronous encounter probability of rich poor precipitation between water source area and water receiving area based on Copula theory in middle route of South to North water transfer project [J].Water Resources and Power,2009,27(6):9-11，112.

[12]王偉,鐘永華,雷曉輝,等.引漢濟(jì)渭工程水源區(qū)與受水區(qū)豐枯遭遇分析 [J].南水北調(diào)與水利科技,2012,10(5):23-26，36.

Wang W,Zhong Y H,Lei X H,et al.Synchronous-asynchronous encounter probability of rich-poor precipitation between water source area and water receiving area of the Hanjiang-to-Weihe River water transfer project [J].South-to-North Water Transfers and Water Science & Technology,2012,10(5):23-26，36.

[13]Zoubeida K B,Andras B.Modeling short duration extreme pr-ecipitation patterns using copula and generalized maximum pseudo-likelihood estimation with censoring [J].Advances in Water Resources,2015,84:1-13.

[14]陳璐,郭生練,張洪剛,等.長(zhǎng)江上游干支流洪水遭遇分析 [J].水科學(xué)進(jìn)展,2011,22(3):323-330.

Chen L,Guo S L,Zhang H G,et al.Flood coincidence probability analysis for the upstream Yangtze River and its tributaries [J].Advances in Water Science,2011,22(3):323-330.

[15]李偉.Copula函數(shù)在多變量洪水聯(lián)合分布中的應(yīng)用研究 [D].武漢：華中科技大學(xué),2013.

Li W.The application and research of multivate flood joint distribution based on Copula theory [D].Wuhan:Huazhong University of Science and Technology,2013.

[16]Fan Y R,Huang W W,Huang G H,et al.Hydrologic risk analysis in the Yangtze River basin through coupling Gaussian mixtures into copulas [J].Advances in Water Resources,2016,88:170-185.

[17]張翔,冉啟香,夏軍,等.基于Copula函數(shù)的水量水質(zhì)聯(lián)合分布函數(shù) [J].水利學(xué)報(bào),2011,42(4):483-489.

Zhang X,Ran Q X,Xia J,et al.Jointed distribution function of water quality and water quantity based on Copula [J].Journal of Hydraulic Engineering,2011,42(4):483-489.

[18]陳仁杰,錢(qián)海雷,闞海東，等.水質(zhì)評(píng)價(jià)綜合指數(shù)法的研究進(jìn)展 [J].環(huán)境與職業(yè)醫(yī)學(xué),2009,26(6):581-584.

Chen R J,Qian H L,Kan H D,et al.A review of comprehensive index methods in water quality assessment [J].Journal of Environmental & Occupational Medicine,2009,26(6):581-584.

[19]Zhang L.Multivariate hydrological frequency analysis and risk mapping [D].Baton Rouge:Louisiana State University,2005.

[20]李凡修,陳武,梅平，等.淺層地下水環(huán)境質(zhì)量評(píng)價(jià)的綜合指數(shù)模型 [J].地下水,2004,26(1):36-37.

Li F X,Chen W,Mei P,et al.Synthetic index model of shallow groundwater environment quality evaluation [J].Ground Water,2004,26(1):36-37.

Joint risk of water quantity and quality in water sources of water diversion project

XU Cong,FENG Minquan

(StateKeyLaboratoryBaseofEco-hydraulicEngineeringinAridArea,Xi’anUniversityofTechnology,Xi’an，Shaanxi710048,China)

【Objective】 This paper studied the joint risk of water quantity and quality in water sources of water diversion project to provide scientific basis for successful construction of water diversion project.【Method】 Water quality was evaluated by synthetic index method.Water quantity data was used together with water quality synthetic index and the most serious pollution indicator to build the Copula joint distribution function.Then,the joint risk of water quantity and quality was investigated using the water diversion project from Qin River to Fen River as an example.【Result】 The water sources of water diversion project from Qin River to Fen River have good water quality,and the overall water quality reaches ribbon Class Ⅳ water requirements.From the joint risk calculated by water quality and quantity synthetic index,the probability that is conducive to diversion is 99.29%,and the probability that is not conducive to diversion is 2.01% during the entire water diversion period.From the joint risk of water contamination by the most serious water quality indicator TN,the probability that is conducive to diversion is 32.16% and the probability that is not conducive to diversion is 69.16%.【Conclusion】 The water sources of water diversion project from Qin River to Fen River is suitable and the control and governance of TN should be strengthened in future.

water diversion project from Qin River to Fen River;synthetic index method for water quality assessment;Copula function;joint risk of water quantity and quality

時(shí)間:2016-08-0909:41DOI：10.13207/j.cnki.jnwafu.2016.09.031

2015-02-06

高等學(xué)校博士學(xué)科點(diǎn)專(zhuān)項(xiàng)科研基金項(xiàng)目(20126118110015)；中央財(cái)政支持地方高校發(fā)展專(zhuān)項(xiàng)資金特色重點(diǎn)學(xué)科項(xiàng)目(106-00X101，106-5X1205)；山西省中部引黃工程管理局計(jì)劃基金項(xiàng)目(2014-9)

徐聰(1990-)，男，河南長(zhǎng)葛人，在讀碩士，主要從事水環(huán)境模擬與預(yù)測(cè)研究。E-mail:1064974686@qq.com

馮民權(quán)(1964-)，男，山西永濟(jì)人，教授，博士，博士生導(dǎo)師，主要從事水環(huán)境模擬與污染控制研究。

E-mail:mqfeng@xaut.edu.cn

X824.04

A

1671-9387(2016)09-0228-07

網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20160809.0941.062.html