一种变电站设备时钟同步的改进方法技术

本发明专利技术公开了一种变电站设备时钟同步的改进方法，步骤一：本地时钟的频偏计算，解得本地时钟频偏η的最大似然估计值为；步骤二：本地时钟的相偏计算，计算N次对时相偏Φk的平均值作为相偏值Φ；步骤三：对本地时钟晶振固有频率进行补偿；步骤四：对本地时钟进行补偿。针对智能变电站网络负荷较重时，网络传输往返延时不对称，传统的精确时钟同步协议(Precision Time Protocol，PTP)单次对时在本地时钟的校正上误差较大，本发明专利技术提供的一种变电站设备时钟同步的改进方法，通过统计学方法在多次对时后计算出本地时钟频偏与相偏误差的最优解进而进行补偿。

【技术实现步骤摘要】
一种变电站设备时钟同步的改进方法
本专利技术涉及一种变电站设备时钟同步的改进方法，属于变电站网络对时

技术介绍
现代化变电站以全站信息数字化、通信平台网络化、信息共享标准化为基本要求，自动完成信息采集、测量、控制、保护、监测和计量等基本功能。但无论测控装置、保护保护装置，还是计量装置的计算处理都要求采样数据应在同一个时间点上采集，以免相位和幅值产生误差。对于过电流保护等保护，因为微机保护的本地晶振时钟的短期稳定性非常高，不会对保护的动作精度造成影响。但对于差动保护和计量，由于微机保护本地晶振时钟并不十分准确，经过长时间的误差累积，会造成跨间隔间不同采样数据的相位误差和幅值误差的逐渐扩大，导致差动保护的误动作和计量的严重误差，因此需要对时技术对本地时钟进行校正。IEEE1588定义了一个能在测量和控制系统中实现高精度时钟同步的协议——精确时钟同步协议(PrecisionTimeProtocol，PTP)。PTP协议的基本原理是主时钟与本地时钟之间进行同步信息包的发送，对信息包的发出时间和接收时间信息进行记录，并且在每一条信息包上打上时间戳标签。PTP协议为了使用和描述时间信息，定义了四种多点传送的信息包，即同步信息包Sync、同步信息包之后的跟随信息包Follow_Up、延时测量信息包Delay_Req和延时测量应答信息包Delay_Resp。同步信息包传递的机制成为“延时—请求响应机制”，如图1所示，主时钟周期性发送包含时钟质量的Sync信息包，紧接着发送Follow_up信息包告知本地时钟上个信息包的实际发送时间Tm1；本地时钟记录Sync信息包的到达时间Ts1，随后在Ts3时刻发送Delay_Req信息包；主时钟记录Delay_Req信息包的到达时间Tm3，并发送信息包Delay_Resp把Tm3告知本地时钟。本地时钟根据这四个时间信息计算出时钟的偏移和传输延迟。假设在时钟的信息通信中认为信息包的往返延迟是相等的，则图1中Tm1、Ts1、Ts3、Tm3四个点的连线是一个等腰梯形。本地时钟便可计算出自身与主时钟之间的信息交换延时Tdelay为Tdelay＝[(Tm3﹣Tm1)﹣(Ts3﹣Ts1)]/2(1.1)本地时钟与主时钟的时钟相位偏移Toffset为Toffset＝[(Ts1﹣Tm1)+(Ts3﹣Tm3)]/2(1.2)本地时钟根据计算出来的相位偏移修改本地时间，从而达到与主时钟同步。PTP协议通过这种“延时—请求响应机制”消除了物理层处理延时的影响，从而进一步提高了对时精度。虽然IEEE1588通过PTP协议消除了上层处理延时与物理层处理延时造成的影响，但它与所有网络对时算法一样存在局限性，PTP协议为计算时钟相偏所进行的四次测量都以对时信息的传输延迟在传输方向上相同为基础，但在实际使用中，不可能绝对满足这个前提。虽然IEEE1588协议引入了透明时钟和边界时钟，但需要交换机支持IEEE1588时钟机制，同时网络传输延时测量的随机性和较低频率，但当网络负载突发性变化时，造成信息包在传输方向上的延时差异变大，因此公式(1.1)与公式(1.2)成立前提的假设——信息包交换的往返延迟相等便不成立，图1中Tm1、Ts1、Ts3、Tm3四个点的连线将是一个不等腰梯形，若简单地通过公式(1.2)对本地时钟相偏进行计算，便会对本地时钟的校正带来较大误差。该方法只对本地时钟的相位偏移进行了补偿，并未对本地时钟的频率偏移进行补偿。而造成相位偏移的根本原因是本地晶振时钟频率并不十分准确，与主时钟频率相比存在误差，经过长时间的误差累积而造成相位偏差。
技术实现思路
目的：为了克服现有技术中存在的不足，本专利技术提供一种变电站设备时钟同步的改进方法。技术方案：为解决上述技术问题，本专利技术采用的技术方案为：一种变电站设备时钟同步的改进方法，包括步骤如下：步骤一：本地时钟的频偏计算，假设本地时钟与主时钟之间由于频偏η的影响，时钟相偏单调递增，即频偏η不变，构建时钟相偏与频偏的对时信息交换的时钟模型；一次对时中对时信息共进行N次交换，其中第k次对时信息交换的T2,k与T3,k均为本地时钟的本地时间，与均为主时钟的本地时间；忽略上层处理延时，物理层处理延时为定值d，信息由主时钟到本地时钟的网络传输延时为随机变量Xk，信息由本地时钟到主时钟的网络传输延时为随机变量Yk；选择T1,1为参考时间，即T1,1为主时钟实际本地时间的零时刻；则与分别为主时钟实际本地时间与的相对时间；表示在时刻的参考时钟相位偏移，其中Φr代表本地时钟的实际相偏；本地时钟接收到第k次对时信息时的相对观测时间为T2,k，发送第k次对时信息时的相对观测时间为T3,k，有：有同理有T3,k＝(T4,k-d-Yk)(1+η)+Φ(2.2)由式(2.1)，用T2,N减去T2,1可得：T2,N-T2,1＝(T1,N-T1,1+XN-X1)(1+η)(2.3)同理，根据式(2.2)可得：T3,N-T3,1＝(T4,N-T4,1+Y1-YN)(1+η)(2.4)将第N个与第1个时间戳的差值分别定义为D1＝T1,N-T1,1，D2＝T2,N-T2,1，D3＝T3,N-T3,1，D4＝T4,N-T4,1，且第N次信息交换与第1次信息交换的网络传输差值定义为DX＝XN-X1与DY＝YN-Y1，则式(2.3)与(2.4)可分别改写为D2＝(D1+DX)(1+η)(2.5)D3＝(D4-DY)(1+η)(2.6)假设网络传输延时和为相互独立但同指数分布的随机变量，指数分布均值为λ，则DX与DY为零均值的拉普拉斯随机变量，方差为2λ2；因此DX与DY的联合概率密度函数为：式(2.5)与(2.6)描述的简化观测模型的似然函数为：令代入式(2.8)得：通过使式(2.8)的似然函数L(η′,λ)取得最大值，则可以得到η′的最大似然估计值，为式中为的从小到大顺序统计量，Ki为与δ(i)对应的D2或D3，当δ(1)＝δ1时，K1＝D2，K2＝D3；当δ(1)＝δ2时，K1＝D3，K2＝D2。解得本地时钟频偏η的最大似然估计值为：步骤二：本地时钟的相偏计算，相偏计算时需要先对本地时钟的观测时间T2,k、T3,k(k＝1…N)进行频率补偿，得到修正后的时间T’2,k、T’3,k：T′2,k＝T2,k-(T1,k+d+Xk)ηMLE(3.1)T′3,k＝T3,k-(T4,k-d-Yk)ηMLE(3.2)由式(3.3)计算出N次对时每次的相偏ΦkΦk＝T′2,k+T′3,k-T1,k-T4,k(3.3)计算N次对时相偏Φk的平均值作为相偏值Φ步骤三：设本地时钟晶振固有频率为f，对其进行补偿后，实际频率f’为f’＝f×(1+ηMLE)(2.12)步骤四：设本地时钟此时输出时间为Ci(t)，对其进行补偿后，实际输出时间C’i(t)为C’i(t)＝Ci(t)+Φ(3.5)有益效果：针对智能变电站网络负荷较重时，网络传输往返延时不对称，传统的精确时钟同步协议(PrecisionTimeProtocol，PTP)单次对时在本地时钟的校正上误差较大，本专利技术提供的一种变电站设备时钟同步的改进方法，通过统计学方法在多次对时后计算出本地时钟频偏与相偏误差的最优解进而进行补偿。附图说明图1为延时—请求响应机制本文档来自技高网...

【技术保护点】
一种变电站设备时钟同步的改进方法，其特征在于：包括步骤如下：步骤一：本地时钟的频偏计算，假设本地时钟与主时钟之间由于频偏η的影响，时钟相偏单调递增，即频偏η不变，构建时钟相偏与频偏的对时信息交换的时钟模型；一次对时中对时信息共进行N次交换，其中第k次对时信息交换的T2,k与T3,k均为本地时钟的本地时间，与均为主时钟的本地时间；忽略上层处理延时，物理层处理延时为定值d，信息由主时钟到本地时钟的网络传输延时为随机变量Xk，信息由本地时钟到主时钟的网络传输延时为随机变量Yk；选择T1,1为参考时间，即T1,1为主时钟实际本地时间的零时刻；则与分别为主时钟实际本地时间与的相对时间；表示在时刻的参考时钟相位偏移，其中Φr代表本地时钟的实际相偏；本地时钟接收到第k次对时信息时的相对观测时间为T2,k，发送第k次对时信息时的相对观测时间为T3,k，有：有同理有T3,k＝(T4,k‑d‑Yk)(1+η)+φ(2.2)由式(2.1)，用T2,N减去T2,1可得：T2,N‑T2,1＝(T1,N‑T1,1+XN‑X1)(1+η)   (2.3)同理，根据式(2.2)可得：T3,N‑T3,1＝(T4,N‑T4,1+Y1‑YN)(1+η)   (2.4)将第N个与第1个时间戳的差值分别定义为D1＝T1,N‑T1,1，D2＝T2,N‑T2,1，D3＝T3,N‑T3,1，D4＝T4,N‑T4,1，且第N次信息交换与第1次信息交换的网络传输差值定义为DX＝XN‑X1与DY＝YN‑Y1，则式(2.3)与(2.4)可分别改写为D2＝(D1+DX)(1+η)   (2.5)D3＝(D4‑DY)(1+η)   (2.6)假设网络传输延时和为相互独立但同指数分布的随机变量，指数分布均值为λ，则DX与DY为零均值的拉普拉斯随机变量，方差为2λ2；因此DX与DY的联合概率密度函数为：式(2.5)与(2.6)描述的简化观测模型的似然函数为：令代入式(2.8)得：通过使式(2.8)的似然函数L(η′,λ)取得最大值，则可以得到η′的最大似然估计值，为式中为的从小到大顺序统计量， Ki为与δ(i)对应的D2或D3，当δ(1)＝δ1时，K1＝D2，K2＝D3；当δ(1)＝δ2时，K1＝D3，K2＝D2。解得本地时钟频偏η的最大似然估计值为：步骤二：本地时钟的相偏计算，相偏计算时需要先对本地时钟的观测时间T2,k、T3,k(k＝1…N)进行频率补偿，得到修正后的时间T’2,k、T’3,k：T′2,k＝T2,k‑(T1,k+d+Xk)ηMLE   (3.1)T′3,k＝T3,k‑(T4,k‑d‑Yk)ηMLE   (3.2)由式(3.3)计算出N次对时每次的相偏ΦkΦk＝T′2,k+T′3,k‑T1,k‑T4,k   (3.3) 计算N次对时相偏Φk的平均值作为相偏值Φ步骤三：设本地时钟晶振固有频率为f，对其进行补偿后，实际频率f’为f’＝f×(1+ηMLE)   (2.12)步骤四：设本地时钟此时输出时间为Ci(t)，对其进行补偿后，实际输出时间C’i(t)为C’i(t)＝Ci(t)+Φ   (3.5)。...

【技术特征摘要】
1.一种变电站设备时钟同步的改进方法，其特征在于：包括步骤如下：步骤一：本地时钟的频偏计算，假设本地时钟与主时钟之间由于频偏η的影响，时钟相偏单调递增，即频偏η不变，构建时钟相偏与频偏的对时信息交换的时钟模型；一次对时中对时信息共进行N次交换，其中第k次对时信息交换的T2,k与T3,k均为本地时钟的本地时间，与均为主时钟的本地时间；忽略上层处理延时，物理层处理延时为定值d，信息由主时钟到本地时钟的网络传输延时为随机变量Xk，信息由本地时钟到主时钟的网络传输延时为随机变量Yk；选择T1,1为参考时间，即T1,1为主时钟实际本地时间的零时刻；则与分别为主时钟实际本地时间与的相对时间；表示在时刻的参考时钟相位偏移，其中Φr代表本地时钟的实际相偏；本地时钟接收到第k次对时信息时的相对观测时间为T2,k，发送第k次对时信息时的相对观测时间为T3,k，有：有同理有T3,k＝(T4,k-d-Yk)(1+η)+Φ(2.2)由式(2.1)，用T2,N减去T2,1可得：T2,N-T2,1＝(T1,N-T1,1+XN-X1)(1+η)(2.3)同理，根据式(2.2)可得：T3,N-T3,1＝(T4,N-T4,1+Y1-YN)(1+η)(2.4)将第N个与第1个时间戳的差值分别定义为D1＝T1,N-T1,1，D2＝T2,N-T2,1，D3＝T3,N-T3,1，D4＝T4,N-T4,1，且第N次信息交换与第1次信息交换的网络传输差值定义为DX＝XN-X1与DY＝YN-Y1，则式(2.3)与(2.4)可分别改写为D2＝(D1+DX)(1+η)(2.5)D3＝(D4-DY)(1+η)(2.6)假设网络传输延时和为相互独立但同指数分布的随机变量，指数分布均值为λ，则DX与DY为零均值的拉普拉斯随机变量，方差为2λ2；因此DX与DY的联合概率密度函数为：式(2.5)与(2.6)描述的简化观测模型的似然函数为：令...

【专利技术属性】
技术研发人员：陈昊，朱超，张兆君，黄祖荣，
申请(专利权)人：江苏省电力公司检修分公司，国家电网公司，国网江苏省电力公司，
类型：发明
国别省市：江苏;32