CN108599137A - A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character - Google Patents

A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character Download PDF

Info

Publication number
CN108599137A
CN108599137A CN201711459115.7A CN201711459115A CN108599137A CN 108599137 A CN108599137 A CN 108599137A CN 201711459115 A CN201711459115 A CN 201711459115A CN 108599137 A CN108599137 A CN 108599137A
Authority
CN
China
Prior art keywords
formula
heat
power
pipeline section
temperature
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201711459115.7A
Other languages
Chinese (zh)
Other versions
CN108599137B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201711459115.7A priority Critical patent/CN108599137B/en
Publication of CN108599137A publication Critical patent/CN108599137A/en
Application granted granted Critical
Publication of CN108599137B publication Critical patent/CN108599137B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28—Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • H02J3/383—
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/50—Photovoltaic [PV] energy
    • Y02E10/56—Power conversion systems, e.g. maximum power point trackers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

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

Abstract

The invention discloses a kind of multipotency streaming system optimizing operation methods considering region heat supply network transient state heat-transfer character, including：Step 1) establishes the region heat supply network transient state heat-transfer character model for considering return pipe net topological structure, pipeline section temperature change dynamic process and pipeline section propagation delay time based on graph theory, which includes mainly six aspect contents：Heat source characteristic, thermic load characteristic, node flow balance, node power fusion characteristics, pipeline section heat-transfer character and pipeline section temporarily/steady state heat transfer Performance criteria；Step 2) establishes multipotency streaming system optimal operation model, with the minimum object function of total operating cost in 4 hours futures of system, the operation constraints that should meet in conjunction with system, determine each equipment in multipotency streaming system it is real-time contribute, the real-time supply and return water temperature etc. of each heat source of heat supply network, realize that the United Dispatching of system and optimization are run.The present invention considers the transient state heat-transfer character of region heat supply network, and the optimized operation state of multipotency streaming system is made more to coincide with actual motion state.

Description

A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character
Technical field
The present invention relates to a kind of multipotency streaming system optimizing operation methods, particularly belong to energy system operation optimisation technique neck Domain.
Background technology
With the increasingly exacerbation of world energy sources crisis and problem of environmental pollution, energy internet, integrated energy system, " mutually The innovation idea such as networking+" wisdom energy continues to bring out, and changes energy supply model, builds the modern energy system of clean and effective It is of great significance.Multipotency streaming system relies on the relevant technologies such as energy source technology, energy conversion and regenerative resource Constantly innovation, the various energy resources networks such as coupling power grid, heat supply network and natural gas grid, is advantageously implemented multipotency cooperative supply and the energy is comprehensive Cascade utilization is closed, one of important trend of energy field is increasingly becoming.
It is one of the important signs that multipotency streaming system is different from traditional energy system, all kinds of energy that various energy resources, which intercouple, Property difference new challenge is proposed to the planning operation of system.Medium of the region heat supply network as transporting heat energy, due to passing There is notable difference with power grid in terms of defeated loss and time scale, attract in the multipotency streaming system of structure multipotency cooperative supply Extensive research concern.
Ripe has been tended to for the research of heat supply network hydraulic regime at present, but the research of thermal drying is still in and explores rank Section.Some scholars are concerned only with the heat-transfer character of heat supply network, using the methods of numerical simulation, Experimental modeling, although can obtain higher The heating power data of precision, but model used is complex, time-consuming for solution, and needs given user's thermic load data and heat Source supply water temperature, it is difficult to be applied in the ultra-short term scheduling of multipotency streaming system.Also some scholars are right towards steady state heat transfer problem A large amount of simplification have been done in the heat-transfer character of heat supply network and network constraint, although so that heat supply network model is easy to solve, work as system call When the time scale of strategy is suitable with heat supply network delay, significant error will be caused.Thus the key of problem is to establish one simultaneous Care for the heat supply network model of accuracy and solving complexity so that the optimized operation state and actual motion state of multipotency streaming system are more While coincideing, energy management Optimized model is easy to solve enough again.
Invention content
The technical problem to be solved by the present invention is to：In view of the deficiencies of the prior art, the invention discloses a kind of consideration areas The multipotency streaming system optimizing operation method of domain heat supply network transient state heat-transfer character, by focusing on opening up for return pipe net in heat supply network Models Sets The delay character that structure, the dynamic process of pipeline section temperature change and pipeline section transmit is flutterred, can reflect system ultra-short term scheduling strategy The dynamic process of lower heat supply network heat transfer so that the optimized operation state of multipotency streaming system is more coincide with actual motion state.Simultaneously It is optimal for target with system economy, in conjunction with the operation constraints that system should meet, determine each equipment in multipotency streaming system It contributes in real time, the real-time supply and return water temperature etc. of each heat source of heat supply network, realizes that the United Dispatching of system and optimization are run.
The present invention uses following technical scheme to solve above-mentioned technical problem：
A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character, includes the following steps：
Step 1)：It is established based on graph theory and considers that return pipe net topological structure, pipeline section temperature change dynamic process and pipeline section pass The region heat supply network transient state heat-transfer character model of defeated time delay, the model include six aspect contents：Heat source characteristic, thermic load characteristic, Node flow balance, node power fusion characteristics, pipeline section heat-transfer character and pipeline section temporarily/steady state heat transfer Performance criteria；
Step 2)：Multipotency streaming system optimal operation model is established, with the operating cost in one rolling scheduling period of system Minimum object function determines going out in real time for each equipment in multipotency streaming system in conjunction with the operation constraints that system should meet The real-time supply and return water temperature of each heat source of power, heat supply network realizes that the United Dispatching of system and optimization are run.
Further, multipotency streaming system optimizing operation method proposed by the invention sees hot water pipe net in step 1) Make fluid network, every pipeline is a branch, and the tie point of heat source, thermic load and pipeline is node；It is temporary to establish region heat supply network The process of state heat-transfer character model is：
Step 101) considers heat source characteristic：
If TSOAnd TSIThe respectively confession of heat source, return water temperature need to ensure the normal heat supply of heat source to heat source Supply and return water temperature is limited, as shown in formula (1)：
In formula, TSO,maxAnd TSO,minThe upper and lower limit of heat source supply water temperature is indicated respectively；TSI,maxAnd TSI,minHeat is indicated respectively The upper and lower limit of source return water temperature；On the other hand, heat source is made of physical equipment, should meet power Climing constant, such as formula (2) It is shown：
In formula, Qs(t) the output thermal power of current scheduling moment heat point source is indicated；Qs(when t+ Δs t) indicates next scheduling The output thermal power of punctum heat source；Δ t indicates scheduling time inter；WithHeat source climbing thermal power is indicated respectively Upper and lower limit；
Step 102) considers thermic load characteristic：
For the heat supply network branch comprising thermic load, the power of taking of confession, the temperature of water return pipeline and thermic load meets formula (3)：
In formula：TLIAnd TLOThe column vector that supply water temperature and return water temperature at respectively all thermic loads are constituted；C is heat Matchmaker's specific heat capacity takes 4.2kJ/ (kg DEG C)；ρ is heating agent density, takes 934.667kg/m3；ηHXFor each heat transmission equipment of heat supply network load side Average efficiency；QL/FFor the column vector for taking power and pipeline flow composition of all thermic loads, as shown in formula (4)：
QL/F=[QL1/fL1,QL2/fL2,…QLk/fLk]TFormula (3)
In formula：QLiAnd fLiIndicate the operating flux of power that i-th of load is taken and pipeline where it respectively, i=1, 2 ..., k, k are load sum；
In order to ensure the heating quality of thermic load, need to limit the supply and return water temperature at thermic load, such as formula (5) institute Show：
In formula, TLIAnd TLOConfession, return water temperature respectively at thermic load；TLI,maxAnd TLI,minIndicate that thermic load supplies respectively The upper and lower limit of coolant-temperature gage；TLO,maxAnd TLO,minThe upper and lower limit of thermic load return water temperature is indicated respectively；
Step 103) considers node flow balance：
For any node in heat supply network, the sum of flow rate of heat medium of inflow is equal to the sum of the flow rate of heat medium of outflow, i.e.,:
In formula,WithThe upper and lower incidence matrix of heat supply network is indicated respectively；F is flow column vector, and upper and lower incidence matrix is determined Shown in justice such as formula (7)：
Shown in the definition such as formula (8) of heat supply network flow column vector：
F=[f1,f2,…,fb]TFormula (7)
In formula：fiIndicate that the flow rate of heat medium of i-th branch, i=1,2 ..., b, b are branch sum, n is node total number；
Step 104) considers node power fusion characteristics：
By law of conservation of energy it is found that for any node in heat supply network, the sum of power of inflow is equal to the power of outflow The sum of, i.e.,：
In formula：TSFAnd TEFThe column vector being made of pipeline section beginning, last temperature and flow, defines such as formula (10)：
In formula：TsiAnd TeiThe initial temperature and terminal temperature of i-th pipeline are indicated respectively；
After power fusion occurs at node, the temperature of the node is equal with the outflow heat medium temperature of the node, i.e.,：
Tni=Ts1,i=Ts2,i=...=Tsk,iFormula (10)
In formula：TniFor the temperature of i-th of node；Ts1,i,Ts2,i,…,Tsk,iIt is connected directly with i-th of node to be all, and And heating agent flows out the initial temperature of the pipeline of the node；
Step 105) considers pipeline section heat-transfer character, when the propagation delay time of pipeline section is much smaller than the time scale of system call strategy When, at each scheduling instance point, heat supply network has reached stable state, and the heat-transfer character of heat supply network is described using steady state heat transfer formula；And When the propagation delay time of pipeline section is suitable with the time scale of system call strategy, the heat transfer of heat supply network is described using transient state heat transfer formula Characteristic, to ensure the accuracy of scheduling strategy；Specifically include step 1051)~1052)：
Step 1051) pipeline section steady state heat transfer characteristic：
For one-dimensional pipeline, known by Fourier law：The thermal power of conduction is directly proportional to temperature gradient and heat transfer area, I.e.：
In formula：QlossThe heat power lost of pipeline is represented, λ indicates the thermal conductivity of pipeline；A indicates the facing area of heat transfer；T Indicate temperature；δ indicates thickness；Negative sign on the right of equal sign indicates that direction of heat flow is always opposite with temperature gradient direction；
Heat power lost of the fluid in pipeline is calculated using multi-layer cylinder wall heat conduction model：
In formula：TaFor ambient temperature, r is the radius of cylinder, riFor the radius of i-th layer of cylinder, λiFor the i-th layer material Thermal conductivity, dl be pipeline infinitesimal length；
Remember that pipeline unit length thermal resistance R is：
Then the thermal power of unit length fluids within pipes loss is：
For a pipeline section, the thermal power of head end is：
Q0=c ρ f (Ts-Tr) formula (15)
In formula：C indicates the specific heat capacity of heating agent；ρ indicates the density of heating agent；F indicates the operating flux of heating agent；TsIndicate heat source Supply water temperature；TrIndicate return water temperature；
Heat supply network transmission thermal power at head end x is equal to the difference of thermal power and heat supply network transmission loss power at heat source, i.e.,：
Fluids within pipes temperature is solved to be distributed as with pipe range：
In view of the coefficient of exponential term x in formula (18) is close to zero, by its Taylor expansion and quadratic term and more high order are omitted , the form after being linearized：
Make formula (19) into matrix form：
TE'=η Ts' formula (19)
In formula：TE' and Ts' it is to be originated by the equivalent pipeline section that temperature and ambient temperature are constituted at pipeline section head end temperature, x respectively Temperature column vector and equivalent pipeline section terminal temperature column vector；η is the matrix being made of pipeline section parameter：
TS'=[Ta-Ts1,Ta-Ts2,…,Ta-Tsb]T
TE'=[Te1-Ts1,Te2-Ts2,…,Teb-Tsb]T
Step 1052) pipeline section transient state heat-transfer character：
It is that the fluid infinitesimal of dl is analyzed to take length in pipeline section, disregards the interaction between fluid infinitesimal, only considers Fluid infinitesimal is by tube wall to extraneous dissipated heat Hloss, from Fourier law：
In formula：T (t) indicate the fluid infinitesimal t moment temperature；
Initial time, the fluid infinitesimal are just sent out from heat source, and the energy having is：
In formula：D is the internal diameter of pipeline section；
The energy that the t moment fluid infinitesimal has subtracts dispersed heat equal to primary power, i.e.,：
The temperature that simultaneous formula (23) and (24) can solve fluid infinitesimal pipeline section Nei changes with time relationship, such as formula (25) shown in：
If the fluid infinitesimal reaches after time t at pipeline section x,：
Formula (26) is substituted into formula (25) temperature of fluid infinitesimal at this time can be obtained and is：
This temperature is the steady temperature at pipeline section x, sees formula (18), i.e., after time t, the temperature at pipeline x is lucky Reach stable state；
For one-dimensional pipeline, heat-transfer character can be by γ1、γ2Two parameter approximations are portrayed；If remembering heating agent infinitesimal from pipeline section The time that head end is flowed into pipeline section end outflow needs is critical duration tc, i.e.,：
In formula：L indicates the length of pipeline section；
Then as time interval Δ t >=γ of system call strategy2tcWhen, in each scheduling instance point tiFlow out pipeline section end Heat medium temperature T (l, ti) all it is up to stable state, therefore T (l, ti) can be acquired by formula (29)：
In formula：Ts(ti-1) indicate to flow into the heat medium temperature at pipeline section beginning in (i-1)-th scheduling instance point；
When the time interval of system call strategy meets γ1tc≤Δt≤γ2tcWhen, in each scheduling instance point tiOutflow Heat medium temperature T (l, the t of pipeline section endi) it will all be in propradation, T (l, ti) uphill process useable linear model approximation carve It draws, therefore can be acquired by formula (30)：
In formula：Tp[Ts(ti-1)] indicate to flow into the steady temperature of pipeline section beginning heating agent in (i-1)-th scheduling instance point；
When the time interval of system call strategy is relatively small, meet 0≤Δ t≤γ1tcWhen, in each scheduling instance point tiFlow out heat medium temperature T (l, the t of pipeline section endi) portrayed with transmission delay model approximation, it can be acquired by formula (31)：
In formula：t0It indicates in initial schedule moment point；Shown in the definition of k such as formula (32)：
K=INT [γ2tc/ Δ t]+1 formula (32)
In formula：INT [] indicates downward bracket function；
Step 106) consider pipeline section temporarily/steady state heat transfer Performance criteria：
By formula (26) it is found that from the heating agent infinitesimal that a pipeline section head end flows into after a dispatching cycle Δ t, if not reaching The end of the pipeline section, then there are transient state heat-transfer characters for the pipeline section, i.e.,：
Formula (33) be heat supply network pipeline section temporarily/steady state heat transfer Performance criteria.
Further, multipotency streaming system optimizing operation method proposed by the invention in step 2), establishes multipotency streaming system The process of optimal operation model is：
Step 201) establishes object function：
For the multipotency streaming system being incorporated into the power networks, the object function of rolling optimization is the fortune made in a rolling scheduling period Row expense is minimum, i.e.,：
Min C=Ce+CgFormula (34)
It is that target determines that next scheduling instance point determines real-time output, the reality of each heat source of heat supply network of each equipment with formula (35) When supply and return water temperature etc.；In formula：C is the operation totle drilling cost in a rolling scheduling period；CeAnd CgRespectively the power purchase of system at Originally and gas cost is bought, sees formula (35)：
In formulaWithElectricity price at respectively t-th of scheduling instance point and Gas Prices,WithRespectively For the combustion gas work(of the power purchase power of power grid, the gas horsepower of gas turbine consumption and gas fired-boiler consumption at t-th of scheduling instance Rate, Δ t are the temporal resolution for rolling the period；
Step 202) establishes constraints, and specific steps include step 2021)~1025)：
Step 2021) establishes power-balance constraint
In formula：PinAnd PoutExpression flows in and out the electrical power of electrical busbar respectively；CinAnd CoutRespectively indicate flow into and Flow out the refrigeration work consumption of air busbar；WithExpression flows in and out the steam thermal power of steam busbar respectively；With Expression flows in and out the hot water thermal power of hot water pipe net respectively；WithThe thermal power and heat that expression heat source is sent out respectively are negative The thermal power that lotus absorbs；
Step 2022) establishes the constraint of heat supply network characteristic
Heat source characteristic meets formula (1) and formula (2), and thermic load characteristic meets formula (3) and formula (5), and node flow balance meets Formula (6), node power fusion characteristics meet formula (9) and formula (11), and pipeline section heat-transfer character meets formula (29), formula (30) and formula (31), pipeline section temporarily/steady state heat transfer Performance criteria meets formula (33)；
Step 2023) establishes non-energy storage device constraint
In formula：Pi、QiAnd CiElectrical power, thermal power and the refrigeration work consumption of a certain equipment output are indicated respectively；Subscript m in tables Show that the lower limit of output power, subscript m ax indicate the upper limit of output power；In addition, should also meet power Climing constant, such as formula (38) It is shown：
In formula：Pi(t)、Qi(t) and Ci(t) electromotive power output, thermal power and the system of current scheduling moment a certain equipment are indicated Cold power；Pi(t+Δt)、Qi(t+ Δ t) and Ci(t+ Δs t) indicates the electromotive power output of next scheduling instance equipment, thermal power And refrigeration work consumption；WithIndicate the electromotive power output of the equipment, the climbing lower limit of thermal power and refrigeration work consumption；WithIndicate the electromotive power output of the equipment, the climbing upper limit of thermal power and refrigeration work consumption；
Step 2024) establishes energy storage device constraint
Storage/exoergic power constraint：
In formula：PES,CAnd PES,DThe charge and discharge power of electric energy storage device is indicated respectively；WithElectric storage device is indicated respectively The lower limit of charge and discharge power；WithThe upper limit of electric storage device charge and discharge power is indicated respectively；QHS,CAnd QHS,DTable respectively Show the charge and discharge power of heat-storing device,WithThe lower limit of the charge and discharge power of heat-storing device is indicated respectively,WithThe upper limit of the charge and discharge power of heat-storing device is indicated respectively；CCS,CAnd CCS,DThe charge and discharge work(that accumulator is set is indicated respectively Rate,The lower limit for the charge and discharge power that accumulator is set is indicated respectively,WithIndicate what accumulator was set respectively The upper limit of charge and discharge power；
Energy storage capacity constrains：
Wimin≤Wi≤WimaxFormula (40)
In formula：WiIndicate the gross energy of the storage of device；WimaxAnd WiminThe upper and lower limit of energy storage capacity is indicated respectively；
Storage/exergonic process constraint：
In formula：W (t) and W (t+1) indicate the energy storage capacity at current time and subsequent time energy storage device respectively, σ be energy oneself Loss rate, PCAnd PDStorage/exoergic power, η are indicated respectivelyCAnd ηDThe efficiency of storage/exoergic is indicated respectively；
Step 2025) establishes state variable constraint
If 0-1 variable Xs (1, t) indicate that the energy accumulating state of t moment energy storage device, X (1, t)=1 indicate that equipment is in accumulation of energy Working condition, X (1, t)=0 indicate that equipment is not in accumulation of energy working condition；Similarly, if 0-1 variable Xs (2, t) indicate t moment The exoergic state of energy storage device cannot be existed simultaneously from the accumulation of energy of energy storage device, exoergic state：
X (1, t)+X (2, t)≤1 formula (42)
If 0-1 variable Xs (3, t) indicate that power purchase state of the t moment multipotency streaming system from power grid, X (3, t)=1 indicate system In the state to power grid power purchase, X (3, t)=0 expression system is not in the state to power grid power purchase；Similarly, if 0-1 variables X (4, t) indicates sale of electricity state of the t moment multipotency streaming system to power grid；It cannot be understood simultaneously to power grid purchase, sale of electricity from system：
X (3, t)+X (4, t)≤1 formula (43)
So far, according to the object function of foundation and the constraints of setting, the reality of each equipment in multipotency streaming system is determined When contribute, the real-time supply and return water temperature of each heat source of heat supply network, realize that the United Dispatching of system and optimization are run.
Further, multipotency streaming system optimizing operation method proposed by the invention, one rolls tune described in step 2) It is 4 hours to spend the period, and the temporal resolution Δ t for rolling the period is 15 minutes.The present invention is using above technical scheme and existing skill Art is compared, and is had the following technical effects：
The multipotency streaming system optimizing operation method proposed by the present invention for considering region heat supply network transient state heat-transfer character, is primarily based on Graph theory establishes region heat supply network transient state heat-transfer character model, and derive judge pipeline section temporarily/foundation of steady state heat transfer characteristic.With it is existing Heat supply network model compare, model emphasis proposed by the invention consider the topological structure of return pipe net, pipeline section temperature change it is dynamic The delay character of state process and pipeline section transmission, uses piecewise-linear techniques in model and describes heat-transfer character, greatly reduce The complexity of model is easy to solve.
Secondly, the present invention is answered with one minimum object function of total operating cost rolled in the period of system in conjunction with system The operation constraints of the satisfaction establishes multipotency streaming system optimal operation model.The Optimized model is MIXED INTEGER linear gauge Model is drawn, ripe solver (such as CPLEX) can be called to be solved, it is possible thereby to determine each equipment in multipotency streaming system It contributes in real time, the real-time supply and return water temperature etc. of each heat source of heat supply network, the United Dispatching and optimization for realizing system are run.The operation side Method considers the transient state heat-transfer character of region heat supply network, and the optimized operation state of multipotency streaming system is made more to be kissed with actual motion state It closes.
In addition, the heat supply network model that the present invention is established is suitable for the two dimension of different topology structure, three-dimensional hot water pipe net, and Can easily the amount of being generalized to adjust, stage by stage change flow matter adjust, intermittent regulation situations such as.
Description of the drawings
Fig. 1 is the flow chart of implementation of the present invention.
Fig. 2 is the system construction drawing in the embodiment of the present invention.
Fig. 3 is the heat supply network structure chart in embodiment in the present invention.
Fig. 4 is the heat supply network schematic diagram in embodiment in the present invention.
Fig. 5 is the heat-net-pipeline heat-transfer character schematic diagram in embodiment in the present invention.
Specific implementation mode
Technical scheme of the present invention is described in further detail below in conjunction with the accompanying drawings：
Those skilled in the art of the present technique are it is understood that unless otherwise defined, all terms used herein (including skill Art term and scientific terminology) there is meaning identical with the general understanding of the those of ordinary skill in fields of the present invention.Also It should be understood that those terms such as defined in the general dictionary should be understood that with in the context of the prior art The consistent meaning of meaning, and unless defined as here, will not be explained with the meaning of idealization or too formal.
System structure in the embodiment of the present invention as shown in Fig. 2, system by gas turbine, gas fired-boiler, waste heat boiler, electricity Refrigeration machine, hot water operated absorption refrigerating machine, heat exchanger and accumulator composition, while photovoltaic apparatus is accessed, the parameter of each equipment is shown in Table 1.In the system, gas turbine, gas fired-boiler, regenerative resource are energy production installations, and utility network is to replenishment system Insufficient electrical energy demands.System is connected with power grid, does not see from power grid power purchase but to power grid sale of electricity, tou power price and Gas Prices Table 3.System provides energy requirement to electricity, heat, cold three kinds of loads simultaneously, and each load is shown in Table 3 with photovoltaic output rolling forecast value.Figure 3 be the heat supply network structure chart in system, and heat supply network schematic diagram is shown in that Fig. 4, heat supply network share 8 water supply lines, 7 water return pipelines and 6 sections Point, parameter are shown in Table 2.Heat source S1 is collectively constituted by gas turbine and waste heat boiler, and heat source S2 is gas fired-boiler；Load QL1, QL2 It is respectively shopping centre, Office Area and residential quarter thermic load with QL3.
Each device parameter in 1 example of table
Heat supply network parameter in 2 example of table
Each load, photovoltaic power generation output forecasting value and tou power price and Gas Prices table in 3 example of table
A kind of multipotency streaming system for consideration region heat supply network transient state heat-transfer character that the embodiment of the present invention proposes optimizes operation side Method, refering to what is shown in Fig. 1, including the following steps：
Step 1)：It is established based on graph theory and considers that return pipe net topological structure, pipeline section temperature change dynamic process and pipeline section pass The region heat supply network transient state heat-transfer character model of defeated time delay, the model include mainly six aspect contents：Heat source characteristic, thermic load are special Property, node flow balance, node power fusion characteristics, pipeline section heat-transfer character and pipeline section temporarily/steady state heat transfer Performance criteria.
Step 2)：Establish multipotency streaming system optimal operation model, in one rolling scheduling period of system it is total operation at This minimum object function determines the real-time of each equipment in multipotency streaming system in conjunction with the operation constraints that system should meet It contributes, the real-time supply and return water temperature etc. of each heat source of heat supply network, realizes that the United Dispatching of system and optimization are run.
In the above-described embodiments, as preference, in the step 1), region heat supply network transient state heat-transfer character model is established Process includes step 101)~106)：
Step 101) considers heat source characteristic：
If TSOAnd TSIThe respectively confession of heat source, return water temperature need to ensure the normal heat supply of heat source to heat source Supply and return water temperature is limited, as shown in formula (1)：
In formula：Subscript indicates the number of heat source.On the other hand, heat source is made of physical equipment, thus is met power and climbed Slope constrains, as shown in formula (2)：
In formula, Qs(t) the output thermal power of current scheduling moment heat point source, unit are indicated：kW；Qs(t+ Δs t) indicates next The output thermal power of a scheduling instance heat point source, unit：kW；Δ t indicates scheduling time inter, takes 15 minutes；Subscript indicates heat The number in source.
Step 102) considers thermic load characteristic：
For the heat supply network branch comprising thermic load, the power of taking of confession, the temperature of water return pipeline and thermic load meets formula (3)：
In formula：TLIAnd TLOSupply water temperature and return water temperature respectively at thermic load；QLPower is taken for thermic load；f The operating flux of pipeline where thermic load；Subscript indicates the number of thermic load.
In order to ensure the heating quality of thermic load, need to limit the supply and return water temperature at thermic load, such as formula (4) institute Show：
In formula, TLIAnd TLOConfession, return water temperature respectively at thermic load；Subscript indicates thermic load number.
Step 103) considers node flow balance：
For any node in heat supply network, the sum of flow rate of heat medium of inflow is equal to the sum of the flow rate of heat medium of outflow, i.e.,：
In formula：It is neutral element that matrix, which does not fill out part,.
Step 104) considers node power fusion characteristics：
By law of conservation of energy it is found that for any node in heat supply network, the sum of power of inflow is equal to the power of outflow The sum of, i.e.,：
In formula：TsiAnd Tei(i=1,2 ..., b) indicates the initial temperature and terminal temperature of i-th pipeline respectively.
After power fusion occurs at node, the temperature of the node is equal with the outflow heat medium temperature of the node, i.e.,：
In formula：TniFor the temperature of i-th of node.
Step 105) considers pipeline section heat-transfer character, and specific steps include step 1051)~1052)：
The heat-transfer character of pipeline section refers to variation characteristic of the heat medium temperature in pipeline section with transmission range and transmission time.Stable state Heat-transfer character refers to that the temperature of heating agent in pipeline section no longer changes over time, and has reached stable state, only consider heat medium temperature with The variation characteristic of transmission range；And the temperature that transient state heat-transfer character describes certain point in pipeline section changes to separately from a stable state The characteristic of one stable state.When the propagation delay time of pipeline section is much smaller than the time scale of system call strategy, in each scheduling instance At point, heat supply network has reached stable state, and steady state heat transfer formula can be used to describe the heat-transfer character of heat supply network；And when the transmission of pipeline section When prolonging suitable with the time scale of system call strategy, transient state heat transfer formula should be used to describe the heat-transfer character of heat supply network, to protect Demonstrate,prove the accuracy of scheduling strategy.
Step 1051) pipeline section steady state heat transfer characteristic：
For one-dimensional pipeline, known by Fourier law：The thermal power of conduction is directly proportional to temperature gradient and heat transfer area. I.e.：
In formula：λ indicates the thermal conductivity of pipeline；A indicates the facing area of heat transfer；T indicates temperature；δ indicates thickness；Equal sign is right The negative sign on side indicates that direction of heat flow is always opposite with temperature gradient direction.
Heat power lost of the fluid in pipeline is calculated using multi-layer cylinder wall heat conduction model：
In formula：TaIt indicates ambient temperature, 0 DEG C is taken in this embodiment.
Remember that pipeline unit length thermal resistance R is：
Then the thermal power of unit length fluids within pipes loss is：
For a pipeline section, the thermal power of head end is：
Q0=c ρ f (Ts-Tr) formula (12)
In formula：C indicates the specific heat capacity of heating agent, takes 4.2kJ/ (kg DEG C)；ρ indicates the density of heating agent, takes 934.667kg/ m3；TsIndicate the supply water temperature of heat source；TrIndicate return water temperature.
Heat supply network transmission thermal power at pipeline section head end x be equal at heat source thermal power and heat supply network transmission loss power it Difference, i.e.,：
Fluids within pipes temperature is solved to be distributed as with pipe range：
In view of the coefficient of exponential term x in formula (14) is close to zero, by its Taylor expansion and quadratic term and more can be omitted High-order term, the form after being linearized：
Make formula (15) into matrix form：
TE'=η Ts' formula (16)
In formula：TE' and Ts' it is to be originated by the equivalent pipeline section that temperature and ambient temperature are constituted at pipeline section head end temperature, x respectively Temperature column vector and equivalent pipeline section terminal temperature column vector；η is the matrix being made of pipeline section parameter：
TS'=[Ta-Ts1,Ta-Ts2,…,Ta-Ts15]T
TE'=[Te1-Ts1,Te2-Ts2,…,Te15-Ts15]T
Step 1052) pipeline section transient state heat-transfer character：
It is that the fluid infinitesimal of dl is analyzed to take length in pipeline section, disregards the interaction between fluid infinitesimal, only considers Fluid infinitesimal is by tube wall to extraneous dissipated heat Hloss, from Fourier law：
In formula：T (t) indicate the fluid infinitesimal t moment temperature.
Initial time, the fluid infinitesimal are just sent out from heat source, and the energy having is：
In formula：D is the internal diameter of pipeline section.
The energy that the t moment fluid infinitesimal has subtracts dispersed heat equal to primary power, i.e.,：
The temperature that simultaneous formula (19) and (20) can solve fluid infinitesimal pipeline section Nei changes with time relationship, such as formula (21) shown in：
If the fluid infinitesimal reaches after time t at pipeline section x,：
Formula (22) is substituted into formula (21) temperature of fluid infinitesimal at this time can be obtained and is：
This temperature is the steady temperature (see formula (14)) at pipeline section x, i.e., after time t, the temperature at pipeline x is lucky Reach stable state.
For one-dimensional pipeline, heat-transfer character can be by γ1、γ2Two parameter approximations are portrayed.Heat-net-pipeline conducts heat in example Characteristic schematic diagram is as shown in figure 5, take：
If the time that note heating agent infinitesimal is flowed into pipeline section end outflow needs from pipeline section head end is critical duration tc, i.e.,：
In formula：L indicates the length of pipeline section.
Then when the time interval of system call strategyWhen, in each scheduling instance point tiFlow out the heat of pipeline section end Matchmaker temperature T (l, ti) all it is up to stable state.Therefore T (l, ti) can be acquired by formula (29)：
In formula：Ts(ti-1) indicate to flow into the heat medium temperature at pipeline section beginning in (i-1)-th scheduling instance point.
When the time interval of system call strategy meetsWhen, in each scheduling instance point tiFlow out pipeline section Heat medium temperature T (l, the t of endi) it will all be in propradation.T(l,ti) uphill process useable linear model approximation portray, because This can be acquired by formula (30)：
In formula：Tp[Ts(ti-1)] indicate to flow into the steady temperature of pipeline section beginning heating agent in (i-1)-th scheduling instance point.
When the time interval of system call strategy is relatively small, meetWhen, in each scheduling instance point tiStream Go out heat medium temperature T (l, the t of pipeline section endi) portrayed with transmission delay model approximation, it can be acquired by formula (28)：
In formula：t0It indicates in initial schedule moment point；Shown in the definition of k such as formula (29)：
In formula：INT [] indicates downward bracket function.
Step 106) consider pipeline section temporarily/steady state heat transfer Performance criteria：
By formula (25) it is found that from the heating agent infinitesimal that a pipeline section head end flows into after a dispatching cycle Δ t, if not reaching The end of the pipeline section, then there are transient state heat-transfer characters for the pipeline section, i.e.,：
Formula (30) be heat supply network pipeline section temporarily/steady state heat transfer Performance criteria.
As preference, in the step 20), the step of establishing multipotency streaming system optimal operation model includes step 201)~202)：
Step 201) establishes object function：
For the multipotency streaming system being incorporated into the power networks, the object function of rolling optimization is the fortune made in a rolling scheduling period Row expense is minimum, i.e.,：
Min C=Ce+CgFormula (31)
It is that target determines that next scheduling instance point determines real-time output, the reality of each heat source of heat supply network of each equipment with formula (31) When supply and return water temperature etc..In formula：C is the operation totle drilling cost in a rolling scheduling period, and it is 4 that the rolling period is taken in the present embodiment Hour, CeAnd CgThe respectively purchases strategies of system and purchase gas cost, are shown in formula (32)：
In formulaWithElectricity price at respectively t-th of scheduling instance point and Gas Prices,WithRespectively For the combustion gas work(of the power purchase power of power grid, the gas horsepower of gas turbine consumption and gas fired-boiler consumption at t-th of scheduling instance Rate, Δ t are temporal resolution, are taken 15 minutes in the present embodiment.
Step 202) establishes constraints, and specific steps include step 2021)~1025)：
Step 2021) establishes power-balance constraint
Step 2022) establishes the constraint of heat supply network characteristic
Heat source characteristic meets formula (1) and formula (2), and thermic load characteristic meets formula (3) and formula (4), and node flow balance meets Formula (5), node power fusion characteristics meet formula (6) and formula (7), and pipeline section heat-transfer character meets formula (26), formula (27) and formula (28), Pipeline section temporarily/steady state heat transfer Performance criteria meets formula (30).In the embodiment of the present invention, 1,3,4,5,6,7,10,11,12 and of number 15 heat-net-pipeline heat-transfer character meets formula (26), and the heat-net-pipeline heat-transfer character of number 2,8,9,13 and 14 meets formula (27).The original state of heat supply network is shown in Table 4.
4 example system initial state of table
Step 2023) establishes non-energy storage device constraint
In formula：The unit of each power is kW.In addition, should also meet power Climing constant, as shown in formula (35)：
In formula：The unit of each power is kW.
Step 2024) establishes energy storage device constraint
Storage/exoergic power constraint：
In formula：PES,CAnd PES,DThe charge and discharge power of accumulator is indicated respectively；The unit of each power is kW.
Energy storage capacity constrains：
200≤WES≤ 900 formulas (37)
In formula：The unit of energy storage capacity is：kJ.
Storage/exergonic process constraint：
In formula：W (t) and W (t+1) indicates the energy storage capacity at current time and subsequent time energy storage device respectively.
Step 2025) establishes state variable constraint
If 0-1 variable Xs (1, t) indicate that the energy accumulating state of t moment energy storage device, X (1, t)=1 indicate that equipment is in accumulation of energy Working condition, X (1, t)=0 indicate that equipment is not in accumulation of energy working condition；Similarly, if 0-1 variable Xs (2, t) indicate t moment The exoergic state of energy storage device.It cannot be existed simultaneously from the accumulation of energy of energy storage device, exoergic state：
X (1, t)+X (2, t)≤1 formula (39)
In the embodiment of the present invention, the object function and constraints of optimal operation model are linear, and are become containing 0-1 Amount, therefore be MIXED INTEGER nonlinear problem, it can be solved, be determined each in multipotency streaming system with the ripe solver such as CPLEX The real-time output of equipment, the real-time supply and return water temperature etc. of each heat source of heat supply network, realize that the United Dispatching of system and optimization are run.
In conclusion the embodiment of the present invention, which is primarily based on graph theory and establishes, considers that return pipe net topological structure, pipeline section temperature become Change the region heat supply network transient state heat-transfer character model of dynamic process and pipeline section propagation delay time, which includes mainly in six aspects Hold：Heat source characteristic, thermic load characteristic, node flow balance, node power fusion characteristics, pipeline section heat-transfer character and pipeline section are temporary/steady State heat-transfer character criterion；Next establishes multipotency streaming system optimal operation model, with total fortune in one rolling scheduling period of system The minimum object function of row cost determines each equipment in multipotency streaming system in conjunction with the operation constraints that system should meet It contributes in real time, the real-time supply and return water temperature of each heat source of heat supply network, realizes that the United Dispatching of system and optimization are run.The operation method is examined The transient state heat-transfer character for having considered region heat supply network makes the optimized operation state of multipotency streaming system more coincide with actual motion state, Heat supply network model used is suitable for the two dimension of different topology structure, three-dimensional hot water pipe net, and can easily the amount of being generalized to adjust, Situations such as changing matter adjusting, the intermittent regulation of flow stage by stage.
The above is only some embodiments of the present invention, it is noted that for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (5)

1. a kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character, which is characterized in that the operation Method includes the following steps：
Step 1)：When establishing consideration return pipe net topological structure, pipeline section temperature change dynamic process and pipeline section transmission based on graph theory The region heat supply network transient state heat-transfer character model prolonged, the model include six aspect contents：Heat source characteristic, thermic load characteristic, node Flow equilibrium, node power fusion characteristics, pipeline section heat-transfer character and pipeline section temporarily/steady state heat transfer Performance criteria；
Step 2)：Multipotency streaming system optimal operation model is established, it is minimum with the operating cost in one rolling scheduling period of system For object function real-time output, the heat of each equipment in multipotency streaming system are determined in conjunction with the operation constraints that system should meet The real-time supply and return water temperature for netting each heat source realizes that the United Dispatching of system and optimization are run.
2. the multipotency streaming system optimizing operation method described in accordance with the claim 1 for considering region heat supply network transient state heat-transfer character, It is characterized in that：In the step 1), regard hot water pipe net as fluid network, every pipeline is a branch, heat source, thermic load Tie point with pipeline is node；The process for establishing region heat supply network transient state heat-transfer character model is：
Step 101) considers heat source characteristic：
If TSOAnd TSIThe respectively confession of heat source, return water temperature are needed to heat source to ensure the normal heat supply of heat source for return water Temperature is limited, as shown in formula (1)：
In formula, TSO,maxAnd TSO,minThe upper and lower limit of heat source supply water temperature is indicated respectively；TSI,maxAnd TSI,minIndicate that heat source returns respectively The upper and lower limit of coolant-temperature gage；On the other hand, heat source is made of physical equipment, should meet power Climing constant, as shown in formula (2)：
In formula, Qs(t) the output thermal power of current scheduling moment heat point source is indicated；Qs(t+ Δs t) indicates next scheduling instance point The output thermal power of heat source；Δ t indicates scheduling time inter；WithThe upper and lower of heat source climbing thermal power is indicated respectively Limit；
Step 102) considers thermic load characteristic：
For the heat supply network branch comprising thermic load, the power of taking of confession, the temperature of water return pipeline and thermic load meets formula (3)：
In formula：TLIAnd TLOThe column vector that supply water temperature and return water temperature at respectively all thermic loads are constituted；C is heating agent specific heat Hold, takes 4.2kJ/ (kg DEG C)；ρ is heating agent density, takes 934.667kg/m3；ηHXFor being averaged for each heat transmission equipment of heat supply network load side Efficiency；QL/FFor the column vector for taking power and pipeline flow composition of all thermic loads, as shown in formula (4)：
QL/F=[QL1/fL1,QL2/fL2,…QLk/fLk]TFormula (3)
In formula：QLiAnd fLiIndicate the operating flux of power that i-th of load is taken and pipeline where it respectively, i=1,2 ..., k, K is load sum；
In order to ensure the heating quality of thermic load, need to limit the supply and return water temperature at thermic load, as shown in formula (5)：
In formula, TLIAnd TLOConfession, return water temperature respectively at thermic load；TLI,maxAnd TLI,minThermic load supply water temperature is indicated respectively Upper and lower limit；TLO,maxAnd TLO,minThe upper and lower limit of thermic load return water temperature is indicated respectively；
Step 103) considers node flow balance：
For any node in heat supply network, the sum of flow rate of heat medium of inflow is equal to the sum of the flow rate of heat medium of outflow, i.e.,:
In formula,WithThe upper and lower incidence matrix of heat supply network is indicated respectively；F is flow column vector, and the definition of upper and lower incidence matrix is such as Shown in formula (7)：
Shown in the definition such as formula (8) of heat supply network flow column vector：
F=[f1,f2,…,fb]TFormula (7)
In formula：fiIndicate that the flow rate of heat medium of i-th branch, i=1,2 ..., b, b are branch sum, n is node total number；
Step 104) considers node power fusion characteristics：
By law of conservation of energy it is found that for any node in heat supply network, the sum of power of inflow is equal to the sum of the power of outflow, I.e.：
In formula：TSFAnd TEFThe column vector being made of pipeline section beginning, last temperature and flow, defines such as formula (10)：
In formula：TsiAnd TeiThe initial temperature and terminal temperature of i-th pipeline are indicated respectively；
After power fusion occurs at node, the temperature of the node is equal with the outflow heat medium temperature of the node, i.e.,：
Tni=Ts1,i=Ts2,i=...=Tsk,iFormula (10)
In formula：TniFor the temperature of i-th of node；Ts1,i,Ts2,i,…,Tsk,iIt is connected directly with i-th of node to be all, and heat Matchmaker flows out the initial temperature of the pipeline of the node；
Step 105) considers pipeline section heat-transfer character, when the propagation delay time of pipeline section is much smaller than the time scale of system call strategy, At each scheduling instance point, heat supply network has reached stable state, and the heat-transfer character of heat supply network is described using steady state heat transfer formula；And when pipe When the propagation delay time of section is suitable with the time scale of system call strategy, the heat transfer that heat supply network is described using transient state heat transfer formula is special Property, to ensure the accuracy of scheduling strategy；Specifically include step 1051)~1052)：
Step 1051) pipeline section steady state heat transfer characteristic：
For one-dimensional pipeline, known by Fourier law：The thermal power of conduction is directly proportional to temperature gradient and heat transfer area, i.e.,：
In formula：QlossThe heat power lost of pipeline is represented, λ indicates the thermal conductivity of pipeline；A indicates the facing area of heat transfer；T is indicated Temperature；δ indicates thickness；Negative sign on the right of equal sign indicates that direction of heat flow is always opposite with temperature gradient direction；
Heat power lost of the fluid in pipeline is calculated using multi-layer cylinder wall heat conduction model：
In formula：TaFor ambient temperature, r is the radius of cylinder, riFor the radius of i-th layer of cylinder, λiFor the heat of the i-th layer material Conductance, dl are the length of pipeline infinitesimal；
Remember that pipeline unit length thermal resistance R is：
Then the thermal power of unit length fluids within pipes loss is：
For a pipeline section, the thermal power of head end is：
Q0=c ρ f (Ts-Tr) formula (15)
In formula：C indicates the specific heat capacity of heating agent；ρ indicates the density of heating agent；F indicates the operating flux of heating agent；TsIndicate the confession of heat source Coolant-temperature gage；TrIndicate return water temperature；
Heat supply network transmission thermal power at head end x is equal to the difference of thermal power and heat supply network transmission loss power at heat source, i.e.,：
Fluids within pipes temperature is solved to be distributed as with pipe range：
In view of the coefficient of exponential term x in formula (18) is close to zero, by its Taylor expansion and quadratic term and more high-order term are omitted, is obtained Form after to linearisation：
Make formula (19) into matrix form：
TE'=η Ts' formula (19)
In formula：TE' and Ts' it is the equivalent pipeline section initial temperature being made of temperature and ambient temperature at pipeline section head end temperature, x respectively Column vector and equivalent pipeline section terminal temperature column vector；η is the matrix being made of pipeline section parameter：
TS'=[Ta-Ts1,Ta-Ts2,…,Ta-Tsb]T
TE'=[Te1-Ts1,Te2-Ts2,…,Teb-Tsb]T
Step 1052) pipeline section transient state heat-transfer character：
It is that the fluid infinitesimal of dl is analyzed to take length in pipeline section, disregards the interaction between fluid infinitesimal, only considers fluid Infinitesimal is by tube wall to extraneous dissipated heat Hloss, from Fourier law：
In formula：T (t) indicate the fluid infinitesimal t moment temperature；
Initial time, the fluid infinitesimal are just sent out from heat source, and the energy having is：
In formula：D is the internal diameter of pipeline section；
The energy that the t moment fluid infinitesimal has subtracts dispersed heat equal to primary power, i.e.,：
The temperature that simultaneous formula (23) and (24) can solve fluid infinitesimal pipeline section Nei changes with time relationship, such as formula (25) institute Show：
If the fluid infinitesimal reaches after time t at pipeline section x,：
Formula (26) is substituted into formula (25) temperature of fluid infinitesimal at this time can be obtained and is：
This temperature is the steady temperature at pipeline section x, sees formula (18), i.e., after time t, the temperature at pipeline x reaches just Stable state；
For one-dimensional pipeline, heat-transfer character can be by γ1、γ2Two parameter approximations are portrayed；If remembering heating agent infinitesimal from pipeline section head end It is critical duration t to be flowed into the time that the outflow of pipeline section end needsc, i.e.,：
In formula：L indicates the length of pipeline section；
Then as time interval Δ t >=γ of system call strategy2tcWhen, in each scheduling instance point tiFlow out the heat of pipeline section end Matchmaker temperature T (l, ti) all it is up to stable state, therefore T (l, ti) can be acquired by formula (29)：
In formula：Ts(ti-1) indicate to flow into the heat medium temperature at pipeline section beginning in (i-1)-th scheduling instance point；
When the time interval of system call strategy meets γ1tc≤Δt≤γ2tcWhen, in each scheduling instance point tiFlow out pipeline section Heat medium temperature T (l, the t of endi) it will all be in propradation, T (l, ti) uphill process useable linear model approximation portray, because This can be acquired by formula (30)：
In formula：Tp[Ts(ti-1)] indicate to flow into the steady temperature of pipeline section beginning heating agent in (i-1)-th scheduling instance point；
When the time interval of system call strategy is relatively small, meet 0≤Δ t≤γ1tcWhen, in each scheduling instance point tiStream Go out heat medium temperature T (l, the t of pipeline section endi) portrayed with transmission delay model approximation, it can be acquired by formula (31)：
In formula：t0It indicates in initial schedule moment point；Shown in the definition of k such as formula (32)：
K=INT [γ2tc/ Δ t]+1 formula (32)
In formula：INT [] indicates downward bracket function；
Step 106) consider pipeline section temporarily/steady state heat transfer Performance criteria：
By formula (26) it is found that from the heating agent infinitesimal that a pipeline section head end flows into after a dispatching cycle Δ t, if not reaching the pipe The end of section, then there are transient state heat-transfer characters for the pipeline section, i.e.,：
Formula (33) be heat supply network pipeline section temporarily/steady state heat transfer Performance criteria.
3. the multipotency streaming system optimizing operation method of region heat supply network transient state heat-transfer character is considered according to claim 2, It is characterized in that：In the step 2), the process for establishing multipotency streaming system optimal operation model is：
Step 201) establishes object function：
For the multipotency streaming system being incorporated into the power networks, the object function of rolling optimization is the running cost made in a rolling scheduling period With minimum, i.e.,：
Min C=Ce+CgFormula (34)
It is that target determines that next scheduling instance point determines the real-time confession of the real-time output, each heat source of heat supply network of each equipment with formula (35) Return water temperature etc.；In formula：C is the operation totle drilling cost in a rolling scheduling period；CeAnd CgRespectively the purchases strategies of system and Gas cost is bought, sees formula (35)：
In formulaWithElectricity price at respectively t-th of scheduling instance point and Gas Prices,WithRespectively t The gas horsepower of the power purchase power of power grid, the gas horsepower of gas turbine consumption and gas fired-boiler consumption, Δ t at a scheduling instance To roll the temporal resolution in period；
Step 202) establishes constraints, and specific steps include step 2021)~1025)：
Step 2021) establishes power-balance constraint
In formula：PinAnd PoutExpression flows in and out the electrical power of electrical busbar respectively；CinAnd CoutIt indicates to flow in and out respectively The refrigeration work consumption of air busbar；WithExpression flows in and out the steam thermal power of steam busbar respectively；WithRespectively Expression flows in and out the hot water thermal power of hot water pipe net；WithThe thermal power and thermic load suction that heat source is sent out are indicated respectively The thermal power of receipts；
Step 2022) establishes the constraint of heat supply network characteristic
Heat source characteristic meets formula (1) and formula (2), and thermic load characteristic meets formula (3) and formula (5), and node flow balance meets formula (6), node power fusion characteristics meet formula (9) and formula (11), and pipeline section heat-transfer character meets formula (29), formula (30) and formula (31), Pipeline section temporarily/steady state heat transfer Performance criteria meets formula (33)；
Step 2023) establishes non-energy storage device constraint
In formula：Pi、QiAnd CiElectrical power, thermal power and the refrigeration work consumption of a certain equipment output are indicated respectively；Subscript m in indicates defeated Go out the lower limit of power, subscript m ax indicates the upper limit of output power；In addition, should also meet power Climing constant, such as formula (38) institute Show：
In formula：Pi(t)、Qi(t) and Ci(t) electromotive power output, thermal power and the refrigeration work(of current scheduling moment a certain equipment are indicated Rate；Pi(t+Δt)、Qi(t+ Δ t) and Ci(t+ Δs t) indicates electromotive power output, thermal power and the system of next scheduling instance equipment Cold power；WithIndicate the electromotive power output of the equipment, the climbing lower limit of thermal power and refrigeration work consumption；WithIndicate the electromotive power output of the equipment, the climbing upper limit of thermal power and refrigeration work consumption；
Step 2024) establishes energy storage device constraint
Storage/exoergic power constraint：
In formula：PES,CAnd PES,DThe charge and discharge power of electric energy storage device is indicated respectively；WithIndicate that electric storage device is filled, put respectively The lower limit of electrical power；WithThe upper limit of electric storage device charge and discharge power is indicated respectively；QHS,CAnd QHS,DHeat accumulation is indicated respectively The charge and discharge power of device,WithThe lower limit of the charge and discharge power of heat-storing device is indicated respectively,WithRespectively Indicate the upper limit of the charge and discharge power of heat-storing device；CCS,CAnd CCS,DThe charge and discharge power that accumulator is set is indicated respectively,The lower limit for the charge and discharge power that accumulator is set is indicated respectively,WithWhat expression accumulator was set respectively fills, puts The upper limit of electrical power；
Energy storage capacity constrains：
Wimin≤Wi≤WimaxFormula (40)
In formula：WiIndicate the gross energy of the storage of device；WimaxAnd WiminThe upper and lower limit of energy storage capacity is indicated respectively；
Storage/exergonic process constraint：
In formula：W (t) and W (t+1) indicates that the energy storage capacity at current time and subsequent time energy storage device, σ are that energy damages certainly respectively Rate, PCAnd PDStorage/exoergic power, η are indicated respectivelyCAnd ηDThe efficiency of storage/exoergic is indicated respectively；
Step 2025) establishes state variable constraint
If 0-1 variable Xs (1, t) indicate that the energy accumulating state of t moment energy storage device, X (1, t)=1 indicate that equipment is in accumulation of energy work State, X (1, t)=0 indicate that equipment is not in accumulation of energy working condition；Similarly, if 0-1 variable Xs (2, t) indicate t moment energy storage The exoergic state of equipment cannot be existed simultaneously from the accumulation of energy of energy storage device, exoergic state：
X (1, t)+X (2, t)≤1 formula (42)
If 0-1 variable Xs (3, t) indicate that power purchase state of the t moment multipotency streaming system from power grid, X (3, t)=1 expression system are in To the state of power grid power purchase, X (3, t)=0 expression system is not in the state to power grid power purchase；Similarly, if 0-1 variable Xs (4, T) sale of electricity state of the t moment multipotency streaming system to power grid is indicated；It cannot be understood simultaneously to power grid purchase, sale of electricity from system：
X (3, t)+X (4, t)≤1 formula (43)
So far, according to the object function of foundation and the constraints of setting, going out in real time for each equipment in multipotency streaming system is determined The real-time supply and return water temperature of each heat source of power, heat supply network realizes that the United Dispatching of system and optimization are run.
4. the multipotency streaming system optimizing operation method described in accordance with the claim 1 for considering region heat supply network transient state heat-transfer character, It is characterized in that：A rolling scheduling period described in step 2) is 4 hours.
5. the multipotency streaming system optimizing operation method described in accordance with the claim 3 for considering region heat supply network transient state heat-transfer character, It is characterized in that：The temporal resolution Δ t for rolling the period is 15 minutes.
CN201711459115.7A 2017-12-28 2017-12-28 Multi-energy flow system optimization operation method considering transient heat transfer characteristics of regional heat supply network Active CN108599137B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711459115.7A CN108599137B (en) 2017-12-28 2017-12-28 Multi-energy flow system optimization operation method considering transient heat transfer characteristics of regional heat supply network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711459115.7A CN108599137B (en) 2017-12-28 2017-12-28 Multi-energy flow system optimization operation method considering transient heat transfer characteristics of regional heat supply network

Publications (2)

Publication Number Publication Date
CN108599137A true CN108599137A (en) 2018-09-28
CN108599137B CN108599137B (en) 2020-03-31

Family

ID=63633236

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711459115.7A Active CN108599137B (en) 2017-12-28 2017-12-28 Multi-energy flow system optimization operation method considering transient heat transfer characteristics of regional heat supply network

Country Status (1)

Country Link
CN (1) CN108599137B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241676A (en) * 2018-10-18 2019-01-18 东南大学 The selection method of heat supply network and building temporal resolution in integrated energy system
CN109523137A (en) * 2018-10-29 2019-03-26 天津大学 Consider the garden comprehensive energy Optimization Scheduling of building thermal load demands response
CN109697308A (en) * 2018-11-30 2019-04-30 天津大学 A kind of natural gas transmission systematic steady state modeling method considering pipe network time-delay characteristics
CN109726906A (en) * 2018-12-21 2019-05-07 国网浙江省电力有限公司电力科学研究院 Co-generation unit dispatching method a few days ago based on the constraint of heat supply network partial differential equation
CN109740242A (en) * 2018-12-29 2019-05-10 重庆大学 Consider that the electric-gas integrated energy system of natural gas thermal procession unifies energy flux computation method
CN109840305A (en) * 2019-03-26 2019-06-04 中冶赛迪技术研究中心有限公司 A kind of steam pipe network waterpower-thermal calculation method and system
CN110762602A (en) * 2019-10-08 2020-02-07 华南理工大学 Energy storage quantitative calculation method for distributed energy heat supply pipe network
CN111209703A (en) * 2020-01-03 2020-05-29 云南电网有限责任公司电力科学研究院 Delay-considered regional steam heat supply network topological structure optimization method and system
CN112362096A (en) * 2020-10-26 2021-02-12 南方电网科学研究院有限责任公司 Method and device for monitoring running state of multi-energy flow, terminal equipment and storage medium
CN113434988A (en) * 2021-06-11 2021-09-24 西安交通大学 Heat transfer characteristic analysis method for heat pipe in nuclear power system
CN113537264A (en) * 2020-04-21 2021-10-22 阿里巴巴集团控股有限公司 Space application state detection method, management method, device and equipment
CN117332725A (en) * 2023-12-01 2024-01-02 中国科学院广州能源研究所 Steam network dynamic calculation method and device of comprehensive energy system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040202457A1 (en) * 2003-04-08 2004-10-14 Lg Electronics Inc. Water heating apparatus and refrigerator having the same
CN103593540A (en) * 2013-11-28 2014-02-19 中南大学 Method for determining position of root of blast furnace cohesive zone through multi-source information fusion
CN106447113A (en) * 2016-10-08 2017-02-22 东南大学 Multi-region integrated energy system operation method based on operation optimization model
CN106998079A (en) * 2017-04-28 2017-08-01 东南大学 A kind of modeling method of combined heat and power Optimal Operation Model

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040202457A1 (en) * 2003-04-08 2004-10-14 Lg Electronics Inc. Water heating apparatus and refrigerator having the same
CN103593540A (en) * 2013-11-28 2014-02-19 中南大学 Method for determining position of root of blast furnace cohesive zone through multi-source information fusion
CN106447113A (en) * 2016-10-08 2017-02-22 东南大学 Multi-region integrated energy system operation method based on operation optimization model
CN106998079A (en) * 2017-04-28 2017-08-01 东南大学 A kind of modeling method of combined heat and power Optimal Operation Model

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YANNAN SUN ET AL.: "The Impact of Uncertain Physical Parameters on HVAC Demand Response", 《IEEE TRANSACTIONS ON SMART GRID》 *
顾伟等: "多区域综合能源系统热网建模及系统运行优化", 《中国电机工程学报》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241676A (en) * 2018-10-18 2019-01-18 东南大学 The selection method of heat supply network and building temporal resolution in integrated energy system
CN109523137A (en) * 2018-10-29 2019-03-26 天津大学 Consider the garden comprehensive energy Optimization Scheduling of building thermal load demands response
CN109523137B (en) * 2018-10-29 2022-11-04 天津大学 Garden comprehensive energy optimization scheduling method considering building thermal load demand response
CN109697308A (en) * 2018-11-30 2019-04-30 天津大学 A kind of natural gas transmission systematic steady state modeling method considering pipe network time-delay characteristics
CN109697308B (en) * 2018-11-30 2023-04-07 天津大学 Natural gas transmission system steady-state modeling method considering delay characteristic of pipe network
CN109726906A (en) * 2018-12-21 2019-05-07 国网浙江省电力有限公司电力科学研究院 Co-generation unit dispatching method a few days ago based on the constraint of heat supply network partial differential equation
CN109740242B (en) * 2018-12-29 2020-11-24 重庆大学 Unified energy flow calculation method of electricity-gas comprehensive energy system considering natural gas thermal process
CN109740242A (en) * 2018-12-29 2019-05-10 重庆大学 Consider that the electric-gas integrated energy system of natural gas thermal procession unifies energy flux computation method
CN109840305A (en) * 2019-03-26 2019-06-04 中冶赛迪技术研究中心有限公司 A kind of steam pipe network waterpower-thermal calculation method and system
CN109840305B (en) * 2019-03-26 2023-07-18 中冶赛迪技术研究中心有限公司 Steam pipe network hydraulic-thermal calculation method and system
CN110762602A (en) * 2019-10-08 2020-02-07 华南理工大学 Energy storage quantitative calculation method for distributed energy heat supply pipe network
CN111209703A (en) * 2020-01-03 2020-05-29 云南电网有限责任公司电力科学研究院 Delay-considered regional steam heat supply network topological structure optimization method and system
CN111209703B (en) * 2020-01-03 2023-10-13 云南电网有限责任公司电力科学研究院 Regional steam heating network topology optimization method and system considering delay
CN113537264A (en) * 2020-04-21 2021-10-22 阿里巴巴集团控股有限公司 Space application state detection method, management method, device and equipment
CN113537264B (en) * 2020-04-21 2022-09-20 阿里巴巴集团控股有限公司 Space application state detection method, management method, device and equipment
CN112362096A (en) * 2020-10-26 2021-02-12 南方电网科学研究院有限责任公司 Method and device for monitoring running state of multi-energy flow, terminal equipment and storage medium
CN112362096B (en) * 2020-10-26 2022-04-12 南方电网科学研究院有限责任公司 Method and device for monitoring running state of multi-energy flow, terminal equipment and storage medium
CN113434988A (en) * 2021-06-11 2021-09-24 西安交通大学 Heat transfer characteristic analysis method for heat pipe in nuclear power system
CN113434988B (en) * 2021-06-11 2022-10-28 西安交通大学 Heat transfer characteristic analysis method for heat pipe in nuclear power system
CN117332725A (en) * 2023-12-01 2024-01-02 中国科学院广州能源研究所 Steam network dynamic calculation method and device of comprehensive energy system

Also Published As

Publication number Publication date
CN108599137B (en) 2020-03-31

Similar Documents

Publication Publication Date Title
CN108599137A (en) A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character
CN106998079B (en) Modeling method of combined heat and power optimization scheduling model
CN106447113B (en) A kind of multizone integrated energy system operation method based on optimal operation model
CN107067116A (en) A kind of multizone electric heating integrated system economic environment combined dispatching method for solving
Liu et al. Coordinated operation of multi-integrated energy system based on linear weighted sum and grasshopper optimization algorithm
CN107808218A (en) Urban energy internet tidal current computing method based on hotspot stress regulation
CN111313429B (en) Reliability assessment method and system for comprehensive energy system
CN106056478B (en) The Interval Power Flow calculation method of heat supply network in a kind of electric-thermal coupled system
CN113190942B (en) Method and device for calculating virtual energy storage capacity of heat supply/cold system and electronic equipment
CN110110913A (en) Large-scale garden integrated energy system energy source station Optimal Configuration Method
CN106339794A (en) Electric-thermal coupling multi-energy flow network node energy price calculation method
CN109389248B (en) Thermoelectric power coordinated scheduling method of comprehensive energy system based on regional heat supply network
CN106844895A (en) A kind of decoupling calculation method of supply of cooling, heating and electrical powers microgrid energy stream
WO2019075872A1 (en) Intra-day rolling scheduling method for thermoelectric coupling multi-energy-flow system
Blacha et al. Dynamic simulation of bidirectional low-temperature networks-a case study to facilitate the integration of renewable energies
CN110245819A (en) A kind of energy efficiency evaluating method of multiple-energy-source interacted system
CN113091124B (en) Regional energy supply capacity-based electric-heat complementary heating system partition regulation and control method
Cui et al. Optimal operation of CCHP microgrids with multiple shiftable loads in different auxiliary heating source systems
Liu et al. A new power flow model for combined heat and electricity analysis in an integrated energy system
CN111783309A (en) Dynamic simulation method of steam heating network based on internal conservation
Dai et al. Modeling and analysis of electrical heating system based on entransy dissipation-based thermal resistance theory
CN115234973A (en) Method for determining dynamic heat supply strategy of heat supply network
CN110119524A (en) A kind of energy resource system Multiple Time Scales interaction optimizing operation method and system
CN115062555A (en) Integrated energy system exergy flow direct calculation method based on unbalanced node exergy
Wang et al. Integrated energy system optimal operation using Data-Driven district heating network model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant