CN106203766A - 一种500kV变电站串抗区域工频电场优化方法 - Google Patents

一种500kV变电站串抗区域工频电场优化方法 Download PDF

Info

Publication number
CN106203766A
CN106203766A CN201610460270.XA CN201610460270A CN106203766A CN 106203766 A CN106203766 A CN 106203766A CN 201610460270 A CN201610460270 A CN 201610460270A CN 106203766 A CN106203766 A CN 106203766A
Authority
CN
China
Prior art keywords
phase sequence
electric field
sequence arrangement
frequency electric
region
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
CN201610460270.XA
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.)
Shenzhen Power Supply Co ltd
Chongqing University
Original Assignee
Shenzhen Power Supply Co ltd
Chongqing University
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 Shenzhen Power Supply Co ltd, Chongqing University filed Critical Shenzhen Power Supply Co ltd
Priority to CN201610460270.XA priority Critical patent/CN106203766A/zh
Publication of CN106203766A publication Critical patent/CN106203766A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • 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

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供一种500kV变电站串抗区域工频电场优化方法,包括确定500kV变电站串抗区域的相关信息,并根据相关信息构建计算模型;针对预设的双回母线相序排列范围中每一相序排列,均采用边界元法对计算模型进行仿真求解运算,得到每一相序排列分别对应的工频电场强度,并筛选出满足一定条件的工频电场强度及其对应的相序排列;获取500kV变电站串抗区域双回母线的当前相序排列,并将当前相序排列修订为筛选出的相序排列。实施本发明,分析及调整串抗区域的双回高压线路的相序排列,从而降低串抗区域的工频电场强度,达到优化电磁环境的目的。

Description

一种500kV变电站串抗区域工频电场优化方法
技术领域
本发明涉及输电线路电磁环境检测技术领域,尤其涉及一种500kV变电站串抗区域工频电场优化方法。
背景技术
随着国民经济的快速发展,电力需求越来越迫切,然而现阶段能源分配并不合理,需求大的负荷中心一般都离能源基地较远。为了快速实现资源优化配置,通常会在能源基地与负荷中心之间铺设高压输电线路,这样有利于经济和环境效益的提升,但是高压输电线路的高电压、大电流特点带来了一系列的电磁环境问题。
国内外学者及研究机构发现,500kV变电站串抗区域的双回高压线路所产生的电磁环境影响也比较严重,不仅受到线路结构的影响,还受到运行情况的影响,因此均作了大量的研究。例如美国邦纳维尔电力公司从1976年开始进行特高压线路的机械结构的研究,并进行了电晕和电场、生态和环境、操作和雷电冲击绝缘等方面的研究;例如日本,前苏联及加拿大等国家均相继对串抗区输电线路的电磁环境影响作了相关研究。
但是,发明人发现,国内外的研究主要侧重于串抗区域线路结构及运行情况对输电线路下工频电磁场的影响研究,而针对相序的影响情况,研究较少。
而随着我国变电站数量的增加,对500kV变电站串抗区域的电磁环境研究也需要进一步的提高。此时,对500kV变电站的串抗区域电磁环境优化设计方法研究就十分必要。
发明内容
本发明实施例所要解决的技术问题在于,提供一种500kV变电站串抗区域工频电场优化方法,分析及调整串抗区域的双回高压线路的相序排列,从而降 低串抗区域的工频电场强度,达到优化电磁环境的目的。
为了解决上述技术问题,本发明实施例提供了一种500kV变电站串抗区域工频电场优化方法,所述方法包括:
a、确定500kV变电站串抗区域的相关信息,并根据所述确定的相关信息构建计算模型;
b、针对预设的双回母线相序排列范围中每一相序排列,均采用边界元法对所述构建的计算模型进行仿真求解运算,得到每一相序排列分别对应的工频电场强度,并在所述得到的每一相序排列分别对应的工频电场强度中,筛选出满足一定条件的工频电场强度及其对应的相序排列;
c、获取所述500kV变电站串抗区域双回母线的当前相序排列,并将所述500kV变电站串抗区域双回母线的当前相序排列修订为所述筛选出的相序排列。
其中,所述步骤a中的“500kV变电站串抗区域的相关信息”包括双回母线的相序排列、双回母线距地的线路高度以及双回母线的相间位。
其中,所述步骤b中的“预设的双回母线相序排列范围”有六种相序排列,包括ABC-ABC、ABC-ACB、ABC-BAC、ABC-BCA、ABC-CAB、ABC-CBA。
其中,所述步骤b中的“每一相序排列分别对应的工频电场强度”通过对每一相序排列在所述构建的计算模型中取双回母线距地的线路高度均为1.5米时,仿真求解而获得的。
其中,所述步骤b中的“筛选出满足一定条件的工频电场强度及其对应的相序排列”具体为:筛选出最小的工频电场强度及其对应的相序排列。
实施本发明实施例,具有如下有益效果:
在本发明实施例中,由于针对每一相序排列均采用边界元法进行仿真求解运算,筛选出满足一定条件的工频电场强度及其对应的相序排列,并将所筛选的相序排列用于调整串抗区域双回高压线路的当前相序排列,从而降低串抗区域的工频电场强度,达到优化电磁环境的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施 例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,根据这些附图获得其他的附图仍属于本发明的范畴。
图1为本发明实施例提供的500kV变电站串抗区域工频电场优化方法的流程图;
图2为本发明实施例提供的500kV变电站串抗区域工频电场优化方法应用场景中电偶极子向量图;
图3为本发明实施例提供的500kV变电站串抗区域工频电场优化系统的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
如图1所示,为本发明实施例中,提供的一种500kV变电站串抗区域工频电场优化方法,所述方法包括:
步骤S1、确定500kV变电站串抗区域的相关信息,并根据所述确定的相关信息构建计算模型;
具体过程为,确定500kV变电站串抗区域的相关信息,包括双回母线的相序排列、双回母线距地的线路高度以及双回母线的相间位置等,并根据上述相关信息构建计算模型。
应当说明的是,双回母线应在串抗区域的上方。
步骤S2、针对预设的双回母线相序排列范围中每一相序排列,均采用边界元法对所述构建的计算模型进行仿真求解运算,得到每一相序排列分别对应的工频电场强度,并在所述得到的每一相序排列分别对应的工频电场强度中,筛选出满足一定条件的工频电场强度及其对应的相序排列;
具体过程为,预先定义双回母线相序排列范围,该相序排列范围包括ABC-ABC、ABC-ACB、ABC-BAC、ABC-BCA、ABC-CAB、ABC-CBA等六种相序排列,并对每一种均采用边界元法对计算模型进行仿真求解运算,得到 每一相序排列分别对应的工频电场强度,并进一步每一相序排列分别对应的工频电场强度中,筛选出满足一定条件的工频电场强度及其对应的相序排列。
在一个实施例中,在计算模型中取双回母线下方距地高度均为1.5米仿真求解后,得到六种双回母线相序排列距地面1.5m处分别对应的工频电场强度,并筛选出最小的工频电场强度及其对应的相序排列。
在本发明实施例中,由于间接边界积分是基于点电荷的基本解和叠加原理建立边界积分方程的,因此边界元计算方法具体如下:
(I)根据单层源的观点,认为带电体表面的电荷按面密度σ(x,y,z)分布,若带电体位于有限空间,且以无限远处为电位参考点,则场域中任意点M的电位按叠加原理可写成,如下述公式(1)所示:
φ M = ∫ s σ 4 πϵr 1 d S - - - ( 1 )
式(1)中,r1为源点和场点之间的距离。
由于导体表面的因此式(1)可以改写为:
式(2)为单层源的间接边界积分方程,适用的条件为各向同性、线性的均匀介质。
(II)对于双层等效源,它和电位之间的关系式可推导为:
如图2所示,电偶极子的偶极距p=qΔl,其中Δl为等量而异号的两点电荷之间的微小距离,其方向从负电荷指向正电荷。
则点M处的电位为
由于Δl很小,可以用微分代替,因此,式(3)可改写为
因此对于面密度为σp的双层面偶极距,产生的电位如公式(5)所示:
当场点与源点重合时,式(5)中便会出现奇点,这种奇点为可去奇点,具体如公式(6)所示:
式(6)中:np为曲面上p点的外法线方向上的单位矢量;为从p点指向M点的单位矢量,是面单元对场点M所张的立体角dΩ,因此,式(6)可改写为如公式(7)所示:
当场点M接近于源点p0时,将曲面上的p0点孤立起来,且引进符号S-P0表示除去P0点后的边界。
当场点M从偶极层的内部接近P0点时,则有
当场点M从偶极层的外部接近P0点时,则有
从式(8)和式(9)可以看出,通过偶极层时,电位将引起突变,具体如式(10)所示:
分层介质的间接边界积分方程:对于电力设备中的绝缘子,其表面的电位是未知的,根据分层介质分界面电通量连续性可以得出其分界面电位连续性条件,具体如式(11)所示:
式(11)中:εd为介质的介电常数;εa为空气的介电常数;εan为介质表面电场位于空气侧的法向分量;εdn为介质表面电场位于介质侧的法向分量。对于绝缘介质与空气分界面,根据式(8)以及立体角的定义可以得出,位于介质表面一侧点Pi的电位的法向导数为
对于位于空气一侧的点Pi的电位的法向导数为
式(12)和式(13)中:en为dS的法向方向;R表示源点到场点的位置矢量。式(12)和式(13)中的法向量方向由介质指向空气。将式(12)和式(13)代入(11)中可以得到介质分界面的间接积分方程为
1 2 ϵ ϵ d + ϵ a ϵ d - ϵ a δ i + 1 4 π ϵ ∫ S δ i R · e n R 3 d S = 0 - - - ( 14 )
悬浮导体的间接边界积分方程:悬浮导体边界面上节点的电位也是未知的,将悬浮导体看作介电常数为无穷大的介质,即将式(14)中的εd→∞,对式(14)取极限得到悬浮导体的边界积分方程如公式(15)所示:
1 2 ϵ 0 δ i + 1 4 π ϵ ∫ S δ i R · e n R 3 d S = 0 - - - ( 15 )
利用仿真软件对该串抗区域的电场进行建模仿真,综合仿真数据和计算结果,得出不同相序排列情况下串抗区的电场大小,再选择出满足一定条件的相序排列。
步骤S3、获取所述500kV变电站串抗区域双回母线的当前相序排列,并将所述500kV变电站串抗区域双回母线的当前相序排列修订为所述筛选出的相序排列。
具体过程为,对500kV变电站串抗区域双回母线的当前相序排列进行修订,从而使得所筛选的相序排列替换成双回母线的当前相序排列,达到减小串抗区域工频电场强度的目的。
如图3所示,为本发明实施例中,提供的一种500kV变电站串抗区域工频电场优化系统,所述系统包括:
计算模型构建单元210,用于确定500kV变电站串抗区域的相关信息,并根据所述确定的相关信息构建计算模型;
相序排列筛选单元220,用于针对预设的双回母线相序排列范围中每一相序排列,均采用边界元法对所述构建的计算模型进行仿真求解运算,得到每一相序排列分别对应的工频电场强度,并在所述得到的每一相序排列分别对应的工频电场强度中,筛选出满足一定条件的工频电场强度及其对应的相序排列;
相序排列修订单元230,用于获取所述500kV变电站串抗区域双回母线的当前相序排列,并将所述500kV变电站串抗区域双回母线的当前相序排列修订为所述筛选出的相序排列。
实施本发明实施例,具有如下有益效果:
在本发明实施例中,由于针对每一相序排列均采用边界元法进行仿真求解运算,筛选出满足一定条件的工频电场强度及其对应的相序排列,并将所筛选的相序排列用于调整串抗区域双回高压线路的当前相序排列,从而降低串抗区域的工频电场强度,达到优化电磁环境的目的。
值得注意的是,上述系统实施例中,所包括的各个系统单元只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可; 另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,所述的存储介质,如ROM/RAM、磁盘、光盘等。
以上所揭露的仅为本发明较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。

Claims (5)

1.一种500kV变电站串抗区域工频电场优化方法,其特征在于,所述方法包括:
a、确定500kV变电站串抗区域的相关信息,并根据所述确定的相关信息构建计算模型;
b、针对预设的双回母线相序排列范围中每一相序排列,均采用边界元法对所述构建的计算模型进行仿真求解运算,得到每一相序排列分别对应的工频电场强度,并在所述得到的每一相序排列分别对应的工频电场强度中,筛选出满足一定条件的工频电场强度及其对应的相序排列;
c、获取所述500kV变电站串抗区域双回母线的当前相序排列,并将所述500kV变电站串抗区域双回母线的当前相序排列修订为所述筛选出的相序排列。
2.如权利要求1所述的方法,其特征在于,所述步骤a中的“500kV变电站串抗区域的相关信息”包括双回母线的相序排列、双回母线距地的线路高度以及双回母线的相间位置。
3.如权利要求1所述的方法,其特征在于,所述步骤b中的“预设的双回母线相序排列范围”有六种相序排列,包括ABC-ABC、ABC-ACB、ABC-BAC、ABC-BCA、ABC-CAB、ABC-CBA。
4.如权利要求1所述的方法,其特征在于,所述步骤b中的“每一相序排列分别对应的工频电场强度”通过对每一相序排列在所述构建的计算模型中取双回母线下方距地的高度为1.5米时,仿真求解而获得的。
5.如权利要求1所述的方法,其特征在于,所述步骤b中的“筛选出满足一定条件的工频电场强度及其对应的相序排列”具体为:筛选出最小的工频电场强度及其对应的相序排列。
CN201610460270.XA 2016-06-23 2016-06-23 一种500kV变电站串抗区域工频电场优化方法 Pending CN106203766A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610460270.XA CN106203766A (zh) 2016-06-23 2016-06-23 一种500kV变电站串抗区域工频电场优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610460270.XA CN106203766A (zh) 2016-06-23 2016-06-23 一种500kV变电站串抗区域工频电场优化方法

Publications (1)

Publication Number Publication Date
CN106203766A true CN106203766A (zh) 2016-12-07

Family

ID=57461360

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610460270.XA Pending CN106203766A (zh) 2016-06-23 2016-06-23 一种500kV变电站串抗区域工频电场优化方法

Country Status (1)

Country Link
CN (1) CN106203766A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107784187A (zh) * 2017-11-22 2018-03-09 重庆大学 一种基于边界元法的变电站电磁环境分布研究方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2535809Y (zh) * 2000-12-29 2003-02-12 重庆德通科技开发有限公司 50Hz工频电场无线遥测仪
US20140052393A1 (en) * 2011-04-29 2014-02-20 Shanxi Electric Power Research Institute Method for determining spatial location of conducting wire and aerial earth wire of power transmission line
CN103810329A (zh) * 2014-01-18 2014-05-21 国家电网公司 根据输电线路功率调整输电线路导线相序排列的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2535809Y (zh) * 2000-12-29 2003-02-12 重庆德通科技开发有限公司 50Hz工频电场无线遥测仪
US20140052393A1 (en) * 2011-04-29 2014-02-20 Shanxi Electric Power Research Institute Method for determining spatial location of conducting wire and aerial earth wire of power transmission line
CN103810329A (zh) * 2014-01-18 2014-05-21 国家电网公司 根据输电线路功率调整输电线路导线相序排列的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
中国电力科学研究院: "《特高压输电技术 交流输电分册》", 30 May 2012, 中国电力出版社 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107784187A (zh) * 2017-11-22 2018-03-09 重庆大学 一种基于边界元法的变电站电磁环境分布研究方法

Similar Documents

Publication Publication Date Title
Acosta et al. Methodology for optimizing the capacity and costs of overhead transmission lines by modifying their bundle geometry
CN106547990B (zh) 特高压直流换流阀塔阀层集成宽频等效电路模型的建模方法
CN104375056B (zh) 基于电压及电流初始行波的变电站电缆出线故障监测方法
CN207541193U (zh) 一种配电网多态复杂故障模拟试验平台
CN110955990B (zh) 一种基于多物理耦合场的变压器绕组瞬态形变量计算方法
CN203772925U (zh) Gis用多输出电子式电压互感器
CN103605836A (zh) 一种高压变电站三维电磁场并行计算方法
CN107846013A (zh) 基于peec法的gis外壳环流和暂态地电位升高建模与分析方法
CN105005665A (zh) 油纸电容式套管末屏的特高频辐射特性仿真方法与系统
CN103474998B (zh) 500kV电缆及架空混合线路过电压优化控制方法
Gao et al. High-impedance arc fault modeling for distribution networks based on dynamic geometry dimension
CN110807243A (zh) 一种计及频变参数的变压器绕组等值电路模型建立方法
CN106203766A (zh) 一种500kV变电站串抗区域工频电场优化方法
CN107679285B (zh) 特高压直流换流阀二端口电路
CN103762066A (zh) 一种gis用多输出电子式电流电压组合互感器
CN102201669A (zh) 一种±660kV阀厅直流金具选型方法
CN105678014B (zh) 一种气体介质的储能特征与放电电压预测方法
Liu et al. The influence of VFTO on SCSR windings with an equivalent PEEC model
CN204028345U (zh) 一种检测电压互感器抵御地电位升高能力的系统
CN106446418A (zh) 一种电缆终端进水缺陷电场分布分析方法及装置
CN203774078U (zh) Gis用多输出电子式电流电压组合互感器
Zhang Optimization of corona ring design used in UHV composite insulator by PSO algorithm
Wang et al. Distribution of electric field and structure optimisation on the surface of a±1100 kV smoothing reactor
CN105514981A (zh) 一种基于建筑信息模型的光伏逆变器接入位置优化方法
CN106777632A (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: 20161207

RJ01 Rejection of invention patent application after publication