CN111447030A - 地震采集器的网络授时精度控制方法、系统及终端设备 - Google Patents

地震采集器的网络授时精度控制方法、系统及终端设备 Download PDF

Info

Publication number
CN111447030A
CN111447030A CN202010244198.3A CN202010244198A CN111447030A CN 111447030 A CN111447030 A CN 111447030A CN 202010244198 A CN202010244198 A CN 202010244198A CN 111447030 A CN111447030 A CN 111447030A
Authority
CN
China
Prior art keywords
value
data packet
client
time
server
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.)
Granted
Application number
CN202010244198.3A
Other languages
English (en)
Other versions
CN111447030B (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.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
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 Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN202010244198.3A priority Critical patent/CN111447030B/zh
Publication of CN111447030A publication Critical patent/CN111447030A/zh
Application granted granted Critical
Publication of CN111447030B publication Critical patent/CN111447030B/zh
Priority to PCT/CN2021/077592 priority patent/WO2021196923A1/zh
Priority to US17/475,369 priority patent/US11716162B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J3/00Time-division multiplex systems
    • H04J3/02Details
    • H04J3/06Synchronising arrangements
    • H04J3/0635Clock or time synchronisation in a network
    • H04J3/0682Clock or time synchronisation in a network by delay compensation, e.g. by compensation of propagation delay or variations thereof, by ranging
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J3/00Time-division multiplex systems
    • H04J3/02Details
    • H04J3/06Synchronising arrangements
    • H04J3/0635Clock or time synchronisation in a network
    • H04J3/0638Clock or time synchronisation among nodes; Internode synchronisation
    • H04J3/0658Clock or time synchronisation among packet nodes
    • H04J3/0661Clock or time synchronisation among packet nodes using timestamps
    • H04J3/0667Bidirectional timestamps, e.g. NTP or PTP for compensation of clock drift and for compensation of propagation delays
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L7/00Arrangements for synchronising receiver with transmitter
    • H04L7/0016Arrangements for synchronising receiver with transmitter correction of synchronization errors
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W56/00Synchronisation arrangements
    • H04W56/004Synchronisation arrangements compensating for timing error of reception due to propagation delay
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/24Recording seismic data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2200/00Details of seismic or acoustic prospecting or detecting in general
    • G01V2200/10Miscellaneous details
    • G01V2200/12Clock synchronization-related issues

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)
  • Computer And Data Communications (AREA)

Abstract

本发明适用于地震数据采集技术领域,提供了一种地震采集器的网络授时精度控制方法、系统及终端设备,包括:采用中断模式使传递数据包;计算最优网络延迟;根据最优网络延迟对网络授时过程中的传递误差进行校正后,由服务器的物理层接收所述数据包,通过中断模式将数据由服务器的物理层发送到服务器的应用层为数据包进行授时。本申请采用中断模式使数据包在物理层和应用层之间传递,实现了数据包时间戳精确加盖,并且由于此算法在客户端和服务器都采用对称算法模式,其产生的执行误差可由于对称性大大减少。针对客户端和服务器之间数据包传递产生的误差通过计算最优网络延时对其进行调整,大大减小授时误差,提高网络授时精度。

Description

地震采集器的网络授时精度控制方法、系统及终端设备
技术领域
本发明涉及地震数据采集领域,尤其涉及一种地震采集器的网络授时精度控制方法、系统、终端设备及计算机可读存储介质。
背景技术
在地球物理领域,特别是地震数据采集领域,其地震采集器的采集的地震数据的时间信息是地震信号非常重要的信息之一,具体到地震数据采集器,其自身的时间服务也是其重要指标之一,现在所有的地震数据采集器其时间授时服务来源主要来自于两个方面,一个是基于GPS的时间服务,具体为从串口,I2C或者SPI等通信接口接收固定格式的报文信息,其信息内包含具体的年月日,时分秒信息,另外配合PPS秒脉冲实现时间的精确授时;另外一种是基于网络的时间授时机制,其具体分为两种,一种是基于通用网络的NTP网络授时协议,其授时机制工作在网络应用层,可以实现在通用的网络硬件下实现时间同步,缺点是授时精度低,其授时精度在ms级;另外一种基于网络的时间同步协议是PTP网络授时协议,其工作在特定的网络链路下,需要特殊的硬件配合才能实现网络授时,其优点是授时精度高,其授时精度可精确到微妙级,缺点是实施成本高,包括必要的硬件成本和软件运行成本,且需要特殊的网络硬件支持,因此在通用的网络硬件条件下不利于实施。
对于地震数据采集领域,常规的地震数据采集在处理地球物理天然地震信号的情况下要求的采样率在500以下,且对于主动源地震采集,其数据采样率要求也在10K以下,因此对于其特殊的应用领域,其地震数据采集器的时钟精度需要控制在微秒到十分之一毫秒之间,且对于基于实时数据传输到地震采集器而言,由于其需要将采集到数据通过网络传输到固定的服务器或者台网中,因此其必然标配正常的网络通信模组,可实现网络应用层上的软件嵌入,同时由于地震数据采集器其部分需要布设到矿井,大楼内部等特殊环境,其环境内不可能接受到GPS信号。
有鉴于此,需要提出一种新的基于NTP的网络协议的授时控制方法,以实现时间的精确授时。
发明内容
鉴于此,本发明实施例提供了一种地震采集器的网络授时精度控制方法、系统及终端设备,该方法可应用于具有网络传输硬件的数据采集器中,可在不改变硬件条件下通过软件算法实现时间授时,并大幅的提高授时精度,且为达到最大兼容性,其算法可在NTPV0版本上运行。
本发明实施例的第一方面提供了一种地震采集器的网络授时精度控制方法,所述网络授时精度控制方法包括:
对客户端接收到的数据包,采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层;
所述客户端向服务器发型所述数据包时,预先设置网络环境约束值;
连续检测不同时刻所述网络环境约束值次对时过程中的GAP值和DELAY 值,其中,GAP值为所述客户端和所述服务器之间的时钟偏移量,DELAY值为数据包在所述客户端和所述服务器之间往返过程中的时间延迟;
对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟;
根据所述最优网络延迟对网络授时过程中的传递误差进行校正后,由所述服务器的物理层接收所述数据包;
通过所述中断模式将所述数据由服务器的物理层发送到服务器的应用层为数据包进行授时。
可选地,所述采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层,包括:
预先设置客户端的应用层的发送数据包的模式及物理层接收数据包的模式为中断模式,设置所述客户端的物理层优先级为最高,应用层中断的优先级为所述客户端的物理层优先级减1;
T1时刻客户端将NTP数据包打包,发送时间字段填充为打包时刻所获取的所述客户端自身时间CT1加1ms,等待中断发生时由应用层将数据包发出;
物理层检测到所述数据包到来时,接收中断,并将此刻获取的所述客户端的自身时间进行存储,对数据包进行标志位后等待处理。
可选地,所述对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟,包括:
首次循环选取网络环境约束值中的DELAY值最小值作为DELAY0,并将 DELAY0以及其对应的GAP0存储到样本库中;
持续统计和测量网络环境约束值中次对时过程中的GAP值和DELAY值,选取所统计的网络环境约束值中的最小值为DELAY1和GAP1;
若DELAY1小于DELAY0,且GAP1在GAP0±(DELAY1-DELAY0)/收敛系数范围内,则所述DELAY1和GAP1更新为新的GPA0和DELAY0,标记此值可用,样本为有效样本;
若有效样本个数大于拟合库样本最大数,则对有效样本中的GAP值和 DELAY值进行卡尔曼滤波,以得到最优解的DELAYD值和GAPD值;
将所得到的最优解的DELAYD值和GAPD值作为最优网络延迟。
可选地,所述网络授时精度控制方法还包括:
若DELAY1不小于DELAY0,或GAP1不在GAP0±(DELAY1-DELAY0)/ 收敛系数范围内,则舍弃此次统计的DELAY值和GAP值。
可选地,所述拟合库样本最大数为64。
本发明实施例的第二方面提供了一种地震采集器的网络授时精度控制系统,所述网络授时精度控制系统包括:
数据处理模块,用于对客户端接收到的数据包,采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层;
最优网络延迟计算模块,用于在所述客户端向服务器发型所述数据包时,预先设置网络环境约束值;连续检测不同时刻所述网络环境约束值次对时过程中的GAP值和DELAY值,其中,GAP值为所述客户端和所述服务器之间的时钟偏移量,DELAY值为数据包在所述客户端和所述服务器之间往返过程中的时间延迟;对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟;
校正授时模块,用于根据所述最优网络延迟对网络授时过程中的传递误差进行校正后,由所述服务器的物理层接收所述数据包;并通过所述中断模式将所述数据由服务器的物理层发送到服务器的应用层为数据包进行授时。
可选地,所述数据处理模块具体用于,
预先设置客户端的应用层的发送数据包的模式及物理层接收数据包的模式为中断模式,设置所述客户端的物理层优先级为最高,应用层中断的优先级为所述客户端的物理层优先级减1;
T1时刻客户端将NTP数据包打包,发送时间字段填充为打包时刻所获取的所述客户端自身时间CT1加1ms,等待中断发生时由应用层将数据包发出;
物理层检测到所述数据包到来时,接收中断,并将此刻获取的所述客户端的自身时间进行存储,对数据包进行标志位后等待处理。
可选地,所述最优网络延迟计算模块在对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟时,具体用于:
首次循环选取网络环境约束值中的DELAY值最小值作为DELAY0,并将 DELAY0以及其对应的GAP0存储到样本库中;
持续统计和测量网络环境约束值中次对时过程中的GAP值和DELAY值,选取所统计的网络环境约束值中的最小值为DELAY1和GAP1;
若DELAY1小于DELAY0,且GAP1在GAP0±(DELAY1-DELAY0)/收敛系数范围内,则所述DELAY1和GAP1更新为新的GPA0和DELAY0,标记此值可用,样本为有效样本;
若有效样本个数大于拟合库样本最大数,则对有效样本中的GAP值和 DELAY值进行卡尔曼滤波,以得到最优解的DELAYD值和GAPD值;
将所得到的最优解的DELAYD值和GAPD值作为最优网络延迟。
本发明实施例的第三方面提供了一种终端设备,包括存储器,处理器及存储在存储器上并可在处理器上运行的计算机程序,上述处理器执行上述计算机程序时实现上述第一方面提及的方法。
本发明实施例的第四方面提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,上述计算机程序被处理器执行时实现上述第一方面提及的方法。
本发明实施例与现有技术相比存在的有益效果是:本申请采用中断模式使数据包在客户端的物理层和应用层之间以及服务器的物理层和应用层之间传递,实现了数据包时间戳精确加盖,并且由于此算法在客户端和服务器都采用对称算法模式,其产生的执行误差可由于对称性大大减少。针对客户端和服务器之间数据包传递产生的误差通过计算最优网络延时对其进行调整,大大减小授时误差,提高网络授时精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的地震采集器的网络授时精度控制方法的流程示意图;
图2为本发明实施例提供数据包在客户端和服务器传递过程的示意图;
图3为本发明实施例提供的地震采集器的网络授时精度控制系统的结构示意图;
图4为本发明实施例提供的终端设备的结构示意图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。
应当理解,当在本说明书和所附权利要求书中使用时,术语“包括”指示所描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
还应当理解,在此本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样,除非上下文清楚地指明其它情况,否则单数形式的“一”、“一个”及“该”意在包括复数形式。
还应当进一步理解,在本发明说明书和所附权利要求书中使用的术语“和/ 或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。
为了说明本发明所述的技术方案,下面通过具体实施例来进行说明。
实施例一
图1是本发明实施例提供的地震采集器的网络授时精度控制方法的流程示意图,该方法可以包括以下步骤:
S101:对客户端接收到的数据包,采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层。
S102:所述客户端向服务器发型所述数据包时,预先设置网络环境约束值。
S103:连续检测不同时刻所述网络环境约束值次对时过程中的GAP值和 DELAY值,其中,GAP值为所述客户端和所述服务器之间的时钟偏移量, DELAY值为数据包在所述客户端和所述服务器之间往返过程中的时间延迟。
S104:对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟。
S105:根据所述最优网络延迟对网络授时过程中的传递误差进行校正后,由所述服务器的物理层接收所述数据包。
S106:通过所述中断模式将所述数据由服务器的物理层发送到服务器的应用层为数据包进行授时。
实施例二
下面结合具体例子对本申请提供的授时精度控制方法进行说明;图2为本发明另一实施例提供数据包在客户端和服务器传递过程的示意图。
在NTPV0协议基础上,采用此协议的客户端服务器模式,通信方式为UDP, 在此协议基础上,其授时源为服务器,客户端为地震采集器(一下简称为客户端),其为被授时单元,在服务器段,默认此时服务器已有精确时间授时,其时间来源可来自于GPS或者其他授时单元,在此不做讨论,但对于服务器的授时过程中的算法限制,其有一部分需要有应用层软件限制,在此特别提出。
另外,涉及到NTP协议中的其他常规名称和通用的同步过程,在本发明中不再赘述。
为达到上述发明特点,发明的技术解决方案如下:
如图2所示,其中在图示中,各个字符标号的含义为:T1:同步阶段客户端应用层打包后发出的时间;T2:同步阶段客户端物理层打包发出的时间;T3:同步阶段服务器物理层收到数据包的时间;T4:同步阶段服务器应用层收到数据包的时间;T5:同步阶段服务器应用层打包发出的时间;T6:同步阶段服务器物理层打包发出的时间;T7:同步阶段客户端物理层收到数据包的时间;T8:同步阶段客户端应用层收到数据包的时间。
进一步的,对于上述的字符标号,便于下面算法阐述,其衍生的字符标号包括:D1:同步阶段客户端从应用层打包到物理层打包所经历的时间间隔;D2: 同步阶段数据包从客户端到服务器之间的网络链路延迟;D3:同步阶段服务器从物理层接收到应用层接收的时间间隔;D4:同步阶段服务器应用层从接收到数据包到发出数据包的时间间隔;D5:V同步阶段服务器从应用层打包到物理层经历的时间间隔;D6:同步阶段数据包从服务器到客户端之间的网络链路延迟;D7: 同步阶段客户端从物理层接收到应用层接收的时间间隔。
进一步的,由于NTP的工作机制,其建立在UDP的基础上,是建立在网络应用层之上的,其授时误差主要来源于两个方面,一个是对于不同的操作系统,其数据包从物理层到应用层的数据不可控,且无法估计;另外一个是在同步过程中,其从客户端到服务器和从服务器到客户端的网络延时不可控,且随着不同的网络环境而改变,因此本算法是针对此上述NTP产生误差的两大途径进行的优化,从而实现授时精度的提高。
进一步的,对于数据包从物理层到网络层以及数据包从网络层到物理层的授时误差,提出一种基于硬件中断的应用层与物理层之间的延时控制算法,具体到客户端,其表述为:1,将物理层数据包发送设置为中断模式,取消轮训模式,并设置发送优先级为最高,从而保证发送时刻此任务不受其他任务打断; 2,T1时刻将NTP数据包打包,并填充必要字段,其发送时间字段填充为打包时刻获取的客户端自身时间CT1加1ms,启动内部定时器CTIMER1定时中断,中断的优先级设置为物理层优先级减1(即为仅次于物理层发送优先级),并且设置TIMER中断的时间为1ms;3,等待CTIMER1中断,当中断发生时,也即是T2时刻启动物理层数据包发生中断,数据包发出;4,设置物理层接收为中断模式,且中断优先级为最高,由于物理层发送和接收不可能同时发生,因此其两者不会发生冲突;5,当数据包到来时刻,自动触发物理层接收中断,此时将此刻获取的客户端自身时间存储到内存变量中,同时设置客户端数据包处理标志位,等待处理;6,应用层轮训执行序列设置NTP数据接收处理标志位为最前端,可保证其轮询队列可在第一时间处理其数据包。
进一步的,由于将客户端的数据包发送和接收时间戳加盖都是基于中断,且设置为优先级为最高,因此在主频设置为10M以上时,其执行后的时间精度可控制在1微妙以内。
进一步的,对于服务器端,其具体执行算法表述为:1,设置物理层接收为中断模式,其中断优先级为最高,其后服务器接收到物理层数据包后可自动触发;2,物理层接收到数据包后触发服务器中断,此时加盖服务器此时的时间戳,将接收时刻的时间存放在内存变量SVAR1中;3,延时,延时间隔设置为2ms; 4,准备发送数据包,此时将变量SVAR1填充到要发送的数据包中,并获取当前时间,将获取的当前时间填充到到SAVAR2中,并将SAVR2中的时间加1ms 后重新赋值给SAVR2;5,启动内部定时器STIMER1,将此定时器的中断优先级设置为最高,并将时间间隔设置为1ms,并设置物理层发送中断为最高优先级; 6,等待STIMER1中断,当中断发送时,启动物理层数据包发送中断,数据包发出。
进一步的,由于采用上述客户端和服务器发送接收算法,此时可实现,数据包时间戳精确加盖,并且由于此算法在客户端和服务器都采用对称算法模式,其产生的执行误差可由于对称性大大减少。
进一步的,针对网络延时产生的误差,采用如下的网络延时最优解估计算法,大大减少网络延时引起的误差。先将此网络延时最优解估计算法表述如下,在表述过程中涉及到的变量包括:1,GAP:客户端和服务器之间的时钟偏移;2,DELAY:从客户端到服务器以及服务器到客户端之间的往返的时间延迟。在此次两个变量的计算方式可从NTP数据包中提取相关字段就行计算,且计算方式为常规操作,在此不再赘述;3,NETCONSUM:为网络环境约束值,此值可由用户自行设定,当网络环境较差时可增加此值,来实现质量约束;4,STRVAR:收敛系数,此系数控制授时的收敛范围,默认设置为2,此系数选取越高,约束收敛效应越高,可控制排除异常值;5,LIBSUM:拟合库样本最大数,此值控制收敛库中最大样本数,大于此数的旧的样本将会被舍弃,通常设置此值为64; 6,FILSUM,滤波器样本数,采用卡尔曼滤波器的样本个数,采用卡尔曼滤波器对其进行平滑滤波的的样点数,通常选取32,可根据需要选择。
现将此网络延时最优解估计算法表述如下:1,设置NETCONSUM为固定常熟,在一般网络环境下,设置此值为16,如遇到网络环境变差情况,可加此值变高;2,统计和测量NETCONSUM次对时过程中的GAP和DELAY值;3,首次循环选取NETCONSUM中的DELAY最小值作为DELAY0,并将DELAY0 以及其对应的GAP0放入到样本库中,标记此值可用,并可用此值用作后续的时间校正计算;4,持续统计和测量NETCONSUM次对时过程中的GAP和 DELAY值,选取NETCONSUM中的最小值为DELAY1和GAP1;5,对DELAY1 和GAP1进行约束收敛判定,如若DELAY1小于DELAY0且GAP1在GAP0 ±(DELAY1-DELAY0)/STRVAR范围内,那么将此DELAY1和GAP1更新为新的GPA0和DELAY0,标记此值可用,并可用此值用作后续的时间校正计算,并且当有效样本大于LIBSUM时采用队列方式将之前最老的样本值替换为最新样本,如若GAP1不在此范围内,则将次值舍弃;6,统计卡尔曼滤波有效样本个数FILCNT,如若有效样本个数大于FILSUM时对其包含样本中的GAP和 DELAY值进行卡尔曼滤波,后得到最优解的DELAYD和GAPD,可应用于最后的时间校正,此时将FILCNT清零,后返回状态4继续执行,如若FILCNT小于FILSUM,则将FILCNT递增后返回状态4执行。
本发明的提供一种基于实时数据传输地震采集器的高精度网络授时算法,该算法可应用于具有网络传输硬件的数据采集器中,可在不改变硬件条件下通过该软件算法实现时间授时,并大幅的提高授时精度,且为达到最大兼容性,其算法可在NTPV0版本上运行,因此在此算法的实现基础不涉及后续高版本的NTP协议,由于NTP的高级版本是向下兼容的,因此用户在使用高级版本的 NTP协议时并不影响本算法的应用。
实施例三
图3是本发明另一实施例提供的地震采集器的网络授时精度控制方法,详述如下:
对于此基于实时数据传输地震采集器的高精度网络授时算法所依赖的客户端和服务器在此实例中均为基于ARM架构的Cortex-M3嵌入式MCU,其选用型号为意法半导体的STM32F2系列,在此此算法并不需要对MCU型号做出具体要求,仅需要此MCU有相应的网络物理层中断寄存器,可配置和响应此寄存器中断。对于客户端和服务器所采用的操作系统,在次算法中也不做强制依赖要求,但需要此操作系统中具备如下两个条件:1,有网络物理层中断挂接函数或者接口,可以在中断发生时执行相应的操作;2对于中断相应序列能做具体的优先级配置,在同时发生不同中断时,可根据不同的优先级来响应不同的中断函数,在此实例中没有采用现成的既有操作系统,仅采用基于中断驱动和状态机的嵌入式软件系统设计,其硬件中断可通过其优先级打断主状态机执行,且主状态机可通过设置不同的状态标志位实现执行序列的执行先后顺序。
进一步的,对于数据包从物理层到网络层以及数据包从网络层到物理层的授时误差,提出的一种基于硬件中断的应用层与物理层之间的延时控制算法,具体到客户端,此客户端MCU芯片采用STM32F207芯片,其芯片具备10/100 Ethernet MAC为物理层,且提供硬件中断向量,其具体表述为:1,将物理层数据包发送设置为中断模式,取消轮训模式,并设置发送优先级为最高,从而保证发送时刻此任务不受其他任务打断,在此实例中设置物理层发送中断主优先级为0,从优先级为0,也即在此芯片中为最高优先级;2,T1时刻将NTP数据包打包,并填充必要字段,在此NTP数据包打包包括NTP服务器地址填充, UDP必要数据包填充,其发送时间字段填充为打包时刻获取的客户端自身时间CT1加1ms,启动内部定时器CTIMER1定时中断,中断的优先级设置为物理层优先级减1(即为仅次于物理层发送优先级),并且设置CTIMER中断的时间为1ms,在此由于其发送时间字段在填充时已经加1ms,并且采用CTIMER中断可实现1ms延时精确控制,因此在CTIMER1发送时刻可以精确抵消1ms时间延迟;3,等待CTIMER1中断,当中断发生时,也即是T2时刻启动物理层数据包发生中断,数据包发出,在此将物理层数据发送放在CTIMER1中断中执行,由于其为同等优先级,其发送过程不会被打断;4,设置物理层接收为中断模式,且中断优先级为最高,在此设置接收主优先级为0,从优先级为0,由于物理层发送和接收不可能同时发生,因此其两者不会发生冲突;5,当数据包到来时刻,自动触发物理层接收中断,此时将此刻获取的客户端自身时间存储到内存变量中,在此实例中变量主要包括数据包中的四个时间变量,同时设置客户端数据包处理标志位,等待处理;6,应用层轮训执行序列设置NTP数据接收处理标志位为最前端,可保证其轮询队列可在第一时间处理其数据包,在此实例中主要通过轮训标志位设置其轮训处理优先级别,将其轮询处理级别设置为最高,可保证其在后续轮训过程中得到第一时间处理。
进一步的,由于将客户端的数据包发送和接收时间戳加盖都是基于中断,且设置为优先级为最高,因此在主频设置为10M以上时,其执行后的时间精度可控制在1微妙以内,在此实例中,其客户端晶振采用12M晶振,通过倍频器后其主时钟频率为48M。
进一步的,对于服务器端,在此实例中同样采用STM32F207作为主芯片,外围网络硬件跟客户端相同,同样具备10/100Ethernet MAC为物理层,且提供硬件中断向量,由于其需要更多的并发操作,设置倍频器倍频后的主频为120M,其具体执行算法表述为:1,设置物理层接收为中断模式,其中断优先级为最高,在此实例中设置其中断主优先级为0,从优先级为0,也即最高优先级,其后服务器接收到物理层数据包后可自动触发;2,物理层接收到数据包后触发服务器中断,此时加盖服务器此时的时间戳,将接收时刻的时间存放在内存变量SVAR1中,在此次时间戳为服务器接收时间,便于后续填充到给客户端返回的数据包字段中;3,延时,延时间隔设置为2ms,在此实例中设置为2ms,由于此延时间隔并非精确要求,且其精确度不会影响到整个的对钟流程的精确性,因此采用普通的延时算法,在此不再赘述;4,准备发送数据包,此时将变量 SVAR1填充到要发送的数据包中,并获取当前时间,将获取的当前时间填充到到SAVAR2中,并将SAVR2中的时间加1ms后重新赋值给SAVR2;5,启动内部定时器STIMER1,将此定时器的中断优先级设置为最高,并将时间间隔设置为1ms,并设置物理层发送中断为最高优先级,在此设置STIMER1和物理层发送中断的中断优先级为主优先级为0,从优先级为0;6,等待STIMER1 中断,当中断发送时,启动物理层数据包发送中断,数据包发出。
进一步的,由于采用上述客户端和服务器发送接收算法,此时可实现,数据包时间戳精确加盖,并且由于此算法在客户端和服务器都采用对称算法模式,其产生的执行误差可由于对称性大大减少。
进一步的,针对网络延时产生的误差,采用如下的网络延时最优解估计算法,大大减少网络延时引起的误差。先将此网络延时最优解估计算法表述如下,在表述过程中涉及到的变量在上节文档中均已介绍,因此在此不再复述。
现将此网络延时最优解估计算法表述如下:1,设置NETCONSUM为固定常数,在一般网络环境下,设置此值为16,如遇到网络环境变差情况,可加此值变高,在此实例中设置此值为16;2,统计和测量NETCONSUM次对时过程中的GAP和DELAY值,在此实例中将此变量作为大小为NETCONSUM的二维数组作为统计变量存入;3,首次循环选取NETCONSUM中的DELAY最小值作为DELAY0,并将DELAY0以及其对应的GAP0放入到样本库中,标记此值可用,并可用此值用作后续的时间校正计算,在此实例中比较上述数组中的最小值存入到DEALY0;4,持续统计和测量NETCONSUM次对时过程中的 GAP和DELAY值,选取NETCONSUM中的最小值为DELAY1和GAP1,在此实例中持续统计和测量的时间间隔为10秒,此间隔为用户自行设定,可根据其对钟频率设定,设定值越小,其对钟频次越快;5,对DELAY1和GAP1进行约束收敛判定,如若DELAY1小于DELAY0且GAP1在GAP0± (DELAY1-DELAY0)/STRVAR范围内,那么将此DELAY1和GAP1更新为新的 GPA0和DELAY0,标记此值可用,并可用此值用作后续的时间校正计算,在此实例中,由于DELAY1和DELAY0相互之间大小关系不确定,因此± (DELAY1-DELAY0)后产生的是范围值,需要用户注意,并且当有效样本大于LIBSUM时采用队列方式将之前最老的样本值替换为最新样本,在此实例中,采用的是先进后出的循环队列结构,因此其可由最新值替换掉最老的变量值,如若GAP1不在此范围内,则将次值舍弃,不影响此次计算;6,统计卡尔曼滤波有效样本个数FILCNT,如若有效样本个数大于FILSUM时对其包含样本中的GAP和DELAY值进行卡尔曼滤波,后得到最优解的DELAYD和GAPD,可应用于最后的时间校正,在此实例中,此DELAYD和GAPD即为可信可用的变量值,用户可应用此变量来进行时钟校正,此时将FILCNT清零,后返回状态4继续执行,如若FILCNT小于FILSUM,则将FILCNT递增后返回状态 4执行,此时还不满足新此卡尔曼滤波系数缓冲,因此不产生新的可信可用的 DELAYD和GAPD。
实施例四
图3为本申请提供的地震采集器的网络授时精度控制系统的结构示意图,为了便于说明,仅示出了与本发明实施例相关的部分。
所述地震采集器的网络授时精度控制系统,包括:
数据处理模块31,用于对客户端接收到的数据包,采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层;
最优网络延迟计算模块32,用于在所述客户端向服务器发型所述数据包时,预先设置网络环境约束值;连续检测不同时刻所述网络环境约束值次对时过程中的GAP值和DELAY值,其中,GAP值为所述客户端和所述服务器之间的时钟偏移量,DELAY值为数据包在所述客户端和所述服务器之间往返过程中的时间延迟;对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟;
校正授时模块33,用于根据所述最优网络延迟对网络授时过程中的传递误差进行校正后,由所述服务器的物理层接收所述数据包;并通过所述中断模式将所述数据由服务器的物理层发送到服务器的应用层为数据包进行授时。
实施例五
图4是本发明实施例四提供的终端设备的结构示意图。如图4所示,该实施例的终端设备4包括:处理器40、存储器41以及存储在所述存储器41中并可在所述处理器40上运行的计算机程序42。所述处理器40执行所述计算机程序42时实现上述方法实施例一中的步骤,例如图1所示的步骤S101至S106。所述处理器40执行所述计算机程序42时实现上述各装置实施例中各模块/单元的功能,例如图3所示模块31至33的功能。
示例性的,所述计算机程序42可以被分割成一个或多个模块/单元,所述一个或者多个模块/单元被存储在所述存储器41中,并由所述处理器40执行,以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述所述计算机程序42在所述终端设备4中的执行过程。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种地震采集器的网络授时精度控制方法,其特征在于,所述网络授时精度控制方法包括:
对客户端接收到的数据包,采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层;
所述客户端向服务器发型所述数据包时,预先设置网络环境约束值;
连续检测不同时刻所述网络环境约束值次对时过程中的GAP值和DELAY值,其中,GAP值为所述客户端和所述服务器之间的时钟偏移量,DELAY值为数据包在所述客户端和所述服务器之间往返过程中的时间延迟;
对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟;
根据所述最优网络延迟对网络授时过程中的传递误差进行校正后,由所述服务器的物理层接收所述数据包;
通过所述中断模式将所述数据由服务器的物理层发送到服务器的应用层为数据包进行授时。
2.根据权利要求1所述的网络授时精度控制方法,其特征在于,所述采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层,包括:
预先设置客户端的应用层的发送数据包的模式及物理层接收数据包的模式为中断模式,设置所述客户端的物理层优先级为最高,应用层中断的优先级为所述客户端的物理层优先级减1;
T1时刻客户端将NTP数据包打包,发送时间字段填充为打包时刻所获取的所述客户端自身时间CT1加1ms,等待中断发生时由应用层将数据包发出;
物理层检测到所述数据包到来时,接收中断,并将此刻获取的所述客户端的自身时间进行存储,对数据包进行标志位后等待处理。
3.根据权利要求1所述的网络授时精度控制方法,其特征在于,所述对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟,包括:
首次循环选取网络环境约束值中的DELAY值最小值作为DELAY0,并将DELAY0以及其对应的GAP0存储到样本库中;
持续统计和测量网络环境约束值中次对时过程中的GAP值和DELAY值,选取所统计的网络环境约束值中的最小值为DELAY1和GAP1;
若DELAY1小于DELAY0,且GAP1在GAP0±(DELAY1-DELAY0)/收敛系数范围内,则所述DELAY1和GAP1更新为新的GPA0和DELAY0,标记此值可用,样本为有效样本;
若有效样本个数大于拟合库样本最大数,则对有效样本中的GAP值和DELAY值进行卡尔曼滤波,以得到最优解的DELAYD值和GAPD值;
将所得到的最优解的DELAYD值和GAPD值作为最优网络延迟。
4.根据权利要求3所述的网络授时精度控制方法,其特征在于,所述网络授时精度控制方法还包括:
若DELAY1不小于DELAY0,或GAP1不在GAP0±(DELAY1-DELAY0)/收敛系数范围内,则舍弃此次统计的DELAY值和GAP值。
5.根据权利要求1至4任一项所述的网络授时精度控制方法,其特征在于,所述拟合库样本最大数为64。
6.一种地震采集器的网络授时精度控制系统,其特征在于,所述网络授时精度控制系统包括:
数据处理模块,用于对客户端接收到的数据包,采用中断模式使所述数据包由客户端的应用层接收所述数据包,并采用中断模式将所述数据包由客户端的应用层发送到客户端的物理层;
最优网络延迟计算模块,用于在所述客户端向服务器发型所述数据包时,预先设置网络环境约束值;连续检测不同时刻所述网络环境约束值次对时过程中的GAP值和DELAY值,其中,GAP值为所述客户端和所述服务器之间的时钟偏移量,DELAY值为数据包在所述客户端和所述服务器之间往返过程中的时间延迟;对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟;
校正授时模块,用于根据所述最优网络延迟对网络授时过程中的传递误差进行校正后,由所述服务器的物理层接收所述数据包;并通过所述中断模式将所述数据由服务器的物理层发送到服务器的应用层为数据包进行授时。
7.根据权利要求6所述的网络授时精度控制系统,其特征在于,所述数据处理模块具体用于,
预先设置客户端的应用层的发送数据包的模式及物理层接收数据包的模式为中断模式,设置所述客户端的物理层优先级为最高,应用层中断的优先级为所述客户端的物理层优先级减1;
T1时刻客户端将NTP数据包打包,发送时间字段填充为打包时刻所获取的所述客户端自身时间CT1加1ms,等待中断发生时由应用层将数据包发出;
物理层检测到所述数据包到来时,接收中断,并将此刻获取的所述客户端的自身时间进行存储,对数据包进行标志位后等待处理。
8.根据权利要求6所述的网络授时精度控制系统,其特征在于,所述最优网络延迟计算模块在对所述GAP值和DELAY值进行约束收敛判定,计算最优网络延迟时,具体用于:
首次循环选取网络环境约束值中的DELAY值最小值作为DELAY0,并将DELAY0以及其对应的GAP0存储到样本库中;
持续统计和测量网络环境约束值中次对时过程中的GAP值和DELAY值,选取所统计的网络环境约束值中的最小值为DELAY1和GAP1;
若DELAY1小于DELAY0,且GAP1在GAP0±(DELAY1-DELAY0)/收敛系数范围内,则所述DELAY1和GAP1更新为新的GPA0和DELAY0,标记此值可用,样本为有效样本;
若有效样本个数大于拟合库样本最大数,则对有效样本中的GAP值和DELAY值进行卡尔曼滤波,以得到最优解的DELAYD值和GAPD值;
将所得到的最优解的DELAYD值和GAPD值作为最优网络延迟。
9.一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至5中任一项所述方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至5中任一项所述方法的步骤。
CN202010244198.3A 2020-03-31 2020-03-31 地震采集器的网络授时精度控制方法、系统及终端设备 Expired - Fee Related CN111447030B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202010244198.3A CN111447030B (zh) 2020-03-31 2020-03-31 地震采集器的网络授时精度控制方法、系统及终端设备
PCT/CN2021/077592 WO2021196923A1 (zh) 2020-03-31 2021-02-24 地震采集器的网络授时精度控制方法、系统及终端设备
US17/475,369 US11716162B2 (en) 2020-03-31 2021-09-15 Method and system for controlling network timing precision of a seismic collector, and terminal device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010244198.3A CN111447030B (zh) 2020-03-31 2020-03-31 地震采集器的网络授时精度控制方法、系统及终端设备

Publications (2)

Publication Number Publication Date
CN111447030A true CN111447030A (zh) 2020-07-24
CN111447030B CN111447030B (zh) 2021-02-19

Family

ID=71650975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010244198.3A Expired - Fee Related CN111447030B (zh) 2020-03-31 2020-03-31 地震采集器的网络授时精度控制方法、系统及终端设备

Country Status (3)

Country Link
US (1) US11716162B2 (zh)
CN (1) CN111447030B (zh)
WO (1) WO2021196923A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113037418A (zh) * 2021-03-17 2021-06-25 广州虎牙科技有限公司 一种网络授时的误差校正方法及相关装置
WO2021196923A1 (zh) * 2020-03-31 2021-10-07 中国科学院地质与地球物理研究所 地震采集器的网络授时精度控制方法、系统及终端设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7363009B2 (en) * 2003-02-03 2008-04-22 Qualcomm Incorporated Method and apparatus for determining propagation delays for use in wide area networks
US7535931B1 (en) * 2003-07-14 2009-05-19 Symmetricom, Inc. Two-way timing protocol
US20090284373A1 (en) * 2007-07-31 2009-11-19 Nec Corporation Disaster noticing system, disaster noticing server, disaster reporting terminal method, and program
CN102647269A (zh) * 2011-02-22 2012-08-22 中国石油集团长城钻探工程有限公司 井下仪器的时间同步方法
US20140064303A1 (en) * 2012-09-04 2014-03-06 Khalifa University of Science, Technology, and Research Methods and devices for clock synchronization
CN108919634A (zh) * 2018-08-13 2018-11-30 中国科学院国家授时中心 一种北斗三频非差非组合观测值时间传递系统及方法
CN109167667A (zh) * 2018-09-18 2019-01-08 武汉海晟科讯科技有限公司 一种基于ptp同步的传感器数据采集传输系统和方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2799004B1 (fr) * 1999-09-27 2001-10-19 Inst Francais Du Petrole Methode et systeme de transmission utilisant un reseau de transmission standard, pour relier entre eux des elements d'un dispositif sismique
US6985499B2 (en) * 2000-04-20 2006-01-10 Symmetricom, Inc. Precise network time transfer
KR101323128B1 (ko) * 2012-03-13 2013-10-30 주식회사 다산네트웍스 마스터의 발신 메시지 개수 자동 조절을 통한 시간 동기화 장치 및 방법
CN104113517A (zh) * 2013-04-22 2014-10-22 华为技术有限公司 时间戳生成方法、装置及系统
US20160094298A1 (en) * 2014-09-25 2016-03-31 Seabed Geosolutions B.V. Wireless data transfer for an autonomous seismic node
KR20170069706A (ko) * 2015-12-11 2017-06-21 (주) 소암컨설턴트 미소지진 계측 시스템 및 이의 시간 동기화 방법
US10005523B2 (en) * 2016-01-19 2018-06-26 Ion Geophysical Corporation Unmanned marine vessel for node deployment and retrieval
MX2018008831A (es) * 2016-01-19 2018-11-09 Ion Geophysical Corp Embarcacion maritima no tripulada para la implementacion y recuperacion de nodos.
US10341083B2 (en) * 2016-09-09 2019-07-02 Huawei Technologies Co., Ltd. System and methods for network synchronization
US10104148B2 (en) * 2017-01-03 2018-10-16 Globalfoundries Inc. Nanosecond accuracy under precision time protocol for ethernet by using high accuracy timestamp assist device
US11050501B2 (en) * 2018-06-14 2021-06-29 Microchip Technology Incorporated Performing PHY-level hardware timestamping and time synchronization in cost-sensitive environments
US11974238B2 (en) * 2018-06-21 2024-04-30 Nokia Technologies Oy Time-synchronized radio bearer for supporting precision timing protocol (PTP) based time sensitive network (TSN) applications
CN108897210A (zh) * 2018-07-11 2018-11-27 成都爱为贝思科技有限公司 一种地震勘探数据无线采集和gps时间同步方法
CN111447030B (zh) * 2020-03-31 2021-02-19 中国科学院地质与地球物理研究所 地震采集器的网络授时精度控制方法、系统及终端设备

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7363009B2 (en) * 2003-02-03 2008-04-22 Qualcomm Incorporated Method and apparatus for determining propagation delays for use in wide area networks
US7535931B1 (en) * 2003-07-14 2009-05-19 Symmetricom, Inc. Two-way timing protocol
US20090284373A1 (en) * 2007-07-31 2009-11-19 Nec Corporation Disaster noticing system, disaster noticing server, disaster reporting terminal method, and program
CN102647269A (zh) * 2011-02-22 2012-08-22 中国石油集团长城钻探工程有限公司 井下仪器的时间同步方法
US20140064303A1 (en) * 2012-09-04 2014-03-06 Khalifa University of Science, Technology, and Research Methods and devices for clock synchronization
CN108919634A (zh) * 2018-08-13 2018-11-30 中国科学院国家授时中心 一种北斗三频非差非组合观测值时间传递系统及方法
CN109167667A (zh) * 2018-09-18 2019-01-08 武汉海晟科讯科技有限公司 一种基于ptp同步的传感器数据采集传输系统和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
许卫卫: "地震观测数据质量控制研究综述", 《地球物理学进展》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021196923A1 (zh) * 2020-03-31 2021-10-07 中国科学院地质与地球物理研究所 地震采集器的网络授时精度控制方法、系统及终端设备
CN113037418A (zh) * 2021-03-17 2021-06-25 广州虎牙科技有限公司 一种网络授时的误差校正方法及相关装置
CN113037418B (zh) * 2021-03-17 2023-06-20 广州虎牙科技有限公司 一种网络授时的误差校正方法及相关装置

Also Published As

Publication number Publication date
WO2021196923A1 (zh) 2021-10-07
US20220006548A1 (en) 2022-01-06
US11716162B2 (en) 2023-08-01
CN111447030B (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
CN111555834B (zh) 基于ieee1588协议的无线网络的时钟同步方法及系统
CN109818826B (zh) 一种网络路径延迟测量方法及其装置与时间同步系统
CN110249372A (zh) 传感系统及时间同步方法
CN111447030B (zh) 地震采集器的网络授时精度控制方法、系统及终端设备
CN1845546A (zh) 面向测量与控制的精确时间同步方法与系统
EP2448168A1 (en) Method and system for bearing time synchronization protocol in optical transport network
US20110222561A1 (en) Systems and methods for providing time synchronization
US11336383B2 (en) Packet scheduling system with desired physical transmission time for packets
CN105846941A (zh) 基于ieee 1588标准的时钟透明传输装置及方法
WO2021008532A1 (zh) 有源天线单元时延对齐方法、装置及有源天线单元
CN103248471A (zh) 基于ptp协议与反射内存网的时钟同步方法
US20220360350A1 (en) Method and apparatus for acquiring timestamp of data stream, storage medium, and electronic apparatus
CN112395231B (zh) 一种多处理器间的数据延迟时间计算方法及装置
CN113419598B (zh) 一种多soc系统时钟同步系统及方法
WO2012071977A1 (zh) 一种机架内的时间同步方法及设备
Akpinar et al. Improved clock synchronization algorithms for the controller area network (CAN)
CN108768573B (zh) 一种时钟同步方法及装置
CN106130711B (zh) 基于pac控制器的ieee1588对时方法及装置
CN102638339A (zh) 一种用于实现精确时间同步的方法和装置
CN103490856B (zh) 用于深空探测的Mark5B格式VLBI数据接收解码与纠错系统及方法
US11907754B2 (en) System to trigger time-dependent action
US20220269223A1 (en) Wireless time servicing method, device, and system
CN102983959B (zh) 在多个mac中实现ptp一步模式和两步模式的方法及装置
TW201729628A (zh) 頻率校正裝置及方法
CN214480655U (zh) 一种兼容可定义确定性通信以太网的嵌入式设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210219

CF01 Termination of patent right due to non-payment of annual fee