Multi-region interconnected power grid resource scheduling method and system based on big data

Abstract

The invention relates to the technical field of resource scheduling, and discloses a multi-region interconnected power grid resource scheduling method and system based on big data, and the method comprises the steps: obtaining the number of municipal areas and the power consumption data of each municipal area, and calculating the power utilization change rate; obtaining substation parameters of the power supply line, and performing power adjustment; and updating the power consumption data in real time, and adjusting the power distribution of each municipal district. Electric quantity distribution can be dynamically optimized by acquiring electric quantity data of a municipal district, substation parameters of a power supply line and a power adjustment mechanism in real time, so that the operation efficiency of an electric power system is effectively improved. A model based on power consumption change rate and power factor adjustment is adopted, accurate adjustment of the power grid load is achieved, and efficient utilization of electric power resources and power supply stability are ensured.

Description

Technical Field

[0001] The present invention relates to the technical field of resource scheduling, and in particular to a method and system for scheduling multi-regional interconnected power grid resources based on big data. Background Art

[0002] With the continuous advancement of urbanization and the rapid development of the economy, the demand for electricity in cities has shown an increasing trend year by year. In order to meet the growing electricity demand in cities, the planning, construction and management of power systems have become particularly important. Especially in large cities and urban areas, the refined management of power distribution has become an important means to improve energy efficiency and reduce energy waste. Existing power dispatching systems mostly rely on traditional load forecasting methods, such as time series analysis and linear regression. However, these methods often cannot provide sufficient accuracy and flexibility when faced with complex changes in electricity demand, the diversity of electricity consumption characteristics, and the mutual influence of many variables in the power system.

[0003] In order to improve the accuracy and timeliness of power distribution, many researchers have tried to apply data-driven intelligent algorithms, such as machine learning and deep learning methods, to optimize the dispatch and load forecasting of power systems. These methods can identify power consumption patterns and trends by analyzing historical data and make dynamic adjustments based on real-time data. However, traditional machine learning algorithms may face problems such as high data dimensions and high computational complexity when processing multi-dimensional power data. Therefore, how to reduce computing resource consumption while ensuring the effectiveness of the algorithm has become one of the main challenges facing current power dispatch optimization technology.

[0004] In terms of power distribution in urban areas, traditional methods usually rely on static power calculations and fixed power distribution strategies, which are often unable to be flexibly adjusted according to real-time changing needs. Unbalanced power distribution can lead to oversupply or undersupply in some areas, resulting in resource waste or power outages. Therefore, how to update the power distribution in urban areas based on real-time data and reasonably adjust the power load in each urban area has become a key issue that needs to be urgently addressed in the power system. Summary of the invention

[0005] In view of the above-mentioned problems, the present invention is proposed.

[0006] Therefore, the technical problem solved by the present invention is: in an urban power system, how to dynamically adjust the power distribution in the urban area according to real-time data to improve the efficiency, stability and energy consumption management of the power system.

[0007] In order to solve the above technical problems, the present invention provides the following technical solutions: a multi-region interconnected power grid resource scheduling method based on big data, comprising:

[0008] Obtain the number of municipal districts and the power consumption data of each municipal district, and calculate the power consumption change rate;

[0009] Obtain substation parameters of the power supply line and adjust power;

[0010] Update power consumption data in real time and adjust power distribution in each city district.

[0011] As a preferred solution of the multi-region interconnected power grid resource dispatching method based on big data described in the present invention, wherein: the obtaining of the number of municipal districts and the power consumption data of each municipal district includes obtaining the specific number of municipal districts and regional division information from the power management system;

[0012] Real-time and historical electricity consumption data of each district are collected from electricity meters, smart grid monitoring equipment and power metering systems in each municipal district.

[0013] As a preferred solution of the multi-regional interconnected power grid resource dispatching method based on big data described in the present invention, wherein: calculating the amount of electricity to be allocated includes obtaining the current month as m and the current date as d;

[0014] Get the total number of days in month m, recorded as dt;

[0015] The number of all municipal districts in the target area, denoted as mn;

[0016] Get the current power consumption corresponding to each city district. The formula is:

[0017] (oj1,oj2...oj mn )

[0018] Among them, oj1 represents the current power consumption of the first district; oj2 represents the current power consumption of the second district; oj mn represents the current power consumption of the mnth municipal district;

[0019] Obtain the electricity consumption of each city district in month m in the past three years as the historical electricity consumption;

[0020] The historical electricity consumption of the first municipal district is expressed as ahj (1) 、bhj (1) and chj (1) ;

[0021] Among them, ahj (1) represents the electricity consumption of the first municipal district in the third year in month m, bhj (1) represents the power consumption of the first municipal district in the second year in month m, chj (1) It represents the electricity consumption in month m of the first year in the first municipal district;

[0022] The historical electricity consumption of the second municipal district is ahj(2) 、bhj (2) and chj (2) ;

[0023] Among them, ahj (2) represents the electricity consumption of the second urban district in the third year in month m, bhj (2) represents the power consumption of the second district in the second year in month m, chj (2) It represents the electricity consumption in month m of the first year of the second municipal district;

[0024] The historical electricity consumption of the mnth municipal district is ahj (mn) 、bhj (mn) and chj (mn) ;

[0025] Among them, ahj (mn) represents the power consumption of the mnth district in the third year in month m, bhj (mn) represents the power consumption of the mnth district in the second year in month m, chj (mn) It represents the electricity consumption of the mth district in the first year and month m;

[0026] According to the historical power consumption, the formula for calculating the power consumption change rate of each municipal district is expressed as:

[0027]

[0028] Wherein, ub(1) represents the electricity consumption change rate of the first city district; ub(2) represents the electricity consumption change rate of the second city district; (mn) represents the electricity consumption change rate of the mnth city district; E i,t represents the electricity consumption of city district i at time t, E i,t-1 It represents the electricity consumption of urban district i in the previous time period.

[0029] As a preferred solution of the multi-regional interconnected power grid resource scheduling method based on big data described in the present invention, the calculation of the power consumption change rate also includes obtaining the total power consumption of each municipal district last month, which is expressed as a reference power formula:

[0030] (hj1,hj2...hj mn )

[0031] Where hj1 represents the total power consumption of the first city district last month; hj2 represents the total power consumption of the second city district last month; hj mn It represents the total electricity consumption of the mnth municipal district last month;

[0032] Calculate the amount of electricity to be distributed in the first municipal district. The formula is:

[0033]

[0034] Calculate the amount of electricity to be distributed in the second municipal district, the formula is:

[0035]

[0036] The formula for calculating the amount of electricity to be distributed in the mnth municipal district is:

[0037]

[0038] Among them, pj1, pj2, …, pj mn represents the amount of electricity to be allocated in the jth municipal district, j represents the number of the municipal district, mn represents the total number of municipal districts, ub (1) represents the electricity consumption change rate of the first city district; ub (2) represents the electricity consumption change rate of the second urban district; ub (mn) represents the rate of change of electricity consumption in the mnth district; dt represents the time interval, oj1 represents the change of electricity consumption in the first district; oj2 represents the change of electricity consumption in the second district; oj mn represents the change in electricity consumption in the mnth municipal district; apj represents the sum of the electricity to be allocated in all municipal districts;

[0039] The active power of the power station in the target area is obtained as PO, and the reactive power is recorded as QO;

[0040] The formula for calculating the power factor PFO of a power station is:

[0041] PFO=PO / QO

[0042] Calculate the total power generation of the power station in month m, denoted as ZJ formula:

[0043] ZJ=PO×tt

[0044] Among them, tt represents the power generation time of the power station in month m; PO represents the active power of the power station in the target area, ZJ represents the total power generation of the power station in month m; QO represents the reactive power of the power station in the target area, and PFO represents the power factor of the power station;

[0045] When ZJ≤apj, the amount of electricity generated by the power station is low, the power of the power station is increased, and the active power and reactive power of the generator are adjusted;

[0046] When ZJ>[apj×(1+α)], the power generated by the power station is sufficient and the power of the power station is not adjusted;

[0047] Wherein, α represents the power loss rate.

[0048] As a preferred solution of the multi-regional interconnected power grid resource scheduling method based on big data described in the present invention, power adjustment includes: (1) 、bhj (1) and chj (1) ;

[0049] Calculate the rate of change of power consumption in the first city district, denoted as ub (1) ;

[0050] Calculate the rate of change in electricity consumption from the third year to the second year in the first municipal district, recorded as up (3-2) The formula is:

[0051] up (3-2) =[bhj (1) / chj (1) ]-1

[0052] Calculate the rate of change in electricity consumption from the second year to the first year in the first city district, recorded as up (2-1) The formula is:

[0053] up (2-1) =[ahj (1) / bhj (1) ]-1

[0054] According to up (3-2) and up (2-1) , construct the state transfer matrix, and use the NumPy library function to iterate the state transfer matrix to obtain the steady-state transfer matrix B1;

[0055] Define a 2×1 order receiving matrix {0, 1}, denoted as matrix B2;

[0056] Use matrix multiplication to calculate the matrix B1 multiplied by the matrix B2 to obtain the matrix B3;

[0057] Calculate the average value of each parameter in matrix B3 as the power consumption change rate of the first district, and get ub (1) ;

[0058] Calculate the power consumption change rate of the 2nd to mnth districts, and get ub (2) ) (mn) ;

[0059] The adjusted generator active power is calculated and recorded as PO1. The formula is:

[0060] PO1=PO+[(apj-ZJ) / tt]

[0061] [(apj-ZJ) / tt]≥(PO×α)

[0062] Among them, PO1 represents the adjusted active power of the generator, up (3-2) Indicates the rate of change of electricity consumption from the third year to the second year in the first municipal district; bhj (1) Indicates the power consumption of the first municipal district in the second or third year; chj (1) Indicates the power consumption of the first municipal district in the third or second year; up (2-1) Indicates the rate of change of electricity consumption from the second year to the first year in the first city district, ahj (1) represents the power consumption of the first municipal district in the first year, B1 represents the steady-state transfer matrix; B2 represents the acceptance matrix, B3 represents the matrix multiplication sum, PO represents the active power of the power station in the target area; apj represents the total amount of power to be allocated, ZJ represents the total power generation of the power station in month m, and tt represents the power generation time of the power station in month m;

[0063] Assuming the reactive power of the power station remains unchanged, calculate the new power factor of the power station, denoted as PFO1, and the formula is expressed as:

[0064] PFO1=PO1 / QO

[0065] Assuming the power factor of the power station remains unchanged, calculate the new reactive power of the power station, denoted as QO1, and the formula is expressed as:

[0066] QO1=PO1 / FPO

[0067] Among them, FPO represents the power factor; PO1 represents the adjusted active power of the generator, QO1 represents the new reactive power of the power station, PFO1 represents the new power factor of the power station, and QO represents the reactive power of the power station in the target area;

[0068] Calculate the rate of change of reactive power of the power station, denoted as b (QO) The formula is:

[0069] b (QO) =(QO1-QO) / QO

[0070] Calculate the rate of change of the power factor of the power station, denoted as b (PFO) The formula is:

[0071] b (PFO) =(PFO1-PFO) / PFO

[0072] Among them, b (PFO) represents the rate of change of the power factor of the power station, PFO1 represents the new power factor of the power station, PFO represents the power factor of the power station, QO represents the reactive power of the power station in the target area, QO1 represents the new reactive power of the power station, PO1 represents the adjusted active power of the generator, b (QO) Indicates the rate of change of reactive power of a power station;

[0073] The adjusted reactive power of the power station is recorded as QO2;

[0074] When b (QO) ≥ b (PFO) When QO2=PO1 / [(1+b (PFO) )×FPO];

[0075] If b (QO) (PFO) QO2 = QO × (1 + b (QO) );

[0076] Where, FPO represents power factor; b (PFO) Indicates the rate of change of the power factor of the power station, b (QO) It indicates the rate of change of reactive power of power station; QO2 indicates the reactive power of power station after adjustment;

[0077] The active power of the power station is adjusted to PO1 and the reactive power is adjusted to QO2.

[0078] As a preferred solution of the multi-regional interconnected power grid resource scheduling method based on big data described in the present invention, the power adjustment also includes obtaining the current active power of the power station and recording it as PP o , reactive power is recorded as QQ o ;

[0079] Calculate the current power factor of the power plant, denoted as PF P The formula is:

[0080] PFp=PP o / QQ o

[0081] The formula for obtaining the amount of electricity to be allocated corresponding to the 1st to mnth municipal districts is as follows:

[0082] (pj1,pj2...pj mn )

[0083] Among them, pj1 represents the amount of electricity to be allocated corresponding to the first municipal district; pj2 represents the amount of electricity to be allocated corresponding to the second municipal district; pj mn Indicates the amount of electricity to be allocated corresponding to the mn-th municipal district;

[0084] Calculate the transformation ratio of the receiving transformer in the first municipal district, denoted as bb1, and the formula is:

[0085] bb1=(PP o ×tt) / pj1

[0086] Calculate the transformation ratio of the receiving transformer in the second municipal district, denoted as bb1, and the formula is: ​

[0087] bb2=(PP o ×tt) / pj2

[0088] Calculate the transformation ratio of the receiving transformer in the mnth municipal district, denoted as bb mn The formula is:

[0089] bb mn =(PP o ×tt) / pj mn

[0090] Among them, PP o represents the current active power of the power station, tt represents the power generation time of the power station in month m, bb mn represents the transformation ratio of the receiving transformer in the mnth municipal district; bb2 represents the transformation ratio of the receiving transformer in the second municipal district, bb1 represents the transformation ratio of the receiving transformer in the first municipal district; pj1 represents the amount of electricity to be distributed corresponding to the first municipal district; pj2 represents the amount of electricity to be distributed corresponding to the second municipal district; pj mn Indicates the amount of electricity to be allocated corresponding to the mn-th municipal district;

[0091] Obtain the substation parameters and the connection relationship of the substations of the mnth power supply line to obtain the substation data mn; adjust the active power and reactive power of each substation on the mnth power supply line according to the substation data mn to complete the power distribution to the mnth municipal district;

[0092] Substation parameter representation includes active power and reactive power of the substation;

[0093] Taking the receiving transformer as the starting point and the starting transformer as the end point, obtain the number of all substations on the first power supply line, recorded as ts;

[0094] Receiving transformer to the first substation;

[0095] Summarize the active power of all substations connected to the first transformer and adjust the power of the first substation;

[0096] The active power of the first substation is obtained and recorded as P1, and the reactive power is recorded as Q1;

[0097] The load power of the first substation is denoted as pl1;

[0098] Determine whether the first substation is connected only to the receiving transformer of the first municipal district and determine the value of pl1;

[0099] When the first substation is connected only to the receiving transformer of the first municipal district, pl1 = (pj1 / tt);

[0100] When the first substation is connected to multiple receiving transformers in municipal districts, the receiving transformers in different municipal districts connected to the first substation are recorded as type A transformers as tr1, and the value of pl1 is calculated;

[0101] Calculate the power factor of the first substation, denoted as pf1, and the formula is:

[0102] pf1=P1 / Q1

[0103] Compare the size of P1 and pl1 to determine whether the active power and reactive power of the first substation are adjusted;

[0104] When P1≥pl1, the active power and reactive power of the first substation are not adjusted;

[0105] When P1<pl1, adjust the active power and reactive power of the first substation;

[0106] The adjusted active power of the first substation is recorded as tP1, and the adjusted reactive power is recorded as tQ1. The formula is expressed as follows:

[0107] tP1=P1×[1+(P1 / pl1)]

[0108] tQ1=tP1 / pf1

[0109] Among them, tP1 represents the adjusted active power of the first substation, tQ1 represents the adjusted reactive power, P1 represents the active power of the first substation, pl1 represents the load power of the first substation, pf1 represents the power factor of the first substation, Q1 represents the reactive power of the first substation, and pj1 represents the amount of electricity to be allocated corresponding to the first municipal district;

[0110] The active power of the first substation is used as the load power of the second substation, and the power of the second substation is adjusted until the transformation ratio of the starting transformer is adjusted;

[0111] The substation connected to the starting transformer is taken as the target substation, and the number of target substations is obtained, which is recorded as ta;

[0112] The formula for obtaining the active power of each target substation is expressed as:

[0113] (Pa1,Pa2...Pa ta )

[0114] Where, Pa1 represents the active power of the first target substation; Pa2 represents the active power of the second target substation; Pa ta represents the active power of the ta-th target substation;

[0115] Get the current active power PP of the power stationo , reactive power QQ o and the power factor PF of the current power station p ;

[0116]

[0117] Where aPa represents the total power of the starting substation; Pa1 represents the active power of the first target substation; Pa2 represents the active power of the second target substation; Pa ta represents the active power of the ta-th target substation;

[0118] When aPa<PP o When ch=PP o / aPa;

[0119] If aPa≥PP o When the effective power of the power station is adjusted to aPa, the active power and reactive power of the generator are adjusted repeatedly, and the reactive power of the power station is adjusted twice to obtain the reactive power of the power station after the secondary adjustment, which is recorded as QO t ;

[0120] Adjust the effective power of the power station to aPa and the ineffective power to QO t , repeat the power distribution to the first municipal district until the total power of the starting substation is less than the effective power of the power station.

[0121] As a preferred solution of the multi-regional interconnected power grid resource dispatching method based on big data described in the present invention, wherein: adjusting the power distribution of each municipal area includes obtaining the transformation ratio of each class A transformer to obtain a formula expressed as:

[0122]

[0123] Among them, cr (1,1) Indicates the transformation ratio of the first class A transformer; cr (1,2) Indicates the transformation ratio of the second class a transformer; cr (1,tr1) It represents the transformation ratio of the Tr1th class A transformer;

[0124] Calculate the combined transformation ratio of all class A transformers, denoted as zb1, and the formula is:

[0125]

[0126] Among them, cr (1,j) represents the transformation ratio of the jth class A transformer, and the value range of j is 1 to tr1;

[0127] The formula for obtaining the amount of electricity to be distributed in the municipal area corresponding to each Class A transformer is expressed as:

[0128]

[0129] Among them, pja1 represents the amount of electricity to be distributed in the city area corresponding to the first type a transformer; pja2 represents the amount of electricity to be distributed in the city area corresponding to the second type a transformer; pja tr1 It indicates the amount of electricity to be distributed in the municipal area corresponding to the tr1th type a transformer;

[0130] The formula for calculating pl1 is expressed as:

[0131]

[0132] Among them, pja p It indicates the amount of electricity to be distributed in the city area corresponding to the p-th type A transformer, and the value range of p is 1 to tr1;

[0133] The active power of the second substation is obtained and recorded as P2, and the reactive power is recorded as Q2;

[0134] The voltage value of the second substation is obtained and recorded as Vi2, the phase angle of the current voltage is recorded as vθ2; the current value is recorded as Ii2;

[0135] The voltage assignment of the first substation is recorded as Vi1, the phase angle of the current voltage is recorded as vθ1; the current assignment is recorded as Ii1;

[0136] Calculate the equivalent admittance from the first substation to the second substation. The formula is:

[0137]

[0138] Among them, Y (1-2) represents the equivalent admittance from the first substation to the second substation, V i1 Represents the voltage of the first substation, V i2 Indicates the voltage of the second substation, I i1 Represents the current of the first substation, I i2 represents the current of the second substation, vθ1 represents the phase angle of the voltage of the first substation, and vθ2 represents the phase angle of the voltage of the second substation;

[0139] Calculate the forward active power transmitted from the second substation to the first substation, denoted as PP (2-1) The formula is expressed as:

[0140] PP (2-1) =(Vi2×Vi1)×Y (1-2) ×cos(νθ2-νθ1)

[0141] Among them, V i1Represents the voltage of the first substation, V i2 Indicates the voltage of the second substation, I i1 Represents the current of the first substation, I i2 represents the current of the second substation, vθ1 represents the phase angle of the voltage of the first substation, and vθ2 represents the phase angle of the voltage of the second substation; Y (1-2) represents the equivalent admittance from the first substation to the second substation, PP (2-1) It represents the positive active power transmitted from the second substation to the first substation;

[0142] When PP (2-1) When ≥P2, the active power and reactive power of the second substation are not adjusted;

[0143] When PP (2-1) When <P2, adjust the active power and reactive power of the second substation and calculate the power factor of the second substation;

[0144] Calculate the power factor of the second substation, denoted as pf2, and the formula is:

[0145] pf2=P2 / Q2

[0146] The adjusted active power of the second substation is recorded as tp2, and the adjusted reactive power is recorded as tQ2. The formula is expressed as follows:

[0147] tp2=P2×[1+(PP (2-1) / P2)]

[0148] tQ2=tp2 / pf2

[0149] Among them, pf2 represents the power factor of the second substation, P2 represents the active power of the substation before adjustment, PP (2-1) It represents the forward active power transmitted from the second substation to the first substation, tQ2 represents the adjusted reactive power of the second substation, and tp2 represents the adjusted active power of the second substation.

[0150] A multi-region interconnected power grid resource dispatching system based on big data, wherein:

[0151] The data acquisition module acquires the number of urban districts in the target area; acquires the current power consumption corresponding to each urban district, and obtains power consumption data;

[0152] The data analysis module obtains the historical power consumption of each city district and calculates the amount of power to be allocated in each city district based on the power consumption data; obtains the active power and reactive power of the power station in the target area, and adjusts the power of the power station based on the amount of power to be allocated in each city district, and enters the power dispatch module;

[0153] The power dispatching module obtains the substation parameters and connection relationship of the power supply line corresponding to each municipal district, and obtains the substation data; according to the amount of electricity to be allocated in each municipal district, adjusts the starting transformer ratio and the receiving transformer ratio on each power supply line, and then calculates the power flow for each substation according to the substation data, adjusts the active power and reactive power of the substation, and allocates power to each municipal district;

[0154] The continuous scheduling module updates power consumption data in real time and adjusts the power consumption of each municipal district.

[0155] A computer device comprises: a memory and a processor; the memory stores a computer program, wherein the processor implements the steps of any one of the methods of the present invention when executing the computer program.

[0156] A computer-readable storage medium stores a computer program, wherein the computer program, when executed by a processor, implements the steps of any one of the methods of the present invention.

[0157] Beneficial effects of the present invention: The multi-regional interconnected power grid resource dispatching method based on big data provided by the present invention can dynamically optimize power distribution by acquiring the power data of the municipal area, the substation parameters of the power supply line and the power adjustment mechanism in real time, thereby effectively improving the operating efficiency of the power system. The model based on the power consumption change rate and power factor adjustment is adopted to achieve accurate regulation of the power grid load, ensure the efficient use of power resources and power supply stability, and optimize the power adjustment process by using the state transfer matrix and matrix operation, which is helpful to achieve real-time adaptive adjustment in complex power systems, reduce energy waste, reduce energy consumption, and improve the overall dispatching capability and reliability of the system, and has strong industrial application value. BRIEF DESCRIPTION OF THE DRAWINGS

[0158] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0159] Figure 1 An overall flow chart of a multi-region interconnected power grid resource scheduling method based on big data provided by the first embodiment of the present invention;

[0160] Figure 2 A system schematic diagram of a multi-region interconnected power grid resource scheduling method based on big data provided by a second embodiment of the present invention;

[0161] Figure 3A schematic diagram of the first power supply line of a multi-regional interconnected power grid resource scheduling method based on big data provided in the first embodiment of the present invention. DETAILED DESCRIPTION

[0162] In order to make the above-mentioned purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are described in detail below in conjunction with the drawings of the specification. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in the art without creative work should fall within the scope of protection of the present invention.

[0163] Example 1, reference Figure 1 , is an embodiment of the present invention, and provides a multi-region interconnected power grid resource scheduling method based on big data, comprising:

[0164] S1: Obtain the number of urban districts and the power consumption data of each urban district, and calculate the power consumption change rate.

[0165] The obtaining of the number of city districts and the power consumption data of each city district includes obtaining the specific number of city districts and regional division information from the power management system.

[0166] Real-time and historical electricity consumption data of each district are collected from electricity meters, smart grid monitoring equipment and power metering systems in each municipal district.

[0167] Get the current month as m and the current day as d.

[0168] Get the total number of days in month m, recorded as dt.

[0169] The number of all urban districts in the target area is denoted as mn.

[0170] Get the current power consumption corresponding to each city district. The formula is:

[0171] (oj1,oj2...oj mn )

[0172] Among them, oj1 represents the current power consumption of the first district; oj2 represents the current power consumption of the second district; oj mn Indicates the current power consumption of the mnth municipal district.

[0173] Obtain the electricity consumption of each municipal district in month m in the past three years as the historical electricity consumption.

[0174] The historical electricity consumption of the first municipal district is expressed as ahj (1) 、bhj (1) and chj (1) .

[0175] Among them, ahj (1) represents the electricity consumption of the first municipal district in the third year in month m, bhj (1) represents the power consumption of the first municipal district in the second year in month m, chj (1) It represents the electricity consumption in month m of the first year in the first urban district.

[0176] The historical electricity consumption of the second municipal district is ahj (2) 、bhj (2) and chj (2) .

[0177] Among them, ahj (2) represents the electricity consumption of the second urban district in the third year in month m, bhj (2) represents the power consumption of the second district in the second year in month m, chj (2) It represents the electricity consumption in month m of the first year in the second urban district.

[0178] The historical electricity consumption of the mnth municipal district is ahj (mn) 、bhj (mn) and chj (mn) .

[0179] Among them, ahj (mn) represents the power consumption of the mnth district in the third year in month m, bhj (mn) represents the power consumption of the mnth district in the second year in month m, chj (mn) It represents the electricity consumption in month m of the first year of the mn-th urban district.

[0180] According to the historical power consumption, the formula for calculating the power consumption change rate of each municipal district is expressed as:

[0181]

[0182] Among them, ub (1) represents the electricity consumption change rate of the first city district; ub (2) represents the electricity consumption change rate of the second urban district; ub (mn) represents the electricity consumption change rate of the mnth city district; E i,t represents the electricity consumption of city district i at time t, E i,t-1 It represents the electricity consumption of urban district i in the previous time period.

[0183] Get the total power consumption of each city district last month, and use the reference power formula as follows:

[0184] (hj1,hj2...hj mn )

[0185] Where hj1 represents the total power consumption of the first city district last month; hj2 represents the total power consumption of the second city district last month; hj mn Represents the total electricity consumption of the mnth urban district last month.

[0186] Calculate the amount of electricity to be distributed in the first municipal district. The formula is:

[0187]

[0188] Calculate the amount of electricity to be distributed in the second municipal district, the formula is:

[0189]

[0190] The formula for calculating the amount of electricity to be distributed in the mnth municipal district is:

[0191]

[0192] Among them, pj1, pj2, …, pj mn represents the amount of electricity to be allocated in the jth municipal district, j represents the number of the municipal district, mn represents the total number of municipal districts, ub (1) represents the electricity consumption change rate of the first city district; ub (2) represents the electricity consumption change rate of the second urban district; ub (mn) represents the rate of change of electricity consumption in the mnth district; dt represents the time interval, oj1 represents the change of electricity consumption in the first district; oj2 represents the change of electricity consumption in the second district; oj mn represents the change in electricity consumption in the mnth municipal district; apj represents the sum of the electricity to be allocated in all municipal districts;

[0193] The active power of the power station in the target area is obtained and recorded as PO, and the reactive power is recorded as QO.

[0194] The formula for calculating the power factor PFO of a power station is:

[0195] PFO=PO / QO

[0196] Calculate the total power generation of the power station in month m, denoted as ZJ formula:

[0197] ZJ=PO×tt

[0198] Among them, tt represents the power generation time of the power station in month m; PO represents the active power of the power station in the target area, ZJ represents the total power generation of the power station in month m; QO represents the reactive power of the power station in the target area, and PFO represents the power factor of the power station.

[0199] When ZJ≤apj, the amount of electricity generated by the power station is low, the power of the power station is increased, and the active power and reactive power of the generator are adjusted.

[0200] When ZJ>[apj×(1+α)], the power station generates sufficient electricity and the power of the power station does not need to be adjusted.

[0201] Wherein, α represents the power loss rate.

[0202] Furthermore, by analyzing the power consumption of each district in the target area, the power consumption change rate of each district can be dynamically calculated, and based on this, the power to be allocated in each district can be predicted, thereby realizing the precise dispatch of the smart grid. By obtaining real-time power consumption data, historical data, and the power consumption change trend of each district, the system can automatically adjust the power output of the generator according to the power demand, power consumption change, and power generation capacity of different districts, ensuring that the power generation capacity is increased when the power is tight, and no unnecessary power adjustment is made when the power is sufficient, thereby improving the stability and power supply efficiency of the power grid. At the same time, by calculating the power factor and the total power generation, the operation mode of the power grid can be further optimized, the power loss can be reduced, and the reasonable allocation and efficient use of power resources can be ensured.

[0203] S2: Obtain the substation parameters of the power supply line and adjust the power.

[0204] According to ahj (1) 、bhj (1) and chj (1) .

[0205] Among them, ahj (1) represents the electricity consumption of the first municipal district in the third year in month m, bhj (1) represents the power consumption of the first municipal district in the second year in month m, chj (1) It represents the electricity consumption in month m of the first year in the first urban district.

[0206] Calculate the rate of change of power consumption in the first city district, denoted as ub (1) .

[0207] Calculate the rate of change in electricity consumption from the third year to the second year in the first municipal district, recorded as up (3-2) The formula is:

[0208] up (3-2) =[bhj (1) / chj (1) ]-1

[0209] Calculate the rate of change in electricity consumption from the second year to the first year in the first city district, recorded as up (2-1) The formula is:

[0210] up(2-1) =[ahj (1) / bhj (1) ]-1

[0211] According to up (3-2) and up (2-1) , construct the state transfer matrix, and use the NumPy library function to iterate the state transfer matrix to obtain the steady-state transfer matrix B1.

[0212] A 2×1 order receiving matrix {0, 1} is defined, denoted as matrix B2.

[0213] Use matrix multiplication to multiply matrix B1 by matrix B2 to obtain matrix B3.

[0214] Calculate the average value of each parameter in matrix B3 as the power consumption change rate of the first district, and get ub (1) .

[0215] Calculate the power consumption change rate of the 2nd to mnth districts, and get ub (2) ) (mn) .

[0216] Calculate the adjusted active power of the generator, denoted as PO1, and the formula is expressed as:

[0217] PO1=PO+[(apj-ZJ) / tt]

[0218] [(apj-ZJ) / tt]≥(PO×α)

[0219] Among them, PO1 represents the adjusted active power of the generator, up (3-2) Indicates the rate of change of electricity consumption from the third year to the second year in the first municipal district; bhj (1) Indicates the power consumption of the first municipal district in the second or third year; chj (1) Indicates the power consumption of the first municipal district in the third or second year; up (2-1) Indicates the rate of change of electricity consumption from the second year to the first year in the first city district, ahj (1) It represents the power consumption of the first urban district in the first year, B1 represents the steady-state transfer matrix; B2 represents the acceptance matrix, B3 represents the matrix multiplication sum, PO represents the active power of the power station in the target area; apj represents the total amount of power to be allocated, ZJ represents the total power generation of the power station in month m, and tt represents the power generation time of the power station in month m.

[0220] Assuming the reactive power of the power station remains unchanged, calculate the new power factor of the power station, denoted as PFO1, and the formula is expressed as:

[0221] PFO1=PO1 / QO

[0222] Assuming the power factor of the power station remains unchanged, calculate the new reactive power of the power station, denoted as QO1, and the formula is expressed as:

[0223] QO1=PO1 / FPO

[0224] Among them, FPO represents the power factor; PO1 represents the adjusted active power of the generator, QO1 represents the new reactive power of the power station, PFO1 represents the new power factor of the power station, and QO represents the reactive power of the power station in the target area.

[0225] Calculate the rate of change of reactive power of the power station, denoted as b (QO) The formula is:

[0226] b (QO) =(QO1-QO) / QO

[0227] Calculate the rate of change of the power factor of the power station, denoted as b (PFO) The formula is:

[0228] b (PFO) =(PFO1-PFO) / PFO

[0229] Among them, b (PFO) represents the rate of change of the power factor of the power station, PFO1 represents the new power factor of the power station, PFO represents the power factor of the power station, QO represents the reactive power of the power station in the target area, QO1 represents the new reactive power of the power station, PO1 represents the adjusted active power of the generator, b (QO) Indicates the rate of change of reactive power of a power station.

[0230] The adjusted reactive power of the power station is denoted as QO2.

[0231] When b (QO) ≥ b (PFO) When QO2=PO1 / [(1+b (PFO) )×FPO];

[0232] If b (QO) (PFO) QO2 = QO × (1 + b (QO) );

[0233] Where, FPO represents power factor; b (PFO) Indicates the rate of change of the power factor of the power station, b (QO) It indicates the rate of change of reactive power of power station; QO2 indicates the reactive power of power station after adjustment.

[0234] The active power of the power station is adjusted to PO1 and the reactive power is adjusted to QO2.

[0235] ​Power adjustment also includes obtaining the current active power of the power station and recording it as PP o , reactive power is recorded as QQ o .

[0236] Calculate the current power factor of the power plant, denoted as PF P The formula is:

[0237] PFp=PP o / QQ o

[0238] The formula for obtaining the amount of electricity to be allocated corresponding to the 1st to mnth municipal districts is as follows:

[0239] (pj1,pj2...pj mn )

[0240] Among them, pj1 represents the amount of electricity to be allocated corresponding to the first municipal district; pj2 represents the amount of electricity to be allocated corresponding to the second municipal district; pj mn Indicates the amount of electricity to be allocated corresponding to the mn-th urban district.

[0241] Calculate the transformation ratio of the receiving transformer in the first municipal district, denoted as bb1, and the formula is:

[0242] bb1=(PP o ×tt) / pj1

[0243] Calculate the transformation ratio of the receiving transformer in the second municipal district, denoted as bb1, and the formula is:

[0244] bb2=(PP o ×tt) / pj2

[0245] Calculate the transformation ratio of the receiving transformer in the mnth municipal district, denoted as bb mn The formula is:

[0246] bb mn =(PP o ×tt) / pj mn

[0247] Among them, PP o represents the current active power of the power station, tt represents the power generation time of the power station in month m, bb mn represents the transformation ratio of the receiving transformer in the mnth municipal district; bb2 represents the transformation ratio of the receiving transformer in the second municipal district, bb1 represents the transformation ratio of the receiving transformer in the first municipal district; pj1 represents the amount of electricity to be distributed corresponding to the first municipal district; pj2 represents the amount of electricity to be distributed corresponding to the second municipal district; pj mn Indicates the amount of electricity to be allocated corresponding to the mn-th urban district.

[0248] The substation parameters and the connection relationship of the substations of the mnth power supply line are obtained to obtain the substation data mn; the active power and reactive power of each substation on the mnth power supply line are adjusted according to the substation data mn to complete the power distribution to the mnth municipal district.

[0249] The substation parameters include the active power and reactive power of the substation.

[0250] like Figure 3 As shown, taking the receiving transformer as the starting point and the starting transformer as the end point, the number of all substations on the first power supply line is obtained, which is recorded as ts.

[0251] Receiving transformer to the 1st substation.

[0252] Summarize the active power of all substations connected to the first transformer and adjust the power of the first substation.

[0253] The active power of the first substation is obtained and recorded as P1, and the reactive power is recorded as Q1.

[0254] The load power of the first substation is denoted as pl1.

[0255] Determine whether the first substation is only connected to the receiving transformer of the first municipal district and determine the value of pl1.

[0256] When the first substation is connected only to the receiving transformer of the first municipal district, pl1 = (pj1 / tt).

[0257] When the first substation is connected to receiving transformers in multiple municipalities, the receiving transformers in different municipalities connected to the first substation are recorded as Class A transformers as tr1, and the value of pl1 is calculated.

[0258] Calculate the power factor of the first substation, denoted as pf1, and the formula is:

[0259] pf1=P1 / Q1

[0260] Compare the sizes of P1 and pl1 to determine whether the active power and reactive power of the first substation are adjusted.

[0261] When P1≥pl1, the active power and reactive power of the first substation are not adjusted.

[0262] When P1<pl1, adjust the active power and reactive power of the first substation.

[0263] The adjusted active power of the first substation is recorded as tP1, and the adjusted reactive power is recorded as tQ1. The formula is expressed as follows:

[0264] tP1=P1×[1+(P1 / pl1)]

[0265] tQ1=tP1 / pf1

[0266] Among them, tP1 represents the adjusted active power of the first substation, tQ1 represents the adjusted reactive power, P1 represents the active power of the first substation, pl1 represents the load power of the first substation, pf1 represents the power factor of the first substation, Q1 represents the reactive power of the first substation, and pj1 represents the amount of electricity to be allocated corresponding to the first municipal district.

[0267] The active power of the first substation is used as the load power of the second substation, and the power of the second substation is adjusted until the transformation ratio of the starting transformer is adjusted.

[0268] The substation connected to the starting transformer is taken as the target substation, and the number of target substations is obtained, which is recorded as ta.

[0269] The formula for obtaining the active power of each target substation is expressed as:

[0270] (Pa1,Pa2...Pa ta )

[0271] Wherein, Pa1 represents the active power of the first target substation; Pa2 represents the active power of the second target substation; Pa ta Represents the active power of the ta-th target substation.

[0272] Get the current active power PP of the power station o , reactive power QQ o and the power factor PF of the current power station p .

[0273]

[0274] Where aPa represents the total power of the starting substation; Pa1 represents the active power of the first target substation; Pa2 represents the active power of the second target substation; Pa ta Represents the active power of the ta-th target substation.

[0275] When aPa<PP o When ch=PP o / aPa.

[0276] If aPa≥PP o When the effective power of the power station is adjusted to aPa, the active power and reactive power of the generator are adjusted repeatedly, and the reactive power of the power station is adjusted twice to obtain the reactive power of the power station after the secondary adjustment, which is recorded as QO t .

[0277] Adjust the effective power of the power station to aPa and the ineffective power to QO t , repeat the power distribution to the first municipal district until the total power of the starting substation is less than the effective power of the power station.

[0278] Furthermore, by gradually adjusting the active power and reactive power of the power station, it is ensured that the power distribution of the power station can meet the power demand of each municipal district, while maintaining the stability of the system and the rationality of the power factor. By dynamically adjusting the power of the generator, it is possible to effectively respond to changes in power consumption in each municipal district, optimize the allocation of power resources, and improve the operating efficiency and reliability of the power supply system. In addition, by adjusting the power factor and load power of the substation, power loss can be reduced, power supply quality can be improved, and power coordination between substations can be ensured, so that the entire power network operates in a steady state and avoids overload and imbalance problems. This process is particularly important for complex power systems, which can optimize power distribution in actual operation and avoid insufficient or wasteful power supply.

[0279] S3: Update power consumption data in real time and adjust power distribution in each city district.

[0280] Adjusting the power distribution of each municipal area includes obtaining the transformation ratio of each class A transformer and obtaining the formula expressed as:

[0281]

[0282] Among them, cr (1,1) Indicates the transformation ratio of the first class A transformer; cr (1,2) Indicates the transformation ratio of the second class a transformer; cr (1,tr1) Represents the transformation ratio of the Tr1th class A transformer.

[0283] Calculate the combined transformation ratio of all class A transformers, denoted as zb1, and the formula is:

[0284]

[0285] Among them, cr (1,j) It represents the transformation ratio of the jth class A transformer, and the value range of j is 1 to tr1.

[0286] The formula for obtaining the amount of electricity to be distributed in the municipal area corresponding to each Class A transformer is expressed as:

[0287]

[0288] Among them, pja1 represents the amount of electricity to be distributed in the city area corresponding to the first type a transformer; pja2 represents the amount of electricity to be distributed in the city area corresponding to the second type a transformer; pja tr1It indicates the amount of electricity to be distributed in the municipal area corresponding to the Tr1th A-type transformer.

[0289] The formula for calculating pl1 is expressed as:

[0290]

[0291] Among them, pja p It indicates the amount of electricity to be distributed in the municipal area corresponding to the p-th type A transformer, and the value range of p is 1 to tr1.

[0292] The active power of the second substation is obtained and recorded as P2, and the reactive power is recorded as Q2.

[0293] The voltage assignment of the second substation is recorded as Vi2, the phase angle of the current voltage is recorded as vθ2; the current assignment is recorded as Ii2.

[0294] The value of the voltage at the first substation is recorded as Vi1, the phase angle of the current voltage is recorded as vθ1; the value of the current is recorded as Ii1.

[0295] Calculate the equivalent admittance from the first substation to the second substation. The formula is:

[0296]

[0297] Among them, Y (1-2) represents the equivalent admittance from the first substation to the second substation, V i1 Represents the voltage of the first substation, V i2 Indicates the voltage of the second substation, I i1 Represents the current of the first substation, I i2 represents the current of the second substation, vθ1 represents the phase angle of the voltage of the first substation, and vθ2 represents the phase angle of the voltage of the second substation.

[0298] Calculate the forward active power transmitted from the second substation to the first substation, denoted as PP (2-1) The formula is expressed as:

[0299] PP (2-1) =(Vi2×Vi10×Y (1-2) ×cos(νθ2-νθ1)

[0300] Among them, V i1 Represents the voltage of the first substation, V i2 Indicates the voltage of the second substation, I i1 Represents the current of the first substation, I i2 represents the current of the second substation, vθ1 represents the phase angle of the voltage of the first substation, and vθ2 represents the phase angle of the voltage of the second substation; Y(1-2) represents the equivalent admittance from the first substation to the second substation, PP (2-1) It represents the forward active power transmitted from the second substation to the first substation.

[0301] When PP (2-1) When ≥P2, the active power and reactive power of the second substation are not adjusted.

[0302] When PP (2-1) When <P2, adjust the active power and reactive power of the second substation and calculate the power factor of the second substation.

[0303] Calculate the power factor of the second substation, denoted as pf2, and the formula is:

[0304] pf2=P2 / Q2

[0305] The adjusted active power of the second substation is recorded as tp2, and the adjusted reactive power is recorded as tQ2. The formula is expressed as follows:

[0306] tp2=P2×[1+(PP (2-1) / P2)]

[0307] tQ2=tp2 / pf2

[0308] Among them, pf2 represents the power factor of the second substation, P2 represents the active power of the substation before adjustment, PP (2-1) It represents the forward active power transmitted from the second substation to the first substation, tQ2 represents the adjusted reactive power of the second substation, and tp2 represents the adjusted active power of the second substation.

[0309] Furthermore, by updating and dynamically adjusting the power distribution in real time, the balance and stability of power supply in each city district can be ensured. By accurately calculating the transformer ratio, the amount of power to be distributed, and the power and admittance parameters of the substation, the load distribution of the power grid can be optimized, and the power factor and operating efficiency of the system can be improved. When there is power imbalance in the substation, adjusting its active and reactive power can not only balance the load of the power grid, but also avoid excessive power loss, ensuring the safe and stable operation of the power system.

[0310] Embodiment 2 is an embodiment of the present invention, which provides a multi-region interconnected power grid resource dispatching system based on big data, including:

[0311] The data acquisition module obtains the number of urban districts in the target area; obtains the current power consumption corresponding to each urban district, and obtains power consumption data.

[0312] The data analysis module obtains the historical power consumption of each urban district and calculates the amount of power to be allocated in each urban district based on the power consumption data; obtains the active power and reactive power of the power station in the target area, and adjusts the power of the power station according to the amount of power to be allocated in each urban district, and enters the power dispatching module.

[0313] The power dispatching module obtains the substation parameters and connection relationship of the power supply line corresponding to each municipal district, and obtains the substation data; according to the amount of electricity to be distributed in each municipal district, the starting transformer ratio and the receiving transformer ratio on each power supply line are adjusted, and then the power flow calculation is performed for each substation according to the substation data, and the active power and reactive power of the substation are adjusted to distribute power to each municipal district.

[0314] The continuous scheduling module updates power consumption data in real time and adjusts the power consumption of each municipal district.

[0315] Embodiment 3, an embodiment of the present invention, is different from the first two embodiments in that:

[0316] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc., which can store program code.

[0317] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in conjunction with such instruction execution systems, devices or apparatuses. For the purposes of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in conjunction with such instruction execution systems, devices or apparatuses.

[0318] More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk case (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be a paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, deciphering or, if necessary, processing in another suitable manner, and then stored in a computer memory.

[0319] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0320] In Example 4, real-time power consumption data of each municipal district is collected, and the amount of electricity to be allocated in each municipal district is calculated through historical data. Then, by obtaining the data of the power station in the target area, the power of the power station is adjusted to meet the needs of each municipal district. In actual application, the active power and reactive power of the power station are optimized and adjusted according to the power demand of each municipal district, and the active power and reactive power of the substation are also adjusted accordingly. In order to ensure the stability and reliability of the power grid, the experiment continuously updates the power consumption data and adjusts the power on this basis.

[0321] Obtain power consumption data and historical data for the target area. Calculate the amount of electricity to be allocated for each municipal district based on real-time power consumption and historical data, and determine the power demand for each municipal district. Adjust the power of the power station and the active and reactive power on the power supply lines in each municipal district based on the amount of electricity to be allocated. Update power consumption data in real time, and adjust the amount of electricity based on the new power consumption data to optimize the power distribution structure and improve the utilization rate of electric energy.

[0322] During the experiment, the differences between the traditional power dispatching method and the method of the present invention were compared by using the optimization algorithm, especially the advantages in power adjustment of power stations, power flow calculation of substations and power distribution. The experimental results are shown in Table 1.

[0323] Table 1 Experimental data table

[0324]

[0325] It can be seen from the table that the power grid resource dispatching method of the present invention is generally implemented in multiple municipal districts, and the power distribution efficiency is generally high, up to 96.25% (municipal district D), which is much higher than the traditional power dispatching method. By optimizing the power of the power station and the active power and reactive power of the substation, the power loss is significantly reduced, and the energy loss rate of the municipal district is generally maintained at a low level (2.1% to 3.0%). Especially in municipal district A, although the current power consumption is as high as 1200kW, after power adjustment, the power distribution efficiency can be effectively improved and energy waste can be minimized.

[0326] Compared with traditional power dispatching methods, the present invention demonstrates strong optimization capabilities in processing power consumption data, calculating the amount of power to be distributed, and adjusting the power of generators and substations. Traditional methods usually fail to adjust the power demand of each municipal district in real time, resulting in an unreasonable power distribution structure and high energy losses. The present invention not only improves the power distribution efficiency but also reduces energy losses by updating power consumption data in real time and accurately adjusting the power of each municipal district. In addition, in terms of power station power adjustment, the optimization dispatching method based on big data can ensure that the power demand of each municipal district is responded to in a timely manner through accurate power flow calculation and power adjustment of substations, thereby avoiding the over-generation or under-generation common in traditional dispatching methods.

[0327] It can be seen from the data that the grid resource dispatching method implemented in the present invention has significant advantages over traditional methods, especially in terms of power distribution efficiency and energy loss control, showing innovation and high practicality. Through further optimization, it is expected that the stability and reliability of grid operation can be greatly improved, providing solutions for power dispatching in a wider range.

[0328] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention can be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should be included in the scope of the claims of the present invention.

Claims

1. A multi-regional interconnected power grid resource scheduling method based on big data, characterized in that: include: Obtain the number of municipal districts and the power consumption data of each municipal district, and calculate the power consumption change rate; Obtain substation parameters of the power supply line and adjust power; Update power consumption data in real time and adjust power distribution in each city district.

2. The multi-regional interconnected power grid resource scheduling method based on big data according to claim 1, characterized in that: The obtaining of the number of city districts and the power consumption data of each city district includes obtaining the specific number of city districts and regional division information from the power management system; Real-time and historical electricity consumption data of each district are collected from electricity meters, smart grid monitoring equipment and power metering systems in each municipal district.

3. The multi-regional interconnected power grid resource scheduling method based on big data as claimed in claim 2, characterized in that: Calculating the electricity consumption change rate includes obtaining the current month as m and the current date as d; Get the total number of days in month m, recorded as dt; The number of all municipal districts in the target area, denoted as mn; Get the current power consumption corresponding to each city district. The formula is: (oj1,oj2...oj mn ) Among them, oj1 represents the current power consumption of the first district; oj2 represents the current power consumption of the second district; oj mn represents the current power consumption of the mnth municipal district; Obtain the electricity consumption of each city district in month m in the past three years as the historical electricity consumption; The historical electricity consumption of the first municipal district is expressed as ahj (1) 、bhj (1) and chj (1) ; Among them, ahj (1) represents the electricity consumption of the first municipal district in the third year in month m, bhj (1) represents the power consumption of the first municipal district in the second year in month m, chj (1) It represents the electricity consumption in month m of the first year in the first municipal district; The historical electricity consumption of the second municipal district is ahj (2) 、bhj (2) and chj (2) ; Among them, ahj (2) represents the electricity consumption of the second urban district in the third year in month m, bhj (2) represents the power consumption of the second district in the second year in month m, chj (2) It represents the electricity consumption in month m of the first year of the second municipal district; The historical electricity consumption of the mnth municipal district is ahj (mn) 、bhj (mn) and chj (mn) ; Among them, ahj (mn) represents the power consumption of the mnth district in the third year in month m, bhj (mn) represents the power consumption of the mnth district in the second year in month m, chj (mn) It represents the electricity consumption of the mth district in the first year and month m; According to the historical power consumption, the formula for calculating the power consumption change rate of each municipal district is expressed as: Wherein, ub(1) represents the electricity consumption change rate of the first city district; ub(2) represents the electricity consumption change rate of the second city district; (mn) represents the electricity consumption change rate of the mnth city district; E i,t represents the electricity consumption of city district i at time t, E i,t-1 It represents the electricity consumption of urban district i in the previous time period.

4. The multi-regional interconnected power grid resource scheduling method based on big data as claimed in claim 3, characterized in that: Calculating the electricity consumption change rate also includes obtaining the total electricity consumption of each city district last month, which is expressed as the reference electricity formula: (hj1,hj2...hj mn ) Where hj1 represents the total power consumption of the first city district last month; hj2 represents the total power consumption of the second city district last month; hj mn It represents the total electricity consumption of the mnth municipal district last month; Calculate the amount of electricity to be distributed in the first municipal district. The formula is: Calculate the amount of electricity to be distributed in the second municipal district, the formula is: The formula for calculating the amount of electricity to be distributed in the mnth municipal district is: Among them, pj1, pj2, …, pj mn represents the amount of electricity to be allocated in the jth municipal district, j represents the number of the municipal district, mn represents the total number of municipal districts, ub (1) represents the electricity consumption change rate of the first city district; ub (2) represents the electricity consumption change rate of the second urban district; ub (mn) represents the rate of change of electricity consumption in the mnth district; dt represents the time interval, oj1 represents the change of electricity consumption in the first district; oj2 represents the change of electricity consumption in the second district; oj mn represents the change in electricity consumption in the mnth municipal district; apj represents the sum of the electricity to be allocated in all municipal districts; The active power of the power station in the target area is obtained as PO, and the reactive power is recorded as QO; The formula for calculating the power factor PFO of a power station is: PFO==PO / QO Calculate the total power generation of the power station in month m, denoted as ZJ formula: ZJ=PO×tt Among them, tt represents the power generation time of the power station in month m; PO represents the active power of the power station in the target area, ZJ represents the total power generation of the power station in month m; QO represents the reactive power of the power station in the target area, and PFO represents the power factor of the power station; When ZJ≤apj, the amount of electricity generated by the power station is low, the power of the power station is increased, and the active power and reactive power of the generator are adjusted; When ZJ>[apj×(1+α)], the power generated by the power station is sufficient and the power of the power station is not adjusted; Wherein, α represents the power loss rate.

5. The multi-regional interconnected power grid resource scheduling method based on big data as claimed in claim 4, characterized in that: Power adjustment includes, according to ahj (1) 、bhj (1) and chj (1) ; Among them, ahj (1) represents the electricity consumption of the first municipal district in the third year in month m, bhj (1) represents the power consumption of the first municipal district in the second year in month m, chj (1) It represents the electricity consumption in month m of the first year in the first municipal district; Calculate the rate of change of power consumption in the first city district, denoted as ub (1) ; Calculate the rate of change in electricity consumption from the third year to the second year in the first municipal district, recorded as up (3-2) The formula is: up (3-2) =[bhj (1) / chj (1) ]-1 Calculate the rate of change in electricity consumption from the second year to the first year in the first city district, recorded as up (2-1) The formula is: up (2-1) =[ahj (1) / bhj (1) ]-1 According to up (3-2) and up (2-1) , construct the state transfer matrix, and use the NumPy library function to iterate the state transfer matrix to obtain the steady-state transfer matrix B1; Define a 2×1 order receiving matrix {0, 1}, denoted as matrix B2; Use matrix multiplication to calculate the matrix B1 multiplied by the matrix B2 to obtain the matrix B3; Calculate the average value of each parameter in matrix B3 as the power consumption change rate of the first city district, and get ub (1) ; Calculate the power consumption change rate of the 2nd to mnth districts, and get ub (2) ) (mn) ; The adjusted generator active power is calculated and recorded as PO1. The formula is: PO1=PO+[(apj-ZJ) / tt] [(apj-ZJ) / tt]≥(PO×α) Among them, PO1 represents the adjusted active power of the generator, up (3-2) Indicates the rate of change of electricity consumption from the third year to the second year in the first municipal district; bhj (1) Indicates the power consumption of the first municipal district in the second or third year; chj (1) Indicates the power consumption of the first municipal district in the third or second year; up (2-1) Indicates the rate of change of electricity consumption from the second year to the first year in the first city district, ahj (1) represents the power consumption of the first municipal district in the first year, B1 represents the steady-state transfer matrix; B2 represents the acceptance matrix, B3 represents the matrix multiplication sum, PO represents the active power of the power station in the target area; apj represents the total amount of power to be allocated, ZJ represents the total power generation of the power station in month m, and tt represents the power generation time of the power station in month m; Assuming the reactive power of the power station remains unchanged, calculate the new power factor of the power station, denoted as PFO1, and the formula is expressed as: PFO1=PO1 / QO Assuming the power factor of the power station remains unchanged, calculate the new reactive power of the power station, denoted as QO1, and the formula is expressed as: QO1=PO1 / FPO Among them, FPO represents the power factor; PO1 represents the adjusted active power of the generator, QO1 represents the new reactive power of the power station, PFO1 represents the new power factor of the power station, and QO represents the reactive power of the power station in the target area; Calculate the rate of change of reactive power of the power station, denoted as b (QO) The formula is: b (QO) =(QO1-QO) / QO Calculate the rate of change of the power factor of the power station, denoted as b (PFO) The formula is: b (PFO) =(PFO1-PFO) / PFO Among them, b (PFO) represents the rate of change of the power factor of the power station, PFO1 represents the new power factor of the power station, PFO represents the power factor of the power station, QO represents the reactive power of the power station in the target area, QO1 represents the new reactive power of the power station, PO1 represents the adjusted active power of the generator, b (QO) Indicates the rate of change of reactive power of a power station; The adjusted reactive power of the power station is recorded as QO2; When b (QO) ≥ b (PFO) When QO2=PO1 / [(1+b (PFO) )×FPO]; If b (QO) (PFO) QO2 = QO × (1 + b (QO) );​ Where, FPO represents power factor; b (PFO) Indicates the rate of change of the power factor of the power station, b (QO) It indicates the rate of change of reactive power of power station; QO2 indicates the reactive power of power station after adjustment; The active power of the power station is adjusted to PO1 and the reactive power is adjusted to QO2.

6. The multi-regional interconnected power grid resource scheduling method based on big data according to claim 5, characterized in that: Power adjustment also includes obtaining the current active power of the power station and recording it as PP o , reactive power is recorded as QQ o ; Calculate the current power factor of the power plant, denoted as PF P The formula is: PFp=PP o / QQ o The formula for obtaining the amount of electricity to be allocated corresponding to the 1st to mnth municipal districts is as follows: (pj1,pj2...pj mn ) Among them, pj1 represents the amount of electricity to be allocated corresponding to the first municipal district; pj2 represents the amount of electricity to be allocated corresponding to the second municipal district; pj mn Indicates the amount of electricity to be allocated corresponding to the mn-th municipal district; Calculate the transformation ratio of the receiving transformer in the first municipal district, denoted as bb1, and the formula is: bb1=(PP o ×tt) / pj1 Calculate the transformation ratio of the receiving transformer in the second municipal district, denoted as bb1, and the formula is: bb2=(PP o ×tt) / pj2 Calculate the transformation ratio of the receiving transformer in the mnth municipal district, denoted as bb mn The formula is: bb mn =(PP o ×tt) / pj mn Among them, PP o represents the current active power of the power station, tt represents the power generation time of the power station in month m, bb mn represents the transformation ratio of the receiving transformer in the mnth municipal district; bb2 represents the transformation ratio of the receiving transformer in the second municipal district, bb1 represents the transformation ratio of the receiving transformer in the first municipal district; pj1 represents the amount of electricity to be distributed corresponding to the first municipal district; pj2 represents the amount of electricity to be distributed corresponding to the second municipal district; pj mn Indicates the amount of electricity to be allocated corresponding to the mn-th municipal district; Obtain the substation parameters and the connection relationship of the substations of the mnth power supply line to obtain the substation data mn; adjust the active power and reactive power of each substation on the mnth power supply line according to the substation data mn to complete the power distribution to the mnth municipal district; Substation parameter representation includes active power and reactive power of the substation; Taking the receiving transformer as the starting point and the starting transformer as the end point, obtain the number of all substations on the first power supply line, recorded as ts; Receiving transformer to the first substation; Summarize the active power of all substations connected to the first transformer and adjust the power of the first substation; The active power of the first substation is obtained and recorded as P1, and the reactive power is recorded as Q1; The load power of the first substation is denoted as pl1; Determine whether the first substation is connected only to the receiving transformer of the first municipal district and determine the value of pl1; When the first substation is connected only to the receiving transformer of the first municipal district, pl1 = (pj1 / tt); When the first substation is connected to multiple receiving transformers in municipal districts, the receiving transformers in different municipal districts connected to the first substation are recorded as type A transformers as tr1, and the value of pl1 is calculated; Calculate the power factor of the first substation, denoted as pf1, and the formula is: pf1=P1 / Q1 Compare the size of P1 and pl1 to determine whether the active power and reactive power of the first substation are adjusted; When P1≥pl1, the active power and reactive power of the first substation are not adjusted; When P1<pl1, adjust the active power and reactive power of the first substation; The adjusted active power of the first substation is recorded as tP1, and the adjusted reactive power is recorded as tQ1. The formula is expressed as follows: tP1=P1×[1+(P1 / pl1)] tQ1=tP1 / pf1 Among them, tP1 represents the adjusted active power of the first substation, tQ1 represents the adjusted reactive power, P1 represents the active power of the first substation, pl1 represents the load power of the first substation, pf1 represents the power factor of the first substation, Q1 represents the reactive power of the first substation, and pj1 represents the amount of electricity to be allocated corresponding to the first municipal district; The active power of the first substation is used as the load power of the second substation, and the power of the second substation is adjusted until the transformation ratio of the starting transformer is adjusted; The substation connected to the starting transformer is taken as the target substation, and the number of target substations is obtained, which is recorded as ta; The formula for obtaining the active power of each target substation is expressed as: (Pa1,Pa2...Pa ta ) Where, Pa1 represents the active power of the first target substation; Pa2 represents the active power of the second target substation; Pa ta represents the active power of the ta-th target substation; Get the current active power PP of the power station o , reactive power QQ o and the power factor PF of the current power station p ; Where aPa represents the total power of the starting substation; Pa1 represents the active power of the first target substation; Pa2 represents the active power of the second target substation; Pa ta represents the active power of the ta-th target substation; When aPa<PP o When ch=PP o / aPa; If aPa≥PP o When the effective power of the power station is adjusted to aPa, the active power and reactive power of the generator are adjusted repeatedly, and the reactive power of the power station is adjusted twice to obtain the reactive power of the power station after the secondary adjustment, which is recorded as QO t ; Adjust the effective power of the power station to aPa and the ineffective power to QO t , repeat the power distribution to the first municipal district until the total power of the starting substation is less than the effective power of the power station.

7. The multi-regional interconnected power grid resource scheduling method based on big data according to claim 6, characterized in that: Adjusting the power distribution of each municipal area includes obtaining the transformation ratio of each class A transformer and obtaining the formula expressed as: Among them, cr (1,1) Indicates the transformation ratio of the first class A transformer; cr (1,2) Indicates the transformation ratio of the second class A transformer; cr (1,tr1) It represents the transformation ratio of the Tr1th class A transformer; Calculate the combined transformation ratio of all class A transformers, denoted as zb1, and the formula is: Among them, cr (1,j) represents the transformation ratio of the jth class A transformer, and the value range of j is 1 to tr1; The formula for obtaining the amount of electricity to be distributed in the municipal area corresponding to each Class A transformer is expressed as: Among them, pja1 represents the amount of electricity to be distributed in the city area corresponding to the first type a transformer; pja2 represents the amount of electricity to be distributed in the city area corresponding to the second type a transformer; pja tr1 It indicates the amount of electricity to be distributed in the municipal area corresponding to the tr1th type a transformer; The formula for calculating pl1 is expressed as: Among them, pja p It indicates the amount of electricity to be distributed in the city area corresponding to the p-th type A transformer, and the value range of p is 1 to tr1; The active power of the second substation is obtained and recorded as P2, and the reactive power is recorded as Q2; The voltage assignment of the second substation is recorded as Vi2, the phase angle of the current voltage is recorded as vθ2; the current assignment is recorded as Ii2; The voltage assignment of the first substation is recorded as Vi1, the phase angle of the current voltage is recorded as vθ1; the current assignment is recorded as Ii1; Calculate the equivalent admittance from the first substation to the second substation. The formula is: Among them, Y (1-2) represents the equivalent admittance from the first substation to the second substation, V i1 Represents the voltage of the first substation, V i2 Indicates the voltage of the second substation, I i1 Represents the current of the first substation, I i2 represents the current of the second substation, vθ1 represents the phase angle of the voltage of the first substation, and vθ2 represents the phase angle of the voltage of the second substation; Calculate the forward active power transmitted from the second substation to the first substation, denoted as PP (2-1) The formula is expressed as: PP (2-1) =(Vi2×Vi1)×Y (1-2) ×cos(νθ2-νθ1 Among them, V i1 Represents the voltage of the first substation, V i2 Indicates the voltage of the second substation, I i1 Represents the current of the first substation, I i2 represents the current of the second substation, vθ1 represents the phase angle of the voltage of the first substation, and vθ2 represents the phase angle of the voltage of the second substation; Y (1-2) represents the equivalent admittance from the first substation to the second substation, PP (2-1) It represents the positive active power transmitted from the second substation to the first substation; When PP (2-1) When ≥P2, the active power and reactive power of the second substation are not adjusted; When PP (2-1) When <P2, adjust the active power and reactive power of the second substation and calculate the power factor of the second substation; Calculate the power factor of the second substation, denoted as pf2, and the formula is: pf2=P2 / Q2 The adjusted active power of the second substation is recorded as tp2, and the adjusted reactive power is recorded as tQ2. The formula is expressed as follows: tp2=P2×[1+(PP (2-1) / P2)] tQ2=tp2 / pf2 Among them, pf2 represents the power factor of the second substation, P2 represents the active power of the substation before adjustment, PP (2-1) It represents the forward active power transmitted from the second substation to the first substation, tQ2 represents the adjusted reactive power of the second substation, and tp2 represents the adjusted active power of the second substation.

8. A multi-regional interconnected power grid resource dispatching system based on big data using the method according to any one of claims 1 to 7, characterized in that: The data acquisition module acquires the number of urban districts in the target area; acquires the current power consumption corresponding to each urban district, and obtains power consumption data; The data analysis module obtains the historical power consumption of each city district and calculates the power to be allocated in each city district based on the power consumption data; Obtain the active power and reactive power of the power station in the target area, and adjust the power of the power station according to the amount of electricity to be allocated in each municipal area, and enter the power dispatch module; The power dispatching module obtains the substation parameters and connection relationship of the power supply line corresponding to each municipal district, and obtains the substation data; according to the amount of electricity to be allocated in each municipal district, adjusts the starting transformer ratio and the receiving transformer ratio on each power supply line, and then calculates the power flow for each substation according to the substation data, adjusts the active power and reactive power of the substation, and allocates power to each municipal district; The continuous scheduling module updates power consumption data in real time and adjusts the power consumption of each municipal district.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the multi-regional interconnected power grid resource scheduling method based on big data are implemented in any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the multi-regional interconnected power grid resource scheduling method based on big data described in any one of claims 1 to 7 are implemented.

Images