电-氢混合储能型多微电网系统有功功率均衡控制方法

李奇; 李蕊睿; 李朔; 蒲雨辰; 孙彩; 陈维荣

doi:10.3969/j.issn.0258-2724.20210506

电-氢混合储能型多微电网系统有功功率均衡控制方法

doi: 10.3969/j.issn.0258-2724.20210506

西南交通大学电气工程学院，四川成都 610031

基金项目: 国家自然科学基金（51977181，52077180）；霍英东教育基金会高等院校青年教师基金（171104）

详细信息

作者简介:
李奇（1984—），男，教授，博士生导师，研究方向为轨道交通新能源技术、综合能源系统运行与控制等，E-mail：liqi0800@163.com

中图分类号: TM732
计量
- 文章访问数: 281
- HTML全文浏览量: 170
- PDF下载量: 65
- 被引次数: 0
出版历程
- 收稿日期: 2021-06-22
- 修回日期: 2022-01-22
- 网络出版日期: 2024-03-13
- 刊出日期: 2022-03-31

Control Method for Active Power in Electric-Hydrogen Hybrid Energy-Storage Microgrids

School of Electrical Engineering, Southwest Jiaotong University, Chengdu 610031, China

摘要

摘要:
在多电-氢微电网并联运行的交流系统中，采用传统控制方法的逆变器会受到输出线路阻抗差异所产生的影响产生较大环流，甚至无法实现有功功率的合理分配，考虑到系统中电压偏差与有功功率间的关系，对多微电网的有功功率分配问题进行分析研究. 首先，构建含光伏、燃料电池、电解槽、蓄电池的电-氢混合储能微电网交流系统模型；其次，根据反下垂控制中有功功率与电压间的关系，构造考虑电压偏差的反下垂控制，并为使额定有功功率自适应调节，提出基于功率跟随控制的反下垂控制方法；最后，在多电-氢微电网并联运行的系统中，对本文所提方法进行RT-LAB半实物实验验证及与其他方法的对比验证. 实验结果显示：基于本文所提控制方法的系统稳定后功率分配的准确度达97.50%，母线电压的精准度达99.86%，环流大小的范围约为[−3.0, 3.0] A，均优于其余方法.
- 电-氢微电网 /
- 混合储能 /
- 功率控制
Abstract:
In an AC system with parallel operation of electricity-hydrogen micro-grids, when adopting traditional control methods, the difference in output line impedance will affect inverters, and large circulating currents occurs such that a reasonable distribution of active power become inaccessible. Given the relationship between the voltage deviation and the active power, the active power distribution between multiple micro-grids is analyzed. Firstly, the AC system model of the electric-hydrogen hybrid energy-storage micro-grid is constructed, including photovoltaic, battery, fuel cell, and electrolyzer. Secondly, following the relationship between the active power and the voltage in the reverse droop control, a reverse droop control based on power following control is proposed to work out voltage deviation and auto-adjust the rated active power. Finally, within a parallel operation of electric-hydrogen micro-grids, the proposed method is compared with other methods and is verified by hardware-in-the-loop experiments on the platform of RT-LAB. The experimental results show that the proposed method outperforms other methods: the accuracy of power distribution after stabilization is 97.50%, the accuracy of the bus voltage is 99.86%, and the circulating current mostly ranges in [−3.0, 3.0] A.
- electric-hydrogen microgrid /
- hybrid energy storage /
- power control

HTML全文

双制式列车能实现市域铁路与市内城轨交通的互联互通、贯通运营，减少换乘. 然而，双制式列车存在车体—钢轨回流、车体电位抬升和轴承电蚀等问题，车体回流过大、车体电位过高会对车载弱电设备造成影响，并恶化列车的电磁环境，而轴承的损伤直接影响列车的行车安全. 这些问题与列车整车接地系统的设计密切相关，因此，研究双制式列车接地系统配置具有重要意义.

目前，针对列车接地系统的研究主要集中在单一制式列车方面，如直流（DC）制或交流（AC）制列车. 文献[1]量化了直流制地铁列车轴承电流的大小及其与牵引回流的关系，并提出一种抑制轴承电流的接地方式. 文献[2]建立了动态工况下的地铁列车车网联合仿真模型，对列车接地回流分布进行仿真分析. 文献[3-5]建立高速动车组接地系统的仿真模型，通过优化列车接地系统，降低车体环流和车体—轴端电位.

对于双制式牵引供电系统，为限制不同供电制式系统之间发生牵引回流交换，建议在双制式系统转换段使用绝缘节与特殊馈电装置和开关装置连接^[6]. 针对过渡区段设置方案，文献[7-8]研究了双制式牵引供电系统复合钢轨电位问题，提出应在交、直转换段钢轨处设置绝缘节；文献[9-10]分析了当双制式线路的直、交流区段由同一主所供电时，交流区段产生的谐波问题会影响直流区段；文献[11]建立了双制式列车通过锚段关节式电分相过程的等值电路模型，并分析了列车通过锚段关节式电分相时的暂态特性.

双制式列车接地系统需兼容交、直牵引回流. 目前，对于双制式列车接地系统的研究较少. 文献[12]介绍了双制式列车牵引传动系统. 文献[13]提出一种双制式列车接地系统方案设计，提出在各节车厢车体与汇流排之间均应设置保护电阻.

目前，对列车接地系统的研究主要以仿真建模为主，缺乏统一的数学模型，这限制了列车接地系统在各区段运行过程性能的深入研究. 为此，本文提出列车接地系统的车—轨链式等效电路模型，并结合牵引供电系统进行车-地一体化联合计算，分析双制式列车运行在直流区段、交流区段的车体环流分布，并对列车接地系统提出改进建议. 该方法也适用于单供电制式或多供电制式列车接地系统的分析.

1. 双制式列车接地系统

双制式列车可在交、直流区段贯通运行. 交、直流区段之间需通过交—直转换段进行分离，其结构如图1所示.

图 1 双制式牵引供电系统

Figure 1. Dual-system traction power supply system

下载: 全尺寸图片幻灯片

为适应直流和交流两种供电制式，双制式列车除了和普通地铁列车一样装有牵引逆变器以外，也装有与牵引变压器和牵引变流器. 列车在交—直转换段完成不同制式的切换.

国内某型双制式列车为6节编组，其接地系统如图2所示. 列车牵引传动和接地系统均首尾对称. 在直流牵引工况下，列车工作接地与保护接地共用同一接地汇流排. 在交流牵引工况下，列车工作接地设置专用汇流排与保护接地的汇流排分开.

图 2 线路列车既有接地系统

Figure 2. Existing train grounding system

下载: 全尺寸图片幻灯片

1.1 车—轨链式电路

列车接地系统可以分为车体、接地汇流排和钢轨三层结构，如图3所示. 图中：Z_ctlj和Z_ct分别为列车车间连接线阻抗和列车车体阻抗，Z_hlp为汇流排间连接线阻抗，R₀为车体与接地汇流排之间的保护电阻，Z_jd1和Z_jd2为接地汇流排与轮对之间连接线阻抗，Z_r1为同一转向架轮对间的钢轨阻抗，Z_r2为不同一转向架间的钢轨阻抗，Z_r3为不同车厢间的钢轨阻抗. 接地汇流排切面包含4个节点（图3中 ①～④）.

图 3 列车接地系统电路模型

Figure 3. Circuit model of train grounding system

下载: 全尺寸图片幻灯片

单个接地汇流排的节点导纳矩阵 ${{\boldsymbol{Y}}_{{\text{cb}}}}$ 如式（1）所示.

$\qquad\qquad\qquad\quad {{\boldsymbol{Y}}_{{\text{cb}}}} = \left( {\begin{array}{*{20}{c}} {1/{R_0}}&{ - 1/{R_0}}&0&0 \\ { - 1/{R_0}}&{1/{R_0} + 1/{Z_{{\text{jd}}2}} + 1/{Z_{{\text{jd}}1}}}&{ - 1/{Z_{{\text{jd}}2}}}&{ - 1/{Z_{{\text{jd}}1}}} \\ 0&{ - 1/{Z_{{\text{jd}}1}}}&{1/{Z_{{\text{jd}}1}} + 1/{Z_{{\text{r}}1}}}&{ - 1/{Z_{{\text{r}}1}}} \\ 0&{ - 1/{Z_{{\text{jd}}2}}}&{ - 1/{Z_{{\text{r}}1}}}&{1/{Z_{{\text{jd}}2}} + 1/{Z_{{\text{r}}1}}} \end{array}} \right).$

(1)

各接地汇流排之间还存在车间连接线、汇流排间连接线等横连线，其中，车间连接线的横连线导纳矩阵Y_lj如式（2）所示.

${{\boldsymbol{Y}}_{{\text{lj}}}} = \left( {\begin{array}{*{20}{c}} {1/{Z_{{\text{ctlj}}}}}&{}&{ - 1/{Z_{{\text{ctlj}}}}}&{} \\ {}&{\boldsymbol{O}}&{}&{\boldsymbol{O}} \\ { - 1/{Z_{{\text{ctlj}}}}}&{}&{1/{Z_{{\text{ctlj}}}}}&{} \\ {}&{\boldsymbol{O}}&{}&{\boldsymbol{O}} \end{array}} \right),$

(2)

式中：O为3阶零矩阵.

由于列车接地汇流排的布置具有对称性，可以根据列车接地汇流排划分切面，形成车—轨链式等效电路，如图4所示. 图中，Y_cbi为第i个接地汇流排节点导纳矩阵，Z_hli为横连线阻抗矩阵，I_cti为注入第i个接地汇流排的牵引电流向量， $i = 1,2, \cdots ,m$ .

图 4 车—轨链式电路

Figure 4. Train–rail chain circuit

下载: 全尺寸图片幻灯片

形成的列车车—轨导纳矩阵为

${{\boldsymbol{Y}}_{{\text{CT}}}} = \left( {\begin{array}{*{20}{c}} {{\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb1}}}}}&{ - {\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1}}&{}&{}&{}&{} \\ { - {\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1}}&{{\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb2}}}} + {\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1}}&{ - {\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1}}&{}&{}&{} \\ {}&{ - {\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1}}&{{\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb3}}}} + {\boldsymbol{Z}}_{{\text{hl3}}}^{ - 1}}&{ - {\boldsymbol{Z}}_{{\text{hl3}}}^{ - 1}}&{}&{} \\ {}&{}& & \vdots & &{} \\ {}&{}&{}&{ - {\boldsymbol{Z}}_{{\text{hl}}m{{-2}}}^{ - 1}}&{{\boldsymbol{Z}}_{{\text{hl}}m{{-2}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb}}m{{-}}1}} + {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}}&{ - {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}} \\ {}&{}&{}&{}&{ - {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}}&{{{\boldsymbol{Y}}_{{\text{cb}}m}} + {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}} \end{array}} \right).$

(3)

对于牵引电流向量，需分别考虑列车处于直流和交流两种工况. 含有牵引电流注入的车—轨切面电流向量如式（4）所示.

$\begin{gathered} {{\boldsymbol{I}}_{{\text{ct}}i}}(t) = \left\{ {\begin{array}{*{20}{c}} {{{\left[ {\begin{array}{*{20}{c}} 0&{{{\dot I}_{{\text{tract}}i}}(t)}&0&0 \end{array}} \right]}^{\mathrm{T}}}\begin{array}{*{20}{c}} , \;\; &{{T_{\mathrm{r}}}(t) = {T_{{\mathrm{r}}0}}}, \end{array}} \\ {{{\left[ {\begin{array}{*{20}{c}} 0&0&{{{\dot I}_{{\text{tract}}i}}(t)}&0 \end{array}} \right]}^{\mathrm{T}}}\begin{array}{*{20}{c}} , \;\; &{{T_{\mathrm{r}}}(t) = {T_{{\mathrm{r}}1}}}, \end{array}} \end{array}} \right. \end{gathered}$

(4)

式中：T_r(t)为时刻t列车的供电制式状态；T_r0为列车处于直流制牵引供电；T_r1为列车处于交流制牵引供电； $\dot{I}_{{\mathrm{tract}}i}(t)$ 为时刻t注入第i个接地汇流排的牵引电流，直流制式下 $\dot{I}_{{\mathrm{tract}}i}(t)$ 仅取实部.

对列车车—轨首、末切面的电流向量进行修正，修正后的车—轨首、末切面的电流向量I_ct1new、I_ctmnew如式（5）所示.

$\begin{gathered} \left\{ {\begin{array}{*{20}{c}} {{{\boldsymbol{I}}_{{\text{ct1new}}}} = {{\boldsymbol{I}}_{{\text{ct1}}}} + {{\left[ {\begin{array}{*{20}{c}} 0&0&{ - {{\dot I}_{{\text{GL}}}}(t)}&0 \end{array}} \right]}^{\mathrm{T}}}}, \\ {{{\boldsymbol{I}}_{{\text{ct}}{{m}}{\text{new}}}} = {{\boldsymbol{I}}_{{\text{ct}}{{m}}}} + {{\left[ {\begin{array}{*{20}{c}} 0&0&0&{ - {{\dot I}_{{\text{GR}}}}(t)} \end{array}} \right]}^{\mathrm{T}}}}, \end{array}} \right. \\ \end{gathered}$

(5)

$\sum\limits_{i = 1}^m {{{\dot I}_{{\text{tract}}i}}(t) = {{\dot I}_{{\text{tr}}}}(t) = {{\dot I}_{{\text{GL}}}}(t) + {{\dot I}_{{\text{GR}}}}(t)},$

(6)

式中： $\dot{I}_{\mathrm{tr}}(t)$ 为列车取流， $\dot{I}_{\mathrm{GL}}(t)$ 为车—轨首切面钢轨电流， $\dot{I}_{\mathrm{GR}}(t)$ 为车—轨末切面钢轨电流. 直流制式下 $\dot{I}_{\mathrm{GL}}(t)$ 、 $\dot{I}_{\mathrm{GR}}(t)$ 仅取实部.

列车车—轨电流向量为

${{\boldsymbol{I}}_{{\text{CT}}}} = {({{\boldsymbol{I}}_{{\text{ct1new}}}},{{\boldsymbol{I}}_{{\text{ct2}}}},{{\boldsymbol{I}}_{{\text{ct3}}}}, \cdots ,{{\boldsymbol{I}}_{{\text{ct}}(m - 2)}},{{\boldsymbol{I}}_{{\text{ct}}(m - 1)}},{{\boldsymbol{I}}_{{\text{ct}}m{\text{new}}}})^{\mathrm{T}}}.$

(7)

形成的车—轨链式等效电路节点电压向量为

${{\boldsymbol{U}}_{{\text{CT}}}} = {\boldsymbol{Y}}_{{\text{CT}}}^{ - 1}{{\boldsymbol{I}}_{{\text{CT}}}}.$

(8)

2. 双制式牵引供电系统建模

2.1 交、直流区段牵引供电系统建模

交、直流区段的牵引网均满足平行多导体结构，故可采用链式电路模型表示. 交、直流区段l长度的牵引网阻抗矩阵 ${{\boldsymbol{Z}}_l}$ 、导纳矩阵 ${{\boldsymbol{Y}}_l}$ 如式（9）所示^[14].

$\left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{Z}}_l} = \sinh (\sqrt {{\boldsymbol{ZY}}} l){{({\boldsymbol{ZY}})}^{ - \frac{1}{2}}}{\boldsymbol{Z}} =\\ \quad {\boldsymbol{Z}}{\displaystyle\sum\limits_{n = 0}^\infty {\dfrac{{{l^{2n + 1}}}}{{(2n + 1)!}}{{({\boldsymbol{ZY}})}^n}} }, \\ {{\boldsymbol{Y}}_l}/2 = {\boldsymbol{Z}}{{({\boldsymbol{ZY}})}^{ - \frac{1}{2}}}\tanh (\sqrt {{\boldsymbol{ZY}}} l/2) = \\ \quad {{\boldsymbol{Z}}_l}^{ - 1}\displaystyle\sum\limits_{n = 1}^\infty {\dfrac{{{l^{2n}}}}{{(2n)!}}{{({\boldsymbol{ZY}})}^n}}, \end{array}} \right.$

(9)

式中：Z、Y分别为牵引网单位阻抗矩阵、导纳矩阵.

在直流区段中，只考虑Z、Y的电阻部分.

牵引变电所可采用诺顿等效电路模型. 交流牵引变电所α、β供电臂的等效电流源 $\dot{I}_{{\mathrm{equ}}\alpha}$ 、 $\dot{I}_{{\mathrm{equ}}\beta }$ 可用式（10）求得^[15].

$\left\{ {\begin{array}{*{20}{l}} {\left[ {\begin{array}{*{20}{c}} {{{\dot{I}}_{{\text{equ}}\alpha }}}&{} \\ {}&{{{\dot{I}}_{{\text{equ}}\beta }}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{Z_\alpha }} \\ {{Z_\beta }} \end{array}} \right] = {\boldsymbol{N}}{{\boldsymbol{E}}_{{\mathrm{ABC}}}} - {Z_\gamma }\left[ {\begin{array}{*{20}{c}} {{{\dot{I}}_\beta }} \\ {{{\dot{I}}_\alpha }} \end{array}} \right], } \\ {{{\boldsymbol{Z}}_{\alpha \beta }} = \left[ {\begin{array}{*{20}{c}} {{ Z_\alpha }}&{{ Z_\gamma }} \\ {{ Z_\gamma }}&{{ Z_\beta }} \end{array}} \right],} \end{array}} \right.$

(10)

式中：N为电压变化矩阵，E_ABC为牵引变压器一次侧电压向量，Z_α、Z_β、 $Z_\gamma$ 为变电所T型等效电路等值阻抗，Z_αβ为变电所等效阻抗矩阵， $\dot{I}_\alpha$ 、 $\dot{I}_\beta$ 分别为变电所α、β供电臂端口电流.

直流牵引变电所使用24脉波整流机组（24 MPR），其等效电流源I_d和端口等效电阻r₀为^[16]

$\left\{ {\begin{array}{*{20}{l}} {{I_{\text{d}}} = \dfrac{{{U_{\text{d}}}}}{{{r_0}}} = \dfrac{1}{{{r_0}}}\dfrac{{24\sqrt 2 }}{{\text{π}} }\sin \dfrac{{\text{π}} }{{24}}{U_2} = 1.41\dfrac{{{U_2}}}{{{r_0}}}}, \\ {{r_0} = \dfrac{3}{{2{\text{π}} }}{X_{\text{c}}}}, \end{array}} \right.$

(11)

式中：U_d为24 MPR等效电压源，U₂为24 MPR二次侧线电压额定值，X_c为24 MPR换流阻抗.

2.2 交、直流区段联合牵引供电计算

在双制式线路上存在不同供电制式的列车运行，故需对交、直流区段进行联合牵引供电计算.

对于牵引网的建模，若线路均为高架区段，交流区段采用直供带回流方式，其牵引网由回流线、接触网、钢轨、综合地线构成，单位牵引网阻抗、导纳为4阶矩阵；直流区段牵引网由接触网、钢轨、排流网构成，单位牵引网阻抗、导纳为3阶矩阵. 为形成统一的牵引网阻抗、导纳矩阵，直流区段需添加虚拟导线使其单位牵引网阻抗、导纳变为4阶矩阵. 为使得系统节点导纳矩阵非奇异，虚拟导线与钢轨之间可设置大电阻横向连接. 双制式线路牵引网模型如所示. 图中， $\dot{I}_{\mathrm{trdc}}$ 、 $\dot{I}_{\mathrm{trac}}$ 分别为列车在直流区段、交流区段取流.

图 5 双制式牵引网模型

Figure 5. Dual-system traction network model

下载: 全尺寸图片幻灯片

根据叠加定理，交、直流区段统一牵引供电计算中节点电压方程组为

$\left\{ {\begin{array}{*{20}{c}} {{{\boldsymbol{U}}_{{\text{Sdc}}}} = {\boldsymbol{G}}_{{\text{S}}j(a + b) \times j(a + b)}^{ - 1}{{\boldsymbol{I}}_{{\mathrm{dcs}}}}}, \\ {{{\boldsymbol{U}}_{{\text{Sac}}}} = {\boldsymbol{Y}}_{{\text{S}}j(a + b) \times j(a + b)}^{ - 1}{{\boldsymbol{I}}_{{\mathrm{acs}}}}}, \end{array}} \right.$

(12)

式中：U_Sdc为系统节点电压向量直流分量，U_Sac为系统节点电压向量交流分量，I_dcs为直流电流向量，I_acs为交流电流向量，G_S为仅直流源作用时系统牵引网节点导纳矩阵，Y_S为仅交流源作用时系统牵引网节点导纳矩阵，j为单位牵引网阻抗导纳矩阵阶数，a、b分别为直流区段、交流区段切面数.

3. 车-地一体的交、直流区段联合牵引供电计算

由双制式牵引系统供电计算可以得到各时刻列车的取流信息以及列车位置相邻切面的钢轨电流信息. 根据这些数据，结合车—轨链式等效电路模型，可以得出列车在各时刻的车体环流分布以及各节车厢车体—轴端电位.

车-地一体的交、直流区段联合牵引供电计算电路模型如图6所示，图6（a）中列车接地系统等值电路见图6（b）. 图中：Z_c为接触网阻抗，Z_r为钢轨阻抗，Z_p为排流网阻抗，Z_fb为回流线阻抗， Z_r-p为钢轨对排流网阻抗，Z_p-e为排流网对地阻抗，Z_gw为贯通地线阻抗，Z_e为贯通地线接地阻抗.

图 6 车—地一体的交、直流区段联合牵引供电计算电路模型

Figure 6. Circuit model of united traction power supply calculation for vehicle-ground integration in AC and DC sections

下载: 全尺寸图片幻灯片

该联合计算流程如图7所示. 该方法不仅适用于双流制接地系统的研究，也可用于单一供电制式或者其他多供电制式列车接地系统的分析.

图 7 供电系统与列车接地系统联合计算流程

Figure 7. Flowchart of united calculation of power supply system and train grounding system

下载: 全尺寸图片幻灯片

4. 列车接地保护电阻优化模型

4.1 目标函数

从车体—轴端电位角度，建立列车接地保护电阻优化模型，尽可能地抑制车体—轴端电位大小.

选择列车运行过程中，列车各车车体—轴端电位最大值作为目标函数f₁(X)，如式（13）所示.

$\left\{ {\begin{array}{*{20}{c}} {\min {f_1}(X) = {U_{{\text{CZ,}}k}}}, \\ {X = \{ {x_k}\} ,\quad {x_k} \in \varepsilon } , \end{array}} \right.$

(13)

式中：U_CZ,k为列车第k节车厢车体—轴端电位最大值，X为列车各车保护电阻值配置集合，x_k为列车第k节车厢的保护电阻值，ε为以10 mΩ为步长的随机离散变量组合.

4.2 约束条件

求解单目标函数式（13），其约束条件为列车保护电阻上、下限值以及通过保护电阻的车体电流上限值，如式（14）所示.

$\left\{ {\begin{array}{*{20}{c}} {{x_{\min }} \leqslant {x_k} \leqslant {x_{\max }}}, \\ {I_{{\text{CT}},k} \leqslant {I_{{\text{CT}},\max }}}, \end{array}} \right.$

(14)

式中：x_max、x_min分别为保护电阻设置的上、下限，I_CT,k为列车第k节车厢保护电阻的电流最大值，I_CT,max为列车保护电阻的电流上限值.

5. 算例分析

为验证算法的有效性，以国内某条双制式轨道交通线路（图8）为例，分析列车处于直流区段和交流区段2种情况下的车体环流分布，并与实测结果进行对比.

图 8 双制式轨道交通线路

Figure 8. Dual-system rail transit line

下载: 全尺寸图片幻灯片

5.1 仿真条件

该线路全长28.22 km，各牵混所以及主所位置见表1. 全线共设置2个直流车站和5个交流车站，所有车站均为高架站. T1和T2为区间牵引所，TS1为带牵引所的直流区段车站，S2～S7为普通车站，K为线路设计起点，TPS1和TPS2为交流区段牵引主变电所. 交—直转换段NS设置在S2和S3站点之间. K—NS为直流区段，采用1500 V柔性接触网授流，通过钢轨回流；采用地面制动电阻吸收列车再生制动能量；TS1处的钢轨电位限制器（OVPD）接地. NS—S7为交流区段，采用25 kV柔性接触网授流，通过回流线、贯通地线以及钢轨回流. 定义K—S7为上行方向， S7—K为下行方向.

表 1 牵引所以及主所位置

Table 1. Location of traction station and main substation km

牵混所编号	位置	主所编号	位置
T1	0.15	TPS1	19.81
TS1	1.74	TPS2	28.21
T2	3.70

下载: 导出CSV

| 显示表格

该型列车车辆参数见表2所示，其接地系统如图2所示. 定义其中1车为头车，方向朝S7；6车为尾车，方向朝TS1. 全天采用K—S7单一运行路线，发车间隔设置为500 s.

表 2 车辆参数

Table 2. Train parameters

参数	数值
列车编组	6A
车重/t	223.88
结构速度/（km·h⁻¹）	120
最大加速度/（m·s⁻²）	1.1 （DC），0.9 （AC）
最大减速度/（m·s⁻²）	1.2

下载: 导出CSV

| 显示表格

列车接地系统中各连接线阻抗以及车体阻抗由现场实际测量列车得出，对应参数见表3.

表 3 列车接地系统模型参数

Table 3. Model parameters of train grounding system

参数	AC	DC
Z_ct/mΩ	1.500	0.027
Z_ctlj/mΩ	0.82	0.82
R₀/mΩ	50.6	50.0
Z_jd1/mΩ	0.89	0.89
Z_jd2/mΩ	1.21	1.21
Z_hlp/mΩ	3.9	3.9
Z_r/（mΩ·km⁻¹）	191.0^[17]	36.4^[18]

下载: 导出CSV

| 显示表格

牵引计算模拟得到的上下行列车功率-时间曲线如图9所示.

图 9 上、下行列车功率曲线

Figure 9. Power curves of up and down trains

下载: 全尺寸图片幻灯片

5.2 实测验证

对该线路某辆列车的一组动力单元的接地回流系统的电流进行同步检测，在列车保护电阻处加装电流传感器，规定流出车体电流方向为正. 监测点的详细布置方案如图10所示，图中，k-p为第k节车厢第p个接地汇流排处的保护电阻，余图同. 不同位置设备的检测信号通过GPS同步授时.

图 10 监测点布置方案

Figure 10. Monitoring point layout

下载: 全尺寸图片幻灯片

5.2.1 直流区段

在直流工况下，列车牵引电流通过各车接地汇流排分散注入至钢轨. 选取列车在上行方向TS1—K区间运行过程，计算并分析列车车体环流分布情况. 该区间列车取流以及左、右切面钢轨电流计算值如图11所示. 列车一组动力单元的车体环流分布计算值和实测值如图12所示.

图 11 列车取流以及钢轨电流（直流）

Figure 11. Train’s current and rail current (DC)

下载: 全尺寸图片幻灯片

在该直流区间中，列车各保护电阻的电流计算值如图13所示. 由图可知，通过列车保护电阻电流的最大值为22.3 A. 该时刻列车处于牵引工况，车体环流路径为4～6车流入车体，1～3车流出车体.

图 12 保护电阻电流计算与实测对比（直流）

Figure 12. Comparison of protection resistance current calculation and measured value (DC)

下载: 全尺寸图片幻灯片

图 13 列车各保护电阻电流（直流）

Figure 13. Protection resistance current of train (DC)

下载: 全尺寸图片幻灯片

5.2.2 交流区段

列车在交流工况下，牵引电流通过3、4车牵引变压器一次侧专用汇流排注入至钢轨. 选取列车在上行方向S5—S4区间运行过程计算列车取流以及左、右切面钢轨电流计算值，如图14所示. 列车一组动力单元的车体环流分布计算值和实测值如图15所示.

图 14 列车取流以及钢轨电流（交流）

Figure 14. Train’s current and rail current (AC)

下载: 全尺寸图片幻灯片

在该交流区间中，列车各保护电阻的电流计算值如图16所示. 通过列车保护电阻电流的最大峰值为9.6A，该时刻列车处于牵引工况，车体环流路径为3～6车流入车体，1～2车流出车体.

图 16 列车各保护电阻电流（交流）

Figure 16. Protection resistance current of train (AC)

下载: 全尺寸图片幻灯片

计算与实测存在误差的原因是由于列车车体环流分布受列车司机操纵设置以及钢轨均流线等横联线位置的影响. 列车1～3车监测位置的保护电阻电流实测值与计算值的总体变化趋势大致相同，说明该车-地一体联合牵引供电计算模型可用于分析列车车体环流分布，验证了列车在交、直流区段时模型的准确性.

图 15 保护电阻电流计算与实测对比（交流）

Figure 15. Comparison of protection resistance current calculation and measured value (AC)

下载: 全尺寸图片幻灯片

5.3 对列车接地系统配置的讨论与分析

在正常工况下，列车接地系统配置需能抑制车体环流以及车体—轴端电位. 对于交流制动车组，车体—轴端电位的最大值不应超过1 V^[19]. 地铁列车的车体—轴端电位和车体环流没有明确的限值. 关于列车保护电阻，TB/T 2977—2016^[20]规定，列车车体与轨道之间的保护电阻值不应大于50 mΩ^[20].

由图12、15可见，在交、直流区段中该型双制式列车的车体环流分布总体上是头尾车（1、6车）最大，中间车（3、4车）最小. 为使列车在交、直流区段运行过程中降低头尾车车体—轴端电位的同时，尽量抑制车体环流，本文对既有的列车接地系统方案提出2种改进方案：方案1为更改1、6车的保护电阻值，方案2为在1、6车增设1处接地汇流排，并更改1、6车的保护电阻值. 改进方案如图17所示.

图 17 列车接地系统改进方案

Figure 17. Improved scheme of train grounding system

下载: 全尺寸图片幻灯片

设置x_min=10 mΩ，x_max=50 mΩ，I_CT,max=60 A，使用2种改进方案重新计算，可得其在交、直流区段运行中的车体电流和车体—轴端电位分布.

方案1的车体—轴端电位以及车体环流分布最大值如图18所示. 当1、6车的保护电阻从50 mΩ降为10 mΩ，列车在直流区段的车体—轴端电位最大值从1.25 V降至1.04 V，而车体电流最大值从22.33 A增至88.01 A. 列车在交流区段的车体—轴端电位最大值从0.89 V降至0.67 V，而车体电流最大值从16.86 A增至45.98 A.

图 18 列车的车体—轴端电位和车体环流分布（方案1）

Figure 18. Distribution of current circulation and potential of train’s body–axle end (scheme 1)

下载: 全尺寸图片幻灯片

方案2的车体—轴端电位以及车体环流分布最大值如图19所示. 在1、6车各增设1处汇流排，1、6车的保护电阻从50 mΩ降为10 mΩ，列车在直流区段的车体—轴端电位最大值从2.12 V降至1.64 V，而车体电流最大值从28.99 A增至107.97 A. 列车在交流区段的车体—轴端电位最大值从0.79 V降至0.50 V，而车体电流最大值从10.95 A增至25.11 A.

图 19 列车的车体—轴端分布和车体环流分布（方案2）

Figure 19. Distribution of current circulation and potential of train’s body–axle end (scheme 2)

下载: 全尺寸图片幻灯片

方案1相较方案2，列车的车体—轴端电位最大值降低36.58%～41.04%，车体环流最大值降低18.49%～22.97%；并且在方案1中，1、6车采用20 mΩ保护电阻接地可使列车的车体—轴端电位最大值为1.15 V，车体电流最大值为50.30 A，达到抑制车体—轴端电位的最优效果.

6. 结　论

1）建立了车地一体化的交、直流区段联合牵引供电计算模型，并将模型计算值与实测数据进行比较，验证模型的正确性.

2）列车车体环流除受自身牵引取流影响外，也受邻车位置和工况、牵引网横联线设置的影响. 在交、直流区段中，列车头尾车的车体电流和车体—轴端电位值最大，中间车的车体环流和车体—轴端电位值最小.

3）建立了列车接地保护电阻优化模型，并以既有的双制式列车接地系统为例，讨论分析了2种改进的接地系统方案. 方案1中头尾车采用20 mΩ保护电阻接地，能达到抑制车体—轴端电位的最优效果.

本文仅对稳态情况下列车接地保护电阻的设置进行了探讨，车体浪涌的暂态模型还需进行进一步研究.

图 1 多微电网并联运行结构示意

Figure 1. Schematic of micro-grid parallel operation

下载: 全尺寸图片幻灯片

图 2 多微网并联运行等效电路

Figure 2. Equivalent circuit of micro-grid parallel operation

下载: 全尺寸图片幻灯片

图 3 系统运行示意

Figure 3. Schematic of system operation

下载: 全尺寸图片幻灯片

图 4 功率跟随控制框图

Figure 4. Block diagram of power following control

下载: 全尺寸图片幻灯片

图 5 基于功率跟随控制的反下垂控制框图

Figure 5. Block diagram of reverse droop control based on power following control

下载: 全尺寸图片幻灯片

图 6 基于功率跟随控制的反下垂控制频谱响应特性

Figure 6. Spectrum response characteristics of reverse droop control based on power following control

下载: 全尺寸图片幻灯片

图 7 RT-LAB实时仿真平台

Figure 7. Real-time simulation platform RT-LAB

下载: 全尺寸图片幻灯片

图 8 传统下垂控制方法实验结果

Figure 8. Experimental results of traditional droop control

下载: 全尺寸图片幻灯片

图 9 基于电压恢复机制的反下垂控制方法实验结果

Figure 9. Experimental results of reverse droop control based on voltage recovery mechanism

下载: 全尺寸图片幻灯片

图 10 RDVD控制方法实验结果

Figure 10. Experimental results of RDVD control method

下载: 全尺寸图片幻灯片

图 11 本文所提控制方法实验结果

Figure 11. Experimental results of proposed control method

下载: 全尺寸图片幻灯片

表 1 系统参数

Table 1. System parameters

参数	取值
逆变器 1 线路阻抗 Z_{line_1}/Ω	0.642 + j0.083
逆变器 2 线路阻抗 Z_{line_2}/Ω	0.963 + j0.124 5
逆变器 3 线路阻抗 Z_{line_3}/Ω	1.284 + j0.166
滤波电感/MH	5
滤波电容/MF	0.5
额定电压/V	220
额定无功功率/kvar	5
逆变器容量比例	1∶1∶1

下载: 导出CSV

表 2 元件参数

Table 2. Component parameters

元件	参数	数值
光伏阵列	环境温度/℃	25
燃料电池	额定功率/W	2000
燃料电池	额定电压/V	24
电解槽	额定功率/W	2000
储氢罐	最大容许压强/MPa	35
	体积/L	1
	初始容量/%	50
蓄电池	最大充放电功率/W	4000
	容量/（A·h）	40
	额定电压/V	500
	初始荷电状态/%	50

下载: 导出CSV

参考文献(16)

[1]	李霞林,郭力,王成山,等. 直流微电网关键技术研究综述[J]. 中国电机工程学报,2016,36(1): 2-17. LI Xialin, GUO Li, WANG Chengshan, et al. Key technologies of DC microgrids: an overview[J]. Proceedings of the CSEE, 2016, 36(1): 2-17.
[2]	LI Q, QIU Y, YANG H, et al. Stability-constrained two-stage robust optimization for integrated hydrogen hybrid energy system[J]. CSEE Journal of Power and Energy Systems, 2021, 7(1): 162-171.
[3]	刘欣博,刘宁,宋晓通,等. 基于交流恒功率负载特性的交直流混合微电网系统大信号稳定性判据[J]. 高电压技术,2021,47(10): 3441-3451. LIU Xinbo, LIU Ning, SONG Xiaotong, et al. Large-signal stability criteria of AC/DC hybrid microgrid based on AC constant power loads[J]. High Voltage Engineering, 2021, 47(10): 3441-3451.
[4]	蔡国伟,孔令国,薛宇,等. 风氢耦合发电技术研究综述[J]. 电力系统自动化,2014,38(21): 127-135. doi: 10.7500/AEPS20131231004 CAI Guowei, KONG Lingguo, XUE Yu, et al. Overview of research on wind power coupled with hydrogen production technology[J]. Automation of Electric Power Systems, 2014, 38(21): 127-135. doi: 10.7500/AEPS20131231004
[5]	EHSAN A, YANG Q. Optimal integration and planning of renewable distributed generation in the power distribution networks: a review of analytical techniques[J]. Applied Energy, 2018, 210: 44-59. doi: 10.1016/j.apenergy.2017.10.106
[6]	李奇,蒲雨辰,韩莹,等. 电-氢孤岛直流微电网的分层能量管理[J]. 西南交通大学学报,2020,55(5): 912-919. LI Qi, PU Yuchen, HAN Ying, et al. Hierarchical energy management for electric-hydrogen island direct current micro-grid[J]. Journal of Southwest Jiaotong University, 2020, 55(5): 912-919.
[7]	米阳,陈鑫,季亮,等. 基于虚拟额定电流的直流微电网分布式储能单元精确电流分配研究[J]. 电网技术,2020,44(3): 823-835. MI Yang, CHEN Xin, JI Liang, et al. Accurate current sharing of distributed energy storage units in DC microgrid based on virtual rated current[J]. Power System Technology, 2020, 44(3): 823-835.
[8]	朱晓荣,韩丹慧,孟凡奇,等. 提高直流微电网稳定性的并网换流器串联虚拟阻抗方法[J]. 电网技术,2019,43(12): 4523-4531. ZHU Xiaorong, HAN Danhui, MENG Fanqi, et al. Grid converter series virtual impedance method for improving DC microgrid stability[J]. Power System Technology, 2019, 43(12): 4523-4531.
[9]	柴秀慧,张纯江,柴建国,等. 改进互联通信荷电状态下垂控制及功率均衡优化[J]. 电工技术学报,2021,36(16): 3365-3374. CHAI Xiuhui, ZHANG Chunjiang, CHAI Jianguo, et al. Improved interconnected communication state of charge droop control and power balance optimization[J]. Transactions of China Electrotechnical Society, 2021, 36(16): 3365-3374.
[10]	王晓寰,王书光,刘聪哲,等. 下垂控制并联系统中功率精确分配的控制策略[J]. 电力电子技术,2019,53(8): 4-7. WANG Xiaohuan, WANG Shuguang, LIU Congzhe, et al. Accurate power sharing strategy for parallel system based on droop control[J]. Power Electronics, 2019, 53(8): 4-7.
[11]	麦倩屏,陈鸣. 用于多微源低压微电网的虚拟阻抗反下垂控制[J]. 电力系统保护与控制,2018,46(1): 96-102. doi: 10.7667/PSPC162039 MAI Qianping, CHEN Ming. P-V/Q-f droop control strategy with virtual impedance for low-voltage microgrid with multiple micro sources[J]. Power System Protection and Control, 2018, 46(1): 96-102. doi: 10.7667/PSPC162039
[12]	白小丹,苗虹,曾成碧,等. 适用于低压微网中逆变器无功均分的改进下垂控制策略[J]. 高电压技术,2020,46(4): 1310-1318. BAI Xiaodan, MIAO Hong, ZENG Chengbi, et al. Improved droop control strategy for reactive power sharing of inverters in low-voltage microgrids[J]. High Voltage Engineering, 2020, 46(4): 1310-1318.
[13]	PU Y C, LI Q, CHEN W R, et al. Hierarchical energy management control for islanding DC microgrid with electric-hydrogen hybrid storage system[J]. International Journal of Hydrogen Energy, 2019, 44(11): 5153-5161. doi: 10.1016/j.ijhydene.2018.10.043
[14]	LI L Y, CHEN W R, HAN Y, et al. A stability enhancement method based on adaptive virtual resistor for electric-hydrogen hybrid DC microgrid grid-connected inverter under weak grid[J]. Electric Power Systems Research, 2021, 191: 106882.1-106882.
[15]	LI Q, WANG T H, DAI C H, et al. Power management strategy based on adaptive droop control for a fuel cell-battery-supercapacitor hybrid tramway[J]. IEEE Transactions on Vehicular Technology, 2018, 67(7): 5658-5670. doi: 10.1109/TVT.2017.2715178
[16]	IPSAKIS D, VOUTETAKIS S, SEFERLIS P, et al. Power management strategies for a stand-alone power system using renewable energy sources and hydrogen storage[J]. International Journal of Hydrogen Energy, 2009, 34(16): 7081-7095. doi: 10.1016/j.ijhydene.2008.06.051