CN106772575A - 一种基于折射波与面波联合反演剩余煤层厚度的方法 - Google Patents

一种基于折射波与面波联合反演剩余煤层厚度的方法 Download PDF

Info

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
Application number
CN201611063733.5A
Other languages
English (en)
Inventor
胡泽安
张平松
许光泉
黄晖
韩必武
刘亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN201611063733.5A priority Critical patent/CN106772575A/zh
Publication of CN106772575A publication Critical patent/CN106772575A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282Application of seismic models, synthetic seismograms

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),直至两种曲线的残差同时达到残差终止条件,计算停止,此时得到最优解的剩余煤厚模型。
CN201611063733.5A 2016-11-28 2016-11-28 一种基于折射波与面波联合反演剩余煤层厚度的方法 Pending CN106772575A (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 中国石油化工股份有限公司 一种基于表层横波速度反演的纵横波联合静校正方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
傅皓淳: "煤田槽波地震勘探中层析成像技术应用研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (9)

* Cited by examiner, † Cited by third party
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) 一种基于折射波与面波联合反演剩余煤层厚度的方法
CN105589100B (zh) 一种微地震震源位置和速度模型同时反演方法
CN109613616B (zh) 一种浅层岩溶隧道地质综合立体探测方法
CN103472483B (zh) 基于真地表或浮动基准面的速度建模方法
CN105093320B (zh) 针对高速结晶盐壳覆盖区层析静校正初至拾取方法
CN107783187B (zh) 一种将测井速度和地震速度结合建立三维速度场的方法
CN101561512A (zh) 一种井间多尺度sirt层析成像的方法
CN111077572A (zh) 一种基于透射槽波频散曲线反演的煤厚定量预测方法
CN104422963B (zh) 一种变偏移距vsp资料时差校正方法
CN110687602A (zh) 浅层地震多波联合勘探方法
CN104374828A (zh) 一种隐患探测的超声波层析成像方法
CN105301647B (zh) 区分灰质泥岩和砂岩的方法
CN104422955B (zh) 一种利用旅行时变化量进行各向异性参数提取的方法
CN113376695B (zh) 一种适用于煤层底板复杂陷落柱的全波形反演方法
JPH077068B2 (ja) 水平鑿井中での非常に高い分解能の地震探査法
CN102877828A (zh) 一种三维多井联合井地ct成像方法
CN117192615A (zh) 基于透射地震尾波的采煤工作面内隐伏地质构造探测方法
CN102466818B (zh) 一种利用井间地震数据对各向异性介质成像的方法
CN111852467A (zh) 一种砂岩铀矿矿体延伸范围的圈定方法及系统
CN104570104A (zh) 一种基于两步法avf的纵横波地震品质因子提取方法
CN104570091B (zh) 一种获得初至波射线的方法和装置
CN111897029A (zh) 一种岩心-测井交互对比确定地下背斜中裂缝空间展布的方法
CN111812714B (zh) 基于折射纵波与高频槽波的煤层纵横波速度求取方法
CN105676283B (zh) 一种利用斜井直达地震波走时计算地层地震波速度的方法
CN110764140A (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