CN110531413B - 一种小断层超前三维可视化建模方法 - Google Patents
一种小断层超前三维可视化建模方法 Download PDFInfo
- Publication number
- CN110531413B CN110531413B CN201910773648.5A CN201910773648A CN110531413B CN 110531413 B CN110531413 B CN 110531413B CN 201910773648 A CN201910773648 A CN 201910773648A CN 110531413 B CN110531413 B CN 110531413B
- Authority
- CN
- China
- Prior art keywords
- fault
- point
- coordinate
- seismic
- advanced
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000000007 visual effect Effects 0.000 title claims abstract description 8
- 239000003245 coal Substances 0.000 claims abstract description 37
- 230000010287 polarization Effects 0.000 claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000012800 visualization Methods 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 101150072344 argA gene Proteins 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 17
- 230000005641 tunneling Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 241000219000 Populus Species 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
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
本发明公开了一种小断层超前三维可视化建模方法,先建立煤巷掘进超前预报观测系统;然后对各个三分量检波器接收到的三分量Love型反射槽波信号,进行时频域极化偏移成像,获取反射槽波在本盘界面上反射点位置及断层走向;在接收到的三分量信号中确定绕射体波信号,然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理,获取断层对盘绕射点的位置,进行收敛成像,并获得断层落差;通过上述得到的反射点及绕射点几何位置关系,获取小断层的倾角参数;最后获得断层性质及倾向,根据上述断层参数即可建立小断层三维空间模型。因此本发明能精确探测出巷道前方断层的三维空间下的各要素信息,从而能实现超前探测小断层的三维可视化。
Description
技术领域
本发明涉及一种矿井地球物理探测方法,具体是一种小断层超前三维可视化建模方法。
背景技术
随着煤矿开采深度的不断增加,三维地震探测断层的难度增加,导致目前对断层的探查程度不够,同时煤巷掘进过程中由于断层的影响造成的煤与瓦斯突出事故和突水频发,给煤矿生产带来巨大的安全隐患,因此积极发展井下探查技术是目前煤矿生产迫切需要。
煤巷掘进超前探测断层主要通过矿井物探方法(反射槽波超前探、瞬变电磁法,直流电法等)和钻探技术,钻探技术虽然能够直观的反应前方的构造情况,但存在“一孔之径”的问题所导致的探测范围的有限,且钻探成本较高,目前地震反射槽波最适合超前探测断层,目前针对回采工作面内断层的探测,前人做了大量工作,而针对掘进迎头前方断层探测较少。专利申请号“201711007670.6”中杨思通提出厚煤层沿底板或顶板掘进煤巷反射槽波超前探测方法,利用Rayleigh槽波进行超前探测,但Rayleigh槽波形成条件较为苛刻,而Love型槽波更具普遍性。专利申请号“201711391055.X”中王保利提出了基于孔巷联合的井下独头巷道槽波超前探测方法及系统,利用偏振滤波对地震记录进行区噪处理,采用散射成像方法对去噪后的地震记录进行偏移成像;综上所述,针对煤巷掘进迎头前方断层研究上取得了一定的成果,但均没有解决如何可视化探测断层的几何要素(如断层面的具体走向、倾向、倾角等)的问题。因此如何发明一种小断层超前三维可视化建模方法是本行业的研究方向。
发明内容
针对上述现有技术存在的问题,本发明提供一种小断层超前三维可视化建模方法,能精确探测出巷道前方断层,在三维空间中的各要素参数信息,从而能实现超前探测断层的三维可视化。
为了实现上述目的,本发明采用的技术方案是:一种小断层超前三维可视化建模方法,该方法的具体步骤为:
步骤一:在煤巷左右帮的中心高度分别布置n个三分量检波器接收地震信号,道间距(即相邻三分量检波器之间的距离)为N米,从迎头反向依此布置;
步骤二:采用CO2炮在煤巷左右帮煤层中各激发一次,分别产生震源S1和震源S2;(震源产生的原理为CO2炮激发后对煤岩体压裂过程中会产生震动)
步骤三:沿煤巷建立三维坐标系,X方向为指向迎头,Y方向为垂直于煤壁方向,Z方向为垂直于煤层底板,以模型边界上煤巷中心为原点建立三维空间坐标系,将上述震源位置与三分量检波器位置归入到空间坐标系中,建立三维煤巷地震超前预报观测系统;
步骤四:分析震源S1和S2的地震记录,选取信噪比高的地震记录进行处理分析;
步骤五:在地震记录中确定来自本盘的反射Love型槽波信号,进行时频域极化偏移成像,获取本盘断层的位置及走向,并获得最大振幅反射点位置;
步骤六:在地震记录中确定来自对盘的绕射体波信号,然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理,获取对盘绕射点的位置,进行收敛成像;
步骤七:根据步骤六中对盘绕射点位置求取断层的落差;
步骤八:结合步骤五、步骤六中反射点位置和绕射点位置,通过几何位置关系求取断层的倾角、性质、倾向;
步骤九:根据步骤五至八获取的断层的位置、走向、倾角、性质及倾向的参数信息,建立断层三维空间模型。
进一步,所述步骤五具体步骤为:
(1)在由t组成的时窗t1范围内的反射Love型数据Ai(t)采用解析信号ui(t)来表示:
(2)根据上述时频特征,构建协方差矩阵:
(3)通过上述协方差矩阵求取的最大特征值对应的特征向量,进行归一化后表示为[x1(t),y1(t),z1(t)];
(4)主极化方向的极化优势倾角表示为:
(5)根据以上参数进行时频域极化成像,其原理如下所示:
其中Ω代表所有地震道,A(t,f)代表标量振幅的叠加函数,Ai代表振幅能量,Ls+Lr代表地震波传播的总路径,v代表地震波速;主极化方向的权重函数Pc表达如下:
(6)通过上述(1)至(5),得到断层本盘界面,计算得出断层界面走向为α1和断层距离坐标原点的水平长度为L;
(7)在断层本盘界面上选取最大振幅反射点位置,设其为A点,三维空间坐标为(xi,yi,zi)。
进一步,所述步骤六的具体过程为:
(1)在地震记录中确定断层对盘断点绕射体波波列,选定时窗t2;
(2)从时窗t2中选取负相位最大振幅采样点,对该采样点的三分量信号进行时频域希尔伯特极化分析,获得主极化方向的方位角及倾角参数;
(3)通过n个检波器之间的组合计算关系,结合上步骤方位角及倾角参数,获得射线交汇点的空间分布;
(4)将上步骤交汇点的位置归入到三维空间坐标系内,然后将三维空间坐标系以1m为单位划分三维空间网格;
(5)统计每个三维空间网格内落入的交汇点,将落入交汇点最多的三维空间网格确定为对盘绕射点位置,设其点为B,坐标为(xj,yj,zj);
进一步,所述步骤七的具体过程为:
(1)将步骤六中对盘断层绕射点B投影到XOZ面上,即可获得B点坐标为(xj,zj);
(2)根据几何位置关系即可获得落差H为:
H=zj
进一步,所述步骤八的具体过程为:
(1)由于断层本盘反射点A的坐标为(xi,yi,zi)以及断层对盘绕射点B坐标为(xj,yj,zj);将A、B两点分别投影到XOZ面上,即坐标分别为(xi,zi),(xj,zj),根据坐标关系求取AB两点连线投影角度即:
θ即为断层倾角。
(2)在XOZ面上,根据几何位置坐标关系,判定断层的性质,具体步骤为:
根据上述A、B两点的位置坐标关系,即可获得断层的性质,可分为以下四种情况:
①当坐标关系满足zj>zi,xj>xi时,为正断层;
②当坐标关系满足zj<zi,xj<xi时,为逆断层;
③当坐标关系满足zj<zi,xj>xi时,为正断层;
④当坐标关系满足zj>zi,xj<xi时,为逆断层;
(3)根据上述步骤(2)得到的断层性质,结合确定的断层走向,即获得断层的倾向。
与现有技术相比,本发明先建立煤巷掘进超前预报观测系统;然后对各个三分量检波器接收到的三分量Love型反射槽波信号,进行时频域极化偏移成像,获取反射槽波在断层本盘界面上反射点位置及走向、倾向的大小;在接收到的三分量信号中确定绕射体波信号,然后对绕射体波信号进行时频域希尔伯特极化成像方法,获取断层对盘绕射点的位置,进行收敛成像,并获得小断层落差;通过上述得到的反射点及绕射点位置关系,通过几何关系,获取小断层的倾角大小;最后根据获取的断层参数,建立小断层三维空间模型。因此本发明能精确探测出巷道前方断层的三维空间下的各要素信息,从而能实现超前探测小断层的三维可视化。
附图说明
图1是本发明中检波器和震源的布设示意图;
图2是图3的俯视图;
图3是本发明探测断层时的剖面示意图;
图4是本发明中震源S1激发后的地震记录信号及反射槽波组和反射体波组图;
图5是本发明中震源S2激发后的地震记录信号及反射槽波组和反射体波组图;
图6是本发明三维坐标系中XOY面探测成像示意图;
图7是本发明三维坐标系中XOZ面探测成像示意图;
图8是本发明在XOZ面中计算断层落差的示意图;
图9是本发明在XOZ面中计算断层面倾角的示意图;
图10是本发明断层性质判断示意图。
图中:1、震源S1,2、三分量检波器,3、巷道,4、断层界面,5、断层本盘,6、断层对盘。
具体实施方式
下面将对本发明做进一步说明。
如图所示,以图1中的右方为前方进行专利描述,本发明的具体步骤为:
步骤一:在煤巷左右帮的中心高度分别布置n个三分量检波器2接收地震信号,道间距(即相邻三分量检波器之间的距离)为N米,从迎头反向依此布置;
步骤二:采用CO2炮在煤巷左右帮煤层中各激发一次,分别产生震源S1和震源S2;(震源产生的原理为CO2炮激发后对煤岩体压裂过程中会产生震动)
步骤三:沿煤巷建立三维坐标系,X方向为指向迎头,Y方向为垂直于煤壁方向,Z方向为垂直于煤层底板,以模型边界上煤巷中心为原点建立三维空间坐标系,将上述震源位置与三分量检波器2位置归入到空间坐标系中,建立三维煤巷地震超前预报观测系统;
步骤四:分析震源S1和S2的地震记录,选取信噪比高的地震记录进行处理分析;
步骤五:在地震记录中确定来自本盘的反射Love型槽波信号,进行时频域极化偏移成像,获取断层本盘5的位置及走向,并获得最大振幅反射点位置;
步骤六:在地震记录中确定来自对盘的绕射体波信号,然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理,获取对盘绕射点的位置,进行收敛成像;
步骤七:根据步骤六中对盘绕射点位置求取断层的落差;
步骤八:结合步骤五、步骤六中反射点位置和绕射点位置,通过几何位置关系求取断层的倾角、性质、倾向;
步骤九:根据步骤五至八获取的断层的位置、走向、倾角、性质及倾向的参数信息,建立断层三维空间模型。
进一步,所述步骤五具体步骤为:
(1)在由t组成的时窗t1范围内的反射Love型数据Ai(t)采用解析信号ui(t)来表示:
(2)根据上述时频特征,构建协方差矩阵:
(3)通过上述协方差矩阵求取的最大特征值对应的特征向量,进行归一化后表示为[x1(t),y1(t),z1(t)];
(4)主极化方向的极化优势倾角表示为:
(5)根据以上参数进行时频域极化成像,其原理如下所示:
其中Ω代表所有地震道,A(t,f)代表标量振幅的叠加函数,Ai代表振幅能量,Ls+Lr代表地震波传播的总路径,v代表地震波速;主极化方向的权重函数Pc表达如下:
(6)通过上述(1)至(5),得到断层本盘5界面,计算得出断层界面4走向为α1和断层距离坐标原点的水平长度为L;
(7)在断层本盘5界面上选取最大振幅反射点位置,设其为A点,三维空间坐标为(xi,yi,zi)。
进一步,所述步骤六的具体过程为:
(1)在地震记录中确定断层对盘断点绕射体波波列,选定时窗t2;
(2)从时窗t2中选取负相位最大振幅采样点,对该采样点的三分量信号进行时频域希尔伯特极化分析,获得主极化方向的方位角及倾角参数;
(3)通过n个检波器之间的组合计算关系,结合上步骤方位角及倾角参数,获得射线交汇点的空间分布;具体步骤为:
①根据任一个检波器的空间位置和步骤(2)获取主极化方向的方位角及倾角参数构建空间矢量,如下:
设检波器R的空间坐标为(xR,yR,zR),检波器R接收的绕射波的方位角、倾角的方向信息为(m,n,p),列出参数方程(t表示时间参数):
计算任意两个主极化方向射线相交的交点位置:
②在三维空间中得出两个相交的主极化方向射线,并将该交点的位置坐标记录;
③将剩下的检波器的空间坐标及检波器所接收的对盘断点绕射波传播方向重复步骤①和②,获得所有交汇点的位置;
④统计出各个交汇点在三维空间的分布特征;
(4)将上步骤交汇点的位置归入到三维空间坐标系内,然后将三维空间坐标系以1m为单位划分三维空间网格;
(5)统计每个三维空间网格内落入的交汇点,将落入交汇点最多的三维空间网格确定为对盘绕射点位置,设其点为B,坐标为(xj,yj,zj);
进一步,所述步骤七的具体过程为:
(1)将步骤六中对盘断层绕射点B投影到XOZ面上,即可获得B点坐标为(xj,zj);
(2)根据几何位置关系即可获得落差H为:
H=zj
进一步,所述步骤八的具体过程为:
(1)由于断层本盘反射点A的坐标为(xi,yi,zi)以及断层对盘绕射点B坐标为(xj,yj,zj);将A、B两点分别投影到XOZ面上,即坐标分别为(xi,zi),(xj,zj),根据坐标关系求取AB两点连线投影角度即:
θ即为断层倾角。
(2)在XOZ面上,根据几何位置坐标关系,判定断层的性质,具体步骤为:
根据上述A、B两点的位置坐标关系,即可获得断层的性质,可分为以下四种情况:
①当坐标关系满足zj>zi,xj>xi时,为正断层;
②当坐标关系满足zj<zi,xj<xi时,为逆断层;
③当坐标关系满足zj<zi,xj>xi时,为正断层;
④当坐标关系满足zj>zi,xj<xi时,为逆断层;
(3)根据上述步骤(2)得到的断层性质,结合确定的断层走向,即获得断层的倾向。
Claims (5)
1.一种小断层超前三维可视化建模方法,其特征在于,该方法的具体步骤为:
步骤一:在煤巷左右帮的中心高度分别布置n个三分量检波器接收地震信号,道间距为N米,从迎头反向依此布置;
步骤二:采用CO2炮在煤巷左右帮煤层中各激发一次,分别产生震源S1和震源S2;
步骤三:沿煤巷建立三维坐标系,X方向为指向迎头,Y方向为垂直于煤壁方向,Z方向为垂直于煤层底板,以模型边界上煤巷中心为原点建立三维空间坐标系,将震源位置与三分量检波器位置归入到空间坐标系中,建立三维煤巷地震超前预报观测系统;
步骤四:分析震源S1和S2的地震记录,选取信噪比高的地震记录进行处理分析;
步骤五:在地震记录中确定来自本盘的反射Love型槽波信号,进行时频域极化偏移成像,获取本盘断层的位置及走向,并获得最大振幅反射点位置;
步骤六:在地震记录中确定来自对盘的绕射体波信号,然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理,获取对盘绕射点的位置,进行收敛成像;
步骤七:根据步骤六中对盘绕射点位置求取断层的落差;
步骤八:结合步骤五、步骤六中反射点位置和绕射点位置,通过几何位置关系求取断层的倾角、性质、倾向;
步骤九:根据步骤五至八获取的断层的位置、走向、倾角、性质及倾向的参数信息,建立断层三维空间模型。
2.根据权利要求1所述的一种小断层超前三维可视化建模方法,其特征在于,所述步骤五具体步骤为:
(1)在由t组成的时窗t1范围内的反射Love型数据Ai(t)采用解析信号ui(t)来表示:
(2)构建协方差矩阵:
(3)通过上述协方差矩阵求取的最大特征值对应的特征向量,进行归一化后表示为[x1(t),y1(t),z1(t)];
(4)主极化方向的极化优势倾角表示为:
(5)根据以上参数进行时频域极化成像,其原理如下所示:
其中Ω代表所有地震道,A(t,f)代表标量振幅的叠加函数,Ai代表振幅能量,Ls+Lr代表地震波传播的总路径,v代表地震波速;主极化方向的权重函数Pc表达如下:
(6)通过上述步骤(1)至(5),得到断层本盘界面,计算得出断层界面走向为α1和断层距离坐标原点的水平长度为L;
(7)在断层本盘界面上选取最大振幅反射点位置,设其为A点,三维空间坐标为(xi,yi,zi)。
3.根据权利要求1所述的一种小断层超前三维可视化建模方法,其特征在于,所述步骤六的具体过程为:
(1)在地震记录中确定断层对盘断点绕射体波波列,选定时窗t2;
(2)从时窗t2中选取负相位最大振幅采样点,对该采样点的三分量信号进行时频域希尔伯特极化分析,获得主极化方向的方位角及倾角参数;
(3)通过n个检波器之间的组合计算关系,结合上步骤方位角及倾角参数,获得射线交汇点的空间分布;
(4)将上步骤交汇点的位置归入到三维空间坐标系内,然后将三维空间坐标系以1m为单位划分三维空间网格;
(5)统计每个三维空间网格内落入的交汇点,将落入交汇点最多的三维空间网格确定为对盘绕射点位置,设其点为B,坐标为(xj,yj,zj)。
4.根据权利要求1所述的一种小断层超前三维可视化建模方法,其特征在于,所述步骤七的具体过程为:
(1)将步骤六中对盘断层绕射点B投影到XOZ面上,即可获得B点坐标为(xj,zj);
(2)根据几何位置关系即可获得落差H为:
H=zj。
5.根据权利要求2所述的一种小断层超前三维可视化建模方法,其特征在于,所述步骤八的具体过程为:
(1)由于断层本盘反射点A的坐标为(xi,yi,zi)以及断层对盘绕射点B坐标为(xj,yj,zj);将A、B两点分别投影到XOZ面上,即坐标分别为(xi,zi),(xj,zj),根据坐标关系求取AB两点连线投影角度即:
θ即为断层倾角;
(2)在XOZ面上,根据几何位置坐标关系,判定断层的性质,具体步骤为:
根据上述A、B两点的位置坐标关系,即可获得断层的性质,分为以下四种情况:
①当坐标关系满足zj>zi,xj>xi时,为正断层;
②当坐标关系满足zj<zi,xj<xi时,为逆断层;
③当坐标关系满足zj<zi,xj>xi时,为正断层;
④当坐标关系满足zj>zi,xj<xi时,为逆断层;
(3)根据上述步骤(2)得到的断层性质,结合获取的走向,即获得断层的倾向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910773648.5A CN110531413B (zh) | 2019-08-21 | 2019-08-21 | 一种小断层超前三维可视化建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910773648.5A CN110531413B (zh) | 2019-08-21 | 2019-08-21 | 一种小断层超前三维可视化建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110531413A CN110531413A (zh) | 2019-12-03 |
CN110531413B true CN110531413B (zh) | 2020-10-30 |
Family
ID=68662461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910773648.5A Active CN110531413B (zh) | 2019-08-21 | 2019-08-21 | 一种小断层超前三维可视化建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110531413B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111965695B (zh) * | 2020-08-17 | 2023-04-25 | 山西潞安环保能源开发股份有限公司五阳煤矿 | 一种基于反射槽波的小断层落差探测方法 |
CN114185082B (zh) * | 2021-12-02 | 2023-04-21 | 中国矿业大学 | 一种基于工作面透射地震观测的煤层下伏陷落柱探测方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2032103B (en) * | 1978-10-20 | 1983-01-19 | Coal Industry Patents Ltd | Method and apparatus for locating flaws in a mineral seam |
US5060204A (en) * | 1990-06-27 | 1991-10-22 | Chevron Research And Technology Company | Method of layer stripping to determine fault plane stress build-up |
JP4818010B2 (ja) * | 2006-06-23 | 2011-11-16 | 常生 山内 | 地震時の建造物変形量に基づく地震マグニチュードの早期予測方法及び地震マグニチュードの早期予測プログラム |
WO2010019314A1 (en) * | 2008-08-11 | 2010-02-18 | Exxonmobil Upstream Research Company | Removal of surface-wave noise in seismic data |
US7978563B2 (en) * | 2009-08-18 | 2011-07-12 | Microseismic, Inc. | Method for passive seismic emission tomography including polarization correction for source mechanism |
CN102681004A (zh) * | 2012-05-14 | 2012-09-19 | 中国矿业大学(北京) | 以掘进机为震源的巷道随掘地震超前探测装置及方法 |
CN104166164B (zh) * | 2014-08-08 | 2017-08-29 | 山东科技大学 | 煤巷掘进地质构造三分量多波反射三维地震超前探测方法 |
CN105572745A (zh) * | 2015-12-10 | 2016-05-11 | 北京中矿大地地球探测工程技术有限公司 | 一种煤矿井下三分量槽波地震勘探方法 |
CN106875471B (zh) * | 2017-01-13 | 2019-12-06 | 山东科技大学 | 煤系含或隔水层三维可视化建模方法 |
CN109765606B (zh) * | 2019-01-28 | 2020-08-04 | 阳泉煤业(集团)有限责任公司 | 一种基于反射槽波的回采工作面隐蔽断层性质的探测方法 |
-
2019
- 2019-08-21 CN CN201910773648.5A patent/CN110531413B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110531413A (zh) | 2019-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Three-dimensional seismic ahead-prospecting method and application in TBM tunneling | |
CN102395902B (zh) | 使用快速面向目标照明计算的地震成像系统及方法 | |
Liu et al. | A new 3D observation system designed for a seismic ahead prospecting method in tunneling | |
US8233351B2 (en) | Analyzing borehole seismic data to locate subsurface diffractors | |
CN102681004A (zh) | 以掘进机为震源的巷道随掘地震超前探测装置及方法 | |
CN102749645B (zh) | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 | |
CN105700010B (zh) | 煤矿联合地震全息预测方法 | |
CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
CN110687602A (zh) | 浅层地震多波联合勘探方法 | |
CN110531413B (zh) | 一种小断层超前三维可视化建模方法 | |
CN101663597B (zh) | 根据二维地震数据进行三维散射体成像 | |
CN105676279A (zh) | 一种采用同心圆等炮检距的地震反射数据采集方法 | |
CN105223607A (zh) | 基于反稳相滤波的绕射波成像方法 | |
CN110632667B (zh) | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 | |
CN103235334B (zh) | 基于地震信号的低抽巷道超前探测断层落差方法 | |
CN110687607B (zh) | 一种斯通利波探测方法及系统 | |
EP2593815B1 (en) | Method for accentuating specular and non-specular seismic events from within shallow subsurface rock formations | |
Ji et al. | Three-component polarization migration of channel waves for prediction ahead of coal roadway | |
Song et al. | Experiment of 3D seismic reflection technique for forward probing on TBM tunnel Face | |
CN102798884B (zh) | 巷道顶板二维地震勘探方法及系统 | |
CN106772637A (zh) | 一种基于槽波超前探测tvsp扫描方法 | |
CN110531416B (zh) | 一种基于三分量反射信号时频域极化参数的断层确定方法 | |
US3354985A (en) | Seismic prospecting method of locating a suspected fault | |
CN102478663A (zh) | 一种三维地震观测系统偏移噪声获取方法及装置 | |
Wang et al. | A polarization migration velocity model building method for geological prediction ahead of the tunnel face |
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 |