1.造纸机干燥部关键工艺参数的一种在线监测方法，其特征在于包括如下步骤：(1)纸张在干燥部干燥过程中遵守物料和能量守恒定律，以此为基础建立以干燥部监测部位处的纸张温度(Tp)和湿度(u)、监测部位处袋区空气的温度(Ta)和湿度(AHa)为未知量的四元偏微分方程组过程模型；

(2)实时在线采集造纸机车速(V)、成纸定量(BWo)、成纸水份(Wo)，实时在线采集纸张在监测部位进口处的湿度(u0)和温度(Tp，0)，实时在线采集监测部位处通风与余热回收子系统的送风风量(ma，sup)、送风温度(Ta，sup)及送风气体湿度(AHa，atm)数据，实时在线采集监测部位处蒸汽冷凝水子系统的蒸汽压力(Ps)数据；

(3)利用步骤(2)中采集的实时在线监测数据，以纸张监测部位进口的湿度(u0)和温度(Tp，0)、监测部位处通风与余热回收子系统的送风温度(Ta，sup)及送风气体湿度(AHa，atm)为初值，以监测部位处蒸汽冷凝水子系统的蒸汽压力(Ps)为边界条件，采用Runge‑Kutta法数值求解步骤(1)中的四元偏微分方程组，即可实时获得监测点处的纸张温度(Tp)和湿度(μ)及监测点处袋区空气的温度(Ta)和湿度(AHa)数据。

2.一种如权利要求1所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的四元偏微分方程组过程模型如下：

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10；

Δ1和Δ2为与干燥面有关的系数，若监测点处的干燥面由烘缸‑纸张‑干网‑空气组成则Δ1＝FRF、Δ2＝1，若监测点处的干燥面由烘缸‑干网‑纸张‑空气或烘缸‑纸张‑空气组成则Δ1＝1、Δ2＝1，若监测点处的干燥面由空气‑纸张‑干网‑空气组成则Δ1＝1+FRF、Δ2＝2，若监测点处的干燥面由空气‑纸张‑空气组成则Δ1＝2、Δ2＝2，FRF表示干网影响因子；Rp和Rs为与供热有关的系数，当纸张直接接受来自蒸汽提供的热量时Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，否则Rp为0，当袋区空气直接接受来自蒸汽提供的热量时Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)，否则Rs为0；

u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。

3.一种如权利要求2所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的监测点为单挂上排供热烘缸由烘缸‑纸张‑干网‑空气组成的干燥界面Phase1，所建立的四元偏微分方程组为

所述的监测点为单挂上排供热烘缸由空气‑纸张‑干网‑空气组成的干燥界面Phase2，所建立的四元偏微分方程组为

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10 ，Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)；

FRF表示干网影响因子，u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。

4.一种如权利要求2所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的监测点为单挂上排不供热烘缸由烘缸‑纸张‑干网‑空气组成的干燥界面Phase1，所建立的四元偏微分方程组为

所述的监测点为单挂上排供不热烘缸由空气‑纸张‑干网‑空气组成的干燥界面Phase2，所建立的四元偏微分方程组为ma，sup(Ca+CvAHa)]

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10 ，Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)；

FRF表示干网影响因子，u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。

5.一种如权利要求2所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的监测点为单挂下排供热烘缸由烘缸‑干网‑纸张‑空气组成的干燥界面Phase1，所建立的四元偏微分方程组为

所述的监测点为单挂下排供热烘缸由空气‑纸张‑干网‑空气组成的干燥界面Phase2，所建立的四元偏微分方程组为

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10 ，Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)；

FRF表示干网影响因子，u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。

6.一种如权利要求2所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的监测点为单挂下排不供热烘缸由烘缸‑干网‑纸张‑空气组成的干燥界面Phase1，所建立的四元偏微分方程组为

所述的监测点为单挂下排不供热烘缸由空气‑纸张‑干网‑空气组成的干燥界面Phase2，所建立的四元偏微分方程组为

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10 ，Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)；

FRF表示干网影响因子，u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。

7.一种如权利要求2所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的监测点为双挂供热烘缸由烘缸‑纸张‑空气组成的干燥界面Phase1和phase3，所建立的四元偏微分方程组为

ma，sup(Ca+CvAHa)]

所述的监测点为双挂供热烘缸由烘缸‑纸张‑干网‑空气组成的干燥界面Phase2，所建立的四元偏微分方程组为

所述的监测点为双挂供热烘缸由空气‑纸张‑空气组成的干燥界面Phase4，所建立的四元偏微分方程组为

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10 ，Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)；

FRF表示干网影响因子，u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。

8.一种如权利要求2所述的造纸机烘缸关键工艺参数的在线监测方法，其特征在于，所述的监测点为双挂不供热烘缸由烘缸‑纸张‑空气组成的干燥界面Phase1和phase3，所建立的四元偏微分方程组为

所述的监测点为双挂不供热烘缸由烘缸‑纸张‑干网‑空气组成的干燥界面Phase2，所建立的四元偏微分方程组为

所述的监测点为双挂不供热烘缸由空气‑纸张‑空气组成的干燥界面Phase4，所建立的四元偏微分方程组为

1.877

其中，mev，a＝{exp[23.1934‑3819.44/(Tp+227.03)]·[1‑exp(‑47.58u ‑

1.0585

0.10085Tpu )]‑exp[23.1934‑3819.44/(Ta+227.03)]}/(Ta+273.15)，ΔHev＝

1.0585 2 1.877 1.0585

46.544292·u .(Tp+273.15) ·{1/[1‑exp(‑47.58u ‑0.10085Tpu )]‑1}+3

(2504.7‑2.4789Tp)·10 ，Rp＝hs‑p(3816.44/(23.1934‑lnPs)‑227.03‑Tp)，Rs＝hs‑a(3816.44/(23.1934‑lnPs)‑227.03‑Ta)；

FRF表示干网影响因子，u表示纸张湿度，Tp表示纸张温度，AHa表示袋区空气绝对湿度，Ta表示袋区空气温度，t表示时间；

Cp，dp表示绝干纸张的比热取定值1423.5J/(kg·℃)，Ca表示绝干空气的比热取定值

1010J/(kg·℃)，Cv表示水蒸气的比热取定值1880J/(kg·℃)，γ0表示在标准大气压下0℃2

时水的汽化潜热取定值2501000J/kg，Cp，w＝0.0139(Tp/℃) ‑1.3129(Tp/℃)+4206.1表示与温度有关的水的比热，ε表示通风机给对应的烘缸组送风量的分配系数，Ly表示纸机幅宽，hs‑p表示蒸汽与纸张之间的对流传热系数，hs‑a表示蒸汽与袋区空气之间的对流传热系数，hp‑a表示纸张与袋区空气之间的对流传热系数，Rv是气体常数取定值461.5J/(kg·℃)，Kp‑a表示纸张‑空气接触面对流传质系数。