水资源短缺风险评估模型建立与分析数学建模

1、摘要一、问题提出二、问题分析三、模型的假设四、模型的建立4.1 模型一4.1.1 1符号说明4.1.2 模型一的建立4.1.3 模型一的求解4. 2模型二4.1.1 符号说明4.1.2 模型二的建立4.1.3 模型二的求解五、模型的优缺点分析六、城市水资源保护利用的建议报告七、参考文献八、感谢信九、附录附录一：涉及水资源短缺风险因子的各项数据表附录二：求解逆矩阵的C+源代码附录三：马氏距离 MATLAB单程序水资源短缺风险评估模型建立与分析摘要目前，水资源问题，特别是水资源短缺问题严重地制约着我国经济发展， 是我们迫切所要解决的问题之一。水资源是人类生产生活最重要的资源，可是 如今，生态环境遭

2、到严重破坏，水体污染严重，水资源开采过度，水资源的保 护和水污染的治理成为现代社会最关注的问题。尤其近年来，我国北方地区水 资源短缺问题日趋严重，已成为人们关注的焦点。本文以北京地区水资源短缺问题为背景，通过建立线性加权均值风险评估 模型和短缺风险模糊综合评估模型，对其水资源短缺情况进行综合风险评估。 同时根据研究分析所得数据，给出了合理使用和保护水资源的建议报告。对于模型一：首先，我们借助软件 SPSS寸所搜集的数据进行筛选处理，按 照影响程度进行加工整合，同时依据评指标选择原则，选取了农业用水量、生 活用水量、工业用水量、地下水资源量、污水排放总量、地表水资源量、入境 水流量、再生水资源量

3、、环境用水量、COD排放总量作为水资源短缺风险的评估因子。其次，随机抽出20组数据分成两大组，作为矩阵 A和B,通过矩阵运 算、加权均值、求解方程组，确定了各项因子的评分系数。在此过程中，我们 通过编写C+取序求解可逆矩阵和借助 MATLAB+算方程组的解，从而确定评估 指标函数，在此基础上确定了水资源短缺评估的临界值，并依次对水资源短缺 风险预警级别进行等级划分。最后，随机选择了2009和2010年的数据进行检验，得出评估模型所求数据与实际情况基本一致。对于模型二：首先，对影响北京水资源短缺的风险因子进行了分析，得到 了一部分重要的风险因子。其次，利用了判别分析法中的Mahalanobis距

4、离法对主要的敏感因子进行筛选，取得了良好的结果，基本符合客观事实。针对水 资源短缺风险评估中各指标的模糊性和不确定性，建立短缺风险模糊综合评估 模型。通过对风险率、脆弱性、可恢复性、事故周期、风险度这五个评价指标 建立综合评估体系，作出了等级划分。通过对风险因子的调控，降低了北京水 资源短缺的风险，并对未来北京市的水资源情况进行预测，对实际生活、生产 具有指导意义。关键字：水资源短缺风险因子线性加权均值风险评估模型短缺风险模糊综合评估模型风险预测建议一、问题提出近年来，我国北方地区水资源短缺问题日趋严重，制约着经济的发展，水 资源短缺已成为人们关注的焦点话题。以北京市为例，北京是世界上水资源严

5、重缺乏的大都市之一，其人均水资 源占有量不足300小，为全国人均的1/8,世界人均的1/30,属重度缺水地区， 附表中所列的数据给出了 1979年至2000年北京市水资源短缺的状况。北京市水 资源短缺已经成为影响和制约首都社会和经济发展的主要因素。政府采取了一 系列措施，如南水北调工程建设，建立污水处理厂，产业结构调整等。但是，气 候变化和经济社会不断发展，水资源短缺风险始终存在。如何对水资源风险的 主要因子进行识别，对风险造成的危害等级进行划分，对不同风险因子采取相 应的有效措施规避风险或减少其造成的危害，这对社会经济的稳定、可持续发 展战略的实施具有重要的意义。根据北京统计年鉴及市政统计资

6、料提供的北京市水资源有关信息，分析 讨论以下问题：(1)以北京市水资源资料为例，分析水资源短缺的风险因子，并对这些风险 因子进行重要性分析；(2)建立一个数学模型对北京市水资源短缺风险进行综合评价，作出风险等 级划分并陈述理由；(3)从用水量、用水结构、水资源存量几个方面对北京市未来五年水资源进 行预测；(4)以北京市有关部门为报告对象，写一份建议报告，至少从水资源短缺成 因、水资源风险控制以及水资源保护几个方面提出建议和对策。二、问题分析2.1 水资源短缺风险评价指标水资源的短缺取决于供水和需水两方面影响，而这两方面都具有随机性和不 确定性。因此，水资源短缺风险也具有随机性和不确定性。在进行

7、风险评价时， 要充分考虑风险的特点以及水资源系统的复杂性，要把存在风险的概率、风险出现的时间、风险造成的损失有多少、风险解除的时间、缺水量的分布等一系列因 素考虑在内。因此难以用某一种指标对其进行全面描述和评价，必须从多方面的指标综合考虑。评价指标选择的原则是：(1)能集中反映缺水地区的缺水风险；(2)能集中反映缺水风险的程度；(3)能反映水资源短缺风险发生后水资源系统的承受能力；(4)代表性好，针对性强，易于量化。依据上述原则，并参考相关文献，选取了影响水资源短缺的主要风险因子， 联系实际，建立了水资源短缺风险评估指标，构造出评估模型。2.2 北京市水资源短缺主要原因(1)我国水资源时空分布

8、不均以及入境水流量少是导致北京水资源紧张的主 要原因；(2)北京人口增长过快，工农业发展迅速，加剧了水资源紧张程度；(3)北京地区生态环境恶化，干旱频率加大；(4)居民生活用水的污染和地表水资源急剧减少；三、模型的假设(1)假设题中所给数据基本真实有效；(2)假设没有重大的自然灾害发生如干旱、洪涝等；(3)假设北京地区人口流动正常；(4)假设南水北调及其它工程正常运行；(5)只考虑模型中的十个主要因子，其它的因子忽略不计。四、模型的建立4.1 模型一4.1.1 符号说明Y：样本数据的评分值C :对应影响因子的评分系数X :各影响因子实际数据Y 0、Y1 ：样本数据评分值的算数平均值XX1: A

9、组、B组原始数据Xi0、X1： A组、B组数据各列的平均值A= X： -X：B= X； - Xj1S=(A+B)Y A0、Y1 ： a组、b组评分值的平均水平Y C :橙色预警的临界值4.1.2模型一的建立水资源短缺风险评估模型是评价水资源短缺程度的评价指标，通过对北京 市水资源资料的调查分析建立这一模型。本模型主要通过对水资源短缺风险因COD排放总量进行了Xo为农业实际用水量。Xi为生活实际用水量。X2为工业实际用水量。X3为地下实际水资源量 X4为污水实际排放总量, X5为地下实际水资源量子：农业用水量、生活用水量、工业用水量、地下水资源量、污水排放总量、 地表水资源量、入境水流量、人均用

10、水量、环境用水量、 研究，它们分别是：1 .农业用水量：C0为农业用水量风险评估系数,2 .生活用水量：Ci为生活用水量风险评估系数,3 .工业用水量：C2为工业用水量风险评估系数,4 .地下水资源量：C3为地下水资源量评估系数,5 .污水排放总量：C4为污水排放总量评估系数,6 .地表水资源量：C5为地下水资源量评估系数,7 .入境水流量：C6为入境水流量评估系数，X6为实际入境水流量。8 .再生水资源量：C7为再生水资源量评估系数，X7为实际再生水资源量9 .环境用水量：C8为环境用水量评估系数，X8为环境实际用水量10 . COD排放总量：C9为COD排放总量评估系数，X9为COD实际排

11、放总以上十项的线性组合Y=C0X0+ClXl+C2X2+C3X3+C4X4+C5X5+C6X6+C7X7+C8X8+C9X9 记为水资源短缺风险评估指标值。（各项水资源短缺风险评估因子具体数据见附 录一）0下面根据收集的水资源短缺风险因子样本资料，求出各评估系数C o从北京市1979年至2001年水资源使用情况数据库中分别随机抽取 s和t个样本分别 记为A、B组。则它们对应的评估指标值为： 组A的评估指标值 组B的评估指标值1 s1 t. -n又做Y =-£ Y , Y =-Z Yi ,即Y为组A的平均值，Y为组B的平均S i 1t i =1值。为使组A与组B之间有明显区别，希望它们

12、平均值之间差距越大越好，而 组内离差平方和越小越好，即L Co,Ci,CzC3c 4c 6c £ ,C,C 9=-(YY1)2(1)£ (Y?-0)2+Z (Y-Y )越大越好,从而建立评估系数 Co,Ci,C2,C3, C4,C5,C6,C7, C8,C9为L(C。温56C9,C ,C ,C ,C ,C ,C ,C )的极大值点。12由微分方程可知（Co，G，C2，C3，C4，C5，C6，C7，C8，C9 ）为方程组:L C0,Ci,C2,C3,C4,C5,C6,C7,C8,C9八二0Ci(i = 0,1，,9)的解4.1.3模型一的求解1.将原始数据写成矩阵组A的矩阵

13、组B的矩阵 各列的平均X0X；) =1-3.2 -4.3 -3.6 -3.7 -3.5-3.1 -2.3 -1.5 -2.3-4.0(XXX»-；)=3.8 -3.1 -2.7 -3.9 -3.6-3.2 -2.9 -2.0 -3.5-2.32做新的矩阵A、B,及两组的算术平均值矩阵S000010-X0-X00Y00Y0X01X1X08X80V 00V 0X11X1X18X8IIIIIIHIIIIIIIX09-X； IY 0 Y0X19X9X90-X：0 y 0X91 一人1 )X98 - X8X99 - X；1-2 0.9-8.8-11.11 1 8 . 9-5 . 1-1 15

14、. 0I 4.82 . 9：3.22 7.56.5 -9.5 -19.04.8 9.65.66.23 . 54.32 89.48.618.44.9 9 . 74.711.73.03 . 52 5:19 .11 0:22 0.24.711.94.95 . 73 . 18.97- 722 .3.73.41 . 15 . 16 . 12 .4-18.7 - 7 . 86.2-11.6- 1A 1 . 3 -14.9-2 1.6-6.1 9 -3 79. 14.16.71.5 -18.51 8 -2 . 28 . 522 .10.29 . 93.97 . 30.49.1333.02.255满足式（1）

15、的（C0,C1,C2,C3,C4,C5,C6,C7,C8,C9）为方程组.51 . 717.49 . 14.18 . 80 . 734842.20.310.12 1 .1 1 .6 . 14.60解得:-C01(元- X011C1X10 一 X11C2X° - X2C3X" X31C4x： -x4S X的解,C5X； -x；C6x： 一 X61C7X - X71C8X； -X；C9 -1X90 _X1-C°l1X00-一 n-X01-0.6 1C1X10-X；-1.2C2X2°-1-X2-0.9C 3X30 -1-X30.2C4X4°-1-x4

16、_ 二0.1C 5=SX50 -X5=S0.1C6X60 -X60.6C7X70 -X70.5C8X8°-X80.8Q 一X90-X9 一i7.7 一Co，Cl，C2，C3，C4，C5，C6，C7，C8，C9 =（其中S通过C+设计程序求得，具体程序见附录二）3.确定评估指标函数Y =6.4X0 5.2X 1 3.8X 2 4.5X 3 2.3X4 1.4X51.6X60.87X72.12X8 0.03X9(2)4.组A、B评估指标值的平均水平YA nCoX； CiXi0 c飘；C9X0 = -95.541Yb1 =CoX0 +GX11+ +CIX8 +C9X9 = -94.879再

17、取出两组数据（每组20个），重复以上过程，最终求得临界值：Yc=-190.425.水资源短缺风险评估等级梯度值划分为:黄色预警I级：得分-147.44,虽然水资源处于短缺状况，但尚且可满 足工农业及生活用水，地下水资源开采较大，地表水、入境水用量较大，污 水处理率较高（再生水用效较多），水资源短缺风险较小。黄色预警II级：147.44得分 -168.93,虽然水资源处于短缺状况， 但尚且可满足工农业及生活用水，地下水资源开采较大，地表水、入境水用 量较大，污水处理率不高，水资源短缺风险也较小。黄色预警III级：-168.93得分 -190.42,水资源处于短缺状况，不能 完全满足工农业及生活用

18、水（有些产业、地区的用水被限制），地下水资源开 采较大，地表水、入境水用量较大，污水处理率不高，水资源短缺风险较 大。橙色预警I级：-190.42得分方-216.91,水资源处于短缺状况，不能 完全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开 采大（有些地区已无地下水可开采），地表水、入境水用量较大，污水处理率 不高，水资源短缺风险也较大。橙色预警II级：-216.91 得分方-233.40,水资源处于短缺状况，不能 完全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开 采大（有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出 现河流干涸、湖水面

19、积不断缩小），污水处理率不高，水资源短缺风险大。红色预警I级：-233.40得分方-259.89,水资源处于短缺状况，不能 完全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开 采大（有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出 现河流干涸、湖水面积不断缩小），污水处理率低，水资源短缺风险大。红色预警II级：-259.89得分 -286.38,水资源处于短缺状况，不能 完全满足工农业及生活用水（有些产业、地区的用水被限制），地下水资源开 采大（有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出 现河流干涸、湖水面积不断缩小，甚至出现旱灾，严重影响了

20、生产、生活） ， 污水处理率低，水资源短缺风险很大。红色预警III级：得分 -286.38,水资源处于短缺状况，不能完全满足 工农业及生活用水（有些产业、地区的用水被强力限制），地下水资源开采大 （有些地区已无地下水可开采），地表水、入境水用量较大（有些地区出现河 流干涸、湖水面积不断缩小，甚至出现旱灾，严重影响了生产、生活、影响 经济发展)，污水处理率低，水资源短缺风险很大。4.2模型二4.2.1 符号说明U评价对象的因素论域V评语论域W 各因素对水资源短缺风险指标的权重R模糊关系矩阵Xg表示横轴自变量(g =1,2,3,4,5,6 )Vg表示纵轴因变量(g =1,2,3,4,5,6 )Cj

21、表示编号为j的第i个因素的值(i =1,2,3,4,5; j =1,2,3|凡30 )rj表示编号为j的第i个因素的值在30年的同一因素值的总和中所占的比 例A各风险因子的的权重分配1表示因素论域U中第i个因素u对应于评语论域 V中第j个等级Vj的 相对隶属度W表示各风险因子对水资源短缺风险指标的权重。4.2.2 模型二的建立(1)水资源短缺的影响因子大概可分为两类：自然因素，即北京市总人口，水资源总量，水位深度，河 流和气候等；社会经济因素，污水排放和处理，农业用水，工业用水，服务业 用水等由于经济发展对水资源需求。(2)判别分析判别分析可用于识别影响水资源短缺风险的敏感因子能够从诸多表明观

22、测 对象特征的自变量中筛选出提供较多信息的变量且使这些变量之间的相关程度 较低线性判别函数的一般形式如下其中为y判别分数，乙,4,为反映研究对象特征的变量，ai, a21an为 各变量的系数也称判别系数。常用的判别分析法有欧式距离法和马氏距离法，由于欧式距离判定法存在这一定的缺点，故常用的判别分析方法是 Mahalanobis距离法，即每步都使得 相距最近的两类间的Mahalanobis距离最大的变量进入判别函数，其计算公式如 下：其中x是某一类中的观测量，Y是另一类，可以求出x与Y的Mahalanobis距离。(3)水资源短缺影响因子筛选 根据上文中提出的水资源短缺风险影响因子， 利用Mah

23、alanobis距离法筛选出水资源短缺风险敏感因子,见表5。从表5中第3 栏可以看出，水资源总量、工业排放总量、农业用水量、生活用水量在步骤1至步骤4中移出模型的概率均小于 0.1,同时在每步中这4个变量均使得最近的两 类间的Mahalanobis距离最大因此，这4个变量是影响北京地区水资源短缺风险 的敏感因子且影响力由大到小为：水资源总量 农业用水 第三产业用水 工业用水步 骤影响因素容许度移出概率最小马氏 距离的平方组问1工业用总量1.0000.0892工业用总量0.6820.0200.1842,5水资源总量0.6820.0000.2361,4工业用总量0.3910.0280.8461,5

24、3水资源总量0.6780.0000.7222,4农业用水量0.4600.0341.2252,5工业用总量0.2510.0356.5501,54水资源总量0.3280.0001.3852,4农业用水量0.1230.0031.2422,5第二产业用水量0.1020.0232.9632,5敏感因子筛选(4)模糊的分析：由于降水量，河流的流量的随机性，北京市供水和需水量 也存在这不确定因素。因此水资源短缺也存在一定的随机性，具有水资源短缺 的风险。所谓水资源短缺风险是指在特定的环境条件下，由于供水水和用水两 方面存在不确定性，使区域水资源系统发生供水短缺的概率以及由此产生的损 失。水资源短缺风险评价指

25、标主要描述水资源短缺的状况，只针对水资源的系 统不包括社会经济系统，是水资源短缺风险却化的基础。选取风险率、脆弱 性、可恢复性、重现期和风险度几个评价指标利用模糊综合评价方法对水资源 短缺风险进行综合评价。(5)评价指标1 .风险率根据风险理论，载荷是造成系统非正常状态的动力，抗力是维护系统正常 的能力。如果把水资源系统的失事状态记为FJRaP),正常状态记为S-（R>P）,那么水资源系统的风险率为其中，Xi为水资源系统状态变量。如果水资源系统的工作状态有长期记录风险 率也可以定义为系统不能正常工作的时间和整个系统工作时间的比值，即： 其中，NS为水资源系统工作的总时间；It为水资源的状

26、态变量。2 .脆弱性脆弱性是描述水资源系统失事损失严重程度的期望值的大小。假定系统第1次失事的损失程度为S ,它相对应的发生概率为 P，那么水资源系统的脆弱 性为：其中，nf为系统失事状态的总次数。3 .重现性周期是两次进入非常状态模式 F之间的时间问隔，也叫平均重现期。用 d（,n读示第n间隔时间的历时，则平均重现期为： 其中，N =N（N）为10,t】时段内属于模式F的事故数目。4 .可恢复性可恢复性是指系统从失事状态恢复到正常状态的可能性。恢复性越高说明 系统从失事状态转入正常状态的时间越短。可以用下面条件概率来定义： 引入整数变量： 这样可得到：5 .风险度用概率分布的数学特征，如标准

27、差吟可以说明风险的大小。.越大,则风险 越大，反之越小。这是因为概率分布越分散，实际结果远离期望值的概率就越大。 4.2.3模型二的求解采用上述定义的的风险率、脆弱性、可恢复性、重复性、风险度作为评价 指标，采用模糊综合评价方法对风险进行评价。设定两个有限论域（即集合）：U = u1,u2,u3,lll,umKDV =»*,111&其中，U表示综合评判的因素所 组成的集合，V表示评价语组成的集合。则模糊综合评判即表示下列的模糊变 换8=人口,式中A为U上的模糊子集。而评判结果 B是V上的模糊子集并且 可以表示人=（獴必|,£ ）, 0 < Z «1；

28、 B=（h，b2，|但% ）,0 <bj <1；。其中九i表示单 因素Ui在总评定因素中所起作用大小的变量，也在一定程度上代表根据单因素 Ui评定等级的能力；bj为等级Vj对综合评定所得模糊子集 B的隶属度，它表示 综合评价结果。关系矩阵R可表示为：式中：。表示因素u的评价对等级w的隶属度，因而矩阵R中R =（'2,111,几） 即为对第i个因素u的单因素评判结果。在评价计算中 A =（九，儿2,1也鼠）代表了各 个因素对综合评判重要性的权系数，因此满足 £ % =1,（i =1,2,|,m）;同时，模糊变换AmR也即退化为普通矩阵计算，即bj =min（1,&#

29、163;加rj ）,i = 1,2| ,m;j =1,2,|,n。上述权系数的确定可用层次分析法 AHP马到。由上述分析可以得到：评价因素集 U对应评语集V,而评判矩阵中rj即为某因素ui对应等级Vj的隶属度，其值课根据各评价因素的实际数值对照各因素的分级指标推求。我们把评语级划为 5个如下等级：水资源短u缺风险 风险脆弱性u2可恢复u3重现性u4风险度u5对于水资源系统的风险率，脆弱性，风险度是越小越好，所以按越小越 优性指标对3,2刀5构建如下的隶属函数：对于水资源系统可恢复性和重现性则是越大越优性指标，对U3,山各评价语级构建如下隶属函数：对于水资源短缺风险评价的因素集U而言，对应一个测

30、定指标向量Y =（*11、24,。4丸）。其中七是U寸于uj的测定值。这样uVi（*j ）就是相对于ui属于Vi的程度。对于因素集U®有一下的模糊关系矩阵；水资源短缺风险评价各因素的权重确定采用层次分析法（AHP）,设权重计算结果为A = （4%,%54,%）,于是可得出综合评判向量：我们选取“加权平均型”的模型进行综合评价，即bj =min（1Z %irj ）nz %rj <1,所以这个模型可以用一般实数加法，即：i 1maxbj 对应的评价语即为水资源短缺评价的评价结果对风险程度的综合评价可以将其划为5级，制作如下的级别评价表格分别 为：低风险，较低风险，中风险，较高风险，

31、高风险。水资源系统风险等级划分风险等级风险级别风险特征1低风险可以忽略的风险2较低风险可以接受的风险3中风险边缘风险4较图风险比较严重的风险5局风险无法承受的风险五、模型的优缺点分析模型一的优缺点：本模型的优点是操作简单，实用性强。通过大量随机数据的验证，模型基本符合实际情况，对水资源短缺风险评估具有一定的帮助作用。另外，在模型 的计算中进行了 C+刑序和MATLABt机实现，对以后水资源短缺风险评估标准 的程序化、软件化具有一定的指导作用。本模型的缺点是考虑的情况较为理想化，有地区局限性。模型二的优缺点：一般风险研究只对个别风险指标进行描述，而对水资源短缺的综合评价则 比较少，而模糊综合评价

32、模型从不同角度对水资源系统进行了比较全面的风 析，与实际情况比较符合。在分析敏感因子简化问题，抓住了主要原因，取得 了良好的结果。本模型的缺点是模型的使用范围存在着一定的局限性。六、城市水资源保护利用的建议报告目前，北京市水资源短缺问题日益严重，已沦为世界上水资源严重匮乏的大 都市之一。水资源短缺问题影响着人们生产生活的正常进行， 制约着北京市经济 的可持续发展。通过对水资源状况的综合分析，我们建立了线性加权均值风险评 估模型和短缺风险模糊综合评价模型，得出了影响北京市水资源短缺的主要风险 因子，划分了水资源短缺风险预警级别，并且对北京市未来水资源短缺情况进行预测。对此，我们提出以下几点浅陋建

33、议，希望对水资源的保护利用有所帮助(1) 强水土保持与综合治理日益严重的水土流失使地面土壤表层不断受到破坏， 造成土壤养分流失，制 约了经济发展，导致农民生活贫困，更为重要的是破坏了流域的生态平衡。因此, 治理水土流失，开发当地水土资源是改善当地生态环境，促进国民经济可持续发展的主要措施， 可通过以下措施加强水土保持及综合治理工作。(1) 通过生物措施工程措施和耕作措施的合理配置，使流域形成由水、林、田、路统一规划，坡面、沟道均有水保工程的立体防护体系。(2) 根据土地资源适宜性评价结果，确定流域经济与生产发展方向将陡坡地段的坡耕地退耕还林，增加林地用地比例，大举植树造林 或封禁治理，提高土地

34、利用率，减少水土流失面积，逐渐使流域的土地利用结构 趋于合理，农林牧各业协调发展，改善生态环境。(2) 强水资源环境保护与污染治理切实加强水资源的环境保护，保护现有植被，减少水土流失及地面径流对工 业园区设置单独的污水处理厂，严禁污水直接排入自然水体将合理利用水资源与 控制水污染相结合，提高水资源重复利用率，近、中、远期分别达到60% 70%、 80%Z上建筑工地废水必须经过处理回用或合格排放，禁止无组织占用道路经营 洗车。通过固体废物控制，优化污染治理，进而推进水资源保护工作。(1)生活垃圾的处理采用定点收集垃圾，集中无害处理，焚烧和填埋并行， 并发展综合利用家庭生活垃圾收集：实施家庭垃圾分

35、类袋装化，根据垃圾的可否 再生利用处理难易程度等特点，由家庭事先进行分类袋装，在居住区设置专用垃 圾收集间和特定集装箱，在商业中心设置垃圾收集站，对商业垃圾进行分类收集， 并纳放家庭垃圾收集系统。(3) 工业垃圾的处理以填埋方式为主。(4) 医疗废弃物必须单独收集运输并用焚烧法处理，避免与一般垃圾相混。(5) 建筑垃圾以填埋方式为主，对施工产生的余泥尽可能就地回填或申报有 关管理部门及时运走，堆放到合适地方，并注意清洁运输，防止运输过程中洒落， 对于建筑固体废物除部分木材和竹料经再加工可再利用外，一般不能重新利用， 需要堆置存放。3 .实行严格的水资源管理制度结合现有水资源管理的经验，首先在水

36、资源管理观念上应实现六个转变：一是在管理理念上，要加快从供水管理向需水管理转变；二是在规划思路上，要把 水资源开发利用优先转变为节约保护优先；三是在保护举措上，要加快从事后治 理向事前预防转变；四是在开发方式上，要加快从过度开发无序开发向合理开发 有序开发转变；五是在用水模式上，要加快从粗放利用向高效利用转变；六是在 管理手段上，要加快从注重行政管理向综合管理转变。其次，实行严格的水资 源管理制度，严格水资源开发利用红线、水功能区限制纳污红线和用水效率控制 红线，通过完善的水资源法规体系，强有力的水资源监管措施，坚决打击违法填湖行为，强化取水资源论证、节水考核、入河排污口设置、水域岸线利用、河

37、 道采砂、水工程建设等方面的管理和执法监督，严格执行水资源管理的各项法规。4 . 加强用水节能教育，推进全面节水工程(1)推动节水型社会建设。认真研究提出关于开展节水型社会建设的指导性文件，适时召开全国节水型社会建设工作会议， 继续开展全国节水型社会建设试 点工作，研究提出水资源宏观分配指标和微观取水定额指标， 推进国家水权制度 建设。(2) 开展城市节水工作。积极开展节水产品研发，加大节水设备和器具的推 广，指导各地加快供水管网改造，降低管网漏失率，推动公共建筑、生活小区、 住宅节水和中水回用设施建设，推进污水处理及再生利用，加快城市供水和污水 处理市场的改革，大力开展雨水回收利用，建筑物安

38、装屋顶虹吸式雨水收集系统， 通过这个收集系统，尽量收集更多的雨水，采用植草砖、渗透砖收集地表径流水， 经处理后做景观绿化用水，并在适当的地方建设地窖存放丰雨期收集到的雨水。(3) 推进农业节水。继续推进农业节水灌溉，推广农业节水灌溉设备应用， 大力推进大中型灌区节水改造，积极开展农业末级渠系节水改造试点。在丘陵、 山区和干旱地区积极开展雨水积蓄利用，支持农村水窖建设，推广旱作农业技术， 发展旱作节水农业，扩大节水作物品种和种植面积，开展农村、集镇生态卫生旱 厕试点。(4) 推进节水技术改造和海水利用。推进高耗水行业节水技术改造，矿井水 资源化利用，推进沿海缺水城市海水淡化和海水直接利用。总之，

39、合理利用并保护水资源是北京市自然环境与自然资源保护中的一项重 大任务，只有妥善解决，才能实现北京市社会稳定，经济的平稳发展。以上所述几点建议，只是我们个人浅显的看法和认识。由于所学知识和阅历 有限，言辞中难免有所不当，望各位领导批评指正！七、参考文献1阮本清,梁瑞驹,陈韶君.一种供水系统的风险分析与评价方法J.水利学报,2000,(9):17.2傅湘,王丽萍,纪昌明.极值统计学在洪灾风险评价中的应用J.水利学报,2001,(7):610.3阮本清，韩宇平，王浩等，水资源短缺风 险的模糊综 合评价J.水利学 报,2005(8).4 Iman Karimi,Eyke Hllermeier,Risk

40、 assessment system ofnatural hazards : A new approach based on fuzzy probability, Fuzzy Sets and Systems.2007,987-999.5杨树滩，更自强.模糊数学在水资源紧缺程度评价中的应用研究J.长江科学院院报,2005,22(1).6林武星，吴.基于自适应向量评估遗传算法的水资源优化配置模型及应用J.南水北调与水利科技,2008,6(3):69-71.7左其亭，吴泽宁，赵伟.水资源系统中的不确定性及风险分析方法 J.干旱区 地理.2003,(2).8胡国华，夏军.风险分析的灰色随机风险率方法

41、研究J.水利学 报,2001,(4):15.9陶涛，付湘，纪昌明.区域水资源供需风险分析的应用研究J.武汉大学学报(工学版),2002,(3).10黄明聪，解建仓，阮本清等.基于支持向量机的水资源短缺风险评价模型及应用J.水利学报，2007,38(3).11冯平.供水系统干旱期的水资源风险管理J.自然资源学报，1998,13(2).12八、感谢信通过此次建模活动，我们不仅强烈的感受到水资源短缺问题的严重性及解 决的必要性，也充分的意识到了数学知识在现实生活中解决问题的作用之大， 团队合作的力量之大。在此，首先，感谢学校能提供此次建模活动，这不仅让我们感受到了数学 知识在解决实际问题时的威力，还

42、让我们的专业知识得到了进一步巩固和提 升，提高了团队合作的意识。其次，我们真心感谢吴红萍老师、杨小东老师对我们的认真指导，无论从数据的搜集还是模型的建立以及论文的成稿，无不体现了他们治学与作研究的 严谨态度，同时被他们的精神所感染。他们是我们学习的榜样，真心感谢老师 能抽出宝贵的时间来指导我们。在此，衷心祝老师工作顺利、身体健康！九、附录附录一：涉及水资源短缺风险因子的各项数据表表1: 1979年至2000年北京市水资源短缺的状况年份总用水 量（亿立 方米）农业用 水（亿立 方米）工业用 水（亿立 方米）第二产业及生活等具它用水（亿立方米）水资源总量（亿方）197942.9224.1814.3

43、74.3738.23198050.5431.8313.774.9426198148.1131.612.214.324198247.2228.8113.894.5236.6198347.5631.611.244.7234.7198440.0521.8414.3764.01739.31198531.7110.1217.24.3938198636.5519.469.917.1827.03198730.959.6814.017.2638.66198842.4321.9914.046.439.18198944.6424.4213.776.4521.55199041.1221.7412.347.0435.

44、86199142.0322.711.97.4342.29199246.4319.9415.5110.9822.44199345.2220.3515.289.5919.67199445.8720.9314.5710.3745.42199544.8819.3313.7811.7730.34199640.0118.9511.769.345.87199740.3218.1211.111.122.25199840.4317.3910.8412.237.7199941.7118.4510.5612.714.22200040.416.4910.5213.3916.86年份降水量 （毫米）气（心均温二）日照时

45、数（时）平均风 速（米/ 秒）平均气 压（百 帕）大风 日数 （日）雨日 数（日）1978664.811.62865.42.61012.835641979718.411.12667.42.51012.233631980380.711.02920.82.51012.729831981393.212.32803.92.51010.815921982544.412.82825.12.61010.526921983489.913.02844.32.41010.3291001984488.811.92767.62.41010.618901985721.011.52511.92.21010.4121041

46、986665.312.12804.12.31010.721961987683.912.32631.92.41010.3231021988673.312.72558.12.41010.817961989442.213.22626.21.91011312.72325.01.91010.6121131991747.912.52536.62.11010.88981992541.512.82712.52.21011.061001993506.713.02669.82.61010.812911994813.213.72470.52.51010.19921995572.513.32

47、519.12.61010.316891996700.912.72418.72.61011.0161031997430.913.12596.52.51012.911761998731.713.12420.72.31012.510931999266.913.12594.02.41012.57862000371.112.82667.22.51012.710832001338.912.92611.72.41012.910782002370.413.22588.42.31012.715842003444.912.92260.22.51013.36932004483.513.52515.42.41012.

48、612942005410.713.22576.12.41012.85792006318.013.42192.72.21012.55862007483.914.02351.12.21012.65782008626.313.42391.42.21012.68100表2:北京市气象情况（1978-2008年）表3:排水及节水附录表年份降项量 （亳米、平均 笠温日照时 券 （时）平均风 速（米20屋气8± （百大风日2（ 及（日）07日数 （日）台匕一H、匕 / n 秒）帕）no1978处理能 664.8（1帝力弘65)42.63291012.83553 641979 级718.4（1方1月

49、干）42.5 319 1012.233 313 631980380.711.02920.82.51012.7298319眼年纪.理 393.21(23立 7122803.92.5104：芍1010.815 98865921982_5444-12 82825 12 61010 526921 DONO*T*T .*TIN.。N0NJ. 1N .U1 U 1 U. U污水厂（万立方米）9843891161二三级（万立方米）9607090302污水处理率(%)78.976.2集中处理率(%)74.569.6污水排放总量（万立方米）132095129820排水管道长度（公里）88818526污水管（公里）44584357再生水利用量（万立方米）6000049501节水量（万立方米）1955916513节水措施（项）245370表4：北京市200