CN106772575A - 一种基于折射波与面波联合反演剩余煤层厚度的方法 - Google Patents
一种基于折射波与面波联合反演剩余煤层厚度的方法 Download PDFInfo
- Publication number
- CN106772575A CN106772575A CN201611063733.5A CN201611063733A CN106772575A CN 106772575 A CN106772575 A CN 106772575A CN 201611063733 A CN201611063733 A CN 201611063733A CN 106772575 A CN106772575 A CN 106772575A
- Authority
- CN
- China
- Prior art keywords
- wave
- coal seam
- velocity
- model
- thickness
- 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.)
- Pending
Links
- 239000003245 coal Substances 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000006185 dispersion Substances 0.000 claims abstract description 13
- 239000011435 rock Substances 0.000 claims abstract description 9
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 9
- 238000004458 analytical method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000010454 slate Substances 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005641 tunneling Effects 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/282—Application of seismic models, synthetic seismograms
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
技术领域
本发明涉及煤层探测方法领域,具体是一种基于折射波与面波联合反演剩余煤层厚度的方法。
背景技术
在矿井生产中,了解煤层结构和厚度的变化情况是合理确定分层开采的必要条件,也是防止煤炭资源损失的关键。长期以来,煤层开采剩余厚度的探测工作一直采用钻探测的方法,不仅劳动强度大,效率低,而且往往不能及时准确地提供煤厚资料。与之相比地震波法具有操作简单、迅速,解释精度高的特点,可以更加有效地进行剩余煤厚的探测。目前巷道探测剩余煤厚的物探方法主要有相遇折射波法,反射波法等。然而,从国内外资料分析可见,由于通常剩余煤厚很薄;且测试场地破坏松动,造成勘探的地震波频率很低,分辨率不高。因此为提高煤层开采利用率,急迫需要一种快捷有效的基于地震信号处理的,能够探测剩余煤厚的方法。
煤层与底板岩石波阻抗差异明显,根据地震波的传播原理,地震波传播过程中一部分能量以折射波的路径传播,另一部分以面波的形态在剩余的煤层中传播。折射波和面波均可单独勘探剩余煤厚,因为可利用联合反演的方法来提高探测剩余煤厚的精度与准确性。
发明内容
本发明的目的是提供一种基于折射波与面波联合反演剩余煤层厚度的方法,以解决现有技术剩余煤层厚度探测方法精度和准确性差的问题。
为了达到上述目的,本发明所采用的技术方案为:
一种基于折射波与面波联合反演剩余煤层厚度的方法,其特征在于:包括以下步骤:
(1)、在巷道底板位置布置多个炮点,炮点从左向右布置,在炮点正前方处布置22个X、Z两分量检波点,接收排列随炮点移动,通过巷道内布置的地震仪获取地震记录,根据地震波传播时间先后顺序,地震记录中可清楚确定直达波、折射波和面波;
(2)、拾取首至波旅行时形成的时距曲线L1,根据折射波法得到初始的两层速度模型,即可以得到出表层煤层纵波速度V1'和下层底板岩石纵波速度V2'以及剩余煤层的厚度H1;
(3)、对数据中的面波信号进行频散处理,得到实测面波信号的频散曲线L2,根据面波法得到初始的两层速度模型,又根据横波速度与纵波速度的经验公式,可得表层煤层纵波速度V1”和下层底板岩石纵波速度V2”以及剩余煤层的厚度H2;
(4)、根据折射波和面波反演模型参数,求取均值,作为新的初始反演模型。模型参数为以此模型计算理论首至波曲线和理论面波频散曲线
(5)、计算L1和曲线L2之间的残差,并设置残差终止条件(一般设置10%以内)。只要其中之一不满足残差终止条件,则调整模型的表层煤层纵波速度参数下层底板岩石纵波速度和剩余煤层的厚度H*;
(6)、调整模型后,重复步骤(4)、(5),直至两种曲线的残差同时达到残差终止条件,计算停止,此时得到最优解的剩余煤厚模型。
本发明提供了一种基于折射波与面波联合反演剩余煤层厚度的方法,利用折射波和面波探测数据,进行重构地质模型。通过迭代的方式,构建一个同时满足折射波和面波特征的地质模型,达到探测剩余煤厚的目标。
附图说明
图1为本发明观测系统布置图。
图2为折射波与面波信号图。
图3为折射波反演流程图,其中:
图3a为直达波和折射波到时拾取,图3b为直达波和折射波的拟合时距曲线;
图3c为折射波反演结果。
图4为面波反演流程图,其中:
图4a为面波信号频散分析与曲线拾取,图4b为面波频散曲线反演结果,
图4c为面波反演结果。
图5为折射波和面波联合反演流程图。
具体实施方式
一种基于折射波与面波联合反演剩余煤层厚度的方法,包括以下步骤:
(1)、如图1所示,在巷道底板位置布置多个炮点,炮点从左向右布置,在炮点正前方(图中向右为掘进方向)处布置22个X、Z两分量检波点,接收排列随炮点移动,通过巷道内布置的地震仪获取地震记录;
以某一相同测点的联合反演为例进行说明。实测的共炮点数据,根据地震波传播时间先后顺序,地震记录中可清楚确定直达波、折射波和面波,如图2所示。
(2)、拾取首至波旅行时(如图3a所示)形成的时距曲线L1(如图3b所示),根据折射波法得到初始的两层速度模型(如图3c所示),即可以得到出表层煤层纵波速度V1'和下层底板岩石纵波速度V2'以及剩余煤层的厚度H1;
(3)、对数据中的面波信号进行频散处理,得到实测面波信号的频散曲线L1,如图4所示,根据面波法得到初始的两层速度模型,又根据横波速度与纵波速度的经验公式,可得表层煤层纵波速度V1”和下层底板岩石纵波速度V2”以及剩余煤层的厚度H2;
(4)、根据折射波和面波反演模型参数,求取均值,作为新的初始反演模型。模型参数为以此模型计算理论首至波曲线和理论面波频散曲线
(5)、计算L1和曲线L2之间的残差,并设置残差终止条件(一般设置10%以内)。只要其中之一不满足残差终止条件,则调整模型的表层煤层纵波速度参数下层底板岩石纵波速度和剩余煤层的厚度H*;
(6)、调整模型后,重复步骤(4)、(5),直至两种曲线的残差同时达到残差终止条件,计算停止,此时得到最优解的剩余煤厚。详见图5,其中图5中剩余煤厚为5m。
Claims (1)
1.一种基于折射波与面波联合反演剩余煤层厚度的方法,其特征在于:包括以下步骤:
(1)、在巷道底板位置布置多个炮点,炮点从左向右布置,在炮点正前方处布置22个X、Z两分量检波点,接收排列随炮点移动,通过巷道内布置的地震仪获取地震记录,根据地震波传播时间先后顺序,地震记录中可清楚确定直达波、折射波和面波;
(2)、拾取首至波旅行时形成的时距曲线L1,根据折射波法得到初始的两层速度模型,即可以得到出表层煤层纵波速度V1'和下层底板岩石纵波速度V2'以及剩余煤层的厚度H1;
(3)、对数据中的面波信号进行频散处理,得到实测面波信号的频散曲线L2,根据面波法得到初始的两层速度模型,又根据横波速度与纵波速度的经验公式,可得表层煤层纵波速度V1”和下层底板岩石纵波速度V2”以及剩余煤层的厚度H2;
(4)、根据折射波和面波反演模型参数,求取均值,作为新的初始反演模型。模型参数为以此模型计算理论首至波曲线和理论面波频散曲线
(5)、计算L1和曲线L2之间的残差,并设置残差终止条件(一般设置10%以内)。只要其中之一不满足残差终止条件,则调整模型的表层煤层纵波速度参数V1 *、下层底板岩石纵波速度和剩余煤层的厚度H*;
(6)、调整模型后,重复步骤(4)、(5),直至两种曲线的残差同时达到残差终止条件,计算停止,此时得到最优解的剩余煤厚模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611063733.5A CN106772575A (zh) | 2016-11-28 | 2016-11-28 | 一种基于折射波与面波联合反演剩余煤层厚度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611063733.5A CN106772575A (zh) | 2016-11-28 | 2016-11-28 | 一种基于折射波与面波联合反演剩余煤层厚度的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106772575A true CN106772575A (zh) | 2017-05-31 |
Family
ID=58901885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611063733.5A Pending CN106772575A (zh) | 2016-11-28 | 2016-11-28 | 一种基于折射波与面波联合反演剩余煤层厚度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772575A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108426545A (zh) * | 2018-03-22 | 2018-08-21 | 天津大学 | 一种应用超声表面波无损检测薄膜厚度的方法 |
CN110687602A (zh) * | 2019-10-31 | 2020-01-14 | 山东电力工程咨询院有限公司 | 浅层地震多波联合勘探方法 |
WO2020029015A1 (zh) * | 2018-08-06 | 2020-02-13 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
CN112731558A (zh) * | 2020-12-16 | 2021-04-30 | 中国科学技术大学 | 一种地震面波与接收函数的联合反演方法及装置 |
CN115932945A (zh) * | 2022-10-14 | 2023-04-07 | 扎赉诺尔煤业有限责任公司 | 掘进巷道剩余煤厚地震多波多分量的探测方法 |
CN116660974A (zh) * | 2023-04-11 | 2023-08-29 | 中国地震局地球物理研究所 | 一种基于结构耦合约束的体波和面波三维联合反演方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
CN102359924A (zh) * | 2011-09-19 | 2012-02-22 | 中国地质大学(北京) | 一种基于多波地震数据的煤岩强度的检测方法 |
CN104730579A (zh) * | 2013-12-18 | 2015-06-24 | 中国石油化工股份有限公司 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
-
2016
- 2016-11-28 CN CN201611063733.5A patent/CN106772575A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
CN102359924A (zh) * | 2011-09-19 | 2012-02-22 | 中国地质大学(北京) | 一种基于多波地震数据的煤岩强度的检测方法 |
CN104730579A (zh) * | 2013-12-18 | 2015-06-24 | 中国石油化工股份有限公司 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
Non-Patent Citations (1)
Title |
---|
傅皓淳: "煤田槽波地震勘探中层析成像技术应用研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108426545A (zh) * | 2018-03-22 | 2018-08-21 | 天津大学 | 一种应用超声表面波无损检测薄膜厚度的方法 |
WO2020029015A1 (zh) * | 2018-08-06 | 2020-02-13 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
US11137511B2 (en) | 2018-08-06 | 2021-10-05 | Southern University Of Science And Technology | Active source surface wave prospecting method, surface wave exploration device and computer-readable storage medium |
CN110687602A (zh) * | 2019-10-31 | 2020-01-14 | 山东电力工程咨询院有限公司 | 浅层地震多波联合勘探方法 |
CN112731558A (zh) * | 2020-12-16 | 2021-04-30 | 中国科学技术大学 | 一种地震面波与接收函数的联合反演方法及装置 |
CN115932945A (zh) * | 2022-10-14 | 2023-04-07 | 扎赉诺尔煤业有限责任公司 | 掘进巷道剩余煤厚地震多波多分量的探测方法 |
CN115932945B (zh) * | 2022-10-14 | 2024-04-02 | 扎赉诺尔煤业有限责任公司 | 掘进巷道剩余煤厚地震多波多分量的探测方法 |
CN116660974A (zh) * | 2023-04-11 | 2023-08-29 | 中国地震局地球物理研究所 | 一种基于结构耦合约束的体波和面波三维联合反演方法 |
CN116660974B (zh) * | 2023-04-11 | 2024-02-23 | 中国地震局地球物理研究所 | 一种基于结构耦合约束的体波和面波三维联合反演方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772575A (zh) | 一种基于折射波与面波联合反演剩余煤层厚度的方法 | |
CN109613616B (zh) | 一种浅层岩溶隧道地质综合立体探测方法 | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
CN103472483B (zh) | 基于真地表或浮动基准面的速度建模方法 | |
CN105093320B (zh) | 针对高速结晶盐壳覆盖区层析静校正初至拾取方法 | |
CN111077572A (zh) | 一种基于透射槽波频散曲线反演的煤厚定量预测方法 | |
CN101561512A (zh) | 一种井间多尺度sirt层析成像的方法 | |
CN110687602A (zh) | 浅层地震多波联合勘探方法 | |
CN109884710A (zh) | 针对激发井深设计的微测井层析成像方法 | |
CN105301647B (zh) | 区分灰质泥岩和砂岩的方法 | |
CN102877828A (zh) | 一种三维多井联合井地ct成像方法 | |
CN107942372A (zh) | 井间及井地联合地震波ct成像方法及装置 | |
CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
CN113376695B (zh) | 一种适用于煤层底板复杂陷落柱的全波形反演方法 | |
JPH077068B2 (ja) | 水平鑿井中での非常に高い分解能の地震探査法 | |
CN117192615A (zh) | 基于透射地震尾波的采煤工作面内隐伏地质构造探测方法 | |
CN102466818B (zh) | 一种利用井间地震数据对各向异性介质成像的方法 | |
CN116626757A (zh) | 一种隧道超前地质预报方法 | |
CN104570104A (zh) | 一种基于两步法avf的纵横波地震品质因子提取方法 | |
Gu et al. | Investigation of fractures using seismic computerized crosshole tomography | |
CN110764136B (zh) | 各向异性纵横波走时线性组合与非线性组合联合定位方法 | |
CN111897029A (zh) | 一种岩心-测井交互对比确定地下背斜中裂缝空间展布的方法 | |
CN110780345A (zh) | 一种隧道超前地震探测地震数据的三维速度分析方法 | |
CN111812714B (zh) | 基于折射纵波与高频槽波的煤层纵横波速度求取方法 | |
CN105676283B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |
|
RJ01 | Rejection of invention patent application after publication |