CN102967780B - 一种变电站智能告警建模方法 - Google Patents

一种变电站智能告警建模方法 Download PDF

Info

Publication number
CN102967780B
CN102967780B CN201210475190.3A CN201210475190A CN102967780B CN 102967780 B CN102967780 B CN 102967780B CN 201210475190 A CN201210475190 A CN 201210475190A CN 102967780 B CN102967780 B CN 102967780B
Authority
CN
China
Prior art keywords
fault
model
substation
electric network
attribute
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
CN201210475190.3A
Other languages
English (en)
Other versions
CN102967780A (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.)
Nari Technology Co Ltd
Original Assignee
Nari Technology 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 Nari Technology Co Ltd filed Critical Nari Technology Co Ltd
Priority to CN201210475190.3A priority Critical patent/CN102967780B/zh
Publication of CN102967780A publication Critical patent/CN102967780A/zh
Application granted granted Critical
Publication of CN102967780B publication Critical patent/CN102967780B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)

Abstract

本发明提供了一种变电站智能告警建模方法,包括:针对变电站继电保护动作元件,建立保护动作元件推理模型;建立变电站的一次设备、二次设备、量测数据、状态数据、保护信息之间的映射关系及继电保护装置的配置信息;通过网络拓扑获得一次设备的接线方式以及与相邻间隔的连接关系;基于专家***知识数据库,针对各类故障的条件、征兆、过程特征和结果信息建立故障推理模型;采用基于XML语言的故障描述语言描述所述保护动作元件推理模型和故障推理模型。本发明针对变电站的拓扑结构、元件参数和故障特点,建立了故障推理模型,具有实现灵活、扩展方便、适应性强等优点,因此基于故障推理模型的进行故障推理具有较高的可实施性。

Description

一种变电站智能告警建模方法
技术领域
本发明属于变电站自动化技术领域,特别涉及一种变电站智能告警建模方法。
背景技术
目前,国内主流厂家大多采用基于规则的专家***来实现变电站故障分析,但是这种方式存在着***维护困难、容错能力较差等缺陷,工程化程度低,不利于推广应用。如果可以针对基于规则的专家***的缺陷,通过对站内各种故障以及专家***的知识库进行分析、总结和提炼,将其抽象为故障推理模型,并基于这些故障推理模型对故障进行分析,就能够有效的解决基于规则的专家***存在的问题。因此,在变电站自动化领域中亟待出现一种实现故障推理模型的智能告警建模方法。
发明内容
变电站故障主要包括电网故障和二次设备故障两类,对于电网故障而言,虽然其过程非常复杂,用于故障判断的量测数据、状态数据信息量大,故障特征及过程状态多,但是经过继电保护装置处理后上送的数据却相对简洁明了,这使得基于继电保护装置上送的信息,再辅以测控装置的SOE以及一次设备在线监测***的监测结果信息,对某种特定故障进行描述,进而生成故障模型成为一种现实的选择。而对于二次设备故障而言,目前其判据主要基于二次设备上传的自检信号以及通信状态,虽然存在着设备种类多、型号多、生产厂家多、特性各异,而且上送的自检信号缺乏统一的标准等不利条件,但是对具体故障而言,故障特征及判据都较为明确,其故障模型与电网故障模型相比较为简单,更容易建模。
根据变电站故障的上述特点,本发明提供了一种变电站智能告警建模方法,用于生成变电站故障推理模型,从而依靠建立的故障推理模型实现变电站智能告警,具有实现灵活、扩展方便、适应性强的优点。
所述变电站智能告警建模方法,其特征在于,包括以下步骤:
步骤1,针对变电站继电保护动作元件,建立保护动作元件推理模型;
步骤2,建立变电站的一次设备、二次设备、量测数据、状态数据、保护信息之间的映射关系及继电保护装置的配置信息;
步骤3,通过网络拓扑获得一次设备的接线方式以及与相邻间隔的连接关系;
步骤4,基于专家***知识数据库,针对各类故障的条件、征兆、过程特征和结果信息建立故障推理模型;
步骤5,采用基于XML语言的故障描述语言描述所述保护动作元件推理模型和故障推理模型。
优选地,所述步骤1中针对变电站继电保护动作元件的电压等级、接线方式、保护对象等条件,针对所述继电保护动作元件的动作原因、动作特征、动作影响和动作后处理等参数建立所述保护动作元件推理模型。
优选地,所述步骤2中,在常规变电站中从变电站监控***中提取相关信息进行配置,在基于IEC61850标准变电站中从IEC61850模型提取相关信息进行配置。
优选地,所述步骤4包括:针对单电气间隔电网故障的条件、特征,基于所述保护动作元件推理模型和专家***知识数据库建立单间隔电网故障模型;针对跨电气间隔电网故障,在单间隔电网故障模型基础上,结合网络拓扑信息建立跨间隔电网故障模型;针对变电站二次设备故障,建立二次设备故障模型;结合所述单间隔电网故障模型、跨间隔电网故障模型和二次设备故障模型,形成具有二次设备故障辅助判断的电网故障模型。
优选地,在步骤5中,所述故障描述语言采用分段式字符串的变电站对象及对象属性命名体系。进一步优选地,所述命名体系用于表示在所述故障推理模型中定义所涉及的设备对象、设备对象属性以及设备对象属性的属性,各分段间以“.”分隔,完整的字符结构为“厂站.间隔.设备.设备属性.属性的属性”。
优选地,所述步骤5中,所述故障描述语言定义了用于故障分析的约束条件。
可见,本发明所述的变电站智能告警建模方法针对变电站的拓扑结构、元件参数和故障特点,建立了故障推理模型,具有实现灵活、扩展方便、适应性强等优点,因此基于故障推理模型的进行故障推理具有较高的可实施性。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明实施例中保护动作元件推理建模示意图。
图2是本发明实施例中通过IEC61850模型获得变电站一次设备与二次设备信息映射关系示意图。
图3为本发明实施例的分段式字符串命名体系结构图。
图4为本发明实施例的故障推理模型生成示意图。
图5为220kV线路典型接线图。
图6为图5所示220kV线路某瞬时故障推理模型采用故障逻辑描述语言描述的示例。
具体实施方式
为了使本技术领域的人员更好地理解本发明的技术方案,并使本发明的上述目的、特征和优点能够更加明显易懂,下面结合实施例及实施例附图对本发明作进一步详细的说明。
变电站故障主要包括电网故障和二次设备故障两类,本发明的变电站智能告警建模方法针对以上两种故障分别建立相应的故障推理模型。其中,对于电网故障而言,基于继电保护装置上送的信息,再辅以测控装置的SOE以及一次设备在线监测***的监测结果信息,对某种特定故障进行描述,进而生成故障模型。而对于二次设备故障而言,主要基于二次设备上传的自检信号以及通信状态,根据故障特征及判据建立相应的模型。
本发明所述的智能告警建模方法包括以下步骤:
步骤1,针对变电站继电保护动作元件,建立保护动作元件推理模型。
对于电网故障而言,虽然其过程非常复杂,用于故障判断的量测数据、状态数据信息量大,故障特征及过程状态多,但是经过继电保护装置处理后上送的数据却相对简洁明了,这使得基于继电保护装置上送的信息,再辅以测控装置的SOE以及一次设备在线监测***的监测结果信息,对某种特定故障进行描述,进而生成故障模型成为一种现实的选择。
如图1所示为保护动作元件推理建模示意图,如图所示,获得变电站继电保护动作元件的电压等级、接线方式、保护对象等条件,针对所述继电保护动作元件的动作原因、动作特征、动作影响和动作后处理等参数建立所述保护动作元件推理模型。
步骤2,建立变电站的一次设备、二次设备、量测数据、状态数据、保护信息之间的映射关系及继电保护装置的配置信息。
建立故障推理模型所需的变电站一次设备、二次设备、各类量测数据、状态数据、保护信息之间的映射关系与继电保护装置的配置信息,在常规变电站中,可从变电站监控***中提取相关信息进行配置,在基于IEC61850标准变电站中,可从变电站的IEC61850模型提取相关信息进行配置。如图2所示,为通过IEC61850模型获得变电站一次设备与二次设备信息映射关系示意图,包括从变电站IEC61850模型的变电站***配置描述文件(SCD)中读取测控IED模型和保护IED模型,并且从变电站规范模型文件SSD中读取变电站一次设备模型,从而建立变电站故障分析建模所需要的映射关系,包括量测数据与一次设备、状态数据与一次设备、保护信息与一次设备之间的映射关系,以及保护IED的配置。
步骤3,通过网络拓扑获得一次设备的接线方式以及与相邻间隔的连接关系。
步骤4,基于专家***知识数据库,针对各类故障的条件、征兆、过程特征和结果信息建立故障推理模型。在本步骤中,对变电站内各种故障以及专家***知识库进行分析、总结,提炼出故障条件、故障征兆、故障过程特征以及故障结果,将其抽象为故障推理模型图4所示为变电站故障推理模型生成示意图。针对单电气间隔电网故障的条件、特征,基于步骤1所述的所述保护动作元件推理模型和专家***知识数据库,建立反映变电站单电气间隔电网故障的故障条件和故障特征的单间隔电网故障模型;针对跨电气间隔电网故障,在单间隔电网故障模型基础上,结合步骤3所述的网络拓扑信息建立跨间隔电网故障模型;针对变电站二次设备故障,建立二次设备故障模型;结合所述单间隔电网故障模型、跨间隔电网故障模型和二次设备故障模型,形成具有二次设备故障辅助判断的电网故障模型。
步骤5,采用基于XML语言的故障描述语言描述所述保护动作元件推理模型和故障推理模型。在故障描述语言中,定义了用于故障分析的电压等级、相别、时间等约束条件。为了便于使用故障描述语言描述变电站故障推理模型的故障条件、故障特征等属性,建立了一套分段式字符串的变电站对象及对象属性命名体系,该命名体系参照了IEC61850标准中Reference的定义。如图3所示,本发明采用的分段式字符串变电站设备及设备属性命名体系结构,不同分段直接以“.”分隔,完整的字符结构为“厂站.间隔.设备.设备属性.属性的属性”,例如“SFac01.GLine01.Br5011.Ia.Quality”。其中需要特别说明的是:
a) “SXXX”字段用于定义厂站信息。该字段以“S”开头,意为Station。当需要建立多厂站故障模型时,需要定义该段内容,采用“S+厂站关键字”形式定义;单厂站模式下填写缺省值“SDefault”,意为当前厂站。
b) “GXXX”字段用于定义指定厂站中的电气间隔(设备组),该字段以“G”开头,意为Group。当模型中需要获得相邻电气间隔信息时,在该字段指定相关的电气间隔,采用“G+电气间隔关键字”形式定义;本间隔内属性填写缺省值“GDefault”,意为本电气间隔。
c) 第三字段用于定义指定电气间隔中的对象(设备)。一个电气间隔包括多种对象(设备),例如断路器、隔离开关、接地刀闸、测控IED、保护IED等,断路器以“Br”开头,采用“Br + 断路器关键字”形式定义;隔离开关以“Sw”开头,采用“Sw + 隔离开关关键字”形式定义;接地刀闸以“ESw”开头,采用“ESw + 接地刀闸关键字”形式定义;测控IED以“Mcd”开头,采用“Mcd + 测控IED关键字”形式定义;保护IED以“Rly”开头,采用“Rly + 保护IED关键字”形式定义,其它设备定义,以此类推。
d) 第四字段用于定义前三段指定对象的属性,例如电压,电流,位置等,不同属性采用不同的字符串定义,例如Ia,Uab,Pos等。
e) 可以扩展使用第五字段,用于定义属性的属性,例如Ia的品质,Uab的品质等。
下面以具体实例说明采用故障逻辑描述语言描述的故障推理模型的示例。图5为220kV线路典型接线图。如图6所示,为图5所示220kV线路某瞬时故障推理模型采用故障逻辑描述语言描述示例。线路发生瞬时故障的典型特征为:线路发生故障时,线路保护装置动作出口跳开断路器,然后线路保护重合闸动作重合断路器,重合闸成功。该故障的约束条件包括:从故障发生到故障切除的时间(通常为7-8秒)、电压等级、开关是否分相等等,基于XML的灵活性,这些约束条件可以通过扩展关键字来定义。模型中的故障条件、故障特征等属性,采用分段式字符串变电站设备及设备属性命名体系。故障模型中一些特征判断需要进行逻辑运算,采用脚本来实现。
可见,本发明所述的变电站智能告警建模方法针对变电站的拓扑结构、元件参数和故障特点,建立了故障推理模型,具有实现灵活、扩展方便、适应性强等优点,因此基于故障推理模型的进行故障推理具有较高的可实施性。
以上所述,仅为本发明的具体实施方式,本发明还可以应用在其它设备中;以上描述中的尺寸和数量均仅为参考性的,本领域技术人员可根据实际需要选择适当的应用尺寸,而不脱离本发明的范围。本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求所界定的保护范围为准。

Claims (1)

1.一种变电站智能告警建模方法,其特征在于,包括以下步骤:
步骤1,针对变电站继电保护动作元件,建立保护动作元件推理模型;根据电压等级、接线方式和保护对象条件,结合所述继电保护动作元件的动作原因、动作特征、动作影响和动作后处理参数建立所述保护动作元件推理模型;
步骤2,建立变电站的一次设备、二次设备、量测数据、状态数据、保护信息之间的映射关系及获取继电保护装置的配置信息;在常规变电站中从变电站监控***中提取相关信息进行配置,在基于IEC61850标准变电站中从IEC61850模型提取相关信息进行配置;
步骤3,通过网络拓扑获得一次设备的接线方式以及与相邻间隔的连接关系;
步骤4,基于专家***知识数据库,针对各类故障的条件、征兆、过程特征和结果信息建立故障推理模型;针对单电气间隔电网故障的条件、特征,基于所述保护动作元件推理模型和专家***知识数据库建立单间隔电网故障模型;针对跨电气间隔电网故障,在单间隔电网故障模型基础上,结合网络拓扑信息建立跨间隔电网故障模型;针对变电站二次设备故障,建立二次设备故障模型;结合所述单间隔电网故障模型、跨间隔电网故障模型和二次设备故障模型,形成具有二次设备故障辅助判断的电网故障模型;
步骤5,采用基于XML语言的故障描述语言描述所述保护动作元件推理模型和故障推理模型;所述故障描述语言采用形式为分段式字符串的变电站对象及对象属性命名体系;所述命名体系用于表示在所述故障推理模型中定义所涉及的设备对象、设备对象属性以及设备对象属性的属性,各分段间以“.”分隔,完整的字符结构为“厂站.间隔.设备.设备属性.属性的属性”;所述故障描述语言定义了用于故障分析的约束条件。
CN201210475190.3A 2012-11-21 2012-11-21 一种变电站智能告警建模方法 Active CN102967780B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210475190.3A CN102967780B (zh) 2012-11-21 2012-11-21 一种变电站智能告警建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210475190.3A CN102967780B (zh) 2012-11-21 2012-11-21 一种变电站智能告警建模方法

Publications (2)

Publication Number Publication Date
CN102967780A CN102967780A (zh) 2013-03-13
CN102967780B true CN102967780B (zh) 2015-04-08

Family

ID=47798060

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210475190.3A Active CN102967780B (zh) 2012-11-21 2012-11-21 一种变电站智能告警建模方法

Country Status (1)

Country Link
CN (1) CN102967780B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103278731B (zh) * 2013-05-20 2016-01-20 国家电网公司 一种基于二次回路信息的故障诊断方法及***
CN103296757B (zh) * 2013-05-21 2014-09-03 国家电网公司 基于多参量辨识的智能变电站二次***故障判断方法
CN103324845B (zh) * 2013-06-13 2016-12-07 国电南瑞科技股份有限公司 一种基于不确定性推理的智能变电站开关状态估计方法
CN103632308B (zh) * 2013-11-01 2017-06-16 广东电网公司电力调度控制中心 基于专家***的涉网电厂继电保护定值风险评估***
CN103679354A (zh) * 2013-12-02 2014-03-26 南京磐能电力科技股份有限公司 基于多层故障推理库和智能推理机的智能告警方法
CN104375020B (zh) * 2014-10-08 2015-11-25 河海大学 一种分布式在线电网故障诊断方法
CN104699801B (zh) * 2015-03-19 2018-04-27 国家电网公司 一种二次设备关联配置的方法及装置
CN105139158A (zh) * 2015-10-10 2015-12-09 国家电网公司 电网异常信息智能告警和辅助决策方法
CN106685704A (zh) * 2016-12-13 2017-05-17 云南电网有限责任公司电力科学研究院 变电站二次设备故障的排查方法及装置
CN106680634B (zh) * 2016-12-30 2019-08-02 国家电网公司 一种智能变电站继电保护安全措施的智能预警方法
CN107733082B (zh) * 2017-10-11 2020-02-28 国网湖南省电力公司 一种继电保护设备装置异常信息收集方法及***
CN108596459B (zh) * 2018-04-13 2021-08-10 广东电网有限责任公司广州供电局 变电站风险告警方法和***
CN110955694A (zh) * 2018-09-26 2020-04-03 中国电力科学研究院有限公司 一种变电站智能告警通用规则库的描述方法及***
CN109683571B (zh) * 2018-12-28 2020-04-21 南京南瑞继保工程技术有限公司 一种水电站监控***的建模方法和装置
CN109633381A (zh) * 2019-02-21 2019-04-16 贵州电网有限责任公司 一种电网故障诊断智能分析方法
CN110674193A (zh) * 2019-10-25 2020-01-10 国网四川省电力公司技能培训中心 一种智能变电站继电保护故障信息建模方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101989766A (zh) * 2009-08-03 2011-03-23 南京南瑞继保电气有限公司 基于保护动作链的电力***故障诊断的方法
CN102035260A (zh) * 2010-12-17 2011-04-27 国电南瑞科技股份有限公司 智能变电站监控***智能告警及故障分析实现方法
CN102123045A (zh) * 2011-01-18 2011-07-13 河南龙源许继科技发展股份有限公司 基于iec61850的高压设备在线监测的建模方法
CN102497027A (zh) * 2011-12-20 2012-06-13 浙江省电力公司 数字化故障录波装置自动建模的方法
CN102497024A (zh) * 2011-12-16 2012-06-13 广东电网公司茂名供电局 基于整数规划的智能告警***
CN102708472A (zh) * 2012-06-20 2012-10-03 江西省电力公司超高压分公司 组合式电网告警信息智能化处理方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101266953B1 (ko) * 2011-02-23 2013-05-30 한전케이디엔주식회사 분산전원을 이용한 건전구간 전원공급 시스템 및 방법

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101989766A (zh) * 2009-08-03 2011-03-23 南京南瑞继保电气有限公司 基于保护动作链的电力***故障诊断的方法
CN102035260A (zh) * 2010-12-17 2011-04-27 国电南瑞科技股份有限公司 智能变电站监控***智能告警及故障分析实现方法
CN102123045A (zh) * 2011-01-18 2011-07-13 河南龙源许继科技发展股份有限公司 基于iec61850的高压设备在线监测的建模方法
CN102497024A (zh) * 2011-12-16 2012-06-13 广东电网公司茂名供电局 基于整数规划的智能告警***
CN102497027A (zh) * 2011-12-20 2012-06-13 浙江省电力公司 数字化故障录波装置自动建模的方法
CN102708472A (zh) * 2012-06-20 2012-10-03 江西省电力公司超高压分公司 组合式电网告警信息智能化处理方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
应用于电网故障诊断的统一信息模型;高湛军等;《电力***自动化》;20100510;第34卷(第09期);第1-3节 *

Also Published As

Publication number Publication date
CN102967780A (zh) 2013-03-13

Similar Documents

Publication Publication Date Title
CN102967780B (zh) 一种变电站智能告警建模方法
CN105429291B (zh) 智能变电站一次接线图和二次***逻辑连接图的关联方法
Higgins et al. Distributed power system automation with IEC 61850, IEC 61499, and intelligent control
CN101551440B (zh) 一种发电机变压器组故障录波分析装置
CN102937685B (zh) 一种变电站用基于实时仿真技术的一体化测试平台
CN103715767B (zh) 智能电网信息集成与展示平台
CN102184209A (zh) 一种基于cim电网模型接口的仿真数据接入方法
CN111478324B (zh) 一种电站告警分析***
CN108241771A (zh) 基于电磁暂态仿真平台的二次设备建模方法
CN108982996B (zh) 一种实景替代式的特高压直流保护现场测试装置
CN103995172B (zh) 变电站gis母线负荷电流在线监测方法
CN103971292B (zh) 基于实时站内图映射的故障准确辨识方法
CN103425847B (zh) 基于电网调度防误***的线路预判改造***及其工作方法
WO2020171814A1 (en) Event selection for power grid disturbance
CN104104151A (zh) 配电终端设备信息的接入方法
CN103455678B (zh) 基于标准化虚端子的二次虚回路图表协同设计方法
CN115169506A (zh) 一种供变电关键设备故障快速诊断方法及***
CN106569096B (zh) 一种配电网单相故障的在线定位方法
CN109149534B (zh) 基于dts虚拟开关的快速诊断电网模型拓扑故障的方法
CN105046011A (zh) 基于分布式电网拓扑计算的电气设备状态快速分析方法
CN105913337B (zh) 一种psasp数据模型转换为rscad数据模型的方法
Kezunovic et al. Merging PMU, operational, and non-operational data for interpreting alarms, locating faults and preventing cascades
CN110489729B (zh) D5000-matpower的电网多灾种耦合连锁故障模型的自动转换方法及***
CN111899127B (zh) 基于层次化的故障波及树建模方法、模型及其应用方法
CN106709158B (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
C14 Grant of patent or utility model
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Xu Zufeng

Inventor after: Zhang Haibin

Inventor after: Zhang Jian

Inventor after: Pei Yulong

Inventor after: Wang Haifeng

Inventor after: Dai Xingya

Inventor after: Zheng Yuping

Inventor before: Xu Zufeng

Inventor before: Zhang Haibin

Inventor before: Zhang Jian

Inventor before: Pei Yulong

Inventor before: Wang Haifeng

Inventor before: Dai Xingya

COR Change of bibliographic data
CB03 Change of inventor or designer information

Inventor after: Xu Zufeng

Inventor after: Wu Xiaozhong

Inventor after: Sun Chaowu

Inventor after: Zheng Yuping

Inventor after: Zhang Haibin

Inventor after: Zhang Jian

Inventor after: Pei Yulong

Inventor after: Wang Haifeng

Inventor after: Dai Xingya

Inventor after: Huang Lai

Inventor before: Xu Zufeng

Inventor before: Zhang Haibin

Inventor before: Zhang Jian

Inventor before: Pei Yulong

Inventor before: Wang Haifeng

Inventor before: Dai Xingya

Inventor before: Zheng Yuping

COR Change of bibliographic data