CN105160058B - 一种用于电网电磁暂态计算的网络边界等值方法 - Google Patents
一种用于电网电磁暂态计算的网络边界等值方法 Download PDFInfo
- Publication number
- CN105160058B CN105160058B CN201510402799.1A CN201510402799A CN105160058B CN 105160058 B CN105160058 B CN 105160058B CN 201510402799 A CN201510402799 A CN 201510402799A CN 105160058 B CN105160058 B CN 105160058B
- Authority
- CN
- China
- Prior art keywords
- network
- bus
- harmonic
- boundary
- node
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 239000000523 sample Substances 0.000 claims abstract description 4
- 230000005611 electricity Effects 0.000 claims description 5
- 230000006698 induction Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000000717 retained effect Effects 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000001453 impedance spectrum Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 6
- 238000013178 mathematical model Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 238000012502 risk assessment Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种用于电网电磁暂态计算的网络边界等值方法,包括以下步骤:S1、预设谐波阻抗观测母线的步骤,S2、预设等值边界的步骤,S3、获取谐波导纳矩阵基础参数的步骤,S4、基础参数预处理的步骤,S5、形成全网络谐波导纳矩阵的步骤,S6、计算从观测母线向全网络看的谐波阻抗的步骤,S7、对预设等值边界进行Ward等值的步骤,S8、求解网络化简后从观测母线向保留子网看的谐波阻抗的步骤,S9、比较等值前后扫描母线集谐波阻抗的步骤,S10、确定最终等值方案的步骤。本发明的有益效果是:快速得出网络中任意母线向系统看的谐波阻抗值,为电磁暂态计算提供保留原始网络阻抗频谱特性和边界节点间完整电气联系的等值结果,从而提高仿真计算的准确度。
Description
技术领域
本发明涉及电力系统分析计算领域,具体地,涉及一种用于电网电磁暂态计算的网络边界等值方法。
背景技术
电磁暂态计算广泛应用于电力系统的设计、运行和研究分析,是电力系统设备选型、风险评估、决策制定和故障分析的重要技术手段。由于所使用的模型复杂、计算量大,因此电磁暂态计算能够仿真的电网规模较为有限,实际应用中必须对网络边界进行等值处理。
电网电磁暂态计算需要对数微秒至数秒间的电磁暂态过程进行模拟,因此在数赫兹到数千赫兹频带内准确的模拟电网元件以及对外部网络进行等值化简,对于提高电磁暂态计算的准确性极为重要。现有电磁暂态计算程序拥有较为完善和准确的电网元件模型,但由于电网结构的复杂性,难以在等值过程中准确地保留外部电网的原始特性。
目前工程实用的网络边界等值方法是对边界节点逐一进行短路电流计算,再根据外部电网向该节点提供的短路电流计算从该节点向外部网络看的工频等值阻抗。上述方法存在以下缺陷:(1)无法准确描述外部网络阻抗的频谱特性。在外部网络较弱或外部网络含有大容量无功补偿装置时,从边界节点向外部网络看的等值阻抗将呈现出明显的非线性频谱特性。若忽略这些特性,将给电磁暂态计算结果带来较大误差;(2)无法正确描述各边界节点间的电气联系。当保留网络处于原始网络的内部而非边缘时,各边界节点间除了在保留网络内的电气联系外,在外部网络内部亦存在电气联系,传统方法割裂了边界节点在外部网络中的电气联系;(3)对每个边界节点单独处理的方式,可能重复计入外部网络提供的短路电流,使计算出的工频等值阻抗不准确。
鉴于电磁暂态计算在电网设计、运行中的重要地位,研究适用于电网电磁暂态计算的网络边界等值方法,对于提高电网电磁暂态计算的准确性具有重要价值。
发明内容
本发明所要解决的技术问题是提供一种提高仿真计算准确度、用于电网电磁暂态计算的网络边界等值方法。
本发明使用的技术术语:
电磁暂态计算:用数值计算方法对从数微秒至数秒之间的电磁暂态过程进行仿真模拟。
网络边界等值:将网络划分为内部网络、外部网络和边界节点三部分,保持内部网络不变,将外部网络用包含边界节点的简化网络替代的过程。
谐波阻抗:支路或网络的阻抗随频率变化的特性。
Ward等值:一种对线性有源网络进行星—网变换以进行网络化简的方法。
节点导纳矩阵:一种由网络节点电压方程导出,用以描述网络节点连接关系、支路导纳参数的对称稀疏矩阵。
本发明解决上述问题所采用的技术方案是:
一种用于电网电磁暂态计算的网络边界等值方法,包括以下步骤:
S1、预设谐波阻抗观测母线的步骤:谐波阻抗观测母线包含以下三类:(1)内部子网的枢纽站点高压\中压母线;(2)电磁暂态计算中待模拟操作的母线;(3)内部子网中关键发电厂的高压母线。
S2、预设等值边界的步骤;
优选的,该步骤中,以包含所有谐波阻抗观测母线的最小子网为核心内部子网,将核心内部子网与外部网络的边界预设为初始等值边界。
S3、获取谐波导纳矩阵基础参数的步骤:依次读取电网指定运行方式下的母线、线路、变压器、发电机、无功补偿、负荷投入情况及对应的工频参数,包括线路电阻RL、电抗XL、电纳BL、串联补偿XC,变压器电阻RT、电抗XT、非标准变比tk,发电机定子电阻Ra、次暂态电抗Xd”、发电机额定容量Sn,负荷有功功率PL、无功功率QL,感性无功补偿对地电抗XR,容性无功补偿对地电抗XRC。
S4、基础参数预处理的步骤:消去XL<1.0e-5p.u.的短接线支路和相关短接线母线,对有功功率小于0的负荷将有功功率置0,去除所有发电机、负荷、线路、变压器及无功补偿未使用的闲置母线,对所有节点间和节点对地的并联支路进行合并,并更新线路、变压器、发电机、无功补偿、负荷工频参数。
S5、形成全网络谐波导纳矩阵的步骤:在fstart至fend范围内,按照设定fstep步长依次形成谐波导纳矩阵,扫描频率f=fstart+hfstep(h=0,1,…(fend-fstart)/fstep),其中,fstart为谐波阻抗扫描起始频率,fend为谐波阻抗扫描结束频率,fstep为扫描频率变化的步长,按照设定fstep步长依次建立元件谐波计算模型,基于元件谐波计算模型,形成全网络的谐波导纳矩阵簇;
优选的,各元件的谐波计算模型为:
对于发电机,串联电阻为串联电抗为其中SB=100MVA为系统基准容量,fb为系统额定频率;
对于无功补偿,并联电感支路感抗为并联电容支路容抗为
对于负荷,并联电阻支路电阻为并联电抗支路感抗并联电容支路容抗
对于输电线路,i,j两节点间串联阻抗为
,
i、j侧并联导纳为
对于变压器,i,j两节点间串联阻抗为i侧并联导纳为j侧并联导纳为其中
基于上述元件谐波计算模型,形成全网络的谐波导纳矩阵簇。
进一步的,其中,非对角元与i、j节点间支路导纳为对角元等于该行所有非对角元的相反数与该节点对地导纳之和,即
,其中,YTeq=YTeq1(or YTeq2)。
S6、计算从观测母线向全网络看的谐波阻抗的步骤:针对每一个扫描频率下的谐波导纳矩阵依次设置谐波阻抗观测母线的注入电流Iδ为1,利用三角分解法求解下列方程,
求得Uδ即为该母线在该指定频率下向系统看的阻抗,上式中,n为电网的母线数,δ为观测母线在导纳矩阵中的序列号。再依次求解其他谐波阻抗观测母线的观测阻抗,直到所有指定的频率都计算完毕为止。最终汇集所有观测母线在所计算频带内的谐波阻抗值,获取待扫描母线集的谐波阻抗阵列Zh,其结构如下:
其中,k为观测母线个数,m为扫描频率点数。
S7、对预设等值边界进行Ward等值的步骤:将网络中的节点集合划分为内部系统节点子集{I}、边界节点子集{B}和外部网络节点子集{E},然后将整个网络的工频节点导纳矩阵按节点集合的划分写成分块矩阵
消去外部系统的节点子集,获取消去外部系统后在边界节点附加的节点导纳矩阵YEQ,经等值处理后得到简化系统。基于节点导纳矩阵中的对应关系,得到化简后等值边界附加支路的工频参数。
S8、求解网络化简后从观测母线向保留子网看的谐波阻抗的步骤:按照前述形成网络谐波导纳矩阵的步骤,形成含保留子网和等值化简网络的谐波导纳矩阵,再计算从观测母线向保留子网看的谐波阻抗,最终形成网络等值化简后待扫描母线集的谐波阻抗阵列Z′h。
S9、比较等值前后扫描母线集谐波阻抗的步骤:将矩阵Zh和Z′h中所有元素均转换成极坐标形式,获得模矩阵Mh、M′h和相角矩阵Ph、P′h。计算等值前后的差异矩阵
其中,
再计算出综合差异指标y;
优选的,差异指标
确定最终等值方案的步骤:当综合差异指标大于阈值ymax,则返回修改预设等值边界,将等值边界向外部网络延拓一级,重复上述步骤,直到综合差异指标满足阈值要求,优选的,所述的综合差异指标阈值ymax=0.05。
利用本发明提出的方法,可以获知每一母线向系统看的谐波阻抗,为谐波、故障分析提供参考,在对网络边界进行等值的过程中不仅保留了外部网络的谐波阻抗特性,也保留了等值边界各节点间的互阻抗联系,从而使等值结果更多的留存了原始网络的频率阻抗特性,提高了网络等值和电磁暂态计算的准确性。
为在网络边界等值过程中准确保留外部网络的谐波阻抗特性和边界节点间的电气联系,该方法通过读取电网的线路、变压器、发电机、无功补偿及负荷参数,形成指定运行方式下网络的谐波导纳矩阵,基于等值前网络谐波导纳矩阵对内部网络指定节点进行谐波阻抗扫描,对预设边界节点进行Ward等值获取包含等值化简网络的保留子网,再基于保留子网谐波导纳矩阵对指定节点进行谐波阻抗扫描,分析比较等值前后指定节点的谐波阻抗,当二者不满足误差条件时,返回修改预设边界、扩大保留子网范围,再对边界进行Ward等值和扫描保留子网指定节点谐波阻抗,直至等值前后指定节点谐波阻抗满足误差条件为止。
综上,本发明的有益效果是:
电磁暂态计算在电力系统的设计、运行和分析研究等领域中具有举足轻重的作用,是电力系统设备选型、风险评估、决策制定和故障分析的重要技术手段。由于所使用的模型复杂、计算量大,因此电磁暂态计算能够仿真的电网规模较为有限,实际应用中必须对网络边界进行等值处理。传统的电磁暂态计算边界等值方法是对等值边界节点逐一进行短路电流计算,再根据外部电网向该节点提供的短路电流计算从该节点向外部网络看的工频等值阻抗,忽略了外部网络阻抗的频谱特性,也无法正确描述各边界节点间的电气联系,从而对电磁暂态计算的准确度造成影响。
本发明基于等值前后的网络谐波导纳矩阵,对内部网络指定母线进行谐波阻抗扫描,可获取指定母线向电网看的谐波阻抗特性,同时也使等值子网保留了原始网络的阻抗频谱特性,解决了弱联系电网的边界等值问题。
本发明对边界节点进行Ward等值,获取包含等值化简网络的保留子网,从而完整的保留了等值边界节点间的电气联系,解决了环网内部子网络的边界等值问题。
本发明对等值前后内部网络指定母线的谐波阻抗模值、相位进行比较,从而可以不依赖于数值仿真结果来验证等值的准确性。
本发明提出的方法解决了在等值子网中保留原始网络阻抗频谱特性的问题和环网内部子网的边界等值问题,并且提供了一种不依赖于数值仿真结果来验证等值准确性的途径。利用本发明提出的方法,可以快速得出网络中任意母线向系统看的谐波阻抗值,可以为电磁暂态计算提供保留原始网络阻抗频谱特性和边界节点间完整电气联系的等值结果,从而提高仿真计算的准确度。
附图说明
图1是发电机支路的数学模型;
图2是无功补偿支路的数学模型;
图3是负荷支路的数学模型;
图4是输电线路的数学模型;
图5是变压器的数学模型;
图6是不同频率下网络的谐波导纳矩阵簇;
图7是Ward等值后的边界节点及其支路连接关系;
图8是实施例所使用的电网结构图。
具体实施方式
下面结合实施例及附图,对本发明作进一步地的详细说明,但本发明的实施方式不限于此。
实施例:
一种用于电网电磁暂态计算的网络边界等值方法,包括以下步骤:
S1、预设谐波阻抗观测母线的步骤:指定电网数据的工作目录和潮流方式号,设置频率扫描的起始频率fstart=50Hz,终止频率fend=2000Hz,扫描步长fstep=10Hz。设置谐波阻抗观测母线,如图8所示的母线YUL、LANCJ和BANGD为谐波阻抗观测母线。
S2、预设等值边界的步骤:设置初始等值边界,如图8所示的线路LANCJ-BAT。
S3、获取谐波导纳矩阵基础参数的步骤:依次读取PSASP程序基础数据库中与指定潮流方式号对应的母线、线路、变压器、发电机、无功补偿、负荷投入情况及工频参数,包括线路电阻RL、电抗XL、电纳BL、串连补偿XC,变压器电阻RT、电抗XT、非标准变比tk,发电机定子电阻Ra、次暂态电抗Xd”、发电机额定容量Sn,负荷有功功率PL、无功功率QL,感性无功补偿对地电抗XR,容性无功补偿对地电抗XRC。
S4、基础参数预处理的步骤:消去XL<1.0e-5p.u.的短接线支路和相关附加节点,当负荷PL<0时,将PL置0,去除所有发电机、负荷、线路、变压器及无功补偿未使用的闲置母线,对所有节点间和节点对地的并联支路进行合并,并更新线路、变压器、发电机、无功补偿、负荷工频参数。
S5、形成全网络谐波导纳矩阵的步骤:在fstart至fend范围内,按照设定fstep步长依次形成谐波导纳矩阵,扫描频率f=fstart+hfstep(h=0,1,…195),基于元件谐波计算模型,形成全网络的谐波导纳矩阵簇。
本实施例各元件的谐波计算模型如图1—图5所示。
对于发电机,如图1所示,串联电阻为串联电抗为其中SB=100MVA为系统基准容量,fb=50Hz为系统额定频率;
对于无功补偿,如图2所示,并联电感支路感抗为并联电容支路容抗为
对于负荷,如图3所示,并联电阻支路电阻为并联电抗支路感抗并联电容支路容抗
对于输电线路,如图4所示,i,j两节点间串联阻抗为
,i、j侧并联导纳为
对于变压器,如图5所示,i,j两节点间串联阻抗为i侧并联导纳为j侧并联导纳为其中
基于元件谐波计算模型,形成全网络的谐波导纳矩阵簇,如图6所示。其中,非对角元与i、j节点间支路导纳为对角元等于该行所有非对角元的相反数与该节点对地导纳之和,即
YLeq+YTeq+sum(Yij),其中,YTeq=YTeq1(or YTeq2)。
S6、计算从观测母线向全网络看的谐波阻抗的步骤:针对每一个扫描频率下的谐波导纳矩阵依次设置谐波阻抗观测母线YUL、LANCJ、BANGD的注入电流Iδ为1,利用三角分解法求解下列方程,
求得Uδ即为该母线在该指定频率下向系统看的阻抗,上式中,n为电网的母线数,δ为观测母线在导纳矩阵中的序列号。再依次求解其他谐波阻抗观测母线的观测阻抗,直到所有指定的频率都计算完毕为止。最终汇集所有待计算观测母线在所计算频带内的谐波阻抗值,获取待扫描母线集的谐波阻抗阵列Zh,其结构如下:
其中,m为扫描频率点数,m=195。
S7、对预设等值边界进行Ward等值的步骤:将网络中的节点集合划分为内部系统节点子集{I}、边界节点子集{B}和外部网络节点子集{E},然后将整个系统的工频节点导纳矩阵按节点集合的划分写成分块矩阵
消去外部系统的节点子集,获取消去外部系统后在边界节点附加的节点导纳矩阵YEQ,经等值处理后的简化系统如图7所示。基于节点导纳矩阵中的对应关系,得到化简后等值边界附加支路的工频参数。
S8、求解网络化简后从观测母线向保留子网看的谐波阻抗的步骤:按照形成网络谐波导纳矩阵的步骤,形成含保留子网和等值化简网络的谐波导纳矩阵,再计算从观测母线向保留子网看的谐波阻抗,最终形成网络等值化简后待扫描母线集的谐波阻抗阵列Z′h。
S9、比较等值前后扫描母线集谐波阻抗的步骤:将矩阵Zh和Z′h中所有元素均转换成极坐标形式,获得模矩阵Mh、M′h和相角矩阵Ph、P′h。计算等值前后的差异矩阵
其中,
再计算出综合差异指标y;
指标。
S10、确定最终等值方案的步骤:设定综合差异指标y的阈值ymax为0.05,当综合差异指标y>0.05,则返回S2修改预设等值边界,将等值边界向外部网络延拓一级,如图8所示的等值边界2、等值边界3,重复步骤S2~S9,直到综合差异指标y<0.05。
如上所述,可较好的实现本发明。
Claims (6)
1.一种用于电网电磁暂态计算的网络边界等值方法,其特征在于,包括以下步骤:
S1、预设谐波阻抗观测母线的步骤:将内部子网的枢纽站点高压\中压母线、电磁暂态计算中待模拟操作的母线、内部子网中关键发电厂的高压母线设置为谐波阻抗观测母线;
S2、预设等值边界的步骤;
S3、获取谐波导纳矩阵基础参数的步骤:依次读取电网指定运行方式下的线路电阻RL、电抗XL、电纳BL、串连补偿XC,变压器电阻RT、电抗XT、非标准变比tk,发电机定子电阻Ra、次暂态电抗Xd”、发电机额定容量Sn,负荷有功功率PL、无功功率QL,感性无功补偿对地电抗XR,容性无功补偿对地电抗XRC;
S4、基础参数预处理的步骤:消去XL<1.0e-5p.u.的短接线支路和相关附加节点,对有功功率小于0的负荷将有功功率置0,去除所有发电机、负荷、线路、变压器及无功补偿未使用的闲置母线,对所有节点间和节点对地的并联支路进行合并,并更新线路、变压器、发电机、无功补偿、负荷工频参数;
S5、形成全网络谐波导纳矩阵的步骤:基于设备工频参数在fstart至fend范围内,按照设定fstep步长依次建立元件谐波计算模型,基于元件谐波计算模型,形成全网络的谐波导纳矩阵簇;fstart为谐波阻抗扫描起始频率,fend为谐波阻抗扫描结束频率,fstep为扫描频率变化的步长;
S6、计算从观测母线向全网络看的谐波阻抗的步骤:针对每一个扫描频率下的谐波导纳矩阵依次设置谐波阻抗观测母线的注入电流Iδ为1,利用三角分解法求解下列方程,
求得Uδ即为该母线在该指定频率下向系统看的阻抗,上式中,n为电网的母线数,δ为观测母线在导纳矩阵中的序列号;再依次求解其他谐波阻抗观测母线的观测阻抗,直到所有指定的频率都计算完毕为止;最终汇集所有待计算观测母线在所计算频带内的谐波阻抗值,获取待扫描母线集的谐波阻抗阵列Zh:
其中,k为观测母线个数,m为扫描频率点数;
S7、对预设等值边界进行Ward等值的步骤:将网络中的节点集合划分为内部系统节点子集{I}、边界节点子集{B}和外部网络节点子集{E},然后将整个系统的工频节点导纳矩阵按节点集合的划分写成分块矩阵
消去外部系统的节点子集,获取消去外部系统的节点子集后在边界节点附加的节点导纳矩阵YEQ,基于节点导纳矩阵中的对应关系,得到化简后等值边界附加支路的工频参数;
S8、求解网络化简后从观测母线向保留子网看的谐波阻抗的步骤:形成含保留子网和等值化简网络的谐波导纳矩阵,再计算从观测母线向保留子网看的谐波阻抗,最终形成网络等值化简后待扫描母线集的谐波阻抗阵列Z′h;
S9、比较等值前后扫描母线集谐波阻抗的步骤:将矩阵Zh和Z′h中所有元素均转换成极坐标形式,获得模矩阵Mh、M′h和相角矩阵Ph、P′h;
计算等值前后的差异矩阵
其中, 再计算出综合差异指标y;
S10、确定最终等值方案的步骤:当综合差异指标y大于阈值ymax,则返回修改预设等值边界,将等值边界向外部网络延拓一级,重复上述步骤,直到综合差异指标满足阈值要求。
2.根据权利要求1所述的一种用于电网电磁暂态计算的网络边界等值方法,其特征在于,所述的元件谐波计算模型为:
对于发电机,串联电阻为串联电抗为其中SB=100MVA为系统基准容量,fb为系统额定频率;对于无功补偿,并联电感支路感抗为并联电容支路容抗为对于负荷,并联电阻支路电阻为并联电抗支路感抗并联电容支路容抗对于输电线路,i,j两节点间串联阻抗为 i、j侧并联导纳为 对于变压器,i,j两节点间串联阻抗为i侧并联导纳为j侧并联导纳为其中扫描频率f=fstart+hfstep(h=0,1,…(fend-fstart)/fstep),fstart为谐波阻抗扫描起始频率,fend为谐波阻抗扫描结束频率,fstep为扫描频率变化的步长。
3.根据权利要求2所述的一种用于电网电磁暂态计算的网络边界等值方法,其特征在于,对于所述的谐波导纳矩阵簇,其中,非对角元与i、j节点间支路导纳为对角元等于该行所有非对角元的相反数与该节点对地导纳之和,即
其中,YTeq=YTeq1(or YTeq2)。
4.根据权利要求1或2所述的一种用于电网电磁暂态计算的网络边界等值方法,其特征在于,所述的综合差异指标
5.根据权利要求1或2所述的一种用于电网电磁暂态计算的网络边界等值方法,其特征在于,所述的综合差异指标阈值ymax=0.05。
6.根据权利要求1或2所述的一种用于电网电磁暂态计算的网络边界等值方法,其特征在于,所述的步骤S2是以包含所有谐波阻抗观测母线的最小子网为核心内部子网,将核心内部子网与外部网络的边界预设为初始等值边界。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510402799.1A CN105160058B (zh) | 2015-07-10 | 2015-07-10 | 一种用于电网电磁暂态计算的网络边界等值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510402799.1A CN105160058B (zh) | 2015-07-10 | 2015-07-10 | 一种用于电网电磁暂态计算的网络边界等值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105160058A CN105160058A (zh) | 2015-12-16 |
CN105160058B true CN105160058B (zh) | 2018-12-25 |
Family
ID=54800914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510402799.1A Active CN105160058B (zh) | 2015-07-10 | 2015-07-10 | 一种用于电网电磁暂态计算的网络边界等值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105160058B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105915303B (zh) * | 2016-04-18 | 2018-03-09 | 华北电力大学(保定) | 一种基于支路追加法的配电网电力线载波通信信道建模法 |
CN109254187A (zh) * | 2017-07-14 | 2019-01-22 | 中国电力工程顾问集团华东电力设计院有限公司 | 一种确定变电站系统短路电流的系统及方法 |
CN110188386B (zh) * | 2019-04-26 | 2023-02-28 | 南方电网科学研究院有限责任公司 | 交流系统谐波阻抗分区方法、装置、设备及存储介质 |
CN112134288B (zh) * | 2020-08-28 | 2022-07-01 | 国网江苏省电力有限公司滨海县供电分公司 | 基于前/后向扫描谐波潮流的谐波污染配电网络重构方法 |
CN116822436B (zh) * | 2023-06-30 | 2024-02-27 | 四川大学 | 一种直流输电送端交流系统的振荡风险灵敏度分析方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006065859A2 (en) * | 2004-12-14 | 2006-06-22 | Brian Costa | Dynamic load control in a customer distribution network |
CN101957872A (zh) * | 2010-07-09 | 2011-01-26 | 南方电网技术研究中心 | 一种交直流电力系统的混合实时仿真方法 |
CN102104252A (zh) * | 2011-02-28 | 2011-06-22 | 浙江大学 | 一种适用于电磁暂态分析的电力系统动态等值方法 |
CN103166237A (zh) * | 2013-02-28 | 2013-06-19 | 四川电力科学研究院 | 一种弱联系受端小电网联络线高抗补偿方法 |
CN103559347A (zh) * | 2013-10-31 | 2014-02-05 | 浙江大学 | 一种大规模交直流电力系统电磁暂态仿真模型的构建方法 |
CN104009475A (zh) * | 2014-06-19 | 2014-08-27 | 重庆大学 | 考虑等值元件全面性与参数物理约束的非拓扑静态等值法 |
CN104698346A (zh) * | 2015-03-31 | 2015-06-10 | 国网内蒙古东部电力有限公司通辽供电公司 | 含源配电网故障分析方法及装置 |
-
2015
- 2015-07-10 CN CN201510402799.1A patent/CN105160058B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006065859A2 (en) * | 2004-12-14 | 2006-06-22 | Brian Costa | Dynamic load control in a customer distribution network |
CN101957872A (zh) * | 2010-07-09 | 2011-01-26 | 南方电网技术研究中心 | 一种交直流电力系统的混合实时仿真方法 |
CN102104252A (zh) * | 2011-02-28 | 2011-06-22 | 浙江大学 | 一种适用于电磁暂态分析的电力系统动态等值方法 |
CN103166237A (zh) * | 2013-02-28 | 2013-06-19 | 四川电力科学研究院 | 一种弱联系受端小电网联络线高抗补偿方法 |
CN103559347A (zh) * | 2013-10-31 | 2014-02-05 | 浙江大学 | 一种大规模交直流电力系统电磁暂态仿真模型的构建方法 |
CN104009475A (zh) * | 2014-06-19 | 2014-08-27 | 重庆大学 | 考虑等值元件全面性与参数物理约束的非拓扑静态等值法 |
CN104009475B (zh) * | 2014-06-19 | 2017-01-18 | 重庆大学 | 考虑等值元件全面性与参数物理约束的非拓扑静态等值法 |
CN104698346A (zh) * | 2015-03-31 | 2015-06-10 | 国网内蒙古东部电力有限公司通辽供电公司 | 含源配电网故障分析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105160058A (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105160058B (zh) | 一种用于电网电磁暂态计算的网络边界等值方法 | |
CN107453357B (zh) | 一种基于分层求解的配电网状态估计方法 | |
CN103728535B (zh) | 一种基于小波变换暂态能量谱的特高压直流输电线路故障测距方法 | |
CN101505061B (zh) | 一种基于时域仿真的戴维南等值参数跟踪的计算方法 | |
CN105529711B (zh) | 一种基于bpa数据的交流系统谐波阻抗扫描方法 | |
CN106250660A (zh) | 一种海上风电场谐振分析方法 | |
CN103018534B (zh) | 确定谐波电压的方法及系统 | |
CN110221237A (zh) | 一种电容式电压互感器谐波电压测量方法和装置 | |
CN105676022A (zh) | 一种长线路风电场并网谐振信息提取方法 | |
CN108828318B (zh) | 一种提取级联型隔离变压器寄生电容的方法 | |
CN103605856B (zh) | 基于分数阶线路模型的输电线路参数估计方法 | |
CN105823929B (zh) | 基于双端异步数据的长距离同塔四回输电线路零序参数精确测量方法 | |
Mendonça et al. | Wind farm and system modelling evaluation in harmonic propagation studies | |
CN104865441B (zh) | 一种串扰频谱计算方法 | |
CN112072692A (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
CN111444592B (zh) | 一种变压器宽频导纳模型及建立方法 | |
CN102709911A (zh) | 一种混合仿真谐波特性的接口设计方法 | |
CN114167133B (zh) | 一种输电网内新建站点的谐波电压评估方法和装置 | |
CN112966364A (zh) | 一种用于特征值计算的光伏电站等值模型建模方法及装置 | |
Bigdeli et al. | Identification of transient model parameters of transformer using genetic algorithm | |
CN103713196B (zh) | 一种交流特高压同塔双回线路相自导纳、相自阻抗参数测量方法 | |
CN116299138A (zh) | 一种10kV电压互感器的亚纳秒脉冲响应特性检测方法 | |
CN107436995B (zh) | 考虑外网对地支路和灵敏度信息的等值三相短路计算方法 | |
Ramzan et al. | Accurate harmonic analysis of distribution systems | |
Eguia et al. | Characterization of network harmonic impedance for grid connection studies of renewable plants |
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 |