CN101869860A - 一种球磨机磨矿产品的粒度分布预测方法 - Google Patents
一种球磨机磨矿产品的粒度分布预测方法 Download PDFInfo
- Publication number
- CN101869860A CN101869860A CN201010191503A CN201010191503A CN101869860A CN 101869860 A CN101869860 A CN 101869860A CN 201010191503 A CN201010191503 A CN 201010191503A CN 201010191503 A CN201010191503 A CN 201010191503A CN 101869860 A CN101869860 A CN 101869860A
- Authority
- CN
- China
- Prior art keywords
- grade
- granularity
- broken
- percent
- product
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种球磨机磨矿产品的粒度分布预测方法,对每一个粒级破碎速率随时间变化快慢不一致的非线性破碎,将磨矿时间进行分段求取每个时间段上各粒级的破碎速率,对磨矿产品粒度进行预测时,按时间顺序先预测第一时间段的产品粒度分布,再以第一时间段的预测结果作为第二时间段的给料,预测第二时间段的产品粒度分布,依此类推预测出最终产品的粒度分布。本发明用于球磨机分批磨矿产品粒度分布预测,相对误差均在±5%内,绝对误差均在±2%内,减少工人的劳动强度,提高工作效率;扩展用于工业球磨机产品粒度预测,85.7%和72.6%的相对误差分别位于±10%和±5%内,91.2%和71.5%的绝对误差分别位于±5%和±3%内,本发明为实现工业磨矿过程的模拟和优化、节能降耗奠定了基础。
Description
技术领域本发明涉及一种球磨机磨矿产品的粒度分布预测方法。
背景技术磨矿历来是选矿过程中耗能最高的环节,其能耗约占整个选矿过程的50%,而且,随着矿物资源的贫乏,一些矿物品位也随之变化频繁,从而导致了其磨矿效率低、选矿过程工艺指标波动大,造成资源和能源的浪费。采用仿真技术对磨矿过程进行产品粒度预测和生产条件优化是一种有效的提高磨矿效率的手段,并已在国外得到了比较广泛的应用。
过程模型是实现过程仿真的基础。球磨机的模型描述了磨机输出粒度分布和输入粒度分布、磨矿条件以及矿物性质之间的关系。
磨矿过程中,某个粒级含量的累积=该粒级的破碎+该粒级的生成。则采用物料平衡及破碎动力学表示的时间连续、粒度离散的球磨机分批磨矿总体平衡模型为:
其中,t为时间;i和j为粒级编号,n为总粒级数,粒级1表示最粗的粒级,n表示最细的粒级;wi(t)为t时刻第i粒级的质量分数;S为破碎速率函数,粒度离散化为n个粒级后对应的S为对角矩阵则Si(t)表示t时刻第i粒级的破碎速率;b为破碎分布函数,粒度离散化为n个粒级后对应的b为下三角矩阵bij表示第j粒级破碎后进入第i粒级的量占第j粒级已破碎量的质量百分比。
如果已知破碎速率函数S和破碎分布函数b,即可采用公式(1)对磨矿产品的粒度进行预测。而bij和Si(t)可以根据实际的磨矿数据,直接计算得到,或者采用式(1)反算得到,然后再用于预测。
在实际的磨矿中,S可能仅与粒径有关而不依赖于时间,也可能既与粒径有关也随着时间的变化而变化。
当Si仅与粒径有关,而不依赖于时间时,则某粒级在磨矿产品中未破碎的量与其在给料中的量之比和时间的关系在半对数坐标上呈一条斜率为-S的直线(如附图1),称该粒级符合一阶破碎动力学。如果在磨碎过程中,第1粒级只有破碎没有生成,用公式表示其一阶破碎动力学为:
w1(t)=w1(0)exp(-S1t) (2)
则其破碎速率为:
研究发现,一些矿物在一定的磨矿条件下遵循一阶破碎的规律,此时可采用式(1)逐个粒级预测产品的粒度分布,也即首先计算第1粒级的粒度分布,然后将第1粒级的粒度分布代入式(1)中,计算第2粒级的粒度分布,从而逐步迭代实现对磨矿产品的粒度分布预测。
当物料的破碎速率随磨矿时间的增加而变化时,称为非一阶破碎(如附图1所示)。此时,式(1)除了粒级间相互关联外,破碎速率和粒度分布同时为时间的函数,难以求解。对非一阶破碎,目前的处理方法有:
1)将不同粒级矿物的破碎速率采用相同的依赖于时间的加速度表示,从而将实际的磨矿时间转化到另一个时间域(称作伪时间,false time)内,即令
Si(t)=k(t)Si(t=0)≡k(t)Si0 (4)
Si0为0时刻第i粒级的破碎速率,是一个常数,k(t)为所有粒级的破碎速率随时间变化的速度。
则所有粒级的破碎速率随着磨矿时间的增加而以相同的速度k(t)变化,令
θ即称为伪时间。
那么破碎速率随时间变化的非一阶破碎就可以在θ域上转化为一阶的(即破碎速率非时变)破碎模型:
进而可逐粒级迭代预测磨矿产品的粒度分布。很明显,该方法在不同粒级的破碎速率随时间变化的速度不同时,预测结果就会产生较大的误差。
2)采用一个依赖于粒径的显函数k(x)和一个依赖于粒度分布的泛函ζ的乘积来表示破碎速率:
S≡k(x)ζ[P(x,z)w(z)] (7)
其中,w(z)是粒径为z的颗粒的量,P(x,z)表示粒径为z的颗粒对粒径为x的颗粒的破碎速率的影响程度。此种方法可以描述不同颗粒之间的相互作用造成的非一阶破碎,特别是细颗粒对其他颗粒的破碎速率的影响,但该方法比较复杂,寻找一个合适的泛函是比较困难的。
发明内容
本发明的目的是提供一种球磨机磨矿产品的粒度分布预测方法。
本发明的球磨机磨矿产品粒度分布预测方法基于分段磨矿的思想,即对时间连续的磨矿,将磨矿时间进行分段,可视每一个时间段的磨矿为一次独立的磨矿,上一个时间段的磨矿产品为下一个时间段的给料,且每一时间段内的破碎速率可线性描述。基于此思想,对不同粒级的破碎速率具有不同加速度的非一阶破碎,将磨矿时间T分为0-t1,t1-t2,t2-t3,…,tr-2-tr-1,t>tr-1等r个段,根据已知的磨矿给料粒度分布和对应产品粒度分布,用0时刻的给料,求取0-t1时间段的破碎速率Si1,然后再以t1时刻的产品作为给料,求取t1-t2时间段的破碎速率Si2,依此类推,求取每一个时间段上各粒级的破碎速率,在每一时间段上各粒级的破碎速率为:
且在所有时间段上,i=n时,Si=0;
其中,wi(t)为第i粒级在t时刻的百分含量。
则在各时间段上第i粒级百分含量的预测公式为:
其中,θ表示时间,bij为破碎分布函数,表示第j粒级破碎后进入第i粒级的量占第j粒级已破碎量的质量百分比。
则由式(8)计算得到破碎速率函数S后,根据式(9)以及磨矿给料和对应产品粒度分布数据,可反算求得破碎分布函数b。
在对磨矿产品的粒度分布进行预测时,已知破碎分布函数b,由0时刻的给料,根据0-t1时间段的破碎速率Si1,预测得到t1时刻的粒度分布,然后再以t1时刻的粒度作为给料,根据t1-t2时间段的破碎速率Si2,预测t2时刻的粒度分布,依此类推,从而分步预测得到最终磨矿产品的粒度分布,而不是以零时刻的给料直接预测最终产品的粒度分布。
将该方法用于球磨机分批磨矿产品粒度分布的预测,预测结果的相对误差均位于±5%以内,绝对误差均位于±2%以内,从而可以代替人工试验,极大的减少了工作人员的劳动强度,提高了工作效率;将该方法扩展用于工业级球磨机连续磨矿产品粒度分布预测,预测精度85.7%的相对误差位于±10%以内,72.6%位于±5%以内,91.2%的绝对误差位于±5%内,71.5%位于±3%内,为实现工业磨矿过程的模拟和优化,达到节能降耗的目的奠定了基础。
该方法适用于任何呈现非一阶破碎的矿物和球磨机磨矿过程。
附图说明
图1为一阶破碎和非一阶破碎示意图,一阶破碎表现为在不同磨矿时间下,某粒级在磨矿产品中的量与在给料中的量之比与时间的关系,在半对数坐标上为一条直线,即斜率S(1)=S(1.7)=S(3),且而非一阶破碎,其破碎速率随着时间的增加而减小,假设0<t1<t2<T,其中,T为总的磨矿时间,则 且S(t1)>S(t2);
图2为分段线性化的非一阶破碎;
图3为球磨机分批磨矿产品粒度分布预测值与实测值对比;
图4为工业球磨机连续磨矿产品粒度分布预测值与实测值对比。
具体实施方式
下面给出本发明中球磨机磨矿产品粒度分布预测方法的具体实施过程。
实施例1
1.采用铜矿石,将球磨机分批磨矿给料和产品分为13个粒级,每一个粒级用其粒径的上限来表示该粒级的粒径,最细的粒级编号为13,最粗的粒级编号为1,则从最细到最粗,13个粒级的粒径为0.045,0.075,0.106,0.15,0.212,0.3,0.425,0.6,0.85,1.18,1.7,2.36,3.35,单位为mm;
2.将磨矿时间分为4个时间段,即0-1min,1-4min,4-8min,>8min,采用式(8)计算每一个粒级在每一时间段上的破碎速率如表1所示;
表1不同时间段上各粒级的破碎速率
粒级的粒径(mm) | 0-1min | 1-4min | 4-8min | >8min |
0.045 | 0 | 0 | 0 | 0 |
0.075 | 0.283 | 0.105 | 0.105 | 0.105 |
0.106 | 0.333 | 0.152 | 0.152 | 0.152 |
0.150 | 0.387 | 0.221 | 0.221 | 0.221 |
0.212 | 0.621 | 0.621 | 0.621 | 0.621 |
0.300 | 0.688 | 0.688 | 0.688 | 0.688 |
0.425 | 0.720 | 0.720 | 0.720 | 0.720 |
粒级的粒径(mm) | 0-1min | 1-4min | 4-8min | >8min |
0.600 | 0.870 | 0.870 | 0.870 | 0.870 |
0.850 | 1.205 | 1.205 | 1.205 | 1.205 |
1.180 | 1.130 | 1.063 | 1.063 | 1.063 |
1.700 | 0.870 | 0.601 | 0.601 | 0.601 |
2.360 | 0.654 | 0.324 | 0.308 | 0.308 |
3.350 | 0.468 | 0.183 | 0.140 | 0.109 |
3.由第2步所得到的S和式(9),采用L.G.Austin的反算法,从磨矿给料及对应产品的粒度分布数据反算求取破碎分布函数b;
4.预测过程为:当磨矿时间小于1min钟时,输入给料,由公式(9)和第一时间段各粒级的破碎速率,可以直接预测磨矿产品的粒度分布;当磨矿时间位于1-4min之间时,则先根据其给料和第一时间段的破碎速率,预测1min的产品粒度分布,然后再以预测得到的1min的粒度分布作为给料,根据第二时间段的破碎速率,预测最终的产品粒度分布;依此类推,从而可以预测4-8min和大于8min时的产品粒度分布。
用上述方法对分批磨矿0.5min,1min,2min,4min和8min的产品粒度分布进行预测,所得结果如图3所示,预测结果的相对误差均位于±5%以内,绝对误差均位于±2%以内;将该方法扩展用于工业级球磨机的磨矿产品粒度分布预测,共分为14个粒级,即0.045,0.075,0.106,0.15,0.212,0.3,0.425,0.6,0.85,1.18,1.7,2.36,3.35,5.6,单位为mm,预测精度相对误差85.7%位于±10%以内,72.6%位于±5%以内,绝对误差91.2%位于±5%内,71.5%位于±3%内,部分预测结果如图4所示。
Claims (1)
1.一种球磨机磨矿产品的粒度分布预测方法,其特征在于:基于球磨机磨矿的总体平衡模型
对不同粒级的破碎速率具有不同加速度的非一阶破碎,将磨矿时间T分为0-t1,t1-t2,t2-t3,…,tr-2-tr-1,t>tr-1共r个段,根据已知的磨矿给料和对应产品粒度分布,用0时刻的给料,求取0-t1时间段的破碎速率Si1,再以t1时刻的产品作为给料,求取t1-t2时间段的破碎速率Si2,依此类推,求取每一个时间段上各粒级的破碎速率,则在每一时间段上各粒级的破碎速率为:
且在所有时间段上,i=n时,Si=0;
则在各时间段上第i粒级百分含量的预测公式为:
其中,t和θ为均表示时间;i和j为粒级编号,n为总粒级数,粒级1表示最粗的粒级,n表示最细的粒级;wi(θ)为θ时刻第i粒级的质量分数;S为破碎速率函数,粒度离散化为n个粒级后对应的S为对角矩阵则Si(t)表示t时刻第i粒级的破碎速率;b为破碎分布函数,粒度离散化为n个粒级后对应的b为下三角矩阵bij表示第j粒级破碎后进入第i粒级的量占第j粒级已破碎量的质量百分比;
根据分段求得的破碎速率函数S,采用式(3)以及磨矿数据,反算求得破碎分布函数b,根据S和b,即可对磨矿产品的粒度分布进行预测;
预测时,已知b,由0时刻的给料,根据0-t1时间段的破碎速率Si1,采用式(3)预测得到t1时刻的粒度分布,然后再以t1时刻的粒度作为给料,根据t1-t2时间段的破碎速率Si2,预测t2时刻的粒度分布,依此类推,分步预测得到最终磨矿产品的粒度分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101915033A CN101869860B (zh) | 2010-06-04 | 2010-06-04 | 一种球磨机磨矿产品的粒度分布预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101915033A CN101869860B (zh) | 2010-06-04 | 2010-06-04 | 一种球磨机磨矿产品的粒度分布预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101869860A true CN101869860A (zh) | 2010-10-27 |
CN101869860B CN101869860B (zh) | 2012-05-23 |
Family
ID=42995056
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101915033A Expired - Fee Related CN101869860B (zh) | 2010-06-04 | 2010-06-04 | 一种球磨机磨矿产品的粒度分布预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101869860B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318330A (zh) * | 2014-10-20 | 2015-01-28 | 东北大学 | 快速蒙特卡洛方法预测球磨机磨矿粒度分布的系统及方法 |
CN106874568A (zh) * | 2017-01-17 | 2017-06-20 | 中南大学 | 一种球磨过程的物料粒径分布预测方法 |
CN108763166A (zh) * | 2018-04-17 | 2018-11-06 | 内蒙古科技大学 | 一种粉碎产品粒度分布的拟合方法 |
CN108918356A (zh) * | 2018-05-18 | 2018-11-30 | 中建材(合肥)粉体科技装备有限公司 | 一种混合物料挤压粉碎后的粒度预测方法及系统 |
CN109085809A (zh) * | 2018-08-29 | 2018-12-25 | 东北大学 | 加长型球磨机粒度分布的分布参数模拟系统 |
CN112816371A (zh) * | 2021-04-22 | 2021-05-18 | 南京兴盛德水泥有限公司 | 一种球磨机水泥生粉磨粒度测量系统及其测量方法 |
CN114112819A (zh) * | 2022-01-25 | 2022-03-01 | 北矿智云科技(北京)有限公司 | 一种测量磨矿粒度的方法及装置 |
CN114372635A (zh) * | 2022-01-12 | 2022-04-19 | 中冶沈勘秦皇岛工程设计研究总院有限公司 | 磨矿分级作业的预测方法、装置、电子设备及存储介质 |
CN114733617A (zh) * | 2022-03-10 | 2022-07-12 | 中南大学 | 一种磨矿分级过程的nmpc-pi控制方法及系统 |
CN114372635B (zh) * | 2022-01-12 | 2024-07-05 | 中冶沈勘秦皇岛工程设计研究总院有限公司 | 磨矿分级作业的预测方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3596839A (en) * | 1969-12-10 | 1971-08-03 | Westinghouse Electric Corp | Slurry particle size determination |
JPH10211445A (ja) * | 1997-01-29 | 1998-08-11 | Mitsubishi Chem Corp | 粉粒体の性状制御方法 |
-
2010
- 2010-06-04 CN CN2010101915033A patent/CN101869860B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3596839A (en) * | 1969-12-10 | 1971-08-03 | Westinghouse Electric Corp | Slurry particle size determination |
JPH10211445A (ja) * | 1997-01-29 | 1998-08-11 | Mitsubishi Chem Corp | 粉粒体の性状制御方法 |
Non-Patent Citations (2)
Title |
---|
《无机盐工业》 20070531 粟海锋,徐颖利,胡雪玲,文衍宣,童张法 重钙超细粉碎工艺与粒度分布的定量关系 22-25 1 第39卷, 第5期 2 * |
《淮南矿业学院学报》 19930630 俞少功,姜利群 煤超细研磨过程的数学模型研究 81-86 1 第13卷, 第2期 2 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318330B (zh) * | 2014-10-20 | 2017-06-13 | 东北大学 | 快速蒙特卡洛方法预测球磨机磨矿粒度分布的系统及方法 |
CN104318330A (zh) * | 2014-10-20 | 2015-01-28 | 东北大学 | 快速蒙特卡洛方法预测球磨机磨矿粒度分布的系统及方法 |
CN106874568A (zh) * | 2017-01-17 | 2017-06-20 | 中南大学 | 一种球磨过程的物料粒径分布预测方法 |
CN108763166B (zh) * | 2018-04-17 | 2021-11-19 | 内蒙古科技大学 | 一种粉碎产品粒度分布的拟合方法 |
CN108763166A (zh) * | 2018-04-17 | 2018-11-06 | 内蒙古科技大学 | 一种粉碎产品粒度分布的拟合方法 |
CN108918356A (zh) * | 2018-05-18 | 2018-11-30 | 中建材(合肥)粉体科技装备有限公司 | 一种混合物料挤压粉碎后的粒度预测方法及系统 |
CN108918356B (zh) * | 2018-05-18 | 2021-07-09 | 中建材(合肥)粉体科技装备有限公司 | 一种混合物料挤压粉碎后的粒度预测方法及系统 |
CN109085809A (zh) * | 2018-08-29 | 2018-12-25 | 东北大学 | 加长型球磨机粒度分布的分布参数模拟系统 |
CN112816371A (zh) * | 2021-04-22 | 2021-05-18 | 南京兴盛德水泥有限公司 | 一种球磨机水泥生粉磨粒度测量系统及其测量方法 |
CN114372635A (zh) * | 2022-01-12 | 2022-04-19 | 中冶沈勘秦皇岛工程设计研究总院有限公司 | 磨矿分级作业的预测方法、装置、电子设备及存储介质 |
CN114372635B (zh) * | 2022-01-12 | 2024-07-05 | 中冶沈勘秦皇岛工程设计研究总院有限公司 | 磨矿分级作业的预测方法、装置、电子设备及存储介质 |
CN114112819A (zh) * | 2022-01-25 | 2022-03-01 | 北矿智云科技(北京)有限公司 | 一种测量磨矿粒度的方法及装置 |
CN114733617A (zh) * | 2022-03-10 | 2022-07-12 | 中南大学 | 一种磨矿分级过程的nmpc-pi控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101869860B (zh) | 2012-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101869860B (zh) | 一种球磨机磨矿产品的粒度分布预测方法 | |
CN101339574B (zh) | 基于回弹补偿的混凝土搅拌叶片模具型面设计系统及方法 | |
CN104573237B (zh) | 一种基于摩擦磨损cae分析的模具优化方法 | |
Wang et al. | Investigation of material deformation in multi-pass conventional metal spinning | |
CN105373683A (zh) | 一种20CrMnTiH钢热变形过程微观组织演化规律的预测方法 | |
CN104156783A (zh) | 计及气象累积效应的电力系统最大日负荷预测系统及方法 | |
CN102672003B (zh) | 不锈钢带钢拉矫机机组的工艺参数设定方法 | |
CN102672059B (zh) | 根据仿真冲压工件厚度确定出模具凹凸模修改型面的方法 | |
CN102601881B (zh) | 一种橡胶硬度的在线质量监控和预测模型更新的方法 | |
CN106157176A (zh) | 一种gpu加速的电力潮流雅可比矩阵的lu分解方法 | |
CN106461522A (zh) | 可塑性材料的评价方法及可塑性材料的塑性加工的评价方法 | |
CN104680012A (zh) | 一种烧结配料计算模型 | |
CN101655490B (zh) | 一种不锈钢盘条抗拉强度的预测方法及其控制方法 | |
CN102142054A (zh) | 一种铝合金蒙皮拉伸成形道次划分方法 | |
CN104932464A (zh) | 一种动态配料控制系统设计方法 | |
CN103878181B (zh) | 一种可提高型钢轧辊使用寿命的轧辊孔型优化方法 | |
Zhang et al. | Enhancing the green efficiency of fundamental sectors in China’s industrial system: A spatial-temporal analysis | |
CN103163033A (zh) | 级配碎石回弹模量数值测定方法 | |
CN109085809B (zh) | 加长型球磨机粒度分布的分布参数模拟系统 | |
CN103272982B (zh) | 金属薄壁件铆接装配的铆钉镦粗方向确定方法 | |
CN104318330A (zh) | 快速蒙特卡洛方法预测球磨机磨矿粒度分布的系统及方法 | |
CN101611922A (zh) | 一种提高加料后物料含水率均匀性的方法 | |
CN102383529B (zh) | 一种确定使用弯矩下钢砼梁刚度折减与配筋率关系的方法 | |
CN102063094B (zh) | 轧钢工序机组间的电能分配优化方法 | |
CN102004834B (zh) | 一种预测大锻件轴向中心线上空洞闭合率的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120523 Termination date: 20150604 |
|
EXPY | Termination of patent right or utility model |