96基于古典概率模型的水资源短缺综合评价及其应用

1基于古典概率的水资源短缺风险综合评价模型及其应用摘要基于北京市水资源严重缺乏现状，本文利用SPSS软件采用主成分析法，找到水资源短缺风险敏感因子：降雨量、人口总数、年污水再生量、工业用水量、第三产业用水量、再生水。基于古典概率模型建立了水资源短缺风险综合评价模型，对水资源发生的概率和水资源缺水影响程度做了定量分析，得到水资源短缺风险综合值。首先，建立反映缺水影响程度的风险度的隶属函数；然后，利用古典概率，计算每年发生水资源短缺风险的概率为0.742，风险度和风险率的乘积反映风险综合程度，即风险综合值。利用k均值聚类法，将风险等级分为五类，得到五个聚类中心：0.01、0.12、0.31、0.46、0.66，整理每个聚类中心所对应的所有风险综合值，将风险等级按风险综合值分为五等：0-0.066为可以忽略的风险，0.066-0.2075为可以接受的风险，0.2075-0.4035为边缘风险，0.4035-0.536为比较严重的风险，0.536-1为无法承受的风险。利用时间序列分析模型，得到2011年及2012年农业用水、工业用水、第三产业及生活等其他用水、水资源总量、年污水再生量、再生水量、南水北调工程调水量七个数据，分别为（单位：亿立方米）：11.1154、3.7719、20.3191、21.74、13.1096、8.7175、2.6；10.7868、3.2148、21.2572、21.74、14.1832、9.6916、2.6，得到2011年和2012年的需水量分别为（单位：亿立方米）：35.2064、35.2588，供水量分别为（单位：亿立方米）：46.1671、48.2148，缺水量（单位：亿立方米）分别为：-10.9607、-12.956。根据建立的水资源短缺风险综合评价模型，得到2011年和2012年风险综合值均为0，风险等级为可以忽略的风险。降低水资源短缺风险，本文从减少需水量，增加供水量这两个角度提出相应建议。以水资源短缺风险敏感因子为主要突破口，采取措施降低风险，本问提出了减少工业用水、第三产业及生活等其他用水的具体措施，增加年污水再生量以及再生水利用量的具体措施，并提出“生态移民”这一重要措施，这也是向北京市水行政主管部门写的建议的总体思路和内容。关键词：水资源短缺风险；古典概率；隶属函数；时间序列分析模型；主成分分析2一问题重述近年来，我国，特别是北方地区水资源短缺问题日趋严重，水资源成为焦点话题。以北京市为例，其人均水资源占有量不足3003m，为全国人均的1/8，世界人均的1/30，属重度缺水地区，北京市水资源短缺已经成为影响和制约首都社会和经济发展的主要因素。现针对如何对水资源风险的主要因子进行识别，对风险造成的危害等级进行划分，对不同风险因子采取相应的有效措施规避风险或减少其造成的危害这三个方面解决下面四个问题：1.评价判定北京市水资源短缺风险的主要风险因子是什么？（影响水资源的因素很多,例如：气候条件、水利工程设施、工业污染、农业用水、管理制度，人口规模等。）2.建立一个数学模型对北京市水资源短缺风险进行综合评价，作出风险等级划分并陈述理由。对主要风险因子,如何进行调控，使得风险降低？3.对北京市未来两年水资源的短缺风险进行预测，并提出应对措施。4.以北京市水行政主管部门为报告对象，写一份建议报告。二问题分析要解决本题，应该针对各小问题逐个分析，逐步解决，最后建立模型解决本题。解决第一问的关键是归纳总结所有影响水资源短缺风险的因子，然后用SPSS软件采用主层次分析法找出其中的敏感因子，即水资源短缺风险的主要风险因子。第二问首先需要建立一个风险评价模型，将每年的风险进行量化，再用聚类分析法将所有风险分成若干等级，针对主要因子产生的原因，进行弥补，使风险降低。第三问是一个时间序列预测模型，根据已有数据，采用线性回归计算未来两年2011年、2012年的农业用水、工业用水、第三产业及生活等其他用水、水资源总量、再生水、年污水再生量、南水北调工程调水量7个值，从而求出风险度，进一步求出风险综合值，完成对风险等级的预测，根据预测结果提出相应应对措施。第四问以水资源短缺风险敏感因子为突破口，提出相应建议。三问题假设1.假设用完的水必须的经过处理后排放，即所产生的污水总量与所处理的污水总量相等，忽略随意乱排污现象。2.假设每年水资源短缺风险概率是一样的。3.年鉴中2001年-2009年再生水和南水北调的数据是从某一年开始的，假设在该年前这两个数据均为0。4.假设南水北调工程调水量在短期内保持不变。3四符号说明1x降水量2x日照时数3x水资源总量4x人口总数5x年污水再生量6x水产品价格指数7x工业用水8x农业用水9x第三产业及生活等其它用水量10x再生水11x南水北调工程调水量g供水量n需水量iq第i年的缺水量minq缺水情况下最少的缺水量maxq最大的缺水量i第i年的风险度if第i年的风险率iw第i年的风险综合值kj第k组风险综合值的第j个数据t年份4五模型的建立与求解1.主要风险因子的确定需要求评价判定北京市水资源短缺风险的主要风险因子，首先需要了解一个概念：水资源短缺风险。水资源短缺风险，泛指在特定的时空环境条件下，由于来水和用水两方面存在不确定性，使区域水资源系统发生供水短缺的可能性以及由此产生的损失。所以，本问可从来水和用水两个方面，从自然因素、社会经济、工程技术三个角度收集可能影响水资源短缺的各个因素，并将其量化，采用层次分析法确定主要风险因子。1.1主层次分析法的引入主成分分析：将原来较多的指标简化为少数几个新的综合指标的多元统计方法。主成分：由原始指标综合形成的几个新指标。依据主成分所含信息量的大小成为第一主成分，第二主成分等等。主成分分析法在医学、心理学、经济学等科学领域以及社会化生产中得到广泛的应用。对于本问，可以通过主成分分析得到主成份，各个变量的系数反应了其贡献率及影响力的大小。1.2收集与水资源短缺风险相关的各个因子确定水资源短缺风险因子，首先建立一个风险指标体系]1[，内容如下：水资源短缺风险指标体系：自然因素：降水量1x，见附表1日照时数2x，见附表1水资源总量3x，见附表2、3社会经济：人口总数4x,见附表4年污水再生量5x（5x=日污水处理能力*365*污水处理率），见附表5水产品价格指数6x,见附表6工业用水量7x,见附表2、3农业用水量8x,见附表2、3第三产业及生活等其它用水量9x,见附表2、35工程技术：再生水10x,见附表3南水北调工程调水量11x,见附表3年污水再生量计算如下：年污水再生量365**污水处理率日污水处理量整理附表1-6得各因素数据统计如表1：表1：各因子统计表(1979-2009)年份降雨量1x（mm)日照时数2x(时）水资源总量3x(亿立方米）人口总数4x(万人）年污水再生量5x（亿立方米）水产品价格指数6x工业用水量7x（亿立方米）1979718.42667.438.23897.10.0856107.214.371980380.72920.826904.30.789131.313.771981393.22803.924919.20.0986100.012.211982544.42825.136.6935.00.0995100.013.891983489.92844.334.7950.00.0931100.011.241984488.82767.639.31965.00.0913106.614.3761985721.02511.938981.00.0913235.317.21986665.32804.127.031028.00.0845129.79.911987683.92631.938.661047.00.0731116.114.011988673.32558.139.181061.00.0702144.914.041989442.22626.221.551075.00.0626114.513.771990697.32325.035.861086.00.0799101.712.341991747.92536.642.291094.00.0723102.011.91992541.52712.522.441102.00.0022101.815.511993506.72669.819.671112.00.0057109.315.281994813.22470.545.421125.00.0876125.614.571995572.52519.130.341251.10.4178111.113.781996700.92418.745.871259.40.4565104.511.7661997430.92596.522.251240.00.4738109.711.11998731.72420.737.71245.60.484594.810.841999266.92594.014.221257.20.538495.510.562000371.12667.216.861363.61.855107.310.522001338.92611.719.21385.12.20896.99.22002370.42588.416.11423.22.97395.77.52003444.92260.218.41456.43.932102.58.42004483.52515.421.41492.75.017107.17.72005410.72576.123.21538.07.379104.76.82006318.02192.724.51581.08.916101.96.22007483.92351.123.81633.09.818108.85.82008626.32391.434.21695.09.475120.15.22009480.62511.821.81755.010.43104.55.2续表1年份农业用水量8x第三产业及生活等其它用水量9x（亿立方米）再生水10x（亿立方米）南水北调工程调水量11x（亿立方米）197924.184.3700198031.834.9400198131.64.300198228.814.5200198331.64.7200198421.844.01700198510.124.3900198619.467.180019879.687.2600198821.996.4007198924.426.4500199021.747.0400199122.77.4300199219.9410.9800199320.359.5900199420.9310.3700199519.3311.7700199618.959.300199718.1211.100199817.3912.200199918.4512.700200016.4913.3900200117.412.300200215.511.600200313.813.62.10200413.513.42.00200513.214.52.60200612.815.33.60200712.416.65.00200812.017.96.00200912.018.36.52.6将所有数据标准化]1[：iiijjisxxy（1）ijx是第i个因子的第j个值，个因子的均值，是第ixiis第i个因子的标准差，整理得表2：8表2：各因子标准值（1979-2009）年份降雨量1x日照时数2x水资源总量3x人口总数4x年污水再生量5x水产品价格指数6x19791.26420.52501.0005-1.3125-0.6134-0.21231980-1.04491.9985-0.3242-1.2834-0.40310.73251981-0.95941.3187-0.5408-1.2230-0.6095-0.494619820.07441.44200.8240-1.1591-0.6092-0.49461983-0.29821.55360.6182-1.0983-0.6111-0.49461984-0.30571