一种基于敏感性分析和差分进化算法的参数自动校准方法

摘要

本发明涉及一种基于敏感性分析和差分进化算法的活性污泥模型的参数自动校准方法，该方法包括：基于活性污泥模型在线或离线仿真模拟平台模拟污水厂工艺得到相应的出水模拟值；基于敏感性分析的校准参数自动选择模块，利用上述模拟平台，根据设定的参数个数或阈值从众多活性污泥模型参数中依次选择确定需要校准的参数；设定关注的出水指标的可接受误差范围，从而建立校准的损失函数；参数自动校准模块，利用上述模拟平台与设定的损失函数，对选择的参数通过误差反馈进行校准，以达到设定的目标值；自动校准激活模块，通过设定误差限或校准周期，对平台上的水厂的参数在必要时进行校准操作。该方法具有设计合理、操作简便、可视化强、性能可靠稳定、能够被广泛应用的特点。

说明书

技术领域

本发明涉及污水处理领域，尤其涉及基于国际水协活性污泥模型的水厂仿真模拟模型的参数校准，并将该方法嵌入水厂仿真在线或者离线模拟平台，周期或按需对水厂仿真模拟模型的参数进行自动校准过程。

背景技术

建立活性污泥工艺模型后，可结合水质和运行数据对此期间的实际运行进行模拟预测。理想情况下，模拟预测值应与实测值一致或基本接近。但实际上，应用活性污泥模型缺省参数及运行数据的模拟预测结果可能与实测值之间存在一定差异，称之为模拟偏差。造成模拟偏差的原因比较复杂，按其来源可分为：1)活性污泥模型缺省参数与实际值之间的差异造成的模拟偏差。活性污泥模型缺省参数在特定条件测定，这一特定条件可能与实际运行条件不符合，从而在此条件下测定的模型参数与实际情况不符；2)进出水物料偏差，包括流量参数不准确、进水水质划分错误等；3)运行数据中的错误造成的模拟偏差，包括系统边界的定义错误(没有包括全部的物料进口和出口)、取样的错误及测量误差等。

应用活性污泥工艺模型的基础在于模拟值与实测值的一致性。因此，当模拟值与实测值之间的差异超出一定的容许范围之后，就必须对工艺模型进行检查和修正。目前对于工艺模型的校正，多数研究集中在利用敏感性分析等手段考察单因素对模型的影响，且按常规顺序对缺省参数逐一修正。这种方法的弊端是无法衡量各缺省参数之间对模型的联动影响，因此往往通过反复多轮校正从而达到模型的准确性，既耗时又会对校正精度造成一定程度影响。

因此，有必要通过敏感性分析和多参数优化结合的方式，对造成模型偏差的关键因素进行联动参数修正。

发明内容

针对现有污水厂仿真模拟软件参数校准耗时耗力，需要比较有经验的人员进行反复校准的缺点，本发明的目的在于：提供一种基于敏感性分析和差分进化算法的参数自动校准方法，省去人工校准参数的繁琐，节约仿真模拟软件的参数率定时间，其具有设计合理、操作简便、可视化强、性能可靠稳定、能够被广泛推广使用的特点。

为了达到上述目的，本发明采用如下技术方案：

该方法基于以国际水协活性污泥模型为基本模型构建的在线或离线仿真平台，模拟水厂运行情况，计算出水指标。对仿真模拟经验较少的人员通过利用敏感性分析，对相应的出水指标针对不同参数依次进行模拟计算，将对该值具有显著影响(超过一定阈值或敏感性排序中的前若干个)的参数予以记录，加入待校准参数中，而对于经验丰富的仿真模拟人员则支持自定义校准参数设定。由于水厂仪表以及模型本身的计算误差，通过经验设定关注的出水指标的可接受误差范围，建立用以校准的损失函数，将该损失函数推送给参数自动校准模块，利用上述模拟平台与设定的损失函数，对选择的参数通过误差反馈进行校准，至达到设定的目标值时止。由于季节、水质等不可抗力因素的影响，校准参数具有时效性，在一段时间后某些或全部参数将不适用于当前时段的水厂运行模拟，因此必要时需要重新校准参数，为达到此目的，加入了自动校准激活模块，该模块通过设定误差限(计算关注的目标值与预测值间的误差)或校准周期，对平台上的水厂的参数在必要时反复进行自动校准操作。其详细步骤如下：

前期条件：以国际水协活性污泥模型为基础的在线或离线仿真模拟平台、对特定水厂确定各项出水指标的误差限(由于不同水厂出水标准、仪表等不同存在差异)、对指定校准参数的人员需要给定各个参数的变化范围否则可选择根据敏感性分析自动选择参数并自动校准参数。

1.使用国际水协活性污泥模型或指定校准参数的直接进入2，否则，通过(1)式遍历各个参数对出水水质的影响程度，通过排序选取前10％的所有参数构成自动校准的参数向量进入2进行优化；

其中

2.假设通过1确定的需要自动校正的参数总个数为N，则随机(按照一致分布)生成第一代的P＝10*N个个体，第i个个体为

3.根据式(2)以及仿真模拟平台进行模拟计算每个个体的损失值，令第一代的最优解为

其中

4.对第G代种群中每个个体

其中1≤r1≠r2≠r3≠i≤P，为避免出现早熟现象使用自适应变异算子F：

F＝F

其中F

5.对第G代个体及其产生的变异个体V

其中j＝1,2,…,N,rand(0,1)表示按照一致分布在区间(0,1)内随机取值。

6.使用贪婪算法选择下一代个体：

7.根据(2)式以及仿真模拟平台模拟计算G+1代个体的损失值，并产生第G+1代的最优解

8.若G>G

9.利用仿真模拟平台对水厂水质数据不断模拟，若80％出出水水质模拟值与真实值的差连续20次超过相应的θ

附图说明

附图1是本发明一种基于敏感性分析和差分进化算法的整体流程示意图。

附图2是本发明一种基于敏感性分析和差分进化算法的某水厂基于稳态模拟后动态模拟TN人工校准模拟值、自动校准模拟值与实际值的对比图。

附图3是本发明一种基于敏感性分析和差分进化算法的某水厂基于稳态模拟后动态模拟TP人工校准模拟值、自动校准模拟值与实际值的对比图。

附图4是本发明一种基于敏感性分析和差分进化算法的某水厂基于稳态模拟后动态模拟NH

附图5是本发明一种基于敏感性分析和差分进化算法的某水厂基于稳态模拟后动态模拟COD人工校准模拟值、自动校准模拟值与实际值的对比图。

具体实施方式

下面结合附图1-5对本发明一种基于敏感性分析和差分进化算法的参数自动校准方法作以详细说明，并结合一具体实施例对该方法的效果作以说明。

一种基于敏感性分析和差分进化算法的参数自动校准方法，该方法基于以国际水协活性污泥模型为基本模型构建的在线或离线仿真平台，模拟水厂运行情况，计算出水指标。对仿真模拟经验较少的人员通过利用敏感性分析，对相应的出水指标针对不同参数依次进行模拟计算，将对该值具有显著影响(超过一定阈值或敏感性排序中的前若干个)的参数予以记录，加入待校准参数中，而对于经验丰富的仿真模拟人员则支持自定义校准参数设定。由于水厂仪表以及模型本身的计算误差，通过经验设定关注的出水指标的可接受误差范围，建立用以校准的损失函数，将该损失函数推送给参数自动校准模块，利用上述模拟平台与设定的损失函数，对选择的参数通过误差反馈进行校准，至达到设定的目标值时止。由于季节、水质等不可抗力因素的影响，校准参数具有时效性，在一段时间后某些或全部参数将不适用于当前时段的水厂运行模拟，因此必要时需要重新校准参数，为达到此目的，加入了自动校准激活模块，该模块通过设定误差限(计算关注的目标值与预测值间的误差)或校准周期，对平台上的水厂的参数在必要时反复进行自动校准操作。

以某基于活性污泥法的污水厂仿真模拟模型进行稳态模拟为例，对该方法步骤进行详细说明：

1.根据仿真模拟人员经验以及敏感性分析结果(根据敏感性分析获得需校准的参数上下界)，影响某污水厂出水氨氮、硝氮和总磷的活性污泥模型(以ASM2d为例其他模型步骤完全相同，仅参数名称因模型而异)动力学系数、进水组分占比同出水水质可接受范围值见下表：

表1用于校准的参数名称、默认值及范围

2.通过1确定的需要自动校正的参数总个数为7，则随机(按照一致分布)生成第一代P＝70个个体，第i个个体为

3.根据式(8)以及仿真模拟平台进行稳态模拟计算每个个体的损失值，初始化

其中

4.对第G代种群中每个个体

其中1≤r1≠r2≠r3≠i≤P，为避免出现早熟现象使用自适应变异算子F：

F＝F

其中F

5.对第G代个体及其产生的变异个体V

在此CR取值为0.8，rand(0,1)表示在区间(0,1)内按照一致分布取随机值。

6.使用贪婪算法选择下一代个体：

7.根据(8)式以及仿真模拟平台稳态模拟计算G+1代个体的损失值，并产生第G+1代的最优解

8.若G>G

表2自动校准模拟值、人工校准模拟值与真实值对比

9.将输出的X