第9卷第4期 2011年8月 南水北调与水利科技 South-to-North Water Diversion and Water Science&Technology Vo1.9 No.4 Aug.2011 6oi:10.3724/SP.d.1201.2011.04102 河流水量分 配系 统的研发与应用 岑栋浩,邵东国,肖 淳,汪 奎,邹一雄 (武汉大学水资源与水电工程科学国家重点实验室,武汉430072) 摘要:目前,水量分配的内涵、原理及方法仍有待完善,现整合水量分配理论,以不同河流普遍适用为目标,通过VB、 ArcEngine及Access技术集成可视化的河流水量分配系统。河流水量分配系统内嵌水量分配指标库及方法库,融 入了层次分析法(AHP)、多目标半结构性模糊优选模型法及基于实数编码加速遗传算法的投影寻踪分类模型 (RAGA_PPC)。系统在府环河流域得到了有效的应用,基于指标库构建分水指标体系,分别使用AHP和RAGA- PPC方法确定指标的权重及行政区的分水比例,分水结果显示,通过本系统AHP和RAGA-PPC方法得到的府环 河流域分水方案与其需水预测结果大致吻合。 关键词:用水总量控制;ArcEngine;VB;府环河流域;层次分析法;投影寻踪分类模型 中图分类号:TV212.4;TP319 文献标识码:A 文章编号:1672-1683(2011)04—0102—03 Development and Application of River Water Allocation System CEN Dong-hao,SHAO Dong-guo,XIAO Chun,WANG Kui,ZOU Yi—xiong (State Key Laboratory of Water Resources&Hydropower Engineering Science, Wuhan University,Wuhan 430072,China) Abstract:Conflict of water supply and demand is acute in our country.The strictest management system of water resources has been built.Ifs important tO allocate water in a river to achieve tota1 quantity control of water to be consumed.At present,the connotation,principle and method are not perfect enough.This article summarized the theory of water allocation,developed a system for water allocation in rivers using technology of VB6.0,ArcEngine and Access with the goal of being suitable for differ— ent rivers.This system contained a series of indices and methods for water allocation.The methods included the analytic hierar— chy process method(AHP),the multi—objective and semi-structural fuzzy optimization model method,and the projection pursuit classification model(PPC),which based on real coding based accelerating genetic algorithm(RAGA).The system had an effec— tive application in Fuhuan River Basin.The indices system for water allocation was built based on indices contained in the sys— tern.AHP and RAGA-PPC had been used in calculating water proportion among different regions.The results showed that the water proportion gotten by AHP and RAGA-PPC was close to the forecast results of water need. Key words:total quantity control of water consumed;ArcEngine;VB;River Basin;AHP;PPC 水量分配是对水资源可利用量或可分配水量向行政区 域进行逐级分配的过程,我国在河流水量分配领域中开展了 不少研究与实践[1 ]。针对水利部制定的“用水总量控制、用 水效率控制以及河流纳污总量控制”三条红线,河流水量分 配作为实现用水总量控制的重要措施,迫切需要在大范 围内实践。 法¨3],进行最优解模拟,代表系统包括MIKEBASIN、Water— ware和WMS等。这些系统功能强大,但使用过程也很复 杂,基于此,提出具有水量分配功能、操作简单、运行可靠的 河流水量分配系统。河流水量分配系统是基于指标体系法 原理,构建水量分配指标体系,计算权重进而得到不同地区 的分水比例。 目前,水资源优化配置系统的应用已很成熟,但缺乏可 视化的河流水量分配系统。水资源优化配置系统与河流水 量分配系统差异很大,水资源优化配置系统应用系统优化方 本文研发的河流水量分配系统内嵌3个子模块:①分水 总量计算模块;②水量分配计算模块,内嵌指标库和计算方 法库,方法库融入层次分析法(简称AHP) 】、多目标半结构 收稿日期:2011—05—13 修回日期:2011—06—17 网络出版时间:2011 07 23 网络出版地址:http://www.cnki.net/kcms/detail/13.1334.TV.20110723.1707.033.html 基金项目:国家水体污染控制与治理科技重大专项(2011ZX07213—5) 作者简介:岑栋浩(1987一),男,广东湛江人,硕士,主要从事水量水质联合及系统开发等方面的研究。E-mail:cendonghao@163.com ・102・ 爱| |焉 菇 岑栋浩等・河流水量分配系统的研发与应用 性模糊优选模型法[5 以及基于实数编码加速遗传算法的投 影寻踪分类模型(简称RAGA-PPC)E7];③水量分配结果显 存水量分配的年度分水方案计算结果。④月份分水表。储 存水量分配的月份分水方案计算结果。 示模块。为验证系统的科学性和有效性,本文以湖北省府环 河流域为例,计算府环河流域的水量分配方案,分析表明,水 量分配结果与现实情况大致吻合。 2水量分配子模块介绍 2.1 河流水量分配原理 河流水量分配是以河流的可利用水量或可分配水量作 为分配对象,结合河流的供需水情况,再综合考虑影响分水 l河流水量分配系统的基本架构 河流水量分配系统以VB 6.0为开发语言,使用Ar— cEngind。 二次开发工具包实现ArcGIS提供的空间查询与分析、 可视化表达等功能,并采用Access 2003作为后台管理数据库。 1.1 系统结构与功能设计 河流水量分配包括分水总量计算子模块、水量分配子模 块以及展示子模块,与系统各项功能一一对应,系统的模块 结构情况见图1。 …………… … …………星萤 蔓 图1 河流水量分配系统模块结构 Fig.1 Model Structure of River Water Allocation System 1.2数据库设计 河流水量分配系统的数据主要涉及地形地貌、行政区域 等基础地理数据和河流水系、降水、河川径流、地下水位、供 需水等水文数据以及水利工程、社会经济等其他数据。为了 便于数据管理,系统以Access数据库作为后台管理数据库, 主要包含4种类型数据表,见图2。①基本属性表。储存 图2数据库框架 Fig.2 Database Structure 流域内各个地区的基本信息(GDP、人口、土地面积等)。② 指标值表。储存各个地区所有指标值。③年度分水表。储 结果的因素,构建水量分配指标体系。目前,国内外有关研 究中普遍采用的水量分配方法,主要分为3类:①应用层次 分析方法和模糊评价方法;②应用群决策模型;③应用德尔 菲法(又称专家决策法)。具体包括综合权重法、层次分析法 (AHP)、模糊物元分析法、灰色关联度法、多目标半结构性模 糊优选模型法等以及最近出现的一种非常新颖的方法—— 基于实数编码加速遗传算法(RAGA)的投影寻踪分类(PPC) 模型。选择其中一种方法,确定指标体系中所有指标的权 重,最终得到不同行政区的分水比例。 2.2水量分配子模块概况 2.2.1指标库操作 不同流域社会经济发展情况及产业结构、用水构成的差 异,导致水量分配指标选择的侧重点各不相同。系统内嵌的 指标库提供常用指标选择及自定义指标输入两种模式,同时 归纳为公平性、有效性及可持续性3大类指标,形成一个流 域普遍适用、综合全面的指标体系构建平台。内嵌的公平性 指标包括人口总数、缺水率、耕地面积、人均用水量等14个 常用指标;有效性指标包括GDP总量、人均GDP、单方水 GDP、工业用水重复利用率等14个常用指标;可持续性指标 包括单位GDP排污量、水质达标率、人口增长率等6个常用 指标。 2.2.2方法库操作 系统内嵌的方法库包含3种常见的水量分配方法: AHP、多目标半结构性模糊优选模型法及RAGA-PPC。 AHP及多目标半结构性模糊优选模型法的原理请参考文献 [4—5],下面重点阐述RAGA-PPC的原理。 2.3 RAGA—PPC的原理 模型的建模过程如下。 步骤1:分水指标集的归一化处理。 步骤2:构造投影指标函数。PPC方法就是把P维数据 {,27 ( , )l 一1,2,…,P) 宗合成以口一{。(1),a(2),口(3),…, n( ))为投影方向的一维投影值Z( )。 z( )一 口( ) ( , )( 1,2,…, ) (1) J一1 综合投影指标值时,要求投影值Z( )的散布特征为局部 投影点尽可能密集,最好凝聚成若干个点团,而在整体上投 影点团之间尽可能散开。因此,投影指标函数可以表达成 Q(&)一SzDz (2) 式中:Sz一投影值Z( )的标准差;De一投影值Z( )的局部密 度,即 /妻[ )一E(z)] Sz ^√——— ,i=i z(广一 (3) ∞镑 罄| % 弹簪 1% g辨 生卷与力:绣 ・103・ Ⅲ 第9卷总第55期・南水北调与水利科技・2011年第4期 Dz一∑∑(R—r( , ))・“(R—r(i, )) (4) 式中:E(Z)一序列{Z( )l 一1,2,…, }的平均值;R一局部密 度的窗口半径,R可以根据试验确定,一般取值0.1 Sz; r(i, )一样本之间的距离,r(i, )一lZ( )一Z( )I;“(£)一单 位阶跃函数,当£≥O时,其值为1;当t<0时,其值为0。 步骤3:优化投影指标函数。 向。,即 . 通过求解投影指标函数最大化问题来估计最佳投影方 ; 清选择标题:磴 …标题 ……… —囊 I肯选掸蚓形:瞵i … —: ~—一l 曳J 墅璧夔剧 ;最大化目标函数 maxQ(a)一Sz・Dz (5) 点堕£萋到标题化置盘l 捌 ; x轴.F————~ Y轴: 芳 一—一 捌剜 约束条件 0 ∑ ( )一1 (6) 这是一个以{a( )l 一1,2,…,P}为优化变量的复杂非 线性优化问题,用传统的优化方法处理比较困难。因此,本 文应用基于实数编码的加速遗传算法(RAGA)来解决其高 维全局寻优问题。 步骤4:RAGA求解。RAGA的计算原理请参考文献E9]。 3 实例应用 府环河是湖北省内仅次于汉江、清江的第3大水系,发源 于大洪山北麓,随州岗镇双门洞,经随州、广水、安陆、孝 感、黄陂等县市于武汉市谌家矶入长江,是长江的一级支流。 干流全长348 km,集水面积14 455 km2(包括环水上游河南省 内面积),省内集水面积14 287 ,支流348条。下面以府 环河流域为研究对象,应用流域水量分配系统,已求得来水 频率P一50%,水平年2020年情况下,府环河流域可分配水 量为248 329万m3,以下是府环河流域水量分配计算过程。 从指标库选取8个公平性指标,11个有效性指标,以及 6个可持续性指标,构建府环河流域水量分配指标体系,图3 为水量分配指标的选择界面。 图3选择水量分配指标 Fig.3 Selection for indices of water allocation 从方法库选择AHP、RAGA-PPC两种水量分配方法, 计算得到2020年(水平年),P=50 频率下的府环河流域水 量分配结果,如图4所示。 比较府环河流域2020年用水预测值和年度分水结果, 统计得到如表1所示的情况。总的来说,对随州市和孝感市 而言,使用AHP和PPC方法的得到的分水结果大于2020 ・104・0咚 图4府环河流域水量分配结果(图形显示方式) Fig.4 Water allocation results in Fuhuan River Basin 年的用水预测值;武汉市的情况刚好相反。对AHP方法而 言,判断矩阵的确定会影响最终的分水结果;在PPC方法的 计算过程中,交叉概率、变异概率及加速次数的确定也会影 响最终的分水结果。 表1水量分配结果 Table 1 Results of Water Allocation 城市 。 揪 AH 箬 果即 果 随州市8O 741 82 445 83 190 孝感市 100 079 104 298 101 3l8 武汉市 67 509 61 586 63 821 4结论 流域水量分配是实现用水总量控制的重要措施之 一,针对我国用水矛盾突出的现状,基于GIS平台的水量分 配系统研发显得十分重要。本文基于VB 6.0、ArcEngine及 Access研发了流域水量分配系统,系统具备分水总量计算、 水量分配及可视化结果显示功能,为水量分配提供一套较完 整的解决方案,并在府环河流域的水量分配过程中得到验 证。系统内嵌水量分配指标库及方法库,可供用户灵活选 择,旨在为不同流域的水量分配构建一个通用平台。 参考文献(References): [15何俊仕,李秀明,尉成海,等.大凌河流域水量分配方法研究 _J].人民黄河,2008,30(4):50—52.(HE Jun shi,LI Xiu—ming, YU Cheng—haiet al,Research of Method on Water Allocation in DalingRiver[J].YellowRiver,2008,30(4):5O 52.(inChinese)) [23赵鸣雁,陈吉宁,程春田,等.水资源配置多准则决策敏感性分 析[J].水电能源科学,2009,(2):35 36.(ZHAO Ming—yan, CHEN Ji—ning,CHENG Chun-tian,et a1.Sensitivity Analysis of Multi Criteria Decision for Water Resources Allocation[J].Wa— ter Resources and Power,2009,(2):35—36.(in Chinese)) [3]陈良程,蔡治国,王晓霖,等.黄河流域水量调度方案编制方法 研究EJ].应用基础与工程科学学报,2003,(2):96—103.(CHEN Liang—cheng,CAI Zhi—guo,WANG Xiao-lin.Methods of Consti— tuting the Yellow River Water Resources Regulation Scheme _J]Journal of Basic Science and Engineering,2003,(2):96— 103.(in Chinese)) (下转第133页) 张宁霞等・松嫩平原碳酸盐渍土的工程特性 有实验结果表明,碳酸盐渍土的抗剪特性与其干密度和 含水量有一定的关系:含水量相同时,干密度越大,黏聚力和 内摩擦角也增大。当干密度相同时,含水量增大,黏聚力和 内摩擦角减小。这是因为含盐量一定时,其抗剪强度与土中 参考文献(References) [1]罗伟甫.盐渍土地区公路工程[M].人民交通出版社,1980. (I UO Wei—fu.Saline Soil Area Project M 1.China Communica— tions Press。l980.(in Chinese)) 含液量关系较大。当含液量增加时,土中的结晶盐被溶解, 同时,颗粒问的连接关系由水胶连结转变为水连结,使土体 之间的联结力下降;含液量降低时,土中液态的盐转变为固 态,在土中起胶结作用,同时,颗粒粒问水连结转变为水胶连 [2]罗金明,邓伟,张晓平.从小尺度探讨苏打盐渍土的特征以及演 变规律l_I1].干旱区资源与环境,2009,23(2):138 140.(LUO Jin ming.DENG Wei,ZHANG Xiao ping.Properties and Evolvement of Carbonated Saline Soil in MicroScale【J I.Journal of Arid I.and Resources and Environment,2009,23(2):1 38— 140.(in Chinese)) 结,抗剪强度明显增加[5_6。 据其他有关盐渍土研究结果显示:在气候干燥的状态下, 苏打型盐渍土与常规土相比,更紧密坚硬,强度较高,力学性 能较好。在吸水试验过程中随着含水量的增加,土体的抗剪 强度迅速降低。并且与普通盐渍土相比,苏打型盐渍土抗剪 强度随含水量增加而降低的更为显著。因此,苏打型盐渍土 在浸水(雨、雪)条件下,吸湿性强,稳定性低,并且不易排水, [3]孙婉.松嫩平原苏打型盐渍土强度挣l生研究[I)].北京:中国地质大学, 2006.(SUNWan.Researchy onStrength ofSongNenPlainSoda划ine SiolED].Beijing:China University of Geoscience.s,2006.(in Chinese)) [4]张国辉,于青春,魏国强.松嫩平原碳酸盐渍土的基本性质研究 lJ].水文地质工程地质.2007(2):37 39.(ZHANC-Guo hui, YU Qing chun,WEI Guo—qiang.Basic Property of Songnen 很难干燥,力学性能迅速恶化,造成路基的严重塌陷 。 Plain Saline alkali Soil[J ̄.Hydrogeology and Engineering Ge ology,2007,(2):37 39.(in Chinese)) 3结论 通过人工配制不同含盐量的盐渍土,测得随着含盐量的 增大碳酸盐渍土的液塑限都有所提高。通过碳酸盐渍土洗 盐前后的颗粒分析试验结果可知:洗盐后,细粒径土颗粒所 占百分比有很大程度的提高,这显示出,盐渍土中的盐的胶 结作用极大的影响了其颗粒组成。对于碳酸盐渍土的击实 特性的讨论,可以看出,含盐量越大,则需要较大的含水量才 能达到最佳压实状态。通过膨胀试验来研究遇水膨胀特性, [5]陈炜韬,李姝,王鹰.含盐量、含盐类别对盐渍土抗剪强度的影 响_J].铁道建筑技术,2005,(6):54 57.(CHEN Wei tao,I I Shu,WANG Ying.The Effect of Salt Content and Salt Type on the Shear Strength of Saline Soil[J].Railway Construction Technology,2005,(6):54 57.(in Chinese)) [6]苑杰,陈晓飞.棕壤盐水入渗水分运移特征研究[J].水电能源 科学,2010,29(8):104 106,168.(YUAN Jie,CHEN Xiao fei. Features of Infiltration Water Migration in Brown Soil with Sa line Water[J].Water Resources and Power,201 0,29(8):104 lO6。168.(in Chinese)) 结果发现:当碳酸盐含量大于0.5 时,膨胀量显著增大,当 含盐量从0.5 增加到2.0 时,膨胀率最大可以达到23 。 [7]梁健伟,房营光,陈松.含盐量对极细颗粒黏土强度影响的试验 研究[J].岩土力学与工程学报,2009,28(2):3823 3827.(Ex— perimental Research On Effect()f Salt Content On Strength of 通过对盐渍土抗剪强度的研究发现:即土的含盐量增加时, 黏聚力有减小的趋势,内摩擦角先减小后增大。 Tiny-Partical Caly[J].Chinese Journal of Rock Mechanics and Engineering,2009,28(2):3823 3827.(in Chinese)) (上接第104页) [4] 袁伟,郭宗楼,袁华.层次方在黑河水量调度决策支持系统 中的应用[J].灌溉排水学报,2005,24(2):58—61.(YUAN Wei, GUO Zong lou,YUAN Hua.The Application of Analytical Hi erarchy Process in Decision Support System Used in Hei River 45.(FU Qiang,JIN Ju—liang,HANG Chuan.Application of Projection Pursuit Model to Optimize Paddy Irrigation Schedule _J].Journal of Hydraulic Engineering,2002,(10):39—45.(in Chinese)) Water Dispatch[J].Journal of Irrigation and Drainage,2005,24 (2):58 61.(in Chinese)) [8]朱仕杰,南卓铜.基于ArcEngine的GIS软件框架建设[J].遥 感技术与应用,2006,21(4):385—390.(ZHU Shijie,NAN [5] 陈守煜,丁晓阳.多目标半结构性决策模糊识别理论在北水南 调中的应用[J].水资源保护,2003,(1):29 33.(CHEN Shou— yu,DING Xiao yang.Application of Multi—objective and Semi Structural Decision making Fuzzy Recognition Theory to The Zhuo tong.Building GIS Framework with ArcEngine[J].Re— mote Sensing Technology and Application,2006,21(4):385 390.(in Chinese)) [9]金菊良,杨晓华,丁晶.基于实数编码的加速遗传算法[J].四川 大学学报(工程科学版),2000,32(4):20 24.(JIN Ju—liang, YANG Xiao—hua,DIN(}J ing.Real Coding Based Acceleration North to South Water Diversion Project lJ.Water Resources Protection,2003,(1):29—33.(in Chinese)) [6] 陈守煜.多目标系统半结构性决策理论与方法[J].辽宁工程技 术大学学报(自然科学版),2001,20(5):578 582.(CHEN Shou yu.Semi structural Decision-making Theory and ApproachFor Genetic Algorithm[J].Journal of Sichuan University(Engineer Science Edition,2000,32(4):20 24.(in Chinese)) [1O]宋关福,钟耳顺.组件式地理信息系统研究与开发_J].中国图 像图形学报,1998,3(4):313 317.(SONG Guan fu,ZHONG Er shun.Research and Development of Components Geo Multi一()bjetive System[J].Journal of Liaoning Technical Uni versity(Natural Science),2001,20(5):578 582.(in Chinese)) [7] 付强,金菊良,粱川.基于实码加速遗传算法的投影寻踪分类模 型在水稻灌溉制度优化中的应用[J].水利学报,2002,(10):39一 graphic Information Systems[J].Journal of Image and Graph— ics,1998,3(4):313-317.(in Chinese)) 水文地质与工程地质 ・133・