王鋆晟，孟文博，文曉哲，王洪超

（吉林大學儀器科學與電氣工程學院，吉林 長春 130026）

七芯鎧裝電纜傳輸特性仿真模型研究

王鋆晟，孟文博，文曉哲，王洪超

（吉林大學儀器科學與電氣工程學院，吉林 長春 130026）

為經濟有效地獲取任意長度七芯鎧裝電纜的傳輸特性，基于多導體傳輸線理論對七芯鎧裝電纜的傳輸特性進行建模仿真分析。首先采用有限元方法計算模型參數即電纜單位長度分布參數，將趨膚效應和鋼制鎧裝鄰近效應的影響計算在內；然后以相似變化思想解模型方程即多導體傳輸線方程，獲得電纜近端及遠端傳輸特性曲線；最后以長度為3km的七芯鎧裝電纜為例，將模型仿真結果與實測數據進行對比，證明模型的有效性。結果表明：外圍纜芯的傳輸特性更加復雜多變且衰減劇烈；中間纜芯的傳輸特性衰減變化較平緩，更適合用于高速通信。

七芯鎧裝電纜；多導體傳輸線；有限元仿真；傳輸特性

0 引 言

在VSP測井及水力壓裂過程監(jiān)測等井中測試領域，對作為傳輸介質的七芯鎧裝電纜的傳輸特性有清晰的認識是進行高質量數據傳輸的前提。測試過程中所使用的七芯鎧裝電纜長達幾千米，信號在電纜上的傳輸屬于長線傳輸，受分布電阻、分布電感、分布電容等參數的影響,電纜可用帶寬很窄,長期以來一直是提升井中數據傳輸速率和質量的瓶頸。此外,電纜的傳輸特性具有非線性，非線性信道對不同形式的信號所產生的影響也不盡相同，很難給出固定參數來精確計算信號傳輸特性[1]。

自20世紀70年代末，Clayton教授[2-4]系統地提出多導體傳輸線理論，并采用這種有效的近似分析方法分析了帶狀電纜、雙絞線電纜信號傳輸過程的衰減、串擾問題，使得多導體傳輸線理論在多芯電纜傳輸特性的分析領域得到了廣泛的應用。如Mazloom等[5]針對瑞典鐵路用懸空信號電纜，進行了建模和試驗研究；Spadacini等[6]借助多導體傳輸線理論對衛(wèi)星供電線路的電磁輻射量進行預測研究；謝彥召等[7]對架空多導體傳輸線纜的電磁脈沖響應進行了計算，為防護設計提供參考。七芯鎧裝電纜相對于常用多芯電纜結構更加復雜，七根絞在一起的屏蔽電纜的外側尚有一層厚厚的鋼制鎧裝。受鐵磁性物體臨近效應的影響，電纜傳輸特性變化情況更加難以預測。目前國內主要采取測試方法獲得七芯鎧裝電纜的傳輸特性[8]。但對千米級七芯鎧裝電纜而言，測試法需針對不同的通信方案和硬件連接狀況分別進行測試，不僅會嚴重延遲研發(fā)進度，還需要大量經費支持。

為了以一種經濟有效的方法獲得任意長度的七芯鎧裝電纜的傳輸特性。參照上述多芯電纜傳輸特性分析方法，本文將七芯鎧裝電纜視作一種準TEM（橫向電磁場，即電場和磁場分量均分布在與傳播方向垂直的平面內，也即沿傳播方向沒有電場和磁場分量）結構，則各纜芯間的電壓和流過纜芯的電流能夠被唯一的定義，繼而可通過建模方法獲得電纜的傳輸特性。

圖1 多導體傳輸線系統

1 理論基礎

本文將七芯鎧裝電纜視作一種準TEM結構，也是一種多導體傳輸線結構，因而可借助多導體傳輸線理論[9]建立模型分析鎧裝電纜的傳輸特性。

分析一個由n+1個導體組成的總長為L的多導體傳輸線結構，n+1個導體互相平行且平行于Z坐標軸（如圖 1 所示）。 在近端（z=0）和遠端（z=L）處的鏈接網絡以廣義戴維南模型表示，包含集中式激勵源。

假設電磁場沿導體傳輸時是橫向電磁場（TEM）結構，以0號導體作為參考導體和電流回路，傳輸線方程為

其中V和I為n×1階向量，包含n個導體的電壓[V（z，t）]i=Vi（z，t）和 n 個導體的電流[I（z，t）]i=Ii（z，t），i代表此分量處于向量中的第i行，也代表此分量是第i個導體的數據，z表征其在z軸上的坐標，t表征時間。同樣其中R、L、C、G分別是單位長度電阻矩陣、單位長度電感矩陣、單位長度電容矩陣和單位長度電導矩陣。各矩陣組成元素如圖2所示，單位長度參數矩陣的形式如下式：

從多導體傳輸線理論可知獲得多導體傳輸線的單位長度分布參數是對多導體傳輸線分析的前提條件，解多導體傳輸線方程是獲得各個導體上的電壓、電流分布的必要過程。因而本文將首先通過有限元法計算七芯鎧裝電纜的單位長度分布參數矩陣R、L、C、G（獲取模型參數），然后通過矩陣轉換方法解多導體傳輸線頻域方程（解模型方程）以獲得千米級電纜的傳輸特性。

圖2 單位長度等效電路圖

2 模型建立

由理論分析可知，對七芯鎧裝電纜的傳輸特性建模分析主要分為兩步。

2.1 有限元法計算模型參數

一般情況下，傳輸線的單位長度參數多采用解析方式獲得。解析方式簡潔迅速，但只能解決簡單結構的問題，對結構更加復雜的多導體傳輸線（如七芯鎧裝電纜）是無能為力的[10]。

本文采用有限元方法計算七芯鎧裝電纜的單位分布參數，不受復雜結構的限制，且能充分考慮到趨膚效應和鄰近效應對單位參數的影響，獲得各對應頻率點下分布參數值。具體思路如下：

一個任意形狀的截面圖形Ω如圖3所示，其中Ω1，Ω2，…，Ωn為 n 個導體（本文中即為纜芯），Ωn+1，Ωn+2，…，Ωn+m為包圍著n個導體的m塊介電常數不同的絕緣區(qū)域（本文中為纜芯外絕緣物質）。

又有差分形式的頻域麥克斯韋方程組為

式中：σ——電導率；

E——電場強度；

B——磁通量密度；

圖3 截面的幾何構造

H——磁場強度；

D——電通量密度；

ρ——電荷密度；

ω——角頻率。

2.1.1 求解電阻、電感

結合式（3）和式（4）可得：

采用有限元方法結合諾埃曼邊界條件n×A=0（n為截面區(qū)域Ω邊界的法向量）在電纜截面區(qū)域內解方程,可解得截面區(qū)域內電流分布狀況。

從能量損耗角度[11-14]計算自電阻值rii：

式中Jz為沿Z軸方向流過電纜橫截面的電流密度矢量。

在計算出自電阻后，假設i導體中電流為Ix，j導體中電流為-Ix，則功率損耗Q與互電阻rij的關系為

自電感Lii和互電感Lij是通過計算整個區(qū)域的磁蓄能Wm獲得：

假設i導體和j導體流過的電流分別是Ii和Ij，其他導體內電流為0，則：

2.1.2 求解電容

自電容和互電容是通過在所有非導體區(qū)域內解給出狄利克雷邊界條件的拉普拉斯方程獲得電勢分布后求得。

計算i導體的自電容Cii時將其表面電位設為Vi，其他導體表面和鎧裝表面電位均為0，則：

We為在所有非導體區(qū)域Ωm=Ωn+1∪Ωn+2∪…∪Ωn+m內的電場蓄能：

互電容Cij則是假設相應的i導體和j導體的表面電位分別為Vi和Vj。

2.1.3 求解電導

電導表征導體周圍介質的導電性能，因鎧裝電纜采用的絕緣介質，電導值很小趨近于0，本文中假設其為0。

2.2 頻域方法解模型方程

按上述方法求得模型參數后，獲得傳輸特性的過程轉化為求解模型方程（多導體傳輸線方程）的過程。本文通過相似變換的方法求解多導體傳輸線方程[15]。

假定隨時間變化的傳輸線的激勵源是正弦波并已處于穩(wěn)態(tài)。因此，傳輸線上的電壓和電流也是與激勵源同頻率的正弦波。則多導體傳輸線方程的頻域形式為

其中n×n單位長度阻抗矩陣 Z=R+jωL，導納矩陣Y=G+jωC。

假設傳輸線是均勻的，即傳輸線橫截面的尺度和傳輸線周圍介質的特性沿傳輸線均是不變的（與Z無關）。將式（16）轉化為二階常微分方程

上式方程是互相耦合的。因為ZY和YZ是滿秩矩陣，即每一組的電壓 Vi（z）和電流 Ii（z）影響著其他所有組的電壓和電流。求解方法的基本思想是采用相似變換對方程解耦。

進行模量變換：

其中n×n非奇異復矩陣Tx和Ty定義為傳輸線上實際的相量電壓V和I與模電壓Vm和Im間的變量變換，則式（17）可變?yōu)?/p>

這里γ2是一個n×n對角矩陣：

若能找到變換矩陣Tx和Ty，則式（19）被成功解耦。

因為Z和Y均為對稱陣，則ZY的轉置即是YZ，有相同的特征值，因此只需找到Tx和Ty中的一個，就能解耦方程。方程（19）解耦后，參照雙導體傳輸線計算思路，可得到解為

其中，指數矩陣定義為

Vm和Im是n×1待定常數向量，它們與前向、反向橫行波有關，即：

利用式（18），可以將求得的模電壓和模電流通過變換獲得傳輸線上實際的電壓和電流：

也可轉化為僅用Im表示，即：

其中Im+、Im-在方程解偶后結合如圖1所示的終端條件，由下式求得：

其中特征阻抗矩陣 Zc=ZTyγ-1Ty-1。 電壓矩陣 VS、VL表征圖1中近端及遠端的電壓量，阻抗矩陣ZS、ZL表征圖1中近端及遠端的端接阻抗。

通過上述過程求得長度為z的七芯鎧裝電纜上的電壓、電流值，即可獲得其傳輸特性。

3 模型驗證及實例應用

按上述建模方法對一長度為3 km的七芯鎧裝電纜的傳輸特性進行建模仿真，首先使用有限元軟件計算七芯鎧裝電纜的單位長度分布參數。

使用三角形有限元對七芯鎧裝電纜的截面進行剖分，在結構尺寸越加精細的部分進行更加精細的劃分，以滿足計算精度要求（見圖4）。圖中黑色部分即為結構交界處及精細結構被進行精細剖分的結果。各部分材質屬性設置如表1所示，纜芯的材質是銅，鎧裝的材質是鋼。

圖4 電纜截面剖分效果圖

表1 電纜仿真材質屬性

如上文建模方法中所述對截面區(qū)域內的電流密度和電勢分布進行了分析（如圖5、圖6所示），可看到電流密度分布和電勢分布受趨膚效應和鄰近效應影響明顯，即鎧裝電纜的單位長度分布參數受頻率影響。通過有限元方法計算電纜單位長度分布參數即可將鄰近效應和趨膚效應的影響納入計算范圍，得到與實際狀況相符的結果。

在使用有限元法獲得電纜單位長度分布參數后，按2.2中所述方法使用Fortran語言編寫程序解模型方程，如圖1所示，在電纜末端加50Ω匹配電阻，分別以中間芯和外圍芯進行信號傳輸仿真，獲得電纜近端和遠端傳輸特性。

對一長度為3 km纏繞在電纜盤上的七芯鎧裝電纜的近端及遠端傳輸特性進行測試，在近端傳輸信號的纜芯處通過一內阻為50 Ω的信號源輸入幅值為5V的正弦激勵，其余纜芯通過50Ω的匹配電阻接地，電纜遠端同仿真模型一樣接以50Ω匹配電阻，逐個頻率點進行測試獲得測試結果。測試結果與仿真數據對比，如圖7～圖9所示。

圖5 相鄰纜芯分別通以1A，-1A（300kHz）電流時電流密度分布

可以看出電纜近端的仿真結果可以與實測數據很好的吻合，隨著頻率的增大，衰減也在逐漸增大，在100kHz附近增益約為-6dB。

經過3km的傳輸后，外圍纜芯遠端的仿真結果也能較好地與實測數據吻合。但可看出幅度增益波動劇烈，在100kHz附近波動可達5dB，因而外圍纜芯不適合用于高速傳輸方案，會產生較高的誤碼率。

圖8 3km電纜外圍芯遠端傳輸特性

圖9 3km電纜中間芯遠端傳輸特性

從圖中可看出，中間纜芯遠端的仿真結果能更好地與測試結果吻合，且通過與圖8對比可看出，中間纜芯衰減更小，幅度曲線更平滑，可以很好地解釋七芯鎧裝電纜應用中的經驗性結論：外圍纜芯的傳輸特性更加復雜多變，且衰減劇烈；比較而言，中間纜芯更適合用于高速通信。

綜上所述，通過文中所用仿真模型可以簡單有效地獲得任意長度七芯鎧裝電纜的傳輸特性，且結果中的誤差是在假設整根七芯鎧裝電纜是均勻的，即電纜任意一處橫截面的尺寸和周圍介質的特性沿電纜軸向均不變的前提下產生的，這與實際情況有一定偏差。若對精度有更高要求，可將長距離鎧裝電纜進行分段，每一段的截面尺寸和周圍介質屬性根據實際情況分別設置，上一段的輸出作為下一段的輸入進行更加精細的處理，即可根據需求得到更高的精度。

4 結束語

本文基于多導體傳輸線理論對七芯鎧裝電纜的傳輸特性進行了建模，采用有限元法計算模型參數，通過矩陣轉換方法解模型方程，獲得了電纜近端與遠端的傳輸特性。并以長度為3km的電纜為例，將模型仿真結果與實測結果進行了對比，證明了仿真模型的有效性。實際應用中鎧裝電纜沿軸向方向截面特性受多因素影響，并非均勻不變的。為進一步提高精度可將電纜截面特性沿軸向變化情況作為模型參數之一納入考慮范圍，進行分段計算，但精度的提升是以計算量的增加為代價的。

[1]張菊茜，盧濤，李群，等.一種基于OFDM技術的900 Kbit/s測井數據傳輸系統[J].測井技術，2009，33（1）：84-88.

[2]CLAYTON R P.Solution of the transmission line equations for lossy conductorsand imperfect earth[J].Proceedings of the Institution of Electrical Engineers，1975，122（2）：177-182.

[3]CLAYTON R P.A brief history of work in transmission lines for EMC applications[J].IEEE Transactions on Electromagnetic Compatibility，2007，49（2）：237-252.

[4]CLAYTON R P.A serious problem in using the jarvis model for the prediction of crosstalk[J].IEEE Electromagnetic Compatibility Magazine，2012，1（3）：46-55.

[5]MAZLOOM Z， THEETHAYI N， THOTTAPPLLIL R.Indirect lightning-induced voltages along a railway catenary-track multiconductor transmission-line system with lumped components[J].IEEE Transactions on Electromagnetic Compatibility，2011，53（2）：537-539.

[6]SPADACINI G， GRASSI F， BELLAN D， et al.Prediction of conducted emissions in satellite power buses[J].International Journal of Aerospace Engineering，2015，41（3）：439-448.

[7]謝彥召，王贊基，王群書，等.架空多導體傳輸線纜的電磁脈沖響應計算[J].清華大學學報（自然科學版），2006，46（4）：449-452.

[8]楊小鋒.測井電纜特性測試技術研究[D].大慶：東北石油大學，2014.

[9]ANTONINIG，ORLANDI A，PIGNARIS A.Review of Clayton R.PaulStudies on Multiconductor Transmission Lines[J].IEEE Transactions on Electromagnetic Compatibility，2013，55（4）：639-646.

[10]DARCHERIF A， RAIZER A,SAKELLARIS J， et al.On the use of the surface impedance concept in shielded and multiconductor cable characterization by the finite element method[J].IEEE Transactions on Magnetics，1992，28（2）：1446-1449.

[11]BONYADI-RAM S， KORDIB， BRIDGES G E.A Full-space conformal mapping for the calculation of series impedance of overhead transmission lines and underground cables[J].Electric Power Systems Research，2012（91）：95-103.

[12]杜學峰.長電纜傳輸研究[D].合肥：中國科學技術大學，2006.

[13]DEMEESTER T， GINSTE D V， DANI?L D Z.Accurate study of the electromagnetic and circuit behavior of finite conducting wedges and interconnects with arbitrary cross-sections[C]∥IEEE 19th Conference on Electrical Performance of Electronic Packaging and Systems.Washington DC：IEEE Computer Society，2011：133-136.

[14]RASOULPOOR M，MIRZAIE M，MIRIMANI S M.Losses distribution in sheathed power cables under non-sinu soidal currents using numerical method[J].Computer Aplications in Engineering Education，2016，24（5）:692-705.

[15]CLAYTON R P.多導體傳輸線分析[M].楊曉憲，鄭濤，譯.2版.北京：中國電力出版社，2013:176-191.

（編輯：李妮）

Simulation research on seven-core armored cable’s transmission characteristics

WANG Junsheng，MENG Wenbo，WEN Xiaozhe，WANG Hongchao
（College of Instrumentation ＆ Electrical Engineering，Jilin University，Changchun 130026，China）

In order to obtain the transmission characteristics of seven-core armored cable of arbitrary length in a cost-effective manner，based on the theory of multiconductor transmission lines，the transmission characteristics of seven-core armored cable were simulated and analyzed.Firstly，the finite elementmethod wasused to calculate the cable unitlength distribution parameters，and the effect of the skin effect and proximity effect caused by steel armor were calculated.Then，the transmission equation of the multiconductor was solved by the similarity transformation idea.Finally，the simulation results of seven-core armored cable with length of 3km were compared with the measured data，which proved the effectiveness of the model.The results show that the transmission characteristics of the outer cable core are more complicated and decay drastically.The transmission characteristics of the inter cable core change more slowly and more suitable for high-speed transmission.

seven-core armored cable；multiconductor transmission line；finite element simulation；transmission characteristics

A

1674-5124（2017）04-0118-07

10.11857/j.issn.1674-5124.2017.04.025

2016-09-03；

2016-10-20

國家自然科學基金項目（41074074）

王鋆晟（1989-），男，河南平頂山市人，碩士研究生，專業(yè)方向為井下通信技術。