CN112327047B - 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
- CN112327047B CN112327047B CN201911231261.3A CN201911231261A CN112327047B CN 112327047 B CN112327047 B CN 112327047B CN 201911231261 A CN201911231261 A CN 201911231261A CN 112327047 B CN112327047 B CN 112327047B
- Authority
- CN
- China
- Prior art keywords
- power
- data
- active
- section
- reactive power
- 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.)
- Active
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
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 acquired 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 curve 0 And forming a through section at an interval of 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 node power data of a bus in a station, wherein the data of each node are not measured at the same time, acquiring all day data of all power nodes on the bus through a power meter and recording the corresponding time of the data, the voltage corresponding to different nodes in sequence and the current corresponding to the different nodes are
V i ,i＝(1,2,…,N)
I j ,j＝(1,2,…,N)
And 2, step: 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 change curve and 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 limits, i.e.
U imin ≤U i ≤U imax
I jmin ≤I j ≤I jmax
When the measured data is not in the interval, the data has errors, if the single point or 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 may 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:
∑P G ＝∑P L ＝∑P D +∑P S +∑P C
Q GC ＝Q G∑ +Q C∑
in the formula, sigma P G Active power, SIG P, from the power plant L The total load of the system, SIG P D -active load of user, ∑ P S -factory active load of power plant, ∑ P C Active loss of the network, Q GC Reactive power, Q, from a reactive power source G∑ Reactive power required by the load, Q C∑ -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, P j The power, P, corresponding to the nodes in turn j ' -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 S 1 ＝{O i ＝(x i ,y i ) I =1, \8230Nis the data set of the first power node on the plane. To a point O 0 Local regression line L 0 Y = ax + b, obtained by minimizing the following quadratic function:
wherein w i Is point O i The non-negative weight of (1); by weighted regression as 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 O 0 X axis and L 0 In a parallel new coordinate system;
second order regression of G * :
Obtained by minimizing the following function:
it is noted that
To G * The mapping of (c) is (0, c); finally, O * Moving to point (0, c) transformed M -1 The 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 curve 0 And forming a through section at an interval of 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 data acquisition on all the power nodes on the bus all day by a power meter and recording the corresponding time of the data, the voltage corresponding to different nodes in sequence and the current are
V i ,i＝(1,2,…,N)
I j ,j＝(1,2,…,N)
And 2, step: 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 change curve and 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.
U imin ≤U i ≤U imax
I jmin ≤I j ≤I jmax
When the measured data is not in the interval, the data has errors, if the single point or 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:
∑P G ＝∑P L ＝∑P D +∑P S +∑P C
Q GC ＝Q G∑ +Q C∑
in the formula, sigma P G Active power, SIG P, from the power plant L The total load of the system, Σ P D Subscriber's active load, SIG P S -plant utility active load, Σ P, of a power plant C Active loss of the network, Q GC Reactive power, Q, from a reactive power source GΣ Reactive power required by the load, Q CΣ 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, P j The nodes in turn corresponding to a power, P j ' -rated power corresponding to node in turn
And 4, step 4: will countThe calculated active and reactive power data are processed by fitting according to the following method, and S is set 1 ＝{O i ＝(x i ,y i ) I =1, \ 8230;, N } is the data set of the first power node on the plane. For a point O 0 Local regression line L 0 Y = ax + b, obtained by minimizing the following quadratic function:
wherein w i Is a point O i Is 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 O 0 X axis and L 0 In a parallel new coordinate system.
Second regression of G * :
Obtained by minimizing the following function:
it is noted that
To G * The mapping of (c) is (0, c). Finally, O * Moving to point (0, c) transformed M -1 The new position obtained. And performing the regression operation on each data point of the first power node, and performing the fitting processing on each remaining node according to the mode.
And 5: selecting time and section to collect data, changing the power values of all nodes into corresponding curves, and selecting section t for the curves 0 And 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 obtained by performing curve fitting by adopting the method, it can be clearly seen that, after adopting the method, in a time interval without data in fig. 1, data in the time interval can also be collected through a curve obtained by time fitting 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
V i ,i＝(1,2,…,N)
I j ,j＝(1,2,…,N)
And 2, step: 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:
determining the maximum value and the minimum value of the voltage and the current of each node according to the load variation curve and 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.
U imin ≤U i ≤U imax
I jmin ≤I j ≤I jmax
When the measured data is not in the interval, the data has errors, if the single point or 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 has 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 3, 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:
∑P G ＝∑P L ＝∑P D +∑P S +∑P C
Q GC ＝Q G∑ +Q C∑
in the formula, sigma P G Active power, sigma P, from the power plant L The total load of the system, sigma P D -active load of user, ∑ P S -factory active load of power plant, ∑ P C -active loss of network, Q GC Reactive power, Q, from a reactive power source G∑ Reactive power required by the load, Q C∑ -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, P j -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 S 1 ＝{O i ＝(x i ,y i ) I =1, \8230Nis the data set of the first power node on the plane. To a point O 0 Local regression line L 0 Y = ax + b, obtained by minimizing the following quadratic function:
wherein w i Is point O i The non-negative weight of (1); by weighted regression, P can be calculated * Is preferably a regression line L * (ii) a All data points are transformed by the transformation M to the origin O 0 X axis and L 0 In a parallel new coordinate system;
second order regression of G * :
This is obtained by minimizing the following function:
it is noted that
To G * The mapping of (a) is (0, c); finally, O * Moving to point (0, c) transformed M -1 The new position obtained; performing the regression operation on each data point of the first power node, and performing fitting processing on each remaining node according to the mode;
and 5: selecting time and section to collect data:
selecting section t for curve 0 And 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 CN112327047A (en) | 2021-02-05 |
| CN112327047B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115276228A (en) * | 2022-07-29 | 2022-11-01 | 南京国电南自电网自动化有限公司 | Method, system, storage device and computing device for synchronously uploading power data |
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 |
|---|
| 基于全节点负荷断面聚类的无功优化资产配置;桂前进等;《合肥工业大学学报(自然科学版)》;20160328;第39卷(第03期);328-332 * |
| 计及新能源电站参与的配电网全天无功计划;赵丰明等;《可再生能源》;20180716;第36卷(第07期);1067-1073 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112327047A (en) | 2021-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102128975B (en) | Voltage stabilization online monitoring phasor data measurement device and phasor measurement method | |
| CN108535543B (en) | Transformer area phase identification method based on synchronous sampling of concentrator and ammeter | |
| CN107454143B (en) | Intelligent inspection method for metering performance of digital electric energy meter in wide-area measurement mode | |
| CN102361353B (en) | Method for aligning section raw data on basis of double time scale delay evaluation | |
| CN110648249B (en) | Annual power balance measuring and calculating method, device and equipment | |
| JP2010193594A (en) | Maximum power generation amount estimating method for photovoltaic power generation system, method for controlling power distribution system, and distribution system control apparatus | |
| CN103812220A (en) | Fault hologram system based on standardization processing technology and method thereof | |
| CN102497019B (en) | Transformer substation-dispatching center two-stage state estimation multi-source data timing and splicing method | |
| CN103884931A (en) | Testing and recording device for load characteristics of transformer substation bus | |
| CN107944680A (en) | A kind of substation's electric energy balance monitoring method based on intelligent electric energy meter | |
| CN112327047B (en) | Method for realizing power same-section data measurement in transformer substation | |
| CN110736957B (en) | Online detection system and method for electric energy metering device based on clock synchronization | |
| CN102901868A (en) | Method for electric energy acquisition system data checking | |
| CN111221794A (en) | Multi-source data fusion method for wide-area measurement control system of power distribution network | |
| 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 | |
| CN103217581B (en) | The method and system of line parameter circuit value identification is realized based on stable state telemetry | |
| Bhatti et al. | Implementation of low cost non-DFT based Phasor Measurement Unit for 50 Hz power system | |
| Watitwa et al. | A review on active distribution system state estimation | |
| CN107767060B (en) | Theoretical line loss calculation system and method for distribution network line | |
| Haji et al. | Practical considerations in the design of distribution state estimation techniques | |
| Kononov et al. | The low cost method of measurements synchronization in the low-voltage network on frequency and the voltage phase | |
| Wu et al. | State Estimation of Distribution Network Considering Data Compatibility | |
| CN112886579B (en) | Power distribution network line parameter identification method considering unmatched PMU and SCADA measurement time scales | |
| 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 |
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 |