一种直流输电线路电场强度及无线电干扰计算的方法技术

本发明专利技术公开了一种直流输电线路电场强度及无线电干扰计算的方法，解决了现有技术和传统方法对于输电线路导线断面具有起伏地势无法直接设置导线镜像电荷的问题；另外，解决了直流输电线路采用高斯牛顿法等线性化方程会遗留原方程不适定的问题。与传统局放方法对比，本发明专利技术提高了试验效率，经过试验验证，具有实用性和可行性。

【技术实现步骤摘要】
一种直流输电线路电场强度及无线电干扰计算的方法
本专利技术属于电力
，尤其涉及在高压、特高压领域，直流输电线路电场强度及无线电干扰计算的方法。
技术介绍
高压或特高压直流输电技术广泛用于西电东送主通道的输电工程，然而随着电压等级的提高，其电磁环境影响越来越受到公众的关注和线路设计的重点考虑。直流线路下方的合成电场、离子电流、直流磁场，以及电晕产生的无线电干扰和可听噪声会对沿线居明的生产和生活产生一定的影响。研究表明，无线电干扰水平是直流输电线路环评的重要电磁环境指标之一，其主要危害是造成输电线路周围的无线电广播、电视等无线电信号干扰。文献《负直流下绞线电晕起始电压分析》、《高海拔条件下钢芯铝绞线的正直流电晕起始电压分析》基于气体放电理论分析了正负直流绞线电晕起始电压，文献《高海拔±800kV直流输电线路电磁环境测量》进行了特高压直流线路无线电干扰特性研究，上述理论分析和运行试验的实际测量均选定输电线路导线垂直断面为平面，结合场域边界条件，利用马克特-门格尔法、逐步镜像法进行无线电干扰分析和计算。然而实际的输电线路导线垂直断面经常处于不平坦地形。基于上述文献中的计算方法不能直接用于复杂导线垂直断面地势下直流线路电场强度和无线电干扰计算。文献《复杂地形情况下高压交流输电线路电磁环境特性分析》、《复杂地势下超高压输电线路的工频电场》基于优化的模拟电荷法开展了复杂地形情况下高压交流输电线路电磁环境特性分析，但该方法是凭经验设置模拟电荷的位置。文献《考虑气象条件的输电导线工频电场计算新方法》、《交叉跨越输电导线三维工频电磁场计算》利用遗传算法分别进行电场逆运算的输电导线弧垂和特高压绝缘子均压环优化，然而传统方法只针对平坦的导线断面，利用设置镜像电荷等效大地表面的感应电荷，针对输电线路导线断面具有起伏地势的情况，无法直接设置导线镜像电荷。
技术实现思路
为了解决上述问题，本专利技术提供的一种直流输电线路电场强度及无线电干扰计算的方法，具体步骤如下：S1：确定模拟电荷空间位置模拟电荷的位置满足式(1)的约束条件；为提高分辨率将分裂子导线和避雷线导线采用10位二进制编码，镜像电荷采用10位二进制编码；式(1)中(xi,yi)和(xj,yj)为分裂子导线和避雷线导线内的模拟电荷坐标，(xk,yk)为大地表面以下的镜像电荷坐标，R1和R2为分裂子导线和避雷线导线的半径，A和N分别为避雷线导线内模拟电荷的数量和总的模拟电荷数量；f(xk,yk)为描述大地表面连续曲线的函数；S2：模拟电荷位置初始种群P(t)的产生分裂子导线和避雷线导线的模拟电荷分别均分布于直径1/4的圆环，各镜像电荷分别位于相对地面垂直的等距离镜像位置；将式(2)作为适应度函数评价个体的优劣：式(2)中M为分裂子导线表面的匹配点数量，Q和分别未知模拟电荷的电量向量和已知的匹配点的电位向量，为分裂子导线、避雷线导线或大地表面第i个匹配点的电位；S3：判断收敛性若|fmax-fmin|<ε，则输出寻优结果，停止迭代，转步骤S4；否则转入S5；其中fmax,fmin分别为当前种群中最优和最差个体的函数值，ε为给定准确度；S4：正则化在式(2)中，为实现对非线性算子F(Q)的近似，利用泰勒公式将F(Q+δQ)在Q处的泰勒展开，当||δQ||足够小时，采用一次近似且忽略其余的高阶小量，得到式(3)；令Q*＝Q+δQ是方程的准确解，则在接近于Q*的Q处可以由式(3)得到如下线性算子方程式(4)：结合线性化处理的式(4)，基于全局正则化且通过优化问题得出一个最佳的δQ如式(5)：式(5)中，W为线性算子，η为信赖域大小，μ(η)为罚函数；对于(5)式所示的优化问题，利用f为极小解的必要条件是f的梯度为零的点，得到式(6)的线性方程：Qk+1＝Qk+δQ式(6)中，为雅可比矩阵的转置；式(2)和式(5)的非线性最小二乘问题的求解归结为其法方程的求解如式(6)；S5：选择、交叉和变异在S3中，若不满足|fmax-fmin|<ε，需进行种群的选择、交叉和变异：其中，选择规则是按照适应度大小对个体进行排序，从中选取i个适应度最大个体，按照下式(7)进行选择概率Ps(xi)：交叉规则采用两位交叉方式，自适应交叉概率Pc为式(8)：式(8)中，favg为群体的平均适应度，f'为两交叉个体较大的适应值，fmax为当前种群中最优个体的函数值；变异规则是将每个个体的每两位二进制编码随机产生一个在(0，1)之间的数，若大于个体的变异率，则该编码由1变成0，或者由0变成1，否则该编码不变异，Pm变异概率为式(9)：式(9)中，favg为群体的平均适应度，f为每个个体的每两位二进制编码随机产生一个在(0，1)之间的数，fmax为当前种群中最优个体的函数值；然后，保留较优个体，用子代种群的适应值较大的个体替代父代种群适应值较小的相同数目个体，即P(t)＝P(t+1)，提高获得最优个体的概率，将当前保留的新种群作为步骤S2的初始种群重新计算；S6：判断式(6)的是则转S12；否则转S7；S7：计算δQ,fk,fk+1,f1k+1通过式(6)的信赖域方法确定适合的δQ：式中Δf是非线性目标值的改变，Δf1是线性化目标值的改变，同时fk+1,fk分别是k+1和k次迭代非线性目标值，f1k+1是是k+1次迭代线性化目标值，具体计算方法如式(11)：式(11)中，fk+1,fk分别是k+1和k次迭代非线性目标值，f1k+1是是k+1次迭代线性化目标值，W为线性算子，Q和分别为未知模拟电荷的电量向量和已知的匹配点的电位向量；S8：判断fk>fk+1，是则转S9；否则将信赖域大小η变为原来的0.9倍，即0.9η，转入S7迭代计算；S9：判断式中，Δf是非线性目标值的改变，Δf1是线性化目标值的改变，τ为信赖域控制参数；若转S10；否则，将信赖域大小η变为原来的1.1倍，即1.1η，转入S7重新迭代计算；S10：判断是否大于最大迭代次数实际工程计算中，为避免大量运算，需要设定最大迭代次数，最大迭代次数小于等于1000次；若大于最大迭代次数，转S5；否则，转S11；S11：按式(12)进行fk+1的k+1迭代非线性目标值计算，Qk+1的k+1迭代，然后转S4；其中，fk+1,fk分别是k+1和k次迭代非线性目标值，Qk+1，Qk分别是k+1和k次迭代未知模拟电荷的电量向量；S12：计算输电线路表面电场强度式(13)中，E为电场强度的大小，其分量分别为Ex和Ey，Li为模拟电荷至计算点的距离；因此，基于式(13)计算直流双极分裂子导线表面各点场强，求解分裂子导线最大表面场强Emax；S13：计算任意地形的无线电干扰水平；结合分裂子导线半径和计算点距离正极的距离等参数，直流线路无线电干扰水平计算的CISPR推荐方法如式(14)：式中G为特高压直流无线电干扰水平，Emax为分裂子导线最大表面场强，r为子分裂子导线半径，n为分裂分裂子导线数，R为参考点至最近分裂子导线的空间距离；ΔEf为干扰频率修正项，ΔEf＝5[1～2(log10f)2]，f为干扰频率，现场可忽略ΔEf；ΔEw为气象修正项，好天气取0dBV/m，坏天气取3dBV/m；S14：计算可听噪声；可听噪声PdB采用美国电力科学研究院(本文档来自技高网...

【技术保护点】
一种直流输电线路电场强度及无线电干扰计算的方法，其特征在于，S1：确定模拟电荷空间位置模拟电荷的位置满足式(1)的约束条件；为提高分辨率将导线和避雷线采用10位二进制编码，镜像电荷采用10位二进制编码；(xi-x0i)2+(yi-y0i)2<R12(i=1...M)(xj-x0j)2+(yj-y0j)2<R22(j=1...A)f(xk,yk)<0,(k=1...N-A)---(1)]]>式(1)中(xi,yi)和(xj,yj)为分裂子导线和避雷线内的模拟电荷坐标，(xk,yk)为大地表面以下的镜像模拟电荷坐标，R1和R2为分裂子导线和避雷线的半径，A和N分别为避雷线内模拟电荷的数量和总的模拟电荷数量；f(xk,yk)为描述大地表面连续曲线的函数；S2：模拟电荷位置初始种群P(t)的产生导线和避雷线导线的模拟电荷分别均分布于直径1/4的圆环，各镜像电荷分别位于相对垂直地面的等距离镜像位置；将式(2)作为适应度函数评价个体的优劣：式(2)中M为导线表面的匹配点数量，Q和分别未知模拟电荷的电量向量和已知的匹配点的电位向量，为导线、避雷线或大地表面第i个匹配点的电位；S3：判断收敛性若|fmax‑fmin|＜ε，则输出寻优结果，停止迭代，转步骤S4；否则转入S5；其中fmax,fmin分别为当前种群中最优和最差个体的函数值，ε为给定准确度；S4：正则化在式(2)中，为实现对非线性算子F(Q)的近似，利用泰勒公式将F(Q+δQ)在Q处的泰勒展开，当||δQ||足够小时，采用一次近似且忽略其余的高阶小量，得到式(3)；F(Q+δQ)=F(Q)+∂F∂QδQ+12∂2F∂Q2(δQ)2+R(Q,δQ)≈F(Q)+∂F∂QδQ---(3)]]>令Q*＝Q+δQ是方程的准确解，则在接近于Q*的Q处可以由式(3)得到如下线性算子方程式(4)：结合线性化处理的式(4)，基于全局正则化且通过优化问题在该区域内得出一个最佳的δQ如式(5)：式(5)中，W为线性算子，η为信赖域大小，μ(η)为罚函数；对于(5)式所示的优化问题，利用f为极小解的必要条件是f的梯度为零的点，得到式(6)的线性方程：Qk+1＝Qk+δQ式(6)中，为雅可比矩阵的转置；式(2)和式(5)的非线性最小二乘问题的求解归结为其法方程的求解如式(6)；S5：选择、交叉和变异在S3中，若不满足|fmax‑fmin|＜ε，需进行种群的选择、交叉和变异：其中，选择规则是按照适应度大小对个体进行排序，从中选取i个适应度最大个体，按照下式(7)进行选择概率Ps(xi)：Ps(xi)=f(xi)/Σif(xi)---(7)]]>交叉规则采用两位交叉方式，自适应交叉概率Pc为式(8)：PC=0.9-0.3×(fmax-f′)fmax-favgf′≥favg0.9f′≥favg---(8)]]>式(8)中，favg为群体的平均适应度，f'为两交叉个体较大的适应值，fmax为当前种群中最优个体的函数值；变异规则是将每个个体的每两位二进制编码随机产生一个在(0，1)之间的数，若大于个体的变异率，则该编码由1变成0，或者由0变成1，否则该编码不变异，Pm变异概率为式(9)：Pm=0.2-0.5×(fmax-f)fmax-favgf≥favg0.2f≥favg---(9)]]>式(9)中，favg为群体的平均适应度，f为每个个体的每两位二进制编码随机产生一个在(0，1)之间的数，fmax为当前种群中最优个体的函数值；然后，保留较优个体，用子代种群的适应值较大的个体替代父代种群适应值较小的相同数目个体，即P(t)＝P(t+1)，提高获得最优个体的概率，将当前保留的新种群作为步骤S2的初始种群重新计算；S6：判断式(6)的则转S12；否则转S7；S7：计算δQ,fk,fk+1,f1k+1通过式(6)的信赖域方法确定适合的δQ：Δf＝fk+1‑fk        (10)Δf1＝f1k+1‑fk式中Δf是非线性目标值的改变，Δf1是线性化目标值的改变，同时fk+1,fk分别是k+1和k次迭代非线性目标值，f1k+1是是k+1次迭代线性化目标值，具体计算方法如式(11)：式(11)中，fk+1,fk分别是k+1和k次迭代非线性目标值，f1k+1是是k+1次迭代线性化目标值，W为线性算子，Q和分别为未知模拟电荷的电量向量和已知的匹配点的电位向量；S8：判断...

【技术特征摘要】
1.一种直流输电线路电场强度及无线电干扰计算的方法，其特征在于，S1：确定模拟电荷空间位置模拟电荷的位置满足式(1)的约束条件；为提高分辨率将分裂子导线和避雷线导线采用10位二进制编码，镜像电荷采用10位二进制编码；式(1)中(xi,yi)和(xj,yj)为分裂子导线和避雷线导线内的模拟电荷坐标，(xk,yk)为大地表面以下的镜像电荷坐标，R1和R2为分裂子导线和避雷线导线的半径，A和N分别为避雷线导线内模拟电荷的数量和总的模拟电荷数量；f(xk,yk)为描述大地表面连续曲线的函数；S2：模拟电荷位置初始种群P(t)的产生分裂子导线和避雷线导线的模拟电荷分别均分布于直径1/4的圆环，各镜像电荷分别位于相对地面垂直的等距离镜像位置；将式(2)作为适应度函数评价个体的优劣：式(2)中M为分裂子导线表面的匹配点数量，Q和分别未知模拟电荷的电量向量和已知的匹配点的电位向量，为分裂子导线、避雷线导线或大地表面第i个匹配点的电位；S3：判断收敛性若|fmax-fmin|<ε，则输出寻优结果，停止迭代，转步骤S4；否则转入S5；其中fmax,fmin分别为当前种群中最优和最差个体的函数值，ε为给定准确度；S4：正则化在式(2)中，为实现对非线性算子F(Q)的近似，利用泰勒公式将F(Q+δQ)在Q处的泰勒展开，当||δQ||足够小时，采用一次近似且忽略其余的高阶小量，得到式(3)；令Q*＝Q+δQ是方程的准确解，则在接近于Q*的Q处可以由式(3)得到如下线性算子方程式(4)：结合线性化处理的式(4)，基于全局正则化且通过优化问题得出一个最佳的δQ如式(5)：式(5)中，W为线性算子，η为信赖域大小，μ(η)为罚函数；对于(5)式所示的优化问题，利用f为极小解的必要条件是f的梯度为零的点，得到式(6)的线性方程：式(6)中，为雅可比矩阵的转置；式(2)和式(5)的非线性最小二乘问题的求解归结为其法方程的求解如式(6)；S5：选择、交叉和变异在S3中，若不满足|fmax-fmin|<ε，需进行种群的选择、交叉和变异：其中，选择规则是按照适应度大小对个体进行排序，从中选取i个适应度最大个体，按照下式(7)进行选择概率Ps(xi)：交叉规则采用两位交叉方式，自适应交叉概率Pc为式(8)：式(8)中，favg为群体的平均适应度，f'为两交叉个体较大的适应值，fmax为当前种群中最优个体的函数值；变异规则是将每个个体的每两位二进制编码随机产生一个在(0，1)之间的数，若大于个体的变异率，则该编码由1变成0，或者由0变成1，否则该编码不变异，Pm变异概率为式(9)：

【专利技术属性】
技术研发人员：沈鑫，闫永梅，曹敏，张林山，
申请(专利权)人：云南电网有限责任公司电力科学研究院，
类型：发明
国别省市：云南;53