CN101315545B - 一种加氢站高效加氢的三级氢气加注系统 - Google Patents
一种加氢站高效加氢的三级氢气加注系统 Download PDFInfo
- Publication number
- CN101315545B CN101315545B CN2008100635841A CN200810063584A CN101315545B CN 101315545 B CN101315545 B CN 101315545B CN 2008100635841 A CN2008100635841 A CN 2008100635841A CN 200810063584 A CN200810063584 A CN 200810063584A CN 101315545 B CN101315545 B CN 101315545B
- Authority
- CN
- China
- Prior art keywords
- filling
- pressure
- hydrogen
- hydrogen storage
- storage vessel
- 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
- 238000005984 hydrogenation reaction Methods 0.000 title claims description 30
- 239000001257 hydrogen Substances 0.000 claims abstract description 183
- 229910052739 hydrogen Inorganic materials 0.000 claims abstract description 183
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 172
- 238000000034 method Methods 0.000 claims abstract description 55
- 150000002431 hydrogen Chemical class 0.000 claims abstract description 12
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 239000007789 gas Substances 0.000 claims description 141
- 238000003860 storage Methods 0.000 claims description 119
- 238000005429 filling process Methods 0.000 claims description 35
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 239000012159 carrier gas Substances 0.000 claims description 2
- 238000001595 flow curve Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 4
- 239000000446 fuel Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000005273 aeration Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004517 catalytic hydrocracking Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005868 electrolysis reaction Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/32—Hydrogen storage
Landscapes
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
本发明涉及一种加氢优化控制方法及其系统，旨在提供一种加氢站高效加氢的三级加注优化控制方法，即加氢站中储气罐分为高中低三组；为了同时实现高取气率和快速加注的目标，提出了多目标优化算法，采用约束法将多目标问题转化为若干个单目标优化问题；每个单目标规划问题属于动态规划范筹，采用由前向后的顺序搜索方法；根据优化结果设计出具体的程序控制流程，为加氢站的加注系统设计出最佳的加注方式。本发明还提供了一种用于前述加氢站高效加氢的三级加注优化控制方法的氢气加注系统，包括：储氢容器、充气阀、压力传感器、氢气分配器、流量计和单向阀；储氢容器包括低压储氢容器、中压储氢容器和高压储氢容器；还包括一个程序控制器。
Description
技术领域
本发明涉及一种加氢系统，更具体的说，是涉及种加氢站高效加氢的三级氢气加注系统。
背景技术
在燃料电池汽车推向市场的过程中，目前至少有3个主要的技术壁垒需要克服：(1)小型、轻量、廉价和高效的燃料电池发动机的集成；(2)高储存量车载储氢罐的设计；(3)氢气加注设施的开发和加氢站的建设。发展至今，后者更为关键。目前，氢加注系统建设主要有二大发展趋势：(1)加注气态氢的高参数化，即将车载容器贮氢压力增加到最高经济压力70MPa，以提高氢能汽车的续驶里程；(2)液氢加注方案，液氢燃料供应系统具有体积密度高、压力低等优点。但液化储氢存在二大难点：一是液化能耗大；二是液氢的蒸发损失比较严重，存在安全隐患。因此，目前国内外的加氢站研究更为关注的是高压气态加氢，由此引出高取气率序列取气和快速加注控制策略研究。
直接与加氢站控制有关的技术进展，目前主要是对加注系统组成及控制硬件进行的研究，如同济大学在汽车用氢燃料加氢站系统配置的研究中，给出了加氢站中氢气的压缩系统、储存系统和加注机等的控制系统设计方法；北京航天试验技术研究所则给出了加气机中使用的阀门、流量计、加气枪、拉断阀和接头等控制硬件的性能要求。中国专利CN200610088913.9针对加氢站提出了一种快速充装氢气的制氢加氢站系统及其方法，其加注系统通过使用40MPa储氢容器及75MPa超高压储氢容器的二级加注，来提高氢气的利用率和充装速度。而且，采用此专利其已研制成功国内第一台氢加注机，可对压力、温度、流量、流速实时计量，自动收费，同时对车载气瓶的压力、温度、应力变化进行监控，可耐压45-75MPa，安全保护装置齐备。但此方法只包含对各种参数进行监控，并没有加注过程的优化控制。中国专利CN200710098598.2则给出了一种车厢内设有依次连通的紧缩型水电解设备、低压储氢罐组、专用氢气压缩机、高压储氢罐组以及氢气售气机等组成的专用车载制氢加氢站。总之，目前国内尚无关于高压加氢过程中同时实现高取气率序列取气和快速加注的优化控制优化方法的报道。
发明内容
本发明的目的在于克服现有技术中的不足，提供一种加氢站高效加氢的三级加注优化控制方法，本发明更进一步的目的在于提供用于该加氢站高效加氢的三级加注优化控制方法的氢气加注系统。
为了解决以上问题，本发明是通过如下技术方案实现。
本发明提供了一种加氢站高效加氢的三级加注优化控制方法，包括以下步骤：
(1)获取低压、中压和高压储氢容器及车载气瓶中的氢气初始状态参数，并设定：
V1：代表第i次加注前的低压储氢容器的容积；
V2：代表第i次加注前的中压储氢容器的容积；
m、V：分别代表车载气瓶额定充装质量、容积；
m0：加注前车载气瓶内的氢气质量；
xi、yi、zi：分别代表第i次加注时低压、中压和高压储氢容器的取气质量；
m1i、m2i、m3i：分别代表第i次加注前低压、中压和高压储氢容器中的氢气质量；
hi：第i次加注时的取气次序标志，对于三级加注来说共有6种加注次序，hi＝1，表示按低-中-高的次序取气；hi＝2，表示按低-高-中的次序取气；hi＝3，表示按中-低-高的次序取气；hi＝4，表示按中-高-低的次序取气；hi＝5，表示按高-低-中的次序取气；hi＝6，表示按高-中-低的次序取气；
n(i)：第i次加注过程下可充装的车载气瓶数量；
nmax：在一定加注时间约束下，可充装的车载气瓶数量的最大值；
x平衡：加注到储氢容器跟车载气瓶压力平衡时向车载气瓶中充装的质量；
tk：约束时间，表示搜索过程中设定的时间最大值，用这个约束时间作为判断标准，来判断搜索到的加注方式所用的加注时间是否小于该约束时间，只有小于该约束时间的加注方式才可用；
P：储氢容器压力，是计算流量系数时设定的一个压力参数，不具体代表哪个储氢容器，泛指储氢容器的压力，即一个储氢容器压力对应于一个流量系数；
q、qf：分别代表等熵流动状态下的质量流量和实际摩擦流动状态下的质量流量；
Cd：流量系数，Cd＝qf/q；
qs、qfs：分别代表等熵流动临界状态下的质量流量与摩擦流动临界状态下的质量流量；
qa：代表加注开始时的质量流量；
ms、mfs：分别代表加注过程中达到等熵流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量、加注过程中达到摩擦流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量；
t1：为先取气的储氢容器向车载气瓶加注质量为x的氢气所用的时间；
x1充：为某时间下，先取气的储氢容器已经为车载气瓶充入的氢气质量；
tmin：代表整个加注过程只充装一个车载气瓶时，所用的最短时间；
tmax：代表不考虑加注时间约束时，氢气利用率最高的加注过程中最慢一次加注所用时间；
Nmax：代表不考虑加注时间约束时，所能充装的最大车载气瓶数；
i：表示加注第i个车载气瓶的标志；
(2)将区间tmin至tmax之间按照等分的原则取r个值，确定r个约束时间tk，对每一个约束时间tk进行一次加注过程的优化，搜索获取各个约束时间下的氢气利用率最高的加注过程，即可充装的车载气瓶数最多的最优加注过程；
(3)开始时i＝1，此时低压、中压、高压储氢容器的状态为第一次加注前的状态，根据以下氢气状态方程得到初始质量m11、m21、m31：
pV＝mRT(1+αp/T)(α＝1.9155×10-6，R＝4124.3)
对第一个车载气瓶加注时，如果按低——中——高的取气次序取气，则低压储氢容器与车载气瓶达到压力平衡时充入的气体质量x1平衡为：
低压储氢容器的取气质量范围为
此时将取气质量区间分为若干等分，x1的实际取值范围为各个等分点构成的集合；
对于低压储氢容器的每一个取气质量值，中压储氢容器进行加注时，达到压力平衡时充入的氢气质量为：
中压储氢容器的取气质量范围为
此时将取气质量区间分为若干等分，y1的实际取值范围为各个等分点构成的集合；
当低压、中压储氢容器的取气质量均已确定时，高压储氢容器的取气质量为z1＝m-x1-y1；对第1个车载气瓶加注时，每一个x1和y1值对应于加注该车载气瓶时的一种组合；
(4)对每一种取气质量组合，判断六种取气次序是否可能实现：
对第二步取气的储氢容器判断其加注完毕后压力是否高于车载气瓶的压力，如果是，则此种取气序列不能实现；否则，判断最后取气的储氢容器加注完毕后压力是否高于车载气瓶额度充装压力，如果是，则该取气次序可能实现，否则，该取气次序不能实现；
(5)对于可能实现的取气次序，计算各自的加注时间，其方法如下：
将加注软管之前的管路假定为直径逐渐减小的等熵喷管，将加注软管假定为有摩擦的直管，实际管路系统中氢气的流动状态为摩擦流动；
分别在等熵流动状态下和摩擦流动状态下计算储氢容器压力为10MPa、20MPa、30MPa、40MPa、50MPa、60MPa、70MPa时车载气瓶压力从0到储氢容器压力之间变化时的流量q和qf，得到流量系数Cd，然后对各储氢容器压力下的流量系数按最小二乘法原则拟合为四次曲线；实际某压力下的流量系数的计算，举以下例子说明：储氢容器压力在0——10MPa，其流量系数取10MPa下的流量系数。其它压力下的流量系数取法相同；
0-70MPa之间，摩擦流动临界状态下的压力比和车载气瓶与储氢容器压力的比，两者相差很小，可以统一取40MPa时临界状态下的压比；计算某储氢容器压力下等熵流动状态下的临界压力比步骤如下：在0-70MPa之间，每隔5MPa取一个压力值，计算出等熵流动临界状态时的压力比，再作出P-
曲线，然后对曲线按最小二乘法进行二次拟合，最后，通过拟合得到的多项式计算等熵流动状态下的质量流量；
当压力比
小于或等于实际管路系统的临界压比
时，质量流量恒定，按摩擦流动方程计算；当压力比
大于摩擦流动临界状态时的压力比
时，流量按如下公式计算：qf＝Cdq；
对于先取气的储氢容器来说，如果加注开始时
充气结束时那么先求出该储氢容器的初始压力和加注结束后的压力，然后求出初始质量流量qa和加注结束时的质量流量qb，如果加注开始时
加注结束时
那么分两段求取时间t1，求出
时的质量流量qfs及加注结束时的质量流量qb，
如果加注开始时
加注结束时
那么分四段求取充气时间，其中，计算中所用的
及ms为初始压力pi下求出的等熵流动状态下的临界压力比；将后面的时间质量流量曲线按充气质量等分的原则，处理为两段直线，其中，qb为等分点x1充＝ms+0.5(x-ms)时的质量流量，qc为加注结束时的质量流量；如果加注开始时
加注结束时
qa及qb分别为加注开始和加注结束时的质量流量；如果加注开始时
加注结束时
那么分三段计算充气时间，
qb为等分点的质量流量，qc为加注结束时的质量流量；如果，加注开始时
那么
qb为加注结束时的质量流量；
后取气的两个储氢容器加注时间的计算方法与先取气的储氢容器相同，对此车载气瓶总的加注时间为分别从三个储氢容器中取气的加注时间之和；
(6)选择加注时间最短的那种取气次序：
判断最短的加注时间是否小于或等于约束时间，即是否满足ti≤tk；如果满足，则令加注标志n(i)＝i，根据氢气的状态方程获取各储氢容器加注后的状态参数作为加注第二个车载气瓶时的初始状态，计算出各组储氢容器取气质量的取值范围，进行第二次加注的运算；如果不满足，则继续第一次加注过程的循环搜索，搜索此加注过程中其他的取气质量组合，作同样的判断；
(7)根据第一次加注完毕获得出的加注第二个车载气瓶时的初始状态以及各储氢容器取气质量的取值范围，进行第二次加注的搜索，其方法与第一次加注的搜索相同；
(8)由前面的加注过程逐步向后搜索直到所有的循环搜索全部完成，对于第i个车载气瓶加注过程中的某种组合如果能够充满车载气瓶并满足时间上的要求，那么就令加注标志n(i)＝i，n(i)的初始值为0；并记录下该次及该次以前的加注过程中各组储氢容器的取气质量及取气次序；最后求出数组n(i)，中的最大值——即nmax，i＝1：Nmax，并由程序控制器根据nmax对应的为各车载气瓶加注时各储氢容器的取气质量及取气次序向阀门控制装置输出控制信号。
本发明还提供了一种用于前述加氢站高效加氢的三级加注优化控制方法的氢气加注系统，包括：储氢容器、充气阀、压力传感器、氢气分配器、流量计和单向阀，充气阀、氢气分配器、单向阀依次连接至车载储气瓶；所述储氢容器是三级压力的储罐，包括低压储氢容器、中压储氢容器和高压储氢容器，分别通过各自的控制阀门cv1、cv2、cv3连接至充气阀；还包括一个程序控制器，压力传感器连接氢气分配器和程序控制器，流量计的测量点位于氢气分配器和单向阀之间，并通过信号线连接至程序控制器，程序控制器通过信号线与低压储氢容器、中压储氢容器和高压储氢容器各自的控制阀门cv1、cv2、cv3连接。
与现有技术相比，本发明的有益效果是：
针对氢加注系统的特点，提出了过程优化控制方法，这是在以前的加注系统设计及研究中所没有的；利用本控制方法可对氢气的利用率和加注时间同时进行优化，能确保高的氢气利用率的同时，较大地缩短加注时间；可以根据不同的加注参数及要求，产生不同的优化加注方式；使加注过程简单、直接并易于操作。
附图说明
图1是加注系统优化控制方法流程图；
图2是加氢站的加注系统流程示意图。
图中附图标记：1储气罐、2充气阀、3压力传感器、4氢气分配器、5流量计、6单向阀、7车载储气瓶、8程序控制器、低压储氢容器控制阀门cv1、中压储氢容器控制阀门cv2、高压储氢容器控制阀门cv3。
具体实施方式
结合附图，下面通过具体实施例对本发明进行详细说明。
以下根据加氢站具体的加注参数，用本控制方法实现的实例，并对优化后的具体加注过程进行了详细说明。本发明中的程序控制器8可选用西门子公司生产的PLC型号为S7-400的产品。
实例1
对加氢站进行三级加注优化计算，其中低压储氢容器容积V1＝0.6m3，压力P1＝42MPa；中压储氢容器容积V2＝0.4m3，压力P2＝45MPa；高压储氢容器容积V3＝0.2m3，压力P3＝70MPa；设定车载气瓶容积V＝0.15m3，充满气后的压力P＝35MPa，每次充气前的初始压力P0＝2MPa。
根据发明内容中的控制方法，进行编程求解，获得控制方法的非劣解。而后，根据过程优化的目的和实际要求从其中选出最优解，编制成具体的实际操作控制步骤。
其中控制方法的非劣解如表1所示。
其中，t为本加注方式下，充满一辆车所用的最长的时间。
由上表可以根据不同的加注要求，来选择不同的加注方式。为了同时达到高的氢气利用率和快速充气的目的，现以充气时间为242s，充气罐数为6罐的方式为最佳。
其中，氢气利用率可以通过以下公式计算：
式中m为车载气瓶规定的充气质量，n为实际可充装的气瓶数量，M为各储氢容器的氢气总质量。氢气质量的计算方法：
P为压力，单位Pa；V为容积，单位m3；R＝4124.3.3，T＝293.15K，α＝1.9155×10-6。
将以上结果对应的加注方式编入程序控制器，通过程控器来实现对加注过程中阀门的切换和优化控制。
根据本控制方法得到的优化加注方式的具体操作步骤如下：
a)开始第一辆车的充气，由程序控制器控制开启阀门2，并打开cv1先利用低压组的储气罐进行充气，当流量计5的流量达到3.0750kg时，低压组储气罐的阀门cv1关闭，中压组储气罐cv2的阀门打开，此时当流量计5测试的流量达到0.2141kg时，切断中压组储气罐cv2的阀门，此时车载储气瓶充满。关闭阀门2完成第一辆车的充气过程；
b)开始第二辆车的充气，由程序控制器控制开启阀门2，并打开cv1先利用低压组的储气罐进行充气，当流量计5的流量达到2.4600kg时，低压组储气罐的阀门cv1关闭，中压组储气罐cv2的阀门打开，此时当流量计5测试的流量达到0.8291kg时，切断中压组储气罐cv2的阀门，此时车载储气瓶充满。关闭阀门2完成第二辆车的充气过程；
c)接下来的车辆按照此方法，根据程序控制器设定好的程序来进行充气。
实例2
加氢站低压组固定容器容积V1＝0.6m3，压力P1＝42MPa，中压组固定容器容积V2＝0.4m3，压力P2＝70MPa，高压组固定容器容积V3＝0.2m3，压力P3＝70MPa，设定车载气瓶容积V＝0.15m3，充满气后的压力P＝35MPa，每次充气前的初始压力P0＝2MPa。
根据本控制方法，可以得到以上几种加注方式。根据具体的加注要求可确定最佳方式，如当充气为8个罐的时候，其氢气利用率为
而当充气为7个罐的时，氢气利用率为
为了同时具有较短的充气时间和高取气率，优先选用充气时间为181s，利用率为57.92％的加注方式，其与平衡加注方式相比可节省充气时间近50％。另外，在利用率均为57.92％的几种加注方式中，可以明显地观察到，所用时间最短，即181s，那种加注方式中，存在加注罐顺序的改变。如在此加注方式中，在充装第6罐的时候不是采用低-中-高的加注顺序，而是直接用中-高充气，到第7罐又恢复到低-中-高的加注顺序。可见，由于存在这种加注次序的自动切换，可以在相同利用率的情况下使其加注时间明显缩短。由于现在给出的实例中，只是系统一轮次加注过程，实际中，通过压缩机补气，采用的是多轮次加注。因而，当加注车辆更多时，采用本方法来改变加注次序并同时实现高利用率和缩短加注时间的效果将会更加突出。
最后，需要注意的是，以上列举的仅是本发明的具体实施例。显然，本发明不限于以上实施例，还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形，均应认为是本发明的保护范围。
Claims (1)
1.一种加氢站高效加氢的三级氢气加注系统，包括：储氢容器、充气阀、压力传感器、氢气分配器、流量计和单向阀，充气阀、氢气分配器、单向阀依次连接至车载储气瓶；其特征在于，所述储氢容器是三级压力的储罐，包括低压储氢容器、中压储氢容器和高压储氢容器，分别通过各自的控制阀门(cv1、cv2、cv3)连接至充气阀；还包括一个程序控制器，压力传感器连接氢气分配器和程序控制器，流量计的测量点位于氢气分配器和单向阀之间，并通过信号线连接至程序控制器，程序控制器通过信号线与低压储氢容器、中压储氢容器和高压储氢容器各自的控制阀门(cv1、cv2、cv3)连接；
所述三级氢气加注系统，其优化的氢气加注控制通过以下方法实现：
(1)获取低压、中压和高压储氢容器及车载气瓶中的氢气初始状态参数，并设定：
V1：代表第i次加注前的低压储氢容器的容积；
V2：代表第i次加注前的中压储氢容器的容积；
m、V：分别代表车载气瓶额定充装质量、容积；
m0：加注前车载气瓶内的氢气质量；
xi、yi、zi：分别代表第i次加注时低压、中压和高压储氢容器的取气质量；
m1i、m2i、m3i：分别代表第i次加注前低压、中压和高压储氢容器中的氢气质量；
hi：第i次加注时的取气次序标志，对于三级加注来说共有6种取气次序，hi＝1，表示按低-中-高的次序取气；hi＝2，表示按低-高-中的次序取气；hi＝3，表示按中-低-高的次序取气；hi＝4，表示按中-高-低的次序取气；hi＝5，表示按高-低-中的次序取气；hi＝6，表示按高-中-低的次序取气；
n(i)：第i次加注过程下可充装的车载气瓶数量；
nmax：在一定加注时间约束下，可充装的车载气瓶数量的最大值；
x平衡：加注到储氢容器跟车载气瓶压力平衡时向车载气瓶中充装的质量；
tk：约束时间，表示搜索过程中设定的时间最大值，用这个约束时间作为判断标准，来判断搜索到的加注方式所用的加注时间是否小于该约束时间，只有小于该约束时间的加注方式才可用；
P：储氢容器压力，是计算流量系数时设定的一个压力参数，不具体代表哪个储氢容器，泛指储氢容器的压力，即一个储氢容器压力对应于一个流量系数；
q、qf：分别代表等熵流动状态下的质量流量和实际摩擦流动状态下的质量流量；
Cd：流量系数，Cd＝qf/q；
qs、qfs：分别代表等熵流动临界状态下的质量流量与摩擦流动临界状态下的质量流量；
qa：代表加注开始时的质量流量：
ms、mfs：分别代表加注过程中达到等熵流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量、加注过程中达到摩擦流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量；
t1：为先取气的储氢容器向车载气瓶加注质量为x的氢气所用的时间；
x1充：为某时间下，先取气的储氢容器已经为车载气瓶充入的氢气质量；
tmin：代表整个加注过程只充装一个车载气瓶时，所用的最短时间；
tmax：代表不考虑加注时间约束时，氢气利用率最高的加注过程中最慢一次加注所用时间；
Nmax：代表不考虑加注时间约束时，所能充装的最大车载气瓶数；
i：表示加注第i个车载气瓶的标志；
(2)将区间tmin至tmax之间按照等分的原则取r个值，确定r个约束时间tk，对每一个约束时间tk进行一次加注过程的优化，搜索获取各个约束时间下的氢气利用率最高的加注过程，即可充装的车载气瓶数最多的最优加注过程；
(3)开始时i＝1，此时低压、中压、高压储氢容器的状态为第一次加注前的状态，根据以下氢气状态方程得到初始质量m11、m21、m31：
pV＝mRT(1+αp/T)，其中α＝1.9155×10-6，R＝4124.3，T＝293.15K
对第一个车载气瓶加注时，如果按低——中——高的取气次序取气，则低压储氢容器与车载气瓶达到压力平衡时充入的气体质量x1平衡为：
低压储氢容器的取气质量范围为此时将取气质量区间分为若干等分，x1的实际取值范围为各个等分点构成的集合；
对于低压储氢容器的每一个取气质量值，中压储氢容器进行加注时，达到压力平衡时充入的氢气质量为：中压储氢容器的取气质量范围为
此时将取气质量区间分为若干等分，y1的实际取值范围为各个等分点构成的集合；
当低压、中压储氢容器的取气质量均已确定时，高压储氢容器的取气质量为z1＝m-x1-y1；对第1个车载气瓶加注时，每一个x1和y1值对应于加注该车载气瓶时的一种组合；
(4)对每一种取气质量组合，判断六种取气次序是否可能实现：
对第二步取气的储氢容器判断其加注完毕后压力是否高于车载气瓶的压力，如果是，则此种取气序列不能实现；否则，判断最后取气的储氢容器加注完毕后压力是否高于车载气瓶额度充装压力，如果是，则该取气次序可能实现，否则，该取气次序不能实现；
(5)对于可能实现的取气次序，计算各自的加注时间，其方法如下：
将加注软管之前的管路假定为直径逐渐减小的等熵喷管，将加注软管假定为有摩擦的直管，实际管路系统中氢气的流动状态为摩擦流动；
分别在等熵流动状态下和摩擦流动状态下计算储氢容器压力为10MPa、20MPa、30MPa、40MPa、50MPa、60MPa、70MPa时车载气瓶压力从0到储氢容器压力之间变化时的流量q和qf，得到流量系数Cd，然后对各储氢容器压力下的流量系数按最小二乘法原则拟合为四次曲线；
0-70MPa之间，摩擦流动临界状态下的压力比和车载气瓶与储氢容器压力的比，两者相差很小，可以统一取40MPa时临界状态下的压比；计算某储氢容器压力下等熵流动状态下的临界压力比步骤如下：在0-70MPa之间，每隔5MPa取一个压力值，计算出等熵流动临界状态时的压力比，再作出P-
曲线，然后对曲线按最小二乘法进行二次拟合，最后，通过拟合得到的多项式计算等熵流动状态下的质量流量；
当压力比
小于或等于实际管路系统的临界压比
时，质量流量恒定，按摩擦流动方程计算；当压力比
大于摩擦流动临界状态时的压力比
时，流量按如下公式计算：qf＝Cdq；
对于先取气的储氢容器来说，如果加注开始时充气结束时
那么先求出该储氢容器的初始压力和加注结束后的压力，然后求出初始质量流量qa和加注结束时的质量流量qb，如果加注开始时
加注结束时
那么分两段求取时间t1，求出时的质量流量qfs及加注结束时的质量流量qb，
如果加注开始时
加注结束时
那么分四段求取充气时间，其中，计算中所用的
及ms为初始压力pi下求出的等熵流动状态下的临界压力比；将后面的时间质量流量曲线按充气质量等分的原则，处理为两段直线，
其中，qb为等分点x1充＝ms+0.5(x-ms)时的质量流量，qc为加注结束时的质量流量；如果加注开始时
加注结束时
qa及qb分别为加注开始和加注结束时的质量流量；如果加注开始时加注结束时
那么分三段计算充气时间，
qb为等分点的质量流量，qc为加注结束时的质量流量；如果，加注开始时
那么
qb为加注结束时的质量流量；
后取气的两个储氢容器加注时间的计算方法与先取气的储氢容器相同，对此车载气瓶总的加注时间为分别从三个储氢容器中取气的加注时间之和；
(6)选择加注时间最短的那种取气次序：
判断最短的加注时间是否小于或等于约束时间，即是否满足ti≤tk；如果满足，则令加注标志n(i)＝i，根据氢气的状态方程获取各储氢容器加注后的状态参数作为加注第二个车载气瓶时的初始状态，计算出各组储氢容器取气质量的取值范围，进行第二次加注的运算；如果不满足，则继续第一次加注过程的循环搜索，搜索此加注过程中其他的取气质量组合，作同样的判断；
(7)根据第一次加注完毕获得出的加注第二个车载气瓶时的初始状态以及各储氢容器取气质量的取值范围，进行第二次加注的搜索，其方法与第一次加注的搜索相同；
(8)由前面的加注过程逐步向后搜索直到所有的循环搜索全部完成，对于第i个车载气瓶加注过程中的某种组合如果能够充满车载气瓶并满足时间上的要求，那么就令加注标志n(i)＝i，n(i)的初始值为0；并记录下该次及该次以前的加注过程中各组储氢容器的取气质量及取气次序；最后求出数组n(i)中的最大值——即nmax，i＝1∶Nmax，并由程序控制器根据nmax对应的为各车载气瓶加注时各储氢容器的取气质量及取气次序向阀门控制装置输出控制信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2008100635841A CN101315545B (zh) | 2008-06-27 | 2008-06-27 | 一种加氢站高效加氢的三级氢气加注系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2008100635841A CN101315545B (zh) | 2008-06-27 | 2008-06-27 | 一种加氢站高效加氢的三级氢气加注系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101315545A CN101315545A (zh) | 2008-12-03 |
| CN101315545B true CN101315545B (zh) | 2010-06-09 |
Family
ID=40106576
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2008100635841A Active CN101315545B (zh) | 2008-06-27 | 2008-06-27 | 一种加氢站高效加氢的三级氢气加注系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101315545B (zh) |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101832467B (zh) * | 2010-06-03 | 2011-09-14 | 浙江大学 | 基于温升与能耗控制的加氢站氢气优化加注方法及其系统 |
| CN102637247A (zh) * | 2012-02-23 | 2012-08-15 | 浙江工商大学 | 一种集成概化与消隐的高维数据匿名处理方法 |
| US9279541B2 (en) * | 2013-04-22 | 2016-03-08 | Air Products And Chemicals, Inc. | Method and system for temperature-controlled gas dispensing |
| JP6001600B2 (ja) * | 2014-06-26 | 2016-10-05 | 株式会社日本自動車部品総合研究所 | ガス供給部を制御するシステムおよびガス充填方法 |
| CN106227907B (zh) * | 2016-05-30 | 2019-06-07 | 重庆大学 | 一种基于演化聚类学习的液氢加注系统安全性评估方法 |
| CN108613012A (zh) * | 2016-12-09 | 2018-10-02 | 中国石油化工股份有限公司 | 一种加氢站气动式氢气压缩储气系统 |
| CN107139943A (zh) * | 2017-05-12 | 2017-09-08 | 中车唐山机车车辆有限公司 | 有轨电车交通系统 |
| CN109945067A (zh) * | 2019-03-21 | 2019-06-28 | 浙江昊凡科技有限公司 | 一种加氢站氢能源快速充装的控制方法及系统 |
| CN112539337B (zh) * | 2019-09-20 | 2022-08-05 | 中国石油化工股份有限公司 | 高压管束车供氢加氢站设备选型方法 |
| CN112539339B (zh) * | 2019-09-20 | 2022-08-23 | 中国石油化工股份有限公司 | 天然气重整供氢加氢站设备选型方法 |
| CN111365607B (zh) * | 2020-03-20 | 2021-11-09 | 重庆大学 | 基于智能预测与控制的车载氢气瓶三级加注方法 |
| CN112902014B (zh) * | 2021-04-07 | 2022-11-18 | 阳光新能源开发股份有限公司 | 一种氢能站及其控制系统和控制方法 |
| CN113357535A (zh) * | 2021-06-10 | 2021-09-07 | 苏州工业园区奥菲特气体设备有限公司 | 分布式加压系统 |
| CN113978316B (zh) * | 2021-10-21 | 2024-04-16 | 潍柴动力股份有限公司 | 一种可续航里程的计算方法、装置及存储介质 |
| CN114484266A (zh) * | 2022-01-27 | 2022-05-13 | 氢华能源技术(武汉)有限公司 | 移动式加氢站 |
| CN114609894B (zh) * | 2022-05-11 | 2022-08-02 | 浙江浙能航天氢能技术有限公司 | 一种基于总线控制的可移动加氢站的控制系统及控制方法 |
| CN115164088B (zh) * | 2022-07-13 | 2024-03-01 | 沈阳工业大学 | 储氢系统及其快速加氢过程的压力和温度优化控制方法 |
| CN116428513B (zh) * | 2023-03-13 | 2024-02-13 | 航天氢能科技有限公司 | 一种基于模型预测控制的加氢站氢气优化加注方法及系统 |
| CN118520895A (zh) * | 2024-07-25 | 2024-08-20 | 山东永安特种装备有限公司 | 基于金属电子标签的工业气瓶信息管理方法、装置及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3754125A (en) * | 1971-01-06 | 1973-08-21 | Du Pont | Vent gas control system |
| US5000924A (en) * | 1987-06-02 | 1991-03-19 | Elsagainternational B.V. | Autoacceleration control for exothermic reactors |
| CN1426985A (zh) * | 2001-12-18 | 2003-07-02 | 石油大学(北京) | 一种乙炔加氢反应器出口乙炔含量的控制方法 |
| CN101033821A (zh) * | 2006-03-10 | 2007-09-12 | 大阳日酸株式会社 | 氢气填充方法及装置 |
| CN101180495A (zh) * | 2005-04-25 | 2008-05-14 | 乔治洛德方法研究和开发液化空气有限公司 | 填充用于气囊系统的气体容器的两步方法及用于两步填充方法的气体填充装置 |
-
2008
- 2008-06-27 CN CN2008100635841A patent/CN101315545B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3754125A (en) * | 1971-01-06 | 1973-08-21 | Du Pont | Vent gas control system |
| US5000924A (en) * | 1987-06-02 | 1991-03-19 | Elsagainternational B.V. | Autoacceleration control for exothermic reactors |
| CN1426985A (zh) * | 2001-12-18 | 2003-07-02 | 石油大学(北京) | 一种乙炔加氢反应器出口乙炔含量的控制方法 |
| CN101180495A (zh) * | 2005-04-25 | 2008-05-14 | 乔治洛德方法研究和开发液化空气有限公司 | 填充用于气囊系统的气体容器的两步方法及用于两步填充方法的气体填充装置 |
| CN101033821A (zh) * | 2006-03-10 | 2007-09-12 | 大阳日酸株式会社 | 氢气填充方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 郑津洋等.高压氢气储运设备及其风险评价.太阳能学报第27卷 第11期.2006,第27卷(第11期),1168-1174. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101315545A (zh) | 2008-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101315545B (zh) | 一种加氢站高效加氢的三级氢气加注系统 | |
| CN201218417Y (zh) | 一种加氢站高效加氢的三级加注优化控制系统 | |
| CN108087717B (zh) | 一种液氢及高压氢联合加氢系统及加注方法 | |
| CN106801787A (zh) | 一种采用预压式气液置换的气瓶充气系统 | |
| CN102155615B (zh) | 具有bog回收功能的lng无泵加气方法的设备 | |
| CN109708000B (zh) | 一种l-ch2型加氢站热管理系统 | |
| Tietze et al. | Bulk storage vessels for compressed and liquid hydrogen | |
| CN201964142U (zh) | 具有bog回收功能的lng无泵加气设备 | |
| CN108644604A (zh) | 低温杜瓦容器以及低温高压储氢系统 | |
| CN101905647A (zh) | 一种机动车辆燃料系统的使用方法 | |
| CN207246782U (zh) | 制氢加氢站 | |
| Li et al. | Investigation on the changes of pressure and temperature in high pressure filling of hydrogen storage tank | |
| CN111692520B (zh) | 一种供应气体的方法 | |
| CN103344292B (zh) | 火箭贮箱推进剂量实时动态测量方法 | |
| CN115388324A (zh) | 一种气、液氢联供加氢站 | |
| CN201946692U (zh) | 车用氢燃料的冷却加注装置 | |
| CN111692521A (zh) | 一种高效灵活的多模式加氢站 | |
| CN216524275U (zh) | 液氢泵驱动的质量法液氢流量标准装置 | |
| CN203868681U (zh) | 一种液化天然气闪蒸气回收装置 | |
| CN216307425U (zh) | 一种现场制氢与外购氢互补的加氢站系统 | |
| CN110715164A (zh) | 一种天然气或氢气加注系统及加注方法 | |
| US11391415B1 (en) | Method for minimizing power demand for hydrogen refueling station | |
| CN213299608U (zh) | 一种多模式加氢站 | |
| Yuan et al. | Effects of structure parameter and material property on thermal performance of on-board hydrogen storage tanks during fast refueling | |
| CN103322410B (zh) | 无动力撬装液化天然气加注站 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |