本发明涉及一种基于分布参数模型的超特高压系统电晕损耗计算方法，属于电力系统线损技术领域。
背景技术：
在现代电力系统中，超/特高压输电线路已成为我国的主力电网。对于一般电压等级的输电线路，可以用集中参数等值电路。但是超/特高压输电线路传输距离长达几百公里甚至上千公里，研究此线路的传输特性时，就需要采用分布参数模型，此时，如再使用集中参数模型就会有很大的偏差。
由于超/特高压输电线路的电压等级高，线路表面易产生放电现象，因而超/特高压输电线路的电晕损耗是不可忽视却十分复杂的现象。据大量实测研究，500kv及以上线路电晕损耗对线路实际运行电压非常敏感，且占线路总损耗相当大的比重，因而在计算系统线路总损耗时电晕损耗不应忽略。
技术实现要素：
目的：为了克服现有计算方法中存在的不足，本发明提供一种基于分布参数模型的超特高压系统电晕损耗计算方法。
技术方案：为解决上述技术问题，本发明采用的技术方案为：
一种基于分布参数模型的超特高压系统电晕损耗计算方法，具体包括如下步骤：
s01：基于长距离交流输电线路的分布参数模型推导出超特高压长距离交流输电线路的π形等值电路参数计算公式；
s02：通过精确考虑超特高压系统长线路上的电压分布推导出线路单位长度上的年平均电晕损耗公式；
s03：基于电晕损耗功率方程进行牛顿-拉夫逊迭代，进而求出每条线路上的电晕损耗和系统的电晕损耗。
在步骤s01中推导出长线路等值电路各参数的计算公式：
基于超特高压输电线路的特性阻抗：
其中：
线路传播常数：
其中：
精确考虑输电线路的分布参数特性：
其中：l为输电线路的长度；
令u＝αl，v＝βl，利用双曲函数
经过推导和化简，令
最终可得超特高压输电线路π形等值电路的参数计算公式：
在步骤s02中推导出线路单位长度上的年平均电晕损耗公式：
在超特高压输电线路的长线模型下，已知线路末端电压、电流时，沿线路任意点的电压、电流计算公式：
其中：
则超特高压系统中长线路每一点的电压表达式：
其中：
单位长度上好天、雪天和雨天的电晕损耗计算公式为：
其中：pf、ps、pr分别表示单位长度上好天、雪天和雨天的电晕损耗(kw/km)；req表示分裂导线等效半径(cm)；em为导线表面最大场强(kv/cm)；e0为以皮克公式计算的临界电场强度(kv/cm)；
中相的表面场强为：
边相的表面场强为：
其中：n为分裂导线根数；r’为导线半径(cm)；d为分裂间距(cm)；dm为几何均距(cm)；u为实际运行电压(kv)；
k2＝1.03k1/1.1；
皮克公式表示为：
其中：m1为导线表面粗糙系数；m2为气象系数；δ为空气的相对密度；
如果设一年中好天为tf，雪天为ts，雨天为tr,同时在线路中均匀产生t+1个点，从首端到末端的每点电压依次为u1、u2、……、ut+1，则线路的电晕损耗即为t个小线段电晕损耗的总和，小线段实际运行电压表示为该段首端和末端的平均电压，线路中n、m两点间的小线段单位长度上电晕损耗的年平均值为：
其中，un、um分别为线路中n、m两点的电压；
那么线路单位长度上电晕损耗的年平均值为：
在步骤s03中求出每条线路上的电晕损耗和系统的电晕损耗：
为简化计算起见，将线路的电晕损耗分为相等的两部分悬挂在输电线路的两端，认为是“有功负荷”，这样推导出的对于一个n节点系统来说，节点i的功率方程为：
将
其中：pi为pq和pv节点的注入有功功率；qi为pq节点的注入无功功率；gij、bij分别为支路电导和电纳；θij为节点i和j之间的电压相角差；lij为线路i-j的长度；ui、uj分别为节点i、j的电压；
牛拉法迭代后可得系统每条线路单位长度年平均电晕损耗值pcij，进而求得每条线路上的电晕损耗pcijlij，最后累加可得系统的电晕损耗。
有益效果：本发明提供的一种基于分布参数模型的超特高压系统电晕损耗计算方法，此计算方法精确考虑了线路上的每一点电压，因而结果更加准确，对超/特高压系统的设计以及分析具有重大的意义。
附图说明
图1是本发明的算例八节点系统图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示，八节点系统其中平衡节点5，pv节点6，其余都为pq节点；系统共有5条线路，其中l1、l2、l3、l4的线长为400km，l5的线长为800km；tf＝7471h，ts＝56h，tr＝1233h，分段数t＝20，空气相对密度δ＝1，导线粗糙系数m1＝0.9，好天时m2＝1.0，雪天时m2＝0.997，雨天时m2＝0.993，500kv线路型号为4*lgj-300，导线半径r’＝1.213cm，分裂间距d＝45cm，三相水平排列，相间距离为13m。
在步骤s01中求出长线路等值电路各参数的值：
线路l1、l2、l3、l4的r＝9.88572ω，x＝109.126ω，b＝0.00321s,线路l5的r＝16.279ω，x＝199.285ω，b＝0.00673s。
在步骤s02中求出线路单位长度上电晕损耗的年平均值：
在步骤s03中求出每条线路上的电晕损耗和系统的电晕损耗：
为简化计算起见，将线路的电晕损耗分为相等的两部分悬挂在输电线路的两端，认为是“有功负荷”，这样推导出的对于一个n节点系统来说，节点i的功率方程为：
将
其中：pi为pq和pv节点的注入有功功率；qi为pq节点的注入无功功率；lij为线路i-j的长度。
牛拉法迭代后可得系统每条线路单位长度年平均电晕损耗值pcij，进而求得每条线路上的电晕损耗pcijlij：线路l1、l2上的电晕损耗均为3.55693mw，线路l3、l4上的电晕损耗均为3.90702mw，线路l5上的电晕损耗为7.94192mw；最后可得系统的电晕损耗为22.8698mw。以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。