6000 1473 95000 1273 4000 1073 3000 873 2000 673 1000 473 273 30 6090120 150 180 Time/min 图5模拟高炉条件下部分反应后复合团块冷抗碎强度的变化 Fig.5 Change of BCB cold crushing strength after partial reaction under simulated BR 3.2模型验证 用表2中的实验条件对复合团块模型进行验证。在模拟每种实验方案下的生物质炭复合团块反应行为 时，由表2和进入反应管的总气体流量确定模型的边界条件。由于氧化铁颗粒的烧结，表3中反应(1-3) 的反应速率中的难于确定，因此，本文采用了试错法。用于确危®值的数据点如图6所示。在每种实 验方案下，数据点以40s的时间步长进行选择。aes的取值为1000m2m、1250m2m3、1500m2m3、1750 m2-m3和2000m2m3。模型中，时刻1的团块质量损失分数串戒 (12计算。 (12) 0.9 Scenario Ill 0.7 0.6 Scenario Il Scenario I 300 600900120015001800 Time /s 图6用于确定as的数据点 Fig.6 Selected data points for determinings 对每个as值，模型预测值和实验测量值之间的一致性通过均方误差(MSE)进行评估，MSE表示为式 (13)。结果表明，在aes=1000m2m3下，最小均方差为0.0053，在ags=1250m2m3下为0.0037，在as= 1500m2m3下为0.0032，在aes=1750m2m3下为0.0031，在ae=2000m2m3时为0.0034。因此可以认为 ag=1750m2m3是最合适的。 MSE=( ∑(ym-vam)/Np (13) 其中，Vm为模拟值，'为实验值，Np为图6中比较数据点总数。图 5 模拟高炉条件下部分反应后复合团块冷抗碎强度的变化 Fig. 5 Change of BCB cold crushing strength after partial reaction under simulated BF conditions 3.2 模型验证 用表 2 中的实验条件对复合团块模型进行验证。在模拟每种实验方案下的生物质炭复合团块反应行为 时，由表 2 和进入反应管的总气体流量确定模型的边界条件。由于氧化铁颗粒的烧结，表 3 中反应（1-3） 的反应速率中的 ags 难于确定，因此，本文采用了试错法。用于确定 ags 值的数据点如图 6 所示。在每种实 验方案下，数据点以 40 s 的时间步长进行选择。ags的取值为 1000 m2 ·m-3、1250 m2 ·m-3、1500 m2 ·m-3、1750 m2 ·m-3和 2000 m2 ·m-3。模型中，时刻 t 的团块质量损失分数由式（12）计算。 t /2 2 m O 1 2 3 C 4 C,0 O,0 0 0 1.0-(4 ( ( ) ) ) / ( + ) d f M R R R M R r drdt m m = + + +    (12) 图 6 用于确定 ags的数据点 Fig. 6 Selected data points for determining ags 对每个ags值，模型预测值和实验测量值之间的一致性通过均方误差（MSE）进行评估，MSE表示为式 （13）。结果表明，在 ags1 1000 m2 ·m-3下，最小均方差为 0.0053，在 ags1 1250 m2 ·m-3下为 0.0037，在 ags 1 1500 m2 ·m-3下为 0.0032，在 ags 1 1750 m2 ·m-3下为 0.0031，在 ags1 2000 m2 ·m-3时为 0.0034。因此可以认为 ags1 1750 m2 ·m-3是最合适的。 2 sim exp ( ( ) ) / NP P i MSE v v N = −  (13) 其中， sim v 为模拟值， exp v 为实验值， NP 为图 6 中比较数据点总数。 录用稿件，非最终出版稿