基于最大樹法的華北地區(qū)節(jié)水型農(nóng)業(yè)分區(qū)

基于最大樹法的華北地區(qū)節(jié)水型農(nóng)業(yè)分區(qū)劉玉鳳 黃介生 伍靖偉摘要：針對(duì)華北地區(qū)節(jié)水農(nóng)業(yè)發(fā)展的區(qū)域差異性和一致性，根據(jù)地形地貌、氣候特征、缺水程度、灌區(qū)類型、農(nóng)業(yè)結(jié)構(gòu)、經(jīng)濟(jì)水平以及農(nóng)業(yè)節(jié)水發(fā)展現(xiàn)狀等建立農(nóng)業(yè)節(jié)水區(qū)劃的分區(qū)指標(biāo)體系，采用模糊聚類法中的最大樹法將華北地區(qū)分為10個(gè)農(nóng)業(yè)節(jié)水類型區(qū)，并根據(jù)各區(qū)特點(diǎn)提出各自的發(fā)展方向。關(guān)鍵詞：節(jié)水分區(qū)；指標(biāo)體系；最大樹法；模糊聚類中圖分類號(hào)：S27文獻(xiàn)標(biāo)識(shí)碼：A0引言華北地區(qū)水資源嚴(yán)重匱乏，作為全國糧、棉、油的主要產(chǎn)區(qū)，由于水資源短缺，地下水超采嚴(yán)重，已引起一系列的環(huán)境問題1，迫切需要大力發(fā)展節(jié)水農(nóng)業(yè)。為指導(dǎo)農(nóng)業(yè)節(jié)水工作更好地發(fā)展，根據(jù)各分區(qū)的特點(diǎn)，因地制宜地制定節(jié)水農(nóng)業(yè)發(fā)展模式，進(jìn)行農(nóng)業(yè)節(jié)水分區(qū)是十分必要的。國內(nèi)關(guān)于節(jié)水農(nóng)業(yè)分區(qū)，已經(jīng)有經(jīng)驗(yàn)法、指標(biāo)法、類型法、重疊法、聚類法、模糊聚類法、共區(qū)優(yōu)選法等方法2。就華北地區(qū)節(jié)水農(nóng)業(yè)分區(qū)而言，黃修橋等人曾在1995年以地貌特征、氣候特征、灌溉程度、缺水程度作為分區(qū)指標(biāo)，采用模糊聚類和傳統(tǒng)經(jīng)驗(yàn)相結(jié)合的方法將華北分為10個(gè)節(jié)水型農(nóng)業(yè)分區(qū)，并對(duì)2000年華北地區(qū)節(jié)水型農(nóng)業(yè)發(fā)展規(guī)模進(jìn)行了預(yù)測3。本文以此為基礎(chǔ)，根據(jù)最新收集到的華北地區(qū)節(jié)水農(nóng)業(yè)現(xiàn)狀等相關(guān)數(shù)據(jù)資料，進(jìn)一步豐富指標(biāo)體系，采用模糊聚類法中的最大樹法進(jìn)行分區(qū)。1華北地區(qū)節(jié)水型農(nóng)業(yè)分區(qū)11分區(qū)基本單元的確定華北地區(qū)有近300個(gè)縣市，若以此作為分析單元，計(jì)算統(tǒng)計(jì)量相當(dāng)大?？紤]節(jié)水型農(nóng)業(yè)分區(qū)中的一些指標(biāo)在地域上的相似性以及不打破各行政區(qū)劃的目的，初步將華北地區(qū)46個(gè)市級(jí)行政區(qū)作為分區(qū)的基本單元。為了方便，用數(shù)字編號(hào)代表分區(qū)基本單元，數(shù)字代號(hào)與分區(qū)基本單元的對(duì)應(yīng)情況如表1。12分區(qū)指標(biāo)的設(shè)置分區(qū)指標(biāo)的設(shè)置應(yīng)滿足全面性、概括性、易于取得的要求，并能反應(yīng)節(jié)水農(nóng)業(yè)的區(qū)域特征。根據(jù)節(jié)水灌溉特點(diǎn)，結(jié)合華北地區(qū)的具體情況，從眾多影響節(jié)水門農(nóng)業(yè)的指標(biāo)中，篩選出七項(xiàng)對(duì)節(jié)水農(nóng)業(yè)地域分異影響較大的定量指標(biāo)，即干燥度、農(nóng)民人均純收入、耕地灌溉率、畝均水資源量、缺水程度、糧食作物及棉花種植比例。另設(shè)置兩項(xiàng)定性指標(biāo)：地形地貌、灌區(qū)類型。為了便于聚類分析，將定性指標(biāo)定量化，量化標(biāo)準(zhǔn)見表2。各項(xiàng)定量指標(biāo)的確定方法如下。121干燥度K干燥度是一定時(shí)段內(nèi)可能蒸發(fā)量與平均降水量的比值。是反映氣候干濕程度的指標(biāo)，它的計(jì)算采用中國氣候資源及區(qū)劃推薦的經(jīng)驗(yàn)公式即：K=016T/P （1）式中：T 為大于10期間的積溫；P 為同期的降雨量。K 與氣候干濕的關(guān)系：K200干旱，15K200半干旱，10K15半濕潤，05K100濕潤，K05十分濕潤4。122農(nóng)民人均純收入、畝均水資源量農(nóng)民人均純收入以2007年到2010年華北地區(qū)各省市的經(jīng)濟(jì)統(tǒng)計(jì)年鑒的統(tǒng)計(jì)數(shù)據(jù)為準(zhǔn)；畝均水資源量以20072010年華北地區(qū)各省市的水資源公報(bào)的統(tǒng)計(jì)數(shù)據(jù)為準(zhǔn)。123耕地灌溉率耕地灌溉率為有效灌溉面積與耕地面積的比值，即：124缺水程度缺水程度是直接反應(yīng)某地區(qū)灌溉水資源量豐缺程度的指標(biāo)。用耕地面積上的可用水量與綜合作物需水量的比值表示，即 （2）式中：Wy為可用灌溉水量，m3/hm2；Ws為綜合作物需水量，m3/hm2。綜合作物需水量按下式計(jì)算： （3）式中：n為作物種類數(shù)；Wi為第i 種作物的需水量，m3/hm2；ai為第i種作物的種植比例，。與缺水程度的關(guān)系為：1為極缺水區(qū)；12為缺水區(qū)；23為微缺水區(qū)；3為不缺水區(qū)。125糧食作物、棉花種植比例用下列公式計(jì)算：13指標(biāo)權(quán)重的確定指標(biāo)權(quán)重的確定很重要，應(yīng)根據(jù)各項(xiàng)指標(biāo)對(duì)農(nóng)業(yè)節(jié)水分區(qū)的重要性，結(jié)合當(dāng)?shù)氐膶?shí)際情況進(jìn)行權(quán)衡。指標(biāo)權(quán)重的確定方法分主觀賦權(quán)法和客觀賦權(quán)法兩大類。層次分析法是主觀賦權(quán)法的一種，運(yùn)用層次分析法確定指標(biāo)權(quán)重思路簡單明了，不需要樣本數(shù)據(jù)，僅憑對(duì)評(píng)價(jià)指標(biāo)內(nèi)涵與外延的理解即可作出判斷，其缺點(diǎn)是存在主觀性，在一定程度上可能會(huì)降低可信度。熵值法是客觀賦權(quán)法的一種，它客觀地反映了指標(biāo)信息熵的效用價(jià)值，確定的指標(biāo)權(quán)重相對(duì)來講有較高的可信度，但缺乏指標(biāo)間的橫向比較。此外，層次分析法確定的指標(biāo)權(quán)重與實(shí)際調(diào)查樣本信息無關(guān)，即只要評(píng)定者不變，指標(biāo)權(quán)重一旦確定就不會(huì)隨樣本的變化而變化。熵值法確定的指標(biāo)權(quán)重基于樣本的信息進(jìn)行計(jì)算，其值隨樣本的變化而變化5?？傮w上說，無論是主觀賦權(quán)法還是客觀賦權(quán)法，均存在著一定的不足。為避免單一方法的缺陷，本文將主、客觀兩種方法獲得的指標(biāo)權(quán)重有機(jī)結(jié)合起來，采用綜合指標(biāo)權(quán)重，即：wj=w1j+（1-）w2j（j1，2，m；01）（4）式中：wj為第j 項(xiàng)指標(biāo)的綜合指標(biāo)權(quán)重；w1j為第j 項(xiàng)指標(biāo)采用層次分析法得到的指標(biāo)權(quán)重；w2j為采用熵值法得到的評(píng)價(jià)指標(biāo)的權(quán)重。取不同的值即可得到相應(yīng)的綜合指標(biāo)權(quán)重，本文取05，把層次分析法和信息熵法所確定的指標(biāo)權(quán)重代入式（4）即可得到各項(xiàng)指標(biāo)的綜合權(quán)重，結(jié)果見表3。14分區(qū)方法采用模糊聚類分析中的最大樹法進(jìn)行農(nóng)業(yè)節(jié)水分區(qū)。模糊聚類分析是數(shù)理統(tǒng)計(jì)研究中“物以聚類”的一種多元分析方法，即用數(shù)學(xué)定量地確定樣品的親疏關(guān)系，從而客觀地分型劃分6。最大樹法是模糊聚類分析方法中較為經(jīng)典的一種方法，做出的樹狀圖能生動(dòng)的反映樣品間的相似關(guān)系，由給定的進(jìn)行切割最大樹，就可以進(jìn)行分類，這樣更為直觀方便，特別適合在計(jì)算機(jī)上進(jìn)行運(yùn)算聚類。其具體步驟如下：（1） 原始矩陣構(gòu)成。以基本單元為樣本，各項(xiàng)指標(biāo)作為描述樣本的特征指標(biāo)，一個(gè)樣本為一行，所有單元構(gòu)成模糊聚類的原始矩陣A。（5）其中：xij=Wjaij（i=1，2，46；j1，2，9）（6）式中：aij表示第i個(gè)基本單元對(duì)應(yīng)的第j 個(gè)分區(qū)指標(biāo)的指標(biāo)值；xij表示考慮指標(biāo)權(quán)重后的基本單元的指標(biāo)值；Wj表示第j 個(gè)分區(qū)指標(biāo)的綜合權(quán)重。（2）數(shù)據(jù)標(biāo)準(zhǔn)化。為了消除指標(biāo)量綱的差別和數(shù)量級(jí)不同的影響，需對(duì)各指標(biāo)值進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化，使每一個(gè)指標(biāo)都統(tǒng)一于某種共同的數(shù)值特性范圍，以便于比較和分析。數(shù)據(jù)標(biāo)準(zhǔn)化的公式為： （7）式中：xij為經(jīng)標(biāo)準(zhǔn)化處理后的各項(xiàng)指標(biāo)值（表略）；xij為原始矩陣中的指標(biāo)值；xj為原始矩陣中指標(biāo)值的平均值；c為原始矩陣中指標(biāo)值的標(biāo)準(zhǔn)差。（3）建立模糊相似關(guān)系矩陣R。算出衡量被分類基本單元間相似程度的統(tǒng)計(jì)量ij（i，j1，2，n），n為被分類基本單元的個(gè)數(shù)，從而確定模糊相似關(guān)系：（8）（9）式中：ij（i1，n；j1，n）為模糊相似關(guān)系矩陣R 的第i行第j 列元素，表示被分類基本單元間的相似程度；c為一特定的正數(shù)，用以保證0ij1，經(jīng)分析確定c0040；xik，xjk分別為數(shù)據(jù)標(biāo)準(zhǔn)化后矩陣中的第i行，第j行k 列元素。（4）模糊聚類分析。采用模糊聚類分析中的最大樹法，首先畫出最大樹，然后選定閾值切割最大樹進(jìn)行動(dòng)態(tài)聚類，具體步驟如下：用Kruskal法畫最大樹7，找出R 中最大值0941，對(duì)應(yīng)的頂點(diǎn)分別是12，18；連邊；標(biāo)權(quán)重，依次按照ij從大到小順序連邊標(biāo)權(quán)重，直到所有的頂點(diǎn)都出現(xiàn)在圖中。如圖1所示，該圖表示各個(gè)單元的相似程度?？橙D中權(quán)重小于073的邊，最大樹分為互不相連的三棵樹，分別是1、2和剩余所有，考慮到1、2的行政獨(dú)特性，將兩者歸為一類，故此時(shí)分為兩類，1、2歸為一類，其余相連的各區(qū)為一類；砍去圖中權(quán)重小于075的邊，則相連的43、44、45歸為一類，其余相連的各點(diǎn)為一類；砍去圖中權(quán)重小于077的邊，則相連的9、13、15歸為一類，其余相連的各點(diǎn)為一類；砍去圖中權(quán)重小于078的邊，則相連的11、10、12、5、18、26、25、42、38、40、41、39、14歸為一類，其余相連的各點(diǎn)為一類；砍去圖中權(quán)重小于084的邊，則相連的11、10、12、5、18、26、25、42、38、40、41、39、14分成兩類，即11、10、12、5、18、26、25為一類，42、38、40、41、39、14為另一類，原來的剩余相連部分成了四類：29、30、31、33、27、36、34為一類，20、35、37、28、32為第二類，19、22、46為第三類，其余部分為第四類；砍去圖中權(quán)重小于085的邊，則剩下相連部分分為兩類：24、17、16、21、23歸為一類，3、4、6、7、8歸為另一類。15分區(qū)命名和結(jié)果分區(qū)命名由地理位置、地貌形態(tài)、氣候特征、缺水程度、耕地灌溉率等因素組成，根據(jù)最大樹法得到的結(jié)果，并考慮地域完整性，將華北地區(qū)劃分為10個(gè)分區(qū)，分區(qū)名稱及有關(guān)指標(biāo)值見表4。華北地區(qū)節(jié)水農(nóng)業(yè)分區(qū)見圖2。2分區(qū)特點(diǎn)及農(nóng)業(yè)節(jié)水發(fā)展方向（1）區(qū)。本區(qū)包括北京和天津兩大直轄市，降雨量相對(duì)較大，屬于華北地區(qū)相對(duì)濕潤、水資源較為豐富的地區(qū)。該地區(qū)社會(huì)經(jīng)濟(jì)發(fā)展水平較高，主要種植小麥、玉米等糧食作物及大豆、花生等經(jīng)濟(jì)作物，農(nóng)田灌溉方式為蓄水灌溉和提水灌溉。此區(qū)菜園、果園、牧草及一些新增作物的種植面積發(fā)展迅速，但田間節(jié)水未得到足夠重視，今后田間灌溉應(yīng)該朝著減少地面無效蒸發(fā)和深層滲漏方向發(fā)展，可以對(duì)地下滴灌、滲灌等節(jié)水灌溉技術(shù)進(jìn)行研究，技術(shù)設(shè)備成熟后進(jìn)行推廣。（2）區(qū)。本區(qū)包括內(nèi)蒙古的呼和浩特、鄂爾多斯、包頭三個(gè)地級(jí)市，降雨量少，大部分屬干旱、半干旱的大陸性季風(fēng)氣候，水資源短缺，水土流失嚴(yán)重。主要種植小麥、玉米等糧食作物及大豆等經(jīng)濟(jì)作物，農(nóng)田灌溉方式為引水灌溉和提水灌溉，耕地次生鹽堿化問題突出，改善農(nóng)作物結(jié)構(gòu)如擴(kuò)大耐鹽堿作物的種植面積將提高農(nóng)業(yè)用水的效率。本區(qū)大多處于山地，坡耕地比重較大，但土層薄，保水保肥能力差，結(jié)合山溪性小溪及地面坡降較大的特點(diǎn)，宜重點(diǎn)發(fā)展自壓噴灌技術(shù)。（3）區(qū)。本區(qū)包括山東省的日照、臨沂、萊蕪三個(gè)地級(jí)市，處于北溫帶半濕潤季風(fēng)氣候區(qū)，降雨量較大，水資源年際年內(nèi)變化劇烈，水資源地區(qū)分布不均，開發(fā)利用難度較大。本區(qū)主要種植小麥、玉米等糧食作物及大豆、花生等經(jīng)濟(jì)作物，農(nóng)田灌溉方式為蓄水灌溉和引水灌溉。本區(qū)降雨量相對(duì)較大，應(yīng)當(dāng)重視農(nóng)業(yè)節(jié)水配套技術(shù)，加強(qiáng)雨水資源的利用。（4）區(qū)。本區(qū)包括山東的德州、濱州、聊城、東營和河北的邯鄲、衡水、邢臺(tái)七個(gè)地級(jí)市，灌區(qū)多處于海河下游平原地區(qū)，降雨量較少，氣候濕潤，水資源短缺。主要種植小麥、玉米等糧食作物和棉花等經(jīng)濟(jì)作物，農(nóng)田灌溉方式主要為引水灌溉和蓄水灌溉。本區(qū)地形平緩，田面平整度不高，宜加大平整土地和田間工程的建設(shè)，本區(qū)棉花、園藝等特色經(jīng)濟(jì)作物的比重較大，宜穩(wěn)步發(fā)展膜下軟管灌溉、棉花膜下滴灌、園藝滴灌等高效節(jié)水技術(shù)。（5）區(qū)。本區(qū)包括河南的安陽、濮陽、新鄉(xiāng)、焦作、鶴壁和山東的菏澤六個(gè)地級(jí)市，降雨量較少，氣候相對(duì)濕潤，水資源總量小且利用率低，缺水嚴(yán)重。主要種植小麥、玉米等糧食作物，農(nóng)田灌溉方式以引水灌溉為主，蓄水灌溉和井灌為輔。本區(qū)地勢平坦，田面平整度不高，應(yīng)重點(diǎn)發(fā)展以精細(xì)平地技術(shù)為核心的水平畦灌和覆膜灌溉技術(shù)。本區(qū)旱地較多，應(yīng)當(dāng)大力推廣旱作節(jié)水技術(shù)體系，推廣積水工程及集雨補(bǔ)灌旱作節(jié)水農(nóng)業(yè)，通過農(nóng)業(yè)工程和農(nóng)藝措施，將自然降水在時(shí)間和空間上進(jìn)行人工再分配，提高利用效率。（6）區(qū)。本區(qū)包括山西的陽泉、長治、晉城、晉中、太原、運(yùn)城、臨汾七個(gè)地級(jí)市，降雨量稀少，蒸發(fā)量大，氣候干燥，區(qū)域大部分屬于干旱或半干旱地區(qū)，水資源量少，缺水嚴(yán)重。主要種植玉米等糧食作物，農(nóng)田灌溉方式主要為引水灌溉。此區(qū)最大的特點(diǎn)就是干旱，今后應(yīng)該從提高天然降水利用率和效率、增強(qiáng)抗旱節(jié)水增產(chǎn)能力出發(fā)，改善不合理的耕作方式和栽培模式，利用各種蓄水保墑技術(shù)，建設(shè)田間土壤水庫，提高田間水分的生產(chǎn)效益。（7）區(qū)。本區(qū)包括山西的忻州、呂梁、大同、朔州和河北的張家口五個(gè)地級(jí)市，降水量較少，蒸發(fā)量大，屬于華北地區(qū)相對(duì)濕潤地區(qū)，水資源短缺，主要種植玉米等糧食作物，農(nóng)田灌溉方式主要為蓄水灌溉。本區(qū)大多處于山地，坡耕地比重較大，但土層薄，保水保肥能力差，結(jié)合山溪性小溪及地面坡降較大的特點(diǎn)，宜重點(diǎn)發(fā)展自壓噴灌技術(shù)。（8）區(qū)。本區(qū)包括河北的秦皇島、承德和內(nèi)蒙古的赤峰三個(gè)地級(jí)市，降水量較少，蒸發(fā)量大，氣候干燥，區(qū)域大部分屬于干旱或半干旱地區(qū)，水資源少，缺水現(xiàn)象非常嚴(yán)重。主要種植玉米等糧食作物和花生等經(jīng)濟(jì)作物，農(nóng)田灌溉方式主要為蓄水灌溉。本區(qū)大多處于山丘區(qū)，改進(jìn)農(nóng)作物的灌溉制度，積極開展渠道防滲，減少渠系滲漏損失，是節(jié)水灌溉的重要方面。（9）區(qū)。本區(qū)包括河北的石家莊、保定、唐山、廊坊、滄州五個(gè)地級(jí)市，降雨量相對(duì)較大，屬于華北地區(qū)相對(duì)濕潤地區(qū)，水資源短缺。主要種植小麥、玉米等糧食作物，農(nóng)田灌溉方式以蓄水灌溉為主、以引水灌溉為輔。本區(qū)大多處于平丘區(qū)，渠灌是主要的灌溉方式之一，節(jié)水潛力較大，調(diào)整渠系，渠道防滲，更新改造老化工程，加強(qiáng)灌溉管理，優(yōu)化配水等綜合技術(shù)改造措施將實(shí)現(xiàn)節(jié)水增產(chǎn)。（10）區(qū)。本區(qū)包括山東省的濟(jì)南、淄博、濰坊、濟(jì)寧、泰安五個(gè)地級(jí)市，降雨量相對(duì)較大，氣候濕潤，水資源年際年內(nèi)變化劇烈，水資源地區(qū)分布不均，開發(fā)利用難度較大。本區(qū)主要種植小麥、玉米等糧食作物及蔬菜等經(jīng)濟(jì)作物，農(nóng)田灌溉方式以引水灌溉為主、蓄水灌溉為輔。本區(qū)降雨量相對(duì)較大但有效利用率較小，今后應(yīng)重視農(nóng)業(yè)節(jié)水配套技術(shù)，加強(qiáng)雨水資源的利用。本區(qū)地勢較為平坦，田面平整度不高，宜重點(diǎn)發(fā)展以精細(xì)平地技術(shù)為核心的水平畦灌和覆膜灌溉技術(shù)。本文注重全面性的考慮，將黃修橋等人之前的四指標(biāo)擴(kuò)展為九指標(biāo)體系，雖然同分為10個(gè)分區(qū)，但具體分區(qū)結(jié)果存在著差異；黃修橋等人對(duì)節(jié)水型農(nóng)業(yè)發(fā)展大體規(guī)模進(jìn)行了數(shù)字上的預(yù)測，但未能提出各區(qū)各自的節(jié)水發(fā)展具體措施，本文盡力避免這一不足，在對(duì)分區(qū)的特點(diǎn)進(jìn)行綜合分析的基礎(chǔ)上，因地制宜地提出具體的節(jié)水發(fā)展措施，用于節(jié)水灌溉實(shí)踐。3結(jié)語根據(jù)華北地區(qū)地形地貌、氣候特征、缺水程度、灌區(qū)類型、農(nóng)業(yè)結(jié)構(gòu)、經(jīng)濟(jì)水平以及農(nóng)業(yè)節(jié)水發(fā)展現(xiàn)狀等組成農(nóng)業(yè)節(jié)水區(qū)劃的分區(qū)指標(biāo)體系，采用模糊聚類法中的最大樹法將華北地區(qū)分為十個(gè)類型區(qū)，并根據(jù)各區(qū)特點(diǎn)提出各自的發(fā)展方向。本文采用了定量定性相結(jié)合的辦法，根據(jù)華北地區(qū)實(shí)際情況定性確定了指標(biāo)體系，然后對(duì)定性指標(biāo)進(jìn)行等級(jí)及評(píng)分，后將定性指標(biāo)定量化，指標(biāo)的權(quán)重的設(shè)置采用主觀和客觀相結(jié)合的辦法，并且對(duì)每個(gè)分區(qū)的綜合指標(biāo)值進(jìn)行了嚴(yán)格的計(jì)算。本文限于資料數(shù)據(jù)以地級(jí)市為單位，故只相應(yīng)地進(jìn)行了一級(jí)區(qū)的劃分，今后可以在完善資料的基礎(chǔ)上進(jìn)行二級(jí)分區(qū)，也可參照各省市已有的分區(qū)資料進(jìn)行更深層次的分區(qū)。在華北地區(qū)未來的節(jié)水灌溉規(guī)劃中，可以參考這個(gè)分區(qū)，相信對(duì)促進(jìn)華北地區(qū)節(jié)水農(nóng)業(yè)發(fā)展的有序化和系統(tǒng)化將發(fā)揮一定的指導(dǎo)作用。參考文獻(xiàn)：1王道波，張廣錄，周曉果華北水資源利用現(xiàn)狀及其宏觀調(diào)控對(duì)策研究J干旱區(qū)資源與環(huán)境，2005，19（2）：46502韓棟節(jié)水灌溉規(guī)劃分區(qū)方法D北京：北京工業(yè)大學(xué)，20053黃修橋，李英能華北地區(qū)節(jié)水型農(nóng)業(yè)分區(qū)和發(fā)展預(yù)測J灌溉排水，1995，14（3）：184水利部農(nóng)村水利司水土資源評(píng)價(jià)與節(jié)水灌溉規(guī)劃M北京：中國水利水電出版社，19885舒衛(wèi)萍，崔遠(yuǎn)來層次分析法在灌區(qū)綜合評(píng)價(jià)中的應(yīng)用J中國農(nóng)村水利水電，2005，（6）：1091116董媛媛模糊聚類分析的兩種方法J井岡山學(xué)院學(xué)報(bào)（自然科學(xué)），2006，27（6）：26287王紅霞，盧文喜河北省節(jié)水灌溉區(qū)劃中模糊聚類分析與應(yīng)用J工程勘察，2006（12）：3235作者簡介：劉玉鳳（1989），女，碩士，研究方向?yàn)樗临Y源。