CN111725799A - 一种基于坐标变换的阻抗计算方法和系统 - Google Patents

一种基于坐标变换的阻抗计算方法和系统 Download PDF

Info

Publication number
CN111725799A
CN111725799A CN201910204280.0A CN201910204280A CN111725799A CN 111725799 A CN111725799 A CN 111725799A CN 201910204280 A CN201910204280 A CN 201910204280A CN 111725799 A CN111725799 A CN 111725799A
Authority
CN
China
Prior art keywords
coordinate system
synchronous generator
virtual synchronous
initial
matrix
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.)
Granted
Application number
CN201910204280.0A
Other languages
English (en)
Other versions
CN111725799B (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 Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201910204280.0A priority Critical patent/CN111725799B/zh
Publication of CN111725799A publication Critical patent/CN111725799A/zh
Application granted granted Critical
Publication of CN111725799B publication Critical patent/CN111725799B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/40—Synchronising a generator for connection to a network or to another generator

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明涉及一种基于坐标变换的阻抗计算方法和系统，包括：根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。本发明提出的阻抗计算方法精准的计算了受扰动后的弱电网线路的阻抗。

Description

一种基于坐标变换的阻抗计算方法和系统
技术领域
本发明涉及逆变器建模领域，具体涉及一种基于坐标变换的阻抗计算方法和系统。
背景技术
新能源发电大量分布于戈壁、山区和海岛等偏远地区，其与主网间的连接强度较弱，电网通常呈现高感抗弱电网特征。新能源并网逆变器接入弱电网易发生谐波振荡等交互稳定性问题。虚拟同步发电机模拟同步发电机的外特性，为大规模新能源友好并网提供新的技术手段，成为当前新能源发电技术的热点。
虚拟同步技术改变传统逆变器控制策略，在控制系统中引入虚拟惯量和阻尼，使得逆变器的输出外特性具有同步发电机类似的旋转惯量、下垂特性和阻尼特性，可为电网提供必要的频率和电压支撑。
由于虚拟同步发电机具备这些优异的性能，虚拟同步技术被推广应用到交-直流微电网的母线电压惯性增强控制、电动汽车快充接口的惯性和阻尼控制以及风力发电变换器和柔性直流输电变换器的惯性控制等领域。
虚拟同步发电机能够主动支撑弱电网的前提是其能够稳定运行，但虚拟同步发电机运行在失去大电网支撑的离网模式时，电网为本地负载供电时会面临许多问题：
一方面，光伏、风电等可再生微源出力的随机性、波动性和负载扰动的不可预测性等因素较大，其会给系统的稳定性控制增加了一定的难度；
另一方面，如果逆变器的参数设计不恰当，将其并入系统时也会引起与同步机并入大电网的过程中发生类似的功率振荡现象，此功率振荡同样严重影响系统的稳定性。
目前，针对虚拟同步发电机的稳定性研究，主流方法是对逆变器系统进行小信号建模，根据Nyquist判据判断稳定性，或是建立系统的阻抗模型，根据阻抗判据判断稳定性。
就目前而言，尚未有相关科学技术手段利用坐标变换对逆变器系统进行小信号建模，通过考虑坐标变换对虚拟同步发电机阻抗计算的影响来减小阻抗计算误差对虚拟同步发电机的稳定性的影响。
发明内容
针对现有技术的不足，本发明的目的是提供一种能精准的计算扰动对弱电网线路阻抗影响的阻抗计算方法，该方法利用坐标变换计算虚拟同步发电机的阻抗，尽可能的减小了阻抗计算误差对虚拟同步发电机的稳定性的影响。
本发明的目的是采用下述技术方案实现的：
本发明提出一种基于坐标变换的阻抗计算方法，其改进之处在于，所述方法包括：
根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。
优选的，所述根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵，包括：
按下式确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure BDA0001998461530000021
Figure BDA0001998461530000022
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
Figure BDA0001998461530000023
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000024
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵。
优选的，所述利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵；
根据所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
进一步的，所述利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵，包括：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA0001998461530000031
Figure BDA0001998461530000032
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000033
式中，
Figure BDA0001998461530000034
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000035
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000036
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；I2为二阶单位矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA0001998461530000037
Figure BDA0001998461530000038
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000039
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
进一步的，所述根据等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000041
Figure BDA0001998461530000042
式中，B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵。
优选的，所述将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000043
Figure BDA0001998461530000044
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000045
式中，
Figure BDA0001998461530000046
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000047
为等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
Figure BDA0001998461530000048
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000049
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Dd为初始dq坐标系下d轴占空比的稳态值；Dq为初始dq坐标系下d轴占空比的稳态值；F为占空比电压系数矩阵；H为占空比电流系数矩阵。
进一步的，利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵，包括：
将初始坐标系下的虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵代入虚拟同步发电机的主电路小信号模型，获得所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，如
下式：Z＝(UdcF+I2)-1(A-UdcH)
式中；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；I2为二阶单位矩阵。
优选的，按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型：
Figure BDA0001998461530000051
式中，
Figure BDA0001998461530000052
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000053
为初始dq坐标系下的占空比扰动矩阵；
Figure BDA0001998461530000054
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的主电路系数矩阵A：
Figure BDA0001998461530000055
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
进一步的，按下式确定所述占空比电压系数矩阵F：
Figure BDA0001998461530000056
式中，Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；
按下式确定所述占空比电流系数矩阵H：
Figure BDA0001998461530000061
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵；E为角度扰动的电压系数矩阵。
进一步的，按下式确定所述虚拟同步发电机的双环控制的电流系数矩阵B；
Figure BDA0001998461530000062
式中，D1为电压调节系数；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；h1为电压环PI控制器取值、h2为电流环PI控制器取值；ω为虚拟同步发电机的同步旋转角速度；Ln为LC滤波器的电感；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的双环控制的电压系数矩阵C；
Figure BDA0001998461530000063
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；C为LC滤波器的电容；
其中，按下式确定所述电压环PI控制器取值h1：
Figure BDA0001998461530000064
式中，kp1为电压环PI控制器比例系数；ki1为电压环PI控制器积分系数；
按下式确定所述电流环PI控制器取值h2：
Figure BDA0001998461530000071
式中，kp2为电流环PI控制器比例系数；kp2为电流环PI控制器积分系数。
进一步的，按下式确定所述角度扰动的电流系数矩阵矩阵D；
Figure BDA0001998461530000072
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
按下式确定所述角度扰动的电压系数矩阵E；
Figure BDA0001998461530000073
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
本发明提出的一种基于坐标变换的阻抗计算系统，其改进之处在于，所述系统包括：
第一确定模块：用于根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
第二确定模块：用于利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
转换模块：用于将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
第三确定模块：用于利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。
优选的，所述第一确定模块，用于：
按下式确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure BDA0001998461530000074
Figure BDA0001998461530000081
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
Figure BDA0001998461530000082
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000083
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵。
优选的，所述第二确定模块，包括：
第一确定子单元：用于利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵；
第二确定子单元：用于根据所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
进一步的，所述第一确定子单元，用于：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA0001998461530000084
Figure BDA0001998461530000085
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000086
式中，
Figure BDA0001998461530000087
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000091
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000092
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；I2为二阶单位矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA0001998461530000093
Figure BDA0001998461530000094
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000095
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
进一步的，所述第二确定子单元，用于：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000096
Figure BDA0001998461530000097
式中，B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵。
优选的，所述转换模块，用于：
按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000098
Figure BDA0001998461530000101
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000102
式中，
Figure BDA0001998461530000103
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000104
为等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
Figure BDA0001998461530000105
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000106
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Dd为初始dq坐标系下d轴占空比的稳态值；Dq为初始dq坐标系下d轴占空比的稳态值；F为占空比电压系数矩阵；H为占空比电流系数矩阵。
进一步的，所述第三确定模块，用于：
将初始坐标系下的虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵代入虚拟同步发电机的主电路小信号模型，获得所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，如
下式：Z＝(UdcF+I2)-1(A-UdcH)
式中；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；I2为二阶单位矩阵。
优选的，按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型：
Figure BDA0001998461530000107
式中，
Figure BDA0001998461530000108
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000109
为初始dq坐标系下的占空比扰动矩阵；
Figure BDA00019984615300001010
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的主电路系数矩阵A：
Figure BDA0001998461530000111
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
进一步的，按下式确定所述占空比电压系数矩阵F：
Figure BDA0001998461530000112
式中，Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；
按下式确定所述占空比电流系数矩阵H：
Figure BDA0001998461530000113
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵；E为角度扰动的电压系数矩阵。
进一步的，按下式确定所述虚拟同步发电机的双环控制的电流系数矩阵B；
Figure BDA0001998461530000114
式中，D1为电压调节系数；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；h1为电压环PI控制器取值、h2为电流环PI控制器取值；ω为虚拟同步发电机的同步旋转角速度；Ln为LC滤波器的电感；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的双环控制的电压系数矩阵C；
Figure BDA0001998461530000121
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；C为LC滤波器的电容；
其中，按下式确定所述电压环PI控制器取值h1：
Figure BDA0001998461530000122
式中，kp1为电压环PI控制器比例系数；ki1为电压环PI控制器积分系数；
按下式确定所述电流环PI控制器取值h2：
Figure BDA0001998461530000123
式中，kp2为电流环PI控制器比例系数；kp2为电流环PI控制器积分系数。
进一步的，按下式确定所述角度扰动的电流系数矩阵矩阵D；
Figure BDA0001998461530000124
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
按下式确定所述角度扰动的电压系数矩阵E；
Figure BDA0001998461530000125
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
与最接近的现有技术相比，本发明具有的有益效果：
本发明提供的技术方案，根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵，本发明提出的阻抗计算方法计算简单，易于操作，适用于离线仿真、在线实时仿真和硬件在环仿真系统，能更精准的计算虚拟同步发电机受扰动后的输出阻抗。
附图说明
图1是本发明提供的基于坐标变换的阻抗计算方法流程图；
图2是本发明实施例1波形图；
图3是本发明实施例2波形图；
图4是本发明实施例3波形图；
图5是本发明实施例4波形图；
图6是本发明提供的基于坐标变换的阻抗计算系统流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
本发明提供了一种基于坐标变换的阻抗计算方法，如图1所示，所述包括：
101.根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；虚拟同步发电机受扰动后在初始dq坐标系下的参量包括：初始dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和初始dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵；
102.利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
103.将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
104.利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。
具体的，所述步骤101，包括：
按下式确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure BDA0001998461530000141
Figure BDA0001998461530000142
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；s为拉普拉斯算子；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
Figure BDA0001998461530000143
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000144
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
进一步的，在本发明最优实施例中，确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure BDA0001998461530000145
的过程可以为：
获取虚拟同步发电机同步控制有功-频率控制原理：
Figure BDA0001998461530000146
式中，P＝(edid+eqiq)，J是转子转动惯量；Tm为虚拟同步发电机的机械转矩；Te为虚拟同步发电机的电磁转矩；Pset为虚拟同步发电机的机械功率、P为虚拟同步发电机的有功功率；D'为阻尼系数；id为初始dq坐标系下LC滤波器经过电感的d轴电流；iq为初始dq坐标系下LC滤波器经过电感的q轴电流；ed初始dq坐标系下虚拟同步发电机的d轴输出电压；eq为初始dq坐标系下虚拟同步发电机的q轴输出电压；ω为虚拟同步发电机的同步旋转角速度；ωo为虚拟同步发电机的标准同步旋转角速度；
将所述虚拟同步发电机同步控制有功-频率控制原理转换到S域内，得：
Figure BDA0001998461530000151
将S域内的虚拟同步发电机同步控制有功-频率控制原理代入至公式
Figure BDA0001998461530000152
中，得：
Figure BDA0001998461530000153
在虚拟同步发电机稳态工作状态下施加小扰动，获取虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure BDA0001998461530000154
Figure BDA0001998461530000155
即：
Figure BDA0001998461530000156
式中，
Figure BDA0001998461530000157
Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
进一步的，所述步骤102，包括：
利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵；
根据所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
具体的，所述利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵，包括：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA0001998461530000158
Figure BDA0001998461530000159
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000161
式中，
Figure BDA0001998461530000162
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000163
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000164
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；I2为二阶单位矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
进一步的，在本发明的最优实施例中，获取所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA0001998461530000165
的过程可以为；
根据坐标变换原理，确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA0001998461530000166
Figure BDA0001998461530000167
式中，
Figure BDA0001998461530000168
为等效dq坐标系下LC滤波器经过电感的d轴电流的稳态值；
Figure BDA0001998461530000169
为等效dq坐标系下LC滤波器经过电感的q轴电流的稳态值；
化简所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA00019984615300001610
得：
Figure BDA00019984615300001611
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA00019984615300001612
Figure BDA0001998461530000171
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000172
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
进一步的，在本发明的最优实施例中，获取所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA0001998461530000173
的过程可以为：
根据坐标变换原理，确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA0001998461530000174
Figure BDA0001998461530000175
其中，
Figure BDA0001998461530000176
为等效dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；
Figure BDA0001998461530000177
为等效dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
化简所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA0001998461530000178
得：
Figure BDA0001998461530000179
将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA00019984615300001710
和电流扰动矩阵
Figure BDA00019984615300001711
代入预先获取的虚拟同步发电机的双环控制小信号模型中，确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
具体的，所述根据等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000181
Figure BDA0001998461530000182
式中，B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵。
进一步的，在本发明的最优实施例中，获取虚拟同步发电机的双环控制公式
Figure BDA0001998461530000183
的过程可以为：
获取虚拟同步发电机的双环控制原理：
Figure BDA0001998461530000184
Em＝E0-D1Q，Q＝(eq id-ed iq)
其中，Em是虚拟同步发电机控制电压；E0是虚拟同步发电机的空载电压；D1是电压调节系数；Q是无功功率；
Figure BDA0001998461530000185
为初始dq坐标系下控制器输出d轴电压；
Figure BDA0001998461530000186
为初始dq坐标系下控制器输出q轴电压；
在虚拟同步发电机稳态工作状态下施加小扰动，获取虚拟同步发电机的双环控制公式：
Figure BDA0001998461530000187
取调制度为
Figure BDA0001998461530000188
将上述虚拟同步发电机的双环控制公式改写为矩阵形式，得：
Figure BDA0001998461530000189
式中，
Figure BDA00019984615300001810
进一步的，所述步骤103，包括：
按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000191
Figure BDA0001998461530000192
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000193
式中，
Figure BDA0001998461530000194
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000195
为等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
Figure BDA0001998461530000196
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000197
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Dd为初始dq坐标系下d轴占空比的稳态值；Dq为初始dq坐标系下d轴占空比的稳态值；F为占空比电压系数矩阵；H为占空比电流系数矩阵。
进一步的，在本发明的最优实施例中，确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000198
的过程可以为：
根据坐标变换原理，确定初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000199
Figure BDA00019984615300001910
式中，
Figure BDA00019984615300001911
为等效dq坐标系下d轴占空比的稳态值；
Figure BDA00019984615300001912
为等效dq坐标系下d轴占空比的稳态值；
化简所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA00019984615300001913
得：
Figure BDA0001998461530000201
将所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000202
代入到所述虚拟同步发电机的主电路小信号模型中，确定所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，其中，
Figure BDA0001998461530000203
进一步的，所述步骤104，包括：
将初始坐标系下的虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵代入虚拟同步发电机的主电路小信号模型，获得所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，如下式：
Z＝(UdcF+I2)-1(A-UdcH)
式中；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；I2为二阶单位矩阵。
具体的，按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型：
Figure BDA0001998461530000204
式中，
Figure BDA0001998461530000205
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000206
为初始dq坐标系下的占空比扰动矩阵；
Figure BDA0001998461530000207
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的主电路系数矩阵A：
Figure BDA0001998461530000208
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
例如，在本发明的最优实施例中，获取初始dq坐标系下虚拟同步发电机的主电路小信号模型的过程可以为：
获取虚拟同步发电机的主电路原理：
Figure BDA0001998461530000211
式中，id为初始dq坐标系下LC滤波器经过电感的d轴电流；iq为初始dq坐标系下LC滤波器经过电感的q轴电流；ed初始dq坐标系下虚拟同步发电机的d轴输出电压；eq为初始dq坐标系下虚拟同步发电机的q轴输出电压；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度；dd为初始dq坐标系下的d轴占空比；dq为初始dq坐标系下的q轴占空比；udc是虚拟同步发电机直流侧电压。
在虚拟同步发电机稳态工作状态下施加小扰动，获取虚拟同步发电机的主电路小信号模型：
Figure BDA0001998461530000212
式中，
Figure BDA0001998461530000213
为初始dq坐标系下LC滤波器经过电感的d轴扰动电流；
Figure BDA0001998461530000214
为初始dq坐标系下LC滤波器经过电感的q轴扰动电流；
Figure BDA0001998461530000215
为初始dq坐标系下虚拟同步发电机的d轴扰动输出电压；
Figure BDA0001998461530000216
为初始dq坐标系下虚拟同步发电机的q轴扰动输出电压；
Figure BDA0001998461530000217
为初始dq坐标系下的d轴扰动占空比；
Figure BDA0001998461530000218
为初始dq坐标系下的q轴扰动占空比；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值。
将上述虚拟同步发电机的主电路小信号模型线性化处理，得：
Figure BDA0001998461530000221
将上述线性化处理后的虚拟同步发电机的主电路小信号模型整理成矩阵形式，得：
Figure BDA0001998461530000222
式中，
Figure BDA0001998461530000223
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000224
为初始dq坐标系下的占空比扰动矩阵；
Figure BDA0001998461530000225
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
将矩阵形式的虚拟同步发电机的主电路小信号模型变化到S域内，得：
Figure BDA0001998461530000226
式中，A为虚拟同步发电机的主电路系数矩阵；
其中，所述虚拟同步发电机的主电路系数矩阵A：
Figure BDA0001998461530000227
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
本发明技术方案涉及的矩阵的具体计算过程如下，按下式确定所述占空比电压系数矩阵F：
Figure BDA0001998461530000228
式中，Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；
按下式确定所述占空比电流系数矩阵H：
Figure BDA0001998461530000231
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵；E为角度扰动的电压系数矩阵。
具体的，按下式确定所述虚拟同步发电机的双环控制的电流系数矩阵B；
Figure BDA0001998461530000232
式中，D1为电压调节系数；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；h1为电压环PI控制器取值、h2为电流环PI控制器取值；ω为虚拟同步发电机的同步旋转角速度；Ln为LC滤波器的电感；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的双环控制的电压系数矩阵C；
Figure BDA0001998461530000233
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；C为LC滤波器的电容；
其中，按下式确定所述电压环PI控制器取值h1：
Figure BDA0001998461530000234
式中，kp1为电压环PI控制器比例系数；ki1为电压环PI控制器积分系数；
按下式确定所述电流环PI控制器取值h2：
Figure BDA0001998461530000241
式中，kp2为电流环PI控制器比例系数；kp2为电流环PI控制器积分系数。
具体的，按下式确定所述角度扰动的电流系数矩阵矩阵D；
Figure BDA0001998461530000242
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
按下式确定所述角度扰动的电压系数矩阵E；
Figure BDA0001998461530000243
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
通常，设计弱电网中线路阻抗参数，需要通过稳定性分析验证阻抗的准确性，验证方法为：令系统阻抗回比矩阵为：
Figure BDA0001998461530000244
式中Z是虚拟同步发电机输出阻抗，Zload是包含线路阻抗的负载阻抗。由广义奈奎斯特判据可知，当且仅当阻抗回比矩阵L的每个特征函数的奈奎斯特曲线在s平面逆时针环绕(-1,j0)的圈数，与虚拟同步发电机输出阻抗Z和负载阻抗Zload在右半平面极点总数相同时，系统是稳定的。负载阻抗不存在右半平面极点，由于虚拟同步发电机单独工作稳定，输出阻抗Z也不存在右半平面极点，因此，阻抗回比矩阵L的每个特征函数的奈奎斯特曲线都不环绕(-1,j0)点，系统是稳定的。
例如，对上述方案进行模拟仿真，虚拟同步发电机主电路参数为：直流侧电压Udc＝800V，LC滤波器电感Ln＝3mH，LC滤波器电容C＝20pF，线路阻抗Lg＝0.5mH，恒功负载P＝5000W。
系统控制参数：PWM开关频率为10kHz，电压环PI控制器参数kp1＝1，ki1＝100,电流环PI控制器参数kp2＝1，ki2＝100,空载电动势E0＝381V,电压调节系数Dq＝0.005,转子转动惯量J＝0.5kg.m2，给定机械功率Pset＝5000W，阻尼系数D＝10.4，给定角速度ω0＝100Πrad/s。
本发明的阻抗计算方法在不同线路阻抗下的仿真结果如图2-5所示：
图2为当虚拟同步发电机在稳态时(Lg＝0.5mH)时Nyquist波形图；
图3为当虚拟同步发电机在失稳时(Lg＝2mH)时Nyquist波形图；
图4为本发明阻抗计算方法在MATLAB/Simulink环境下做出的实施例的线路阻抗突变(第0.3秒时，线路阻抗Lg由0.5mH→2mH)工况下仿真结果中虚拟同步发电机的输出电压和电流波形；
图5为本发明阻抗计算方法在MATLAB/Simulink环境下做出的实施例的线路阻抗突变第(0.7秒时，Lg由2mH→0.5mH)工况下仿真结果中虚拟同步发电机的输出电压和电流波形。
以上所述为本发明的一种实施方式，基于这种实施方式可在MATLAB/Simulink下进行本发明的虚拟同步发电机阻抗计算的仿真，基于这种实施方式亦可以在dSPACE或者RT-LAB等类似的实时仿真器中进行虚拟同步发电机运行在不同工况下的仿真。
本发明提供了一种基于坐标变换的阻抗计算系统，如图6所示，所述系统包括：
第一确定模块：用于根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
第二确定模块：用于利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
转换模块：用于将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
第三确定模块：用于利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。
具体的，所述第一确定模块，用于：
按下式确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure BDA0001998461530000251
Figure BDA0001998461530000252
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
Figure BDA0001998461530000261
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000262
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵。
具体的，所述第二确定模块，包括：
第一确定子单元：用于利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵；
第二确定子单元：用于根据所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
具体的，所述第一确定子单元，用于：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure BDA0001998461530000263
Figure BDA0001998461530000264
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000265
式中，
Figure BDA0001998461530000266
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000267
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000268
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；I2为二阶单位矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure BDA0001998461530000271
Figure BDA0001998461530000272
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000273
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
具体的，所述第二确定子单元，用于：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000274
Figure BDA0001998461530000275
式中，B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵。
具体的，所述转换模块，用于：
按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure BDA0001998461530000276
Figure BDA0001998461530000277
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure BDA0001998461530000278
式中，
Figure BDA0001998461530000279
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure BDA0001998461530000281
为等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
Figure BDA0001998461530000282
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000283
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Dd为初始dq坐标系下d轴占空比的稳态值；Dq为初始dq坐标系下d轴占空比的稳态值；F为占空比电压系数矩阵；H为占空比电流系数矩阵。
具体的，所述第三确定模块，用于：
将初始坐标系下的虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵代入虚拟同步发电机的主电路小信号模型，获得所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，如下式：Z＝(UdcF+I2)-1(A-UdcH)
式中；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；I2为二阶单位矩阵。
具体的，按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型：
Figure BDA0001998461530000284
式中，
Figure BDA0001998461530000285
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure BDA0001998461530000286
为初始dq坐标系下的占空比扰动矩阵；
Figure BDA0001998461530000287
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的主电路系数矩阵A：
Figure BDA0001998461530000288
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
具体的，按下式确定所述占空比电压系数矩阵F：
Figure BDA0001998461530000289
式中，Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；
按下式确定所述占空比电流系数矩阵H：
Figure BDA0001998461530000291
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵；E为角度扰动的电压系数矩阵。
具体的，按下式确定所述虚拟同步发电机的双环控制的电流系数矩阵B；
Figure BDA0001998461530000292
式中，D1为电压调节系数；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；h1为电压环PI控制器取值、h2为电流环PI控制器取值；ω为虚拟同步发电机的同步旋转角速度；Ln为LC滤波器的电感；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的双环控制的电压系数矩阵C；
Figure BDA0001998461530000293
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；C为LC滤波器的电容；
其中，按下式确定所述电压环PI控制器取值h1：
Figure BDA0001998461530000301
式中，kp1为电压环PI控制器比例系数；ki1为电压环PI控制器积分系数；
按下式确定所述电流环PI控制器取值h2：
Figure BDA0001998461530000302
式中，kp2为电流环PI控制器比例系数；kp2为电流环PI控制器积分系数。
具体的，按下式确定所述角度扰动的电流系数矩阵矩阵D；
Figure BDA0001998461530000303
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
按下式确定所述角度扰动的电压系数矩阵E；
Figure BDA0001998461530000304
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (22)

1.一种基于坐标变换的阻抗计算方法，其特征在于，所述方法包括：
根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。
2.如权利要求1所述的方法，其特征在于，所述根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵，包括：
按下式确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure FDA0001998461520000011
Figure FDA0001998461520000012
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
Figure FDA0001998461520000013
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000014
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵。
3.如权利要求1所述的方法，其特征在于，所述利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵；
根据所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
4.如权利要求3所述的方法，其特征在于，所述利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵，包括：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure FDA0001998461520000021
Figure FDA0001998461520000022
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure FDA0001998461520000023
式中，
Figure FDA0001998461520000024
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure FDA0001998461520000025
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000026
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；I2为二阶单位矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure FDA0001998461520000027
Figure FDA0001998461520000028
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure FDA0001998461520000031
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
5.如权利要求4所述的方法，其特征在于，所述根据等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure FDA0001998461520000032
Figure FDA0001998461520000033
式中，B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵。
6.如权利要求1所述的方法，其特征在于，所述将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵，包括：
按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure FDA0001998461520000034
Figure FDA0001998461520000035
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure FDA0001998461520000036
式中，
Figure FDA0001998461520000037
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure FDA0001998461520000038
为等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
Figure FDA0001998461520000041
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000042
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Dd为初始dq坐标系下d轴占空比的稳态值；Dq为初始dq坐标系下d轴占空比的稳态值；F为占空比电压系数矩阵；H为占空比电流系数矩阵。
7.如权利要求6所述的方法，其特征在于，利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵，包括：
将初始坐标系下的虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵代入虚拟同步发电机的主电路小信号模型，获得所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，如下式：
Z＝(UdcF+I2)-1(A-UdcH)
式中；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；I2为二阶单位矩阵。
8.如权利要求7所述的方法，其特征在于，按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型：
Figure FDA0001998461520000043
式中，
Figure FDA0001998461520000044
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000045
为初始dq坐标系下的占空比扰动矩阵；
Figure FDA0001998461520000046
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的主电路系数矩阵A：
Figure FDA0001998461520000047
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
9.如权利要求6或7所述的方法，其特征在于，按下式确定所述占空比电压系数矩阵F：
Figure FDA0001998461520000051
式中，Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；
按下式确定所述占空比电流系数矩阵H：
Figure FDA0001998461520000052
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵；E为角度扰动的电压系数矩阵。
10.如权利要求9所述的方法，其特征在于，按下式确定所述虚拟同步发电机的双环控制的电流系数矩阵B；
Figure FDA0001998461520000053
式中，D1为电压调节系数；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；h1为电压环PI控制器取值、h2为电流环PI控制器取值；ω为虚拟同步发电机的同步旋转角速度；Ln为LC滤波器的电感；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的双环控制的电压系数矩阵C；
Figure FDA0001998461520000054
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；C为LC滤波器的电容；
其中，按下式确定所述电压环PI控制器取值h1：
Figure FDA0001998461520000061
式中，kp1为电压环PI控制器比例系数；ki1为电压环PI控制器积分系数；
按下式确定所述电流环PI控制器取值h2：
Figure FDA0001998461520000062
式中，kp2为电流环PI控制器比例系数；kp2为电流环PI控制器积分系数。
11.如权利要求9所述的方法，其特征在于，按下式确定所述角度扰动的电流系数矩阵矩阵D；
Figure FDA0001998461520000063
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
按下式确定所述角度扰动的电压系数矩阵E；
Figure FDA0001998461520000064
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
12.一种基于坐标变换的阻抗计算系统，其特征在于，所述系统包括：
第一确定模块：用于根据虚拟同步发电机受扰动后在初始dq坐标系下的参量确定虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
第二确定模块：用于利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
转换模块：用于将所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵转换为初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
第三确定模块：用于利用所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵确定虚拟同步发电机受扰动后的输出阻抗矩阵。
13.如权利要求12所述的系统，其特征在于，所述第一确定模块，用于：
按下式确定所述虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵
Figure FDA0001998461520000071
Figure FDA0001998461520000072
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
Figure FDA0001998461520000073
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000074
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵。
14.如权利要求12所述的系统，其特征在于，所述第二确定模块，包括：
第一确定子单元：用于利用所述偏差角小信号矩阵确定等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵；
第二确定子单元：用于根据所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵和电压扰动矩阵确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵。
15.如权利要求14所述的系统，其特征在于，所述第一确定子单元，用于：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电流扰动矩阵
Figure FDA0001998461520000075
Figure FDA0001998461520000076
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure FDA0001998461520000081
式中，
Figure FDA0001998461520000082
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure FDA0001998461520000083
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000084
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；I2为二阶单位矩阵；D为角度扰动的电流系数矩阵；E为角度扰动的电压系数矩阵；
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的电压扰动矩阵
Figure FDA0001998461520000085
Figure FDA0001998461520000086
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure FDA0001998461520000087
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
16.如权利要求15所述的系统，其特征在于，所述第二确定子单元，用于：
按下式确定所述等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure FDA0001998461520000088
Figure FDA0001998461520000089
式中，B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵。
17.如权利要求12所述的系统，其特征在于，所述转换模块，用于：
按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵
Figure FDA0001998461520000091
Figure FDA0001998461520000092
将虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵代入上式中，得：
Figure FDA0001998461520000093
式中，
Figure FDA0001998461520000094
为虚拟同步发电机受扰动后的等效dq坐标系与所述初始dq坐标系的偏差角小信号矩阵；
Figure FDA0001998461520000095
为等效dq坐标系下虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵；
Figure FDA0001998461520000096
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000097
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；Dd为初始dq坐标系下d轴占空比的稳态值；Dq为初始dq坐标系下d轴占空比的稳态值；F为占空比电压系数矩阵；H为占空比电流系数矩阵。
18.如权利要求17所述的系统，其特征在于，所述第三确定模块，用于：
将初始坐标系下的虚拟同步发电机的主电路小信号模型中的占空比扰动矩阵代入虚拟同步发电机的主电路小信号模型，获得所述虚拟同步发电机受扰动后的输出阻抗矩阵Z，如下式：
Z＝(UdcF+I2)-1(A-UdcH)
式中；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；I2为二阶单位矩阵。
19.如权利要求12所述的系统，其特征在于，按下式确定所述初始dq坐标系下虚拟同步发电机的主电路小信号模型：
Figure FDA0001998461520000101
式中，
Figure FDA0001998461520000102
为初始dq坐标系下LC滤波器经过电感的电流扰动矩阵；
Figure FDA0001998461520000103
为初始dq坐标系下的占空比扰动矩阵；
Figure FDA0001998461520000104
为初始dq坐标系下虚拟同步发电机的电压扰动矩阵；A为虚拟同步发电机的主电路系数矩阵；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的主电路系数矩阵A：
Figure FDA0001998461520000105
式中，s为拉普拉斯算子；Ln为LC滤波器的电感；Rn为LC滤波器的寄生电阻；ω为虚拟同步发电机的同步旋转角速度。
20.如权利要求17或18所述的系统，其特征在于，按下式确定所述占空比电压系数矩阵F：
Figure FDA0001998461520000106
式中，Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；
按下式确定所述占空比电流系数矩阵H：
Figure FDA0001998461520000107
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；Dd为初始dq坐标系下的d轴占空比的稳态值；Dq为初始dq坐标系下的q轴占空比的稳态值；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；B为虚拟同步发电机的双环控制的电流系数矩阵；C为虚拟同步发电机的双环控制的电压系数矩阵；E为角度扰动的电压系数矩阵。
21.如权利要求20所述的系统，其特征在于，按下式确定所述虚拟同步发电机的双环控制的电流系数矩阵B；
Figure FDA0001998461520000111
式中，D1为电压调节系数；Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值；h1为电压环PI控制器取值、h2为电流环PI控制器取值；ω为虚拟同步发电机的同步旋转角速度；Ln为LC滤波器的电感；Udc为虚拟同步发电机直流侧电压的稳态值；
按下式确定所述虚拟同步发电机的双环控制的电压系数矩阵C；
Figure FDA0001998461520000112
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；C为LC滤波器的电容；
其中，按下式确定所述电压环PI控制器取值h1：
Figure FDA0001998461520000113
式中，kp1为电压环PI控制器比例系数；ki1为电压环PI控制器积分系数；
按下式确定所述电流环PI控制器取值h2：
Figure FDA0001998461520000114
式中，kp2为电流环PI控制器比例系数；kp2为电流环PI控制器积分系数。
22.如权利要求20所述的系统，其特征在于，按下式确定所述角度扰动的电流系数矩阵矩阵D；
Figure FDA0001998461520000115
式中，Id为初始dq坐标系下LC滤波器经过电感的d轴电流的稳态值；Iq为初始dq坐标系下LC滤波器经过电感的q轴电流的稳态值；s为拉普拉斯算子；ωo为虚拟同步发电机的标准同步旋转角速度；D'为阻尼系数；J为虚拟同步发电机的转子转动惯量；
按下式确定所述角度扰动的电压系数矩阵E；
Figure FDA0001998461520000121
式中，Ed为初始dq坐标系下虚拟同步发电机的d轴输出电压的稳态值；Eq为初始dq坐标系下虚拟同步发电机的q轴输出电压的稳态值。
CN201910204280.0A 2019-03-18 2019-03-18 一种基于坐标变换的阻抗计算方法和系统 Active CN111725799B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910204280.0A CN111725799B (zh) 2019-03-18 2019-03-18 一种基于坐标变换的阻抗计算方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910204280.0A CN111725799B (zh) 2019-03-18 2019-03-18 一种基于坐标变换的阻抗计算方法和系统

Publications (2)

Publication Number Publication Date
CN111725799A true CN111725799A (zh) 2020-09-29
CN111725799B CN111725799B (zh) 2023-11-03

Family

ID=72562280

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910204280.0A Active CN111725799B (zh) 2019-03-18 2019-03-18 一种基于坐标变换的阻抗计算方法和系统

Country Status (1)

Country Link
CN (1) CN111725799B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106712100A (zh) * 2017-03-02 2017-05-24 江苏大学 一种基于扰动观测的虚拟并网同步逆变器控制方法及系统
CN107994620A (zh) * 2017-12-28 2018-05-04 东南大学 柔性环网控制器双端虚拟电机控制方法
CN108258735A (zh) * 2018-03-21 2018-07-06 广东电网有限责任公司电力科学研究院 一种在虚拟同步机中引入虚拟阻抗的仿真控制方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106712100A (zh) * 2017-03-02 2017-05-24 江苏大学 一种基于扰动观测的虚拟并网同步逆变器控制方法及系统
CN107994620A (zh) * 2017-12-28 2018-05-04 东南大学 柔性环网控制器双端虚拟电机控制方法
CN108258735A (zh) * 2018-03-21 2018-07-06 广东电网有限责任公司电力科学研究院 一种在虚拟同步机中引入虚拟阻抗的仿真控制方法及系统

Also Published As

Publication number Publication date
CN111725799B (zh) 2023-11-03

Similar Documents

Publication Publication Date Title
CN107968591B (zh) 基于预测公共点电压的并网逆变器虚拟惯性功率解耦控制方法
CN111541274B (zh) 一种基于虚拟同步发电机特性的孤岛微电网控制方法
CN107732978B (zh) 一种基于虚拟同步发电机改进功率分配策略的逆变器系统
CN106356884B (zh) 一种基于虚拟同步机的光伏并网控制方法、装置及系统
CN110504711B (zh) 一种基于新能源同步机的新能源并网控制系统和方法
Adzic et al. Maximum power search in wind turbine based on fuzzy logic control
CN111193291B (zh) 一种适用于不平衡条件下复合型虚拟同步机控制方法
CN110739678B (zh) 一种并网换流器串联虚拟阻抗的控制方法
CN108599264B (zh) 一种基于虚拟同步发电机控制的频率电压无差调节方法
CN107565599B (zh) 一种基于vsg的风电并网半实物仿真系统
CN115313419B (zh) 一种基于构网型变流器的直驱风电系统控制方法及系统
CN109638881B (zh) 电网强度自适应优化的储能逆变器虚拟同步方法及系统
CN107887919B (zh) 一种分布式储能逆变器协调控制方法及系统
CN110768239B (zh) 一种基于p-u下垂特性的虚拟直流电机控制方法
CN110752762B (zh) 一种并网换流器并联虚拟电容的控制方法
CN106712115A (zh) 一种不用电流反馈的虚拟同步发电机控制器
CN111049178A (zh) 一种直驱永磁风电机组经vsc-hvdc并网稳定控制分析方法
Cheng et al. Coordinated derived current control of DFIG’s RSC and GSC without PLL under unbalanced grid voltage conditions
Li et al. A grid frequency support control strategy of the three phase cascaded H-bridge based photovoltaic generation system
CN105099320A (zh) 永磁直驱风力发电机输出有功功率控制方法及装置
Kenne et al. Adaptive PI control strategy for a self-excited induction generator driven by a variable speed wind turbine
CN111725799A (zh) 一种基于坐标变换的阻抗计算方法和系统
Pradhan et al. Adaptive position observer for multimode wind-BES based microgrid interfaced to distribution network
Noussi et al. Towards a global nonlinear control strategy for DFIG-based wind turbine in a high wind energy penetrated system
Yu et al. Fuzzy-PI and feedforward control strategy of DFIG wind turbine

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