CN109933883A - 一种基于多元线性回归的阻性电流增长率计算方法 - Google Patents

一种基于多元线性回归的阻性电流增长率计算方法 Download PDF

Info

Publication number
CN109933883A
CN109933883A CN201910175756.2A CN201910175756A CN109933883A CN 109933883 A CN109933883 A CN 109933883A CN 201910175756 A CN201910175756 A CN 201910175756A CN 109933883 A CN109933883 A CN 109933883A
Authority
CN
China
Prior art keywords
linear regression
phase angle
multiple linear
current
growth rate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910175756.2A
Other languages
English (en)
Other versions
CN109933883B (zh
Inventor
韦海荣
吴有文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Jiangsu Electric Power Co ltd Innovation And Innovation Center
Original Assignee
NANJING DAONANENG TECHNOLOGY CO LTD
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NANJING DAONANENG TECHNOLOGY CO LTD filed Critical NANJING DAONANENG TECHNOLOGY CO LTD
Priority to CN201910175756.2A priority Critical patent/CN109933883B/zh
Publication of CN109933883A publication Critical patent/CN109933883A/zh
Application granted granted Critical
Publication of CN109933883B publication Critical patent/CN109933883B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明涉及一种基于多元线性回归的阻性电流增长率计算方法,包括:采样、多元线性回归分析、相角差转换、健康状态判断等步骤。本发明提供的基于多元线性回归的阻性电流增长率计算方法,能够间接去除环境因素对三相泄露电流相角差的影响,定量的表示其与环境因素的关系;能够针对任意布局下的避雷器阻性电流增长率计算;能够预测不同环境下的三相泄漏电流的相角差;能避免因为天气原因而导致误判避雷器的健康状况。

Description

一种基于多元线性回归的阻性电流增长率计算方法
技术领域
本发明涉及一种基于多元线性回归的阻性电流增长率计算方法,属于电力***安全控制技术领域。
背景技术
金属氧化物避雷器(Metal Oxide Surge Arrester,简记为MOA)以其优异的非线性特性和大通流容量等优点,成为电力***过电压保护的主要设备,能有效抑制电力***中的操作过电压和大气过电压。
泄漏电流中的阻性电流能有效的反映避雷器健康状况,泄漏电流幅值对阻性电流并不是很敏感,但泄漏电流相加却能间接且有效的反映出阻性电流的增长率。但是泄漏电流的相角很容易受到外界耦合电流的影响,那么去除外界耦合电流的影响就显得尤为重要,而外界耦合电流受很多环境因素的影响,例如温度、湿度、光照强度和大气压强等。
回归分析是研究变量之间相关关系的数学工具,回归分析可以准确的计量各个因素之间的相关程度与回归拟合程度的高低,借助回归分析的方法,在衡量相角与环境各因素之间的关系时,能用解析式显示的表达泄漏电流与环境之间的关系。
发明内容
本发明要解决技术问题是:克服上述技术的缺点,提供一种明确了温度、湿度、光照强度和大气压强对三相泄漏电流相角差的影响,能够将实际环境下测得的相角差转化为标准状态下的相角的电流增长率计算方法。
为了解决上述技术问题,本发明提出的技术方案是:一种基于多元线性回归的阻性电流增长率计算方法,包括如下步骤:
(1)采样:采集不同温度、湿度、光照强度和大气压强下的三相避雷器的泄漏电流相角差;
(2)利用多元线性回归分析,将上步骤中的温度、湿度、光照强度和大气压强作为自变量,三相避雷器的泄漏电流相角差作为因变量,拟合曲线;
(3)通过测量得到当前的温度、湿度、光照强度和大气压强,以及泄漏电流相角差的实测值,根据拟合曲线计算得到泄漏电流相角差的估计值,将泄漏电流相角差的实测值与估计值之间的差作为对标准状态下泄漏电流相角差的补偿,得到优化后的标准状态下三相泄漏电流相角差,并根据优化后的标准状态下泄漏电流相角差计算阻性电流增长率;
(4)根据已知的阻性电流增长率与避雷器健康状态的关系,得到当前避雷器的健康状况。
上述方案进一步的改进在于:所述步骤(2)中的多元线性回归方法为:设温度为T、湿度为RH、光照强度为L和大气压强为P,Δ∠XY是三相避雷器中的任意两项的泄漏电流相角差;假定它们之间有如下线性关系:Δ∠XY=β01T+β2RH+β3L+β4P;
其中Δ∠XY是可测量的随机变量,βi(i=0,…,4)是回归系数;
取n组不同的自变量(Ti,RHi,Li,Pi)(i=1,2,…,n),每组自变量对应一个测量得到的泄漏电流相角差;根据上述线性关系使用最小二乘法求解多元线性回归方程。
上述方案进一步的改进在于:对线性回归方程进行显著性检验;
对于给定的显著性水平α,查F分布得到临界值Fα(k,n-k-1),并根据F检验计算检验统计量F值,当F>Fα(k,n-k-1)时认为Δ∠XYi组成的矩阵与T、RH、L和P中任意一个因素组成的矩阵之间存在线性相关关系,所建立的线性回归方程式显著的;否则所建立的线性回归方程式不显著的。
上述方案进一步的改进在于:根据求得的多元线性回归方程得到确定系数R,确定系数的表达式:
其中yi为原始数据,为预测数据,为原始数据的均值。
上述方案进一步的改进在于:进行多重共线性诊断;采用基于方差因子膨胀法的多重共线性诊断;根据自变量关于其余自变量的多元线性回归,计算模型的判定系数,定义为第i个自变量的方差膨胀因子:
基于方差膨胀因子的多重共线性诊断规则为:VIFi<5,认为不存在共线性;5≤VIFi≤10,认为存在中等程度共线性;VIFi>10,认为共线性严重;若自变量全无依赖关系时,那么线性回归拟合完成。
上述方案进一步的改进在于:所述步骤(2)中的多元线性回归方法为多项式回归拟合。
上述方案进一步的改进在于:进行多次拟合曲线,并选用拟合优度高的拟合曲线。
上述方案进一步的改进在于:在所述步骤(1)之前,还包括:确定与设备发生缺陷相关联的关联因素。
本发明提供的基于多元线性回归的阻性电流增长率计算方法,能够间接去除环境因素对三相泄露电流相角差的影响,定量的表示其与环境因素的关系;能够针对任意布局下的避雷器阻性电流增长率计算;能够预测不同环境下的三相泄漏电流的相角差;能避免因为天气原因而导致误判避雷器的健康状况。由于存在外界干扰因素,因此传统方案会考虑避雷器布局因素。本发明,通过环境条件和实测值进行相关性拟合,将外界干扰隐含的体现在回归系数中,因此,不需要考虑布局。
附图说明
图1为本发明一个优选的实施例的流程图。
图2是任意布局下的单个避雷器受外部耦合电容影响模型图。
具体实施方式
实施例
本实施例针对任意布局下的三相避雷器排列方式,以A相为例,建立单个避雷器受外部耦合电容影响的模型,如图2所示。图中为电网三相电压矢量;C和R为自身相避雷器的电容与非线性电阻;Cba、Cca为B相与C相对A相避雷器的耦合电容;Ci(i=1,…,n)为其他邻近相避雷器对A相的耦合电容;为耦合电容相对应的电压矢量。根据电路原理可知,所有耦合电容所产生的耦合电流都可以用一个耦合电容和相应的电压矢量来等效。由于电网的标准,电压矢量几乎是不变的,又耦合电容受环境因素的影响严重,所以借助回归分析衡量耦合电容受环境因素的影响。
本实施例的基于多元线性回归的阻性电流增长率计算方法(流程图见图1),包括如下步骤:
(1)采样:采集不同温湿度、光照强度和大气压强下的三相避雷器的泄漏电流相角差;
(2)利用多元线性回归分析,将上步骤中的温湿度、光照强度和大气压强作为自变量,三相避雷器的泄漏电流相角差作为因变量,拟合曲线;
(3)通过测量得到当前的温度、湿度、光照强度和大气压强,以及泄漏电流相角差的实测值,根据拟合曲线计算得到泄漏电流相角差的估计值,将泄漏电流相角差的实测值与估计值之间的差作为对标准状态下泄漏电流相角差的补偿,得到优化后的标准状态下三相泄漏电流相角差,并根据优化后的标准状态下泄漏电流相角差计算阻性电流增长率;
(4)根据国标阻性电流增长率与避雷器健康状态的关系,得到当前避雷器的健康状况。
在本实施例方案中回归分析采用两种方法:线性回归分析和多项式回归分析。在进行回归分析前,先对数据进行适当的分段。
线性回归分析法原理如下,以A相与B相的泄漏电流相角差为例,设影响泄露电流的各自变量参数温度T、湿度RH、光照强度L和大气压强P为影响泄露电流I的四个因素,假定他们之间有如下线性关系:
Δ∠AB=β01T+β2RH+β3L+β4P;
其中I是可测量的随机变量,βi(i=0,…,4)是未知参数,称为回归系数。
取n组不同的(Ti,RHi,Li,Pi)(i=1,2,…,n),每组自变量对应一个测量得到的泄漏电流相角差;根据上述线性关系使用最小二乘法求解多元线性回归方程。
利用最小二乘法求多元线性回归方程后,必须对线性回归方程进行显著性检验。对于给定的显著性水平α,可以查F分布得到临界值Fα(k,n-k-1)。所以当F>Fα(k,n-k-1)时认为Δ∠ABi组成的矩阵与T、RH、L和P中任意一个因素组成的矩阵之间存在线性相关关系,所建立的线性回归方程式显著的;否则所建立的线性回归方程式不显著的。
根据求得的回归方程得到确定系数R,确定系数R越接近1表明自变量对因变量解释能力越好。确定系数的表达式如下
其中yi为原始数据,为预测数据,为原始数据的均值,
在显著性分析中,多重共线性指的是多元回归模型的自变量间存在近似的线性关系,它的存在使得估计的精确性大幅度降低,估计值稳定性变差,甚至会导致一些回归系数通不过显著性检验,正负号倒置,使得无法从回归方程中得到合理的物理意义的解释,降低了回归方程的应用价值。多重共线性诊断的方法有很多,本实施例采用基于方差因子膨胀法的多重共线性诊断。
本实施例采用的基于方差膨胀因子的多重共线性诊断,考虑到自变量关于其余自变量的多元线性回归,计算模型的判定系数,定义为第i个自变量的方差膨胀因子:
当自变量有依赖于其他自变量的线性关系时,接近于1,VIFi接近于无穷大;反之,接近于0,VIFi接近于1。VIFi越大说明线性依赖关系越严重,即存在共线性。通常情况下,基于方差膨胀因子的多重共线性诊断规则为:VIFi<5,认为不存在共线性(或共线性较弱);5≤VIFi≤10,认为存在中等程度共线性;V/Fi>10,认为共线性严重。若自变量全无依赖关系时,那么线性回归拟合完成。
另一种拟合方式为多项式回归拟合,此处以单变量多项式回归模型为例,具体实际中可以采用多变量,视情况而定。
作为单变量多项式回归模型的例子,考虑
y=β01x+β2x2+ε;
该模型称为单变量的二阶模型。有时也称之为二次模型,这是由于y的期望值为
E(y)=β01x+β2x2
通常将β0称为均值系数,β1称为线性效应系数,而将β2称为二次效应系数。
与多元线性回归相似,进行最小二乘法求解,并得到确定系数,最后进行多重共线性分析。
通过以上两种当中拟合效果最好即确定系数最接近1的确定拟合的曲线,到此采样阶段就完成了。接下来就是实际应用阶段,利用实际测得的环境变量与拟合的曲线,将不同环境下的三相电流相角差转化为标准温度下的相角差,并利用单相阻性电流增长率与三相相角差之间的关系得到单相阻性电流增长率,与国家标准对照进行判断避雷器的健康状况。
以A相为例,单相阻性电流增长率与三相泄漏电流相角差的关系如下:A相相角变大,意味着AB之间的相角变大,AC之间的相角变小,那么与健康状况下的相角差之间的变化就是A相的变化量,通过对单相阻性电流增长率导致单相相角的变化量为20%/1°,通过A相相角的变化得到A相阻性电流增长率
本发明不局限于上述实施例。凡采用等同替换形成的技术方案,均落在本发明要求的保护范围。

Claims (8)

1.一种基于多元线性回归的阻性电流增长率计算方法,其特征在于,包括如下步骤:
(1)采样:采集不同温度、湿度、光照强度和大气压强下的三相避雷器的泄漏电流相角差;
(2)利用多元线性回归分析,将上步骤中的温度、湿度、光照强度和大气压强作为自变量,三相避雷器的泄漏电流相角差作为因变量,拟合曲线;
(3)通过测量得到当前的温度、湿度、光照强度和大气压强,以及泄漏电流相角差的实测值,根据拟合曲线计算得到泄漏电流相角差的估计值,将泄漏电流相角差的实测值与估计值之间的差作为对标准状态下泄漏电流相角差的补偿,得到优化后的标准状态下三相泄漏电流相角差,并根据优化后的标准状态下泄漏电流相角差计算阻性电流增长率;
(4)根据已知的阻性电流增长率与避雷器健康状态的关系,得到当前避雷器的健康状况。
2.根据权利要求1所述的基于多元线性回归的阻性电流增长率计算方法,其特征在于,所述步骤(2)中的多元线性回归方法为:设温度为T、湿度为RH、光照强度为L和大气压强为P,Δ∠XY是三相避雷器中的任意两项的泄漏电流相角差;假定它们之间有如下线性关系:Δ∠XY=β01T+β2RH+β3L+β4P;
其中Δ∠XY是可测量的随机变量,βi(i=0,…,4)是回归系数;
取n组不同的自变量(Ti,RHi,Li,Pi)(i=1,2,…,n),每组自变量对应一个测量得到的泄漏电流相角差;根据上述线性关系使用最小二乘法求解多元线性回归方程。
3.根据权利要求2所述的基于多元线性回归的阻性电流增长率计算方法,其特征在于:对线性回归方程进行显著性检验;
对于给定的显著性水平α,查F分布得到临界值Fα(k,n-k-1),并根据F检验计算检验统计量F值,当F>Fα(k,n-k-1)时认为Δ∠XYi组成的矩阵与T、RH、L和P中任意一个因素组成的矩阵之间存在线性相关关系,则所建立的线性回归方程式显著的;否则所建立的线性回归方程式不显著的。
4.根据权利要求3所述的基于多元线性回归的阻性电流增长率计算方法,其特征在于:根据求得的多元线性回归方程得到确定系数R,确定系数的表达式:
其中yi为原始数据,为预测数据,为原始数据的均值。
5.根据权利要求3所述的基于多元线性回归的阻性电流增长率计算方法,其特征在于:进行多重共线性诊断;采用基于方差因子膨胀法的多重共线性诊断;根据自变量关于其余自变量的多元线性回归,计算模型的判定系数,定义为第i个自变量的方差膨胀因子:
基于方差膨胀因子的多重共线性诊断规则为:VIFi<5,认为不存在共线性;5≤VIFi≤10,认为存在中等程度共线性;VIFi>10,认为共线性严重;若自变量全无依赖关系时,那么线性回归拟合完成。
6.根据权利要求1所述的基于多元线性回归的阻性电流增长率计算方法,其特征在于,所述步骤(2)中的多元线性回归方法为多项式回归拟合。
7.根据权利要求1所述的基于多元线性回归的阻性电流增长率计算方法,其特征在于,进行多次拟合曲线,并选用拟合优度高的拟合曲线。
8.根据权利要求1所述的基于条件概率的电力设备缺陷事件预警方法,其特征在于:在所述步骤(1)之前,还包括:确定与设备发生缺陷相关联的关联因素。
CN201910175756.2A 2019-03-08 2019-03-08 一种基于多元线性回归的阻性电流增长率计算方法 Active CN109933883B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910175756.2A CN109933883B (zh) 2019-03-08 2019-03-08 一种基于多元线性回归的阻性电流增长率计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910175756.2A CN109933883B (zh) 2019-03-08 2019-03-08 一种基于多元线性回归的阻性电流增长率计算方法

Publications (2)

Publication Number Publication Date
CN109933883A true CN109933883A (zh) 2019-06-25
CN109933883B CN109933883B (zh) 2023-01-24

Family

ID=66986437

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910175756.2A Active CN109933883B (zh) 2019-03-08 2019-03-08 一种基于多元线性回归的阻性电流增长率计算方法

Country Status (1)

Country Link
CN (1) CN109933883B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112504348A (zh) * 2020-12-11 2021-03-16 厦门汇利伟业科技有限公司 一种融合环境因素的物体状态显示方法和***
CN113064013A (zh) * 2021-03-31 2021-07-02 广东电网有限责任公司 一种避雷器状态监测方法、装置及***
CN116680621A (zh) * 2023-08-02 2023-09-01 国网江苏省电力有限公司 一种避雷器相角差预测方法、装置、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6076048A (en) * 1997-09-26 2000-06-13 Betzdearborn, Inc. System and method for least squares filtering based leak flow estimation/detection using exponentially shaped leak profiles
CN105514990A (zh) * 2015-12-10 2016-04-20 中国电力科学研究院 综合经济性与安全性的输电线路利用率提高平台及方法
CN109406926A (zh) * 2018-09-10 2019-03-01 国网江苏省电力有限公司连云港供电分公司 环境因素对避雷器泄漏电流影响的计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6076048A (en) * 1997-09-26 2000-06-13 Betzdearborn, Inc. System and method for least squares filtering based leak flow estimation/detection using exponentially shaped leak profiles
CN105514990A (zh) * 2015-12-10 2016-04-20 中国电力科学研究院 综合经济性与安全性的输电线路利用率提高平台及方法
CN109406926A (zh) * 2018-09-10 2019-03-01 国网江苏省电力有限公司连云港供电分公司 环境因素对避雷器泄漏电流影响的计算方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112504348A (zh) * 2020-12-11 2021-03-16 厦门汇利伟业科技有限公司 一种融合环境因素的物体状态显示方法和***
CN113064013A (zh) * 2021-03-31 2021-07-02 广东电网有限责任公司 一种避雷器状态监测方法、装置及***
CN116680621A (zh) * 2023-08-02 2023-09-01 国网江苏省电力有限公司 一种避雷器相角差预测方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN109933883B (zh) 2023-01-24

Similar Documents

Publication Publication Date Title
CN109933883A (zh) 一种基于多元线性回归的阻性电流增长率计算方法
CN102830319B (zh) 一种氧化锌避雷器绝缘状态带电检测装置及方法
CN104142421B (zh) 变电设备绝缘在线监测***及其工作方法
CN116879662A (zh) 基于数据分析的变压器故障检测方法
CN110389289A (zh) 绝缘故障检测方法及控制装置
CN107656181B (zh) 一种定量诊断油纸绝缘老化程度的方法
CN106597350B (zh) 一种关口电能表运行误差评估方法及***
CN102288804B (zh) 一种基于正交变换的避雷器阻性电流计算方法
CN113238142A (zh) 用于集成电路的方法和***
CN114034997A (zh) 一种基于多参量的绝缘子劣化程度预测方法及***
Papliński et al. Application of leakage current parameters for technical diagnostics of surge arresters
CN105891625A (zh) 基于能量流的电力***次同步振荡扰动源辨别方法
CN108828292B (zh) 考虑接地排影响的氧化锌避雷器阻性电流在线监测方法
CN112083299B (zh) 一种基于卡尔曼滤波的直流***绝缘故障预测方法
CN114660387A (zh) 基于泄漏电流传感器和bp神经网络算法的避雷器监控方法
CN110988444A (zh) 氧化锌避雷器分压装置、过电压检测避雷装置及检测方法
CN112989587A (zh) 一种电容式电压互感器劣化成因在线分析方法及***
CN109406926A (zh) 环境因素对避雷器泄漏电流影响的计算方法
CN1797014A (zh) 避雷器阻性电流测试的检测方法
CN111965511A (zh) 一种晶闸管电路的故障检测方法、***及介质
CN115575880B (zh) 一种电力互感器计量误差状态在线检测方法
CN109298033B (zh) 一种电化学气体传感器及其诊断方法
CN114325176B (zh) 一种氧化锌避雷器电阻阀片受潮老化的性能评估方法
CN106053984B (zh) 确定防滑刹车控制装置高温故障分布的方法
CN109324223A (zh) 一种高原型氧化锌避雷器在线监测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20231108

Address after: Room 1019, 10th Floor, No. 20 Beijing West Road, Gulou District, Suzhou City, Jiangsu Province, 210024

Patentee after: State Grid Jiangsu Electric Power Co.,Ltd. innovation and Innovation Center

Address before: Room 1406, Building 1, No. 251, Heyan Road, Maigaoqiao Street, Qixia District, Nanjing City, Jiangsu Province, 210000

Patentee before: NANJING DAONANENG TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right
CP02 Change in the address of a patent holder

Address after: 210024 Room 1019, 10th Floor, No. 20, Beijing West Road, Gulou District, Nanjing, Jiangsu Province

Patentee after: State Grid Jiangsu Electric Power Co.,Ltd. innovation and Innovation Center

Address before: Room 1019, 10th Floor, No. 20 Beijing West Road, Gulou District, Suzhou City, Jiangsu Province, 210024

Patentee before: State Grid Jiangsu Electric Power Co.,Ltd. innovation and Innovation Center

CP02 Change in the address of a patent holder