一种基于切比雪夫多项式的地表电位快速评估方法技术

本发明专利技术属于特高压直流输电工程接地极造成大范围地表电位干扰的分析领域，具体涉及一种基于切比雪夫多项式的地表电位快速评估方法，首先采集数据，将采集数据代入切比雪夫多项式快速评估地表电位的数值。与现有技术相比，本发明专利技术具有如下的有益效果：大幅降低了直流接地极造成广大区域地表电位分布的计算难度，切比雪夫多项式快速评估地表电位的数值计算方法有助于形成计算软件。算方法有助于形成计算软件。算方法有助于形成计算软件。

【技术实现步骤摘要】
一种基于切比雪夫多项式的地表电位快速评估方法


[0001]本专利技术属于特高压直流输电工程接地极造成大范围地表电位干扰的分析领域，具体涉及一种直流输电发生单极大地返回运行方式时，直流接地极造成广大区域地表电位分布的计算方法，利用切比雪夫多项式快速评估地表电位的数值计算方法。

技术介绍

[0002]运行经验表明，直流接地极造成大范围的电网直流偏磁危害其根源在于直流接地极自身造成的不均匀地表电位分布。电网直流偏磁危害具体表现为：电网由于通过输电线路连接了相互距离较远的众多变电站，而变电站内高压变压器中性点是直接接地运行的；大地与电力网之间形成了一个串联的电流通道，直流输电的部分入地电流会从大地中被“抽取”到电力系统内；由于有直流电流作用于变压器绕组上，造成变压器的铁芯饱和，从而产生变压器的半周饱和现象，也就是直流偏磁危害。直流偏磁危害具体表现为变压器的振动和异响，以及局部温升和谐波等问题。直流偏磁危害重危害电力系统安全运行。
[0003]直流接地极造成大范围的地表电位V可写为如下形式：
[0004][0005]其中，r为地表与直流接地极的距离，ρ1为首层大地的电阻率，I为直流输电系统的大地返回电流，J0为第一类零阶贝塞尔函数，λ为积分系数，k1为表征水平分层大地结构的等效电阻率函数，可有如下的递推公式来求解：
[0006][0007]其中，水平分层大地结构示意图如图1所示，d为分层的厚度，ρ为分层的电阻率，μ
i
为第i个分层的反射系数。
[0008]式(1)属于广义无穷积分，J0随rλ的增大呈现震荡衰减的特性，传统的数值计算方法，无法准确求解式(1)。如果k1可以进行无穷指数级数展开，式(1)则可以运用经典镜像法求解，但无法针对分层较多的情况作求解。造成这个方面困难的原因是分层较多的情况下k1的无穷指数级数理论表达过于复杂，另外是计算时间太长，无法用于实际。
[0009]目前有学者研究使用复镜像法求解直流接地极造成广大区域地表电位分布，但复镜像法算法复杂，不便于推广应用。有鉴于此，本专利技术提出一种直流输电发生单极大地返回运行方式时直流接地极造成广大区域地表电位分布的计算方法，通过引入切比雪夫多项式
快速评估地表电位，为直流偏磁风险评估提供新思路和新手段。

技术实现思路

[0010]本专利技术的上述技术问题主要是通过下述技术方案得以解决的：
[0011]一种基于切比雪夫多项式的地表电位快速评估方法，其特征在于，包括
[0012]步骤1，采集数据，包括地表与直流接地极的距离r、首层大地的电阻率ρ1、直流输电系统的大地返回电流I、积分系数λ、表征水平分层大地结构的等效电阻率函数k1、分层的厚度d、分层的电阻率ρ、第i个分层的反射系数μ
i
；
[0013]步骤2，设首层的衰减系数E为：
[0014]E＝exp(
‑
2λd1)
ꢀꢀ
(1)
[0015]对于第i层大地，其衰减系数为
[0016][0017]步骤3，公式(1)、公式(2)代入公式(3)，
[0018][0019]即得：
[0020][0021]使用级数展开，即得
[0022][0023]步骤4，公式(5)代入公式(6)
[0024][0025]有
[0026][0027]李谱希茨分公式
[0028][0029]代入公式(7)有
[0030][0031]输出计算结果,输出的结果是地表与直流接地极距离为r处的地表电位；因为输出的结果即为电表电位，最终可以得到距离接地极不同位置的地表电位分布，从而对其进行评估；
[0032]其中，x、y分别为满足李普希茨条件的自变量和因变量，只要满足条件，该等式就成立；
[0033]李普希茨条件：
[0034]设y＝f(x)为定义在[a,b]上的函数，k为一正常数，若对于[a,b]中任意两点x1，x2，都有
[0035]|f(x1)
‑
f(x2)|≤k|x1‑
x2|
[0036]则称f(x)在区间[a,b]上满足李普希茨条件；r为地表与直流接地极的距离，ρ1为首层大地的电阻率，I为直流输电系统的大地返回电流，J0为第一类零阶贝塞尔函数，λ为积分系数，k1为表征水平分层大地结构的等效电阻率函数，d为分层的厚度，ρ为分层的电阻率，μ
i
为第i个分层的反射系数,A
i
是将公式(4)使用级数展开之后的系数，如公式(5)所示；V(r)是指地表与直流接地极的距离为r处的地表电位。
[0037]在上述的一种基于切比雪夫多项式的地表电位快速评估方法，步骤3中需要引入移位运算、系数展开和截断误差判定三大环节来求解式(5)；
[0038]步骤3.1，移位运算；切比雪夫多项式的值域范围为[
‑
1，1]，而系数k1的取值范围为[0，1]；因此切比雪夫多项式要作移位运算，移位切比雪夫多项式T
j
表达式取为
[0039][0040]步骤3.2，系数展开；进行k1的移位切比雪夫多项式展开，具体表达式取为
[0041][0042]其中，c
j
为T
j
的系数，c
j
的具体表达式为
[0043][0044]其中，ξ
b
为T
b
系数方程的根，ξ
b
的具体表达式为
[0045][0046]步骤3.3，截断误差判定；m阶移位切比雪夫多项式表达式的误差表达式可写为
[0047][0048]式(11)近似取误差主项
[0049]ε
m
(E)≈c
m+1
T
m+1
(E)
ꢀꢀ
(15)
[0050]E值域误差的总误差为
[0051][0052]式(13)积分的理论结果为
[0053][0054]故可按如下接待的流程确定公式(8)中的最高阶数m。
[0055]在上述的一种基于切比雪夫多项式的地表电位快速评估方法，步骤3中，确定公式(8)中的最高阶数m包括如下步骤：
[0056]步骤1、设定允许误差ε＝1e
‑
6，令m＝1；
[0057]步骤2、按照公式(8)计算c
j
和ε
tot
[0058]步骤3、判断ε
tot
与ε的大小关系：
[0059]当ε
tot
＞ε时，令m＝m+1，重复步骤2、3；
[0060]当ε
tot
＜ε时，直接输出c
j
和ε
tot
。
[0061]与现有技术相比，本专利技术具有如下的有益效果：大幅降低了直流接地极造成广大区域地表电位分布的计算难度，切比雪夫多项式快速评估地表电位的数值计算方法有助本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于切比雪夫多项式的地表电位快速评估方法，其特征在于，包括步骤1，采集数据，包括地表与直流接地极的距离r、首层大地的电阻率ρ1、直流输电系统的大地返回电流I、积分系数λ、表征水平分层大地结构的等效电阻率函数k1、分层的厚度d、分层的电阻率ρ、第i个分层的反射系数μ
i
；步骤2，设首层的衰减系数E为：E＝exp(
‑
2λd1)
ꢀꢀꢀꢀꢀꢀꢀ
(1)对于第i层大地，其衰减系数为步骤3，公式(1)、公式(2)代入公式(3)，即得：使用级数展开，即得步骤4，公式(5)代入公式(6)有李谱希茨分公式
代入公式(7)有输出计算结果,输出的结果是地表与直流接地极距离为r处的地表电位；因为输出的结果即为电表电位，最终可以得到距离接地极不同位置的地表电位分布，从而对其进行评估；其中，x、y分别为满足李普希茨条件的自变量和因变量，只要满足条件，该等式就成立；李普希茨条件：设y＝f(x)为定义在[a,b]上的函数，k为一正常数，若对于[a,b]中任意两点x1，x2，都有|f(x1)
‑
f(x2)|≤k|x1‑
x2|则称f(x)在区间[a,b]上满足李普希茨条件；r为地表与直流接地极的距离，ρ1为首层大地的电阻率，I为直流输电系统的大地返回电流，J0为第一类零阶贝塞尔函数，λ为积分系数，k1为表征水平分层大地结构的等效电阻率函数，d为分层的厚度，ρ为分层的电阻率，μ
i
为第i个分层的反射系数,A
i
是将公式(4)使用级数展开之后的系数，如公式(5)所示；V(r)是指地表与直流接地极的距离为r处的地表电位。2.根据权利要求1所述的一种基于切比雪夫多项式的地表电位快速评估方法，其特征在...

【专利技术属性】
技术研发人员：齐桓若，陈晨，闫向阳，康祎龙，王顺然，刁旭，宋景博，李凯，张金凤，薛文杰，翟育新，宋晓帆，郭放，
申请(专利权)人：国网河南省电力公司经济技术研究院，
类型：发明
国别省市：