GRID-FORMING ENERGY STORAGE SYSTEM AND CONTROL METHOD THERFOR

Abstract

A grid-forming energy storage system includes a data interface and transfer module, a state analysis and prediction module, an event response and strategy adjustment module, and an energy scheduling and management module. In the present disclosure, the energy management efficiency in the smart grid is significantly improved through real-time data formatting and efficient information transfer, real-time transmission and accuracy of data are ensured by using a JSON formatted data interface and a TCP/IP protocol, the data processing flow is optimized in conjunction with Apache Kafka, and the throughput capacity is improved and the response time is shortened, such that the grid manages resources more efficiently. Through state analysis and demand prediction by Spark Streaming, the grid responds in real time and predicts energy demand changes to reduce the energy waste.

Claims

1. A grid-forming energy storage system, comprising: a data interface and transfer module, configured to perform JSON formatting on key parameters of a voltage, a current, a temperature, and a remaining capacity of an energy storage unit, according to data of an energy storage unit, wherein the data is transmitted to a central control configuration through a TCP/IP protocol to obtain an encoded data stream, and the encoded data stream is inputted to an Apache Kafka message queue to acquire a queue data stream; a state analysis and prediction module, configured to input the queue data stream to a real-time data processing engine based on Spark Streaming to analyze data of the voltage, the current, the temperature, and the capacity, so as to obtain an analysis result of state time series, and dynamically estimate an energy demand of each unit to acquire a demand prediction result according to the analysis result of state time series; an event response and strategy adjustment module, configured to identify an energy demand peak and an emergent bulk power consumption event according to the demand prediction result, match an event pattern, execute event processing by using an Esper engine to obtain an event response strategy, and adjust energy supply by means of the event response decision to acquire an optimized energy distribution plan; and an energy scheduling and management module, configured to receive the optimized energy distribution plan, adjust a charge-discharge strategy according to demand, record an energy flow condition, and analyze a charge-discharge efficiency of energy storage and adaptability to obtain a charge-discharge scheduling data result of energy storage of the smart grid.

2. The grid-forming energy storage system according to claim 1, wherein the step of acquiring the queue data stream comprises: based on the data of the energy storage unit, converting each key parameter value of the energy storage unit into an attribute of a JSON object through a formatting function, generating a formatted data packet by using the following formula: $\mathrm{Ja}=\underset{kK}{.Math.}f\left(p_k\right)w_{k}N$ in the formula, K denotes the key parameter, p.sub.k denotes a real-time numerical value of a parameter k, f denotes a data formatting function, w.sub.k denotes a weight of the parameter k, N denotes the total number of parameters, and Ja denotes the formatted data packet; by means of a transmitting function, transmitting the formatted data packet through the TCP/IP protocol, generating an encoded data stream by using the following formula: $\mathrm{Ea}=\frac{1}{\sqrt{n}}\stackrel{n}{\underset{i=1}{.Math.}}s\left({\mathrm{Ja}}_i\right)$ in the formula, Ja.sub.i denotes the i.sup.th formatted data packet, n denotes the quantity of the data packets, s denotes the transmitting function, and Ea denotes the encoded data stream; and inputting the encoded data stream to the Apache Kafka message queue, and by means of a queue function, generating a queue data stream by using the following formula: $\mathrm{Qa}=q\left(\mathrm{Ea}\right)\left(1+\frac{1}{L}\right)$ in the formula, Ea denotes the encoded data stream, q denotes the queue function, L denotes a load factor, and Qa denotes the queue data stream.

3. The grid-forming energy storage system according to claim 2, wherein the step of acquiring the analysis result of state time series comprises: performing real-time processing on the queue data stream through a Spark Streaming engine, adjusting a processing speed and responsiveness of the data stream by applying a transfer function and combining an adjustment parameter, and generating a processed data stream by using the following formula: $\mathrm{St}=g\left(\mathrm{Qa}\right).Math.\left(+\frac{{}_0}{1+e^{-k.Math.\mathrm{len}\left(\mathrm{Qa}\right)}}\right)$ in the formula, is a baseline adjustment factor, .sub.0 and k are dynamic adjustment parameters, len(Qa) is configured to calculate a queue length, St denotes the processed data stream, Qa denotes an inputted queue data stream, and g is a real-time processed transfer function; extracting key electrical features from the processed data stream, comprising the voltage, current, temperature, and capacity, and performing weighting in conjunction with a feature weight, and generating a data stream feature set by using the following formula: $\mathrm{Ft}=\frac{1}{Z}\left(\underset{k\left\{V,I,T,C\right\}}{.Math.}h\left(k,\mathrm{St}\right).Math.{}_k\right)$ in the formula, h is a feature extraction function, .sub.k is the feature weight, Z is a normalization constant, Ft denotes the data stream feature set, and V, I, T, C are respectively the voltage, the current, the temperature, and the capacity; and performing time series analysis on the data stream feature set, by means of an analytic function in conjunction with a time factor, generating an analysis result of state time series by using the following formula: $\mathrm{At}=m\left(\mathrm{Ft}\right).Math.{\left(\frac{t}{+t}\right)}^{}$ in the formula, t denotes a time variable, is a decay constant, is a time-dependent intensity factor, At denotes the analysis result of state time series, Ft denotes the data stream feature set, and m denotes the time series analytic function.

4. The grid-forming energy storage system according to claim 3, wherein the step of acquiring the demand prediction result specifically comprises: integrating energy data of each unit in the analysis result of state time series, combining historical data with current data by means of an enhanced integration function, and generating integrated data by using the following formula: $\mathrm{Bo}=f\left(\mathrm{At}\right){\left(\stackrel{n}{\underset{i=1}{.Math.}}\frac{D_i}{\sqrt{n}}\right)}^2$ in the formula, D.sub.i denotes a historical energy demand of the i.sup.th unit, n denotes the total number of the units, At denotes the analysis result of state time series, f denotes the enhanced integration function, and Bo denotes the integrated data; applying a prediction model to the integrated data to dynamically calculate a demand in a future time period, and generating a predicted energy demand by using the following formula: $\mathrm{Po}=p\left(\mathrm{Bo}\right).Math.{\left(\frac{t}{+t}\right)}^{1/}$ in the formula, is a time delay constant, tis a time interval from a time when data is integrated to a predicted time, Po denotes the predicted energy demand, Bo denotes the integrated data, and p denotes the demand prediction model; and outputting the predicted energy demand in a formatted manner, matching the predicted energy demand to an energy unit management configuration, and generating a demand prediction result by using the following formula: $\mathrm{Do}=\mathrm{format}\left(\mathrm{Po}\right){}^2$ in the formula, format is a formatting function, is a formatting parameter, Do denotes a demand prediction result, and Po denotes the predicted energy demand.

5. The grid-forming energy storage system according to claim 4, wherein the step of acquiring the event response decision specifically comprises: defining the energy demand peak and identifying the emergent bulk power consumption event by using the demand prediction result, and by applying a threshold identification method, generating an event pattern by using the following formula: $\mathrm{Eq}=\left\{x\mathrm{Do}|x>.Math.\left(1+\frac{1}{\log \left(e+n\right)}\right)\right\}$ in the formula, is a dynamically adjusted demand peak threshold, n denotes the quantity of data points in a time series, Eq denotes the event pattern, Do denotes the demand prediction result, and x denotes energy demand data of a single data point; applying the event pattern to the Esper engine for pattern matching to identify a key energy event, and generating a matched event by using the following formula: $\mathrm{Mq}=\mathrm{match}\left(\mathrm{Eq}\right)\left(\frac{{.Math.}_{i=1}^n{}_i}{}\right)$ in the formula, .sub.i is an importance weight of each event, is the sum of the event weights, Mq denotes the matched event, Eq denotes the event pattern, and match denotes a pattern match function of the Esper engine; and based on the matched event, executing a decision logic and responding to the identified event, generating an event response decision by using the following formula: $\mathrm{Wq}=\mathrm{decision}\left(\mathrm{Mq}\right)\left(\frac{1}{1+e^{-k.Math.m}}\right)$ in the formula, k is the adjustment coefficient, m denotes the quantity of the matched events, Wq denotes the event response decision, Mq denotes the matched event, and decision denotes a decision function applied to the matched event.

6. The grid-forming energy storage system according to claim 5, wherein the step of acquiring the optimized energy distribution plan specifically comprises: extracting data from the event response decision, analyzing and identifying the energy supply strategy needed to be adjusted, and generating analyzed energy data by using the following formula: $\mathrm{Ar}=\left(\stackrel{n}{\underset{i=1}{.Math.}}{\mathrm{Wq}}_i.Math.{}_i\right)n^{1.5}$ in the formula, Wq.sub.i denotes data of the i.sup.th event response decision, .sub.i denotes the adjustment coefficient, and Ar denotes the analyzed energy data; adjusting the energy supply strategy according to the analyzed energy data, and matching the current energy demand with the predicted energy demand, and generating adjusted supply data by using the following formula: $\mathrm{Vr}=m\left(\mathrm{Ar}\right).Math.\sin \left(\frac{.Math.\mathrm{Ar}}{\max \left(\mathrm{Ar}\right)}\right)$ in the formula, m is a function that adjusts the energy supply, Vr denotes the adjusted supply data, max(Ar) denotes the maximum value in the analyzed data, and Ar denotes the analyzed energy data; and converting the adjusted supply data into an operation strategy, and forming the optimized energy distribution plan by using the following formula: $\mathrm{Hr}=\mathrm{strategize}\left(\mathrm{Vr}\right).Math.\left(1+\frac{\log \left(1+\mathrm{Vr}\right)}{\log \left(2\right)}\right)$ in the formula, strategize is a function that adjusts and converts the supply into a strategy, Hr denotes the optimized energy distribution plan, Vr denotes the adjusted supply data, log(1+Vr) denotes a logarithmic transformation of small value adjustment, and log(2) denotes a base of logarithm.

7. The grid-forming energy storage system according to claim 6, wherein the step of acquiring the charge-discharge scheduling data result of energy storage of the smart grid specifically comprises: receiving the optimized energy distribution plan, and forming a summarized data view in conjunction with real-time energy consumption data by using the following formula: $\mathrm{Iu}=\frac{{.Math.}_{i=1}^n\left({\mathrm{Hr}}_i{C}_i\right)}{\sqrt{n}}$ in the formula, Iu denotes the summarized data, Hr.sub.i denotes the i.sup.th strategy result data point, C.sub.i denotes corresponding energy consumption data, and n is the quantity of the data points; monitoring energy flow by means of the summarized data view, and matching the supply with the demand, adjusting energy charge and discharge, and generating an energy flow state by using the following formula: $\mathrm{Fu}={\left(\stackrel{n}{\underset{i=1}{.Math.}}{\mathrm{Iu}}_i{}_i\right)}^{\frac{1}{v}}$ in the formula, .sub.i is a weight factor of each data point, v is a global parameter that adjusts a nonlinear response, Fu denotes the energy flow state, and Iu denotes the summarized data; and based on the energy flow state, analyzing an efficiency and adaptability of energy distribution, and generating the charge-discharge scheduling data result of energy storage of the smart grid by using the following formula: $\mathrm{Eu}=\left(\frac{{.Math.}_{i=1}^n{\mathrm{Fu}}_i}{{}_i}\right).Math.$ in the formula, .sub.i is an efficiency benchmark of each state, is a coefficient for adaptive analysis, Eu denotes the charge-discharge scheduling data result of energy storage of the smart grid, and Fu denotes the energy flow state.

8. A grid-forming energy storage control method, executed through the grid-forming energy storage system according to claim 1, comprising the following steps: collecting the data of voltage, current, temperature, and remaining capacity of the energy storage unit through a control interface, formatting the data into JSON and controlling a data transmission process through the TCP/IP protocol, and optimizing a size and a transmission frequency of a data packet and matching a network condition to generate an encoded data stream; transmitting the data to Apache Kafka for control processing by means of the encoded data stream, and optimizing reliability of message processing by controlling sequencing and error detection of the data stream in the message queue to obtain a queue data stream; performing data control analysis on the queue data stream in a Spark Streaming environment, continuously monitoring the voltage, current, temperature, and capacity, analyzing a time series variation of parameters, and controlling and identifying an abnormal pattern to generate the analysis result of state time series; based on the analysis result of state time series, controlling a calculation of demand prediction, defining and identifying the energy peak value and the emergent bulk power consumption event, controlling real-time response processing of the event, and adjusting an output configuration of the energy storage unit to acquire an event response decision; and controlling energy supply and output of energy storage by using the event response decision, adjusting the charge-discharge strategy according to a demand, controlling an energy flow and consumption condition of a network, and analyzing and adjusting the energy distribution strategy and optimizing energy storage control configuration performance to establish the charge-discharge scheduling data result of energy storage of the smart grid.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0062] FIG. 1 is a flow chart of a system of the present disclosure.

[0063] FIG. 2 is a flow chart of a queue data stream in the present disclosure.

[0064] FIG. 3 is a flow chart of an analysis result of state time series in the present disclosure.

[0065] FIG. 4 is a flow chart of a demand prediction result in the present disclosure.

[0066] FIG. 5 is a flow chart of an event response decision in the present disclosure.

[0067] FIG. 6 is a flow chart of an optimized energy distribution plan in the present disclosure.

[0068] FIG. 7 is a flow chart of a charge-discharge scheduling data result of energy storage of a smart grid in the present disclosure.

DETAILED DESCRIPTIONS OF THE EMBODIMENTS

[0069] In order to make objectives, technical solutions, and advantages of the present disclosure clearer, the present disclosure will be further described in detail below in conjunction with drawings and embodiments. It should be understood that the specific embodiments described herein are merely used to explain the present disclosure, rather than limiting the present disclosure.

[0070] In the description of the present disclosure, it should be noted that an orientation relationship or a position relationship indicated by the terms length, width, upper, lower, front, back, left, right, vertical, horizontal, top, bottom, inside, outside, and the like is an orientation relationship or a position relationship based on the drawing and is merely for ease of describing present disclosure and simplifying the description, rather than indicating or implying that a specified apparatus or element necessarily has a specific orientation or is constructed and operated in a specific orientation. Therefore, the terms should not be construed as a limitation on the present disclosure. In the description of the present disclosure, unless otherwise specified, a plurality of means two or more than two.

Embodiment I

[0071] Referring to FIG. 1, the present disclosure provides a technical solution: a grid-forming energy storage system includes: [0072] a data interface and transfer module, configured to perform JSON formatting on key parameters of a voltage, a current, a temperature, and a remaining capacity of an energy storage unit according to data of the energy storage unit, where the data is transmitted to a central control configuration through a TCP/IP protocol to obtain an encoded data stream, and the encoded data stream is inputted to an Apache Kafka message queue to acquire a queue data stream; [0073] a state analysis and prediction module, configured to input the queue data stream to a real-time data processing engine based on Spark Streaming to analyze data of the voltage, the current, the temperature, and the capacity, so as to obtain an analysis result of state time series, and dynamically estimate an energy demand of each unit according to the analysis result of state time series to acquire a demand prediction result; [0074] an event response and strategy adjustment module, configured to identify an energy demand peak and an emergent bulk power consumption event according on the demand prediction result, match an event pattern, execute event processing by using an Esper engine to obtain an event response strategy, and adjust energy supply by means of the event response decision to acquire an optimized energy distribution plan; and [0075] an energy scheduling and management module, configured to receive the optimized energy distribution plan, adjust a charge-discharge strategy according to demand, record an energy flow condition, and analyze a charge-discharge efficiency of energy storage and adaptability to obtain a charge-discharge scheduling data result of energy storage of the smart grid.

[0076] The encoded data stream includes data serialization, transmission optimization, and protocol compatibility. The queue data stream includes a message queue length, a cache strategy, and transmission stability. The analysis result of state time series includes voltage trend analysis, current pattern identification, and temperature anomaly detection. The demand prediction result includes consumption estimation, demand responsiveness, and prediction accuracy. The event response decision includes strategy formation speed, event trigger accuracy, and decision execution efficiency. The optimized energy distribution plan includes supply adjustment range, output balance, and system response speed. The charge-discharge scheduling data result of energy storage of the smart grid includes monitoring system efficiency, data report precision, and scheduling adaptability.

[0077] Referring to FIG. 2, the step of acquiring the queue data stream specifically includes: [0078] converting each key parameter value of the energy storage unit into an attribute of a JSON object through a formatting function according to the data of the energy storage unit, generating a formatted data packet by using the following formula:

[00019]$\mathrm{Ja}=\underset{kK}{.Math.}f\left(p_k\right)w_{k}N$ [0079] in the formula, K denotes the key parameter, p.sub.k denotes a real-time numerical value of a parameter k, f denotes a data formatting function, w.sub.k denotes a weight of the parameter k, N denotes the total number of parameters, and Ja denotes the formatted data packet; [0080] by means of a transmitting function, transmitting the formatted data packet through the PCP/IP protocol, generating an encoded data stream by using the following formula:

[00020]$\mathrm{Ea}=\frac{1}{\sqrt{n}}\underset{i=1}{\overset{n}{.Math.}}s\left({\mathrm{Ja}}_i\right)$ [0081] in the formula, Ja.sub.i denotes an i.sup.th formatted data packet, n denotes the quantity of the data packets, s denotes the transmitting function, and Ea denotes the encoded data stream; and [0082] inputting the encoded data stream to the Apache Kafka message queue, and by means of a queue function, generating a queue data stream by using the following formula:

[00021]$\mathrm{Qa}=q\left(\mathrm{Ea}\right)\left(1+\frac{1}{L}\right)$ [0083] in the formula, Ea denotes the encoded data stream, q denotes the queue function, L denotes a load factor, and Qa denotes the queue data stream. A formatted data packet formula:

[00022]$\mathrm{Ja}=\underset{kK}{.Math.}f\left(p_k\right)w_{k}N$ [0084] K: a set of all key parameters, for example, voltage (V), current (I), temperature (T), and remaining capacity (C)); [0085] p.sub.k: a real-time numerical value of each parameter, for example, V=3.7V, I=2.0 A, T=35 C., C=75%; [0086] w.sub.k: a weight of each parameter, assuming that the weights are equal here, and each of the weights is 1; [0087] N: the total number of the parameters, which is 4 here; [0088] data preparation: V=3.7, 1=2.0, T=35, C=0.75, w.sub.V=w.sub.I=w.sub.T=w.sub.C=1, N=4; [0089] Ja is calculated as follows:

[00023]$\mathrm{Ja}=\frac{\left(3.71\right)+\left(2.01\right)+\left(351\right)+\left(0.751\right)}{4}=\frac{41.45}{4}=10\text{.3625}.$

[0090] This result (10.3625) denotes a parameter value after weighted averaging for further data processing.

[0091] An encoded data stream formula:

[00024]$\mathrm{Ea}=\frac{1}{\sqrt{n}}{.Math.}_{i=1}^{n}s\left(J{a}_i\right).$ [0092] Ja.sub.i: a formatted data packet formula, which is 10.3625 in the above example; [0093] n: the quantity of the data packets, assuming that 3 data packets exist; [0094] It is assumed: Ja.sub.1=Ja.sub.2=Ja.sub.3=10.3625, n=3; [0095] Ea is calculated as follows:

[00025]$\mathrm{Ea}=\frac{1}{\sqrt{3}}{.Math.}_{i=1}^{3}10.3625=\frac{1}{1.732}31.0875=17.943.$

[0096] Result explanation: this result (17.943) denotes a processed weighted data stream adapted to network transmission.

[0097] A queue data stream formula:

[00026]$\mathrm{Qa}=q\left(\mathrm{Ea}\right)\left(1+\frac{1}{L}\right).$ [0098] Ea: an encoded data stream, which is 17.943 in the above example; [0099] L: a system load factor, assuming that it is 5; [0100] it is assumed: Ea=17.943, L=5; [0101] Qa is calculated as follows:

[00027]$\mathrm{Qa}=17.943\left(1+\frac{1}{5}\right)=17.9431.2=21.532.$

[0102] This result (21.532) denotes a processed data stream inputted to Apache Kafka for subsequent message queue processing.

[0103] Referring to FIG. 3, the step of acquiring the analysis result of state time series specifically includes: [0104] performing real-time processing on the queue data stream through a Spark Streaming engine, adjusting a processing speed and responsiveness of the data stream by applying a transfer function and combining an adjustment parameter, and [0105] generating a processed data stream by using the following formula:

[00028]$St=g\left(\mathrm{Qa}\right).Math.\left(+\frac{{}_0}{1+e^{-k.Math.\mathrm{len}\left(\mathrm{Qa}\right)}}\right)$ [0106] in the formula, is a baseline adjustment factor, .sub.0 and k are dynamic adjustment parameters, len(Qa) is configured to calculate a queue length, St denotes the processed data stream, Qa denotes an inputted queue data stream, and g is a real-time processed transfer function; [0107] extracting key electrical features from the processed data stream, including the voltage, the current, the temperature, and the capacity, and performing weighting in conjunction with a feature weight, and [0108] generating a data stream feature set by using the following formula:

[00029]$\mathrm{Ft}=\frac{1}{Z}\left(\underset{k\left\{V,I,T,C\right\}}{.Math.}h\left(k,\mathrm{St}\right).Math.{}_k\right)$ [0109] in the formula, h is a feature extraction function, .sub.k is the feature weight, Z is a normalization constant, Ft denotes the data stream feature set, and V, I, T, C are respectively the voltage, the current, the temperature, and the capacity; and [0110] performing time sequence analysis on the data stream feature set, by means of an analytic function in conjunction with a time factor, generating an analysis result of state time series by using the following formula:

[00030]$\mathrm{At}=m\left(\mathrm{Ft}\right).Math.{\left(\frac{t}{+t}\right)}^{}$ [0111] in the formula, t denotes a time variable, is a decay constant, is a time-dependent intensity factor, At denotes the analysis result of state time series, Ft denotes the data stream feature set, and m denotes the time sequence analytic function.

[0112] A processed data stream formula:

[00031]$\mathrm{St}=g\left(\mathrm{Qa}\right).Math.\left(+\frac{{}_0}{1+e^{-k.Math.\mathrm{len}\left(Qa\right)}}\right).$ [0113] g(Qa) denotes a real-time processing transfer function applied to the inputted queue data stream Qa; [0114] : a baseline adjustment factor, configured to control a basic rate of processing; [0115] .sub.0: a dynamic adjustment parameter, configured to adjust a change of response with the data stream length; [0116] k: a coefficient to adjust sensitivity, which affects the application intensity of .sub.0; [0117] len(Qa): a length of Qa, i.e., the quantity of data items in the data stream; [0118] it is assumed that len(Qa)=100 items of data, g(Qa)=200 (denotes the processed data base value), =1.5, .sub.0=0.5, k=0.1; [0119] e.sup.k.Math.len (Qa)=e.sup.0.1.Math.100=e.sup.100.000045 is calculated; [0120] an adjustment item

[00032]$\frac{0.5}{1+0.000045}0.4999775$ is calculated; [0121] a total adjustment factor

[00033]$+\frac{{}_0}{1+e^{-k.Math.\mathrm{len}\left(\mathrm{Qa}\right)}}=1.5+0.49997752.0$ is calculated [0122] a final result St=200.Math.2.0=400 is calculated; [0123] the finally processed data stream St=400 indicates that the adjusted processing intensity is two times that of a certain value, indicating that the processing engine has strong responsiveness to a large data flow.

[0124] A data stream feature set formula:

[00034]$\mathrm{Ft}=\frac{1}{Z}\left(\underset{k\left\{V,I,T,C\right\}}{.Math.}h\left(k,\mathrm{St}\right).Math.{}_k\right)$ [0125] h(k, St): a feature extraction function, for data items k (voltage V, current I, temperature T, and capacity C) in St; [0126] .sub.k: a weight of each feature; [0127] Z: a normalization constant, ensuring that the weights are reasonably distributed; [0128] it is assumed: h(V, St)=50, h(I, St)=60, h(T, St)=40, h(C, St)=70, .sub.V=0.2, .sub.I=0.3, .sub.T=0.1, .sub.C=0.4, Z=2; [0129] the sum of weighted features is calculated as follows:

[00035]$\underset{k\left\{V,I,T,C\right\}}{.Math.}h\left(k,\mathrm{St}\right).Math.{}_k=500.2+600.3+400.1+700.4=10+18+4+28=60$ [0130] Ft is calculated as follows:

[00036]$\mathrm{Ft}=\frac{1}{2}60=30.$

[0131] The final feature set Ft=30 denotes a weighted and normalized result after comprehensively considering all features, for subsequent time series analysis.

[0132] A formula of the analysis results of state time series:

[00037]$\mathrm{At}=m\left(\mathrm{Ft}\right).Math.{\left(\frac{t}{+t}\right)}^{}$ [0133] m(Ft): a time sequence analytic function for processing Ft; [0134] t: a current time point; [0135] : a time decay constant; [0136] : a time-dependent intensity factor; [0137] it is assumed that m(Ft)=3Ft=90 (it is assumed that the base value of the time sequence analysis is three times of Ft), t=5, =10, 8=2; [0138] a time factor is calculated as follows:

[00038]${\left(\frac{t}{+t}\right)}^{}={\left(\frac{5}{10+5}\right)}^2={\left(\frac{5}{15}\right)}^2={\left(\frac{1}{3}\right)}^2=\frac{1}{9};$ [0139] a final result At is calculated as follows:

[00039]$\mathrm{At}=90\frac{1}{9}=10.$

[0140] The final analysis result of state time series At=10 denotes that under the influence of a given time point, an analysis result indicates weakening of feature influence, reflecting a significant influence of the time factor on the result.

[0141] Referring to FIG. 4, the step of acquiring the prediction result specifically includes: [0142] integrating energy data of each unit in the analysis result of state time series, combining historical data with current data by means of an enhanced integration function, and [0143] generating integrated data by using the following formula:

[00040]$\mathrm{Bo}=f\left(\mathrm{At}\right){\left(\underset{i=1}{\overset{n}{.Math.}}\frac{D_i}{\sqrt{n}}\right)}^2$ [0144] in the formula, D.sub.i denotes a historical energy demand of the i.sup.th unit, n denotes the total number of the units, At denotes the analysis result of state time series, f denotes the enhanced integration function, and Bo denotes the integrated data; [0145] applying a prediction model to the integrated data to dynamically calculate a demand in a future time period, and generating a predicted energy demand by using the following formula:

[00041]$\mathrm{Po}=p\left(\mathrm{Bo}\right).Math.{\left(\frac{t}{+t}\right)}^{1/}$ [0146] in the formula, is a time delay constant, t is a time interval from a time when data is integrated to a predicted time, Po denotes the predicted energy demand, Bo denotes the integrated data, and p denotes the demand prediction model; and [0147] outputting the predicted energy demand in a formatted manner, matching the predicted energy demand to an energy unit management configuration, and [0148] generating a demand prediction result by using the following formula Do=format(Po).sup.2, in the formula, format is a formatting function, is a formatting parameter, Do denotes a demand prediction result, and Po denotes the predicted energy demand.

[0149] An integrated data stream formula:

[00042]$\mathrm{Bo}=f\left(\mathrm{At}\right){\left(\underset{i=1}{\overset{n}{.Math.}}\frac{D_i}{\sqrt{n}}\right)}^2$ [0150] At: an analysis result of state time series, where it is assumed that the analysis result of state time series is 1000, representing an energy state indicator at a certain time point; [0151] f(At): an integration function, which may be assumed as f(At)=At0.1, indicating scaling of the state time sequence result to facilitate combination with historical data; [0152] D.sub.i: a historical energy demand of the i.sup.th unit, where it is assumed that the energy demand of three units is D.sub.1=500, D.sub.2=300, D.sub.3=400; [0153] n: the total number of the units, which is 3 here; [0154] f(At): is calculated as follows: [0155] f(At)=10000.1=100, [0156] the square of the average energy demand is calculated as follows:

[00043]$\underset{i=1}{\overset{n}{.Math.}}\frac{D_i}{\sqrt{n}}=\left(\frac{500+300+400}{\sqrt{3}}\right)=\frac{1200}{\sqrt{3}}692\text{.82}$${\left({.Math.}_{i=1}^n\frac{D_i}{\sqrt{n}}\right)}^2=692.82^{2}480000.$

[0157] The finally integrated data Bo:

[00044]$\mathrm{Bo}=100480000=48000000.$

[0158] This integrated data denotes a comprehensive energy demand prediction indicator after considering the historical energy demand and the current state.

[0159] An energy prediction demand formula:

[00045]$\mathrm{Po}=p\left(\mathrm{Bo}\right).Math.{\left(\frac{t}{+t}\right)}^{1/}$ [0160] p(Bo): a demand prediction model function, which is assumed as p(Bo)=Bo0.05, indicating that the integrated data is converted into a specific energy demand prediction; [0161] : a time decay constant, which is assumed as 5; [0162] t: a time interval, from a time when the data is integrated to a predicted time, which is assumed as 10 h; [0163] p(Bo) is calculated as follows:

[00046]$p\left(\mathrm{Bo}\right)=480000000.05=2400000,$${\left(\frac{t}{+t}\right)}^{1/}={\left(\frac{10}{5+10}\right)}^{1/5}0.896.$

[0164] A final prediction demand Po:

[00047]$\mathrm{Po}=24000000.8962150400.0$

[0165] This prediction demand denotes a prediction energy demand taking the time decay into consideration.

[0166] A demand prediction result formula:

[00048]$\mathrm{Do}=\mathrm{format}\left(\mathrm{Po}\right){}^2$ [0167] format(Po): a formatting function, which is assumed to be simplified to format(Po)=Po, without format change; [0168] : a formatting parameter, which affects the output format and is assumed as 10; [0169] a formatted prediction demand:

[00049]$\mathrm{format}\left(\mathrm{Po}\right)=2150400.0$

[0170] A final output is calculated as follows:

[00050]$\mathrm{Do}=2150400{10}^2=2.1504{10}^7$

[0171] This result denotes the formatted energy demand prediction, ready to be used to report to the management system.

[0172] Referring to FIG. 5, the step of acquiring the event response decision specifically includes: [0173] defining the energy demand peak and identifying the emergent bulk power consumption event by using the demand prediction result, and by applying a threshold identification method, generating an event pattern by using the following formula:

[00051]$\mathrm{Eq}=\left\{x\mathrm{Do}.Math.x>.Math.\left(1+\frac{1}{\log \left(e+n\right)}\right)\right\}$ [0174] in the formula, is a dynamically adjusted demand peak threshold, n denotes the quantity of data points in a time sequence, Eq denotes the event pattern, Do denotes the demand prediction result, and x denotes energy demand data of a single data point; [0175] applying the event pattern to the Esper engine for pattern matching to identify a key energy event, and generating a matched event by using the following formula:

[00052]$\mathrm{Mq}=\mathrm{match}\left(\mathrm{Eq}\right)\left(\frac{{.Math.}_{i=1}^n{}_i}{}\right)$ [0176] in the formula, .sub.i is an importance weight of each event, is the sum of the event weights, Mq denotes the matched event, Eq denotes the event pattern, and match denotes a pattern match function of the Esper engine; and [0177] executing a decision logic and responding to the identified event according to the matched event to generate an event response decision by using the following formula:

[00053]$\mathrm{Wq}=\mathrm{decision}\left(\mathrm{Mq}\right)\left(\frac{1}{1+e^{-k.Math.m}}\right)$ [0178] in the formula, k is the adjustment coefficient, m denotes the quantity of the matched events, Wq denotes the event response decision, Mq denotes the matched event, and decision denotes a decision function applied to the matched event.

[0179] An event pattern formula:

[00054]$\mathrm{Eq}=\left\{x\mathrm{Do}.Math.x>.Math.\left(1+\frac{1}{\log \left(e+n\right)}\right)\right\}$ [0180] Do: energy demand prediction data in every 10 minute in a hour, with a unit of megawatt (MW), Do=[10, 15, 20, 25, 30, 35]; [0181] : a set demand peak threshold is 20 MW; [0182] n: the quantity of the data points is 6 (i.e., six ten-minute intervals in 1 h); [0183] a dynamic threshold adjustment factor is calculated as follows:

[00055]$1+\frac{1}{\log \left(e+6\right)}1+\frac{1}{\log \left(8.71828\right)}1+0.13031.1303;$ [0184] then a new threshold is calculated as follows: =201.130322.606; [0185] energy demands meeting conditions are screened according to the new threshold:

[00056]$\mathrm{Eq}=\left\{x\left[10,15,20,25,30,35\right]|x>22.606\right\}=\left\{25,30,35\right\};$ [0186] a result Eq denotes that during observation, the energy demands at three time points exceed the adjusted peak threshold, which affects the stability and operation of the energy system.

[0187] An event pattern formula:

[00057]$\mathrm{Mq}=\mathrm{match}\left(\mathrm{Eq}\right)\left(\frac{{.Math.}_{i=1}^n{}_i}{}\right).$

[0188] It is assumed that match (Eq) simply returns the quantity of matched events, which is 3 here (three events with the energy demands exceeding the threshold); [0189] .sub.i: a weight of each event is assumed as 1; [0190] =n=6 (the sum of weights of all events); [0191] a weight ratio is calculated as follows:

[00058]$\frac{{.Math.}_{i=1}^3{}_i}{}=\frac{3}{6}=0.5.$

[0192] A final result is calculated as follows:

[00059]$\mathrm{Mq}=30.5=1.5.$

[0193] A result Mq=1.5 denotes that after weight adjustment, the system identifies a key event with 1.5 average influence units, which is beneficial for a system decision maker to further analyze and respond.

[0194] An event response decision formula:

[00060]$\mathrm{Wq}=\mathrm{decision}\left(\mathrm{Mq}\right)\left(\frac{1}{1+e^{-k.Math.m}}\right).$ [0195] It is assumed that decision (Mq) returns a basic weight of an event decision, and the basic weight is 2 here (it is assumed that two strategies need to be copied with); [0196] k=0.5, m=1.5 (a matched event obtained from the previous step); [0197] An exponential response is calculated as follows:

[00061]$\frac{1}{1+e^{-k.Math.m}}=\frac{1}{1+e^{-0.51.5}}=\frac{1}{1+e^{-0.75}}\frac{1}{1+0.4724}0.679;$ [0198] a final result is calculated as follows:

[00062]$\mathrm{Wq}=20.6791.358.$

[0199] A result Wq1.358 denotes that a strategy response intensity with about 1.358 units is needed after adjustment to process the identified event according to a decision weight of the system in the current state.

[0200] Referring to FIG. 6, the step of acquiring the optimized energy distribution plan specifically includes: [0201] extracting data from the event response decision, analyzing and identifying the energy supply strategy needed to be adjusted, and [0202] generating analyzed energy data by using the following formula:

[00063]$\mathrm{Ar}=\left(\underset{i=1}{\overset{n}{.Math.}}{\mathrm{Wq}}_i.Math.{}_i\right)n^{1.5}$ [0203] in the formula, Wq.sub.i denotes data of the i.sup.th event response decision, .sub.i denotes the adjustment coefficient, and Ar denotes the analyzed energy data; [0204] adjusting the energy supply strategy according to the analyzed energy data, and matching the current energy demand with the predicted energy demand, and [0205] generating adjusted supply data by using the following formula:

[00064]$\mathrm{Vr}=m\left(\mathrm{Ar}\right).Math.\sin \left(\frac{.Math.\mathrm{Ar}}{\max \left(Ar\right)}\right)$ [0206] in the formula, m is a function that adjusts the energy supply, Vr denotes the adjusted supply data, max(Ar) denotes the maximum value in the analyzed data, and Ar denotes the analyzed energy data; and converting the adjusted supply data into an operation strategy, and [0207] forming the optimized energy distribution plan by using the following formula:

[00065]$\mathrm{Hr}=\mathrm{strategize}\left(\mathrm{Vr}\right).Math.\left(1+\frac{\log \left(1+\mathrm{Vr}\right)}{\log \left(2\right)}\right)$ [0208] in the formula, strategize is a function that adjusts and converts the supply into a strategy, Hr denotes the optimized energy distribution plan, Vr denotes the adjusted supply data, log(1+Vr) denotes a logarithmic transformation of small value adjustment, and log(2) denotes a base of logarithm.

[0209] An analyzed energy data formula:

[00066]$\mathrm{Ar}=\left(\underset{i=1}{\overset{n}{.Math.}}{\mathrm{Wq}}_i.Math.{}_i\right)n^{1.5}$ [0210] Wq.sub.i: a value of the i.sup.th decision data point, which is assumed as an emergency level score of response; [0211] .sub.i: an adjustment coefficient of the i.sup.th data point, reflecting the relative importance of the data point; [0212] n: the total number of the data points, indicating the number of decision points taken into consideration; [0213] it is assumed that 3 decision data points exist, i.e., n=3, the decision data points Wq are respectively 10, 15, and 20, and the adjustment coefficients are respectively 1.2, 0.8, and 1.0; [0214] a weighted value of each decision point is calculated as follows:

[00067]${\mathrm{Wq}}_1.Math.{}_1=101.2=12,$${\mathrm{Wq}}_2.Math.{}_2=150.8=12,$${\mathrm{Wq}}_3.Math.{}_4=201.0=20;$ [0215] all weighted values are summated:

[00068]$\underset{i=1}{\overset{n}{.Math.}}{\mathrm{Wq}}_i.Math.{}_i=12+12+20=44$ [0216] divided by n.sup.1.5:

[00069]$n^{1.5}=3^{1.5}=5.196.$

[0217] Ar is calculated as follows:

[00070]$\mathrm{Ar}=\frac{44}{5.196}8.47.$

[0218] The obtained Ar8.47 denotes the weighted average decision emergency level.

[0219] An adjusted supply data formula:

[00071]$\mathrm{Vr}=m\left(\mathrm{Ar}\right).Math.\sin \left(\frac{.Math.\mathrm{Ar}}{\max \left(Ar\right)}\right)$ [0220] m(Ar): an energy adjustment function, which is assumed as m(x)=2x; [0221] Ar: analyzed energy data: [0222] max(Ar): the maximum value of Ar, which is assumed as 10; [0223] Ar8.47 (calculated from the previous step) is used; [0224] m(Ar) is calculated as follows:

[00072]$m\left(\mathrm{Ar}\right)=m\left(8.47\right)=28.47=16.94.$

[0225] A sinusoidal part is calculated as follows:

[00073]$\sin \left(\frac{.Math.\mathrm{Ar}}{\max \left(\mathrm{Ar}\right)}\right)=\sin \left(\frac{.Math.8.47}{10}\right)=\sin \left(2.654\right)0.445.$

[0226] A final result:

[00074]$\mathrm{Vr}=16.940.4457.54.$

[0227] The obtained Vr7.54 denotes the adjusted energy supply index.

[0228] An optimized energy distribution plan formula:

[00075]$\mathrm{Hr}=\mathrm{strategize}\left(\mathrm{Vr}\right)\left(1+\frac{\log \left(1+\mathrm{Vr}\right)}{\log \left(2\right)}\right)$ [0229] strategize(Vr): a function that converts energy supply into a specific strategy, which is assumed as strategize(Vr)=x.sup.2; [0230] Vr: the adjusted supply data; [0231] log(2): a base of Natural Logarithm 2, which is approximately equal to 0.693; [0232] Vr7.54 (calculated from the previous step) is used: [0233] a strategy conversion function is calculated as follows:

[00076]$\mathrm{strategize}\left(\mathrm{Vr}\right)=\mathrm{strategize}\left(7.54\right)={\left(7.54\right)}^2=56.8516;$ [0234] a logarithmic part is calculated as follows:

[00077]$\frac{\log \left(1+\mathrm{Vr}\right)}{\log \left(2\right)}=\frac{\log \left(1+7.54\right)}{\log \left(2\right)}=\frac{\log \left(8.54\right)}{0.693}3.41.$

[0235] A final result:

[00078]$\mathrm{Hr}=56.8516\left(1+3.41\right)294.44$

[0236] The obtained Hr249.44 denotes a comprehensive effect index of the finally formed adjustment strategy.

[0237] Referring to FIG. 7, the step of acquiring the charge-discharge scheduling data result of energy storage of the smart grid specifically includes: [0238] receiving the optimized energy distribution plan, and forming a summarized data view in conjunction with real-time energy consumption data by using the following formula:

[00079]$\mathrm{Iu}=\frac{{.Math.}_{i=1}^n\left({\mathrm{Hr}}_i{C}_i\right)}{\sqrt{n}}$ [0239] in the formula, Iu denotes the summarized data, Hr.sub.i denotes the i.sup.th strategy result data point, C.sub.i denotes corresponding energy consumption data, and n is the quantity of the data points; [0240] monitoring energy flow by means of the summarized data view, and matching the supply with the demand, adjusting energy charge and discharge, and [0241] generating an energy flow state by using the following formula:

[00080]$\mathrm{Fu}={\left(\stackrel{n}{\underset{i=1}{.Math.}}{\mathrm{Iu}}_i{}_i\right)}^{\frac{1}{v}}$ [0242] in the formula, .sub.i is a weight factor of each data point, v is a global parameter that adjusts a nonlinear response, Fu denotes the energy flow state, and Iu denotes the summarized data; and [0243] analyzing an efficiency and adaptability of energy distribution according to the energy flow state, and generating the charge-discharge scheduling data result of energy storage of the smart grid by using the following formula:

[00081]$\mathrm{Eu}=\left(\frac{{.Math.}_{i=1}^n{\mathrm{Fu}}_i}{{}_i}\right).Math.$ [0244] in the formula, .sub.i is an efficiency benchmark of each state, is a coefficient for adaptive analysis, Eu denotes the charge-discharge scheduling data result of energy storage of the smart grid, and Fu denotes the energy flow state.

[0245] A summarized data formula:

[00082]$\mathrm{Iu}=\frac{{.Math.}_{i=1}^n\left({\mathrm{Hr}}_i{C}_i\right)}{\sqrt{n}}$ [0246] Hr.sub.i denotes a data point of the i.sup.th optimized energy distribution plan, assuming that the data is from calculation of the previous step or a real-time system feedback; [0247] C.sub.i denotes the i.sup.th energy consumption data point according to the feedback of monitoring the system in real time; [0248] n denotes the total number of the data points for standardized data processing; [0249] it is assumed that 3 strategy result data points exist: Hr=[2, 3, 5], the corresponding energy consumption data is C=[100, 150, 200], so n=3; [0250] Iu is calculated as follows:

[00083]$\mathrm{Iu}=\frac{2100+3150+5200}{\sqrt{3}}=\frac{200+450+1000}{1.732}\frac{1650}{1.732}952.81.$

[0251] This result (952.81) denotes the integrated data index under a circumstance that the strategy adjustment and the energy consumption data are given.

[0252] An energy flow state formula:

[00084]$\mathrm{Fu}={\left(\stackrel{n}{\underset{i=1}{.Math.}}{\mathrm{Iu}}_i{}_i\right)}^{\frac{1}{v}}$ [0253] Iu: an integrated data point obtained from the first step; [0254] .sub.i: a weight factor of each data point, which is assumed as a standard weight acquired from the historical data or an expert system here; [0255] v: a global parameter, configured to adjust nonlinear response; [0256] Iu=952.81 in the previous example is continuously used, assuming .sub.i=[1, 1.5, 2], v=2; [0257] Fu is calculated as follows:

[00085]$\mathrm{Fu}={\left(952.811+952.811.5+952.812\right)}^{\frac{1}{2}}={\left(952.814.5\right)}^{\frac{1}{2}}65.48.$

[0258] This result (65.48) denotes an adjustment state index of an energy flow.

[0259] A charge-discharge scheduling data result formula of energy storage of a smart grid:

[00086]$\mathrm{Eu}=\left(\frac{{.Math.}_{i=1}^n{\mathrm{Fu}}_i}{{}_i}\right).Math.$ [0260] Fu: an energy flow state data point obtained from the previous step; [0261] .sub.i: an efficiency benchmark of each state, assuming that it is estimated according to an industrial standard or previous performance; [0262] : a global adaptability coefficient, which enhances the depth of the system adaptability analysis;

[00087]$\mathrm{Fu}=65.48\mathrm{is}\mathrm{used},\mathrm{assuming}{}_i=\left[1,1,1\right],=1.1;$ [0263] Eu is calculated as follows:

[00088]$\mathrm{Eu}=\left(\frac{65.48+65.48+65.48}{1+1+1}\right).Math.1.1=\left(65.48\right).Math.1.1=72.028.$

[0264] This result (72.028) denotes a comprehensive estimation index of the energy distribution efficiency and the system adaptability.

[0265] A grid-forming energy storage control method executed through the grid-forming energy storage system includes the following steps: [0266] S1: collecting the data of voltage, current, temperature, and remaining capacity of the energy storage unit through a control interface, formatting the data into JSON and controlling a data transmission process through the TCP/IP protocol, and optimizing a size and a transmission frequency of a data packet and matching a network condition to generate an encoded data stream; [0267] S2: transmitting the data to Apache Kafka for control processing by means of the encoded data stream, and optimizing reliability of message processing by controlling sequencing and error detection of the data stream in the message queue to obtain a queue data stream; [0268] S3: performing data control analysis on the queue data stream in a Spark Streaming environment, continuously monitoring the voltage, the current, the temperature, and the capacity, analyzing a time series variation of parameters, and controlling and identifying an abnormal pattern to generate the analysis result of state time series; [0269] S4: controlling a calculation of demand prediction according to the analysis result of state time series, defining and identifying the energy peak value and the emergent bulk power consumption event, controlling real-time response processing of the event, and adjusting an output configuration of the energy storage unit to acquire an event response decision; and [0270] S5: controlling energy supply and output of energy storage by using the event response decision, adjusting the charge-discharge strategy according to a demand, controlling an energy flow and consumption condition of a network, and analyzing and adjusting the energy distribution strategy and optimizing energy storage control configuration performance to establish the charge-discharge scheduling data result of energy storage of the smart grid.

[0271] The above is merely preferred embodiments of the present disclosure, rather than limiting the present disclosure in other forms. Those skilled in the art may alter or modify the technical content disclosed above into equivalent embodiments of equivalent variations and apply them to other fields. However, any simply amendment, equivalent change and modification made on the above embodiments according to the technical essence of the present disclosure without departing from the content of the technical solution of the present disclosure still fall within the protection scope of the technical solution of the present disclosure.