CN113176458A - Low-voltage transformer area household relation identification method aiming at incomplete data - Google Patents

Low-voltage transformer area household relation identification method aiming at incomplete data Download PDF

Info

Publication number
CN113176458A
CN113176458A CN202110250758.0A CN202110250758A CN113176458A CN 113176458 A CN113176458 A CN 113176458A CN 202110250758 A CN202110250758 A CN 202110250758A CN 113176458 A CN113176458 A CN 113176458A
Authority
CN
China
Prior art keywords
low
user
identified
voltage
phase
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
Application number
CN202110250758.0A
Other languages
Chinese (zh)
Other versions
CN113176458B (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.)
Shenzhen Polytechnic
Original Assignee
Shenzhen Polytechnic
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 Shenzhen Polytechnic filed Critical Shenzhen Polytechnic
Priority to CN202110250758.0A priority Critical patent/CN113176458B/en
Publication of CN113176458A publication Critical patent/CN113176458A/en
Application granted granted Critical
Publication of CN113176458B publication Critical patent/CN113176458B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/06—Measuring real component; Measuring reactive component
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a low-voltage transformer area household relationship identification method aiming at incomplete data. The method comprises the following steps: calculating the active current time sequence data of each user to be identified and each phase bus at the low-voltage side of the low-voltage transformer area; classifying users to be identified based on the voltage time sequence curve correlation coefficient between the users to be identified to form an electric meter class library; solving a phase-family relation identification optimization model based on a kirchhoff current law to obtain a phase-family relation of a first low-voltage transformer area; and correcting the phase-to-user relationship of the first low-voltage transformer area based on the correlation characteristics between the users to be identified and the low-voltage side three-phase buses of the low-voltage transformer area to obtain a final phase-to-user relationship identification result of the low-voltage transformer area. The method can solve the problem of identification of the mutual relationship of the low-voltage distribution areas when the user data collected by the power grid company is incomplete, does not need to increase an acquisition terminal, and has the characteristics of low cost and high engineering application value.

Description

Low-voltage transformer area household relation identification method aiming at incomplete data
Technical Field
The invention relates to the technical field of electric power low-voltage distribution networks, in particular to a low-voltage transformer area household relationship identification method aiming at incomplete data.
Background
The current intelligent development of the low-voltage distribution network is limited by the loss or inaccuracy of the physical topological connection information of the low-voltage distribution network. The lack of accurate low-voltage topological relation can lead to the problems of difficult three-phase unbalance management, abnormal line loss statistics, untimely power failure rush repair and the like. Meanwhile, the intelligent electric energy meter and the low-voltage meter reading technology gradually realize full coverage in a low-voltage transformer area in recent years. The power grid company can remotely acquire the power utilization data and the power grid operation data of the user. The identification method based on data driving has the advantages of small modification amount, large input-output ratio and the like, and becomes an important technical direction for solving the problems of household variable relation verification and ammeter phase sequence identification of low-voltage distribution area topology identification. However, due to the problems of communication quality, manual error, electricity stealing and the like, the user data collected by the power grid company is not the complete electricity utilization data of the low-voltage distribution area.
The method is characterized in that a low-voltage distribution network phase-line-user identification method fusing user clustering, quadratic programming and a probability distribution model is provided by Thangjie and the like based on electric meter voltage and current data, but the method does not consider the situation that electric meter data collected by a concentrator is incomplete (Thangjie, Chua Yongzhi, Zhouyi, and the like, a data-driven low-voltage distribution network line-user relation identification method [ J ] power system automation 2020,44(11): 127-. Data incompleteness has universality in an actual low-voltage distribution network, so that the research on the low-voltage distribution network household identification method considering the data incompleteness is of great significance.
Disclosure of Invention
The invention aims to solve the problem of identification of the household relationship of the low-voltage distribution area when the user data is incomplete, and is beneficial to improving the operation benefit of a power grid enterprise and the customer satisfaction index.
According to the method, firstly, a phase-family relation identification optimization model based on the kirchhoff current law is used for determining a phase-family relation I of a low-voltage transformer area, then a multidimensional correction link is constructed based on correlation characteristics among users and between the users and a low-voltage side three-phase bus of the low-voltage transformer area, and the final phase-family relation of the low-voltage transformer area is determined. Compared with other identification methods, the method is based on kirchhoff current law and voltage correlation characteristics, a multidimensional correction link of the identification result is constructed, the problem of identification of the user relationship of the low-voltage transformer area when user data is incomplete can be solved, and the method has the characteristics of strong practicability, convenience in operation, reduction of labor cost of power companies and efficiency improvement.
The purpose of the invention is realized by at least one of the following technical solutions.
A low-voltage transformer area household relationship identification method aiming at incomplete data comprises the following steps:
s1, acquiring time sequence data of voltage and active power of a user to be identified and a low-voltage side three-phase bus of a low-voltage distribution room where the user to be identified is located, and calculating active current time sequence data of each user to be identified and each phase bus of the low-voltage side of the low-voltage distribution room;
s2, classifying users to be identified based on the voltage time sequence curve correlation coefficient among the users to be identified to form an electric meter class library, summing the active current values of the electric meters in the same class, and taking the value as the current value of the class;
s3, solving a kirchhoff current law-based phase-family relation identification optimization model to obtain a first low-voltage transformer area phase-family relation;
and S4, correcting the phase-to-user relationship of the first low-voltage transformer area based on the correlation characteristics between the users to be identified and the low-voltage side three-phase buses of the low-voltage transformer area to obtain a final phase-to-user relationship identification result of the low-voltage transformer area.
Further, in step S1, the active current timing sequence data of each user to be identified and each phase bus at the low-voltage side of the low-voltage transformer area are calculated as follows:
Figure BDA0002965961120000021
wherein, Ii,t、Pi,tAnd ui,tRespectively representing the active current value, the active power and the voltage of the user i to be identified at the moment t;
Figure BDA0002965961120000022
and
Figure BDA0002965961120000023
respectively represents the low-voltage side of the distribution transformer at the time t
Figure BDA0002965961120000024
The active current value, active power and voltage of the phase bus; i 1,2, …, M, T1, 2, …, T,
Figure BDA0002965961120000025
m represents a meter reading directory of a low-voltage transformer area and comprisesT represents the total number of time segments of the data acquisition cycle.
Further, step S2 specifically includes the following steps:
s2.1, calculating a voltage curve correlation coefficient matrix R between users to be identified included in a meter reading directory of the low-voltage transformer area, wherein the u-th row element in the matrix R is a voltage time sequence curve correlation coefficient between the users to be identified u and all the users to be identified, and the method specifically comprises the following steps:
Figure BDA0002965961120000026
m represents the total number of users to be identified in a meter reading directory of the low-voltage transformer area; theta is a user to be identified included in the meter reading catalog of the low-voltage transformer area; r isuvThe voltage time sequence correlation coefficient representing the time sequence correlation coefficient between the user u to be identified and the user v to be identified is as follows:
Figure BDA0002965961120000031
in the formula (I), the compound is shown in the specification,
Figure BDA0002965961120000032
and
Figure BDA0002965961120000033
voltage values of a user u to be identified and a user v to be identified at a time T are respectively, u and v are respectively equal to theta, and T is 1,2, … and T;
s2.2, based on the matrix R, classifying each user to be identified and the user to be identified with the maximum voltage curve correlation except the user to be identified into one class, and obtaining M double-table classifications in total;
s2.3, merging the classifications including the same user to obtain a user class library omega including N classesclaAnd summing the active current values of the electric meters under the same category, and taking the value as the current value of the category.
Further, in step S3, the objective function of the identification optimization model based on the phase-user relationship of kirchhoff' S current law is the minimum sum of absolute values of the difference between the sum of the active current flowing out from the head end of the low-voltage side bus of each phase and the active current belonging to the category of the electric meter of the phase in the data acquisition period:
Figure BDA0002965961120000034
in the formula, T represents the total time period number of the data acquisition period;
Figure BDA0002965961120000035
indicating the low-voltage side of the distribution transformation at time t
Figure BDA0002965961120000036
An active current value of the phase bus;
Figure BDA0002965961120000037
a binary variable representing the attribution relationship between the electric meter class g and each phase of the target transformer area
Figure BDA0002965961120000038
When equal to 1, it indicates that the electric meter class g belongs to
Figure BDA0002965961120000039
Otherwise, the electric meter class g does not belong to
Figure BDA00029659611200000310
Phase (1); k is the total number of the categories of the electric meters;
the constraint condition of the phase-user relationship identification optimization model based on kirchhoff current law needs to ensure that each ammeter class cannot belong to a plurality of phase lines at the same time, as shown in the following formula:
Figure BDA00029659611200000311
further, in step S3, the optimization model is identified by solving the correlation based on kirchhoff' S current law with CPLEX.
Further, in step S4, the specific steps are as follows:
s4.1, screening the first low-voltage transformer area correlation theta1Identifying the wrong suspected user by the middle phase sequence;
s4.2, correcting results based on the correlation characteristics between the user to be identified and the low-voltage side buses of the low-voltage distribution area;
s4.3, correcting results based on the correlation characteristics among the electric meters;
and S4.4, determining the inter-house relationship of the low-voltage transformer area by combining the correlation characteristics between the user to be identified and the low-voltage side bus of the low-voltage transformer area and the correlation characteristics between the users to be identified, and finishing the inter-house relationship identification of the low-voltage transformer area.
Further, in step S4.1, the specific steps are as follows:
s4.1.1, calculating the first low-voltage transformer area household relation theta1In the method, the voltage curve correlation coefficients between users to be identified belonging to the same phase sequence are calculated according to the following formula:
Figure BDA0002965961120000041
in the formula ui,t、uj,tRespectively representing the voltage values of a user i to be identified and a user j to be identified at the moment t, wherein i and j are equal to [1,2, …, M]M is the total number of users to be identified, T is 1, 2.
S4.1.2, calculating the average value of the voltage curve correlation coefficient of each user to be identified and the users to be identified belonging to the same phase sequence,
Figure BDA0002965961120000042
in the formula (I), the compound is shown in the specification,
Figure BDA0002965961120000043
is attributed to
Figure BDA0002965961120000044
The user i to be identified of the phase is associated with the user i to be identified in the same phase sequenceIdentifying the mean value, r, of the correlation coefficient of the voltage curve of the userigFor the voltage curve correlation coefficient of the user i to be identified and the user c to be identified,
Figure BDA0002965961120000045
is attributed to
Figure BDA0002965961120000046
The set of users to be identified of the facies,
Figure BDA0002965961120000047
is attributed to
Figure BDA0002965961120000048
The total number of the users to be identified;
s4.1.3, averaging the average value of the voltage curve correlation coefficients of the users to be identified belonging to the same phase sequence and the users to be identified belonging to the phase sequence,
Figure BDA0002965961120000049
in the formula (I), the compound is shown in the specification,
Figure BDA00029659611200000410
is attributed to
Figure BDA00029659611200000411
The mean value of the correlation coefficient mean values of voltage curves between the phase users to be identified;
if it is
Figure BDA00029659611200000412
Is less than
Figure BDA00029659611200000413
The user i to be identified is set as the suspected user with the wrong phase sequence identification.
Further, in step S4.2, the specific steps are as follows:
s4.2.1, respectively calculating voltage time sequence curve correlation coefficients of all suspected users in the step S4.1 and the low-voltage side three-phase bus of the low-voltage transformer area, and taking the bus phase sequence with the maximum voltage time sequence curve correlation coefficient value as a first phase sequence result of the suspected users;
s4.2.2, replace the correlation theta of the first low-voltage transformer area with the first phase sequence result of the suspected user in step S4.2.11The phase sequence of the suspected user in the second low-voltage transformer area obtains the phase-to-user relationship theta of the second low-voltage transformer area2。
Further, in step S4.3, the specific steps are as follows:
s4.3.1, relationship theta between the first low-voltage transformer area and the user1Removing the suspected user in the step S4.1 to obtain the relation theta of the third low-voltage transformer area with the user3；
S4.3.2, calculating the relationship theta between each suspected user and the third low-voltage transformer area3The voltage curve correlation coefficients of all users to be identified in each phase are obtained, so that each suspected user can obtain 3 voltage curve correlation coefficient vectors corresponding to A, B, C three phases, all elements in the 3 vectors are respectively averaged, and the phase sequence to which the voltage curve correlation coefficient vector with the largest average value belongs is taken as a second phase sequence result of the suspected user;
s4.3.3, replacing the correlation theta of the first low-voltage transformer area with the second order result of the suspected user in step S4.3.21The phase sequence of the suspected user in (1) is obtained to obtain the phase-user relationship theta of the fourth low-voltage transformer area4。
Further, in step S4.4, the specific steps are as follows:
s4.4.1, screening out the correlation theta of the second low-voltage transformer area2The relation theta with the fourth low-voltage transformer area4Obtaining a user set zeta by users to be identified with inconsistent middle phase sequences;
s4.4.2, averaging the voltage time sequence curves of the users to be identified, sorting the users to be identified according to the voltage average value from large to small, and setting a threshold coefficient tau to be in the range of 0,0.5]Extracting the top of the sequencing result of the user to be identified
Figure BDA0002965961120000051
A user to be identifiedForming a set d as a user set close to the head end of the low-voltage side of the low-voltage transformer area, wherein M is the total number of users to be identified;
s4.4.3, the phase sequence of the users to be identified in the user set belonging to the user set near the head end of the low-voltage side of the low-voltage distribution area in the user set zeta is in accordance with the second low-voltage distribution area phase relationship theta2If the user phase sequence relation cannot be updated, the user phase sequence relation is updated according to the fourth low-voltage transformer area phase sequence relation theta4And updating the user phase sequence relation to obtain the final low-voltage transformer area phase relationship.
The invention has the beneficial effects that:
(1) the identification result multidimensional correction link is constructed based on the kirchhoff current law and the voltage correlation characteristic, is suitable for identifying the phase relationship of the low-voltage transformer area when user data is incomplete, and is beneficial to improving the accuracy of identifying the phase relationship of the low-voltage transformer area in practical application;
(2) the acquisition terminal does not need to be added in the low-voltage distribution network, so the method has the characteristics of low cost and small engineering quantity.
Drawings
FIG. 1 is a flow chart of a method for identifying a correlation between low-voltage transformer areas with incomplete data;
fig. 2 is a schematic diagram of a network connection in a real cell.
Detailed Description
The following description of the embodiments of the present invention is provided in connection with the accompanying drawings and examples.
Example (b):
a method for identifying a subscriber relationship in a low-voltage distribution area with incomplete data, as shown in fig. 1, includes the following steps:
s1, acquiring voltage and active power time sequence data of the user to be identified and the low-voltage side three-phase bus of the low-voltage distribution room where the user to be identified is located, and calculating active current time sequence data of each user to be identified and each phase bus of the low-voltage side of the low-voltage distribution room, wherein the active current time sequence data are as follows:
Figure BDA0002965961120000061
wherein, Ii,t、Pi,tAnd ui,tRespectively representing the active current value, the active power and the voltage of the user i to be identified at the moment t;
Figure BDA0002965961120000062
and
Figure BDA0002965961120000063
respectively represents the low-voltage side of the distribution transformer at the time t
Figure BDA0002965961120000064
The active current value, active power and voltage of the phase bus; i 1,2, …, M, T1, 2, …, T,
Figure BDA0002965961120000065
m represents the total number of users to be identified in the meter reading directory of the low-voltage transformer area, and T represents the total time period number of the data acquisition period.
In this embodiment, as shown in fig. 2, the low-voltage transformer area includes 29 households of a-phase electric meters, 32 households of B-phase electric meters, 15 households of C-phase electric meters, and 8 households of three-phase users. In fig. 2, S indicates a single-phase user, and T indicates a three-phase user. And collecting 2-day voltage time sequence data of the low-voltage distribution transformer low-voltage side bus and all the electric meters, wherein 6 tables affected by communication fail to collect information, and the 7 tables are S1, S2, S7, S8, S9 and S10 respectively.
S2, classifying users to be identified based on the voltage time sequence curve correlation coefficient between the users to be identified to form an electric meter class library, summing the active current values of the electric meters in the same class, and taking the value as the current value of the class, wherein the method specifically comprises the following steps:
s2.1, calculating a voltage curve correlation coefficient matrix R between users to be identified included in a meter reading directory of the low-voltage transformer area, wherein the u-th row element in the matrix R is a voltage time sequence curve correlation coefficient between the users to be identified u and all the users to be identified, and the method specifically comprises the following steps:
Figure BDA0002965961120000066
m represents the total number of users to be identified in a meter reading directory of the low-voltage transformer area; theta is a user to be identified included in the meter reading catalog of the low-voltage transformer area; r isuvThe voltage time sequence correlation coefficient representing the time sequence correlation coefficient between the user u to be identified and the user v to be identified is as follows:
Figure BDA0002965961120000067
in the formula (I), the compound is shown in the specification,
Figure BDA0002965961120000068
and
Figure BDA0002965961120000069
voltage values of a user u to be identified and a user v to be identified at a time T are respectively, u and v are respectively equal to theta, and T is 1,2, … and T;
s2.2, based on the matrix R, classifying each user to be identified and the user to be identified with the maximum voltage curve correlation except the user to be identified into one class, and obtaining M double-table classifications in total;
s2.3, merging the classifications including the same user to obtain a user class library omega including N classesclaAnd summing the active current values of the electric meters under the same category, and taking the value as the current value of the category.
In this embodiment, the obtained user category library is shown in table 1.
TABLE 1 user Category Bank
User categories User categories User categories
T11,S17,S18 T12,S11,S12 S3,S19,S20,S21
S25,S26,S27,S31 S13,S14,S15,S16 S6,T13,S23,S24
S28,S29,S30 S32,S33,T62,T72,T82,S45,S47 S34,T33,T43,T53
S35,S39,S40,S41,S42,S43 S36,T22,T32,T42,T52 S37,T23
S38,T21,T31,T41,T51 S51,S54,S55,S56,S57 T63,T73,T83,S46
T61,T71,T81,S44 S64,S65 S49,S50
S58,S59,S60,S61 S66,S67,S68,S69 S52,S53,S70,S71
S62,S63 S4,S5,S22,S48 S72,S73,S74,S75,S76
S3, solving a kirchhoff current law-based phase-family relation identification optimization model to obtain a first low-voltage transformer area phase-family relation;
the objective function of the phase-user relationship identification optimization model based on the kirchhoff current law is that the sum of absolute values of the difference values of the sum of the active current flowing out from the head end of the low-voltage side bus of each phase and the active current belonging to the class of the ammeter of the phase in a data acquisition period is minimum:
Figure BDA0002965961120000071
in the formula, T represents the total time period number of the data acquisition period;
Figure BDA0002965961120000072
indicating the low-voltage side of the distribution transformation at time t
Figure BDA0002965961120000073
An active current value of the phase bus;
Figure BDA0002965961120000074
a binary variable representing the attribution relationship between the electric meter class g and each phase of the target transformer area
Figure BDA0002965961120000075
When equal to 1, it indicates that the electric meter class g belongs to
Figure BDA0002965961120000076
Otherwise, the electric meter class g does not belong to
Figure BDA0002965961120000077
Phase (1); k is the total number of the categories of the electric meters;
the constraint condition of the phase-user relationship identification optimization model based on kirchhoff current law needs to ensure that each ammeter class cannot belong to a plurality of phase lines at the same time, as shown in the following formula:
Figure BDA0002965961120000078
and solving a phase-user relationship identification optimization model based on the kirchhoff current law through CPLEX.
In this embodiment, the relationship between the first low-voltage transformer areas and the subscriber is shown in the following table.
TABLE 2 Low Voltage distribution area household relationship one
Figure BDA0002965961120000081
S4, correcting the first low-voltage transformer area phase-to-phase relationship based on the correlation characteristics between the users to be identified and the low-voltage side three-phase buses of the low-voltage transformer area to obtain a final low-voltage transformer area phase-to-phase relationship identification result, and the method specifically comprises the following steps:
s4.1, screening the first low-voltage transformer area correlation theta1The method comprises the following steps of:
s4.1.1, calculating the first low-voltage transformer area household relation theta1In the method, the voltage curve correlation coefficients between users to be identified belonging to the same phase sequence are calculated according to the following formula:
Figure BDA0002965961120000082
in the formula ui,t、uj,tRespectively representing the voltage values of a user i to be identified and a user j to be identified at the moment t, wherein i and j are equal to [1,2, …, M]M is the total number of users to be identified, T is 1, 2.
S4.1.2, calculating the average value of the voltage curve correlation coefficient of each user to be identified and the users to be identified belonging to the same phase sequence,
Figure BDA0002965961120000083
in the formula (I), the compound is shown in the specification,
Figure BDA0002965961120000084
is attributed to
Figure BDA0002965961120000085
Mean value r of the correlation coefficients of the voltage curves of the phase subscriber i to be identified and the subscribers to be identified belonging to the same phase sequenceigFor the voltage curve correlation coefficient of the user i to be identified and the user c to be identified,
Figure BDA0002965961120000086
is attributed to
Figure BDA0002965961120000087
The set of users to be identified of the facies,
Figure BDA0002965961120000088
is attributed to
Figure BDA0002965961120000089
The total number of the users to be identified;
s4.1.3, averaging the average value of the voltage curve correlation coefficients of the users to be identified belonging to the same phase sequence and the users to be identified belonging to the phase sequence,
Figure BDA0002965961120000091
in the formula (I), the compound is shown in the specification,
Figure BDA0002965961120000092
is attributed to
Figure BDA0002965961120000093
The mean value of the correlation coefficient mean values of voltage curves between the phase users to be identified;
if it is
Figure BDA0002965961120000094
Is less than
Figure BDA0002965961120000095
The user i to be identified is set as the suspected user with the wrong phase sequence identification.
S4.2, based on the correlation characteristic correction result between the user to be identified and the low-voltage side bus of the low-voltage area, the method specifically comprises the following steps:
s4.2.1, respectively calculating voltage time sequence curve correlation coefficients of all suspected users in the step S4.1 and the low-voltage side three-phase bus of the low-voltage transformer area, and taking the bus phase sequence with the maximum voltage time sequence curve correlation coefficient value as a first phase sequence result of the suspected users;
s4.2.2, replace the correlation theta of the first low-voltage transformer area with the first phase sequence result of the suspected user in step S4.2.11The phase sequence of the suspected user in the second low-voltage transformer area obtains the phase-to-user relationship theta of the second low-voltage transformer area2。
S4.3, correcting results based on the correlation characteristics among the electric meters, and specifically comprising the following steps:
s4.3.1, relationship theta between the first low-voltage transformer area and the user1Removing the suspected user in the step S4.1 to obtain the relation theta of the third low-voltage transformer area with the user3；
S4.3.2, calculating the relationship theta between each suspected user and the third low-voltage transformer area3The voltage curve correlation coefficients of all users to be identified in each phase are obtained, so that each suspected user can obtain 3 voltage curve correlation coefficient vectors corresponding to A, B, C three phases, all elements in the 3 vectors are respectively averaged, and the phase sequence to which the voltage curve correlation coefficient vector with the largest average value belongs is taken as a second phase sequence result of the suspected user;
s4.3.3, replacing the correlation theta of the first low-voltage transformer area with the second order result of the suspected user in step S4.3.21The phase sequence of the suspected user in (1) is obtained to obtain the phase-user relationship theta of the fourth low-voltage transformer area4。
S4.4, determining the low-voltage transformer area household relation by combining the correlation characteristic between the user to be identified and the low-voltage side bus of the low-voltage transformer area and the correlation characteristic between the users to be identified, and completing the identification of the low-voltage transformer area household relation, wherein the specific steps are as follows:
s4.4.1, screening out the correlation theta of the second low-voltage transformer area2The relation theta with the fourth low-voltage transformer area4Obtaining a user set zeta by users to be identified with inconsistent middle phase sequences;
s4.4.2, averaging the voltage time sequence curves of the users to be identified, sorting the users to be identified according to the voltage average value from large to small, and setting a threshold coefficient tau to be in the range of 0,0.5]Extracting the top of the sequencing result of the user to be identified
Figure BDA0002965961120000096
Forming a set d by the users to be identified as a user set close to the head end of the low-voltage side of the low-voltage transformer area, wherein M is the total number of the users to be identified;
s4.4.3, the phase sequence of the users to be identified in the user set belonging to the user set near the head end of the low-voltage side of the low-voltage distribution area in the user set zeta is in accordance with the second low-voltage distribution area phase relationship theta2If the user phase sequence relation cannot be updated, the user phase sequence relation is updated according to the fourth low-voltage transformer area phase sequence relation theta4And updating the user phase sequence relation to obtain the final low-voltage transformer area phase relationship.
In this embodiment, the relationship between the subscriber identities of the low-voltage distribution area is finally shown in the following table.
TABLE 3 end user phase sequence results
Figure BDA0002965961120000101
As can be seen from fig. 2, the results shown in table 3 correctly reflect the correlation of the station area.
To sum up, the above embodiments illustrate the effectiveness of the low voltage transformer area household relationship identification method based on the voltage correlation characteristic provided by the present invention.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents and are intended to be included in the scope of the present invention.

Claims (10)

1. A low-voltage transformer area household relationship identification method aiming at incomplete data is characterized by comprising the following steps:
s1, acquiring time sequence data of voltage and active power of a user to be identified and a low-voltage side three-phase bus of a low-voltage distribution room where the user to be identified is located, and calculating active current time sequence data of each user to be identified and each phase bus of the low-voltage side of the low-voltage distribution room;
s2, classifying users to be identified based on the voltage time sequence curve correlation coefficient among the users to be identified to form an electric meter class library, summing the active current values of the electric meters in the same class, and taking the value as the current value of the class;
s3, solving a kirchhoff current law-based phase-family relation identification optimization model to obtain a first low-voltage transformer area phase-family relation;
and S4, correcting the phase-to-user relationship of the first low-voltage transformer area based on the correlation characteristics between the users to be identified and the low-voltage side three-phase buses of the low-voltage transformer area to obtain a final phase-to-user relationship identification result of the low-voltage transformer area.
2. The method for identifying the phase-user relationship of the low-voltage transformer area with incomplete data according to claim 1, wherein in step S1, the active current timing sequence data of each user to be identified and each phase bus at the low-voltage side of the low-voltage transformer area are calculated as follows:
Figure FDA0002965961110000011
wherein, Ii,t、Pi,tAnd ui,tRespectively representing the active current value, the active power and the voltage of the user i to be identified at the moment t;
Figure FDA0002965961110000012
and
Figure FDA0002965961110000013
respectively represents the distribution low voltage at the time tSide wall
Figure FDA0002965961110000014
The active current value, active power and voltage of the phase bus; i 1,2, …, M, T1, 2, …, T,
Figure FDA0002965961110000015
m represents the total number of users to be identified in the meter reading directory of the low-voltage transformer area, and T represents the total time period number of the data acquisition period.
3. The method for identifying the subscriber relationship of the low voltage transformer area with incomplete data as claimed in claim 1, wherein step S2 specifically includes the following steps:
s2.1, calculating a voltage curve correlation coefficient matrix R between users to be identified included in a meter reading directory of the low-voltage transformer area, wherein the u-th row element in the matrix R is a voltage time sequence curve correlation coefficient between the users to be identified u and all the users to be identified, and the method specifically comprises the following steps:
Figure FDA0002965961110000016
m represents the total number of users to be identified in a meter reading directory of the low-voltage transformer area; theta is a user to be identified included in the meter reading catalog of the low-voltage transformer area; r isuvThe voltage time sequence correlation coefficient representing the time sequence correlation coefficient between the user u to be identified and the user v to be identified is as follows:
Figure FDA0002965961110000021
in the formula (I), the compound is shown in the specification,
Figure FDA0002965961110000022
and
Figure FDA0002965961110000023
respectively being a user u to be identified and a user v to be identified at the momentT, u and v ∈ Θ, T being 1,2, …, T;
s2.2, based on the matrix R, classifying each user to be identified and the user to be identified with the maximum voltage curve correlation except the user to be identified into one class, and obtaining M double-table classifications in total;
s2.3, merging the classifications including the same user to obtain a user class library omega including K classesclaAnd summing the active current values of the electric meters under the same category, and taking the value as the current value of the category.
4. The method for identifying the phase relationship of the low-voltage transformer area with incomplete data according to claim 1, wherein in step S3, the sum of absolute values of the difference between the active current flowing from the head end of the low-voltage side bus of each phase and the sum of the active currents belonging to the electric meter categories of the phase is the smallest based on the objective function of the phase relationship identification optimization model based on kirchhoff' S current law:
Figure FDA0002965961110000024
in the formula, T represents the total time period number of the data acquisition period;
Figure FDA0002965961110000025
indicating the low-voltage side of the distribution transformation at time t
Figure FDA0002965961110000026
An active current value of the phase bus;
Figure FDA0002965961110000027
a binary variable representing the attribution relationship between the electric meter class g and each phase of the target transformer area
Figure FDA0002965961110000028
When equal to 1, it indicates that the electric meter class g belongs to
Figure FDA0002965961110000029
Otherwise, the electric meter class g does not belong to
Figure FDA00029659611100000210
Phase (1); k is the total number of the categories of the electric meters;
the constraint condition of the phase-user relationship identification optimization model based on kirchhoff current law needs to ensure that each ammeter class cannot belong to a plurality of phase lines at the same time, as shown in the following formula:
Figure FDA00029659611100000211
5. the method for identifying phase relationships in low-voltage transformer areas with incomplete data according to claim 4, wherein in step S3, the phase relationship identification optimization model based on kirchhoff current law is solved through CPLEX.
6. The method for identifying the subscriber relationship of the low-voltage transformer area with incomplete data as claimed in claim 1, wherein in step S4, the specific steps are as follows:
s4.1, screening the first low-voltage transformer area correlation theta1Identifying the wrong suspected user by the middle phase sequence;
s4.2, correcting results based on the correlation characteristics between the user to be identified and the low-voltage side buses of the low-voltage distribution area;
s4.3, correcting results based on the correlation characteristics among the electric meters;
and S4.4, determining the inter-house relationship of the low-voltage transformer area by combining the correlation characteristics between the user to be identified and the low-voltage side bus of the low-voltage transformer area and the correlation characteristics between the users to be identified, and finishing the inter-house relationship identification of the low-voltage transformer area.
7. The method for identifying the subscriber relationship in the low-voltage transformer area with incomplete data as claimed in claim 6, wherein in step S4.1, the specific steps are as follows:
s4.1.1, calculating the first low-voltage transformer area household relation theta1In the method, the voltage curve correlation coefficients between users to be identified belonging to the same phase sequence are calculated according to the following formula:
Figure FDA0002965961110000031
in the formula ui,t、uj,tRespectively representing the voltage values of a user i to be identified and a user j to be identified at the moment t, wherein i and j are equal to [1,2, …, M]M is the total number of users to be identified, T is 1, 2.
S4.1.2, calculating the average value of the voltage curve correlation coefficient of each user to be identified and the users to be identified belonging to the same phase sequence,
Figure FDA0002965961110000032
in the formula (I), the compound is shown in the specification,
Figure FDA0002965961110000033
is attributed to
Figure FDA0002965961110000034
Mean value r of the correlation coefficients of the voltage curves of the phase subscriber i to be identified and the subscribers to be identified belonging to the same phase sequenceicFor the voltage curve correlation coefficient of the user i to be identified and the user c to be identified,
Figure FDA0002965961110000035
is attributed to
Figure FDA0002965961110000036
The set of users to be identified of the facies,
Figure FDA0002965961110000037
is attributed to
Figure FDA0002965961110000038
The total number of the users to be identified;
s4.1.3, averaging the average value of the voltage curve correlation coefficients of the users to be identified belonging to the same phase sequence and the users to be identified belonging to the phase sequence,
Figure FDA0002965961110000039
in the formula (I), the compound is shown in the specification,
Figure FDA00029659611100000310
is attributed to
Figure FDA00029659611100000311
The mean value of the correlation coefficient mean values of voltage curves between the phase users to be identified;
if it is
Figure FDA00029659611100000312
Is less than
Figure FDA00029659611100000313
The user i to be identified is set as the suspected user with the wrong phase sequence identification.
8. The method for identifying the subscriber relationship in the low-voltage transformer area with incomplete data as claimed in claim 6, wherein in step S4.2, the specific steps are as follows:
s4.2.1, respectively calculating voltage time sequence curve correlation coefficients of all suspected users in the step S4.1 and the low-voltage side three-phase bus of the low-voltage transformer area, and taking the bus phase sequence with the maximum voltage time sequence curve correlation coefficient value as a first phase sequence result of the suspected user;
s4.2.2, replace the correlation theta of the first low-voltage transformer area with the first phase sequence result of the suspected user in step S4.2.11The phase sequence of the suspected user in the second low-voltage transformer area obtains the phase-to-user relationship theta of the second low-voltage transformer area2。
9. The method for identifying the subscriber relationship in the low-voltage transformer area with incomplete data as claimed in claim 6, wherein in step S4.3, the specific steps are as follows:
s4.3.1, relationship theta between the first low-voltage transformer area and the user1Removing the suspected user in the step S4.1 to obtain the relation theta of the third low-voltage transformer area with the user3；
S4.3.2, calculating the relationship theta between each suspected user and the third low-voltage transformer area3Obtaining 3 voltage curve correlation coefficient vectors corresponding to A, B, C three phases by each suspected user, respectively averaging all elements in the 3 vectors, and taking the phase sequence to which the voltage curve correlation coefficient vector with the largest average value belongs as a second phase sequence result of the suspected user;
s4.3.3, replacing the correlation theta of the first low-voltage transformer area with the second order result of the suspected user in step S4.3.21The phase sequence of the suspected user in (1) is obtained to obtain the phase-user relationship theta of the fourth low-voltage transformer area4。
10. The low-voltage transformer area user relationship identification method for incomplete data according to any one of claims 1 to 9, characterized in that in step S4.4, the specific steps are as follows:
s4.4.1, screening out the correlation theta of the second low-voltage transformer area2The relation theta with the fourth low-voltage transformer area4Obtaining a user set zeta by users to be identified with inconsistent middle phase sequences;
s4.4.2, averaging the voltage time sequence curves of the users to be identified, sorting the users to be identified according to the voltage average value from large to small, and setting a threshold coefficient tau to be in the range of 0,0.5]Extracting the top of the sequencing result of the user to be identified
Figure FDA0002965961110000041
The users to be identified form a set d as a user set close to the head end of the low-voltage side of the low-voltage area, wherein M is the user set to be identifiedThe total number of users;
s4.4.3, the phase sequence of the users to be identified in the user set belonging to the user set near the head end of the low-voltage side of the low-voltage distribution area in the user set zeta is in accordance with the second low-voltage distribution area phase relationship theta2If the user phase sequence relation cannot be updated, the user phase sequence relation is updated according to the fourth low-voltage transformer area phase sequence relation theta4And updating the user phase sequence relation to obtain the final low-voltage transformer area phase relationship.
CN202110250758.0A 2021-03-08 2021-03-08 Low-voltage transformer area household relation identification method aiming at incomplete data Active CN113176458B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110250758.0A CN113176458B (en) 2021-03-08 2021-03-08 Low-voltage transformer area household relation identification method aiming at incomplete data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110250758.0A CN113176458B (en) 2021-03-08 2021-03-08 Low-voltage transformer area household relation identification method aiming at incomplete data

Publications (2)

Publication Number Publication Date
CN113176458A true CN113176458A (en) 2021-07-27
CN113176458B CN113176458B (en) 2022-08-30

Family

ID=76921878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110250758.0A Active CN113176458B (en) 2021-03-08 2021-03-08 Low-voltage transformer area household relation identification method aiming at incomplete data

Country Status (1)

Country Link
CN (1) CN113176458B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111505446A (en) * 2020-05-25 2020-08-07 广州市奔流电力科技有限公司 Method, device and equipment for identifying subscriber relationship of platform area vacant house subscriber
CN115577993A (en) * 2022-12-09 2023-01-06 江苏瑞电智芯信息科技有限公司 Transformer area household change identification method based on time sequence matching
CN117131355A (en) * 2023-06-19 2023-11-28 广州番禺职业技术学院 Method and device for identifying looks at home

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7095786B1 (en) * 2003-01-11 2006-08-22 Neo Magic Corp. Object tracking using adaptive block-size matching along object boundary and frame-skipping when object motion is low
CN111505446A (en) * 2020-05-25 2020-08-07 广州市奔流电力科技有限公司 Method, device and equipment for identifying subscriber relationship of platform area vacant house subscriber
CN111505443A (en) * 2020-05-13 2020-08-07 广州市奔流电力科技有限公司 Low-voltage transformer area line-to-user relationship identification method and device and computer equipment
CN111541301A (en) * 2020-04-27 2020-08-14 杭州沃瑞电力科技有限公司 Power distribution network line variation relation checking method based on fluctuation rate
CN111817291A (en) * 2020-06-16 2020-10-23 华南理工大学 Low-voltage transformer area household relation identification method based on voltage correlation characteristics
CN111881189A (en) * 2020-08-04 2020-11-03 广东电网有限责任公司计量中心 Method and system for verifying station area user-to-user relationship
CN112132283A (en) * 2020-10-20 2020-12-25 江苏其厚智能电气设备有限公司 Non-signal injection type user variable topological relation identification method based on genetic algorithm
CN112162136A (en) * 2020-09-14 2021-01-01 佰聆数据股份有限公司 Electricity larceny prevention method and system applied to power distribution network fusion terminal
CN112380496A (en) * 2020-09-18 2021-02-19 青岛鼎信通讯股份有限公司 Pearson correlation coefficient-based anti-electricity-stealing analysis method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7095786B1 (en) * 2003-01-11 2006-08-22 Neo Magic Corp. Object tracking using adaptive block-size matching along object boundary and frame-skipping when object motion is low
CN111541301A (en) * 2020-04-27 2020-08-14 杭州沃瑞电力科技有限公司 Power distribution network line variation relation checking method based on fluctuation rate
CN111505443A (en) * 2020-05-13 2020-08-07 广州市奔流电力科技有限公司 Low-voltage transformer area line-to-user relationship identification method and device and computer equipment
CN111505446A (en) * 2020-05-25 2020-08-07 广州市奔流电力科技有限公司 Method, device and equipment for identifying subscriber relationship of platform area vacant house subscriber
CN111817291A (en) * 2020-06-16 2020-10-23 华南理工大学 Low-voltage transformer area household relation identification method based on voltage correlation characteristics
CN111881189A (en) * 2020-08-04 2020-11-03 广东电网有限责任公司计量中心 Method and system for verifying station area user-to-user relationship
CN112162136A (en) * 2020-09-14 2021-01-01 佰聆数据股份有限公司 Electricity larceny prevention method and system applied to power distribution network fusion terminal
CN112380496A (en) * 2020-09-18 2021-02-19 青岛鼎信通讯股份有限公司 Pearson correlation coefficient-based anti-electricity-stealing analysis method
CN112132283A (en) * 2020-10-20 2020-12-25 江苏其厚智能电气设备有限公司 Non-signal injection type user variable topological relation identification method based on genetic algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CUIZHE KUANG,ET AL: "Multi-energy Flow Regulation and Operation Method of Distribution Network Based on Smart City Energy Internet", 《2020 INTERNATIONAL CONFERENCE ON COMMUNICATIONS, INFORMATION SYSTEM AND COMPUTER ENGINEERING (CISCE)》 *
尹飞等: "基于二次聚类的低压台区负荷特性分析", 《电力技术交流》 *
张勇等: "基于AHP-GCA及多元线性回归模型的电压合格率预测", 《2015年全国智能电网用户端能源管理学术年会论文集》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111505446A (en) * 2020-05-25 2020-08-07 广州市奔流电力科技有限公司 Method, device and equipment for identifying subscriber relationship of platform area vacant house subscriber
CN115577993A (en) * 2022-12-09 2023-01-06 江苏瑞电智芯信息科技有限公司 Transformer area household change identification method based on time sequence matching
CN117131355A (en) * 2023-06-19 2023-11-28 广州番禺职业技术学院 Method and device for identifying looks at home
CN117131355B (en) * 2023-06-19 2024-06-07 广州番禺职业技术学院 Method and device for identifying looks at home

Also Published As

Publication number Publication date
CN113176458B (en) 2022-08-30

Similar Documents

Publication Publication Date Title
CN113176458B (en) Low-voltage transformer area household relation identification method aiming at incomplete data
WO2021253806A1 (en) Voltage association characteristic-based method for identifying phase sequence-user relation in low-voltage transformer area
CN110376469B (en) Method and device for identifying relationship between phase sequence and line subscriber of transformer area
CN111505433B (en) Low-voltage transformer area indoor variable relation error correction and phase identification method
CN112329806B (en) Ammeter clustering method for topology identification of low-voltage transformer area
CN113078630B (en) Low-voltage distribution network topology identification method based on real-time measurement data
CN111476427A (en) Low-voltage distribution area topology identification method and identification device
CN112966219A (en) Method, system, equipment and medium for identifying relationship between household meter and meter box
Zhou et al. Consumer phase identification under incomplete data condition with dimensional calibration
CN111723339A (en) Method for identifying low-voltage hitching of transformer area based on trend similarity and distance measurement
CN110264032A (en) A kind of acquisition methods of direct current distribution Plan mode result
CN112198375B (en) Phase identification method, device, equipment and storage medium for single-phase user
CN111835006A (en) Low-voltage transformer area topology identification method based on voltage curve and least square
CN112636335B (en) Low-voltage distribution network user transformer relation verification method and system based on voltage data decorrelation selection and correlation analysis
CN112508254B (en) Method for determining investment prediction data of transformer substation engineering project
CN112085065A (en) Low-voltage user-to-home phase identification method based on voltage and active power reading
Yi et al. Improving correlation-based consumer phase identification for incomplete data
CN111245095A (en) Topology identification method of low-voltage distribution network topology identification system
CN113806899B (en) Power distribution network topological relation identification method and device and mobile terminal
CN116819188A (en) Low-voltage distribution network phase sequence identification method based on Bayesian probability theory
CN116257973A (en) Particle swarm optimization-based low-voltage power network line impedance and loss calculation method and system
CN115796711A (en) Auxiliary decision-making method and device for loss reduction of distribution line
CN112000917B (en) Low-voltage user meter box identification method integrating known phase difference and address information
CN112615375A (en) Low-voltage transformer area multi-level branch line user identification method integrated with intelligent switch
CN113075468A (en) Phase determination method based on low-voltage user

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