CN111260158A - Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation - Google Patents

Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation Download PDF

Info

Publication number
CN111260158A
CN111260158A CN202010117662.2A CN202010117662A CN111260158A CN 111260158 A CN111260158 A CN 111260158A CN 202010117662 A CN202010117662 A CN 202010117662A CN 111260158 A CN111260158 A CN 111260158A
Authority
CN
China
Prior art keywords
time
gas
interest
power
main body
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
Application number
CN202010117662.2A
Other languages
Chinese (zh)
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.)
Economic and Technological Research Institute of State Grid Sichuan Electric Power Co Ltd
Original Assignee
Economic and Technological Research Institute of State Grid Sichuan 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 Economic and Technological Research Institute of State Grid Sichuan Electric Power Co Ltd filed Critical Economic and Technological Research Institute of State Grid Sichuan Electric Power Co Ltd
Priority to CN202010117662.2A priority Critical patent/CN111260158A/en
Publication of CN111260158A publication Critical patent/CN111260158A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q30/00Commerce
    • G06Q30/02Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201Market modelling; Market analysis; Collecting market data
    • G06Q30/0206Price or cost determination based on market factors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S50/00Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/14Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards

Landscapes

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

Abstract

The invention discloses a market multi-interest main body transaction behavior modeling method based on internal interactive cooperation, which is characterized in that an energy hub basic framework of a market multi-interest main body is constructed based on the basic principle of hub input and output; respectively carrying out independent modeling analysis on multi-benefit subjects interactively cooperating in the energy concentrator, namely carrying out modeling analysis on an electric benefit subject, a gas benefit subject and a cogeneration benefit subject; establishing a coupling matrix to represent the incidence relation among all benefit subjects; and analyzing and calculating the obtained optimal conclusion through the coupling matrix, and guiding the reasonable distribution of the multi-benefit subjects according to the optimal conclusion. According to the invention, through quantifying the energy flow of the power connecting line and the natural gas pipeline between the system and the main network, the electric heating load variation trend is respectively analyzed, and an optimal conclusion is obtained under the constraint of an objective function.

Description

Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation
Technical Field
The invention belongs to the technical field of multi-energy main body cooperation management, and particularly relates to a market multi-benefit main body transaction behavior modeling method based on internal interactive cooperation.
Background
Energy is the basis for human survival and development and is the life line of national economy. Renewable energy is a great measure for improving the energy structure of China, promoting environmental protection and maintaining the sustainable development of society and economy. In the process of considering energy development, electric energy is no longer in a single energy form, and multiple energy forms such as gas heat greatly influence the proportion of electric energy in an energy structure. However, the development of gas-heat energy has far influence on the environment beyond electric energy, so that a multi-energy system formed by reasonably planning and designing electricity, heat and gas and optimizing operation provides an important solution for realizing the change of the energy field. When each energy source participates in the market, the market behavior of each energy source is restricted by other energy source forms. Each energy form is an independent benefit main body, and the operation and optimization of the power distribution network are inevitably influenced by other energy forms under the framework of mainly electric energy, so that the influence on the electricity benefit main body under the condition that a distribution network area taking an intelligent park and the like as a typical demonstration area is considered, and the coupling of various energy forms is considered. The existing literature is used for researching a combined dispatching model of energy systems such as electricity-gas, electricity-heat and the like. The multi-energy system has complex space-time characteristics, and has different research methods for different time scales and space scales.
On the spatial scale, the research of the multi-energy system can follow 'multi-energy element-multi-energy element integration-multi-energy network', and the space is divided into regional multi-energy systems which mainly adopt electric heating coupling and electric coupling; on the time scale, the research of the multi-energy system can be divided into dynamic problems on short time scales such as seconds and minutes, and dynamic problems on time scales such as hours, days and months. The main body of the multi-energy system includes not only single elements such as a thermal power plant, an electric boiler, etc., but also a combination of these elements, for example, a commercial building may include various elements such as a distributed natural gas power generation device, a distributed photovoltaic device, an electric boiler, a heat pump, etc. Therefore, it is necessary to study the integration characteristics of the device based on the characteristics of the device, and in academia, the integration of the multi-energy device is generally called energy hub. The main body of the multi-energy system comprises single elements such as a thermal power plant, an electric boiler and the like, and more, the single elements are higher-level packaging and combination of the elements, the elements are mutually coupled, the connection relationship is complex, and the analysis is often difficult. In fact, the most critical to a system is its input-output characteristics, i.e., the "external port" characteristics. For a multi-energy system, the problem to be solved is how to equate the coupling conversion relationship of internal elements and establish an external end port model for mapping the input and output of electricity, gas and heat. Under the concept of the energy hub, energy hubs are produced. The method can improve the flexibility of the energy system and can be widely applied to the problems of system planning, optimal power flow, demand side management, optimal operation and the like, the concept of optimal operation is mostly appeared in the field of power system analysis, and the method is popularized to the optimal operation of a multi-energy hybrid system. On the basis, the multi-energy-form participation market behavior is defined as participation of a plurality of interest-related bodies in the market, the power distribution network region is mainly researched, and how the multi-interest body internally couples and interactively cooperates to guide the market behavior is influenced, and the economy of the system is influenced.
Therefore, it is necessary to provide a modeling method for trading behavior of market multi-interest main body based on internal interactive collaboration.
Disclosure of Invention
The invention aims to provide a market multi-benefit-subject trading behavior modeling method based on internal interactive cooperation, aims to construct three different benefit subjects based on an energy hub and taking coordination and complementation of electricity, heat and gas as an example in a power distribution network area, and provides the market multi-benefit-subject trading behavior modeling method based on the internal interactive cooperation. Through quantifying the energy flow of the power connecting line and the natural gas pipeline between the system and the main network, including the output characteristics of each device in the hub in different time periods, a corresponding coupling matrix model is established, and the multi-dimensional characteristics of three stakeholders are modeled and analyzed. Meanwhile, the market behaviors of multi-benefit subjects are guided by considering the demand response characteristics of loads in multiple periods, analyzing the change trend of the connection electric load and the heat load of the energy hub and the electric purchasing condition and calculating the total electric purchasing cost. And finally, taking 3 typical parks in the multi-interest main body hub system as an example, the superiority of the distribution network multi-interest main body market trading mechanism based on the internal interactive cooperation of the energy hub is demonstrated.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the market multi-interest main body trading behavior modeling method based on internal interactive collaboration comprises the following steps:
step 1: constructing an energy hub basic framework of a market multi-interest main body based on the input and output basic principle of the hub;
step 2: respectively carrying out independent modeling analysis on the multi-benefit main bodies interactively cooperated in the energy hub in the step 1, namely carrying out modeling analysis on an electric benefit main body, a gas benefit main body and a cogeneration benefit main body;
and step 3: establishing a coupling matrix to represent the incidence relation among all the benefit bodies in the step 2;
and 4, step 4: estimating the photovoltaic and fan output according to historical data, measuring the tie line power of the system and calculating the net exchange power; measuring gas flow data of a natural gas pipeline, and setting a target function to calculate the gas and electricity acquisition cost; analyzing the change trend of the thermoelectric load according to the response change; and (4) analyzing and calculating the obtained optimal conclusion through the coupling matrix in the step (3), and guiding the reasonable distribution of the multi-benefit subjects according to the optimal conclusion.
Further, in step 1, the basic framework of the energy hub is as follows:
Figure BDA0002391992110000021
Figure BDA0002391992110000022
in the formula (1), the vector
Figure BDA0002391992110000023
And
Figure BDA0002391992110000024
representing inputs and outputs of different energies of a single hub, matrix lambdamnUsed as a coupling matrix describing the relationship between the input and the output, where the value is related to the conversion efficiency of the power supply device.
Further, the modeling analysis of the electric interest body in the step 2 specifically comprises:
the electric interest main body comprises photovoltaic, wind power generation, comprehensive dispatching of energy storage equipment and load demand response;
for photovoltaic and wind power generation, a random planning method is selected, and the uncertainty of a distributed power supply is described by using a probability density distribution function, wherein the photovoltaic output obeys beta distribution, and the fan output obeys two-parameter Weibu distribution;
for the comprehensive scheduling of the energy storage device, the energy storage device can be flexibly called in a scheduling period, and a dynamic model of the energy storage device can be expressed as follows:
Figure BDA0002391992110000031
Figure BDA0002391992110000032
in the formula (3), x is a binary variable and represents the charge-discharge state of the energy storage device at the ith day t, 1 is in a working state, 0 is in a non-working state,
Figure BDA0002391992110000033
in order to charge and discharge power in real time,
Figure BDA0002391992110000034
is a charge-discharge power limit; in the formula (4), Ei,0,Ei,24Representing the stored energy at the beginning and end of the energy storage device; es,t,jReal-time charge;
Figure BDA0002391992110000035
limiting the real-time charge capacity of the equipment;
Figure BDA0002391992110000036
the maximum charge capacity of the equipment; ei,tRepresenting the electric quantity stored by the energy storage equipment at adjacent moments, α and β are charge and discharge efficiencies, and are both set to be 0.9;
for the load demand response, the PBDR is selected as a research object in the load demand response project, the multi-interest subject is built, and the user is guided to carry out load distribution by responding the time-of-use electricity price signal in real time; the dual electricity price mechanism of the time-of-use electricity price and the real-time electricity price is selected, a user responds to the time-of-use electricity price and can participate in the PBDR project according to the real-time electricity price, and load demand calculation after PBDR is realized can be represented as follows:
Figure BDA0002391992110000037
in the formula:
Figure BDA0002391992110000038
in order to implement the load after PBDR,
Figure BDA0002391992110000039
initial load before PBDR implementation; t is a total scheduling period, and s and T are sub-periods; e.g. of the typestThe value is the price elastic coefficient, and s-t is the self elastic coefficient, otherwise, the value is the cross elastic coefficient;
Figure BDA00023919921100000310
to implement the real-time electricity rates after PBDR,
Figure BDA00023919921100000311
initial time-of-use electricity prices before PBDR is implemented; the real-time electricity price can be calculated according to the following formula:
Figure BDA0002391992110000041
in the formula:
Figure BDA0002391992110000042
is the total load demand for a day;
Figure BDA0002391992110000043
is the average load of a day.
Further, the modeling analysis on the gas benefit subject in the step 2 specifically comprises:
the interest subject of the gas boiler takes the gas boiler as a considered object, and the model of the interest subject of the gas boiler is represented by the coupling relation of the heating power and the gas consumption in unit time:
Figure BDA0002391992110000044
in the formula, muGB,thRepresents the electrical conversion efficiency;
Figure BDA0002391992110000045
the heat supply power of the gas-fired boiler is increased,
Figure BDA0002391992110000046
is the amount of natural gas consumed per unit time.
Further, the modeling analysis of the cogeneration interest body in the step 2 specifically comprises:
the combined heat and power benefit agent takes natural gas as a consideration object, and the power supply and heat supply quantity of the combined heat and power benefit agent is expressed as follows:
Figure BDA0002391992110000047
in the formula:
Figure BDA0002391992110000048
the electric quantity and the heating quantity of the cogeneration unit generated at the first day t,
Figure BDA0002391992110000049
is the consumption of natural gas per unit time, muCHP,ele,μCHP,thThe gas-electricity conversion efficiency and the gas-heat conversion efficiency are respectively.
Further, the coupling matrix in step 3 is specifically:
by modeling the coupling relationship among the internal devices of the multi-interest main body, the analysis of the coupling relationship among the multi-interest main body can be further expressed as follows:
Figure BDA00023919921100000410
Figure BDA00023919921100000411
Figure BDA00023919921100000412
in the middle of the above formula, the compound has the following structure,
Figure BDA00023919921100000413
the total electric load demand at the time t of the first day;
Figure BDA00023919921100000414
the total heat load demand at the time t of the first day;
Figure BDA00023919921100000415
storing power for the heat storage device in real time;
Figure BDA00023919921100000416
purchasing electric power from the side of the power grid;
Figure BDA00023919921100000417
outputting power for the distributed power supply;
Figure BDA00023919921100000418
is the total demand response;
Figure BDA00023919921100000419
the method comprises the following steps of (1) generating electric quantity and heating quantity for a cogeneration unit at t time of the first day;
Figure BDA00023919921100000420
the heat is generated by the gas boiler;
Figure BDA00023919921100000421
outputting the total electric quantity for the concentrator;
Figure BDA0002391992110000051
outputting the total heat for the concentrator;
Figure BDA0002391992110000052
inputting the electric quantity of a power grid for the concentrator;
Figure BDA0002391992110000053
converting heat into electricity;
Figure BDA0002391992110000054
generating heat for a steam boiler; mu.sCHP,eleμGB,thAnd muCHP,thThe electric, electric and gas-heat conversion efficiencies are respectively 0.3, 0.4 and 0.9.
The beneficial technical effects of the invention are as follows: (1) the method comprises the steps of constructing an internal coupling interaction model of a multi-benefit main body based on an energy hub, establishing a definite coupling matrix, quantifying the energy flow of a power connecting line and a natural gas pipeline between a system and a main network, analyzing the variation trend of electric heating loads respectively, and obtaining an optimal conclusion under the constraint of an objective function.
(2) An electricity price type demand response model is introduced, an electricity price mechanism is fully responded, distributed energy storage is matched, a typical region of a power distribution network is taken as an example, and the specific influence of a market mechanism on an electric heating load under the coupling condition of a multi-benefit subject is researched.
Drawings
FIG. 1 is a block diagram of a multi-stakeholder hub architecture according to an embodiment of the invention.
Fig. 2 is a schematic diagram illustrating an electrical load variation of the hub according to an embodiment of the present invention.
Fig. 3 is a schematic diagram illustrating a thermal load variation of the hub according to an embodiment of the present invention.
Fig. 4 is a schematic diagram showing the change of the purchased power under two scenarios according to the embodiment of the present invention.
Fig. 5 is a schematic diagram showing the change of the gas purchase amount under two scenarios according to the embodiment of the invention.
Fig. 6 is a schematic diagram showing the electrical load change before and after the PBDR is implemented in the hub 1 according to the embodiment of the present invention.
Fig. 7 is a schematic diagram showing the electrical load change before and after the PBDR implementation of the hub 2 according to the embodiment of the present invention.
Fig. 8 is a schematic diagram showing the electrical load change before and after the PBDR implementation of the hub 3 according to the embodiment of the present invention.
Fig. 9 shows a probability density diagram of photovoltaic contribution for an embodiment of the present invention.
Fig. 10 shows a probability distribution diagram of photovoltaic contribution for an embodiment of the present invention.
FIG. 11 is a probability density diagram of fan output according to an embodiment of the present invention.
FIG. 12 is a schematic illustration of a probability distribution of fan output according to an embodiment of the present invention.
FIG. 13 is a flowchart illustrating the overall steps of an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to fig. 1 to 13 of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example (b):
as shown in fig. 13, the modeling method for trading behavior of multi-benefit agent in market based on internal interactive collaboration comprises the following steps:
step 1: constructing an energy hub basic framework of a market multi-interest main body based on the input and output basic principle of the hub;
step 2: respectively carrying out independent modeling analysis on the multi-benefit main bodies interactively cooperated in the energy hub in the step 1, namely carrying out modeling analysis on an electric benefit main body, a gas benefit main body and a cogeneration benefit main body;
and step 3: establishing a coupling matrix to represent the incidence relation among all the benefit bodies in the step 2;
and 4, step 4: estimating the photovoltaic and fan output according to historical data, measuring the tie line power of the system and calculating the net exchange power; measuring gas flow data of a natural gas pipeline, and setting a target function to calculate the gas and electricity acquisition cost; analyzing the change trend of the thermoelectric load according to the response change; and (4) analyzing and calculating the obtained optimal conclusion through the coupling matrix in the step (3), and guiding the reasonable distribution of the multi-benefit subjects according to the optimal conclusion.
The energy hub is not only a connection center between various energy forms and different types of loads, but also comprises the transmission, conversion and storage of energy of different forms, and has good control characteristics. The problem of weak correlation between the power grid and the micro-grid is solved to a certain extent. As a virtual control platform, the system can monitor the running states of the micro-grid and the power distribution network in real time, respond to the change of distributed energy output and optimize various load requirements in time. The diversity of distributed energy sources dictates that multi-stakeholders must include a variety of energy supply equipment. The basic framework of a multi-benefit-subject hub is designed by integrating benefit subjects such as electric power, heat, gas and the like into research objects, and is shown in figure 1.
For a single multi-benefit agent hub (MSH), the multi-port network should include photovoltaic, wind turbines, cogeneration units, gas boilers, transformers, and energy storage systems, among other functional units. For a single input-output energy conversion device, the coupling relationship between the input and output is typically described by a coupling matrix in the following specific form:
Figure BDA0002391992110000061
Figure BDA0002391992110000062
in the formula (1), the vector
Figure BDA0002391992110000063
And
Figure BDA0002391992110000064
representing inputs and outputs of different energies of a single hub, matrix lambdamnUsed as a coupling matrix describing the relationship between the input and the output, where the value is related to the conversion efficiency of the power supply device.
A. Electric benefits subject
For the power benefit agent, the comprehensive scheduling and load demand response of photovoltaic, wind power generation and energy storage equipment are mainly considered. In order to solve the uncertainty of photovoltaic and fan output, a stochastic programming method is selected, and the uncertainty of a distributed power supply is described by using a probability density distribution function (CDF), wherein the photovoltaic output obeys beta distribution, and the fan output obeys wibu distribution of two parameters. The detailed figures are shown in detail in fig. 9-12.
(1) Energy storage model
The energy storage device can be flexibly called in the scheduling period, and the dynamic model can be expressed as follows:
Figure BDA0002391992110000071
Figure BDA0002391992110000072
in the formula (3), x is a binary variable and represents the charge-discharge state of the energy storage device at the ith day t, 1 is in a working state, 0 is in a non-working state,
Figure BDA0002391992110000073
in order to charge and discharge power in real time,
Figure BDA0002391992110000074
is a charge-discharge power limit; in the formula (4), Ei,0,Ei,24Representing the stored energy at the beginning and end of the energy storage device; es,t,jReal-time charge;
Figure BDA0002391992110000075
limiting the real-time charge capacity of the equipment;
Figure BDA0002391992110000076
the maximum charge capacity of the equipment; ei,tRepresenting the electric quantity stored by the energy storage equipment at adjacent moments, α and β are charge and discharge efficiencies, and are both set to be 0.9;
for the load demand response, the PBDR is selected as a research object in the load demand response project, the multi-interest subject is built, and the user is guided to carry out load distribution by responding the time-of-use electricity price signal in real time; the dual electricity price mechanism of the time-of-use electricity price and the real-time electricity price is selected, a user responds to the time-of-use electricity price and can participate in the PBDR project according to the real-time electricity price, and load demand calculation after PBDR is realized can be represented as follows:
Figure BDA0002391992110000077
in the formula:
Figure BDA0002391992110000078
in order to implement the load after PBDR,
Figure BDA0002391992110000079
initial load before PBDR implementation; t is a total scheduling period, and s and T are sub-periods; e.g. of the typestThe value is the price elastic coefficient, and s-t is the self elastic coefficient, otherwise, the value is the cross elastic coefficient;
Figure BDA00023919921100000710
to implement the real-time electricity rates after PBDR,
Figure BDA00023919921100000711
initial time-of-use electricity prices before PBDR is implemented; the real-time electricity price can be calculated according to the following formula:
Figure BDA00023919921100000712
in the formula:
Figure BDA00023919921100000713
is the total load demand for a day;
Figure BDA00023919921100000714
is the average load of a day.
B. Gas benefit agent
The gas boiler benefit body (GBS) is mainly typified by a gas boiler. They usually only switch gas and warm air. The model is represented by the coupling relation between the heat supply power and the gas consumption in unit time:
Figure BDA0002391992110000081
in the formula, muGB,thRepresents the electrical conversion efficiency;
Figure BDA0002391992110000082
the heat supply power of the gas-fired boiler is increased,
Figure BDA0002391992110000083
is the amount of natural gas consumed per unit time.
C. Combined heat and power benefit agent
For the benefit subject of cogeneration, natural gas is used as the main energy of the thermal power generating unit, and has the unique advantages of high power generation efficiency, comprehensive utilization of energy and the like. The power supply and heat supply are expressed as:
Figure BDA0002391992110000084
in the formula:
Figure BDA0002391992110000085
the electric quantity and the heating quantity of the cogeneration unit generated at the first day t,
Figure BDA0002391992110000086
is the consumption of natural gas per unit time, muCHP,ele,μCHP,thThe gas-electricity conversion efficiency and the gas-heat conversion efficiency are respectively.
D. Hub internal coupling relationship
By modeling the coupling relationship among the internal devices of the multi-interest main body, the analysis of the coupling relationship among the multi-interest main body can be further expressed as follows:
Figure BDA0002391992110000087
Figure BDA0002391992110000088
Figure BDA0002391992110000089
in the upper type,
Figure BDA00023919921100000810
The total electric load demand at the time t of the first day;
Figure BDA00023919921100000811
the total heat load demand at the time t of the first day;
Figure BDA00023919921100000812
storing power for the heat storage device in real time;
Figure BDA00023919921100000813
purchasing electric power from the side of the power grid;
Figure BDA00023919921100000814
outputting power for the distributed power supply;
Figure BDA00023919921100000815
is the total demand response;
Figure BDA00023919921100000816
the method comprises the following steps of (1) generating electric quantity and heating quantity for a cogeneration unit at t time of the first day;
Figure BDA00023919921100000817
the heat is generated by the gas boiler;
Figure BDA00023919921100000818
outputting the total electric quantity for the concentrator;
Figure BDA00023919921100000819
outputting the total heat for the concentrator;
Figure BDA00023919921100000820
inputting the electric quantity of a power grid for the concentrator;
Figure BDA00023919921100000821
converting heat into electricity;
Figure BDA00023919921100000822
generating heat for a steam boiler; mu.sCHP,eleμGB,thAnd muCHP,thThe electric, electric and gas-heat conversion efficiencies are respectively 0.3, 0.4 and 0.9.
Simulation analysis
Three hub hubs are taken as an example, the rationality of the scheme is verified, three different parks are included, and the influence of multi-dimensional characteristics of multi-benefit subjects on gas and power purchasing optimization under different conditions is discussed. The scheduling period is about 25 days, i is 25, T is 24 hours, the time step is 1 hour, and the internal parameters of each hub are identical. The thermal load and electrical load values of different hubs are shown in the figure, and a probability density function curve of an output is obtained by inputting historical data. See figures 2-4 for details. In addition, other power limits are shown in table 2, and the time-of-use electricity rates are shown in table 3. For better results, two basic scenarios were set for comparison. Firstly, the system does not contain a demand response item, and only carries out gas-electricity consumption analysis without relating to the movement of user-side resources; the other scene is to introduce a demand corresponding item, because the user's electricity utilization behavior can be adjusted through the real-time response to the electricity price signal, the user side resource can be fully adjusted, and the peak-load and valley-load adjusting effect is better under the stimulation of the time-of-use electricity price. The resulting pattern can be contrasted sharply.
Figure BDA0002391992110000091
TABLE 2 electric quantity and Natural gas quantity Transmission constraints
Figure BDA0002391992110000092
TABLE 3 time-of-use price for buying and selling electricity
For case 1, which contains all of the stakeholders in the hub, but no demand response items, the CDF curve for the gas and electricity purchases in this case is shown in fig. 5. As can be clearly seen from the figure, the purchase amount of the natural gas is between 0kW and 9000kW, the purchase amount reaches a peak value around 7500kW, and the probability value reaches 1. Meanwhile, the electricity purchasing power is between 0 and 5000kw, the electricity purchasing quantity reaches a peak value at about 4500kw, and the probability value is 1. Due to the uncertainty of the distributed power generation, when the photovoltaic and the fan can not meet the basic load requirement of urban domestic electricity, the requirements of electricity and heat load can be adjusted by purchasing electricity and gas, and the constraint of the limit of power transmission power needs to be met.
For case 2, we can see that the rising trend of the CDF curve is similar to case 1. But the amount of electricity and gas purchased is significantly reduced due to the implementation of the PBDR project. This is because the PBDR reduces the power purchase from the distribution network in the garden to some extent by actively guiding the power utilization behavior of the user, and has good peak load regulation and valley filling capability. On the basis of case 1, the gas supply pressure of the natural gas system is also reduced. The improvement of the energy supply proportion of the power system is beneficial to improving the economic operation level of the whole system and relieving the power supply pressure of the system.
FIGS. 6-8 show the change in internal load capacity of a single hub before and after the implementation of a PBDR project. Therefore, after the PBDR project is implemented, the load curve has good variation trend in different peak-valley periods, and the good characteristic that PBDR has peak clipping and valley filling is reflected. Different hub load values are different, but the curves have similar basic trends, which shows that the advantages of PBDR in system operation are particularly obvious. Further, since the reduction in the electric load affects the purchase amount of natural gas to some extent, the dependence on natural gas is also reduced, and thus the cost influence due to the internal coupling characteristics of the multi-benefit agent is easily obtained.
The gas and electricity procurement costs for the hub are shown in table 1 below. It can be seen by comparison that the gas and electricity procurement costs are higher for case 1 than for case 2. The main reason is that the implementation of the PBDR project actively guides the electricity utilization behavior of users according to the time-of-use electricity price, and reduces the running cost of the whole system and the power supply pressure of each concentrator.
Figure BDA0002391992110000101
TABLE 1 comparison of cost calculations for each part under different conditions
The method mainly researches multi-dimensional feature modeling of multi-interest main bodies in a park, and the multi-dimensional feature modeling comprises multi-period operation characteristics, internal coupling characteristics and space characteristics, provides the concept of a multi-interest related person concentrator system, and constructs a multi-energy coupling model in a concentrator. How multi-benefit agents interact with the power distribution network through internal interactive cooperation is analyzed. Finally, the effectiveness of the trading behavior of the distribution network multi-interest main body based on internal interactive cooperation is illustrated through a calculation example. Meanwhile, the PBDR project is used as a comparison object, so that the implementation of the PBDR can improve the economy of system operation and relieve the power supply pressure of a distribution network.
In the description of the present invention, it is to be understood that the terms "counterclockwise", "clockwise", "longitudinal", "lateral", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc., indicate orientations or positional relationships based on those shown in the drawings, and are used for convenience of description only, and do not indicate or imply that the referenced devices or elements must have a particular orientation, be constructed and operated in a particular orientation, and thus, are not to be considered as limiting.

Claims (6)

1. The market multi-interest main body trading behavior modeling method based on internal interactive collaboration is characterized by comprising the following steps of:
step 1: constructing an energy hub basic framework of a market multi-interest main body based on the input and output basic principle of the hub;
step 2: respectively carrying out independent modeling analysis on the multi-benefit main bodies interactively cooperated in the energy hub in the step 1, namely carrying out modeling analysis on an electric benefit main body, a gas benefit main body and a cogeneration benefit main body;
and step 3: establishing a coupling matrix to represent the incidence relation among all the benefit bodies in the step 2;
and 4, step 4: estimating the photovoltaic and fan output according to historical data, measuring the tie line power of the system and calculating the net exchange power; measuring gas flow data of a natural gas pipeline, and setting a target function to calculate the gas and electricity acquisition cost; analyzing the change trend of the thermoelectric load according to the response change; and (4) analyzing and calculating the obtained optimal conclusion through the coupling matrix in the step (3), and guiding the reasonable distribution of the multi-benefit subjects according to the optimal conclusion.
2. The internal interactive collaboration based modeling method for trading behavior of multi-interest main body in market according to claim 1, wherein in step 1, the energy hub basic framework is as follows:
Figure FDA0002391992100000011
Figure FDA0002391992100000012
in the formula (1), the vector
Figure FDA0002391992100000013
And
Figure FDA0002391992100000014
representing inputs and outputs of different energies of a single hub, matrix lambdamnUsed as a coupling matrix describing the relationship between the input and the output, where the value is related to the conversion efficiency of the power supply device.
3. The market multi-interest-principal trading behavior modeling method based on internal interactive collaboration as claimed in claim 2, wherein the modeling analysis of the electric interest principal in the step 2 is specifically as follows:
the electric interest main body comprises photovoltaic, wind power generation, comprehensive dispatching of energy storage equipment and load demand response;
for photovoltaic and wind power generation, a random planning method is selected, and the uncertainty of a distributed power supply is described by using a probability density distribution function, wherein the photovoltaic output obeys beta distribution, and the fan output obeys two-parameter Weibu distribution;
for the comprehensive scheduling of the energy storage device, the energy storage device can be flexibly called in a scheduling period, and a dynamic model of the energy storage device can be expressed as follows:
Figure FDA0002391992100000015
Figure FDA0002391992100000021
in the formula (3), x is a binary variable and represents the charge-discharge state of the energy storage device at the ith day t, 1 is in a working state, 0 is in a non-working state,
Figure FDA0002391992100000022
in order to charge and discharge power in real time,
Figure FDA0002391992100000023
is a charge-discharge power limit; in the formula (4), Ei,0,Ei,24Representing the stored energy at the beginning and end of the energy storage device; es,t,jReal-time charge;
Figure FDA0002391992100000024
limiting the real-time charge capacity of the equipment;
Figure FDA0002391992100000025
the maximum charge capacity of the equipment; ei,tRepresenting the electric quantity stored by the energy storage equipment at adjacent moments, α and β are charge and discharge efficiencies, and are both set to be 0.9;
for the load demand response, the PBDR is selected as a research object in the load demand response project, the multi-interest subject is built, and the user is guided to carry out load distribution by responding the time-of-use electricity price signal in real time; the dual electricity price mechanism of the time-of-use electricity price and the real-time electricity price is selected, a user responds to the time-of-use electricity price and can participate in the PBDR project according to the real-time electricity price, and load demand calculation after PBDR is realized can be represented as follows:
Figure FDA0002391992100000026
in the formula:
Figure FDA0002391992100000027
in order to implement the load after PBDR,
Figure FDA0002391992100000028
initial load before PBDR implementation; t is a total scheduling period, and s and T are sub-periods; e.g. of the typestThe value is the price elastic coefficient, and s-t is the self elastic coefficient, otherwise, the value is the cross elastic coefficient;
Figure FDA0002391992100000029
to implement the real-time electricity rates after PBDR,
Figure FDA00023919921000000210
initial time-of-use electricity prices before PBDR is implemented; the real-time electricity price can be calculated according to the following formula:
Figure FDA00023919921000000211
in the formula:
Figure FDA00023919921000000212
is the total load demand for a day;
Figure FDA00023919921000000213
is the average load of a day.
4. The internal interactive collaboration-based market multi-interest entity trading behavior modeling method according to claim 3, wherein the modeling analysis on the gas interest entity in the step 2 specifically comprises the following steps:
the interest subject of the gas boiler takes the gas boiler as a considered object, and the model of the interest subject of the gas boiler is represented by the coupling relation of the heating power and the gas consumption in unit time:
Figure FDA0002391992100000031
in the formula, muGB,thRepresents the electrical conversion efficiency;
Figure FDA0002391992100000032
the heat supply power of the gas-fired boiler is increased,
Figure FDA0002391992100000033
is the amount of natural gas consumed per unit time.
5. The internal interactive collaboration-based market multi-interest entity trading behavior modeling method according to claim 4, wherein the modeling analysis on the cogeneration interest entity in the step 2 is specifically as follows:
the combined heat and power benefit agent takes natural gas as a consideration object, and the power supply and heat supply quantity of the combined heat and power benefit agent is expressed as follows:
Figure FDA0002391992100000034
in the formula:
Figure FDA0002391992100000035
the electric quantity and the heating quantity of the cogeneration unit generated at the first day t,
Figure FDA0002391992100000036
is the consumption of natural gas per unit time, muCHP,ele,μCHP,thThe gas-electricity conversion efficiency and the gas-heat conversion efficiency are respectively.
6. The internal interactive collaboration-based market multi-interest entity trading behavior modeling method according to claim 5, wherein the coupling matrix in the step 3 is specifically:
by modeling the coupling relationship among the internal devices of the multi-interest main body, the analysis of the coupling relationship among the multi-interest main body can be further expressed as follows:
Figure FDA0002391992100000037
Figure FDA0002391992100000038
Figure FDA0002391992100000039
in the middle of the above formula, the compound has the following structure,
Figure FDA00023919921000000310
the total electric load demand at the time t of the first day;
Figure FDA00023919921000000311
the total heat load demand at the time t of the first day;
Figure FDA00023919921000000312
storing power for the heat storage device in real time;
Figure FDA00023919921000000313
purchasing electric power from the side of the power grid;
Figure FDA00023919921000000314
outputting power for the distributed power supply;
Figure FDA00023919921000000315
is the total demand response;
Figure FDA00023919921000000316
the method comprises the following steps of (1) generating electric quantity and heating quantity for a cogeneration unit at t time of the first day;
Figure FDA00023919921000000317
the heat is generated by the gas boiler;
Figure FDA00023919921000000318
outputting the total electric quantity for the concentrator;
Figure FDA00023919921000000319
outputting the total heat for the concentrator;
Figure FDA00023919921000000320
inputting the electric quantity of a power grid for the concentrator;
Figure FDA00023919921000000321
converting heat into electricity;
Figure FDA00023919921000000322
generating heat for a steam boiler; mu.sCHP,eleμGB,thAnd muCHP,thThe electric, electric and gas-heat conversion efficiencies are respectively 0.3, 0.4 and 0.9.
CN202010117662.2A 2020-02-25 2020-02-25 Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation Pending CN111260158A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010117662.2A CN111260158A (en) 2020-02-25 2020-02-25 Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010117662.2A CN111260158A (en) 2020-02-25 2020-02-25 Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation

Publications (1)

Publication Number Publication Date
CN111260158A true CN111260158A (en) 2020-06-09

Family

ID=70952762

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010117662.2A Pending CN111260158A (en) 2020-02-25 2020-02-25 Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation

Country Status (1)

Country Link
CN (1) CN111260158A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111861551A (en) * 2020-07-08 2020-10-30 海宁同兴智慧园区管理有限公司 Intelligent monitoring and management system for electricity consumption of park
CN117131982A (en) * 2023-08-21 2023-11-28 四川大学 User side energy optimization method for guaranteeing reliable load supply under extremely high temperature event

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106532709A (en) * 2016-10-25 2017-03-22 天津大学 Distribution-grid-reconstruction-based optimal power flow calculation method for integrated communication energy system
CN107611966A (en) * 2017-09-20 2018-01-19 天津大学 A kind of active power distribution network evaluation of power supply capability method for considering difference reliability
CN110245810A (en) * 2019-07-09 2019-09-17 国网山东省电力公司聊城供电公司 A kind of Regional Energy internet coordination optimization operation method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106532709A (en) * 2016-10-25 2017-03-22 天津大学 Distribution-grid-reconstruction-based optimal power flow calculation method for integrated communication energy system
CN107611966A (en) * 2017-09-20 2018-01-19 天津大学 A kind of active power distribution network evaluation of power supply capability method for considering difference reliability
CN110245810A (en) * 2019-07-09 2019-09-17 国网山东省电力公司聊城供电公司 A kind of Regional Energy internet coordination optimization operation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
FRANZISKA ADAMEK 等: "On Decisive Storage Parameters for Minimizing Energy Supply Costs in Multicarrier Energy Systems", 《IEEE TRANSACTIONS ON SUSTAINABLE ENERGY》 *
倪伟 等: "基于机会约束规划的能源集线器系统气电购置优化建模", 《电网技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111861551A (en) * 2020-07-08 2020-10-30 海宁同兴智慧园区管理有限公司 Intelligent monitoring and management system for electricity consumption of park
CN111861551B (en) * 2020-07-08 2023-12-05 海宁同兴智慧园区管理有限公司 Intelligent monitoring management system for power consumption of park
CN117131982A (en) * 2023-08-21 2023-11-28 四川大学 User side energy optimization method for guaranteeing reliable load supply under extremely high temperature event
CN117131982B (en) * 2023-08-21 2024-04-05 四川大学 User side energy optimization method for guaranteeing reliable load supply under extremely high temperature event

Similar Documents

Publication Publication Date Title
Khalili et al. Optimal performance of microgrid in the presence of demand response exchange: A stochastic multi-objective model
Tan et al. Multi-scenario operation optimization model for park integrated energy system based on multi-energy demand response
WO2022166114A1 (en) Electro-hydrogen energy system scheduling method thst considers flexible hydrogen requirements
Ju et al. Multi-objective electro-thermal coupling scheduling model for a hybrid energy system comprising wind power plant, conventional gas turbine, and regenerative electric boiler, considering uncertainty and demand response
CN111260158A (en) Market multi-interest main body transaction behavior modeling method based on internal interactive cooperation
CN116911533B (en) Multi-microgrid energy sharing method for regional comprehensive energy system
CN111832807B (en) Multi-micro-grid coordination optimization scheduling method considering load characteristics and demand response
Deng et al. Research on economic dispatch of integrated energy system based on improved krill swarm algorithm
CN112883571A (en) Comprehensive energy system planning method considering dynamic multi-energy demand response
Liu et al. Optimal dispatch strategy of virtual power plants using potential game theory
Zhang et al. Bi-layer economic scheduling for integrated energy system based on source-load coordinated carbon reduction
Wu et al. Dynamic pricing and energy management of hydrogen-based integrated energy service provider considering integrated demand response with a bi-level approach
Yunshou et al. A collaborative demand-controlled operation strategy for a multi-energy system
CN115879613A (en) Agricultural industry park load optimization scheduling method considering carbon transaction
CN111967639A (en) Power distribution network optimization method and system considering wind-solar power output and price demand response
CN115099550A (en) Market interaction profit model calculation device considering flexibility of comprehensive energy system
CN116957139B (en) Multi-comprehensive-energy microgrid optimal operation method and system considering carbon transaction among microgrids
CN114444797A (en) Community energy Internet double-layer distributed interactive optimization method considering carbon emission reduction
Khajehzadeh et al. Novel control framework for optimal scheduling in microgrid with demand response support under contingency events
Zhuang et al. Flexibility exploitation with nonlinear integrated demand response for multi-energy system against load estimation mistake
Ma et al. Low-carbon economic operation of energy hub integrated with linearization model and nodal energy-carbon price
CN111754076B (en) Method and equipment suitable for evaluating electric wholesale market mode
CN116128163B (en) Comprehensive energy optimization method and device considering green hydrogen production and storage and user satisfaction
CN112001613A (en) Benefit distribution strategy based on multi-micro-energy-source network system cooperative game decision mechanism
CN116341268A (en) Cold-hot electric coupling type microgrid group three-layer main body distributed optimization method

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200609