水资源短缺风险综合评价模型_数学建模竞赛参赛论文(编辑修改稿)

水资源短缺风险综合评价模型_数学建模竞赛参赛论文(编辑修改稿)内容摘要：

二、 北京没有天然湖泊，其大量用水靠水库提供，但北京的密云水库、官厅水库、怀柔水库、海子水库四座大型水库的水源入境量都有不同程度的衰减，近年来官厅水库尤为严重。 由于 各年 入境水量 不稳定 ， 所以 ， 应该大力改进修建水库 ， 扩大蓄水量 以保证水资源的合理利用。 三、 根据所给数据显示，北京市近年来的工、农业用水虽然都呈现下降趋势，但第三产业用水却呈现明显增长趋势。 为此，必须从各方面狠抓节约用水，以节约用水支持社会经济发展。 首先，要调整北京的工业结构布局，使高新技术产业成为首都经济发展的主导。 其次， 要 科学调整农业种植结构，建节水灌溉工程，并 加强城镇生活用水的管理。 最后，实行用水总量控制和定额管理相结合的制度，尽快 制定行业综合用水定额、居民生活用水定额，在实践中不断完善和改进 用水方案 ， 并研究科学可行的评估方法对各个用水单位进行用水评估 ，且对评估结果进行一定的奖惩措施。 四、 污水按性质分为 自然污染水和人为污染水。 北京当前对水体危害较大的是人为污染。 北运河是北京天然河道 五大水系中一条常年有水的河流， 但由于 城市污水污染、工业废水污染、农业回流水污染、固体废物污染等使得 北运河将在七年内 才能 完成综合治理。 对于污水处理问题，可以加大污水处理厂的建设， 也可以投入建设中水设施，以洗浴、盥洗等日常杂排水为水源，经过处理达到中水水质标准后，可以用于冲厕、洗车、绿化等。 因此，建立科学合理的水资源利用效率指标体系，研究科学可行的评估方法是很有必要的，对于节水 型社会建设具有重要的社会意义。 针对北京市的特殊背景和情况，上述建议为地方和行业水资源利用效率、调整产业结构布局提供了信息支撑。 以上报告，如无不妥，请批转各地执行。 通化师范学院 二 0一一年五月五日 模型 评价 与改进 11 本文提出的是使用模糊综合评价 的方法解决 水资源短缺问题，与其它计算方法相比，最主要的特点是： 模型的 优点 ： 本模型基于多个指标对北京市水资源短缺风险进行了综合评价，具有 一定的可操作性和实用性。 本模型成功地 分析了北京市水资源短缺 风险 问题 ，其中采用了计算机数据对实际理论情况的分析， 分析数据的 结果 表明此 模型 的建立 有效合理。 本模型与实际联系较为密切，综合实际情况进行 了 合理分析， 为北京市未来水资源短缺风险的预测提供了有效的理论依据，也为解决水资源短缺问题提供了切实可行的方案。 模型的 缺点 ： 由于现实生活中 水资源系统的 不确定性， 所以 ， 还有很多难以预料 的不 定 因素的影响 ，这使得运算结果 具有 一定的误差 ；同时，我们的理论体系还不够完善 ，需要进一步的改进。 模型的改进 ： 在计算概率 时 还可以应用 Logistic 回归方法进一步提高模型的精度， 由于时间关系，一些改进的思路还来不及实现。 模型推广 水资源是人类赖以生存和发展的重要资源之一，是不可或缺、不可代替的特殊资源。 通过定性分析影响水资源短缺风险因子，可以了解影响水资源短缺的因素，从而，从这些因素入手制定合理的方案降低水资源短缺现状。 模型推广之后，可以对北京市未来水资源短缺风险进行预测并提 出相应的应对措施 ，具有一定的实用性 和可操作性。 也可以借此方法来研究我国其他缺水地区水资源短缺风险。 参考文献 [1] 北京市统计局 .北京 2020 年统计年鉴 [M].中国统计 ,2020. [2] 刘卫国 .MATLAB 程序设计与应用 .北京：高等教育出版社， 2020. [3] 冯平 .供水系统干旱期的水资源风险管理 [J].自然资源学报， 1998,13（ 2）： 139—144. [4] 叶其孝 .大学生数学建模竞赛辅导教材 .长沙：湖南教育出版社， 1993. [5] 姜启源 .数学模型（第二版） .高等教育出版社， 1993. 12 附录 附表 1979 年至 2020 年北京市水资源短缺的状况 年份 总用水量(亿立方米 ) 农业用水(亿立方米 ) 工业用水(亿立方米 ) 第三产业及生活等其它用水 (亿立方米 ) 水资源总量（亿方） 1979 1980 26 1981 24 1982 1983 1984 1985 38 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2020 2020 2020 2020 2020 2020 2020 2020 2020 13 程序 1（计算 水 资源风险率 衡量 指标 ） clc clear M=[ ]。 %总用水量 G=[ 26 24 38 ]。 %水资源总量 C= M G。 D= M G M (C)=1。 M (D)=0 It=M %水资源系统状态量 Tsx=0。 for k= M Tsx=Tsx+k end clc Tsx %水资源状态总量 T=30 %水资源系统工作的总历时 U1=Tsx/T %水资源风险率衡量指标 运行结果： Tsx = 26 T = 30 U1 = %总用水量与水资源总量分布图 % M=[ ]。 14 plot(M,39。 ^39。 ) %总用水量分布图 % plot(G,39。 *39。 ) %水资源总量分布图 % 程序 2（计算 水 资源 脆弱性衡量 指标 ） clc clear M=[ ]。 %总用水量 G=[ 26 24 38 ]。 %水资源总量 C= M G。 D= M G M (C)=1。 M (D)=0 T1=0。 for k= M 15 T1= T1+k end clc T1 M=[ ]。 %总用水量 G=[ 26 24 38 ]。 %水资源总量 QZt1=MG B=QZt10。 QZt1(B)=0。 clc QZt1 M=[ ]。 %总用水量 G=[ 26 24 38 ]。 %水资源总量 QZt2=GM B=QZt20。 QZt2(B)=0。 clc QZt2 QZt=QZt1+QZt2 QZ=sum(QZt1+QZt2) clc QZt %第 t年水资源缺水量 QZ %水资源短缺总量 T1 %水资源系统失事总时间 X=QZ/T1 %失事状态下缺水量年均值 SZt=G。 %第 t年水资源总量 SZ=sum(SZt)。 %水资源总量 T=30 %水资源系统工作的总历时 Y=SZ/T %水资源总量年均值 16 U2=X/Y %水资源脆弱性衡量指标 运行结果： QZt = Columns 1 through 6 Columns 7 through 12 Columns 13 through 18 Columns 19 through 24 Columns 25 through 30 QZ = T1 = 26 X = T = 30 Y = U2 = %第 t年水资源缺水量与第 t年水资源总量 % plot(QZt,39。 g39。 ) %第 t年水资源缺水 plot(SZt) 17 %第 t年水资源总量 程序 3（计算 水 资源 重现期衡量 指标 ） clc clear M=[ ]。 %总用水量 G=[ 26 24 38 ]。 %。