CN105160047A - 基于ybco超导带材的电阻型超导限流器数字建模仿真方法 - Google Patents
基于ybco超导带材的电阻型超导限流器数字建模仿真方法 Download PDFInfo
- Publication number
- CN105160047A CN105160047A CN201510204556.7A CN201510204556A CN105160047A CN 105160047 A CN105160047 A CN 105160047A CN 201510204556 A CN201510204556 A CN 201510204556A CN 105160047 A CN105160047 A CN 105160047A
- Authority
- CN
- China
- Prior art keywords
- superconducting tape
- resistance
- ybco
- superconducting
- electric current
- 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
Landscapes
- Superconductors And Manufacturing Methods Therefor (AREA)
Abstract
一种基于YBCO超导带材的电阻型超导限流器数字建模仿真方法,在建立YBCO超导带材等效结构模型、YBCO超导带材等效电路模型、YBCO超导带材热传导模型和电阻型超导限流器电路模型的基础上,根据给定限流器的电路参数、超导带材的结构参数和初始运行条件,计算线路电流和超导带材电流;根据YBCO超导带材的热传导模型,计算超导带材的温度;根据YBCO超导带材等效电路模型、超导带材温度和电流,计算超导带材的电阻;根据电阻型超导限流器的电路模型,计算线路电流和超导带材电流,实现电阻型超导限流器的仿真计算。本发明可用于多单元和复杂结构的超导限流器。
Description
技术领域
本发明涉及一种超导限流器数字建模仿真方法。
背景技术
随着国民经济的快速发展,社会对电力的需求不断增加,带动了电力***的不断发展,单机和发电厂容量、变电所容量、城市和工业中心负荷不断增加,就使得电力***之间互联,各级电网中的短路电流水平不断提高,短路故障对电力***及其相连的电气设备的破坏性也越来越大。然而,大电网的暂态稳定性问题比较突出,其中最重要的原因之一是由于常规电力技术缺乏行之有效的短路故障电流限制技术。
电阻型超导限流器是一种具有远大发展前途的限流器,它利用超导体从超导态到正常态的转变,即从无阻状态到有阻状态的变化来限流,具有自动故障辨识和自动响应故障的特点。电阻型超导限流器的限流故障是非线性变化过程,其变化过程涉及到电磁场和热场的相互耦合。对于电阻型超导限流器的数字建模仿真方法,在文献“ThermalandElectricalAnalysisofCoatedConductorUnderACOver-Current(IEEETRANSACTIONSONAPPLIEDSUPERCONDUCTIVITY,VOL.17,NO.2,JUNE2007)”中有过YBCO超导带材的电磁场和热场相互耦合关系进了建模研究,分析了YBCO超导带材在经历超导态、磁通流阻态和正常态过程中,所产生的电阻及超导的温度等。然而,该数字建模仿真方法简单地把YBCO超导带材在超导态、磁通流阻态和正常态的电阻分段处理,使得在故障限流的暂态过程中,建模无法平滑过渡,经常产生奇异点,函数无法收敛;同时,由于参数多,给分析问题带来了极大的不便。
从横截面来看,基于YBCO超导带材从上到下呈层状结构,分为6层,上表面铜层,依次为银层、YBCO层、缓冲层、哈氏合金基底层和下表面的铜层,如图1所示。当超导带材处于超导态时,YBCO层的电阻为零,电流通过YBCO层而导通;当超导带材处于正常态时,YBCO层因失超而产生电阻,电流超导带材的各层中按电阻分配而导通。
对于基于YBCO超导带材的电阻型超导限流器,不可避免地需要进行超导体从超导态、磁通流阻态和正常态的电磁场和热场的耦合关系分析,只有准确分析此过程,才能准确计算出YBCO超导带材所产生的电阻,从而实现限流器的设计。
发明内容
本发明的目的是克服已有技术的不足,提出一种基于YBCO超导带材的电阻型超导限流器数字建模仿真方法。本发明不但能够保证建模的准确性,而且避免了超导限流器暂态分析过程中的函数收敛问题,提高了***仿真效率。
本发明采用的技术方案:
本发明在建立YBCO超导带材等效结构模型、YBCO超导带材等效电路模型、YBCO超导带材热传导模型和电阻型超导限流器电路模型的基础上,根据给定电阻型超导限流器的电路参数、超导带材的结构参数和初始运行条件,在超导无感线圈电阻Rsc为零的前提下,计算线路电流的初始值和超导带材电流的初始值;根据YBCO超导带材的热传导模型计算超导带材的温度;根据YBCO超导带材等效电路模型、超导带材温度和电流,计算超导带材的电阻;根据电阻型超导限流器的电路模型,计算并反馈线路电流和超导带材电流,实现电阻型超导限流器的建模仿真。
本发明建模仿真方法的具体步骤如下:
步骤1.建立YBCO超导带材等效结构模型;
建立YBCO超导带材等效结构模型时,为便于分析YBCO超导带材各层在失超过程中的电阻情况,采用以下简化方法:
因YBCO超导带材的缓冲层的厚度很小,电阻很大,可以忽略。同时,又因YBCO超导带材的上表面铜层和下表面铜层的物理特性相同,可以合二为一进行分析。因此,所建立的YBCO超导带材等效结构模型简化分为4层:表面铜层、银层、YBCO层和哈氏合金基底层。表面铜层包括上表面铜层和下表面铜层。
步骤2.建立电阻型超导限流器的电路模型,给定电阻型超导限流器的电路参数:线路等效电感Ls、线路等效电阻r和负载电阻Rload,以及交流电源电压Us,在超导无感线圈电阻Rsc等于零的前提下,计算线路电流初始值I0(t)和超导带材电流初始值Is0(t);
电阻型超导限流器包括交流电源Us、线路等效感抗X、线路等效电阻r、断路器Br、超导无感线圈Rsc和负载电阻Rload。交流电源Us、线路等效感抗X、线路等效电阻r、断路器Br、超导无感线圈电阻Rsc和负载电阻Rload依次串联,交流电源Us和负载电阻Rload的一端接地;断路器Br和超导无感线圈电阻Rsc连接在第一连接点A上,超导无感线圈电阻Rsc和负载电阻Rload连接在第二连接点B上。
根据全电路欧姆定律,当超导无感线圈电阻Rsc为零时,电阻型超导限流器的电流与电压关系为:
Us=I0X+I0r+I0Rload(1)
式中:Rload为负载电阻,I0为线路电流,Us为交流电源电压,X为线路等效感抗,r为线路等效电阻。
其中,交流电源电压Us、线路等效感抗X分别表示为:
X=j2πfLs(3)
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,j为虚函数符号。
线路电流的初始值为:
超导无感线圈的电阻为Rsc(T),超导无感线圈一般由m根超导带材并联组成,m≥1,根据分流定律:
Is0(t)=I0(t)/m(5)
式中,I0(t)和Is0(t)分别为线路电流初始值和超导带材的电流初始值。
步骤3.建立YBCO超导带材的热传导模型,给定超导带材的结构参数:超导带材的宽度w、厚度d、长度le,所述的厚度d包括表面铜层厚度d1、银层的厚度d2、YBCO层的厚度d3和哈氏合金基底层厚度d4;给定超导无感线圈的初始运行条件:工作温度Top,计算超导带材的温度T;
(1)YBCO超导带材直接在液氮浸泡环境中冷却,根据热平衡方程,沿超导带材的长度方向,一维热传导方程为:
式中,Kcm(T)为热传导系数,Vcm和Ccm分别为YBCO超导带材的密度与比热容,gj(T)和Wcool(T)分别是超导带材的焦耳热和散失热量。
YBCO超导带材的密度Vcm与比热容Ccm满足式:
其中,ν1,ν2,ν3,ν4分别为表面铜层、银层、YBCO层和哈氏合金基底层的密度。d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,d2为银层的厚度,d3为YBCO层的厚度,d4为哈氏合金基底层的厚度,d为YBCO超导带材的厚度。C1,C2,C3,C4分别为表面铜层、银层、YBCO层和哈氏合金基底层的比热容。各种密度和比热容均可查阅手册获得,比如《超导电力技术基础》(科学出版社,2011年)。
其中,K1,K2,K3,K4分别为表面铜层,银层,YBCO层和哈氏合金基底层的热导率。各种材料的热导率均可查阅手册获得,比如《超导电力技术基础》(科学出版社,2011年)。
(2)根据焦耳定律,超导带材的焦耳热
其中,US(t)为超导带材的电压,IS(t)为超导带材的电流。
(3)超导带材浸泡在液氮中,根据经验公式,散失热量完全由液氮带走,超导带材的散失热量Wcool(T)为:
Wcool(T)=hA(T)(10)
其中,A为超导带材与液氮的接触面积,即超导带材的表面积;h为液氮传热系数,与超导带材和液氮的温度差ΔT(T-Top)有关,根据实验可知,液氮传热系数h的取值。对应于不同的温差ΔT(T-Top),热传递过程有对流、核沸腾、过渡态和膜沸腾4种状态,不同的状态对应不同的液氮传热系数h。液氮传热系数h的与拟合结果:
步骤4.建立YBCO超导带材等效电路模型,依据超导带材温度T和超导带材电流Is(t),计算YBCO超导带材的电阻rsc;
YBCO超导带材的等效电路为4个电阻并联结构。第一电阻r1为表面铜层电阻,第二电阻r2为银层的电阻,第三电阻r3为YBCO层的电阻,第四电阻r4为哈氏合金基底层的电阻。
根据所建立的YBCO超导带材等效电路模型,并根据超导带材温度和电流,根据欧姆定律和电路原理,计算YBCO超导带材的等效电路的电阻:
(1)第一电阻r1为表面铜层电阻,为铜材料制作。第一电阻r1是超导带材温度T的函数。根据欧姆定律可得:
其中,ρ1(T)为铜的电阻率,是超导带材温度T的函数,w为YBCO超导带材的宽度,d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,le为超导带材的长度。
(2)第二电阻r2为银层电阻,为银材料制作。第二电阻r2是超导带材温度T的函数。根据欧姆定律可得:
其中,ρ2(T)为银的电阻率,w为YBCO超导带材的宽度,d2为银层的厚度,le为超导带材的长度。
(3)第三电阻r3为YBCO层电阻。根据超导带材的电阻变化规律,第三电阻r3的电阻率ρ3(T)是超导带材温度T的函数,通过对YBCO超导带材的特性测试,并采用数值拟合的方法可得第三电阻r3的电阻率ρ3(T):
ρ3(T)=ρ31(T)+ρ32(T)(14)
其中,ρ31(T)和ρ32(T)由分段函数拟合而成:
其中,Jc(T)是超导带材临界电流密度Jc随温度T变化的函数:
Jc(T)=Jc0[(Tc-T)/(Tc-Top)]1.5(17)
其中,Jc为超导带材临界电流密度。Jc0=2×106A/cm2,为77K下的临界电流密度;Tc=92K,为YBCO的临界温度,Top为工作温度,在液氮池中为77K。参数n1=3;n2=20;γ=2,均为根据超导带材的特性而得到的拟合参数。针对ρ31(T)和ρ32(T)为零的情况,在实际计算过程中,可以取为一个非常小的数据,如10-19Ω··cm等,避免运算出错。
根据欧姆定律,超导带材的电流Is表示为:
IS(T)=J(T)/(dw)(18)
式中:d为超导带材的厚度,w为YBCO超导带材的宽度。
其中,超导带材的厚度d表示为:
d=(d1+d2+d3+d4)(19)
式中,d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,d2为银层的厚度,d3为YBCO层的厚度,d4为哈氏合金基底层的厚度。
YBCO层的第三电阻r3,根据欧姆定律可得:
式中:d3为YBCO层厚度,w为YBCO超导带材的宽度,Is为超导带材电流,le为超导带材的长度,T为超导带材的温度。
(4)第四电阻r4是哈氏合金基底层电阻,第四电阻r4是超导带材温度T的函数。根据欧姆定律可得:
式中,ρ4(T)为银的电阻率,w为YBCO超导带材的宽度,d4为哈氏合金基底层的厚度。
(5)按照YBCO超导带材等效电路的电阻并联结构,根据全电路欧姆定律,超导带材的电阻rsc为:
式中:r1、r2、r3、r4分别为第一电阻r1、第二电阻r2、第三电阻r3、第四电阻r4的阻值,Is为超导带材电流,T为超导带材的温度。
在电网稳态运行时,超导带材工作在超导态,即J<JC(T)时,超导带材的YBCO层的电阻为零,电流都通过YBCO层而导通,不会对电网造成电压降。当电网发生故障时,电网电流增大,超导带材失超而产生电阻,电流将在超导带材的各层之间分配。
步骤5.根据电阻型超导限流器的电路模型、超导带材的电阻rsc,计算线路电流I(t)和超导带材电流Is(t),输出线路电流I(t)和超导带材电流Is(t),并反馈超导带材的温度T给步骤3、反馈超导带材电流Is(t)给步骤4,实现***循环建模和仿真;
根据全电路欧姆定律,电阻型超导限流器的电流与电压关系为:
Us=IX+Ir+IRsc+IRload(23)
式中:Rload为负载电阻,Rsc为超导无感线圈电阻,I为线路电流,Us为交流电源电压,X为线路等效感抗,r为线路等效电阻。
其中,交流电源电压Us、线路等效感抗X分别表示为:
X=j2πfLs(25)
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,j为虚函数符号。
超导无感线圈的电阻为Rsc(T),超导无感线圈一般由m根超导带材并联组成,m≥1,根据分流定律:
Rsc(T)=rsc(T)/m(26)
I(t)=mIs(t)(27)
式中,rsc(T)为超导带材的电阻、I(t)和Is(t)分别为线路电流和超导带材的电流。
根据全电路欧姆定律,把式(24)-(27)代入式(23),得到电阻型超导限流器的电流与电压关系:
其中,相角θ为:
因此,通过超导带材的电流:
超导带材的电流峰值为:
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,m为超导带材的并联根数,r为线路等效电阻,Rload为负载电阻,rsc为超导带材的电阻。
在电网稳态时,电阻型超导限流器工作在超导态,超导无感线圈的电阻为Rsc(T)=0,根据式(30),线路电流为:
在电网发生短路故障时,负载电阻Rload减小为零,同时,电阻型超导限流器工作的正常态而产生电阻,因此,根据式(23),线路电流为:
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,m为超导带材的并联根数,r为线路等效电阻,Rload为负载电阻,rsc为超导带材的电阻。线路等效电感Ls和线路等效电阻r表示了电网短路故障的程度,而超导无感线圈的电阻rsc/m则体现了电阻型超导限流器的限流能力,当电阻rsc/m增大时,线路电流I(t)就减小。
本发明的主要优点:
1.本发明通过对YBCO超导带材结构和电路的简化和等效,使得超导带材在过流失超过程中的暂态电阻变化关系更加明确,便于计算。
2.本发明通过建立YBCO超导带材的热传导模型,准确分析了限流暂态过程中,超导带材温度变化的原因和热扩散的特征。
3.本发明给出了YBCO超导带材的电路和热传导建模所需的多种参数,使得超导带材建模更加简单。
4.本发明所建立的电阻型超导限流器的电路的数字模型,从超导带材的结构和热传导模型出发,解决了超导限流器在限流过程中的复杂的电磁场-热场耦合问题。
5.本发明所建立的电阻型超导限流器的电路的数字模型,使得建模方法简单化,为包含电阻型超导限流器的限流技术的建模和控制,提供了一种简单、易于操作的方法。
附图说明
图1为YBCO超导带材结构图;
图2为本发明的电阻型超导限流器的仿真计算流程图;
图3为本发明YBCO超导带材等效结构模型;
图4为本发明的电阻型超导限流器电路原理图;
图5为本发明的液氮传热系数与带材温度和运行温度的温差的关系图;
图6为本发明YBCO超导带材等效电路模型。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
本发明包括建立YBCO超导带材等效结构模型、等效电路模型和热传导模型、电阻型超导限流器的电路原理的数字建模仿真等步骤,如图2所示。
本发明的具体步骤如下:
步骤1.建立YBCO超导带材等效结构模型;
建立YBCO超导带材等效结构模型时,为便于分析YBCO超导带材各层在失超过程中的电阻情况,采用以下简化方法:
因YBCO超导带材的缓冲层的厚度很小,电阻很大,可以忽略。同时,又因YBCO超导带材的上表面铜层和下表面铜层的物理特性相同,可以合二为一进行分析。因此,所建立的YBCO超导带材等效结构模型简化分为4层:表面铜层、银层、YBCO层和哈氏合金基底层。表面铜层包括上表面铜层和下表面铜层,如图3所示。
步骤2.建立电阻型超导限流器的电路模型,给定电阻型超导限流器的电路参数:线路等效电感Ls、线路等效电阻r和负载电阻Rload,以及交流电源电压Us,在超导无感线圈电阻Rsc等于零的前提下,计算线路电流初始值I0(t)和超导带材电流初始值Is0(t);
电阻型超导限流器的电路原理图如图4所示,包括交流电源Us、线路等效感抗X、线路等效电阻r、断路器Br、超导无感线圈Rsc和负载电阻Rload。交流电源Us、线路等效感抗X、线路等效电阻r、断路器Br、超导无感线圈电阻Rsc和负载电阻Rload依次串联,交流电源Us和负载电阻Rload的一端接地;断路器Br和超导无感线圈电阻Rsc连接在第一连接点A上,超导无感线圈电阻Rsc和负载电阻Rload连接在第二连接点B上。
根据全电路欧姆定律,当超导无感线圈电阻Rsc为零时,电阻型超导限流器的电流与电压关系为:
Us=I0X+I0r+I0Rload(1)
式中:Rload为负载电阻,I0为线路电流,Us为交流电源电压,X为线路等效感抗,r为线路等效电阻。
其中,交流电源电压Us、线路等效感抗X分别表示为:
X=j2πfLs(3)
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,j为虚函数符号。
线路电流的初始值为:
超导无感线圈的电阻为Rsc(T),超导无感线圈一般由m(m≥1)根超导带材并联组成,根据分流定律:
Is0(t)=I0(t)/m(5)
式中,I0(t)和Is0(t)分别为线路电流初始值和超导带材的电流初始值。
步骤3.建立YBCO超导带材的热传导模型,给定超导带材的结构参数:超导带材的宽度w、厚度d和长度le,所述的后度d包括了包括表面铜层厚度d1、银层的厚度d2、YBCO层的厚度d3和哈氏合金基底层厚度d4;给定超导无感线圈的初始运行条件:工作温度Top,计算超导带材的温度T;
(1)YBCO超导带材直接在液氮浸泡环境中冷却,根据热平衡方程,沿超导带材的长度方向,一维热传导方程为:
式中,Kcm(T)为热传导系数,Vcm和Ccm分别为YBCO超导带材的密度与比热容,gj(T)和Wcool(T)分别是超导带材的焦耳热和散失热量。
YBCO超导带材的密度Vcm与比热容Ccm满足式:
其中,ν1,ν2,ν3,ν4分别为表面铜层、银层、YBCO层和哈氏合金基底层的密度。d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,d2为银层的厚度,d3为YBCO层的厚度,d4为哈氏合金基底层的厚度,d为YBCO超导带材的厚度。C1,C2,C3,C4分别为表面铜层、银层、YBCO层和哈氏合金基底层的比热容。各种密度和比热容均可查阅手册获得。
其中,K1,K2,K3,K4分别为表面铜层,银层,YBCO层和哈氏合金基底层的热导率。各种材料的热导率均可查阅手册获得。
(2)根据焦耳定律,超导带材的焦耳热
其中,US(t)为超导带材的电压,IS(t)为超导带材的电流。
(3)超导带材浸泡在液氮中,根据经验公式,散失热量完全由液氮带走,超导带材的散失热量Wcool(T)为:
Wcool(T)=hA(T)(10)
其中,A为超导带材与液氮的接触面积,即超导带材的表面积;h为液氮传热系数,与超导带材和液氮的温度差ΔT(T-Top)有关,根据实验可知,液氮传热系数h的取值,如图5所示。对应于不同的温差ΔT(T-Top),热传递过程有对流、核沸腾、过渡态和膜沸腾4种状态,不同的状态对应不同的液氮传热系数h。液氮传热系数h的与拟合结果:
步骤4.建立YBCO超导带材等效电路模型,依据超导带材温度T和超导带材电流Is(t),计算YBCO超导带材的电阻rsc;
YBCO超导带材的等效电路为4个电阻并联结构。第一电阻r1为表面铜层电阻,第二电阻r2为银层的电阻,第三电阻r3为YBCO层的电阻,第四电阻r4为哈氏合金基底层的电阻,如图6所示。
根据所建立的YBCO超导带材等效电路模型,并根据超导带材温度和电流,根据欧姆定律和电路原理,计算YBCO超导带材的等效电路的电阻:
(1)第一电阻r1为表面铜层电阻,为铜材料制作。第一电阻r1是超导带材温度T的函数。根据欧姆定律可得:
其中,ρ1(T)为铜的电阻率,是超导带材温度T的函数,w为YBCO超导带材的宽度,d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,le为超导带材的长度。
(2)第二电阻r2为银层电阻,为银材料制作。第二电阻r2是超导带材温度T的函数。根据欧姆定律可得:
其中,ρ2(T)为银的电阻率,w为YBCO超导带材的宽度,d2为银层的厚度,le为超导带材的长度。
(3)第三电阻r3为YBCO层电阻。根据超导带材的电阻变化规律,第三电阻r3的电阻率ρ3(T)是超导带材温度T的函数,通过对YBCO超导带材的特性测试,并采用数值拟合的方法可得第三电阻r3的电阻率ρ3(T):
ρ3(T)=ρ31(T)+ρ32(T)(14)
其中,ρ31(T)和ρ32(T)由分段函数拟合而成:
其中,Jc(T)是超导带材临界电流密度Jc随温度T的函数:
Jc(T)=Jc0[(Tc-T)/(Tc-Top)]1.5(17)
其中,Jc为超导带材临界电流密度。Jc0=2×106A/cm2,为77K下的临界电流密度;Tc=92K,为YBCO的临界温度,Top为工作温度,在液氮池中为77K。参数n1=3;n2=20;γ=2,均为根据超导带材的特性而得到的拟合参数。针对ρ31(T)和ρ32(T)为零的情况,在实际计算过程中,可以取为一个非常小的数据,如10-19Ω··cm等,避免运算出错。
根据欧姆定律,超导带材的电流Is表示为:
IS(T)=J(T)/(dw)(18)
式中:d为超导带材的厚度,w为YBCO超导带材的宽度。
其中,超导带材的厚度d表示为:
d=(d1+d2+d3+d4)(19)
式中,d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,d2为银层的厚度,d3为YBCO层的厚度,d4为哈氏合金基底层的厚度。
YBCO层的第三电阻r3,根据欧姆定律可得:
式中:d3为YBCO层厚度,w为YBCO超导带材的宽度,Is为超导带材电流,le为超导带材的长度,T为超导带材的温度。
(4)第四电阻r4是哈氏合金基底层电阻,第四电阻r4是超导带材温度T的函数。根据欧姆定律可得:
式中,ρ4(T)为银的电阻率w为YBCO超导带材的宽度,d4为哈氏合金基底层的厚度。
(5)按照YBCO超导带材等效电路的电阻并联结构,根据全电路欧姆定律,超导带材的电阻rsc为:
式中:r1、r2、r3、r4分别为第一电阻r1、第二电阻r2、第三电阻r3、第四电阻r4的阻值,Is为超导带材电流,T为超导带材的温度。
在电网稳态运行时,超导带材工作在超导态,即J<JC(T)时,超导带材的YBCO层的电阻为零,电流都通过YBCO层而导通,不会对电网造成电压降。当电网发生故障时,电网电流增大,超导带材失超而产生电阻,电流将在超导带材的各层之间分配。
步骤5.根据电阻型超导限流器的电路模型、超导带材的电阻rsc,计算线路电流I(t)和超导带材电流Is(t);输出线路电流I(t)和超导带材电流Is(t),并反馈超导带材的温度T给步骤3、反馈超导带材电流Is(t)给步骤4,实现***循环建模和仿真。
根据全电路欧姆定律,电阻型超导限流器的电流与电压关系为:
Us=IX+Ir+IRsc+IRload(23)
式中:Rload为负载电阻,Rsc为超导无感线圈电阻,I为线路电流,Us为交流电源电压,X为线路等效感抗,r为线路等效电阻。
其中,交流电源电压Us、线路等效感抗X分别表示为:
X=j2πfLs(25)
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,j为虚函数符号。
超导无感线圈的电阻为Rsc(T),超导无感线圈一般由m根超导带材并联组成,m≥1,根据分流定律:
Rsc(T)=rsc(T)/m(26)
I(t)=mIs(t)(27)
式中,rsc(T)为超导带材的电阻、I(t)和Is(t)分别为线路电流和超导带材的电流。
根据全电路欧姆定律,把式(24)-(27)代入式(23),得到电阻型超导限流器的电流与电压关系:
其中,相角θ为:
因此,通过超导带材的电流:
超导带材的电流峰值为:
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,m为超导带材的并联根数,r为线路等效电阻,Rload为负载电阻,rsc为超导带材的电阻。
在电网稳态时,电阻型超导限流器工作在超导态,超导无感线圈的电阻为Rsc(T)=0,根据式(30),线路电流为:
在电网发生短路故障时,负载电阻Rload减小为零,同时,电阻型超导限流器工作的正常态而产生电阻,因此,根据式(23),线路电流为:
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,m为超导带材的并联根数,r为线路等效电阻,Rload为负载电阻,rsc为超导带材的电阻。线路等效电感Ls和线路等效电阻r表示了电网短路故障的程度,而超导无感线圈的电阻rsc/m则体现了电阻型超导限流器的限流能力,当电阻rsc/m增大时,线路电流I(t)减小。
本发明基于YBCO超导带材的电阻型超导限流器数字建模仿真方法所构造的超导带材结构模型和超导带材等效电路模型,简化了超导带材的结构,全面描述了超导带材在由超导态转化为正常态过程的电阻变化规律;通过YBCO超导带材的热传导建模的建立,分析了超导带材各个组分在不同温度下的比热容、电导率等的变化规律,以及带材的散热情况;通过电阻型超导限流器的电路模型的建立,提供了简单而有效的电路和限流器模型和分析方法。通过建模提供了一种有效的基于YBCO超导带材的电阻型超导限流器的故障分析方法,提高了工作效率。通过本发明所提供的建模方法,不仅解决了基于YBCO超导带材的电阻型超导限流器精确分析方法,而且,也可用于多单元和复杂结构的超导限流器的研究,推动了超导限流器在电网中的应用。
Claims (2)
1.一种基于YBCO超导带材的电阻型超导限流器数字建模仿真方法,其特征在于:所述建模仿真方法在建立YBCO超导带材等效结构模型、YBCO超导带材等效电路模型、YBCO超导带材热传导模型和电阻型超导限流器电路模型的基础上,根据给定电阻型超导限流器的电路参数、超导带材的结构参数和初始运行条件,在超导无感线圈电阻Rsc为零的前提下,计算线路电流的初始值和超导带材电流的初始值;根据YBCO超导带材的热传导模型计算超导带材的温度;根据YBCO超导带材等效电路模型、超导带材温度和电流,计算超导带材的电阻;根据电阻型超导限流器的电路模型,计算并反馈线路电流和超导带材电流,实现电阻型超导限流器的建模和仿真。
2.按照权利要求1所述的基于YBCO超导带材的电阻型超导限流器数字建模仿真方法,其特征在于:所述的建模仿真方法的具体步骤如下:
步骤1.建立YBCO超导带材等效结构模型;
建立YBCO超导带材等效结构模型时,为便于分析YBCO超导带材各层在失超过程中的电阻情况,忽略YBCO超导带材缓冲层的厚度,将YBCO超导带材的上表面铜层和下表面铜层合二为一进行分析;因此,所建立的YBCO超导带材等效结构模型简化分为4层:表面铜层、银层、YBCO层和哈氏合金基底层;表面铜层包括上表面铜层和下表面铜层;
步骤2.建立电阻型超导限流器的电路模型,给定电阻型超导限流器的电路参数:线路等效电感Ls、线路等效电阻r和负载电阻Rload,以及交流电源电压Us,在超导无感线圈电阻Rsc等于零的前提下,计算线路电流初始值I0(t)和超导带材电流初始值Is0(t);
电阻型超导限流器包括交流电源U s 、线路等效感抗X、线路等效电阻r、断路器Br、超导无感线圈Rsc和负载电阻Rload;交流电源Us、线路等效感抗X、线路等效电阻r、断路器Br、超导无感线圈电阻Rsc和负载电阻Rload依次串联,交流电源Us和负载电阻Rload的一端接地;断路器Br和超导无感线圈电阻Rsc连接在第一连接点A上,超导无感线圈电阻Rsc和负载电阻Rload连接在第二连接点B上;
根据全电路欧姆定律,当超导无感线圈电阻Rsc为零时,电阻型超导限流器的电流与电压关系为:
Us=I0X+I0r+I0Rload(1)
式中:Rload为负载电阻,I0为线路电流,Us为交流电源电压,X为线路等效感抗,r为线路等效电阻;
其中,交流电源电压Us、线路等效感抗X分别表示为:
X=j2πfLs(3)
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,j为虚函数符号;
线路电流初始值为:
超导无感线圈的电阻为Rsc(T),超导无感线圈一般由m根超导带材并联组成,m≥1,根据分流定律:
Is0(t)=I0(t)/m(5)
式中,I0(t)为线路电流初始值,Is0(t)为超导带材的电流初始值;
步骤3.建立YBCO超导带材的热传导模型,给定超导带材的结构参数:超导带材的宽度w、包括表面铜层厚度d1,银层厚度d2,YBCO层厚度d3,哈氏合金基底层厚度d4的厚度d,以及长度le,给定超导无感线圈的初始运行条件:工作温度Top,计算超导带材的温度T;
(1)YBCO超导带材直接在液氮浸泡环境中冷却,根据热平衡方程,沿超导带材的长度方向,一维热传导方程为:
式中,Kcm(T)为热传导系数,Vcm和Ccm分别为YBCO超导带材的密度与比热容,gj(T)和Wcool(T)分别是超导带材的焦耳热和散失热量;
YBCO超导带材的密度Vcm与比热容Ccm满足式:
其中,ν1,ν2,ν3,ν4分别为表面铜层、银层、YBCO层和哈氏合金基底层的密度,d1为表面铜层厚度,为上表面铜层和下表面铜层的厚度之和,d2为银层的厚度,d3为YBCO层的厚度,d4为哈氏合金基底层的厚度,d为YBCO超导带材的厚度;C1,C2,C3,C4分别为表面铜层、银层、YBCO层和哈氏合金基底层的比热容;
其中,K1,K2,K3,K4分别为表面铜层,银层,YBCO层和哈氏合金基底层的热导率;
(2)根据焦耳定律,超导带材的焦耳热:
其中,US(t)为超导带材的电压,IS(t)为超导带材的电流;
(3)超导带材浸泡在液氮中,根据经验公式,散失热量完全由液氮带走,超导带材的散失热量Wcool(T)为:
Wcool(T)=hA(T)(10)
其中,A为超导带材与液氮的接触面积,即超导带材的表面积;h为液氮传热系数,与超导带材和液氮的温度差ΔT(T-Top)有关,根据实验可知,液氮传热系数h的取值;对应于不同的温差ΔT(T-Top),热传递过程有对流、核沸腾、过渡态和膜沸腾4种状态,不同的状态对应不同的液氮传热系数h;液氮传热系数h的与拟合结果:
步骤4.建立YBCO超导带材等效电路模型,依据超导带材温度T和超导带材电流Is(t),计算YBCO超导带材的电阻rsc;
YBCO超导带材的等效电路为4个电阻并联结构;第一电阻r1为表面铜层电阻,第二电阻r2为银层的电阻,第三电阻r3为YBCO层的电阻,第四电阻r4为哈氏合金基底层的电阻;
根据所建立的YBCO超导带材等效电路模型,并根据超导带材温度和电流,根据欧姆定律和电路原理,计算YBCO超导带材的等效电路的电阻:
(1)第一电阻r1为表面铜层电阻,为铜材料制作;第一电阻r1是超导带材温度T的函数;根据欧姆定律可得:
其中,ρ1(T)为铜的电阻率,是超导带材温度T的函数,w为YBCO超导带材的宽度,d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,le为超导带材的长度;
(2)第二电阻r2为银层电阻,为银材料制作;第二电阻r2是超导带材温度T的函数;根据欧姆定律可得:
其中,ρ2(T)为银的电阻率,w为YBCO超导带材的宽度,d2为银层的厚度,le为超导带材的长度;
(3)第三电阻r3为YBCO层电阻;根据超导带材的电阻变化规律,第三电阻r3的电阻率ρ3(T)是超导带材温度T的函数,通过对YBCO超导带材的特性测试,并采用数值拟合的方法可得第三电阻r3的电阻率ρ3(T):
ρ3(T)=ρ31(T)+ρ32(T)(14)
其中,ρ31(T)和ρ32(T)由分段函数拟合而成:
其中,Jc(T)是超导带材临界电流密度Jc随温度T的函数:
Jc(T)=Jc0[(Tc-T)/(Tc-Top)]1.5(17)
其中,Jc为超导带材临界电流密度;Jc0=2×106A/cm2,为77K下的临界电流密度;Tc=92K,为YBCO的临界温度,Top为工作温度,在液氮池中为77K;参数n1=3;n2=20;γ=2,均为根据超导带材的特性而得到的拟合参数;根据欧姆定律,超导带材的电流Is表示为:
IS(T)=J(T)/(dw)(18)
式中:d为超导带材的厚度,w为YBCO超导带材的宽度;
其中,超导带材的厚度d表示为:
d=(d1+d2+d3+d4)(19)
式中,d1为表面铜层厚度,即上表面铜层和下表面铜层的厚度之和,d2为银层的厚度,d3为YBCO层的厚度,d4为哈氏合金基底层的厚度;
YBCO层的第三电阻r3,根据欧姆定律可得:
式中:d3为YBCO层厚度,w为YBCO超导带材的宽度,Is为超导带材电流,le为超导带材的长度,T为超导带材的温度;
(4)第四电阻r4是哈氏合金基底层电阻,第四电阻r4是超导带材温度T的函数,根据欧姆定律可得:
式中,ρ4(T)为银的电阻率w为YBCO超导带材的宽度,d4为哈氏合金基底层的厚度;
(5)按照YBCO超导带材等效电路的电阻并联结构,根据全电路欧姆定律,超导带材的电阻rsc为:
式中:r1、r2、r3、r4分别为第一电阻r1、第二电阻r2、第三电阻r3、第四电阻r4的阻值,Is为超导带材电流,T为超导带材的温度;
在电网稳态运行时,超导带材工作在超导态,即J<JC(T)时,超导带材的YBCO层的电阻为零,电流都通过YBCO层而导通,不会对电网造成电压降;当电网发生故障时,电网电流增大,超导带材失超而产生电阻,电流将在超导带材的各层之间分配;
步骤5.根据电阻型超导限流器的电路模型、超导带材的电阻rsc,计算线路电流I(t)和超导带材电流Is(t);输出线路电流I(t)和超导带材电流Is(t),并反馈超导带材的温度T给步骤3,反馈超导带材电流Is(t)给步骤4,实现***循环建模和仿真;
根据全电路欧姆定律,电阻型超导限流器的电流与电压关系为:
Us=IX+Ir+IRsc+IRload(23)
式中:Rload为负载电阻,Rsc为超导无感线圈电阻,I为线路电流,Us为交流电源电压,X为线路等效感抗,r为线路等效电阻;
其中,交流电源电压Us、线路等效感抗X分别表示为:
X=j2πfLs(25)
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,j为虚函数符号;
超导无感线圈的电阻为Rsc(T),超导无感线圈一般由m根超导带材并联组成,m≥1,根据分流定律:
Rsc(T)=rsc(T)/m(26)
I(t)=mIs(t)(27)
式中,rsc(T)为超导带材的电阻、I(t)和Is(t)分别为线路电流和超导带材的电流;
根据全电路欧姆定律,把式(24)-(27)代入式(23),得到电阻型超导限流器的电流与电压关系:
其中,相角θ为:
因此,通过超导带材的电流:
超导带材的电流峰值为:
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,m为超导带材的并联根数,r为线路等效电阻,Rload为负载电阻,rsc为超导带材的电阻;
在电网稳态时,电阻型超导限流器工作在超导态,超导无感线圈的电阻为Rsc(T)=0,根据式(30),线路电流为:
在电网发生短路故障时,负载电阻Rload减小为零,同时,电阻型超导限流器工作的正常态而产生电阻,因此,根据式(23),线路电流为:
式中,U0为交流电源电压有效值,f为交流电源频率,Ls为线路等效电感,m为超导带材的并联根数,r为线路等效电阻,Rload为负载电阻,rsc为超导带材的电阻;线路等效电感Ls和线路等效电阻r表示了电网短路故障的程度,而超导无感线圈的电阻rsc/m则体现了电阻型超导限流器的限流能力,当电阻rsc/m增大时,线路电流I(t)减小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510204556.7A CN105160047B (zh) | 2015-04-27 | 2015-04-27 | 基于ybco超导带材的电阻型超导限流器数字建模仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510204556.7A CN105160047B (zh) | 2015-04-27 | 2015-04-27 | 基于ybco超导带材的电阻型超导限流器数字建模仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105160047A true CN105160047A (zh) | 2015-12-16 |
CN105160047B CN105160047B (zh) | 2018-04-03 |
Family
ID=54800903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510204556.7A Expired - Fee Related CN105160047B (zh) | 2015-04-27 | 2015-04-27 | 基于ybco超导带材的电阻型超导限流器数字建模仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105160047B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107478938A (zh) * | 2017-09-28 | 2017-12-15 | 南方电网科学研究院有限责任公司 | 电阻型直流限流器、电阻型直流限流器的测试***及方法 |
CN108763712A (zh) * | 2018-05-22 | 2018-11-06 | 深圳供电局有限公司 | 一种高温超导corc仿真建模方法 |
CN108923401A (zh) * | 2018-08-27 | 2018-11-30 | 广东电网有限责任公司 | 一种超导限流器故障电阻和限流电流的计算方法及装置 |
CN109712755A (zh) * | 2018-12-29 | 2019-05-03 | 上海超导科技股份有限公司 | 适用于电阻型高温超导限流器的带材 |
CN110161329A (zh) * | 2019-04-19 | 2019-08-23 | 华中科技大学 | 一种二代高温超导带材的失超电阻仿真预测方法和*** |
CN110441602A (zh) * | 2019-08-05 | 2019-11-12 | 南京尚众微电子科技有限公司 | 用于表征铁基超导体电子向列相的复电阻抗测量方法 |
CN111244920A (zh) * | 2020-03-09 | 2020-06-05 | 广东电网有限责任公司电力科学研究院 | 一种高压大容量电阻型超导限流器仿真建模方法及设备 |
CN112018741A (zh) * | 2020-09-10 | 2020-12-01 | 广东电网有限责任公司电力科学研究院 | 基于重合闸的超导交流限流器的设计方法和装置 |
CN112906268A (zh) * | 2021-02-07 | 2021-06-04 | 中国电力科学研究院有限公司 | 一种ybco高温超导单元失超电阻率的计算方法及*** |
CN113064108A (zh) * | 2021-03-15 | 2021-07-02 | 西安交通大学 | 一种测量超导带在强制对流冷却下失超及恢复特性的装置 |
CN113076642A (zh) * | 2021-03-31 | 2021-07-06 | 上海超导科技股份有限公司 | 根据线圈设计定制化超导带材的方法、超导带材及线圈 |
CN114839488A (zh) * | 2022-05-07 | 2022-08-02 | 上海超导科技股份有限公司 | 超导带材耐过流冲击判定***和方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1619913A (zh) * | 2003-09-30 | 2005-05-25 | Lg产电株式会社 | 电阻型超导故障限流器 |
KR20100042782A (ko) * | 2008-10-17 | 2010-04-27 | 한국산업기술대학교산학협력단 | 안정화층이 형성되지 않은 초전도 선재의 임계전류 측정장치 |
CN101741076A (zh) * | 2010-01-15 | 2010-06-16 | 中国科学院电工研究所 | 一种基于ybco高温超导带材的短路故障限流器 |
-
2015
- 2015-04-27 CN CN201510204556.7A patent/CN105160047B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1619913A (zh) * | 2003-09-30 | 2005-05-25 | Lg产电株式会社 | 电阻型超导故障限流器 |
KR20100042782A (ko) * | 2008-10-17 | 2010-04-27 | 한국산업기술대학교산학협력단 | 안정화층이 형성되지 않은 초전도 선재의 임계전류 측정장치 |
CN101741076A (zh) * | 2010-01-15 | 2010-06-16 | 中国科学院电工研究所 | 一种基于ybco高温超导带材的短路故障限流器 |
Non-Patent Citations (2)
Title |
---|
秦玥 等: "电阻型超导限流器对电力***暂态稳定的影响分析", 《华东电力》 * |
陈妍君 等: "电阻型超导限流器仿真模型及其对10kV配电网的影响", 《电力自动化设备》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107478938A (zh) * | 2017-09-28 | 2017-12-15 | 南方电网科学研究院有限责任公司 | 电阻型直流限流器、电阻型直流限流器的测试***及方法 |
CN107478938B (zh) * | 2017-09-28 | 2023-04-28 | 南方电网科学研究院有限责任公司 | 电阻型直流限流器、电阻型直流限流器的测试***及方法 |
CN108763712A (zh) * | 2018-05-22 | 2018-11-06 | 深圳供电局有限公司 | 一种高温超导corc仿真建模方法 |
US11190010B2 (en) | 2018-08-27 | 2021-11-30 | Guangdong Power Grid Co. Ltd. | Annular bearer network and service bearing implementation method therefor |
CN108923401A (zh) * | 2018-08-27 | 2018-11-30 | 广东电网有限责任公司 | 一种超导限流器故障电阻和限流电流的计算方法及装置 |
CN109712755A (zh) * | 2018-12-29 | 2019-05-03 | 上海超导科技股份有限公司 | 适用于电阻型高温超导限流器的带材 |
CN110161329A (zh) * | 2019-04-19 | 2019-08-23 | 华中科技大学 | 一种二代高温超导带材的失超电阻仿真预测方法和*** |
CN110161329B (zh) * | 2019-04-19 | 2020-08-04 | 华中科技大学 | 一种二代高温超导带材的失超电阻仿真预测方法和*** |
CN110441602A (zh) * | 2019-08-05 | 2019-11-12 | 南京尚众微电子科技有限公司 | 用于表征铁基超导体电子向列相的复电阻抗测量方法 |
CN110441602B (zh) * | 2019-08-05 | 2021-06-18 | 南京尚众微电子科技有限公司 | 用于表征铁基超导体电子向列相的复电阻抗测量方法 |
CN111244920A (zh) * | 2020-03-09 | 2020-06-05 | 广东电网有限责任公司电力科学研究院 | 一种高压大容量电阻型超导限流器仿真建模方法及设备 |
CN112018741A (zh) * | 2020-09-10 | 2020-12-01 | 广东电网有限责任公司电力科学研究院 | 基于重合闸的超导交流限流器的设计方法和装置 |
CN112906268A (zh) * | 2021-02-07 | 2021-06-04 | 中国电力科学研究院有限公司 | 一种ybco高温超导单元失超电阻率的计算方法及*** |
CN112906268B (zh) * | 2021-02-07 | 2024-03-29 | 中国电力科学研究院有限公司 | 一种ybco高温超导单元失超电阻率的计算方法及*** |
CN113064108A (zh) * | 2021-03-15 | 2021-07-02 | 西安交通大学 | 一种测量超导带在强制对流冷却下失超及恢复特性的装置 |
CN113076642A (zh) * | 2021-03-31 | 2021-07-06 | 上海超导科技股份有限公司 | 根据线圈设计定制化超导带材的方法、超导带材及线圈 |
CN114839488A (zh) * | 2022-05-07 | 2022-08-02 | 上海超导科技股份有限公司 | 超导带材耐过流冲击判定***和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105160047B (zh) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105160047A (zh) | 基于ybco超导带材的电阻型超导限流器数字建模仿真方法 | |
Ye et al. | Modeling and simulation of high temperature resistive superconducting fault current limiters | |
CN106897487B (zh) | 一种高温超导磁体非线性分析的建模方法 | |
CN104834781A (zh) | 一种基于平波电抗器多次换向失败时的暂态温度场仿真方法 | |
CN108763712A (zh) | 一种高温超导corc仿真建模方法 | |
De Sousa et al. | Simulation of a superconducting fault current limiter: A case study in the Brazilian power system with possible recovery under load | |
Zermeno et al. | Simulation of an HTS synchronous superconducting generator | |
CN103840442B (zh) | 一种短路故障限流器 | |
CN103616578A (zh) | 电力***中高温超导电缆***的多状态可靠性确定方法 | |
Liang et al. | Modeling and experiment of the current limiting performance of a resistive superconducting fault current limiter in the experimental system | |
CN103986151B (zh) | 一种交流受端电网直流受电规模的计算方法 | |
Xiang et al. | DC fault study of a point-to-point HVDC system integrating offshore wind farm using high-temperature superconductor DC cables | |
CN105245188A (zh) | 光伏逆变器能耗特性在线预测方法及其装置 | |
Nemdili et al. | Electrothermal modeling of coated conductor for a resistive superconducting fault-current limiter | |
CN101446598A (zh) | 一种变截面电流引线 | |
WO2024139415A1 (zh) | 一种应用于超导电缆数字孪生***的仿真方法及平台 | |
CN101446611A (zh) | 一种使用变截面电流引线测量高温超导体临界电流特性的装置 | |
Park et al. | Experimental and numerical analysis of high resistive coated conductor for conceptual design of fault current limiter | |
Yasui et al. | Temperature and pressure distribution simulations of 3-km-long high-temperature superconducting power cable system with fault current for 66-kV-class transmission lines | |
Sun et al. | Accuracy analysis of static voltage stability indices based on power flow model | |
Chen et al. | Numerical study of normal‐zone evolution and stability of composite superconductors | |
Sharifi et al. | Electro-thermal modeling and optimization algorithm of resistive superconducting fault current limiters | |
CN108388744A (zh) | 一种稳态变压器热效应的仿真方法及*** | |
Gardiner et al. | Operation of power electronic converters at cryogenic temperatures for utility energy conditioning applications | |
CN209015528U (zh) | 一种并联型线缆模拟装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 |
Granted publication date: 20180403 Termination date: 20190427 |
|
CF01 | Termination of patent right due to non-payment of annual fee |