CN108805449A - Cooperative game method towards integrated energy system cost sharing and distribution of income - Google Patents
Cooperative game method towards integrated energy system cost sharing and distribution of income Download PDFInfo
- Publication number
- CN108805449A CN108805449A CN201810596949.0A CN201810596949A CN108805449A CN 108805449 A CN108805449 A CN 108805449A CN 201810596949 A CN201810596949 A CN 201810596949A CN 108805449 A CN108805449 A CN 108805449A
- Authority
- CN
- China
- Prior art keywords
- main body
- income
- alliance
- participating
- matrix
- 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.)
- Pending
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims abstract description 64
- 238000005553 drilling Methods 0.000 claims abstract description 20
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 20
- 230000005611 electricity Effects 0.000 claims description 23
- 238000010248 power generation Methods 0.000 claims description 10
- 230000035945 sensitivity Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000001737 promoting effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 230000002269 spontaneous effect Effects 0.000 claims 1
- 238000004146 energy storage Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000571 coke Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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/06313—Resource planning in a project environment
-
- 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)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The cooperative game method towards integrated energy system cost sharing and distribution of income that the invention discloses a kind of can analyze each fairness for participating in main body cooperation result, and this method includes:It obtains in integrated energy system and participates in main information;According to main body quantity is participated in, first alliance's matrix is established;According to the participation type of subject and free generalization energy resource system and the correspondence for participating in main body, first alliance matrix is converted, second alliance's matrix is obtained;The controllable capacity that main body is participated according to integrated energy system, establishes contributrion margin column vector;Increased totle drilling cost and increased total revenue under main body modality for co-operation, which are participated in, according to integrated energy system calculates cooperation surplus;Each earnings pattern for participating in main body is established according to Shapley allocation rules and second alliance matrix, the contributrion margin column vector and cooperation surplus;According to each cost for participating in main body and the income, net profit is calculated;Key parameter sensitivity analysis finally is carried out to net profit.
Description
Technical field
The present invention relates to industrial park comprehensive energy technical fields, and in particular to towards integrated energy system cost sharing and
The cooperative game method of distribution of income.
Background technology
As the propulsion of China's energy technology revolution process and power system reform are deepened, in integrated energy system comprehensively
It participates in main body and tends to diversification, and subject of operation has the chance directly to user's sale of electricity.Industrial park integrated energy system one
As with natural gas cold, heat and power triple supply system (CCHP) be core energy conversion device, while adaptation to local conditions configure distributed photovoltaic
Other power supplys such as electricity generation system, distributed energy storage system, therefore also expedited the emergence of the subject of operation of diversification.In integrated energy system
The game behaviors participated between main body and operation mode have a significant impact to the overall efficiency of each side's income and system more.Excellent
While changing configuration garden resource, fairly and reasonably income is distributed, overhead cost is each coke for participating in main body and being paid close attention in game
Point.
It is existing to study the non-cooperation operation mode for being more focused on multi-energy system, to each participation main body cooperation operation mode
Game behavior and cost sharing and income distribution method research are less.Generally be integrated in the integrated energy system of industrial park with
A variety of secondary energy sources workload demands such as hot and cold, electric based on industrial user, have that demand is big, part throttle characteristics is complicated, reliability
It is required that the features such as high, and power grid, cold heat net need to be mutually coordinated in garden, to achieve the purpose that improve efficiency of energy utilization.Industry
Garden integrated energy system takes the operation mode of cooperative game, can increase each interests for participating in main body, be improved more for user
Good energy supply service, improves industrial park whole economic efficiency and social benefit.
Therefore a kind of method that can be analyzed each fairness for participating in main body cooperation result is urgently provided.In order to
While improving integrated energy system overall operation benefit, ensure the justice of each cost sharing and distribution of income for participating in main body
Rationally.
Invention content
The cooperative game method towards integrated energy system cost sharing and distribution of income that the present invention provides a kind of, can
The method that each fairness for participating in main body cooperation result is analyzed.In order to improve integrated energy system overall operation effect
While beneficial, ensure it is each participate in main body cost sharing and distribution of income it is fair and reasonable.
A kind of cooperative game method towards integrated energy system cost sharing and distribution of income includes in the present invention:
It obtains in integrated energy system and participates in main information, including participating in type of subject, quantity information and arbitrary
Integrated energy system and the correspondence for participating in main body;
According to the participation main body quantity, first alliance's matrix is established, first alliance matrix is (2n- 1) × n rank squares
Battle array;
According to the participation type of subject and free generalization energy resource system and the correspondence for participating in main body, to institute
It states first alliance's matrix to be converted, obtains second alliance's matrix, second alliance matrix is (2n- 1) × n rank 0-1 matrixes,
Every a line in matrix represents an alliance, and the serial number of each row represents a participation main body;" 0 " element representation in matrix is not
Cooperative game is participated in, " 1 " element, which represents, participates in cooperative game;Matrix constitutes a binary number per a line each element, and from the
Two rows start, and the binary number that each row element is constituted is equal to previous row element composition binary number and adds 1;
The controllable capacity that main body is participated according to integrated energy system, establishes contributrion margin column vector;
Increased totle drilling cost and increased total revenue under main body modality for co-operation, which are participated in, according to integrated energy system calculates cooperation
It is remaining;
According to Shapley allocation rules and second alliance matrix, the contributrion margin column vector and cooperation surplus
Establish each earnings pattern for participating in main body;
According to each cost for participating in main body and the income, net profit is calculated;
Key parameter sensitivity analysis is carried out to the net profit, and draws sensitivity analysis song in same reference axis
Line can be intersected in same point if respectively participating in the sensitivity curves of main body, then proving that the distribution method is fair.
Optionally,
The controllable capacity for participating in main body described in step according to integrated energy system, establishes contributrion margin column vector, wherein joining
The contributrion margin of alliance is expressed as with main body;
Wherein, υ (A) is the overall utility that all participation main bodys constitute alliance, if υ (A) represents income, the calculating of f (q)
As a result indicate that q-th of contributrion margin for participating in main body to alliance's income, υ (A) indicate as the same when cost;Sq, SkIt is expressed as
The q controllable resources capacity that main body is respectively participated in kth;
Optionally,
Increased totle drilling cost and increased total revenue under main body modality for co-operation are participated according to integrated energy system described in step
Further include before calculating cooperation surplus:
It obtains compared to each summation for participating in main body additional income before cooperation;
The increased totle drilling cost of main body needs is respectively participated in after obtaining cooperation.
Optionally,
It obtains and includes compared to each summation for participating in main body additional income before cooperation:Obtain increased peak clipping income, CCHP
Running optimizatin income and promoted photovoltaic from consumption rate and increased income, wherein peak clipping income are expressed as B1=Δ Qhl·Ph, B1For
Increased peak clipping income, Δ QhlFor increased peak clipping electricity, PhFor local peak value electricity price;
Optimization operation income is expressed as B2=Δ Qmove·(Ph-Pl), wherein B2Optimize operation income for CCHP;ΔQmoveFor
The CCHP generated energy of crest segment is moved to by paddy section;PlFor local valley electricity price;
Distributed photovoltaic power generation project total installed capacity is expressed as B3=CPV·HPV·Δβ·(Pre_load-Pre_PV), wherein CPV
For distributed photovoltaic power generation project total installed capacity;HPVFor distributed photovoltaic power generation annual utilization hours;Δ β is distributed photovoltaic power generation
It generates power for their own use increased ratio;Pre_loadFor local time-of-use tariffs mean value;Pre_PVFor photovoltaic online electricity price;
Increased total revenue is expressed as B under cooperation operation modeAll=B1+B2+B3。
Optionally,
The increased totle drilling cost of main body needs is respectively participated in after obtaining cooperation described in step, it is wherein increased under cooperation operation mode
Totle drilling cost is expressed as
Wherein, CAllFor increased totle drilling cost under cooperation operation mode;M indicates the electric wiring quantity for needing to be transformed in total;
LiFor the length of i-th distribution line;CiFor the unit length cost on i-th line road;N is the hot and cold tube netting index amount of transformation, WjFor
The construction length of the cold and hot pipe network of j-th strip;Chc_jFor the cold and hot pipe network unit length cost of j-th strip;CsysFor building for garden electrical secondary system
It is set as this.
Optionally,
According to Shapley allocation rules and second alliance matrix, contributrion margin column vector and described described in step
Cooperation surplus establishes each earnings pattern for participating in main body, and the earnings pattern is shown as:
In formula, φi(υ) indicates game person aiThe income of gained;E (S) indicates the total revenue of alliance S;S-{aiIndicate in S
Game person's set of ai is removed, π (S) indicates the probability that S occurs in coalition structure;| S | indicate game person's number in alliance S
Amount, n are the sum of game person in feasible zone;CSFor the second alliance battle array;CctrFor contributrion margin column vector;P is cooperation surplus.
Participating in main body cooperation interpretation of result device the present invention also provides a kind of integrated energy system includes:
Acquiring unit participates in main information for obtaining in integrated energy system;
Alliance's matrix establishes unit, for according to the participation main body quantity, establishing first alliance's matrix, described first
Alliance's matrix is (2n-1) × n rank matrixes;
Alliance's matrix conversion unit, for according to the participation type of subject and free generalization energy resource system and the ginseng
With the correspondence of main body, first alliance matrix is converted, obtains second alliance's matrix, second alliance matrix
It is (2n- 1) × n rank 0-1 matrixes;
Contributrion margin column vector establishes unit, and the controllable capacity for participating in main body according to integrated energy system establishes side
Contribute column vector in border;
Cooperation surplus computing unit, for according to integrated energy system participate in main body modality for co-operation under increased totle drilling cost and
Increased total revenue calculates cooperation surplus;
Earnings pattern establishes unit, for according to Shapley allocation rules and second alliance matrix, marginal tribute
It offers column vector and the cooperation surplus establishes each earnings pattern for participating in main body;
Net profit computing unit, for according to each cost for participating in main body and the income, calculating net profit;
Key parameter sensitivity analysis unit, for carrying out key parameter sensitivity analysis to the net profit, and same
Sensitivity analysis curve is drawn in one reference axis, same point can be intersected in if respectively participating in the sensitivity curves of main body, then demonstrate,proving
The bright distribution method is fair.
Optionally,
Described device further includes:
First acquisition unit, for obtaining compared to each summation for participating in main body additional income before cooperation;
Second acquisition unit, for respectively participating in the increased totle drilling cost of main body needs after obtaining cooperation.
Optionally, first acquisition unit includes:
Subelement is obtained, for obtaining increased peak clipping income, CCHP running optimizatins income and promoting photovoltaic from consumption rate
And increased income;
Computation subunit, for being dissolved certainly to the increased peak clipping income, CCHP running optimizatins income and promotion photovoltaic
Rate and increased income carry out sum operation.
The invention has the advantages that:
It obtains in integrated energy system first and participates in main information;According to main body quantity is participated in, first alliance's matrix is established;
According to the participation type of subject and free generalization energy resource system and the correspondence for participating in main body, to described first
Alliance's matrix is converted, and second alliance's matrix is obtained;The controllable capacity of main body is participated according to integrated energy system, establishes marginal tribute
Offer column vector;Increased totle drilling cost and increased total revenue under main body modality for co-operation, which are participated in, according to integrated energy system calculates cooperation
It is remaining;It is built according to Shapley allocation rules and second alliance matrix, the contributrion margin column vector and cooperation surplus
Vertical each earnings pattern for participating in main body;According to each cost for participating in main body and the income, net profit is calculated;Finally to net profit
Key parameter sensitivity analysis is carried out, and draws sensitivity analysis curve in same reference axis, if respectively participating in the sensitivity of main body
Linearity curve can be intersected in same point, then proving that the distribution method is fair.
Description of the drawings
Fig. 1 is a kind of towards the implementation of the cooperative game method of integrated energy system cost sharing and distribution of income in the present invention
Example flow chart;
Fig. 2 is a kind of integrated energy system participation main body cooperation interpretation of result device embodiment structural representation in the present invention
Figure;
Fig. 3 is a kind of towards the implementation of the cooperative game method of integrated energy system cost sharing and distribution of income in the present invention
The cooperative game relation schematic diagram of main body is respectively participated in example;
Fig. 4 is a kind of towards the implementation of the cooperative game method of integrated energy system cost sharing and distribution of income in the present invention
A key parameter sensitivity analysis result schematic diagram in example;
Fig. 5 is a kind of towards the implementation of the cooperative game method of integrated energy system cost sharing and distribution of income in the present invention
Another key parameter sensitivity analysis result schematic diagram in example.
Specific implementation mode
The cooperative game method towards integrated energy system cost sharing and distribution of income that the present invention provides a kind of, can
The method that each fairness for participating in main body cooperation result is analyzed.In order to improve integrated energy system overall operation effect
While beneficial, ensure it is each participate in main body cost sharing and distribution of income it is fair and reasonable.
A kind of cooperative game method towards integrated energy system cost sharing and distribution of income includes in the present invention:
101, obtain integrated energy system in participate in main information, including participate in type of subject, quantity information and
Free generalization energy resource system and the correspondence for participating in main body;
In the present embodiment, device, which is configured as obtaining in integrated energy system, participates in main information, including participation
Type of subject, quantity information and free generalization energy resource system and the correspondence for participating in main body;Such as:In industrial park
3 participation main bodys, respectively grid company, genco and sale of electricity company are shared, number consecutively is 1,2 and 3 respectively, can also be remembered
It is set A={ 1,2,3 }.
102, according to the participation main body quantity, first alliance's matrix is established, first alliance matrix is (2n-1)×n
Rank matrix;
In the present embodiment, system is configured as, according to the participation main body quantity, establishing first alliance's matrix;Such as 3
A participation main body then has 23There is kind of a combination.
103, type of subject and free generalization energy resource system and the correspondence for participating in main body are participated according to described,
First alliance matrix is converted, second alliance's matrix is obtained, second alliance matrix is (2n- 1) × n ranks 0-1
Matrix.
Every a line in matrix represents an alliance, and the serial number of each row represents a participation main body;" 0 " member in matrix
Element expression is not involved in cooperative game, and " 1 " element, which represents, participates in cooperative game;Matrix constitutes a binary system per a line each element
Number, and since the second row, the binary number that each row element is constituted adds 1 equal to previous row element composition binary number.
Such as the corresponding second alliance matrix of 3 participation main bodys is;
104, the controllable capacity that main body is participated according to integrated energy system, establishes contributrion margin column vector;
In the present embodiment, system is configured as respectively participating in the controllable resources of main body according to industrial park integrated energy system
Capacity is expressed as the contributrion margin of alliance wherein participating in main body
Wherein, υ (A) is the overall utility that all participation main bodys constitute alliance, if υ (A) represents income, the calculating of f (q)
As a result indicate that q-th of contributrion margin for participating in main body to alliance's income, υ (A) indicate as the same when cost;Sq, SkIt is expressed as
The q controllable resources capacity that main body is respectively participated in kth;
Such as shown in table 2.The contributrion margin collection of formation is combined into { 0.5357,0.2857,0.1786 }.
Participate in main body | Controllable resources type | Controllable capacity | Contributrion margin |
Grid company | Battery energy storage | 6MW | 0.5357 |
Genco | Battery energy storage+CCHP | 2MW+3MW | 0.2857 |
Sale of electricity company | Battery energy storage | 2MW | 0.1786 |
1 industrial park integrated energy system of table respectively participates in main body controllable resources situation table
105, increased totle drilling cost and increased total revenue under main body modality for co-operation is participated according to integrated energy system to calculate
Cooperation surplus;
In the present embodiment, system is configured as participating in increased assembly under main body modality for co-operation according to integrated energy system
This and increased total revenue calculate cooperation surplus and are modeled to the cooperative game behavior for respectively participating in main body in integrated energy system,
Establish the computation model of cooperation surplus.Obtained by respectively being participated in the integrated energy system of industrial park in embodiment after main body cooperation
Extra returns include increased peak clipping income, CCHP running optimizatins income and promote photovoltaic from consumption rate and increased income, point
Not can by formula (1)~(3) calculate and obtain, involved in formula to calculating parameter refer to table 3.
B1=Δ Qhl·Ph (1)
In formula:B1For increased peak clipping income;ΔQhlFor increased peak clipping electricity;PhFor local peak value electricity price.
B2=Δ Qmove·(Ph-Pl) (2)
In formula:B2Optimize operation income for CCHP;ΔQmoveTo be moved to the CCHP generated energy of crest segment by paddy section;PlTo work as
The valley electricity price on ground.
B3=CPV·HPV·Δβ·(Pre_load-Pre_PV) (3)
In formula:CPVFor distributed photovoltaic power generation project total installed capacity;HPVFor distributed photovoltaic power generation annual utilization hours;Δβ
It generates power for their own use increased ratio for distributed photovoltaic power generation;Pre_loadFor local time-of-use tariffs mean value;Pre_PVFor photovoltaic online
Electricity price.
Increased total revenue B under cooperation operation modeAllIt can be calculated by formula (4):
BAll=B1+B2+B3 (4)
Increased cost of investment includes garden two after industrial park integrated energy system respectively participates in main body cooperation in embodiment
Subsystem construction cost, electric wiring improvement cost and cold and hot pipe network improvement cost can be calculated by formula (5), involved in calculating
To parameter refer to table 2 and table 3.
In formula, CAllFor increased totle drilling cost under cooperation operation mode;M indicates the electric wiring quantity for needing to be transformed in total;
LiFor the length of i-th distribution line;CiFor the unit length cost on i-th line road;N is the hot and cold tube netting index amount of transformation, WjFor
The construction length of the cold and hot pipe network of j-th strip;Chc_jFor the cold and hot pipe network unit length cost of j-th strip;CsysFor building for garden electrical secondary system
It is set as this.
Parameter | Numerical value | Parameter | Numerical value |
Peak value electricity price | 100.39 point/kilowatt hour | Flat section electricity price | 60.84 point/kilowatt hour |
Paddy section electricity price | 30.42 point/kilowatt hour | Photovoltaic online electricity price | 0.75 yuan/kilowatt hour |
Unit price is averagely transformed in electric wiring | 2500000 yuan/km | Electrical secondary system construction cost | 20000000 yuan |
Unit price is averagely transformed in cold and hot pipe network | 500000 yuan/km |
2 industrial park relevant rudimentary data of table
Parameter | Numerical value | Parameter | Numerical value |
Length is transformed in electric wiring | 10 kms | Increase peak clipping electricity | 294.7MWh |
Cold and hot pipe network improving length | 10 kms | CCHP translates electricity | 6000000 kilowatt hours |
3 industrial park integrated energy system of table participates in main body cooperation surplus more and calculates relevant parameter
106, according to Shapley allocation rules and second alliance matrix, the contributrion margin column vector and cooperation
Residue establishes each earnings pattern for participating in main body;The earnings pattern is shown as:
In formula, φi(υ) indicates game person aiThe income of gained;E (S) indicates the total revenue of alliance S;S-{aiIndicate in S
Remove aiGame person set, π (S) indicate in coalition structure S occur probability;| S | indicate game person's number in alliance S
Amount, n are the sum of game person in feasible zone;CSFor the second alliance battle array;CctrFor contributrion margin column vector;P is cooperation surplus.
107, according to each cost for participating in main body and the income, net profit is calculated;
System is configured as, according to each cost for participating in main body and the income, calculating net profit, tool in the present embodiment
Body operation method is not specifically limited, and checkout result refers to table 4
Grid company | Genco | Sale of electricity company | It is total | |
Overhead cost (ten thousand yuan) | 133.93 | 71.43 | 44.64 | 250.00 |
Distribute income (ten thousand yuan) | 242.92 | 129.56 | 80.97 | 453.45 |
Net profit (ten thousand yuan) | 108.99 | 58.13 | 36.33 | 203.45 |
Table 4 respectively participates in main body cost sharing and distribution of income result
108, key parameter sensitivity analysis is carried out to the net profit, and sensitivity analysis is drawn in same reference axis
Curve can be intersected in same point if respectively participating in the sensitivity curves of main body, then proving that the distribution method is fair.
Refer to Fig. 4 and Fig. 5.
Participating in main body cooperation interpretation of result device the present invention also provides a kind of integrated energy system includes:
Acquiring unit 201 participates in main information for obtaining in integrated energy system;
Alliance's matrix establishes unit 202, for according to the participation main body quantity, establishing first alliance's matrix, described the
One alliance's matrix is (2n-1) × n rank matrixes;
Alliance's matrix conversion unit 203, for according to the participation type of subject and free generalization energy resource system and institute
The correspondence for participating in main body is stated, first alliance matrix is converted, second alliance's matrix, second alliance are obtained
Matrix is (2n- 1) × n rank 0-1 matrixes;
Contributrion margin column vector establishes unit 204, and the controllable capacity for participating in main body according to integrated energy system is established
Contributrion margin column vector;
Cooperation surplus computing unit 205, for participating in increased assembly under main body modality for co-operation according to integrated energy system
This and increased total revenue calculate cooperation surplus;
Earnings pattern establishes unit 206, for according to Shapley allocation rules and second alliance matrix, limit
Contribution column vector and the cooperation surplus establish each earnings pattern for participating in main body;
Net profit computing unit 207, for according to each cost for participating in main body and the income, calculating net profit;
Key parameter sensitivity analysis unit 208 is used for net profit progress key parameter sensitivity analysis, and
Sensitivity analysis curve is drawn in same reference axis, and same point can be intersected in if respectively participating in the sensitivity curves of main body, then
Prove that the distribution method is fair.
The device further comprises:
First acquisition unit 209, for obtaining compared to each summation for participating in main body additional income before cooperation;
Second acquisition unit 210, for respectively participating in the increased totle drilling cost of main body needs after obtaining cooperation.
Wherein, first acquisition unit 209 includes:
Subelement 2091 is obtained, is disappeared certainly for obtaining increased peak clipping income, CCHP running optimizatins income and promoting photovoltaic
Receive rate and increased income;
Computation subunit 2092 is used for the increased peak clipping income, CCHP running optimizatins income and promotes photovoltaic certainly
Consumption rate and increased income carry out sum operation.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art can easily think of the change or the replacement in the technical scope of present disclosure, should all contain
Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (9)
1. a kind of cooperative game method towards integrated energy system cost sharing and distribution of income, which is characterized in that including:
It obtains in integrated energy system and participates in main information, including participation type of subject, quantity information and free generalization
Energy resource system and the correspondence for participating in main body;
According to the participation main body quantity, first alliance's matrix is established, first alliance matrix is (2n- 1) × n rank matrixes;
According to the participation type of subject and free generalization energy resource system and the correspondence for participating in main body, to described the
One alliance's matrix is converted, and second alliance's matrix is obtained, and second alliance matrix is (2n- 1) × n rank 0-1 matrixes, matrix
In every a line represent an alliance, the serial number of each row represents a participation main body;" 0 " element representation in matrix is not involved in
Cooperative game, " 1 " element, which represents, participates in cooperative game;Matrix constitutes a binary number per a line each element, and from the second row
Start, the binary number that each row element is constituted is equal to previous row element composition binary number and adds 1;
The controllable capacity that main body is participated according to integrated energy system, establishes contributrion margin column vector;
Increased totle drilling cost and increased total revenue under main body modality for co-operation, which are participated in, according to integrated energy system calculates cooperation surplus;
It is established according to Shapley allocation rules and second alliance matrix, the contributrion margin column vector and cooperation surplus
Each earnings pattern for participating in main body;
According to each cost for participating in main body and the income, net profit is calculated;
Key parameter sensitivity analysis is carried out to the net profit, and draws sensitivity analysis curve in same reference axis, if
Each sensitivity curves for participating in main body can be intersected in same point, then proving that the distribution method is fair.
2. the cooperative game method according to claim 1 towards integrated energy system cost sharing and distribution of income,
It is characterized in that,
The controllable capacity for participating in main body described in step according to integrated energy system, establishes contributrion margin column vector, wherein participating in master
Body is expressed as the contributrion margin of alliance;
Wherein, υ (A) is the overall utility that all participation main bodys constitute alliance, if υ (A) represents income, the result of calculation of f (q)
Indicate that q-th of contributrion margin for participating in main body to alliance's income, υ (A) indicate as the same when cost;Sq, Sk are expressed as q-th
The controllable resources capacity of main body is respectively participated in kth.
3. the cooperative game method according to claim 1 towards integrated energy system cost sharing and distribution of income,
It is characterized in that,
Increased totle drilling cost and increased total revenue under main body modality for co-operation are participated in described in step according to integrated energy system to calculate
Further include before cooperation surplus:
It obtains compared to each summation for participating in main body additional income before cooperation;
The increased totle drilling cost of main body needs is respectively participated in after obtaining cooperation.
4. the cooperative game method according to claim 3 towards integrated energy system cost sharing and distribution of income,
It is characterized in that,
It obtains and includes compared to each summation for participating in main body additional income before cooperation:Obtain increased peak clipping income, CCHP operations
Optimize income and promoted photovoltaic from consumption rate and increased income, wherein peak clipping income are expressed as B1=Δ Qhl·Ph, B1To increase
Peak clipping income, Δ QhlFor increased peak clipping electricity, PhFor local peak value electricity price;
Optimization operation income is expressed as B2=Δ Qmove·(Ph-Pl), wherein B2Optimize operation income for CCHP;ΔQmoveFor by paddy
Section is moved to the CCHP generated energy of crest segment;PlFor local valley electricity price;
Distributed photovoltaic power generation project total installed capacity is expressed as B3=CPV·HPV·Δβ·(Pre_load-Pre_PV), wherein CPVTo divide
Cloth photovoltaic generation project total installed capacity;HPVFor distributed photovoltaic power generation annual utilization hours;Δ β is that distributed photovoltaic power generation is spontaneous
Personal increased ratio;Pre_loadFor local time-of-use tariffs mean value;Pre_PVFor photovoltaic online electricity price;
Increased total revenue is expressed as B under cooperation operation modeAll=B1+B2+B3。
5. the cooperative game method according to claim 3 towards integrated energy system cost sharing and distribution of income,
It is characterized in that,
The increased totle drilling cost of main body needs, wherein increased assembly under cooperation operation mode are respectively participated in after obtaining cooperation described in step
Originally it is expressed as
Wherein, CAllFor increased totle drilling cost under cooperation operation mode;M indicates the electric wiring quantity for needing to be transformed in total;LiFor
The length of i-th distribution line;CiFor the unit length cost on i-th line road;N is the hot and cold tube netting index amount of transformation, WjFor jth
The construction length of the cold and hot pipe network of item;Chc_jFor the cold and hot pipe network unit length cost of j-th strip;CsysFor the construction of garden electrical secondary system
Cost.
6. the cooperative game method according to claim 1 towards integrated energy system cost sharing and distribution of income,
It is characterized in that,
According to Shapley allocation rules and second alliance matrix, the contributrion margin column vector and cooperation described in step
Residue establishes each earnings pattern for participating in main body, and the earnings pattern is shown as:
In formula, φi(e) game person a is indicatediThe income of gained;E (S) indicates the total revenue of alliance S;S-{aiIndicate to remove a in Si
Game person set, π (S) indicate in coalition structure S occur probability;| S | indicate game person's quantity in alliance S, n is can
The sum of game person in row domain;CSFor the second alliance battle array;CctrFor contributrion margin column vector;P is cooperation surplus.
7. a kind of integrated energy system participates in main body cooperation interpretation of result device, which is characterized in that including:
Acquiring unit participates in main information for obtaining in integrated energy system;
Alliance's matrix establishes unit, for according to the participation main body quantity, establishing first alliance's matrix, first alliance square
Battle array is (2n-1) × n rank matrixes;
Alliance's matrix conversion unit, for being led with the participation according to the participation type of subject and free generalization energy resource system
The correspondence of body converts first alliance matrix, obtains second alliance's matrix, and second alliance matrix is
(2n-1) × n rank 0-1 matrixes;
Contributrion margin column vector establishes unit, the controllable capacity for participating in main body according to integrated energy system, establishes marginal tribute
Offer column vector;
Cooperation surplus computing unit, for participating in increased totle drilling cost and increase under main body modality for co-operation according to integrated energy system
Total revenue calculate cooperation surplus;
Earnings pattern establishes unit, for being arranged according to Shapley allocation rules and second alliance matrix, contributrion margin
Cooperation surplus described in vector sum establishes each earnings pattern for participating in main body;
Net profit computing unit, for according to each cost for participating in main body and the income, calculating net profit;
Key parameter sensitivity analysis unit, for carrying out key parameter sensitivity analysis to the net profit, and in same seat
Sensitivity analysis curve is drawn in parameter, same point can be intersected in if respectively participating in the sensitivity curves of main body, then proving should
Distribution method is fair.
8. integrated energy system according to claim 7 participates in main body cooperation interpretation of result device, which is characterized in that
Described device further includes:
First acquisition unit, for obtaining compared to each summation for participating in main body additional income before cooperation;
Second acquisition unit, for respectively participating in the increased totle drilling cost of main body needs after obtaining cooperation.
9. integrated energy system according to claim 8 participates in main body cooperation interpretation of result device, which is characterized in that first
Acquiring unit includes:
Subelement is obtained, is increased from consumption rate for obtaining increased peak clipping income, CCHP running optimizatins income and promoting photovoltaic
The income added;
Computation subunit, for the increased peak clipping income, CCHP running optimizatins income and promote photovoltaic from consumption rate and
Increased income carries out sum operation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810596949.0A CN108805449A (en) | 2018-06-11 | 2018-06-11 | Cooperative game method towards integrated energy system cost sharing and distribution of income |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810596949.0A CN108805449A (en) | 2018-06-11 | 2018-06-11 | Cooperative game method towards integrated energy system cost sharing and distribution of income |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108805449A true CN108805449A (en) | 2018-11-13 |
Family
ID=64088445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810596949.0A Pending CN108805449A (en) | 2018-06-11 | 2018-06-11 | Cooperative game method towards integrated energy system cost sharing and distribution of income |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108805449A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109784714A (en) * | 2019-01-11 | 2019-05-21 | 上海电机学院 | A kind of station P2G operation mode under multiagent based on Shapley value method is joint |
CN110110008A (en) * | 2019-05-08 | 2019-08-09 | 北京理工大学 | A kind of shared motivational techniques of the block chain medical data based on Charolais cattle |
CN112257951A (en) * | 2020-11-02 | 2021-01-22 | 国网安徽省电力有限公司合肥供电公司 | Optimized operation method of comprehensive energy system and power distribution company based on cooperative game |
CN112669052A (en) * | 2020-12-04 | 2021-04-16 | 国网辽宁省电力有限公司经济技术研究院 | Flexible income distribution method for park comprehensive energy system participating in peak shaving |
WO2022011968A1 (en) * | 2020-07-17 | 2022-01-20 | 广西电网有限责任公司电力科学研究院 | Multi-agent investment proportion optimization method and system based on cooperative game |
-
2018
- 2018-06-11 CN CN201810596949.0A patent/CN108805449A/en active Pending
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109784714A (en) * | 2019-01-11 | 2019-05-21 | 上海电机学院 | A kind of station P2G operation mode under multiagent based on Shapley value method is joint |
CN110110008A (en) * | 2019-05-08 | 2019-08-09 | 北京理工大学 | A kind of shared motivational techniques of the block chain medical data based on Charolais cattle |
CN110110008B (en) * | 2019-05-08 | 2020-12-25 | 北京理工大学 | Block chain medical data sharing excitation method based on charapril value |
WO2022011968A1 (en) * | 2020-07-17 | 2022-01-20 | 广西电网有限责任公司电力科学研究院 | Multi-agent investment proportion optimization method and system based on cooperative game |
CN112257951A (en) * | 2020-11-02 | 2021-01-22 | 国网安徽省电力有限公司合肥供电公司 | Optimized operation method of comprehensive energy system and power distribution company based on cooperative game |
CN112257951B (en) * | 2020-11-02 | 2023-12-19 | 国网安徽省电力有限公司合肥供电公司 | Comprehensive energy system and power distribution company optimized operation method based on cooperative game |
CN112669052A (en) * | 2020-12-04 | 2021-04-16 | 国网辽宁省电力有限公司经济技术研究院 | Flexible income distribution method for park comprehensive energy system participating in peak shaving |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108805449A (en) | Cooperative game method towards integrated energy system cost sharing and distribution of income | |
Zhang et al. | Unit commitment model in smart grid environment considering carbon emissions trading | |
Jin et al. | Game theoretical analysis on capacity configuration for microgrid based on multi-agent system | |
Ding et al. | Economic and environmental assessment of multi-energy microgrids under a hybrid optimization technique | |
CN103455852A (en) | Power transmission and distribution cost allocation method based on DEA cooperative game | |
Jia et al. | Joint operation and transaction mode of rural multi microgrid and distribution network | |
CN105375472B (en) | A kind of intelligent distribution network low-carbon performance evaluation method | |
CN104915790A (en) | Time of use price optimization method for promoting wind power generation | |
CN116911533B (en) | Multi-microgrid energy sharing method for regional comprehensive energy system | |
CN108155673B (en) | Power system optimal scheduling method considering uncertainty of combined power generation at power generation side | |
Chupradit et al. | Modeling and Optimizing the Charge of Electric Vehicles with Genetic Algorithm in the Presence of Renewable Energy Sources | |
CN115864376A (en) | Multi-time scale power supply conservation analysis method and system | |
CN105226707A (en) | A kind of methodology based on Shapley value wind-electricity integration system fixed cost of power transmission | |
Yang et al. | Optimal scheduling of electro-thermal system considering refined demand response and source-load-storage cooperative hydrogen production | |
Yang et al. | A two-stage operation optimization model for isolated integrated energy systems with concentrating solar power plant considering multi-energy and multi-type demand response | |
CN113988471A (en) | Multi-objective optimization method for micro-grid operation | |
Liu et al. | Analysis of flexible energy use behavior of rural residents based on two-stage questionnaire: A case study in Xi’an, China | |
CN108062721A (en) | Industrial park integrated energy system income distribution method, device and system | |
Tang et al. | Combined electricity-gas-heat energy internet scheduling with power-to-gas and renewable power uncertainty | |
CN115663832A (en) | Source network and storage interactive game operation strategy for multi-energy complementary power distribution system | |
CN115906456A (en) | Hydrogen-containing energy IES scheduling optimization model considering response uncertainty of demand side | |
CN115358559A (en) | Comprehensive energy system source storage capacity optimal configuration method considering comprehensive demand response | |
CN114336607A (en) | Power grid economic operation domain generation method and system considering source-load bilateral uncertainty | |
Eghlimi et al. | Decision model for cross-border electricity trade considering renewable energy sources | |
Na et al. | Coordination and optimization model of wind power grid integrated system considering demand response |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181113 |
|
RJ01 | Rejection of invention patent application after publication |