A power big data analysis method and system based on cloud computing
Through cloud -based electricity big data analysis methods, the problems of lack of systemic lack of systematic, data quality and human factors in the existing power analysis system have been solved, and the efficiency of power analysis is improved and the power supply is optimized.
Patent Information
- Application Number
- CN202411533468.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-31
- Publication Date
- 2025-01-24
- Estimated Expiration
- 2044-10-31
AI Technical Summary
The existing power analysis methods and systems have systematic shortcomings, data quality problems and lack of human factors, resulting in low power analysis efficiency.
The power big data analysis method based on cloud computing is used to obtain the power consumption data and transmission data of the municipal area, calculate the priority coefficient and equivalent power parameter values, adjust the transformer ratio, and realize the automation and dynamic adjustment of power distribution.
It improves the data processing capabilities and analysis accuracy of the power system, optimizes power supply, improves system stability and operating efficiency, reduces power supply costs and enhances power supply reliability.
Smart Images

Figure CN119047797B_ABST
Abstract
Description
Technical Field
[0001] The present invention discloses a power big data analysis method and system based on cloud computing, and relates to the field of power analysis. Background Art
[0002] The existing power analysis methods and analysis systems have the following deficiencies:
[0003] Lack of systematicity: In the existing power system, there is a lack of systematic planning in the application stage of cloud computing technology. Specifically, the various information systems are not highly correlated with each other, and multiple databases cannot share resource information with each other, which hinders the orderly progress of the information construction of the power system.
[0004] Data quality issues: The existing system acquires too much raw data without any correlation, resulting in large errors in the process of data analysis and data fitting. At the same time, too much raw data will also lead to problems such as data distortion, incompleteness, and inaccuracy during the acquisition process, reducing the credibility of the system analysis results.
[0005] Human factors: There is currently a lack of professional talent in the power system, which will increase the difficulty and cost of implementing cloud computing technology applications. Summary of the invention
[0006] In view of the shortcomings of the prior art, the purpose of the present invention is to provide a cloud computing-based power big data analysis method and system, aiming to solve the problem of low power analysis efficiency.
[0007] In order to achieve the above object, the present invention is implemented by the following technical solution: A power big data analysis method based on cloud computing includes:
[0008] Step S1: Obtain the number of all municipal districts in the target area; obtain the number of residents, current power consumption and historical power consumption of each municipal district to obtain power consumption data; calculate the priority coefficient of each municipal district based on the power consumption data to obtain an allocation sequence;
[0009] Step S2: Obtain the no-load electromotive force, rated voltage and reactance of the power generation end, obtain the receiving voltage and receiving power of each city district, and obtain the power transmission data; construct a mathematical relationship between the power generation end, the power transmission end and the power receiving end according to the power transmission data, and calculate the equivalent power parameter value of the power supply line of each city district according to the mathematical relationship;
[0010] Step S3: Obtain the series-parallel relationship of the transformers between each power supply line to obtain the connection data of the transformers; based on the allocation sequence, adjust the transformer ratio of each power supply line according to the connection data and equivalent power parameter value of the transformer, and allocate power to each city district;
[0011] Step S4: Continuously obtain the current power consumption of each city district, and update and adjust the power distribution of each city district.
[0012] Furthermore, the specific steps of step S1 are as follows:
[0013] Step S11: Obtain the number of all city districts in the target area, denoted as mn;
[0014] Get the number of residents corresponding to each city district and record it as pn1, pn2~pn mn , the current power consumption is recorded as oj1, oj2~oj mn ;
[0015] Among them, pn1 and oj1 represent the number of residents and current power consumption in the first city district respectively; pn2 and oj2 represent the number of residents and current power consumption in the second city district respectively; and so on. mn and oj mn They represent the number of residents and current power consumption in the mnth municipal district respectively;
[0016] Step S12: Extract oj1~oj mn The maximum value is recorded as ojm, the minimum value is recorded as ojl, and the mode is recorded as ojz; calculate oj1~oj mn The average value of is denoted as aoj;
[0017] The relay power is recorded as ojc; compare ojz with aoj to determine the value of ojc;
[0018] If ojz>aoj, then the value of ojc is ojz;
[0019] If ojz≤aoj, then the value of ojc is aoj;
[0020] Step S13: Calculate the difference between ojm and ojl, denoted as ooj;
[0021] Assume that the current power consumption of the i-th district is oj i , the value range of i is: 1~mn;
[0022] The priority allocation coefficient corresponding to the i-th municipal district is denoted as α(i); the calculation formula of α(i) is defined as:
[0023] Judge oj i Whether it is within the interval [ojc, ojm];
[0024] If , then the calculation formula of α(i) is as follows:
[0025] α(i)=[(oj i -ojc) / (ojm-ojc)]×0.5;
[0026] If not, the calculation formula of α(i) is as follows:
[0027] α(i)=[(oj i -ojl) / (ojc-ojl)]×0.5;
[0028] Step S14: oj1~oj mn Substituting into the calculation formula of α(i), we can obtain the priority allocation coefficients corresponding to the 1st to the mnth municipal districts: α(1), α(2)~α(mn);
[0029] Among them, α(1) represents the priority allocation coefficient of the first municipal district; α(2) represents the priority allocation coefficient of the second municipal district; and so on, α(mn) represents the priority allocation coefficient of the mnth municipal district;
[0030] Obtain the electricity consumption of each municipal district in the past three years as historical electricity consumption;
[0031] The electricity consumption of the first city district in three years is: ahj(1), bhj(1) and chj(1); among them, ahj(1) represents the electricity consumption of the first city district in the third year, bhj(1) represents the electricity consumption of the first city district in the second year, and chj(1) represents the electricity consumption of the first city district in the first year;
[0032] The electricity consumption of the second city district in three years is: ahj(2), bhj(2) and chj(2); among them, ahj(2) represents the electricity consumption of the second city district in the third year, bhj(2) represents the electricity consumption of the second city district in the second year, and chj(2) represents the electricity consumption of the second city district in the first year;
[0033] By analogy, the power consumption of the mn-th municipal district in three years is: ahj(mn), bhj(mn) and chj(mn); among them, ahj(mn) represents the power consumption of the mn-th municipal district in the third year, bhj(mn) represents the power consumption of the mn-th municipal district in the second year, and chj(mn) represents the power consumption of the mn-th municipal district in the first year;
[0034] Step S15: Calculate the electricity consumption change rate of each city district according to the historical power consumption to obtain ub(1), ub(2)~ub(mn);
[0035] 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 city district; and so on, ub(mn) represents the electricity consumption change rate of the mnth city district;
[0036] Step S16: Calculate the priority coefficient of the first municipal district, denoted as ff(1); the calculation formula of ff(1) is: ff(1)=pn1×α(1)×ub(1);
[0037] Calculate the priority coefficient of the second municipal district, denoted as ff(2); the calculation formula of ff(2) is: ff(2) = pn2×α(2)×ub(2);
[0038] Similarly, the priority coefficient of the mnth district is calculated, denoted as ff(mn); the calculation formula of ff(mn) is: ff(mn) = pn mn ×α(mn)×ub(mn);
[0039] The priority coefficients ff(1) to ff(mn) corresponding to the 1st to mnth municipal districts are summarized to obtain the allocation sequence of the municipal districts, and then the process goes to step S2.
[0040] Furthermore, the specific steps of step S15 are as follows:
[0041] Step S151: Calculate the power consumption change rate of the first city district according to ahj(1), bhj(1) and chj(1), denoted as ub(1);
[0042] 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) =[bhj(1) / chj(1)]-1;
[0043] 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) ; up (2-1) The calculation formula is: up (2-1) =[ahj(1) / bhj(1)]-1;
[0044] Step S152: According to up (3-2) 、up (2-1) , ahj(1) and bhj(1), construct the state transfer matrix, denoted as matrix B; the mathematical expression of matrix B is as follows:
[0045] ; Where n represents the number of iterations of matrix B, × represents matrix multiplication; n is a natural number, and its initial value is 1; each time matrix B is iterated, n is increased by 1 on the original basis;
[0046] Use the NumPy library function to iterate the matrix B and obtain the steady-state transfer matrix, which is recorded as matrix B1;
[0047] Define a 1×2 order receiving matrix {1, 0}, 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 element in the matrix B3 as the estimated power consumption of the first municipal district, denoted as jj(1);
[0048] Step S153: Calculate the power consumption change rate ub(1) of the first city district, where ub(1) is calculated as follows: ub(1)=[jj(1) / ahj(1)]-1;
[0049] Step S154: Repeat the same steps of calculating the power consumption change rate of the first urban district, calculate the power consumption change rates of the second to mn-th urban districts, and obtain ub(2)~ub(mn).
[0050] Furthermore, the specific steps of step S2 are as follows:
[0051] Step S21: Obtain the no-load electromotive force of the generating end, recorded as Eq(o), the rated voltage as Un(o), and the reactance as Xd(o); obtain the current frequency of the generating end, recorded as f;
[0052] Calculate the power generation period, recorded as T; the calculation formula of T is: T = 1 / f;
[0053] Step S22: Get the current time, recorded as dt; calculate the remainder of dt and T, recorded as cd; the calculation formula of cd is: cd = dt%T; where % represents the modulus symbol;
[0054] Taking the current as the positive direction, calculate the phase angle between voltage and current, which is recorded as θ. The calculation formula of θ is: θ = (cd / T) × 2π; where π represents pi;
[0055] Calculate the active power of the power generation equipment, recorded as PG (o); the calculation formula of PG (o) is:
[0056] PG(o)=[Eq(o)×Un(o)] / Xd(o)×sin(θ);
[0057] Step S23: Obtain the transformer voltage connected to the power grid corresponding to each municipal district (within the target area) as the received voltage; and the transformer power as the received power;
[0058] The receiving voltages corresponding to the 1st to mnth municipal districts are recorded as: Ui1, Ui2~Ui mn ;
[0059] The received power is recorded as: Pi1, Pi2~Pi mn ;
[0060] Among them, Ui1 and Pi1 represent the receiving voltage and receiving power of the first city district respectively; Ui2 and Pi2 represent the receiving voltage and receiving power of the second city district respectively; and so on, Ui mn and Pi mn represent the receiving voltage and receiving power of the mn-th municipal district respectively;
[0061] Step S24: according to the data in steps S21 to S23, the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance from the power generation end to each municipal district are calculated in sequence;
[0062] Step S25: Summarize the equivalent reactances corresponding to the 1st to mnth municipal districts: XT1, XT2~XT mn ; Equivalent susceptance: BT1, BT2~BT mn , equivalent resistance: RT1, RT2~RT mn , equivalent conductance: GT1, GT2~GT mn , enter step S3 as the equivalent power parameter value.
[0063] Furthermore, the specific steps of step S24 are as follows:
[0064] Step S241: the power generation end reaches the first municipal district;
[0065] Step S2411: Calculate the receiving current of the first city district, denoted as Ii1;
[0066] The calculation formula of Ii1 is: Ii1=Pi1 / Ui1;
[0067] Calculate the equivalent reactance from the generator to the first municipal district, denoted as XT1,
[0068] The calculation formula of XT1 is: XT1 = (Un (o) - Ui1) / Ii1;
[0069] Step S2412: Calculate the equivalent susceptance from the power generation end to the first municipal district, denoted as BT1;
[0070] The calculation formula of BT1 is: BT1 = [3 × (Ii1) 2 ] 1 / 2 / Ui1;
[0071] Calculate the equivalent resistance from the generator end to the first municipal district, denoted as RT1;
[0072] The calculation formula of RT1 is: RT1 = [PG (o) × (Ui1) 2 ] / [(Pi1) 2 ];
[0073] Calculate the equivalent conductance from the power generation end to the first municipal district, denoted as GT1;
[0074] The calculation formula of GT1 is: GT1 = (Pi1) / [(Ui1) 2 ];
[0075] Step S242: the power generation end reaches the second municipal district;
[0076] Step S2421: Calculate the receiving current of the second city district, recorded as Ii2;
[0077] The calculation formula of Ii2 is: Ii2=Pi2 / Ui2;
[0078] Calculate the equivalent reactance from the generator to the second municipal district, denoted as XT2,
[0079] The calculation formula of XT2 is: XT2 = (Un (o) - Ui2) / Ii2;
[0080] Step S2422: Calculate the equivalent susceptance from the power generation end to the second municipal district, denoted as BT2;
[0081] The calculation formula of BT2 is: BT2 = [3 × (Ii2) 2 ] 2 / 2 / Ui2;
[0082] Calculate the equivalent resistance from the generator end to the second municipal district, recorded as RT2;
[0083] The calculation formula of RT2 is: RT2 = [PG (o) × (Ui2) 2 ] / [(Pi2) 2 ];
[0084] Calculate the equivalent conductance from the generating end to the second municipal district, denoted as GT2;
[0085] The calculation formula of GT2 is: GT2 = (Pi2) / [(Ui2) 2 ];
[0086] Step S243: Similarly, the power generation end reaches the mnth municipal district;
[0087] Step S2431: Calculate the receiving current of the mnth city district, denoted as Ii mn ;
[0088] Ii mn The calculation formula is: Ii mn =Pi mn / Ui mn ;
[0089] Calculate the equivalent reactance from the generator to the mnth municipal district, denoted as XT mn ,
[0090] XT mn The calculation formula is: XT mn = (Un (o) - Ui mn ) / Ii mn ;
[0091] Step S2432: Calculate the equivalent susceptance from the power generation end to the mnth municipal district, denoted as BT mn ;
[0092] BT mn The calculation formula is: BT mn =[3×(Ii mn ) 2 ] mn / 2 / Ui mn ;
[0093] Calculate the equivalent resistance from the generator end to the mnth municipal district, denoted as RT mn ;
[0094] RT mn The calculation formula is: RT mn =[PG(o)×(Ui mn ) 2 ] / [(Pi mn ) 2 ];
[0095] Calculate the equivalent conductance from the power generation end to the mnth municipal district, denoted as GT mn ;
[0096] GT mn The calculation formula is: GT mn =(Pi mn ) / [(Ui mn ) 2 ].
[0097] Furthermore, the specific steps of step S3 are as follows:
[0098] Step S31: constructing a power transmission relationship map for each city district according to the equivalent power parameter value;
[0099] The first municipal district: {PG (o), Un (o)}—{XT1, BT1, RT1, GT1}—{Pi1, Ui1};
[0100] The second municipal district: {PG (o), Un (o)}—{XT2, BT2, RT2, GT2}—{Pi2, Ui2};
[0101] And so on, the mnth district: {PG(o), Un(o)}—{XT mn , BT mn , RTmn , GT mn}—{Pi mn , Ui mn};
[0102] Step S32: Obtain the priority coefficients corresponding to the 1st to mnth municipal districts: ff(1), ff(2) to ff(mn);
[0103] Calculate the value of ff(1):ff(2):~:ff(mn), recorded as ap1:ap2:~:ap mn ;(i.e. ff(1):ff(2):~:ff(mn)=ap1:ap2:~:ap mn ) Among them, ap1~ap mn For positive numbers, round down;
[0104] Extraction of ap1~ap mn The minimum value of apl is recorded as apl, the received power corresponding to apl is recorded as Pil, and the initial transformer ratio from the power generation end to each municipal district is calculated;
[0105] Step S321: Calculate the standard transformation ratio of the initial transformer, denoted as fb;
[0106] The calculation formula of fb is: fb = [PG (o) / Pil] 1 / 2 ;
[0107] Step S322: Calculate the initial transformer ratio for power transmission from the power generation end to the first municipal district, denoted as B (o-1) ; B (o-1) The calculation formula is: (o-1) = (apl / ap1) × fb;
[0108] Calculate the initial transformer ratio for power transmission from the generator to the second municipal district, denoted as B (o-2) ; B (o-2) The calculation formula is: (o-2) = (apl / ap2) × fb;
[0109] By analogy, calculate the initial transformer ratio for power transmission from the power generation end to the mnth municipal district, denoted as B (o-mn) ; B (o-mn) The calculation formula is: (o-mn) = (apl / ap mn )×fb;
[0110] Step S33: adjusting the transformer ratio of the first power supply line according to the power transmission relationship mapping of the first municipal district;
[0111] Step S34: Repeat the same steps of adjusting the transformation ratio of the transformer of the first power supply line to adjust the transformation ratios of the transformers of the second to mnth power supply lines.
[0112] Furthermore, the specific steps of step S33 are as follows:
[0113] Step S331: For the first power supply line, obtain the number of transformers other than the initial transformer and the receiving transformer, recorded as tr; obtain the length of the power transmission line from the initial transformer to the receiving transformer, recorded as lr;
[0114] Step S332: the first power supply line, the initial transformer to the first transformer;
[0115] Step S333: The first power supply line, the first transformer to the second transformer;
[0116] Determine whether the second transformer is connected in parallel with other transformers;
[0117] If not connected in parallel, proceed to step S334;
[0118] If in parallel, then obtain the transformer information connected in parallel with the second transformer and proceed to step S335;
[0119] Step S334: the second transformer is not connected in parallel with other transformers, and the adjusted transformation ratio of the second transformer is calculated;
[0120] Step S335: the second transformer is connected in parallel with other transformers, and the adjusted transformation ratio of the second transformer is calculated;
[0121] Step S3351: other transformers connected in parallel with the two transformers are used as parallel transformers;
[0122] Count the number of transformers in parallel, denoted as oh;
[0123] Get the current transformation ratio of the 1st to ohth parameter transformer and record it as: OB1, OB2~OB oh ;
[0124] Among them, OB1 represents the current transformation ratio of the first parallel transformer; OB2 represents the current transformation ratio of the second parallel transformer; and so on. oh Indicates the current transformation ratio of the ohth parallel transformer;
[0125] Step S3352: Using the transformation ratio DQ of the second transformer (2) As the first term, calculate DQ (2) :OB1:OB2:~:OB oh The ratio is recorded as: uu1:uu2:~:uu oh :uu (oh+1); (ie DQ (2) :OB1:OB2:~:OB oh =uu1:uu2:~:uu oh :uu (oh+1) )
[0126] Calculate uu1:uu2:~:uu oh :uu (oh+1) The inverse ratio is expressed as: nn1:nn2:~:nn oh :nn (oh+1) ;
[0127] Step S3353: Repeat the same steps of calculating the adjusted transformation ratio of the second transformer on the premise that the second transformer is not connected in parallel with other transformers, and determine QQ (2) The value of QQ (2) Indicates the transformation ratio of the second transformer after adjustment;
[0128] Calculate nn1 + nn2 + ~ + nn oh +nn (oh+1) ; and is recorded as ann;
[0129] The transformation ratio of the second transformer after secondary adjustment is denoted as cQQ (2) ;cQQ (2) The calculation formula is:
[0130] QQQ (2) =(QQ (2) )×(nn1 / ann);
[0131] Step S336: Repeat the same steps of step S332 to step S335 to adjust the transformation ratios of the 3rd to trth transformers in sequence.
[0132] Furthermore, the specific steps of step S332 are as follows:
[0133] Step S3321: Obtain the length of the transmission line from the initial transformer to the first transformer, denoted as lr (o-1) ;
[0134] Calculate the distribution coefficient of the first transformer, denoted as β (1) β (1) The calculation formula is: (1) =lr (o-1) / lr;
[0135] Calculate the voltage drop from the initial transformer to the first transformer, denoted as ΔU (o-1) ; ΔU (o-1) The calculation formula is: ΔU (o-1) =(β (1) )×(Un(o)-Ui1);
[0136] Step S3322: Obtain the actual receiving voltage of the first transformer, recorded as DU (1) , actual power is recorded as DP (1) ; Get the current transformation ratio of the first transformer and record it as DQ (1) ;
[0137] Calculate the actual receiving current of the first transformer, recorded as DI (1) DI (1) The calculation formula is: DI (1) =DP (1) / DU (1) ;
[0138] The actual reactance of the first transformer is denoted by DX (1) ,
[0139] DX (1) The calculation formula is: DX (1) =[ΔU (o-1) ] / DI (1) ;
[0140] The actual susceptance of the first transformer is denoted as DB (1) ;
[0141] DB (1) The calculation formula is: DB (1) =[3×(DI (1) ) 2 ] 1 / 2 / DU (1) ;
[0142] The actual resistance of the first transformer is denoted as DR (1) ;
[0143] DR (1) The calculation formula is: DR (1) =[PG(o)×(DU (1) ) 2 ] / [mn×(DP (1) ) 2 ];
[0144] The actual conductance of the first transformer is denoted as DG (1) ;
[0145] DG (1) The calculation formula is: DG (1) =(DP (1) ) / [(DU (1) ) 2 ];
[0146] Step S3323: Calculate DX (1)The offset coefficient with equivalent susceptance is denoted as zx (1) ;
[0147] zx (1) The calculation formula is: zx (1) =1-[DX (1) / (β (1) ×XT1)];
[0148] Calculate DB (1) The deviation coefficient from the equivalent reactance is denoted as zb (1) ;
[0149] zb (1) The calculation formula is: (1) =1-[DB (1) / (β (1) ×BT1)];
[0150] Calculating DR (1) The deviation coefficient from the equivalent reactance is denoted as zr (1) ;
[0151] z (1) The calculation formula is: (1) =1-[DR (1) / (β (1) ×RT1)];
[0152] Calculate DG (1) The deviation coefficient from the equivalent reactance is denoted as zg (1) ;
[0153] z (1) The calculation formula is: (1) =1-[β (1) / (DG (1) ×GT1)];
[0154] Step S3324: Calculate zx (1) +zb (1) +zr (1) +zg (1) The sum of is denoted as zz (1) ;
[0155] Assume that the transformation ratio of the first transformer after adjustment is QQ (1) , judge zz (1) Size DeterminationQQ (1) Calculation formula of
[0156] Ifzz (1) ≥0, then QQ (1) The calculation formula is: =[(2-zz (1) ) 1 / 2 ]×DQ (1) ;
[0157] Ifzz (1) <0, then QQ (1) The calculation formula is: =[(1-zz (1) ) 1 / 2 ]×DQ (1) ;
[0158] Among them, QQ (1) is a positive integer or a rational number.
[0159] Furthermore, the specific steps of step S334 are as follows:
[0160] Step S3341: Obtain the length of the transmission line from the first transformer to the second transformer, denoted as lr (1-2) ;
[0161] Calculate the distribution coefficient of the second transformer, denoted as β (2) β (2) The calculation formula is: (2) =lr (1-2) / lr;
[0162] Calculate the voltage drop from the first to the second transformer, denoted as ΔU (1-2) ; ΔU (1-2) The calculation formula is: ΔU (1-2) =(β (2) )×(Un(o)-Ui1);
[0163] Step S3342: Obtain the actual receiving voltage of the second transformer, recorded as DU (2 , actual power is recorded as DP (2) ; Get the current transformation ratio of the second transformer and record it as DQ (2) ;
[0164] Calculate the actual receiving current of the first transformer, recorded as DI (2) DI (2) The calculation formula is: DI (2) =DP (2) / DU (2) ;
[0165] The actual reactance of the second transformer is denoted by DX (2) ,
[0166] DX (2) The calculation formula is: DX (2) =[ΔU (1-2) ] / DI (2) ;
[0167] The actual susceptance of the second transformer is denoted as DB (2) ;
[0168] DB (2) The calculation formula is: DB (2) =[3×(DI (2) ) 2 ] 1 / 2 / DU (2) ;
[0169] The actual resistance of the second transformer is denoted as DR (2) ;
[0170] DR (2) The calculation formula is: DR (2) =[PG(o)×(DU (2) ) 2 ] / [mn×(DP (2) ) 2 ];
[0171] The actual conductance of the second transformer is denoted as DG (1) ;
[0172] DG (2) The calculation formula is: DG (2) =(DP (2) ) / [(DU (2) ) 2 ];
[0173] Step S3343: Calculate DX (2) The offset coefficient with equivalent susceptance is denoted as zx (2) ;
[0174] zx (2) The calculation formula is: zx (2) =1-[DX (2) / (β (2) ×XT1)];
[0175] Calculate DB (1) The deviation coefficient from the equivalent reactance is denoted as zb (2) ;
[0176] zb (2) The calculation formula is: (2) =1-[DB (2) / (β (2) ×BT1)];
[0177] Calculating DR (2) The deviation coefficient from the equivalent reactance is denoted as zr (2) ;
[0178] z (2) The calculation formula is: (2) =1-[DR (2) / (β (2)×RT1)];
[0179] Calculate DG (1) The deviation coefficient from the equivalent reactance is denoted as zg (2) ;
[0180] z (2) The calculation formula is: (2) =1-[β (2) / (DG (2) ×GT1)];
[0181] Step S3344: Calculate zx (2) +zb (2) +zr (2) +zg (2) The sum of is denoted as zz (2) ;
[0182] Assume that the transformation ratio of the second transformer after adjustment is QQ (2) , judge zz (2) Size DeterminationQQ (2) Calculation formula of
[0183] Ifzz (2) ≥0, then QQ (2) The calculation formula is: =[(2-zz (2) ) 1 / 2 ]×DQ (2) ;
[0184] Ifzz (2) <0, then QQ (1) The calculation formula is: =[(1-zz (2) ) 1 / 2 ]×DQ (2) ;
[0185] Among them, QQ (2) is a positive integer or a rational number.
[0186] A power big data analysis system based on cloud computing includes:
[0187] Power consumption analysis module: used to obtain the number of all urban districts in the target area; obtain the number of residents, current power consumption and historical power consumption of each urban district to obtain power consumption data; calculate the priority coefficient of each urban district based on the power consumption data to obtain the allocation sequence;
[0188] Transmission analysis module: used to obtain the no-load electromotive force, rated voltage and reactance of the power generation end, obtain the receiving voltage and receiving power of each municipal district, and obtain the transmission data; construct the mathematical relationship between the power generation end, the power transmission end and the power receiving end according to the transmission data, and calculate the equivalent power parameter value of the power supply line in each municipal district according to the mathematical relationship;
[0189] Power transmission control module: used to obtain the series-parallel relationship of the transformers between each power supply line and obtain the connection data of the transformers; based on the allocation sequence, according to the connection data and equivalent power parameter values of the transformers, the transformer ratio of each power supply line is adjusted to allocate power to each municipal area;
[0190] Continuous update module: used to continuously obtain the current power consumption of each city district, update and adjust the power distribution of each city district.
[0191] Compared with the prior art, the present invention has the following beneficial effects:
[0192] Large-scale data processing capability: The present invention can provide large-scale data processing for the power system; through data statistics, fitting and calculation capabilities, the present invention analyzes and adjusts the data at the power generation end, transmission end and power consumption end in the power system, which can meet the power system's needs for big data processing and power distribution.
[0193] Intelligent analysis and prediction
[0194] Optimizing power supply: The present invention helps the power system predict load and optimize scheduling by analyzing and summarizing data in the power system, thereby improving the stability and operating efficiency of the system; power companies can understand the power demand at any time and the power supply obtained from different sources, optimize power supply capacity, reduce power supply costs and enhance power supply reliability.
[0195] Improve the efficiency of power resource utilization: The present invention uses data processing technology to build a resource visualization platform for power enterprises and realize unified management and allocation of power resources; at the same time, the present invention adds a lot of automation content, which is friendly to industry practitioners. BRIEF DESCRIPTION OF THE DRAWINGS
[0196] Other features, objects and advantages of the present invention will become more apparent from the detailed description of non-limiting embodiments made with reference to the following drawings:
[0197] Figure 1 It is a schematic diagram of the method of the present invention;
[0198] Figure 2 It is a schematic diagram of the system of the present invention;
[0199] Figure 3 It is an equivalent schematic diagram of the power transmission process of the present invention;
[0200] Figure 4 This is a schematic diagram of the power supply circuit of the present invention. DETAILED DESCRIPTION
[0201] 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.
[0202] Embodiment 1
[0203] See also Figure 1 , a power big data analysis method based on cloud computing includes:
[0204] Step S1: Obtain the number of all municipal districts in the target area; obtain the number of residents, current power consumption and historical power consumption of each municipal district (in the target area) to obtain power consumption data; calculate the priority coefficient of each municipal district based on the power consumption data to obtain an allocation sequence;
[0205] It should be noted that the “target area” in the present invention means: the city-level area where the present invention (a power big data analysis method and system based on cloud computing) is used to perform power analysis;
[0206] The specific steps of step S1 are as follows:
[0207] Step S11: Obtain the number of all city districts in the target area, denoted as mn;
[0208] Get the number of residents corresponding to each city district and record it as pn1, pn2~pn mn , the current power consumption (unit: degree or kilowatt-hour) is recorded as oj1, oj2~oj mn ;
[0209] Among them, pn1 and oj1 represent the number of residents and current power consumption in the first city district respectively; pn2 and oj2 represent the number of residents and current power consumption in the second city district respectively; and so on. mn and oj mn They represent the number of residents and current power consumption in the mnth municipal district respectively;
[0210] Step S12: Extract oj1~oj mn The maximum value is recorded as ojm, the minimum value is recorded as ojl, and the mode is recorded as ojz; calculate oj1~oj mn The average value of is denoted as aoj;
[0211] The relay power is recorded as ojc; compare ojz with aoj to determine the value of ojc;
[0212] If ojz>aoj, then the value of ojc is ojz;
[0213] If ojz≤aoj, then the value of ojc is aoj;
[0214] Step S13: Calculate the difference between ojm and ojl, recorded as ooj; (i.e., ojm-ojl=ooj)
[0215] Assume that the current power consumption of the ith district (in the target area) is oj i , the value range of i is: 1~mn;
[0216] The priority allocation coefficient corresponding to the i-th municipal district is denoted as α(i); the calculation formula of α(i) is defined as:
[0217] Judge oj i Whether it is within the interval [ojc, ojm];
[0218] If , then the calculation formula of α(i) is as follows:
[0219] α(i)=[(oj i -ojc) / (ojm-ojc)]×0.5;
[0220] If not, the calculation formula of α(i) is as follows:
[0221] α(i)=[(oj i -ojl) / (ojc-ojl)]×0.5;
[0222] Step S14: oj1~oj mn Substituting into the calculation formula of α(i), we can obtain the priority allocation coefficients corresponding to the 1st to mnth municipal districts (within the target area): α(1), α(2)~α(mn);
[0223] Among them, α(1) represents the priority allocation coefficient of the first municipal district; α(2) represents the priority allocation coefficient of the second municipal district; and so on, α(mn) represents the priority allocation coefficient of the mnth municipal district;
[0224] Obtain the electricity consumption of each municipal district (within the target area) in the past three years as historical electricity consumption;
[0225] The electricity consumption of the first city district in three years is: ahj(1), bhj(1) and chj(1); among them, ahj(1) represents the electricity consumption of the first city district in the third year, bhj(1) represents the electricity consumption of the first city district in the second year, and chj(1) represents the electricity consumption of the first city district in the first year;
[0226] The electricity consumption of the second city district in three years is: ahj(2), bhj(2) and chj(2); among them, ahj(2) represents the electricity consumption of the second city district in the third year, bhj(2) represents the electricity consumption of the second city district in the second year, and chj(2) represents the electricity consumption of the second city district in the first year;
[0227] By analogy, the power consumption of the mn-th municipal district in three years is: ahj(mn), bhj(mn) and chj(mn); among them, ahj(mn) represents the power consumption of the mn-th municipal district in the third year, bhj(mn) represents the power consumption of the mn-th municipal district in the second year, and chj(mn) represents the power consumption of the mn-th municipal district in the first year;
[0228] Step S15: Calculate the electricity consumption change rate of each city district (within the target area) based on the historical electricity consumption to obtain ub(1), ub(2)~ub(mn);
[0229] Among them, ub(1) represents the electricity consumption change rate of the first district in the target area; ub(2) represents the electricity consumption change rate of the second district in the target area; and so on, ub(mn) represents the electricity consumption change rate of the mnth district in the target area;
[0230] Step S151: Calculate the power consumption change rate of the first city district according to ahj(1), bhj(1) and chj(1), denoted as ub(1);
[0231] 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) =[bhj(1) / chj(1)]-1;
[0232] 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) ; up (2-1) The calculation formula is: up (2-1) =[ahj(1) / bhj(1)]-1;
[0233] Step S152: According to up (3-2) 、up (2-1) , ahj(1) and bhj(1), construct the state transfer matrix, denoted as matrix B; the mathematical expression of matrix B is as follows:
[0234] ; Where n represents the number of iterations of matrix B, × represents matrix multiplication; n is a natural number, and its initial value is 1; each time matrix B is iterated, n is increased by 1 on the original basis;
[0235] Use the NumPy library function to iterate the matrix B and obtain the steady-state transfer matrix, which is recorded as matrix B1;
[0236] Define a 1×2 order receiving matrix {1, 0}, 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 element in the matrix B3 as the estimated power consumption of the first municipal district, denoted as jj(1);
[0237] Step S153: Calculate the power consumption change rate ub(1) of the first city district, where ub(1) is calculated as follows: ub(1)=[jj(1) / ahj(1)]-1;
[0238] Step S154: Repeat the same steps of calculating the power consumption change rate of the first city district (i.e., steps S151 to S153), calculate the power consumption change rates of the second to mn-th city districts, and obtain ub(2) to ub(mn);
[0239] Step S16: Calculate the priority coefficient of the first municipal district, denoted as ff(1); the calculation formula of ff(1) is: ff(1)=pn1×α(1)×ub(1);
[0240] Calculate the priority coefficient of the second municipal district, denoted as ff(2); the calculation formula of ff(2) is: ff(2) = pn2×α(2)×ub(2);
[0241] Similarly, the priority coefficient of the mnth district is calculated, denoted as ff(mn); the calculation formula of ff(mn) is: ff(mn) = pn mn ×α(mn)×ub(mn);
[0242] The priority coefficients ff(1) to ff(mn) corresponding to the 1st to mnth municipal districts are summarized to obtain the allocation sequence of the municipal districts, and then the process goes to step S2.
[0243] Step S2: Obtain the no-load electromotive force, rated voltage and reactance of the generating end, obtain the receiving voltage and receiving power of each municipal district (within the target area), and obtain the transmission data; construct a mathematical relationship between the generating end, the transmission end and the receiving end according to the transmission data, and calculate the equivalent power parameter value of the power supply line of each municipal district (within the target area) according to the mathematical relationship;
[0244] It should be noted that the "generating end" in the present invention means: a large-scale generating set or a large-scale power grid provided by a power bureau or a power plant for supplying power to a target area; the "power supply line" means: the shortest power line from the generating end to the city-level district receiving power from the transformer; if there are multiple "shortest power lines" in the city-level district, the city-level district is: the "shortest power line with the most transformers" from the generating end to the city-level district receiving power from the transformer;
[0245] The specific steps of step S2 are as follows:
[0246] Step S21: Obtain the no-load electromotive force of the generating end, recorded as Eq(o), the rated voltage as Un(o), and the reactance as Xd(o); obtain the current frequency (AC) of the generating end, recorded as f;
[0247] Calculate the power generation period, recorded as T; the calculation formula of T is: T = 1 / f;
[0248] Step S22: Get the current time, recorded as dt; calculate the remainder of dt and T, recorded as cd; the calculation formula of cd is: cd = dt%T; where % represents the modulus symbol;
[0249] Taking the current as the positive direction, calculate the phase angle between voltage and current, which is recorded as θ. The calculation formula of θ is: θ = (cd / T) × 2π; where π represents pi;
[0250] Calculate the active power of the power generation equipment, recorded as PG (o); the calculation formula of PG (o) is:
[0251] PG(o)=[Eq(o)×Un(o)] / Xd(o)×sin(θ);
[0252] Step S23: Obtain the transformer voltage connected to the power grid corresponding to each municipal district (within the target area) as the received voltage; and the transformer power as the received power;
[0253] The receiving voltages corresponding to the 1st to mnth municipal districts are recorded as: Ui1, Ui2~Ui mn ;
[0254] The received power is recorded as: Pi1, Pi2~Pi mn ;
[0255] Among them, Ui1 and Pi1 represent the receiving voltage and receiving power of the first city district respectively; Ui2 and Pi2 represent the receiving voltage and receiving power of the second city district respectively; and so on, Ui mn and Pi mn represent the receiving voltage and receiving power of the mn-th municipal district respectively;
[0256] Step S24: Please refer to Figure 3 , according to the data in step S21 to step S23, the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance from the power generation end to each municipal district are calculated in sequence;
[0257] Step S241: the power generation end reaches the first municipal district;
[0258] Step S2411: Calculate the receiving current of the first city district, denoted as Ii1;
[0259] The calculation formula of Ii1 is: Ii1=Pi1 / Ui1;
[0260] Calculate the equivalent reactance from the generator to the first municipal district, denoted as XT1,
[0261] The calculation formula of XT1 is: XT1 = (Un (o) - Ui1) / Ii1;
[0262] Step S2412: Calculate the equivalent susceptance from the power generation end to the first municipal district, denoted as BT1;
[0263] The calculation formula of BT1 is: BT1 = [3 × (Ii1) 2 ] 1 / 2 / Ui1;
[0264] Calculate the equivalent resistance from the generator end to the first municipal district, denoted as RT1;
[0265] The calculation formula of RT1 is: RT1 = [PG (o) × (Ui1) 2 ] / [(Pi1) 2 ];
[0266] Calculate the equivalent conductance from the power generation end to the first municipal district, denoted as GT1;
[0267] The calculation formula of GT1 is: GT1 = (Pi1) / [(Ui1) 2 ];
[0268] Step S242: the power generation end reaches the second municipal district;
[0269] Step S2421: Calculate the receiving current of the second city district, recorded as Ii2;
[0270] The calculation formula of Ii2 is: Ii2=Pi2 / Ui2;
[0271] Calculate the equivalent reactance from the generator to the second municipal district, denoted as XT2,
[0272] The calculation formula of XT2 is: XT2 = (Un (o) - Ui2) / Ii2;
[0273] Step S2422: Calculate the equivalent susceptance from the power generation end to the second municipal district, denoted as BT2;
[0274] The calculation formula of BT2 is: BT2 = [3 × (Ii2) 2 ] 2 / 2 / Ui2;
[0275] Calculate the equivalent resistance from the generator end to the second municipal district, recorded as RT2;
[0276] The calculation formula of RT2 is: RT2 = [PG (o) × (Ui2) 2] / [(Pi2) 2 ];
[0277] Calculate the equivalent conductance from the power generation end to the second municipal district, denoted as GT2;
[0278] The calculation formula of GT2 is: GT2 = (Pi2) / [(Ui2) 2 ];
[0279] Step S243: Similarly, the power generation end reaches the mnth municipal district;
[0280] Step S2431: Calculate the receiving current of the mnth city district, denoted as Ii mn ;
[0281] Ii mn The calculation formula is: Ii mn =Pi mn / Ui mn ;
[0282] Calculate the equivalent reactance from the generator to the mnth municipal district, denoted as XT mn ,
[0283] XT mn The calculation formula is: XT mn = (Un (o) - Ui mn ) / Ii mn ;
[0284] Step S2432: Calculate the equivalent susceptance from the power generation end to the mnth municipal district, denoted as BT mn ;
[0285] BT mn The calculation formula is: BT mn =[3×(Ii mn ) 2 ] mn / 2 / Ui mn ;
[0286] Calculate the equivalent resistance from the generator end to the mnth municipal district, denoted as RT mn ;
[0287] RT mn The calculation formula is: RT mn =[PG(o)×(Ui mn ) 2 ] / [(Pi mn ) 2 ];
[0288] Calculate the equivalent conductance from the power generation end to the mnth municipal district, denoted as GT mn ;
[0289] GT mn The calculation formula is: GT mn =(Pi mn ) / [(Ui mn ) 2 ];
[0290] Step S25: Summarize the equivalent reactances corresponding to the 1st to mnth municipal districts: XT1, XT2~XT mn ; Equivalent susceptance: BT1, BT2~BT mn , equivalent resistance: RT1, RT2~RT mn , equivalent conductance: GT1, GT2~GT mn , as the equivalent power parameter value, enter step S3;
[0291] Among them, XT1, BT1, RT1 and GT1 represent the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance of the first city district, respectively; XT2, BT2, RT2 and GT2 represent the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance of the second city district, respectively; and so on, XT mn , BT mn , RT mn and GT mn , respectively represent the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance of the mnth municipal district.
[0292] Step S3: Obtain the series-parallel relationship of the transformers between each power supply line (in the target area) to obtain the connection data of the transformers; based on the allocation sequence, adjust the total transformer ratio of each power supply line (in the target area) according to the connection data and equivalent power parameter values of the transformers, and allocate power to each municipal district (in the target area);
[0293] The specific steps of step S3 are as follows:
[0294] Step S31: constructing a power transmission relationship map for each city district according to the equivalent power parameter value;
[0295] The first municipal district: {PG (o), Un (o)}—{XT1, BT1, RT1, GT1}—{Pi1, Ui1};
[0296] The second municipal district: {PG (o), Un (o)}—{XT2, BT2, RT2, GT2}—{Pi2, Ui2};
[0297] And so on, the mnth district: {PG(o), Un(o)}—{XT mn , BT mn , RT mn , GT mn}—{Pimn , Ui mn};
[0298] Step S32: Obtain the priority coefficients corresponding to the 1st to mnth municipal districts: ff(1), ff(2) to ff(mn);
[0299] Calculate the value of ff(1):ff(2):~:ff(mn), recorded as ap1:ap2:~:ap mn ;(i.e. ff(1):ff(2):~:ff(mn)=ap1:ap2:~:ap mn ) Among them, ap1~ap mn For positive numbers, round down;
[0300] Extraction of ap1~ap mn The minimum value of apl is recorded as apl, the received power corresponding to apl is recorded as Pil, and the initial transformer ratio from the power generation end to each municipal district is calculated;
[0301] Step S321: Calculate the standard transformation ratio of the initial transformer, denoted as fb;
[0302] The calculation formula of fb is: fb = [PG (o) / Pil] 1 / 2 ;
[0303] Step S322: Calculate the initial transformer ratio for power transmission from the power generation end to the first municipal district, denoted as B (o-1) ; B (o-1) The calculation formula is: (o-1) = (apl / ap1) × fb;
[0304] Calculate the initial transformer ratio for power transmission from the generator to the second municipal district, denoted as B (o-2) ; B (o-2) The calculation formula is: (o-2) = (apl / ap2) × fb;
[0305] By analogy, calculate the initial transformer ratio for power transmission from the power generation end to the mnth municipal district, denoted as B (o-mn) ; B (o-mn) The calculation formula is: (o-mn) = (apl / ap mn )×fb;
[0306] Step S33: Please refer to Figure 4 , according to the transmission relationship mapping of the first municipal district, adjust the transformer ratio of the first power supply line;
[0307] Step S331: For the first power supply line, obtain the number of transformers other than the initial transformer and the receiving transformer, recorded as tr; obtain the length of the power transmission line from the initial transformer to the receiving transformer, recorded as lr;
[0308] Step S332: the first power supply line, the initial transformer to the first transformer;
[0309] Step S3321: Obtain the length of the transmission line from the initial transformer to the first transformer, denoted as lr (o-1) ;
[0310] Calculate the distribution coefficient of the first transformer, denoted as β (1) β (1) The calculation formula is: (1) =lr (o-1) / lr;
[0311] Calculate the voltage drop from the initial transformer to the first transformer, denoted as ΔU (o-1) ; ΔU (o-1) The calculation formula is: ΔU (o-1) =(β (1) )×(Un(o)-Ui1);
[0312] Step S3322: Obtain the actual receiving voltage of the first transformer, recorded as DU (1) , actual power is recorded as DP (1) ; Get the current transformation ratio of the first transformer and record it as DQ (1) ;
[0313] Calculate the actual receiving current of the first transformer, recorded as DI (1) DI (1) The calculation formula is: DI (1) =DP (1) / DU (1) ;
[0314] The actual reactance of the first transformer is denoted by DX (1) ,
[0315] DX (1) The calculation formula is: DX (1) =[ΔU (o-1) ] / DI (1) ;
[0316] The actual susceptance of the first transformer is denoted as DB (1) ;
[0317] DB (1) The calculation formula is: DB (1) =[3×(DI (1) ) 2 ]1 / 2 / DU (1) ;
[0318] The actual resistance of the first transformer is denoted as DR (1) ;
[0319] DR (1) The calculation formula is: DR (1) =[PG(o)×(DU (1) ) 2 ] / [mn×(DP (1) ) 2 ];
[0320] The actual conductance of the first transformer is denoted as DG (1) ;
[0321] DG (1) The calculation formula is: DG (1) =(DP (1) ) / [(DU (1) ) 2 ];
[0322] Step S3323: Calculate DX (1) The offset coefficient with equivalent susceptance is denoted as zx (1) ;
[0323] zx (1) The calculation formula is: zx (1) =1-[DX (1) / (β (1) ×XT1)];
[0324] Calculate DB (1) The deviation coefficient from the equivalent reactance is denoted as zb (1) ;
[0325] zb (1) The calculation formula is: (1) =1-[DB (1) / (β (1) ×BT1)];
[0326] Calculating DR (1) The deviation coefficient from the equivalent reactance is denoted as zr (1) ;
[0327] z (1) The calculation formula is: (1) =1-[DR (1) / (β (1) ×RT1)];
[0328] Calculate DG (1) The deviation coefficient from the equivalent reactance is denoted as zg (1) ;
[0329] z (1) The calculation formula is: (1) =1-[β (1) / (DG (1) ×GT1)];
[0330] Step S3324: Calculate zx (1) +zb (1) +zr (1) +zg (1) The sum of is denoted as zz (1) ;
[0331] Assume that the transformation ratio of the first transformer after adjustment is QQ (1) , judge zz (1) Size DeterminationQQ (1) Calculation formula of
[0332] Ifzz (1) ≥0, then QQ (1) The calculation formula is: =[(2-zz (1) ) 1 / 2 ]×DQ (1) ;
[0333] Ifzz (1) <0, then QQ (1) The calculation formula is: =[(1-zz (1) ) 1 / 2 ]×DQ (1) ;
[0334] Among them, QQ (1) is a positive integer or a rational number;
[0335] Step S333: The first power supply line, the first transformer to the second transformer;
[0336] Determine whether the second transformer is connected in parallel with other transformers;
[0337] If not connected in parallel, proceed to step S334;
[0338] If in parallel, then obtain the transformer information connected in parallel with the second transformer and proceed to step S335;
[0339] Step S334: the second transformer is not connected in parallel with other transformers, and the adjusted transformation ratio of the second transformer is calculated;
[0340] Step S3341: Obtain the length of the transmission line from the first transformer to the second transformer, denoted as lr (1-2) ;
[0341] Calculate the distribution coefficient of the second transformer, denoted as β (2)β (2) The calculation formula is: (2) =lr (1-2) / lr;
[0342] Calculate the voltage drop from the first to the second transformer, denoted as ΔU (1-2) ; ΔU (1-2) The calculation formula is: ΔU (1-2) =(β (2) )×(Un(o)-Ui1);
[0343] Step S3342: Obtain the actual receiving voltage of the second transformer, recorded as DU (2 , actual power is recorded as DP (2) ; Get the current transformation ratio of the second transformer and record it as DQ (2) ;
[0344] Calculate the actual receiving current of the first transformer, recorded as DI (2) DI (2) The calculation formula is: DI (2) =DP (2) / DU (2) ;
[0345] The actual reactance of the second transformer is denoted by DX (2) ,
[0346] DX (2) The calculation formula is: DX (2) =[ΔU (1-2) ] / DI (2) ;
[0347] The actual susceptance of the second transformer is denoted as DB (2) ;
[0348] DB (2) The calculation formula is: DB (2) =[3×(DI (2) ) 2 ] 1 / 2 / DU (2) ;
[0349] The actual resistance of the second transformer is denoted as DR (2) ;
[0350] DR (2) The calculation formula is: DR (2) =[PG(o)×(DU (2) ) 2 ] / [mn×(DP (2) ) 2 ];
[0351] The actual conductance of the second transformer is denoted as DG(1) ;
[0352] DG (2) The calculation formula is: DG (2) =(DP (2) ) / [(DU (2) ) 2 ];
[0353] Step S3343: Calculate DX (2) The offset coefficient with equivalent susceptance is denoted as zx (2) ;
[0354] zx (2) The calculation formula is: zx (2) =1-[DX (2) / (β (2) ×XT1)];
[0355] Calculate DB (1) The deviation coefficient from the equivalent reactance is denoted as zb (2) ;
[0356] zb (2) The calculation formula is: (2) =1-[DB (2) / (β (2) ×BT1)];
[0357] Calculating DR (2) The deviation coefficient from the equivalent reactance is denoted as zr (2) ;
[0358] z (2) The calculation formula is: (2) =1-[DR (2) / (β (2) ×RT1)];
[0359] Calculate DG (1) The deviation coefficient from the equivalent reactance is denoted as zg (2) ;
[0360] z (2) The calculation formula is: (2) =1-[β (2) / (DG (2) ×GT1)];
[0361] Step S3344: Calculate zx (2) +zb (2) +zr (2) +zg (2) The sum of is denoted as zz (2) ;
[0362] Assume that the transformation ratio of the second transformer after adjustment is QQ(2) , judge zz (2) Size DeterminationQQ (2) Calculation formula of
[0363] Ifzz (2) ≥0, then QQ (2) The calculation formula is: =[(2-zz (2) ) 1 / 2 ]×DQ (2) ;
[0364] Ifzz (2) <0, then QQ (1) The calculation formula is: =[(1-zz (2) ) 1 / 2 ]×DQ (2) ;
[0365] Among them, QQ (2) is a positive integer or a rational number;
[0366] Step S335: the second transformer is connected in parallel with other transformers, and the adjusted transformation ratio of the second transformer is calculated;
[0367] Step S3351: other transformers connected in parallel with the two transformers are used as parallel transformers;
[0368] Count the number of transformers in parallel, denoted as oh;
[0369] Get the current transformation ratio of the 1st to ohth parameter transformer and record it as: OB1, OB2~OB oh ;
[0370] Among them, OB1 represents the current transformation ratio of the first parallel transformer; OB2 represents the current transformation ratio of the second parallel transformer; and so on. oh Indicates the current transformation ratio of the ohth parallel transformer;
[0371] Step S3352: Using the transformation ratio DQ of the second transformer (2) As the first term, calculate DQ (2) :OB1:OB2:~:OB oh The ratio is recorded as: uu1:uu2:~:uu oh :uu (oh+1) ; (ie DQ (2) :OB1:OB2:~:OB oh =uu1:uu2:~:uu oh :uu (oh+1) )
[0372] Calculate uu1:uu2:~:uu oh :uu (oh+1)The inverse ratio is expressed as: nn1:nn2:~:nn oh :nn (oh+1) ;
[0373] Step S3353: Repeat the same steps (steps S3341 to S3344) of calculating the adjusted transformation ratio of the second transformer on the premise that the second transformer is not connected in parallel with other transformers, and determine QQ (2) The value of QQ (2) Indicates the transformation ratio of the second transformer after adjustment;
[0374] Calculate nn1 + nn2 + ~ + nn oh +nn (oh+1) ; and is written as ann;
[0375] The transformation ratio of the second transformer after secondary adjustment is denoted as cQQ (2) ;cQQ (2) The calculation formula is:
[0376] QQQ (2) =(QQ (2) )×(nn1 / ann);
[0377] Step S336: repeat the same steps of step S332 to step S335, and adjust the transformation ratios of the 3rd to trth transformers in sequence;
[0378] Step S34: Repeat the same steps of adjusting the transformation ratio of the transformer of the first power supply line to adjust the transformation ratios of the transformers of the second to mnth power supply lines.
[0379] It should be noted that, because the transformer is an inductive element, the theoretical sodium under natural conditions after multiple groups of transformers are connected in series and parallel in "Electrical Machinery" should be 0; but in the actual long-distance power supply process, because the power supply cable is thicker and longer, the self-inductance and mutual inductance generated by the cable during power transmission cannot be ignored, resulting in a decrease in the electric energy received by the transformer; therefore, in the calculation process of the present invention, the electrical susceptance generated by the self-inductance and mutual inductance of the cable is aggregated into the transformer and regarded as "generated" by the transformer; this processing method is only for the convenience of calculation, which is inconsistent with objective physical facts, and is specially explained again.
[0380] Step S4: continuously obtain the current power consumption of each city district (within the target area), and update and adjust the power distribution of each city district (within the target area).
[0381] Embodiment 2
[0382] See also Figure 2A power big data analysis system based on cloud computing includes: a power consumption analysis module, a power transmission analysis module, a power transmission control module, a continuous update module, a database and a server; wherein the power consumption analysis module, the power transmission analysis module, the power transmission control module and the continuous update module are connected to the database and the server respectively.
[0383] Power consumption analysis module: used to obtain the number of all municipal districts in the target area; obtain the number of residents, current power consumption and historical power consumption of each municipal district (in the target area) to obtain power consumption data; calculate the priority coefficient of each municipal district based on the power consumption data to obtain the allocation sequence;
[0384] Transmission analysis module: used to obtain the no-load electromotive force, rated voltage and reactance of the power generation end, obtain the receiving voltage and receiving power of each municipal district (within the target area), and obtain the transmission data; construct the mathematical relationship between the power generation end, the power transmission end and the power receiving end according to the transmission data, and calculate the equivalent power parameter value of the power supply line of each municipal district (within the target area) according to the mathematical relationship;
[0385] Power transmission control module: used to obtain the series-parallel relationship of transformers between each power supply line (in the target area) and obtain the connection data of the transformers; based on the allocation sequence, adjust the total transformer ratio of each power supply line (in the target area) according to the connection data and equivalent power parameter values of the transformers, and allocate power to each municipal district (in the target area);
[0386] Continuous update module: used to continuously obtain the current power consumption of each city district (within the target area), update and adjust the power distribution of each city district (within the target area).
[0387] The above formulas are all dimensionless and numerical calculations. 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 technical personnel in this field according to actual conditions. For example, if there are weight coefficients and proportional coefficients, their set sizes are 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, it is sufficient as long as it does not affect the proportional relationship between the parameter and the quantized value.
[0388] 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 protection scope of the present invention is not limited thereto. Although the present invention is described in detail with reference to the above-described embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-described embodiments within the technical scope disclosed by the present invention, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.
Claims
1. A power big data analysis method based on cloud computing, characterized in that: The method comprises: Step S1: Obtain the number of all municipal districts in the target area; obtain the number of residents, current power consumption and historical power consumption of each municipal district to obtain power consumption data; calculate the priority coefficient of each municipal district based on the power consumption data to obtain an allocation sequence; Step S2: Obtain the no-load electromotive force, rated voltage and reactance of the power generation end, obtain the receiving voltage and receiving power of each city district, and obtain the power transmission data; construct a mathematical relationship between the power generation end, the power transmission end and the power receiving end according to the power transmission data, and calculate the equivalent power parameter value of the power supply line of each city district according to the mathematical relationship; Step S3: Obtain the series-parallel relationship of the transformers between each power supply line to obtain the connection data of the transformers; based on the allocation sequence, adjust the transformer ratio of each power supply line according to the connection data and equivalent power parameter value of the transformer, and allocate power to each city district; Step S4: continuously obtain the current power consumption of each city district, and update and adjust the power distribution of each city district; The specific steps of step S3 are as follows: Step S31: constructing a power transmission relationship map for each city district according to the equivalent power parameter value; The first municipal district: {PG (o), Un (o)}—{XT1, BT1, RT1, GT1}—{Pi1, Ui1}; The second municipal district: {PG (o), Un (o)}—{XT2, BT2, RT2, GT2}—{Pi2, Ui2}; And so on, the mnth district: {PG(o), Un(o)}—{XT mn , BT mn , RT mn , GT mn }—{Pi mn , Ui mn }; Step S32: Obtain the priority coefficients corresponding to the 1st to mnth municipal districts: ff(1), ff(2) to ff(mn); Calculate the value of ff(1):ff(2):~:ff(mn), recorded as ap1:ap2:~:ap mn ; Among them, ap1~ap mn For positive numbers, round down; Extraction of ap1~ap mn The minimum value of apl is recorded as apl, the received power corresponding to apl is recorded as Pil, and the initial transformer ratio from the power generation end to each municipal district is calculated; Step S321: Calculate the standard transformation ratio of the initial transformer, denoted as fb; The calculation formula of fb is: fb = [PG (o) / Pil] 1 / 2 ; Step S322: Calculate the initial transformer ratio for power transmission from the power generation end to the first municipal district, denoted as B (o-1) ; B (o-1) The calculation formula is: (o-1) = (apl / ap1) × fb; Calculate the initial transformer ratio for power transmission from the generator to the second municipal district, denoted as B (o-2) ; B (o-2) The calculation formula is: (o-2) = (apl / ap2) × fb; By analogy, calculate the initial transformer ratio for power transmission from the power generation end to the mnth municipal district, denoted as B (o-mn) ; B (o-mn) The calculation formula is: (o-mn) = (apl / ap mn )×fb; Step S33: adjusting the transformer ratio of the first power supply line according to the power transmission relationship mapping of the first municipal district; Step S34: Repeat the same steps of adjusting the transformation ratio of the transformer of the first power supply line to adjust the transformation ratios of the transformers of the second to mnth power supply lines.
2. The power big data analysis method based on cloud computing according to claim 1 is characterized in that: The specific steps of step S1 are as follows: Step S11: Obtain the number of all city districts in the target area, denoted as mn; Get the number of residents corresponding to each city district and record it as pn1, pn2~pn mn , the current power consumption is recorded as oj1, oj2~oj mn ; Among them, pn1 and oj1 represent the number of residents and current power consumption in the first city district respectively; pn2 and oj2 represent the number of residents and current power consumption in the second city district respectively; and so on. mn and oj mn They represent the number of residents and current power consumption in the mnth municipal district respectively; Step S12: Extract oj1~oj mn The maximum value is recorded as ojm, the minimum value is recorded as ojl, and the mode is recorded as ojz; calculate oj1~oj mn The average value of is denoted as aoj; The relay power is recorded as ojc; compare ojz with aoj to determine the value of ojc; If ojz>aoj, then the value of ojc is ojz; If ojz≤aoj, then the value of ojc is aoj; Step S13: Calculate the difference between ojm and ojl, denoted as ooj; Assume that the current power consumption of the i-th district is oj i , the value range of i is: 1~mn; The priority allocation coefficient corresponding to the i-th municipal district is denoted as α(i); the calculation formula of α(i) is defined as: Judge oj i Whether it is within the interval [ojc, ojm]; If , then the calculation formula of α(i) is as follows: α(i)=[(oj i -ojc) / (ojm-ojc]]×0.5; If not, the calculation formula of α(i) is as follows: α(i)=[(oj i 𝑖𝑖𝑗𝑙) / (𝑖𝑗𝑐-𝑖𝑗𝑙)]×0.5; Step S14: oj1~oj mn Substituting into the calculation formula of α(i), we can obtain the priority allocation coefficients corresponding to the 1st to the mnth municipal districts: α(1), α(2)~α(mn); Among them, α(1) represents the priority allocation coefficient of the first municipal district; α(2) represents the priority allocation coefficient of the second municipal district; and so on, α(mn) represents the priority allocation coefficient of the mnth municipal district; Obtain three years of electricity consumption in each municipal district as historical electricity consumption; The electricity consumption of the first city district in three years is: ahj(1), bhj(1) and chj(1); among them, ahj(1) represents the electricity consumption of the first city district in the third year, bhj(1) represents the electricity consumption of the first city district in the second year, and chj(1) represents the electricity consumption of the first city district in the first year; The electricity consumption of the second city district in three years is: ahj(2), bhj(2) and chj(2); among them, ahj(2) represents the electricity consumption of the second city district in the third year, bhj(2) represents the electricity consumption of the second city district in the second year, and chj(2) represents the electricity consumption of the second city district in the first year; By analogy, the power consumption of the mn-th municipal district in three years is: ahj(mn), bhj(mn) and chj(mn); among them, ahj(mn) represents the power consumption of the mn-th municipal district in the third year, bhj(mn) represents the power consumption of the mn-th municipal district in the second year, and chj(mn) represents the power consumption of the mn-th municipal district in the first year; Step S15: Calculate the electricity consumption change rate of each city district according to the historical power consumption to obtain ub(1), ub(2)~ub(mn); 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 city district; and so on, ub(mn) represents the electricity consumption change rate of the mnth city district; Step S16: Calculate the priority coefficient of the first municipal district, denoted as ff(1); the calculation formula of ff(1) is: ff(1)=pn1×α(1)×ub(1); Calculate the priority coefficient of the second municipal district, denoted as ff(2); the calculation formula of ff(2) is: ff(2) = pn2×α(2)×ub(2); Similarly, the priority coefficient of the mnth district is calculated, denoted as ff(mn); the calculation formula of ff(mn) is: ff(mn) = pn mn ×α(mn)×ub(mn); The priority coefficients ff(1) to ff(mn) corresponding to the 1st to mnth municipal districts are summarized to obtain the allocation sequence of the municipal districts, and then the process goes to step S2.
3. The power big data analysis method based on cloud computing according to claim 1 is characterized in that: The specific steps of step S2 are as follows: Step S21: Obtain the no-load electromotive force of the generating end, recorded as Eq(o), the rated voltage as Un(o), and the reactance as Xd(o); obtain the current frequency of the generating end, recorded as f; Calculate the power generation period, recorded as T; the calculation formula of T is: T = 1 / f; Step S22: Get the current time, recorded as dt; calculate the remainder of dt and T, recorded as cd; the calculation formula of cd is: cd = dt%T; where % represents the modulus symbol; Taking the current as the positive direction, calculate the phase angle between voltage and current, which is recorded as θ. The calculation formula of θ is: θ = (cd / T) × 2π; where π represents pi; Calculate the active power of the power generation equipment, recorded as PG (o); the calculation formula of PG (o) is: PG(o)=[Eq(o)×Un(o)] / Xd(o)×sin(θ); Step S23: Obtain the transformer voltage connected to the power supply grid corresponding to each municipal district as the received voltage; and the transformer power as the received power; The receiving voltages corresponding to the 1st to mnth municipal districts are recorded as: Ui1, Ui2~Ui mn ; The received power is recorded as: Pi1, Pi2~Pi mn ; Among them, Ui1 and Pi1 represent the receiving voltage and receiving power of the first city district respectively; Ui2 and Pi2 represent the receiving voltage and receiving power of the second city district respectively; and so on, Ui mn and Pi mn Respectively represent the receiving voltage and receiving power of the mn-th municipal district; Step S24: according to the data in steps S21 to S23, the equivalent reactance, equivalent susceptance, equivalent resistance and equivalent conductance from the power generation end to each municipal district are calculated in sequence; Step S25: Summarize the equivalent reactances corresponding to the 1st to mnth municipal districts: XT1, XT2~XT mn ; Equivalent susceptance: BT1, BT2~BT mn , equivalent resistance: RT1, RT2~RT mn , equivalent conductance: GT1, GT2~GT mn , enter step S3 as the equivalent power parameter value.
4. According to the cloud computing-based power big data analysis method of claim 2, it is characterized in that: The specific steps of step S15 are as follows: Step S151: Calculate the power consumption change rate of the first city district according to ahj(1), bhj(1) and chj(1), 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) ; up (3-2) The calculation 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) ; up (2-1) The calculation formula is: up (2-1) =[ahj(1) / bhj(1)]-1; Step S152: According to up (3-2) 、up (2-1) , ahj(1) and bhj(1), construct the state transfer matrix, denoted as matrix B; the mathematical expression of matrix B is as follows: ; Where n represents the number of iterations of matrix B, × represents matrix multiplication; n is a natural number, and its initial value is 1; each time matrix B is iterated, n is increased by 1 on the original basis; Use the NumPy library function to iterate the matrix B and obtain the steady-state transfer matrix, which is recorded as matrix B1; Define a 1×2 order receiving matrix {1, 0}, 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 element in the matrix B3 as the estimated power consumption of the first municipal district, denoted as jj(1); Step S153: Calculate the power consumption change rate ub(1) of the first city district, where ub(1) is calculated as follows: ub(1)=[jj(1) / ahj(1)]-1; Step S154: Repeat the same steps of calculating the power consumption change rate of the first urban district, calculate the power consumption change rates of the second to mn-th urban districts, and obtain ub(2)~ub(mn).
5. The power big data analysis method based on cloud computing according to claim 3 is characterized in that: The specific steps of step S24 are as follows: Step S241: the power generation end reaches the first municipal district; Step S2411: Calculate the receiving current of the first city district, denoted as Ii1; The calculation formula of Ii1 is: Ii1=Pi1 / Ui1; Calculate the equivalent reactance from the generator to the first municipal district, denoted as XT1, The calculation formula of XT1 is: XT1 = (Un (o) - Ui1) / Ii1; Step S2412: Calculate the equivalent susceptance from the power generation end to the first municipal district, denoted as BT1; The calculation formula of BT1 is: BT1 = [3 × (Ii1) 2 ] 1 / 2 / Ui1; Calculate the equivalent resistance from the generator end to the first municipal district, denoted as RT1; The calculation formula of RT1 is: RT1 = [PG (o) × (Ui1) 2 ] / [(Pi1) 2 ]; Calculate the equivalent conductance from the power generation end to the first municipal district, denoted as GT1; The calculation formula of GT1 is: GT1 = (Pi1) / [(Ui1) 2 ]; Step S242: the power generation end reaches the second municipal district; Step S2421: Calculate the receiving current of the second city district, denoted as Ii2; The calculation formula of Ii2 is: Ii2=Pi2 / Ui2; Calculate the equivalent reactance from the generator to the second municipal district, denoted as XT2, The calculation formula of XT2 is: XT2 = (Un (o) - Ui2) / Ii2; Step S2422: Calculate the equivalent susceptance from the power generation end to the second municipal district, denoted as BT2; The calculation formula of BT2 is: BT2 = [3 × (Ii2) 2 ] 2 / 2 / Ui2; Calculate the equivalent resistance from the generator end to the second municipal district, recorded as RT2; The calculation formula of RT2 is: RT2 = [PG (o) × (Ui2) 2 ] / [(Pi2) 2 ]; Calculate the equivalent conductance from the power generation end to the second municipal district, denoted as GT2; The calculation formula of GT2 is: GT2 = (Pi2) / [(Ui2) 2 ]; Step S243: Similarly, the power generation end reaches the mnth municipal district; Step S2431: Calculate the receiving current of the mnth city district, denoted as Ii mn ; Ii mn The calculation formula is: Ii mn =Pi mn / Ui mn ; Calculate the equivalent reactance from the generator to the mnth municipal district, denoted as XT mn , XT mn The calculation formula is: XT mn = (Un (o) - Ui mn ) / Ii mn ; Step S2432: Calculate the equivalent susceptance from the power generation end to the mnth municipal district, denoted as BT mn ; BT mn The calculation formula is: BT mn =[3×(Ii mn ) 2 ] mn / 2 / Ui mn ; Calculate the equivalent resistance from the generator end to the mnth municipal district, denoted as RT mn ; RT mn The calculation formula is: RT mn =[PG(o)×(Ui mn ) 2 ] / [(Pi mn ) 2 ]; Calculate the equivalent conductance from the power generation end to the mnth municipal district, denoted as GT mn ; GT mn The calculation formula is: GT mn =(Pi mn ) / [(Ui mn ) 2 ].
6. The power big data analysis method based on cloud computing according to claim 1 is characterized in that: The specific steps of step S33 are as follows: Step S331: For the first power supply line, obtain the number of transformers other than the initial transformer and the receiving transformer, recorded as tr; obtain the length of the power transmission line from the initial transformer to the receiving transformer, recorded as lr; Step S332: initial transformer to the first transformer; Step S333: from the first transformer to the second transformer; Determine whether the second transformer is connected in parallel with other transformers; If not connected in parallel, proceed to step S334; If in parallel, then obtain the transformer information connected in parallel with the second transformer and proceed to step S335; Step S334: the second transformer is not connected in parallel with other transformers, and the adjusted transformation ratio of the second transformer is calculated; Step S335: the second transformer is connected in parallel with other transformers, and the adjusted transformation ratio of the second transformer is calculated; Step S3351: other transformers connected in parallel with the two transformers are used as parallel transformers; Count the number of transformers in parallel, denoted as oh; Get the current transformation ratio of the 1st to ohth parameter transformer and record it as: OB1, OB2~OB oh ; Among them, OB1 represents the current transformation ratio of the first parallel transformer; OB2 represents the current transformation ratio of the second parallel transformer; and so on. oh Indicates the current transformation ratio of the ohth parallel transformer; Step S3352: Using the transformation ratio DQ of the second transformer (2) As the first term, calculate DQ (2) :OB1:OB2:~:OB oh The ratio is recorded as: uu1:uu2:~:uu oh :uu (oh+1) ; Calculate uu1:uu2:~:uu oh :uu (oh+1) The inverse ratio is expressed as: nn1:nn2:~:nn oh :nn (oh+1) ; Step S3353: Repeat the same steps of calculating the adjusted transformation ratio of the second transformer on the premise that the second transformer is not connected in parallel with other transformers, and determine QQ (2) The value of QQ (2) Indicates the transformation ratio of the second transformer after adjustment; Calculate nn1 + nn2 + ~ + nn oh +nn (oh+1) The word "ann" is written as "ann"; The transformation ratio of the second transformer after secondary adjustment is denoted as cQQ (2) ;cQQ (2) The calculation formula is: cQQ (2) =(QQ (2) )×(nn1 / ann); Step S336: Repeat the same steps of step S332 to step S335 to adjust the transformation ratios of the 3rd to trth transformers in sequence.
7. The power big data analysis method based on cloud computing according to claim 6 is characterized in that: The specific steps of step S332 are as follows: Step S3321: Obtain the length of the transmission line from the initial transformer to the first transformer, denoted as lr (o-1) ; Calculate the distribution coefficient of the first transformer, denoted as β (1) β (1) The calculation formula is: (1) =lr (o-1) / lr; Calculate the voltage drop from the initial transformer to the first transformer, denoted as ΔU (o-1) ; ΔU (o-1) The calculation formula is: ΔU (o-1) =(β (1) )×(Un(o)-Ui1); Step S3322: Obtain the actual receiving voltage of the first transformer, recorded as DU (1) , actual power is recorded as DP (1) ; Get the current transformation ratio of the first transformer and record it as DQ (1) ; Calculate the actual receiving current of the first transformer, recorded as DI (1) DI (1) The calculation formula is: DI (1) =DP (1) / DU (1) ; The actual reactance of the first transformer is denoted by DX (1) , DX (1) The calculation formula is: DX (1) =[ΔU (o-1) ] / DI (1) ; The actual susceptance of the first transformer is denoted as DB (1) ; DB (1) The calculation formula is: DB (1) =[3×(DI (1) ) 2 ] 1 / 2 / DU (1) ; The actual resistance of the first transformer is denoted as DR (1) ; DR (1) The calculation formula is: DR (1) =[PG(o)×(DU (1) ) 2 ] / [mn×(DP (1) ) 2 ]; The actual conductance of the first transformer is denoted as DG (1) ; DG (1) The calculation formula is: DG (1) =(DP (1) ) / [(DU (1) ) 2 ]; Step S3323: Calculate DX (1) The offset coefficient with equivalent susceptance is denoted as zx (1) ; zx (1) The calculation formula is: zx (1) =1-[DX (1) / (β (1) ×XT1)]; Calculate DB (1) The deviation coefficient from the equivalent reactance is denoted as zb (1) ; zb (1) The calculation formula is: (1) =1-[DB (1) / (β (1) ×BT1)]; Calculating DR (1) The deviation coefficient from the equivalent reactance is denoted as zr (1) ; z (1) The calculation formula is: (1) =1-[DR (1) / (β (1) ×RT1)]; Calculate DG (1) The deviation coefficient from the equivalent reactance is denoted as zg (1) ; z (1) The calculation formula is: (1) =1-[β (1) / (DG (1) ×GT1)]; Step S3324: Calculate zx (1) +zb (1) +zr (1) +zg (1) The sum of is denoted as zz (1) ; Assume that the transformation ratio of the first transformer after adjustment is QQ (1) , judge zz (1) Size DeterminationQQ (1) Calculation formula of If zz (1) ≥0, then QQ (1) The calculation formula is: =[(2-zz (1) ) 1 / 2 ]×DQ (1) ; If zz (1) <0, then QQ (1) The calculation formula is: =[(1-zz (1) ) 1 / 2 ]×DQ (1) ; Among them, QQ (1) is a positive integer or a rational number.
8. The power big data analysis method based on cloud computing according to claim 6 is characterized in that: The specific steps of step S334 are as follows: Step S3341: Obtain the length of the transmission line from the first transformer to the second transformer, denoted as lr (1-2) ; Calculate the distribution coefficient of the second transformer, denoted as β (2) β (2) The calculation formula is: (2) =lr (1-2) / lr; Calculate the voltage drop from the first to the second transformer, denoted as ΔU (1-2) ; ΔU (1-2) The calculation formula is: ΔU (1-2) =(β (2) )×(Un(o)-Ui1); Step S3342: Obtain the actual receiving voltage of the second transformer, recorded as DU (2 , actual power is recorded as DP (2) ; Get the current transformation ratio of the second transformer and record it as DQ (2) ; Calculate the actual receiving current of the first transformer, recorded as DI (2) DI (2) The calculation formula is: DI (2) =DP (2) / DU (2) ; The actual reactance of the second transformer is denoted by DX (2) , DX (2) The calculation formula is: DX (2) =[ΔU (1-2) ] / DI (2) ; The actual susceptance of the second transformer is denoted as DB (2) ; DB (2) The calculation formula is: DB (2) =[3×(DI (2) ) 2 ] 1 / 2 / DU (2) ; The actual resistance of the second 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 second transformer is denoted as DG (1) ; DG (2) The calculation formula is: DG (2) =(DP (2) ) / [(DU (2) ) 2 ]; Step S3343: Calculate DX (2) The offset coefficient with equivalent susceptance is denoted as zx (2) ; zx (2) The calculation formula is: zx (2) =1-[DX (2) / (β (2) ×XT1)]; Calculate DB (1) The deviation coefficient from the equivalent reactance is denoted as zb (2) ; zb (2) The calculation formula is: (2) =1-[DB (2) / (β (2) ×BT1)]; Calculating DR (2) The deviation coefficient from the equivalent reactance is denoted as zr (2) ; z (2) The calculation formula is: (2) =1-[DR (2) / (β (2) ×RT1)]; Calculate DG (1) The deviation coefficient from the equivalent reactance is denoted as zg (2) ; z (2) The calculation formula is: (2) =1-[β (2) / (DG (2) ×GT1)]; Step S3344: Calculate zx (2) +zb (2) +zr (2) +zg (2) The sum of is denoted as zz (2) ; Assume that the transformation ratio of the second transformer after adjustment is QQ (2) , judge zz (2) Size DeterminationQQ (2) Calculation formula of If zz (2) ≥0, then QQ (2) The calculation formula is: =[(2-zz (2) ) 1 / 2 ]×DQ (2) ; If zz (2) <0, then QQ (1) The calculation formula is: =[(1-zz (2) ) 1 / 2 ]×DQ (2) ; Among them, QQ (2) is a positive integer or a rational number.
9. A power big data analysis system based on cloud computing, applicable to a power big data analysis method based on cloud computing as claimed in any one of claims 1 to 8, characterized in that: The system comprises: Power consumption analysis module: used to obtain the number of all urban districts in the target area; obtain the number of residents, current power consumption and historical power consumption of each urban district to obtain power consumption data; calculate the priority coefficient of each urban district based on the power consumption data to obtain the allocation sequence; Transmission analysis module: used to obtain the no-load electromotive force, rated voltage and reactance of the power generation end, obtain the receiving voltage and receiving power of each municipal district, and obtain the transmission data; construct the mathematical relationship between the power generation end, the power transmission end and the power receiving end according to the transmission data, and calculate the equivalent power parameter value of the power supply line in each municipal district according to the mathematical relationship; Power transmission control module: used to obtain the series-parallel relationship of the transformers between each power supply line and obtain the connection data of the transformers; based on the allocation sequence, according to the connection data and equivalent power parameter values of the transformers, the transformer ratio of each power supply line is adjusted to allocate power to each municipal area; Continuous update module: used to continuously obtain the current power consumption of each city district, update and adjust the power distribution of each city district.
Citation Information
Patent Citations
Wide-area distributed integrated topology analysis method for large power grid
CN104156835A
Urban intelligent regulation and control electrical power system
CN109119992A