現(xiàn)代水網(wǎng)建設(shè)的區(qū)域生態(tài)效應(yīng)可變模糊評價(jià)

高峰 齊真 王好芳 趙然杭 王興菊 白珊

摘要：現(xiàn)代水網(wǎng)建設(shè)旨在解決區(qū)域水資源短缺、地下水超采以及生態(tài)惡化等問題，為了評估現(xiàn)代水網(wǎng)建設(shè)的生態(tài)功能，對現(xiàn)代水網(wǎng)建設(shè)產(chǎn)生的生態(tài)效應(yīng)進(jìn)行評價(jià)。以青州市為研究區(qū)域，以2013年為基準(zhǔn)年，2020年和2025年為規(guī)劃年，在50%、75%和95%保證率下，通過青州市現(xiàn)代水網(wǎng)進(jìn)行水資源優(yōu)化配置和調(diào)度，根據(jù)配置和調(diào)度結(jié)果，選取5個(gè)生態(tài)效應(yīng)評價(jià)指標(biāo)，構(gòu)建生態(tài)效應(yīng)評價(jià)指標(biāo)體系，采用可變模糊評價(jià)法，對青州市現(xiàn)代水網(wǎng)建設(shè)產(chǎn)生的生態(tài)效應(yīng)進(jìn)行評價(jià)。結(jié)果表明，青州市現(xiàn)代水網(wǎng)建設(shè)后，生態(tài)狀況等級比建設(shè)前提升了1～3級，生態(tài)環(huán)境明顯改善。因此，區(qū)域現(xiàn)代水網(wǎng)建設(shè)和生態(tài)調(diào)度對生態(tài)環(huán)境改善有重要作用。

關(guān)鍵詞：生態(tài)效應(yīng)；可變模糊評價(jià)；現(xiàn)代水網(wǎng)建設(shè)；水資源；青州市

中圖分類號：TV213.9文獻(xiàn)標(biāo)志碼：A文章編號：16721683（2018）03010208

Variable fuzzy evaluation of regional ecological effects of modern water network construction

GAO Feng1，QI Zhen2，WANG Haofang2，ZHAO Ranhang2，WANG Xingju2，BAI Shan2

（1.Water Resources Research Institute of Shandong Province，Jinan 250013，China；

2.School of Civil Engineering，Shandong University，Jinan 250061，China）

Abstract：The modern water network construction is intended to solve the problems of water shortage，groundwater overdraft，and ecological deterioration.This paper evaluates the ecological effects of modern water network construction so as to assess the ecological function of modern water network construction.We took Qingzhou City as the research area，and took 2013 as the base year and 2020 and 2025 as the planning years.At 50%，75% and 95% assurance rates，we conducted optimal allocation and dispatching of water resources through Qingzhou modern water network.According to the allocation and dispatching results，we selected 5 evaluation indexes of ecological effects，and constructed an index system of ecological effects evaluation.Using variable fuzzy evaluation method，we evaluated the ecological effects of the construction of modern water network in Qingzhou.The results showed that the ecological level improved by 13 grades after the modern water network was constructed，and the ecological environment was obviously improved.Therefore，the regional modern water network construction and ecological dispatching play an important role in improving the ecological environment.

Key words：ecological effect；variable fuzzy evaluation；modern water network construction；water resources；Qingzhou City

生態(tài)效應(yīng)是指人類活動(dòng)引發(fā)的生態(tài)系統(tǒng)的變化和響應(yīng)，按照性質(zhì)可以分為正效應(yīng)和負(fù)效應(yīng)[1]。水利工程生態(tài)效應(yīng)是指人類通過修建各類水利工程對生態(tài)系統(tǒng)造成的結(jié)構(gòu)和功能的影響，包括生態(tài)破壞及生態(tài)修復(fù)兩種效果[23]。針對缺水地區(qū)，進(jìn)行區(qū)域水網(wǎng)建設(shè)是生態(tài)修復(fù)的有效途徑。20世紀(jì)以來，美國、前蘇聯(lián)、澳大利亞及中國等國家，通過跨流域或區(qū)域調(diào)水，對水資源進(jìn)行重新分配，以緩解區(qū)域水資源短缺等問題，促進(jìn)社會(huì)經(jīng)濟(jì)的發(fā)展[4]。近年來，隨著生態(tài)問題的日益嚴(yán)重，水網(wǎng)建設(shè)等調(diào)水工程產(chǎn)生的生態(tài)效應(yīng)逐漸受到關(guān)注[5]。

生態(tài)效應(yīng)具有一定的模糊性，一般是通過構(gòu)建指標(biāo)體系，采用相應(yīng)評價(jià)方法進(jìn)行評價(jià)。目前，常用的評價(jià)方法主要有模糊綜合評價(jià)法[6]、模糊模式識(shí)別法[7]、可變模糊評價(jià)法[813]、集對評價(jià)法[14]等，其中，可變模糊評價(jià)法能夠科學(xué)合理地確定各樣本指標(biāo)的相對隸屬度和相對隸屬函數(shù)，進(jìn)而確定各樣本的級別特征值。在水利行業(yè)，該方法廣泛應(yīng)用于河流健康評價(jià)、生態(tài)安全評價(jià)和風(fēng)險(xiǎn)評估等領(lǐng)域[1519]。因此，本文以青州市為研究區(qū)域，以2013年為基準(zhǔn)年，2020年和2025年為規(guī)劃年，在50%、75%和95%保證率下，通過研究區(qū)域現(xiàn)代水網(wǎng)的水資源優(yōu)化配置和調(diào)度，應(yīng)用可變模糊集理論[89]和陳守煜提出的可變模糊評價(jià)法[1013]對基于現(xiàn)代水網(wǎng)建設(shè)的配置和調(diào)度結(jié)果，進(jìn)行區(qū)域生態(tài)效應(yīng)評價(jià)，為區(qū)域水網(wǎng)的運(yùn)行調(diào)度提供科學(xué)依據(jù)。

第16卷 總第96期·南水北調(diào)與水利科技·2018年6月高峰等·現(xiàn)代水網(wǎng)建設(shè)的區(qū)域生態(tài)效應(yīng)可變模糊評價(jià)1研究方法與研究區(qū)域

1.1研究方法

設(shè)U為論域，u為U中的任意元素。設(shè)區(qū)間[a，b]為實(shí)軸上模糊可變集合V的吸引域區(qū)間，相對差異度DA（u）∈（0，1]；區(qū)間[c，d]為實(shí)軸上包含[a，b]的V的范圍域區(qū)間，相對差異度DA（u）∈[-1，1]；設(shè)M為[a，b]中相對差異度DA（u）=1的點(diǎn)值；設(shè)x為[c，d]中的任意點(diǎn)值，如圖1所示[10]。

建立評價(jià)指標(biāo)特征值矩陣X=（xij）、指標(biāo)標(biāo)準(zhǔn)值矩陣Y=（yih）以及點(diǎn)值Mih的矩陣。其中，i為指標(biāo)；j為樣本；h為級別，h=1，2，…，k，k為級別數(shù) 。

根據(jù)標(biāo)準(zhǔn)值矩陣Y，建立吸引域矩陣Iab和范圍域矩陣Icd，判斷樣本特征值xij的位置。根據(jù)式（1）計(jì)算相對差異度矩陣[DA（xij）h]，根據(jù)式（2）計(jì)算相對隸屬度矩陣[μA（xij）h]。

DA（u）=x-aM-aβ；x∈[a，M]

DA（u）=-x-ac-aβ；x∈[c，a]

DA（u）=x-bM-bβ；x∈[M，b]

DA（u）=-x-bd-aβ；x∈[b，d]

DA（u）=-1；x

瘙 綋 （c，d）（1）

μA（u）=（1+DA（u））/2（2）

式中：β為非負(fù)指數(shù)，通?？扇ˇ?1。

根據(jù)指標(biāo)重要性排序一致性定理[11]，確定評價(jià)指標(biāo)的權(quán)向量ωi。

應(yīng)用公式

ju′h=11+∑mi=1[ωi（1-μA（xij）h）]p∑mi=1（ωiμA（xij）h）pαp（3）

式中：ju′h為非歸一化的綜合相對隸屬度；ωi為指標(biāo)權(quán)重；α為模型優(yōu)化準(zhǔn)則參數(shù)，取1或2；m為識(shí)別指標(biāo)數(shù)；p為距離參數(shù)，取1或2，其中，p=1為海明距離，p=2為歐氏距離。

在α和p的不同組合情況下，分別求得樣本的相對隸屬度，并通過歸一化公式

juh=ju′h∑ch=1ju′h（4）

得到樣本級別特征值矩陣[jUh]。式中，juh為歸一化的綜合相對隸屬度。

應(yīng)用公式

H=（1，2，…，k）·U（5）

計(jì)算樣本級別特征值H。

1.2研究區(qū)域概況

青州市地處山東半島中部，水資源短缺，人均水資源量僅為全國平均水平的1/6。受氣候和地理位置影響，青州市河川徑流量分布不均。西南部河川徑流量占全市的85.4%；東北北部平原區(qū)僅占全市的14.6%。而青州市東北部又是工業(yè)區(qū)，水資源短缺嚴(yán)重，出現(xiàn)了工農(nóng)業(yè)用水大量擠占生態(tài)環(huán)境用水，造成了東北部地下水嚴(yán)重超采，形成了漏斗區(qū)，對青州市的生態(tài)環(huán)境造成了嚴(yán)重影響[20]。

為此，青州市于2012年開始進(jìn)行了現(xiàn)代水網(wǎng)建設(shè)，見圖2，以青州市“五庫”（即仁河水庫、黑虎山水庫、七一水庫、壩溝子水庫和趙家莊水庫）和“六河”（即淄河、仁河、大石河、南陽河、彌河和北陽河）為骨架，使庫河聯(lián)通、庫庫聯(lián)通以及河河聯(lián)通，實(shí)現(xiàn)青州市水系貫通，豐枯互補(bǔ)。通過水資源合理配置和調(diào)度，將西南部豐富的水資源輸送到東北部平原區(qū)，以解決青州市東北部缺水及地下水漏斗區(qū)問題[2122]。

2可變模糊評價(jià)

本文以2013年為基準(zhǔn)年，2020年和2025年為規(guī)劃年，在50%、75%和95%保證率下，通過青州市現(xiàn)代水網(wǎng)建設(shè)的水資源優(yōu)化配置和調(diào)度方案，對青州市現(xiàn)代水網(wǎng)建設(shè)產(chǎn)生的生態(tài)效應(yīng)進(jìn)行評價(jià)。

2.1水資源優(yōu)化調(diào)度方案

基于青州市現(xiàn)代水網(wǎng)建設(shè)，以2013年為基準(zhǔn)年，2020年和2025年為規(guī)劃年，確定50%、75%和95%保證率下的各水源可供水量及各用水戶需水量，以社會(huì)、經(jīng)濟(jì)和生態(tài)綜合效益最大的為目標(biāo)，建立多目標(biāo)的水資源優(yōu)化配置模型，見式（6）和式（7）。

目標(biāo)函數(shù)：minf1=∑Kk=1（Dk-Qk）

maxf2=∑Kk=1（bkQk-ckQk）

minf3=（De-QeDe2（6）

約束條件：Qkmin≤Qk≤Dk

∑Kk=1Qk≤Q

Qk≥0（7）

式中：f1為社會(huì)效益；Dk為第k個(gè)用戶的需水量（萬m3）；Qk為第k個(gè)用水戶的供水量（萬m3）；K為用水戶數(shù)量；f2為經(jīng)濟(jì)效益；bk為第k個(gè)用水戶的供水效益系數(shù)（元/ m3），ck為第k個(gè)用水戶的供水費(fèi)用系數(shù)（元/ m3）；f3為生態(tài)效益；De生態(tài)需水量；Qe為生態(tài)供水量（萬m3）；Qkmin為第k個(gè)用水戶的供水下限，Q為可供水量。

利用理想點(diǎn)法求解該模型，得到青州市水資源優(yōu)化配置方案，然后，在近水近用、優(yōu)水優(yōu)用等原則下，既考慮水網(wǎng)工程的供水能力和用水實(shí)際需求，制定青州市水資源優(yōu)化調(diào)度方案，其中，青州市水資源優(yōu)化調(diào)度方案中生態(tài)部分詳見表1和表2。

2.2指標(biāo)體系構(gòu)建

為了有效地評價(jià)現(xiàn)代水網(wǎng)建設(shè)產(chǎn)生的生態(tài)效應(yīng)，根據(jù)實(shí)地考察發(fā)現(xiàn)，青州市生態(tài)系統(tǒng)脆弱或惡化是由于水資源短缺，可供水量不能有效地滿足生態(tài)用水需求造成的。因此，本文結(jié)合青州市生態(tài)狀況，以科學(xué)性、代表性、簡明和可操作性的原則，重點(diǎn)從水量方面考慮河道內(nèi)和河道外生態(tài)用水狀況，選取生態(tài)效應(yīng)指標(biāo)，構(gòu)建生態(tài)效應(yīng)評價(jià)指標(biāo)體系，以評價(jià)青州市現(xiàn)代水網(wǎng)建設(shè)對青州市東北部生態(tài)狀況的改善情況。各指標(biāo)選取如下。

（1）生態(tài)水資源滿足程度（%）。

生態(tài)水資源滿足程度指河道內(nèi)外總的生態(tài)供水量與需水量的比值（以百分比表示）。該指標(biāo)反映青州市東北部平原區(qū)生態(tài)用水的整體滿足程度。

（2）地下水開采率（%）。

地下水開采率是指地下水開采量與可開采量的比值（以百分比表示）。可開采量指每年的地下水補(bǔ)給量。當(dāng)?shù)叵滤_采量大于補(bǔ)給量時(shí)，就會(huì)造成地下水超采，超采嚴(yán)重時(shí)會(huì)導(dǎo)致地下漏斗區(qū)的產(chǎn)生。因此，地下水開采率能夠反映青州市東北部平原區(qū)地下水開采對地下水漏斗區(qū)惡化或緩解程度。

（3）河道內(nèi)生態(tài)水資源滿足程度（%）。

河道內(nèi)生態(tài)水資源滿足程度是指河道內(nèi)生態(tài)供水流量與生態(tài)需水流量的比值（以百分比表示）。青州市東北部平原區(qū)最主要的兩條河流為北陽河和南陽河。因此，河道內(nèi)生態(tài)水資源滿足程度主要是指北陽河和南陽河河道內(nèi)的生態(tài)用水滿足情況，反映北陽河和南陽河河道內(nèi)的生態(tài)狀況。

（4）河道外生態(tài)水資源滿足程度（%）。

河道外生態(tài)水資源滿足程度是指河道外生態(tài)供水量與需水量的比值（以百分比表示）。主要指青州市東北部綠化用水、澆灑路面用水等滿足情況。

（5）生態(tài)供水保證率（%）。

生態(tài)供水保證率是指在用來供水的河道內(nèi)，流量大于基本生態(tài)流量的天數(shù)比例。彌河是青州市最大的過境河流，水資源較為豐富，是重要的客水資源，從彌河調(diào)水可有效緩解青州市東北部的缺水狀態(tài)，對地下漏斗區(qū)的修復(fù)具有重要作用。但從彌河調(diào)水時(shí)應(yīng)首先要保證彌河河道內(nèi)的基本生態(tài)用水，因此，選用該指標(biāo)來度量彌河可向青州市東北部調(diào)度水資源的能力。

根據(jù)青州市2013年水資源開發(fā)利用情況、現(xiàn)代水網(wǎng)建設(shè)的水資源優(yōu)化調(diào)度方案（表1和表2）以及各水平年彌河水資源狀況，得到各評價(jià)指標(biāo)特征值，詳見表3。

評價(jià)等級劃分

評價(jià)等級劃分是依據(jù)研究區(qū)的實(shí)際生態(tài)狀況，結(jié)合水資源規(guī)劃等相關(guān)資料，統(tǒng)計(jì)各指標(biāo)的標(biāo)準(zhǔn)差、平均值、中位數(shù)、最大值和最小值等統(tǒng)計(jì)參數(shù)特征值，并參考相關(guān)文獻(xiàn)[2325]，將青州市各評價(jià)指標(biāo)分為5個(gè)等級，等級1情況屬優(yōu)，等級2情況屬良，等級3情況屬中等，等級4情況屬較差，等級5情況屬很差。各評價(jià)指標(biāo)劃分標(biāo)準(zhǔn)詳見表4。

根據(jù)表1和表2建立青州市評價(jià)指標(biāo)特征值矩陣X=（xij）與指標(biāo)標(biāo)準(zhǔn)值矩陣Y=（yih）。

X=312416

127199338

20125

656050

758350=（xij）

Y=

90～10080～9060～8030～600～30

0～5050～7575～100100～200200～500

90～10080～9060～8030～600～30

90～10080～9060～8030～600～30

95～10085～9570～8550～700～50=（yih）

其中，i=1，2，…，5為指標(biāo)號；j=1，2，3為水平年號；h=1，2，…，5為級別編號。

參照指標(biāo)標(biāo)準(zhǔn)值矩陣Y，建立青州市生態(tài)狀況評價(jià)指標(biāo)可變集合的吸引域矩陣Iab與范圍域矩陣Icd以及點(diǎn)值Mih的矩陣。Iab=[90，100][80，90][60，80][30，60][0，30]

[0，50][50，75][75，100][100，200][200，500]

[90，100][80，90][60，80][30，60][0，30]

[90，100][80，90][60，80][30，60][0，30]

[95，100][85，95][70，85][50，70][0，50]

Icd=[80，100][60，100][30，90][0，80][0，60]

[0，75][0，100][50，200][75，500][100，500]

[80，100][60，100][30，90][0，80][0，60]

[80，100][60，100][30，90][0，80][0，60]

[85，100][70，100][50，95][0，85][0，70]

M=1008570450

06590150500

1008570450

1008570450

1009080600

根據(jù)公式（1）計(jì)算相對差異度矩陣[DA（xij）h]。

[DA（xi1）h]=-1.00-1.00-0.970.07-0.03

-1.00-1.00-0.270.54-0.73

-1.00-1.00-1.000.330.33

-1.00-0.75-0.500.25-1.00

-1.00-0.67-0.500.33-1.00

[DA（xi2）h]=-1.00-1.00-0.11-0.200.20

-1.00-1.00-0.990.02-0.01

-1.00-1.00-1.00-0.600.60

-1.00-1.000.000.00-1.00

-1.00-0.130.40-0.87-1.00

[DA（xi3）h]=-1.00-1.00-1.00-0.470.47

-1.00-1.00-1.00-0.460.46

-1.00-1.00-1.00-0.830.83

-1.00-1.00-0.330.67-0.67

-1.00-1.00-1.000.000.00

根據(jù)公式（2）計(jì)算指標(biāo)j對h級的相對隸屬度矩陣[μA（xij）h]。

[μA（xi1）h]=0.000.000.020.540.49

0.000.000.370.770.14

0.000.000.000.340.67

0.000.130.750.380.00

0.000.170.750.340.00

[μA（xi2）h]=0.000.000.000.400.60

0.000.000.010.510.50

0.000.000.000.200.80

0.000.000.500.500.00

0.000.440.700.070.00

[μA（xi3）h]=0.000.000.000.270.74

0.000.000.000.270.73

0.000.000.000.090.92

0.000.000.340.840.17

0.000.000.000.500.50

根據(jù)指標(biāo)重要性排序一致性定理[11]，得到指標(biāo)重要性排序一致性標(biāo)度矩陣F。

排序

F=0.51111

00.5111

000.511

0000.51

00000.5①

②

③

④

⑤

按矩陣F關(guān)于重要性的排序，以排序?yàn)棰俚闹笜?biāo)x1逐一地與排序?yàn)棰?、③、④、⑤的指?biāo)，進(jìn)行二元比較。指標(biāo)x1與指標(biāo)x2相比，處于“略為”；指標(biāo)x1與指標(biāo)x3相比，處于“明顯”；指標(biāo)x1與指標(biāo)x3相比，處于“十分”；指標(biāo)x1與指標(biāo)x5相比，處于“非?！?，根據(jù)表5得到非歸一化的評價(jià)指標(biāo)權(quán)向量為

ω′1=（1.000，0.667，0.429，0.250，0.176）

分析與討論

表4和表5結(jié)果表明，不同權(quán)重下，各樣本評價(jià)等級一致，可作為最終結(jié)果；青州市現(xiàn)代水網(wǎng)建設(shè)后，相較于2013年，2020年生態(tài)狀況在50%、75%和95%保證率下分別由4級、4～5級和5級提高到2級、2級和3～4級，2025年生態(tài)狀況分別提高到1～2級、1～2級和4級。

其他的評價(jià)方法通常選擇一套權(quán)重以及建立一套評價(jià)模型，往往不能保證結(jié)果的可靠性和準(zhǔn)確性。本文采用可變模糊評價(jià)法，選取兩套權(quán)重，并在基本模型的基礎(chǔ)上通過變換參數(shù)來變化模型（通過變換參數(shù)α和p，將模型變化為線性和非線性），將穩(wěn)定結(jié)果作為最終結(jié)果，對青州市現(xiàn)代水網(wǎng)建設(shè)的生態(tài)效應(yīng)進(jìn)行評價(jià)，提高了結(jié)果的可靠性和可信度。

3結(jié)論

以青州市為研究區(qū)域，進(jìn)行現(xiàn)代水網(wǎng)建設(shè)的區(qū)域生態(tài)效應(yīng)評價(jià)。結(jié)果表明，青州市現(xiàn)代水網(wǎng)建設(shè)后，在不同保證率下，生態(tài)狀況評價(jià)等級均有所提高，生態(tài)環(huán)境得到明顯改善，生態(tài)效應(yīng)顯著?？勺兡：u價(jià)法能夠客觀確定現(xiàn)代水網(wǎng)建設(shè)前后的生態(tài)狀況等級，結(jié)果較為可靠準(zhǔn)確，為定量評價(jià)生態(tài)效應(yīng)提供了參考依據(jù)，對區(qū)域水生態(tài)文明建設(shè)具有積極意義。

隨著水網(wǎng)建成，在實(shí)際運(yùn)行管理中，建議增加水質(zhì)和生物多樣性等資料，使該指標(biāo)體系更加完善，評價(jià)結(jié)果更加科學(xué)合理。

參考文獻(xiàn)（References）：

[1]馬鐵民，曾容，徐菲.遼河流域典型水庫生態(tài)效應(yīng)評價(jià)[J].東北水利水電，2008（10）：3740.（MA T M，ZENG R，XU F.Ecological effect evaluation of typical reservoir in Liao river basin[J].Water Resources & Hydropower of Northeast，2008（10）：3740.（in Chinese）） DOI：10.3969/j.issn.10020624.2008.10.015.

[2]常本春，耿雷華，劉翠善，等.水利水電工程的生態(tài)效應(yīng)評價(jià)指標(biāo)體系[J].水利水電科技進(jìn)展，2006（6）：1115.（CHANG B C，GENG L H，LIU C S，et al.Index system for ecological effect evaluation of hydropower projects[J].Advances in Science and Technology of Water Resources，2006（6）：1115.（in Chinese）） DOI：10.3880/j.issn.10067647.2006.06.003.

[3]孫宗鳳，董增川.水利工程的生態(tài)效應(yīng)分析[J].水利水電技術(shù)，2004，35（4）：58.（SUN Z F，DONG Z C.Analysis on ecological effect of water conservancy project[J].Water Resources and Hydropower Engineering，2004，35（4）：58.（in Chinese）） DOI：10.3969/j.issn.10000860.2004.04.002.

[4]楊立信，劉國緯.國外調(diào)水工程[M].北京：中國水利水電出版社，2003.（YANG L X，LIU G W.Foreign water transfer project[M].Beijing：China Water & Power Press，2003.（in Chinese））.

[5]呂學(xué)研，吳時(shí)強(qiáng)，張?jiān)?，?調(diào)水引流工程生態(tài)與環(huán)境效應(yīng)研究進(jìn)展[J].水資源與水工程學(xué)報(bào)，2015，26（4）：3845.（LV X Y，WU S Q，ZHANG Y，et al.Progress of ecology and environment effect of water diversion project[J].Journal of Water Resources & Water Engineering，2015，26（4）：3845.（in Chinese）） DOI：10.11705/j.issn.1672643X.2015.04.08.

[6]錢挹清.應(yīng)用模糊綜合評判法進(jìn)行東莞市水資源規(guī)劃宏觀經(jīng)濟(jì)社會(huì)效益評價(jià)[J].珠江現(xiàn)代建設(shè)，2006，12（6）：13，8.（QIAN Y Q.Evaluation of macroeconomy and social benefit of water resources planning in Dongguan City by Fuzzy Comprehensive Evaluation Method[J].Modern Construction of Zhujiang River，2006，12（6）：13，8.（in Chinese））.

[7]陳守煜.區(qū)域水資源可持續(xù)利用評價(jià)理論模型與方法[J].中國工程科學(xué)，2001，3（2）：3338.（CHEN S Y.Theoretical model and method of evaluation of regional water resources sustainable utilization[J].Engineering Sciences，2001，3（2）：3338.（in Chinese）） DOI：10.3969/j.issn.10091742.2001.02.008.

[8]ZADEH L A.Fuzzy Sets[J].Information and Control，1965，8（3）：338353.

[9]陳守煜.工程可變模糊集理論與模型模糊水文水資源學(xué)數(shù)學(xué)基礎(chǔ)[J].大連理工大學(xué)學(xué)報(bào)，2005，45，（2）：308312.（CHEN S Y.The theory and model of engineering variable fuzzy set[J].Journal of Dalian University of Technology，2005，45，（2）：308312.（in Chinese）） DOI：10.3321/j.issn：10008608.2005.02.031.

[10]陳守煜，胡吉敏.可變模糊評價(jià)法及在水資源承載能力評價(jià)中的應(yīng)用[J].水利學(xué)報(bào)，2006，37（3）：264 271.（CHEN S Y，HU J M.Variable fuzzy assessment method and its application in assessing water resources carrying capacity[J].Journal of Hydraulic Engineering，2006，37（3）：264 271.（in Chinese）） DOI：10.3321/j.issn：05599350.2006.03.003.

[11]陳守煜.工程模糊集理論與應(yīng)用[M].北京：國防工業(yè)出版社，1998.（CHEN S Y.Theory and application of engineering fuzzy sets[M].Beijing：National Defense Industry Press，1998.（in Chinese））

[12]陳守煜.開放的復(fù)雜水資源巨系統(tǒng)概念及決策技術(shù)[J].大連理工大學(xué)學(xué)報(bào)，1993，33（5）：591597.（CHEN S Y.Open complex water resources giant system concept and decision technology[J].Journal of Dalian University of Technology，1993，33（5）：591597.（in Chinese））

[13]陳守煜.水資源與防洪系統(tǒng)可變模糊集理論與方法[M].大連：大連理工大學(xué)出版社，2005.（CHEN S Y.The theory and method of variable fuzzy set for water resources and flood control system[M].Dalian：Dalian University of Technology Press，2005.（in Chinese））

[14]王文圣，金菊良，丁晶，等.水資源系統(tǒng)評價(jià)新方法集對評價(jià)法[J].中國科學(xué)E輯：技術(shù)科學(xué)，2009，39（9）：15291534.（WANG W S，JIN J L，DING J，et al.New method of water resources system evaluation set pair evaluation method[J].Science in China（Series E）：Technical Science，2009，39（9）：15291534.（in Chinese））

[15]田林鋼，靳聰聰.水景觀建設(shè)對城市河流生態(tài)健康的影響[J].水力發(fā)電，2014，40（4）：48.（TIAN L G，JIN C C.Influence of the construction of water landscape on urban river ecology health[J].Water Power，2014，40（4）：48.（in Chinese）） DOI：10.3969/j.issn.05599342.2014.04.002.

[16]陸志強(qiáng)，李吉鵬，章耕耘，等.基于可變模糊評價(jià)模型的東山灣生態(tài)系統(tǒng)健康評價(jià)[J].生態(tài)學(xué)報(bào)，2015，35（14）：49074919.（LU Z Q，LI J P，ZHANG G Y，et al.Ecosystem health assessment based on variable fuzzy evaluation model in Dongshan Bay，F(xiàn)ujian，China[J].Acta Ecologica Sinica，2015，35（14）：49074919.（in Chinese）） DOI：10.5846/stxb201312052899.

[17]許雪青，劉忠波，張晉，等.基于可變模糊評價(jià)法的渤海海洋生態(tài)環(huán)境綜合評價(jià)[J].海洋湖沼通報(bào)，2014，（3）：5157.（XU X Q，LIU Z B，ZHANG J，et al.Assessment of the marine ecological environment in Bohai Sea based on a variable fuzzy sets method[J].Transaction of Oceanology and Limnology，2014，（3）：5157.（in Chinese））

[18]曹永強(qiáng)，朱明明，張亮亮，等.基于可變模糊評價(jià)法的大連市水資源承載力分析[J].水利水運(yùn)工程學(xué)報(bào)，2016（4）：4046.（CAO Y Q，ZHU M M，ZHANG L L，et al.Analysis of carrying capacity of water resources in Dalian based on variable fuzzy assessment method[J].HydroScience and Engineering，2016（4）：4046.（in Chinese）） DOI：10.16198/j.cnki.1009640X.2016.04.006.

[19]劉玉玉，許士國.可變模糊評價(jià)方法在渾河上游河段健康評價(jià)中的應(yīng)用[J].水利水電科技進(jìn)展，2013，33（3）：6467，88.（LIU Y Y，XU S G.River health assessment in the upper Hun River based on variable fuzzy model[J].Advances in Science and Technology of Water Resources，2013，33（3）：6467，88.（in Chinese）） DOI：10.3880/j.issn.10067647.2013.03.014.

[20]趙敏敏，張興玨，閆永鑾，等.山東省青州市利用亞行貸款地下水漏斗區(qū)域綜合治理示范工程可行性研究報(bào)告[R].2015.（ZHAO M M，ZHANG X Y，YAN Y L，et al.The report on the feasibility of demonstration project of comprehensive management of groundwater funnel in Qingzhou City，Shandong Province by ADB Loan[R].2015（in Chinese））

[21]鄭良勇，齊春三，宋煒.濰坊市現(xiàn)代水網(wǎng)構(gòu)建研究[J]水利規(guī)劃與設(shè)計(jì)，2005（11）：35，47.（ZHENG L Y，QI C S，SONG W.Study on modern water network construction in Weifang City[J].Water Resources Planning and Design，2005（11）：35，47.（in Chinese）） DOI：10.3969/j.issn.16722469.2015.11.002.

[22]李莉，彭慧，王傳榮，等.山東省青州市現(xiàn)代水網(wǎng)建設(shè)規(guī)劃[R].2012.（LI L，PENG H，WANG C R，et al.Planning of modern water network in Qingzhou City，Shandong Province[R].2012.（in Chinese））

[23]董哲仁.河流生態(tài)系統(tǒng)研究的理論框架[J].水利學(xué)報(bào)，2009，2（2）：129137.（DONG Z R.The theoretical framework of rivers ecosystem research[J].Journal of Hydraulic Engineering，2009，2（2）：129137.（in Chinese）） DOI：10.3321/j.issn：05599350.2009.02.001.

[24]山東省水資源綜合利用中長期規(guī)劃[Z].2016.（Medium and long term plan for comprehensive utilization of water resources in Shandong Province[Z].2006.（in Chinese））

[25]青州統(tǒng)計(jì)年鑒2013[Z].山東省青州市統(tǒng)計(jì)局，2013.（Qingzhou statistical yearbook in 2013[Z].Statistics Bureau of Qingzhou City，Shandong Province，2013.（in Chinese））