CN115081758A - Calculation transfer demand response system oriented to coordination data center and power grid - Google Patents

Calculation transfer demand response system oriented to coordination data center and power grid Download PDF

Info

Publication number
CN115081758A
CN115081758A CN202211002763.0A CN202211002763A CN115081758A CN 115081758 A CN115081758 A CN 115081758A CN 202211002763 A CN202211002763 A CN 202211002763A CN 115081758 A CN115081758 A CN 115081758A
Authority
CN
China
Prior art keywords
data center
cluster
ith
center cluster
energy storage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211002763.0A
Other languages
Chinese (zh)
Other versions
CN115081758B (en
Inventor
汤瑞欣
陈忠颖
李俊华
黄佳允
梁志豪
欧韦宜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhaoqing Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Zhaoqing Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhaoqing Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Zhaoqing Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority to CN202211002763.0A priority Critical patent/CN115081758B/en
Publication of CN115081758A publication Critical patent/CN115081758A/en
Application granted granted Critical
Publication of CN115081758B publication Critical patent/CN115081758B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46Multiprogramming arrangements
    • G06F9/50Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5061Partitioning or combining of resources
    • G06F9/5072Grid computing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46Multiprogramming arrangements
    • G06F9/50Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5094Allocation of resources, e.g. of the central processing unit [CPU] where the allocation takes into account power or heat criteria
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q30/00Commerce
    • G06Q30/02Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201Market modelling; Market analysis; Collecting market data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q30/00Commerce
    • G06Q30/06Buying, selling or leasing transactions
    • G06Q30/0601Electronic shopping [e-shopping]
    • G06Q30/0611Request for offers or quotes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q40/00Finance; Insurance; Tax strategies; Processing of corporate or income taxes
    • G06Q40/04Trading; Exchange, e.g. stocks, commodities, derivatives or currency exchange
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00006Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment
    • H02J13/00016Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment using a wired telecommunication network or a data transmission bus
    • H02J13/00017Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment using a wired telecommunication network or a data transmission bus using optical fiber
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/007Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
    • H02J3/0075Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2310/00The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/50The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads
    • H02J2310/56The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads characterised by the condition upon which the selective controlling is based
    • H02J2310/62The condition being non-electrical, e.g. temperature
    • H02J2310/64The condition being economic, e.g. tariff based load management

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Finance (AREA)
  • General Business, Economics & Management (AREA)
  • Accounting & Taxation (AREA)
  • Power Engineering (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Tourism & Hospitality (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Data Mining & Analysis (AREA)
  • Educational Administration (AREA)
  • Mathematical Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Medical Informatics (AREA)
  • Computing Systems (AREA)
  • Technology Law (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application discloses calculation transfer demand response system for coordination data center and power grid, comprising: the system comprises a data center A capable of participating in power market demand response, a data center cluster B incapable of participating in power market demand response, a declaration calculation module, a scheduling center, a power trading center, an intelligent pricing module and an execution calculation module. The system coordinates the operation of the power grid and the communication network, provides the constraints and characteristics of the communication network and the power grid of the data center with different nodes, transfers the load to the nodes with distributed power generation or energy storage by using a calculation task transfer method, and realizes the scheme that the data center participates in the power grid demand response. Furthermore, in consideration of the fact that the optimal transfer data volume needs to be obtained by adopting an optimization algorithm when the data center transfers the calculation task, the invention provides an optimization algorithm reconstruction model, which greatly shortens the calculation time while maintaining the accuracy of the optimal solution, thereby playing a role in reducing the energy consumption of the data center.

Description

Calculation transfer demand response system oriented to coordination data center and power grid
Technical Field
The application relates to the technical field of electric power, in particular to a calculation transfer demand response system facing a coordination data center and a power grid.
Background
5G communication is an application in the latest generation mobile networks. It makes possible a new type of network designed to connect almost all people and everything together, including machines, living beings and devices. The 5G wireless technology type mass data processing is accelerated, meanwhile, the reliability and the usability are greatly improved, and the user experience quality is obviously improved. Meanwhile, the higher the 5G communication performance, the faster the data is processed, and the higher the power consumption required for the data. Studies have shown that in 5G communication networks, the energy consumption of one data center is about 3-4 times that of a data center in 4G communication networks. Therefore, as the data age develops, data becomes more and more, and 5G communication may cause the computing tasks required to be processed by the data center to increase exponentially, so that a low-voltage distribution network has a remarkable power load increase phenomenon.
In order to improve the energy efficiency of the 5G data center, researchers focus on energy consumption reduction and electric energy utilization. Research on reducing energy consumption is mainly focused on detailed operation and management inside data centers; in terms of power utilization, there is a major focus on how to utilize the backup battery in demand response without considering the internal operational details of the data center. In fact, the operational details inside the data center are with the potential to participate in demand response, in addition to the backup battery. For example, during a peak or power gap period, the data center under a certain node in the power grid needs to be under-loaded, and the computing task of the data center of the node can be partially transferred to the data center or the base station of other nodes for processing, so as to provide a new direction for the data center of the node to participate in the power market.
In summary, under the conditions of limited computing resources and short power consumption, a solution capable of effectively exploring potential resources is very important, so that the quality of user experience is ensured, and resource allocation and potential resources are optimized and utilized fully and reasonably.
Disclosure of Invention
The application provides a calculation transfer demand response system facing a coordination data center and a power grid, which is used for fully and reasonably optimizing resource allocation and utilizing potential resources.
In view of the above, a first aspect of the present application provides a compute-shift demand response system for coordinating a data center and a power grid, the system including: the system comprises a data center A which can participate in power market demand response, a data center cluster B which cannot participate in power market demand response, a declaration calculation module, a dispatching center, a power trading center, an intelligent pricing module and an execution calculation module;
the declaration calculation module is used for receiving the CPU calculation task and the required electric energy of the data center A, and receiving the electric energy, the distributed generation data and the energy storage battery data required by the CPU calculation task of the data center cluster B; establishing a first objective function and a mathematical model of declaration constraint by taking the lowest power consumption generated by processing a calculation task by the data center A in a declaration time period as a target, reconstructing the mathematical model, and solving and processing the reconstructed model by adopting an improved high-efficiency Lagrange reconstruction multiplier method to obtain a time sequence virtual output boundary of the data center A and the power consumption generated by processing the calculation task transferred by the data center A in the data center cluster B;
the intelligent pricing module is used for receiving the time sequence virtual output boundary, the historical market clearing result and the historical execution result, and meanwhile, determining the declared price of the data center A which is optimal for the data center A to participate in the current power market demand response according to the historical market clearing result, the historical execution result and the time sequence virtual output boundary by adopting a Q learning algorithm, so that the optimal time sequence reported price of the current participation in the power market demand response is output;
the scheduling center is used for providing energy demand information of the operation day, and the energy demand information comprises: time series demand capacity and time series demand price;
the electric power trading center is used for receiving the optimal time sequence reported electricity price, the time sequence demand capacity and the time sequence demand price, calculating and outputting a market clearing result, wherein the market clearing result comprises: winning capacity in time sequence of each virtual power plant;
the execution calculation module is used for receiving the bid winning capacity in the time sequence and receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; and establishing a second objective function and a mathematical model of execution constraint in execution time at the lowest cost of a middle scalar quantity which meets the virtual output of the data center A in the spot market, solving an optimal objective cost function value by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving the execution resource combination and the time sequence output capacity of the data center cluster B in an execution stage.
Optionally, the declaration calculation module specifically includes:
the objective function module is used for receiving the CPU calculation task of the data center A and the required electric energy, establishing the first objective function by taking the lowest electricity consumption generated by the data center A for processing the calculation task in the declaration time period as a target, and inputting the first objective function into the optimization algorithm reconstruction module as a data model;
the declaration constraint module is used for receiving the CPU calculation task and the required electric energy of the data center A, and receiving the electric energy, the distributed generation data and the energy storage battery data required by the CPU calculation task of the data center cluster B; establishing power utilization constraint of the data center cluster B in a declaration time period, charging and discharging quantity upper and lower limit constraints of an energy storage battery of the ith data center in the data center cluster B in unit time, energy storage battery capacity constraint of the ith data center in the data center cluster B, calculation task processing time delay constraint and transmission time delay constraint among data, and inputting the constraints into an optimization algorithm reconstruction module as a data model;
an optimization algorithm reconstruction module, configured to receive the mathematical models of the objective function module and the declaration constraint module, perform mathematical model reconstruction at the same time, reconstruct the mathematical model containing inequality constraints and equality constraints into a mathematical model containing only equality constraints, obtain a reconstruction model, solve the optimal data amount of the data center a migrated to the ith data center in the data center cluster B and the independent variable of the first objective function by using an improved high-efficiency lagrangian reconstruction multiplier method, and solve a time sequence virtual output boundary of the data center a and an electricity consumption generated by a CPU of each data center in the data center cluster B due to processing of the calculation task transferred by the data center a;
and reconstructing inequality constraints in an index constraint function into the first objective function in a form of penalty terms by the improved efficient Lagrange reconstruction multiplier method.
Optionally, the power usage constraint comprises: reporting the electricity consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task within a time period; declaring electric power generated by distributed power generation of the ith data center in the data center cluster B within a time period; and reporting the charge and discharge amount of the energy storage battery of the ith data center in the data center cluster B in the time period.
Optionally, the intelligent pricing module specifically includes:
the initialization module is used for setting an initialization Q table, wherein the initialization Q table is a two-dimensional table with M rows and 2+ N columns, M represents the number of times that the data center A participates in the historical demand response, N represents the action behavior selection of the data center, the first column of the initialization Q table is used for reporting the electricity price for the historical time sequence of the data center A, the second column of the initialization Q table is used for responding to the historical demand participation of the data center A, and the third column to 2+ N column of the initialization Q table are used for setting future reward expectation values obtained by the data center A for different actions made in time;
a behavior decision module for receiving the time sequence virtual output boundary, the historical market clearing result and the historical execution result, selecting an action Q for the intelligent agent, and acquiring a corresponding state according to the initialized Q table, so as to pass through
Figure 996861DEST_PATH_IMAGE001
Determining an optimal time sequence reporting electricity price by a strategy, and outputting the optimal time sequence reporting electricity price, wherein a historical market clearing result is a scalar in a time sequence, and a historical execution result is a time sequence output capacity and execution resource combination;
and the updating module is used for receiving the electricity price reported by the optimal time sequence and updating the initialized Q table by using a Bellman equation.
Optionally, the executing calculation module specifically includes:
the execution constraint module is used for receiving the bid winning capacity in the time sequence and receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; establishing power utilization constraint of the ith data center of the data center cluster B in execution time, distributed generation power constraint of the ith data center of the data center cluster B, upper and lower limit constraint of charging and discharging power of an energy storage battery of the ith data center of the data center cluster B, capacity constraint of the energy storage battery of the ith data center of the data center cluster B and power consumption constraint of a CPU (central processing unit) of the ith data center of the data center cluster B in execution time period, wherein the power consumption constraint is generated by processing calculation tasks and is used as a data model to be input to an optimization algorithm module;
the target cost function module is used for receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; establishing the second objective function with the lowest cost of the medium scalar quantity which meets the virtual output of the data center A in the spot market, and inputting the second objective function serving as a data model into an optimization algorithm module;
the optimization algorithm module is used for receiving the mathematical models of the target cost function module and the execution constraint module, solving an optimal target cost function value by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving an execution resource combination and a time sequence output capacity of the data center cluster B in an execution stage;
and reconstructing inequality constraints in the index constraint function into the second objective function in a form of penalty terms by using the improved efficient Lagrange reconstruction multiplier method.
Optionally, the first objective function of the reconstruction model is:
Figure 477521DEST_PATH_IMAGE002
in the formula (I), the compound is shown in the specification,
Figure 721027DEST_PATH_IMAGE003
the representative reconstruction model aims at minimizing the sum of the power consumption and each penalty term generated by the CPU of the data center A in the declaration time interval due to the processing of the calculation tasks;
Figure 250229DEST_PATH_IMAGE004
an argument representing the first objective function is a total data volume of the reporting period migrated to the data center cluster B by the data center A;
Figure 944515DEST_PATH_IMAGE005
the total data volume is the total data volume when the data center A does not carry out calculation task transfer;
Figure 48607DEST_PATH_IMAGE006
calculating a CPU period required for a unit data volume for the data center A;
Figure 828344DEST_PATH_IMAGE007
calculating the energy consumption of each CPU cycle for the data center A;
Figure 161236DEST_PATH_IMAGE008
Figure 975608DEST_PATH_IMAGE009
penalty coefficients with values larger than 1 are respectively;
Figure 486486DEST_PATH_IMAGE010
an odd number greater than 1 for the amplification factor;
Figure 753520DEST_PATH_IMAGE011
an upper boundary penalty term of a jth inequality constraint of an ith data center in the data center cluster B;
Figure 890103DEST_PATH_IMAGE012
a lower boundary penalty term of a jth inequality constraint of an ith data center in the data center cluster B;
Figure 558982DEST_PATH_IMAGE013
and
Figure 739296DEST_PATH_IMAGE014
respectively an upper boundary penalty item and a lower boundary penalty item of the distributed power generation amount constraint of the ith data center in the data center cluster B in unit time;
Figure 431309DEST_PATH_IMAGE015
and
Figure 168320DEST_PATH_IMAGE016
respectively an upper boundary penalty item and a lower boundary penalty item of the charge and discharge amount constraint of the energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 400629DEST_PATH_IMAGE017
and
Figure 830473DEST_PATH_IMAGE018
respectively limiting upper and lower boundary punishment items for the capacity of the energy storage battery of the ith data center in the data center cluster B;
Figure 9782DEST_PATH_IMAGE019
an upper bound penalty term for declaring a data capacity constraint of the ith data center in the data center cluster B within the time period;
Figure 284905DEST_PATH_IMAGE020
an upper boundary penalty item is constrained for the channel transmission delay between the data center A and the ith data center in the data center cluster B;
Figure 115327DEST_PATH_IMAGE021
the distributed power generation amount of the ith data center in the data center cluster B in unit time is obtained;
Figure 388177DEST_PATH_IMAGE022
the distributed power generation amount is the maximum value of the distributed power generation amount of the ith data center in the data center cluster B in unit time;
Figure 117098DEST_PATH_IMAGE023
to declare the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period,
Figure 884328DEST_PATH_IMAGE024
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is calculated;
Figure 382306DEST_PATH_IMAGE025
the maximum charge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 91636DEST_PATH_IMAGE026
the maximum discharge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 42274DEST_PATH_IMAGE027
the minimum capacity of an energy storage battery of the ith data center in the data center cluster B is obtained;
Figure 111730DEST_PATH_IMAGE028
the initial capacity of an energy storage battery of the ith data center in the data center cluster B before the declaration time period is obtained;
Figure 464214DEST_PATH_IMAGE029
the maximum capacity of an energy storage battery of the ith data center in the data center cluster B is obtained;
Figure 344445DEST_PATH_IMAGE030
reporting the data volume of the ith data center in the data center cluster B which is not subjected to calculation task transfer in the time period, namely the calculation task to be processed by the ith data center;
Figure 467866DEST_PATH_IMAGE031
the data volume of the ith data center migrated to the data center cluster B for the data center A;
Figure 888483DEST_PATH_IMAGE032
reporting the maximum data capacity of the ith data center in the data center cluster B in the time period;
Figure 298736DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 599136DEST_PATH_IMAGE034
computing the maximum transmission time delay of the ith data center in the cluster B of the data centers from the data center A to the computing task in the reporting time period;
Figure 462050DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
Optionally, the equation for reporting electricity usage for data center cluster B over time is constrained by:
Figure 686357DEST_PATH_IMAGE036
in the formula (I), the compound is shown in the specification,
Figure 436270DEST_PATH_IMAGE037
reporting the total power consumption of the ith data center in the data center cluster B to the market in the day ahead for a reporting time period;
Figure 720621DEST_PATH_IMAGE038
reporting the electricity consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the time period;
Figure 70831DEST_PATH_IMAGE039
generating power for distributed generation of an ith data center in the data center cluster B within the declaration period;
Figure 98830DEST_PATH_IMAGE023
and reporting the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period.
Optionally, a variable
Figure 467363DEST_PATH_IMAGE038
Is constrained by:
Figure 594719DEST_PATH_IMAGE040
in the formula (I), the compound is shown in the specification,
Figure 180028DEST_PATH_IMAGE041
calculating the CPU period required by unit data volume for the ith data center in the data center cluster B in the declaration time period;
Figure 746138DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 985490DEST_PATH_IMAGE043
calculating channel consumption power when an ith data center in a data center cluster B and a data center A perform task transmission in unit time;
Figure 533014DEST_PATH_IMAGE044
calculating the time required for transferring the task from the data center A to the ith data center in the data center cluster B within the declaration time period;
Figure 920133DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 227618DEST_PATH_IMAGE045
reporting the channel capacity between the data center A and the ith data center in the data center cluster B in a time period, namely the maximum average information rate which can be transmitted by the channel between the data center A and the ith data center;
Figure 806629DEST_PATH_IMAGE046
transferring the computing task from the data center A to the ith data center in the data center cluster B to obtain the maximum power consumption of the channel;
Figure 603684DEST_PATH_IMAGE004
reporting the total data volume of the data center A migrated to the data center cluster B in the time period;
Figure 150203DEST_PATH_IMAGE035
is a number ofThe number of data centers included in data center cluster B.
Optionally, a variable
Figure 323695DEST_PATH_IMAGE039
The equation constraints are:
Figure 521327DEST_PATH_IMAGE047
in the formula (I), the compound is shown in the specification,
Figure 426966DEST_PATH_IMAGE039
to declare the power generated by distributed generation of the ith data center in data center cluster B over the period of time,
Figure 523098DEST_PATH_IMAGE021
distributed power generation per unit time for the ith data center in the data center cluster B,
Figure 920188DEST_PATH_IMAGE048
to declare a time period.
Optionally, a variable
Figure 723059DEST_PATH_IMAGE023
Is constrained by:
Figure 861916DEST_PATH_IMAGE049
in the formula (I), the compound is shown in the specification,
Figure 632295DEST_PATH_IMAGE023
to declare the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period,
Figure 819694DEST_PATH_IMAGE024
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is set;
Figure 804968DEST_PATH_IMAGE050
and the charging and discharging efficiency of the energy storage battery of the ith data center in the data center cluster B is obtained.
Optionally, the solving of the time-sequence virtual output boundary of the data center a and the power consumption generated by the CPU of each data center in the data center cluster B due to processing the calculation task transferred by the data center a specifically includes:
calculating a time sequence virtual output boundary through a time sequence virtual output boundary calculation formula, and calculating power consumption through a power consumption calculation formula;
the time sequence virtual output boundary calculation formula is as follows:
Figure 803142DEST_PATH_IMAGE051
the formula for calculating the electricity consumption is as follows:
Figure 811549DEST_PATH_IMAGE052
in the formula (I), the compound is shown in the specification,
Figure 51906DEST_PATH_IMAGE053
a timing virtual force boundary;
Figure 891686DEST_PATH_IMAGE054
the electricity consumption is used;
Figure 310029DEST_PATH_IMAGE004
data volume for migrating the computing tasks in the data center A to the data center cluster B;
Figure 287956DEST_PATH_IMAGE006
calculating a CPU period required for a unit data volume for the data center A;
Figure 145054DEST_PATH_IMAGE007
calculating the energy consumption of each CPU cycle for the data center A;
Figure 777023DEST_PATH_IMAGE031
the data volume of the ith data center migrated to the data center cluster B for the data center A;
Figure 428585DEST_PATH_IMAGE041
calculating the CPU period required by unit data volume for the ith data center in the data center cluster B in the declaration time period;
Figure 660852DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 993744DEST_PATH_IMAGE043
calculating channel consumption power when an ith data center in a data center cluster B and a data center A perform task transmission in unit time;
Figure 542537DEST_PATH_IMAGE044
the time required for the transfer of the computing task from data center a to the ith data center in data center cluster B within the declaration period.
Optionally, the second objective function is:
Figure 318994DEST_PATH_IMAGE055
in the formula (I), the compound is shown in the specification,
Figure 523711DEST_PATH_IMAGE056
Figure 722611DEST_PATH_IMAGE057
distributed generating capacity per unit time and charging and discharging capacity of an energy storage battery per unit time of the ith data center in the data center cluster B in the execution time are respectively also independent variables of the second objective function;
Figure 312861DEST_PATH_IMAGE058
Figure 509487DEST_PATH_IMAGE059
respectively the distributed power generation cost and the energy storage battery power generation cost of the ith data center in the data center cluster B;
Figure 949302DEST_PATH_IMAGE060
Figure 686314DEST_PATH_IMAGE061
distributed generation capital recovery coefficients and battery capital recovery coefficients for the ith data center in the data center cluster B, respectively;
Figure 147383DEST_PATH_IMAGE062
Figure 764178DEST_PATH_IMAGE063
respectively representing the distributed power generation and the sticking rate of the energy storage battery of the ith data center in the data center cluster B;
Figure 677907DEST_PATH_IMAGE064
Figure 218610DEST_PATH_IMAGE065
the service life of distributed power generation equipment and energy storage battery equipment of the ith data center in the data center cluster B respectively;
Figure 550496DEST_PATH_IMAGE066
Figure 823346DEST_PATH_IMAGE067
the installation costs of distributed power generation and energy storage batteries of the ith data center in the data center cluster B are respectively set;
Figure 739218DEST_PATH_IMAGE068
Figure 818033DEST_PATH_IMAGE069
respectively, of the ith data center in data center cluster BThe operation and maintenance costs of the distributed power generation and energy storage batteries;
Figure 253693DEST_PATH_IMAGE070
is an execution time period;
Figure 445247DEST_PATH_IMAGE050
charging and discharging efficiency of an energy storage battery of the ith data center in the data center cluster B;
Figure 661464DEST_PATH_IMAGE035
the number of data centers included in data center cluster B;
Figure 481653DEST_PATH_IMAGE071
the aim is that the sum of the distributed power generation cost of the data center cluster B and the power generation cost of the energy storage battery in the execution time period is the lowest.
Optionally, the power usage constraint is:
Figure 21087DEST_PATH_IMAGE072
in the formula (I), the compound is shown in the specification,
Figure 901319DEST_PATH_IMAGE073
total power consumption declared to the day-ahead market for the ith data center in the data center cluster B in the execution period;
Figure 27669DEST_PATH_IMAGE074
the power consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the execution period;
Figure 572920DEST_PATH_IMAGE075
electric power generated for distributed power generation of the ith data center in the data center cluster B within the execution period;
Figure 45489DEST_PATH_IMAGE076
for the ith data in the data center cluster B in the execution periodThe charge and discharge capacity of the core energy storage cell.
Optionally, the distributed generation power constraint is:
Figure 831043DEST_PATH_IMAGE077
in the formula (I), the compound is shown in the specification,
Figure 756273DEST_PATH_IMAGE056
the distributed power generation amount of the ith data center in the data center cluster B in the execution stage in unit time is calculated;
Figure 666067DEST_PATH_IMAGE022
the maximum value of the distributed power generation amount of the ith data center in the data center cluster B in unit time.
Optionally, the charging and discharging power upper and lower limits of the energy storage battery are constrained as follows:
Figure 930826DEST_PATH_IMAGE078
in the formula (I), the compound is shown in the specification,
Figure 215177DEST_PATH_IMAGE057
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in the execution stage in unit time is calculated;
Figure 814655DEST_PATH_IMAGE025
the maximum charging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 514757DEST_PATH_IMAGE026
the maximum discharge amount of the energy storage battery of the ith data center in the data center cluster B in unit time.
Optionally, the energy storage battery capacity constraint is:
Figure 696340DEST_PATH_IMAGE079
in the formula (I), the compound is shown in the specification,
Figure 840008DEST_PATH_IMAGE027
Figure 677514DEST_PATH_IMAGE029
respectively obtaining the minimum capacity and the maximum capacity of an energy storage battery of the ith data center in the data center cluster B;
Figure 430575DEST_PATH_IMAGE080
the initial capacity of an energy storage battery of the ith data center in the data center cluster B before the execution stage is obtained;
Figure 466664DEST_PATH_IMAGE076
the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the execution period is shown.
Optionally, the power usage constraint is:
Figure 30501DEST_PATH_IMAGE081
in the formula (I), the compound is shown in the specification,
Figure 417620DEST_PATH_IMAGE082
the data volume of the ith data center in the data center cluster B, which is not subjected to the calculation task transfer, is obtained;
Figure 43265DEST_PATH_IMAGE083
the data volume of the ith data center in the data center cluster B migrated to the data center A in the execution stage;
Figure 137123DEST_PATH_IMAGE041
calculating a CPU period required by unit data volume for the ith data center in the data center cluster B;
Figure 934178DEST_PATH_IMAGE042
ith data in data center cluster BThe center calculates the energy consumption of each CPU period;
Figure 729964DEST_PATH_IMAGE084
calculating the channel consumption power in unit time when the ith data center in the data center cluster B and the data center A in the execution stage perform calculation task transmission;
Figure 903457DEST_PATH_IMAGE085
the time required for transferring the calculation task from the data center A to the ith data center in the data center cluster B in the execution stage is calculated;
Figure 851821DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 508193DEST_PATH_IMAGE086
the channel capacity between the data center A and the ith data center in the data center cluster B in the execution stage, namely the maximum average information rate which can be transmitted by the channel between the data center A and the ith data center in the data center cluster B;
Figure 604325DEST_PATH_IMAGE046
transferring the computing task from the data center A to the ith data center in the data center cluster B to obtain the maximum power consumption of the channel;
Figure 988033DEST_PATH_IMAGE087
to perform the total amount of data that data center a migrates to data center cluster B during the phase,
Figure 305750DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
Optionally, the solving of the execution resource combination and the time-series capacity of the data center cluster B in the execution stage specifically includes:
substituting the execution resource combination into a time sequence output capacity calculation formula to calculate and obtain the sum time sequence output capacity of the data center cluster B in the execution stage;
wherein, the time sequence output capacity calculation formula is as follows:
Figure 444608DEST_PATH_IMAGE088
in the formula (I), the compound is shown in the specification,
Figure 965719DEST_PATH_IMAGE089
time-series capacity for the data center cluster B in the execution phase;
Figure 900920DEST_PATH_IMAGE035
the number of data centers included in data center cluster B;
Figure 886194DEST_PATH_IMAGE056
Figure 133636DEST_PATH_IMAGE057
the distributed power generation amount per unit time and the charging and discharging amount per unit time of the energy storage battery of the ith data center in the data center cluster B in the execution time are also independent variables of the second objective function, namely the execution resource combination.
According to the technical scheme, the method has the following advantages:
the application provides a calculation transfer demand response system facing a coordination data center and a power grid, which not only considers communication networks among data centers of different nodes, but also combines processing calculation tasks of the data centers with operation of the power grid. The system coordinates the operation of the power grid and the communication network simultaneously, provides the constraints and characteristics of the communication network and the power grid of the data center with different nodes, transfers the load to the nodes with distributed power generation or energy storage by utilizing a calculation task transfer method, realizes the participation of the data center in a power grid demand response scheme, and plays a positive role in peak clipping and valley filling for the power system in the peak period of power utilization.
In consideration of the fact that an optimization algorithm is needed to obtain the optimal transfer data volume when a data center transfers a calculation task, although the traditional Lagrange multiplier method based on the KKT condition can solve the optimal transfer data volume, the traditional Lagrange multiplier method based on the KKT condition has more constraint conditions, so that the calculation time is longer, and the energy consumption is higher when a CPU of the data center performs the solution.
According to the invention, through the characteristics of a communication network and data transferability, the data center which can participate in market demand response transfers calculation tasks to the data center which has different nodes and a small power supply, thereby playing a guiding role in participating in the demand response market by a communication data type heavy load aggregator, playing a role in peak clipping and valley filling, reducing the high load operation cost of a power system, improving the social energy efficiency level, and having a substantial guiding role in promoting the communication data amount type load aggregator to participate in the demand response scheme of the power market. Under the condition of power supply shortage and power gap, or when the transformer, the line, the feeder line and the like exist in a local area and have heavy overload risks, the pressure of the power system is relieved, potential resources are explored, optimal configuration of response resources is realized, and a positive effect is provided for stable and reliable operation of the power system.
Drawings
Fig. 1 is a schematic structural diagram of an embodiment of a computation-based transfer demand response system for a coordination data center and a power grid provided in an embodiment of the present application;
FIG. 2 is a flow diagram of a declaration calculation module provided in an embodiment of the present application;
FIG. 3 is a flow diagram of an intelligent pricing module provided in an embodiment of the present application;
fig. 4 is a flowchart of a module for performing computations provided in an embodiment of the present application.
Detailed Description
In order to make those skilled in the art better understand the technical solutions of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Referring to fig. 1, in an embodiment of the present application, a computing transfer demand response system for coordinating a data center and a power grid includes: the system comprises a data center A capable of participating in power market demand response, a data center cluster B incapable of participating in power market demand response, a declaration calculation module, a scheduling center, a power trading center, an intelligent pricing module and an execution calculation module;
it should be noted that the data center a and the data center cluster B are located in different nodes of the power grid, but belong to the same data manager, and the data manager can schedule calculation tasks among different data centers according to needs, thereby indirectly transferring power consumption among the data centers.
The data center A is a load aggregator or a large electricity consumption user in an electricity market, consumed electric energy is from a power grid, an energy storage or small power supply is not arranged, and electric energy required by processing of a CPU (central processing unit) of the data center on a calculation task is mainly considered in electric load consumption. In an early declaration stage, the data center A has direct control capability and other conditions which are respectively aggregated to form a virtual power plant, the virtual power plant can be used as a demand side response of the virtual power plant participating in the power market, and a data manager can transfer a calculation task of the data center A to a data center cluster B with an energy storage or small power supply through a declaration calculation module. In the later execution calculation stage, the data center a is used for receiving a time sequence bid amount, which is a market clearing result issued by the power trading center, and transmitting the time sequence bid amount to the data center cluster B.
The data center cluster B is a load aggregation business cluster or a power consumption large user cluster of different nodes in the power grid, and the data center cluster B is connected with distributed power generation and energy storage, so that the power source of the data center cluster B comprises the power grid, self distributed power generation and self energy storage batteries. The data center cluster B is only used as a load in the power market, the power consumption purchased from the power grid is fixed as the power consumption declared in the market at the day before, and the power center cluster B does not serve as a demand response of the virtual power plant participating in the power market. In an early declaration stage, the data center cluster B is used as input data of a declaration calculation module, and the data volume of the data center cluster B which can maximally accommodate the data center A is obtained through an optimization algorithm. In the later execution calculation stage, the data center cluster B is used for receiving the time sequence bid-winning capacity obtained by the data center A from the power trading center and transmitting the time sequence bid-winning capacity to the execution calculation module, so that the optimal execution scheme of the data center cluster B is obtained.
The reporting and calculating module is used for receiving the CPU calculating task and the required electric energy of the data center A, and receiving the electric energy, the distributed generation data and the energy storage battery data required by the CPU calculating task of the data center cluster B; establishing a first objective function and a mathematical model of declaration constraint by taking the lowest power consumption generated by processing a calculation task by a data center A in a declaration time period as a target, reconstructing the mathematical model, and solving and processing the reconstructed model by adopting an improved high-efficiency Lagrange reconstruction multiplier method to obtain a time sequence virtual output boundary of the data center A and the power consumption generated by processing the calculation task transferred by the data center A in a data center cluster B;
as shown in fig. 2, in this embodiment, the declaration calculation module specifically includes an objective function module, a declaration constraint module, and an optimization algorithm reconstruction module.
Wherein:
the objective function module is used for receiving a CPU calculation task of the data center A and required electric energy, establishing an objective function by taking the lowest electricity consumption generated by the data center A for processing the calculation task in a reporting time interval (a time interval when the data center A participates in demand response) as a target, and inputting the objective function into the optimization algorithm reconstruction module as a data model;
it should be noted that the objective function is as follows:
Figure 391311DEST_PATH_IMAGE090
(1)
Figure 444717DEST_PATH_IMAGE003
representing the occurrence of computing tasks being processed by the CPU of data centre A within the declaration periodThe lowest electricity consumption is taken as a target;
Figure 222180DEST_PATH_IMAGE004
representing the objective function independent variable, and reporting the total data volume (bit) of the data center A transferred to the data center cluster B in the time period;
Figure 391256DEST_PATH_IMAGE005
the total data volume (bit) when the data center A does not carry out calculation task transfer is obtained;
Figure 683697DEST_PATH_IMAGE006
a CPU cycle required to calculate a unit data amount (1 bit) for the data center a;
Figure 478477DEST_PATH_IMAGE007
the energy consumption per CPU cycle is calculated for data center a.
The reporting constraint module is used for receiving the CPU calculation task and the required electric energy of the data center A, and receiving the electric energy, the distributed generation data and the energy storage battery data required by the CPU calculation task of the data center cluster B; establishing power utilization constraint of the data center cluster B in a declaration time period, charge and discharge quantity upper and lower limit constraint of an energy storage battery of the ith data center in the data center cluster B in unit time, energy storage battery capacity constraint of the ith data center in the data center cluster B, calculation task processing time delay constraint and transmission time delay constraint among data, and inputting the data serving as a data model into an optimization algorithm reconstruction module;
it should be noted that the declaration constraint module is configured to receive the CPU computation task and the required power from the data center a, the energy consumption, the distributed power generation data, the data related to the energy storage battery required by the CPU computation task from the data center cluster B, and the channel data between the data center a and each data center cluster B. For the data center cluster B, corresponding electric quantity is purchased to the power grid in the market at the day before in a contract form, and certain punishment is carried out if the quantity of the contract is violated. Therefore, the data manager reports the virtual output boundary to the electric power trading center in consideration of the data center AIn time, the power consumption purchased from each node of the data center cluster B to the power grid needs to be considered and fixed as the total power consumption declared in the market at present
Figure 172764DEST_PATH_IMAGE037
. The power utilization constraint of the data center cluster B in the reporting time is shown as an equation (2):
Figure 11276DEST_PATH_IMAGE091
(2)
Figure 994275DEST_PATH_IMAGE037
reporting the total power consumption of the ith data center in the data center cluster B to the day-ahead market in the reporting time period;
Figure 389485DEST_PATH_IMAGE038
reporting the electricity consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the time period;
Figure 623764DEST_PATH_IMAGE039
generating power for distributed generation of an ith data center in the data center cluster B within the declaration period;
Figure 649488DEST_PATH_IMAGE023
and reporting the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period. Wherein
Figure 916522DEST_PATH_IMAGE038
The formula (3) is shown in the formula,
Figure 36793DEST_PATH_IMAGE039
the formula (4) is shown in the formula,
Figure 643355DEST_PATH_IMAGE023
the formula (6) is shown in the formula;
Figure 590714DEST_PATH_IMAGE092
(3)
Figure 17147DEST_PATH_IMAGE030
reporting the data volume (namely the calculation task to be processed) of the ith data center in the data center cluster B which is not subjected to calculation task transfer in the time period;
Figure 206689DEST_PATH_IMAGE031
a data volume (bit) for migrating the data center A to the ith data center in the data center cluster B;
Figure 667757DEST_PATH_IMAGE041
calculating a CPU period required by a unit data volume (1 bit) for the ith data center in the data center cluster B in the declaration time period;
Figure 832022DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 759133DEST_PATH_IMAGE043
calculating channel power consumption when the ith data center in the data center cluster B and the data center A perform calculation task transmission in unit time;
Figure 299836DEST_PATH_IMAGE044
calculating the time required for transferring the task from the data center A to the ith data center in the data center cluster B within the declaration time period;
Figure 615411DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 216157DEST_PATH_IMAGE045
for the data center A and the data center cluster B in the declaration periodi, the channel capacity among the data centers, namely the maximum average information rate which can be transmitted by the channels among the data centers;
Figure 866450DEST_PATH_IMAGE046
transferring the computing task from the data center A to the channel maximum power consumption of the ith data center in the data center cluster B;
Figure 148526DEST_PATH_IMAGE004
total data volume (bit) for migrating the data center A to the data center cluster B in the declaration time period;
Figure 646504DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
Figure 840987DEST_PATH_IMAGE093
(4)
Figure 57205DEST_PATH_IMAGE094
(5)
Formula (4) represents the power generated by the distributed power generation of the ith data center in the data center cluster B within the declaration time period, and formula (5) is the upper and lower limit constraints of the distributed power generation of the ith data center in the data center cluster B within unit time. Wherein
Figure 877393DEST_PATH_IMAGE039
To declare the power generated by distributed generation of the ith data center in data center cluster B over the period of time,
Figure 416828DEST_PATH_IMAGE021
the distributed power generation amount of the ith data center in the data center cluster B in unit time,
Figure 93797DEST_PATH_IMAGE048
in order to declare a time period,
Figure 469415DEST_PATH_IMAGE022
the maximum value of the distributed power generation amount of the ith data center in the data center cluster B in unit time.
Figure 155611DEST_PATH_IMAGE095
(6)
Figure 48087DEST_PATH_IMAGE096
(7)
Figure 99220DEST_PATH_IMAGE097
(8)
Formula (6) represents the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the declaration time period; formula (7) is the charge and discharge amount upper and lower limit constraints of the energy storage battery of the ith data center in the data center cluster B in unit time; equation (8) is the energy storage battery capacity constraint for the ith data center in data center cluster B.
Wherein
Figure 24450DEST_PATH_IMAGE023
To declare the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period,
Figure 435709DEST_PATH_IMAGE024
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is calculated;
Figure 762785DEST_PATH_IMAGE050
for the charging and discharging efficiency of the energy storage battery of the ith data center in the data center cluster B, the energy storage battery has an electric energy conversion process, and the mutual conversion between the electric energy and the chemical energy can cause the electric energy loss, so that the coefficient is obtained;
Figure 984819DEST_PATH_IMAGE048
to declare a time period;
Figure 820182DEST_PATH_IMAGE025
the maximum charge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 848181DEST_PATH_IMAGE026
the maximum discharge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 967447DEST_PATH_IMAGE027
the minimum capacity of an energy storage battery of the ith data center in the data center cluster B is obtained;
Figure 609650DEST_PATH_IMAGE028
the initial capacity of the energy storage battery of the ith data center in the data center cluster B before the reporting time period is reported;
Figure 509472DEST_PATH_IMAGE029
and the maximum capacity of the energy storage battery of the ith data center in the data center cluster B.
The data transmission between the ith data center in the data center cluster B and the data center A needs to meet certain service quality, including the calculation task processing delay constraint and the transmission delay constraint between data. See formula (9) and formula (10):
Figure 13266DEST_PATH_IMAGE098
(9)
Figure 49355DEST_PATH_IMAGE099
(10)
equation (9) is a data capacity constraint of the ith data center in the data center cluster B in the declaration time period, and represents that the calculation task undertaken in the declaration time period does not exceed the upper limit of the calculation capacity, and once the upper limit of the capacity is exceeded, more time is needed for the calculationProcessing the computational tasks, resulting in a degradation of quality of service; equation (10) is the channel transmission delay constraint between data center a and the ith data center in data center cluster B, which is the maximum transmission delay that can be received by the user. Wherein the content of the first and second substances,
Figure 360994DEST_PATH_IMAGE030
reporting the data volume (namely the calculation task to be processed) of the ith data center in the data center cluster B which is not subjected to calculation task transfer in the time period;
Figure 482534DEST_PATH_IMAGE031
a data volume (bit) for migrating the data center A to the ith data center in the data center cluster B;
Figure 790019DEST_PATH_IMAGE032
reporting the maximum data capacity of the ith data center in the data center cluster B in the time period;
Figure 133144DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 930199DEST_PATH_IMAGE034
and reporting the maximum transmission time delay of the calculation task transferred from the data center A to the ith data center in the data center cluster B in the time period.
The optimization algorithm reconstruction module is used for receiving the mathematical models of the target function module and the declaration constraint module, reconstructing the mathematical models containing inequality constraints and equality constraints into the mathematical models only containing equality constraints to obtain reconstructed models, solving the optimal data volume of the data center A transferred to the ith data center in the data center cluster B and the target function independent variable by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving the time sequence virtual output boundary of the data center A and the power consumption generated by the CPU of each data center in the data center cluster B due to the processing of the calculation task transferred by the data center A;
and reconstructing inequality constraints in the index constraint function into the target function in a form of penalty terms by using the improved efficient Lagrange reconstruction multiplier method.
It should be noted that the optimization algorithm reconstruction module is used for receiving the mathematical models of the objective function module and the declaration constraint module, and solving the optimal mathematical model by adopting the improved high-efficiency Lagrange reconstruction multiplier method
Figure 476718DEST_PATH_IMAGE031
And
Figure 384631DEST_PATH_IMAGE004
thus, the time-series virtual output boundary is further calculated and the value is output. The improved high-efficiency Lagrange reconstruction multiplier method reconstructs inequality constraints in the index constraint function into the target function in a form of penalty terms, and the calculation speed is higher than that of a typical Lagrange multiplier method based on a KKT condition and has a similar global optimal value. The first objective function of its reconstructed mathematical model (reconstruction model) is as follows:
Figure 83728DEST_PATH_IMAGE002
(11)
Figure 989367DEST_PATH_IMAGE100
(12)
Figure 85499DEST_PATH_IMAGE101
(13)
Figure 984054DEST_PATH_IMAGE102
(14)
Figure 114821DEST_PATH_IMAGE103
(15)
equation (11) is a first objective function of the reconstructed model; formula (12), formula (13)Equation (14) and equation (15) are the electricity consumption equation constraint and the variable of the data center cluster B in the declaration time respectively
Figure 191361DEST_PATH_IMAGE038
Is constrained by an equality of (1), a variable
Figure 509210DEST_PATH_IMAGE039
Equality constraint, variable
Figure 709991DEST_PATH_IMAGE023
Is constrained by the equation of (a).
Figure 632948DEST_PATH_IMAGE003
The representative reconstruction model aims at minimizing the sum of the power consumption and each penalty term generated by the CPU of the data center A in the declaration time interval due to the processing of the calculation tasks;
Figure 129657DEST_PATH_IMAGE004
an argument representing the first objective function is a total data volume (bit) of the reporting period migrated to the data center cluster B by the data center a;
Figure 934802DEST_PATH_IMAGE005
the total data volume bit when the data center A does not carry out calculation task transfer is obtained;
Figure 925892DEST_PATH_IMAGE006
calculating a CPU period required by 1bit of unit data volume for the data center A;
Figure 765672DEST_PATH_IMAGE007
calculating the energy consumption of each CPU cycle for the data center A;
Figure 934747DEST_PATH_IMAGE008
Figure 430450DEST_PATH_IMAGE009
penalty coefficients with values greater than 1, respectively;
Figure 21969DEST_PATH_IMAGE010
an odd number greater than 1 for the amplification factor;
Figure 168785DEST_PATH_IMAGE011
an upper boundary penalty term of a jth inequality constraint of an ith data center in the data center cluster B;
Figure 820346DEST_PATH_IMAGE012
a lower boundary penalty term of a jth inequality constraint of an ith data center in the data center cluster B;
Figure 537767DEST_PATH_IMAGE013
and
Figure 932976DEST_PATH_IMAGE014
respectively an upper boundary penalty item and a lower boundary penalty item of the distributed power generation amount constraint of the ith data center in the data center cluster B in unit time;
Figure 432834DEST_PATH_IMAGE015
and
Figure 192980DEST_PATH_IMAGE016
respectively an upper boundary penalty item and a lower boundary penalty item of the charge and discharge amount constraint of the energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 460013DEST_PATH_IMAGE017
and
Figure 845864DEST_PATH_IMAGE018
respectively limiting upper and lower boundary punishment items for the capacity of the energy storage battery of the ith data center in the data center cluster B;
Figure 514743DEST_PATH_IMAGE019
an upper bound penalty term for declaring a data capacity constraint of the ith data center in the data center cluster B within the time period;
Figure 445789DEST_PATH_IMAGE020
an upper boundary penalty item is constrained for the channel transmission delay between the data center A and the ith data center in the data center cluster B;
Figure 200119DEST_PATH_IMAGE037
reporting the total power consumption of the ith data center in the data center cluster B to the market in the day ahead for a reporting time period;
Figure 891125DEST_PATH_IMAGE038
reporting the electricity consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the time period;
Figure 352194DEST_PATH_IMAGE039
generating power for distributed generation of the ith data center in the data center cluster B within the declaration period;
Figure 516459DEST_PATH_IMAGE023
and reporting the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period.
Figure 945035DEST_PATH_IMAGE030
Reporting the data volume (namely the calculation task to be processed) of the ith data center in the data center cluster B which is not subjected to calculation task transfer in the time period;
Figure 423421DEST_PATH_IMAGE031
a data volume (bit) for migrating the data center A to the ith data center in the data center cluster B;
Figure 801313DEST_PATH_IMAGE041
calculating a CPU period required by a unit data volume (1 bit) for the ith data center in the data center cluster B in the declaration time period;
Figure 110982DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 574324DEST_PATH_IMAGE043
calculating channel consumption power when an ith data center in a data center cluster B and a data center A perform task transmission in unit time;
Figure 856401DEST_PATH_IMAGE044
calculating the time required for transferring the task from the data center A to the ith data center in the data center cluster B within the declaration time period;
Figure 88799DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 47397DEST_PATH_IMAGE045
reporting the channel capacity between the data center A and the ith data center in the data center cluster B in a time period, namely the maximum average information rate which can be transmitted by the channel between the data center A and the ith data center;
Figure 935718DEST_PATH_IMAGE046
transferring the computing task from the data center A to the channel maximum power consumption of the ith data center in the data center cluster B;
Figure 818223DEST_PATH_IMAGE004
total data volume (bit) for migrating the data center A to the data center cluster B in the declaration period;
Figure 859123DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
Wherein the content of the first and second substances,
Figure 801671DEST_PATH_IMAGE039
for reporting in the data within the time periodThe distributed generation of electricity by the ith data center in the core cluster B,
Figure 177289DEST_PATH_IMAGE021
distributed power generation per unit time for the ith data center in the data center cluster B,
Figure 863485DEST_PATH_IMAGE048
in order to declare a time period,
Figure 257426DEST_PATH_IMAGE022
the maximum value of the distributed power generation amount per unit time of the ith data center in the data center cluster B.
Wherein the content of the first and second substances,
Figure 308559DEST_PATH_IMAGE023
to declare the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period,
Figure 233790DEST_PATH_IMAGE024
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is calculated;
Figure 143583DEST_PATH_IMAGE050
for the charging and discharging efficiency of the energy storage battery of the ith data center in the data center cluster B, the energy storage battery has an electric energy conversion process, and the mutual conversion between electric energy and chemical energy can cause electric energy loss, so that the coefficient is obtained;
Figure 408343DEST_PATH_IMAGE025
the maximum charge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 427114DEST_PATH_IMAGE026
the maximum discharge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 26592DEST_PATH_IMAGE027
the minimum capacity of an energy storage battery of the ith data center in the data center cluster B is set;
Figure 54590DEST_PATH_IMAGE028
the initial capacity of an energy storage battery of the ith data center in the data center cluster B before the declaration time period is obtained;
Figure 173856DEST_PATH_IMAGE029
and the maximum capacity of the energy storage battery of the ith data center in the data center cluster B.
Figure 317524DEST_PATH_IMAGE032
Reporting the maximum data capacity of the ith data center in the data center cluster B in the time period;
Figure 951768DEST_PATH_IMAGE034
and reporting the maximum transmission time delay of the calculation task transferred from the data center A to the ith data center in the data center cluster B in the time period.
Therefore, the optimization algorithm reconstruction module reconstructs the mathematical model containing inequality constraint and equality constraint into the mathematical model only containing equality constraint, and the reconstruction model can be solved by adopting a Lagrange multiplier method to solve the optimal mathematical model
Figure 455561DEST_PATH_IMAGE031
And
Figure 757230DEST_PATH_IMAGE004
then, the time-series virtual output boundary of the data center a and the power consumption of the CPU of each data center in the data center cluster B due to processing the calculation task transferred by the data center a can be obtained by equations (16) and (17).
Figure 570334DEST_PATH_IMAGE104
(16)
Figure 629556DEST_PATH_IMAGE105
(17)
Wherein, the formula (16) is a time sequence virtual output boundary of the data center A
Figure 999358DEST_PATH_IMAGE053
Represents; equation (17) is the power consumption of the CPU in the ith data center in the data center cluster B for processing the calculation task transferred by the data center A
Figure 841019DEST_PATH_IMAGE054
And (4) showing.
Figure 575756DEST_PATH_IMAGE004
Data volume (bit) for migration of computing tasks in data center A to data center cluster B;
Figure 184592DEST_PATH_IMAGE006
calculating a CPU period required by a unit data volume (1 bit) for the data center A;
Figure 279456DEST_PATH_IMAGE007
the energy consumption per CPU cycle is calculated for data center a.
Figure 290137DEST_PATH_IMAGE031
A data volume (bit) for migrating the data center A to the ith data center in the data center cluster B;
Figure 195777DEST_PATH_IMAGE041
calculating a CPU period required by a unit data volume (1 bit) for the ith data center in the data center cluster B in the declaration time period;
Figure 980324DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 691928DEST_PATH_IMAGE043
calculating channel consumption power when an ith data center in a data center cluster B and a data center A perform task transmission in unit time;
Figure 494799DEST_PATH_IMAGE044
the time required for the transfer of the computing task from data center a to the ith data center in data center cluster B within the declaration period.
The intelligent pricing module is used for receiving the time sequence virtual output boundary, the historical market clearing result and the historical execution result, and meanwhile, determining the declared price of the data center A which is the best in response to the current power market demand by adopting a Q learning algorithm according to the historical market clearing result, the historical execution result and the time sequence virtual output boundary, so that the time sequence reported price of the current power market demand response is output;
it should be noted that the machine learning method adopted by the intelligent pricing module is a Q learning algorithm, and an optimal action strategy is found by using a Q function, and the maximum reported price of electricity that can be profitable in the future is obtained. Q is the sum of Q (s,
Figure 820607DEST_PATH_IMAGE106
) I.e. in s-state, take action
Figure 404035DEST_PATH_IMAGE106
The expected value can be obtained, and the environment feeds back corresponding rewards according to the actions of the intelligent agent. The main idea of the Q learning algorithm is therefore to combine the state(s) and the action(s) ((s))
Figure 857013DEST_PATH_IMAGE106
) Forming a Q table to store Q values that will represent the maximum future reward value to be achieved by taking the corresponding action given the best strategy
Figure 262193DEST_PATH_IMAGE001
And the strategy selects the behavior action with the maximum Q value in the Q table.
As shown in fig. 3, the intelligent pricing module of this embodiment specifically includes an initialization module, a behavior decision module, and an update module.
Wherein:
the initialization module is used for setting an initialization Q table, the initialization Q table is a two-dimensional table with M rows and 2+ N columns, wherein M represents the number of times that the data center A historically participates in demand response, N represents action behavior selection of the data center, the first column of the initialization Q table is used for reporting the electricity price of the historical time sequence of the data center A, the second column of the initialization Q table is used for responding to the demand history of the data center A, and the third column to the 2+ N column of the initialization Q table are future reward expected values obtained by the data center A for different actions made within the time;
it should be noted that the initialization module is an initialization Q table, which is a two-dimensional table with M rows and 2+ N columns, where M represents the number of times that the data center a has historically participated in the demand response, N represents the selection of the action behavior of the data center, the first column of the Q table is the report price of electricity for the historical time sequence of the data center a, the second column of the Q table is the time that the data center a has historically participated in the demand response, and the third to 2+ N columns of the Q table are the expected future rewards obtained by the data center a for different actions performed within the time, and the value is 0 in the initial state (the Q table needs to be initialized only when the data center a first participates in the demand response). The corresponding relation of historical data corresponding to the time sequence reported electricity price and participation time of the data center A in the Q table determines N action states, and corresponding actions taken in each action state have corresponding reward values, namely reward values R fed back by the environment to the action. The calculation formula of the reward value R is shown as formula (18):
Figure 571952DEST_PATH_IMAGE107
(18)
wherein the content of the first and second substances,
Figure 580359DEST_PATH_IMAGE108
reporting the electricity price to the time sequence declared by the data center A to the power market;
Figure 633766DEST_PATH_IMAGE087
a total data volume (bit) for migrating the data center A to the data center cluster B in the execution phase;
Figure 660497DEST_PATH_IMAGE006
calculating a CPU period required by a unit data volume (1 bit) for the data center A;
Figure 78840DEST_PATH_IMAGE007
calculating the energy consumption of each CPU cycle for the data center A;
Figure 636860DEST_PATH_IMAGE035
the number of data centers included in data center cluster B;
Figure 916794DEST_PATH_IMAGE058
Figure 611080DEST_PATH_IMAGE059
the distributed power generation cost and the energy storage battery power generation cost of the ith data center in the data center cluster B are respectively.
A behavior decision module for receiving the time sequence virtual output boundary, the historical market clearing result and the historical execution result, selecting an action Q for the intelligent agent, and acquiring a corresponding state according to the initialized Q table, thereby passing through
Figure 465904DEST_PATH_IMAGE001
Determining an optimal time sequence reporting electricity price by a strategy, and outputting the time sequence reporting electricity price, wherein a historical market clearing result is a scalar in the time sequence, and a historical execution result is a time sequence output capacity and execution resource combination;
it should be noted that the behavior decision module is configured to receive historical reported electricity prices, time sequence virtual output boundaries, and historical market clearing results (time sequence bid amount, time sequence output capacity, and execution resource combination) of the data center a; by passing
Figure 245641DEST_PATH_IMAGE001
Policy determination of optimalityReporting the electricity price by the time sequence and outputting the reported electricity price by the time sequence. The behavior decision module will select an action for the agent
Figure 827801DEST_PATH_IMAGE106
This action corresponds to the S state in the currently existing Q table. Decision module adoption
Figure 579856DEST_PATH_IMAGE001
The strategy, ε, is the agent's exploration rate, and is typically set to a value of 0 to 1.
Figure 402319DEST_PATH_IMAGE001
The strategy will cause the model to generate a random number and compare this random number to epsilon. If the random number is larger than the epsilon value, the action of the intelligent agent is random exploration; if the random number is less than or equal to epsilon, the action of the agent selects the action with the maximum Q value in the corresponding state in the existing Q table. The action value is the time sequence reported electricity price of the data center A. The Q value is calculated by the formula (19):
Figure 354838DEST_PATH_IMAGE109
(19)
wherein the content of the first and second substances,
Figure 553738DEST_PATH_IMAGE110
is at the same timesBehavior of actions in State
Figure 160300DEST_PATH_IMAGE111
The new value of the Q value is obtained,
Figure 153664DEST_PATH_IMAGE112
for the current value of Q, the value of Q,
Figure 94944DEST_PATH_IMAGE113
is in the next state
Figure 769639DEST_PATH_IMAGE114
The maximum Q value obtained for all of the next actions,
Figure 293024DEST_PATH_IMAGE115
for the forward return to current value,
Figure 411284DEST_PATH_IMAGE116
in order for the agent to feed back the learning rate,
Figure 590592DEST_PATH_IMAGE117
the ratio of reward attenuation of the agent to the long term, and R is the reward value of the environment to the behavior feedback.
And the updating module is used for receiving the time sequence reported electricity price and updating the initialized Q table by using a Bellman equation.
It should be noted that the update module is configured to receive the time sequence reported electricity price output by the behavior decision module. For actions selected in a behavior decision module
Figure 865716DEST_PATH_IMAGE106
Executing the action returns a new state
Figure 696138DEST_PATH_IMAGE114
And corresponding feedback values. The new Q table is updated using the Bellman equation, see equation (18).
Figure 31304DEST_PATH_IMAGE118
(20)
Wherein the content of the first and second substances,
Figure 697909DEST_PATH_IMAGE119
selecting actions for the s state
Figure 776723DEST_PATH_IMAGE106
The new value of Q is obtained in the following way,
Figure 960186DEST_PATH_IMAGE112
is operated in s state
Figure 669516DEST_PATH_IMAGE106
Current Q value of,
Figure 620155DEST_PATH_IMAGE120
The reward value for the context to the behavioral feedback,
Figure 689611DEST_PATH_IMAGE121
is the next new state
Figure 42095DEST_PATH_IMAGE114
All possible actions
Figure 922326DEST_PATH_IMAGE122
The maximum Q value of (2).
Figure 360261DEST_PATH_IMAGE116
The body learning rate is fed back intelligently;
Figure 469293DEST_PATH_IMAGE117
the larger the value, the more important the benefit of the agent to the distant period is, for the rate of decay of the agent to the distant period. The new Q value will be saved in the Q table to provide learning data for the next action of the behavior decision module.
The dispatching center is used for providing energy demand information of the operation days, and the energy demand information comprises: time series demand capacity and time series demand price;
the electric power trading center is used for receiving the time sequence reported electricity price, the time sequence required capacity and the time sequence required price, calculating and outputting a market clearing result, wherein the market clearing result comprises the following steps: winning capacity in time sequence of each virtual power plant;
the execution calculation module is used for receiving the bid winning capacity in the time sequence and receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; and establishing a second objective function and a mathematical model of execution constraint in execution time with the lowest cost of a middle scalar quantity which meets the virtual output of the data center A in the spot market, solving an optimal objective cost function value by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving the execution resource combination and the time sequence output capacity of the data center cluster B in the execution stage.
As shown in fig. 4, in this embodiment, the execution calculation module specifically includes an execution constraint module, an objective cost function module, and an optimization algorithm module.
Wherein:
the target cost function module is used for receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; establishing a second objective function with the lowest cost of the medium scalar quantity which meets the virtual output of the data center A in the spot market, and inputting the second objective function serving as a data model into the optimization algorithm module;
it should be noted that the objective cost function module is configured to receive energy consumption, distributed power generation data, and energy storage battery related data required by the CPU computation task of the data center cluster B. Computing tasks by transfer to satisfy the medium scalar for data center A's virtual contribution in the spot market
Figure 879546DEST_PATH_IMAGE087
And in the data cluster B, adjusting distributed power generation and battery energy storage of the data center cluster B to realize the optimal execution scheme with the lowest bid amount cost, wherein a second objective function is as follows:
Figure 992996DEST_PATH_IMAGE123
(21)
in the formula (I), the compound is shown in the specification,
Figure 105177DEST_PATH_IMAGE056
Figure 329485DEST_PATH_IMAGE057
distributed generating capacity in unit time and charging and discharging capacity of an energy storage battery in unit time of the ith data center in the data center cluster B in the execution time are respectively also independent variables of the second objective function;
Figure 328665DEST_PATH_IMAGE058
Figure 613016DEST_PATH_IMAGE059
respectively the distributed power generation cost and the energy storage battery power generation cost of the ith data center in the data center cluster B;
Figure 711028DEST_PATH_IMAGE060
Figure 676710DEST_PATH_IMAGE061
distributed generation capital recovery coefficients and battery capital recovery coefficients for the ith data center in the data center cluster B, respectively;
Figure 982927DEST_PATH_IMAGE062
Figure 172600DEST_PATH_IMAGE063
respectively representing the distributed power generation and the sticking rate of the energy storage battery of the ith data center in the data center cluster B;
Figure 760838DEST_PATH_IMAGE064
Figure 326949DEST_PATH_IMAGE065
the service life of distributed power generation equipment and energy storage battery equipment of the ith data center in the data center cluster B respectively;
Figure 566300DEST_PATH_IMAGE066
Figure 113825DEST_PATH_IMAGE067
the installation costs of distributed power generation and energy storage batteries of the ith data center in the data center cluster B are respectively set;
Figure 500944DEST_PATH_IMAGE068
Figure 808429DEST_PATH_IMAGE069
the operation and maintenance costs of the distributed power generation and energy storage battery of the ith data center in the data center cluster B are respectively set;
Figure 699024DEST_PATH_IMAGE070
is an execution time period;
Figure 181565DEST_PATH_IMAGE050
charging and discharging efficiency of an energy storage battery of the ith data center in the data center cluster B;
Figure 728084DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
The execution constraint module is used for receiving the bid winning capacity in the time sequence and receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; establishing power utilization constraint of the ith data center of the data center cluster B in execution time, distributed generation power constraint of the ith data center of the data center cluster B, energy storage battery charging and discharging power upper and lower limit constraint of the ith data center of the data center cluster B, energy storage battery capacity constraint of the ith data center of the data center cluster B and power consumption constraint of a CPU (central processing unit) in the ith data center of the data center cluster B in execution time period, wherein the power consumption constraint is generated by processing calculation tasks and is input to the optimization algorithm module as a data model;
it should be noted that the execution constraint module is configured to receive a scalar in a time sequence, that is, a medium bid amount of virtual output of the data center a in the spot market, energy consumption, distributed power generation data, data related to an energy storage battery required by a CPU computation task from the data center cluster B, and channel data between the data center a and each data center cluster B. For the data center cluster B, corresponding electric quantity is purchased to the power grid in the market at the day before in a contract form, and certain punishment is carried out if the quantity of the contract is violated. Therefore, during the execution period, the condition that the power consumption purchased from each node of the data center cluster B to the power grid is fixed as the total power consumption declared in the market at present
Figure 901576DEST_PATH_IMAGE073
. Intra-execution-time data center groupThe electricity utilization constraint of the ith data center of the set B is shown in a formula (22); the distributed generation power constraint of the ith data center in the data center cluster B is shown as a formula (23); the upper and lower limits of the charging and discharging power of the energy storage battery of the ith data center in the data center cluster B are constrained as shown in a formula (24); the energy storage battery capacity constraint of the ith data center in the data center cluster B is shown in a formula (25); the constraint of the electricity consumption generated by the CPU in the ith data center in the data center cluster B for processing the calculation task in the execution period is shown as a formula (26);
Figure 99208DEST_PATH_IMAGE124
(22)
Figure 4847DEST_PATH_IMAGE125
(23)
Figure 100979DEST_PATH_IMAGE126
(24)
Figure 500999DEST_PATH_IMAGE127
(25)
Figure 366187DEST_PATH_IMAGE128
(26)
in the formula (I), the compound is shown in the specification,
Figure 442727DEST_PATH_IMAGE073
total power consumption declared to the day-ahead market for the ith data center in the data center cluster B in the execution period;
Figure 26155DEST_PATH_IMAGE074
the power consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the execution period;
Figure 462821DEST_PATH_IMAGE075
electric power generated for distributed power generation of the ith data center in the data center cluster B within the execution period;
Figure 385778DEST_PATH_IMAGE076
the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the execution period is shown.
Figure 695537DEST_PATH_IMAGE056
The distributed power generation amount of the ith data center in the data center cluster B in the execution stage in unit time is calculated;
Figure 475184DEST_PATH_IMAGE022
the maximum value of the distributed power generation amount of the ith data center in the data center cluster B in unit time.
Figure 528591DEST_PATH_IMAGE057
The charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in the execution stage in unit time is calculated;
Figure 306054DEST_PATH_IMAGE025
the maximum charging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 973665DEST_PATH_IMAGE026
the maximum discharge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 266106DEST_PATH_IMAGE027
Figure 60886DEST_PATH_IMAGE029
respectively obtaining the minimum capacity and the maximum capacity of an energy storage battery of the ith data center in the data center cluster B;
Figure 443588DEST_PATH_IMAGE080
the initial capacity of an energy storage battery of the ith data center in the data center cluster B before the execution stage is obtained;
Figure 95150DEST_PATH_IMAGE082
the data volume (namely the calculation task to be processed) of the ith data center in the data center cluster B, which is not subjected to calculation task transfer;
Figure 78149DEST_PATH_IMAGE083
the data volume (bit) of the ith data center migrated to the data center cluster B by the data center A in the execution stage;
Figure 660309DEST_PATH_IMAGE041
calculating a CPU period required by a unit data volume (1 bit) for the ith data center in the data center cluster B;
Figure 209102DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B;
Figure 234827DEST_PATH_IMAGE084
calculating the channel consumption power in unit time when the ith data center in the data center cluster B and the data center A in the execution stage perform calculation task transmission;
Figure 187346DEST_PATH_IMAGE085
the time required for transferring the calculation task from the data center A to the ith data center in the data center cluster B in the execution stage is calculated;
Figure 386246DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 727229DEST_PATH_IMAGE086
for the execution of the intra-phase data center A andthe channel capacity between the ith data centers in the data center cluster B, namely the maximum average information rate which can be transmitted by the channels between the ith data centers;
Figure 986172DEST_PATH_IMAGE046
transferring the computing task from the data center A to the channel maximum power consumption of the ith data center in the data center cluster B;
Figure 927452DEST_PATH_IMAGE087
to perform the total amount of data (bit) that data center a migrates to data center cluster B during the phase,
Figure 664464DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
The optimization algorithm module is used for receiving the mathematical models of the target cost function module and the execution constraint module, solving an optimal target cost function value by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving execution resource combination and time sequence output capacity of the data center cluster B in an execution stage;
and reconstructing inequality constraints in the index constraint function into the target function in a form of penalty terms by using the improved efficient Lagrange reconstruction multiplier method.
It should be noted that the optimization algorithm module is configured to receive the mathematical models of the objective cost function module and the execution constraint module, and solve the optimal objective cost function value by using an improved and efficient lagrangian reconstruction multiplier method
Figure 125532DEST_PATH_IMAGE129
And at the same time obtaining the independent variable
Figure 243792DEST_PATH_IMAGE130
I.e. the combination of resources of the execution phase. The time-series output capacity can be obtained by the formula (27):
Figure 219838DEST_PATH_IMAGE131
(27)
wherein the content of the first and second substances,
Figure 698224DEST_PATH_IMAGE089
the time-ordered total capacity of data center cluster B within the execution phase is calculated.
In summary, the embodiment of the present invention provides a demand response system for computing transfer oriented to a coordination data center and a power grid, and a method for computing task transfer is adopted by a data center with different nodes, so as to participate in demand response transaction. In order to enable a data center manager participating in demand response to obtain the maximum benefit, the system also provides an intelligent pricing model based on machine learning, and the benefit maximization is realized. The calculation transfer demand response system for the coordination data center and the power grid has the functions of relieving the pressure of the power system, exploring potential resources, realizing the optimal configuration of response resources and providing a positive effect for the stable and reliable operation of the power system under the conditions of short power supply and power gaps or heavy overload risks of transformers, lines, feeders and the like in local areas.
The terms "first," "second," "third," "fourth," and the like in the description of the application and the above-described figures, if any, are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the application described herein are, for example, capable of operation in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
It should be understood that in the present application, "at least one" means one or more, "a plurality" means two or more. "and/or" for describing an association relationship of associated objects, indicating that there may be three relationships, e.g., "a and/or B" may indicate: only A, only B and both A and B are present, wherein A and B may be singular or plural. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship. "at least one of the following" or similar expressions refer to any combination of these items, including any combination of the singular or plural items. For example, at least one (one) of a, b, or c, may represent: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", wherein a, b, c may be single or plural.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. For example, the above-described apparatus embodiments are merely illustrative, and for example, a division of a unit is merely a logical division, and an actual implementation may have another division, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
Units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit may be implemented in the form of hardware, or may also be implemented in the form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be substantially implemented or contributed to by the prior art, or all or part of the technical solution may be embodied in a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a portable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions in the embodiments of the present application.

Claims (18)

1. A compute-shift demand response system oriented to coordination of a data center and a power grid, comprising: the system comprises a data center A which can participate in power market demand response, a data center cluster B which cannot participate in power market demand response, a declaration calculation module, a dispatching center, a power trading center, an intelligent pricing module and an execution calculation module;
the declaration calculation module is used for receiving the CPU calculation task and the required electric energy of the data center A, and receiving the electric energy, the distributed generation data and the energy storage battery data required by the CPU calculation task of the data center cluster B; establishing a first objective function and a mathematical model of declaration constraint by taking the lowest power consumption generated by processing a calculation task by the data center A in a declaration time period as a target, reconstructing the mathematical model, and solving and processing the reconstructed model by adopting an improved high-efficiency Lagrange reconstruction multiplier method to obtain a time sequence virtual output boundary of the data center A and the power consumption generated by processing the calculation task transferred by the data center A in the data center cluster B;
the intelligent pricing module is used for receiving the time sequence virtual output boundary, the historical market clearing result and the historical execution result, and meanwhile, determining the declared price of the data center A which is optimal for the data center A to participate in the current power market demand response according to the historical market clearing result, the historical execution result and the time sequence virtual output boundary by adopting a Q learning algorithm, so that the optimal time sequence reported price of the current participation in the power market demand response is output;
the dispatching center is used for providing energy demand information of the operation day, and the energy demand information comprises: time series demand capacity and time series demand price;
the electric power trading center is used for receiving the optimal time sequence reported electricity price, the time sequence demand capacity and the time sequence demand price, calculating and outputting a market clearing result, wherein the market clearing result comprises: winning capacity in time sequence of each virtual power plant;
the execution calculation module is used for receiving the bid winning capacity in the time sequence and receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; and establishing a second objective function and a mathematical model of execution constraint in execution time at the lowest cost of a middle scalar quantity which meets the virtual output of the data center A in the spot market, solving an optimal objective cost function value by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving the execution resource combination and the time sequence output capacity of the data center cluster B in an execution stage.
2. The system for computing transfer demand response oriented to coordination of a data center and a power grid according to claim 1, wherein the declaration computing module specifically comprises:
the objective function module is used for receiving the CPU calculation task of the data center A and the required electric energy, establishing the first objective function by taking the lowest electricity consumption generated by the data center A for processing the calculation task in the declaration time period as a target, and inputting the first objective function into the optimization algorithm reconstruction module as a data model;
the declaration constraint module is used for receiving the CPU calculation task and the required electric energy of the data center A, and receiving the electric energy, the distributed generation data and the energy storage battery data required by the CPU calculation task of the data center cluster B; establishing power utilization constraint of the data center cluster B in a declaration time period, charge and discharge quantity upper and lower limit constraint of an energy storage battery of the ith data center in the data center cluster B in unit time, energy storage battery capacity constraint of the ith data center in the data center cluster B, calculation task processing time delay constraint and transmission time delay constraint among data, and inputting the constraints into an optimization algorithm reconstruction module as a data model;
an optimization algorithm reconstruction module, configured to receive the mathematical models of the objective function module and the declaration constraint module, perform mathematical model reconstruction at the same time, reconstruct the mathematical model including inequality constraint and equality constraint into a mathematical model including equality constraint only, obtain a reconstruction model, solve the optimal data amount of the data center a migrated to the ith data center in the data center cluster B and the argument of the first objective function by using an improved high-efficiency lagrangian reconstruction multiplier method, and solve a time sequence virtual output boundary of the data center a and power consumption generated by a CPU of each data center in the data center cluster B due to processing of the calculation task transferred by the data center a;
and reconstructing inequality constraints in an index constraint function into the first objective function in a form of penalty terms by the improved efficient Lagrange reconstruction multiplier method.
3. The coordinated data center and grid oriented computing transfer demand response system of claim 2, wherein the electricity usage constraint comprises: reporting the electricity consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task within a time period; declaring electric power generated by distributed power generation of the ith data center in the data center cluster B within a time period; and reporting the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in a time period.
4. The system for calculating and transferring demand response oriented to coordination of a data center and a power grid as claimed in claim 1, wherein the intelligent pricing module specifically comprises:
the initialization module is used for setting an initialization Q table, wherein the initialization Q table is a two-dimensional table with M rows and 2+ N columns, M represents the number of times that the data center A participates in the historical demand response, N represents the action behavior selection of the data center, the first column of the initialization Q table is used for reporting the electricity price for the historical time sequence of the data center A, the second column of the initialization Q table is used for responding to the historical demand participation of the data center A, and the third column to 2+ N column of the initialization Q table are used for setting future reward expectation values obtained by the data center A for different actions made in time;
a behavior decision module for receiving the time sequence virtual output boundary, the historical market clearing result and the historical execution result, selecting an action Q for the intelligent agent, and acquiring a corresponding state according to the initialized Q table, so as to pass through
Figure 313152DEST_PATH_IMAGE001
Determining an optimal time sequence reporting electricity price by a strategy, and outputting the optimal time sequence reporting electricity price, wherein a historical market clearing result is a scalar in a time sequence, and a historical execution result is a time sequence output capacity and execution resource combination;
and the updating module is used for receiving the reporting electricity price of the optimal time sequence and updating the initialization Q table by using a Bellman equation.
5. The system for computing transfer demand response oriented to coordination of a data center and a power grid according to claim 1, wherein the executing the computing module specifically comprises:
the execution constraint module is used for receiving the bid winning capacity in the time sequence and receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; establishing power utilization constraint of the ith data center of the data center cluster B in execution time, distributed generation power constraint of the ith data center of the data center cluster B, energy storage battery charging and discharging power upper and lower limit constraint of the ith data center of the data center cluster B, energy storage battery capacity constraint of the ith data center of the data center cluster B and power consumption constraint generated by a CPU (central processing unit) processing a calculation task of the ith data center of the data center cluster B in execution time period, and inputting the power consumption constraint into an optimization algorithm module as a data model;
the target cost function module is used for receiving electric energy, distributed generation data and energy storage battery data required by the CPU calculation task of the data center cluster B; establishing a second objective function with the lowest cost of the intermediate scalar satisfying the virtual output of the data center A in the spot market, and inputting the second objective function serving as a data model into an optimization algorithm module;
the optimization algorithm module is used for receiving the mathematical models of the objective cost function module and the execution constraint module, solving an optimal objective cost function value by adopting an improved high-efficiency Lagrange reconstruction multiplier method, and solving execution resource combination and time sequence output capacity of a data center cluster B in an execution stage;
and reconstructing inequality constraints in the index constraint function into the second objective function in a form of penalty terms by using the improved efficient Lagrange reconstruction multiplier method.
6. The coordinated data center and power grid oriented compute transfer demand response system of claim 1, wherein the first objective function of the reconstructed model is:
Figure 794075DEST_PATH_IMAGE002
in the formula (I), the compound is shown in the specification,
Figure 573812DEST_PATH_IMAGE003
the representative reconstruction model aims at minimizing the sum of the power consumption and each penalty term generated by the CPU of the data center A in the declaration time interval due to the processing of the calculation tasks;
Figure 31338DEST_PATH_IMAGE004
an argument representing the first objective function is a total data volume of the reporting period migrated to the data center cluster B by the data center A;
Figure 908027DEST_PATH_IMAGE005
the total data volume is the total data volume when the data center A does not carry out calculation task transfer;
Figure 730490DEST_PATH_IMAGE006
calculating a CPU period required for a unit data volume for the data center A;
Figure 59840DEST_PATH_IMAGE007
calculating the energy consumption of each CPU cycle for the data center A;
Figure 825451DEST_PATH_IMAGE008
Figure 556647DEST_PATH_IMAGE009
penalty coefficients with values greater than 1, respectively;
Figure 612328DEST_PATH_IMAGE010
is an amplification factor, with an odd number greater than 1;
Figure 428974DEST_PATH_IMAGE011
is the j inequality of the ith data center in the data center cluster BAn upper bound penalty term for the bundle;
Figure 729768DEST_PATH_IMAGE012
a lower boundary penalty term of a jth inequality constraint of an ith data center in the data center cluster B;
Figure 253153DEST_PATH_IMAGE013
and
Figure 745314DEST_PATH_IMAGE014
respectively an upper boundary penalty item and a lower boundary penalty item of the distributed power generation amount constraint of the ith data center in the data center cluster B in unit time;
Figure 49256DEST_PATH_IMAGE015
and
Figure 324380DEST_PATH_IMAGE016
respectively an upper boundary penalty item and a lower boundary penalty item of the charge and discharge amount constraint of the energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 30168DEST_PATH_IMAGE017
and
Figure 365334DEST_PATH_IMAGE018
respectively limiting upper and lower boundary punishment items for the capacity of the energy storage battery of the ith data center in the data center cluster B;
Figure 655108DEST_PATH_IMAGE019
an upper bound penalty term for declaring a data capacity constraint of the ith data center in the data center cluster B within the time period;
Figure 796239DEST_PATH_IMAGE020
an upper boundary penalty item is constrained for the channel transmission delay between the data center A and the ith data center in the data center cluster B;
Figure 294217DEST_PATH_IMAGE021
the distributed power generation amount of the ith data center in the data center cluster B in unit time is obtained;
Figure 128180DEST_PATH_IMAGE022
the distributed power generation amount is the maximum value of the distributed power generation amount of the ith data center in the data center cluster B in unit time;
Figure 141136DEST_PATH_IMAGE023
to declare the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period,
Figure 23641DEST_PATH_IMAGE024
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is calculated;
Figure 939907DEST_PATH_IMAGE025
the maximum charging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 882455DEST_PATH_IMAGE026
the maximum discharge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time;
Figure 382706DEST_PATH_IMAGE027
the minimum capacity of an energy storage battery of the ith data center in the data center cluster B is obtained;
Figure 865640DEST_PATH_IMAGE028
the initial capacity of the energy storage battery of the ith data center in the data center cluster B before the reporting time period is reported;
Figure 338210DEST_PATH_IMAGE029
as a data centerThe maximum capacity of an energy storage battery of the ith data center in the cluster B;
Figure 513976DEST_PATH_IMAGE030
reporting the data volume of the ith data center in the data center cluster B which is not subjected to calculation task transfer in the time period, namely the calculation task to be processed by the ith data center;
Figure 439207DEST_PATH_IMAGE031
the data volume of the ith data center migrated to the data center cluster B for the data center A;
Figure 218508DEST_PATH_IMAGE032
reporting the maximum data capacity of the ith data center in the data center cluster B in the time period;
Figure 342322DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 626672DEST_PATH_IMAGE034
calculating the maximum transmission time delay of the ith data center in the cluster B for transferring the task from the data center A to the data center in the declaration time period;
Figure 101516DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
7. The coordinated data center and power grid oriented compute transfer demand response system of claim 6, wherein the equation constraint declaring the electricity usage of data center cluster B over time is:
Figure 129515DEST_PATH_IMAGE036
in the formula (I), the compound is shown in the specification,
Figure 373414DEST_PATH_IMAGE037
reporting the total power consumption of the ith data center in the data center cluster B to the market in the day ahead for a reporting time period;
Figure 563087DEST_PATH_IMAGE038
reporting the electricity consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the time period;
Figure 26692DEST_PATH_IMAGE039
generating power for distributed generation of the ith data center in the data center cluster B within the declaration period;
Figure 592803DEST_PATH_IMAGE023
and reporting the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period.
8. The coordinated data center and grid oriented compute transfer demand response system of claim 7, wherein variables are variables
Figure 956788DEST_PATH_IMAGE038
Is constrained by:
Figure 379679DEST_PATH_IMAGE040
in the formula (I), the compound is shown in the specification,
Figure 766798DEST_PATH_IMAGE041
calculating the CPU period required by unit data volume for the ith data center in the data center cluster B in the declaration time period;
Figure 198916DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 89512DEST_PATH_IMAGE043
calculating channel consumption power when an ith data center in a data center cluster B and a data center A perform task transmission in unit time;
Figure 447419DEST_PATH_IMAGE044
calculating the time required for transferring the task from the data center A to the ith data center in the data center cluster B within the declaration time period;
Figure 56254DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 292064DEST_PATH_IMAGE045
reporting the channel capacity between the data center A and the ith data center in the data center cluster B in a time period, namely the maximum average information rate which can be transmitted by the channel between the data center A and the ith data center;
Figure 365062DEST_PATH_IMAGE046
transferring the computing task from the data center A to the channel maximum power consumption of the ith data center in the data center cluster B;
Figure 333018DEST_PATH_IMAGE004
reporting the total data volume of the data center A migrated to the data center cluster B in the time period;
Figure 491467DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
9. The coordinated data center and grid oriented compute transfer demand response system of claim 8, wherein variables are variables
Figure 766853DEST_PATH_IMAGE039
The equation constraints are:
Figure 632040DEST_PATH_IMAGE047
in the formula (I), the compound is shown in the specification,
Figure 833215DEST_PATH_IMAGE039
to declare the power generated by distributed generation of the ith data center in data center cluster B over the period of time,
Figure 416643DEST_PATH_IMAGE021
distributed power generation per unit time for the ith data center in the data center cluster B,
Figure 728675DEST_PATH_IMAGE048
to declare a time period.
10. The coordinated data center and grid oriented compute transfer demand response system of claim 9, wherein variables are variables
Figure 776266DEST_PATH_IMAGE023
Is constrained by:
Figure 86024DEST_PATH_IMAGE049
in the formula (I), the compound is shown in the specification,
Figure 723460DEST_PATH_IMAGE023
to declare the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the time period,
Figure 776867DEST_PATH_IMAGE024
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in unit time is calculated;
Figure 678964DEST_PATH_IMAGE050
and the charging and discharging efficiency of the energy storage battery of the ith data center in the data center cluster B is obtained.
11. The system according to claim 2, wherein the solving of the time-series virtual output boundary of the data center a and the power consumption of the CPUs of the data centers in the data center cluster B due to processing the calculation task transferred by the data center a specifically comprises:
calculating a time sequence virtual output boundary through a time sequence virtual output boundary calculation formula, and calculating power consumption through a power consumption calculation formula;
the time sequence virtual output boundary calculation formula is as follows:
Figure 159624DEST_PATH_IMAGE051
the calculation formula of the power consumption is as follows:
Figure 514381DEST_PATH_IMAGE052
in the formula (I), the compound is shown in the specification,
Figure 433796DEST_PATH_IMAGE053
a timing virtual force boundary;
Figure 128082DEST_PATH_IMAGE054
the electricity consumption is;
Figure 343425DEST_PATH_IMAGE004
data volume for migrating the computing tasks in the data center A to the data center cluster B;
Figure 388742DEST_PATH_IMAGE006
for data center ACalculating the CPU period required by unit data amount;
Figure 846268DEST_PATH_IMAGE007
calculating the energy consumption of each CPU cycle for the data center A;
Figure 395061DEST_PATH_IMAGE031
the data volume of the ith data center migrated to the data center cluster B for the data center A;
Figure 545420DEST_PATH_IMAGE041
calculating the CPU period required by unit data volume for the ith data center in the data center cluster B in the declaration time period;
Figure 874770DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle for the ith data center in the data center cluster B in the declaration time period;
Figure 73670DEST_PATH_IMAGE043
calculating channel power consumption when the ith data center in the data center cluster B and the data center A perform calculation task transmission in unit time;
Figure 37821DEST_PATH_IMAGE044
the time required for the transfer of the computing task from data center a to the ith data center in data center cluster B within the declaration period.
12. The coordinated data center and grid oriented compute transfer demand response system of claim 5, wherein the second objective function is:
Figure 296764DEST_PATH_IMAGE055
in the formula (I), the compound is shown in the specification,
Figure 113411DEST_PATH_IMAGE056
Figure 850423DEST_PATH_IMAGE057
distributed generating capacity in unit time and charging and discharging capacity of an energy storage battery in unit time of the ith data center in the data center cluster B in the execution time are respectively also independent variables of the second objective function;
Figure 498442DEST_PATH_IMAGE058
Figure 928286DEST_PATH_IMAGE059
respectively the distributed power generation cost and the energy storage battery power generation cost of the ith data center in the data center cluster B;
Figure 468114DEST_PATH_IMAGE060
Figure 8817DEST_PATH_IMAGE061
distributed generation capital recovery coefficients and battery capital recovery coefficients for the ith data center in the data center cluster B, respectively;
Figure 714605DEST_PATH_IMAGE062
Figure 112088DEST_PATH_IMAGE063
respectively representing the distributed power generation and the sticking rate of the energy storage battery of the ith data center in the data center cluster B;
Figure 841009DEST_PATH_IMAGE064
Figure 982141DEST_PATH_IMAGE065
the service life of distributed power generation equipment and energy storage battery equipment of the ith data center in the data center cluster B respectively;
Figure 58548DEST_PATH_IMAGE066
Figure 564616DEST_PATH_IMAGE067
the installation costs of distributed power generation and energy storage batteries of the ith data center in the data center cluster B are respectively set;
Figure 843151DEST_PATH_IMAGE068
Figure 787973DEST_PATH_IMAGE069
the operation and maintenance costs of the distributed power generation and energy storage battery of the ith data center in the data center cluster B are respectively set;
Figure 140457DEST_PATH_IMAGE070
is an execution time period;
Figure 145322DEST_PATH_IMAGE050
charging and discharging efficiency of an energy storage battery of the ith data center in the data center cluster B;
Figure 583257DEST_PATH_IMAGE035
the number of data centers included in data center cluster B;
Figure 567655DEST_PATH_IMAGE071
the lowest sum of the distributed power generation cost of the data center cluster B and the power generation cost of the energy storage battery in the execution period is taken as a target.
13. The coordinated data center and power grid oriented compute transfer demand response system of claim 5, wherein the power usage constraint is:
Figure 40225DEST_PATH_IMAGE072
in the formula (I), the compound is shown in the specification,
Figure 950412DEST_PATH_IMAGE073
total power consumption declared to the day-ahead market for the ith data center in the data center cluster B in the execution period;
Figure 937960DEST_PATH_IMAGE074
the power consumption generated by the CPU of the ith data center in the data center cluster B for processing the calculation task in the execution period;
Figure 162268DEST_PATH_IMAGE075
electric power generated for distributed power generation of the ith data center in the data center cluster B within the execution period;
Figure 551661DEST_PATH_IMAGE076
the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the execution period is shown.
14. The coordinated data center and grid oriented compute transfer demand response system of claim 5, wherein the distributed generation power constraint is:
Figure 836012DEST_PATH_IMAGE077
in the formula (I), the compound is shown in the specification,
Figure 809390DEST_PATH_IMAGE056
the distributed power generation amount of the ith data center in the data center cluster B in the execution stage in unit time is calculated;
Figure 571810DEST_PATH_IMAGE022
the maximum distributed power generation amount per unit time of the ith data center in the data center cluster B.
15. The system for calculating and transferring demand response oriented to the coordination data center and the power grid according to claim 5, wherein the upper and lower limits of the charge and discharge power of the energy storage battery are constrained as follows:
Figure 815710DEST_PATH_IMAGE078
in the formula (I), the compound is shown in the specification,
Figure 333279DEST_PATH_IMAGE057
the charging and discharging amount of an energy storage battery of the ith data center in the data center cluster B in the execution stage in unit time is calculated;
Figure 233101DEST_PATH_IMAGE025
the maximum charge capacity of an energy storage battery of the ith data center in the data center cluster B in unit time is obtained;
Figure 861529DEST_PATH_IMAGE026
the maximum discharge amount of the energy storage battery of the ith data center in the data center cluster B in unit time.
16. The coordinated data center and grid oriented compute transfer demand response system of claim 5, wherein the energy storage battery capacity constraints are:
Figure 897618DEST_PATH_IMAGE079
in the formula (I), the compound is shown in the specification,
Figure 87553DEST_PATH_IMAGE027
Figure 536989DEST_PATH_IMAGE029
respectively, of the ith data center in data center cluster BEnergy storage battery minimum and maximum capacities;
Figure 641211DEST_PATH_IMAGE080
the initial capacity of an energy storage battery of the ith data center in the data center cluster B before the execution stage is obtained;
Figure 859703DEST_PATH_IMAGE076
the charging and discharging amount of the energy storage battery of the ith data center in the data center cluster B in the execution period is shown.
17. The coordination data center and grid oriented compute transfer demand response system of claim 5, wherein the power usage constraint is:
Figure 656758DEST_PATH_IMAGE081
in the formula (I), the compound is shown in the specification,
Figure 327910DEST_PATH_IMAGE082
the data volume of the ith data center in the data center cluster B, which is not subjected to the calculation task transfer, is obtained;
Figure 68114DEST_PATH_IMAGE083
the data volume of the ith data center in the data center cluster B migrated to the data center A in the execution stage;
Figure 78796DEST_PATH_IMAGE041
calculating a CPU period required by unit data volume for the ith data center in the data center cluster B;
Figure 109069DEST_PATH_IMAGE042
calculating the energy consumption of each CPU cycle by the ith data center in the data center cluster B;
Figure 205200DEST_PATH_IMAGE084
calculating the channel consumption power in unit time when the ith data center in the data center cluster B and the data center A in the execution stage perform calculation task transmission;
Figure 713542DEST_PATH_IMAGE085
the time required for transferring the calculation task from the data center A to the ith data center in the data center cluster B in the execution stage is calculated;
Figure 906626DEST_PATH_IMAGE033
the optical fiber transmission speed of the channel between each node;
Figure 45483DEST_PATH_IMAGE086
the channel capacity between the data center A and the ith data center in the data center cluster B in the execution stage, namely the maximum average information rate which can be transmitted by the channel between the data center A and the ith data center in the data center cluster B;
Figure 192693DEST_PATH_IMAGE046
transferring the computing task from the data center A to the channel maximum power consumption of the ith data center in the data center cluster B;
Figure 442409DEST_PATH_IMAGE087
to perform the total amount of data that data center a migrates to data center cluster B during the phase,
Figure 489999DEST_PATH_IMAGE035
the number of data centers included in data center cluster B.
18. The system of claim 5, wherein the solving of the execution resource combination and the time-series output capacity of the data center cluster B in the execution phase specifically comprises:
substituting the execution resource combination into a time sequence output capacity calculation formula to calculate and obtain the sum time sequence output capacity of the data center cluster B in the execution stage;
wherein, the time sequence output capacity calculation formula is as follows:
Figure 862075DEST_PATH_IMAGE088
in the formula (I), the compound is shown in the specification,
Figure 932799DEST_PATH_IMAGE089
time-series capacity for the data center cluster B in the execution phase;
Figure 48523DEST_PATH_IMAGE035
the number of data centers included in data center cluster B;
Figure 888303DEST_PATH_IMAGE056
Figure 929815DEST_PATH_IMAGE057
the distributed power generation amount per unit time and the charging and discharging amount of the energy storage battery per unit time of the ith data center in the data center cluster B in the execution time are also independent variables of the second objective function, namely the execution resource combination.
CN202211002763.0A 2022-08-22 2022-08-22 Calculation transfer demand response system oriented to coordination data center and power grid Active CN115081758B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211002763.0A CN115081758B (en) 2022-08-22 2022-08-22 Calculation transfer demand response system oriented to coordination data center and power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211002763.0A CN115081758B (en) 2022-08-22 2022-08-22 Calculation transfer demand response system oriented to coordination data center and power grid

Publications (2)

Publication Number Publication Date
CN115081758A true CN115081758A (en) 2022-09-20
CN115081758B CN115081758B (en) 2023-01-03

Family

ID=83243983

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211002763.0A Active CN115081758B (en) 2022-08-22 2022-08-22 Calculation transfer demand response system oriented to coordination data center and power grid

Country Status (1)

Country Link
CN (1) CN115081758B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150278968A1 (en) * 2009-10-23 2015-10-01 Viridity Energy, Inc. Facilitating revenue generation from data shifting by data centers
CN111461920A (en) * 2020-03-23 2020-07-28 青海大学 Smart grid-oriented data center demand response method and system
CN113078665A (en) * 2021-04-07 2021-07-06 国网上海能源互联网研究院有限公司 Method and system for optimizing energy storage system based on data center demand response
WO2021135332A1 (en) * 2019-12-31 2021-07-08 东北大学 Edge computing-based multi-agent load regulation and control method
CN113177655A (en) * 2021-03-25 2021-07-27 天津大学 Comprehensive energy system multi-main-body operation optimization method and device based on reinforcement learning
CN114662319A (en) * 2022-03-25 2022-06-24 华北电力大学 Construction method of active power distribution network planning model considering data center

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150278968A1 (en) * 2009-10-23 2015-10-01 Viridity Energy, Inc. Facilitating revenue generation from data shifting by data centers
WO2021135332A1 (en) * 2019-12-31 2021-07-08 东北大学 Edge computing-based multi-agent load regulation and control method
CN111461920A (en) * 2020-03-23 2020-07-28 青海大学 Smart grid-oriented data center demand response method and system
CN113177655A (en) * 2021-03-25 2021-07-27 天津大学 Comprehensive energy system multi-main-body operation optimization method and device based on reinforcement learning
CN113078665A (en) * 2021-04-07 2021-07-06 国网上海能源互联网研究院有限公司 Method and system for optimizing energy storage system based on data center demand response
CN114662319A (en) * 2022-03-25 2022-06-24 华北电力大学 Construction method of active power distribution network planning model considering data center

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HUANG,XIN等: "Optimal Dispatch and Control Strategy of Integrated Energy System Considering Multiple P2H to Provide Integrated Demand Response", 《FRONTIERS IN ENERGY RESEARCH》 *
王晴等: "考虑数据中心需求响应的城市电网阻塞管理", 《电网技术》 *

Also Published As

Publication number Publication date
CN115081758B (en) 2023-01-03

Similar Documents

Publication Publication Date Title
US10317924B2 (en) Collaborative load balancing within a community of energy nodes
CN108599373B (en) Cascade analysis method for transmission and distribution coordination scheduling target of high-proportion renewable energy power system
Elkazaz et al. A hierarchical and decentralized energy management system for peer-to-peer energy trading
CN113610311B (en) Comprehensive energy service provider cooperation operation optimization method considering carbon emission reduction under double-layer cooperation architecture
CN113393125A (en) Comprehensive energy system cooperative scheduling method based on source-load bilateral interactive game
Zhang et al. Optimization strategy for power sharing and low-carbon operation of multi-microgrid IES based on asymmetric nash bargaining
CN114256836A (en) Capacity optimization configuration method for new energy power station shared energy storage
Shi et al. Research on power to hydrogen optimization and profit distribution of microgrid cluster considering shared hydrogen storage
EP2251822A1 (en) A method for load balancing in a power system
CN116342166A (en) Energy game regulation and control method and equipment based on multi-region sharing
Liu et al. Optimal dispatch strategy of virtual power plants using potential game theory
CN115102231A (en) Wind and light storage station optimization control method and system under multi-scale electricity-carbon mode
Lasemi et al. Energy cost optimization of globally distributed internet data centers by copula-based multidimensional correlation modeling
Dai et al. An equilibrium model of the electricity market considering the participation of virtual power plants
CN115081758B (en) Calculation transfer demand response system oriented to coordination data center and power grid
CN112884381A (en) P2P energy-consuming market planning method considering supply and demand uncertainty
CN116432862A (en) Multi-main-body game optimization method and device for renewable energy micro-grid
CN116436048A (en) Multi-target-driven micro-grid group cloud energy storage optimal configuration method and device
Samadi et al. Energy blockchain for demand response and distributed energy resource management
CN115018554A (en) Virtual power plant secondary frequency modulation method considering comprehensive energy resources
Saukh et al. Modelling of market equilibrium on the basis of Smart Grid market system decomposition
Zhang et al. Joint optimal operation and bidding strategy of scenic reservoir group considering energy storage sharing
CN113177655A (en) Comprehensive energy system multi-main-body operation optimization method and device based on reinforcement learning
Minhas et al. Optimal control of power cost and consumer satisfaction using smart grid intelligent energy management system
Spasova et al. Energy exchange strategy for local energy markets with heterogenous renewable sources

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant