一种高浓磨浆系统游离度的软测量方法.pdf

摘要
申请专利号：	CN201410768527.9	申请日：	2014.12.12
公开号：	CN104459089A	公开日：	2015.03.25
当前法律状态：	授权	有效性：	有权
法律详情：	授权\|\|\|实质审查的生效IPC(主分类):G01N33/34申请日:20141212\|\|\|公开
IPC分类号：	G01N33/34	主分类号：	G01N33/34
申请人：	东北大学
发明人：	周平; 赵梓丞; 王宏; 代凌旖
地址：	110819辽宁省沈阳市和平区文化路3号巷11
优先权：
专利代理机构：	沈阳亚泰专利商标代理有限公司21107	代理人：	许宇来
PDF下载：	PDF下载

内容摘要

一种高浓磨浆系统游离度的软测量方法涉及化学机械制浆生产过程高浓磨浆系统游离度的测量技术领域，具体地说是涉及一种基于粒子群算法优化ε-支持向量回归机(PSO-ε-SVR)的高浓磨浆系统游离度的软测量方法。本发明提供一种高浓磨浆系统游离度的软测量方法，该方法能够实现对高浓磨浆机出浆游离度的实时测量，提高人工测量的效率和精度，为工业实时控制和实现节能降耗提供技术支持。本发明方法由硬件平台及测量软件组成，其中硬件平台核心包括高浓磨浆系统三盘磨浆机、压力传感器、震动传感器、加速度传感器、流量传感器、分料螺旋等。

权利要求书

1.  一种高浓磨浆系统游离度的软测量方法，包括以下步骤：(1)辅助变量选择和模型输入变量确定，(2)PSO-ε-SVR软测量模型的训练和软测量实现。
(1)辅助变量选择和模型输入变量确定
需要软测量的高浓磨浆系统纸浆关键参数为游离度(ml)。根据生产工艺及机理以及可测可观和变量间相关性分析(因子分析、协方差分析)，选择i个软测量辅助变量(i＝9)，包括：
设定产量u₁(adt/d)；
磨浆机功率u₂(Mw)；
传动侧稀释水量u₃(l/min)；
非传动侧稀释水量u₄(l/min)；
非传动侧震动加速度u₅(％)；
非传动侧震动u₆(mm/s)；
非传动侧磨盘间隙u₇(mm)；
传动侧磨盘压力u₈(bar)；
非传动侧磨盘压力u₉(bar)。
确定上述9个软测量辅助变量为软测量模型的输入变量，游离度(ml)为软测量模型的输出变量。
(2)PSO-ε-SVR软测量模型的训练和软测量实现
(A)开始：变量初始化；
(B)选择进行软测量模型训练，转至(C)，读取软测量训练样本；选择高浓磨浆系统游离度软测量，转至(K)，读取软测量模型最优参数(惩罚因子C和核函数参数δ)_；
(C)读取软测量训练样本：从数据库导入或输入训练集T＝{(x₁,y₁),L,(x_l,y_l)}∈(Rⁿ×γ)^l。其中，软测量模型输入：x_i∈Rⁿ，i＝1,2,L,l，软测量模型输出：y_i∈γ＝R，i＝1,2,L,l，l为训练集样本数量；
(D)数据预处理：
采用尖峰滤波算法将高浓磨浆系统的尖峰异常数据进行剔除；9个辅助变量存在着量纲差，对每个变量分别进行中心标准化处理，即每类样本数据减去对应样本的平均值，然后除以其样本方差，即：
计算均值：

计算方差：

标准化：

其中，x_ij为第i类变量的第j个训练样本，l为训练样本数，为第i类变量的均值，为中心标准化后的训练样本，表示第i类变量的方差，表示第i类变量的标准差；
(E)确定ε-SVR软测量模型的初始化参数，包括：惩罚因子C和误差容限值ε，核函数Φ(g)的类型及其核函数参数δ。惩罚因子C取20，误差容限值ε取0.01，核函数取高斯径向基函数K(x,x')＝exp(-||x-x'||²/δ²)，其中核函数参数δ＝1；
(F)确定ε-SVR待优化参数惩罚因子C的搜索范围为[0,100]，核函数参数δ的搜索范围为[0,50]，粒子群个体的数目40，粒子的最大迭代寻优次数600，粒子的惩罚因子C迭代速度的最大值为25，粒子的核函数参数δ迭代速度的最大值为25，粒子的初始位置、初始速度、个体极值以及全局极值随机生成；
(G)ε-SVR模型训练：
ε-SVR计算高浓磨浆系统游离度的回归超平面方程y＝(ωgx)+b，其中ω为其法方向，b为其截距。
设模型样本的第j个目标输出为y_j，则ε-SVR将拟合问题等价转化为如下所述的二次规划问题：

其中，C为惩罚因子，ε为误差容限值，ξ_i和(简记为)为松弛变量。
求解式(4)中二次规划问题的解(ω,b,ξ)后，得到回归超平面方程：
y＝(ωgx)+b   (5)
将Lagrange函数引入上述问题，导出其对偶问题：

解得Lagrange乘子向量
计算b：选择位于开区间..中的α^(*)的分量α_j或若选到的是α_j，则有：

若选择的是则有：

构造回归方程：

其中，y为训练样本i的输出。(x_igx_k)是核函数Φ(g)，这里取高斯径向基函数K(x,x')＝exp(-||x-x'||²/δ²)，是对应的拉格朗日乘子的各个分量。
(H)计算适应度：
设定优化的目标函数，并计算适应度，对每组惩罚因子C和核函数参数δ的效果进行评价，通过相应的误差函数计算适应度。对于误差大的粒子，认为其适应度小，粒子的适应度f_p表示为：
f_p＝1/(e_p+1)   (10)
其中，e_p是回归方程的误差函数，即：

其中，是回归方程的预测输出，y_i是回归方程的目标输出；
更新每个粒子的速度和位置

其中，表示待更新粒子的速度，表示待更新粒子的位置，p_id表示待更新粒子的个体最优值，p_gd表示整个粒子群的全局最优值。w为惯性权重，c₁、c₂为加速度常数，rand₁和rand₂为[0,1]范围内的随机数；
对粒子p，如果更新后的适应度值大于原个体的最优值，则更新粒子的个体最优值p_id，即：
p_id＝f_p   (13)
如果粒子p的个体最优值p_id大于原来粒子群的全局最优值p_gd，则更新原来的全局最优值p_gd：
p_gd＝p_id   (14) 。

2.  根据权利要求1所述一种高浓磨浆系统游离度的软测量方法，其特征在于还包括以下步骤：
(I)建模效果评价：
判定建模效果是否满足要求。若满足，结束寻优，得到惩罚因子C和核函数参数δ的最优值，进入步骤(J)，保存ε-SVR惩罚因子C和核函数参数δ的最优参数；否则返回步骤(F)，继续迭代寻优，直到达到最大寻优代数，结束寻优；
(J)保存ε-SVR惩罚因子C和核函数参数δ的最优参数；
(K)读取ε-SVR惩罚因子C和核函数参数δ的最优参数；
(L)读取软测量模型的输入数据；
(M)将输入变量数据中心标准化处理后，使用ε-SVR惩罚因子C和核函数参数δ的最优参数进行高浓磨浆系统游离度在线软测量；
(N)游离度软测量结果显示：在软测量系统人机界面上显示游离度软测量的结果；
(O)游离度软测量结果保存：将本次软测量的输入输出数据保存到相应的数据库，供查询、修正使用；
(P)判断软测量是否结束，结束转至步骤(Q)，否则转至步骤(L)继续进行游离度软测量；
(Q)结束：完成游离度软测量。

3.  根据权利要求2所述一种高浓磨浆系统游离度的软测量方法，其特征在于还包括以下步骤：
定期将游离度离线间歇式数据输入到训练集中，更新游离度ε-支持向量回归机软测量模型。

说明书

一种高浓磨浆系统游离度的软测量方法
技术领域
本发明涉及化学机械制浆生产过程高浓磨浆系统游离度的测量技术领域，具体地说是涉及一种基于粒子群算法优化ε-支持向量回归机(PSO-ε-SVR)的高浓磨浆系统游离度的软测量方法。
背景技术
造纸过程由制浆和造纸两大环节串联而成。制浆环节的主要功能是从植物纤维原材料生产出具有特定形态的纤维，而造纸环节的功能主要是以特定形态的纤维为原料生产出各种纸制产品。制浆造纸工业在我国是重点能耗行业之一，生产过程中消耗大量热能、电能和水能。特别是制浆过程的能耗及其所产生的纸浆的质量直接关系到后续造纸环节的能耗及产品质量。
按现有工艺，制浆分为机械浆、化学浆和化学机械浆三大类。其中，化学机械浆生产过程主要包括：木片预处理，高浓磨浆，粗浆漂白，粗浆洗涤，低浓磨浆，渣浆磨浆，纤维回收，浆料浓缩脱水等过程。高浓磨浆作为化学机械浆生产的关键步骤，其用电量占到整个工厂用电的35％以上，决定着纸浆质量和生产能耗。
然而，高浓磨浆过程具有多变量、强耦合和非线性的特征，使得高浓磨浆过程的机理分析、建模存在很大困难。目前，研究的假设性强，得到的高浓磨浆模型缺乏通用性。此外，研究主要集中在低浓磨浆过程、单盘磨浆机和磨盘改进上，对于高浓磨浆系统三盘磨浆机的建模和分析还没有公开报道。
纤维网滤水性(游离度)是检测纸浆通过铜网的滤水程度快慢的一个指标。作为评价纸浆质量的关键指标之一，该指标不仅关系着最终纸张的质量(透光性，吸水性等)，同时与整个生产过程的能耗密切相关。为此，使用游离度来评价纤维质量十分合理。
公开资料表明，目前还没有传感器对于高浓磨浆系统游离度直接进行测量，主要依靠离线间歇人工测量。该测量方法存在着较大的时滞性和误差，进而无法实时对制浆生产进行控制，阻碍实现节能降耗和纸浆质量控制。为此，提出一种基于粒子群算法优化ε-支持向量机(PSO-ε-SVR)的高浓磨浆系统游离度软测量的技术。本发明可通过常规在线测量方法提供的模型输入数据，给出输入数据对应的在线估计值，为制浆生产过程的优化操作和优化运行提供关键质量指标，进而实现质量控制和节能降耗。
发明内容
为了克服高浓磨浆过程中游离度无法直接测量，离线人工测量时效性差、精度低等缺点，本发明提供一种高浓磨浆系统游离度的软测量方法，该方法能够实现对高浓磨浆机出浆游离度的实时测量，提高人工测量的效率和精度，为工业实时控制和实现节能降耗提供技术支持。
本发明方法由硬件平台及测量软件组成，其中硬件平台核心包括高浓磨浆系统三盘磨浆机、压力传感器、震动传感器、加速度传感器、流量传感器、分料螺旋等。
一种最优参数的支持向量回归机高浓磨浆系统游离度软测量包括以下步骤：(1)辅助变量选择和模型输入变量确定，(2)PSO-ε-SVR软测量模型的训练和软测量实现。
(1)辅助变量选择和模型输入变量确定
需要软测量的高浓磨浆系统纸浆关键参数为游离度(ml)。根据生产工艺及机理以及可测可观和变量间相关性分析(因子分析、协方差分析)，选择i个软测量辅助变量(i＝9)，包括：
设定产量u₁(adt/d)；
磨浆机功率u₂(Mw)；
传动侧稀释水量u₃(l/min)；
非传动侧稀释水量u₄(l/min)；
非传动侧震动加速度u₅(％)；
非传动侧震动u₆(mm/s)；
非传动侧磨盘间隙u₇(mm)；
传动侧磨盘压力u₈(bar)；
非传动侧磨盘压力u₉(bar)。
确定上述9个软测量辅助变量为软测量模型的输入变量，游离度(ml)为软测量模型的输出变量。
(2)PSO-ε-SVR软测量模型的训练和软测量实现
(A)开始：变量初始化；
(B)选择进行软测量模型训练，转至(C)，读取软测量训练样本；选择高浓磨浆系统游离度软测量，转至(K)，读取软测量模型最优参数(惩罚因子C和核函数参数δ)；
(C)读取软测量训练样本：从数据库导入或输入训练集T＝{(x₁,y₁),L,(x_l,y_l)}∈(Rⁿ×Υ)^l。其中，软测量模型输入：x_i∈Rⁿ，i＝1,2,L,l，软测量模型输出：y_i∈Υ＝R，i＝1,2,L,l，l为训练集样本数量；
(D)数据预处理：
采用尖峰滤波算法将高浓磨浆系统的尖峰异常数据进行剔除；由于选取的9个辅助变量存在着量纲差，对每个变量分别进行中心标准化处理，即每类样本数据减去对应样本的平均值，然后除以其样本方差，即：
计算均值：
x&OverBar;i=1lΣj=1lxij---(1)]]>
计算方差：
σxi2=1l-1Σj=1l(xij-x&OverBar;i)---(2)]]>
标准化：

其中，x_ij为第i类变量的第j个训练样本，l为训练样本数，为第i类变量的均值，为中心标准化后的训练样本，表示第i类变量的方差，表示第i类变量的标准差；
(E)确定ε-SVR软测量模型的初始化参数，包括：惩罚因子C和误差容限值ε，核函数Φ(g)的类型及其核函数参数δ。惩罚因子C取20，误差容限值ε取0.01，核函数取高斯径向基函数K(x,x′)＝exp(-||x-x'||²/δ²)，其中核函数参数δ＝1；
(F)确定ε-SVR待优化参数惩罚因子C的搜索范围为[0,100]，核函数参数δ的搜索范围为[0,50]，粒子群个体的数目40，粒子的最大迭代寻优次数600，粒子的惩罚因子C迭代速度的最大值为25，粒子的核函数参数δ迭代速度的最大值为25，粒子的初始位置、初始速度、个体极值以及全局极值随机生成；
(G)ε-SVR模型训练：
ε-SVR计算高浓磨浆系统游离度的回归超平面方程y＝(ωgx)+b，其中ω为其法方向，b为其截距。
设模型样本的第i个目标输出为y_i，则ε-SVR将拟合问题等价转化为如下所述的二次规划问题：
minw,b,ξ(*)12||ω||2+CΣi=1l(ξi+ξi*)s.t.((ωgxi)+b)-yi≤ϵ+ξi,i=1,L,l;yj-((ωgxi)+b)≤ϵ+ξi*,i=1,K,l;ξj*&GreaterEqual;0,i=1,L,l.---(4)]]>
其中，C为惩罚因子，ε为误差容限值，ξ_i和(简记为)为松弛变量。
求解式(4)中二次规划问题的解后，得到回归超平面方程：
y＝(ωgx)+b   (5)
将Lagrange函数引入上述问题，导出其对偶问题：
minα(*)&Element;R2l12Σi,j=1l(αi*-αi)(αj*-αj)(xigxj)+ϵΣi=1l(αi*+αi)-Σi=1lyi(αi*-αi)=0,s,t.Σi=1l(αi*-αi)=0,0≤αi(*)≤C,i=1,K,l.---(6)]]>
解得Lagrange乘子向量α(*)=(α1,α1*,L,αl,αl*)T.]]>
计算b：选择位于开区间(0,C)中的α^(*)的分量α_j或若选到的是α_j，则有：
b=yj-Σi=1l(αi*-αi)(xigxj)+ϵ---(7)]]>
若选择的是则有：
b=yk-Σi=1l(αi*-αi)(xigxk)-ϵ---(8)]]>
构造回归方程：
y=Σi=1l(αi*-αi)(xigxk)+b---(9)]]>
其中，y为训练样本输出。(x_igx_k)是核函数Φ(g)，这里取高斯径向基函数K(x,x')＝exp(-||x-x'||²/δ²)，是对应的拉格朗日乘子的各个分量。
(H)计算适应度：
设定优化的目标函数，并计算适应度，对每组惩罚因子C和核函数参数δ的效果进行评价，通过相应的误差函数计算适应度。对于误差大的粒子，认为其适应度小，粒子的适应度f_p表示为：
f_p＝1/(e_p+1)   (10)
其中，e_p是回归方程的误差函数，即：
ep=1lΣi=1l(y^i-yi)2---(11)]]>
其中，是回归方程的预测输出，y_i是回归方程的目标输出；
更新每个粒子的速度和位置
vidk=wvidk-1+c1rand1(pid-xidk-1)+c2rand2(pgd-xidk-1)xidk=xidk-1+vidk-1---(12)]]>
其中，表示待更新粒子的速度，表示待更新粒子的位置，p_id表示待更新粒子的个体最优值，p_gd表示整个粒子群的全局最优值。w为惯性权重，c₁、c₂为加速度常数，rand₁和rand₂为[0,1]范围内的随机数；
对粒子p，如果更新后的适应度值大于原个体的最优值，则更新粒子的个体最优值p_id，即：
p_id＝f_p   (13)
如果粒子p的个体最优值p_id大于原来粒子群的全局最优值p_gd，则更新原来的全局最优值p_gd：
p_gd＝p_id   (14)
作为一种优选方案：本发明还包括以下步骤：
(I)建模效果评价：
判定建模效果是否满足要求。若满足，结束寻优，得到惩罚因子C和核函数参数δ的最优值，进入步骤(J)，保存ε-SVR惩罚因子C和核函数参数δ的最优参数；否则返回步骤(F)，继续迭代寻优，直到达到最大寻优代数，结束寻优；
(J)保存ε-SVR惩罚因子C和核函数参数δ的最优参数；
(K)读取ε-SVR惩罚因子C和核函数参数δ的最优参数；
(L)读取软测量模型的输入数据；
(M)将输入变量数据中心标准化处理后，使用ε-SVR惩罚因子C和核函数参数δ的最优参数进行高浓磨浆系统游离度在线软测量；
(N)游离度软测量结果显示：在软测量系统人机界面上显示游离度软测量的结果；
(O)游离度软测量结果保存：将本次软测量的输入输出数据保存到相应的数据库，供查询、修正使用；
(P)判断软测量是否结束，结束转至步骤(Q)，否则转至步骤(L)继续进行游离度软测量；
(Q)结束：完成游离度软测量。
作为另一种优选方案：本发明还包括以下步骤：
定期将游离度离线间歇式数据输入到训练集中，更新游离度ε-支持向量回归机软测量模型。
本发明的一种高浓磨浆系统的游离度软测量方法具有以下技术效果：
本发明针对目前高浓磨浆过程中无法直接测量游离度的现状，在没有设备测量游离度的前提下，选取对高浓磨浆系统游离度有直接影响的9个操作变量和检测变量作为软测量模型的输入数据，实现对高浓磨浆系统游离度的软测量，其测量精度高，为实现整个造纸过程的节能降耗提供理论依据和模型支持。
附图说明：
图1是本发明选择的研究对象：高浓磨浆系统三盘磨浆机的结构及13个主要操作变量和检测变量
图2是本发明的一种基于粒子群优化ε-支持向量回归机(PSO-ε-SVR)的高浓磨浆系统游离度软测量算法的流程图
图3是本发明选择的高浓磨浆系统制浆生产过程
高浓磨浆过程中，经预处理后的木片经过传动侧分料螺旋和非传动侧分料螺旋，传动侧负载感应螺旋和非传动侧负载感应螺旋进入高浓磨浆机，经过高浓磨机磨浆后的浆料进入旋风分离器从而实现气浆分离，便形成了粗浆。
图4是本发明寻得最优参数对游离度进行软测量的效果图
具体实施方式
如图1所示，高浓磨浆系统三盘磨浆机由中间一个具有双齿面结构的定盘与两侧两个定盘构成双磨室结构。三盘磨浆机在高转速时运行平稳，不会出现动盘偏斜等问题，在不提高转速及增大磨盘直径的情况下，磨盘面积增加一倍，具有提高产量、改进磨浆质量、方便热能回收等优点。
本发明提供一种高浓磨浆系统三盘磨浆机游离度的软测量方法。该方法能够实现对高浓磨浆机出浆游离度的实时测量，提高人工测量的效率和精度，为工业实时控制和实现节能降耗提供技术支持。
参照图1、图2，一种基于粒子群算法优化ε-支持向量回归机(PSO-ε-SVR)的高浓磨浆系统游离度软测量的方法，其特征在于有常规测量系统、数据采集器、PSO-ε-SVR软测量软件以及运行软测量软件的计算机系统构成。流量计、压力机、震动计、加速度计等常规测量仪表安装于高浓磨浆系统的各个相应位置。数据采集器连接常规测量系统，并通过通信总线连接运行软测量软件的计算机系统。常规测量系统主要包括如下常规测量仪表：
二个流量计(春辉LW JY-LWGY)，分别用于在线测量高浓磨浆系统传动侧稀释水量，非传动侧稀释水量；
一个振动加速度传感器(YMK 151A100)，用于在线测量高浓磨浆系统非传动侧的震动；
一个加速度计(PCB SA51SC)，用于在线测量高浓磨浆系统非传动侧震动加速度的大小；
一个电涡流位移传感器(航振HZ891XL)，用于在线测量高浓磨浆系统非传动侧定盘与动盘之间的间隙；
二个压力计(GP50211)，用于在线测量高浓磨浆系统传动侧磨盘压力和非传动侧磨盘压力。
此外，常规测量系统还包括如下两个测量仪：
一个高浓磨浆系统设定产量测量仪：通过传动侧分料螺旋和非传动侧分料螺旋转速设定来设定高浓磨浆系统产量；
一个高浓磨浆系统磨浆机功率测量仪：通过各子系统能耗计量得到高浓磨浆系统磨浆机功率。
本发明实现方法包括：(1)辅助变量选择和模型输入变量确定；(2)PSO-ε-SVR软测量模型的训练和软测量实现。
(1)辅助变量选择和模型输入变量确定
需要软测量的高浓磨浆系统纸浆关键参数为游离度(ml)。根据生产工艺及机理以及可测可观和变量间相关性分析(因子分析、协方差分析)，选择i个(i＝9)软测量辅助变量包括：
设定产量u₁(adt/d)；
磨浆机功率u₂(Mw)；
传动侧稀释水量u₃(l/min)；
非传动侧稀释水量u₄(l/min)；
非传动侧震动加速度u₅(％)；
非传动侧震动u₆(mm/s)；
非传动侧磨盘间隙u₇(mm)；
传动侧磨盘压力u₈(bar)；
非传动侧磨盘压力u₉(bar)。
确定上述9个软测量辅助变量为为软测量模型的输入变量，游离度(ml)作为软测量模型的输出变量。
(2)PSO-ε-SVR软测量模型的训练和软测量实现
(A)开始：变量初始化；
(B)选择进行软测量模型训练，转至(C)，读取软测量训练样本；选择高浓磨浆系统游离度软测量，转至(K)，读取软测量模型最优参数(惩罚因子C和核函数参数δ)；
(C)读取软测量训练样本：从数据库导入或输入训练集T＝{(x₁,y₁),L,(x_l,y_l)}∈(Rⁿ×Υ)^l。其中，软测量模型输入：x_i∈Rⁿ，i＝1,2,L,l，软测量模型输出：y_i∈Υ＝R，i＝1,2,L,l，l为训练集样本数量；
(D)数据预处理：
采用尖峰滤波算法将高浓磨浆系统的尖峰异常数据进行剔除；由于选取的9个辅助变量存在着量纲差，对每个变量分别进行中心标准化处理，即每类样本数据减去对应样本的平均值，然后除以其样本方差，即：
计算均值：
x&OverBar;i=1lΣj=1lxij---(15)]]>
计算方差：
σxi2=1l-1Σj=1l(xij-x&OverBar;i)---(16)]]>
标准化：

其中，x_ij为第i类变量的第j个训练样本，l为训练样本数，为第i类变量的均值，为中心标准化后的训练样本，表示第i类变量的方差，表示第i类变量的标准差；
(E)确定ε-SVR软测量模型的初始化参数，包括：惩罚因子C和误差容限值ε，核函数Φ(g)的类型及其核函数参数δ。惩罚因子C取20，误差容限值ε取0.01，核函数取高斯径向基函数K(x,x')＝exp(-||x-x'||²/δ²)，其中核函数参数δ＝1；
(F)确定ε-SVR待优化参数惩罚因子C的搜索范围为[0,100]，核函数参数δ的搜索范围为[0,50]，粒子群个体的数目40，粒子的最大迭代寻优次数600，粒子的惩罚因子C迭代速度的最大值为25，粒子的核函数参数δ迭代速度的最大值为25，粒子的初始位置、初始速度、个体极值以及全局极值随机生成；
(G)ε-SVR模型训练：
ε-SVR计算高浓磨浆系统游离度的回归超平面方程y＝(ωgx)+b，其中ω为其法方向，b为其截距。
设模型样本的第i个目标输出为y_i，则ε-SVR将拟合问题等价转化为如下所述的二次规划问题：
minw,b,ξ(*)12||ω||2+CΣi=1l(ξi+ξi*)s.t.((ωgxi)+b)-yi≤ϵ+ξi,i=1,L,l;yj-((ωgxi)+b)≤ϵ+ξi*,i=1,K,l;ξj*&GreaterEqual;0,i=1,L,l.---(18)]]>
其中，C为惩罚因子，ε为误差容限值，ξ_i和(简记为)为松弛变量。
求解上述二次规划问题的解后，得到回归方程：
y＝(ωgx)+b   (19)
将Lagrange函数引入上述问题，导出其对偶问题：
minα(*)&Element;R2l12Σi,j=1l(αi*-αi)(αj*-αj)(xigxj)+ϵΣi=1l(αi*+αi)-Σi=1lyi(αi*-αi)=0,s,t.Σi=1l(αi*-αi)=0,0≤αi(*)≤C,i=1,K,l.---(20)]]>
解得Lagrange乘子向量α(*)=(α1,α1*,L,αl,αl*)T.]]>
计算b：选择位于开区间(0,C)中的α^(*)的分量α_j或若选到的是α_j，则有：
b=yj-Σi=1l(αi*-αi)(xigxj)+ϵ---(21)]]>
若选择的是则有：
b=yk-Σi=1l(αi*-αi)(xigxk)-ϵ---(22)]]>
构造回归方程：
y=Σi=1l(αi*-αi)(xigxk)+b---(23)]]>
其中，y为训练样本输出。(x_igx_k)是核函数Φ(g)，这里取高斯径向基函数K(x,x')＝exp(-||x-x'||²/δ²)，是对应的拉格朗日乘子的各个分量。
(H)计算适应度：
设定优化的目标函数，并计算适应度，对每组惩罚因子C和核函数参数δ的效果进行评价，通过相应的误差函数计算适应度。对于误差大的粒子，认为其适应度小，粒子的适应度f_p表示为：
f_p＝1/(e_p+1)   (24)
其中，e_p是回归方程的误差函数，即：
ep=1lΣj=1l(y^j-yj)2---(25)]]>
其中，是回归方程的预测输出，y_i是回归方程的目标输出；
更新每个粒子的速度和位置
vidk=wvidk-1+c1rand1(pid-xidk-1)+c2rand2(pgd-xidk-1)xidk=xidk-1+vidk-1---(26)]]>
其中，表示待更新粒子的速度，表示待更新粒子的位置，p_id表示待更新粒子的个体最优值，p_gd表示整个粒子群的全局最优值。w为惯性权重，c₁为加速度常数，rand₁和rand₂为[0,1]范围内的随机数；
对粒子p，如果更新后的适应度值大于原个体的最优值，则更新粒子的个体最优值p_id，即：
p_id＝f_p   (27)
如果粒子p的个体最优值p_id大于原来粒子群的全局最优值p_gd，则更新原来的全局最优值p_gd：
p_gd＝p_id   (28)
(I)建模效果评价：
判定建模效果是否满足要求。若满足，结束寻优，得到惩罚因子C和核函数参数δ的最优值，进入步骤(J)，保存ε-SVR惩罚因子C和核函数参数δ的最优参数；否则返回步骤(F)，继续迭代寻优，直到达到最大寻优代数，结束寻优；
(J)保存ε-SVR惩罚因子C和核函数参数δ的最优参数；
(K)读取ε-SVR惩罚因子C和核函数参数δ的最优参数；
(L)读取软测量模型的输入数据；
(M)将输入变量数据中心标准化处理后，使用ε-SVR惩罚因子C和核函数参数δ的最优参数进行高浓磨浆系统游离度在线软测量；
(N)游离度软测量结果显示：在软测量系统人机界面上显示游离度软测量的结果；
(O)游离度软测量结果保存：将本次软测量的输入输出数据保存到相应的数据库，供查询、修正使用；
(P)判断软测量是否结束，结束转至步骤(Q)，否则转至步骤(L)继续进行游离度软测量；
(Q)结束：完成游离度软测量。
作为优选的一种方案：所述软测量方法还包括以下步骤：
定期将游离度离线间歇式数据输入到训练集中，更新游离度ε-支持向量回归机软测量模型。
图4为一段时间内，高浓磨浆系统游离度软测量的效果，可以看出游离度的预测值与其实际值基本一致，误差比较小，且变化趋势基本一致。此外，本发明方法速度快、精度高、泛化能力强且有严格的数学解释，相比于其他方法有较高的优越性。因此本发明是一种具有很高实用价值的、低成本的高浓磨浆系统游离度计量手段。
可以理解的是，以上关于本发明的具体描述，仅用于说明本发明而并非受限于本发明实施例所描述的技术方案，本领域的普通技术人员应当理解，仍然可以对本发明进行修改或等同替换，以达到相同的技术效果；只要满足使用需要，都在本发明的保护范围之内。

资源描述

《一种高浓磨浆系统游离度的软测量方法.pdf》由会员分享，可在线阅读，更多相关《一种高浓磨浆系统游离度的软测量方法.pdf（18页珍藏版）》请在专利查询网上搜索。

1、(10)申请公布号 (43)申请公布日 (21)申请号 201410768527.9 (22)申请日 2014.12.12 G01N 33/34(2006.01) (71)申请人东北大学地址 110819 辽宁省沈阳市和平区文化路 3 号巷 11 (72)发明人周平赵梓丞王宏代凌旖 (74)专利代理机构沈阳亚泰专利商标代理有限公司 21107 代理人许宇来 (54) 发明名称一种高浓磨浆系统游离度的软测量方法 (57) 摘要一种高浓磨浆系统游离度的软测量方法涉及化学机械制浆生产过程高浓磨浆系统游离度的测量技术领域，具体地说是涉及一种基于粒子群算法优化-支持向量回归。

2、机(PSO-SVR)的高浓磨浆系统游离度的软测量方法。本发明提供一种高浓磨浆系统游离度的软测量方法，该方法能够实现对高浓磨浆机出浆游离度的实时测量，提高人工测量的效率和精度，为工业实时控制和实现节能降耗提供技术支持。本发明方法由硬件平台及测量软件组成，其中硬件平台核心包括高浓磨浆系统三盘磨浆机、压力传感器、震动传感器、加速度传感器、流量传感器、分料螺旋等。 (51)Int.Cl. (19)中华人民共和国国家知识产权局 (12)发明专利申请权利要求书4页说明书10页附图3页 (10)申请公布号 CN 104459089 A (43)申请公布日 2015.0。

3、3.25 CN 104459089 A 1/4 页 2 1. 一种高浓磨浆系统游离度的软测量方法，包括以下步骤： (1) 辅助变量选择和模型输入变量确定， (2)PSO-SVR 软测量模型的训练和软测量实现。 (1) 辅助变量选择和模型输入变量确定需要软测量的高浓磨浆系统纸浆关键参数为游离度(ml)。根据生产工艺及机理以及可测可观和变量间相关性分析 ( 因子分析、协方差分析 )，选择 i 个软测量辅助变量 (i 9)，包括：设定产量 u1(adt/d) ；磨浆机功率 u2(Mw) ；传动侧稀释水量 u3(l/min) ；非传动侧稀释水量 u4(l/min) ；非传。

4、动侧震动加速度 u5( ) ；非传动侧震动 u6(mm/s) ；非传动侧磨盘间隙 u7(mm) ；传动侧磨盘压力 u8(bar) ；非传动侧磨盘压力 u9(bar)。确定上述9个软测量辅助变量为软测量模型的输入变量，游离度(ml)为软测量模型的输出变量。 (2)PSO-SVR 软测量模型的训练和软测量实现 (A) 开始：变量初始化； (B) 选择进行软测量模型训练，转至 (C)，读取软测量训练样本；选择高浓磨浆系统游离度软测量，转至 (K)，读取软测量模型最优参数 ( 惩罚因子 C 和核函数参数 )； (C) 读取软测量训练样本：从数据。

5、库导入或输入训练集 T (x1,y1),L,(xl,yl) (Rn)l。其中，软测量模型输入： xi Rn， i 1,2,L,l，软测量模型输出： yi R， i 1,2,L,l， l 为训练集样本数量； (D) 数据预处理：采用尖峰滤波算法将高浓磨浆系统的尖峰异常数据进行剔除； 9 个辅助变量存在着量纲差，对每个变量分别进行中心标准化处理，即每类样本数据减去对应样本的平均值，然后除以其样本方差，即：计算均值：计算方差：标准化：权利要求书 CN 104459089 A 2 2/4 页 3 其中， xij为第 i 类变量的第。

6、j 个训练样本， l 为训练样本数，为第 i 类变量的均值，为中心标准化后的训练样本，表示第 i 类变量的方差，表示第 i 类变量的标准差； (E)确定-SVR软测量模型的初始化参数，包括：惩罚因子C和误差容限值，核函数 (g) 的类型及其核函数参数。惩罚因子 C 取 20，误差容限值取 0.01，核函数取高斯径向基函数 K(x,x) exp(-|x-x|2/2)，其中核函数参数 1 ； (F) 确定 -SVR 待优化参数惩罚因子 C 的搜索范围为 0,100，核函数参数的搜索范围为 0,50，粒子群个体的数目 40，粒子的最大迭代寻优次数 600，粒子的。

7、惩罚因子 C 迭代速度的最大值为 25，粒子的核函数参数迭代速度的最大值为 25，粒子的初始位置、初始速度、个体极值以及全局极值随机生成； (G)-SVR 模型训练： -SVR 计算高浓磨浆系统游离度的回归超平面方程 y (gx)+b，其中为其法方向， b 为其截距。设模型样本的第j个目标输出为yj，则-SVR将拟合问题等价转化为如下所述的二次规划问题：其中， C 为惩罚因子，为误差容限值， i和( 简记为) 为松弛变量。求解式 (4) 中二次规划问题的解 (,b,) 后，得到回归超平面方程： y (gx)+b (5) 将 Lagrange 函数引入上。

8、述问题，导出其对偶问题：解得 Lagrange 乘子向量计算 b ：选择位于开区间中的 (*)的分量 j或若选到的是 j，则有：若选择的是则有：权利要求书 CN 104459089 A 3 3/4 页 4 构造回归方程：其中， y为训练样本i的输出。 (xigxk)是核函数(g)，这里取高斯径向基函数K(x,x) exp(-|x-x|2/2)，是对应的拉格朗日乘子的各个分量。 (H) 计算适应度：设定优化的目标函数，并计算适应度，对每组惩罚因子 C 和核函数参数的效果进行评价，通过相应的误差函数计算适应度。对于误差大的粒子，认为其适应度小，。

9、粒子的适应度 fp表示为： fp 1/(ep+1) (10) 其中， ep是回归方程的误差函数，即：其中，是回归方程的预测输出， yi是回归方程的目标输出；更新每个粒子的速度和位置其中，表示待更新粒子的速度，表示待更新粒子的位置， pid表示待更新粒子的个体最优值， pgd表示整个粒子群的全局最优值。w 为惯性权重， c1、 c2为加速度常数， rand1 和 rand2为 0,1 范围内的随机数；对粒子 p，如果更新后的适应度值大于原个体的最优值，则更新粒子的个体最优值 pid，即： pid fp (13) 如果粒子 p 的个体最优值 pid大于原来粒子群的全。

10、局最优值 pgd，则更新原来的全局最优值 pgd： pgd pid (14) 。 2. 根据权利要求 1 所述一种高浓磨浆系统游离度的软测量方法，其特征在于还包括以下步骤： (I) 建模效果评价：判定建模效果是否满足要求。若满足，结束寻优，得到惩罚因子 C 和核函数参数的最优值，进入步骤 (J)，保存 -SVR 惩罚因子 C 和核函数参数的最优参数；否则返回步骤 (F)，继续迭代寻优，直到达到最大寻优代数，结束寻优； (J) 保存 -SVR 惩罚因子 C 和核函数参数的最优参数； (K) 读取 -SVR 惩罚因子 C 和核函数参数的最优参数；。

11、权利要求书 CN 104459089 A 4 4/4 页 5 (L) 读取软测量模型的输入数据； (M)将输入变量数据中心标准化处理后，使用-SVR惩罚因子C和核函数参数的最优参数进行高浓磨浆系统游离度在线软测量； (N) 游离度软测量结果显示：在软测量系统人机界面上显示游离度软测量的结果； (O) 游离度软测量结果保存：将本次软测量的输入输出数据保存到相应的数据库，供查询、修正使用； (P) 判断软测量是否结束，结束转至步骤 (Q)，否则转至步骤 (L) 继续进行游离度软测量； (Q) 结束：完成游离度软测量。 3. 根据权利要求 2 所述一种高。

12、浓磨浆系统游离度的软测量方法，其特征在于还包括以下步骤：定期将游离度离线间歇式数据输入到训练集中，更新游离度 - 支持向量回归机软测量模型。权利要求书 CN 104459089 A 5 1/10 页 6 一种高浓磨浆系统游离度的软测量方法技术领域 0001 本发明涉及化学机械制浆生产过程高浓磨浆系统游离度的测量技术领域，具体地说是涉及一种基于粒子群算法优化 - 支持向量回归机 (PSO-SVR) 的高浓磨浆系统游离度的软测量方法。背景技术 0002 造纸过程由制浆和造纸两大环节串联而成。制浆环节的主要功能是从植物纤维原材料生产出具有特定形态的纤维，而造纸环。

13、节的功能主要是以特定形态的纤维为原料生产出各种纸制产品。制浆造纸工业在我国是重点能耗行业之一，生产过程中消耗大量热能、电能和水能。特别是制浆过程的能耗及其所产生的纸浆的质量直接关系到后续造纸环节的能耗及产品质量。 0003 按现有工艺，制浆分为机械浆、化学浆和化学机械浆三大类。其中，化学机械浆生产过程主要包括：木片预处理，高浓磨浆，粗浆漂白，粗浆洗涤，低浓磨浆，渣浆磨浆，纤维回收，浆料浓缩脱水等过程。高浓磨浆作为化学机械浆生产的关键步骤，其用电量占到整个工厂用电的 35以上，决定着纸浆质量和生产能耗。 0004 然而，高浓磨浆过程具有多变量、。

14、强耦合和非线性的特征，使得高浓磨浆过程的机理分析、建模存在很大困难。目前，研究的假设性强，得到的高浓磨浆模型缺乏通用性。此外，研究主要集中在低浓磨浆过程、单盘磨浆机和磨盘改进上，对于高浓磨浆系统三盘磨浆机的建模和分析还没有公开报道。 0005 纤维网滤水性 ( 游离度 ) 是检测纸浆通过铜网的滤水程度快慢的一个指标。作为评价纸浆质量的关键指标之一，该指标不仅关系着最终纸张的质量 ( 透光性，吸水性等 )，同时与整个生产过程的能耗密切相关。为此，使用游离度来评价纤维质量十分合理。 0006 公开资料表明，目前还没有传感器对于高浓磨浆系统游离度直接进行测量，主要。

15、依靠离线间歇人工测量。该测量方法存在着较大的时滞性和误差，进而无法实时对制浆生产进行控制，阻碍实现节能降耗和纸浆质量控制。为此，提出一种基于粒子群算法优化 -支持向量机(PSO-SVR)的高浓磨浆系统游离度软测量的技术。本发明可通过常规在线测量方法提供的模型输入数据，给出输入数据对应的在线估计值，为制浆生产过程的优化操作和优化运行提供关键质量指标，进而实现质量控制和节能降耗。发明内容 0007 为了克服高浓磨浆过程中游离度无法直接测量，离线人工测量时效性差、精度低等缺点，本发明提供一种高浓磨浆系统游离度的软测量方法，该方法能够实现对高浓磨浆机出浆游离度的实时。

16、测量，提高人工测量的效率和精度，为工业实时控制和实现节能降耗提供技术支持。 0008 本发明方法由硬件平台及测量软件组成，其中硬件平台核心包括高浓磨浆系统三盘磨浆机、压力传感器、震动传感器、加速度传感器、流量传感器、分料螺旋等。说明书 CN 104459089 A 6 2/10 页 7 0009 一种最优参数的支持向量回归机高浓磨浆系统游离度软测量包括以下步骤： (1) 辅助变量选择和模型输入变量确定， (2)PSO-SVR 软测量模型的训练和软测量实现。 0010 (1) 辅助变量选择和模型输入变量确定 0011 需要软测量的高浓磨浆系统纸浆关键参数为游离度 (m。

17、l)。根据生产工艺及机理以及可测可观和变量间相关性分析 ( 因子分析、协方差分析 )，选择 i 个软测量辅助变量 (i 9)，包括： 0012 设定产量 u1(adt/d) ； 0013 磨浆机功率 u2(Mw) ； 0014 传动侧稀释水量 u3(l/min) ； 0015 非传动侧稀释水量 u4(l/min) ； 0016 非传动侧震动加速度 u5( ) ； 0017 非传动侧震动 u6(mm/s) ； 0018 非传动侧磨盘间隙 u7(mm) ； 0019 传动侧磨盘压力 u8(bar) ； 0020 非传动侧磨盘压力 u9(bar)。 0021 确定上述9个软测量辅助变量为软。

18、测量模型的输入变量，游离度(ml)为软测量模型的输出变量。 0022 (2)PSO-SVR 软测量模型的训练和软测量实现 0023 (A) 开始：变量初始化； 0024 (B) 选择进行软测量模型训练，转至 (C)，读取软测量训练样本；选择高浓磨浆系统游离度软测量，转至 (K)，读取软测量模型最优参数 ( 惩罚因子 C 和核函数参数 ) ； 0025 (C) 读取软测量训练样本：从数据库导入或输入训练集 T (x1,y1),L,(xl,yl) (Rn)l。其中，软测量模型输入： xi Rn， i 1,2,L,l，软测量模型。

19、输出： yi R， i 1,2,L,l， l 为训练集样本数量； 0026 (D) 数据预处理： 0027 采用尖峰滤波算法将高浓磨浆系统的尖峰异常数据进行剔除；由于选取的 9 个辅助变量存在着量纲差，对每个变量分别进行中心标准化处理，即每类样本数据减去对应样本的平均值，然后除以其样本方差，即： 0028 计算均值： 0029 0030 计算方差： 0031 0032 标准化： 0033 说明书 CN 104459089 A 7 3/10 页 8 0034 其中， xij为第 i 类变量的第 j 个训练样本， l 为训练样本数，为第 i 类变量的均值，为。

20、中心标准化后的训练样本，表示第 i 类变量的方差，表示第 i 类变量的标准差； 0035 (E)确定-SVR软测量模型的初始化参数，包括：惩罚因子C和误差容限值，核函数 (g) 的类型及其核函数参数。惩罚因子 C 取 20，误差容限值取 0.01，核函数取高斯径向基函数 K(x,x ) exp(-|x-x|2/2)，其中核函数参数 1 ； 0036 (F) 确定 -SVR 待优化参数惩罚因子 C 的搜索范围为 0,100，核函数参数的搜索范围为 0,50，粒子群个体的数目 40，粒子的最大迭代寻优次数 600，粒子的惩罚因子 C迭代速度的最大值为25，粒。

21、子的核函数参数迭代速度的最大值为25，粒子的初始位置、初始速度、个体极值以及全局极值随机生成； 0037 (G)-SVR 模型训练： 0038 -SVR 计算高浓磨浆系统游离度的回归超平面方程 y (gx)+b，其中为其法方向， b 为其截距。 0039 设模型样本的第i个目标输出为yi，则-SVR将拟合问题等价转化为如下所述的二次规划问题： 0040 0041 其中， C 为惩罚因子，为误差容限值， i和( 简记为) 为松弛变量。 0042 求解式 (4) 中二次规划问题的解后，得到回归超平面方程： 0043 y (gx)+b (5) 0044 将 Lagrang。

22、e 函数引入上述问题，导出其对偶问题： 0045 0046 解得 Lagrange 乘子向量 0047 计算 b ：选择位于开区间 (0,C) 中的 (*)的分量 j或若选到的是 j，则有： 0048 说明书 CN 104459089 A 8 4/10 页 9 0049 若选择的是则有： 0050 0051 构造回归方程： 0052 0053 其中， y为训练样本输出。 (xigxk)是核函数(g)，这里取高斯径向基函数K(x,x) exp(-|x-x|2/2)，是对应的拉格朗日乘子的各个分量。 0054 (H) 计算适应度： 0055 设定优化的目标函数，并计算适应度。

23、，对每组惩罚因子 C 和核函数参数的效果进行评价，通过相应的误差函数计算适应度。对于误差大的粒子，认为其适应度小，粒子的适应度 fp表示为： 0056 fp 1/(ep+1) (10) 0057 其中， ep是回归方程的误差函数，即： 0058 0059 其中，是回归方程的预测输出， yi是回归方程的目标输出； 0060 更新每个粒子的速度和位置 0061 0062 其中，表示待更新粒子的速度，表示待更新粒子的位置， pid表示待更新粒子的个体最优值， pgd表示整个粒子群的全局最优值。w 为惯性权重， c1、 c2为加速度常数， rand1和 rand2为 0,1 。

24、范围内的随机数； 0063 对粒子 p，如果更新后的适应度值大于原个体的最优值，则更新粒子的个体最优值 pid，即： 0064 pid fp (13) 0065 如果粒子 p 的个体最优值 pid大于原来粒子群的全局最优值 pgd，则更新原来的全局最优值 pgd： 0066 pgd pid (14) 0067 作为一种优选方案：本发明还包括以下步骤： 0068 (I) 建模效果评价： 0069 判定建模效果是否满足要求。若满足，结束寻优，得到惩罚因子C和核函数参数的最优值，进入步骤 (J)，保存 -SVR 惩罚因子 C 和核函数参数的最优参数；否则返回。

25、步骤 (F)，继续迭代寻优，直到达到最大寻优代数，结束寻优； 0070 (J) 保存 -SVR 惩罚因子 C 和核函数参数的最优参数； 0071 (K) 读取 -SVR 惩罚因子 C 和核函数参数的最优参数；说明书 CN 104459089 A 9 5/10 页 10 0072 (L) 读取软测量模型的输入数据； 0073 (M) 将输入变量数据中心标准化处理后，使用 -SVR 惩罚因子 C 和核函数参数的最优参数进行高浓磨浆系统游离度在线软测量； 0074 (N) 游离度软测量结果显示：在软测量系统人机界面上显示游离度软测量的结果； 0075 (O) 游。

26、离度软测量结果保存：将本次软测量的输入输出数据保存到相应的数据库，供查询、修正使用； 0076 (P) 判断软测量是否结束，结束转至步骤 (Q)，否则转至步骤 (L) 继续进行游离度软测量； 0077 (Q) 结束：完成游离度软测量。 0078 作为另一种优选方案：本发明还包括以下步骤： 0079 定期将游离度离线间歇式数据输入到训练集中，更新游离度 - 支持向量回归机软测量模型。 0080 本发明的一种高浓磨浆系统的游离度软测量方法具有以下技术效果： 0081 本发明针对目前高浓磨浆过程中无法直接测量游离度的现状，在没有设备测量游离度的前提下，选取对。

27、高浓磨浆系统游离度有直接影响的 9 个操作变量和检测变量作为软测量模型的输入数据，实现对高浓磨浆系统游离度的软测量，其测量精度高，为实现整个造纸过程的节能降耗提供理论依据和模型支持。附图说明： 0082 图1是本发明选择的研究对象：高浓磨浆系统三盘磨浆机的结构及13个主要操作变量和检测变量 0083 图 2 是本发明的一种基于粒子群优化 - 支持向量回归机 (PSO-SVR) 的高浓磨浆系统游离度软测量算法的流程图 0084 图 3 是本发明选择的高浓磨浆系统制浆生产过程 0085 高浓磨浆过程中，经预处理后的木片经过传动侧分料螺旋和非传动侧分料螺旋，传动侧负载感应螺。

28、旋和非传动侧负载感应螺旋进入高浓磨浆机，经过高浓磨机磨浆后的浆料进入旋风分离器从而实现气浆分离，便形成了粗浆。 0086 图 4 是本发明寻得最优参数对游离度进行软测量的效果图具体实施方式 0087 如图 1 所示，高浓磨浆系统三盘磨浆机由中间一个具有双齿面结构的定盘与两侧两个定盘构成双磨室结构。三盘磨浆机在高转速时运行平稳，不会出现动盘偏斜等问题，在不提高转速及增大磨盘直径的情况下，磨盘面积增加一倍，具有提高产量、改进磨浆质量、方便热能回收等优点。 0088 本发明提供一种高浓磨浆系统三盘磨浆机游离度的软测量方法。该方法能够实现对高浓磨浆机出浆游离度的实时测量。

29、，提高人工测量的效率和精度，为工业实时控制和实现节能降耗提供技术支持。 0089 参照图 1、图 2，一种基于粒子群算法优化 - 支持向量回归机 (PSO-SVR) 的高说明书 CN 104459089 A 10 6/10 页 11 浓磨浆系统游离度软测量的方法，其特征在于有常规测量系统、数据采集器、 PSO-SVR 软测量软件以及运行软测量软件的计算机系统构成。流量计、压力机、震动计、加速度计等常规测量仪表安装于高浓磨浆系统的各个相应位置。数据采集器连接常规测量系统，并通过通信总线连接运行软测量软件的计算机系统。常规测量系统主要包括如下常规测量仪表： 00。

30、90 二个流量计 ( 春辉 LW JY-LWGY)，分别用于在线测量高浓磨浆系统传动侧稀释水量，非传动侧稀释水量； 0091 一个振动加速度传感器 (YMK 151A100)，用于在线测量高浓磨浆系统非传动侧的震动； 0092 一个加速度计 (PCB SA51SC)，用于在线测量高浓磨浆系统非传动侧震动加速度的大小； 0093 一个电涡流位移传感器 ( 航振 HZ891XL)，用于在线测量高浓磨浆系统非传动侧定盘与动盘之间的间隙； 0094 二个压力计 (GP50211)，用于在线测量高浓磨浆系统传动侧磨盘压力和非传动侧磨盘压力。 0095 此外，常规测量系统还。

31、包括如下两个测量仪： 0096 一个高浓磨浆系统设定产量测量仪：通过传动侧分料螺旋和非传动侧分料螺旋转速设定来设定高浓磨浆系统产量； 0097 一个高浓磨浆系统磨浆机功率测量仪：通过各子系统能耗计量得到高浓磨浆系统磨浆机功率。 0098 本发明实现方法包括： (1) 辅助变量选择和模型输入变量确定； (2)PSO-SVR 软测量模型的训练和软测量实现。 0099 (1) 辅助变量选择和模型输入变量确定 0100 需要软测量的高浓磨浆系统纸浆关键参数为游离度(ml)。根据生产工艺及机理以及可测可观和变量间相关性分析 ( 因子分析、协方差分析 )，选择 i 个 (i 。

32、9) 软测量辅助变量包括： 0101 设定产量 u1(adt/d) ； 0102 磨浆机功率 u2(Mw) ； 0103 传动侧稀释水量 u3(l/min) ； 0104 非传动侧稀释水量 u4(l/min) ； 0105 非传动侧震动加速度 u5( ) ； 0106 非传动侧震动 u6(mm/s) ； 0107 非传动侧磨盘间隙 u7(mm) ； 0108 传动侧磨盘压力 u8(bar) ； 0109 非传动侧磨盘压力 u9(bar)。 0110 确定上述9个软测量辅助变量为为软测量模型的输入变量，游离度(ml)作为软测量模型的输出变量。 0111 (2)PSO-SVR 软测量模型的。

33、训练和软测量实现 0112 (A) 开始：变量初始化； 0113 (B) 选择进行软测量模型训练，转至 (C)，读取软测量训练样本；选择高浓磨浆系说明书 CN 104459089 A 11 7/10 页 12 统游离度软测量，转至 (K)，读取软测量模型最优参数 ( 惩罚因子 C 和核函数参数 ) ； 0114 (C) 读取软测量训练样本：从数据库导入或输入训练集 T (x1,y1),L,(xl,yl) (Rn)l。其中，软测量模型输入： xi Rn， i 1,2,L,l，软测量模型输出： yi R， i 1,2,L,l。

34、， l 为训练集样本数量； 0115 (D) 数据预处理： 0116 采用尖峰滤波算法将高浓磨浆系统的尖峰异常数据进行剔除；由于选取的 9 个辅助变量存在着量纲差，对每个变量分别进行中心标准化处理，即每类样本数据减去对应样本的平均值，然后除以其样本方差，即： 0117 计算均值： 0118 0119 计算方差： 0120 0121 标准化： 0122 0123 其中， xij为第 i 类变量的第 j 个训练样本， l 为训练样本数，为第 i 类变量的均值，为中心标准化后的训练样本，表示第 i 类变量的方差，表示第 i 类变量的标准差； 0124 (E)确定。

35、-SVR软测量模型的初始化参数，包括：惩罚因子C和误差容限值，核函数 (g) 的类型及其核函数参数。惩罚因子 C 取 20，误差容限值取 0.01，核函数取高斯径向基函数 K(x,x) exp(-|x-x|2/2)，其中核函数参数 1 ； 0125 (F) 确定 -SVR 待优化参数惩罚因子 C 的搜索范围为 0,100，核函数参数的搜索范围为 0,50，粒子群个体的数目 40，粒子的最大迭代寻优次数 600，粒子的惩罚因子 C迭代速度的最大值为25，粒子的核函数参数迭代速度的最大值为25，粒子的初始位置、初始速度、个体极值以及全局极值随机生成； 0。

36、126 (G)-SVR 模型训练： 0127 -SVR 计算高浓磨浆系统游离度的回归超平面方程 y (gx)+b，其中为其法方向， b 为其截距。 0128 设模型样本的第i个目标输出为yi，则-SVR将拟合问题等价转化为如下所述的二次规划问题：说明书 CN 104459089 A 12 8/10 页 13 0129 0130 其中， C 为惩罚因子，为误差容限值， i和( 简记为) 为松弛变量。 0131 求解上述二次规划问题的解后，得到回归方程： 0132 y (gx)+b (19) 0133 将 Lagrange 函数引入上述问题，导出其对偶问题： 0134。

37、 0135 解得 Lagrange 乘子向量 0136 计算 b ：选择位于开区间 (0,C) 中的 (*)的分量 j或若选到的是 j，则有： 0137 0138 若选择的是则有： 0139 0140 构造回归方程： 0141 0142 其中， y为训练样本输出。 (xigxk)是核函数(g)，这里取高斯径向基函数K(x,x) exp(-|x-x|2/2)，是对应的拉格朗日乘子的各个分量。 0143 (H) 计算适应度： 0144 设定优化的目标函数，并计算适应度，对每组惩罚因子 C 和核函数参数的效果进行评价，通过相应的误差函数计算适应度。对于误差大的粒子，认为其适。

38、应度小，粒子的适应度 fp表示为： 0145 fp 1/(ep+1) (24) 0146 其中， ep是回归方程的误差函数，即：说明书 CN 104459089 A 13 9/10 页 14 0147 0148 其中，是回归方程的预测输出， yi是回归方程的目标输出； 0149 更新每个粒子的速度和位置 0150 0151 其中，表示待更新粒子的速度，表示待更新粒子的位置， pid表示待更新粒子的个体最优值， pgd表示整个粒子群的全局最优值。w 为惯性权重， c1为加速度常数， rand1 和 rand2为 0,1 范围内的随机数； 0152 对粒子 p，如果更新后。

39、的适应度值大于原个体的最优值，则更新粒子的个体最优值 pid，即： 0153 pid fp (27) 0154 如果粒子 p 的个体最优值 pid大于原来粒子群的全局最优值 pgd，则更新原来的全局最优值 pgd： 0155 pgd pid (28) 0156 (I) 建模效果评价： 0157 判定建模效果是否满足要求。若满足，结束寻优，得到惩罚因子C和核函数参数的最优值，进入步骤 (J)，保存 -SVR 惩罚因子 C 和核函数参数的最优参数；否则返回步骤 (F)，继续迭代寻优，直到达到最大寻优代数，结束寻优； 0158 (J) 保存 -SVR 惩罚因子。

40、 C 和核函数参数的最优参数； 0159 (K) 读取 -SVR 惩罚因子 C 和核函数参数的最优参数； 0160 (L) 读取软测量模型的输入数据； 0161 (M) 将输入变量数据中心标准化处理后，使用 -SVR 惩罚因子 C 和核函数参数的最优参数进行高浓磨浆系统游离度在线软测量； 0162 (N) 游离度软测量结果显示：在软测量系统人机界面上显示游离度软测量的结果； 0163 (O) 游离度软测量结果保存：将本次软测量的输入输出数据保存到相应的数据库，供查询、修正使用； 0164 (P) 判断软测量是否结束，结束转至步骤 (Q)，否则转至步骤 (。

41、L) 继续进行游离度软测量； 0165 (Q) 结束：完成游离度软测量。 0166 作为优选的一种方案：所述软测量方法还包括以下步骤： 0167 定期将游离度离线间歇式数据输入到训练集中，更新游离度 - 支持向量回归机软测量模型。 0168 图 4 为一段时间内，高浓磨浆系统游离度软测量的效果，可以看出游离度的预测值与其实际值基本一致，误差比较小，且变化趋势基本一致。此外，本发明方法速度快、精度高、泛化能力强且有严格的数学解释，相比于其他方法有较高的优越性。因此本发明是一种说明书 CN 104459089 A 14 10/10 页 15 具有很高。

42、实用价值的、低成本的高浓磨浆系统游离度计量手段。 0169 可以理解的是，以上关于本发明的具体描述，仅用于说明本发明而并非受限于本发明实施例所描述的技术方案，本领域的普通技术人员应当理解，仍然可以对本发明进行修改或等同替换，以达到相同的技术效果；只要满足使用需要，都在本发明的保护范围之内。说明书 CN 104459089 A 15 1/3 页 16 图 1 说明书附图 CN 104459089 A 16 2/3 页 17 图 2 图 3 说明书附图 CN 104459089 A 17 3/3 页 18 图 4 说明书附图 CN 104459089 A 18 。

展开阅读全文