CN109581130B - Voltage curve similarity calculation method considering voltage adjustment measure influence - Google Patents

Voltage curve similarity calculation method considering voltage adjustment measure influence Download PDF

Info

Publication number
CN109581130B
CN109581130B CN201910049892.7A CN201910049892A CN109581130B CN 109581130 B CN109581130 B CN 109581130B CN 201910049892 A CN201910049892 A CN 201910049892A CN 109581130 B CN109581130 B CN 109581130B
Authority
CN
China
Prior art keywords
voltage
similarity
distribution
curves
distribution transformer
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
Application number
CN201910049892.7A
Other languages
Chinese (zh)
Other versions
CN109581130A (en
Inventor
唐泽洋
周鲲鹏
曹侃
万磊
余笑东
蔡德福
万黎
王文娜
饶渝泽
邢杰
周鸿雁
刘焱
邱丹
余飞
忻俊慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Hubei Electric Power Co Ltd Electric Power Research Institute
State Grid Corp of China SGCC
State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Hubei Electric Power Co Ltd Electric Power Research Institute
State Grid Corp of China SGCC
State Grid Hubei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Hubei Electric Power Co Ltd Electric Power Research Institute, State Grid Corp of China SGCC, State Grid Hubei Electric Power Co Ltd filed Critical State Grid Hubei Electric Power Co Ltd Electric Power Research Institute
Priority to CN201910049892.7A priority Critical patent/CN109581130B/en
Publication of CN109581130A publication Critical patent/CN109581130A/en
Application granted granted Critical
Publication of CN109581130B publication Critical patent/CN109581130B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/26Arrangements for eliminating or reducing asymmetry in polyphase networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/50Arrangements for eliminating or reducing asymmetry in polyphase networks

Abstract

The invention provides a voltage curve similarity calculation method considering the influence of voltage adjustment measures, which comprises the following steps of: A. obtaining distribution transformation related data, wherein the distribution transformation related data comprises: the corresponding relation between the 10kV line and the distribution transformer, the name of the distribution transformer, the address code of the distribution transformer and the three-phase voltage of the distribution transformer; B. b, preprocessing the distribution transformer outlet voltage data obtained in the step A; C. b, identifying a voltage catastrophe point according to the distribution transformer outlet voltage data preprocessed in the step B, and segmenting a voltage curve; D. calculating the similarity of the segmented curves according to the voltage curve segmentation result of the step C; E. and D, comprehensively judging the similarity of the voltage curves of the two transformer areas according to the similarity result of the sectional curves calculated in the step D. The method provided by the invention can be used for correctly calculating the voltage curve similarity when voltage regulation measures exist.

Description

Voltage curve similarity calculation method considering voltage adjustment measure influence
Technical Field
The invention relates to the technical field of distribution transformer operation, in particular to a voltage curve similarity calculation method considering the influence of voltage adjustment measures.
Background
When the distribution network operates, in order to reduce active network loss, balance load and improve reliability, the network structure of the distribution network needs to be adjusted, relevant information in the power consumption information acquisition system may not be updated in time after the network structure is adjusted, and generally, relevant data are verified and updated in a manual mode, which consumes time and labor.
The applicant of the present patent proposes a distribution network topology structure verification and maintenance method (publication number: CN107508297A), which realizes automatic verification and maintenance of a distribution network topology structure by analyzing distribution-transformation outlet voltage curve similarity, wherein the similarity calculation method is a correlation coefficient. However, practical operation experience shows that some transformer areas are provided with voltage adjusting equipment, when the voltage adjusting equipment is put into use, the voltage at the outlet of the distribution transformer is increased or decreased steeply, and at the moment, the similarity degree of voltage curves of the two transformer areas is directly calculated, so that the calculation result is inaccurate.
Chinese patent publication No. CN108564485A (entitled "low voltage transformer area user phase identification method based on voltage curve similarity analysis") proposes a low voltage transformer area user phase identification method based on voltage curve similarity analysis of a smart meter, and determines the similarity of voltage curves by calculating DTW distance. The chinese patent publication No. CN104092481A (entitled: a method for distinguishing between a distribution area and a phase by voltage characteristics) compares the similarity between the voltage curve of a node to be determined and the voltage curve recorded by a concentrator to determine whether the node belongs to the distribution area, wherein the similarity calculation method is a correlation coefficient method. None of the above patents relate to a voltage curve similarity calculation method that is affected by voltage adjustment measures.
Disclosure of Invention
The invention provides a voltage curve similarity calculation method considering the influence of voltage adjustment measures, and the similarity of two voltage curves of a transformer area can be correctly calculated by the method when the voltage adjustment measures exist.
The technical scheme adopted by the invention is as follows:
a voltage curve similarity calculation method considering the influence of voltage adjustment measures comprises the following steps:
A. obtaining distribution transformation related data, wherein the distribution transformation related data comprises: the corresponding relation between the 10kV line and the distribution transformer, the name of the distribution transformer, the address code of the distribution transformer and the three-phase voltage of the distribution transformer;
B. b, preprocessing the distribution transformer outlet voltage data obtained in the step A;
C. b, identifying a voltage catastrophe point according to the distribution transformer outlet voltage data preprocessed in the step B, and segmenting a voltage curve;
D. calculating the similarity of the segmented curves according to the voltage curve segmentation result of the step C;
E. and D, comprehensively judging the similarity of the voltage curves of the two transformer areas according to the similarity result of the sectional curves calculated in the step D.
Further, the preprocessing of the distribution transformer outlet voltage data in the step B specifically includes:
the production management system records the three-phase voltage value of a platform area in one day, records data of one point at regular intervals and N points all day,the A phase voltage value is Ua1,Ua2……UaNAnd the phase voltage value of B is Ub1,Ub2……UbNC phase voltage value is Uc1,Uc2……UcNCalculating the distribution transformer outlet voltage U without considering the influence of three-phase load unbalance by adopting the following formulai
Wherein U isai、Ubi、UciThe voltage of the ith point A, B, C is 1 … N, UiIn order to not consider the distribution transformer outlet voltage at the ith point when the influence of three-phase load unbalance is not considered, U can be obtained by an iteration methodi
Further, the voltage discontinuity identification method in the step C comprises:
1) calculating the voltage difference value of the two curves at the corresponding moment; distribution transformer outlet voltage of transformer area 1 is UT1iThe distribution transformer outlet voltage of the transformer area 2 is UT2iWhere i is 1 … … 96. Voltage difference value U of two station areas at corresponding time1i
U1i=UT1i-UT2i,i=1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations1iChange amount of (Δ U)1i
ΔU1i=U1(i+1)-U1i,i=1……95
3) Carrying out normal test on the variable quantity through Kolmogorov-Smirnov, and if the variable quantity is subjected to normal distribution, indicating that no mutation point exists; if the normal distribution is not obeyed, the mutation points are shown;
4) and judging abnormal points through the box body diagram, and identifying the positions of the voltage catastrophe points.
Further, the step D of calculating the similarity of the piecewise curves specifically includes:
assuming that there are two mutation points, the positions of the mutation points are respectively k1,k2Since the voltage curve has 96 points, it is based on two protrusionsThe position of the variable point can divide the voltage curve into three sections of [0, k ]1),[k1,k2),[k2,96]Respectively calculating the similarity S of the three voltage curves1,S2,S3
Further, step E comprehensively judges the similarity between the two voltage curves of the transformer area, specifically: comparing the similarity S of the piecewise curves calculated in the step D1,S2,S3And taking the maximum value as the similarity of the voltage curves of the two transformer areas.
According to the distribution transformer outlet voltage historical operation data-based distribution transformer outlet voltage distribution transformer outlet.
Drawings
FIG. 1 is a schematic flow chart of a voltage curve similarity calculation method in consideration of the influence of voltage adjustment measures according to the present invention;
FIG. 2 is a graph of typical two stage exit voltages (after pre-processing by step B) according to the present invention;
FIG. 3 is a Q-Q diagram of voltage discontinuity detection in accordance with the present invention;
fig. 4 is a box diagram of voltage discontinuity detection according to the present invention.
Detailed Description
The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings.
Fig. 1 is a schematic flow chart of a voltage curve similarity calculation method considering the influence of a voltage adjustment measure according to the present invention, which includes the following steps:
A. obtaining distribution transformation related data, wherein the distribution transformation related data comprises: the corresponding relation between the 10kV line and the distribution transformer, the name of the distribution transformer, the address code of the distribution transformer and the three-phase voltage of the distribution transformer;
B. b, preprocessing the distribution transformer outlet voltage data acquired in the step A
The production management system records three-phase voltage values of a platform area in one day, data of N points are recorded at regular intervals in total every day, and the voltage value of the phase A is Ua1,Ua2……UaNAnd the phase voltage value of B is Ub1,Ub2……UbNC phase voltage value is Uc1,Uc2……UcNCalculating the distribution transformer outlet voltage U without considering the influence of three-phase load unbalance by adopting the following formulai
Wherein U isai、Ubi、UciThe voltage of the ith point A, B, C is 1 … N, UiIn order to not consider the distribution transformer outlet voltage at the ith point when the influence of three-phase load unbalance is not considered, U can be obtained by an iteration methodi
C. B, identifying a voltage catastrophe point according to the distribution transformer outlet voltage data preprocessed in the step B, and segmenting a voltage curve;
1) calculating the voltage difference value of the two curves at the corresponding moment; distribution transformer outlet voltage of transformer area 1 is UT1iThe distribution transformer outlet voltage of the transformer area 2 is UT2iWhere i is 1 … … 96. Voltage difference value U of two station areas at corresponding time1i
U1i=UT1i-UT2i,i=1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations1iChange amount of (Δ U)1i
ΔU1i=U1(i+1)-U1i,i=1……95
3) Carrying out normal test on the variable quantity through Kolmogorov-Smirnov, and if the variable quantity is subjected to normal distribution, indicating that no mutation point exists; if the normal distribution is not obeyed, the mutation points are shown;
4) and judging abnormal points through the box body diagram, and identifying the positions of the voltage catastrophe points.
D. According to the voltage curve segmentation result of the step C, calculating the similarity of the segmentation curves
Assuming that there are two mutation points, the positions of the mutation points are respectively k1,k2Because the voltage curve has 96 points, the voltage curve can be divided into three sections of [0, k ] according to the positions of two catastrophe points1),[k1,k2),[k2,96]Respectively calculating the similarity S of the three voltage curves1,S2,S3
E. And D, comprehensively judging the similarity of the voltage curves of the two transformer areas according to the similarity result of the sectional curves calculated in the step D. Comparing the similarity S of the piecewise curves calculated in the step D1,S2,S3And taking the maximum value as the similarity of the voltage curves of the two transformer areas.
The technical scheme and effect of the invention are explained in detail by a specific embodiment as follows:
in the step A, the distribution and transformation outlet voltage data of two transformer areas (a new country transformer area and a next transformer area) supplied by a 10kV public friend pavilion line are obtained from a related system, and the distribution and transformation outlet voltage value of the next transformer area is shown as the following table:
distribution transformer outlet voltage value of channel region in 12018 years, 6 months and 1 day
In step B, the distribution transformer outlet voltage data is preprocessed, and the distribution transformer outlet voltage is calculated without considering the influence of the three-phase load imbalance as shown in fig. 2.
In the step C, 1) calculating a voltage difference value of the two curves at the corresponding moment; distribution transformer outlet voltage of transformer area 1 is UT1iThe distribution transformer outlet voltage of the transformer area 2 is UT2iWhere i is 1 … … 96. Voltage difference value U of two station areas at corresponding time1i
U1i=UT1i-UT2i,i=1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations1iChange amount of (Δ U)1i
ΔU1i=U1(i+1)-U1i,i=1……95
3) The results of the normal test on the variation by Kolmogorov-Smirnov and the normal test on the variation of the voltage difference value are shown in the following table, and it can be seen that the p value is 0<0.05, which indicates that the variation of the voltage difference value does not obey positive-too-distribution, and indicates that the outlet voltage of one of the two platform regions has a mutation value. The results of the Q-Q plot shown in fig. 3 also show non-normal distribution compliance.
TABLE 2 results of the normality test
4) And judging abnormal points through the box body diagram, and identifying the positions of the voltage catastrophe points. It is found from the box diagram shown in fig. 4 that there are two abrupt change points, which are the 19 th and 78 th data points (corresponding to the 20 th and 79 th points of the original data), and the voltage curve is divided into three segments according to the positions of the abrupt change points, wherein 1-19 are the first segment, 20-78 are the second segment, and 79-96 are the third segment.
In step D, the positions of the two mutation points are respectively k1=20,k2When the voltage curve is 79, the voltage curve can be divided into three segments [0,20 ], [20,79 ], [79,96 ] according to the positions of the two discontinuities]Respectively calculating the similarity S of the three voltage curves1,S2,S3Wherein the magnitude of the similarity is represented by calculating a correlation coefficient.
TABLE 3 correlation coefficient segmentation calculation results
Calculation results Correlation coefficient
Integral body 0.077
First stage 0.76
Second section 0.90
Third stage 0.93
In the step E, according to the calculation result of the step D, the similarity of the third segment is 0.93 at most in the segment similarity, so that the similarity of the two curves is 0.93.
If the similarity of the voltage outlet curves of the two transformer areas is not calculated by adopting the method provided by the invention, the similarity of the voltage outlet curves of the two transformer areas obtained by direct calculation is 0.077, the calculation result is not accurate, and the calculation result is accurate by adopting the method provided by the invention.
The above description is only an embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention.

Claims (4)

1. A voltage curve similarity calculation method considering the influence of voltage adjustment measures is characterized by comprising the following steps:
A. obtaining distribution transformation related data, wherein the distribution transformation related data comprises: the corresponding relation between the 10kV line and the distribution transformer, the name of the distribution transformer, the address code of the distribution transformer and the three-phase voltage of the distribution transformer;
B. b, preprocessing the distribution transformer outlet voltage data obtained in the step A;
C. b, identifying a voltage catastrophe point according to the distribution transformer outlet voltage data preprocessed in the step B, and segmenting a voltage curve;
D. calculating the similarity of the segmented curves according to the voltage curve segmentation result of the step C;
E. d, comprehensively judging the similarity of the voltage curves of the two transformer areas according to the similarity result of the sectional curves calculated in the step D;
the step C of identifying the voltage catastrophe point specifically comprises the following steps:
1) calculating the voltage difference value of the two curves at the corresponding moment; distribution transformer outlet voltage of transformer area 1 is UT1iThe distribution transformer outlet voltage of the transformer area 2 is UT2iWhere i is 1 … … 96, the voltage difference U between two stations at the corresponding time1i
U1i=UT1i-UT2i,i=1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations1iChange amount of (Δ U)1i
ΔU1i=U1(i+1)-U1i,i=1……95
3) Carrying out normal test on the variable quantity through Kolmogorov-Smirnov, and if the variable quantity is subjected to normal distribution, indicating that no mutation point exists; if the normal distribution is not obeyed, the mutation points are shown;
4) and judging abnormal points through the box body diagram, and identifying the positions of the voltage catastrophe points.
2. The voltage curve similarity calculation method considering the influence of the voltage adjustment measure according to claim 1, wherein: the step B of preprocessing the distribution transformer outlet voltage data obtained in the step A specifically comprises the following steps:
the production management system records three-phase voltage values of a platform area in one day, data of N points are recorded at regular intervals in total every day, and the voltage value of the phase A is Ua1,Ua2……UaNAnd the phase voltage value of B is Ub1,Ub2……UbNC phase voltage value is Uc1,Uc2……UcNThe shadow without considering the three-phase load unbalance is calculated by the following formulaTime-of-response distribution transformer outlet voltage Ui
Wherein U isai、Ubi、UciThe voltage of the ith point A, B, C is 1 … N, UiIn order to not consider the distribution transformer outlet voltage at the ith point when the influence of three-phase load unbalance is not considered, U can be obtained by an iteration methodi
3. The voltage curve similarity calculation method considering the influence of the voltage adjustment measure according to claim 1, wherein: step D, calculating the similarity of the piecewise curves specifically includes:
assuming that there are two mutation points, the positions of the mutation points are respectively k1,k2Because the voltage curve has 96 points, the voltage curve is divided into three sections of [0, k ] according to the positions of two catastrophe points1),[k1,k2),[k2,96]Respectively calculating the similarity S of the three voltage curves1,S2,S3
4. A voltage curve similarity calculation method considering the influence of voltage adjustment measures according to claim 3, characterized in that: step E, comprehensively judging the similarity of the voltage curves of the two transformer areas, which specifically comprises the following steps: comparing the similarity S of the piecewise curves calculated in the step D1,S2,S3And taking the maximum value as the similarity of the voltage curves of the two transformer areas.
CN201910049892.7A 2019-01-18 2019-01-18 Voltage curve similarity calculation method considering voltage adjustment measure influence Active CN109581130B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910049892.7A CN109581130B (en) 2019-01-18 2019-01-18 Voltage curve similarity calculation method considering voltage adjustment measure influence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910049892.7A CN109581130B (en) 2019-01-18 2019-01-18 Voltage curve similarity calculation method considering voltage adjustment measure influence

Publications (2)

Publication Number Publication Date
CN109581130A CN109581130A (en) 2019-04-05
CN109581130B true CN109581130B (en) 2020-10-16

Family

ID=65917014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910049892.7A Active CN109581130B (en) 2019-01-18 2019-01-18 Voltage curve similarity calculation method considering voltage adjustment measure influence

Country Status (1)

Country Link
CN (1) CN109581130B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112564110B (en) * 2021-02-24 2021-07-02 国网江西省电力有限公司电力科学研究院 Transformer area low-voltage treatment method and system

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104092481B (en) * 2014-07-17 2016-04-20 江苏林洋能源股份有限公司 A kind of by voltage characteristic differentiation Tai Qu and phase method for distinguishing
JP6392155B2 (en) * 2015-03-27 2018-09-19 株式会社東芝 Power distribution monitoring and control device
CN105184070A (en) * 2015-08-31 2015-12-23 华南理工大学 Voltage integral method based segmentation fitting method for calculating volt-time characteristic curve
CN107508297B (en) * 2017-08-10 2019-07-26 国家电网公司 A kind of verification of distribution topological structure and maintaining method
CN108199375B (en) * 2018-01-31 2021-01-22 天津大学 Intelligent power distribution network topology identification method based on synchronous phasor measurement
CN108564485A (en) * 2018-04-16 2018-09-21 国网河南省电力公司电力科学研究院 Low-voltage platform area user's phase recognition methods based on voltage curve similarity analysis
CN108600362A (en) * 2018-04-19 2018-09-28 南京林洋电力科技有限公司 A method of managing taiwan area archives automatically based on voltage characteristic

Also Published As

Publication number Publication date
CN109581130A (en) 2019-04-05

Similar Documents

Publication Publication Date Title
Wang et al. Multi-criteria building energy performance benchmarking through variable clustering based compromise TOPSIS with objective entropy weighting
US20150316620A1 (en) Method for identifying a system anomaly in a power distribution system
CN109581130B (en) Voltage curve similarity calculation method considering voltage adjustment measure influence
CN109596929B (en) Voltage curve similarity judgment method considering clock asynchronous influence
Bin Othman et al. Determination of transmission reliability margin using parametric bootstrap technique
WO2020098728A1 (en) Cluster load prediction method and apparatus, and storage medium
CN107508297A (en) A kind of distribution topological structure verification and maintaining method
Mutanen et al. Customer classification and load profiling based on AMR measurements
CN107203947A (en) A kind of digital substation relay protection evaluation mode of system condition method
CN106682763B (en) Power load optimization prediction method for large amount of sample data
Chen et al. Change detection of electric customer behavior based on AMR measurements
Ding et al. Individual nonparametric load estimation model for power distribution network planning
Grigoraş et al. Load estimation for distribution systems using clustering techniques
CN108875972A (en) A kind of power grid O&amp;M expense prediction technique and forecasting system
CN109598435B (en) Power distribution network cable running state evaluation method and system
CN111539657A (en) Typical electricity consumption industry load characteristic classification and synthesis method combined with user daily electricity consumption curve
CN110991527A (en) Similarity threshold determination method considering voltage curve average fluctuation rate
CN111061821A (en) Low-voltage distribution network topology verification method and system based on improved k value clustering algorithm
Canale et al. A comparative analysis among Standard Load Profiles for Natural Gas consumption simulation at urban scale
CN109784777B (en) Power grid equipment state evaluation method based on time sequence information fragment cloud similarity measurement
CN110264107B (en) Large data technology-based abnormal diagnosis method for line loss rate of transformer area
CN111091223A (en) Distribution transformer short-term load prediction method based on Internet of things intelligent sensing technology
Sarjan et al. Loss estimation in power distribution networks with limited information using clustering techniques
CN108665090B (en) Urban power grid saturation load prediction method based on principal component analysis and Verhulst model
CN112712203A (en) Method and system for predicting daily maximum load of 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