CN106099940B - 变电站内变压器分接头档位对母线电压影响的计算方法 - Google Patents
变电站内变压器分接头档位对母线电压影响的计算方法 Download PDFInfo
- Publication number
- CN106099940B CN106099940B CN201610692725.0A CN201610692725A CN106099940B CN 106099940 B CN106099940 B CN 106099940B CN 201610692725 A CN201610692725 A CN 201610692725A CN 106099940 B CN106099940 B CN 106099940B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- substation
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种变电站内变压器分接头档位对母线电压影响的计算方法，属于电力系统中电压自动控制技术领域。本方法在实时的潮流计算基础上，求出采用准稳态潮流构造变压器绕组有载分头变压比发生变化与变压器绕组支路潮流变化的线性方程，求出采用准稳态潮流构造注入潮流变化量与母线电压变化量的线性方程，计算变压器绕组有载调压分头变压比发生变化后对电网模型中所有母线的电压幅值变化值，找出与变压器分接头同属于一个变电站的母线电压变化量。本方法给出了变电站变压器分接头档位对母线电压影响的灵敏度结果，在变电站电压自动控制的过程中，对变压器分接头档位发生变化后各条母线的电压进行预估，保证变电站变压器档位控制策略可靠性。
Description
技术领域
本发明涉及一种变电站内变压器分接头档位对母线电压影响的计算方法，属于电力系统中电压自动控制技术领域。
背景技术
自动电压控制(以下简称AVC，Automatic Voltage Control)系统是电力系统最重要的自动控制系统之一。其中变电站的自动电压控制能力与效果直接影响了电力系统自动电压控制的整体控制结果。变电站里均配有一定容量的离散无功设备(电容器，电抗器)来完成自动电压控制的目标。但由于变电站均配有多台电容器，电抗器，且在选择投入或切除具体某台电容器或电抗器时，需考虑投入(切除)不同的电容电抗器对不同母线电压的效果不同，且每台电容器，电抗器都有不同的动作时间，动作次数的限制，以及投入(切除)无功设备的相关顺序约束等。
现有的变电站控制系统中在接到需投入(切除)电容电抗器的控制要求后，利用传统经验来决定具体由哪台无功设备来执行控制任务，缺乏量化，可靠的决策标准，难以满足复杂，特殊情况下的控制要求，并且其固化的评判标准不具备应对在不同控制策略下，对离散无功设备的各项性能指标的不同要求。
变电站内投入(切除)电容电抗器应具有如下原则：设备操作后，通过计算该设备对母线的电压灵敏度，可以算出投入(切除)该设备后，各条母线的电压变化量，所投入(切除)的设备应使母线电压越线数减少或消除电压越线现象；优先选择对受控母线的无功影响大，即灵敏度大的设备；优先投入(切除)越线母线上的无功设备。在考虑主变功率因数越限策略时，根据设备对主变绕组无功灵敏度，可以算出投入(切除)该设备后，各绕组首端无功变化量，所投入(切除)的设备应使主变功率因越限数减少或消除功率因数越限现象，且投入(切除)的设备应满足主变绕组无功倒送阀值约束的要求；优先选择对受控母线的无功影响大，即灵敏度大的设备。
在上述变电站控制过程中涉及到控制灵敏度的计算。孙宏斌，张伯明，相年德在《准稳态的灵敏度分析方法》(中国电机工程学报，1999年4月V19N4，pp.9-13)中提出了准稳态灵敏度方法，与常规的静态的灵敏度分析方法不同，准稳态灵敏度方法考虑了电力系统准稳态的物理响应，计及系统控制前后新旧稳态间的总变化，有效提高了灵敏度分析的精度。该方法基于电力系统的PQ解耦模型，当发电机安装有自动电压调节器(AVR)时，可认为该发电机节点为PV节点；而当发电机装有自动无功功率调节(AQR)或自动功率因数调节(APFR)时，可认为该发电机节点与普通负荷节点相同均为PQ节点。此外，将负荷电压静特性考虑成节点电压的一次或二次曲线。这样所建立的潮流模型就自然地将这些准稳态的物理响应加以考虑，从而基于潮流模型计算出的灵敏度即为准稳态的灵敏度。
电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的数据，计算母线的电压、各元件的功率及网损，并对电网各处的运行状态进行评估。再根据计算得到的数据对电网系统的运行进行监测和优化，从而提高供电方案或运行方式的合理性、可靠性和经济性。对于正在运行的电力系统，通过潮流计算可以评估当前系统中母线的电压、支路的功率等参数是否超限；如果出现异常，就应采取措施，调整运行方式。
极坐标下的电网潮流方程可以表示为：
其中，其中Vi为母线i电压的幅值，θi为母线i电压的相角，Vj为母线j电压的幅值，θj为母线j电压的相角，I表示通过支路与母线i连接的母线j的集合；gij为连接母线i和母线j的支路的电导，bij为支路的电纳。PGi为母线i上的有功发电，PDi为母线i上的有功负荷，QGi为母线i的无功发电，QDi为母线i上的无功负荷。
从上述潮流方式可以进一步得到支路潮流计算表达式如下：
其中Pij为连接母线i及其相邻的母线j的支路上流动的有功潮流，Pij表示从母线i流向母线j的潮流，Pji表示从母线j流向母线j的潮流；支路包括连接母线的线路支路或变压器绕组支路。Qij为支路无功潮流，yc表示线路的半导纳。
潮流计算可以在电网结构参数确定的情况下确定电网的稳态运行状态。电网中的节点依给定的变量不同可划分为PQ节点、PV节点和Vθ节点，要计算的是电网状态变量即节点电压幅值和相角。潮流方程是一组高阶非线性代数方程组，要用迭代法求解。牛顿-拉夫逊潮流算法是具有二阶收敛性的潮流算法，因此得到了广泛的应用。但由于该方法的雅可比矩阵是待求状态变量函数，所以在迭代过程中要重新形成雅可比矩阵并进行高斯消去法求解，每次迭代的计算量较大。由于它是各种潮流计算方法的基础，因此在电网分析中有其重要的地位。
张伯明、陈寿孙等在《高等电力网络分析》(清华大学出版社，1996年，190-191页)中对牛顿-拉夫逊潮流算法进行了说明，该算法中的修正公式如下：
而式中雅可比矩阵的各个元素则分别为
为求取这些偏导数，可将Pi、Qi分别展开如下
计及
j≠i时，由于对特定的j，只有该特定节点的δj，从而特定的δij＝δi-δj是变量，可得
相似的，由于对特定的j，只有该特定节点的Uj是变量，可得
j＝i时，由于δi是变量，所有δij＝δi-δj都是变量，可得
相似的，由于Ui是变量，可得
根据式(5)-(8)，当Ui、Uj已知时，可以计算得到式(1)中雅各比矩阵的各元素，形成常数雅各比矩阵。
在AVC控制策略计算过程中，无功电压灵敏度计算目前有两种方法：一种是在B”矩阵中增广PV节点的传统灵敏度矩阵，一种是考虑无功控制设备响应的准稳态灵敏度矩阵。常规的灵敏度分析仅仅依赖于电力网络方程的线性化，不考虑电力设备(如：发电机、负荷等)对各种控制操作和扰动的准稳态的物理响应，这在一般的静态的电力网络分析中是可行的，但要服务于控制决策就无法实用了。基于上述考虑，在AVC的二级电压控制中，采用准稳态控制灵敏度分析方法，重点考虑发电机准稳态的无功电压外特性，在优化计算中考虑安装有AVC机组的准稳态响应，保障控制决策的精度和可靠性。
在AVC的变电站控制中，其控制对象为主变低压侧电容器、电抗器等无功设备，与发电机无功的连续调节不同，这些无功设备的投切具有离散特点，其投切后会造成变电站高、中、低三侧母线电压的阶跃性变化。为了保证AVC控制过程中各级母线电压的合格，要求无功设备对变电站高、中、低三侧母线电压的灵敏度非常精确，才能准确地对控制效果进行预估并生成正确的控制策略。
发明内容
本发明的目的是提出一种变电站内变压器分接头档位对母线电压影响的计算方法，通过计算给出变电站无功设备对母线电压影响的灵敏度结果，从而在变电站电压自动控制的过程中，对低压侧无功设备动作后，变电站内各母线的电压进行预估，保证变电站无功设备控制策略可靠性。
本发明提出的变电站内变压器分接头档位对母线电压影响的计算方法，包括以下步骤：
(1)设电网模型中的支路包括连接母线的线路或连接母线的变压器绕组，电网模型中支路无功潮流计算表达式如下：
其中，Pab表示电网模型中支路a-b首端a流向电网模型中支路a-b末端b的有功潮流，Qab表示从电网模型中支路a-b首端a流向电网模型中支路a-b末端b的无功潮流，Va为电网模型中支路a-b首端a电压的幅值，Vb为电网模型中支路a-b末端b电压的幅值，Gab为电网模型中支路a-b的电导，Bab为电网模型中支路a-b的电纳，yc表示电网模型中支路a-b的半导纳,θab为电网模型中支路a-b首端a电压相角与末端b电压相角的差值；
(2)由于电网模型中变电站内变压器绕组支路m-n的半导纳yc＝0，变电站内变压器绕组支路m-n的电导值Gmn远小于变电站内变压器绕组支路m-n的电纳值Bmn，因此上述支路无功潮流计算表达式中的Gmn忽略，当变压器绕组有载调压分接头变压比为t不为1时，根据上述步骤(1)的电网模型中支路无功潮流计算表达式，根据下式求出电网模型中变电站内变压器绕组支路m-n的首端m有功功率Pm和无功功率Qmn以及变压器绕组支路m-n的末端有功功率Pnm和无功功率Qnm对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量：
变电站内变压器绕组支路m-n首端m流向变电站内变压器绕组支路m-n末端n的有功功率Pmn和无功功率Qmn对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量为：
变电站内变压器绕组支路m-n末端n流向变电站内变压器绕组支路m-n首端m的有功功率Pnm和无功功率Qnm对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量为：
(3)根据上述步骤(2)中的变电站内变压器绕组的有载调压分接头变压比t发生变化时，对应变电站内变压器绕组支路m-n的首端m和末端n的有功潮流和无功潮流的变化量，构建一个变电站内变压器绕组支路m-n的潮流方程如下：
其中，为雅可比矩阵，H、M、N、L分别为雅可比矩阵的子阵，H表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化时，电网模型中有功功率对所有母线电压相角影响的偏导矩阵，L表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中有功功率对所有母线电压幅值影响的偏导矩阵，M表示在电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中无功功率对所有电压相角的影响的偏导矩阵，N表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中无功功率对电压幅值影响的偏导矩阵，Δθ'表示电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中所有母线电压相角变化量，ΔU′为电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中所有母线电压相角变化量，ΔP'表示网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，变电站内变压器绕组支路端点有功功率变化量，ΔQ'表示网模型中变电站内变压器绕组的有载调压分接头变压比发生变化时，变电站内变压器绕组支路端点无功功率变化量；
(4)利用上述步骤(3)的变电站内变压器绕组支路m-n的潮流方程，计算得到电网模型中变电站内变压器绕组支路m-n的有载调压分接头变压比t发生一个档位变化时，变电站内变压器绕组支路m-n首端m和末端n的有功功率变化量ΔP'和无功功率变化量ΔQ'，有功功率变化量ΔP'和无功功率变化量ΔQ'的数值来自上述步骤(2)，其它元素为0，即首先将步骤(3)的潮流方程的右边向量置为:
再求解上述步骤(3)的潮流方程，得到变电站内变压器绕组支路m-n的有载调压分接头变压比t发生一个档位变化后，电网模型中所有母线电压幅值变化量ΔU′和电压相角变化量Δθ'；
(5)设电网模型中第x条母线与变电站内变压器绕组支路m-n同属于一个变电站，则从ΔU′中得到第x条母线对应的电压幅值变化量ΔUx′，电网模型中变电站内变压器绕组支路m-n有载调压分接头档位t发生变化时，对变电站内第x条母线的电压幅值影响灵敏度为：
Stx＝ΔUx′Δtstep
其中，Δtstep为变电站内变压器绕组支路m-n的有载调压分接头单位增加1档，对应的变压比t的值，t为常数，取值范围为1－100；
(6)遍历电网模型中所有变电站内变压器绕组支路的有载调压分接头，返回步骤(4)，完成变电站内变压器分接头档位对母线电压影响的计算。
本发明提出的变电站内变压器分接头档位对母线电压影响的计算方法，在实时的潮流计算结果基础上，首先求出采用准稳态潮流构造注入潮流变化量与母线电压变化量的线性方程，然后计算出在某一条母线上注入无功后对电网模型中母线电压的变化幅值，即该母线对电网模型中其它母线的无功电压灵敏度，由此可以得出变电站无功设备所连接的母线对电网模型中其它母线的电压灵敏度，即变电站无功设备对母线无功电压灵敏度。通过本发明计算方法，可以给出变电站无功设备对母线电压影响的灵敏度即无功电压灵敏度，可以提高变电站自动电压控制的可靠性。
附图说明
图1是本发明计算方法涉及的电网模型中，变电站内变压器绕组支路m-n模型。
图2是本发明方法涉及的电网模型中的变电站模型。
具体实施方式
本发明提出的一种变电站内变压器分接头档位对母线电压影响灵敏度的计算方法，结合实施例详细说明如下：
该方法包括：
(1)设电网模型中的支路包括连接母线的线路或连接母线的变压器绕组，电网模型中支路无功潮流计算表达式如下：
其中，Pab表示电网模型中支路a-b首端a流向电网模型中支路a-b末端b的有功潮流，Qab表示从电网模型中支路a-b首端a流向电网模型中支路a-b末端b的无功潮流，Va为电网模型中支路a-b首端a电压的幅值，Vb为电网模型中支路a-b末端b电压的幅值，Gab为电网模型中支路a-b的电导，Bab为电网模型中支路a-b的电纳，yc表示电网模型中支路a-b的半导纳,θab为电网模型中支路a-b首端a电压相角与末端b电压相角的差值；本发明实施例考虑一个有高、中、低三侧母线的变电站，具体连接关系见附图2，变电站内一条变压器绕组支路为m-n具体连接关系见附图1。以下通过具体的例子来分析灵敏度计算的过程。
(2)由于电网模型中变电站内变压器绕组支路m-n的半导纳yc＝0，变电站内变压器绕组支路m-n的电导值Gmn远小于变电站内变压器绕组支路m-n的电纳值Bmn，上述支路无功潮流计算表达式中的Gmn忽略，当变压器绕组有载调压分接头变压比为t不为1时，根据上述步骤(1)的电网模型中支路无功潮流计算表达式，根据下式求出电网模型中变电站内变压器绕组支路m-n的首端m有功功率Pm和无功功率Qmn以及变压器绕组支路m-n的末端有功功率Pnm和无功功率Qnm对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量：
变电站内变压器绕组支路m-n首端m流向变电站内变压器绕组支路m-n末端n的有功功率Pmn和无功功率Qmn对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量为：
变电站内变压器绕组支路m-n末端n流向变电站内变压器绕组支路m-n首端m的有功功率Pnm和无功功率Qnm，对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量为：
(3)根据上述步骤(2)中的变电站内变压器绕组的有载调压分接头变压比t发生变化时，对应变电站内变压器绕组支路m-n的首端m和末端n的有功潮流和无功潮流的变化量，构建一个变电站内变压器绕组支路m-n的潮流方程如下：
其中，为雅可比矩阵，其中H、M、N、L分别为雅可比矩阵的子阵，H表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化时，电网模型中有功功率对所有母线电压相角影响的偏导矩阵，L表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中有功功率对所有母线电压幅值影响的偏导矩阵，M表示在电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中无功功率对所有电压相角的影响的偏导矩阵，N表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中无功功率对电压幅值影响的偏导矩阵，Δθ'表示电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中所有母线电压相角变化量，ΔU′为电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中所有母线电压相角变化量，ΔP'表示网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，变电站内变压器绕组支路端点有功功率变化量，ΔQ'表示网模型中变电站内变压器绕组的有载调压分接头变压比发生变化时，变电站内变压器绕组支路端点无功功率变化量；
(4)利用上述步骤(3)的变电站内变压器绕组支路m-n的潮流方程，计算得到电网模型中变电站内变压器绕组支路m-n的有载调压分接头变压比t发生一个档位变化时，变电站内变压器绕组支路m-n首端m和末端n的有功功率变化量ΔP'和无功功率变化量ΔQ'，有功功率变化量ΔP'和无功功率变化量ΔQ'的数值来自上述步骤(2)，其它元素为0，即首先将步骤(3)的潮流方程的右边向量置为:
求解上述步骤(3)的潮流方程，得到变电站内变压器绕组支路m-n的有载调压分接头变压比t发生一个档位变化后，电网模型中所有母线电压幅值变化量ΔU′和电压相角变化量Δθ'；
(5)设电网模型中第x条母线与变电站内变压器绕组支路m-n同属于一个变电站，则从ΔU′中得到第x条母线对应的电压幅值变化量ΔUx′，电网模型中变电站内变压器绕组支路m-n有载调压分接头档位t发生变化时，对变电站内第x条母线的电压幅值影响灵敏度为：
Stx＝ΔUx′Δtstep
其中，Δtstep为变电站内变压器绕组支路m-n的有载调压分接头单位增加1档，对应的变压比t的值Δtstep＝1，则电网模型中变电站内变压器绕组支路m-n对变电站内第x条母线的电压幅值影响灵敏度为0.447kv/档。
(6)遍历电网模型中所有变电站内变压器绕组支路的有载调压分接头，返回步骤(4)，完成变电站内变压器分接头档位对母线电压影响的计算。
Claims (1)
1.一种变电站内变压器分接头档位对母线电压影响的计算方法，该方法包括以下步骤：
(1)设电网模型中的支路包括连接母线的线路或连接母线的变压器绕组，电网模型中支路无功潮流计算表达式如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>V</mi>
<mi>a</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>G</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>a</mi>
</msub>
<msub>
<mi>V</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mi>a</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>a</mi>
</msub>
<msub>
<mi>V</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，Pab表示电网模型中支路a-b首端a流向电网模型中支路a-b末端b的有功潮流，Qab表示从电网模型中支路a-b首端a流向电网模型中支路a-b末端b的无功潮流，Va为电网模型中支路a-b首端a电压的幅值，Vb为电网模型中支路a-b末端b电压的幅值，Gab为电网模型中支路a-b的电导，Bab为电网模型中支路a-b的电纳，yc表示电网模型中支路a-b的半导纳,θab为电网模型中支路a-b首端a电压相角与末端b电压相角的差值；
(2)由于电网模型中变电站内变压器绕组支路m-n的半导纳yc＝0，变电站内变压器绕组支路m-n的电导值Gmn远小于变电站内变压器绕组支路m-n的电纳值Bmn，因此在支路无功潮流计算表达式中Gmn可以忽略不计，当变压器绕组有载调压分接头变压比为t不为1时，根据上述步骤(1)的电网模型中支路无功潮流计算表达式，根据下式求出电网模型中变电站内变压器绕组支路m-n的首端m有功功率Pmn和无功功率Qmn以及变压器绕组支路m-n的末端有功功率Pnm和无功功率Qnm对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量：
变电站内变压器绕组支路m-n首端m流向变电站内变压器绕组支路m-n末端n的有功功率Pmn和无功功率Qmn对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量为：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>m</mi>
</msub>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>V</mi>
<mi>m</mi>
</msub>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
变电站内变压器绕组支路m-n末端n流向变电站内变压器绕组支路m-n首端m的有功功率Pnm和无功功率Qnm对变压器绕组有载调压分接头变压比t的有功潮流、无功潮流变化量为：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>V</mi>
<mi>m</mi>
</msub>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
<msub>
<mi>B</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<msup>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>B</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>m</mi>
</msub>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
<msub>
<mi>B</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
(3)根据上述步骤(2)中的变电站内变压器绕组的有载调压分接头变压比t发生变化时，对应变电站内变压器绕组支路m-n的首端m和末端n的有功潮流和无功潮流的变化量，构建一个变电站内变压器绕组支路m-n的潮流方程如下：
<mrow>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>H</mi>
</mtd>
<mtd>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>J</mi>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;&theta;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;U</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;P</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;Q</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，为雅可比矩阵，H、M、N、L分别为雅可比矩阵的子阵，H表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化时，电网模型中有功功率对所有母线电压相角影响的偏导矩阵，L表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中有功功率对所有母线电压幅值影响的偏导矩阵，M表示在电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中无功功率对所有电压相角的影响的偏导矩阵，N表示电网当前潮流状态下，电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中无功功率对电压幅值影响的偏导矩阵，Δθ'表示电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中所有母线电压相角变化量，ΔU′为电网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，电网模型中所有母线电压相角变化量，ΔP'表示网模型中变电站内变压器绕组的有载调压分接头变压比发生变化后，变电站内变压器绕组支路端点有功功率变化量，ΔQ'表示网模型中变电站内变压器绕组的有载调压分接头变压比发生变化时，变电站内变压器绕组支路端点无功功率变化量；
(4)利用上述步骤(3)的变电站内变压器绕组支路m-n的潮流方程，计算得到电网模型中变电站内变压器绕组支路m-n的有载调压分接头变压比t发生一个档位变化时，变电站内变压器绕组支路m-n首端m和末端n的有功功率变化量ΔP'和无功功率变化量ΔQ'，有功功率变化量ΔP'和无功功率变化量ΔQ'的数值来自上述步骤(2)，其它元素为0，即首先将步骤(3)的潮流方程的右边向量置为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;P</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>&Delta;Q</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
再求解上述步骤(3)的潮流方程，得到变电站内变压器绕组支路m-n的有载调压分接头变压比t发生一个档位变化后，电网模型中所有母线电压幅值变化量ΔU′和电压相角变化量Δθ'；
(5)设电网模型中第x条母线与变电站内变压器绕组支路m-n同属于一个变电站，则从ΔU′中得到第x条母线对应的电压幅值变化量ΔU′x，电网模型中变电站内变压器绕组支路m-n有载调压分接头档位t发生变化时，对变电站内第x条母线的电压幅值影响灵敏度为：
Stx＝ΔU′xΔtstep
其中，Δtstep为变电站内变压器绕组支路m-n的有载调压分接头单位增加1档，对应的变压比t的值，t为常数，取值范围为1－100；
(6)遍历电网模型中所有变电站内变压器绕组支路的有载调压分接头，返回步骤(4)，完成变电站内变压器分接头档位对母线电压影响的计算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610692725.0A CN106099940B (zh) | 2016-08-19 | 2016-08-19 | 变电站内变压器分接头档位对母线电压影响的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610692725.0A CN106099940B (zh) | 2016-08-19 | 2016-08-19 | 变电站内变压器分接头档位对母线电压影响的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106099940A CN106099940A (zh) | 2016-11-09 |
| CN106099940B true CN106099940B (zh) | 2018-04-27 |
Family
ID=58069741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610692725.0A Active CN106099940B (zh) | 2016-08-19 | 2016-08-19 | 变电站内变压器分接头档位对母线电压影响的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106099940B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110134709B (zh) * | 2019-04-12 | 2023-06-09 | 云南电网有限责任公司曲靖供电局 | 电网数据查询方法及装置 |
| CN113690893B (zh) * | 2021-08-26 | 2022-08-02 | 广东电网有限责任公司广州供电局 | 一种含调压器微电网的潮流计算方法和系统 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN203690924U (zh) * | 2014-02-08 | 2014-07-02 | 广西电网公司河池供电局 | 一种可控电抗器协调控制装置 |
| CN104917186B (zh) * | 2015-06-16 | 2017-12-08 | 傅辉明 | 基于电压无功解耦判据实现最优目标控制的解耦协调方法 |
| CN105140925A (zh) * | 2015-07-21 | 2015-12-09 | 河南行知专利服务有限公司 | 一种电力系统自动电压控制方法及装置 |
-
2016
- 2016-08-19 CN CN201610692725.0A patent/CN106099940B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106099940A (zh) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103269069B (zh) | 一种低压输电系统网络重构的方法 | |
| US20060282239A1 (en) | Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system | |
| CN110718918B (zh) | 一种基于hvac和lcc-vsc hvdc混合系统的潮流计算方法 | |
| CN105552920B (zh) | 一种基于馈线负荷分支电压灵敏度的含分布式电源的10kV馈线电压调整控制方法 | |
| CN106099939A (zh) | 一种变电站无功设备对母线电压影响灵敏度的计算方法 | |
| CN105375513B (zh) | 一种基于实时在线等值的110千伏风电场自动电压控制方法 | |
| Seifi et al. | Power flow study and comparison of FACTS: Series (SSSC), Shunt (STATCOM), and Shunt-Series (UPFC) | |
| Mohamed et al. | Line-wise power flow and voltage collapse | |
| CN105162144A (zh) | 一种配电网电压无功优化的线性逼近求解方法 | |
| CN104934973B (zh) | 一种考虑温度的电力线路损耗计算方法 | |
| CN103346573B (zh) | 基于黄金分割云粒子群优化算法的风电系统无功规划方法 | |
| CN105048473B (zh) | 变电站自动电压控制中低压无功设备电压灵敏度统计方法 | |
| CN106953320A (zh) | 一种基于upfc的输电线路过载控制方法 | |
| CN106099940B (zh) | 变电站内变压器分接头档位对母线电压影响的计算方法 | |
| CN102904265B (zh) | 一种基于区间潮流的变电站电压无功控制方法和系统 | |
| CN104700205B (zh) | 一种改变电网网络拓扑结构及选择并联补偿装置的方法 | |
| CN106229995A (zh) | 基于风电场抗台风运行模式下的备用电源并联电抗器参数优化方法 | |
| CN104300546B (zh) | 一种基于广域量测信息的电压稳定约束无功优化方法 | |
| Ram et al. | Voltage stability analysis using L-index under various transformer tap changer settings | |
| CN107451743B (zh) | 一种考虑源荷运行特性的有源配电网孤岛运行方法 | |
| CN105071397A (zh) | 风电外送的多类型无功补偿设备的协调无功电压控制方法 | |
| CN106159946B (zh) | 变电站内无功设备对变压器绕组无功功率影响的计算方法 | |
| Rao et al. | Enhancement of power quality by using unified power quality conditioner with pid and fuzzy logic controller | |
| CN107769187A (zh) | 一种适用于电压型直流潮流控制器的优化配置方法 | |
| CN109830964B (zh) | 一种微电网动态功率平衡方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191129 Address after: Yunda economic and Technological Development Zone in Yunnan province Kunming city 650217 West Road No. 105 Patentee after: Electric Power Research Institute of Yunnan Power System Ltd Address before: Guandu District of Yunnan city of Kunming Province Yun Road 650000 No. 105 cloud electric science and Technology Park Co-patentee before: Beijing King Star Hi-Tech System Control Co., Ltd. Patentee before: Electric Power Research Institute of Yunnan Power System Ltd |
|
| TR01 | Transfer of patent right |