CN112861469A - 配电变压器模型参数估计方法、装置、设备和存储介质 - Google Patents
配电变压器模型参数估计方法、装置、设备和存储介质 Download PDFInfo
- Publication number
- CN112861469A CN112861469A CN202110255432.7A CN202110255432A CN112861469A CN 112861469 A CN112861469 A CN 112861469A CN 202110255432 A CN202110255432 A CN 202110255432A CN 112861469 A CN112861469 A CN 112861469A
- Authority
- CN
- China
- Prior art keywords
- distribution transformer
- transfer function
- function value
- low
- value
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012546 transfer Methods 0.000 claims abstract description 102
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 30
- 230000004044 response Effects 0.000 claims abstract description 29
- 230000002068 genetic effect Effects 0.000 claims abstract description 28
- 230000006870 function Effects 0.000 claims description 123
- 238000004804 winding Methods 0.000 claims description 114
- 239000003990 capacitor Substances 0.000 claims description 29
- 238000010586 diagram Methods 0.000 claims description 25
- 238000009413 insulation Methods 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 22
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 9
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000005094 computer simulation Methods 0.000 abstract description 7
- 238000013459 approach Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Genetics & Genomics (AREA)
- Physiology (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
本申请公开了配电变压器模型参数估计方法、装置、设备和存储介质，该配电变压器模型参数估计方法包括获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值；使用遗传算法估计所述配电变压器的模型参数，利用该方法，无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。本发明提供的配电变压器模型参数估计装置、设备和存储介质，同样能够无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。
Description
技术领域
本发明属于配电变压器技术领域，特别是涉及配电变压器模型参数估计方法、装置、设备和存储介质。
背景技术
配电变压器是配电网的重要资产，当其出现包括停电和着火在内的灾难性后果时，会给电网和用户带来严重的经济利益损失。配电变压器包括铁芯、绕组、绝缘装置在内的组成部分，可以等效为一个电阻、电容和电感连接组成的电路，任何机械变形或绝缘装置的退化都会导致这些参数的变化。因此，建立一个可靠的配电变压器中低频模型对配电变压器的运行状态的分析是必不可少的，工频变压器的模型参数可以通过一些实验测量来确定，如开路和短路测试，中低频频变压器的模型参数一般根据变压器的精确物理参数进行计算得到，但由于厂家生产保密和运行年限过久，部分参数可能无法获得，因此亟需一种应用方便，无需变压器本体参数即可估算等效电路参数的估计方法。
发明内容
为解决上述问题，本发明提供了配电变压器模型参数估计方法、装置、设备和存储介质，无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。
本发明提供的配电变压器模型参数估计方法包括：
获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；
绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值；
使用遗传算法估计所述配电变压器的模型参数。
优选的，在上述配电变压器模型参数估计方法中，所述获取配电变压器的中低频频率响应信号，获得传递函数值的测量值包括：
向所述配电变压器的低压侧的其中一相绕组端子注入一个低振幅、可变频率的正弦波电压，并在高压侧的同相绕组端子处捕获响应信号，得到以中低频频率响应信号表示的所述传递函数值的测量值为：
其中，Vin为所述低压侧的输入电压幅值，Vout为所述高压侧的输出电压幅值。
优选的，在上述配电变压器模型参数估计方法中，所述绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值包括：
根据配电变压器的电磁传递特性，绘制出配电变压器的等效电路图，对图中的各个节点运用基尔霍夫定律，得到配电变压器的等效电路的状态空间微分方程为：
式中，
表示对变量vk取对时间的微分；
其中，GH0为高压绕组与接地箱体间的绝缘电导，CH0为高压绕组与接地箱体间的绝缘电容，GH为高压绕组电导，CH为高压绕组电容，RH为高压绕组电阻，LH为高压绕组电感，GHL为高低压绕组间的绝缘电导，CHL为高低压绕组间的绝缘电容，GL0为低压绕组与接地铁芯间的绝缘电导，CL0为低压绕组与接地铁芯间的绝缘电容，GL为低压绕组电导，CL为低压绕组电容，RL为低压绕组电阻，LL为低压绕组电感，Vj-1、Vj、Vj+1、Vj+2分别为节点j-1、j、j+1、j+2的电压，Ij-1、Ij、Ij+1分别为节点j-1、j、j+1流出的电流，Vi-1、Vi、Vi+1、Vi+2分别为节点i-1、i、i+1、i+2的电压，Ii-1、Ii、Ii+1分别为节点i-1、i、i+1流出的电流；
结合所述等效电路图和所述状态空间微分方程，在某一输入电压和电流频率下各元件参数确定后，利用MATLAB软件计算配电变压器的传递函数值的计算值。
优选的，在上述配电变压器模型参数估计方法中，所述使用遗传算法估计所述配电变压器的模型参数包括：
所述遗传算法的目标函数为测量得到的传递函数值的测量值与传递函数值的计算值之间的误差最小，即
min J＝J1+J2
其中，J为遗传算法目标函数，J1和J2分别为全频率点子目标函数和谐振频率点子目标函数，(|TF|m)measured和(|TF|m)estimated分别为在频率点m处测量得到的传递函数值和估计的传递函数值，n为频率点的总数，(|TF|mr)measured和(|TF|mr)estimated分别为在谐振频率点mr处测量的传递函数值的测量值和传递函数值的计算值，z为谐振频率点总数。
本发明提供的配电变压器模型参数估计装置包括：
传递函数值测量值获得单元，用于获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；
传递函数值计算值计算单元，用于绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值；
模型参数估计单元，用于使用遗传算法估计所述配电变压器的模型参数。
优选的，在上述配电变压器模型参数估计装置中，所述传递函数值测量值获得单元具体用于：
向所述配电变压器的低压侧的其中一相绕组端子注入一个低振幅、可变频率的正弦波电压，并在高压侧的同相绕组端子处捕获响应信号，得到以中低频频率响应信号表示的所述传递函数值的测量值为：
其中，Vin为所述低压侧的输入电压幅值，Vout为所述高压侧的输出电压幅值。
优选的，在上述配电变压器模型参数估计装置中，所述传递函数值计算值计算单元具体用于：
根据配电变压器的电磁传递特性，绘制出配电变压器的等效电路图，对图中的各个节点运用基尔霍夫定律，得到配电变压器的等效电路的状态空间微分方程为：
式中，
表示对变量vk取对时间的微分；
其中，GH0为高压绕组与接地箱体间的绝缘电导，CH0为高压绕组与接地箱体间的绝缘电容，GH为高压绕组电导，CH为高压绕组电容，RH为高压绕组电阻，LH为高压绕组电感，GHL为高低压绕组间的绝缘电导，CHL为高低压绕组间的绝缘电容，GL0为低压绕组与接地铁芯间的绝缘电导，CL0为低压绕组与接地铁芯间的绝缘电容，GL为低压绕组电导，CL为低压绕组电容，RL为低压绕组电阻，LL为低压绕组电感，Vj-1、Vj、Vj+1、Vj+2分别为节点j-1、j、j+1、j+2的电压，Ij-1、Ij、Ij+1分别为节点j-1、j、j+1流出的电流，Vi-1、Vi、Vi+1、Vi+2分别为节点i-1、i、i+1、i+2的电压，Ii-1、Ii、Ii+1分别为节点i-1、i、i+1流出的电流；
结合所述等效电路图和所述状态空间微分方程，在某一输入电压和电流频率下各元件参数确定后，利用MATLAB软件计算配电变压器的传递函数值的计算值。
优选的，在上述配电变压器模型参数估计装置中，所述模型参数估计单元具体用于：
所述遗传算法的目标函数为测量得到的传递函数值的测量值与传递函数值的计算值之间的误差最小，即
min J＝J1+J2
其中，J为遗传算法目标函数，J1和J2分别为全频率点子目标函数和谐振频率点子目标函数，(|TF|m)measured和(|TF|m)estimated分别为在频率点m处测量得到的传递函数值和估计的传递函数值，n为频率点的总数，(|TF|mr)measured和(|TF|mr)estimated分别为在谐振频率点mr处测量的传递函数值的测量值和传递函数值的计算值，z为谐振频率点总数。
本发明提供的一种计算机设备包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如上面任一项所述配电变压器模型参数估计方法的步骤。
本发明提供的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上面任一项所述配电变压器模型参数估计方法的步骤。
通过上述描述可知，本发明提供的上述配电变压器模型参数估计方法，由于包括获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值；使用遗传算法估计所述配电变压器的模型参数，因此能够无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。本发明提供的上述配电变压器模型参数估计装置、设备和存储介质具有与上述方法相同的优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本发明提供的配电变压器模型参数估计方法的实施例的示意图；
图2为本申请采用的配电变压器的等效电路图；
图3为本发明提供的配电变压器模型参数估计装置的实施例的示意图；
图4为本发明提供的一种计算机设备的实施例的示意图。
具体实施方式
本发明的核心是提供配电变压器模型参数估计方法、装置、设备和存储介质，无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供的配电变压器模型参数估计方法的实施例如图1所示，图1为本发明提供的配电变压器模型参数估计方法的实施例的示意图，该方法可以包括如下步骤：
S1：获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；
具体的，可以包括向配电变压器的低压侧的其中一相绕组端子注入一个低振幅、可变频率的正弦波电压，具体的频率可以但不限于为1Hz至5Hz，并在高压侧的同相绕组端子处捕获响应信号，得到以中低频频率响应信号表示的传递函数值的测量值为：
其中，Vin为低压侧的输入电压幅值，Vout为高压侧的输出电压幅值。
S2：绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算配电变压器的传递函数值的计算值；
具体的，可以根据配电变压器的电磁传递特性，绘制出配电变压器的等效电路图，可以参考图2，图2为本申请采用的配电变压器的等效电路图，根据配电变压器的电磁传递特性，在进行配电变压器的等效时，高压绕组模型用高压绕组的电阻RH和电感LH串联后与电导GH、电容CH并联表示；低压绕组模型用低压绕组的电阻RL和电感LL串联后与电导GL、电容CL并联表示；高低压绕组间的绝缘模型由绝缘电容CHL和绝缘电导GHL并联而成；低压绕组与接地铁芯间的绝缘模型由低压绕组与接地铁芯间的绝缘电导GL0和绝缘电容CL0的二分之一并联而成；高压绕组与接地箱体间的绝缘模型由高压绕组与接地箱体间的绝缘电导GH0和绝缘电容CH0的二分之一并联而成。在连接方式上：A相高压绕组模型两端分别接入节点j-1和j，A相低压绕组模型两端分别接入节点i-1和i，节点j-1和i-1间通过高低压绕组间的绝缘模型相连，A相低压绕组与接地铁芯间的绝缘模型一端接入节点j-1、另一端接地，A相高压绕组与接地箱体间的绝缘模型一端接入节点i-1、另一端接地；B、C相的各模型连接方式与A相相同，只是节点编号分别增加1和2。
对图中的各个节点运用基尔霍夫定律，得到配电变压器的等效电路的状态空间微分方程为：
式中，
表示对变量vk取对时间的微分；
其中，GH0为高压绕组与接地箱体间的绝缘电导，CH0为高压绕组与接地箱体间的绝缘电容，GH为高压绕组电导，CH为高压绕组电容，RH为高压绕组电阻，LH为高压绕组电感，GHL为高低压绕组间的绝缘电导，CHL为高低压绕组间的绝缘电容，GL0为低压绕组与接地铁芯间的绝缘电导，CL0为低压绕组与接地铁芯间的绝缘电容，GL为低压绕组电导，CL为低压绕组电容，RL为低压绕组电阻，LL为低压绕组电感，Vj-1、Vj、Vj+1、Vj+2分别为节点j-1、j、j+1、j+2的电压，Ij-1、Ij、Ij+1分别为节点j-1、j、j+1流出的电流，Vi-1、Vi、Vi+1、Vi+2分别为节点i-1、i、i+1、i+2的电压，Ii-1、Ii、Ii+1分别为节点i-1、i、i+1流出的电流；
结合等效电路图和状态空间微分方程，在某一输入电压和电流频率下各元件参数确定后，利用MATLAB软件计算配电变压器的传递函数值的计算值。
S3：使用遗传算法估计配电变压器的模型参数。
需要说明的是，该遗传算法的核心在于：根据中低频输入下测量得到的配电变压器传递函数值的计算值，采用遗传算法估算各元件参数，代入至推导得到传递函数值计算程序中，使计算得到的传递函数值逐渐逼近匹配上述传递函数值的测量值，最终得到满足使用精度要求的变压器模型参数值。
通过上述描述可知，本发明提供的上述配电变压器模型参数估计方法的实施例中，由于包括获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算配电变压器的传递函数值的计算值；使用遗传算法估计配电变压器的模型参数，因此能够无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。
在上述配电变压器模型参数估计方法的一个具体实施例中，使用遗传算法估计配电变压器的模型参数可以包括如下方式：
遗传算法的目标函数为测量得到的传递函数值的测量值与传递函数值的计算值之间的误差最小，即
min J＝J1+J2
其中，J为遗传算法目标函数，J1和J2分别为全频率点子目标函数和谐振频率点子目标函数，(|TF|m)measured和(|TF|m)estimated分别为在频率点m处测量得到的传递函数值和估计的传递函数值，n为频率点的总数，(|TF|mr)measured和(|TF|mr)estimated分别为在谐振频率点mr处测量的传递函数值的测量值和传递函数值的计算值，z为谐振频率点总数。
具体而言，该算法的计算流程可以包括如下4个步骤：
1)初始化：创建配电变压器等效电路元件电气参数的随机初始种群；
2)评估：根据上一步定义的参数和图2的等效电路，计算估计的传递函数响应值，并与实际测量得到的响应值进行比较，评估上述公式中的目标函数；
3)更新：应用选择和繁殖算子(精英论、交叉和变异)创建一个新的种群，并将最优异的个体复制到新创建的种群中；
4)终止：重复上述更新过程，直到目标函数达到全局最小值或迭代次数达到最大值。
本发明提供的配电变压器模型参数估计装置的实施例如图3所示，图3为本发明提供的配电变压器模型参数估计装置的实施例的示意图，该装置包括：
传递函数值测量值获得单元301，用于获取配电变压器的中低频频率响应信号，获得传递函数值的测量值，具体的，该传递函数值测量值获得单元可以用于向配电变压器的低压侧的其中一相绕组端子注入一个低振幅、可变频率的正弦波电压，具体的频率可以但不限于为1Hz至5Hz，并在高压侧的同相绕组端子处捕获响应信号，得到以中低频频率响应信号表示的传递函数值的测量值为：
其中，Vin为低压侧的输入电压幅值，Vout为高压侧的输出电压幅值；
传递函数值计算值计算单元302，用于绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算配电变压器的传递函数值的计算值，根据配电变压器的电磁传递特性，绘制出配电变压器的等效电路图，对图中的各个节点运用基尔霍夫定律，得到配电变压器的等效电路的状态空间微分方程为：
式中，
表示对变量vk取对时间的微分；
其中，GH0为高压绕组与接地箱体间的绝缘电导，CH0为高压绕组与接地箱体间的绝缘电容，GH为高压绕组电导，CH为高压绕组电容，RH为高压绕组电阻，LH为高压绕组电感，GHL为高低压绕组间的绝缘电导，CHL为高低压绕组间的绝缘电容，GL0为低压绕组与接地铁芯间的绝缘电导，CL0为低压绕组与接地铁芯间的绝缘电容，GL为低压绕组电导，CL为低压绕组电容，RL为低压绕组电阻，LL为低压绕组电感，Vj-1、Vj、Vj+1、Vj+2分别为节点j-1、j、j+1、j+2的电压，Ij-1、Ij、Ij+1分别为节点j-1、j、j+1流出的电流，Vi-1、Vi、Vi+1、Vi+2分别为节点i-1、i、i+1、i+2的电压，Ii-1、Ii、Ii+1分别为节点i-1、i、i+1流出的电流；
结合等效电路图和状态空间微分方程，在某一输入电压和电流频率下各元件参数确定后，利用MATLAB软件计算配电变压器的传递函数值的计算值；
模型参数估计单元303，用于使用遗传算法估计配电变压器的模型参数，需要说明的是，该遗传算法的核心在于：根据中低频输入下测量得到的配电变压器传递函数值的计算值，采用遗传算法估算各元件参数，代入至推导得到传递函数值计算程序中，使计算得到的传递函数值逐渐逼近匹配上述传递函数值的测量值，最终得到满足使用精度要求的变压器模型参数值。
在上述配电变压器模型参数估计装置的一个具体实施例中，模型参数估计单元具体可以用于：
遗传算法的目标函数为测量得到的传递函数值的测量值与传递函数值的计算值之间的误差最小，即
min J＝J1+J2
其中，J为遗传算法目标函数，J1和J2分别为全频率点子目标函数和谐振频率点子目标函数，(|TF|m)measured和(|TF|m)estimated分别为在频率点m处测量得到的传递函数值和估计的传递函数值，n为频率点的总数，(|TF|mr)measured和(|TF|mr)estimated分别为在谐振频率点mr处测量的传递函数值的测量值和传递函数值的计算值，z为谐振频率点总数。
具体而言，该算法的计算流程可以包括如下4个步骤：
1)初始化：创建配电变压器等效电路元件电气参数的随机初始种群；
2)评估：根据上一步定义的参数和图2的等效电路，计算估计的传递函数响应值，并与实际测量得到的响应值进行比较，评估上述公式中的目标函数；
3)更新：应用选择和繁殖算子(精英论、交叉和变异)创建一个新的种群，并将最优异的个体复制到新创建的种群中；
4)终止：重复上述更新过程，直到目标函数达到全局最小值或迭代次数达到最大值。
本发明提供的一种计算机设备的实施例如图4所示，图4为本发明提供的一种计算机设备的实施例的示意图，该设备包括：
存储器401，用于存储计算机程序；
处理器402，用于执行所述计算机程序时实现如上面任一项所述配电变压器模型参数估计方法的步骤。
本发明提供的一种计算机可读存储介质的实施例中，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上面任一项所述配电变压器模型参数估计方法的步骤。
利用上述装置、设备和存储介质，无需变压器本体参数即可方便的估计出等效电路的参数，便于配电变压器仿真建模和内部故障的识别。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.配电变压器模型参数估计方法，其特征在于，包括：
获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；
绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值；
使用遗传算法估计所述配电变压器的模型参数。
2.根据权利要求1所述的配电变压器模型参数估计方法，其特征在于，所述获取配电变压器的中低频频率响应信号，获得传递函数值的测量值包括：
向所述配电变压器的低压侧的其中一相绕组端子注入一个低振幅、可变频率的正弦波电压，并在高压侧的同相绕组端子处捕获响应信号，得到以中低频频率响应信号表示的所述传递函数值的测量值为：
其中，Vin为所述低压侧的输入电压幅值，Vout为所述高压侧的输出电压幅值。
3.根据权利要求1所述的配电变压器模型参数估计方法，其特征在于，所述绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值包括：
根据配电变压器的电磁传递特性，绘制出配电变压器的等效电路图，对图中的各个节点运用基尔霍夫定律，得到配电变压器的等效电路的状态空间微分方程为：
式中，
表示对变量vk取对时间的微分；
其中，GH0为高压绕组与接地箱体间的绝缘电导，CH0为高压绕组与接地箱体间的绝缘电容，GH为高压绕组电导，CH为高压绕组电容，RH为高压绕组电阻，LH为高压绕组电感，GHL为高低压绕组间的绝缘电导，CHL为高低压绕组间的绝缘电容，GL0为低压绕组与接地铁芯间的绝缘电导，CL0为低压绕组与接地铁芯间的绝缘电容，GL为低压绕组电导，CL为低压绕组电容，RL为低压绕组电阻，LL为低压绕组电感，Vj-1、Vj、Vj+1、Vj+2分别为节点j-1、j、j+1、j+2的电压，Ij-1、Ij、Ij+1分别为节点j-1、j、j+1流出的电流，Vi-1、Vi、Vi+1、Vi+2分别为节点i-1、i、i+1、i+2的电压，Ii-1、Ii、Ii+1分别为节点i-1、i、i+1流出的电流；
结合所述等效电路图和所述状态空间微分方程，在某一输入电压和电流频率下各元件参数确定后，利用MATLAB软件计算配电变压器的传递函数值的计算值。
4.根据权利要求1所述的配电变压器模型参数估计方法，其特征在于，所述使用遗传算法估计所述配电变压器的模型参数包括：
所述遗传算法的目标函数为测量得到的传递函数值的测量值与传递函数值的计算值之间的误差最小，即
min J＝J1+J2
其中，J为遗传算法目标函数，J1和J2分别为全频率点子目标函数和谐振频率点子目标函数，(|TF|m)measured和(|TF|m)estimated分别为在频率点m处测量得到的传递函数值和估计的传递函数值，n为频率点的总数，(|TF|mr)measured和(|TF|mr)estimated分别为在谐振频率点mr处测量的传递函数值的测量值和传递函数值的计算值，z为谐振频率点总数。
5.配电变压器模型参数估计装置，其特征在于，包括：
传递函数值测量值获得单元，用于获取配电变压器的中低频频率响应信号，获得传递函数值的测量值；
传递函数值计算值计算单元，用于绘制配电变压器的等效电路，结合配电变压器状态空间微分方程，计算所述配电变压器的传递函数值的计算值；
模型参数估计单元，用于使用遗传算法估计所述配电变压器的模型参数。
6.根据权利要求5所述的配电变压器模型参数估计装置，其特征在于，所述传递函数值测量值获得单元具体用于：
向所述配电变压器的低压侧的其中一相绕组端子注入一个低振幅、可变频率的正弦波电压，并在高压侧的同相绕组端子处捕获响应信号，得到以中低频频率响应信号表示的所述传递函数值的测量值为：
其中，Vin为所述低压侧的输入电压幅值，Vout为所述高压侧的输出电压幅值。
7.根据权利要求5所述的配电变压器模型参数估计方法，其特征在于，所述传递函数值计算值计算单元具体用于：
根据配电变压器的电磁传递特性，绘制出配电变压器的等效电路图，对图中的各个节点运用基尔霍夫定律，得到配电变压器的等效电路的状态空间微分方程为：
式中，
表示对变量vk取对时间的微分；
其中，GH0为高压绕组与接地箱体间的绝缘电导，CH0为高压绕组与接地箱体间的绝缘电容，GH为高压绕组电导，CH为高压绕组电容，RH为高压绕组电阻，LH为高压绕组电感，GHL为高低压绕组间的绝缘电导，CHL为高低压绕组间的绝缘电容，GL0为低压绕组与接地铁芯间的绝缘电导，CL0为低压绕组与接地铁芯间的绝缘电容，GL为低压绕组电导，CL为低压绕组电容，RL为低压绕组电阻，LL为低压绕组电感，Vj-1、Vj、Vj+1、Vj+2分别为节点j-1、j、j+1、j+2的电压，Ij-1、Ij、Ij+1分别为节点j-1、j、j+1流出的电流，Vi-1、Vi、Vi+1、Vi+2分别为节点i-1、i、i+1、i+2的电压，Ii-1、Ii、Ii+1分别为节点i-1、i、i+1流出的电流；
结合所述等效电路图和所述状态空间微分方程，在某一输入电压和电流频率下各元件参数确定后，利用MATLAB软件计算配电变压器的传递函数值的计算值。
8.根据权利要求5所述的配电变压器模型参数估计装置，其特征在于，所述模型参数估计单元具体用于：
所述遗传算法的目标函数为测量得到的传递函数值的测量值与传递函数值的计算值之间的误差最小，即
min J＝J1+J2
其中，J为遗传算法目标函数，J1和J2分别为全频率点子目标函数和谐振频率点子目标函数，(|TF|m)measured和(|TF|m)estimated分别为在频率点m处测量得到的传递函数值和估计的传递函数值，n为频率点的总数，(|TF|mr)measured和(|TF|mr)estimated分别为在谐振频率点mr处测量的传递函数值的测量值和传递函数值的计算值，z为谐振频率点总数。
9.一种计算机设备，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至4任一项所述配电变压器模型参数估计方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至4任一项所述配电变压器模型参数估计方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110255432.7A CN112861469A (zh) | 2021-03-09 | 2021-03-09 | 配电变压器模型参数估计方法、装置、设备和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110255432.7A CN112861469A (zh) | 2021-03-09 | 2021-03-09 | 配电变压器模型参数估计方法、装置、设备和存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112861469A true CN112861469A (zh) | 2021-05-28 |
Family
ID=75995082
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110255432.7A Pending CN112861469A (zh) | 2021-03-09 | 2021-03-09 | 配电变压器模型参数估计方法、装置、设备和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112861469A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018010028A1 (fr) * | 2016-07-14 | 2018-01-18 | HYDRO-QUéBEC | Détection d'anomalies dans un réseau électrique |
| CN108920751A (zh) * | 2018-05-24 | 2018-11-30 | 西安交通大学 | 基于拓扑优化的电力变压器绕组变形状态逆向求解方法 |
| CN109061407A (zh) * | 2018-06-28 | 2018-12-21 | 云南电网有限责任公司保山供电局 | 基于FDS的ZnO阀片等效电路模型的参数估计方法 |
| CN110826018A (zh) * | 2020-01-09 | 2020-02-21 | 四川大学 | 基于有理函数逼近的一次信号恢复方法及系统 |
| CN110826255A (zh) * | 2019-12-02 | 2020-02-21 | 国网河北省电力有限公司电力科学研究院 | 一种电力装备中绕组线圈部件的建模方法 |
-
2021
- 2021-03-09 CN CN202110255432.7A patent/CN112861469A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018010028A1 (fr) * | 2016-07-14 | 2018-01-18 | HYDRO-QUéBEC | Détection d'anomalies dans un réseau électrique |
| CN108920751A (zh) * | 2018-05-24 | 2018-11-30 | 西安交通大学 | 基于拓扑优化的电力变压器绕组变形状态逆向求解方法 |
| CN109061407A (zh) * | 2018-06-28 | 2018-12-21 | 云南电网有限责任公司保山供电局 | 基于FDS的ZnO阀片等效电路模型的参数估计方法 |
| CN110826255A (zh) * | 2019-12-02 | 2020-02-21 | 国网河北省电力有限公司电力科学研究院 | 一种电力装备中绕组线圈部件的建模方法 |
| CN110826018A (zh) * | 2020-01-09 | 2020-02-21 | 四川大学 | 基于有理函数逼近的一次信号恢复方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Abu-Siada et al. | Estimating power transformer high frequency model parameters using frequency response analysis | |
| Rahimpour et al. | Transfer function method to diagnose axial displacement and radial deformation of transformer windings | |
| CN112485509B (zh) | 基于非线性宽频模型的暂态过电压测量装置与方法 | |
| CN110175351B (zh) | 一种换流变压器建模方法 | |
| Ragavan et al. | Construction of physically realizable driving-point function from measured frequency response data on a model winding | |
| CN117233687B (zh) | 一种基于历史数据的cvt初始误差评估方法、介质及终端 | |
| CN108490379B (zh) | 一种基于自激振荡波的变压器绕组波过程校验方法 | |
| Samarawickrama et al. | Impulse generator optimum setup for transient testing of transformers using frequency-response analysis and genetic algorithm | |
| Karimifard et al. | Determination of axial displacement extent based on transformer winding transfer function estimation using vector‐fitting method | |
| Karimifard et al. | Localization of winding radial deformation and determination of deformation extent using vector fitting‐based estimated transfer function | |
| CN110826255B (zh) | 一种电力装备中绕组线圈部件的建模方法 | |
| Orosz et al. | Modeling the high frequency behavior of the Rogowski-coil passive L/r integrator current transducer with analytical and finite element method | |
| Bagheri et al. | Transformer winding parameter identification based on frequency response analysis using hybrid wavelet transform (WT) and simulated annealing algorithm (SA) and compare with genetic algorithm (GA) | |
| CN112861469A (zh) | 配电变压器模型参数估计方法、装置、设备和存储介质 | |
| Høidalen et al. | Analysis of gray box modelling of transformers | |
| Abu-Siada | High frequency transformer modelling using state space representation for FRA studies | |
| Mwaniki et al. | Characterizing power transformer frequency responses using bipolar pseudo-random current impulses | |
| Mikulović et al. | Calculation of transients in transformer winding and determination of winding parameters | |
| Albert et al. | Comparing two topology transformer hysteresis models with power transformer measurements | |
| Stevanović et al. | Behavioral circuit modeling of single-and three-phase chokes with multi-resonances | |
| CN113496096A (zh) | 避雷器暂态电压获取方法、装置、计算机设备和存储介质 | |
| Saied | New solution technique for the frequency and transient response of transformer windings with all inter-turn mutual inductances and capacitances included | |
| Yu et al. | A circuit model in a wide frequency range for power transformer and analysis of its characteristics | |
| Bigdeli et al. | Optimized modeling of transformer in transient state with genetic algorithm | |
| Januário et al. | Modelling of the impedance frequency response of transformers at no load and under nominal conditions |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230425 Address after: 310018 No. 11 street, Hangzhou economic and Technological Development Zone, Hangzhou, Zhejiang 91 Applicant after: HANGZHOU ELECTRIC EQUIPMENT MANUFACTURING Co.,Ltd. Applicant after: State Grid Zhejiang Electric Power Co., Ltd. Hangzhou Yuhang District Power Supply Co. Applicant after: HANGZHOU POWER SUPPLY COMPANY, STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd. Address before: 310018 No. 11 street, Hangzhou economic and Technological Development Zone, Zhejiang 91 Applicant before: HANGZHOU ELECTRIC EQUIPMENT MANUFACTURING Co.,Ltd. Applicant before: STATE GRID ZHEJIANG HANGZHOU YUHANG POWER SUPPLY Co. Applicant before: HANGZHOU POWER SUPPLY COMPANY, STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd. |
|
| TA01 | Transfer of patent application right |