CN106953331A - 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 - Google Patents
一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 Download PDFInfo
- Publication number
- CN106953331A CN106953331A CN201710333984.9A CN201710333984A CN106953331A CN 106953331 A CN106953331 A CN 106953331A CN 201710333984 A CN201710333984 A CN 201710333984A CN 106953331 A CN106953331 A CN 106953331A
- Authority
- CN
- China
- Prior art keywords
- node
- harmonic
- value
- phase
- current
- 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.)
- Pending
Links
- 238000001228 spectrum Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000004075 alteration Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
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/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
-
- 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/01—Arrangements for reducing harmonics or ripples
-
- 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/26—Arrangements for eliminating or reducing asymmetry in polyphase 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/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/40—Arrangements for reducing harmonics
-
- 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/50—Arrangements for eliminating or reducing asymmetry in polyphase networks
Abstract
本发明涉及一种考虑不确定性和三相不平衡的电力系统谐波潮流算法,包括:选择谐波次数h;计算非线性负荷节点所接谐波源注入系统的第h次谐波电流值;建立系统元件包括各线路及变压器等的第h次谐波模型;初始化各节点谐波电压与相角;设置迭代次数,每次迭代过程,计算各线性负荷节点的等效注入谐波电流向值;从末节点开始,通过对支路谐波电流的求和计算,获得各条支路每相的谐波电流值;从馈线始端至馈线末端依次计算节点谐波电压值;迭代终止判据为各节点每相的谐波电压相对于上一次迭代的数值偏差小于允许值;根据所获基波电压值以及各次谐波电压的值求取各节点各相谐波电压畸变率。
Description
技术领域
本发明涉及电力系统谐波分析领域,特别涉及基于不确定数学理论的电力系统三相不平衡谐波分布计算方法。
技术背景
随着电力市场的出现以及分布式电源的接入,在电力系统实际运行中会存在大量的不确定因素,如分布式电源出力的不确定性、变动负荷的不确定性、电力原件模型参数的不确定性等,这些因素都会对系统的谐波水平造成影响,进而导致谐波分布具有不确定性,传统的确定性谐波分布计算方法在对具有不确定性信息的系统进行分析计算时,表现出一定的局限性。不确定谐波分布计算在电力系统的谐波水平评估中已有了一些应用,但由于问题本身难度较大,它一直是谐波分布计算中的一个难点。
目前针对不确定谐波分布的计算主要包括概率谐波分布方法和模糊谐波分布方法。对于具备统计性质的随机信息,可以用概率的理论来表示负荷的不确定性;对于不具备统计性质的不确定因素,则可以采用模糊数学的方法来处理。但这些详细的信息并不容易获得,对于负荷预测的结果不精确但在一定范围内准确时,采用区间数学来求解非常适用。
发明内容
本发明的目的是提供谐波分布分析过程中对不确定信息的处理方法,在谐波源详细的概率分布无法获得的情况下,能以区间表征节点谐波的水平,提供基于不确定性数学理论的三相不平衡谐波分布计算方法。本发明采用的技术方案是:
一种考虑不确定性和三相不平衡的电力系统谐波潮流算法,包括下列步骤:
(1)选择谐波次数h;
(2)根据基波计算结果,计算非线性负荷节点所接谐波源注入系统的第h次谐波电流值
(3)计及谐波影响,建立系统元件包括各线路及变压器等的第h次谐波模型;
(4)初始化各节点谐波电压与相角;
(5)将迭代次数k置零,k=0;
(6)计算各线性负荷节点的等效注入谐波电流向值:
式中:是线性负荷节点j的第h次谐波电压第k次迭代的值;
是线性负荷节点j的第h次谐波等效阻抗;
(7)从末节点开始,通过对支路谐波电流的求和计算,获得各条支路每相的谐波电流值:
式中:是节点j-1到节点j第i相的谐波线路电流值;
是节点j第i相的谐波注入电流值;
Ωj是节点j所连支路的集合;
(8)从馈线始端至馈线末端依次计算节点谐波电压值:
式中:是节点j的谐波电压值矩阵;
是节点j-1谐波电压值矩阵;
是节点j-1到节点j的线路阻抗矩阵;
(9)判断终止,若满足迭代误差要求则停止,否则返回(6)继续迭代,k=k+1,迭代终止判据为各节点每相的谐波电压相对于上一次迭代的数值偏差小于允许值ε:
(10)判断谐波计算是否结束,若未结束则返回(1)继续选择谐波次数,否则,执行下一步;
(11)根据所获基波电压值以及各次谐波电压的值求取各节点各相谐波电压畸变率:
式中:是节点j第i相的电压谐波畸变率;
是节点j处第i相基波电压值;
是节点j处第i相的第h次谐波电压值。
其中步骤(2)可以如下:
1)将节点注入功率统一表示成不确定量的形式:
式中为节点j第i相负荷值,i=a,b,c;
为节点j第i相的额定功率;
为节点j第i相的噪声元;
为噪声元系数;
2)经三相基波潮流算法,求得各节点电压值进而求取非线性负荷节点的电流值:
式中为节点j第i相的注入电流值;
为节点j第i相的电压值;
3)根据基波电流值和谐波源典型频谱,求取各次谐波电流的幅值和相角,得到谐波源的恒流源不确定性模型:
式中I1-spectrum为基波电流幅值;
Ih-spectrum为第h次谐波电流幅值;
为节点j第i相的基波电流相角值;
θ1-spectrum为基波相角值;
θh-spectrum为第h次谐波相角值;
为节点j第i相第h次谐波电流的相角值;
|·|为模值求取运算。
本发明利用区间数学,在稳态条件下,计及负荷的不确定性,以区间来表示谐波分布情况,获得配网中各节点谐波电压的波动范围,方便对非测量点的谐波水平进行合理估计,建立基于不确定性数学理论的三相不平衡谐波分布计算方法,为治理配网电能质量问题提供合理的依据。
附图说明
图1基于不确定性数学理论的三相不平衡电力系统谐波分布计算流程图。
图2 33节点算例接线图。
具体实施方式
本发明的考虑不确定性和三相不平衡的电力系统谐波潮流算法具体步骤如下
1.根据需求获取基波潮流计算的结果;
2.选择谐波次数h;
3.根据基波计算结果,结合测量、仿真、典型频谱等方式,计算非线性负荷节点所接谐波源注入系统的第h次谐波电流值方法如下:
1)将节点注入功率统一表示成不确定量的形式:
式中为节点j第i相负荷值,i=a,b,c;
为节点j第i相的额定功率;
为节点j第i相的噪声元;
为噪声元系数;
2)经三相基波潮流算法,求得各节点电压值进而求取非线性负荷节点的电流值:
式中为节点j第i相的注入电流值;
为节点j第i相的电压值;
3)根据基波电流值和谐波源典型频谱,求取各次谐波电流的幅值和相角,得到谐波源的恒流源不确定性模型:
式中I1-spectrum为基波电流幅值;
Ih-spectrum为第h次谐波电流幅值;
为节点j第i相的基波电流相角值;
θ1-spectrum为基波相角值;
θh-spectrum为第h次谐波相角值;
为节点j第i相第h次谐波电流的相角值;
|·|为模值求取运算;
4.计及谐波影响,更新网络参数,建立系统元件包括各线路及变压器等的第h次谐波模型;
5.初始化各节点每相第h次谐波电压的幅值和相角;
6.将迭代次数k置零,k=0;
7.计算线性负荷的等效注入谐波电流值:
式中:是线性负荷节点j的第h次谐波电压第k次迭代的值;
是线性负荷节点j的第h次谐波等效阻抗;
8.从末节点开始,通过对支路谐波电流的求和计算,获得各条支路每相的始端谐波电流值:
式中:是节点j-1到节点j第i相的谐波线路电流值;
是节点j第i相的谐波注入电流值;
Ωj是节点j所连支路的集合;
9.从馈线始端至馈线末端依次计算节点谐波电压值:
式中:是节点j的谐波电压值矩阵;
是节点j-1谐波电压值矩阵;
是节点j-1到节点j的线路阻抗矩阵;
10.判断终止,若满足迭代误差要求则停止,否则返回7继续迭代,k=k+1。迭代终止判据为各节点每相的谐波电压相对于上一次迭代的数值偏差小于允许值:
11.判断谐波计算是否结束,若结束则继续计算THD(%),否则返回(2)继续选择谐波次数;
12.计算谐波畸变率,根据所获基波电压复仿射值以及各次谐波电压的复仿射值求取各节点各相谐波电压畸变率:
式中:是节点j第i相的电压谐波畸变率值;
是节点j处第i相基波电压值;
是节点j处第i相的第h次谐波电压值。
13.输出计算结果。
以某33节点电力系统为例进一步说明本发明,表1给出算例中接入的两种谐波源电流频谱,表2为此算例非线性负荷节点信息,系统其它数据与33节点算例相同,不再赘述。本发明采用有名值进行计算(误差精度为0.01),输出的结果以有名值的形式表示。
表3为33节点配电系统三相不平衡谐波分布各节点各相THD(%)的计算结果。
表1
表2
表3
计算结果表明,本发明提出的基于不确定性理论的三相不平衡谐波分布计算方法能够获得系统不确定谐波分布的区间解,解决了确定性谐波分布不能解决的问题。
Claims (2)
1.一种考虑不确定性和三相不平衡的电力系统谐波潮流算法,包括下列步骤:
(1)选择谐波次数h;
(2)根据基波计算结果,计算非线性负荷节点所接谐波源注入系统的第h次谐波电流值
(3)计及谐波影响,建立系统元件包括各线路及变压器等的第h次谐波模型;
(4)初始化各节点谐波电压与相角;
(5)将迭代次数k置零,k=0;
(6)计算各线性负荷节点的等效注入谐波电流向值:
式中:是线性负荷节点j的第h次谐波电压第k次迭代的值;
是线性负荷节点j的第h次谐波等效阻抗;
(7)从末节点开始,通过对支路谐波电流的求和计算,获得各条支路每相的谐波电流值:
式中:是节点j-1到节点j第i相的谐波线路电流值;
是节点j第i相的谐波注入电流值;
Ωj是节点j所连支路的集合;
(8)从馈线始端至馈线末端依次计算节点谐波电压值:
式中:是节点j的谐波电压值矩阵;
是节点j-1谐波电压值矩阵;
是节点j-1到节点j的线路阻抗矩阵;
(9)判断终止,若满足迭代误差要求则停止,否则返回(6)继续迭代,k=k+1,迭代终止判据为各节点每相的谐波电压相对于上一次迭代的数值偏差小于允许值ε:
(10)判断谐波计算是否结束,若未结束则返回(1)继续选择谐波次数,否则,执行下一步;
(11)根据所获基波电压值以及各次谐波电压的值求取各节点各相谐波电压畸变率:
式中:是节点j第i相的电压谐波畸变率;
是节点j处第i相基波电压值;
是节点j处第i相的第h次谐波电压值。
2.根据权利要求1所述的算法,其特征在于,步骤2)如下:
1)将节点注入功率统一表示成不确定量的形式:
式中为节点j第i相负荷值,i=a,b,c;
为节点j第i相的额定功率;
为节点j第i相的噪声元;
为噪声元系数;
2)经三相基波潮流算法,求得各节点电压值进而求取非线性负荷节点的电流值:
式中为节点j第i相的注入电流值;
为节点j第i相的电压值;
3)根据基波电流值和谐波源典型频谱,求取各次谐波电流的幅值和相角,得到谐波源的恒流源不确定性模型:
式中I1-spectrum为基波电流幅值;
Ih-spectrum为第h次谐波电流幅值;
为节点j第i相的基波电流相角值;
θ1-spectrum为基波相角值;
θh-spectrum为第h次谐波相角值;
为节点j第i相第h次谐波电流的相角值;
|·|为模值求取运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710333984.9A CN106953331A (zh) | 2017-05-12 | 2017-05-12 | 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710333984.9A CN106953331A (zh) | 2017-05-12 | 2017-05-12 | 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106953331A true CN106953331A (zh) | 2017-07-14 |
Family
ID=59478545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710333984.9A Pending CN106953331A (zh) | 2017-05-12 | 2017-05-12 | 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106953331A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108471119A (zh) * | 2018-04-24 | 2018-08-31 | 长沙理工大学 | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 |
CN111274701A (zh) * | 2020-01-20 | 2020-06-12 | 福州大学 | 一种采用区间监测数据降维回归的谐波源仿射建模方法 |
CN112485523A (zh) * | 2020-11-25 | 2021-03-12 | 云南电网有限责任公司电力科学研究院 | 一种判定谐波电压测量失真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236691A (zh) * | 2013-04-16 | 2013-08-07 | 天津大学 | 基于复仿射数学理论的三相不平衡潮流计算方法 |
CN103972886A (zh) * | 2014-05-12 | 2014-08-06 | 河海大学 | 一种主动配电网谐波潮流分析方法 |
-
2017
- 2017-05-12 CN CN201710333984.9A patent/CN106953331A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236691A (zh) * | 2013-04-16 | 2013-08-07 | 天津大学 | 基于复仿射数学理论的三相不平衡潮流计算方法 |
CN103972886A (zh) * | 2014-05-12 | 2014-08-06 | 河海大学 | 一种主动配电网谐波潮流分析方法 |
Non-Patent Citations (2)
Title |
---|
王守相等: "配电系统三相不确定谐波潮流的复仿射计算方法", 《电力系统自动化》 * |
田洪等: "基于蒙特卡洛模拟的配电网三相潮流算法", 《现代电力》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108471119A (zh) * | 2018-04-24 | 2018-08-31 | 长沙理工大学 | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 |
CN108471119B (zh) * | 2018-04-24 | 2021-07-06 | 长沙理工大学 | 含智慧社区的配电网三相不平衡动态潮流模型预测控制方法 |
CN111274701A (zh) * | 2020-01-20 | 2020-06-12 | 福州大学 | 一种采用区间监测数据降维回归的谐波源仿射建模方法 |
CN111274701B (zh) * | 2020-01-20 | 2022-06-07 | 福州大学 | 一种采用区间监测数据降维回归的谐波源仿射建模方法 |
CN112485523A (zh) * | 2020-11-25 | 2021-03-12 | 云南电网有限责任公司电力科学研究院 | 一种判定谐波电压测量失真方法 |
CN112485523B (zh) * | 2020-11-25 | 2024-04-09 | 云南电网有限责任公司电力科学研究院 | 一种判定谐波电压测量失真方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Smith et al. | A review of iterative harmonic analysis for AC-DC power systems | |
CN104638678B (zh) | 一种多风机并网谐振分析方法 | |
US10164434B2 (en) | Methods and systems for power injection or extraction in a power network | |
CN103969552B (zh) | 一种分布式发电系统的谐波源定位分析方法 | |
CN103236691B (zh) | 基于复仿射数学理论的三相不平衡潮流计算方法 | |
CN109494724B (zh) | 基于lu分解的大电网戴维南等值参数在线辨识方法 | |
CN111082427A (zh) | 一种基于全纯函数的微电网潮流计算方法 | |
CN110556831A (zh) | 一种多机多节点发电系统的稳定性评估方法及系统 | |
CN109066694A (zh) | 含线间潮流控制器的电力系统的多目标潮流优化方法 | |
CN108054757A (zh) | 一种内嵌无功和电压的n-1闭环安全校核方法 | |
CN104616061B (zh) | 基于小波包对数能量熵及遗传算法优化的孤岛检测方法 | |
Sharma et al. | An iterative approach for analysis of harmonic penetration in the power transmission networks | |
CN106953331A (zh) | 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 | |
CN110034567A (zh) | 基于传播特性的配电网短路型电压暂降特征量计算方法 | |
CN113158573B (zh) | 一种基于深度学习的小水电群区域配网最优潮流获取方法 | |
CN104794531A (zh) | 基于网损等值负荷模型的改进直流最优潮流方法 | |
CN109921426A (zh) | 基于cv-kde的风电并网系统概率潮流计算方法 | |
Liu et al. | Frequency‐coupled impedance model based subsynchronous oscillation analysis for direct‐drive wind turbines connected to a weak AC power system | |
CN108075480A (zh) | 一种交直流系统的状态估计方法及系统 | |
CN103956735A (zh) | 一种分布式发电系统的谐波潮流分析方法 | |
CN108418242A (zh) | 一种基于相似度同调的双馈风电机动态等值方法 | |
CN110188480A (zh) | 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析系统和方法 | |
CN113224788B (zh) | 含有风电接入的交直流混联电网静态电压稳定评估方法 | |
CN107658880B (zh) | 基于关联矩阵运算的快速分解法系数矩阵计算方法 | |
CN110048428B (zh) | 基于概率守恒原理的电力系统概率潮流计算方法 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170714 |
|
WD01 | Invention patent application deemed withdrawn after publication |