CN110632667A - 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 - Google Patents
一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 Download PDFInfo
- Publication number
- CN110632667A CN110632667A CN201910993568.0A CN201910993568A CN110632667A CN 110632667 A CN110632667 A CN 110632667A CN 201910993568 A CN201910993568 A CN 201910993568A CN 110632667 A CN110632667 A CN 110632667A
- Authority
- CN
- China
- Prior art keywords
- point
- seismic source
- shock wave
- seismic
- points
- 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
- 230000035939 shock Effects 0.000 title claims abstract description 46
- 238000001514 detection method Methods 0.000 title claims abstract description 17
- 230000005284 excitation Effects 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 30
- 239000003245 coal Substances 0.000 claims abstract description 16
- 230000005641 tunneling Effects 0.000 claims abstract description 6
- 238000003384 imaging method Methods 0.000 claims abstract description 4
- 230000010287 polarization Effects 0.000 claims description 36
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005516 engineering process Methods 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/02—Generating seismic energy
-
- 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
-
- 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
本发明公开了一种基于冲击波激发激震条件下的隐伏陷落柱超前探测方法,采用冲击波作为激发震源,将冲击波激发震源放置到底板钻孔最深处,通过改变不同孔深的方式分别进行激发,从而产生多个震源,且多个震源与检波器之间的偏移距均不相同,每次激发后获得一次共炮点道集,然后将各个共炮点道集叠加进行处理,获取各个震源产生的绕射波在隐伏陷落柱上最大振幅绕射点位置;最后通过各个最大振幅绕射点位置进行曲线拟合,将曲线显示在坐标系的XOZ面上,该曲线的位置即为隐伏陷落柱的侧界面位置;本发明通过冲击波激发震源进行多次激发,从而利用绕射波信号对煤层下方的陷落柱进行探测成像,不仅精度高,而且方便快捷,对煤巷安全掘进提供数据支撑。
Description
技术领域
本发明涉及一种隐伏陷落柱的超前探测方法,具体是一种基于冲击波激震条件下的隐伏陷落柱超前探测方法。
背景技术
巷道掘进是煤矿安全事故的频发阶段,有些巷道掘进过程中会发生突水事故,使煤矿受到巨大的经济损失或人员伤亡,突水事故的主要原因是底板隐伏陷落柱导通奥灰水形成。目前针对掘进迎头前方煤层的底板隐伏陷落柱的探测方法有限,因此亟需发展井下探测技术,确保工作面安全掘进。
专利号“201610323984.6”中提出一种基于绕射波的回采工作面小型陷落柱的探测方法,利用工作面煤层内的陷落柱透射方向的绕射波进行探查,然而超前探测是利用反射方向的绕射波,两者差异大,因此无法采用;故如何对隐伏陷落柱进行超前探测是本行业的研究方向。
发明内容
针对上述现有技术存在的问题,本发明提供一种基于冲击波激震条件下的隐伏陷落柱超前探测方法,通过冲击波激震产生多次激发,从而利用底板绕射波信号对底板以下的陷落柱进行探测成像,不仅精度高,而且使用方便,对煤巷安全掘进提供数据支撑。
为了实现上述目的,本发明采用的技术方案是:一种基于冲击波激发震源条件下的隐伏陷落柱超前探测方法,该方法的具体步骤为:
步骤一:在巷道底板距离迎头10m位置施工一个激发钻孔,钻孔长度为N米,钻孔与巷道底板垂直;
步骤二:将冲击波激发震源送入钻孔最深处;
步骤三:在巷道底板布设n个安装孔,相邻安装孔的间距为3米,n个安装孔与激发钻孔的位置处在同一直线上,在每个安装孔内均布置一个三分量检波器;各个三分量检波器与地震记录仪连接,建立煤巷掘进超前预报观测***;
步骤四:冲击波激发震源从钻孔最深处进行第一次激发,然后冲击波激发震源由下往上每偏移3m激发一次,直至冲击波激发震源距离钻孔孔口3m时,停止激发并退出钻孔,共激发n次,形成震源S1、S2、S3…Sn,各个三分量检波器实时接收各个冲击波激发震源产生的地震波并传递给地震记录仪;
步骤五:沿巷道建立三维坐标系,其中X方向为指向迎头方向,Y方向为垂直于煤壁且垂直于迎头方向,Z方向为垂直于巷道顶板且垂直于迎头方向,以模型左边界上巷道面的中心为原点建立三维坐标系,将各个震源位置与三分量检波器位置归入上述三维坐标系中;
步骤六:每个震源产生一个共炮点道集,共产生n个共炮点道集,对每个共炮点道集进行分析,选取其中来自前方的绕射波信号,对绕射波信号进行时频域希尔伯特极化成像,计算煤层下方隐伏陷落柱上绕射点的位置;
步骤七:获取每个共炮点道集中最大振幅绕射点,记为A1(x1,y1,z1),A2(x2,y2,z2),A3(x3,y3,z3)…An(xn,yn,zn);
步骤八:将所有最大振幅绕射点反映在三维坐标系的XOZ面上,其坐标分别为A1(x1,z1)、A2(x2,z2),A3(x3,z3)…….An(xn,zn);
步骤九:将步骤八中的各个绕射点进行曲线拟合,得出的曲线即为隐伏陷落柱的侧界面。
进一步,所述步骤六具体为:
(1)在震源S1形成的共炮点道集中确定来自底板的绕射波信号,选定带有绕射波信号的时窗t;
(2)通过对绕射波信号进行时频域希尔伯特极化分析,获得绕射波的方位角及倾角参数;
(3)通过检波器组合计算关系结合极化参数并三维网格化后,获得射线交汇点的空间分布,具体步骤为:
Ⅰ、根据任一个检波器的空间位置和步骤(2)获取相应的极化参数构建空间矢量,如下:
设检波器R的空间坐标为(xR,yR,zR),检波器R接收的绕射波的方位角、倾角的主极化方向信息为(m,l,p),列出该检波器的主极化方向参数方程,t表示时间参数:
计算任意两个检波器主极化方向射线相交的交点位置:
上述方程可用以下矩阵形式表示:
VT=P
其中:
Ⅱ、在三维空间中若能得出两个相交的主极化方向射线,则记录该交点的位置坐标;地震记录由于噪声干扰和极化方法误差的影响,若未能得出两个相交的主极化方向射线,则使目标函数的解J达到最小:
J=||T-V-1P||2
采用最小二乘法求解,即:
T=(VTV)-1VTP
从而得出两个未相交的主极化方向射线的最小距离,并记录该最小距离中点的位置坐标;
Ⅲ、重复步骤Ⅰ和Ⅱ将各个检波器的主极化方向射线两两分组,使每组的两个主极化方向射线相交或距离最小;并得出各组的交点或中点的位置坐标,这些坐标则为所有绕射点的位置坐标;
(4)对震源S2、S3…Sn形成的三分量共炮点道集重复步骤(1)至(3),获得震源S2、S3…Sn各自所有绕射点的位置坐标。
进一步,所述步骤七的具体过程为:
(1)统计震源S1形成的各个绕射点在三维空间的分布特征,将各组得出的交点或中点的位置坐标放入在建立的三维空间坐标系内,然后将三维空间坐标系以1m为单位划分三维空间网格;
(2)统计每个三维空间网格内落入的交点或中点,将落入交点或中点数量最多的三维空间网格确定为最大振幅绕射点位置,设其点为A1,坐标为(x1,y1,z1);
(3)对震源S2、S3…Sn形成的各个绕射点重复步骤(1)和(2),分别得出各个震源对应的最大振幅绕射点位置分别为A2(x2,y2,z2),A3(x3,y3,z3),A4(x4,y4,z4)…An(xn,yn,zn)。
进一步,所述步骤九的具体过程为:
(1)各个最大振幅绕射点在XOZ面上的坐标分别为A1(x1,z1)、A2(x2,z2),A3(x3,z3)…….An(xn,zn);
(2)设拟合多项式曲线为:
z=a0+a1x+…+akxk
其中a0、a1…ak为求解参数;
(3)将所有最大振幅绕射点坐标分别代入到步骤(2)中的拟合多项式,构成如下矩阵:
(4)设置k=1,k=2,k=3,k=4,分别代入步骤(3)中矩阵,利用最小二乘法进行求解,即可对应的获得四组参数;
(5)为了确定最优解,将四组a0、a1…ak的参数值代入拟合曲线,分别归入到三维坐标系的XOZ面上;
与现有技术相比,本发明采用冲击波激震仪作为震源,将冲击波震源放置到底板钻孔内,利用冲击波在钻孔内能多次激发的特点,通过改变不同孔深的方式分别进行激发,从而产生多个震源,且多个震源与检波器之间的偏移距均不相同,每次激发后地震观测***获取一次共炮点道集,然后将各个共炮点道集进行处理,获取各个震源产生的绕射波在隐伏陷落柱上最大振幅反射点位置;最后通过各个最大振幅反射点位置拟合曲线,将曲线显示在坐标系的XOZ面上,该曲线的位置即为隐伏陷落柱的反射界面;本发明通过冲击波激发震源产生多次激发,震源新型且多次循环激发效率快,解决了现有钻探、三维地震勘探精度低的难题。
附图说明
图1是本发明中观测***的布置示意图;
图2是采用本发明得出的反射点成像散点分布图。
图中:1、三分量检波器,2、巷道,3、迎头,4、钻孔,5、隐伏陷落柱。
具体实施方式
下面将对本发明做进一步说明。
如图所示,本发明的具体步骤为:
步骤一:在巷道2底板距离迎头10m位置施工一个激发钻孔,钻孔长度为N米,钻孔4与巷道2底板垂直;
步骤二:将冲击波激发震源送入钻孔4最深处;
步骤三:在巷道2底板布设n个安装孔,相邻安装孔的间距为3米,n个安装孔与激发钻孔的位置处在同一直线上,在每个安装孔内均布置一个三分量检波器1;各个三分量检波器1与地震记录仪连接,建立煤巷掘进超前预报观测***;
步骤四:冲击波激发震源从钻孔4最深处进行第一次激发,然后冲击波激发震源由下往上每偏移3m激发一次,直至冲击波激发震源距离钻孔孔口3m时,停止激发并退出钻孔4,共激发n次,形成震源S1、S2、S3…Sn,各个三分量检波器1实时接收各个冲击波激发震源产生的地震波并传递给地震记录仪;
步骤五:沿巷道建立三维坐标系,其中X方向为指向迎头方向,Y方向为垂直于煤壁且垂直于迎头方向,Z方向为垂直于巷道顶板且垂直于迎头方向,以模型左边界上巷道面的中心为原点建立三维坐标系,将各个震源位置与三分量检波器1位置归入上述三维坐标系中;
步骤六:每个震源产生一个共炮点道集,共产生n个共炮点道集,对每个共炮点道集进行分析,选取其中来自前方的绕射波信号,对绕射波信号进行时频域希尔伯特极化成像,计算煤层下方隐伏陷落柱5上绕射点的位置;
步骤七:获取每个共炮点道集中最大振幅绕射点,记为A1(x1,y1,z1),A2(x2,y2,z2),A3(x3,y3,z3)…An(xn,yn,zn);
步骤八:将所有最大振幅绕射点反映在三维坐标系的XOZ面上,其坐标分别为A1(x1,z1)、A2(x2,z2),A3(x3,z3)…….An(xn,zn);
步骤九:将步骤八中的各个绕射点进行曲线拟合,得出的曲线即为隐伏陷落柱5的侧界面。
进一步,所述步骤六具体为:
(1)在震源S1形成的共炮点道集中确定来自底板的绕射波信号,选定带有绕射波信号的时窗t;
(2)通过对绕射波信号进行时频域希尔伯特极化分析,获得绕射波的方位角及倾角参数;
(3)通过检波器组合计算关系结合极化参数并三维网格化后,获得射线交汇点的空间分布,具体步骤为:
Ⅰ、根据任一个检波器的空间位置和步骤(2)获取相应的极化参数构建空间矢量,如下:
设检波器R的空间坐标为(xR,yR,zR),检波器R接收的绕射波的方位角、倾角的主极化方向信息为(m,l,p),列出该检波器的主极化方向参数方程,t表示时间参数:
计算任意两个检波器主极化方向射线相交的交点位置:
上述方程可用以下矩阵形式表示:
VT=P
其中:
Ⅱ、在三维空间中若能得出两个相交的主极化方向射线,则记录该交点的位置坐标;地震记录由于噪声干扰和极化方法误差的影响,若未能得出两个相交的主极化方向射线,则使目标函数的解J达到最小:
J=||T-V-1P||2
采用最小二乘法求解,即:
T=(VTV)-1VTP
从而得出两个未相交的主极化方向射线的最小距离,并记录该最小距离中点的位置坐标;
Ⅲ、重复步骤Ⅰ和Ⅱ将各个检波器的主极化方向射线两两分组,使每组的两个主极化方向射线相交或距离最小;并得出各组的交点或中点的位置坐标,这些坐标则为所有绕射点的位置坐标;
(4)对震源S2、S3…Sn形成的三分量共炮点道集重复步骤(1)至(3),获得震源S2、S3…Sn各自所有绕射点的位置坐标。
进一步,所述步骤七的具体过程为:
(1)统计震源S1形成的各个绕射点在三维空间的分布特征,将各组得出的交点或中点的位置坐标放入在建立的三维空间坐标系内,然后将三维空间坐标系以1m为单位划分三维空间网格;
(2)统计每个三维空间网格内落入的交点或中点,将落入交点或中点数量最多的三维空间网格确定为最大振幅绕射点位置,设其点为A1,坐标为(x1,y1,z1);
(3)对震源S2、S3…Sn形成的各个绕射点重复步骤(1)和(2),分别得出各个震源对应的最大振幅绕射点位置分别为A2(x2,y2,z2),A3(x3,y3,z3),A4(x4,y4,z4)…An(xn,yn,zn)。
进一步,所述步骤九的具体过程为:
(1)各个最大振幅绕射点在XOZ面上的坐标分别为A1(x1,z1)、A2(x2,z2),A3(x3,z3)…….An(xn,zn);
(2)设拟合多项式曲线为:
z=a0+a1x+…+akxk
其中a0、a1…ak为求解参数;
(3)将所有最大振幅绕射点坐标分别代入到步骤(2)中的拟合多项式,构成如下矩阵:
(4)设置k=1,k=2,k=3,k=4,分别代入步骤(3)中矩阵,利用最小二乘法进行求解,即可对应的获得四组参数;
(5)为了确定最优解,将四组a0、a1…ak的参数值代入拟合曲线,分别归入到三维坐标系的XOZ面上;
Claims (4)
1.一种基于冲击波激发震源条件下的隐伏陷落柱超前探测方法,其特征在于,该方法的具体步骤为:
步骤一:在巷道底板距离迎头10m位置施工一个激发钻孔,钻孔长度为N米,钻孔与巷道底板垂直;
步骤二:将冲击波激发震源送入钻孔最深处;
步骤三:在巷道底板布设n个安装孔,相邻安装孔的间距为3米,n个安装孔与激发钻孔的位置处在同一直线上,在每个安装孔内均布置一个三分量检波器;各个三分量检波器与地震记录仪连接,建立煤巷掘进超前预报观测***;
步骤四:冲击波激发震源从钻孔最深处进行第一次激发,然后冲击波激发震源由下往上每偏移3m激发一次,直至冲击波激发震源距离钻孔孔口3m时,停止激发并退出钻孔,共激发n次,形成震源S1、S2、S3...Sn,各个三分量检波器实时接收各个冲击波激发震源产生的地震波并传递给地震记录仪;
步骤五:沿巷道建立三维坐标系,其中X方向为指向迎头方向,Y方向为垂直于煤壁且垂直于迎头方向,Z方向为垂直于巷道顶板且垂直于迎头方向,以模型左边界上巷道面的中心为原点建立三维坐标系,将各个震源位置与三分量检波器位置归入上述三维坐标系中;
步骤六:每个震源产生一个共炮点道集,共产生n个共炮点道集,对每个共炮点道集进行分析,选取其中来自前方的绕射波信号,对绕射波信号进行时频域希尔伯特极化成像,计算煤层下方隐伏陷落柱上绕射点的位置;
步骤七:获取每个共炮点道集中最大振幅绕射点,记为A1(x1,y1,z1),A2(x2,y2,z2),A3(x3,y3,z3)…An(xn,yn,zn);
步骤八:将所有最大振幅绕射点反映在三维坐标系的XOZ面上,其坐标分别为A1(x1,z1)、A2(x2,z2),A3(x3,z3).......An(xn,zn);
步骤九:将步骤八中的各个绕射点进行曲线拟合,得出的曲线即为隐伏陷落柱的侧界面。
2.根据权利要求1所述的一种基于冲击波激发震源条件下的隐伏陷落柱超前探测方法,其特征在于,所述步骤六具体为:
(1)在震源S1形成的共炮点道集中确定来自底板的绕射波信号,选定带有绕射波信号的时窗t;
(2)通过对绕射波信号进行时频域希尔伯特极化分析,获得绕射波的方位角及倾角参数;
(3)通过检波器组合计算关系结合极化参数并三维网格化后,获得射线交汇点的空间分布,具体步骤为:
I、根据任一个检波器的空间位置和步骤(2)获取相应的极化参数构建空间矢量,如下:
设检波器R的空间坐标为(xR,yR,zR),检波器R接收的绕射波的方位角、倾角的主极化方向信息为(m,l,p),列出该检波器的主极化方向参数方程,t表示时间参数:
计算任意两个检波器主极化方向射线相交的交点位置:
上述方程可用以下矩阵形式表示:
VT=P
其中:
II、在三维空间中若能得出两个相交的主极化方向射线,则记录该交点的位置坐标;地震记录由于噪声干扰和极化方法误差的影响,若未能得出两个相交的主极化方向射线,则使目标函数的解J达到最小:
J=||T-V-1P||2
采用最小二乘法求解,即:
T=(VTV)-1VTP
从而得出两个未相交的主极化方向射线的最小距离,并记录该最小距离中点的位置坐标;
III、重复步骤I和II将各个检波器的主极化方向射线两两分组,使每组的两个主极化方向射线相交或距离最小;并得出各组的交点或中点的位置坐标,这些坐标则为所有绕射点的位置坐标;
(4)对震源S2、S3…Sn形成的三分量共炮点道集重复步骤(1)至(3),获得震源S2、S3…Sn各自所有绕射点的位置坐标。
3.根据权利要求1所述的一种基于冲击波激发震源条件下的隐伏陷落柱超前探测方法,其特征在于,所述步骤七的具体过程为:
(1)统计震源S1形成的各个绕射点在三维空间的分布特征,将各组得出的交点或中点的位置坐标放入在建立的三维空间坐标系内,然后将三维空间坐标系以1m为单位划分三维空间网格;
(2)统计每个三维空间网格内落入的交点或中点,将落入交点或中点数量最多的三维空间网格确定为最大振幅绕射点位置,设其点为A1,坐标为(x1,y1,z1);
(3)对震源S2、S3…Sn形成的各个绕射点重复步骤(1)和(2),分别得出各个震源对应的最大振幅绕射点位置分别为A2(x2,y2,z2),A3(x3,y3,z3),A4(x4,y4,z4)…An(xn,yn,zn)。
4.根据权利要求1所述的一种基于冲击波激发震源条件下的隐伏陷落柱超前探测方法,其特征在于,所述步骤九的具体过程为:
(1)各个最大振幅绕射点在XOZ面上的坐标分别为A1(x1,z1)、A2(x2,z2),A3(x3,z3).......An(xn,zn);
(2)设拟合多项式曲线为:
z=a0+a1x+…+akxk
其中a0、a1…ak为求解参数;
(3)将所有最大振幅绕射点坐标分别代入到步骤(2)中的拟合多项式,构成如下矩阵:
(4)设置k=1,k=2,k=3,k=4,分别代入步骤(3)中矩阵,利用最小二乘法进行求解,即可对应的获得四组参数;
(5)为了确定最优解,将四组a0、a1…ak的参数值代入拟合曲线,分别归入到三维坐标系的XOZ面上;
(7)比较的大小,所有最大振幅绕射点与曲线的距离平方和最短代表曲线拟合最好,因此最小值对应的为最优拟合曲线,该曲线对应位置即为隐伏陷落柱的侧界面位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910993568.0A CN110632667B (zh) | 2019-10-18 | 2019-10-18 | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910993568.0A CN110632667B (zh) | 2019-10-18 | 2019-10-18 | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110632667A true CN110632667A (zh) | 2019-12-31 |
CN110632667B CN110632667B (zh) | 2020-12-18 |
Family
ID=68976692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910993568.0A Active CN110632667B (zh) | 2019-10-18 | 2019-10-18 | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110632667B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113296411A (zh) * | 2021-05-26 | 2021-08-24 | 广东电网有限责任公司 | 跟踪高频噪声功率增益调节控制策略参数的方法及*** |
CN113687308A (zh) * | 2021-09-07 | 2021-11-23 | 哈尔滨工程大学 | 一种基于弯曲波的冰上震源定位方法 |
CN114185082A (zh) * | 2021-12-02 | 2022-03-15 | 中国矿业大学 | 一种基于工作面透射地震观测的煤层下伏陷落柱探测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4935904A (en) * | 1989-06-08 | 1990-06-19 | Western Atlas International, Inc. | Method for suppressing coherent noise in seismic data |
WO2005073688A1 (en) * | 2004-01-21 | 2005-08-11 | California Institute Of Technology | Extreme event performance evaluation of a structure using real-time hystersis monitoring |
CN104166164A (zh) * | 2014-08-08 | 2014-11-26 | 山东科技大学 | 煤巷掘进地质构造三分量多波反射三维地震超前探测方法 |
US20150063066A1 (en) * | 2013-09-05 | 2015-03-05 | William A. Burnett | Creating Seismic Images Using Expanded Image Gathers |
CN104594901A (zh) * | 2014-12-08 | 2015-05-06 | 太原理工大学 | 一种工作面通过火成岩侵入区的方法 |
CN105911588A (zh) * | 2016-05-16 | 2016-08-31 | 中国矿业大学 | 一种基于绕射波的回采工作面小型陷落柱探测方法 |
US20180164452A1 (en) * | 2016-12-14 | 2018-06-14 | Pgs Geophysical As | Generating Pseudo Pressure Wavefields Utilizing a Warping Attribute |
CN108845358A (zh) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | 断层及构造异常体识别方法及装置 |
CN109387834A (zh) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | 地下水库坝体损伤探测方法及电子设备 |
CN109541690A (zh) * | 2018-11-30 | 2019-03-29 | 中铁第四勘察设计院集团有限公司 | 一种浅层介质结构面松散程度评价方法 |
-
2019
- 2019-10-18 CN CN201910993568.0A patent/CN110632667B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4935904A (en) * | 1989-06-08 | 1990-06-19 | Western Atlas International, Inc. | Method for suppressing coherent noise in seismic data |
WO2005073688A1 (en) * | 2004-01-21 | 2005-08-11 | California Institute Of Technology | Extreme event performance evaluation of a structure using real-time hystersis monitoring |
US20150063066A1 (en) * | 2013-09-05 | 2015-03-05 | William A. Burnett | Creating Seismic Images Using Expanded Image Gathers |
CN104166164A (zh) * | 2014-08-08 | 2014-11-26 | 山东科技大学 | 煤巷掘进地质构造三分量多波反射三维地震超前探测方法 |
CN104594901A (zh) * | 2014-12-08 | 2015-05-06 | 太原理工大学 | 一种工作面通过火成岩侵入区的方法 |
CN105911588A (zh) * | 2016-05-16 | 2016-08-31 | 中国矿业大学 | 一种基于绕射波的回采工作面小型陷落柱探测方法 |
US20180164452A1 (en) * | 2016-12-14 | 2018-06-14 | Pgs Geophysical As | Generating Pseudo Pressure Wavefields Utilizing a Warping Attribute |
CN109387834A (zh) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | 地下水库坝体损伤探测方法及电子设备 |
CN108845358A (zh) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | 断层及构造异常体识别方法及装置 |
CN109541690A (zh) * | 2018-11-30 | 2019-03-29 | 中铁第四勘察设计院集团有限公司 | 一种浅层介质结构面松散程度评价方法 |
Non-Patent Citations (1)
Title |
---|
赵禄顺 等: ""岩溶陷落柱绕射成像研究"", 《能源技术与管理》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113296411A (zh) * | 2021-05-26 | 2021-08-24 | 广东电网有限责任公司 | 跟踪高频噪声功率增益调节控制策略参数的方法及*** |
CN113687308A (zh) * | 2021-09-07 | 2021-11-23 | 哈尔滨工程大学 | 一种基于弯曲波的冰上震源定位方法 |
CN113687308B (zh) * | 2021-09-07 | 2023-12-19 | 哈尔滨工程大学 | 一种基于弯曲波的冰上震源定位方法 |
CN114185082A (zh) * | 2021-12-02 | 2022-03-15 | 中国矿业大学 | 一种基于工作面透射地震观测的煤层下伏陷落柱探测方法 |
CN114185082B (zh) * | 2021-12-02 | 2023-04-21 | 中国矿业大学 | 一种基于工作面透射地震观测的煤层下伏陷落柱探测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110632667B (zh) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110632667B (zh) | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 | |
US10197691B2 (en) | Acoustic anisotropy and imaging by means of high resolution azimuthal sampling | |
Liu et al. | A new 3D observation system designed for a seismic ahead prospecting method in tunneling | |
CN105676279B (zh) | 一种采用同心圆等炮检距的地震反射数据采集方法 | |
CN112485823B (zh) | 高效综合超前地质预报方法 | |
CN106443766A (zh) | 三维隧道地震超前探测方法 | |
CN109765606B (zh) | 一种基于反射槽波的回采工作面隐蔽断层性质的探测方法 | |
Zhao et al. | Inversion of seepage channels based on mining-induced microseismic data | |
CN110687602A (zh) | 浅层地震多波联合勘探方法 | |
CN104181581A (zh) | 基于任意排布的地震波地下工程空间观测的***及方法 | |
CN112578428A (zh) | 一种基于巷道垂向虚拟测线的散射多波超前探测方法 | |
CN114839672A (zh) | 一种超前水平钻探随钻随测的隧道快速超前地质预报方法 | |
Nie et al. | TBM-mounted seismic ahead-prospecting for fast detecting anomalous geology ahead of tunnel face | |
CN114280669A (zh) | 一种基于折射波周期振幅衰减的薄煤带探测方法及*** | |
CN113126144A (zh) | 一种隧道三维超前地质预报***和方法 | |
CN117192615A (zh) | 基于透射地震尾波的采煤工作面内隐伏地质构造探测方法 | |
CN110850472B (zh) | 一种基于冲击波激发震源的可变偏移距超前探测断层方法 | |
CN110531413B (zh) | 一种小断层超前三维可视化建模方法 | |
CN111852467B (zh) | 一种砂岩铀矿矿体延伸范围的圈定方法及*** | |
CN111158050B (zh) | 数据采集***、方法及隧道地震波超前预报方法 | |
CN116626757A (zh) | 一种隧道超前地质预报方法 | |
CN114185082B (zh) | 一种基于工作面透射地震观测的煤层下伏陷落柱探测方法 | |
CN215180925U (zh) | 一种隧道三维超前地质预报数据采集*** | |
CN112684504B (zh) | 一种基于全散射模型的城市地下空洞快速探测方法 | |
CN111025383B (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 |