CN109446236A - 基于随机分布的水泥粒径分布预测方法 - Google Patents
基于随机分布的水泥粒径分布预测方法 Download PDFInfo
- Publication number
- CN109446236A CN109446236A CN201811218359.0A CN201811218359A CN109446236A CN 109446236 A CN109446236 A CN 109446236A CN 201811218359 A CN201811218359 A CN 201811218359A CN 109446236 A CN109446236 A CN 109446236A
- Authority
- CN
- China
- Prior art keywords
- cement
- particle size
- distribution
- function
- weight
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002245 particle Substances 0.000 title claims abstract description 75
- 238000009826 distribution Methods 0.000 title claims abstract description 73
- 238000013277 forecasting method Methods 0.000 title claims abstract description 10
- 239000004568 cement Substances 0.000 claims abstract description 43
- 238000000227 grinding Methods 0.000 claims abstract description 42
- 239000000843 powder Substances 0.000 claims abstract description 13
- 230000002159 abnormal effect Effects 0.000 claims abstract description 4
- 238000012216 screening Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 29
- 238000012549 training Methods 0.000 claims description 17
- 238000012360 testing method Methods 0.000 claims description 10
- 238000005315 distribution function Methods 0.000 claims description 8
- 239000008187 granular material Substances 0.000 claims description 8
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 2
- 238000011897 real-time detection Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000003062 neural network model Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010298 pulverizing process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Dispersion Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Crushing And Grinding (AREA)
- Disintegrating Or Milling (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于水泥细度预测研究技术领域,具体是一种基于随机分布的水泥粒径分布预测方法。包括以下步骤,S100:获取水泥粉磨***中磨机喂料量、磨机工作电流、磨机出入口温度、磨机压差、循环提升机电流、选粉机转速、循环风机电流以及对应时刻水泥粒径概率分布,并将所有参数信号存储为历史数据集;S200:建立水泥粒径概率分布密度函数的基函数表示模型;S300:筛除异常数据,依照分类赋值样本权重,构成新数据样本;S400:建立输入变量与前n‑1个权值向量之间的非线性预测模型,预测下一时刻水泥粒径分布;S500:通过模型输出误差值,更新水泥粒径概率分布密度函数的基函数表示模型参数。本发明可实时检测水泥粉磨***中的水泥粒径分布。
Description
技术领域
本发明属于水泥细度预测研究技术领域,具体是一种基于随机分布的水泥粒径分布预测方法。
背景技术
水泥是一种粉体产品,由不同大小的颗粒按照一定的比例构成。水泥粉磨***是将水泥熟料研磨成符合国家质量标准的水泥产品的过程。它是生产水泥过程的最后一个环节,也是至关重要的环节。目前,水泥质量标准的评价以水泥细度为主,水泥细度对水泥强度、凝结时间、准稠度用水量等有着很大的影响。水泥细度主要包含筛余、比表面积、颗粒级配三项指标,而在企业生产中通常采用筛余或比表面积来进行水泥细度的测定和评判。实际上,筛余和比表面积只能部分体现水泥的颗粒大小情况,并不能完全表现水泥真实的质量信息。
水泥细度检测分为在线和离线检测,在线检测的主要设备是粒度监测仪,尽管可以快速、实时地对水泥细度进行测量,但是该类设备价格昂贵且维修成本高,并不适合大多数中小规模企业。而离线监测通过每一个小时到现场采样,将样本在实验室进行分析,但是由于离线检测取样量小,导致其检测结果的代表性不够。另外,间隔时间长会导致最后测定的水泥比表面积和细度跟实际生产不一致。软测量技术的出现很好解决了以上问题,目前针对水泥粒径细度的软测量技术主要通过神经网络、模糊控制、最小二乘法等建立预测模型。这些软测量技术均以水泥细度的百分比为预测结果,不能完全反应水泥实际质量性能。另外,也未考虑水泥粒径分布具有非高斯特性,导致预测精度不高。
发明内容
本发明为了解决上述技术问题,提供一种基于随机分布的水泥粒径分布预测方法。
本发明采取以下技术方案:一种基于随机分布的水泥粒径分布预测方法,包括以下步骤,
S100:获取水泥粉磨***中磨机喂料量、磨机工作电流、磨机出入口温度、磨机压差、循环提升机电流、选粉机转速、循环风机电流以及对应时刻水泥粒径概率分布,并将所有参数信号存储为历史数据集;
S200:建立水泥粒径概率分布密度函数的基函数表示模型;
S300:筛除异常数据,依照分类赋值样本权重,构成新数据样本;
S400:建立输入变量与前n-1个权值向量之间的非线性预测模型,预测下一时刻水泥粒径分布;
S500:通过模型输出误差值,更新水泥粒径概率分布密度函数的基函数表示模型参数。
所述S100的具体方法如下,历史数据集包括磨机喂料量Wmw、磨机电流Imd、磨机入口温度Tmrw、磨机出口温度Tmcw、磨机压差Pmy、循环提升机电流Imw、选粉机转速Vxf和循环风机电流Ixf,数据集为P组数据,而同一组下包含K时刻的各个数据,具体表示为Ui(K)=[Wi mw(k),Ii md(k),Ti mrw(k),Ti mcw(k),Pi my(k),Ii xt(k),Vi xz(k),Ii xf(k)]K*8i=1,2,3.....P,利用实验室离线分析获取相应时刻的水泥粒径概率分布函数Pi k(D,Ui),i=1,2,3...P,其中D为采集的水泥颗粒直径信息,Ui为输入参数。
所述S200具体方法如下,
S201:从S100获得的水泥粒径概率分布Pi K(D,Ui)利用公式1求出水泥粒径分布概率分布密度函数
其中d1、d2为水泥颗粒直径;
S202:建立水泥粒径概率分布函数的基函数表达式,采用高斯形式的RBF网络作为基函数,其表达式为:
其中D为采集的水泥颗粒直径信息,μi,σi为第i个网络节点函数的中心值和宽度;
根据RBF网络逼近原理,确定基函数的权值表示,将水泥粒径分布的概率密度函数利用基函数加权和的形式进行表示,其表达式为:
其中,C(D)=[B1(D),B2(D),…,Bn-1(D)],Bn(D)为第n个基函数,V(k)=[w1(k),w2(k),…,wn-1(k)]T,wn(k)为第n个基函数对应的权值,e0(D,k)为对不同情况、不同时刻的概率分布密度函数逼近的误差;
S203:确定每个基函数的权值,第n个权值ωn(k)可用权值向量V(k)的非线性函数h(V(k))表示为:
其中,
忽略逼近误差,结合式3和4可以得到:
两边左乘[CT(D)Bn(D)]T,并在区间[Dmin Dmax]上进行积分,当
矩阵非奇异时,可以转化得到:
利用上式6求出水泥粒径概率分布密度函数各个基函数的权值
V(k)。
所述S300具体方法如下,
S301:将磨机喂料量Wmw、磨机电流Imd、磨机入口温度Tmrw、
磨机出口温度Tmcw、磨机压差Pmy、循环提升机电流Imw、选粉机转
速Vxf、循环风机电流Ixf、前一时刻变量U(K-1)和水泥粒径分布概率
密度权值V(K)组合成m个训练样本X;
X=[Wi mw(k),Ii md(k),Ti mrw(k),Ti mcw(k),Pi my(k),Ii xt(k),Vi xz(k),Ii xf(k),Ui(K-1),Vi(k)]K*(8+8+m)i=1,2,3,...m,m∈P
S302:初始给定各均值向量μj以及分类数G,计算样本Xi与μj之
间的距离;
根据距离最近的均值向量确定xi的族标记:λi=argmini∈{1,2,3...k}dji,将样本xi划入相应的族:
S303:计算新的均值向量
S304:判断是否μ′j=μj,重复S302至S303,直至满足本条件;
S305:依据对训练样本分类结果,按照式子进行对相同分类的样本赋权重值;
nk:同类的样本个数
NP:训练样本总数
由此得到带有权重的训练样本X′=WX。
所述S400具体方法如下,
S401:选择随机权神经网络,其网络模型表示为:
其中,ωj=[ωj1,ωj2,…,ωjm]T为m个输入节点连接第j个隐含单元的输入权重,βj=[βj1,βj2,…,βj(n-1)]T为第j个隐含层连接输出节点的输出权重,bj是第j个隐含单元的偏置,对模型进行训练,随机给定一组输入层权值和偏置,利用采集的n组历史数据对模型进行训练,其目标函数为:
其中Xi为第i个训练样本,Vi为Xi对应的粒径分布权重。
通过求取H矩阵的广义逆得到最优的输出权值其公式为:
S402:建立基函数权值与输入变量之间的关系:利用测试集Xp-m检测该模型:
S403:输出测试集数据水泥粒径概率分布gp-m通过式12可得:
gp-m(D,k+1)=C(D)Vp-m(k+1)+Bn(T)wn(k+1) (12)
其中gp-m为预测的水泥粒径分布概率函数。
所述S500具体方法如下,
S501:定义每一组Xi的性能指标函数Zn:
Zn是关于基函数参数的函数,其中K为每组总采样时间个数,γp-m(D,U(k))为测试集Xp-m的水泥粒径实际输出概率分布密度函数,通过式14和式15更新基函数参数μ和σ;
其中,μn=[μ1 μ2 ... μn],σn=[σ1 σ2 ... σn],n为定义的基函数参数个数;
S502:通过式16可变换为概率分布函数,由此可预测下一时刻的水泥粒径分布:
与现有技术相比,本发明采用随机分布算法预测水泥粒径的概率分布密度函数,克服了水泥粒径分布具有随机性,避免了筛余、比表面积、颗粒级配的缺点,完全具体地反应出被检测水泥的质量信息。本发明通过采用k-means分类算法对数据分类、赋值权重,以及采用性能指标函数对基函数的参数优化,提高了该模型的可靠性、准确度。同时,本发明可实时检测水泥粉磨***中的水泥粒径分布。
附图说明
图1为***工作框图;
图2为***实施实例。
实施方式
基于随机分布预测算法的水泥粉磨质量监测在实施方式上,可以采用建模过程和实时过程来实现:
建模过程包括以下步骤:
S100:将所有参数信号存储为历史数据集,建立历史数据集可以采用如下方式实现:
历史数据需要采集磨机喂料量、磨机电流、磨机入口温度、磨机出口温度、磨机压差、循环提升机电流、选粉机转速、循环风机电流。在送料前,把物料放到称重台进行称重获取磨机喂料量,固定选取P组喂料量。对该磨机、循环提升机、循环风机设备安装电流互感器记录设备工作电流,在磨机3/1入口处和3/1出口处安装温度传感器记录温度。安装转速测量设备在选粉机上以获取选粉机转速。对获取的数据先后进行采样、过滤、平滑处理,获取K个时间点的标准数据集。利用实验室离线分析获取此时对应的水泥粒径分布。
S200:建立水泥粒径概率分布密度函数的基函数表示模型;通过基函数对已知的水泥粒径概率分布进行表示,具体采用如下方式实现:
S201:从S100获得的水泥粒径概率分布Pi K(D,Ui)利用公式1求出水泥粒径分布概率分布密度函数
其中d1、d2为水泥颗粒直径;
S202:建立水泥粒径概率分布函数的基函数表达式,采用高斯形式的RBF网络作为基函数,其表达式为:
其中D为采集的水泥颗粒直径信息,μi,σi为第i个网络节点函数的中心值和宽度;
根据RBF网络逼近原理,确定基函数的权值表示,将水泥粒径分布的概率密度函数利用基函数加权和的形式进行表示,其表达式为:
其中,C(D)=[B1(D),B2(D),…,Bn-1(D)],Bn(D)为第n个基函数,V(k)=[w1(k),w2(k),…,wn-1(k)]T,wn(k)为第n个基函数对应的权值,e0(D,k)为对不同情况、不同时刻的概率分布密度函数逼近的误差;
S203:确定每个基函数的权值,第n个权值ωn(k)可用权值向量V(k)的非线性函数h(V(k))表示为:
其中,
忽略逼近误差,结合式3和4可以得到:
两边左乘[CT(D) Bn(D)]T,并在区间[Dmin Dmax]上进行积分,当
矩阵非奇异时,可以转化得到:
利用上式6求出水泥粒径概率分布密度函数各个基函数的权值V(k)。
S300:筛除异常数据,依照分类赋值权重,构成新训练数据样本,具体采用如下实现方式:
S301:将磨机喂料量Wmw、磨机电流Imd、磨机入口温度Tmrw、磨机出口温度Tmcw、磨机压差Pmy、循环提升机电流Imw、选粉机转速Vxf、循环风机电流Ixf、前一时刻变量U(K-1)和水泥粒径分布概率密度权值V(K)组合成m个训练样本X;
X=[Wi mw(k),Ii md(k),Ti mrw(k),Ti mcw(k),Pi my(k),Ii xt(k),Vi xz(k),Ii xf(k),Ui(K-1),Vi(k)]K*(8+8+m)i=1,2,3,...m,m∈P
S302:初始给定各均值向量μj以及分类数G,计算样本Xi与μj之间的距离;
根据距离最近的均值向量确定xi的族标记:λi=argmini∈{1,2,3...k}dji,将样本xi划入相应的族:
S303:计算新的均值向量
S304:判断是否μj′=μj,重复S302至S303,直至满足本条件;
S305:依据对训练样本分类结果,按照式子进行对相同分类的样本赋权重值,本实施例选择通过K-means聚类算法,对数据样本进行分类,把样本划分为K个数据样本:
nk:同类的样本个数
NP:训练样本总数
由此得到带有权重的训练样本X′=WX。
S400:建立输入变量与前n-1个权值向量之间的非线性模型,将输入变量与水泥粒径概率分布密度建立联系,具体采用如下方式实现:
本实施例选择通过建立随机权值神经网络模型,建立输入变量与前n-1个权值之间的关系。首先从采集的历史数据中选择磨机喂料量Wmw、磨机电流Imd、磨机入口温度Tmrw、磨机出口温度Tmcw、磨机压差Pmy、循环提升机电流Imw、选粉机转速Vxf、循环风机电流Ixf,、水泥粒径概率分布P与此时前n-1个基函数权值向量共同组合构成随机权神经网络的输入变量,随机权神经网络模型形式为:
其中,ωj=[ωj1,ωj2,…,ωjm]T为m个输入节点连接第j个隐含单元的输入权重,βj=[βj1,βj2,…,βj(n-1)]T为第j个隐含层连接输出节点的输出权重,bj是第j个隐含单元的偏置,对模型进行训练,随机给定一组输入层权值和偏置,利用采集的n组历史数据对模型进行训练,其目标函数为:
其中Xi为第i个训练样本,Vi为Xi对应的粒径分布权重。
通过求取H矩阵的广义逆得到最优的输出权值其公式为:
S402:建立基函数权值与输入变量之间的关系:利用测试集Xp-m检测该模型:
S403:输出测试集数据水泥粒径概率分布gp-m通过式12可得:
gp-m(D,k+1)=C(D)Vp-m(k+1)+Bn(T)wn(k+1) (12)
其中gp-m为预测的水泥粒径分布概率函数。
利用采集的多组历史数据及S200的基函数表示方法,计算此时水泥概率分布密度函数的权重Vi,i=1,2,…,M,将此信息作为神经网络的输出,通过最小化网络目标函数得到最优的网络输出权值得到训练好的随机权神经网络模型。
S500:通过模型输出误差值,更新水泥粒径概率密度函数的基函数表示模型参数,具体采用如下实现方式:
S501:本实施例定义一个以基函数参数μ和σ为自变量的性能指标函数Zn其函数模型为:
该函数将S400建立的非线性预测模型输出的概率密度函数g(X)与实际值γ(X)做差求期望,更新的基函数参数为:
S502:通过式16可变换为概率分布函数,由此可预测下一时刻的水泥粒径分布:
利用更新后的基函数参数,优化整个模型,得到更好的水泥粒径概率分布检测模型。
通过上述5个步骤完成历史水泥粒径概率分布的基函数表示方法并确定输入变量与前n-1个权值向量的非线性模型,更新基函数参数,完成建模过程。
实时过程时主要包括以下内容:
利用S100所述方法实时采集水泥粉磨工艺参数信息并存储;
利用S400所述训练的随机权神经网络预测下一时刻的水泥粒径概率分布,其形式为:
γ(D,k+1)=C(D)V(k+1)+Bn(D)wn(k+1)。
本专利提出的基于随机分布预测算法的水泥粉磨颗粒预测整体流程如附图1所示,***实施实例如图2所示。本发明通过随机分布算法对水泥粉磨过程中的水泥粒径概率分布密度进行估计,建立粉磨工艺流程输入变量与水泥粒径概率分布之间的关系,实现了更具体清晰的表示水泥质量在线监测。本发明不仅提高在线监测的正确度,而且可以通过降低不必要的能源损耗为企业降低生产成本同时提高水泥质量性能。
通过上述具体实施方式可以实现基于随机分布预测算法的水泥粒径分布预测。但是本发明并不限于上述实施方式,在不脱离本发明精神和范围的任何修改或局部替换,均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种基于随机分布的水泥粒径分布预测方法,其特征在于:包括以下步骤,
S100:获取水泥粉磨***中磨机喂料量、磨机工作电流、磨机出入口温度、磨机压差、循环提升机电流、选粉机转速、循环风机电流以及对应时刻水泥粒径概率分布,并将所有参数信号存储为历史数据集;
S200:建立水泥粒径概率分布密度函数的基函数表示模型;
S300:筛除异常数据,依照分类赋值样本权重,构成新数据样本;
S400:建立输入变量与前n-1个权值向量之间的非线性预测模型,预测下一时刻水泥粒径分布;
S500:通过模型输出误差值,更新水泥粒径概率分布密度函数的基函数表示模型参数。
2.根据权利要求1所述的基于随机分布的水泥粒径分布预测方法,其特征在于:所述S100的具体方法如下,历史数据集包括磨机喂料量Wmw、磨机电流Imd、磨机入口温度Tmrw、磨机出口温度Tmcw、磨机压差Pmy、循环提升机电流Imw、选粉机转速Vxf和循环风机电流Ixf,数据集为P组数据,而同一组下包含K时刻的各个数据,具体表示为
Ui(K)=[Wi mw(k),Ii md(k),Ti mrw(k),Ti mcw(k),Pi my(k),Ii xt(k),Vi xz(k),Ii xf(k)]K*8i=1,2,3.....P,利用实验室离线分析获取相应时刻的水泥粒径概率分布函数Pi k(D,Ui),i=1,2,3...P,其中D为采集的水泥颗粒直径信息,Ui为输入参数。
3.根据权利要求2所述的基于随机分布的水泥粒径分布预测方法,其特征在于:所述S200具体方法如下,
S201:从S100获得的水泥粒径概率分布Pi K(D,Ui)利用公式1求出水泥粒径分布概率分布密度函数γi K(D,Ui);
其中d1、d2为水泥颗粒直径;
S202:建立水泥粒径概率分布函数的基函数表达式,采用高斯形式的RBF网络作为基函数,其表达式为:
其中D为采集的水泥颗粒直径信息,μi,σi为第i个网络节点函数的中心值和宽度;
根据RBF网络逼近原理,确定基函数的权值表示,将水泥粒径分布的概率密度函数利用基函数加权和的形式进行表示,其表达式为:
其中,C(D)=[B1(D),B2(D),…,Bn-1(D)],Bn(D)为第n个基函数,V(k)=[w1(k),w2(k),…,wn-1(k)]T,wn(k)为第n个基函数对应的权值,e0(D,k)为对不同情况、不同时刻的概率分布密度函数逼近的误差;
S203:确定每个基函数的权值,第n个权值ωn(k)可用权值向量V(k)的非线性函数h(V(k))表示为:
其中,
忽略逼近误差,结合式3和4可以得到:
两边左乘[CT(D) Bn(D)]T,并在区间[Dmin Dmax]上进行积分,当矩阵非奇异时,可以转化得到:
利用上式6求出水泥粒径概率分布密度函数各个基函数的权值V(k)。
4.根据权利要求3所述的基于随机分布的水泥粒径分布预测方法,其特征在于:所述S300具体方法如下,
S301:将磨机喂料量Wmw、磨机电流Imd、磨机入口温度Tmrw、磨机出口温度Tmcw、磨机压差Pmy、循环提升机电流Imw、选粉机转速Vxf、循环风机电流Ixf、前一时刻变量U(K-1)和水泥粒径分布概率密度权值V(K)组合成m个训练样本X;
X=[Wi mw(k),Ii md(k),Ti mrw(k),Ti mcw(k),Pi my(k),Ii xt(k),Vi xz(k),Ii xf(k),Ui(K-1),Vi(k)]K*(8+8+m)i=1,2,3,...m,m∈P
S302:初始给定各均值向量μj以及分类数G,计算样本Xi与μj之间的距离;
dij=||xi-μj||2 (7)
根据距离最近的均值向量确定xi的族标记:λi=argmini∈{1,2,3...k}dji,将样本xi划入相应的族:
S303:计算新的均值向量
S304:判断是否μ′j=μj,重复S302至S303,直至满足本条件;
S305:依据对训练样本分类结果,按照式子进行对相同分类的样本赋权重值;
nk:同类的样本个数;
NP:训练样本总数;
由此得到带有权重的训练样本X′=WX。
5.根据权利要求4所述的基于随机分布的水泥粒径分布预测方法,其特征在于:所述S400具体方法如下,
S401:选择随机权神经网络,其网络模型表示为:
其中,ωj=[ωj1,ωj2,…,ωjm]T为m个输入节点连接第j个隐含单元的输入权重,βj=[βj1,βj2,…,βj(n-1)]T为第j个隐含层连接输出节点的输出权重,bj是第j个隐含单元的偏置,对模型进行训练,随机给定一组输入层权值和偏置,利用采集的n组历史数据对模型进行训练,其目标函数为:
其中Xi为第i个训练样本,Vi为Xi对应的粒径分布权重;
通过求取H矩阵的广义逆得到最优的输出权值其公式为:
S402:建立基函数权值与输入变量之间的关系:利用测试集Xp-m检测该模型:
S403:输出测试集数据水泥粒径概率分布gp-m通过式12可得:
gp-m(D,k+1)=C(D)Vp-m(k+1)+Bn(T)wn(k+1) (12)
其中gp-m为预测的水泥粒径分布概率函数。
6.根据权利要求5所述的基于随机分布的水泥粒径分布预测方法,其特征在于:所述S500具体方法如下,
S501:定义每一组Xi的性能指标函数Zn:
Zn是关于基函数参数的函数,其中K为每组总采样时间个数,γp-m(D,U(k))为测试集Xp-m的水泥粒径实际输出概率分布密度函数,通过式14和式15更新基函数参数μ和σ;
其中,μn=[μ1 μ2...μn],σn=[σ1 σ2...σn],n为定义的基函数参数个数;
S502:通过式16可变换为概率分布函数,由此可预测下一时刻的水泥粒径分布:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811218359.0A CN109446236B (zh) | 2018-10-18 | 2018-10-18 | 基于随机分布的水泥粒径分布预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811218359.0A CN109446236B (zh) | 2018-10-18 | 2018-10-18 | 基于随机分布的水泥粒径分布预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109446236A true CN109446236A (zh) | 2019-03-08 |
CN109446236B CN109446236B (zh) | 2021-12-21 |
Family
ID=65547385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811218359.0A Active CN109446236B (zh) | 2018-10-18 | 2018-10-18 | 基于随机分布的水泥粒径分布预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109446236B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110222825A (zh) * | 2019-06-10 | 2019-09-10 | 燕山大学 | 一种水泥成品比表面积预测方法及*** |
CN110322014A (zh) * | 2019-07-10 | 2019-10-11 | 燕山大学 | 一种基于bp神经网络的水泥成品比表面积预测方法 |
CN110322077A (zh) * | 2019-07-10 | 2019-10-11 | 燕山大学 | 基于卷积神经网络的水泥原料立磨生料细度指标预测方法 |
CN110987737A (zh) * | 2019-12-23 | 2020-04-10 | 华中科技大学 | 一种基于光散射响应的气溶胶粒谱、浓度测量方法 |
CN112712861A (zh) * | 2021-01-07 | 2021-04-27 | 北京明略软件***有限公司 | 模型构建方法、装置、设备及计算机可读介质 |
CN112765868A (zh) * | 2020-12-29 | 2021-05-07 | 哈尔滨工业大学 | 一种早期水泥石毛细水中可冻水含量的预判方法 |
CN112949183A (zh) * | 2021-03-04 | 2021-06-11 | 上海交通大学宁波人工智能研究院 | 一种水泥生料立磨***异常工况检测***及方法 |
CN113405956A (zh) * | 2021-06-15 | 2021-09-17 | 中建材(合肥)粉体科技装备有限公司 | 一种粒度分析仪检测数据的在线校正方法、***、设备 |
CN114112819A (zh) * | 2022-01-25 | 2022-03-01 | 北矿智云科技(北京)有限公司 | 一种测量磨矿粒度的方法及装置 |
CN116046618A (zh) * | 2023-02-27 | 2023-05-02 | 一夫科技股份有限公司 | 一种超细硫酸钙的粉磨细度监测方法及*** |
CN117174219A (zh) * | 2023-11-02 | 2023-12-05 | 济南大学 | 基于隐变量模型的水泥强度条件概率分布估计方法及*** |
CN117174218A (zh) * | 2023-11-02 | 2023-12-05 | 济南大学 | 一种基于流模型的水泥强度条件概率分布估计方法及*** |
CN118060052A (zh) * | 2024-04-18 | 2024-05-24 | 山东埃尔派粉体科技股份有限公司 | 一种超细气流粉碎分级***及粒径调控方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7495436B2 (en) * | 2004-03-18 | 2009-02-24 | Baker Hughes Incorporated | Rock properties prediction, categorization, and recognition from NMR echo-trains using linear and nonlinear regression |
CN101482750A (zh) * | 2009-02-06 | 2009-07-15 | 北京矿冶研究总院 | 湿法冶金合成过程草酸钴粒度预测方法 |
CN104915505A (zh) * | 2015-06-18 | 2015-09-16 | 东北大学 | 一种高浓磨浆***输出纤维形态分布pdf建模方法 |
CN105020705A (zh) * | 2015-03-04 | 2015-11-04 | 内蒙古瑞特优化科技股份有限公司 | 循环流化床锅炉燃烧性能实时优化控制方法及*** |
CN106056243A (zh) * | 2016-05-27 | 2016-10-26 | 东北大学 | 一种高浓磨浆***输出纤维形态分布的控制***及方法 |
US20160370492A1 (en) * | 2015-02-26 | 2016-12-22 | Halliburton Energy Services, Inc. | Methods and systems employing nmr-based prediction of pore throat size distributions |
CN106651030A (zh) * | 2016-12-21 | 2017-05-10 | 重庆邮电大学 | 一种改进的rbf神经网络热点话题用户参与行为预测方法 |
CN106874568A (zh) * | 2017-01-17 | 2017-06-20 | 中南大学 | 一种球磨过程的物料粒径分布预测方法 |
-
2018
- 2018-10-18 CN CN201811218359.0A patent/CN109446236B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7495436B2 (en) * | 2004-03-18 | 2009-02-24 | Baker Hughes Incorporated | Rock properties prediction, categorization, and recognition from NMR echo-trains using linear and nonlinear regression |
CN101482750A (zh) * | 2009-02-06 | 2009-07-15 | 北京矿冶研究总院 | 湿法冶金合成过程草酸钴粒度预测方法 |
US20160370492A1 (en) * | 2015-02-26 | 2016-12-22 | Halliburton Energy Services, Inc. | Methods and systems employing nmr-based prediction of pore throat size distributions |
CN105020705A (zh) * | 2015-03-04 | 2015-11-04 | 内蒙古瑞特优化科技股份有限公司 | 循环流化床锅炉燃烧性能实时优化控制方法及*** |
CN104915505A (zh) * | 2015-06-18 | 2015-09-16 | 东北大学 | 一种高浓磨浆***输出纤维形态分布pdf建模方法 |
CN106056243A (zh) * | 2016-05-27 | 2016-10-26 | 东北大学 | 一种高浓磨浆***输出纤维形态分布的控制***及方法 |
CN106651030A (zh) * | 2016-12-21 | 2017-05-10 | 重庆邮电大学 | 一种改进的rbf神经网络热点话题用户参与行为预测方法 |
CN106874568A (zh) * | 2017-01-17 | 2017-06-20 | 中南大学 | 一种球磨过程的物料粒径分布预测方法 |
Non-Patent Citations (2)
Title |
---|
FANG MINGJIAN 等: "Study on air fine particles pollution prediction of main traffic route using artificial neural network", 《2011 INTERNATIONAL CONFERENCE ON COMPUTER DISTRIBUTED CONTROL AND INTELLIGENT ENVIRONMENTAL MONITORING》 * |
裘国华 等: "基于GM-RBF神经网络组合模型的水泥强度预测方法", 《科技导报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110222825A (zh) * | 2019-06-10 | 2019-09-10 | 燕山大学 | 一种水泥成品比表面积预测方法及*** |
CN110222825B (zh) * | 2019-06-10 | 2021-04-09 | 燕山大学 | 一种水泥成品比表面积预测方法及*** |
CN110322014A (zh) * | 2019-07-10 | 2019-10-11 | 燕山大学 | 一种基于bp神经网络的水泥成品比表面积预测方法 |
CN110322077A (zh) * | 2019-07-10 | 2019-10-11 | 燕山大学 | 基于卷积神经网络的水泥原料立磨生料细度指标预测方法 |
CN110322077B (zh) * | 2019-07-10 | 2022-08-02 | 燕山大学 | 基于卷积神经网络的水泥原料立磨生料细度指标预测方法 |
CN110987737A (zh) * | 2019-12-23 | 2020-04-10 | 华中科技大学 | 一种基于光散射响应的气溶胶粒谱、浓度测量方法 |
CN110987737B (zh) * | 2019-12-23 | 2021-02-05 | 华中科技大学 | 一种基于光散射响应的气溶胶粒谱、浓度测量方法 |
CN112765868B (zh) * | 2020-12-29 | 2022-11-25 | 哈尔滨工业大学 | 一种早期水泥石毛细水中可冻水含量的预判方法 |
CN112765868A (zh) * | 2020-12-29 | 2021-05-07 | 哈尔滨工业大学 | 一种早期水泥石毛细水中可冻水含量的预判方法 |
CN112712861A (zh) * | 2021-01-07 | 2021-04-27 | 北京明略软件***有限公司 | 模型构建方法、装置、设备及计算机可读介质 |
CN112949183A (zh) * | 2021-03-04 | 2021-06-11 | 上海交通大学宁波人工智能研究院 | 一种水泥生料立磨***异常工况检测***及方法 |
CN112949183B (zh) * | 2021-03-04 | 2024-01-16 | 上海交通大学宁波人工智能研究院 | 一种水泥生料立磨***异常工况检测***及方法 |
CN113405956A (zh) * | 2021-06-15 | 2021-09-17 | 中建材(合肥)粉体科技装备有限公司 | 一种粒度分析仪检测数据的在线校正方法、***、设备 |
CN113405956B (zh) * | 2021-06-15 | 2023-07-28 | 中建材(合肥)粉体科技装备有限公司 | 一种粒度分析仪检测数据的在线校正方法、***、设备 |
CN114112819A (zh) * | 2022-01-25 | 2022-03-01 | 北矿智云科技(北京)有限公司 | 一种测量磨矿粒度的方法及装置 |
CN116046618A (zh) * | 2023-02-27 | 2023-05-02 | 一夫科技股份有限公司 | 一种超细硫酸钙的粉磨细度监测方法及*** |
CN116046618B (zh) * | 2023-02-27 | 2024-02-27 | 一夫科技股份有限公司 | 一种超细硫酸钙的粉磨细度监测方法及*** |
CN117174219A (zh) * | 2023-11-02 | 2023-12-05 | 济南大学 | 基于隐变量模型的水泥强度条件概率分布估计方法及*** |
CN117174218A (zh) * | 2023-11-02 | 2023-12-05 | 济南大学 | 一种基于流模型的水泥强度条件概率分布估计方法及*** |
CN117174218B (zh) * | 2023-11-02 | 2024-01-26 | 济南大学 | 一种基于流模型的水泥强度条件概率分布估计方法及*** |
CN117174219B (zh) * | 2023-11-02 | 2024-02-02 | 济南大学 | 基于隐变量模型的水泥强度条件概率分布估计方法及*** |
CN118060052A (zh) * | 2024-04-18 | 2024-05-24 | 山东埃尔派粉体科技股份有限公司 | 一种超细气流粉碎分级***及粒径调控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109446236B (zh) | 2021-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109446236A (zh) | 基于随机分布的水泥粒径分布预测方法 | |
CN110458195B (zh) | 一种基于多判据融合的异常数据甄别方法 | |
CN104778837B (zh) | 一种道路交通运行态势多时间尺度预测方法 | |
Yao et al. | A support vector machine with the tabu search algorithm for freeway incident detection | |
CN108492555A (zh) | 一种城市路网交通状态评价方法及装置 | |
CN110322014A (zh) | 一种基于bp神经网络的水泥成品比表面积预测方法 | |
CN109858679A (zh) | 一种结合人机物的反窃电稽查监控***及其工作方法 | |
CN111369070A (zh) | 一种基于包络线聚类的多模融合光伏功率预测方法 | |
CN111079906A (zh) | 基于长短时记忆网络的水泥成品比表面积预测方法及*** | |
CN113762329A (zh) | 一种大型轧机状态预测模型的构建方法及构建*** | |
CN109675935A (zh) | 一种变控制限的ipca轧制过程在线故障诊断方法 | |
CN101863088A (zh) | 一种橡胶混炼过程中门尼粘度的预报方法 | |
Deng et al. | Bead geometry prediction for multi-layer and multi-bead wire and arc additive manufacturing based on XGBoost | |
CN113591374B (zh) | 基于卷积简单循环单元网络的水泥成品比表面积预测方法 | |
CN110210648A (zh) | 基于灰色长短期记忆网络的管制空域战略流量预测方法 | |
CN109342279B (zh) | 基于粉磨机理和神经网络的混合软测量方法 | |
CN113627735A (zh) | 工程建设项目安全风险的预警方法及*** | |
CN110245702A (zh) | 机械磨损颗粒识别方法、装置、电子设备及存储介质 | |
Lee et al. | An empirical experiment on deep learning models for predicting traffic data | |
CN110837857A (zh) | 工业用电负荷预测方法、***及其存储介质 | |
Lu et al. | Operating performance assessment based on semi-supervised cluster generative adversarial networks for gold flotation process | |
CN106874568A (zh) | 一种球磨过程的物料粒径分布预测方法 | |
CN109508820A (zh) | 基于差异化建模的校园用电量预测建模方法 | |
CN117251814A (zh) | 一种高速公路充电桩电量损耗异常的分析方法 | |
CN112541296A (zh) | 一种基于pso-lssvm的so2预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |