发明名称 一种电‑热耦合多能流网络节点能价计算方法
摘要 本发明涉及一种电‐热耦合多能流网络节点能价方法,属于综合能源系统中的能源市场领域。本方法考虑了电热耦合网络中电价和热价的耦合关系,打破了热电联供机组联合成本定价只能解耦计算的现状,还原了电热耦合系统中的市场互动机制,相较于现有的对于供电、供热系统分别进行定价的方法,提高了边际成本计算的精确性。同时,节点能价作为实时价格中的一种,弥补了热网中实时热价的空白,为供热网络进行阻塞管理和网损分摊提供了一定的思路。
申请公布号 CN106339794A 申请公布日期 2017.01.18
申请号 CN201610674418.X 申请日期 2016.08.16
申请人 清华大学 发明人 孙宏斌;郭庆来;王彬;邓莉荣;潘昭光
分类号 G06Q10/06(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/06(2012.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 一种电‐热耦合多能流网络节点能价计算方法,其特征在于该方法包括以下步骤:(1)建立一个电‑热耦合多能流网络最优潮流的目标函数,如下:<maths num="0001"><math><![CDATA[<mrow><mi>min</mi><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>N</mi><mrow><mi>C</mi><mi>H</mi><mi>P</mi></mrow></msub></mrow></munder><mrow><mo>(</mo><msub><mi>&alpha;</mi><mrow><mn>5</mn><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>p</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><msub><mi>q</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>&alpha;</mi><mrow><mn>4</mn><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>q</mi><mrow><mi>G</mi><mi>i</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>&alpha;</mi><mrow><mn>3</mn><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>p</mi><mrow><mi>G</mi><mi>i</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>&alpha;</mi><mrow><mn>2</mn><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>q</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>&alpha;</mi><mrow><mn>1</mn><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>p</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>&alpha;</mi><mrow><mn>0</mn><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001079585300000011.GIF" wi="1390" he="143" /></maths>其中:i为电‑热耦合多能流网络中热电机组的编号,N<sub>CHP</sub>为热电机组的数量,α<sub>0,i</sub>~α<sub>5,i</sub>为第i台热电机组的燃料费用系数,α<sub>0,i</sub>~α<sub>5,i</sub>从第i台热电机组的产品说明书中获取,p<sub>Gi</sub>、q<sub>Gi</sub>分别为第i台热电机组的电出力和热出力;(2)建立一个电‐热耦合多能流网络最优潮流的约束条件,包括:(2‐1)电‐热耦合多能流网络中的热电联供机组可行域约束,如下:<maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>a</mi><mi>l</mi></msub><msub><mi>p</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>b</mi><mi>l</mi></msub><msub><mi>q</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>&GreaterEqual;</mo><msub><mi>c</mi><mi>l</mi></msub></mrow></mtd><mtd><mrow><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><msub><mi>N</mi><mi>l</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001079585300000012.GIF" wi="558" he="55" /></maths>其中:l为组成热电联供机组可行域的线段编号,N<sub>l</sub>为线段总数,a<sub>l</sub>,b<sub>l</sub>,c<sub>l</sub>为各条线段参数,a<sub>l</sub>,b<sub>l</sub>,c<sub>l</sub>从热电联供机组的产品说明书中获取,<img file="FDA0001079585300000016.GIF" wi="41" he="41" />表示任意,<img file="FDA0001079585300000017.GIF" wi="56" he="50" />表示对于任意的l都有上述的线性约束;(2‐2)电‑热耦合多能流网络中电网采用直流潮流模型,直流潮流模型中的有功功率平衡约束如下:<maths num="0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>B</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>&theta;</mi><mi>j</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001079585300000013.GIF" wi="662" he="143" /></maths>其中:P<sub>i</sub>为电网中节点i的注入有功功率,θ<sub>j</sub>为电网中节点j的电压相角,B<sub>ij</sub>为电网节点导纳矩阵Y的第i行、第j列元素的虚部,电网节点导纳矩阵Y从电‑热耦合多能流系统的能量管理系统中获取;(2‐3)电‑热耦合多能流网络中电网的传输功率约束: P<sub>ij</sub>=‑B<sub>ij</sub>(θ<sub>i</sub>‑θ<sub>j</sub>)≤P<sub>ijmax</sub>其中:P<sub>ij</sub>为电网中节点i到节点j的线路传输有功功率,P<sub>ijmax</sub>为电网中节点i到节点j的线路传输有功功率上限;(2‑4)电‑热耦合多能流网络中热网的热力平衡方程约束:<maths num="0004"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><msub><mi>A</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><msub><mi>G</mi><mi>c</mi></msub><msubsup><mi>A</mi><mrow><mi>a</mi><mn>1</mn></mrow><mi>T</mi></msubsup><mi>T</mi><mo>-</mo><msub><mi>A</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><msub><mi>G</mi><mi>c</mi></msub><msub><mi>T</mi><mi>c</mi></msub></mrow>]]></math><img file="FDA0001079585300000014.GIF" wi="590" he="69" /></maths> G<sub>0</sub>=A<sub>a</sub>·G<sub>c</sub><maths num="0005"><math><![CDATA[<mrow><msub><mi>T</mi><mi>c</mi></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mi>A</mi><mrow><mi>a</mi><mn>1</mn></mrow><mi>T</mi></msubsup><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&lambda;</mi><mi>L</mi><mo>/</mo><msub><mi>c</mi><mi>p</mi></msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>+</mo><msub><mi>T</mi><mi>a</mi></msub></mrow>]]></math><img file="FDA0001079585300000015.GIF" wi="846" he="75" /></maths>其中:Q为热网节点的注入热功率列向量,G<sub>c</sub>为热网管段流量热当量对角阵,G<sub>0</sub>为热网节点流量热当量列向量,T为热网节点流出温度列向量,T<sub>c</sub>为热网管段出口温度列向量,T<sub>a</sub>为热网管段所在的环境温度,c<sub>p</sub>为热网介质水的比热容,取值为4182焦耳/(千克·摄氏度),L为热网管段长度,λ为热网管段单位长度的导热系数列向量,L和λ从电‐热耦合多能流系统的能量管理系统中获取,<img file="FDA0001079585300000021.GIF" wi="49" he="47" />为热网管段的介质水质量流量,A<sub>a1</sub>为热网管段的起点增广关联矩阵,A<sub>a2</sub>为热网管段的终点增广关联矩阵,A<sub>a</sub>为热网管段的增广关联矩阵,A<sub>a1</sub>、A<sub>a2</sub>和A<sub>a</sub>从电‐热耦合多能流系统的能量管理系统中获取,T为矩阵转置;(2‑5)电‑热耦合多能流网络中热网的运行约束:<maths num="0006"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>G</mi><mrow><mi>c</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>G</mi><mi>c</mi></msub><mo>&le;</mo><msub><mi>G</mi><mrow><mi>c</mi><mi>max</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>c</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>T</mi><mi>c</mi></msub><mo>&le;</mo><msub><mi>T</mi><mrow><mi>c</mi><mi>max</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mi>min</mi></msub><mo>&le;</mo><mi>T</mi><mo>&le;</mo><msub><mi>T</mi><mi>max</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001079585300000022.GIF" wi="446" he="175" /></maths>其中:G<sub>cmin</sub>和G<sub>cmax</sub>分别为热网管段流量热当量的运行上限和下限,T<sub>cmin</sub>和T<sub>cmax</sub>为分别热网管段出口温度的上限和下限,T<sub>min</sub>和T<sub>max</sub>分别为热网节点流出温度的上限和下限,T<sub>cmin</sub>、T<sub>cmax</sub>、T<sub>min</sub>和T<sub>max</sub>分别从电‐热耦合多能流系统的能量管理系统中获取;(3)利用拉格朗日乘子法,将步骤(1)中的目标函数和步骤(2)中的约束条件构成一个拉格朗日函数,利用原对偶内点法求解该拉格朗日函数,得到电力系统节点功率平衡方程和热力系统节点功率平衡方程所对应的拉格朗日乘子,即分别为电‐热耦合多能流网络各节点的节点电价和节点热价。
地址 100084 北京市海淀区清华园1号