CN111965695A - 一种基于反射槽波的小断层落差探测方法 - Google Patents
一种基于反射槽波的小断层落差探测方法 Download PDFInfo
- Publication number
- CN111965695A CN111965695A CN202010827051.7A CN202010827051A CN111965695A CN 111965695 A CN111965695 A CN 111965695A CN 202010827051 A CN202010827051 A CN 202010827051A CN 111965695 A CN111965695 A CN 111965695A
- Authority
- CN
- China
- Prior art keywords
- wave
- reflected
- fault
- seismic source
- frequency
- 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 23
- 239000003245 coal Substances 0.000 claims abstract description 36
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000005520 cutting process Methods 0.000 claims abstract description 4
- 239000000839 emulsion Substances 0.000 claims abstract description 4
- 239000002360 explosive Substances 0.000 claims abstract description 4
- 238000001228 spectrum Methods 0.000 claims description 18
- 238000000605 extraction Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000005065 mining Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 4
- 238000000034 method Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000002679 ablation Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000005641 tunneling 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/282—Application of seismic models, synthetic seismograms
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
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
本发明公开了一种基于反射槽波的小断层落差探测方法,在煤巷单巷道内探测断层区域的煤层中间高度等间隔布置n个横波检波器,采用矿用乳化***在探测起始点激发一次,产生震源;分析震源记录,确定反射槽波位置,进行直达波切除处理;对切除直达波后的震源记录进行主振幅提取,获取X分量反射槽波振幅;对X分量反射槽波进行滤波,获取滤波后X分量反射槽波能量;对比分析获得断层的能量变化规律,根据数据分析断层位置;提取反射槽波的子波类型,确定其频率大小;根据所得的信息建立不同落差断层的三维空间模型,得到不同落差断层反射槽波的能量变化规律,达到探测小断层落差的目的。
Description
技术领域
本发明涉及地震数据分析技术领域,尤其是一种基于反射槽波的小断层 落差探测方法。
背景技术
煤矿生产普遍使用大型采煤机,大型回采面越来越多,工作面回采安全 问题成为煤矿发展日益关注的焦点。
槽波勘探对断层、陷落柱等有害地质体能够起到一定的勘查效果,比如 识别煤层的走向,但是在预估断层落差大小方面仍有待探索和研究。槽波勘 探利用煤层对地震波的制导作用,对煤层内各种地质构造进行有效勘查。断 层是煤巷掘进过程中潜伏的安全隐患,可能造成瓦斯突出或突水事故,因此 探寻一种可行的断层分析方法十分必要。目前,我国对煤层中断层的研究有 很多,但是对断层落差,尤其是中小断层落差的探测分析研究较少。
发明内容
针对上述问题,本发明提供一种基于反射槽波的小断层落差探测方法, 用于探测煤矿中的中小断层落差。
一种基于反射槽波的小断层落差探测方法,包括以下步骤:
步骤1,在煤巷单巷道内探测断层区域建立三维煤巷单巷道观测***,观 测***主要由震源、等间隔分布的多个检波器、以及为震源和检波器提供位 置坐标的三维直角空间坐标系构成;
步骤2,分析震源记录,确定反射槽波位置,进行直达波切除处理;
步骤3,对切除直达波后的震源记录进行主振幅提取,获取X分量反射 槽波振幅;
步骤4,对X分量反射槽波进行滤波,获取滤波后X分量反射槽波能量;
步骤5,对比分析获得断层的能量变化规律,根据数据分析断层位置;
步骤6,提取反射槽波的子波类型,确定其频率大小;
步骤7,根据步骤5和步骤6所得的信息建立不同落差断层的三维空间模 型,得到不同落差断层反射槽波的能量变化规律,达到探测小断层落差的目 的。
进一步的,所述步骤1具体包括以下步骤:
步骤A1,在煤巷单巷道内探测断层区域的煤层中间高度等间隔布置n个 横波检波器;
步骤A2,采用矿用乳化***在探测起始点激发一次,产生震源;
步骤A3,沿煤巷建立三维直角空间坐标系,X方向指向待探测区域,Y 方向垂直于煤壁,Z方向垂直于煤层顶板;
步骤A4,确定震源、n个检波器在三维直角空间坐标中的坐标位置。
进一步的,所述步骤4包括以下步骤:
步骤B1,获得反射槽波频谱信息;
步骤B2,截取较高能量频段的频宽;
步骤B3,对截取的反射槽波进行频带宽度为50Hz的带通滤波;
步骤B4,对偏移距进行处理得到反射槽波入射角,定义为横坐标;
步骤B5,对滤波后的反射槽波数据进行最大振幅提取。
进一步的,所述步骤6包括以下步骤:
步骤C4,根据频率能量分布确定使用频率大小。
本发明通过建立三维煤巷单巷道观测***获取相关信息,根据所得的信 息建立不同落差断层的三维空间模型,得到不同落差断层反射槽波的能量变 化规律,达到探测小断层落差的目的。
附图说明
图1是震源和检波器布置示意图;
图2是反射槽波示意图;
图3是切波前的模拟数据图;
图4是切波后的模拟数据图;
图5是反射槽波子波形态图;
图6是反射槽波示意图;
图7是断层三维空间模型;
图8是不同落差模型提取能量对比图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。本发明的 实施例是为了示例和描述起见而给出的,而并不是无遗漏的或者将本发明限 于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显而易 见的。选择和描述实施例是为了更好说明本发明的原理和实际应用,并且使 本领域的普通技术人员能够理解本发明从而设计适于特定用途的带有各种修 改的各种实施例。
实施例1
一种基于反射槽波的小断层落差探测方法,包括以下步骤:
1、在煤巷单巷道内探测断层区域建立三维煤巷单巷道观测***,观测系 统主要由震源、等间隔分布的多个检波器、以及为震源和检波器提供位置坐 标的三维直角空间坐标系构成,如图1所示。
三维直角空间坐标系沿煤巷建立,X方向指向待探测区域,Y方向垂直 于煤壁,Z方向垂直于煤层顶板。震源采用矿用乳化***在探测起始点激发产 生。三维观测***观测对象包括围岩和煤层,围岩的速度和煤层的体波速度 完全参照现实测试的煤层和岩石速度参数设置,煤层和围岩网格大小一致。
2、分析震源记录,确定反射槽波位置,如图2所示,进行直达波切除处 理。这里的切取直达波指的是包括直达槽波在内的直达波完全切除,获得只 留下反射槽波数据。图3和图4分别是切波前后模拟数据图。
3、对切除直达波后的震源记录进行主振幅提取,获取X分量反射槽波振 幅。振幅提取是将时频-空间域地震记录转换到时-频域,在时-频域定量提取 可表征各地质层位反射波波形总体特征的主特征参数。由信号分析理论可知, 复合信号在时-频域分辨率最高,所以在时-频域中提取的信号特征主参数比在 时间域中提取的精度要高。能量提取是按道数提取反射槽波的能量,提取反 射槽波能量分析时提取每道的最大值。
时频-空间域的能量函数E(t,f,x)=|ε(t,f,x)|
由以上各式可得时频-空间域的主能量函数E(t0,f0,x)=Emax(t0,f0,x),其中 f0和t0分别为地震信号的主频率和群延迟,t0可称为地震信号主能量的延迟时 间。对E0开平方可得到反射波主振幅的值
4、对X分量反射槽波进行滤波,获取滤波后X分量反射槽波能量,具 体包括以下步骤:
⑴获得反射槽波频谱信息;
⑵截取较高能量频段的频宽;
⑶对截取的反射槽波进行频带宽度为50Hz的带通滤波;
(4)对偏移距进行处理得到反射槽波入射角,定义为横坐标;
⑸对滤波后的反射槽波数据进行最大振幅提取。
5、对比分析获得断层的能量变化规律,把观测***的道数换算成反射波 的入射角度,成图进行分析,确定不同落差断层的能量信息,分析断层位置;
6、提取反射槽波的子波类型,确定其频率大小,具体包括以下步骤:
(4)根据频率能量分布确定使用频率大小。
三维模拟子波形态(参照图5)为从实际信号中提取的子波形态,子波的 频率为实际的槽波频谱分析以后得到的最大能量最高的主频,子波类别采用 可以产生love槽波的双向平行的震源。模型选取子波类型为single子波,频 率120Hz,模型边界采用多层pml边界,确保没有边界反射的干扰,有限差 分采用多梯度差分,在煤层中细化差分网格,在模型边界采用大网格差分, 节省模型运算时间。
7、反射槽波示意图如图6所示,根据步骤5和步骤6所得的信息建立不 同落差断层的三维空间模型,如图7所示,得到不同落差断层反射槽波的能 量变化规律,如图8所示。从图8可以看到振幅大小与入射角有着较强的关 联性,而入射角与断层落差存在较强的关联性,从而不同落差断层反射槽波 的能量变化规律通过达到探测小断层落差的目的。
显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实 施例。基于本发明中的实施例,本领域及相关领域的普通技术人员在没有作 出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范 围。
Claims (5)
1.一种基于反射槽波的小断层落差探测方法,其特征在于,包括以下步骤:
步骤1,在煤巷单巷道内探测断层区域建立三维煤巷单巷道观测***,观测***主要由震源、等间隔分布的多个检波器、以及为震源和检波器提供位置坐标的三维直角空间坐标系构成;
步骤2,分析震源记录,确定反射槽波位置,进行直达波切除处理;
步骤3,对切除直达波后的震源记录进行主振幅提取,获取X分量反射槽波振幅;
步骤4,对X分量反射槽波进行滤波,获取滤波后X分量反射槽波能量;
步骤5,对比分析获得断层的能量变化规律,根据数据分析断层位置;
步骤6,提取反射槽波的子波类型,确定其频率大小;
步骤7,根据步骤5和步骤6所得的信息建立不同落差断层的三维空间模型,得到不同落差断层反射槽波的能量变化规律,达到探测小断层落差的目的。
2.根据权利要求1所述的基于反射槽波的小断层落差探测方法,其特征在于,所述步骤1具体包括以下步骤:
步骤A1,在煤巷单巷道内探测断层区域的煤层中间高度等间隔布置n个横波检波器;
步骤A2,采用矿用乳化***在探测起始点激发一次,产生震源;
步骤A3,沿煤巷建立三维直角空间坐标系,X方向指向待探测区域,Y方向垂直于煤壁,Z方向垂直于煤层顶板;
步骤A4,确定震源、n个检波器在三维直角空间坐标中的坐标位置。
4.根据权利要求1或2所述的基于反射槽波的小断层落差探测方法,其特征在于,所述步骤4包括以下步骤:
步骤B1,获得反射槽波频谱信息;
步骤B2,截取较高能量频段的频宽;
步骤B3,对截取的反射槽波进行频带宽度为50Hz的带通滤波;
步骤B4,对偏移距进行处理得到反射槽波入射角,定义为横坐标;
步骤B5,对滤波后的反射槽波数据进行最大振幅提取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010827051.7A CN111965695B (zh) | 2020-08-17 | 2020-08-17 | 一种基于反射槽波的小断层落差探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010827051.7A CN111965695B (zh) | 2020-08-17 | 2020-08-17 | 一种基于反射槽波的小断层落差探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111965695A true CN111965695A (zh) | 2020-11-20 |
CN111965695B CN111965695B (zh) | 2023-04-25 |
Family
ID=73388192
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010827051.7A Active CN111965695B (zh) | 2020-08-17 | 2020-08-17 | 一种基于反射槽波的小断层落差探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111965695B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4653031A (en) * | 1979-10-19 | 1987-03-24 | Coal Industry (Patents) Limited | Mapping faults in a geological seam |
US20080123469A1 (en) * | 2006-11-28 | 2008-05-29 | Magnitude Spas | System and method for seismic pattern recognition |
WO2013152468A1 (zh) * | 2012-04-13 | 2013-10-17 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN106526673A (zh) * | 2016-11-10 | 2017-03-22 | 淮南矿业(集团)有限责任公司 | 一种反射法槽波地震勘探巷道超前探测断层试验方法 |
WO2018129844A1 (zh) * | 2017-01-10 | 2018-07-19 | 中国科学院地质与地球物理研究所 | 一种地震绕射波分离方法和装置 |
CN108761532A (zh) * | 2018-05-08 | 2018-11-06 | 河北煤炭科学研究院 | 煤矿工作面断层预测方法 |
CN110531418A (zh) * | 2019-08-21 | 2019-12-03 | 徐州工程学院 | 一种基于希尔伯特极化成像的断点三维精细定位方法 |
CN110531413A (zh) * | 2019-08-21 | 2019-12-03 | 中国矿业大学 | 一种小断层超前三维可视化建模方法 |
-
2020
- 2020-08-17 CN CN202010827051.7A patent/CN111965695B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4653031A (en) * | 1979-10-19 | 1987-03-24 | Coal Industry (Patents) Limited | Mapping faults in a geological seam |
US20080123469A1 (en) * | 2006-11-28 | 2008-05-29 | Magnitude Spas | System and method for seismic pattern recognition |
WO2013152468A1 (zh) * | 2012-04-13 | 2013-10-17 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN106526673A (zh) * | 2016-11-10 | 2017-03-22 | 淮南矿业(集团)有限责任公司 | 一种反射法槽波地震勘探巷道超前探测断层试验方法 |
WO2018129844A1 (zh) * | 2017-01-10 | 2018-07-19 | 中国科学院地质与地球物理研究所 | 一种地震绕射波分离方法和装置 |
CN108761532A (zh) * | 2018-05-08 | 2018-11-06 | 河北煤炭科学研究院 | 煤矿工作面断层预测方法 |
CN110531418A (zh) * | 2019-08-21 | 2019-12-03 | 徐州工程学院 | 一种基于希尔伯特极化成像的断点三维精细定位方法 |
CN110531413A (zh) * | 2019-08-21 | 2019-12-03 | 中国矿业大学 | 一种小断层超前三维可视化建模方法 |
Non-Patent Citations (2)
Title |
---|
乔勇虎等: "含小断层煤层Rayleigh型槽波波场和频散分析", 《地球物理学报》 * |
李刚等: "煤矿井下工作面内隐伏断层透射槽波探测技术", 《煤田地质与勘探》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111965695B (zh) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10302788B2 (en) | Method and apparatus for identifying collapsed coal column | |
CN104678428B (zh) | 隧道掘进机破岩震源和主动源三维地震联合超前探测*** | |
Schwenk et al. | Surface-wave methods for anomaly detection | |
Senkaya et al. | A semi-automatic approach to identify first arrival time: the cross-correlation technique (CCT) | |
NO20130024A1 (no) | Fremgangsmate for analysering av seismiske data | |
CN105223608A (zh) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 | |
CN107255831A (zh) | 一种叠前频散属性的提取方法 | |
D’hour et al. | Detection of subtle hydromechanical medium changes caused by a small-magnitude earthquake swarm in NE Brazil | |
Girard et al. | Automated processing strategies for ambient seismic data | |
CN108828664B (zh) | 一种多次波识别方法及装置 | |
Rector et al. | Tomographic imaging of surface waves: A case study from the Phoenix Mine, Battle Mountain, Nevada | |
Preiswerk et al. | Effects of geometry on the seismic wavefield of Alpine glaciers | |
Pitarka et al. | Using dense array waveform correlations to build a velocity model with stochastic variability | |
Poiata et al. | Complexity of deep low‐frequency earthquake activity in Shikoku (Japan) imaged from the analysis of continuous seismic data | |
Li et al. | High resolution processing of 3D seismic data for thin coal seam in Guqiao coal mine | |
Noguchi et al. | FDM simulation of an anomalous later phase from the Japan Trench subduction zone earthquakes | |
WO2013057484A2 (en) | Seismic survey data collection and processing | |
CN111965695A (zh) | 一种基于反射槽波的小断层落差探测方法 | |
Alaei | Improved fault imaging by integration of frequency decomposition and fault attributes, example from mid Norwegian sea | |
Tang et al. | Using Green's theorem to satisfy data requirements of multiple removal methods: The impact of acquisition design | |
CN103760599B (zh) | 一种微小断层检测方法及断层检测装置 | |
Nasıf | Processing and joint interpretation of multi-resolution marine seismic datasets | |
JP4069177B2 (ja) | 断層破砕帯判定方法及び判定支援装置 | |
Petrosino et al. | Shallow shear-wave velocity structure of Solfatara volcano (Campi Flegrei, Italy), from inversion of Rayleigh-wave dispersion curves | |
Chen | Application of Radon Transform in in-seam Wave Advanced Detection Technology |
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 |