长江大学学报（自然科学版）理工２０１２年１０月第９卷第ｌｏ期 Ｊｏｕｒｎａｌ ｏｆ Ｙａｎｇｔｚｅ Ｕｎｉｖｅｒｓｉｔｙ（Ｎａｔ Ｓｃｉ Ｅｄｉｔ）Ｓｃｉ＆Ｅｎｇ Ｏｃｔ．２０１２，Ｖｏ１．９ Ｎｏ．１０ ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７３—１４０９（Ｎ）．２０１２．１０．００８ 城市水资源质量系统动力学仿真与评价研究 毛敏华 （南水北调中线干线建设管理局机电物资部，北京１ ００３８２） ［摘要］在分析城市水资源质量影响因素的基础上，建立了城市水资源质量系统动力学仿真模型，模型由 工业用水子系统、城市居民用水子系统及水污染处理子系统构成。通过ｗ市实际统计数据对模型调试检验 后，以工业产值增长率、人口增长率和污水处理率为决策变量，以３种发展方案对ｗ市未来２Ｏ年的ＣＯＤ 排放量仿真，仿真结果说明城市发展规划需要协调发展，特别是对人口增长率及水污染处理率不能忽视。 ［关键词］水资源质量；系统动力学；仿真；评价 ［中图分类号］ＴＶ２１ ３．４ ［文献标识码］Ａ ［文章编号］１６７３—１４０９（２０１２）】０一Ｎ０２２—０４ 水是城市经济协调健康发展的战略性基础资源，也是生态环境系统维持良性循环的控制性要素，城 市经济发展和水安全之间的互动影响关系越来越受到重视［１］。城市飞速发展的同时给区域内用水安全带 来了更大的压力，保证城市水资源质量已是一个很重要的课题。 在对城市水资源质量的研究中，大部分主要针对居民生活废水的净化再利用环节进行研究，李梅将 系统动力学应用于污水再生利用的系统分析和预测 ；徐志嫱等用系统动力学方法建立了分散式污水再 生回用系统模型ｌ＿３ ；黄微等在污水再生利用规划中运用系统动力学模型针对石家庄特点设计多情景，并 以非线性优化方法得到具体再生水水源和回用对象的配水方案 ］。以上的研究中将城市居民生活用水造 成的污染作为城市水资源质量单一影响因素进行分析，但实际情况中工业的发展也是城市的重要组成部 分；与此同时，在水资源承载力、需水量预测、水资源供需平衡分析等方面也开展了广泛的研究，但存 在系统结构不完整、反馈回路断裂的问题影响系统模拟效果，且未对水体质量作深入分析。 系统动力学（ｓｙｓｔｅｍ ｄｙｎａｍｉｃｓ，ＳＤ）是一种将结构、功能和历史相结合，以反馈控制理论为基础， 借助计算机进行模拟仿真而定量地研究高阶次、非线性、多重反馈复杂时变系统的系统分析理论和方 法，是研究具有复杂反馈关系系统的重要方法。其本质是带时滞的一阶微分方程组，可作为实际系统， 特别是社会、经济、生态等复杂大系统的“实验室”，是研究水资源系统的重要方法之一［５］。 下面，笔者针对城市水资源质量变化影响因素，应用系统动力学方法对其进行仿真建模分析，通过 已有的经济水文数据检验调试模型的精度后，选取工业产值增长率、人口增长率和污水处理率为决策变 量，通过３种不同的城市发展方案对ｗ市未来２０年水资源质量进程仿真，通过仿真结果的分析对城市 未来发展规划提出相应的政策建议。 １城市水资源质量ＳＤ模型的建立 城市水资源的主要污染源为居民生活污水及工业污水，经过污水处理设施后，这部分水再次进入城 市内部循环，建立系统动力学模型首先需要对建模的目的及模型边界、假设条件等进行确定。笔者建立 的城市水资源质量模型，拟解决以下几个问题：①更加直观的表现城市水资源质量的主要影响因素，特 别是经济发展导致居民生活用水及工业用水的增加，同时经济的发展也会使城市水污染处理能力的增 加，通过建模较为直观的表述上述因素之前的相互作用、相互影响的关系；②通过调节模型决策变量来 制定多种方案，模拟不同方案下水资源质量的变化情况；③定量描述Ｗ市未来２Ｏ年水资源质量的动态 变化情况；④将各方案对比分析，确定最优方案，为决策者提供依据。 由于数学模型是现实情况的抽象变现，故允许存在一定的误差，笔者建立的城市水资源质量系统动力学 模型的边界为对象城市范围内，将城市的水资源相对独立起来进行研究，且假设城市污水处理能力不足。 ［收稿日期］２０１２—０５—２５ ［作者简介］毛敏华（１９６８一），男，１９９０年大学毕业，高级工程师，现主要从事长距离输水的安全及可靠性方面的研究工作。 第９卷第１Ｏ期 毛敏华等：城市水资源质量系统动力学仿真与评价研究 建模的主要目的及模型的边界、相关假设条件确定后，运用系统动力学软件Ｖｅｎｓｉｍ ＰＩ Ｅ得到城市 水资源质量主要因果关系回路，按照城市主要污水产生系统分为居民生活污水及工业产生污水，具体的 城市水资源质量主要因果关系回路如图１所示。图１中共有２个正（＋）反馈回路：工业产值、城市人 口一社会总产值一污水处理投资一污水处理量一ＣＯＤ排放量一水资源质量一工业产值、城市人口，在一定时期 内，随着工业产值及城市人口的增加，社会总产值必将随之增加，于是政府对污水治理投资加大， ｃＯＤ排放下降，从而使水资源质量上升，反过来进一步促进工业产值及城市人口的增加；２个负（一） 反馈回路：工业产值、城市人口一总污水量一直接排放量－ＣＯＤ排放量一水资源质量，主要表现为工业产 值和城市人口使得总体污水排放量增加，尽管污水处理能力也在上升，但总体结果还是导致ＣＯＤ排放 量上升，并最终导致水质恶化＿６］。 ２城市水资源质量ＳＤ分析 城市水资源质量主要因果关 系回路如图１所示，由图１可知 城市水资源污染源主要分为工业 污染和居民生活污染，同时将水 污染处理作为单独的子系统，将 量 工业用水、城市居民用水和水污 染处理分别作为独立的子系统作 详细的系统动力学分析 ］。 ２．１工业用水子系统 工业用水子系统即城市中工 图１城市水资源质量主要因果关系回路 业导致的污水量及相应的污水处 理系统。该方法是通过调查工业万元产值取水量的现状和历史变化趋势，推测目前或将来为实现某一工 业产值目标所需的工业用水量，并通过单位产值工业污水产生率研究相应时间段内工业污水的产生量。该 子系统的核心变量是工业产值，由于工业用水的统计数据相对较少，该研究通过设定表函数，结合城市相 关发展规划指标的方法来表示工业产值增长率，确定工业用水量的变化状况。该系统ＳＤ结构如图２所示。 ２．２城市居民用水子系统 选定人口数量作为流位变量，人口变化率作为流率变量，总人口变化率可以通过多年的统计数据， 通过ＶＥＮＳＩＭ软件的表函数、分段函数相结合的方法来确定。该系统ＳＤ结构图如图３所示。 工业产值 万元产值用水量 量 ／／ 增长 人均用水额度 ． 城市生活需水量 图２工业用水子系统 图３城市居民用水子系统 ２．３水污染处理子系统 工业用水子系统及城市居民 生活用水子系统产生的污水通过 城市污水处理系统处理后排放 （见图４），由于城市内污水处理 能力有限，故有一部分废水直接 排放，处理后污水与直接排放污 水的区别在于ｃＯＤ浓度的不同， 同时随着工业产值与城市人口的 图４水污染处理子系统 长江大学学报（自然科学版）理工＊水文与水资源 ２０１２年１Ｏ月 增加，导致社会总产值的上升，政府加大污水处理力度使得污水处理率呈动态变化趋势。 ３模型的调试与检验 根据建模的目的，模型运行之前应对模型结构和功能进行全面的检验，以保障模型的有效性和真实 性。系统动力学模型的适用性与一致性的检验过程包括以下４组：结构的合适性检验、模型行为的适用 性检验、模型结构与实际系统的一致性检验、模型行为与实际系统的一致性检验 。 将Ｗ市２００６～２０１０年的历史参数输 入模型进行模型仿真，然后与历史实际发 生的行为、数据进行对比。选取人口数量 和工业产值进行历史验证。验证结果如表 表１模型的运行结果检验 １所示。由表１可知，人口数量和工业产 值的仿真结果与历史值基本符合，其中人 口数量的最大误差不超过６ ，工业产值 额最大误差不超过２ ，认为模型模拟结 果与实际值拟合较好，可运用该模型进行仿真。 ４城市水资源ＳＤ仿真与结果分析 已建立的城市水资源质量ＳＤ模型通过实际数据的仿真检验保证了模型的精度符合要求，系统动力 学又称为政策实验室，决策变量的选取引用多指标综合评价法中指标的选取原则，结合ｗ市水资源质 量影响因素的具体情况，选取工业产值增长率、人口增长率、污水处理率３个量作为决策变量。通过改 变决策变量的值得到不同方案下相关指标的模拟结果。再通过对不同方案的模拟结果分析方案可取之处 及需要改进的地方。 选取３个方案进行仿真模拟，其中方案１为现实情况，方案２及方案３为假设情况，具体３种情况 的决策变量值如表２所示。这３种方案变现出了３种不同的城市发展理念，方案２主要强调保持工业产 值的高速增长，同时对人口增长率的控制较方案１更为严格；方案３强调工业产值与人口增长率的稳步 发展，虽然工业产值不如其他方案增加的快，但对人口增长控制最为严格，同时方案３对污水处理的投 入更大，使污水处理率维持在一个高水平。通过ＳＤ模型仿真后，得出城市资源水质量动态变化情况， 即ＣｏＤ排放量仿真结果如表３所示。 表２各方案决策变量取值 １．８５ １．９Ｏ １．９０ １．９Ｏ 以上仿真结果表明，针对ｗ市水资源质量情况， 工业用水量与城市居民用水量产生的污水是影响水质的 主要原因，仿真方案中，虽然方案２的工业增长值最 大，但由于对人口增长控制不严及污水处理率不高， ＣＯＤ排放量均大于方案３，同时该方案中还有很多指． 年份 表３各方案模拟仿真结果 方案１ ＣＯＤ排放量／１０ ｔ 方案２ 方案３ 标有待改进，说明在城市的发展规划中，不能只关注工 业产值的增长，人口的控制及污水处理水平的提升也不容忽视。 第９卷第１ｏ期 毛敏华等：城市水资源质量系统动力学仿真与评价研究 ・２５・ ５结 语 以城市水资源质量为研究对象，用系统动力学方法建模。通过Ｗ市实际统计数据对模型调试与检 验后仿真，将模型仿真精度校准到规定范围后，选取城市发展规划中工业产值增长率、人口增长率及水 污染治理率为决策变量，了解３个参数对Ｃ０Ｄ排放量产生的影响，为城市发展相关政策分析起参考作用。 仿真结果表明，只有工业产值增长率、人口增长率和水污染治理率协调发展才能保障城市水资源质量。 ［参考文献］ ［１］ 钱正英，张光斗．中国可持续发展水资源战略研究综合报告及各专题报告［Ｍ］．北京：中国水利水电出版社，２００１． ［２］ 李梅．城市污水再生回用系统分析及模拟预测［Ｄ］．西安：西安建筑科技大学，２００３． ［３］ 徐志嫱，黄廷林．分散式污水再生回用系统优化分析［Ｊ］．西安建筑科技大学学报，２００５，３７（２）：１６４—１６８． ［４］ 黄徽，杜鹏飞，曾思育．基于系统动力学模型的石家庄污水再生利用规划口］．清华大学学报（自然科学版），２０１０，５０（３）：３９１—３９５． ［５］ 王其藩．系统动力学［Ｍ］．上海：上海财经大学出版有限公司，２００９． Ｅ６］ 王薇，雷学东，余新晓，等．基于ＳＤ模型的水资源承载力计算理论研究——以青海共和盆地水资源承载力研究为例［Ｊ］．水资源与 水工程学报，２００５，１６（３）：１１—１５． ［７３ 惠泱河，蒋晓辉，黄强，等．水资源承载力评价指标体系研究［Ｊ］．水土保持通报，２００１，２１（１）：３０—３４． ［８］ 陈红，回燕斌．辽中南城市群水资源承载力分析［Ｊ］．科技情报开发与经济，２００７，１７（１ｏ）：１５８—１６０． ［编辑］ 洪云飞 （上接第２１页） ３）光催化氧化对渗滤液的ＤＯＭ中的绝大部分物质具有较明显降解效果，但对ＤＯＭ不同组分的 作用效率和规律差异较大，并且光催化氧化对ＨＩＡ组分的降解并不明显。可见ｒ，ＨＩＡ组分比较稳定， 难以光催化降解，是制约光催化氧化处理效率提高的主要组分。 ［参考文献］ ［１］Ｌｕ Ｆ，Ｚｈａｎｇ Ｈ，Ｃｈａｎｇ Ｃ Ｈ，ｅｔ ａ１．Ｄｉｓｓｏｌｖｅｄ ｏｒｇａｎｉｃ ｍａｔｔｅｒ ａｎｄ ｅｓｔｒｏｇｅｎｉｃ ｐｏｔｅｎｔｉａｌ ｏｆ ｌａｎｄｆｉｌｌ ｌｅａｃｈａｔｅ［Ｊ］．Ｃｈｅｍｏｓｐｈｅｒｅ，２００８。 ７２（９）：１３８ｉ－１３８６． Ｅ２］Ｈｕｏ Ｓ Ｉ ，ⅪＢ Ｄ，Ｙｕ Ｈ Ｃ，ｅｔ ａ１．Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｄｉｓｓｏｌｖｅｄ ｏｒｇａｎｉｃ ｍａｔｔｅｒ（ＤＯＭ）ｉｎ ｌｅａｃｈａｔｅ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｌａｎｄｆｉｌｌ ａｇｅｓ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅｓ－Ｃｈｉｎａ，２００８，２０（４）：４９２—４９８． ［３］Ｈａｎ Ｎ，Ｔｈｏｍｐｓｏｎ Ｍ Ｌ．Ｃｏｐｐｅｒ—ｂｉｎｄｉｎｇ ａｂｉｌｉｔｙ ｏｆ ｄｉｓｓｏｌｖｅｄ ｏｒｇａｎｉｃ ｍａｔｔｅｒ ｄｅｒｉｖｅｄ ｆｒｏｍ ａｎａｅｒｏｂｉｃａｌｌｙ ｄｉｇｅｓｔｅｄ ｂｉｏｓｏｌｉｄｓ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｑｕａｌｉｔｙ，１９９９，２８（３）：９３９ ９４４． ［４］ｄｅ Ｍｏｒａｉｓ Ｊ Ｌ，Ｚａｍｏｒａ Ｐ Ｕｓｅ ｏｆ ａｄｖａｎｃｅｄ ｏｘｉｄａｔｉｏｎ ｐｒｏｃｅｓｓｅｓ ｔｏ ｉｍｐｒｏｖｅ ｔｈｅ ｂｉｏｄｅｇｒａｄａｂｉｌｉｔｙ ｏｆ ｍａｔｕｒｅ ｌａｎｄｆｉｌｌ ｌｓｅｃｈａｔｅｓ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｈａｚａｒｄｏｕｓ Ｍａｔｅｒｉａｌｓ，２００５，１２３（３）：１８１—１８６． ［５］Ｃｈｏ Ｓ Ｐ，Ｈｏｎｇ Ｓ Ｃ，Ｈｏｎｇ Ｓ Ｌ Ｓｔｕｄｙ ｏｆ ｔｈｅ ｅｎｄ ｐｏｉｎｔ ｏｆ ｐｈｏｔｏｃａｔｌａｙｔｉｃ ｄｅｇｒａｄａｔｉｏｎ ｏｆ ｌａｎｄｆｉｌｌ ｌｅａｃｈａｔｅ ｃｏｎｔａｉｉｎｎｇ ｒｅｆｒａｃｔｏｒｙ ｍａｔｔｅｒ口］．Ｃｈｅｍｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｊｏｕｍａｌ，２００４，９８（３）：２４５—２５３． ［６］Ｌｏｕ Ｚ Ｙ，Ｚｈａｏ Ｙ Ｃ，Ｙｕａｎ Ｔ，ｅｔ ａ１．Ｎａｔｕｒａｌ ａｔｔｅｎｕａｔｉｏｎ ａｎｄ ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ ｏｆ ｃｏｎｔａｍｉｎａｎｔｓ ｃｏｍｐｏｓｉｔｉｏｎ ｉｎ ｌａｎｄｆｉｌｌ ｌｅａｃｈａｔｅ ｕｎｄｅｒ ｄｉｆｆｅｒｅｎｔ ｄｉｓｐｏｓｉｎｇ ａｇｅｓ口］．Ｓｃｉｅｎｃｅ ｏｆ ｔｈｅ Ｔｏｔａｌ Ｅｎｖｉｒｏｎｍｅｎｔ，２００９，４０７（１０）：３３８５—３３９１． ［７］Ｔｈｅｉｓ Ｔ Ｌ，Ｙｏｕｎｇ Ｔ Ｃ，Ｈｕａｎｇ Ｍ，ｅｔ ａ１．Ｌｅａｃｈａｔｅ Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｎｄ Ｃｏｍｐｏｓｉｔｉｏｎ ｏｆ Ｃｙａｎｉｄｅ－Ｂｅａｒｉｎｇ Ｗａｓｔｅｓ ｆｒｏｍ Ｍａｎｕｆａｃｔｕｒｅｄ－Ｇａｓ Ｐｌａｎｔｓ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅ＆Ｔｅｃｈｎｏｌｏｇｙ，１９９４，２８（１）：９９—１０６． ［８］Ｇａｕ Ｓ Ｈ，Ｃｈｏｗ Ｊ ｎ Ｌａｎｄｆｉｌｌ ｌｅａｃｈａｔｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｎｄ ｍｏｄｅｌｉｎｇ ｏｆ ｍｕｎｉｃｉｐａｌ ｓｏｌｉｄ ｗａｓｔｅｓ ｃｏｍｂｉｎｅｄ ｗｉｔｈ ｉｎｃｉｎｅｒａｔｅｄ ｒｅｓｉｄｕａｌｓ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｈａｚａｒｄｏｕｓ Ｍａｔｅｒｉａｌｓ，１９９８，５８（１）：２４９—２５９． ［９］Ｅｄｚｗａｌｄ Ｊ Ｋ，Ｔｏｂｉａｓｏｎ Ｊ Ｅ．Ｅｎｈａｎｃｅｄ ｃｏａｇｕｌａｔｉｏｎ：ＵＳ ｒｅｑｕｉｒｅｍｅｎｔｓ ａｎｄ ａ ｂｒｏａｄｅｒ ｖｉｅｗ［Ｊ］．Ｗａｔｅｒ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ，１９９９， ４０（９）ｔ ６３－７０． ［１Ｏ］Ｗａｎｇ Ｌ ｓ，Ｈｕ Ｈ Ｙ，Ｗａｎｇ Ｃ．Ｅｆｆｅｃｔ ｏｆ ａｍｍｏｎｉａ ｎｉｔｒｏｇｅｎ ａｎｄ ｄｉｓｓｏｌｖｅｄ ｏｒｇａｎｉｃ ｍａｔｔｅｒ ｆｒａｃｔｉｏｎｓ ｏｎ ｔｈｅ ｇｅｎｏｔｏｘｉｃｉｔｙ ｏｆ ｗａｓｔｅｗａｔｅｒ ｅｆｆｌｕｅｎｔ ｄｕｒｉｎｇ ｃｈｌｏｒｉｎｅ ｄｉｓｉｎｆｅｃｔｉｏｎ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅ＆Ｔｅｃｈｎｏｌｏｇｙ，２００７，４１（１）：１６０—１６５． ［１１］Ｈｅ Ｐ Ｊ，Ｘｕｅ Ｊ Ｆ，Ｓｈａｏ Ｉ Ｍ，ｅｔ ａ１．Ｄｉｓｓｏｌｖｅｄ ｏｒｇａｎｉｃ ｍａｔｔｅｒ（ＤＯＭ）ｉｎ ｒｅｃｙｃｌｅｄ ｌｅａｃｈａｔｅ ｏｆ ｂｉｏｒｅａｃｔｏｒ ｌａｎｄｆｉｌｌ口］．Ｗａｔｅｒ Ｒｅｓｅａｒｃｈ， ２００６，４０（７）：１４６５—１４７３． ［１２］Ｆａｎ Ｈ Ｊ，Ｓｈｕ Ｈ Ｙ，Ｙａｎｇ Ｈ Ｓ，ｅｔ ａ１．Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｌａｎｄｆｉｌｌ ｌｅａｃｈａｔｅｓ ｉｎ ｃｅｎｔｒａｌ Ｔａｉｗａｎ［Ｊ］．Ｓｃｉｅｎｃｅ ｏｆ ｔｈｅ Ｔｏｔａｌ Ｅｎｖｉｒｏｎｍｅｎｔ， ２００６，３６（３）：２５—３７． ［１３］王怡中．二氧化钛悬浆体系中８种染料的太阳光催化氧化降解［Ｊ］．催化学报，２０００，２１（４）：３２７—３３１． ［编辑］ 李启栋