(12)发明专利申请
(10)申请公布号 CN 110912140 A(43)申请公布日 2020.03.24
(21)申请号 201911284405.1(22)申请日 2019.12.13
(71)申请人 福州大学
地址 350108 福建省福州市闽侯县福州大
学城乌龙江北大道2号福州大学(72)发明人 江岳文 费迎阳
(74)专利代理机构 福州元创专利商标代理有限
公司 35100
代理人 钱莉 蔡学俊(51)Int.Cl.
H02J 3/06(2006.01)H02J 3/24(2006.01)
权利要求书4页 说明书10页 附图1页
(54)发明名称
一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法(57)摘要
本发明涉及一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，包括以下步骤：步骤S1：分别根据压缩机与电解槽的内部电压、电流关系得到压缩机功率模型与电解槽功率模型；步骤S2：求取制氢系统节点注入功率，并以此确定功率失配量方程，求取功率失配量对节点电压的一阶导数；步骤S3：采用牛顿法，建立直流配电网的雅可比矩阵，进行迭代计算，当迭代计算收敛后得到电解槽的功率，进而得到制氢系统功率值。本发明使制氢系统的内部电流参与到迭代计算中，及时修正节点电压量，与配电网运行相互协调，能够得到更准确的潮流分布结果。
CN 110912140 ACN 110912140 A
权 利 要 求 书
1/4页
1.一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，其特征在于，包括以下步骤：
步骤S1：分别根据压缩机与电解槽的内部电压、电流关系得到压缩机功率模型与电解槽功率模型；
步骤S2：求取制氢系统节点注入功率，并以此确定功率失配量方程，求取功率失配量对节点电压的一阶导数；
步骤S3：采用牛顿法，建立直流配电网的雅可比矩阵，进行迭代计算，当迭代计算收敛后得到电解槽的功率，进而得到制氢系统功率值。
2.根据权利要求1所述的一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，其特征在于，步骤S1具体为：
根据电解槽的内部电流电压关系得到所述电解槽功率模型：
其中，
式中，为电解槽功率，U0为驱动电解水反应的最小电压，为电解槽注入电流，R为
普适气体常数，T为电解槽工作温度，z为半反应电荷转移化学计量系数，F为法拉第常数，αa
为阳极电荷转移系数，A为电解槽横截面积，βtmem为膜的厚度，σ1的计算公式如上式所示，mem为膜的电导，joa为阳极交换电流密度，jL为极限电流密度，为氢气分压；
根据压缩机的进出口压强得到所述压缩机功率模型：
式中，Pcom,p为压缩机消耗的有功功率，ngas表示压缩机出口的气体流速，P1表示压缩机入口压强，P2表示压缩机出口压强，m表示压缩机运行在多变过程下的指数。
3.根据权利要求1所述的一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，其特征在于，步骤S2具体包括以下步骤：
步骤S21：求取制氢系统的节点注入功率：
式中，Pi为第i个节点的注入功率，PG表示节点发电机注入功率，Pload表示节点负荷功率，N为电解槽数目；
根据求取的制氢系统节点注入功率建立失配量方程：
2
CN 110912140 A
权 利 要 求 书
2/4页
ΔPi＝Pi-UiIi＝PG-Pload-UiIi；式中，ΔPi为第i个节点的功率失配量，Ui表示节点i的电压，Ii表示节点i的注入电流；其中Ii的表达式为：
式中，Yii表示节点i自导纳，Yij表示节点i与其相连节点互导纳，n表示直流配电网络中节点数目；由于电解槽工作电压低，通过直流换流器进行电压转换，YiiUE为该直流换流器的接入对制氢系统接入节点注入电流的修正量；
步骤S22：求取功率失配量对节点电压的一阶导数：
其中，
4.根据权利要求3所述的一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，其特征在于，步骤S2还包括以下步骤：
步骤S23：令压缩机消耗的有功功率Pcom,p为恒定值，其对节点电压的导数为0，令节点发电机注入功率PG同样为恒定值，其对节点电压的导数为0，因此仅需考虑电解槽，将步骤S22中的一阶导数公式转换为：
其中，
其中，电解槽注入电流为：
3
CN 110912140 A
权 利 要 求 书
3/4页
则，
式中，Icom为压缩机负荷注入节点电流，Icom＝-Pcom,p/Ui。
5.根据权利要求1所述的一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，其特征在于，步骤S3具体包括以下步骤：
步骤S31：形成直流配电网的雅可比矩阵，设制氢系统负荷接入直流配电网的节点i，配电网中含有n个节点，当不考虑负荷模型时，直流配电网络节点电压方程为：
Pi＝UiIi；
建立此时的雅可比矩阵元素为：
当考虑负荷模型后，建立的雅可比矩阵为：
步骤S32：采用牛顿拉夫逊法进行迭代计算，第k+1次迭代计算后，节点i电压新值为：式中，为节点i第k次迭代计算得到的节点电压值，到的电压修正量；
根据节点电压方程求出更新后的节点注入电流：
为节点i第k+1次迭代计算得
4
CN 110912140 A
权 利 要 求 书
4/4页
其中，Yij为节点i与其相连节点的互导纳，根据公式得到电解槽注入电流：
步骤S33：根据电解槽功率模型，更新功率值，若迭代计算没有收敛，则返回步骤S31再次求取雅可比矩阵，进入下一次迭代，若迭代计算收敛，则进入步骤S34；
步骤S34：得到当前的电解槽注入电流，并根据该电解槽注入电流求取制氢系统的功率值；
步骤S35：进一步得到制氢系统的功率值Pload：
6.一种计及电解水制氢系统精细化模型的直流配电网潮流计算系统，其特征在于，包括存储器、处理器以及存储于存储器上并能够被所述处理器运行的计算机程序，当处理器运行该计算机程序时，实现如权利要求1-5任一项所述的方法步骤。
5
CN 110912140 A
说 明 书
1/10页
一种计及电解水制氢系统精细化模型的直流配电网潮流计算
方法
技术领域
[0001]本发明涉及电力系统技术领域，特别是一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法。
背景技术
[0002]电解水制氢系统作为一种可控负荷，理想情况下，可在额定功率的0-100％范围内运行，对于提高功率波动较大的风电、光伏发电利用率具有重大意义。针对功率波动较大的风电、光伏发电并网难消纳问题，研究者提出了非并网理论，应用风电的同时，绕开电网这一瓶颈，使风电、光伏发电与蓄电池以及制氢系统相互配合运行，提高风电利用率。然而，风电、光伏发电受风力、光照等自然因素影响较大，电压值、功率值波动较大，因此，需获悉不同状态下电解槽内的潮流值，对电解槽的工作状态进行评估。
[0003]当前文献将制氢系统作为恒电流负荷或恒功率负荷参与到系统潮流计算及控制策略的研究中，无法观察电解槽内温度、压力等物理参数的变化对配电网络潮流分布的影响。
发明内容
[0004]有鉴于此，本发明的目的是提出一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，使制氢系统的内部电流参与到迭代计算中，及时修正节点电压量，与配电网运行相互协调，能够得到更准确的潮流分布结果。[0005]本发明采用以下方案实现：一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，包括以下步骤：[0006]步骤S1：分别根据压缩机与电解槽的内部电压、电流关系得到压缩机功率模型与电解槽功率模型；[0007]步骤S2：求取制氢系统节点注入功率，并以此确定功率失配量方程，求取功率失配量对节点电压的一阶导数；[0008]步骤S3：采用牛顿法，建立直流配电网的雅可比矩阵，进行迭代计算，当迭代计算收敛后得到电解槽的功率，进而得到制氢系统功率值。[0009]进一步地，步骤S1具体为：
[0010]根据电解槽的内部电流电压关系得到所述电解槽功率模型：
[0011][0012]
其中，
6
CN 110912140 A
说 明 书
2/10页
[0013]
[0014][0015][0016]
式中，为电解槽功率，U0为驱动电解水反应的最小电压，为电解槽注入电流，
R为普适气体常数，T为电解槽工作温度，z为半反应电荷转移化学计量系数，F为法拉第常
数，αA为电解槽横截面积，βtmem为膜的厚a为阳极电荷转移系数，1的计算公式如上式所示，度，σjoa为阳极交换电流密度，jL为极限电流密度，为氢气分压；mem为膜的电导，
[0017][0018][0019]
根据压缩机的进出口压强得到所述压缩机功率模型：
式中，Pcom,p为压缩机消耗的有功功率，ngas表示压缩机出口的气体流速，P1表示压缩机入口压强，P2表示压缩机出口压强，m表示压缩机运行在多变过程下的指数。[0020]进一步地，步骤S2具体包括以下步骤：[0021]步骤S21：求取制氢系统的节点注入功率：
[0022]
式中，Pi为第i个节点的注入功率，PG表示节点发电机注入功率，Pload表示节点负荷功率，N为电解槽数目；
[0024]根据求取的制氢系统节点注入功率建立失配量方程：[0025]ΔPi＝Pi-UiIi＝PG-Pload-UiIi；[0026]式中，ΔPi为第i个节点的功率失配量，Ui表示节点i的电压，Ii表示节点i的注入电流；其中Ii的表达式为：
[0027][0028]
[0023]
式中，Yii表示节点i自导纳，Yij表示节点i与其相连节点互导纳，n表示直流配电网络中节点数目；由于电解槽工作电压低，通过直流换流器进行电压转换，UEYii为该直流换流器的接入对制氢系统接入节点注入电流的修正量；[0029]步骤S22：求取功率失配量对节点电压的一阶导数：
[0030]
[0031][0032]
其中，
7
CN 110912140 A[0033]
说 明 书
3/10页
[0034][0035]
进一步地，为了降低算法的复杂度，步骤S2还包括步骤S23：令压缩机消耗的有功
功率Pcom,p为恒定值，其对节点电压的导数为0，令节点发电机注入功率PG同样为恒定值，其对节点电压的导数为0，因此仅需考虑电解槽，将步骤S22中的一阶导数公式转换为：
[0036]
[0037][0038][0039]
其中，
[0040]
[0041][0042][0043][0044][0045][0046]
其中，电解槽注入电流为：
则，
[0047][0048]
式中，Icom为压缩机负荷注入节点电流，Icom＝-Pcom,p/Ui。[0049]进一步地，步骤S3具体包括以下步骤：[0050]步骤S31：形成直流配电网的雅可比矩阵，设制氢系统负荷接入直流配电网的节点i，配电网中含有n个节点，当不考虑负荷模型时，直流配电网络节点电压方程为：
[0051]
[0052]
Pi＝UiIi；
8
CN 110912140 A[0053][0054][0055][0056]
说 明 书
4/10页
建立此时的雅可比矩阵元素为：
当考虑负荷模型后，建立的雅可比矩阵为：
[0057]
[0058]
步骤S32：采用牛顿拉夫逊法进行迭代计算，第k+1次迭代计算后，节点i电压新值
为：
[0059][0060]
式中，为节点i第k次迭代计算得到的节点电压值，为第k+1次迭代计算得到
的电压修正量；
[0061]根据节点电压方程求出更新后的节点注入电流：
[0062][0063]
其中，Yij为节点i与其相连节点互导纳，根据公式得到电解槽注入电
流：
[0064][0065]
步骤S33：根据电解槽功率模型，更新功率值，若迭代计算没有收敛，则返回步骤
S31再次求取雅可比矩阵，进入下一次迭代，若迭代计算收敛，则进入步骤S34；[0066]步骤S34：得到当前的电解槽注入电流，并根据该电解槽注入电流求取制氢系统的功率值；
[0067]步骤S35：进一步得到制氢系统的功率值Pload：
[0068]
较佳的，本发明还可以根据上文的计算结果进一步得到电解槽的工作状态值，其包括电解池内每个电解槽端电压，电解槽内电流密度，以及电解池制氢效率。
[0070]本发明还提供了一种计及电解水制氢系统精细化模型的直流配电网潮流计算系统，包括存储器、处理器以及存储于存储器上并能够被所述处理器运行的计算机程序，当处理器运行该计算机程序时，实现如上文所述的方法步骤。[0071]与现有技术相比，本发明有以下有益效果：本发明在进行潮流计算时抛弃传统的将制氢系统看作电流恒定或功率恒定的负荷的方法，考虑负荷模型，使制氢系统的内部电
9
[0069]
CN 110912140 A
说 明 书
5/10页
流参与到迭代计算中，及时修正节点电压量，与配电网运行相互协调，能够使得得到的潮流计算结果更加准确，即所求取的制氢系统功率值与电解槽的工作状态值更加准确。附图说明
[0072]图1为本发明实施例的方法流程示意图。
[0073]图2为本发明实施例的两端直流配电网结构示意图。
具体实施方式
[0074]下面结合附图及实施例对本发明做进一步说明。[0075]应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
[0076]需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。[0077]如图1所示，本实施例提供了一种计及电解水制氢系统精细化模型的直流配电网潮流计算方法，包括以下步骤：[0078]步骤S1：分别根据压缩机与电解槽的内部电压、电流关系得到压缩机功率模型与电解槽功率模型；[0079]步骤S2：求取制氢系统节点注入功率，并以此确定功率失配量方程，求取功率失配量对节点电压的一阶导数；[0080]步骤S3：采用牛顿法，建立直流配电网的雅可比矩阵，进行迭代计算，当迭代计算收敛后得到电解槽的功率，进而得到制氢系统功率值。[0081]在本实施例中，步骤S1具体为：
[0082]根据电解槽的内部电流电压关系得到所述电解槽功率模型：
[0083][0084][0085]
其中，
[0086][0087][0088]
式中，为电解槽功率，U0为驱动电解水反应的最小电压，为电解槽注入电流，
R为普适气体常数，T为电解槽工作温度，z为半反应电荷转移化学计量系数，F为法拉第常
10
CN 110912140 A
说 明 书
6/10页
数，αA为电解槽横截面积，βtmem为膜的厚a为阳极电荷转移系数，1的计算公式如上式所示，度，σjoa为阳极交换电流密度，jL为极限电流密度，为氢气分压；mem为膜的电导，
[0089][0090][0091]
根据压缩机的进出口压强得到所述压缩机功率模型：
式中，Pcom,p为压缩机消耗的有功功率，ngas表示压缩机出口的气体流速，P1表示压缩机入口压强，P2表示压缩机出口压强，m表示压缩机运行在多变过程下的指数。[0092]在本实施例中，步骤S2具体包括以下步骤：[0093]步骤S21：求取制氢系统的节点注入功率：
[0094]
式中，Pi为第i个节点的注入功率，PG表示节点发电机注入功率，Pload表示节点负荷功率，N为电解槽数目；
[0096]根据求取的制氢系统节点注入功率建立失配量方程：[0097]ΔPi＝Pi-UiIi＝PG-Pload-UiIi；[0098]式中，ΔPi为第i个节点的功率失配量，Ui表示节点i的电压，Ii表示节点i的注入电流；其中Ii的表达式为：
[0099][0100]
[0095]
式中，Yii表示节点i自导纳，Yij表示节点i与其相连节点互导纳，n表示直流配电网络中节点数目；由于电解槽工作电压低，通过直流换流器进行电压转换，UEYii为该直流换流器的接入对制氢系统接入节点注入电流的修正量；[0101]步骤S22：求取功率失配量对节点电压的一阶导数：
[0102]
[0103][0104][0105]
其中，
[0106][0107]
在本实施例中，为了降低算法的复杂度，步骤S2还包括步骤S23：令压缩机消耗的有功功率Pcom,p为恒定值，其对节点电压的导数为0，令节点发电机注入功率PG同样为恒定值，其对节点电压的导数为0，因此仅需考虑电解槽，将步骤S22中的一阶导数公式转换为：
11
CN 110912140 A
说 明 书
7/10页
[0108]
[0109][0110][0111]
其中，
[0112]
[0113][0114][0115][0116][0117][0118]
其中，电解槽注入电流为：
则，
[0119][0120]
式中，Icom为压缩机负荷注入节点电流，Icom＝-Pcom,p/Ui。
[0121]在本实施例中，步骤S3具体包括以下步骤：[0122]步骤S31：形成直流配电网的雅可比矩阵，设制氢系统负荷接入直流配电网的节点i，配电网中含有n个节点，当不考虑负荷模型时，直流配电网络节点电压方程为：
[0123]
[0124][0125][0126][0127][0128]
Pi＝UiIi；
建立此时的雅可比矩阵元素为：
当考虑负荷模型后，建立的雅可比矩阵为：
12
CN 110912140 A
说 明 书
8/10页
[0129]
[0130]
步骤S32：采用牛顿拉夫逊法进行迭代计算，第k+1次迭代计算后，节点i电压新值
为：
[0131][0132]
式中，为节点i第k次迭代计算得到的节点电压值，为节点i第k+1次迭代计
算得到的电压修正量；
[0133]根据节点电压方程求出更新后的节点注入电流：
[0134][0135]
其中，Yij为节点i与其相连节点互导纳，根据公式得到电解槽注入电
流：
[0136][0137]
步骤S33：根据电解槽功率模型，更新功率值，若迭代计算没有收敛，则返回步骤
S31再次求取雅可比矩阵，进入下一次迭代，若迭代计算收敛，则进入步骤S34；[0138]步骤S34：得到当前的电解槽注入电流，并根据该电解槽注入电流求取制氢系统的功率值；
[0139]步骤S35：进一步得到制氢系统的功率值Pload：
[0140]
较佳的，本实施例还可以根据上文的计算结果进一步得到电解槽的工作状态值，其包括电解池内每个电解槽端电压，电解槽内电流密度，以及电解池制氢效率。
[0142]本实施例还提供了一种计及电解水制氢系统精细化模型的直流配电网潮流计算系统，包括存储器、处理器以及存储于存储器上并能够被所述处理器运行的计算机程序，当处理器运行该计算机程序时，实现如上文所述的方法步骤。[0143]特别的，本实施例采用如图2所示的两端直流配电网结构，电压等级为10kV。程序运行环境为Matlab2017b，CPU为Inter(R)Core(TM)i5-3235M，2.6GHz。光伏发电PDG1＝0.7MW,储能供电PDG2＝0.7MW。VSC1采用定Udc、定Qs控制，VSC1中变压器阻抗标幺值Zt＝0.015+0.112j，滤波器导纳标幺值Bf＝0.0877，电抗器以及换流器损耗等效阻抗标幺值Zc＝0.0001+0.6428。DC/DC1与DC/DC2变换器分别位于直流节点1与4处，皆采用不调压控制方式，占空比D＝0.285，IPOS拓扑结构中变换器子模块数量N＝6。变换器位于直流节点2处，采用不调压控制方式，拓扑结构中，连接直流负荷制氢系统，制氢系统中电解槽与压缩机参数分别见表1、表2。高压侧电压基准值UBH＝10kV，低压侧电压基准值UBL＝0.38kV，功率基准值SB＝100MVA。交直流子系统收敛精度皆为0.0001。
13
[0141]
CN 110912140 A[0144]
说 明 书
9/10页
表1电解槽参数值
[0145]
[0146]
表2压缩机参数取值
[0147]
[0148]
如图2所示，分布式电源以及直流负载均通过控制直流变换器的变比连接到低压
侧。网络潮流计算结果影响制氢系统内部电流，进而影响功率，计算负荷功率与节点功率差值，根据功率失配量对节点电压进行修正。虽然考虑制氢系统内部电压电流关系减缓了收敛速度，但是负荷功率在迭代过程中可以随着节点电压的变化而变化，便于观察配网中电压以及功率波动下制氢系统的运行状况，对于发展非并网风电理论，解决无法上网的风电及光伏发电利用率问题具有重要意义。
[0149]本实施例的直流潮流计算结果以及制氢系统的潮流计算结果分别如表3与表4所示。
[0150]表3直流潮流计算结果
14
CN 110912140 A
说 明 书
10/10页
[0151]
[0152]
表4制氢系统潮流计算结果
[0153]
[0154]
其中，电解池由多个电解槽串联，每个电解槽电压1.778V，电解池内部电流
有名值为1.99KA，消耗功率0.745MW，制氢系统与压缩机消耗功率总和为
0.0077，有名值0.77MW。本实施例的计算结果皆在合理范围内。
[0155]以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。
15
CN 110912140 A
说 明 书 附 图
1/1页
图1
图2
16
因篇幅问题不能全部显示，请点此查看更多更全内容