CN112232984B - Distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method - Google Patents
Distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method Download PDFInfo
- Publication number
- CN112232984B CN112232984B CN202011102939.0A CN202011102939A CN112232984B CN 112232984 B CN112232984 B CN 112232984B CN 202011102939 A CN202011102939 A CN 202011102939A CN 112232984 B CN112232984 B CN 112232984B
- Authority
- CN
- China
- Prior art keywords
- energy
- comprehensive
- energy system
- power
- iii
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000005457 optimization Methods 0.000 title claims abstract description 19
- 238000011156 evaluation Methods 0.000 claims abstract description 22
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 230000008878 coupling Effects 0.000 claims description 10
- 238000010168 coupling process Methods 0.000 claims description 10
- 238000005859 coupling reaction Methods 0.000 claims description 10
- 239000005431 greenhouse gas Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000002910 solid waste Substances 0.000 claims description 9
- 230000003749 cleanliness Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012423 maintenance Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000007797 corrosion Effects 0.000 claims description 4
- 238000005260 corrosion Methods 0.000 claims description 4
- 238000005336 cracking Methods 0.000 claims description 4
- 230000006378 damage Effects 0.000 claims description 4
- 238000004140 cleaning Methods 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000011109 contamination Methods 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 abstract description 5
- 239000002699 waste material Substances 0.000 abstract description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 24
- 239000003345 natural gas Substances 0.000 description 12
- 239000010410 layer Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009440 infrastructure construction Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 239000002918 waste heat Substances 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013278 delphi method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/82—Energy audits or management systems therefor
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Artificial Intelligence (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a distributed data central computing power and energy flow fused comprehensive energy system optimization scheduling method, which comprises the following steps: s1, establishing a comprehensive energy system mathematical model integrating distributed data central computing power and energy flow; s2, establishing a comprehensive energy system operation evaluation index system comprising three primary indexes of economy, safety and cleanness and a plurality of secondary indexes; s3, determining the comprehensive weight of each index by adopting a comprehensive evaluation method; s4, constructing an optimized dispatching model of the comprehensive energy system by taking the lowest operation cost, the highest safety and the lowest pollution emission as three objective functions; and S5, accessing the optimized dispatching model of the comprehensive energy system into the mathematical model of the comprehensive energy system to obtain an optimal dispatching method and a result. The distributed data center is incorporated into the comprehensive energy system, and the comprehensive optimization scheduling is carried out on the other energy sources such as computing power, electric power, heating power and the like on the whole, so that the energy waste is reduced, the consumption of clean energy is increased, and the flexibility of the system is improved.
Description
Technical Field
The invention belongs to the field of comprehensive energy systems, and particularly relates to a distributed data central computing power and energy flow fused comprehensive energy system optimization scheduling method.
Background
The distributed data center is used as an upgrading scheme of a centralized data center and is an information infrastructure of a next generation ultra-high speed network. At present, the construction of a new generation of information infrastructure, namely 'new infrastructure', is enhanced in China, and the information infrastructure construction is a main edition block, such as communication network infrastructure represented by 5G, internet of things, industrial Internet and satellite Internet, new technical infrastructure represented by artificial intelligence, cloud computing, block chains and the like, computing infrastructure represented by a data center and an intelligent computing center and the like. Distributed data centers have become an integral part of future industrial park infrastructure construction. However, the introduction of the distributed data center will inevitably increase the power consumption of the area greatly, and will cause a great impact on the planning and operation optimization of the comprehensive energy system in the area. Meanwhile, the power calculating equipment generates a large amount of heat during operation, and how to utilize the heat is one of the focuses of people.
Although many technologies for providing heating or hot water supply service for local users by using waste heat of computing equipment as a heat source appear at the present stage, in the technologies, a distributed data center is often independent of a local regional comprehensive energy system in the aspects of design, planning, operation and control, and cannot be comprehensively considered, and situations of untimely regulation and control and unmatched supply and demand often appear.
The distributed data center is incorporated into the comprehensive energy system, so that unified optimization scheduling is performed on other clean energy such as computing power, electric power, heating power and the like on the whole, the problems of untimely regulation and control, mismatching of supply and demand and the like can be well solved, energy waste is reduced, and the flexibility of the system is improved.
Disclosure of Invention
In order to solve the problems in the prior art and optimally schedule the integrated energy system integrating distributed data center computing power and energy flow, the invention discloses an optimized scheduling method of the integrated energy system integrating distributed data center computing power and energy flow.
The invention specifically adopts the following technical scheme:
a distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method comprises the following specific steps:
s1, establishing a comprehensive energy system mathematical model integrating distributed data central computing power and energy flow;
s2, establishing a comprehensive energy system operation evaluation index system comprising three primary indexes of economy, safety and cleanliness and a plurality of secondary indexes;
s3, determining the comprehensive weight of each index by adopting a comprehensive evaluation method;
s4, constructing an optimized scheduling model of the comprehensive energy system by taking the lowest operation cost, the highest safety and the lowest pollution emission as three objective functions;
and S5, accessing the optimized dispatching model of the comprehensive energy system into the mathematical model of the comprehensive energy system to obtain an optimal dispatching method and a result.
Further, the step S1 of establishing a mathematical model of the integrated energy system with integrated distributed data central computing power and energy flow includes:
the comprehensive energy system is modeled based on graph theory, the joints of equipment such as a power plant, a solar photovoltaic panel, a distributed data center, a transformer, a heat source, a heating power station, a user terminal and the like are abstracted into connection nodes, lines between the two nodes are abstracted into edges, each piece of equipment in the system is simplified into a physical model according to actual system data, and a topological structure of the comprehensive energy system is established.
Further, for a single energy conversion device, the relationship between the input IN and the output OUT is:
OUT=μ·IN (1)
where μ represents a coupling coefficient between the input and the output;
for an integrated energy system comprising a plurality of energy transforming devices and a plurality of forms of energy, the relationship between input and output may be described by a matrix:
in the formula, subscripts P, G, H, C represent electric power, clean energy, thermal power, and computational power, respectively; mu.s PG The coupling coefficient between the input and the output when the electric power is used as the input and the clean energy is used as the output; mu.s GP The coupling coefficient between the input and the output is determined by taking clean energy as input and taking electric power as output.
Further, the operation evaluation system of the comprehensive energy system in the step S2 mainly comprises three primary indexes of economy, safety and cleanness;
wherein the economic indicator I comprises energy material cost I 1 Equipment operation maintenance cost I 2 Energy loss cost I 3 Three secondary indexes;
safety index II includes equipment failure rate II 1 And probability of fatigue cracking of heat supply network pipeline II 2 Corrosion failure rate II 3 Three secondary indexes;
cleanliness indices III include solid waste discharge III 1 Clean energy utilization rate III 2 Greenhouse gas emission III 3 Three secondary indexes.
Further, in the step S3, the comprehensive weight of each index is determined by using a comprehensive evaluation method, wherein an entropy weight method is used as an objective evaluation method, and an analytic hierarchy process is used as a subjective evaluation method. The entropy weight method obtains the objective weight of each index asThe analytic hierarchy process obtains the subjective weight of each index asThen its composite weight isi=I,I 1 ,I 2 ,I 3 ,II,II 1 ,II 2 ,II 3 ,III,III 1 ,III 2 ,III 3 。
Further, the economic indicator I represents a function x (I) as:
the safety index II represents a function x (II) of:
the cleanliness index III represents a function x (III) of:
in the formulae (3), (4) and (5), w i The weight of each index is represented, and x (i) represents each influence index.
Further, the influence indexes are as follows:
x(I 1 )、x(I 2 ) And x (I) 3 ) Respectively as follows:
in the formula (I), the compound is shown in the specification,respectively represent the prices for purchasing electric power, clean energy (natural gas, wind energy, solar energy, etc.); z P ,Z G Representing the amount of purchased electricity, clean energy (natural gas, wind, solar, etc.); gamma ray P ,γ G ,γ H ,γ C Respectively representing the unit cost of operation and maintenance of electric equipment, clean energy (natural gas, wind energy, solar energy and the like), thermal equipment and power calculation equipment, and converting the unit cost into unit power; p P ,P G ,P H ,P C Respectively representing electric power equipment and clean energy (day)Natural gas, wind energy, solar energy, etc.) equipment, thermal equipment, computational power equipment; delta phi P ,ΔΦ G ,ΔΦ H ,ΔΦ C Respectively represents the power loss, the loss of clean energy (natural gas, wind energy, solar energy and the like), the thermal loss and the computational power loss in the operation process.
Due to influence on equipment failure rate II 1 And probability of fatigue cracking of heat supply network pipeline II 2 Corrosion failure rate II 3 The three indexes have many factors, mainly including running time, medium temperature, performance coefficients of equipment, external environment factors and artificial destruction probability, so x (II) 1 )、x(II 2 ) And x (II) 3 ) The complex functions of (a) are:
in the formula, f () can be determined according to practical conditions based on experience, T represents the running time, p represents the bearing pressure of the equipment, T represents the temperature of a medium, alpha represents each performance coefficient of the equipment, beta represents an external environment factor, and lambda represents the probability of artificial damage.
x(III 1 )、x(III 2 ) And x (III) 3 ) Respectively as follows:
wherein epsilon represents the degree of contamination of the solid waste; s Solid body ,S Cleaning of ,S General (1) Respectively representing the solid waste discharge amount, the clean energy utilization amount and the total energy utilization amount; p Total output of The total generated heat and calculated force of the system are represented and converted into electric power; l is Greenhouse gases Indicating greenhouse gas emissions.
Further, in the step S4, the lowest operation cost, the highest safety, and the lowest pollution emission are used as three objective functions to construct an optimized scheduling model of the integrated energy system, which specifically includes:
an objective function f (x, w) of the comprehensive energy system optimization scheduling model:
constraint conditions of the model:
H=H demand for ,C=C Demand for (11)
In the formula, P, T, I and Y respectively represent the power, temperature, medium flow (current, heat flow, clean energy (natural gas, wind energy, solar energy and the like) flow and calculation flow) of each device in the comprehensive energy system and the bearable pressure of the devices; the subscripts min and max represent the minimum and maximum values of the index, respectively; h and C represent total heat and total computing power generated by the system; h Demand for ,C Demand for Is a constant representing the user's thermal and computational demands of the dispatch.
Further, in step S5, the comprehensive energy system optimized scheduling model is accessed to the mathematical model thereof to obtain an optimal scheduling method and result, and the specific method is as follows:
solving the optimal solution set of the model by adopting a non-dominated sorting genetic algorithm (NSGA-II) with an elite strategy to obtain a Pareto optimal solution set: f (x, w) = { f 1 (x,w),f 2 (x,w),…,f n (x, w) }, when y = { y = 1 ,y 2 ,…,y n Y is an operation parameter set of each system device in the scheduling process;
and selecting a solution closest to the optimal scheme from the Pareto optimal solution set by adopting a TOPSIS method to obtain an optimal scheduling result, wherein the weights of the three target functions are w respectively I ,w II ,w III 。
The beneficial effects of the invention are:
according to the comprehensive energy system mathematical model integrating the calculation capacity and the energy flow of the distributed data center, the calculation capacity of the distributed data center is contained in the comprehensive energy system to form one of the generalized energy flow carriers, so that on one hand, the waste heat emitted by the distributed data center can be better recovered, the energy waste is reduced, on the other hand, the consumption of clean energy in the comprehensive energy system can be increased by introducing the distributed data center, and the flexibility of the system is improved. Meanwhile, a comprehensive evaluation method is adopted, a comprehensive evaluation system of the comprehensive energy system is established, reasonable combination of the subjective and objective weight method is achieved, and the obtained result is more accurate.
Drawings
FIG. 1 is the main steps of the process of the present invention;
FIG. 2 is a diagram of a comprehensive evaluation system of the comprehensive energy system;
FIG. 3 is a schematic diagram of an integrated energy system energy flow with distributed data centric power and energy flow fusion;
FIG. 4 is a schematic diagram of a primary energy conversion process;
FIG. 5 is a flow chart of an optimization algorithm.
Detailed Description
The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings and embodiments. The technical scope of the present invention is not limited to the contents of the specification, and the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments, and various changes and modifications can be made by workers within the scope of the technical spirit of the present invention without departing from the technical scope of the present invention, and thus the technical scope thereof must be determined by the scope of the claims.
The following further illustrates the optimized scheduling method of the present invention by embodiments:
as shown in fig. 1, a distributed data center computing power and energy flow integrated comprehensive energy system optimization scheduling method specifically includes the following steps:
the step S1 of establishing a comprehensive energy system mathematical model integrating distributed data central computing power and energy flow comprises the following specific steps:
as shown in fig. 3, the integrated energy system is modeled based on graph theory, the joints of devices such as a power plant, a solar photovoltaic panel, a distributed data center, a transformer, a heat source, a heating power station, a user terminal and the like are abstracted into connection nodes, lines between the two nodes are abstracted into edges, each device in the system is simplified into a physical model according to actual system data, and a topological structure of the integrated energy system is established.
For a single energy conversion device, the relationship between the input IN and the output OUT is:
OUT=μ·IN (1)
where μ represents a coupling coefficient between the input and the output;
as shown in fig. 4, the integrated energy system including a plurality of energy transforming devices and a plurality of energy forms can describe the relationship between the input and the output by a matrix:
wherein the subscript P, G, H, C represents electric power, clean energy, thermal power, computational power, μ PG The coupling coefficient between the input and the output when the electric power is used as the input and the clean energy is used as the output; mu.s GP The coupling coefficient between the input and output is determined by taking clean energy as input and taking electric power as output.
As shown in fig. 2, the operation evaluation system of the integrated energy system in step S2 mainly includes three primary indexes of economy, safety and cleanliness;
wherein the economic indicator I comprises energy material cost I 1 Equipment operating and maintenance cost I 2 Energy loss cost I 3 Three secondary indexes;
safety index II includes the failure rate of the equipment II 1 And the fatigue cracking probability II of the heat supply network pipeline 2 Corrosion failure rate II 3 Three secondary indexes;
cleanliness indices III include solid waste discharge III 1 Clean energy utilization rate III 2 Greenhouse gas emission III 3 Three secondary indexes.
Wherein, the step S3 adopts a comprehensive evaluation method to determine the comprehensive weight of each index, and the comprehensive evaluation method comprises an objective evaluation method and a subjective evaluation methodAnd a valence method, wherein an entropy weight method is adopted as an objective evaluation method, and an analytic hierarchy process is adopted as a subjective evaluation method. Obtaining objective weight of each index by entropy weight methodObtaining subjective weight of each index by using analytic hierarchy processThen its composite weight isi=I,I 1 ,I 2 ,I 3 ,II,II 1 ,II 2 ,II 3 ,III,III 1 ,III 2 ,III 3 。
The specific method comprises the following steps:
an entropy weight method is used to determine the objective (quantitative) weight of the selected indicator.
Firstly, normalization processing is carried out on each secondary index data, and the method comprises the following steps:
for example, the jth secondary index of the ith scheme is x ij Normalized to x i ' j The following two equations may be used for normalization:
If the index is a forward index, selecting a first formula; if the indicator is a negative indicator, then the second formula is selected. Where max (x) j ) Is the maximum value of the j index, min (x) j ) Is the minimum value of the j-th index. Obtaining a comprehensive evaluation index matrix X ' = (X ' consisting of m scheme n indexes after the treatment is finished ' ij ) m×n ,i=1,2,…,m;j=1,2,…,n。
When the second-level indexes are normalized, the qualitative indexes in the second-level indexes are converted into quantitative indexes by selecting 4 grades of index sets {20%,40%,60% and 80% } which respectively represent poor, good and good. Wherein the grade of the qualitative index can be determined according to expert opinions and empirical formulas.
Then, the entropy and the weight of each index are calculated, and the method comprises the following steps:
1) Calculating the specific gravity of the jth index of the ith user:
2) Calculating the information entropy of the j index:
And determining the subjective (qualitative) weight of the evaluation index by combining an analytic hierarchy process with a Delphi method.
Firstly, a hierarchical structure model is established, and the structural layer comprises a target layer, a criterion layer and a scheme layer.
And (3) constructing a pair comparison matrix, starting from the second layer, using the pair comparison matrix and the scale from 1 to 9, wherein the construction method of each layer is as follows:
wherein, a i,j Represents the ithThe relative importance of the indicator to the factor of the previous layer relative to the jth indicator.
Then, calculating a single-rank-order weight vector and performing consistency check, namely:
and calculating the maximum eigenvalue and the corresponding eigenvector of each pair of comparison matrixes, and performing consistency check by using the consistency index, the random consistency index and the consistency ratio. If the test is passed, the feature vector is normalized to be a weight vector; if not, the comparison matrix needs to be reconstructed. Wherein:
the consistency index is defined as:
The random consistency index is RI, and the random consistency index is obtained by table look-up.
The consistency ratio is:
the identity test was considered passed when CR < 0.1.
And finally, calculating a total sorting vector and carrying out consistency check, wherein the total sorting consistency check method is consistent with the single-layer sorting consistency check method. After passing the test, the subjective (qualitative) weight can be obtained:
the objective (quantitative) weight and the subjective (qualitative) weight of each index are combined to obtain the comprehensive weight of each index.
w=(w 1 ,w 2 ,...,w j ,…,w n ) (8)
Wherein w isComprehensive weight vector matrix, w j Is the integrated weight vector of the jth index,
the economic indicator I represents a function x (I) as follows:
the safety index II represents a function x (II) of:
the cleanliness index III represents a function x (III) of:
in the formulae (10), (11) and (12), w i The weight of each index is represented, and x (i) represents each influence index.
The above-mentioned various influence indexes are:
x(I 1 )、x(I 2 ) And x (I) 3 ) Respectively as follows:
in the formula (I), the compound is shown in the specification,respectively represent the prices for purchasing electricity, clean energy (natural gas, wind energy, solar energy, etc.); z P ,Z G Representing the amount of purchased electricity, clean energy (natural gas, wind, solar, etc.); gamma ray P ,γ G ,γ H ,γ C Respectively representing power equipment and clean energy (natural gas and wind energy)Solar energy, etc.) equipment, thermodynamic equipment, computational power equipment, converted to unit power; p P ,P G ,P H ,P C Respectively representing the working power of electric equipment, clean energy (natural gas, wind energy, solar energy and the like) equipment, thermal equipment and computing power equipment; delta phi P ,ΔΦ G ,ΔΦ H ,ΔΦ C Respectively represents the power loss, the loss of clean energy (natural gas, wind energy, solar energy and the like), the thermal loss and the computing power loss in the operation process.
x(II 1 )、x(II 2 ) And x (II) 3 ) Respectively as follows:
in the formula, f () can be determined according to practical conditions according to experience (in the case of the scheme, coefficients of influencing factors are obtained through multiple experiments by adopting a least square method), T represents running time, p represents bearing pressure of equipment, T represents medium temperature, alpha represents each performance coefficient of the equipment, beta represents an external environment factor, and lambda represents artificial destruction probability.
x(III 1 )、x(III 2 ) And x (III) 3 ) Respectively as follows:
wherein epsilon represents the degree of contamination of the solid waste; s Solid body ,S Cleaning of ,S General assembly Respectively representing the solid waste discharge amount, the clean energy utilization amount and the total energy utilization amount; p Total yield The total generated heat and calculated force of the system are expressed and converted into electric power; l is Greenhouse gases Indicating greenhouse gas emissions.
S4, constructing an optimized dispatching model of the comprehensive energy system by taking the lowest operation cost, the highest safety and the lowest pollution emission as three objective functions, wherein the specific method comprises the following steps:
an objective function f (x, w) of the comprehensive energy system optimization scheduling model:
constraint conditions of the model:
H=H demand for ,C=C Demand for (18)
In the formula, P, T, I and Y respectively represent power, temperature, medium flow (current, heat flow, clean energy (natural gas, wind energy, solar energy and the like) flow and computational power flow) of each device in the comprehensive energy system and other important performance coefficients; the subscripts min and max represent the minimum and maximum values of the index, respectively; h and C represent total heat and total computing power generated by the system; h Demand for ,C Demand for Is a constant representing the user's thermal and computational demands of the dispatch.
As shown in fig. 5, in step S5, the comprehensive energy system optimized scheduling model is accessed to the mathematical model thereof to obtain an optimal scheduling method and result, and the specific method includes:
solving the optimal solution set of the model by adopting a non-dominated sorting genetic algorithm (NSGA-II) with an elite strategy to obtain a Pareto optimal solution set: f (x, w) = { f 1 (x,w),f 2 (x,w),...,f n (x, w) }, when y = { y = 1 ,y 2 ,...,y n Y is an operation parameter set of each system device in the scheduling process;
selecting a solution closest to the optimal scheme from the Pareto optimal solution set by adopting a TOPSIS method to obtain an optimal scheduling result, wherein the weights of the three target functions are w respectively I ,w II ,w III 。
Claims (7)
1. A distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method is characterized by comprising the following steps:
s1, establishing a comprehensive energy system mathematical model integrating distributed data central computing power and energy flow;
s2, establishing a comprehensive energy system operation evaluation index system comprising three primary indexes of economy, safety and cleanness and a plurality of secondary indexes;
s3, determining the comprehensive weight of each index by adopting a comprehensive evaluation method;
s4, constructing an optimized dispatching model of the comprehensive energy system by taking the lowest operation cost, the highest safety and the lowest pollution emission as three objective functions;
s5, accessing the optimized scheduling model of the comprehensive energy system into the mathematical model of the comprehensive energy system to obtain an optimal scheduling method and a result;
the step S1 specifically includes: for a single energy conversion device, the relationship between the input IN and the output OUT is:
OUT=μ·IN (1)
where μ represents a coupling coefficient between the input and the output;
for an integrated energy system comprising a plurality of energy transforming devices and a plurality of forms of energy, the relationship between input and output may be described by a matrix:
in the formula, subscripts P, G, H, C represent electric power, clean energy, thermal power, and computational power, respectively; mu.s PG The coupling coefficient between the input and the output when the electric power is used as the input and the clean energy is used as the output; mu.s GP When the clean energy is used as input and the electric power is used as output, the coupling coefficient between the input and the output is the coupling coefficient;
the step S4 specifically comprises the following steps:
an objective function f (x, w) of the comprehensive energy system optimization scheduling model:
constraint conditions of the model:
H=H demand for ,C=C Demand for (5)
In the formula, P, T, I and Y respectively represent the power, temperature, medium flow and bearable pressure of each device in the comprehensive energy system; the subscripts min and max represent the minimum and maximum values of the index, respectively; h and C represent total heat and total computing power generated by the system; h Demand for ,C Demand for Is a constant representing the user's thermal and computational demands of the dispatch.
2. The distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method according to claim 1, wherein the step S1 specifically comprises:
the method comprises the steps of modeling the comprehensive energy system based on graph theory, abstracting the connection position of equipment into connection nodes, abstracting lines between the two nodes into edges, simplifying each equipment in the system into a physical model according to actual system data, and establishing a topological structure of the comprehensive energy system.
3. The distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method according to claim 1, wherein in step S2, of the three primary indexes:
the economic indicator I comprises energy material cost I 1 Equipment operation maintenance cost I 2 Energy loss cost I 3 Three secondary indexes;
safety index II includes equipment failure rate II 1 And the fatigue cracking probability II of the heat supply network pipeline 2 Corrosion failure rate II 3 Three secondary indexes;
cleanliness indices III include solid waste discharge III 1 Clean energy utilization rate III 2 Greenhouse gas emission III 3 Three secondary indexes.
4. The distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method according to claim 3, wherein the step S3 specifically comprises: obtaining objective weight of each index by entropy weight methodObtaining subjective weight of each index by using analytic hierarchy processThen its composite weight isi=I,I 1 ,I 2 ,I 3 ,II,II 1 ,II 2 ,II 3 ,III,III 1 ,III 2 ,III 3 。
5. The distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method according to claim 4,
the economic indicator I represents a function x (I) as follows:
the safety index II represents a function x (II) of:
the cleanliness index III represents a function x (III) of:
in the formulae (6), (7) and (8), w i The weight of each index is represented, and x (i) represents each influence index.
6. The method for optimizing and scheduling a distributed data center computing power and energy flow integrated comprehensive energy system according to claim 5, wherein the influence indexes are as follows:
x(I 1 )、x(I 2 ) And x (I) 3 ) Respectively as follows:
in the formula (I), the compound is shown in the specification,respectively representing the prices of electric power and clean energy; z P ,Z G Representing the amount of purchased electric power and clean energy; gamma ray P ,γ G ,γ H ,γ C Respectively representing the unit cost of operation and maintenance of electric power equipment, clean energy equipment, thermal power equipment and power calculation equipment, and converting the unit cost into unit power; p P ,P G ,P H ,P C Respectively representing the working power of electric equipment, clean energy equipment, thermal equipment and force calculating equipment; delta phi P ,ΔΦ G ,ΔΦ H ,ΔΦ C Respectively representing power loss, clean energy loss, heat loss and computational power loss in the operation process;
x(II 1 )、x(II 2 ) And x (II) 3 ) Respectively as follows:
in the formula, T represents the running time, p represents the bearing pressure of the equipment, T represents the medium temperature, alpha represents each performance coefficient of the equipment, beta represents the external environment factor, and lambda represents the artificial destruction probability;
x(III 1 )、x(III 2 ) And x (III) 3 ) Respectively as follows:
wherein epsilon represents the degree of contamination of the solid waste; s Solid body ,S Cleaning of ,S General assembly Respectively representing the solid waste discharge amount, the clean energy utilization amount and the total energy utilization amount; p Total yield The total generated heat and calculated force of the system are expressed and converted into electric power; l is Greenhouse gases Indicating greenhouse gas emissions.
7. The distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method according to claim 1, wherein the step S5 specifically comprises:
adopting a non-dominated sorting genetic algorithm NSGA-II with an elite strategy to solve the optimal solution set of the model to obtain a Pareto optimal solution set: f (x, w) = { f 1 (x,w),f 2 (x,w),...,f n (x, w) }, when y = { y = 1 ,y 2 ,...,y n Y is an operation parameter set of each system device in the scheduling process;
selecting a solution closest to the optimal scheme from the Pareto optimal solution set by adopting a TOPSIS method to obtain an optimal scheduling result, wherein the weights of the three target functions are w respectively I ,w II ,w III 。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011102939.0A CN112232984B (en) | 2020-10-15 | 2020-10-15 | Distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011102939.0A CN112232984B (en) | 2020-10-15 | 2020-10-15 | Distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112232984A CN112232984A (en) | 2021-01-15 |
CN112232984B true CN112232984B (en) | 2023-02-17 |
Family
ID=74113744
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011102939.0A Active CN112232984B (en) | 2020-10-15 | 2020-10-15 | Distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112232984B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113158486B (en) * | 2021-05-08 | 2023-02-24 | 国网河北省电力有限公司经济技术研究院 | User side comprehensive energy network optimization control method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110363353A (en) * | 2019-07-16 | 2019-10-22 | 厦门大学 | The optimization design and dispatching method and system of a kind of Distributed Integration energy resource system |
CN110417063A (en) * | 2019-07-31 | 2019-11-05 | 广东电网有限责任公司 | A kind of integrated energy system Optimization Scheduling a few days ago |
CN110571863A (en) * | 2019-08-06 | 2019-12-13 | 国网山东省电力公司经济技术研究院 | Distributed power supply maximum acceptance capacity evaluation method considering flexibility of power distribution network |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110029142A1 (en) * | 2010-07-02 | 2011-02-03 | David Sun | System tools that provides dispatchers in power grid control centers with a capability to make changes |
-
2020
- 2020-10-15 CN CN202011102939.0A patent/CN112232984B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110363353A (en) * | 2019-07-16 | 2019-10-22 | 厦门大学 | The optimization design and dispatching method and system of a kind of Distributed Integration energy resource system |
CN110417063A (en) * | 2019-07-31 | 2019-11-05 | 广东电网有限责任公司 | A kind of integrated energy system Optimization Scheduling a few days ago |
CN110571863A (en) * | 2019-08-06 | 2019-12-13 | 国网山东省电力公司经济技术研究院 | Distributed power supply maximum acceptance capacity evaluation method considering flexibility of power distribution network |
Non-Patent Citations (5)
Title |
---|
《Energy Conscious Scheduling for Distributed Computing Systems under Different Operating Conditions》;Young Choon Lee et al.;《IEEE Transactions on Parallel and Distributed Systems》;20101230;1374-1381 * |
《Sizing domestic air-source heat pump systems with thermal storage under varying electrical load shifting strategies》;Marini D et al.;《Applied Energy》;20190918;1-14 * |
《基于区块链的能源电力供需网调度优化模型》;胡伟 等;《系统管理学》;20191206;1134-1142 * |
基于区块链的虚拟电厂模型研究;邵炜晖等;《计算机科学》;20180215(第02期);32-38 * |
考虑综合评估的孤立微网优化调度;郭宁等;《电力建设》;20160401(第04期);67-73 * |
Also Published As
Publication number | Publication date |
---|---|
CN112232984A (en) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Stochastic-weighted robust optimization based bilayer operation of a multi-energy building microgrid considering practical thermal loads and battery degradation | |
Zheng et al. | Coordinated scheduling strategy to optimize conflicting benefits for daily operation of integrated electricity and gas networks | |
CN107301472A (en) | Distributed photovoltaic planing method based on scene analysis method and voltage-regulation strategy | |
WO2019237316A1 (en) | Knowledge-transfer-based modeling method for blast furnace coal gas scheduling system | |
CN109193628B (en) | Comprehensive energy system energy management method based on consistency | |
Liu et al. | Integrated optimal dispatching strategy considering power generation and consumption interaction | |
CN112701687A (en) | Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response | |
CN106786801A (en) | One kind is based on micro-capacitance sensor operation method in a balanced way of bidding | |
CN114638502A (en) | Power distribution network line reinforcement strategy considering flexibility of demand side resources | |
CN112348276A (en) | Comprehensive energy system planning optimization method based on multiple elements and three levels | |
CN115860413A (en) | Grid-connected combined cooling heating and power micro-grid system economic scheduling method based on load demand response and double-layer adjustable robust optimization | |
CN113779792A (en) | Affine-based comprehensive energy system optimal configuration method | |
CN117236629A (en) | Hierarchical cooperative scheduling method for electric-carbon coupled multi-energy system | |
CN112418488A (en) | Comprehensive energy system scheduling method and device based on two-stage energy optimization | |
CN112232984B (en) | Distributed data center computing power and energy flow fused comprehensive energy system optimization scheduling method | |
Li et al. | Intraday multi-objective hierarchical coordinated operation of a multi-energy system | |
Nappu et al. | Economic emission dispatch for thermal power plant in Indonesia | |
Jiang et al. | Scenario-transformation-based optimal sizing of hybrid hydrogen-battery storage for multi-timescale islanded microgrids | |
CN116362421B (en) | Energy supply distribution prediction system and method based on comprehensive overall analysis of energy sources | |
CN115860241B (en) | SALR-based distributed scheduling method and device for interconnected multi-region comprehensive energy system | |
CN116187538B (en) | Energy scheduling method and device | |
CN112564182A (en) | Distributed robust opportunity constraint-based unified scheduling method for multi-energy combined system | |
Qiu et al. | Local integrated energy system operational optimization considering multi‐type uncertainties: A reinforcement learning approach based on improved TD3 algorithm | |
Liu et al. | Influence Evaluation of Integrated Energy System on the Unit Commitment in Power System | |
CN109472716A (en) | A kind of production capacity end power distribution method multifactor based on the considerations of genetic algorithm |
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 |