CN112346120A - 一种槽波频散曲线的变异度计算方法及异常体的识别方法 - Google Patents

一种槽波频散曲线的变异度计算方法及异常体的识别方法 Download PDF

Info

Publication number
CN112346120A
CN112346120A CN202011202939.8A CN202011202939A CN112346120A CN 112346120 A CN112346120 A CN 112346120A CN 202011202939 A CN202011202939 A CN 202011202939A CN 112346120 A CN112346120 A CN 112346120A
Authority
CN
China
Prior art keywords
wave
dispersion curve
frequency
calculating
groove
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
Application number
CN202011202939.8A
Other languages
English (en)
Other versions
CN112346120B (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.)
Beijing Sikaiweike Physical Geography Information Technology Co ltd
Institute of Geographic Sciences and Natural Resources of CAS
Original Assignee
Beijing Sikaiweike Physical Geography Information Technology Co ltd
Institute of Geographic Sciences and Natural Resources of CAS
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 Beijing Sikaiweike Physical Geography Information Technology Co ltd, Institute of Geographic Sciences and Natural Resources of CAS filed Critical Beijing Sikaiweike Physical Geography Information Technology Co ltd
Priority to CN202011202939.8A priority Critical patent/CN112346120B/zh
Publication of CN112346120A publication Critical patent/CN112346120A/zh
Application granted granted Critical
Publication of CN112346120B publication Critical patent/CN112346120B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

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所示，从而不能准确拾取旅行时，从而使得通过拾取旅行时层析成像得到的煤层内的异常体信息不准确。
申请内容
本申请实施例通过提供一种槽波频散曲线的变异度计算方法及异常体的识别方法，用以解决现有技术中当槽波射线经过变异体时，槽波频散曲线连续性变差，从而不能准确拾取旅行时，从而使得通过拾取旅行时层析成像得到的煤层内的异常体信息不准确的问题。
为了解决上述问题，第一方面，本发明实施例提供了一种槽波频散曲线的变异度计算方法，包括：获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据；根据地震数据计算槽波实际频散曲线及槽波频谱曲线；根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线；根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度。
可选地，根据地震数据计算实际槽波频散曲线及槽波频谱曲线，包括：将地震数据进行窄带滤波，得到槽波实际频散曲线；将地震数据进行傅里叶变换，得到槽波频谱曲线。
可选地，根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度，包括：根据槽波实际频散曲线确定槽波实际频散曲线中各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，第一频率小于第二频率；根据槽波频谱曲线中预设振幅值对应的最大频率确定第三频率；根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点的加权因子；根据槽波实际频散曲线及槽波理论频散曲线确定槽波实际频散曲线中各断点对应的第二频率对应的槽波理论速度值；根据各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，第三频率，各断点的加权因子及各断点对应的第二频率对应的槽波理论速度值计算槽波实际频散曲线的变异度。
可选地，根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点的加权因子，包括：根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点对应的频谱第一振幅值及频谱第二振幅值；根据各断点对应的频谱第一振幅值及频谱第二振幅值确定槽波实际频散曲线中各断点的加权因子。
第二方面，本发明实施例提供了一种异常体的识别方法，通过如第一方面或第一方面任意实施方式中的槽波频散曲线的变异度计算方法计算槽波频散曲线的变异度，方法包括：获取待测区域位置数据、各槽波射线路径数据及各地震数据；根据各地震数据计算各槽波频散曲线的变异度；根据待测区域位置数据对待测区域进行网格划分；根据网格划分结果、各槽波射线路径数据计算各网格区域中各槽波的路径系数；根据各槽波频散曲线的变异度及各网格区域中各槽波的路径系数计算各网格区域的槽波异度值；根据各网格区域的槽波异度值对待测区域的异常体进行识别。
可选地，根据网格划分结果、各槽波射线路径数据计算各网格区域中各槽波的路径系数，包括：根据网格划分结果、各槽波射线路径数据确定各网格区域中各射线路径经过的长度及各网格的对角线长度；根据各网格区域中各射线路径经过的长度及各网格的对角线长度计算各网格区域中各槽波的路径系数。
可选地，根据各槽波频散曲线的变异度及各网格区域中各槽波的路径系数计算各网格区域的槽波异度值，包括：根据各槽波频散曲线的变异度与各网格区域中各槽波的路径系数计算各网格区域中各槽波的异度值；将各网格区域中各槽波的异度值进行累加得到各网格区域的槽波异度值。
可选地，根据各网格区域的槽波异度值对待测区域的异常体进行识别，包括：根据各网格区域的槽波异度值及槽波异度值预设等级对应的预设阈值确定各网格区域的槽波异度值所属的等级；根据各网格区域的槽波异度值所属的等级对待测区域的异常体进行识别。
第三方面，本发明实施例提供了一种计算机，包括：至少一个处理器；以及与至少一个处理器通信连接的存储器；其中，存储器存储有可被一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器执行如第一方面或第一方面任意实施方式中的槽波频散曲线的变异度计算方法或如第二方面或第二方面任意实施方式中的异常体的识别方法。
第四方面，本发明实施例提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机指令，计算机指令用于使计算机执行如第一方面或第一方面任意实施方式中的槽波频散曲线的变异度计算方法或如第二方面或第二方面任意实施方式中的异常体的识别方法。
本发明实施例提供的槽波频散曲线的变异度计算方法及异常体的识别方法，通过获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据，根据地震数据计算槽波实际频散曲线及槽波频谱曲线，根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线，根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度，从而可以对频散曲线的变异度进行定量计算，从而可以根据待测区域内各频散曲线的变异度对待测区域的异常体进行成像，从而识别待测区域的异常体，由于可以准确地定量待测区域内各频散曲线的变异度，从而根据待测区域内各频散曲线的变异度成像后得到的待测区域的异常体信息更加准确，从而可以解决现有技术中当槽波射线经过变异体时，槽波频散曲线连续性变差，从而不能准确拾取旅行时，使得通过拾取旅行时层析成像得到的煤层内的异常体信息不准确的问题。
上述说明仅是本申请技术方案的概述，为了能够更清楚了解本申请的技术手段，而可依照说明书的内容予以实施，并且为了让本申请的上述和其它目的、特征和优点能够更明显易懂，以下特举本申请的具体实施方式。
附图说明
图1为本发明实施例中一种槽波频散曲线的变异度计算方法的流程示意图；
图2为本发明实施例中槽波频谱曲线图；
图3为本发明实施例中槽波实际频散曲线、槽波频谱曲线、槽波理论频散曲线混合图；
图4为本发明实施例中一种异常体的识别方法的流程示意图；
图5为本发明实施例中激发点与接收点的射线路径图；
图6为本发明实施例中待测区域网格化图。
图7为本发明实施例中一种计算机的硬件结构示意图。
具体实施方式
本发明实施例提供了一种槽波频散曲线的变异度计算方法，如图1所示，包括：
S101.获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据；具体地，地震数据为在煤层内激发地震波时记录的数据，包括激发点与接收点的距离、激发地震波的开始时刻，在接收点接收到的槽波数据等。
S102.根据地震数据计算槽波实际频散曲线及槽波频谱曲线；具体地，可以将地震数据进行窄带滤波，得到槽波实际频散曲线，将地震数据进行傅里叶变换，并进行归一化处理，得到槽波频谱曲线。
S103.根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线；具体地，可以根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据解频散方程，得到槽波理论频散曲线。
S104.根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度。具体地，可定义槽波实际频散曲线变异度计算表达式为：
Figure BDA0002756034040000061
其中，d为变异度，R为接收点(k＝1,2，…，K)，S为激发点(j＝1,2，……，J)，v为槽波速度，f为频率，n为断点，断点总数为N、i为第几个断点，wi为第i个断点的加权因子，
Figure BDA0002756034040000062
Figure BDA0002756034040000063
为第i个断点对应的频谱第一振幅值，
Figure BDA0002756034040000064
为第i个断点对应的频谱第二振幅值，wi<1，V为围岩的横波速度，F为槽波频谱曲线中预设振幅值对应的最大频率fb，即第三频率，f2i-1为第i个断点对应的第一频率，f2i为第i个断点对应的第二频率，v2i为第i个断点对应的第二频率f2i对应的槽波实际速度值，v′2i为第i个断点对应的第二频率f2i对应的槽波理论速度值。槽波实际频散曲线变异度计算表达式等号右边第一项表示速度(时间)间断性大小，第二项表示频率间断性大小。根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线变异度计算表达式中各变量，可计算得到槽波实际频散曲线的变异度。槽波实际频散曲线变异度计算表达式表明，断点数多，变异度大，并且d值大，变异度大，地质问题突出，d值小，变异度小，地质问题越减弱。
本发明实施例提供的槽波频散曲线的变异度计算方法，通过获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据，根据地震数据计算槽波实际频散曲线及槽波频谱曲线，根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线，根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度，从而可以对频散曲线的变异度进行定量计算，从而可以根据待测区域内各频散曲线的变异度对待测区域的异常体进行成像，从而识别待测区域的异常体，由于可以准确地定量待测区域内各频散曲线的变异度，从而根据待测区域内各频散曲线的变异度成像后得到的待测区域的异常体信息更加准确，从而可以解决现有技术中当槽波射线经过变异体时，槽波频散曲线连续性变差，从而不能准确拾取旅行时，使得通过拾取旅行时层析成像得到的煤层内的异常体信息不准确的问题。
在可选的实施例中，根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度，具体包括：根据槽波实际频散曲线确定槽波实际频散曲线中各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，第一频率小于第二频率；根据槽波频谱曲线中预设振幅值对应的最大频率确定第三频率；根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点的加权因子；根据槽波实际频散曲线及槽波理论频散曲线确定槽波实际频散曲线中各断点对应的第二频率对应的槽波理论速度值；根据各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，第三频率，各断点的加权因子及各断点对应的第二频率对应的槽波理论速度值计算槽波实际频散曲线的变异度。
具体地，槽波频谱曲线示例性的可如图2所示，在图2中，可得到相应的主频fm，主频振幅值。预设振幅值可设定为主频振幅值的0.707倍，从而可以得到预设振幅值对应的频率fa和fb，其中fa小于fb，从而可以根据fb确定第三频率。并且可以确定fa和fb两者之差为频带宽度。将槽波频谱曲线及槽波理论频散曲线投影至槽波实际频散曲线中，示例性的可如图3所示。在图3中，根据预设断点分类表可确定槽波实际频散曲线各断点的断点类型，可确定槽波实际频散曲线的断点总数N，各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，例如，第1个断点n1对应的第一频率为f1，对应的第二频率为f2，第二频率f2对应的槽波实际速度值为v2。根据槽波实际频散曲线及槽波频谱曲线可确定槽波实际频散曲线中各断点对应的频谱第一振幅值和频谱第二振幅值，从而可确定加权因子。根据槽波实际频散曲线及槽波理论频散曲线可确定槽波实际频散曲线中各断点对应的第二频率对应的槽波理论速度值，例如，第1个断点n1对应的第二频率f2对应的槽波理论速度值为v′ 2。
断点分类表中可包括Ⅰ、Ⅱ、Ⅲ类断点。Ⅰ断点，曲线一致，只是频率f断开，如图3中第1个断点所示。Ⅱ类断点，曲线上下错动，但频率f不间断，如图3中第2个断点所示。Ⅲ类断点，曲线上下错动，且频率f间断，如图3中第3个断点所示。根据槽波实际频散曲线变异度计算表达式，以及各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，第三频率，各断点的加权因子及各断点对应的第二频率对应的槽波理论速度值可以计算槽波实际频散曲线的变异度；其中，对于Ⅰ断点，有v′2i＝v2i，此时槽波实际频散曲线变异度计算表达式中第一项为0，第二项不为0。对于Ⅱ类断点，有v’2i≠v2i，此时槽波实际频散曲线变异度计算表达式中第一项不为0，第二项为0。对于Ⅲ类断点，有v’2i≠v2i，此时槽波实际频散曲线变异度计算表达式中第一项不为0，第二项不为0，从而断点数越多，变异度越大。
通过根据各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，第三频率，各断点的加权因子及各断点对应的第二频率对应的槽波理论速度值计算槽波实际频散曲线的变异度，考虑了槽波实际频散曲线中速度的间断性大小和频率间断性大小，从而可以根据速度的间断性大小和频率间断性大小对槽波实际频散曲线的变异度进行定量。
在可选的实施例中，根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点的加权因子，包括：根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点对应的频谱第一振幅值及频谱第二振幅值；根据各断点对应的频谱第一振幅值及频谱第二振幅值确定槽波实际频散曲线中各断点的加权因子。
具体地，如图3所示，各断点对应的频谱第一振幅值为各断点对应的第一频率对应的频谱第一振幅值，各断点对应的频谱第二振幅值为各断点对应的第二频率对应的频谱第二振幅值，选取各断点对应的频谱第一振幅值和频谱第二振幅值中较大的值作为各断点的加权因子。
通过根据槽波实际频散曲线及槽波频谱曲线确定槽波实际频散曲线中各断点对应的频谱第一振幅值及频谱第二振幅值；根据各断点对应的频谱第一振幅值及频谱第二振幅值确定槽波实际频散曲线中各断点的加权因子，从而考虑了槽波实际频散曲线各断点的振幅，从而可以使计算到的槽波实际频散曲线变异度更加准确。
本发明实施例还提供了一种异常体的识别方法，通过如上述任意实施方式中的槽波频散曲线的变异度计算方法计算槽波频散曲线的变异度，如图4所示，异常体的识别方法包括：
S401.获取待测区域位置数据、各槽波射线路径数据及各地震数据；具体地，待测区域位置数据包括待测区域的长度、宽度、形状等。槽波射线路径数据如图5所示，可以包括激发点Sj位置、接收点Rk位置、激发点与接收点之间的射线路径。地震数据为在煤层内激发地震波时记录的数据，包括激发点与接收点的距离、激发地震波的开始时刻，在接收点接收到的槽波数据等。还可以包括煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据。
S402.根据各地震数据计算各槽波频散曲线的变异度；具体地，可以通过上述任意实施方式中的槽波频散曲线的变异度计算方法计算各槽波频散曲线的变异度，具体的实施步骤不再赘述。
S403.根据待测区域位置数据对待测区域进行网格划分；具体地，对待测区域进行网格划分时，网格的长度、宽度可以不相等。
S404.根据网格划分结果、各槽波射线路径数据计算各网格区域中各槽波的路径系数；具体地，示例性的，如图6所示，可以根据各槽波射线路径数据将各激发点与各接收点之间的射线路径投影至划分网格后的待测区域内，并根据各射线路径经过各网格的长度l计算各网格区域中各槽波的路径系数g。
S405.根据各槽波频散曲线的变异度及各网格区域中各槽波的路径系数计算各网格区域的槽波异度值；具体地，如图6所示，可以将各槽波频散曲线的变异度
Figure BDA0002756034040000101
均匀地投影到其射线路径经过的网格区域内。将各网格区域中各槽波频散曲线的变异度乘以其对应的槽波的路径系数，并进行累加，可以计算到各网格区域的槽波异度值。
S406.根据各网格区域的槽波异度值对待测区域的异常体进行识别。具体地，槽波异度值越大，异常体的可能性越大。可以将各网格区域中的槽波异度值按照最高槽波异度值进行归一化，使其归一化槽波异度值区间为[0,1]，然后通过异常体分布区对应的预设阈值对异常体进行识别，得到异常体信息。
本发明实施例提供的异常体的识别方法，通过获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据，根据地震数据计算槽波实际频散曲线及槽波频谱曲线，根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线，根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度，从而可以对频散曲线的变异度进行定量计算，从而可以根据待测区域内各频散曲线的变异度对待测区域的异常体进行成像，从而识别待测区域的异常体，由于可以准确地定量待测区域内各频散曲线的变异度，从而根据待测区域内各频散曲线的变异度成像后得到的待测区域的异常体信息更加准确，从而可以解决现有技术中当槽波射线经过变异体时，槽波频散曲线连续性变差，从而不能准确拾取旅行时，使得通过拾取旅行时层析成像得到的煤层内的异常体信息不准确的问题。
在可选的实施例中，根据网格划分结果、各槽波射线路径数据计算各网格区域中各槽波的路径系数，包括：根据网格划分结果、各槽波射线路径数据确定各网格区域中各射线路径经过的长度及各网格的对角线长度；根据各网格区域中各射线路径经过的长度及各网格的对角线长度计算各网格区域中各槽波的路径系数。
具体地，如图6所示，可以根据各槽波射线路径数据将各激发点与各接收点之间的射线路径投影至划分网格后的待测区域内，然后根据各网格的长度、宽度可以计算各网格的对角线长度p以及各网格区域中各射线路径经过的长度l。根据各网格区域中各射线路径经过的长度l及各网格的对角线长度p的比值，可以计算得到各网格区域中各槽波的路径系数g。
在可选的实施例中，根据各槽波频散曲线的变异度及各网格区域中各槽波的路径系数计算各网格区域的槽波异度值，包括：根据各槽波频散曲线的变异度与各网格区域中各槽波的路径系数计算各网格区域中各槽波的异度值；将各网格区域中各槽波的异度值进行累加得到各网格区域的槽波异度值。
具体地，根据各槽波频散曲线的变异度与各网格区域中各槽波的路径系数的乘积可以计算各网格区域中各槽波的异度值。通过将各网格区域中各槽波的异度值进行累加得到各网格区域的槽波异度值，可以计算得到各网格区域中所以的槽波频散曲线对应的槽波异度值，从而使得各网格区域的槽波异度值更加准确。
在可选的实施例中，根据各网格区域的槽波异度值对待测区域的异常体进行识别，具体包括：根据各网格区域的槽波异度值及槽波异度值预设等级对应的预设阈值确定各网格区域的槽波异度值所属的等级；根据各网格区域的槽波异度值所属的等级对待测区域的异常体进行识别。
具体地，可以对槽波异度值进行等级设定，可以分为高值区及低值区。高值区代表异常体分布区，低值区代表无异常体分布区。高值区对应的预设阈值的区间可以为[0.7,1]，低值区对应的预设阈值的区间可以为[0,0.7)，但并不以此为限，预设阈值的区间可以根据实际情况进行具体设定。根据各网格区域的槽波异度值及槽波异度值预设等级对应的预设阈值可以确定各网格区域的槽波异度值所属的等级，并可以通过槽波异度值预设等级对应的成像图像将各网格区域进行成像，例如，高值区用红色图像标识，低值区用蓝色图像标识。从而根据各网格区域的成像可以判断各网格区域是否为异常体分布区，从而可以根据异常体分布区确定异常体。
通过根据各网格区域的槽波异度值及槽波异度值预设等级对应的预设阈值确定各网格区域的槽波异度值所属的等级，根据各网格区域的槽波异度值所属的等级对待测区域的异常体进行识别，可以比较直观、形象地识别待测区域的异常体。
基于与前述实施例中一种槽波频散曲线的变异度计算方法或一种异常体的识别方法同样的发明构思，本发明还提供一种计算机，其上存储有计算机程序，该程序被处理器执行时实现前文所述一种槽波频散曲线的变异度计算方法或一种异常体的识别方法的任一方法的步骤。
其中，在图7中，总线架构(用总线300来代表)，总线300可以包括任意数量的互联的总线和桥，总线300将包括由处理器302代表的一个或多个处理器和存储器304代表的存储器的各种电路链接在一起。总线300还可以将诸如外围设备、稳压器和功率管理电路等之类的各种其他电路链接在一起，这些都是本领域所公知的，因此，本文不再对其进行进一步描述。总线接口306在总线300和接收器301和发送器303之间提供接口。接收器301和发送器303可以是同一个元件，即收发机，提供用于在传输介质上与各种其他装置通信的单元。
处理器302负责管理总线300和通常的处理，而存储器304可以被用于存储处理器302在执行操作时所使用的数据。
基于与前述实施例中一种槽波频散曲线的变异度计算方法或一种异常体的识别方法同样的发明构思，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现以下步骤：
获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据；根据地震数据计算槽波实际频散曲线及槽波频谱曲线；根据煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线；根据槽波实际频散曲线、槽波频谱曲线及槽波理论频散曲线计算槽波实际频散曲线的变异度。
或该程序被处理器执行时实现以下步骤：
获取待测区域位置数据、各槽波射线路径数据及各地震数据；根据各地震数据计算各槽波频散曲线的变异度；根据待测区域位置数据对待测区域进行网格划分；根据网格划分结果、各槽波射线路径数据计算各网格区域中各槽波的路径系数；根据各槽波频散曲线的变异度及各网格区域中各槽波的路径系数计算各网格区域的槽波异度值；根据各网格区域的槽波异度值对待测区域的异常体进行识别。
在具体实施过程中，该程序被处理器执行时，还可以实现上述实施例中的任一方法步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程信息处理设备的处理器以产生一个机器，使得通过计算机或其他可编程信息处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程信息处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程信息处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (10)

1.一种槽波频散曲线的变异度计算方法，其特征在于，包括：
获取地震数据、煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据；
根据所述地震数据计算槽波实际频散曲线及槽波频谱曲线；
根据所述煤层厚度数据、煤层的横波速度数据及围岩的横波速度数据计算槽波理论频散曲线；
根据所述槽波实际频散曲线、所述槽波频谱曲线及所述槽波理论频散曲线计算所述槽波实际频散曲线的变异度。
2.根据权利要求1所述的槽波频散曲线的变异度计算方法，其特征在于，所述根据所述地震数据计算实际槽波频散曲线及槽波频谱曲线，包括：
将所述地震数据进行窄带滤波，得到槽波实际频散曲线；
将所述地震数据进行傅里叶变换，得到槽波频谱曲线。
3.根据权利要求1所述的槽波频散曲线的变异度计算方法，其特征在于，所述根据所述槽波实际频散曲线、所述槽波频谱曲线及所述槽波理论频散曲线计算所述槽波实际频散曲线的变异度，包括：
根据所述槽波实际频散曲线确定所述槽波实际频散曲线中各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，所述第一频率小于第二频率；
根据所述槽波频谱曲线中预设振幅值对应的最大频率确定第三频率；
根据所述槽波实际频散曲线及所述槽波频谱曲线确定所述槽波实际频散曲线中各断点的加权因子；
根据所述槽波实际频散曲线及所述槽波理论频散曲线确定所述槽波实际频散曲线中各断点对应的第二频率对应的槽波理论速度值；
根据所述各断点对应的第一频率、第二频率及第二频率对应的槽波实际速度值，所述第三频率，所述各断点的加权因子及所述各断点对应的第二频率对应的槽波理论速度值计算所述槽波实际频散曲线的变异度。
4.根据权利要求3所述的槽波频散曲线的变异度计算方法，其特征在于，所述根据所述槽波实际频散曲线及所述槽波频谱曲线确定所述槽波实际频散曲线中各断点的加权因子，包括：
根据所述槽波实际频散曲线及所述槽波频谱曲线确定所述槽波实际频散曲线中各断点对应的频谱第一振幅值及频谱第二振幅值；
根据所述各断点对应的频谱第一振幅值及频谱第二振幅值确定所述槽波实际频散曲线中各断点的加权因子。
5.一种异常体的识别方法，其特征在于，通过如权利要求1-4任意一项所述的槽波频散曲线的变异度计算方法计算槽波频散曲线的变异度，所述方法包括：
获取待测区域位置数据、各槽波射线路径数据及各地震数据；
根据各所述地震数据计算各槽波频散曲线的变异度；
根据所述待测区域位置数据对所述待测区域进行网格划分；
根据网格划分结果、各所述槽波射线路径数据计算各网格区域中各槽波的路径系数；
根据各所述槽波频散曲线的变异度及各网格区域中各所述槽波的路径系数计算各网格区域的槽波异度值；
根据所述各网格区域的槽波异度值对待测区域的异常体进行识别。
6.根据权利要求5所述的异常体的识别方法，其特征在于，所述根据网格划分结果、各所述槽波射线路径数据计算各网格区域中各槽波的路径系数，包括：
根据网格划分结果、各所述槽波射线路径数据确定各网格区域中各射线路径经过的长度及各网格的对角线长度；
根据各网格区域中各所述射线路径经过的长度及所述各网格的对角线长度计算各网格区域中各所述槽波的路径系数。
7.根据权利要求5所述的异常体的识别方法，其特征在于，所述根据各所述槽波频散曲线的变异度及各网格区域中各所述槽波的路径系数计算各网格区域的槽波异度值，包括：
根据各所述槽波频散曲线的变异度与各网格区域中各所述槽波的路径系数计算各网格区域中各槽波的异度值；
将各网格区域中各槽波的异度值进行累加得到各网格区域的槽波异度值。
8.根据权利要求5所述的异常体的识别方法，其特征在于，所述根据所述各网格区域的槽波异度值对待测区域的异常体进行识别，包括：
根据所述各网格区域的槽波异度值及槽波异度值预设等级对应的预设阈值确定所述各网格区域的槽波异度值所属的等级；
根据所述各网格区域的槽波异度值所属的等级对待测区域的异常体进行识别。
9.一种计算机，其特征在于，包括：
至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行如权利要求1-4任意一项所述的槽波频散曲线的变异度计算方法或如权利要求5-8任意一项所述的异常体的识别方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行如权利要求1-4任意一项所述的槽波频散曲线的变异度计算方法或如权利要求5-8任意一项所述的异常体的识别方法。
CN202011202939.8A 2020-11-02 2020-11-02 一种槽波频散曲线的变异度计算方法及异常体的识别方法 Active CN112346120B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011202939.8A CN112346120B (zh) 2020-11-02 2020-11-02 一种槽波频散曲线的变异度计算方法及异常体的识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011202939.8A CN112346120B (zh) 2020-11-02 2020-11-02 一种槽波频散曲线的变异度计算方法及异常体的识别方法

Publications (2)

Publication Number Publication Date
CN112346120A true CN112346120A (zh) 2021-02-09
CN112346120B CN112346120B (zh) 2021-05-18

Family

ID=74355323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011202939.8A Active CN112346120B (zh) 2020-11-02 2020-11-02 一种槽波频散曲线的变异度计算方法及异常体的识别方法

Country Status (1)

Country Link
CN (1) CN112346120B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987095A (zh) * 2021-02-25 2021-06-18 中国科学院地理科学与资源研究所 一种地质断层探测方法及装置
CN113156504A (zh) * 2021-04-09 2021-07-23 中国科学院地理科学与资源研究所 一种地震波速度的确定方法及装置
CN114509813A (zh) * 2022-03-03 2022-05-17 中国科学院地理科学与资源研究所 基于槽波确定煤层厚度的方法、装置及电子设备
CN116166988A (zh) * 2023-04-25 2023-05-26 山东科技大学 一种含煤地层异常构造识别分类方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040230389A1 (en) * 2003-04-11 2004-11-18 Westerngeco L.L.C., A Delaware Limited Liability Company Determination of waveguide parameters
CN102788991A (zh) * 2012-07-16 2012-11-21 中煤科工集团西安研究院 基于z分量的透射槽波高密度快速探测方法
CN111077572A (zh) * 2019-12-17 2020-04-28 安徽理工大学 一种基于透射槽波频散曲线反演的煤厚定量预测方法
CN111551989A (zh) * 2020-05-20 2020-08-18 中国科学院地理科学与资源研究所 一种透射槽波成像方法、设备及计算机可读存储介质

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040230389A1 (en) * 2003-04-11 2004-11-18 Westerngeco L.L.C., A Delaware Limited Liability Company Determination of waveguide parameters
CN102788991A (zh) * 2012-07-16 2012-11-21 中煤科工集团西安研究院 基于z分量的透射槽波高密度快速探测方法
CN111077572A (zh) * 2019-12-17 2020-04-28 安徽理工大学 一种基于透射槽波频散曲线反演的煤厚定量预测方法
CN111551989A (zh) * 2020-05-20 2020-08-18 中国科学院地理科学与资源研究所 一种透射槽波成像方法、设备及计算机可读存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵云佩: "槽波反射法在探测采空区中的应用", 《世界核地质科学》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987095A (zh) * 2021-02-25 2021-06-18 中国科学院地理科学与资源研究所 一种地质断层探测方法及装置
CN112987095B (zh) * 2021-02-25 2021-10-15 中国科学院地理科学与资源研究所 一种地质断层探测方法及装置
CN113156504A (zh) * 2021-04-09 2021-07-23 中国科学院地理科学与资源研究所 一种地震波速度的确定方法及装置
CN114509813A (zh) * 2022-03-03 2022-05-17 中国科学院地理科学与资源研究所 基于槽波确定煤层厚度的方法、装置及电子设备
CN116166988A (zh) * 2023-04-25 2023-05-26 山东科技大学 一种含煤地层异常构造识别分类方法
CN116166988B (zh) * 2023-04-25 2024-02-02 山东科技大学 一种含煤地层异常构造识别分类方法

Also Published As

Publication number Publication date
CN112346120B (zh) 2021-05-18

Similar Documents

Publication Publication Date Title
CN112346120B (zh) 一种槽波频散曲线的变异度计算方法及异常体的识别方法
DE102016222821B4 (de) Radarvorrichtung für ein Fahrzeug und Zielobjektmessverfahren dafür
Lederman et al. A data fusion approach for track monitoring from multiple in-service trains
DE102009043882B4 (de) Vorrichtung und Verfahren zum Abschwächen von Mehrwegsignalen
EP2930533B1 (en) Systems for a code carrier divergence high-pass filter monitor
US11867858B2 (en) Seismic observation device, seismic observation method, and recording medium
DE102013227199B4 (de) Gerät und Verfahren für das Entfernen von Rauschen eines Ultraschallsystems
CN107993001A (zh) 一种风险评估的可视化方法、装置及存储介质
DE10236900A1 (de) Verfahren und System zum Bitfehlerratentesten in minimaler Zeit
CN102436742A (zh) 浮动车系统交通信息服务水平评价方法及装置
EP2239598B1 (de) Verfahren und Vorrichtung zur Abschätzung des Integritätsrisikos in einem Satellitennavigationssystem
CN104880713A (zh) Gbas完好性评估处理方法和系统
DE102017203543B4 (de) Verfahren zum Empfangen und Überwachen eines Signals sowie eine Vorrichtung zum Empfangen und Überwachen von Signalen
US20220268591A1 (en) Method, Computer Program, and Device for Processing Data Recorded by a Motor Vehicle
CN115510998A (zh) 交易异常值检测方法及装置
CN115269389A (zh) 一种项目质量确定方法、装置、电子设备及存储介质
DE102015209939A1 (de) Verfahren und Vorrichtung zum Bewerten einer Signifikanz von mittels eines an einem Fahrzeug angebrachten Ultraschallsensors empfangenen Ultraschallsignalen
CN109781229A (zh) 车辆超载检测记录方法、装置、计算机设备及存储介质
CN104766318A (zh) 一种水环境遥感监测校正检验方法及装置
CN104363614A (zh) 一种小区的覆盖方向异常的检测方法及装置
Rui et al. Inaccuracy in pipeline compressor station construction cost estimation
EP1835299A1 (de) Unterbrechungsdetektion bei einer Trägerphasen-Messung
CN114491912A (zh) 管道的应力腐蚀开裂开挖检测位置的确定方法及设备
DE102015103632B4 (de) Verfahren zum Überprüfen einer Auswirkung einer Veränderung von Sensordaten auf ein Fahrerassistenzsystem und programmierbare Hardwarekomponente
CN110485222B (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
GR01 Patent grant
GR01 Patent grant