CN114325823A - 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 - Google Patents
一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 Download PDFInfo
- Publication number
- CN114325823A CN114325823A CN202111666298.6A CN202111666298A CN114325823A CN 114325823 A CN114325823 A CN 114325823A CN 202111666298 A CN202111666298 A CN 202111666298A CN 114325823 A CN114325823 A CN 114325823A
- Authority
- CN
- China
- Prior art keywords
- microseismic
- data
- iot
- time window
- effective
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000011435 rock Substances 0.000 title claims abstract description 68
- 238000012544 monitoring process Methods 0.000 title claims abstract description 33
- 239000013598 vector Substances 0.000 claims abstract description 52
- 230000005059 dormancy Effects 0.000 claims abstract description 36
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 238000007906 compression Methods 0.000 claims abstract description 8
- 230000006835 compression Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims description 38
- 230000001133 acceleration Effects 0.000 claims description 36
- 238000006243 chemical reaction Methods 0.000 claims description 24
- 238000004891 communication Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 22
- 238000000605 extraction Methods 0.000 claims description 13
- 230000003321 amplification Effects 0.000 claims description 9
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 9
- 238000013144 data compression Methods 0.000 claims description 8
- 238000012806 monitoring device Methods 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 5
- 239000007924 injection Substances 0.000 claims description 5
- 230000003213 activating effect Effects 0.000 claims description 2
- 238000004806 packaging method and process Methods 0.000 claims description 2
- 230000001902 propagating effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 14
- 230000002265 prevention Effects 0.000 abstract description 4
- 239000003990 capacitor Substances 0.000 description 15
- 101100236764 Caenorhabditis elegans mcu-1 gene Proteins 0.000 description 12
- 230000007958 sleep Effects 0.000 description 11
- 101100434411 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) ADH1 gene Proteins 0.000 description 10
- 101150102866 adc1 gene Proteins 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 4
- HCUOEKSZWPGJIM-YBRHCDHNSA-N (e,2e)-2-hydroxyimino-6-methoxy-4-methyl-5-nitrohex-3-enamide Chemical compound COCC([N+]([O-])=O)\C(C)=C\C(=N/O)\C(N)=O HCUOEKSZWPGJIM-YBRHCDHNSA-N 0.000 description 3
- 101001109689 Homo sapiens Nuclear receptor subfamily 4 group A member 3 Proteins 0.000 description 3
- 101000598778 Homo sapiens Protein OSCP1 Proteins 0.000 description 3
- 101001067395 Mus musculus Phospholipid scramblase 1 Proteins 0.000 description 3
- 102100022673 Nuclear receptor subfamily 4 group A member 3 Human genes 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000006378 damage Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000087 stabilizing effect Effects 0.000 description 3
- 239000007822 coupling agent Substances 0.000 description 2
- 238000005034 decoration Methods 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- HJVCHYDYCYBBQX-HLTLHRPFSA-N (2s,3s,4e,6e,8s,9s)-3-amino-9-methoxy-2,6,8-trimethyl-10-phenyldeca-4,6-dienoic acid Chemical compound OC(=O)[C@@H](C)[C@@H](N)/C=C/C(/C)=C/[C@H](C)[C@@H](OC)CC1=CC=CC=C1 HJVCHYDYCYBBQX-HLTLHRPFSA-N 0.000 description 1
- 102100024348 Beta-adducin Human genes 0.000 description 1
- 101000689619 Homo sapiens Beta-adducin Proteins 0.000 description 1
- 101710096655 Probable acetoacetate decarboxylase 1 Proteins 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- QVFWZNCVPCJQOP-UHFFFAOYSA-N chloralodol Chemical compound CC(O)(C)CC(C)OC(O)C(Cl)(Cl)Cl QVFWZNCVPCJQOP-UHFFFAOYSA-N 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000006837 decompression Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000368 destabilizing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- PBLZLIFKVPJDCO-UHFFFAOYSA-N omega-Aminododecanoic acid Natural products NCCCCCCCCCCCC(O)=O PBLZLIFKVPJDCO-UHFFFAOYSA-N 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000002618 waking effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Arrangements For Transmission Of Measured Signals (AREA)
Abstract
本发明提供一种基于NB‑IOT的岩体破裂失稳微震信号无线监测方法与装置，属于地质灾害防治技术领域，所述方法包括：采集岩体破裂产生的微震信号，使用长短时窗幅值比值法拾取有效微震波形数据，从有效微震波形数据中提取微震特征向量，根据微震特征向量执行不同的压缩策略得到压缩数据包，使用休眠唤醒算法对NB‑IOT模块进行调度，在需要发送数据的时候唤醒，而无需发送数据时进行休眠；在NB‑IOT模块处于唤醒状态时，将数据包发送至云平台，云平台对所接收的数据包进行分析，若满足预警条件则发出预警信息。本发明结合NB‑IOT技术解决岩体破裂失稳微震信号无线监测的数据远距离传输问题，同时也延长NB‑IOT模块工作时间，对于地质防灾减灾具有重要的实用价值。
Description
技术领域
本发明涉及地质灾害防治技术领域，尤其涉及一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法与装置。
背景技术
危岩失稳崩塌具有较高的突发性和较强的破坏性，强大的冲击力直接引起下方建筑物的垮塌、破坏，严重影响公路、铁路交通的正常运营，造成巨大的财产损失，而且近些年，频发旅游景区危岩崩塌和工程施工现场人员伤亡事件，对人民群众的生命造成了严重的威胁。微震(Microseism,MS)是频率小于100赫兹的低频声信号，是岩体在受到外界扰动应力及温度等因素的影响下，岩体内部出现应力集中，引起岩体微观裂隙的产生、扩展、贯通过程中伴随的低能量的弹性波或应力波。同时，在失稳崩塌演变的过程中，无论是从塌方演变初期——微裂纹发育阶段，还是后期——宏观破裂阶段，均包含着丰富的低频段微震弹性波，微震信号监测设备能够有效探测到其弹性波的大小及发生部位，并通过对危岩体破裂微震信号进行分析，揭示塌方发生的可能性及破坏剧烈程度，采取有效的防治或提前规避手段降低其所造成危害大小。到目前为止，微震监测技术已经成为了众多演示工程灾害监测和预报的重要手段，如煤矿的开采、水电站的建设、公路旁危岩体的监测等场景。传统的微震监测技术，主要采用有线的方式，传感器需要和采集仪等设备使用线缆连接起来，采集仪和服务器之间也主要通过有线光缆进行传输，这对一些较为偏远的工程实际来说，不仅需要架设电缆，还需要架设通讯光缆，使其使用代价较大。而传统的无线微震监测方式，不仅受限于无线通讯技术速率和延迟，无法做到实时传输和检测；同时也由于无线传输设备的功耗过大而导致工作时间过短不利于在室外环境使用。
NB-IOT(Narrow Band Internet of Things)是物联网领域一种新兴的技术，支持低功耗设备在广域网的蜂窝数据连接，也被称作低功耗广域网(LPWAN)。相比起其他不同的物联网无线传输技术如LoRa、ZigBee等，其优点在于可直接部署于GSM网络、UMTS网络或LTE网络，能实现平滑升级，因此其具有覆盖范围广阔的优点，在野外、山区等地方均能被覆盖到；NB-IOT具有低功耗的优点，特别适合于一些不能经常更换电池的场合，如安置于高山荒野偏远地区；此外，NB-IOT还具有低成本的优点，有利于在生产环境中大规模应用。
然而，现有的NB-IOT技术结合到岩体破裂微震信号无线监测这一具体的任务上，尚未形成一个有效的解决方案，尚存在诸多待解决的问题。首当其冲的是通信模块的功耗问题，由于其安装在偏远的野外地区，经常更换电池是不切实际的做法，只能选用低功耗的通信模块。然而对于微震信号传输这一任务来说，区别于传统的物联网应用的低频率数据传输场景，其数据传输的频率呈现了短时高频率的特点，这给通信模块的长时间工作带来了不小的挑战。除此之外，岩体破裂时产生的微震信号复杂多样，存在一些值得进一步分析的微震信号，如何将有分析价值的微震信号传输到云平台也是一个具有现实意义的问题。由于NB-IOT的窄带宽的特点，其在发送容量较大的数据包的时候存在延迟较高、阻塞信道甚至发送失败的现象。
发明内容
本发明的目的在于提供一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法与装置，解决现有岩体破裂失稳微震信号有线监测方案中因通讯线缆导致使用场合受限的技术问题。以及现有无线传输装置因功耗过高而导致无法安装在偏远地区的问题。实现了将岩体破裂失稳产生的微震信号通过NB-IOT技术远程地发送到云端，确保了信号发送的及时性同时延长了NB-IOT模块的工作时间。
为了实现上述目的，本发明采用的技术方案如下：
一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，所述方法包括如下步骤：
步骤1：在待监测岩体上安装一个或者一个以上的数据采集模块，采集岩体破裂产生的微震信号，使用长短时窗幅值比值法拾取有效微震波形数据；
步骤2：创建数据发送缓冲队列，从有效微震波形数据中提取微震特征向量，根据微震特征向量执行不同的压缩策略得到压缩数据包，将压缩数据包放置于数据发送缓冲队列尾部，同时发出一个携带了压缩数据包信息的信号量；
步骤3：使用休眠唤醒算法对NB-IOT模块进行调度，使其在需要发送数据的时候被唤醒，而不需要发送数据的时候进行休眠，节省电量以延长工作时间；
步骤4：当NB-IOT模块处于唤醒状态时，与云平台模块建立通信链路，将所述的压缩数据包发送至云平台模块；
步骤5：云平台模块对接收到的数据包进行分析，若满足预警条件则发出预警信息。
进一步地，步骤1中，采集岩体破裂产生的微震信号的具体过程为：
步骤1.1.1：当岩体产生微观裂隙时会释放出低能量的弹性波并向四周传播，引发岩体的微小震动，数据采集模块的加速度传感器能够将岩体微小震动的加速度值转换为电压值得到微震信号，转换公式如式子(1)所示：
式中，d表示压电常数，m表示加速度传感器内部压电元件的质量，C表示压电元件两端的电容，在加速度传感器型号确定的情况下其d、m和C均为确定常量，因此电压值V与加速度a成正比，所采集的微震信号的电压值越大，表示岩体震动的加速度越大；
步骤1.1.2：使用所述数据采集模块的A/D转换单元将加速度传感器输出的模拟信号转换为数字化的幅值并传输至数据处理模块作进一步处理。
进一步地，步骤1中，使用长短时窗幅值比值法拾取有效微震波形数据的具体过程为：
步骤1.2.1：创建微震数据包暂存队列、有效波形暂存队列、短时窗数组和长时窗数组，所述的短时窗数组和长时窗数组的最大长度分别为M和N，设有效信号提取阈值为R；所述的M、N和R根据监测现场环境决定；
步骤1.2.2：初始化长时窗数组，逐个读取从所述数据采集模块传输过来的幅值，将其按顺序保存于长时窗数组内，直至长时窗数组的元素个数达到N；
步骤1.2.3：逐个读取从所述数据采集模块传输过来的幅值，将其按顺序保存于短时窗数组内，当短时窗数组内的元素达到M时，分别计算短时窗数组和长时窗数组内的平均电压幅值Amplitude1和Amplitude2，接着将长时窗数组的元素整体往前搬移M个位置，将短时窗数组内的元素依次添加到长时窗数组末尾的M个位置上，然后求得Amplitude1和Amplitude2的比值r＝Amplitude1/Amplitude2，若r大于等于R则表示开始拾取有效微震波形数据，接着转步骤1.2.4：；否则，清空短时窗数组内的元素并重复执行步骤1.2.3；
步骤1.2.4：将短时窗数组内的元素依次添加到有效波形暂存队列尾部中；接着清空短时窗数组内的元素，然后继续逐个读取从所述数据采集模块传输过来的幅值，当短时窗数组内的元素达到M后，分别计算短时窗数组和长时窗数组的平均电压幅值Amplitude1和Amplitude2并求得二者的比值r＝Amplitude1/Amplitude2，若r小于R则表示结束拾取有效的微震波形数据，转步骤1.2.5，否则，将长时窗数组的元素整体往前搬移M个位置，将短时窗数组内的元素依次添加到长时窗数组末尾的M个位置上，重复执行步骤1.2.4；
步骤1.2.5：将短时窗数组内的元素按顺序添加到有效波形暂存队列尾部中，获取当前的时间戳，将有效波形暂存队列的所有元素和时间戳打包成一个微震数据包，并将其放置于微震数据包暂存队列尾部，表示当前拾取了一个有效微震波形数据，接着清空有效波形暂存队列，然后重复执行步骤1.2.3，开始拾取下一个有效微震波形数据。
进一步地，步骤2中，从有效微震波形数据中提取微震特征向量的具体过程为：
创建一个微震波形数据处理数组，从所述的微震数据包暂存队列的头部中取出一个微震数据包，将其中的有效微震波形数据逐元素复制到微震波形数据处理数组中；从微震波形数据处理数组的第一个元素开始遍历，到最后一个元素结束遍历，计算得到触发位置、结束位置、信号持续时间、上升时间、最大幅值和最大幅值的位置共六个特征值，将这六个特征值组成微震特征向量，六个特征值的具体定义为：
触发位置：表示数组中第一个大于阈值的幅值对应的数组下标；
结束位置：表示数组中最小幅值对应的数组下标；
信号持续时间：表示结束位置与触发位置的差值；
最大幅值：表示数组所有幅值中最大的幅值；
最大幅值的位置：表示数组中最大幅值所对应的数组下标；
上升时间：表示最大幅值的位置与触发位置的差值。
进一步地，步骤2中，根据微震特征向量执行不同的压缩策略得到压缩数据包的具体过程为：
若所述微震特征向量的最大幅值大于或等于设定的阈值Vthreshold，则表示需要将有效微震波形数据发送至云平台备份，将有效微震波形数据进行压缩以适应NB-IOT低带宽的场合，首先将微震波形数据处理数组的所有浮点数以逗号分隔首尾相连写成字符串，其次采用gzip算法对该字符串进行压缩得到压缩字符串；接着将微震特征向量、压缩字符串和布尔值标志位Flag编码成JSON格式的数据，然后用Hpack算法对其进行压缩得到压缩数据包；所述的布尔值标志位Flag，其值被设置为真，表示该压缩数据包含有有效微震波形数据和微震特征向量；
若所述微震特征向量中的最大幅值小于设定的阈值Vthreshold，则表示无需发送有效微震波形数据到云平台备份，只需发送微震特征向量到云平台保存；首先将微震特征向量和布尔值标志位Flag编码成JSON格式的数据，然后使用Hpack算法对其进行压缩得到压缩数据包；所述的布尔值标志位Flag，其值被设置为假，表示该压缩数据包仅含有微震特征向量。
进一步地，步骤3中，使用休眠唤醒算法对NB-IOT模块进行调度的具体过程为：
步骤3.1：系统初始化时，所述数据处理模块的休眠唤醒控制单元向NB-IOT模块按顺序发送激活AT指令、注网AT指令和休眠AT指令，使NB-IOT模块在被激活后进行注网操作，随后切换至低功耗模式，进入休眠状态；
步骤3.2：创建一个计数器Counter和一个计时器Timer，所述计数器Counter的溢出条件是其数值达到Cm，所述的计时器Timer的溢出条件是其时间达到Tm；
步骤3.3：重置计数器Counter和计时器Timer；休眠唤醒控制单元等待接收信号量，当其首次接收到信号量时，若信号量的布尔值标志位Flag为假，则计数器Counter的数值加1，且计时器Timer同步开启，然后进入步骤3.4，否则，进入步骤3.5；
步骤3.4：休眠唤醒控制单元继续等待接收信号量，每次接收到信号量时，若信号量的布尔值标志位Flag为真，则进入步骤3.5，否则，计数器Counter的数值加1，若计数器Counter满足溢出条件或计时器Timer满足溢出条件时，进入步骤3.5；否则重复执行步骤3.4；
步骤3.5：首先，所述休眠唤醒控制单元向NB-IOT模块发送唤醒AT指令以唤醒NB-IOT模块，接着，执行数据发送流程从数据发送缓冲队列中逐个取出压缩数据包，并发送压缩数据包到云平台，当数据发送缓冲队列为空时，所述休眠唤醒控制单元向NB-IOT模块发送休眠AT指令使其进入休眠状态，最后，进入步骤3.3。
一种基于NB-IOT的岩体破裂失稳微震信号无线监测装置，包括一个或者一个以上的数据采集模块、数据处理模块、NB-IOT模块和云平台模块，一个或者一个以上的数据采集模块与数据处理模块有线或者无线连接，数据处理模块经NB-IOT模块与云平台模块无线连接；
一个或者一个以上的数据采集模块用于获取岩体破裂所产生的微震信号，并将所述的微震信号转换成需要处理的数字化的电压信号，数据处理模块用于拾取有效微震波形数据，并从中提取微震特征向量，对微震特征向量和有效微震波形数据的压缩以满足NB-IOT低带宽的特点，同时对NB-IOT模块的休眠和唤醒状态之间的调度，NB-IOT模块用于实现与云平台模块建立通信链路，将压缩数据包发送至云平台模块，云平台模块用于与NB-IOT模块进行远程连接并建立通信链路，接收并存储来自所述的NB-IOT模块的数据，对数据进行分析，若满足预警条件则发送预警信息。
进一步地，数据采集模块包括加速度传感器单元、放大电路单元和A/D转换单元，加速度传感器单元经放大电路单元与A/D转换单元连接，加速度传感器单元用于采集岩体破裂失稳产生的微震信号，放大电路单元用于将所述加速度传感器单元采集到的微震信号进行放大，A/D转换单元用于将所述放大电路单元输出的微震信号转换为数字化的电压信号。
进一步地，数据处理模块包括有效波形拾取单元、特征提取单元、数据压缩单元和休眠唤醒控制单元，有效波形拾取单元的输入端与A/D转换单元连接，有效波形拾取单元的输出端与特征提取单元连接，特征提取单元的输出端与数据压缩单元连接，休眠唤醒控制单元与NB-IOT模块连接；
有效波形拾取单元用于使用长短时窗幅值比值法从所述A/D转换单元传输过来的数据中提取有效微震波形数据，特征提取单元用于从有效微震波形数据中提取微震特征向量，数据压缩单元用于压缩微震特征向量和有效微震波形数据，减少其所占的存储空间以便能够在NB-IOT的窄带宽环境下进行数据传输，休眠唤醒控制单元用于使用休眠唤醒算法对所述NB-IOT模块进行调度，使其在需要发送数据时被唤醒，而无需发送数据时休眠，节省电量以延长NB-IOT模块的工作时间。
本发明由于采用了上述技术方案，具有以下有益效果：
(1)本发明结合NB-IOT技术构建一个用于岩体破裂微震信号无线监测方法及装置，充分发挥NB-IOT技术的低成本、广覆盖、低功耗的优势来弥补现有微震信号监测系统工作时间过短、通信效率低下以及难以远距离传输的缺点，克服传统有线监测方式使用场景受限的难题和传统无线监测方式通讯效率低、延迟高且无法做到实时监测的难题，实现了微震监测数据的高效传输和实时监测。
(2)本发明提出的NB-IOT模块休眠唤醒算法，进一步延长了NB-IOT模块的工作时长，尤其适用于微震信号的短时高频率发送的场景，使得NB-IOT模块能够安装在偏远的山区地方而不用时常更换电池。
(3)本发明提出的微震数据压缩方法，使得压缩后的微震数据能够在NB-IOT的窄带宽环境下传输到云平台中，这有助于保存与岩体破裂失稳相关的、有分析价值的微震数据，为微震数据的后续处理和分析创造条件。
附图说明
图1是本发明实施例中基于NB-IOT的岩体破裂微震信号无线监测方法的流程图；
图2是本发明实施例中使用长短时窗幅值比值法拾取有效微震波形数据的流程图；
图3是本发明实施例中使用长短时窗幅值比值法拾取的有效微震波形数据的波形图；
图4是本发明实施例中使用休眠唤醒算法对NB-IOT模块进行调度的流程图；
图5是本发明实施例中的基于NB-IOT的岩体破裂微震信号无线监测装置的实现电路原理图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举出优选实施例，对本发明进一步详细说明。然而，需要说明的是，说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解，即便没有这些特定的细节也可以实现本发明的这些方面。
如图1-4所示，一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，图1为该方法的流程图，具体包括步骤101-105：
步骤101：在待监测的岩体不同位置安装数据采集模块，采集由岩体破裂失稳而产生的微震信号，使用长短时窗幅值比值法拾取有效微震波形数据。
在本发明实施例中，在待监测的岩体上安装数据采集模块，对待监测岩体进行实时监测；首先在待监测岩体表面钻出加速度传感器的安装孔，直径略大于加速度传感器直径，然后将加速度传感器深埋入安装孔内，并采用适当的耦合剂将待监测岩体与加速度传感器嵌为一体，耦合剂具有可塑性、速凝性、均一性等特点，将耦合界面内的空气、水分排除实现传感器与岩壁的直接接触。
在本发明实施例中，当岩体产生微观裂隙时会释放出低能量的弹性波并向四周传播，引发岩体的微小震动，数据采集模块的加速度传感器能够将岩体微小震动的加速度值转换为电压值，转换公式如式子(1)所示：
式中，d表示压电常数，m表示加速度传感器内部压电元件的质量，C表示压电元件两端的电容，在加速度传感器型号确定的情况下其d、m和C均为确定常量，因此电压值V与加速度a成正比，所采集的微震信号的电压值越大，表示岩体震动的加速度越大。
在本发明实施例中，数据采集模块的加速度传感器与数据采集模块的A/D转换单元使用导线相连接，使用A/D转换单元将加速度传感器输出的模拟信号转换为数字化的幅值并传输至数据处理模块作进一步处理。
在本发明实施例中，使用长短时窗幅值比值法拾取有效微震波形数据，具体地，图2为该方法的具体流程图，包括步骤A1～A10；图3为使用该方法拾取的有效微震波形数据的波形图：
A1：创建微震数据包暂存队列、有效波形暂存队列、短时窗数组和长时窗数组。
具体的，短时窗数组和长时窗数组的最大长度分别为M和N，优选的，M取为50，N取为1000；设有效信号提取阈值为R，优选的，R设置为1.5。
A2：初始化长时窗数组。
具体的，逐个读取从数据采集模块传输过来的电压幅值，将这些电压幅值按顺序保存在长时窗数组内，当长时窗数组的元素的个数达到N时，结束初始化过程。
A3：填充短时窗数组。
具体的，首先清空短时窗数组的所有元素，然后逐个读取从数据采集模块传输过来的电压幅值，将这些电压幅值按顺序保存在短时窗数组内，当短时窗数组的元素个数达到M时，结束填充短时窗数组的过程。
A4：计算短时窗数组和长时窗数组的平均电压幅值的比值r。
具体的，累加短时窗数组所有元素的值，然后除以M，得到短时窗数组的平均电压幅值Amplitude1；累加长时窗数组所有元素的值，然后除以N，得到长时窗数组的平均电压幅值Amplitude2；将长时窗数组的所有元素整体往前搬移M个位置，将短时窗数组内的元素依次添加到长时窗数组末尾的M个位置上；然后求得Amplitude1和Amplitude2的比值r＝Amplitude1/Amplitude2。
A5：判断比值r是否大于信号提取阈值R，并根据判断结果进行跳转。
具体的，若r≥R，则表示需要开始拾取有效微震波形，转步骤A6；若r＜R，则表示不需要拾取有效微震波形，转步骤A3；
A6：拾取有效微震波形数据。
具体的，将短时窗数组内的元素依次添加到有效波形暂存队列尾部，然后清空短时窗数组内的所有元素。
A7：填充短时窗数组。
具体的，逐个读取从数据采集模块传输过来的电压幅值，将这些电压幅值按顺序添加到短时窗数组内，当短时窗数组的元素达到M后，停止填充短时窗数组的流程。
A8：计算短时窗数组和长时窗数组的平均电压幅值的比值r。
具体的，累加短时窗数组所有元素的值，然后除以M，得到短时窗数组的平均电压幅值Amplitude1；累加长时窗数组所有元素的值，然后除以N，得到长时窗数组的平均电压幅值Amplitude2；将长时窗数组的所有元素整体往前搬移M个位置，将短时窗数组内的元素依次添加到长时窗数组末尾的M个位置上；然后求得Amplitude1和Amplitude2的比值r＝Amplitude1/Amplitude2。
A9：判断比值r是否大于等于R，并根据判断结果进行跳转。
具体的，若r≥R，则表示当前的微震波形仍属于有效波形，需要继续保持有效波形的拾取流程，转步骤A6；若r＜R，则表示结束拾取有效微震波形数据的流程，转步骤A10。
A10：结束拾取有效微震波形数据的流程。
具体的，将短时窗数组内的元素按顺序添加到有效波形暂存队列尾部中，获取当前的时间戳，将有效波形暂存队列的所有元素和时间戳打包成一个微震数据包，并将其放置于微震数据包暂存队列尾部，表示当前拾取了一个有效微震波形数据；接着清空有效波形暂存队列，然后转步骤A3开始拾取下一个有效微震波形数据。
步骤102：创建数据发送缓冲队列，从有效微震波形数据中提取微震特征向量，根据微震特征向量执行不同的压缩策略得到压缩数据包，将压缩数据包放置于数据发送缓冲队列尾部，同时发出一个携带了压缩数据包信息的信号量。
在本发明实施例中，从有效微震波形数据中提取微震特征向量；首先，创建一个微震波形数据处理数组；其次，从所述的微震数据包暂存队列的头部取出一个微震数据包；接着，将其中的有效微震波形数据逐元素复制到微震波形数据处理数组中；然后，从微震波形数据处理数组的第一个元素开始遍历，到最后一个元素结束遍历，计算得到触发位置、结束位置、信号持续时间、上升时间、最大幅值和最大幅值的位置共六个特征值；最后，将这六个特征值组成微震特征向量。
在本发明实施例中，微震特征向量中的触发位置表示数组中第一个大于阈值的幅值对应的数组下标；结束位置表示数组中最小幅值对应的数组下标；信号持续时间表示结束位置与触发位置的差值；最大幅值表示数组所有幅值中最大的幅值；最大幅值的位置表示数组中最大幅值所对应的数组下标；上升时间表示最大幅值的位置与触发位置的差值。
在本发明实施例中，判断微震特征向量中的最大幅值是否大于设定的阈值Vthreshold，如果大于或等于该阈值，则表示同时需要将有效微震波形数据发送至云端备份，将有效微震波形数据进行压缩以适应NB-IOT低带宽的场合：首先将微震波形数据处理数组的所有浮点数以逗号相隔首尾相连写成字符串，其次采用gzip算法对该字符串进行压缩得到压缩字符串；接着将微震特征向量、压缩字符串和布尔值标志位Flag编码成JSON格式的数据，然后采用Hpack算法对该JSON格式的数据进行压缩，得到含有有效微震波形数据和微震特征向量的压缩数据包；所述的布尔值标志位Flag，其值被设置为真，表示该压缩数据包含有有效微震波形数据和微震特征向量；
如果微震特征向量中的最大幅值小于设定的阈值Vthreshold，那么就无需发送有效微震波形数据到云端备份，只需发送微震特征向量到云平台保存；首先将微震特征向量和布尔值标志位Flag编码成JSON格式的数据包，然后使用Hpack算法对该JSON格式的数据包进行压缩，得到仅含有微震特征向量的压缩数据包；所述的布尔值标志位Flag，其值被设置为假，表示该压缩数据包仅含有微震特征向量。
步骤103：为了延长NB-IOT模块的工作时间并且能够满足实时发送微震数据的需求，使用休眠唤醒算法对NB-IOT模块进行调度。
在本发明实施例中，使用休眠唤醒算法对NB-IOT模块进行调度，具体地，图4为该方法的具体流程图，包括步骤B1～B11：
B1：NB-IOT模块初始化。
具体的，数据处理模块的休眠唤醒控制单元通过串口向NB-IOT模块发送激活AT指令，该AT指令的形式为：“AT+SM＝LOCK”；接着发送注网AT指令，该AT指令的形式为：“AT+CEREG＝1”；然后发送休眠AT指令，该AT指令的形式为：“AT+CFUN＝0”。
B2：分别创建一个计数器Counter和一个计时器Timer。
具体的，计数器Counter的溢出条件是其数值达到Cm；而计时器Timer的溢出条件是其时间达到Tm。
B3：重置计数器和计时器。
具体的，计数器Counter的数值设置为0；计时器Timer的时间设置为0秒，每过一秒，时间都会自增一秒。
B4：休眠唤醒控制单元等待接收信号量。
B5：判断信号量的布尔值标志位Flag是否为真，并根据判断结果进行跳转。
具体的，若布尔值标志位Flag为假，说明产生了一个仅含有微震特征向量的压缩数据包，该压缩数据包没有必要马上发送到云平台分析，所以此时暂时不唤醒NB-IOT模块，而是利用批处理的思想来等待更多的仅含有微震特征向量的压缩数据包的到来，再集中发送数据到云平台，转步骤B6；若布尔值标志位Flag为真，说明产生了一个包含有效微震波形数据和微震特征向量的压缩数据包，该压缩数据包需要马上发送至云平台进行分析，转步骤B11。
B6：计数器的数值加1，同时开启计时器。
B7：休眠唤醒控制单元等待接收信号量。
B8：判断信号量的布尔值标志位Flag是否为真，并根据判断结果进行跳转。
具体的，若布尔值标志位Flag为真，则转步骤B11；若布尔值标志位Flag为假，转步骤B9。
B9：计数器的数值加1。
B10：判断计数器是否满足溢出条件或计时器是否满足溢出条件，并根据判断结果进行跳转。
具体的，若计数器Counter的数值等于Cm或计时器Timer的时间等于Tm，则表示需要批量发送目前数据发送缓冲队列存在的压缩数据包到云平台，转步骤B11；否则，转步骤B7。
B11：唤醒NB-IOT模块，执行数据发送流程来发送数据。
具体的，首先，所述休眠唤醒控制单元向NB-IOT模块发送唤醒AT指令以唤醒NB-IOT模块；接着，执行数据发送流程从数据发送缓冲队列中逐个取出压缩数据包，并发送压缩数据包到云平台，当数据发送缓冲队列为空时，所述休眠唤醒控制单元向NB-IOT模块发送休眠AT指令使其进入休眠状态；最后，进入步骤B3。
步骤104：在NB-IOT模块处于唤醒状态时，与云平台模块建立通信链路，将所述压缩数据包发送至云平台。
在本发明实施例中，步骤104具体包括了子步骤C1和C2：
C1：休眠唤醒控制单元通过串口向NB-IOT模块发送连接云平台AT指令，该指令为“AT+QMTOPEN＝0,"iot-as-mqtt.cn-shanghai.aliyuncs.com",1883”，指令包含了云平台的网络域名和端口号；其次，当NB-IOT模块接收到该指令后，会向云平台发起一个基于MQTT协议的连接请求，云平台从请求中取出NB-IOT的身份信息和密钥进行验证，若验证通过则会与NB-IOT模块建立一个稳定的MQTT连接；最后，当NB-IOT模块与云平台建立通信链路后，便可以开始执行数据发送流程将数据发送到云平台；
C2：从数据发送缓冲队列的头部取出一个压缩数据包，将其发送至通信链路中，当接收到云平台返回的携带成功标志位的响应数据后，表示云平台已经成功接收到压缩数据包；然后从数据发送缓冲队列内取出下一份压缩数据包，重复发送过程，直至数据发送缓冲队列为空。
步骤105：云平台实时接收NB-IOT模块发送的数据包，对其进行分析预警操作，然后将数据保存到数据库中。
在本发明实施例中，步骤105包括子步骤D1和D2：
D1：云平台开启通信端口监听，若接收到来自NB-IOT模块的握手协议数据包，则对其鉴权信息进行鉴权，若通过则返回携带成功标志位的响应数据包并打开通信链路；否则，返回携带失败标志位的响应数据包；当通信链路打开后，连续地读取由NB-IOT模块发送过来的数据包，每读取到一个数据包后，将数据包放进队列中等待处理并向NB-IOT模块返回携带成功标志位的响应数据包；
D2：连续地从队列中取出数据包，并对每一个数据包进行分析、预警和存储操作；
首先，使用Hpack算法解压缩得到布尔值标志位Flag和微震特征向量；读取布尔值标志位Flag，若其值为真，说明解压缩还会得到压缩字符串，对其使用gzip算法解压缩得到有效微震波形数据；将微震特征向量和有效微震波形数据都保存到数据库中，同时云平台模块发出一个预警信息推送到管理员的手机中，该预警信息携带了该有效微震波形数据及时间戳、监测岩体所在的位置等；该预警信息用于辅助管理员对监测岩体的健康状态做出判断。
若布尔值标志位Flag的值为假，说明该压缩数据包只有微震特征向量，而不含有有效微震波形数据，因此将有微震特征向量保存到数据库中即可。
以上详细介绍了一种基于NB-IOT的岩体破裂失稳微震信号无线监测的方法流程，该方法也可以通过相应的装置实现，下面详细介绍该装置的结构和功能。
本发明实施例提供一种基于NB-IOT的岩体破裂失稳微震信号无线监测的装置，包括：
数据采集模块，用于获取岩体破裂所产生的微震信号，并将所述的微震信号转换成需要处理的数字化的电压信号；
数据处理模块，用于拾取有效微震波形数据，并从中提取微震特征向量，用于实现对微震特征向量和有效微震波形数据的压缩以满足NB-IOT低带宽的特点，还用于实现对NB-IOT模块的休眠和唤醒状态之间的调度；
NB-IOT模块，用于实现与云平台模块建立通信链路，将压缩数据包发送至云平台模块；
云平台模块，用于与NB-IOT模块进行远程连接并建立通信链路，接收并存储来自所述的NB-IOT模块的数据，对数据进行分析，若满足预警条件则发送预警信息；
在本发明实施例中，所述数据采集模块包括：
加速度传感器单元，用于采集岩体破裂失稳产生的微震信号；
放大电路单元，用于将所述加速度传感器单元采集到的微震信号进行放大；
A/D转换单元，用于将所述放大电路单元输出的微震信号转换为数字化的电压信号；
在本发明实施例中，所述数据处理模块包括：
有效波形拾取单元，用于使用长短时窗幅值比值法从所述A/D转换单元传输过来的数据中提取有效微震波形数据；
特征提取单元，用于从有效微震波形数据中提取微震特征向量；
数据压缩单元，用于压缩微震特征向量和有效微震波形数据，减少其所占的存储空间以便能够在NB-IOT的窄带宽环境下进行数据传输；
休眠唤醒控制单元，用于使用休眠唤醒算法对所述NB-IOT模块进行调度，使其在需要发送数据时被唤醒，而无需发送数据时休眠，节省电量以延长NB-IOT模块的工作时间。
本发明实施例提供一种基于NB-IOT的岩体破裂失稳微震信号无线监测装置的硬件实现电路，该电路实现了所述的数据采集模块、数据处理模块和NB-IOT无线传输模块，请参见图5，具体为：
所述数据采集模块包括微震传感器SENSOR1，其VCC引脚接电源，GND引脚接地，信号输出端OUT接电容C1的一端，电容C1的另一端接电阻R1的一端，电阻R1的另一端接运放器AR1的同相输入端，电阻R2的一端接运放器AR1的同相输入端，电阻R2的另一端接电容C2的一端，电容C2的另一端接电源，电阻R3和电阻R4的一端接运放器AR1的反相输入端，电阻R3的另一端接地，电阻R4的另一端接运放器的输出端，电阻R5接运放器AR1的输出端，电阻R5的另一端接运放器AR2的同相输入端，运放器AR2的反相输入端接其输出端，运放器AR2的输出端接A/D转换单元的IN0引脚；或非门NOR1的输出端接ADC1的OE引脚，或非门NOR2的输出端接ADC1的START引脚和ALE引脚，ADC1的VCC引脚和VREF(+)引脚接电源，ADC1的GND引脚和VREF(-)引脚接地。优选的，A/D转换单元可以使用ADC0809芯片。
所述数据处理模块包括微型处理器MCU1，其VCC引脚接电源，GND引脚接地，MCU1的PC0至PC7引脚分别与所述的A/D转换单元ADC1的2-1至2-8引脚相连接，MCU1的PB0引脚接ADC1的或非门NOR1的一个输入端，MCU1的PB1引脚接ADC1的或非门NOR1的另一个输入端，MCU1的PB2引脚接ADC1的或非门NOR2的一个输入端，MCU1的PB3引脚、PB4引脚和PB5引脚分别接ADC1的ADDA引脚、ADDB引脚和ADDC引脚，MCU1的PB6引脚与ADC1的CLOCK引脚相连，MCU1的PB7引脚与ADC1的EOC引脚相连。优选的，微型处理器MCU1可以采用STM32F1系列的芯片。
所述NB-IOT模块包括NB-IOT通信模块NB1，NB1的VBAT引脚接电容C3的一端，NB1的GND引脚接电容C3的另一端，电容C4、电容C5和稳压二极管D1与电容C3并联，稳压二极管D1的一端接电源，稳压二极管的另一端接地；NB1的RESET引脚接稳压二极管D2的一端，稳压二极管D2的另一端接地，NB1的RESET引脚还接了按键开关KEY1的一端，按键开关KEY1的另一端接地；NB1的SIM_GND引脚接电容C6的一端，电容C6的另一端接SIM卡连接器SCC的VCC引脚，NB1的SIM_GND引脚还接了SIM卡连接器SCC的GND引脚，NB1的SIM_VDD引脚接了SIM卡连接器SCC的VCC引脚，NB1的SIM_VDD引脚还接了电阻R9的一端，电阻R9的另一端接了NB1的SIM_DATA引脚，电阻R9的另一端还接了电阻R8的一端，电阻R8另一端还接了SIM卡连接器SCC的IO引脚，NB1的SIM_RST引脚接了电阻R6的一端，电阻R6的另一端还接了SIM卡连接器SCC的RST引脚，NB1的SIM_CLK引脚接电阻R7的一端，电阻R7的另一端接SIM卡连接器SCC的CLK引脚；NB1的NETLIGHT引脚接电阻R14的一端，电阻R14的另一端接晶体管Q3的基极，晶体管Q3的发射极接电阻R15的一端，电阻R15的另一端接地以及晶体管Q3的基极，晶体管Q3的集电极接电阻R16的一端，电阻R16的另一端接发光二极管LED1的一端，发光二极管LED1的另一端接NB1的VDD引脚；NB1的TXD引脚接晶体管Q1的发射极，晶体管Q1的集电极接MCU1的PA1引脚，晶体管Q1的集电极还接了电阻R10的一端，电阻R10的另一端接MCU1的VDD引脚，晶体管Q1的基极接电阻R11的一端，电阻R11的另一端接NB1的VDD_EXT引脚，晶体管Q1的基极还接了电容C7的一端，电容C7的另一端接NB1的VDD_EXT引脚，NB1的RXD引脚接电阻R12的一端，电阻R12的另一端接NB1的VDD_EXT引脚，电阻R12的另一端还接了电阻R13和电容C8的一端，电阻R13和电容C8的另一端接晶体管Q2的基极，NB1的RXD引脚还接了晶体管Q2的集电极，晶体管Q2的发射极接MCU1的PA0引脚。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的方法可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
此外，在本申请各个实施例中的各个单元可以集成在一个模块中，也可以是单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (9)
1.一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，其特征在于，所述方法包括如下步骤：
步骤1：在待监测岩体上安装一个或者一个以上的数据采集模块，采集岩体破裂产生的微震信号，使用长短时窗幅值比值法拾取有效微震波形数据；
步骤2：创建数据发送缓冲队列，从有效微震波形数据中提取微震特征向量，根据微震特征向量执行不同的压缩策略得到压缩数据包，将压缩数据包放置于数据发送缓冲队列尾部，同时发出一个携带压缩数据包信息的信号量；
步骤3：使用休眠唤醒算法对NB-IOT模块进行调度，使其在需要发送数据的时候被唤醒，而不需要发送数据的时候进行休眠，节省电量以延长工作时间；
步骤4：当NB-IOT模块处于唤醒状态时，与云平台模块建立通信链路，将所述的压缩数据包发送至云平台模块；
步骤5：云平台模块对接收到的数据包进行分析，若满足预警条件则发出预警信息。
2.根据权利要求1所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，其特征在于，步骤1中，采集岩体破裂产生的微震信号的具体过程为：
步骤1.1.1：当岩体产生微观裂隙时会释放出低能量的弹性波并向四周传播，引发岩体的微小震动，数据采集模块的加速度传感器能够将岩体微小震动的加速度值转换为电压值得到微震信号，转换公式如式子(1)所示：
式中，d表示压电常数，m表示加速度传感器内部压电元件的质量，C表示压电元件两端的电容，在加速度传感器型号确定的情况下其d、m和C均为确定常量，因此电压值V与加速度a成正比，所采集的微震信号的电压值越大，表示岩体震动的加速度越大；
步骤1.1.2：使用所述数据采集模块的A/D转换单元将加速度传感器输出的模拟信号转换为数字化的幅值并传输至数据处理模块作进一步处理。
3.根据权利要求1所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，其特征在于，步骤1中，使用长短时窗幅值比值法拾取有效微震波形数据的具体过程为：
步骤1.2.1：创建微震数据包暂存队列、有效波形暂存队列、短时窗数组和长时窗数组，所述的短时窗数组和长时窗数组的最大长度分别为M和N，设有效信号提取阈值为R；所述的M、N和R根据监测现场环境决定；
步骤1.2.2：初始化长时窗数组，逐个读取从所述数据采集模块传输过来的幅值，将其按顺序保存于长时窗数组内，直至长时窗数组的元素个数达到N；
步骤1.2.3：逐个读取从所述数据采集模块传输过来的幅值，将其按顺序保存于短时窗数组内，当短时窗数组内的元素达到M时，分别计算短时窗数组和长时窗数组内的平均电压幅值Amplitude1和Amplitude2，接着将长时窗数组的元素整体往前搬移M个位置，将短时窗数组内的元素依次添加到长时窗数组末尾的M个位置上，然后求得Amplitude1和Amplitude2的比值r＝Amplitude1/Amplitude2，若r大于等于R则表示开始拾取有效微震波形数据，接着转步骤1.2.4：；否则，清空短时窗数组内的元素并重复执行步骤1.2.3；
步骤1.2.4：将短时窗数组内的元素依次添加到有效波形暂存队列尾部中；接着清空短时窗数组内的元素，然后继续逐个读取从所述数据采集模块传输过来的幅值，当短时窗数组内的元素达到M后，分别计算短时窗数组和长时窗数组的平均电压幅值Amplitude1和Amplitude2并求得二者的比值r＝Amplitude1/Amplitude2，若r小于R则表示结束拾取有效的微震波形数据，转步骤1.2.5，否则，将长时窗数组的元素整体往前搬移M个位置，将短时窗数组内的元素依次添加到长时窗数组末尾的M个位置上，重复执行步骤1.2.4；
步骤1.2.5：将短时窗数组内的元素按顺序添加到有效波形暂存队列尾部中，获取当前的时间戳，将有效波形暂存队列的所有元素和时间戳打包成一个微震数据包，并将其放置于微震数据包暂存队列尾部，表示当前拾取了一个有效微震波形数据，接着清空有效波形暂存队列，然后重复执行步骤1.2.3，开始拾取下一个有效微震波形数据。
4.根据权利要求1所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，其特征在于，步骤2中，从有效微震波形数据中提取微震特征向量的具体过程为：
创建一个微震波形数据处理数组，从所述的微震数据包暂存队列的头部中取出一个微震数据包，将其中的有效微震波形数据逐元素复制到微震波形数据处理数组中；从微震波形数据处理数组的第一个元素开始遍历，到最后一个元素结束遍历，计算得到触发位置、结束位置、信号持续时间、上升时间、最大幅值和最大幅值的位置共六个特征值，将这六个特征值组成微震特征向量，六个特征值的具体定义为：
触发位置：表示数组中第一个大于阈值的幅值对应的数组下标；
结束位置：表示数组中最小幅值对应的数组下标；
信号持续时间：表示结束位置与触发位置的差值；
最大幅值：表示数组所有幅值中最大的幅值；
最大幅值的位置：表示数组中最大幅值所对应的数组下标；
上升时间：表示最大幅值的位置与触发位置的差值。
5.根据权利要求1所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，其特征在于，步骤2中，根据微震特征向量执行不同的压缩策略得到压缩数据包的具体过程为：
若所述微震特征向量的最大幅值大于或等于设定的阈值Vthreshold，则表示需要将有效微震波形数据发送至云平台备份，将有效微震波形数据进行压缩以适应NB-IOT低带宽的场合，首先将微震波形数据处理数组的所有浮点数以逗号分隔首尾相连写成字符串，其次采用gzip算法对该字符串进行压缩得到压缩字符串；接着将微震特征向量、压缩字符串和布尔值标志位Flag编码成JSON格式的数据，然后用Hpack算法对其进行压缩得到压缩数据包；所述的布尔值标志位Flag，其值被设置为真，表示该压缩数据包含有有效微震波形数据和微震特征向量；
若所述微震特征向量中的最大幅值小于设定的阈值Vthreshold，则表示无需发送有效微震波形数据到云平台备份，只需发送微震特征向量到云平台保存；首先将微震特征向量和布尔值标志位Flag编码成JSON格式的数据，然后使用Hpack算法对其进行压缩得到压缩数据包；所述的布尔值标志位Flag，其值被设置为假，表示该压缩数据包仅含有微震特征向量。
6.根据权利要求1所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测方法，其特征在于，步骤3中，使用休眠唤醒算法对NB-IOT模块进行调度的具体过程为：
步骤3.1：系统初始化时，所述数据处理模块的休眠唤醒控制单元向NB-IOT模块按顺序发送激活AT指令、注网AT指令和休眠AT指令，使NB-IOT模块在被激活后进行注网操作，随后切换至低功耗模式，进入休眠状态；
步骤3.2：创建一个计数器Counter和一个计时器Timer，所述计数器Counter的溢出条件是其数值达到Cm，所述的计时器Timer的溢出条件是其时间达到Tm；
步骤3.3：重置计数器Counter和计时器Timer；休眠唤醒控制单元等待接收信号量，当其首次接收到信号量时，若信号量的布尔值标志位Flag为假，则计数器Counter的数值加1，且计时器Timer同步开启，然后进入步骤3.4，否则，进入步骤3.5；
步骤3.4：休眠唤醒控制单元继续等待接收信号量，每次接收到信号量时，若信号量的布尔值标志位Flag为真，则进入步骤3.5，否则，计数器Counter的数值加1，若计数器Counter满足溢出条件或计时器Timer满足溢出条件时，进入步骤3.5；否则重复执行步骤3.4；
步骤3.5：首先，所述休眠唤醒控制单元向NB-IOT模块发送唤醒AT指令以唤醒NB-IOT模块，接着，执行数据发送流程从数据发送缓冲队列中逐个取出压缩数据包，并发送压缩数据包到云平台，当数据发送缓冲队列为空时，所述休眠唤醒控制单元向NB-IOT模块发送休眠AT指令使其进入休眠状态，最后，进入步骤3.3。
7.一种基于NB-IOT的岩体破裂失稳微震信号无线监测装置，其特征在于：包括一个或者一个以上的数据采集模块、数据处理模块、NB-IOT模块和云平台模块，一个或者一个以上的数据采集模块与数据处理模块有线或者无线连接，数据处理模块经NB-IOT模块与云平台模块无线连接；
一个或者一个以上的数据采集模块用于获取岩体破裂所产生的微震信号，并将所述的微震信号转换成需要处理的数字化的电压信号，数据处理模块用于拾取有效微震波形数据，并从中提取微震特征向量，对微震特征向量和有效微震波形数据的压缩以满足NB-IOT低带宽的特点，同时对NB-IOT模块的休眠和唤醒状态之间的调度，NB-IOT模块用于实现与云平台模块建立通信链路，将压缩数据包发送至云平台模块，云平台模块用于与NB-IOT模块进行远程连接并建立通信链路，接收并存储来自所述的NB-IOT模块的数据，对数据进行分析，若满足预警条件则发送预警信息。
8.根据权利要求7所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测装置，其特征在于：数据采集模块包括加速度传感器单元、放大电路单元和A/D转换单元，加速度传感器单元经放大电路单元与A/D转换单元连接，加速度传感器单元用于采集岩体破裂失稳产生的微震信号，放大电路单元用于将所述加速度传感器单元采集到的微震信号进行放大，A/D转换单元用于将所述放大电路单元输出的微震信号转换为数字化的电压信号。
9.根据权利要求8所述的一种基于NB-IOT的岩体破裂失稳微震信号无线监测装置，其特征在于：数据处理模块包括有效波形拾取单元、特征提取单元、数据压缩单元和休眠唤醒控制单元，有效波形拾取单元的输入端与A/D转换单元连接，有效波形拾取单元的输出端与特征提取单元连接，特征提取单元的输出端与数据压缩单元连接，休眠唤醒控制单元与NB-IOT模块连接；
有效波形拾取单元用于使用长短时窗幅值比值法从所述A/D转换单元传输过来的数据中提取有效微震波形数据，特征提取单元用于从有效微震波形数据中提取微震特征向量，数据压缩单元用于压缩微震特征向量和有效微震波形数据，减少其所占的存储空间以便能够在NB-IOT的窄带宽环境下进行数据传输，休眠唤醒控制单元用于使用休眠唤醒算法对所述NB-IOT模块进行调度，使其在需要发送数据时被唤醒，而无需发送数据时休眠，节省电量以延长NB-IOT模块的工作时间。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111666298.6A CN114325823B (zh) | 2021-12-30 | 2021-12-30 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
| PCT/CN2022/095165 WO2023123844A1 (zh) | 2021-12-30 | 2022-05-26 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111666298.6A CN114325823B (zh) | 2021-12-30 | 2021-12-30 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114325823A true CN114325823A (zh) | 2022-04-12 |
| CN114325823B CN114325823B (zh) | 2022-07-15 |
Family
ID=81021631
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111666298.6A Active CN114325823B (zh) | 2021-12-30 | 2021-12-30 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN114325823B (zh) |
| WO (1) | WO2023123844A1 (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023123844A1 (zh) * | 2021-12-30 | 2023-07-06 | 广西大学 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117289344B (zh) * | 2023-11-24 | 2024-01-30 | 北京科技大学 | 基于震源空间分布的煤岩失稳破坏快速判识方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130336091A1 (en) * | 2012-06-18 | 2013-12-19 | Halliburton Energy Services, Inc. | Statistics-based seismic event detection |
| CN103513293A (zh) * | 2013-10-12 | 2014-01-15 | 广西大学 | 隧道地质综合超前预报专家系统及其实现方法 |
| CN103913772A (zh) * | 2014-04-02 | 2014-07-09 | 西南石油大学 | 基于储层地质力学参数的微地震事件正演模拟方法 |
| CN105956526A (zh) * | 2016-04-22 | 2016-09-21 | 山东科技大学 | 基于多尺度排列熵的低信噪比微震事件辨识方法 |
| CN107561579A (zh) * | 2017-08-31 | 2018-01-09 | 北京市政建设集团有限责任公司 | 一种隧道施工微震监测系统及监测方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105487114B (zh) * | 2015-12-08 | 2018-02-02 | 中南大学 | 一种微震信号p波初至点综合拾取方法 |
| EP3361289A1 (en) * | 2017-02-08 | 2018-08-15 | Shell International Research Maatschappij B.V. | Method, sensor and system for wireless seismic networking |
| CN208125931U (zh) * | 2018-04-10 | 2018-11-20 | 成都东方监控技术有限公司 | 震灾信息监测应急器 |
| CN111083677B (zh) * | 2019-12-31 | 2023-04-28 | 山东普赛通信科技股份有限公司 | 一种低功耗下行通信系统及方法 |
| CN112130175A (zh) * | 2020-09-21 | 2020-12-25 | 中国地质环境监测院 | 一种地质监测系统及方法 |
| CN112731518B (zh) * | 2020-12-09 | 2022-06-07 | 中国科学院地质与地球物理研究所 | 一种人工智能实时微地震监测节点 |
| CN113050159B (zh) * | 2021-03-23 | 2021-11-16 | 中国矿业大学 | 一种煤岩水力压裂裂缝微震定位及扩展机理监测方法 |
| CN114325823B (zh) * | 2021-12-30 | 2022-07-15 | 广西大学 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
-
2021
- 2021-12-30 CN CN202111666298.6A patent/CN114325823B/zh active Active
-
2022
- 2022-05-26 WO PCT/CN2022/095165 patent/WO2023123844A1/zh unknown
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130336091A1 (en) * | 2012-06-18 | 2013-12-19 | Halliburton Energy Services, Inc. | Statistics-based seismic event detection |
| CN103513293A (zh) * | 2013-10-12 | 2014-01-15 | 广西大学 | 隧道地质综合超前预报专家系统及其实现方法 |
| CN103913772A (zh) * | 2014-04-02 | 2014-07-09 | 西南石油大学 | 基于储层地质力学参数的微地震事件正演模拟方法 |
| CN105956526A (zh) * | 2016-04-22 | 2016-09-21 | 山东科技大学 | 基于多尺度排列熵的低信噪比微震事件辨识方法 |
| CN107561579A (zh) * | 2017-08-31 | 2018-01-09 | 北京市政建设集团有限责任公司 | 一种隧道施工微震监测系统及监测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 刘玉春，等: "含水煤岩单轴压缩微震信号特征试验研究", 《中国安全生产科学技术》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023123844A1 (zh) * | 2021-12-30 | 2023-07-06 | 广西大学 | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114325823B (zh) | 2022-07-15 |
| WO2023123844A1 (zh) | 2023-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114325823B (zh) | 一种基于nb-iot的岩体破裂失稳微震信号无线监测方法与装置 | |
| CN104281247B (zh) | 执行唤醒控制的方法及装置 | |
| CN101470013B (zh) | 一种架空输电线路舞动监测方法及装置 | |
| CN106246476B (zh) | 海上风力发电机叶片故障智能监测系统及方法 | |
| CN109555563A (zh) | 一种煤岩动力灾害危险的电磁辐射和地音监测预警方法 | |
| CN201181992Y (zh) | 基于移动终端的实时温湿度监测与控制装置 | |
| CN103052146B (zh) | 一种低功耗埋入式无线传感器网络节点装置 | |
| CN201134152Y (zh) | 一种架空输电线路舞动监测装置及系统 | |
| CN102129752A (zh) | 室外设施人为破坏在线监测报警系统 | |
| CN102466811A (zh) | 一种微地震监测系统和方法 | |
| CN108761522A (zh) | 搭载于隧道掘进机上的地震波超前探测系统及方法 | |
| CN111862569B (zh) | 一种智能局放监测数据传输的算法 | |
| D’Errico et al. | Structural health monitoring and earthquake early warning on 5g urllc network | |
| CN104295316A (zh) | 一种地音监测方法与预警系统 | |
| CN102506988A (zh) | 输电线路在线舞动及风偏监测系统及其监测方法 | |
| CN104554086B (zh) | 一种能够累计工程机械工作时间的车辆追踪器 | |
| CN208140639U (zh) | 一种基于声发射技术的风机叶片监测装置 | |
| CN104730150A (zh) | 多通道声发射监测系统 | |
| CN207742067U (zh) | 一种水利工程渗流实时监测系统 | |
| CN109443521A (zh) | 一种矿用地音传感器 | |
| CN205562823U (zh) | 一种岩爆预测系统 | |
| CN204904215U (zh) | 一种多功能无线空中鼠标系统 | |
| CN204630657U (zh) | 一种电缆管道振动监测节点 | |
| CN211654300U (zh) | 一种矿井下安全生产语音合成装置 | |
| CN208969247U (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 |