阅读说明：本技术 一种多智能家庭的分布式智能能量管理方法 (Distributed intelligent energy management method for multiple intelligent families ) 是由 胡晓松 邓昕晨 邓忠伟 李佳承 胡兵 徐立军 王小娟 于 2021-09-16 设计创作，主要内容包括：本发明涉及一种多智能家庭的分布式智能.能量管理方法,属于智能电网技术领域,包括以下步骤：S1：建立智能家庭各部件的数学模型；建立社区中各智能家庭之间的拓扑结构,以实现其能量交易；制定各智能家庭之间的能量交易规则；S2：针对各智能家庭中的电动汽车建立相应的动力电池老化模型,包括循环老化模型和日历老化模型；S3：针对各智能家庭建立其日常电能成本最低的优化目标,设计各部件相应的约束；S4：对原优化问题转换成凸优化问题,保证所得解为全局最优；S5：设计各智能家庭之间能量交易的一致性辅助变量,利用Consensus ADMM算法对各智能家庭的优化问题进行分解协调至最优,使整个社区的电能成本最低。(The invention relates to a distributed intelligence energy management method for multiple intelligent families, which belongs to the technical field of intelligent power grids and comprises the following steps: s1: establishing a mathematical model of each part of the intelligent home; establishing a topological structure among intelligent families in a community to realize energy trading; formulating an energy transaction rule among all intelligent families; s2: establishing corresponding power battery aging models including a circulating aging model and a calendar aging model for the electric automobiles in each intelligent family; s3: establishing an optimization target with the lowest daily electric energy cost aiming at each intelligent home, and designing corresponding constraints of each part; s4: converting the original optimization problem into a convex optimization problem, and ensuring that the obtained solution is globally optimal; s5: and designing a consistency auxiliary variable of energy transaction among the intelligent families, and decomposing and coordinating the optimization problem of each intelligent family to be optimal by using a Consensus ADMM algorithm so as to minimize the electric energy cost of the whole community.)

1. A distributed intelligence energy management method for multiple intelligent families is characterized by comprising the following steps: the method comprises the following steps:

s1: establishing a mathematical model of each part of the intelligent home; establishing a topological structure among intelligent families in a community to realize energy trading; formulating an energy transaction rule among all intelligent families;

s2: establishing corresponding power battery aging models including a circulating aging model and a calendar aging model for the electric automobiles in each intelligent family;

s3: establishing an optimization target with the lowest daily electric energy cost aiming at each intelligent home, and designing corresponding constraints of each part;

s4: converting the original optimization problem into a convex optimization problem, and ensuring that the obtained solution is globally optimal;

s5: and designing a consistency auxiliary variable of energy transaction among the intelligent families, and decomposing and coordinating the optimization problem of each intelligent family to be optimal by using a Consensus ADMM algorithm so as to minimize the electric energy cost of the whole community.

2. The method for distributed intelligence energy management of a multi-intelligent home of claim 1, wherein: in the step S1, each part of the intelligent home comprises a power utilization and generation prediction model, an energy storage system and an electric automobile; the step S1 specifically includes the following steps:

s11: establishing a mathematical model of each component of the intelligent home, wherein the mathematical model of the energy storage system is as follows:

wherein t is the time when the system operates,is the energy contained in the energy storage system, P_i ^ch(t) and P_i ^dis(T) is the charge and discharge power, eta is the charge and discharge efficiency, and delta T is the operation step length;

the mathematical model of the power battery of the electric automobile is as follows:

wherein, P_i ^ev(t) is the output power of the power battery at the time t,the open rate voltage of the battery related to the SOC of the battery at different moments, I (t) is charging and discharging current, R is internal resistance of the battery, and N is the number of single batteries in the electric automobile;

s12: establishing a topological structure among intelligent families in a community, and designing circuit connection among the families;

s13: and formulating energy transaction rules and prices among the intelligent families.

3. The method for distributed intelligence energy management of a multi-intelligent home of claim 1, wherein: the step S2 specifically includes the following steps:

s21: the method comprises the following steps of selecting an electric automobile power battery aging model, wherein the electric automobile power battery aging model comprises a cyclic aging model and a calendar aging model, the cyclic aging of a battery is related to a current multiplying power c, a temperature T and an ampere hour throughput Ah, the calendar aging is related to a battery state of charge (SOC), a temperature T and a storage time T, and the expressions are respectively as follows:

wherein the content of the first and second substances,andrespectively, capacity loss due to cyclic charge and discharge and battery capacity loss due to storage aging, B_cycAnd B_calExponential pre-factors, Ea, for the cyclic aging model and the calendar aging model, respectively_cycAnd Ea_calActivation energies of a cyclic aging model and a calendar aging model, respectively, and are related to a current rate and a battery SOC, respectively, R is an ideal gas constant, z_cycAnd z_calIs a dimensionless constant;

s22: the environmental temperature of the electric automobile is kept at 25 ℃, a cyclic aging model and a calendar aging model of the battery are converted into a cyclic aging capacity loss and a calendar aging capacity loss in unit time, which are respectively functions related to charge and discharge power and the SOC of the battery, and the expressions are as follows:

wherein, P_bFor the charge and discharge power of the battery, t is the time of system operation, N is the number of battery cells of the electric automobile, E_bThe rated energy of the battery;

s23: the two formulas in step S22 are respectively fitted with piecewise linear functions, as follows:

wherein, a_0,1、a_1,1、a_0,2、a_1,2、a_0,3、a_1,3Respectively representing the coefficients of a linear piecewise function fitted according to the cyclic aging variation of the battery per unit time, b_0,1、b_1,1、b_0,2、b_1,2、b_0,3、b_1,3Respectively, represent linear piecewise function coefficients fitted according to the calendar aging variation of the battery per unit time.

4. The method for distributed intelligence energy management of a multi-intelligent home of claim 1, wherein: the step S3 specifically includes the following steps:

establishing an objective function and related constraints of the minimum daily electric energy cost of each intelligent home, including energy trading among the intelligent homes, and optimizing the grid power utilization function of the ith intelligent home at the moment tRate of changePower sold to the gridCharging and discharging power of energy storage batteryAndpower purchased to Smart Home jAnd power of saleThe expression is as follows:

wherein the objective function is the cost of the family i in the time domain T, including the expenditure of the power gridRevenue for selling power to the gridCost of energy storage agingAging cost of power battery of electric automobileCost of purchasing power to family jAnd revenue of vending electric energy

5. The method for distributed intelligence energy management of a multi-intelligent home of claim 1, wherein: the step S4 specifically includes an equation constrained convex form and a model convex form of the electric vehicle power battery, which is as follows:

wherein, I_tAnd c and d are coefficients of linear functions for the current of the power battery at the moment t.

6. The method for distributed intelligence energy management of a multi-intelligent home of claim 1, wherein: the step S5 specifically includes the following steps:

s51: design uniformity auxiliary variablesAndcoordinating a plurality of intelligent families to enable the community to reach a consensus as a whole, wherein the consistency auxiliary variable is the electric energy buying and selling between the ith family and the jth family, and the following formula is adopted:

wherein k is the number of iterations;

s52: decomposing the optimization problem with the lowest energy cost of the whole community into a plurality of sub-problems by using a Consensus ADMM algorithm, and leading the energy cost of the community to be the lowest by coordinating the whole to reach the agreement, wherein the expression of the sub-optimization problem is as follows:

wherein N is_iP is a parameter predefined in the multiplier method for the total number of households in which electric energy is traded with the household i,for scaled dual variables, λ_ijIs a lagrange multiplier, k is the number of iterations,a dual variable scaled by a consistency constraint for the purchase of power by family i to family j at time t,a dual variable scaled by a consistency constraint for selling electric energy from family i to family j at time t;

s53: in each iteration process, after the calculation of the sub-optimization problem is completed, the scaled dual variable is updated, and the expression is as follows:

s54: in order to ensure that the Consensus ADMM algorithm coordinates each sub-optimization problem to be optimal, the following convergence standard is established:

wherein the content of the first and second substances,for the residual error of the original optimization problem,for dual problem residual, epsilon₁And ε₂Set threshold value, Z_ij,kA consistency assist variable for energy trading between family i and family j for low k iterations.

Technical Field

The invention belongs to the technical field of smart power grids, and relates to a distributed intelligent energy management method for multiple intelligent families.

Background

The community containing a plurality of intelligent families is an indispensable component of the intelligent power grid, can provide energy autonomous functional modules in local areas, and is the development trend of the intelligent power grid. The intelligent home comprises a controllable load device, an uncontrollable load device, an Electric Vehicle (EV), Renewable Energy Sources (RESs) and an Energy Storage System (ESS), and high Energy autonomy can be realized by managing the EV, so that the Energy cost is greatly reduced. At present, distributed management of a plurality of intelligent families becomes a research hotspot, the intelligent families are both energy suppliers and energy consumers, and the multilateral trading mechanism can promote the development of energy autonomous communities.

Traditional centralized management methods require a high bandwidth communication infrastructure and all communication information must be collected and processed by a central controller. Furthermore, for large communities, the related optimization problem may be computationally difficult to solve. Distributed algorithms can decompose the entire optimization problem into smaller scale local problems in a distributed manner, requiring only limited information to be exchanged with neighboring users. The distributed management mode has the advantages of expandability, flexibility, robustness and the like, and the parallel computing capacity can accelerate the convergence speed.

Disclosure of Invention

In view of this, the present invention provides an advanced multi-smart home distributed energy management method, which ensures extensibility, flexibility, robustness and high efficiency of management.

In order to achieve the purpose, the invention provides the following technical scheme:

the energy management method comprises the following steps:

s1: establishing a mathematical model of each part of the intelligent home; establishing a topological structure among intelligent families in a community to realize energy trading; formulating an energy transaction rule among all intelligent families;

s2: establishing corresponding power battery aging models including a circulating aging model and a calendar aging model for the electric automobiles in each intelligent family;

s3: establishing an optimization target with the lowest daily electric energy cost aiming at each intelligent home, and designing corresponding constraints of each part;

s4: converting the original optimization problem into a convex optimization problem, and ensuring that the obtained solution is globally optimal;

s5: and designing a consistency auxiliary variable of energy transaction among the intelligent families, and decomposing and coordinating the optimization problem of each intelligent family to be optimal by using a Consensus ADMM algorithm so as to minimize the electric energy cost of the whole community.

Further, in the step S1, each component of the smart home includes a power consumption and generation prediction model, an energy storage system, and an electric vehicle; the step S1 specifically includes the following steps:

s11: the method comprises the following steps of establishing mathematical models of all parts of the intelligent home to realize intelligent control of all the parts in the intelligent home and reduce the electric energy cost, wherein the mathematical models of the energy storage system are as follows:

wherein t is the time when the system operates,is the energy contained in the energy storage system, P_i ^ch(t) and P_i ^dis(T) is the charge and discharge power, eta is the charge and discharge efficiency, and delta T is the operation step length.

The mathematical model of the power battery of the electric automobile is as follows:

wherein, P_i ^ev(t) is the output power of the power battery at the time t,the open rate voltage of the battery is related to the SOC of the battery at different moments, I (t) is charging and discharging current, R is internal resistance of the battery, and N is the number of single batteries in the electric automobile.

S12: establishing a topological structure among intelligent families in a community, and designing circuit connection among the intelligent families to realize energy trading;

s13: energy trading rules and prices among the intelligent families are formulated, redundant power generation is promoted to be sold to users with high power utilization requirements, electric energy purchased to a power grid is reduced, and therefore the power utilization cost of the whole community is reduced.

Further, the step S2 specifically includes the following steps:

s21: the method comprises the steps of selecting an electric automobile power battery aging model, wherein the electric automobile power battery aging model comprises a cyclic aging model and a calendar aging model, the cyclic aging of a battery is related to a current multiplying power c, a temperature T and an ampere hour throughput Ah, the calendar aging is related to a battery State of Charge (SOC), the temperature T and a storage time T, and the expressions are as follows:

wherein the content of the first and second substances,andrespectively, capacity loss due to cyclic charge and discharge and battery capacity loss due to storage aging, B_cycAnd B_calExponential pre-factors, Ea, for the cyclic aging model and the calendar aging model, respectively_cycAnd Ea_calActivation energies of a cyclic aging model and a calendar aging model, respectively, and are related to a current rate and a battery SOC, respectively, R is an ideal gas constant, z_cycAnd z_calIs a dimensionless constant;

s22: the environmental temperature of the electric automobile is kept at 25 ℃, a cyclic aging model and a calendar aging model of the battery are converted into a cyclic aging capacity loss and a calendar aging capacity loss in unit time, which are respectively functions related to charge and discharge power and the SOC of the battery, and the expressions are as follows:

wherein, P_bFor the charge and discharge power of the battery, t is the time of system operation, N is the number of battery cells of the electric automobile, E_bThe rated energy of the battery;

s23: the two formulas in step S22 are respectively fitted with piecewise linear functions, as follows:

wherein, a_0,1、a_1,1、a_0,2、a_1,2、a_0,3、a_1,3Respectively representing the coefficients of a linear piecewise function fitted according to the cyclic aging variation of the battery per unit time, b_0,1、b_1,1、b_0,2、b_1,2、b_0,3、b_1,3Respectively, represent linear piecewise function coefficients fitted according to the calendar aging variation of the battery per unit time.

Further, the step S3 specifically includes the following steps:

establishing an objective function and related constraints of the minimum daily electric energy cost of each intelligent home, including energy trading among the intelligent homes, and optimizing the power consumption of the power grid of the ith intelligent home at the moment tPower sold to the gridCharging and discharging power of energy storage batteryAndpower purchased to Smart Home jAnd power of saleThe expression is as follows:

wherein the objective function is the cost of the family i in the time domain T, including the expenditure of the power gridRevenue for selling power to the gridCost of energy storage agingAging cost of power battery of electric automobileCost of purchasing power to family jAnd revenue of vending electric energy

Further, the step S4 specifically includes an equation constrained convex form and a model convex form of the electric vehicle power battery, as follows:

wherein, I_tAnd c and d are coefficients of linear functions for the current of the power battery at the moment t.

Further, the step S5 specifically includes the following steps:

s51: design uniformity auxiliary variablesAndcoordinating a plurality of intelligent families to enable the community to reach a consensus as a whole, wherein the consistency auxiliary variable is the electric energy buying and selling between the ith family and the jth family, and the following formula is adopted:

s52: decomposing the optimization problem with the lowest energy cost of the whole community into a plurality of sub-problems by using a Consensus ADMM algorithm, and leading the energy cost of the community to be the lowest by coordinating the whole to reach the agreement, wherein the expression of the sub-optimization problem is as follows:

wherein N is_iP is a parameter predefined in the multiplier method for the total number of households in which electric energy is traded with the household i,for scaled dual variables, λ_ijIs a lagrange multiplier, k is the number of iterations,for family i to family jthe consistency of the purchased power at time t constrains the scaled dual variables,a dual variable scaled by a consistency constraint for selling electric energy from family i to family j at time t;

s53: in each iteration process, after the calculation of the sub-optimization problem is completed, the scaled dual variable is updated, and the expression is as follows:

s54: in order to ensure that the Consensus ADMM algorithm coordinates each sub-optimization problem to be optimal, the following convergence standard is established:

wherein the content of the first and second substances,for the residual error of the original optimization problem,for dual problem residual, epsilon₁And ε₂Set threshold value, Z_ij,kA consistency assist variable for energy trading between family i and family j for low k iterations.

The invention has the beneficial effects that: the method utilizes electrochemical impedance spectrum analysis to establish a fractional order equivalent circuit model with an inductive element for accurately simulating the high nonlinearity of the battery, utilizes Grunnwald-Letnikov (GL) definition in fractional order calculus to carry out modeling to obtain an improved fractional order model, and then combines a fractional order extended Kalman algorithm to realize the high-precision SOC estimation of the battery.

The advantages of the invention include the following:

(1) traditional centralized energy management approaches for multiple smart homes require a high bandwidth communication infrastructure because all the detailed information needs to be collected by the central controller. In contrast, in the distributed energy management approach proposed by the present invention, the small-scale local problem of a single subsystem is solved in a distributed manner, requiring only limited information to be exchanged with neighboring subsystems and/or a central controller during the optimization process.

(2) In the distributed energy management method proposed by the present invention, a single subsystem does not need to communicate its confidential information to other subsystems and the central controller.

(3) Compared with centralized management, the distributed management method provided by the invention has higher expandability and flexibility in system change, and particularly takes the topological structures of the power grid and communication infrastructure into consideration, so that the distributed management method is more dynamic in the smart power grid.

(4) Distributed management is more robust than centralized management. In centralized optimization, the entire system may be interrupted while the central controller is offline, whereas the distributed management method proposed by the present invention can be solved asynchronously by a single local controller, i.e. some local controllers may continue their normal function due to the loss of other local controllers. Furthermore, once the lost local controller is brought back online, accurate results can ultimately be obtained.

(5) Compared with a centralized management method, the distributed management method has higher efficiency. In the face of a large system or community, the centralized management method is low in calculation efficiency of large data and even cannot be solved due to the limitation of calculation capacity. The distributed management method provided by the invention can decompose a large optimization problem into a plurality of small problems to be solved respectively, so that the calculation efficiency is greatly improved, and the calculation complexity is equivalent to a small optimization problem by virtue of the parallel calculation capability.

Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention may be realized and attained by the means of the instrumentalities and combinations particularly pointed out hereinafter.

Drawings

For the purposes of promoting a better understanding of the objects, aspects and advantages of the invention, reference will now be made to the following detailed description taken in conjunction with the accompanying drawings in which:

FIG. 1 is a flow chart of a distributed energy management method for a multi-smart home in accordance with the present invention;

FIG. 2 is an EV power battery model;

FIG. 3 is a schematic diagram of a multi-smart-home topology of the present invention;

FIG. 4 is a fitting graph after adjustment of an aging model of an EV power battery;

FIG. 5 is a schematic diagram of the distributed management algorithm Consensus ADMM message communication in the present invention;

fig. 6 is a schematic diagram of the iterative process of Consensus ADMM in the present invention.

Detailed Description

The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention in a schematic way, and the features in the following embodiments and examples may be combined with each other without conflict.

Wherein the showings are for the purpose of illustrating the invention only and not for the purpose of limiting the same, and in which there is shown by way of illustration only and not in the drawings in which there is no intention to limit the invention thereto; to better illustrate the embodiments of the present invention, some parts of the drawings may be omitted, enlarged or reduced, and do not represent the size of an actual product; it will be understood by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted.

The same or similar reference numerals in the drawings of the embodiments of the present invention correspond to the same or similar components; in the description of the present invention, it should be understood that if there is an orientation or positional relationship indicated by terms such as "upper", "lower", "left", "right", "front", "rear", etc., based on the orientation or positional relationship shown in the drawings, it is only for convenience of description and simplification of description, but it is not an indication or suggestion that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and therefore, the terms describing the positional relationship in the drawings are only used for illustrative purposes, and are not to be construed as limiting the present invention, and the specific meaning of the terms may be understood by those skilled in the art according to specific situations.

Referring to fig. 1, the distributed intelligent energy management method for multiple intelligent homes according to the present invention includes the following steps:

step S1: establishing mathematical models of all parts (a power utilization and generation prediction model, an energy storage system and an electric automobile) of the intelligent home; establishing a topological structure among intelligent families in a community to realize energy trading; formulating an energy transaction rule among all intelligent families;

step S2: establishing corresponding power battery aging models including a circulating aging model and a calendar aging model for the electric automobiles in each intelligent family;

step S3: establishing an optimization target with the lowest daily electric energy cost aiming at each intelligent home, and designing corresponding constraints of each part;

step S4: converting the original optimization problem into a convex optimization problem, and ensuring that the obtained solution is globally optimal;

step S5: and designing a consistency auxiliary variable of energy transaction among the intelligent families, and decomposing and coordinating the optimization problem of each intelligent family to be optimal by using a Consensus ADMM algorithm so as to minimize the electric energy cost of the whole community.

In an embodiment of the present invention, a distributed intelligent energy management method for multiple intelligent homes specifically includes the following steps:

the step S1 specifically includes steps S11 to S13:

step S11: the mathematical model of each part of the intelligent family (electricity utilization and power generation prediction model, energy storage system and electric automobile) is established to realize the intelligent control of each part in the intelligent family, reduce its electric energy cost, wherein the mathematical model of energy storage system and electric automobile power battery respectively as follows:

wherein t is the time when the system operates,is the energy contained in the energy storage system, P_i ^ch(t) and P_i ^dis(T) is the charge and discharge power, eta is the charge and discharge efficiency, and delta T is the operation step length. Wherein, P_i ^ev(t) is the output power of the power battery at the time t,the open rate voltage of the battery is related to the SOC of the battery at different moments, I (t) is charging and discharging current, R is internal resistance of the battery, and N is the number of single batteries in the electric automobile.

Please refer to fig. 2 for a mathematical model of the power battery of the electric vehicle.

Step S12: establishing a topological structure among intelligent households in a community, and designing circuit connection among the households to realize energy trading, wherein a specific topological structure schematic diagram of the multiple intelligent households refers to fig. 3;

step S13: energy trading rules and prices among the intelligent families are formulated, redundant power generation is promoted to be sold to users with high power utilization requirements, electric energy purchased to a power grid is reduced, and therefore the power utilization cost of the whole community is reduced.

The step S2 specifically includes S21 to S23:

step S21: the method comprises the steps of selecting an electric automobile power battery aging model, wherein the electric automobile power battery aging model comprises a cyclic aging model and a calendar aging model, the cyclic aging of a battery is related to a current multiplying power c, a temperature T and an ampere hour throughput Ah, the calendar aging is related to a battery State of Charge (SOC), the temperature T and a storage time T, and the expressions are as follows:

wherein the content of the first and second substances,andrespectively, capacity loss due to cyclic charge and discharge and battery capacity loss due to storage aging, B_cycAnd B_calExponential pre-factors, Ea, for the cyclic aging model and the calendar aging model, respectively_cycAnd Ea_calActivation energies of a cyclic aging model and a calendar aging model, respectively, and are related to a current rate and a battery SOC, respectively, R is an ideal gas constant, z_cycAnd z_calAre dimensionless constants.

Step S22: assuming that the ambient temperature of the electric vehicle is constantly 25 ℃, converting the cyclic aging model and the calendar aging model of the battery into a cyclic aging capacity loss and a calendar aging capacity loss per unit time, which are functions related to the charge and discharge power and the battery SOC, respectively, and expressing the following:

wherein, P_bFor the charge and discharge power of the battery, t is the time of system operation, N is the number of battery cells of the electric automobile, E_bThe rated energy of the battery.

Step S23: fitting the two formulas with piecewise linear functions respectively as follows:

wherein, a_0,1、a_1,1、a_0,2、a_1,2、a_0,3、a_1,3Respectively representing the coefficients of a linear piecewise function fitted according to the cyclic aging variation of the battery per unit time, b_0,1、b_1,1、b_0,2、b_1,2、b_0,3、b_1,3Respectively, represent linear piecewise function coefficients fitted according to the calendar aging variation of the battery per unit time. See fig. 4 for the effect after fitting.

Step S3 specifically includes:

establishing an objective function and related constraints of the minimum daily electric energy cost of each intelligent home, including energy trading among the intelligent homes, and optimizing the power consumption of the power grid of the ith intelligent home at the moment tPower sold to the gridCharging and discharging power of energy storage batteryAndpower purchased to Smart Home jAnd power of saleThe expression is as follows:

wherein the objective function is the cost of the family i in the time domain T, including the expenditure of the power gridRevenue for selling power to the gridCost of energy storage agingAging cost of power battery of electric automobileCost of purchasing power to family jAnd revenue of vending electric energy

Step S4 specifically includes:

converting the original optimization problem into a convex optimization problem, and enabling the obtained solution to be a global optimal solution which comprises a convex form of equality constraint and a model convex form of an electric automobile power battery, wherein the model convex form comprises the following formula:

wherein, I_tAnd c and d are coefficients of linear functions for the current of the power battery at the moment t.

The step S5 specifically includes S51 to S54:

step S51: design uniformity auxiliary variablesAndcoordinating a plurality of intelligent families to enable the community to reach a consensus as a whole, wherein the consistency auxiliary variable is the electric energy buying and selling between the ith family and the jth family, and the following formula is adopted:

step S52: the Consensus ADMM algorithm is used to decompose the optimization problem with the lowest energy cost in the whole community into a plurality of sub-problems, the energy cost in the community is minimized by coordinating the whole to reach the Consensus, the data communication mode refers to fig. 5, and the expression of the sub-optimization problem is as follows:

wherein N is_iP is a parameter predefined in the multiplier method for the total number of households in which electric energy is traded with the household i,for scaled dual variables, λ_ijIs a lagrange multiplier, k is the number of iterations,a dual variable scaled by a consistency constraint for the purchase of power by family i to family j at time t,a dual variable scaled by a consistency constraint for selling electric energy from family i to family j at time t;

step S53: referring to fig. 6, in each iteration process, after the sub-optimization problem is calculated, the scaled dual variable is updated, and the expression is as follows:

step S54: in order to ensure that the Consensus ADMM algorithm coordinates each sub-optimization problem to be optimal, the following convergence standard is established:

wherein the former is the residual error of the original optimization problem, and the latter is the residual error of the dual problem, epsilon₁And ε₂Set threshold value, Z_ij,kA consistency assist variable for energy trading between family i and family j for low k iterations.

Finally, the above embodiments are only intended to illustrate the technical solutions of the present invention and not to limit the present invention, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions, and all of them should be covered by the claims of the present invention.