【毕业学位论文】（Word原稿）北京市水资源短缺风险关系数学模型-统计教育学

北京市水资源短缺风险关系 数学模型 北京航空航天大学 摘要 题目要求 分析北京市资源短缺风险，对其进行综合评定及相关预测，并提出合理的调控措施，降低风险。 鉴于影响因素众多， 但相关 程度各不相同，首先我们借用 学软件对各项影响因素进行主成成分分析，找出其中相关系数较大的因素，即为主要风险因子。 其次， 为了更好的刻画水资源利用情况，我们引进一个概念 水资源短缺程度，并将其定量表示为： Z=（用水量 总水量） /用水量。 水资源短缺 带来的风险， 将直接反映在经济增长上，因此 我们主要研究水资源利用情况给经 济 带来的影响 。 利用所给数据在 模拟出 历年来 工业 、农业、第三产业的经济增长状况与 Z 值的关系，并拟合出对应的函数表达式。结合图形，以经济增长量为参考基准，对风险因子进行等级划分。 根据已知的数据并对未来几年的水资源短缺程度和对经济的影响程度进行预测，指出未来可能遇到的水资源风险， 最后提出了合理的调控措施以降低风险，同时，对模型进行了一定的评价与讨论。 关键词： 水资源短缺风险、经济增长、主要因子、风险预测 1 问题 描述 题目中提供了 1978 2008 年间北京市的水资源利用情况， 其中包括总用水量、工 业用水、农业用水、第三产业用水和生活用水、水资源总量等信息。 要求根据所给数据及其他相关资料解决如下问题： 1） 从众多影响北京市水资源利用的因素中 分离 判定 出北京市水资源短缺风险的主要风险因子； 2） 建立合适的数学模型对北京市水资源短缺风险进行综合评价， 作出风险等级划分 ； 3） 对主要风险因子进行具体分析 ,提出合理的调控措施，使得风险降低。 4） 对北京市未来两年水资源的短缺风险进行预测，并提出应对措施。 以北京市水行政主管部门为报告对象，写一份建议报告。 2 问题分析 水资源短缺风险，泛指在特定的时空环境条件下 ，由于来水和用水两方面存在不确定性，使区域水资源系统发生供水短缺的可能性以及由此产生的损失。 我们要判定水资源短缺风险，需先对水资源短缺程度进行分析，进而由短缺所造成的经济损失来判定水资源短缺风险。 我们可以以水资源短缺程度为目标函数，对可能影响水资源短缺程度的因子进行主要因子分析，从而得出影响水资源短缺的主要因子，并根据北京市的经济增长数据对风险因子的作用进行评估。对未来的预测需要对已知的数据信息进行处理，在没有更多外在影响的情况下，可以得出未来几年的水资源短缺程度。 根据预测的信息，我们可以对主要风险因子 进行调控，从而可以向北京市水行政主管部门写一份建议报告。 3 假设与符号 本假设 1）北京的变化是连续与光滑的 ； 2）没有火灾、地震等重大事情发生 ； 3）北京的政策没有大的变动 ； 4）其他因素的影响是随机与对称的 ； 号说明 4 模型的建立与求解 要风险因子 年份 总用水量(亿立方米 ) 农业用水(亿立方米 ) 工业用水 (亿立方米 ) 第三产业及生活等其它用水(亿立方米 ) 水资源总量（亿方） 年降水量（毫米） 人口规模 (万人） 水资源短缺程度 (%) 1979 980 6 981 4 982 35 983 50 984 65 985 8 721 981 986 028 987 047 988 061 989 075 990 086 991 094 992 102 993 112 994 125 995 996 997 240 998 999 000 001 002 003 004 005 538 006 18 1581 007 633 008 2 695 产业经济 增长率 与水资源短缺程度 之间的关系如下表所示： 年份 农业产值增长率 工业产值增长率 第三产业增长 率 总产业增长率 目标函数 Z 1980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 000 001 002 003 004 005 006 007 008 据生活经验，我们筛选出上述 六 个影响水资源利用的因素进行主成成分分析。 我们采用 主成分 分析法来分析影响水资源短缺程度的主要因子，在进行主因子分析的过程中，我们是对每一个参数进行过参数标注化处理，从而参数的单位不会对结果造成影响。 g ; n,m= N=m; y=:,N); y=(y)/(y)y); :,1: 1,); 1,); i=2:i*1,); i*1,); ; pc 00* pc*; ; 主成分分析 ); 这是主成分分析所得的权重图。其中横坐标表示实际因子进行空间转化后的因子，纵坐标表示百分比。图中只列出了对目标参数影响之和超过 99%的几个参数。其中的直方图表示，每个 目标参数的影响权重，其中的折线表示，前几个 主成分分析的结果如下： 表 PC 业用水 业用水 三产业用水 资源 降雨量 口规模 中横坐标表示主因子分析中空间转后的因素，从坐标表示每个实际具体因子在空间转换后因素上的权重。 表 主因子分析中 2 3 4 5 6 百分比 % 中目标参数就是我们的缺水度。 由上述两 个表格，我们可以得到实际因子对缺水度影响的权重 项目 农业用水 工业用水 第三产业用水 水资源总量 年降水量 人口规模 相关系数 表格可以看出， 影响水资源短缺程度的主要 风险 因子有农业用水， 工业用水、第三产业及人口规模等 ， 这些因素与我们在感性上的认识是相同的。 险等级划分 用 合曲线如下： 图示中，横坐标为缺水程度，纵坐标为各产业经济增长值。 同时，我们对每一个产业进行风险划分，当水资源短缺率为 0 以 下时，我们认为风险很低，当经济增长率岁水资源短缺的变化发生很大 变化时，为风险较大的区域。同时对风险低和风险高之间进行细分，认为风险较低和风险较高。 说明：在绘图过程中大部分点相对集中，我们略去少数偏离较远的点。因为这些点所代表的增长率不仅与水有关，还在很大程度上受到了其他非相关因素，如经济政策、自然条件等的影响，因此应该不予考虑。 表 1 历年来总产值增长率 年份 1980 1981 1982 1983 1984 1985 1986 1987 增长率 水度 Z 份 1988 1989 1990 1991 1992 1993 1994 1995 增长率 水度 Z 份 1996 1997 1998 1999 2000 2001 2002 2003 增长率 水度 Z 份 2004 2005 2006 2007 2008 增长率 缺水度 Z 年份 1980 1981 1982 1983 1984 1985 1986 1987 1 历年来农业生产总值增长率 年份 1980 1981 1982 1983 1984 1985 1986 1987 增长率 水度 Z 份 1988 1989 1990 1991 1992 1993 1994 1995 增长率 年份 1980 1981 1982 1983 1984 1985 1986 1987 份 1988 1989 1990 1991 1992 1993 1994 1995 份 1996 1997 1998 1999 2000 2001 2002 2003 份 2004 2005 2006 2007 2008 Z 图 2 农业 产值 缺水度 缺水度 Z 份 1996 1997 1998 1999 2000 2001 2002 2003 增长率 水度 Z 份 2004 2005 2006 2007 2008 增长率 缺水度 Z 表 1 历年来工业生产总值增长率 年份 1980 1981 1982 1983 1984 1985 1986 1987 增长率 水度 Z 份 1988 1989 1990 1991 1992 1993 1994 1995 增长率 水度 Z 份 1996 1997 1998 1999 2000 2001 2002 2003 增长率 水度 Z 份 2004 2005 2006 2007 2008 增长率 缺水度 Z 图 3 工业 缺水度 表 1 历年来第三产业生产总值增长率 年份 农业用水 工业用水 第三产业用水 人口规模 缺水程度 1999 000 001 002 003 004 005 538 006 581 007 633 3 第三产业 缺水度 从图 1 中可以很明显地看出缺水对于经济总量 影响，当 Z 值小于 0 时，曲线基本呈 持平趋势，几乎与横轴平行，说明水分充足时，不对经济增长造成不良影响；而一旦 Z 值大于 0 时，总产值 增长率 明显呈下降趋势，表明缺水会对经济增长形成相应的阻碍作用，且缺水越严重，阻碍作用越强烈。 然而，单独观察图 2、 3、 4，又可 看出，缺水度对各项产值的影响有着很大区别。分析如下： 比较各图可知，缺水度对于农业产值影响最大，当 Z 值大于 0 时 ，农业产值增长率立即呈现 明显的下降趋势。之所 以 出现这种情况，是因为农业上需要大量的水进行灌溉，而 由于我国 地广人多的基本情况， 先进的 农业灌溉技术 ，如喷灌、滴灌等 并没有完全普及 ，因此在灌溉作物时，大部分水得不到合理利用，也无法回收。因此，农业对水的需求量是非常大的。一旦水资源出现匮乏，必然会使农业产值的增长严重滞后。 相比较而言，工业对于水资源的依赖程度要相对轻一些。当缺水度 Z60%时，工业产值基本 保 持不变的增长率。联系我们国家对于排污治污的各项严格规定就不难解释这个问题了。工业用水一般都经过治理达标才可排放，以达到可重复利用的目的。 第三产业产值增长情况则介于上述两者之间。 f(x) = a*b*x) + c*d*x) 农业 a = b = c = d = 业 a = b = c = d = 第三产业 a = 5720 b = c = d = 总产值 a = b = c = d = 未来四年水资源风险的预测 灰度预测程序如下： x) g if x (:,1)=1 x=x; n=x); z=0; i=1:n z=z+x (i,:); i,:)=z; i=2:n y ()=x (i,:); i=1:c (i,:)=i,:)+i+1,:); j=1:e (j,:)=1; i=1: (i,1)=c (i,:); B (i,2)=e (i,:); B*B)*B*y; i=1:n+1 i,:)=(x (1,:)2,:)/,:)*,:)*(+2,:)/,:); 1,:)=1,:); i=1:i+1,:)=i+1,:)i,:); 00*x; 1、对农业用水进行灰度预测 表 业用 水量灰度预测 年份 实际农业用水量 预测值 绝对误差 相对误差 % 1979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 000 001 002 003 004 005 006 007 008 12 预测的结果可以看出误差很大 ，从而 应对最近几年的进行预测： 表 次预测 年份 实际农业用水量 预测值 绝对误差 相对误差 % 2001 002 003 004 005 006 007 008 12 预测结果课果可以看出，现在的预测程序是可信的，对农业用水进行预测可得： 表 预测结果 年份 2009 2010 2011 2012 预测值 、对工业用水量的灰色预测 表 业用水量灰度预测结果 年份 实际工业用水量 预测值 绝对误差 相对误差 % 1979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 000 001 002 003 004 005 006 007 008 以看出最后几年的相对误差很大，这是 由于科学技术的进步 ，从而 每年的用水量会有很大的不同 。 因而不能对改革开放三十年以来所有的进行灰色预测 ，于是对最近十年进行灰色预测： 表 次与预测 年份 实际工业用水量 预测值 绝对误差 相对误差 % 1999 000 001 002 003 004 005 006 007 008 从预测结果课果可以看出，现在的预测程序是可信的，对 工 业用水进行预测可得： 表 测结果 年份 2009 2010 2011 2012 预测值 、对第三产业用水量的灰色预测 表 三产业用水量灰色预测 年份 实际第三产业 预测值 绝对误差 相对误差 % 1979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 000 001 002 003 004 005 006 007 008 以看出最后几年的相对误差很大， 这是由于三十年来，中国社会巨大的变化造成的，从而 不能对改革开放三十年以来所有的 数据 进行灰色预测 ，于是 对最近七 年进行灰色预测： 表 次预测 年份 实际第三产业 预测值 绝 对误差 相对误差 % 2002 003 004 005 006 007 008 预测结果课果可以看出，现在的预测程序是可信的，对第三产业用水进行预测可得： 表 测结果 年份 2009 2010 2011 2012 预测值 、对人口规模的灰色预测 表 口规模灰色预测 年份 实际人口规模 预测值 绝对误差 相对误差 % 1979 980 981 982 935 983 950 984 965 985 981 986 1028 987 1047 988 1061 989 1075