CN109510224B - Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy - Google Patents

Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy Download PDF

Info

Publication number
CN109510224B
CN109510224B CN201811368811.1A CN201811368811A CN109510224B CN 109510224 B CN109510224 B CN 109510224B CN 201811368811 A CN201811368811 A CN 201811368811A CN 109510224 B CN109510224 B CN 109510224B
Authority
CN
China
Prior art keywords
power
photovoltaic
energy storage
energy
network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811368811.1A
Other languages
Chinese (zh)
Other versions
CN109510224A (en
Inventor
刘学智
陈思捷
严正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201811368811.1A priority Critical patent/CN109510224B/en
Publication of CN109510224A publication Critical patent/CN109510224A/en
Application granted granted Critical
Publication of CN109510224B publication Critical patent/CN109510224B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin

Landscapes

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

Abstract

The invention relates to a capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy, which comprises the following steps: (1) acquiring basic data of a distributed energy system and a thermal power and power load prediction data curve accurate to the hour of a year; (2) configuring the capacity of energy conversion equipment in the distributed energy system by taking the minimum annual total investment operation cost of the distributed energy system as an optimization target; (3) determining an optimized output curve of the distributed energy system; (4) determining a power net load curve; (5) collecting basic data of a photovoltaic energy storage system; (6) and optimally configuring the battery energy storage rated power, the battery energy storage rated capacity and the charge and discharge strategy in the photovoltaic energy storage system by taking the minimum annual total investment operation cost of the photovoltaic energy storage system as an optimization target according to the basic data and the electric net load curve of the photovoltaic energy storage system. Compared with the prior art, the method is rapid and comprehensive, and the result is accurate and reliable.

Description

Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy
Technical Field
The invention relates to a capacity configuration and operation optimization method for a comprehensive energy network, in particular to a capacity configuration and operation optimization method combining photovoltaic energy storage and distributed energy.
Background
The comprehensive energy system considers the synergistic effect of various energy sources such as cold/heat/electricity/gas and the like, meets the energy requirements of social energy terminal users on power supply, heat supply, cold supply and the like by integrating energy conversion and storage among the various energy network systems, plans and operates the various energy network systems as a whole, and improves the operation efficiency of the system so as to achieve the purposes of energy conservation and emission reduction. Under the framework of a comprehensive energy system, the cooperative optimization of a plurality of energy systems of electric power, heat and gas is realized by further introducing high-capacity energy storage, the large space-time range optimal configuration capacity of the energy system is improved, and the problems of renewable energy consumption, peak regulation and the like can be effectively solved.
In the aspect of distributed energy conversion equipment, as different geographical climate resource conditions have certain complexity, unified planning for integrating and accommodating various conditions is lacked, and various energy storage and conversion equipment are systematically integrated. The heat pump heating is widely used in Japan, northern Europe and Germany use high proportion of regional heating, and great amount of gas boilers are used for direct heating in the UK by means of developed natural gas pipe network infrastructures. Study doctor academic paper of zhou ji of qinghua university: in the modeling and optimization of the distributed energy system with multi-energy coordination, a set of energy system equipment optimization design model based on superstructure linearization modeling is provided for the distributed energy system with multi-energy coordination, but energy flow calculation (power flow, thermal flow and gas flow) of a multi-energy network is not considered. A double-layer optimization planning design method of a combined cooling heating and power micro-grid system is disclosed in the patent application No.: CN201310661953.8 discloses a two-stage modeling and planning method for coupling characteristics of planning design and operation optimization of a micro-grid, which can realize interactive optimization of an outer-layer device type, a capacity optimization module and an inner-layer operation strategy optimization module, but does not consider modeling of a thermal power grid and a gas grid, and the double-layer planning and solving speed is very slow, often reaching several hours or more.
In the aspect of configuration and operation of the photovoltaic-energy storage system, the capacity design and operation strategy of battery energy storage are key. Too small a capacity does not effectively take up photovoltaic power, and too large a capacity incurs too high an investment cost for the battery, so there is an optimum balance point between the capacity and price of the battery. The benefits of photovoltaic and battery energy storage coupling come from: 1. and 2, the difference between the grid electricity price when the redundant photovoltaic power is stored and the photovoltaic outgoing electricity price when the redundant photovoltaic power is stored and the grid electricity price when the part of power is released. The heavy use of energy conversion equipment such as heat pumps, CHP, air conditioners, electric vehicles and the like affects the power curve of the load, thereby affecting the configuration and operation of the photovoltaic energy storage system. At present, a great deal of research is carried out on photovoltaic-energy storage system configuration and operation scheduling at home and abroad, and the problem of maximum profit is discussed from an individual perspective by focusing on an individual user taking a house or a building as a unit. But research on the optimal configuration and operation of the energy conversion equipment and the photovoltaic energy storage system is not considered integrally.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a capacity configuration and operation optimization method combining photovoltaic energy storage and distributed energy sources.
The purpose of the invention can be realized by the following technical scheme:
a capacity configuration and operation optimization method combining photovoltaic energy storage and distributed energy sources comprises the following steps:
(1) acquiring basic data of a distributed energy system and a thermal power and power load prediction data curve accurate to the hour of a year;
(2) according to basic data of the distributed energy system and a thermal power and power load prediction data curve, the capacity of energy conversion equipment in the distributed energy system is configured by taking the minimum annual total investment and operation cost of the distributed energy system as an optimization target;
(3) determining an optimized output curve of the distributed energy system according to the capacity configuration result of the energy conversion equipment;
(4) determining an electric net load curve according to the electric load prediction data curve and the distributed energy system operation optimization output curve;
(5) collecting basic data of a photovoltaic energy storage system;
(6) and optimally configuring the battery energy storage rated power, the battery energy storage rated capacity and the charge and discharge strategy in the photovoltaic energy storage system by taking the minimum annual total investment operation cost of the photovoltaic energy storage system as an optimization target according to the basic data and the electric net load curve of the photovoltaic energy storage system.
The step (2) is specifically as follows:
(21) establishing an optimization objective function of the distributed energy system:
Min Ctotal-DER=LFn·CCapex-DER+COpex-DER
wherein, Ctotal-DERFor the annual total investment and operation cost of the distributed energy system, CCapex-DERInitial total investment cost of energy conversion equipment in distributed energy system,COpex-DERAnnual energy conversion equipment operating charge, LF, for a distributed energy systemnIn order to be a coefficient of capital recovery,
Figure GDA0003248035430000031
d is the discount rate, and n is the service life of the energy conversion equipment;
(22) establishing a constraint function of the distributed energy system, wherein the constraint function comprises an equality constraint function and an inequality constraint function, and the equality constraint function is as follows:
Figure GDA0003248035430000032
in the formula, PiIs the active power of the power of node i, ViIs the voltage of node i, VjIs the voltage of node j, NeIs the number of nodes, G, of the power systemijIs the conductance of line ij, BijIs the susceptance, θ, of line ijijIs the voltage angle difference of nodes i and j, QiElectric reactive power, C, for node ipIs the specific heat capacity of water, AhIs a heat power network association matrix and is characterized in that,
Figure GDA00032480354300000311
for mass flow, TsFor temperature of the water supply, ToFor the outlet water temperature, phi is the thermal power vector consumed or provided by each thermal node, BhIs a heat power network loop incidence matrix, KhIs the coefficient of resistance of the heat distribution network pipe, CsAs a water supply network coefficient matrix, bsIs a constant vector, CrAs a coefficient matrix of the backwater network, TrTo the return water temperature, brIs a constant vector, AgIs a gas network incidence matrix, vgIs node air pressure, vqTo node the gas flow rate, BgIs a gas network loop incidence matrix, KgIs the resistance coefficient of the gas network pipeline, and k is an exponential constant;
the inequality constraint function includes:
(a) and (3) limiting the upper limit and the lower limit of the capacity of the energy conversion equipment:
Figure GDA0003248035430000033
Figure GDA0003248035430000034
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000035
for the value of the power active power provided by the conversion device at peak load,
Figure GDA0003248035430000036
the lower limit of the active power of the electricity supplied to the conversion equipment,
Figure GDA0003248035430000037
the upper limit of the electrical active power provided to the conversion equipment,
Figure GDA0003248035430000038
for converting the thermal power supplied by the device at load peaks,
Figure GDA0003248035430000039
the lower limit of the thermal power is,
Figure GDA00032480354300000310
is the upper limit of thermodynamic power;
(b) the power grid inequality constrains:
Vimin≤Vi≤Vimax,i=1,…,Ne
Figure GDA0003248035430000041
Figure GDA0003248035430000042
Figure GDA0003248035430000043
in the formula, ViIs the voltage of node i, ViminLower voltage limit, V, of node iimaxIs the upper voltage limit of node i, NeIs the number of nodes, PgeniFor the active power output of the ith generator,
Figure GDA0003248035430000044
the lower active power output limit of the ith generator,
Figure GDA0003248035430000045
upper limit of active power output, Q, of the ith generatorgeniIs the reactive power output of the ith generator,
Figure GDA0003248035430000046
the lower limit of reactive power output of the ith generator,
Figure GDA0003248035430000047
is the upper limit of reactive power output of the ith generator, NgeNumber of generators, SkIs the electrical power of the kth branch,
Figure GDA0003248035430000048
the lower limit of the electric power of the kth branch,
Figure GDA0003248035430000049
is the upper electric power limit of the kth branch, NleThe number of branches;
(c) the thermal power grid inequality constrains:
Figure GDA00032480354300000410
Ts_min≤Ts≤Ts_max
Tr_min≤Tr≤Tr_max
in the formula (I), the compound is shown in the specification,
Figure GDA00032480354300000411
the flow rate of each pipeline of the heating power network,
Figure GDA00032480354300000412
the lower limit of the flow of each pipeline of the heating power network,
Figure GDA00032480354300000413
is the upper limit of the flow of each pipeline of the heating power network, TsSupply of water temperature, T, to each node of the heating networks_minLower limit of water supply temperature, T, for each node of heating power networks_maxUpper limit of water supply temperature, T, for each node of heating power networkrFor the return water temperature, T, of each node of the heating power networkr_minThe lower limit of the return water temperature of each node of the heating power network, Tr_maxThe upper limit of the return water temperature of each node of the heating power network is set;
(d) the inequality constraint of the gas network is as follows:
pgmin≤pg≤pgmax
vg_min≤vg≤vg_max
in the formula, pgFor the gas pressure, p, at each node of the gas networkgminIs the lower limit of gas pressure of each node of the gas network, pgmaxIs the upper limit of gas pressure of each node of the gas network, vgIs the gas flow of each pipeline of the gas network, vg_minIs the lower limit of the gas flow of each pipeline of the gas network, vg_maxThe upper limit of the gas flow of each pipeline of the gas network;
(23) and optimizing and solving to obtain the capacity of the energy conversion equipment in the distributed energy system.
Initial total investment cost C of energy conversion equipment in distributed energy systemCapex-DERThe method specifically comprises the following steps:
Figure GDA00032480354300000414
wherein the content of the first and second substances,
Figure GDA00032480354300000415
in order to convert the initial investment cost of the equipment i,
Figure GDA00032480354300000416
the electrical power provided by the inverter device i at peak load,
Figure GDA00032480354300000417
the V-shaped index represents the operator "OR" for the thermal power supplied by the switching device i at load peaks.
Annual energy conversion equipment operating charge C in distributed energy systemOpex-DERThe method specifically comprises the following steps:
Figure GDA0003248035430000051
wherein, CeFor purchase of electricity price, Pimport(t) buying electric power to the upper grid for the t hour of the year, CgAs gas price, vgtotal(t) Total gas consumption Rate in the t hour of the year, CcarbonIs the carbon price, ξeIs the carbon emission intensity of the powergIs the carbon emission intensity of natural gas, CO&MFor operating maintenance costs, T is the number of hours of a year.
The step (6) is specifically as follows:
(61) establishing an optimization objective function of the photovoltaic energy storage system:
Min Ctotal-PVB=CCapex-PVB+COpex-PVB
wherein, Ctotal-PVBThe total annual investment and operation cost of the photovoltaic energy storage system, CCapex-PVBFor the annual investment costs of the photovoltaic energy storage system, COpex-PVBThe annual running cost of the photovoltaic energy storage system is saved;
(62) establishing a constraint function of the photovoltaic energy storage system, comprising:
(a) charge and discharge power constraint:
Figure GDA0003248035430000052
Figure GDA0003248035430000053
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000054
charging power for the t hour on day d,
Figure GDA0003248035430000055
the discharge power at the t hour on day d,
Figure GDA0003248035430000056
d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24;
(b) energy balance constraint:
Figure GDA0003248035430000057
in the formula, EBESS(d,t)EBESS(t) is the amount of energy stored in the battery at the time of day d and t, EBESS(d,t-1)EBESS(t-1) is the electric quantity of the battery stored energy at the t-1 moment on the day d, etacEfficiency of charging the battery, ηdFor the battery discharge efficiency, Δ t is the time step;
(c) battery storage capacity constraint:
Figure GDA0003248035430000058
in the formula, EBESS(d, t) is the electric quantity stored by the battery at the time t on the day d,
Figure GDA0003248035430000059
storing energy rating, SOC, for a batteryminAllowed by ratio of residual capacity to capacity of batteryMinimum value, SOCmaxThe minimum value allowed by the ratio of the energy storage residual capacity to the capacity of the battery is obtained;
(63) and optimizing and solving to obtain the battery energy storage rated power, the battery energy storage rated capacity and the charging and discharging power at each moment in the photovoltaic energy storage system.
Annual investment cost C of photovoltaic energy storage systemCapex-PVBThe method specifically comprises the following steps:
Figure GDA0003248035430000061
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000062
the rated capacity of the energy storage of the battery,
Figure GDA0003248035430000063
the rated power is stored for the battery,
Figure GDA0003248035430000064
in order to obtain the rated output power of the photovoltaic,
Figure GDA0003248035430000065
the capacity investment cost per kWh for battery energy storage,
Figure GDA0003248035430000066
the power investment cost per kW for battery energy storage,
Figure GDA0003248035430000067
is the unit power investment cost of the photovoltaic.
Annual operating charge C of photovoltaic energy storage systemOpex-PVBThe method specifically comprises the following steps:
Figure GDA0003248035430000068
in the formula, PPV→load(d, t) is the power of the photovoltaic supply load at the tth day, PPV→bess(d, t) isPower of photovoltaic charging of battery at t hour on day d, PPV→grid(d, t) is the power of photovoltaic externally transmitted to a superior power grid in the tth hour on the day d, PPV(d, t) is the generated power of the photovoltaic,
Figure GDA0003248035430000069
charging power for the t hour on day d,
Figure GDA00032480354300000610
discharge power at day d, t hour, Ce(t) is the power grid electricity price at the moment t,
Figure GDA00032480354300000611
is the electricity price of the photovoltaic power generation outgoing power grid,
Figure GDA00032480354300000612
d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24 for the subsidy price per kW of photovoltaic power generation.
Compared with the prior art, the invention has the following advantages:
the invention provides a rapid, comprehensive and extensible capacity configuration and operation optimization method for photovoltaic energy storage and distributed energy combination, aiming at the mutual relation of energy storage and conversion equipment influence load curves, a whole system design and operation model of a comprehensive energy network is established, the whole optimization is realized, the influence of the deployment of a distributed energy system (DER) on net load change is calculated according to the configuration result of the conversion equipment in the DER, on the basis, the design and operation of a photovoltaic energy storage system (PVB) are researched, and economic incentive is provided for owners of a photovoltaic cell energy storage system.
Drawings
FIG. 1 is a block flow diagram of a method for capacity allocation and operation optimization of photovoltaic energy storage in conjunction with distributed energy resources in accordance with the present invention;
FIG. 2 is a block diagram of capacity allocation and operation optimization of the distributed energy system of the integrated energy grid according to the present invention;
FIG. 3 is a block diagram illustrating capacity allocation and operation optimization of the energy storage and conversion device according to the present invention;
FIG. 4 is a peak daily curve of power and heating load in the example;
FIG. 5 is an exploded view of the cost of the energy conversion device planning result in three scenarios of the embodiment;
FIG. 6 is a daily net load curve every half hour for the example;
FIG. 7 is a diagram showing the relationship between the design of the capacity of the battery for storing energy and the price under three scenarios in the embodiment;
fig. 8 is a power balance diagram every half hour for the three scenarios in the example.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments. Note that the following description of the embodiments is merely a substantial example, and the present invention is not intended to be limited to the application or the use thereof, and is not limited to the following embodiments.
Examples
As shown in fig. 1, a capacity allocation and operation optimization method for photovoltaic energy storage and distributed energy resource combination includes the following steps:
(1) acquiring basic data of a distributed energy system and a thermal power and power load prediction data curve accurate to the hour of a year;
(2) according to basic data of the distributed energy system and a thermal power and power load prediction data curve, the capacity of energy conversion equipment in the distributed energy system is configured by taking the minimum annual total investment and operation cost of the distributed energy system as an optimization target;
(3) determining an optimized output curve of the distributed energy system according to the capacity configuration result of the energy conversion equipment;
(4) determining an electric net load curve according to the electric load prediction data curve and the distributed energy system operation optimization output curve;
(5) collecting basic data of a photovoltaic energy storage system;
(6) and optimally configuring the battery energy storage rated power, the battery energy storage rated capacity and the charge and discharge strategy in the photovoltaic energy storage system by taking the minimum annual total investment operation cost of the photovoltaic energy storage system as an optimization target according to the basic data and the electric net load curve of the photovoltaic energy storage system.
Specifically, the basic data of the distributed energy system (DER) in step (1) include investment cost, service life, efficiency, fuel price, carbon strength, carbon price, grid electricity price of CHP power generation, load factor and the like of DER equipment (such as CHP, heat pump and boiler). The time scale of the collected thermal and power load peak daily curve is hours or half-hour or 15 minutes or 5 minutes.
And (2) realizing capacity allocation and operation optimization of the distributed energy system, wherein a block diagram of the capacity allocation and operation optimization of the distributed energy system of the comprehensive energy network is shown in fig. 2, the distribution situation of cold/heat/electricity/gas energy and the load prediction situation in the selected area are comprehensively considered, the minimum total annual investment and operation cost in a planning period is taken as a target, and an optimization framework is shown in fig. 2. The investment cost is determined by the reduction of the total equipment investment to the annual investment cost within its operating period. The model obtains the optimal configuration and operation of the energy conversion equipment with multi-energy coordination by minimizing the annual cost. The input of the optimization problem is the current state parameters of a cold/heat/electricity/gas network, the future new energy, the load output and the like in the planned area, and the output is the results of the capacity, the investment, the operation cost and the like of the energy conversion equipment. The optimization model is used for determining the installation type, the capacity and the operation strategy of the energy conversion equipment, and the multi-energy flow optimization problem is solved by constructing multi-energy flow analysis into equality constraint of the multi-energy flow optimization model and calling an interior point method. The large use of energy conversion equipment such as heat pumps, CHP, air conditioners, electric vehicles and the like affects the power load curve, thereby affecting the configuration and operation of the photovoltaic energy storage system.
Therefore, the step (2) is specifically:
(21) the model obtains the optimal planning configuration of the multifunctional cooperative comprehensive energy system conversion equipment by minimizing the annual total investment and operation cost, and realizes the lowest investment cost CAPEX and operation cost OPEX of the whole energy system. The model optimizes the production and conversion of electric power and heat of each network energy conversion device, so that an optimization objective function of the distributed energy system is established:
Min Ctotal-DER=LFn·CCapex-DER+COpex-DER
wherein, Ctotal-DERFor the annual total investment and operation cost of the distributed energy system, CCapex-DERFor the initial total investment cost, C, of energy conversion equipment in a distributed energy systemOpex-DERAnnual energy conversion equipment operating charge, LF, for a distributed energy systemnIn order to be a coefficient of capital recovery,
Figure GDA0003248035430000081
d is the discount rate, and n is the service life of the energy conversion equipment;
initial total investment cost C of energy conversion equipment in distributed energy systemCapex-DERThe method specifically comprises the following steps:
Figure GDA0003248035430000082
wherein the content of the first and second substances,
Figure GDA0003248035430000083
in order to convert the initial investment cost of the equipment i,
Figure GDA0003248035430000084
the electrical power provided by the inverter device i at peak load,
Figure GDA0003248035430000085
the V-shaped index represents the operator "OR" for the thermal power supplied by the switching device i at load peaks.
Annual energy conversion equipment operating charge C in distributed energy systemOpex-DERThe method specifically comprises the following steps:
Figure GDA0003248035430000086
wherein, CeFor purchase of electricity price, Pimport(t) buying electric power to the upper grid for the t hour of the year, CgAs gas price, vgtotal(t) is the tth hour of the yearTotal gas consumption rate, CcarbonIs the carbon price, ξeIs the carbon emission intensity of the powergIs the carbon emission intensity of natural gas, CO&MFor the operation and maintenance cost, T is the hour of one year, and T is 8760 h.
In the embodiment, the operation cost is optimized based on the load of each hour in the whole year, and then the operation cost in the whole year is obtained through accumulation, and the calculation time of the method is too slow, so that the peak load of the power and the heat in the use year and the introduced load rate are used. In this case, the calculated running cost uses only the peak load, not the 8760 hour load. This approach greatly reduces the complexity and computation time of the problem. If the number of running points is 8760 hours, the running time is shortened by 8760 times. The equation is therefore simplified as follows:
Figure GDA0003248035430000091
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000092
for the operating costs of the system at peak load, etalfThe load factor is used.
(22) Establishing constraint functions of the distributed energy system, including equality constraint functions and inequality constraint functions,
the equality constraint function establishes an energy flow analysis model of independent networks such as a natural gas network and a thermal power network through analogy with power system analysis. The power flow calculation of the power system is used for analyzing the voltage and the power of each node of the power line, so that power supply capacity check, line loss analysis and the like are performed. The heat flow calculation is used for analyzing the temperature, flow, pressure and state change conditions of the heat transfer medium, so that the heat energy loss and the flow pressure loss are analyzed. And (3) calculating and analyzing the flow, pressure and state change conditions of the fuel gas in the transmission and distribution gas pipeline by using the gas flow, thereby analyzing the gas flow pressure loss. Performing combined modeling analysis on each network energy flow and energy conversion equipment of the comprehensive energy system, and considering coupling elements for connecting a power network, a heating power network and a gas network, such as a Combined Heat and Power (CHP), a heat pump, a gas boiler and the like, wherein the network couplingThe resultant part is characterized by a multi-vector conversion efficiency matrix and a permutation matrix. The conversion device model (P, phi, v) in the multi-energy flow joint equation is linked by the conversion efficiency matrixq) And realizing the coupling of the electric/heat/gas network. The multiple energy flow power equality constraint of the comprehensive energy system for balancing the electric power, the heat power, the fuel gas and the load is a combined equation equality of the comprehensive energy system.
Specifically, the equality constraint function is:
Figure GDA0003248035430000093
in the formula, PiIs the active power of the power of node i, ViIs the voltage of node i, VjIs the voltage of node j, NeIs the number of nodes, G, of the power systemijIs the conductance of line ij, BijIs the susceptance, θ, of line ijijIs the voltage angle difference of nodes i and j, QiElectric reactive power, C, for node ipIs the specific heat capacity of water, AhIs a heat power network association matrix and is characterized in that,
Figure GDA0003248035430000101
for mass flow, TsFor temperature of the water supply, ToFor the outlet water temperature, phi is the thermal power vector consumed or provided by each thermal node, BhIs a heat power network loop incidence matrix, KhIs the coefficient of resistance of the heat distribution network pipe, CsAs a water supply network coefficient matrix, bsIs a constant vector, CrAs a coefficient matrix of the backwater network, TrTo the return water temperature, brIs a constant vector, AgIs a gas network incidence matrix, vgIs node air pressure, vqTo node the gas flow rate, BgIs a gas network loop incidence matrix, KgIs the resistance coefficient of the gas network pipeline, and k is an exponential constant;
the inequality constraint function includes:
(a) and (3) limiting the upper limit and the lower limit of the capacity of the energy conversion equipment:
Figure GDA0003248035430000102
Figure GDA0003248035430000103
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000104
for the value of the power active power provided by the conversion device at peak load,
Figure GDA0003248035430000105
the lower limit of the active power of the electricity supplied to the conversion equipment,
Figure GDA0003248035430000106
the upper limit of the electrical active power provided to the conversion equipment,
Figure GDA0003248035430000107
for converting the thermal power supplied by the device at load peaks,
Figure GDA0003248035430000108
the lower limit of the thermal power is,
Figure GDA0003248035430000109
is the upper limit of thermodynamic power;
(b) the power grid inequality constrains:
Vimin≤Vi≤Vimax,i=1,…,Ne
Figure GDA00032480354300001010
Figure GDA00032480354300001011
Figure GDA00032480354300001012
in the formula, ViIs the voltage of node i, ViminLower voltage limit, V, of node iimaxIs the upper voltage limit of node i, NeIs the number of nodes, PgeniFor the active power output of the ith generator,
Figure GDA00032480354300001013
the lower active power output limit of the ith generator,
Figure GDA00032480354300001014
upper limit of active power output, Q, of the ith generatorgeniIs the reactive power output of the ith generator,
Figure GDA00032480354300001015
the lower limit of reactive power output of the ith generator,
Figure GDA00032480354300001016
is the upper limit of reactive power output of the ith generator, NgeNumber of generators, SkIs the electrical power of the kth branch,
Figure GDA00032480354300001017
the lower limit of the electric power of the kth branch,
Figure GDA00032480354300001018
is the upper electric power limit of the kth branch, NleThe number of branches;
(c) the thermal power grid inequality constrains:
Figure GDA00032480354300001019
Ts_min≤Ts≤Ts_max
Tr_min≤Tr≤Tr_max
in the formula (I), the compound is shown in the specification,
Figure GDA00032480354300001020
the flow rate of each pipeline of the heating power network,
Figure GDA00032480354300001021
the lower limit of the flow of each pipeline of the heating power network,
Figure GDA00032480354300001022
is the upper limit of the flow of each pipeline of the heating power network, TsSupply of water temperature, T, to each node of the heating networks_minLower limit of water supply temperature, T, for each node of heating power networks_maxUpper limit of water supply temperature, T, for each node of heating power networkrFor the return water temperature, T, of each node of the heating power networkr_minThe lower limit of the return water temperature of each node of the heating power network, Tr_maxThe upper limit of the return water temperature of each node of the heating power network is set;
(d) the inequality constraint of the gas network is as follows:
pgmin≤pg≤pgmax
vg_min≤vg≤vg_max
in the formula, pgFor the gas pressure, p, at each node of the gas networkgminIs the lower limit of gas pressure of each node of the gas network, pgmaxIs the upper limit of gas pressure of each node of the gas network, vgIs the gas flow of each pipeline of the gas network, vg_minIs the lower limit of the gas flow of each pipeline of the gas network, vg_maxThe upper limit of the gas flow of each pipeline of the gas network;
(23) and optimizing and solving to obtain the capacity of the energy conversion equipment in the distributed energy system.
And (3) determining an optimized output curve of the distributed energy system according to the capacity configuration result of the energy conversion equipment on the basis.
And (4) subtracting the output of the distributed energy from the electric load by the electric net load in the electric net load curve.
And (5) the basic data of the photovoltaic energy storage system comprise a typical daily curve of photovoltaic power generation, PVB investment cost, service life, efficiency, power grid time-sharing electricity price, photovoltaic on-grid electricity price, photovoltaic peak value and annual utilization hours.
After the capacity configuration and operation optimization of the distributed energy system is completed, the capacity configuration and operation optimization of the photovoltaic energy storage system is performed, and fig. 3 is a block diagram of the capacity configuration and operation optimization of the energy storage and conversion equipment.
Specifically, the capacity configuration and operation optimization of the photovoltaic energy storage system are realized in the step (6), which mainly comprises the establishment of an optimization objective function of the photovoltaic energy storage system and the establishment of a constraint function of the photovoltaic energy storage system, and before that, the power flow modeling of the photovoltaic energy storage system is firstly carried out.
Specifically, the power flow of the photovoltaic energy storage system is modeled as:
photovoltaic power PPVIs directly absorbed by load PPV→loadOr storing P by batteryPV→bessOr delivered to higher-level grid PPV→gridSatisfying the equation:
PPV→load(d,t)+PPV→bess(d,t)+PPV→grid(d,t)=PPV(d,t),
PPV→load(d, t) is the power of the photovoltaic supply load at the tth day, PPV→bess(d, t) is the power of the photovoltaic cell to charge the cell at the tth day on day d, PPV→grid(d, t) is the power of photovoltaic externally transmitted to a superior power grid in the tth hour on the day d, PPVAnd (d, t) is the photovoltaic power generation power, d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24.
The power balance equation of the load expresses that the load power is equal to the sum of the photovoltaic power generation, the energy storage charge and discharge and the power of the superior power grid:
Figure GDA0003248035430000111
Figure GDA0003248035430000121
charging power for the t hour on day d,
Figure GDA0003248035430000122
the discharge power at the t hour on day d;
thus, it is possible to obtain:
Figure GDA0003248035430000123
the photovoltaic power flow is dependent on the photovoltaic versus load power as shown in table 1.
(a) If the photovoltaic power is less than the load, the photovoltaic is fully supplied to the load.
(b) If the photovoltaic power is larger than the load and the redundant part does not exceed the energy storage charging limit, the charging power P of the photovoltaic to the batteryPV→bessEqual to this remaining part PPV-PloadElse PPV→bessEqual to the maximum charging power of the battery
Figure GDA0003248035430000124
(c) If the photovoltaic cell has P remaining after charging the battery at maximum powerPV-Pload>PPV→bessAnd then the power P of the photovoltaic outgoing superior power gridPV→gridIf the difference is not between the two, otherwise PPV→grid=0。
TABLE 1 photovoltaic Power output flow Direction vs. load relationship
Figure GDA0003248035430000125
Thus, the step (6) includes the steps of:
(61) establishing an optimization objective function of the photovoltaic energy storage system:
Min Ctotal-PVB=CCapex-PVB+COpex-PVB
wherein, Ctotal-PVBThe total annual investment and operation cost of the photovoltaic energy storage system, CCapex-PVBFor the annual investment costs of the photovoltaic energy storage system, COpex-PVBThe annual running cost of the photovoltaic energy storage system is saved;
(62) establishing a constraint function of the photovoltaic energy storage system, comprising:
(a) charge and discharge power constraint:
Figure GDA0003248035430000126
Figure GDA0003248035430000127
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000128
charging power for the t hour on day d,
Figure GDA0003248035430000129
the discharge power at the t hour on day d,
Figure GDA00032480354300001210
d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24;
(b) energy balance constraint:
Figure GDA00032480354300001211
in the formula, EBESS(d,t)EBESS(t) is the amount of energy stored in the battery at the time of day d and t, EBESS(d,t-1)EBESS(t-1) is the electric quantity of the battery stored energy at the t-1 moment on the day d, etacEfficiency of charging the battery, ηdFor the battery discharge efficiency, Δ t is the time step;
(c) battery storage capacity constraint:
Figure GDA0003248035430000131
in the formula, EBESS(d, t) is the electric quantity stored by the battery at the time t on the day d,
Figure GDA0003248035430000132
storing energy rating, SOC, for a batteryminFor storing the minimum value allowed by the ratio of the residual capacity to the capacity of the battery, SOCmaxThe minimum value allowed by the ratio of the energy storage residual capacity to the capacity of the battery is obtained;
(63) and optimizing and solving to obtain the battery energy storage rated power, the battery energy storage rated capacity and the charging and discharging power at each moment in the photovoltaic energy storage system.
Annual investment cost C of photovoltaic energy storage systemCapex-PVBThe method specifically comprises the following steps:
Figure GDA0003248035430000133
in the formula (I), the compound is shown in the specification,
Figure GDA0003248035430000134
the rated capacity of the energy storage of the battery,
Figure GDA0003248035430000135
the rated power is stored for the battery,
Figure GDA0003248035430000136
in order to obtain the rated output power of the photovoltaic,
Figure GDA0003248035430000137
the capacity investment cost per kWh for battery energy storage,
Figure GDA0003248035430000138
the power investment cost per kW for battery energy storage,
Figure GDA0003248035430000139
is the unit power investment cost of the photovoltaic.
Annual operating charge C of photovoltaic energy storage systemOpex-PVBThe method specifically comprises the following steps:
Figure GDA00032480354300001310
in the formula, PPV→load(d, t) is the power of the photovoltaic supply load at the tth day, PPV→bess(d, t) is the power of the photovoltaic cell to charge the cell at the tth day on day d, PPV→grid(d, t) is the power of photovoltaic externally transmitted to a superior power grid in the tth hour on the day d, PPV(d, t) is the generated power of the photovoltaic,
Figure GDA00032480354300001311
charging power for the t hour on day d,
Figure GDA00032480354300001312
discharge power at day d, t hour, Ce(t) is the power grid electricity price at the moment t,
Figure GDA00032480354300001313
is the electricity price of the photovoltaic power generation outgoing power grid,
Figure GDA00032480354300001314
d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24 for the subsidy price per kW of photovoltaic power generation.
The capacity configuration and operation optimization of the photovoltaic energy storage system are completed, the internal yield rate IRR, the photovoltaic absorption rate SCR and the self-load self-supporting rate SSR of the PVB system can be obtained, and the economic benefit of the photovoltaic system depends on the local absorption of the photovoltaic and is not the income delivered to a superior power grid. The battery energy storage increases matching between the photovoltaic and the load, and improves the self-consistency ratio (SCR) and the self-sufficiency ratio (SSR) of the photovoltaic.
The photovoltaic consumption rate SCR is defined as the amount of electricity E consumed locally by the photovoltaicPV,usedTotal electric quantity E generated by photovoltaicPV,genThe ratio of:
Figure GDA0003248035430000141
the self-sufficiency SSR of the load is defined as the amount E of electricity consumed locally by the photovoltaicPV,usedTotal amount of electricity consumed by load EloadThe ratio of:
Figure GDA0003248035430000142
in the embodiment, campus data of the university of manchester in the united kingdom are adopted, the campus data comprise a 6.6kV power distribution network, a heat power network and a gas network, the gas network is adopted in an east area A of the Oxford road, and the steam heat power network is adopted in a west area B of the Oxford road. The peak daily curves of thermal and electrical load are shown in figure 4. The energy conversion device configuration for the three energy supply scenarios is shown in table 2:
table 2 energy conversion device configuration for three energy supply scenarios
Figure GDA0003248035430000143
The energy price, carbon strength and carbon emission price are shown in table 3:
TABLE 3 energy price, carbon strength and carbon emission price
Figure GDA0003248035430000144
The parameters and prices of the energy conversion devices are shown in table 4:
TABLE 4 parameters and prices of energy conversion devices
Figure GDA0003248035430000145
Factors influencing the economical efficiency of the photovoltaic energy storage system are shown in table 5, the higher the grid peak-valley difference electricity price is, the lower the photovoltaic on-grid electricity price is, the better the energy storage benefit is
TABLE 5 factors affecting photovoltaic energy storage system economics
Figure GDA0003248035430000151
The solution of the optimization problem is accelerated by introducing the load rate into the objective function and parallel computing, and the computing time is less than 1 minute. The planned capacity and cost decomposition of the energy conversion device in the three scenarios is shown in fig. 5, fig. 5(a) is an exploded view of the planned capacity and cost of the energy conversion device in 2016, and fig. 5(b) is a decomposition of the planned capacity and cost of the energy conversion device in 2030. The planning results quantify the impact of energy and carbon prices on energy conversion equipment planning. CHP scenario 2 has advantages over other options without the objective function taking into account carbon prices. Scenario 2 still has advantages if the objective function considers carbon prices using 2016 financial data. However, as the grid carbon strength decreased significantly and the carbon price increased in 2030, the option of heat pump scenario 3 dominates. The results show that when the intensity of the carbon emission of the power grid is reduced to be below 130gCO2/kW, the total cost of the heat pump is lower than that of CHP.
The load peak day operation simulation was performed based on the planned capacities of the boiler, CHP and heat pump, and the net load curve obtained is shown in fig. 6. Where net load refers to the original load plus the power consumed by the heat pump or minus the power generated by the CHP.
The relationship between the optimal configuration capacity of the battery energy storage and the price in the scenario 3 is shown in fig. 7, and fig. 7(a), 7(b), and 7(c) correspond to scenario 1, scenario 2, and scenario 3 in this order.
When the battery energy storage price is reduced to be lower than 150$/kWh, the battery energy storage installation capacity is greatly increased, and the benefit of the photovoltaic energy storage system is greatly improved.
In the CHP scenario, the photovoltaic absorption rate SCR is very low at 0.34, and the self-supporting rate SSR is very high; in the heat pump scenario, the opposite is true, SCR 0.99 ~ 1.
The higher the photovoltaic absorption rate, the higher the Internal Return Rate (IRR) of the photovoltaic energy storage system. In the CHP scenario, most of the photovoltaic is delivered to the upper grid, and therefore the IRR of the photovoltaic energy storage system is low.
The gains of photovoltaic → stored energy (P _ (PV → bess)) come from: 1 storing the difference between the grid electricity price when the redundant photovoltaic power is stored and the photovoltaic outgoing electricity price, and 2 storing the difference between the grid electricity price when the redundant photovoltaic power is stored and the grid electricity price when the part of power is released. In the two-stage time-of-use electricity price, because the electricity transferred to the photovoltaic through the stored energy has no electricity price difference, and the profit of photovoltaic → stored energy comes from the difference between the electricity price of the power grid and the electricity price of photovoltaic delivery, the photovoltaic consumption rate SCR is in direct proportion to the system profit. The use of a heat pump increases the electrical load and therefore increases the photovoltaic rate of absorption SCR and thus the yield of the photovoltaic energy storage system. In the three-part time-of-use electricity price, the photovoltaic power profit absorbed by the heat pump may be smaller than the profit stored in the part of the photovoltaic power and released at the evening peak price, so the use of the heat pump does not necessarily increase the profit of the photovoltaic energy storage system.
Using the predicted 2025 year battery price 96$/kWh, the power balance every half hour for the 3 scenarios is shown in fig. 8, where fig. 8(a), 8(b), and 8(c) correspond in sequence to scenario 1, scenario 2, and scenario 3. The histogram shows the photovoltaic, battery charge and discharge, upper grid, CHP and heat pump power balance. The operation strategy of the battery energy storage is to charge the power grid at the valley price and with redundant photovoltaic power, and discharge the power grid at the peak price and at the evening without photovoltaic power.
In the reference scenario 1, the redundant photovoltaic power is used for storing and charging energy for the battery in the daytime, and the photovoltaic power is released in the evening.
In CHP scenario 2, most of the excess photovoltaic power during the day is sent out to the upper grid, since the power for the load is already supplied by the CHP during the evening.
In heat pump scenario 3, most of the redundant photovoltaic power is consumed by the heat pump during the day, and the battery energy storage is mainly charged by using the midnight valley price.
The energy storage arbitrage of scenario 1 mainly comes from the difference between the grid electricity price and the photovoltaic delivery electricity price. The arbitrage of scenario 3 mainly comes from the peak-to-valley difference of the grid electricity prices.
The above embodiments are merely examples and do not limit the scope of the present invention. These embodiments may be implemented in other various manners, and various omissions, substitutions, and changes may be made without departing from the technical spirit of the present invention.

Claims (5)

1. A capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy sources is characterized by comprising the following steps:
(1) acquiring basic data of a distributed energy system and a thermal power and power load prediction data curve accurate to the hour of a year;
(2) according to basic data of the distributed energy system and a thermal power and power load prediction data curve, the capacity of energy conversion equipment in the distributed energy system is configured by taking the minimum annual total investment and operation cost of the distributed energy system as an optimization target;
(3) determining an optimized output curve of the distributed energy system according to the capacity configuration result of the energy conversion equipment;
(4) determining an electric net load curve according to the electric load prediction data curve and the distributed energy system operation optimization output curve;
(5) collecting basic data of a photovoltaic energy storage system;
(6) according to basic data and an electric net load curve of the photovoltaic energy storage system, optimally configuring battery energy storage rated power, battery energy storage rated capacity and a charging and discharging strategy in the photovoltaic energy storage system by taking the minimum annual total investment operation cost of the photovoltaic energy storage system as an optimization target;
the step (2) is specifically as follows:
(21) establishing an optimization objective function of the distributed energy system:
Min Ctotal-DER=LFn·CCapex-DER+COpex-DER
wherein, Ctotal-DERFor the annual total investment and operation cost of the distributed energy system, CCapex-DERFor the initial total investment cost, C, of energy conversion equipment in a distributed energy systemOpex-DERAnnual energy conversion equipment operating charge, LF, for a distributed energy systemnIn order to be a coefficient of capital recovery,
Figure FDA0003248035420000011
d is the discount rate, and n is the service life of the energy conversion equipment;
(22) establishing constraint functions of the distributed energy system, including equality constraint functions and inequality constraint functions,
wherein the equality constraint function is:
Figure FDA0003248035420000021
in the formula, PiIs the active power of the power of node i, ViIs the voltage of node i, VjIs the voltage of node j, NeIs the number of nodes, G, of the power systemijIs the conductance of line ij, BijIs the susceptance, θ, of line ijijIs the voltage angle difference of nodes i and j, QiElectric reactive power, C, for node ipIs the specific heat capacity of water, AhIs a heat power network association matrix and is characterized in that,
Figure FDA00032480354200000213
for mass flow, TsFor temperature of the water supply, ToFor the outlet water temperature, phi is the thermal power vector consumed or provided by each thermal node, BhIs a heat power network loop incidence matrix, KhIs the coefficient of resistance of the heat distribution network pipe, CsAs a water supply network coefficient matrix, bsIs a constant vector, CrAs a coefficient matrix of the backwater network, TrTo the return water temperature, brIs a constant vector, AgIs a gas network incidence matrix, vgIs node air pressure, vqTo node the gas flow rate, BgIs a gas network loop incidence matrix, KgIs the resistance coefficient of the gas network pipeline, and k is an exponential constant;
the inequality constraint function includes:
(a) and (3) limiting the upper limit and the lower limit of the capacity of the energy conversion equipment:
Figure FDA0003248035420000022
Figure FDA0003248035420000023
in the formula (I), the compound is shown in the specification,
Figure FDA0003248035420000024
for the value of the power active power provided by the conversion device at peak load,
Figure FDA0003248035420000025
the lower limit of the active power of the electricity supplied to the conversion equipment,
Figure FDA0003248035420000026
the upper limit of the electrical active power provided to the conversion equipment,
Figure FDA0003248035420000027
for converting the thermal power supplied by the device at load peaks,
Figure FDA0003248035420000028
the lower limit of the thermal power is,
Figure FDA0003248035420000029
is the upper limit of thermodynamic power;
(b) the power grid inequality constrains:
Vimin≤Vi≤Vimax,i=1,…,Ne
Figure FDA00032480354200000210
Figure FDA00032480354200000211
Figure FDA00032480354200000212
in the formula, ViIs the voltage of node i, ViminIs a nodeLower voltage limit of i, VimaxIs the upper voltage limit of node i, NeIs the number of nodes, PgeniFor the active power output of the ith generator,
Figure FDA0003248035420000031
the lower active power output limit of the ith generator,
Figure FDA0003248035420000032
upper limit of active power output, Q, of the ith generatorgeniIs the reactive power output of the ith generator,
Figure FDA0003248035420000033
the lower limit of reactive power output of the ith generator,
Figure FDA0003248035420000034
is the upper limit of reactive power output of the ith generator, NgeNumber of generators, SkIs the electrical power of the kth branch,
Figure FDA0003248035420000035
the lower limit of the electric power of the kth branch,
Figure FDA0003248035420000036
is the upper electric power limit of the kth branch, NleThe number of branches;
(c) the thermal power grid inequality constrains:
Figure FDA0003248035420000037
Ts_min≤Ts≤Ts_max
Tr_min≤Tr≤Tr_max
in the formula (I), the compound is shown in the specification,
Figure FDA0003248035420000038
the flow rate of each pipeline of the heating power network,
Figure FDA0003248035420000039
the lower limit of the flow of each pipeline of the heating power network,
Figure FDA00032480354200000310
is the upper limit of the flow of each pipeline of the heating power network, TsSupply of water temperature, T, to each node of the heating networks_minLower limit of water supply temperature, T, for each node of heating power networks_maxUpper limit of water supply temperature, T, for each node of heating power networkrFor the return water temperature, T, of each node of the heating power networkr_minThe lower limit of the return water temperature of each node of the heating power network, Tr_maxThe upper limit of the return water temperature of each node of the heating power network is set;
(d) the inequality constraint of the gas network is as follows:
pgmin≤pg≤pgmax
vg_min≤vg≤vg_max
in the formula, pgFor the gas pressure, p, at each node of the gas networkgminIs the lower limit of gas pressure of each node of the gas network, pgmaxIs the upper limit of gas pressure of each node of the gas network, vgIs the gas flow of each pipeline of the gas network, vg_minIs the lower limit of the gas flow of each pipeline of the gas network, vg_maxThe upper limit of the gas flow of each pipeline of the gas network;
(23) optimizing and solving to obtain the capacity of energy conversion equipment in the distributed energy system;
the step (6) is specifically as follows:
(61) establishing an optimization objective function of the photovoltaic energy storage system:
Min Ctotal-PVB=CCapex-PVB+COpex-PVB
wherein, Ctotal-PVBThe total annual investment and operation cost of the photovoltaic energy storage system, CCapex-PVBFor the annual investment costs of the photovoltaic energy storage system, COpex-PVBThe annual running cost of the photovoltaic energy storage system is saved;
(62) establishing a constraint function of the photovoltaic energy storage system, comprising:
(a) charge and discharge power constraint:
Figure FDA00032480354200000311
Figure FDA00032480354200000312
in the formula (I), the compound is shown in the specification,
Figure FDA0003248035420000041
charging power for the t hour on day d,
Figure FDA0003248035420000042
the discharge power at the t hour on day d,
Figure FDA0003248035420000043
d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24;
(b) energy balance constraint:
Figure FDA0003248035420000044
in the formula, EBESS(d,t)EBESS(t) is the amount of energy stored in the battery at the time of day d and t, EBESS(d,t-1)EBESS(t-1) is the electric quantity of the battery stored energy at the t-1 moment on the day d, etacEfficiency of charging the battery, ηdFor the battery discharge efficiency, Δ t is the time step;
(c) battery storage capacity constraint:
Figure FDA0003248035420000045
in the formula, EBESS(d, t) is day d, time tThe amount of electricity that the battery stores,
Figure FDA0003248035420000046
storing energy rating, SOC, for a batteryminFor storing the minimum value allowed by the ratio of the residual capacity to the capacity of the battery, SOCmaxThe minimum value allowed by the ratio of the energy storage residual capacity to the capacity of the battery is obtained;
(63) and optimizing and solving to obtain the battery energy storage rated power, the battery energy storage rated capacity and the charging and discharging power at each moment in the photovoltaic energy storage system.
2. The method of claim 1, wherein the initial total investment cost C of the energy conversion equipment in the distributed energy system isCapex-DERThe method specifically comprises the following steps:
Figure FDA0003248035420000047
wherein the content of the first and second substances,
Figure FDA0003248035420000048
in order to convert the initial investment cost of the equipment i,
Figure FDA0003248035420000049
the electrical power provided by the inverter device i at peak load,
Figure FDA00032480354200000410
the V-shaped index represents the operator "OR" for the thermal power supplied by the switching device i at load peaks.
3. The method of claim 1, wherein the annual energy conversion equipment operating cost C in the distributed energy system isOpex-DERThe method specifically comprises the following steps:
Figure FDA00032480354200000411
wherein, CeFor purchase of electricity price, Pimport(t) buying electric power to the upper grid for the t hour of the year, CgIn order to be the price of the gas,
Figure FDA00032480354200000412
the total gas consumption rate in the t hour of the year, CcarbonIs the carbon price, ξeIs the carbon emission intensity of the powergIs the carbon emission intensity of natural gas, CO&MFor operating maintenance costs, T is the number of hours of a year.
4. The method for optimizing capacity allocation and operation of photovoltaic energy storage and distributed energy resource combination of claim 1, wherein the annual investment cost C of the photovoltaic energy storage systemCapex-PVBThe method specifically comprises the following steps:
Figure FDA0003248035420000051
in the formula (I), the compound is shown in the specification,
Figure FDA0003248035420000052
the rated capacity of the energy storage of the battery,
Figure FDA0003248035420000053
the rated power is stored for the battery,
Figure FDA0003248035420000054
in order to obtain the rated output power of the photovoltaic,
Figure FDA0003248035420000055
the capacity investment cost per kWh for battery energy storage,
Figure FDA0003248035420000056
the power investment cost per kW for battery energy storage,
Figure FDA0003248035420000057
is the unit power investment cost of the photovoltaic.
5. The method of claim 1, wherein the photovoltaic energy storage and distributed energy combined capacity allocation and operation optimization method is characterized in that the photovoltaic energy storage system is operated at an annual cost COpex-PVBThe method specifically comprises the following steps:
Figure FDA0003248035420000058
in the formula, PPV→load(d, t) is the power of the photovoltaic supply load at the tth day, PPV→bess(d, t) is the power of the photovoltaic cell to charge the cell at the tth day on day d, PPV→grid(d, t) is the power of photovoltaic externally transmitted to a superior power grid in the tth hour on the day d, PPV(d, t) is the generated power of the photovoltaic,
Figure FDA0003248035420000059
charging power for the t hour on day d,
Figure FDA00032480354200000510
discharge power at day d, t hour, Ce(t) is the power grid electricity price at the moment t,
Figure FDA00032480354200000511
is the electricity price of the photovoltaic power generation outgoing power grid,
Figure FDA00032480354200000512
d is more than or equal to 1 and less than or equal to 365, and t is more than or equal to 1 and less than or equal to 24 for the subsidy price per kW of photovoltaic power generation.
CN201811368811.1A 2018-11-16 2018-11-16 Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy Active CN109510224B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811368811.1A CN109510224B (en) 2018-11-16 2018-11-16 Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811368811.1A CN109510224B (en) 2018-11-16 2018-11-16 Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy

Publications (2)

Publication Number Publication Date
CN109510224A CN109510224A (en) 2019-03-22
CN109510224B true CN109510224B (en) 2021-11-09

Family

ID=65748884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811368811.1A Active CN109510224B (en) 2018-11-16 2018-11-16 Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy

Country Status (1)

Country Link
CN (1) CN109510224B (en)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110289611A (en) * 2019-03-29 2019-09-27 国网浙江海宁市供电有限公司 A kind of active distribution network planing method promoting distributed generation resource digestion capability
CN110165698B (en) * 2019-04-01 2022-05-24 广西电网有限责任公司 Wind power plant smooth grid-connection method for realizing prospective error asset conversion
CN109919399B (en) * 2019-04-15 2021-07-27 中国科学院电工研究所 Day-ahead economic dispatching method and system for comprehensive energy system
CN110322068A (en) * 2019-07-03 2019-10-11 上海交通大学 A kind of integrated energy system energy transition equipment configuration method and configuration system
CN110610031B (en) * 2019-08-22 2023-04-18 广东电网有限责任公司 Design method of photovoltaic heat pump energy storage comprehensive energy system
CN112819191A (en) * 2019-11-18 2021-05-18 国网能源研究院有限公司 Garden comprehensive energy system operation optimization analysis method oriented to engineering application
CN111244993B (en) * 2020-01-21 2021-11-30 国网湖南省电力有限公司 Capacity optimization configuration method for energy storage participating in power grid peak shaving application
CN111697617B (en) * 2020-03-10 2022-04-05 广东电网有限责任公司广州供电局 High-photovoltaic-permeability power distribution network optimized operation method and system
CN113393010A (en) * 2020-03-12 2021-09-14 中关村海华信息技术前沿研究院 Data processing method, system and storage medium for control system
CN111697604B (en) * 2020-06-19 2021-12-28 中国电力科学研究院有限公司 Configuration method, system and equipment for integrating multiple stations
CN112085377A (en) * 2020-09-04 2020-12-15 浙江大学 Energy internet distributed photovoltaic configuration method under carbon emission
CN112084463B (en) * 2020-09-04 2024-03-15 浙江大学 Distributed photovoltaic configuration method for power system under carbon emission
CN112200347A (en) * 2020-09-08 2021-01-08 国网天津市电力公司电力科学研究院 Method for identifying weak links and expanding and transforming weak links of park comprehensive energy system
CN113051752B (en) * 2021-03-22 2023-04-07 西安热工研究院有限公司 Method for determining optimal heat source of high-pressure air energy storage system electrically coupled with coal
CN114294708B (en) * 2021-12-03 2023-03-14 中国电建集团河北省电力勘测设计研究院有限公司 Method for adjusting heat storage of long-distance heat supply pipe network
CN114301097A (en) * 2022-01-04 2022-04-08 珠海格力电器股份有限公司 Photovoltaic system and electric energy control method thereof
CN114819508B (en) * 2022-03-28 2024-03-29 上海交通大学 Comprehensive energy system distributed photovoltaic maximum admittance capacity calculation method and system
CN115065078B (en) * 2022-05-06 2023-07-04 合肥工业大学 Energy storage capacity configuration method and system in micro-grid environment
CN116151436B (en) * 2023-01-09 2023-08-18 广州汇锦能效科技有限公司 Household-user-oriented photovoltaic building energy planning method and system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103490410A (en) * 2013-08-30 2014-01-01 江苏省电力设计院 Micro-grid planning and capacity allocation method based on multi-objective optimization
CN107104433A (en) * 2017-05-15 2017-08-29 国网江苏省电力公司电力科学研究院 A kind of light-preserved system participates in the acquisition methods of power distribution network Optimal Operation Strategies
CN108695868A (en) * 2018-06-26 2018-10-23 上海交通大学 Power distribution network energy storage addressing constant volume method based on electric power electric transformer

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012129675A1 (en) * 2011-03-31 2012-10-04 Energent Incorporated A computer implemented electrical energy hub management system and method
US20180046160A1 (en) * 2016-08-12 2018-02-15 Robert Bosch Llc Methods and systems for providing photovoltaic plant power feed-in

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103490410A (en) * 2013-08-30 2014-01-01 江苏省电力设计院 Micro-grid planning and capacity allocation method based on multi-objective optimization
CN107104433A (en) * 2017-05-15 2017-08-29 国网江苏省电力公司电力科学研究院 A kind of light-preserved system participates in the acquisition methods of power distribution network Optimal Operation Strategies
CN108695868A (en) * 2018-06-26 2018-10-23 上海交通大学 Power distribution network energy storage addressing constant volume method based on electric power electric transformer

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
考虑分时电价的风光储联合优化调度研究;杨晓萍等;《西安理工大学学报》;20161230(第04期);全文 *
高光伏渗透率配电网中分布式储能系统的优化配置方法;陶琼等;《高电压技术》;20160718(第07期);全文 *

Also Published As

Publication number Publication date
CN109510224A (en) 2019-03-22

Similar Documents

Publication Publication Date Title
CN109510224B (en) Capacity allocation and operation optimization method combining photovoltaic energy storage and distributed energy
Ma et al. The optimal structure planning and energy management strategies of smart multi energy systems
Lin et al. Decentralized solution for combined heat and power dispatch through benders decomposition
CN109919478B (en) Comprehensive energy microgrid planning method considering comprehensive energy supply reliability
CN109659927B (en) Energy storage capacity configuration method of comprehensive energy microgrid considering energy storage participation degree
Mellouk et al. Design and energy management optimization for hybrid renewable energy system-case study: Laayoune region
Torreglosa et al. Control based on techno-economic optimization of renewable hybrid energy system for stand-alone applications
Liu et al. Optimal coordinated operation of a multi-energy community considering interactions between energy storage and conversion devices
CN108154309B (en) Energy internet economic dispatching method considering multi-load dynamic response of cold, heat and electricity
Chen et al. Optimal design and performance assessment for a solar powered electricity, heating and hydrogen integrated energy system
CN111242806A (en) Planning method of electric-thermal-hydrogen multi-energy system considering uncertainty
CN109993345A (en) A kind of isolated operation towards garden is provided multiple forms of energy to complement each other system dynamic economic dispatch method
CN107392791B (en) Distributed photovoltaic and gas-electricity hybrid capacity planning method and system for multi-energy complementary system
Vyas et al. A Review of Algorithms for Control and Optimization for Energy Management of Hybrid Renewable Energy Systems
Xue et al. Optimal capacity allocation method of integrated energy system considering renewable energy uncertainty
Veremiichuk et al. Energy sources selection for industrial enterprise combined power supply system
Liu et al. Influence evaluation of integrated energy system on the unit commitment in power system
CN110224397B (en) User-side battery energy storage cost benefit analysis method under wind and light access background
Meng et al. Economic optimization operation approach of integrated energy system considering wind power consumption and flexible load regulation
Kermani et al. Optimal operation of a real power hub based on PV/FC/GenSet/BESS and demand response under uncertainty
CN114139837A (en) Regional multi-system double-layer distributed optimization scheduling method considering double-layer carbon emission optimization distribution model
Moghaddam et al. A multi-slack optimization model for scheduling energy hubs in smart Grids
CN109345030B (en) Multi-microgrid comprehensive energy system thermoelectric energy flow distribution type optimization method and device
Erdemir et al. Development and comparative evaluation of integrated solar‐driven off‐grid energy systems with hydrogen production and storage options for sustainable buildings
Guo et al. Optimal configuration of multi-energy storage for load aggregators considering user behavior

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 200030 Dongchuan Road, Minhang District, Minhang District, Shanghai

Applicant after: Shanghai Jiaotong University

Address before: 200030 Huashan Road, Shanghai, No. 1954, No.

Applicant before: Shanghai Jiaotong University

GR01 Patent grant
GR01 Patent grant