Evaluation Method of Railway Schemes Along Rivers in Risk Areas of Moraine-Dammed Lake Outburst
-
摘要:
为建立青藏高原南缘冰碛湖溃决危险区铁路方案风险评价方法，首先，对青藏高原南缘74例冰碛湖溃决灾害实例溃决机制进行统计分析；在此基础上，针对冰碛湖溃决影响因素众多、关系复杂且具有高度不确定性的特点，利用多态贝叶斯网络建立了冰碛湖溃决概率预测模型；然后，以冰碛湖溃决危险区沿河线路工程为承灾体，综合冰碛湖溃决危险性评估与沿河线河谷地貌特征提出了各类风险区线路长度计算程式，建立了针对冰碛湖溃决危险性的铁路选线方案评价方法；最后，以中尼铁路跨喜马拉雅段的樟木、吉隆局部走向方案为例，说明评价方法的作业程式. 研究结果表明：樟木方案在冰碛湖分布区的线路长度小于吉隆方案，从选线角度属于线路短直方案，但樟木方案各风险区段线路总长超出吉隆方案约45%；从冰碛湖溃决风险角度，吉隆方案优于樟木方案.
Abstract:This study is aimed at risk assessments of railway schemes in moraine-dammed lake outburst risk areas in the southern Tibet Plateau. Based on the mechanism analysis of 74 outburst events in the study area, a multi-state Bayesian-Network prediction model is established considering the complex relationships among the influencing factors and high uncertainties of the outburst events. Taking the railway line engineering along rivers in moraine-dammed lake outburst risk areas as the hazard bearing body and combining the risk assessment of moraine-dammed lakes with geomorphological characteristics of the river valley, a calculation program of the railway scheme length in various risk areas and an evaluation system for railway schemes considering moraine-dammed lake outburst risk are proposed. In addition, the operation procedure of the proposed method is demonstrated taking the Zhangmu and Gyirong schemes of the China-Nepal railway as examples The results show that the Zhangmu scheme has a smaller route length in the moraine lake distribution area than the Gyirong scheme, and belongs to the short straight scheme from the perspective of route selection; however, the total length of the Zhangmu scheme in each risk area exceeds that of the Gyirong scheme by about 45%. Therefore, the Gyirong scheme is the better selection when considering the moraine lake outburst risk.
-
Key words:
- moraine-dammedlake outburst /
- B-Ns /
- China–Nepal railway /
- risk assessment
-
2020年，中国在第75届联合国大会上提出，力争2030年前二氧化碳排放达到峰值，努力争取2060年前实现碳中和[1]. 在双碳目标和全球能源革命的背景下，新能源产业将迎来高质量、跨越式发展. “多能互补、源网荷储协同”的综合能源发展规划已经明确，清洁低碳、安全高效的能源体系正在加快构建[2]. 光伏以其取之不尽、用之不竭的优点受到青睐.
近年来，我国铁路交通发展迅猛，路网运输能力和效率显著提升，其能耗问题也日益显现[3]. 光伏接入铁路系统作为一项新型的铁路节能技术，在国内外的铁路车站、信号系统以及站舍屋顶等非牵引领域已大量应用[4-5]，而在牵引供电领域的研究及应用相对较少. 考虑到牵引能耗占铁路系统总能耗比例最大的现实情况，将光伏发电技术应用在牵引供电领域是降低铁路系统能耗、实现铁路低碳化运行的关键[6]. 能够使铁路沿线土地及太阳光资源得到充分利用，可有效减少光伏发电传输、转换过程中的能量损耗，响应国家“清洁电能就近消纳”的号召[7]，显著缓解外部电网的供电压力. 然而，在把光伏接入牵引供电系统的过程中存在诸多困难，主要有光伏电能消纳、电能质量影响、列车再生制动能量返送等[8]. 文献[9]介绍了德国铁路牵引供电系统中2种将光伏发电直接接入牵引供电系统的应用形式，但是，此方法无法解决光伏接入对电能质量的影响问题.
铁路功率调节器（RPC）是一种能够实现牵引供电系统无功、负序、谐波综合补偿的电力电子装置. RPC通过控制牵引供电系统中两供电臂的能量双向流动，实现动态补偿，有效对牵引供电系统电能质量进行改善[10-11]. 考虑光伏出力的特性以及列车负荷的多元化和不确定性，在RPC 的直流环节引入储能部分，及时对功率进行超额吸收和缺额补偿，平抑光伏功率波动. 由于电网供电是单向的，其自身不具有储存能量的功能，因此，引入储能部分能够实现对列车再生制动能量的回收再利用. 文献[12]提出了一种基于RPC的光伏并网功率调节器及其控制策略，在实现光伏有效并网的同时实现了对牵引供电系统无功、负序和谐波的综合补偿，但是缺少储能环节，难以在实际情况中加以应用. 文献[13]为有效回收电气化铁路交流传动型机车产生的大量再生制动电能，同时兼顾改善牵引供电系统电能质量的目的，基于RPC的结构，提出一种背靠背混合储能系统及其控制方法. 文献[14]同样针对高速铁路牵引供电系统中再生制动能量不能得到有效利用的问题，提出一种基于RPC的牵引供电系统储能方案及其控制方法，有效回收利用高速铁路牵引供电系统产生的再生能量，但是未考虑光伏接入系统的情况.
上述文献对RPC的光伏接入、电能质量补偿、再生制动回收等方面的应用进行了研究，但是均聚焦在单一方面，未能针对RPC的特性进行充分利用.
综上，本文利用RPC作为光伏并网接口，同时接入储能部分，提出一种基于RPC的光储接入牵引供电系统的协调控制方法，根据RPC补偿方案拓扑结构给出各类典型运行工况能量协调控制策略，进行仿真分析验证所提出方法的可行性、有效性.
1. 基于RPC的光储系统接入牵引供电系统
1.1 系统结构
基于RPC的光储系统接入牵引供电系统结构如图1所示，包括光储发电系统、RPC以及牵引供电系统3部分，图中：$ {\dot I_{\rm{A}}}、{\dot I_{\rm{B}}}、{\dot I_{\rm{C}}}$分别为电力系统的三相电流，${\dot I_\alpha }、{\dot I_\beta }$分别为两供电臂（$\alpha、\beta $相）等效牵引电流，${\dot I_{\alpha {\rm{L}}}}、{\dot I_{\beta {\rm{L}}}}$分别为两供电臂负荷电流，${\dot I_{{\rm{c}}\alpha} }、{\dot I_{{\rm{c}}\beta} }$分别为两供电臂补偿电流，${\dot U_\alpha }、{\dot U_\beta }$分别为牵引变压器副边电压， ${L_\alpha }$和${L_\beta }$分别为两供电臂电感. 在光储发电系统中，光伏阵列和储能系统分别通过各自的DC/DC变换器接入RPC的直流电容两侧，为牵引供电系统提供电能. 本文使用蓄电池作为储能部分. 其中，光伏阵列作为光储发电系统的主要电能来源，储能系统作为辅助，在平抑光伏功率波动的同时，还能对机车的再生制动能量进行回收.
图 1 基于RPC的牵引供电系统综合补偿结构图Figure 1. Comprehensive compensation structure of traction power supply system based on RPCRPC作为光伏接入牵引供电系统的接口并且负担电能补偿功能，考虑应用稳定性，采用基于全桥逆变结构，其输出端口基于电感${L_\alpha }$和${L_\beta }$实现L型滤波，如图2所示. RPC包括2个背靠背结构的电压源变流器和一个直流电容${C_{\rm{d}}}$，连接在两供电臂中间，实现两供电臂之间有功功率的转移，同时也能对无功和谐波进行补偿，从而达到无功、负序、谐波综合补偿的多种目标.
牵引供电系统中，来自电力系统的220 kV三相工频交流电通过V/V降压变压器转换成27.5 kV的单相工频交流电，为α相和β相两供电臂上运行的电力机车供电.
1.2 电能质量补偿机理
为了实现RPC对牵引供电系统电能质量问题的动态补偿，本文采用瞬时鉴相检测法来检测负序电流的参考信号. 假设牵引供电系统原边的电网A相单位电压$ {u_{\rm{A}}}\left( t \right) = \sin {\omega t} $，那么2个牵引供电臂的单位电压分别表示为$ {e_\alpha }\left( t \right) = \sin \left( {\omega t - {\text{π}}/6} \right)$和$ {e_\beta }\left( t \right) = \sin \left( {\omega t - {\text{π}} /2} \right)$. 用傅里叶级数表示2个供电臂的负载电流$ i_{\alpha {\rm{L}}}$和$ i_{\beta {\rm{L}}}$分别为[15]
{iαL=Iαpsinωt+Iαqcosωt+∞∑n=2Iαnsin(nωt+φαn),iβL=Iβpsinωt+Iβqcosωt+∞∑n=2Iβnsin(nωt+φβn), (1) 式中：${I_{\alpha {\rm{p}}}}$和${I_{\beta {\rm{p}}}}$分别为$\alpha、\beta $相两供电臂负载电流的有功分量幅值，${I_{\alpha {\rm{q}}}}$和${I_{\beta {\rm{q}}}}$分别为两供电臂负载电流的无功分量的幅值，${I_{\alpha n}}{\rm{sin}}\left( {n\omega t + {\varphi _{\alpha n}}} \right)$为α相的n次谐波电流， $I_{\beta n} \sin(n\omega t + \varphi_\beta n )$为β相n次谐波电流，${\varphi _{\alpha n}}$和${\varphi _{\beta n}}$分别为左、右供电臂n次谐波电流相位.
根据鉴相检测原理[16]，先将两牵引供电臂负载电流分别乘以同步信号，相加之后通过低通滤波器（lower pass filter，LPF）可以提取出基波有功电流幅值的1/2，提取后的电流${I_{\rm{p}}}$表示为
Ip=12(Iαp+Iβp). (2) 根据负序补偿原理，计算得到补偿后的两臂期望电流分量为
{iα(t)=Ipsin(ωt−π6)+Iptanπ6cos(ωt−π6),iβ(t)=Ipsin(ωt−π2)+Iptanπ6cos(ωt−π2). (3) 为实现RPC的补偿功能，需在两臂负荷电流的基础上注入补偿电流，使牵引供电系统中两臂等效电流达到式（3）表示的期望值，故RPC的两侧补偿电流指令值${i_{{\text{c}}\alpha }}$和${i_{{\text{c}}\beta }}$可以表示为
{icα(t)=iα(t)−iαL(t),icβ(t)=iβ(t)−iβL(t). (4) 式（4）中包含了谐波电流分量，因此，不但能够补偿负序，而且能够起到抑制谐波的作用.
当光储接入后，此时流过变流器的电流不仅有补偿电流，还包括光储提供给牵引供电系统的有功分量. 为实现光储的有效并网的应用目标，又要保留RPC对无功、负序、谐波的补偿功能，所以需要改变变流器的期望指令. 光伏储能系统接入的目的是给牵引供电系统提供有功功率，换言之就是将光储所输出的有功功率输送给两供电臂进行利用，即在${i_\alpha } $和$ {i_\beta }$上叠加相应的光储有功电流分量. 叠加后得到两相电流的期望值$ {i_{{\text{d}}\alpha }} $和$ {i_{{\text{d}}\beta }} $分别为
{idα=iα+Iαhsin(ωt+π),idβ=iβ+Iβhsin(ωt+π3), (5) Iαh=Iβh=√2(Ppv+Pbat)Uα+Uβ, (6) 式中：${I_{\alpha {\rm{h}}}}$、${I_{\beta {\rm{h}}}}$为光储输出的有功正分量幅值，${P_{{\rm{pv}}}}$为光伏的输出功率，${P_{{\rm{bat}}}}$为蓄电池的输出或吸收的功率，${U_\alpha } 、 {U_\beta }$分别为牵引变压器副边两相电压有效值.
此时，只需将$ {i_{{\rm{d}}\alpha }} $、$ {i_{{\rm{d}}\beta }} $与$ {i_{\alpha {\rm{L}}}} $、$ {i_{\beta {\rm{L}}}}$相减可得到新的参考信号.
2. 协调控制方法
2.1 铁路功率调节器控制策略
在不考虑光储接入的情况下，若要保证RPC的正常工作，此时需要控制RPC的变流器使两供电臂进行功率交换，同时维持直流环节的电压稳定. 直流侧的误差通过PI控制器后，分别与α相和β相电压的同步信号$ {\rm{sin}}\left( {\omega t - {{30}^ \circ }} \right)$和$ {\rm{sin}}\left( {\omega t - {{90}^ \circ }} \right)$相乘，得到RPC两变流器的有功指令，然后，再跟负序和谐波参考信号相叠加就得到了最后的实际指令$i_{\alpha {\rm{c}}}^*$、$i_{\beta {\rm{c}}}^*$. 两变流器采用PI控制[10]方法，使其具有较快的响应速度.
当光储接入后，为使并网后的电流实现平衡，需要对RPC原有的控制策略进行改进. 改进部分的控制策略如图3的α相和β相变流器控制所示，由于加入了光伏和蓄电池储能，故需要将这两部分所提供的功率平均地分配到RPC两侧. 图中：${P_{{\rm{c}}{\rm{max}}}}$和${P_{{\rm{d}}{\rm{max}}}}$为蓄电池的最大充、放电功率，前者取负值；${P_{{\rm{lim}}}}$为光伏的功率限值；$ {u_\alpha }、{u_\beta }$分别为两供电臂的电压实测值；$P_{{\rm{ES}}} $为光储部分为两臂提供的总功率； $U_{{\rm{dc}}}^*$为直流侧电压给定值；$ {U_{{\text{dc}}}} $为直流侧电压测量值，两者比较计算得到直流侧的误差通过PI调节后，分别与电力系统A相和B相电压的同步信号相乘，叠加到两相电流期望值上，再跟${i_{\alpha {\rm{L}}}}$和${i_{\beta {\rm{L}}}}$相减，得到新指令，其他部分控制方法不变. 此时，负序和谐波依旧得到有效的补偿.
2.2 光储荷系统控制策略
光伏和储能接入牵引供电系统是未来绿色铁路发展的一个重要方向，光伏发电能够给牵引供电系统负荷供电，实现铁路沿线光伏电能的就近消纳，同时储能系统能够灵活储存和释放电量，容纳列车的再生制动能量，防止其返送电能对牵引供电系统电能质量的负面影响，更好地实现节能减排. 所以需要良好的协调控制策略来控制光伏和储能的能量流动，使光储系统实时调整电能分配，能够更好地满足牵引所需用能.
光储荷系统的功率控制流程如图4所示，图中：蓄电池的荷电状态（SOC）上、下限分别为${\rm{SO}}{{\rm{C}}_{{\rm{max}}}}$和${\rm{SO}}{{\rm{C}}_{{\rm{min}}}}$， ${P_\alpha }、{P_\beta }$为两臂功率. 首先，需要获得两供电臂的电压电流实际测量值，计算出两臂的功率 ${P_\alpha }、 {P_\beta }$. 考虑到RPC能够控制两供电臂间功率流动，此处可以考虑采用两供电臂的功率和$P = {P_\alpha } + {P_\beta }$来判断，大于0、小于0、等于0分别表示此时牵引供电系统局部处于牵引、制动或空载工况，以此简化工况类型并且便于分析. 对蓄电池采取限值管理策略防止过充过放损害使用寿命，以其荷电状态SOC为依据，监测并判断蓄电池是否正常进行充放电工作.
图 4 光储荷系统功率控制流程图Figure 4. Flow chart of power control for “PV + battery + load” system工况1：局部牵引工况下，当列车牵引需求量大于光伏发电量，即$P \geqslant {P_{{\rm{PV}}}}$，若蓄电池SOC已到达下限，蓄电池无法提供电能，光储系统中仅有光伏提供牵引用能，即${P_{{\rm{ES}}}} = {P_{{\rm{PV}}}}$.
工况2：局部牵引工况下，若$P \geqslant {P_{{\rm{PV}}}}$且蓄电池${\rm{SOC}} \geqslant {\rm{SO}}{{\rm{C}}_{{\rm{min}}}}$，则光伏和蓄电池一起提供电能. 当$P \geqslant {P_{{\rm{PV}}}} + {P_{{\rm{d}}{\rm{max}}}}$,则蓄电池以最大放电功率进行放电尽可能满足牵引用能，此时${P_{{\rm{ES}}}} = {P_{{\rm{PV}}}} + {P_{{\rm{d}}{\rm{max}}}}$.
工况3：局部牵引工况下，若$P \geqslant {P_{{\rm{PV}}}}$且${\rm{SOC}} \geqslant {{\rm{SOC}}_{{\rm{min}}}}$，如果$P < {P_{{\rm{PV}}}} + {P_{{\rm{d}}{\rm{max}}}}$，则蓄电池放电功率为${P_{{\rm{bat}}}} = P - {P_{{\rm{PV}}}}$，${P_{{\rm{ES}}}} = P$.
以上3种工况中，蓄电池均面临释放自身电能的外界需求，此时仅需${{\rm{SOC}}_{{\rm{min}}}}$作为标准进行判断，无需考虑蓄电池过充问题.
工况4：局部牵引工况下，若牵引需求量小于光伏所发电能$P < {P_{{\rm{PV}}}}$，则光伏单独为列车供电且有余电. 若此时蓄电池SOC到达其上限值$（{\rm{SOC}} \leqslant {{\rm{SOC}}_{{\rm{max}}}}）$，此时蓄电池无法储存多余的光伏发电电量，此时只能对光伏的输出功率进行限制，即${P_{{\rm{PV}}}} = {P_{{\rm{lim}}}}$，${P_{{\rm{ES}}}} = {P_{{\rm{lim}}}}$.
工况5：局部牵引工况下，若$P < {P_{{\rm{PV}}}}$且蓄电池${\rm{SOC}} \leqslant {\rm{SOC}}{_{{\rm{max}}}}$，则光伏单独为列车供电，光伏发电余量输送到蓄电池进行存储. 若$P \leqslant{P_{{\rm{PV}}}} - \left| {{P_{{\rm{c}}{\rm{max}}}}} \right|$，则蓄电池以最大充电功率为目标进行充电.
工况6：局部牵引工况下，若$P < {P_{{\rm{PV}}}}$且蓄电池${\rm{SOC}} \leqslant {{\rm{SOC}}_{{\rm{max}}}}$，如果$P > {P_{{\rm{PV}}}} - \left| {{P_{{\rm{c}}{\rm{max}}}}} \right|$，则蓄电池充电，为${P_{{\rm{bat}}}} = P - {P_{{\rm{PV}}}}$，此时${P_{{\rm{ES}}}} = P$.
工况7：局部制动工况下，若此时蓄电池SOC到达其上限值，则蓄电池受到过充限制不能继续充电，此时不能回收再生制动能量及光伏发电的电能.
工况8：局部制动工况下，如果蓄电池${\rm{SOC}} \leqslant {\rm{SO}}{{\rm{C}}_{{\rm{max}}}}$，能够正常充电，并且$\left| P \right| + {P_{{\rm{PV}}}} \geqslant {P_{{\rm{c}}{\rm{max}}}}$，则蓄电池以最大充电功率充电以尽可能回收光伏及再生制动能量，此时${P_{{\rm{bat}}}} = {P_{{\rm{c}}{\rm{max}}}}$.
工况9：局部制动工况下，如果蓄电池${\rm{SOC}} \leqslant {\rm{SO}}{{\rm{C}}_{{\rm{max}}}}$，且$\left| P \right| + {P_{{\rm{PV}}}} <{P_{{\rm{c}}{\rm{max}}}}$，表明此时蓄电池容量充足，能够完全容纳光伏发电电量、回收列车再生制动能量，此时${P_{{\rm{bat}}}} = \left| P \right| + {P_{{\rm{PV}}}}$.
工况4～9中，蓄电池均面临储存电能的外界需求，此时仅需${\rm{SO}}{{\rm{C}}_{{\rm{max}}}}$作为标准进行判断，无需考虑蓄电池过放问题.
工况10：局部的空载工况情况，即无列车经过或列车惰行情况下，此时两供电臂上无功率需求也无再生制动能量反送，若${\rm{SO}}{{\rm{C}}_{{\rm{min}}}} \leqslant {\rm{SOC}} \leqslant {\rm{SO}}{{\rm{C}}_{{\rm{max}}}}$，表明此时蓄电池能够正常工作，可以对光伏发电电能进行吸收.
工况11：局部的空载工况情况，若蓄电池SOC已达充放电上下限，此时蓄电池不能对光伏发电电能进行回收.
3. 仿真分析
为验证所提协调控制方法的有效性，本文在RT-Lab实时仿真机上搭建了基于RPC的光储系统接入牵引供电系统的实时仿真模型，并分别在RPC的电能质量补偿效果、光储荷系统功率分配以及实际工况下系统运行仿真结果. 仿真主要参数如表1.
表 1 仿真模型主要参数Table 1. Main parameters of simulation model子系统 参数 数值 牵引变压器 变比 220/27.5 降压变压器 变比 27.5/1 变流器 直流侧电压/V 4000 RPC 直流电容/F 0.02 RPC 交流电感/H 0.0001 光伏装机 容量/MW 3.265 蓄电池 容量/MW 6 3.1 RPC电能质量补偿效果
为验证RPC系统对电能质量的补偿效果，设定2种局部牵引工况.
工况A1：左侧供电臂负载功率为5 MW，右侧供电臂负载功率为3 MW.
工况A2：左侧供电臂负载功率为3 MW，右侧供电臂空载.
2种工况下的补偿前、后三相侧电流如图5所示，在第1 s时接入RPC进行补偿. 从图5中可以看到：接入RPC补偿前，三相电流各相幅值差异明显；1.0 s时接入RPC补偿后，经过0.1 s左右的时间后实现了稳定，此时三相电流基本对称，负序有了明显的改善，验证了RPC的补偿功能和动态响应性能.
图 5 补偿前、后三相侧电流变化Figure 5. Three-phase side current changes before and after compensation3.2 光储荷系统功率流分配结果
考虑到铁路沿线光照强度不同、列车的制动牵引等工况不断变化，修改光照条件和供电臂运行工况实现对光伏发电和列车负荷的仿真，验证在光伏发电及列车负荷变化下，光储荷系统协调控制策略分配功率的显著效果，如图6所示. 图中：PL为左侧供电臂负载功率，PR为右侧供电臂负载功率.
工况B1 （0～1.0 s）：两臂均处于牵引工况，此时设定左臂牵引功率为3.000 MW，右臂牵引功率为2.000 MW. 在0～0.5 s时间段内，光伏强度较弱，此时仅能输出1.806 MW，为满足两供电臂上的列车负载需求，蓄电池输出功率约为3.200 MW. 在0.5～1.0 s时间段内，由于光照强度增强，光伏输出逐步增长为3.265 MW，此时蓄电池相应减少能量释放，此时输出功率约为1.730 MW. 整个过程中，两臂的负载所需电能均由光储系统提供.
工况B2 （1.0～2.0 s）：左侧供电臂制动右侧供电臂牵引，此时设定左臂再生制动功率为1.500 MW，右臂牵引功率为1.000 MW. 在1.0～1.5 s时间段内，光照强度较弱，光伏输出1.806 MW，因此，蓄电池吸收列车再生制动能量及光伏发电能量，蓄电池储能功率约为2.300 MW. 在1.5～2.0 s时间段内，光伏强度增大，光伏输出增长为3.265 MW，此时蓄电池吸收功率为3.760 MW. 整个过程中蓄电池保持正常充电状态.
工况B3 （2.0～3.0 s）：两臂均处于空载或惰行状态，此时光储系统与牵引供电系统没有能量的交换，光伏发电功率直接送到储能部分进行储存，2.5 s时光照强度突变，光伏发电量增加，蓄电池继续储能，整个过程中除去未稳定时刻，蓄电池吸收功率恒等于光伏发电功率.
工况B4 （3.0～4.0 s）：两臂均处于制动工况，此时设定左臂再生制动功率为0.500 MW，右臂再生制动功率为1.000 MW. 在3.0～3.5 s时间段内，光伏输出1.806 MW，蓄电池吸收光伏发电能量及两供电臂上的再生制动能量，功率为3.300 MW，在3.5～4.0 s时间段内，光伏输出突变3.265 MW，蓄电池储能功率相应增加约为4.770 MW. 整个过程中，两臂的再生制动能量与光伏所发电能均由蓄电池吸收.
综上所述，该控制策略能够在不同工况及光照条件下，实现光储系统按照牵引供电系统所需功率进行合理的能量分配与转移，同时保证了再生制动能量的有效回收.
3.3 实际工况下系统运行结果
为分析实际工况下光储和RPC接入后对牵引供电系统电能质量的影响，截取了某牵引变电所某时间段的左、右供电臂实际负荷情况并压缩成50.0 s，如图7所示.
图7中：负荷大于0表示处于牵引工况，小于0表示处于制动工况，等于0表示惰行或者空载状态. 由图7可以看到，牵引工况下最大负荷可以达到20 MW左右. 根据仿真情况并通过计算，50.0 s时间段内再生制动能量占总牵引能量的9.2%，当通过RPC接入光伏储能系统后，再生制动能量的回收率达到86.7%. 在较为复杂和多变的实际负荷情况下，依旧能够达到较高的回收率，有效实现制动能量回收，有助于实现能量的再利用，减少反送电能对牵引供电系统电能质量的负面影响，节省铁路运行成本.
将当地实际光照强度的两组数据输入到光伏仿真模型中得到光伏输出功率，如图8所示，图8（a）对应的光照强度变化相对剧烈，所以造成第一组光伏功率波动较大，而图8（b）对应的光照强度变化相对较小，故光伏输出功率较为平滑.
仿真得到光储 + RPC接入前三相电流的负序情况，如图9. 可以看到此时负序电流较大，波动明显.
图 9 光储 + RPC接入前负序电流Figure 9. Negative sequence current before PV and battery + RPC access图10为光储 + RPC接入后三相电流的负序情况. 通过比较接入前、后的负序电流波形图可以看到：接入后的负序电流相比接入前的负序电流幅值减小，接入后的负序电流保持在15 A以内，基本保持在5 A以下，接入前后最大降幅能够达到55 A左右，整体负序补偿效果十分显著；当左右两臂负荷差距较大时，负序电流较大，符合理论分析结果. 针对2组不同光照的数据进行仿真，在光照强度相差较大的情况下，2组数据的负序电流相差不大且较为相似，可见，光照强度等环境因素的变化对电流负序的影响不大. 采用本文的控制方法能较好地改善负序的同时对光伏发电电量进行消纳. 且补偿前后高压侧的功率因数也有了一定的改善，从原来的平均功率因数0.863提升到0.918，整体提高了0.055，符合国家标准.
图 10 光储 + RPC接入后2种光照下的负序电流Figure 10. Negative sequence current under two kinds of illumination after PV and battery + RPC access为了验证谐波的补偿效果，截取了某时间段内某牵引变电所供电臂的1～50次谐波，并将其注入到仿真模型中. 通过谐波分析可以发现，3、5、7、9、11、31、33、35、37、39次的谐波电流含有率（HRI）相对较高，如表2所示. 针对光储 + RPC接入前、后的2种情况，比较各次谐波电流含有率的变化情况. 仿真结果如图11.
表 2 各次谐波电流含有率Table 2. Current ratio of each harmonic谐波次数/次 含有率/% 谐波次数/次 含有率/% 3 9.725 31 5.108 5 2.683 33 3.230 7 2.750 35 9.415 9 1.750 37 3.388 11 1.575 39 4.708 
图 11 光储 + RPC接入前后谐波电流含有率变化Figure 11. Changes in current ratio of each harmonic before and after PV and battery + RPC access对比图11可以看出：光储 + RPC接入后谐波电流的含有率有了很大的改善，对于高次谐波的改善更加明显. 其中： A相电流总谐波畸变率从19.16%降低至5.66%；B相电流总谐波畸变率从15.29%降低至3.01%；C相电流总谐波畸变率从18.82%降低至4.69%. 验证了RPC的谐波抑制效果，表明该控制策略有利于提高系统稳定性及电能质量.
4. 结 论
本文提出了光储通过铁路功率调节器接入牵引供电系统的控制策略，并根据预设工况和实际工况验证了该系统及策略的有效性和优越性，可以得出以下结论：
1） 光储通过铁路功率调节器接入牵引供电系统，通过合理的功率控制方式能够提高光伏的利用率，减小牵引供电系统对电力系统的依赖，并且能够回收再生制动能量，有利于提高能源利用效率.
2） 通过建立系统模型，基于预设工况下RPC接入前后三相电流由不对称变为基本对称. 在两供电臂牵引、制动、空载工况下，各部分功率流都进行了合理有效的转移和利用，符合预期期望.
3） 在实际工况及不同的光照情况下，接入后负序电流都有了极大的改善，功率因数也有了一定的提高，兼顾了较高的再生制动能量回收率，实现较好的节能效果. 本文所采用的数据时间较短，考虑到实际应用情况，还需要收集更多数据进行仿真验证对模型进行优化.
4） 基于本文的仿真结果进行分析，就改善的效果而言，电能质量各项指标都有了一定程度的提高，但就经济性而言，如何配置光储的容量以及蓄电池的种类参数选取等仍需继续研究，也是今后研究的重点.
-
图 1 青藏高原南缘冰川、冰湖及冰碛湖溃决灾害分布
Figure 1. Distributions of glaciers, glacier lakes and moraine-dammed lake outburst events in the southern Tibet Plateau
图 3 基于B-Ns的冰湖溃决危险性评估程式
Figure 3. Risk assessment procedure for moraine-dammed lake outburst based on B-Ns method
图 4 各类风险区线路长度计算流程
Figure 4. Calculation process of the railway line length in different risk areas
图 5 吉隆、波曲流域危险冰碛湖与线路方案示意
Figure 5. Railway schemes and distribution of dangerous moraine-dammed lakes in the Gyirong and Poiqu river
表 1 青藏高原南缘冰碛湖主要溃决机制统计
Table 1. Main outburst mechanisms of the moraine-dammed lakes in the southern Tibet Plateau
冰碛湖主要溃决机制 案例数 占比/% （冰崩/冰滑坡）涌浪漫顶 46 74.2 （降雨/冰川融水）洪水漫顶 12 19.4 管涌 5 8.1 冰碛物坍塌 5 8.1 外界动力激发 1 1.6 
下载: 导出CSV
表 2 根节点因素多区间状态划分标准
Table 2. Division of risk factors in multi-interval safety states for root nodes
二级指标变量 二级指标状态值 0 1 x1/（°） < 8 > 8 x2 < 0.11 > 0.11 x3/km > 1 < 1 x4 相对不发育 相对较发育 x5 < 0.23 > 0.23 x6/m > 60 < 60 x7 固结 松散 x8 未满水 满水 x9 < 0.20 > 0.20 x10/（°） < 20 > 20 x11 无/少量 丰富 x12 < 0.25g > 0.25g
下载: 导出CSV
表 3 冰碛湖溃决根节点的先验概率
Table 3. Prior probabilities of root nodes of moraine-dammed lake outburst
节点 状态数量 先验概率 状态 1 状态 2 x1 2 0.524 0.476 x2 0.286 0.714 x3 0.625 0.375 x4 0.192 0.808 x5 0.632 0.368 x6 0.240 0.760 x7 0.823 0.177 x8 0.423 0.577 x9 0.577 0.423 x10 0.712 0.288 x11 0.481 0.519 x12 0.900 0.100
下载: 导出CSV
表 4 中间节点y1的条件概率
Table 4. Conditional probabilities of node y1
x1 x2 x3 x4 x5 y1 0 1 2 0 0 0 0 0 1.0 0 0 0 0 0 0 1 1.0 0 0 0 0 0 1 0 1.0 0 0 0 0 0 1 1 0.2 0.8 0 0 0 1 0 0 1.0 0 0 0 0 1 0 1 0.2 0.8 0 $\vdots $ $\vdots $ $\vdots $ $\vdots $ $\vdots $ $\vdots $ $\vdots $ $\vdots $ 1 1 1 0 1 0 0.2 0.8 1 1 1 1 0 0 0.2 0.8 1 1 1 1 1 0 0 1.0
下载: 导出CSV
表 5 铁路方案冰碛湖溃决危险性评价指标
Table 5. Evaluation indexes of railway schemes considering moraine-dammed lake outburst risk
指标名称 比选因素 流域指标 冰碛湖总数量/个 冰碛湖总面积/km2 历史冰碛湖溃决案例数/次 线路指标 冰碛湖分布区线路总长度/km 溃决高风险区段线路长度/km 溃决次高风险区段线路长度/km 溃决中风险区段线路长度/km
下载: 导出CSV
表 6 吉隆、樟木线路方案比选表
Table 6. Comparison of the railway schemes in the Gyirong and Poiqu river basins
指标类型 比选因素 吉隆 樟木 流域指标 冰碛湖总数量/个 130 121 冰碛湖总面积/ km2 5.8 17.3 历史冰碛湖溃决次数/次 2 6 线路指标 冰碛湖分布区
线路长度/ km107.9 93.5 冰碛湖溃决高危险区段
线路长度/ km21.2 38.3 冰碛湖溃决次高危险区
段线路长度/ km54.1 132.0 冰碛湖溃决中危险区段
线路长度/ km10.5 54.1
下载: 导出CSV
-
[1] 王欣, 刘时银, 丁永建. 中国喜马拉雅山冰碛湖溃决灾害评价方法与应用研究[M]. 北京: 科学出版社, 2016. [2] ZHANG C, HUANG Y D, YAO L K, et al. Prediction of the maximum wave elevation in moraine-dammed lakes during resonant earthquake excitation[J]. Journal of Engineering Mechanics, 2020, 146(2): 04019118.1-04019118.13. [3] 徐道明,冯清华. 西藏喜马拉雅山区危险冰湖及其溃决特征[J]. 地理学报,1989,44(3): 343-351,385. doi: 10.3321/j.issn:0375-5444.1989.03.011XU Daoming, FENG Qinghua. Dangerous glacial lake and outburst features in Xizang Himalayas[J]. Acta Geographica Sinica, 1989, 44(3): 343-351,385. doi: 10.3321/j.issn:0375-5444.1989.03.011 [4] 吕儒仁, 唐邦兴, 朱平一, 等. 西藏泥石流与环境[M]. 成都: 成都科技大学出版社, 1999. [5] HUGGEL C, KÄÄB A, HAEBERLI W, et al. Remote sensing based assessment of hazards from glacier lake outbursts: a case study in the Swiss Alps[J]. Canadian Geotechnical Journal, 2002, 39(2): 316-330. doi: 10.1139/t01-099 [6] 陈晓清,崔鹏,杨忠,等. 喜马拉雅山中段波曲流域近期冰湖溃决危险性分析与评估[J]. 冰川冻土,2007,29(4): 509-516.CHEN Xiaoqing, CUI Peng, YANG Zhong, et al. Risk assessment of glacial lake outburst in the poiqu river basin of Tibet Autonomous Region[J]. Journal of Glaciology and Geocryology, 2007, 29(4): 509-516. [7] 舒有锋,王钢城,庄树裕,等. 基于粗糙集的权重确定方法在我国喜马拉雅山地区典型冰碛湖溃决危险性评价中的应用[J]. 水土保持通报,2010,30(5): 109-114. doi: 10.13961/j.cnki.stbctb.2010.05.008SHU Youfeng, WANG Gangcheng, ZHUANG Shuyu, et al. A weight determination method based on rough set for hazard assessment of typical moraine-dammed lake outburst in Himalayan region[J]. Bulletin of Soil and Water Conservation, 2010, 30(5): 109-114. doi: 10.13961/j.cnki.stbctb.2010.05.008 [8] 王欣,刘时银,郭万钦,等. 我国喜马拉雅山区冰碛湖溃决危险性评价[J]. 地理学报,2009,64(7): 782-790. doi: 10.3321/j.issn:0375-5444.2009.07.002WANG Xin, LIU Shiyin, GUO Wanqin, et al. Hazard assessment of moraine-dammed lake outburst floods in the Himalayas, China[J]. Acta Geographica Sinica, 2009, 64(7): 782-790. doi: 10.3321/j.issn:0375-5444.2009.07.002 [9] 马德仲,周真,于晓洋,等. 基于模糊概率的多状态贝叶斯网络可靠性分析[J]. 系统工程与电子技术,2012,34(12): 2607-2611. doi: 10.3969/j.issn.1001-506X.2012.12.35MA Dezhong, ZHOU Zhen, YU Xiaoyang, et al. Reliability analysis of multi-state Bayesian networks based on fuzzy probability[J]. Systems Engineering and Electronics, 2012, 34(12): 2607-2611. doi: 10.3969/j.issn.1001-506X.2012.12.35 [10] 林鹏智,陈宇. 基于贝叶斯网络的梯级水库群漫坝风险分析[J]. 工程科学与技术,2018,50(3): 46-53.LIN Pengzhi, CHEN Yu. Risk analysis of dam overtopping for cascade reservoirs based on Bayesian network[J]. Advanced Engineering Sciences, 2018, 50(3): 46-53. [11] 程尊兰,田金昌,张正波,等. 藏东南冰湖溃决泥石流形成的气候因素与发展趋势[J]. 地学前缘,2009,16(6): 207-214. doi: 10.3321/j.issn:1005-2321.2009.06.023CHENG Zunlan, TIAN Jinchang, ZHANG Zhengbo, et al. Debris flow induced by glacial-lake break in Southeast Tibet[J]. Earth Science Frontiers, 2009, 16(6): 207-214. doi: 10.3321/j.issn:1005-2321.2009.06.023 [12] 姚晓军,刘时银,孙美平,等. 20世纪以来西藏冰湖溃决灾害事件梳理[J]. 自然资源学报,2014,29(8): 1377-1390. [13] NIE Y, LIU Q, WANG J, et al. An inventory of historical glacial lake outburst floods in the Himalayas based on remote sensing observations and geomorphological analysis[J]. Geomorphology, 2018, 308: 91-106. [14] ICIMOD. Glacial lakes and glacial lake outburst floods in Nepal[R]. Kathmandu: ICIMOD, 2011. [15] KOMORI J, KOIKE T, YAMANOKUCHI T, et al. Glacial lake outburst events in the Bhutan Himalayas[J]. Global Environmental Research, 2012, 16(1): 59-70. [16] BHAMBRI R, HEWITT K, KAWISHWAR P, et al. Ice-dams, outburst floods, and movement heterogeneity of glaciers, Karakoram[J]. Global and Planetary Change, 2019, 180: 100-116. doi: 10.1016/j.gloplacha.2019.05.004 [17] GURUNG D R, KHANAL N R, BAJRACHARYA S R, et al. Lemthang Tsho glacial lake outburst flood (GLOF) in Bhutan: cause and impact[J]. Geoenvironmental Disasters, 2017, 4(17): 1-13. doi: 10.1186/s40677-017-0080-2 [18] AI H Z, YAO L K, ZHOU Y L. Laboratory investigations of earthquake-and landslide-induced composite surges[J]. Journal of Mountain Science, 2017, 14(8): 1537-1549. doi: 10.1007/s11629-016-4339-y [19] 刘宁, 程尊兰, 崔鹏. 堰塞湖及其风险控制[M]. 北京: 科学出版社, 2013. [20] WASHAKH R, CHEN N, WANG T, et al. GLOF risk assessment model in the Himalayas: a case study of a Hydropower project in the upper Arun river[J]. Water, 2019, 11(9): 1-23. doi: 10.3390/w11091839 期刊类型引用(10)
1. 王浩震. 波动载荷下旋挖钻机动力头液压系统控制技术. 液压气动与密封. 2024(02): 86-91 . 
百度学术2. 张彩婷，李化龙，吕阳伟，孙赟. 基于改进NSGA-Ⅲ算法的抽油机井群控变频启动技术. 化工自动化及仪表. 2024(02): 319-324+337 . 百度学术3. 刘淑妍，张爱莉. 稻麦联合收割机变频驱动液压控制中强扰动问题解决方法. 液压气动与密封. 2024(04): 87-92 . 百度学术4. 刘桓龙. 电液混合动力系统关键技术及能量管理研究综述. 西南交通大学学报. 2024(03): 600-614 . 
本站查看5. 孙宝林，崔凯，张千朋，韩朋刚，只升震. 基于PLC的船舶电站发电机组自动启动控制方法. 自动化应用. 2024(12): 111-113 . 百度学术6. 邓亚伟，王康，戴明明，杜海红. 基于电力模拟信号的井式电炉启动控制方法研究. 工业加热. 2024(08): 16-21 . 百度学术7. 高爽，张忠，容慧娴. 基于电流源型变频器的电机长线驱动方法. 自动化技术与应用. 2024(10): 176-180 . 百度学术8. 南亮生，陈璐，冯铁成，虞继敏. 基于改进EMD阈值的伺服电机驱动控制仿真. 计算机仿真. 2024(09): 346-350 . 百度学术9. 王燕. 基于载波频率的高频电流传感器均流控制系统设计. 计算机测量与控制. 2022(08): 111-115+121 . 百度学术10. 柯万宇，雷力，汤前进，彭瑞. 数控加工装备多驱动系统漏电流抑制方法. 机械与电子. 2022(10): 49-53 . 百度学术其他类型引用(1)
-
下载:
下载:
百度学术
点击查看大图
计量
- 文章访问数: 368
- HTML全文浏览量: 178
- PDF下载量: 29
- 被引次数: 11