1.本发明属于超导电缆输电领域，具体涉及一种超导电缆输电线路序参数的计算方法、系统、设备及介质。
背景技术：2.超导电缆具有低损耗和紧凑大容量的特征，是电网的潜在变革技术，国内外开展了大量的基础研究与并网试验示范。2000年美国南线公司并网运行世界第一条30m长、12.5kv/1.25ka的三相统包冷绝缘高温超导电缆，2014 年德国投运1公里长、10kv/1ka三相同轴超导电缆。目前，欧、美、日和韩等密集开展了数十条示范工程，系统验证超导电缆的技术性能。国内，2003 年英纳超导并网33.5米、35kv/2ka热绝缘高温超导电缆，初步探索了热绝缘超导电缆的并网技术；2013年上缆所在宝钢投运三相分体冷绝缘超导电缆，标志着国内基本掌握超导电缆制造技术；2021年，上海和深圳分别投运的超导电缆输电示范工程，标志着超导电缆在国内进入实际并网阶段。
3.典型的冷绝缘超导电缆包括内支撑、电缆导体、低温恒温器、电气绝缘层和电缆屏蔽层与外护套等部分。冷绝缘超导电缆与常规电缆最大的不同是其屏蔽层，一般采用与导体层相同的超导带材同轴绕制而成，运行时屏蔽层也处在低温环境中，称为超导屏蔽层。冷绝缘超导电缆屏蔽层的作用是屏蔽超导电缆导体输电时产生的交变磁场，屏蔽交变磁场一方面防止在低温恒温器等金属导体中产生感应电流，避免涡流损耗，减少制冷功率，降低安全风险；另一方面减少相间的间隔距离，节约电缆敷设占用的管廊空间，提高空间的利用效率。
4.与常规电缆相比，超导屏蔽层有两方面的影响。首先，在线路内部和外部的非对称故障下，传输的三相电流之和不为零，因此屏蔽层入地电流不为零，如图1所示，这样会改变电网系统的零序电流。其次，超导屏蔽层的电阻很小，且与导电层耦合较好，一般情况下超导屏蔽层的感应电流为传输电流的90％左右，感应电流的存在减少了电缆线路电抗。这两方面的影响都与超导电缆的线路序参数相关，因此，有必要建立考虑冷绝缘超导屏蔽层结构的超导电缆线路序参数分析方法。
技术实现要素：5.本发明的目的在于提供一种超导电缆输电线路序参数的计算方法、系统、设备及介质，以克服现有技术存在的缺陷，本发明建立一种计算冷绝缘超导输电线路参数的分析方法，具有重要的理论价值和实践意义。
6.为达到上述目的，本发明采用如下技术方案：
7.超导电缆输电线路序参数的计算方法，包括：
8.获取超导电缆的结构参数和电气物性参数，所述电气物性参数包括磁导率、电导率和电阻；
9.基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算
屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和超导电缆的电阻计算导电层和屏蔽层传输电流；
10.通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数。
11.进一步地，所述超导电缆的二维等效磁场模型如式(1)所示：
[0012][0013]
式中：μ为磁场的磁导率；je为外部激励电流的密度矢量；js为感应电流密度矢量；a为矢量磁位；为旋度计算符号。
[0014]
进一步地，所述基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量，具体为：
[0015]
依据超导电缆的结构参数，采用棱边有限元法求解二维等效磁场模型磁场区域的整体插值函数如式(2)所示：
[0016][0017]
式中：a为矢量磁位，在此利用整体插值函数表示矢量磁位；nn(x,y,z) 为单元插值函数；n为第n个单元的单元编号；n1为总单元数；an表示每个单元用插值函数表示的矢量磁位；
[0018]
依据格林定理对式(2)求解，得到相应磁场区域的伽辽金加权余量方程：
[0019][0020]
式中：nm，nn分别表示第m个单元和第n个单元的单元插值函数；m表示第m个单元的单元编号；j为电流密度；
[0021]
二维等效磁场模型中若已知电流分布，将伽辽金加权余量方程离散形成代数方程组，求解代数方程组得到超导电缆横截面空间的矢量磁位a；
[0022]
利用矢量磁位a求得磁感应强度b，如式(4)所示：
[0023][0024]
利用电磁感应定律求解空间感应电场e，如式(5)所示：
[0025][0026]
通过欧姆定律获得感应电流密度矢量js，如式(6)所示：
[0027]js
＝σe
ꢀꢀꢀ
(6)
[0028]
式中，σ为电导率。
[0029]
进一步地，所述等效电路模型包括a，b和c三相电源，a，b和c三相电源分别通过电阻rs与超导电缆的导电层组成回路，形成导电层回路，电阻rs包括电源内阻和负载电阻；a、b和c三相的屏蔽层在两端短接后串入电阻rg，形成屏蔽层回路，电阻rg为接地电阻。
[0030]
5.根据权利要求4所述的超导电缆输电线路序参数的计算方法，其特征在于，所述导电层回路的电路方程如式(7)所示：
[0031][0032]
式中：ea、eb和ec分别是导电层回路a、b和c相的电源电势；ia、ib和ic分别是导电层回路a、b和c相的传输电流；ua、ub和uc分别是a、 b和c相的超导电缆导电层两端的电压降，由二维等效磁场模型获得；
[0033]
所述屏蔽层回路的电路方程如式(8)所示：
[0034][0035]
式中：ua、ub和uc分别是屏蔽层回路a、b和c三相感应电动势；ia、ib和ic分别是屏蔽层回路a、b和c三相的传输电流，由二维等效磁场模型获得； ig是屏蔽层回路入地电流，当ia，ib和ic三相对称时，入地电流ig为0。
[0036]
进一步地，所述通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数，具体包括：
[0037]
1)已知某一时刻的电源电势值，在等效电路模型中求解超导电缆导电层和屏蔽层传输电流；
[0038]
2)将导电层和屏蔽层传输电流代入二维等效磁场模型，通过求解式(1) 和式(4)获得整体空间磁场的矢量磁位a和磁感应强度b，然后通过式(5) 和式(6)得到导电层和屏蔽层电压降u和感应电流密度矢量js；
[0039]
3)将导电层和屏蔽层的电压降u代入等效电路模型，通过迭代步骤1) 和步骤2)，直到导电层电流和屏蔽层电流满足收敛条件后，进入步骤4)；
[0040]
4)将导电层电流和屏蔽层电流进行拟合，获得幅值和相位；
[0041]
5)将电压和电流的向量值相除得到序参数，所述向量值由幅值和相位构成。
[0042]
进一步地，所述超导电缆输电线路序参数计算包括零序阻抗计算以及正序阻抗和负序阻抗计算；
[0043]
所述正序阻抗的计算方法如下：
[0044]
设定激励电源电势ea、eb和ec为三相对称，相位互差120
°
，通过对a、 b和c相区域的电流密度进行积分，获得ia、ib和ic随时间变化的曲线，将 ia、ib和ic随时间变化的曲线，截取稳态部分，按照正弦曲线进行拟合，如式(9)所示：
[0045]
i1＝a
1 sin(b1t+c1)
ꢀꢀꢀ
(9)
[0046]
式中：t为时间；a1表示正序电流幅值；b1为正序角频率；c1为正序相位，拟合获得的i1为导电层回路a、b和c相的传输电流；
[0047]
假设超导电缆的正序阻抗为：
[0048]
z1＝r1+jx1ꢀꢀꢀ
(10)
[0049]
式中，z1为正序阻抗；r1为正序电阻；x1为正序电抗；j为虚数；
[0050]
将式(9)和式(10)带入式(7)得到：
[0051][0052]
将拟合后的式(9)带入式(11)，则正序阻抗公式为：
[0053][0054]
所述负序阻抗与正序阻抗相等；
[0055]
所述零序阻抗的计算方法如下：
[0056]
设定激励电源电势ea、eb和ec的幅值相等、相位相同进行仿真计算，积分获得ia、ib和ic随时间变化的曲线，截取稳态部分后拟合，如式(13) 所示：
[0057]
i0＝a
0 sin(b0t+c0)
ꢀꢀꢀ
(13)
[0058]
式中，a0表示零序电流幅值，b0为零序角频率，c0为零序相位，拟合获得的i0即为a、b或c相的零序电流；
[0059]
令超导电缆的零序阻抗为：
[0060]
z0＝r0+jx0ꢀꢀꢀ
(14)
[0061]
式中，z0为零序阻抗；r0为零序电阻；x0为零序电抗；
[0062]
假设超导电缆两端的电压为u0，则：
[0063]
u0＝e
a-i0rsꢀꢀꢀ
(15)
[0064]
零序阻抗为：
[0065][0066]
将拟合后的式(13)带入式(16)并结合式(15)，计算获得零序阻抗为：
[0067][0068]
超导电缆输电线路序参数的计算系统，包括：
[0069]
参数获取模块：用于获取超导电缆的结构参数和电气物性参数，所述电气物性参数包括磁导率、电导率和电阻；
[0070]
中间量计算模块：用于基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和超导电缆的电阻计算导电层和屏蔽层传输电流；
[0071]
序参数计算模块：用于通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数。
[0072]
进一步地，所述超导电缆的二维等效磁场模型如式(1)所示：
[0073][0074]
式中：μ为磁场的磁导率；je为外部激励电流的密度矢量；js为感应电流密度矢量；a为矢量磁位；为旋度计算符号。
[0075]
进一步地，所述基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量，具体为：
[0076]
依据超导电缆的结构参数，采用棱边有限元法求解二维等效磁场模型磁场区域的整体插值函数如式(2)所示：
[0077][0078]
式中：a为矢量磁位，在此利用整体插值函数表示矢量磁位；nn(x,y,z) 为单元插值函数；n为第n个单元的单元编号；n1为总单元数；an表示每个单元用插值函数表示的矢量磁位；
[0079]
依据格林定理对式(2)求解，得到相应磁场区域的伽辽金加权余量方程：
[0080][0081]
式中：nm，nn分别表示第m个单元和第n个单元的单元插值函数；m表示第m个单元的单元编号；j为电流密度；
[0082]
二维等效磁场模型中若已知电流分布，将伽辽金加权余量方程离散形成代数方程组，求解代数方程组得到超导电缆横截面空间的矢量磁位a；
[0083]
利用矢量磁位a求得磁感应强度b，如式(4)所示：
[0084][0085]
利用电磁感应定律求解空间感应电场e，如式(5)所示：
[0086][0087]
通过欧姆定律获得感应电流密度矢量js，如式(6)所示：
[0088]js
＝σe
ꢀꢀꢀ
(6)
[0089]
式中，σ为电导率。
[0090]
进一步地，所述等效电路模型包括a，b和c三相电源，a，b和c三相电源分别通过电阻rs与超导电缆的导电层组成回路，形成导电层回路，电阻rs包括电源内阻和负载电阻；a、b和c三相的屏蔽层在两端短接后串入电阻rg，形成屏蔽层回路，电阻rg为接地电阻。
[0091]
进一步地，所述导电层回路的电路方程如式(7)所示：
[0092][0093]
式中：ea、eb和ec分别是导电层回路a、b和c相的电源电势；ia、ib和ic分别是导电层回路a、b和c相的传输电流；ua、ub和uc分别是a、 b和c相的超导电缆导电层两端的电压降，由二维等效磁场模型获得；
[0094]
所述屏蔽层回路的电路方程如式(8)所示：
[0095][0096]
式中：ua、ub和uc分别是屏蔽层回路a、b和c三相感应电动势；ia、ib和ic分别是屏蔽层回路a、b和c三相的传输电流，由二维等效磁场模型获得； ig是屏蔽层回路入地电流，当ia，ib和ic三相对称时，入地电流ig为0。
[0097]
进一步地，所述通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数，具体包括：
[0098]
1)已知某一时刻的电源电势值，在等效电路模型中求解超导电缆导电层和屏蔽层传输电流；
[0099]
2)将导电层和屏蔽层传输电流代入二维等效磁场模型，通过求解式(1) 和式(4)获得整体空间磁场的矢量磁位a和磁感应强度b，然后通过式(5) 和式(6)得到导电层和屏蔽层电压降u和感应电流密度矢量js；
[0100]
3)将导电层和屏蔽层的电压降u代入等效电路模型，通过迭代步骤1) 和步骤2)，直到导电层电流和屏蔽层电流满足收敛条件后，进入步骤4)；
[0101]
4)将导电层电流和屏蔽层电流进行拟合，获得幅值和相位；
[0102]
5)将电压和电流的向量值相除得到序参数，所述向量值由幅值和相位构成。
[0103]
进一步地，所述超导电缆输电线路序参数计算包括零序阻抗计算以及正序阻抗和负序阻抗计算；
[0104]
所述正序阻抗的计算方法如下：
[0105]
设定激励电源电势ea、eb和ec为三相对称，相位互差120
°
，通过对a、 b和c相区域的电流密度进行积分，获得ia、ib和ic随时间变化的曲线，将 ia、ib和ic随时间变化的曲线，截取稳态部分，按照正弦曲线进行拟合，如式(9)所示：
[0106]
i1＝a
1 sin(b1t+c1)
ꢀꢀꢀ
(9)
[0107]
式中：t为时间；a1表示正序电流幅值；b1为正序角频率；c1为正序相位，拟合获得的i1为导电层回路a、b和c相的传输电流；
[0108]
假设超导电缆的正序阻抗为：
[0109]
z1＝r1+jx1ꢀꢀꢀ
(10)
[0110]
式中，z1为正序阻抗；r1为正序电阻；x1为正序电抗；j为虚数；
[0111]
将式(9)和式(10)带入式(7)得到：
[0112][0113]
将拟合后的式(9)带入式(11)，则正序阻抗公式为：
[0114]
[0115]
所述负序阻抗与正序阻抗相等；
[0116]
所述零序阻抗的计算方法如下：
[0117]
设定激励电源电势ea、eb和ec的幅值相等、相位相同进行仿真计算，积分获得ia、ib和ic随时间变化的曲线，截取稳态部分后拟合，如式(13) 所示：
[0118]
i0＝a
0 sin(b0t+c0)
ꢀꢀꢀ
(13)
[0119]
式中，a0表示零序电流幅值，b0为零序角频率，c0为零序相位，拟合获得的i0即为a、b或c相的零序电流；
[0120]
令超导电缆的零序阻抗为：
[0121]
z0＝r0+jx0ꢀꢀꢀ
(14)
[0122]
式中，z0为零序阻抗；r0为零序电阻；x0为零序电抗；
[0123]
假设超导电缆两端的电压为u0，则：
[0124]
u0＝e
a-i0rsꢀꢀꢀ
(15)
[0125]
零序阻抗为：
[0126][0127]
将拟合后的式(13)带入式(16)并结合式(15)，计算获得零序阻抗为：
[0128][0129]
一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述超导电缆输电线路序参数的计算方法的步骤。
[0130]
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述超导电缆输电线路序参数的计算方法的步骤。
[0131]
与现有技术相比，本发明具有以下有益的技术效果：
[0132]
与常规电缆不同，超导输电的零电阻特性和导电层/屏蔽层强耦合特性，使得超导输电的线路参数特性有新的特点，本发明基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和超导电缆的电阻计算导电层和屏蔽层传输电流，然后通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导输电线路的正序、负序和零序线路参数，能够为计算系统短路电流、继电保护整定、电力系统潮流计算和选择合理运行方式等工作提供重要依据，具有重要的理论价值和实践意义。
附图说明
[0133]
说明书附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
[0134]
图1为三芯超导电缆系统接线示意图；
[0135]
图2为本发明方法流程示意图；
[0136]
图3为本发明系统结构示意图；
[0137]
图4为冷绝缘三相超导电缆磁场计算模型示意图；
[0138]
图5为三相超导电缆的电路模型示意图；
[0139]
图6为电缆正序电抗与屏蔽层电阻的关系图；
[0140]
图7为电缆零序电抗与屏蔽层电阻的关系图。
具体实施方式
[0141]
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
[0142]
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0143]
本发明基于超导电缆时域场路耦合电磁模型进行线路序参数计算分析，如图2，其中等效电路模型计算导电层和屏蔽层传输电流，二维等效磁场模型计算屏蔽层感应电场和感应电流密度矢量，由屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆的线路序参数。
[0144]
本发明针对超导电缆的时域场路耦合电磁模型有如下假设：
[0145]
1)本发明的时域场路耦合电磁模型不考虑超导的非线性电阻对线路阻抗的影响，侧重考查线路的电抗分析方法。
[0146]
2)忽略入地电流的磁场效应，即不考虑入地电流与超导导电层与超导屏蔽层间的互感。
[0147]
3)忽略超导电缆绕制参数的影响，采用电缆横截面的二维模型作为磁场仿真对象。
[0148]
二维等效磁场模型
[0149]
基于二维等效磁场模型分析超导电缆导电层与屏蔽层的电磁感应关系，如图4所示，采用棱边有限元法以及矢量磁位a，建立超导电缆的二维等效磁场模型：
[0150][0151]
式中：μ为磁场的磁导率，表征媒介的磁性；je为外部激励电流的密度矢量，表征激励电流的大小和方向；js为感应电流密度矢量，表征感应电流的大小和方向；为旋度计算符号。
[0152]
依据超导电缆的结构参数，采用棱边有限元法求解二维等效磁场模型磁场区域的整体插值函数为
[0153][0154]
式中：a为矢量磁位，在此利用整体插值函数表示矢量磁位；nn(x,y,z) 为单元插值函数；n为第n个单元的单元编号；n1为总单元数；an表示每个单元用插值函数表示的矢量磁位。
[0155]
依据格林定理对式(2)求解，可得相应磁场区域的伽辽金加权余量方程：
[0156][0157]
式中：nm，nn分别表示第m个单元和第n个单元的单元插值函数；m表示第m个单元的单元编号；j为电流密度。
[0158]
二维等效磁场模型中若已知电流分布，将伽辽金加权余量方程离散形成代数方程组，求解可得超导电缆横截面空间的矢量磁位a。
[0159]
利用矢量磁位a求得磁感应强度b：
[0160][0161]
利用电磁感应定律求解空间感应电场e
[0162][0163]
获得屏蔽层的感应电场e后，通过欧姆定律可以获得感应电流密度矢量 js；
[0164]js
＝σe
ꢀꢀꢀ
(6)
[0165]
式中σ为电导率。
[0166]
等效电路模型
[0167]
等效电路模型如图5所示，三相电源通过电阻rs与超导电缆的导电层组成回路，称为导电层回路，电阻rs包括电源内阻和负载电阻；a、b和c三相的屏蔽层在两端短接后串入电阻rg，称为屏蔽层回路，电阻rg主要包括地电阻。在本发明的仿真计算中，电缆长度l为1km。
[0168]
导电层回路的电路方程为
[0169][0170]
式中：ea、eb和ec分别是导电层回路a、b和c相的电源电势；ia、ib和ic分别是a、b和c相回路的传输电流；ua、ub和uc分别是a、b和c 相的超导电缆导电层两端的电压降，由二维等效磁场模型获得。
[0171]
同理，屏蔽层回路的电路方程为
[0172][0173]
式中：ua、ub和uc分别是屏蔽层回路a、b和c三相感应电动势；ia、 ib和ic分别是屏蔽层回路a、b和c三相的传输电流，通过磁场模型获得；ig是屏蔽层回路入地电流，显然当ia，ib和ic对称三相时，入地电流ig为0。
[0174]
时域场路耦合电磁模型的求解包括以下几个步骤：
[0175]
1)已知某一时刻的电源电势值，在等效电路模型中求解超导电缆导电层和屏蔽层传输电流。
[0176]
2)将导电层和屏蔽层传输电流代入二维等效磁场模型，通过求解式(1) 和式(4)获得整体空间磁场的矢量磁位a和磁感应强度b，然后通过式(5) 和式(6)得到导电层和屏蔽层电压降u(单位长度下与与空间感应电场e 数值相等)和感应电流密度矢量js。
[0177]
3)将导电层和屏蔽层的电压降u代入等效电路模型，通过迭代步骤1) 和步骤2)，直到导电层电流和屏蔽层电流满足收敛条件后，进入步骤4)；
[0178]
4)将导电层电流和屏蔽层电流进行拟合，获得幅值和相位；
[0179]
5)将电压和电流的向量值相除得到序参数，所述向量值由幅值和相位构成。
[0180]
正序参数计算方法
[0181]
超导电缆序参数的提取方法参考工程中线路序参数的实测方法。工程上通过短接线路一端的a、b和c三相，另一端注入三相电压源，测量线路的电压、电流和功率来获得线路的序参数。
[0182]
在时域场路耦合电磁模型中，基本思路是通过三相电压源的激励，获得超导电缆a、b和c相的电流，然后计算出超导电缆的等效正负序和等效零序参数。
[0183]
对电缆而言，正序阻抗与负序阻抗相等，后面以正序阻抗为例进行说明。
[0184]
在正序参数计算中，设定激励电源电势ea、eb和ec为三相对称，相位互差120
°
，经过场路耦合模型仿真计算后得到磁场分布，通过对a、b和c 相区域的电流密度进行积分，获得ia、ib和ic随时间变化的曲线。将ia、ib和ic随时间的曲线，截取稳态部分，按照正弦曲线进行拟合：
[0185]
i1＝a
1 sin(b1t+c1)
ꢀꢀꢀ
(9)
[0186]
式中：t为时间；a1表示正序电流幅值；b1为正序角频率；c1为正序相位，拟合获得的i1为导电层回路a、b和c相的传输电流。
[0187]
假设超导电缆的正序阻抗为
[0188]
z1＝r1+jx1ꢀꢀꢀ
(10)
[0189]
式中，z1为正序阻抗；r1为正序电阻；x1为正序电抗；j为虚数。
[0190]
将式(9)和式(10)带入式(7)可得：
[0191][0192]
将拟合后的式(9)带入式(11)，则阻抗公式为
[0193][0194]
零序参数计算方法
[0195]
在零序参数计算中，设定激励电源电势ea、eb和ec为幅值相等，相位相同进行仿真计算，积分获得ia、ib和ic随时间变化的曲线，截取稳态部分后拟合：
[0196]
i0＝a
0 sin(b0t+c0)
ꢀꢀꢀ
(13)
[0197]
同理a0表示零序电流幅值，b0为零序角频率，c0为零序相位，拟合获得的i0即为a、b或c相的零序电流。
[0198]
令超导电缆的零序阻抗为
[0199]
z0＝r0+jx0ꢀꢀꢀ
(14)
[0200]
式中，z0为零序阻抗；r0为零序电阻；x0为零序电抗；
[0201]
假设超导电缆两端的电压为u0，则
[0202]
u0＝e
a-i0rsꢀꢀꢀ
(15)
[0203]
零序参数为
[0204][0205]
将拟合后的式(13)带入式(16)结合式(15)，计算获得的零序阻抗为
[0206][0207]
本发明还提供一种超导电缆输电线路序参数的计算系统，如图3，包括：
[0208]
参数获取模块：用于获取超导电缆的结构参数和电气物性参数，所述电气物性参数包括磁导率、电导率和电阻；
[0209]
中间量计算模块：用于基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和超导电缆的电阻计算导电层和屏蔽层传输电流；
[0210]
序参数计算模块：用于通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数。
[0211]
实施例
[0212]
基于本发明的方法，根据表1参数结合有限元分析方法和电路分析方法，可以计算获得超导输电线路的序参数，如表2。
[0213]
表1 超导电缆的结构参数和电气物性参数
[0214][0215][0216]
当屏蔽层电阻取值0时，求取超导输电线路的序参数，然后通过改变屏蔽层电阻取值，获得电缆正序电抗与屏蔽层电阻的关系(如图6)以及电缆零序电抗与屏蔽层电阻的关系(如图7)。
[0217]
超导电缆的正序阻抗参数和零序阻抗参数的准确性关系到电网的安全稳定运行，是计算系统短路电流、继电保护整定、电力系统潮流计算和选择合理运行方式等工作的依据，本发明通过计算超导输电线路的正序、负序和零序线路参数，能够为计算系统短路电流、继电保护整定、电力系统潮流计算和选择合理运行方式等工作提供重要依据，具有重要的理论价值和实践意义。
[0218]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0219]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流
程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0220]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0221]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0222]
最后应当说明的是：以上实施例仅用于说明本发明的技术方案而非对其保护范围的限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解:本领域技术人员阅读本发明后依然可对发明的具体实施方式进行种种变更、修改或者等同替换，但这些变更、修改或者等同替换，均在发明待批的权利要求保护范围之内。
技术特征：1.超导电缆输电线路序参数的计算方法，其特征在于，包括：获取超导电缆的结构参数和电气物性参数，所述电气物性参数包括磁导率、电导率和电阻；基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和超导电缆的电阻计算导电层和屏蔽层传输电流；通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数。2.根据权利要求1所述的超导电缆输电线路序参数的计算方法，其特征在于，所述超导电缆的二维等效磁场模型如式(1)所示：式中：μ为磁场的磁导率；j
e
为外部激励电流的密度矢量；j
s
为感应电流密度矢量；a为矢量磁位；为旋度计算符号。3.根据权利要求2所述的超导电缆输电线路序参数的计算方法，其特征在于，所述基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量，具体为：依据超导电缆的结构参数，采用棱边有限元法求解二维等效磁场模型磁场区域的整体插值函数如式(2)所示：式中：a为矢量磁位，在此利用整体插值函数表示矢量磁位；n
n
(x,y,z)为单元插值函数；n为第n个单元的单元编号；n1为总单元数；a
n
表示每个单元用插值函数表示的矢量磁位；依据格林定理对式(2)求解，得到相应磁场区域的伽辽金加权余量方程：式中：n
m
，n
n
分别表示第m个单元和第n个单元的单元插值函数；m表示第m个单元的单元编号；j为电流密度；二维等效磁场模型中若已知电流分布，将伽辽金加权余量方程离散形成代数方程组，求解代数方程组得到超导电缆横截面空间的矢量磁位a；利用矢量磁位a求得磁感应强度b，如式(4)所示：利用电磁感应定律求解空间感应电场e，如式(5)所示：通过欧姆定律获得感应电流密度矢量j
s
，如式(6)所示：j
s
＝σe
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)式中，σ为电导率。
4.根据权利要求3所述的超导电缆输电线路序参数的计算方法，其特征在于，所述等效电路模型包括a，b和c三相电源，a，b和c三相电源分别通过电阻r
s
与超导电缆的导电层组成回路，形成导电层回路，电阻r
s
包括电源内阻和负载电阻；a、b和c三相的屏蔽层在两端短接后串入电阻r
g
，形成屏蔽层回路，电阻r
g
为接地电阻。5.根据权利要求4所述的超导电缆输电线路序参数的计算方法，其特征在于，所述导电层回路的电路方程如式(7)所示：式中：e
a
、e
b
和e
c
分别是导电层回路a、b和c相的电源电势；i
a
、i
b
和i
c
分别是导电层回路a、b和c相的传输电流；u
a
、u
b
和u
c
分别是a、b和c相的超导电缆导电层两端的电压降，由二维等效磁场模型获得；所述屏蔽层回路的电路方程如式(8)所示：式中：u
a
、u
b
和u
c
分别是屏蔽层回路a、b和c三相感应电动势；i
a
、i
b
和i
c
分别是屏蔽层回路a、b和c三相的传输电流，由二维等效磁场模型获得；i
g
是屏蔽层回路入地电流，当i
a
，i
b
和i
c
三相对称时，入地电流i
g
为0。6.根据权利要求5所述的超导电缆输电线路序参数的计算方法，其特征在于，所述通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数，具体包括：1)已知某一时刻的电源电势值，在等效电路模型中求解超导电缆导电层和屏蔽层传输电流；2)将导电层和屏蔽层传输电流代入二维等效磁场模型，通过求解式(1)和式(4)获得整体空间磁场的矢量磁位a和磁感应强度b，然后通过式(5)和式(6)得到导电层和屏蔽层电压降u和感应电流密度矢量j
s
；3)将导电层和屏蔽层的电压降u代入等效电路模型，通过迭代步骤1)和步骤2)，直到导电层电流和屏蔽层电流满足收敛条件后，进入步骤4)；4)将导电层电流和屏蔽层电流进行拟合，获得幅值和相位；5)将电压和电流的向量值相除得到序参数，所述向量值由幅值和相位构成。7.根据权利要求6所述的超导电缆输电线路序参数的计算方法，其特征在于，所述超导电缆输电线路序参数计算包括零序阻抗计算以及正序阻抗和负序阻抗计算；所述正序阻抗的计算方法如下：设定激励电源电势e
a
、e
b
和e
c
为三相对称，相位互差120
°
，通过对a、b和c相区域的电流密度进行积分，获得i
a
、i
b
和i
c
随时间变化的曲线，将i
a
、i
b
和i
c
随时间变化的曲线，截取稳态部分，按照正弦曲线进行拟合，如式(9)所示：i1＝a
1 sin(b1t+c1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)式中：t为时间；a1表示正序电流幅值；b1为正序角频率；c1为正序相位，拟合获得的i1为
导电层回路a、b和c相的传输电流；假设超导电缆的正序阻抗为：z1＝r1+jx1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)式中，z1为正序阻抗；r1为正序电阻；x1为正序电抗；j为虚数；将式(9)和式(10)带入式(7)得到：将拟合后的式(9)带入式(11)，则正序阻抗公式为：所述负序阻抗与正序阻抗相等；所述零序阻抗的计算方法如下：设定激励电源电势e
a
、e
b
和e
c
的幅值相等、相位相同进行仿真计算，积分获得i
a
、i
b
和i
c
随时间变化的曲线，截取稳态部分后拟合，如式(13) 所示：i0＝a0sin(b0t+c0)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)式中，a0表示零序电流幅值，b0为零序角频率，c0为零序相位，拟合获得的i0即为a、b或c相的零序电流；令超导电缆的零序阻抗为：z0＝r0+jx0ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)式中，z0为零序阻抗；r0为零序电阻；x0为零序电抗；假设超导电缆两端的电压为u0，则：u0＝e
a-i0r
s
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(15)零序阻抗为：将拟合后的式(13)带入式(16)并结合式(15)，计算获得零序阻抗为：8.超导电缆输电线路序参数的计算系统，其特征在于，包括：参数获取模块：用于获取超导电缆的结构参数和电气物性参数，所述电气物性参数包括磁导率、电导率和电阻；中间量计算模块：用于基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和
超导电缆的电阻计算导电层和屏蔽层传输电流；序参数计算模块：用于通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数。9.根据权利要求8所述的超导电缆输电线路序参数的计算系统，其特征在于，所述超导电缆的二维等效磁场模型如式(1)所示：式中：μ为磁场的磁导率；j
e
为外部激励电流的密度矢量；j
s
为感应电流密度矢量；a为矢量磁位；为旋度计算符号。10.根据权利要求9所述的超导电缆输电线路序参数的计算系统，其特征在于，所述基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量，具体为：依据超导电缆的结构参数，采用棱边有限元法求解二维等效磁场模型磁场区域的整体插值函数如式(2)所示：式中：a为矢量磁位，在此利用整体插值函数表示矢量磁位；n
n
(x,y,z)为单元插值函数；n为第n个单元的单元编号；n1为总单元数；a
n
表示每个单元用插值函数表示的矢量磁位；依据格林定理对式(2)求解，得到相应磁场区域的伽辽金加权余量方程：式中：n
m
，n
n
分别表示第m个单元和第n个单元的单元插值函数；m表示第m个单元的单元编号；j为电流密度；二维等效磁场模型中若已知电流分布，将伽辽金加权余量方程离散形成代数方程组，求解代数方程组得到超导电缆横截面空间的矢量磁位a；利用矢量磁位a求得磁感应强度b，如式(4)所示：利用电磁感应定律求解空间感应电场e，如式(5)所示：通过欧姆定律获得感应电流密度矢量j
s
，如式(6)所示：j
s
＝σe
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)式中，σ为电导率。11.根据权利要求10所述的超导电缆输电线路序参数的计算系统，其特征在于，所述等效电路模型包括a，b和c三相电源，a，b和c三相电源分别通过电阻r
s
与超导电缆的导电层组成回路，形成导电层回路，电阻r
s
包括电源内阻和负载电阻；a、b和c三相的屏蔽层在两端短接后串入电阻r
g
，形成屏蔽层回路，电阻r
g
为接地电阻。12.根据权利要求11所述的超导电缆输电线路序参数的计算系统，其特征在于，所述导电层回路的电路方程如式(7)所示：
式中：e
a
、e
b
和e
c
分别是导电层回路a、b和c相的电源电势；i
a
、i
b
和i
c
分别是导电层回路a、b和c相的传输电流；u
a
、u
b
和u
c
分别是a、b和c相的超导电缆导电层两端的电压降，由二维等效磁场模型获得；所述屏蔽层回路的电路方程如式(8)所示：式中：u
a
、u
b
和u
c
分别是屏蔽层回路a、b和c三相感应电动势；i
a
、i
b
和i
c
分别是屏蔽层回路a、b和c三相的传输电流，由二维等效磁场模型获得；i
g
是屏蔽层回路入地电流，当i
a
，i
b
和i
c
三相对称时，入地电流i
g
为0。13.根据权利要求12所述的超导电缆输电线路序参数的计算系统，其特征在于，所述通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数，具体包括：1)已知某一时刻的电源电势值，在等效电路模型中求解超导电缆导电层和屏蔽层传输电流；2)将导电层和屏蔽层传输电流代入二维等效磁场模型，通过求解式(1)和式(4)获得整体空间磁场的矢量磁位a和磁感应强度b，然后通过式(5)和式(6)得到导电层和屏蔽层电压降u和感应电流密度矢量j
s
；3)将导电层和屏蔽层的电压降u代入等效电路模型，通过迭代步骤1)和步骤2)，直到导电层电流和屏蔽层电流满足收敛条件后，进入步骤4)；4)将导电层电流和屏蔽层电流进行拟合，获得幅值和相位；5)将电压和电流的向量值相除得到序参数，所述向量值由幅值和相位构成。14.根据权利要求13所述的超导电缆输电线路序参数的计算系统，其特征在于，所述超导电缆输电线路序参数计算包括零序阻抗计算以及正序阻抗和负序阻抗计算；所述正序阻抗的计算方法如下：设定激励电源电势e
a
、e
b
和e
c
为三相对称，相位互差120
°
，通过对a、b和c相区域的电流密度进行积分，获得i
a
、i
b
和i
c
随时间变化的曲线，将i
a
、i
b
和i
c
随时间变化的曲线，截取稳态部分，按照正弦曲线进行拟合，如式(9)所示：i1＝a1sin(b1t+c1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)式中：t为时间；a1表示正序电流幅值；b1为正序角频率；c1为正序相位，拟合获得的i1为导电层回路a、b和c相的传输电流；假设超导电缆的正序阻抗为：z1＝r1+jx1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)式中，z1为正序阻抗；r1为正序电阻；x1为正序电抗；j为虚数；将式(9)和式(10)带入式(7)得到：
将拟合后的式(9)带入式(11)，则正序阻抗公式为：所述负序阻抗与正序阻抗相等；所述零序阻抗的计算方法如下：设定激励电源电势e
a
、e
b
和e
c
的幅值相等、相位相同进行仿真计算，积分获得i
a
、i
b
和i
c
随时间变化的曲线，截取稳态部分后拟合，如式(13)所示：i0＝a0sin(b0t+c0)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)式中，a0表示零序电流幅值，b0为零序角频率，c0为零序相位，拟合获得的i0即为a、b或c相的零序电流；令超导电缆的零序阻抗为：z0＝r0+jx0ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)式中，z0为零序阻抗；r0为零序电阻；x0为零序电抗；假设超导电缆两端的电压为u0，则：u0＝e
a-i0r
s
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(15)零序阻抗为：将拟合后的式(13)带入式(16)并结合式(15)，计算获得零序阻抗为：15.一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述超导电缆输电线路序参数的计算方法的步骤。16.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述超导电缆输电线路序参数的计算方法的步骤。
技术总结本发明公开了超导电缆输电线路序参数的计算方法、系统、设备及介质，获取超导电缆的结构参数和电气物性参数，所述电气物性参数包括磁导率、电导率和电阻；基于超导电缆的二维等效磁场模型和超导电缆的结构参数、磁导率及电导率计算屏蔽层感应电场和感应电流密度矢量；基于超导电缆的等效电路模型和超导电缆的电阻计算导电层和屏蔽层传输电流；通过屏蔽层感应电场和感应电流密度矢量以及导电层和屏蔽层传输电流耦合计算超导电缆输电线路序参数。本发明建立一种计算冷绝缘超导输电线路参数的分析方法，具有重要的理论价值和实践意义。具有重要的理论价值和实践意义。具有重要的理论价值和实践意义。
技术研发人员：张会明 张宏杰 陈盼盼 焦婷 李红雷 郑健 谢伟 魏本刚 丘明
受保护的技术使用者：国网上海市电力公司 国家电网有限公司
技术研发日：2022.03.30
技术公布日：2022/7/5
