CN110826018A - 基于有理函数逼近的一次信号恢复方法及系统 - Google Patents
基于有理函数逼近的一次信号恢复方法及系统 Download PDFInfo
- Publication number
- CN110826018A CN110826018A CN202010022734.5A CN202010022734A CN110826018A CN 110826018 A CN110826018 A CN 110826018A CN 202010022734 A CN202010022734 A CN 202010022734A CN 110826018 A CN110826018 A CN 110826018A
- Authority
- CN
- China
- Prior art keywords
- transfer function
- function
- fitting
- current voltage
- voltage transformer
- 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
Links
Images
Classifications
-
- 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/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
Abstract
本发明公开了基于有理函数逼近的一次信号恢复方法及系统,包括:对直流电压互感器一次侧输入不同频率的电压;在直流电压互感器二次侧测量同频电压;求二次侧测量值与一次侧测量值之比,获得不同频率下的互感器传递函数值;根据厂商给定的等效电路图确定直流电压互感器拟合传递函数的阶数;将拟合传递函数与测量值利用有理函数逼近原理构建误差函数;利用最优化算法求解误差函数的最优解和拟合传递函数的系数;利用二次侧电压信号和拟合传递函数恢复一次侧电压信号。本发明通过测量传递函数值拟合得到的互感器传递函数与二次侧电压信号能准确恢复一次侧电压信号,仿真分析验证了所提方法在宽频域下具有较高的准确性和良好的抗噪能力。
Description
技术领域
本发明涉及直流电压互感器研究领域,具体地,涉及基于有理函数逼近的一次信号恢复方法及系统。
背景技术
特高压直流输电是目前常用的大容量远距离输电技术,随着国内特高压直流输电网络的建设,直流输电引起的谐波不稳定以及非线性负载的增加,电网中的谐波问题变得越发严重,直流电压互感器(DC voltage transformer,DCTV)能有效测量一次侧的直流电压,但难以测量系统中的谐波电压。另一方面,系统发生故障时DCTV二次侧电压难以准确反映系统的暂态情况,可能导致保护装置延时动作或误动作,严重时影响电力系统的安全运行。而对DCTV宽频域传递特性的准确求取将有效提高DCTV测量一次侧电压的准确性。
目前对于特高压直流互感器的宽频域传递特性的建模方法研究相对滞后,主要集中在互感器等效电路仿真建模及理论研究。其一,利用测量得到的散射参数和矢量匹配法对电压传输函数进行有理函数逼近得到互感器传递函数,为了降低模型复杂度,这种方法只能测量一定频率范围内的数据,降低了模型精度。其二,利用谐波源法得到互感器的传递函数,其传递特性不具有泛化能力,仅能反映该谐波频率点的测量特征。最常用的方法是对阻容式直流电压互感器进行集中参数建模和分布参数建模,由于二次分压系统中隔离放大器截断频率难以获取,使二次分压系统的传递函数难以确定,只能根据厂家给定的元件额定参数推算直流电压互感器传递函数,导致产生较大的误差。其主要原因有三点:第一,分压电阻的自热效应和环境温度会导致额定电阻阻值变化;第二,安装电阻的绝缘支架有泄露电流会使高压臂的实际阻值发生变化;第三,在潮湿的环境中,湿气侵入电阻内部,通电时会发生电腐蚀,导致电阻值变化。
综上所述,目前缺乏一种不依赖直流电压互感器内部参数可以准确求取互感器的宽频域传递特性的方法,导致二次侧电压信号不能准确恢复一次侧电压信号。
发明内容
传统方法利用额定参数推算直流电压互感器传递函数,但由于互感器长时间工作后元件腐蚀老化,参数改变,导致推算得到的传递函数不能准确反映直流电压互感器实际的传递特性。针对这一问题,本发明通过将DCTV等效为一个二端口系统,通过扫频法获取不同频率下直流电压互感器的传递函数值,同时根据等效电路图确定拟合传递函数阶数,利用有理函数逼近原理构造拟合传递函数与测量值的误差函数,通过采用粒子群优化算法求取该误差函数最小值来确定拟合传递函数的系数,将所求系数代入拟合传递函数,便可准确得到DCTV的传递函数。利用二次侧电压信号和DCTV传递函数可以准确地恢复一次侧电压信号。
本发明通过下述技术方案实现:
基于有理函数逼近的一次信号恢复方法,包括以下步骤:
A、将直流电压互感器视为一个二端口系统,通过扫频法求得直流电压互感器在不同频率下的传递函数值,完成采样过程;
B、根据直流电压互感器的等效电路图确定直流电压互感器传递函数阶数并构建拟合传递函数,将采样值与拟合传递函数利用有理函数逼近原理构建误差函数,利用优化算法求误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的逼近拟合;
C、根据逼近拟合得到的DCTV传递函数和二次侧电压信号实现一次侧电压信号的恢复。
进一步地,基于有理函数逼近的一次信号恢复方法,所述步骤A具体包括:
A1、在直流电压互感器一次侧施加不同频率的电压信号;
A2、在直流电压互感器二次侧测量同频电压信号,二次侧电压值和一次侧电压值之比便得到该频率点的传递函数值,实现传递函数值采样。
进一步地,基于有理函数逼近的一次信号恢复方法,所述步骤B具体包括:
B2、综合考虑幅频误差和相频误差,利用有理函数逼近原理构造拟合传递函数和采样值的误差函数为:
(3)
B3、误差函数J展开后为非线性函数,对非线性函数最优化求解有多种算法,如遗传算法(GA)、模退火算法(SA)、粒子群算法(PSO)和梯度下降法(GD)等,利用合适的优化算法求解误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的拟合。
进一步的,基于有理函数逼近的一次信号恢复方法,所述步骤C具体包括:
本发明还提供了一种基于有理函数逼近的一次信号恢复系统,所述系统包括:
采样单元,用于通过扫频法求得直流电压互感器在不同频率下的传递函数值,获得采样值;
处理单元,用于获得直流电压互感器传递函数阶数,并构建拟合传递函数,基于采样值与拟合传递函数构建误差函数,求解误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的逼近拟合;
恢复单元,用于根据逼近拟合得到的直流电压互感器传递函数和二次侧电压信号完成一次侧电压信号的恢复。
本发明提供的一个或多个技术方案,至少具有如下技术效果或优点:
现有方法通常误差较大或者得到的传递特性不具有泛化能力,仅能反映该谐波频率点的测量特征。本发明通过测量DCTV不同频率的传递函数值拟合得到DCTV的传递函数,只依赖于实测数据自身的属性,反映的是系统自身特性,得到的传递函数能准确地反映DCTV的宽频域传递特性,利用拟合传递函数和二次侧电压信号能够更加准确地恢复一次侧电压信号。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定;
图1为本发明流程示意图;
图2为本发明直流电压互感器等效电路模型图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在相互不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述范围内的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
本领域技术人员应理解的是,在本发明的揭露中,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系是基于附图所示的方位或位置关系,其仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此上述术语不能理解为对本发明的限制。
可以理解的是,术语“一”应理解为“至少一”或“一个或多个”,即在一个实施例中,一个元件的数量可以为一个,而在另外的实施例中,该元件的数量可以为多个,术语“一”不能理解为对数量的限制。
实施例一
请参考图1,本发明实施例提供了基于有理函数逼近的一次信号恢复方法,包括以下步骤:
A、将直流电压互感器视为一个二端口系统,通过扫频法求得直流电压互感器在不同频率下的传递函数值,完成采样过程;
B、根据直流电压互感器的等效电路图确定直流电压互感器传递函数阶数并构建拟合传递函数,将采样值与拟合传递函数利用有理函数逼近原理构建误差函数,利用优化算法求误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的逼近拟合;
C、根据逼近拟合得到的DCTV传递函数和二次侧电压信号实现一次侧电压信号的恢复。
其中,在本发明实施例中,基于有理函数逼近的一次信号恢复方法所述步骤A具体包括:
A1、在直流电压互感器一次侧施加不同频率的电压信号;
A2、在直流电压互感器二次侧测量同频电压信号,二次侧电压值和一次侧电压值之比便得到该频率点的传递函数值,实现传递函数值采样。
其中,在本发明实施例中,基于有理函数逼近的一次信号恢复方法,所述步骤B具体包括:
B2、综合考虑幅频误差和相频误差,利用有理函数逼近原理构造拟合传递函数和采样值的误差函数为:
B3、误差函数J展开后为非线性函数,对非线性函数最优化求解有多种算法,如遗传算法(GA)、模退火算法(SA)、粒子群算法(PSO)和梯度下降法(GD)等,利用合适的优化算法求解误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的拟合。
其中,在本发明实施例中,基于有理函数逼近的一次信号恢复方法,所述步骤C具体包括:
根据图2等效电路图,可以确定拟合传递函数分子分母阶数都为2,构建拟合传递函数:
通过仿真计算,得到的拟合传递函数与实际传递函数在角频率ω a=0到ω b=63000rad/s的每个频率对应的平均相对误差为2.1908e-16%。
恢复一次侧电压信号,具体方法为:
当DCTV处于谐波状态下时,分别就不同频率谐波输入,进行一次侧电压信号的还原。
表1
表2
表2为本发明实施例中元件参数。
其中,在本发明实施例中,本方法在实际应用中,本方法包括:
步骤一,对直流电压互感器一次侧输入不同频率的电压;步骤二,在直流电压互感器二次侧测量同频电压;步骤三,求二次侧测量值与一次侧测量值之比,获得不同频率下的互感器传递函数值;步骤四,根据厂商给定的等效电路图确定直流电压互感器拟合传递函数的阶数;步骤五,将拟合传递函数与测得的测量值利用有理函数逼近原理构建误差函数;步骤六,利用优化算法求解误差函数的最优解和拟合传递函数的系数;步骤七,利用二次侧电压信号和拟合传递函数恢复一次侧电压信号。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.基于有理函数逼近的一次信号恢复方法,其特征在于,所述方法包括:
A、通过扫频法求得直流电压互感器在不同频率下的传递函数值,基于传递函数值获得采样值;
B、获得直流电压互感器传递函数阶数,并构建拟合传递函数,基于采样值与拟合传递函数构建误差函数,求解误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的逼近拟合;
C、根据逼近拟合得到的直流电压互感器传递函数和二次侧电压信号完成一次侧电压信号的恢复。
2.根据权利要求1所述的基于有理函数逼近的一次信号恢复方法,其特征在于,将采样值与拟合传递函数利用有理函数逼近原理构建误差函数。
3.根据权利要求1所述的基于有理函数逼近的一次信号恢复方法,其特征在于,所述步骤A具体包括:
A1、在直流电压互感器一次侧施加不同频率的电压信号;
A2、在直流电压互感器二次侧测量同频电压信号,基于二次侧电压值和一次侧电压值之比,获得该频率点的传递函数值,实现传递函数值采样。
4.根据权利要求1-3中任意一个所述的基于有理函数逼近的一次信号恢复方法,其特征在于,所述步骤B具体包括:
B2、综合考虑幅频误差和相频误差,利用有理函数逼近原理构造拟合传递函数和采样值的误差函数为:
B3、求解误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的拟合。
5.根据权利要求1-3中任意一个所述的基于有理函数逼近的一次信号恢复方法,其特征在于,利用优化方法求解误差函数的最优解及拟合传递函数的系数,优化方法包括:遗传算法或模拟退火算法或粒子群算法或蚁群算法或免疫克隆选择算法或禁忌搜索算法或人群搜索算法或细菌觅食算法或填充函数法或梯度下降法。
7.一种基于有理函数逼近的一次信号恢复系统,其特征在于,所述系统包括:
采样单元,用于通过扫频法求得直流电压互感器在不同频率下的传递函数值,基于传递函数值获得采样值;
处理单元,用于获得直流电压互感器传递函数阶数,并构建拟合传递函数,基于采样值与拟合传递函数构建误差函数,求解误差函数的最优解及拟合传递函数的系数,完成对直流电压互感器传递函数的逼近拟合;
恢复单元,用于根据逼近拟合得到的直流电压互感器传递函数和二次侧电压信号完成一次侧电压信号的恢复。
8.根据权利要求7所述的基于有理函数逼近的一次信号恢复系统,其特征在于,所述采样单元具体用于:
在直流电压互感器一次侧施加不同频率的电压信号;
在直流电压互感器二次侧测量同频电压信号,基于二次侧电压值和一次侧电压值之比,获得该频率点的传递函数值,实现传递函数值采样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010022734.5A CN110826018B (zh) | 2020-01-09 | 2020-01-09 | 基于有理函数逼近的一次信号恢复方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010022734.5A CN110826018B (zh) | 2020-01-09 | 2020-01-09 | 基于有理函数逼近的一次信号恢复方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110826018A true CN110826018A (zh) | 2020-02-21 |
CN110826018B CN110826018B (zh) | 2020-04-14 |
Family
ID=69546514
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010022734.5A Expired - Fee Related CN110826018B (zh) | 2020-01-09 | 2020-01-09 | 基于有理函数逼近的一次信号恢复方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110826018B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111610482A (zh) * | 2020-06-10 | 2020-09-01 | 国网四川省电力公司电力科学研究院 | 一种换流站直流电压互感器宽频传递特性测试系统及方法 |
CN112861469A (zh) * | 2021-03-09 | 2021-05-28 | 杭州电力设备制造有限公司 | 配电变压器模型参数估计方法、装置、设备和存储介质 |
CN115078820A (zh) * | 2022-08-19 | 2022-09-20 | 石家庄科林电气股份有限公司 | 低压智能断路器的保护电流互感器饱和处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103840674A (zh) * | 2014-03-04 | 2014-06-04 | 国家电网公司 | 一种校准直流互感器的自适应电流发生系统及其发生方法 |
US20140268926A1 (en) * | 2013-03-14 | 2014-09-18 | General Electric Company | High voltage direct current (hvdc) converter system and method of operating the same |
CN106772189A (zh) * | 2016-11-18 | 2017-05-31 | 云南电网有限责任公司电力科学研究院 | 一种直流电压互感器同步校准方法及系统 |
CN107017654A (zh) * | 2017-05-31 | 2017-08-04 | 国网四川省电力公司电力科学研究院 | 一种抑制交流暂时过电压的特高压直流输电控制装置 |
-
2020
- 2020-01-09 CN CN202010022734.5A patent/CN110826018B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140268926A1 (en) * | 2013-03-14 | 2014-09-18 | General Electric Company | High voltage direct current (hvdc) converter system and method of operating the same |
CN103840674A (zh) * | 2014-03-04 | 2014-06-04 | 国家电网公司 | 一种校准直流互感器的自适应电流发生系统及其发生方法 |
CN106772189A (zh) * | 2016-11-18 | 2017-05-31 | 云南电网有限责任公司电力科学研究院 | 一种直流电压互感器同步校准方法及系统 |
CN107017654A (zh) * | 2017-05-31 | 2017-08-04 | 国网四川省电力公司电力科学研究院 | 一种抑制交流暂时过电压的特高压直流输电控制装置 |
Non-Patent Citations (5)
Title |
---|
QIN SHU 等: ""Estimate Utility Harmonic Impedance via the Cor-relation of Harmonic Measurements in Different Time Intervals"", 《IEEE》 * |
舒勤 等: ""多直流馈入受端城市电网直流落点处系统侧谐波阻抗估计方法"", 《高电压技术》 * |
董丽梅 等: ""基于遗传算法的EMD电力信号去噪"", 《计算机仿真》 * |
郝巍 等: ""直流输电引起的谐波不稳定及其相关问题"", 《电力系统自动化》 * |
黄银龙 等: ""一种新型的高压直流输电系统直流侧谐波电压测量方法"", 《电工技术学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111610482A (zh) * | 2020-06-10 | 2020-09-01 | 国网四川省电力公司电力科学研究院 | 一种换流站直流电压互感器宽频传递特性测试系统及方法 |
CN111610482B (zh) * | 2020-06-10 | 2022-09-13 | 国网四川省电力公司电力科学研究院 | 一种换流站直流电压互感器宽频传递特性测试系统及方法 |
CN112861469A (zh) * | 2021-03-09 | 2021-05-28 | 杭州电力设备制造有限公司 | 配电变压器模型参数估计方法、装置、设备和存储介质 |
CN115078820A (zh) * | 2022-08-19 | 2022-09-20 | 石家庄科林电气股份有限公司 | 低压智能断路器的保护电流互感器饱和处理方法 |
CN115078820B (zh) * | 2022-08-19 | 2022-11-18 | 石家庄科林电气股份有限公司 | 低压智能断路器的保护电流互感器饱和处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110826018B (zh) | 2020-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110826018B (zh) | 基于有理函数逼近的一次信号恢复方法及系统 | |
US20170030958A1 (en) | Transformer parameter estimation using terminal measurements | |
Apostolopoulos et al. | Accurate fault location algorithm for double-circuit series compensated lines using a limited number of two-end synchronized measurements | |
CN103401238B (zh) | 一种基于总体测辨法的电力负荷建模方法 | |
WO2012034528A1 (zh) | 一种计算电网一次时间常数的方法 | |
Radojevic et al. | Smart overhead lines autoreclosure algorithm based on detailed fault analysis | |
Guo et al. | Global-sensitivity-based theoretical analysis and fast prediction of traveling waves with respect to fault resistance on HVDC transmission lines | |
CN111141995B (zh) | 一种基于比幅原理的线路双端稳态量测距方法和系统 | |
Prabakar et al. | Use of traveling wave signatures in medium-voltage distribution systems for fault detection and location | |
CN112415273A (zh) | 一种双回非全线平行输电线路零序参数精确测量方法 | |
CN103487778B (zh) | 一种模拟量校正方法及基于该方法的单cpu低压保护装置 | |
Saber | A new fault location technique for four‐circuit series‐compensated transmission lines | |
Karimifard et al. | Determination of axial displacement extent based on transformer winding transfer function estimation using vector‐fitting method | |
CN110024249B (zh) | 用于检测电力传输线路中的故障的方法和使用该方法的保护系统 | |
Goklani et al. | Quantification of minimum unbalance required for accurate estimation of sequence parameters of transmission line using PMU data | |
JP2020076671A (ja) | 対地静電容量測定装置および対地静電容量測定方法 | |
Lin et al. | A hybrid method for current transformer saturation detection and compensation in smart grid | |
CN109713641B (zh) | 一种发电机定子接地保护系统及参数测量方法 | |
Gani et al. | A closed-form mathematical model and method for fast fault location on a Low voltage dc feeder using single-ended measurements | |
Filomena et al. | Extended impedance-based fault location formulation for unbalanced underground distribution systems | |
Banks et al. | Parameter estimation of a distribution transformer model using pseudo-random impulse sequence perturbation | |
CN110726887A (zh) | 一种就地化电力系统电压暂降检测方法 | |
Kamath et al. | Modeling of transformer characteristics using fractional order transfer functions | |
Milojevic et al. | A robust synchrophasor-based method for determining three-phase line impedance | |
Kang et al. | Parameter identification algorithm for fault location based on distributed transmission line model |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210122 Address after: 610000, No. 24, south section of Ring Road, Sichuan, Chengdu Patentee after: SICHUAN University Patentee after: STATE GRID SICHUAN ELECTRIC POWER Research Institute Address before: 610000, No. 24, south section of Ring Road, Sichuan, Chengdu Patentee before: SICHUAN University |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200414 Termination date: 20220109 |
|
CF01 | Termination of patent right due to non-payment of annual fee |