CN109412164A - A kind of electric-thermal association system trend processing method - Google Patents

A kind of electric-thermal association system trend processing method Download PDF

Info

Publication number
CN109412164A
CN109412164A CN201811639858.7A CN201811639858A CN109412164A CN 109412164 A CN109412164 A CN 109412164A CN 201811639858 A CN201811639858 A CN 201811639858A CN 109412164 A CN109412164 A CN 109412164A
Authority
CN
China
Prior art keywords
trend
node
load
thermal medium
electric
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
CN201811639858.7A
Other languages
Chinese (zh)
Other versions
CN109412164B (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.)
Zhejiang University ZJU
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Zhejiang University ZJU
State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, State Grid Jiangsu Electric Power Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201811639858.7A priority Critical patent/CN109412164B/en
Publication of CN109412164A publication Critical patent/CN109412164A/en
Application granted granted Critical
Publication of CN109412164B publication Critical patent/CN109412164B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • 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/46Controlling of the sharing of output between the generators, converters, or transformers

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of electric-thermal association system trend processing methods.For electric-thermal association system, the known basic network data of following power grid and heat supply network are obtained by the sensor in power grid or by tool acquisition before calculating, voltage V after the trend of each node in the power grid after trend iteration is obtained according to known basic network data processingi' and heat supply network in each node trend after thermal medium flow into temperature, and continuous trend handles iteration, after each trend iteration, calculate the maximum value that trend iteration front and back heat supply network interior joint thermal medium flows into the maximum value of temperature difference and the difference of each node voltage of power grid: and judge that realization completes trend and handles until reaching requirement.Electric-thermal association system trend processing of the invention can carry out calculation processing more quickly, committed memory is small, and power flow solutions accuracy is good, can optimize operation for the building of electric-thermal association system, accurate data is provided to dispatching of power netwoks, improves the applicability of direct trend method.

Description

A kind of electric-thermal association system trend processing method
Technical field
The invention belongs to the Operation of Electric Systems and control technology field of the form containing various energy resources, have been specifically related to one kind Electric-thermal association system trend processing method,
Background technique
Currently, with the large-scale application of cogeneration of heat and power technology, electric-thermal association system effectively improves energy benefit as one kind It is rapidly developed with the form of efficiency.The system that is coupled to of the two different form of energy determines that operating parameter and state become again There are many new challenges when amount.
Existing trend processing method is only capable of the trend of processing system for distribution network of power, can not handle two kinds of energy couplings of electric-thermal The case where conjunction, thus it is existing lacked it is a kind of can be for electric-thermal association system trend processing method.
Summary of the invention
For the problems in above-mentioned background technology, the invention proposes a kind of electric-thermal association system trend processing method, energy Suitable for the processing of electric-thermal association system trend, solve the processing of electric-thermal association system trend in the prior art inaccuracy, when Between the long technical problem accounted for more than memory.
The technical scheme is that
Electric-thermal association system of the present invention includes three parts of power grid, heat supply network and cogeneration units, thermoelectricity connection It produces unit and is separately connected heat supply network and power grid, cogeneration units provide thermal energy and electric energy, foundation to heat supply network and power grid simultaneously respectively The power parameter and power output mode of cogeneration units proportionally carry out the output distribution of electric energy and thermal energy;Heat supply network is divided into heat supply Net and backheat net, heating network and backheat net are both connected between cogeneration units and heating equipment, the heat of cogeneration units Medium is delivered to heating equipment through heating network as thermal energy, and the thermal medium of heating equipment is delivered to thermoelectricity through heating network as thermal energy Coproduction unit;There are load bus in power grid and heat supply network, the route or pipeline being connected between each node are branch, negative in power grid Lotus node is electric load node, and the load bus in heat supply network is heating power load bus, is saved with higher level's supply transformer as source Point, higher level's supply transformer and cogeneration units convey to power grid and provide electric energy;
The method of the present invention including the following steps:
1) it is obtained known to following power grid and heat supply network by the sensor in power grid before calculating or by tool acquisition Basic network data, comprising:
Electric-thermal association system: total node number n, route and duct length Lij, wherein i, j indicate the ordinal number of node, i, j ∈n;
Power grid: route unit impedance Z, voltage rating VN, electric load node power consumption power PLoad.i
Heat supply network: the rated temperature T of the heat transfer coefficient λ of pipeline unit length, heating networks.N, backheat net rated temperature To.N、 The thermal energy consumption power φ of heating power load busLoad.i, heating network rated temperature Ts.N, backheat net rated temperature To.N;Simultaneously The thermal medium of heating power load bus flows out temperature To.Load.iThe as rated temperature T of backheat neto.N, i.e. To.Load.i=To.N;Thermoelectricity The thermal medium of coproduction unit flows out temperature To.CHPThe as rated temperature T of heating networks.N, i.e. Ts.N=To.CHP
Cogeneration units: power output the power proportions k, φ of electric energy and thermal energyCHP=k × PCHP, wherein PCHPFor cogeneration of heat and power The electric energy of unit goes out activity of force, φCHPGo out activity of force for the thermal energy of cogeneration units;
2) according to the following formula, temperature T is flowed into using thermal medium before the trend of heating power load buss.Load.i, heating power load The thermal medium of node flows out temperature To.Load.i, heating power load bus thermal energy consumption power φLoad.iWith the specific heat capacity of thermal medium Cp, calculate thermal medium flow m before the trend of each heating power load busq.Load.i
3) then by thermal medium flow m before the trend of each load bus of above-mentioned steps calculating gainedq.Load.iComposition is load The thermal medium traffic matrix m of nodeq.Load:
It is calculated by the following formula the thermic load traffic matrix m for obtaining heat supply network branch again:
M=Amq.Load (2)
Wherein, A is node branch incidence matrix, and m indicates the thermic load traffic matrix of heat supply network branch, mijTo be flowed through in heat supply network The thermic load flow of branch between node i and node j;
4) by following formula, the temperature T out of backheat net thermal medium from each branch is obtainedend.H.ij:
Wherein, Tstart.H.ijIndicate the inflow temperature of branch thermal medium in backheat net, λ indicates the heat transfer of pipeline unit length Coefficient, mH.ijIndicate the flow of branch thermal medium between backheat net interior joint i and node j, LijIndicate heat supply network interior joint i and node The length of bypass line between j, e indicate natural Exponents e, are the truth of a matter of natural logrithm;
5) by following formula, the thermal medium for obtaining the node to cross in backheat net there are pipeline flows out temperature Tout.H:
Wherein, Tout.HWith Tin.HRespectively indicate the thermal medium outflow temperature and stream of the node to cross in backheat net there are pipeline Enter temperature;mout.HWith min.HIt respectively indicates the thermal medium outflow flow of the node to cross in backheat net there are pipeline and flows into stream Amount;
6) temperature T is flowed out with the thermal medium of the node to cross in backheat net there are pipelineout.HAs cogeneration units Thermal medium flows into temperature Ts.CHP, i.e. Ts.CHP=Tout.H;Flow is flowed out with the thermal medium of the node to cross in backheat net there are pipeline mout.HThermal medium as cogeneration units flows through flow mq.CHP, i.e. mq.CHP=mout.H;Then pass through following formula, according to The thermal medium of cogeneration units flows into temperature Ts.CHP, cogeneration units thermal medium flow through flow mq.CHPAnd cogeneration of heat and power The thermal medium of unit flows out temperature To.CHP, thermal energy goes out activity of force φ before obtaining the trend of cogeneration unitsCHP:
φCHP=Cp·mq.CHP·(Ts.CHP-To.CHP) (5)
7) by following formula, according to the power output power proportions k and cogeneration of heat and power of the electric energy of cogeneration units and thermal energy Thermal energy goes out activity of force φ before the trend of unitCHP, obtain electric energy power output power P before the trend of cogeneration unitsCHP:
8) by following formula, connection cogeneration units (CHP) is obtained respectively and without connection cogeneration units (CHP) in the case of, the outflow electric current I of each electric load node in power gridi:
If each electric load node does not connect cogeneration units (CHP) in power grid, according to power load each in power grid The power consumption power P of lotus nodeLoad.iWith voltage V before the trend of each electric load nodeiEach electric load node is calculated Outflow electric current Ii:
Wherein, ()*The conjugation of representing matrix calculates;
If each electric load node is connected with cogeneration units (CHP) in power grid, according to electric load each in power grid The power consumption power P of nodeLoad.i, by before the trend of cogeneration units electric energy contribute power PCHPWith each electric load node Trend before voltage ViThe outflow electric current I of each electric load node is calculatedi:
9) the outflow electric current I of each electric load node obtained by step 8)iForm the outflow electric current of electric load node Matrix I:
By following formula, electricity is obtained according to the outflow current matrix I of electric load node and node branch incidence matrix A The current matrix B of net branch:
B=AT·I (9)
Wherein, BijThe electric current of branch road between power grid interior joint i and node j;
10) it is acquired by the outflow current matrix I of electric load node and node branch incidence matrix A with direct trend method The tidal current voltage V of each electric load node in power gridi′;
11) by following formula, for an electric load node for being connected with cogeneration units, according to the power load Voltage V after the trend of lotus nodei', the outflow electric current I of the electric load nodeiWith the power consumption power of the electric load node PLoad.i, obtain electric energy power output power P after the trend of cogeneration units 'CHP:
P′CHP=PLoad.i-Ii·Vi′ (10)
Wherein, PLoad.iIndicate the power consumption power of each electric load node in power grid;
12) by following formula, according to the power output power proportions k and cogeneration of heat and power of the electric energy of cogeneration units and thermal energy After the trend of unit electric energy power output power P 'CHP, thermal energy goes out activity of force φ ' after obtaining the trend of cogeneration unitsCHP:
φ′CHP=k × P 'CHP (11)
13) by following formula, go out activity of force φ ' using thermal energy after the trend of cogeneration unitsCHP, cogeneration of heat and power machine The thermal medium of group flows into temperature Ts.CHPTemperature T is flowed out with the thermal medium of cogeneration unitso.CHP, obtain cogeneration units Thermal medium flows through flow m ' after trendq.CHP:
14) flow m ' is flowed through with thermal medium after the trend of cogeneration unitsq.CHPThermal medium flowing out stream as heating network Measure mG, temperature T is flowed out with the thermal medium of cogeneration unitso.CHPThermal medium as heating network flows out temperature TG;Then pass through Following formula calculates the thermal medium outflow temperature T for obtaining the node to cross in heating network there are pipelineout.G:
Wherein, Tout.GWith Tin.GRespectively indicate the thermal medium outflow temperature and stream of the node to cross in heating network there are pipeline Enter temperature;mout.GWith min.GIt respectively indicates the thermal medium outflow flow of the node to cross in heating network there are pipeline and flows into stream Amount;
Temperature T is flowed out with the thermal medium of the node i to cross in heating network there are pipelineout.GAfter the trend of the node i Thermal medium flows into temperature T 's.Load.i, flow T is flowed out with the thermal medium of the node i to cross in heating network there are pipelineout.GAs Thermal medium flow m ' after the trend that the heating equipment of the node i flows throughq.Load.i
15) by following formula, the inflow temperature T of the thermal medium from each branch in heating network is calculatedstart.G.ijWith outflow Temperature Tend.G.ij:
Wherein, Tstart.G.ijIndicate the inflow temperature of each branch thermal medium in heating network, λ indicates the biography of pipeline unit length Hot coefficient, mH.ijIndicate the flow of branch thermal medium between backheat net interior joint i and node j, LijIndicate heat supply network interior joint i and section The length of bypass line between point j, e indicate natural Exponents e, are the truth of a matter of natural logrithm;
16) by following formula, the maximum value of heat supply network interior joint thermal medium inflow temperature difference before and after trend iteration is calculated ΔTs.maxWith the maximum value Δ V of the difference of each node voltage of power gridi.max:
ΔTs.max=max (| T 's.Load.i-Ts.Load.i|) (15)
ΔVi.max=max (| Vi′-Vi|) (16)
Wherein, T 's.Load.iThermal medium flows into temperature after indicating the trend of the node i to cross in heating network there are pipeline, Ts.Load.iThermal medium flows into temperature before indicating the trend of the node i to cross in heating network there are pipeline;ViIndicate electric load section Voltage before the trend of point i, voltage after the trend of V ' expression electric load node i;
17) step (2)~(16) are constantly iteratively repeated and carry out trend iteration, after each trend iteration, in the following ways It carries out judging whether trend iteration restrains;
If restraining after this trend iteration, the processing result after exporting this trend iteration obtains this trend iteration Voltage V after the trend of each node in power grid afterwardsi' and heat supply network in each node trend after thermal medium flow into temperature T 's.Load.i
If not converged after this trend iteration, using the processing result after this trend iteration as next trend iteration before Numerical value before trend, by voltage V after the trend of each node in the power grid after this trend iterationi' as next trend iteration when Voltage V before the trend of each node in power gridi, thermal medium after the trend of each node in the heat supply network after this trend iteration is flowed into temperature Spend T 's.Load.iThermal medium flows into temperature T before the trend of each node in heat supply network when as next trend iterations.Load.i, carry out down Trend iterative processing.
In the step 17), whether trend iteration convergent to be sentenced and holds back foundation are as follows: heat supply network node thermal medium before and after trend iteration Flow into the maximum value Δ T of temperature differences.maxWith the maximum value Δ V of grid nodes difference in voltagei.maxWhether it is respectively less than and is equal to 10-5:
ΔTs.max,ΔVi.max≤10-5 (17)
If meeting above-mentioned formula, restrained after this trend iteration;Otherwise it is not yet restrained after this trend iteration.
In the step 17), whether trend iteration convergent to be sentenced and holds back foundation are as follows: it is more than 100 that trend, which iterates to calculate number k, It is secondary, then judge that the electric-thermal association system trend does not restrain, terminates electric-thermal association system Load flow calculation.
By voltage V before the trend of electric load node each in power gridiIt is initially set to voltage rating VN, by heating power each in heat supply network Thermal medium flows into temperature T before the trend of load buss.Load.iIt is initially set to the rated temperature T of heating networks.N
Beneficial effects of the present invention:
Electric-thermal association system trend processing of the invention can carry out calculation processing more quickly, and committed memory is small, trend knot Fruit accuracy is good, can optimize operation for the building of electric-thermal association system, provides accurate data to dispatching of power netwoks, and it is accurate to provide Network state.
The present invention makes full use of coupling (heating power node with power node have identical relevance) of the heat supply network with power grid, heat Net and power grid this feature of node branch incidence matrix A having the same, greatly reduce the EMS memory occupation of program.
The improved electric-thermal association system trend processing of the present invention improves the applicability of direct trend method, and is being promoted The high efficiency and high robust of the direct trend of original power system are still maintained on the basis of applicability.
Detailed description of the invention
Fig. 1 shows simple electric-thermal association system structural schematic diagrams of the invention.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples.
The electric-thermal association system that the present invention is embodied is specific as follows:
Heat supply network transmits energy by achieving the purpose that using heat supply pipeline conveying thermal medium (steam or hot water).Work as heat supply network When using cogeneration units etc. as heat supply network heat source, therrmodynamic system, electric system thus coupling are formd increasingly complex Electric-thermal association system.
Electric-thermal association system includes three parts of power grid, heat supply network and cogeneration units, and cogeneration units connect respectively Heat supply network and power grid are connect, cogeneration units provide thermal energy and electric energy to heat supply network and power grid simultaneously respectively, according to cogeneration units Power parameter and power output mode proportionally carry out the output distribution of electric energy and thermal energy;Heat supply network is divided into heating network and backheat net, Heating network and backheat net are both connected between cogeneration units and heating equipment, and the thermal medium of cogeneration units is as thermal energy It is delivered to heating equipment through heating network, the thermal medium of heating equipment is delivered to cogeneration units through heating network as thermal energy;I.e. Heating network is the piping network that thermal medium is delivered to heating equipment by cogeneration units, and backheat net is that heat is situated between by heating equipment Matter is delivered to the piping network of cogeneration units;Also, heating network and backheat net have same paths, but pipeline is different.
Heat supply network trend process is as follows: cogeneration units generate thermal energy, are sent to the thermal medium of thermal energy by heating network Remaining thermal medium is back to cogeneration units by backheat net after heating equipment uses thermal energy by each heating equipment, is completed Heat supply network heat supply.
There are load bus in power grid and heat supply network, the route or pipeline being connected between each node are branch, and route is electric power Route, pipeline are heat distribution pipeline, and the load bus in power grid is electric load node, and the load bus in heat supply network is heating power load Node becomes source node with higher level's supply transformer, and higher level's supply transformer and cogeneration units, which are conveyed to power grid, to be provided Electric energy;The equipment for collecting, distributing and transmitting energy (heat energy/electric energy) is known as node, wherein be known as the equipment for using energy Higher level's supply transformer is become source node by load bus, and the route or pipeline being connected between each node are branch.
It is of the invention based on direct trend method in power grid and heat supply network topological structure having the same, born in power grid and heat supply network Connection structure of the lotus node in two networks is identical, utilizes the matrix and heating power system constructed in trend method direct in electric system Discharge model and temperature model in system, the algorithm of simultaneous solution electric-thermal association system trend.The trend of electric-thermal association system Calculate the voltage V for calculating each load bus in power gridi, in heat supply network each load bus thermal medium flow mqWith inflow temperature Ts.Load.i, cogeneration units electric energy PCHPWith thermal energy φCHPActivity of force out.
The embodiment of the present invention and its implementation process are as follows:
1) it is obtained known to following power grid and heat supply network by the sensor in power grid before calculating or by tool acquisition Basic network data, comprising:
Electric-thermal association system: total node number n, route and duct length Lij, wherein i, j indicate the ordinal number of node, i, j ∈n;
Power grid: route unit impedance Z, voltage rating VN, electric load node power consumption power PLoad.i
Heat supply network: the rated temperature T of the heat transfer coefficient λ of pipeline unit length, heating networks.N, backheat net rated temperature To.N、 The thermal energy consumption power φ of heating power load busLoad.i, heating network rated temperature Ts.N, backheat net rated temperature To.N;Simultaneously The thermal medium of heating power load bus flows out temperature To.Load.iThe as rated temperature T of backheat neto.N, i.e. To.Load.i=To.N;Thermoelectricity The thermal medium of coproduction unit flows out temperature To.CHPThe as rated temperature T of heating networks.N, i.e. Ts.N=To.CHP
Cogeneration units: power output the power proportions k, φ of electric energy and thermal energyCHP=k × PCHP, wherein PCHPFor cogeneration of heat and power The electric energy of unit goes out activity of force, φCHPGo out activity of force for the thermal energy of cogeneration units;
For example following electric-thermal association system is embodied:
Electric-thermal association system as shown in Figure 1 is made of four nodes.Wherein, there are four nodes, on node 1 represents Grade supply transformer is source node;Node 2 and node 3 are load bus, and node 4 is cogeneration units.Node 1, node It is connected between 2 and node 3 by grid branch (i.e. power circuit), node 2 and node 3 pass through heat supply network branch (i.e. thermal pipe Road) it is connected to cogeneration units (CHP).
And obtain the known basic network data of power grid and heat supply network are as follows:
Power grid:
Power node number: 3;
Line parameter circuit value: 3.7 Ω of route resistance per unit length/km, 0.369 Ω of route unit length reactance/km
Line length: route 12:0.01km;Route 23:0.01km
Power node voltage rating: 220V;
The power of electric load node electrical equipment: power node 2: active power 500kW;Power node 3: active power 400kW。
Heat supply network:
Heating power number of nodes: 3;
Pipe parameter: the heat transfer coefficient of pipeline unit length: 0.321Wm-1·℃-1
Duct length: pipeline 23:0.015km;Pipeline 34:0.0175km;
Heating power heating network rated temperature: 70 DEG C;
Heating power backheat net rated temperature: 30 DEG C;
The power of heating power load bus heating equipment: heating power node 2:0.0805kW: heating power node 3:0.107kW.
Cogeneration units:
Electric-thermal goes out force mode: electricity production power is 1.3 times of heat production power.
2) during electric-thermal association system trend, for the first time in iterative processing, by the trend of electric load node each in power grid Preceding voltage ViIt is initially set to voltage rating VN, thermal medium before the trend of heating power load bus each in heat supply network is flowed into temperature Ts.Load.i It is initially set to the rated temperature T of heating networks.N
According to the following formula, temperature T is flowed into using thermal medium before the trend of heating power load buss.Load.i, heating power load section The thermal medium of point flows out temperature To.Load.i, heating power load bus thermal energy consumption power φLoad.iWith the specific heat capacity C of thermal mediump, Calculate thermal medium flow m before the trend of each heating power load busq.Load.i
3) then by thermal medium flow m before the trend of each load bus of above-mentioned steps calculating gainedq.Load.iComposition is load The thermal medium traffic matrix m of nodeq.Load:
It is calculated by the following formula the thermic load traffic matrix m for obtaining heat supply network branch again:
M=Amq.Load (2)
Wherein, A is node branch incidence matrix, represents the relevance between node and branch, and m indicates heat supply network branch Thermic load traffic matrix, mijFor the thermic load flow for flowing through the branch between node i and node j in heat supply network;
In example electric-thermal association system as shown in Figure 1, the thermic load traffic matrix of heat supply network branch isNode Branch incidence matrix isThe thermal medium traffic matrix of load bus isTherefore haveI.e.
4) by following formula, the temperature T out of backheat net thermal medium from each branch is obtainedend.H.ij:
Wherein, Tstart.H.ijIndicate the inflow temperature of branch thermal medium in backheat net, λ indicates the heat transfer of pipeline unit length Coefficient, mH.ijIndicate the flow of branch thermal medium between backheat net interior joint i and node j, LijIndicate heat supply network interior joint i and node The length of bypass line between j, e indicate natural Exponents e, are the truth of a matter of natural logrithm;
In example electric-thermal association system as shown in Figure 1, in backheat net, between load bus 2 and load bus 3 Branch 23, the inflow temperature T of thermal mediumstart.H.23Temperature T is flowed out for the thermal medium of load bus 2o.Load.2, i.e., Tstart.H.23=To.Load.2.Its interior thermal medium flow mH.23It is acquired by formula (2), i.e. mH.23=mq.Load.2.And it is asked according to formula (3) 23 end of branch and the heat medium temperature T before crossing with node 3 in backheat net outend.H.23
5) by following formula, the thermal medium for obtaining the node to cross in backheat net there are pipeline flows out temperature Tout.H:
Wherein, Tout.HWith Tin.HRespectively indicate the thermal medium outflow temperature and stream of the node to cross in backheat net there are pipeline Enter temperature;mout.HWith min.HIt respectively indicates the thermal medium outflow flow of the node to cross in backheat net there are pipeline and flows into stream Amount;
The thermal medium of above-mentioned cogeneration units flows into temperature Ts.CHPFlow is flowed through with the thermal medium of cogeneration units mq.CHPPartial parameters as cogeneration units.
In example electric-thermal association system as shown in Figure 1, there are pipelines to cross for node 3, and heat is situated between in the backheat net of the node There are two the channels that mass flow enters: the thermal medium of thermal medium and 3 heating equipment of node that branch 23 flows out flowed through.
Therefore: two thermal mediums inflow flows are respectively as follows: m in the backheat net of node 3in.H.1=m23、min.H.2= mq.Load.3, corresponding two thermal mediums flow into temperature and are respectively as follows: Tin.H.1=Tend.23、Tin.H.2=To.Load.3.Meanwhile it saving Thermal medium only has an efflux channel in the backheat net of point 3, flows out flow mout.HTo flow into the sum of flow, it may be assumed that mout.H=∑ min.H.Accordingly, the thermal medium that can find out node 3 according to each known parameters and formula (4) flows out temperature Tout.H
To sum up, thermal medium flows out temperature T in the backheat net of node 3out.HAs the thermal medium of backheat net flows out temperature TH, section The thermal medium of point 3 flows out flow mout.HAs the thermal medium of backheat net flows out flow mH
According to step 3)~step 5), the calculating of all trend parameters in backheat net is completed.
6) temperature T is flowed out with the thermal medium of the node to cross in backheat net there are pipelineout.HAs cogeneration units Thermal medium flows into temperature Ts.CHP, i.e. Ts.CHP=Tout.H;Flow is flowed out with the thermal medium of the node to cross in backheat net there are pipeline mout.HThermal medium as cogeneration units flows through flow mq.CHP, i.e. mq.CHP=mout.H;Then pass through following formula, according to The thermal medium of cogeneration units flows into temperature Ts.CHP, cogeneration units thermal medium flow through flow mq.CHPAnd cogeneration of heat and power The thermal medium of unit flows out temperature To.CHP, thermal energy goes out activity of force φ before obtaining the trend of cogeneration unitsCHP:
φCHP=Cp·mq.CHP·(Ts.CHP-To.CHP) (5)
7) by following formula, according to the power output power proportions k and cogeneration of heat and power of the electric energy of cogeneration units and thermal energy Thermal energy goes out activity of force φ before the trend of unitCHP, obtain electric energy power output power P before the trend of cogeneration unitsCHP:
8) by following formula, connection cogeneration units (CHP) is obtained respectively and without connection cogeneration units (CHP) in the case of, the outflow electric current I of each electric load node in power gridi:
If each electric load node does not connect cogeneration units (CHP) in power grid, according to power load each in power grid The power consumption power P of lotus nodeLoad.iWith voltage V before the trend of each electric load nodeiEach electric load node is calculated Outflow electric current Ii:
Wherein, ()*The conjugation of representing matrix calculates;
If each electric load node is connected with cogeneration units (CHP) in power grid, according to electric load each in power grid The power consumption power P of nodeLoad.i, by before the trend of cogeneration units electric energy contribute power PCHPWith each electric load node Trend before voltage ViThe outflow electric current I of each electric load node is calculatedi:
9) the outflow electric current I of each electric load node obtained by step 8)iForm the outflow electric current of electric load node Matrix I:
By following formula, electricity is obtained according to the outflow current matrix I of electric load node and node branch incidence matrix A The current matrix B of net branch:
B=AT·I (9)
Wherein, BijThe electric current of branch road between power grid interior joint i and node j;
In example electric-thermal association system as shown in Figure 1, the current matrix of grid branch isLoad bus Flowing out current matrix isIt is formedI.e.
10) it is acquired by the outflow current matrix I of electric load node and node branch incidence matrix A with direct trend method The tidal current voltage V of each electric load node in power gridi′;
Direct trend method uses J.H.Teng, " A direct approach for distribution system Load flow solutions, " in IEEE Trans.Power Delivery, vol.18, pp.882-887, July 2003. Direct trend method.The outflow current matrix I of electric load node and node branch incidence matrix A are input to direct trend method In can obtain voltage V after trendi′。
11) by following formula, for an electric load node for being connected with cogeneration units, according to the power load Voltage V after the trend of lotus nodei', the outflow electric current I of the electric load nodeiWith the power consumption power of the electric load node PLoad.i, obtain electric energy power output power P after the trend of cogeneration units 'CHP:
P′CHP=PLoad.i-Ii·Vi′ (10)
Wherein, PLoad.iIndicate the power consumption power of each electric load node in power grid;
12) by following formula, according to the power output power proportions k and cogeneration of heat and power of the electric energy of cogeneration units and thermal energy After the trend of unit electric energy power output power P 'CHP, thermal energy goes out activity of force φ ' after obtaining the trend of cogeneration unitsCHP:
φ′CHP=k × P 'CHP (11)
13) by following formula, go out activity of force φ ' using thermal energy after the trend of cogeneration unitsCHP, cogeneration of heat and power machine The thermal medium of group flows into temperature Ts.CHPTemperature T is flowed out with the thermal medium of cogeneration unitso.CHP, obtain cogeneration units Thermal medium flows through flow m ' after trendq.CHP:
14) flow m ' is flowed through with thermal medium after the trend of cogeneration unitsq.CHPThermal medium flowing out stream as heating network Measure mG, temperature T is flowed out with the thermal medium of cogeneration unitso.CHPThermal medium as heating network flows out temperature TG;Then pass through Following formula calculates the thermal medium outflow temperature T for obtaining the node to cross in heating network there are pipelineout.G:
Wherein, Tout.GWith Tin.GRespectively indicate the thermal medium outflow temperature and stream of the node to cross in heating network there are pipeline Enter temperature;mout.GWith min.GIt respectively indicates the thermal medium outflow flow of the node to cross in heating network there are pipeline and flows into stream Amount;
Temperature T is flowed out with the thermal medium of the node i to cross in heating network there are pipelineout.GAfter the trend of the node i Thermal medium flows into temperature T 's.Load.i, flow T is flowed out with the thermal medium of the node i to cross in heating network there are pipelineout.GAs Thermal medium flow m ' after the trend that the heating equipment of the node i flows throughq.Load.i
In example electric-thermal association system as shown in Figure 1, there are pipelines to cross for node 3, and heat is situated between in the heating network of the node The channel that mass flow enters is the thermal medium of cogeneration units, i.e. min.G=m 'q.CHP, corresponding thermal medium flows into temperature and is Tin.G=To.CHP.Meanwhile there are two the channels that thermal medium flows out in the heating network of node 3: the thermal medium and section that branch 23 flows into The thermal medium of 3 heating equipments of point flowed through, and the sum of the two outflow flow is inflow flow min.G.Therefore: the heating network of node 3 In two thermal mediums flow into flows and be respectively as follows: mout.G.1=m23、mout.G.2=min.G-mout.G.1.Accordingly, according to each known parameters Thermal medium in the heating network of node 3, which is found out, with formula (13) flows out temperature Tout.G
To sum up, thermal medium flow m ' after the trend of calculating 3 heating equipment of posterior nodal point flowed throughq.Load.3The as confession of node 2 Thermal medium flows to the outflow flow m of 3 heating equipment of node in heat supply networkout.G.2, it may be assumed that m 'q.Load.3=mout.G.2;Calculate posterior nodal point 3 Trend after thermal medium flow into temperature T 's.Load.3Thermal medium flows out temperature T as in the heating network of node 3out.G, i.e. T 's.Load.3 =Tout.G
15) by following formula, the inflow temperature T of the thermal medium from each branch in heating network is calculatedstart.G.ijWith outflow Temperature Tend.G.ij:
Wherein, Tstart.G.ijIndicate the inflow temperature of each branch thermal medium in heating network, λ indicates the biography of pipeline unit length Hot coefficient, mH.ijIndicate the flow of branch thermal medium between backheat net interior joint i and node j, LijIndicate heat supply network interior joint i and section The length of bypass line between point j, e indicate natural Exponents e, are the truth of a matter of natural logrithm;
In example electric-thermal association system as shown in Figure 1, for heating network branch 23, the inflow temperature of thermal medium Tstart.G.ijTemperature T is flowed out for thermal medium in the heating network of node 3out.G, i.e. Tstart.G.ij=Tout.G.Its interior thermal medium flow mG.23The outflow flow m of branch 23 is flowed to for thermal medium in the heating network of node 3out.G.1, i.e. mG.23=mout.G.1.The heat of node i Medium flows into temperature T 's.Load.iFor the heat medium temperature T at the end heating network branch iend.G.ij, and heat supply is found out according to formula (14) The heat medium temperature T at 23 end of net branchend.G.23, calculate the thermal medium inflow temperature T ' of posterior nodal point 2s.Load.2As heating network branch The heat medium temperature T at 23 end of roadend.G.23, i.e. T 's.Load.2=Tend.G.23
According to step (13)~step (15), the calculating of all trend parameters in heating network is completed.
16) by following formula, the maximum value of heat supply network interior joint thermal medium inflow temperature difference before and after trend iteration is calculated ΔTs.maxWith the maximum value Δ V of the difference of each node voltage of power gridi.max:
ΔTs.max=max (| T 's.Load.i-Ts.Load.i|) (15)
ΔVi.max=max (| Vi′-Vi|) (16)
Wherein, T 's.Load.iThermal medium flows into temperature after indicating the trend of the node i to cross in heating network there are pipeline, Ts.Load.iThermal medium flows into temperature before indicating the trend of the node i to cross in heating network there are pipeline;ViIndicate electric load section Voltage before the trend of point i, voltage after the trend of V ' expression electric load node i;
17) step (2)~(17) are constantly iteratively repeated and carry out trend iteration, after each trend iteration, record iteration time Number k, and whether trend restrains after judging current iteration.It carries out judging whether trend iteration restrains in the following ways;
If restraining after this trend iteration, the processing result after exporting this trend iteration obtains this trend iteration Voltage V after the trend of each node in power grid afterwardsi' and heat supply network in each node trend after thermal medium flow into temperature T 's.Load.i
If not converged after this trend iteration, using the processing result after this trend iteration as next trend iteration before Numerical value before trend, by voltage V after the trend of each node in the power grid after this trend iterationi' as next trend iteration when Voltage V before the trend of each node in power gridi, thermal medium after the trend of each node in the heat supply network after this trend iteration is flowed into temperature Spend T 's.Load.iThermal medium flows into temperature T before the trend of each node in heat supply network when as next trend iterations.Load.i, carry out down Trend iterative processing.
Heat supply network node thermal medium flows into the maximum value Δ T of temperature difference before and after trend iterations.maxWith grid nodes voltage it The maximum value Δ V of differencei.maxWhether it is respectively less than and is equal to 10-5:
ΔTs.max,ΔVi.max≤10-5 (17)
If meeting above-mentioned formula, restrained after this trend iteration;Otherwise it is not yet restrained after this trend iteration.
After embodiment carries out the processing of electric-thermal association system trend, total the number of iterations is carried out 8 times, obtains final result are as follows: electricity Voltage is followed successively by 211.2405V and 213.0379V, heating power heating network heating power node 2 after the trend of power node 2 and power node 3 Temperature is flowed into successively with thermal medium after the trend of power node 3 are as follows: 69.4444 DEG C and 69.5389.

Claims (5)

1. a kind of electric-thermal association system trend processing method, it is characterised in that:
The electric-thermal association system includes three parts of power grid, heat supply network and cogeneration units, cogeneration units difference Heat supply network and power grid are connected, cogeneration units provide thermal energy and electric energy to heat supply network and power grid simultaneously respectively, according to cogeneration of heat and power machine The power parameter and power output mode of group proportionally carry out the output distribution of electric energy and thermal energy;Heat supply network is divided into heating network and backheat Net, heating network and backheat net are both connected between cogeneration units and heating equipment, the thermal medium conduct of cogeneration units Thermal energy is delivered to heating equipment through heating network, and the thermal medium of heating equipment is delivered to cogeneration of heat and power machine through heating network as thermal energy Group;There are load bus in power grid and heat supply network, the route or pipeline being connected between each node are branch, the load bus in power grid For electric load node, the load bus in heat supply network is heating power load bus, becomes source node, higher level with higher level's supply transformer Supply transformer and cogeneration units, which are conveyed to power grid, provides electric energy;
The method of the present invention including the following steps:
1) the known basis of following power grid and heat supply network is obtained by the sensor in power grid before calculating or by tool acquisition Network data, comprising:
Electric-thermal association system: total node number n, route and duct length Lij, wherein i, j indicate the ordinal number of node, i, j ∈ n;
Power grid: route unit impedance Z, voltage rating VN, electric load node power consumption power PLoad.i
Heat supply network: the rated temperature T of the heat transfer coefficient λ of pipeline unit length, heating networks.N, backheat net rated temperature To.N, heating power The thermal energy consumption power φ of load busLoad.i, heating network rated temperature Ts.N, backheat net rated temperature To.N;Heating power simultaneously The thermal medium of load bus flows out temperature To.Load.iThe as rated temperature T of backheat neto.N, i.e. To.Load.i=To.N;Cogeneration of heat and power The thermal medium of unit flows out temperature To.CHPThe as rated temperature T of heating networks.N, i.e. Ts.N=To.CHP
Cogeneration units: power output the power proportions k, φ of electric energy and thermal energyCHP=k × PCHP, wherein PCHPFor cogeneration units Electric energy go out activity of force, φCHPGo out activity of force for the thermal energy of cogeneration units;
2) according to the following formula, temperature T is flowed into using thermal medium before the trend of heating power load buss.Load.i, heating power load bus Thermal medium flow out temperature To.Load.i, heating power load bus thermal energy consumption power φLoad.iWith the specific heat capacity C of thermal mediump, meter Calculate thermal medium flow m before the trend of each heating power load busq.Load.i
3) then by thermal medium flow m before the trend of each load bus of above-mentioned steps calculating gainedq.Load.iComposition is load bus Thermal medium traffic matrix mq.Load:
It is calculated by the following formula the thermic load traffic matrix m for obtaining heat supply network branch again:
M=Amq.Load (2)
Wherein, A is node branch incidence matrix, and m indicates the thermic load traffic matrix of heat supply network branch, mijTo flow through node in heat supply network The thermic load flow of branch between i and node j;
4) by following formula, the temperature T out of backheat net thermal medium from each branch is obtainedend.H.ij:
Wherein, Tstart.H.ijIndicating the inflow temperature of branch thermal medium in backheat net, λ indicates the heat transfer coefficient of pipeline unit length, mH.ijIndicate the flow of branch thermal medium between backheat net interior joint i and node j, LijIt indicates between heat supply network interior joint i and node j The length of bypass line, e indicate natural Exponents e, are the truth of a matter of natural logrithm;
5) by following formula, the thermal medium for obtaining the node to cross in backheat net there are pipeline flows out temperature Tout.H:
Wherein, Tout.HWith Tin.HIt respectively indicates the thermal medium outflow temperature of the node to cross in backheat net there are pipeline and flows into temperature Degree;mout.HWith min.HIt respectively indicates the thermal medium outflow flow of the node to cross in backheat net there are pipeline and flows into flow;
6) temperature T is flowed out with the thermal medium of the node to cross in backheat net there are pipelineout.HHeat as cogeneration units is situated between Mass flow enters temperature Ts.CHP, i.e. Ts.CHP=Tout.H;Flow m is flowed out with the thermal medium of the node to cross in backheat net there are pipelineout.H Thermal medium as cogeneration units flows through flow mq.CHP, i.e. mq.CHP=mout.H;Then by following formula, according to thermoelectricity The thermal medium of coproduction unit flows into temperature Ts.CHP, cogeneration units thermal medium flow through flow mq.CHPAnd cogeneration units Thermal medium flow out temperature To.CHP, thermal energy goes out activity of force φ before obtaining the trend of cogeneration unitsCHP:
φCHP=Cp·mq.CHP·(Ts.CHP-To.CHP) (5)
7) by following formula, according to the power output power proportions k and cogeneration units of the electric energy of cogeneration units and thermal energy Trend before thermal energy go out activity of force φCHP, obtain electric energy power output power P before the trend of cogeneration unitsCHP:
8) by following formula, connection cogeneration units (CHP) is obtained respectively and without connection cogeneration units (CHP) feelings Under condition, the outflow electric current I of each electric load node in power gridi:
If each electric load node does not connect cogeneration units (CHP) in power grid, according to electric load section each in power grid The power consumption power P of pointLoad.iWith voltage V before the trend of each electric load nodeiThe stream of each electric load node is calculated Electric current I outi:
Wherein, ()*The conjugation of representing matrix calculates;
If each electric load node is connected with cogeneration units (CHP) in power grid, according to electric load node each in power grid Power consumption power PLoad.i, by before the trend of cogeneration units electric energy contribute power PCHPWith the tide of each electric load node Voltage V before flowingiThe outflow electric current I of each electric load node is calculatedi:
9) the outflow electric current I of each electric load node obtained by step 8)iForm the outflow current matrix I of electric load node:
By following formula, power grid branch is obtained according to the outflow current matrix I of electric load node and node branch incidence matrix A The current matrix B on road:
B=AT·I (9)
Wherein, BijThe electric current of branch road between power grid interior joint i and node j;
10) power grid is acquired with direct trend method by the outflow current matrix I of electric load node and node branch incidence matrix A In each electric load node tidal current voltage Vi′;
11) by following formula, for an electric load node for being connected with cogeneration units, according to the electric load section Voltage V after the trend of pointi', the outflow electric current I of the electric load nodeiWith the power consumption power of the electric load node PLoad.i, obtain electric energy power output power P after the trend of cogeneration units 'CHP:
P′CHP=PLoad.i-Ii·Vi′ (10)
Wherein, PLoad.iIndicate the power consumption power of each electric load node in power grid;
12) by following formula, according to the power output power proportions k and cogeneration units of the electric energy of cogeneration units and thermal energy Trend after electric energy power output power P 'CHP, thermal energy goes out activity of force φ ' after obtaining the trend of cogeneration unitsCHP:
φ′CHP=k × P 'CHP (11)
13) by following formula, go out activity of force φ ' using thermal energy after the trend of cogeneration unitsCHP, cogeneration units Thermal medium flows into temperature Ts.CHPTemperature T is flowed out with the thermal medium of cogeneration unitso.CHP, obtain the trend of cogeneration units Thermal medium flows through flow m ' afterwardsq.CHP:
14) flow m ' is flowed through with thermal medium after the trend of cogeneration unitsq.CHPThermal medium as heating network flows out flow mG, Temperature T is flowed out with the thermal medium of cogeneration unitso.CHPThermal medium as heating network flows out temperature TG;Then pass through following public affairs Formula calculates the thermal medium outflow temperature T for obtaining the node to cross in heating network there are pipelineout.G:
Wherein, Tout.GWith Tin.GIt respectively indicates the thermal medium outflow temperature of the node to cross in heating network there are pipeline and flows into temperature Degree;mout.GWith min.GIt respectively indicates the thermal medium outflow flow of the node to cross in heating network there are pipeline and flows into flow;
Temperature T is flowed out with the thermal medium of the node i to cross in heating network there are pipelineout.GIt is situated between as heat after the trend of the node i Mass flow enters temperature T 's.Load.i, flow T is flowed out with the thermal medium of the node i to cross in heating network there are pipelineout.GAs the section Thermal medium flow m ' after the trend that the heating equipment of point i flows throughq.Load.i
15) by following formula, the inflow temperature T of the thermal medium from each branch in heating network is calculatedstart.G.ijWith outflow temperature Tend.G.ij:
Wherein, Tstart.G.ijIndicate the inflow temperature of each branch thermal medium in heating network, λ indicates the heat transfer system of pipeline unit length Number, mH.ijIndicate the flow of branch thermal medium between backheat net interior joint i and node j, LijIndicate heat supply network interior joint i and node j Between bypass line length, e indicate natural Exponents e, be natural logrithm the truth of a matter;
16) by following formula, the maximum value Δ of heat supply network interior joint thermal medium inflow temperature difference before and after trend iteration is calculated Ts.maxWith the maximum value Δ V of the difference of each node voltage of power gridi.max:
ΔTs.max=max (| T 's.Load.i-Ts.Load.i|) (15)
ΔVi.max=max (| Vi′-Vi|) (16)
Wherein, T 's.Load.iThermal medium flows into temperature, T after indicating the trend of the node i to cross in heating network there are pipelines.Load.i Thermal medium flows into temperature before indicating the trend of the node i to cross in heating network there are pipeline;ViIndicate the tide of electric load node i Flow preceding voltage, voltage after the trend of V ' expression electric load node i;
17) step (2)~(16) are constantly iteratively repeated and carry out trend iteration, after each trend iteration, are carried out in the following ways Judge whether trend iteration restrains;
If being restrained after this trend iteration, the processing result after exporting this trend iteration, after obtaining this trend iteration Voltage V after the trend of each node in power gridi' and heat supply network in each node trend after thermal medium flow into temperature T 's.Load.i
If not converged after this trend iteration, using the processing result after this trend iteration as the trend before next trend iteration Preceding numerical value, by voltage V after the trend of each node in the power grid after this trend iterationi' as next trend iteration when power grid In each node trend before voltage Vi, thermal medium after the trend of each node in the heat supply network after this trend iteration is flowed into temperature T′s.Load.iThermal medium flows into temperature T before the trend of each node in heat supply network when as next trend iterations.Load.i, carry out next Secondary trend iterative processing.
2. a kind of electric-thermal association system trend processing method according to claim 1, it is characterised in that:
In the step 17), whether trend iteration convergent to be sentenced and holds back foundation are as follows: heat supply network node thermal medium flows into before and after trend iteration The maximum value Δ T of temperature differences.maxWith the maximum value Δ V of grid nodes difference in voltagei.maxWhether it is respectively less than and is equal to 10-5:
ΔTs.max,ΔVi.max≤10-5 (17)
If meeting above-mentioned formula, restrained after this trend iteration;Otherwise it is not yet restrained after this trend iteration.
3. a kind of electric-thermal association system trend processing method according to claim 1, it is characterised in that:
In the step 17), whether trend iteration convergent to be sentenced and holds back foundation are as follows: it is more than 100 times that trend, which iterates to calculate number k, then Judge that the electric-thermal association system trend does not restrain, terminates electric-thermal association system Load flow calculation.
4. a kind of electric-thermal association system trend processing method according to claim 1, it is characterised in that:
By voltage V before the trend of electric load node each in power gridiIt is initially set to voltage rating VN, by heating power load each in heat supply network Thermal medium flows into temperature T before the trend of nodes.Load.iIt is initially set to the rated temperature T of heating networks.N
5. a kind of electric-thermal association system trend processing method according to claim 1, it is characterised in that:
The node is the equipment that will collect, and distribute and transmit energy (heat energy/electric energy).
CN201811639858.7A 2018-12-29 2018-12-29 Power flow processing method of electric-thermal combined system Active CN109412164B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811639858.7A CN109412164B (en) 2018-12-29 2018-12-29 Power flow processing method of electric-thermal combined system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811639858.7A CN109412164B (en) 2018-12-29 2018-12-29 Power flow processing method of electric-thermal combined system

Publications (2)

Publication Number Publication Date
CN109412164A true CN109412164A (en) 2019-03-01
CN109412164B CN109412164B (en) 2020-06-23

Family

ID=65462765

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811639858.7A Active CN109412164B (en) 2018-12-29 2018-12-29 Power flow processing method of electric-thermal combined system

Country Status (1)

Country Link
CN (1) CN109412164B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110120684A (en) * 2019-05-15 2019-08-13 浙江大学 It is a kind of meter and electric-thermal Demand Side Response cogeneration units configuration method
CN113515853A (en) * 2021-06-08 2021-10-19 国网江苏省电力有限公司连云港供电分公司 Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106056478A (en) * 2016-06-12 2016-10-26 清华大学 Interval power flow calculating method of heat supply network in electrical-thermal coupled system
CN107871058A (en) * 2017-11-24 2018-04-03 清华-伯克利深圳学院筹备办公室 Tidal current computing method, device, equipment and the storage medium of electric heating association system
CN108229753A (en) * 2018-01-30 2018-06-29 清华大学 A kind of distributed scheduling method of coupled thermomechanics system
CN106067677B (en) * 2016-05-27 2018-07-06 清华大学 A kind of electric-thermal couples multipotency streaming system static security analysis method
CN108599266A (en) * 2018-03-21 2018-09-28 浙江大学 A kind of Demand Side Response dispatching method under the hot multipotency stream coupling of electric-gas-

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106067677B (en) * 2016-05-27 2018-07-06 清华大学 A kind of electric-thermal couples multipotency streaming system static security analysis method
CN106056478A (en) * 2016-06-12 2016-10-26 清华大学 Interval power flow calculating method of heat supply network in electrical-thermal coupled system
CN107871058A (en) * 2017-11-24 2018-04-03 清华-伯克利深圳学院筹备办公室 Tidal current computing method, device, equipment and the storage medium of electric heating association system
CN108229753A (en) * 2018-01-30 2018-06-29 清华大学 A kind of distributed scheduling method of coupled thermomechanics system
CN108599266A (en) * 2018-03-21 2018-09-28 浙江大学 A kind of Demand Side Response dispatching method under the hot multipotency stream coupling of electric-gas-

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
谢华宝 等: "电热联合系统潮流计算", 《供用电》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110120684A (en) * 2019-05-15 2019-08-13 浙江大学 It is a kind of meter and electric-thermal Demand Side Response cogeneration units configuration method
CN113515853A (en) * 2021-06-08 2021-10-19 国网江苏省电力有限公司连云港供电分公司 Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation

Also Published As

Publication number Publication date
CN109412164B (en) 2020-06-23

Similar Documents

Publication Publication Date Title
CN109242365A (en) A kind of Interval Power Flow calculation method of electric-thermal interconnection integrated energy system
CN110175311B (en) Optimized power flow calculation method based on multi-energy coupling model
CN109492325B (en) Flow analysis method of multi-energy coupling system based on extended energy concentrator
CN108921727A (en) Consider the regional complex energy resource system reliability estimation method of thermic load dynamic characteristic
Awad et al. Integrated optimal power flow for electric power and heat in a microgrid
CN111061986B (en) Thermoelectric comprehensive energy system tide calculation method with multiple operation modes
CN111191182B (en) Static combined heat and power flow hybrid solving method based on linearized thermal model
CN111428351B (en) Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method
CN101882785B (en) Flow analysis system and method for intelligent power distribution networks containing distributed generation
CN109412164A (en) A kind of electric-thermal association system trend processing method
CN113255105B (en) Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure
CN107257130B (en) Low-voltage distribution network loss calculation method based on regional measurement decoupling
CN112035984B (en) Collaborative planning method for comprehensive energy system of electricity-gas-storage area
CN111859590B (en) Method for constructing heat transport thermal resistance model
CN104600697A (en) Quasi-direct current optimal power flow method considering temperature influence
CN114493258A (en) Power grid carbon footprint tracking method based on electrical subdivision
CN109344480A (en) Consider the electric heating association system tidal current computing method of vapor network
CN108594658A (en) Multi-objective optimization model for probability maximum load margin of electric-gas coupling system and solving method thereof
CN115062555A (en) Integrated energy system exergy flow direct calculation method based on unbalanced node exergy
Man et al. State estimation for integrated energy system containing electricity, heat and gas
CN111310310B (en) Thermodynamic system static state trend rapid decoupling calculation method for quantity adjustment
CN112052543A (en) Bottom-preserving net rack search modeling method based on mixed integer second-order cone programming
Yang et al. A power flow analysis method for the integrated electricity-heat system in distribution network based on forward/backward iterations
CN113515830B (en) Heat supply pipeline network topology transformation-based heat supply network model optimization method
CN111598320B (en) Micro-energy system scheduling method

Legal Events

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