Lightning strike probability-based coordinated power control method for energy storage system and load

Abstract

A lightning strike probability-based coordinated power control method for an energy storage system and a load adjusts output power of an energy storage system and a regulatable load based on a real-time predicted lightning strike probability. This reduces a correlation between a power grid in a lightning region and an external power grid, minimizes impact of a lightning strike if any on the external power grid, and decreases an economic loss due to the lightning strike in the lightning region. The present disclosure establishes a relationship between the output power of the energy storage system, the regulated power of the regulatable load, and a lightning strike probability in a same region, and then regulates the power of the energy storage system and the regulatable load in real time based on the lightning strike probability.

Claims

1. A lightning strike probability-based coordinated power control method for an energy storage system and a load, wherein the method comprises: step 1: collecting a predicted lightning probability value for a current time interval ΔT, and performing step 2 when the predicted lightning probability value p.sub.lightning is >μ.sub.1, wherein μ.sub.1 represents a minimum predicted lightning probability value at which an energy storage system and a load can perform power regulation, and the time interval ΔT represents a lightning probability prediction cycle, which is equivalent to a regulation cycle of the energy storage system and in units of hour; step 2: determining a short-term scheduling cost loss of the energy storage system and a compensation cost function of regulated power of a regulatable load; step 3: establishing a power control model of the energy storage system based on a lightning strike probability; and step 4: establishing a power control model of the regulatable load based on the lightning strike probability, adjusting the regulated power of the regulatable load, and then returning to step 1 to collect the predicted lightning probability value for a next time interval and perform coordinated power control on the energy storage system and the load.

2. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 1, wherein in step 2, the short-term scheduling cost loss of the energy storage system satisfies the following relational expression: $J_{\mathrm{bat}}=\frac{{\mathrm{CC}}_{\mathrm{bat}}\times \Delta T\times P_{\mathrm{ba}}\times {\eta }_c}{N_{\mathrm{cycle}}}+{\mathrm{CC}}_b\times P_{\mathrm{ba}}^2\times \Delta T$ wherein J.sub.bat represents the short-term scheduling cost loss of the energy storage system, which is in units of Yuan; CC.sub.b represents a loss cost caused by a charging and discharging current, with a unit price of Yuan/(kW).sup.2 h; N.sub.cycle represents a quantity of cycles of the energy storage system, which is in units of times; P.sub.ba represents actual power of the energy storage system, which is in units of kW; CC.sub.bat represents an investment cost of the energy storage system, which is in units of Yuan/kWh; ΔT represents the lightning probability prediction cycle, which is equivalent to the regulation cycle of the energy storage system and in units of hour; and η.sub.c represents efficiency of the energy storage system.

3. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 2, wherein in step 2, the compensation cost function of the regulated power of the regulatable load is expressed as follows: $\{\begin{array}{c}{C}_{k,t}^{+}={a_k^{+}\left(\Delta P_{k,t}^{+}\right)}^2+b_k^{+}\Delta P_{k,t}^{+}-b_k^{+}\Delta P_{k,t}^{+}{\vartheta }_k^{+}\\ C_{k,t}^{-}={a_k^{-}\left(\Delta P_{k,t}^{-}\right)}^2+b_k^{-}\Delta P_{k,t}^{-}-b_k^{-}\Delta P_{k,t}^{-}{\vartheta }_k^{-}\end{array}$ wherein C.sub.k,t.sup.+ and C.sub.k,t.sup.− respectively represent a compensation cost caused by increasing power of the regulatable load and a compensation cost caused by decreasing the power of the regulatable load, which are in units of Yuan; a.sub.k.sup.+ represents a quadratic term coefficient of the compensation cost caused by increasing the power of the regulatable load; b.sub.k.sup.+ represents a single term coefficient of the compensation cost caused by increasing the power of the regulatable load; a.sub.k.sup.− represents a quadratic term coefficient of the compensation cost caused by decreasing the power of the regulatable load; b.sub.k.sup.− represents a single term coefficient of the compensation cost caused by decreasing the power of the regulatable load; and respectively represent parameters indicating willingness of a user for load adjustment when the power of the regulatable load is increased and decreased, both are greater than 0 and less than 1, and a greater value leads to stronger willingness of the user for increasing or reducing the power of the load; and ΔP.sub.k,t.sup.+ and ΔP.sub.k,t.sup.− respectively represent power variations when the power of the regulatable load is increased and decreased.

4. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 3, wherein step 3 comprises the following substeps: step 3.1: determining a possible loss of the system due to a lightning strike when power regulation of the energy storage system is not considered; step 3.2: determining a power loss of the system due to the lightning strike after the energy storage system preforms power regulation; step 3.3: determining optimal output power P*.sub.ba of the energy storage system; and step 3.4: adjusting output power of the energy storage system based on an operation constraint of the energy storage system.

5. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 4, wherein in step 3.1, the possible loss of the system due to the lightning strike when power regulation of the energy storage system is not considered satisfies the following relational expression:
Y.sub.1=P.sub.load×price×p.sub.lightning×ΔT wherein P.sub.load represents a power loss caused to the load by the lightning strike, which is in units of kW; and price represents a unit price of a lost electric quantity, which is in units of Yuan/kWh.

6. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 5, wherein in step 3.2, the power loss of the energy storage system due to the lightning strike after the energy storage system preforms power regulation satisfies the following relational expression: $Y_2=\left(P_{\mathrm{load}}-P_{\mathrm{ba}}\right)\times \mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T+J_{\mathrm{bat}}\left(P_{\mathrm{load}}-P_{\mathrm{ba}}\right)\times \mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T+\frac{{\mathrm{CC}}_{\mathrm{bat}}\times {\eta }_c\times \Delta T\times P_{\mathrm{ba}}}{N_{\mathrm{cycle}}}+{\mathrm{CC}}_b\times \Delta T\times {\left(P_{\mathrm{ba}}\right)}^2$ wherein P.sub.ba represents the output power of the energy storage system, which is in units of kW; and price represents the unit price of the lost electric quantity, which is in units of Yuan/kWh.

7. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 6, wherein in step 3.3, a difference between the loss when the energy storage system does not perform power regulation and the loss after the energy storage system performs power regulation satisfies the following relational expression: $\Delta Y=Y_2-Y_1={\mathrm{CC}}_b\times \Delta T\times {\left(P_{\mathrm{ba}}\right)}^2+\left(\frac{{\mathrm{CC}}_{\mathrm{bat}}\times {\eta }_c\times \Delta T}{N_{\mathrm{cycle}}}-\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T\right)\times P_{\mathrm{ba}}$ based on the difference ΔY between the loss when the energy storage system does not perform power regulation and the loss after the energy storage system performs power regulation, it is determined that the obtained optimal output power of the energy storage system is P*.sub.ba and satisfies the following relational expression: $P_{\mathrm{ba}}^{*}=\frac{\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T\times N_{\mathrm{cycle}}-{\mathrm{CC}}_{\mathrm{ba}}\times {\eta }_c\times \Delta T}{2\times {\mathrm{CC}}_b\times \Delta T\times N_{\mathrm{cycle}}}$ in other words, when ΔY is a minimum value, the optimal output power of the energy storage system is P*.sub.ba.

8. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 7, wherein in step 3.4, based on the operation constraint of the energy storage system, the output power of the energy storage system is adjusted as follows: $P_{\mathrm{ba}}=\{\begin{array}{cc}{P}_{\mathrm{ba}}^{*}& \left({\mathrm{SOC}}_{\min }\le \mathrm{SOC}\left(t\right)\right)\&\left(P_{\mathrm{ba}}^{*}\le P_{\mathrm{ba}\max }\right)\\ 0& \mathrm{SOC}\left(t\right)\le {\mathrm{SOC}}_{\min }\\ P_{\mathrm{ba}\max }& P_{\mathrm{ba}}^{*}\ge P_{\mathrm{ba}\max }\end{array}$ wherein the operation constraint of the energy storage system is as follows: $\{\begin{array}{c}{\mathrm{SOC}}_{\min }\le \mathrm{SOC}\left(t\right)\le {\mathrm{SOC}}_{\max }\\ P_{\mathrm{ba}\min }\le P_{\mathrm{ba}}\le P_{\mathrm{ba}\max }\end{array}$ wherein SOC.sub.min and SOC.sub.max respectively represent minimum and maximum allowable states of charge for operation of the energy storage system; SOC(t) represents an actual state of charge of the energy storage system; P.sub.bamin and P.sub.bamax respectively represent minimum and maximum allowable power for operation of the energy storage system, P.sub.bamin<0, and P.sub.bamax>0; and P.sub.ba represents the actual power of the energy storage system.

9. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 4, wherein step 4 comprises the following substeps: step 4.1: determining a possible loss of the system due to the lightning strike when the energy storage system participates in power regulation and the load does not perform power regulation; step 4.2: determining a possible loss of the energy storage system due to the lightning strike after both the energy storage system and the regulatable load preform power regulation; step 4.3: determining an optimal power regulation amount ΔP*.sub.l of the load; and step 4.4: regulating the regulated power of the regulatable load.

10. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 8, wherein step 4 comprises the following substeps: step 4.1: determining a possible loss of the system due to the lightning strike when the energy storage system participates in power regulation and the load does not perform power regulation; step 4.2: determining a possible loss of the energy storage system due to the lightning strike after both the energy storage system and the regulatable load preform power regulation; step 4.3: determining an optimal power regulation amount ΔP*.sub.l of the load; and step 4.4: regulating the regulated power of the regulatable load.

11. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 9, wherein in step 4.1, the possible loss of the system due to the lightning strike when the energy storage system participates in power regulation and the load does not perform power regulation satisfies the following relational expression:
X.sub.1=(P.sub.load−P*.sub.ba)×price×p.sub.lightning×ΔT wherein P.sub.load represents the power loss caused to the load by the lightning strike, which is in units of kW; price represents the unit price of the lost electric quantity, which is in units of Yuan/kWh; and P*.sub.ba represents the optimal output power of the energy storage system.

12. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 10, wherein in step 4.1, the possible loss of the system due to the lightning strike when the energy storage system participates in power regulation and the load does not perform power regulation satisfies the following relational expression:
X.sub.1=(P.sub.load−P*.sub.ba)×price×p.sub.lightning×ΔT wherein P.sub.load represents the power loss caused to the load by the lightning strike, which is in units of kW; price represents the unit price of the lost electric quantity, which is in units of Yuan/kWh; and P*.sub.ba represents the optimal output power of the energy storage system.

13. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 11, wherein in step 4.2, the possible loss of the energy storage system due to the lightning strike after both the energy storage system and the regulatable load preform power regulation satisfies the following relational expression:
X.sub.2=(P.sub.load−P*.sub.ba−ΔP.sub.l)×price×p.sub.lightning×ΔT+a.sub.k.sup.−(ΔP.sub.l).sup.2+b.sub.k.sup.−ΔP.sub.l−b.sub.k.sup.−ΔP.sub.l, wherein ΔP.sub.l represents a power regulation amount of the load.

14. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 12, wherein in step 4.2, the possible loss of the energy storage system due to the lightning strike after both the energy storage system and the regulatable load preform power regulation satisfies the following relational expression:
X.sub.2=(P.sub.load−P*.sub.ba−ΔP.sub.l)×price×p.sub.lightning×ΔT+a.sub.k.sup.−(ΔP.sub.l).sup.2+b.sub.k.sup.−ΔP.sub.l−b.sub.k.sup.−ΔP.sub.l, wherein ΔP.sub.l represents a power regulation amount of the load.

15. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 13, wherein in step 4.3, a function of a difference between the possible loss due to the lightning strike when the regulatable load does not perform regulation and the possible loss due to the lightning strike after the regulatable load performs regulation satisfies the following relational expression:
ΔX=X.sub.2−X.sub.1=a.sub.k.sup.−(ΔP.sub.l).sup.2+(b.sub.k.sup.−−b.sub.k.sup.−−price×p.sub.lightning×ΔT)ΔP.sub.l wherein a.sub.k.sup.− represents the quadratic term coefficient of the compensation cost caused by decreasing the power of the regulatable load; b.sub.k.sup.− represents the single term coefficient of the compensation cost caused by decreasing the power of the regulatable load; represents the parameter indicating the willingness of the user for load adjustment when the power of the regulatable load is decreased; and based on the function of the difference between the possible loss due to the lightning strike when the regulatable load does not perform regulation and the possible loss due to the lightning strike after the regulatable load performs regulation, the optimal power regulation amount ΔP*.sub.l of the load can be obtained and is expressed as follows: $\Delta P_l^{*}=\frac{b_k^{-}{\vartheta }_k^{-}+\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T-b_k^{-}}{2a_k^{-}}.$

16. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 14, wherein in step 4.3, a function of a difference between the possible loss due to the lightning strike when the regulatable load does not perform regulation and the possible loss due to the lightning strike after the regulatable load performs regulation satisfies the following relational expression:
ΔX=X.sub.2−X.sub.1=a.sub.k.sup.−(ΔP.sub.l).sup.2+(b.sub.k.sup.−−b.sub.k.sup.−−price×p.sub.lightning×ΔT)ΔP.sub.l wherein a.sub.k.sup.− represents the quadratic term coefficient of the compensation cost caused by decreasing the power of the regulatable load; b.sub.k.sup.− represents the single term coefficient of the compensation cost caused by decreasing the power of the regulatable load; represents the parameter indicating the willingness of the user for load adjustment when the power of the regulatable load is decreased; and based on the function of the difference between the possible loss due to the lightning strike when the regulatable load does not perform regulation and the possible loss due to the lightning strike after the regulatable load performs regulation, the optimal power regulation amount ΔP*.sub.l of the load can be obtained and is expressed as follows: $\Delta P_l^{*}=\frac{b_k^{-}{\vartheta }_k^{-}+\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T-b_k^{-}}{2a_k^{-}}.$

17. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 15, wherein in step 4.4, a method for regulating the regulated power of the regulatable load is as follows: setting a maximum power regulation amount of the load to ΔP.sub.lmax; and if ΔP*.sub.l≤ΔP.sub.lmax, determining that the regulated power of the regulatable load is ΔP*.sub.l; or if ΔP*.sub.l>ΔP.sub.lmax, determining that the regulated power of the regulatable load is ΔP.sub.lmax.

18. The lightning strike probability-based coordinated power control method for the energy storage system and the load according to claim 16, wherein in step 4.4, a method for regulating the regulated power of the regulatable load is as follows: setting a maximum power regulation amount of the load to ΔP.sub.lmax; and if ΔP*.sub.l≤ΔP.sub.lmax, determining that the regulated power of the regulatable load is ΔP*.sub.l; or if ΔP*.sub.l>ΔP.sub.lmax, determining that the regulated power of the regulatable load is ΔP.sub.lmax.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The sole FIGURE is a flowchart of a lightning strike probability-based coordinated power control method for an energy storage system and a load according to the present disclosure.

DETAILED DESCRIPTION

(2) The present disclosure is further described below with reference to the accompanying drawings. The following embodiments are only used for describing the technical solutions of the present disclosure more clearly, and are not intended to limit the protection scope of the present disclosure.

(3) The sole FIGURE is a flowchart of a lightning strike probability-based coordinated power control method for an energy storage system and a load according to the present disclosure. The method specifically includes the following content.

(4) Step 1: Collect a predicted lightning probability value p.sub.lightning for a current time interval ΔT, and perform step 2 when p.sub.lightning>μ.sub.1, where μ.sub.1 represents a minimum predicted lightning probability value at which an energy storage system and a load can perform power regulation.

(5) The time interval ΔT represents a lightning probability prediction cycle, which is equivalent to a regulation cycle of the energy storage system and in units of hour.

(6) When the predicted lightning probability value p.sub.lightning is ≤μ.sub.1, it indicates that there is no lightning; when the predicted lightning probability value p.sub.lightning is greater than μ.sub.1 and less than P.sub.2, it indicates that there is basically no lightning; when the predicted lightning probability value p.sub.lightning is equal to or greater than μ.sub.2 and less than μ.sub.3, it indicates that there may be lightning; when the predicted lightning probability value p.sub.lightning is ≥μ.sub.3, it indicates that there is lightning. Preferably, μ.sub.1 is 20%, μ.sub.2 is 50%, and μ.sub.3 is 70%.

(7) Step 2: Determine a short-term scheduling cost loss of the energy storage system and a compensation cost function of regulated power of a regulatable load.

(8) A regulation cost of the energy storage system depends on a quantity of cycles of the energy storage system. In addition, a discharging current also increases internal resistance of a battery. Based on the above two aspects, the short-term scheduling cost loss of the energy storage system can be determined according to the following formula:

(9) $J_{\mathrm{bat}}=\frac{C{C}_{\mathrm{bat}}\times \Delta T\times P_{ba}\times {\eta }_c}{N_{\mathrm{cycle}}}+C{C}_b\times P_{ba}^2\times \Delta T$

(10) In the above formula, J.sub.bat represents the short-term scheduling cost loss of the energy storage system, which is in units of Yuan; CC.sub.b represents a loss cost caused by a charging and discharging current, with a unit price of Yuan/(kW).sup.2 h; N.sub.cycle represents the quantity of cycles of the energy storage system, which is in units of times; P.sub.ba represents actual power of the energy storage system, which is in units of kW; CC.sub.bat represents an investment cost of the energy storage system, which is in units of Yuan/kWh; ΔT represents the lightning probability prediction cycle, which is equivalent to the regulation cycle of the energy storage system and in units of hour; and η.sub.c represents efficiency of the energy storage system.

(11) The compensation cost function of the regulated power of the regulatable load satisfies the following relational expression:

(12) The compensation cost function of the regulated power of the regulatable load is expressed as follows:

(13) 0$\{\begin{array}{c}{C}_{k,t}^{+}={a_k^{+}\left(\Delta P_{k,t}^{+}\right)}^2+b_k^{+}\Delta P_{k,t}^{+}-b_k^{+}\Delta P_{k,t}^{+}{\vartheta }_k^{+}\\ C_{k,t}^{-}={a_k^{-}\left(\Delta P_{k,t}^{-}\right)}^2+b_k^{-}\Delta P_{k,t}^{-}-b_k^{-}\Delta P_{k,t}^{-}{\vartheta }_k^{-}\end{array}.$

(14) In the above function, C.sub.k,t.sup.+ and C.sub.k,t.sup.− respectively represent a compensation cost caused by increasing power of the regulatable load and a compensation cost caused by decreasing the power of the regulatable load, which are in units of Yuan;

(15) a.sub.k.sup.+ represents a quadratic term coefficient of the compensation cost caused by increasing the power of the regulatable load;

(16) b.sub.k.sup.+ represents a single term coefficient of the compensation cost caused by increasing the power of the regulatable load;

(17) a.sub.k.sup.− represents a quadratic term coefficient of the compensation cost caused by decreasing the power of the regulatable load;

(18) b.sub.k.sup.− represents a single term coefficient of the compensation cost caused by decreasing the power of the regulatable load;

(19) + and respectively represent parameters indicating willingness of a user for load adjustment when the power of the regulatable load is increased and decreased, both are greater than 0 and less than 1, and a greater value leads to stronger willingness of the user for increasing or reducing the power of the load;

(20) those skilled in the art may set the willingness parameters based on actual industry characteristics, user characteristics and user load characteristics, where for an industry or a user who has no demand for power regulation of the system and the load, values of and are 0; for an industry such as a steelwork, the power of the regulatable load is generally not regulated based on system scheduling; for an industry using an electrolytic aluminum load, because a load power regulation capacity of the industry is greater than that of a common industry, the values of and are 0.5; for another industry or user who can regulate load power at any time, the values of and are greater than 0.5 and less than 1; and in this embodiment, the values of and are 0.8 for rural residents because level-3 loads are basically used for the rural residents, and short-term power outage can be carried out; and

(21) ΔP.sub.k,t.sup.+ and ΔP.sub.k,t.sup.− respectively represent power variations when the power of the regulatable load is increased and decreased.

(22) Step 3: Establish a power control model of the energy storage system based on a lightning strike probability.

(23) In order to ensure user satisfaction in electricity use, an energy storage device is preferentially used for power regulation to bear partial load power in a region with a lightning probability.

(24) Step 3.1: Determine a possible loss of the system due to a lightning strike when power regulation of the energy storage system is not considered.

(25) When the energy storage system does not participate in regulation, the possible loss due to the lightning strike is as follows:
Y.sub.1=P.sub.load×price×p.sub.lightning×ΔT

(26) In the above formula, P.sub.load represents a power loss caused to the load by the lightning strike, which is in units of kW; price represents a unit price of a lost electric quantity, which is in units of Yuan/kWh; and p.sub.lightning represents the lightning strike probability.

(27) Step 3.2: Determine a possible power loss of the system due to the lightning strike after the energy storage system preforms power regulation.

(28) After the energy storage system performs power regulation, the energy storage system bears partial power, and in this case, the possible loss due to the lightning strike is as follows:

(29) $Y_2=\left(P_{\mathrm{load}}-P_{\mathrm{ba}}\right)\times \mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T+J_{\mathrm{bat}}\left(P_{\mathrm{load}}-P_{\mathrm{ba}}\right)\times \mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T+\frac{{\mathrm{CC}}_{\mathrm{bat}}\times {\eta }_c\times \Delta T\times P_{\mathrm{ba}}}{N_{\mathrm{cycle}}}+{\mathrm{CC}}_b\times \Delta T\times {\left(P_{\mathrm{ba}}\right)}^2$

(30) In the above formula, P.sub.ba represents output power of the energy storage system, which is in units of kW; and price represents the unit price of the lost electric quantity, which is in units of Yuan/kWh.

(31) Step 3.3: Determine optimal output power P*.sub.ba of the energy storage system.

(32) An expected difference between the loss when the energy storage system does not perform power regulation and the loss after the energy storage system performs power regulation satisfies the following relational expression:

(33) $\Delta Y=Y_2-Y_1={\mathrm{CC}}_b\times \Delta T\times {\left(P_{\mathrm{ba}}\right)}^2+\left(\frac{{\mathrm{CC}}_{\mathrm{bat}}\times {\eta }_c\times \Delta T}{N_{\mathrm{cycle}}}-\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T\right)\times P_{\mathrm{ba}}$

(34) Apparently, when ΔY<0, in other words, when ΔY is a minimum value and a reduced loss is a maximum value, the optimal output power of the energy storage system is P*.sub.ba and can be determined based on the expected difference ΔY between the loss when the energy storage system does not perform power regulation and the loss after the energy storage system performs power regulation, which satisfies the following relational expression:

(35) $P_{\mathrm{ba}}^{*}=\frac{\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T\times N_{\mathrm{cycle}}-{\mathrm{CC}}_{\mathrm{ba}}\times {\eta }_c\times \Delta T}{2\times {\mathrm{CC}}_b\times \Delta T\times N_{\mathrm{cycle}}}$

(36) Step 3.4: Adjust the output power of the energy storage system based on an operation constraint of the energy storage system.

(37) $\{\begin{array}{c}{\mathrm{SOC}}_{\min }\le \mathrm{SOC}\left(t\right)\le {\mathrm{SOC}}_{\max }\\ P_{\mathrm{ba}\min }\le P_{\mathrm{ba}}\le P_{\mathrm{ba}\max }\end{array}$

(38) In the above operation constraint of the energy storage system, SOC.sub.min and SOC.sub.max respectively represent minimum and maximum allowable states of charge for operation of the energy storage system; SOC(t) represents an actual state of charge of the energy storage system; P.sub.bamin and P.sub.bamax respectively represent minimum and maximum allowable power for operation of the energy storage system, P.sub.bamin<0, and P.sub.bamax>0; and P.sub.ba represents the actual power of the energy storage system, and in this specification, P.sub.ba>0 because discharging is needed.

(39) After the above constraint is considered, the output power of the energy storage system can be obtained according to the following formula:

(40) $P_{\mathrm{ba}}=\{\begin{array}{cc}{P}_{\mathrm{ba}}^{*}& \left({\mathrm{SOC}}_{\min }\le \mathrm{SOC}\left(t\right)\right)\&\left(P_{\mathrm{ba}}^{*}\le P_{\mathrm{ba}\max }\right)\\ 0& \mathrm{SOC}\left(t\right)\le {\mathrm{SOC}}_{\min }\\ P_{\mathrm{ba}\max }& P_{\mathrm{ba}}^{*}\ge P_{\mathrm{ba}\max }\end{array}$

(41) Step 4: Establish a power control model of the regulatable load based on the lightning strike probability, adjust the regulated power of the regulatable load, and then return to step 1 to collect the predicted lightning probability value for a next time interval and perform coordinated power control on the energy storage system and the load.

(42) As a virtual energy storage system, the regulatable load can also have an effect similar to an effect of the energy storage system. On a premise that the energy storage system participates in regulation, the load power can also be regulated to further improve the initiative of the system, reduce a correlation between a region with the lightning probability and an external power grid, and reduce a possible loss caused by lightning.

(43) Step 4.1: Determine a possible loss of the system due to the lightning strike when the energy storage system participates in power regulation and the load does not perform power regulation.

(44) When the battery participates in regulation and the load does not participate in regulation, the power loss due to the lightning strike can be determined according to the following formula:
X.sub.1=(P.sub.load−P*.sub.ba)×price×p.sub.lightning×ΔT

(45) Step 4.2: Determine a possible loss of the system due to the lightning strike after both the energy storage system and the regulatable load preform power regulation.

(46) After the load participates in regulation, the possible loss due to the lightning strike can be calculated according to the following formula:
X.sub.2=(P.sub.load−P*.sub.ba−ΔP.sub.l)×price×p.sub.lightning×ΔT+a.sub.k.sup.−(ΔP.sub.l).sup.2+b.sub.k.sup.−ΔP.sub.l−b.sub.k.sup.−ΔP.sub.l,

(47) In the above formula, ΔP.sub.l represents a power regulation amount of the load.

(48) Step 4.3: Determine an optimal power regulation amount ΔP*.sub.l of the load.

(49) A function of a difference between the possible loss due to the lightning strike when the regulatable load does not perform regulation and the possible loss due to the lightning strike after the regulatable load performs regulation can be expressed as follows:
ΔX=X.sub.2−X.sub.1=a.sub.k.sup.−(ΔP.sub.l).sup.2+(b.sub.k.sup.−−b.sub.k.sup.−−−price×p.sub.lightning×ΔT)ΔP.sub.l

(50) In the above expression, a.sub.k.sup.− represents the quadratic term coefficient of the compensation cost caused by decreasing the power of the regulatable load;

(51) b.sub.k.sup.− represents the single term coefficient of the compensation cost caused by decreasing the power of the regulatable load; and

(52) represents the parameter indicating the willingness of the user for load adjustment when the power of the regulatable load is decreased.

(53) Apparently, based on the function of the difference between the possible loss due to the lightning strike when the regulatable load does not perform regulation and the possible loss due to the lightning strike after the regulatable load performs regulation, the optimal power regulation amount of the load can be obtained, and may be expressed as follows:

(54) $\Delta P_l^{*}=\frac{b_k^{-}{\vartheta }_k^{-}+\mathrm{price}\times p_{\mathrm{lightning}}\times \Delta T-b_k^{-}}{2a_k^{-}}$

(55) In this case, after the regulatable load participates in regulation, the possible loss caused by lightning can be further reduced.

(56) Step 4.4: Adjust the regulated power of the regulatable load based on a power regulation threshold of the load.

(57) A maximum power regulation amount of the load is ΔP.sub.lmax. if ΔP*.sub.l≤ΔP.sub.lmax, the regulated power of the regulatable load is ΔP*.sub.l; if ΔP*.sub.l>ΔP.sub.lmax, the regulated power of the regulatable load is ΔP.sub.lmax.

(58) In this case, optimal regulated power of the regulatable load can be determined.

(59) After the current regulation cycle, a next regulation cycle is followed, and time t correspondingly increases by ΔT.

(60) The applicant of the present disclosure has made a detailed description of the implementation examples of the present disclosure with reference to the accompanying drawings in the specification. However, those skilled in the art should understand that the above implementation examples are only preferred implementation solutions of the present disclosure, and the detailed description is only to help readers better understand the spirit of the present disclosure, rather than to limit the protection scope of the present disclosure. On the contrary, any improvement or modification based on the spirit of the present disclosure shall fall within the protection scope of the present disclosure.