井地一体微震监测系统监测单元采集数据授时对齐方法 |
|||||||
申请号 | CN202210586305.X | 申请日 | 2022-05-27 | 公开(公告)号 | CN114895353B | 公开(公告)日 | 2023-03-10 |
申请人 | 中国矿业大学; 徐州弘毅科技发展有限公司; | 发明人 | 巩思园; 葛庆; 窦林名; | ||||
摘要 | 一种井地一体微震监测系统监测单元采集数据授时对齐方法,利用井地一体微震监测系统获得地面及井下监测单元采集的带有时间戳的矿震 波形 数据段,分别计算两段波形数据中的所有相邻GPS授时时间戳的差值和相邻GPS授时时间戳内的 采样 点数,等间距补齐所有相邻GPS授时时间戳内未打时间戳的采样点的采样时刻,计算两段矿震波形数据的平均采样 频率 ,分别等间距补齐两段采样数据中首、尾GPS授时时间戳之外的采样点采样时刻;获得所有采样点上的采样时刻后,按照固定 采样频率 对采样时刻重新取样,采用线性插值公式计算新采样时刻下的波形数据,进而对齐两段波形数据的采样时刻。本 发明 能准确对齐地面监测单元和井下监测单元采集的震动 波数 据。 | ||||||
权利要求 | 1.一种井地一体微震监测系统监测单元采集数据授时对齐方法,其特征在于,包括如下步骤: |
||||||
说明书全文 | 井地一体微震监测系统监测单元采集数据授时对齐方法技术领域背景技术[0002] 煤炭是我国主体能源,其安全开采至关重要,近年来随着许多煤矿进入深部开采,由冲击地压灾害带来的安全问题日益突出。目前,微震监测技术作为冲击地压危险分析和预警所依赖的重要技术,在国内煤矿已获得广泛应用,取得了较显著的监测预警效果。但受限于井下近水平巷道的传感器安装条件,使得单一使用井下布设传感器的方式无法突破大幅提高微震监测精度的瓶颈。 [0003] 为达到对冲击地压精准探测的目的,通过在地面安装一定数量的传感器,形成井地联合观测模式,可在空间上立体包围监测区域。为方便移动地面监测点和确保井下监测数据质量,地面监测单元采用无线通讯方式,井下监测单元采用有线通讯方式,两种监测方式下分别采用独立的GPS进行采样数据的非连续授时,同时,监测单元也存在采样频率不稳定的情况,导致汇集到监测中心的两类监测单元的授时数据之间存在授时时间无法对齐的问题,极大限制了微震监测精度的提升。 发明内容[0004] 针对上述现有技术存在的问题,本发明提供一种井地一体微震监测系统监测单元采集数据授时对齐方法,该方法能够准确对齐井地一体微震监测系统中地面监测单元和井下监测单元采集的震动波数据,保证后续矿震定位和能量求解的准确性,提升微震监测精度。 [0005] 为了实现上述目的,本发明采用的技术方案是:一种井地一体微震监测系统监测单元采集数据授时对齐方法,包括如下步骤: [0006] (1)分别取由井地一体微震监测系统地面无线监测单元和井下有线监测单元独立采集的带有多个GPS授时时间戳的矿震波形数据一和矿震波形数据二;矿震波形数据一和矿震波形数据二为长度不同的两段数据段; [0007] (2)分别计算矿震波形数据一和矿震波形数据二中所有相邻GPS授时时间戳之间的差值以及相邻GPS授时时间戳内包含的采样点数后,等间距补齐相邻GPS授时时间戳内未打时间戳的采样点的采样时刻; [0008] (3)分别计算矿震波形数据一和矿震波形数据二的平均采样频率,分别等间距补齐矿震波形数据一和矿震波形数据二中的首、尾GPS授时时间戳之外的采样点的采样时刻; [0009] (4)根据步骤(2)和步骤(3)得到的所有采样点上的采样时刻,按照统一设置的采样频率对采样时刻进行重新取样; [0010] (5)根据步骤(2)和步骤(3)得到的所有采样点上的采样时刻和波形数据,采用线性插值公式计算步骤(4)中新采样时刻下的矿震波形数据,对齐矿震波形数据一和矿震波形数据二的采样时刻。 [0011] 进一步地,所述步骤(1)中,地面无线监测单元和井下有线监测单元分别采用独立的高精度GPS进行授时,其采样频率分别为fs和fu;所取的矿震波形数据一和矿震波形数据二分别为 和 其中,i=1,2,...,l,l为地面监测单元的采样长度;j=1,2,...,p,p为井下监测单元的采样长度;所述GPS授时时间戳为两个独立的GPS分别对矿震波形数据一和矿震波形数据二上的部分采样点打时间戳,其中, 的部分采样点 上的时间戳形成′ sGPS授时序列为 其中,i =1,2,...,m;index为在 中有时间戳的位置序号,m为地面无线监测单元所监测的矿震波形数据一中带有GPS时间戳的采样点个数; 的部分采样点′ u 上的时间戳形成的GPS授时序列为 其中,j =1,2,...,n;index 为在 中有时间戳的位置序号,n为井下有线监测单元所监测的矿震波形数据二中带有GPS时间戳的采样点个数。 [0012] 进一步地,所述步骤(2)中,所述等间距补齐相邻GPS授时时间戳内未打时间戳的采样点的采样时刻的步骤为: [0013] 201)分别计算GPS授时序列 中相邻时间戳的差值和其内包含的采样点数分别为: [0014] 时间差: [0015] 采样点数: [0016] 202)利用时间差和采样点数,采用公式: [0017] [0018] 等间距计算 上采样点 之后第k个点的采样时刻,直至计算至采样点停止计算; [0019] 采用公式: [0020] [0021] 等间距计算 上采样点 之后第k个点的采样时刻,直至计算至采样点停止计算; [0022] 203)重复步骤201)和202),直至补齐矿震波形数据一和矿震波形数据二中所有相邻GPS授时时间戳内未打时间戳的采样点的采样时刻。 [0023] 进一步地,所述步骤(3)中,等间距补齐矿震波形数据一和矿震波形数据二中的首、尾GPS授时时间戳之外的采样点的采样时刻的步骤为: [0024] 301)分别计算矿震波形数据一和矿震波形数据二的平均采样频率 和 为: [0025] [0026] [0027] 302)采用算式 等间距补齐地面无线监测单元所监测的矿震波形数据一首个GPS授时时间戳之前第k个采样点的采样时刻,直至补充至矿震波形数据一的第1个采样点,停止补充;采用算式 等间距补齐地面无线监测单元所监测的矿震波形数据一最后一个GPS授时时间戳之后第k个采样点的采样时刻,直至补充至矿震波形数据一的最后1个采样点,停止补充; [0028] 采用算式 等间距补齐井下有线监测单元所监测的矿震波形数据二首个GPS授时时间戳之前第k个采样点的采样时刻,直至补充至矿震波形数据二的第1个采样点,停止补充;采用算式 等间距补齐井下有线监测单元所监测的矿震波形数据二最后一个GPS授时时间戳之后第k个采样点的采样时刻,直至补充至矿震波形数据二的最后1个采样点,停止补充。 [0029] 进一步地,所述步骤(4)中,所述的采样频率统一设置为f,补齐后 对应的采样时刻为 补齐后 对应的采样时刻为 对采样时刻进行重新取样的步骤为: [0030] 401)分别找到采样时刻 和 的最小值 和最大值 [0031] 402)确定重新取样后采样时刻 的起始值 结束值采样时刻 的起始值 结束值 式中mod为取模运算; [0032] 403)确定重新取样后的 为 [0033] [0034] 其中, 重新取样后的 为 [0035] [0036] 其中, [0037] 进一步地,所述步骤(5)中,根据采样时刻 和 计算对应矿震波形数据的步骤为: [0038] 501)找到满足 的下标位置index; [0039] 其中,k=1,2,...,z,index=1,2,...,z‑1,z=l或p,x=u或s; [0040] 502)利用下标位置index,采用线性插值公式 [0041] [0042] 计算得到波形数据 [0043] 本发明利用井地一体微震监测系统获得地面及井下监测单元采集的带有时间戳的矿震波形数据段,分别计算两段矿震波形数据中的所有相邻GPS授时时间戳的差值和相邻GPS授时时间戳内包含的采样点数,然后等间距补齐所有相邻GPS授时时间戳内未打时间戳的采样点的采样时刻,通过统计计算两段矿震波形数据的平均采样频率,分别等间距补齐两段采样数据中首、尾GPS授时时间戳之外的采样点的采样时刻;在获得所有采样点上的采样时刻后,按照统一固定采样频率对采样时刻进行重新取样,之后采用线性插值公式计算新采样时刻下的波形数据,进而对齐两段波形数据的采样时刻。本发明解决了井地一体微震监测系统中地面监测单元和井下监测单元采集的GPS授时波形数据之间无法对齐和同步的问题,有效提高了后续矿震震源参数求解精度,保证了矿震定位和能量求解的准确性,提升了微震监测精度。另外,本发明适用范围广,且易于电脑编程。附图说明 [0044] 图1是本发明的工作流程图; [0045] 图2是实施例中地面监测单元监测的矿震波形数据和GPS时间戳信息; [0046] 图3是实施例中井下监测单元监测的矿震波形数据和GPS时间戳信息; [0047] 图4是实施例中对地面监测单元补齐后的所有采样点采样时刻; [0048] 图5是实施例中对井下监测单元补齐后的所有采样点采样时刻; [0049] 图6是实施例中井下和地面监测单元监测波形授时对齐后,以时间为横轴、幅值为纵轴的波形图。 具体实施方式[0050] 下面将结合附图对本发明作进一步说明。 [0051] 如图1所示,一种井地一体微震监测系统监测单元采集数据授时对齐方法,包括如下步骤: [0052] (1)分别取由井地一体微震监测系统地面无线监测单元和井下有线监测单元独立采集的带有多个GPS授时时间戳的矿震波形数据一和矿震波形数据二;矿震波形数据一和矿震波形数据二为长度不同的两段数据段; [0053] (2)分别计算矿震波形数据一和矿震波形数据二中所有相邻GPS授时时间戳之间的差值以及相邻GPS授时时间戳内包含的采样点数后,等间距补齐相邻GPS授时时间戳内未打时间戳的采样点的采样时刻; [0054] (3)分别计算矿震波形数据一和矿震波形数据二的平均采样频率,分别等间距补齐矿震波形数据一和矿震波形数据二中的首、尾GPS授时时间戳之外的采样点的采样时刻; [0055] (4)根据步骤(2)和步骤(3)得到的所有采样点上的采样时刻,按照统一设置的采样频率对采样时刻进行重新取样; [0056] (5)根据步骤(2)和步骤(3)得到的所有采样点上的采样时刻和波形数据,采用线性插值公式计算步骤(4)中新采样时刻下的矿震波形数据,对齐矿震波形数据一和矿震波形数据二的采样时刻。 [0057] 进一步地,所述步骤(1)中,地面无线监测单元和井下有线监测单元分别采用独立的高精度GPS进行授时,其采样频率分别为fs和fu;所取的矿震波形数据一和矿震波形数据二分别为 和 其中,i=1,2,...,l,l为地面监测单元的采样长度;j=1,2,...,p,p为井下监测单元的采样长度;所述GPS授时时间戳为两个独立的GPS分别对矿震波形数据一和矿震波形数据二上的部分采样点打时间戳,其中, 的部分采样点 上的时间戳形s成GPS授时序列为 其中,i′=1,2,...,m;index为在 中有时间戳的位置序号,m为地面无线监测单元所监测的矿震波形数据一中带有GPS时间戳的采样点个数; 的部分采样u 点 上的时间戳形成的GPS授时序列为 其中,j′=1,2,...,n;index为在 中有时间戳的位置序号,n为井下有线监测单元所监测的矿震波形数据二中带有GPS时间戳的采样点个数。 [0058] 进一步地,所述步骤(2)中,所述等间距补齐相邻GPS授时时间戳内未打时间戳的采样点的采样时刻的步骤为: [0059] 201)分别计算GPS授时序列 中相邻时间戳的差值和其内包含的采样点数分别为: [0060] 时间差: [0061] 采样点数: [0062] 202)利用时间差和采样点数,采用公式: [0063] [0064] 等间距计算 上采样点 之后第k个点的采样时刻,直至计算至采样点停止计算; [0065] 采用公式: [0066] [0067] 等间距计算 上采样点 之后第k个点的采样时刻,直至计算至采样点停止计算; [0068] 203)重复步骤201)和202),直至补齐矿震波形数据一和矿震波形数据二中所有相邻GPS授时时间戳内未打时间戳的采样点的采样时刻。 [0069] 进一步地,所述步骤(3)中,等间距补齐矿震波形数据一和矿震波形数据二中的首、尾GPS授时时间戳之外的采样点的采样时刻的步骤为: [0070] 301)分别计算矿震波形数据一和矿震波形数据二的平均采样频率 和 为: [0071] [0072] [0073] 302)采用算式 等间距补齐地面无线监测单元所监测的矿震波形数据一首个GPS授时时间戳之前第k个采样点的采样时刻,直至补充至矿震波形数据一的第1个采样点,停止补充;采用算式 等间距补齐地面无线监测单元所监测的矿震波形数据一最后一个GPS授时时间戳之后第k个采样点的采样时刻,直至补充至矿震波形数据一的最后1个采样点,停止补充; [0074] 采用算式 等间距补齐井下有线监测单元所监测的矿震波形数据二首个GPS授时时间戳之前第k个采样点的采样时刻,直至补充至矿震波形数据二的第1个采样点,停止补充;采用算式 等间距补齐井下有线监测单元所监测的矿震波形数据二最后一个GPS授时时间戳之后第k个采样点的采样时刻,直至补充至矿震波形数据二的最后1个采样点,停止补充。 [0075] 进一步地,所述步骤(4)中,所述的采样频率统一设置为f,补齐后 对应的采样时刻为 补齐后 对应的采样时刻为 对采样时刻进行重新取样的步骤为: [0076] 401)分别找到采样时刻 和 的最小值 和最大值 [0077] 402)确定重新取样后采样时刻 的起始值 结束值采样时刻 的起始值 结束值 式中mod为取模运算; [0078] 403)确定重新取样后的 为 [0079] [0080] 其中, [0081] 重新取样后的 为 [0082] [0083] 其中, [0084] 进一步地,所述步骤(5)中,根据采样时刻 和 计算对应矿震波形数据的步骤为: [0085] 501)找到满足 的下标位置index; [0086] 其中,k=1,2,...,z,index=1,2,...,z‑1,z=l或p,x=u或s; [0087] 502)利用下标位置index,采用线性插值公式 [0088] [0089] 计算得到波形数据 [0090] 实施例: [0091] (1)如图2所示,取井地一体微震监测系统地面无线监测单元以采样频率fs=500Hz采集的带有GPS授时时间戳的矿震波形数据 其中, 的某些 采样点 上的时间戳形成的序列为 [0092] 如图3所示,取井地一体微震监测系统井下有线监测单元以采样频率fu=500Hz采集的带有GPS授时时间戳的矿震波形数据 的某些采样点上的时间戳形成的序列为 [0093] (2)等间距补齐相邻GPS授时时间戳内未打时间戳的采样点的采样时刻的步骤为: [0094] 201)计算波形数据 的GPS授时 相邻时间戳的差值 其内包含的采样点数为 计算波形数 据 的GPS授时 相邻时间戳的差值 其内包含的采样点数 为 [0095] 202)利用时间差和采样点数,采用式 [0096] [0097] 等间距计算 上采样点 之后第k个点的采样时刻,直至计算至采样点[0098] 采用式 [0099] [0100] 等间距计算 上采样点 之后第k个点的采样时刻,直至计算至采样点 [0101] 203)重复步骤201)和202),直至补齐两段波形数据中所有相邻GPS授时时间戳内未打时间戳的采样点的采样时刻。 [0102] (3)等间距补齐波形数据 和 中首、尾GPS授时时间戳之外的采样点的采样时刻的步骤为: [0103] 301)分别计算波形数据 和 的平均采样频率 和 为: [0104] [0105] [0106] 302)采用式 等间距补齐地面无线监测波形数据 首个GPS授时时间戳之前第k个采样点的采样时刻,直至波形数据的第1个采样点;采用式等间距补齐地面无线监测波形数据 最后一个GPS授时时间 戳之后第k个采样点的采样时刻,直至波形数据的最后1个采样点;采用式等间距补齐井下有线监测波形数据 首个GPS授时 时间戳之前第k个采样点的采样时刻,直至波形数据的第1个采样点;采用式等间距补齐井下有线监测波形数据 最后一个GPS 授时时间戳之后第k个采样点的采样时刻,直至波形数据的最后1个采样点。 [0107] (4)如图4和图5所示,利用步骤(2)和步骤(3)得到补齐后所有采样点上的采样时刻 和 按照统一固定采样频率f=500Hz对采样时刻进行重新取样的步骤为: [0108] 401)找到采样时刻 的最小值 和最大值 找到采样时刻 的最小值 和最大值 [0109] 402)确定重新取样后采样时刻 的起始值为: [0110] [0111] 结束值为: [0112] [0113] 采样时刻 的起始值为: [0114] [0115] 结束值为: [0116] [0117] 403)确定重新取样后的 为: [0118] [0119] [0120] 重新取样后的 为: [0121] [0122] [0123] (5)采用线性插值公式计算步骤(4)中新采样时刻 和计算对应波形数据进而对齐采样时刻的步骤为: [0124] 步骤501、找到满足 的下标位置index;k=1,2,...,l或p,l=6000,p=6505,index=1,2,...,5999,x=u或s; [0125] 步骤502、利用下标位置index,采用线性插值公式: [0126] [0127] 计算得到如图6所示的波形数据 或p,l=6000,p=6505,x=u或s。 |