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 PDFInfo
- 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
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 86
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000009466 transformation Effects 0.000 claims abstract description 41
- 238000004804 winding Methods 0.000 claims description 137
- 239000011159 matrix material Substances 0.000 claims description 48
- 238000004590 computer program Methods 0.000 claims description 15
- 238000012544 monitoring process Methods 0.000 abstract description 15
- 238000007405 data analysis Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 3
- 238000000819 phase cycle Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009666 routine test Methods 0.000 description 1
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
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:
in the above formula, K is the transformer transformation ratio of the target transformer,impedance of the target transformerThe sum of impedance parameters after the primary winding and the secondary winding of the target transformer are converted to the primary side,andthe currents of the primary side and the secondary side of the target transformer are respectively,andthe 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 measurementAndand the primary and secondary voltages of the target transformerAnd
optionally, the detailed steps of step 3) include:
3.1A) setting an intermediate matrix variable according to the parameter identification model:
in the above formula, K is the transformer transformation ratio of the target transformer,is the impedance of the target transformer and,andthe currents of the primary side and the secondary side of the target transformer are respectively,andvoltage 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 phasorsWherein 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:
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 methodTo obtainThereby obtaining the transformation ratio K and the impedance of the transformerThe 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, impedanceAs 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:
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,respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,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 phasorsThe 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 measurementAnd the current of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer
Optionally, the detailed steps of step 3) include:
3.1B) setting an intermediate matrix variable according to the parameter identification model:
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,respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,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 phasorsThe amplitude of (a) of (b) is,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndangle difference therebetween, ξ1~ξ4Respectively represents K1,K2,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 phasorsWherein 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:
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 equationTo obtainThereby obtaining the impedance of the primary winding of the target transformerImpedance of low voltage windingAnd 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, whereinThe voltage and the current of the primary side and the secondary side of the two-winding transformer are respectively,respectively converted into the voltage and current values of the primary side, Gm+jBmFor the admittance parameter of the excitation branch,in order to be the exciting current,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:
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:
in the above formula, K is the transformer transformation ratio of the target transformer,impedance of the target transformerThe sum of impedance parameters after the primary winding and the secondary winding of the target transformer are converted to the primary side,andthe currents of the primary side and the secondary side of the target transformer are respectively,andthe 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 transformerAndand the primary and secondary voltages of the target transformerAnd
in this embodiment, the detailed steps of step 3) include:
3.1A) setting an intermediate matrix variable according to the parameter identification model:
in the above formula, K is the transformer transformation ratio of the target transformer,is the impedance of the target transformer and,andthe currents of the primary side and the secondary side of the target transformer are respectively,andvoltage 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 phasorsWherein 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:
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 methodTo obtainThereby obtaining the transformation ratio K and the impedance of the transformerThe 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, impedanceKnown 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:
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,respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,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 phasorsThe 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 transformerAnd the current phasor of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer
In this embodiment, the detailed steps of step 3) include:
3.1B) setting an intermediate matrix variable according to the parameter identification model:
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,respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,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 phasorsThe amplitude of (a) of (b) is,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndangle difference therebetween, ξ1~ξ4Respectively represents K1,K2,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 phasorsWherein 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:
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 equationTo obtainThereby obtaining the impedance of the primary winding of the target transformerImpedance of low voltage windingAnd 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:
in the above formula, K is the transformer transformation ratio of the target transformer,impedance of the target transformerThe sum of impedance parameters after the primary winding and the secondary winding of the target transformer are converted to the primary side,andthe currents of the primary side and the secondary side of the target transformer are respectively,andthe voltages of the primary side and the secondary side of the target transformer are respectively.
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:
in the above formula, K is the transformer transformation ratio of the target transformer,is the impedance of the target transformer and,andthe currents of the primary side and the secondary side of the target transformer are respectively,andvoltage 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 phasorsWherein 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:
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 methodTo obtainThereby obtaining the transformation ratio K and the impedance of the transformerThe 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, impedanceAs 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:
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,respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,respectively a high-voltage winding and a medium-voltage winding of a target transformerThe voltage phasors of the group and low-voltage windings,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 phasorsThe 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 measurementAnd the current of the high-voltage winding, the medium-voltage winding and the low-voltage winding of the target transformer
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:
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,respectively impedance phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer,voltage phasors of a high-voltage winding, a medium-voltage winding and a low-voltage winding of the target transformer are respectively,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 phasorsThe amplitude of (a) of (b) is,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndthe difference in the angle between the two,representing the phasor of the currentAndangle difference therebetween, ξ1~ξ4Respectively represents K1,K2,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 phasorsWherein 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:
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 thetaTo obtainThereby obtaining the impedance of the primary winding of the target transformerImpedance of low voltage windingAnd 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.
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)
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)
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 |
-
2019
- 2019-11-15 CN CN201911120033.9A patent/CN110910001B/en active Active
Patent Citations (4)
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)
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)
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 |