CN102876825A - 一种高炉炉缸侵蚀线位置的判定方法 - Google Patents
一种高炉炉缸侵蚀线位置的判定方法 Download PDFInfo
- Publication number
- CN102876825A CN102876825A CN2012103362679A CN201210336267A CN102876825A CN 102876825 A CN102876825 A CN 102876825A CN 2012103362679 A CN2012103362679 A CN 2012103362679A CN 201210336267 A CN201210336267 A CN 201210336267A CN 102876825 A CN102876825 A CN 102876825A
- Authority
- CN
- China
- Prior art keywords
- temperature
- blast furnace
- cupola well
- grid
- hearth
- 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.)
- Pending
Links
Images
Landscapes
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
一种高炉炉缸侵蚀线位置的判定方法，属于高炉长寿技术领域。该方法提供一种比现有方法计算量更小且精度高的确定高炉炉缸侵蚀线位置的方法。该方法依据稳态导热方程，直接用稳定网格的有限元或有限差分的方法计算炉缸砖衬中的温度场，在迭代计算中根据网格节点的温度确定砖衬导热系数，如果温度高于1450℃，将采用铁水的有效导热系数，计算确定炉缸砖衬中的温度场，根据温度场确定1150℃侵蚀线的位置，这种确定侵蚀线位置的方法只需计算一次温度场，因此较之于以往的方法，计算量大幅缩小，计算周期短，能满足高炉操作的各种要求，计算稳定性大幅提高，且将铁水与砖衬的分界温度设为1450℃，更符合高炉实际。
Description
所属技术领域
本发明属于高炉长寿技术领域，用于确定冶金行业炼铁高炉炉缸中1150℃侵蚀线的位置，为高炉操作提供基础信息。
背景技术
高炉长寿和高效稳定生产的关键之一是高炉炉缸状况能得到及时、准确的监控，为此国内外开发了大量高炉炉缸监控系统，主要是监控炉缸中1150℃侵蚀线的位置。主流方式是依据炉缸炉衬中的热电偶温度信息通过传热模型的计算确定1150℃侵蚀线的位置，其中，主要的确定方式是利用边界元和有限元的方式逆解算，即利用非线性优化的方式假定1150℃侵蚀线的位置，然后以侵蚀线位置为边界条件用边界元或有限元的方法计算温度场，然后得到热电偶处的温度计算值，温度计算值和实测值方差最小的1150℃侵蚀线的假定位置为最终得到的侵蚀线位置。显然，该方式需要多次试探侵蚀线的位置，并计算温度场，计算量将很大，导致计算时间长、计算网格不稳定并且计算精度却无法得到保证。
针对这个问题，已有采用稳定网格正解算方式确定1150℃侵蚀线的位置。即直接用原始砖衬的位置数据设置网格用有限元或有限差分的方法计算炉缸砖衬中的温度场，根据温度场确定1150℃侵蚀线的位置，温度高于1150℃的区域视为铁水区，温度低于1150℃的区域视为砖衬区，在铁水区将导热系数换为铁水的导热系数，然后再次计算温度场并重复上述过程，直至温度场稳定下来，此时根据温度场确定的1150℃侵蚀线的位置为最终得到的侵蚀线位置。这种方式计算网格稳定且容易设定，计算精度有保证，计算量也大幅减小，但同样需要重复计算温度场，导致计算量仍然不小，且将1150℃视为铁水区和砖衬的分界点不尽合理，因为铁水导热率大大高于砖衬的导热率，实际高炉中砖衬与铁水接触处温度应接近铁水的主体温度。
发明内容
本发明所要解决的技术问题是：提供一种较之于现有技术判定高炉炉缸侵蚀线位置计算量更小且精度高的方法，确定高炉炉缸侵蚀线位置，为高炉操作提供稳定和计算周期短的炉缸基础信息。
本发明解决其技术问题采用的技术方案是：一种高炉炉缸侵蚀线位置的判定方法，其特征在于：依据稳态导热方程，采用稳定网格的有限元或有限差分法计算炉缸砖衬中的温度场，在迭代计算中根据网格节点的温度确定导热系数，进而得到炉缸砖衬中的温度场，再根据温度场判定1150℃侵蚀线的位置，具体步骤包括：
a、将炉缸分为炉缸侧壁区、炉底中心区和炉缸拐角区三部分，其中所述炉缸侧壁区的传热方向垂直于炉壁，所述炉底中心区的传热方向垂直于炉底，将三部分区域分别定义为长方形，各自都包含铁水和砖衬；
b、读入高炉炉缸中预设的热电偶温度，砖衬及铁水的导热系数及各区域的网格位置数据；
c、分区域进行稳态导热计算，在炉缸侧壁区和炉底中心区为由热电偶温度递推区域内网格上各节点的温度的一维计算，所述炉缸侧壁区的递推方向为由外而内，所述炉底中心区的递推方向为自下而上，在炉缸拐角区以炉缸侧壁区和炉底中心区所得到的边界温度作为边界条件，采用有限元或有限差分法二维计算本区域内网格上各节点的温度；
d、将分区域的计算结果合并构成炉缸的温度场，并插值计算得到高炉炉缸1150℃侵蚀线的位置。
进一步地，步骤c中各区域网格上的导热系数根据温度选定，当温度高于1450℃时选用铁水的导热系数，否则选用对应温度砖衬的导热系数。
进一步地，步骤b后还包括选定并去除高炉已损坏部分的热电偶温度。
进一步地，所述铁水与砖衬的分界温度为1450℃。
本发明判定方法的研制和应用，使得确定高炉炉缸侵蚀线位置只需计算一次温度场，因此和已有技术方法相比，计算量大幅缩小，计算周期短，能满足高炉操作的各种要求，计算稳定性大幅提高，可在一台电脑上实现监控多个高炉的炉缸状况；且将铁水与砖衬的分界温度设为1450℃，更符合高炉实际。
附图说明
图1是本发明判定方法针对炉缸分区计算的分区示意图。
图2是本发明判定方法的流程示意图。
具体实施方式
为克服现有技术在高炉长寿监控方面的不足，概括来看：本发明依据稳态导热方程，直接用稳定网格的有限元或有限差分的方法计算炉缸砖衬中的温度场，在迭代计算中根据网格节点的温度确定砖衬导热系数，如果温度高于1450℃，将采用铁水的有效导热系数，计算确定炉缸砖衬中的温度场，根据温度场确定1150℃侵蚀线的位置，这种确定侵蚀线位置的方法只需计算一次温度场，因此和以前的方法比，计算量大幅缩小，计算周期短，能满足高炉操作的各种要求，计算稳定性大幅提高，且将铁水与砖衬的分界温度设为1450℃，更符合高炉实际。
以下便结合附图，对本发明判定方法做如下详细说明，以使本领域技术人员能充分理解本发明创新实质及其应用价值。如图1所示的是该判定方法针对炉缸分区计算的分区示意图。从图1中可以清楚地看到，作为判定方法的开端步骤，需要对炉缸作划分，主要分为炉缸侧壁区1、炉底中心区2和炉缸拐角区3三部分。其中炉缸侧壁区1的传热方向6垂直于炉壁，炉底中心区2的传热方向7垂直于炉底，而炉缸拐角区3的传热方向无法直接以比照的方式简单确定，只可在炉缸侧壁区1和炉底中心区2各自的传热方向之间，选择一个波动范围。本发明中，将三部分区域分别定义为长方形，各自都包含铁水4和砖衬5。各区域中根据热电偶的布置安排计算网格，在计算中网格不再改变，因此称为稳定网格。由此便完成了温度场确定的基础。如图1箭头6、7所示，对于炉缸侧壁区1和炉底中心区2可采用稳定网格递推计算的一维计算方式，递推方向由下往上或由外向内。而对于相对难以确定的炉缸拐角区3，以炉缸侧壁区1和炉底中心区2所得到的边界温度作为边界条件，也采用稳定网格的有限元或有限差分的方法进行二维计算，得到炉缸拐角区砖衬的温度变化呈现。
上述分区域计算中，网格上导热系数依据温度选定对应温度的砖衬导热系数，如果温度高于1450℃，将采用铁水的有效导热系数。
最后，将三个区域的温度场合并为高炉炉缸的温度场，并依据炉缸温度场插值计算确定1150℃侵蚀线的位置。
上述该种判定方法，根据高炉炉缸处热电偶温度信息，可很快计算出1150℃侵蚀线的位置，计算方法简便，计算精度有保证，可在一台电脑上监控多个高炉的炉缸状况。
从明确的判定实施步骤来看：Step1：读入高炉热电偶温度、砖衬及铁水的导热系数及各区域的网格位置数据等。
Step2：剔除不合理的温度数据。由于高炉内热电偶经常损坏，需要对热电偶的温度数据进行整理，剔除掉不合理的温度数据。判定规则为：1、接近热面的热电偶温度应高于接近冷面的热电偶温度。2、热电偶的最大温度不应超过铁水温度。
Step3：分区域进行稳态导热计算，其中，炉缸侧壁区1和炉底中心区2为一维传热计算；炉缸拐角区3为二维传热计算。采用有限元或有限差分法计算区域内网格各节点的温度。在上述迭代计算中，递推方向均为从下往上，从外往内，根据网格节点温度确定砖衬耐材的导热系数，如果温度高于1450℃，将采用铁水的有效导热系数。另外，炉缸侧壁区1和炉底中心区2得到的边界温度，将作为炉缸拐角区3的边界条件，用于二维传热计算。计算得到温度场数据。三个区域传热控制方程分别为：
（1）‘炉缸侧壁区1’的传热控制方程
（1）
式中，q r 为炉缸区耐材中的传热量。
（2）‘炉底中心区2’的传热控制方程
式中，qZ为炉缸区耐材中的传热量。
（3）‘炉缸拐角区3’的传热控制方程
（3）
式中，λ为耐材的导热系数。
Step4：将炉缸侧壁区1、炉底中心区2和炉缸拐角区3三个区域的温度场计算结果合并为炉缸区的温度场。温度场数据包括：各个风口截面的温度场，炉底炉缸各层耐火材料的温度场。
Step5：由炉缸区的温度场插值计算得到高炉炉缸1150℃侵蚀线的位置。
需要特别重视的是：高炉炉缸热电偶工作环境很恶劣，损坏的很多，为保证计算结果的正确性，必须剔除不合理的数据，以提高侵蚀线判定的准确度。
综上，本发明判定方法的研制和应用，使得确定高炉炉缸侵蚀线位置只需计算一次温度场。因此和已有技术方法相比，计算量大幅缩小，计算周期短，能满足高炉操作的各种要求，计算稳定性大幅提高，可在一台电脑上实现监控多个高炉的炉缸状况；且将铁水与砖衬的分界温度设为1450℃，更符合高炉实际。
然而上述有关实施例的详细描述旨在便于本发明技术方案的理解，并非以此限制本发明专利申请的保护范围，但凡在上述实施例基础上经等效替换所形成的技术方案，均应视作为本发明权利要求书所要求保护的范畴。
Claims (4)
1.一种高炉炉缸侵蚀线位置的判定方法，其特征在于：依据稳态导热方程，采用稳定网格的有限元或有限差分法计算炉缸砖衬中的温度场，在迭代计算中根据网格节点的温度确定导热系数，进而得到炉缸砖衬中的温度场，再根据温度场判定1150℃侵蚀线的位置，具体步骤包括：
a、将炉缸分为炉缸侧壁区、炉底中心区和炉缸拐角区三部分，其中所述炉缸侧壁区的传热方向垂直于炉壁，所述炉底中心区的传热方向垂直于炉底，将三部分区域分别定义为长方形，各自都包含铁水和砖衬；
b、读入高炉炉缸中预设的热电偶温度，砖衬及铁水的导热系数及各区域的网格位置数据；
c、分区域进行稳态导热计算，在炉缸侧壁区和炉底中心区为由热电偶温度递推区域内网格上各节点的温度的一维计算，所述炉缸侧壁区的递推方向为由外而内，所述炉底中心区的递推方向为自下而上，在炉缸拐角区以炉缸侧壁区和炉底中心区所得到的边界温度作为边界条件，采用有限元或有限差分法二维计算本区域内网格上各节点的温度；
d、将分区域的计算结果合并构成炉缸的温度场，并插值计算得到高炉炉缸1150℃侵蚀线的位置。
2.根据权利要求1所述一种高炉炉缸侵蚀线位置的判定方法，其特征在于：步骤c中，各区域网格上的导热系数根据温度选定，当温度高于1450℃时选用铁水的导热系数，否则选用对应温度砖衬的导热系数。
3.根据权利要求1所述的一种高炉炉缸侵蚀线位置的判定方法，其特征在于：步骤b后还包括选定并去除高炉已损坏部分的热电偶温度。
4.根据权利要求1所述一种高炉炉缸侵蚀线位置的判定方法，其特征是：所述铁水与砖衬的分界温度为1450℃。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2012103362679A CN102876825A (zh) | 2012-09-13 | 2012-09-13 | 一种高炉炉缸侵蚀线位置的判定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2012103362679A CN102876825A (zh) | 2012-09-13 | 2012-09-13 | 一种高炉炉缸侵蚀线位置的判定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102876825A true CN102876825A (zh) | 2013-01-16 |
Family
ID=47478335
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2012103362679A Pending CN102876825A (zh) | 2012-09-13 | 2012-09-13 | 一种高炉炉缸侵蚀线位置的判定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102876825A (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104313224A (zh) * | 2014-10-16 | 2015-01-28 | 攀钢集团西昌钢钒有限公司 | 一种无陶瓷杯钒钛磁铁矿高炉的炉缸浸蚀结厚报警方法 |
| CN104460574A (zh) * | 2014-04-10 | 2015-03-25 | 江苏永钢集团有限公司 | 一种可实时立体监视高炉状态的系统和方法 |
| CN105005632A (zh) * | 2014-04-17 | 2015-10-28 | 宝山钢铁股份有限公司 | 多层耐火砖炉墙结构的高炉炉缸侵蚀预测方法 |
| CN105241576A (zh) * | 2015-09-10 | 2016-01-13 | 马鞍山市安工大工业技术研究院有限公司 | 一种基于分布式光纤的高炉热风炉内衬侵蚀建模系统及方法 |
| CN108090293A (zh) * | 2017-12-26 | 2018-05-29 | 国家超级计算天津中心 | 一种高炉炉缸炉底侵蚀包络面确定方法 |
| CN110688376A (zh) * | 2019-09-27 | 2020-01-14 | 中冶赛迪重庆信息技术有限公司 | 一种温度数据清洗方法、系统及设备 |
| CN110826172A (zh) * | 2019-09-23 | 2020-02-21 | 南京航空航天大学 | 一种高炉炉缸横截面侵蚀边界二维快速计算方法 |
| CN111737899A (zh) * | 2020-06-22 | 2020-10-02 | 中冶赛迪重庆信息技术有限公司 | 一种面向高炉炉缸长寿的热应力计算方法及系统 |
| CN113971359A (zh) * | 2021-10-28 | 2022-01-25 | 中冶南方工程技术有限公司 | 一种炉缸侵蚀的识别方法、电子设备以及存储介质 |
| CN114139412A (zh) * | 2021-10-28 | 2022-03-04 | 中冶南方工程技术有限公司 | 一种炉缸侵蚀的评估方法、电子设备以及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR940008056B1 (ko) * | 1991-12-30 | 1994-09-01 | 포항종합제철 주식회사 | 고로노저 연와 잔존 두께 및 응고층 두께 추정방법 |
| JPH07278632A (ja) * | 1994-04-11 | 1995-10-24 | Nippon Steel Corp | 高炉炉底状況の推定方法 |
| JPH0978113A (ja) * | 1995-09-18 | 1997-03-25 | Sumitomo Metal Ind Ltd | 高炉炉底構造の設計方法 |
-
2012
- 2012-09-13 CN CN2012103362679A patent/CN102876825A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR940008056B1 (ko) * | 1991-12-30 | 1994-09-01 | 포항종합제철 주식회사 | 고로노저 연와 잔존 두께 및 응고층 두께 추정방법 |
| JPH07278632A (ja) * | 1994-04-11 | 1995-10-24 | Nippon Steel Corp | 高炉炉底状況の推定方法 |
| JPH0978113A (ja) * | 1995-09-18 | 1997-03-25 | Sumitomo Metal Ind Ltd | 高炉炉底構造の設計方法 |
Non-Patent Citations (3)
| Title |
|---|
| 吴俐俊: "基于边界元法的高炉炉底炉缸侵蚀模型", 《上海交通大学学报》 * |
| 李家新等: "高炉炉底侵蚀线的计算", 《包头钢铁学院学报》 * |
| 陈良玉等: "高炉炉缸炉底内衬侵蚀形貌的计算和应用", 《中国科技论文在线》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104460574A (zh) * | 2014-04-10 | 2015-03-25 | 江苏永钢集团有限公司 | 一种可实时立体监视高炉状态的系统和方法 |
| CN105005632A (zh) * | 2014-04-17 | 2015-10-28 | 宝山钢铁股份有限公司 | 多层耐火砖炉墙结构的高炉炉缸侵蚀预测方法 |
| CN105005632B (zh) * | 2014-04-17 | 2017-11-28 | 宝山钢铁股份有限公司 | 多层耐火砖炉墙结构的高炉炉缸侵蚀预测方法 |
| CN104313224B (zh) * | 2014-10-16 | 2016-05-11 | 攀钢集团西昌钢钒有限公司 | 一种无陶瓷杯钒钛磁铁矿高炉的炉缸浸蚀结厚报警方法 |
| CN104313224A (zh) * | 2014-10-16 | 2015-01-28 | 攀钢集团西昌钢钒有限公司 | 一种无陶瓷杯钒钛磁铁矿高炉的炉缸浸蚀结厚报警方法 |
| CN105241576A (zh) * | 2015-09-10 | 2016-01-13 | 马鞍山市安工大工业技术研究院有限公司 | 一种基于分布式光纤的高炉热风炉内衬侵蚀建模系统及方法 |
| CN105241576B (zh) * | 2015-09-10 | 2017-12-15 | 马鞍山市安工大工业技术研究院有限公司 | 一种基于分布式光纤的高炉热风炉内衬侵蚀建模方法 |
| CN108090293B (zh) * | 2017-12-26 | 2021-04-30 | 国家超级计算天津中心 | 一种高炉炉缸炉底侵蚀包络面确定方法 |
| CN108090293A (zh) * | 2017-12-26 | 2018-05-29 | 国家超级计算天津中心 | 一种高炉炉缸炉底侵蚀包络面确定方法 |
| CN110826172A (zh) * | 2019-09-23 | 2020-02-21 | 南京航空航天大学 | 一种高炉炉缸横截面侵蚀边界二维快速计算方法 |
| CN110688376A (zh) * | 2019-09-27 | 2020-01-14 | 中冶赛迪重庆信息技术有限公司 | 一种温度数据清洗方法、系统及设备 |
| CN111737899A (zh) * | 2020-06-22 | 2020-10-02 | 中冶赛迪重庆信息技术有限公司 | 一种面向高炉炉缸长寿的热应力计算方法及系统 |
| CN111737899B (zh) * | 2020-06-22 | 2022-05-31 | 中冶赛迪重庆信息技术有限公司 | 一种面向高炉炉缸长寿的热应力计算方法及系统 |
| CN113971359A (zh) * | 2021-10-28 | 2022-01-25 | 中冶南方工程技术有限公司 | 一种炉缸侵蚀的识别方法、电子设备以及存储介质 |
| CN114139412A (zh) * | 2021-10-28 | 2022-03-04 | 中冶南方工程技术有限公司 | 一种炉缸侵蚀的评估方法、电子设备以及存储介质 |
| WO2023070693A1 (zh) * | 2021-10-28 | 2023-05-04 | 中冶南方工程技术有限公司 | 一种炉缸侵蚀的评估方法、电子设备以及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102876825A (zh) | 一种高炉炉缸侵蚀线位置的判定方法 | |
| Xu et al. | A new type of two-supply, one-return, triple pipe-structured heat loss model based on a low temperature district heating system | |
| CN103866061B (zh) | 一种高炉炉缸侵蚀偏差监测方法 | |
| CN102095516B (zh) | 一种废杂铜冶炼炉温测量方法 | |
| CN106916915B (zh) | 一种在线监测高炉软熔带根部位置的方法 | |
| CN105005632A (zh) | 多层耐火砖炉墙结构的高炉炉缸侵蚀预测方法 | |
| CN106874648B (zh) | 一种高炉高热负荷区域操作炉型计算方法 | |
| CN112528572B (zh) | 一种基于overset模型的低温碳化炉丝束加热过程三维模拟方法 | |
| CN108090293B (zh) | 一种高炉炉缸炉底侵蚀包络面确定方法 | |
| CN105087882A (zh) | 一种立式淬火炉热处理阶段划分方法 | |
| CN110781566A (zh) | 一种炉缸凝铁层计算方法、系统、存储介质及电子终端 | |
| Ding et al. | Multi‐Objective Optimization of Slab Heating Process in Walking Beam Reheating Furnace Based on Particle Swarm Optimization Algorithm | |
| KR102531803B1 (ko) | 용광로의 내화물 라이닝의 마모를 모니터링하기 위한 방법 | |
| CN104611483A (zh) | 一种基于冷却壁热负荷和炉衬测厚的炉温在线检测方法 | |
| Li et al. | 3D simulation of borosilicate glass all‐electric melting furnaces | |
| Tang et al. | Modeling of steel slab reheating process in a walking beam reheating furnace | |
| CN103088176A (zh) | 高炉炉壳气隙的侦测方法 | |
| CN104313224B (zh) | 一种无陶瓷杯钒钛磁铁矿高炉的炉缸浸蚀结厚报警方法 | |
| Ahmed et al. | A Truncated transient slab model for a reheating furnace | |
| JPH0978113A (ja) | 高炉炉底構造の設計方法 | |
| CN108917987A (zh) | 混凝土内部温度预测方法、系统与装置 | |
| Feng et al. | A new method to optimize furnace designs using daily flow rates to maximize energy savings in the steady production | |
| Dong et al. | Mathematical model for fast computation of erosion profile in submerged arc furnace with freeze lining | |
| Guojun et al. | Three-dimensional zone method mathematic model on conventional reheating furnace | |
| Zhang et al. | The study of heat-mass transfer characteristics and optimization on electric arc furnace with the novel fence electrode |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130116 |