CN106324490A - 变压器有载分接开关机械故障诊断方法 - Google Patents

变压器有载分接开关机械故障诊断方法 Download PDF

Info

Publication number
CN106324490A
CN106324490A CN201610632448.4A CN201610632448A CN106324490A CN 106324490 A CN106324490 A CN 106324490A CN 201610632448 A CN201610632448 A CN 201610632448A CN 106324490 A CN106324490 A CN 106324490A
Authority
CN
China
Prior art keywords
characteristic vector
load tap
changer
vibration
load
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.)
Pending
Application number
CN201610632448.4A
Other languages
English (en)
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.)
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610632448.4A priority Critical patent/CN106324490A/zh
Publication of CN106324490A publication Critical patent/CN106324490A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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/327Testing of circuit interrupters, switches or circuit-breakers
    • G01R31/3271Testing of circuit interrupters, switches or circuit-breakers of high voltage or medium voltage devices
    • G01R31/3275Fault detection or status indication
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/025Measuring arrangements

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Protection Of Transformers (AREA)

Abstract

本发明涉及一种变压器有载分接开关机械故障诊断方法,其技术特点是包括:将振动检测探头贴敷于有载分接开关的箱壁上,捕获分接开关操作过程中产生的振动信号;通过小波包分析技术对振动信号进行能量频带分解,提取信号的频谱特征,从而构成状态特征向量,然后通过信号采集和存储形成特征向量库;利用欧几里德距离函数分析特征向量的变化趋势,对有载分接开关的机械状态进行诊断和评估。本发明能够实时捕获分接开关操作过程中产生的振动信号,提取振动信号中的特征量,并与历时数据进行比对,判断是否出现故障以及故障发展趋势,节约了检修成本,提高了检测效率,可及时有效地发现早期机械故障,避免故障严重化。

Description

变压器有载分接开关机械故障诊断方法
技术领域
本发明属于变压器高压试验技术领域,尤其是一种变压器有载分接开关机械故障诊断方法。
背景技术
有载分接开关是一种为变压器在负载变化时提供恒定电压的开关设备,是变压器中唯一频繁操作的机械与电气结合一体的设备。随着调压次数的增多,有载分接开关的机械故障和电气故障也相应增加,很容易出现烧毁造成变压器停电事故。据有关统计,有载分接开关的机械故障占总故障数量的90%以上。目前,对有载分接开关机械故障的诊断方法主要是停电检修。当有载分接开关运行一定年限或调压次数较多时,需要将有载分接开关的切换芯子从变压器中取出进行检修维护,检查切换机构和组部件的机械状态是否良好。有载分接开关的停电检修往往周期较长,难以及时发现早期的机械故障,经常在停电检修前即发生故障损坏,并且停电检修影响变压器正常运行,需要耗费大量人力、物力和财力。
发明内容
本发明的目的在于克服现有技术的不足,提供一种变压器有载分接开关机械故障诊断方法,解决停电检修存在的检测周期长、浪费人力财力等问题。
本发明解决其技术问题是采取以下技术方案实现的:
一种变压器有载分接开关机械故障诊断方法,包括以下步骤:
步骤1、将振动检测探头贴敷于有载分接开关的箱壁上,捕获分接开关操作过程中产生的振动信号;
步骤2、通过小波包分析技术对振动信号进行能量频带分解,提取信号的频谱特征,从而构成状态特征向量,然后通过信号采集和存储形成特征向量库;
步骤3、利用欧几里德距离函数分析特征向量的变化趋势,对有载分接开关的机械状态进行诊断和评估。
所述步骤2的具体实现方法为:
对振动信号进行j层小波包分解,得到2j个频率段,定义在某一特定频率段的小波能量为该段小波系数的平方和Em
E m = Σ k | C m ( k ) | 2
得到总能量Etot的表达式为:
E t o t = Σ m Σ k | C m ( k ) | 2 = Σ m E m
定义某频段的相对小波包能量pm
p m = E m E t o t
定义相对小波包能量特征向量Ki
K i = ( p 0 , p 1 , ... , p 2 j - 1 )
上式中,m=0,1,……,2j-1,
最后根据特征向量建立特征向量库。
所述步骤3的具体实现方法为:
采用待监测的特征向量KT与标准特征向量库中特征向量Ki的距离d(KT,Ki)作为模式判别的依据,如果:
d(KT,Ki)=min{d(KT,Ki),...,d(KT,Kn)}
则将待监测向量的状态归入正常状态,采用如下欧几里德距离函数计算距离:
d ( K T , K i ) = Σ k = 0 2 j - 1 ( p T , k - p i , k )
式中,j为小波包分解的层数;
在状态检测过程中,取如下两组数据:一组是待检测时刻的变压器有载分接开关振动数据,用来建立测试数据库,另一组是变压器有载分接开关正常工作时的振动数据,用来建立标准数据库;当欧几里德距离小于某一阈值时,表明测试数据库中特征向量和标准数据库中的特征向量表征同一状态,为正常工作状态,当距离大于某一阈值时,表明该时刻有载分接开关可能出现机械故障。
本发明的优点和积极效果是:
本发明将振动检测探头贴敷于有载分接开关的箱壁上,捕获分接开关操作过程中产生的振动信号,通过数据处理,提取振动信号中的特征量,并与历时数据进行比对,判断是否出现故障以及故障发展趋势,本发明不需要变压器停电,也不需要将切换芯子从变压器中取出,节约了检修成本,提高了检测效率,可以及时有效的发现早期机械故障,避免故障严重化。
具体实施方式
以下结合一个实施例对本发明做进一步详述:
一种变压器有载分接开关机械故障诊断方法的设计原理为:由于在有载分接开关出现机械故障后,其分接转换过程的振动特性将会有明显的改变,主要表现在不同频率段的振动信号具有不同的衰减或增强现象,因此,首先通过小波包分析技术对振动信号进行能量频带分解,提取信号的频谱特征,从而构成状态特征向量,然后通过信号采集和存储形成特征向量库,分析特征向量的变化趋势即可对有载分接开关的机械状态进行诊断和评估。
本发明变压器有载分接开关机械故障诊断方法包括以下步骤:
步骤1、将振动检测探头贴敷于有载分接开关的箱壁上,捕获分接开关操作过程中产生的振动信号。
步骤2、通过小波包分析技术对振动信号进行能量频带分解,提取信号的频谱特征,从而构成状态特征向量,然后通过信号采集和存储形成特征向量库。
在本步骤中,通过对振动信号进行j层小波包分解,可以得到2j个频率段,定义在某一特定频率段的小波能量为该段小波系数的平方和:
E m = Σ k | C m ( k ) | 2 , ( m = 0 , 1 , ... , 2 j - 1 ) - - - ( 1 )
由此可得到总能量的表达式为:
E t o t = Σ m Σ k | C m ( k ) | 2 = Σ m E m , ( m = 0 , 1 , ... , 2 j - 1 ) - - - ( 2 )
定义某频段的相对小波包能量:
p m = E m E t o t , ( m = 0 , 1 , ... , 2 j - 1 ) - - - ( 3 )
定义相对小波包能量特征向量:
K i = ( p 0 , p 1 , ... , p 2 j - 1 ) - - - ( 4 )
根据特征数据得到表征其状态的特征向量之后,就可以建立特征向量库了。特征向量库具有开放性,而且随着检测数据的增加而增加。
步骤3、利用欧几里德距离函数分析特征向量的变化趋势,即可对有载分接开关的机械状态进行诊断和评估。
在特征向量库建立以后,就可以进行模式识别的工作了。采用待监测的特征向量KT与标准特征向量库中特征向量Ki的距离d(KT,Ki)作为模式判别的依据。若
d(KT,Ki)=min{d(KT,Ki),...,d(KT,Kn)}(5)
也就是说KT与Ki最接近,则把待监测向量的状态归入正常状态。距离的计算采用最常见的欧几里德距离函数,计算公式为:
d ( K T , K i ) = Σ k = 0 2 j - 1 ( p T , k - p i , k ) - - - ( 6 )
式中,j为小波包分解的层数。在状态检测中,一般取两组数据,一组是待检测时刻的变压器有载分接开关振动数据,用来建立测试数据库,一组是变压器有载分接开关正常工作时的振动数据,用来建立标准数据库。
当距离矩阵中的欧几里德距离小于某一阈值时,表明测试数据库中特征向量和标准数据库中的特征向量表征同一状态,为正常工作状态,因此欧几里德距离越小越好;当距离大于某一阈值时,表明该时刻有载分接开关可能出现机械故障,可作为故障判断的依据。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。

Claims (3)

1.一种变压器有载分接开关机械故障诊断方法,其特征在于包括以下步骤:
步骤1、将振动检测探头贴敷于有载分接开关的箱壁上,捕获分接开关操作过程中产生的振动信号;
步骤2、通过小波包分析技术对振动信号进行能量频带分解,提取信号的频谱特征,从而构成状态特征向量,然后通过信号采集和存储形成特征向量库;
步骤3、利用欧几里德距离函数分析特征向量的变化趋势,对有载分接开关的机械状态进行诊断和评估。
2.根据权利要求1所述的变压器有载分接开关机械故障诊断方法,其特征在于:所述步骤2的具体实现方法为:
对振动信号进行j层小波包分解,得到2j个频率段,定义在某一特定频率段的小波能量为该段小波系数的平方和Em
E m = Σ k | C m ( k ) | 2
得到总能量Etot的表达式为:
E t o t = Σ m Σ k | C m ( k ) | 2 = Σ m E m
定义某频段的相对小波包能量pm
p m = E m E t o t
定义相对小波包能量特征向量Ki
K i = ( p 0 , p 1 , . . . , p 2 j - 1 )
上式中,m=0,1,……,2j-1,
最后根据特征向量建立特征向量库。
3.根据权利要求1所述的变压器有载分接开关机械故障诊断方法,其特征在于:所述步骤3的具体实现方法为:
采用待监测的特征向量KT与标准特征向量库中特征向量Ki的距离d(KT,Ki)作为模式判别的依据,如果:
d(KT,Ki)=min{d(KT,Ki),...,d(KT,Kn)}
则将待监测向量的状态归入正常状态,采用如下欧几里德距离函数计算距离:
d ( K T , K i ) = Σ k = 0 2 j - 1 ( p T , k - p i , k )
式中,j为小波包分解的层数;
在状态检测过程中,取如下两组数据:一组是待检测时刻的变压器有载分接开关振动数据,用来建立测试数据库,另一组是变压器有载分接开关正常工作时的振动数据,用来建立标准数据库;当欧几里德距离小于某一阈值时,表明测试数据库中特征向量和标准数据库中的特征向量表征同一状态,为正常工作状态,当距离大于某一阈值时,表明该时刻有载分接开关可能出现机械故障。
CN201610632448.4A 2016-08-03 2016-08-03 变压器有载分接开关机械故障诊断方法 Pending CN106324490A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610632448.4A CN106324490A (zh) 2016-08-03 2016-08-03 变压器有载分接开关机械故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610632448.4A CN106324490A (zh) 2016-08-03 2016-08-03 变压器有载分接开关机械故障诊断方法

Publications (1)

Publication Number Publication Date
CN106324490A true CN106324490A (zh) 2017-01-11

Family

ID=57739698

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610632448.4A Pending CN106324490A (zh) 2016-08-03 2016-08-03 变压器有载分接开关机械故障诊断方法

Country Status (1)

Country Link
CN (1) CN106324490A (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107273585A (zh) * 2017-05-25 2017-10-20 国网山东省电力公司青岛供电公司 一种有载分接开关故障检测方法及装置
CN107907320A (zh) * 2017-12-01 2018-04-13 国网重庆市电力公司电力科学研究院 变压器有载分接开关的机械状态特征的诊断方法及装置
CN108362485A (zh) * 2018-01-15 2018-08-03 国网江苏省电力公司检修分公司特高压交直流运检中心 一种换流变分接开关典型机械故障模拟系统及故障模拟方法
CN108535638A (zh) * 2018-01-25 2018-09-14 国网浙江省电力有限公司电力科学研究院 基于多层滤波器的有载分接开关机械状态监测方法及系统
CN108593095A (zh) * 2018-04-26 2018-09-28 盐城博鸣信息科技有限公司 一种换流变压器分接开关的振动故障信号采集处理方法
CN109443717A (zh) * 2018-10-16 2019-03-08 国网江苏省电力有限公司南京供电分公司 一种有载分接开关机械故障在线监测方法
CN109443528A (zh) * 2018-11-16 2019-03-08 国网江苏省电力有限公司盐城供电分公司 一种基于振动信号分析的变压器故障诊断系统及其诊断方法
CN109443190A (zh) * 2018-11-20 2019-03-08 武汉拓清科技有限公司 基于暂态行波的变压器绕组变形在线监测方法及装置
CN109506761A (zh) * 2018-06-12 2019-03-22 国网四川省电力公司乐山供电公司 一种变压器表面振动特征提取方法
CN109541455A (zh) * 2018-12-03 2019-03-29 国网江苏省电力有限公司南京供电分公司 一种基于s变换时频谱svd降噪的oltc冲击特性提取方法
CN109856530A (zh) * 2018-12-25 2019-06-07 国网江苏省电力有限公司南京供电分公司 一种有载分接开关在线监测故障诊断方法
CN109900465A (zh) * 2019-03-12 2019-06-18 西安交通大学 一种基于声音信号的有载分接开关故障诊断方法
CN110132566A (zh) * 2019-05-28 2019-08-16 河海大学 一种基于模糊聚类的oltc故障诊断方法
CN110187268A (zh) * 2019-06-14 2019-08-30 国网江苏省电力有限公司 变压器有载分接开关数据管理的方法、系统和终端设备
CN114068209A (zh) * 2021-11-01 2022-02-18 中国南方电网有限责任公司超高压输电公司检修试验中心 一种分接开关档位选择方法、系统、存储介质及选择器
EP4163616A1 (en) * 2021-10-06 2023-04-12 Hitachi Energy Switzerland AG Monitoring of a power device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06117917A (ja) * 1992-10-02 1994-04-28 Shinko Electric Co Ltd 機械振動系振動応答特性計測方法およびその計測方法を用いた計測装置
CN101661080A (zh) * 2009-09-23 2010-03-03 河海大学 电力变压器有载分接开关触头脱落故障诊断方法
CN101825489A (zh) * 2010-01-29 2010-09-08 浙江大学 电力变压器有载分接开关振动信号分离方法
TW201100773A (en) * 2009-06-23 2011-01-01 Univ Southern Taiwan Tech Envelope signals derived from application of Morlet wavelet in gain of mechanical vibration signals and method for defect diagnosis of machines using exponential decay frequency of envelope signals
CN202066583U (zh) * 2011-03-17 2011-12-07 中国电力科学研究院 一种变压器有载分接开关的振动测量系统
CN105510812A (zh) * 2015-10-12 2016-04-20 国家电网公司 一种换流变分接开关状态监测系统及非平稳信号频谱估计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06117917A (ja) * 1992-10-02 1994-04-28 Shinko Electric Co Ltd 機械振動系振動応答特性計測方法およびその計測方法を用いた計測装置
TW201100773A (en) * 2009-06-23 2011-01-01 Univ Southern Taiwan Tech Envelope signals derived from application of Morlet wavelet in gain of mechanical vibration signals and method for defect diagnosis of machines using exponential decay frequency of envelope signals
CN101661080A (zh) * 2009-09-23 2010-03-03 河海大学 电力变压器有载分接开关触头脱落故障诊断方法
CN101825489A (zh) * 2010-01-29 2010-09-08 浙江大学 电力变压器有载分接开关振动信号分离方法
CN202066583U (zh) * 2011-03-17 2011-12-07 中国电力科学研究院 一种变压器有载分接开关的振动测量系统
CN105510812A (zh) * 2015-10-12 2016-04-20 国家电网公司 一种换流变分接开关状态监测系统及非平稳信号频谱估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
胡昌华等: "《基于MATLAB6.X的系统分析与设计-小波分析》", 31 January 2004, 西安电子科技大学出版社 *
虞和济 等: "《设备故障诊断工程》", 30 June 2001, 冶金工业出版社 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107273585A (zh) * 2017-05-25 2017-10-20 国网山东省电力公司青岛供电公司 一种有载分接开关故障检测方法及装置
CN107273585B (zh) * 2017-05-25 2021-04-13 国网山东省电力公司青岛供电公司 一种有载分接开关故障检测方法及装置
CN107907320B (zh) * 2017-12-01 2020-04-10 国网重庆市电力公司电力科学研究院 变压器有载分接开关的机械状态特征的诊断方法及装置
CN107907320A (zh) * 2017-12-01 2018-04-13 国网重庆市电力公司电力科学研究院 变压器有载分接开关的机械状态特征的诊断方法及装置
CN108362485A (zh) * 2018-01-15 2018-08-03 国网江苏省电力公司检修分公司特高压交直流运检中心 一种换流变分接开关典型机械故障模拟系统及故障模拟方法
CN108362485B (zh) * 2018-01-15 2020-10-20 国网江苏省电力有限公司检修分公司 一种换流变分接开关典型机械故障模拟系统及故障模拟方法
CN108535638A (zh) * 2018-01-25 2018-09-14 国网浙江省电力有限公司电力科学研究院 基于多层滤波器的有载分接开关机械状态监测方法及系统
CN108593095A (zh) * 2018-04-26 2018-09-28 盐城博鸣信息科技有限公司 一种换流变压器分接开关的振动故障信号采集处理方法
CN109506761B (zh) * 2018-06-12 2021-08-27 国网四川省电力公司乐山供电公司 一种变压器表面振动特征提取方法
CN109506761A (zh) * 2018-06-12 2019-03-22 国网四川省电力公司乐山供电公司 一种变压器表面振动特征提取方法
CN109443717A (zh) * 2018-10-16 2019-03-08 国网江苏省电力有限公司南京供电分公司 一种有载分接开关机械故障在线监测方法
CN109443528A (zh) * 2018-11-16 2019-03-08 国网江苏省电力有限公司盐城供电分公司 一种基于振动信号分析的变压器故障诊断系统及其诊断方法
CN109443190A (zh) * 2018-11-20 2019-03-08 武汉拓清科技有限公司 基于暂态行波的变压器绕组变形在线监测方法及装置
CN109541455A (zh) * 2018-12-03 2019-03-29 国网江苏省电力有限公司南京供电分公司 一种基于s变换时频谱svd降噪的oltc冲击特性提取方法
CN109856530A (zh) * 2018-12-25 2019-06-07 国网江苏省电力有限公司南京供电分公司 一种有载分接开关在线监测故障诊断方法
CN109856530B (zh) * 2018-12-25 2021-11-02 国网江苏省电力有限公司南京供电分公司 一种有载分接开关在线监测故障诊断方法
CN109900465A (zh) * 2019-03-12 2019-06-18 西安交通大学 一种基于声音信号的有载分接开关故障诊断方法
CN109900465B (zh) * 2019-03-12 2020-07-28 西安交通大学 一种基于声音信号的有载分接开关故障诊断方法
CN110132566A (zh) * 2019-05-28 2019-08-16 河海大学 一种基于模糊聚类的oltc故障诊断方法
CN110187268A (zh) * 2019-06-14 2019-08-30 国网江苏省电力有限公司 变压器有载分接开关数据管理的方法、系统和终端设备
EP4163616A1 (en) * 2021-10-06 2023-04-12 Hitachi Energy Switzerland AG Monitoring of a power device
WO2023057514A1 (en) * 2021-10-06 2023-04-13 Hitachi Energy Switzerland Ag Monitoring of a power device
CN114068209A (zh) * 2021-11-01 2022-02-18 中国南方电网有限责任公司超高压输电公司检修试验中心 一种分接开关档位选择方法、系统、存储介质及选择器
CN114068209B (zh) * 2021-11-01 2023-12-15 中国南方电网有限责任公司超高压输电公司检修试验中心 一种分接开关档位选择方法、系统、存储介质及选择器

Similar Documents

Publication Publication Date Title
CN106324490A (zh) 变压器有载分接开关机械故障诊断方法
CN104808109B (zh) 基于录波数据的高压输电线路故障识别方法和系统
CN103499753B (zh) 快速识别高压直流输电系统暂态故障的智能系统
CN111198309B (zh) 一种海上平台干式变压器状态监测与故障预警方法
WO2016058435A1 (zh) 变压器外部故障下绕组状态评估方法
CN104849628B (zh) 一种电缆附件故障分析检测方法
CN106908674A (zh) 一种基于多状态量预测的变压器状态评估方法
CN110146269A (zh) 基于eemd的多尺度模糊熵的oltc故障诊断方法
CN105676121A (zh) 一种220kV及以上瓷柱式断路器闪络故障评估方法
CN103278715B (zh) 电力设备试验方法
CN103776622B (zh) 一种高压断路器机械故障快速诊断系统
CN103513125A (zh) 一种220kv级以上变压器一体化智能诊断系统及其方法
CN102567783A (zh) 并联混合型电能质量调节器专家故障分析与诊断方法
CN106443433B (zh) 基于分合闸线圈电流的断路器状态监测系统及监测方法
CN108664741A (zh) 基于时间序列模型特征的变电站故障检测方法
CN109507580A (zh) 一种高压断路器电气寿命监测方法及系统
CN103499770B (zh) 基于广域电压的电网故障定位方法
CN103578042B (zh) 一种电力变压器可靠度层次化评估方法
CN103630810A (zh) 一种并联补偿电容器绝缘状态诊断方法
CN104360194A (zh) 一种智能电网故障诊断方法
CN102288881B (zh) 变压器油纸绝缘尖刺放电缺陷严重程度诊断方法
CN102147440A (zh) 变压器油纸绝缘沿面局部放电缺陷严重程度诊断方法
CN106338673A (zh) 一种变压器绕组故障特征参量提取方法
CN108152733B (zh) 一种用于高压断路器机械特性的带电检测分合闸时间方法
CN112345885A (zh) 一种高压直流输电线路故障检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20170111

RJ01 Rejection of invention patent application after publication