CN109359266B - 感应电动机暂态响应的求解方法、系统、设备及存储介质 - Google Patents
感应电动机暂态响应的求解方法、系统、设备及存储介质 Download PDFInfo
- Publication number
- CN109359266B CN109359266B CN201811252348.4A CN201811252348A CN109359266B CN 109359266 B CN109359266 B CN 109359266B CN 201811252348 A CN201811252348 A CN 201811252348A CN 109359266 B CN109359266 B CN 109359266B
- Authority
- CN
- China
- Prior art keywords
- induction motor
- negative sequence
- equation
- sequence
- power consumption
- 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
- 230000006698 induction Effects 0.000 title claims abstract description 244
- 230000001052 transient effect Effects 0.000 title claims abstract description 169
- 230000004044 response Effects 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000004364 calculation method Methods 0.000 claims abstract description 90
- 238000004088 simulation Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 4
- 230000000875 corresponding effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 12
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 239000000126 substance Substances 0.000 description 3
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
-
- 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种感应电动机暂态响应的求解方法,包括:获取感应电动机的参数数据;根据参数数据计算出近似负序阻抗Zin2;再结合电网系统的电网数据,确定出负序电磁转矩Te‑,负序有功消耗Pe‑及负序无功消耗Qe‑;联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出暂态响应;不对称故障期间三阶机电暂态方程中的总电磁转矩Te计入负序电磁转矩Te‑,有功消耗计算方程计入负序有功消耗Pe‑,无功消耗计算方程计入负序无功消耗Qe‑。应用本申请的方案,有效地提高了不对称扰动下计算出的暂态响应的准确性且计算简便。本申请还公开了一种感应电动机暂态响应的求解系统、设备及存储介质,具有相应效果。
Description
技术领域
本发明涉及电网分析技术领域,特别是涉及一种感应电动机暂态响应的求解方法、系统、设备及存储介质。
背景技术
负荷特性对电力系统的稳定性有着重要影响,感应电动机在电力系统的负荷中的占比高达60%~70%,是影响电力系统稳定性的重要因素。
通常采用数字仿真分析大电网的稳定性,因此,数字仿真能否准确模拟感应电动机的暂态响应,直接影响电力系统稳定性分析结果的可信度,进而影响相应的控制策略的正确性。
目前,分析时广泛采用的是感应电动机的三阶机电暂态模型。该模型的适用条件是电动机母线电压为基波正序,当电力系统发生对称故障时,模型的计算结果较为准确,但发生不对称故障时,由于该模型忽略了负序分量的影响,但此时感应电动机的负序定子电流一般较大,忽略负序电流会影响暂态无功和电压响应的计算精度,此外,感应电动机转子上会出现制动性质的负序转矩,忽略这一转矩可能影响转子滑差及电动机有功功耗的计算精度。而感应电动机的五阶电磁暂态模型虽然能够较为准确地反映电动机在对称和不对称扰动下的真实暂态响应,但由于五阶电磁暂态模型要求电力网络也用微分方程进行描述,这会大幅增加计算负担,无法应用在大电网的稳定性分析中。
综上所述,如何在不大幅增加计算负担的前提下,有效地提高不对称扰动下计算出的感应电动机的暂态响应的准确性,是目前本领域技术人员急需解决的技术问题。
发明内容
本发明的目的是提供一种感应电动机暂态响应的求解方法、系统、设备及存储介质,在不大幅增加计算负担的前提下,有效地提高不对称扰动下计算出的感应电动机的暂态响应的准确性。
为解决上述技术问题,本发明提供如下技术方案:
一种感应电动机暂态响应的求解方法,包括:
获取感应电动机的参数数据;
根据所述参数数据计算出所述感应电动机的近似负序阻抗Zin2;
联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出所述感应电动机的暂态响应;
其中,不对称故障期间所述感应电动机三阶机电暂态方程中的总电磁转矩Te为计入所述负序电磁转矩Te-的电磁转矩,所述有功消耗计算方程为计入负序有功消耗Pe-的方程,所述无功消耗计算方程为计入负序无功消耗Qe-的方程。
通过所述电网系统的结构建立负序网络方程;
其中,所述Rs为所述参数数据中的所述感应电动机的定子电阻,所述X′为所述感应电动机的转子短路电抗,满足:X′=Xs+XrXm/(Xr+Xm);所述Xr为所述参数数据中的所述感应电动机的转子电抗,所述Xs为所述参数数据中的所述感应电动机的定子电抗,所述Xm为所述参数数据中的所述感应电动机的励磁电抗。
优选的,所述根据所述参数数据计算出所述感应电动机的近似负序阻抗Zin2,包括:
通过忽略滑差变化构建所述感应电动机的负序稳态等值电路;
根据所述参数数据以及所述负序稳态等值电路,计算出所述感应电动机的负序等值阻抗,并将所述负序等值阻抗作为计算出的所述感应电动机的近似负序阻抗Zin2。
优选的,计算出的所述近似负序阻抗Zin2表示为:
其中,所述Rs为所述参数数据中的所述感应电动机的定子电阻,所述Xr为所述参数数据中的所述感应电动机的转子电抗,所述Xm为所述参数数据中的所述感应电动机的励磁电抗,所述Xs为所述参数数据中的所述感应电动机的定子电抗,所述Rr为所述参数数据中的所述感应电动机的转子电阻。
优选的,确定出的所述负序电磁转矩Te-表示为:
其中,real表示取相量的实部,*表示相量的共轭;
确定出的所述负序有功消耗Pe-表示为:
确定出的所述负序无功消耗Qe-表示为:
其中,imag表示取相量的虚部。
优选的,计入所述负序电磁转矩Te-的所述总电磁转矩Te表示为:
计入所述负序有功消耗Pe-的所述有功消耗计算方程表示为:
计入所述负序无功消耗Qe-的所述无功消耗计算方程表示为:
一种感应电动机暂态响应的求解系统,包括:
参数数据获取模块,用于获取感应电动机的参数数据;
负序阻抗计算模块,用于根据所述参数数据计算出所述感应电动机的近似负序阻抗Zin2;
暂态响应计算模块,用于联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出所述感应电动机的暂态响应;其中,不对称故障期间所述感应电动机三阶机电暂态方程中的总电磁转矩Te为计入所述负序电磁转矩Te-的电磁转矩,所述有功消耗计算方程为计入负序有功消耗Pe-的方程,所述无功消耗计算方程为计入负序无功消耗Qe-的方程。
一种感应电动机暂态响应的求解设备,包括:
存储器,用于存储暂态响应求解程序;
处理器,用于执行所述暂态响应求解程序以实现上述任一项所述的感应电动机暂态响应的求解方法的步骤。
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述任一项所述的感应电动机暂态响应的求解方法的步骤。
应用本发明实施例所提供的技术方案,在三阶机电暂态模型的基础上考虑了负序分量。具体的,不对称故障期间感应电动机三阶机电暂态方程中的总电磁转矩Te为计入负序电磁转矩Te-的电磁转矩,有功消耗计算方程为计入负序有功消耗Pe-的方程,无功消耗计算方程为计入负序无功消耗Qe-的方程。也就使得计算出的暂态响应更为准确,并且,由于仍是在三阶机电暂态模型的基础上进行暂态分量的计算,使得不会提高计算的复杂度。而负序分量获得的过程也较为简便,具体的,只需要获取感应电动机的参数数据以及感应电动机所处的电网系统的电网数据,便可以计算出感应电动机的负序定子电压相量负序定子电流相量/>以及负序转子暂态电势相量/>进而计算出感应电动机的负序电磁转矩Te-,负序有功消耗Pe-以及负序无功消耗Qe-;这一过程中并不需要将电力网络用微分方程的形式进行表述,仅需获得电网系统的电网数据即可,使得计算简便。因此,本申请的方案,在不大幅增加计算负担的前提下,有效地提高了不对称扰动下计算出的感应电动机的暂态响应的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中一种感应电动机暂态响应的求解方法的实施流程图;
图2为本发明一种具体实施方式中感应电动机的近似负序稳态等值电路的结构示意图;
图3为本发明一种具体实施方式中电网系统的结构示意图;
图4为对称电压跌落下IEEE-2型感应电动机的现有的三阶机电暂态模型与电磁暂态仿真的比较示意图;
图5为对称电压跌落下IEEE-5型感应电动机的现有的三阶机电暂态模型与电磁暂态仿真的比较示意图;
图6为不对称电压跌落下IEEE-2型感应电动机的现有的三阶机电暂态模型与电磁暂态仿真的比较示意图;
图7为不对称电压跌落下IEEE-5型感应电动机的现有的三阶机电暂态模型与电磁暂态仿真的比较示意图;
图8为不对称电压跌落下IEEE-2型感应电动机的本申请的一种具体方案与电磁暂态仿真的比较示意图;
图9为不对称电压跌落下IEEE-5型感应电动机的本申请的一种具体方案与电磁暂态仿真的比较示意图;
图10为本发明中一种感应电动机暂态响应的求解系统的结构示意图;
图11为本发明中一种感应电动机暂态响应的求解设备的结构示意图。
具体实施方式
本发明的核心是提供一种不对称故障下考虑负序分量影响的感应电动机暂态响应的求解方法,能够在不大幅增加计算负担的前提下,有效地提高计算出的感应电动机的暂态响应的准确性。
为了使本技术领域的人员更好地理解本发明方案,下面结合附图和具体实施方式对本发明作进一步的详细说明。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参考图1,图1为本发明中一种感应电动机暂态响应的求解方法的实施流程图,该方法包括以下步骤:
步骤S101:获取感应电动机的参数数据。
该步骤中的参数数据表示的是感应电动机的基础数据,并且转子侧参数已经归算至定子侧。具体的,可以包括感应电动机的定子电阻Rs,定子电抗Xs,转子电抗Xr,励磁电抗Xm,转子电阻Rr,惯性时间常数Tj以及负载转矩Tm。
步骤S102:根据参数数据计算出感应电动机的近似负序阻抗Zin2。
具体的,可以结合感应电动机的参数数据计算出感应电动机的近似负序阻抗Zin2。
在本发明的一种具体实施方式中,步骤S102可以具体包括以下两个步骤:
步骤一:通过忽略滑差变化构建感应电动机的负序稳态等值电路;
步骤二:根据参数数据以及负序稳态等值电路,计算出感应电动机的负序等值阻抗,并将负序等值阻抗作为计算出的感应电动机的近似负序阻抗Zin2。
具体实施时,申请人考虑到由于滑差较小,并且其变化对负序阻抗的计算值影响较小,因此忽略滑差变化构建负序稳态等值电路后,可以方便地计算出近似负序阻抗Zin2。
可参阅图2,为该种实施方式中的感应电动机的近似负序稳态等值电路的结构示意图,通过该近似负序稳态等值电路可以近似替代感应电动机,从而结合感应电动机的参数数据计算出相应的负序等值阻抗,并作为计算出的感应电动机的近似负序阻抗Zin2,可以表示为:
其中的Rs即为参数数据中的定子电阻,Xr为参数数据中的转子电抗,Xm为参数数据中的励磁电抗,Xs为参数数据中的定子电抗,Rr为参数数据中的转子电阻。
在计算出感应电动机的近似负序阻抗Zin2之后,可以执行步骤S103的操作。
申请人考虑到在发生不对称故障时,感应电动机的定子端出现负序电压在电动机转子中激发出负序暂态电势/>产生相应的负序定子电流/>并且,考虑到/>的数值较小且变化率较大,可以认为/>在故障发生时刻即达到其暂态终值,在故障消失时刻即降低为零。因此,故障期间的/>可用感应电机的稳态等值电路来计算。
具体的,可以通过负序网络方程确定出负序定子电压相量负序网络方程可以表示为:Y2V2=I2。其中的Y2为计入了近似负序阻抗Zin2的负序节点导纳矩阵,V2为节点负序电压相量,I2为节点负序注入电流相量。之后,再结合近似负序阻抗Zin2计算出负序定子电流相量/>以及负序转子暂态电势相量/>
其中的Rs即为定子电阻,X′为感应电动机的转子短路电抗,满足:X′=Xs+XrXm/(Xr+Xm);Xr为转子电抗,Xm为励磁电抗,Xs为定子电抗。
具体的,确定出的负序电磁转矩Te-可以表示为:
确定出的负序有功消耗Pe-可以表示为:
确定出的负序无功消耗Qe-可以表示为:
其中,real表示取相量的实部,*表示相量的共轭,imag表示取相量的虚部。
步骤S105:联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出感应电动机的暂态响应;
其中,不对称故障期间感应电动机三阶机电暂态方程中的总电磁转矩Te为计入负序电磁转矩Te-的电磁转矩,有功消耗计算方程为计入负序有功消耗Pe-的方程,无功消耗计算方程为计入负序无功消耗Qe-的方程。
本申请的方案中,在计算暂态响应时考虑进了负序分量的影响,但仍是建立在感应电动机三阶机电暂态模型的基础上进行暂态响应的计算,因此降低了计算的负担。
具体的,正序电力网络方程的一般形式可以为:YV=I,其中,Y为节点导纳矩阵,V为节点电压相量,I为节点注入电流相量。
感应电动机定子电压方程可以表示为:其中,/>为感应电动机的转子正序暂态电势相量,/>为感应电动机的正序定子电流相量,/>为感应电动机的正序定子电压相量,X′为感应电动机的转子短路电抗,Rs为定子电阻。
感应电动机三阶机电暂态方程可以表示为:
其中,E′d和E′q为的d轴和q轴分量;T′0表示感应电动机定子开路、转子回路暂态时间常数,可由感应电动机的参数数据获得:T′0=(Xr+Xm)/(ωsRr);ωs为同步转速,ωs=2πf,f为电网的额定频率,我国为50Hz。X为感应电动机的转子开路电抗,X=Xs+Xm;s为感应电动机的滑差;Ids和Iqs为/>的d轴和q轴分量。
计入负序有功消耗Pe-的有功消耗计算方程可以表示为:
计入负序无功消耗Qe-的无功消耗计算方程可以表示为:
联立上述正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程之后,可以计算出感应电动机的暂态响应。具体的,可以参照现有的机电暂态仿真方法进行求解,即先用梯形积分公式把微分方程差分化,再联立所有方程,用牛顿法求解,典型的计算步长为0.01s。
应用本发明实施例所提供的感应电动机暂态响应的求解方法,在三阶机电暂态模型的基础上考虑了负序分量。具体的,不对称故障期间感应电动机三阶机电暂态方程中的总电磁转矩Te为计入负序电磁转矩Te-的电磁转矩,有功消耗计算方程为计入负序有功消耗Pe-的方程,无功消耗计算方程为计入负序无功消耗Qe-的方程。也就使得计算出的暂态响应更为准确,并且,由于仍是在三阶机电暂态模型的基础上进行暂态分量的计算,使得不会提高计算的复杂度。而负序分量获得的过程也较为简便,具体的,只需要获取感应电动机的参数数据以及感应电动机所处的电网系统的电网数据,便可以计算出感应电动机的负序定子电压相量负序定子电流相量/>以及负序转子暂态电势相量/>进而计算出感应电动机的负序电磁转矩Te-,负序有功消耗Pe-以及负序无功消耗Qe-;这一过程中并不需要将电力网络用微分方程的形式进行表述,仅需获得电网系统的电网数据即可,使得计算简便。因此,本申请的方案,在不大幅增加计算负担的前提下,有效地提高了不对称扰动下计算出的感应电动机的暂态响应的准确性。
在本发明的一种具体实施方式中,步骤S103可以具体包括:
通过电网系统的结构建立负序网络方程;
可以根据电网系统的具体结构建立负序网络方程,例如根据图3的电网系统,建立的负序网络方程可以为:
需要说明的是,本实施例中便于描述以图3的两节点系统为例,在实际应用中,电网系统可以为其他结构,相适应地进行负序网络方程的建立即可,并不影响本发明的实施。
在图3的实施方式中,假定电网数据中的电源电压Eeq为115kV,输电回路等值阻抗Zeq=(3.4+j8.0425)Ω,f=50Hz,负序网络方程中的可以根据故障信息由对称分量法获知。将Zin2代入负序网络方程中,便可以计算出负序定子电压相量/>进而求得/>以及/>
以表1的IEEE-2和IEEE-5这两组感应电动机的参数数据作为算例,验证图3的实施方式的暂态响应。
表一:IEEE-2和IEEE-5的感应电动机的参数数据
类型 | R<sub>s</sub>(p.u.) | X<sub>s</sub>(p.u.) | R<sub>r</sub>(p.u.) | X<sub>r</sub>(p.u.) | X<sub>m</sub>(p.u.) | T<sub>j</sub>(s) | A |
IEEE-2 | 0.013 | 0.067 | 0.009 | 0.170 | 3.80 | 3.0 | 0.8 |
IEEE-5 | 0.077 | 0.107 | 0.079 | 0.098 | 2.22 | 1.48 | 0.46 |
(p.u.)表示相应数据为标幺值,负载转矩Tm=Aωm 2,ωm即为标幺转速,在进行有名值的计算时,该种实施方式中,采用的数据是:感应电动机的额定电压Un=110kV,额定电流In=152.5A。
设定电源在0.2s时发生持续时间为0.2s的对称和不对称电压跌落,计算感应电动机的暂态响应,验证本申请的有效性。电压跌落数据如下:对称电压跌落期间,Eeqa=Eeqb=Eeqc=0.p.u.,等价于电力系统发生了三相接地故障,不对称电压跌落期间,Eeqa=0.p.u.,Eeqb和Eeqc保持不变,等价于电力系统发生了A相接地故障。考虑到电磁暂态仿真是电力系统最精确的仿真手段,该种实施方式中便用其检验本申请方案的有效性和精度。
可参阅图4,为对称电压跌落下,采用IEEE-2的参数数据时,现有技术中的三阶机电暂态模型计算得到的感应电动机的暂态响应,与PSCAD/EMTDC电磁暂态仿真软件计算结果的比较。具体的,坐标系C41至坐标系C44依次示出了感应电动机的定子电压、滑差、有功消耗以及无功消耗。图4中的机电暂态即为三阶机电暂态模型计算得到的结果,电磁暂态即为电磁暂态仿真软件的计算结果。图5则为采用IEEE-5的参数数据的相应比较结果。
由图4及图5可以看出,受扰后由于定子电压下降,电动机电磁转矩下降,转速降低,滑差增大,功率消耗降低。扰动消失后,转速、功率等经历暂态过程后逐渐恢复正常。两种方法的计算结果吻合良好,这表明现有技术中的三阶机电暂态模型在对称故障下求得的感应电动机暂态响应较为准确。
可参阅图6以及图7,图6为不对称电压跌落下,采用IEEE-2的参数数据时,现有技术中的三阶机电暂态模型计算得到的感应电动机的暂态响应,与PSCAD/EMTDC电磁暂态仿真软件计算结果的比较。图6中的坐标系C61至坐标系C66依次示出了感应电动机的A、B、C相的定子电压、滑差、有功消耗以及无功消耗。图7则为采用IEEE-5的参数数据的相应比较结果。可以看出,图6中的无功功耗和定子电压的误差较大。图7中,滑差、有功消耗、无功消耗以及定子电压的计算精度都不足。表明在不对称扰动下,不计负序分量影响的传统的三阶机电暂态模型计算得到的感应电动机暂态响应并不准确。
可参阅图8以及图9,图8为不对称电压跌落下,采用IEEE-2的参数数据时,本申请的方案计算得到的感应电动机的暂态响应,与PSCAD/EMTDC电磁暂态仿真软件计算结果的比较。图8中的坐标系C81至坐标系C86依次示出了感应电动机的A、B、C相的定子电压、滑差、有功消耗以及无功消耗,图中的新方法即为本申请的方案。图9则为采用IEEE-5的参数数据的相应比较结果。可以看出,本申请的方案的计算结果与电磁暂态仿真结果吻合良好,验证了本申请方案的有效性。
相应于上面的方法实施例,本发明实施例还提供了一种感应电动机暂态响应的求解系统,下文描述的感应电动机暂态响应的求解系统与上文描述的感应电动机暂态响应的求解方法可相互对应参照。
参见图10所示,为本发明中一种感应电动机暂态响应的求解系统的结构示意图,包括:
参数数据获取模块201,用于获取感应电动机的参数数据;
负序阻抗计算模块202,用于根据参数数据计算出感应电动机的近似负序阻抗Zin2;
暂态响应计算模块205,用于联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出感应电动机的暂态响应;其中,感应电动机三阶机电暂态方程中的总电磁转矩Te为计入负序电磁转矩Te-的电磁转矩,有功消耗计算方程为计入负序有功消耗Pe-的方程,无功消耗计算方程为计入负序无功消耗Qe-的方程。
在本发明的一种具体实施方式中,负序相量计算模块203,具体用于:
通过电网系统的结构建立负序网络方程;
其中,Rs为参数数据中的感应电动机的定子电阻,X′为感应电动机的转子短路电抗,满足:X′=Xs+XrXm/(Xr+Xm);Xr为参数数据中的感应电动机的转子电抗,Xs为参数数据中的感应电动机的定子电抗,Xm为参数数据中的感应电动机的励磁电抗。
在本发明的一种具体实施方式中,负序阻抗计算模块202,具体用于:
通过忽略滑差变化构建感应电动机的负序稳态等值电路;
根据参数数据以及负序稳态等值电路,计算出感应电动机的负序等值阻抗,并将负序等值阻抗作为计算出的感应电动机的近似负序阻抗Zin2。
在本发明的一种具体实施方式中,计算出的近似负序阻抗Zin2表示为:
其中,Rs为参数数据中的感应电动机的定子电阻,Xr为参数数据中的感应电动机的转子电抗,Xm为参数数据中的感应电动机的励磁电抗,Xs为参数数据中的感应电动机的定子电抗,Rr为参数数据中的感应电动机的转子电阻。
在本发明的一种具体实施方式中,确定出的负序电磁转矩Te-表示为:
其中,real表示取相量的实部,*表示相量的共轭;
确定出的负序有功消耗Pe-表示为:
确定出的负序无功消耗Qe-表示为:
其中,imag表示取相量的虚部。
在本发明的一种具体实施方式中,计入负序电磁转矩Te-的总电磁转矩Te表示为:
计入负序有功消耗Pe-的有功消耗计算方程表示为:
计入负序无功消耗Qe-的无功消耗计算方程表示为:
相应于上面的方法和系统实施例,本发明实施例还提供了一种感应电动机暂态响应的求解设备以及一种计算机可读存储介质,参阅图11,该设备可以包括:
存储器301,用于存储暂态响应求解程序;
处理器302,用于执行暂态响应求解程序以实现上述任一实施例中的感应电动机暂态响应的求解方法的步骤。
计算机可读存储介质上存储有计算机程序,计算机程序被处理器执行时实现上述任一实施例中的感应电动机暂态响应的求解方法的步骤。这里所说的计算机可读存储介质包括随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的技术方案及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (8)
1.一种感应电动机暂态响应的求解方法,其特征在于,包括:
获取感应电动机的参数数据;
根据所述参数数据计算出所述感应电动机的近似负序阻抗Z in2;
联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出所述感应电动机的暂态响应;
其中,不对称故障期间所述感应电动机三阶机电暂态方程中的总电磁转矩T e为计入所述负序电磁转矩T e-的电磁转矩,所述有功消耗计算方程为计入负序有功消耗P e-的方程,所述无功消耗计算方程为计入负序无功消耗Q e-的方程;
通过所述电网系统的结构建立负序网络方程;
根据所述参数数据计算出所述感应电动机的近似负序阻抗Z in2,包括:
通过忽略滑差变化构建所述感应电动机的负序稳态等值电路;
根据所述参数数据以及所述负序稳态等值电路,计算出所述感应电动机的负序等值阻抗,并将所述负序等值阻抗作为计算出的所述感应电动机的近似负序阻抗Z in2;
联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出所述感应电动机的暂态响应,包括:
联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,基于机电暂态仿真方法,计算出所述感应电动机的暂态响应。
6.一种感应电动机暂态响应的求解系统,其特征在于,包括:
参数数据获取模块,用于获取感应电动机的参数数据;
负序阻抗计算模块,用于根据所述参数数据计算出所述感应电动机的近似负序阻抗Z in2;
负序相量计算模块,用于通过所述感应电动机所处的电网系统的电网数据,所述近似负序阻抗Z in2以及所述参数数据,计算出所述感应电动机的负序定子电压相量,负序定子电流相量/>以及负序转子暂态电势相量/>;
暂态响应计算模块,用于联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出所述感应电动机的暂态响应;其中,不对称故障期间所述感应电动机三阶机电暂态方程中的总电磁转矩T e为计入所述负序电磁转矩T e-的电磁转矩,所述有功消耗计算方程为计入负序有功消耗P e-的方程,所述无功消耗计算方程为计入负序无功消耗Q e-的方程;
通过所述电网系统的结构建立负序网络方程;
根据所述参数数据计算出所述感应电动机的近似负序阻抗Z in2,包括:
通过忽略滑差变化构建所述感应电动机的负序稳态等值电路;
根据所述参数数据以及所述负序稳态等值电路,计算出所述感应电动机的负序等值阻抗,并将所述负序等值阻抗作为计算出的所述感应电动机的近似负序阻抗Z in2;
联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,计算出所述感应电动机的暂态响应,包括:
联立正序电力网络方程,感应电动机定子电压方程,感应电动机三阶机电暂态方程,有功消耗计算方程以及无功消耗计算方程,基于机电暂态仿真方法,计算出所述感应电动机的暂态响应。
7.一种感应电动机暂态响应的求解设备,其特征在于,包括:
存储器,用于存储暂态响应求解程序;
处理器,用于执行所述暂态响应求解程序以实现如权利要求1至5任一项所述的感应电动机暂态响应的求解方法的步骤。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至5任一项所述的感应电动机暂态响应的求解方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811252348.4A CN109359266B (zh) | 2018-10-25 | 2018-10-25 | 感应电动机暂态响应的求解方法、系统、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811252348.4A CN109359266B (zh) | 2018-10-25 | 2018-10-25 | 感应电动机暂态响应的求解方法、系统、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109359266A CN109359266A (zh) | 2019-02-19 |
CN109359266B true CN109359266B (zh) | 2023-03-28 |
Family
ID=65346730
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811252348.4A Active CN109359266B (zh) | 2018-10-25 | 2018-10-25 | 感应电动机暂态响应的求解方法、系统、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109359266B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110501614A (zh) * | 2019-08-31 | 2019-11-26 | 广东电网有限责任公司 | 适用于电磁暂态仿真的电动机负荷暂态过程计算方法 |
CN112886570A (zh) * | 2021-01-18 | 2021-06-01 | 湖南大学 | 一种高速铁路牵引供电系统综合负荷等效建模方法 |
CN112928806B (zh) * | 2021-01-28 | 2022-09-13 | 湖南炬神电子有限公司 | 一种超小型电源实现方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB491934A (en) * | 1937-03-12 | 1938-09-12 | Charles Leslie Faudell | Improvements in or relating to electric oscillation generators |
CN101599644A (zh) * | 2009-07-09 | 2009-12-09 | 华北电力大学 | 以电网正常或非正常运行信息辨识等值综合负荷模型 |
CN106356877A (zh) * | 2016-09-20 | 2017-01-25 | 广东工业大学 | 一种针对三相负载不平衡的无源补偿装置及方法 |
CN107591839A (zh) * | 2017-09-25 | 2018-01-16 | 江苏大学 | 基于励磁状态分析的虚拟同步发电机低电压穿越控制方法 |
-
2018
- 2018-10-25 CN CN201811252348.4A patent/CN109359266B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB491934A (en) * | 1937-03-12 | 1938-09-12 | Charles Leslie Faudell | Improvements in or relating to electric oscillation generators |
CN101599644A (zh) * | 2009-07-09 | 2009-12-09 | 华北电力大学 | 以电网正常或非正常运行信息辨识等值综合负荷模型 |
CN106356877A (zh) * | 2016-09-20 | 2017-01-25 | 广东工业大学 | 一种针对三相负载不平衡的无源补偿装置及方法 |
CN107591839A (zh) * | 2017-09-25 | 2018-01-16 | 江苏大学 | 基于励磁状态分析的虚拟同步发电机低电压穿越控制方法 |
Non-Patent Citations (1)
Title |
---|
电网故障下双馈感应发电机暂态特性解析与短路计算模型;付元欢;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20180731;第1-82页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109359266A (zh) | 2019-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109359266B (zh) | 感应电动机暂态响应的求解方法、系统、设备及存储介质 | |
CN104485684B (zh) | 一种含双馈风电机组的电力系统故障电流计算方法 | |
CN102510263B (zh) | 基于抛载试验和数值差分的同步发电机实用参数辨识方法 | |
CN107783006B (zh) | 水轮发电机转子绕组匝间短路故障检测方法 | |
CN107069800A (zh) | 一种双馈风电场短路电流模型的建立方法 | |
CN103500269A (zh) | 一种双馈异步发电机组暂态短路电流的计算方法 | |
CN107451364B (zh) | 一种dfig风电场等值参数的辨识方法 | |
CN107121604B (zh) | 一种不对称短路电流直流分量衰减时间常数获取方法 | |
CN107453401A (zh) | 一种双馈风力发电机参数辨识方法 | |
WO2021036279A1 (zh) | 适用于电磁暂态仿真的电动机负荷暂态过程计算方法 | |
Kopse et al. | Applying a wide-area measurement system to validate the dynamic model of a part of European power system | |
CN111141995B (zh) | 一种基于比幅原理的线路双端稳态量测距方法和系统 | |
CN104820895B (zh) | 电压源-阻抗形式的双馈风力发电机等效模型 | |
CN102916647A (zh) | 一种在线测量异步电机转子时间常数的方法及装置 | |
CN112907075B (zh) | 一种电力系统综合负荷模型参数辨识方法 | |
CN110082596B (zh) | 基于高频暂态量的距离保护故障阻抗计算方法及系统 | |
CN106897514B (zh) | 一种全功率变换型新能源场站的短路电流计算模型的建立方法 | |
CN101788608B (zh) | 一种独立三绕组变压器电抗参数的估计方法 | |
CN111562464A (zh) | 中性点不接地系统的参数在线测试方法、存储介质及装置 | |
CN117081066A (zh) | 考虑大扰动的感应电机负荷暂态电压失稳分析方法及系统 | |
CN104267243B (zh) | 同步发电机内电势和电抗参数的测量方法及装置 | |
WO2022110290A1 (zh) | 一种基于变机械负载转矩的感应电动机转子滑差计算方法 | |
CN106936359A (zh) | 一种机电暂态同步机模型及初始q轴饱和参数的确定方法 | |
Cari et al. | A constrained minimization approach for the estimation of parameters of transient generator models | |
CN106842021B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |