CN110929980A - Planning method and device for regional comprehensive energy system - Google Patents

Planning method and device for regional comprehensive energy system Download PDF

Info

Publication number
CN110929980A
CN110929980A CN201910958835.0A CN201910958835A CN110929980A CN 110929980 A CN110929980 A CN 110929980A CN 201910958835 A CN201910958835 A CN 201910958835A CN 110929980 A CN110929980 A CN 110929980A
Authority
CN
China
Prior art keywords
energy
load center
energy station
station
load
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
CN201910958835.0A
Other languages
Chinese (zh)
Other versions
CN110929980B (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.)
Tsinghua University
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
Tsinghua University
Electric Power Research Institute of State Grid Shandong Electric Power 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 Tsinghua University, Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN201910958835.0A priority Critical patent/CN110929980B/en
Publication of CN110929980A publication Critical patent/CN110929980A/en
Application granted granted Critical
Publication of CN110929980B publication Critical patent/CN110929980B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • 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

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a regional integrated energy system planning method and a device, wherein the method comprises the steps of calculating a regional integrated energy system partition load center, determining the grade and the quantity of energy stations, calculating the coordinates of the energy stations, matching the energy stations with the load center, and calculating the service capacity of the energy stations. According to the invention, the energy stations are divided into different levels, and the energy stations are planned according to the service radius, the service capacity and the service requirement of the load center and the lowest overall construction cost, so that the number, the service scale and the positions of the energy stations are obtained, the optimal layout of a plurality of distributed energy stations is realized, and the construction requirement of the distributed energy stations according to local conditions is met.

Description

Planning method and device for regional comprehensive energy system
Technical Field
The invention relates to the technical field of energy planning, in particular to a planning method and a planning device of a regional comprehensive energy system comprising a plurality of distributed energy stations.
Background
The existing energy system planning is usually based on the supply and demand relationship, and the special planning of each energy source is respectively carried out. In the actual process, the special plans are influenced interactively, so that the problems of repeated load calculation, occupied land conflict, unreasonable energy structure and the like can be caused. In energy internet planning, planning of all energy systems is integrated, and how to realize optimal matching of energy supply and demand under interaction influence and constraint relation of all energy is a key problem of comprehensive energy system planning and design.
The existing planning work usually adopts a top-down method, namely: the method comprises the steps of planning the generation and transmission of main energy in a large range, and planning the energy supply in a local area. For an integrated energy system with interactive influence of various energy sources, the method is difficult to carry out.
In order to solve the problem of difficult top-down planning in energy systems, a new planning method is urgently needed.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a method and an apparatus for planning a regional integrated energy system, so as to solve the problem that the current energy system planning cannot meet the distributed energy supply requirement of the regional integrated energy system.
The invention provides a planning method of a regional comprehensive energy system, which comprises the following steps:
s110: establishing load partitions of the regional comprehensive energy system according to the planned regional plot division, establishing centers of all the load partitions as load centers, and calculating energy demand of all the load centers;
s120: determining the number and the grade of the energy stations according to the energy demand of each load center, and determining the position of each energy station according to the position of each load center;
s130: judging whether a feasible energy pipeline path exists from the load center to each energy station, if so, further judging whether the energy pipeline path meets the requirement of the service radius of one energy station, and if so, selecting the energy station to serve the load center;
s140: if the requirement of the service radius is not met, judging whether the requirement of the service radius of other energy stations is met, and selecting a proper energy station to serve the load center according to the requirement of the service radius and the requirement of the maximum service capacity of the energy station;
s150: for the load centers which are not matched with the proper energy stations, the steps S120 to S140 are circulated until all the load centers are matched with the proper energy stations;
s160: and after each load center is matched with a proper energy station, calculating the required service capacity of each energy station.
In addition, it is preferable that the calculation formula of the energy demand of the load center is:
Qj＝εejQej+εhjQhj+εcjQcj
wherein Q isjEnergy demand for jth load center;
Qej、Qhj、Qcjthe annual power consumption, the heat consumption and the cold consumption of the load center are respectively;
εej、εhj、εcjrespectively, the unified conversion coefficients of electricity, heat and cold.
In addition, it is preferable that in the step S120, the service-by-service is performedThe capacity divides the energy station into A, B, C three energy station grades, and the unit cost of each grade energy station is cA、cB、cCFixed construction cost αA、αB、αC(ii) a The service radius of each grade of energy station is rA、rB,、rCGiven an initial number s of three types of energy stationsA＝5、s＝8、sC＝8；
Wherein, in the process of determining the number and the type of the energy stations according to the energy demand of each load center,
the first step is as follows: defining an association decision matrix DM of the load center and the energy station, the matrix DM being(s)A+sB+sC) The x n dimension, a matrix element of 1 indicates that the load center is to be served by the corresponding energy station, and a matrix element of 0 indicates that the load center is not to be served by the corresponding energy station; wherein the content of the first and second substances,
the sum of each column of the matrix DM must be 1, and the decision variables are the elements in the matrix DM: dmij，dmij∈{0,1}，i＝1...m，j＝1...n，m＝sA+sB+sC；
The second step is that: establishing an objective function, wherein the objective function is that the total cost of the energy station is the least, and the specific formula is as follows:
min(C·DM·Q+α·D)
wherein the content of the first and second substances,
Figure BDA0002228257990000031
Q＝[Q1,Q2,...Qn]T，
Figure BDA0002228257990000032
D＝[d1,d2,...dm]T，
Figure BDA0002228257990000033
and (3) establishing an equality constraint:
Figure BDA0002228257990000034
n, the equality constraint indicates that a certain load center has and only one energy station to serve it;
establishing inequality constraints:
Figure BDA0002228257990000035
1, wherein SmaxiMaximizing building service capacity for energy stations
The third step: and solving the problem established in the first step by using branch-and-bound routing, and solving an association decision matrix of the load center and the energy station and A, B, C the number of three energy station grades.
In addition, it is preferable that the requirement formula of the service capability of each energy station is as follows:
Figure BDA0002228257990000036
wherein Q isjRepresenting the energy demand of the jth load center;
dmijindicating the association of the jth load center with the ith energy station,
when dm ij1 means that the jth load center is served by the ith energy station,
when dm ij0 means that the jth load center is not associated with the 2 nd energy station.
In addition, it is preferable that the position of each energy station is determined by using the formula that the coordinate (X) of the ith energy stationi,Yi) Comprises the following steps:
Figure BDA0002228257990000041
wherein the content of the first and second substances,
Figure BDA0002228257990000042
representing the service capability requirements of the energy station.
The invention also provides a planning device of the regional comprehensive energy system, which comprises:
the load center energy demand acquisition unit is used for establishing load partitions of the regional comprehensive energy system according to the planned regional plot division, establishing the centers of the load partitions as load centers, and calculating the energy demand of the load centers;
the energy station quantity and type determining unit is used for determining the quantity and the grade of the energy stations according to the energy demand of each load center;
the energy station position acquisition unit is used for determining the position of each energy station according to the position of each load center;
the matching unit of the load center and the energy stations is used for judging whether a feasible energy pipeline path exists from the load center to each energy station or not, further judging whether the energy pipeline path meets the requirement of the service radius of one energy station or not if the feasible energy pipeline path exists, and selecting the energy stations to serve the load center if the feasible energy pipeline path meets the requirement of the service radius;
if the requirement of the service radius is not met, judging whether the requirement of the service radius of other energy stations is met, and selecting a proper energy station to serve the load center according to the requirement of the service radius and the requirement of the maximum service capacity of the energy station;
for the load centers which are not matched with the proper energy stations, the matching mode is circulated until each load center is matched with the proper energy station;
and the energy station service capacity calculating unit is used for calculating the required service capacity of each energy station.
In addition, preferably, the calculation formula of the energy demand of the load center energy demand acquisition unit is as follows:
Qj＝εejQej+εhjQhj+εcjQcj
wherein Q isjEnergy demand for jth load center;
Qej、Qhj、Qcjthe annual power consumption, the heat consumption and the cold consumption of the load center are respectively;
εej、εhj、εcjrespectively, the unified conversion coefficients of electricity, heat and cold.
Furthermore, it is preferable that the energy station number and type determination unit divides the energy stations into A, B, C three energy station classes according to service capacity, and the cost per unit of each class of energy station is cA、cB、cCFixed construction cost αA、αB、αC(ii) a The service radius of each grade of energy station is rA、rB,、rCGiven an initial number s of three types of energy stationsA＝5、s＝8、sC＝8；
Wherein, in the process of determining the number and the type of the energy stations according to the energy demand of each load center,
the first step is as follows: defining an association decision matrix DM of the load center and the energy station, the matrix DM being(s)A+sB+sC) The x n dimension, a matrix element of 1 indicates that the load center is to be served by the corresponding energy station, and a matrix element of 0 indicates that the load center is not to be served by the corresponding energy station; wherein the content of the first and second substances,
the sum of each column of the matrix DM must be 1, and the decision variables are the elements in the matrix DM: dmij，dmij∈{0,1}，i＝1...m，j＝1...n，m＝sA+sB+sC；
The second step is that: establishing an objective function, wherein the objective function is that the total cost of the energy station is the least, and the specific formula is as follows:
min(C·DM·Q+α·D)
wherein the content of the first and second substances,
Figure BDA0002228257990000051
Q＝[Q1,Q2,...Qn]T，
Figure BDA0002228257990000052
D＝[d1,d2,...dm]T，
Figure BDA0002228257990000053
the third step: and solving the problem established in the first step by using branch-and-bound routing, and solving an association decision matrix of the load center and the energy station and A, B, C the number of three energy station grades.
In addition, it is preferable that the requirement formula of the service capability of each energy station of the energy station service capability calculation unit is as follows:
Figure BDA0002228257990000061
wherein Q isjRepresenting the energy demand of the jth load center;
dmijindicating the association of the jth load center with the ith energy station,
when dm ij1 means that the jth load center is served by the ith energy station,
when dm ij0 means that the jth load center is not associated with the 2 nd energy station.
In addition, it is preferable that the energy station position acquiring unit determines the position of each energy station by using the formula that the coordinate (X) of the ith energy stationi,Yi) Comprises the following steps:
Figure BDA0002228257990000062
wherein the content of the first and second substances,
Figure BDA0002228257990000063
representing the service capability requirements of the energy station.
According to the technical scheme, the planning method and the planning device of the regional comprehensive energy system divide the energy stations into different levels, plan the energy stations according to the service radius, the service capacity and the service requirement of the load center and the lowest overall construction cost, and obtain the number, the service scale and the positions of the energy stations. The planning method breaks through the limitation that network constraint or capacity constraint is only considered in site selection and volume fixing of the comprehensive energy station, takes the energy station close to the load center as a principle, realizes the economic target of energy system construction, ensures that the energy system planning scheme has more reasonability and adaptability, ensures reliable energy supply of the load center, improves the utilization rate of energy station facilities, and reduces the investment of the energy system.
To the accomplishment of the foregoing and related ends, one or more aspects of the invention comprise the features hereinafter fully described. The following description and the annexed drawings set forth in detail certain illustrative aspects of the invention. These aspects are indicative, however, of but a few of the various ways in which the principles of the invention may be employed. Further, the present invention is intended to include all such aspects and their equivalents.
Drawings
Other objects and results of the present invention will become more apparent and more readily appreciated as the same becomes better understood by reference to the following description taken in conjunction with the accompanying drawings. In the drawings:
fig. 1 is a schematic flow chart of a planning method of a regional integrated energy system according to an embodiment of the present invention;
FIG. 2 is a schematic view of a load center distribution according to an embodiment of the present invention;
fig. 3 is a diagram illustrating a result of energy station planning according to an embodiment of the present invention;
fig. 4 is a schematic diagram of a planning apparatus of a regional integrated energy system according to an embodiment of the present invention.
The same reference numbers in all figures indicate similar or corresponding features or functions.
Detailed Description
In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of one or more embodiments. It may be evident, however, that such embodiment(s) may be practiced without these specific details.
At present, the problem that the planning of an energy system can not meet the distributed energy supply requirement of a regional comprehensive energy system is solved, and a bottom-up planning method is adopted for the regional comprehensive energy system comprising a plurality of energy stations. For the regional integrated energy system as an example, the user-oriented differentiated requirements should be combined with the local resource conditions of the user to perform the optimal design of each local energy supply system, on the basis, the energy interaction requirements of each local system and the outside are determined, and then the regional centralized energy supply device and the network are optimally configured and distributed. The invention provides an optimized planning method of a regional integrated energy system with a plurality of energy stations aiming at the distributed energy supply requirement of the regional integrated energy system.
Compared with the current comprehensive energy system planning method, the regional comprehensive energy system planning method realizes the optimized layout of a plurality of distributed energy stations and meets the construction requirements of the distributed energy stations according to local conditions.
Specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
In order to illustrate the method for planning the regional integrated energy system provided by the present invention, fig. 1 shows a flow of the method for planning the regional integrated energy system according to an embodiment of the present invention.
As shown in fig. 1, the method for planning the regional integrated energy system according to the present invention includes:
s110: establishing load partitions of the regional comprehensive energy system according to the planned regional plot division, establishing centers of all the load partitions as load centers, and calculating energy demand of all the load centers;
s120: determining the number and the grade of the energy stations according to the energy demand of each load center, and determining the position of each energy station according to the position of each load center;
s130: judging whether a feasible energy pipeline path exists from the load center to each energy station, if so, further judging whether the energy pipeline path meets the requirement of the service radius of one energy station, and if so, selecting the energy station to serve the load center;
s140: if the requirement of the service radius is not met, judging whether the requirement of the service radius of other energy stations is met, and selecting a proper energy station to serve as the load center according to the requirement of the service radius and the requirement of the maximum service capacity of the energy station;
s150: for the load centers which are not matched with the proper energy stations, the steps S120 to S140 are circulated until all the load centers are matched with the proper energy stations;
s160: and after each load center is matched with a proper energy station, calculating the required service capacity of each energy station.
In step S110, the regional integrated energy system partition load center is calculated. Establishing regional comprehensive energy system load partitions according to planned regional plot division, and counting the electricity, heat and cold loads of each partition in a planned year; the center of each plot was established as the load center, and a load center profile was plotted.
Wherein, the energy demand of the jth load center is set as QjThe position coordinate of which is (x)j,yj)；QjThe calculation method is Q for the comprehensive demand of cooling, heating and power of the load centerj＝εejQej+εhjQhj+εcjQcjWherein Q isej、Qhj、QcjAnnual electricity, heat and cold consumption, epsilon, of load centersej、εhj、εcjRespectively, the unified conversion coefficients of electricity, heat and cold.
In an embodiment of the present invention, take εej＝1，εhj＝0.35、εcj0.25. The loads of the plots in the examples are shown in the following table:
Figure BDA0002228257990000081
Figure BDA0002228257990000091
in step S120, the rank and number of energy stations are determined. In the embodiment of the invention, the energy station is divided into A, B, C energy station grades, and the unit cost of each grade energy station is cA、cB、cCFixed construction cost αA、αB、αC(ii) a The service radius of each grade of energy station is rA、rB,、rCGiven an initial number s of three types of energy stationsA＝5、s＝8、sC8. The specific method comprises the following steps:
1) defining an associated decision moment DM between the load center and the energy station, the matrix being(s)A+sB+sC) The xn dimension, a matrix element of 1 indicates that the load center is to be served by the corresponding energy station, and an element of 0 indicates that the load center is not to be served by the corresponding energy station. The sum of each column of the matrix must be 1, i.e. there is one and only one interval to supply it to any one load. The decision variables are the elements in the matrix DM: dmij，dmij∈{0,1}，i＝1...m，j＝1...n，m＝sA+sB+sC；
2) Establishing an objective function, wherein the objective function is that the total cost of the energy station is the least, and the specific formula is as follows:
min(C·DM·Q+α·D)
wherein the content of the first and second substances,
Figure BDA0002228257990000101
Q＝[Q1,Q2,...Qn]T，
Figure BDA0002228257990000102
D＝[d1,d2,...dm]T，
Figure BDA0002228257990000103
and (3) establishing an equality constraint:
Figure BDA0002228257990000104
n, which indicates that there is one and only one energy station serving a certain load center.
Establishing inequality constraints:
Figure BDA0002228257990000105
1, wherein SmaxiAnd the maximum construction service capacity of the energy station is provided.
3) Solving a binary Integer Linear Programming (ILP) problem formed by the established objective function constraints, and solving an association decision matrix DM of the load center and the energy station, and A, B, C the number of three energy station levels. In this embodiment, the preliminarily selected energy station includes: 2A-level energy stations, 3B-level energy stations and 3C-level energy stations.
4) And calculating coordinates of each energy station. Theoretical site coordinates (X) of the ith energy station on the principle that the energy station is as close as possible to the load centeri,Yi) Comprises the following steps:
Figure BDA0002228257990000106
in step S130, it is determined whether there is a feasible energy pipeline path from the load center to each energy station, and the service radius of the energy station is satisfied, and if the radius requirement is satisfied, the energy station is selected to serve. In this embodiment, the determinable energy stations are shown in the following table:
Figure BDA0002228257990000111
in step S140, for the load center that does not satisfy the service radius requirement, it is determined whether the service radius of the other energy station is satisfied, and the service capability constraint of the energy station is satisfied, and the nearest energy station that satisfies the above condition is found to be accessed. The determination of the energy station in this embodiment is shown in the following table:
load center serial number Energy station Load center serial number Energy station
1 Class a energy station 1 16 C-class energy station 1
2 Class a energy station 1 17
3 Class a energy station 1 18
4 Class a energy station 1 19 C-class energy station 1
5 20 Class B energy station 4
6 21 Class A energy station 2
7 Class B energy station 1 22 Class A energy station 2
8 Class B energy station 1 23
9 Class B energy station 3 24 Class A energy station 2
10 Class B energy station 3 25 C-class energy station 2
11 26
12 Class a energy station 1 27 Class B energy station 2
13 Class a energy station 1 28 Class B energy station 2
14 C-class energy station 1 29 Class B energy station 2
15 C-class energy station 3 30 Class B energy station 2
In step S150, for the load center with no available energy station, the process returns to step S120 to plan a new energy station to serve. In an embodiment of the invention, the load centers 5, 6, 11, 17, 18, 20, 23, 26 do not obtain an energy station plan. Returning to the steps S120, S130, and S140 to plan a new energy station for the load center, where the planning result is: the load centers 5, 6, 17, 11, 20 plan a class a energy station (labeled class a energy station 3), the load centers 18, 23 plan a class B energy station (labeled class B energy station 4), and the load center 26 plans a class C energy station (labeled class C energy station 4). To this end, all load centers are planned with energy stations.
In step S160, the service capability requirement of each energy station is calculated:
Figure BDA0002228257990000121
wherein Q isjRepresenting the energy demand of the jth load center; dmijShowing the association relationship between the jth load center and the ith energy station when dm ij1 means that the jth load center is driven by the ith energy sourceStation service, when dm ij0 means that the jth load center is not associated with the 2 nd energy station.
According to the planning result, the service capability requirement of each energy station in this embodiment is as follows:
Figure BDA0002228257990000122
Figure BDA0002228257990000131
as shown in fig. 3, the positions of the energy stations are plotted on the distribution diagram of the load centers, each energy station is represented by "■", and the service relationship between the energy station and each load center is represented by "- - -".
Corresponding to the method, the invention also provides a planning device of the regional integrated energy system, and fig. 4 shows a logical structure of the planning device of the regional integrated energy system according to the embodiment of the invention.
As shown in fig. 4, the planning apparatus 400 for regional integrated energy system according to the present invention includes:
a load center energy demand obtaining unit 410, configured to establish load partitions of the regional integrated energy system according to the planned regional parcel division, establish the centers of the load partitions as load centers, and calculate energy demand of each load center;
the energy station number and type determining unit 420 is used for determining the number and the level of the energy stations according to the energy demand of each load center;
an energy station position obtaining unit 430, configured to determine positions of the energy stations according to positions of the load centers;
a load center and energy station matching unit 440, configured to determine whether there is a feasible energy pipeline path from the load center to each energy station, further determine whether the energy pipeline path meets a requirement of a service radius of one of the energy stations if there is a feasible energy pipeline path, and select the energy station to serve the load center if the requirement of the service radius is met;
if the requirement of the service radius is not met, judging whether the requirement of the service radius of other energy stations is met, and selecting a proper energy station to serve the load center according to the requirement of the service radius and the requirement of the maximum service capacity of the energy station;
for the load centers which are not matched with the proper energy stations, the matching mode is circulated until each load center is matched with the proper energy station;
the energy station service capability calculating unit 450 is configured to calculate a required service capability of each energy station.
The calculation formula of the energy demand of the load center energy demand acquisition unit load center is as follows:
Qj＝εejQej+εhjQhj+εcjQcj
wherein Q isjEnergy demand for jth load center;
Qej、Qhj、Qcjthe annual power consumption, the heat consumption and the cold consumption of the load center are respectively;
εej、εhj、εcjrespectively, the unified conversion coefficients of electricity, heat and cold.
The energy station quantity and type determining unit divides the energy stations into A, B, C three energy station levels according to service capacity, and the unit cost of each level of energy stations is cA、cB、cCFixed construction cost αA、αB、αC(ii) a The service radius of each grade of energy station is rA、rB,、rCGiven an initial number s of three types of energy stationsA＝5、s＝8、sC＝8；
Wherein, in the process of determining the number and the type of the energy stations according to the energy demand of each load center,
the first step is as follows: defining an association decision matrix DM of said load center and said energy station, said matrixDM is(s)A+sB+sC) The x n dimension, a matrix element of 1 indicates that the load center is to be served by the corresponding energy station, and a matrix element of 0 indicates that the load center is not to be served by the corresponding energy station; wherein the content of the first and second substances,
the sum of each column of the matrix DM must be 1, and the decision variables are the elements in the matrix DM: dmij，dmij∈{0,1}，i＝1...m，j＝1...n，m＝sA+sB+sC；
The second step is that: establishing an objective function, wherein the objective function is that the total cost of the energy station is the least, and the specific formula is as follows:
min(C·DM·Q+α·D)
wherein the content of the first and second substances,
Figure BDA0002228257990000141
Q＝[Q1,Q2,...Qn]T，
Figure BDA0002228257990000142
D＝[d1,d2,...dm]T，
Figure BDA0002228257990000143
and (3) establishing an equality constraint:
Figure BDA0002228257990000151
n, the equality constraint indicates that a certain load center has and only one energy station to serve it;
establishing inequality constraints:
Figure BDA0002228257990000152
1, wherein SmaxiThe maximum construction service capacity of the energy station is provided;
the third step: and solving the problem established in the first step by using branch-and-bound routing, and solving an association decision matrix of the load center and the energy station and A, B, C the number of three energy station grades.
The requirement formula of the service capability of each energy station of the energy station service capability calculation unit is as follows:
Figure BDA0002228257990000153
wherein Q isjRepresenting the energy demand of the jth load center;
dmijindicating the association of the jth load center with the ith energy station,
when dm ij1 means that the jth load center is served by the ith energy station,
when dm ij0 means that the jth load center is not associated with the 2 nd energy station.
Wherein the energy station position obtaining unit determines the position of each energy station by adopting the following formula, and the coordinate (X) of the ith energy stationi,Yi) Comprises the following steps:
Figure BDA0002228257990000154
wherein the content of the first and second substances,
Figure BDA0002228257990000155
representing the service capability requirements of the energy station.
According to the method and the device for planning the regional comprehensive energy system, the energy stations are divided into different levels, and the energy stations are planned according to the service radius, the service capacity and the service requirement of the load center and the total construction cost to obtain the number, the service scale and the position of the energy stations. The planning method breaks through the limitation that network constraint or capacity constraint is only considered in site selection and volume fixing of the comprehensive energy station, takes the energy station close to the load center as a principle, realizes the economic target of energy system construction, ensures that the energy system planning scheme has more reasonability and adaptability, ensures reliable energy supply of the load center, improves the utilization rate of energy station facilities, and reduces the investment of the energy system.
The planning method and apparatus of the regional integrated energy system according to the present invention are described above by way of example with reference to the accompanying drawings. However, it should be understood by those skilled in the art that various modifications can be made to the method and apparatus for planning a regional integrated energy system of the present invention without departing from the scope of the present invention. Therefore, the scope of the present invention should be determined by the contents of the appended claims.

Claims (10)

1. A method for planning a regional integrated energy system includes:
s110: establishing load partitions of the regional comprehensive energy system according to the planned regional plot division, establishing centers of all the load partitions as load centers, and calculating energy demand of all the load centers;
s120: determining the number and the grade of the energy stations according to the energy demand of each load center, and determining the position of each energy station according to the position of each load center;
s130: judging whether a feasible energy pipeline path exists from the load center to each energy station, if so, further judging whether the energy pipeline path meets the requirement of the service radius of one energy station, and if so, selecting the energy station to serve the load center;
s140: if the requirement of the service radius is not met, judging whether the requirement of the service radius of other energy stations is met, and selecting a proper energy station to serve the load center according to the requirement of the service radius and the requirement of the maximum service capacity of the energy station;
s150: for the load centers which are not matched with the proper energy stations, the steps S120 to S140 are circulated until all the load centers are matched with the proper energy stations;
s160: and after each load center is matched with a proper energy station, calculating the required service capacity of each energy station.
2. The method for planning the regional integrated energy system according to claim 1,
the calculation formula of the energy demand of the load center is as follows:
Qj＝εejQej+εhjQhj+εcjQcj
wherein Q isjEnergy demand for jth load center;
Qej、Qhj、Qcjthe annual power consumption, the heat consumption and the cold consumption of the load center are respectively;
εej、εhj、εcjrespectively, the unified conversion coefficients of electricity, heat and cold.
3. The method for planning the regional integrated energy system according to claim 1,
in the step S120, the energy station is divided into A, B, C three energy station levels according to service capacity, and the unit cost of each level energy station is cA、cB、cCFixed construction cost αA、αB、αC(ii) a The service radius of each grade of energy station is rA、rB,、rCGiven an initial number s of three types of energy stationsA＝5、s＝8、sC＝8；
Wherein, in the process of determining the number and the type of the energy stations according to the energy demand of each load center,
the first step is as follows: defining an association decision matrix DM of the load center and the energy station, the matrix DM being(s)A+sB+sC) The x n dimension, a matrix element of 1 indicates that the load center is to be served by the corresponding energy station, and a matrix element of 0 indicates that the load center is not to be served by the corresponding energy station; wherein the content of the first and second substances,
the sum of each column of the matrix DM must be 1, and the decision variables are the elements in the matrix DM: dmij，dmij∈{0,1}，i＝1...m，j＝1...n，m＝sA+sB+sC；
The second step is that: establishing an objective function, wherein the objective function is that the total cost of the energy station is the least, and the specific formula is as follows:
min(C·DM·Q+α·D)
wherein the content of the first and second substances,
Figure FDA0002228257980000021
Q＝[Q1,Q2,...Qn]T，
Figure FDA0002228257980000022
D＝[d1,d2,...dm]T，
Figure FDA0002228257980000023
and (3) establishing an equality constraint:
Figure FDA0002228257980000024
the equality constraint indicates that a certain load center has only one energy station serving it;
establishing inequality constraints:
Figure FDA0002228257980000025
wherein S ismaxiThe maximum construction service capacity of the energy station is provided;
the third step: and solving the problem established in the first step by using branch-and-bound routing, and solving an association decision matrix of the load center and the energy station and A, B, C the number of three energy station grades.
4. The method for planning the regional integrated energy system according to claim 3,
the requirement formula of the service capacity of each energy station is as follows:
Figure FDA0002228257980000031
wherein Q isjRepresenting the energy demand of the jth load center;
dmijindicating the association of the jth load center with the ith energy station,
when dmij1 means that the jth load center is served by the ith energy station,
when dmij0 means that the jth load center is not associated with the 2 nd energy station.
5. The method for planning the regional integrated energy system according to claim 4,
the position of each energy station is determined by the following formula, wherein the coordinate (X) of the ith energy stationi,Yi) Comprises the following steps:
Figure FDA0002228257980000032
wherein the content of the first and second substances,
Figure FDA0002228257980000033
representing the service capability requirements of the energy station.
6. An apparatus for planning a regional integrated energy system, comprising:
the load center energy demand acquisition unit is used for establishing load partitions of the regional comprehensive energy system according to the planned regional plot division, establishing the centers of the load partitions as load centers, and calculating the energy demand of the load centers;
the energy station quantity and type determining unit is used for determining the quantity and the grade of the energy stations according to the energy demand of each load center;
the energy station position acquisition unit is used for determining the position of each energy station according to the position of each load center;
the matching unit of the load center and the energy stations is used for judging whether a feasible energy pipeline path exists from the load center to each energy station or not, further judging whether the energy pipeline path meets the requirement of the service radius of one energy station or not if the feasible energy pipeline path exists, and selecting the energy stations to serve the load center if the feasible energy pipeline path meets the requirement of the service radius;
if the requirement of the service radius is not met, judging whether the requirement of the service radius of other energy stations is met, and selecting a proper energy station to serve the load center according to the requirement of the service radius and the requirement of the maximum service capacity of the energy station;
for the load centers which are not matched with the proper energy stations, the matching mode is circulated until each load center is matched with the proper energy station;
and the energy station service capacity calculating unit is used for calculating the required service capacity of each energy station.
7. The regional integrated energy system planning apparatus of claim 6,
the calculation formula of the energy demand of the load center energy demand acquisition unit load center is as follows:
Qj＝εejQej+εhjQhj+εcjQcj
wherein Q isjEnergy demand for jth load center;
Qej、Qhj、Qcjthe annual power consumption, the heat consumption and the cold consumption of the load center are respectively;
εej、εhj、εcjrespectively, the unified conversion coefficients of electricity, heat and cold.
8. The regional integrated energy system planning apparatus of claim 6,
the energy station quantity and type determining unit divides the energy stations into A, B, C three energy station levels according to service capacity, and the unit cost of each level of energy stations is cA、cB、cCFixed buildingThe manufacturing cost is αA、αB、αC(ii) a The service radius of each grade of energy station is rA、rB,、rCGiven an initial number s of three types of energy stationsA＝5、s＝8、sC＝8；
Wherein, in the process of determining the number and the type of the energy stations according to the energy demand of each load center,
the first step is as follows: defining an association decision matrix DM of the load center and the energy station, the matrix DM being(s)A+sB+sC) The x n dimension, a matrix element of 1 indicates that the load center is to be served by the corresponding energy station, and a matrix element of 0 indicates that the load center is not to be served by the corresponding energy station; wherein the content of the first and second substances,
the sum of each column of the matrix DM must be 1, and the decision variables are the elements in the matrix DM: dmij，dmij∈{0,1}，i＝1...m，j＝1...n，m＝sA+sB+sC；
The second step is that: establishing an objective function, wherein the objective function is that the total cost of the energy station is the least, and the specific formula is as follows:
min(C·DM·Q+α·D)
wherein the content of the first and second substances,
Figure FDA0002228257980000051
Q＝[Q1,Q2,...Qn]T，
Figure FDA0002228257980000052
D＝[d1,d2,...dm]T，
Figure FDA0002228257980000053
and (3) establishing an equality constraint:
Figure FDA0002228257980000054
the equality constraint indicates that a certain load center has and only hasOne energy station serves it;
establishing inequality constraints:
Figure FDA0002228257980000055
wherein S ismaxiThe maximum construction service capacity of the energy station is provided;
the third step: and solving the problem established in the first step by using branch-and-bound routing, and solving an association decision matrix of the load center and the energy station and A, B, C the number of three energy station grades.
9. The regional integrated energy system planning apparatus of claim 8,
the requirement formula of the service capability of each energy station of the energy station service capability calculation unit is as follows:
Figure FDA0002228257980000056
wherein Q isjRepresenting the energy demand of the jth load center;
dmijindicating the association of the jth load center with the ith energy station,
when dmij1 means that the jth load center is served by the ith energy station,
when dmij0 means that the jth load center is not associated with the 2 nd energy station.
10. The regional integrated energy system planning apparatus of claim 9,
the energy station position acquisition unit determines the position of each energy station by adopting the following formula, namely the coordinate (X) of the ith energy stationi,Yi) Comprises the following steps:
Figure FDA0002228257980000061
wherein the content of the first and second substances,
Figure FDA0002228257980000062
representing the service capability requirements of the energy station.
CN201910958835.0A 2019-10-10 2019-10-10 Planning method and device for regional comprehensive energy system Active CN110929980B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910958835.0A CN110929980B (en) 2019-10-10 2019-10-10 Planning method and device for regional comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910958835.0A CN110929980B (en) 2019-10-10 2019-10-10 Planning method and device for regional comprehensive energy system

Publications (2)

Publication Number Publication Date
CN110929980A true CN110929980A (en) 2020-03-27
CN110929980B CN110929980B (en) 2022-06-24

Family

ID=69848819

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910958835.0A Active CN110929980B (en) 2019-10-10 2019-10-10 Planning method and device for regional comprehensive energy system

Country Status (1)

Country Link
CN (1) CN110929980B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112700175A (en) * 2021-01-21 2021-04-23 潜能恒信能源技术股份有限公司 Construction method and system of underground energy and overground demand coupled distributed energy station

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030105868A1 (en) * 2001-12-04 2003-06-05 Kimbrel Tracy J. Dynamic resource allocation using known future benefits
CN102751722A (en) * 2012-06-18 2012-10-24 上海交通大学 Grid network optimization method based on shortest feasible path and historical flow data
CN102880921A (en) * 2012-10-16 2013-01-16 山东电力集团公司电力科学研究院 Method for optimizing site selection of electric vehicle charging stations
CN110046791A (en) * 2019-02-28 2019-07-23 中国电力科学研究院有限公司 A kind of energy source station scene building method and system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030105868A1 (en) * 2001-12-04 2003-06-05 Kimbrel Tracy J. Dynamic resource allocation using known future benefits
CN102751722A (en) * 2012-06-18 2012-10-24 上海交通大学 Grid network optimization method based on shortest feasible path and historical flow data
CN102880921A (en) * 2012-10-16 2013-01-16 山东电力集团公司电力科学研究院 Method for optimizing site selection of electric vehicle charging stations
CN110046791A (en) * 2019-02-28 2019-07-23 中国电力科学研究院有限公司 A kind of energy source station scene building method and system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112700175A (en) * 2021-01-21 2021-04-23 潜能恒信能源技术股份有限公司 Construction method and system of underground energy and overground demand coupled distributed energy station
CN112700175B (en) * 2021-01-21 2024-03-01 潜能恒信能源技术股份有限公司 Method and system for constructing distributed energy station coupling underground energy and overground demand

Also Published As

Publication number Publication date
CN110929980B (en) 2022-06-24

Similar Documents

Publication Publication Date Title
CN105096006A (en) Method for optimizing a routing of an intelligent ammeter distributing vehicle
CN109492950B (en) Prediction method capable of meeting space saturation load of large area range based on GIS technology
CN108921327A (en) Shelf method for carrying, apparatus and system applied to goods to people's system
CN110443397A (en) A kind of order allocator
CN110717226A (en) Regional comprehensive energy system layout planning method considering energy network topological characteristics
CN110929980B (en) Planning method and device for regional comprehensive energy system
CN101159048A (en) Oil products delivery cistern car scheduling system and method thereof
CN106096815A (en) A kind of power distribution network operational efficiency Impact analysis method
CN106874598A (en) A kind of electric power meter is directly with the storage dispensing Network Optimization Design method under pattern
CN108876050A (en) A kind of setting of the main processing procedure of iron and steel enterprise's contract and automatic switching method
Guo et al. Joint optimization of LORA and spares stocks considering corrective maintenance time
CN101437049B (en) Method for calculating service provider number and time length based on probability under service architecture
CN106228243B (en) The station track metro depot Hui Ku arrangement method
CN104869154A (en) Distributed resource scheduling method for balancing resource credibility and user satisfaction
CN101819652A (en) Multistage workflow modeling method based on serial number
CN105825308A (en) Reservoir connectivity analysis method based on multi-target analysis
Yang et al. Two-stage robust optimization strategy for spatially-temporally correlated data centers with data-driven uncertainty sets
CN114297249B (en) New energy project optimal selection ordering method considering power grid demand degree
Moradipari et al. Pricing differentiated services in an electric vehicle public charging station network
Asl et al. A mathematical model and solving method for multi-depot and multi-level vehicle routing problem with fuzzy time windows
Huang et al. An intelligent energy management mechanism for electric vehicles
Chen et al. Using scheduled ordering to improve the performance of distribution supply chains
CN110210631A (en) A kind of formulating method of high-speed railway comprehensive maintenance scheme
Zhou et al. A genetic algorithm approach on reverse logistics optimization for product return distribution network
Wang et al. Optimization approach to station location of car sharing system

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