CN112487641A - 基于多尺度建模的煤化学链气化工艺设计方法及其应用 - Google Patents
基于多尺度建模的煤化学链气化工艺设计方法及其应用 Download PDFInfo
- Publication number
- CN112487641A CN112487641A CN202011363689.6A CN202011363689A CN112487641A CN 112487641 A CN112487641 A CN 112487641A CN 202011363689 A CN202011363689 A CN 202011363689A CN 112487641 A CN112487641 A CN 112487641A
- Authority
- CN
- China
- Prior art keywords
- coal
- clg
- simulation
- gasification process
- modeling
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 134
- 230000008569 process Effects 0.000 title claims abstract description 90
- 239000003245 coal Substances 0.000 title claims abstract description 79
- 239000000126 substance Substances 0.000 title claims abstract description 34
- 238000002309 gasification Methods 0.000 title claims abstract description 32
- 238000013461 design Methods 0.000 title claims abstract description 25
- 238000004088 simulation Methods 0.000 claims abstract description 66
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 33
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 30
- 239000001301 oxygen Substances 0.000 claims abstract description 30
- 239000007789 gas Substances 0.000 claims abstract description 27
- 238000000329 molecular dynamics simulation Methods 0.000 claims abstract description 25
- 238000004519 manufacturing process Methods 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 20
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 20
- 230000007246 mechanism Effects 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 239000013256 coordination polymer Substances 0.000 claims abstract description 8
- 239000000446 fuel Substances 0.000 claims abstract description 8
- 239000012530 fluid Substances 0.000 claims abstract description 7
- 239000012634 fragment Substances 0.000 claims abstract description 6
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims description 23
- 241000183024 Populus tremula Species 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 229910052799 carbon Inorganic materials 0.000 claims description 5
- 238000007792 addition Methods 0.000 claims description 4
- 238000000137 annealing Methods 0.000 claims description 4
- 229910052739 hydrogen Inorganic materials 0.000 claims description 4
- 238000002156 mixing Methods 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 229910052717 sulfur Inorganic materials 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 2
- 239000002994 raw material Substances 0.000 claims description 2
- 238000005265 energy consumption Methods 0.000 abstract description 2
- 238000000197 pyrolysis Methods 0.000 description 18
- QPLDLSVMHZLSFG-UHFFFAOYSA-N Copper oxide Chemical compound [Cu]=O QPLDLSVMHZLSFG-UHFFFAOYSA-N 0.000 description 12
- 238000012546 transfer Methods 0.000 description 12
- 239000000571 coke Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 239000007787 solid Substances 0.000 description 6
- 238000005094 computer simulation Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000000969 carrier Substances 0.000 description 3
- 239000010949 copper Substances 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 230000005610 quantum mechanics Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 239000005751 Copper oxide Substances 0.000 description 2
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N Iron oxide Chemical compound [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 description 2
- 125000004429 atom Chemical group 0.000 description 2
- 125000004432 carbon atom Chemical group C* 0.000 description 2
- 229910000431 copper oxide Inorganic materials 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000005411 Van der Waals force Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002802 bituminous coal Substances 0.000 description 1
- 238000006664 bond formation reaction Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000006229 carbon black Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010351 charge transfer process Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N iron(III) oxide Inorganic materials O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000324 molecular mechanic Methods 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000006068 polycondensation reaction Methods 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000013341 scale-up Methods 0.000 description 1
- 239000011269 tar Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Industrial Gases (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明涉及清洁生产领域领域，具体涉及基于多尺度建模的煤化学链气化工艺设计方法及其应用。建立煤模型并进行全局能量优化，采用ReaxFF‑MD模拟对CP进行机理分析，通过分析碎片数量得到最优CP温度；采用CFD对CLG过程进行模拟，对燃料反应器进行建模，确定载氧体的种类和最佳停留时间；基于CP机理分析和CLG流体动力学模拟结果，对CP‑CLG系统进行稳态模拟，并通过灵敏度分析得到影响合成气生产的关键参数。通过多尺度建模获得高纯度的合成气，考虑清洁生产各项性能的综合优化，能同时实现清洁生产、低能耗、高生产率以及多反应器动态稳定性，保证了工艺安全。
Description
技术领域
本发明涉及清洁生产领域，具体涉及基于多尺度建模的煤化学链气化工艺设计方法及其应用。
背景技术
公开该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不必然被视为承认或以任何形式暗示该信息构成已经成为本领域一般技术人员所公知的现有技术。
煤炭是我国能源结构中占主导地位的化石燃料之一。但传统的煤炭利用技术能源效率低，污染问题严重。因此，洁净煤技术已成为近年来研究和开发的重点。煤热解过程(CP)是一种更清洁的煤转化方式，可生产高附加值的热解气、焦和焦油。化学链气化(CLG)是另一种利用载氧体实现高效能源转化和内部CO2分离的清洁技术。因此，CP和CLG过程的集成可以比单个过程同时提高生产效率和减少环境污染。
CP是一个复杂的过程，涉及无数的反应途径，目前尚缺乏对其机理的详细描述。分子动力学(MD)模拟是一种适用于大分子系统的有效建模方法。反力场(ReaxFF)是为了弥补MD与量子力学(QM)之间的差距而发展起来的。它适用于原子数不超过一百万个的各种化学系统，不受量子力学中原子数的限制。因此，ReaxFF-MD为CP等复杂反应体系的MD模拟提供了一种有前景的方法。
但是CP和CLG过程复杂，在模拟过程中，处理对象的复杂程度不同直接影响建模过程中的参数选择、难度和模拟结果。发明人发现，现有技术虽然有公开CP和CLG的建模方法，但单一尺度(仅模拟CP过程或仅模拟CLG过程)的建模分析无法研究CP和CLG过程参数在CP-CLG完整过程中的相互配合作用。而使用一种建模构思模拟CP-CLG完整过程的参数选择，因为没有针对CP和CLG反应过程的不同特点和重点参数进行模拟，无法兼顾CP和CLG过程的不同特点，因此仅使用一种建模构思构建的CP-CLG过程参数的选择有待完善，反应产物的稳定性有待提高。
还有一些建模方法处理对象理想化，与实际生产过程中的煤构型和反应过程有较大偏差，有的侧重研究炭黑产率，无法兼顾合成气质量，有的合成气质量流量虽然大，但是处理对象理想化导致缺乏实际参考，无法进行动态研究，因此不能满足工业合成气的生产要求。
发明内容
针对上述现有技术中存在的问题，本发明采用基于多尺度建模的煤化学链气化工艺设计方法。以CP-CLG流程模拟为核心，利用MD模拟对CP过程进行机理研究，再利用CFD(计算流体力学)模拟对CLG过程进行建模和分析，为CP-CLG稳态模拟提供了最优的操作参数，互相配合，实现CP-CLG稳态运行。并在CP-CLG稳态模拟基础上，建立了多反应器系统的动态控制方案，并通过添加不同扰动测试了高纯度合成气的稳定生产。
具体地，本发明是通过如下所述的技术方案实现：
本发明第一方面，本发明提供一种基于多尺度建模的煤化学链气化工艺设计方法，包括：
建立煤模型并进行全局能量优化，采用ReaxFF-MD模拟对CP进行机理分析，通过分析碎片数量得到最优CP温度；
采用CFD对CLG过程进行模拟，对燃料反应器进行建模，确定载氧体的种类和最佳停留时间；
基于CP机理分析和CLG流体动力学模拟结果，对CP-CLG系统进行稳态模拟，并通过灵敏度分析得到影响合成气生产的关键参数。
本发明第二方面，提供一种基于多尺度建模的煤化学链气化工艺设计方法在煤化学链气化工艺中的应用。
本发明第三方面，提供一种煤化学链气化生产方法，采用基于多尺度建模的煤化学链气化工艺设计方法进行模拟。
本发明的一个或多个实施例具有以下有益效果：
(1)选用含有更多元素，碳原子数更多的煤构型进行研究，真实地模拟实际生产过程中煤的CP-CLG过程，该建模方法能对碳原子数少的煤构型同样有效。
(2)CP机理研究主要包括煤模型的全局优化和热解过程的分析，过程复杂，涉及很多个复杂反应，如果通过CFD进行模拟的话很难建模，缺少大量的反应动力学数据，并且很难得到准确的结果。为了提高模拟准确度，选择MD模拟CP。模拟结果表明，热解温度需在1100K以下进行，与实际运行参数相同，说明CP模拟过程准确度高。
(3)CLG过程的CFD模拟主要是通过反应器设计建模，以载氧体为分散相，以水蒸气为连续相的两者混合程度来确定的停留时间。选择了CuO作为载氧体，确定其停留时间为2s。实验发现，Aspen Plus软件模拟CP和CLG过程，无法兼顾CP和CLG过程不同的特点，模拟获得的参数契合度低，直接表现为合成气质量低，产率低。与Aspen Plus软件模拟CLG过程相比，MD模拟CP过程获得的参数能与CFD模拟CLG过程获得的参数更好的契合，形成稳定的模拟体系，合成气质量更高。
(4)利用Aspen Plus软件对CP-CLG系统进行稳态模拟，在该建模方法处理更复杂的煤构型获得的工艺参数运行条件下，合成气质量流量达到508.5kg/h，能满足工业合成气的生产要求。
(5)在该建模方法处理更复杂的煤构型获得的工艺参数运行条件下，通过动态仿真对CP-CLG系统的控制结构进行了验证，证明了在10％的扰动下合成气生产是稳定进行的。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。以下，结合附图来详细说明本发明的实施方案，其中：
图1为本发明实施例1中最低能量煤模型；
图2为本发明实施例1中MD模拟结果(C1-C3，CO2)；
图3为本发明实施例1中MD模拟结果(H2，CO)；
图4为本发明实施例1中CLG过程装置尺寸图；
图5为本发明实施例1中CuO载氧体传热性能图；
图6为本发明实施例1中Fe2O3载氧体传热性能图；
图7为本发明实施例1中CP-CLG稳态流程图；
图8为本发明实施例1中CP-CLG动态控制方案图；
图9为本发明实施例1中添加10％蒸汽流量扰动CP-CLG过程H2动态曲线图
图10为本发明实施例1中添加10％蒸汽流量扰动CP-CLG过程CO动态曲线图
图11为本发明实施例1中添加10％温度扰动CP-CLG过程H2动态扰动曲线图
图12为本发明实施例1中添加10％温度扰动CP-CLG过程CO动态扰动曲线图。
具体实施方式
下面结合具体实施例，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。下列实施例中未注明具体条件的实验方法，通常按照常规条件或按照制造厂商所建议的条件。
除非另行定义，文中所使用的所有专业与科学用语与本领域熟练人员所熟悉的意义相同。此外，任何与所记载内容相似或均等的方法及材料皆可应用于本发明方法中。文中所述的较佳实施方法与材料仅作示范之用。
针对现有技术研究CP-CLG系统存在的处理对象简单和单一建模效果差的问题，本发明提出一种基于多尺度建模与控制的煤化学链气化工艺设计方法。
综合MD和CFD的多尺度建模适用于定量、综合地描述、预测和理解复杂机制。AspenDynamics是一种流行的动态仿真工具，它将多个学科合并到动态仿真中。因此，本发明利用它通过对控制系统进行动态响应测试。本发明采用多尺度建模方法建立了CP-CLG仿真系统，然后针对关键工艺参数对合成气纯度的影响，设计了多反应器联合动态控制系统。
通过多尺度建模获得高纯度的合成气，考虑清洁生产各项性能的综合优化，能同时实现清洁生产、低能耗、高生产率以及多反应器动态稳定性，保证了工艺安全。
为了实现上述技术效果，本发明的第一方面提供如下技术方案：
一种基于多尺度建模与控制的煤化学链气化工艺设计方法，包括：
建立煤模型并进行全局能量优化，采用ReaxFF-MD模拟对CP进行机理分析，通过分析碎片数量得到最优CP温度；
采用CFD对CLG过程进行模拟，对燃料反应器进行建模，确定载氧体的种类和最佳停留时间；
基于CP机理分析和CLG流体动力学模拟结果，对CP-CLG系统进行稳态模拟，并通过灵敏度分析得到影响合成气生产的关键参数。
CP模拟的难点在于模型结构复杂，运行速度慢，不进行退火模拟很难收敛，因此重点和难点在于结构优化。CLG模拟的难点在于反应动力学参数的确立以及编写的准确性，一旦编写出现错误很难得到准确的结果。
而且CP过程复杂，涉及很多个复杂反应，如果通过CFD进行模拟的话很难建模，缺少大量的反应动力学数据，并且很难得到准确的结果，因此选择MD模拟。载氧体的选择和停留时间的确定只能通过CFD模拟确定，所以目标不同，选择的方法不同。而如何使两个过程模拟的参数很好契合，共同形成完整的CP-CLG过程，在稳态模拟过程中实现较好的效果，成为技术难题之一。
为了获得不同的参数，所以选择了不同的模拟方法。结合MD和CFD模拟数据，更好地建立CP与CLG过程的联系。使CP与CLG过程相互促进，煤模型在CP确定的最佳热解温度下热解，得到产物焦炭和焦油，以焦炭为反应底物，结合反应器结构特点，以及传热性能，确定载氧剂种类和载氧剂停留时间，这些参数的确定为接下来的CP-CLG过程模拟提供了坚实的基础。
由于MD模拟的只是煤分子的结构单元，不符合实际企业，得到的仅仅是热解温度，不足够支撑Aspen Plus稳态模拟。所以本发明完善了煤的建模，匹配了实际煤的结构，并且得到的煤热解温度与实际企业运行温度一致，达到了很好的精度。除此之外，本发明用CFD模拟对载氧体进行选择，并确定其停留时间，使Aspen Plus稳态模拟更加全面，结果更令人信服。
为了增加建模的真实性，在本发明的一个或多个实施例中，所述煤模型为含C、H、O、N、S元素的模型。
为了真实反映复杂煤在化学链气化过程中的反应过程，在本发明的一个或多个实施例中，所述煤模型选自碳原子数为500-700的煤分子构型。
在本发明的一个或多个实施例中，通过退火循环模拟，所述煤模型选自3种分子式相加为C1992H1755N33O74S6的煤分子构型，组装成密度为1.3g/m3的立方三维分子模型。
在本发明的一个或多个实施例中，CFD模拟载氧体停留时间主要是通过反应器设计建模，以载氧体为分散相，以水蒸气为连续相的两者混合程度来确定的停留时间。
本发明的核心构思是进行CP-CLG的Aspen Plus稳态模拟，CFD是对CLG的模拟，目的是指导CP-CLG的Aspen Plus稳态模拟。CFD对CLG的模拟能指导Aspen Plus模拟，主要体现在载氧体的选择和停留时间的确定，旨在为Aspen Plus模拟提供更多的参数支持，避免了Aspen Plus模拟过程参数不完整的情况出现，对Aspen Plus模拟的结果准确性提供了保障。
在本发明的一个或多个实施例中，在通过Aspen Dynamics软件验证CP-CLG系统的抗干扰能力步骤中，直接采用CP产物作为原料。
在本发明的一个或多个实施例中，燃料反应器中CuO与水蒸气反应，压力和温度分别设定为3.0-4.0MPa和1220-1300K，优选为3.5MPa和1223K；
或者，在空气反应器Cu与氧的反应中，压力和温度分别设置为1.5-2.5MPa和1350-1450K，优选为2MPa和1400K。
在本发明的一个或多个实施例中，基于关键参数的选择，通过Aspen Dynamics软件验证CP-CLG系统的抗干扰能力，提出CP-CLG系统的控制方案；
优选地，验证抗干扰能力过程中，干扰条件为在2h时加入10％的蒸汽流量，观察合成气中H2和CO纯度的变化。
经过本发明的模拟可以通过合成气纯度这个指标来判断，本发明的效果比任何单独的过程都要好。
在本发明的第二方面，本发明提供一种基于多尺度建模与控制的煤化学链气化工艺设计方法在煤化学链气化工艺中的应用。
在本发明的第三方面，本发明提供一种基于多尺度建模与控制的煤化学链气化工艺设计方法进行模拟。
实施例1
一种基于多尺度建模与控制的煤化学链气化清洁生产工艺设计，对CP-CLG系统的建立，其步骤具体如下：
(1)建立煤模型并进行了全局的能量优化，并采用ReaxFF-MD模拟对CP进行机理分析，通过分析碎片数量得到了最优的CP温度。
首先建立煤的结构单元来表示分子骨架。然后用可调氢对煤分子结构单元进行加氢，得到合理的初始构型。最后，以能量为目标函数，对初始构型进行分子力学模拟，优化其几何结构。具体参数如下：任务类型为几何优化，优化质量为中等，能源优化方法为智能。用电荷平衡法计算了库仑力和范德瓦尔斯力。能量收敛质量中等，能量收敛最大阈值为0.001kcal·mol-1，力收敛最大阈值为
最大迭代次数为500次。通过退火循环模拟，可以得到能量较低的5种煤分子构型(C1,C2,C3,C4,C5)，其总势能如表1所示。为保证烟煤模型中C、H、O、N、S的合理组成，将分子式为C1992H1755N33O74S6的C1结构3块组装成密度为1.3g/m3的立方三维分子模型。在单元中加入周期性边界条件，以保持粒子数目不变。采用镜像法计算了边界原子间的力，使边界力均匀，消除了边界效应。在优化了单元的几何形状后，通过ReaxFF MD模拟得到煤模型的最低能量几何形状，如图1所示。在此模型基础上，ReaxFF MD方法模拟了CP过程，因为ReaxFF MD可以根据键序概念描述复杂反应分子体系的键形成和电荷转移过程。ReaxFF MD模拟中收集煤模型的细胞体积、原子类型、原子序数和原子坐标信息进行热解。热解初始温度为450K，之后在NVT条件(等粒子数、等体积、等温度、的环境)下以20K/ps的速率升高至2450K。主要运行条件为:步长为0.25fs，总模拟时间为100ps，设置键截断系数为0.3。温度从450K升高到2450K时碎片数量变化如图2和图3所示。可以看出，在500k后，由于C-O键等弱键的断裂，煤分子开始分解为C1-C3、CO2等各种元素分子气体。大约1100k后，由于半焦缩聚反应，碎片中出现大量的H2分子。同时，C-C键开始断裂，随着温度的升高，C1-C3碎片、CO2碎片和CO碎片也增加。因此，为了在后续CLG过程中获得更清洁的合成气，应将热解温度设置在1100k以下，以抑制CP中过量出现H2。
表1. 5种煤分子构型的能量组成
(2)通过对比不同氧载体的传热性能，选择CLG过程中最佳的氧载体。然后采用CFD模拟对燃料反应器进行建模，确定载氧体的最佳停留时间。
采用CFD方法对其进行模拟。与传统的气化技术不同，煤CLG系统通常由FR和AR组成，氧气载体在两者之间循环。如图4所示，该过程随FR的几何形状和设备尺寸而变化。FR几何形状越复杂，设备尺寸越大，CFD的收敛性就越差。从FR底部开始的煤与下落的高温载氧颗粒充分接触。FR内部循环流场的存在使煤在FR内停留较长时间，有利于气化反应的全面进行。CFD是分析流体动力学和化学动力学之间相互作用的有用方法，可用于优化操作和扩大单元操作规模。与MD适用于复杂反应机理研究不同，CFD可以耦合气固分离领域中的质量、热量和动量传递现象。本实施例CLG的CFD模拟，包括载氧体的选择、FR的CFD建模和两相(载氧体和水蒸气)混合。
如式(1)和(2)所示，固体的传热速度主要取决于密度(冲压)、定压热容量(Cp)和热导率(k)等物理性质。将固体传热物理场与传质耦合，研究了氧化铜和氧化铁的传热性能。CFD模拟得到的CuO和Fe2O3的温度分布如图5和图6所示。可以看出CuO的温度在2s内从1200K下降到1000K，而Fe2O3的温度下降同样需要3s。使用时间越短，传热率越高。因此，从传热角度出发，选择CuO作为氧载体。
(3)基于CP机理分析和CLG流体动力学模拟结果，对CP-CLG系统进行了稳态模拟，并通过灵敏度分析得到影响合成气生产的关键参数。
通过以上讨论，结合MD和CFD模拟数据，更好地建立CP与CLG过程的联系。使CP与CLG过程相互促进，煤模型在CP确定的最佳热解温度下热解，得到产物焦炭和焦油，以焦炭为反应底物，结合反应器结构特点，以及传热性能，确定载氧剂种类和载氧剂停留时间，这些参数的确定为接下来的CP-CLG过程模拟提供了坚实的基础。
整个CP-CLG体系主要由热解、FR、AR组成，进行CP-CLG模拟的前提是定义组分。在Aspen Plus数据库中，煤、灰分和煤焦被视为非常规组分。性质方法选用Peng Robinson-Boston Mathias(PR-BM)方程，因为它可以得到煤等复杂组分的具体性质。另外，ReaxFF MD模拟的热解温度应设置在1100K以下，因此热解压力和温度分别为1MPa和1000K。在CLG中，选用CuO作为氧载体，通过CFD模拟确定CuO的停留时间为2s。其它工艺操作参数如温度、压力、设备尺寸等均由某煤化工企业实际确定。所建立的CP-CLG系统流程图如图7所示，进行CP与CLG的耦合过程。在稳态CP模拟中，采用RYIELD单元模拟热解和干燥反应器。干燥反应器的压力和温度分别设置为0.1MPa和400k。1号分离装置是通过除去大部分水来获得干煤。经过热解反应器，焦炭、热煤气和焦油通过2号分离器单元进行分离。然后焦炭进入焦炭单元生产焦炭。在分解过程中，分解的焦炭进入FR(燃料反应器)与水蒸汽接触，生成合成气。此过程中所有的反应均列于式(3)和(4)。AR和FR反应器采用RGIBBS表示，允许计算含固体成分的反应。在FR中CuO与水蒸汽反应如式(5)所示，经过压力和温度分别设定为3.5MPa和1223K的FR反应器后，合成气进一步清洗，通过4号分离器装置得到H2和CO。AR(空气反应器)反应器中Cu与氧的反应如式(6)所示，其压力和温度分别设置为2MPa和1400K。最后，氧化铜通过6号分离器从AR循环到FR。
COAL→C+O2+H2+H2O+CH4+N2+S+ASH (3)
C+H2O→CO+H2 (4)
C+CuO→CO+Cu (5)
2Cu+O2→2CuO (6)
(4)基于关键参数的选择，通过Aspen Dynamics软件验证CP-CLG系统的抗干扰能力，提出了CP-CLG系统的控制方案。
在上述过程仿真的基础上，还需要进行动态仿真来验证产品在面对外界干扰时的稳定性，探讨CP-CLG系统中有关合成气质量的稳定性控制。因此，对包含多种物质的复杂动力学模型进行简化如下：该动力学模型直接采用CP产物作为原料，原因是其热解反应复杂且动力学参数不严格。表2列出了动态过程所需的所有新增阀门的压降和管式反应器的尺寸。表3列出了所有的控制器参数，包括9个主要参数的控制器，如操纵变量、控制变量、积分时间、增益。
表2.CP-CLG动态过程所需的所有新增阀门的压降和管式反应器的尺寸
| Tube Reactor | R01 | R02 |
| Length,m | 1 | 0.8 |
| Diameter,m | 0.5 | 0.4 |
| Valve | V1,V3,V4,V5 | V2,V6,V7,V8,VI,VII |
| Pressure drop,MPa | 0.1 | 0.2 |
| Valid phase | Vapor-Liquid | Vapor only |
表3.CP-CLG动态过程中控制器参数
基于以上设计，对CP-CLG系统进行动态仿真，如图8所示。图9-12显示了在2h时加入10％的蒸汽流量或改变温度后，合成气中H2和CO纯度的变化。图9和图10显示了在2h时加入10％的蒸汽流量后，合成气中H2和CO纯度的变化。结果表明，合成气中H2的纯度逐渐提高到0.281，合成气中CO的纯度逐渐接近0.38。在R01温度升高10％的2小时后，合成气中H2和CO的纯度变化如图11和图12所示，可以看出，合成气中H2的纯度逐渐下降到0.245，CO的纯度逐渐接近0.41。图9-12说明了图8所示控制方案的有效性。本申请的控制结构更为复杂，使用不同建模构思分别模拟CP和CLG过程，结果可靠性更高。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，包括：
建立煤模型并进行全局能量优化，采用ReaxFF-MD模拟对CP进行机理分析，通过分析碎片数量得到最优CP温度；
采用CFD对CLG过程进行模拟，对燃料反应器进行建模，确定载氧体的种类和最佳停留时间；
基于CP机理分析和CLG流体动力学模拟结果，对CP-CLG系统进行稳态模拟，并通过灵敏度分析得到影响合成气生产的关键参数。
2.根据权利要求1所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，所述煤模型为含C、H、O、N、S元素的模型。
3.根据权利要求1所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，所述煤模型选自碳原子数为500-700的煤分子构型。
4.根据权利要求3所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，通过退火循环模拟，所述煤模型选自3种分子式相加为C1992H1755N33O74S6的煤分子构型，组装成密度为1.3g/m3的立方三维分子模型。
5.根据权利要求1所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，CFD模拟载氧体停留时间主要是通过反应器设计建模，以载氧体为分散相，以水蒸气为连续相的两者混合程度来确定的停留时间。
6.根据权利要求1所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，在通过Aspen Dynamics软件验证CP-CLG系统的抗干扰能力步骤中，直接采用CP产物作为原料。
7.根据权利要求1所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，燃料反应器中CuO与水蒸气反应，压力和温度分别设定为3.0-4.0MPa和1220-1300K，优选为3.5MPa和1223K；
或者，在空气反应器Cu与氧的反应中，压力和温度分别设置为1.5-2.5MPa和1350-1450K，优选为2MPa和1400K。
8.根据权利要求1所述的基于多尺度建模的煤化学链气化工艺设计方法，其特征在于，基于关键参数的选择，通过Aspen Dynamics软件验证CP-CLG系统的抗干扰能力，提出CP-CLG系统的控制方案；
优选地，验证抗干扰能力过程中，干扰条件为在2h时加入10％的蒸汽流量，观察合成气中H2和CO纯度的变化。
9.权利要求1至8中任一项所述的基于多尺度建模的煤化学链气化工艺设计方法在煤化学链气化工艺中的应用。
10.一种煤化学链气化生产方法，其特征在于，采用权利要求1至8中任一项所述的基于多尺度建模的煤化学链气化工艺设计方法进行模拟。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011363689.6A CN112487641B (zh) | 2020-11-27 | 2020-11-27 | 基于多尺度建模的煤化学链气化工艺设计方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011363689.6A CN112487641B (zh) | 2020-11-27 | 2020-11-27 | 基于多尺度建模的煤化学链气化工艺设计方法及其应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112487641A true CN112487641A (zh) | 2021-03-12 |
| CN112487641B CN112487641B (zh) | 2023-04-07 |
Family
ID=74936722
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011363689.6A Active CN112487641B (zh) | 2020-11-27 | 2020-11-27 | 基于多尺度建模的煤化学链气化工艺设计方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112487641B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114722736A (zh) * | 2022-04-08 | 2022-07-08 | 青岛科技大学 | 一种基于实验与多尺度模拟的煤化学链气化全面优化方法 |
| CN115508401A (zh) * | 2021-06-23 | 2022-12-23 | 中国石油化工股份有限公司 | 分解放热性化学品热稳定微观特征的量化表征方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002356681A (ja) * | 2001-05-30 | 2002-12-13 | Mitsubishi Materials Corp | 石炭ガス化および燃焼をシミュレーションするためのモデリング方法およびシミュレーションシステム |
| CN101324908A (zh) * | 2008-08-01 | 2008-12-17 | 重庆大学 | 一种固体氧化物燃料电池的多尺度模拟方法 |
| US20150073188A1 (en) * | 2012-03-01 | 2015-03-12 | The Trustees Of Princeton University | Processes for producing synthetic hydrocarbons from coal, biomass, and natural gas |
| CN105138809A (zh) * | 2015-10-21 | 2015-12-09 | 中国石油大学(华东) | 一种无定形态烟煤结构孔隙特性的分析方法 |
| CN106960070A (zh) * | 2016-12-28 | 2017-07-18 | 山东科技大学 | 一种基于有限元‑离散元ct重构煤体的渗流模拟方法 |
| CN107325846A (zh) * | 2017-06-28 | 2017-11-07 | 大连理工大学 | 一种基于低阶煤梯级利用的煤热解化学链气化耦合工艺 |
| US20190003704A1 (en) * | 2015-12-22 | 2019-01-03 | Institutt For Energiteknikk | Sustainable Oxygen Carriers for Chemical Looping Combustion with Oxygen Uncoupling and Methods for Their Manufacture |
| CN111777035A (zh) * | 2020-06-01 | 2020-10-16 | 东南大学 | 一种基于铜铁基复合载氧体的自热式褐煤化学链气化方法 |
-
2020
- 2020-11-27 CN CN202011363689.6A patent/CN112487641B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002356681A (ja) * | 2001-05-30 | 2002-12-13 | Mitsubishi Materials Corp | 石炭ガス化および燃焼をシミュレーションするためのモデリング方法およびシミュレーションシステム |
| CN101324908A (zh) * | 2008-08-01 | 2008-12-17 | 重庆大学 | 一种固体氧化物燃料电池的多尺度模拟方法 |
| US20150073188A1 (en) * | 2012-03-01 | 2015-03-12 | The Trustees Of Princeton University | Processes for producing synthetic hydrocarbons from coal, biomass, and natural gas |
| CN105138809A (zh) * | 2015-10-21 | 2015-12-09 | 中国石油大学(华东) | 一种无定形态烟煤结构孔隙特性的分析方法 |
| US20190003704A1 (en) * | 2015-12-22 | 2019-01-03 | Institutt For Energiteknikk | Sustainable Oxygen Carriers for Chemical Looping Combustion with Oxygen Uncoupling and Methods for Their Manufacture |
| CN106960070A (zh) * | 2016-12-28 | 2017-07-18 | 山东科技大学 | 一种基于有限元‑离散元ct重构煤体的渗流模拟方法 |
| CN107325846A (zh) * | 2017-06-28 | 2017-11-07 | 大连理工大学 | 一种基于低阶煤梯级利用的煤热解化学链气化耦合工艺 |
| CN111777035A (zh) * | 2020-06-01 | 2020-10-16 | 东南大学 | 一种基于铜铁基复合载氧体的自热式褐煤化学链气化方法 |
Non-Patent Citations (4)
| Title |
|---|
| KARTIKEYA MAHALATKAR ET AL.: ""CFD simulation of a chemical-looping fuel reactor utilizing solid fuel"", 《CHEMICAL ENGINEERING SCIENCE》 * |
| KAVITHA G. MENON ET AL.: ""CFD simulation of fuel reactor for chemical looping combustion of Indian coal"", 《FUEL》 * |
| ZHE CUI ET AL.: ""Multi-scale Cleaner Analysis of Chemical Looping Gasification Coupled Coal Pyrolysis Process"", 《CHEMICAL ENGINEERING TRANSACTIONS》 * |
| 钱亚男: ""外热式内构件固定床/移动床煤热解反应器的数值模拟"", 《中国博士学位论文全文数据库 工程科技I辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115508401A (zh) * | 2021-06-23 | 2022-12-23 | 中国石油化工股份有限公司 | 分解放热性化学品热稳定微观特征的量化表征方法 |
| CN114722736A (zh) * | 2022-04-08 | 2022-07-08 | 青岛科技大学 | 一种基于实验与多尺度模拟的煤化学链气化全面优化方法 |
| CN114722736B (zh) * | 2022-04-08 | 2024-08-23 | 青岛科技大学 | 一种基于实验与多尺度模拟的煤化学链气化全面优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112487641B (zh) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Walspurger et al. | Sorption enhanced methanation for substitute natural gas production: Experimental results and thermodynamic considerations | |
| Yang et al. | Development of a full CO2 capture process based on the rectisol wash technology | |
| Yu et al. | Multi-objective optimization of ANN-based PSA model for hydrogen purification from steam-methane reforming gas | |
| CN112487641B (zh) | 基于多尺度建模的煤化学链气化工艺设计方法及其应用 | |
| Guo et al. | Thermodynamic and economic analysis of polygeneration system integrating atmospheric pressure coal pyrolysis technology with circulating fluidized bed power plant | |
| Yang et al. | Progress in coal chemical technologies of China | |
| Liu et al. | A proposed coal-to-methanol process with CO2 capture combined Organic Rankine Cycle (ORC) for waste heat recovery | |
| Wang | Systems engineering | |
| Srinivas et al. | Thermodynamic equilibrium model and exergy analysis of a biomass gasifier | |
| Al-Zareer et al. | Influence of selected gasification parameters on syngas composition from biomass gasification | |
| Yang et al. | Optimal design and exergy analysis of biomass-to-ethylene glycol process | |
| Xie et al. | Energy consumption and exergy analysis of MEA-based and hydrate-based CO2 separation | |
| CN102175345A (zh) | 多喷嘴对置式水煤浆气化炉炉膛温度的软测量方法 | |
| Liu et al. | Conceptual design of the coal to synthetic natural gas (SNG) process based on BGL gasifier: modeling and techno-economic analysis | |
| Yi et al. | CO2 capture and use in a novel coal-based polygeneration system | |
| Yu et al. | Design and economic evaluation of a coal-to-synthetic natural gas process | |
| Cui et al. | Multi-scale modeling and control of chemical looping gasification coupled coal pyrolysis system for cleaner production of synthesis gas | |
| Cao et al. | Performance modeling of integrated chemical looping air separation and IGCC with CO2 capture | |
| Yang et al. | Multi-aspect comparative analyses of two innovative methanol and power cogeneration systems from two different sources | |
| Zhu et al. | Applying the wavelet neural network to estimate hydrogen dissolution in underground sodium chloride solutions | |
| Adu et al. | Optimization and energy assessment of technological process for CO2 capture system of natural gas and coal combustion | |
| Zhang et al. | Machine learning assisted rediscovery of methane storage and separation in porous carbon from material literature | |
| Wodołażski et al. | Modelling and process integration study of dimethyl ether synthesis from syngas derived from biomass gasification: Flowsheet simulation | |
| Zhang et al. | Advanced process integration and machine learning-based optimization to enhance techno-economic-environmental performance of CO2 capture and conversion to methanol | |
| Tian et al. | Molecular dynamics and steady-state process simulation coupled research on CCLG system design for cleaner production of EG |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |