基于传热反问题的高炉炉衬单层侵蚀模型

工程科学学报，第39卷，第4期：574580,2017年4月 Chinese Journal of Engineering,Vol.39,No.4:574-580,April 2017 DOI:10.13374/j.issn2095-9389.2017.04.013:http://journals.ustb.edu.cn 基于传热反问题的高炉炉衬单层侵蚀模型 苏福永区，温治 北京科技大学能源与环境工程学院，北京100083 ☒通信作者，E-mail:sfyong(@ustb.edu.cm 摘要基于传热反问题，建立了高炉炉衬侵蚀过程的数学模型，确定了模型的边界条件，并采用共轭梯度法将反问题分解 为三个问题：正问题、灵敏度问题和伴随问题进行求解.通过不同形状函数的反演结果证明了其可行性，并分别研究初始猜 测形状曲线、测点数等对反演结果的影响.研究结果表明，初始猜测曲线的选取对反演结果影响很小，充分说明该方法不受 初始猜测曲线的限制，具有较好的通用性.而测点数的选取对反演结果有一定的影响，测点数越多，曲线特征被捕捉的越好 但在保证得到曲线特征的前提下，较少的测点数也能得到比较满意的反演结果，平均相对误差控制在3%以内 关键词炉村侵蚀：边界形状：导热反问题：共轭梯度法 分类号TF549 Blast furnace layer lining erosion problem based on inverse the heat conduction model SU Fu-yong,WEN Zhi School of Energy and Environmental Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:sfyong@ustb.edu.cn ABSTRACT A blast furnace lining mathematical model was established based on the inverse heat transfer problem.After determi- ning the boundary conditions of the model,this inverse heat transfer problem is divided into three problems which are the direct prob- lem,the sensitivity problem and the adjoint problem,and these were solved using the conjugate gradient method.The feasibility of this model was proved by the inversion results of different shape functions and then it was discussed that the initial guess shape and number of measurement points effect on the inversion results.The results show that the accuracy of the inverse solution is independent of the the initial guess shape,but the number of measurement points has some impact on these results,whereby the more points are meas- ured,the better the curve features are captured.An accurate inverse solution can be obtained with fewer measurement points and an average relative error within 3%,even though the arrangement of more points can achieve a slightly better solution. KEY WORDS lining erosion:boundary shape:inverse heat conduction problem:conjugate gradient method 高炉炼铁是钢铁生产中的一个重要环节，由于高 导热反问题(inverse heat conduction problem,H- 炉炉衬时时刻刻受到高温、炉气、炉尘、炉渣、铁水等的 CP)B阿产生于许多实际测量比较困难的传热问题. 热应力、化学作用和剧烈冲刷，结果导致炉衬严重侵 自Solz囿于1960年发表了第一篇关于导热反问题的 蚀Ⅲ.然而，如何能在正常生产条件下，通过监测和检 论文以来，许多科学工作者对传热反问题的数值解法 测与高炉炉衬侵蚀厚度相关的其他热工参数来及时预 进行了研究，提出了各种有效的求解方法，如Tikhonov 报高炉炉衬内的侵蚀情况，继而采取相应的措施来延 与Arsenin)的正则化法(1977年)、Beck的序列化法 长高炉炉衬寿命是摆在我们面前亟待解决的重大 (1985年)、Blackwell的时间推进算法图(1981年)、 课题四， Weber的空间推进算法(1981年)、Elden等的Fou- 收稿日期：2016-08-11 基金项目：中央高校基本科研业务费专项资金资助项目(FRF-TP-15O22A3)

工程科学学报，第 39 卷，第 4 期: 574--580，2017 年 4 月 Chinese Journal of Engineering，Vol． 39，No． 4: 574--580，April 2017 DOI: 10． 13374 /j． issn2095--9389． 2017． 04． 013; http: / /journals． ustb． edu． cn 基于传热反问题的高炉炉衬单层侵蚀模型 苏福永，温 治 北京科技大学能源与环境工程学院，北京 100083 通信作者，E-mail: sfyong@ ustb． edu． cn 摘 要 基于传热反问题，建立了高炉炉衬侵蚀过程的数学模型，确定了模型的边界条件，并采用共轭梯度法将反问题分解 为三个问题: 正问题、灵敏度问题和伴随问题进行求解． 通过不同形状函数的反演结果证明了其可行性，并分别研究初始猜 测形状曲线、测点数等对反演结果的影响． 研究结果表明，初始猜测曲线的选取对反演结果影响很小，充分说明该方法不受 初始猜测曲线的限制，具有较好的通用性． 而测点数的选取对反演结果有一定的影响，测点数越多，曲线特征被捕捉的越好． 但在保证得到曲线特征的前提下，较少的测点数也能得到比较满意的反演结果，平均相对误差控制在 3% 以内． 关键词 炉衬侵蚀; 边界形状; 导热反问题; 共轭梯度法 分类号 TF549 Blast furnace layer lining erosion problem based on inverse the heat conduction model SU Fu-yong ，WEN Zhi School of Energy and Environmental Engineering，University of Science and Technology Beijing，Beijing 100083，China Corresponding author，E-mail: sfyong@ ustb． edu． cn ABSTＲACT A blast furnace lining mathematical model was established based on the inverse heat transfer problem． After determi￾ning the boundary conditions of the model，this inverse heat transfer problem is divided into three problems which are the direct prob￾lem，the sensitivity problem and the adjoint problem，and these were solved using the conjugate gradient method． The feasibility of this model was proved by the inversion results of different shape functions and then it was discussed that the initial guess shape and number of measurement points effect on the inversion results． The results show that the accuracy of the inverse solution is independent of the the initial guess shape，but the number of measurement points has some impact on these results，whereby the more points are meas￾ured，the better the curve features are captured． An accurate inverse solution can be obtained with fewer measurement points and an average relative error within 3% ，even though the arrangement of more points can achieve a slightly better solution． KEY WOＲDS lining erosion; boundary shape; inverse heat conduction problem; conjugate gradient method 收稿日期: 2016--08--11 基金项目: 中央高校基本科研业务费专项资金资助项目( FＲF--TP--15--022A3) 高炉炼铁是钢铁生产中的一个重要环节，由于高 炉炉衬时时刻刻受到高温、炉气、炉尘、炉渣、铁水等的 热应力、化学作用和剧烈冲刷，结果导致炉衬严重侵 蚀［1］． 然而，如何能在正常生产条件下，通过监测和检 测与高炉炉衬侵蚀厚度相关的其他热工参数来及时预 报高炉炉衬内的侵蚀情况，继而采取相应的措施来延 长高炉炉衬寿命是摆在我们面前亟待解决的重大 课题［2］． 导热反问题( inverse heat conduction problem，IH￾CP) ［3--5］产生于许多实际测量比较困难的传热问题． 自 Stolz［6］于 1960 年发表了第一篇关于导热反问题的 论文以来，许多科学工作者对传热反问题的数值解法 进行了研究，提出了各种有效的求解方法，如 Tikhonov 与 Arsenin ［7］的正则化法( 1977 年) 、Beck 的序列化法 ( 1985 年) 、Blackwell 的 时 间 推 进 算 法［8］ ( 1981 年) 、 Weber 的空间推进算法［9］( 1981 年) 、Eldén 等的 Fou-

苏福永等：基于传热反问题的高炉炉衬单层侵蚀模型 575 rier正则化法@(2000年)、目前国内外流行的共轭梯 度法-四等，这些方法极大地推动了导热反问题在实 右侧壁面 际工程和科学实验测量中的应用.共轭梯度法 左侧壁面 To (conjugate gradient)是介于最速下降法与牛顿法之间 的一个方法，它仅需利用一阶导数信息，但克服了最速 0- 下降法收敛慢的缺点，又避免了牛顿法需要存储和计 r=f 算Hesse矩阵并求逆的缺点，具有步收敛性和稳定性 ●- 高的特点 热电偶 郑坤灿等提出了基于传热反问题的高炉炉衬 侵蚀厚度的在线动态实时预报方法，并提出了详细的 技术方案和实施策略，为高炉炉衬侵蚀厚度的在线实 图1二维计算区域示意图 时预报提供了新的解决思路，但其模型主要基于高炉 Fig.1 Two-dimensional schematic diagram of the calculation area 炉衬表面温度的测量值，容易受到外部环境的干扰 1.1正问题 张伟等闭分析了当前高炉炉衬检测技术，介绍了冲击 在假设已知内壁侵蚀曲线的条件下，通过求解二 回波法的基本原理和测试特点，并应用此方法对炉缸 维稳态导热问题来获得整个区域的温度场分布.正问 炉衬等进行了厚度测试研究及应用，其试验结果表明 题的数值求解采用有限差分的方法，半隐式差分格式， 冲击回波法能为高炉炉衬厚度监测提供一个新的有效 离散方程采用三对角矩阵算法(ridiagonal matrix algo- 手段，但目前该方法还并未进行大规模应用.鞠贵冬 ithm,TDMA)进行计算.其控制方程及边界条件 等网基于计算流体力学(C℉D)方法，建立了一种考虑 如下： 冷却水进出口温差的冷却壁数值模型，并利用热态试 验研究确立了更为准确的数值模拟边界条件，全面模 =0 C2: 拟了冷却壁运行时的各项指标。左海滨等网建立了 -k x=b; 高炉冷却壁稳态传热模型，模拟了不同材质不同工况 下的温度场分布，并分析了渣皮对冷却壁体起温度的 T=To, x=f(y）: 影响和保护作用.石琳等利用非稳态热传导反问 0 在y=0和在y=L 题数学模型研究了高炉监测热电偶数据与高炉炉衬侵 ay 独厚度的关系，采用了整个时间域上的空间距离步进 T(xy =T, i=1,2,…,m 迭代算法和离散软化过滤法，模拟计算得出高炉炉衬 (1) 剩余厚度随时间的变化关系，为高炉炉衬侵蚀状况的 式中：k为材料热导率，W·m1·K;9为外壁面热流 离线模拟研究做出了贡献. 强度，W·m2;T。为内壁表面温度，K;T为测点i处的 本文基于高炉炉衬侵蚀问题，试图依靠炉衬内部 实际测得的温度值，K:2为计算区域内的坐标合集 的热电偶的测量数据来反推炉衬侵蚀轮廓线，即求解 1.2反问题 边界形状导热反问题 反问题是通过若干测量点处温度值来反推内壁边 1数学模型 界形状，即测量点实际测得的温度值与猜测边界形状 x=∫(y)通过式(1)求解正问题计算得到的测点处温 本文建立的数学模型主要是求解边界形状导热反 度值差的平方的最小值： 问题（即高炉炉衬的内部侵蚀形状曲线）.计算区域如 J(f(y))=T:-T2= 图1所示，研究的是二维无内热源的稳态导热问题 三rW-T1 (2) 左侧内壁面上温度保持恒定值T。,而右侧外壁面上 有一个恒定热流值9。，上下壁面为绝热条件，初始形 式中：T(xy:为测点i处计算得出的温度值，K; 状曲线函数为x=∫(y),计算区域坐标合集为2.在 T(x,y)为测点i处实际测得的温度值，K 左侧壁面或接近左侧壁面处布置若干个热电偶，通 1.3灵敏度问题 过此热电偶测出的温度值来反演内壁面真实的侵蚀 灵敏度问题就是指当边界形状∫(y)有一个增量 曲线 △f(y)时，测点处温度T(x,y)也会有一相应的变化量 本文采用共轭梯度法求解边界形状导热反问题， △T(x,y),其变化量值为多少可以通过求解灵敏度问 计算过程分解为三个问题的求解，分别为正问题、灵敏 题方程，具体方法就是在式(1)中，用T+△T代替原来 度问题和伴随问题. 的T,而f(y)+△f(y)代替原来的f(y),然后减去式

苏福永等: 基于传热反问题的高炉炉衬单层侵蚀模型 rier 正则化法［10］( 2000 年) 、目前国内外流行的共轭梯 度法［11--12］等，这些方法极大地推动了导热反问题在实 际工程和科学实验测量中的应用． 共轭梯度法［13--15］ ( conjugate gradient) 是介于最速下降法与牛顿法之间 的一个方法，它仅需利用一阶导数信息，但克服了最速 下降法收敛慢的缺点，又避免了牛顿法需要存储和计 算 Hesse 矩阵并求逆的缺点，具有步收敛性和稳定性 高的特点． 郑坤灿等［16］提出了基于传热反问题的高炉炉衬 侵蚀厚度的在线动态实时预报方法，并提出了详细的 技术方案和实施策略，为高炉炉衬侵蚀厚度的在线实 时预报提供了新的解决思路，但其模型主要基于高炉 炉衬表面温度的测量值，容易受到外部环境的干扰． 张伟等［17］分析了当前高炉炉衬检测技术，介绍了冲击 回波法的基本原理和测试特点，并应用此方法对炉缸 炉衬等进行了厚度测试研究及应用，其试验结果表明 冲击回波法能为高炉炉衬厚度监测提供一个新的有效 手段，但目前该方法还并未进行大规模应用． 鞠贵冬 等［18］基于计算流体力学( CFD) 方法，建立了一种考虑 冷却水进出口温差的冷却壁数值模型，并利用热态试 验研究确立了更为准确的数值模拟边界条件，全面模 拟了冷却壁运行时的各项指标． 左海滨等［19］建立了 高炉冷却壁稳态传热模型，模拟了不同材质不同工况 下的温度场分布，并分析了渣皮对冷却壁体起温度的 影响和保护作用． 石琳等［20］利用非稳态热传导反问 题数学模型研究了高炉监测热电偶数据与高炉炉衬侵 蚀厚度的关系，采用了整个时间域上的空间距离步进 迭代算法和离散软化过滤法，模拟计算得出高炉炉衬 剩余厚度随时间的变化关系，为高炉炉衬侵蚀状况的 离线模拟研究做出了贡献． 本文基于高炉炉衬侵蚀问题，试图依靠炉衬内部 的热电偶的测量数据来反推炉衬侵蚀轮廓线，即求解 边界形状导热反问题． 1 数学模型 本文建立的数学模型主要是求解边界形状导热反 问题( 即高炉炉衬的内部侵蚀形状曲线) ． 计算区域如 图 1 所示，研究的是二维无内热源的稳态导热问题． 左侧内壁面上温度保持恒定值 T0，而右侧外壁面上 有一个恒定热流值 q0，上下壁面为绝热条件，初始形 状曲线函数为 x = f( y) ，计算区域坐标合集为 Ω． 在 左侧壁面或接近左侧壁面处布置若干个热电偶，通 过此热电偶测出的温度值来反演内壁面真实的侵蚀 曲线． 本文采用共轭梯度法求解边界形状导热反问题， 计算过程分解为三个问题的求解，分别为正问题、灵敏 度问题和伴随问题． 图 1 二维计算区域示意图 Fig． 1 Two-dimensional schematic diagram of the calculation area 1. 1 正问题 在假设已知内壁侵蚀曲线的条件下，通过求解二 维稳态导热问题来获得整个区域的温度场分布． 正问 题的数值求解采用有限差分的方法，半隐式差分格式， 离散方程采用三对角矩阵算法( tridiagonal matrix algo￾rithm，TDMA) 进 行 计 算． 其 控 制 方 程 及 边 界 条 件 如下:   ( x k T  ) x +   ( y k T  ) y = 0， Ω; － k T x = q0， x = b; T = T0， x = f( y) ; T y = 0， 在 y = 0 和在 y = L; T( xi，yi ) = Tc i ， i = 1，2，…，            m． ( 1) 式中: k 为材料热导率，W·m － 1·K － 1 ; q0 为外壁面热流 强度，W·m － 2 ; T0为内壁表面温度，K; Tc i 为测点 i 处的 实际测得的温度值，K; Ω 为计算区域内的坐标合集． 1. 2 反问题 反问题是通过若干测量点处温度值来反推内壁边 界形状，即测量点实际测得的温度值与猜测边界形状 x = f( y) 通过式( 1) 求解正问题计算得到的测点处温 度值差的平方的最小值: J( f( y) ) = ‖Ti － Tc i ‖2 = ∑ m i = 1 ［Ti ( xi，yi ) － Tc i ( xi，yi ) ］2 ． ( 2) 式中: Ti ( xi，yi ) 为测 点 i 处 计 算 得 出 的 温 度 值，K; Tc i ( xi，yi ) 为测点 i 处实际测得的温度值，K． 1. 3 灵敏度问题 灵敏度问题就是指当边界形状 f( y) 有一个增量 Δf( y) 时，测点处温度 T( x，y) 也会有一相应的变化量 ΔT( x，y) ，其变化量值为多少可以通过求解灵敏度问 题方程，具体方法就是在式( 1) 中，用 T + ΔT 代替原来 的 T，而 f( y) + Δf( y) 代替原来的 f( y) ，然后减去式 · 575 ·

576 工程科学学报，第39卷，第4期 (1)可得以下灵敏度问题方程： 三(T,-T)4T: =0,C2: 0n= (9) ayay =0, x=b; 三m- ax 把以上整理归纳可得出共轭梯度法求解步骤 47=4yr, x=f(y):(3) 如下. (1)先任意猜测一个内壁形状初始值f(y),通过 aAT =0, y=0: 求解正问题，即式(1)，得出整个温度分布场，从而获 ay 得测点处的计算值T(xy). =0, l ay y=L. (2)根据T(x,y)和T,(x,y),判断式(2)的值 1.4 伴随问题 是否满足收敛停止条件，即J((y))<e.若满足，说 伴随问题即泛函求导问题，是为了得到目标函数 明初始猜测值∫(y)即为所得形状函数.否则，进行下 的梯度方向J(》].为了导出伴随问题方程，可把 一步. 式(1)乘上拉格朗日算子A(x,y)(即伴随函数)，并对 (3)求解伴随问题，即式(4)，得出整个区域的伴 其区域积分，然后加入到式(2)中，再用T+△T代替 随函数值，从而利用正问题的温度场分布解，通过式 T,f(y)+△(y)代替f(y),经过一系列变形，得到变分 (5)求出目标函数的梯度值J(f(y)). △了(y)].最后利用灵敏度问题的边界条件，使 (4)第一次求解时，可直接利用式(6)求出搜索方 △J(]趋近于0，得到伴随问题的方程： 向.否则，先利用式(7)求出共轭系数，再通过式(6) 求出搜索方向. [A+A=0 C2: ax oy (5)使△f=d。,通过求解灵敏度问题，即式(3)得 出温度增量△T. =2(T-T)8(y-y)，x=b: dy (6)利用式(9)算出搜索步长，再利用式(8)计算 入=0， x=f(y）: (4) 得到新的形状函数∫(y).之后，利用式(1)求出新的 a以=0， y=0: T(xy),回到第二步. ay a以=0， 2 模型验证及仿真结果分析 ay y=L. 2.1模型验证条件 解得伴随函数入(x,y)的相应值后，利用下式求出 本文以高炉炉衬侵蚀为研究背景，取二维炉衬为 目标函数的导数： 研究区域，炉衬内壁面温度T。=1150℃，即高炉炉渣 J(V(y))=-a证 (5) ax dx ls-f) 熔点值.炉衬初期没有侵蚀时，x轴方向上的厚度b= 1.5迭代搜索方法 0.5m,y轴方向上的长度L=1m其测点均匀布置在 共轭梯度法的基本思想是把共轭性和最速下降法 外壁面处，个数为M.此研究区域二维尺寸为0.5m× 相结合，利用已知点处的梯度构造一组共轭方向，并沿 1m,网格数选取50×100，网格大小为1cm,这就意味 此组方向进行搜索，从而求出目标函数的极小值.此 着形状函数最小变化量为1cm.网格数的选取是综合 方向称为搜索方向，公式如下： 考虑了形状函数的变化量和程序运行时间得出的.计 「d=-J, n=1; 算中热流密度9。取值50000W·m2,初始曲线全部采 (6) 用的是形状函数x=f(y）)=0. ldn=-Jg+B.d.-,n≥2. 2.2计算结果分析 式中：d为第n次的搜索方向，第一次搜索方向直接取 (1)第一种真实形状函数方程： 负梯度方向：J:为梯度方向：B。为共轭系数.共轭系数 x=f(y）=0.35-0.1sin(2πy）,0≤y≤1.(10) Bn可由下式计算： 在真实的形状函数下，右侧固定位置上的热电偶 (J)dy 实测值见表1. B. Jo ，n≥2. (7) 利用这11个测点的实测值，可以通过共轭梯度 (J)2dy 法，把求解边界形状反问题转化成不断求解正问题、灵 形状函数(y)的变化过程可由下式得到： 敏度问题及伴随问题的循环，在满足停止收敛标准下， fni=fn+ad。 (8) 得到最终的反演曲线，它与真实侵蚀边界形状函数的 式中，a为搜索步长.可通过下式求得： 比较如图2所示.反演得到的曲线是由相邻两个反演

工程科学学报，第 39 卷，第 4 期 ( 1) 可得以下灵敏度问题方程:   ( x k ΔT  ) x +   ( y k ΔT  ) y = 0， Ω; － k ΔT x = 0， x = b; ΔT = Δf T x ， x = f( y) ; ΔT y = 0， y = 0; ΔT y = 0， y =              L． ( 3) 1. 4 伴随问题 伴随问题即泛函求导问题，是为了得到目标函数 的梯度方向 J'［f( y) ］． 为了导出伴随问题方程，可把 式( 1) 乘上拉格朗日算子 λ( x，y) ( 即伴随函数) ，并对 其区域积分，然后加入到式( 2) 中，再用 T + ΔT 代替 T，f( y) + Δf( y) 代替 f( y) ，经过一系列变形，得到变分 ΔJ［f ( y) ］． 最后 利 用 灵 敏 度 问 题 的 边 界 条 件，使 ΔJ［f( y) ］趋近于 0，得到伴随问题的方程:  2 λ x 2 +  2 λ y 2 = 0， Ω; － λ y = 2( Ti － Tc i ) δ( y － yi ) ， x = b; λ = 0， x = f( y) ; λ y = 0， y = 0; λ y = 0， y =            L． ( 4) 解得伴随函数 λ( x，y) 的相应值后，利用下式求出 目标函数的导数: J'( f( y) ) = － λ x T x x = f( y) ． ( 5) 1. 5 迭代搜索方法 共轭梯度法的基本思想是把共轭性和最速下降法 相结合，利用已知点处的梯度构造一组共轭方向，并沿 此组方向进行搜索，从而求出目标函数的极小值． 此 方向称为搜索方向，公式如下: d1 = － J'1， n = 1; dn = － J' { n + βn dn － 1， n≥2． ( 6) 式中: dn为第 n 次的搜索方向，第一次搜索方向直接取 负梯度方向; J'n为梯度方向; βn为共轭系数． 共轭系数 βn可由下式计算: βn = ∫ L 0 ( J'n ) 2 dy ∫ L 0 ( J'n －1 ) 2 dy ，n≥2． ( 7) 形状函数 f( y) 的变化过程可由下式得到: fn + 1 = fn + αn dn ． ( 8) 式中，αn为搜索步长． 可通过下式求得: αn = ∑ m i = 1 ( Ti － Tc i ) ΔTi ∑ m i = 1 ( Ti － Tc i ) 2 ． ( 9) 把以上整理归纳可得出共轭梯度法求解步骤 如下． ( 1) 先任意猜测一个内壁形状初始值 f( y) ，通过 求解正问题，即式( 1) ，得出整个温度分布场，从而获 得测点处的计算值 Ti ( xi，yi ) ． ( 2) 根据 Tc i ( xi，yi ) 和 Ti ( xi，yi ) ，判断式( 2) 的值 是否满足收敛停止条件，即 J( f( y) ) ＜ ε． 若满足，说 明初始猜测值 f( y) 即为所得形状函数． 否则，进行下 一步． ( 3) 求解伴随问题，即式( 4) ，得出整个区域的伴 随函数值，从而利用正问题的温度场分布解，通过式 ( 5) 求出目标函数的梯度值 J'( f( y) ) ． ( 4) 第一次求解时，可直接利用式( 6) 求出搜索方 向． 否则，先利用式( 7) 求出共轭系数，再通过式( 6) 求出搜索方向． ( 5) 使 Δf = dn，通过求解灵敏度问题，即式( 3) 得 出温度增量 ΔT． ( 6) 利用式( 9) 算出搜索步长，再利用式( 8) 计算 得到新的形状函数 f( y) ． 之后，利用式( 1) 求出新的 Ti ( xi，yi ) ，回到第二步． 2 模型验证及仿真结果分析 2. 1 模型验证条件 本文以高炉炉衬侵蚀为研究背景，取二维炉衬为 研究区域，炉衬内壁面温度 T0 = 1150 ℃，即高炉炉渣 熔点值． 炉衬初期没有侵蚀时，x 轴方向上的厚度 b = 0. 5 m，y 轴方向上的长度 L = 1 m． 其测点均匀布置在 外壁面处，个数为 M． 此研究区域二维尺寸为 0. 5 m × 1 m，网格数选取 50 × 100，网格大小为 1 cm，这就意味 着形状函数最小变化量为 1 cm． 网格数的选取是综合 考虑了形状函数的变化量和程序运行时间得出的． 计 算中热流密度 q0取值 50000 W·m － 2，初始曲线全部采 用的是形状函数 x = f( y) = 0． 2. 2 计算结果分析 ( 1) 第一种真实形状函数方程: x = f( y) = 0. 35 － 0. 1sin ( 2πy) ，0≤y≤1． ( 10) 在真实的形状函数下，右侧固定位置上的热电偶 实测值见表 1． 利用这 11 个测点的实测值，可以通过共轭梯度 法，把求解边界形状反问题转化成不断求解正问题、灵 敏度问题及伴随问题的循环，在满足停止收敛标准下， 得到最终的反演曲线，它与真实侵蚀边界形状函数的 比较如图2所示． 反演得到的曲线是由相邻两个反演 · 675 ·

苏福永等：基于传热反问题的高炉炉村单层侵蚀模型 577 表1第一种形状函数条件下测温点坐标及实测值 值连线而得的，相当于由11个点构成的10段连接线. Table1 Temperartures and coordinates of measurement points for shape 图2中可以看出，11个测点反推出的曲线反演值基本 function I 重合在真实的侵蚀边界上，符合度极高，同时真实值与 测点编号 测点x坐标/m测点y坐标/m温度实测值/℃ 反演值每个对应点的相对误差保持在了1%左右，平 0.5 0 230.50 均相对误差为0.83%，控制在3%以内 2 0.5 0.1 203.76 (2)第二种真实形状函数方程： 0.5 0.2 175.17 =f0)）=0.45-0.4, 0≤y<0.5: 0.5 0.3 203.29 10.25+0.4(y-0.5),0.5≤y≤1. (11) 5 0.5 0.4 308.40 此形状函数的特点是靠近上下壁面处侵蚀程度 6 0.5 0.5 487.64 强，中间侵蚀程度弱.在真实形状函数下，右侧固定位 7 0.5 0.6 707.81 置上的热电偶实测值见表2. 8 0.5 0.7 869.31 利用11个测点的实测值，通过共轭梯度法进行求 9 0.5 0.8 872.59 解，在满足停止收敛标准下，得到最终的反演曲线，它 0 0.5 0.9 748.91 与真实侵蚀边界形状函数的比较如图3所示.反演得 11 0.5 1.0 670.73 到的曲线是由相邻两个反演值连线而得的，相当于由 1.0 (a) 温度 1150 一直实值 0.8 1052 ·反演值 0.8 954.0 856.0 0.6 0.6 758.0 660.0 562.0 0.4 0.4 464.0 366.0 268.0 0.2 03 70.0 0.10.20.30.405 00.10.20.30.40.5 x/m x/m 图2第一种形状函数条件下反演结果.()形状曲线图：(b)温度云图 Fig.2 Real and calculated shapes for shape function 1:(a)shape curve:(b)temperature nephogram 表2第二种形状函数条件下测温点坐标及实测值 上，同时真实值与反演值之间的相对误差基本都控制 Table 2 Temperartures and coordinates of measurement points for shape 在了1%以内，整体的平均相对误差保持在0.73%，控 function 2 制在3%以内. 测点编号测点x坐标/m测点y坐标/m温度实测值/℃ 2.3测点数的影响 1 0.5 0 848.79 保持其他参数不变，即炉衬外壁面热流强度值9。= 0.5 0.1 723.19 50000W·m2,初始猜测曲线为x=∫(y）=0,导热率 0.5 0.2 544.80 k=10W·m·K,同样的选取上面2种具有代表性 4 0.5 0.3 388.90 的形状函数来进行计算，将测点数M从11增加到21. 5 0.5 0.4 279.94 在真实侵蚀边界曲线为第一种形状函数条件下， 6 0.5 0.5 240.15 求解正问题得到整个二维区域内的温度值，同时也得 0.5 0.6 279.94 到了测点处的实际值.利用这些测点温度值，通过共 8 0.5 0.7 388.90 轭梯度法得到最终反演侵蚀曲线，其结果见图4.如果 9 0.5 0.8 544.80 测点数为11，则反演曲线就是由11个点之间连结而 o 0.5 0.9 723.19 成的10小段直线构成，若为21个，则反演曲线由21 11 0.5 1.0 848.79 个点之间连结而成的20小段直线构成. 从图4可以得知，当测点数M=11和21时，其反 11个点构成的10段连接线.图3中可以看出，11个测 演结果基本没有差别，与真实曲线的吻合度也极高，相 点反推出的曲线反演值基本重合在真实的侵蚀边界 对平均误差分别为0.83%和0.81%，都控制在了3%

苏福永等: 基于传热反问题的高炉炉衬单层侵蚀模型 表 1 第一种形状函数条件下测温点坐标及实测值 Table 1 Temperartures and coordinates of measurement points for shape function 1 测点编号 测点 x 坐标/m 测点 y 坐标/m 温度实测值/℃ 1 0. 5 0 230. 50 2 0. 5 0. 1 203. 76 3 0. 5 0. 2 175. 17 4 0. 5 0. 3 203. 29 5 0. 5 0. 4 308. 40 6 0. 5 0. 5 487. 64 7 0. 5 0. 6 707. 81 8 0. 5 0. 7 869. 31 9 0. 5 0. 8 872. 59 10 0. 5 0. 9 748. 91 11 0. 5 1. 0 670. 73 值连线而得的，相当于由 11 个点构成的 10 段连接线． 图 2 中可以看出，11 个测点反推出的曲线反演值基本 重合在真实的侵蚀边界上，符合度极高，同时真实值与 反演值每个对应点的相对误差保持在了 1% 左右，平 均相对误差为 0. 83% ，控制在 3% 以内． ( 2) 第二种真实形状函数方程: x = f( y) = 0. 45 － 0. 4y， 0≤y ＜ 0. 5; {0. 25 + 0. 4( y － 0. 5) ， 0. 5≤y≤1． ( 11) 此形状函数的特点是靠近上下壁面处侵蚀程度 强，中间侵蚀程度弱． 在真实形状函数下，右侧固定位 置上的热电偶实测值见表 2． 利用 11 个测点的实测值，通过共轭梯度法进行求 解，在满足停止收敛标准下，得到最终的反演曲线，它 与真实侵蚀边界形状函数的比较如图 3 所示． 反演得 到的曲线是由相邻两个反演值连线而得的，相当于由 图 2 第一种形状函数条件下反演结果 . ( a) 形状曲线图; ( b) 温度云图 Fig． 2 Ｒeal and calculated shapes for shape function 1: ( a) shape curve; ( b) temperature nephogram 表 2 第二种形状函数条件下测温点坐标及实测值 Table 2 Temperartures and coordinates of measurement points for shape function 2 测点编号 测点 x 坐标/m 测点 y 坐标/m 温度实测值/℃ 1 0. 5 0 848. 79 2 0. 5 0. 1 723. 19 3 0. 5 0. 2 544. 80 4 0. 5 0. 3 388. 90 5 0. 5 0. 4 279. 94 6 0. 5 0. 5 240. 15 7 0. 5 0. 6 279. 94 8 0. 5 0. 7 388. 90 9 0. 5 0. 8 544. 80 10 0. 5 0. 9 723. 19 11 0. 5 1. 0 848. 79 11 个点构成的10 段连接线． 图3 中可以看出，11 个测 点反推出的曲线反演值基本重合在真实的侵蚀边界 上，同时真实值与反演值之间的相对误差基本都控制 在了 1% 以内，整体的平均相对误差保持在 0. 73% ，控 制在 3% 以内． 2. 3 测点数的影响 保持其他参数不变，即炉衬外壁面热流强度值 q0 = 50000 W·m － 2，初始猜测曲线为 x = f( y) = 0，导热率 k = 10 W·m － 1·K － 1，同样的选取上面 2 种具有代表性 的形状函数来进行计算，将测点数 M 从 11 增加到 21． 在真实侵蚀边界曲线为第一种形状函数条件下， 求解正问题得到整个二维区域内的温度值，同时也得 到了测点处的实际值． 利用这些测点温度值，通过共 轭梯度法得到最终反演侵蚀曲线，其结果见图 4． 如果 测点数为 11，则反演曲线就是由 11 个点之间连结而 成的 10 小段直线构成，若为 21 个，则反演曲线由 21 个点之间连结而成的 20 小段直线构成． 从图 4 可以得知，当测点数 M = 11 和 21 时，其反 演结果基本没有差别，与真实曲线的吻合度也极高，相 对平均误差分别为 0. 83% 和 0. 81% ，都控制在了 3% · 775 ·

·578· 工程科学学报，第39卷，第4期 1.0 1.0m (a) b 温度℃ 1150 一真实值 0.8 1059 ·反演值 0.8 968.0 877.0 786.0 0.6 0.6 695.0 604.0 0.4 0.4 513.0 422.0 331.0 02 0.2 240.0 0.10.20.3040.5 000.10.2,0.30405 x/m x/m 图3第二种形状函数条件下反演结果.()形状曲线图：(b)温度云图 Fig.3 Real and calculated shapes for shape function 2:(a)shape curve:(b)temperature nephogram 1.0 1.0 一真实值 一真实值 0.8 M=11 0.8 ■M=11 M=21 4M-21 0.6 0.6 0.4 0.4 0.2 0.2 0.10.20.30.40.5 00 0.1 020.30.405 x/m x/m 图4第一种形状函数条件下测点数M对反演结果的影响 图5第二种形状函数条件下测点数M对反演结果的影响 Fig.4 Real and calculated shapes for shape function 1 (M=11 Fig.5 Real and calculated shapes for shape function 2 (M 11 and M 21) and M 21) 以内. 式中：9o为炉衬外壁面热流强度值，W·m2:k为材料热 在真实侵蚀边界曲线为第二种形状函数条件下， 导率，W·m·K:T。为内壁表面温度，K:T为测点i 其反演侵蚀曲线结果见图5.可以看出，当测点数M= 处的实际测得的温度值，K:b为炉衬初期没有侵蚀时x 11和21时，其反演结果也基本差别不大，与真实曲线 轴方向上的厚度 的吻合度也极高，两者相对平均误差分别为0.73%和 在真实侵蚀边界曲线为第一种形状函数条件下， 0.66%,都控制在了3%以内. 求解正问题得到了测点处的真实值，即所谓的实际测 细细比较，测点数M从11变为21时，其反演结果 得温度值.之后分别选取不同的初始猜测曲线，通过 的准确性略有提升，这是因为测点数目的增加在一定 共轭梯度法，获得了最终的反演侵蚀曲线，其结果见图 程度上更能捕捉曲线的特征. 6.图6中，反演值1是在初始猜测曲线为x=f(y)=0 2.4初始猜测曲线的影响 情况下得到的最终结果，反演值2的初始猜测曲线为 保持其他参数都不变，即炉衬外壁面热流强度值 式(12).可以看出，两者的反演结果基本没有差别，与 6=50000W·m2,测点数M=11,导热率k=10W· 真实曲线的吻合度极高，相对平均误差分别为0.83% m·K改变初始的猜测曲线，反演值1的初始猜测 和0.66%，所有点都控制在了3%以内. 曲线取的是x=∫(y)=0,而反演值2的初始猜测曲线 在真实侵蚀边界曲线为第二种形状函数条件下， 函数如下： 求解正问题得到了测点处的真实值，即所谓的实际测 本=/八)=6-k-I (12) 得温度值.之后分别选取不同的初始猜测曲线，通过 9 共轭梯度法，反推出了最终的侵蚀曲线，其结果见图

工程科学学报，第 39 卷，第 4 期 图 3 第二种形状函数条件下反演结果． ( a) 形状曲线图; ( b) 温度云图 Fig． 3 Ｒeal and calculated shapes for shape function 2: ( a) shape curve; ( b) temperature nephogram 图 4 第一种形状函数条件下测点数 M 对反演结果的影响 Fig． 4 Ｒeal and calculated shapes for shape function 1 ( M = 11 and M = 21) 以内． 在真实侵蚀边界曲线为第二种形状函数条件下， 其反演侵蚀曲线结果见图 5． 可以看出，当测点数 M = 11 和 21 时，其反演结果也基本差别不大，与真实曲线 的吻合度也极高，两者相对平均误差分别为 0. 73% 和 0. 66% ，都控制在了 3% 以内． 细细比较，测点数 M 从11 变为21 时，其反演结果 的准确性略有提升，这是因为测点数目的增加在一定 程度上更能捕捉曲线的特征． 2. 4 初始猜测曲线的影响 保持其他参数都不变，即炉衬外壁面热流强度值 q0 = 50000 W·m － 2，测点数 M = 11，导热率 k = 10 W· m － 1·K － 1 ． 改变初始的猜测曲线，反演值 1 的初始猜测 曲线取的是 x = f( y) = 0，而反演值 2 的初始猜测曲线 函数如下: x = f( y) = b － k T0 － Tc i q0 ． ( 12) 图 5 第二种形状函数条件下测点数 M 对反演结果的影响 Fig． 5 Ｒeal and calculated shapes for shape function 2 ( M = 11 and M = 21) 式中: q0为炉衬外壁面热流强度值，W·m － 2 ; k 为材料热 导率，W·m － 1·K － 1 ; T0为内壁表面温度，K; Tc i 为测点 i 处的实际测得的温度值，K; b 为炉衬初期没有侵蚀时 x 轴方向上的厚度． 在真实侵蚀边界曲线为第一种形状函数条件下， 求解正问题得到了测点处的真实值，即所谓的实际测 得温度值． 之后分别选取不同的初始猜测曲线，通过 共轭梯度法，获得了最终的反演侵蚀曲线，其结果见图 6． 图 6 中，反演值 1 是在初始猜测曲线为 x = f( y) = 0 情况下得到的最终结果，反演值 2 的初始猜测曲线为 式( 12) ． 可以看出，两者的反演结果基本没有差别，与 真实曲线的吻合度极高，相对平均误差分别为 0. 83% 和 0. 66% ，所有点都控制在了 3% 以内． 在真实侵蚀边界曲线为第二种形状函数条件下， 求解正问题得到了测点处的真实值，即所谓的实际测 得温度值． 之后分别选取不同的初始猜测曲线，通过 共轭梯度法，反推出了最终的侵蚀曲线，其结果见图 · 875 ·