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
power
node
load
heat
chp
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

本发明公开了一种电‑热联合系统潮流处理方法。针对电‑热联合系统,在计算前通过电网中的传感器或者通过工具采集获得以下电网与热网的已知基础网络数据,根据已知基础网络数据处理获得潮流迭代后的电网中各节点的潮流后电压Vi′和热网中各节点的潮流后热介质流入温度,并且不断潮流处理迭代,每次潮流迭代后,计算潮流迭代前后热网中节点热介质流入温度之差的最大值和电网各节点电压之差的最大值:并判断直至达到要求,实现完成了潮流处理。本发明的电‑热联合系统潮流处理能更快速地进行计算处理,占用内存小,潮流结果准确性好,对于电‑热联合系统的构建能优化运行,向电网调度提供准确数据,提升了直接潮流法的适用性。The invention discloses a power flow processing method for an electric-heat combined system. For the combined power-heat system, the following known basic network data of the power grid and heat network are obtained through sensors in the power grid or through tools before calculation, and the power flow of each node in the power grid after power flow iteration is obtained by processing the known basic network data. The post-voltage V i ′ and the post-power flow temperature of the heat medium in each node in the heat network, and the power flow processing iteration is performed continuously. After each power flow iteration, the maximum value of the difference between the heat medium inflow temperature of the nodes in the heat network before and after the power flow iteration and the power grid are calculated. The maximum value of the difference between the voltages of each node: and judge until the requirements are met, and the power flow processing is completed. The power flow processing of the combined electricity-heating system of the present invention can perform calculation processing more quickly, occupies less memory, and has good accuracy of power flow results, can optimize the operation for the construction of the combined electricity-heating system, provide accurate data for grid scheduling, and improve direct Applicability of the tide 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.一种电-热联合系统潮流处理方法,其特征在于:1. an electric-heat combined system power flow processing method is characterized in that: 所述的电-热联合系统包括电网、热网和热电联产机组的三个部分,热电联产机组分别连接热网和电网,热电联产机组向热网和电网分别同时提供热能与电能,依据热电联产机组的功率参数和出力方式按照比例进行电能与热能的出力分配;热网分为供热网和回热网,供热网和回热网均连接在热电联产机组和用热设备之间,热电联产机组的热介质作为热能经供热网输送至用热设备,用热设备的热介质作为热能经供热网输送至热电联产机组;电网和热网中存在负荷节点,各节点之间相连的线路或管道为支路,电网中的负荷节点为电力负荷节点,热网中的负荷节点为热力负荷节点,以上级供电变压器成为源节点,上级供电变压器与热电联产机组均向电网输送提供电能;The power-heat combined system includes three parts: a power grid, a heat grid and a combined heat and power unit. The combined heat and power unit is connected to the heat grid and the power grid, respectively, and the combined heat and power unit provides thermal energy and electric energy to the heat grid and the power grid at the same time, respectively. According to the power parameters and output mode of the cogeneration unit, the output of electric energy and heat energy is distributed in proportion; the heat network is divided into a heating network and a heat recovery network, and the heating network and the heat recovery network are connected to the cogeneration unit and the heat consumption network. Between equipments, the heat medium of the cogeneration unit is transported as heat energy to the heating equipment through the heating network, and the heat medium of the heating equipment is transported as heat energy to the cogeneration unit through the heating network; there are load nodes in the power grid and the heating network. , the line or pipeline connected between each node is a branch, the load node in the power grid is an electric load node, the load node in the heat network is a thermal load node, the upper-level power supply transformer becomes the source node, and the upper-level power supply transformer and cogeneration All units supply power to the grid; 本发明方法包括以下几个步骤:The method of the present invention comprises the following steps: 1)在计算前通过电网中的传感器或者通过工具采集获得以下电网与热网的已知基础网络数据,包括:1) Before the calculation, the known basic network data of the following power grids and heat grids, including: 电-热联合系统:总节点数n、线路与管道长度Lij,其中i,j均表示节点的序数,i,j∈n;Electric-heat combined system: the total number of nodes n, the length of lines and pipelines L ij , where i, j both represent the ordinal numbers of nodes, i, j∈n; 电网:线路单位阻抗Z、额定电压VN、电力负荷节点的电能消耗功率PLoad.iPower grid: line unit impedance Z, rated voltage V N , power consumption power P Load.i of power load nodes; 热网:管道单位长度的传热系数λ、供热网的额定温度Ts.N、回热网的额定温度To.N、热力负荷节点的热能消耗功率φLoad.i、供热网的额定温度Ts.N、回热网的额定温度To.N;同时热力负荷节点的热介质流出温度To.Load.i即为回热网的额定温度To.N,即To.Load.i=To.N;热电联产机组的热介质流出温度To.CHP即为供热网的额定温度Ts.N,即Ts.N=To.CHPHeating network: the heat transfer coefficient λ per unit length of the pipeline, the rated temperature T sN of the heating network, the rated temperature T oN of the regenerative heating network, the thermal energy consumption power φ Load.i of the thermal load node, and the rated temperature T sN of the heating network , the rated temperature T oN of the regenerative network; at the same time, the heat medium outflow temperature T o.Load.i of the thermal load node is the rated temperature T oN of the regenerative network, that is, T o.Load.i =T oN ; cogeneration The heat medium outflow temperature T o.CHP of the unit is the rated temperature T sN of the heating network, that is, T sN =T o.CHP ; 热电联产机组:电能与热能的出力功率比例k,φCHP=k×PCHP,其中PCHP为热电联产机组的电能出力功率,φCHP为热电联产机组的热能出力功率;Cogeneration unit: output power ratio k of electric energy and thermal energy, φ CHP = k×P CHP , where P CHP is the electrical output power of the cogeneration unit, and φ CHP is the thermal output power of the cogeneration unit; 2)根据以下公式,利用热力负荷节点的潮流前热介质流入温度Ts.Load.i、热力负荷节点的热介质流出温度To.Load.i、热力负荷节点的热能消耗功率φLoad.i和热介质的比热容Cp,计算出每个热力负荷节点的潮流前热介质流量mq.Load.i2) According to the following formula, the inflow temperature T s.Load.i of the thermal medium before the power flow of the thermal load node, the outflow temperature T o.Load.i of the thermal medium of the thermal load node, and the thermal energy consumption power φ Load.i of the thermal load node are used and the specific heat capacity C p of the heat medium, calculate the heat medium flow m q.Load.i before the power flow of each heat load node; 3)接着由上述步骤计算所得各负荷节点的潮流前热介质流量mq.Load.i组成为负荷节点的热介质流量矩阵mq.Load3) Next, the heat medium flow m q.Load.i before the power flow of each load node calculated by the above steps is composed of the heat medium flow matrix m q.Load of the load node: 再通过以下公式计算获得热网支路的热负荷流量矩阵m:Then calculate the heat load flow matrix m of the heat network branch by the following formula: m=A·mq.Load (2)m=A m q. Load (2) 其中,A为节点支路关联矩阵,m表示热网支路的热负荷流量矩阵,mij为热网中流经节点i与节点j之间的支路的热负荷流量;Among them, A is the node branch correlation matrix, m represents the heat load flow matrix of the heat network branch, m ij is the heat load flow through the branch between node i and node j in the heat network; 4)通过以下公式,获得回热网从各支路中热介质的出温度Tend.H.ij4) Obtain the outlet temperature T end.H.ij of the heat medium from each branch of the regenerative network by the following formula: 其中,Tstart.H.ij表示回热网中支路热介质的流入温度,λ表示管道单位长度的传热系数,mH.ij表示回热网中节点i和节点j之间支路热介质的流量,Lij表示热网中节点i和节点j之间支路管道的长度,e表示自然指数e,为自然对数的底数;Among them, T start.H.ij represents the inflow temperature of the branch heat medium in the regenerative network, λ represents the heat transfer coefficient per unit length of the pipeline, and m H.ij represents the branch heat between node i and node j in the regenerative network. The flow rate of the medium, Li ij represents the length of the branch pipe between node i and node j in the heat network, e represents the natural exponent e, which is the base of the natural logarithm; 5)通过以下公式,获得回热网中存在管道交汇的节点的热介质流出温度Tout.H5) Obtain the outflow temperature T out.H of the heat medium at the node where the pipes intersect in the heat recovery network by the following formula: 其中,Tout.H与Tin.H分别表示回热网中存在管道交汇的节点的热介质流出温度和流入温度;mout.H与min.H分别表示回热网中存在管道交汇的节点的热介质流出流量和流入流量;Among them, T out.H and T in.H represent the outflow temperature and inflow temperature of the heat medium at the nodes where the pipes intersect in the regenerative network , respectively; The outflow and inflow of the thermal medium of the node; 6)以回热网中存在管道交汇的节点的热介质流出温度Tout.H作为热电联产机组的热介质流入温度Ts.CHP,即Ts.CHP=Tout.H;以回热网中存在管道交汇的节点的热介质流出流量mout.H作为热电联产机组的热介质流经流量mq.CHP,即mq.CHP=mout.H;然后通过以下公式,根据热电联产机组的热介质流入温度Ts.CHP、热电联产机组的热介质流经流量mq.CHP和热电联产机组的热介质流出温度To.CHP,获得热电联产机组的潮流前热能出力功率φCHP6) Take the heat medium outflow temperature T out.H of the node where the pipes intersect in the heat recovery network as the heat medium inflow temperature T s.CHP of the cogeneration unit, that is, T s.CHP =T out.H ; The outflow flow m out.H of the heat medium at the nodes where the pipes intersect in the network is taken as the flow m q.CHP of the heat medium of the cogeneration unit, that is, m q.CHP = m out.H ; then through the following formula, according to the thermoelectricity The inflow temperature T s.CHP of the heat medium of the cogeneration unit, the flow rate m q.CHP of the heat medium of the cogeneration unit and the outflow temperature of the heat medium T o.CHP of the heat and power unit, and the flow front of the cogeneration unit is obtained. Thermal output power φ CHP : φCHP=Cp·mq.CHP·(Ts.CHP-To.CHP) (5)φ CHP =C p · m q . CHP · (T s. CHP -T o. CHP ) (5) 7)通过以下公式,根据热电联产机组的电能与热能的出力功率比例k与热电联产机组的潮流前热能出力功率φCHP,获得热电联产机组的潮流前电能出力功率PCHP7) According to the following formula, according to the output power ratio k of the electric energy and thermal energy of the cogeneration unit and the pre-tidal flow thermal output power φ CHP of the cogeneration unit, the pre-tidal electric power output P CHP of the cogeneration unit is obtained: 8)通过以下公式,分别获得连接热电联产机组(CHP)和没有连接热电联产机组(CHP)情况下,电网中各电力负荷节点的流出电流Ii8) Obtain the outflow current I i of each power load node in the power grid when the combined heat and power unit (CHP) is connected and the combined heat and power unit (CHP) is not connected by the following formula: 若电网中各电力负荷节点没有连接热电联产机组(CHP),则根据电网中各电力负荷节点的电能消耗功率PLoad.i与各电力负荷节点的潮流前电压Vi计算得到各电力负荷节点的流出电流IiIf each power load node in the power grid is not connected to a combined heat and power unit (CHP), then each power load node is calculated according to the power consumption power P Load.i of each power load node in the power grid and the pre-power flow voltage V i of each power load node The outflow current I i : 其中,()*表示矩阵的共轭计算;Among them, () * represents the conjugate calculation of the matrix; 若电网中各电力负荷节点连接有热电联产机组(CHP),则根据电网中各电力负荷节点的电能消耗功率PLoad.i、由热电联产机组的潮流前电能出力功率PCHP与各电力负荷节点的潮流前电压Vi计算得到各电力负荷节点的流出电流IiIf each power load node in the power grid is connected to a combined heat and power unit ( CHP ), then according to the power consumption power P Load . The outflow current I i of each power load node is obtained by calculating the pre-flow voltage V i of the load node: 9)由步骤8)获得的各电力负荷节点的流出电流Ii组成电力负荷节点的流出电流矩阵I:9) The outflow current I i of each power load node obtained in step 8) forms the outflow current matrix I of the power load node: 通过以下公式,根据电力负荷节点的流出电流矩阵I与节点支路关联矩阵A获得电网支路的电流矩阵B:Through the following formula, the current matrix B of the grid branch is obtained according to the outflow current matrix I of the power load node and the node branch correlation matrix A: B=AT·I (9)B=A T · I (9) 其中,Bij为电网中节点i和节点j之间支路上的电流;Among them, B ij is the current on the branch between node i and node j in the power grid; 10)由电力负荷节点的流出电流矩阵I与节点支路关联矩阵A运用直接潮流法求得电网中各电力负荷节点的潮流电压Vi′;10) Using the direct power flow method to obtain the power flow voltage V i ′ of each power load node in the power grid from the outflow current matrix I of the power load node and the node branch correlation matrix A; 11)通过以下公式,针对连接有热电联产机组的一个电力负荷节点,根据该电力负荷节点的潮流后电压Vi′、该电力负荷节点的流出电流Ii与该电力负荷节点的电能消耗功率PLoad.i,获得热电联产机组的潮流后电能出力功率P′CHP11) Through the following formula, for a power load node connected with a cogeneration unit, according to the post-power flow voltage V i ′ of the power load node, the outflow current I i of the power load node and the power consumption power of the power load node P Load.i , obtain the post-tidal power output power P′ CHP of the cogeneration unit: P′CHP=PLoad.i-Ii·Vi′ (10)P' CHP =P Load.i -I i ·V i ' (10) 其中,PLoad.i表示电网中各电力负荷节点的电能消耗功率;Among them, P Load.i represents the electric energy consumption power of each electric load node in the power grid; 12)通过以下公式,根据热电联产机组的电能与热能的出力功率比例k与热电联产机组的潮流后电能出力功率P′CHP,获得热电联产机组的潮流后热能出力功率φ′CHP12) Through the following formula, according to the output power ratio k of the electric energy and thermal energy of the cogeneration unit and the electric power output P′ CHP of the cogeneration unit after the tidal current, the thermal output power φ′ CHP of the cogeneration unit after the tidal current is obtained: φ′CHP=k×P′CHP (11)φ′ CHP = k×P′ CHP (11) 13)通过以下公式,利用热电联产机组的潮流后热能出力功率φ′CHP、热电联产机组的热介质流入温度Ts.CHP和热电联产机组的热介质流出温度To.CHP,获得热电联产机组的潮流后热介质流经流量m′q.CHP13) Through the following formula, using the post-tidal heat output power φ′ CHP of the cogeneration unit, the inflow temperature T s.CHP of the heat medium of the cogeneration unit and the outflow temperature T o.CHP of the heat medium of the cogeneration unit to obtain After the power flow of the cogeneration unit, the heat medium flows through the flow m' q.CHP : 14)以热电联产机组的潮流后热介质流经流量m′q.CHP作为供热网的热介质流出流量mG,以热电联产机组的热介质流出温度To.CHP作为供热网的热介质流出温度TG;然后通过以下公式,计算获得供热网中存在管道交汇的节点的热介质流出温度Tout.G14) The flow rate m′ q.CHP of the heat medium after the power flow of the cogeneration unit is taken as the outflow flow m G of the heat medium of the heating network, and the outflow temperature of the heat medium T o.CHP of the heat and power unit is taken as the heat supply network The outflow temperature of the heat medium T G ; then the outflow temperature T out.G of the heat medium at the node where the pipes intersect in the heating network is calculated by the following formula: 其中,Tout.G与Tin.G分别表示供热网中存在管道交汇的节点的热介质流出温度和流入温度;mout.G与min.G分别表示供热网中存在管道交汇的节点的热介质流出流量和流入流量;Among them, T out.G and T in.G respectively represent the outflow temperature and inflow temperature of the heat medium at the nodes where the pipes intersect in the heating network; m out.G and min.G respectively represent the pipes intersecting in the heating network. The outflow and inflow of the thermal medium of the node; 以供热网中存在管道交汇的节点i的热介质流出温度Tout.G作为该节点i的潮流后热介质流入温度T′s.Load.i,以供热网中存在管道交汇的节点i的热介质流出流量Tout.G作为该节点i的用热设备流经的潮流后热介质流量m′q.Load.iTake the outflow temperature T out.G of the heat medium at the node i where the pipes intersect in the heating network as the inflow temperature T′ s.Load.i of the heat medium after the power flow of the node i, and take the node i where the pipes intersect in the heating network The heat medium outflow flow T out.G is taken as the heat medium flow m′ q.Load.i after the power flow of the heat-using equipment at the node i; 15)通过以下公式,计算供热网中从各支路中热介质的流入温度Tstart.G.ij与流出温度Tend.G.ij15) Calculate the inflow temperature T start.G.ij and outflow temperature T end.G.ij of the heat medium from each branch in the heating network by the following formula: 其中,Tstart.G.ij表示供热网中各支路热介质的流入温度,λ表示管道单位长度的传热系数,mH.ij表示回热网中节点i和节点j之间支路热介质的流量,Lij表示热网中节点i和节点j之间支路管道的长度,e表示自然指数e,为自然对数的底数;Among them, T start.G.ij represents the inflow temperature of the heat medium in each branch in the heating network, λ represents the heat transfer coefficient per unit length of the pipeline, and m H.ij represents the branch between node i and node j in the heating network. The flow rate of the heat medium, L ij represents the length of the branch pipe between node i and node j in the heat network, e represents the natural exponent e, which is the base of the natural logarithm; 16)通过以下公式,计算潮流迭代前后热网中节点热介质流入温度之差的最大值ΔTs.max和电网各节点电压之差的最大值ΔVi.max16) Calculate the maximum value ΔT s.max of the difference between the inflow temperature of the heat medium at the nodes in the thermal network before and after the power flow iteration and the maximum value ΔV i.max of the difference between the voltages of each node in the power grid by the following formula: ΔTs.max=max(|T′s.Load.i-Ts.Load.i|) (15)ΔT s.max =max(|T′ s.Load.i -T s.Load.i |) (15) ΔVi.max=max(|Vi′-Vi|) (16)ΔV i.max =max(|V i ′-V i |) (16) 其中,T′s.Load.i表示供热网中存在管道交汇的节点i的潮流后热介质流入温度,Ts.Load.i表示供热网中存在管道交汇的节点i的潮流前热介质流入温度;Vi表示电力负荷节点i的潮流前电压,V′表示电力负荷节点i的潮流后电压;Among them, T' s.Load.i represents the inflow temperature of the heat medium after the power flow at the node i where the pipelines intersect in the heating network, and T s.Load.i represents the heat medium before the power flow at the node i where the pipelines intersect in the heating network Inflow temperature; V i represents the pre-power flow voltage of power load node i, and V′ represents the post-power flow voltage of power load node i; 17)不断迭代重复步骤(2)~(16)进行潮流迭代,每次潮流迭代后,采用以下方式进行判断潮流迭代是否收敛;17) Repeat steps (2) to (16) continuously for power flow iteration. After each power flow iteration, use the following method to judge whether the power flow iteration has converged; 若本次潮流迭代后收敛,则输出本次潮流迭代后的处理结果,获得本次潮流迭代后的电网中各节点的潮流后电压Vi′和热网中各节点的潮流后热介质流入温度T′s.Load.iIf the power flow is converged after the current power flow iteration, the processing result after the power flow iteration is output, and the post-power flow voltage V i ′ of each node in the power grid and the post-power flow temperature of the thermal medium in the heat network after this power flow iteration are obtained. T's.Load.i ; 若本次潮流迭代后未收敛,将本次潮流迭代后的处理结果作为下次潮流迭代前的潮流前数值,将本次潮流迭代后的电网中各节点的潮流后电压Vi′作为下次潮流迭代时的电网中各节点的潮流前电压Vi,将本次潮流迭代后的热网中各节点的潮流后热介质流入温度T′s.Load.i作为下次潮流迭代时的热网中各节点的潮流前热介质流入温度Ts.Load.i,进行下一次潮流迭代处理。If the power flow does not converge after this power flow iteration, the processing result after this power flow iteration is taken as the pre-power flow value before the next power flow iteration, and the post-power flow voltage V i ′ of each node in the power grid after this power flow iteration is taken as the next power flow iteration The pre-power flow voltage V i of each node in the power grid during the power flow iteration, and the post-power flow thermal medium inflow temperature T′ s.Load.i of each node in the thermal network after the current power flow iteration is used as the thermal network in the next power flow iteration The inflow temperature T s.Load.i of the heat medium before the power flow of each node in , and the next power flow iteration processing is performed. 2.根据权利要求1所述的一种电-热联合系统潮流处理方法,其特征在于:2. a kind of electric-heat combined system power flow processing method according to claim 1 is characterized in that: 所述步骤17)中,潮流迭代是否收敛的判敛依据为:潮流迭代前后热网节点热介质流入温度之差的最大值ΔTs.max与电网节点电压之差的最大值ΔVi.max是否均小于等于10-5In the step 17), the basis for judging whether the power flow iteration has converged is: whether the maximum value ΔT s . are less than or equal to 10 -5 : ΔTs.max,ΔVi.max≤10-5 (17)ΔT s.max ,ΔV i.max ≤10 -5 (17) 如果满足上述公式,则本次潮流迭代后收敛;否则本次潮流迭代后尚未收敛。If the above formula is satisfied, the power flow has converged after this iteration; otherwise, it has not converged after the current power flow iteration. 3.根据权利要求1所述的一种电-热联合系统潮流处理方法,其特征在于:3. a kind of electric-heat combined system power flow processing method according to claim 1 is characterized in that: 所述步骤17)中,潮流迭代是否收敛的判敛依据为:潮流迭代计算次数k超过100次,则判断该电-热联合系统潮流不收敛,结束电-热联合系统潮流计算。In the step 17), the convergence basis for determining whether the power flow iteration is convergent is: if the number of power flow iteration calculations k exceeds 100, it is judged that the power flow of the combined electric-thermal system does not converge, and the power flow calculation of the combined electric-thermal system is terminated. 4.根据权利要求1所述的一种电-热联合系统潮流处理方法,其特征在于:4. a kind of electric-heat combined system power flow processing method according to claim 1, is characterized in that: 将电网中各电力负荷节点的潮流前电压Vi初始设为额定电压VN,将热网中各热力负荷节点的潮流前热介质流入温度Ts.Load.i初始设为供热网的额定温度Ts.NThe pre-load flow voltage V i of each power load node in the power grid is initially set as the rated voltage V N , and the pre-load flow temperature T s.Load.i of the thermal medium in the heating network is initially set as the rated voltage of the heating network. Temperature T sN . 5.根据权利要求1所述的一种电-热联合系统潮流处理方法,其特征在于:5. a kind of electric-heat combined system power flow processing method according to claim 1, is characterized in that: 所述的节点为将汇集、分配和传送能量(热能/电能)的设备。Said nodes are devices that will collect, distribute and transmit energy (thermal/electrical).
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 清华-伯克利深圳学院筹备办公室 Power flow calculation method, device, equipment and storage medium of electric heating combined 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 清华-伯克利深圳学院筹备办公室 Power flow calculation method, device, equipment and storage medium of electric heating combined 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
Ayele et al. An extended energy hub approach for load flow analysis of highly coupled district energy networks: Illustration with electricity and heating
CN106339772B (en) Combined heat-electricity optimal scheduling method based on heat storage efficiency of heating pipe network
CN109242365A (en) A kind of Interval Power Flow calculation method of electric-thermal interconnection integrated energy system
CN111061986B (en) Thermoelectric comprehensive energy system tide calculation method with multiple operation modes
CN111815042A (en) Optimal scheduling method of electric-heat integrated energy system considering refined heat network model
Awad et al. Integrated optimal power flow for electric power and heat in a microgrid
CN109492325B (en) Flow analysis method of multi-energy coupling system based on extended energy concentrator
CN111428351B (en) Power flow calculation method for electric-thermal integrated energy system based on forward push-back method
CN113255105B (en) Power flow calculation method for electric and thermal integrated energy system with bidirectional coupling network structure
CN104113061B (en) A kind of distribution network three-phase power flow method containing distributed power source
CN111555285B (en) Energy flow decoupling analysis and calculation method for distributed combined cooling heating and power comprehensive energy system
CN111191182A (en) Hybrid solution method of combined power flow of static heat and electricity based on linearized thermodynamic model
CN107257130B (en) Loss calculation method of low-voltage distribution network based on regional measurement decoupling
CN109389248A (en) Integrated energy system thermoelectric (al) power coordinated scheduling method based on region heat supply network
CN111859590B (en) Method for constructing heat transport thermal resistance model
CN109412164A (en) A kind of electric-thermal association system trend processing method
CN111310310B (en) Thermodynamic system static state trend rapid decoupling calculation method for quantity adjustment
CN110688744A (en) An Asynchronous Distributed State Estimation Method for Thermoelectric Coupling Networks
CN113324283B (en) Central heating network equivalence method considering temperature semi-dynamic characteristics
CN112052543B (en) A search modeling method for bottom-guaranteed grids based on mixed integer second-order cone programming
CN108418211A (en) A Power Supply Path Optimization Method Based on Dijlstra Algorithm and Genetic Algorithm
CN102437643B (en) Method for determining available power transmission capacity of each power transmission section in multiple sections
CN116611706A (en) Dynamic carbon emission factor measuring and calculating method based on multi-energy main body
CN114549232A (en) Hybrid energy flow calculation method for electricity-heat comprehensive energy system
CN114566972A (en) Composite power flow tracking line loss optimization method considering new energy fluctuation uncertainty

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