目前由于国内造纸产业布局、原料结构、企业规模不合理，生产工艺、技术设备的过时，碱回收率的低下和污染物排放标准偏低等问题，造纸废水治理已成为工业污染防治的焦点、热点和难点.造纸行业废水排放量大、色度高、纤维悬浮物多、耗氧量大及可生化性差，是浓度较高的有机废水，排出的废水对整个水体造成了严重的污染，对生态环境造成了极大的破坏，对人们的健康和生态环境的保护构成极大的威胁.所以，造纸行业废水的治理是当前水处理领域急需解决的问题之一，造纸行业排放限值的制定对该行业污染物的减排具有关键性的作用.随着清河流域经济的迅速发展，清河流域面临的污染问题日益严峻，由于缺少成熟的污染物排放控制体系，水质污染治理任务十分紧迫，通过清河流域总量控制管理是解决此问题的一种有效手段.本文选取清河流域沿岸的典型造纸厂，结合其排污特点，采用基于Petri网理论的OPMSE软件，对造纸厂水污染物BAT(best available technology)处理工艺的削减能力进行了仿真计算，探讨研究制定造纸行业污染物排放限值的方法，从而严格控制造纸行业排放的水污染物对清河流域水质的影响，以期改善清河流域水质.

自欧盟执行委员会提出IPPC指令至今，欧盟已有能源、钢铁、有色、化工和造纸等27个行业的最佳可行技术(BAT)参考文件(BREF)被相继颁布实施，此外6个跨行业的BREFs也提出了相应的BAT^[1].

Jensen^[6]引入了着色Petri网，将标识赋予多种颜色以示区分，将系统的层次关系通过不同颜色来表示，减少了Petri网模型的状态空间数量，增加了建模的灵活性.Wang等^[7]将Petri网与面向对象的思想相融合，提出了面向对象Petri网(OOPN)，从而使建模的灵活性和模型的模块化能力得到提高.Meng^[8]在面向对象Petri网的基础上对企业的可重构制造系统进行了探究，并将时间约束与分层次建模相结合.面向对象技术的引入使模型的可重用性得到了保障，便于企业应对快速变化的市场需求.

在柔性制造系统方面，Petri网被用于建立FMS模型，成为当今研究热点之一^[2].以Petri网的初始标识作为树根，通过变迁的触发所产生的标识作为树枝，通过变迁的逐步触发得到最终的标识，从而搜索Petri网模型的可达树^[3]；利用A^*算法来搜索Petri网的可达树，通过估价函数来决定应该触发的变迁，以此类推，直到出现最终标识，可避免Petri网模型整个可达树的产生^[4]；对全局搜索采用回溯搜索方法，对局部搜索采用A^*算法，从而避免在寻优过程中陷入局部收敛的问题^[5].

使用OPMSE的分析技术可以得到系统的结构特征和性能特征，借助OPMSE强大的描述能力可以描述一个复杂系统，结合其他技术可以形成系统的、不同层次的描述，OPMSE使用图形化的建模方式，具有严密的理论基础和强大的描述能力^[9].

在面向对象Petri网的理论基础上，通过使用OPMSE对开原造纸厂水污染物的BAT处理工艺进行模拟仿真.主要研究了基于OPMSE的仿真算法，在结合开原造纸厂BAT技术的基础上，通过设定BAT处理工艺中各个处理工段对污染物去除率的区间范围，对开原造纸厂水污染物的BAT处理工艺进行多次仿真计算.在使用OPMSE进行仿真计算时，用区间计算方法替代传统的单一数值计算方法，即对区间范围内的数字进行随机选取并代入计算，然后经过多次仿真计算，并将所有计算结果作为样本，估计总体均值的置信区间，使得出的结果更具有科学性.在经过仿真计算之后，根据计算结果分析开原造纸厂经过BAT削减处理后的水污染物排放质量浓度，从而确定造纸行业水污染物直接排放限值.主要分为以下几步：

1) 典型污染行业具体企业调研、统计.

2) 仿真系统的分析.仿真系统包括仿真软件、仿真模型、数据库等方面.首先，采用OPMSE软件，利用编程加图示逻辑的过程建立仿真模型.数据则采用独立计算方式，即单个模型仅采用本模型所需的数据.

3) 仿真系统的建模.根据开原造纸厂具体的生产流程，结合具体的排污特性，确定污水处理BAT工艺流程，通过计算机绘制流程图、编程等确定污染物处理工艺的仿真计算过程，从而得到开原造纸厂经过BAT削减处理后的污染物排放质量浓度.

4) 仿真系统的计算和数据分析.基于OPMSE，通过区间计算的方法，进行模拟仿真，并对计算出的数据进行分析，依据仿真结果得出结论，从而制定造纸行业水污染物的直接排放限值.

纸张生产过程中，废水主要来源：制备工段废液，蒸煮废液，筛选工段废液，漂白工段废液.在该工艺模拟流程图中，主要对污水COD，BOD₅，NH₃-N三种物质的质量浓度进行模拟.由于实际生产过程中，每日所产生污染物的量不同，导致污水的初始污染指标不尽相同，经历史监测数据查找，模拟中设置污水的COD质量浓度为1 000~5 000 mg/L，BOD₅质量浓度为500~1 500 mg/L，NH₃-N质量浓度为5~15 mg/L；历史监测数据显示污水经处理后COD的质量浓度为15.17~30.55 mg/L，BOD₅的质量浓度为5.12~8.19 mg/L，NH₃-N 的质量浓度为4.03~7.95 mg/L.造纸废水处理工艺基本流程为：粗格栅过滤—细格栅过滤—沉淀池—冷却池—均衡池—调制池—上流式厌氧污泥反应器—曝气池—二沉池—混合反应池—终沉池.废水处理工艺模拟流程图，即造纸生产工艺中废水处理工段细化如图 1所示.

图1(Figure 1)

图1 造纸厂废水处理工艺 Fig. 1 The waste water treatment process in paper mill

由于实际处理过程中，处理效率不可能为一个固定的数值，因此，模拟中设定沉淀池对NH₃-N的处理效率为0~5%；上流式厌氧污泥反应器对COD的处理效率为75%~80%，对BOD₅的处理效率为85%~88%；设定曝气池对COD的处理效率为65%~70%，对BOD₅的处理效率为92%~97%，对NH₃-N的处理效率为30%~35%；设定混合反应池对COD的处理效率为80%~85%.对于含有多个随机变量的系统而言，其统计结果的方法是非常前沿的.选取开原造纸厂其中某一次仿真，从开始到结束如表 1所示.

表1(Table 1)

表1 单次模拟仿真结果 Table 1 A single simulation results mg· L-1 初始质量浓度 仿真最终结果 COD BOD5 NH3-N COD BOD5 NH3-N 2 118.1 620.5 6.94 69.51 12.73 3.30

表1 单次模拟仿真结果 Table 1 A single simulation results

通过对清河流域的开原造纸厂分析研究，对基于OPMSE的BAT水处理组合进行了模拟仿真，且通过模拟仿真得到了最终排放污染物COD，BOD₅，NH₃-N的质量浓度.图 2~图 4给出了计算仿真系统的结果.

图2(Figure 2)

图2 出水COD质量浓度仿真结果 Fig. 2 Simulation results of mass concentration of COD

图3(Figure 3)

图3 出水BOD5质量浓度仿真结果 Fig. 3 Simulation results of mass concentration of BOD5

图4(Figure 4)

图4 出水NH3-N质量浓度仿真结果 Fig. 4 Simulation results of mass concentration of NH3-N

对于开原造纸厂100次的仿真结果，对其均值、标准差、最大值以及最小值进行了统计，结果见表 2.

表2(Table 2)

表2 造纸厂BAT水处理计算仿真模型的仿真数据分析表 Table 2 Data analysis of paper mill's BAT water treatment calculation simulation model mg·L-1 指标 仿真次数 均值 标准差 最大值 最小值 COD 100 37.965 14.74 70.47 12.79 BOD5 100 7.53 3.33 14.51 2.62 NH3-N 100 6.81 2.01 10.22 3.33

表2 造纸厂BAT水处理计算仿真模型的仿真数据分析表 Table 2 Data analysis of paper mill's BAT water treatment calculation simulation model

根据表 2的数值，可以得出开原造纸厂不同污染指标在用正态分布统计时不同置信水平的置信区间，如表 3所示.

表3(Table 3)

表3 造纸厂BAT水处理出水污染物质量浓度置信区间 Table 3 Confidence interval of pollutants mass concentration of paper mill’s BAT water treatment mg·L-1 置信水平/% 90 95 99 COD (35.52，40.41) (35.04，40.89) (34.09，41.84) BOD5 (6.97，8.08) (6.87，8.19) (6.65，8.40) NH3-N (6.48，7.15) (6.41，7.21) (6.28，7.34)

表3 造纸厂BAT水处理出水污染物质量浓度置信区间 Table 3 Confidence interval of pollutants mass concentration of paper mill’s BAT water treatment

从表 3可以看出，BAT水处理组合对该厂的污染物去除的总体效果较好，出水水质满足《辽宁省地方标准污水综合排放标准》(DB21/1627—2008)直接排放的水质要求^[10].表 2仿真结果中的最小值指出了通过最大化处理效率所能达到的水质，即通过人为提高水处理的效率和去除率之后，所能达到的出水污染指标质量浓度.建议所有造纸行业在采用该BAT处理组合的基础上，提高各处理方法的工作效率，达到更好的污染物减排效果和最低的出水污染指标质量浓度.表 4给出了理想状况下最佳和最差的出水水质的污染指标质量浓度.

表4(Table 4)

表4 BAT水处理出水最佳及最差水质 Table 4 Best and worst water quality of BAT water treatment mg·L-1 污染指标 COD BOD5 NH3-N 最佳质量浓度 12.79 2.62 3.33 最差质量浓度 70.47 14.51 10.22

表4 BAT水处理出水最佳及最差水质 Table 4 Best and worst water quality of BAT water treatment

OPMSE模型是含有多个随机变量的仿真系统，纵观历史监测数据与模拟仿真结果，可以验证该模型灵敏度较高、稳定性较强，通过输入变量的改变，能及时得到相应的反馈，较好地控制该模型，进一步证明了该模型较为实用.

历史监测数据均在99%置信区间范围内，根据仿真结果的置信区间、历史监测数据以及现有的标准进行对比^[10]，结果见表 5.

表5(Table 5)

表5 结果与现有排放限值的对比 Table 5 Comparison between results and existing emission limits mg·L-1 污染指标 COD BOD5 NH3-N 99%置信区间 (34.09，41.84) (6.65，8.40) (6.28，7.34) 历史监测数据 (15.17，30.55) (5.12，8.19) (4.03，7.95) 直接排放标准 50 10 8

表5 结果与现有排放限值的对比 Table 5 Comparison between results and existing emission limits

综合上述仿真结果、历史监测数据以及和现有标准的对比，初步拟定基于BAT的清河流域造纸行业污染物排放限值一份，见表 6.

表6(Table 6)

表6 基于BAT的清河流域造纸行业污染物排放限值 Table 6 Emission limits of the paper-making industry in Qinghe basin based on BAT mg·L-1 污染指标 COD BOD5 NH3-N 直接排放限值 45 9 8

表6 基于BAT的清河流域造纸行业污染物排放限值 Table 6 Emission limits of the paper-making industry in Qinghe basin based on BAT

1) 对BAT削减能力的仿真计算得出了造纸行业污水经过BAT处理之后COD,BOD₅,NH₃-N的质量浓度正态分布均值在不同置信水平下的置信区间.

2) 经过仿真可以得到BAT处理之后的最佳出水质量浓度和最差出水质量浓度，最佳出水的各指标质量浓度可以为污染物削减潜力和企业减排的空间给出指导意义，而最差出水质量浓度也可以为污染源预警提供指导.

3) 仿真结果表明，造纸行业BAT对于BOD₅的去除效果较好，COD的去除效果相对较差，企业在现有基础上应该适当改进处理工艺来提高处理工艺对COD的去除.根据仿真结果并结合现有的排放限值，拟定清河流域造纸行业污染物排放限值，直接排放限值为ρ(COD)=45 mg/L,ρ(BOD₅)=9 mg/L,ρ(NH₃-N) =8 mg/L.