CN106019374B - 基于反射槽波频散相似度的断层成像方法 - Google Patents

基于反射槽波频散相似度的断层成像方法 Download PDF

Info

Publication number
CN106019374B
CN106019374B CN201610517975.0A CN201610517975A CN106019374B CN 106019374 B CN106019374 B CN 106019374B CN 201610517975 A CN201610517975 A CN 201610517975A CN 106019374 B CN106019374 B CN 106019374B
Authority
CN
China
Prior art keywords
wave
imaging
frequency
point
slot wave
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
Application number
CN201610517975.0A
Other languages
English (en)
Other versions
CN106019374A (zh
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.)
Xian Research Institute Co Ltd of CCTEG
Original Assignee
Xian Research Institute Co Ltd of CCTEG
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 Xian Research Institute Co Ltd of CCTEG filed Critical Xian Research Institute Co Ltd of CCTEG
Priority to CN201610517975.0A priority Critical patent/CN106019374B/zh
Publication of CN106019374A publication Critical patent/CN106019374A/zh
Application granted granted Critical
Publication of CN106019374B publication Critical patent/CN106019374B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/65—Source localisation, e.g. faults, hypocenters or reservoirs

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

基于反射槽波频散相似度的断层成像方法
技术领域
本发明属于煤矿井下物探技术领域，具体涉及一种基于反射槽波频散相似度的断层成像方法。
背景技术
理论的反射槽波成像方法为延时求和偏移法。其思想是依据理论频散曲线对各道信号进行相位补偿，然而在计算出延迟时间后将各道信号延时时刻的振幅叠加起来。实际资料处理中，由于频散曲线的理论值与实际值有一定差异，相位补偿的效果并不理想。因此常用方法是依据中心频率进行带通滤波后求包络，再对包络进行延时叠加。当反射槽波较强能形成单独包络时，采用此方法能够获得一定的成像效果。然而，由于井下噪声较强且直达槽波波列较长，反射槽波波列经常被淹没在噪声或直达槽波波列中，无法形成包络，成像效果较差。已有的反射槽波法对断层成像的效果并不理想，难以满足实际需要。
事实上，频散性是槽波区别于其他弹性波的重要特性之一。已有方法无论是相位补偿还是求包络，都没有利用频散性，甚至将其视为有害信息而加以消除。针对这一问题，我们研究发现反射对槽波频散的影响并不大，发生反射前后槽波在“速度-频率”域内具有很强的相似性。因此可以利用这一性质，首先提取直达槽波的“速度-频率”谱作为模板，在成像过程中，根据成像点的位置将数据转换到“速度-频率”域内以后计算其与直达槽波“速度-频率”谱的相似性。如果成像点是反射点，则相似性就高。通过计算每一个成像点处反射槽波与直达槽波“速度-频率”谱的相似度，就可以对断层成像。
发明内容
有鉴于此，本发明的主要目的在于提供一种基于反射槽波频散相似度的断层成像方法。
为达到上述目的，本发明的技术方案是这样实现的：
本发明实施例提供一种基于反射槽波频散相似度的断层成像方法，通过以下步骤实现：
步骤一：时频分解，用广义S变换求出各道地震数据的时频谱
其中sn(t)为第n道地震数据的时域形式，Sn(τ,f)为其时频谱，n∈[1,N]，N为总道数；步骤二：求直达槽波频散曲线，依据各道的炮点位置(xsn,ysn)与检波点位置(xrn,yrn)，求出直达槽波的传播距离Ln
再由Ln对Sn(τ,f)进行坐标转换，可求出各道直达槽波的频散曲线，叠加后获得总体数据中直达槽波的频散曲线C(v,f)
步骤三：消除直达槽波，根据Ln将C(v,f)从“速度-频率”域转换回“时间-频率”域后，与Sn(τ,f)相减，就可以得到消除了直达槽波后的各道时频谱S’n(τ,f)
步骤四：成像，首先选定一个成像区域，区域的一个点P，坐标为(xp,yp)，假设P点是一个反射点，则槽波从第n道数据的炮点(xsn,ysn)出发，在P点反射后，被(xrn,yrn)处的检波器接收到，则传播路径的总长度为
如果P点是反射点，则根据Lpn求出的频散曲线Cpn(v,f)应该与C(v,f)一致，因此，将所有道Cpn(v,f)与C(v,f)的相关系数的和Rp作为该成像点的值，
当该点确实为反射点时，Rp较大，而当它不是反射点时，Rp较小，对成像区域内所有成像点进行上述处理，得到整个区域的成像结果。
上述方案中，所述步骤一具体为：
1)令n＝1，读取第n道地震数据sn(t)；
2)用广义S变换求sn(t)对应的时频谱Sn(τ,f)；
3)保存Sn(τ,f)；
4)令n＝n+1，处理下一道直到所有道处理完毕。
上述方案中，所述步骤二具体为：
1)令n＝1，读取第n道时频谱Sn(τ,f)；
2)从第n道道头中读取炮点位置(xsn,ysn)与检波点位置(xrn,yrn)；
3)计算直达槽波的传播距离Ln；
4)时频谱Sn(τ,f)转换为速度-频率谱Cn(v,f)；
5)令n＝n+1，处理下一道直到所有道处理完毕；
6)将所有道的Cn(v,f)叠加得到直达槽波的频散曲线C(v,f)。
上述方案中，所述步骤三具体为：
1)读取直达槽波频散曲线C(v,f)；
2)令n＝1，读取第n道对应的直达槽波传播距离Ln；
3)读取第n道时频谱Sn(τ,f)；由Ln将C(v,f)转换到“时间-频率”域后与Sn(τ,f)相减，得到消除了直达槽波后时频谱S’n(τ,f)；
4)保存S’n(τ,f)；
5)令n＝n+1，处理下一道直到所有道处理完毕。
上述方案中，所述步骤四具体为：
1)选定成像范围；
2)令x＝1，y＝1，从第1个点开始对成像范围内的每一个点逐个处理；
3)令n＝1，从道头中读取炮点位置(xsn,ysn)与检波点位置(xrn,yrn)；
4)计算槽波从炮点到反射点，再到检波点的传播距离Lpn；
5)读取S’n(τ,f)；
6)求当以(x,y)为坐标的点P为反射点时，第n道反射槽波的频散曲线Cpn(v,f)；
7)求Cpn(v,f)与C(v,f)的相关系数Rpn；
8)令n＝n+1，处理下一道直到所有道处理完毕；
9)将所有道的Rpn叠加后的值Rp作为P点对应的成像值；
10)令x＝x+1或y＝y+1，求下一个点的Rp直到成像区域内所有点成像完毕；
11)最终显示成像结果。
与现有技术相比，本发明的有益效果：
本发明充分利用了槽波的频散性质，依据频散性减少了其他噪声对成像的影响，减少了画弧现象，大大提高了断层成像的精度。
附图说明
图1为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的总流程图；
图2为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的时频分解流程图；
图3为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的求直达槽波频散曲线流程图；
图4为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的直达槽波压制流程图；
图5为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的反射槽波成像流程图；
图6为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的断层位置与观测***布置图；
图7为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的第20炮模型数据；
图8为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的直达槽波频散曲线；
图9为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的按直达波路径求出第30道的速度-频率谱；
图10为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的压制直达槽波后结果；
图11为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的按反射路径求出的速度-频率谱；
图12为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的已有方法的断层成像结果；
图13为本发明实施例提供一种基于反射槽波频散相似度的断层成像方法的本发明的断层成像结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明实施例提供一种基于反射槽波频散相似度的断层成像方法，通过以下步骤实现：
步骤一：用广义S变换求出各道地震数据的时频谱
其中sn(t)为第n道地震数据的时域形式，Sn(τ,f)为其时频谱。n∈[1,N]，N为总道数；
具体地，如图2所示，令n＝1，读取第n道地震数据sn(t)；用广义S变换求sn(t)对应的时频谱Sn(τ,f)；保存Sn(τ,f)；令n＝n+1，处理下一道直到所有道处理完毕。
步骤二：依据各道的炮点位置(xsn,ysn)与检波点位置(xrn,yrn)，求出直达槽波的传播距离Ln
具体地，如图3所示，令n＝1，读取第n道时频谱Sn(τ,f)；从第n道道头中读取炮点位置(xsn,ysn)与检波点位置(xrn,yrn)；计算直达槽波的传播距离Ln；时频谱Sn(τ,f)转换为速度-频率谱Cn(v,f)；令n＝n+1，处理下一道直到所有道处理完毕；将所有道的Cn(v,f)叠加得到直达槽波的频散曲线C(v,f)。
步骤三：由Ln对Sn(τ,f)进行坐标转换，可求出各道直达槽波的频散曲线，叠加后获得总体数据中直达槽波的频散曲线C(v,f)
再根据Ln将C(v,f)从“速度-频率”域转换回“时间-频率”域后，与Sn(τ,f)相减，就可以得到消除了直达槽波后的各道时频谱S’n(τ,f)
具体地，如图4所示，读取直达槽波频散曲线C(v,f)；令n＝1，读取第n道对应的直达槽波传播距离Ln；读取第n道时频谱Sn(τ,f)；由Ln将C(v,f)转换到“时间-频率”域后与Sn(τ,f)相减，得到消除了直达槽波后时频谱S’n(τ,f)；保存S’n(τ,f)；令n＝n+1，处理下一道直到所有道处理完毕。步骤四：成像，首先选定一个成
像区域，区域的一个点P，坐标为(xp,yp)。
假设P点是一个反射点，则槽波从第n道数据的炮点(xsn,ysn)出发，在P点反射后，被(xrn,yrn)处的检波器接收到，则传播路径的总长度为
如果P点是反射点，则根据Lpn求出的频散曲线Cpn(v,f)应该与C(v,f)一致，因此，将所有道Cpn(v,f)与C(v,f)的相关系数的和Rp作为该成像点的值，
当该点确实为反射点时，Rp较大，而当它不是反射点时，Rp较小，对成像区域内所有成像点进行上述处理，得到整个区域的成像结果。
具体地，如图5所示，读取直达槽波频散曲线C(v,f)；选定成像范围；令x＝1，y＝1，从第1个点开始对成像范围内的每一个点逐个处理；令n＝1，从道头中读取炮点位置(xsn,ysn)与检波点位置(xrn,yrn)；计算槽波从炮点到反射点，再到检波点的传播距离Lpn；读取S’n(τ,f)；求当以(x,y)为坐标的点P为反射点时，第n道反射槽波的频散曲线Cpn(v,f)；求Cpn(v,f)与C(v,f)的相关系数Rpn；令n＝n+1，处理下一道直到所有道处理完毕；将所有道的Rpn叠加后的值Rp作为P点对应的成像值；令x＝x+1或y＝y+1，求下一个点的Rp直到成像区域内所有点成像完毕；最终显示成像结果。
下面以理论模型为例说明本发明效果：模型由顶板、煤层、底板三层
组成。围岩纵波速度3800m/s，横波速度2000m/s，密度2.4g/cm3；煤层厚度5m，纵波速度1800m/s，横波速度1100m/s，密度1.4g/cm3。模型中包含一条斜向断层，落差5m。断层位置与观测***布置如图6所示。图7为根据三维弹性波正演模拟方法得到的第20炮正演模拟的结果，从中可以分辨出直达槽波与断层的反射槽波。图8由所有炮叠加得到的直达槽波“速度-频率”谱，其中黑色线条为根据物性参数计算出的理论频散曲线。图9是第20炮的第30道数据时频谱转换到“速度-频率谱”的结果，压缩其中的直达槽波，结果如图10所示。根据反射的几何关系，第20炮第20道数据的反射槽波应来自于断层上的点(x＝429,y＝137)。如果将这一点作为反射点，求第20炮第30道数据的“速度-频率谱”，结果如图11所示。可以看出反射槽波与直达槽波的频散曲线具有很强的相似性。本发明正是利用这一相似性对断层成像的。图12是用已有方法成像的结果。图13是用本发明方法成像的结果，可以看出与已有方法相比，新方法的成像结果减少了画弧效应，具有更高的分辨率，能够更准确反映断层的发育情况。断层的像在距离模型右边界20m范围内有能量减弱的现象，这是因为能够在这段范围形成有效反射的炮点已经超出了模型范围，所以这段断层无法成像是正常现象。
以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。

Claims (5)

1.一种基于反射槽波频散相似度的断层成像方法，其特征在于：通过以下步骤实现：
步骤一：时频分解，用广义S变换求出各道地震数据的时频谱
其中sn(t)为第n道地震数据的时域形式，Sn(τ,f)为其时频谱，n∈[1,N]，N为总道数；
步骤二：求直达槽波频散曲线，依据各道的炮点位置(xsn,ysn)与检波点位置(xrn,yrn)，求出直达槽波的传播距离Ln
再由Ln对Sn(τ,f)进行坐标转换，可求出各道直达槽波的频散曲线，叠加后获得总体数据中直达槽波的频散曲线C(v,f)
步骤三：消除直达槽波，根据Ln将C(v,f)从“速度-频率”域转换回“时间-频率”域后，与Sn(τ,f)相减，就可以得到消除了直达槽波后的各道时频谱S’n(τ,f)
步骤四：成像，首先选定一个成像区域，区域的一个点P，坐标为(xp,yp)，假设P点是一个反射点，则槽波从第n道数据的炮点(xsn,ysn)出发，在P点反射后，被(xrn,yrn)处的检波器接收到，则传播路径的总长度为
如果P点是反射点，则根据Lpn求出的频散曲线Cpn(v,f)应该与C(v,f)一致，因此，将所有道Cpn(v,f)与C(v,f)的相关系数的和Rp作为该成像点的值，
当该点确实为反射点时，Rp较大，而当它不是反射点时，Rp较小，对成像区域内所有成像点进行上述处理，得到整个区域的成像结果。
2.根据权利要求1所述的基于反射槽波频散相似度的断层成像方法，其特征在于，所述步骤一具体为：
1)令n＝1，读取第n道地震数据sn(t)；
2)用广义S变换求sn(t)对应的时频谱Sn(τ,f)；
3)保存Sn(τ,f)；
4)令n＝n+1，处理下一道直到所有道处理完毕。
3.根据权利要求1或2所述的基于反射槽波频散相似度的断层成像方法，其特征在于，所述步骤二具体为：
1)令n＝1，读取第n道时频谱Sn(τ,f)；
2)从第n道道头中读取炮点位置(xsn,ysn)与检波点位置(xrn,yrn)；
3)计算直达槽波的传播距离Ln；
4)时频谱Sn(τ,f)转换为速度-频率谱Cn(v,f)；
5)令n＝n+1，处理下一道直到所有道处理完毕；
6)将所有道的Cn(v,f)叠加得到直达槽波的频散曲线C(v,f)。
4.根据权利要求3所述的基于反射槽波频散相似度的断层成像方法，其特征在于，所述步骤三具体为：
1)读取直达槽波频散曲线C(v,f)；
2)令n＝1，读取第n道对应的直达槽波传播距离Ln；
3)读取第n道时频谱Sn(τ,f)；
4)由Ln将C(v,f)转换到“时间-频率”域后与Sn(τ,f)相减，得到消除了直达槽波后时频谱S’n(τ,f)；
5)保存S’n(τ,f)；
6)令n＝n+1，处理下一道直到所有道处理完毕。
5.根据权利要求4所述的基于反射槽波频散相似度的断层成像方法，其特征在于，所述步骤四具体为：
1)选定成像范围；
2)令x＝1，y＝1，从第1个点开始对成像范围内的每一个点逐个处理；
3)令n＝1，从道头中读取炮点位置(xsn,ysn)与检波点位置(xrn,yrn)；
4)计算槽波从炮点到反射点，再到检波点的传播距离Lpn；
5)读取S’n(τ,f)；
6)求当以(x,y)为坐标的点P为反射点时，第n道反射槽波的频散曲线Cpn(v,f)；
7)求Cpn(v,f)与C(v,f)的相关系数Rpn；
8)令n＝n+1，处理下一道直到所有道处理完毕；
9)将所有道的Rpn叠加后的值Rp作为P点对应的成像值；
10)令x＝x+1或y＝y+1，求下一个点的Rp直到成像区域内所有点成像完毕；
11)最终显示成像结果。
CN201610517975.0A 2016-07-04 2016-07-04 基于反射槽波频散相似度的断层成像方法 Active CN106019374B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610517975.0A CN106019374B (zh) 2016-07-04 2016-07-04 基于反射槽波频散相似度的断层成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610517975.0A CN106019374B (zh) 2016-07-04 2016-07-04 基于反射槽波频散相似度的断层成像方法

Publications (2)

Publication Number Publication Date
CN106019374A CN106019374A (zh) 2016-10-12
CN106019374B true CN106019374B (zh) 2018-06-26

Family

ID=57106812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610517975.0A Active CN106019374B (zh) 2016-07-04 2016-07-04 基于反射槽波频散相似度的断层成像方法

Country Status (1)

Country Link
CN (1) CN106019374B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107229597A (zh) * 2017-05-31 2017-10-03 成都理工大学 同步挤压广义s变换信号时频分解与重构方法
CN107831530B (zh) * 2017-10-25 2019-03-26 山东科技大学 厚煤层沿底板或顶板掘进煤巷反射槽波超前探测方法
CN109581500B (zh) * 2018-12-18 2020-06-30 东华理工大学 一种反射地震记录频变速度分析方法
CN109765606B (zh) * 2019-01-28 2020-08-04 阳泉煤业（集团）有限责任公司 一种基于反射槽波的回采工作面隐蔽断层性质的探测方法
CN113325467B (zh) * 2021-06-08 2023-10-24 中煤科工集团西安研究院有限公司 一种基于槽波频散特征的微地震震源定位方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102788991B (zh) * 2012-07-16 2015-05-06 中煤科工集团西安研究院有限公司 基于z分量的透射槽波高密度快速探测方法
CN102788995A (zh) * 2012-08-02 2012-11-21 中煤科工集团西安研究院 以切割震动为地震信号的煤矿工作面探测方法
US20140153365A1 (en) * 2012-11-30 2014-06-05 Chevron U.S.A. Inc. System and method for producing local images of subsurface targets
CN105158794B (zh) * 2015-08-21 2017-07-14 山东科技大学 一种基于四通道冗余滤波器组的槽波弱信号提取方法
CN105549087B (zh) * 2015-12-10 2018-03-06 北京中矿大地地球探测工程技术有限公司 一种煤矿井下槽波地震勘探的走时和振幅联合反演方法

Also Published As

Publication number Publication date
CN106019374A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
CN106019374B (zh) 基于反射槽波频散相似度的断层成像方法
CN106932824B (zh) 陆地地震勘探资料的降维自适应层间多次波压制方法
CN102939546B (zh) 用于在地震处理中的局部属性匹配的***和方法
CN102636811B (zh) 一种海上二维地震资料中多次波的消除方法
CN101957455B (zh) 三维保幅叠前时间偏移方法
CN104391327B (zh) 一种海上斜井井间地震叠前逆时深度偏移成像方法
CN105093292B (zh) 一种地震成像的数据处理方法和装置
Nakamura et al. FDM simulation of seismic‐wave propagation for an aftershock of the 2009 Suruga Bay earthquake: effects of ocean‐bottom topography and seawater layer
CN104483705B (zh) 一种三维剩余静校正方法
MX2012008214A (es) Metodo para procesar datos sismicos marinos.
CN103645497B (zh) 一种基于出射角的多分量波场分离方法
CN104820242B (zh) 一种面向叠前反演的道集振幅分频补偿方法
CN107817523B (zh) 绕射波偏移速度的分析方法及装置
CN102162859A (zh) 一种斜井井间地震波场的成像方法
CN104570116A (zh) 基于地质标志层的时差分析校正方法
CN104280777A (zh) 一种压制陆上地震资料多次波干扰的方法
CN110018517A (zh) 一种多尺度地面微地震逆时干涉定位方法
CN105652322A (zh) 多分量地震数据的t-f-k域极化滤波方法
CN102053260B (zh) 获得地震纵波的方位速度的方法及处理地震数据的方法
CN101609167B (zh) 基于起伏地表的井间地震波动方程叠前深度偏移成像方法
CN103018777A (zh) 一种能够消除逆时偏移噪声的成像方法及装置
CN103149588A (zh) 一种利用井震标定计算vti各向异性参数的方法及***
CN110118993A (zh) 绕射波成像方法及装置
CN106950600A (zh) 一种近地表散射面波的去除方法
CN104216018A (zh) 一种基于小道距层析反演的表层结构调查方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant