专利名称:一种预测多相混合均匀性的时间序列模型和方法
技术领域:
本发明涉及属于化学工程技术领域,具体地说是一种适应于化工领域所有流体混 合均勻性预测的方法。
背景技术:
时间序列分析是根据系统观测得到的时间序列数据,通过曲线拟合和参数估计来 建立数学模型的理论和方法。它一般采用曲线拟合和参数估计方法(如非线性最小二乘 法)进行。时间序列分析常用在国民经济宏观控制、区域综合发展规划、企业经营管理、市 场潜量预测、气象预报、水文预报、地震前兆预报、农作物病虫灾害预报、环境污染控制、生 态平衡、天文学和海洋学等方面。时间序列建模基本步骤是①用观测、调查、统计、抽样等 方法取得被观测系统时间序列动态数据。②根据动态数据作相关图,进行相关分析,求自相 关函数。相关图能显示出变化的趋势和周期,并能发现跳点和拐点。跳点是指与其他数据 不一致的观测值。如果跳点是正确的观测值,在建模时应考虑进去,如果是反常现象,则应 把跳点调整到期望值。拐点则是指时间序列从上升趋势突然变为下降趋势的点。如果存在 拐点,则在建模时必须用不同的模型去分段拟合该时间序列,例如采用门限回归模型。③辨 识合适的随机模型,进行曲线拟合,即用通用随机模型去拟合时间序列的观测数据。对于短 的或简单的时间序列,可用趋势模型和季节模型加上误差来进行拟合。对于平稳时间序列, 可用通用ARMA模型(自回归滑动平均模型)及其特殊情况的自回归模型、滑动平均模型或 组合-ARMA模型等来进行拟合。当观测值多于50个时一般都采用ARMA模型。对于非平稳 时间序列则要先将观测到的时间序列进行差分运算,化为平稳时间序列,再用适当模型去 拟合这个差分序列时间序列分析主要用于①系统描述。根据对系统进行观测得到的时间 序列数据,用曲线拟合方法对系统进行客观的描述。②系统分析。当观测值取自两个以上 变量时,可用一个时间序列中的变化去说明另一个时间序列中的变化,从而深入了解给定 时间序列产生的机理。③预测未来。一般用ARMA模型拟合时间序列,预测该时间序列未来 值。④决策和控制。根据时间序列模型可调整输入变量使系统发展过程保持在目标值上, 即预测到过程要偏离目标时便可进行必要的控制。DPS数据处理系统提供给用户一套较完 整的时间序列建模分析、进行预测预报的工具,包括平稳无趋势时间序列分析预测、有趋势 的时间序列预测、具季节性周期的时间序列预测以及差分自回归滑动平均(ARIMA)建模分 析、预测等时间序列分析和建模技术。然而上述的方法都是针对非混沌的时间序列,对混沌 的则预测模型复杂且效果差。特别是在化工领域混合均勻性方面则更差。而对整个过程实 验验证则费时费力,而且低成本,还不能对非正常操作进行检测和避免,甚至对事故产生的 避免。
发明内容
本发明的目的是克服了上述方法的不足,提供一种具有较高的应用价值的、简便 可行的,适应性强,预测效果好的用于化学工业中流体混合均勻性的预测模型和方法。
本发明预测多相混合均勻性的时间序列模型的建立方的技术方案是该方法采用混沌理论以及混沌时间序列预测方法,具体步骤如下(1)利用电子断层成像技术或高速摄像机来获得多相搅拌混合实时图样;(2)用编写好的程序计算所获得的实时图样的第0维贝蒂数;(3)用(2)的方法每隔一定时间获得之初的相关数据,从而得到相应时间序列;(4)利用(3)得到的时间序列应用相重构法计算时间延迟变量和嵌入维数以及带 宽;(5)再对(3)中的时间序列采用一步向前法和局部加权线性预测,从而就有完整的时间 序列;(6)再用实验方法来验证我们的预测误差。针对透明或半透明流体该方法通过高速摄像机来获得混合的实时图样;针对不透 明流体采用电子断层成像技术EPT来获得混合的的实时图样。实验中针对均相流体可以加入示踪粒子。之前每隔一定时间t截取搅拌反应器横截面或者纵截面来获得混合图样,贝蒂数 的计算可以用chomp免费软件获得,从而获得相应时间序列。时间延迟变量的选取采用最大互信息法,嵌入维数则采用伪最小邻近法,带宽则 采用0. 05的置信水平来选取。采用高斯核函数作为权重,公式如下K(Guass) = e(_(而 1獻仙)2)距离用By correlation来衡量,如下 采用均方根百分比误差来衡量预测效果,且该方法的预测误差控制在5%以内,公 式如下在算代数拓扑学中,第0维贝蒂数的含义是区域中连通成分的个数,简单的说就 是区域中块的个数。借助此理论知识对流体混合的复杂相图进行计算,定量给出该混合区 域中块的个数,表征了混合均勻性。即在统计学中非参数时间序列预测一步向前局部加权 线性法是用相重构的思想进行的。本发明利用该理论来预测多相混合均勻性。具体方法 是(1)利用电子断层成像技术或高速摄像机来获得多相搅拌混合实时图样;(2)用编写好 的程序计算所获得的实时图样的第0维贝蒂数;(3)用(2)的方法每隔一定时间获得之初 的相关数据,从而得到相应时间序列;(4)利用(3)得到的时间序列应用相重构法计算时间 延迟变量和嵌入维数以及带宽;(5)再对(3)中的时间序列采用一步向前法和局部加权线 性预测,从而就有完整的时间序列(6)再用实验方法来验证我们的预测误差。若流体为均勻相,即分辨不出来各相的界面,则可以加入相应的示踪粒子。本发明的有益效果是1该方法补充了经典的时间序列方法如ARMA模型,曲线拟合和参数估计方法(如
M_非线性最小二乘法)的不足;2、该方法简单可行,它是一种非参数预测方法;3、可以对搅拌反应器里的不正常操作进行检测并避免危险的产生,且省时省钱。本发明应用在对所有流体混合均勻性的预测,该方法简单方便,且具有很高的实 用价值,对化工实验中预测混合均勻性及理论上指导搅拌反应器的设计及搅拌的操作,提 供了一种可靠实用的预测方法,且可有效及时的避免事故的发生。
图1为本发明的第0维贝蒂数时间序列实际跟预测对比图。图2为本发明的第0维贝蒂数时间序列实际跟预测对比图。图3为本发明的第0维贝蒂数时间序列实际跟预测对比图。
具体实施例方式具体步骤如下(1)利用电子断层成像技术或高速摄像机来获得多相搅拌混合实 时图样;(2)用编写好的程序计算所获得的实时图样的第0维贝蒂数;(3)用(2)的方法每 隔一定时间获得之初的相关数据,从而得到相应时间序列;(4)利用(3)得到的时间序列应 用相重构法计算时间延迟变量和嵌入维数以及带宽;(5)再对(3)中的时间序列采用一步 向前法和局部加权线性预测,从而就有完整的时间序列(6)再用实验方法来验证我们的预 测误差。针对透明或半透明流体该方法通过高速摄像机来获得混合的实时图样;针对不透 明流体采用电子断层成像技术FPT来获得混合的的实时图样。实验中针对均相流体可以加 入示踪粒子。之前每隔一定时间t截取搅拌反应器横截面或者纵截面来获得混合图样,贝 蒂数的计算可以用chomp免费软件获得,从而获得相应时间序列。时间延迟变量的选取采 用最大互信息法,嵌入维数则采用伪最小邻近法,带宽则采用0. 05的置信水平来选取。采 用高斯核函数作为权重,公式如下K(Guass) = e{-(disX3nce'hf)距离用By correlation来衡量,如下D(By correlation)=丨1 ,。采用均方根百分比误差来衡量预测效果,且
该方法的预测误差控制在5 %以内,公式如下
1 I"RMSPE =- Tiyl - y* f 100% 。
mean \实施例1 某造纸厂用一种以木薯根碎纸碎木屑混合物为原材料制作纸浆的方法,一开始就 使用混合原材料,其中木薯根的成分多少都可以,但比重在5-50%比较合适,其他的物质为 碎纸及碎木屑,比重为95-50%,以上比例是重量比。现在实验室条件下进行混合均勻判断 去适量原料放在搅拌器中加水混合。
5
通过电子断层成像技术 EPT (Electrical Process Tomography)或 CT(Computerized Tomography)(针对不透明流体)来获得混合的实时图样,保存数据, 然后借助计算机计算贝蒂数的程序进行计算,将第0维贝蒂数的时间序列计算出来,测量 1-9000,预测 9001-9100 (如图 1),经计算 RMSPE = 4. 89% 实施例2 在一种工业制备催化剂方法中将铜、锌、铝的可溶性盐类共溶于去离子水,将模板 剂在另一个容器中溶于去离子水,然后将这两种水溶液在强搅拌下充分混合,生成一种均 勻的混合溶液,现在实验室条件下进行混合实验通过粒子测速仪(半透明流体)、来获得混 合的实时图样,保存数据,然后借助计算机计算贝蒂数的程序进行计算,将第0维贝蒂数的 时间序列计算出来,测量1-9000,预测9001-9100。(如图2),经计算RMSPE = 4. 51%0实施例3 化学实验中配置某试剂,将NaCl、(NH4) 2S04、、NH4C1固体放到磁力搅拌器中加水 进行混合,由于放在水中属于均相,故需加入示踪粒子来监控,其混合状态。通过粒子测速 仪(针对透明流体)来获得混合的实时图样保存数据,然后借助计算机计算贝蒂数的程序 进行计算,将第0维贝蒂数的时间序列计算出来,测量1-9000,预测9001-9100。(如图3), 经计算 RMSPE = 5. 00%。
权利要求
一种预测多相混合均匀性的时间序列模型的建立方法,其特征在于该方法采用混沌理论以及混沌时间序列预测方法,具体步骤如下(1)利用电子断层成像技术或高速摄像机来获得多相搅拌混合实时图样;(2)用编写好的程序计算所获得的实时图样的第0维贝蒂数;(3)用(2)的方法每隔一定时间获得之初的相关数据,从而得到相应时间序列;(4)利用(3)得到的时间序列应用相重构法计算时间延迟变量和嵌入维数以及带宽;(5)再对(3)中的时间序列采用一步向前法和局部加权线性预测,从而就有完整的时间序列;(6)再用实验方法来验证我们的预测误差。
2.根据权利要求1所述的预测多相混合均勻性的时间序列模型的建立方法,其特征在 于针对透明或半透明流体该方法通过高速摄像机来获得混合的实时图样;针对不透明流 体采用电子断层成像技术EPT来获得混合的的实时图样。
3.根据权利要求1所述预测多相混合均勻性的时间序列模型的建立方法,其特征在 于实验中针对均相流体可以加入示踪粒子。
4.根据权利要求1所述的预测多相混合均勻性的时间序列模型的建立方法,其特征在 于之前每隔一定时间t截取搅拌反应器横截面或者纵截面来获得混合图样,贝蒂数的计 算可以用chomp免费软件获得,从而获得相应时间序列。
5.根据权利要求1所述的预测多相混合均勻性的时间序列模型的建立方法,其特征 在于时间延迟变量的选取采用最大互信息法,嵌入维数则采用伪最小邻近法,带宽则采用 0. 05的置信水平来选取。
6.根据权利要求1所述的预测多相混合均勻性的时间序列模型的建立方法,其特征在 于采用高斯核函数作为权重,公式如下K{Guass) = 彻—)2)距离用By correlation来衡量,如下
7.根据权利要求1所述的预测多相混合均勻性的时间序列模型和方法,其特征在于 采用均方根百分比误差来衡量预测效果,且该方法的预测误差控制在5%以内,公式如下RMSPE =So 厂父)2100%
全文摘要
本发明是一种预测多相混合均匀性的时间序列模型的建立方法。本发明主要应用于化学工程实验中预测流体混合效果及理论上指导及校正实验的设计。具体是(1)利用电子断层成像技术或高速摄像机来获得多相搅拌混合实时图样;(2)用编写好的程序计算所获得的实时图样的第0维贝蒂数;(3)用(2)的方法每隔一定时间获得之初的相关数据,从而得到相应时间序列;(4)利用(3)得到的时间序列应用相重构法计算时间延迟变量和嵌入维数以及带宽;(5)再对(3)中的时间序列采用一步向前法和局部加权线性预测,从而就有完整的时间序列;(6)再用实验方法来验证我们的预测误差。本发明应用在对所有流体混合效果的预测,该方法简单方便,且具有很高的实用价值。
文档编号G01N33/00GK101929994SQ201010261128
公开日2010年12月29日 申请日期2010年8月24日 优先权日2010年8月24日
发明者徐建新, 朱道飞, 王仕博, 王 华, 范国锋 申请人:昆明理工大学