水资源短缺风险评估模型建立与分析数学建模

1、目录摘 要 1一、问题提出 错误 ! 未定义书签。二、问题分析 错误 !未定义书签。三、模型的假设 错误 !未定义书签。四、模型的建立 错误 !未定义书签。4.1 模型一 错. 误 ! 未定义书签。4.1 1 符号说明 错误!未定义书签。4.1.2 模型一的建立 错误 !未定义书签。4.1.3 模型一的求解 错误 !未定义书签。4. 2模型二 错误!未定义书签。4.2.1 符号说明 错误 !未定义书签。4.2.2 模型二的建立 错误 !未定义书签。4.2.3 模型二的求解 错误 !未定义书签。五、模型的优缺点分析 错误 !未定义书签。六、城市水资源保护利用的建议报告 错误 !未定义书签。七、参

2、考文献 错误 !未定义书签。八、感谢信 错误 !未定义书签。九、附 录 错误 ! 未定义书签。附录一：涉及水资源短缺风险因子的各项数据表 错误 ! 未定义书签。附录二：求解逆矩阵的C+源代码错误!未定义书签。附录三：马氏距离 MATLAB军程序错误!未定义书签。水资源短缺风险评估模型建立与分析摘要目前，水资源问题， 特别是水资源短缺问题严重地制约着我国经济发展， 是 我们迫切所要解决的问题之一。 水资源是人类生产生活最重要的资源， 可是如今， 生态环境遭到严重破坏， 水体污染严重， 水资源开采过度， 水资源的保护和水污 染的治理成为现代社会最关注的问题。 尤其近年来， 我国北方地区水资源短缺问

3、 题日趋严重，已成为人们关注的焦点。本文以北京地区水资源短缺问题为背景， 通过建立线性加权均值风险评估模 型和短缺风险模糊综合评估模型， 对其水资源短缺情况进行综合风险评估。 同时 根据研究分析所得数据，给出了合理使用和保护水资源的建议报告。对于模型一：首先，我们借助软件 SPSS寸所搜集的数据进行筛选处理，按 照影响程度进行加工整合， 同时依据评指标选择原则， 选取了农业用水量、 生活 用水量、工业用水量、地下水资源量、污水排放总量、地表水资源量、入境水流 量、再生水资源量、环境用水量、COD放总量作为水资源短缺风险的评估因子。 其次，随机抽出20组数据分成两大组，作为矩阵 A和B,通过矩阵

4、运算、加权 均值、求解方程组，确定了各项因子的评分系数。在此过程中，我们通过编写 C+S序求解可逆矩阵和借助 MATLA计算方程组的解，从而确定评估指标函数， 在此基础上确定了水资源短缺评估的临界值， 并依次对水资源短缺风险预警级别 进行等级划分。 最后， 随机选择了 2009和2010年的数据进行检验， 得出评估模 型所求数据与实际情况基本一致。对于模型二： 首先，对影响北京水资源短缺的风险因子进行了分析， 得到了 一部分重要的风险因子。 其次，利用了判别分析法中的 Mahalanobis 距离法对主 要的敏感因子进行筛选， 取得了良好的结果， 基本符合客观事实。 针对水资源短 缺风险评估中

5、各指标的模糊性和不确定性， 建立短缺风险模糊综合评估模型。 通 过对风险率、脆弱性、可恢复性、事故周期、风险度这五个评价指标建立综合评 估体系，作出了等级划分。 通过对风险因子的调控， 降低了北京水资源短缺的风 险，并对未来北京市的水资源情况进行预测，对实际生活、生产具有指导意义。关键字： 水资源短缺 风险因子 线性加权均值风险评估模型 短缺风险模糊综合评估模型 风险预测 建议、问题提出近年来，我国北方地区水资源短缺问题日趋严重， 制约着经济的发展， 水资 源短缺已成为人们关注的焦点话题。以北京市为例， 北京是世界上水资源严重缺乏的大都市之一， 其人均水资源 占有量不足3OOm3,为全国人均的

6、1/8,世界人均的1/30，属重度缺水地区，附 表中所列的数据给出了 1979年至2000年北京市水资源短缺的状况。 北京市水资 源短缺已经成为影响和制约首都社会和经济发展的主要因素。 政府采取了一系列 措施, 如南水北调工程建设 , 建立污水处理厂 , 产业结构调整等。但是，气候变 化和经济社会不断发展， 水资源短缺风险始终存在。 如何对水资源风险的主要因 子进行识别， 对风险造成的危害等级进行划分， 对不同风险因子采取相应的有效 措施规避风险或减少其造成的危害， 这对社会经济的稳定、 可持续发展战略的实 施具有重要的意义。根据北京统计年鉴 及市政统计资料提供的北京市水资源有关信息， 分析

7、讨论以下问题：（ 1 ）以北京市水资源资料为例， 分析水资源短缺的风险因子， 并对这些风险 因子进行重要性分析；（ 2）建立一个数学模型对北京市水资源短缺风险进行综合评价，作出风险等级划分并陈述理由；（3）从用水量、 用水结构、水资源存量几个方面对北京市未来五年水资源进 行预测；（4）以北京市有关部门为报告对象， 写一份建议报告， 至少从水资源短缺成 因、水资源风险控制以及水资源保护几个方面提出建议和对策。二、问题分析2.1 水资源短缺风险评价指标水资源的短缺取决于供水和需水两方面影响， 而这两方面都具有随机性和不 确定性。因此，水资源短缺风险也具有随机性和不确定性。在进行风险评价时， 要充分

8、考虑风险的特点以及水资源系统的复杂性， 要把存在风险的概率、 风险出 现的时间、 风险造成的损失有多少、 风险解除的时间、 缺水量的分布等一系列因 素考虑在内。 因此难以用某一种指标对其进行全面描述和评价， 必须从多方面的 指标综合考虑。评价指标选择的原则是：（1）能集中反映缺水地区的缺水风险；能集中反映缺水风险的程度；（3）能反映水资源短缺风险发生后水资源系统的承受能力；（4）代表性好，针对性强，易于量化。依据上述原则，并参考相关文献，选取了影响水资源短缺的主要风险因子， 联系实际，建立了水资源短缺风险评估指标，构造出评估模型。2.2北京市水资源短缺主要原因（1）我国水资源时空分布不均以及入

9、境水流量少是导致北京水资源紧张的主 要原因；（2）北京人口增长过快，工农业发展迅速，加剧了水资源紧张程度；（3）北京地区生态环境恶化，干旱频率加大；（4）居民生活用水的污染和地表水资源急剧减少；三、模型的假设（1）假设题中所给数据基本真实有效；（2）假设没有重大的自然灾害发生如干旱、洪涝等；（3）假设北京地区人口流动正常；（4）假设南水北调及其它工程正常运行；（5）只考虑模型中的十个主要因子，其它的因子忽略不计。四、模型的建立4.1模型一4.1.1符号说明丫：样本数据的评分值C :对应影响因子的评分系数X :各影响因子实际数据丫0、丫1 :样本数据评分值的算数平均值X：、X1 : A组、B组原

10、始数据Xi0、X1 : A组、B组数据各列的平均值A= X0 Xi0B= X； X1S=（A+B）3Y、Y1 : a组、B组评分值的平均水平YC :橙色预警的临界值4.1.2模型一的建立水资源短缺风险评估模型是评价水资源短缺程度的评价指标，通过对北京市 水资源资料的调查分析建立这一模型。本模型主要通过对水资源短缺风险因子： 农业用水量、生活用水量、工业用水量、地下水资源量、污水排放总量、地表水资源量、入境水流量、它们分别是：人均用水量、环境用水量、COD排放总量进行了研究,1.农业用水量：Co为农业用水量风险评估系数,Xo为农业实际用水量2.生活用水量：Ci为生活用水量风险评估系数,Xi为生活

11、实际用水量3.工业用水量：C2为工业用水量风险评估系数,X2为工业实际用水量#X3为地下实际水资源量X4为污水实际排放总量X5为地下实际水资源量4. 地下水资源量：C3为地下水资源量评估系数,5. 污水排放总量：C4为污水排放总量评估系数,6. 地表水资源量：C5为地下水资源量评估系数,7. 入境水流量：C6为入境水流量评估系数，X6为实际入境水流量。8. 再生水资源量：C7为再生水资源量评估系数，X7为实际再生水资源量9. 环境用水量：C8为环境用水量评估系数，X8为环境实际用水量10. COD排放总量：C9为COD排放总量评估系数，X9为COD实际排放总 量。以上十项的线性组合Y=C0X0

12、+C1X1+C2X2+C3X3+C4X4+C5X5+C6X6+C7X7+C8X 8+C9X9记为水资源短缺风险评估指标值。（各项水资源短缺风险评估因子具体数据见附 录一）。下面根据收集的水资源短缺风险因子样本资料，求出各评估系数 C。从北 京市1979年至2001年水资源使用情况数据库中分别随机抽取 s和t个样本分别 记为A、B组。则它们对应的评估指标值为：组A的评估指标值丫0CX 00C1X 01C2X 02rCgX O8C9X09丫1cx 10C1X11C2X 12+CgX 18C9X 19丫2CX 20C1X21C2X 22* *CgX 28C9X 29YsC0X s0C1X S1C2X

13、 s2C8Xs8C9X s9组B的评估指标值丫0C0X 00C1X 01C2X 02“C8X O8C9X 09丫1C0X 10GX11C2X 12*C8X 18C9X 19丫2C0X 20C1X21C2X 22 *屮 * +C8X 28C9X 29YtC0Xt 0C1Xt1C2Xt2” C8X t8C9Xt9又做Y01 $Y，Y11 t Y，即丫0为组A的平均值，Y1为组B的平均S i 1t i 1值。为使组A与组B之间有明显区别，希望它们平均值之间差距越大越好，而组内离差平方和越小越好，即L Co ,Ci,C2, C3, C4,C5, C6, C7, Cg, C9(YY 一1)2(1)(YY

14、 _0)2i 1(Y1i 1Y1)25#越大越好，从而建立评估系数 Co,G,C2,C3,C4,C5,C6,C7,C8,C9为L Co1,23456C9,C ,C ,C ,C ,C ,C ,C 的极大值点。12由微分方程可知 C,G,C2,C3,C4,C5,C6,C7,C8,C9为方程组L C0 ,C1, C2, C3 ,C4 ,C5, C6, C7, C8, C9 Ci4.1.3模型一的求解1.将原始数据写成矩阵组A的矩阵0000 001 08009010X0 X011 180190900 0X91X9809924.421.722.728.819.918.918.117.418.516.51

15、3.911.812.011.813.411.713.914.914.815.113.812.311.915.515.314.613.811.811.110.813.714.613.115.114.615.613.219.315.918.38.39.18.58.49.09.29.310.812.513.05.86.36.57.36.65.77.111.35.48.11.21.40.91.71.51.11.91.91.12.02.62.12.12.01.91.66.07.05.57.70.60.70.40.51.51.32.44.25.64.10.060.030.030.020.060.070.0

16、50.040.010.01各列的平均(X0 X1 x8-X0) = -3.2 -4.3-3.6-3.7 -3.5 -3.1-2.3-1.5 -2.3-4.0(XX 一1 xX_ 19) = -3.8-3.1-2.7 -3.9 -3.6 -3.2-2.9-2.0 -3.5-2.324.231.831.628.831.612.010.813.012.812.414.413.812.213.911.215.714.714.816.518.68.39.18.58.49.07.85.36.18.27.61.81.31.11.21.61.61.92.12.02.60.30.80.60.61.10.020.

17、030.060.040.07组B的矩阵X00X01X08X；9X1。x；1 -X；8X；9x9。x9r-X；8X；921.810.119.59.721.912.713.913.713.814.214.417.210.014.014.018.516.221.415.117.89.29.310.812.513.06.07.612.86.87.21.01.61.81.82.23.65.06.06.56.71.62.73.23.63.90.080.050.060.050.082做新的矩阵A、B，及两组的算术平均值矩阵 SX00X0x0X01X10 X 0X08X；x0X09xAX?0X0x0X11X1

18、0X 0X18xfx 0X19X90x9。X0x0X91X10x0X98X：x0X99X020.927.52825.128.118.77.8187.4188.86.59.49.18.96.211.612.211.510.311.19.58.610.27.711.314.98.511.710.118.919.018.420.222.121.618.522.917.421.75.14.84.94.75.5-6.179.310.29.1119.69.711.911.19.19.914.39.111.15.05.64.74.95.14.13.93.34.16.14.86.211.75.76.16.77

19、.37.58.84.62.93.53.03.12.41.50.41.70.703.24.33.53.73.43.02.251.442.253.811111 111X00 X0 X01 X1 X08X8 X09X9BX；。X0 x11X11x：X8 X；9x9X；。X1 X；1 X；X：X8 x99X120.618.62024.916.415.715.215.41514.210.18.79.37.99.98.51112.911.312.810.09.29.211.611.811.410.99.87.68.517.517.715.81918.118.816.121.319.420.64.565.8

20、4.55.56.06.48.89.010.79.69.49.211.210.18.910.013.38.910.45.04.53.65.65.04.34.83.94.64.36.45.210.85.95.44.88.99.09.0103.22.42.33.42.01.90.52.22.11.83.73.12.73.93.43.12.91.93.52.3且S (A B)240.646.1485045.524.42333.422.432.218.815.218.71718.818.122.625.122.823.121.219.721.819.522.725.818.329.318.624.836

21、.436.734.239.240.240.424.644.236.842.319.610.810.79.211.012.113.418.119.219.8220.619.018.923.121.218.019.927.68.021.510.110.18.310.510.18.48.77.28.710.411.211.222.511.611.511.516.216.517.814.66.15.95.36.54.43.40.13.92.81.86.97.36.27.53.86.15.13.46.756.17满足式（1 的 Co,G,C2,C3,C4,C5,C6,C7,C8,C9 为方程组C0C1C

22、2C3C4X：C5的解,X31（2）8C8C9X?1X6X 1X；C6X001x；00.601X10X111.2X20x20.9X30x30.21X40x4；0.1S 10=S 1X5X50.1X60X60.6X70x70.5010.8X8X 8011.7X9X9C 0C1C2C 3C4C5C6C7C8C9解得:C0, C1,C2 , C3 ,C4, C5, C6 ,C7 ,C8 , C9= 6.4 5.2 3.8 4.5 2.3 1.4 1.6 0.87 2.12 0.03（其中S1通过C+设计程序求得，具体程序见附录二）3. 确定评估指标函数Y 6.4X05.2X；3.8X24.5X32.

23、3X41.4X51.6X60.87X72.12X80.03X9组A、B评估指标值的平均水平（2）#（2）#YA QX0 C1X100 0-0X8 C9X0= 95.541（2）#YB CoX； gX；一 C8X8 C9X；= 94.8799#再取出两组数据（每组20个），重复以上过程，最终求得临界值：Yc190.424. 水资源短缺风险评估等级梯度值划分为：黄色预警I级：得分147.44,虽然水资源处于短缺状况，但尚且可满 足工农业及生活用水，地下水资源开采较大，地表水、入境水用量较大，污水 处理率较高（再生水用效较多），水资源短缺风险较小。黄色预警II级：147.44得分168.93,虽然水

24、资源处于短缺状况， 但尚且可满足工农业及生活用水， 地下水资源开采较大，地表水、入境水用量 较大，污水处理率不高，水资源短缺风险也较小。黄色预警III级：168.93得分190.42,水资源处于短缺状况，不能 完全满足工农业及生活用水（有些产业、地区的用水被限制） ，地下水资源开 采较大，地表水、入境水用量较大，污水处理率不高，水资源短缺风险较大。橙色预警I级：190.42得分216.91,水资源处于短缺状况，不能完 全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开采 大（有些地区已无地下水可开采），地表水、入境水用量较大，污水处理率不 高，水资源短缺风险也较大。橙色预警II

25、级：216.91得分233.40,水资源处于短缺状况，不能完 全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开采 大（有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出现 河流干涸、湖水面积不断缩小），污水处理率不高，水资源短缺风险大。红色预警I级：233.40得分259.89,水资源处于短缺状况，不能完 全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开采 大（有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出现 河流干涸、湖水面积不断缩小），污水处理率低，水资源短缺风险大。红色预警II级：259.89得分286.38,水资源处于短缺

26、状况，不能 完全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开采大（有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出 现河流干涸、湖水面积不断缩小，甚至出现旱灾，严重影响了生产、生活），污水处理率低，水资源短缺风险很大。红色预警III级：得分V 286.38,水资源处于短缺状况，不能完全满足工 农业及生活用水（有些产业、地区的用水被强力限制），地下水资源开采大（有 些地区已无地下水可开采），地表水、入境水用量较大（有些地区出现河流干 涸、湖水面积不断缩小，甚至出现旱灾，严重影响了生产、生活、影响经济发 展），污水处理率低，水资源短缺风险很大。4.2模型二4.2.

27、1符号说明U评价对象的因素论域V评语论域W 各因素对水资源短缺风险指标的权重R模糊关系矩阵X表示横轴自变量g 123,4,5,6Vg表示纵轴因变量 g 1,2,3,4,5,6Cj表示编号为j的第i个因素的值i 1,2,3,4,5; j 1,2,3,30rj表示编号为j的第i个因素的值在30年的同一因素值的总和中所占的比 例A各风险因子的的权重分配rj表示因素论域U中第i个因素Ui对应于评语论域V中第j个等级Vj的 相对隶属度W表示各风险因子对水资源短缺风险指标的权重。4.2.2模型二的建立（1）水资源短缺的影响因子大概可分为两类：自然因素，即北京市总人口，水资源总量，水位深度，河 流和气候等；

28、社会经济因素，污水排放和处理，农业用水，工业用水，服务业用 水等由于经济发展对水资源需求。（2）判别分析判别分析可用于识别影响水资源短缺风险的敏感因子能够从诸多表明观测 对象特征的自变量中筛选出提供较多信息的变量且使这些变量之间的相关程度 较低线性判别函数的一般形式如下y ai 1 a2 2 , , a* n其中为y判别分数，1， 2,n为反映研究对象特征的变量，a1， a2, an为各变量的系数也称判别系数常用的判别分析法有欧式距离法和马氏距离法， 由于欧式距离判定法存在这 一定的缺点，故常用的判别分析方法是 Mahalanobis距离法，即每步都使得相距 最近的两类间的Mahalanobi

29、s距离最大的变量进入判别函数，其计算公式如下 ：d2 x,Yyii Xyi其中x是某一类中的观测量,Y是另一类,，可以求出x与丫的Mahalanobis距离。（3）水资源短缺影响因子筛选 根据上文中提出的水资源短缺风险影响因子， 利用Mahalanobis距离法筛选出水资源短缺风险敏感因子,见表5。从表5中第3 栏可以看出，水资源总量、工业排放总量、农业用水量、生活用水量在步骤1至步骤4中移出模型的概率均小于 0.1,同时在每步中这4个变量均使得最近的两 类间的Mahalanobis距离最大因此,这4个变量是影响北京地区水资源短缺风险 的敏感因子且影响力由大到小为：水资源总量 农业用水第三产业

30、用水 工业用水步 骤影响因素容许度移出概率最小马氏 距离的平方组间1工业用总量1.0000.0892工业用总量0.6820.0200.1842,5水资源总量0.6820.0000.2361,4工业用总量0.3910.0280.8461,53水资源总量0.6780.0000.7222,4农业用水量0.4600.0341.2252,5工业用总量0.2510.0356.5501,54水资源总量0.3280.0001.3852,4农业用水量0.1230.0031.2422,5第三产业用水量0.1020.0232.9632,5敏感因子筛选（4）模糊的分析：由于降水量，河流的流量的随机性，北京市供水和需水

31、量 也存在这不确定因素。因此水资源短缺也存在一定的随机性， 具有水资源短缺的 风险。所谓水资源短缺风险是指在特定的环境条件下，由于供水水和用水两方面 存在不确定性，使区域水资源系统发生供水短缺的概率以及由此产生的损失。水资源短缺风险评价指标主要描述水资源短缺的状况，只针对水资源的系统不包括 社会经济系统，是水资源短缺风险却化的基础。选取风险率、脆弱性、可恢复性、 重现期和风险度几个评价指标利用模糊综合评价方法对水资源短缺风险进行综 合评价。（5）评价指标1. 风险率根据风险理论，载荷是造成系统非正常状态的动力，抗力是维护系统正常的 能力。如果把水资源系统的失事状态记为 F R，正常状态记为S

32、R那么水资源系统的风险率为r p Rp x F其中，人为水资源系统状态变量。如果水资源系统的工作状态有长期记录风险率也可以定义为系统不能正常工作的时间和整个系统工作时间的比值，即：1 NSa hNS t i其中，NS为水资源系统工作的总时间；lt为水资源的状态变量0,系统正常工作 j S1,系统失事 t F2. 脆弱性脆弱性是描述水资源系统失事损失严重程度的期望值的大小。假定系统第1 次失事的损失程度为，它相对应的发生概率为，那么水资源系统的脆弱性为：NFE SRSjt其中，NF为系统失事状态的总次数3. 重现性周期是两次进入非常状态模式 F之间的时间间隔，也叫平均重现期。用d ,n表示第n间

33、隔时间的历时，则平均重现期为:,n其中，NN 为0,t时段内属于模式F的事故数目4. 可恢复性可恢复性是指系统从失事状态恢复到正常状态的可能性。恢复性越高说明系统从失事状态转入正常状态的时间越短。可以用下面条件概率来定义:P Xt 1 F,Xt S引入整数变量:P Xt i F1,Xt F0,Xt S1, XtF,Xt S0, 其他这样可得到:NFNF乙/t 1t 15. 风险度用概率分布的数学特征，如标准差C，可以说明风险的大小。（7 -越大，则风险越大，反之越小。这是因为概率分布越分散,实际结果远离期望值的概率就越大。14#1/2D XXii 11/22E X / n 1423模型二的求解

34、采用上述定义的的风险率、脆弱性、可恢复性、重复性、风险度作为评价指 标，采用模糊综合评价方法对风险进行评价。设定两个有限论域（即集合）：Uq,U2，u3,um和Vv,v2, v3,vm其中，U表示综合评判的因素所组成的集合，V表示评价语组成的集合。则模糊综合评判即表示下列的模糊变换 B AR ，式中A为U上的模糊子集。而评判结果B是V上的模糊子集，并且可 以表示A 1, 2,，m , 0 i 1; B Db，,bm ,0 bj 1;。其中i表示单因 素Ui在总评定因素中所起作用大小的变量，也在一定程度上代表根据单因素Ui评定等级的能力；bj为等级Vj对综合评定所得模糊子集B的隶属度，它表示综

35、合评价结果。关系矩阵R可表示为：r11r12 r1nr21广 【22OnRr31 r323nrm1rm2式中：rj表示因素的评价对等级Vi的隶属度，因而矩阵F中 R mg,扁即 为对第i个因素的单因素评判结果。在评价计算中 A 1, 2,，m代表了各个 因素对综合评判重要性的权系数，因此满足i 1, i 1,2，,m ;同时，模糊变换A R也即退化为普通矩阵计算，即bj min 1, 飞，i 1,2,m;j 1,2，n。上述权系数的确定可用层次分析法 AHP得到。由上述分析可以得到：评价因素集U对应评语集V,而评判矩阵中用即为某因素Ui对应等级Vj的隶属度，其值课根据各评价因素的实际数值对照各

36、因素的 分级指标推求。我们把评语级划为 5个如下等级：水资源短 缺风险风险脆弱性u2可恢复U3重现性U4风险度U50.2000.2000.2000.2000.20020.200 0.4000.200 0.4000.200 0.4000.200 0.4000.200 0.40030.400 0.6000.400 0.6DGD.400 0.6000.400 0.6000.400 0.60040.600 0.8000.600 0.800).600 0.8000.600 0.8000.600 0.80050.8000.8000.8000.8000.800对于水资源系统的风险率，脆弱性，风险度是越小越好

37、，所以按越小越 优性指标对u1,u2,u5构建如下的隶属函数：1,ai1ai1ai1ai2ai21, ai1ai2,ai1ai2ai1V2ai3420,ai2ai3ai20,ai315，ai10,ai2ai2,ai1a?Vj1,ai2ai3ai3aai4ai3,ai3 ai2,ai4ai2ai3V5a3a；3ai4ai31,ai3ai3ai4ai4V4ai31,ai3ai3 ai2ai3ai4ai4对于水资源系统可恢复性和重现性则是越大越优性指标，对U3 , U4各评价语级构建如下隶属函数:1,ai13i2ai2ai2,ai1ai2ai1,ai13i23i2ai1V31,ai3ai2ai3ai

38、3ai3ai4,ai4,ai4,9i3V5a47ai4ai31,ai4对于水资源短缺风险评价的因素集ai17ai11,ai1ai2V2ai3ai29i3 :3i2 ,a 3,ai2ai3,9i33i2V41,ai4a3ai47ai4U而言，对应一个测定指标向量丫 11, 12, 13, 14, 15 其中j是U寸于Uj的测定值。这样Uvij就是相对于Ui属于Vi的程度。对于因素集U便有一下的模糊关系矩阵17v1 11v2 11v3 11v1 12RUv1 13v2 12v2 13v3 12v3 13v412v512v413v513v4 11v5 11v4 14v5 14v1 14v2 14v3

39、 14v1 15v2 15v3 15v4 15v5 15水资源短缺风险评价各因素的权重确定采用层次分析法(AHP),设权重计算18结果为A 1, 2, 3，4, 5，于是可得出综合评判向量B A?R我们选取“加权平均型”的模型进行综合评价，即bj min 1, 百n山1，所以这个模型可以用一般实数加法，即：i 1nbj舟 j 1,2,3，mi 1max bij对应的评价语即为水资源短缺评价的评价结果。对风险程度的综合评价可以将其划为5级，制作如下的级别评价表格分别 为：低风险，较低风险，中风险，较高风险，高风险。水资源系统风险等级划分风险等级风险级别风险特征1低风险可以忽略的风险2较低风险可以

40、接受的风险3中风险边缘风险4较高风险比较严重的风险5咼风险无法承受的风险五、模型的优缺点分析模型一的优缺点：本模型的优点是操作简单，实用性强。通过大量随机数据的验证，模型基本 符合实际情况，对水资源短缺风险评估具有一定的帮助作用。 另外，在模型的计 算中进行了 C+程序和MATLABk机实现，对以后水资源短缺风险评估标准的程 序化、软件化具有一定的指导作用。本模型的缺点是考虑的情况较为理想化，有地区局限性。模型二的优缺点：一般风险研究只对个别风险指标进行描述， 而对水资源短缺的综合评价则比 较少，而模糊综合评价模型从不同角度对水资源系统进行了比较全面的风析， 与 实际情况比较符合。 在分析敏感

41、因子简化问题， 抓住了主要原因， 取得了良好的 结果。本模型的缺点是模型的使用范围存在着一定的局限性。六、城市水资源保护利用的建议报告 目前，北京市水资源短缺问题日益严重，已沦为世界上水资源严重匮乏的大 都市之一。 水资源短缺问题影响着人们生产生活的正常进行， 制约着北京市经济 的可持续发展。 通过对水资源状况的综合分析， 我们建立了线性加权均值风险评 估模型和短缺风险模糊综合评价模型， 得出了影响北京市水资源短缺的主要风险 因子，划分了水资源短缺风险预警级别， 并且对北京市未来水资源短缺情况进行 预测。对此，我们提出以下几点浅陋建议，希望对水资源的保护利用有所帮助。1. 加强水土保持与综合治

42、理 日益严重的水土流失使地面土壤表层不断受到破坏， 造成土壤养分流失， 制 约了经济发展，导致农民生活贫困， 更为重要的是破坏了流域的生态平衡。 因此， 治理水土流失，开发当地水土资源是改善当地生态环境， 促进国民经济可持续 发展的主要措施， 可通过以下措施加强水土保持及综合治理工作。(1) 通过生物措施工程措施和耕作措施的合理配置，使流域形成由水、林、 田、路统一规划，坡面、沟道均有水保工程的立体防护体系。(2) 根据土地资源适宜性评价结果，确定流域经济与生产发展方 向将陡坡地段的坡耕地退耕还林，增加林地用地比例，大举植树造林 或封禁治理，提高土地利用率，减少水土流失面积， 逐渐使流域的土地

43、利用结构 趋于合理，农林牧各业协调发展，改善生态环境。2. 加强水资源环境保护与污染治理切实加强水资源的环境保护， 保护现有植被， 减少水土流失及地面径流对工 业园区设置单独的污水处理厂， 严禁污水直接排入自然水体将合理利用水资源与 控制水污染相结合， 提高水资源重复利用率， 近、中、远期分别达到 60%、70% 、 80%以上建筑工地废水必须经过处理回用或合格排放，禁止无组织占用道路经营 洗车。通过固体废物控制，优化污染治理，进而推进水资源保护工作。(1) 生活垃圾的处理采用定点收集垃圾，集中无害处理，焚烧和填埋并行， 并发展综合利用家庭生活垃圾收集 : 实施家庭垃圾分类袋装化， 根据垃圾的

44、可否 再生利用处理难易程度等特点， 由家庭事先进行分类袋装， 在居住区设置专用垃 圾收集间和特定集装箱， 在商业中心设置垃圾收集站， 对商业垃圾进行分类收集， 并纳放家庭垃圾收集系统。(2) 工业垃圾的处理以填埋方式为主。(3) 医疗废弃物必须单独收集运输并用焚烧法处理，避免与一般垃圾相混。(4) 建筑垃圾以填埋方式为主，对施工产生的余泥尽可能就地回填或申报有 关管理部门及时运走， 堆放到合适地方， 并注意清洁运输，防止运输过程中洒落， 对于建筑固体废物除部分木材和竹料经再加工可再利用外，一般不能重新利用， 需要堆置存放。3. 实行严格的水资源管理制度结合现有水资源管理的经验，首先在水资源管理

45、观念上应实现六个转变 : 一 是在管理理念上， 要加快从供水管理向需水管理转变 ; 二是在规划思路上， 要把 水资源开发利用优先转变为节约保护优先 ; 三是在保护举措上， 要加快从事后治 理向事前预防转变 ; 四是在开发方式上， 要加快从过度开发无序开发向合理开发 有序开发转变 ; 五是在用水模式上， 要加快从粗放利用向高效利用转变 ; 六是在 管理手段上， 要加快从注重行政管理向综合管理转变。其次，实行严格的水资 源管理制度， 严格水资源开发利用红线、 水功能区限制纳污红线和用水效率控制 红线，通过完善的水资源法规体系，强有力的水资源监管措施， 坚决打击违法 填湖行为，强化取水资源论证、节水

46、考核、入河排污口设置、水域岸线利用、河 道采砂、水工程建设等方面的管理和执法监督， 严格执行水资源管理的各项法规。4. 加强用水节能教育，推进全面节水工程(1) 推动节水型社会建设。认真研究提出关于开展节水型社会建设的指导性 文件，适时召开全国节水型社会建设工作会议， 继续开展全国节水型社会建设试 点工作，研究提出水资源宏观分配指标和微观取水定额指标， 推进国家水权制度 建设。(2) 开展城市节水工作。积极开展节水产品研发，加大节水设备和器具的推 广，指导各地加快供水管网改造，降低管网漏失率，推动公共建筑、生活小区、 住宅节水和中水回用设施建设， 推进污水处理及再生利用， 加快城市供水和污水

47、处理市场的改革，大力开展雨水回收利用， 建筑物安装屋顶虹吸式雨水收集系统， 通过这个收集系统， 尽量收集更多的雨水， 采用植草砖、 渗透砖收集地表径流水， 经处理后做景观绿化用水，并在适当的地方建设地窖存放丰雨期收集到的雨水。(3) 推进农业节水。继续推进农业节水灌溉，推广农业节水灌溉设备应用， 大力推进大中型灌区节水改造，积极开展农业末级渠系节水改造试点。在丘陵、 山区和干旱地区积极开展雨水积蓄利用， 支持农村水窖建设， 推广旱作农业技术， 发展旱作节水农业， 扩大节水作物品种和种植面积， 开展农村、 集镇生态卫生旱 厕试点。(4) 推进节水技术改造和海水利用。推进高耗水行业节水技术改造，矿

48、井水 资源化利用，推进沿海缺水城市海水淡化和海水直接利用。总之，合理利用并保护水资源是北京市自然环境与自然资源保护中的一项重 大任务，只有妥善解决，才能实现北京市社会稳定，经济的平稳发展。以上所述几点建议，只是我们个人浅显的看法和认识。由于所学知识和阅历 有限，言辞中难免有所不当，望各位领导批评指正！七、参考文献1 阮本清 , 梁瑞驹 , 陈韶君 . 一种供水系统的风险分析与评价方法 J. 水利学 报,2000,(9):17.2 傅湘 , 王丽萍 , 纪昌明 . 极值统计学在洪灾风险评价中的应用 J. 水利学 报,2001,(7):610.3 阮本清,韩宇平,王浩等, 水资源短缺风 险的模糊综 合评价 J. 水利学 报,2005(8).4 Iman Karimi,Eyke Hllermeier,Risk assessment system ofnatural hazards : A new approach based on fuzzy probability, Fuzzy Sets and Systems.2007,987-999.5 杨树滩, 夏自强.模糊数学在水资源紧缺程度评价中的应用研究 J. 长江科学 院院报 ,2005,22(