主权项 |
一种电‐热耦合多能流网络节点能价计算方法,其特征在于该方法包括以下步骤:(1)建立一个电‑热耦合多能流网络最优潮流的目标函数,如下:<maths num="0001"><math><![CDATA[<mrow><mi>min</mi><munder><mo>Σ</mo><mrow><mi>i</mi><mo>∈</mo><msub><mi>N</mi><mrow><mi>C</mi><mi>H</mi><mi>P</mi></mrow></msub></mrow></munder><mrow><mo>(</mo><msub><mi>α</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>α</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>α</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>α</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>α</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>α</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>≥</mo><msub><mi>c</mi><mi>l</mi></msub></mrow></mtd><mtd><mrow><mo>∀</mo><mi>l</mi><mo>∈</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>Σ</mi><mrow><mi>j</mi><mo>∈</mo><mi>i</mi></mrow></munder><msub><mi>B</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>θ</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>×</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>λ</mi><mi>L</mi><mo>/</mo><msub><mi>c</mi><mi>p</mi></msub><mover><mi>m</mi><mo>·</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>≤</mo><msub><mi>G</mi><mi>c</mi></msub><mo>≤</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>≤</mo><msub><mi>T</mi><mi>c</mi></msub><mo>≤</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>≤</mo><mi>T</mi><mo>≤</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)中的约束条件构成一个拉格朗日函数,利用原对偶内点法求解该拉格朗日函数,得到电力系统节点功率平衡方程和热力系统节点功率平衡方程所对应的拉格朗日乘子,即分别为电‐热耦合多能流网络各节点的节点电价和节点热价。 |