CN114172401B - 基于降阶观测器的npc三电平逆变器多类故障诊断方法 - Google Patents

基于降阶观测器的npc三电平逆变器多类故障诊断方法 Download PDF

Info

Publication number
CN114172401B
CN114172401B CN202111344797.3A CN202111344797A CN114172401B CN 114172401 B CN114172401 B CN 114172401B CN 202111344797 A CN202111344797 A CN 202111344797A CN 114172401 B CN114172401 B CN 114172401B
Authority
CN
China
Prior art keywords
phase
matrix
fault
inverter
sensor
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
CN202111344797.3A
Other languages
English (en)
Other versions
CN114172401A (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202111344797.3A priority Critical patent/CN114172401B/zh
Publication of CN114172401A publication Critical patent/CN114172401A/zh
Application granted granted Critical
Publication of CN114172401B publication Critical patent/CN114172401B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS 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
    • H02M7/00Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/42Conversion of dc power input into ac power output without possibility of reversal
    • H02M7/44Conversion of dc power input into ac power output without possibility of reversal by static converters
    • H02M7/48Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/53Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
    • H02M7/537Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
    • H02M7/5387Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration
    • H02M7/53871Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration with automatic control of output voltage or current
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/40Testing power supplies
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS 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
    • H02M1/00Details of apparatus for conversion
    • H02M1/08Circuits specially adapted for the generation of control voltages for semiconductor devices incorporated in static converters
    • H02M1/088Circuits specially adapted for the generation of control voltages for semiconductor devices incorporated in static converters for the simultaneous control of series or parallel connected semiconductor devices
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Inverter Devices (AREA)

Abstract

本发明公开了基于降阶观测器的NPC三电平逆变器多类故障诊断方法,属于故障诊断技术领域。该方法包括以下步骤:建立混合逻辑动态模型,利用矩阵变换进行解耦,将传感器故障与系统状态分解并建立降阶观测器,定义电流形态因子、故障诊断自适应阈值,进行故障诊断,再进行故障定位。本发明采用矩阵变换将传感器故障与系统状态进行分解,使得可以同时诊断逆变器开关管故障和传感器故障;本发明采用的电流形态因子和自适应阈值提高故障诊断的准确性和鲁棒性;本发明通过三相输出电流和传感器故障估计值的特征,来定义开关管故障检测特征值、传感器故障检测特征值和故障定位特征值对NPC三电平逆变器进行多类故障诊断。

Description

基于降阶观测器的NPC三电平逆变器多类故障诊断方法
技术领域
本发明涉及故障诊断技术领域,尤其涉及基于降阶观测器的NPC三电平逆变器多类故障诊断方法。
背景技术
在光伏发电系统中,逆变器发挥着重要的作用。NPC三电平逆变器具有损耗低、输出电压和电流波形谐波含量低、器件应力小、输出容量大等诸多优点,被广泛运用在光伏微电网中。尽管NPC三电平逆变器具有上述众多优点,但由于光伏发电系统所处环境复杂多变,逆变器中功率器件和传感器很容易发生故障,严重威胁光伏发电系统的安全,而NPC三电平逆变器结构复杂,具有过多电力电子器件,导致了运行时故障率的增大,不仅降低了系统运行的可靠性,而且增加了系统的运行和维护成本。为了保证NPC三电平逆变器在实际运行中的可靠性,对逆变器的开关管故障诊断的快速性和准确性提出了更高的要求。
NPC三电平逆变器的开关管故障主要可以分为短路故障和开路故障,开关管的短路故障由保护电路保护,当系统发生短路故障时保护电路迅速断开,最终会将开关管短路故障转换为开路故障,鉴于短路故障时间短,会迅速转换为开路故障,故只考虑对NPC三电平逆变器开关管开路故障进行诊断;电流传感器故障则分为增益故障、卡死故障、开路故障等,其中开路故障会导致逆变器控制系统获取不到参考电流信号,增益故障、卡死故障会导致逆变器控制系统获取到发生畸变的参考电流信号,均会导致NPC逆变器的输出电流发生严重畸变,进而导致整个逆变器系统的崩溃,因此对电流传感器的故障诊断也尤为重要。
目前对逆变器的故障诊断技术大多都是只针对开关管故障,大致可分为以下几种:
1、基于特征提取的方法。此类方法主要利用主成分分析等方法对故障的主要成分进行提取分析,并用智能分类器对故障进行诊断,如基于小波变换的方法和基于瞬时频率的方法等,具体的相关论文及专利如《A Diagnosis Algorithm for Multiple Open-Circuited Faults of Microgrid Inverters Based on Main Fault ComponentAnalysis》、《一种基于小波分析和SVM的逆变器故障诊断方法》(申请公布号CN 105095566A)、《基于瞬时频率的NPC三电平逆变器开路故障诊断方法》(申请公布号CN111077471 A)等,此类方法存在信号处理复杂度大、诊断周期长等问题。
2、基于知识的方法。其基本的理论思想是通过模拟人的思维方式去实现逆变器的故障诊断。如基于神经网络的方法和基于支持向量机的方法等,具体的相关专利文献如《一种NPC三电平光伏逆变器开路故障诊断方法》(申请公布号CN108229544A)、《基于优化支持向量机的三电平逆变器开路故障诊断方法.》(申请公布号CN110068776A)等,此类方法存在诊断计算量大、建立知识库比较难且知识库维护难度大的问题。
3、基于解析模型的方法。此类方法主要思想是建立逆变器的数学模型,将估计的系统输出与测量信息进行比较得到残差,对残差进行分析以实现变换器的故障诊断,相关专利文献如《基于区间滑模观测器的逆变器开路故障诊断方法》、《一种基于模型的多电平逆变器参数性故障诊断方法》(申请公布号CN108649600A)等,此类方法存在数学模型要求高,鲁棒性差的问题。
综上所述,现有的技术存在信号处理复杂度大、诊断周期长、计算量大、建立知识库比较难且知识库维护难度大、数学模型要求高和鲁棒性差等问题。
发明内容
本发明的目的是提供一种基于降阶观测器的NPC三电平逆变器多类故障诊断方法,解决上述现有技术中存在的问题。具体的,利用矩阵变换进行解耦,将传感器故障与系统状态分解并建立降阶观测器,避免传感器故障对观测器输出电流的影响;利用自适应阈值取代传统的固定阈值,减小故障诊断时间,提高故障诊断的鲁棒性。
为实现上述目的,本发明提供了一种基于降阶观测器的NPC三电平逆变器多类故障诊断方法,该方法涉及的NPC三电平逆变器的拓扑结构包括直流电源、两个相同的支撑电容、主逆变电路、三个相同的传感器、三个相同的电感、三个相同的电阻和控制模块;所述直流电源的直流电压记为Udc,两个所述的支撑电容分别记为支撑电容C1和支撑电容C2,支撑电容C1和支撑电容C2串联后并联在直流电源的直流正母线Q1和直流负母线Q2之间;
所述主逆变电路分为三相桥臂,三相桥臂均与直流电源并联,将三相桥臂记为k相桥臂,k表示相序,k=a,b,c;在三相桥臂中,每相桥臂由四个开关管串联组成,即主逆变电路共包含12个开关管,将12个开关管记为V,σ表示开关管的序号,σ=1,2,3,4;在三相桥臂的每相桥臂中,开关管Vk1、开关管Vk2、开关管Vk3、开关管Vk4依次串联,开关管Vk2和开关管Vk3的连接点记为主逆变电路的输出点ψk,k=a,b,c;
将所述三个相同的传感器记为传感器Sk,三个相同的电感记为电感Lk、三个相同的电阻记为Rk,k=a,b,c,所述传感器Sk的一端与主逆变电路的输出点ψk相连,另一端与电感Lk相连,电感Lk的另一端与电阻Rk相连,电阻Rk的另一端接地;
所述控制模块的输入端分别连接传感器Sa、传感器Sb、传感器Sc,所述控制模块的输出端分别连接12个开关管V
所述基于降阶观测器的NPC三电平逆变器多类故障诊断方法包括以下步骤:
步骤1,将NPC三电平逆变器记为逆变器,建立逆变器的混合逻辑动态模型,并计算k相相电压Uk的估计值k=a,b,c;
NPC三电平逆变器的混合逻辑动态模型的表达式为:
其中,为k相端电压的估计值,ξk为k相桥臂的开关函数,k=a,b,c;
k相相电压Uk的估计值的表达式为:
步骤2,通过传感器Sa、传感器Sb、传感器Sc检测逆变器的三相输出电流ia,ib,ic,采样三相输出相电压Ua,Ub,Uc,建立逆变器在故障状态下的系统状态方程,其表达式为:
其中,x为逆变器输出电流, 为x的导数,U为逆变器输出相电压,A为系数矩阵1,/>R为电阻Ra的电阻值,L为电感La的电感值,B为系数矩阵2,/>C为系数矩阵3,C=I3,I3为三阶的单位矩阵,F为故障系数矩阵,D为扰动系数矩阵,y为逆变器系统输出信号,fx为开关管故障信号,fs为传感器故障信号,d为逆变器扰动信号;
步骤3,对故障状态下的系统状态方程进行增广得到增广状态方程,其表达式为:
其中,为逆变器输出电流x的增广,/> 为逆变器输出电流x的增广/>的导数,/>为增广系数矩阵1,/>03×3为三阶的零矩阵,/>为增广系数矩阵2,/> 为增广系数矩阵3,/> 为增广系数矩阵4,/> 为增广故障系数矩阵,/>f为增广故障信号,/>
步骤4,定义一次变换矩阵T,其中/>为增广系数矩阵4的正交补的转置;定义一次变换电流/> 并将一次变换电流/>展开为/> 为一次变换电流/>的展开式1,/>为一次变换电流/>的展开式2;
将步骤3得到的增广状态方程进行第一次变换,得到一次变换状态方程为,其表达式为:
其中,为一次变换电流/>的导数,/>为一次变换系数矩阵1,/> 为一次变换系数矩阵2,/> 为一次变换系数矩阵3,/>T-1为一次变换矩阵T的逆矩阵;
步骤5,将一次变换矩阵T的逆矩阵T-1展开为 为一次变换矩阵T的逆矩阵T-1的展开式1,/>为一次变换矩阵T的逆矩阵T-1的展开式2,/>为一次变换矩阵T的逆矩阵T-1的展开式3,/>为一次变换矩阵T的逆矩阵T-1的展开式4;定义二次变换矩阵V,并将二次变换矩阵V的逆矩阵V-1展开为/> 为二次变换矩阵V的逆矩阵V-1的展开式1,/>为二次变换矩阵V的逆矩阵V-1的展开式2,/>为二次变换矩阵V的逆矩阵V-1的展开式3,/>为二次变换矩阵V的逆矩阵V-1的展开式4;
将步骤4得到的一次变换状态方程中第一个等式两边同时左乘二次变换矩阵V的逆矩阵V-1,得到二次变换状态方程,其表达式为:
其中,为二次变换系数矩阵1,/> 为二次变换系数矩阵2,并将/>展开为/> 为二次变换系数矩阵2的展开式1,/>为二次变换系数矩阵2的展开式2,/>为二次变换系数矩阵2的展开式3,/>为二次变换系数矩阵2的展开式4;/>为二次变换系数矩阵3,/>并将/>展开为/> 为二次变换系数矩阵3的展开式1,/>为二次变换系数矩阵3的展开式2;/>为二次变换系数矩阵4,/>并将/>展开为/> 为二次变换系数矩阵4的展开式1,/>为二次变换系数矩阵4的展开式2;
将二次变换状态方程展开可得到展开后的二次变换状态方程,其表达式为:
步骤6,对步骤5得到的展开后的二次变换状态方程进行第三次变换,具体步骤如下:
步骤6.1,存在正定矩阵P和增益矩阵J使矩阵不等式成立,使用LMI工具箱求出正定矩阵P和增益矩阵J,并将正定矩阵P展开为/>P1为正定矩阵P的展开式1,P2为正定矩阵P的展开式2,P3为正定矩阵P的展开式3,P4为正定矩阵P的展开式4;
步骤6.2,将展开后的二次变换状态方程方程中第一个等式两边同时左乘三次变换矩阵Z,得到三次变换状态方程,其表达式为:
其中,为P1的逆矩阵;
步骤7,设计广义降阶观测器,表达式为:
其中,为中间变量,/>为中间变量/>的导数,/>为/>的估计值,/>为/>的估计值,/>为逆变器输出电流估计值,/> 为三相输出电流ia,ib,ic的估计值,估计值,/>为传感器故障估计值,/>为逆变器扰动估计值;
步骤8,设计故障诊断自适应阈值Tth
定义K相电流形态因子Mk,其表达式为:
其中,e为自然常数,ln()为以自然常数e为底的对数函数,|ik|rms为k相输出电流ik的绝对值的均方根值;
将三相输出电流ia,ib,ic表达为k相输出电流ik,将三相输出电流ia,ib,ic的估计值表达为k相输出电流ik的估计值/>,定义K相故障诊断自适应阈值Tthk
其中,||为绝对值函数,为k相输出电流ik的估计值/>的绝对值的均方根值;
步骤9,对逆变器故障进行诊断,具体步骤如下:
步骤9.1,比较k相电流形态因子Mk和k相故障诊断自适应阈值Tthk,并进行以下判断:
若Ma<Ttha且Mb<Tthb且Mc<Tthc,则判定逆变器正常工作,故障诊断结束;
若三相中任意一相的Mk≥Tthk,将Mk≥Tthk的相记为g相,并判定逆变器g相发生故障,将g相对应的传感器故障估计值记为g相传感器故障估计值g或等于a或等于b或等于c;
步骤9.2,定义g相故障定位阈值Tg,Tg=ig max×5%,其中,ig为g相输出电流,ig max为g相输出电流ig的最大值,比较g相传感器故障估计值和g相故障定位阈值Tg,并进行以下判断:
,则判定逆变器g相开关管发生故障,进入步骤9.3;
,则判定逆变器g相传感器发生故障,进入步骤9.4;
步骤9.3,定义g相开关管故障检测特征量fg和g相故障定位特征量wgwg=sign(|ig|-Tg),其中,/>为g相输出电流ig的估计值,sign()为符号函数;
按以下条件进行开关管故障定位:
当fg=1,wg=1,则开关管Vg1发生开路故障;
当fg=1,wg=-1,则开关管Vg2发生开路故障;
当fg=-1,wg=-1,则开关管Vg3发生开路故障;
当fg=-1,wg=1,则开关管Vg4发生开路故障;
步骤9.4,定义g相传感器故障检测特征量vg和g相故障定位特征量wgwg=sign(|ig|-Tg),其中,/>为g相输出电流ig的导数,sign()为符号函数;
按以下条件进行传感器故障分类:
当vg=0,wg=1,则传感器Sg发生卡死故障;
当vg=0,wg=-1,则传感器Sg发生开路故障;
当vg≠0,wg=1,则传感器Sg发生增益故障。
优选地,步骤1所述k相桥臂的开关函数ξk按照以下方式确定:
规定电流从NPC三电平逆变器流向电感Lk为正,电流从电感Lk流向NPC三电平逆变器为负,定义逻辑变量μk,μk=1表示k相电流为正、μk=0表示k相电流为负;将开关管V的开关信号记为δ,k相桥臂的开关函数ξk的表达式如下:
其中,符号“-”表示逻辑非。
由于采用了上述故障诊断方法,与现有技术相比,本发明的有益效果为:
1、将传感器故障与系统状态分解并建立降阶观测器,避免传感器故障对观测器输出电流的影响,实现了逆变器开关管故障和传感器故障的同时诊断;
2、选取自适应阈值来故障诊断,对扰动具有抗干扰性,提高了故障检测的准确性和鲁棒性;
3、故障诊断过程不需要增加额外的传感器,降低了故障检测的成本。
附图说明
图1是本发明实施例中NPC三电平逆变器的拓扑图;
图2是本发明基于降阶观测器的NPC三电平逆变器多类故障诊断方法的流程图;
图3是本发明实施例中开关管Va1发生开路故障时a相输出电流ia、a相电流估计值、a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;
图4是本发明实施例中开关管Va1发生开路故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相开关管故障检测特征值fa和a相故障定位特征值wa的仿真波形图;
图5是本发明实施例中传感器Sa发生开路故障时a相输出电流ia、a相电流估计值、a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;
图6是本发明实施例中传感器Sa发生开路故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相传感器故障检测特征值va和a相故障定位特征值wa的仿真波形图;
图7是本发明实施例中传感器Sa发生卡死故障时a相输出电流ia、a相电流估计值、a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;
图8是本发明实施例中传感器Sa发生卡死故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相传感器故障检测特征值va和a相故障定位特征值wa的仿真波形图;
图9是本发明实施例中传感器Sa发生增益故障时a相输出电流ia、a相电流估计值、a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;
图10是本发明实施例中传感器Sa发生增益故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相传感器故障检测特征值va和a相故障定位特征值wa的仿真波形图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整的描述。
图1是本发明实施例中的NPC三电平逆变器的拓扑图。由图可见本发明方法涉及的NPC三电平逆变器的拓扑结构包括直流电源、两个相同的支撑电容、主逆变电路、三个相同的传感器、三个相同的电感、三个相同的电阻和控制模块;所述直流电源的直流电压记为Udc,两个所述的支撑电容分别记为支撑电容C1和支撑电容C2,支撑电容C1和支撑电容C2串联后并联在直流电源的直流正母线Q1和直流负母线Q2之间。
所述主逆变电路分为三相桥臂,三相桥臂均与直流电源并联,将三相桥臂记为k相桥臂,k表示相序,k=a,b,c;在三相桥臂中,每相桥臂由四个开关管串联组成,即主逆变电路共包含12个开关管,将12个开关管记为V,σ表示开关管的序号,σ=1,2,3,4;在三相桥臂的每相桥臂中,开关管Vk1、开关管Vk2、开关管Vk3、开关管Vk4依次串联,开关管Vk2和开关管Vk3的连接点记为主逆变电路的输出点ψk,k=a,b,c。
将所述三个相同的传感器记为传感器Sk,三个相同的电感记为电感Lk、所述三个相同的电阻记为Rk,k=a,b,c,所述传感器Sk的一端与主逆变电路的输出点ψk相连,另一端与电感Lk相连,电感Lk的另一端与电阻Rk相连,电阻Rk的另一端接地。
所述控制模块的输入端分别连接传感器Sa、传感器Sb、传感器Sc,所述控制模块的输出端分别连接12个开关管V
在本实施例中Udc=220V。
在图1中,点O为支撑电容C1和支撑电容C2的公共节点。从图1可见,在三相桥臂中,每相桥臂还包括两个二极管,即三相桥臂共包括6个二极管,将六个二极管记为Dkh,h表示二极管的序号,h=1,2。具体的,二极管Dk1的阳极连接中性点O,二极管Da1的阴极连接开关管Va1的集电极,二极管Db1的阴极连接开关管Vb1的集电极,二极管Dc1的阴极连接开关管Vc1的集电极;二极管Dk2的阴极连接中性点O,二极管Da2的阳极连接开关管Va3的发射极,二极管Db2的阳极连接开关管的发射极,二极管Dc2的阳极连接开关管Vc3的发射极。
图2是本发明基于降阶观测器的NPC三电平逆变器多类故障诊断方法的流程图,由图2可见,基于降阶观测器的NPC三电平逆变器多类故障诊断方法包括以下步骤:
步骤1,将NPC三电平逆变器记为逆变器,建立逆变器的混合逻辑动态模型,并计算k相相电压Uk的估计值,k=a,b,c;
NPC三电平逆变器的混合逻辑动态模型的表达式为:
其中,为k相端电压的估计值,ξk为k相桥臂的开关函数,k=a,b,c;
k相相电压Uk的估计值的表达式为:
在本实施例中,所述k相桥臂的开关函数ξk按照以下方式确定:
规定电流从NPC三电平逆变器流向电感Lk为正,电流从电感Lk流向NPC三电平逆变器为负,定义逻辑变量μk,μk=1表示k相电流为正、μk=0表示k相电流为负;将开关管V的开关信号记为δ,k相桥臂的开关函数ξk的表达式如下:
其中,符号“-”表示逻辑非。
步骤2,通过传感器Sa、传感器Sb、传感器Sc检测逆变器的三相输出电流ia,ib,ic,采样三相输出相电压Ua,Ub,Uc,建立逆变器在故障状态下的系统状态方程,其表达式为:
其中,x为逆变器输出电流, 为x的导数,U为逆变器输出相电压,A为系数矩阵1,/>R为电阻Ra的电阻值,L为电感La的电感值,B为系数矩阵2,/>C为系数矩阵3,C=I3,I3为三阶的单位矩阵,F为故障系数矩阵,D为扰动系数矩阵,y为逆变器系统输出信号,fx为开关管故障信号,fs为传感器故障信号,d为逆变器扰动信号。
在本实施例中,R=10Ω,L=80mH,d=0.01sin(100πt)。
步骤3,对故障状态下的系统状态方程进行增广得到增广状态方程,其表达式为:
其中,为逆变器输出电流x的增广,/> 为逆变器输出电流x的增广/>的导数,/>为增广系数矩阵1,/>03×3为三阶的零矩阵,/>为增广系数矩阵2,/> 为增广系数矩阵3,/> 为增广系数矩阵4,/> 为增广故障系数矩阵,/>f为增广故障信号,/>
步骤4,定义一次变换矩阵T,其中/>为增广系数矩阵4的正交补的转置;定义一次变换电流/> 并将一次变换电流/>展开为/> 为一次变换电流/>的展开式1,/>为一次变换电流/>的展开式2;
将步骤3得到的增广状态方程进行第一次变换,得到一次变换状态方程为,其表达式为:
其中,为一次变换电流/>的导数,/>为一次变换系数矩阵1,/> 为一次变换系数矩阵2,/> 为一次变换系数矩阵3,/>T-1为一次变换矩阵T的逆矩阵。
步骤5,将一次变换矩阵T的逆矩阵T-1展开为 为一次变换矩阵T的逆矩阵T-1的展开式1,/>为一次变换矩阵T的逆矩阵T-1的展开式2,/>为一次变换矩阵T的逆矩阵T-1的展开式3,/>为一次变换矩阵T的逆矩阵T-1的展开式4;定义二次变换矩阵V,并将二次变换矩阵V的逆矩阵V-1展开为/> 为二次变换矩阵V的逆矩阵V-1的展开式1,/>为二次变换矩阵V的逆矩阵V-1的展开式2,/>为二次变换矩阵V的逆矩阵V-1的展开式3,/>为二次变换矩阵V的逆矩阵V-1的展开式4;
将步骤4得到的一次变换状态方程中第一个等式两边同时左乘二次变换矩阵V的逆矩阵V-1,得到二次变换状态方程,其表达式为:
其中,为二次变换系数矩阵1,/> 为二次变换系数矩阵2,并将/>展开为/> 为二次变换系数矩阵2的展开式1,/>为二次变换系数矩阵2的展开式2,/>为二次变换系数矩阵2的展开式3,/>为二次变换系数矩阵2的展开式4;/>为二次变换系数矩阵3,/>并将/>展开为/> 为二次变换系数矩阵3的展开式1,/>为二次变换系数矩阵3的展开式2;/>为二次变换系数矩阵4,/>并将/>展开为/> 为二次变换系数矩阵4的展开式1,/>为二次变换系数矩阵4的展开式2;
将二次变换状态方程展开可得到展开后的二次变换状态方程,其表达式为:
步骤6,对步骤5得到的展开后的二次变换状态方程进行第三次变换,具体步骤如下:
步骤6.1,存在正定矩阵P和增益矩阵J使矩阵不等式成立,使用LMI工具箱求出正定矩阵P和增益矩阵J,并将正定矩阵P展开为/>,P1为正定矩阵P的展开式1,P2为正定矩阵P的展开式2,P3为正定矩阵P的展开式3,P4为正定矩阵P的展开式4;
在本实施例中,正定矩阵P和增益矩阵J的具体数字如下:
步骤6.2,将展开后的二次变换状态方程方程中第一个等式两边同时左乘三次变换矩阵Z,得到三次变换状态方程,其表达式为:
其中,为P1的逆矩阵。
步骤7,设计广义降阶观测器,表达式为:
其中,为中间变量,/>为中间变量/>的导数,/>为/>的估计值,/>为/>的估计值,/>为逆变器输出电流估计值,/> 为三相输出电流ia,ib,ic的估计值,估计值,/>为传感器故障估计值,/>为逆变器扰动估计值。
步骤8,设计故障诊断自适应阈值Tth
定义K相电流形态因子Mk,其表达式为:
其中,e为自然常数,ln()为以自然常数e为底的对数函数,|ik|rms为k相输出电流ik的绝对值的均方根值;
将三相输出电流ia,ib,ic表达为k相输出电流ik,将三相输出电流ia,ib,ic的估计值表达为k相输出电流ik的估计值/>定义K相故障诊断自适应阈值Tthk
其中,||为绝对值函数,为k相输出电流ik的估计值/>的绝对值的均方根值。
步骤9,对逆变器故障进行诊断,具体步骤如下:
步骤9.1,比较k相电流形态因子Mk和k相故障诊断自适应阈值Tthk,并进行以下判断:
若Ma<Ttha且Mb<Tthb且Mc<Tthc,则判定逆变器正常工作,故障诊断结束;
若三相中任意一相的Mk≥Tthk,将Mk≥Tthk的相记为g相,并判定逆变器g相发生故障,将g相对应的传感器故障估计值记为g相传感器故障估计值g或等于a或等于b或等于c;
步骤9.2,定义g相故障定位阈值Tg,Tg=ig max×5%,其中,ig为g相输出电流,ig max为g相输出电流ig的最大值,比较g相传感器故障估计值和g相故障定位阈值Tg,并进行以下判断:
则判定逆变器g相开关管发生故障,进入步骤9.3;
则判定逆变器g相传感器发生故障,进入步骤9.4;
步骤9.3,定义g相开关管故障检测特征量fg和g相故障定位特征量wgwg=sign(|ig|-Tg),其中,/>为g相输出电流ig的估计值,sign()为符号函数;
按以下条件进行开关管故障定位:
当fg=1,wg=1,则开关管Vg1发生开路故障;
当fg=1,wg=-1,则开关管Vg2发生开路故障;
当fg=-1,wg=-1,则开关管Vg3发生开路故障;
当fg=-1,wg=1,则开关管Vg4发生开路故障;
步骤9.4,定义g相传感器故障检测特征量vg和g相故障定位特征量wgwg=sign(|ig|-Tg),其中,/>为g相输出电流ig的导数,sign()为符号函数;
按以下条件进行传感器故障分类:
当vg=0,wg=1,则传感器Sg发生卡死故障;
当vg=0,wg=-1,则传感器Sg发生开路故障;
当vg≠0,wg=1,则传感器Sg发生增益故障。
通过仿真对本发明进行了验证。
图3是本发明实施例中开关管Va1发生开路故障时a相输出电流ia、a相电流估计值a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;由该图可见,在0.264秒后a相输出电流ia发生较大的变化,a相电流估计值/>不变,a相传感器故障估计值/>不超过故障定位阈值Ta
图4是本发明实施例中开关管Va1发生开路故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相开关管故障检测特征值fa和a相故障定位特征值wa的仿真波形图;由该图可见,在0.268秒后,a相电流形态因子Ma超过a相故障诊断自适应阈值Ttha,此时a相开关管故障检测特征值fa=1、a相故障定位特征值wa=1。
图5是本发明实施例中传感器Sa发生开路故障时a相输出电流ia、a相电流估计值a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;由该图可见,在0.264秒后a相输出电流ia恒定为0,a相电流估计值/>不变,a相传感器故障估计值/>超过故障定位阈值Ta
图6是本发明实施例中传感器Sa发生开路故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相传感器故障检测特征值va和a相故障定位特征值wa的仿真波形图;由该图可见,在0.268秒后,a相电流形态因子Ma超过a相故障诊断自适应阈值Ttha,此时a相传感器故障检测特征值va=0、a相故障定位特征值wa=-1。
图7是本发明实施例中传感器Sa发生卡死故障时a相输出电流ia、a相电流估计值、a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;由该图可见,在0.264秒后a相输出电流ia恒定为2,a相电流估计值/>不变,a相传感器故障估计值/>超过故障定位阈值Ta
图8是本发明实施例中传感器Sa发生卡死故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相传感器故障检测特征值va和a相故障定位特征值wa的仿真波形图;由该图可见,在0.268秒后,a相电流形态因子Ma超过a相故障诊断自适应阈值Ttha,此时a相传感器故障检测特征值va=0、a相故障定位特征值wa=1。
图9是本发明实施例中传感器Sa发生增益故障时a相输出电流ia、a相电流估计值、a相传感器故障估计值/>和a相故障定位阈值Ta的仿真波形图;由该图可见,在0.264秒后a相输出电流ia发生较大变化,a相电流估计值/>不变,a相传感器故障估计值/>超过故障定位阈值Ta
图10是本发明实施例中传感器Sa发生增益故障时a相电流形态因子Ma和a相故障诊断自适应阈值Ttha、a相传感器故障检测特征值va和a相故障定位特征值wa的仿真波形图;由该图可见,在0.268秒后,a相电流形态因子Ma超过a相故障诊断自适应阈值Ttha,此时a相传感器故障检测特征值va≠0、a相故障定位特征值wa=1。

Claims (2)

1.一种基于降阶观测器的NPC三电平逆变器多类故障诊断方法,该方法涉及的NPC三电平逆变器的拓扑结构包括直流电源、两个相同的支撑电容、主逆变电路、三个相同的传感器、三个相同的电感、三个相同的电阻和控制模块;所述直流电源的直流电压记为Udc,两个所述的支撑电容分别记为支撑电容C1和支撑电容C2,支撑电容C1和支撑电容C2串联后并联在直流电源的直流正母线Q1和直流负母线Q2之间;
所述主逆变电路分为三相桥臂,三相桥臂均与直流电源并联,将三相桥臂记为k相桥臂,k表示相序,k=a,b,c;在三相桥臂中,每相桥臂由四个开关管串联组成,即主逆变电路共包含12个开关管,将12个开关管记为V,σ表示开关管的序号,σ=1,2,3,4;在三相桥臂的每相桥臂中,开关管Vk1、开关管Vk2、开关管Vk3、开关管Vk4依次串联,开关管Vk2和开关管Vk3的连接点记为主逆变电路的输出点ψk,k=a,b,c;
将所述三个相同的传感器记为传感器Sk、三个相同的电感记为电感Lk、三个相同的电阻记为Rk,k=a,b,c,所述传感器Sk的一端与主逆变电路的输出点ψk相连,另一端与电感Lk相连,电感Lk的另一端与电阻Rk相连,电阻Rk的另一端接地;
所述控制模块的输入端分别连接传感器Sa、传感器Sb、传感器Sc,所述控制模块的输出端分别连接12个开关管V
其特征在于,所述基于降阶观测器的NPC三电平逆变器多类故障诊断方法包括以下步骤:
步骤1,将NPC三电平逆变器记为逆变器,建立逆变器的混合逻辑动态模型,并计算k相相电压Uk的估计值k=a,b,c;
NPC三电平逆变器的混合逻辑动态模型的表达式为:
其中,为k相端电压的估计值,ξk为k相桥臂的开关函数,k=a,b,c;
k相相电压Uk的估计值的表达式为:
步骤2,通过传感器Sa、传感器Sb、传感器Sc检测逆变器的三相输出电流ia,ib,ic,采样三相输出相电压Ua,Ub,Uc,建立逆变器在故障状态下的系统状态方程,其表达式为:
其中,x为逆变器输出电流, 为x的导数,U为逆变器输出相电压,/>A为系数矩阵1,/>R为电阻Ra的电阻值,L为电感La的电感值,B为系数矩阵2,C为系数矩阵3,C=I3,I3为三阶的单位矩阵,F为故障系数矩阵,D为扰动系数矩阵,y为逆变器系统输出信号,fx为开关管故障信号,fs为传感器故障信号,d为逆变器扰动信号;
步骤3,对故障状态下的系统状态方程进行增广得到增广状态方程,其表达式为:
其中,为逆变器输出电流x的增广,/> 为逆变器输出电流x的增广/>的导数,/>为增广系数矩阵1,/>03×3为三阶的零矩阵,/>为增广系数矩阵2, 为增广系数矩阵3,/> 为增广系数矩阵4,/> 为增广故障系数矩阵,/>f为增广故障信号,/>
步骤4,定义一次变换矩阵T,其中/>为增广系数矩阵4的正交补的转置;定义一次变换电流/>并将一次变换电流/>展开为/> 为一次变换电流/>的展开式1,/>为一次变换电流/>的展开式2;
将步骤3得到的增广状态方程进行第一次变换,得到一次变换状态方程为,其表达式为:
其中,为一次变换电流/>的导数,/>为一次变换系数矩阵1,/> 为一次变换系数矩阵2,/> 为一次变换系数矩阵3,/>T-1为一次变换矩阵T的逆矩阵;
步骤5,将一次变换矩阵T的逆矩阵T-1展开为T1 -1为一次变换矩阵T的逆矩阵T-1的展开式1,T2 -1为一次变换矩阵T的逆矩阵T-1的展开式2,T3 -1为一次变换矩阵T的逆矩阵T-1的展开式3,T4 -1为一次变换矩阵T的逆矩阵T-1的展开式4;定义二次变换矩阵V,并将二次变换矩阵V的逆矩阵V-1展开为/>V1 -1为二次变换矩阵V的逆矩阵V-1的展开式1,V2 -1为二次变换矩阵V的逆矩阵V-1的展开式2,V3 -1为二次变换矩阵V的逆矩阵V-1的展开式3,V4 -1为二次变换矩阵V的逆矩阵V-1的展开式4;
将步骤4得到的一次变换状态方程中第一个等式两边同时左乘二次变换矩阵V的逆矩阵V-1,得到二次变换状态方程,其表达式为:
其中,为二次变换系数矩阵1,/> 为二次变换系数矩阵2,/>并将/>展开为/> 为二次变换系数矩阵2的展开式1,/>为二次变换系数矩阵2的展开式2,/>为二次变换系数矩阵2的展开式3,/>为二次变换系数矩阵2的展开式4;为二次变换系数矩阵3,/>并将/>展开为/> 为二次变换系数矩阵3的展开式1,/>为二次变换系数矩阵3的展开式2;/>为二次变换系数矩阵4,/>并将/>展开为/> 为二次变换系数矩阵4的展开式1,/>为二次变换系数矩阵4的展开式2;
将二次变换状态方程展开可得到展开后的二次变换状态方程,其表达式为:
步骤6,对步骤5得到的展开后的二次变换状态方程进行第三次变换,具体步骤如下:
步骤6.1,存在正定矩阵P和增益矩阵J使矩阵不等式成立,使用LMI工具箱求出正定矩阵P和增益矩阵J,并将正定矩阵P展开为/>P1为正定矩阵P的展开式1,P2为正定矩阵P的展开式2,P3为正定矩阵P的展开式3,P4为正定矩阵P的展开式4;
步骤6.2,将展开后的二次变换状态方程方程中第一个等式两边同时左乘三次变换矩阵Z,得到三次变换状态方程,其表达式为:
其中,P1 -1为P1的逆矩阵;
步骤7,设计广义降阶观测器,表达式为:
其中,为中间变量,/>为中间变量/>的导数,/>为/>的估计值,/>为/>的估计值,/>为逆变器输出电流估计值,/> 为三相输出电流ia,ib,ic的估计值,/>为传感器故障估计值,/>为逆变器扰动估计值;
步骤8,设计故障诊断自适应阈值Tth
定义K相电流形态因子Mk,其表达式为:
其中,e为自然常数,ln()为以自然常数e为底的对数函数,|ik|rms为k相输出电流ik的绝对值的均方根值;
将三相输出电流ia,ib,ic表达为k相输出电流ik,将三相输出电流ia,ib,ic的估计值表达为k相输出电流ik的估计值/>定义K相故障诊断自适应阈值Tthk
其中,| |为绝对值函数,为k相输出电流ik的估计值/>的绝对值的均方根值;所述k相输出电流ik的估计值/>为三相输出电流估计值
步骤9,对逆变器故障进行诊断,具体步骤如下:
步骤9.1,比较k相电流形态因子Mk和k相故障诊断自适应阈值Tthk,并进行以下判断:
若Ma<Ttha且Mb<Tthb且Mc<Tthc,则判定逆变器正常工作,故障诊断结束;
若三相中任意一相的Mk≥Tthk,将Mk≥Tthk的相记为g相,并判定逆变器g相发生故障,将g相对应的传感器故障估计值记为g相传感器故障估计值g或等于a或等于b或等于c;
步骤9.2,定义g相故障定位阈值Tg,Tg=igmax×5%,其中,ig为g相输出电流,igmax为g相输出电流ig的最大值,比较g相传感器故障估计值和g相故障定位阈值Tg,并进行以下判断:
则判定逆变器g相开关管发生故障,进入步骤9.3;
则判定逆变器g相传感器发生故障,进入步骤9.4;
步骤9.3,定义g相开关管故障检测特征量fg和g相故障定位特征量wgwg=sign(|ig|-Tg),其中,/>为g相输出电流ig的估计值,sign()为符号函数;
按以下条件进行开关管故障定位:
当fg=1,wg=1,则开关管Vg1发生开路故障;
当fg=1,wg=-1,则开关管Vg2发生开路故障;
当fg=-1,wg=-1,则开关管Vg3发生开路故障;
当fg=-1,wg=1,则开关管Vg4发生开路故障;
步骤9.4,定义g相传感器故障检测特征量vg和g相故障定位特征量wgwg=sign(|ig|-Tg),其中,/>为g相输出电流ig的导数,sign()为符号函数;
按以下条件进行传感器故障分类:
当vg=0,wg=1,则传感器Sg发生卡死故障;
当vg=0,wg=-1,则传感器Sg发生开路故障;
当vg≠0,wg=1,则传感器Sg发生增益故障。
2.根据权利要求1所述基于降阶观测器的NPC三电平逆变器多类故障诊断方法,其特征在于,步骤1所述k相桥臂的开关函数ξk按照以下方式确定:
规定电流从NPC三电平逆变器流向电感Lk为正,电流从电感Lk流向NPC三电平逆变器为负,定义逻辑变量μk,μk=1表示k相电流为正、μk=0表示k相电流为负;将开关管V的开关信号记为δ,k相桥臂的开关函数ξk的表达式如下:
其中,符号“-”表示逻辑非。
CN202111344797.3A 2021-11-12 2021-11-12 基于降阶观测器的npc三电平逆变器多类故障诊断方法 Active CN114172401B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111344797.3A CN114172401B (zh) 2021-11-12 2021-11-12 基于降阶观测器的npc三电平逆变器多类故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111344797.3A CN114172401B (zh) 2021-11-12 2021-11-12 基于降阶观测器的npc三电平逆变器多类故障诊断方法

Publications (2)

Publication Number Publication Date
CN114172401A CN114172401A (zh) 2022-03-11
CN114172401B true CN114172401B (zh) 2023-08-29

Family

ID=80478930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111344797.3A Active CN114172401B (zh) 2021-11-12 2021-11-12 基于降阶观测器的npc三电平逆变器多类故障诊断方法

Country Status (1)

Country Link
CN (1) CN114172401B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117491934B (zh) * 2023-12-29 2024-03-12 锦浪科技股份有限公司 一种传感器故障检测方法、系统及光伏系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109756169A (zh) * 2017-11-06 2019-05-14 操纵技术Ip控股公司 永磁dc驱动转向系统的电流传感器故障缓解
CN110649846A (zh) * 2019-09-29 2020-01-03 湖南科技大学 基于滑模观测器的双馈感应发电机故障检测方法
CN111077471A (zh) * 2019-12-30 2020-04-28 合肥工业大学 基于瞬时频率的npc三电平逆变器开路故障诊断方法
CN113534035A (zh) * 2021-06-30 2021-10-22 合肥工业大学 新能源电动汽车驱动系统电流传感器微小故障诊断方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109905048B (zh) * 2017-12-08 2021-01-26 台达电子企业管理(上海)有限公司 三电平电路中点电压平衡控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109756169A (zh) * 2017-11-06 2019-05-14 操纵技术Ip控股公司 永磁dc驱动转向系统的电流传感器故障缓解
CN110649846A (zh) * 2019-09-29 2020-01-03 湖南科技大学 基于滑模观测器的双馈感应发电机故障检测方法
CN111077471A (zh) * 2019-12-30 2020-04-28 合肥工业大学 基于瞬时频率的npc三电平逆变器开路故障诊断方法
CN113534035A (zh) * 2021-06-30 2021-10-22 合肥工业大学 新能源电动汽车驱动系统电流传感器微小故障诊断方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
永磁同步电机无传感器控制技术研究;庞晴晴;《中国优秀硕士学位论文全文数据库》;全文 *

Also Published As

Publication number Publication date
CN114172401A (zh) 2022-03-11

Similar Documents

Publication Publication Date Title
Zhang et al. Fault diagnosis of wind turbine power converter considering wavelet transform, feature analysis, judgment and BP neural network
Huang et al. Multiple open-circuit fault diagnosis based on multistate data processing and subsection fluctuation analysis for photovoltaic inverter
CN111562517B (zh) 一种npc三电平逆变器开关管开路故障诊断方法
Devi et al. Detection of transmission line faults using discrete wavelet transform
Jain et al. Fault classification of double circuit transmission line using artificial neural network
Manohar et al. Enhancing resilience of PV-fed microgrid by improved relaying and differentiating between inverter faults and distribution line faults
Yuan et al. Open-circuit fault diagnosis of NPC inverter based on improved 1-D CNN network
CN102567783A (zh) 并联混合型电能质量调节器专家故障分析与诊断方法
CN113075585B (zh) 一种npc三电平逆变器复合开路故障诊断方法
Kou et al. Fault diagnosis for power electronics converters based on deep feedforward network and wavelet compression
CN114172401B (zh) 基于降阶观测器的npc三电平逆变器多类故障诊断方法
CN113610119B (zh) 基于卷积神经网络识别输电线路发展性故障的方法
CN112731193B (zh) 基于滑模观测器的npc逆变器多类故障诊断方法
CN113702767A (zh) 基于小波滑窗能量的孤岛直流微电网故障诊断方法
An et al. IGBT open circuit fault diagnosis method for a modular multilevel converter based on PNN-MD
CN112731103B (zh) 一种双级矩阵变换器的故障诊断方法
Patil et al. Improved fault detection and location scheme for photovoltaic system
CN112928772B (zh) 一种三相两电平光伏逆变器直流过压故障成因定位方法
Li et al. Fault diagnosis and fault-tolerant control of photovoltaic micro-inverter
Wang et al. A Data-Driven Open-Circuit Fault Diagnosis Method for Three-Phase T-type Three-level Inverters
CN113315100A (zh) 一种基于卷积神经网络的微电网保护方法与系统
Zhan et al. Fault diagnosis of DAB converter based on wavelet packet-RBF neural network
CN111398787A (zh) 复杂工况下的三相电压型pwm整流电路故障诊断方法
CN114914888B (zh) 一种单相三电平整流器多故障估计方法
CN115267474B (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