CN105930666B - 一种气流床气化炉炉膛温度的在线软测量方法 - Google Patents
一种气流床气化炉炉膛温度的在线软测量方法 Download PDFInfo
- Publication number
- CN105930666B CN105930666B CN201610267960.3A CN201610267960A CN105930666B CN 105930666 B CN105930666 B CN 105930666B CN 201610267960 A CN201610267960 A CN 201610267960A CN 105930666 B CN105930666 B CN 105930666B
- Authority
- CN
- China
- Prior art keywords
- represent
- gasification
- vaporizer
- furnace
- temperature
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000002309 gasification Methods 0.000 claims abstract description 124
- 238000006243 chemical reaction Methods 0.000 claims abstract description 56
- 239000003245 coal Substances 0.000 claims abstract description 39
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000013461 design Methods 0.000 claims abstract description 5
- 239000006200 vaporizer Substances 0.000 claims description 80
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 45
- 239000000126 substance Substances 0.000 claims description 44
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 34
- 238000002485 combustion reaction Methods 0.000 claims description 33
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 32
- 238000000197 pyrolysis Methods 0.000 claims description 31
- 239000007789 gas Substances 0.000 claims description 27
- 238000012546 transfer Methods 0.000 claims description 24
- 239000002817 coal dust Substances 0.000 claims description 22
- 239000002245 particle Substances 0.000 claims description 22
- 229910052757 nitrogen Inorganic materials 0.000 claims description 21
- 239000007787 solid Substances 0.000 claims description 20
- 229960004424 carbon dioxide Drugs 0.000 claims description 19
- 239000007790 solid phase Substances 0.000 claims description 19
- 239000011280 coal tar Substances 0.000 claims description 18
- 229910052739 hydrogen Inorganic materials 0.000 claims description 17
- 239000001257 hydrogen Substances 0.000 claims description 17
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims description 17
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 16
- 239000012071 phase Substances 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 229910001868 water Inorganic materials 0.000 claims description 16
- 239000001569 carbon dioxide Substances 0.000 claims description 15
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 14
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 14
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 14
- 239000001301 oxygen Substances 0.000 claims description 14
- 229910052760 oxygen Inorganic materials 0.000 claims description 14
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 claims description 13
- 229910000037 hydrogen sulfide Inorganic materials 0.000 claims description 13
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 6
- 239000008187 granular material Substances 0.000 claims description 5
- 229910002090 carbon oxide Inorganic materials 0.000 claims description 4
- 238000001816 cooling Methods 0.000 claims description 4
- 239000004408 titanium dioxide Substances 0.000 claims description 4
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000001754 furnace pyrolysis Methods 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 239000012159 carrier gas Substances 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 claims description 2
- 238000010025 steaming Methods 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 7
- 238000013178 mathematical model Methods 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 4
- 238000009434 installation Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000013480 data collection Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 239000002956 ash Substances 0.000 description 16
- 238000010791 quenching Methods 0.000 description 9
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000003786 synthesis reaction Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- -1 steam Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000001035 drying Methods 0.000 description 3
- 239000000428 dust Substances 0.000 description 3
- 239000002893 slag Substances 0.000 description 3
- 241000790917 Dioxys <bee> Species 0.000 description 2
- 238000000227 grinding Methods 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 238000006902 nitrogenation reaction Methods 0.000 description 2
- 238000004073 vulcanization Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 125000003636 chemical group Chemical group 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000003250 coal slurry Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000010881 fly ash Substances 0.000 description 1
- 150000004820 halides Chemical class 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C10—PETROLEUM, GAS OR COKE INDUSTRIES; TECHNICAL GASES CONTAINING CARBON MONOXIDE; FUELS; LUBRICANTS; PEAT
- C10J—PRODUCTION OF PRODUCER GAS, WATER-GAS, SYNTHESIS GAS FROM SOLID CARBONACEOUS MATERIAL, OR MIXTURES CONTAINING THESE GASES; CARBURETTING AIR OR OTHER GASES
- C10J3/00—Production of combustible gases containing carbon monoxide from solid carbonaceous fuels
- C10J3/72—Other features
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Organic Chemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Industrial Gases (AREA)
Abstract
一种气流床气化炉炉膛温度的在线软测量方法，涉及气化炉运行参数检测。包括数据采集和计算输出，在接受采集的数据后，通过相应的数学计算式及守恒方程，计算并输出炉内的温度数据给现场操作人员，为其判断气化炉内气化过程及装置的运行状态提供参考。采用数学机理建模的方法建立气化炉的动态数学模型，根据数据采集系统接受煤质数据、装置参数及运行数据，通过在线求解数学模型的方法获得Shell气化炉炉内温度的软测量值。可以为操作人员提供及时与有效的炉内温度的信息，指导其更好得完成气化反应过程的优化与运行操作，克服由于Shell气化装置自身设计缺陷对生产过程及产品质量的影响。
Description
技术领域
本发明涉及气化炉运行参数检测，尤其是涉及根据工业现场采集的运行数据通过事先建立好的模型在线求解炉内重要变量的一种气流床气化炉炉膛温度的在线软测量方法。
背景技术
Shell气流床煤气化技术由荷兰Shell国际石油公司于20世纪50年代开始研发，目前已在德国、荷兰、意大利等多个国家建成并投入使用，是已经商业化的第二代煤气化技术中最有竞争力的技术之一。
Shell煤气化技术的工艺流程主要包括磨煤干燥、加压输送、气化、除渣、除灰、湿洗等单元，煤场的原料煤首先由带式输送机送入磨煤与干燥单元中研磨和干燥，之后依次通过煤粉加压与输送单元内的常压煤粉仓、加压煤粉仓及给料仓，由高压氮气或二氧化碳将煤粉输送至气化炉煤烧嘴，在气化炉内与来自空分的高压氧气发生一系列反应，产生合成气、飞灰与炉渣。其中，气化产物中的炉渣从气化炉底部以熔渣的形式流入渣池，在渣池中被激冷分散后排出气化系统；气化产物中的合成气与飞灰由气化炉顶部送出，在激冷段内被激冷压缩机送来的激冷气激冷至800～900℃，然后经合成气冷却器(废锅)进一步冷却后送入干法除灰单元；合成气中夹带的熔融飞灰在之前的激冷段中由于温度的陡然降低会发生固化，这些固化的飞灰在除灰系统中大部分被高温高压过滤器分离出合成气并送入飞灰贮罐；经过干法除灰的合成气被分为两股，一股回到激冷压缩机作为激冷气，另一股进入湿洗单元，待合成气中的卤化物被除去且飞灰含量进一步降低(<1mg/m3)之后再被分为两股，一股送回激冷压缩机作激冷气，另一股作为最后的粗合成气给下游工段使用。
在目前Shell气化炉工业运行的过程中，由于炉内状况复杂、炉温较高等特点，导致气化反应器内包括温度在内的众多变量无法直接测量，工业过程中的实际测量装置均处于气化工艺中低温段，为此操作人员无法确定炉内的反应状况与反应温度的实时数据，致使气化炉无法达到精细优化操作，实际表现为其开车运行周期短、故障频发等现象。鉴于此种情况，在现有的生产条件下以低成本获得Shell炉内温度的实时数据就成为工业过程中亟待解决的技术问题。
目前，已有学者提出了与Shell气化炉同属气流床气化技术类别的多喷嘴对置式水煤浆气化炉炉膛温度的软测量方法(参见中国专利申请号：CN201110001570.9)，该方法基于BP人工神经网络建立炉内模型，实现炉膛内温度信息的在线预测。
发明内容
本发明的目的在于，提供无需另外增加气化炉测点，在推算气化炉炉炉膛温度的同时考虑到了炉内的化学反应和传热特性，使软测量结果能更实时准确地反映炉膛内温度的真实值，帮助现场操作人员控制气化过程的一种气流床气化炉炉膛温度的在线软测量方法。
本发明包括以下步骤：
1)数据采集；
2)计算输出；
在接受采集的数据后，通过相应的数学计算式及守恒方程，计算并输出炉内的温度数据给现场操作人员，为其判断气化炉内气化过程及装置的运行状态提供参考。
在步骤1)中，所述数据包括入炉煤质数据、气化炉结构参数和气化炉运行实时数据等；所述入炉煤质数据可由相关人员通过煤质分析获得；所述气化炉结构参数包括气化炉高度、气化炉直径、气化炉水冷壁厚度等，所述气化炉结构参数可通过气化炉设计图纸获得；所述气化炉运行实时数据可通过工厂DCS系统采集，所述气化炉运行实时数据包括进煤量、进氧量、蒸汽量、载气(二氧化碳或氮气)量以及气化炉压力。
在步骤2)中，所述计算并输出炉内的温度数据包括炉内的化学反应和炉内向外传热的计算，以及炉温的输出整个软测量，具体方法如下：
(1)根据计算精度与计算速度的要求确定气化室个数；
(2)接收数据采集部分传来的本时刻气化炉运行数据；
(3)根据进料参数计算本时刻进入气化炉热解燃烧室内的煤粉发生热解及燃烧反应的质量与能量平衡，其求解过程按下述表达式描述：
(Cg,i,0,Cp,j,0,T0,u0,ρN0,dp,0,εp,0,εg,0,P)
＝f(Tfeeding,x,mfeeding,x,dp,feeding,ρcoal)
式中Cg,i,0代表热解燃烧室内气相组分i的摩尔浓度，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，Cp,j,0代表热解燃烧室内固相组分j的质量浓度，其中j包括煤焦与灰分两种物质，T0代表热解燃烧室温度，u0代表热解燃烧室内所有物质的流速，ρN0代表热解燃烧室内固体颗粒数，dp,0代表热解燃烧室内煤焦颗粒直径，εp,0代表热解燃烧室内固相物质体积分数，εg,0代表热解燃烧室内气相物质体积分数，P代表气化炉操作压力，Tfeeding,x代表进入气化炉的x组分的温度，其中x包括煤粉、氧气、蒸汽、氮气、二氧化碳五种物质，mfeeding,x代表进入气化炉的x组分的质量流量，其中的x仍然包括煤粉、氧气、蒸汽、氮气、二氧化碳五种物质，dp,feeding代表进入气化炉的煤粉颗粒的直径，ρcoal代表进入气化炉的煤粉密度。
(4)计算本时刻各气化室内物质发生气化反应时的质量与能量源项，其中的气化反应速率(均相、非均相)由相关反应的动力学方程计算：
式中，代表第n个气化室内煤焦的非均相反应质量源项，代表第n个气化室内相关非均相反应m的反应速率，代表第n个气化室内灰分的非均相反应质量源项，xAsh代表灰分在煤粉中的质量分数，代表第n个气化室内气相组分i的非均相反应质量源项，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，xchar代表煤焦在煤粉中的质量分数，MC代表煤焦的摩尔质量，υi,m代表组分i在化学反应m中的化学计量系数，代表第n个气化室内气相组分i的均相反应质量源项，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，代表第n个气化室内相关均相反应m的反应速率，代表第n个气化室内所有固相组分的非均相反应能量源项，ΔHm代表化学反应m的反应焓，代表第n个气化室内所有气相组分的均相反应能量源项。
(5)计算本时刻各气化室与炉壁的传热源项，包括气相对流传热与气固辐射传热两部分。
Qconv,n＝Awα(Tn-Twall)
式中，Qconv,n代表第n个气化室内的对流传热量，Aw代表对流传热面积，α代表对流传热系数，Twall代表炉壁温度，Tn代表第n个气化室内的温度，Qrad,n代表第n个气化室内的辐射传热量，Ap代表辐射传热面积，σ代表玻尔兹曼常数，ep代表炉内固体颗粒黑度，ewall代表炉壁黑度。
(6)建立本时刻各气化室的质量、能量守恒的动态方程式，并计算流速与颗粒数：
气相质量守恒：
固相质量守恒：
颗粒数：
气固两相能量守恒：
流速：
式中，t代表时间，εg,n代表第n个气化室内气相物质体积分数，Cg,i,n代表第n个气化室内气相组分i的摩尔浓度，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，Ln代表第n个气化室的长度，un-1代表第n-1个气化室内的物质流速，εg,n-1代表第n-1个气化室内的气相物质体积分数，Cg,i,n-1代表第n-1个气化室内气相组分i的摩尔浓度，i仍然包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，un代表第n个气化室内的物质流速，εg,n代表第n个气化室内的气相物质体积分数，Cp,j,n代表第n个气化室内固相组分j的摩尔浓度，其中j包括煤焦与灰分两种物质，Cp,j,n-1代表第n-1个气化室内固相组分j的摩尔浓度，j仍然包括煤焦与灰分两种物质，ρNn代表第n个气化室内的固体颗粒数，ρNn-1代表第n-1个气化室内的固体颗粒数，Cpi,n代表第n个气化室内气相组分i的比热容，Cpj,n代表第n个气化室内固相组分j的比热容，Cpi,n-1代表第n-1个气化室内气相组分i的比热容，Cpj,n-1代表第n-1个气化室内固相组分j的比热容，Tn-1代表第n-1个气化室的温度，Tn代表第n个气化室的温度。其余质量、能量与传热源项 Qconv,n、Qrad,n的含义已在步骤四与步骤五中作过介绍，此处不再赘述。
(7)求解每个守恒方程关于时间的积分，输出本时刻气化炉各气化室温度的计算结果，再重复步骤(2)～(7)以进行下一时刻炉内温度的计算。
本发明基于实际应用的需要，根据气化炉内反应将气化炉分为一个热解燃烧室与多个气化室串联的模型结构，且热解燃烧室与各气化室皆视为全混流反应器。气化室的具体个数根据实际计算精度与计算速度的需要由用户自行设定，若要求计算精度较高则设定较多的气化室个数，若要求计算速度较快则设定较少的气化室个数。
经工业现场实验证明，本发明计算结果精度较高，计算速度也较快，具有很好的可行性。此外，本发明根据气化过程机理建立模型，能够得出炉内多个位置温度的实时数据，可以为工业现场操作人员控制气化过程提供参考，实现装置稳定生产，降低设备损耗。同时，还能够对深入研究气化过程机理与进一步控制气化过程提供帮助。
附图说明
图1为本发明的方法示意图。
图2为本发明对Shell气流床气化炉模型结构划分。
图3为本发明的操作流程示意图。
图4为本发明实施例在10h内的工业测量值与模拟值的对照图。
具体实施方式
本发明在实际应用时通过数据采集部分获得煤质数据、气化炉结构参数以及包括进煤量、进氧量在内的实时运行数据，据此依次计算该时间段内的炉内煤粉与气化剂发生各化学反应及传热的质量与能量守恒方程，最终得到该时间段内炉温的数据，再进行下一时间段内的计算。操作人员可根据每一时刻计算得到的温度数据判断气化炉内气化过程及装置的运行状态，这一过程如图1所示。
本发明将气化炉划分为如图2所示的一个热解燃烧室与多个气化室串联的模型结构，且热解燃烧室与各气化室皆视为全混流反应器。气化室的具体个数根据实际计算精度与计算速度的需要由用户自行设定，若要求计算精度较高则设定较多的气化室个数，若要求计算速度较快则设定较少的气化室个数。
图3展示了本发明的大致操作流程，以下结合一气化炉实例进一步说明本发明的具体实施方式。
实例选取河南能源化工集团下属某气化装置，其日处理煤量约2000t，反应器高度为6.8米，直径为2.96m，反应过程产生大量的热，热量通过水冷壁强制循环移除，反应器设计中压蒸汽产量为110t/h，装置进料煤粉采用CO2输送，设计正常操作压力为3.8～4.0MPa，气化温度约为1400～1600℃，碳转化率高达99％，产品气体中不含重烃，甲烷含量极低，合成气中有效气(CO+H2)含量可达85％以上。
数据采集部分包括入炉煤质数据、气化炉结构参数和包括进煤量、进氧量在内的气化炉10h内运行的实时数据，采集数据如表1～3所示。
表1煤质组成数据(％)
表2气化炉结构参数
表3 10h工业过程操作参数一览表
本方法具体步骤如下：
步骤一：根据计算精度与计算速度的要求确定气化室个数。
步骤二：接收数据采集部分传来的本时刻气化炉运行数据。
步骤三：根据进料参数计算本时刻进入气化炉热解燃烧室内的煤粉发生热解及燃烧反应的质量与能量平衡，其求解过程按下述表达式描述：
(Cg,i,0,Cp,j,0,T0,u0,ρN0,dp,0,εp,0,εg,0,P)
＝f(Tfeeding,x,mfeeding,x,dp,feeding,ρcoal)
式中Cg,i,0代表热解燃烧室内气相组分i的摩尔浓度，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，Cp,j,0代表热解燃烧室内固相组分j的质量浓度，其中j包括煤焦与灰分两种物质，T0代表热解燃烧室温度，u0代表热解燃烧室内所有物质的流速，ρN0代表热解燃烧室内固体颗粒数，dp,0代表热解燃烧室内煤焦颗粒直径，εp,0代表热解燃烧室内固相物质体积分数，εg,0代表热解燃烧室内气相物质体积分数，P代表气化炉操作压力，Tfeeding,x代表进入气化炉的x组分的温度，其中x包括煤粉、氧气、蒸汽、氮气、二氧化碳五种物质，mfeeding,x代表进入气化炉的x组分的质量流量，其中的x仍然包括煤粉、氧气、蒸汽、氮气、二氧化碳五种物质，dp,feeding代表进入气化炉的煤粉颗粒的直径，ρcoal代表进入气化炉的煤粉密度。
步骤四：计算本时刻各气化室内物质发生气化反应时的质量与能量源项，其中的气化反应速率(均相、非均相)由相关反应的动力学方程计算：
式中，代表第n个气化室内煤焦的非均相反应质量源项，代表第n个气化室内相关非均相反应m的反应速率，代表第n个气化室内灰分的非均相反应质量源项，xAsh代表灰分在煤粉中的质量分数，代表第n个气化室内气相组分i的非均相反应质量源项，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，xchar代表煤焦在煤粉中的质量分数，MC代表煤焦的摩尔质量，υi,m代表组分i在化学反应m中的化学计量系数，代表第n个气化室内气相组分i的均相反应质量源项，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，代表第n个气化室内相关均相反应m的反应速率，代表第n个气化室内所有固相组分的非均相反应能量源项，ΔHm代表化学反应m的反应焓，代表第n个气化室内所有气相组分的均相反应能量源项。
步骤五：计算本时刻各气化室与炉壁的传热源项，包括气相对流传热与气固辐射传热两部分。
Qconv,n＝Awα(Tn-Twall)
式中，Qconv,n代表第n个气化室内的对流传热量，Aw代表对流传热面积，α代表对流传热系数，Twall代表炉壁温度，Tn代表第n个气化室内的温度，Qrad,n代表第n个气化室内的辐射传热量，Ap代表辐射传热面积，σ代表玻尔兹曼常数，ep代表炉内固体颗粒黑度，ewall代表炉壁黑度。
步骤六：建立本时刻各气化室的质量、能量守恒的动态方程式，并计算流速与颗粒数：
气相质量守恒：
固相质量守恒：
颗粒数：
气固两相能量守恒：
流速：
式中，t代表时间，εg,n代表第n个气化室内气相物质体积分数，Cg,i,n代表第n个气化室内气相组分i的摩尔浓度，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，Ln代表第n个气化室的长度，un-1代表第n-1个气化室内的物质流速，εg,n-1代表第n-1个气化室内的气相物质体积分数，Cg,i,n-1代表第n-1个气化室内气相组分i的摩尔浓度，i仍然包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，un代表第n个气化室内的物质流速，εg,n代表第n个气化室内的气相物质体积分数，Cp,j,n代表第n个气化室内固相组分j的摩尔浓度，其中j包括煤焦与灰分两种物质，Cp,j,n-1代表第n-1个气化室内固相组分j的摩尔浓度，j仍然包括煤焦与灰分两种物质，ρNn代表第n个气化室内的固体颗粒数，ρNn-1代表第n-1个气化室内的固体颗粒数，Cpi,n代表第n个气化室内气相组分i的比热容，Cpj,n代表第n个气化室内固相组分j的比热容，Cpi,n-1代表第n-1个气化室内气相组分i的比热容，Cpj,n-1代表第n-1个气化室内固相组分j的比热容，Tn-1代表第n-1个气化室的温度，Tn代表第n个气化室的温度。其余质量、能量与传热源项 Qconv,n、Qrad,n的含义已在步骤四与步骤五中作过介绍，此处不再赘述。
步骤七：求解每个守恒方程关于时间的积分，输出本时刻气化炉各气化室温度的计算结果，再重复步骤二至步骤七以进行下一时刻炉内温度的计算。
上述步骤涉及微分方程的求解，计算过程通过XD-APC组态软件实现。图4为气化炉正常运行情况下，求解得到的第四个气化室温度与工业现场测量值的10h数据对比图。由于气化炉内温度无法直接测量，因此选择的工业现场测量值为进料氧煤比以及某一测温点温度的变化情况，该测温点位于合成气输气管上，合成气经激冷之后流入输气管并通过该点，由于激冷气温度变化范围不大，该点温度主要受气化炉内合成气温度的影响，因此从理论上讲，该点温度以及氧煤比的变化趋势应与气化室温度变化趋势一致。
从图中可以看到，气化室温度的模拟值(b)与氧煤比(c)、测温点温度(a)两个测量值的变化情况非常吻合，这一结果证明了本发明方法能够满足炉内温度实时软测量的要求。
本发明采用数学机理建模的方法建立气化炉的动态数学模型，根据数据采集系统接受煤质数据、装置参数及运行数据，通过在线求解数学模型的方法获得Shell气化炉炉内温度的软测量值。可以为操作人员提供及时与有效的炉内温度的信息，指导其更好得完成气化反应过程的优化与运行操作，克服由于Shell气化装置自身设计缺陷对生产过程及产品质量的影响。
Claims (6)
1.一种气流床气化炉炉膛温度的在线软测量方法，其特征在于包括以下步骤：
1)数据采集；所述数据包括入炉煤质数据、气化炉结构参数和气化炉运行实时数据；
2)计算输出；
在接受采集的数据后，通过相应的数学计算式及守恒方程，计算并输出炉内的温度数据给现场操作人员，为其判断气化炉内气化过程及装置的运行状态提供参考；
所述计算并输出炉内的温度数据包括炉内的化学反应和炉内向外传热的计算，以及炉温的输出整个软测量，具体方法如下：
(1)根据计算精度与计算速度的要求确定气化室个数；
(2)接收数据采集部分传来的本时刻气化炉运行数据；
(3)根据进料参数计算本时刻进入气化炉热解燃烧室内的煤粉发生热解及燃烧反应的质量与能量平衡，其求解过程按下述表达式描述：
(Cg,i,0,Cp,j,0,T0,u0,ρN0,dp,0,εp,0,εg,0,P)
＝f(Tfeeding,x,mfeeding,x,dp,feeding,ρcoal)
式中Cg,i,0代表热解燃烧室内气相组分i的摩尔浓度，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，Cp,j,0代表热解燃烧室内固相组分j的质量浓度，其中j包括煤焦与灰分两种物质，T0代表热解燃烧室温度，u0代表热解燃烧室内所有物质的流速，ρN0代表热解燃烧室内固体颗粒数，dp,0代表热解燃烧室内煤焦颗粒直径，εp,0代表热解燃烧室内固相物质体积分数，εg,0代表热解燃烧室内气相物质体积分数，P代表气化炉操作压力，Tfeeding,x代表进入气化炉的x组分的温度，其中x包括煤粉、氧气、蒸汽、氮气和二氧化碳五种物质，mfeeding,x代表进入气化炉的x组分的质量流量，其中的x仍然包括煤粉、氧气、蒸汽、氮气和二氧化碳五种物质，dp,feeding代表进入气化炉的煤粉颗粒的直径，ρcoal代表进入气化炉的煤粉密度；
(4)计算本时刻各气化室内物质发生气化反应时的质量与能量源项，其中的气化反应速率由相关反应的动力学方程计算：
式中，代表第n个气化室内煤焦的非均相反应质量源项，代表第n个气化室内相关非均相反应m的反应速率，代表第n个气化室内灰分的非均相反应质量源项，xAsh代表灰分在煤粉中的质量分数，代表第n个气化室内气相组分i的非均相反应质量源项，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，xchar代表煤焦在煤粉中的质量分数，MC代表煤焦的摩尔质量，υi,m代表组分i在化学反应m中的化学计量系数，代表第n个气化室内气相组分i的均相反应质量源项，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，代表第n个气化室内相关均相反应m的反应速率，代表第n个气化室内所有固相组分的非均相反应能量源项，ΔHm代表化学反应m的反应焓，代表第n个气化室内所有气相组分的均相反应能量源项；
(5)计算本时刻各气化室与炉壁的传热源项，包括气相对流传热与气固辐射传热两部分；
Qconv,n＝Awα(Tn-Twall)
式中，Qconv,n代表第n个气化室内的对流传热量，Aw代表对流传热面积，α代表对流传热系数，Twall代表炉壁温度，Tn代表第n个气化室内的温度，Qrad,n代表第n个气化室内的辐射传热量，Ap代表辐射传热面积，σ代表玻尔兹曼常数，ep代表炉内固体颗粒黑度，ewall代表炉壁黑度；
(6)建立本时刻各气化室的质量和能量守恒的动态方程式，并计算流速与颗粒数：
气相质量守恒：
固相质量守恒：
颗粒数：
气固两相能量守恒：
流速：
式中，t代表时间，εg,n代表第n个气化室内气相物质体积分数，Cg,i,n代表第n个气化室内气相组分i的摩尔浓度，其中i包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，Ln代表第n个气化室的长度，un-1代表第n-1个气化室内的物质流速，εg,n-1代表第n-1个气化室内的气相物质体积分数，Cg,i,n-1代表第n-1个气化室内气相组分i的摩尔浓度，i仍然包括二氧化碳、水、氢气、一氧化碳、甲烷、氮气以及硫化氢七种物质，un代表第n个气化室内的物质流速，εg,n代表第n个气化室内的气相物质体积分数，Cp,j,n代表第n个气化室内固相组分j的摩尔浓度，其中j包括煤焦与灰分两种物质，Cp,j,n-1代表第n-1个气化室内固相组分j的摩尔浓度，j仍然包括煤焦与灰分两种物质，ρNn代表第n个气化室内的固体颗粒数，ρNn-1代表第n-1个气化室内的固体颗粒数，Cpi,n代表第n个气化室内气相组分i的比热容，Cpj,n代表第n个气化室内固相组分j的比热容，Cpi,n-1代表第n-1个气化室内气相组分i的比热容，Cpj,n-1代表第n-1个气化室内固相组分j的比热容，Tn-1代表第n-1个气化室的温度，Tn代表第n个气化室的温度；
(7)求解每个守恒方程关于时间的积分，输出本时刻气化炉各气化室温度的计算结果，再重复步骤(2)～(7)以进行下一时刻炉内温度的计算。
2.如权利要求1所述一种气流床气化炉炉膛温度的在线软测量方法，其特征在于所述入炉煤质数据由相关人员通过煤质分析获得。
3.如权利要求1所述一种气流床气化炉炉膛温度的在线软测量方法，其特征在于所述气化炉结构参数包括气化炉高度、气化炉直径和气化炉水冷壁厚度。
4.如权利要求1所述一种气流床气化炉炉膛温度的在线软测量方法，其特征在于所述气化炉结构参数通过气化炉设计图纸获得。
5.如权利要求1所述一种气流床气化炉炉膛温度的在线软测量方法，其特征在于所述气化炉运行实时数据通过工厂DCS系统采集。
6.如权利要求1所述一种气流床气化炉炉膛温度的在线软测量方法，其特征在于所述气化炉运行实时数据包括进煤量、进氧量、蒸汽量、载气量以及气化炉压力，所述载气为二氧化碳或氮气。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610267960.3A CN105930666B (zh) | 2016-04-27 | 2016-04-27 | 一种气流床气化炉炉膛温度的在线软测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610267960.3A CN105930666B (zh) | 2016-04-27 | 2016-04-27 | 一种气流床气化炉炉膛温度的在线软测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105930666A CN105930666A (zh) | 2016-09-07 |
| CN105930666B true CN105930666B (zh) | 2019-01-22 |
Family
ID=56836424
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610267960.3A Expired - Fee Related CN105930666B (zh) | 2016-04-27 | 2016-04-27 | 一种气流床气化炉炉膛温度的在线软测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105930666B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108956693B (zh) * | 2018-04-28 | 2020-11-27 | 中国科学院广州能源研究所 | 一种氢气火焰温度测量方法 |
| CN109726493B (zh) * | 2019-01-04 | 2021-03-23 | 厦门大学 | 一种渣油气化炉炉膛内气体组分浓度在线软测量方法 |
| CN109920486B (zh) * | 2019-01-10 | 2021-03-16 | 江苏理工学院 | 一种基于Shell语言提高分子动力学批量建模效率的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6865440B1 (en) * | 2004-02-26 | 2005-03-08 | Sodick Co., Ltd. | Numeric control apparatus for machine tools |
| CN101038277B (zh) * | 2007-04-19 | 2010-12-08 | 东北大学 | 基于最小二乘-支持向量机的制粉过程煤粉细度软测量方法 |
| CN103150480A (zh) * | 2013-03-19 | 2013-06-12 | 贵州师范大学 | 一种宏观气固反应动力学模型 |
| CN104087344A (zh) * | 2014-05-05 | 2014-10-08 | 董金国 | 一种气流床气化炉炉膛温度测量方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7079908B2 (en) * | 2003-01-31 | 2006-07-18 | Delphi Technologies,Inc. | Horizontally-structured modeling for analysis |
-
2016
- 2016-04-27 CN CN201610267960.3A patent/CN105930666B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6865440B1 (en) * | 2004-02-26 | 2005-03-08 | Sodick Co., Ltd. | Numeric control apparatus for machine tools |
| CN101038277B (zh) * | 2007-04-19 | 2010-12-08 | 东北大学 | 基于最小二乘-支持向量机的制粉过程煤粉细度软测量方法 |
| CN103150480A (zh) * | 2013-03-19 | 2013-06-12 | 贵州师范大学 | 一种宏观气固反应动力学模型 |
| CN104087344A (zh) * | 2014-05-05 | 2014-10-08 | 董金国 | 一种气流床气化炉炉膛温度测量方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| "水煤浆气化炉炉温智能软测量建模应用研究";李杰;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20120615(第6期);第B017-30页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105930666A (zh) | 2016-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Petersen et al. | Experimental investigation and modeling of gasification of sewage sludge in the circulating fluidized bed | |
| CN105930666B (zh) | 一种气流床气化炉炉膛温度的在线软测量方法 | |
| Kolb et al. | Tackling the challenges in modelling entrained-flow gasification of low-grade feedstock | |
| Tagutchou et al. | A new experimental Continuous Fixed Bed Reactor to characterise wood char gasification | |
| CN102799748B (zh) | 煤气化炉的控制方法 | |
| CN1011076B (zh) | 量热系统 | |
| Herdel et al. | Experimental investigations in a demonstration plant for fluidized bed gasification of multiple feedstock’s in 0.5 MWth scale | |
| Golub et al. | Experimental study into the influence of straw content in fuel on parameters of generator gas | |
| Walawender et al. | Steam gasification of pure cellulose. 1. Uniform temperature profile | |
| Moldenhauer et al. | Experimental investigation of chemical-looping combustion and chemical-looping gasification of biomass-based fuels using steel converter slag as oxygen carrier | |
| CN109726493B (zh) | 一种渣油气化炉炉膛内气体组分浓度在线软测量方法 | |
| Punnarapong et al. | Development and performance evaluation of a biomass gasification system for ceramic firing process | |
| Kolbitsch | First fuel tests at a novel 100 kWth dual fluidized bed steam gasification pilot plant | |
| Foong et al. | Coal gasification in a spouted bed | |
| CA2975325C (en) | Standpipe-fluid bed hybrid system for char collection, transport, and flow control | |
| Chowdhury et al. | Modelling and simulation of a downdraft rice husk gasifier | |
| Perry et al. | Continuous gasification of pulverized coal with oxygen and steam by the vortex principle | |
| Lewis et al. | Strongly Caking Coal Gasified in a Stirred-Bed Producer | |
| CN103601375B (zh) | 并流蓄热式石灰窑燃烧膛内气体、固体温度的确定方法 | |
| Foster | Production of Water Gas from Pulverixed | |
| Engelbrecht | Characterisation and fluidised bed gasification of selected high-ash South African coals | |
| US1549160A (en) | Gas apparatus | |
| CN113984635A (zh) | 一种耐火砖层的蚀损速率监测装置和监测方法 | |
| CN109453608A (zh) | 一种回收锅炉余热的装置及制备方法 | |
| Kumararaja | OPERATION OF PACKED BED DOWNDRAFT BIOMASS GASIFIER WITH DIFFERENT AIR SUPPLY ARRANGEMENTS |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190122 |