CN112327047A - Method for realizing power same-section data measurement in transformer substation - Google Patents
Method for realizing power same-section data measurement in transformer substation Download PDFInfo
- Publication number
- CN112327047A CN112327047A CN201911231261.3A CN201911231261A CN112327047A CN 112327047 A CN112327047 A CN 112327047A CN 201911231261 A CN201911231261 A CN 201911231261A CN 112327047 A CN112327047 A CN 112327047A
- Authority
- CN
- China
- Prior art keywords
- power
- data
- active
- section
- node
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/06—Arrangements for measuring electric power or power factor by measuring current and voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/001—Measuring real or reactive component; Measuring apparent energy
- G01R21/002—Measuring real component
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/001—Measuring real or reactive component; Measuring apparent energy
- G01R21/003—Measuring reactive component
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
A method for realizing power co-section data measurement in a transformer substation is characterized in that all power nodes are subjected to simultaneous interval time point data acquisition, and a data fitting technology is applied to construct a numerical curve of all the power nodes, so that the time co-section measurement of the power data is realized. The method comprises the following steps: acquiring power data of a power node of an in-station bus, determining a power change range, determining a maximum value and a minimum value of power according to historical normal data and a load change curve, processing error data, and calculating active reactive power; further determining the error of the collected data according to the active and reactive power balance; fitting the calculated active power and reactive power; selecting time and section to collect data, and selecting section t for curve0And forming a through section at an interval delta t, and collecting data.
Description
Technical Field
The invention relates to the fields of power system operation, power system stability analysis and data acquisition and processing, in particular to a method for realizing power same-section data measurement in a transformer substation.
Background
The power system is a time-dependent system, and is based on the waveform of a time axis no matter the voltage, the current, the phase angle and the power angle change. In recent years, supercritical and ultra-supercritical units are successively connected to the grid for operation, large-area power grids are interconnected, and an ultra-high voltage transmission technology is developed. The safe and stable operation of the power grid provides new requirements for power automation equipment, particularly for time synchronization and the requirement of operation of a relay protection device, an automation device, a safety and stability control system, an Energy Management System (EMS), a production information management system and the like based on a unified time reference so as to meet the requirements of event sequence recording (SOE), fault recording and real-time data acquisition time consistency, ensure the accuracy of line fault location, phasor and power angle dynamic monitoring, unit and power grid parameter verification and the level of power grid accident analysis and stability control, and improve the operation efficiency and reliability. The popularization and application of the digital power technology in the future have higher requirements on time synchronization.
At present, only the clock synchronization between stations and between each device in the stations is technically researched, and the conventional time synchronization technology mainly comprises the following 5 schemes of realizing the satellite time synchronization scheme (1) and receiving the standard time output by a Global Positioning System (GPS) or a Beidou system as the system time; (2) the wired transmission time synchronization scheme adopts the time coding IRIGB or utilizes serial port transmission and the like to realize the precise synchronization of time in short distance; (3) network Time Protocol (NTP) or Simple Network Time Protocol (SNTP), where NTP or SNTP is a scheme for performing high-precision time synchronization using ethernet packets; (4) a precision clock time synchronization (PTP) protocol is a scheme for improving the time transmission precision by adding time information in a physical layer protocol of a network; (5) the time synchronization scheme based on SDH uses SDH as a transmission medium and adopts a bidirectional time synchronization protocol to carry out time calibration so as to realize accurate time synchronization. The method cannot ensure that the collection of the numerical values of all the power nodes on the bus is realized at the same time point during the power measurement in the transformer substation.
Disclosure of Invention
The invention provides a method for realizing power co-section data measurement in a transformer substation, which is characterized in that time interval time point data collection is carried out on all power nodes at the same time, and a numerical curve of all the power nodes is constructed by using a data fitting technology, so that time co-section measurement is carried out on the power data.
The technical scheme of the invention is as follows:
a method for realizing power same-section data measurement in a transformer substation comprises the following steps:
step 1: acquiring power data of power nodes of buses in a station, wherein the data of each node are not measured at the same time, acquiring all day data of all the power nodes on the buses through a power meter and recording the corresponding time of the data, the voltage corresponding to different nodes in sequence and the current are
Vi,i=(1,2,…,N)
Ij,j=(1,2,…,N)
Step 2: determining a power variation range;
step 2.1: determining the maximum value and the minimum value of power according to historical normal data and a load change curve:
and determining the maximum value and the minimum value of the voltage and the current of each node according to the load variation curve and the historical data of each power node under the normal condition.
Step 2.2: and (3) processing error data:
the line voltage and current magnitudes must not exceed certain ranges, i.e.
Uimin≤Ui≤Uimax
Ijmin≤Ij≤Ijmax
When the measured data is not in the interval, the data has errors, if the single point or the two points exceed the normal range, the fault of the measuring equipment is presumed, and the error data is deleted; if the multiple points exceed the normal range, supposing that the line is possible to have an accident, and starting wave recording;
step 2.3: calculating active and reactive power:
and calculating the active and reactive power of each point according to the voltage and current collected by each point:
in the formula, P is active power, Q is reactive power, U is voltage, I is current,-the phase difference of the voltage and the current;
and step 3: and further determining the error of the acquired data according to the active and reactive power balance, wherein the active power balance and the reactive power balance are as follows:
∑PG=∑PL=∑PD+∑PS+∑PC
QGC=QG∑+QC∑
in the formula, sigma PGActive power, SIG P, from the power plantLThe total load of the system, SIG PDSubscriber's active load, SIG PS-factory active load of power plant, ∑ PCActive loss of the network, QGCReactive power, Q, from a reactive power sourceG∑Reactive power required by the load, QC∑-reactive losses in the network;
satisfying the following formula indicates that the power balance is such that the calculated data is error free:
in the formula, PjThe power, P, corresponding to the nodes in turnj' -the nodes correspond to the rated power in turn;
and 4, step 4: fitting the calculated active power and reactive power according to the following method, and setting S1={Oi=(xi,yi) I 1, …, N is the data set of the first power node on the plane. For a point O0Local regression line L0Y ═ ax + b, obtained by minimizing the following quadratic function:
wherein wiIs point OiThe non-negative weight of (1); tong (Chinese character of 'tong')By the weighted regression mentioned above, P can be calculated*Is preferably a regression line L*(ii) a All data points are transformed by the transformation M to the origin O0X axis and L0In a parallel new coordinate system;second regression of G*:Obtained by minimizing the following function:
it is noted thatTo G*The mapping of (a) is (0, c); finally, O*Moving to point (0, c) transformed M-1The new position obtained; the regression operation is performed on each data point of the first power node, and each remaining node is fitted in the manner described above.
And 5: selecting time and section to collect data:
selecting section t for curve0And forming a through section at an interval delta t, and collecting data.
By utilizing a curve fitting technology and power balance of a power grid, the fault of the bus can be judged according to data, and the error range of the data can be estimated. The method comprises the steps of firstly, recording all-day data of each power node and time corresponding to the data, determining a power change range for each node according to historical data under normal conditions, setting the direction of inflow and outflow of the power nodes, judging data errors according to power balance of a power system and rated power of the nodes, changing the data of each node into a curve on a two-dimensional plane by curve fitting, and setting time intervals to realize same-section measurement on the power data of each node visually. The invention has the beneficial effects that: the active and reactive power is calculated by using the collected data, the error range is determined, the same data at different moments are determined by curve fitting of the optimal mathematical curve, the same-section data collection of the data is realized, the generated same-section data is visual, and the error is small.
Drawings
Fig. 1 is a graph of power load variation in an electric power system.
Fig. 2 is a flow chart of error data processing.
Detailed Description
The following further describes the embodiments of the present invention with reference to the drawings.
The method for realizing the power same-section data measurement in the transformer substation comprises the following steps:
step 1: collecting voltage and current data of power nodes of a bus in a station, wherein the data of each node are not measured at the same time, performing all-day data collection on all the power nodes on the bus through a power meter and recording the corresponding time of the data, the voltage corresponding to different nodes is sequentially obtained, and the current is
Vi,i=(1,2,…,N)
Ij,j=(1,2,…,N)
Step 2: determining a power variation range:
step 2.1: determining the maximum value and the minimum value according to historical normal data and a load change curve:
and determining the maximum value and the minimum value of the voltage and the current of each node according to the load variation curve and the historical data of each power node under the normal condition.
Step 2.2: and (3) processing error data:
the line voltage and current magnitudes must not exceed certain ranges, i.e.
Uimin≤Ui≤Uimax
Ijmin≤Ij≤Ijmax
When the measured data is not in the interval, the data has errors, if the single point or the two points exceed the normal range, the fault of the measuring equipment is presumed, and the error data is deleted; if the multiple points exceed the normal range, the accident of the line is presumed to occur, and the recording record is started.
Step 2.3: calculating active and reactive power:
calculating active and reactive power according to the collected voltage and current:
in the formula, P is active power, Q is reactive power, U is voltage, I is current,-the phase difference of the voltage and the current.
And step 3: and further determining the error of the calculated data according to the active and reactive power balance, wherein the active power balance and the reactive power balance are as follows:
∑PG=∑PL=∑PD+∑PS+∑PC
QGC=QG∑+QC∑
in the formula, sigma PGActive power, SIG P, from the power plantLTotal load of the system, Σ PDSubscriber's active load, SIG PS-plant utility active load, Σ P, of a power plantCActive loss of the network, QGCReactive power, Q, from a reactive power sourceGΣReactive power required by the load, QCΣReactive losses in the network
Satisfying the following formula indicates that the power balance, i.e. the measured values are not erroneous:
in the formula, PjThe power, P, corresponding to the nodes in turnj' -rated power corresponding to node in turn
And 4, step 4: fitting the calculated active and reactive power data according to the following method, and setting S1={Oi=(xi,yi) I 1, …, N is the data set of the first power node on the plane. For a point O0Local regression line L0Y ═ ax + b, obtained by minimizing the following quadratic function:
wherein wiIs point OiIs not a negative weight. By the weighted regression described above, P can be calculated*Is preferably a regression line L*. All data points are transformed by the transformation M to the origin O0X axis and L0In a parallel new coordinate system.Second regression of G*:Obtained by minimizing the following function:
it is noted thatTo G*The mapping of (c) is (0, c). Finally, O*Moving to point (0, c) transformed M-1The new position obtained. The regression operation is performed on each data point of the first power node, and each remaining node is fitted in the manner described above.
And 5: selecting time and section to collect data, the power values of all nodes will become corresponding curves,selecting section t for curve0And forming a through section at an interval delta t, and collecting data.
The method of the invention is described by taking data of any two points on a 110KV substation bus within one hour as an example:
and a node A:
and the node B:
taking active power as an example, fig. 1 is an image of two nodes without adopting the method, and fig. 2 is an image after curve fitting is performed by adopting the method, it can be clearly seen that, after adopting the method, in a time interval without data in fig. 1, data of the time interval can also be collected through a curve fitted by time in fig. 2, so as to form data acquisition on the same section.
The above description is only exemplary of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (1)
1. A method for realizing power same-section data measurement in a transformer substation is characterized by comprising the following steps:
step 1: acquiring power data of power nodes of buses in a station, wherein the data of each node are not measured at the same time, acquiring all day data of all the power nodes on the buses through a power meter and recording the corresponding time of the data, the voltage corresponding to different nodes in sequence and the current are
Vi,i=(1,2,…,N)
Ij,j=(1,2,…,N)
Step 2: determining a power variation range;
step 2.1: determining the maximum value and the minimum value of power according to historical normal data and a load change curve:
and determining the maximum value and the minimum value of the voltage and the current of each node according to the load variation curve and the historical data of each power node under the normal condition.
Step 2.2: and (3) processing error data:
the line voltage and current magnitudes must not exceed certain ranges, i.e.
Uimin≤Ui≤Uimax
Ijmin≤Ij≤Ijmax
When the measured data is not in the interval, the data has errors, if the single point or the two points exceed the normal range, the fault of the measuring equipment is presumed, and the error data is deleted; if the multiple points exceed the normal range, supposing that the line is possible to have an accident, and starting wave recording;
step 2.3: calculating active and reactive power:
and calculating the active and reactive power of each point according to the voltage and current collected by each point:
in the formula, P is active power, Q is reactive power, U is voltage, I is current,-the phase difference of the voltage and the current;
and step 3: and further determining the error of the acquired data according to the active and reactive power balance, wherein the active power balance and the reactive power balance are as follows:
∑PG=∑PL=∑PD+∑PS+∑PC
QGC=QG∑+QC∑
in the formula, sigma PGActive power, SIG P, from the power plantLThe total load of the system, SIG PDSubscriber's active load, SIG PS-factory active load of power plant, ∑ PCActive loss of the network, QGCReactive power, Q, from a reactive power sourceG∑Reactive power required by the load, QC∑-reactive losses in the network;
satisfying the following formula indicates that the power balance is such that the calculated data is error free:
in the formula, Pj-power, P ', which nodes correspond in turn'j-the nodes in turn correspond to a nominal power;
and 4, step 4: fitting the calculated active power and reactive power according to the following method, and setting S1={Oi=(xi,yi) I 1, …, N is the data set of the first power node on the plane. For a point O0Local regression line L0Y ═ ax + b, obtained by minimizing the following quadratic function:
wherein wiIs point OiThe non-negative weight of (1); by the weighted regression described above, P can be calculated*Is preferably a regression line L*(ii) a All data points are transformed by the transformation M to the origin O0X axis and L0In a parallel new coordinate system;second regression ofObtained by minimizing the following function:
it is noted thatTo G*The mapping of (a) is (0, c); finally, O*Moving to point (0, c) transformed M-1The new position obtained; the regression operation is performed on each data point of the first power node, and each remaining node is fitted in the manner described above.
And 5: selecting time and section to collect data:
selecting section t for curve0And forming a through section at an interval delta t, and collecting data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911231261.3A CN112327047B (en) | 2019-12-05 | 2019-12-05 | Method for realizing power same-section data measurement in transformer substation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911231261.3A CN112327047B (en) | 2019-12-05 | 2019-12-05 | Method for realizing power same-section data measurement in transformer substation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112327047A true CN112327047A (en) | 2021-02-05 |
CN112327047B CN112327047B (en) | 2022-11-15 |
Family
ID=74319799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911231261.3A Active CN112327047B (en) | 2019-12-05 | 2019-12-05 | Method for realizing power same-section data measurement in transformer substation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112327047B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024021559A1 (en) * | 2022-07-29 | 2024-02-01 | 南京国电南自电网自动化有限公司 | Power data synchronous uploading method and system, storage and computing device |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101291061A (en) * | 2008-05-16 | 2008-10-22 | 南京南瑞继保电气有限公司 | Status estimating method for dynamic process of electrical power system |
JP2010002386A (en) * | 2008-06-23 | 2010-01-07 | Central Res Inst Of Electric Power Ind | Fault locator, fault localization method, and fault localization program |
CN101924393A (en) * | 2010-08-11 | 2010-12-22 | 徐文远 | Method and system for estimating load model parameters of electric system on line |
CN102611106A (en) * | 2012-03-31 | 2012-07-25 | 凯里供电局 | Maximum load supply capability evaluation method of medium-voltage power distribution network for loop power supply |
CN103606921A (en) * | 2013-12-10 | 2014-02-26 | 上海交通大学 | Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence |
CN104463357A (en) * | 2014-11-27 | 2015-03-25 | 国家电网公司 | Method for evaluating random intermittent DG optimized integration based on random optimal power flow |
CN104794206A (en) * | 2015-04-23 | 2015-07-22 | 国网山东省电力公司 | Transformer substation data quality evaluation system and method |
WO2015176564A1 (en) * | 2014-05-22 | 2015-11-26 | 袁志贤 | Method of transformer internal fault diagnosis based on monitoring on virtual parameter changes |
CN108110754A (en) * | 2017-12-29 | 2018-06-01 | 山东大学 | A kind of strategy for security correction method for considering to coordinate between poly zoarium |
CN108304976A (en) * | 2018-03-06 | 2018-07-20 | 西安交通大学 | A kind of electric system fining load prediction and analysis method |
CN109119993A (en) * | 2018-07-27 | 2019-01-01 | 国网辽宁省电力有限公司电力科学研究院 | One kind being based on nine area's optimal control policy analysis method of distribution network system various dimensions |
CN109416379A (en) * | 2016-04-22 | 2019-03-01 | 戴普赛斯股份公司 | The method for determining the mutual voltage sensitivity coefficient between the multiple measuring nodes of power network |
CN109728646A (en) * | 2018-12-29 | 2019-05-07 | 国网北京市电力公司 | A kind of distribution system and method having synchronous phasor measurement function |
CN109752629A (en) * | 2017-11-07 | 2019-05-14 | 中国电力科学研究院有限公司 | A kind of power grid measurement problem intelligent diagnosing method and system |
-
2019
- 2019-12-05 CN CN201911231261.3A patent/CN112327047B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101291061A (en) * | 2008-05-16 | 2008-10-22 | 南京南瑞继保电气有限公司 | Status estimating method for dynamic process of electrical power system |
JP2010002386A (en) * | 2008-06-23 | 2010-01-07 | Central Res Inst Of Electric Power Ind | Fault locator, fault localization method, and fault localization program |
CN101924393A (en) * | 2010-08-11 | 2010-12-22 | 徐文远 | Method and system for estimating load model parameters of electric system on line |
CN102611106A (en) * | 2012-03-31 | 2012-07-25 | 凯里供电局 | Maximum load supply capability evaluation method of medium-voltage power distribution network for loop power supply |
CN103606921A (en) * | 2013-12-10 | 2014-02-26 | 上海交通大学 | Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence |
WO2015176564A1 (en) * | 2014-05-22 | 2015-11-26 | 袁志贤 | Method of transformer internal fault diagnosis based on monitoring on virtual parameter changes |
CN104463357A (en) * | 2014-11-27 | 2015-03-25 | 国家电网公司 | Method for evaluating random intermittent DG optimized integration based on random optimal power flow |
CN104794206A (en) * | 2015-04-23 | 2015-07-22 | 国网山东省电力公司 | Transformer substation data quality evaluation system and method |
CN109416379A (en) * | 2016-04-22 | 2019-03-01 | 戴普赛斯股份公司 | The method for determining the mutual voltage sensitivity coefficient between the multiple measuring nodes of power network |
CN109752629A (en) * | 2017-11-07 | 2019-05-14 | 中国电力科学研究院有限公司 | A kind of power grid measurement problem intelligent diagnosing method and system |
CN108110754A (en) * | 2017-12-29 | 2018-06-01 | 山东大学 | A kind of strategy for security correction method for considering to coordinate between poly zoarium |
CN108304976A (en) * | 2018-03-06 | 2018-07-20 | 西安交通大学 | A kind of electric system fining load prediction and analysis method |
CN109119993A (en) * | 2018-07-27 | 2019-01-01 | 国网辽宁省电力有限公司电力科学研究院 | One kind being based on nine area's optimal control policy analysis method of distribution network system various dimensions |
CN109728646A (en) * | 2018-12-29 | 2019-05-07 | 国网北京市电力公司 | A kind of distribution system and method having synchronous phasor measurement function |
Non-Patent Citations (2)
Title |
---|
桂前进等: "基于全节点负荷断面聚类的无功优化资产配置", 《合肥工业大学学报(自然科学版)》 * |
赵丰明等: "计及新能源电站参与的配电网全天无功计划", 《可再生能源》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024021559A1 (en) * | 2022-07-29 | 2024-02-01 | 南京国电南自电网自动化有限公司 | Power data synchronous uploading method and system, storage and computing device |
Also Published As
Publication number | Publication date |
---|---|
CN112327047B (en) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108535543B (en) | Transformer area phase identification method based on synchronous sampling of concentrator and ammeter | |
CN102128975B (en) | Voltage stabilization online monitoring phasor data measurement device and phasor measurement method | |
RU2577245C2 (en) | Merging unit and method of operating merging unit | |
CN107454143B (en) | Intelligent inspection method for metering performance of digital electric energy meter in wide-area measurement mode | |
CN106199183A (en) | A kind of PMU realizing sub-synchronous oscillation on-line identification alarm and method | |
CN102497019B (en) | Transformer substation-dispatching center two-stage state estimation multi-source data timing and splicing method | |
CN107944680A (en) | A kind of substation's electric energy balance monitoring method based on intelligent electric energy meter | |
CN102944860A (en) | Method and device for comparing running errors of optical current transformer and electromagnetic current transformer in real time | |
CN103884931A (en) | Testing and recording device for load characteristics of transformer substation bus | |
CN110736957B (en) | Online detection system and method for electric energy metering device based on clock synchronization | |
CN112327047B (en) | Method for realizing power same-section data measurement in transformer substation | |
CN102901868A (en) | Method for electric energy acquisition system data checking | |
Ahmad Khan et al. | PTP‐based time synchronisation of smart meter data for state estimation in power distribution networks | |
CN111221794A (en) | Multi-source data fusion method for wide-area measurement control system of power distribution network | |
Watitwa et al. | A review on active distribution system state estimation | |
CN102902894A (en) | Method for evaluating the data quality and estimating the angle error of PMU (Phasor Measurement Unit) of control center based on difference comparison | |
Bhatti et al. | Implementation of low cost non-DFT based Phasor Measurement Unit for 50 Hz power system | |
CN103217581B (en) | The method and system of line parameter circuit value identification is realized based on stable state telemetry | |
Kononov et al. | The low cost method of measurements synchronization in the low-voltage network on frequency and the voltage phase | |
CN103236916A (en) | Dynamic compensation method for access network delay in SV (Sampled Value) networking of digitalized relay protection device | |
CN110555244B (en) | Synchronous line loss calculation system and method for power distribution network | |
CN114172262A (en) | Intelligent substation sampling data quality comprehensive evaluation method and system | |
Khan et al. | Three phase state estimation in power distribution networks by integrating IEEE-1588 with smart meters | |
CN102411085A (en) | Alternating current collecting and merging device and alternating current collecting method | |
CN111709612A (en) | Power distribution network state estimation method considering collected historical data |
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 |