一种煤灰流动温度的预测方法
【专利摘要】本发明提供一种煤灰流动温度的预测方法，该方法包括将待测煤灰成分含量代入煤灰流动温度预测模型，得到待测煤灰流动温度的预测值。其中，通过以下步骤构建煤灰流动温度预测模型：分析入厂煤或混配煤的煤灰成分以及煤灰流动温度，将获得的多组分析数据构成样本点数据集，其中每组分析数据包括煤灰成分含量和对应的煤灰流动温度；基于所述样本点数据集，将煤灰成分含量作为自变量并且将煤灰流动温度作为因变量，采用偏最小二乘回归法建立煤灰流动温度预测模型。本发明解决了煤灰流动温度传统测定方法分析过程复杂、分析周期长等局限性。此外，本发明克服了多重相关性在系统建模中的不良作用，缩小了模型预测的误差，提高了模型的稳健性。
【专利说明】一种煤灰流动温度的预测方法
【技术领域】
[0001] 本发明涉及软测量【技术领域】，具体来说，涉及一种煤灰流动温度的预测方法。
【背景技术】
[0002] 煤灰流动温度是壳牌（Shell)气化炉温度调节控制的重要指标之一，Shell气化 炉的操作温度必须超过煤灰流动温度1〇〇?150°C，且要求入炉煤的煤质必须相对稳定，煤 灰流动温度波动在l〇〇°C以内。炉温过高或过低将会出现渣屏、渣口堵渣等现象。另外，煤 灰流动温度直接关系到SGC (合成气冷却器）入口温度的调整，是影响气化炉负荷提高的重 要因素。因此，快速准确地获得入炉煤、原料煤、工艺混配煤的煤灰流动温度，对于Shell气 化炉的长周期稳定运行具有重要意义。
[0003] 现阶段，可以按照国标《GB/T219-2008煤灰熔融性的测定方法》来进行煤灰流动温 度分析。煤样须经过破碎、研磨、筛分、混匀、干燥、灰化、灰锥制作、灰锥干燥、仪器分析等过 程，单个煤样的分析周期为8?10小时，当样品较多时，分析周期将会超过12小时，从而无 法满足Shell气化炉的生产需要。
[0004] 根据煤灰成分经验关系式对灰熔点进行预测是较早应用的数学方法，但通过经验 关系式计算灰熔点的误差较大，平均相对预测误差超过7%。此外，随着数学及计算机软件 技术的不断发展，近年来，BP神经网络、支持向量机、遗传算法等建模技术已在灰熔点预测 中得到广泛应用，预测准确度也大幅提高。但灰熔点预测模型大多以煤灰软化温度为输出 变量，且均停留在实验室研究阶段，在工业化生产中尚未见报道。
[0005] 综上所述，当前还没有一种能够快速且准确地预测煤灰流动温度的技术，以满足 Shell气化炉实际生产的需要。
【发明内容】
[0006] 本发明的目的是提供一种准确性高、分析速度快的煤灰流动温度预测方法。根据 本发明的一个实施例，提供一种构建煤灰流动温度模型的方法，包括：
[0007] 步骤1)、分析入厂煤或混配煤的煤灰成分以及煤灰流动温度，将获得的多组分析 数据构成样本点数据集，其中每组分析数据包括煤灰成分含量和对应的煤灰流动温度；
[0008] 步骤2)、基于所述样本点数据集，将煤灰成分含量作为自变量并且将煤灰流动温 度作为因变量，采用偏最小二乘回归法建立煤灰流动温度预测模型。
[0009] 在一个实施例中，在步骤2)中包括：
[0010] 步骤21)、在所述样本点数据集中逐步提取自变量和因变量的主成分，对所述煤灰 流动温度预测模型进行拟合，使得所述煤灰流动温度预测模型的拟合度以及主成分对自变 量和因变量的累计解释能力达到最大；
[0011] 步骤22)、检查所述煤灰流动温度预测模型的拟合度以及主成分对自变量和因变 量的累计解释能力；如果不满足预定精度要求，则查找并移除所述样本点数据集中的离群 点，返回步骤21)重新对该模型进行拟合，进一步提高所述煤灰流动温度预测模型的拟合 度。
[0012] 在一个实施例中，用如下公式表示所述煤灰流动温度预测模型：
[0013] Y = e+a1X1+a2X2+. . . +anXn,
[0014] 其中，Y表不煤灰流动温度，X1JfXn分别表不n个煤灰成分含量，为通过 偏最小二乘回归法得到的回归系数，e为通过偏最小二乘回归法得到的参数。
[0015] 在一个实施例中，所述煤灰成分包括 Si02、A1203、Fe203、Ti0 2、CaO、MgO、S03、P205、 K2O 和 Na2O。
[0016] 在一个实施例中，在步骤I)中，采用灰熔点分析仪分析入厂煤或混配煤的煤灰流 动温度；并且采用X-荧光光谱仪、加载煤灰成分标准工作曲线来分析入厂煤或者混配煤的 煤灰成分。
[0017] 在一个实施例中，在步骤1)之前还包括采用下列步骤建立煤灰成分标准工作曲 线.
[0018] 选择各元素含量具有预定梯度分布的多个煤灰标准样品，经化学法确定所述多个 煤灰标准样品中煤灰成分含量；
[0019] 将所述多个煤灰标准样品分别压制成煤灰标准样片；
[0020] 利用X-荧光光谱仪对煤灰标准样片进行扫描，根据各煤灰成分脉冲信号强度建 立煤灰成分标准工作曲线。
[0021] 根据本发明的一个实施例，还提供一种煤灰流动温度的预测方法，包括：
[0022] 将待测煤灰成分含量代入根据上述方法所构建的煤灰流动温度预测模型，得到待 测煤灰流动温度的预测值。
[0023] 本发明具有以下优点及突出性效果：
[0024] ⑴准确度高：截止目前为止，使用本发明提供的煤灰流动温度预测方法已累计报 出煤灰流动温度预测数据1658批次，与实测数据相比较，流动温度的平均相对预测误差和 最大相对预测误差分别为1. 04%和3. 20%，满足实际生产的需要。
[0025] ⑵分析速度快：使用本发明的预测方法来预测煤灰流动温度可使入厂煤灰熔点分 析时间缩短至4小时，入炉煤灰熔点分析时间缩短至3小时，极大地提高了分析效率，为气 化炉工艺参数的快速调整提供了有效的数据支撑。
[0026] ⑶预测方法先进：将煤灰成分的分析数据输入煤灰流动温度预测模型中即可快速 预测出煤灰的流动温度，是一种先进的软测量技术，在工业化生产使用过程中尚未见有相 关报道。
[0027] ⑷节能降耗优：利用本发明的预测方法来预测煤灰流动温度可节省在实测过程中 消耗的物资和人力成本。
[0028] (5)应用范围广：使用本发明提供的方法可构建出各种煤质煤灰流动温度的预测模 型，适用于电厂锅炉运行、粉煤灰电解铝和混配煤工艺生产等领域。
【专利附图】
【附图说明】
[0029] 图1是根据本发明一个实施例的煤灰流动温度预测方法的流程图；
[0030] 图2是X-荧光光谱仪的系统结构示意图；
[0031] 图3是根据本发明一个实施例的用于标识离群点的t/t2平面图；以及
[0032] 图4是采用本发明提供的方法得到的煤灰流动温度预测值与实测值的比较示意 图。
【具体实施方式】
[0033] 下面结合附图和【具体实施方式】对本发明加以说明。应当理解，此处所描述的具体 实施例仅用以解释本发明，并不用于限定本发明。
[0034] 根据本发明的一个实施例，提供一种煤灰流动温度的预测方法。参考图1，该方法 包括以下步骤：
[0035] 第1步：建立煤灰成分标准工作曲线
[0036] 首先，选择各元素含量有一定梯度分布的多个煤灰标准样品，经化学法反复测定 定标，确定煤灰中各成分的含量，并且将这些煤灰标准样品压制成标准样片。
[0037] 接着，调整X-荧光光谱仪的工作条件，利用X-荧光光谱仪对上述煤灰标准样片进 行扫描，检查每个元素通道的分析角度、确定背景、计算峰和背景的测量时间、检查每个通 道的脉冲高度分布和脉冲处理条件，将煤灰标准样品各成分含量输入X-荧光光谱仪软件 中，根据各煤灰成分脉冲信号强度建立煤灰成分标准工作曲线。
[0038] 图2示出了本发明使用的X-荧光光谱仪的系统结构，该X-荧光光谱仪包括X-射 线管、滤光片、准直器、分光晶体、检测器、分析器（未示出）等部件。其测量原理为：x-射 线管发射的高能一次X射线经滤光片过滤后照射样品，样品中的元素被激发而发出特定波 长的X-荧光，X-荧光经准直器后将发散光汇聚成一束平行光，光束强度和分辨性提高后通 过晶体的衍射散布在空间不同的方位上，并由检测器跟踪测量，进行光电转换和能量分辨， 信号输入多道分析器（图中未示出）进行定性和定量分析。
[0039] 在一个优选的实施例中，可选择八个或者更多个煤灰标准样品来建立煤灰成分标 准工作曲线。
[0040] 在一个实施例中，煤灰标准样品的成分可包括Si02、A120 3、Fe203、Ti02、CaO、MgO、 S03、P205、K2O和Na2O,下文将以这十个成分为例展开描述。在其他实施例中，煤灰标准样品 的成分还可以包括MnO2和BaO。
[0041] 第2步：筛选用于建模煤灰流动温度预测模型的样本点。
[0042] 本步骤可分为以下子步骤来执行：
[0043] 2. 1、使用X-荧光光谱仪，加载煤灰成分标准曲线对入厂煤、混配煤煤灰成分进行 分析；并且根据GB/T219-2008煤灰熔融性测定方法，使用灰熔点分析仪分析上述煤的煤灰 流动温度。
[0044] 2. 2、从上面分析得到的实测数据中筛选出多组分析数据，其中每组分析数据包括 煤灰各成分（例如上文给出的十个成分：Si02、A1203、Fe 203、Ti02、CaO、MgO、S03、P20 5、K20、 Na2O)含量和对应的煤灰流动温度的实测数据。
[0045] 在一个优选的实施例中，可以从上面分析得到的多组分析数据中选择煤灰成分的 含量能够覆盖其各个数值范围的代表性数据作为样本点。例如，对于其中某一成分，选择其 含量在5 % -10 %范围内的5个实测数据、选择其含量在10 % -15 %范围内的5个实测数据， 以此类推。
[0046] 在一个优选的实施例中，可以选择至少100组有代表性的煤灰成分含量和煤灰流 动温度的分析数据作为样本点。
[0047] 第3步：建立煤灰流动温度预测模型。
[0048] 本发明采用偏最小二乘回归法来建立煤灰流动温度预测模型，偏最小二乘回归法 作为一种新的多元统计分析法，现已成功应用于分析化学，如紫外光谱、气相色谱和电分 析化学等。偏最小二乘回归法是一种数学优化技术，它通过最小化误差的平方和找到一组 数据的最佳函数匹配，用最简的方法求得一些绝对不可知的真值，而令误差平方之和为最 小。具体来说，假设有P个因变量h，y2, .... yp与m个自变量X1, x2, .... Xm。为研究因变 量和自变量的统计关系，观测n个样本点，由此构成了自变量与因变量的数据表X = Ix1, X2,--X1J和Y = {yi，y2,--yp}。偏最小二乘回归法分别在X与Y中提取出主成分h和 U1 (其中L是X1, x2,--xm的线性组合，U1是y2,--y p的线形组合）。在提取这两个 主成分时，为了回归分析的需要，有以下两个要求：
[0049] (1)&和U1应尽可能大地携带他们各自数据表中的变异信息；
[0050] (2) h和U1的相关程度能够达到最大。
[0051] 这两个要求表明，h和U1应尽可能好地代表数据表X和Y，同时自变量的主成分心 对因变量的主成分U1又有最强的解释能力。
[0052] 在第一个主成分h和U1被提取后，偏最小二乘回归法分别实施X对h的回归以 及Y对U1的回归。如果回归方程已经达到预期的精度，则算法终止；否则，将利用X被h解 释后的残余信息以及Y被^解释后的残余信息进行第二轮的主成分提取。如此往复，直到 能达到一个较满意的精度为止。若最终对X提取了 r个主成分t2, ....&，偏最小二乘 回归法将建立yp ... yp与tp t2,--^的的回归，然后再表示--yp与原自变量的回 归方程式，即建模的偏最小二乘回归方程式。
[0053] 在本步骤中，以煤灰成分含量为自变量，以煤灰流动温度为因变量，应用偏最小二 乘回归法来建立煤灰流动温度预测模型。例如，用所选的多组样本点数据中的煤灰成分 Si02、A1203、Fe203、Ti0 2、CaO、MgO、S03、P205、K20、Na 2O 的含量作为自变量 X，以灰熔点分析仪 实测的煤灰流动温度Ft为因变量Y，建立单因变量的偏最小二乘回归模型作为煤灰流动温 度预测模型。包括以下子步骤：
[0054] 3. 1、在上述样本点中逐步提取自变量X和因变量Y的主成分，对煤灰流动温度预 测模型进行拟合，使得拟合得到的预测模型的整体拟合度、主成分对自变量和因变量的累 计解释能力达到最大。
[0055] 在一个实施例中，当提取到主成分tn+1时，如果对自变量X的累计解释能力、对因 变量Y的累计解释能力或预测模型整体拟合度不再增加或者增加不明显，则提取到主成分 、时，煤灰流动温度预测模型的预测能力最佳（即此时拟合得到的预测模型的整体拟合度、 主成分对自变量和因变量的累计解释能力达到最大）。
[0056] 3. 2、当不满足预期的精度要求时，对模型进行精制，使拟合度以及主成分对自变 量和因变量的累计解释能力进一步提高，达到预期的要求。
[0057] 模型精制的一个例子可包括：查找样本点中的离群点，移除样本点中的离群点，并 且在移除了离群样本点的样本点数据集中再次逐步提取自变量X和因变量Y的主成分，重 新对煤灰流动温度预测模型进行拟合。
[0058] 3. 3、得到煤灰流动温度预测模型。
[0059] 在一个实施例中，煤灰流动温度预测模型可以表示为：
[0060] Y = e+a1X1+a2X2+. . . +anXn
[0061] 其中，Y表示流动温度，X^Xn分别表示n个煤灰成分含量，ai-a n为通过上述 偏最小二乘回归法得到的回归系数，e为通过上述偏最小二乘回归法得到的参数。
[0062] 第4步：对于流动温度未知的待测煤灰，使用X-荧光光谱仪对其进行扫描，根据煤 灰成分标准工作曲线计算待测煤灰各成分（例如，上文描述的十个成分）的含量。接着，将 待测煤灰的各成分含量代入煤灰流动温度预测模型中即可得出待测煤灰的流动温度。
[0063] 下文以褐煤为例，结合上述预测方法，具体描述褐煤煤灰流动温度的预测过程。
[0064] 1、建立褐煤煤灰成分标准工作曲线
[0065] 1)、选择各元素含量有一定梯度分布的八个褐煤煤灰标准样品，通过化学法确定 煤灰中 Si02、A1203、Fe203、Ti0 2、CaO、MgO、S03、P205、K20、Na 2O 十个成分的含量。将这八个煤 灰标准样品分别粉碎研磨制成200目粉样，然后与粘结剂混合，用硼酸作边，在30T压力下 保压30S，压制成标准样片。
[0066] 2)、调整X-荧光光谱仪的工作条件：光路介质选择真空、真空锁定时间为16s、延 迟时间为2s、准直器光栏孔径为27mm、样品类型为压片粉末、样品直径为50mm，各元素分析 条件的设置如表1所示。
[0067]表 1
【权利要求】
1. 一种构建煤灰流动温度模型的方法，包括： 步骤1)、分析入厂煤或混配煤的煤灰成分以及煤灰流动温度，将获得的多组分析数据 构成样本点数据集，其中每组分析数据包括煤灰成分含量和对应的煤灰流动温度； 步骤2)、基于所述样本点数据集，将煤灰成分含量作为自变量并且将煤灰流动温度作 为因变量，采用偏最小二乘回归法建立煤灰流动温度预测模型。
2. 根据权利要求1所述的方法，在步骤2)中包括： 步骤21)、在所述样本点数据集中逐步提取自变量和因变量的主成分，对所述煤灰流动 温度预测模型进行拟合，使得所述煤灰流动温度预测模型的拟合度以及主成分对自变量和 因变量的累计解释能力达到最大； 步骤22)、检查所述煤灰流动温度预测模型的拟合度以及主成分对自变量和因变量的 累计解释能力；如果不满足预定精度要求，则查找并移除所述样本点数据集中的离群点，返 回步骤21)重新对所述煤灰流动温度预测模型进行拟合。
3. 根据权利要求1或2所述的方法，其中，用如下公式表示所述煤灰流动温度预测模 型： Y = e+a1X1+a2X2+. . . +anXn 其中，Y表示煤灰流动温度，Xi、Xn分别表示n个煤灰成分含量，&1一 an为通过偏最 小二乘回归法得到的回归系数，e为通过偏最小二乘回归法得到的参数。
4. 根据权利要求1或2所述的方法，其中，所述煤灰成分包括Si02、A1203、Fe 203、Ti02、 CaO、MgO、S03、P205、K20 和 Na20。
5. 根据权利要求1或2所述的方法，在步骤1)中，采用灰熔点分析仪分析入厂煤或混 配煤的煤灰流动温度。
6. 根据权利要求1或2所述的方法，在步骤1)中，采用X-荧光光谱仪、加载煤灰成分 标准工作曲线来分析入厂煤或者混配煤的煤灰成分。
7. 根据权利要求6所述的方法，其中，在步骤1)之前还包括采用下列步骤建立煤灰成 分标准工作曲线： 选择各元素含量具有预定梯度分布的多个煤灰标准样品，经化学法确定所述多个煤灰 标准样品中煤灰成分含量； 将所述多个煤灰标准样品分别压制成煤灰标准样片； 利用X-荧光光谱仪对煤灰标准样片进行扫描，根据各煤灰成分脉冲信号强度建立煤 灰成分标准工作曲线。
8. 根据权利要求7所述的方法，其中，所述多个煤灰标准样品的个数为8个或者更多 个。
9. 根据权利要求1或2所述的方法，其中，所述样本点数据集包括100组或者更多组分 析数据。
10. -种煤灰流动温度的预测方法，包括： 将待测煤灰成分含量代入根据权利要求1-9中任何一个所述的方法所构建的煤灰流 动温度预测模型，得到待测煤灰流动温度的预测值。
【文档编号】G06F17/50GK104408215SQ201410385010
【公开日】2015年3月11日 申请日期:2014年8月7日 优先权日:2014年8月7日
【发明者】马继生, 高秀娟, 张志明, 丁志刚, 王昌盛, 姜新其, 廖家势, 李宁, 黄晓伟, 张美林 申请人:大唐内蒙古多伦煤化工有限责任公司