CN107870362A - 一种基于样条插值的爆破振动响应预测方法 - Google Patents
一种基于样条插值的爆破振动响应预测方法 Download PDFInfo
- Publication number
- CN107870362A CN107870362A CN201711074727.4A CN201711074727A CN107870362A CN 107870362 A CN107870362 A CN 107870362A CN 201711074727 A CN201711074727 A CN 201711074727A CN 107870362 A CN107870362 A CN 107870362A
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- domain amplitude
- spectrum
- frequency
- waveform
- 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
- 238000005422 blasting Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000001228 spectrum Methods 0.000 claims abstract description 102
- 239000002360 explosive Substances 0.000 claims abstract description 12
- 230000000694 effects Effects 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 239000011435 rock Substances 0.000 claims abstract description 10
- 239000002689 soil Substances 0.000 claims abstract description 9
- 238000000605 extraction Methods 0.000 claims abstract description 6
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 238000004880 explosion Methods 0.000 claims description 14
- 238000012544 monitoring process Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000003814 drug Substances 0.000 claims description 2
- 239000002245 particle Substances 0.000 description 10
- 210000005252 bulbus oculi Anatomy 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000001845 vibrational spectrum Methods 0.000 description 1
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
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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)
Abstract
本发明公开了一种基于样条插值的爆破振动响应预测方法，包括：采集爆破荷载作用下岩土介质振动的实测数据；计算各测点的比例药量，并对采集的各测点波形进行波形滤波；对实测时域波形进行快速傅里叶变换，得到波形频域幅值谱和相位谱；提取波形频域幅值谱和相位谱在各时刻下的频域幅值和相位值；对每个频率成分在各个测点位置波形对应的频域幅值和相位值进行三次样条插值，得到拟合曲线；将预测点的比例药量代入拟合曲线中，得出各频率下的频域幅值和相位值，线性插值得到预测频域幅值谱和相位谱；得到时域振动预测波形。采用本发明能够预测一定范围内任意位置的完整爆破振动时间历程，并且预测所需准备的资料易采集，便于工程实际应用。
Description
技术领域
本发明属于爆破振动频谱预测领域，尤其涉及一种基于样条插值的爆破振动响应预测方法。
背景技术
在岩土爆破过程中，只有20％—30％的能量作用于岩体破坏上，其余的能量以振动、灰尘、飞石、冲击波、噪声等形式消散，其中，爆破振动被认为是爆破施工中最为主要的负面效应。爆破振动是国内外工程爆破领域学者的研究重点之一，过大的爆破振动会引起建筑物变形、山体滑坡、洞室坍塌等灾害事故。爆前评估和预测爆破振动强度，控制爆破振动处于安全范围内是保证工程爆破成功实施的重要前提。
爆破地震波在介质中的传播是一个复杂的力学过程，药量、爆心距、地质地形条件、岩土体介质动力学特性等因素都会影响爆破振动效应，现有的预测方法已逐渐难以满足日益提高的爆破振动安全要求。如经验公式法，虽然便于工程技术人员掌握，能够快速预测目标点的质点峰值速度，但是由于所得预测结果多数过于保守，经济性较差，并且无法预测介质质点的爆破振动全历程；动力分析法虽然理论性强，但是缺乏一定的经济性和简便性，难于在工程应用中大量推广。
现今质点峰值速度和主频已经成为实际工程中最为主要的爆破地震的安全评定依据，我国GB6722-2014安全评定标准中根据建筑物的结构形式和构筑质量划分了不同的安全质点峰值速度和频率范围，综合体现了爆破场地和建筑物基底场地的特性以及建筑物类型和其重要程度。因此在评估爆破振动效应时，常常同时研究质点峰值速度和频率这两个因素，为建筑物的安全防护提供理论依据。
发明内容
为了解决上述问题，本发明提出一种基于样条插值的爆破振动响应预测方法，对一定范围内任意位置的完整爆破振动时间历程进行预测，并且预测所需准备的资料易采集，能够快速获得预测结果，便于应用于工程实际。
一种基于样条插值的爆破振动响应预测方法，具体步骤如下：
(1)采集爆破荷载作用下岩土介质振动的实测数据；
(2)计算各测点的比例药量，并对步骤(1)中采集的各测点波形进行波形滤波；
(3)对实测时域波形进行快速傅里叶变换，得到波形频域幅值谱和相位谱；
(4)提取步骤(3)中得到的波形频域幅值谱和相位谱在各时刻下的频域幅值和相位值，并根据比例药量大小排列组合；
(5)对每个频率成分在各个测点位置波形对应的频域幅值和相位值进行三次样条插值，得到拟合曲线；
(6)将预测点的比例药量代入步骤(5)的拟合曲线中，得出各频率下的频域幅值和相位值，按照频率大小进行排列组合，线性插值得到预测频域幅值谱和相位谱；
(7)对步骤(6)得到的预测频域幅值谱和相位谱进行快速傅里叶逆变换，得到时域振动预测波形。
优选地，所述比例药量的计算公式如下：ρ＝Q1/3/R，其中Q为爆破总药量，R为爆心距，即爆炸源到测点的距离。
优选地，所述采集爆破荷载作用下岩土介质振动的实测数据是指对n个监测点进行采样率为fs，采样时长为T的爆破振动信号采样，采样点总量N＝fsT，对每一个已知的监测点计算其比例药量：ρi＝Qi 1/3/Ri(i＝0,1,2...,n)，并计算预测点的比例药量：
优选地，步骤(3)中通过对实测时域波形快速傅里叶变换，得到间隔频率为dHz,频率范围为0～mdHz的频域幅值谱和相位谱，共有n组不同测点的频域幅值谱和相位谱，每组数据共有m个数据点，采样数据所对应的频率为fj＝jd(j＝0,1,2...,m)。
优选地，提取步骤(3)中得到的频域幅值谱和相位谱在各时刻下的频域幅值和相位值并根据比例药量大小排列组合，此处和中右上标i表示第i(i＝1,2,.....,n)个测点，右下标j表示fj(j＝1,2,.....m)频率。
优选地，定义三次样条函数Sr(x)中自变量x为比例药量，即x＝Q1/3/R，因变量为频域幅值A和相位值θ，即三次样条函数Sr(x)＝Sr(Q1/3/R)＝A和Sr(x)＝Sr(Q1/3/R)＝θ，选用自然边界作为约束条件，即S″r(x0)＝S″r(xn-1)＝0，通过三次样条函数拟合可得到频率fj下的拟合函数其中Q为爆破总药量，R为爆心距，A为频域幅值，θ为相位值。
优选地，拟合函数中令fj＝fpredict，求得各频率下的频域幅值和相位，按照频率大小得到离散的预测频域幅值谱和相位谱，再通过线性插值得到连续的预测频域幅值谱和相位谱。
优选地，将预测的频域幅值谱和相位谱进行快速傅里叶逆变换，即得到预测的时域波形图。
有益效果：本发明公开了一种基于样条插值的爆破振动响应预测方法，能够预测一定范围内任意位置的完整爆破振动时间历程，并且预测所需准备的资料易采集，便于工程实际应用。
附图说明
图1为基于样条插值的爆破振动响应预测方法的技术路线图；
图2为监测数据时域波形图；
图3为各比例药量下频域幅值谱；
图4为各比例药量下相位谱；
图5为频域幅值谱与比例药量关系图；
图6为相位谱与比例药量关系图；
图7为三次样条完成的频域幅值谱与比例药量关系图；
图8为三次样条完成的相位谱与比例药量关系图；
图9为离散化的预测频域幅值谱；
图10为离散化的预测相位谱；
图11为连续的预测频域幅值谱；
图12为连续的预测相位谱；
图13为预测时域波形图；
图14为实测点信息表；
图15为实测点波形图；
图16为实测点频域幅值谱图；
图17为实测点相位谱图；
图18为各频率下的频域幅值与比例药量的三次样条插值函数拟合曲线图；
图19为各频率下的相位与比例药量的三次样条插值函数拟合曲线图；
图20为实测频域幅值谱与预测频域幅值谱对比图；
图21为实测相位谱与预测相位谱对比图；
图22为实测时域波形与预测时域波形对比图；
图23为实测波形与预测波形频域部分对比表；
图24为实测质点峰值振动速度预测结果对比表。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
一种基于样条插值的爆破振动响应预测方法，具体步骤如下：
(1)采集爆破荷载作用下岩土介质振动的实测数据；
(2)计算各测点的比例药量，并对步骤(1)中采集的各测点波形进行波形滤波；
(3)对实测时域波形进行快速傅里叶变换，得到波形频域幅值谱和相位谱；
(4)提取步骤(3)中得到的波形频域幅值谱和相位谱在各时刻下的频域幅值和相位值，并根据比例药量大小排列组合；
(5)对每个频率成分在各个测点位置波形对应的频域幅值和相位值进行三次样条插值，得到拟合曲线；
(6)将预测点的比例药量代入步骤(5)的拟合曲线中，得出各频率下的频域幅值和相位值，按照频率大小进行排列组合，线性插值得到预测频域幅值谱和相位谱；
(7)对步骤(6)得到的预测频域幅值谱和相位谱进行快速傅里叶逆变换，得到时域振动预测波形。
优选地，所述比例药量的计算公式如下：ρ＝Q1/3/R，其中Q为爆破总药量，R为爆心距，即爆炸源到测点的距离。
优选地，所述采集爆破荷载作用下岩土介质振动的实测数据是指对n个监测点进行采样率为fs，采样时长为T的爆破振动信号采样，采样点总量N＝fsT，对每一个已知的监测点计算其比例药量：ρi＝Qi 1/3/Ri(i＝0,1,2...,n)，并计算预测点的比例药量：
优选地，步骤(3)中通过对实测时域波形快速傅里叶变换，得到间隔频率为dHz,频率范围为0～mdHz的频域幅值谱和相位谱，共有n组不同测点的频域幅值谱和相位谱，每组数据共有m个数据点，采样数据所对应的频率为fj＝jd(j＝0,1,2...,m)。
优选地，提取步骤(3)中得到的频域幅值谱和相位谱在各时刻下的频域幅值和相位值并根据比例药量大小排列组合，此处和中右上标i表示第i(i＝1,2,.....,n)个测点，右下标j表示fj(j＝1,2,.....m)频率。
优选地，定义三次样条函数Sr(x)中自变量x为比例药量，即x＝Q1/3/R，因变量为频域幅值A和相位值θ，即三次样条函数Sr(x)＝Sr(Q1/3/R)＝A和Sr(x)＝Sr(Q1/3/R)＝θ，选用自然边界作为约束条件，即S″r(x0)＝S″r(xn-1)＝0，通过三次样条函数拟合可得到频率fj下的拟合函数其中Q为爆破总药量，R为爆心距，A为频域幅值，θ为相位值。
优选地，拟合函数中令fj＝fpredict，求得各频率下的频域幅值和相位，按照频率大小得到离散的预测频域幅值谱和相位谱，再通过线性插值得到连续的预测频域幅值谱和相位谱。
优选地，将预测的频域幅值谱和相位谱进行快速傅里叶逆变换，即得到预测的时域波形图。
如图1所示，为一种基于样条插值的爆破振动响应预测方法的技术路线图，包括：采集爆破荷载作用下岩土介质振动的实测数据并计算各测点的比例药量，对各测点波形进行滤波；对实测时域波形进行快速傅里叶变换，得到频域幅值谱和相位谱；根据各测点的比例药量，提取每个频率成分在各个测点位置波形对应的频谱值，对每个频率成分在各个测点位置波形对应的频谱进行三次样条插值，得到拟合曲线；计算预测点的比例药量并带入拟合曲线，提取相应的各频率下的幅值和相位，按照频率大小进行排列组合，得到预测频域幅值谱和相位谱；对预测频域幅值和相位谱进行傅里叶逆变换，得到时域振动预测波形。
如图2所示，为监测数据时域波图，对n个监测点进行采样率为fs，采样时长为T的爆破振动信号采样，采样点总量N＝fsT。对每一个已知的监测点计算其ρi＝Qi 1/3/Ri(i＝0,1,2...,n)，并计算预测点的
如图3和图4所示，分别为各比例药量下的频域幅值谱和相位谱，通过快速傅里叶变换，得到间隔频率为dHz,频率范围为0～mdHz的频域幅值谱和相位谱，共有n组不同测点的频域幅值谱和相位谱，每组数据共有m个数据点，采样数据所对应的频率为fi＝jd(j＝0,1,2...,m)j＝0,1,2...,m。
如图5和图6所示，分别为频域幅值谱与比例药量关系图和相位谱与比例药量关系图，提取相应的频域幅值谱和相位谱在各时刻下的频域幅值和相位值并根据比例药量大小排列组合。
如图7和图8所示，分别为三次样条插值完成的频域幅值谱与比例药量关系图和三次样条完成的相位谱与比例药量关系图，其中定义三次样条函数Sr(x)中自变量x＝Q1/3/R，因变量为频域幅值A和相位值θ，即三次样条函数Sr(x)＝Sr(Q1/3/R)＝A和Sr(x)＝Sr(Q1/3/R)＝θ，选用自然边界作为约束条件，即S″r(x0)＝S″r(xn-1)＝0。其中Q为爆破总药量，R为爆心距，A为频域幅值，θ为相位值。
如图9和图10所示，分别为离散化的预测频域幅值谱和相位谱，如图11和图12所示，分别为连续的预测频域幅值谱和相位谱，拟合函数中令fj＝fpredict，可求得各频率下的频域幅值和相位，通过排列组合便可得到离散的预测频域幅值谱和相位谱。再通过线性插值便可得到连续的预测频域幅值谱和相位谱。
如图13所示，为预测时域波形图，预测的频域幅值谱和相位谱进行快速傅里叶逆变换即得到预测的时域波形图。
如图14所示，为河北衡水某大厦爆破测点的实测点信息表，Ve为质点峰值速度，5个测点共5组监测数据，选取测点四(P4)作为预测点，其余4组测点数据作为数据样本。
如图15所示，为5个测点实测的竖向振动波形。
如图16和图17所示，为通过对各个测点实测的竖向振动波形快速傅里叶变换获得的各测点频域幅值谱和相位谱图。
如图18和图19所示，频率间隔设置为20Hz，构建各频率下的频谱值、相位与比例药量拟合曲线图。
如图20和图21所示，根据预测点的比例药量，提取各频率下相应的频域幅值和相位，得到预测频域幅值谱和相位谱。
如图22所示，为实测时域波形与预测时域波形对比图。
如图23所示，为实测波形与预测波形频域部分对比表，从表中质点峰值速度、炸药量与爆心距关系式拟合结果为Ve＝Q0.99155R-0.44569。
如图24所示，为基于三次样条插值法的预测质点峰值振动速度和经验公式法结果对比表，预测质点峰值速度的相对误差为20.5％，低于经验公式法的33.4％。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的两种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种基于样条插值的爆破振动响应预测方法，其特征在于，具体步骤如下：
(1)采集爆破荷载作用下岩土介质振动的实测数据；
(2)计算各测点的比例药量，并对步骤(1)中采集的各测点波形进行波形滤波；
(3)对实测时域波形进行快速傅里叶变换，得到波形频域幅值谱和相位谱；
(4)提取步骤(3)中得到的波形频域幅值谱和相位谱在各时刻下的频域幅值和相位值，并根据比例药量大小排列组合；
(5)对每个频率成分在各个测点位置波形对应的频域幅值和相位值进行三次样条插值，得到拟合曲线；
(6)将预测点的比例药量代入步骤(5)的拟合曲线中，得出各频率下的频域幅值和相位值，按照频率大小进行排列组合，线性插值得到预测频域幅值谱和相位谱；
(7)对步骤(6)得到的预测频域幅值谱和相位谱进行快速傅里叶逆变换，得到时域振动预测波形。
2.根据权利要求书1所述的一种基于样条插值的爆破振动响应预测方法，其特征在于：所述比例药量的计算公式如下：ρ＝Q1/3/R，其中Q为爆破总药量，R为爆心距，即爆炸源到测点的距离。
3.根据权利要求书2所述的一种基于样条插值的爆破振动响应预测方法，其特征在于：所述采集爆破荷载作用下岩土介质振动的实测数据是指对n个监测点进行采样率为fs，采样时长为T的爆破振动信号采样，采样点总量N＝fsT，对每一个已知的监测点计算其比例药量：ρi＝Qi 1/3/Ri(i＝0,1,2...,n)，并计算预测点的比例药量：
4.根据权利要求书3所述的一种基于样条插值的爆破振动响应预测方法，其特征在于：步骤(3)中通过对实测时域波形快速傅里叶变换，得到间隔频率为dHz,频率范围为0～mdHz的频域幅值谱和相位谱，共有n组不同测点的频域幅值谱和相位谱，每组数据共有m个数据点，采样数据所对应的频率为fj＝jd(j＝0,1,2...,m)。
5.根据权利要求书4述的一种基于样条插值的爆破振动响应预测方法，其特征在于：提取步骤(3)中得到的频域幅值谱和相位谱在各时刻下的频域幅值和相位值并根据比例药量大小排列组合，此处和中右上标i表示第i(i＝1,2,.....,n)个测点，右下标j表示fj(j＝1,2,.....m)频率。
6.根据权利要求书5所述的一种基于样条插值的爆破振动响应预测方法，其特征在于：定义三次样条函数Sr(x)中自变量x为比例药量，即x＝Q1/3/R，因变量为频域幅值A和相位值θ，即三次样条函数Sr(x)＝Sr(Q1/3/R)＝A和Sr(x)＝Sr(Q1/3/R)＝θ，选用自然边界作为约束条件，即S″r(x0)＝S″r(xn-1)＝0，通过三次样条函数拟合可得到频率fj下的拟合函数其中Q为爆破总药量，R为爆心距，A为频域幅值，θ为相位值。
7.根据权利要求书6所述的一种基于样条插值的爆破振动响应预测方法，其特征在于：拟合函数中令fj＝fpredict，求得各频率下的频域幅值和相位，按照频率大小得到离散的预测频域幅值谱和相位谱，再通过线性插值得到连续的预测频域幅值谱和相位谱。
8.根据权利要求书7所述的一种基于样条插值的爆破振动响应预测方法，其特征在于：将预测的频域幅值谱和相位谱进行快速傅里叶逆变换，即得到预测的时域波形图。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711074727.4A CN107870362B (zh) | 2017-11-06 | 2017-11-06 | 一种基于样条插值的爆破振动响应预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711074727.4A CN107870362B (zh) | 2017-11-06 | 2017-11-06 | 一种基于样条插值的爆破振动响应预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107870362A true CN107870362A (zh) | 2018-04-03 |
| CN107870362B CN107870362B (zh) | 2019-08-06 |
Family
ID=61753433
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711074727.4A Expired - Fee Related CN107870362B (zh) | 2017-11-06 | 2017-11-06 | 一种基于样条插值的爆破振动响应预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107870362B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109752085A (zh) * | 2018-11-09 | 2019-05-14 | 北方爆破科技有限公司 | 一种数码电子雷管逐孔起爆振动预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002055090A (ja) * | 2000-08-09 | 2002-02-20 | Ohbayashi Corp | 液状化判定方法 |
| CN102135445A (zh) * | 2010-06-30 | 2011-07-27 | 中国铁道科学研究院铁道建筑研究所 | 爆破振动预测方法 |
| CN102721978A (zh) * | 2012-06-19 | 2012-10-10 | 中国地震灾害防御中心 | 获得地震动样本的方法和装置 |
| CN104281785A (zh) * | 2014-10-31 | 2015-01-14 | 中国矿业大学(北京) | 一种基于置信水平预测爆破振动速度的方法及装置 |
| CN106126904A (zh) * | 2016-06-21 | 2016-11-16 | 四川大学 | 多层建筑物的吸收爆破振动能量舒适性评价方法 |
-
2017
- 2017-11-06 CN CN201711074727.4A patent/CN107870362B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002055090A (ja) * | 2000-08-09 | 2002-02-20 | Ohbayashi Corp | 液状化判定方法 |
| CN102135445A (zh) * | 2010-06-30 | 2011-07-27 | 中国铁道科学研究院铁道建筑研究所 | 爆破振动预测方法 |
| CN102721978A (zh) * | 2012-06-19 | 2012-10-10 | 中国地震灾害防御中心 | 获得地震动样本的方法和装置 |
| CN104281785A (zh) * | 2014-10-31 | 2015-01-14 | 中国矿业大学(北京) | 一种基于置信水平预测爆破振动速度的方法及装置 |
| CN106126904A (zh) * | 2016-06-21 | 2016-11-16 | 四川大学 | 多层建筑物的吸收爆破振动能量舒适性评价方法 |
Non-Patent Citations (2)
| Title |
|---|
| 杨仁树等: "基于HHT方法的爆炸应变波时频分析", 《振动与冲击》 * |
| 楼矗云等: "基于样条插值的爆破振动时域波形预测方法", 《工程爆破》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109752085A (zh) * | 2018-11-09 | 2019-05-14 | 北方爆破科技有限公司 | 一种数码电子雷管逐孔起爆振动预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107870362B (zh) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xiao et al. | ISRM suggested method for in situ microseismic monitoring of the fracturing process in rock masses | |
| CN101770038B (zh) | 矿山微震源智能定位方法 | |
| CN103306722B (zh) | 一种冲击危险区的微震多维信息综合区域探测评价方法 | |
| CN102253410B (zh) | 全数字化工程双源面波勘探数据采集与处理系统 | |
| CN106840381A (zh) | 一种爆破振动衰减效应的测试方法 | |
| CN107315862A (zh) | 一种建立明挖基坑工程勘察与模拟参数关系的方法 | |
| CN103336297B (zh) | 微破裂向量扫描方法 | |
| CN105956526A (zh) | 基于多尺度排列熵的低信噪比微震事件辨识方法 | |
| CN108508482A (zh) | 一种地下裂缝地震散射响应特征模拟方法 | |
| CN107727737A (zh) | 一种开挖卸荷作用下地下工程未知地质构造活化识别方法 | |
| CN108733942A (zh) | 隧道围岩可靠性的评估方法及装置 | |
| CN105676268A (zh) | 一种基于声音信号波形变化特征的应变型岩爆预警方法 | |
| CN106199700A (zh) | 一种地下水封储油洞库微震监测方法及系统 | |
| CN103323530A (zh) | 利用爆破作业震动波检测巷道群围岩稳定性的方法 | |
| Ma et al. | Intelligent location of microseismic events based on a fully convolutional neural network (FCNN) | |
| CN104297788A (zh) | 基于波形起振趋势线斜率的矿山微震和爆破信号识别方法 | |
| CN109190136A (zh) | 面向地表沉陷动态预计的数值模型岩体力学参数加权反演方法 | |
| CN102073061A (zh) | 使用数字地听仪的地听信息高密度记录系统 | |
| CN114943149A (zh) | 一种隧道内岩爆损伤岩体体积的计算方法 | |
| CN107870362A (zh) | 一种基于样条插值的爆破振动响应预测方法 | |
| Xu et al. | Optimal design of microseismic monitoring networking and error analysis of seismic source location for rock slope | |
| Chiaraluce et al. | The shallow boreholes at The AltotiBerina near fault Observatory (TABOO; northern Apennines of Italy) | |
| Uski et al. | Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region, Fennoscandian Shield | |
| CN205861210U (zh) | 爆破振动全天候远程智能化监测系统 | |
| Dal Moro et al. | Gaussian-filtered Horizontal Motion (GHM) plots of non-synchronous ambient microtremors for the identification of flexural and torsional modes of a building |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Liu Mingqing Inventor after: Xie Deqiang Inventor after: Kong Qingjie Inventor after: Zhang Zhixiang Inventor after: Zhao Xueyi Inventor after: Wu Qian Inventor after: Liu Tingting Inventor before: Xie Deqiang Inventor before: Kong Qingjie Inventor before: Zhang Zhixiang Inventor before: Zhao Xueyi Inventor before: Wu Qian Inventor before: Liu Tingting |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190806 |