王 子 茹， 羅 寶 力， 牛 云 格

(大連理工大學(xué) 建設(shè)工程學(xué)部，遼寧 大連 116024)

0 引 言

水是基礎(chǔ)性的自然資源和戰(zhàn)略性的經(jīng)濟(jì)資源，是生態(tài)環(huán)境的控制要素.隨著信息技術(shù)的飛速發(fā)展，水利信息化建設(shè)在水資源管理中的重要性日益凸顯.計(jì)算機(jī)技術(shù)、科學(xué)可視化滲透在整個(gè)水資源管理系統(tǒng)中，為實(shí)現(xiàn)水資源的動(dòng)態(tài)監(jiān)測(cè)、水資源的調(diào)查和評(píng)價(jià)、水資源可持續(xù)利用和優(yōu)化配置以及提高水資源的管理水平，提供了前所未有的技術(shù)保障和決策支撐手段.

近年來(lái)，利用計(jì)算機(jī)可視化技術(shù)進(jìn)行水資源管理方面已開(kāi)展了大量的研究工作，文獻(xiàn)[1、2]應(yīng)用地理信息和可視化技術(shù)，建設(shè)了市級(jí)水資源管理系統(tǒng)，實(shí)現(xiàn)了區(qū)域水資源的動(dòng)態(tài)監(jiān)測(cè)以及水資源的評(píng)價(jià)、分區(qū)等；文獻(xiàn)[3]以VB為開(kāi)發(fā)平臺(tái)，以區(qū)域水資源經(jīng)濟(jì)凈效益最大化為目標(biāo)函數(shù)，建立了多目標(biāo)規(guī)劃模型，得到水資源利用的最佳方案；文獻(xiàn)[4]將GIS應(yīng)用于地下水模擬、評(píng)價(jià)和管理，為地下水信息化工作引進(jìn)了新的技術(shù)手段.

上述研究無(wú)疑推動(dòng)了我國(guó)信息技術(shù)在水資源管理方面的應(yīng)用和發(fā)展，為數(shù)字水利提供了先進(jìn)的技術(shù)支持.但就已有研究成果中，應(yīng)用陳守煜于2005年創(chuàng)建的可變模糊集理論[5～7]，對(duì)區(qū)域水資源配置方案進(jìn)行評(píng)價(jià)優(yōu)選，并將其以可視化的形式表達(dá)尚未見(jiàn)到.本文將VC＋＋可視化編程技術(shù)和GIS結(jié)合，采用文獻(xiàn)[8]的評(píng)價(jià)方法，建立區(qū)域水資源優(yōu)化配置決策方案評(píng)價(jià)可視化系統(tǒng)，并以大連市水資源優(yōu)化配置為例，對(duì)多種方案進(jìn)行優(yōu)選，為水資源合理配置提供決策支持.

1 水資源優(yōu)化配置決策方案評(píng)價(jià)模型簡(jiǎn)介[8]

設(shè)xi(i＝1，2，…，m)為已知待評(píng)對(duì)象u的指標(biāo)特征值，依據(jù)m個(gè)指標(biāo)c個(gè)級(jí)別的指標(biāo)標(biāo)準(zhǔn)矩陣M＝(Mih)(i＝1，2，…，m；h＝1，2，…，c)進(jìn)行綜合評(píng)價(jià).

根據(jù)文獻(xiàn)[8]，評(píng)價(jià)對(duì)象u指標(biāo)i對(duì)級(jí)別h的相對(duì)隸屬度公式為

又根據(jù)文獻(xiàn)[8]評(píng)價(jià)對(duì)象u對(duì)級(jí)別h的多指標(biāo)綜合相對(duì)隸屬度模型為

α為模型優(yōu)化準(zhǔn)則參數(shù)，當(dāng)α＝1時(shí)，優(yōu)化準(zhǔn)則相當(dāng)于最小一乘方；α＝2時(shí)，優(yōu)化準(zhǔn)則相當(dāng)于最小二乘方.p為距離參數(shù)，當(dāng)p＝1時(shí)，為海明距離；p＝2時(shí)，為歐氏距離.兩者可組成模型的4種組合:(1)α＝1，p＝1；(2)α＝1，p＝2；(3)α＝2，p＝1；(4)α＝2，p＝2.

應(yīng)用文獻(xiàn)[6]中的級(jí)別特征值公式

應(yīng)用文獻(xiàn)[7]中歸屬等級(jí)判斷準(zhǔn)則

2 系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)

2.1 系統(tǒng)的總體結(jié)構(gòu)和功能

基于Windows操作平臺(tái)，利用VC＋＋中的MFC框架和GIS功能組件 MapObjects2.0集成，以大連市水資源配置方案為依據(jù)，應(yīng)用可變模糊評(píng)價(jià)模型，開(kāi)發(fā)建立了水資源優(yōu)化配置可變模糊評(píng)價(jià)可視化系統(tǒng)，系統(tǒng)總體框架如圖1所示.該系統(tǒng)由方案可變模糊評(píng)價(jià)子系統(tǒng)、信息查詢(xún)子系統(tǒng)和GIS水系圖子系統(tǒng)組成.系統(tǒng)的主要功能如下:

圖1 系統(tǒng)的總體結(jié)構(gòu)Fig.1 The general structure of the system

(1)方案可變模糊評(píng)價(jià)子系統(tǒng)

根據(jù)陳守煜最新提出的基于辯證法三大規(guī)律之一的質(zhì)量互變規(guī)律的質(zhì)量互變定理的水資源可持續(xù)利用可變模糊集模型，采用模塊化和參數(shù)化的設(shè)計(jì)理念，通過(guò)建立方案參數(shù)、標(biāo)準(zhǔn)區(qū)間、相對(duì)隸屬度、確定權(quán)重以及評(píng)價(jià)結(jié)果等菜單，開(kāi)發(fā)出了基于可變模糊的水資源優(yōu)化配置評(píng)價(jià)子系統(tǒng).并以大連市水資源優(yōu)化配置方案為例，對(duì)大連市2020年水資源優(yōu)化配置方案進(jìn)行評(píng)價(jià)優(yōu)選.

(2)信息查詢(xún)子系統(tǒng)

以VC和MapObjects相結(jié)合的形式開(kāi)發(fā)了水資源優(yōu)化配置可變模糊評(píng)價(jià)系統(tǒng)中信息查詢(xún)子系統(tǒng).主要功能是:

①可變模糊評(píng)價(jià)信息查詢(xún)(包括水資源優(yōu)化配置方案信息查詢(xún)、指標(biāo)標(biāo)準(zhǔn)矩陣查詢(xún)、可變模糊評(píng)價(jià)模型查詢(xún)、評(píng)價(jià)等級(jí)查詢(xún)等).可以根據(jù)分析要求，查詢(xún)到不同參數(shù)下的級(jí)別特征值.

②大連市各行政區(qū)水資源狀況的空間屬性查詢(xún)，主要是水資源管理業(yè)務(wù)查詢(xún)(水資源基礎(chǔ)資料:地表水、地下水、海水、中水、外調(diào)水等).

③系統(tǒng)考慮到MapObjects組件的功能特點(diǎn)，將地理信息與屬性信息有機(jī)地結(jié)合在一起，實(shí)現(xiàn)圖文互動(dòng)的功能，更直觀地查詢(xún)水資源狀況.

(3)GIS水系圖顯示子系統(tǒng)

主要實(shí)現(xiàn)GIS的基本地圖操作功能，包括地圖的縮放、漫游和全圖顯示等，可使決策者對(duì)大連市轄區(qū)范圍內(nèi)的水資源分配情況有更加直觀、準(zhǔn)確的了解.

2.2 系統(tǒng)主界面設(shè)計(jì)

本系統(tǒng)基于MFC應(yīng)用程序向?qū)9]，以CFormView類(lèi)為基類(lèi)，建立具有參數(shù)化功能的水資源優(yōu)化配置方案可變模糊評(píng)價(jià)系統(tǒng).通過(guò)插入MapObjects的Map控件來(lái)進(jìn)行GIS圖形顯示，系統(tǒng)主界面如圖2所示.

圖2 系統(tǒng)主界面Fig.2 The system main interface

2.3 VC＋＋和MapObjects的集成

2.3.1 生成水系圖 在ArcMap下繪制大連市水系圖[10]，將GIS組件MapObjects與VC相結(jié)合，將繪制好的水系圖通過(guò)Map控件添加到VC環(huán)境中，實(shí)現(xiàn)系統(tǒng)的集成，實(shí)現(xiàn)大連市水系圖在VC環(huán)境下的可視化顯示，并且生成的水系圖具有平移、縮放、漫游和地圖標(biāo)注等功能.具體做法如下:

(1)空間數(shù)據(jù)的建立

在地理信息系統(tǒng)中，空間數(shù)據(jù)與屬性數(shù)據(jù)通過(guò)拓?fù)潢P(guān)系來(lái)建立聯(lián)系.本系統(tǒng)以空間數(shù)據(jù)為研究基礎(chǔ).首先啟動(dòng)ArcCatalog，在ArcCatalog目錄窗口中，選擇New／Shapefile…，進(jìn)入新建Shapefile對(duì)話框.根據(jù)空間數(shù)據(jù)劃分的情況，鍵入Name內(nèi)容和選擇Point、Polyline、Polygon等要素類(lèi)，以及進(jìn)行要素屬性表(Feature Attribute Table)的編輯.然后，啟動(dòng)ArcMap，選擇主菜單File／Add Data，加載新建的Shapefile，進(jìn)行水系圖的編輯.

(2)利用 MFC AppWizard創(chuàng)建應(yīng)用程序

本系統(tǒng)開(kāi)發(fā)一個(gè)基于MFC的單文檔(SDI)應(yīng)用程序.選擇 MFC AppWizard(exe)來(lái)創(chuàng)建應(yīng)用程序.

(3)添加 MapObjects控件

選 擇 菜 單 Project→ Add to Project→Components and Controls Gallery， 加 入MapObjects控件模板，將Map控件添加到Form中，并定義變量.

(4)加載圖層

通過(guò)設(shè)置Map Control控件的屬性對(duì)話框來(lái)加載要顯示的圖層數(shù)據(jù)[11].利用Class Wizard向CWaterMapView中添加 WM＿SIZE事件，在OnSize方法的函數(shù)體中添加代碼.

(5)地圖注記

MapObjects中用于添加注記的對(duì)象是CMoLabelRender對(duì)象和CMoLabelPlacer對(duì)象，這兩個(gè)對(duì)象是常用的注記標(biāo)注工具.本文采用CMoLabelPlacer對(duì)象，可以自動(dòng)添加圖層上對(duì)象的標(biāo)注.

2.3.2 實(shí)現(xiàn)水資源優(yōu)化配置情況的查詢(xún) 以VC和MapObjects相結(jié)合的形式開(kāi)發(fā)了水資源優(yōu)化配置可變模糊評(píng)價(jià)系統(tǒng)中的信息查詢(xún)子系統(tǒng)，實(shí)現(xiàn)了大連市各行政區(qū)的水資源狀況的空間屬性查詢(xún)以及水資源配置方案的信息查詢(xún)等.

查詢(xún)方式有兩種:菜單選擇和GIS地圖選擇.通過(guò)菜單下拉項(xiàng)選擇查詢(xún)目標(biāo)，操作簡(jiǎn)單；地圖選擇則是在GIS地圖上選擇需要查詢(xún)的行政區(qū)域的水資源狀況，鼠標(biāo)點(diǎn)擊水系圖上任意圖層，則該圖層閃亮顯示，彈出一個(gè)窗口，顯示該圖層的屬性信息，為該區(qū)域的水資源狀況.實(shí)現(xiàn)查詢(xún)區(qū)域閃爍效果使用的是MapObjects組件的FlashShape()方法[12].

2.4 系統(tǒng)應(yīng)用——以大連市為例

2.4.1 數(shù)據(jù)資料 根據(jù)大連市的社會(huì)、經(jīng)濟(jì)和環(huán)境等因素，篩選出與水資源可持續(xù)利用和大連社會(huì)經(jīng)濟(jì)發(fā)展有關(guān)的9個(gè)評(píng)價(jià)指標(biāo)[13]進(jìn)行評(píng)價(jià):工業(yè)用水比例x1(%)，中水回用成本x2(106元)，外調(diào)水供水成本x3(106元)，地下水開(kāi)發(fā)利用程度x4(%)，中水回用率x5(%)，污水處理率x6(%)，非傳統(tǒng)水源占總供水比例x7(%)，本地淡水資源利用率x8(%)，外調(diào)水量占總供水量比例x9(%).其中，外調(diào)水量占總供水量比例是外調(diào)水利用量／總供水量.考慮現(xiàn)在的成本有所提高，故將評(píng)價(jià)指標(biāo)標(biāo)準(zhǔn)值x2、x3做了相應(yīng)的調(diào)整.評(píng)價(jià)指標(biāo)及各級(jí)標(biāo)準(zhǔn)值見(jiàn)表1.采用文獻(xiàn)[13]的有關(guān)數(shù)據(jù)，大連市2020年城鎮(zhèn)水資源優(yōu)化配置方案評(píng)價(jià)指標(biāo)值如表2所示.

表1 指標(biāo)分級(jí)標(biāo)準(zhǔn)值Tab.1 Index classification standard values

表2 大連市水資源配置方案評(píng)價(jià)指標(biāo)值Tab.2 Water resources allocation plan evaluation index values in Dalian

采用文獻(xiàn)[13]的指標(biāo)權(quán)重向量為

2.4.2 可變模糊評(píng)價(jià)子系統(tǒng)的實(shí)現(xiàn) 可變模糊評(píng)價(jià)子系統(tǒng)由方案參數(shù)輸入、指標(biāo)分級(jí)標(biāo)準(zhǔn)值、指標(biāo)相對(duì)隸屬度計(jì)算、指標(biāo)權(quán)重、評(píng)價(jià)結(jié)果等系統(tǒng)菜單組成.評(píng)價(jià)時(shí)，在相應(yīng)的對(duì)話框中輸入?yún)?shù)，即可得到該方案的評(píng)價(jià)結(jié)果(各種對(duì)話框略).

以方案1為例，通過(guò)可變模糊評(píng)價(jià)系統(tǒng)對(duì)其進(jìn)行評(píng)價(jià)，最終將4組參數(shù)下的評(píng)價(jià)結(jié)果匯總，得出方案1的評(píng)價(jià)等級(jí)，如圖3所示.

圖3 評(píng)價(jià)結(jié)果匯總Fig.3 Integration of evaluation result

2.4.3 評(píng)價(jià)結(jié)果分析 應(yīng)用可變模糊評(píng)價(jià)方法，采用可變模糊集中的4組模型分別對(duì)3套方案進(jìn)行評(píng)判，評(píng)價(jià)結(jié)果見(jiàn)表3.

表3 大連市水資源優(yōu)化配置方案評(píng)價(jià)等級(jí)Tab.3 Rank of water resources optimizing allocation plans in Dalian

根據(jù)判斷準(zhǔn)則式(5)，通過(guò)分析評(píng)價(jià)結(jié)果可知，3套方案的評(píng)價(jià)等級(jí)均為2級(jí)偏1級(jí).但考慮到水資源系統(tǒng)是一個(gè)復(fù)雜的非線性系統(tǒng)，又考慮到在所選取的評(píng)價(jià)指標(biāo)中有線性指標(biāo)(例如中水回用成本等)，故采用α＝1，p＝1和α＝2，p＝1兩種參數(shù)組合的平均值作為最后的評(píng)價(jià)結(jié)果，如表4所示.由此可見(jiàn)，方案2的評(píng)價(jià)結(jié)果與其他方案相比較，相對(duì)于最優(yōu)級(jí)的隸屬度較高，故方案2的配置方案為最優(yōu)方案.

表4 兩種參數(shù)組合下的評(píng)價(jià)結(jié)果Tab.4 Evaluation result of two kinds of parameters combination

2.4.4 系統(tǒng)查詢(xún)

(1)菜單選擇查詢(xún)信息:通過(guò)菜單下拉項(xiàng)選擇查詢(xún)目標(biāo)，圖4為指標(biāo)值查詢(xún).

(2)選擇GIS地圖查詢(xún)空間屬性:在GIS地圖上可以任意選擇需要查詢(xún)的行政區(qū)域的水資源狀況，點(diǎn)擊鼠標(biāo)左鍵即可實(shí)現(xiàn)水資源狀況的空間屬性查詢(xún).圖5為金州以南地區(qū)的空間屬性查詢(xún).

圖4 指標(biāo)值查詢(xún)Fig.4 Inquirement of index values

圖5 空間屬性查詢(xún)Fig.5 Inquirement of spatial attribute

3 結(jié) 語(yǔ)

利用可視化技術(shù)，以Windows系統(tǒng)為操作平臺(tái)，通過(guò)VC＋＋6.0中的MFC框架結(jié)合GIS技術(shù)，開(kāi)發(fā)出了大連市水資源優(yōu)化配置方案評(píng)價(jià)可視化系統(tǒng).根據(jù)文獻(xiàn)[13、14]中的相關(guān)數(shù)據(jù)，系統(tǒng)實(shí)現(xiàn)了運(yùn)用可變模糊集理論及模型對(duì)大連市城鎮(zhèn)水資源優(yōu)化配置方案的評(píng)價(jià)優(yōu)選，得出最優(yōu)方案.文獻(xiàn)[15]同樣說(shuō)明了可變模糊集理論、模型與方法的科學(xué)性.實(shí)例表明，該系統(tǒng)運(yùn)用于水資源優(yōu)化配置方案評(píng)價(jià)，具有復(fù)雜問(wèn)題簡(jiǎn)單化、計(jì)算方便快捷、人機(jī)界面交互性好的特點(diǎn)，為水資源優(yōu)化配置方案評(píng)價(jià)可視化提供了一個(gè)新平臺(tái).

致謝:本文得到了陳守煜教授的悉心指導(dǎo).

[1]高芳琴，潘崇倫，邱紹偉.上海市水資源綜合管理系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)[J].水利信息，2010(3):65－69

[2]劉秀花，白峰青.可視化技術(shù)在咸陽(yáng)市水資源管理系統(tǒng)中的應(yīng)用[J].干旱地區(qū)農(nóng)業(yè)研究，2005，23(2):191－194

[3]孟曉路，梁秀娟，王 靜，等.吉林中部地區(qū)水資源管理決策支持系統(tǒng)的開(kāi)發(fā)[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版)，2006，11(36):39－42

[4]李門(mén)樓，胡 成，陳植華.河北平原區(qū)域地下水資源決策支持系統(tǒng)設(shè)計(jì)與開(kāi)發(fā)[J].中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，2002，27(2):222－226

[5]陳守煜.工程可變模糊集理論與模型——模糊水文水資源學(xué)數(shù)學(xué)基礎(chǔ)[J].大連理工大學(xué)學(xué)報(bào)，2005，45(2):308－312(CHEN Shou－yu.Theory and model of engineering variable fuzzy set — Mathematical basis for fuzzy hydrology and water resources[J].Journal of Dalian University of Technology，2005，45(2):308－312)

[6]陳守煜.工程模糊集理論與應(yīng)用[M].北京:國(guó)防工業(yè)出版社，1998

[7]陳守煜.可變模糊集理論與模型及其應(yīng)用[M].大連:大連理工大學(xué)出版社，2009

[8]陳守煜.基于可變模糊集的辯證法三大規(guī)律數(shù)學(xué)定理及其應(yīng)用[J].大連理工大學(xué)學(xué)報(bào)，2010，50(5):838－844(CHEN Shou－yu.Three mathematical theorems of dialectics based on variable fuzzy sets and their application [J].Journal of Dalian University of Technology，2010，50(5):838－844)

[9]陳志泊，張海燕，王春玲.Visual C＋＋程序設(shè)計(jì)[M].北京:中國(guó)鐵道出版社，2005

[10]宋小冬，鈕心毅.地理信息系統(tǒng)實(shí)習(xí)教程(ArcGIS 9.x)[M].北京:科學(xué)出版社，2007

[11]韓 鵬.地理信息系統(tǒng)開(kāi)發(fā)——MapObjects方法[M].武漢:武漢大學(xué)出版社，2004

[12]齊 紅.基于VC和MapObjects的森林資源管理系統(tǒng)的開(kāi)發(fā)應(yīng)用研究[D].哈爾濱:東北林業(yè)大學(xué)，2005

[13]黃 旭.大連市水資源合理配置及其方案綜合評(píng)價(jià)研究[D].大連:大連理工大學(xué)，2008

[14]大連理工大學(xué)，大連市水務(wù)局.大連市水資源合理配置及初始水權(quán)問(wèn)題研究技術(shù)報(bào)告[R].大連:大連理工大學(xué)，大連市水務(wù)局，2007

[15]陳守煜，王子茹.基于對(duì)立統(tǒng)一質(zhì)量互變定理的水資源系統(tǒng)可變模糊評(píng)價(jià)新方法[J].水利學(xué)報(bào)，2011，42(3):253－262