10kv busbar voltage sensitivity calculation method based on voltage historical data

Abstract

The present application relates to the field of power system automation and control, specifically a 10kV busbar voltage sensitivity calculation method based on voltage historical data, including: reading the distribution network equipment model, reading the mapping relationship between the distribution network feeder and the 10kV busbar, analyzing the distribution network area and the 10kV busbar affiliation, and through the D5000 platform commercial library interface, collecting the historical voltage data of each area and the 96-point historical voltage data of the 10kV busbar, calculating the original sensitivity of the 10kV busbar voltage to the area voltage, and checking the original sensitivity, using the linear regression analysis method to extract the sensitivity of the busbar voltage change to the area voltage change; the sensitivity calculation method can accurately quantify the influence of the busbar voltage on the voltage of each area, thereby realizing dynamic adjustment based on data driving, and significantly improving the efficiency and response speed of voltage regulation.

Description

TECHNICAL FIELD

[0001] The application relates to the field of power system automation and control technology, in particular to a 10kV busbar-to-voltage sensitivity calculation method based on voltage historical data. BACKGROUND

[0002] Voltage sensitivity is a key parameter in circuit design, which can determine the stability of voltage in the running process. For the new power system with a large number of distributed new energy sources connected to the distribution network, the voltage regulation of the distribution network usually depends on fixed regulation strategies and empirical rules, and lacks real-time analysis of the accurate relationship between the bus voltage and the voltage of the distribution area, so the influence of the bus voltage on the voltage of each distribution area cannot be accurately quantified. SUMMARY

[0003] The application provides a 10kV busbar-to-voltage sensitivity calculation method based on voltage historical data, which can accurately calculate the qualified rate of the voltage of each distribution area under the busbar based on the real-time topological structure of the transmission network and the distribution network, and accurately quantify the influence of the bus voltage on the voltage of each distribution area through the sensitivity calculation and sensitivity analysis method, thereby solving the problems in the background art.

[0004] The technical scheme of the application is as follows:

[0005] The 10kV busbar-to-voltage sensitivity calculation method based on voltage historical data comprises the following steps:

[0006] S1: reading the distribution network equipment model and analyzing the topological relationship between the distribution network area and the feeder;

[0007] S2: reading the mapping relationship between the distribution network feeder and the 10kV busbar, and analyzing the subordinate relationship between the distribution network area and the 10kV busbar;

[0008] S3: collecting historical voltage data of each distribution area and 96-point historical voltage data of the 10kV busbar through the D5000 platform commercial library interface; the 96-point refers to the power load data measured every 15 minutes in a day, and there are 96 data points;

[0009] S4: analyzing the distribution network area voltage historical data, performing data verification on the historical data, and generating effective voltage historical data of the distribution area;

[0010] S5: analyzing the 10kV busbar and the voltage historical data of the distribution area, calculating the original sensitivity of the 10kV busbar voltage to the voltage of the distribution area, and checking the original sensitivity;

[0011] S6: using a linear regression analysis method to extract the sensitivity of the bus voltage change to the voltage change of the distribution area.

[0012] Further, in S1, the topology relationship analysis includes the following steps:

[0013] S1.1: Read the distribution network model file through the interface of the distribution automation platform (a commonly used distribution automation system in the current power system);

[0014] S1.2: Parse the distribution network model file, and generate the topology relationship between the distribution network area and the feeder according to the device node topology in the model.

[0015] Further, in S2, the distribution network area and 10kV busbar affiliation analysis process is as follows:

[0016] S2.1: Read the transmission network model file through the D5000 platform interface; D5000 platform is an important technical support platform for power grid operation control and dispatching production management, which integrates the traditional basic application systems of the dispatching center and consists of four key levels: hardware layer, operation program layer, support system layer, and application layer;

[0017] S2.2: Parse the transmission network model file, generate the topology relationship between the 10kV busbar and the feeder according to the device node topology in the model, and generate the distribution network area and 10kV busbar affiliation according to the topology relationship between the distribution network area and the feeder.

[0018] Further, in S4, the specific generation process of the effective voltage historical data of the area is as follows:

[0019] S4.1: According to the distribution network area and 10kV busbar affiliation, parse the historical voltage data of the distribution network area under the 10kV busbar;

[0020] S4.2: According to the upper limit of the distribution network area voltage 10.7kV and the lower limit of the voltage 9.7kV, calculate whether the area voltage is in a reasonable range, and count the voltage distribution of the area under the busbar;

[0021] S4.3: Exclude the area voltage outside the reasonable range to generate the effective voltage historical data of the area.

[0022] Further, in S5, the original sensitivity checking step is as follows:

[0023] S5.1: According to the 10kV busbar historical voltage data and the historical voltage data of the distribution network area under the busbar, calculate the voltage change amount of the 96-point adjacent time historical data of the busbar and the voltage change amount of the 96-point adjacent time historical data of the distribution network area in turn;

[0024] S5.2: Set a 96-point sensitivity list, set the ratio of the voltage change amount of the distribution network area to the busbar voltage at zero time to 1.0, and store it to the 96-point sensitivity list;

[0025] S5.3: Based on the voltage change of the distribution network area and the voltage change of the bus, calculate 95 voltage sensitivities in sequence and store them in a 96-point sensitivity list;

[0026] S5.4: Based on the physical characteristics of the power system, the sensitivity of the 10kV bus voltage b change to the voltage of the distribution network area should be within a reasonable range of 0.8 to 1.2. The 96 points of sensitivity that are not within the reasonable range due to the influence of distributed power sources and load fluctuations in the distribution network are eliminated, and the eliminated points are assigned a default value of 1.0, so as to obtain the 96 points of sensitivity of the 10kV bus voltage change to the voltage of the distribution network area at different times.

[0027] Furthermore, in S6, the sensitivity steps for bus voltage changes to transformer area voltage changes are as follows:

[0028] S6.1: Calculate the average voltage avgLoadVolt of the transformer area based on the historical effective voltage data of the transformer area;

[0029] S6.2: Calculate the average voltage sensitivity avgSens of the transformer area based on the sensitivity of 96 points in the transformer area;

[0030] S6.3: Calculate the covariance using the least squares method In the formula, and Let represent the voltage sensitivity and voltage of the i-th transformer area, respectively;

[0031] S6.4: Calculate the variance using the least squares method ;

[0032] S6.5: Calculate the slope of the voltage sensitivity of the transformer substation to the voltage substation using the least squares method. , cut off ;

[0033] S6.6: Complete the sensitivity extraction of bus voltage change to transformer area voltage change.

[0034] Compared with the prior art, the present invention has the following beneficial effects:

[0035] 1. The application is based on the real-time topology of the power transmission network and the distribution network, accurately calculates the qualified rate of the voltage of each substation under the bus, and uses a sensitivity calculation method to accurately quantify the influence of the bus voltage on the voltage of each substation, thereby realizing dynamic adjustment based on data driving, significantly improving the efficiency and response speed of voltage adjustment; sensitivity analysis can monitor the influence of bus voltage on substation voltage in real time, accurately identify substation with greater influence, provide data support, so that the voltage qualified rate can be improved as a whole, and voltage abnormalities can be avoided. The method aims to accurately and effectively adjust the substation voltage of the distribution network, thereby accurately and efficiently improving the overall level of the substation voltage qualified rate of the distribution network, and ultimately promoting the stable and reliable operation of the power system. BRIEF DESCRIPTION OF DRAWINGS

[0036] Fig. 1 is the flow chart of the voltage sensitivity calculation method of the application;

[0037] Fig. 2 is the flow chart of the linear regression analysis of the application. DETAILED DESCRIPTION

[0038] The embodiments of the application will be further described in detail below in combination with the drawings and examples. The following examples are used to illustrate the application, but cannot be used to limit the scope of the application.

[0039] The application provides a 10kV bus to substation voltage sensitivity calculation method based on voltage historical data, and the following Figs. 1-2 , a detailed description of an embodiment of the method of the application is provided, including the following steps:

[0040] S1: read the distribution network equipment model, analyze the topological relationship between the distribution network substation and the feeder;

[0041] S1.1: read the distribution network model file through the interface of the distribution automation platform (which is the commonly used distribution automation system in the current power system);

[0042] S1.2: analyze the topological relationship between the distribution network substation and the feeder according to the device node topology in the model.

[0043] S2: read the mapping relationship between the distribution network feeder and the 10kV bus, and analyze the subordination relationship between the distribution network substation and the 10kV bus;

[0044] S2.1: read the transmission network model file through the D5000 platform interface; the D5000 platform is an important technical support platform for power grid operation control and dispatching production management, which integrates the traditional basic application systems of the dispatching center and is composed of four key levels: hardware layer, operation program layer, support system layer and application layer;

[0045] S2.2: Analyze the power grid model file, generate the topology relationship of 10kV bus and feeder according to the device node topology in the model, and generate the relationship between the distribution network area and the 10kV bus according to the topology relationship between the distribution network area and the feeder.

[0046] S3: Collect and read the historical voltage data of each area and the 96-point historical voltage data of the 10kV bus through the D5000 platform commercial library interface; 96 points refer to the power load data measured every 15 minutes in a day, a total of 96 data points;

[0047] S4: Analyze the historical voltage data of the distribution network area, and perform data verification on the historical data to generate effective voltage historical data of the area;

[0048] S4.1: According to the relationship between the distribution network area and the 10kV bus, analyze the historical voltage data of the distribution network area under the 10kV bus;

[0049] S4.1.1: In the analysis process, first preprocess the historical voltage data to remove outliers and ensure the accuracy and integrity of the data;

[0050] S4.1.2: The determination process of the outliers is as follows:

[0051] Calculate the mean:

[0052]

[0053] Calculate the standard deviation:

[0054]

[0055] Then calculate the Z-score (Z = (x - x) / s), when |Z| > 3, the data point is an outlier; In the formula, x represents the historical voltage data value, j represents the index of the historical voltage data value, and n represents the number of historical voltage data values;

[0056] S4.2: According to the upper limit of the distribution network area voltage 10.7kV and the lower limit of the voltage 9.7kV, calculate whether the area voltage is in a reasonable range, and count the voltage distribution of the area under the bus;

[0057] S4.3: Remove the area voltage outside the reasonable range to generate effective voltage historical data of the area.

[0058] S5: Analyze the 10kV bus and area voltage historical data, calculate the original sensitivity of the 10kV bus voltage to the area voltage, and check the original sensitivity;

[0059] S5.1: According to the 10kV bus historical voltage data and the historical voltage data of the distribution network area under the bus;​​

[0060] S5.2: Calculate the voltage variation of the adjacent time point of the bus 96 points in sequence;

[0061] S5.2.1: For each pair of adjacent time points, calculate the voltage variation. The specific calculation formula is:

[0062]

[0063] Wherein, is the voltage value of the kth time point, is the voltage value of the k+1th time point;

[0064] S5.3: Calculate the voltage variation of the adjacent time point of the distribution network area 96 points in sequence, and the calculation process is the same as above;

[0065] S5.4: Set a 96 point sensitivity list, set the ratio of the voltage variation of the distribution network area to the bus voltage at zero time point to 1.0, and store it to the 96 point sensitivity list;

[0066] S5.5: According to the voltage variation of the distribution network area and the voltage variation of the bus, calculate the voltage sensitivity of 95 points in sequence, and store it to the 96 point sensitivity list;

[0067] S5.6: According to the physical characteristics of power system, the sensitivity of 10kV bus voltage b change to distribution network area voltage should be in the reasonable range of 0.8 to 1.2, and the 96 point sensitivity affected by the fluctuation of distributed power and load outside the reasonable range is removed, and the removed point is assigned a default value of 1.0, and the 96 point sensitivity of 10kV bus voltage change to distribution network area voltage at different time is obtained.

[0068] S6: Using linear regression analysis method, extract the sensitivity of bus voltage change to area voltage change;

[0069] S6.1: According to the effective voltage history data of the area, calculate the average voltage avgLoadVolt of the area;

[0070] S6.2: According to the 96 point sensitivity of the area, calculate the average voltage sensitivity avgSens of the area;

[0071] S6.3: Calculate the covariance by least square method ;

[0072] S6.4: Calculate the variance by least square method ;

[0073] S6.5: Calculate the slope of the area voltage to the area voltage sensitivity by least square method , intercept ;

[0074] S6.6: Complete the sensitivity extraction of bus voltage change to transformer area voltage change.

[0075] The embodiments of the present invention are given for the purposes of illustration and description. Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for calculating the voltage sensitivity of a 10kV busbar to a transformer substation based on historical voltage data, characterized in that, Includes the following steps: S1: Read the power distribution network equipment model and analyze the topological relationship between the power distribution network area and the feeder; S2: Read the mapping relationship between the distribution network feeder and the 10kV bus, and analyze the subordinate relationship between the distribution network area and the 10kV bus; S3: Collect historical voltage data of each transformer substation and 96 historical voltage data of the 10kV busbar through the commercial library interface of the D5000 platform. S4: Analyze the historical voltage data of the distribution network area, verify the historical data, and generate the effective voltage historical data of the area. S5: Analyze historical data of 10kV bus and transformer substation voltage, calculate the original sensitivity of 10kV bus voltage to transformer substation voltage, and verify the original sensitivity. S6: Linear regression analysis is used to extract the sensitivity of bus voltage changes to transformer area voltage changes; In S5, the original sensitivity verification steps are as follows: S5.1: Based on the historical voltage data of the 10kV bus and the historical voltage data of the distribution network area below the bus, calculate the voltage change of the historical data of the bus at 96 adjacent times and the voltage change of the distribution network area at 96 adjacent times in sequence. S5.2: Set a 96-point sensitivity list, set the ratio of the change in voltage of the distribution network area to the change in bus voltage at zero point to 1.0, and store it in the 96-point sensitivity list; S5.3: Based on the voltage change of the distribution network area and the voltage change of the bus, calculate 95 voltage sensitivities in sequence and store them in a 96-point sensitivity list; S5.4: Based on the physical characteristics of the power system, the sensitivity of the 10kV bus voltage b change to the voltage of the distribution network area should be within a reasonable range of 0.8 to 1.

2. The 96 points of sensitivity that are not within the reasonable range due to the influence of distributed power sources and load fluctuations in the distribution network are eliminated, and the eliminated points are assigned a default value of 1.

0. The 96 points of sensitivity of the 10kV bus voltage change to the voltage of the distribution network area at different times are obtained. In S6, the sensitivity steps of bus voltage change to transformer area voltage change are as follows: S6.1: Calculate the average voltage avgLoadVolt of the transformer area based on the historical effective voltage data of the transformer area; S6.2: Calculate the average voltage sensitivity avgSens of the transformer area based on the sensitivity of 96 points in the transformer area; S6.3: Calculate the covariance using the least squares method ; S6.4: Calculate the variance using the least squares method ; In the formula, and Let represent the voltage sensitivity and voltage of the i-th transformer area, respectively; S6.5: Calculate the slope of the voltage sensitivity of the transformer substation to the voltage substation using the least squares method. , cut off ; S6.6: Complete the sensitivity extraction of bus voltage change to transformer area voltage change.

2. The method for calculating the voltage sensitivity of a 10kV busbar to a transformer substation based on historical voltage data as described in claim 1, characterized in that: In S1, topological relationship analysis includes the following steps: S1.1: Read the distribution network model file through the distribution automation platform interface; S1.2: Parse the distribution network model file and generate the topological relationship between the distribution network area and the feeder based on the topology of the equipment nodes in the model.

3. The method for calculating the voltage sensitivity of a 10kV busbar to a transformer substation based on historical voltage data as described in claim 2, characterized in that: In S2, the analysis process of the subordinate relationship between the distribution network area and the 10kV bus is as follows: S2.1: Read the transmission network model file through the D5000 platform interface; S2.2: Parse the transmission network model file, generate the topological relationship between the 10kV bus and the feeder based on the equipment node topology in the model, and generate the subordinate relationship between the distribution network area and the 10kV bus based on the topological relationship between the distribution network area and the feeder.

4. The method for calculating the voltage sensitivity of a 10kV busbar to a transformer substation based on historical voltage data as described in claim 3, characterized in that: In S4, the specific process for generating historical effective voltage data for transformer substations is as follows: S4.1: Based on the subordinate relationship between the distribution network area and the 10kV bus, analyze the historical voltage data of the distribution network area under the 10kV bus; S4.2: Based on the upper limit of voltage of 10.7kV and the lower limit of voltage of distribution network area, calculate whether the voltage of distribution area is within a reasonable range, and count the voltage distribution of distribution area under bus. S4.3: Remove voltages from transformer substations that are outside the reasonable range and generate historical data on the effective voltage of the substations.

Images