![点击分享此内容可以赚币 分享](/master/images/share_but.png)
双导体有损频变均匀传输线的电磁暂态时域仿真模型研究_刘刚.pdf
《双导体有损频变均匀传输线的电磁暂态时域仿真模型研究_刘刚.pdf》由会员分享,可在线阅读,更多相关《双导体有损频变均匀传输线的电磁暂态时域仿真模型研究_刘刚.pdf(11页珍藏版)》请在咨信网上搜索。
1、第 47 卷 第 2 期 电 网 技 术 Vol.47 No.2 2023 年 2 月 Power System Technology Feb.2023 文章编号:1000-3673(2023)02-0823-10 中图分类号:TM 721 文献标志码:A 学科代码:47040 双导体有损频变均匀传输线的电磁暂态 时域仿真模型研究 刘刚1,温晓芳1,郝世缘1,纪锋2(1华北电力大学电力工程系,河北省 保定市 071003;2先进输电技术国家重点实验室(全球能源互联网研究院有限公司),北京市 昌平区 102209)Study on Electromagnetic Transient Time-d
2、omain Simulation Model of Lossy Frequency-dependent Uniform Transmission Line With Two-conductor LIU Gang1,WEN Xiaofang1,HAO Shiyuan1,JI Feng2(1.Department of Electrical Engineering,North China Electric Power University,Baoding 071003,Hebei Province,China;2.State Key Laboratory of Advanced Power Tra
3、nsmission Technology(Global Energy Interconnection Research Institute Co.,Ltd.,)Changping District,Beijing 102209,China)ABSTRACT:The characteristic impedance and the transfer function are two important parameters in the electromagnetic transient time-domain simulation model of the multi-conductor tr
4、ansmission line.In this paper,taking the lossy frequency-dependent uniform double-conductor transmission line as the research object,the coupling amount of the transmission line is decoupled by using the phase-mode transformation.Then,in the mode domain,the time-domain equivalent circuit of the char
5、acteristic impedance is obtained by using the Pade approximation and the Foster circuit model.Simultaneously,the approximate function of the transfer function is got by using the Laplace inverse transformation and the Prony approximation.Based on the equivalent circuit and the approximation function
6、,the electromagnetic transient time-domain simulation model of the double-conductor transmission line is established by adopting the MATLAB programming.Finally,this time-domain simulation model is used to achieve the voltage responses at the beginning and the terminal ends of the transmission line u
7、nder the step excitation.The simulation results are compared with those using the PSCAD/EMTDC software to verify the time-domain model in this paper.KEY WORDS:Bergeron model;Pade approximation;Prony approximation;frequency-dependent parameters;electromagnetic transient 摘要:特征阻抗和传输函数是多导体传输线电磁暂态时域仿 基金项
8、目:先进输电技术国家重点实验室开放基金项目(GEIRI SKL2018001)。Project Supported by State Key Laboratory of Advanced Power Transmission Technology(GEIRISKL2018001)真模型中的 2 个重要参数。该文以有损频变均匀双导体传输线为研究对象,采用相模变换实现传输线耦合量的解耦,在模态域利用 Pade 近似和 Foster 电路模型获得特征阻抗的时域等效电路。同时利用拉普拉斯反变换和 Prony 近似得到传输函数的近似函数。该文基于等效电路和近似函数通过MATLAB 编程建立了双导体传输线
9、的电磁暂态时域仿真模型。最后,用该时域仿真模型获得了阶跃激励下传输线首末端的电压响应,并与 PSCAD/EMTDC 软件仿真结果进行对比分析,验证了所提时域模型的有效性。关键词:Bergeron 模型;Pade 近似;Prony 近似;频变参数;电磁暂态 DOI:10.13335/j.1000-3673.pst.2021.2318 0 引言 当电力系统发生短路、过电压或其它非正常运行工况时,传输线上的电压电流会在短时间内快速变化,可能会对电气设备甚至整个电力系统造成危害1,为了避免此类危害的发生,需要准确分析传输线和电气设备上的电压电流动态变化过程,因此建立准确的传输线暂态时域模型至关重要2-
10、3。最初提出的传输线时域模型为 Bergeron 模型,但 Bergeron 模型未考虑损耗分布特性和线路参数频变特性的影响4。而实际传输线的距离很长,并且线路参数是随频率变化的5-7,因此,很多学者开始研究频变传输线电磁暂态模型。Budner 提出了时域瞬态解的频率相关线路模型,他在导纳线路模型中使用了加权函数的概念8。但是,加权函数高度振荡且较难精准拟合。之后,学者们开始在算法改进和拟合精度等方面展开研究。Snelson 提出了行824 刘刚等:双导体有损频变均匀传输线的电磁暂态时域仿真模型研究 Vol.47 No.2 波法9,引入一个常数将电流和电压在时域中联系起来(其表达式与 Berg
11、eron 模型中对波过程解释的表达式类似),然后转换到频域中求解。Meyer 和Dommel 在 Snelson 的基础上,提出了前、反行波权函数法10,使计算结果更加准确,但是该方法含有多个卷积积分,计算复杂。为解决这一问题,Semlyen 通过指数函数来拟合线路的冲击响应和阶跃响应,并用插值法把复杂的卷积运算转化为递归公式的求解11。J.Marti 基于 Semlyen 思想,建立了更加有效的 J.Marti 模型,该模型采用渐进拟合法在频域中拟合特征阻抗和传输函数,将模拟滤波技术应用于频变传输线的暂态计算中12,该思想受到广泛应用,但由于拟合的有理式零极点被局限于实数域,在一定程度上降低
12、了拟合的精确度。后来,Gustavsen 提出了矢量匹配法,该方法减少了计算量并使拟合精度进一步提高13,其不足为不同的初始迭代极点会导致不同的有理式;因此,有学者将矢量匹配法和遗传算法结合起来完成拟合14。这种方法精度较高且计算速度较快,但过程稍显复杂。针对以上研究现状,为简化特征阻抗和传输函数的近似方法,本文将分别用 Pade 近似和 Prony近似实现特征阻抗和传输函数的拟合,该方法只需在已知线路参数的前提下,采用 Pade 近似得到特征阻抗的 Foster 等效电路,采用 Prony 近似得到传输函数的近似函数,然后将其应用于 Bergeron 模型,最终通过 MATLAB 编程获得适
13、用于时域仿真的有损频变均匀传输线电磁暂态模型。本文所提方法无需迭代,更简单易懂,易于程序实现。为了验证方法的有效性,论文分别用所提方法建立的时域传输线模型和 PSCAD/EMTDC 对双导体传输线的阶跃响应进行仿真,并对比分析了二者的结果。1 相模变换 为了分析方便,本文以双导体传输线为例,其电路图如图 1 所示。0t dyRdyRdyuLRLR 图 1 双导体传输线电路图 Fig.1 Circuit of two-conductor transmission line 图中:udy表示电压源;Rdy表示电源内阻;l表示传输线长度;RL表示负载电阻;L、C、R、G分别表示传输线的串联电感分布参
14、数矩阵、并联电容分布参数矩阵、串联电阻分布参数矩阵、并联电导分布参数矩阵。电路的时域电报方程为 abababxtuiRiL(1)abababxtiuGuC(2)式中:aaababbbuiui ui,。电力系统中,传输线一般为循环对称结构,因此有:zzhzhzzzhzhz0=00=0RLLRLLGCCGCCRLGC,(3)式(1)(2)变换到复频域,可得如下方程:()ababsxURL I(4)()ababsxIGC U(5)式中:aaababbbUIUIUI,。从式(3)知,R,L,G,C 均为对称矩阵。导线之间的电磁联系导致了互电感和互电容的产生,使 L 矩阵和 C 矩阵中存在耦合元素,为方
15、便计算,需采用相模变换实现解耦,相模变换示意图如图 2所示。图 2 相模变换示意图 Fig.2 Schematic diagram of phase-mode transformation 以 ab 代表 2 个导体,代表 2 个模量。从两相自然坐标系 ab 到静止坐标系 的 Clarke变换为=ab T,111=112T(6)矩阵 T 的 2 个列向量是一组标准正交基,其转置矩阵即为自身的逆矩阵,故对应的 Clarke 反变换为 1=ab T,1T111=112TT(7)由于该变换与频率无关,所以可将式(7)代入式(4)(5),整理得到:1()sxUT RL TI(8)1()sxIT GC
16、T U(9)在循环对称结构下,T(R+sL)T1和 T(G+sC)T1第 47 卷 第 2 期 电 网 技 术 825 均是对角矩阵15-16。这样就利用 Clarke 变换实现了解耦,从而可以应用单导体传输线理论解决双导体问题。2 复频域 Bergeron 等效电路 复频域下单导体传输线首末端的电压和电流可以用传输矩阵17表示为()()()()()()()()(k,)(k,)eeee()(m,)22(m,)1eeee()22s ls ls ls ls ls ls ls lUsIsZ sUsIsZ s(10)式中:k 代表线路首端;m 代表线路末端;U(k,s),I(k,s)分别表示线路首端的
17、电压和电流,传输函数G(s)=e(s)l;传播常数()()()sRsL GsC;特征阻抗()()()RsLZ sGsC;U(m,s)、I(m,s)分别表示线路末端的电压和电流。用式第一行减去式(10)第二行乘以 Z(s)并整理得到线路首端电压:()(k,)(k,)()(m,)(m,)e()()s lUsIs Z sUsIsZ sZ s(11)类似地,线路末端电压 ()(m,)(m,)()(k,)(k,)e()()s lUsIs Z sUsIsZ sZ s(12)根据式(11)(12)可得复频域下 Bergeron 等效电路如图 3 所示。图 3 复频域 Bergeron 等效电路 Fig.3
18、Bergeron equivalent circuit in complex frequency domain 图 3 中,Iks和 Ims为线路首末端的历史电流源,如式(13)(14)所示。()k(m,)(m,)e()s lsUsIIsZ s (13)()m(k,)(k,)e()s lsUsIIsZ s (14)从图中可知,求得单导体传输线首末端电压电流的前提是特征阻抗和传输函数的准确拟合。3 基于 Pade 近似的特征阻抗拟合 本文对于特征阻抗拟采用 Pade 近似,Pade 近似是一种有理分式逼近法,以尽量快的速度与泰勒级数展开式相匹配。如对1212()1xx做近似,Pade 近似采用
19、5 阶的有理分式得到的结果与标准值的精确度仅差在 108数量级,而泰勒级数展开式需要取11 项才能达到相同精度。此外,Pade 近似往往比截断的泰勒级数准确,而且当泰勒级数不收敛时,Pade 近似往往仍可行,所以多应用于基于计算机的科学计算中。其基本思路为将原函数近似为有理多项式相除的形式18-19,本文中Pade近似函数用Zd(s)表示:1121d112()()1nnnnnnny sy sy syZ sZsx sx sx s(15)式(15)两端同乘以分母并整理得 1121112()()()()nnnnnnny sy sy syZ s x sZ s x sZ s x sZ s(16)式中有
20、2n+1 个待定系数,选择 2n+1 个频点构成以下方程组,求解即可。MBZ (17)其中,111111212121212121=1()()1()()nnnnnnnnnnssZ s sZ s sssZ ssZ ssM 1+11 nnnyyyxxB T121=()()nZ sZ sZ 至此,可求出式(16)中的 x、y 值,也就得到了特征阻抗的有理近似函数,之后,将其展开为部分分式和的形式,如式(18)所示。d11211112()=1nnnjnnnnjjnZsry sy sy sykspx sx sx s(18)首先将式(18)中的 k 等效为一个电阻,然后将每一项分式都等效为一个 RC 并联模
21、块。这样就得到了特征阻抗的Foster模型(即一系列RC并联模块的串联),如图 4 所示,其中 Cj和 Rj由有理近似函数的极点和留数确定。单个 RjCj并联模块如图 5 所示。根据图 5,求出复频域下该电路的阻抗,并将其与部分分式建立等式关系,如式(19)。826 刘刚等:双导体有损频变均匀传输线的电磁暂态时域仿真模型研究 Vol.47 No.2 图 4 Foster 等效模型 Fig.4 Foster equivalent model 图 5 RjCj并联模块 Fig.5 RjCj parallel module 1=11jjjjjjjjjjRsCCrZspRssCC R(19)为方便计算
22、,去掉极点和留数的虚部,但拟合效果仍比较理想(从图 9 和图 10 中可知)。经整理得Foster 模型中的电容和电阻值如式(20)所示。1/jjjjjCrRrp (20)4 基于 Prony 近似的传输函数拟合 如果用 Pade 近似来拟合传输函数,会得到一些虚部不可忽略的极点和留数,在之后进行历史电流源的计算时难以处理,所以本文换一种思路,采用 Prony 近似拟合传输函数。Prony 近似采用指数函数的线性组合来描述等间距采样数据,该算法通过求解常系数线性方程组来得到近似函数,避免了非线性方程组的求解,广泛应用于电力系统响应信号分析、电磁暂态过程以及各种控制系统的协调设计研究中20,其基
23、本思路如下:设原函数为 g(t),近似函数为 gd(t),将 gd(t)由2n 个等间距采样点表示21-22 d1()ejnp tjjg tr(21)可知式(21)中有 2n 个参数,分别为留数 rj和极点 pj,每个采样点对应原函数中一个时刻的值:()0,1,21ig iTgin,(22)式中 T 为采样周期,第一个采样点为原点。令 e 1,jp Tjzjn,(23)结合式,得:1=0,1,21niijjjgr zin,(24)此时引入变量 i,有:01()1nniijnijzzz,(25)根据式(24)和(25),可得:00110()=()0 0,1,1nnnk ik iijjiiijnn
24、kijjijjigr zr zzkn,(26)进而:10 0,1,1nk iik niggkn,(27)由式(27)求出 i;将 i代入式(25),求出 zj;根据式(23),可得 pj;将已求出的值代入式(24),求得rj,如此可得近似函数 gd(t)。为简化计算,先对传播常数做适当变形,具体推导过程见附录 A。2222()()1()()24RsL GsCRGRCGLLC sLCL C(28)令|2RCGLLC,1()2RGcLC,传播常数可表示为 22()()sLCsc (29)令llLCv,波速1=vLC,则传输函数为 22()()es cG s (30)表 1 为拉普拉斯变换对23。表
25、 1 拉普拉斯变换对 Table 1 Laplace transform pairs 时域形式 复频域形式 22122e()e()()ctcttIttt 1e()()*ejnp tcjjttr eiptir 22()es c 1(e)enjcsjjrsp iirsp 第 47 卷 第 2 期 电 网 技 术 827 由表 1 可知传输函数的拉普拉斯反变换为 22122()e()e()()ctctg ttIttt (31)令 22122()e()()ctI tIttt(32)式(31)中含有贝塞尔函数221()It,故先对 I(t)作 Prony 近似,得到其近似函数:d1()ejnp tjjI
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 导体 有损 均匀 传输线 电磁 时域 仿真 模型 研究
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。