深圳供電有限公司 洪錳 華南理工大學(xué)電力學(xué)院 盧穎

1 配網(wǎng)電纜實時導(dǎo)體溫度計算方法

1.1 配網(wǎng)電纜暫態(tài)熱路模型

結(jié)合傅立葉傳熱定律，配網(wǎng)三芯電纜可以建立暫態(tài)熱路模型[1-2]。配網(wǎng)三芯電纜結(jié)構(gòu)示意圖如圖1所示。配網(wǎng)三芯電纜暫態(tài)熱路模型如圖2所示。

圖1 配網(wǎng)三芯電纜結(jié)構(gòu)示意圖

圖2 配網(wǎng)三芯電纜暫態(tài)熱路模型

其中：Qc為單個線芯電纜導(dǎo)體及內(nèi)屏蔽層產(chǎn)生的損耗；Qd為單個線芯絕緣及絕緣屏蔽層的介質(zhì)損耗；λ1為單個線芯金屬套及屏蔽的損耗因素；λ2為單個線芯鎧裝層的損耗因素；θ1為電纜導(dǎo)體溫度；θ2為絕緣外表面及金屬屏蔽層溫度；θ3為鎧裝層溫度；θ0為電纜外護套表面溫度；T1為單個線芯絕緣層熱阻；T2為填充層及內(nèi)護套熱阻；T3為外護套熱阻；T4為電纜外部熱源至電纜表面介質(zhì)熱阻；C1'為導(dǎo)體及內(nèi)屏蔽熱容；C2'為絕緣層及絕緣屏蔽層熱容；C3'為金屬套屏蔽層熱容；C4'為填充層及內(nèi)護套熱容；C5'為鎧裝層熱容；C6'為外護套熱容。

將暫態(tài)熱路模型中同一節(jié)點上的熱流源、并聯(lián)熱容合并，可簡化為配網(wǎng)三芯電纜簡化暫態(tài)熱路模型，如圖3所示。

圖3 配網(wǎng)三芯電纜簡化暫態(tài)熱路模型

基于熱電類比理論，對圖3中的配網(wǎng)三芯電纜簡化暫態(tài)熱路模型列寫節(jié)點方程，通過配網(wǎng)三芯電纜各層結(jié)構(gòu)的等效熱參數(shù)與電纜溫度初值，得出任意時刻的電纜實時導(dǎo)體溫度。

配網(wǎng)三芯電纜簡化暫態(tài)熱路模型列寫節(jié)點方程[3]：

1.2 初值問題描述

基于熱電類比理論可以列出電纜暫態(tài)熱路模型的微分方程組，在給定電纜導(dǎo)體溫度的初值作為邊界條件的前提下，實現(xiàn)對電纜導(dǎo)體溫度的求解。而對于實際運行中的存量電纜線路，由于導(dǎo)體上高壓的存在，無法實現(xiàn)對實時導(dǎo)體溫度的直接測量。因而計算初值的設(shè)置和選取會在電纜暫態(tài)熱路方程求解時逐步累積產(chǎn)生較大的誤差，進而影響導(dǎo)體溫度精確計算。電纜負荷往往是連續(xù)波動變化，當(dāng)前時刻的導(dǎo)體溫度是電纜歷史運行狀態(tài)累積作用的結(jié)果。

2 導(dǎo)體溫度初值問題收斂性的證明

2.1 熱路模型的簡化

對于待解決的導(dǎo)體溫度初值問題，新建電纜線路可使用新型內(nèi)置光纖的電纜直接測量導(dǎo)體溫度，但對于存量電纜線路，一般只能通過電纜表皮的溫度進行導(dǎo)體溫度的推算，針對研究對象即為存量電纜線路。實際中大多數(shù)的電力電纜結(jié)構(gòu)復(fù)雜，其暫態(tài)熱路模型對應(yīng)的微分方程是高階微分方程，不便于計算。

為了簡化計算，在計及電纜各部分損耗的基礎(chǔ)上，利用調(diào)和平均法，將電纜的各層損耗歸算到導(dǎo)體損耗，并將導(dǎo)體外的各層參數(shù)歸算為同一層，則電纜的暫態(tài)熱路模型簡化為一個只有絕緣層和導(dǎo)體的一階T 型熱容熱阻傳熱模型。此時，電纜熱路模型中的微分方程組被簡化成為一個一階熱容熱阻傳熱的微分方程：

式中：P1為電纜損耗；CT為簡化后的等效熱容；TT為簡化后的等效熱阻；θ1為導(dǎo)體溫度；θ2為電纜表皮溫度；τ為時間。

為了證明經(jīng)過一段時間后，導(dǎo)體溫度在取任意初值和任意負荷的情況下，仍然能夠收斂于一個定值，需要對微分方程進行求解，得出電纜導(dǎo)體溫度關(guān)于時間的表達式，證明電纜導(dǎo)體溫升的過程中，導(dǎo)體計算溫度與真實溫度的差值隨著時間的增大而減小。相較于直接求解電纜暫態(tài)模型的微分方程組，求解該一階微分方程更為簡便。

2.2 收斂理論的推導(dǎo)證明

為了求解上述常系數(shù)微分方程，得出電纜導(dǎo)體溫度關(guān)于時間的表達式，需要將有關(guān)時間的變量和常數(shù)項進行分離。導(dǎo)體的交流電阻Rs是一個與溫度有關(guān)的函數(shù)，直接代入不能對微分方程進行參數(shù)的分離和求解，本質(zhì)上Rs與溫度呈線性關(guān)系，為了實現(xiàn)參數(shù)的分離，將Rs等效為一個只與溫度有關(guān)的線性函數(shù)：

結(jié)合電纜熱源的計算公式，可得［2］：

對上式的兩端分別進行積分，可得：

式中：C為微分方程解中待定參數(shù)。

將電纜導(dǎo)體溫度的初值作為邊界條件帶入即可得到其數(shù)值：

將式（6）代入式（5）可得：

由式（7）可知，導(dǎo)體溫度關(guān)于時間的表達式中含有導(dǎo)體溫度的初值，因此，需要分析電纜導(dǎo)體溫度的初值與指數(shù)項的乘積以及導(dǎo)體溫度的初值對導(dǎo)體溫度收斂性的影響，即分析此時θ1表達式中所含的指數(shù)項的值與零的關(guān)系。若式（8）成立，則導(dǎo)體溫度可收斂。

對式（8）整理可得：

基于傳熱學(xué)理論可以對電纜暫態(tài)熱路模型中各參數(shù)進行計算，進而對電纜導(dǎo)體溫度的收斂性進行精確的數(shù)學(xué)分析。

結(jié)合電纜導(dǎo)體交流電阻的計算公式，式（9）可以轉(zhuǎn)變?yōu)槭剑?0）[2]。基于式（10）可得滿足初值收斂條件下電纜負荷有效值應(yīng)符合的條件：

為了分析電纜負荷對于收斂速度的影響，需要利用數(shù)量級進行估算。對于不同型號的電纜，其電纜結(jié)構(gòu)參數(shù)的數(shù)量級相差不大，由估算可知TT≈10，α20（1+Ys+Yp）≈104，在20℃條件下電纜導(dǎo)體直流電阻R0≈10-5，將各參數(shù)的數(shù)量級代入式（11）進行估算，可得I1≈104。對于實際運行中的配網(wǎng)三芯電纜，其正常運行時負荷的有效值最大不超過103的數(shù)量級，則可證明式（11）恒成立，即配網(wǎng)三芯電纜導(dǎo)體溫度在任意初值的情況下均能收斂于一個定值。

3 收斂速度與導(dǎo)體溫度初值的確定方法

取實際配網(wǎng)三芯電纜運行時電流不超過103的數(shù)量級進行估算，代入下式（12）：

式（12）中：z為電纜負荷和交流電阻對該指數(shù)項的影響率，計算結(jié)果小于1%，在工程計算中可以忽略。

因此，在工程上，電纜導(dǎo)體溫度的表達式可簡化為：

為了求解實際導(dǎo)體溫度與假定導(dǎo)體溫度之間的誤差經(jīng)過多長時間才能收斂到工程上可以接受的范圍之內(nèi)，假設(shè)θ1為實際導(dǎo)體溫度，θ'1為假定導(dǎo)體溫度。定義此時真實溫度和假設(shè)溫度的誤差Δθ，則Δθ的表達式為：

由式（14）可知，此時導(dǎo)體實際溫度和假定溫度之間的誤差是以指數(shù)進行衰減，即經(jīng)過一段時間后，導(dǎo)體溫度的初值對導(dǎo)體溫度的終值的影響逐漸消失。定義工程上實際導(dǎo)體溫度和假定導(dǎo)體溫度的誤差在5%以內(nèi)可以接受，則：

基于式（15）可以推導(dǎo)得到使導(dǎo)體溫度假定值精度滿足工程實際需求的收斂時間求解公式：

由式（16）可知，導(dǎo)體溫度假定初值和真實初值間的誤差越大，收斂時間越長，雖然任意選擇導(dǎo)體溫度初值最終會收斂于其真實導(dǎo)體溫度，但是選擇合適的初值將有利于收斂過程的進行。

選定零時刻導(dǎo)體溫度初值的方法包括以下三種。

一是取零度與電纜能夠長期工作的最高溫度的中間值。

二是根據(jù)載流量大小預(yù)設(shè)，令導(dǎo)體初值溫度與電纜滿載狀態(tài)下導(dǎo)體溫度之比等于初始負荷電流與載流量之比。

三是選擇零時刻電纜負荷的穩(wěn)態(tài)導(dǎo)體計算溫度作為初值。

第三種選取方式較為接近零時刻導(dǎo)體的真實溫度。

4 結(jié)論

對導(dǎo)體動態(tài)溫度的準(zhǔn)確測算是實現(xiàn)配網(wǎng)三芯電纜動態(tài)增容的前提?，F(xiàn)有研究缺乏溫度初值選取對電纜動態(tài)熱評估結(jié)果精度影響方面的分析。針對現(xiàn)有存量配網(wǎng)三芯電纜實時熱評估是導(dǎo)體溫度初值往往無法直接通過測量獲取的問題。

基于傳熱學(xué)原理建立了配網(wǎng)三芯電纜的暫態(tài)熱路模型，并利用調(diào)和平均法將多階熱路模型簡化為一階熱路模型。通過對一階微分方程的求解，以比較數(shù)量級的方式從理論方面論證了電纜導(dǎo)體溫度初值無論如何選取，在經(jīng)歷了一段時間的迭代計算后，導(dǎo)體溫度計算結(jié)果均能收斂于真實值，為利用歷史監(jiān)測數(shù)據(jù)實現(xiàn)導(dǎo)體溫度的實時監(jiān)測提供了堅實的理論支撐。