水资源短缺风险综合评价摘要本文首先对北京市的2001-2009年得水资源短缺状况进行了调查，在综合考虑系统属性等风险过程后，利用层次分析模型对北京市各缺水影响因子的权重进行了定量分析，并基于致险因子承险因子及损害程度等影响因子构建了水资源系统风险的评价指标及模型：该指标体系由4层次共20个指标构成，能更好的表征风险的产生和构成；该模型包括参数计算与风险分级，能简便计算风险级别的划分。
其次，本文在综合考虑水资源呈现能力后得出结论，北京市能应对水资源系统风险，但是仍受约束性风险限制，可通过开源节流，调整产业结构及规划水资源管理来应付.然后我们对北京市2001-2009年水资源总量，地表水资源以及地下水资源量进行了调查。
运用Matlab处理系统对历年降水量进行了拟合，用origin处理系统对万元GDP水耗做出了拟合由此得出了缺水量波动性较大的结论。
最后本文对所建模型进行了升级及完善，采用灰色模型的建立改进方法，通过对无偏GM(1,1)模型的求解，得出2010年和2011年度致险率(RBI)、承险率(RSI)、脆弱性(CI)以及风险(ωD Ri)、风险损失(DI)的值，并由此得出结论，北京市未来两年的水资源短缺风险分别为28.40%和30.50%，正在呈上升趋势，逐年增高，不过基本上还在约束性风险级别内，为此建议管理机构还是要约束水资源使用来防范风险，通过推荐高效水资源系统管理，促进水资源优化配置进程等途径来促进水资源系统恢复，有效地减弱风险发生及潜在损害。
关键词：层次分析法水资源短缺风险多元回归拟合无偏GM（1,1）1.问题的重述水资源，是指可供人类直接利用，能够不断更新的天然水体。
主要包括陆地上的地表水和地下水。
近年来，我国、特别是北方地区水资源短缺问题日趋严重，水资源成为焦点话题。
水资源系统风险是由于天然来水的波动、地下水持续保障能力不足、供水条件落后以及水资源社会经济承载负担过重等因素综合作用的结果, 对社会、经济、环境存在潜在损害。
目前北京是世界上水资源严重缺乏的大都市之一，其人均水资源占有量不足300m3，为全国人均的1/8，世界人均的1/30，属重度缺水地区，北京市水资源短缺已经成为影响和制约首都社会和经济发展的主要因素。
政府采取了一系列措施, 如南水北调工程建设, 建立污水处理厂,产业结构调整等。
但是，气候变化和经济社会不断发展，水资源短缺风险始终存在。
如何对水资源风险的主要因子进行识别，对风险造成的危害等级进行划分，对不同风险因子采取相应的有效措施规避风险或减少其造成的危害，这对社会经济的稳定、可持续发展战略的实施具有重要的意义。
根据《北京2009统计年鉴》及市政统计资料提供了北京市水资源的有关信息。
利用这些资料和我们自己获得的其他资料，讨论得知北京市水资源短缺现在主要的主要原因是水资源供应小于水资源需求的矛盾，而如果想要解决这个问题，就必须要从影响因子来着手，所以问题细分为：北京市水资源短缺的主要因子是什么？各因子对于风险程度的贡献是多少？北京市已经到什么样的风险程度了？针对于主要因子我们怎样应对才能降低风险，从而做到有效调控？北京是未来几年又将面临什么程度的水资源短缺风险，又该如何应对呢？最后，通过建模等一系列过程进行分析检验并得出结论，且向北京市主管部门写一份建议报告。
2.问题的分析北京水资源人均占有量在世界各国首都中排名百位之后。
自上世纪七十年代以来，随着人口的大量增加和经济的发展，缺水成为北京面临的严重问题之一，近几年每年缺水均在4亿立方米左右。
地下水资源开采量逐年剧增，尽管目前对地下水开采进行了限制，地下水位有所上升，但仍处于超采状态。
地下水的超采会形成漏斗区，到目前为止，已经形成以朝阳区为中心，西到石景山、东至顺义、南至南苑、北到昌平约1600平方公里的漏斗区，引起地面沉降。
由于水位不断下降，致使井越打越深，形成恶性循环。
为应对这种情况，多年来，北京通过各种方式保证供水安全，除了通过调整产业结构加大节水力度，多次提高水价，强力推行农业、工业和城市节水，关、停、转移高耗水企业外，还独创了地表水、地下水、再生水、过境水、雨洪水和外调水的六水联调模式，对水资源进行合理调配，以此提高城市的供水能力。
尽管方法尽施，可是但是究竟是什么导致了北京市的如此现状呢，有没有什么原因呢，该怎样解决呢，经过本小组成员的查找资料和激烈的讨论，我们认为从水资源系统结构来看, 风险来源于系统属性和过程对潜在危害的抵抗乏力。
系统本身的输入主体短缺、过程波动及输出脆弱程度是导致系统风险产生的重要原因, 他们是水资源系统风险的致险因子; 水资源系统对致险因子进行反馈, 引导系统对潜在风险进行抵抗从而削减风险产生及危害, 我们将这一种反馈及抵抗性质称为承险能力; 系统在致险与承险因子相互作用下, 当致险压力大于承险能力时, 风险就产生了。
所以风险因子分为致险因子和承险因子，然后致险因子和承险因子又会细分为很多条。
不同的因子给风险所能带来的贡献的大小是不一样的，这个可以通过建模的方式利用Excel 或者Matlab ，origin 软件进行运算得出，之后通过制定一个风险等级将北京市的现状表现出来，并对作出巨大贡献的因子进行合理的调控。
至于北京市水资源短缺未来两年的预测，我们可以用回归的思想，再用灰色理论进一步改进，即可对北京市进行预测了。
3. 模型的假设与符号说明假设1：收集的北京市水资源各个数据都实际数值相差不大； 假设2：各个影响因素不会因突发事件发生突变； 假设3：建模收集数据真实可靠；假设4：建模中涉及主观分析的结论基本与事实相符；aij ——1—9标度理论得出的第i 项较第j 项的相对重要值 μij ——测度判断值ωD ——准则层D 下的相对权向量 CI —— 脆弱性ωSi ——系统风险的发生及传递对系统损害率在相关评价指标体系上的指标权重的重分配值Pi——风险潜在发生概率 Ri——系统风险因子的指标值 DI——风险损失 RBI ——致险率 RSI——承险率 ωD R i ——风险 b——回归系数；bint——回归系数的区间估计； r——残差；rint ——置信区间；stats——用于检验回归模型的统计量，有三个数值：相关系数r2、F 值、与F 对应的概率p 。
其中相关系数r2越接近1，说明回归方程越显著；F > F1-α（k ，n-k-1）时拒绝H0，F 越大，说明回归方程越显著；与F 对应的概率p α<时拒绝H0，回归模型成立。
α——显著性水平，一般为0.05或0.01，本论文中为0.05.4. 模型建立与求解4.1 针对于问题1：我们对北京市全年水资源总量与全年供水(用水)总量进行了调查和比较如表1所示。
表1[7]-[9]利用Origin软件处理的三位直观图如图1。
图1如图所示北京市全年供水总量均高于水资源总量，由于2008年北京市举办奥运会，为了保证北京市水资源能够及时供应，所以产生了供水量与水资源总量相差不大的局面。
其余年份数据反映北京市各年均呈缺水事态。
我们通过具体的系统属性，系统潜在损害指标来对系统表征层进行表征。
通过对水文水资源循环机理的研究，充分考虑水资源系统风险产生和传递机制，经过层层挑选，从而得到一下20个评价指标(如图2)[1]。
图2 水资源系统风险评价指标体系图说明： Ａ 代表目标层，即水资源系统风险评价要达到的目标， 即水资源系统风险求算。
Ｂ 代表风险属性层，即表征潜在危害排除系统承险能力而产生和传递的过程。
风险属性层有短缺性、波动性、脆弱性及承险性。
系统承险能力通过水资源系统的承险性来表征． Ｃ 代表风险属性系统表征层，即水资源系统面临风险的综合表征层， 但不能明确地表征系统的风险， 因此需要通过系统潜在风险损害可能及损害程度来细化表征。
本指标系统将风险属性层通过１２个指标来体现。
Ｄ 代表评价指标层，即通过具体的系统属性、系统潜在损害指标来对系统表征层进行表征。
通过对水文水资源循环机理的研究， 充分考虑水资源系统风险产生和传递机制， 经过层层挑选， 从而得到以下２０个评价指标。
这些指标不但从水量、水质方面对水资源系统风险进行了表征， 还兼顾了社会、经济、生态环境的效应， 并综合考虑了人对风险的积极适应以及人的应急性处理等主观能动的作用。
我们认为从水资源系统结构来看, 风险来源于系统属性和过程对潜在危害的抵抗乏力。
Kaplans 等从定量角度对风险进行了定义[2]，系统本身的输入主体短缺、过程波动及输出脆弱程度是导致系统风险产生的重要原因, 他们是水资源系统风险的致险因子; 水资源系统对致险因子进行反馈, 引导系统对潜在风险进行抵抗从而削减风险产生及危害, 我们将这一种反馈及抵抗性质称为承险能力; 系统在致险与承险因子相互作用下, 当致险压力大于承险能力时, 风险就产生了。
水资源系统风险的要素还包括损害程度, 致险压力、承险能力及损害程度，综合作用下的风险过程, 如图3所示[1]。
图3水资源系统风险过程及属性表征系统风险因子(如图2所示)可归结为致险因子和承险因子, 前者是指引起系统变化的因子, 包括系统结构的变化和外界干扰, 致使风险发生的概率为致险率; 而后者是指系统充分反馈或在历史事故后自我调节、自我适应而达到的能够应对危害的要素, 系统对风险削减能力为承险率。
下面对致险因子以及承险因子进行展开讨论。
水资源系统致险因子：①短缺性: 指水资源系统在自身运行过程中输入主体容易受到损害的性质,表征系统输入主体抵抗风险的不完备性。
短缺性体现在系统运行的供需不满足性以及系统已经受到损害的程度。
具体来看, 水资源系统的短缺性体现在使用短缺性、蓄水短缺性和环境短缺性三个方面, 即水资源缺水率、地下水超采和水体污染造成损害性。
②波动性: 相对于水资源系统多年正常运作的稳定程度, 波动性是指水资源系统因为系统波动或要素波动造成的系统不平衡运作的性质。
系统波动性来源于系统多年的不平衡性和系统输入输出的变动性。
波动可以用平均状态和极值差异来表征, 同时系统输入输出的稳定性也非常重要, 因此, 水资源系统波动性由多年波动、极值波动和水源波动来表征。
③脆弱性: 表征系统面临风险的潜在损害度, 即系统潜在输出抵抗风险的脆弱程度。
脆弱性指标能够衡量因风险产生而引起的损害程度, 体现在水资源- 社会经济耦合系统中, 主要是引起社会、经济、生态应对风险能力的下降。
社会损失体现在人均潜在利用能力的损失, 经济损失表征对生产活动造成的损害, 而生态损害体现在生态环境的破坏上。
水资源系统承险因子：水资源系统本身是一个动态的开放系统, 能通过自身的反馈调节来应对风险, 系统本身形成了一套承担风险发生和阻止风险破坏的承险因子体系。
水资源系统本身的资源禀赋、系统内部对风险事件形成的适应性以及风险发生时通过人为调度的应急性, 都是系统应对风险的有力保障, 是水资源系统承险因子。
