一种基于灵敏度分析的区域地区直流波动源定位方法技术

本发明专利技术提供一种基于灵敏度分析的区域地区直流波动源定位方法，该方法通过获取区域中变电站变压器中性点直流电流INi和线波动源的经度X、维度Y及电流I；利用获取的参数来建立的线波动源模型；再计算区域电网中所有中性点接地的变电站对所有线直流波动源的灵敏度系数，并形成区域交流电网变电站直流信号灵敏度列表；然后根据灵敏度列表确定选择监测站点的选择方法，并对其进行蒙特卡洛仿真；最后通过建立的线波动源模型对区域直流波动源进行定位。

【技术实现步骤摘要】
一种基于灵敏度分析的区域地区直流波动源定位方法
本专利技术涉及输变电
，更具体地，涉及一种基于灵敏度分析的区域地区直流波动源定位方法。
技术介绍
大地是不良导体，一般情况下是没有大幅电流流过的，但最近多项研究发现，如地磁暴、高空核弹爆炸、城市轨道交通的杂散电流、直流输电系统大地返回运行等情况下均会引起高幅值的电流从地中流过，引发大范围地表电位的明显波动，在电力系统接地运行的变压器中性点上产生直流电流或者是近似于直流电流的甚低频电流这个电流不仅危及交流电网的安全运行，也会对地下金属管道和工业用电设备产生众多不利影响，还会干扰地震监测台网的正常工作。这些影响不仅会造成大地电位波动源附近地区的电力、埋地金属设备损坏，地震监测台网误报，还会产生一系列连锁反应，造成重大经济损失，给现代社会带来巨大的负面影响。因此，有必要分析研究不同波动干扰源的类型，确定其位置和强度，排除交流电网中的隐患，从而对电力系统和地下管网进行必要修缮，确保电力系统和管网系统的安全稳定运行。目前在我国的多个地区，多个接地运行的电力变压器在排除了直流接地极以及轨道交通的影响后仍然检测出了噪声异常的情况。经变压器中性点电流监测的直流分量测试分析发现，主变噪声是随直流分量幅值的变化高低起伏，即主变出现直流偏磁导致噪声异常，因此可推测出其它直流干扰源的影响。输电网络由于其广阔的地理分布和低阻网络特性，地电位波动能够在电力系统内产生特征电流分布，所以输电网络是一个能够很好地发现地电位波动的监测网络。对变压器中性点电流进行有效监测，利用检测值反演地电位波动源的强度和位置十分重要。
技术实现思路
本专利技术提供一种基于灵敏度分析的区域地区直流波动源定位方法，该方法为干扰源的排查与抑制地电位波动造成的不利影响实施提供依据。为了达到上述技术效果，本专利技术的技术方案如下：一种基于灵敏度分析的区域地区直流波动源定位方法，包括以下步骤：S1：获取区域中变电站变压器中性点直流电流INi和线波动源的经度X、维度Y及电流I；S2：利用S1中的参数来建立的线波动源模型；S3：计算区域电网中所有中性点接地的变电站对所有线直流波动源的灵敏度系数，并形成区域交流电网变电站直流信号灵敏度列表；S4：根据灵敏度列表确定选择监测站点的选择方法，并对其进行蒙特卡洛仿真；S5：通过建立的线波动源模型对区域直流波动源进行定位。进一步地，所述步骤S2的过程是：令变压器T1的中性点直流电流为IN1，当考虑背景噪声时，T1的中性点直流电流与线波动源电流的关系满足：IN1＝M1(X,Y)I+JN1为求取经度X，纬度Y和电流I这三个未知数，结合三个不同位置变压器的中性点电流测量值，形成如式方程：互电阻函数M1、M2和M3是非线性函数，JN1、JN2、JN3是噪声信号，由于交流系统的参数是确定的，变化的只是X，Y和I，波动源定位视为最小二乘优化问题，即：使用迭代形式的优化方法就可以求解上式的目标函数最优解，通过更多的监测点来实现波动源定位，上式变成：式中，k为监测变压器的数目。进一步地，所述步骤S3的过程是：式中，η为灵敏度系数，m和n分别为变电站和波动源的编号，IN为变压器中性点电流，I为波动源电流。进一步地，所述步骤S4的具体过程是：1)、设置I为非零值，令抽样样本数目为N；2)、在交流电网区域内按均匀分布的形式产生N个随机的X和Y；3)、在完成蒙特卡洛随机仿真的N个样本后，对所有站点的M值进行加权平均：式中，ki为加权系数，如(Xi，Yi)落入重点监控区域，加权系数应适当增加，对于一般情况有ki＝1。与现有技术相比，本专利技术技术方案的有益效果是：本专利技术通过获取区域中变电站变压器中性点直流电流INi和线波动源的经度X、维度Y及电流I；利用获取的参数来建立的线波动源模型；再计算区域电网中所有中性点接地的变电站对所有线直流波动源的灵敏度系数，并形成区域交流电网变电站直流信号灵敏度列表；然后根据灵敏度列表确定选择监测站点的选择方法，并对其进行蒙特卡洛仿真；最后通过建立的线波动源模型对区域直流波动源进行定位。本专利技术方法从逆向思维来考虑波动源与电力系统内特征电流分布响应的问题，进一步完善了直流分布计算。通过公式的推导，在理论上证明了测量数据反演地电位波动源的可行性；充分利用变压器中性点监测值，构建电位波动测量网络，为干扰源的排查与抑制地电位波动造成的不利影响实施提供依据；模型简洁，结合蒙特卡洛算法，可达到较高的计算精度。附图说明图1为本专利技术方法流程图。具体实施方式附图仅用于示例性说明，不能理解为对本专利的限制；为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。下面结合附图和实施例对本专利技术的技术方案做进一步的说明。实施例1如图1所示，一种基于灵敏度分析的区域地区直流波动源定位方法，包括以下步骤：S1：获取区域中变电站变压器中性点直流电流INi和线波动源的经度X、维度Y及电流I；S2：利用S1中的参数来建立的线波动源模型；S3：计算区域电网中所有中性点接地的变电站对所有线直流波动源的灵敏度系数，并形成区域交流电网变电站直流信号灵敏度列表；S4：根据灵敏度列表确定选择监测站点的选择方法，并对其进行蒙特卡洛仿真；S5：通过建立的线波动源模型对区域直流波动源进行定位。A、建立线波动源模型：地电位波动源的几何形态可分为单点、直线、曲线、面波动源、及其他形式的点线形式组合。线形态的地电位波动源主要由城市轨道交通引起，我国目前多个城市内变电站内变压器受轨道交通杂散电流的影响。GIC和HEMP会形成面形态的地电位波动源，需确定经度和纬度方向上的地下电场强度。线波动源待求解的参数包括：(1)电流；(2)位置，即经度和纬度。B、波动源参数推导变压器中性点特征电流响应：由直流电流分布原理，对于所有变压器中性点电流列向量IN，有IN＝MI(1)M＝SgΓ1Γ2Φdc(2)Γ1＝[S-Γ2ΦgGgAg]-1(3)上式中，M为比例系数矩阵，M的元素为非线性函数，与大地电阻率参数，变压器位置，电网拓扑/电网网络参数，变电站接地电阻等因素相关；S和Sg分别为交流电网节点电导矩阵和变压器母线–中性点的电导矩阵，由系统参数决定；Φg和Φdc分别是变电站–变电站和干扰源–变电站之间的互阻矩阵，与干扰源位置和变电站位置相关，其表达式是位置的非线性函数，与大地电阻率有关；Gg是变电站的接地电导矩阵；Ag是变电站节点的关联矩阵，与系统的拓扑结构有关，；E是单位矩阵；Γ1和Γ2均为中间矩阵。C、根据监测系统获取直流电流参数。利用线定位原理对区域直流波动源进行定位。假设：(1)线的波动源由经度X，纬度Y和电流I三个未知因素共同确定确定。(2)电力系统中某些变电站的变压器中性点直流电流INi已知。假设变压器T1的中性点直流电流为IN1，当考虑背景噪声时，T1的中性点直流电流与线波动源电流的关系满足IN1＝M1(X,Y)I+JN1(5)为求取经度X，纬度Y和电流I这三个未知数，结合三个不同位置变压器的中性点电流测量值，形成如式方程：虽然互电阻函数M1、M2和M3是非线性函数，但由于交流系统的参数是确定的，变化的只是X，Y和I，所以波动源定位可视为最小二乘本文档来自技高网...

【技术保护点】
一种基于灵敏度分析的区域地区直流波动源定位方法，其特征在于，包括以下步骤：S1：获取区域中变电站变压器中性点直流电流INi和线波动源的经度X、维度Y及电流I；S2：利用S1中的参数来建立的线波动源模型；S3：计算区域电网中所有中性点接地的变电站对所有线直流波动源的灵敏度系数，并形成区域交流电网变电站直流信号灵敏度列表；S4：根据灵敏度列表确定选择监测站点的选择方法，并对其进行蒙特卡洛仿真；S5：通过建立的线波动源模型对区域直流波动源进行定位。

【技术特征摘要】
1.一种基于灵敏度分析的区域地区直流波动源定位方法，其特征在于，包括以下步骤：S1：获取区域中变电站变压器中性点直流电流INi和线波动源的经度X、维度Y及电流I；S2：利用S1中的参数来建立的线波动源模型；S3：计算区域电网中所有中性点接地的变电站对所有线直流波动源的灵敏度系数，并形成区域交流电网变电站直流信号灵敏度列表；S4：根据灵敏度列表确定选择监测站点的选择方法，并对其进行蒙特卡洛仿真；S5：通过建立的线波动源模型对区域直流波动源进行定位。2.根据权利要求1所述的基于灵敏度分析的区域地区直流波动源定位方法，其特征在于，所述步骤S2的过程是：令变压器T1的中性点直流电流为IN1，当考虑背景噪声时，T1的中性点直流电流与线波动源电流的关系满足：IN1＝M1(X,Y)I+JN1为求取经度X，纬度Y和电流I这三个未知数，结合三个不同位置变压器的中性点电流测量值，形成如式方程：互电阻函数M1、M2和M3是非线性函数，JN1、JN2、JN3是噪声信号，由于交流系统的参数是确定的，变化的只是X，Y和I，波动源定位视为最小二乘优化问题，即：

【专利技术属性】
技术研发人员：李雷，王俊波，罗容波，章涛，邢立坤，
申请(专利权)人：广东电网有限责任公司佛山供电局，
类型：发明
国别省市：广东,44