CN113435065B - 一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 - Google Patents
一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 Download PDFInfo
- Publication number
- CN113435065B CN113435065B CN202110908462.3A CN202110908462A CN113435065B CN 113435065 B CN113435065 B CN 113435065B CN 202110908462 A CN202110908462 A CN 202110908462A CN 113435065 B CN113435065 B CN 113435065B
- Authority
- CN
- China
- Prior art keywords
- transfer function
- converter
- modeling
- frequency band
- frequency
- 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 33
- 238000012546 transfer Methods 0.000 claims abstract description 64
- 238000005070 sampling Methods 0.000 claims abstract description 33
- 238000010586 diagram Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 12
- 230000007474 system interaction Effects 0.000 abstract description 3
- 230000003993 interaction Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation 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/20—Design optimisation, verification or simulation
-
- 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
- H02J1/00—Circuit arrangements for dc mains or dc distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M3/00—Conversion of dc power input into dc power output
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明提出一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,基于扫频结果选取主导特征频段,并在该频段内选取采样点;基于二阶传递函数进行“黑箱”DC/DC变换器一次建模;基于全频段一致性建立最终的等效模型。本发明的方法适用于实际系统变换器的有效建模,在变换器电路结构、零极点个数等关键电路信息未知的情况下,仅根据扫频结果,获得变换器的等效模型;等效模型用于系统稳定性分析时,与系统实际传递函数吻合度较高,适用于有建模或分析需求的含“黑箱”DC/DC变换器的直流配电网场景,对于掌握变换器接入对系统稳定性的影响,分析系统交互特性,预测可能出现的谐振点,改善电网电能质量等具有重要意义。
Description
技术领域
本发明属数据辨识技术领域,特别是涉及一种基于“黑箱”DC/DC变换器的等效阻抗建模方法。
背景技术
随着“双碳”目标的设定及柔性直流输电、电力电子变压器等技术的快速发展,直流配电网展现出其供电容量大、线路损耗小、更适于分布式电源接入等优势。分布式电源接入所带来的降低系统阻尼,与配电网中的恒功率负载发生交互作用,可能引起宽频振荡、谐振,降低系统的稳定性等问题。在实际系统中,由于商业竞争、用户信息保护等因素,部分变换器的电路参数、拓扑结构、控制方式等信息不能完全公开,导致部分DC/DC变换器呈现“黑箱”状态,不利于系统建模和稳定性分析。合理、准确地建立变换器的等效模型,能为提升系统稳定性、研究交互作用规律提供科学依据。
目前“黑箱”DC/DC变换器建模方法较少,现有建模方法有适用于所有直流功率变换装置的标准通用电路模型,通过一次全频带、融合扫频测量即可得到传递函数模型等。大部分的建模方法只分析装置本身的谐振特性,无法评估装置接入对系统稳定性的影响,无法建立全频段等值阻抗模型。
建模的目的在于分析配电网的交互作用和各模块之间的相互影响,为分布式能源接入后深入分析直流配电网变换器之间的稳定性和谐振分析、电能质量研究提供模型支撑。因此,适用、合理的“黑箱”DC/DC变换器的等效阻抗建模方法对分析直流配电网交互特性、提高配电网稳定性、改善直流配电网电能质量等具有重要意义。
鉴于此,本发明提出了一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,克服了现有配电网中大量电路拓扑或参数等信息未知的“黑箱”DC/DC变换器建模的问题,可用于实际系统的DC/DC变换器等效阻抗建模。
发明内容
为了解决上述问题,本发明提出了一种基于“黑箱”DC/DC变换器的等效阻抗建模方法。运用本发明中的等效阻抗建模方法,能够在电路拓扑、零极点个数等关键信息不完整的条件下获取DC/DC变换器等效传递函数及频域特征量。
为实现此目的,本发明采用如下技术方案。
一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,包括:
A、基于扫频结果选取主导特征频段,并基于该频段选取采样点;
B、基于二阶传递函数进行“黑箱”DC/DC变换器一次建模;
C、基于全频段一致性建立最终的等效模型。
进一步,所述步骤A包括:
步骤A1:根据扫频结果绘制Bode图,选取Bode图中存在波峰或幅值相位变化明显的频段,得到主导特征频段;
步骤A2:根据Bode图中的主导特征频段,采用穷举法,对所述主导特征频段内所有采样点进行组合,每个组合包括三个采样点。
进一步,所述步骤B包括:
步骤B1、基于二阶传递函数的采样点解的计算
根据二阶传递函数公式确定所选采样点对应的唯一解,实现对“黑箱”DC/DC变换器选取采样点的频段的一次建模,公式表示如下:
其中,A、B、C、D、E为传递函数分子和分母的系数,是求解的对象,s为传递函数的自变量;
步骤B2、基于最小误差的二阶传递函数参数的确定
计算G(s)频段误差,重复步骤A2和B1,选择出与扫频结果最为接近的传递函数,即对所有的采样点组合进行传递函数系数求解,得到在选择频段范围内与建模对象扫频结果误差最小的传递函数。
更进一步,所述步骤B1中,将二阶传递函数公式分离实部和虚部,形成五元一次方程组进行求解,从而得出每组采样点对应的唯一解,
Re[]为实部,Im[]为虚部,1,2,3为所选择的三个采样点。
进一步,所述步骤C包括:
步骤C1、基于一次建模的新建模对象的获取
将扫频结果除以所述步骤B1所得二阶传递函数,得到新的建模对象,公式表示如下:
式中,Q(s)为实际传递函数,在计算过程中用扫频结果代替;Gx(s)为第x次建模得到的二阶传递函数;
步骤C2、对于四阶以上系统,对比所得传递函数与扫频结果Bode图是否全频段一致,多次重复步骤A、B、C1,直至在整个频域内,所得传递函数Bode图与扫频结果基本一致,满足精度需求,得到最终模型,公式表示如下:
附图说明
图1是本发明具体实施方式中基于“黑箱”DC/DC变换器的等效阻抗建模方法的流程图;
图2是本发明具体实施方式中一次建模的采样点位置图;
图3是本发明具体实施方式中一次建模Bode图;
图4是本发明具体实施方式中二次建模Bode图;
图5是本发明具体实施方式中三次建模Bode图;
图6是本发明具体实施方式中与实际传递函数对比Bode图;
图7是本发明具体实施方式中与实际传递函数对比Bode图。
具体实施方式
下面结合附图表,对本发明作详细说明。
以下公开详细的推理分析方法和示范分析例。然而,此处公开的具体推理及分析过程细节仅仅是出于描述示范分析例的目的。
然而,应该理解,本发明不局限于公开的具体示范实施例,而是覆盖落入本公开范围内的所有修改、等同物和替换物。在对全部附图的描述中,相同的附图标记表示相同的元件。
通过本发明的一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,能够简洁、有效地表述“黑箱”DC/DC变换器的频域特性,得出频域特征量,可用于系统稳定性及交互特性的研究。
如图1所示,本发明实施方式中“黑箱”DC/DC变换器的等效阻抗建模方法包括以下步骤:
步骤A、基于扫频结果选取主导特征频段,并基于该频段选取采样点。
步骤A1:根据扫频结果绘制Bode图,观察Bode图,选取图中存在波峰或幅值相位变化明显的频段,确定该频段为主导特征频段。
步骤A2:根据Bode图中的主导特征频段,在该频段内选择采样点。采样点的选择是影响建模准确性的重要因素,采用穷举法,对特征频段内所有采样点进行组合,每个组合包括三个采样点。
步骤B、基于二阶传递函数进行“黑箱”DC/DC变换器一次建模;
步骤B1、基于二阶传递函数的采样点解的计算
在控制理论中,一定条件下,高阶系统可以降维为二阶系统,而不失其主导运动特征。因此,本发明以二阶传递函数的形式对“黑箱”DC/DC变换器选取采样点的频段进行建模。将二阶传递函数公式分离实部和虚部,形成五元一次方程组进行求解,从而得出每组采样点对应的唯一解,
实现对“黑箱”DC/DC变换器选取采样点的频段的一次建模。公式表示如下:
其中,A、B、C、D、E为传递函数分子和分母的系数,是求解的对象,s为传递函数的自变量。
对于三个采样点的一个组合,其中:
Re[]为实部,Im[]为虚部,1,2,3为所选择的三个采样点。
步骤B2、基于最小误差的二阶传递函数参数的确定
由于受高阶部分干扰,经上述步骤等效的传递函数与实际传递函数仍存在一定的差别。因此,计算G(s)频段误差,重复步骤A2、B1,选择与扫频结果最为接近的传递函数,即对所有的采样点组合进行传递函数系数求解,得到在选择频段范围内与建模对象扫频结果误差最小的传递函数。
步骤C、基于全频段一致性建立最终的等效模型。
步骤C1、基于一次建模的新建模对象的获取
对于二阶以上的传递函数,二阶模型只能对其中部分频段进行较为准确地建模,较难实现完整Bode图的等效。因此,为了建立全频段等值阻抗模型,在所建立二阶模型的基础上,需将扫频结果除以上述步骤得到的二阶传递函数,得到新的建模对象。每次新建模对象应为扫频结果除以之前得到的所有二阶模型之积,公式表示如下:
式中,Q(s)为实际传递函数,在计算过程中用扫频结果代替;Gx(s)为第x次建模得到的二阶传递函数。
步骤C2、对于四阶以上系统,需多次重复步骤A、B、C1,对比所得传递函数与扫频结果Bode图是否全频段一致,直至全轨迹基本重合,进而得到与扫频结果全频段一致的最终模型。公式表示如下:
以下通过一个具体实施方式来说明本发明。
提取某直流配电网中分布式电源的光伏输出端等效阻抗Zpv扫频结果,观察其Bode图确定Zpv的主要特征频段为20rad/s-200rad/s,并在该频段内选择采样点。根据二次传递函数公式,确定A,B,C,D,E的值。通过对比,所建传递函数在主要特征频段内误差最小时,采样点位置如图2所示,求得传递函数如下:
如图3G1(s)的Bode图所示,G1(s)仅在低频段及所选特征频段与Zpv轨迹保持一致,无法体现其他频段变化规律。选取二次建模对象Zpv/G1(s),主要特征频段为700rad/s-5000rad/s,重复以上步骤,得到误差最小时的传递函数G2(s),实际传递函数为G(s)=G1(s)*G2(s),观察图4的Bode图可知,G(s)与Zpv轨迹仍存在差异;因此,需进行三次建模,步骤同上,求得G3(s),实际传递函数为G(s)=G1(s)*G2(s)*G3(s),如图5所示,此时G(s)与Zpv在全频段变化基本一致。此外,列出了G(s)和Zpv的频域特征量,如表1所示。
表1 Zpv等效模型特征量分析数据
利用本发明所提“黑箱”建模方法形成等效仿真模型,得和/>与实际传递函数对比图,如图6、7,对比发现,二者吻合度较高,本发明能够较好地体现系统的频域特性。
本发明所提出的基于“黑箱”DC/DC变换器的等效阻抗建模方法,主要优点在于在变换器电路结构及参数等关键电路信息未知的情况下,仅根据扫频结果,得出等效传递函数及频域特征量;所建立等效模型与实际输出结果一致性较好,能够较好地体现系统的输出特性,为系统稳定性、谐振等问题的分析提供了模型基础。本发明无需提供除扫频结果以外的任何信息,有效避免了其他建模方法需要提供零极点个数、拓扑结构等信息的限制,有效解决了在上述信息未知情况下仅能根据扫频结果,分析装置本身的谐振特性,数据量大且无法评估装置接入对系统交互特性的影响等问题,适用于有建模或分析需求的含“黑箱”DC/DC变换器的直流配电网场景。
通过上述说明,本发明基于“黑箱”DC/DC变换器的等效阻抗建模方法的基本功能得到了阐述。本发明的一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,实现了对DC/DC变换器的有效建模,克服了现有建模方法中的不足,对于掌握变换器接入对系统稳定性的影响,分析系统交互特性,预测可能出现的谐振点,改善电网电能质量等具有重要意义。
需要说明的是,上述实施方式仅为本发明较佳的实施方式,不能将其理解为对本发明保护范围的限制,在未脱离本发明构思的前提下,对本发明所做的任何微小变化与修饰均属于本发明的保护范围。
Claims (2)
1.一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,其特征在于,包括以下步骤:
步骤A、基于扫频结果选取主导特征频段,并基于该频段选取采样点,包括:
步骤A1:根据扫频结果绘制Bode图,选取Bode图中存在波峰或幅值相位变化明显的频段,得到主导特征频段;
步骤A2:根据Bode图中的主导特征频段,采用穷举法,对所述主导特征频段内所有采样点进行组合,每个组合包括三个采样点;
步骤B、基于二阶传递函数进行“黑箱”DC/DC变换器一次建模,包括:
步骤B1、基于二阶传递函数的采样点解的计算
根据二阶传递函数公式确定所选采样点对应的唯一解,实现对“黑箱”
DC/DC变换器选取采样点的频段的一次建模,公式表示如下:
其中,A、B、C、D、E为传递函数分子和分母的系数,是求解的对象,s为传递函数的自变量;
步骤B2、基于最小误差的二阶传递函数参数的确定
计算G(s)频段误差,重复步骤A2和B1,选择出与扫频结果最为接近的传递函数,即对所有的采样点组合进行传递函数系数求解,得到在选择频段范围内与建模对象扫频结果误差最小的传递函数;
步骤C、基于全频段一致性建立最终的等效模型,包括:
步骤C1、基于一次建模的新建模对象的获取
将扫频结果除以所述步骤B1所得二阶传递函数,得到新的建模对象,公式表示如下:
式中,Q(s)为实际传递函数,在计算过程中用扫频结果代替;Gx(s)为第x次建模得到的二阶传递函数;
步骤C2、对于四阶以上系统,对比所得传递函数与扫频结果Bode图是否全频段一致,多次重复步骤A、B、C1,直至在整个频域内,所得传递函数Bode图与扫频结果基本一致,满足精度需求,得到最终模型,公式表示如下:
2.根据权利要求1所述的一种基于“黑箱”DC/DC变换器的等效阻抗建模方法,其特征在于:所述步骤B1中,将二阶传递函数公式分离实部和虚部,形成五元一次方程组进行求解,从而得出每组采样点对应的唯一解,
Re[]为实部,Im[]为虚部,1,2,3为所选择的三个采样点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110908462.3A CN113435065B (zh) | 2021-08-09 | 2021-08-09 | 一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110908462.3A CN113435065B (zh) | 2021-08-09 | 2021-08-09 | 一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113435065A CN113435065A (zh) | 2021-09-24 |
CN113435065B true CN113435065B (zh) | 2024-03-15 |
Family
ID=77763047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110908462.3A Active CN113435065B (zh) | 2021-08-09 | 2021-08-09 | 一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113435065B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114880863B (zh) * | 2022-05-20 | 2024-03-01 | 合肥工业大学 | 分布式可再生能源集群阻抗聚合模型自适应分频降阶方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107317321A (zh) * | 2017-05-15 | 2017-11-03 | 北京交通大学 | Dc‑dc变换器并联系统的稳定性分析方法 |
CN109980634A (zh) * | 2019-02-28 | 2019-07-05 | 南京工程学院 | 一种基于谐波传递函数矩阵的dc/dc变换器稳定性分析方法 |
CN110471310A (zh) * | 2019-08-13 | 2019-11-19 | 广西电网有限责任公司电力科学研究院 | 一种功率型数模混合仿真系统稳定性改进方法 |
CN110535348A (zh) * | 2019-09-03 | 2019-12-03 | 华北电力大学 | 直流换流器高频阻抗建模方法 |
CN113054640A (zh) * | 2021-03-03 | 2021-06-29 | 湖南大学 | 一种基于阻抗分解的直流变换器并联系统稳定性判据方法 |
-
2021
- 2021-08-09 CN CN202110908462.3A patent/CN113435065B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107317321A (zh) * | 2017-05-15 | 2017-11-03 | 北京交通大学 | Dc‑dc变换器并联系统的稳定性分析方法 |
CN109980634A (zh) * | 2019-02-28 | 2019-07-05 | 南京工程学院 | 一种基于谐波传递函数矩阵的dc/dc变换器稳定性分析方法 |
CN110471310A (zh) * | 2019-08-13 | 2019-11-19 | 广西电网有限责任公司电力科学研究院 | 一种功率型数模混合仿真系统稳定性改进方法 |
CN110535348A (zh) * | 2019-09-03 | 2019-12-03 | 华北电力大学 | 直流换流器高频阻抗建模方法 |
CN113054640A (zh) * | 2021-03-03 | 2021-06-29 | 湖南大学 | 一种基于阻抗分解的直流变换器并联系统稳定性判据方法 |
Non-Patent Citations (2)
Title |
---|
Black-Box Modelling of Low-Switching-Frequency Power Inverters for EMC Analyses in Renewable Power Systems;Luis Arnedo 等;《2009 Twenty-Fourth Annual IEEE Applied Power Electronics Conference and Exposition》;1476-1481 * |
基于频域法的DC/DC变换器传递函数设计方法研究;高嵩 等;《系统仿真学报》;第34卷(第2期);334-341 * |
Also Published As
Publication number | Publication date |
---|---|
CN113435065A (zh) | 2021-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107453357B (zh) | 一种基于分层求解的配电网状态估计方法 | |
CN108667048B (zh) | 新能源并网系统振荡稳定性的频域判稳方法及装置 | |
KR101219545B1 (ko) | 전력계통에서의 최적화 기법을 적용한 파라미터 추정 방법 | |
CN109802406B (zh) | 一种分析柔性直流输电接入系统谐振稳定性的方法 | |
CN102305891A (zh) | 一种电力系统低频振荡在线监测方法 | |
CN113078630B (zh) | 一种基于实时量测数据的低压配电网拓扑辨识方法 | |
CN114252103B (zh) | 融合电站运行故障预测方法 | |
CN105938578A (zh) | 一种基于聚类分析的大规模光伏电站等值建模方法 | |
CN113435065B (zh) | 一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 | |
CN112966219A (zh) | 识别户表与表箱关系的方法及系统、设备、介质 | |
CN115453193B (zh) | 基于pqm、ttu和sm量测数据协同的配电网谐波状态估计方法 | |
CN105160058B (zh) | 一种用于电网电磁暂态计算的网络边界等值方法 | |
CN113114497B (zh) | 一种低压台区拓扑自动识别方法 | |
CN102856896A (zh) | 一种直流输电损耗的在线分析方法 | |
CN115036924A (zh) | 一种考虑多源数据不确定性的配电网分布式状态估计方法 | |
CN115877312A (zh) | 一种基于台区电能量守恒的电能表信息化评价校准模型 | |
CN110707693A (zh) | 一种基于ami全量测点分区的集合卡尔曼滤波动态状态估计方法 | |
CN113466624A (zh) | 一种多端混合直流输电线路故障区域的检测方法及系统 | |
CN113985339A (zh) | 智能电表的误差诊断方法及系统、设备、存储介质 | |
CN115912349A (zh) | 一种基于压缩感知和图卷积神经网络的宽频振荡扰动源定位方法 | |
CN112487748B (zh) | 一种考虑时频域特性的锂离子电池分数阶模型建立方法 | |
CN114421465A (zh) | 基于元件拓扑的电网参数辨识校验方法、装置及存储介质 | |
CN110865328B (zh) | 基于ami的智能电表相位识别、拓扑辨识及阻抗估计方法 | |
CN115347553B (zh) | 基于线性回归初值估计的配电网线路参数辨识方法 | |
CN106597848B (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 |