CN109581130B  Voltage curve similarity calculation method considering voltage adjustment measure influence  Google Patents
Voltage curve similarity calculation method considering voltage adjustment measure influence Download PDFInfo
 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
Links
 238000004364 calculation method Methods 0.000 title claims abstract description 22
 230000000875 corresponding Effects 0.000 claims abstract description 13
 230000001131 transforming Effects 0.000 claims abstract description 10
 238000007781 preprocessing Methods 0.000 claims abstract description 7
 230000011218 segmentation Effects 0.000 claims abstract description 6
 230000035772 mutation Effects 0.000 claims description 14
 238000010586 diagram Methods 0.000 claims description 7
 230000002159 abnormal effect Effects 0.000 claims description 4
 238000004519 manufacturing process Methods 0.000 claims description 3
 230000003203 everyday Effects 0.000 claims description 2
 230000003247 decreasing Effects 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000006467 substitution reactions Methods 0.000 description 1
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
 G06F17/10—Complex mathematical operations
 G06F17/15—Correlation function computation including computation of convolution operations

 H—ELECTRICITY
 H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
 H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
 H02J3/00—Circuit arrangements for ac mains or ac distribution networks
 H02J3/26—Arrangements for eliminating or reducing asymmetry in polyphase networks

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
 Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
 Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
 Y02E40/50—Arrangements 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 threephase 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
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 distributiontransformation 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 threephase 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 threephase 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 U_{a1}，U_{a2}……U_{aN}And the phase voltage value of B is U_{b1}，U_{b2}……U_{bN}C phase voltage value is U_{c1}，U_{c2}……U_{cN}Calculating the distribution transformer outlet voltage U without considering the influence of threephase load unbalance by adopting the following formula_{i}：
Wherein U is_{ai}、U_{bi}、U_{ci}The voltage of the ith point A, B, C is 1 … N, U_{i}In order to not consider the distribution transformer outlet voltage at the ith point when the influence of threephase load unbalance is not considered, U can be obtained by an iteration method_{i}。
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 U_{T1i}The distribution transformer outlet voltage of the transformer area 2 is U_{T2i}Where i is 1 … … 96. Voltage difference value U of two station areas at corresponding time_{1i}；
U_{1i}＝U_{T1i}U_{T2i}，i＝1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations_{1i}Change amount of (Δ U)_{1i}；
ΔU_{1i}＝U_{1(i+1)}U_{1i}，i＝1……95
3) Carrying out normal test on the variable quantity through KolmogorovSmirnov, 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 k_{1}，k_{2}Since 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}),[k_{1},k_{2}),[k_{2},96]Respectively calculating the similarity S of the three voltage curves_{1}，S_{2}，S_{3}。
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 D_{1}，S_{2}，S_{3}And 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 databased 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 preprocessing by step B) according to the present invention;
FIG. 3 is a QQ 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 threephase voltage of the distribution transformer;
B. b, preprocessing the distribution transformer outlet voltage data acquired in the step A
The production management system records threephase 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 U_{a1}，U_{a2}……U_{aN}And the phase voltage value of B is U_{b1}，U_{b2}……U_{bN}C phase voltage value is U_{c1}，U_{c2}……U_{cN}Calculating the distribution transformer outlet voltage U without considering the influence of threephase load unbalance by adopting the following formula_{i}：
Wherein U is_{ai}、U_{bi}、U_{ci}The voltage of the ith point A, B, C is 1 … N, U_{i}In order to not consider the distribution transformer outlet voltage at the ith point when the influence of threephase load unbalance is not considered, U can be obtained by an iteration method_{i}。
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 U_{T1i}The distribution transformer outlet voltage of the transformer area 2 is U_{T2i}Where i is 1 … … 96. Voltage difference value U of two station areas at corresponding time_{1i}；
U_{1i}＝U_{T1i}U_{T2i}，i＝1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations_{1i}Change amount of (Δ U)_{1i}；
ΔU_{1i}＝U_{1(i+1)}U_{1i}，i＝1……95
3) Carrying out normal test on the variable quantity through KolmogorovSmirnov, 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 k_{1}，k_{2}Because 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 points_{1}),[k_{1},k_{2}),[k_{2},96]Respectively calculating the similarity S of the three voltage curves_{1}，S_{2}，S_{3}。
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 D_{1}，S_{2}，S_{3}And 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 threephase 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 U_{T1i}The distribution transformer outlet voltage of the transformer area 2 is U_{T2i}Where i is 1 … … 96. Voltage difference value U of two station areas at corresponding time_{1i}；
U_{1i}＝U_{T1i}U_{T2i}，i＝1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations_{1i}Change amount of (Δ U)_{1i}；
ΔU_{1i}＝U_{1(i+1)}U_{1i}，i＝1……95
3) The results of the normal test on the variation by KolmogorovSmirnov 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 positivetoodistribution, and indicates that the outlet voltage of one of the two platform regions has a mutation value. The results of the QQ plot shown in fig. 3 also show nonnormal 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 119 are the first segment, 2078 are the second segment, and 7996 are the third segment.
In step D, the positions of the two mutation points are respectively k_{1}＝20，k_{2}When 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 curves_{1}，S_{2}，S_{3}Wherein 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 threephase 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 U_{T1i}The distribution transformer outlet voltage of the transformer area 2 is U_{T2i}Where i is 1 … … 96, the voltage difference U between two stations at the corresponding time_{1i}：
U_{1i}＝U_{T1i}U_{T2i}，i＝1……96
2) Then calculating the variation of the voltage difference; voltage difference U between two stations_{1i}Change amount of (Δ U)_{1i}：
ΔU_{1i}＝U_{1(i+1)}U_{1i}，i＝1……95
3) Carrying out normal test on the variable quantity through KolmogorovSmirnov, 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 threephase 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 U_{a1}，U_{a2}……U_{aN}And the phase voltage value of B is U_{b1}，U_{b2}……U_{bN}C phase voltage value is U_{c1}，U_{c2}……U_{cN}The shadow without considering the threephase load unbalance is calculated by the following formulaTimeofresponse distribution transformer outlet voltage U_{i}：
Wherein U is_{ai}、U_{bi}、U_{ci}The voltage of the ith point A, B, C is 1 … N, U_{i}In order to not consider the distribution transformer outlet voltage at the ith point when the influence of threephase load unbalance is not considered, U can be obtained by an iteration method_{i}。
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 k_{1}，k_{2}Because the voltage curve has 96 points, the voltage curve is divided into three sections of [0, k ] according to the positions of two catastrophe points_{1}),[k_{1},k_{2}),[k_{2},96]Respectively calculating the similarity S of the three voltage curves_{1}，S_{2}，S_{3}。
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 D_{1}，S_{2}，S_{3}And taking the maximum value as the similarity of the voltage curves of the two transformer areas.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201910049892.7A CN109581130B (en)  20190118  20190118  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)  20190118  20190118  Voltage curve similarity calculation method considering voltage adjustment measure influence 
Publications (2)
Publication Number  Publication Date 

CN109581130A CN109581130A (en)  20190405 
CN109581130B true CN109581130B (en)  20201016 
Family
ID=65917014
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201910049892.7A Active CN109581130B (en)  20190118  20190118  Voltage curve similarity calculation method considering voltage adjustment measure influence 
Country Status (1)
Country  Link 

CN (1)  CN109581130B (en) 
Families Citing this family (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN112564110B (en) *  20210224  20210702  国网江西省电力有限公司电力科学研究院  Transformer area lowvoltage treatment method and system 
Family Cites Families (7)
Publication number  Priority date  Publication date  Assignee  Title 

CN104092481B (en) *  20140717  20160420  江苏林洋能源股份有限公司  A kind of by voltage characteristic differentiation Tai Qu and phase method for distinguishing 
JP6392155B2 (en) *  20150327  20180919  株式会社東芝  Power distribution monitoring and control device 
CN105184070A (en) *  20150831  20151223  华南理工大学  Voltage integral method based segmentation fitting method for calculating volttime characteristic curve 
CN107508297B (en) *  20170810  20190726  国家电网公司  A kind of verification of distribution topological structure and maintaining method 
CN108199375B (en) *  20180131  20210122  天津大学  Intelligent power distribution network topology identification method based on synchronous phasor measurement 
CN108564485A (en) *  20180416  20180921  国网河南省电力公司电力科学研究院  Lowvoltage platform area user's phase recognition methods based on voltage curve similarity analysis 
CN108600362A (en) *  20180419  20180928  南京林洋电力科技有限公司  A method of managing taiwan area archives automatically based on voltage characteristic 

2019
 20190118 CN CN201910049892.7A patent/CN109581130B/en active Active
Also Published As
Publication number  Publication date 

CN109581130A (en)  20190405 
Similar Documents
Publication  Publication Date  Title 

Wang et al.  Multicriteria 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&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)  Lowvoltage 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 technologybased abnormal diagnosis method for line loss rate of transformer area  
CN111091223A (en)  Distribution transformer shortterm 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 