High-capacity industrial load safety quantification and peak regulation method based on energy supply safety index

Abstract

The invention discloses a high-capacity industrial load safety quantification and peak regulation method based on an energy supply safety index, and belongs to the technical field of electrical engineering. The method comprises the following steps: firstly, selecting and determining an industrial load interruption priority, formulating an industrial load energy supply safety rating by applying the priority, and determining an industrial load interruption priority system; then, on the basis of industrial load energy supply safety rating, industrial load safety indexes are designed, an industrial load peak regulation strategy is constructed, evaluation of industrial load energy supply safety is achieved, and clear guidance is provided for safety management of industrial load energy supply; and finally, the improved parrot optimization algorithm is applied to industrial load peak regulation, and the algorithm can find an optimal peak regulation scheme in a short time so as to realize optimization of industrial load safety indexes. According to the method, power resource distribution is optimized, key production links are guaranteed preferentially when power supply is insufficient, and the stability and safety of whole industrial production are maintained.

Description

Technical Field

[0001] The present invention belongs to the field of electrical engineering technology, and in particular relates to a large-capacity industrial load safety quantification and peak regulation method based on energy supply safety indicators. Background Art

[0002] Industrial load refers to the total amount of electricity consumed by various mechanical equipment, production facilities, and auxiliary systems during the industrial production process, including the electricity consumption required for production equipment operation, lighting, data processing, etc. Quantifying the energy supply security of industrial loads can identify and evaluate the risks and impacts that may be caused to production when the power supply is interrupted or unstable during the production process. By quantifying the safety boundaries of industrial loads, companies can develop effective preventive measures and emergency plans to reduce potential economic losses and safety accidents, while improving energy efficiency and the reliability of production systems. In addition, this approach helps to optimize the allocation of power resources and ensure that key production links are given priority when power supply is tight, thereby maintaining the stability and safety of the entire industrial production.

[0003] The core goal of the production process safety boundary quantification technology of industrial loads is to determine the minimum and maximum power demand ranges within which industrial loads can operate safely under different circumstances. This requires accurate measurement and analysis of power consumption in the industrial production process, taking into account the supply and demand balance of the power grid and the dynamics of the power market. During periods of tight power supply or peak demand, industrial enterprises may need to adjust their production plans to avoid exceeding the carrying capacity of the power grid. At the same time, through coordination with power grid operators, industrial enterprises can participate in demand response programs to obtain economic compensation or priority power supply rights by reducing or postponing some non-critical power consumption.

[0004] Quantifying the safety margins of industrial load energy supply is a multidisciplinary field that requires close cooperation between industrial enterprises and grid operators, using advanced technologies and data analysis methods to achieve optimization of power supply and risk management.

[0005] Therefore, there is an urgent need for a new technical solution in the existing technology to comprehensively consider large-capacity industrial loads to ensure the stable operation of industrial production and power systems. Summary of the invention

[0006] The technical problem to be solved by the present invention is to provide a large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators, which is used to solve the technical problem that industrial load safety indicators are difficult to quantify in the prior art.

[0007] The technical solution adopted by the present invention is to provide a large-capacity industrial load safety quantification and peak regulation method based on energy supply safety indicators, comprising the following steps:

[0008] S1: Construct the comprehensive weight of energy supply of industrial enterprises based on the sequential relationship analysis method and the indicator importance correlation method;

[0009] S2: The TOPSIS method based on comprehensive weights is used to determine the load interruption priority of industrial enterprises in the industrial park;

[0010] S3: Agreement on industrial load safety index: Construct an industrial load energy supply safety quantification method and an industrial load safety index optimization model based on industrial load safety index;

[0011] S4: Optimization of industrial load safety indicators based on the improved Parrot optimization algorithm.

[0012] The determination of the comprehensive weight of energy supply of industrial enterprises in step S1 includes the following steps:

[0013] S1.1: Evaluation index formulation: The energy supply characteristics of industrial enterprises u 1 , energy characteristics u 2 , load interruption willingness u 3 As an enterprise factor indicator, energy consumption 1 、Enterprise power reactive power incentives 2 As the energy supply factor, the enterprise factor index is used as the enterprise energy supply comprehensive weight input index, and the energy supply factor is used as the TOPSIS evaluation index;

[0014] S1.2: Use the G1 method to calculate the subjective weight of each evaluation indicator;

[0015] S1.3: Use the CRITIC method to calculate the objective weight of each evaluation indicator;

[0016] S1.4: Calculate the final comprehensive weight using the correction factor.

[0017] The specific calculation method of the comprehensive weight in step S1.4 is as follows:

[0018] (1) The evaluation indicators are arranged in order from small to large according to their objective weight;

[0019] (2) Calculate the difference coefficient and correction coefficient

[0020]

[0021] Among them, D c is the difference coefficient, C f is the correction coefficient, where the target weight vector Cw j (Cw 1 , Cw 2 , ..., Cw n ) have been rearranged in ascending order;

[0022] (3) Based on the calculated correction coefficient C f , calculate the comprehensive weight w of each evaluation index j :

[0023] w j =(1-C f )×Gw j +Cw j ×C f

[0024] Among them, Gw j Represents the subjective weight of the j-th evaluation indicator.

[0025] The method for determining the industrial load interruption priority in step S2 is specifically as follows:

[0026] S2.1: Positive index processing: unify the trend of the quality of each index, set all indicators to "the bigger the better", and for very small indicators, use the following formula to "positive" them

[0027]

[0028] Among them, x ij is the jth indicator after the i-th solution is positively transformed, x imax 、x imin is the maximum and minimum value of indicator j under scheme i, x ijk is the index before positive transformation, i represents the i-th scheme, j represents the j-th index, and the energy consumption volume s 1 、Enterprise power reactive power incentives 2 As an evaluation indicator of TOPSIS;

[0029] S2.2: Use the following calculation method to unify the dimensions of each indicator:

[0030]

[0031] Among them, z is the normalization matrix, z ij represents the element in the i-th row and j-th column of z, where the rows of z represent the solutions to be evaluated and the columns of z represent different evaluation indicators;

[0032] S2.3: Calculation of ideal and non-ideal solutions: The maximum value in each column of the matrix z with normalized indicators and unified dimensions is taken as the ideal solution, and the minimum value in each column is taken as the non-ideal solution;

[0033] Define the ideal solution z b The calculation is as follows:

[0034] z b =[z 1max z 2max ...zjmax ]

[0035] Define non-ideal solution z w as follows:

[0036] z w =[z 1min z 2min ...z jmin ]

[0037] Among them, z 1max 、z 1min Represents the maximum and minimum values ​​of the first indicator;

[0038] S2.4: Calculate the distance score of each solution: Calculate the distance between each solution and the ideal solution and the non-ideal solution

[0039] The ideal solution distance is calculated as follows:

[0040]

[0041] The distance of the non-ideal solution is as follows:

[0042]

[0043] Among them, d ib d iw are the distance between solution i and the ideal solution and the distance between solution i and the non-ideal solution, respectively, j is the comprehensive weight of indicator j;

[0044] S2.5: Scheme scoring construction, calculation of industrial enterprise scoring S i as follows:

[0045]

[0046] Rating: S i The larger the value, the better the solution. This is how we score each solution. i Will be used for subsequent load regulation; S of industrial enterprise production load i , the larger the load interrupt priority, S i The smaller the load interruption priority, the higher the priority, and industrial enterprises with lower load interruption priority will have priority in using electricity.

[0047] The industrial load safety index I in step S3 z It is the industrial load power consumption characteristic index I ch Industrial load power supply guarantee index I sa The weighted sum of industrial load safety index I z :

[0048] I z =wc I ch +w s I sa

[0049] Among them, w c is the weight of the power consumption characteristic index, which is set to 0.5 here, w s It is the weight of the electricity security indicator, which is set to 0.5 here.

[0050] The industrial load power supply guarantee index I sa :Based on the assessment of the interruptibility priority of industrial loads, set the industrial load power supply guarantee index I sa for:

[0051]

[0052] Among them, P i.t is the electricity consumption of the i-th industrial enterprise at time t, ΔP i.t is the load change of the i-th industrial enterprise at time t, ΔP i.tmax is the maximum load that the i-th industrial enterprise can change at time t, N is the number of industrial enterprises in the industrial park, S i Prioritize industrial enterprises;

[0053] ΔP i.t =P ia.t +P ib.t

[0054] ΔP i.tmax =max(P ia.tpmax ,|P ia.tnmax |)+|P ib.tmax |

[0055] Among them, P ia.t , P ib.t are the time-shifted load and reduced load of industrial enterprise i at time t; P ia.tpmax , P ia.tnmax , P ib.tmax are the maximum positive value, negative value and maximum load reduction of industrial enterprise i at time t respectively;

[0056] Time-shiftable load P ia.t The constraints are as follows:

[0057] P ia.tnmax ≤P ia.t ≤P ia.tpmax

[0058]

[0059] Where n is the number of time periods in a day. The above formula shows that the time-shiftable load has upper and lower bound constraints and load transfer constraints. When P ia.t When P is negative, it means that industrial enterprise i transfers load to other time periods at time t. ia.t The timing indicates that industrial enterprise i accepts the load transferred from other periods at time t;

[0060] Reducible load P ib.t There are transfer upper and lower bound constraints:

[0061] P ib.tmax ≤P ib.t ≤0

[0062] From the above formula, we can see that P ib.t If it is not a positive number, it means that the load that can be reduced is the load reduced from the original load;

[0063] Setting industrial load power supply guarantee index I sa The smaller the value, the better. i The higher the value, the more important the power supply of the industrial enterprise is. Therefore, when adjusting the industrial load, adjust S i Smaller industrial enterprise loads will receive a smaller industrial load power guarantee index I sa .

[0064] The industrial load power consumption characteristic index I ch The calculation method is as follows:

[0065]

[0066] Among them, I 1 is the maximum industrial load, I 2 is the peak-valley rate of industrial load, I 3 is the peak load duration, I 4 is the industrial load factor, I 1a is the maximum industrial load before peak load regulation, I 2a is the peak-to-valley ratio of industrial load before peak regulation, I 3a is the peak load duration before peak load regulation, I 4a is the industrial load factor before peak regulation;

[0067] Industrial load characteristic index I ch The design of the above-mentioned fractional form can eliminate the adverse effects caused by the different dimensions of different indicators and enhance the balance of each indicator. ch The smaller the better;

[0068] Industrial load maximum I 1 :

[0069] I 1 =max(P t )

[0070] Among them, P t is the industrial load time series, max(P t ) represents the maximum value of industrial load in the dispatching cycle. The maximum value of industrial load can evaluate the power supply capacity and distribution facility capacity required by industrial enterprises. t It is defined as the total electricity consumption of all industrial enterprises in the industrial park during period t, and is calculated as follows:

[0071]

[0072] Industrial load peak-valley rate I 2 :

[0073]

[0074] Among them, min(P t ) represents the minimum value of industrial load within the scheduling period. By calculating the peak-valley rate of industrial load, the load stability of industrial enterprises can be evaluated;

[0075] Peak load duration I 3 :

[0076]

[0077] I 3 =∑t i

[0078] Among them, t i is a binary variable for determining peak load, P t Defined as peak load, I 3 is the total period of peak load for industrial load, P line.max It is the maximum power provided by the upper power grid to the industrial park. The peak load duration shows the duration of the industrial load exceeding a certain value in a specific time period.

[0079] Industrial load factor I 4 :

[0080]

[0081] The industrial load factor is defined as the ratio of average load to maximum load and is used to reflect the stability of industrial load.

[0082] In step S4, the Parrot optimization algorithm is improved and applied to the optimization of industrial load safety index, the population size and the number of iterations are set, the time-shiftable load and the load that can be reduced are defined as decision variables, and the upper and lower bounds of the load changes of each enterprise are counted as the boundary constraints of the decision variables. The objective function is the industrial load safety index I z, the fitness value is the objective function value.

[0083] The improved parrot optimization algorithm in step S4 comprises the following steps:

[0084] S4.1: Population parameter definition;

[0085] S4.2: Initialize the population;

[0086] S4.3: Enter the iteration loop;

[0087] S4.4: Randomly select iteration mode;

[0088] S4.5: Use the Cauchy mutation factor to perturb the individuals in the population;

[0089] S4.6: Individual decision variable modification;

[0090] S4.7: Determine whether the iteration is completed. If the iteration is not completed, repeat steps 3 to 7;

[0091] S4.8: Complete the iteration and output the optimal industrial load safety index and the optimal load change value of each enterprise.

[0092] Through the above design scheme, the present invention can bring the following beneficial effects:

[0093] Quantifying the safety margin of industrial loads during the production process can evaluate the priority of industrial enterprise load interruptions on the energy consumption and supply sides and quantify the energy consumption characteristics of industrial enterprises. By implementing industrial load peak regulation plans, while optimizing load indicators, the willingness of industrial enterprises to interrupt loads can be fully considered, which is conducive to ensuring the energy safety of industrial enterprises, reducing grid load fluctuations and grid overload risks, weakening the peak-to-valley difference of the grid, and improving the stability, reliability and safety of the grid. BRIEF DESCRIPTION OF THE DRAWINGS

[0094] Figure 1 It is a flow chart of the safety boundary quantification of the production process of industrial loads of the large-capacity industrial load safety quantification and peak regulation method based on the energy supply safety index of the present invention;

[0095] Figure 2 It is a flow chart of the improved Parrot optimization algorithm of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention;

[0096] Figure 3 Example 1 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention is a load change diagram of enterprise 1;

[0097] Figure 4 The load variation diagram of enterprise 2 in Example 1 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0098] Figure 5 This is a load variation diagram of enterprise 3 in Example 1 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention;

[0099] Figure 6 This is a load variation diagram of enterprise 4 in Example 1 of a large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention;

[0100] Figure 7 The load variation diagram of an industrial park in Example 1 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention;

[0101] Figure 8 The load variation diagram of enterprise 1 in Example 2 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0102] Fig. 9 It is a load variation diagram of enterprise 2 of embodiment 2 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0103] Fig.10 It is a load variation diagram of enterprise 3 in Example 2 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0104] Fig.11 It is a load variation diagram of enterprise 4 in Example 2 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0105] Fig.12 The load variation diagram of an industrial park in Example 2 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention;

[0106] Fig.13 It is a load variation diagram of enterprise 1 in Example 3 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0107] Fig.14 It is a load variation diagram of Enterprise 2 in Example 3 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0108] Fig.15 It is a load variation diagram of enterprise 3 of embodiment 3 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index of the present invention;

[0109] Fig.16 This is an industrial park load variation diagram of Example 3 of the large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators of the present invention. DETAILED DESCRIPTION

[0110] The present invention is further described below with reference to the accompanying drawings and specific embodiments:

[0111] Example 1: Example 1 uses spring load data of an industrial park in Shandong Province, which includes four industrial enterprises. The selected industrial enterprises are steel, electrolytic aluminum, cement and textile enterprises. The following steps are followed to perform industrial load interruption priority negotiation and industrial load safety index optimization.

[0112] S01: Evaluation index formulation: The energy supply characteristics of industrial enterprises u 1 , energy characteristics u 2 , load interruption willingness u 3 As an enterprise factor indicator. Energy consumption 1 、Enterprise power reactive power incentives 2 As a factor of the energy supply side. The enterprise factor index is used as the comprehensive weight input index of the energy supply of the enterprise; the energy supply side factor is used as the evaluation index of TOPSIS.

[0113] S02: Use the G1 method to calculate the subjective weight of each evaluation index. The G1 method is also called the sequential relationship analysis method, which is a commonly used subjective weighting method. The calculation process of this method is simple and fast, and it does not require consistency verification like the hierarchical analysis method. However, the G1 method has some limitations, such as strong dependence on the judgment matrix and high requirements for consistency verification. The G1 method process is as follows:

[0114] 1. Construct evaluation indicators: Determine the evaluation indicators based on the evaluation object.

[0115] 2. Data standardization: To eliminate the dimensional impact between different indicators, the data of each indicator is standardized. Commonly used standardization methods include range standardization and normalization.

[0116] 3. Construct a judgment matrix: Construct a judgment matrix. The elements in the matrix represent the relative importance of indicators, which are usually determined by expert scoring or other methods.

[0117] 4. Calculate the weight vector: Use the judgment matrix to calculate the weight vector of each indicator.

[0118] 5. Perform consistency check: Perform consistency check on the judgment matrix to ensure the rationality of the expert scoring. If the consistency check fails, the judgment matrix needs to be readjusted.

[0119] 6. Calculate the comprehensive score of the indicator: Combine the standardized indicator data with the weight vector to calculate the comprehensive score of each evaluation object.

[0120] 7. Sorting: Sort the evaluation objects according to their comprehensive scores to determine their positions in the evaluation system.

[0121] S03: Use the CRITIC method to calculate the objective weight of each evaluation indicator. The CRITIC method is also called the indicator importance correlation method, which can be used to deal with the objective weighting problem under multiple standards. This method mainly determines the objective weight of the indicator by analyzing the objective laws between the indicators. However, the results of the CRITIC method are greatly affected by data distribution and sample size. The CRITIC method process is as follows:

[0122] 1. Determine the evaluation indicators: Select the evaluation indicators.

[0123] 2. Data standardization: Eliminate the impact of different indicator dimensions and ensure that the data has the same order of magnitude.

[0124] 3. Calculate standard deviation: Calculate the standard deviation for each indicator.

[0125] 4. Calculate the correlation coefficient: Calculate the correlation coefficient between indicators to measure the conflict between indicators.

[0126] 5. Determine the conflict index of the indicator: Calculate the conflict index of each indicator based on the standard deviation and correlation coefficient.

[0127] 6. Determine the distinguishing ability of the indicators: Calculate the distinguishing ability index based on the correlation coefficient between the indicators. The distinguishing ability index is usually expressed as 1 minus the absolute value of the correlation coefficient.

[0128] 7. Calculate CRITIC weight: Combine the conflict index and the discrimination index to calculate the CRITIC weight of each indicator. The CRITIC weight is the product of the conflict index and the discrimination index.

[0129] 8. Indicator ranking: Indicators are ranked according to CRITIC weights. Indicators with larger weights are more important.

[0130] 9. Result analysis: Analyze the indicator importance ranking results and identify key indicators and secondary indicators.

[0131] S04: Calculate the final comprehensive weight of the enterprise factor index through the correction coefficient. The comprehensive weight of each industrial enterprise is shown in the following table:

[0132] Table 1 Comprehensive weights of enterprises in Example 1

[0133] Number of enterprises Enterprise 1 weight Enterprise 2 Weight Enterprise 3 weight Enterprise 4 weight Comprehensive weight 0.18 0.32 0.13 0.37

[0134] S05: Use the approximate ideal solution (TOPSIS) based on the energy supply factors to determine the priority of industrial load interruption. Calculate the score S of the industrial enterprise i, the larger the load interrupt priority, S i The smaller the load interruption priority, the higher the load interruption priority. Industrial enterprises with lower load interruption priority will be given priority to use electricity. Based on the score, the load interruption priority order of each enterprise is as follows:

[0135] Table 2 Example 1 Enterprise Priority Order

[0136]

[0137] From the results in the above table, it can be seen that when load adjustment is carried out, the power supply of enterprise 4 will be guaranteed first, that is, the load adjustment of enterprise 4 will be avoided as much as possible.

[0138] S06: Construction of industrial load power supply guarantee index: Industrial load power supply guarantee index I sa :Based on the assessment of the interruptibility priority of industrial loads, the power supply guarantee indicators for industrial loads are set as follows:

[0139]

[0140] Among them, P i.t is the electricity consumption of the i-th industrial enterprise at time t. ΔP i.t is the load change of the i-th industrial enterprise at time t. ΔP i.tmax is the maximum load that the i-th industrial enterprise can change at time t. N is the number of industrial enterprises in the industrial park. i Score the priority of industrial enterprises.

[0141] ΔP i.tmax With ΔP i.t The calculation method is as follows:

[0142] ΔP i.t =P ia.t +P ib.t

[0143] ΔP i.tmax =max(P ia.tpmax ,|P ia.tnmax |)+|P ib.tmax |

[0144] Among them, P ia.t , P ib.t are the time-shifted load and reduced load of industrial enterprise i at time t; P ia.tpmax , P ia.tnmax , P ib.tmax are respectively the maximum positive value, negative value and maximum load reduction of industrial enterprise i at time t.

[0145] Time-shiftable load P ia.t The constraints are as follows:

[0146] P ia.tnmax ≤P ia.t ≤P ia.tpmax

[0147]

[0148] Where n is the number of time periods in a day. The above formula shows that the time-shiftable load has upper and lower bounds and load transfer constraints. ia.t When P is negative, it means that industrial enterprise i transfers load to other time periods at time t. ia.t It is the timing that industrial enterprise i accepts the load transferred from other periods at time t.

[0149] Reducible load P ib.t There are transfer upper and lower bound constraints:

[0150] P ib.tmax ≤P ib.t ≤0

[0151] From the above formula, we can see that P ib.t It is not a positive number, which means that the load that can be reduced is the load reduced from the original load.

[0152] Setting industrial load power supply guarantee index I sa The smaller the value, the better. i The higher the value, the more important the power supply of the industrial enterprise is. Therefore, when adjusting the industrial load, adjust S i Smaller industrial enterprise loads will receive a smaller industrial load power guarantee index I sa .

[0153] S07: Construction of industrial load electricity consumption characteristic indicators: Industrial load electricity consumption characteristic indicators are indicators formulated based on the characteristics of the industrial load itself. They can evaluate the impact of industrial load on the upper-level power grid and reflect the effect of industrial load management.

[0154] Industrial load power consumption characteristic index I ch The calculation method is as follows:

[0155]

[0156] Among them, I 1 is the maximum industrial load, I 2 is the peak-valley rate of industrial load, I 3 is the peak load duration, I 4 is the industrial load factor. 1a is the maximum industrial load before peak load regulation, I 2a is the peak-to-valley ratio of industrial load before peak regulation, I 3a is the peak load duration before peak load regulation, I 4a is the industrial load factor before peak regulation.

[0157] Industrial load characteristic index I ch The above-mentioned fractional form can eliminate the adverse effects of different indicators and enhance the balance of each indicator. ch The smaller the better.

[0158] Industrial load maximum I 1 :

[0159] I 1 =max(P t )

[0160] Among them, P t is the industrial load time series, max(P t ) represents the maximum value of industrial load in the dispatching cycle. The maximum value of industrial load can evaluate the power supply capacity and distribution facility capacity required by industrial enterprises. t It is defined as the total electricity consumption of all industrial enterprises in the industrial park during period t, and is calculated as follows:

[0161]

[0162] Industrial load peak-valley rate I 2 :

[0163]

[0164] Among them, min(P t ) represents the minimum value of industrial load in the scheduling cycle. By calculating the peak-to-valley rate of industrial load, the load stability of industrial enterprises can be evaluated.

[0165] Peak load duration I 3 :

[0166]

[0167] I 3 =∑t i

[0168] Among them, t i is a binary variable for determining peak load. t Defined as peak load. 3 is the total period of time when the industrial load is at peak load. line.max It is the maximum power provided by the upper power grid to the industrial park. The peak load duration shows the duration that the industrial load exceeds a certain value in a specific time period.

[0169] Industrial load factor I 4 :

[0170]

[0171] The industrial load factor is defined as the ratio of the average load to the maximum load. The industrial load factor is often used to reflect the stability of industrial load.

[0172] S08: Constructing Industrial Load Safety Index I z As shown below:

[0173] I z =w c I ch +w s I sa

[0174] Among them, w c is the weight of the power consumption characteristic index, which is set to 0.5 here, w s is the weight of the power supply guarantee indicator, which is set to 0.5 here. z It is an industrial load safety indicator.

[0175] From the above formula, we can see that adjusting the industrial load can indirectly improve the power consumption characteristic index, but it is not conducive to the power consumption guarantee index. Therefore, appropriate adjustment of the industrial load can achieve a compromise solution that takes into account both goals, improve the industrial load safety index, and thus ensure the power safety of the industrial park.

[0176] S09: Improve the Parrot optimization algorithm and apply it to the optimization of industrial load safety indicators. First, define the parameters: set the population size and the number of iterations. Define the time-shiftable load and the load that can be reduced as decision variables, and count the upper and lower bounds of the load changes of each enterprise as the boundary constraints of the decision variables. The objective function is the industrial load safety index I z The fitness value is the objective function value (industrial load safety index I z )

[0177] S10: Initialize the population: Initialize the population in the following way:

[0178] x i.0 =lb+rand(ub-lb)

[0179] Among them, x i.0 is the first generation individual, lb is the lower bound of the variable, ub is the upper bound of the variable, and rand is a random number in the range of 0 to 1.

[0180] S11: Enter the iteration loop: calculate the fitness of each individual, and record the decision variables and fitness of the best individual of each generation and the global best individual.

[0181] S12: Randomly select an iteration mode. Each individual in the population controls the variables according to the iteration mode and randomly selects one of the four iteration modes of foraging, staying, communicating, and fear of strangers to iterate.

[0182] Foraging updates are as follows:

[0183]

[0184] Among them, x i.t+1 , x i.t is the decision variable of individual i in the population at the t+1th, tth iteration. best , x mean.t is the global optimal individual and the average position of the population at the t-th iteration, t is the t-th iteration, ger is the total number of iterations. Levy is the Levy coefficient, which is calculated as follows:

[0185]

[0186] μ~N(0,1)

[0187] v~N(0,1)

[0188]

[0189] Among them, μ and ν are random numbers uniformly distributed between 0 and 1. β is the Levy coefficient, which is a constant in the range of 0 to 2.

[0190] The dwell behavior is updated as follows:

[0191] x i.t+1 =x i.t +x best ·Levy+rand

[0192] The communication behavior is updated as follows:

[0193]

[0194] Wherein, P is a random number in the range of 0 to 1. Different update methods are selected by judging the range of the random number P.

[0195] The fear behavior is updated as follows:

[0196]

[0197] S13: To enhance the performance of the algorithm, the Cauchy mutation factor is used to perturb the individuals in the population. The Cauchy mutation strategy is as follows:

[0198]

[0199] Among them, Cauchy(0,1) is the Cauchy distribution, rate is the mutation probability, and rand is a random number between 0 and 1.

[0200] S14: Individual decision variable correction. If the decision variable exceeds the upper and lower limits of the change, the variable is corrected

[0201] S15: Determine whether the iteration is completed. If not, repeat Figure 2 Steps 3 to 7 in the tutorial.

[0202] S16: Complete the iteration and output the optimal industrial load safety index and the optimal load change value of each enterprise. The result is as follows Figure 3 to Figure 7 shown.

[0203] The code flow for improving the parrot optimization algorithm is as follows:

[0204]

[0205]

[0206] Example 2: The experiment selected summer load data of an industrial park in Hebei Province. The industrial park includes 4 industrial enterprises, all of which are cement production enterprises. The specific implementation steps are similar to those of Example 1. The comprehensive weight of each industrial enterprise is shown in Table 3, the priority order of each enterprise is shown in Table 4, and the load scheduling diagram is shown in Table 4. Figure 8 to Figure 12 shown.

[0207] Table 3 Comprehensive weights of enterprises in Example 2

[0208] Number of enterprises Enterprise 1 weight Enterprise 2 Weight Enterprise 3 weight Enterprise 4 weight Comprehensive weight 0.47 0.15 0.21 0.17

[0209] Table 4 Enterprise priority order in Example 2

[0210]

[0211] Example 3: The experiment selected load data of an industrial zone in Liaoning. The load data used autumn load data. The industrial zone includes three industrial enterprises. Enterprises 1 and 2 are both steel enterprises, and enterprise 3 is an automobile manufacturer. The specific implementation steps are similar to those in Example 1, but only TOPSIS is used to evaluate the energy supply index, without considering the comprehensive weight of the enterprise-side index. The priority of industrial enterprises is shown in Table 5, and the load scheduling diagram is shown in Figure 13 to Figure 16 shown.

[0212] Table 5 Enterprise priority order in Example 3

[0213]

[0214] The implementation methods of the present invention are not limited to the above-mentioned embodiments, and any other changes, modifications, substitutions, combinations, and simplifications made without departing from the spirit and principle of the present invention should be equivalent replacement methods and are included in the protection scope of the present invention.

Claims

1. A large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety indicators, characterized by: The following steps are involved: S1: Construct the comprehensive weight of energy supply of industrial enterprises based on the sequential relationship analysis method and the indicator importance correlation method; S2: The TOPSIS method based on comprehensive weights is used to determine the load interruption priority of industrial enterprises in the industrial park; S3: Agreement on industrial load safety index: Construct an industrial load energy supply safety quantification method and an industrial load safety index optimization model based on industrial load safety index; S4: Optimization of industrial load safety indicators based on the improved Parrot optimization algorithm.

2. The large-capacity industrial load safety quantification and peak-shaving method based on energy supply safety index according to claim 1 is characterized by: The determination of the comprehensive weight of energy supply of industrial enterprises in step S1 includes the following steps: S1.1: Evaluation index formulation: The energy supply characteristics u1, energy consumption characteristics u2, and load interruption willingness u3 of industrial enterprises are used as enterprise factor indicators, energy consumption volume s1 and enterprise power consumption reactive power incentive s2 are used as energy supply factors, enterprise factor indicators are used as comprehensive weight input indicators of enterprise energy supply, and energy supply factors are used as TOPSIS evaluation indicators; S1.2: Use the G1 method to calculate the subjective weight of each evaluation indicator; S1.3: Use the CRITIC method to calculate the objective weight of each evaluation indicator; S1.4: Calculate the final comprehensive weight using the correction factor.

3. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 2 is characterized in that: The specific calculation method of the comprehensive weight in step S1.4 is as follows: (1) The evaluation indicators are arranged in order from small to large according to their objective weight; (2) Calculate the difference coefficient and correction coefficient Among them, D c is the difference coefficient, C f is the correction coefficient, where the target weight vector Cw j (Cw1, Cw2, ..., Cw n ) have been rearranged in ascending order; (3) Based on the calculated correction coefficient C f , calculate the comprehensive weight w of each evaluation index j : w j =(1-C f )×Gw j +Cw j ×C f Among them, Gw j Represents the subjective weight of the jth evaluation indicator.

4. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 1 is characterized in that: The method for determining the industrial load interruption priority in step S2 is specifically as follows: S2.1: Positive index processing: unify the trend of the quality of each index, set all indicators to "the bigger the better", and for very small indicators, use the following formula to "positive" them Among them, x ij is the jth indicator after the i-th solution is positively transformed, x imax 、x imin is the maximum and minimum value of indicator j under scheme i, x ijk is the index before positive transformation, i represents the i-th scheme, j represents the j-th index, and the energy consumption volume s1 and the enterprise power reactive power incentive s2 are used as the evaluation indicators of TOPSIS; S2.2: Use the following calculation method to unify the dimensions of each indicator: Among them, z is the normalization matrix, z ij represents the element in the i-th row and j-th column of z, where the rows of z represent the solutions to be evaluated and the columns of z represent different evaluation indicators; S2.3: Calculation of ideal and non-ideal solutions: The maximum value in each column of the matrix z with normalized indicators and unified dimensions is taken as the ideal solution, and the minimum value in each column is taken as the non-ideal solution; Define the ideal solution z b The calculation is as follows: With b =[of 1max With 2max ...With jmax ] Define the non-ideal solution z w as follows: With w =[of 1min With 2min ...With jmin ] Among them, z 1max 、z 1min Represents the maximum and minimum values ​​of the first indicator; S2.4: Calculate the distance score of each solution: Calculate the distance of each solution from the ideal solution and the non-ideal solution. The distance calculation method for the ideal solution is as follows: The distance of the non-ideal solution is as follows: Among them, d ib ,d iw are the distance between solution i and the ideal solution and the distance between solution i and the non-ideal solution, respectively, j is the comprehensive weight of indicator j; S2.5: Scheme scoring construction, calculation of industrial enterprise scoring S i as follows: Rating: S i The larger the value, the better the solution. This is how we score each solution. i Will be used for subsequent load regulation; S of industrial enterprise production load i , the larger the load interrupt priority, S i The smaller the load interruption priority, the higher the priority, and industrial enterprises with lower load interruption priority will have priority in using electricity.

5. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 1 is characterized in that: The industrial load safety index I in step S3 z It is the industrial load power consumption characteristic index I ch Industrial load power supply guarantee index I sa The weighted sum of industrial load safety index I z : AND z =in c AND ch +in s AND sa Among them, w c is the weight of the power consumption characteristic index, which is set to 0.5 here, w s It is the weight of the electricity security indicator, which is set to 0.5 here.

6. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 5 is characterized by: The industrial load power supply guarantee index I sa :Based on the assessment of the interruptibility priority of industrial loads, set the industrial load power supply guarantee index I sa for: Among them, P i.t is the electricity consumption of the i-th industrial enterprise at time t, ΔP i.t is the load change of the i-th industrial enterprise at time t, ΔP i.tmax is the maximum load that the i-th industrial enterprise can change at time t, N is the number of industrial enterprises in the industrial park, S i Prioritize industrial enterprises; ΔP i.t =P ia.t +P ib.t ΔP i.tmax =max(P ia.tpmax ,|P ia.tnmax |)+|P ib.tmax | Among them, P ia.t , P ib.t are the time-shifted load and reduced load of industrial enterprise i at time t; P ia.tpmax , P ia.tnmax , P ib.tmax are the maximum positive value, negative value and maximum load reduction of industrial enterprise i at time t respectively; Time-shiftable load P ia.t There are the following constraints: P ia.tnmax ≤P ia.t ≤P ia.tpmax Where n is the number of time periods in a day. The above formula shows that the time-shiftable load has upper and lower bound constraints and load transfer constraints. When P ia.t When P is negative, it means that industrial enterprise i transfers load to other time periods at time t. ia.t The timing indicates that industrial enterprise i accepts the load transferred from other periods at time t; Reducible load P ib.t There are transfer upper and lower bound constraints: P ib.tmax ≤P ib.t ≤0 From the above formula, we can see that P ib.t If it is not a positive number, it means that the load that can be reduced is the load reduced from the original load; Setting industrial load power supply guarantee index I sa The smaller the value, the better. i The higher the value, the more important the power supply of the industrial enterprise is. Therefore, when adjusting the industrial load, adjust S i Smaller industrial enterprise loads will receive a smaller industrial load power guarantee index I sa .

7. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 5 is characterized in that: The industrial load power consumption characteristic index I ch The calculation method is as follows: Among them, I1 is the maximum industrial load, I2 is the peak-to-valley rate of industrial load, I3 is the duration of peak load, I4 is the industrial load factor, I 1a is the maximum industrial load before peak load regulation, I 2a is the peak-to-valley ratio of industrial load before peak regulation, I 3a is the peak load duration before peak load regulation, I 4a is the industrial load factor before peak regulation; Industrial load characteristic index I ch The design of the above-mentioned fractional form can eliminate the adverse effects caused by the different dimensions of different indicators and enhance the balance of each indicator. ch The smaller the better; Maximum industrial load I1: I1=max(P t ) Among them, P t is the industrial load time series, max(P t ) represents the maximum value of industrial load in the dispatching cycle. The maximum value of industrial load can evaluate the power supply capacity and distribution facility capacity required by industrial enterprises. t It is defined as the total electricity consumption of all industrial enterprises in the industrial park during period t, and is calculated as follows: Industrial load peak-valley rate I2: Among them, min(P t ) represents the minimum value of industrial load within the scheduling period. By calculating the peak-valley rate of industrial load, the load stability of industrial enterprises can be evaluated; Peak load duration I3: I3=∑t i Among them, t i is a binary variable for determining peak load, P t Defined as peak load, I3 is the total period of industrial load as peak load, P line.max It is the maximum power provided by the upper power grid to the industrial park. The peak load duration shows the duration of the industrial load exceeding a certain value in a specific time period. Industrial load factor I4: The industrial load factor is defined as the ratio of average load to maximum load and is used to reflect the stability of industrial load.

8. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 6 is characterized in that: In step S4, the Parrot optimization algorithm is improved and applied to the optimization of industrial load safety index, the population size and the number of iterations are set, the time-shiftable load and the load that can be reduced are defined as decision variables, and the upper and lower bounds of the load changes of each enterprise are counted as the boundary constraints of the decision variables. The objective function is the industrial load safety index I z , the fitness value is the objective function value.

9. The method for safe quantification and peak load regulation of large-capacity industrial loads based on energy supply safety indicators according to claim 8 is characterized in that: The improved parrot optimization algorithm in step S4 comprises the following steps: S4.1: Population parameter definition; S4.2: Initialize the population; S4.3: Enter the iteration loop; S4.4: Randomly select iteration mode; S4.5: Use the Cauchy mutation factor to perturb the individuals in the population; S4.6: Individual decision variable modification; S4.7: Determine whether the iteration is completed. If the iteration is not completed, repeat steps 3 to 7; S4.8: Complete the iteration and output the optimal industrial load safety index and the optimal load change value of each enterprise.

Images