CN105548644A - 一种双线并绕的罗氏线圈及信号处理方法 - Google Patents

一种双线并绕的罗氏线圈及信号处理方法 Download PDF

Info

Publication number
CN105548644A
CN105548644A CN201510875045.8A CN201510875045A CN105548644A CN 105548644 A CN105548644 A CN 105548644A CN 201510875045 A CN201510875045 A CN 201510875045A CN 105548644 A CN105548644 A CN 105548644A
Authority
CN
China
Prior art keywords
coil
luo
current
double
tap
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
CN201510875045.8A
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
North China Electric Power University
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Liaoning Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
North China Electric Power University
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Liaoning 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, North China Electric Power University, Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Liaoning Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510875045.8A priority Critical patent/CN105548644A/zh
Publication of CN105548644A publication Critical patent/CN105548644A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R15/00Details of measuring arrangements of the types provided for in groups G01R17/00 - G01R29/00, G01R33/00 - G01R33/26 or G01R35/00
    • G01R15/14Adaptations providing voltage or current isolation, e.g. for high-voltage or high-current networks
    • G01R15/18Adaptations providing voltage or current isolation, e.g. for high-voltage or high-current networks using inductive devices, e.g. transformers
    • G01R15/181Adaptations providing voltage or current isolation, e.g. for high-voltage or high-current networks using inductive devices, e.g. transformers using coils without a magnetic core, e.g. Rogowski coils

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Instrument Details And Bridges, And Automatic Balancing Devices (AREA)

Abstract

本发明涉及电子技术领域,尤其涉及一种双线并绕的罗氏线圈及信号处理方法,具体是电子式电流互感器中的罗氏线圈绕制方法和原始信号的检测方法。本发明包括由双线并绕的罗氏线圈,采用双路漆包线紧挨、平行、并均匀缠绕在罗氏线圈骨架上,形成环形线圈。本发明不用改变线圈结构,仅在绕制方法和连接方式上做了改进,即可获得良好的性能。由于采用双路漆包线或其它表层绝缘的导线并绕在传统罗氏线圈骨架上,获得两路差动信号,经此差动信号的检测,灵敏度提高了一倍,还具有抗共模干扰的特点。可应用于电站的计量仪表、继电保护,可控硅整流,变频调速,电工领域的测量仪器等场合。在电力系统等领域具有广阔的应用前景。

Description

一种双线并绕的罗氏线圈及信号处理方法
技术领域
本发明涉及电子技术领域,尤其涉及一种双线并绕的罗氏线圈及信号处理方法,具体是电子式电流互感器中的罗氏线圈绕制方法和原始信号的检测方法。
背景技术
罗氏线圈(RogwsikiCoil,罗哥夫斯基线圈)又叫电流测量线圈、微分电流传感器,作为一种新型电子式电流互感器测量装置,主要特点是被测电流与线圈无直接接触,具有测量范围宽、无磁饱和、体积小等优点,不含铁磁性材料,无磁滞效应,几乎为零的相位误差;无磁饱和象,因而测量范围可从数安培到数百千安的电流;结构简单,并且和被测电流之间没有直接的电路联系;响应频带宽0.1Hz-1MHz。与带铁芯的传统互感器相比,罗氏线圈具有测量范围宽,精度高,稳定可靠,响应频带宽。
相关研究结果表明,电子设备对温度变化十分敏感,无论是高温或是低温都会对其性能造成影响。电子式互感器的工作环境在室外,四季温差较大,这样会对互感器的测量精度造成影响。实验表明,在温度变化时,其输出电压在一定的范围内浮动。
对于小电流测量,感应出的电压信号很微弱(为几毫伏至几十毫伏),极易受干扰而难以精确测量;而且其精度都很低,最高达到0.2级。
发明内容
为了解决上述现有技术中存在的问题,本发明提供了一种双线并绕的罗氏线圈及信号处理方法,目的是可以改善电流互感器对温度的变化引起测量精度的影响,同时提高了检测灵敏度。
为了实现上述发明目的,本发明的技术方案是这样实现的:
一种双线并绕的罗氏线圈及信号处理方法,包括由双线并绕的罗氏线圈,采用双路漆包线紧挨、平行、并均匀缠绕在罗氏线圈骨架上,形成环形线圈;
信号处理方法如下:利用前置差动信号电压放大器对e(t)处理;差动信号的处理采用差动式放大器,集放大、积分和相位调整为一体;通过调整比值a,可以调整差动放大器的放大倍数,通过调整比值b,可以微调整放大器输入和输出之间的相位。
所述的双路漆包线用表层绝缘的导线代替。
所述的骨架为非铁磁性材料。
所述的线圈包括:线圈A和线圈B,和四根抽头:抽头A1,抽头A2,抽头B1,抽头B2分别为线圈A的头和尾;抽头B1和抽头B2分别为线圈B的头和尾;线圈A的头与线圈B的尾联接作为参考点,即抽头A1与抽头B2联接;在线圈A的尾抽头A2和线圈B的头抽头B1获得两路差动信号。
所述的骨架为矩形结构的线圈时,其线圈端电压和线圈电流的计算方法为:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流;
对于截面为矩形的罗氏线圈:
M = μ 0 N h 2 π l n ( b a ) - - - ( 3 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,ab为线圈横截面的内外径,h为线圈截面高度。
所述的骨架为圆形结构时,其线圈端电压和线圈电流的计算方法为:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流;
对于截面为圆形的罗氏线圈:
M = μ 0 N S l - - - ( 4 )
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,S为每匝线圈截面面积,l(2πr)为线圈环形骨架芯平均周长。
所述的线圈,其中线圈B1、线圈B2输出电压为e1(t),线圈A1、线圈A2输出电压为e2(t),o(t)为前置放大器输出,比值a和电阻值R确定前置放大器的放大倍数,比值b和电阻值R1可微调输入/输出之间的相差;
计算方法如下:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流;
对于截面为矩形的罗氏线圈:
M = μ 0 N h 2 π l n ( b a ) - - - ( 3 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,a、b为线圈横截面的内外径,h为线圈截面高度;
对于截面为园的罗氏线圈:
M = μ 0 N S l - - - ( 4 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,S为每匝线圈截面面积,l(2πr)为线圈环形骨架芯平均周长。
本发明的优点及有益效果是:
本发明的优点是不用改变线圈结构,仅在绕制方法和连接方式上做了改进,即可获得良好的性能。又由于采用双路漆包线或其它表层绝缘的导线并绕在传统罗氏线圈骨架上,获得两路差动信号,经此差动信号的检测,使灵敏度提高了一倍。本发明还具有抗共模干扰的特点。本发明罗氏线圈可应用于电站的计量仪表、继电保护,可控硅整流,变频调速,电工领域的测量仪器等场合。在电力系统等领域具有广阔的应用前景,成为互感器发展的重要方向之一。
以下结合附图和具体实施例,对本发明作进一步详细的说明,但不受本实施例所限。
附图说明
图1是骨架为矩形结构的线圈示意图;
图2是骨架为圆形结构的线圈示意图;
图3是本发明线圈结构抽头连接示意图;
图4是本发明信号检测电路图;
图5是本发明线圈与前置放大器连接示意图;
图6是本发明差动式放大、积分、移相电路图,R5=10Ω;
图7是输入输出相位比较图;
图8是R5=100kΩ,电路图;
图9是输入输出相位比较图;
图10是双线并绕罗氏线圈连接后的等效电路;
图11是是罗氏线圈连接前置放大器示意图。
图1中:双路漆包线1,骨架2,线圈A,线圈B,抽头A1,抽头A2,抽头B1,抽头B2,线圈A1,线圈B1、,线圈A2,线圈B2。
具体实施方式
如图1-图3所示,本发明是一种双线并绕的罗氏线圈及信号处理方法,所述的双线并绕的罗氏线圈是基于在传统罗氏线圈骨架上,采用双路漆包线1或其它表层绝缘的导线,紧挨、平行、并均匀缠绕在骨架2上,形成环形线圈。所述的骨架2为非铁磁性材料。这样获两组线圈A和线圈B,和四根抽头:抽头A1,抽头A2,抽头B1,抽头B2分别为线圈A的头和尾;抽头B1和抽头B2分别为线圈B的头和尾。线圈A的头与线圈B的尾联接作为参考点,即抽头A1与抽头B2联接。在线圈A的尾抽头A2和线圈B的头抽头B1获得两路差动信号。差动信号的检测,灵敏度提高了一倍。差动式积分放大电路。e(t)为罗氏线圈输出的差动电压。
如图1所示,图1是骨架为矩形结构的线圈示意图;其线圈端电压和线圈电流的计算方法为:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流。
对于截面为矩形的罗氏线圈:
M = μ 0 N h 2 π l n ( b a ) - - - ( 3 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,a、b为线圈横截面的内外径,h为线圈截面高度。
图2是骨架为圆形结构的线圈示意图;其线圈端电压和线圈电流的计算方法为:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流。
对于截面为圆形的罗氏线圈:
M = μ 0 N S l - - - ( 4 )
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,S为每匝线圈截面面积,l(2πr)为线圈环形骨架芯平均周长。
图3是本发明线圈结构抽头连接示意图。
图4是本发明信号检测电路图。一种基于集放大、积分和相位调整为一体的差动式前置放大器。
本发明的信号处理方法如下:
利用前置差动信号电压放大器对e(t)处理。
差动信号的处理采用差动式放大器,集放大、积分和相位调整为一体。通过调整比值a,可以调整差动放大器的放大倍数,通过调整比值b,可以微调整放大器输入和输出之间的相位。如图5所示,图5是差动式罗氏线圈与差动式前置放大器连接示意图。
输出信号是电流对时间的微分。通过一个对输出的电压信号进行积分的电路,就可以真实还原测量电流,如图4所示,图4是本发明信号检测电路图。
将几十毫伏级的原始信号放大到伏级。差动式放大器可以抑制零点漂移,由于它是一个对称电路,可使漂移信号相互抵消。由温度变化等因素引起两路的输出漂移电压必然是大小相等,极性相同,即为共模信号。具有抗共模干扰的特点。
图4和图5中,线圈B1、线圈B2输出电压为e1(t),线圈A1、线圈A2输出电压为e2(t),o(t)为前置放大器输出,比值a和电阻值R确定前置放大器的放大倍数,比值b和电阻值R1可微调输入/输出之间的相差。
本发明的计算方法如下:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流。
对于截面为矩形的罗氏线圈:图1所示。
M = μ 0 N h 2 π l n ( b a ) - - - ( 3 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,a、b为线圈横截面的内外径,h为线圈截面高度。
对于截面为园的罗氏线圈:图2所示。
M = μ 0 N S l - - - ( 4 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,S为每匝线圈截面面积,l(2πr)为线圈环形骨架芯平均周长。
以下是检测电路与输入输出波形:
图6为差动式放大、积分、移相电路图,当前置放大器输入V1=20mV;V2=20mV的差动信号;其它元件参数如图6所示。当R5=10Ω时,在示波器XSC2可观察输入/输出波形。图5为输入/输出相位比较,红为输入,黑为输出。
图7是输入输出相位比较图,可以看出输出相位超前输入90﹣Δ。当前置放大器输入V1=20mV;V2=20mV的差动信号;其它元件参数如图8所示。当R5=100kΩ时,在示波器XSC2可观察输入/输出波形。图7为输入/输出相位比较,红为输入,黑为输出。
图9是输入输出相位比较图,可以看出输出相位超前输入90﹢Δ。改变R5阻值,可微调输入/输出相位差。
图10是双线并绕罗氏线圈连接后的等效电路。e1(t)、e2(t)为被测电流在两组双线并绕线圈中的感应电动势;L1、L2和r1、r2分别表示罗氏线圈的自感和内阻;C1、C2线圈的分布;Rt1和Rt2线圈端接电阻。
图11是罗氏线圈连接前置放大器,ein(t)=e′1(t)-e′2(t)为罗氏线圈输出电压;调整a和电阻R1确定前置放大器的放大倍数;调整b和电阻R确定输出u0(t)与输入ein(t)之间的相位差。

Claims (7)

1.一种双线并绕的罗氏线圈及信号处理方法,其特征是:包括由双线并绕的罗氏线圈,采用双路漆包线(1)紧挨、平行、并均匀缠绕在罗氏线圈骨架(2)上,形成环形线圈;
信号处理方法如下:
利用前置差动信号电压放大器对e(t)处理;
差动信号的处理采用差动式放大器,集放大、积分和相位调整为一体;通过调整比值a,可以调整差动放大器的放大倍数,通过调整比值b,可以微调整放大器输入和输出之间的相位。
2.根据权利要求1所述的一种双线并绕的罗氏线圈及信号处理方法,其特征是:所述的双路漆包线(1)用表层绝缘的导线代替。
3.根据权利要求1所述的一种双线并绕的罗氏线圈及信号处理方法,其特征是:所述的骨架(2)为非铁磁性材料。
4.根据权利要求1所述的一种双线并绕的罗氏线圈及信号处理方法,其特征是:所述的线圈包括:线圈A和线圈B,和四根抽头:抽头A1,抽头A2,抽头B1,抽头B2分别为线圈A的头和尾;抽头B1和抽头B2分别为线圈B的头和尾;线圈A的头与线圈B的尾联接作为参考点,即抽头A1与抽头B2联接;在线圈A的尾抽头A2和线圈B的头抽头B1获得两路差动信号。
5.根据权利要求1所述的一种双线并绕的罗氏线圈及信号处理方法,其特征是:所述的骨架(2)为矩形结构的线圈时,其线圈端电压和线圈电流的计算方法为:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流;
对于截面为矩形的罗氏线圈:
M = μ 0 N h 2 π l n ( b a ) - - - ( 3 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,ab为线圈横截面的内外径,h为线圈截面高度。
6.根据权利要求1所述的一种双线并绕的罗氏线圈及信号处理方法,其特征是:所述的骨架(2)为圆形结构时,其线圈端电压和线圈电流的计算方法为:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流;
对于截面为圆形的罗氏线圈:
M = μ 0 N S l - - - ( 4 )
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,S为每匝线圈截面面积,l(2πr)为线圈环形骨架芯平均周长。
7.根据权利要求1所述的一种双线并绕的罗氏线圈及信号处理方法,其特征是:所述的线圈,其中线圈B1、线圈B2输出电压为e1(t),线圈A1、线圈A2输出电压为e2(t),o(t)为前置放大器输出,比值a和电阻值R确定前置放大器的放大倍数,比值b和电阻值R1可微调输入/输出之间的相差;
计算方法如下:
e ( t ) = - M d i d t - - - ( 1 ) ;
i ( t ) = - 1 M ∫ e ( t ) d t - - - ( 2 ) ;
线圈的端电压e(t)与一次电流的时间变化率成正比,与罗氏线圈和载流导体之间的互感M成正比,i(t)载流导体电流;
对于截面为矩形的罗氏线圈:
M = μ 0 N h 2 π l n ( b a ) - - - ( 3 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,a、b为线圈横截面的内外径,h为线圈截面高度;
对于截面为园的罗氏线圈:
M = μ 0 N S l - - - ( 4 ) ;
μ0=4π×10-7H/m为真空磁导率,N为线圈匝数,S为每匝线圈截面面积,l(2πr)为线圈环形骨架芯平均周长。
CN201510875045.8A 2015-12-02 2015-12-02 一种双线并绕的罗氏线圈及信号处理方法 Pending CN105548644A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510875045.8A CN105548644A (zh) 2015-12-02 2015-12-02 一种双线并绕的罗氏线圈及信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510875045.8A CN105548644A (zh) 2015-12-02 2015-12-02 一种双线并绕的罗氏线圈及信号处理方法

Publications (1)

Publication Number Publication Date
CN105548644A true CN105548644A (zh) 2016-05-04

Family

ID=55827967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510875045.8A Pending CN105548644A (zh) 2015-12-02 2015-12-02 一种双线并绕的罗氏线圈及信号处理方法

Country Status (1)

Country Link
CN (1) CN105548644A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107768124A (zh) * 2017-11-08 2018-03-06 国网辽宁省电力有限公司电力科学研究院 一种带差动放大电路的双绕线罗氏线圈电子式电流互感器
CN108414812A (zh) * 2018-02-06 2018-08-17 国网安徽省电力有限公司 一种基于罗氏线圈的电子式电流互感器及其特性分析方法
CN112730944A (zh) * 2020-12-19 2021-04-30 华中科技大学 一种基于罗氏线圈的电流测量方法及装置
CN114814702A (zh) * 2022-01-25 2022-07-29 武汉大学 双环骨架罗氏线圈电流传感器灵敏度的快速计算方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107768124A (zh) * 2017-11-08 2018-03-06 国网辽宁省电力有限公司电力科学研究院 一种带差动放大电路的双绕线罗氏线圈电子式电流互感器
CN108414812A (zh) * 2018-02-06 2018-08-17 国网安徽省电力有限公司 一种基于罗氏线圈的电子式电流互感器及其特性分析方法
CN112730944A (zh) * 2020-12-19 2021-04-30 华中科技大学 一种基于罗氏线圈的电流测量方法及装置
CN112730944B (zh) * 2020-12-19 2022-03-11 华中科技大学 一种基于罗氏线圈的电流测量方法及装置
CN114814702A (zh) * 2022-01-25 2022-07-29 武汉大学 双环骨架罗氏线圈电流传感器灵敏度的快速计算方法
CN114814702B (zh) * 2022-01-25 2024-05-14 武汉大学 双环骨架罗氏线圈电流传感器灵敏度的快速计算方法

Similar Documents

Publication Publication Date Title
CN103235170B (zh) 差分式D-dot电压传感器
CN103575960A (zh) 巨磁阻效应电流传感器
CN106771477B (zh) 大口径高灵敏度的高压直流电缆泄漏电流检测传感器
CN201689125U (zh) 零磁通电流传感器
CN103308743B (zh) 直流电流计量装置
CN105548644A (zh) 一种双线并绕的罗氏线圈及信号处理方法
CN104155499A (zh) 基于组合式pcb型罗氏线圈的电流测量装置及方法
CN110927428B (zh) 一种宽量程宽频高精度磁平衡式电流测量装置
CN102129059A (zh) 电流互感器校验用5ka零磁道式直流电流比较仪
CN103630724B (zh) 基于带间隙铁芯线圈的自校准谐波电流互感器
CN203838218U (zh) 一种电子式电流互感器
CN104851580A (zh) 基于磁位计补偿的带间隙铁芯式罗氏线圈互感器
CN110824229A (zh) 一种单磁芯多绕组磁平衡式电流检测装置
CN105510673A (zh) 一种直流电流测量装置
CN204464036U (zh) 基于磁位计补偿的带间隙铁芯式罗氏线圈互感器
CN210803568U (zh) 一种交流电流传感器
CN101872006B (zh) 带升压器的电子式周期性非正弦波基准电压互感器
Fubin et al. Effects of geometrical parameters on the performance of Rogowski coil for current measuring
CN107768124A (zh) 一种带差动放大电路的双绕线罗氏线圈电子式电流互感器
CN110456142A (zh) 磁调制式直流漏电流传感器
CN204666709U (zh) 半磁芯电流传感器
CN113341211B (zh) 一种磁通门电流传感器及其测量方法
CN108226610A (zh) 一种用于交直流小电流的钳型测量装置
CN114675064A (zh) 一种高精度双量程磁通门电流传感器
CN102998522B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20160504

RJ01 Rejection of invention patent application after publication