CN109255461B - 一种氢气资源的优化方法及优化系统 - Google Patents
一种氢气资源的优化方法及优化系统 Download PDFInfo
- Publication number
- CN109255461B CN109255461B CN201710563606.XA CN201710563606A CN109255461B CN 109255461 B CN109255461 B CN 109255461B CN 201710563606 A CN201710563606 A CN 201710563606A CN 109255461 B CN109255461 B CN 109255461B
- Authority
- CN
- China
- Prior art keywords
- hydrogen
- model
- optimization
- unit
- module
- 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.)
- Active
Links
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 615
- 239000001257 hydrogen Substances 0.000 title claims abstract description 615
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 title claims abstract description 597
- 238000005457 optimization Methods 0.000 title claims abstract description 238
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000011084 recovery Methods 0.000 claims abstract description 156
- 238000004364 calculation method Methods 0.000 claims abstract description 76
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims abstract description 71
- 238000004088 simulation Methods 0.000 claims abstract description 64
- 239000007789 gas Substances 0.000 claims description 62
- 238000004519 manufacturing process Methods 0.000 claims description 50
- 238000001179 sorption measurement Methods 0.000 claims description 46
- 238000013461 design Methods 0.000 claims description 40
- 239000012528 membrane Substances 0.000 claims description 38
- 239000000203 mixture Substances 0.000 claims description 36
- 238000000926 separation method Methods 0.000 claims description 35
- 238000004422 calculation algorithm Methods 0.000 claims description 30
- 239000000463 material Substances 0.000 claims description 30
- 238000004517 catalytic hydrocracking Methods 0.000 claims description 21
- 239000012535 impurity Substances 0.000 claims description 19
- 238000011160 research Methods 0.000 claims description 19
- 150000002431 hydrogen Chemical class 0.000 claims description 18
- 239000006227 byproduct Substances 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 17
- 229930195733 hydrocarbon Natural products 0.000 claims description 16
- 150000002430 hydrocarbons Chemical class 0.000 claims description 16
- 239000004215 Carbon black (E152) Substances 0.000 claims description 15
- 239000000047 product Substances 0.000 claims description 15
- 238000002407 reforming Methods 0.000 claims description 13
- 239000003054 catalyst Substances 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 12
- 239000003463 adsorbent Substances 0.000 claims description 9
- 238000006477 desulfuration reaction Methods 0.000 claims description 9
- 230000023556 desulfurization Effects 0.000 claims description 9
- 238000004230 steam cracking Methods 0.000 claims description 9
- 238000004523 catalytic cracking Methods 0.000 claims description 6
- 239000003245 coal Substances 0.000 claims description 6
- 238000004939 coking Methods 0.000 claims description 6
- 238000002485 combustion reaction Methods 0.000 claims description 6
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 claims description 5
- 239000005977 Ethylene Substances 0.000 claims description 5
- 238000005868 electrolysis reaction Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 description 32
- 239000003921 oil Substances 0.000 description 26
- 238000005984 hydrogenation reaction Methods 0.000 description 23
- 238000012545 processing Methods 0.000 description 18
- 239000000126 substance Substances 0.000 description 15
- 230000008901 benefit Effects 0.000 description 11
- 239000012071 phase Substances 0.000 description 9
- 150000001875 compounds Chemical class 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 239000003350 kerosene Substances 0.000 description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 239000002994 raw material Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 239000003502 gasoline Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 150000001336 alkenes Chemical class 0.000 description 2
- 125000003118 aryl group Chemical group 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 239000002283 diesel fuel Substances 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- JRZJOMJEPLMPRA-UHFFFAOYSA-N olefin Natural products CCCCCCCC=C JRZJOMJEPLMPRA-UHFFFAOYSA-N 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 238000012946 outsourcing Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000012466 permeate Substances 0.000 description 2
- WGACMNAUEGCUHG-VYBOCCTBSA-N (2s)-2-[[(2s)-2-[[(2s)-2-acetamidopropanoyl]amino]propanoyl]amino]-n-[(2s)-6-amino-1-[[(2s)-1-[(2s)-2-[[(2s)-1-[[(2s)-5-amino-1-[[(2s)-1-[[(2s)-1-[[(2s)-6-amino-1-[[(2s)-1-amino-3-(4-hydroxyphenyl)-1-oxopropan-2-yl]amino]-1-oxohexan-2-yl]amino]-3-hydroxy- Chemical compound CC(=O)N[C@@H](C)C(=O)N[C@@H](C)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CCCCN)C(=O)N[C@@H](CCCN=C(N)N)C(=O)N1CCC[C@H]1C(=O)N[C@@H](CO)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CCCN=C(N)N)C(=O)N[C@@H](CO)C(=O)N[C@@H](CCCCN)C(=O)N[C@H](C(N)=O)CC1=CC=C(O)C=C1 WGACMNAUEGCUHG-VYBOCCTBSA-N 0.000 description 1
- 230000003197 catalytic effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000004821 distillation Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- -1 flow rate Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 108010074544 myelin peptide amide-12 Proteins 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 239000012465 retentate Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
本发明公开了一种氢气资源的优化方法及优化系统。该方法包括以下步骤：划分氢气系统的各单元；建立各单元模型，再建立整个氢系统的非线性的数学模拟模型；对非线性的数学模拟模型进行模拟求解；确定整个氢系统的目标函数；建立用氢单元数学规划模型；判断氢回收单元各流股是否有作为用氢单元的供氢；建立氢回收单元数学规划模型；确定供氢单元的氢气公用工程量；判断整个氢系统的目标函数是否达到计算条件。本发明方法能够避免复杂系统优化求解常遇到的难收敛、局部最优等问题，可在复杂氢气系统下有效提高氢气系统优化求解的计算效率，获得氢气系统最优解，得到最佳匹配方案，降低炼厂氢耗、提高氢气回收效率、减少氢公用工程量。
Description
技术领域
本发明涉及氢气资源优化领域，具体涉及一种炼厂氢气资源的优化方法及优化系统。
背景技术
随着加工原油的重质化、劣质化趋势加剧及产品质量升级要求加强，炼厂对氢气的需求量越来越大，氢气已成为炼油厂原料成本中仅次于原油成本的第二大成本要素。氢气资源优化利用技术是企业实现降本增效的重要手段，已受到企业愈来愈多的重视。
目前，氢气资源优化研究主要有夹点分析法及数学规划法，其特点如下：
（1）夹点分析法。开展氢气网络夹点分析，首先需得到“氢气体积分数-流股流量”关系的氢源复合曲线及氢阱复合曲线，通过比较氢源复合曲线及氢阱复合曲线在同一氢气纯度下的面积得到氢气剩余曲线图，即为氢气网络的夹点图。氢夹点分析技术优点是简单直观，可以在氢网络设计前得到氢气公用工程量最小需求目标。缺点是夹点分析法仅考虑流股的流量、氢纯度，未考虑压力影响，造成结果是优化的结果可能存在过度优化而无法在实际生产中实施，另一缺点是夹点分析法对于大规模实际问题使用起来比较困难、繁琐；
（2）数学规划法。基本原理是将氢气系统优化过程转化成数学建模方程的求解过程，优化目的即为目标函数，根据炼厂实际情况设置约束条件，最后利用易于收敛的求算方法获取目标函数的最优解。优点是处理复杂氢网络优化问题能省却图解法夹点分析的繁琐计算过程。缺点：一是由于氢气系统涉及面广，联立求解方程组较多，因而优化求解的难度较大，此外，诸多设备模型存在非凸非连续的特征，且包含很强的约束条件，造成了优化求解的难度，不容易得到最优解或陷入局域最优解；二是复杂系统下的约束条件不易设定，约束条件增多，计算量增大，优化求解困难，有时难以获得最优解，反之，约束条件减少，优化求解结果的可实施性降低。
例如，CN106485341A涉及一种炼油厂氢气系统优化方法。该方法将炼油厂氢气系统收益最大作为目标，对约束条件进行优选，形成了包括轻烃资源优化在内的新的数学规划算法；与现有氢气系统的数学规划算法相比，更切合实际，可获得更大收益以满足炼油厂节氢、增效的需求。缺点：一将整个氢气系统作为优化研究对象，涉及面广，约束限制多，造成优化求解的计算量偏大，有时不易获得系统最优解；二是在优化求解过程中，需要反复调整待优化参数以确保最大收益函数的解趋于最大值，这种在复杂氢气系统中反复试差调整待优化参数的做法必定会降低优化求解的计算效率。
因此，如何将现有技术问题加以解决，提供一种新的高效准确的氢气资源优化研究方法，即为本领域技术人员的研究方向所在。
发明内容
针对现有技术中的缺陷，本发明提供了一种氢气资源的优化方法及优化系统。本发明方法能够避免复杂系统优化求解常遇到的难收敛、局部最优等问题，可在复杂氢气系统下有效提高氢气系统优化求解的计算效率，获得氢气系统最优解，得到最佳匹配方案，降低炼厂氢耗、提高氢气回收效率、减少氢公用工程量。
本发明提供了一种氢气资源的优化方法，包括以下步骤：
步骤S1：将氢气系统划分为供氢单元、用氢单元、氢回收单元，确定各单元的设计参数及操作参数；
步骤S2：分别建立供氢单元模型、用氢单元模型、氢回收单元模型、瓦斯系统模型，各模型间按实际物料关系建立连接，进而建立整个氢系统的非线性的数学模拟模型；
步骤S3：对非线性的数学模拟模型进行模拟求解，得到模拟运算结果；
步骤S4：确定整个氢系统的目标函数；
步骤S5：建立用氢单元数学规划模型，然后优化求解；其中，所述用氢单元数学规划模型的优化求解可将步骤S3模拟计算结果作为优化运算的初始可行解，在优化变量的取值范围内，采用优化算法进行求解，得到求解优化结果；
步骤S6：判断氢回收单元各流股是否有作为用氢单元的供氢，如有作为用氢单元的供氢，则返回步骤S5；未有作为用氢单元的供氢，则执行步骤S7；
步骤S7：建立氢回收单元数学规划模型，优化求解；所述氢回收单元优化求解是将步骤S3模拟结果为优化运算初始可行解，步骤S5用氢单元优化结果为进料条件，在优化变量取值范围内，采用优化算法进行求解，得到求解优化结果；
步骤S8：确定供氢单元的氢气公用工程量；
步骤S9：判断整个氢系统的目标函数是否达到计算条件；达到计算条件，则将优化求解结果作为氢气系统模型的运行参数，未达到计算条件，则返回步骤S5。
在所述步骤S1中，所述供氢单元包括水蒸汽裂解制氢、煤制氢、重整副产氢、乙烯装置副产氢、离子膜电解水副产氢、焦化装置、催化裂化装置中的一种或几种装置，所述用氢单元包括加氢精制装置、加氢裂化装置中的一种或几种装置；所述氢回收单元包括变压吸附装置、膜分离装置中的一种或几种装置；
所述设计参数包括装置的设计数据、催化剂及吸附剂的物化性质；所述操作参数包括装置的操作条件，物料平衡，以及流股的流量、压力、组成。
在所述步骤S2中，所述供氢单元模型包括简易供氢模型或详细供氢模型；其中，详细供氢模型包括水蒸汽裂解制氢模型、煤制氢模型，所述详细供氢模型为考虑氢气生产过程中存在的各个反应器的生产工艺过程模型，用于模拟动力学、热效应及催化剂活性，基于变化的进料组成、流量预测产物收率、组成；所述简易供氢模型包括重整副产氢模型及其他供氢装置模型，所述简易供氢模型为仅考虑装置的供氢流量、组成的模型；
所述用氢单元模型包括加氢精制处理模型和加氢裂化模型；
所述氢回收单元模型包括变压吸附装置模型、膜分离装置模型、轻烃回收装置、脱硫装置简易模型；
所述瓦斯系统模型采用均相管网模型，用于描述氢气系统中部分含氢流股排放至瓦斯管线或直送至燃烧炉的状况。
在所述步骤S4中，氢气系统目标目标函数为年总成本最小、总操作成本最少或氢气公用工程量最小，优选为氢气公用工程量最小。
在所述步骤S5中，在用氢单元数学规划模型进行优化求解时，以新氢流量、循环氢流量为优化变量，以最小化总用氢量为目标函数值，以满足装置最小用氢需求、设备运行的要求为约束条件；所述用氢单元的优化求解采用优化算法进行求解，得到求解优化结果，所述优化算法优选为序列二次规划算法或序贯线性规划算法；所述用氢单元优化求解结果通过物料连接关系传递给下游装置模型。
在所述步骤S6中，含氢流股是否可作为用氢单元新氢的标准为：含氢流股氢纯度大于用氢装置反应器入口氢纯度，CH4体积分数小于CH4设定值，气体杂质含量小于气体杂质设定值， CH4设定值和气体杂质设定值需根据装置工艺要求进行设定。
在所述步骤S7中，氢回收单元采用基于协同优化理论的建模方法，将整个氢回收系统作为系统级优化研究对象，从而建立整体系统优化模型；将变压吸附装置、膜分离装置作为氢回收系统下的各子系统优化研究对象，从而建立子系统优化模型；所述系统级优化对象与各子系统通过欲回收的含氢流股相耦合。
（a）整体系统优化模型如下：
整体系统优化模型的目标函数：
式中，F表示氢气回收系统的总回收效益；F i (X i )表示各氢回收装置的回收效益；X i 表示模型决策变量；
整体系统优化模型的约束条件：
氢气回收系统与子系统i等式一致性约束为：
……
……
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量，ε表示预设松弛因子，N表示变压吸附单元或膜分离单元中包含的变压吸附装置或膜分离装置的数量；
（b）子系统优化模型如下：
子系统优化模型的目标函数：
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量；
子系统优化模型的约束条件：
式中： 
、
——子系统i的等式不等式约束。
在所述步骤S8中，由用氢单元、氢回收单元优化结果经物料衡算确定氢气公用工程量。
在所述步骤S9中，所述达到计算条件：若两次目标函数计算相差满足收敛要求，则结束系统计算；若两次目标函数计算相差未达到收敛要求，且迭代次数仍未达到设计上限，则返回步骤S5再进行优化计算；若两次目标函数计算相差未满足收敛要求，但迭代次数达到上限，则停止计算。
本发明还提供了一种氢气资源的优化系统，包括以下模块：划分模块、数学模拟模型模块、模拟求解模块、第一确定模块、用氢单元数学规划模型模块、第一判断模块、氢回收单元数学规划模型模块、第二确定模块和第二判断模块；
划分模块，用于将氢气系统划分为供氢单元、用氢单元、氢回收单元，确定各单元的设计参数及操作参数；
数学模拟模型模块，用于分别建立供氢单元模型、用氢单元模型、氢回收单元模型、瓦斯系统模型，各模型间按实际物料关系建立连接，进而建立整个氢系统的非线性的数学模拟模型；
模拟求解模块，用于对非线性的数学模拟模型进行模拟求解，得到模拟运算结果；
第一确定模块，用于确定整个氢系统的目标函数；
用氢单元数学规划模型模块，用于建立用氢单元数学规划模型，然后优化求解；其中，所述用氢单元数学规划模型的优化求解可将模拟求解模块的模拟计算结果作为优化运算的初始可行解，在优化变量的取值范围内，采用优化算法进行求解，得到求解优化结果；
第一判断模块，用于判断氢回收单元各流股是否有作为用氢单元的供氢，如有作为用氢单元的供氢，则返回用氢单元数学规划模型模块；未有作为用氢单元的供氢，则执行氢回收单元数学规划模型模块；
氢回收单元数学规划模型模块，用于建立氢回收单元数学规划模型，优化求解；所述氢回收单元优化求解是将模拟求解模块的模拟结果为优化运算初始可行解，以用氢单元数学规划模型模块的用氢单元优化结果为进料条件，在优化变量取值范围内，采用优化算法进行求解，得到求解优化结果；
第二确定模块，用于确定供氢单元的氢气公用工程量；
第二判断模块，用于判断整个氢系统的目标函数是否达到计算条件；达到计算条件，则将优化求解结果作为氢气系统模型的运行参数，未达到计算条件，则返回用氢单元数学规划模型模块。
在所述划分模块中，所述供氢单元包括水蒸汽裂解制氢、煤制氢、重整副产氢、乙烯装置副产氢、离子膜电解水副产氢、焦化装置、催化裂化装置中的一种或几种装置，所述用氢单元包括加氢精制装置、加氢裂化装置中的一种或几种装置；所述氢回收单元包括变压吸附装置、膜分离装置中的一种或几种装置；
所述设计参数包括装置的设计数据、催化剂及吸附剂的物化性质；所述操作参数包括装置的操作条件，物料平衡，以及流股的流量、压力、组成。
所述数学模拟模型模块中，所述供氢单元模型包括简易供氢模型或详细供氢模型；其中，详细供氢模型包括水蒸汽裂解制氢模型、煤制氢模型，所述详细供氢模型为考虑氢气生产过程中存在的各个反应器的生产工艺过程模型，用于模拟动力学、热效应及催化剂活性，基于变化的进料组成、流量预测产物收率、组成；所述简易供氢模型包括重整副产氢模型及其他供氢装置模型，所述简易供氢模型为仅考虑装置的供氢流量、组成的模型；
所述用氢单元模型包括加氢精制处理模型和加氢裂化模型；
所述氢回收单元模型包括变压吸附装置模型、膜分离装置模型、轻烃回收装置、脱硫装置简易模型；
所述瓦斯系统模型采用均相管网模型，用于描述氢气系统中部分含氢流股排放至瓦斯管线或直送至燃烧炉的状况。
在所述第一确定模块中，氢气系统的目标目标函数为年总成本最小、总操作成本最少或氢气公用工程量最小，优选为氢气公用工程量最小。
所述用氢单元数学规划模型模块中，在用氢单元数学规划模型进行优化求解时，以新氢流量、循环氢流量为优化变量，以最小化总用氢量为目标函数值，以满足装置最小用氢需求、设备运行的要求为约束条件；所述用氢单元的优化求解采用优化算法进行求解，得到求解优化结果，所述优化算法优选为序列二次规划算法或序贯线性规划算法；所述用氢单元优化求解结果通过物料连接关系传递给下游装置模型。
所述第一判断模块中，含氢流股是否可作为用氢单元新氢的标准为：含氢流股氢纯度大于用氢装置反应器入口氢纯度，CH4体积分数小于CH4设定值，气体杂质含量小于气体杂质设定值， CH4设定值和气体杂质设定值需根据装置工艺要求进行设定。
所述氢回收单元数学规划模型模块中，氢回收单元采用基于协同优化理论的建模方法，将整个氢回收系统作为系统级优化研究对象，从而建立整体系统优化模型；将变压吸附装置、膜分离装置作为氢回收系统下的各子系统优化研究对象，从而建立子系统优化模型；所述系统级优化对象与各子系统通过欲回收的含氢流股相耦合。
（a）整体系统优化模型如下：
整体系统优化模型的目标函数：
式中，F表示氢气回收系统的总回收效益；F i (X i )表示各氢回收装置的回收效益；X i 表示模型决策变量；
整体系统优化模型的约束条件：
氢气回收系统与子系统i等式一致性约束为：
……
……
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量，ε表示预设松弛因子，N表示变压吸附单元或膜分离单元中包含的变压吸附装置或膜分离装置的数量；
（b）子系统优化模型如下：
子系统优化模型的目标函数：
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量；
子系统优化模型的约束条件：
式中： 
、
——子系统i的等式不等式约束。
所述第二确定模块中，由用氢单元、氢回收单元优化结果经物料衡算确定氢气公用工程量。
第二判断模块，所述达到计算条件：若两次目标函数计算相差满足收敛要求，则结束系统计算；若两次目标函数计算相差未达到收敛要求，且迭代次数仍未达到设计上限，则返回用氢单元数学规划模型模块再进行优化计算；若两次目标函数计算相差未满足收敛要求，但迭代次数达到上限，则停止计算。
与现有技术相比，本发明的氢气资源的优化方法及系统具有如下优点：
（1）本发明的氢气资源的优化方法，先建立整个氢系统的非线性的数学模拟模型、确定整个氢系统目标函数，将整个氢气系统的数学规划模型的求解问题，简化为具有序贯连接关系的各子单元优化求解，然后在各子单元优化求解过程中与整个氢系统非线性的数学模拟模型建立联系，从而不再简单地将整个氢气系统作为优化研究对象，而是充分研究考虑了氢气系统内部单元的逻辑关联，将整个氢气系统的优化研究问题简化为耗氢单元、氢回收单元、供氢单元各子系统的优化研究，避免了复杂系统优化求解常遇到的难收敛、局部最优等问题，可在氢气大系统下有效提高氢气系统优化求解的计算效率，获得氢气系统最优解，得到最佳匹配方案、操作方案，显著降低炼厂氢耗、提高氢气回收效率、减少氢公用工程量，进而提高企业经济效益，提升企业综合实力。
（2）本发明从系统角度出发，考虑了供氢单元、耗氢单元、氢回收单元间及各单元内部装置间的协同优化，充分发挥了整个氢气系统优化节氢的潜力，可有效提高氢气系统的操作水平，最小化装置用氢，最大化氢回收，提高氢气系统运行效率，降低氢气公用工程消耗。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明一实施例提供的炼厂氢气资源的优化方法的流程图；
图2为一炼厂氢气系统的氢气管网的示意图；
图3为蒸汽裂解制氢装置模型结构的示意图；
图4为炼厂氢系统在优化后氢气管网的示意图；
图5是本发明另一实施例提供的炼厂氢气资源的优化系统的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明一实施例提供了一种炼厂氢气资源的优化方法，参见图1所示的流程图，该方法包括以下步骤：
步骤101：将氢气系统划分为供氢单元、用氢单元、氢回收单元，确定各单元的设计参数及操作参数。
在本步骤中，供氢单元主要涉及水蒸汽裂解制氢、煤制氢、重整副产氢、乙烯装置副产氢、离子膜电解水副产氢，可以是其中的一类装置或几类装置。供氢单元还可以包括焦化装置、催化裂化装置。
用氢单元主要涉及汽油加氢、航煤加氢、柴油加氢、蜡油加氢、渣油加氢、润滑油加氢等各类炼油加氢精制装置及加氢裂化装置，可以是其中的一类装置或几类装置。氢回收单元主要涉及变压吸附装置、膜分离装置，可以是其中一类装置或几类装置。氢回收装置还包括与上述氢回收装置相关联的轻烃回收装置、脱硫装置。
设计参数主要涉及装置的基础设计数据、催化剂及吸附剂的基础物化性质。操作参数主要涉及装置的操作条件、物料平衡及流股的流量、压力、组成。
例如，图2为一炼厂氢气系统示意图，表1为炼厂氢气系统划分的三单元装置构成，表2为部分装置设计参数及操作参数汇总。
表1 氢气系统三单元装置构成
| 项目 | 装置 |
| 供氢单元 | 1#天然气制氢装置、2#天然气制氢装置、连续重整装置、外购高纯氢、延迟焦化、催化裂化 |
| 用氢单元 | 1#、2#柴油加氢装置、航煤加氢装置、Szorb装置、加氢裂化装置、蜡油加氢装置 |
| 氢回收单元 | 1#制氢装置PSA、2#制氢装置PSA、连续重整氢PSA、轻烃回收装置 |
表2 部分装置设计参数及操作参数汇总
| 蜡油加氢 | 加氢裂化 | 2#柴油加氢 | 1#柴油加氢 | 航煤加氢 | S Zorb | |
| 氢油比 | ≥500 | ≥700 | ≥500 | ≥450 | 约45 | mol≥0.26 |
| 氢分压，MPa | ≥8 | ≥12.0（循环氢纯度≥80%） | ≥6.4 | 入口压力7.0 | 入口压力2.0 | 0.61 |
| 反应器设计/操作压力，MPa | 12/11.0 | 15.7/14.0 | 9.3/9 | 8.0/7.1 | 2.3/2.0 | 2.38/2.7 |
| 新氢压缩机入/出口压力，MPa | 2→12 | 2→15.7 | 2→9.3 | 2→7.2 | 1.4→2.3 | 1.8→3.5 |
| 新氢压缩机设计/最大气量，Nm<sup>3</sup>/h | 34000/40000 | 38500/45000 | 27835/33000 | 10000/15000 | 1500/10000 | 3500/8000 |
| 新氢压缩机类型 | 容积式 | 往复式 | 往复式 | 往复式 | 往复式 | 往复式 |
| 循环氢入/出口压力，MPa | 9.5→12 | 13.3→15.3 | 7.1→9.1 | 6.59→7.83 | - | 2.2→3.5 |
| 压缩机设计/最大气量，Nm<sup>3</sup>/h | 200000/250000 | 370000/480000 | 207000/250000 | 40000/50000 | - | - |
| 压缩机类型 | 离心式 | 离心式 | 离心式 | 离心式 | - | - |
| 热高分设计/操作温度，℃ | - | 280/225 | 280/220 | 280/230 | /129 | |
| 热高分设计/操作压力，MPa | - | 15/13 | 8.8/8.1 | 6.8/6.5 | - | /2.72 |
| 热低设计/操作温度，℃ | - | 300/240 | 280/220 | 280/225 | - | - |
| 热低分设计/操作压力，MPa | - | 3.3/2.7 | 3.0/2.4 | 2.5/2.2 | - | - |
| 冷高分设计/操作温度，℃ | 140/49 | 150/55 | 150/50 | 150/50 | 150/45 | - |
| 冷高分设计/操作压力，MPa | 11.5/9.5 | 14.9/13 | 8.65/7.45 | 6.8/6.3 | 2.5/1.95 | - |
| 冷低分设计/操作温度，℃ | 100/51 | 170/40 | 150/50 | 150/50 | 150/40 | - |
| 冷低分设计/操作压力，MPa | 2.5/2 | 3.2/2.6 | 3/2.4 | 3/2.2 | 2.5/1.7 | - |
注：S Zorb的氢油比按摩尔氢油比，其它均按体积氢油比。
步骤102：分别建立供氢单元模型、用氢单元模型、氢回收单元模型、瓦斯系统模型，各模型间按实际物料关系建立连接，进而建立整个氢系统的非线性的数学模拟模型。
（1）供氢单元模型
在步骤102中，供氢单元模型包括简易供氢模型或详细供氢模型；其中，详细供氢模型包括水蒸汽裂解制氢模型、煤制氢模型，详细供氢模型为考虑氢气生产过程中存在的各个反应器的生产工艺过程模型，用于模拟动力学、热效应及催化剂活性，基于变化的进料组成、流量预测产物收率、组成。
简易供氢模型包括重整副产氢模型、蒸汽裂解制氢装置模型及其他供氢装置模型，简易供氢模型为仅考虑装置的供氢流量、组成的模型，不涉及产氢工艺过程的模拟。
蒸汽裂解制氢装置模型结构示意如图3所示，转化炉管内反应可描述如下：
上述反应对应的动力学方程分别为：
式中：
中温变换反应器中反应可如下描述：
（2）用氢单元模型
用氢单元模型为严格反应动力学模型，精确描述原料经反应器后的油品质量变化及轻端组分变化过程。油品质量主要通过油品馏程数据、相对密度表示。轻端组分包括H2、NH3、H2S、N2、O2、H2O、CO、CO2、C1~C5轻烃等。用氢单元模型能模拟当前工况下的生产状况，也能预测反应器进料性质变化或工艺条件变化后的产物分布情况。工艺条件主要涉及反应器温度、压力、氢油比及高低压分离器的温度、压力。
用氢单元模型仅涉及新氢压缩机、反应器、高低压分离器、循环氢压缩机模拟计算，不涉及后续分离工艺装置的核算。反应器入口进料形式为新氢、循环氢、原料油混合物，进料性质主要涉及气相物料的组成、流量、杂质含量及原料油性质。反应器出口为转化后的生成油与轻端组分的混合物，经高低压分离器分离成循环氢、低分气、干气及液相油品。
用氢单元模型可以包括加氢精制处理模型和加氢裂化模型。其中，加氢精制处理模型主要涉及汽油加氢、航煤加氢等各类别加氢精制处理装置模型，加氢精制处理模型要呈现以下几种转化反应，包括加氢脱硫反应、加氢脱氮反应、加氢脱氧反应、烯烃加氢饱和反应及芳烃加氢饱和反应。加氢裂化模型主要涉及各类别加氢裂化装置模型，加氢裂化模型要呈现以下几种转化反应，包括加氢脱硫反应、加氢脱氮反应、烃类裂化反应。
下面以加氢精制处理模型、加氢裂化模型为例进行详细阐述：
对于加氢精制处理模型，以测定给定加氢处理反应器给定反应机理i的反应动力学公式可表示如下：
式中：
用以上方程式分别对加氢脱硫反应、加氢脱氮反应、加氢脱氧反应、烯烃加氢饱和反应及芳烃加氢饱和反应求解，对各反应氢耗速率求和，即可求得该装置总化学氢耗速率。
对于加氢裂化装置模型，除必须考虑上述的脱硫、脱氮等反应外，还必须考虑烃类裂化消耗的氢。模型将反应器进出口油品按每15℃划分虚拟组分，采用下述反应动力学模型测定其相应参数，进而建立加氢裂化装置的反应动力学模型。
式中：
对于轻烃组分的预测，可采用下式计算，
式中：
轻烃组分分布如表3所示。
表3 用氢装置反应产生轻烃组成分布
| 组成 | 摩尔分率 |
| CH<sub>4</sub> | x<sub>1</sub> |
| C<sub>2</sub>H<sub>6</sub> | x<sub>2</sub> |
| C<sub>2</sub>H<sub>4</sub> | x<sub>3</sub> |
| C<sub>3</sub>H<sub>8</sub> | x<sub>4</sub> |
| C<sub>3</sub>H<sub>6</sub> | x<sub>5</sub> |
| C<sub>4</sub>H<sub>10</sub> | x<sub>6</sub> |
| C<sub>4</sub>H<sub>8</sub> | x<sub>7</sub> |
| C<sub>5</sub>H<sub>12</sub> | x<sub>8</sub> |
| C<sub>5</sub>H<sub>10</sub> | x<sub>9</sub> |
| ∑ | 1.0 |
（3）氢回收单元模型
氢回收单元模型分为变压吸附装置模型、膜分离装置模型，此外还包括轻烃回收装置模型、脱硫装置简易模型。
本发明方法中，变压吸附装置模型要能够模拟当前工况下的分离状况，还要能预测进料流量、组成变化后的吸附分离情况。本发明方法中，在变压吸附装置模型中，从吸附床底至顶，按装填吸附剂类别不同划分不同的微元段，假设压力及温度恒定不变，流动模型采用轴向分散活塞流模型，由总物质平衡计算吸附引起的流速变化，传质速率方程采用线性推动力模型（LDF），吸附平衡方程采用扩展的兰格缪尔模型描述。变压吸附装置数学模型如下：
微元体积中气体组分i的传质平衡方程为：
其中，D L 表示床层轴向扩散系数；C i 表示组分i气相总浓度；
表示气流速度；ρ p 表示吸附压力P下气相密度；ε表示分子吸附床孔隙率；
表示组分i吸附平衡浓度，Z表示微元段轴向距离；t表示时间；
总传质平衡方程为：
其中，D L 表示床层轴向扩散系数；C表示床层气相浓度；ρ p 表示吸附压力P下气相密度，C i 表示组分i气相总浓度；v表示气流速度； 
表示组分i吸附平衡浓度；ε表示分子吸附床孔隙率；Z表示微元段轴向距离，t表示时间；
传质速率方程：
其中，
表示气固传质系数；
表示组分i吸附床气相浓度；
表示组分i吸附平衡浓度；
吸附平衡方程为：
其中，θ i 表示待吸附的混合气体中，某层吸附剂上气体组分i的覆盖率；p i 表示待吸附的混合气体中，气体组分i的分压；B i 表示气体组分i在该层吸附剂上的兰格缪尔吸附常数；q i 表示i组分的吸附量，q max,i 表示混合气分中i组分的最大吸附量。
本发明方法中，膜分离装置模型要能够模拟当前工况下的分离状况，还要能预测进料流量、组成变化后的吸附分离情况。本发明方法中，在膜分离装置模型中，忽略膜两侧流体的流动阻力，假设原料侧气体组成呈线性变化，渗透侧为全混合形式，组分i气体渗透量的数学模型如下：
式中：
本发明方法中，轻烃回收装置模型可以采用固定特定组分回收率的方法进行简化处理，回收率可人为设定。
本发明方法中，脱硫装置模型可以采用固定脱硫率、脱CO2率的方法进行简化处理，脱硫率、脱CO2率由人为设定。
瓦斯系统模型采用均相管网模型，用于描述氢气系统中部分含氢流股排放至瓦斯管线或直送至燃烧炉的状况。
本发明方法中，氢气系统模型由上述建立的供氢单元模型、用氢单元模型、氢回收单元模型及瓦斯系统模型，依据实际物流关系，采用物料连接线连接而成。物料连接线仅代表两物理装置物料的连接关系，优选的，可根据需要能进行管线压降、流速的计算。
步骤103：对非线性的数学模拟模型进行模拟求解，得到模拟运算结果。
步骤104：确定整个氢系统的目标函数。
本发明方法中，步骤104中氢气系统目标函数为年总成本最小、总操作成本最少或氢气公用工程量最小，优选的为氢气公用工程量最小为优化目标函数。其中，年总成本=设备年均投资费用+设备操作费用+外购成本+公用工程消耗-外卖收益-燃料价值。总操作成本=∑F i ×P i ，即总操作成本等于流入设备的气体总量×单位气体流量的操作费用。氢气公用工程量指水蒸汽裂解制氢及煤制氢的量。
以氢气系统氢气公用工程量最小为目标函数进行系统优化，氢系统物料平衡关系如下：
氢气公用工程量+副产氢+回收氢=用氢装置用氢+排放至瓦斯系统的氢
欲使氢气系统氢气公用工程量最小，那么必须满足以下条件，即用氢单元耗氢最小、氢回收单元回收氢气最大。
步骤105：建立用氢单元数学规划模型，然后优化求解；其中，所述用氢单元数学规划模型的优化求解可将步骤103模拟计算结果作为优化运算的初始可行解，在优化变量的取值范围内，采用优化算法进行求解，得到求解优化结果。
在步骤105中，用氢单元优化求解采用序列二次规划算法、序贯线性规划算法。用氢单元优化求解结果通过物料连接关系传递给下游装置模型。
在步骤105中，用氢单元数学优化模型的目标函数可以为总用氢最小，调节参数为新氢流量、循环氢流量（即氢分压、氢油比）参数，约束条件为满足装置最小用氢需求（最小氢油比、氢分压、装置耗氢）、关键设备运行要求。
例如，用氢装置优化需满足：反应器入口氢分压不小于工艺要求值，氢油比不小于工艺要求值，新氢中纯氢补给量不小于装置化学氢耗，循环氢氢纯度不小于工艺要求值，反应器进口需满足物料平衡。
式中：
步骤106：判断氢回收单元各流股是否有作为用氢单元的供氢，如有作为用氢单元的供氢，则返回步骤105；未有作为用氢单元的供氢，则执行步骤107；
在步骤106中，含氢流股是否可作为用氢单元新氢的标准如下：含氢流股氢纯度大于某用氢装置反应器入口氢纯度，CH4体积分数小于某一设定值，杂质含量（CO、O2、N2等）小于某一设定值，上述的设定值需根据装置工艺要求进行设定。
例如用氢单元新氢的标准可设定如下：
若含氢流股氢纯度大于88.0%，CH4体积分数小于4%，杂质含量（CO、O2、N2等）小于50ppm，可将该流股送至对氢气质量要求不太高的用氢装置，如汽油加氢、航煤加氢、Szorb等装置；若含氢流股氢纯度＞93%，CH4体积分数小于2%，杂质含量（CO、O2、N2等）小于50ppm，除可将该流股上述装置外，还可送至柴油加氢、蜡油加氢及部分加氢裂化装置。
步骤107：建立氢回收单元数学规划模型，优化求解；氢回收单元优化求解是将步骤103模拟结果为优化运算初始可行解，步骤105用氢单元优化结果为进料条件，在优化变量取值范围内，采用优化算法进行求解，得到求解优化结果。
在步骤107中，对于氢回收单元数学规划模型的建立，采用基于协同优化理论的建模方法，将整个氢回收系统作为系统级优化研究对象，将变压吸附装置、膜分离装置作为氢回收系统下的各子系统优化研究对象；系统级优化对象与各子系统通过欲回收的含氢流股相耦合。
本发明方法中，系统级优化求解、子系统优化求解采用常规的优化算法，可以是遗传算法、多目标优化求解法、既约梯度法等。
（a）整个氢回收系统级优化模型
对于整个氢回收系统，目标函数为最大化氢气回收，约束条件为各子系统优化目标函数，整体系统优化模型的目标函数：
式中，
表示氢气回收系统的总回收效益；
表示各氢回收装置的回收效益；
表示模型决策变量；
整体系统优化模型的约束条件：
氢气回收系统与子系统i等式一致性约束为：
……
……
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量，ε表示预设松弛因子。可以理解的是，这里的子系统i表示变压吸附单元或膜分离单元。N表示变压吸附单元或膜分离单元中包含的变压吸附装置或膜分离装置的数量。
（b）子系统优化模型
本发明方法中，子系统为各变压吸附装置、膜分离装置，子系统优化目标函数为各装置子系统变量与系统级变量偏差最小，约束条件为各系统的相应的等式及不等式约束，子系统i的优化模型如下：
子系统优化模型的目标函数：
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量。
子系统优化模型的约束条件：
式中： 
、
为子系统i的等式不等式约束。
步骤108：确定供氢单元的氢气公用工程量；
在步骤108中，须由用氢单元、氢回收单元优化结果经物料衡算确定氢气公用工程量。例如，供氢单元氢气公用工程量可以通过下述关系式确定：
氢公用工程＝用氢装置用氢＋排放至瓦斯系统的氢－副产氢－回收氢
步骤109：判断整个氢系统的目标函数是否达到计算条件；达到计算条件，则将优化求解结果作为氢气系统模型的运行参数，未达到计算条件，则返回步骤105。
本发明方法中，步骤109达到计算条件主要涉及判断目标函数计算是否收敛及迭代次数是否达到上限：若两次目标函数计算相差满足一定要求，则结束系统计算；若两次目标函数计算相差未达到设计要求，且迭代次数仍未达到设计上限，则返回步骤105再进行优化计算；若两次目标函数计算相差满足一定要求，但迭代次数达到上限，则停止计算。系统级迭代最大次数优选为5~50。
图4为炼厂氢系统在优化后氢气管网的示意图，表4为优化前后氢回收系统整体性能对比。从结果可看出，优化后系统多回收氢气6425Nm3/h，进而减少制氢装置制氢6425Nm3/h，减少了能源消耗及CO2排放，具有显著的经济效益。
表4 优化前后氢回收系统整体性能对比
| 优化前 | 优化后 | |
| 总供纯氢，Nm<sup>3</sup>/h | 82987 | 75970 |
| 装置总用氢，Nm<sup>3</sup>/h | 63138 | 62546 |
| 排放至瓦斯系统的纯氢，Nm<sup>3</sup>/h | 19849 | 13424 |
| PSA回收氢气，Nm<sup>3</sup>/h | - | 6425 |
| 系统氢气利用率，% | 76.08 | 82.33 |
如图5所示，本发明还提供了一种炼厂氢气资源的优化系统，包括以下模块：划分模块21、数学模拟模型模块22、模拟求解模块23、第一确定模块24、用氢单元数学规划模型模块25、第一判断模块26、氢回收单元数学规划模型模块27、第二确定模块28、第二判断模块29；
划分模块21：将氢气系统划分为供氢单元、用氢单元、氢回收单元，确定各单元的设计参数及操作参数；
数学模拟模型模块22：分别建立供氢单元模型、用氢单元模型、氢回收单元模型、瓦斯系统模型，各模型间按实际物料关系建立连接，进而建立整个氢系统的非线性的数学模拟模型；
模拟求解模块23：对非线性的数学模拟模型进行模拟求解，得到模拟运算结果；
第一确定模块24：确定整个氢系统的目标函数；
用氢单元数学规划模型模块25：建立用氢单元数学规划模型，然后优化求解；其中，用氢单元数学规划模型的优化求解可将模拟求解模块23的模拟计算结果作为优化运算的初始可行解，在优化变量的取值范围内，采用优化算法进行求解，得到求解优化结果；
第一判断模块26：判断氢回收单元各流股是否有作为用氢单元的供氢，如有作为用氢单元的供氢，则返回用氢单元数学规划模型模块25；未有作为用氢单元的供氢，则执行氢回收单元数学规划模型模块27；
氢回收单元数学规划模型模块27：建立氢回收单元数学规划模型，优化求解；氢回收单元优化求解是将模拟求解模块23的模拟结果为优化运算初始可行解，以用氢单元数学规划模型模块25的用氢单元优化结果为进料条件，在优化变量取值范围内，采用优化算法进行求解，得到求解优化结果；
第二确定模块28：确定供氢单元的氢气公用工程量；
第二判断模块29：判断整个氢系统的目标函数是否达到计算条件；达到计算条件，则将优化求解结果作为氢气系统模型的运行参数，未达到计算条件，则返回用氢单元数学规划模型模块25。
在划分模块21中，供氢单元包括水蒸汽裂解制氢、煤制氢、重整副产氢、乙烯装置副产氢、离子膜电解水副产氢、焦化装置、催化裂化装置中的一种或几种装置，用氢单元包括加氢精制装置、加氢裂化装置中的一种或几种装置；氢回收单元包括变压吸附装置、膜分离装置中的一种或几种装置；
设计参数包括装置的设计数据、催化剂及吸附剂的物化性质；操作参数包括装置的操作条件、物料平衡以及流股的流量、压力、组成。
数学模拟模型模块22中，供氢单元模型包括简易供氢模型或详细供氢模型；其中，详细供氢模型包括水蒸汽裂解制氢模型、煤制氢模型，详细供氢模型为考虑氢气生产过程中存在的各个反应器的生产工艺过程模型，用于模拟动力学、热效应及催化剂活性，基于变化的进料组成、流量预测产物收率、组成；简易供氢模型包括重整副产氢模型及其他供氢装置模型，简易供氢模型为仅考虑装置的供氢流量、组成的模型；
用氢单元模型包括加氢精制处理模型和加氢裂化模型；
氢回收单元模型包括变压吸附装置模型、膜分离装置模型、轻烃回收装置、脱硫装置简易模型；
瓦斯系统模型采用均相管网模型，用于描述氢气系统中部分含氢流股排放至瓦斯管线或直送至燃烧炉的状况。
第一确定模块24中氢气系统目标目标函数为年总成本最小、总操作成本最少或氢气公用工程量最小，优选为氢气公用工程量最小。
用氢单元数学规划模型模块25中，在用氢单元数学规划模型进行优化求解时，以新氢流量、循环氢流量为优化变量，以最小化总用氢量为目标函数值，以满足装置最小用氢需求、设备运行的要求为约束条件；用氢单元优化求解采用序列二次规划算法、序贯线性规划算法，得到求解优化结果。
用氢单元数学规划模型模块25中，用氢单元优化求解结果通过物料连接关系传递给下游装置模型。
第一判断模块26中，含氢流股是否可作为用氢单元新氢的标准为：含氢流股氢纯度大于某用氢装置反应器入口氢纯度，CH4体积分数小于CH4设定值，气体杂质含量小于气体杂质设定值， CH4设定值和气体杂质设定值需根据装置工艺要求进行设定。
氢回收单元数学规划模型模块27中，氢回收单元采用基于协同优化理论的建模方法，将整个氢回收系统作为系统级优化研究对象，从而建立整体系统优化模型；将变压吸附装置、膜分离装置作为氢回收系统下的各子系统优化研究对象，从而建立子系统优化模型；系统级优化对象与各子系统通过欲回收的含氢流股相耦合。
（a）整体系统优化模型如下：
整体系统优化模型的目标函数：
式中，
表示氢气回收系统的总回收效益；
表示各氢回收装置的回收效益；
表示模型决策变量；
整体系统优化模型的约束条件：
氢气回收系统与子系统i等式一致性约束为：
……
……
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量，ε表示预设松弛因子。可以理解的是，这里的子系统i表示变压吸附单元或膜分离单元。N表示变压吸附单元或膜分离单元中包含的变压吸附装置或膜分离装置的数量；
（b）子系统优化模型如下：
子系统优化模型的目标函数：
其中，
表示子系统i实际回收处理含氢流股r j 的流量，是子系统i传递给系统级的一个常量；
表示氢气回收系统中预定的子系统i回收处理含氢流股r j 中的流量，是系统级变量。
子系统优化模型的约束条件：
式中： 
、
为子系统i的等式不等式约束。
第二确定模块28中，由用氢单元、氢回收单元优化结果经物料衡算确定氢气公用工程量。
第二判断模块29，第二判断模块，所述达到计算条件：若两次目标函数计算相差满足收敛要求，则结束系统计算；若两次目标函数计算相差未达到收敛要求，且迭代次数仍未达到设计上限，则返回用氢单元数学规划模型模块再进行优化计算；若两次目标函数计算相差未满足收敛要求，但迭代次数达到上限，则停止计算。
需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，可以通过任何合适的方式进行任意组合，其同样落入本发明所公开的范围之内。另外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。
以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
Claims (18)
1.一种氢气资源的优化方法，其特征在于，包括以下步骤：
步骤S1：将氢气系统划分为供氢单元、用氢单元、氢回收单元，确定各单元的设计参数及操作参数；
步骤S2：分别建立供氢单元模型、用氢单元模型、氢回收单元模型、瓦斯系统模型，各模型间按实际物料关系建立连接，进而建立整个氢系统的非线性的数学模拟模型；
步骤S3：对非线性的数学模拟模型进行模拟求解，得到模拟运算结果；
步骤S4：确定整个氢系统的目标函数；
步骤S5：建立用氢单元数学规划模型，然后优化求解；其中，所述用氢单元数学规划模型的优化求解可将步骤S3模拟计算结果作为优化运算的初始可行解，在优化变量的取值范围内，采用优化算法进行求解，得到求解优化结果；
步骤S6：判断氢回收单元各流股是否有作为用氢单元的供氢，如有作为用氢单元的供氢，则返回步骤S5；未有作为用氢单元的供氢，则执行步骤S7；
步骤S7：建立氢回收单元数学规划模型，优化求解；所述氢回收单元优化求解是将步骤S3模拟结果为优化运算初始可行解，步骤S5用氢单元优化结果为进料条件，在优化变量取值范围内，采用优化算法进行求解，得到求解优化结果；
步骤S8：确定供氢单元的氢气公用工程量；
步骤S9：判断整个氢系统的目标函数是否达到计算条件；达到计算条件，则将优化求解结果作为氢气系统模型的运行参数，未达到计算条件，则返回步骤S5。
2.根据权利要求1所述的方法，其特征在于：在所述步骤S1中，所述供氢单元包括水蒸汽裂解制氢、煤制氢、重整副产氢、乙烯装置副产氢、离子膜电解水副产氢、焦化装置、催化裂化装置中的一种或几种装置，所述用氢单元包括加氢精制装置、加氢裂化装置中的一种或几种装置；所述氢回收单元包括变压吸附装置、膜分离装置中的一种或几种装置；
所述设计参数包括装置的设计数据、催化剂及吸附剂的物化性质；所述操作参数包括装置的操作条件，物料平衡，以及流股的流量、压力和组成。
3.根据权利要求1所述的方法，其特征在于：在所述步骤S2中，所述供氢单元模型包括简易供氢模型或详细供氢模型；其中，详细供氢模型包括水蒸汽裂解制氢模型、煤制氢模型，所述详细供氢模型为考虑氢气生产过程中存在的各个反应器的生产工艺过程模型，用于模拟动力学、热效应及催化剂活性，基于变化的进料组成和流量预测产物收率和组成；所述简易供氢模型包括重整副产氢模型及其他供氢装置模型，所述简易供氢模型为仅考虑装置的供氢流量和组成的模型；
所述用氢单元模型包括加氢精制处理模型和加氢裂化模型；
所述氢回收单元模型包括变压吸附装置模型、膜分离装置模型、轻烃回收装置、脱硫装置简易模型；
所述瓦斯系统模型采用均相管网模型，用于描述氢气系统中部分含氢流股排放至瓦斯管线或直送至燃烧炉的状况。
4.根据权利要求1所述的方法，其特征在于：在所述步骤S4中，氢气系统目标函数为年总成本最小、总操作成本最少或氢气公用工程量最小。
5.根据权利要求1所述的方法，其特征在于：在所述步骤S5中，在用氢单元数学规划模型进行优化求解时，以新氢流量、循环氢流量为优化变量，以最小化总用氢量为目标函数值，以满足装置最小用氢需求、设备运行的要求为约束条件；所述用氢单元的优化求解采用优化算法进行求解，得到求解优化结果；所述用氢单元优化求解结果通过物料连接关系传递给下游装置模型。
6.根据权利要求1所述的方法，其特征在于：在所述步骤S6中，含氢流股可作为用氢单元新氢的标准为：含氢流股氢纯度大于用氢装置反应器入口氢纯度，CH4体积分数小于CH4设定值，气体杂质含量小于气体杂质设定值， CH4设定值和气体杂质设定值需根据装置工艺要求进行设定。
7.根据权利要求1所述的方法，其特征在于：在所述步骤S7中，氢回收单元采用基于协同优化理论的建模方法，将整个氢回收系统作为系统级优化研究对象，从而建立整体系统优化模型；将变压吸附装置、膜分离装置作为氢回收系统下的各子系统优化研究对象，从而建立子系统优化模型；所述系统级优化对象与各子系统通过欲回收的含氢流股相耦合。
8.根据权利要求1所述的方法，其特征在于：在所述步骤S8中，由用氢单元优化结果和氢回收单元优化结果经物料衡算确定氢气公用工程量。
9.根据权利要求1所述的方法，其特征在于：在所述步骤S9中，所述达到计算条件：若两次目标函数计算相差满足收敛要求，则结束系统计算；若两次目标函数计算相差未达到收敛要求，且迭代次数仍未达到设计上限，则返回步骤S5再进行优化计算；若两次目标函数计算相差未满足收敛要求，但迭代次数达到上限，则停止计算。
10.一种氢气资源的优化系统，其特征在于，包括以下模块：划分模块、数学模拟模型模块、模拟求解模块、第一确定模块、用氢单元数学规划模型模块、第一判断模块、氢回收单元数学规划模型模块、第二确定模块和第二判断模块；
划分模块，用于将氢气系统划分为供氢单元、用氢单元、氢回收单元，确定各单元的设计参数及操作参数；
数学模拟模型模块，用于分别建立供氢单元模型、用氢单元模型、氢回收单元模型、瓦斯系统模型，各模型间按实际物料关系建立连接，进而建立整个氢系统的非线性的数学模拟模型；
模拟求解模块，用于对非线性的数学模拟模型进行模拟求解，得到模拟运算结果；
第一确定模块，用于确定整个氢系统的目标函数；
用氢单元数学规划模型模块，用于建立用氢单元数学规划模型，然后优化求解；其中，所述用氢单元数学规划模型的优化求解可将模拟求解模块的模拟计算结果作为优化运算的初始可行解，在优化变量的取值范围内，采用优化算法进行求解，得到求解优化结果；
第一判断模块，用于判断氢回收单元各流股是否有作为用氢单元的供氢，如有作为用氢单元的供氢，则返回用氢单元数学规划模型模块；未有作为用氢单元的供氢，则执行氢回收单元数学规划模型模块；
氢回收单元数学规划模型模块，用于建立氢回收单元数学规划模型，优化求解；所述氢回收单元优化求解是将模拟求解模块的模拟结果为优化运算初始可行解，以用氢单元数学规划模型模块的用氢单元优化结果为进料条件，在优化变量取值范围内，采用优化算法进行求解，得到求解优化结果；
第二确定模块，用于确定供氢单元的氢气公用工程量；
第二判断模块，用于判断整个氢系统的目标函数是否达到计算条件；达到计算条件，则将优化求解结果作为氢气系统模型的运行参数，未达到计算条件，则返回用氢单元数学规划模型模块。
11.根据权利要求10所述的系统，其特征在于：在所述划分模块中，所述供氢单元包括焦化装置、催化裂化装置中的一种或几种装置，所述用氢单元包括加氢精制装置、加氢裂化装置中的一种或几种装置；所述氢回收单元包括变压吸附装置、膜分离装置中的一种或几种装置；
所述设计参数包括装置的设计数据、催化剂及吸附剂的物化性质；所述操作参数包括装置的操作条件，物料平衡，以及流股的流量、压力和组成。
12.根据权利要求10所述的系统，其特征在于：所述数学模拟模型模块中，所述供氢单元模型包括简易供氢模型或详细供氢模型；其中，详细供氢模型包括水蒸汽裂解制氢模型、煤制氢模型，所述详细供氢模型为考虑氢气生产过程中存在的各个反应器的生产工艺过程模型，用于模拟动力学、热效应及催化剂活性，基于变化的进料组成和流量预测产物收率和组成；所述简易供氢模型包括重整副产氢模型及其他供氢装置模型，所述简易供氢模型为仅考虑装置的供氢流量和组成的模型；
所述用氢单元模型包括加氢精制处理模型和加氢裂化模型；
所述氢回收单元模型包括变压吸附装置模型、膜分离装置模型、轻烃回收装置、脱硫装置简易模型；
所述瓦斯系统模型采用均相管网模型，用于描述氢气系统中部分含氢流股排放至瓦斯管线或直送至燃烧炉的状况。
13.根据权利要求10所述的系统，其特征在于：在所述第一确定模块中，氢气系统的目标函数为年总成本最小、总操作成本最少或氢气公用工程量最小。
14.根据权利要求10所述的系统，其特征在于：所述用氢单元数学规划模型模块中，在用氢单元数学规划模型进行优化求解时，以新氢流量、循环氢流量为优化变量，以最小化总用氢量为目标函数值，以满足装置最小用氢需求、设备运行的要求为约束条件；所述用氢单元的优化求解采用优化算法进行求解，得到求解优化结果；所述用氢单元优化求解结果通过物料连接关系传递给下游装置模型。
15.根据权利要求10所述的系统，其特征在于：所述第一判断模块中，含氢流股可作为用氢单元新氢的标准为：含氢流股氢纯度大于用氢装置反应器入口氢纯度，CH4体积分数小于CH4设定值，气体杂质含量小于气体杂质设定值， CH4设定值和气体杂质设定值需根据装置工艺要求进行设定。
16.根据权利要求10所述的系统，其特征在于：所述氢回收单元数学规划模型模块中，氢回收单元采用基于协同优化理论的建模方法，将整个氢回收系统作为系统级优化研究对象，从而建立整体系统优化模型；将变压吸附装置、膜分离装置作为氢回收系统下的各子系统优化研究对象，从而建立子系统优化模型；所述系统级优化对象与各子系统通过欲回收的含氢流股相耦合。
17.根据权利要求10所述的系统，其特征在于：所述第二确定模块中，由用氢单元优化结果和氢回收单元优化结果经物料衡算确定氢气公用工程量。
18.根据权利要求10所述的系统，其特征在于：第二判断模块，所述达到计算条件：若两次目标函数计算相差满足收敛要求，则结束系统计算；若两次目标函数计算相差未达到收敛要求，且迭代次数仍未达到设计上限，则返回用氢单元数学规划模型模块再进行优化计算；若两次目标函数计算相差未满足收敛要求，但迭代次数达到上限，则停止计算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710563606.XA CN109255461B (zh) | 2017-07-12 | 2017-07-12 | 一种氢气资源的优化方法及优化系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710563606.XA CN109255461B (zh) | 2017-07-12 | 2017-07-12 | 一种氢气资源的优化方法及优化系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109255461A CN109255461A (zh) | 2019-01-22 |
| CN109255461B true CN109255461B (zh) | 2022-03-08 |
Family
ID=65050562
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710563606.XA Active CN109255461B (zh) | 2017-07-12 | 2017-07-12 | 一种氢气资源的优化方法及优化系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109255461B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110046474B (zh) * | 2019-05-28 | 2023-08-22 | 中国恩菲工程技术有限公司 | 模拟变温变压吸附器的方法及系统 |
| CN112749823A (zh) * | 2019-10-30 | 2021-05-04 | 中国石油化工股份有限公司 | 存储器、基于回收氢最大化的氢网优化方法、装置和设备 |
| CN112749459A (zh) * | 2019-10-30 | 2021-05-04 | 中国石油化工股份有限公司 | 存储器、基于设备改建的氢网用氢优化方法、装置和设备 |
| CN112750503A (zh) * | 2019-10-30 | 2021-05-04 | 中国石油化工股份有限公司 | 存储器、基于氢利用最大化的氢网优化方法、装置和设备 |
| CN114069661B (zh) * | 2021-10-15 | 2023-07-28 | 国网浙江省电力有限公司嘉善县供电公司 | 一种适于波动性输入的电解水制氢系统性能优化方法 |
| CN114757126B (zh) * | 2022-05-06 | 2024-03-22 | 大连理工大学 | 一种基于随机夹点的氢网络重构方法 |
| CN117132066A (zh) * | 2023-08-31 | 2023-11-28 | 珠江水利委员会珠江水利科学研究院 | 一种复杂河网区的水安全协同调控方法及系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104891439A (zh) * | 2015-06-11 | 2015-09-09 | 四川天一科技股份有限公司 | 一种重整气提高氢气回收率的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8082549B2 (en) * | 2006-11-08 | 2011-12-20 | Board Of Regents, The University Of Texas System | System, method and apparatus for allocating resources by constraint selection |
-
2017
- 2017-07-12 CN CN201710563606.XA patent/CN109255461B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104891439A (zh) * | 2015-06-11 | 2015-09-09 | 四川天一科技股份有限公司 | 一种重整气提高氢气回收率的方法 |
Non-Patent Citations (4)
| Title |
|---|
| 数学规划法优化氢网络的应用研究;黄风林等;《化学工程》;20150315(第03期);第70-73页 * |
| 炼厂氢气系统与轻烃回收系统的集成优化;张龙等;《石油与天然气化工》;20151225;第44卷(第6期);第43-48页 * |
| 炼厂氢资源优化利用技术进展;孟凡忠等;《化工进展》;20150915;第66-69页 * |
| 考虑轻烃收益的炼油厂氢气资源优化策略;孟凡忠等;《炼油技术与工程》;20150915(第09期);第61-63页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109255461A (zh) | 2019-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109255461B (zh) | 一种氢气资源的优化方法及优化系统 | |
| CN110020444B (zh) | 一种氢气资源的优化方法及系统 | |
| Marques et al. | State-of-the-art review of targeting and design methodologies for hydrogen network synthesis | |
| US20100152900A1 (en) | Optimizing refinery hydrogen gas supply, distribution and consumption in real time | |
| CN108108572B (zh) | 一种用于蜡油加氢裂化过程建模和优化方法 | |
| Fonseca et al. | Hydrogen distribution network optimization: a refinery case study | |
| CN106485341B (zh) | 一种炼油厂氢气系统的优化方法及其应用 | |
| Jiao et al. | Improved optimization methods for refinery hydrogen network and their applications | |
| Zhao et al. | Integrated model of refining and petrochemical plant for enterprise-wide optimization | |
| CN110187635B (zh) | 用于连续重整装置的实时优化方法和设备 | |
| CN113764046A (zh) | 一种提高催化柴油加氢转化高附加值产品收率的方法 | |
| CN104463343A (zh) | 一种预测催化裂化轻质油产率的方法 | |
| US20220299952A1 (en) | Control system with optimization of neural network predictor | |
| Jagannath et al. | Retrofit design of hydrogen network in refineries: mathematical model and global optimization | |
| da Silva et al. | Integration of hydrogen network design to the production planning in refineries based on multi-scenarios optimization and flexibility analysis | |
| Jia | Refinery hydrogen network optimisation with improved hydroprocessor modelling | |
| CN107977736B (zh) | 基于热值平衡的燃料气系统的优化方法及系统 | |
| CN111863145B (zh) | 一种低温煤焦油全馏分加氢裂化集总动力学模型建模方法 | |
| Zhang et al. | Simultaneous optimization of energy and materials based on heat exchanger network simulation for diesel hydrotreating units | |
| CN111383719B (zh) | 用于炼油过程中氢气与瓦斯的协同优化方法 | |
| Jiao et al. | An optimization method for the refinery hydrogen network and its application | |
| CN208136178U (zh) | 一种提高加氢裂化过程中航煤、汽油收率的系统 | |
| Sadighi et al. | Modeling and optimizing a vacuum gas oil hydrocracking plant using an artificial neural network | |
| Wu et al. | Stepwise optimization of hydrogen network integrated sulfur compound removal kinetics and a fluid catalytic cracker | |
| Alhajri | Integration of hydrogen and CO2 management within refinery planning |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20231010 Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Patentee after: CHINA PETROLEUM & CHEMICAL Corp. Patentee after: Sinopec (Dalian) Petrochemical Research Institute Co.,Ltd. Address before: 100020 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Patentee before: CHINA PETROLEUM & CHEMICAL Corp. Patentee before: DALIAN RESEARCH INSTITUTE OF PETROLEUM AND PETROCHEMICALS, SINOPEC Corp. |