CN116402017B - 一种用于电力系统频率动态分析的双馈风机简化建模方法 - Google Patents
一种用于电力系统频率动态分析的双馈风机简化建模方法 Download PDFInfo
- Publication number
- CN116402017B CN116402017B CN202310643779.8A CN202310643779A CN116402017B CN 116402017 B CN116402017 B CN 116402017B CN 202310643779 A CN202310643779 A CN 202310643779A CN 116402017 B CN116402017 B CN 116402017B
- Authority
- CN
- China
- Prior art keywords
- doubly
- equivalent
- fed
- fan
- fed fan
- 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 100
- 238000004458 analytical method Methods 0.000 title claims abstract description 19
- 230000033001 locomotion Effects 0.000 claims abstract description 23
- 230000001052 transient effect Effects 0.000 claims abstract description 10
- 230000008569 process Effects 0.000 claims description 40
- 230000001360 synchronised effect Effects 0.000 claims description 40
- 230000010354 integration Effects 0.000 claims description 12
- 230000005284 excitation Effects 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 4
- 230000008093 supporting effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
-
- 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
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Power Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种用于电力系统频率动态分析的双馈风机简化建模方法,属于电力系统频率控制技术领域;包括以下步骤:获取双馈风机的一次、二次系统参数和稳态运行点;定义双馈风机的等效内电势;列写双馈风机等效内电势与机端电压、电机定子侧向电网注入的电流之间的关系;获得双馈风机等效转速的表达式;组合双馈风机内部各模块,得到线性化的双馈风机等效转子运动模型;获得双馈风机参与系统惯性中心频率动态的模型;获得双馈风机参与扰动功率分配的模型。本发明解析了双馈风机的等效暂态电抗和等效惯量,反映双馈风机与电网之间的电气距离和对电网的惯量支撑作用,对电网调度运行人员有重要参考价值。
Description
技术领域
本发明涉及电力系统频率控制技术领域,尤其是涉及一种用于电力系统频率动态分析的双馈风机简化建模方法。
背景技术
随着“碳达峰碳中和”目标的持续推进,大规模新能源接入给电力系统频率稳定带来了严峻挑战。一方面,新能源通过电力电子设备接入电网,无法主动为系统提供惯量支撑,系统惯量水平降低,扰动后频率变化率(RoCoF)显著增加,频率失稳风险增加;另一方面,新能源规模化集中接入导致系统惯量空间分布不均匀,系统频率具有空间分布特性,局部地区频率动态特性恶化严重,局部故障可能引起全网频率连锁崩溃。
新能源对系统频率动态的影响仍是一个未知的问题,无法定量解释新能源对系统惯性中心(COI)频率动态和频率空间分布的作用机制。特别是广泛应用的双馈风力发电机(DFIG),其动态特性由双馈式感应电机和多种控制回路共同决定,并且电机定子与电网直接相连,在所有新能源中具有最复杂的动态特性。DFIG的详细模型阶数高且微分代数方程耦合,只能用于仿真分析电力系统频率动态的过程,研究DFIG对系统频率动态作用机制的关键点在于建立合适的DFIG简化模型。
现有DFIG简化建模的技术路线主要包括两类。其一,基于系统COI频率动态的DFIG简化建模。比如,考虑DFIG机械转子运动、MPPT控制和其他动态模块,采用虚拟惯量(VI)控制的DFIG可以建模为简化传递函数或者进一步提取为等效惯量常数。DFIG简化传递函数可以描述机端频率与其有功功率输出之间的关系,并在COI坐标系下频率响应模型(SFR)中计及DFIG动态特性。然而,基于系统COI的技术路线无法考虑系统频率空间分布特性。其二,支撑系统频率空间分布研究的DFIG简化建模。比如,将DFIG频率时间尺度特性建模为戴维南等值电路方程,能够描述DFIG与系统之间的电气距离。然而,现有技术中DFIG内电势与转子电流有关,无法解析扰动瞬间DFIG分配的扰动功率以及DFIG内电势对电网频率的贡献。
因此,需要一种用于电力系统频率动态分析的DFIG简化建模方法,能够同时考虑DFIG对COI频率的贡献和对频率空间分布特性的影响,支撑新能源电力系统频率动态精确分析与控制。
发明内容
本发明的目的是提供一种用于电力系统频率动态分析的双馈风机简化建模方法,解决上述背景技术中提到的技术问题。
为实现上述目的,本发明提供了一种用于电力系统频率动态分析的双馈风机简化建模方法,包括以下步骤:
S1、获取双馈风机的参数,包括双馈风机一次系统参数、双馈风机二次系统参数和双馈风机的稳态运行点;
S2、定义双馈风机的等效内电势;
S3、列写双馈风机等效内电势与机端电压、电机定子侧向电网注入的电流之间的关系,将该关系建模为相量形式的方程,构成双馈风机的等效电路模型;
S4、获得双馈风机等效转速的表达式;
S5、根据线性化的双馈风机等效转速表达式,组合双馈风机内部各模块,将各模块的输入统一为有功功率,得到非线性形式的双馈风机等效转子运动模型,对其线性化处理得到线性化的双馈风机等效转子运动模型;
S6、获得双馈风机参与系统惯性中心频率动态的模型;
S7、获得双馈风机参与扰动功率分配的模型。
优选的,所述S2中,双馈风机等效内电势表示为系统频率动态时间尺度下的状态变量的代数组合,如下所示:
;
式中,表示双馈风机等效内电势,/>表示等效内电势幅值,/>表示等效内电势相角,/>表示虚数单位,/>表示代数组合,/>表示锁相环的输出相角,/>表示风机的机械转子转速,/>表示有功控制回路积分环节的状态变量。
优选的,所述S3中,双馈风机相量形式的等效电路模型如下式所示:
;
;
;
;
;
式中,是双馈风机机端电压,/>是双馈风机定子向电网注入的有功,/>是双馈风机的等效电阻,/>是等效暂态电抗,/>是双馈电机的额定转速,/>是定子电抗,/>是励磁电抗,/>是有功外环控制的比例控制参数,/>是有功外环控制的积分控制参数,/>是机端电压控制指令,/>是最大功率跟踪控制系数,/>是扰动前双馈风机电机转子q轴电流;
和/>分别是双馈风机向电网注入的总有功和总无功,满足/>和/>,式中/>、/>分别表示双馈电机定子的有功和无功。
优选的,所述S4中,线性化的双馈风机等效转速表达式:
;
;
;
式中,和/>是线性化系数,/>是电网额定频率,/>和/>分别是风机的机械转速和有功控制积分环节的状态变量在扰动前的值,/>是锁相环输出频率。
优选的,所述S5中,线性化的双馈风机等效转子运动模型为:
;
;
;
;
;
;
;
;
;
;
;
;
;
式中,等效转速分量、/>和/>分别对应线性化的双馈风机等效转速表达式中的/>、/>和/>,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>和/>对应的等效机械功率,/>是线性化后的双馈风机定子有功功率;
等效转子运动模型中等效惯量常数为:
;
;
;
式中,是锁相环的比例控制参数,/>是锁相环的积分控制参数,/>是风力机的机械转动惯量;
常数项为:
;
;
;
高阶分量为:
;
。
优选的,对于最大功率跟踪控制的线性化,最大功率跟踪的线性化系数表达式为:
;
对于双馈风机内部有功潮流方程的线性化,有功潮流方程的线性化系数表达式为:
;
式中,和/>分别是双馈风机等效内电势和等效暂态电抗在扰动前的数值;
对于风机空气动力学模型的线性化,考虑到风机空气动力学模型的多样性,若其模型为:
;
;
;
;
式中,是风机的输出转矩,/>是叶尖速比,/>是一个中间变量,/>是风速,/>是额定风速,/>是风机特性/>的额定值,/>、/>、/>、/>、/>是空气动力学参数;
风机空气动力学模型的线性化系数表达式为:
;
式中,是稳态风速,/>是额定叶尖速比。
优选的,所述S6中,双馈风机参与系统惯性中心频率动态的模型类比同步机系统获得,惯性中心频率动态和系统总惯量/>通过下式计算:
;
;
式中,g和w分别表示同步机和双馈风机的序号,是同步机的转子转速,/>和分别是同步机和双馈风机的额定容量,/>是同步机的惯量常数,/>是双馈风机的总等效惯量;/>是系统中第w台双馈风机的等效转速;
的表达式如下:
;
式中,和/>分别是由线性化系数/>和/>归一化的等效惯量分量。
优选的,所述S7中,双馈风机参与扰动功率分配的模型类比同步机系统获得,系统中第w个双馈风机拾取的扰动功率用双馈风机和扰动位置之间的同步功率系数来描述,如下所示:
;
;
式中,是同步机的同步功率系数,/>是扰动前扰动点电压的相角,/>是双馈风机等效内电势与扰动点之间的等效电抗,/>是系统扰动功率,/>是扰动前双馈风机等效内电势的幅值,/>是扰动前扰动点电压的幅值,/>是扰动前双馈风机等效内电势的相角。
本发明所述的一种用于电力系统频率动态分析的双馈风机简化建模方法的有益效果在于:
本发明所提方法类比同步机二阶经典模型建立了双馈风机的简化模型,该模型的优势在于模型能够直接用于电力系统频率动态分析,包括双馈风机参与系统惯性中心频率动态、双馈风机参与扰动功率分配等方面,模型中解析了双馈风机的等效暂态电抗和等效惯量,能够反映双馈风机与电网之间的电气距离和对电网的惯量支撑作用,对电网调度运行人员有重要参考价值,具有实际应用潜力。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的流程图;
图2为本发明相量形式的双馈风机等效电路模型;
图3为本发明根据等效转速组合各模块后的双馈风机动力学模型;
图4为本发明非线性形式的双馈风机等效转子运动模型;
图5为本发明线性化后的双馈风机等效转子运动模型;
图6为本发明双馈风机的等效转子运动模型;
图7为本发明双馈风场-同步机测试系统拓扑示意图;
图8为本发明模型和双馈风机详细模型的输出有功对比图;
图9为本发明模型和双馈风机详细模型的并网点频率对比图;
图10为双馈风机和同步机的输出有功图;
图11为风机和同步机的(等效)转速以及系统惯性中心频率对比图。
具体实施方式
以下通过附图和实施例对本发明的技术方案作进一步说明。
图1为本发明的流程图。如图1所示,一种用于电力系统频率动态分析的双馈风机简化建模方法,包括以下步骤:
S1、获取双馈风机的参数,包括双馈风机一次系统参数、双馈风机二次系统参数和双馈风机的稳态运行点。
双馈风机一次系统参数:双馈电机的额定转速、励磁电抗/>、定子电抗/>和定子电阻/>,风力机的机械转动惯量/>、空气动力学参数(/>、/>、/>、/>、/>)、额定风速/>、风力机特性/>的额定值/>。
双馈风机二次系统参数:有功外环控制的比例控制参数和积分控制参数/>,锁相环的比例控制参数/>和积分控制参数/>,最大功率跟踪控制系数/>。
双馈风机的稳态运行点:风速,无功功率运行点/>,机端电压控制指令/>。
S2、定义双馈风机的等效内电势。
双馈风机等效内电势表示为系统频率动态时间尺度下的状态变量的代数组合,状态变量包括锁相环(PLL)的输出相角、风机的机械转子转速/>和有功控制回路积分环节的状态变量/>。本发明类比同步机的二阶经典模型建立双馈风机的简化模型,双馈风机的等效内电势可以表示为上述三个状态变量的组合:
如下所示:
;
式中,表示双馈风机等效内电势,/>表示等效内电势幅值,/>表示等效内电势相角,/>表示代数组合,/>表示锁相环的输出相角,/>表示风机的机械转子转速,表示有功控制回路积分环节的状态变量。
S3、列写双馈风机等效内电势与机端电压、电机定子侧向电网注入的电流之间的关系,将该关系建模为相量形式的方程,构成双馈风机的等效电路模型。模型中等效电抗参数用于表示双馈风机等效内电势与电网之间的电气距离。
在该步骤中,以锁相环输出相角为坐标参考系时,双馈电机dq轴形式的电路方程为:
;
;
式中,、/>和/>分别是锁相环相角参考系下的机端电压、定子电流和转子电流,下标d、q分别表示d轴和q轴;/>为机械转子额定转速;/>和/>分别是励磁电抗和定子电抗;/>是定子电阻。
在上式基础上,考虑最大功率跟踪(MPPT)控制和有功外环控制环节,等效内电势与机端电压/>、电机定子侧向电网注入的电流/>之间的关系如下式给出,等式左侧与等效内电势/>有关,该式即为dq轴分量形式的双馈风机等效电路模型。
;
;
式中,和/>分别是有功外环控制的比例控制参数和积分控制参数;/>是最大功率跟踪控制系数;/>是扰动前双馈风机电机转子q轴电流。
为了解决存在导致方程dq轴不对称的问题,构造与双馈风机电机电流有关的恒等式:
;
式中,表示共轭转置运算,/>表示共轭运算;/>、/>和/>分别是双馈电机定子的复功率、有功和无功;/>和/>分别为锁相环相角坐标系下双馈电机的机端电压相量和电流相量,其中,锁相环相角坐标系下机端电压相角/>为0,即/>。
上式的虚部和实部分别提取如下:
;
;
上式中,常数0和1被分解为双馈电机dq轴的电流,常数1乘以项后可引入d轴方程,常数0与/>的乘积可引入q轴方程。dq轴分量形式的双馈风机等效电路模型可以重新表述如下:
;
;
此时,上式dq轴对称,可以写成相量形式,将相量形式方程旋转锁相环相角后转换到电网电压相位坐标系,双馈电机相量形式的等效电路模型下式所示。如图2所示。
;
式中,双馈风机的等效电阻、等效暂态电抗/>、等效内电势幅值/>和等效内电势相角/>如下:
;
;
;
;
式中,是双馈风机机端电压,/>是双馈风机定子向电网注入的有功,/>是机端电压控制指令。
和/>分别是双馈风机向电网注入的总有功和总无功,满足/>和/>。
S4、获得双馈风机等效转速的表达式。
双馈风机等效转速为等效内电势相角/>的导数,需要对三个状态变量分别求偏导数,即锁相环的输出相角/>、风机的机械转子转速/>和有功控制回路积分环节的状态变量/>。
等效转速表达式如下所示:
;
式中,锁相环相角替换为锁相环输出频率/>。
将上式线性化处理得到线性化的双馈风机等效转速表达式:
;
;
;
式中,和/>是线性化系数,/>是电网额定频率,/>和/>分别是风机的机械转速和有功控制积分环节的状态变量在扰动前的值。
S5、根据线性化的双馈风机等效转速表达式,组合双馈风机内部各模块,将各模块的输入统一为有功功率,得到非线性形式的双馈风机等效转子运动模型,对其线性化处理得到线性化的双馈风机等效转子运动模型。
根据双馈风机等效转速表达式,组合双馈风机内部各模块,得到如图3所示的框图。图3中,s是拉普拉斯算子,和/>分别是锁相环的比例控制参数和积分控制参数,是风机机械转动惯量,/>是机端电压控制指令,v是风速,/>是机端电压的相角,/>和/>分别是风机的输出转矩和电磁转矩,/>是锁相环相角坐标系下双馈风机机端电压的q轴分量,/>是锁相环积分控制的状态变量,/>和/>分别是是双馈风机输出有功功率实际值和控制指令。图3中/>是风机的空气动力学模型,具体表达式如下所示:
;
;
;
;
式中,、/>、/>、/>和/>均为风机的空气动力学参数,/>是额定风速,/>是描述风力机特性的变量,/>是其额定值,/>是叶尖速比,/>是一个中间变量。
类比同步机的二阶经典模型,图3中各模块的输入需要转化有功功率,首先,对于机械转子运动模块,电磁转矩/>可以表示为:
;
对于有功外环控制的积分环节,双馈风机总有功可以表示为:
;
对于锁相环,锁相环输出相角可以通过双馈风机内部有功潮流方程计算得到,如下:
;
式中,和/>分别是双馈风机等效内电势和等效暂态电抗在扰动前的数值,如下所示:/>
;
;
式中,和/>分别为双馈风机扰动前总有功功率和无功功率。
此时,图3所示框图可以变为图4,图4即为非线性形式的双馈风机等效转子运动模型。进一步,对图4中各非线性模块线性化,包括风机的空气动力学模型:
;
式中,是空气动力学模型的线性化系数;/>是稳态风速,/>是额定叶尖速比。
最大功率跟踪控制:
;
;
式中,是最大功率跟踪的线性化系数。
有功外环控制的积分环节:
;
锁相环:
;
双馈风机的有功潮流方程:
;
;
式中,是有功潮流方程的线性化系数。
结合上述线性化方程,对图4拉普拉斯变换并做简单代数运算处理,得到线性化后的双馈风机等效转子运动模型,由图5表示。图5引入的等效转速分量、/>和/>对应线性化等效转速表达式中的/>、/>和/>。
类比同步机的二阶经典模型,图5可以经代数运算分解为图6。图6对应的数学方程如下:
;
;
;
;
;
;
;
;
;
;
;
;
;
式中,等效转速分量、/>和/>分别对应线性化的双馈风机等效转速表达式中的/>、/>和/>,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>和/>对应的等效机械功率,/>是线性化后的双馈风机定子有功功率;
等效转子运动模型中等效惯量常数为:
;
;
;
式中,是锁相环的比例控制参数,/>是锁相环的积分控制参数,/>是风力机的机械转动惯量;/>
常数项为:
;
;
;
高阶分量为:
;
。
S6、获得双馈风机参与系统惯性中心频率动态的模型。
双馈风机简化模型与同步发电机的二阶经典模型形式相同,参与系统惯性中心频率动态的模型与同步机相同。惯性中心频率动态和系统总惯量/>通过下式计算:
;
;
式中,g和w分别表示同步机和双馈风机的序号,是同步机的转子转速,/>和分别是同步机和双馈风机的额定容量,/>是同步机的惯量常数,/>是双馈风机的总等效惯量;/>是系统中第w台双馈风机的等效转速;
的表达式如下:
;
式中,和/>分别是由线性化系数/>和/>归一化的等效惯量分量。
S7、获得双馈风机参与扰动功率分配的模型。
由于双馈风机的等效内电势建模为状态变量的组合,双馈风机在系统频率动态时间尺度上的动态类型类似于电压源。系统扰动功率在各台发电机之间的分配是频率空间尺度变化的前提,从本发明所提模型看,由于扰动之后双馈风机等效内电势相角不会突变,双馈风机可以主动拾取扰动不平衡功率。与同步机类似,系统中第w个双馈风机拾取的扰动功率可以用双馈风机和扰动位置之间的同步功率系数(SPC)/>来描述,如下所示:/>
;
;
式中,是同步机的同步功率系数,/>是扰动前扰动点电压的相角,/>是双馈风机等效内电势与扰动点之间的等效电抗,/>是系统扰动功率,/>是扰动前双馈风机等效内电势的幅值,/>是扰动前扰动点电压的幅值,/>是扰动前双馈风机等效内电势的相角。
实施例
如图7所示,300 MVA同步发电机和60台1.5 MW双馈风机组成的风电场通过长距离输电线路相连。该系统验证一种用于电力系统频率动态分析的双馈风机简化建模方法的步骤如下:
1、获取双馈风机的参数:双馈风机一次系统参数:双馈电机的额定转速p.u.、励磁电抗/>p.u.、定子电抗/>p.u.和定子电阻/>p.u.,风力机的机械转动惯量/>s、空气动力学参数(/>、/>、/>、/>、)、额定风速/>m/s、风力机特性/>的额定值/>。双馈风机二次系统参数:有功外环控制的比例控制参数/>和积分控制参数/>,锁相环的比例控制参数/>和积分控制参数/>,最大功率跟踪控制系数/>。双馈风机的稳态运行点:风速/>m/s,无功功率运行点/>,机端电压控制指令/>p.u.。
2、对于双馈风机的等效电路模型,计算等效电阻、等效暂态电抗的数值:
;
;
3、对于双馈风机的等效转子运动模型,首先,计算转子转速线性化系数的数值:
;
;
接着,计算双馈风机等效惯量常数的数值:
;
;
;
常数项数值为:
;
;
;
高阶分量数值为:
;
;
4、在双馈风机参与系统惯性中心频率动态分析中,双馈风机的总等效惯量为:
;
5、在双馈风机参与扰动功率分配分析中,双馈风机分配扰动功率比例为:
;
提供双馈风机的详细模型作为比较,即,将双馈风机的详细模型替换为本发明所提出的简化模型,并在相同的扰动条件下重复仿真。当风场互联点(PCC)的负荷增加10 MW时,风场有功功率和系统频率响应分别在图8和图9中给出。可以看出,本发明所提简化模型的有功功率输出与详细模型非常接近,与详细模型的频率曲线基本相同。采用下式定义的误差指标来量化动态过程的相似性,在惯量响应时间尺度(0+时刻至1秒),有功功率输出和频率的误差指标分别为1.95%和0.10%,在一次调频时间尺度(1秒-18秒),误差指标分别为5.63%和0.85%,进一步验证了本发明所提模型的准确性。
;
式中,和/>分别是误差评价的结束时间和开始时间;/>和/>分别是详细模型和发明所提模型的曲线;/>为这两个模型的曲线差的最大值。
发明所提模型用于电力系统频率动态分析时,扰动后有功功率和频率动态分别如图10和图11所示。在双馈风机参与系统惯性中心频率动态分析中,根据所提模型计算的惯性中心频率响应曲线如图10所示,惯性中心坐标系下系统总惯量为1.14秒。需要注意,尽管同步机1秒的惯量常数和双馈风机1.59秒的等效惯量常数数值近似,但同步机的额定容量远大于双馈风机,因此惯性中心频率和系统总惯量主要受到同步机的影响,可以看到,惯性中心频率与同步机转子转速更接近,系统总惯量与同步机惯量常数更接近。在双馈风机参与扰动功率分配分析中,如图11可以看到,双馈风机实际分配的扰动功率与根据发明所提模型理论计算值(6.13%)基本吻合。
因此,本发明采用上述一种用于电力系统频率动态分析的双馈风机简化建模方法,解析了双馈风机的等效暂态电抗和等效惯量,能够反映双馈风机与电网之间的电气距离和对电网的惯量支撑作用,对电网调度运行人员有重要参考价值。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。
Claims (1)
1.一种用于电力系统频率动态分析的双馈风机简化建模方法,其特征在于,包括以下步骤:
S1、获取双馈风机的参数,包括双馈风机一次系统参数、双馈风机二次系统参数和双馈风机的稳态运行点;
S2、定义双馈风机的等效内电势;
S3、列写双馈风机等效内电势与机端电压、电机定子侧向电网注入的电流之间的关系,将该关系建模为相量形式的方程,构成双馈风机的等效电路模型;
S4、获得双馈风机等效转速的表达式;
S5、根据线性化的双馈风机等效转速表达式,组合双馈风机内部各模块,将各模块的输入统一为有功功率,得到非线性形式的双馈风机等效转子运动模型,对其线性化处理得到线性化的双馈风机等效转子运动模型;
S6、获得双馈风机参与系统惯性中心频率动态的模型;
S7、获得双馈风机参与扰动功率分配的模型;
所述S2中,双馈风机等效内电势表示为系统频率动态时间尺度下的状态变量的代数组合,如下所示:
;
式中,表示双馈风机等效内电势,/>表示等效内电势幅值,/>表示等效内电势相角,/>表示虚数单位,/>表示代数组合,/>表示锁相环的输出相角,/>表示风机的机械转子转速,/>表示有功控制回路积分环节的状态变量;
所述S3中,双馈风机相量形式的等效电路模型如下式所示:
;
;
;
;
;
式中,是双馈风机机端电压,/>是双馈风机定子向电网注入的有功,/>是双馈风机的等效电阻,/>是等效暂态电抗,/>是双馈电机的额定转速,/>是定子电抗,/>是励磁电抗,/>是有功外环控制的比例控制参数,/>是有功外环控制的积分控制参数,/>是机端电压控制指令,/>是最大功率跟踪控制系数,/>是扰动前双馈风机电机转子q轴电流,是定子电阻;
和/>分别是双馈风机向电网注入的总有功和总无功,满足/>和/>,式中/>、/>分别表示双馈电机定子的有功和无功;
所述S4中,线性化的双馈风机等效转速表达式:
;
;
;
式中,和/>是线性化系数,/>是电网额定频率,/>和/>分别是风机的机械转速和有功控制积分环节的状态变量在扰动前的值,/>是锁相环输出频率,t是时间;
所述S5中,线性化的双馈风机等效转子运动模型为:
;
;
;
;
;
;
;
;
;
;
;
;
;
式中,等效转速分量、/>和/>分别对应线性化的双馈风机等效转速表达式中的/>、/>和/>,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>中的连续分量和非连续分量,/>和/>分别是/>和/>对应的等效机械功率,/>是线性化后的双馈风机定子有功功率,/>是锁相环输出频率等效惯量常数,/>是风机的机械转子转速等效惯量常数,/>是有功控制回路积分环节的状态变量等效惯量常数;
等效转子运动模型中等效惯量常数为:
;
;
;
式中,是锁相环的比例控制参数,/>是锁相环的积分控制参数,/>是风力机的机械转动惯量;
常数项为:
;
;
;
式中,是锁相环输出频率常数,/>是风机的机械转子转速常数,/>是有功控制回路积分环节的状态变量常数;
高阶分量为:
;
;
是风机的机械转子转速常数的高阶分量传递函数,/>是有功控制回路积分环节的高阶分量传递函数,s是拉普拉斯算子;
对于最大功率跟踪控制的线性化,最大功率跟踪的线性化系数表达式为:
;
对于双馈风机内部有功潮流方程的线性化,有功潮流方程的线性化系数表达式为:
;
式中,和/>分别是双馈风机等效内电势和等效暂态电抗在扰动前的数值;
对于风机空气动力学模型的线性化,考虑到风机空气动力学模型的多样性,若其模型为:
;
;
;
;
式中,是风机的输出转矩,/>是叶尖速比,/>是一个中间变量,/>是风速,/>是额定风速,/>是风机特性/>的额定值,/>、/>、/>、/>、/>是空气动力学参数;
风机空气动力学模型的线性化系数表达式为:
;
式中,是稳态风速,/>是额定叶尖速比;
所述S6中,双馈风机参与系统惯性中心频率动态的模型类比同步机系统获得,惯性中心频率动态和系统总惯量/>通过下式计算:
;
;
式中,和/>分别表示同步机和双馈风机的序号,/>是同步机的转子转速,/>和分别是同步机和双馈风机的额定容量,/>是同步机的惯量常数,/>是双馈风机的总等效惯量;/>是系统中第w台双馈风机的等效转速;
的表达式如下:
;
式中,和/>分别是由线性化系数/>和/>归一化的等效惯量分量;
所述S7中,双馈风机参与扰动功率分配的模型类比同步机系统获得,系统中第w个双馈风机拾取的扰动功率用双馈风机和扰动位置之间的同步功率系数/>来描述,如下所示:
;
;
式中,是同步机的同步功率系数,/>是扰动前扰动点电压的相角,/>是双馈风机等效内电势与扰动点之间的等效电抗,/>是系统扰动功率,/>是扰动前双馈风机等效内电势的幅值,/>是扰动前扰动点电压的幅值,/>是扰动前双馈风机等效内电势的相角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310643779.8A CN116402017B (zh) | 2023-06-02 | 2023-06-02 | 一种用于电力系统频率动态分析的双馈风机简化建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310643779.8A CN116402017B (zh) | 2023-06-02 | 2023-06-02 | 一种用于电力系统频率动态分析的双馈风机简化建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116402017A CN116402017A (zh) | 2023-07-07 |
CN116402017B true CN116402017B (zh) | 2023-08-22 |
Family
ID=87008976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310643779.8A Active CN116402017B (zh) | 2023-06-02 | 2023-06-02 | 一种用于电力系统频率动态分析的双馈风机简化建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116402017B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117350089B (zh) * | 2023-12-06 | 2024-04-02 | 山东大学 | 一种双馈风电场全动态阻抗模型的构建方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108418242A (zh) * | 2018-01-23 | 2018-08-17 | 华南理工大学 | 一种基于相似度同调的双馈风电机动态等值方法 |
CN112928781A (zh) * | 2021-03-23 | 2021-06-08 | 南方电网科学研究院有限责任公司 | 双馈风机暂态稳定控制方法、系统、计算机设备和介质 |
CN115622149A (zh) * | 2022-09-06 | 2023-01-17 | 华中科技大学 | 双馈风机参与一次调频的系统频率响应建模方法和系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105591403B (zh) * | 2016-02-01 | 2018-07-10 | 华南理工大学 | 基于逻辑开关控制的双馈风机多回路切换控制系统 |
CN112671006A (zh) * | 2020-12-18 | 2021-04-16 | 国网江苏省电力有限公司经济技术研究院 | 一种评估海上风电场柔性直流输电系统谐振稳定性的方法 |
-
2023
- 2023-06-02 CN CN202310643779.8A patent/CN116402017B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108418242A (zh) * | 2018-01-23 | 2018-08-17 | 华南理工大学 | 一种基于相似度同调的双馈风电机动态等值方法 |
CN112928781A (zh) * | 2021-03-23 | 2021-06-08 | 南方电网科学研究院有限责任公司 | 双馈风机暂态稳定控制方法、系统、计算机设备和介质 |
CN115622149A (zh) * | 2022-09-06 | 2023-01-17 | 华中科技大学 | 双馈风机参与一次调频的系统频率响应建模方法和系统 |
Non-Patent Citations (1)
Title |
---|
基于相似度同调的双馈风电场聚合研究;朱林;王贝;陈达;张健;;电力建设(第07期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116402017A (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tang et al. | Modeling of DFIG-based wind turbine for power system transient response analysis in rotor speed control timescale | |
BOUDJEMA et al. | Second order sliding mode control of a dual-rotor wind turbine system by employing a matrix converter | |
El Mourabit et al. | Nonlinear backstepping control for PMSG wind turbine used on the real wind profile of the Dakhla‐Morocco city | |
Bhattarai et al. | Parametrically robust dynamic speed estimation based control for doubly fed induction generator | |
CN101119095A (zh) | 基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统 | |
CN116402017B (zh) | 一种用于电力系统频率动态分析的双馈风机简化建模方法 | |
CN106786673B (zh) | 双馈风机串补输电系统次同步谐振的抑制方法及装置 | |
CN114006383B (zh) | 含vsg的双馈风机并网系统建模与次同步振荡分析的方法 | |
CN106786672B (zh) | 抑制串补输电系统次同步谐振的方法及系统 | |
Abdeen et al. | A recent analytical approach for analysis of sub-synchronous resonance in doubly-fed induction generator-based wind farm | |
CN110970925A (zh) | 基于双馈风机快速有功功率调制改善系统阻尼及建模方法 | |
Laghridat et al. | A novel adaptive active disturbance rejection control strategy to improve the stability and robustness for a wind turbine using a doubly fed induction generator | |
Djoudi et al. | Adaptive sensorless SM-DPC of DFIG-based WECS under disturbed grid: Study and experimental results | |
CN113131522A (zh) | 双馈风力发电机的虚拟惯量控制及稳定性分析方法 | |
Douadi et al. | Improvement performances of active and reactive power control applied to DFIG for variable speed wind turbine using sliding mode control and FOC | |
CN116054190A (zh) | 一种双馈风电场控制参数对次同步振荡影响的分析方法 | |
CN110417047B (zh) | 基于复转矩系数分析双馈风机ssci阻尼特性的方法 | |
CN112994111A (zh) | 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法 | |
Shanoob et al. | Wind turbine transient response and fault ride‐through improvements with optimal control | |
Yan et al. | Transient modelling of doubly‐fed induction generator based wind turbine on full operation condition and rapid starting period based on low voltage ride‐through testing | |
CN113708367B (zh) | 一种基于一致性算法的电力系统分布式协同控制方法 | |
CN112886611B (zh) | 一种直驱风机并网系统的次同步振荡抑制方法 | |
CN112994113B (zh) | 一种双馈风力发电机序阻抗建模方法及系统 | |
Oualah et al. | Super-twisting sliding mode control for brushless doubly fed reluctance generator based on wind energy conversion system | |
Xing et al. | On-line Identification of Equivalent Inertia for DFIG Wind Turbines Based on Extended Kalman Filters |
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 |