一种基于转移矩阵的多回电缆稳态温升获取方法
【专利摘要】本发明涉及一种基于转移矩阵的多回电缆稳态温升获取方法，用以获取多回电力电缆的线芯温升，包括以下步骤：1)获取多回电缆稳态温升模型的转移矩阵A，并构建多回电缆稳态温升模型：2)设定初始温升矩阵T0并获取初始温升矩阵T0对应的初始热流量矩阵Q0；3)根据多回电缆稳态温升模型和初始热流量矩阵Q0，获取下一步的温升矩阵T1；4)判断初始温升矩阵T0与下一步的温升矩阵T1对应的所有元素间的最大差值是否大于收敛阈值，若是，则采用T1取代T0，并返回步骤2)，若否，则判定当前对应的温升矩阵为稳态温升矩阵。与现有技术相比，本发明具有计算简单、准确、计算效率高、算法先进等优点。
【专利说明】
-种基于转移矩阵的多回电缆稳态溫升获取方法
技术领域
[0001] 本发明设及电力电缆运行技术领域，尤其是设及一种基于转移矩阵的多回电缆稳 态溫升获取方法。
【背景技术】
[0002] 由于电力电缆运行的特殊性，一般不可能通过直接测量去获得电力电缆线忍溫 度，因此技术人员提出了多种方法去计算电力电缆线忍溫度，均为基于数值解和试验结果 的工程公式或近似公式。如基于IEC60287标准是一种经典的计算电力电缆线忍溫度的计算 方法、基于数值计算的载流量求解方法。然而类似研究多针对单回电缆，对于实际运行中的 多回电缆多采用数值法求解，考虑到工况的复杂性，所需计算量巨大，具体实施中效率很 低。
【发明内容】
[0003] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计算简单、准 确、计算效率高、算法先进的基于转移矩阵的多回电缆稳态溫升获取方法。
[0004] 本发明的目的可W通过W下技术方案来实现：
[0005] -种基于转移矩阵的多回电缆稳态溫升获取方法，用W获取多回电力电缆的线忍 溫升，包括W下步骤：
[0006] 1)根据热场叠加原理，将多回电缆的共同作用离散为多根电缆单独作用的组合， 获取多回电缆稳态溫升模型的转移矩阵A，并构建多回电缆稳态溫升模型，即：
[0007]
[000引 Q=[qi 化…qi]T [0009] T=[ti t2 …ti]T
[0010]其中，A为转移矩阵，Q为热流量矩阵，T为溫升矩阵，ai,I为第i根电缆的自发热影响 系数，曰1,1-1为第i根电缆与第i-1根电缆的相互发热影响系数，且= 为第i根电 缆的热流量，ti为第i根电缆的溫升；
[0011] 2)设定初始溫升矩阵To并获取初始溫升矩阵To对应的初始热流量矩阵Qo;
[0012] 3)根据多回电缆稳态溫升模型和初始热流量矩阵输，获取下一步的溫升矩阵Tl;
[0013] 4)判断初始溫升矩阵To与下一步的溫升矩阵Tl对应的所有元素间的最大差值是否 大于收敛阔值，若是，则采用Tl取代To,并返回步骤2)，若否，则判定当前对应的溫升矩阵为 稳态溫升矩阵。
[0014] 所述的步骤2)中，初始热流量矩阵Qo内第i根电缆的初始热流量qio的计算式为：
[0015]
[0016] 其中，Ii为第i根电缆的电流量，R为第i根电缆在(TC时的直流电阻，k为直流电阻R 的溫度系数，k功考虑损耗的折算系数，tio为初始溫升矩阵To中第i根电缆的初始溫升。
[0017] 所述的步骤4)中，收敛阔值为0.化。
[0018] 所述的步骤2)中，初始溫升矩阵To中第i根电缆的初始溫升为当前的环境溫度。 [0019]所述的步骤1)中，多回电缆稳态溫升模型的转移矩阵A中的参数曰11-1和ail通过有 限元、有限差分或边界元计算的方法获得：
[0020] 11)由于转移矩阵A仅与周围介质的导热系数、与环境的传热系数等热学因素有 关，而与电缆电流无关，因此在边界条件确定的条件下，转移矩阵A在运行中可视为不变，将 A ? Q = T展开为
[0021 ] ai,i*qi+ai,2*q2+.......+ai,i*qi = ti
[0022 ] B2, i*qi+a2,2*Q2+.......+曰2, i*qi = t2
[0023] ....................................
[0024] ai-i, i*qi+ai-i, 2*Q2+.......+ai-i, i*qi = ti-i
[0025] ai, i*qi+ai, 2*Q2+.......+ai, i*qi = ti
[00%]根据上式可知，若能获得足够多的Q矩阵与对应的T矩阵，即可通过求解上式的方 程组来获得ai,i，ai,2.........的值，从而形成转移矩阵A，对于"足够多"的定义为设 计工况正交化，且方程个数不小于未知数个数；
[0027] 12)Q矩阵与对应的T矩阵的获得
[002引借助通用的数值计算工具（如ANSYS或ANSOFT等）或专用计算软件（如COMSOL等）， 建立如附图1的模型，设定某一工况的Q矩阵，即可得到对应的T矩阵，计算结果等位线如图2 所示，通过此过程，即可获得一组对应的Q与T矩阵；
[0029] 13)变换工况，重复步骤12)若干次(满足足够多的要求），即可获得一定数据的新 的Q矩阵与对应的T矩阵；
[0030] 14)建立A ? Q = T的展开式所示方程组后求解，即可得到转移矩阵A。
[0031 ]与现有技术相比，本发明具有W下优点：
[0032] -、计算简单、准确:此模型的确定不依赖于电缆本身发热量或电流大小，只与电 缆周围材料的热特性相关，而一般运行溫度范围内该类特性可认为基本不变，运样在变换 电缆电流时就无需重复有限元或其他数值计算，直接通过简单的矩阵和迭代即可获得满意 的结果。
[0033] 二、计算效率高:在若干次有限元计算后，截面的热学特性得W全面掌握。后续计 算只需简单的矩阵和迭代，利用计算器等普通工具或笔算，即可获得满意的结果。
[0034] =、算法先进:模型本身基本与电缆损耗无关，只反映截面的热学特性，物理意义 清晰，为后续的分析与改进提供了直接依据，特别是在此基础上可探索暂态溫升计算的快 速算法。
【附图说明】
[0035] 图I为本发明实施例中的电缆截面示意图。
[0036] 图2为某一工况的计算结果等位线图。
【具体实施方式】
[0037] 下面结合附图和具体实施例对本发明进行详细说明。
[003引实施例；
[0039] 利用热场的叠加原理，将多回电缆的共同作用离散为多根电缆单独作用的组合， 将两两之间的相互影响W转移系数来描述，进而形成集总参数构成的转移矩阵，即可实现 多工况的多回电缆稳态溫升的快速计算。
[0040] 下面W6根独立的单忍电缆为例进行说明。
[0041 ] 转移矩阵A:
[0042]
[0043] 其中ai,i为自身发热，ai,2为电缆1对电缆2的溫升影响，根据对偶原则，曰2,1 =曰1,2， 其余类似，即转移矩阵为对称的。
[0044] 溫升矩阵T:
[0045]
[0046] g阵Q:
[0047]
[004引多回电缆稳态溫升模型为：
[0049]
[0050] 此模型的确定不依赖于电缆本身发热量或电流大小，只与电缆周围材料的热特性 相关，而一般运行溫度范围内该类特性可认为基本不变，运样在变换电缆电流时就无需重 复有限元或其他数值计算，直接通过简单的矩阵和迭代即可获得满意的结果。
[0051] 本方法的主要步骤包括：
[0052] (1)转移矩阵的求解
[0053] 1)有限元计算
[0054] 有限元计算需要考虑所选计算工况的正交性与计算工况的数量，运取决于同截面 电缆的回路数量，如图1所示。
[0055] 图中Al -A6为电缆截面，载流量任意，边界条件1为对流散热系数为15W/m2*K，溫度 为30摄氏度;边界2、3、4均设为溫度30摄氏度。考虑到高压电力电缆往往包含多层结构，而 且有些结构层很薄。由于电缆是一个圆柱的轴对称结构，各个方向热阻相同，多层的电缆结 构可W采用调和平均法进行等效，将多层电缆中导体外各层结构等效为一层等效外护层， 上例中调和导热系数设为23.3W/m2体，±壤换热系数为7.3W/m2体。
[0056] 在保证计算工况正交化的条件下，重复若干次计算，即可获得，W下结果汇总。
[0057] 表1有限元计算结果
[00581
[0
[0060] 2)转移矩阵计算
[0061 ] 由表1求解转移矩阵A:
[0062]
[0063] 对比矩阵中的对角元素，基本相等，也侧证了本方法的可行性与有效性。
[0064] 转移矩阵的验算：
[0065] 设定初始热流量为Q=[10000;0;1391;2661;9071;3123]，通过有限元法与转移矩 阵法分别求得溫升矩阵如表2所示。
[0066]表2有限元与转移矩阵计算对比 「mA7l
[00側~由表2可见，基于转移矩阵与基于有限元的计算结果基本一致，证明了转移矩阵法' 的正确性。
[0069] (2)载流量的求解
[0070] 获得转移矩阵后，利用式（1)，即可获得热流量矩阵Q与溫升矩阵T之间的关系。在 限定Q或T后，即可得到对应的T或Q。考虑到热流量为溫度的函数，一般还需要通过一定的迭 代来求取。
[0071 ] 具体步骤如下：
[0072] 1)假定环境溫度To下的热流量，巧
，其中Ii为第i回电缆的电流 量，R为第i回电缆在(TC的直流电阻，k为电阻的溫度系数，ki为考虑满流等损耗的折算系 数，其余各回电缆均如此，形成热流量矩阵Qo。
[0073] 2)利用转移矩阵，求解得到溫升矩阵Ti。
[0074] 3)如溫升矩阵To与溫升矩阵Tl对应的各元素间最大差异大于0.化，利用Tl代替To， 形成新的热流量矩阵。
[0075] 4)如此重复，直至溫升矩阵中对应的各元素间最大差异小于0.化，认为计算收敛， 此时的溫升即为稳态溫升。
[0076] 表3为迭代求解过程数据。
[0077] 表3迭代过程数据
[007引
【主权项】
1. 一种基于转移矩阵的多回电缆稳态温升获取方法，用以获取多回电力电缆的线芯温 升，其特征在于，包括以下步骤： 1) 根据热场叠加原理，将多回电缆的共同作用离散为多根电缆单独作用的组合，获取 多回电缆稳态温升模型的转移矩阵A，并构建多回电缆稳态温升模型，即： A · Q = TQ=[qi Q2 ··· qi]T T=[ti t2 ··· ti]T 其中，A为转移矩阵，Q为热流量矩阵，T为温升矩阵，a1;1为第i根电缆的自发热影响系 数，aM-i为第i根电缆与第i-Ι根电缆的相互发热影响系数，且为第i根电缆 的热流量，为第i根电缆的温升； 2) 设定初始温升矩阵To并获取初始温升矩阵To对应的初始热流量矩阵Q0; 3) 根据多回电缆稳态温升模型和初始热流量矩阵Qo，获取下一步的温升矩阵T1; 4) 判断初始温升矩阵To与下一步的温升矩阵T1对应的所有元素间的最大差值是否大于 收敛阈值，若是，则采用!^取代To,并返回步骤2)，若否，则判定当前对应的温升矩阵为稳态 温升矩阵。2. 根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在 于，所述的步骤2)中，初始热流量矩阵Qo内第i根电缆的初始热流量q l()的计算式为： 《川=尺(1 人? 其中，I i为第i根电缆的电流量，R为第i根电缆在〇 °C时的直流电阻，k为直流电阻R的温 度系数，h为考虑损耗的折算系数，tlQ为初始温升矩阵To中第i根电缆的初始温升。3. 根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在 于，所述的步骤4)中，收敛阈值为0.1 K。4. 根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在 于，所述的步骤2)中，初始温升矩阵To中第i根电缆的初始温升为当前的环境温度。5. 根据权利要求1所述的一种基于转移矩阵的多回电缆稳态温升获取方法，其特征在 于，所述的步骤1)中，多回电缆稳态温升模型的转移矩阵A中的参数B 11-JPa11通过有限元、 有限差分或边界元计算的方法获得。
【文档编号】G06F17/50GK106021676SQ201610317319
【公开日】2016年10月12日
【申请日】2016年5月13日
【发明人】傅晨钊, 司文荣, 李红雷, 姚周飞
【申请人】国网上海市电力公司, 华东电力试验研究院有限公司