CN110910001A - Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system - Google Patents

Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system Download PDF

Info

Publication number
CN110910001A
CN110910001A CN201911120033.9A CN201911120033A CN110910001A CN 110910001 A CN110910001 A CN 110910001A CN 201911120033 A CN201911120033 A CN 201911120033A CN 110910001 A CN110910001 A CN 110910001A
Authority
CN
China
Prior art keywords
transformer
voltage winding
parameter identification
target
sequence
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
CN201911120033.9A
Other languages
Chinese (zh)
Other versions
CN110910001B (en
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 Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power 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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201911120033.9A priority Critical patent/CN110910001B/en
Publication of CN110910001A publication Critical patent/CN110910001A/en
Application granted granted Critical
Publication of CN110910001B publication Critical patent/CN110910001B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a transformer parameter identification method, a system and a medium based on a wide-area synchronous phasor measurement system, wherein the transformer parameter identification method comprises the steps of establishing a transformer transformation ratio and impedance parameter identification model aiming at a target transformer; performing synchronous phasor measurement through a wide-area synchronous phasor measurement system to obtain synchronous measurement phasor of the target transformer; and combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method. Different from the traditional transformer on-line monitoring, the invention does not need to additionally install a transformer state sensor or monitoring equipment, but can utilize the measured data of the wide-area synchronous phasor measurement installed in the transformer substation with the voltage class of 220kV or above to identify the parameters of the transformer through data analysis, and is also suitable for the identification of the wide-area synchronous phasor measurement system installed in the transformer substation with the voltage class of below 220kV for the parameters of the transformer.

Description

Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system
Technical Field
The invention relates to the technical field of on-line monitoring of transformers, in particular to a transformer parameter identification method, a system and a medium based on a wide-area synchronous phasor measurement system.
Background
The transformer is an important component of the power system, plays a pivotal role in electric energy transmission, and has a serious influence on the power supply reliability and the system operation stability of the power system due to faults. How to timely and accurately discover and timely eliminate potential defects of a transformer, eliminate potential safety hazards of the transformer, improve the utilization rate of equipment and prolong the service life of the transformer is an important subject of power system research. In order to find out faults and hidden dangers of the transformer in time and accurately master the working state and the running condition of the transformer, the state of the transformer needs to be monitored and tracked in time. The methods for monitoring and tracking the state of the transformer are mainly classified into the following three methods: online monitoring, power failure test and live detection.
On-line monitoring is to collect the data of the equipment operation under the condition that the monitoring device and the equipment normally operate, to continuously or regularly monitor the working state of the equipment, generally needs to install a special on-line monitoring sensor and a special on-line monitoring device, and has large investment. The power failure test has been the important means that the electric power industry carries out fault detection all the time, but the routine test project of power failure is mostly low-voltage, and the test result can not reflect the state of equipment under actual voltage operation completely, and there is certain deviation and defect in the data to the transformer can age gradually at the operation in-process, has some potential trouble to discover in the power failure test, has delayed the electric power and has overhauld. In addition, power failure maintenance needs to stop the equipment, and the operation of the power system is influenced. The live detection is to detect live-line running electrical equipment in a short time through a special test instrument and an instrument device, and can detect and find potential faults without power failure. The three methods have advantages and disadvantages, and the online monitoring generally needs to install a special sensor and monitoring equipment, so that the investment is large; the power failure test will affect the normal operation of the power system; although the live detection does not need power failure, the live detection only provides a discontinuous detection means, and the running state of power equipment such as a transformer cannot be grasped in time.
Disclosure of Invention
The technical problems to be solved by the invention are as follows: the invention provides a transformer parameter identification method, a system and a medium based on a wide area synchronous phasor measurement system, which are different from the traditional transformer on-line monitoring, and the method does not need to additionally install a transformer state sensor or monitoring equipment, but can utilize the measurement data of the wide area synchronous phasor measurement installed in a transformer substation with the voltage class of 220kV or above to identify the parameters of the transformer through data analysis, and the method is also suitable for the identification of the wide area synchronous phasor measurement system installed in the transformer substation with the voltage class of below 220kV for the parameters of the transformer.
In order to solve the technical problems, the invention adopts the technical scheme that:
a transformer parameter identification method based on a wide area synchronous phasor measurement system comprises the following implementation steps:
1) establishing a parameter identification model of transformer transformation ratio and impedance aiming at a target transformer;
2) performing synchronous phasor measurement through a wide-area synchronous phasor measurement system to obtain synchronous measurement phasor of the target transformer;
3) and combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
Optionally, the target transformer is a two-winding transformer, and the functional expression of the parameter identification model of the transformer transformation ratio and the impedance established in step 1) is as follows:
Figure BDA0002275195830000021
Figure BDA0002275195830000022
in the above formula, K is the transformer transformation ratio of the target transformer,
Figure BDA0002275195830000023
impedance of the target transformer
Figure BDA0002275195830000024
The sum of impedance parameters after the primary winding and the secondary winding of the target transformer are converted to the primary side,
Figure BDA0002275195830000025
and
Figure BDA0002275195830000026
the currents of the primary side and the secondary side of the target transformer are respectively,
Figure BDA0002275195830000027
and
Figure BDA0002275195830000028
the voltages of the primary side and the secondary side of the target transformer are respectively.
Optionally, the synchronous measurement phasor obtained in step 2) includes obtaining currents of the primary side and the secondary side of the target transformer by synchronous phasor measurement
Figure BDA0002275195830000029
And
Figure BDA00022751958300000210
and the primary and secondary voltages of the target transformer
Figure BDA00022751958300000211
And
Figure BDA00022751958300000212
optionally, the detailed steps of step 3) include:
3.1A) setting an intermediate matrix variable according to the parameter identification model:
Figure BDA00022751958300000213
in the above formula, K is the transformer transformation ratio of the target transformer,
Figure BDA00022751958300000214
is the impedance of the target transformer and,
Figure BDA00022751958300000215
and
Figure BDA00022751958300000216
the currents of the primary side and the secondary side of the target transformer are respectively,
Figure BDA00022751958300000217
and
Figure BDA00022751958300000218
voltage of primary and secondary sides of the target transformer, ξ1And ξ2Errors of a primary side and a secondary side of the target transformer are respectively obtained; and establishing a matrix equation shown as the following formula according to the parameter identification model:
yk＝xkθ+ξk
3.2A) forming a phasor sequence from the obtained synchronous measured phasors
Figure BDA00022751958300000219
Wherein k is the sequence number, k is 1,2, …, n, and n is the maximum value of the sequence number; establishing a matrix equation as shown in the following formula:
Figure BDA00022751958300000220
in the above formula, y1～ynIs ykSequence of (a) x1～xnIs ykSequence of (iii) ξ1～ξnIs ξkThe sequence of (a);
3.3A) the matrix equation set up in step 3.2A) is expressed as y ═ Φ θ + ξ, and a least squares estimate of the matrix variable θ is obtained for this equation using the least squares method
Figure BDA00022751958300000221
To obtain
Figure BDA00022751958300000222
Thereby obtaining the transformation ratio K and the impedance of the transformer
Figure BDA00022751958300000223
The best estimate of (c); wherein y represents the sequence y1～ynThe constructed matrix, phi, representing the sequence x1～xnFormed matrix, ξ denotes the sequence ξ1～ξnA matrix is formed.
Optionally, the target transformer is a three-winding transformer, and the low-voltage winding tap transformation ratio K of the three-winding transformer3Known, medium voltage winding resistance R2Reactance X2Impedance, impedance
Figure BDA0002275195830000031
As known, the functional expression of the parameter identification model for the transformer transformation ratio and the impedance established in step 1) is shown as follows:
Figure BDA0002275195830000032
Figure BDA0002275195830000033
in the above formula, K1，K2，K3Respectively the transformation ratios of a high-voltage winding, a medium-voltage winding and a low-voltage winding of a target transformer,
Figure BDA0002275195830000034
respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,
Figure BDA0002275195830000035
voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,
Figure BDA0002275195830000036
current phasors, I, of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer respectively1，I2，I3Are respectively the current phasors
Figure BDA0002275195830000037
The amplitude of (c).
Optionally, the step 2) of obtaining the synchronous measurement phasor includes obtaining voltages of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer through the synchronous phasor measurement
Figure BDA0002275195830000038
And the current of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer
Figure BDA0002275195830000039
Optionally, the detailed steps of step 3) include:
3.1B) setting an intermediate matrix variable according to the parameter identification model:
Figure BDA00022751958300000310
Figure BDA00022751958300000311
in the above formula, K1，K2，K3Respectively the transformation ratios of a high-voltage winding, a medium-voltage winding and a low-voltage winding of a target transformer,
Figure BDA00022751958300000312
respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,
Figure BDA00022751958300000313
voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,
Figure BDA00022751958300000314
current phasors, I, of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer respectively1，I2，I3Are respectively the current phasors
Figure BDA00022751958300000315
The amplitude of (a) of (b) is,
Figure BDA00022751958300000316
representing the phasor of the current
Figure BDA00022751958300000317
And
Figure BDA00022751958300000318
the difference in the angle between the two,
Figure BDA00022751958300000319
representing the phasor of the current
Figure BDA00022751958300000320
And
Figure BDA00022751958300000321
the difference in the angle between the two,
Figure BDA00022751958300000322
representing the phasor of the current
Figure BDA00022751958300000323
And
Figure BDA00022751958300000324
the difference in the angle between the two,
Figure BDA00022751958300000325
representing the phasor of the current
Figure BDA00022751958300000326
And
Figure BDA00022751958300000327
angle difference therebetween, ξ1～ξ4Respectively represents K1，K2，
Figure BDA00022751958300000328
Error corresponding to the estimated value; according to the parameter identification model, a matrix equation shown as the following formula is established:
yk＝xkθ+ξk
3.2B) forming a phasor sequence from the obtained synchronous measurement phasors
Figure BDA00022751958300000329
Wherein k is the sequence number, k is 1,2, …, n, and n is the maximum value of the sequence number; establishing a matrix equation as shown in the following formula:
Figure BDA0002275195830000041
in the above formula, y1～ynIs ykSequence of (a) x1～xnIs ykSequence of (iii) ξ1～ξnIs ξkThe sequence of (a);
3.3B) expressing the matrix equation set up in step 3.2A) as y ═ Φ θ + ξ, and using the least squares method to solve the least squares estimate of the matrix variable θ for this equation
Figure BDA0002275195830000042
To obtain
Figure BDA0002275195830000043
Thereby obtaining the impedance of the primary winding of the target transformer
Figure BDA0002275195830000044
Impedance of low voltage winding
Figure BDA0002275195830000045
And the transformation ratio K of the primary winding of the target transformer1Transformation ratio K of medium voltage winding2The best estimate of (c); wherein y represents the sequence y1～ynThe constructed matrix, phi, representing the sequence x1～xnFormed matrix, ξ denotes the sequence ξ1～ξnA matrix is formed.
In addition, the invention also provides a transformer parameter identification system based on the wide-area synchronous phasor measurement system, which comprises the following steps:
the model initialization program unit is used for establishing a parameter identification model of the transformer transformation ratio and the impedance aiming at the target transformer;
the synchronous measurement phasor acquisition program unit is used for carrying out synchronous phasor measurement through the wide-area synchronous phasor measurement system to obtain the synchronous measurement phasor of the target transformer;
and the transformer parameter solving program unit is used for combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
In addition, the present invention also provides a transformer parameter identification system based on a wide area synchrophasor measurement system, which includes a computer device programmed or configured to execute the steps of the transformer parameter identification method based on a wide area synchrophasor measurement system, or a computer program programmed or configured to execute the transformer parameter identification method based on a wide area synchrophasor measurement system is stored on a memory of the computer device.
In addition, the present invention also provides a computer readable storage medium having a computer program stored thereon that is programmed or configured to perform the transformer parameter identification method based on the wide area synchronous phasor measurement system.
Compared with the prior art, the invention has the following advantages:
1. different from the traditional transformer on-line monitoring, the invention does not need to additionally install a transformer state sensor or monitoring equipment, but can utilize the measured data of the wide-area synchronous phasor measurement already installed in the transformer substation with the voltage class of 220kV or above to identify the parameters of the transformer through data analysis.
2. The method is also suitable for identifying the transformer parameters by the wide-area synchronous phasor measurement system installed in the transformer substation with the voltage class below 220 kV.
3. The invention combines the phase sequence with the parameter identification model and solves the transformer parameter by adopting the least square method, thereby further realizing the early warning treatment of the identification result of the transformer parameter according to the requirement, for example, once the identified parameter is found to have continuous abnormity, the early warning is sent out in time.
Drawings
FIG. 1 is a schematic diagram of the basic flow of the process of the present invention.
Fig. 2 is an equivalent model of a two-winding transformer according to an embodiment of the invention.
Fig. 3 is a simplified model of a two-winding transformer according to an embodiment of the present invention.
Fig. 4 is an equivalent model of a three-winding transformer according to a second embodiment of the present invention.
Fig. 5 is a simplified model of a three-winding transformer according to a second embodiment of the present invention.
Detailed Description
The first embodiment is as follows:
as shown in fig. 1, the implementation steps of the transformer parameter identification method based on the wide area synchronous phasor measurement system in the embodiment include:
1) establishing a parameter identification model of transformer transformation ratio and impedance aiming at a target transformer;
2) performing synchronous phasor Measurement through a Wide-area Measurement System (WAMS) to obtain a synchronous Measurement phasor of the target transformer;
3) and combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
The common transformers are divided into two-winding transformers and three-winding transformers, and the target transformer in this embodiment is a two-winding transformer, and its equivalent model is shown in fig. 2, wherein
Figure BDA0002275195830000051
The voltage and the current of the primary side and the secondary side of the two-winding transformer are respectively,
Figure BDA0002275195830000052
respectively converted into the voltage and current values of the primary side, Gm+jBmFor the admittance parameter of the excitation branch,
Figure BDA0002275195830000053
in order to be the exciting current,
Figure BDA0002275195830000054
and K is the sum of impedance parameters after the primary winding and the secondary winding are converted to the primary side, and is the transformation ratio of the transformer. The admittance to ground of the transformer excitation branch is considered to be generally small, so that identification is not needed. Regardless of the excitation branch, the pi-type equivalent circuit of the two-winding transformer is further simplified as shown in fig. 3, so that the following relationships between variables can be obtained:
Figure BDA0002275195830000055
Figure BDA0002275195830000056
on the basis, a parameter identification model for establishing the transformer transformation ratio and the impedance can be deduced. In this embodiment, the functional expression of the parameter identification model of the transformer transformation ratio and the impedance established in step 1) is shown as follows:
Figure BDA0002275195830000057
Figure BDA0002275195830000058
in the above formula, K is the transformer transformation ratio of the target transformer,
Figure BDA0002275195830000059
impedance of the target transformer
Figure BDA00022751958300000510
The sum of impedance parameters after the primary winding and the secondary winding of the target transformer are converted to the primary side,
Figure BDA00022751958300000511
and
Figure BDA00022751958300000512
the currents of the primary side and the secondary side of the target transformer are respectively,
Figure BDA00022751958300000513
and
Figure BDA00022751958300000514
the voltages of the primary side and the secondary side of the target transformer are respectively.
(Wide-area Measurement System,WAMS)
In this embodiment, the synchronous phasor measurement obtained in step 2) includes synchronous phasor measurement to obtain currents of the primary side and the secondary side of the target transformer
Figure BDA00022751958300000515
And
Figure BDA00022751958300000516
and the primary and secondary voltages of the target transformer
Figure BDA00022751958300000517
And
Figure BDA00022751958300000518
in this embodiment, the detailed steps of step 3) include:
3.1A) setting an intermediate matrix variable according to the parameter identification model:
Figure BDA0002275195830000061
in the above formula, K is the transformer transformation ratio of the target transformer,
Figure BDA0002275195830000062
is the impedance of the target transformer and,
Figure BDA0002275195830000063
and
Figure BDA0002275195830000064
the currents of the primary side and the secondary side of the target transformer are respectively,
Figure BDA0002275195830000065
and
Figure BDA0002275195830000066
voltage of primary and secondary sides of the target transformer, ξ1And ξ2Errors of a primary side and a secondary side of the target transformer are respectively obtained; and establishing a matrix equation shown as the following formula according to the parameter identification model:
yk＝xkθ+ξk
3.2A) forming a phasor sequence from the obtained synchronous measured phasors
Figure BDA0002275195830000067
Wherein k is the sequence number, k is 1,2, …, n, and n is the maximum value of the sequence number; establishing a matrix equation as shown in the following formula:
Figure BDA0002275195830000068
in the above formula, y1～ynIs ykSequence of (a) x1～xnIs ykSequence of (iii) ξ1～ξnIs ξkThe sequence of (a);
3.3A) the matrix equation set up in step 3.2A) is expressed as y ═ Φ θ + ξ, and a least squares estimate of the matrix variable θ is obtained for this equation using the least squares method
Figure BDA0002275195830000069
To obtain
Figure BDA00022751958300000610
Thereby obtaining the transformation ratio K and the impedance of the transformer
Figure BDA00022751958300000611
The best estimate of (c); wherein y represents the sequence y1～ynThe constructed matrix, phi, representing the sequence x1～xnFormed matrix, ξ denotes the sequence ξ1～ξnA matrix is formed.
On the basis of combining the phase sequence with the parameter identification model and solving the transformer parameters by adopting a least square method, the early warning treatment can be further realized according to the identification result of the transformer parameters, for example, once the identified parameters are found to be continuous abnormal, the early warning is timely sent out. For example, if the identified parameter is abnormal (e.g. deviated from the normal value by more than 20%) for a long time (e.g. two consecutive hours), an alarm message that the transformer equipment may be abnormal is issued. Further, the equipment abnormity is required to be checked and eliminated through live detection or power failure overhaul.
In addition, this embodiment still provides a transformer parameter identification system based on wide area synchrophasor measurement system, includes:
the model initialization program unit is used for establishing a parameter identification model of the transformer transformation ratio and the impedance aiming at the target transformer;
the synchronous measurement phasor acquisition program unit is used for carrying out synchronous phasor measurement through the wide-area synchronous phasor measurement system to obtain the synchronous measurement phasor of the target transformer;
and the transformer parameter solving program unit is used for combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
In addition, the present embodiment further provides a transformer parameter identification system based on the wan synchrophasor measurement system, which includes a computer device programmed or configured to execute the steps of the transformer parameter identification method based on the wan synchrophasor measurement system according to the present embodiment, or a computer program programmed or configured to execute the transformer parameter identification method based on the wan synchrophasor measurement system according to the present embodiment is stored in a memory of the computer device. In addition, the present embodiment also provides a computer readable storage medium, which stores a computer program programmed or configured to execute the transformer parameter identification method based on the wide area synchronous phasor measurement system according to the present embodiment.
Example two:
the present embodiment is the same as the basic method of the first embodiment, and the main difference is that the target transformer in the present embodiment is a three-winding transformer, and the specific implementation manner of combining the phase sequence with the parameter identification model and solving the transformer parameter by using the least square method in the corresponding step 3) is different. Needless to say, those skilled in the art can also apply the transformer parameter identification method, system and medium based on the wide area synchronous phasor measurement system of the present invention to transformers with more windings according to the first embodiment and the technical teaching of this embodiment, which are not described herein again.
In this embodiment, the target transformer is a three-winding transformer, and the low-voltage winding tap transformation ratio K of the three-winding transformer3Known, medium voltage winding resistance R2Reactance X2Impedance, impedance
Figure BDA0002275195830000071
Known as such, in step 1)The functional expression of the parameter identification model of the transformer transformation ratio and the impedance is established as follows:
Figure BDA0002275195830000072
Figure BDA0002275195830000073
in the above formula, K1，K2，K3Respectively the transformation ratios of a high-voltage winding, a medium-voltage winding and a low-voltage winding of a target transformer,
Figure BDA0002275195830000074
respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,
Figure BDA0002275195830000075
voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,
Figure BDA0002275195830000076
current phasors, I, of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer respectively1，I2，I3Are respectively the current phasors
Figure BDA0002275195830000077
The amplitude of (c).
In this embodiment, the synchronous phasor measurement obtained in step 2) includes synchronous phasor measurement to obtain voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer
Figure BDA0002275195830000078
And the current phasor of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer
Figure BDA0002275195830000079
In this embodiment, the detailed steps of step 3) include:
3.1B) setting an intermediate matrix variable according to the parameter identification model:
Figure BDA00022751958300000710
Figure BDA0002275195830000081
in the above formula, K1，K2，K3Respectively the transformation ratios of a high-voltage winding, a medium-voltage winding and a low-voltage winding of a target transformer,
Figure BDA0002275195830000082
respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,
Figure BDA0002275195830000083
voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,
Figure BDA0002275195830000084
current phasors, I, of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer respectively1，I2，I3Are respectively the current phasors
Figure BDA0002275195830000085
The amplitude of (a) of (b) is,
Figure BDA0002275195830000086
representing the phasor of the current
Figure BDA0002275195830000087
And
Figure BDA0002275195830000088
the difference in the angle between the two,
Figure BDA0002275195830000089
representing the phasor of the current
Figure BDA00022751958300000810
And
Figure BDA00022751958300000811
the difference in the angle between the two,
Figure BDA00022751958300000812
representing the phasor of the current
Figure BDA00022751958300000813
And
Figure BDA00022751958300000814
the difference in the angle between the two,
Figure BDA00022751958300000815
representing the phasor of the current
Figure BDA00022751958300000816
And
Figure BDA00022751958300000817
angle difference therebetween, ξ1～ξ4Respectively represents K1，K2，
Figure BDA00022751958300000818
Error corresponding to the estimated value; and establishing a matrix equation shown as the following formula according to the parameter identification model:
yk＝xkθ+ξk
3.2B) forming a phasor sequence from the obtained synchronous measurement phasors
Figure BDA00022751958300000819
Wherein k is the sequence number, k is 1,2, …, n, and n is the maximum value of the sequence number; establishing a matrix equation as shown in the following formula:
Figure BDA00022751958300000820
in the above formula, y1～ynIs ykSequence of (a) x1～xnIs ykSequence of (iii) ξ1～ξnIs ξkThe sequence of (a);
3.3B) expressing the matrix equation set up in step 3.2A) as y ═ Φ θ + ξ, and using the least squares method to solve the least squares estimate of the matrix variable θ for this equation
Figure BDA00022751958300000821
To obtain
Figure BDA00022751958300000822
Thereby obtaining the impedance of the primary winding of the target transformer
Figure BDA00022751958300000823
Impedance of low voltage winding
Figure BDA00022751958300000824
And the transformation ratio K of the primary winding of the target transformer1Transformation ratio K of medium voltage winding2The best estimate of (c); wherein y represents the sequence y1～ynThe constructed matrix, phi, representing the sequence x1～xnFormed matrix, ξ denotes the sequence ξ1～ξnA matrix is formed.
In addition, this embodiment still provides a transformer parameter identification system based on wide area synchrophasor measurement system, includes:
the model initialization program unit is used for establishing a parameter identification model of the transformer transformation ratio and the impedance aiming at the target transformer;
the synchronous measurement phasor acquisition program unit is used for carrying out synchronous phasor measurement through the wide-area synchronous phasor measurement system to obtain the synchronous measurement phasor of the target transformer;
and the transformer parameter solving program unit is used for combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
In addition, the present embodiment further provides a transformer parameter identification system based on the wan synchrophasor measurement system, which includes a computer device programmed or configured to execute the steps of the transformer parameter identification method based on the wan synchrophasor measurement system according to the present embodiment, or a computer program programmed or configured to execute the transformer parameter identification method based on the wan synchrophasor measurement system according to the present embodiment is stored in a memory of the computer device. In addition, the present embodiment also provides a computer readable storage medium, which stores a computer program programmed or configured to execute the transformer parameter identification method based on the wide area synchronous phasor measurement system according to the present embodiment.
In summary, the present invention provides a transformer parameter identification method based on a wide-area synchronous phasor measurement system, which implements real-time dynamic identification of impedance and transformation ratio parameters of two-winding and three-winding transformers by real-time measurement data of the wide-area synchronous phasor measurement system, including current and voltage value synchronous measurement values. According to the invention, the early warning of the abnormal state of the transformer can be realized by monitoring the impedance parameter and the change of the transformation ratio parameter of the transformer in real time, and a new perception means is provided for the online state monitoring and evaluation of the transformer. The invention utilizes the synchronous measurement data of external voltage and current of transformation operation acquired by the wide-area synchronous phasor system installed in the transformer substation with the voltage class of 220kV or above, and identifies the parameters of impedance, transformation ratio and the like of the transformer in real time by a data analysis method.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention. It should be noted that modifications and embellishments within the scope of the invention may occur to those skilled in the art without departing from the principle of the invention, and are considered to be within the scope of the invention.

Claims (10)

1. A transformer parameter identification method based on a wide area synchronous phasor measurement system is characterized by comprising the following implementation steps:
1) establishing a parameter identification model of transformer transformation ratio and impedance aiming at a target transformer;
2) performing synchronous phasor measurement through a wide-area synchronous phasor measurement system to obtain synchronous measurement phasor of the target transformer;
3) and combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
2. The method according to claim 1, wherein the target transformer is a two-winding transformer, and the functional expression of the parameter identification model of transformer transformation ratio and impedance established in step 1) is as follows:
Figure FDA0002275195820000011
Figure FDA0002275195820000012
in the above formula, K is the transformer transformation ratio of the target transformer,
Figure FDA0002275195820000013
impedance of the target transformer
Figure FDA0002275195820000014
The sum of impedance parameters after the primary winding and the secondary winding of the target transformer are converted to the primary side,
Figure FDA0002275195820000015
and
Figure FDA0002275195820000016
the currents of the primary side and the secondary side of the target transformer are respectively,
Figure FDA0002275195820000017
and
Figure FDA0002275195820000018
the voltages of the primary side and the secondary side of the target transformer are respectively.
3. The method as claimed in claim 2, wherein the step of obtaining the measured phasor in step 2) comprises obtaining primary and secondary currents of the target transformer by the measured phasor
Figure FDA0002275195820000019
And
Figure FDA00022751958200000110
and the primary and secondary voltages of the target transformer
Figure FDA00022751958200000111
And
Figure FDA00022751958200000112
4. the method for identifying transformer parameters based on the wide-area synchronous phasor measurement system according to claim 2, wherein the detailed steps of step 3) include:
3.1A) setting an intermediate matrix variable according to the parameter identification model:
Figure FDA00022751958200000113
in the above formula, K is the transformer transformation ratio of the target transformer,
Figure FDA00022751958200000114
is the impedance of the target transformer and,
Figure FDA00022751958200000115
and
Figure FDA00022751958200000116
the currents of the primary side and the secondary side of the target transformer are respectively,
Figure FDA00022751958200000117
and
Figure FDA00022751958200000118
voltage of primary and secondary sides of the target transformer, ξ1And ξ2Errors of a primary side and a secondary side of the target transformer are respectively obtained; and establishing a matrix equation as shown in the following formula according to the parameter identification model
yk＝xkθ+ξk
3.2A) forming a phasor sequence from the obtained synchronous measured phasors
Figure FDA00022751958200000119
Wherein k is the sequence number, k is 1,2, …, n, and n is the maximum value of the sequence number; establishing a matrix equation as shown in the following formula:
Figure FDA00022751958200000120
in the above formula, y1～ynIs ykSequence of (a) x1～xnIs ykSequence of (iii) ξ1～ξnIs ξkThe sequence of (a);
3.3A) the matrix equation set up in step 3.2A) is expressed as y ═ Φ θ + ξ, and a least squares estimate of the matrix variable θ is obtained for this equation using the least squares method
Figure FDA0002275195820000021
To obtain
Figure FDA0002275195820000022
Thereby obtaining the transformation ratio K and the impedance of the transformer
Figure FDA0002275195820000023
The best estimate of (c); wherein y represents the sequence y1～ynThe constructed matrix, phi, representing the sequence x1～xnFormed matrix, ξ denotes the sequence ξ1～ξnA matrix is formed.
5. The method as claimed in claim 1, wherein the target transformer is a three-winding transformer, and the low-voltage winding tap transformation ratio K of the three-winding transformer is defined as K3Known, medium voltage winding resistance R2Reactance X2Impedance, impedance
Figure FDA0002275195820000024
As known, the functional expression of the parameter identification model for the transformer transformation ratio and the impedance established in step 1) is shown as follows:
Figure FDA0002275195820000025
Figure FDA0002275195820000026
in the above formula, K1，K2，K3Respectively the transformation ratios of a high-voltage winding, a medium-voltage winding and a low-voltage winding of a target transformer,
Figure FDA0002275195820000027
respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,
Figure FDA0002275195820000028
respectively a high-voltage winding and a medium-voltage winding of a target transformerThe voltage phasors of the group and low-voltage windings,
Figure FDA0002275195820000029
current phasors, I, of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer respectively1，I2，I3Are respectively the current phasors
Figure FDA00022751958200000210
The amplitude of (c).
6. The method for identifying transformer parameters based on wide-area synchronous phasor measurement system according to claim 5, wherein the synchronous phasor measurement obtained in step 2) comprises the step of obtaining the voltages of the high voltage winding, the medium voltage winding and the low voltage winding of the target transformer through synchronous phasor measurement
Figure FDA00022751958200000211
And the current of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer
Figure FDA00022751958200000212
7. The method for identifying transformer parameters based on the wide-area synchronous phasor measurement system according to claim 5, wherein the detailed steps of step 3) include:
3.1B) setting an intermediate matrix variable according to the parameter identification model:
Figure FDA00022751958200000213
Figure FDA00022751958200000214
in the above formula, K1，K2，K3Respectively the transformation ratios of a high-voltage winding, a medium-voltage winding and a low-voltage winding of a target transformer,
Figure FDA00022751958200000215
respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,
Figure FDA00022751958200000216
voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,
Figure FDA00022751958200000217
current phasors, I, of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer respectively1，I2，I3Are respectively the current phasors
Figure FDA00022751958200000218
The amplitude of (a) of (b) is,
Figure FDA00022751958200000219
representing the phasor of the current
Figure FDA00022751958200000220
And
Figure FDA00022751958200000221
the difference in the angle between the two,
Figure FDA0002275195820000031
representing the phasor of the current
Figure FDA0002275195820000032
And
Figure FDA0002275195820000033
the difference in the angle between the two,
Figure FDA0002275195820000034
representing the phasor of the current
Figure FDA0002275195820000035
And
Figure FDA0002275195820000036
the difference in the angle between the two,
Figure FDA0002275195820000037
representing the phasor of the current
Figure FDA0002275195820000038
And
Figure FDA0002275195820000039
angle difference therebetween, ξ1～ξ4Respectively represents K1，K2，
Figure FDA00022751958200000310
Error corresponding to the estimated value; according to the parameter identification model, a matrix equation shown as the following formula is established:
yk＝xkg+ξk
3.2B) forming a phasor sequence from the obtained synchronous measurement phasors
Figure FDA00022751958200000311
Wherein k is the sequence number, k is 1,2, …, n, and n is the maximum value of the sequence number; establishing a matrix equation as shown in the following formula:
Figure FDA00022751958200000312
in the above formula, y1～ynIs ykSequence of (a) x1～xnIs ykSequence of (iii) ξ1～ξnIs ξkThe sequence of (a);
3.3B) represent the matrix equation set up in step 3.2A) as y ═ Φ θ + ξ, and use a minimum of two for this equationLeast squares estimation of multiplication to obtain matrix variable theta
Figure FDA00022751958200000313
To obtain
Figure FDA00022751958200000314
Thereby obtaining the impedance of the primary winding of the target transformer
Figure FDA00022751958200000315
Impedance of low voltage winding
Figure FDA00022751958200000316
And the transformation ratio K of the primary winding of the target transformer1Transformation ratio K of medium voltage winding2The best estimate of (c); wherein y represents the sequence y1～ynThe constructed matrix, phi, representing the sequence x1～xnFormed matrix, ξ denotes the sequence ξ1～ξnA matrix is formed.
8. The utility model provides a transformer parameter identification system based on wide area synchrophasor measurement system which characterized in that includes:
the model initialization program unit is used for establishing a parameter identification model of the transformer transformation ratio and the impedance aiming at the target transformer;
the synchronous measurement phasor acquisition program unit is used for carrying out synchronous phasor measurement through the wide-area synchronous phasor measurement system to obtain the synchronous measurement phasor of the target transformer;
and the transformer parameter solving program unit is used for combining the phasor sequence with the parameter identification model and solving the transformer parameters by adopting a least square method.
9. A transformer parameter identification system based on a wide area synchronous phasor measurement system, comprising a computer device, wherein the computer device is programmed or configured to perform the steps of the transformer parameter identification method based on a wide area synchronous phasor measurement system according to any one of claims 1 to 7, or a computer program programmed or configured to perform the transformer parameter identification method based on a wide area synchronous phasor measurement system according to any one of claims 1 to 7 is stored in a memory of the computer device.
10. A computer-readable storage medium, wherein the computer-readable storage medium stores thereon a computer program programmed or configured to execute the transformer parameter identification method based on the wide-area synchrophasor measurement system according to any one of claims 1 to 7.
CN201911120033.9A 2019-11-15 2019-11-15 Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system Active CN110910001B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911120033.9A CN110910001B (en) 2019-11-15 2019-11-15 Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911120033.9A CN110910001B (en) 2019-11-15 2019-11-15 Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system

Publications (2)

Publication Number Publication Date
CN110910001A true CN110910001A (en) 2020-03-24
CN110910001B CN110910001B (en) 2022-03-22

Family

ID=69817617

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911120033.9A Active CN110910001B (en) 2019-11-15 2019-11-15 Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system

Country Status (1)

Country Link
CN (1) CN110910001B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117054798A (en) * 2023-10-12 2023-11-14 广州市扬新技术研究有限责任公司 Method and device for transformer health monitoring by utilizing parameter identification

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261297A (en) * 2008-04-17 2008-09-10 沈阳工业大学 Electric power transformer windings parameter on-line real-time identification device and method
CN103116097A (en) * 2013-01-25 2013-05-22 中国电力科学研究院 Device parameter online identification method based on multi-section hybrid measurement information
US20170023627A1 (en) * 2015-07-24 2017-01-26 Power Standards Lab Inc Method and Apparatus for Precision Phasor Measurements Through a Medium-voltage Distribution Transformer
CN106645960A (en) * 2016-11-11 2017-05-10 国家电网公司 Distribution network line mountain fire early warning method based on WAMS line parameter dynamic identification

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261297A (en) * 2008-04-17 2008-09-10 沈阳工业大学 Electric power transformer windings parameter on-line real-time identification device and method
CN103116097A (en) * 2013-01-25 2013-05-22 中国电力科学研究院 Device parameter online identification method based on multi-section hybrid measurement information
US20170023627A1 (en) * 2015-07-24 2017-01-26 Power Standards Lab Inc Method and Apparatus for Precision Phasor Measurements Through a Medium-voltage Distribution Transformer
CN106645960A (en) * 2016-11-11 2017-05-10 国家电网公司 Distribution network line mountain fire early warning method based on WAMS line parameter dynamic identification

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BOGDAN VICOL: "On-line overhead transmission line And transformer parameters identification based on PMU measurements", 《2014 INTERNATIONAL CONFERENCE AND EXPOSITION ON ELECTRICAL AND POWER ENGINEERING (EPE)》 *
牛胜锁: "面向配电网的广域测量关键技术研究及***开发", 《中国博士学位论文全文数据库-工程科技Ⅱ辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117054798A (en) * 2023-10-12 2023-11-14 广州市扬新技术研究有限责任公司 Method and device for transformer health monitoring by utilizing parameter identification
CN117054798B (en) * 2023-10-12 2023-12-22 广州市扬新技术研究有限责任公司 Method and device for transformer health monitoring by utilizing parameter identification

Also Published As

Publication number Publication date
CN110910001B (en) 2022-03-22

Similar Documents

Publication Publication Date Title
US10132853B2 (en) Wide area fault detection method using PMU data
CN103308823B (en) A kind of power distribution network single-phase disconnection phase to phase fault localization method
US7961112B2 (en) Continuous condition monitoring of transformers
CN104020370A (en) Transformer internal fault diagnosis method based on virtual parameter change monitoring
CN105653764A (en) Method for online estimating and pre-warning grid safety accident risk levels
CN104836223A (en) Power grid parameter error and bad data coordinated identification and estimation method
EP2680017A1 (en) A method of early detection of feeder lines with a high-ohm ground fault in compensated power networks
EP2715907B1 (en) Assessment of power systems
CN107248736B (en) Online identification method for positive sequence parameters of power distribution network line
CN111413589B (en) Power distribution network single-phase short circuit fault positioning method based on grey target decision
CN105447774A (en) Power grid safety accident risk grade online evaluation and early warning method and apparatus
CN104502691A (en) Online monitoring and analysis device for energy efficiency of distribution transformer
CN107797023B (en) Three-phase voltage unbalance source positioning method, device and system
CN110910001B (en) Transformer parameter identification method, system and medium based on wide-area synchronous phasor measurement system
CN108548999A (en) Cable insulation state evaluating method
KR102260550B1 (en) Facility health monitoring method by measuring the electric circuit constant inside the power facility in operation
CN104111381A (en) Dielectric loss on-line monitoring device for 35kV high voltage parallel connection power capacitor group
Ananthan et al. Model-based approach integrated with fault circuit indicators for fault location in distribution systems
CN107831378B (en) Device and method for detecting compensation effect of arc suppression coil
CN113156247A (en) Early warning method and device for low-frequency oscillation of power system
WO2020191825A1 (en) Phase-relationship-based method for single-phase disconnection fault identification in power distribution network
Lertwanitrot et al. Application of magnitude and phase angle to boundary area-based algorithm for unbalance relay protection scheme in 115-kV capacitor bank
CN105203877A (en) Single-loop power transmission line zero sequence parameter identification method capable of eliminating measurement error influences
Popov et al. Applications of synchrophasor measurement to improve the reliability and efficiency of power systems
KR20210053846A (en) Facility health monitoring method by measuring the electric circuit constant inside the power facility in operation

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