CN106021676B - 一种基于转移矩阵的多回电缆稳态温升获取方法 - Google Patents

Abstract

本发明涉及一种基于转移矩阵的多回电缆稳态温升获取方法，用以获取多回电力电缆的线芯温升，包括以下步骤：1)获取多回电缆稳态温升模型的转移矩阵A，并构建多回电缆稳态温升模型：2)设定初始温升矩阵T₀并获取初始温升矩阵T₀对应的初始热流量矩阵Q₀；3)根据多回电缆稳态温升模型和初始热流量矩阵Q₀，获取下一步的温升矩阵T₁；4)判断初始温升矩阵T₀与下一步的温升矩阵T₁对应的所有元素间的最大差值是否大于收敛阈值，若是，则采用T₁取代T₀，并返回步骤2)，若否，则判定当前对应的温升矩阵为稳态温升矩阵。与现有技术相比，本发明具有计算简单、准确、计算效率高、算法先进等优点。

Description

一种基于转移矩阵的多回电缆稳态温升获取方法

技术领域

本发明涉及电力电缆运行技术领域，尤其是涉及一种基于转移矩阵的多回电缆稳态温升获取方法。

背景技术

由于电力电缆运行的特殊性，一般不可能通过直接测量去获得电力电缆线芯温度，因此技术人员提出了多种方法去计算电力电缆线芯温度，均为基于数值解和试验结果的工程公式或近似公式。如基于IEC60287标准是一种经典的计算电力电缆线芯温度的计算方法、基于数值计算的载流量求解方法。然而类似研究多针对单回电缆，对于实际运行中的多回电缆多采用数值法求解，考虑到工况的复杂性，所需计算量巨大，具体实施中效率很低。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计算简单、准确、计算效率高、算法先进的基于转移矩阵的多回电缆稳态温升获取方法。

本发明的目的可以通过以下技术方案来实现：

一种基于转移矩阵的多回电缆稳态温升获取方法，用以获取多回电力电缆的线芯温升，包括以下步骤：

1)根据热场叠加原理，将多回电缆的共同作用离散为多根电缆单独作用的组合，获取多回电缆稳态温升模型的转移矩阵A，并构建多回电缆稳态温升模型，即：

A·Q＝T

Q＝[q₁ q₂ … q_i]^T

T＝[t₁ t₂ … t_i]^T

其中，A为转移矩阵，Q为热流量矩阵，T为温升矩阵，a_i,i为第i根电缆的自发热影响系数，a_i,i-1为第i根电缆与第i-1根电缆的相互发热影响系数，且a_i,i-1＝a_i-1,i，q_i为第i根电缆的热流量，t_i为第i根电缆的温升；

2)设定初始温升矩阵T₀并获取初始温升矩阵T₀对应的初始热流量矩阵Q₀；

3)根据多回电缆稳态温升模型和初始热流量矩阵Q₀，获取下一步的温升矩阵T₁；

4)判断初始温升矩阵T₀与下一步的温升矩阵T₁对应的所有元素间的最大差值是否大于收敛阈值，若是，则采用T₁取代T₀，并返回步骤2)，若否，则判定当前对应的温升矩阵为稳态温升矩阵。

所述的步骤2)中，初始热流量矩阵Q₀内第i根电缆的初始热流量q_i0的计算式为：

其中，l_i为第i根电缆的电流量，R为第i根电缆在0℃时的直流电阻，k为直流电阻R的温度系数，k₁为考虑损耗的折算系数，t_i0为初始温升矩阵T₀中第i根电缆的初始温升。

所述的步骤4)中，收敛阈值为0.1K。

所述的步骤2)中，初始温升矩阵T₀中第i根电缆的初始温升为当前的环境温度。

所述的步骤1)中，多回电缆稳态温升模型的转移矩阵A中的参数a_i,i-1和a_i,i通过有限元、有限差分或边界元计算的方法获得：

11)由于转移矩阵A仅与周围介质的导热系数、与环境的传热系数等热学因素有关，而与电缆电流无关，因此在边界条件确定的条件下，转移矩阵A在运行中可视为不变，将A·Q＝T展开为

a_1,1*q₁+a_1,2*q₂+.......+a_1,i*q_i＝t₁

a_2,1*q₁+a_2,2*q₂+.......+a_2,i*q_i＝t₂

a_i-1,1*q₁+a_i-1,2*q₂+.......+a_i-1,i*q_i＝t_i-1

a_i,1*q₁+a_i,2*q₂+.......+a_i,i*q_i＝t_i

根据上式可知，若能获得足够多的Q矩阵与对应的T矩阵，即可通过求解上式的方程组来获得a_1,1，a_1,2………a_i-1,1，a_i,i的值，从而形成转移矩阵A，对于“足够多”的定义为设计工况正交化，且方程个数不小于未知数个数；

12)Q矩阵与对应的T矩阵的获得

借助通用的数值计算工具(如ANSYS或ANSOFT等)或专用计算软件(如COMSOL等)，建立如附图1的模型，设定某一工况的Q矩阵，即可得到对应的T矩阵，计算结果等位线如图2所示，通过此过程，即可获得一组对应的Q与T矩阵；

13)变换工况，重复步骤12)若干次(满足足够多的要求)，即可获得一定数据的新的Q矩阵与对应的T矩阵；

14)建立A·Q＝T的展开式所示方程组后求解，即可得到转移矩阵A。

与现有技术相比，本发明具有以下优点：

一、计算简单、准确：此模型的确定不依赖于电缆本身发热量或电流大小，只与电缆周围材料的热特性相关，而一般运行温度范围内该类特性可认为基本不变，这样在变换电缆电流时就无需重复有限元或其他数值计算，直接通过简单的矩阵和迭代即可获得满意的结果。

二、计算效率高：在若干次有限元计算后，截面的热学特性得以全面掌握。后续计算只需简单的矩阵和迭代，利用计算器等普通工具或笔算，即可获得满意的结果。

三、算法先进：模型本身基本与电缆损耗无关，只反映截面的热学特性，物理意义清晰，为后续的分析与改进提供了直接依据，特别是在此基础上可探索暂态温升计算的快速算法。

附图说明

图1为本发明实施例中的电缆截面示意图。

图2为某一工况的计算结果等位线图。

具体实施方式

下面结合附图和具体实施例对本发明进行详细说明。

实施例：

利用热场的叠加原理，将多回电缆的共同作用离散为多根电缆单独作用的组合，将两两之间的相互影响以转移系数来描述，进而形成集总参数构成的转移矩阵，即可实现多工况的多回电缆稳态温升的快速计算。

下面以6根独立的单芯电缆为例进行说明。

转移矩阵A：

其中a_1,1为自身发热，a_1,2为电缆1对电缆2的温升影响，根据对偶原则，a_2,1＝a_1,2，其余类似，即转移矩阵为对称的。

温升矩阵T：

热流量矩阵Q：

多回电缆稳态温升模型为：

此模型的确定不依赖于电缆本身发热量或电流大小，只与电缆周围材料的热特性相关，而一般运行温度范围内该类特性可认为基本不变，这样在变换电缆电流时就无需重复有限元或其他数值计算，直接通过简单的矩阵和迭代即可获得满意的结果。

本方法的主要步骤包括：

(1)转移矩阵的求解

1)有限元计算

有限元计算需要考虑所选计算工况的正交性与计算工况的数量，这取决于同截面电缆的回路数量，如图1所示。

图中A1-A6为电缆截面，载流量任意，边界条件1为对流散热系数为15W/m2*K，温度为30摄氏度；边界2、3、4均设为温度30摄氏度。考虑到高压电力电缆往往包含多层结构，而且有些结构层很薄。由于电缆是一个圆柱的轴对称结构，各个方向热阻相同，多层的电缆结构可以采用调和平均法进行等效，将多层电缆中导体外各层结构等效为一层等效外护层，上例中调和导热系数设为23.3W/m2*K，土壤换热系数为7.3W/m2*K。

在保证计算工况正交化的条件下，重复若干次计算，即可获得，以下结果汇总。

表1有限元计算结果

2)转移矩阵计算

由表1求解转移矩阵A：

对比矩阵中的对角元素，基本相等，也侧证了本方法的可行性与有效性。

转移矩阵的验算：

设定初始热流量为Q＝[10000；0；1391；2661；9071；3123]，通过有限元法与转移矩阵法分别求得温升矩阵如表2所示。

表2有限元与转移矩阵计算对比

计算结果	T1	T2	T3	T4	T5	T6
							有限元法计算	29.08	4.91	10.15	12.38	25.99	14.34
转移矩阵法计算	29.03	4.89	10.25	12.38	25.79	14.34
							误差	0.05	0.02	0.10	0.00	0.20	0.00

由表2可见，基于转移矩阵与基于有限元的计算结果基本一致，证明了转移矩阵法的正确性。

(2)载流量的求解

获得转移矩阵后，利用式(1)，即可获得热流量矩阵Q与温升矩阵T之间的关系。在限定Q或T后，即可得到对应的T或Q。考虑到热流量为温度的函数，一般还需要通过一定的迭代来求取。

具体步骤如下：

1)假定环境温度T₀下的热流量，取其中l_i为第i回电缆的电流量，R为第i回电缆在0℃的直流电阻，k为电阻的温度系数，k₁为考虑涡流等损耗的折算系数，其余各回电缆均如此，形成热流量矩阵Q₀。

2)利用转移矩阵，求解得到温升矩阵T₁。

3)如温升矩阵T₀与温升矩阵T₁对应的各元素间最大差异大于0.1K，利用T₁代替T₀，形成新的热流量矩阵。

4)如此重复，直至温升矩阵中对应的各元素间最大差异小于0.1K，认为计算收敛，此时的温升即为稳态温升。

表3为迭代求解过程数据。

表3迭代过程数据

Claims (5)

1.一种基于转移矩阵的多回电缆稳态温升获取方法，用以获取多回电力电缆的线芯温升，其特征在于，包括以下步骤：

1)根据热场叠加原理，将多回电缆的共同作用离散为多根电缆单独作用的组合，获取多回电缆稳态温升模型的转移矩阵A，并构建多回电缆稳态温升模型，即：

A·Q＝T

Q＝[q₁ q₂ … q_i]^T

T＝[t₁ t₂ … t_i]^T

其中，A为转移矩阵，Q为热流量矩阵，T为温升矩阵，a_i,i为第i根电缆的自发热影响系数，a_i,i-1为第i根电缆与第i-1根电缆的相互发热影响系数，且a_i,i-1＝a_i-1,i，q_i为第i根电缆的热流量，t_i为第i根电缆的温升；

2)设定初始温升矩阵T₀并获取初始温升矩阵T₀对应的初始热流量矩阵Q₀；

3)根据多回电缆稳态温升模型和初始热流量矩阵Q₀，获取下一步的温升矩阵T₁；

4)判断初始温升矩阵T₀与下一步的温升矩阵T₁对应的所有元素间的最大差值是否大于收敛阈值，若是，则采用T₁取代T₀，并返回步骤2)，若否，则判定当前对应的温升矩阵为稳态温升矩阵。

2.根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在于，所述的步骤2)中，初始热流量矩阵Q₀内第i根电缆的初始热流量q_i0的计算式为：

其中，l_i为第i根电缆的电流量，R为第i根电缆在0℃时的直流电阻，k为直流电阻R的温度系数，k₁为考虑损耗的折算系数，t_i0为初始温升矩阵T₀中第i根电缆的初始温升。

3.根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在于，所述的步骤4)中，收敛阈值为0.1K。

4.根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在于，所述的步骤2)中，初始温升矩阵T₀中第i根电缆的初始温升为当前的环境温度。

5.根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在于，所述的步骤1)中，多回电缆稳态温升模型的转移矩阵A中的参数a_i,i-1和a_i,i通过有限元、有限差分或边界元计算的方法获得。