发明名称 特高压直流输电线路行波保护雷击动作边界值确定方法
摘要 本发明公开了特高压直流输电线路行波保护雷击动作边界值确定方法。该方法按一定顺序排列,在线路上不同位置进行不同雷电流幅值的雷击仿真时得到的数据后,根据行波保护判据实际算法依次计算出各情况下判据,再以此判断行波保护动作情况,最后确定出不同位置下能令行波保护动作的最小雷电流幅值。本发明特高压直流输电线路行波保护雷击动作边界值确定方法是基于小步长仿真数据的特高压直流输电线路行波保护电压电流判据的计算方法,实现了按线路位置确定行波保护雷击动作边界值。
申请公布号 CN104569552A 申请公布日期 2015.04.29
申请号 CN201510011605.5 申请日期 2015.01.09
申请人 南方电网科学研究院有限责任公司;华南理工大学 发明人 李书勇;郭琦;饶宏;崔柳;黄立滨;蔡泽祥;梁益
分类号 G01R19/00(2006.01)I;G01R19/165(2006.01)I 主分类号 G01R19/00(2006.01)I
代理机构 广州市华学知识产权代理有限公司 44245 代理人 罗观祥
主权项 特高压直流输电线路行波保护雷击动作边界值确定方法,其特征在于包括如下步骤:1)按照到行波保护安装处的距离由近到远,将特高压直流输电线路作用点分别标记为N<sub>1</sub>、N<sub>2</sub>…N<sub>i</sub>…N<sub>n</sub>,则各个作用点到行波保护安装处的距离为D<sub>1</sub>、D<sub>2</sub>…D<sub>i</sub>…D<sub>n</sub>;2)在各作用点进行雷电流幅值分别为I<sub>1</sub>、I<sub>2</sub>…I<sub>k</sub>…I<sub>m</sub>的特高压直流输电线路雷击仿真,5kA≤I<sub>i</sub>≤35kA,幅值间隔1‐5kA,i为1、2…m;记录每次仿真时线路末端的电压U<sub>dl</sub>和线路末端安装的平波电抗器后侧的电流I<sub>dl</sub>,点N<sub>i</sub>的仿真数据记集合D<sub>i</sub>={I<sub>1</sub>(U<sub>dl1</sub>,I<sub>dl1</sub>)…I<sub>m</sub>(U<sub>dlm</sub>,I<sub>dlm</sub>)},I按照幅值从小到大排列;I<sub>1</sub>(U<sub>dl1</sub>,I<sub>dl1</sub>)为作用点N<sub>1</sub>处发生雷电流幅值分别为I<sub>1</sub>时的线路末端的电压U<sub>dl1</sub>和线路末端安装的平波电抗器后侧的电流I<sub>dl1</sub>;如此类推;3)计算作用点N<sub>i</sub>在发生雷电流幅值为I<sub>k</sub>情况下的行波保护电压变化率<img file="FDA0000654344800000011.GIF" wi="147" he="127" />电压变化量ΔU<sub>k</sub>、整流侧电流变化量ΔI<sub>dl_R_k</sub>及逆变侧电流变化量ΔI<sub>dl_I_k</sub>,并判断其出口情况;所述行波保护电压变化率<img file="FDA0000654344800000012.GIF" wi="118" he="134" />通过如下步骤得到:(1)用当前时刻电压采样值<img file="FDA0000654344800000013.GIF" wi="122" he="77" />与前150个时刻电压采样值<img file="FDA0000654344800000014.GIF" wi="166" he="77" />按照式6进行差分,得到当前时刻的电压变化率<img file="FDA0000654344800000015.GIF" wi="132" he="140" />其中下标T<sub>i</sub>表示时刻,ΔT取150μs;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><msub><mi>dU</mi><msub><mi>T</mi><mi>i</mi></msub></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><msub><mi>U</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mi>i</mi></msub></mrow></msub><mo>-</mo><msub><mi>U</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>-</mo><mn>150</mn></mrow></msub></mrow></msub></mrow><mi>&Delta;T</mi></mfrac></mrow>]]></math><img file="FDA0000654344800000016.GIF" wi="496" he="139" /></maths>                     式6(2)在当前T<sub>i</sub>时刻下,取T<sub>i</sub>以前的T<sub>i‑150</sub>、T<sub>i‑300</sub>、T<sub>i‑450</sub>、T<sub>i‑600</sub>、T<sub>i‑750</sub>时刻以及T<sub>i</sub>时刻的电压变化率计算值<img file="FDA0000654344800000017.GIF" wi="1166" he="141" />进行比较,取其中的最大值记为<img file="FDA0000654344800000018.GIF" wi="242" he="141" />作为当前时刻电压变化率;(3)重复以上步骤(1)至(2),直到所有时刻电压变化率<img file="FDA0000654344800000019.GIF" wi="94" he="133" />构成了计算作用点N<sub>i</sub>在发生雷电流幅值为I<sub>k</sub>情况下的行波保护电压变化率<img file="FDA0000654344800000021.GIF" wi="137" he="133" />所述电压变化量ΔU<sub>k</sub>通过如下方法确定:A)取当前时刻的线路末压电压采样值<img file="FDA0000654344800000022.GIF" wi="120" he="85" />与其以前的T<sub>i‑150</sub>、T<sub>i‑300</sub>、T<sub>i‑450</sub>、T<sub>i‑600</sub>、T<sub>i‑750</sub>这5个时刻的线路末端电压采样值进行比较,则比较<img file="FDA0000654344800000023.GIF" wi="799" he="85" /><img file="FDA0000654344800000024.GIF" wi="414" he="79" />然后取其最大值记为<img file="FDA0000654344800000025.GIF" wi="245" he="76" />B)将<img file="FDA0000654344800000026.GIF" wi="214" he="75" />与<img file="FDA0000654344800000027.GIF" wi="160" he="76" />作差后的值记为当前时刻的电压变化量<img file="FDA0000654344800000028.GIF" wi="140" he="77" />根据式7将<img file="FDA00006543448000000222.GIF" wi="108" he="75" />与电压变化量整定值ΔU<sub>set</sub>进行比较;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&Delta;U</mi><msub><mi>T</mi><mi>i</mi></msub></msub><mo>=</mo><mi>Max</mi><msub><mi>U</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mi>i</mi></msub></mrow></msub><mo>-</mo><msub><mi>U</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>-</mo><mn>5</mn></mrow></msub></mrow></msub></mrow>]]></math><img file="FDA0000654344800000029.GIF" wi="548" he="81" /></maths>                          式7C)根据比较结果,若<img file="FDA00006543448000000210.GIF" wi="114" he="75" />大于整定值ΔU<sub>set</sub>,则电压变化量判据出口;D)下一个时刻的<img file="FDA00006543448000000211.GIF" wi="127" he="76" />值计算则以令将上述步骤中计算所用的数据改为T下标为i=i+1时的值;如此重复以上步骤A)至C),得到的是不同时刻的<img file="FDA00006543448000000212.GIF" wi="112" he="83" />值,所有<img file="FDA00006543448000000213.GIF" wi="114" he="80" />构成的图像即为电压变化量判据ΔU<sub>k</sub>;基于采样时序,整流侧电流变化量判据ΔI<sub>dl_R</sub>及逆变侧电流变化量判据ΔI<sub>dl_I</sub>的计算及出口依据为:所述整流侧电流变化量ΔI<sub>dl_R_k</sub>及逆变侧电流变化量判据ΔI<sub>dl_I</sub>通过如下步骤得到:a)取当前时刻以前的第1200个时刻的电流采样值<img file="FDA00006543448000000214.GIF" wi="162" he="85" />作为当前时刻的电流采样值;b)使其与它当前时刻前1350时刻的电流采样值<img file="FDA00006543448000000215.GIF" wi="158" he="78" />作如式8的差分,获得当前时刻的电流变化率<img file="FDA00006543448000000216.GIF" wi="124" he="139" />ΔT取150μs;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfrac><msub><mi>dI</mi><msub><mi>T</mi><mi>i</mi></msub></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>-</mo><mn>1200</mn></mrow></msub></mrow></msub><mo>-</mo><msub><mi>I</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>-</mo><mn>1350</mn></mrow></msub></mrow></msub></mrow><mi>&Delta;T</mi></mfrac></mrow>]]></math><img file="FDA00006543448000000217.GIF" wi="507" he="142" /></maths>                  式8c)对<img file="FDA00006543448000000218.GIF" wi="102" he="140" />进行如式9的定积分处理后获得整流侧电流变化量判据<img file="FDA00006543448000000219.GIF" wi="189" he="75" />然后将<img file="FDA00006543448000000220.GIF" wi="160" he="74" />进行如式10的定积分处理后获得逆变侧电流变化量判据<img file="FDA00006543448000000221.GIF" wi="182" he="72" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&Delta;I</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mi>i</mi></msub><mo>_</mo><mi>R</mi></mrow></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mn>2</mn></mrow><mn>1.99</mn></msubsup><mfrac><msub><mi>dI</mi><msub><mi>T</mi><mi>i</mi></msub></msub><mi>d&tau;</mi></mfrac><mi>dt</mi></mrow>]]></math><img file="FDA0000654344800000031.GIF" wi="439" he="141" /></maths>                     式9<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&Delta;I</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mi>i</mi></msub><mo>_</mo><mi>I</mi></mrow></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mn>2</mn></mrow><mn>1.99</mn></msubsup><msub><mi>&Delta;I</mi><mrow><mi>dl</mi><mo>_</mo><msub><mi>T</mi><mi>i</mi></msub><mo>_</mo><mi>R</mi></mrow></msub><mi>dt</mi></mrow>]]></math><img file="FDA0000654344800000032.GIF" wi="522" he="111" /></maths>                     式10d)分别将<img file="FDA0000654344800000033.GIF" wi="161" he="70" />与其定值ΔI<sub>dl_Rset</sub>比价,将<img file="FDA0000654344800000034.GIF" wi="156" he="69" />与其定值ΔI<sub>dl_Iset</sub>比较,大于整定值的那一侧保护出口;e)下一个时刻的<img file="FDA0000654344800000035.GIF" wi="158" he="77" />值及<img file="FDA0000654344800000036.GIF" wi="160" he="76" />值计算时,将上述步骤中计算所用的数据T下标i进行i=i+1替换进行同样的计算;如此重复步骤a)至d),所有<img file="FDA0000654344800000037.GIF" wi="150" he="74" />构成的图像即为整流侧电流变化量判据ΔI<sub>dl_R</sub>;所有<img file="FDA0000654344800000038.GIF" wi="156" he="70" />构成的图像即为整流侧电流变化量判据ΔI<sub>dl_I</sub>;4)若存在<img file="FDA0000654344800000039.GIF" wi="208" he="136" />大于行波保护电压变化率定值<img file="FDA00006543448000000310.GIF" wi="130" he="129" />的情况,则认为该数据下的电压变化率<img file="FDA00006543448000000311.GIF" wi="118" he="146" />满足判据定值;若此时判据<img file="FDA00006543448000000312.GIF" wi="124" he="115" />ΔU<sub>k</sub>、ΔI<sub>dlk_I</sub>和ΔI<sub>dlk_R</sub>均不出口,则更换下一个I<sub>k</sub>对应的数据,重复步骤3),直到首先出现某个I<sub>k</sub>能令<img file="FDA00006543448000000313.GIF" wi="131" he="116" />ΔU<sub>k</sub>、ΔI<sub>dlk_I</sub>和ΔI<sub>dlk_R</sub>同时达到保护定值时,该值为特高压直流输电线路位置D<sub>i</sub>处的行波保护雷击动作边界值<img file="FDA00006543448000000314.GIF" wi="97" he="76" />5)更换下一个点N<sub>i</sub>的仿真数据集合D<sub>i</sub>,重复步骤3)、4),直到所有点N全部计算完毕为止;6)得到输出线路各位置行波保护雷击动作边界值<img file="FDA00006543448000000315.GIF" wi="643" he="80" /><img file="FDA00006543448000000316.GIF" wi="295" he="75" />
地址 510080 广东省广州市越秀区东风东路水均岗6、8号