水资源短缺风险综合评价优秀论文

水资源短缺风险综合评价摘要:水是人类赖以生存的资源，近些年来随着社会经济的迅猛发展和全球气候的快速变化，水资源短缺风险问题日益突出。政府部门希望可以得到关于水资源短缺风险的综合评价，并根据评价结果提出有效地应对措施以降低风险。本问提出了利用层次分析法确定主要风险因子，然后根据缺水风险的可能性和水资源短缺的程度对水资源短缺风险进行综合评价，并作出风险等级划分，再对未来两年的风险进行预测，最后给北京水行政主管部门写一份建议报告。问题一：我们先从供水，用水等角度考虑确定出六种风险因子（农业用水，工业用水，第三产业及生活等其他用水，水资源总量，降水量，常住人口），然后从客观分析和主观分析两方面考虑，使用层次分析法确定出六种风险因子的权重，最后我们选取权重大的的风险因子作为主要风险因子。问题二：我们认为水资源短缺风险包含水资源短缺程度所决定的。水资源短缺程度越大，那么风险就越高。再根据评价结果划分出风险的五个等级为：低风险、较低风险、中风险、较高风险、高风险，最后给出对于主要风险因子应该如何进行调控使得风险达到最小。在问题三：首先需要根据历史数据拟合和预测出影响风险的六个风险因子各自的值，然后运用问题二中德模型对北京未来两年水资源的短缺风险进行预测。问题四：以北京市水行政主管部门为对象写一份建议报告，在建议中分析北京市水资源短缺存在的风险，并给出相应的应对措施，以供政府部门作出科学的决策。关键词：层次分析法，Matlab，函数拟合，风险等级划分。问题重述:水资源，是指可供人类直接利用，能够不断更新的天然水体。主要包括陆地上的地表水和地下水。风险，是指某一特定危险情况发生的可能性和后果的组合。水资源短缺风险，泛指在特定的时空环境条件下，由于来水和用水两方面存在不确定性，使区域水资源系统发生供水短缺的可能性以及由此产生的损失。近年来，我国、特别是北方地区水资源短缺问题日趋严重，水资源成为焦点话题。以北京市为例，北京是世界上水资源严重缺乏的大都市之一，其人均水资源占有量不足300m3《北京2009统计年鉴》及市政统计资料提供了北京市水资源的有关信息。利用这些资料和你自己可获得的其他资料，讨论以下问题：评价判定北京市水资源短缺风险的主要风险因子是什么？影响水资源的因素很多,例如：气候条件、水利工程设施、工业污染、农业用水、管理制度，人口规模等。2建立一个数学模型对北京市水资源短缺风险进行综合评价，作出风险等级划分并陈述理由。对主要风险因子,如何进行调控，使得风险降低？3对北京市未来两年水资源的短缺风险进行预测，并提出应对措施。4以北京市水行政主管部门为报告对象，写一份建议报告。基本假设:题中和我们所查到的数据都是真实可靠的。在未来两年内不会发生地震，干旱，水灾等严重自然灾害。每一年的相同风险因子之间存在一定的关系（可能是函数，也可能是走势）。符号说明：Ti:i为1,2,3,4,5,6时分别指的是农业用水，工业用水，第三产业及生活等其他用水，水资源总量，降水量，常住人口得加权系数。Yj：代表每一年的水资源总量，j可以去1,2……30，分别代表从1979年到2008年Zj：代表每一年的用水资源总量，j可以去1,2……30，分别代表从1979年到2008年Mj：代表每一年的缺水程度，j可以去1,2……30，分别代表从1979年到2008年。Mj1：代表修正后缺水程度，j可以去1,2……30，分别代表从1979年到2008年。Yt：t取1,2,3分别代表农业，工业，第三产业的用水量。Xt：t取1,2,3分别代表农业，工业，第三产业的相应年份。5.问题分析：问题一：我们认为影响水资源的风险因子有农业用水，工业用水，第三产业及生活等其他用水，水资源总量，降水量，常住人口。我们运用层次分析法得出个因数的比例系数。最后以比例系数大的作为主要的风险因子。问题二：对于风险的评估我们分为两个方面。一是发生缺水可能性的大小，发生的可能性越大这说明风险越大。二是指缺水的程度，缺水程度越大着风险越大。综合这两方面的大小，再结合问题一的比例系数。我们可以得出缺水风险的大小。再进行相应的等级划分。在问题三：首先需要根据历史数据拟合和预测出影响风险的六个风险因子各自的值，然后运用问题二中德模型对北京未来两年水资源的短缺风险进行预测。问题四：以北京市水行政主管部门为对象写一份建议报告，在建议中分析北京市水资源短缺存在的风险，并给出相应的应对措施，以供政府部门作出科学的决策。6模型的建立和求解：6.1问题一模型的建立与解答：6.1.1构造判断矩阵：根据每两个指标的相互重要性，确定出每两者之间重要性的比值。得到的矩阵如下：6.1.2判断矩阵的计算：我们用matlab软件求解出最大特征根max=6.1857，所对应的特征向量A=(0.3965,0.1560,0.0748,0.4006,0.1925)。6.1.3一致性检验：计算出CI==0.1857/5=0.037。RI=1.24。(CI/RI)<0.1。通过一致性检验。将加权系数归一化得：A=(0.1978,0.0778,0.0373,0.1998,0.00960,0.3913)。所以T1=0.1978T2=0.0778T3=0.0373T4=0.1998T5=0.0096T6=0.3913所以得出农业用水，水资源总量和人口数量是主要的风险因子。6.2问题二模型的建立和解答：6.2.1缺水程度的计算。Mj=（Zj-Yj）/Zj。可以得到每年的缺水程度如下：年份缺水程度年份缺水程度年份缺水程度年份缺水程度19790.1092731987-0.2491119950.32397520030.48603419800.48555619880.0765971996-0.1464620040.38150319810.50114319890.51724919970.44816520050.32753619820.22490519900.12791819980.06752420060.28571419830.2703951991-0.0061919990.65907520070.31609219840.01847719920.51669220000.58267320080.0256411985-0.1983619930.56501520010.50642719860.26046519940.0098120020.5346826.2.2缺水程度的修正：我们将缺水程度全部转化为正值且原缺水程度越大的的得到的修正后的缺水程度值也越大。转化公式如下：Mj1=（Mj-Mjmin）/（Mjmax-Mjmin）。转化后的缺水程度数据如下：年份修正后的缺水程度年份修正后的缺水程度年份修正后的缺水程度年份修正后的缺水程度19790.394619880.358619970.767820060.588919800.808919890.843819980.348620070.622319810.826119900.415119991.000020080.302519820.521919910.267520000.915919830.572019920.843220010.831919840.294619930.896420020.863019850.055919940.285120030.809519860.561119950.631020040.69441987019960.113020050.63496.2.3风险程度的划分:我们将风险划分为低风险、较低风险、中风险、较高风险、高风险五个等级。每种风险的范围低风险[0-0.2)较低风险[0.2-0.4)中风险[0.4-0.6)较高风险[0.6-0.8)高风险[0.8-1.0)6.2.4得出上面各年份的风险程度:年份风险等级年份风险等级年份风险等级年份风险等级1979较低风险1988较低风险1997较高风险2006中风险1980高风险1989高风险1998较低风险2007较高风险1981高风险1990中风险1999高风险2008较低风险1982中风险1991较低风险2000高风险1983中风险1992高风险2001高风险1984较低风险1993高风险2002高风险1985低风险1994较低风险2003高风险1986中风险1995较高风险2004较高风险1987低风险1996低风险2005较高风险6.2.5对于主要风险因子的调控：考虑到影响北京缺水的主要风险因子有农业用水，水资源总量，人口数目。农业用水的建议：我们可以适当的改变一下农业用水的模式，例如：改进田间灌水方法与技术：对量大面广的地面灌溉区，平整土地，缩小地块，改进沟畦尺寸规格，控制入沟畦流量；在低压管道输水灌区，采用“小白龙”薄膜塑料管等简易田间管道灌溉方法；在水田地区，推行水稻“浅、湿、晒”控制灌溉技术；在大田旱作物地区，推行浇关键水非充分灌溉技术；在不具备发展常规灌溉的地方，推行半机械化、畜力、人力等多种形式的注水抗旱播种保苗技术；在经济发达地区、城市郊区、经济作物种植区，逐步推广喷灌、滴灌、微喷灌等技术。

耕作栽培与生物节水技术：采用深翻松土，耙耱保墒，增施有机肥，改善土壤结构等方法提高土壤蓄水能力；推广少耕、免耕方法，秸杆覆盖、地膜覆盖、间作套种等方法，增强土壤保墒能力；施用“旱地龙”等生物抗旱剂和土壤保水剂，减少土壤和农作物水分无效蒸发。以达到节约用水的目的。采用南水北调以达到加大农业灌溉面积，满足农业灌溉用水的目的。2.水资源总量：继续开展污水资源化、雨洪利用的研究和应用。把城市污水排放规划管理、污水处理厂建设、再生污水利用三个环节综合起来，全面规划考虑，实现污水资源化。收集和利用城市雨洪，既可防治雨洪灾害，缓解城市雨洪压力；同时又增加了可用水资源，并可通过回灌补给蓄养地下水。以此来加强水资源的储存总量。3.人口总量：对于全国在以后增量的各种资源中，不能再大量集中在中心城市，而要向中小城市倾斜，重点发展中小城市，缩小中小城市在教育、医疗等方面与中心城市的差距，通过建立高屋建瓴的公共政策，使得人才、教育、医疗等资源不再集中于大城市。“围魏救赵”，分散资源，提升其他城市的品质，这样就能减少北京等大城市的人口，一次来减小人口过大带来的压力。6.3问题三：6.3.1我们采用拟合的方法得出2009和2010年的总用水量和水资源总量。对于总用水量我们采用直接拟合和间接拟合两种方法。然后进行对比和数据的修正。得出符合实际的数据。6.3.2总用水量的拟合：直接拟合：图形如下：间接拟合：图形如下：农业用水：

Y1=0.0006*x1.^2+（-24.4521*x1）+2.4888；工业用水：Y2=-0.0199*x2.^2+（78.9436*x2）+-7.839；第三产业用水：Y3=0.004*x3.^2+（-15.5358*x3）+1.5047；总用水量=农业用水+工业用水+第三产业用水。6.3.3水资源总量的拟合：图形如下：运用第二问的模型我们可以得出2009和2010年的风险程度和风险等级为：修正后的风险程度风险等级2009年0.711较高风险2010年0.378较低风险6.4问题四：关于如何降低北京市水资源风险的报告北京市水行政主管部门：你们好！我们通过对北京市农业用水，工业用水，第三产业及生活用水，水资源总量，降水量和人口规模的分析研究，确定出人口规模，水资源总量和农业用水为造成水资源短缺风险的主要风险因子。然后我们从风险造成的后果进行水资源短缺综合评价，最后将风险等级划分为五级，可以看出我们应该尽量控制人口和节约农业用水，加强水资源总量的储备，以达到降低北京市缺水风险的程度。一、加强水资源的合理调度，完善了相关体制和机构建设：1．加强地下水环境保护，建立完善的地下水动态监测系统。一是严格执行地下水保护条例，坚决控制地下水的开采。今后，在地下水超采区、海水入侵区、城市管网区和地下水水质受到污染地区一律停止审批地下水，在其他地区也要严格限制开采地下水。对非法盗采地下水，要坚决予以打击，特别要打击在工业园区、开发区和管网区内偷打地下井行为，更不能违规审批新的地下水取水工程。省水利厅将加大监察力度，发现一起，查处一起，并向全省通报，坚决杜绝私打乱采地下水现象，保护好地下水资源。二是采取强制性措施，大力推进地下水超采区、海侵区和管网区的综合治理。对目前已经开采的地下水，逐步削减开采量；对能够利用水利工程和城市自来水替代的，有计划地关闭地下水井。2．污水资源的再利用。第一，工业用水的回用，经处理的污水可用于非饮用水的其他用水途径。第二，居民及社区的非饮用水回用，对居民来说，污水回用可以用于冲洗马桶，浇花园，洗车，清洗等用途。第三，间接饮用水循环，经过适当的处理可成为可饮用水。第四，直接用水的循环利用。3．应建立外来水源、如南水北调工程。本地水源相互协调的供水网络，实现本地地表水源与外来水源的联合调蓄、地下水与地表水的联合调蓄，提高北京城市供水安全保证程度，支持城市可持续发展。4．为保证城市供水安全，应科学地适度增加地下水开采量，合理开发利用。对已确定的应急供水水源地应尽快投入勘探和开发工作，对其它地区继续开展调查工作，寻找新的后备应急水源。二、使有关单位及市民注意节约用水，提高用水效率：1．通过宣传加强民众和单位的节水意识，另一方面有关部门需要采取一些有力的甚至是强制性的节水措施。例如，继续努力提高用水效率，缩小与世界平均用水水平的差距；鼓励清洁生产、节水生产，实行定额用水制度；还有，要根据各地的不同情况，适时、适度地提高水价，逐步改变水价格背离其价值的情况。2．节约用水是当务之急也是长远发展战略方针，在优先保证城市生活和重点工业供水的前提下，在无法满足需水时，适度压缩农业用水。加强工业、农业节水力度，调整产业结构，大力发展节水型工业、农业。把节约用水纳入城市发展规划，纳入产品结构调整计划和技术及企业改造计划，使在城市和工业部门中逐步做到计划用水、合理用水和科学用水。3.可以通过行政法的手段禁止容易出现跑、冒、滴、漏的用水装置的生产和使用，现今我国已有相应的标准。再如对于漏水的控制，毫无疑问，漏水企业将承担更多水费支出，但这并不足够，还必须通过行政执法部门，对用水单位的用水管道和装置进行定期检查或抽查，发现跑、冒、滴、漏水现象及时予以纠正和处罚。7.模型的检验：我们查阅相关数据得出2009和2010年的相应数据，与我们得出的数据进行比较如下：7.1拟合的数据：总用水量水资源总量2009年35217.2查阅到的数据：总用水量水资源总量2009年35.521.80.38597.3数据误差分析：缺水程度误差=|拟合数据得到的缺水程度-查阅数据得到的缺水程度|/查阅数据得到的缺水程度。缺水程度误差分析2009年3.6%所以本模型具有很好的准确性，能够较好的预测未来几年的缺水情况。8.模型有缺点的分析8.1模型的优点：本模型能够比较直观的反应出隔年的缺水情况（即缺水的风险等级）。本模型对未来几年的缺水情况能够较为准确的进行预测。8.2模型的缺点：本模型没有考虑到每一年发生缺水情况的概率，此处需要加以改进。9.参考文献：[1]姜启源，数学模型，北京：高等教育出版社，1998[2]西北工业大学数学建模指导委员会，数学建模简明教程，200810.附录：10.1程序：10.1.1问题一：10.1.2问题三：农业用水：a=polyfit(x,y,2);

Warning:Polynomialisbadlyconditioned.AddpointswithdistinctX

values,reducethedegreeofthepolynomial,ortrycentering

andscalingasdescribedinHELPPOLYFIT.

>Inpolyfitat80

>>x1=linspace(1979,2008,30);

>>y1=a(1)*x1.^2+a(2)*x1+a(3);

>>plot(x,y,'o')

>>holdon

>>plot(x1,y1,'r')

>>legend('数据点','拟合函数')工业用水：a=polyfit(x,y,2);

Warning:Polynomialisbadlyconditioned.AddpointswithdistinctX

values,reducethedegreeofthepolynomial,ortrycentering

andscalingasdescribedinHELPPOLYFIT.

>Inpolyfitat80

>>x1=linspace(1979,2008,30);

>>y1=a(1)*x1.^2+a(2)*x1+a(3);

>>plot(x,y,'o')

>>holdon

>>plot(x1,y1,'r')

>>legend('数据点','拟合函数')第三产业：a=polyfit(x,y,2);

Warning:Polynomialisbadlyconditioned.AddpointswithdistinctX

values,reducethedegreeofthepolynomial,ortrycentering

andscalingasdescribedinHELPPOLYFIT.

>Inpolyfitat80

>>x1=linspace(1979,2008,30);

>>y1=a(1)*x1.^2+a(2)*x1+a(3);

>>plot(x,y,'o')

>>holdon

>>plot(x1,y1,'r')

总用水量：xx=1979:1:2008;

xx=1979:0.05:2008;

yy=spline(x,y,xx);

plot(x,y,'o',xx,yy);

gridon

水资源总量：xx=1979:1:2008;

xx=1979:0.05:2008;

yy=spline(x,y,xx);

plot(x,y,'o',xx,yy);

gridon

10.2数据：农业用水工业用水第三产业及生活等其他用水总用水量（亿立方）水资源总量（亿立方）水资源短缺（亿立方）1979年24.1814.374.3742.9238.23-4.691980年31.8313.774.9450.5426-24.541981年31.612.214.348.1124-24.111982年28.8113.894.5247.2236.6-10.621983年31.611.244.7247.5634.7-12.861984年21.8414.3764.01740.0539.31-0.741985年10.1217.24.3931.71386.291986年19.469.917.1836.5527.03-9.521987年9.6814.017.2630.9538.667.711988年21.9914.046.442.4339.18-3.251989年24.4213.776.4544.6421.55-23.091990年21.7412.347.0441.1235.86-5.261991年22.711.97.4342.0342.290.261992年19.9415.5110.9846.4322.44-23.991993年20.3515.289.5945.2219.67-25.551994年20.9314.5710.3745.8745.42-0.451995年19.3313.7811.7744.8830.34-14.541996年18.9511.769.340.0145.875.861997年18.1222.25-18.071998年17.3910.8412.240.4337.7-2.731999年18.4510.5612.741.7114.22-27.492000年16.4910.5213.3940.416.86-23.542001年38.919.2-19.72002年15.57.511.634.616.1-18.52003年13.88.413.635.818.4-17.42004年13.57.713.434.621.4-13.22005年13.26.814.534.523.2-11.32006年34.324.5-9.82007年12.45.816.634.823.8-112008年125.217.935.134.2-0.92009年12.05.218.335.521.8-13.7年份常住人口(万人)按性别分按城乡分城镇人口乡村人口人口出生率(‰)人口死亡率(‰)人口自然增长率(‰)男女197887244342947939312.936.126.81197989745544251038713.675.927.75198090445844652138315.566.39.26198598150048158639515.455.759.71986102852450462140715.824.4711.351987104752552263741017.295.411.891988106153452765041114.435.089.351989107553853766443112.845.537.491990108654554167341313.045.757.23199110945475476834118.035.822.211992110255454869241199311125595537074059.356.163.19199411255645617254008.965.763.2199512516276248154367.925.122.8199612596396208294308.025.342.68199712406296118264147.916.021.891998124663161583840865.30.7199912576366218544020001364711653105830200113837216621080302002142374368011183056.65.730.872003145676169511513055.065.15-0.09年份户数总人口非农业人口1949－－－1950－204.3－1951－222.0－1952－251.9－1953－287.1－195462.5310.4－195566.1321.0－195681.0383.8301.5195786.4401.1－1958116.6552.3－195985.2426.6－196085.7470.3－196192.1456.1396.7196292.8448.6387.7196394.1462.0398.5196495.2472.2407.0196595.5476.3410.41966－465.0397.01967－471．0403．01968－465.0394．01969－443.0371.01970－439.1366．01971101.9445.1371.21972104.3452.1378.01973106.4457.5381.41974109．0464.8389.41975111．0465.7389.91976113.6468.9392.81977117.7473.9396.81978122.1486.0409.31979128.4509.5434.11980135.9530.9454.81981143.5543.0466.51982148.8554.5476.81983155.1567.1489.01984160.4575.9498.51985167.5586.4510.4基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实现单片机嵌入式以太网防盗报警系统基于51单片机的嵌入式Internet系统的设计与实现单片机监测系统在挤压机上的应用MSP430单片机在智能水表系统上的研究与应用