A cloud computing-based intelligent control system for power grid operation

Abstract

The present invention provides a cloud computing-based intelligent control system for power grid operation, which belongs to the field of power control and solves the problem of low power resource control efficiency. Specifically, the data acquisition module acquires historical power consumption data and estimates the expected power consumption of each municipal district. The data analysis module acquires transmission data and calculates the equivalent transmission parameters corresponding to each municipal district. The power dispatching module adjusts each power supply line and the working state of each substation according to the expected power and the equivalent transmission parameters. The continuous dispatching module continuously acquires the historical power consumption and transmission data of each municipal district and adjusts the power distribution of each municipal district. The present invention adjusts the power consumption state of each municipal district in the target area by acquiring, processing and analyzing the power consumption information and power grid information of the target area.

Description

Technical Field

[0001] The present invention discloses a cloud computing-based intelligent control system for power grid operation, which relates to the field of power control. Background Art

[0002] The existing power grid control system has the following deficiencies:

[0003] System complexity and cost: The construction of existing systems involves multiple links, including sensor deployment, data collection, algorithm development, and system integration. The construction and maintenance costs of the entire system are relatively high. Especially for old power grids, the investment required for intelligent transformation is very huge. In addition, the complexity of intelligent systems also requires a high level of technical expertise and talent support.

[0004] Personnel adaptability and training issues: The operation of intelligent control systems requires professionals with high technical levels; however, traditional power industry workers may lack experience in operating and maintaining intelligent systems, so a lot of training and skills improvement are needed.

[0005] System stability: As electricity demand increases, the complexity and scale of the power grid are also expanding. Existing systems need to process large amounts of real-time data and dispatch and optimize various aspects of the power grid. However, the complexity of the system and the limitations of processing capacity may lead to response delays, which in turn affect the stability and operational efficiency of the power grid. Summary of the Invention

[0006] In view of the shortcomings of the existing technology, the purpose of the present invention is to provide a cloud computing-based intelligent control system for power grid operation, aiming to solve the problem of low efficiency in power resource control.

[0007] To achieve the above objectives, the present invention is implemented through the following technical solutions: A cloud computing-based intelligent control system for power grid operation includes:

[0008] Data acquisition module: used to obtain the number of all urban districts in the current month and the target area; based on the current month, obtain the corresponding historical electricity consumption of each urban district and estimate the expected electricity consumption of each urban district;

[0009] Data analysis module: used to obtain the no-load electromotive force, rated voltage, reactance voltage phase, and current phase of the power supply end, and obtain the received voltage and received power of each municipal district as transmission data; based on the transmission data, the expected power corresponding to the power supply end to each municipal district, as well as the equivalent reactance, equivalent susceptance, equivalent resistance, and equivalent conductance, are calculated to obtain the equivalent transmission parameters;

[0010] Power dispatch module: used to adjust the corresponding receiving substation ratio of each municipal district according to the expected power; obtain the substation connection relationship of the power supply line corresponding to each municipal district, use the receiving substation as the starting point and the power supply end as the end point, calculate the reverse power flow for each municipal district based on the equivalent transmission parameters, adjust the working status of each substation, and allocate power to each municipal district;

[0011] Continuous scheduling module: used to continuously obtain historical electricity consumption and transmission data of each urban district and adjust the power distribution of each urban district.

[0012] Furthermore, the workflow of the data analysis module is as follows:

[0013] Process A1: Get the current month, recorded as m; get the number of districts in the target area, recorded as mn;

[0014] Get the electricity consumption of the 1st to mnth municipal districts in month m in the past three years as historical electricity consumption, and get:

[0015] The electricity consumption of the 1st to mnth municipal districts in the third year: (1,3) ~he (mn,3) ;

[0016] The electricity consumption of the first to mnth municipal districts in the second year: (1,2) ~he (mn,2) ;

[0017] Electricity consumption of the first to mnth municipal districts in the first year: (1,1) ~he (mn,1) ;

[0018] Process A2: Calculate the expected electricity consumption of the first city district, denoted as Ee (1) ;

[0019] Process A21: 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) ; up (3-2) The calculation formula is: up (3-2) =(he (1,2) -he (1,3) ) / he (1,3) ;

[0020] The rate of change in electricity consumption from the second year to the first year is recorded as up (2-1) ; up (2-1) The calculation formula is: up (2-1) =(he (1,1) -he (1,2) ) / he (1,2) ;

[0021] Process A22: Obtain the electricity consumption of the first municipal district in month (m-1), recorded as hhe;

[0022] Get the current date, recorded as d; get the number of days in month m, recorded as dm;

[0023] Get the cumulative electricity consumption of the first municipal district on day m, denoted as ne;

[0024] Process A23: According to up (3-2) 、up (2-1) , hhe and [(dm×ne) / d] construct the state transfer matrix, and use the NumPy library function to iterate and obtain the steady-state transfer matrix, which is recorded as matrix CB (1) ;

[0025] Define a 1×2 order receiving matrix {1 / 2, 1 / 2}, denoted as matrix CB (2) ; Calculate the matrix CB (1) Multiply by matrix CB (2) , we get the matrix CB (3) ; Calculate the matrix CB (3) The average value of each parameter in is used as the expected electricity consumption of the first urban district in the mth month, denoted as Ee (1) ;

[0026] Process A3: Repeat the calculation of the expected electricity consumption Ee of the first municipal district (1) The same process is used to calculate the expected electricity consumption of the 2nd to mnth districts, and Ee (1) ~Ee (mn) , and enter the data analysis module.

[0027] Furthermore, the workflow of the data analysis module is as follows:

[0028] Process B1: Obtain the no-load electromotive force of the target area's municipal power grid power supply terminal, denoted as Eq (o) , rated voltage is recorded as Un (o) , reactance is denoted as Xd (o) ;

[0029] The voltage phase at the power supply end is obtained and recorded as φ (U) , the current phase is denoted as φ (I) ;

[0030] Calculate the power factor angle at the power supply end, recorded as θ; the calculation formula for θ is: θ = φ (U) -φ (I) ;

[0031] Calculate the active power at the power supply end, recorded as PG (o) PG (o) The calculation formula is:

[0032] PG (o) =[Eq (o) ×Un (o) ] / Xd (o) ×cos(θ);

[0033] Process B2: Get the expected electricity consumption corresponding to the 1st, 2nd, and up to the mnth municipal districts: Ee (1) 、Ee (2) ~Ee (mn) ;

[0034] Calculate Ee (1) 、Ee (2) ~Ee (mn) The ratio of , we get: bl (1) , bl (2) ~bl (mn) ;

[0035] Get the receiving voltage corresponding to the 1st, 2nd, and up to the mnth municipal districts, recorded as: Ui (1) 、Ui (2) ~Ui (mn) ;

[0036] Get the receiving power corresponding to the 1st, 2nd, and up to the mnth municipal districts, recorded as: Pj (1) 、Pj (2) ~Pj (mn) ;

[0037] Process B3: Calculate bl (1) ~bl (mn) The sum of is denoted as abl; calculate Pj (1) ~Pj (mn) The sum is denoted as aPj;

[0038] Calculate the expected power of the first district, denoted as Pi (1) ;Pi (1) The calculation formula is:

[0039] Pi (1) =(bl (1) / abl)×aPj;

[0040] Calculate the expected power of the second city district, recorded as Pi (2) ;Pi (2) The calculation formula is:

[0041] Pi (2) =(bl (2) / abl)×aPj;

[0042] Similarly, calculate the expected power of the mnth district, recorded as Pi (mn) ;Pi (mn)The calculation formula is:

[0043] Pi (mn) =(bl (mn) / abl)×aPj;

[0044] Process B4: According to Ui (1) ~Ui (mn) and Pi (1) ~Pi (mn) , calculate the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance corresponding to the power supply end to each municipal district in turn;

[0045] Process B5: Summarize the equivalent reactances corresponding to the 1st to mnth municipal districts: Xi (mn) ~Xi (mn) , equivalent susceptance: Bi (mn) ~Bi (mn) , equivalent resistance: Ri (mn) ~Ri (mn) , equivalent conductance: Gi (mn) ~Gi (mn) , as equivalent transmission parameters, and enter the power dispatch module.

[0046] Furthermore, the specific process of process B4 is as follows:

[0047] Process B41: Power supply to the first municipal district;

[0048] Process B411: Calculate the receiving current of the first municipal district, denoted as Ii (1) ;

[0049] Ii (1) The calculation formula is: Ii (1) =Pi (1) / Ui (1) ;

[0050] Calculate the equivalent reactance from the power supply end to the first municipal district, denoted as Xi (1) ,

[0051] Xi (1) The calculation formula is: Xi (1) =(Un (o) -Ui (1) ) / Ii (1) ;

[0052] Process B412: Calculate the equivalent susceptance from the power supply end to the first municipal district, denoted as Bi (1) ;Bi (1) The calculation formula is: Bi (1) ={[3×(Ii (1) ) 2 ] 1 / 2} / Ui(1) ;

[0053] Calculate the equivalent resistance from the power supply end to the first municipal district, recorded as Ri (1) ;

[0054] Ri (1) The calculation formula is: Ri (1) =[PG(o)×(Ui (1) ) 2 ] / [(Pi (1) ) 2 ]; calculate the equivalent conductance from the power supply end to the first municipal district, denoted as Gi (1) ;

[0055] Gi (1) The calculation formula is: Gi (1) =(Pi (1) ) / [(Ui (1) ) 2 ];

[0056] Process B42: power supply to the second municipal district;

[0057] Process B421: Calculate the receiving current of the second city district, denoted as Ii (2) ;

[0058] Ii (2) The calculation formula is: Ii (2) =Pi (2) / Ui (2) ;

[0059] Calculate the equivalent reactance from the power supply end to the second municipal district, denoted as Xi (2) ,

[0060] Xi (2) The calculation formula is: Xi (2) =(Un (o) -Ui (2) ) / Ii (2) ;

[0061] Process B422: Calculate the equivalent susceptance from the power supply end to the second municipal district, denoted as Bi (2) ;Bi (2) The calculation formula is: Bi (2) ={[3×(Ii (2) ) 2 ] 1 / 2} / Ui (2) ;

[0062] Calculate the equivalent resistance from the power supply end to the second municipal district, recorded as Ri (2) ;

[0063] Ri(2) The calculation formula is: Ri (2) =[PG(o)×(Ui (2) ) 2 ] / [(Pi (2) ) 2 ]; calculate the equivalent conductance from the power supply end to the second municipal district, denoted as Gi (2) ;

[0064] Gi (2) The calculation formula is: Gi (2) =(Pi (2) ) / [(Ui (2) ) 2 ];

[0065] Process B43: Similarly, calculate the equivalent reactance Xi from the power supply end to the mnth municipal district (mn) , the equivalent electrical susceptance is denoted as Bi (mn) , equivalent resistance Ri (mn) , equivalent conductance Gi (mn) .

[0066] Furthermore, the workflow of the power dispatch module is as follows:

[0067] Process C1: Construct the transmission relationship map of each municipal district based on equivalent transmission parameters;

[0068] The first municipal district: {PG (o) ,Un (o)}—{Xi (1) 、Bi (1) 、Ri (1) , Gi (1)}—{Pi (1) , Ui (1)};

[0069] Second municipal district: {PG (o) ,Un (o)}—{Xi (2) 、Bi (2) 、Ri (2) , Gi (2)}—{Pi (2) , Ui (2)};

[0070] And so on, the mnth district: {PG (o) ,Un (o)}—{Xi (mn) 、Bi (mn) 、Ri (mn) , Gi (mn)}—{Pi (mn) , Ui (mn)};

[0071] Process C2: The power lines that transmit electricity from the power supply end to the 1st, 2nd, and up to the mnth municipal districts are denoted as Line 1, Line 2 to Line mn;

[0072] Get the length of line 1, recorded as ll (1) ; Get the length of line 2, recorded as ll (2) ; Similarly, get the length of line mn, recorded as ll (mn) ;

[0073] Count the number of substations on line 1, denoted as tr (1) ; The number of substations on line 2, denoted as tr (2) Similarly, the number of substations on line mn is denoted as tr (mn) ;

[0074] Get the active power PG of the power supply end (o) , get the expected power of the first, second, and mnth districts: Pi (1) 、Pi (2) ~Pi (mn) ;

[0075] Process C3: Adjust the transformation ratio of receiving substations in the 1st to mnth municipal districts;

[0076] Calculate the transformation ratio of the receiving substation in the first municipal district, denoted as b (o-1) ; b (o-1) The calculation formula is: (o-1) =Pi (1) / PG (o) ;

[0077] Calculate the transformation ratio of the receiving substation in the second municipal district, denoted as b (o-2) ; b (o-2) The calculation formula is: (o-2) =Pi (2) / PG (o) ;

[0078] Similarly, the transformation ratio of the receiving substation in the mnth municipal district is denoted as b (o-mn) ; b (o-mn) The calculation formula is: (o-mn) =Pi (mn) / PG (o) ;

[0079] Among them, b (o-1) ~b (o-mn) All are rational numbers (i.e., finite decimals or infinite repeating decimals);

[0080] According to b (o-1) ~b (o-mn)Adjust the transformation ratio of the receiving substations corresponding to the 1st to mnth municipal districts;

[0081] Process C4: Adjust the working status of each substation on line 1 according to the data in processes C1 to C3;

[0082] Process C5: Repeat the same process of adjusting the working status of each substation on line 1 to adjust line 2 to line mn.

[0083] Furthermore, the specific process of process C4 is as follows:

[0084] Process C41: Obtain the actual receiving voltage of the receiving substation and record it as dU (1) , actual power is recorded as dP (1) According to Pi (1) and Ui (1) , calculate the active power Po of the receiving substation (1) and transmission voltage Uo (1) ;

[0085] Process C42: Get the length of the transmission line from the receiving substation to the first substation, denoted as lr (o-1) ;

[0086] According to Po (1) and Uo (1) , calculate the transmission loss from the first substation to the receiving substation, and get the active power Po of the first substation (2) and transmission voltage Uo (2) ;

[0087] Process C43: Get the actual receiving voltage of the first substation and record it as dU (1) , the actual power is recorded as dP (1) ;

[0088] Determine whether the output terminal of the first substation is connected to other substations, and (2) Adjust the working status of the first substation;

[0089] Process C44: Repeat the same process of calculating the output power and transmission voltage of the first substation and adjusting the working state of the first substation to adjust the second to third substations on line 1. (1) The working status of a substation.

[0090] Furthermore, the specific process of the process C41 is as follows:

[0091] Process C411: Calculate the actual receiving current of the receiving substation, recorded as dI (1) ;dI (1) The calculation formula is: dI (1) =dP (1) / dU (1) ;

[0092] Calculate the actual reactance of the receiving substation, denoted as dX (1) ;

[0093] dX (1) The calculation formula is: dX (1) =[dU (1) -Ui (1) ] / dI (1) Among them, Ui (1) Indicates the receiving voltage of the first municipal district;

[0094] Calculate the actual susceptance of the receiving substation, expressed in dB (1) ;

[0095] dB (1) The calculation formula is: dB (1) ={[3×(dI (1) ) 2 ] 1 / 2} / dU (1) ;

[0096] Calculate the actual resistance of the receiving substation, denoted as dR (1) ;

[0097] dR (1) The calculation formula is: dR (1) =[PG (o) ×(dU (1) ) 2 ] / [mn×(dP (1) ) 2 ];

[0098] Calculate the actual conductance of the receiving substation, denoted as dG (1) ;

[0099] dG (1) The calculation formula is: dG (1) =(dP (1) ) / [(dU (1) ) 2 ]

[0100] Process C412: The equivalent admittance of the receiving substation is recorded as dY (1) ;dY (1) The mathematical expression is: dY (1) =dG (1) +(dB (1) ×j); where j represents an imaginary unit;

[0101] The equivalent impedance of the receiving substation is denoted as dZ (1) ;dZ (1)The mathematical expression is: dZ (1) =dR (1) +(dX (1) ×j); where j represents an imaginary unit;

[0102] Process C413: Pi (1) and Ui (1) For the final state of the receiving substation, reverse power flow calculation is performed on the receiving substation to calculate the active power Po of the receiving substation (1) and transmission voltage Uo (1) ;

[0103] Process C4131: Calculate the load circuit current of the receiving substation, recorded as Iu; the calculation formula for Iu is: Iu = Pi (1) / Ui (1) ;

[0104] Calculate the power consumed by the equivalent impedance of the receiving substation, denoted as Pz (1) =Iu×(Iu * )×dZ (1) ; Wherein, * represents conjugation;

[0105] Calculate the transmission voltage Uo at the receiving substation (1) , Uo (1) The calculation formula is: Uo (1) =(Iu×dZ (1) )+Ui (1) ;

[0106] Process C4132: Calculate the power consumed by the equivalent admittance of the receiving substation, denoted as Py (1) =(Uo (1) ) 2 ×(dY (1) * ); where * represents conjugation;

[0107] Calculate the active power Po of the receiving substation (1) , Po (1) The calculation formula is: Po (1) =Pz (1) +Py (1) +Pi (1) .

[0108] Furthermore, the specific process of the process C42 is as follows:

[0109] Process C421: Calculate the allocation coefficient from the first substation to the receiving substation, denoted as β (1) β (1) The calculation formula is: (1) =lr (o-1) / lr;

[0110] Calculate the equivalent admittance dY of the transmission line from the first substation to the receiving substation (2,1) and dY (2,2) ;dY (2,1) and dY (2,2) The mathematical formula is: dY (2,1) =dY (2,2) =[(β (1) ×Gi (1) )+j×(β (1) ×Bi (1) )] / 2; where j represents the imaginary unit;

[0111] Calculate the equivalent impedance dZ of the transmission line from the first substation to the receiving substation (2) ;dZ (2) The mathematical formula is: dZ (2) =(β (1) ×Ri (1) )+j×(β (1) ×Xi (1) );

[0112] Process C422: Calculate the equivalent admittance dY (2,1) The power consumed is recorded as the equivalent admittance Py (2,1) ;Py (2,1) The calculation formula is: Py (2,1) =(Uo (1) ) 2 ×(dY (2,1) * ); where * represents conjugation;

[0113] Calculate the power of node A, denoted as S (1) ;S (1) The calculation formula is: S (1) =Py (2,1) +Po (1) ;

[0114] Process C423: Calculate the transmission voltage Uo (2) , Uo (2) The calculation formula is: Uo (2) =Uo (1) +[(S (1) / Uo (1) ) * ]×dZ (2) ; Wherein, * represents conjugation;

[0115] Calculate the equivalent impedance dZ (2) The power consumed is denoted as Pz (2) ;Pz (2) The calculation formula is: Pz (2) =(S(1) / Uo (2) ) 2 ×dY (2,2) ;

[0116] Calculate the power of node B, denoted as S (2) ;S (2) The calculation formula is: S (2) =Pz (2) +S (1) ;

[0117] Process C424: Calculate the active power Po of the first substation (2) , Po (2) The calculation formula is:

[0118] Po (2) =S (2) +Uo (2) ×[(Uo (2) ×dY (2,2) ) * ].

[0119] Furthermore, the specific process of process C43 is as follows:

[0120] Process C431: The output terminal of the first substation is not connected to other substations; obtain the reactive power of the first substation and record it as Qo (2) , and adjust the working status of the first substation;

[0121] Process C432: The output terminal of the first substation is connected to other substations;

[0122] The other substations connected to the output of the first substation are regarded as Class A substations. The power supply and voltage of each Class A substation are obtained. The length of the transmission line connecting each Class A substation to the first substation is obtained. The same process of calculating the output power and transmission voltage of the first substation is repeated to calculate the new output power aPo of the first substation. (2) , calculate the new output power aUo of the first substation (2) ;

[0123] Po (2) Replaced by aPo (2) , Uo (2) Replaced with aUo (2) , re-execute process C431 and adjust the working status of the first substation.

[0124] Furthermore, the specific process of the process C431 is as follows:

[0125] Process C4311: Calculate the actual received current of the first substation, denoted as dI (2);dI (2) The calculation formula is: dI (2) =dP (2) / dU (2) ;

[0126] The actual reactance of the first transformer is denoted as dX (2) ,

[0127] dX (2) The calculation formula is: dX (2) =[dU (2) ] / dI (2) ;

[0128] The actual susceptance of the first transformer, expressed in dB (2) ;

[0129] dB (2) The calculation formula is: dB (2) =[3×(dI (2) ) 2 ] 1 / 2 / dU (2) ;

[0130] The actual resistance of the first transformer is denoted as dR (2) ;

[0131] dR (2) The calculation formula is: dR (2) =[PG(o)×(dU (2) ) 2 ] / [mn×(dP (2) ) 2 ];

[0132] The actual conductance of the first transformer is denoted as dG (2) ;

[0133] dG (2) The calculation formula is: dG (2) =(dP (2) ) / [(dU (2) ) 2 ];

[0134] Process C4312: Calculate dX (2) The offset coefficient of the equivalent susceptance is denoted as zx (2) ;

[0135] zx (2) The calculation formula is: zx (2) =1-[dX (1) / (β (1) ×Xi (1) )];

[0136] Calculating dB(2) The deviation coefficient from the equivalent reactance is denoted as zb (2) ;

[0137] zb (2) The calculation formula is: (2) =1-[dB (2) / (β (1) ×Bi (1) )];

[0138] Calculating dR (2) The offset coefficient with the equivalent reactance is denoted as zr (2) ;

[0139] zr (2) The calculation formula is: (2) =1-[dR (2) / (β (1) ×Ri (1) )];

[0140] Calculate dG (2) The deviation coefficient from the equivalent reactance is denoted as zg (2) ;

[0141] zg (2) The calculation formula is: (2) =1-[β (1) / (dG (1) ×Gi (1) )];

[0142] Process C4313: Calculate zx (2) +zb (2) +zr (2) +zg (2) The sum of , recorded as zz (2) ;

[0143] Calculate the power factor of the first substation and record it as pQ. The calculation formula of pQ is: pQ=Po (2) / Qo(2 ) ;

[0144] Assume that the active power after adjustment of the first transformer is PP (1) , judge zz (1) Size determination PP (1) Calculation formula:

[0145] Ifzz (1) ≥0, then PP (1) The calculation formula is: =[(2-zz (1) ) 1 / 2 ]×Po (2) ;

[0146] Ifzz (1)<0, then PP (1) The calculation formula is: =[(1-zz (1) ) 1 / 2 ]×Po (2) ;

[0147] Process C4314: Let the reactive power Qo of the first substation be (2) The power factor after the change is recorded as pQa; the calculation formula of pQa is: pQa=PP (1) / Qo (2) ;

[0148] Let the power factor pQ of the first substation remain unchanged, and record the reactive power after the change as Qa; the calculation formula of pQa is: Qa=PP (1) / pQ;

[0149] Calculate Qa relative to Qo (2) The rate of change of Qo (2) ); b(Qo (2) ) is calculated as: b(Qo (2) )=(Qa-Qo (2) ) / Qo (2) ;

[0150] Calculate the rate of change of pQa relative to pQ, recorded as b(pQ); the calculation formula for b(pQ) is: b(pQ) = (pQ - pQa) / pQ;

[0151] Process C4315: Comparison b(Qo (2) ) and the size of b(pQ);

[0152] If b(Qo (2) )≤b(pQ), then the active power of the first substation is adjusted to PP (1) , reactive power is adjusted to Qa;

[0153] If b(Qo (2) )>b(pQ), then the active power of the first substation is adjusted to PP (1) , reactive power remains unchanged;

[0154] Process C4316: Repeat the same process to calculate the active power and transmission voltage of the receiving substation, according to Po (2) and Uo (2) Calculate the active power and transmission voltage at the second substation.

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

[0156] Improve grid dispatching efficiency: This invention can optimize the grid dispatching strategy through algorithms, thereby improving dispatching efficiency while ensuring safe operation. By analyzing historical and real-time data, the intelligent control system can optimize the dispatching of power generation, transmission, and distribution, reduce resource waste, and improve energy utilization.

[0157] Automated management of power grids: This invention can achieve automated management of power grids and reduce human intervention. Through the automated control system, the power grid can autonomously perform status monitoring, scheduling optimization, and fault handling, thereby reducing the error rate of manual operations and minimizing accidents and losses caused by human factors.

[0158] Optimizing the configuration of power resources: The present invention can optimize the configuration of power resources while ensuring the safe operation of the power grid; the system rationally allocates power generation resources according to the actual demand of the power grid load, avoiding certain generator sets from operating in a low-load state for a long time, thereby improving the operating efficiency and economy of the generator sets. BRIEF DESCRIPTION OF THE DRAWINGS

[0159] Other features, objects and advantages of the present invention will become more apparent upon reading the detailed description of non-limiting embodiments with reference to the following drawings:

[0160] Figure 1 Schematic diagram of the system of the present invention;

[0161] Figure 2 This is a schematic diagram of an equivalent model of a receiving substation according to the present invention;

[0162] Figure 3 Schematic diagram of the equivalent model of the power transmission line of the present invention. DETAILED DESCRIPTION

[0163] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0164] See also Figure 1 ,A power grid operation intelligent control system based on cloud computing includes: a data acquisition module, a data analysis module, a power dispatching module, and a continuous dispatching module;

[0165] Data acquisition module: used to obtain the number of all urban districts in the current month and the target area; based on the current month, obtain the corresponding historical electricity consumption of each urban district and estimate the expected electricity consumption of each urban district;

[0166] It should be noted that the “target area” in the present invention refers to the municipal area where the present invention (a cloud computing-based intelligent control system for power grid operation) is used to dispatch power resources;

[0167] Process A: The workflow of the data analysis module is as follows:

[0168] Process A1: Get the current month, recorded as m; get the number of districts in the target area, recorded as mn;

[0169] Obtain the electricity consumption of the first to mnth municipal districts (within the target area) in the past three years in month m as historical electricity consumption, and obtain:

[0170] The electricity consumption of the 1st to mnth municipal districts in the third year: (1,3) ~he (mn,3) ;

[0171] The electricity consumption of the first to mnth municipal districts in the second year: (1,2) ~he (mn,2) ;

[0172] Electricity consumption of the first to mnth municipal districts in the first year: (1,1) ~he (mn,1) ;

[0173] Process A2: Calculate the expected electricity consumption of the first city district, denoted as Ee (1) ;

[0174] Process A21: 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) ; up (3-2) The calculation formula is: up (3-2) =(he (1,2) -he (1,3) ) / he (1,3) ;

[0175] The rate of change in electricity consumption from the second year to the first year is recorded as up (2-1) ; up (2-1) The calculation formula is: up (2-1) =(he (1,1) -he (1,2) ) / he (1,2) ;

[0176] Process A22: Obtain the electricity consumption of the first municipal district in month m-1, recorded as hhe;

[0177] Get the current date, recorded as d (i.e., today is the d-day of month m); get the number of days in month m, recorded as dm;

[0178] Get the cumulative electricity consumption of the first municipal district on day m, denoted as ne;

[0179] Process A23: According to up (3-2) 、up (2-1), hhe and [(dm×ne) / d] construct the state transfer matrix, and use the NumPy library function to iterate and obtain the steady-state transfer matrix, which is recorded as matrix CB (1) ;

[0180] Define a 1×2 order receiving matrix {1 / 2, 1 / 2}, denoted as matrix CB (2) ; Calculate the matrix CB (1) Multiply by matrix CB (2) , we get the matrix CB (3) ; Calculate the matrix CB (3) The average value of each parameter in is used as the expected electricity consumption of the first urban district in the mth month, denoted as Ee (1) ;

[0181] It should be noted that if month m is January, then month m-1 is December of the previous year; for example, if month m is January 2024, then month m-1 is December 2023;

[0182] Process A3: Repeat the calculation of the expected electricity consumption Ee of the first municipal district (1) The same process is used to calculate the expected electricity consumption of the 2nd to mnth districts, and Ee (1) ~Ee (mn) , and enter the data analysis module.

[0183] Data analysis module: used to obtain the no-load electromotive force, rated voltage, reactance voltage phase, and current phase of the power supply end, and obtain the received voltage and received power of each municipal district as transmission data; based on the transmission data, the expected power corresponding to the power supply end to each municipal district, as well as the equivalent reactance, equivalent susceptance, equivalent resistance, and equivalent conductance, are calculated to obtain the equivalent transmission parameters;

[0184] Process B: The workflow of the data analysis module is as follows:

[0185] Process B1: Obtain the no-load electromotive force of the power supply end of the municipal power grid (or national power grid) in the target area and record it as Eq (o) , rated voltage is recorded as Un (o) , reactance is denoted as Xd (o) ;

[0186] The voltage phase at the power supply end is obtained and recorded as φ (U) , the current phase is denoted as φ (I) ;

[0187] Calculate the power factor angle at the power supply end, recorded as θ; the calculation formula for θ is: θ = φ (U) -φ (I) ;

[0188] Calculate the active power at the power supply end, recorded as PG (o) PG (o)The calculation formula is:

[0189] PG (o) =[Eq (o) ×Un (o) ] / Xd (o) ×cos(θ);

[0190] Process B2: Obtain the expected electricity consumption corresponding to the first, second, and mnth municipal districts (within the target area): Ee (1) 、Ee (2) ~Ee (mn) ;

[0191] Calculate Ee (1) 、Ee (2) ~Ee (mn) The ratio of , we get: bl (1) , bl (2) ~bl (mn) ;

[0192] Get the receiving voltage corresponding to the 1st, 2nd, and up to the mnth municipal districts, recorded as: Ui (1) 、Ui (2) ~Ui (mn) ;

[0193] Get the receiving power corresponding to the 1st, 2nd, and up to the mnth municipal districts, recorded as: Pj (1) 、Pj (2) ~Pj (mn) ;

[0194] It should be noted that the "received voltage" in the present invention means: the voltage of the receiving substation, and the "received power" means: the active power of the receiving substation;

[0195] "Receiving substation" means: a substation connected to a municipal area;

[0196] Process B3: Calculate bl (1) ~bl (mn) The sum of is denoted as abl; calculate Pj (1) ~Pj (mn) The sum of is denoted as aPj;

[0197] Calculate the expected power of the first district, denoted as Pi (1) ;Pi (1) The calculation formula is:

[0198] Pi (1) =(bl (1) / abl)×aPj;

[0199] Calculate the expected power of the second city district, recorded as Pi (2) ;Pi(2) The calculation formula is:

[0200] Pi (2) =(bl (2) / abl)×aPj;

[0201] Similarly, calculate the expected power of the mnth district, recorded as Pi (mn) ;Pi (mn) The calculation formula is:

[0202] Pi (mn) =(bl (mn) / abl)×aPj;

[0203] Process B4: According to Ui (1) ~Ui (mn) and Pi (1) ~Pi (mn) , calculate the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance corresponding to the power supply end (target area) to each municipal district in turn;

[0204] Process B41: Power supply to the first municipal district;

[0205] Process B411: Calculate the receiving current of the first municipal district, denoted as Ii (1) ;

[0206] Ii (1) The calculation formula is: Ii (1) =Pi (1) / Ui (1) ;

[0207] Calculate the equivalent reactance from the power supply end to the first municipal district, denoted as Xi (1) ,

[0208] Xi (1) The calculation formula is: Xi (1) =(Un (o) -Ui (1) ) / Ii (1) ;

[0209] Process B412: Calculate the equivalent susceptance from the power supply end to the first municipal district, denoted as Bi (1) ;Bi (1) The calculation formula is: Bi (1) ={[3×(Ii (1) ) 2 ] 1 / 2} / Ui (1) ;

[0210] Calculate the equivalent resistance from the power supply end to the first municipal district, recorded as Ri (1) ;

[0211] Ri (1) The calculation formula is: Ri (1) =[PG(o)×(Ui (1) ) 2 ] / [(Pi (1) ) 2 ]; calculate the equivalent conductance from the power supply end to the first municipal district, denoted as Gi (1) ;

[0212] Gi (1) The calculation formula is: Gi (1) =(Pi (1) ) / [(Ui (1) ) 2 ];

[0213] Process B42: power supply to the second municipal district;

[0214] Process B421: Calculate the receiving current of the second city district, denoted as Ii (2) ;

[0215] Ii (2) The calculation formula is: Ii (2) =Pi (2) / Ui (2) ;

[0216] Calculate the equivalent reactance from the power supply end to the second municipal district, denoted as Xi (2) ,

[0217] Xi (2) The calculation formula is: Xi (2) =(Un (o) -Ui (2) ) / Ii (2) ;

[0218] Process B422: Calculate the equivalent susceptance from the power supply end to the second municipal district, denoted as Bi (2) ;Bi (2) The calculation formula is: Bi (2) ={[3×(Ii (2) ) 2 ] 1 / 2} / Ui (2) ;

[0219] Calculate the equivalent resistance from the power supply end to the second municipal district, recorded as Ri (2) ;

[0220] Ri (2) The calculation formula is: Ri (2) =[PG(o)×(Ui (2) ) 2 ] / [(Pi (2)) 2 ]; calculate the equivalent conductance from the power supply end to the second municipal district, denoted as Gi (2) ;

[0221] Gi (2) The calculation formula is: Gi (2) =(Pi (2) ) / [(Ui (2) ) 2 ];

[0222] Process B43: and so on, the power supply is provided to the mnth municipal district;

[0223] Process B431: Calculate the receiving current of the mnth district, denoted as Ii (mn) ;

[0224] Ii (mn) The calculation formula is: Ii (mn) =Pi (mn) / Ui (mn) ;

[0225] Calculate the equivalent reactance from the power supply end to the mnth municipal district, denoted as Xi (mn) ,

[0226] Xi (mn) The calculation formula is: Xi (mn) =(Un (o) -Ui (mn) ) / Ii (mn) ;

[0227] Process B432: Calculate the equivalent susceptance from the power supply end to the mnth municipal district, denoted as Bi (mn) ;

[0228] Bi (mn) The calculation formula is: Bi (mn) ={[3×(Ii (mn) ) 2 ] 1 / 2} / Ui (mn) ;

[0229] Calculate the equivalent resistance from the power supply end to the mnth municipal district, recorded as Ri (mn) ;

[0230] Ri (mn) The calculation formula is: Ri (mn) =[PG(o)×(Ui (mn) ) 2 ] / [(Pi (mn) ) 2 ];

[0231] Calculate the equivalent conductance from the power supply end to the mnth municipal district, denoted as Gi(mn) ;

[0232] Gi (mn) The calculation formula is: Gi (mn) =(Pi (mn) ) / [(Ui (mn) ) 2 ];

[0233] Process B5: Summarize the equivalent reactances corresponding to the 1st to mnth municipal districts: Xi (mn) ~Xi (mn) , equivalent susceptance: Bi (mn) ~Bi (mn) , equivalent resistance: Ri (mn) ~Ri (mn) , equivalent conductance: Gi (mn) ~Gi (mn) , as equivalent transmission parameters, and enter the power dispatch module.

[0234] Power dispatch module: used to adjust the corresponding receiving substation ratio of each municipal district according to the expected power; obtain the substation connection relationship of the power supply line corresponding to each municipal district, use the receiving substation as the starting point and the power supply end as the end point, calculate the reverse power flow for each municipal district based on the equivalent transmission parameters, adjust the working status of each substation, and allocate power to each municipal district;

[0235] Process C: The workflow of the power dispatch module is as follows:

[0236] Process C1: Construct the transmission relationship map of each municipal district based on equivalent transmission parameters;

[0237] The first municipal district: {PG (o) ,Un (o)}—{Xi (1) 、Bi (1) 、Ri (1) , Gi (1)}—{Pi (1) , Ui (1)};

[0238] Second municipal district: {PG (o) ,Un (o)}—{Xi (2) 、Bi (2) 、Ri (2) , Gi (2)}—{Pi (2) , Ui (2)};

[0239] And so on, the mnth district: {PG (o) ,Un (o)}—{Xi(mn) 、Bi (mn) 、Ri (mn) , Gi (mn)}—{Pi (mn) , Ui (mn)};

[0240] Process C2: The power lines that transmit electricity from the power supply end to the 1st, 2nd, and up to the mnth municipal districts are denoted as Line 1, Line 2 to Line mn;

[0241] Get the length of line 1, recorded as ll (1) ; Get the length of line 2, recorded as ll (2) ; Similarly, get the length of line mn, recorded as ll (mn) ;

[0242] Count the number of substations on line 1, denoted as tr (1) ; The number of substations on line 2, denoted as tr (2) Similarly, the number of substations on line mn is denoted as tr (mn) ;

[0243] Get the active power PG of the power supply end (o) , get the expected power of the first, second, and mnth districts: Pi (1) 、Pi (2) ~Pi (mn) ;

[0244] Process C3: Adjust the transformation ratio of receiving substations in the 1st to mnth municipal districts;

[0245] Calculate the transformation ratio of the receiving substation in the first municipal district, denoted as b (o-1) ; b (o-1) The calculation formula is: (o-1) =Pi (1) / PG (o) ;

[0246] Calculate the transformation ratio of the receiving substation in the second municipal district, denoted as b (o-2) ; b (o-2) The calculation formula is: (o-2) =Pi (2) / PG (o) ;

[0247] Similarly, the transformation ratio of the receiving substation in the mnth municipal district is denoted as b (o-mn) ; b (o-mn) The calculation formula is: (o-mn) =Pi (mn) / PG (o) ;

[0248] Among them, b (o-1)~b (o-mn) All are rational numbers (i.e., finite decimals or infinite repeating decimals);

[0249] According to b (o-1) ~b (o-mn) Adjust the transformation ratio of the receiving substations corresponding to the 1st to mnth municipal districts;

[0250] Process C4: Adjust the working status of each substation on line 1 according to the data in processes C1 to C3;

[0251] Process C41: Get the actual receiving voltage of the receiving substation (on line 1) and record it as dU (1) , the actual power is recorded as dP (1) According to Pi (1) and Ui (1) , calculate the active power Po of the receiving substation (receiving electric energy) (1) and transmission voltage Uo (1) ;

[0252] Process C411: Calculate the actual receiving current of the receiving substation, recorded as dI (1) ;dI (1) The calculation formula is: dI (1) =dP (1) / dU (1) ;

[0253] Calculate the actual reactance of the receiving substation, denoted as dX (1) ;

[0254] dX (1) The calculation formula is: dX (1) =[dU (1) -Ui (1) ] / dI (1) Among them, Ui (1) Indicates the receiving voltage of the first municipal district;

[0255] Calculate the actual susceptance of the receiving substation, expressed in dB (1) ;

[0256] dB (1) The calculation formula is: dB (1) ={[3×(dI (1) ) 2 ] 1 / 2} / dU (1) ;

[0257] Calculate the actual resistance of the receiving substation, denoted as dR (1) ;

[0258] dR (1) The calculation formula is: dR (1)=[PG (o) ×(dU (1) ) 2 ] / [mn×(dP (1) ) 2 ];

[0259] Calculate the actual conductance of the receiving substation, denoted as dG (1) ;

[0260] dG (1) The calculation formula is: dG (1) =(dP (1) ) / [(dU (1) ) 2 ]

[0261] Process C412: The equivalent admittance of the receiving substation is recorded as dY (1) ;dY (1) The mathematical expression is: dY (1) =dG (1) +(dB (1) ×j); where j represents an imaginary unit;

[0262] The equivalent impedance of the receiving substation is denoted as dZ (1) ;dZ (1) The mathematical expression is: dZ (1) =dR (1) +(dX (1) ×j); where j represents an imaginary unit;

[0263] Process C413: Please refer to Figure 2 , with Pi (1) and Ui (1) For the final state of the receiving substation, reverse power flow calculation is performed on the receiving substation to calculate the active power Po of the receiving substation (1) and transmission voltage Uo (1) ;

[0264] Process C4131: Calculate the load circuit current of the receiving substation, recorded as Iu; the calculation formula for Iu is: Iu = Pi (1) / Ui (1) ;

[0265] Calculate the power (complex power) consumed by the equivalent impedance of the receiving substation, denoted as Pz (1) =Iu×(Iu * )×dZ (1) ; Wherein, * represents conjugation;

[0266] Calculate the transmission voltage Uo at the receiving substation (1) , Uo (1) The calculation formula is: Uo (1)=(Iu×dZ (1) )+Ui (1) ;

[0267] Process C4132: Calculate the power (complex power) consumed by the equivalent admittance of the receiving substation, denoted as Py (1) =(Uo (1) ) 2 ×(dY (1) * ); where * represents conjugation;

[0268] Calculate the active power Po of the receiving substation (1) , Po (1) The calculation formula is: Po (1) =Pz (1) +Py (1) +Pi (1) ;

[0269] Process C42: Please refer to Figure 3 , obtain the length of the transmission line from the receiving substation to the first substation, denoted as lr ( o-1);

[0270] According to Po (1) and Uo (1) , calculate the transmission loss from the first substation to the receiving substation, and get the active power Po of the first substation (2) and transmission voltage Uo (2) ;

[0271] Process C421: Calculate the allocation coefficient from the first substation to the receiving substation, denoted as β (1) β (1) The calculation formula is: (1) =lr ( o-1) / lr;

[0272] Calculate the equivalent admittance dY of the transmission line from the first substation to the receiving substation (2,1) and dY (2,2) ;dY (2,1) and dY (2,2) The mathematical formula is: dY (2,1) =dY (2,2) =[(β (1) ×Gi (1) )+j×(β (1) ×Bi (1) )] / 2; where j represents the imaginary unit;

[0273] Calculate the equivalent impedance dZ of the transmission line from the first substation to the receiving substation (2) ;dZ (2) The mathematical formula is: dZ(2) =(β (1) ×Ri (1) )+j×(β (1) ×Xi (1) );

[0274] Process C422: Calculate the equivalent admittance dY (2,1) The power consumed is recorded as the equivalent admittance Py (2,1) ;Py (2,1) The calculation formula is: Py (2,1) =(Uo (1) ) 2 ×(dY (2,1) * ); where * represents conjugation;

[0275] Calculate the power of node A, denoted as S (1) ;S (1) The calculation formula is: S (1) =Py (2,1) +Po(1);

[0276] Process C423: Calculate the transmission voltage Uo (2) , Uo (2) The calculation formula is: Uo (2) =Uo (1) +[(S (1) / Uo (1) ) * ]×dZ (2) ; Wherein, * represents conjugation;

[0277] Calculate the equivalent impedance dZ (2) The power consumed is denoted as Pz (2) ;Pz (2) The calculation formula is: Pz (2) =(S (1) / Uo (2) ) 2 ×dY (2 ,2);

[0278] Calculate the power of node B, denoted as S (2) ;S (2) The calculation formula is: S (2) =Pz (2) +S (1) ;

[0279] Process C424: Calculate the active power Po of the first substation (2) , Po (2) The calculation formula is:

[0280] Po (2) =S (2) +Uo(2) ×[(Uo (2) ×dY (2 ,2)) * ];

[0281] Process C43: Get the actual receiving voltage of the first substation and record it as dU (1) , the actual power is recorded as dP (1) ;

[0282] Determine whether the output terminal of the first substation is connected to other substations, and (2) Adjust the working status of the first substation;

[0283] Process C431: The output terminal of the first substation is not connected to other substations; obtain the reactive power of the first substation and record it as Qo (2) , and adjust the working status of the first substation;

[0284] Process C4311: Calculate the actual received current of the first substation, denoted as dI (2) ;dI (2) The calculation formula is: dI (2) =dP (2) / dU (2) ;

[0285] The actual reactance of the first transformer is denoted as dX (2) ,

[0286] dX (2) The calculation formula is: dX (2) =[dU (2) ] / dI (2) ;

[0287] The actual susceptance of the first transformer, expressed in dB (2) ;

[0288] dB (2) The calculation formula is: dB (2) =[3×(dI (2) ) 2 ] 1 / 2 / dU (2) ;

[0289] The actual resistance of the first transformer is denoted as dR (2) ;

[0290] dR (2) The calculation formula is: dR (2) =[PG(o)×(dU (2) ) 2 ] / [mn×(dP (2) ) 2 ];

[0291] The actual conductance of the first transformer is denoted as dG (2) ;

[0292] dG (2) The calculation formula is: dG (2) =(dP (2) ) / [(dU (2) ) 2 ];

[0293] Process C4312: Calculate dX (2) The offset coefficient of the equivalent susceptance is denoted as zx (2) ;

[0294] zx (2) The calculation formula is: zx (2) =1-[dX (1) / (β (1) ×Xi (1) )];

[0295] Calculating dB (2) The offset coefficient with equivalent reactance is denoted as zb (2) ;

[0296] zb (2) The calculation formula is: (2) =1-[dB (2) / (β (1) ×Bi (1) )];

[0297] Calculating dR (2) The offset coefficient with the equivalent reactance is denoted as zr (2) ;

[0298] zr (2) The calculation formula is: (2) =1-[dR (2) / (β (1) ×Ri (1) )];

[0299] Calculate dG (2) The deviation coefficient from the equivalent reactance is denoted as zg (2) ;

[0300] zg (2) The calculation formula is: (2) =1-[β (1) / (dG (1) ×Gi (1) )];

[0301] Process C4313: Calculate zx (2) +zb (2) +zr(2) +zg (2) The sum of , recorded as zz (2) ;

[0302] Calculate the power factor of the first substation and record it as pQ. The calculation formula of pQ is: pQ=Po (2) / Qo (2) ;

[0303] Assume that the active power after adjustment of the first transformer is PP (1) , judge zz (1) Size determination PP (1) Calculation formula:

[0304] Ifzz (1) ≥0, then PP (1) The calculation formula is: =[(2-zz (1) ) 1 / 2 ]×Po (2) ;

[0305] Ifzz (1) <0, then PP (1) The calculation formula is: =[(1-zz (1) ) 1 / 2 ]×Po (2) ;

[0306] Process C4314: Let the reactive power Qo of the first substation be (2) The power factor after the change is recorded as pQa; the calculation formula of pQa is: pQa=PP (1) / Qo (2) ;

[0307] Let the power factor pQ of the first substation remain unchanged, and record the changed reactive power as Qa; the calculation formula of pQa is: Qa=PP (1) / pQ;

[0308] Calculate Qa relative to Qo (2) The rate of change of Qo (2) ); b(Qo (2) ) is calculated as: b(Qo (2) )=(Qa-Qo (2) ) / Qo (2) ;

[0309] Calculate the rate of change of pQa relative to pQ, recorded as b(pQ); the calculation formula for b(pQ) is: b(pQ) = (pQ - pQa) / pQ;

[0310] Process C4315: Comparison b(Qo (2) ) and the size of b(pQ);

[0311] If b(Qo (2) )≤b(pQ), then the active power of the first substation is adjusted to PP (1) , reactive power is adjusted to Qa;

[0312] If b(Qo (2) )>b(pQ), then the active power of the first substation is adjusted to PP (1) , reactive power remains unchanged;

[0313] Process C4316: Repeat the same process (i.e., process C41) to calculate the active power and transmission voltage of the receiving substation, according to Po (2) and Uo (2) Calculate the active power and transmission voltage of the second substation;

[0314] Process C432: The output terminal of the first substation is connected to other substations;

[0315] The other substations connected to the output end of the first substation are regarded as Class A substations. The power supply and supply voltage of each Class A substation are obtained. The length of the transmission line connecting each Class A substation to the first substation is obtained. The same process of calculating the output power and transmission voltage of the first substation is repeated (i.e., process C42) to calculate the new output power aPo of the first substation. (2) , calculate the new output power aUo of the first substation (2) ;

[0316] Po (2) Replaced by aPo (2) , Uo (2) Replaced with aUo (2) , re-execute process C431 and adjust the working status of the first substation;

[0317] Process C44: Repeat the same process (i.e., process C42 to process C43) of calculating the output power and transmission voltage of the first substation and adjusting the working state of the first substation, and adjust the second to the tr (1) Working status of each substation;

[0318] Process C5: Repeat the same process of adjusting the working status of each substation on line 1 to adjust line 2 to line mn.

[0319] Continuous scheduling module: used to continuously obtain historical electricity consumption and transmission data of each urban district and adjust the power distribution of each urban district.

[0320] The above formulas are all dimensionless and calculated by taking their numerical values. The formula is a formula for the most recent real situation obtained by collecting a large amount of data and performing software simulation. The preset parameters in the formula are set by technicians in this field according to actual conditions. For example, if there are weight coefficients and proportional coefficients, the size of the settings is to quantify each parameter to obtain a specific value, which is convenient for subsequent comparison. Regarding the size of the weight coefficient and the proportional coefficient, as long as it does not affect the proportional relationship between the parameter and the quantized value, it is fine.

[0321] Finally, it should be noted that the above-described embodiments are only specific implementations of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the above-described embodiments, those skilled in the art should understand that any person skilled in the art can modify or easily conceive of changes to the technical solutions described in the above-described embodiments within the technical scope disclosed by the present invention, or replace some of the technical features therein with equivalents. Such modifications, changes, or replacements do not deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A cloud computing-based intelligent control system for power grid operation, characterized in that: The system comprises: Data acquisition module: used to obtain the number of all municipal districts in the current month and the target area; according to the current month, obtain the historical electricity consumption corresponding to each municipal district, and estimate the expected electricity consumption of each municipal district; Data analysis module: used to obtain the no-load electromotive force, rated voltage, reactance voltage phase and current phase of the power supply end, and obtain the receiving voltage and receiving power of each municipal district as transmission data; calculate the expected power corresponding to the power supply end to each municipal district and the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance according to the transmission data, and obtain the equivalent transmission parameters; Power dispatching module: used to adjust the transformation ratio of the receiving substation corresponding to each municipal district according to the expected power; obtain the substation connection relationship of the power supply line corresponding to each municipal district, take the receiving substation as the starting point and the power supply end as the end point, calculate the reverse power flow for each municipal district according to the equivalent transmission parameters, adjust the working status of each substation, and allocate power to each municipal district; Continuous scheduling module: used to continuously obtain the historical power consumption and transmission data of each urban district and adjust the power distribution of each urban district.

2. The cloud computing-based intelligent control system for power grid operation according to claim 1, characterized in that: The workflow of the data analysis module is as follows: Process A1: Get the current month, recorded as m; get the number of districts in the target area, recorded as mn; Get the electricity consumption of the 1st to mnth municipal districts in month m in the past three years as the historical electricity consumption, and get: Electricity consumption of the 1st to mnth municipal districts in the third year: (1,3) ~he (mn,3) ; Electricity consumption of the first to mnth municipal districts in the second year: (1,2) ~he (mn,2) ; Electricity consumption in the first year of the 1st to mnth municipal districts: (1,1) ~he (mn,1) ; Process A2: Calculate the expected electricity consumption of the first city district, denoted as Ee (1) ; Process A21: Calculate the rate of change in electricity consumption from the third year to the second year in the first city district, recorded as up (3-2) ; up (3-2) The calculation formula is: up (3-2) =(he (1,2) -he (1,3) ) / he (1,3) ; The rate of change in electricity consumption from the second year to the first year is recorded as up (2-1) ; up (2-1) The calculation formula is: up (2-1) =(he (1,1) -he (1,2) ) / he (1,2) ; Process A22: Obtain the electricity consumption of the first municipal district in month (m-1), recorded as hhe; Get the current date, recorded as d; get the number of days in month m, recorded as dm; Get the cumulative electricity consumption of the first municipal district on day m of month d, recorded as ne; Process A23: According to up (3-2) 、up (2-1) , hhe and [(dm×ne) / d] construct the state transfer matrix, and iterate with the NumPy library function to obtain the steady-state transfer matrix, denoted as matrix CB (1) ; Define a 1×2 order receiving matrix {1 / 2, 1 / 2}, denoted as matrix CB (2) ; Calculate the matrix CB (1) Multiply by matrix CB (2) , we get the matrix CB (3) ; Calculate the matrix CB (3) The average value of each parameter in is taken as the expected electricity consumption of the first municipal district in the mth month, denoted as Ee (1) ; Process A3: Repeat the calculation of the expected electricity consumption Ee of the first municipal district (1) The same process is used to calculate the expected electricity consumption of the 2nd to mnth municipal districts, and Ee is obtained. (1) ~Ee (mn) , and enter the data analysis module.

3. The cloud computing-based intelligent control system for power grid operation according to claim 2 is characterized in that: The workflow of the data analysis module is as follows: Process B1: Obtain the no-load electromotive force of the power supply end of the municipal power grid in the target area, denoted as Eq (o) , rated voltage is recorded as Un (o) , reactance is denoted as Xd (o) ; The voltage phase at the power supply end is obtained and recorded as φ (U) , the current phase is denoted as φ (I) ; Calculate the power factor angle at the power supply end, recorded as θ; the calculation formula for θ is: θ = φ (U) -φ (I) ; Calculate the active power at the power supply end, denoted as PG (o) PG (o) The calculation formula is: PG (o) =[Eq (o) ×Un (o) ] / Xd (o) ×cos(θ); Process B2: Get the expected electricity consumption corresponding to the 1st, 2nd, and up to the mnth municipal districts: Ee (1) ,Ee (2) ~Ee (mn) ; Calculate Ee (1) ,Ee (2) ~Ee (mn) The ratio of , we get: bl (1) ,bl (2) ~bl (mn) ; Get the receiving voltage corresponding to the 1st, 2nd, and up to the mnth municipal districts, denoted as: Ui (1) 、Ui (2) ~Ui (mn) ; Get the receiving power corresponding to the 1st, 2nd, and up to the mnth municipal districts, denoted as: Pj (1) 、Pj (2) ~Pj (mn) ; Process B3: Calculate bl (1) ~bl (mn) The sum of is denoted as abl; calculate Pj (1) ~Pj (mn) The sum of is denoted as aPj; Calculate the expected power of the first city district, denoted as Pi (1) ;Pi (1) The calculation formula is: Pi (1) =(bl (1) / abl)×aPj; Calculate the expected power of the second city district, denoted as Pi (2) ;Pi (2) The calculation formula is: Pi (2) =(bl (2) / abl)×aPj; Similarly, calculate the expected power of the mnth city district, denoted as Pi (mn) ;Pi (mn) The calculation formula is: Pi (mn) =(bl (mn) / abl)×aPj; Process B4: According to Ui (1) ~Ui (mn) and Pi (1) ~Pi (mn) , calculate the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance corresponding to the power supply end to each municipal district in turn; Process B5: Summarize the equivalent reactances corresponding to the 1st to mnth municipal districts: Xi (mn) ~Xi (mn) , equivalent susceptance: Bi (mn) ~Bi (mn) , equivalent resistance: Ri (mn) ~Ri (mn) , equivalent conductance: Gi (mn) ~Gi (mn) , as equivalent transmission parameters, and enter the power dispatching module.

4. The cloud computing-based intelligent control system for power grid operation according to claim 3 is characterized in that: The specific process of process B4 is as follows: Process B41: supplying power to the first municipal district; Process B411: Calculate the receiving current of the first municipal district, denoted as Ii (1) ; Ii (1) The calculation formula is: Ii (1) =Pi (1) / Ui (1) ; Calculate the equivalent reactance from the power supply end to the first municipal district, denoted as Xi (1) , Xi (1) The calculation formula is: Xi (1) =(Un (o) -Ui (1) ) / Ii (1) ; Process B412: Calculate the equivalent susceptance from the power supply end to the first municipal district, denoted as Bi (1) ; Bi (1) The calculation formula is: Bi (1) ={[3×(Ii (1) ) 2 ] 1 / 2 } / Ui (1) ; Calculate the equivalent resistance from the power supply end to the first municipal district, denoted as Ri (1) ; Ri (1) The calculation formula is: Ri (1) =[PG(o)×(Ui (1) ) 2 ] / [(Pi (1) ) 2 ]; Calculate the equivalent conductance from the power supply end to the first municipal district, denoted as Gi (1) ; Gi (1) The calculation formula is: Gi (1) =(Pi (1) ) / [(Ui (1) ) 2 ]; Process B42: supplying power to the second municipal district; Process B421: Calculate the receiving current of the second municipal district, denoted as Ii (2) ; Ii (2) The calculation formula is: Ii (2) =Pi (2) / Ui (2) ; Calculate the equivalent reactance from the power supply end to the second municipal district, denoted as Xi (2) , Xi (2) The calculation formula is: Xi (2) =(Un (o) -Ui (2) ) / Ii (2) ; Process B422: Calculate the equivalent susceptance from the power supply end to the second municipal district, denoted as Bi (2) ; Bi (2) The calculation formula is: Bi (2) ={[3×(Ii (2) ) 2 ] 1 / 2 } / Ui (2) ; Calculate the equivalent resistance from the power supply end to the second municipal district, denoted as Ri (2) ; Ri (2) The calculation formula is: Ri (2) =[PG(o)×(Ui (2) ) 2 ] / [(Pi (2) ) 2 ]; Calculate the equivalent conductance from the power supply end to the second municipal district, denoted as Gi (2) ; Gi (2) The calculation formula is: Gi (2) =(Pi (2) ) / [(Ui (2) ) 2 ]; Process B43: Similarly, calculate the equivalent reactance Xi from the power supply end to the mnth municipal district (mn) , the equivalent electrical susceptance is denoted as Bi (mn) , equivalent resistance Ri (mn) , equivalent conductance Gi (mn) .

5. The cloud computing-based intelligent control system for power grid operation according to claim 3, characterized in that: The workflow of the power dispatch module is as follows: Process C1: Construct the transmission relationship map of each municipal area based on equivalent transmission parameters; The first municipal district: {PG (o) , (o) }—{Xi (1) 、Bi (1) 、Ri (1) , Gi (1) }—{Pi (1) , Ui (1) }; Second Municipal District: {PG (o) , (o) }—{Xi (2) 、Bi (2) 、Ri (2) , Gi (2) }—{Pi (2) , Ui (2) }; And so on, the mnth district: {PG (o) , (o) }—{Xi (mn) 、Bi (mn) 、Ri (mn) , Gi (mn) }—{Pi (mn) , Ui (mn) }; Process C2: The power lines that transmit power from the power supply end to the 1st, 2nd, and even the mnth municipal districts are recorded as line 1, line 2 to line mn; Get the length of line 1, denoted as ll (1) ; Get the length of line 2, recorded as ll (2) ; And so on, get the length of line mn, recorded as ll (mn) ; Count the number of substations on line 1, denoted as tr (1) ; The number of substations on line 2, denoted as tr (2) ; Similarly, the number of substations on line mn is denoted as tr (mn) ; Get the active power PG of the power supply end (o) , get the expected power of the 1st, 2nd, and up to the mnth districts: Pi (1) 、Pi (2) ~Pi (mn) ; Process C3: Adjust the transformation ratio of receiving substations in the 1st to mnth municipal districts; Calculate the transformation ratio of the receiving substation in the first municipal district, denoted as b (o-1) ; b (o-1) The calculation formula is: (o-1) =Pi (1) / PG (o) ; Calculate the transformation ratio of the receiving substation in the second municipal district, denoted as b (o-2) ; b (o-2) The calculation formula is: (o-2) =Pi (2) / PG (o) ; By analogy, the transformation ratio of the receiving substation in the mnth municipal district is denoted as b (o-mn) ; b (o-mn) The calculation formula is: (o-mn) =Pi (mn) / PG (o) ; Among them, b (o-1) ~b (o-mn) are all rational numbers; According to b (o-1) ~b (o-mn) Adjust the transformation ratio of the receiving substations corresponding to the 1st to mnth municipal districts; Process C4: According to the data in processes C1 to C3, the working status of each substation on line 1 is adjusted; Process C5: Repeat the same process of adjusting the working status of each substation on line 1 to adjust line 2 to line mn.

6. The cloud computing-based intelligent control system for power grid operation according to claim 5, characterized in that: The specific process of process C4 is as follows: Process C41: Obtain the actual receiving voltage of the receiving substation and record it as dU (1) , actual power is recorded as dP (1) According to Pi (1) and Ui (1) , calculate the active power Po of the receiving substation (1) And the transmission voltage Uo (1) ; Process C42: Get the length of the transmission line from the receiving substation to the first substation, denoted as lr (o-1) ; According to Po (1) and Uo (1) , calculate the transmission loss from the first substation to the receiving substation, and get the active power Po of the first substation (2) And the transmission voltage Uo (2) ; Process C43: Get the actual receiving voltage of the first substation and record it as dU (1) , actual power is recorded as dP (1) ; Determine whether the output terminal of the first substation is connected to other substations, and (2) Adjust the working status of the first substation; Process C44: Repeat the same process of calculating the output power and transmission voltage of the first substation and adjusting the working state of the first substation to adjust the second to third substations on line 1. (1) The working status of a substation.

7. The cloud computing-based intelligent control system for power grid operation according to claim 6, characterized in that: The specific process of process C41 is as follows: Process C411: Calculate the actual receiving current of the receiving substation, denoted as dI (1) ;dI (1) The calculation formula is: dI (1) =dP (1) / dU (1) ; Calculate the actual reactance of the receiving substation, denoted as dX (1) ; dX (1) The calculation formula is: dX (1) =[dU (1) -Ui (1) ] / dI (1) ; Among them, Ui (1) Indicates the receiving voltage of the first municipal district; Calculate the actual susceptance of the receiving substation, expressed in dB (1) ; dB (1) The calculation formula is: dB (1) ={[3×(dI (1) ) 2 ] 1 / 2 } / dU (1) ; Calculate the actual resistance of the receiving substation, denoted as dR (1) ; dR (1) The calculation formula is: dR (1) =[PG (o) ×(dU (1) ) 2 ] / [mn×(dP (1) ) 2 ]; Calculate the actual conductance of the receiving substation, denoted as dG (1) ; dG (1) The calculation formula is: dG (1) =(dP (1) ) / [(dU (1) ) 2 ] Process C412: The equivalent admittance of the receiving substation is recorded as dY (1) ;dY (1) The mathematical expression is: dY (1) =dG (1) +(dB (1) ×j); where j represents an imaginary unit; The equivalent impedance of the receiving substation is denoted as dZ (1) ;dZ (1) The mathematical expression of dZ is: (1) =dR (1) +(dX (1) ×j); where j represents an imaginary unit; Process C413: Pi (1) and Ui (1) For the final state of the receiving substation, reverse power flow calculation is performed on the receiving substation to calculate the active power Po of the receiving substation (1) And the transmission voltage Uo (1) ; Process C4131: Calculate the load circuit current of the receiving substation, recorded as Iu; the calculation formula for Iu is: Iu = Pi (1) / Ui (1) ; Calculate the power consumed by the equivalent impedance of the receiving substation, denoted as Pz (1) =Iu×(Iu * )×dZ (1) ; Wherein, * represents conjugation; Calculate the transmission voltage Uo at the receiving substation (1) , Uo (1) The calculation formula is: Uo (1) =(Iu×dZ (1) )+Ui (1) ; Process C4132: Calculate the power consumed by the equivalent admittance of the receiving substation, denoted as Py (1) =(Uo (1) ) 2 ×(dY (1) * ); where * represents conjugation; Calculate the active power Po of the receiving substation (1) , Po (1) The calculation formula is: Po (1) =Pz (1) +Py (1) +Pi (1) .

8. The cloud computing-based intelligent control system for power grid operation according to claim 6, characterized in that: The specific process of process C42 is as follows: Process C421: Calculate the allocation coefficient from the first substation to the receiving substation, denoted as β (1) β (1) The calculation formula is: (1) =lr (o-1) / lr; Calculate the equivalent admittance dY of the transmission line from the first substation to the receiving substation (2,1) and dY (2,2) ;dY (2,1) and dY (2,2) The mathematical formula is: dY (2,1) =dY (2,2) =[(β (1) ×Gi (1) )+j×(β (1) ×Bi (1) )] / 2; where j represents an imaginary unit; Calculate the equivalent impedance dZ of the transmission line from the first substation to the receiving substation (2) ;dZ (2) The mathematical formula is: dZ (2) =(β (1) ×Ri (1) )+j×(β (1) ×Xi (1) ); Process C422: Calculate the equivalent admittance dY (2,1) The power consumed is recorded as the equivalent admittance Py (2,1) ;Py (2,1) The calculation formula is: (2,1) =(Uo (1) ) 2 ×(dY (2,1) * ); where * represents conjugation; Calculate the power of node A, denoted as S (1) ; S (1) The calculation formula is: S (1) =Py (2,1) +Po (1) ; Process C423: Calculate the transmission voltage Uo (2) , Uo (2) The calculation formula is: Uo (2) =Uo (1) +[(S (1) / Uo (1) ) * ]×dZ (2) ; Wherein, * represents conjugation; Calculate the equivalent impedance dZ (2) The power consumed is denoted as Pz (2) ;Pz (2) The calculation formula is: Pz (2) =(S (1) / Uo (2) ) 2 ×D (2,2) ; Calculate the power of node B, denoted as S (2) ; S (2) The calculation formula is: S (2) =Pz (2) +S (1) ; Process C424: Calculate the active power Po of the first substation (2) , Po (2) The calculation formula is: Po (2) =S (2) +Uo (2) ×[(Uo (2) ×dY (2,2) ) * ]。 9. The cloud computing-based intelligent control system for power grid operation according to claim 6, characterized in that: The specific process of process C43 is as follows: Process C431: The output end of the first substation is not connected to other substations; obtain the reactive power of the first substation and record it as Qo (2) , and adjust the working status of the first substation; Process C432: The output terminal of the first substation is connected to other substations; Take the other substations connected to the output end of the first substation as Class A substations, obtain the power supply power and power supply voltage of each Class A substation, obtain the length of the transmission line connecting each Class A substation to the first substation, repeat the same process of calculating the output power and transmission voltage of the first substation, and calculate the new output power aPo of the first substation (2) , calculate the new output power aUo of the first substation (2) ; Po (2) Replaced by aPo (2) , will Uo (2) Replaced by aUo (2) , re-execute process C431 and adjust the working status of the first substation.

10. The cloud computing-based intelligent control system for power grid operation according to claim 9, characterized in that: The specific process of process C431 is as follows: Process C4311: Calculate the actual received current of the first substation, denoted as dI (2) ;dI (2) The calculation formula is: dI (2) =dP (2) / dU (2) ; The actual reactance of the first transformer is denoted by dX (2) , dX (2) The calculation formula is: dX (2) =[dU (2) ] / dI (2) ; The actual susceptance of the first transformer, expressed in dB (2) ; dB (2) The calculation formula is: dB (2) =[3×(dI (2) ) 2 ] 1 / 2 / dU (2) ; The actual resistance of the first transformer is denoted as dR (2) ; dR (2) The calculation formula is: dR (2) =[PG(o)×(dU (2) ) 2 ] / [mn×(dP (2) ) 2 ]; The actual conductance of the first transformer is denoted as dG (2) ; dG (2) The calculation formula is: dG (2) =(dP (2) ) / [(dU (2) ) 2 ]; Process C4312: Calculate dX (2) The offset coefficient with equivalent susceptance is denoted as zx (2) ; zx (2) The calculation formula is: zx (2) =1-[dX (1) / (β (1) ×Xi (1) )]; Calculating dB (2) The deviation coefficient from the equivalent reactance is denoted as zb (2) ; zb (2) The calculation formula is: (2) =1-[dB (2) / (β (1) ×Bi (1) )]; Calculate dR (2) The deviation coefficient from the equivalent reactance is denoted as zr (2) ; z (2) The calculation formula is: (2) =1-[dR (2) / (β (1) ×Ri (1) )]; Calculate dG (2) The deviation coefficient from the equivalent reactance is denoted as zg (2) ; z (2) The calculation formula is: (2) =1-[β (1) / (dG (1) ×Gi (1) )]; Process C4313: Calculate zx (2) +zb (2) +zr (2) +zg (2) The sum of is denoted as zz (2) ; Calculate the power factor of the first substation and record it as pQ. The calculation formula of pQ is: pQ = Po (2) / Qo (2) ; Assume that the active power after adjustment of the first transformer is PP (1) , judge zz (1) Size determination PP (1) Calculation formula of If zz (1) ≥0, then PP (1) The calculation formula is: =[(2-zz (1) ) 1 / 2 ]×Po (2) ; If zz (1) <0, then PP (1) The calculation formula is: =[(1-zz (1) ) 1 / 2 ]×Po (2) ; Process C4314: Let the reactive power Qo of the first substation be (2) The power factor after the change is recorded as pQa; the calculation formula of pQa is: pQa = PP (1) / Qo (2) ; Let the power factor pQ of the first substation remain unchanged, and record the reactive power after the change as Qa; the calculation formula of pQa is: Qa=PP (1) / pQ; Calculate Qa relative to Qo (2) The rate of change of Qo (2) ); b(Qo (2) ) is calculated as follows: b(Qo (2) )=(Qa-Qo (2) ) / Qo (2) ; Calculate the rate of change of pQa relative to pQ, recorded as b(pQ); The calculation formula of b(pQ) is: b(pQ) = (pQ - pQa) / pQ; Process C4315: Comparison b(Qo (2) ) and the size of b(pQ); If b(Qo (2) )≤b(pQ), then the active power of the first substation is adjusted to PP (1) , reactive power is adjusted to Qa; If b(Qo (2) )>b(pQ), then the active power of the first substation is adjusted to PP (1) , reactive power remains unchanged; Process C4316: Repeat the same process of calculating the active power and transmission voltage of the receiving substation, according to Po (2) and Uo (2) Calculate the active power and transmission voltage of the second substation.

Images