技术领域本发明涉及用于污染物的直接或分布数字控制系统领域，具体为一种用于湿法脱硫系统运行工况分析的定性方法。
背景技术：
目前，污染源自动监控系统研究与建设还主要处于污染源“末端监控”阶段，“末端监控”是指直接从并且只从企业的排污口进行数据采集与监控。由于数据采集器误差和人为造假等因素存在，“末端监控”不能保证数据的准确性和真实性，不能说清楚排污总量。虽然已经有研究注意到了现有污染源自动监控系统的数据准确性和真实性问题，并分析了其中的原因，但还未能提出一个完整完善的解决方案。
技术实现要素：
为了克服现有技术的缺陷，提供一种的，本发明公开了一种用于湿法脱硫系统运行工况分析的定性方法。本发明通过如下技术方案达到发明目的：一种用于湿法脱硫系统运行工况分析的定性方法，应用标准化模块以定性与定量相结合方式分析企业末端数据真实性的过程，定性分析能够实现对排污异常数据的报警，其特征是：按如下步骤依次实施：Ⅰ.工况采集数据进入标准化模块，首先对工况采集数据进行定性分析，判断其可用性(只有正常工况数据才能通过模型算出正确的末端结果数)，定性分析的过程：过滤异常数据，对其中的信号采集异常导致的畸变数据进行剔除；利用工况参数有效性校验方法判别工况设备运行是否正常；Ⅱ.通过定性分析后的正常工况数据，直接输出在线监测结果；A.数据预处理A.1数据缺失判断单点数据缺失判断，用数据质量点判断。如果数据质量点为timeout，则为数据缺失；A.2满屏跳计算一段时间的波动程度指标index和平均波动与均值ratio的比值(测试时间长度暂定为10分钟，测试间隔为1分钟)公式为：波动程度指标index=Σ|xi-xi-1|n(F‾),]]>ratio=(Fx‾-1)×100%,]]>其中，n是点的个数，是点的平均振幅，即把数据分成长度为m个点的若干段，求其平均振幅，F＝|xmax-xmin|，当index>0.05且ratio>20％，认为是满屏跳，其数值初步确定，经过测试可能会修改；A.3数据跳变/吹扫处理即滤波实时数据因为外界环境、吹扫等原因导致的数据跳变，系统将对判定并剔除。出现数据跳变的点主要包含入口烟气流量、出口烟气流量、入口SO2浓度、出口SO2浓度、脱硫效率。目前采用的判定方法为：采用1hour计算一次，每次取数据片段为本小时和前1h作为数据样本，进行2h数据的方差计算，数据小周期初步设定为5分钟，每1分钟进行前5分钟的均值计算，再设定方差范围，如果数据点或者数据片段超过设定的方差范围，则认为数据跳变，此时，还需要对数据进行范围比较，如果数据在有效性范围内，则不进行处理，认为正常，如果不在范围内，则需对数据进行剔除；如果数据跳变长度小于1min，则直接剔除，不报警；如果数据跳变长度大于1min，剔除，并对数据进行差值计算，输出数据跳变起始和结束时间；A.4恒定值判断目前有两种形式的恒定值：一种是原始测点故障，在DCS中已经是恒定值，这种情况，采用判断数据长周期差异的方式进行数据判断；一种是数据库插值导致的恒定值，这种情况，采用取数据点状态判断；在实际应用中，因不知道恒定值到底是哪种形式，则结合两种方式，首先查询判断时段数据的质量点，如果是timeout点，则直接输出数据缺失，如果数据质量点为good，则进入第二步判断，采用数据长周期差异性的方式进行判断；长周期差异性方法：计算一段时间的方差，测试时间定为10min，测试间隔为1min，若方差小于0.01，则说明这段时间内的数值是恒定值；(其数值初步确定，对于不同机组，会有相应的数据驯化和调整。)B.定性判断B.1数据有效性范围校验数据有效性范围基于脱硫工艺设计、锅炉和脱硫性能试验、脱硫运行经验，确定了不同装机容量下石灰石湿法脱硫工艺的关键因子合理性范围，数据有效性范围校验1hour检验一次，检验该小时的实时值；实时值与范围比较，如果不在范围内，则进行标记，如果累计有15min超过范围，认为该小时数据超限，该小时数据不可用；B.2关联度判断关联度判断进行三步进行计算，首先对实时数据进行去噪处理，然后进行参数长周期关联度计算，第三步进行陡升陡降时的精细化关联度计算；首先对数据进行去噪处理以减少采样值波动的影响，采用去噪的方法为：每秒钟计算一次，每次计算前5min的均值，存储为最新样本；长周期计算是对一天数据大趋势的检验，检验两个因子之间大趋势是否一致；当大趋势不在关联度范围内，再进行数据样本中陡升陡降片段的抓取，确定陡升陡降片段之后，按照时间周期1hour，计算该段时间参数之间的关联度；陡升陡降抓取的方式是，对因子小周期均值进行计算，计算该因子在本周期与上一周期的均值是否大于10％。如果大于10％，则认为陡升陡降；去噪处理之后，再进行关联度的计算，首先进行两个参数的关联度计算，每个因子都和其他因子做关联度计算，然后建立各因子之间的关联度矩阵，引进了对相关程度的度量——Pearson相关系数来定量分析两变量的线性相关性，Pearson相关系数又叫相关系数或者线性相关系数，用字母r表示，由两个变量样本取值得到，是一个描述线性相关强度的量。其中-1<r<1，|r|表明两变量间相关的程度，r>0表示正相关，r<0表示负相关，r＝0表示零相关。|r|越接近于1，表明两变量相关程度越高，它们之间的关系越密切，相关系数用r公式为：rxy=Σi=1n(xi-x‾)(yi-y‾)Σi=1n(xi-x‾)2Σi=1n(yi-y‾)2,]]>采用关联度矩阵的方式对烟气侧各因子进行定性判断，每个因子都和其他因子做关联度计算，然后建立如下表所示的关联度矩阵：关联度矩阵首先确定每一组参数关联度的权重，然后结合单因子与其他因子关联度计算的结果，用权重计算的方式，进行结果计算，最终确定是哪个因子出现异常；这里面需要确定的是，我们首先认定机组负荷是正常的，以此为基本样本，这样做的目的是避免了因多个因子异常而导致的最终结果误判。但这种做法的前提是，机组负荷测点正常。测点异常的情况包括主机数据中断、恒定值。关联度矩阵的计算方法，以燃煤量为例，燃煤量关联度矩阵结果＝a×A1+b×B1+c×B2+d×B3+e×B4+f×B5其中，A1为机组负荷与燃煤量的关联度值，以此类推，a为机组负荷与燃煤量关联度的权重；1.负荷—燃煤量：负荷变大，燃煤量变大：正相关，2.负荷—机组总送风量：正相关，3.负荷—入口烟气流量：正相关，4.负荷—增压风机电流：正相关，5.负荷—引风机电流：正相关，6.燃煤量—机组总送风量：正相关，7.燃煤量—入口烟气流量：正相关，8.燃煤量—增压风机电流：正相关，9.燃煤量—引风机电流：正相关，10.机组总送风量—入口烟气流量：正相关，11.机组总送风量—增压风机电流：正相关，12.机组总送风量—引风机电流：正相关，13.入口烟气流量—出口烟气流量：正相关，且出口烟气流量>入口烟气流量，14.入口烟气流量—增压风机电流：正相关，15.入口烟气流量—引风机电流：正相关，16.增压风机电流—风机电流：正相关，B.3逻辑判断当目标因子出现不随其他因子变化时，再进入逻辑关系判断，对其关联的因子变化进行判断，确定其变化趋势，最终异常原因进行定性，B.3.1脱硫效率脱硫效率有三类异常，一类是跳变，用数据有效性范围进行判断；第二类是与出入口浓度无关，检验方式是用出入口浓度计算脱硫效率，然后与实测脱硫效率进行比对；第三类是脱硫效率为恒定值，此时出口浓度直接判定为作假。脱硫效率＝(入口SO2浓度-出口SO2浓度)/入口SO2浓度×100％，B.3.2pHB.3.2.1pH定性判断原因由于吸收塔脱硫能力计算原理主要依赖于物料平衡和化学反应平衡，物料平衡计算主要依赖于钙硫比，而在实际应用中供浆流量是间断性的，而且浆液密度计损坏率高，因此使用钙硫比进行计算吸收塔脱硫能力时误差较大。化学反应平衡计算主要依赖于液气比、吸收塔pH、原烟气SO2浓度，而这几个参数中问题比较多的就是pH，但相对供浆流量和供浆密度而言，其计算精度高。但是pH是影响脱硫效率的一个关键因素，需要重视。比如入口浓度2000mg/m3，液气比为15的条件下，不同pH下的脱硫效率如下表所示：pH4.64.74.84.95.05.15.25.35.45.55.65.75.8效率88.989.890.591.391.992.593.093.493.794.094.294.394.4B.3.2.2pH异常分类pH值存在三类异常，一是跳变偏低或者偏高；二是恒定值；三是不随其他因子变化，跳变偏低或者偏高用数据有效性范围进行校验，恒定值采用数据质量点结合数据差异性进行判断，不随其他因子变化主要存在当出口浓度升高时pH值不变和出口烟气流量升高时pH值不变；当pH存在两个测点即能获得pH1和pH2两个值时，当pH1和pH2均正常时，对pH1和pH2进行平均作为使用值；当pH1或者pH2中有一个异常，一个正常时，则使用正常的这个值；当pH1和pH2都异常时，则对pH值进行补遗；pH补遗方案如下：当pH异常<72hour时，用异常前1hour进行补遗；B.4燃煤量B.4.1出口烟气流量烟气流量趋势问题主要有以下四种情况：趋势一致、烟气流量不随负荷变化、烟气流量陡升/陡降、限值；通过因子关联度判断出趋势不一致时，再确定是烟气流量未升降还是负荷未升降，或者是反向升降，采用均值的方式确定时间点和变化趋势，陡升陡降判断：计算一段时间内的均值，测试时间长度10min和测试间隔长度1min，设本次测定的均值为m1，上次测定的均值为m2，波动幅度为t，当|m1/m2-1|>t(t＝0.1)，则认为该时间发生了陡升陡降；对于限值判断，首先确定关联度不一致，负荷上升、烟气流量不变，然后再根据三阶拟合公式检验出口流量是否不负荷逻辑关系，三阶拟合公式是拟合负荷和烟气流量之间的；三阶拟合公式：Y＝a+bX+cX2+dX3，其中：Y为负荷，X为烟气流量，a、b、c、d为系数；三阶拟合公式在系统预设阶段做的，选取企业运行较好的一段时间作为数据样本进行公式拟合，确定系数a、b、c、d。B.4.2原烟气流量模型校验实测原烟气流量经过定性判断后，确定为正常的，需要再和模型进行绝对值比对；如果异常，直接输出模型结果；在核算模型原烟气流量需要用的参数包含全硫份、空干基灰分、空干基水份、固定碳、燃煤量，而其中燃煤量需要进行定性判断，定性判断目前采用关联度矩阵的形式进行判断，如果正常，直接使用，如果异常，则使用机组负荷(热电厂使用主蒸汽流量)进行推算燃煤量，推算燃煤量的度电煤耗采用CEMS补数规范，小于1day的，用前一天和后一天的均值核算，大于1day的，用前720hour的有效数据进行补遗；模型原烟气流量与实测原烟气流量进行比对，如果实测原烟气流量在模型原烟气流量的±20％以内，则认为实测正确，输出实测原烟气流量，如果实测原烟气流量不在模型原烟气流量的±20％以内，则认为实测原烟气流量绝对值有误，输出模型原烟气流量，并标记实测原烟气流量与机组负荷逻辑不符；B.4.3净烟气流量模型校验实测净烟气流量首先根据定性判断，其是否正常，如果异常，则直接使用模型值，如果正常，则与模型净烟气流量进行比对，如果实测净烟气流量在模型净烟气流量的±20％以内，则认为实测正确，输出实测净烟气流量，如果实测净烟气流量不在模型净烟气流量的±20％以内，则认为实测净烟气流量绝对值有误，输出模型净烟气流量，并标记实测净烟气流量与机组负荷逻辑不符；B.5出口SO2浓度在脱硫运行时，出口SO2浓度从表观上分为两大类，分别为表观正常、表观异常，表观异常再细划分包含四类，浓度偏低、恒定值、限值、浓度偏高；浓度偏低情况：需要区分本身的使用煤的硫份低还是脱硫效果好。对于像伊敏这样的厂，使用的煤硫份为0.09，在不使用脱硫时出口浓度已经达标，此种情况，再经过脱硫，出口浓度降低为20mg/m3是正确的。鉴于此种情况，我们使用上报硫份对预判做分阶，当上报硫份>0.4时，出口浓度初步认定小于40mg/m3，为偏低，直接使用模型数据；如果>40mg/m3，则启用模型对出口浓度进行校验；当上报硫份<0.4时，直接使用模型对出口浓度进行校验；备注1：这里40mg/m3是按照入口浓度为2000mg/m3，脱硫效率为98％时定的。实际应用中，入口浓度采用实测浓度作为基准值，给定脱硫效率为98％，从而确定出口浓度在(入口浓度×2％)以下认为出口实测浓度偏低，直接标记。设限值：一般企业排出的气体中SO2的浓度超过一定数值会被处罚，因此当出口烟气SO2浓度的实测值超过该数值后，企业会把该数值改为等于或稍小于该数值，我们称该数值为限值。正因为如此，限值一般是样本点的最大值。但并不是所有的最大值都是限值。一般情况下，当最大值(允许有一定误差，即小于最大值程度很小的数也认为是最大值)出现的点的个数或持续的时间超过一定时间，认为该值就是限值，设定方法如下：1)求出该段时间点最大值；2)设定一个误差值，当该点与最大值的差值小于误差值，则认为该点为疑似限值处理；3)疑似限值的个数超过600，即10min，则认为这段时间内该数据作了限值处理；恒定值：使用长周期差异性进行确定；浓度超标情况：分为旁路打开和旁路未打开两种情况，旁路打开，增压风机停，则判定为脱硫停运导致的超标，如果增压风机开，则继续判断循环浆液泵是否开，开几台，如果全没开，则说明脱硫停运，如果有部分开或者全开，判定为部分脱硫；旁路未打开，则根据模型判断入口浓度是否升高，如果入口浓度升高，则理论上循环浆液泵电流、石灰石供给流量、石膏排除泵流量3个都应该升高，如果3个都不变或降低，则输出脱硫条件不满足；如果有部分升高，那么没升高的参数判定为单点异常参数；如果入口浓度不变或者降低，理论上循环浆液泵电流、石灰石供给流量、石膏排除泵流量3个都应该降低，如果3个都不变或者升高，则说明出口作假，如果有任何一个降低，则输出脱硫条件不满足，并判断参数不变的因子为单点异常因子；B.5.1原烟气SO2浓度模型校验实测SO2浓度检验方法采用上报硫份推算的方式进行检验，首先对上报硫份进行正太分布统计，确定硫份的分布情况，如果上报的新硫份不在范围内，则系统会提示上报硫份与以往硫份差异较大，请再次核对，但继续使用该硫份作为计算硫份；根据全硫份、空干基水份、空干基灰分、固定碳根据工艺基础表格查出模型原烟气SO2浓度，再与实测原烟气SO2浓度进行比对，如果实测原烟气SO2浓度>模型的100％～20％，则认为实测是对的，输出实测值；如果实测原烟气SO2浓度<模型的100％～20％，则使用模型原烟气SO2浓度，并输出硫份与实测原烟气SO2浓度逻辑不符。这里需要注意的是：浓度比较，由于上报的硫份和实测浓度都是企业给定的，我们优先考虑相信上报的硫份，但当两个之间存在逻辑不符时，从整体角度考虑，我们使用的原则是实测浓度高，使用实测浓度，模型浓度高使用模型浓度。B.5.2净烟气SO2浓度模型校验净烟气SO2浓度检验和核定的方式是采用原烟气SO2浓度和脱硫效率，原烟气SO2浓度核定在上节已经确定；脱硫效率核定方法主要是采用的液气比和钙硫比的方式核算塔区的去除能力，使用的是核算仪原烟气SO2浓度、核算仪原烟气流量、循环浆液泵额定流量、循环浆液泵电流、吸收塔pH来核算的；首先计算液气比，然后根据pH、液气比、原烟气SO2浓度查工艺基础表格进行确定脱硫效率；液气比计算方法：通过脱硫效率和原烟气SO2浓度计算出模型净烟气SO2浓度，然后与实测净烟气SO2浓度进行范围比对，如果实测净烟气SO2浓度>模型的100％～20％，则认为实测是对的，输出实测值；如果实测净烟气SO2浓度<模型的100％～20％，则使用模型净烟气SO2浓度；这里使用模型净烟气SO2浓度时，需要确定异常原因，两种原因导致的，一种是原烟气SO2浓度，一种是脱硫效率，如果原烟气浓度模型高于实测，则此处输出原因为原烟气浓度高，脱硫不能达到预期要求；如果原烟气浓度模型小于等于实测，则说明脱硫效率模型小于实测，此时输出塔区反应条件不足。本发明属于环保管理领域的污染源自动监控系统范畴，应用于火电厂的湿法脱硫系统运行情况(工况、过程)分析，为环保管理提供支持。本发明的基本原理是从电厂的DCS系统中采集机组、FGD、CEMS数据，通过工况前端进行数据的采集、存储和传输，再上传环保厅。另核算仪采集前端工况数据、直采CEMS数据，与企业服务系统交换企业上报数据，进行过程工况数据的检验定性，如果过程数据异常，则启动定量进行核定，确定理论排放数据，核算仪需要进行的是工况的排口数据与直采的排口数据进行比对，确定排口数据是否在DCS中被更改。核算仪定性和定量结果通过环保专网上传至中心平台。本发明应用标准化模块以定性与定量相结合方式分析企业末端数据真实性的过程，定性分析能够实现对排污异常数据的报警。附图说明图1是定性分析流程图；图2是恒定值判断的流程图；图3是pH值和脱硫效率的函数关系图；图4是本发明中当存在两个pH值测点时的处理流程图；图5是pH补遗方案的流程图；图6是烟气流量趋势判断时的三阶拟合公式确定系数的处理流程图；图7是本发明中原烟气流量的核算流程图；图8是本发明中净烟气流量的核算流程图；图9是原烟气SO2浓度核算方法的流程图；图10是净烟气SO2浓度核算方法的流程图。具体实施方式以下通过具体实施例进一步说明本发明。实施例1一种用于湿法脱硫系统运行工况分析的定性方法，具体的流程图如图1所示，按如下步骤依次实施：Ⅰ.工况采集数据进入标准化模块，首先对工况采集数据进行定性分析，判断其可用性(只有正常工况数据才能通过模型算出正确的末端结果数)，定性分析的过程：过滤异常数据，对其中的信号采集异常导致的畸变数据进行剔除；利用工况参数有效性校验方法判别工况设备运行是否正常；Ⅱ.通过定性分析后的正常工况数据，直接输出在线监测结果；本产品将标准化模块开发成主要污染物智能核算仪，应用于主要污染物智能核算系统建设中，实现利用具备“真实性、准确性”的数据服务于环保业务管理。A.数据预处理A.1数据缺失判断单点数据缺失判断，用数据质量点判断。如果数据质量点为timeout，则为数据缺失；A.2满屏跳计算一段时间的波动程度指标index和平均波动与均值ratio的比值(测试时间长度暂定为10分钟，测试间隔为1分钟)公式为：波动程度指标index=Σ|xi-xi-1|n(F‾),]]>ratio=(Fx‾-1)×100%,]]>其中，n是点的个数，是点的平均振幅，即把数据分成长度为m个点的若干段，求其平均振幅，F＝|xmax-xmin|，当index>0.05且ratio>20％，认为是满屏跳(其数值初步确定，经过测试可能会修改)；A.3数据跳变/吹扫处理即滤波实时数据因为外界环境、吹扫等原因导致的数据跳变，系统将对判定并剔除。出现数据跳变的点主要包含入口烟气流量、出口烟气流量、入口SO2浓度、出口SO2浓度、脱硫效率。目前采用的判定方法为：采用1hour计算一次，每次取数据片段为本小时和前1hour作为数据样本，进行2h数据的方差计算，数据小周期初步设定为5min，每1min进行前5min的均值计算，再设定方差范围，如果数据点或者数据片段超过设定的方差范围，则认为数据跳变，此时，还需要对数据进行范围比较，如果数据在有效性范围内，则不进行处理，认为正常，如果不在范围内，则需对数据进行剔除；如果数据跳变长度小于1min，则直接剔除，不报警；如果数据跳变长度大于1min，剔除，并对数据进行差值计算，输出数据跳变起始和结束时间；A.4恒定值判断目前有两种形式的恒定值：一种是原始测点故障，在DCS中已经是恒定值，这种情况，采用判断数据长周期差异的方式进行数据判断；一种是数据库插值导致的恒定值，这种情况，采用取数据点状态判断；在实际应用中，因不知道恒定值到底是哪种形式，则结合两种方式，首先查询判断时段数据的质量点，如果是timeout点，则直接输出数据缺失，如果数据质量点为good，则进入第二步判断，采用数据长周期差异性的方式进行判断；如图2所示；长周期差异性方法：计算一段时间的方差，测试时间定为10min，测试间隔为1min，若方差小于0.01，则说明这段时间内的数值是恒定值；其数值初步确定，对于不同机组，会有相应的数据驯化和调整。B.定性判断B.1数据有效性范围校验数据有效性范围基于脱硫工艺设计、锅炉和脱硫性能试验、脱硫运行经验，确定了不同装机容量下石灰石湿法脱硫工艺的关键因子合理性范围，数据有效性范围校验1hour检验一次，检验该小时的实时值；实时值与范围比较，如果不在范围内，则进行标记，如果累计有15min超过范围，认为该小时数据超限，该小时数据不可用；B.2关联度判断关联度判断进行三步进行计算，首先对实时数据进行去噪处理，然后进行参数长周期关联度计算，第三步进行陡升陡降时的精细化关联度计算；首先对数据进行去噪处理以减少采样值波动的影响，采用去噪的方法为：每秒钟计算一次，每次计算前5min的均值，存储为最新样本；长周期计算是对一天数据大趋势的检验，检验两个因子之间大趋势是否一致；当大趋势不在关联度范围内，再进行数据样本中陡升陡降片段的抓取，确定陡升陡降片段之后，按照时间周期1hour，计算该段时间参数之间的关联度；陡升陡降抓取的方式是，对因子小周期均值进行计算，计算该因子在本周期与上一周期的均值是否大于10％。如果大于10％，则认为陡升陡降；去噪处理之后，再进行关联度的计算，首先进行两个参数的关联度计算，每个因子都和其他因子做关联度计算，然后建立各因子之间的关联度矩阵，引进了对相关程度的度量——Pearson相关系数来定量分析两变量的线性相关性，Pearson相关系数又叫相关系数或者线性相关系数，用字母r表示，由两个变量样本取值得到，是一个描述线性相关强度的量。其中-1<r<1，|r|表明两变量间相关的程度，r>0表示正相关，r<0表示负相关，r＝0表示零相关。|r|越接近于1，表明两变量相关程度越高，它们之间的关系越密切，相关系数用r公式为：rxy=Σi=1n(xi-x‾)(yi-y‾)Σi=1n(xi-x‾)2Σi=1n(yi-y‾)2,]]>采用关联度矩阵的方式对烟气侧各因子进行定性判断，每个因子都和其他因子做关联度计算，然后建立关联度矩阵，如下表所示：关联度矩阵首先确定每一组参数关联度的权重，然后结合单因子与其他因子关联度计算的结果，用权重计算的方式，进行结果计算，最终确定是哪个因子出现异常；这里面需要确定的是，我们首先认定机组负荷是正常的，以此为基本样本，这样做的目的是避免了因多个因子异常而导致的最终结果误判。但这种做法的前提是，机组负荷测点正常。测点异常的情况包括主机数据中断、恒定值。关联度矩阵的计算方法，以燃煤量为例，燃煤量关联度矩阵结果＝a×A1+b×B1+c×B2+d×B3+e×B4+f×B5其中，A1为机组负荷与燃煤量的关联度值，以此类推；a为机组负荷与燃煤量关联度的权重。1.负荷—燃煤量：负荷变大，燃煤量变大，正相关，2.负荷—机组总送风量：正相关，3.负荷—入口烟气流量：正相关，4.负荷—增压风机电流：正相关，5.负荷—引风机电流：正相关，6.燃煤量—机组总送风量：正相关，7.燃煤量—入口烟气流量：正相关，8.燃煤量—增压风机电流：正相关，9.燃煤量—引风机电流：正相关，10.机组总送风量—入口烟气流量：正相关，11.机组总送风量—增压风机电流：正相关，12.机组总送风量—引风机电流：正相关，13.入口烟气流量—出口烟气流量：正相关，且出口烟气流量>入口烟气流量，14.入口烟气流量—增压风机电流：正相关，15.入口烟气流量—引风机电流：正相关，16.增压风机电流—风机电流：正相关。B.3逻辑判断当目标因子出现不随其他因子变化时，再进入逻辑关系判断，对其关联的因子变化进行判断，确定其变化趋势，最终异常原因进行定性，B.3.1脱硫效率脱硫效率有三类异常，一类是跳变，用数据有效性范围进行判断；第二类是与出入口浓度无关，检验方式是用出入口浓度计算脱硫效率，然后与实测脱硫效率进行比对；第三类是脱硫效率为恒定值，此时出口浓度直接判定为作假。脱硫效率＝(入口SO2浓度-出口SO2浓度)/入口SO2浓度×100％，B.3.2pHB.3.2.1pH定性判断原因由于吸收塔脱硫能力计算原理主要依赖于物料平衡和化学反应平衡，物料平衡计算主要依赖于钙硫比，而在实际应用中供浆流量是间断性的，而且浆液密度计损坏率高，因此使用钙硫比进行计算吸收塔脱硫能力时误差较大。化学反应平衡计算主要依赖于液气比、吸收塔pH、原烟气SO2浓度，而这几个参数中问题比较多的就是pH，但相对供浆流量和供浆密度而言，其计算精度高。但是pH是影响脱硫效率的一个关键因素，需要重视。比如入口浓度2000mg/m3，液气比为15的条件下，不同pH下的脱硫效率如下表所示：pH4.64.74.84.95.05.15.25.35.45.55.65.75.8效率88.989.890.591.391.992.593.093.493.794.094.294.394.4趋势如图3所示。B.B.2pH异常分类pH值存在三类异常，一是跳变偏低或者偏高；二是恒定值；三是不随其他因子变化，跳变偏低或者偏高用数据有效性范围进行校验，恒定值采用数据质量点结合数据差异性进行判断，不随其他因子变化主要存在当出口浓度升高时pH值不变和出口烟气流量升高时pH值不变；当pH存在两个测点即能获得pH1和pH2两个值时，处理流程如图4所示，当pH1和pH2均正常时，对pH1和pH2进行平均作为使用值；当pH1或者pH2中有一个异常，一个正常时，则使用正常的这个值；当pH1和pH2都异常时，则对pH值进行补遗；pH补遗方案如图5所示，具体为：当pH异常<72hour时，用异常前1hour进行补遗；B.4燃煤量B.4.1出口烟气流量烟气流量趋势问题主要有以下四种情况：趋势一致、烟气流量不随负荷变化、烟气流量陡升/陡降、限值，通过因子关联度判断出趋势不一致时，再确定是烟气流量未升降还是负荷未升降，或者是反向升降，采用均值的方式确定时间点和变化趋势，陡升陡降判断：计算一段时间内的均值，测试时间长度10min和测试间隔长度1min，设本次测定的均值为m1，上次测定的均值为m2，波动幅度为t，当|m1/m2-1|>t(t＝0.1)，则认为该时间发生了陡升陡降；对于限值判断，首先确定关联度不一致，负荷上升、烟气流量不变，然后再根据三阶拟合公式检验出口流量是否不负荷逻辑关系，三阶拟合公式是拟合负荷和烟气流量之间的；三阶拟合公式：Y＝a+bX+cX2+dX3，其中：Y为负荷，X为烟气流量，a、b、c、d为系数；三阶拟合公式在系统预设阶段做的，选取企业运行较好的一段时间作为数据样本进行公式拟合，确定系数a、b、c、d。对于不同情况的处理方式如图6所示；B.4.2原烟气流量模型校验原烟气流量的核算流程图如图7所示；实测原烟气流量经过定性判断后，确定为正常的，需要再和模型进行绝对值比对；如果异常，直接输出模型结果；在核算模型原烟气流量需要用的参数包含全硫份、空干基灰分、空干基水份、固定碳、燃煤量，而其中燃煤量需要进行定性判断，定性判断目前采用关联度矩阵的形式进行判断，如果正常，直接使用，如果异常，则使用机组负荷(热电厂使用主蒸汽流量)进行推算燃煤量，推算燃煤量的度电煤耗采用CEMS补数规范，小于1day的，用前一天和后一天的均值核算，大于1day的，用前720hour的有效数据进行补遗；模型原烟气流量与实测原烟气流量进行比对，如果实测原烟气流量在模型原烟气流量的±20％以内，则认为实测正确，输出实测原烟气流量，如果实测原烟气流量不在模型原烟气流量的±20％以内，则认为实测原烟气流量绝对值有误，输出模型原烟气流量，并标记实测原烟气流量与机组负荷逻辑不符；B.4.3净烟气流量模型校验净烟气流量的核算方法如图8所示；实测净烟气流量首先根据定性判断，其是否正常，如果异常，则直接使用模型值，如果正常，则与模型净烟气流量进行比对，如果实测净烟气流量在模型净烟气流量的±20％以内，则认为实测正确，输出实测净烟气流量，如果实测净烟气流量不在模型净烟气流量的±20％以内，则认为实测净烟气流量绝对值有误，输出模型净烟气流量，并标记实测净烟气流量与机组负荷逻辑不符；B.5出口SO2浓度在脱硫运行时，出口SO2浓度从表观上分为两大类，分别为表观正常、表观异常，表观异常再细划分包含四类，浓度偏低、恒定值、限值、浓度偏高；浓度偏低情况：需要区分本身的使用煤的硫份低还是脱硫效果好。对于像伊敏这样的厂，使用的煤硫份为0.09，在不使用脱硫时出口浓度已经达标，此种情况，再经过脱硫，出口浓度降低为20mg/m3是正确的。鉴于此种情况，我们使用上报硫份对预判做分阶，当上报硫份>0.4时，出口浓度初步认定小于40mg/m3，为偏低，直接使用模型数据；如果>40mg/m3，则启用模型对出口浓度进行校验；当上报硫份<0.4时，直接使用模型对出口浓度进行校验；备注1：这里40mg/m3是按照入口浓度为2000mg/m3，脱硫效率为98％时定的。实际应用中，入口浓度采用实测浓度作为基准值，给定脱硫效率为98％，从而确定出口浓度在(入口浓度×2％)以下认为出口实测浓度偏低，直接标记。设限值：一般企业排出的气体中SO2的浓度超过一定数值会被处罚，因此当出口烟气SO2浓度的实测值超过该数值后，企业会把该数值改为等于或稍小于该数值，我们称该数值为限值。正因为如此，限值一般是样本点的最大值。但并不是所有的最大值都是限值。一般情况下，当最大值(允许有一定误差，即小于最大值程度很小的数也认为是最大值)出现的点的个数或持续的时间超过一定时间，认为该值就是限值，设定方法如下：1)求出该段时间点最大值；2)设定一个误差值，当该点与最大值的差值小于误差值，则认为该点为疑似限值处理；3)疑似限值的个数超过600，即10min，则认为这段时间内该数据作了限值处理；恒定值：使用长周期差异性进行确定；浓度超标情况：分为旁路打开和旁路未打开两种情况，旁路打开，增压风机停，则判定为脱硫停运导致的超标，如果增压风机开，则继续判断循环浆液泵是否开，开几台，如果全没开，则说明脱硫停运，如果有部分开或者全开，判定为部分脱硫；旁路未打开，则根据模型判断入口浓度是否升高，如果入口浓度升高，则理论上循环浆液泵电流、石灰石供给流量、石膏排除泵流量3个都应该升高，如果3个都不变或降低，则输出脱硫条件不满足；如果有部分升高，那么没升高的参数判定为单点异常参数；如果入口浓度不变或者降低，理论上循环浆液泵电流、石灰石供给流量、石膏排除泵流量3个都应该降低，如果3个都不变或者升高，则说明出口作假，如果有任何一个降低，则输出脱硫条件不满足，并判断参数不变的因子为单点异常因子；B.5.1原烟气SO2浓度模型校验原烟气SO2浓度的核算方法流程如图9所示。实测SO2浓度检验方法采用上报硫份推算的方式进行检验，首先对上报硫份进行正太分布统计，确定硫份的分布情况，如果上报的新硫份不在范围内，则系统会提示上报硫份与以往硫份差异较大，请再次核对，但继续使用该硫份作为计算硫份；根据全硫份、空干基水份、空干基灰分、固定碳根据工艺基础表格查出模型原烟气SO2浓度，再与实测原烟气SO2浓度进行比对，如果实测原烟气SO2浓度>模型的100％～20％，则认为实测是对的，输出实测值；如果实测原烟气SO2浓度<模型的100％～20％，则使用模型原烟气SO2浓度，并输出硫份与实测原烟气SO2浓度逻辑不符。这里需要注意的是：浓度比较，由于上报的硫份和实测浓度都是企业给定的，我们优先考虑相信上报的硫份，但当两个之间存在逻辑不符时，从整体角度考虑，我们使用的原则是实测浓度高，使用实测浓度，模型浓度高使用模型浓度。B.5.2净烟气SO2浓度模型校验净烟气SO2浓度核算方法的流程如图10所示。净烟气SO2浓度检验和核定的方式是采用原烟气SO2浓度和脱硫效率，原烟气SO2浓度核定在上节已经确定；脱硫效率核定方法主要是采用的液气比和钙硫比的方式核算塔区的去除能力，使用的是核算仪原烟气SO2浓度、核算仪原烟气流量、循环浆液泵额定流量、循环浆液泵电流、吸收塔pH来核算的；首先计算液气比，然后根据pH、液气比、原烟气SO2浓度查工艺基础表格进行确定脱硫效率；液气比计算方法：通过脱硫效率和原烟气SO2浓度计算出模型净烟气SO2浓度，然后与实测净烟气SO2浓度进行范围比对，如果实测净烟气SO2浓度>模型的100％～20％，则认为实测是对的，输出实测值；如果实测净烟气SO2浓度<模型的100％～20％，则使用模型净烟气SO2浓度；这里使用模型净烟气SO2浓度时，需要确定异常原因，两种原因导致的，一种是原烟气SO2浓度，一种是脱硫效率，如果原烟气浓度模型高于实测，则此处输出原因为原烟气浓度高，脱硫不能达到预期要求；如果原烟气浓度模型小于等于实测，则说明脱硫效率模型小于实测，此时输出塔区反应条件不足。某机组某段时间内脱硫停运时出口SO2浓度设定限值，限值在2057mg/m3左右。通过定性分析中逻辑判断本例为出口SO2浓度的四类异常中设定限值，即当最大值(允许有一定误差，即小于最大值程度很小的数也认为是最大值)出现的点的个数或持续的时间超过一定时间，我们就认为该值就是限值。方法如下：1)求出该段时间点最大值(为2057mg/m3)；2)设定一个误差值，当该点与最大值的差值小于误差值，则认为该点为疑似限值处理；3)疑似限值的个数超过600，即10分钟，则认为这段时间内该数据作了限值处理。某机组某段时间内出口烟气流量设限值为2000km3/h。而限值解除后，其平均负荷在450MW，其出口烟气流量均值在1200km3/h，其流量明显偏低。通过定性分析中逻辑判断本例为出口烟气流量的四种异常中设定限值。