水资源短缺风险综合评价模型_数学建模竞赛参赛论文.doc

封一答卷编号（参赛学校填写）：THSY06答卷编号（竞赛组委会填写）：论文题目：B水资源短缺风险综合评价模型 组 别：本科生参赛队员信息(必填)： 姓 名专业班级及学号联系电话参赛队员1数学与应用数学200906010242参赛队员2数学与应用数学200906010220参赛队员3数学与应用数学200906010145参赛学校：报名序号：（可以不填）30封二答卷编号（参赛学校填写）：THSY06答卷编号（竞赛组委会填写）：评阅情况（学校评阅专家填写）：学校评阅1.学校评阅2.学校评阅3. 评阅情况（省赛评阅专家填写）：省赛评阅1.省赛评阅2.省赛评阅水资源短缺风险综合评价模型摘要 本文通过对水资源系统影响因素的分析，确立了水资源短缺的主要风险因子：风险率、脆弱性、重现期、可恢复性、风险度，并以这些因子为衡量指标建立了数学模型对北京市水资源短缺风险进行了综合评价。在问题1中，建立了衡量指标模型、概率模型，运用方差、标准差对风险因子进行了检验。在问题2中，综合问题1的风险因子建立了一个描述北京市水资源短缺风险程度的模糊综合评价模型，在模型建立过程中运用了分类、划分方法将水资源短缺风险因子分为5个等级，并用层次分析法确定了风险因子的权系数，构造隶属函数确立了水资源短缺风险综合评价模型，最终划分了水资源短缺风险的5个等级。在问题3中，运用图形分析法，根据北京市水资源情况图像及问题2中建立的水资源短缺风险综合评价模型对北京市未来两年的水资源短缺风险进行了预测，预测结果表明：尽管加大再生水利用量、南水北调工程在一定程度上缓解了北京市水资源短缺的紧张局面，但未来两年北京市水资源短缺风险仍将处于高风险水平。问题4则对北京市水行政主管部门提出了解决水资源短缺的可行性方案。在整个模型建立过程中还用到了EXCEL、MATLAB等工具。关键词 北京、水资源短缺风险、模糊数学、概率、风险因子问题重述改革开放以来，我国工业和城镇生活用水持续增长，由于气候条件、水利工程设施、工业污染、人口规模等因素，近年来我国北方地区水资源短缺问题日益严重，而北京人均水资源占有量不足300立方米，为全国人均的1/8，世界人均的1/30，已成为世界上水资源严重缺乏的大都市之一。现根据北京2009统计年鉴及市政统计资料提供的有关1979年至2000年北京市水资源短缺的状况表（见附录）讨论如下问题：1、北京市水资源短缺风险的主要风险因子。（影响水资源的因素，如：气候条件、水利工程设施、工业污染、农业用水、管理制度，人口规模等。）2、对北京市水资源短缺风险进行综合评价并根据综合评价标准作出风险等级划分，并考虑如何对主要风险因子进行调控才能使得风险降低？3、对北京市未来两年水资源的短缺风险进行预测，并提出应对措施。4、以北京市水行政主管部门为报告对象，写一份建议报告。问题分析 近年来，受气候变化和经济社会不断发展的影响，水资源短缺问题日趋严重，对水资源短缺风险的研究已引起了广泛的重视。要对北京市水资源短缺风险进行综合评价，首先要明确北京市水资源短缺风险的影响因子及其影响程度的大小；其次，要建立模型考察各风险因子对水资源短缺风险的具体影响。在影响北京市水资源短缺风险的各个主要因子中又有其它因素对其影响，在这里我们忽略了影响较小的因素。判断北京市水资源短缺风险问题要综合考虑，既要考虑各因素对风险因子的影响，又要考虑各风险因子对水资源短缺风险的影响，然后再对其进行综合评价。根据此题特征，我们很容易就想到了模糊概率的相关知识，通过计算概率来实现风险的综合评价。为此，先构造各风险因子的评价指标，然后，再根据各风险因子来构造水资源短缺风险的评价函数。最后，根据模型函数进行风险等级划分并提出相应的解决方案。符号说明符号说明S水资源失事状态Z水资源正常状态M总用水量G水资源总量Xt水资源系统状态变量It水资源系统状态量ui（i=1、2、3、4、5）分别表示衡量风险率、脆弱性、重现期、可恢复性、风险度的指标T水资源系统工作总历时T1水资源系统失事总时间QZt第t年水资源短缺量SZt第t年水资源总量X失事状态下缺水量年均值Y水资源总量年均值Zt水资源可恢复状态量THF水资源可恢复状态总量TSX水资源状态总量D(x)衡量风险度的方差E(x)衡量风险度的标准差d(q,n)第n与n+1 次失事时间间隔Vi(i=1,2,3,4,5)水资源短缺风险等级U、V定义的两个有限域i（i=1,2,3,4,5）各风险的权系数bi（i=1,2,3,4,5）水资源短缺风险等级判别指标A、B分别为U、V的模糊子集RU变换关系矩阵Uvi()Ui对Vi的隶属函数模型假设1、 北京市人均年用水量一定；2、 北京市单位产业用水量一定；3、 不考虑水资源自身净化能力；4、 假设水资源短缺风险的各风险因子间影响程度可以忽略；5、 在计算过程中忽略微小数值对结果的影响；6、 忽略信息的有限性和不完全性对模型的影响。模型建立、求解与分析水资源系统是一个复杂的大系统，广泛存在着随机性和模糊性，由于随机性是因果率的破缺、模糊性是排中率的破缺，所以在水资源短缺风险评价模型的设计中应同时考虑这两种因素的影响。为了缓解北京市水资源短缺问题，建立数学模型对北京市水资源短缺风险进行综合评价。根据北京市水资源短缺风险的主要风险因子：风险率、脆弱性、重现期、可恢复性、风险度建立如下模型：问题1 水资源是人类赖以生存的能源。影响水资源的因素有很多，例如：气候条件、自然灾害、水利工程、工农业污染、污水利用率、管理制度、人口数量及规模、入境水量等，而这些因素都影响着水资源系统处于失事状态时所经历的时间、损失的严重程度、失事的次数、恢复到正常状态经历的时间等。综合以上原因可以将水资源的风险因子概括为：1、风险率 考虑到水资源系统的不确定性，如果供水系统处于失事状态， 且水资源系统的工作有长期的记录，风险率可以定义为水资源系统不能正常工作的时间与整个工作历时之比，即水资源系统风险率衡量指标为： 其中2、脆弱性脆弱性是指水资源系统处于失事状态的平均损失严重程度。为定量描述水资源系统的脆弱性，我们假设水资源系统第t年水资源短缺量为QZt水资源总量为SZt，T1为水资源系统失事总时间,则水资源系统脆弱性衡量指标为：上式中,失事状态下缺水量年均值X，水资源总量年均值Y为： 3、重现期重现期是指再次进入失事状态S时所经历的时间。则水资源短缺重现期衡量指标为：其中，重现期4、可恢复性可恢复性即水资源系统从失事状态到正常状态的可能性。系统的可恢复性越高，表明该系统从失事状态恢复到正常状态的能力越强，反之，则越弱。根据已知条件它可以用如下指标来衡量： 其中，5、风险度标准差反映了变量的稳定与波动、集中与离散的程度，标准差越小，波动越小，反之，则波动越大。利用标准差的特征可以定义以下指标来衡量风险度的大小：其中，方差在类似分析中，一般假定P1=P2=PT=1/(T-1) 标准差，问题2针对问题1中的各种风险因子，建立一个基于模糊综合评价的水资源短缺风险模型。如果北京市年水资源总量小于年需水量，即水资源处于短缺状态。则基于水资源系统的模糊不确定性，可以构造一个合适的模糊综合评价模型：设U=U1,U2,U3,U4,U5和V=V1,V2,V3,V4,V5为两个有限域，A= 1 ，2 ，3 ，4 ，5,0i1,B=b1,b2,b3,b4,b5分别为U、V的模糊子集,则模糊综合评价可表示为下列模糊变换： B=ARU由上述假设可知，关系矩阵RU与水资源短缺的主要风险因子有关。关系矩阵RU中的元素即为因素Ui对应等级Vi的隶属度，其值可根据各级评价因子的实际数值对照各因子的分级指标推求。下面我们将评语分为5个等级：表一：各风险因子指标分级水资源短缺风险等级风险率指标U1脆弱性指标U2可恢复性指标U3重现期指标U4风险度指标U5V1低0.200.200.809.000.20V2较低0.20-0.400.20-0.400.60-0.806.00-9.000.20-0.60V3中0.40-0.600.40-0.600.40-0.603.00-6.000.60-1.00V4较高0.60-0.800.60-0.800.20-0.401.00-3.001.00-2.00V5高0.800.800.201.002.00由于水资源风险率、脆弱性、风险度是越小越优性指标 ，而可恢复性、重现期是越大越优性指标，所以，对于U1、U2、U3、U4、U5各级评语构造如下隶属函数： 从而，得到RU的具体向量表示：利用层次分析法确定水资源短缺风险各因子的权重系数，即A= 1 ，2 ，3 ，4 ，5=1/2，1/4,1/8,1/12,1/24于是可得到如下综合评判向量：在综合评判中，我们可以将上述的向量转化为一般实数的加法，即最后，选取maxbj对应的评语作为水资源短缺风险的综合评价结果。 根据上述评价标准，我们可以将水资源短缺风险分为如下级别：表二：水资源短缺风险级别评价水资源短缺风险等级风险级别水资源短缺风险的特征V1(1b)低风险可以忽略的风险V2(0.9b1)较低风险可以接受的风险V3(0.8b0.9)中风险边缘风险V4(0.7b0.8)较高风险不可接受风险V5(b G;D= M G;D= M GM (C)=1;M (D)=0T1=0;for k= MT1= T1+kendclcT1M=42.92 50.54 48.11 47.22 47.56 40.05 31.71 36.55 30.95 42.43 44.64 41.12 42.03 46.43 45.22 45.87 44.88 40.01 40.32 40.43 41.71 40.4 38.9 34.6 35.8 34.6 34.5 34.3 34.8 35.1; %总用水量 G= 38.23 26 24 36.6 34.7 39.31 38 27.03 38.66 39.18 21.55 35.86 42.29 22.44 19.67 45.42 30.34 45.87 22.25 37.7 14.22 16.86 19.2 16.1 18.4 21.4 23.2 24.5 23.8 34.2; %水资源总量QZt1=M-GB=QZt10;QZt1(B)=0;clcQZt1M=42.92 50.54 48.11 47.22 47.56 40.05 31.71 36.55 30.95 42.43 44.64 41.12 42.03 46.43 45.22 45.87 44.88 40.01 40.32 40.43 41.71 40.4 38.9 34.6 35.8 34.6 34.5 34.3 34.8 35.1; %总用水量 G= 38.23 26 24 36.6 34.7 39.31 38 27.03 38.66 39.18 21.55 35.86 42.29 22.44 19.67 45.42 30.34 45.87 22.25 37.7 14.22 16.86 19.2 16.1 18.4 21.4 23.2 24.5 23.8 34.2; %水资源总量QZt2=G-MB=QZt2 G;D= M GM (C)=1;M (D)=0T1=0;for k= MT1= T1+kendclcT1It=M %水资源系统状态量for i=1:29dqnm(i)=It(i+1)-It(i)endclcdqnmdqnn=dqnm G;D= M G;D= M GM (C)=1;M (D)=0Tsx=0;for k= MTsx=Tsx+kendclcTsx %水资源状态总量if Tsx=0; U4=0;else U4=Thf/Tsx;endU4 %水资源可恢复性衡量指标运行结果：Tsx = 26U4 = 0.1538 %水资源可恢复状态量%plot(Zt,+)%水资源可恢复状态量程序5（计算水资源风险度衡量指标） clcclear M=42.92 50.54 48.11 47.22 47.56 40.05 31.71 36.55 30.95 42.43 44.64 41.12 42.03 46.43 45.22 45.87 44.88 40.01 40.32 40.43 41.71 40.4 38.9 34.6 35.8 34.6 34.5 34.3 34.8 35.1; %总用水量 G= 38.23 26 24 36.6 34.7 39.31 38 27.03 38.66 39.18 21.55 35.86 42.29 22.44 19.67 45.42 30.34 45.87 22.25 37.7 14.22 16.86 19.2 16.1 18.4 21.4 23.2 24.5 23.8 34.2; %水资源总量QZt1=M-GB=QZt10;QZt1(B)=0;clcQZt2=G-MB=QZt20;QZt2(B)=0;clcQZt=QZt1+QZt2 %第t年水资源缺水量SZt=G; %第t年水资源总量SZ=sum(SZt); %水资源总量T=30 %水资源系统工作的总历时Y=SZ/T %水资源总量年均值QZc=QZt-Y %第t年缺水量与水资源总量年均值的差值QZcc=QZc(1)2 QZc(2)2 QZc(3)2 QZc(4)2 QZc(5)2 QZc(6)2 QZc(7)2 QZc(8)2 QZc(9)2 QZc(10)2 QZc(11)2 QZc(12)2 QZc(13)2 QZc(14)2 QZc(15)2 QZc(16)2 QZc(17)2 QZc(18)2 QZc(19)2 QZc(20)2 QZc(21)2 QZc(22)2 QZc(23)2 QZc(24)2 QZc(25)2 QZc(26)2 QZc(27)2 QZc(28)2 QZc(29)2 QZc(30)2;QZz=sum(QZcc);Dx=QZz/(T-1) %衡量风险度的方差Ex=sqrt (Dx) %衡量风险度的标准差U5=Ex/Y %水资源风险度衡量指标运行结果：QZt = Columns 1 through 6 4.6900 24.5400 24.1100 10.6200 12.8600 0.7400 Columns 7 through 12 6.2900 9.5200 7.7100 3.2500 23.0900 5.2600 Columns 13 through 18 0.2600 23.9900 25.5500 0.4500 14.5400 5.8600 Columns 19 through 24 18.0700 2.7300 27.4900 23.5400 19.7000 18.5000 Columns 25 through 30 17.4000 13.2000 11.3000 9.8000 11.0000 0.9000T = 30Y = 29.2327QZc= Columns 1 through 6 -24.5427 -4.6927 -5.1227 -18.6127 -16.3727 -28.4927 Columns 7 through 12 -22.9427 -19.7127 -21.5227 -25.9827 -6.1427 -23.9727 Columns 13 through 18 -28.9727 -5.2427 -3.6827 -28.7827 -14.6927 -23.3727 Columns 19 through 24 -11.1627 -26.5027 -1.7427 -5.6927 -9.5327 -10.7327 Columns 25 through 30 -11.8327 -16.0327 -17.9327 -19.4327 -18.2327 -28.3327Dx = 362.5391Ex = 19.0405U5 = 0.6513程序6（计算水资源短缺风险级别）U1=0.8667; %水资源风险率衡量指标U2=0.4695; %水资源脆弱性衡量指标U3=0.1600; %水资源重现期衡量指标U4=0.1538; %水资源可恢复性衡量指标U5=0.6513; %水资源风险度衡量指标U=U1 U2 U3 U4 U5 %风险因子衡量指标的集合miuv1phi11=0.2miuv2phi11=0.3miuv3phi11=0.5miuv4phi11=0.7miuv5phi11=0.8 %风险率、脆弱性、重现期、可恢复性、风险度低风险的指标miuv1phi12=0.2miuv2phi12=0.3miuv3phi12=0.5miuv4phi12=0.7miuv5phi12=0.8 %风险率、脆弱性、重现期、可恢复性、风险度较低风险的指标miuv1phi13=0.8miuv2phi13=0.7miuv3phi13=0.5miuv4phi13=0.3miuv5phi13=0.2 %风险率、脆弱性、重现期、可恢复性、风险度中风险的指标miuv1phi14=9miuv2phi14=7.5miuv3phi14=4.5miuv4phi14=2miuv5phi14=1 %风险率、脆弱性、重现期、可恢复性、风险度较高风险的指标miuv1phi15=0.2miuv2phi15=0.4miuv3phi15=0.8miuv4phi15=1.5miuv5phi15=2 %风险率、脆弱性、重现期、可恢复性、风险度高风险的指标W1=miuv1phi11 miuv2phi11 miuv3phi11 miuv4phi11 miuv5phi11 %风险率U1衡量指标W2=miuv1phi12 miuv2phi12 miuv3phi12 miuv4phi12 miuv5phi12 %脆弱性U2衡量指标W3=miuv1phi13 miuv2phi13 miuv3phi13 miuv4phi13 miuv5phi13 %可恢复性U3衡量指标W4=miuv1phi14 miuv2phi14 miuv3phi14 miuv4phi14 miuv5phi14 %重现期U4衡量指标W5=miuv1phi15 miuv2phi15 miuv3phi15 miuv4phi15 miuv5phi15 %风险度U5衡量指标clcRu0=W1;W2;W3;W4;W5 %变换关系矩阵if U1=miuv1phi11miuv1phi11=1;elseif U1=miuv1phi11miuv2phi11=(miuv2phi11-U1)/(miuv2phi11-miuv1phi11);elseif U1=miuv2phi11 miuv3phi11=0;elseif U1=miuv3phi11 miuv4phi11=0;elseif U1=miuv1phi11 miuv5phi11=0;endW1=miuv1phi11 miuv2phi11 miuv3phi11 miuv4phi11 miuv5phi11W2=miuv1phi12 miuv2phi12 miuv3phi12 miuv4phi12 miuv5phi12W3=miuv1phi13 miuv2phi13 miuv3phi13 miuv4phi13 miuv5phi13W4=miuv1phi14 miuv2phi14