水泥生产过程数值模拟研究.docx
- 文档编号:14492466
- 上传时间:2023-06-23
- 格式:DOCX
- 页数:18
- 大小:504.23KB
水泥生产过程数值模拟研究.docx
《水泥生产过程数值模拟研究.docx》由会员分享,可在线阅读,更多相关《水泥生产过程数值模拟研究.docx(18页珍藏版)》请在冰点文库上搜索。
水泥生产过程数值模拟研究
水泥生产过程数值模拟研究
史才军曹张
(湖南大学土木工程学院,长沙410082)
摘要:
水泥生产工艺是一个非常复杂的物理变化和化学变化的过程,为了研究这一过程,常采用建立数学模型进行模拟的方
法。
对模拟水泥生产的模型进行了评述,并分析了其优点和缺点。
Mujumdar的RoCKS计算模型是目前最为完备的整体计算模
型,但实际的生产情况与计算模型仍存在一定的差异。
对其中单个工艺环节建立的数学模型有一定的局限性,但其计算结果更
加准确,对工艺的改进和进一步建立完整的数学模型意义很大。
TakakoYuko建立的颗粒模型和DeliangShi建立的热力学模
型,同样都是研究颗粒物理化学活动和能量传递的微观机理,准确的把握微观机理,才能准确全面的研究水泥生产工艺。
文章最后提出一种控制温室气体的排放的有效途径:
二氧化碳的捕捉和储藏技术。
关键词:
水泥;生产;模拟;数学模型;二氧化碳
SIMULATIONOFTHECEMENTPRODUCTION
Abstract:
Keywords:
随着我国基础建设的高速发展,水泥产量和消耗量连续几年来以惊人的速度增长。
到
2009年,水泥产量达到16.3亿吨,约占世界总产量的60%每年水泥生产电耗超过1300亿千瓦时,煤耗占全国煤炭总产量的10%年消耗石灰石约14亿吨,排放CQ约10亿吨以上。
[1-3]
水泥生产工艺分为湿法、半干法和干法,基于干法生产的优点,目前正广泛使用。
⑷水泥的
生产过程是一个很复杂的物理化学过程,考虑到节能和减排等问题,我们必须要研究和优化
生产过程。
典型的工艺流程图如图1,含有CaCQSi02、AI2Q、FezQ等成分的生料陆续的通过预热炉、煅烧炉、回转窑、冷却塔,最终形成水泥熟料。
在预热炉里,原料通过来自于煅烧炉的高温气体得到预热,同时原料充分混合。
在煅烧炉,生料充分吸收热量,部分煅烧分解,热
量来源一部分来自于燃料的燃烧,另外利用来自于回转窑的排放气体和冷却器排放的第三级
高温气体,同时煅烧炉煅烧炉的排放气体也用于预热炉的预热。
在回转窑,生料继续煅烧分
解,部分熔融状态,同时形成C2S、C3S、CA、GAF,在这阶段,热量主要来源于燃料的燃烧,煤粉可以通过第二级气体一起从反方向进入窑体内。
高温熟料进入冷却塔,熟料得到冷却,热量可以通过气体进行回收,回收气体分别进入回转窑和煅烧炉循环利用,部分气体直接排
入大气中。
到目前为止,对该过程很多环节的控制方法仍在研究中。
计算模型很多,按照研究对象
考虑分为单个模型和集成模型,按照研究的方法来看分为热力学分析模型和微观机理模型,按模型的特点分为一维、二维和三维模型,按照研究的目的主要有生产过程的优化、工艺更
新等。
1预热炉
预热分解炉是新型干法水泥生产线上预分解技术的重要设备,其属于高温气-固多相反
应器,炉内为气、固两相流动、煤粉燃烧、生料化学反应、同时伴有旋流、回流的复杂流场。
不仅使传热速度大大加快,而且尽可能提高生料的分解率。
1.1Mujumdar模型[5]
Mujumdar等人,由烘干模型引入到旋风式预热炉的计算中,提出了相应的数学模型。
在该模型中主要考虑能量在气-固之间的传递,不考虑石灰石的分解和煤粉的燃烧。
如图2为i
级预热器,假定从i-1、i级预热器出来的固体质量分别为Mn、M,i,温度为Ti-i、T,从i、
i+1级预热器排出的气体带走的固体质量分别为Me,i、Me,i+1,温度为Ti、Ti+1,每级排出气体
质量为M。
物料平衡方程为:
Ms,i」Mse,i1=Ms,iMse,i
(1)
Mse,i二(1一m,p)Ms,i
(2)
在上式中,nm,p表示i级预热器的质量效率。
i级预热器的能量平衡方程满足:
Ms,i4Cp,sTc,iJMse,i1Cp,sTc,i1MgCp,gTc,iMs,iCp,sTc,iMse,iCp,sTc,i
Tc,ihycAcyi(Tc,i_Tiw,i)
上式中Cp,s和Cp,g表示固体和气体的热容量,符号
这种模型的基本原理是质量守衡和能量守衡,计算的连续性很好,输出结果便于下一阶段的
计算。
其缺点主要有:
(1)这里所指的能量简单的指的是该系统的输入能量,如图1中煅烧炉
的尾气排入预热炉的输入能量,事实上很多预热器内部燃烧产生热;
(2)这里的平衡主要指
的是物理意义上的平衡,没有涉及到物料的分解等化学反应;(3)该模型对预热炉的几何尺
寸仅从外形进行考虑,实验证明模型的几何尺寸和进风口的设置会对结果产生影响。
1.2CFD模型
炉内的气、固两相流动的过程是流体力学范畴,可以近似假定满足理想流体理论,通过
dynamics)模型进行计算,假定模型满足Navier-Stokes
利用建立CFD(Computationalfluent方程[6]:
该模型内部边界条件复杂,可以考虑建立气相湍流流动模型,利用计算流体力学中标准
k-八RNGk-;模型和雷诺应力模型(RSM)求解。
李相国、叶旭初、D.K.Fidaros等就应用
商业软件对分解炉内的情况分别进行冷态流场和燃料燃烧分析,[7-9]得出炉内的流动和温度
分布情况。
该方法优点是直观,方便。
缺点是:
(1)该计算模型的边界条件相对独立,其计
算结果不利于下一步的使用;
(2)模型中主要考虑的是理想流体,炉内的气、固两相流动不
是完全意义上的理想流体;(3)不同的求解模型的计算结果有一定的差异。
2煅烧炉
煅烧炉的数学模型主要是Mujumdar等提出的,问如图3。
Mujumdar等假定物料的部分分解发生在煅烧炉中。
首先作如下假设:
(1)气相完全混合均匀;
(2)原料和煤颗粒看成统一粒径的离散相颗粒,且颗粒在炉内停留时间相同。
得到气相质量平衡方程:
dmgcl
—mgin—mgout+[mp,cin一mp,cout]'Np+[mp,Li^—mp,Lout]'Npdt
气相能量平衡方程:
(6)
该方法的计算精度很大程度上依赖于上面两条假设与实际的吻合程度。
煅烧炉的几何尺寸对
气相的均匀混合影响很大。
3回转窑
水泥回转窑是水泥熟料生产线的主要设备。
回转窑由筒体、支承装置、带挡轮支承装置、
传动装置、活动窑头、窑尾密封装置、喷煤管装置等部件组成。
回转窑的窑体与水平呈一定
的倾斜,整个窑体由托轮装置支承,并有控制窑体上下窜动的挡轮装置,传动系统除设置主
传动外,还设置了在主电源中源断时仍能使窑体转动,防止窑体弯曲变形的辅助传动装置,窑头、窑尾密封装置采用了先进的技术,保证了密封的可靠性。
3.1微观分析模型
3.1.1流化床窑技术
通过分析,I.Maki等人提出熟料的化学组成和颗粒尺寸对我们估计熟料的形成机理很有作用。
[10]在此基础上,TakakoYuko等利用微观分析,[11]比较了喷动床窑、流化床窑和回转窑的差别。
方法是,分别在喷动床窑、流化床窑中取样,用JISR5202方法比较不同窑体
的熟料颗粒中CaOSiO2、f-CaO、AI2Q、F&C4、MgON@OKOSO等的含量。
并借助电子显微镜拍摄颗粒照片和XR技术分析熟料中的矿物含量,从微观机理上解释熟料颗粒的形成过程。
JISR5202方法得到的数据比较:
(1)同样粒径的颗粒,流化床窑熟料颗粒中f-CaO的含量明显低于喷动床窑熟料颗粒;
(2)随着粒径增大,CaO勺含量增大,而SiO2的含量减少;
⑶同样粒径的颗粒,流化床窑熟料颗粒中N@OKOSO的含量明显低于喷动床窑熟料颗
粒。
分析结果表明:
(1)C3S在大颗粒熟料中含量较多,而CS在小颗粒熟料中含量较多;
(2)流化床窑熟料颗粒中有害氧化物含量较少,熟料品质更好。
3.1.2影响颗粒间热传递的因素
DeliangShi等利用离散单元法DEMCF[以及热传导计算模拟回转窑内的热量传递情况。
网并用表1实验参数进行实验,研究对比当颗粒材料为玻璃、铁、铝时的各自不同的传热特性,总结出当介质热传导率低或在真空内,热传导系数高的材料传热快;但当介质热传
导率提高后,热传导系数低的材料传热效率大大提高。
表1实验参数
Parametersoftheexperiment
Glass
Steel
Aluminum
Particle
Density(kg/m3)
2700
7900
2700
Diameter(m)
3.0*103
3.0*103
3.0*103
Poissonratio
0.33
0.29
0.33
Young'smoduIus(GPa)
69
193
110
Frictioncoefficient
0.3
0.3
0.3
Thermalconductivity(W/(mK))
1.7
15
180
Heatcapacity(J/(kgK))
800.0
477.0
900.0
Air
Density(kg/m3)
1.23
Viscisitykg/(ms)
1.8*10-5
Thermalconductivity(W/(mK))
0.0261
Heatcapacity(J/(kgK))
1000.7
A.A.BOATENG等利用热力学模型,[13]通过研究研究回转窑内热传导情况和材料颗粒的运动轨迹,得出回转窑转速和内部热传导情况之间的规律。
微观分析方法的缺点是着眼点在局部或某一方面,不能够整体了解到窑内的物料平衡和
能量平衡状况。
但微观分析模型是深入了解窑内的情况必不可少的环节。
TakakoYuko等的
微观分析方法从根本上研究了熟料的生成过程机理,该分析方法有助于我们分析窑内环境和
熟料的组成成分之间的关系。
DeliangShi等的物料热传递分析,AABOATENG等的热力学
模型分析,都没有考虑实际窑内出现的熔融状态,提出了研究颗粒之间热传递的研究方法和
影响热传递效率的主要因素,其模型和分析方法对窑内的整体分析起到借鉴作用。
3.2宏观分析模型
3.2.1Saeman模型
Saeman在1951年建立了回转窑的数学模型,通过该模型可以计算出物料在回转窑内停
留的时间与相关参数的关系如下式Q为体积速率。
[14]
3.2.2Mujumdar模型
3.2.3
基于Sherritt等人提出的关于回转窑轴向混合情况的估计,阿结果表明对于工业水泥
dx
化学反应由下式计算:
该模型是一个典型的一维数学模型,明确的边界条件,可以得到明确的解。
3.2.3回转窑内热传递模型
M.D.Heydenrych等人在基于假设的基础上提出一种新的计算模型。
[18]该模型主要考虑
计算结果同Jauharietal的实验结果进行比较,昭发现该模型对于快速反应系统误差较
大。
即:
k
「k<2模型适用
Z
\k
k>2误差较大
原因是因为,在快速反应系统中,分散或扩散对结果的影响不可忽略。
324CFD模型
E.Mastorakos应用轴对称CFD模型,蒙特卡罗方法建立辐射模型,有限体积法建立能
量方程,同时考虑能量平衡、化学反应。
借助FLOW-3D商业软件进行求解。
[20]其结果与实
验数据有一定的偏差,窑外温度偏差普遍达10%^上,个别数据偏差40%根据计算模型结果
预测,回转窑向外界损失热量为10%左右,用以承担生料分解和熟料生成的热量占40以上。
CFD模型计算方法的计算精度很大程度取决于模型的设置,对于回转窑的数值模拟,CFD
模型还有待进一步改进。
4冷却塔
从窑内出来的水泥熟料,温度很高,处于高温熔融的状态,在急冷急热作用下,来不及结晶而变成玻璃体状的物质。
玻璃体含量多时,可提高熟料质量,易磨性好,使熟料强度提高。
基于NsoforandAdebiyi实验确定可以适用于200-1000摄氏度的气固热传递系数,[21]Mujumdar等人建立了冷却系统的数学模型如图6。
冋其中假定水泥熟料稳定的流量速度进入
该系统,空气由垂直方向即y方向进入,气体按照第二级气体、第三级气体和尾气分别排入
回转窑、煅烧炉和大气中。
这三者各自所占比例按照模型图上各自x方向所占的单元数来确
定。
模型的边界条件如图所示。
质量平衡方程可以写为:
d叫(X)=0
dx
能量平衡方程可以写为:
-「1(1-;)us,xCp,sTs).:
(;(1-;)us,yCp,sTs)
-x:
y
=虫(1-;)ks:
:
Ts/G去(1-;)ks汀s/:
:
y]
:
x:
:
y
以上两式中hc,c为熟料和空气之间的热传递系数,NsoforandAdebiyi通过实验确定
Nusseltnumber计算的经验公式:
Nu=8.749.346(1-;)0.2Re0'2Pr0'33(30vNlk60;50vNlk120)(12)
需要指出的是,对于大多冷却系统,其雷诺数Reynoldsnumbers比NsoforandAdebiyi实验环境要高很多。
如何得到适用高温状态下传递系数,仍需要进一步的实验研究。
5集成模型
要总体考虑水泥生产工艺,需要把前面独立的四部分模型联系在一起,组合成一个整体
系统。
[5,17]Mujumdar等人建立了整体计算模型RoCKS,如图7示。
模型一共包括前处理、调
用子模型、迭代计算、后处理等部分。
前处理主要功能是:
输入需要的数据、数据一致性检测和初始条件假定。
数据输入主要有:
(1)生料从预热炉进入;
(2)空气从冷却设备进入;(3)
煤由煅烧炉和回转窑进入;(4)所有材料特性和必须的参数设置。
通过计算,Mujumdar等人得到了预热炉级数及生料分解情况与单位产量能耗之间的关系,回转窑转速、冷却过程及煤的特性对熟料生成的影响。
虽然所得结果细节上与实际生产
有差距,但分析所得到一些规律与实际生产仍然是基本吻合。
以上提到1-4的单一部分的计算模型相对整体集成模型来说较为简单,虽然不能够整体的把握整个工艺流程,但其计算结果仍然有助于对水泥生产过程的研究和优化。
如何在计算
数学模型中考虑更多的因素,以使计算模型与实际生产的差距越来越小,仍是需要解决的问
题。
特别是对于水泥熟料的微观分析,有助于进一步把握水泥生成机理,个人认为意义重大。
整体集成模型的研究目前基本停留在一维研究阶段,如何对整个过程进行更准确的模拟,随
着计算流体力学CFD模型的引入,为下一步的研究工作提供了参考。
图8二氧化碳捕捉的三种工艺示意图
6二氧化碳的捕捉储藏以及纯氧燃烧
6.1二氧化碳的捕捉和储藏
鉴于温室气体CQ排放量太大,对我们地球环境所造成的巨大压力,从可持续发展战略的角度出发,西方学者首先提出CQ的捕捉和储存的设想。
剑桥大学XiLiang曾指出CQ
的捕捉和储存是化石能源安全利用的唯一途径。
主要原理是利用一定的工艺方法把工业尾
气中的CQ捕捉起来,在满足工业CQ需求后,剩余的CQ经过压缩以及脱水干燥等方法储存起来。
如图8示为CQ捕捉的三种不同工艺,即燃烧前除碳、燃烧后捕捉和纯氧燃烧工艺。
图9是CQ的捕捉和储存基本流程图,捕捉后对CQ进行储存。
[22-24]
6.2Q2/CQ2燃烧
纯氧燃烧为捕捉提供了更大的可行性。
C2/CQ2燃烧,有关文献又称为纯氧燃烧,[25,26]目
前主要国外对该领域的研究很多,特别是西欧已进入实验阶段,预计不久的未来会进入生产
实践阶段。
[27-30]与传统水泥生产流程不同,在燃烧前,将空气中的Q和Na进行分离,纯氧
燃烧一方面减少因N产生的热量损失,另一方面提高CQ在尾气中的含量,更利于CQ的提
取。
CQ的捕捉和储存要想全面得到推广应用,目前除一些技术性问题需要解决外,就是经济上需要降低成本。
[31]
7结语
目前为止,还没有一种广为接受的准确的模拟水泥生产工艺的数学模型。
水泥生产工艺
是一个非常复杂的物理变化和化学变化的过程,为了研究这一过程,常采用建立数学模型进
行模拟的方法。
精确的了解水泥熟料的形成过程,需要准确的建立与实际情况相吻合的数学
模型。
进一步控制熟料的特性、研究如何在生产过程中节能、如何控制温室气体的排放等,这对于水泥工业的发展和可持续发展战略有着重要的意义。
控制温室气体的排放的有效途径
是实行二氧化碳的捕捉和储藏技术。
水泥生产是一个连续工艺的过程,Mujumdar等人建立的RoCKS计算模型是目前最为完
备的整体计算模型,但实际的生产情况与计算模型仍存在一定的差异;建立CFD模型进行
计算,是另外一种可以尝试的计算方法,但它也存在模型过于理想化,数值计算量大。
对其中单个工艺环节建立的数学模型有很大的局限性,但因为研究对象相对单一,考虑
的因素更多,其计算结果更加准确,对工艺的改进和进一步建立完整的数学模型意义很大。
TakakoYuko建立的颗粒模型和DeliangShi建立的热力学模型,同样都是研究颗粒物理
化学活动和能量传递的微观机理,准确的把握微观机理,才能准确全面的研究水泥生产工艺。
References:
[1]DaisukeHayashi,MatthiasKrey.CQ2emissionreductionpotentialoflarge-scaleenergy
efficiencymeasuresinheavyindustryinChina,India,Brazil,IndonesiaandSouthAfrica[J].HamburgInstituteofInternationalEconomics(HWWI).2005.
[2]蒋家超,李明,赵由才等编.工业领域温室气体减排与控制技术[M].化学工业出版社,
2009.
[3]白冰,李小春,刘延锋,张勇.中国CQ2集中排放源调查及其分布特征.岩石力学与
工程学报.2006
(2):
2918-2923
[4]马保国主编.新型干法水泥生产工艺.化学工业出版社.2007.
⑸KaustubhS.Mujumdar,K.V.Ganesh,SaritaB.Kulkarni,VivekV.RanadeRotaryCememnt
KilnSimulator(RoCKS):
lntegratedmodelingofpre-heater,calciner,kilnandclinkercooler[J].ChemicalEngineeringScienee62,(2007)2590-2607
⑹严敬等编.工程流体力学[M].重庆大学出版社.2007.
[7]李相国•预分解系统内流动、燃烧与分解的研究及数值模拟•武汉理工大学.2006.
[8]叶旭初,李祥东,胡道和.喷旋管道式分解炉内燃烧、分解过程的CFD莫拟[J].南京工
业大学学报.2006,28
(1):
62-66
[9]D.K.Fidaros,C.A.Baxevanou,C.D.Dritselis,N.S.Vlachos.Numericalmodellingofflowandtransportprocessesinacalcinerforcementproduction[J].PowderTechnology2007(171):
81-95
[10]I.Maki,T.Tanioka,T.Hibino.FormationofPortlandclinkernodulesinrotarykilnsandfinetexturesofconstituentphases[J].theSecondCem86
(1)(1989)3-10
[11]TakakoYuko,TatsuoIkabata,TatsushiAkiyama,TakanoriYamamoto,NorimitsuKurumada.Newclinkerformationprocessbythefluidzedbedkilnsystem[J].CementandConcreteResearch30,(2000)1113-1120
[12]Deliangshi,WatsonL.Vargas,J.J.McCarthy.Heattransferinrotarykilnswithinterstitialgaser[J].ChemicalEngineeringScienee63,(2008)4506-4516
[13]A.A.BOATENGP.V.BARR.Athermalmodelfortherotarykilnincludingheattransferwithinthebed[J].Int.J.HeatMassTransfer,Vol39,(1996)2131-2147
[14]Saeman,W.C..Passageofsolidsthroughrotarykilns:
factorsaffectingtimepassage[J].ChemicalEngineeringProgress47,508
[15]L.G.Ndiaye,S.Caillat,A.Chinnayya,D.Gambier,B.Baudoin.Applicationofthedynamicmodelofsaemantoanindustryrotarykilnincinerator:
Numericalandexperimentalresults[J].WasteManagement30,(2010)1188-1195
[16]Sherrit,R.G.,Chaouki,J.,Mehrotra,A.K.andBehie,L.A..Axialdispersioninthethree-dimensionalmixingofparticlesinarotatingdrumreactor[J].ChemicalEngineeringScienee,58
(2)401-415
[17]K.S.Mujumdar,V.V.Ranade.Simulationofrotarycementkilnsusingaone-dimensionalmodel[J].ChemicalEngineeringResearchandDesign,84(A3):
165-177
[18]M.D.Heydenrych,P.Greeff,Mehrotra,A.B,M,Heesink,G.F.Versteeg.Masstran
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 水泥 生产过程 数值 模拟 研究