[0001] 本实用新型涉及海上风能利用技术领域，特别是指一种利用海上风电制氢的地下油气藏储氢系统。

[0002] 在海上风能利用和氢能开发的过程中，由于海上风能具有逐时性、不稳定的特点，海上风能产生的电能在时间维度上具有明显的峰谷特征，不仅不适合于并网供电，更容易造成弃风弃电的能源浪费现象。对海上风力生产的电能进行储蓄调峰是解决上述问题的有效方法。另外，氢能不仅是一种清洁无污染、高能量密度的能源，还是一种良好的储能载体，将海上风电以氢能的形式进行储蓄，能实现能源高效利用和稳定输出，以及环境有效保护。相比起传统的储气罐储气以及目前借助有机液体的海上储氢方式，地下储气库储气空间大，能够满足海上风电制氢的大量储存需求。具体来看，海底油气藏储气不仅安全性高，还能节省大量的储气材料，是一种理想的氢气地下储存方式。此外，相比起抽水蓄能和氢气液化存储，海底废弃油气藏储氢还能节省大量的工程建设和经济成本支出。

[0003] 常规的能源系统的效率或产能可根据已知的初始值直接计算得到。然而，对于海上风电制氢地下油气藏储氢系统，由于发电机组生产的电能除了供给海水电解制氢设备外，还供给注气压缩(采气提纯)设备使用。这其中存在一个用电分配比的问题，且该分配比未知，意味着系统的产能和储(采)气量无法直接计算求得，需要通过迭代循环过程计算求解。迭代求解法是一种不断用变量的旧值递推新值，最终得到满足实际条件的解的过程。利用迭代求解法，以天为最小时间步长，能够解决由于系统用电分配百分比未知造成的计算困难，有助于方便、快速、准确地调控海上风电制氢地下油气藏储氢系统的日氢气产能与日储(采)气量。

[0004] 系统产能和储(采)气数量是项目方关心的两个重要参数。然而，目前针对能源系统的计算方法基本只关注于其中一个指标，无法全面地反映系统的调控能力和产出能力。一套综合考虑系统日氢气产能和日储(采)气量的计算方法，能为项目方提供判断系统调控能力和产出能力优劣的参考指标。
实用新型内容

[0005] 本实用新型要解决的技术问题是提供一种利用海上风电制氢的地下油气藏储氢系统。

[0006] 该系统包括海上风电机组、海水电解装置、氢气压缩设备、氢气解压设备、氢气提纯设备、注/采气井，油气藏储库和氢气输出设备，氢气输出设备包括运氢船和氢气管网，海上风电机组通过海底电缆与海水电解装置、氢气压缩设备和氢气提纯设备相连，海水电解装置、氢气压缩设备、氢气解压设备、氢气提纯设备设置在海洋底部，海水电解装置通过输氢管道连接氢气压缩设备，氢气压缩设备通过输氢管道连接缓冲罐，缓冲罐和氢气解压设备之间设置三通阀，氢气解压设备通过输氢管道连接氢气提纯设备，氢气提纯设备连接氢气管网，三通阀设置在注/采气井的井入口，注/采气井开在海床上，海床底部设置油气藏储库，油气藏储库连通注/采气井，油气藏储库上部设置氢气储库上覆盖层岩石，油气藏储库下部设置氢气储库底层岩石，用以维持储库的压力和良好的气密性。

[0007] 其中，海上风电机组按照矩形或圆形排列，受风电机叶片长度以及尾流效应的影响，各风电机之间的半径距离为800-1000米，各风电机的型号相同。

[0008] 海水电解装置包括海水脱盐预处理装置和电解水制氢电解槽，海水电解装置通过海底电缆从海上风电机组输入电能。

[0009] 该系统的调控计算方法，包括步骤如下：

[0010] (1)从选定海域的历史数据中获取连续的τ个离海平面h米日均风速数据，构成日均风速数据集V，然后根据海上风力发电机的风速-输出功率函数求得该海域的单台风力发电机日输出电功率集P，再根据风力发电机的日输出电功率集P、风力发电机的数量和工作时间求得风力发电场日发电量集W；

[0011] (2)引入系统日用电分配百分比集Ι，设置电解海水制氢的日用电占比为ηi，i表示第i天，压缩储气或采气提纯的日用电占比为1-ηi，并为ηi赋初值，其中，0<ηi<1，日用电分配百分比集Ι＝[η1,η2,......,ητ]；

[0012] (3)根据风力发电机日发电量集W以及电解海水制氢的用电分配百分比集Ι，两者对应元素相乘求得电解制氢日耗电量集We，并结合法拉第定律求得系统日电解产氢量集Y；

[0013] (4)对日电解产氢量集Y内的τ个日电解产氢量元素求平均值，获得系统日氢气产能ysd，然后将日电解产氢量集Y中的元素yi减去日氢气产能ysd，并判断：若结果非负则第i天进行压缩储气流程，若结果为负则第i天进行采气提纯流程；

[0014] (5)根据实际流程，计算压缩储气流程的日能耗wc,i或采气提纯流程的日能耗wp,i，并判断：若总能耗值we,i+wc,i或we,i+wp,i与风电机日发电量wi相差小于ε，则当前用电分配比满足要求，输出日氢气产能ysd和日储气或日采气量yi-ysd，并求出系统的平均氢气储产比ζ，即为表征系统氢气产出能力与调控能力的耦合指标，其中，we,i为第i天的电解制氢日耗电量集；若总能耗值we,i+wc,i或we,i+wp,i与风电机日发电量wi相差不小于ε，则进入步骤(6)进一步判断；

[0015] (6)若总能耗值we,i+wc,i或we,i+wp,i大于等于风电机日发电量wi，则当前用电分配比减小为原来的50％，若总能耗值we,i+wc,i或we,i+wp,i小于风电机日发电量wi，则当前用电分配比增大为原来的150％，调整用电分配比完成后，返回步骤(4)进行迭代计算。

[0016] 步骤(1)中风力发电场日发电量集W计算方法如下：

[0017] W＝n×T*P＝[w1,w2,......,wτ]

[0018] 其中，n是风电场所含风电机数量，单位为台，T是风电机的日工作时间集，wi是单台风电机第i日的日发电量，单位为千瓦时，其中，i＝1,2,……,τ。

[0019] 步骤(3)中系统日电解产氢量集Y计算方法如下：

[0020]

[0021] 其中，yi是该海域第i日的日电解产氢量，单位为千克/天，i＝1,2,……,τ；We是电解制氢日耗电量集，m是海水电解制氢产能，单位为千瓦时/千克。

[0022] 系统每天的产氢量会因当天风速不同而改变，然而氢气每日的市场需求是稳定的。为了系统日氢气产能保持稳定，步骤(4)根据τ个日电解产氢量求出平均日电解产氢量，作为日氢气产能。需要注意的是，在步骤(4)所述，求得的日氢气产能只是迭代循环过程的中间值，而不是最终值。若当日电解产氢量大于日氢气产能，多余产量进行压缩储气；若当日电解产氢量小于日氢气产能，不足产量进行采气提纯补充。

[0023] 步骤(5)中通过迭代过程调整系统的用电分配百分比，最终获得系统最优用电分配百分比，即在该分配百分比下，电解制氢能耗与注气压缩(采气提纯)能耗之和与风电机输入电能相差小于ε。

[0024] 步骤(5)中压缩储气流程的日能耗wc,i计算过程如下：

[0025]

[0026] 式中：wc,i是第i日的日压缩储气能耗，单位为焦/天；γ是氢气的比热容比，γ＝1.41，无量纲；f0是压缩系统进口氢气的压力，单位为兆帕；ρ0是压缩系统进口氢气的密度，单位为千克/立方米；fwh是压缩系统出口压力，同时是注采井井口压力，单位为兆帕；

[0027] 采气提纯流程的日能耗wp,i计算过程如下：

[0028]

[0029] 式中：wp,i是第i日的日采气提纯能耗，单位为千焦/小时；T0是采气温度，单位为开；R是通用气体常数，R＝8.314千焦/(千克·开)，β1、β2、βF分别是氢气提纯分离过程中进料气、解析气和产物气中氢气的质量分数。

[0030] 步骤(5)中ε为迭代终止的判断值，根据实际情况设置；

[0031] 系统的平均氢气储产比ζ计算过程如下：

[0032]

[0033] 式中：ysd是系统日氢气产能，单位为千克/天；yi-ysd是系统第i日的日储气量或采气量，单位为千克/天；ζ是系统平均氢气储产比，无量纲。

[0034] 步骤(5)中若所得系统的平均氢气储产比ζ＞1，则表示该系统调控能力较为突出，且ζ数值越大，系统调控能力越强；若ζ＜1，则表示系统产出能力更强，且ζ数值越接近于0，系统产出能力越强。

[0035] 本实用新型的上述技术方案的有益效果如下：

[0036] 上述方案中，能够有效解决目前海上风电输出不稳定、浪费严重和储能成本高、储量低的痛点，充分利用地下空间资源，推广应用高能量密度无环境污染的氢能源，实现能源高效利用和稳定输出，以及环境有效保护。

[0037] 所述方法中，利用迭代求解法，以天为最小时间步长，有效避免计算过程中系统用电分配百分比未知导致无法直接计算的困难，能够方便、快速、准确地调控海上风电制氢地下储氢系统的日氢气稳定输出量和日氢气储(采)气量。

[0038] 此外，针对该系统产能和储(采)气量这两个项目方关心的关键指标，本实用新型结合该系统日稳定产能和日氢气储(采)气量，求解该系统平均氢气储产比，来综合反映该系统的产出能力和调控能力。

[0042] 为使本实用新型要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。

[0043] 本实用新型提供一种利用海上风电制氢的地下油气藏储氢系统。

[0044] 如图1所示，该系统包括海上风电机组1、海水电解装置3、氢气压缩设备5、氢气解压设备8、氢气提纯设备9、注/采气井17，油气藏储库16和氢气输出设备，氢气输出设备包括运氢船10和氢气管网12，氢气管网12布置在陆地11上，海上风电机组1通过海底电缆2与海水电解装置3、氢气压缩设备5和氢气提纯设备9相连，海水电解装置3、氢气压缩设备5、氢气解压设备8、氢气提纯设备9设置在海洋13底部，海水电解装置3通过输氢管道4连接氢气压缩设备5，氢气压缩设备5通过输氢管道4连接缓冲罐6，缓冲罐6和氢气解压设备8之间设置三通阀7，氢气解压设备8通过输氢管道4连接氢气提纯设备9，氢气提纯设备9连接氢气管网12，三通阀7设置在注/采气井17的井入口，注/采气井17开在海床14上，海床14底部设置油气藏储库16，油气藏储库16连通注/采气井17，油气藏储库16上部设置氢气储库上覆盖层岩石15，油气藏储库16下部设置氢气储库底层岩石18。

[0045] 如图2所示，该系统调控计算的方法，包括步骤如下：

[0046] (1)从选定海域的历史数据中获取连续的τ个离海平面h米日均风速数据，构成日均风速数据集V，然后根据海上风力发电机的风速-输出功率函数求得该海域的单台风力发电机日输出电功率集P，再根据风力发电机的日输出电功率集P、风力发电机的数量和工作时间求得风力发电场日发电量集W；

[0047] (2)引入系统日用电分配百分比集Ι，设置电解海水制氢的日用电占比为ηi，i表示第i天，压缩储气或采气提纯的日用电占比为1-ηi，并为ηi赋初值，其中，0<ηi<1，日用电分配百分比集Ι＝[η1,η2,......,ητ]；

[0048] (3)根据风力发电机日发电量集W以及电解海水制氢的用电分配百分比集Ι，两者对应元素相乘求得电解制氢日耗电量集We，并结合法拉第定律求得系统日电解产氢量集Y；

[0049] (4)对日电解产氢量集Y内的τ个日电解产氢量元素求平均值，获得系统日氢气产能ysd，然后将日电解产氢量集Y中的元素yi减去日氢气产能ysd，并判断：若结果非负则第i天进行压缩储气流程，若结果为负则第i天进行采气提纯流程；

[0050] (5)根据实际流程，计算压缩储气流程的日能耗wc,i或采气提纯流程的日能耗wp,i，并判断：若总能耗值we,i+wc,i或we,i+wp,i与风电机日发电量wi相差小于ε，则当前用电分配比满足要求，输出日氢气产能ysd和日储气或日采气量yi-ysd，并求出系统的平均氢气储产比ζ，即为表征系统氢气产出能力与调控能力的耦合指标，其中，we,i为第i天的电解制氢日耗电量集；若总能耗值we,i+wc,i或we,i+wp,i与风电机日发电量wi相差不小于ε，则进入步骤(6)进一步判断；

[0051] (6)若总能耗值we,i+wc,i或we,i+wp,i大于等于风电机日发电量wi，则当前用电分配比减小为原来的50％，若总能耗值we,i+wc,i或we,i+wp,i小于风电机日发电量wi，则当前用电分配比增大为原来的150％，调整用电分配比完成后，返回步骤(4)进行迭代计算。

[0052] 下面结合具体实施例予以说明。

[0053] 在实际设计中，该系统的能量输入来源于海上风能，为了实现紧凑的系统结构和高效的能量传输，体现就地取能用能的能源系统设计原则，海上风电机组、海水电解装置、氢气压缩设备、氢气解压设备和氢气提纯设备都需要在海上进行安装，具体地，安装在海床平面上。同样地，油气藏型地下储氢库是位于海床地下的油气藏，注(采)气井是在海床地下油气藏上钻井建设的气井。

[0054] 海上风电机组按照矩形或圆形排列；其中，受风电机叶片长度以及尾流效应的影响，各风电机之间的半径距离通常在800米至1000米。

[0055] 海水电解装置还包括用于海水脱盐预处理装置和用于电解水制氢的电解槽。其中，海水电解装置通过海底电缆与海上风电机组相连，通过海底电缆从海上风电机组输入电能。

[0056] 氢气压缩设备还包括用于稳定气压的缓冲罐。其中，氢气压缩设备通过海底电缆从海上风电机组输入电能。其中，氢气压缩设备的入口与海水电解装置出口相连。

[0057] 其中，注/采气井还包括用于调节和引导氢气流向的三通阀；其中，三通阀的三个口分别与氢气压缩设备的出口，注/采气井井口和氢气解压设备入口相连接。

[0058] 其中，氢气提纯设备通过海底电缆从海上风电机组输入电能；其中，氢气提纯设备的入口与氢气解压设备的出口相连。

[0059] 其中，油气藏型地下储氢库还包括了用于密封的上覆盖层岩石和底层岩石，用以维持储库的压力和良好的气密性；其中，储气库的储气区域与注/采气井井底相连。

[0060] 其中，氢气输出设备还包括用于氢气海上运输的运氢船和用于氢气陆上运输的氢气管网；其中，氢气输出设备与氢气提纯设备的出口相连。

[0061] 在具体应用中，系统按如下过程运行：

[0062] (1)首先，海上风流过海上风电机组1的叶片使其转动，带动风电机内的涡轮切割磁场发电。

[0063] (2)海上风电机组1生产的电能通过海底电缆2输出到海水电解装置3、氢气压缩设备5和氢气提纯设备9，为上述装置进行供电。

[0064] (3)海水电解装置3内，海水经过脱盐预处理后发生电解反应，生成氢气。

[0065] (4)当日额定的氢气产量依次经过氢气压缩设备5、三通阀7、氢气解压设备8和氢气提纯设备9，最后进入氢气输出设备向用户输出。

[0066] (5)当日多余的氢气产量依次经过氢气压缩设备5、三通阀7、注/采气井17进入油气藏储库16进行储存。

[0067] (6)当日不足的氢气产量依次经过油气藏储库16、注/采气井17、三通阀7、氢气解压设备8和氢气提纯设备9，最后进入氢气输出设备向用户输出，实现采气补给。

[0068] 在具体应用中，该系统按如下过程进行调控：

[0069] (1)从选定海域的历史数据中获取连续的τ个离海平面h米日均风速数据集V。其中V列于下方：

[0070] V＝[v1,v2,......,vτ]   (1)

[0071] 其中，vi(i＝1,2,……,τ)是选定海域的历史数据中第i日的离海平面h米日均风速，米/秒；

[0072] (2)根据海上风力发电机的风速-输出功率函数求得该海域的单台风力发电机日输出电功率集P。其中，单台风力发电机日输出电功率集P列于下方：

[0073] P＝[p1,p2,......,pτ]   (2)

[0074] 其中，pi(i＝1,2,……,τ)是单台风力发电机第i日的日输出电功率，千瓦；

[0075] 其中，风速-输出功率函数列于下方：

[0076]

[0077] 其中，vcut-in是风电机切入风速，vcut-out是风电机切出风速，vr是风电机额定风速，米/秒，pr是风电机额定输出功率，千瓦。

[0078] 其中，风电机型号可根据该海域最低日均风速进行选择，选择的原则是该海域最低日均风速不小于风电机切入风速。

[0079] (3)根据风力发电机的日输出电功率，数量和工作时间可求得风力发电场日发电量集W。其中W列于下方：

[0080] W＝n×T*P＝[w1,w2,......,wτ]   (4)

[0081] 其中，n是风电场所含风电机数量，台，T是风电机的日工作时间集，wi(i＝1,2,……,τ)是单台风电机第i日的日发电量，千瓦时，其中T列于下方：

[0082] T＝[t1,t2,......,tτ]   (5)

[0083] 其中，ti(i＝1,2,……,τ)是单台风电机第i日的日工作时间，小时，“*”表示两个等维向量之间对应元素相乘；

[0084] (4)引入用电分配百分比，其中海水电解槽的用电分配百分比为η，相应地当日压缩注入氢气(采出提纯氢气)的用电分配百分比为1-η，并设置电解海水制氢的用电占比和压缩注入氢气/采出提纯氢气的用电分配比初始值各为50％。其中，海水电解的日用电分配百分比集Ι列于下方：

[0085] Ι＝[η1,η2,......,ητ]   (6)

[0086] 其中，压缩注入氢气/采出提纯氢气的用电分配百分比集[1]-Ι列于下方：

[0087] [1]-Ι＝[(1-η1),(1-η2),......,(1-ητ)]   (7)

[0088] 其中，[1]表示与Ι等维的元素全为1的矩阵。

[0089] (5)根据风力发电机日发电量以及电解海水制氢的用电分配比，结合法拉第定律求得日电解产氢量集Y。其中，海水电解制氢产能表达式列于下方：

[0090]

[0091] 其中，F是法拉第常数，F＝9.65×104库伦/摩尔；n是氢气所含元素的化合价绝对值之和，n＝2；Mmol是氢气的摩尔质量，Mmol＝2克/摩尔；U是输入电压，V，工程上一般为2.0伏；δ是考虑海水预处理能耗所计入的能耗系数，δ>1；m是海水电解制氢产能，千瓦时/千克。

[0092] 其中，Y列于下方：

[0093]

[0094] 其中，yi(i＝1,2,……,τ)是该海域第i日的日电解产氢量，千克/天，We是日电解制氢电耗集；

[0095] (6)对日电解产氢量集Y内的τ个日电解产氢量元素yi求平均值，获得系统日氢气产能ysd。其中，ysd列于下方：

[0096] ysd＝Average(Y)   (10)

[0097] 其中，“Average”表示对向量中所有元素求算术平均值；ysd是系统日氢气产能，千克/天；

[0098] (7)将日电解产氢量yi减去日氢气产能ysd，并判断：若结果非负则该日氢气产能盈余需要进行压缩储气流程，若结果为负则该日氢气产能亏损需要进行采气提纯流程；

[0099] (8)计算第i日的日压缩储气(采气提纯)流程的能耗wc,i(wp,i)，并判断：若系统总能耗we,i+wc,i(wp,i)与所分配用电wi相差小于ε，则当前用电分配比满足要求，进入步骤(10)，否则进入下一判断(9)。其中，井口压力fwh，压缩储气能耗wc,i和采气提纯能耗wp,i列于下方：

[0100]

[0101] 其中:fwh＝(1-α)fwf   (12)

[0102]

[0103]

[0104] 其中：wc,i是第i日的日压缩储气能耗，焦/天；γ是氢气的比热容比，γ＝1.41，无量纲；f0是压缩系统进口氢气的压力，兆帕；ρ0是压缩系统进口氢气的密度，千克/立方米；fwh是压缩系统出口压力，同时是注采井井口压力，兆帕；α是注采气损失百分比；fg是地下储库压力，fg,0是地下储库原始压力，兆帕；μ是储层氢气粘度，帕秒；L是储层厚度，米；K是储层岩石渗透率，毫达西；A是储库面积，平方米；Z是氢气在储层压力下的压缩因子，无量纲；R是通用气体常数，R＝8.314千焦/(千克·开)；Tg是储层温度，开；Vg是储库体积，立方米；

[0105] 其中：

[0106]

[0107] 其中：wp,i是第i日的日采气提纯能耗，千焦/小时；T0是采气温度，开；β1、β2、βF分别是氢气提纯分离过程中进料气、解析气和产物气中氢气的质量分数；

[0108] (9)若系统总能耗we,i+wc,i(wp,i)大于所分配用电wi，则当前用电分配比减小为原来的50％，若系统总能耗we,i+wc,i(wp,i)小于所分配用电wi，则当前用电分配比增大为原来的150％，调整用电分配比完成后，返回步骤(5)进行迭代计算；

[0109] (10)迭代循环计算结束，输出系统日氢气产能ysd，即为系统日氢气产能。输出系统日储(采)气量yi-ysd。计算系统平均氢气储产比ζ：

[0110]

[0111] 其中：ysd是系统日氢气产能，千克/天；yi-ysd是系统第i日的日注气量，千克/天；ζ是系统平均氢气储产比，无量纲。

[0112] 以上所述是本实用新型的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本实用新型所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本实用新型的保护范围。