毕业设计（论文）-水资源短缺风险综合评价.doc

《毕业设计（论文）-水资源短缺风险综合评价.doc》由会员分享，可在线阅读，更多相关《毕业设计（论文）-水资源短缺风险综合评价.doc（32页珍藏版）》请在三一文库上搜索。

1、水资源短缺风险综合评价 摘要：针对北京市水资源短缺的农业用水等九大主要因素，用熵权法得出人口数量和降水量的影响因素较大。本文选取区域水资源短缺风险程度的风险率、脆弱性、可恢复性、重现期和风险度作为评价指标，构建了模糊综合评价模型，结论表明北京市水资源短缺现处于高风险状态，并建立多元线性回归和灰色系统GM模型，预测北京市未来两年水资源短缺仍将持续处于高风险状态。根据所建模型及预测结果向相关部门提出控制在京人口以及合理分配农业、工业、生态用水量来缓解北京水资源短缺现状。一、问题重述1.1问题的提出水资源，是指可供人类直接利用，能够不断更新的天然水体。主要包括陆地上的地表水和地下水。风险，是指某一特

2、定危险情况发生的可能性和后果的组合。水资源短缺风险，泛指在特定的时空环境条件下，由于来水和用水两方面存在不确定性，使区域水资源系统发生供水短缺的可能性以及由此产生的损失。近年来，我国、特别是北方地区水资源短缺问题日趋严重，水资源成为焦点话题。以北京市为例，北京是世界上水资源严重缺乏的大都市之一，其人均水资源占有量不足300m3，为全国人均的1/8，世界人均的1/30，属重度缺水地区，附表中所列的数据给出了1979年至2000年北京市水资源短缺的状况。北京市水资源短缺已经成为影响和制约首都社会和经济发展的主要因素。政府采取了一系列措施, 如南水北调工程建设, 建立污水处理厂,产业结构调整等。但是

3、，气候变化和经济社会不断发展，水资源短缺风险始终存在。如何对水资源风险的主要因子进行识别，对风险造成的危害等级进行划分，对不同风险因子采取相应的有效措施规避风险或减少其造成的危害，这对社会经济的稳定、可持续发展战略的实施具有重要的意义。北京2009统计年鉴及市政统计资料提供了北京市水资源的有关信息。利用这些资料和你自己可获得的其他资料，讨论以下问题：1 评价判定北京市水资源短缺风险的主要风险因子是什么？影响水资源的因素很多,例如：气候条件、水利工程设施、工业污染、农业用水、管理制度，人口规模等。2建立一个数学模型对北京市水资源短缺风险进行综合评价， 作出风险等级划分并陈述理由。对主要风险因子,

4、如何进行调控，使得风险降低？ 3 对北京市未来两年水资源的短缺风险进行预测，并提出应对措施。 4 以北京市水行政主管部门为报告对象，写一份建议报告。1.2 问题简要分析为了能客观准确评价判定北京市水资源短缺风险的主要风险因子，首先找出可能造成北京水资源短缺的多个平行因素，依据历年的耗水指标，根据熵值法确定各因素的权重系数，找出熵值最小（即权重最大）的指标就是水资源短缺风险最主要的风险因子。在实际解决北京水资源短缺的问题时就可以从该主要因素入手。要在风险识别和风险分析的基础上进行水资源短缺风险做综合评价和风险等级划分，把损失概率、损失程度以及其它因素综合起来考虑，分析该风险的影响，寻求风险对策并

5、分析该对策的影响，为风险决策创造条件。通过对定义的风险率、脆弱性、可恢复性、重现期、风险度作为水资源短缺风险的评价指标的计算，采用模糊综合评判方法对水资源短缺风险进行评价，判断风险等级。接着用层次分析法作出评价指标的权重，通过加权平均型模型的建立，计算出风险性的综合评分。对北京市未来两年水资源的短缺风险进行预测，建立灰色系统理论模型和多元线性回归模型两种不同的模型对水资源短缺风险分别进行预测，在假设数据相同计算无误的基础上，该两种模型对未来两年北京水资源的短缺风险的预测结果应该是相同的。二、基本假设和符号说明2.1基本假设1、影响北京水资源短缺的多个因素在未来没有突变情况发生（政府政策的干预，

6、自然灾害的影响等）。2、假设影响北京水资源短缺的多个因素相互独立，即这个指标对风险的影响与其他指标对风险的影响没有关系。3、2.2 符号说明序号符号含义1评价指标的熵权2风险率性能衡量指标3脆弱性性能衡量指标4重现期性能衡量指标5可恢复性性能衡量指标6风险度性能衡量指标7风险性的综合评分值8水资源差值时间序列9水资源差值时间序列级比10模糊关系矩阵11相对于因素 而言属于 vi 的程度12残差13级比偏差值三 模型建立与求解3.1问题一 -P。在本题中，我们在提供的气候条件、水利工程设施、工业用水、农业用水、第三产业及其他用水、人口规模几个影响因子基础上还增加了生态用水量、降水量、园林绿化覆盖

7、率、废水治理设施数、工业排放总量等影响因素。通过对1979年至2009年各个因素数据的分析，构建模糊物元，我们采用了运用熵值法来确定权重系数。由于熵值越小，表明指标值的变异程度越大，提供的信息量越多，在综合评价中所起的作用越大，其权重也越大，故我们最后可以通过各个因素的熵值大小，比较得出各个因子的权重大小，从而得到影响北京水资源短缺的主要风险因子。3.1.1 构建模糊物元模糊物元及复合模糊物元，在物元分析中所描述的事物M及其特征C和量值x组成物元R=(M，C,X)或者R=(M，C，C(M)，同时把事物的名称、特征和量值称为物元三要素。如果物元模型中的量值x具有模糊性，便称其为模糊物元。事物M有

8、n个特征及其相应的量值，则称R为n维模糊物元。m个事物的n维物元组合在一起便构成m个事物的n维复合模糊物元，即 式中：为m个事物的n个模糊特征的复合物元；为第i个事物(i=1,2,m)；为第j个特征(j=1,2,n);为第i个事物第j个特征对应的模糊量值。3.1.2 熵值法确定权重系数在确定评价指标的权重时，往往多采用主观确定权重的方法，如AHP方法等。这样就会造成评价结果可能由于人的主观因素而形成偏差。在信息论中，熵值反映了信息无序化程度，其值越小，系统无序度越小，故可用信息熵评价所获系统信息的有序度及其效用，即由评价指标值构成的判断矩阵来确定指标权重，它能尽量消除各指标权重计算的人为干扰，

9、使评价结果更符合实际。其计算步骤如下：（1）构建m个事物n个评价指标的判R=()nm(i=1,2,n;j=1,2,m)。（2）将判断矩阵归一化处理，得到归一化判断矩阵B 式中：、分别为同指标下不同事物中最满意者或最不满意者(越小越满意或越大越满意)。（3）根据熵的定义，m个评价事物n个评价指标，可以确定评价指标的熵为： 为使有意义，当 =0时，根据风险评价的实际意义，可以理解为一较大的数值，与相乘趋于0，故可认为=0。但当=1，也等于0，这显然与熵所反映的信息无序化程度相悖，不切合实际，故需对加以修正，将其定义为:(4)计算评价指标的熵权W ,且满足3.1.3 权重的求解及主要风险因子的确定我

10、们通过3.1.2中涉及的熵值求权重的方法，运用MATLAB编程得到了各个因素的熵值，通过比较可以得到影响北京市水资源短缺的主要风险因子。熵值法的主要运算结果见下表： 表一：各个风险因子归一化后的年份农业用水生态用水工业用水其他用水人口数量降水量园林覆盖率治理设施工业污染20000.36370.25350.4110.24040.17210.19910.14360.37180.306820010.39350.25260.360.193280.18290.25680.14760.31980.296420020.33130.260.29430.15260.20220.2390.19310.30380.

11、28520030.27550.26190.3310.22720.2190.25640.22880.28040.280920040.26560.26190.30210.22040.23740.29770.23480.26790.266720050.25580.27130.26730.24080.26040.31910.25550.26380.272320060.24270.31810.24410.25090.28210.27880.25730.21790.195620070.22960.4210.22870.25770.30850.22740.26730.25580.165620080.2164

12、0.46780.20550.28480.340.31930.27720.18590.157920090.19680.50510.20550.30520.34420.39820.28720.20590.1534 表二：各个风险因子的熵以及熵权农业用水生态用水工业用水其他用水人口数量降水量园林覆盖率治理设施工业污染H0.13890.14980.14120.15730.15940.15920.15560.14130.1249W0.11220.11080.11190.10980.10960.10960.11010.11190.1141 由上表可知：农业用水、生态用水，工业用水、其他用水，人口数量、降水

13、量、园林覆盖率、治理设施、工业污染的熵权值分别为：0.1122、0.1108、0.1119、0.1098、0.1096、0.1096、0.1101、0.1119、0.1141。由于熵值越小影响越大，其权重也越大的原则，我们可以很清楚的知道：在影响北京水资源短缺的众多因素中，北京市人口总数以及北京市降雨量对其影响程度较大。3.2 问题二 在问题二中我们要对北京市水资源短缺风险进行综合分析评价，并作出风险等级划分。水资源的短缺取决于供水和需水两方面影响，而这两方面都具有随机性和不确定性。因此水资源短缺风险也具有随机性和不确定性。我们在进行风险评价时，充分考虑风险特点以及水资源的复杂性，把存在风险的

14、概率、风险出现的时间、风险造成的损失多少、风险解除时间、缺水量的分布等一系列因素考虑在内。故我们从多方面的指标综合考虑评价北京市现水资源短缺风险等级。我们在评价指标选择中坚持的原则：（1）能集中反映缺水的风险程度；（2）能集中反映水资源短缺风险发生后水资源系统的承受能力；（3）代表性好，针对性强，易于量化。依据上述原则，我们选取水资源的风险率、脆弱性、可恢复性、事故周期、风险度作为水资源系统水资源短缺风险的评价指标。在水资源短缺风险评价等级模型建立过程中，我们采用模糊概率理论，建立了模糊概率模型，通过模型我们可以清晰地看到北京市现水资源风险等级。3.2.1 水资源短缺风险评价指标1. 风险率

15、根据风险理论，荷载是使系统“失事”的驱动力，而抗力则是对象抵御“失事”的能力。如果把水资源系统的失事状态记为 F()，正常状态记为 S()=PF式中：为水资源系统状态变量如果水资源系统的工作状态有长期的记录，风险率也可以定义为水资源系统不能正常工作的时间与整个工作历时之比，即式中：NS 为水资源系统工作的总历时；It是水资源系统的状态变量。 2 脆弱性 脆弱性是描述水资源系统失事损失平均严重程度的重要指标。为了定量表示系统的脆弱性，假定系统第 i 次失事的损失程度为 Si，其相应的发生概率为 Pi，那么系统的脆弱性可表达为：式中：NF为系统失事的总次数。例如，在供水系统的风险分析中，可以用缺水

16、量来描述系统缺水失事的损失程度。类似洪水分析，假定 P1=P2=PNF=1/NF，即不同缺水量的缺水事件是同频率的，这样上式可写为:式中：为第 i 次缺水的缺水量。上式说明干旱的期望缺水量可以用来表示供水系统的脆弱性。为了消除需水量不同的影响，一般采用相对值，即式中：是第 i 次干旱缺水期的需水量。3重现期 事故周期是两次进入失事模式 F之间的时间间隔，也叫平均重现期。用 表示第 n间隔时间的历时，则平均重现期为式中：N=N()是 0 到 t 时段内属于模式 F 的事故数目。4可恢复性 恢复性是描述系统从事故状态返回到正常状态的可能性。系统的恢复性越高，表明该系统能更快地从事故状态转变为正常运

17、行状态。它可以由如下的条件概率来定义：上式亦可用全概率公式改写为引入整数变量及 这样，由全概率公式可得 记 则有：从上式可以看出，当 TF=0，即水资源系统在整个历时一直处于正常工作状态，则=1；而当 TFS=0，即水资源系统一直处于失事状态(TF=NS)，则=。一般来讲，01。这表明水资源系统有时会处于失事状态，但有可能恢复正常状态，而且失事的历时越长，恢复性越小，也就是说水资源系统在经历了一个较长时期的失事之后，转为正常状态是比较困难的。5风险度 用概率分布的数学特征，如标准差或半标准差-，可以说明风险的大小。和-越大，则风险越大，反之越小。这是因为概率分布越分散，实际结果远离期望值的概率

18、就越大 。 用、-比较风险大小虽简单，概念明确，但-为某一物理量的绝对量，当两个比较方案的期望值相差很大时，则可比性差，同时比较结果可能不准确。为了克服用-可比性差的不足，可用其相对量作为比较参数，该相对量定义为风险度，即标准差与期望值的比值(也称变差系数)。这里值得说明的是：风险度不同于风险率，前者的值可大于 1，而后者只能小于或等于 1。3.2.2 水资源短缺风险的模糊综合评判模型的建立风险评价是在风险识别和风险分析的基础上，把损失概率、损失程度以及其它因素综合起来考虑，分析该风险的影响，寻求风险对策并分析该对策的影响，为风险决策创造条件。本文采用上述定义的风险率、脆弱性、可恢复性、重现期

19、、风险度作为水资源短缺风险的评价指标，采用模糊综合评判方法对水资源短缺风险进行评价。设给定 2 个有限论域 和其中，U 代表综合评判的因素所组成的集合；V 代表评语所组成的集合。则模糊综合评判即表示下列的模糊变换 B=AR，式中 A 为 U 上的模糊子集。而评判结果 B 是 V 上的模糊子集，并且可表示为 ，01；，01。其中表示单因素 在总评定因素中所起作用大小的变量，也在一定程度上代表根据单因素评定等级的能力；为等级对综合评定所得模糊子集 B 的隶属度，它表示综合评判的结果。关系矩阵 R 可表示为 式中：表示因素 的评价对等级的隶属度，因而矩阵R中第即为对第i个因素 的单因素评判结果。在评

20、价计算中代表了各个因素对综合评判重要性的权系数，因此满；同时，模糊变换 AR 也即退化为普通矩阵计算，即上述权系数的确定可用层次分析法(AHP)得到。由上述分析可以看出，评价因素集对应评语集而评判矩阵中即为某因素对应等级的隶属度，其值可根据各评价因素的实际数值对照各因素的分级指标推求。我们将评语级分为 5 个级别，各评价因素分级指标见表 1。 表1 各评价因素分级指标水资源短缺风险(风险率)(脆弱性) (可恢复性)(重现期)(风险度) (低) 0.200 0.200 0.800 9.000 0.200 (较低)0.2000.4000.2000.4000.6010.8006.0019.0000.

21、2010.600(中)0.4010.6000.4010.6000.4010.6003.0016.0000.6011.000(较高)0.6010.8000.6010.8000.2000.4001.0003.0001.0012.000(高) 0.800 0.800 0.200 1 2.000由于水资源风险率、脆弱性、风险度是越小越优性指标，所以对于 , 各评语级可构造如下隶属函数： 由于水资源可恢复性和重现期是越大越优性指标，所以对于 、 各评语级可构造如下隶属函数： 对于水资源短缺风险评价的因素集U而言，对应一个测定指标向。 其中 是 U 相对于 的测定值。这样 便表示相对于因素 而言属于 vi

22、 的程度。对于因素集 U，便有下面的模糊关系矩阵： 水资源短缺风险评价各因素的权重确定采用层次分析法(AHP)，设权重计算结果为，于是可得出综合评判向量 在综合评判中，我们选取“加权平均型”的 M (,) 模型，即 由于该模型实际上蜕化为一般的实数加法，即 选取与对应的评语为区域水资源短缺风险的评判结果。为了比较直观的说明风险程度，我们将其分成 5 级，分别叫做低风险、较低风险、中风险、较高风险和高风险，风险各级别按综合分值评判，其评判标准和各级别风险的特征见下表2。表 2 水资源系统水资源短缺风险级别评价水资源短缺风险评价等级风险级别水资源系统的风险特征低风险可以忽略的风险较低风险可以接受的

23、风险中风险边缘风险较高风险不可接受风险高风险灾变风险，系统受到严重破坏3.2.3 模糊综合评价模型的求解 根据五种评价因素各自的定义，我们对北京市水资源总量和总用水量从1979年到2009年的数据进行分析，得到北京市这31年水资源综合的风险率、脆弱性、可恢复性、重现性、风险度的综合性能数值，具体见下表4： 表4：北京市1979年至2009年水资源评价因素综合性能数值风险率脆弱性可恢复性重现性风险度87.1%0.3360.10710.15395.05从表中数据可知道：=0.871，=0.336，=0.1071，=0.1539，=5.05，从而得到测定指标向量。根据各个因素的隶属函数得到因素对应等

24、级的隶属度，从而得到模糊关系矩阵：在计算水资源短缺风险评价各因素的权重时，我们采用层次分析法(AHP)，权重计算结果为：。 由模糊关系矩阵和风险评价各因素的权重，我们可以得到综合评价向量： 我们再根据“加权平均型”的 M (,) 模型，即 ，选取与对应的评语为区域水资源短缺风险的评判结果，具体评判结果见下表： 表6：北京市水资源短缺风险综合评价分值综合评价北京市0.062720.21000.10840.88820.9076高风险由北京市水资源短缺风险综合评价分值可知道：北京市水资源短缺风险已经达到了高风险程度，需要采取及时有效的方法进行改善。3.2.4 对主要风险因子的调控 由熵值法判断权重知

25、影响北京水资源短缺的主要因素为北京市人口，北京市降水量。人口因素为人为可控因素，降水量为自然因素为非可控因素，下面主要讨论改变可控因素的措施。人口因素：从历史和今后的走向看，虽然2002-2009一期南水北调等水利工程对京供水量的增加，但随着北京市人口的快速增加，2010年调查显示常驻人口十年内增加600余万远高于全国平均水平，加之高流动人口。水的承载力并不高，人口对水资源的承载压力是十分巨大，水利工程只能解决一时的问题，难以从根本上解决水资源紧张的问题。解决方法一:适当控制人口增长1调整功能，北京是政治中心和文化中心,北京的多项城市功能高度叠加在这些节点和区域，由此造成北京地区人口高度集中,

26、 要调控北京人口，一个重要的手段，就是要调整北京的功能定位，至少在人口规模上不能够进一步扩充，使市人口得到控制。2区域互动，推动北京周边城市和区域的均衡发展，将会在很大程度上分担区域流动人口的压力，通过自身功能作用的发挥，影响和辐射周边城市和区域，进而缩小与周边城市在经济发展上的差距，形成有区域经济特点的发展模式，减弱周边地区人口进京的动力。3高端发展，就是走科技创新之路，大力发展高端产业。通过调整经济结构和经济规模改变人口结构和人口规模，加大对高端人才引进力度的同时，减少对低端劳动力的需求，达到减少进京人口的目的。4产业转移，经济的高度密集是人口高度集中的根本原因。通过产业转移，调整经济结构

27、在空间上的分布，引导人口形成合理的空间分布，是调控人口的重要手段。解决方法二：采用资源集约型的经济增长方式和结构，提高水的利用率，大力推行节约用水，减少人均用水量。降水量因素：通过人类活动如植树造林、修建水库、人工造湖等，一般都是通过改变人为条件而间接增加降水， 以及必要时进行人工降水作业增加北京市的降雨量。3.3 问题三 在对北京市未来两年水资源的短缺风险进行预测，我们建立了两个模型来进行预测多元线性回归模型和灰色系统GM（1,1）模型。多元线性回归从多个方面考虑各种因素对水资源短缺风险的影响，通过线性拟合得到各个因素影响下供求差值的走势图，；灰色系统GM(1,1)模型直接通过水资源供求差值

28、入手，对数据进行相关处理，从而得到我们所需的预测模型。最后还通过残差和级比偏差来对建立的模型进行检验。通过检验，所建模型均可用来对北京市未来两年水资源短缺风险进行预测。多元线性回归模型（一）多元线性回归模型的概念在许多实际问题重中，我们所研究的因变量的变动可能不仅与一个解释变量有关。因此，有 必要考虑线性模型的更一般形式，即多元线性回归模型： 在这个模型中，Y由X1,X2,X3,XK所解释，有K+1个未知参数。这里，“斜率”j含义是其他变量不变的情况下，Xj改变一个单位对因变量所产生的影响。回到一般模型 即对于n组观测值，有 .其矩阵形式为：其中 二多元线性回归模型的估计多元线性回归模型的估计

29、与双变量线性模型类似，仍采用最小二乘法。当然，计算要复杂得多，通常要借助计算机。理论推导需借助矩阵代数。下面给出最小二乘法应用于多元线性回归模型的假设条件、估计结果及所得到的估计量的性质。1. 假设条件（1） E（ut）=0， t=1,2,n（2） E（uiuj）=0， ij （3） E（ut2）=，t=1,2,n（4） Xjt是非随机量，j=1,2,k t=1,2,n除上面4条外，在多个解释变量的情况下，还有两个条件需要满足：（5） （K+1）n； 即观测值的数目要大于带估计的参数的个数（要有足够数量的数据来拟合回归线）。（6）各解释变量之间不存在严格的线性关系。上述假设条件可用矩阵表示为一

30、下四个条件：（1）（2）由于显然， 仅当E（uiuj）=0， ij E（ut2）=，t=1,2,n这两个条件成立时才成立，因此，此条件相当前面条件（2），（3）两条，即各期扰动项互不相关，并具有常数方差。（3）是一个非随机元素矩阵。（4）Rank（）=（K+1）n.-相当于前面（5）、（6）两条即矩阵的秩=（K+1）n 当然，为了后面区间估计和假设检验的需要，还要加上一条：（5）2.最小二乘估计我们的模型是： 问题是选择，使得残差平方和最小。残差为：要使残差平方和为最小，则应有： ， ，.， 我们得到如下K+1个方程（即正规方程）： 按矩阵形式，上述方程组可表示为： 即 上述结果，亦可从矩阵表

31、示的模型 出发， 完全用矩阵代数推导出来。 残差可用矩阵表示为： 其中： 残差平方和注意到上式中所有项都是标量，且故令 用矩阵微分法，我们可得到 与采用标量式推导所得结果相同。由上述结果，我们有三多元线性拟合由20002009年北京农业用水，生态用水，工业用水，其他用水，人后数量，降水量，园林覆盖率，治理设施，工业污染九个影响因子的数据。如下表1：年份农业用水（亿立方米）生态用水（亿立方米）工业用水（亿立方米）其他用水（亿立方米）人口数量（万人）降水量（毫米）园林覆盖率（%）治理设施（套）工业污染（万吨）200016.490.9110.5213.391363.8266.936.30607140

32、01200117.40.99.212.01385.1371.136.5058113643.0200215.50.987.510.81423.2338.938.7857313250200313.818.4131456.4370.440.5756813107200413.517.712.81492.7444.940.8755512617200513.21.16.813.41538.0483.541.9155312813200612.81.66.213.71581.0410.742.0053010170200712.42.75.813.91633.0318.042.5054991342008123.

33、25.214.71695.0483.943.005148867200911.43.65.215.31703.2626.343.505248713如下表2：2000200120022003200420052006200720082009水资源总量(亿立方米)16.8619.216.118.421.423.224.523.834.221.8总用水量(亿立方米)40.439.834.635.834.634.534.334.835.135.5以以上9个因素为自变量，以北京市“总用水量”和“水资源总量”的差Y作为因变量，建立一个多元线性回归模型。由表2得到因变量Y的数组：y=23.54 20.6 18.

34、5 17.4 13.2 11.3 9.8 11 0.9 13.7由表1得到自变量x1，x2，x3，x4，x5，x6，x7，x8，x9的9个数组：x1=16.49,17.4,15.5,13.8,13.5,13.2,12.8,12.4,12,11.4x2=0.91,0.9,0.98,1,1,1.1,1.6,2.7,3.2,3.6x3=10.52,9.2,7.5,8.4,7.7,6.8,6.2,5.8,5.2,5.2x4=13.39,12.0,10.8,13,12.8,13.4,13.7,13.9,14.7,15.3x5=1363.6,1385.1,1423.2,1456.4,1492.7,1538

35、.0,1581.0,1633.0,1695.0,1703.2x6=266.9,371.1,338.9,370.4,444.9,483.5,410.7,318.0,483.9,626.3x7=36.30,36.50,38.78,40.57,40.87,41.91,42.00,42.50,43.00,43.50x8=607,581,573,568,555,553,530,549,514,524x9=14001,13643,13250,13107,12617,12813,10170,9134,8867,8713将y矩阵进行转置得到 增添一组常数项x0=0,0,0,0,0,0,0,0,0,0 将x=x

36、0,x1,x2,x3,x4,x5,x6,x7,x8,x9转置得到由模型，用矩阵微分法得到，则所以通过matlab进行矩阵运算得到即得到多元拟合线性方程四一元线性拟合回归将如表1所示的9个因子数据组成的是个向量组x1,x2,x3,x4,x5,x6,x7,x8,x9,与因变量y的一元线性拟合，得到以上9个因子与缺水量y的一元线性关系。因为9个因子的单位各不相同，所以统一将各个自变量的进行归一化式中：、分别表示为同因素下最小值和最小值。得到9个行向量如下：b1=0.8483,1,0.6833,0.4,0.35,0.3,0.2333,0.1667,0;b2=0.0037,0,0.0296,0.037,

37、0.037,0.0741,0.2593,0.6667,1;b3=1,0.7519,0.4323,0.6105,0.4699,0.3008,0.188,0.1128,0;b4=0.5756,0.2667,0,0.4889,0.4444,0.5778,0.6444,0.6889,1;b5=0,0.0627,0.1749,0.2727,0.3796,0.513,0.6396,0.7927,1;b6=0,0.2899,0.2003,0.288,0.4953,0.6027,0.4001,0.1422,1;b7=0,0.0278,0.3444,0.5931,0.6347,0.7792,0.7917,0.8

38、611,1;b8=1,0.7204,0.6344,0.5086,0.4409,0.4194,0.172,0.3763,0.1075;b9=1,0.9323,0.858,0.8309,0.7383,0.7753,0.2755,0.0796,0;再由有matlab程序：m=polyfit(x1,y,2);y1=m(1).*x1.2+m(2).*x1+m(3);holdplot(xi,y,c*,xi,yi,c);得到9个因素与因变量y缺水量的一元线性关系，得到散点图以及一元线性拟合图。如下图所示： 注：x1-x9颜色分别为r-*红色-农业用水，y-*黄色-生态用水，g-*绿色-工业用水，b-*蓝色-

39、第三产业用水，c-*青色-浅蓝色-人口数量，m-*品红-降水量，k-*黑色-园林绿化覆盖率，m-o红色-废水治理设施数，b-o蓝色-工业排放废水总量。*Matlab程序在附录3.根据图形所示，得出如下结论1. *红色上升曲线（农业用水）：农业用水增多，因变量y缺水量也随之成上升趋势，所以农业用水为其较主要的缺水主导因素。2. *黄色曲线（生态用水）：随着生态用水增加，缺水量减小，说明早年生态用水不为起主导因素，但从07年往后，生态用水逐渐成为其缺水的主导因素之一 。3. *绿色曲线（工业用水）：随着工业用水逐渐增多，缺水量急剧增加，说明工业用水量逐渐成为北京缺水的主要因素。4. *蓝色曲线（第

40、三产业用水）：随着第三产业用水增多，缺水量不断减少，所以第三产业用水不是北京缺水的主导因素。5. *青色曲线（人口数量）：在07年以前，人口数量对y的影响不大，但在07年以后，人口数量越多，缺水量越多，与缺水量y成正比，逐渐成为主导因素。6. *品红曲线（降水量）：降水量少的时候，缺水量y较大。降水量多的时候，缺水量y较小。X6-降水量与缺水量成反比。所以降水量为其主导因素。7. *黑色曲线（园林覆盖绿化率）：随着园林覆盖绿化率逐渐增大，缺水量成下降趋势，所以，园林覆盖率不是北京缺水的主导因素。8. O品红曲线（废水治理设施数）：随着治理套数增多，缺水量增多，说明废水治理设施数的多少不为北京缺

41、水的主导原因。9. O蓝色曲线（工业排放废水总量）:随着工业排放量增多，缺水量也增多，说明工业、废水排放量是北京缺水的主导原因。由以上分析：影响北京缺水的主导因素分别为：农业用水，生态用水，工业用水，人口数量，降水量，工业排放废水总量。根据表1数据，从2000年-2009年，农业用水、工业用水、工业排放废水总量均成下降趋势，但是下降幅度不大，降水量逐渐增多，增多幅度较大，对于缺水有所缓解。但是，生态用水和人口数量逐渐增大幅度较大，成为北京缺水最主要的原因。所以，因为北京市作为中国首都，既是经济中心也是政治中心，生态环境是必须要大面积发展的方面生态用水必然继续成大幅度上升趋势，另外农民工进京，较

42、多人才“北漂”现象，使外来人口增加，必然会导致北京市缺水持续高风险状态。3.3.2 灰色系统GM预测模型（一）灰色系统模型建立方法灰色预测是指利用GM模型对系统行为特征的发展变化规律进行估计预测，同时也可以对行为特征的异常情况发生的时刻进行估计计算，以及对在特定时区内发生事件的未来时间分布情况做出研究等等。这些工作实质上是将“随机过程”当作“灰色过程”， “随机变量”当作“灰变量”，并主要以灰色系统理论中的 GM(1,1)模型来进行处理。灰色预测的方法设已知参考数据列，做1次累加（AGO）生成数列：其中， （k=1，2,.,n）,求均值数列： 则。于是建立灰微分方程为相应的白化微分方程为 记，则由最小二乘法，求使得达到最小值。于是求解白化微分方程，得：（二）灰色系统预测模型 在本题中，我们主要根据北京市在2000年至2009年水资源总量与总用水量之间的供求差值来建立灰色系统预测模型，北京市水资源供求差值具体数字见下表：年份供求差值（亿立方米）年份供求差值（亿立方