A harmonic responsibility division method and device, electronic equipment and storage medium

Abstract

The application discloses a harmonic responsibility division method and device, electronic equipment and a storage medium, and is used to solve the technical problem that the existing harmonic responsibility division method is dependent on physical modeling, is limited by data acquisition, is difficult to apply in actual engineering, and leads to poor harmonic responsibility division effect in actual application. The application comprises the following steps: obtaining harmonic voltage at a PCC point of a power distribution network and average active power of each user; constructing a first time sequence of the harmonic voltage and a second time sequence of the average active power; constructing a first time sequence feature matrix of the harmonic voltage by using the first time sequence; constructing a second time sequence feature matrix of the average active power of each user by using the second time sequence; calculating a correlation index of the first time sequence feature matrix and each second time sequence feature matrix; and calculating a harmonic responsibility proportion of each user according to the correlation index.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of harmonic analysis, and in particular to a harmonic responsibility division method and device, electronic equipment and a storage medium. BACKGROUND

[0002] With the wide access of distributed new energy, the disturbance uncertainty of distribution network increases, the coupling characteristics of multi-harmonic sources are enhanced, the harmonic pollution is intensified and increasingly complex, which brings severe challenges to high-quality power supply. Therefore, the power department needs to continuously strengthen the research on responsibility division under the multi-harmonic source scene, reasonably assess the harmonic responsibility, and continuously improve the supervision level of harmonics.

[0003] At present, most of the harmonic responsibility division methods focus on calculating accurate equivalent harmonic impedance based on the physical model of the power grid as the theoretical basis, and then establishing a responsibility evaluation index based on the vector projection of harmonic voltage. However, these methods need to include low-latency synchronous vector time-domain recording data with harmonic phase. However, limited by the existing communication architecture and terminal device computing power of the distribution network, the power quality monitoring devices in the actual distribution network can only collect and analyze part of the power quality steady-state data of the distribution network, making it difficult for the above-mentioned model-based methods to be applied in actual engineering, resulting in poor harmonic responsibility division effect in actual application. SUMMARY

[0004] The present application provides a harmonic responsibility division method, device, electronic equipment and storage medium, which is used to solve the technical problem that the existing harmonic responsibility division method relies on physical modeling and is limited by data collection, making it difficult to be applied in actual engineering, resulting in poor harmonic responsibility division effect in actual application.

[0005] The present application provides a harmonic responsibility division method, comprising:

[0006] Obtaining the harmonic voltage at the PCC point of the distribution network and the average active power of each user;

[0007] Constructing a first time series of the harmonic voltage and a second time series of the average active power;

[0008] Constructing a first time series feature matrix of the harmonic voltage using the first time series;

[0009] Constructing a second time series feature matrix of the average active power of each user using the second time series;

[0010] Calculating the correlation index of the first time series feature matrix and each second time series feature matrix;

[0011] According to the correlation index, calculating the harmonic responsibility proportion of each user.

[0012] Optionally, the step of constructing the first time sequence feature matrix of the harmonic voltage by using the first time sequence comprises:

[0013] dividing the first time sequence into a plurality of first time periods;

[0014] extracting first feature parameters from the first time periods;

[0015] constructing the first time sequence feature matrix of the harmonic voltage by using the first feature parameters.

[0016] Optionally, the step of dividing the first time sequence into a plurality of first time periods comprises:

[0017] determining trend turning points of the first time sequence;

[0018] dividing the first time sequence into a plurality of first time periods based on the trend turning points.

[0019] Optionally, the first feature parameters comprise fitting slopes and time spans; the step of extracting first feature parameters from the first time periods comprises:

[0020] calculating fitting slopes of the first time periods;

[0021] obtaining the number of monitoring data points on the first time periods and the current time of each monitoring data point;

[0022] calculating the time span of the first time period according to the current time of each monitoring data point.

[0023] Optionally, the step of constructing the first time sequence feature matrix of the harmonic voltage by using the first feature parameters comprises:

[0024] generating the inclination angle of each first time period by using the fitting slope;

[0025] calculating the time span ratio of each first time period by using the time span and the length of the first time sequence;

[0026] generating the first time sequence feature matrix of the harmonic voltage by using the inclination angle and the time span ratio of each first time period.

[0027] Optionally, the step of calculating the harmonic responsibility proportion of each user according to the correlation index comprises:

[0028] calculating the total average active power of all the users according to the average active power of each user;

[0029] The average active power, the total average active power and the correlation index are used to calculate a harmonic responsibility index of each user;

[0030] The harmonic responsibility index is normalized to obtain a harmonic responsibility proportion of each user.

[0031] The application further provides a harmonic responsibility division device, comprising:

[0032] A harmonic voltage and average active power acquisition module is configured to acquire harmonic voltage at a PCC point of a power distribution network and average active power of each user;

[0033] A time sequence construction module is configured to construct a first time sequence of the harmonic voltage and a second time sequence of the average active power;

[0034] A first time sequence feature matrix construction module is configured to construct a first time sequence feature matrix of the harmonic voltage by using the first time sequence;

[0035] A second time sequence feature matrix construction module is configured to construct a second time sequence feature matrix of the average active power of each user by using the second time sequence;

[0036] A correlation index calculation module is configured to calculate a correlation index of the first time sequence feature matrix and each second time sequence feature matrix;

[0037] A harmonic responsibility proportion calculation module is configured to calculate a harmonic responsibility proportion of each user according to the correlation index.

[0038] Optionally, the first time sequence feature matrix construction module comprises:

[0039] A first time period division sub-module is configured to divide the first time sequence into a plurality of first time periods;

[0040] A first feature parameter extraction sub-module is configured to extract a first feature parameter from the first time period;

[0041] A first time sequence feature matrix construction sub-module is configured to construct the first time sequence feature matrix of the harmonic voltage by using the first feature parameter.

[0042] The application further provides an electronic device, which comprises a processor and a memory:

[0043] The memory is configured to store program code and transmit the program code to the processor;

[0044] The processor is configured to execute the harmonic responsibility division method according to the instructions in the program code.

[0045] The application further provides a computer readable storage medium for storing program code for performing the harmonic responsibility division method according to any one of the above.

[0046] It can be seen from the above technical solutions that the application has the following advantages: the application discloses a harmonic responsibility division method, which comprises the following steps: obtaining harmonic voltage at a PCC point of a power distribution network and average active power of each user; constructing a first time sequence of the harmonic voltage and a second time sequence of the average active power; constructing a first time sequence feature matrix of the harmonic voltage by using the first time sequence; constructing a second time sequence feature matrix of the average active power of each user by using the second time sequence; calculating a correlation index of the first time sequence feature matrix and each second time sequence feature matrix; and calculating a harmonic responsibility proportion of each user according to the correlation index.

[0047] The application obtains the harmonic voltage and the average active power of the user from the PCC point and generates the corresponding time sequence feature matrix, so that the harmonic responsibility proportion of the user is calculated by analyzing the correlation between the time sequence matrices, without physical modeling and without collecting low-time-delay synchronous phasor time-domain recording data of the harmonic phase, thereby reducing the data collection pressure and quickly and accurately realizing the harmonic responsibility division. BRIEF DESCRIPTION OF DRAWINGS

[0048] In order to more clearly illustrate the technical solutions in the embodiments of the application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description only constitute some embodiments of the application, and for those skilled in the art, other drawings can also be obtained from these drawings without creative labor.

[0049] Figure 1 A step flow chart of a harmonic responsibility division method provided by an embodiment of the application;

[0050] Figure 2 A schematic diagram of a multi-harmonic source equivalent circuit provided by an embodiment of the application;

[0051] Figure 3 A step flow chart of a harmonic responsibility division method provided by another embodiment of the application;

[0052] Figure 4 A structural block diagram of a harmonic responsibility division device provided by an embodiment of the application. DETAILED DESCRIPTION

[0053] The embodiment of the present application provides a harmonic responsibility division method, device, electronic equipment and storage medium, and is used for solving the technical problem that the existing harmonic responsibility division method is dependent on physical modeling, is limited by the limitation of data acquisition, is difficult to be applied in actual engineering, and results in poor harmonic responsibility division effect in actual application.

[0054] In order to make the application purpose, features and advantages of the present application more obvious and easy to understand, the technical solutions in the embodiments of the present application will be clearly and completely described below in combination with the drawings in the embodiments of the present application. Obviously, the following described embodiments are only a part of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the protection scope of the present application.

[0055] Please refer to Figure 1 , Figure 1 The embodiment of the present application provides a harmonic responsibility division method.

[0056] The harmonic responsibility division method provided by the present application can specifically include the following steps:

[0057] Step 101, obtaining the harmonic voltage at the PCC point of the power distribution network and the average active power of each user;

[0058] In actual application, the similarity of the harmonic source user active power and the harmonic voltage change rule is high, and the non-harmonic user data has no obvious similarity, so in the embodiment of the present application, the harmonic responsibility division of the user can be performed based on the harmonic voltage and the active power of the user.

[0059] In the embodiment of the present application, the multi-harmonic source equivalent circuit is as shown in Figure 2 The steady-state harmonic voltage V of the PCC point of the power distribution network can be obtained by the power quality monitoring device, and the average active power P of each user at the PCC point can be obtained by the power consumption monitoring device.

[0060] Step 102, constructing a first time sequence of the harmonic voltage and a second time sequence of the average active power;

[0061] After the harmonic voltage of the PCC point is collected, the first time sequence of the harmonic voltage and the second time sequence of the average active power of the user can be constructed respectively.

[0062] The first time sequence represents the change of the harmonic voltage with time; and the second time sequence represents the change of the average active power of the user with time.

[0063] Step 103, constructing a first time sequence feature matrix of the harmonic voltage by using the first time sequence;

[0064] After the first time sequence is acquired, a characteristic parameter of the first time sequence can be extracted to construct a first time sequence characteristic matrix of the harmonic voltage.

[0065] In step 104, a second time sequence characteristic matrix of the average active power of each user is constructed by using the second time sequence.

[0066] Similarly, after the second time sequence is acquired, a time parameter of the second time sequence can be extracted to construct a second time sequence characteristic matrix of the average active power of each user.

[0067] In step 105, a correlation index of the first time sequence characteristic matrix and each second time sequence characteristic matrix is calculated.

[0068] In the embodiment of the present application, the average active power of different users is different, which can cause a difference in the similarity of the harmonic voltage curve. Therefore, the similarity of the harmonic voltage curve and the average active power curve of the user can represent the interference degree of the user to the harmonic voltage. Therefore, in the embodiment of the present application, the correlation index of the first time sequence characteristic matrix and each second time sequence characteristic matrix can be calculated to analyze the harmonic responsibility proportion of each user through the correlation index.

[0069] In step 106, the harmonic responsibility proportion of each user is calculated according to the correlation index.

[0070] In the embodiment of the present application, after the correlation index of the average active power of each user and the harmonic voltage is calculated, the harmonic responsibility proportion of each user can be calculated according to the correlation index.

[0071] The present application obtains the harmonic voltage and the average active power of the user from the PCC point, and generates a corresponding time sequence characteristic matrix, so as to calculate the harmonic responsibility proportion of the user by analyzing the correlation between the time sequence matrices, without physical modeling and without collecting low-time-delay synchronous phasor time-domain recording data of the harmonic phase, thereby reducing the data acquisition pressure and quickly and accurately realizing the harmonic responsibility division.

[0072] Please refer to Figure 3 , Figure 3 A step flowchart of a harmonic responsibility division method provided by another embodiment of the present application can include the following steps.

[0073] In step 301, the harmonic voltage at the PCC point of the power distribution network and the average active power of each user are acquired.

[0074] In step 302, a first time sequence of the harmonic voltage and a second time sequence of the average active power are constructed.

[0075] Steps 301-302 are the same as steps 101-102, and details can be referred to the description of steps 101-102, which will not be repeated here.

[0076] Step 303, dividing the first time sequence into a plurality of first time periods;

[0077] In the embodiments of the present application, in order to capture the harmonic data change rule of different time intervals in the first time sequence, the first time sequence can be divided into a plurality of first time periods according to the change rule of the harmonic data.

[0078] In one example, the step of dividing the first time sequence into a plurality of first time periods can specifically include the following sub-steps:

[0079] S31, determining a trend turning point of the first time sequence;

[0080] S32, dividing the first time sequence into a plurality of first time periods based on the trend turning point.

[0081] In a specific implementation, the change trend of the first time sequence can be divided into three categories: rising, stable, and falling, so the first time sequence can be divided into several time periods by finding the trend turning point of the time sequence.

[0082] The trend turning point of the first time sequence is determined according to the following formula: The trend turning point of the first time sequence can be calculated by the following formula:

[0083]

[0084] Or,

[0085]

[0086] Wherein, x i-1 , x i+1 are power quality monitoring data (here referring to harmonic voltage) of previous and subsequent time.

[0087] Step 304, extracting a first feature parameter from the first time period;

[0088] After the segmentation of the first time sequence is completed, the first feature parameter can be extracted from the first time period.

[0089] In one example, the first feature parameter includes a fitting slope and a time span; the step of extracting the first feature parameter from the first time period can include the following sub-steps:

[0090] S41, calculating a fitting slope of the first time period;

[0091] ​S42, obtain the number of monitoring data points on the first time period, and the current time of each monitoring data point;

[0092] S43, calculate the time span of the first time period according to the current time of each monitoring data point.

[0093] In a specific implementation, after completing the segmentation of the first time sequence, least square fitting can be used for each segmented first time period to perform segmented linear representation, representing the main local trend characteristics and corresponding time characteristics of the harmonic voltage and average active power in the time scale. The expression of the segmented linear representation of the first time period is:

[0094] y n (t)=k n t+b n ,n=1,2,…,N

[0095] In the formula, k n is the fitting slope of the data in the nth first time period, representing the main local trend characteristics in the first time period, and k n is larger, the more intense the data trend of the first time period is, and the definition interval of k n is (-∞, +∞); b n is the fitting intercept of the data in the nth first time period; and N is the number of segments of the first time sequence.

[0096] The time span s n is the time characteristic of the first time sequence, and is expressed as:

[0097] s n =t l -t1

[0098] Wherein, t i (i∈[1, l]) is the time of the monitoring data point x i on the first time period; and l is the number of monitoring data points in the nth first time period. The larger the time span s n , the stronger the description of the first time period on the time sequence pattern is, and the definition interval of s n is (0, +∞).

[0099] Step 305, constructing a first time sequence characteristic matrix of the harmonic voltage by using the first characteristic parameter;

[0100] After obtaining the first characteristic parameter, a first time sequence characteristic matrix of the harmonic voltage can be constructed by using the first characteristic parameter.

[0101] In an example, the step of constructing a first time sequence characteristic matrix of the harmonic voltage by using the first characteristic parameter can include the following sub-steps:

[0102] S51, generating the inclination angle of each first time period by using the fitted slope;

[0103] S52, calculating the time span ratio of each first time period by using the time span and the length of the first time sequence;

[0104] S53, generating the first time sequence feature matrix of the harmonic voltage by using the inclination angle and the time span ratio of each first time period.

[0105] In a specific implementation, since the fitted slope k n , the time span s n , the value range of which is infinite, in order to avoid the time sequence feature being close, and k n , s n difference is too large to produce misjudgment, k n is converted into the inclination angle a n , s n is converted into the time span ratio r n relative to the length S of the first time sequence, that is:

[0106]

[0107]

[0108] After conversion, the harmonic voltage of the PCC point is represented as the first time sequence feature matrix, that is:

[0109] [(a1, r1) (a2, r2) … (a N , r N )]

[0110] Step 306, constructing the second time sequence feature matrix of the average active power of each user by using the second time sequence;

[0111] Similarly, in the embodiment of the application, the second time sequence feature matrix of the average active power of each user can be constructed by using the second time sequence according to the specific description of step 305. Here, no longer be described.

[0112] Step 307, calculating the correlation index of the first time sequence feature matrix and each second time sequence feature matrix;

[0113] In the embodiment of the application, since the time sequence local features of the harmonic voltage and the average active power of each user are different, the key trend turning point The number of segments N obtained by dividing the time series is likely to be different, resulting in a difference in the number of columns of the time series feature matrix. The T-TD based on DTW supports pattern matching of unequal length time series and can prevent the pathological bending phenomenon that may occur when DTW distance matches. The method sets the bending path limit for calculating the minimum cumulative feature distance with the tilt angle α n , the time span ratio r n as a characteristic parameter, implements pattern matching on the time series, and measures the correlation between the PCC point harmonic voltage and the average active power time series feature matrix of different users.

[0114] Specifically, let the first time series feature matrix of the PCC point harmonic voltage be V, and the second time series feature matrix of the average active power of a certain user be P, i.e.:

[0115] V=[(α1,r1)(α2,r2)…(α N ,r N )]=[v1,v2,…,v N ]

[0116] P=[(α1,r1)(α2,r2)…(α N' ,r N' )]=[p1,p2,…,p N' ]

[0117] Wherein, the T-TD definition of the PCC point harmonic voltage V and the average active power P of a certain user is:

[0118]

[0119]

[0120] Wherein, X[1:N] represents a subsequence composed of the first 1 to N column vectors of the time series feature matrix X; represents an empty sequence; represents the base feature distance between v u and p v , representing the difference between the main local trend characteristics and the time series of the u-th time segment of the PCC point harmonic voltage V and the v-th time segment of the average active power P of a certain user. In addition, in order to avoid pathological bending phenomenon, the maximum bending limit is set by the following formula, i.e.:

[0121]

[0122] Wherein, represents the upward rounding operation; δ·N represents the maximum bending radius, and the parameter δ can generally be taken as 0.2.

[0123] Step 308, calculating the harmonic responsibility proportion of each user according to the correlation index.

[0124] In the embodiment of the present application, after calculating the correlation index of the average active power and the harmonic voltage of each user, the harmonic responsibility proportion of each user can be calculated according to the correlation index.

[0125] In one example, the step of calculating the harmonic responsibility proportion of each user according to the correlation index can include the following sub-steps:

[0126] S81, calculating the total average active power of all users according to the average active power of each user;

[0127] S82, calculating the harmonic responsibility index of each user by using the average active power, the total average active power and the correlation index;

[0128] S83, normalizing the harmonic responsibility index to obtain the harmonic responsibility proportion of each user.

[0129] In the specific implementation, the smaller the T-TD of the PCC point harmonic voltage V and the average active power P of a certain user, the more similar the change rules of the two time series, the stronger the correlation, which indicates that the harmonic pollution emitted by the user has a greater impact on the PCC point harmonic voltage and the harmonic responsibility is greater.

[0130] Therefore, according to this rule, the correlation index C k is constructed as follows:

[0131]

[0132] C k The greater the C , the greater the harmonic pollution caused by the user to the PCC point.

[0133] Considering the relationship between the average active power of the user and the harmonic responsibility of the user, a sequence Z k is constructed, and the element Z k of the sequence is calculated as follows:

[0134]

[0135] In the formula, C k represents the correlation index of the user k; P k is the average active power value corresponding to the user k and the correlation coefficient, M is the total sum of the analyzed users, and the sum of the average active power of all users is represented by ΣP.

[0136] In order to intuitively compare the harmonic responsibility of each user, Z kThe normalization is performed so that the harmonic responsibility of each user is within the interval [0, 1], and the sum of the harmonic responsibilities of all users is 1. The normalization process is shown as follows:

[0137]

[0138]

[0139] wherein Z max and Z min are the maximum and minimum values of the harmonic responsibility index Z k of all users, respectively; H k represents the amount of negative values removed from Z k ; and F k represents the normalized harmonic responsibility index of user k, indicating the proportion of the harmonic responsibility of user k in the total harmonic responsibility of all analyzed users M, i.e., the harmonic responsibility proportion of the user. Thus, the harmonic responsibility division is completed.

[0140] The present application obtains the harmonic voltage and the average active power of users from the PCC point, and generates the corresponding time sequence feature matrix, so as to calculate the harmonic responsibility proportion of the user by analyzing the correlation between the time sequence matrices, without physical modeling and without collecting low-time-delay synchronous phasor time-domain recording data of the harmonic phase, thereby reducing the data collection pressure and quickly and accurately realizing the harmonic responsibility division.

[0141] Referring to Figure 4 , Figure 4 FIG. 1 is a structural block diagram of a harmonic responsibility division device provided by an embodiment of the present application.

[0142] The embodiment of the present application provides a harmonic responsibility division device, which comprises:

[0143] A harmonic voltage and average active power acquisition module 401 is configured to acquire the harmonic voltage at the PCC point of a power distribution network and the average active power of each user.

[0144] A time sequence construction module 402 is configured to construct a first time sequence of the harmonic voltage and a second time sequence of the average active power.

[0145] A first time sequence feature matrix construction module 403 is configured to construct a first time sequence feature matrix of the harmonic voltage by using the first time sequence.

[0146] A second time sequence feature matrix construction module 404 is configured to construct a second time sequence feature matrix of the average active power of each user by using the second time sequence.

[0147] A correlation index calculation module 405 is configured to calculate the correlation index of the first time sequence feature matrix and each second time sequence feature matrix.

[0148] The harmonic responsibility proportion calculation module 406 is configured to calculate the harmonic responsibility proportion of each user according to the correlation index.

[0149] In the embodiment of the present application, the first time sequence feature matrix construction module 403 comprises:

[0150] The first time period division sub-module is configured to divide the first time sequence into a plurality of first time periods.

[0151] The first feature parameter extraction sub-module is configured to extract a first feature parameter from the first time period.

[0152] The first time sequence feature matrix construction sub-module is configured to construct a first time sequence feature matrix of the harmonic voltage by using the first feature parameter.

[0153] In the embodiment of the present application, the first time period division sub-module comprises:

[0154] The trend turning point determination unit is configured to determine a trend turning point of the first time sequence.

[0155] The first time period division unit is configured to divide the first time sequence into a plurality of first time periods based on the trend turning point.

[0156] In the embodiment of the present application, the first feature parameter comprises a fitting slope and a time span; and the first feature parameter extraction sub-module comprises:

[0157] The data fitting intercept acquisition unit is configured to acquire a data fitting intercept of the first time period.

[0158] The fitting slope calculation unit is configured to calculate a fitting slope of the first time period by using the data fitting intercept.

[0159] The monitoring data acquisition unit is configured to acquire the number of monitoring data points on the first time period and the current time point of each monitoring data point.

[0160] The time span calculation unit is configured to calculate the time span of the first time period according to the current time point of each monitoring data point.

[0161] In the embodiment of the present application, the first time sequence feature matrix construction sub-module comprises:

[0162] The inclination angle generation unit is configured to generate an inclination angle of each first time period by using the fitting slope.

[0163] The time span ratio calculation unit is configured to calculate the time span ratio of each first time period by using the time span and the length of the first time sequence.

[0164] The first time sequence feature matrix generating unit is configured to generate a first time sequence feature matrix of harmonic voltage by using the tilt angle and the time span ratio of each first time period.

[0165] In the embodiment of the present application, the harmonic responsibility proportion calculation module 406 comprises:

[0166] The total average active power calculation sub-module is configured to calculate total average active power of all users according to the average active power of each user.

[0167] The harmonic responsibility index calculation sub-module is configured to calculate the harmonic responsibility index of each user by using the average active power, the total average active power and the correlation index.

[0168] The harmonic responsibility proportion calculation sub-module is configured to normalize the harmonic responsibility index to obtain the harmonic responsibility proportion of each user.

[0169] The embodiment of the present application also provides an electronic device, which comprises a processor and a memory:

[0170] The memory is configured to store program code and transmit the program code to the processor.

[0171] The processor is configured to execute the harmonic responsibility division method according to the instructions in the program code.

[0172] The embodiment of the present application also provides a computer readable storage medium, which is configured to store program code, and the program code is configured to execute the harmonic responsibility division method.

[0173] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the system, device and unit described above can refer to the corresponding process in the foregoing method embodiments, and will not be described here.

[0174] Each embodiment in the specification is described in a progressive manner, and each embodiment focuses on the difference from other embodiments, and the same and similar parts between each embodiment can be referred to.

[0175] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, device or computer program product. Therefore, the embodiments of the present application can adopt a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware aspects. Moreover, the embodiments of the present application can adopt the form of a computer program product implemented on one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program code.

[0176] The embodiments of the present application are described with reference to the flowchart illustrations and / or block diagrams of the methods, terminal devices (systems) and computer program products according to the embodiments of the present application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing terminal devices to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal devices, create means for implementing the functions specified in the flowchart illustrations and / or block diagrams. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams.

[0177] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing terminal devices to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart illustrations and / or block diagrams. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams.

[0178] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal devices, such that a series of operational steps are carried out on the computer or other programmable terminal devices to produce a computer implemented process so that the instructions executed on the computer or other programmable terminal devices provide steps for implementing the functions specified in the flowchart illustrations and / or block diagrams. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams.

[0179] Although preferred embodiments of the present application have been described, those skilled in the art will be able to make additional modifications and variations to these embodiments without departing from the scope of the present application. Accordingly, the appended claims are intended to encompass all such modifications and variations as falling within the scope of the present application.

[0180] Finally, it needs to be pointed out that in this document, relational terms such as first and second and the like can only be used to distinguish one entity or action from another entity or action, without necessarily requiring or implying that there is any such actual relationship or order between these entities or actions. Moreover, the terms "comprises", "comprising", or any other variations thereof are intended to cover non-exclusive inclusions, so that a process, method, article, or terminal device including a list of elements does not only include those elements, but also includes other elements not explicitly listed, or further includes elements inherent in such a process, method, article, or terminal device. Without more limitations, the element defined by the statement "comprises a" does not exclude the presence of additional identical elements in the process, method, article, or terminal device including the element.

[0181] The above description and the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit the present application; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacements for part of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method of harmonic responsibility division, characterized by, The method comprises the following steps: obtaining harmonic voltage at a PCC point of a power distribution network and average active power of each user; constructing a first time sequence of the harmonic voltage and a second time sequence of the average active power; constructing a first time sequence feature matrix of the harmonic voltage by using the first time sequence; constructing a second time sequence feature matrix of the average active power of each user by using the second time sequence; calculating a correlation index of the first time sequence feature matrix and each second time sequence feature matrix; calculating harmonic responsibility proportion of each user according to the correlation index; wherein the step of constructing the first time sequence feature matrix of the harmonic voltage by using the first time sequence comprises: dividing the first time sequence into a plurality of first time periods; extracting a first feature parameter from the first time period; constructing the first time sequence feature matrix of the harmonic voltage by using the first feature parameter; wherein the first feature parameter comprises a fitting slope and a time span; the step of extracting the first feature parameter from the first time period comprises: calculating the fitting slope of the first time period; obtaining the number of monitoring data points on the first time period and the current time of each monitoring data point; calculating the time span of the first time period according to the current time of each monitoring data point; wherein the step of constructing the first time sequence feature matrix of the harmonic voltage by using the first feature parameter comprises: generating an inclination angle of each first time period by using the fitting slope; calculating a time span ratio of each first time period by using the time span and the length of the first time sequence; generating the first time sequence feature matrix of the harmonic voltage by using the inclination angle and the time span ratio of each first time period; wherein the step of calculating harmonic responsibility proportion of each user according to the correlation index comprises: calculating total average active power of all users according to the average active power of each user; calculating a harmonic responsibility index of each user by using the average active power, the total average active power and the correlation index; normalizing the harmonic responsibility index to obtain the harmonic responsibility proportion of each user.

2. The method of claim 1, wherein, The step of dividing the first time sequence into a plurality of first time periods comprises: determining a trend turning point of the first time sequence; dividing the first time sequence into a plurality of first time periods based on the trend turning point.

3. A harmonic responsibility dividing device characterized by comprising: The method comprises the following steps: harmonic voltage and average active power obtaining module, used for obtaining harmonic voltage at a PCC point of a power distribution network and average active power of each user; time sequence constructing module, used for constructing a first time sequence of the harmonic voltage and a second time sequence of the average active power; first time sequence feature matrix constructing module, used for constructing a first time sequence feature matrix of the harmonic voltage by using the first time sequence; second time sequence feature matrix constructing module, used for constructing a second time sequence feature matrix of the average active power of each user by using the second time sequence; calculating a correlation index of the first time sequence feature matrix and each second time sequence feature matrix; calculating harmonic responsibility proportion of each user according to the correlation index; wherein the step of constructing the first time sequence feature matrix of the harmonic voltage by using the first time sequence comprises: dividing the first time sequence into a plurality of first time periods; extracting a first feature parameter from the first time period; constructing the first time sequence feature matrix of the harmonic voltage by using the first feature parameter; wherein the first feature parameter comprises a fitting slope and a time span; the step of extracting the first feature parameter from the first time period comprises: calculating the fitting slope of the first time period; obtaining the number of monitoring data points on the first time period and the current time of each monitoring data point; calculating the time span of the first time period according to the current time of each monitoring data point; wherein the step of constructing the first time sequence feature matrix of the harmonic voltage by using the first feature parameter comprises: generating an inclination angle of each first time period by using the fitting slope; calculating a time span ratio of each first time period by using the time span and the length of the first time sequence; generating the first time sequence feature matrix of the harmonic voltage by using the inclination angle and the time span ratio of each first time period; wherein the step of calculating harmonic responsibility proportion of each user according to the correlation index comprises: calculating total average active power of all users according to the average active power of each user; calculating a harmonic responsibility index of each user by using the average active power, the total average active power and the correlation index; normalizing the harmonic responsibility index to obtain the harmonic responsibility proportion of each user. The step of dividing the first time sequence into a plurality of first time periods comprises: determining a trend turning point of the first time sequence; dividing the first time sequence into a plurality of first time periods based on the trend turning point. The method comprises the following steps: harmonic voltage and average active power obtaining module, used for obtaining harmonic voltage at a PCC point of a power distribution network and average active power of each user; time sequence constructing module, used for constructing a first time sequence of the harmonic voltage and a second time sequence of the average active power; first time sequence feature matrix constructing module, used for constructing a first time sequence feature matrix of the harmonic voltage by using the first time sequence; second time sequence feature matrix constructing module, used for constructing a second time sequence feature matrix of the average active power of each user by using the second time sequence; The correlation index calculation module is configured to calculate a correlation index of the first time sequence feature matrix and each second time sequence feature matrix. The harmonic responsibility proportion calculation module is configured to calculate a harmonic responsibility proportion of each user according to the correlation index. The first time sequence feature matrix construction module includes: The first time period division sub-module is configured to divide the first time sequence into a plurality of first time periods. The first feature parameter extraction sub-module is configured to extract a first feature parameter from the first time period. The first time sequence feature matrix construction sub-module is configured to construct the first time sequence feature matrix of the harmonic voltage by using the first feature parameter. The first feature parameter includes a fitting slope and a time span. The data fitting intercept acquisition unit is configured to acquire a data fitting intercept of the first time period. The fitting slope calculation unit is configured to calculate a fitting slope of the first time period by using the data fitting intercept. The monitoring data acquisition unit is configured to acquire a number of monitoring data points in the first time period and a current time point of each monitoring data point. The time span calculation unit is configured to calculate a time span of the first time period according to the current time point of each monitoring data point. The first time sequence feature matrix construction sub-module includes: The inclination angle generation unit is configured to generate an inclination angle of each first time period by using the fitting slope. The time span ratio calculation unit is configured to calculate a time span ratio of each first time period by using the time span and a length of the first time sequence. The first time sequence feature matrix generation unit is configured to generate the first time sequence feature matrix of the harmonic voltage by using the inclination angle and the time span ratio of each first time period. The harmonic responsibility proportion calculation module includes: The total average active power calculation sub-module is configured to calculate a total average active power of all users according to average active powers of the users. The harmonic responsibility index calculation sub-module is configured to calculate a harmonic responsibility index of each user by using the average active power, the total average active power, and the correlation index. The harmonic responsibility proportion calculation sub-module is configured to normalize the harmonic responsibility index to obtain the harmonic responsibility proportion of each user.

4. An electronic device, comprising: The device includes a processor and a memory: The memory is configured to store program code and transmit the program code to the processor. The processor is configured to execute the harmonic responsibility division method according to the instructions in the program code.

5. A computer readable storage medium, characterized in that, The computer readable storage medium is configured to store program code for executing the harmonic responsibility division method. The computer readable storage medium is configured to store program code for executing the harmonic responsibility division method.

Images