CN114154293B - 基于实时数据的直驱风电场等值建模方法和装置 - Google Patents

基于实时数据的直驱风电场等值建模方法和装置 Download PDF

Info

Publication number
CN114154293B
CN114154293B CN202111220530.3A CN202111220530A CN114154293B CN 114154293 B CN114154293 B CN 114154293B CN 202111220530 A CN202111220530 A CN 202111220530A CN 114154293 B CN114154293 B CN 114154293B
Authority
CN
China
Prior art keywords
matrix
wind turbine
impedance characteristic
turbine generator
time sequence
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
Application number
CN202111220530.3A
Other languages
English (en)
Other versions
CN114154293A (zh
Inventor
王清玲
黄伟煌
李力
彭发喜
陈辉祥
陈怡静
李岩
许树楷
曹彦朝
赵晓斌
周彦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China South Power Grid International Co ltd
Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
China South Power Grid International Co ltd
Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China South Power Grid International Co ltd, Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN202111220530.3A priority Critical patent/CN114154293B/zh
Publication of CN114154293A publication Critical patent/CN114154293A/zh
Application granted granted Critical
Publication of CN114154293B publication Critical patent/CN114154293B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Control Of Eletrric Generators (AREA)
  • Wind Motors (AREA)

Abstract

本申请提出了一种基于实时数据的直驱风电场等值建模方法，涉及电力系统风电场建模技术领域，其中，该方法包括：获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性曲线的参数；使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。采用上述方案的本申请可以高效、快速地获取直驱馈风电场的等值阻抗，为风电场接入电网的影响分析提供重要依据。

Description

基于实时数据的直驱风电场等值建模方法和装置
技术领域
本申请涉及电力系统风电场建模技术领域，尤其涉及一种基于实时数据的直驱风电场等值建模方法和装置。
背景技术
随着千万千瓦级风电基地的陆续建成，大规模风电机组集中并网给电力系统安全稳定运行带来巨大挑战。构建能够准确描述大规模风电场整体特性的等值模型，是研究高比例风电系统运行和控制的基础，风电场详细模型的等值是风电场动态等值的重要内容。风电场的动态等值建模已成为分析大型风电场并网特征的重要研究手段，目前对单个永磁直驱发电机(PMSG)的特性研究已经非常广泛，但是对风电场的整体特性还鲜有介绍。
发明内容
本申请旨在至少在一定程度上解决相关技术中的技术问题之一。
为此，本申请的第一个目的在于提出一种基于实时数据的直驱风电场等值建模方法，解决了现有方法对风电场的整体特性研究较少的问题，高效、精确地获取了风电系统在宽频域范围的阻抗特性，为大规模风电机组接入电网的影响分析提供了重要依据。
本申请的第二个目的在于提出一种基于实时数据的直驱风电场等值建模装置。
本申请的第三个目的在于提出一种非临时性计算机可读存储介质。
为达上述目的，本申请第一方面实施例提出了一种基于实时数据的直驱风电场等值建模方法，包括：获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性曲线的参数；使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。
可选地，在本申请的一个实施例中，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线，包括以下步骤：
根据电流瞬时值时间序列和风电组机端电压瞬时值时间序列分别生成第一矩阵和第二矩阵，其中，第一矩阵为电流瞬时值时间序列矩阵，第二矩阵为风电组机端电压瞬时值时间序列矩阵；
将第一矩阵分解生成第一目标矩阵和第二目标矩阵，将第二矩阵分解生成第三目标矩阵和第四目标矩阵；
对第一目标矩阵和第二目标矩阵、第三目标矩阵和第四目标矩阵进行计算生成单台风电机组阻抗特性曲线。
可选地，在本申请的一个实施例中，第一矩阵表示为：
其中，N为时间序列长度，当N＝3n，n为正整数，L＝(2/3)×N，当N＝3n+1，n为正整数，L＝(2/3)×(N-1)，当N＝3n+2，n为正整数，L＝(2/3)×(N-2)，电流瞬时值时间序列为IN＝(i(1),i(2),…i(L),…i(N))，第一矩阵的每一行为一个电流瞬时值时间序列，
第二矩阵表示为：
其中，N为时间序列长度，当N＝3n，n为正整数，L＝(2/3)×N，当N＝3n+1，n为正整数，L＝(2/3)×(N-1)，当N＝3n+2，n为正整数，L＝(2/3)×(N-2)，风电组机端电压瞬时值时间序列为UN＝(u(1),u(2),Lu(L),Lu(N))，第二矩阵的每一行为一个风电组机端电压瞬时值时间序列。
可选地，在本申请的一个实施例中，将第一矩阵奇异值分解、将第二矩阵奇异值分解表示为：
DI(N-L)(L+1)和DU(N-L)(L+1)是对角阵，移去和的第一列和最后一列，分别形成矩阵和
生成的第一目标矩阵表示为：
生成的第二目标矩阵表示为：
第三目标矩阵表示为：
第四目标矩阵表示为：
其中，K、D、M矩阵为奇异值分解后的矩阵，并进行数据变化，没有具体含义。
可选地，在本申请的一个实施例中，单台风电机组阻抗特性曲线表示为：
其中，P表示电流和电压中含有的振荡分量模态个数，表示第w个模态电压拟合值的幅值，w＝1…P，表示第w个模态电流拟合值的幅值，ωP表示模态频率， 分别表示第w个模态电压和模态电流拟合值的相位，eαpt表示第P个模态随时间衰减的参数，分别由第一中间矩阵、第二中间矩阵生成，
第一中间矩阵表示为：
其中，表示第w个模态电流拟合值的幅值，表示第w个模态电流拟合值的相位，i(N)表示电流瞬时值时间序列中的第N个，
第二中间矩阵表示为：
其中，表示第w个模态电压拟合值的幅值，w＝1…P，表示第w个模态电压，u(N)表示风电组机端电压瞬时值时间序列中的第N个，
为特征值，为特征值,对降序排列，形成对降序排列，形成取令其中，为实部，表示第P个模态分量的衰减率，为虚部，表示第P个模态分量的角频率。
可选地，在本申请的一个实施例中，多机风电机组由若干单台风电机组组成，使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数，具体为对多机风电机组包含的各台风电机组进行降序排列，表示为：
其中，对于含有J台风电机组汇集的风电场，其阻抗特性曲线为：ZJ(t)，对各台风电机组ωP_J进行降序排列，J表示含有J台风电机组汇集的风电场，表示为J台风机第P个模态分量的角频率的等效值，表示为J台风机第P个模态分量的衰减率的等效值，αP_j表示第j台风机的第P个模态分量的衰减率。
可选地，在本申请的一个实施例中，风电场等值阻抗特性曲线表示为：
其中，P表示电压及电流中含有的振荡模态分量个数，J表示风电机组台数，t表示时间，w表示第w个，w＝1,…P，表示表示第j台风机的w个模态分量的电压幅值，表示第j台风机的w个模态分量的电流幅值，和分别表示J台风机第w振荡模态频率和衰减率的拟合值，表示分别表示第j台风机第w个模态的相位电压和相位电流拟合值。
为达上述目的，本申请第二方面实施例提出了一种基于实时数据的直驱风电场等值建模装置，包括：单台风电机组阻抗特性曲线模块、多机风电机组阻抗特性参数优化模块、风电场等值阻抗特性曲线计算模块，其中，
单台风电机组阻抗特性曲线模块，用于获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；
多机风电机组阻抗特性参数优化模块，用于使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数；
风电场等值阻抗特性曲线计算模块，用于使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。
为了实现上述目的，本申请第三方面实施例提出了一种非临时性计算机可读存储介质，当存储介质中的指令由处理器被执行时，能够执行基于实时数据的直驱风电场等值建模方法。
本申请实施例的基于实时数据的直驱风电场等值建模方法、基于实时数据的直驱风电场等值建模装置和非临时性计算机可读存储介质，解决了现有方法对风电场的整体特性研究较少的问题，高效、精确地获取了风电系统在宽频域范围的阻抗特性，为大规模风电机组接入电网的影响分析提供了重要依据。
本申请附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本申请的实践了解到。
附图说明
本申请上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
图1为本申请实施例一所提供的一种基于实时数据的直驱风电场等值建模方法的流程图；
图2为本申请实施例二所提供的一种基于实时数据的直驱风电场等值建模装置的结构示意图。
具体实施方式
下面详细描述本申请的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本申请，而不能理解为对本申请的限制。
下面参考附图描述本申请实施例的基于实时数据的直驱风电场等值建模方法和装置。
图1为本申请实施例一所提供的一种基于实时数据的直驱风电场等值建模方法的流程图。
如图1所示，该基于实时数据的直驱风电场等值建模方法包括以下步骤：
步骤101，获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；
步骤102，使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性曲线的参数；
步骤103，使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。
本申请实施例的基于实时数据的直驱风电场等值建模方法，通过获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性曲线的参数；使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。由此，能够解决现有方法对风电场的整体特性研究较少的问题，高效、精确地获取了风电系统在宽频域范围的阻抗特性，为大规模风电机组接入电网的影响分析提供了重要依据。
进一步地，在本申请实施例中，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线，包括以下步骤：
根据电流瞬时值时间序列和风电组机端电压瞬时值时间序列分别生成第一矩阵和第二矩阵，其中，第一矩阵为电流瞬时值时间序列矩阵，第二矩阵为风电组机端电压瞬时值时间序列矩阵；
将第一矩阵分解生成第一目标矩阵和第二目标矩阵，将第二矩阵分解生成第三目标矩阵和第四目标矩阵；
对第一目标矩阵和第二目标矩阵、第三目标矩阵和第四目标矩阵进行计算生成单台风电机组阻抗特性曲线。
进一步地，在本申请实施例中，第一矩阵表示为：
其中，N为时间序列长度，当N＝3n，n为正整数，L＝(2/3)×N，当N＝3n+1，n为正整数，L＝(2/3)×(N-1)，当N＝3n+2，n为正整数，L＝(2/3)×(N-2)，电流瞬时值时间序列为IN＝(i(1),i(2),…i(L),…i(N))，第一矩阵的每一行为一个电流瞬时值时间序列，
第二矩阵表示为：
其中，N为时间序列长度，当N＝3n，n为正整数，L＝(2/3)×N，当N＝3n+1，n为正整数，L＝(2/3)×(N-1)，当N＝3n+2，n为正整数，L＝(2/3)×(N-2)，风电组机端电压瞬时值时间序列为UN＝(u(1),u(2),Lu(L),Lu(N))，第二矩阵的每一行为一个风电组机端电压瞬时值时间序列。
进一步地，在本申请实施例中，将第一矩阵奇异值分解、将第二矩阵奇异值分解表示为：
DI(N-L)(L+1)和DU(N-L)(L+1)是对角阵，移去和的第一列和最后一列，分别形成矩阵和
生成的第一目标矩阵表示为：
生成的第二目标矩阵表示为：
第三目标矩阵表示为：
第四目标矩阵表示为：
其中，K、D、M矩阵为奇异值分解后的矩阵，并进行数据变化，没有具体含义。
进一步地，在本申请实施例中，单台风电机组阻抗特性曲线表示为：
其中，P表示电流和电压中含有的振荡分量模态个数，表示第w个模态电压拟合值的幅值，w＝1…P，表示第w个模态电流拟合值的幅值，ωP表示模态频率， 分别表示第w个模态电压和模态电流拟合值的相位，eαpt表示第P个模态随时间衰减的参数，分别由第一中间矩阵、第二中间矩阵生成，
第一中间矩阵表示为：
其中，表示第w个模态电流拟合值的幅值，表示第w个模态电流拟合值的相位，i(N)表示电流瞬时值时间序列中的第N个，
第二中间矩阵表示为：
其中，表示第w个模态电压拟合值的幅值，w＝1…P，表示第w个模态电压，u(N)表示风电组机端电压瞬时值时间序列中的第N个，
为特征值，为特征值,对降序排列，形成对降序排列，形成取令其中，为实部，表示第P个模态分量的衰减率，为虚部，表示第P个模态分量的角频率。
进一步地，在本申请实施例中，多机风电机组由若干单台风电机组组成，使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数，具体为对多机风电机组包含的各台风电机组进行降序排列，表示为：
其中，对于含有J台风电机组汇集的风电场，其阻抗特性曲线为：ZJ(t)，对各台风电机组ωP_J进行降序排列，J表示含有J台风电机组汇集的风电场，表示为J台风机第P个模态分量的角频率的等效值，表示为J台风机第P个模态分量的衰减率的等效值，αP_j表示第j台风机的第P个模态分量的衰减率。
对于含有J台风电机组汇集的风电场，其阻抗特性曲线为：zJ(t)。
进一步地，在本申请实施例中，风电场等值阻抗特性曲线表示为：
其中，P表示电压及电流中含有的振荡模态分量个数，J表示风电机组台数，t表示时间，w表示第w个，w＝1,…P，表示表示第j台风机的w个模态分量的电压幅值，表示第j台风机的w个模态分量的电流幅值，和分别表示J台风机第w振荡模态频率和衰减率的拟合值，表示分别表示第j台风机第w个模态的相位电压和相位电流拟合值。
图2为本申请实施例二所提供的一种基于实时数据的直驱风电场等值建模装置的结构示意图。
如图2所示，该基于实时数据的直驱风电场等值建模装置，包括：单台风电机组阻抗特性曲线模块10、多机风电机组阻抗特性参数优化模块20、风电场等值阻抗特性曲线计算模块30，其中，
单台风电机组阻抗特性曲线模块10，用于获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；
多机风电机组阻抗特性参数优化模块20，用于使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数；
风电场等值阻抗特性曲线计算模块30，用于使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。
本申请实施例的基于实时数据的直驱风电场等值建模装置，包括：单台风电机组阻抗特性曲线模块、多机风电机组阻抗特性参数优化模块、风电场等值阻抗特性曲线计算模块，其中，单台风电机组阻抗特性曲线模块，用于获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据风电机组机端电压瞬时值时间序列和电流瞬时值时间序列计算单台风电机组阻抗特性曲线；多机风电机组阻抗特性参数优化模块，用于使用单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数；风电场等值阻抗特性曲线计算模块，用于使用多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成风电场。由此，能够解决现有方法对风电场的整体特性研究较少的问题，高效、精确地获取了风电系统在宽频域范围的阻抗特性，为大规模风电机组接入电网的影响分析提供了重要依据。
为了实现上述实施例，本申请还提出了一种非临时性计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述实施例的基于实时数据的直驱风电场等值建模方法。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本申请的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本申请的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(RAM)，只读存储器(ROM)，可擦除可编辑只读存储器(EPROM或闪速存储器)，光纤装置，以及便携式光盘只读存储器(CDROM)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
应当理解，本申请的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如，如果用硬件来实现和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(PGA)，现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
此外，在本申请各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器，磁盘或光盘等。尽管上面已经示出和描述了本申请的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本申请的限制，本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (6)

1.一种基于实时数据的直驱风电场等值建模方法，其特征在于，包括以下步骤：
获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据所述风电机组机端电压瞬时值时间序列和所述电流瞬时值时间序列计算单台风电机组阻抗特性曲线；
使用所述单台风电机组阻抗特性曲线优化多机风电机组阻抗特性曲线的参数；
使用所述多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成所述风电场；
其中，所述根据所述风电机组机端电压瞬时值时间序列和所述电流瞬时值时间序列计算单台风电机组阻抗特性曲线，包括以下步骤：
根据所述电流瞬时值时间序列和所述风电组机端电压瞬时值时间序列分别生成第一矩阵和第二矩阵，其中，所述第一矩阵为电流瞬时值时间序列矩阵，所述第二矩阵为风电组机端电压瞬时值时间序列矩阵；
将所述第一矩阵分解生成第一目标矩阵和第二目标矩阵，将所述第二矩阵分解生成第三目标矩阵和第四目标矩阵；
对所述第一目标矩阵和第二目标矩阵、所述第三目标矩阵和第四目标矩阵进行计算生成所述单台风电机组阻抗特性曲线；
所述多机风电机组由若干单台风电机组组成，所述使用所述单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数，具体为对所述多机风电机组包含的各台风电机组进行降序排列，表示为：
其中，对于含有J台风电机组汇集的风电场，其阻抗特性曲线为：ZJ(t)，对各台风电机组ωP_J进行降序排列，J表示含有J台风电机组汇集的风电场，表示为J台风机第P个模态分量的角频率的等效值，表示为J台风机第P个模态分量的衰减率的等效值，αP_j表示第j台风机的第P个模态分量的衰减率；
所述风电场等值阻抗特性曲线表示为：
其中，P表示电压及电流中含有的振荡模态分量个数，J表示风电机组台数，t表示时间，w表示第w个，w＝1,…P，表示表示第j台风机的w个模态分量的电压幅值，表示第j台风机的w个模态分量的电流幅值，和分别表示J台风机第w振荡模态频率和衰减率的拟合值，表示分别表示第j台风机第w个模态的相位电压和相位电流拟合值。
2.如权利要求1所述的基于实时数据的直驱风电场等值建模方法，其特征在于，所述第一矩阵表示为：
其中，N为时间序列长度，当N＝3n，n为正整数，L＝(2/3)×N，当N＝3n+1，n为正整数，L＝(2/3)×(N-1)，当N＝3n+2，n为正整数，L＝(2/3)×(N-2)，所述电流瞬时值时间序列为IN＝(i(1),i(2),…i(L),…i(N))，所述第一矩阵的每一行为一个电流瞬时值时间序列，
所述第二矩阵表示为：
其中，N为时间序列长度，当N＝3n，n为正整数，L＝(2/3)×N，当N＝3n+1，n为正整数，L＝(2/3)×(N-1)，当N＝3n+2，n为正整数，L＝(2/3)×(N-2)，所述风电组机端电压瞬时值时间序列为UN＝(u(1),u(2),Lu(L),Lu(N))，所述第二矩阵的每一行为一个风电组机端电压瞬时值时间序列。
3.如权利要求1所述的基于实时数据的直驱风电场等值建模方法，其特征在于，将所述第一矩阵奇异值分解、所述将所述第二矩阵奇异值分解表示为：
DI(N-L)(L+1)和DU(N-L)(L+1)是对角阵，移去和的第一列和最后一列，分别形成矩阵和
生成的所述第一目标矩阵表示为：
生成的所述第二目标矩阵表示为：
所述第三目标矩阵表示为：
所述第四目标矩阵表示为：
其中，K、D、M矩阵为奇异值分解后的矩阵，并进行数据变化，没有具体含义。
4.如权利要求3所述的基于实时数据的直驱风电场等值建模方法，其特征在于，所述单台风电机组阻抗特性曲线表示为：
其中，P表示电流和电压中含有的振荡分量模态个数，表示第w个模态电压拟合值的幅值，w＝1…P，表示第w个模态电流拟合值的幅值，ωP表示模态频率， 分别表示第w个模态电压和模态电流拟合值的相位，eαpt表示第P个模态随时间衰减的参数，分别由第一中间矩阵、第二中间矩阵生成，
所述第一中间矩阵表示为：
其中，表示第w个模态电流拟合值的幅值，表示第w个模态电流拟合值的相位，i(N)表示电流瞬时值时间序列中的第N个，
所述第二中间矩阵表示为：
其中，表示第w个模态电压拟合值的幅值，w＝1…P，表示第w个模态电压，u(N)表示风电组机端电压瞬时值时间序列中的第N个，
为特征值，为特征值,对降序排列，形成对降序排列，形成取令其中，为实部，表示第P个模态分量的衰减率，为虚部，表示第P个模态分量的角频率。
5.一种基于实时数据的直驱风电场等值建模装置，其特征在于，包括单台风电机组阻抗特性曲线模块、多机风电机组阻抗特性参数优化模块、风电场等值阻抗特性曲线计算模块，其中，
所述单台风电机组阻抗特性曲线模块，用于获取风电机组机端电压瞬时值时间序列和电流瞬时值时间序列，根据所述风电机组机端电压瞬时值时间序列和所述电流瞬时值时间序列计算单台风电机组阻抗特性曲线；
所述多机风电机组阻抗特性参数优化模块，用于使用所述单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数；
所述风电场等值阻抗特性曲线计算模块，用于使用所述多机风电机组阻抗特性曲线计算风电场等值阻抗特性曲线，其中，若干多机风电机组汇集组成所述风电场；
其中，其中，所述根据所述风电机组机端电压瞬时值时间序列和所述电流瞬时值时间序列计算单台风电机组阻抗特性曲线，包括以下步骤：
根据所述电流瞬时值时间序列和所述风电组机端电压瞬时值时间序列分别生成第一矩阵和第二矩阵，其中，所述第一矩阵为电流瞬时值时间序列矩阵，所述第二矩阵为风电组机端电压瞬时值时间序列矩阵；
将所述第一矩阵分解生成第一目标矩阵和第二目标矩阵，将所述第二矩阵分解生成第三目标矩阵和第四目标矩阵；
对所述第一目标矩阵和第二目标矩阵、所述第三目标矩阵和第四目标矩阵进行计算生成所述单台风电机组阻抗特性曲线；
所述多机风电机组由若干单台风电机组组成，所述使用所述单台风电机组阻抗特性曲线优化多机风电机组阻抗特性参数，具体为对所述多机风电机组包含的各台风电机组进行降序排列，表示为：
其中，对于含有J台风电机组汇集的风电场，其阻抗特性曲线为：ZJ(t)，对各台风电机组ωP_J进行降序排列，J表示含有J台风电机组汇集的风电场，表示为J台风机第P个模态分量的角频率的等效值，表示为J台风机第P个模态分量的衰减率的等效值，αP_j表示第j台风机的第P个模态分量的衰减率；
所述风电场等值阻抗特性曲线表示为：
其中，P表示电压及电流中含有的振荡模态分量个数，J表示风电机组台数，t表示时间，w表示第w个，w＝1,…P，表示表示第j台风机的w个模态分量的电压幅值，表示第j台风机的w个模态分量的电流幅值，和分别表示J台风机第w振荡模态频率和衰减率的拟合值，表示分别表示第j台风机第w个模态的相位电压和相位电流拟合值。
6.一种非临时性计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-4中任一所述的基于实时数据的直驱风电场等值建模方法。
CN202111220530.3A 2021-10-20 2021-10-20 基于实时数据的直驱风电场等值建模方法和装置 Active CN114154293B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111220530.3A CN114154293B (zh) 2021-10-20 2021-10-20 基于实时数据的直驱风电场等值建模方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111220530.3A CN114154293B (zh) 2021-10-20 2021-10-20 基于实时数据的直驱风电场等值建模方法和装置

Publications (2)

Publication Number Publication Date
CN114154293A CN114154293A (zh) 2022-03-08
CN114154293B true CN114154293B (zh) 2024-07-16

Family

ID=80462472

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111220530.3A Active CN114154293B (zh) 2021-10-20 2021-10-20 基于实时数据的直驱风电场等值建模方法和装置

Country Status (1)

Country Link
CN (1) CN114154293B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118332321B (zh) * 2024-06-13 2024-09-03 山东大学 一种基于事件数据驱动的风电场阻抗获取方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107769227A (zh) * 2017-10-27 2018-03-06 江苏大学 一种适用于次同步研究的风电场等值建模方法
CN109188188A (zh) * 2018-07-10 2019-01-11 国网浙江省电力有限公司杭州供电公司 基于电压监控的多端柔性直流输电线路故障单端判别方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102882229B (zh) * 2012-09-21 2015-06-17 北京金风科创风电设备有限公司 风电场动态电压自动控制系统
CN106786788B (zh) * 2015-11-20 2019-07-19 中国电力科学研究院 一种计及风电场脱网约束的电网风电送出能力评估方法
CN109521693B (zh) * 2018-12-03 2023-05-16 中国电力科学研究院有限公司 一种风电场电磁暂态实时仿真方法和装置
CN111611696A (zh) * 2020-05-12 2020-09-01 北方工业大学 微电网系统的非线性建模方法
CN113098022B (zh) * 2021-04-09 2024-06-25 国网陕西省电力公司电力科学研究院 风电场并网点无功调节方法、装置、设备及存储介质

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107769227A (zh) * 2017-10-27 2018-03-06 江苏大学 一种适用于次同步研究的风电场等值建模方法
CN109188188A (zh) * 2018-07-10 2019-01-11 国网浙江省电力有限公司杭州供电公司 基于电压监控的多端柔性直流输电线路故障单端判别方法

Also Published As

Publication number Publication date
CN114154293A (zh) 2022-03-08

Similar Documents

Publication Publication Date Title
Rouibah et al. A low-cost monitoring system for maximum power point of a photovoltaic system using IoT technique
CN1254004C (zh) 检测交流电动机定子匝故障的方法和装置
Liu et al. Theoretical analysis on the short-circuit current of inverter-interfaced renewable energy generators with fault-ride-through capability
CN114154293B (zh) 基于实时数据的直驱风电场等值建模方法和装置
Muntean et al. Hardware in the loop wind turbine emulator
Hackl et al. Full-and reduced-order state-space modeling of wind turbine systems with permanent magnet synchronous generator
CN117313293B (zh) 一种直驱风电场小信号等值建模方法、系统、终端及介质
CN116760071B (zh) 一种高压箱储能方法及系统
CN110867850B (zh) 一种发电机转速及风电机组参数计算方法、风电机组模型
CN109782068B (zh) 基于实时数字仿真器的宽频阻抗扫描方法及系统
CN116306236A (zh) 一种基于gru-lstm-fc组合网络的风电场动态等值建模方法及系统
Sun et al. Design and optimization of control parameters based on direct-drive permanent magnet synchronous generator for wind power system
CN109904863A (zh) 一种附加阻尼控制器、水轮机调速器及水轮机控制系统
Selim et al. A comprehensive modeling for wind turbine generator and DSTATCOM into backward/forward load flow algorithm
Nadour et al. Advanced backstepping control of a wind energy conversion system using a doubly-fed induction generator
CN115589031A (zh) 永磁直驱风机构网型控制方法、装置、终端及存储介质
CN114154291A (zh) 基于实时数据的半直驱风电场等值建模方法和装置
CN210015441U (zh) 适用于adpss双馈风机闭环试验的并行异构仿真平台
Katić et al. Wavelet transform for voltage dips detection in a microgrid with distributed generation
Chen et al. Nonlinear maximum power point tracking control of wind turbine based on two-mass model without anemometer
Sahin et al. Optimized fractional order control of a cascaded synchronous buck–boost converter for a wave-uc hybrid energy system
Nguyen et al. Implementation and Validation for Multitasks of a Cost-Effective Scheme Based on ESS and Braking Resistors in PMSG Wind Turbine Systems
Vargas et al. Reduced-order extended harmonic domain modeling of switched networks
CN117494591B (zh) 一种风电电能质量评估方法及装置
CN113777403A (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