CN111064540A - 一种矿井微地震监测时钟同步装置及方法 - Google Patents
一种矿井微地震监测时钟同步装置及方法 Download PDFInfo
- Publication number
- CN111064540A CN111064540A CN201911406667.0A CN201911406667A CN111064540A CN 111064540 A CN111064540 A CN 111064540A CN 201911406667 A CN201911406667 A CN 201911406667A CN 111064540 A CN111064540 A CN 111064540A
- Authority
- CN
- China
- Prior art keywords
- clock
- control station
- signal acquisition
- field control
- acquisition unit
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0635—Clock or time synchronisation in a network
- H04J3/0638—Clock or time synchronisation among nodes; Internode synchronisation
- H04J3/0658—Clock or time synchronisation among packet nodes
- H04J3/0661—Clock or time synchronisation among packet nodes using timestamps
- H04J3/0667—Bidirectional timestamps, e.g. NTP or PTP for compensation of clock drift and for compensation of propagation delays
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0635—Clock or time synchronisation in a network
- H04J3/0638—Clock or time synchronisation among nodes; Internode synchronisation
- H04J3/0644—External master-clock
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Electric Clocks (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种矿山微地震监测***的时钟同步方法,基于GPS授时对地面中央记录***获得时钟一致性匹配,井下现场控制站则采用原子时钟获得时钟一致性匹配,基于现场总线发送同步时钟信号和时间戳协议,由于信号采集单元和现场控制站均采用独立的晶振作为时钟,时间(x,y)间的相互关系可用一系列成对的数据,采用最小二乘法一阶求出现场控制站、信号采集单元时钟矫正初始值,然后利用二分法根据给定的精度要求对信号采集单元时钟频率多次校准,本发明的有益效果:可以满足矿山大规模微地震长时间连续监测的时钟同步要求,较现有技术,克服持续长时间采集时钟误差累积大,具有实现成本低、同步精度高、可靠性好和适用面广的特点。
Description
技术领域
本发明涉及一种矿井微地震监测***的时钟同步方法。具体地,本发明基于GPS时钟或原子钟对现场控制站授时,并对采集单元时钟频率进行校准的同步机制及方法。
背景技术
矿山开采过程中,由于矿山井巷或者采场周围矿体和围岩应力被破坏而发生弹性变形,超过了岩体能承受高压的脆性岩体极限平衡状态,向自由空间突然释放巨大能量,因此而产生的以急剧、猛烈的破坏为特征的动力现象称之矿山动力灾害,如矿山冲击地压、冒顶、矿震、水害突出、煤与瓦斯突出等。微地震监测***就是一种及时的动力灾害监测手段,具有远距离、动态、三维、实时监测的特点,能够给出震动后的各种信息,并根据震源情况确定破裂尺度和性质,具有不损伤岩体、劳动强度小、时间和空间连续,还可以根据震源情况进一步分析破裂尺度、强度和性质等优点,能够为评价范围内的矿山动力灾害隐患提供依据。
矿井微地震监测***是由地面中央控制***及井下现场控制站构成的网络体系,微地震信号采集时间上的误差将成为影响地震数据成像的主要因素,如何保证设备之间的时钟同步是需要解决的关键技术之一。现有地震勘探采集记录***的时钟同步方法,采用单次采集前的同步校准机制,即在开始采集前,从上至下发送同步码字,收到同步码字后,所有采集站在相同的时刻开始采集,从而保证能够在同步校准之后的短时间内实现全网同步采集和传输。但这种机制不适合长时间矿井微地震持续采集,这是由于矿井微地震现场控制站和信号采集单元采用独立晶振作为时钟源,长时间连续采集会产生时钟漂移会造成时钟误差累积,无法满足矿山安全微地震监测时钟一致性需求。
发明内容
本发明目的在于克服现有技术的不足,提供一种矿井微地震监测时钟同步装置及方法,实现中央记录***、现场控制站与基于现场总线技术的信号采集单元全局时钟同步,实现矿井微地震持续长时间采集时钟同步要求。
为实现以上目的,本发明采用如下技术方案。
首先,矿井微地震监测***包括用于地面中央记录***、矿山通讯***、现场控制站、信号采集单元、电源工作站、现场总线和时钟匹配装置。其中现场控制站、信号采集单元、电源工作站、现场总线和时钟匹配装置构建成矿山井下微地震采集子***。
所述的时钟匹配装置为高精度的GPS授时或原子钟,主要为地面中央记录***、现场控制站提供同步时间脉冲和精准时间戳协议,对于地面的监测工作站、超级计算中心可利用GPS授时获得时钟一致性匹配,井下采掘工作面现场控制站采用原子时钟获得时钟一致性匹配,现场控制站通过现场总线实现不同信号采集单元的时钟一致性匹配,即发送同步时钟信号和时间戳协议。
所述基于现场总线技术的信号采集单元和现场控制站均采用独立的晶振作为时钟,时间(x,y)间的相互关系可用一系列成对的数据(x1,y1),(x2,y2)……(xm,ym),设一次系数方程:
y=a0+a1x (1)
其中:a0、a1是任意实数。为确定a0和α1,采用最小二乘法构建函数
φ=∑(yi-yj)2 (2)
式中yi为实测值,yj=a0+a1xi为理论值,把(1)代入(2)中得:
φ=∑(yi-a0-a1xi)2 (3)
为得到∑(yi-yj)2最小值,对函数φ对a0、a1求偏导数,令这两个偏导数等于零。
亦即:
∑yi=ma0+(∑xi)a1 (6)
∑xiyi=(∑xi)a0+(∑xi)a1 (7)
得到的两个关于a0、a1为未知数的两个方程组,解这两个方程组得:
把a0、a1代入试(1)中得到一元线性回归方程。由于回归的关联式不可能通过每个回归数据点(x1,y1),(x2,y2)……(xm,ym),借助相关系数“R”,统计量“m”,进行判断关联式优劣;“R”越趋近于1越好;“m”的绝对值越大越好,并得到剩余标准偏差值。
在式(10)中,m为样本容量;xi、yi分别为任意一组实验数据x、y的数值。借助最小二乘法一阶拟合求解出晶振频率变化情况下的现场控制站、信号采集单元矫正时钟。
考虑到矿井微地震采集记录***长时间连续数据采集时,温度、老化、激励电平、晶体单元感受到的加速度变化和运算的拟合精度等因素的影响,GPS或原子时钟、现场控制站、信号采集单元时间不是单纯线性的关系。最小二乘法一阶拟合后,授时同步装置、现场控制站、信号采集单元时间误差比较大,而且由于授时同步装置、现场控制站、信号采集单元时间不是单纯线性的关系,用一阶拟合参数会存在明显的累计误差。故在上述计算基础上引入二分法。
将方程y=f(x)有根区间[a0,b0]分为两个小区间(称二分),取a0、b0的中点x1=(a0+b0)/2,并计算f(x1)的值;如果f(x1)与f(a0)同号,则方程的根必在[x1,b0]区间,反之,f(x1)与f(a0)异号,则根在[a0,x1]区间。通过这样的方法找出并确定新的有根区间[a1,b1],然后再将新的有根区间二分为两个小区间,如此继续下去,就可得到一个有根区间
直到某区间[ak,bk]的长度小于给定的精度ε,即根的绝对误差小于ε,变可以得到一个满意的近似值。
|x*-xk|≤(bk-ak)≤ε (12)
上述的同步机制过程为概述为最小二乘法一阶拟合求出拟合参数,求出的拟合参数作为现场控制站、信号采集单元时钟矫正初始值,然后利用二分法根据给定的精度要求进而多次矫正。
本发明的有益效果:基于GPS或原子时钟高精准的同步脉冲和时间戳信息,最小二乘法和二分法对现场控制站、信号采集单元时钟频率进行校准的同步机制,可以满足矿山大规模微地震长时间连续实时监测***的时间同步要求。较现有技术,克服持续长时间采集时钟误差累积大,具有实现成本低、同步精度高、可靠性好和适用面广的特点,从而获得同步效果较好的微地震采集数据。
附图说明
图1为矿山区域微地震监测***示意图。
图2时钟同步控制原理框图。
图1中:1、地面中央记录***;2、矿井调度室;3、电源工作站;4、现场控制站;5、时钟匹配装置;6、信号采集单元;7、通讯网络;8、矿山通讯***;9、现场总线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面通过附图和实施例对本发明的技术方案做进一步的详细描述。
首先,根据拟监测矿井区域的地质条件和开采情况布置微地震观测***示意图,如图1所示,其可通过以下步骤实现。
第一步:根据矿井地质条件和开采情况设计并布置微地震监测***,包括地面中央记录***、矿井调度室和现场采集子***。
第二步:采用GPS授时为地面中央记录***、现场控制站提供同步时间脉冲和时间戳信息,井下现场采集子***采用原子钟获得时间脉冲一致性匹配和时间戳信息,通过现场总线自动实现不同信号采集单元的时钟一致性匹配,即发送同步时钟信号和时间戳协议。
第三步:地面中央记录***、矿井调度室和现场控制站接收时钟匹配装置的时间戳协议,并解析时间参数,并将得到的结果更新自身***时间,基于通讯网络或矿山通讯***向现场控制站发送自身修正时钟,保持时钟一致。
第四步:现场控制站基于现场总线向信号采集单元发送时间戳信息,信号采集单元接收时间戳信息,结合信号采集单元独立晶振时钟的相互关系形成一系列成对的数据(x1,y1),(x2,y2)……(xm,ym),其中yi为信号采集单元时钟实测值,xi为现场控制站时间戳信息,根据公式(1)构建相互关系。
第五步:信号采集单元持续接收一段时间现场控制站时间戳信息,如100个时间戳信息,采用最小二乘法构建函数,根据公式(9)解得相互关系未知数a1;进一步根据公式(8)解得相互关系未知数a0,从而得到信号采集单元时钟矫正初始值。
第六步:利用公式(10)评价为了判断时间戳信息和信号采集单元关联式的好坏,并得到剩余标准偏差值±σ。
第七步:由于微地震信号采集单元长时间连续数据采集时,温度、老化、激励电平、网络延时单元、晶体单元感受到的加速度变化和运算的拟合精度等因素的影响,由于时间GPS、FCU时间不是单纯线性的关系,在最小二乘法的基础上再引入二分法。由此,根据矫正初始值计算ym+1-(a0+a1xm+1),若大于0,则对标准偏差区间+σ之间,并进一步计算ym+1-(a0+a1xm+1+σ/2),如此继续,就可得到一个小于给定的精度ε的近似值。反之,小于0,则对标准偏差区间-σ之间,并进一步计算ym+1-(a0+a1xm+1-σ/2),如此同样可得到一个小于给定的精度ε的近似值。
此时信号采集单元得到最终的修正时间,GPS时钟地面中央记录***、原子钟授时、现场控制站和微地震信号采集单元时间精度误差小于300ns,满足矿井微地震数据采集时钟同步要求。
以上给出了具体的实施方式,但本发明不局限于所描述的实施方式。本发明的基本思路在于上述基本方案,对本领域普通技术人员而言,根据本发明设计出各种变形的模型、公式、参数并不需要花费创造性劳动。在不脱离本发明的原理和精神的情况下对实施方式进行的变化、修改、替换和变型都应涵盖在本发明的保护范围之内。
Claims (2)
1.一种矿井微地震监测时钟同步装置及方法,其特征在于:包括用于地面中央记录***、矿山通讯***、现场控制站、信号采集单元、电源工作站、现场总线和时钟匹配装置,信号采集单元和现场控制站均采用独立的晶振作为时钟,时间(x,y)间的相互关系可用一系列成对的数据(x1,y1),(x2,y2)……(xm,ym),采用最小二乘法构建函数得到一元线性回归方程拟合参数,求解出晶振频率变化情况下的现场控制站、信号采集单元时钟的矫正初始值,然后利用二分法多次迭代矫正获得给定的同步精度。
2.根据权利要求1所述一种矿井微地震监测时钟同步装置及方法,其特征在于:时钟匹配装置主要为地面中央记录***、现场控制站提供同步时间脉冲,对于地面的监测工作站、超级计算中心可利用GPS授时获得时钟一致性匹配,井下采掘工作面则采用原子时钟获得时钟一致性匹配,通过现场总线自动实现不同信号采集单元的时钟一致性匹配,即发送同步时钟信号和时间戳协议。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911406667.0A CN111064540A (zh) | 2019-12-31 | 2019-12-31 | 一种矿井微地震监测时钟同步装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911406667.0A CN111064540A (zh) | 2019-12-31 | 2019-12-31 | 一种矿井微地震监测时钟同步装置及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111064540A true CN111064540A (zh) | 2020-04-24 |
Family
ID=70305264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911406667.0A Pending CN111064540A (zh) | 2019-12-31 | 2019-12-31 | 一种矿井微地震监测时钟同步装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111064540A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112713955A (zh) * | 2020-12-22 | 2021-04-27 | 福州华虹智能科技股份有限公司 | 一种基于矿井的精准时间同步方法 |
CN113609971A (zh) * | 2021-08-04 | 2021-11-05 | 广州威拓电子科技有限公司 | 微地震观测设备的检查方法、装置、设备及存储介质 |
CN114739445A (zh) * | 2022-01-27 | 2022-07-12 | 厦门万宾科技有限公司 | 一种城市级排水管网增强扫描方法及*** |
CN115900635A (zh) * | 2023-03-09 | 2023-04-04 | 四川省交通勘察设计研究院有限公司 | 一种隧道变形数据监测方法、装置与*** |
WO2023061513A1 (zh) * | 2022-03-30 | 2023-04-20 | 河南理工大学 | 一种互联式矿井突水灾害微地震监测*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100322267A1 (en) * | 2007-11-30 | 2010-12-23 | Mitsubishi Electric Corporation | Time synchronization system and time synchronization method in train |
CN107395305A (zh) * | 2017-06-21 | 2017-11-24 | 浙江理工大学 | 面向实时Wi‑Fi网络的时间同步方法 |
CN107613458A (zh) * | 2017-09-07 | 2018-01-19 | 西安电子科技大学 | 一种tdoa条件下最优联合时间同步与定位的定位方法 |
CN108693553A (zh) * | 2017-04-06 | 2018-10-23 | 河南理工大学 | 基于互联网与超级计算机的矿山微震监测*** |
-
2019
- 2019-12-31 CN CN201911406667.0A patent/CN111064540A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100322267A1 (en) * | 2007-11-30 | 2010-12-23 | Mitsubishi Electric Corporation | Time synchronization system and time synchronization method in train |
CN108693553A (zh) * | 2017-04-06 | 2018-10-23 | 河南理工大学 | 基于互联网与超级计算机的矿山微震监测*** |
CN107395305A (zh) * | 2017-06-21 | 2017-11-24 | 浙江理工大学 | 面向实时Wi‑Fi网络的时间同步方法 |
CN107613458A (zh) * | 2017-09-07 | 2018-01-19 | 西安电子科技大学 | 一种tdoa条件下最优联合时间同步与定位的定位方法 |
Non-Patent Citations (1)
Title |
---|
周晓铭,陶俊才,秦海鸥: "基于二分法的曲线等误差圆弧拟合算法及其应用", 《南昌大学学报(理科版)》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112713955A (zh) * | 2020-12-22 | 2021-04-27 | 福州华虹智能科技股份有限公司 | 一种基于矿井的精准时间同步方法 |
CN112713955B (zh) * | 2020-12-22 | 2023-09-19 | 福州华虹智能科技股份有限公司 | 一种基于矿井的精准时间同步方法 |
CN113609971A (zh) * | 2021-08-04 | 2021-11-05 | 广州威拓电子科技有限公司 | 微地震观测设备的检查方法、装置、设备及存储介质 |
CN113609971B (zh) * | 2021-08-04 | 2024-07-05 | 广州威拓电子科技有限公司 | 微地震观测设备的检查方法、装置、设备及存储介质 |
CN114739445A (zh) * | 2022-01-27 | 2022-07-12 | 厦门万宾科技有限公司 | 一种城市级排水管网增强扫描方法及*** |
CN114739445B (zh) * | 2022-01-27 | 2023-12-15 | 厦门万宾科技有限公司 | 一种城市级排水管网增强扫描方法及*** |
WO2023061513A1 (zh) * | 2022-03-30 | 2023-04-20 | 河南理工大学 | 一种互联式矿井突水灾害微地震监测*** |
CN115900635A (zh) * | 2023-03-09 | 2023-04-04 | 四川省交通勘察设计研究院有限公司 | 一种隧道变形数据监测方法、装置与*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111064540A (zh) | 一种矿井微地震监测时钟同步装置及方法 | |
WAN et al. | An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence | |
CN101611330B (zh) | 进行井间观测的技术和*** | |
AU2009311440B2 (en) | Practical autonomous seismic recorder implementation and use | |
De Martino et al. | GPS time series at Campi Flegrei caldera (2000-2013) | |
Langbein | Deformation of the Long Valley Caldera, California: Inferences from measurements from 1988 to 2001 | |
CN108415066B (zh) | 一种隧道施工地质灾害预报方法 | |
Puglisi et al. | Ground deformation patterns on Mount Etna, 1992 to 1994, inferred from GPS data | |
Aloisi et al. | The 24 December 2018 eruptive intrusion at Etna volcano as revealed by multidisciplinary continuous deformation networks (CGPS, borehole strainmeters and tiltmeters) | |
CN103471544A (zh) | 一种基于北斗的高精度位移形变监测应用*** | |
CN110609321B (zh) | 一种基于速度模型数据库的微震源定位方法 | |
CN104992054B (zh) | 基于时间序列二维化的电离层垂直总电子含量预报方法 | |
Avallone et al. | Coseismic displacement waveforms for the 2016 August 24 Mw 6.0 Amatrice earthquake (central Italy) carried out from High-Rate GPS data | |
Odbert et al. | Chapter 11 volcano geodesy at the Soufrière Hills Volcano, Montserrat: A review | |
Kuzmin | Paradoxes of the comparative analysis of ground-based and satellite geodetic measurements in recent geodynamics | |
CN104166159A (zh) | 四维微地震监测的裂缝形态处理方法和*** | |
Butyrin et al. | Digital seismic logger “Ermak-5”. Experience of development and implementation | |
Crosson et al. | The southern Italy earthquake of 23 November 1980: an unusual pattern of faulting | |
Murray et al. | Coseismic displacements: 1992 landers, California, earthquake | |
Xu et al. | Optimal design of microseismic monitoring networking and error analysis of seismic source location for rock slope | |
WO2022001368A1 (zh) | 雷达及雷达*** | |
Crowell et al. | Validation of peak ground velocities recorded on very-high rate GNSS against NGA-WEST2 ground motion models | |
Bonaccorso et al. | Continuous tilt monitoring: Lesson learned from 20 years experience at Mt. Etna | |
Alparone et al. | High precision locations of multiplets on south-eastern flank of Mt. Etna (Italy): reconstruction of fault plane geometry | |
Abramov et al. | Key elements of construction and exploitation of a land-based seismic-gravimetric complex |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200424 |
|
RJ01 | Rejection of invention patent application after publication |