CN109597124A - 一种工作面应力集中区槽波地震精细探测方法 - Google Patents
一种工作面应力集中区槽波地震精细探测方法 Download PDFInfo
- Publication number
- CN109597124A CN109597124A CN201811633855.2A CN201811633855A CN109597124A CN 109597124 A CN109597124 A CN 109597124A CN 201811633855 A CN201811633855 A CN 201811633855A CN 109597124 A CN109597124 A CN 109597124A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity
- working face
- stress
- rock
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 239000003245 coal Substances 0.000 claims abstract description 35
- 239000011435 rock Substances 0.000 claims abstract description 35
- 239000006185 dispersion Substances 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims description 14
- 238000004587 chromatography analysis Methods 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims description 2
- 239000004568 cement Substances 0.000 claims 1
- 239000000843 powder Substances 0.000 abstract description 2
- 230000005856 abnormality Effects 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000004575 stone Substances 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/01—Measuring or predicting earthquakes
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
技术领域
本发明属于煤矿地质灾害预测技术领域，具体涉及一种工作面应力集中区槽波地震精细探测新方法。
背景技术
随着煤炭资源的大量开采，浅部资源逐渐枯竭，煤层开采深度逐年增大，面临的地应力也随之增大。对于地质构造复杂的矿井，由于受多期构造运动叠加影响，还会产生较大的构造应力。高应力对巷道的掘进与维护、硐室的开挖和工作面的回采都会造成很大的影响，甚至阻碍正常的生产。同时，矿井采掘过程中遇到这些高应力区也极易发生巷道变形、顶板垮落，造成压架事故，损坏设备，甚至人员伤亡，造成巨大的经济损失。
目前，针对工作面高应力区地质灾害，现有技术中还没有一种有效的预测预报方法。
发明内容
本发明的目的在于提供一种工作面应力集中区槽波地震精细探测方法，其能更好的为煤矿地质灾害进行预测预报。
其技术解决方案包括：
一种工作面应力集中区槽波地震精细探测方法，依次包括以下步骤：
a、通过煤岩体加载条件下纵波和横波波速测试，以获得应力—应变、波速—应变全过程变化曲线；绘制应力—波速变化全过程曲线，并分阶段进行线性和非线性公式，得拟合应力—波速拟合公式；
b、根据地质资料及井田地应力分布特征，计算工作面正常区域地应力大小；
c、基于步骤a所述的应力—波速拟合公式，根据步骤b所述的工作面正常区域地应力大小，预计工作面正常区煤岩体的纵波和横波波速；
d、根据工作面煤岩体物性参数，建立三层对称介质模型，通过数值计算分析槽波理论频散曲线特征；
e、收集现场采集的透射槽波数据并对其进行时频分析，绘制实测时频曲线；
f、将实测时频曲线与理论频散曲线进行对比分析，选取槽波频散曲线稳定段某一频率，拾取该频率对应的槽波旅行时；
g、利用观测***和槽波旅行时信息，并采用层析成像技术计算和绘制槽波速度云图，圈定高速异常区，结合地质条件确定应力集中区及范围。
作为本发明的一个优选方案，步骤a中，所述的煤岩体加载条件下纵波和横波波速测试，是借助岩石力学实验机和声波测试***，加载时采用位移控制法，轴向位移每变化0.01mm采集1次声波数据，直到岩石试件发生破坏或设定的极限变形为止，来获得应力—应变、波速—应变全过程变化曲线。
作为本发明的另一个优选方案，步骤a中，所述的绘制应力—波速变化全过程曲线，并分阶段进行线性和非线性公式拟合步骤，是根据应力、波速与应变的关系，借助应变中间变量，转化为应力与波速的关系，将岩石加载破坏全过程的五个阶段对应的波速响应特征进行应力与波速的线性和非线性公式拟合。
进一步的，步骤d中，根据工作面煤岩体物性参数，建立三层对称介质模型，通过数值计算分析槽波理论频散曲线特征，是指基于工作面煤岩体的纵波、横波波速及密度等参数测试结果，建立三层对称槽波形成的介质模型，模型中L型槽波群速度vL和频率f的关系式如式(1)、(2)所示：
式中，vx为槽波相速度，m/s；vsc为煤层横波波速，m/s；vsr为岩层层横波波速，m/s；d为1/2煤层厚度m；ρc为煤层密度，g/m3；ρr为岩层密度，g/m3；
通过数值计算，带入煤岩体物性参数，得出工作面内槽波理论频散曲线，确定频散曲线稳定段的槽波频率范围。
进一步的，步骤f中，所述的选取槽波频散曲线稳定段某一频率，拾取该频率对应的槽波旅行时，是通过对比实测时频曲线与理论频散曲线，并结合多炮多道记录的频散曲线特征，选取合适的频率，由于频率与旅行时为一一对应关系，进而可进行旅行时记录的提取。
进一步的，所述的采用层析成像技术计算和绘制槽波速度云图，圈定高速异常区，是根据观测***与槽波旅行时相关信息，通过层析成像方法反演计算出槽波速度，并绘制测区的速度云图，将云图中波速超过理论速度的范围视为高速区，结合工作面地质特征综合圈定高速异常区。
与现有技术相比，本发明带来了以下有益技术效果：
本发明提供的工作面应力集中区槽波地震精细探测方法，由于采用现场实测数据与理论成果进行对比分析，通过对槽波频散曲线稳定段进行旅行时提取，操作方便、简洁，人为误差小，能更精准的圈定异常范围，提高探测精度，更好的为煤矿地质灾害进行预测预报。
附图说明
下面结合附图对本发明做进一步说明：
图1是本发明一种工作面应力集中区槽波地震精细探测方法工艺流程图；
图2是具体实施方式中煤样全应力—应变过程波速变化曲线；
图3是具体实施方式中三层对称介质模型示意图；
图4是具体实施方式中不同煤层波速槽波理论频散曲线图；
图5是具体实施方式中实测槽波频散曲线图；
图6是具体实施方式中槽波速度层析成像成果图。
具体实施方式
本发明提出了一种工作面应力集中区槽波地震精细探测方法，为了使本发明的优点、技术方案更加清楚、明确，下面结合具体实施例对本发明做详细说明。
如图1所示，本发明工作面应力集中区槽波地震精细探测方法，具体包括以下步骤：
第一步、首先对煤岩体在加载条件下进行纵波和横波波速测试，借助岩石力学实验机和声波测试***，加载时采用位移控制法，轴向位移每变化0.01mm采集1次声波数据，直到岩石试件发生破坏或设定的极限变形为止，以获得应力—应变、波速—应变全过程变化曲线；绘制煤样全应力—波速变化全过程曲线，如图2所示，并分阶段进行线性和非线性公式，得拟合应力—波速拟合公式；
第二步、根据地质资料及井田地应力分布特征，具体的，根据地质资料、煤岩体岩石物理力学参数，计算工作面正常区域地应力大小；
第三步、基于应力—波速拟合公式，根据工作面地应力大小，预计工作面正常区煤岩体的纵波和横波波速；具体的，将计算出的工作面地应力结果带入基于应力—波速拟合公式，从而得出工作面正常区纵波和横波波速；
第四步、根据工作面煤岩体物性参数，建立三层对称介质模型，其示意图如图3所示，通过数值计算分析槽波理论频散曲线特征；具体的，是指基于工作面煤岩体的纵波、横波波速及密度等参数测试结果，建立三层对称槽波形成的介质模型。模型中L型槽波群速度vL和频率f的关系式：
式中，vx为槽波相速度，m/s；vsc为煤层横波波速，m/s；vsr为岩层层横波波速，m/s；d为1/2煤层厚度m；ρc为煤层密度，g/m3；ρr为岩层密度，g/m3；
通过数值计算，带入煤岩体物性参数，得出工作面内槽波理论频散曲线，确定频散曲线稳定变化段的槽波频率范围；
第五步、将现场采集的透射槽波数据进行时频分析，绘制实测槽波旅行时与频率关系的时频曲线；
第六步、将实测时频曲线与理论频散曲线进行对比分析，选取槽波频散曲线稳定段某一频率，拾取该频率对应的槽波旅行时；通过对比实测时频曲线与理论频散曲线，并结合多炮多道记录的频散曲线特征，选取200Hz的频率，并对该频率对应的旅行时记录进行提取；
第七步、利用观测***和槽波旅行时信息，并采用层析成像技术计算和绘制槽波速度云图，圈定高速异常区，结合地质条件确定应力集中区及范围。具体的，根据观测***与槽波旅行时信息，通过层析成像方法反演计算出槽波速度，并绘制测区的速度云图，将云图中波速超过1000m/s的范围视为高速区，结合工作面地质特征综合圈定高速异常区。
图4是上述不同煤层波速槽波理论频散曲线图，图5是具体实施方式中实测槽波频散曲线图，图6是具体实施方式中槽波速度层析成像成果图。结合上述说明及图1-图6可知，本发明能更精准的圈定异常范围，提高探测精度，更好的为煤矿地质灾害进行预测预报。
本发明未述及的部分借鉴现有技术即可实现。
需要说明的是：在本说明书的教导下本领域技术人员所做出的任何等同方式，或明显变型方式均应在本发明的保护范围内。
Claims (7)
1.一种工作面应力集中区槽波地震精细探测方法，其特征在于，依次包括以下步骤：
a、通过煤岩体加载条件下纵波和横波波速测试，以获得应力—应变、波速—应变全过程变化曲线；绘制应力—波速变化全过程曲线，并分阶段进行线性和非线性公式，得拟合应力—波速拟合公式；
b、根据地质资料及井田地应力分布特征，计算工作面正常区域地应力大小；
c、基于步骤a所述的应力—波速拟合公式，根据步骤b所述的工作面正常区域地应力大小，预计工作面正常区煤岩体的纵波和横波波速；
d、根据工作面煤岩体物性参数，建立三层对称介质模型，通过数值计算分析槽波理论频散曲线特征；
e、收集现场采集的透射槽波数据并对其进行时频分析，绘制实测时频曲线；
f、将实测时频曲线与理论频散曲线进行对比分析，选取槽波频散曲线稳定段某一频率，拾取该频率对应的槽波旅行时；
g、利用观测***和槽波旅行时信息，并采用层析成像技术计算和绘制槽波速度云图，圈定高速异常区，结合地质条件确定应力集中区及范围。
2.根据权利要求1所述的一种工作面应力集中区槽波地震精细探测方法，其特征在于：步骤a中，所述的煤岩体加载条件下纵波和横波波速测试，是借助岩石力学实验机和声波测试***，加载时采用位移控制法，轴向位移每变化0.01mm采集1次声波数据，直到岩石试件发生破坏或设定的极限变形为止，来获得应力—应变、波速—应变全过程变化曲线。
3.根据权利要求1所述的一种工作面应力集中区槽波地震精细探测方法，其特征在于：步骤a中，所述的绘制应力—波速变化全过程曲线，并分阶段进行线性和非线性公式拟合步骤，是根据应力、波速与应变的关系，借助应变中间变量，转化为应力与波速的关系，将岩石加载破坏全过程的五个阶段对应的波速响应特征进行应力与波速的线性和非线性公式拟合。
4.根据权利要求1所述的一种工作面应力集中区槽波地震精细探测方法，其特征在于：步骤d中，根据工作面煤岩体物性参数，建立三层对称介质模型，通过数值计算分析槽波理论频散曲线特征，是指基于工作面煤岩体的纵波、横波波速及密度等参数测试结果，建立三层对称槽波形成的介质模型，模型中L型槽波群速度vL和频率f的关系式如式(1)、(2)所示：
式中，vx为槽波相速度，m/s；vsc为煤层横波波速，m/s；vsr为岩层层横波波速，m/s；d为1/2煤层厚度m；ρc为煤层密度，g/m3；ρr为岩层密度，g/m3；
通过数值计算，带入煤岩体物性参数，得出工作面内槽波理论频散曲线，确定频散曲线稳定段的槽波频率范围。
5.根据权利要求1所述的一种工作面应力集中区槽波地震精细探测方法，其特征在于：步骤f中，所述的选取槽波频散曲线稳定段某一频率，拾取该频率对应的槽波旅行时，是通过对比实测时频曲线与理论频散曲线，并结合多炮多道记录的频散曲线特征，选取合适的频率，由于频率与旅行时为一一对应关系，进而可进行旅行时记录的提取。
6.根据权利要求1所述的一种工作面应力集中区槽波地震精细探测方法，其特征在于：所述的采用层析成像技术计算和绘制槽波速度云图，圈定高速异常区，是根据观测***与槽波旅行时相关信息，通过层析成像方法反演计算出槽波速度，并绘制测区的速度云图，将云图中波速超过理论速度的范围视为高速区，结合工作面地质特征综合圈定高速异常区。
7.根据权利要求6所述的一种工作面应力集中区槽波地震精细探测方法，其特征在于：所述的速度云图中波速超过1000m/s的范围视为高速区。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811633855.2A CN109597124B (zh) | 2018-12-29 | 2018-12-29 | 一种工作面应力集中区槽波地震精细探测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811633855.2A CN109597124B (zh) | 2018-12-29 | 2018-12-29 | 一种工作面应力集中区槽波地震精细探测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109597124A true CN109597124A (zh) | 2019-04-09 |
| CN109597124B CN109597124B (zh) | 2021-06-04 |
Family
ID=65964601
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811633855.2A Active CN109597124B (zh) | 2018-12-29 | 2018-12-29 | 一种工作面应力集中区槽波地震精细探测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109597124B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110118991A (zh) * | 2019-05-16 | 2019-08-13 | 中国矿业大学 | 一种基于微震损伤重构的采动应力评估方法 |
| CN111208567A (zh) * | 2020-01-07 | 2020-05-29 | 中国科学院地理科学与资源研究所 | 一种矿层成像方法、设备及计算机可读存储介质 |
| CN111476891A (zh) * | 2020-02-13 | 2020-07-31 | 山西晋煤集团技术研究院有限责任公司 | 一种井斜坐标数据批量处理的方法 |
| CN111551989A (zh) * | 2020-05-20 | 2020-08-18 | 中国科学院地理科学与资源研究所 | 一种透射槽波成像方法、设备及计算机可读存储介质 |
| WO2021008062A1 (zh) * | 2019-07-17 | 2021-01-21 | 山东科技大学 | 一种槽波信号多径谱分析方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009236874A (ja) * | 2008-03-28 | 2009-10-15 | General Environmental Technos Co Ltd | 二酸化炭素の地中浸透モニタリング方法 |
| CN106610503A (zh) * | 2016-11-29 | 2017-05-03 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿采掘过程中的全方位槽波地震探测方法 |
| CN106680868A (zh) * | 2017-01-05 | 2017-05-17 | 南京君诺鑫成网络科技有限公司 | 一种槽波地震多属性参数层析方法 |
| CN108490488A (zh) * | 2018-03-19 | 2018-09-04 | 中煤科工集团西安研究院有限公司 | 一种基于时窗能量比的槽波地震散射成像方法及装置 |
-
2018
- 2018-12-29 CN CN201811633855.2A patent/CN109597124B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009236874A (ja) * | 2008-03-28 | 2009-10-15 | General Environmental Technos Co Ltd | 二酸化炭素の地中浸透モニタリング方法 |
| CN106610503A (zh) * | 2016-11-29 | 2017-05-03 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿采掘过程中的全方位槽波地震探测方法 |
| CN106680868A (zh) * | 2017-01-05 | 2017-05-17 | 南京君诺鑫成网络科技有限公司 | 一种槽波地震多属性参数层析方法 |
| CN108490488A (zh) * | 2018-03-19 | 2018-09-04 | 中煤科工集团西安研究院有限公司 | 一种基于时窗能量比的槽波地震散射成像方法及装置 |
Non-Patent Citations (5)
| Title |
|---|
| M.GE ET AL.: "Void detection at an anthracite mine using an in-seam seismic method", 《INTERNATIONAL JOURNAL OF COAL GEOLOGY》 * |
| 乐勇等: "槽波地震勘探技术在工作面小构造探测中的应用", 《煤田地质与勘探》 * |
| 巩思园等: "冲击倾向煤岩纵波波速与应力关系试验研究", 《采矿与安全工程学报》 * |
| 杨真等: "0.9m薄煤层SH型槽波频散特征及波形模式", 《地球物理学报》 * |
| 王伟等: "槽波层析成像方法在煤田勘探中的应用-以河南义马矿区为例", 《地球物理学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110118991A (zh) * | 2019-05-16 | 2019-08-13 | 中国矿业大学 | 一种基于微震损伤重构的采动应力评估方法 |
| CN110118991B (zh) * | 2019-05-16 | 2020-06-23 | 中国矿业大学 | 一种基于微震损伤重构的采动应力评估方法 |
| WO2021008062A1 (zh) * | 2019-07-17 | 2021-01-21 | 山东科技大学 | 一种槽波信号多径谱分析方法及*** |
| CN111208567A (zh) * | 2020-01-07 | 2020-05-29 | 中国科学院地理科学与资源研究所 | 一种矿层成像方法、设备及计算机可读存储介质 |
| CN111476891A (zh) * | 2020-02-13 | 2020-07-31 | 山西晋煤集团技术研究院有限责任公司 | 一种井斜坐标数据批量处理的方法 |
| CN111551989A (zh) * | 2020-05-20 | 2020-08-18 | 中国科学院地理科学与资源研究所 | 一种透射槽波成像方法、设备及计算机可读存储介质 |
| CN111551989B (zh) * | 2020-05-20 | 2021-02-09 | 中国科学院地理科学与资源研究所 | 一种透射槽波成像方法、设备及计算机可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109597124B (zh) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109597124A (zh) | 一种工作面应力集中区槽波地震精细探测方法 | |
| CN102879805B (zh) | 一种基于钻孔与地面相结合的地震波空间探测方法 | |
| CN104331745B (zh) | 油气藏内天然裂缝的分期、分成因预测评价方法 | |
| CN112746837A (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集***及方法 | |
| Sil et al. | Dynamic site characterization and correlation of shear wave velocity with standard penetration test ‘N’values for the city of Agartala, Tripura state, India | |
| CN105629325B (zh) | 前陆盆地冲积扇精细刻画与预测方法 | |
| CN111273375B (zh) | 一种应用于缺水地区浅埋地下工程的地质勘察方法 | |
| CN106772554A (zh) | 一种用于复杂地形条件下的多道瞬态面波勘探方法 | |
| CN105277982A (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
| CN214576965U (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集*** | |
| CN107728205B (zh) | 一种地层压力预测方法 | |
| CN111812707A (zh) | 一种基于微动和面波勘探的地铁沿线场地内孤石探测方法 | |
| Ansal et al. | A preliminary microzonation study for the town of Dinar | |
| CN108931812B (zh) | 一种基于透射法槽波埃里相判识构造煤的方法 | |
| Brückl et al. | Seismic monitoring of deep-seated mass movements | |
| McCann et al. | Application of cross-hole seismic measurements in site investigation surveys | |
| Pegah et al. | Evaluating the overconsolidation ratios and peak friction angles of granular soil deposits using noninvasive seismic surveying | |
| Raines et al. | Refraction microtremor (ReMi) to determine the shear-wave velocity structure of the near surface and its application to aid detection of a backfilled mineshaft | |
| Aloisio et al. | Experimental dispersion curves of non-penetrable soils from direct dynamic measurements using the seismic dilatometer (SDMT) | |
| Wang et al. | A Rock Mass Strength Prediction Method Integrating Wave Velocity and Operational Parameters Based on the Bayesian Optimization Catboost Algorithm | |
| Naveen et al. | Evaluation of dynamic properties of municipal solid waste sites by geophysical tests | |
| Zhang et al. | P-wave velocity tomography and acoustic emission characteristics of sandstone under uniaxial compression | |
| Bianchi et al. | A Study of Progressive Failure in Porous Rocks Using Numerical and Experimental Modeling | |
| CN110068864B (zh) | 一种用于探测隧道上覆地层空洞及不密实的方法 | |
| Ladak et al. | Application of active and passive seismic methods for determining the shear wave velocity profile at hard rock site in Eastern Canada |
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 |