CN109980648B - 交直流混联潮流计算的方法及其装置、存储介质及终端 - Google Patents
交直流混联潮流计算的方法及其装置、存储介质及终端 Download PDFInfo
- Publication number
- CN109980648B CN109980648B CN201910309632.9A CN201910309632A CN109980648B CN 109980648 B CN109980648 B CN 109980648B CN 201910309632 A CN201910309632 A CN 201910309632A CN 109980648 B CN109980648 B CN 109980648B
- Authority
- CN
- China
- Prior art keywords
- converter
- voltage
- direct current
- control mode
- active power
- 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
- 238000000034 method Methods 0.000 title claims abstract description 106
- 238000004364 calculation method Methods 0.000 claims abstract description 132
- 239000011159 matrix material Substances 0.000 claims abstract description 96
- 238000002940 Newton-Raphson method Methods 0.000 claims abstract description 42
- 238000002347 injection Methods 0.000 claims description 95
- 239000007924 injection Substances 0.000 claims description 95
- 230000009466 transformation Effects 0.000 claims description 43
- 230000008569 process Effects 0.000 claims description 24
- 230000008859 change Effects 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000010791 quenching Methods 0.000 claims description 2
- 238000010304 firing Methods 0.000 claims 1
- 238000012804 iterative process Methods 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 abstract description 4
- 239000002609 medium Substances 0.000 description 14
- 239000000243 solution Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000009365 direct transmission Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000012120 mounting media Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000171 quenching effect Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/133—Arrangements for measuring electric power or power factor by using digital technique
- G01R21/1331—Measuring real or reactive component, measuring apparent energy
-
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/001—Measuring real or reactive component; Measuring apparent energy
- G01R21/003—Measuring reactive component
-
- 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- 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
- H02J5/00—Circuit arrangements for transfer of electric power between ac networks and dc networks
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明实施例公开了一种交直流混联潮流计算的方法及其装置、存储介质及终端，根据换流站的控制方式完成直流网络的状态计算，将直流网络与交流网络的连接点等效为功率节点；使用牛顿‑拉夫逊法计算潮流。本发明实施例既克服了交替迭代法的交替迭代引起的收敛性变差问题，又避免了联立求解法所引起的初值选取、雅可比矩阵规模扩大问题；具有好的收敛性和迭代过程内存占用少双重优点。
Description
本发明要求在2018年09月07日提交中国专利局、申请号为201811045461.5的中国专利申请的优先权，该申请的全部内容通过引用结合在本发明中。
技术领域
本发明实施例涉及电力系统潮流计算技术领域，尤其涉及一种交直流混联潮流计算的方法及其装置、存储介质及终端。
背景技术
潮流计算是研究电力系统稳态运行情况的一种基本计算。从电的产生到被负荷消耗的过程中流经输配电线中各节点的电压和/或功率等，可通过潮流计算的方式获得。当前，交直流混联电力系统潮流的计算方法主要有交替迭代法和联立求解法两大类。
交替迭代法的主要优点是主迭代中未改变原节点导纳矩阵和雅可比矩阵，只需要稍微修改节点功率平衡方程，易于与原有的潮流算法相结合而编程实现；而交替迭代法的缺点是新增元件装置的控制变量只在子迭代中被修正，而在主迭代过程中控制变量值保持子迭代中修正后的给定值不变，两部分迭代过程相交互所造成的差异使得全算法收敛特性变差，甚至发生数值振荡或者发散，导致算法不收敛，不再具有传统牛顿—拉夫逊法的二阶收敛特性。
联立求解法的优点是保留了传统潮流算法的收敛特性。联立求解法能够将求解系统运行状态变量的方程组与求解新增元件控制变量的方程组进行统一的联立迭代求解，具有传统牛顿—拉夫逊法的收敛特性；并且，联立求解法与原电网潮流计算相比，增加了新的状态变量和控制目标方程或内部约束方程，需要对原雅可比矩阵进行修改扩充。新增的控制变量要考虑初始值的选取，而牛顿—拉夫逊法的求解又对变量的初始值依赖较强，因此联立求解法也存在收敛速度减慢，收敛可靠性变差的问题。同时，新增控制目标方程与经典潮流方程的表达式差异较大，可能出现修正方程病态的情况。
发明内容
针对上述存在问题，本发明实施例提供一种交直流混联潮流计算的方法及其装置、存储介质及终端，能够解决现有技术中潮流计算的收敛可靠性差、易出现修正方程病态的技术问题。
第一方面，本发明实施例提供了一种交直流混联潮流计算的方法，包括：
对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；
根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率；
获取各所述换流器的控制方式；
根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量；
根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量和所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算。
可选的，所述根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率，包括：
根据所述直流网络的结构，获取与所述换流器对应的节点参数；
根据所述节点参数，构建方程组，如下：
其中，Vdk为与所述换流器对应的节点的直流电压，Pdk为与所述换流器对应的节点的有功功率，Idk为流入换流站k的直流电流，Gkj为对应节点k和j之间的导纳矩阵元素，Vdj为与换流器j连接的直流母线的电压，nc为所述直流网络中换流器的个数；
根据所述方程组，求解与所述换流器对应的节点的直流电压和有功功率。
可选的，在所述交直流混联电力系统的直流网络包括分层结构时，所述直流网络串联侧单个换流器输出的有功功率与所述换流器承受的电压占比成正比。
可选的，在所述交直流混联电力系统的直流网络包括分层结构时，满足如下关系：
其中，Idi1和Idi2分别为流经分层结构下换流站的高电压换流器的电流和低电压换流器的电流；Id表示流经换流站整体的电流；Vdr为直流网络送端电压；Vdi1和Vdi2分别表示分层结构下高、低电压换流器的直流电压；Rd为直流线路的电阻；
所述根据所述方程组，求解与所述换流器对应的节点的直流电压和有功功率，具体为：
其中，kidk为换流器k在分层结构中的电压占比，Pd为直流注入换流站的有功功率，Vd为换流站所连接节点的直流电压，Pidk为分层结构下换流器k所输出的有功功率，Vidk为分层结构下换流器k所承受的直流电压。
可选的，所述换流器的控制方式包括第一类控制方式和第二类控制方式；
所述第一类控制方式包括恒定有功功率、恒定直流电压和恒定直流电流；
所述第二类控制方式包括恒定变压器变比和恒定换相角。
可选的，所述换流器的控制方式为所述第一类控制方式时，所述根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量，具体为：
其中，Idk为流入换流器k的直流电流，Pdk为所述有功功率，Vdk为所述直流电压，
为换流器功率因素，Qdk为所述无功注入量。
可选的，所述换流器的控制方式为第二类控制方式，且为恒定换相角的控制方式时，所述根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量，具体为：
其中，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，Pidk为直流注入交流节点i的有功功率，θd为换流器的控制角，Xc为换相电阻，ky为换流器常数，Qdk为所述无功注入量。
可选的，所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，所述根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量，具体为：
其中，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，Va为与换流器相连节点的电压幅值，kT为变压器变比，ky为换流器常数。
可选的，所述根据任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算，包括：
根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量；
根据所述换流器的控制方式以及所述有功功率的不平衡量和所述无功注入量的不平衡量，建立所述潮流计算的雅可比矩阵；其中，在所述换流器的控制方式为第二类控制方式，且为恒定换相角的控制方式时，与所述换流器对应的节点的雅可比矩阵参数仅由交流网络参数求出；所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，与所述换流器对应的节点的雅可比矩阵参数由交流网络参数求出后进行修正；
根据所述雅可比矩阵，通过牛顿-拉夫逊法进行潮流计算。
可选的，在所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，所述根据所述雅可比矩阵，通过牛顿-拉夫逊法进行潮流计算，还包括：
对所述雅可比矩阵元素Lii进行修正，如下：
其中，i为与所述换流器连接的交流网路的节点；Vi为对应节点i的电压幅值，Gij和Bij为导纳矩阵的实部与虚部，Va为与所述换流器相连节点的电压幅值，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，kT为变压器变比，ky为换流器常数，θij为节点i的控制角，H、N、L为所述雅可比矩阵的分块矩阵，ΔP为所述有功功率不平衡量，ΔQ为所述无功注入量的不平衡量，Δθ和ΔV为迭代过程中的变量的修正量。
可选的，所述通过牛顿-拉夫逊法进行潮流计算，还包括：
判断所述潮流计算量是否满足收敛条件；
若是则完成潮流计算；
若否，则重新获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量。
第二方面，本发明实施例还提供了一种交直流混联潮流计算的装置，包括：
电阻获取模块，用于对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；
直流电压和有功功率获取模块，用于根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率；
控制方式获取模块，用于获取各所述换流器的控制方式；
无功注入量计算模块，用于根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量；
潮流计算模块，用于根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流母线的电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算。
第三方面，本发明实施例还提供了一种存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述交直流混联潮流计算的方法。
第四方面，本发明实施例还提供了一种终端，显示屏、存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述交直流混联潮流计算的方法。
本发明实施例提供的一种交直流混联潮流计算的方法及其装置、存储介质和终端，通过分析交直流混联电力系统中各换流器的控制方式，获得相应的参数，以通过牛顿-拉夫逊法进行潮流计算，能够避免因初值选取规模大、雅可比矩阵规模大，不利于计算的问题，从而使得交直流混联潮流计算的具有较好的收敛性，降低计算的复杂性，进一步提高计算速率、降低成本。
另外，本发明为了解决上述问题，还提出了一种计及直流分层结构交直流混联潮流计算的方法及系统，本发明通过将直流网络与交流网络的连接点等效为交流网络中的功率节点，使得潮流计算具有更好的收敛速度、收敛可靠性。
为了实现上述目的，本发明采用如下技术方案：
一种计及直流分层结构交直流混联潮流计算的方法，根据换流站的控制方式完成直流网络的状态计算，将直流网络与交流网络的连接点等效为功率节点；使用牛顿—拉夫逊法计算潮流。
具体的，包括以下步骤：
对直流网络求解电导矩阵，得到每两个换流器间的电阻，或者与分层结构的直流电网连接点间的电阻值；
分析各换流器的控制方式，得到每个节点的电压和有功功率；
根据所得到的控制方式和各节点的电压和有功功率，计算对交流电网的无功注入量；
使用牛顿-拉夫逊法完成潮流计算，得到计算结果。
进一步的，依据换流器控制方式确定节点部分参数，继而构造方程组，以求解节点电压和有功功率为目标完成计算，具体方程组为：
其中：Idi1，Idi2分别表示流经分层结构下换流站的高、低电压换流器的电流；Id表示流经换流站整体的电流；Vdr表示直流网络送端电压；Vdi1、Vdi2分别表示分层结构下高、低电压换流器的直流电压；Rd表示直流线路的电阻。
进一步的，串联侧单个换流器输出的有功功率与所述换流器承受的电压占比成正比。
进一步的，依据节点换流器的控制方式，即恒定换相角控制控制方式下使用式(2)，恒定变压器变比控制方式下使用式(3)完成对交流电网无功注入量的计算，有功注入量的计算通过直流网络状态计算得到；
其中：Vdc为换流站所连接的直流网络节点电压；θd为换流器的控制角，即整流器的触发延迟角、逆变器的熄弧超前角；kT为变压器变比；Xc为换相电阻；考虑到换相角的影响，引入了变量ky；
为换流器从交流系统吸收有功功率(整流器为吸收，逆变器为发出)和无功功率对应的功率因素角；Va为换流器连接的交流网络的电压幅值。
进一步的，涉及到分层结构时，根据分层的换流器电压占比完成功率节点的等效。
进一步的，计算对交流电网的无功注入量的具体过程包括：
(1)如果一个节点对应的换流器控制方式为恒定换相角，则使用公式(2)
计算无功注入量，接着转至步骤(3)；否则转至步骤(2)；
(2)如果一个节点对应的换流器控制方式为恒定变压器变比，则使用公式(3)计算无功注入量，并计算无功注入量对所对应交流电压的导数；
(3)如果存在分层结构，依据
计算各层对交流电网连接点的功率影响；
kidk为换流器k在分层结构中的电压占比，Pd为直流注入换流站的有功功率，Vd为换流站所连接节点的直流电压，Pidk为分层结构下换流器k所输出的有功功率，Vidk为分层结构下换流器k所承受的直流电压。
使用牛顿-拉夫逊法完成潮流计算的具体过程包括：
(a)设定交流网络初始值，求解潮流功率方程不平衡量；
(b)构造雅可比矩阵，流器控制方式为恒定换相角对应的节点所对应的雅可比矩阵参数仅由交流网络参数求出；换流器控制方式为恒定变压器变比对应的节点所对应的雅可比矩阵参数由交流网络参数求出后进行修正；
(c)完成交流网络参数修正，检查收敛条件，达到条件结束迭代，否则转步骤(a)。
更进一步的，修正的方式为：
其中：-Vi∑j∈i，j≠iVj(Gij sinθij-Bij cosθij)+2Vi 2Bii为传统纯交流潮流计算中雅可比矩阵元素L的计算式，Vi为对应节点i的电压幅值，Gij和Bij为导纳矩阵的实部与虚部，Va换流器相连节点的电压幅值，数值上与Vi一致。
一种计及直流分层结构交直流混联潮流计算的系统，运行于处理器上，被配置为执行以下指令：
对直流网络求解电导矩阵，得到每两个换流器间的电阻，或者与分层结构的直流电网连接点间的电阻值；
分析各换流器的控制方式，得到每个节点的电压和有功功率；
根据所得到的控制方式和各节点的电压和有功功率，计算对交流电网的无功注入量；
使用牛顿-拉夫逊法完成潮流计算，得到计算结果。
与现有技术相比，本发明的有益效果为：
1.本发明既克服了交替迭代法的交替迭代引起的收敛性变差问题，又避免了联立求解法所引起的初值选取、雅可比矩阵规模扩大问题；具有好的收敛性和迭代过程内存占用少双重优点。
2.本发明对于既有的纯交流潮流计算程序改动量小，节约软件更新成本。
3.本发明的技术思路完全适用于当前电网新型器件带来的新的网络组成的潮流计算，易于相关软件形成标准化处理。
附图说明
图1是本发明实施例提供的一种交直流混联潮流计算的方法的流程图；
图2是本发明实施例提供的一种分层结构的直流网络的电路结构示意图；
图3是本发明实施例提供的一种获取直流电压和有功功率的方法的流程图；
图4是本发明实施例体用的一种通过牛顿-拉夫逊法进行潮流计算方法的流程图；
图5是本发明实施例提供的又一种通过牛顿-拉夫逊法进行潮流计算方法的流程图；
图6是本发明实施例提供的一种交直流混联潮流计算的装置的结构框图；
图7为本发明实施例提供的一种终端的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
本发明实施例提供了一种交直流混联潮流计算的方法，该方法可适用于对交直流混联电力系统进行潮流计算的情况。本发明实施例提供的交直流混联潮流计算的方法可以由交直流混联潮流计算的装置来执行，该装置可由软件和/或硬件来实现。图1是本发明实施例提供的一种交直流混联潮流计算的方法的流程图。如图1，一种交直流混联潮流计算的方法包括：
S110、对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻。
具体的，交直流混联电力系统包括直流网络和交流电网，直流网络与交流网络通过换流器相互连接，换流器能够将交流电网中的交流电信号转换为直流电信号后输入直流网络，换流器还能够使直流网络中的直流电信号转换为交流电信号后输入交流网络，因此交直流混联电力系统的换流器对交直流混联电力系统的稳定运行具有重要贡献。对于多端的直流网络，可通过求解该直流网络的电导矩阵，即可获得该直流网络中任意两个换流器之间的电阻。或者在交直流混联电力系统的直流网络存在分层结构时，还可通过求解该直流网络的电导矩阵，获得各分层结构的所述直流网络的连接点间的电阻。
S120、根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率。
具体的，直流网络参数通常包括电容、电感等参数值，但是在对交直流混联电力系统的稳定运行进行研究和分析时，通常只考虑直流网络的电阻特性，由直流网络的节点的导纳矩阵Gd表示直流网络：
直流网络的节点的注入电流Id可表示为：
Id＝GdVd
式中Id为直流网络的节点注入的直流电流，Vd为直流网络的节点注入的直流电压。相应的，与换流器对应的节点的直流电压Idk和有功功率Pdk可通过下述方程组进行计算：
其中，Vdk为与所述换流器对应的节点的直流电压，Pdk为与所述换流器对应的节点的有功功率，Idk为流入换流站k的直流电流，Gkj为对应节点k和j之间的导纳矩阵元素，Vdj为与换流器j连接的直流母线的电压，nc为所述直流网络中换流器的个数。
此外，在交直流混联电力系统的直流网络包括分层结构时，直流网络串联侧单个换流器输出的有功功率与所述换流器承受的电压占比成正比。图2是本发明实施例提供的一种分层结构的直流网络的电路结构示意图。如图2，在交直流混联电力系统的直流网络包括分层结构时，直流网络的节点的注入电流Id还应满足如下关系：
其中，Idi1和Idi2分别为流经分层结构下换流站的高电压换流器的电流和低电压换流器的电流；Id表示流经换流站整体的电流；Vdr为直流网络送端电压；Vdi1和Vdi2分别表示分层结构下高、低电压换流器的直流电压；Rd为直流线路的电阻。相应的，与换流器i对应的节点的直流电压Vidk和有功功率Pidk，具体为：
其中，kidk为换流器k在分层结构中的电压占比，Pd为直流注入换流站的有功功率，Vd为换流站所连接节点的直流电压，Pidk为分层结构下换流器k所输出的有功功率，Vidk为分层结构下换流器k所承受的直流电压。
S130、获取各所述换流器的控制方式。
具体的，对于传统的换相换流器，每个换流器有两个独立的控制变量。假定与换流器i相关的变压器抽头可以无缝调节，那么变压器的匝数比kti可以被线性控制。因此，与换流器i连接的直流母线的有功功率Pdci、直流电压Vdci和直流电流Idci可被认为是第一类控制方式中的控制变量，该第一类控制方式可定义为D轴控制方式；而与换流器i相关的变压器的变压器变比(匝数比)kti和换流器i的控制角θi可被认为是第二类控制方式中的控制变量，该第二类控制方式可定义为E轴控制方式。相应的，换流器的控制方式可以分为第一类控制方式和第二类控制方式。
S140、根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量。
具体的，交直流混联电力系统的各换流器的控制方式不同，其对交流电网的无功注入量计算方式不同。换流器的控制方式可以分为第一类控制方式和第二类控制方式。其中，第一类控制方式包括恒定有功功率、恒定直流电压和恒定直流电流等，第二类控制方式包括恒定变压器变比和恒定换相角等，这里的换相角为换流器的控制角的一个确定值。
当换流器的控制方式为所述第一类控制方式时，根据直流电压、有功功率以及换流器的控制方式，计算换流器对交流电网的无功注入量，具体为：
其中，Idk为流入换流器k的直流电流，Pdk为有功功率，Vdk为直流电压，
为换流器功率因素，Qdk为无功注入量。
当换流器的控制方式为第二类控制方式时，需将上述计算无功注入量的方程组与换流器的基本方程结合计算对应的无功注入量。换流器的基本方程如下：
其中，其中，
为直流输电电压的标幺值，
为直流输电电流的标幺值，
为交流母线的线电压的标幺值，
是注入换流器的基频交流电流，kT为变压器变比，θd为换流器的控制角，即整流器的触发延迟角、逆变器的熄弧超前角，
为为换流器从交流系统吸收有功功率(整流器为吸收，逆变器为发出)和无功功率对应的功率因素角，
为换相电阻的标幺值，ky为换流器常数，考虑到换相角的影响，简化分析，近似取常数0.995。
当换流器的控制方式为第二类控制方式，且为恒定换相角的控制方式时，根据直流电压、有功功率以及换流器的控制方式，计算与换流器对交流电网的无功注入量，依据上式可以得出：
其中，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，Pidk为直流注入交流节点i的有功功率，θd为换流器的控制角，Xc为换相电阻，ky为换流器常数，Qdk为所述无功注入量。
当换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，根据直流电压、有功功率以及换流器的控制方式，计算与换流器对交流电网的无功注入量，依据上式可以得出：
其中，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，Va为与换流器相连节点的电压幅值，kT为变压器变比，ky为换流器常数。
S150、根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量和所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算。
具体的，现有的牛顿-拉夫逊法进行潮流计算，通常选取相应的初始值进行迭代计算，此时潮流计算的结果与初始值的选取有关，会因初始值选取规模大，而增加迭代次数，此时雅可比矩阵规模大，不利于计算。通过获取任意两个换流器之间的电阻或连接点间的电阻，再由换流器的控制方式计算直流电压、有功功率、所述无功注入量，通过牛顿-拉夫逊法完成交直流混联电力系统的潮流计算。
本发明实施例通过分析交直流混联电力系统中各换流器的控制方式，获得相应的参数，以通过牛顿-拉夫逊法进行潮流计算，能够避免因初值选取规模大、雅可比矩阵规模大，不利于计算的问题，从而使得交直流混联潮流计算的具有较好的收敛性，降低计算的复杂性，进一步提高计算速率、降低成本。
可选的，在上述实施例的基础上，对获取直流电压和有功功率的具体方法进行了优化。图3是本发明实施例提供的一种获取直流电压和有功功率的方法的流程图。如图3，根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率，具体包括：
S131、根据所述直流网络的结构，获取与所述换流器对应的节点参数；
S132、根据所述节点参数，构建方程组，如下：
其中，Vdk为与所述换流器对应的节点的直流电压，Pdk为与所述换流器对应的节点的有功功率，Idk为流入换流站k的直流电流，Gkj为对应节点k和j之间的导纳矩阵元素，Vdj为与换流器j连接的直流母线的电压，nc为所述直流网络中换流器的个数；
S133、根据所述方程组，求解与所述换流器对应的节点的直流电压和有功功率。
具体的，交直流混联电力系统的直流网络的结构不同，其换流器的对应的节点的直流电压和有功功率计算方式不同。对于一般的直流网络的换流器计算对应节点的直流电压和有功功率可由节点参数构建方程组：
从而根据上述方程组进行直流电压和有功功率的计算。
而对于分层结构的直流网络，其串联侧单个换流器输出的有功功率与所述换流器承受的电压占比成正比。可由下式：
进行直流电压和有功功率的计算。
可选的，在上述实施例的基础上，对通过牛顿-拉夫逊法进行潮流计算的方法进行了优化。图4是本发明实施例体用的一种通过牛顿-拉夫逊法进行潮流计算方法的流程图。如图4，根据任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算，具体包括：
S1511、根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量；
S1512、根据所述换流器的控制方式以及所述有功功率的不平衡量和所述无功注入量的不平衡量，建立所述潮流计算的雅可比矩阵；其中，在所述换流器的流器的控制方式为第二类控制控制方式，且为恒定换相角的控制方式时，与所述换流器对应的节点的雅可比矩阵参数仅由交流网络参数求出；所述换流器的控制方式为第二类控制变量的变压器变比的控制方式时，与所述换流器对应的节点的雅可比矩阵参数由交流网络参数求出后进行修正；
S1513、根据所述雅可比矩阵，通过牛顿-拉夫逊法进行潮流计算。
具体的，在交直流混联电力系统中不同换流器具有不同的控制方式，且换流器对应节点的有功功率和无功注入量与该换流器的控制方式相关。例如，通常潮流计算迭代过程中，不平衡方程为：
其中，Pidk和Qidk均是标量，为正值，其±号的选择原则为整流侧选择正，逆变侧选择负；Pis和Qis为系统发电机与负荷节点总的注入功率；δij为节点i与j之间的相角差；Gij和Bij为导纳矩阵对应元素的实部与虚部。
构建雅可比矩阵如下：
其中，H、N、L为所述雅可比矩阵的分块矩阵，ΔP为所述有功功率不平衡量，ΔQ为所述无功注入量的不平衡量，Δθ和ΔV为迭代过程中的变量的修正量。
当节点对应的换流器控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，无功注入量计算公式如下：
此时，无功注入量的不平衡量应通过下式计算得出：
将所计算的无功注入量代入至雅可比矩阵中，通过牛顿-拉夫逊法完成交直流混联电力系统的潮流计算量。
可选的，在上述实施例的基础上，对通过牛顿-拉夫逊法进行潮流计算的方法可进行进一步的优化。图5是本发明实施例提供的又一种通过牛顿-拉夫逊法进行潮流计算方法的流程图。如图5，通过牛顿-拉夫逊法进行潮流计算，包括：
S1521、根据所述电阻、所述直流母线的电压、所述有功功率和所述无功注入量，获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量；
S1522、根据所述换流器的控制方式以及所述有功功率的不平衡量和所述无功注入量的不平衡量，建立所述潮流计算的雅可比矩阵；
S1523、在所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，对雅可比矩阵元素Lii进行修正；
S1524、判断所述潮流计算量是否满足收敛条件；若是则转至S1525；若否则转至S1521，重新获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量。
S1525、完成潮流计算。
具体的，雅可比矩阵元素Lii的计算公式如下：
而在换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，还需对雅可比矩阵元素Lii进行修正，如下：
其中，i为与换流器连接的交流网路的节点；Vi为对应节点i的电压幅值，Gij和Bij为导纳矩阵的实部与虚部，Va为与换流器相连节点的电压幅值，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，kT为变压器变比，ky为换流器常数，θij为节点i的控制角，H、N、L为雅可比矩阵的分块矩阵，ΔP为有功功率不平衡量，ΔQ为无功注入量的不平衡量，Δθ和ΔV为迭代过程中的变量的修正量。
修正后需对潮流计算量的收敛性进行验证，在潮流计算量满足收敛条件时，结束潮流计算的迭代过程，输出相应的结果；而在潮流计算量不满足收敛条件时，需要重新计算潮流计算过程中的有功功率的不平衡量和无功注入量的不平衡量，直至潮流计算量满足收敛条件。
本发明实施例通过分析交直流混联电力系统中各换流器的控制方式，获得相应的参数，以通过牛顿-拉夫逊法进行潮流计算，能够避免因初值选取规模大、雅可比矩阵规模大，不利于计算的问题，从而使得交直流混联潮流计算的具有较好的收敛性，降低计算的复杂性，进一步提高计算速率、降低成本。
本发明实施例还提供了一种交直流混联潮流计算的装置，该装置可适用于对交直流混联电力系统进行潮流计算的情况。本实施的交直流混联潮流计算的装置可由软件和/或硬件来实现。图6是本发明实施例提供的一种交直流混联潮流计算的装置的结构框图。如图6，交直流混联潮流计算的装置包括电阻获取模块61、控制方式获取模块62、直流电压和有功功率获取模块63、无功注入量计算模块64和潮流计算模块65。
所述电阻获取模块61，用于对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；
所述控制方式获取模块62，用于获取各所述换流器的控制方式；
所述直流电压和有功功率获取模块63，用于根据所述换流器的控制方式，获取与所述换流器对应的节点的直流电压和有功功率；
所述无功注入量计算模块64，用于根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量；
所述潮流计算模块65，用于根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流母线的电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算。
本发明实施例通过分析交直流混联电力系统中各换流器的控制方式，获得相应的参数，以通过牛顿-拉夫逊法进行潮流计算，能够避免因初值选取规模大、雅可比矩阵规模大，不利于计算的问题，从而使得交直流混联潮流计算的具有较好的收敛性，降低计算的复杂性，进一步提高计算速率、降低成本。
本发明实施例还提供了一种包含计算机可执行指令的存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明实施例体用的交直流混联潮流计算的方法，该方法包括：
对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；
根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率；
获取各所述换流器的控制方式；
根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量；
根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量和所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算。
存储介质——任何的各种类型的存储器设备或存储设备。术语“存储介质”旨在包括：安装介质，例如CD-ROM、软盘或磁带装置；计算机系统存储器或随机存取存储器，诸如DRAM、DDR RAM、SRAM、EDO RAM，兰巴斯(Rambus)RAM等；非易失性存储器，诸如闪存、磁介质(例如硬盘或光存储)；寄存器或其它相似类型的存储器元件等。存储介质可以还包括其它类型的存储器或其组合。另外，存储介质可以位于程序在其中被执行的第一计算机系统中，或者可以位于不同的第二计算机系统中，第二计算机系统通过网络(诸如因特网)连接到第一计算机系统。第二计算机系统可以提供程序指令给第一计算机用于执行。术语“存储介质”可以包括可以驻留在不同位置中(例如在通过网络连接的不同计算机系统中)的两个或更多存储介质。存储介质可以存储可由一个或多个处理器执行的程序指令(例如具体实现为计算机程序)。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的交直流混联潮流计算的方法的操作，还可以执行本发明任意实施例所提供的交直流混联潮流计算的方法中的相关操作。
本发明实施例还提供了一种终端，该终端中可集成本发明实施例提供的交直流混联潮流计算的。图7为本发明实施例提供的一种终端的结构示意图。如图7所示，该终端可以包括：显示器(图中未示出)、存储器101、中央处理器(Central Processing Unit，CPU)102(又称处理器，以下简称CPU)、电路板(图中未示出)和电源电路(图中未示出)。所述CPU102和所述存储器101设置在所述电路板上；所述电源电路，用于为所述终端的各个电路或器件供电；所述存储器101，用于存储计算机程序；所述CPU102读取并执行所述存储器101中存储的计算机程序。所述CPU102在执行所述计算机程序时实现以下步骤：对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率；获取各所述换流器的控制方式；根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量；根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量和所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算。
应该理解的是，图示终端100仅仅是终端的一个范例，并且终端100可以具有比图中所示出的更多的或者更少的部件，可以组合两个或更多的部件，或者可以具有不同的部件配置。图中所示出的各种部件可以在包括一个或多个信号处理和/或专用集成电路在内的硬件、软件、或硬件和软件的组合中实现。终端800例如可以为计算机。
本发明实施例提供的终端，实现并行执行本发明任意实施例提供的交直流混联潮流计算的方法的操作，可以有效地对交直流混联电力系统进行潮流计算。
上述实施例中提供的交直流混联潮流计算的装置、存储介质及终端可执行本发明任意实施例所提供的交直流混联潮流计算的方法，具备执行该方法相应的功能模块和有益效果。未在上述实施例中详尽描述的技术细节，可参见本发明任意实施例所提供的交直流混联潮流计算的方法。
本发明实施例还提供了一种计及直流分层结构交直流混联潮流计算的方法，具体包括：
步骤1：对直流网络求解电导矩阵，得到每两个换流器间的电阻，或者与分层结构的直流电网连接点间的电阻值。
步骤2：分析各换流器的控制方式，通过下式
得到每个节点的电压和有功功率。
其中：Idk为流入换流站k的直流电流，Gkj为对应节点k，j之间的导纳矩阵元素。
步骤3：根据步骤2所得到的控制方式和各节点的电压和有功功率，计算对交流电网的无功注入量。
所述步骤3包括：
步骤3.1：如果一个节点对应的换流器控制方式为恒定换相角，则依据式(2)
计算无功注入量，接着转至步骤3.3；否则转至步骤3.2。
步骤3.2：如果一个节点对应的换流器控制方式为恒定变压器变比，则依据式(3)计算无功注入量，并计算依据下式
完成无功注入量对所对应交流电压的导数。
步骤3.3：如果存在分层结构，依据式(4)计算各层对交流电网连接点的功率影响。
步骤4：使用牛顿-拉夫逊法完成潮流计算。
所述步骤4包括：
步骤4.1：设定交流网络初始值，求解潮流功率方程不平衡量。
步骤4.2：构造雅可比矩阵如式
其中：J为雅可比矩阵；当换流器控制方式为恒定换相角时，对应的节点所对应的Lii仅由交流网络参数求出；当换流器控制方式为恒定变压器变比时，对应的节点所对应的Lii除了由交流网络参数求出，需要进一步根据(5)修正；ΔP，ΔQ为功率方程不平衡量，Δθ，ΔV为迭代过程中的变量的修正量。
步骤4.3：求得修正量，完成交流网络参数修正，检查收敛条件，达到条件结束迭代，否则转步骤4.1。
步骤5：输出结果。
具体的，直流网络建模
1)传统直流网络模型
在直流线路参数中，包含电容、电感等参数值。但潮流计算考虑的是稳态情况，故直流线路整体表现为电阻特性。使用节点导纳矩阵表示直流网络：
节点注入电流可以表示为：
Id＝GdVd
式中Id为直流节点注入电流，Vd为直流电压。
换流器基本方程：
式中，小标*表示标幺值，Vdci和Idci分别是直流输电电压和电流；Vi∠δsi是交流母线的线电压向量；Ici是注入换流器的基频交流电流；nti是换流器所含桥的个数；kTi为变压器变比；θi为换流器的控制角，即整流器的触发延迟角、逆变器的熄弧超前角；Xci为换相电阻；考虑到换相角的影响，简化分析，近似取常数kγ＝0.995；
为换流器从交流系统吸收有功功率(整流器为吸收，逆变器为发出)和无功功率对应的功率因素角。
2)直流分层结构的参与
简单的分层接入直流输电结构如图2所示；
分层接入方式直流输电的节点存在串联耦合，如图2中的直流节点1和直流节点2满足如下关系：
其中，各变量含义与式(1)一致。
串联侧单个换流器输出的有功功率与该换流器承受的电压占比成正比，即：
式中：kidk为换流器k在分层结构中的电压占比，Pd为直流注入换流站的有功功率，Vd为换流站所连接节点的直流电压，Pidk为分层结构下换流器k所输出的有功功率，Vidk为分层结构下换流器k所承受的直流电压。
3)换流站控制策略
对于传统的换相换流器，每个换流器有两个独立的控制变量。假定变压器抽头可以无缝调节，那么匝数比kT可以被线性控制。因此，直流母线有功功率Pdc、直流电压Vdc和直流电流Idc被定义为D轴控制变量；变压器变比kT和换流器的控制角θ被称为E轴控制变量。
表1换流器控制策略
D轴控制：对于一个直流网络而言，必有一端换流器D轴控制为电压控制模式，而对于其他端换流器D轴控制无论是恒Pdc还是恒Idc，在已知直流网络电阻的情况下得到Gkj，继而依据式(6)计算得到各端换流器的电压值和有功功率。
2)E轴控制：E轴控制存在两类：
①换流器选择恒定换相角：
抽出功率可以表示为：
其中：
为换流器功率因素，Vdk和Idk为换流器所连接直流节点的电压和电流。
与式(9)联立可得:
②换流器选择恒定变压器变比：
将式(10)式(9)联立可得:
3.潮流计算
在潮流计算迭代过程中，不平衡方程为：
其中Pidc，Qidc均是标量，为正值，±号的选择：整流侧选择正，逆变侧选择负；Pis和Qis为系统发电机与负荷节点总的注入功率；δij为节点i与j之间的相角差；Gij和Bij为导纳矩阵对应元素的实部与虚部。
1)换流器E轴控制选择恒定换相角时，原潮流计算所用的雅可比矩阵无需变化。
2)换流器E轴控制选择恒定变压器变比时
在潮流计算迭代过程中，不平衡量与(11)一致。
雅可比矩阵做如下修改：
式中i对应的是换流器连接的交流节点；
为传统纯交流潮流计算中雅可比矩阵元素L的计算式，Vi为对应节点i的电压幅值，Gij和Bij为导纳矩阵的实部与虚部，Va换流器相连节点的电压幅值，数值上与Vi一致。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整、相互结合和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (18)
1.一种交直流混联潮流计算的方法，其特征在于，包括：
对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；
根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率；
获取各所述换流器的控制方式；所述换流器的控制方式包括第一类控制方式和第二类控制方式；所述第一类控制方式包括恒定有功功率、恒定直流电压和恒定直流电流；所述第二类控制方式包括恒定变压器变比和恒定换相角；
根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量；
根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，获取所述有功功率的不平衡量和所述无功注入量的不平衡量；
根据所述换流器的控制方式以及所述有功功率的不平衡量和所述无功注入量的不平衡量，建立所述潮流计算的雅可比矩阵；其中，在所述换流器的控制方式为第二类控制方式，且为恒定换相角的控制方式时，与所述换流器对应的节点的雅可比矩阵参数仅由交流网络参数求出；所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，与所述换流器对应的节点的雅可比矩阵参数由交流网络参数求出后进行修正；对所述雅可比矩阵元素Lii进行修正，如下：
其中，i为与所述换流器连接的交流网路的节点；Vi为对应节点i的电压幅值，Gij和Bij为导纳矩阵的实部与虚部，Va为与所述换流器相连节点的电压幅值，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，kT为变压器变比，ky为换流器常数，θij为节点i的控制角，H、N、L为所述雅可比矩阵的分块矩阵，ΔP为所述有功功率不平衡量，ΔQ为所述无功注入量的不平衡量，Δθ和ΔV为迭代过程中的变量的修正量；
根据所述雅可比矩阵，判断所述潮流计算量是否满足收敛条件；
若是则完成潮流计算；
若否，则重新获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量。
2.根据权利要求1所述的方法，其特征在于，所述根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率，包括：
根据所述直流网络的结构，获取与所述换流器对应的节点参数；
根据所述节点参数，构建方程组，如下：
其中，Vdk为与所述换流器对应的节点的直流电压，Pdk为与所述换流器对应的节点的有功功率，Idk为流入换流站k的直流电流，Gkj为对应节点k和j之间的导纳矩阵元素，Vdj为与换流器j连接的直流母线的电压，nc为所述直流网络中换流器的个数；
根据所述方程组，求解与所述换流器对应的节点的直流电压和有功功率。
3.根据权利要求2所述的方法，其特征在于，在所述交直流混联电力系统的直流网络包括分层结构时，所述直流网络串联侧单个换流器输出的有功功率与所述换流器承受的电压占比成正比。
4.根据权利要求3所述的方法，其特征在于，在所述交直流混联电力系统的直流网络包括分层结构时，满足如下关系：
其中，Idi1和Idi2分别为流经分层结构下换流站的高电压换流器的电流和低电压换流器的电流；Id表示流经换流站整体的电流；Vdr为直流网络送端电压；Vdi1和Vdi2分别表示分层结构下高、低电压换流器的直流电压；Rd为直流线路的电阻；
所述根据所述方程组，求解与所述换流器对应的节点的直流电压和有功功率，具体为：
其中，kidk为换流器k在分层结构中的电压占比，Pd为直流注入换流站的有功功率，Vd为换流站所连接节点的直流电压，Pidk为分层结构下换流器k所输出的有功功率，Vidk为分层结构下换流器k所承受的直流电压。
5.根据权利要求1所述的方法，其特征在于，所述换流器的控制方式为所述第一类控制方式时，所述根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算所述换流器对交流电网的无功注入量，具体为：
其中，Idk为流入换流器k的直流电流，Pdk为所述有功功率，Vdk为所述直流电压，
为换流器功率因素，Qdk为所述无功注入量。
6.根据权利要求1所述的方法，其特征在于，所述换流器的控制方式为第二类控制方式，且为恒定换相角的控制方式时，所述根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量，具体为：
其中，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，Pidk为直流注入交流节点i的有功功率，θd为换流器的控制角，Xc为换相电阻，ky为换流器常数，Qdk为所述无功注入量。
7.根据权利要求1所述的方法，其特征在于，所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，所述根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量，具体为：
其中，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，Va为与换流器相连节点的电压幅值，kT为变压器变比，ky为换流器常数。
8.一种交直流混联潮流计算的装置，其特征在于，包括：
电阻获取模块，用于对交直流混联电力系统的直流网络求解电导矩阵，获取所述直流网络中任意两个换流器之间的电阻，或者获取各分层结构的所述直流网络的连接点间的电阻；
直流电压和有功功率获取模块，用于根据所述直流网络的结构，获取与所述换流器对应的节点的直流电压和有功功率；
控制方式获取模块，用于获取各所述换流器的控制方式；所述换流器的控制方式包括第一类控制方式和第二类控制方式；所述第一类控制方式包括恒定有功功率、恒定直流电压和恒定直流电流；所述第二类控制方式包括恒定变压器变比和恒定换相角；
无功注入量计算模块，用于根据所述直流电压、所述有功功率以及所述换流器的控制方式，计算与所述换流器对交流电网的无功注入量；
潮流计算模块，用于根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算；
其中，根据任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，通过牛顿-拉夫逊法进行潮流计算，包括：
根据所述任意两个换流器之间的电阻或所述连接点间的电阻以及所述直流电压、所述有功功率、所述无功注入量以及所述换流器的控制方式，获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量；
根据所述换流器的控制方式以及所述有功功率的不平衡量和所述无功注入量的不平衡量，建立所述潮流计算的雅可比矩阵；其中，在所述换流器的控制方式为第二类控制方式，且为恒定换相角的控制方式时，与所述换流器对应的节点的雅可比矩阵参数仅由交流网络参数求出；所述换流器的控制方式为第二类控制方式，且为恒定变压器变比的控制方式时，与所述换流器对应的节点的雅可比矩阵参数由交流网络参数求出后进行修正；对所述雅可比矩阵元素Lii进行修正，如下：
其中，i为与所述换流器连接的交流网路的节点；Vi为对应节点i的电压幅值，Gij和Bij为导纳矩阵的实部与虚部，Va为与所述换流器相连节点的电压幅值，Vdk为直流输电电压，Pdk为流入换流器k的有功功率，kT为变压器变比，ky为换流器常数，θij为节点i的控制角，H、N、L为所述雅可比矩阵的分块矩阵，ΔP为所述有功功率不平衡量，ΔQ为所述无功注入量的不平衡量，Δθ和ΔV为迭代过程中的变量的修正量；
根据所述雅可比矩阵，判断所述潮流计算量是否满足收敛条件；
若是则完成潮流计算；
若否，则重新获取所述潮流计算过程中的所述有功功率的不平衡量和所述无功注入量的不平衡量。
9.一种存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1～7中任一所述的交直流混联潮流计算的方法。
10.一种终端，包括显示屏、存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1～7中任一所述的交直流混联潮流计算的方法。
11.一种计及直流分层结构交直流混联潮流计算的方法，其特征是：根据换流站的控制方式完成直流网络的状态计算，将直流网络与交流网络的连接点等效为功率节点；使用牛顿—拉夫逊法计算潮流；
其中，使用牛顿-拉夫逊法完成潮流计算的具体过程包括：
(a)设定交流网络初始值，求解潮流功率方程不平衡量；
(b)构造雅可比矩阵，流器控制方式为恒定换相角对应的节点所对应的雅可比矩阵参数仅由交流网络参数求出；换流器控制方式为恒定变压器变比对应的节点所对应的雅可比矩阵参数由交流网络参数求出后进行修正；修正的方式为：
其中：-Vi∑j∈i，j≠iVj(Gijsinθij-Bijcosθij)+2Vi 2Bij为传统纯交流潮流计算中雅可比矩阵元素L的计算式，Vi为对应节点i的电压幅值；Gij和Bij为导纳矩阵的实部与虚部；Va为与所述换流器相连节点的电压幅值；Vdc为所述换流器的直流电压；Pdc为所述换流器的有功功率；kT为变压器变比；ky为换流器常数；θij为节点i的控制角；
(c)完成交流网络参数修正，检查收敛条件，达到条件结束迭代，否则转步骤(a)。
12.一种计及直流分层结构交直流混联潮流计算的方法，其特征是：包括以下步骤：
对直流网络求解电导矩阵，得到每两个换流器间的电阻，或者与分层结构的直流电网连接点间的电阻值；
分析各换流器的控制方式，计算得到每个节点的电压和有功功率；
根据所得到的控制方式和各节点的电压和有功功率，计算对交流电网的无功注入量；
使用牛顿-拉夫逊法完成潮流计算，得到计算结果；
其中，使用牛顿-拉夫逊法完成潮流计算的具体过程包括：
(a)设定交流网络初始值，求解潮流功率方程不平衡量；
(b)构造雅可比矩阵，流器控制方式为恒定换相角对应的节点所对应的雅可比矩阵参数仅由交流网络参数求出；换流器控制方式为恒定变压器变比对应的节点所对应的雅可比矩阵参数由交流网络参数求出后进行修正；修正的方式为：
其中：-Vi∑j∈i，j≠iVj(Gijsinθij-Bijcosθij)+2Vi 2Bij为传统纯交流潮流计算中雅可比矩阵元素L的计算式；Vi为对应节点i的电压幅值；Gij和Bij为导纳矩阵的实部与虚部；Va为与所述换流器相连节点的电压幅值；Vdc为所述换流器的直流电压；Pdc为所述换流器的有功功率；kT为变压器变比；ky为换流器常数；θij为节点i的控制角；
(c)完成交流网络参数修正，检查收敛条件，达到条件结束迭代，否则转步骤(a)。
13.根据权利要求11或12所述的一种计及直流分层结构交直流混联潮流计算的方法，其特征是：依据换流器控制方式确定节点部分参数，继而构造方程组
其中：Idk为流入换流站k的直流电流；Gkj为对应节点k，j之间的导纳矩阵元素；Vdk为与所述换流器对应的节点的直流电压；Pdk为与所述换流器对应的节点的有功功率，Vdj为与换流器j连接的直流母线的电压；
以求解节点电压和有功功率为目标完成计算，具体方程组为：
其中：Idi1，Idi2分别表示流经分层结构下换流站的高、低电压换流器的电流；Id表示流经换流站整体的电流；Vdr表示直流网络送端电压；Vdi1、Vdi2分别表示分层结构下高、低电压换流器的直流电压；Rd表示直流线路的电阻。
14.根据权利要求11或12所述的一种计及直流分层结构交直流混联潮流计算的方法，其特征是：串联侧单个所述换流器输出的有功功率与所述换流器承受的电压占比成正比。
15.根据权利要求11或12所述的一种计及直流分层结构交直流混联潮流计算的方法，其特征是：依据节点换流器的控制方式，即恒定换相角控制方式下使用
恒定变压器变比控制方式下使用
完成对交流电网无功注入量的计算，有功注入量的计算通过直流网络状态计算得到；其中：Pidc为所述换流器的抽出功率；Pdc为所述换流器的有功功率；Vdc为所述换流器的直流电压；Vdc为换流站所连接的直流网络节点电压；θd为换流器的控制角，即整流器的触发延迟角、逆变器的熄弧超前角；kT为变压器变比；Xc为换相电阻；考虑到换相角的影响，引入了变量ky；
为换流器从交流系统吸收有功功率和无功功率对应的功率因素角；Va为换流器连接的交流网络的电压幅值。
16.根据权利要求11或12所述的一种计及直流分层结构交直流混联潮流计算的方法，其特征是：涉及到分层结构时，根据分层的换流器电压占比完成功率节点的等效。
17.根据权利要求11或12所述的一种计及直流分层结构交直流混联潮流计算的方法，其特征是：计算对交流电网的无功注入量的具体过程包括：
(1)如果一个节点对应的换流器控制方式为恒定换相角，则用
计算无功注入量，接着转至步骤(3)；否则转至步骤(2)；
其中：Pidc为所述换流器的抽出功率；Pdc为所述换流器的有功功率；Vdc为所述换流器的直流电压；Vdc为换流站所连接的直流网络节点电压；θd为换流器的控制角，即整流器的触发延迟角、逆变器的熄弧超前角；kt为变压器变比；Xc为换相电阻；考虑到换相角的影响，引入了变量ky；
为换流器从交流系统吸收有功功率和无功功率对应的功率因素角；
(2)如果一个节点对应的换流器控制方式为恒定变压器变比，用
计算无功注入量，并计算无功注入量对所对应交流电压的导数；
Va为换流器连接的交流网络的电压幅值；
(3)如果存在分层结构，依据
计算各层对交流电网连接点的功率影响；
式中：kidk为换流器k在分层结构中的电压占比，Pd为直流注入换流站的有功功率，Vd为换流站所连接节点的直流电压，Pidk为分层结构下换流器k所输出的有功功率，Vidk为分层结构下换流器k所承受的直流电压。
18.一种计及直流分层结构交直流混联潮流计算的系统，其特征是：运行于处理器上，被配置为执行以下指令：
对直流网络求解电导矩阵，得到每两个换流器间的电阻，或者与分层结构的直流电网连接点间的电阻值；
分析各换流器的控制方式，计算得到每个节点的电压和有功功率；
根据所得到的控制方式和各节点的电压和有功功率，计算对交流电网的无功注入量；
使用牛顿-拉夫逊法完成潮流计算，得到计算结果；
其中，使用牛顿-拉夫逊法完成潮流计算的具体过程包括：
(a)设定交流网络初始值，求解潮流功率方程不平衡量；
(b)构造雅可比矩阵，流器控制方式为恒定换相角对应的节点所对应的雅可比矩阵参数仅由交流网络参数求出；换流器控制方式为恒定变压器变比对应的节点所对应的雅可比矩阵参数由交流网络参数求出后进行修正；修正的方式为：
其中：-Vi∑j∈i，j≠iVj(Gijsinθij-Bijcosθij)+2Vi 2Bij为传统纯交流潮流计算中雅可比矩阵元素L的计算式；Vi为对应节点i的电压幅值；Gij和Bij为导纳矩阵的实部与虚部；Va为与所述换流器相连节点的电压幅值；Vdc为所述换流器的直流电压；Pdc为所述换流器的有功功率；kT为变压器变比；ky为换流器常数；θij为节点i的控制角；
(c)完成交流网络参数修正，检查收敛条件，达到条件结束迭代，否则转步骤(a)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US16/538,673 US20200081044A1 (en) | 2018-09-07 | 2019-08-12 | Power flow calculation method and device for ac-dc interconnected power system, storage medium and terminal |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811045461.5A CN109462237A (zh) | 2018-09-07 | 2018-09-07 | 计及直流分层结构交直流混联潮流计算的方法及系统 |
| CN2018110454615 | 2018-09-07 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109980648A CN109980648A (zh) | 2019-07-05 |
| CN109980648B true CN109980648B (zh) | 2023-04-07 |
Family
ID=65606584
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811045461.5A Pending CN109462237A (zh) | 2018-09-07 | 2018-09-07 | 计及直流分层结构交直流混联潮流计算的方法及系统 |
| CN201910309632.9A Active CN109980648B (zh) | 2018-09-07 | 2019-04-17 | 交直流混联潮流计算的方法及其装置、存储介质及终端 |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811045461.5A Pending CN109462237A (zh) | 2018-09-07 | 2018-09-07 | 计及直流分层结构交直流混联潮流计算的方法及系统 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20200081044A1 (zh) |
| CN (2) | CN109462237A (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110020510B (zh) * | 2019-06-06 | 2019-09-10 | 西南交通大学 | 一种基于Floquet理论的动车组网侧单相脉冲整流器的稳定性分析方法 |
| CN110718918B (zh) * | 2019-09-12 | 2021-01-05 | 重庆大学 | 一种基于hvac和lcc-vsc hvdc混合系统的潮流计算方法 |
| CN110912140B (zh) * | 2019-12-13 | 2023-04-28 | 福州大学 | 计及电解水制氢系统精细化模型的配电网潮流计算方法 |
| CN111416358B (zh) * | 2020-04-17 | 2021-10-08 | 清华大学 | 交直流混联系统的潮流分析装置和潮流分析方法 |
| CN112421638B (zh) * | 2020-11-19 | 2022-07-22 | 国家电网公司西南分部 | 一种基于潮流贡献度的直流近区资产划分方法 |
| CN114629122A (zh) * | 2020-12-14 | 2022-06-14 | 国网河南省电力公司经济技术研究院 | 一种交直流混合配电网线性潮流计算方法 |
| CN112636381B (zh) * | 2020-12-21 | 2022-10-21 | 中国科学院电工研究所 | 一种手拉手型的交直流配用电系统运行稳定分析方法 |
| CN112865076B (zh) * | 2021-01-12 | 2022-05-20 | 天津大学 | 面向电力电子化柔性配电网的改进增广节点方程建模方法 |
| CN115995818A (zh) * | 2023-03-22 | 2023-04-21 | 国网湖北省电力有限公司经济技术研究院 | 一种多能互补交直流混合配电网潮流计算方法 |
| CN118399417A (zh) * | 2024-04-23 | 2024-07-26 | 山东大学 | 含lcc-vsc串并联的混合直流潮流计算方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103543357A (zh) * | 2013-10-23 | 2014-01-29 | 华北电力大学 | 一种换流阀换相失败的预判方法 |
| CN103904657A (zh) * | 2014-03-24 | 2014-07-02 | 国家电网公司 | 一种基于支持向量机的地区电网无功电压控制方法 |
| CN105759081A (zh) * | 2016-01-21 | 2016-07-13 | 中国电力科学研究院 | 一种高压直流输电系统避雷器老化试验电压波形的产生方法 |
| CN106532711A (zh) * | 2016-11-30 | 2017-03-22 | 大连海事大学 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
| CN106549384A (zh) * | 2016-12-09 | 2017-03-29 | 国网江苏省电力公司经济技术研究院 | 一种含upfc电力系统的通用潮流计算方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2688191A1 (en) * | 2012-07-17 | 2014-01-22 | ABB Research Ltd. | Multi terminal HVDC control |
-
2018
- 2018-09-07 CN CN201811045461.5A patent/CN109462237A/zh active Pending
-
2019
- 2019-04-17 CN CN201910309632.9A patent/CN109980648B/zh active Active
- 2019-08-12 US US16/538,673 patent/US20200081044A1/en not_active Abandoned
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103543357A (zh) * | 2013-10-23 | 2014-01-29 | 华北电力大学 | 一种换流阀换相失败的预判方法 |
| CN103904657A (zh) * | 2014-03-24 | 2014-07-02 | 国家电网公司 | 一种基于支持向量机的地区电网无功电压控制方法 |
| CN105759081A (zh) * | 2016-01-21 | 2016-07-13 | 中国电力科学研究院 | 一种高压直流输电系统避雷器老化试验电压波形的产生方法 |
| CN106532711A (zh) * | 2016-11-30 | 2017-03-22 | 大连海事大学 | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 |
| CN106549384A (zh) * | 2016-12-09 | 2017-03-29 | 国网江苏省电力公司经济技术研究院 | 一种含upfc电力系统的通用潮流计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| 分层接入方式的特高压直流输电逆变侧;卢东斌 等;《中国电机工程学报》;20160405;第36卷(第7期);第1808-1816页 * |
| 基于模块化多电平换流器的交直流;韦延方 等;《电源学报》;20160331;第14卷(第2期);第112-120页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109462237A (zh) | 2019-03-12 |
| CN109980648A (zh) | 2019-07-05 |
| US20200081044A1 (en) | 2020-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109980648B (zh) | 交直流混联潮流计算的方法及其装置、存储介质及终端 | |
| Li et al. | Optimal power flow in stand-alone DC microgrids | |
| US10222404B2 (en) | General load flow calculation method for power systems with unified power flow controller | |
| US20060282239A1 (en) | Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system | |
| CN110932282B (zh) | 一种增广直角坐标下的基于vsc内部修正方程矩阵与交替迭代法的潮流计算方法 | |
| CN115036962B (zh) | 柔性直流输电系统潮流计算及交直流混联潮流计算方法 | |
| CN106953329B (zh) | 含直流系统的复杂拓扑有源配电网潮流计算方法 | |
| CN108711846B (zh) | 一种交直流系统长期电压稳定模型预测控制方法 | |
| CN110556853B (zh) | 一种为电磁暂态仿真提供初值的计算方法及系统 | |
| CN108347057B (zh) | 一种lcc-mmc混合直流电网交替迭代潮流计算方法 | |
| CN108988344B (zh) | 一种电力系统的潮流计算优化方法 | |
| Dong et al. | Output control method of microgrid VSI control network based on dynamic matrix control algorithm | |
| CN109149665A (zh) | 用于柔性直流交流电网联合仿真的多速率仿真方法及系统 | |
| CN109149583A (zh) | 有源配电网简洁前推回代潮流计算方法 | |
| CN112670992B (zh) | 含能量路由器的配电网稳定性分析和失稳校正方法及系统 | |
| CN109861226A (zh) | 一种综合优化谐波稳定性和阻尼损耗的lcl滤波器设计方法 | |
| CN109787211A (zh) | 一种基于vsc的多端直流配电网潮流计算方法及系统 | |
| Nguyen et al. | Power flow solution methods for combined transmission and distribution systems | |
| Abdel-Akher et al. | Development of unbalanced three-phase distribution power flow analysis using sequence and phase components | |
| Jiménez et al. | Unbalanced three-phase power flow studies of distribution systems with plug-in electric vehicles | |
| CN112526396A (zh) | 多回线接地故障分析方法、装置、电子设备及存储介质 | |
| CN112039061A (zh) | 基于电磁暂态仿真的潮流计算方法 | |
| Benato et al. | Power flow solution by a complex admittance matrix method | |
| CN113890069B (zh) | 基于可变均衡系数的级联型储能系统控制方法 | |
| CN111509761A (zh) | 一种与锁相环结合的hvdc动态相量计算方法及系统 |
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 |