专利名称: |
一种颗粒系粒径分布的测量方法 |
摘要: |
本发明公开了一种颗粒系粒径分布的测量方法,包括以下步骤,利用近场散射技术,建立散射光强度分布与颗粒群粒径的积分方程关系式并离散化;构造泛函e,通过最小化泛函e来确定最佳奇异值截止值;采用L曲线法选取正则化参数α,结合最佳奇异值截止值构造滤波因子,得到粒径分布列向量的近似解形式;采用Chahine迭代算法进行迭代计算,得到收敛解,实现颗粒群粒径分布的测量。本方法引入的改进滤波因子极大程度地发挥模型的有效信息作用,对于粒径峰值间隔较近的颗粒系,可以实现高准确性和稳定性的粒径分布测量。 |
专利类型: |
发明专利 |
国家地区组织代码: |
江苏;32 |
申请人: |
东南大学 |
发明人: |
张彪;王晨;姚鸿熙;许传龙 |
专利状态: |
有效 |
申请日期: |
2019-03-04T00:00:00+0800 |
发布日期: |
2019-06-07T00:00:00+0800 |
申请号: |
CN201910162213.7 |
公开号: |
CN109856019A |
代理机构: |
南京苏高专利商标事务所(普通合伙) |
代理人: |
柏尚春 |
分类号: |
G01N15/02(2006.01);G;G01;G01N;G01N15 |
申请人地址: |
211100 江苏省南京市江宁区东南大学路2号 |
主权项: |
1.颗粒系粒径分布的测量方法,其特征在于,包括以下步骤: 步骤一:利用近场散射技术,将CCD相机拍摄到的颗粒系散斑图像处理成包含颗粒粒度信息的散射光强度分布,建立其与颗粒群粒径的积分方程关系式如下: IS(θ)=∫IMie(θ,D)N(D)dD (1) 式中,IS(θ)为散射角为θ时包含颗粒粒度信息的散射光强度,IMie(θ,D)是直径为D的颗粒在散射角为θ时的Mie散射光强,N(D)是颗粒粒径的频率分布; 将散射角θ等分为m份,粒径D等分为n组,离散化处理式(1)得到线性方程组并归一化成矩阵形式: I=AX 其中,I=(i1,i2,...,im)为散射光强角分布列向量,X=(x1,x2,…,xn)为粒径分布列向量,系数矩阵A的元素ai,j表示等效粒径为xi的颗粒在散射角θ时的散射光强; 步骤二:对系数矩阵A进行SVD分解: 其中,ui、vi是A的左右奇异值矩阵列向量,λi是递减排列的奇异值; 建立广义逆矩阵 其中,λk是最佳奇异值截止值; 步骤三:构造泛函e如式(2),通过最小化泛函e来确定最佳奇异值截止值λk; 式中,R为分辨率矩阵,R=A*A;E为单位阵;cov[X]为协方差矩阵,cov[X]=σ2A*(A*)T,σ为观测噪声的方差; 步骤四:采用L曲线法选取正则化参数α,结合最佳奇异值截止值λk构造滤波因子,得到粒径分布列向量的近似解形式: 步骤五:用SIRT算法对式(3)进行非负约束,然后作为Chahine迭代算法的初始值进行迭代计算,得到收敛解,实现颗粒群粒径分布的测量。 2.根据权利要求1所述的颗粒粒径分布测量方法,其特征在于:步骤四中,所述滤波因子如式(4): 其中,τi为滤波因子。 3.根据权利要求1所述的颗粒粒径分布测量方法,其特征在于:步骤一中,所述近场散射技术处理近场散斑图像的具体方法是:将采集到的散斑图像按灰度值矩阵存储,并归一化处理得到无量纲的近场散斑强度差分信号,接着通过快速傅里叶变化得到近场散斑强度的二维功率谱,将二维功率谱图像沿半径按等比分环后得到一维远场散射平均光强分布,即为包含颗粒粒度信息的散射光强度IS(θ)。 |
所属类别: |
发明专利 |