CN111259548B - 一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统 - Google Patents
一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统 Download PDFInfo
- Publication number
- CN111259548B CN111259548B CN202010045691.2A CN202010045691A CN111259548B CN 111259548 B CN111259548 B CN 111259548B CN 202010045691 A CN202010045691 A CN 202010045691A CN 111259548 B CN111259548 B CN 111259548B
- Authority
- CN
- China
- Prior art keywords
- model
- frequency shift
- domain
- identification
- time
- 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.)
- Expired - Fee Related
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 29
- 230000001052 transient effect Effects 0.000 title claims abstract description 16
- 230000005540 biological transmission Effects 0.000 claims abstract description 22
- 230000003993 interaction Effects 0.000 claims abstract description 4
- 239000000126 substance Substances 0.000 claims description 8
- 230000004044 response Effects 0.000 abstract description 7
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000000034 method Methods 0.000 description 21
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 239000003990 capacitor Substances 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统,系统包括辨识模型和频移模型,辨识模型用于等效光伏电站,移模型用于等效交流电网,辨识模型和频移模型之间采用多域传输线模型作为接口模型,并反应辨识模型和频移模型之间的宽频带相互作用关系。本申请提供了面向大规模光伏电站及交流电网的基于辨识模型和频移模型的混合仿真系统,在保证所提出的模型动态响应与实际光伏并网系统一致的同时,可以极大程度上简化电磁分析,加快运算收敛速率,增大收敛精度。
Description
技术领域
本发明涉及电力系统电磁暂态分析技术领域,尤其涉及大规模光伏并网的电磁暂态仿真系统。
背景技术
在研究大规模光伏电站及交流电网的运行时,随着交流电网规模的增大,节点电压方程中系数矩阵的维数将会急剧增加。且由于光伏电站内部具有数量较多的逆变器,包括详细的控制及保护系统,其复杂的拓扑结构使得对大规模光伏电站的电磁暂态建模及分析变得极其繁琐。
更重要的是,在实际的工程环境中,光伏电站或每一个光伏组件内的详细参数及其动态响应通常是无法获知的。
此外,参见图1,原初始模型中光伏电站的控制系统分为电压外环控制以及电流(有功及无功电流)内环控制。直流电压外环的主要作用在于稳定以及调节直流电压,而电流内环旨在实现对id,iq的无静差控制,从而对逆变器输出有功以及无功功率进行调节。控制结构及流程较为复杂,容易造成过程误差。
因此,对于光伏电站及交流电网的运行研究欠缺简化且有效的验证方法。
发明内容
为了克服现有技术的不足,本发明的目的在于提供一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统,其能解决相关问题。
设计原理:一种面向大规模光伏电站及交流电网的基于辨识和频移模型的混合仿真方法,方法拟将仿真系统分为两部分:一、采用频移模型(SFP)来反映交流电网动态特性,二、采用简化后的辨识模型来等效替代光伏电站复杂模型,两者之间的接口特性通过所提出的多域传输线模型(MD-TLM)来体现。通过仿真实验验证该方法的有效性和正确性,为光伏电站及交流电网运行的研究提供一种新的有效途径。
技术方案:本发明的目的采用以下技术方案实现。
一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统,所述仿真系统包括辨识模型和频移模型,所述辨识模型用于等效光伏电站,所述频移模型用于等效交流电网,所述辨识模型和频移模型之间采用多域传输线模型作为接口模型,并反应辨识模型和频移模型之间的宽频带相互作用关系。
所述辨识模型的关键参数包括电容C、电压外环控制的比例环节参数Kp、电流内环控制的时间常数Tp,Tq。
优选的,变量在频移模型中以相量形式表现为:
所述多域传输线模型的建立取决于传输线时延τ与时间步长Δt之间的关系,
当τ>Δt时,多域传输线模型用双层诺顿等效电路来表示为:
当τ<Δt时,在频域下,多域传输线模型通过线性插值的方式来获得频域下电压以及电流相量,即:
其中p,q满足:
且MD-TLM模型表达式为:
将上述结果代入该模型表达式中则有:
相比现有技术,本发明的有益效果在于:面向大规模光伏电站及交流电网基于辨识模型和频移模型的混合仿真系统,在保证所提出的模型动态响应与实际光伏并网系统一致的同时,可以极大程度上简化电磁分析,加快运算收敛速率,增大收敛精度。
附图说明
图1是本发明涉及的现有光伏电站原复杂模型的控制系统结构图;
图2是本发明涉及的光伏并网系统的简化控制辨识模型;
图3是本发明涉及的基于微分进化算法的优化迭代求解流程图;
图4是本发明涉及的多域传输线模型拓扑结构图;
图5是本发明涉及的光伏辨识子系统与频移子系统拓扑结构图;
图6是本发明涉及的光伏辨识与频移子系统的具体实例图;
图7是是本发明涉及的主要实现步骤的流程图;
图8是本发明所提出的模型与实测数据外特性对比结果图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
仿真系统参见图2-图8,具体如下。
(1)给出光伏并网系统的简化控制模型,如图2,是替代了光伏并网系统中控制系统的模型,频移模型体现在交流电网中,而该辨识模型是在光伏逆变器模型中得以体现。两子系统接口特性通过多域传输线体现。
具体而言,则是将电压外环控制中采用P控制器来代替PI控制器,从而减少控制参量,其中P控制器中的比例参量为KP。电流内环控制中等效为一延时环节,有功及无功电流环节中一阶惯性环节中的时间常数分别为Td,Tq。
与简化模型相对应的状态空间方程组表示为:
其中,udc表示直流侧电容两端电压,iPV为光伏电池板等效输出电流,P表示系统中有功功率(不考虑功率损耗),udc_ref、id_ref、iq_ref分别表示直流侧电压udc,电网侧d,q轴电流id、iq参考值,C为直流侧电容,Kp为电压外环控制的比例环节参数,Td,Tq则分别代表电流内环控制的时间常数。
(2)采集原详细仿真模型的输入及动态响应信息,并将根据该输入确定辨识模型中的关键参数,从而使得辨识模型与初始模型的外特性保持一致。将所搭建的辨识模型代替原模型中逆变器的控制系统这一环节,在简化模型的同时亦能正确跟踪并输出功率以及电流等所需电信息。
在简化的辨识模型中,为了使得与原模型动态响应一致,需辨识出四个关键参数的值,即电容C,电压外环控制的比例环节参数Kp,电流内环控制的时间常数Td,Tq。
辨识过程中所涉及到的具体步骤如下:
①:计算外环动态方程
由dudc/dt=[iPV-P/udc]/C可得:
udc(t)=[Pref/udc(t-Δt)-P(t-Δt)/udc(t-Δt)]/CΔt+udc(t-Δt) (2)
其中,Δt为计算步长,t为当前时刻,X(t-Δt)(X可为P、udc)是上一步长下的计算结果。
②:计算id_ref
由式id_ref=kp(udc_ref-udc)即可得。
③:计算id
根据did/dt=(id_ref-id)/Td可得:
id(t)=id(t-Δt)+(id_ref-id(t-Δt))/TdΔt (3)
④:计算iq
由于iq的参考值一般设为0,即iq_ref=0,则根据diq/dt=(iq_ref-id)/Tq有:
iq(t)=iq(t-Δt)+(iq_ref-iq(t-Δt))/TqΔt (4)
图3展示了采用微分进化算法进行迭代求解,从而得到四个待辨识参数最优值的具体流程。该算法的核心思想在于通过计算id、iq,不断进行迭代优化使得两电流值与实测电网侧d,q轴电流值相差最小,则此时所对应的C,Kp,电Td,Tq即为辨识最优值。
图3中,设置三相短路故障即在仿真软件PSCAD中即可设置,如设置三相接地,并设置动态开关,指定故障持续时间;P的参考值为PSCAD运行所得稳定值,Q参考值已有说明,一般设为0,id_ref,iq_ref可通过电压和功率的初始值来进行计算,一般id-ref值为0。其中udc_ref可通过动态方程组中特性方程进行计算。所得的优化模型即为图2所示的简化辨识控制模型。
(3)在交流电网中针对电气元件建立所对应的频移(SFP)模型,从而使得在仿真时可以取得较大的仿真步长,且同时产生瞬时及宽频带相量波形,从而有效反应高频动态特性。
在一个典型的电力系统中,变量通常具有带通特性,其频带以基频ωs为中心,在频移模型中以相量形式表现为:
则对于交流电网中每一个电气元件而言,其动态方程在频移模型中为:
采用隐式梯形积分法对式(7)进行离散化处理可得:
将频移模型中的信号进行坐标变换转换到时域中,该信号形式可表示为:
在式(9)中,该动态响应信号的实部和虚部为:
(4)搭建SFP以及光伏辨识子系统间的接口模型,即所提出的多域传输线(MD-TLM)模型,该模型的建立取决于传输线时延τ与时间步长Δt之间的关系,模型拓扑结构如图4所示。
①:当τ>Δt时,MD-TLM可用双层诺顿等效电路来表示:
②:当τ<Δt时,在频域下,可通过线性插值的方式来获得电压以及电流相量,即:
对于MD-TLM模型有:
将式(13)代入式(12)则有:
(5)将频域以及时域中的变量进行相互变换从而获取多域传输线模型中的关键参数值,并将该多域传输线作为连接两子系统的接口模型,光伏辨识子系统与频移子系统的整体拓扑结构如图5所示,图6则展示了一相应的应用实例。
当从时域坐标系变换到频域下时:
其中,H代表Hilbert变换形式。
另外,当从频域坐标系变换到时域下时:
参见图7,一种基于在线辨识的大规模光伏并网的电磁暂态仿真方法,具体包括以下步骤。
(1)简化光伏逆变器控制系统模型:将光伏电池板简化为恒定电流源;电流内环PI控制简化为一阶惯性环节;电压外环的PI控制简化为P控制,仅含比例环节。
与简化模型相对应的状态空间方程组表示为:
为了提高光伏并网系统的功率因数,一般令Qref=0,即eq=0.。
由式(18)派克变换可知:此时iq_ref=0。
(2)在原复杂仿真模型中设置一对称三相接地故障以产生较大扰动,采集在该故障下原初始模型中光伏控制系统输入给逆变器的电流,电压以及功率等电信息,并根据该输入信息,基于优化算法以及辨识模型中的空间状态方程,确定辨识模型中的关键参数,即电容C,电压外环控制的比例环节参数Kp,电流内环控制的时间常数Tp,Tq,使得辨识模型与初始模型的外特性保持一致。
(3)将交流电网看作一子系统,对其中所含电气元件建立动态方程,并采用隐式梯形积分对其进行离散化处理,进而将信号或变量进行时域与频域坐标下的转换,从而建立对应的等效频移模型。
(4)在反映时域与频域之间的相互作用及关系时,建立多域传输线模型,即时域光伏辨识系统与频域交流电网系统间的接口模型,并根据传输线时延τ与时间步长Δt之间的大小关系,建立不同条件下的数学模型。在此基础上,进行时域与频域坐标系下变量的转换,进而获取多域传输线模型中的参数值搭建接口模型。从而连接两子系统进行立电磁暂态联合仿真,具体结果如图8所示。图中表明在本发明所提出的新型电磁暂态仿真方法下光伏并网系统的动态特性与原始系统的动态特性几乎保持一致,具有高度的可替代性。
本发明的优点及创新点如下。
相比于传统建模方法仿真技术,本发明采用辨识模型以及频移模型替代原复杂模型进行电磁暂态仿真的全新方法。其优点以及创新点在于:
(1)在该发明方法中,光伏电站的辨识模型保持着与实测数据动态响应一致的动态特性。
(2)该发明方法可以在优化算法的基础上,通过光伏电站模型的输入信息辨识出电站关键控制参数,并结合所给出的简化辨识模型,极大简化原复杂模型,降低了工程计算复杂程度。
(3)该方法中,交流电网采用频移模型来进行仿真,从而使得仿真步长可以增大,进而提高仿真效率以及精度,且同时产生瞬时及宽频带相量波形,从而有效反应高频动态特性。
(4)该方法中,辨识与频移两子模型间的接口特性通过多域传输线模型得以体现,可以有效反映两子模型间的宽频带相互作用。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统,其特征在于:所述仿真系统包括辨识模型和频移模型,所述辨识模型用于等效光伏电站,所述频移模型用于等效交流电网,所述辨识模型和频移模型之间采用多域传输线模型作为接口模型,并反映辨识模型和频移模型之间的宽频带相互作用关系;
所述多域传输线模型的建立取决于传输线时延τ与时间步长Δt之间的关系,
当τ>Δt时,多域传输线模型用双层诺顿等效电路来表示为:
当τ<Δt时,在频域下,多域传输线模型通过线性插值的方式来获得频域下电压以及电流相量,即:
其中p,q满足:
且MD-TLM模型表达式为:
将上述结果代入该模型表达式中则有:
2.根据权利要求1所述的电磁暂态仿真系统,其特征在于:所述辨识模型的关键参数包括电容C、电压外环控制的比例环节参数Kp、电流内环控制的时间常数Td,Tq。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010045691.2A CN111259548B (zh) | 2020-01-16 | 2020-01-16 | 一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010045691.2A CN111259548B (zh) | 2020-01-16 | 2020-01-16 | 一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111259548A CN111259548A (zh) | 2020-06-09 |
CN111259548B true CN111259548B (zh) | 2022-09-02 |
Family
ID=70948856
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010045691.2A Expired - Fee Related CN111259548B (zh) | 2020-01-16 | 2020-01-16 | 一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111259548B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102594215A (zh) * | 2012-03-16 | 2012-07-18 | 云南电力试验研究院(集团)有限公司电力研究院 | 一种光伏电站模型参数辨识方法 |
CN103208815A (zh) * | 2013-04-02 | 2013-07-17 | 清华大学 | 光伏发电系统并网逆变器的d-q轴参数辨识方法 |
CN104333326A (zh) * | 2014-11-19 | 2015-02-04 | 东南大学 | 一种光伏逆变器控制器外特性测试的试验方法 |
CN109657276A (zh) * | 2018-11-16 | 2019-04-19 | 清华大学 | 电磁暂态仿真方法及系统 |
CN110516314A (zh) * | 2019-08-02 | 2019-11-29 | 上海交通大学 | 用于大规模风电场的多速率多模态联合仿真方法 |
-
2020
- 2020-01-16 CN CN202010045691.2A patent/CN111259548B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102594215A (zh) * | 2012-03-16 | 2012-07-18 | 云南电力试验研究院(集团)有限公司电力研究院 | 一种光伏电站模型参数辨识方法 |
CN103208815A (zh) * | 2013-04-02 | 2013-07-17 | 清华大学 | 光伏发电系统并网逆变器的d-q轴参数辨识方法 |
CN104333326A (zh) * | 2014-11-19 | 2015-02-04 | 东南大学 | 一种光伏逆变器控制器外特性测试的试验方法 |
CN109657276A (zh) * | 2018-11-16 | 2019-04-19 | 清华大学 | 电磁暂态仿真方法及系统 |
CN110516314A (zh) * | 2019-08-02 | 2019-11-29 | 上海交通大学 | 用于大规模风电场的多速率多模态联合仿真方法 |
Non-Patent Citations (1)
Title |
---|
多频段–动态相量法电磁暂态仿真研究;刘博宁等;《中国电机工程学报》;20191005;第39卷(第19期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111259548A (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20190067939A1 (en) | Multi-time-scale digital/analog hybrid simulation system and method for power distribution network and storage medium | |
CN111541262B (zh) | 模型预测定交流电压控制下mmc频率耦合阻抗建模方法 | |
Chen et al. | Loop frame of reference based three-phase power flow for unbalanced radial distribution systems | |
CN109802406B (zh) | 一种分析柔性直流输电接入系统谐振稳定性的方法 | |
CN111654052B (zh) | 基于动态相量法的柔直换流器建模装置及方法 | |
CN113300383A (zh) | 一种机电暂态建模方法、系统、设备及存储介质 | |
CN110350515A (zh) | 一种适用于故障电流分析的柔性直流电网建模方法 | |
CN111146807A (zh) | 极坐标下的变流器并网系统小干扰稳定性判断方法 | |
CN105449674A (zh) | 航空12脉波逆变器的广义状态空间平均建模方法 | |
CN116937696B (zh) | 一种基于光伏发电系统的自适应等效建模方法 | |
CN111259548B (zh) | 一种基于在线辨识的大规模光伏并网的电磁暂态仿真系统 | |
CN111969639B (zh) | 电力电子化电网级联型换流器多时间尺度暂态建模方法 | |
CN108984847B (zh) | 一种基于分频阻抗补偿的实时数字混合仿真接口方法 | |
CN111628495B (zh) | 含vsc-mmc交直流电网系统的电磁暂态仿真方法 | |
CN111725818B (zh) | 一种三相弱电网变流器群并网仿真方法及仿真终端 | |
CN111987734A (zh) | 一种基于轨迹灵敏度的暂态过电压两阶段优化控制方法 | |
CN112467788A (zh) | 一种减少光伏模型预测控制系统低电压穿越时稳态误差的方法 | |
CN116706927B (zh) | 一种四桥臂逆变器最优电压支撑参考电流计算方法 | |
Wei et al. | Coordinate control of grid-side energy storage power station with consideration of communication delays in grid-connected micro-grids | |
Wei et al. | Segmented transmission delay based decoupling for parallel simulation of a distribution network | |
CN110190627B (zh) | 一种并网逆变器基于辅助环节的电网适应性控制算法 | |
CN115603332A (zh) | 一种柔性直流电网的稳定性分析方法及装置 | |
CN116470522B (zh) | 一种svg通用电磁暂态模型的控制参数识别方法及装置 | |
Shi et al. | An Aggregate Measurement based Identification Model for Large Scale Photovoltaic Stations and Applications on Real time Platform HYPERSIM | |
Bai et al. | A real-time co-simulation of PV power generation system using transmission line model interface |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220902 |