专利名称：：一种加氢站高效加氢的三级加注优化控制方法及其系统的制作方法
技术领域：
：本发明涉及一种加氢优化控制方法及其系统，更具体的说，是涉及加氢站高效加氢的三级加注优化控制方法及其系统。
背景技术：
：在燃料电池汽车推向市场的过程中，目前至少有3个主要的技术壁垒需要克服(1)小型、轻量、廉价和高效的燃料电池发动机的集成；(2)高储存量车载储氢罐的设计；(3)氢气加注设施的开发和加氢站的建设。发展至今，后者更为关键。目前，氢加注系统建设主要有二大发展趋势(1)加注气态氢的高参数化，即将车载容器贮氢压力增加到最高经济压力70MP"，以提高氢能汽车的续驶里程；(2)液氢加注方案，液氢燃料供应系统具有体积密度高、压力低等优点。但液化储氢存在二大难点一是液化能耗大；二是液氢的蒸发损失比较严重,存在安全隐患。因此，目前国内外的加氢站研究更为关注的是高压气态加氢，由此引出高取气率序列取气和快速加注控制策略研究。直接与加氢站控制有关的技术进展，目前主要是对加注系统组成及控制硬件进行的研究，如同济大学在汽车用氢燃料加氢站系统配置的研究中，给出了加氢站中氢气的压縮系统、储存系统和加注机等的控制系统设计方法；北京航天试验技术研究所则给出了加气机中使用的阀门、流量计、加气枪、拉断阀和接头等控制硬件的性能要求。中国专利CN200610088913.9针对加氢站提出了一种快速充装氢气的制氢加氢站系统及其方法，其加注系统通过使用40M/^储氢容器及75M尸a超高压储氢容器的二级加注，来提高氢气的利用率和充装速度。而且，采用此专利其已研制成功国内第一台氢加注机，可对压力、温度、流量、流速实时计量，自动收费，同时对车载气瓶的压力、温度、应力变化进行监控，可耐压45-75MPa，安全保护装置齐备。但此方法只包含对各种参数进行监控，并没有加注过程的优化控制。中国专利CN200710098598.2则给出了一种车厢内设有依次连通的紧縮型水电解设备、低压储氢罐组、专用氢气压縮机、高压储氢罐组以及氢气售气机等组成的专用车载制氢加氢站。总之，目前国内尚无关于高压加氢过程中同时实现高取气率序列取气和快速加注的优化控制优化方法的报道。
发明内容本发明的目的在于克服现有技术中的不足，提供一种加氢站高效加氢的三级加注优化控制方法，本发明更进一步的目的在于提供用于该加氢站高效加氢的三级加注优化控制方法的氢气加注系统。为了解决以上问题，本发明是通过如下技术方案实现。本发明提供了一种加氢站高效加氢的三级加注优化控制方法，包括以下步骤(1)获取低压、中压和高压储氢容器及车载气瓶中的氢气初始状态参数，并设定F1:代表第i次加注前的低压储氢容器的容积；^:代表第i次加注前的中压储氢容器的容积；w、K:分别代表车载气瓶额定充装质量、容积；m。加注前车载气瓶内的氢气质量；X。^、分别代表第i次加注时低压、中压和高压储氢容器的取气质量；W,、W2,、W3i:分别代表第i次加注前低压、中压和高压储氢容器中的氢气质量；第Z'次加注时的取气次序标志，对于三级加注来说共有6种加注次序，/7,=1，表示按低一中一高的次序取气；A,=2，表示按低一高一中的次序取气；/^=3,表示按中一低一高的次序取气；/^=4，表示按中一高一低的次序取气；&=5，表示按高一低一中的次序取气；/z,=6，表示按高一中一低的次序取气；:第i次加注过程下可充装的车载气瓶数量；"max:在一定加注时间约束下，可充装的车载气瓶数量的最大值；:加注到储氢容器跟车载气瓶压力平衡时向车载气瓶中充装的质量；":约束时间，表示搜索过程中设定的时间最大值，用这个约束时间作为判断标准，来判断搜索到的加注方式所用的加注时间是否小于该约束时间，只有小于该约束时间的加注方式才可用；尸储氢容器压力，是计算流量系数时设定的一个压力参数，不具体代表哪个储氢容器，泛指储氢容器的压力，即一个储氢容器压力对应于一个流量系数；g、分别代表等熵流动状态下的质量流量和实际摩擦流动状态下的质量流量；Crf:流量系数，Q^^；&、？/,:分别代表等熵流动临界状态下的质量流量与摩擦流动临界状态下的质量流量；代表加注开始时的质量流量；9:代表加注过程中车载气瓶与先取气的储氢容器的压力比；^、9、分别代表先取气的储氢容器加注过程中达到等熵流动临界状态时及摩擦流动临界状态时车载气瓶与储氢容器的压力比；m,、分别代表加注过程中达到等熵流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量、加注过程中达到摩擦流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量；f,:为先取气的储氢容器向车载气瓶加注质量为x的氢气所用的时间；为某时间下，先取气的储氢容器已经为车载气瓶充入的氢气质量；tmin:代表整个加注过程只充装一个车载气瓶时，所用的最短时间；tmax:代表不考虑加注时间约束时，氢气利用率最高的加注过程中最慢一次加注所用时间；Nmax:代表不考虑加注时间约束时，所能充装的最大车载气瓶数；i:表示加注第i个车载气瓶的标志；(2)将区间tmin至tn^之间按照等分的原则取r个值，确定r个约束时间^，对每一水约束时间^进行一次加注过程的优化，搜索获取各个约束时间下的氢气利用率最高的加注过程，即可充装的车载气瓶数最多的最优加注过程；(3)开始时/=1，此时低压、中压、高压储氢容器的状态为第一次加注前的状态，根据以下氢气状态方程得到初始质量W,、m21、m31:/7F=wir(l+"j5/r)(a=1.9155xl0—6，i=4124.3)对第一个车载气瓶加注时，如果按低——中——高的取气次序取气，则低压储氢容器与车载气瓶达到压力平衡时充入的气体质量Xl￥fft为、平衡='低压储氢容器的取气质量范围为A"O,min(x,平衡，此时将取气质量区间分为若干等分，Xl的实际取值范围为各个等分点构成的集合；对于低压储氢容器的每一个取气质量值，中压储氢容器进行加注时，达到压力平衡时充入的氢气质量为^=^^^，中压储氢容器的取气质量范围为Ae[O,min"平衡，w-xjj;此时将取气质量区间分为若干等分，yi的实际取值范围为各个等分点构成的集合；当低压、中压储氢容器的取气质量均巳确定时，高压储氢容器的取气质量为^二m-x,-力；对第l个车载气瓶加注时，每一个xjQyi值对应于加注该车载气瓶时的一种组合；(4)对每一种取气质量组合，判断六种取气次序是否可能实现对第二步取气的储氢容器判断其加注完毕后压力是否高于车载气瓶的压力，如果是，则此种取气序列不能实现；否则，判断最后取气的储氢容器加注完毕后压力是否高于车载气瓶额度充装压力，如果是，则该取气次序可能实现，否则，该取气次序不能实现；(5)对于可能实现的取气次序，计算各自的加注时间，其方法如下将加注软管之前的管路假定为直径逐渐减小的等熵喷管，将加注软管假定为有摩擦的直管，实际管路系统中氢气的流动状态为摩擦流动；分别在等熵流动状态下和摩擦流动状态下计算储氢容器压力为lOMPa、20MPa、30MPa、40MPa、50MPa、60MPa、70MPa时车载气瓶压力从0到储氢容器压力之间变化时的流量《和^，得到流量系数(^，然后对各储氢容器压力下的流量系数按最小二乘法原则拟合为四次曲线；实际某压力下的流量系数的计算，举以下例子说明储氢容器压力在O——lOMPa，其流量系数取lOMPa下的流量系数。其它压力下的流量系数取法相同；0-70MPa之间，摩擦流动临界状态下的压力比和车载气瓶与储氢容器压力的比，两者相差很小，可以统一取40MPa时临界状态下的压比；计算某储氢容器压力下等熵流动状态下的临界压力比步骤如下在0-70MPa之间，每隔5MPa取一个压力值，计算出等熵流动临界状态时的压力比，再作出尸一(Ps曲线，然后对曲线按最小二乘法进行二次拟合，最后，通过拟合得到的多项式计算等熵流动状态下的质量流量；当压力比9小于或等于实际管路系统的临界压比cPfc时，质量流量恒定，按摩擦流动方程计算；当压力比cp大于摩擦流动临界状态时的压力比9j寸，流量按如下公式计算&=C^;对于先取气的储氢容器来说，如果加注开始时9〈cpfs，充气结束时9<cpfs，那么先求出该储氢容器的初始压力和加注结束后的压力，然后求出初始质量流量^和加注结束时的质量流量&，t广^^;如果加注开始时(p〈9f,，加注结束时(^〈cp〈A，那么分两段求取时间。，qa+qb求出9=^时的质量流量9&及加注结束时的质量流量&，、=^_+^^;如果加注开始qa+qftqfs+qb时(p〈(Pfs，加注结束时(p^Ps,那么分四段求取充气时间，其中，计算中所用的(Ps及W,为初始压力A下求出的等熵流动状态下的临界压力比；将后面的时间质量流量曲线按充气质量等分的原则，处理为两段直线，tl=~^+2(m,m&)+^1+^1，其中，^为等分点qa+qfsqfs+qsqs+qbqb+qcx^=m,+0.5(x-时的质量流量，《e为加注结束时的质量流量；如果加注开始时(pfs〈(p<(ps，加注结束时q^〈9〈化，、=-，&及qb分别为加注开始和加注结束时的质qa+qb量流量；如果加注开始时9fs〈cp〈9s，加注结束时(P〉化，那么分三段计算充气时间，t严^^+^^+^^，^为等分点的质量流量，^为加注结束时的质量流量；如果，qa+qsqs+qbqb+qc加注开始时cp〉(Ps，那么t,"^，^为加注结束时的质量流量；qa+qb后取气的两个储氢容器加注时间的计算方法与先取气的储氢容器相同，对此车载气瓶总的加注时间为分别从三个储氢容器中取气的加注时间之和；(6)选择加注时间最短的那种取气次序判断最短的加注时间是否小于或等于约束时间，即是否满足《^^;如果满足，则令加注标志w(/)=/，根据氢气的状态方程获取各储氢容器加注后的状态参数作为加注第二个车载气瓶时的初始状态，计算出各组储氢容器取气质量的取值范围，进行第二次加注的运算；如果不满足，则继续第一次加注过程的循环搜索，搜索此加注过程中其他的取气质量组合，作同样的判断；(7)根据第一次加注完毕获得出的加注第二个车载气瓶时的初始状态以及各储氢容器取气质量的取值范围，进行第二次加注的搜索，其方法与第一次加注的搜索相同；(8)由前面的加注过程逐步向后搜索直到所有的循环搜索全部完成，对于第i个车载气瓶加注过程中的某种组合如果能够充满车载气瓶并满足时间上的要求，那么就令加注标志"(/)=/，w(O的初始值为O;并记录下该次及该次以前的加注过程中各组储氢容器的取气质量及取气次序；最后求出数组"(O，中的最大值——即"，，z'=l:iVmax，并由程序控制器根据M皿对应的为各车载气瓶加注时各储氢容器的取气质量及取气次序向阀门控制装置输出控制信号。本发明还提供了一种用于前述加氢站高效加氢的三级加注优化控制方法的氢气加注系统，包括储氢容器、充气阀、压力传感器、氢气分配器、流量计和单向阀，充气阀、氢气分配器、单向阀依次连接至车载储气瓶；所述储氢容器是三级压力的储罐，包括低压储氢容器、中压储氢容器和高压储氢容器，分别通过各自的控制阀门cvl、cv2、cv3连接至充气阀；还包括一个程序控制器，压力传感器连接氢气分配器和程序控制器，流量计的测量点位于氢气分配器和单向阀之间，并通过信号线连接至程序控制器，程序控制器通过信号线与低压储氢容器、中压储氢容器和高压储氢容器各自的控制阀门cvl、cv2、cv3连接。与现有技术相比，本发明的有益效果是针对氢加注系统的特点，提出了过程优化控制方法，这是在以前的加注系统设计及研究中所没有的；利用本控制方法可对氢气的利用率和加注时间同时进行优化，能确保高的氢气利用率的同时，较大地縮短加注时间；可以根据不同的加注参数及要求，产生不同的优化加注方式；使加注过程简单、直接并易于操作。图1是加注系统优化控制方法流程图；图2是加氢站的加注系统流程示意图。图中附图标记1储气罐、2充气阀、3压力传感器、4氢气分配器、5流量计、6单向阀、7车载储气瓶、8程序控制器、低压储氢容器控制阀门cvl、中压储氢容器控制阀门cv2、高压储氢容器控制阀门cv3。具体实施例方式结合附图，下面通过具体实施例对本发明进行详细说明。以下根据加氢站具体的加注参数，用本控制方法实现的实例，并对优化后的具体加注过程进行了详细说明。本发明中的程序控制器8可选用西门子公司生产的PLC型号为S7-400的A口广叩o实例1对加氢站进行三级加注优化计算，其中低压储氢容器容积^=0.6m3，压力^:42M尸a;中压储氢容器容积F"2=0.4m3，压力P2=45MP";高压储氢容器容积^=0.2w3,压力P3=70Mi^;设定车载气瓶容积7=0.15m3，充满气后的压力尸=3571/尸"，每次充气前的初始压力P。=2MPa。根据
发明内容中的控制方法，进行编程求解，获得控制方法的非劣解。而后，根据过程优化的目的和实际要求从其中选出最优解，编制成具体的实际操作控制步骤。其中控制方法的非劣解如表1所示。<table>tableseeoriginaldocumentpage11</column></row><table><table>tableseeoriginaldocumentpage12</column></row><table><table>tableseeoriginaldocumentpage13</column></row><table>其中，f为本加注方式下，充满一辆车所用的最长的时间。由上表可以根据不同的加注要求，来选择不同的加注方式。为了同时达到高的氢气利用率和快速充气的目的，现以充气时间为242s，充气罐数为6罐的方式为最佳。其中，氢气利用率可以通过以下公式计算//=，=55.12%。M式中m为车载气瓶规定的充气质量，n为实际可充装的气瓶数量，M为各储氢容器的氢气总质量。氢气质量的计算方法附=PK;P为压力，单位尸fl;r为容积，单位"？；i=4124.3.3，r=293.15《，a=1.9155x10-6。将以上结果对应的加注方式编入程序控制器，通过程控器来实现对加注过程中阀门的切换和优化控制。根据本控制方法得到的优化加注方式的具体操作步骤如下a)开始第一辆车的充气，由程序控制器控制开启阀门2，并打开cvl先利用低压组的储气罐进行充气，当流量计5的流量达到3.0750紐时，低压组储气罐的阀门cvl关闭，中压组储气罐cv2的阀门打开，此时当流量计5测试的流量达到0.2141&时，切断中压组储气罐cv2的阀门，此时车载储气瓶充满。关闭阀门2完成第一辆车的充气过程；b)开始第二辆车的充气，由程序控制器控制开启阀门2，并打开cvl先利用低压组的储气罐进行充气，当流量计5的流量达到2.4600紐时，低压组储气罐的阀门cvl关闭，中压组储气罐cv2的阀门打开，此时当流量计5测试的流量达到0.8291&时，切断中压组储气罐cv2的阀门，此时车载储气瓶充满。关闭阀门2完成第二辆车的充气过程；c)接下来的车辆按照此方法，根据程序控制器设定好的程序来进行充气。实例2加氢站低压组固定容器容积V,二0.6m3，压力P!二42MPa，中压组固定容器容积V2=0.4m3，压力P2二70MPa，高压组固定容器容积V3=0.2m3，压力P3=70MPa，设定车载气瓶容积V二0.15m3，充满气后的压力P=35MPa，每次充气前的初始压力PG=2MPa。<table>tableseeoriginaldocumentpage13</column></row><table><table>tableseeoriginaldocumentpage14</column></row><table><table>tableseeoriginaldocumentpage15</column></row><table>根据本控制方法，可以得到以上几种加注方式。根据具体的加注要求可确定最佳方式，如当充气为8个罐的时候，其氢气利用率为;7=^^=65.47%，而当充气为7个罐的时，氢气利用率为77=^1=57.92%。为了同时具有较短的充气时间和高取气率，优先选用充气时间为181s，利用率为57.92%的加注方式，其与平衡加注方式相比可节省充气时间近50%。另外，在利用率均为57.92%的几种加注方式中，可以明显地观察到，所用时间最短，即181s，那种加注方式中，存在加注罐顺序的改变。如在此加注方式中，在充装第6罐的时候不是采用低-中-高的加注顺序，而是直接用中-高充气，到第7罐又恢复到低-中-高的加注顺序。可见，由于存在这种加注次序的自动切换，可以在相同利用率的情况下使其加注时间明显縮短。由于现在给出的实例中，只是系统一轮次加注过程，实际中，通过压縮机补气，釆用的是多轮次加注。因而，当加注车辆更多时，采用本方法来改变加注次序并同时实现高利用率和缩短加注时间的效果将会更加突出。最后，需要注意的是，以上列举的仅是本发明的具体实施例。显然，本发明不限于以上实施例，还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形，均应认为是本发明的保护范围。权利要求1、一种加氢站高效加氢的三级加注优化控制方法，包括以下步骤(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代表加注开始时的质量流量；id="icf0001"file="A2008100635840002C1.tif"wi="4"he="4"top="224"left="27"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>代表加注过程中车载气瓶与先取气的储氢容器的压力比；id="icf0002"file="A2008100635840002C2.tif"wi="14"he="4"top="232"left="27"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>分别代表先取气的储氢容器加注过程中达到等熵流动临界状态时及摩擦流动临界状态时车载气瓶与储氢容器的压力比；ms、mfs分别代表加注过程中达到等熵流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量、加注过程中达到摩擦流动临界状态时先取气的储氢容器为车载气瓶充入的氢气质量；t1为先取气的储氢容器向车载气瓶加注质量为x的氢气所用的时间；x1充为某时间下，先取气的储氢容器已经为车载气瓶充入的氢气质量；tmin代表整个加注过程只充装一个车载气瓶时，所用的最短时间；tmax代表不考虑加注时间约束时，氢气利用率最高的加注过程中最慢一次加注所用时间；Nmax代表不考虑加注时间约束时，所能充装的最大车载气瓶数；i表示加注第i个车载气瓶的标志；(2)将区间tmin至tmax之间按照等分的原则取r个值，确定r个约束时间tk，对每一个约束时间tk进行一次加注过程的优化，搜索获取各个约束时间下的氢气利用率最高的加注过程，即可充装的车载气瓶数最多的最优加注过程；(3)开始时i＝1，此时低压、中压、高压储氢容器的状态为第一次加注前的状态，根据以下氢气状态方程得到初始质量m11、m21、m31pV＝mRT(1+αp/T)(α＝1.9155×10-6，R＝4124.3)对第一个车载气瓶加注时，如果按低——中——高的取气次序取气，则低压储氢容器与车载气瓶达到压力平衡时充入的气体质量x1平衡为id="icf0003"file="A2008100635840003C1.tif"wi="34"he="10"top="149"left="27"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>低压储氢容器的取气质量范围为id="icf0004"file="A2008100635840003C2.tif"wi="40"he="6"top="150"left="124"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>此时将取气质量区间分为若干等分，x1的实际取值范围为各个等分点构成的集合；对于低压储氢容器的每一个取气质量值，中压储氢容器进行加注时，达到压力平衡时充入的氢气质量为id="icf0005"file="A2008100635840003C3.tif"wi="45"he="10"top="182"left="63"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>中压储氢容器的取气质量范围为id="icf0006"file="A2008100635840003C4.tif"wi="48"he="6"top="198"left="19"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>此时将取气质量区间分为若干等分，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-id="icf0007"file="A2008100635840004C1.tif"wi="3"he="3"top="103"left="66"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>曲线，然后对曲线按最小二乘法进行二次拟合，最后，通过拟合得到的多项式计算等熵流动状态下的质量流量；当压力比id="icf0008"file="A2008100635840004C2.tif"wi="2"he="3"top="119"left="44"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>小于或等于实际管路系统的临界压比id="icf0009"file="A2008100635840004C3.tif"wi="4"he="3"top="119"left="117"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>时，质量流量恒定，按摩擦流动方程计算；当压力比id="icf0010"file="A2008100635840004C4.tif"wi="2"he="3"top="127"left="48"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>大于摩擦流动临界状态时的压力比id="icf0011"file="A2008100635840004C5.tif"wi="3"he="3"top="127"left="116"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>时，流量按如下公式计算qf＝Cdq；对于先取气的储氢容器来说，如果加注开始时id="icf0012"file="A2008100635840004C6.tif"wi="14"he="4"top="134"left="112"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>充气结束时id="icf0013"file="A2008100635840004C7.tif"wi="14"he="4"top="134"left="152"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>那么先求出该储氢容器的初始压力和加注结束后的压力，然后求出初始质量流量qa和加注结束时的质量流量qb，<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>x</mi></mrow><mrow><msub><mi>q</mi><mi>a</mi></msub><mo>+</mo><msub><mi>q</mi><mi>b</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math>id="icf0014"file="A2008100635840004C8.tif"wi="18"he="10"top="156"left="35"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>如果加注开始时id="icf0015"file="A2008100635840004C9.tif"wi="13"he="4"top="159"left="87"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>加注结束时id="icf0016"file="A2008100635840004C10.tif"wi="20"he="4"top="159"left="124"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>那么分两段求取时间t1，求出id="icf0017"file="A2008100635840004C11.tif"wi="9"he="3"top="177"left="27"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>时的质量流量qfs及加注结束时的质量流量qb，<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mrow><mn>2</mn><mi>m</mi></mrow><mi>fs</mi></msub><mrow><msub><mi>q</mi><mi>a</mi></msub><mo>+</mo><msub><mi>q</mi><mi>fs</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>m</mi><mi>fs</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>q</mi><mi>fs</mi></msub><mo>+</mo><msub><mi>q</mi><mi>b</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math>id="icf0018"file="A2008100635840004C12.tif"wi="35"he="10"top="173"left="127"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>如果加注开始时id="icf0019"file="A2008100635840004C13.tif"wi="12"he="4"top="190"left="23"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>加注结束时id="icf0020"file="A2008100635840004C14.tif"wi="12"he="3"top="190"left="61"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>那么分四段求取充气时间，其中，计算中所用的id="icf0021"file="A2008100635840004C15.tif"wi="3"he="3"top="191"left="167"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>及ms为初始压力pi下求出的等熵流动状态下的临界压力比；将后面的时间质量流量曲线按充气质量等分的原则，处理为两段直线，<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mrow><mn>2</mn><mi>m</mi></mrow><mi>fs</mi></msub><mrow><msub><mi>q</mi><mi>a</mi></msub><mo>+</mo><msub><mi>q</mi><mi>fs</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><msub><mi>m</mi><mi>s</mi></msub><mo>-</mo><msub><mi>m</mi><mi>fs</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>q</mi><mi>fs</mi></msub><mo>+</mo><msub><mi>q</mi><mi>s</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>m</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><msub><mi>q</mi><mi>s</mi></msub><mo>+</mo><msub><mi>q</mi><mi>b</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>m</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><msub><mi>q</mi><mi>b</mi></msub><mo>+</mo><msub><mi>q</mi><mi>c</mi></msub></mrow></mfrac><mo>,</mo></mrow>]]></math>id="icf0022"file="A2008100635840004C16.tif"wi="71"he="10"top="212"left="78"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中，qb为等分点x1充＝ms+0.5(x-ms)时的质量流量，qc为加注结束时的质量流量；如果加注开始时id="icf0023"file="A2008100635840004C17.tif"wi="19"he="4"top="242"left="19"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>加注结束时id="icf0024"file="A2008100635840004C18.tif"wi="20"he="4"top="242"left="64"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/><mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>x</mi></mrow><mrow><msub><mi>q</mi><mi>a</mi></msub><mo>+</mo><msub><mi>q</mi><mi>b</mi></msub></mrow></mfrac><mo>,</mo></mrow>]]></math>id="icf0025"file="A2008100635840004C19.tif"wi="19"he="10"top="240"left="90"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>qa及qb分别为加注开始和加注结束时的质量流量；如果加注开始时id="icf0026"file="A2008100635840004C20.tif"wi="23"he="4"top="256"left="70"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>加注结束时id="icf0027"file="A2008100635840004C21.tif"wi="12"he="4"top="256"left="120"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>那么分三段计算充气时间，<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mrow><mn>2</mn><mi>m</mi></mrow><mi>s</mi></msub><mrow><msub><mi>q</mi><mi>a</mi></msub><mo>+</mo><msub><mi>q</mi><mi>s</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>m</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><msub><mi>q</mi><mi>s</mi></msub><mo>+</mo><msub><mi>q</mi><mi>b</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>m</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><msub><mi>q</mi><mi>b</mi></msub><mo>+</mo><msub><mi>q</mi><mi>c</mi></msub></mrow></mfrac><mo>,</mo></mrow>]]></math>id="icf0028"file="A2008100635840004C22.tif"wi="49"he="10"top="267"left="20"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>qb为等分点的质量流量，qc为加注结束时的质量流量；如果，加注开始时id="icf0029"file="A2008100635840005C1.tif"wi="11"he="4"top="32"left="40"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>那么<mathsid="math0006"num="0006"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>x</mi></mrow><mrow><msub><mi>q</mi><mi>a</mi></msub><mo>+</mo><msub><mi>q</mi><mi>b</mi></msub></mrow></mfrac><mo>,</mo></mrow>]]></math>id="icf0030"file="A2008100635840005C2.tif"wi="19"he="10"top="29"left="64"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>qb为加注结束时的质量流量；后取气的两个储氢容器加注时间的计算方法与先取气的储氢容器相同，对此车载气瓶总的加注时间为分别从三个储氢容器中取气的加注时间之和；(6)选择加注时间最短的那种取气次序判断最短的加注时间是否小于或等于约束时间，即是否满足ti≤tk；如果满足，则令加注标志n(i)＝i，根据氢气的状态方程获取各储氢容器加注后的状态参数作为加注第二个车载气瓶时的初始状态，计算出各组储氢容器取气质量的取值范围，进行第二次加注的运算；如果不满足，则继续第一次加注过程的循环搜索，搜索此加注过程中其他的取气质量组合，作同样的判断；(7)根据第一次加注完毕获得出的加注第二个车载气瓶时的初始状态以及各储氢容器取气质量的取值范围，进行第二次加注的搜索，其方法与第一次加注的搜索相同；(8)由前面的加注过程逐步向后搜索直到所有的循环搜索全部完成，对于第i个车载气瓶加注过程中的某种组合如果能够充满车载气瓶并满足时间上的要求，那么就令加注标志n(i)＝i，n(i)的初始值为0；并记录下该次及该次以前的加注过程中各组储氢容器的取气质量及取气次序；最后求出数组n(i)，中的最大值——即nmax，i＝1Nmax，并由程序控制器根据nmax对应的为各车载气瓶加注时各储氢容器的取气质量及取气次序向阀门控制装置输出控制信号。2、一种用于权利要求1所述加氢站高效加氢的三级加注优化控制方法的氢气加注系统，包括储氢容器、充气阔、压力传感器、氢气分配器、流量计和单向阀，充气阀、氢气分配器、单向阀依次连接至车载储气瓶；其特征在于，所述储氢容器是三级压力的储罐，包括低压储氢容器、中压储氢容器和高压储氢容器，分别通过各自的控制阀门(cvl)、(cv2)、(cv3)连接至充气阀；还包括一个程序控制器，压力传感器连接氢气分配器和程序控制器，流量计的测量点位于氢气分配器和单向阀之间，并通过信号线连接至程序控制器，程序控制器通过信号线与低压储氢容器、中压储氢容器和高压储氢容器各自的控制阀门(cvl)、(cv2)、(cv3)连接。全文摘要本发明涉及一种加氢优化控制方法及其系统，旨在提供一种加氢站高效加氢的三级加注优化控制方法，即加氢站中储气罐分为高中低三组；为了同时实现高取气率和快速加注的目标，提出了多目标优化算法，采用约束法将多目标问题转化为若干个单目标优化问题；每个单目标规划问题属于动态规划范筹，采用由前向后的顺序搜索方法；根据优化结果设计出具体的程序控制流程，为加氢站的加注系统设计出最佳的加注方式。本发明还提供了一种用于前述加氢站高效加氢的三级加注优化控制方法的氢气加注系统，包括储氢容器、充气阀、压力传感器、氢气分配器、流量计和单向阀；储氢容器包括低压储氢容器、中压储氢容器和高压储氢容器；还包括一个程序控制器。文档编号G05B13/02GK101315545SQ20081006358公开日2008年12月3日申请日期2008年6月27日优先权日2008年6月27日发明者叶建军,萍唐,平徐,健杨,磊赵,郑津洋申请人:浙江大学