CN109188504B - 一种地下管廊探测方法 - Google Patents

一种地下管廊探测方法 Download PDF

Info

Publication number
CN109188504B
CN109188504B CN201810898904.9A CN201810898904A CN109188504B CN 109188504 B CN109188504 B CN 109188504B CN 201810898904 A CN201810898904 A CN 201810898904A CN 109188504 B CN109188504 B CN 109188504B
Authority
CN
China
Prior art keywords
wave
acquisition
seismic
pipe gallery
point
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
CN201810898904.9A
Other languages
English (en)
Other versions
CN109188504A (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.)
WUHAN ENGINEERING SCIENCE & TECHNOLOGY INSTITUTE
Original Assignee
WUHAN ENGINEERING SCIENCE & TECHNOLOGY INSTITUTE
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 WUHAN ENGINEERING SCIENCE & TECHNOLOGY INSTITUTE filed Critical WUHAN ENGINEERING SCIENCE & TECHNOLOGY INSTITUTE
Priority to CN201810898904.9A priority Critical patent/CN109188504B/zh
Publication of CN109188504A publication Critical patent/CN109188504A/zh
Application granted granted Critical
Publication of CN109188504B publication Critical patent/CN109188504B/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
    • 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/301—Analysis for determining seismic cross-sections or geostructures

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

本发明涉及一种城市地下管廊探测方法，包括如下步骤：步骤1在城市地下管廊的待探测区域上方的地面上沿管廊走向均匀布置多个地震波采集装置组，地震波采集装置组中的每个地震波采集装置为采集区域一个采样点；步骤2采集每个采集区域各采集点的地震波采集装置组所接收的地震波信号；步骤3利用每个采集区域各采集点的地震波信号进行互相关计算；步骤4计算对应的经验格林函数；步骤5提取面波频散曲线；步骤6对频散曲线的基模和高阶模进行分解；步骤7得到相应的待探测区域的速度剖面图；步骤8确定管廊位置。本发明能简便高效的对城市地下管廊进行探测。

Description

一种地下管廊探测方法
技术领域
本发明涉及工程地球物理勘探技术领域，具体涉及一种城市地下管廊探测方法。
背景技术
城市地下管廊是城市建设和发展的瓶颈，是现代化城市正常运行的基本保证。随着城市化进程的不断加快，在我国现使用的城市地下管廊中，不论是在管线材质、铺设手段、还是管线数量上都发生着翻天覆地的变化。由于地理因素和历史等原因，城市地下管廊“老”、“密”、“乱”的状况越来越严重，落后的管廊管理手段以及现有管廊探测技术和城市高速发展之间的矛盾也日益尖锐。而另一方面，在城市化高速发展的今天，城市管线探测的环境越来越苛刻。地面建筑物越来越多，地下设施越来越密集。空间、地面及地下的各类干扰多而严重，各类干扰体(源)形式多样。
常见的城市地下管廊探测方法包括电磁感应法、地质雷达法、瞬态瑞雷面波法、地震映像法、高密度电法等。但这些方法都存在着各自的缺点：电磁感应法和地质雷达法对金属管线探测效果较好，对非金属管线特别是现在广泛采用的PE管探测效果较差，瞬态瑞雷面波法和地震映像法容易受周边人为活动影响，高密度电法容易受地下环境影响，均难以获取较为准确的管廊位置信息，不能满足实际生产需要。
发明内容
本发明的目的在于提供一种城市地下管廊探测方法，该方法能简便高效的对城市地下管廊进行探测。
为解决上述技术问题，本发明公开的一种城市地下管廊探测方法，其特征在于，它包括如下步骤：
步骤1：根据已有的城市地下管廊信息，确定待探测城市地下管廊的走向方向，在城市地下管廊的待探测区域上方的地面上沿管廊走向均匀布置多个地震波采集装置组，每个地震波采集装置组的设置位置作为一个采集区域，每个地震波采集装置组中的多个地震波采集装置在地面上沿管廊走向的垂直方向均匀布置，地震波采集装置组中的每个地震波采集装置为采集区域一个采样点；
步骤2：开启地震波采集装置，采集每个采集区域各采集点的地震波采集装置组所接收的地震波信号；
步骤3：利用每个采集区域各采集点的地震波信号进行互相关计算，计算每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数，具体公式如下：
Figure GDA0002188445380000021
其中，A1(t)表示一个最侧面采集点的地震波信号，An(t)表示第n个采集点的地震波信号，n由一个最侧面采集点外的第二个采集点算起，N为A1(t)和An(t)两个地震波信号的离散数据点的点数个数之和，m表示变数，
Figure GDA0002188445380000022
表示将式中m从0取到N-1，然后求和，t表示时刻，地震波信号为离散信号，A1(m)表示一个最侧面采集点地震波信号的第m个离散点的数值，An(m+t)表示第n个采集点地震波信号的第m+t个离散点的数值；
步骤4：利用步骤3中得到的地震波互相关函数，计算对应的经验格林函数，具体操作如下：
每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数A1n(t)、经验格林函数
Figure GDA0002188445380000023
和真实格林函数G1n(t)之间的关系表示为：
Figure GDA0002188445380000024
其中，其中，
Figure GDA0002188445380000025
表示最侧面采集点的检波器作为震源，第n个采集点的检波器接收信号时的经验格林函数，G1n(t)表示最侧面采集点的检波器作为震源，第n个采集点的检波器接收信号时的真实格林函数，
Figure GDA0002188445380000031
表示第n个采集点的检波器作为震源，最侧面采集点的检波器接收信号时的经验格林函数，Gn1(t)表示第n个采集点的检波器作为震源，最侧面采集点的检波器接收信号时的真实格林函数；
由地震波互相关函数对时间微分即可得到经验格林函数；
Figure GDA0002188445380000032
An1(-t)表示每个采集区域中第n个采集点与一个最侧面采集点之间的地震波互相关函数；
Figure GDA0002188445380000033
An(m)表示第n个采集点地震波信号的第m个离散点的数值，A1(m+t)表示一个最侧面采集点地震波信号的第m+t个离散点的数值；
步骤5：利用步骤4获得的经验格林函数，提取面波频散曲线，其操作如下：
将步骤4得到的
Figure GDA0002188445380000034
组合成起来：
Figure GDA0002188445380000035
并对
Figure GDA0002188445380000036
进行二维傅里叶变换，得到频率波数谱F(f，k)，其中f表示频率，k表示波数；
根据波数的定义
Figure GDA0002188445380000037
其中，vs为横波波速；
将频率波数域记录变换到频率波速域：对频率波数域记录利用：
Figure GDA0002188445380000038
重新采样，得到频率波速谱，在频率波速域上的任意频率，从频率波速谱中取得的极值点，即可得到该频率对应相速度值，该频率指代上述频率波速域上的任意频率，从而得到该面波记录的频散曲线，该面波记录即为步骤2得到的地震波信号；
步骤6：由于待探测地下存在管廊，所以步骤5获得的频散曲线具有多模性，对频散曲线的基模和高阶模进行分解，可以得到频散曲线的基模部分和各高阶模部分；
步骤7：由于面波的特性，各个不同波长波的速度对应着地下相应深度位置的波速，于是根据步骤6所得的若干组频散曲线，可以得到相应的待探测区域的速度剖面图；
步骤8：对城市地下管廊的待探测区域进行探测，并与步骤7中得到的几幅剖面图进行对比分析，舍去不一致的剖面图，留下一致的剖面图，从剖面图中找出地下波速异常区，该地下波速异常区即对应地下管廊位置。
本发明的有益效果：
本发明采用的地震波方法，所以能分辨出电法，如电磁感应法和地质雷达法等难以分辨的PE管；同时，采用被动源方法，便不会有主动源法，如瞬态瑞雷面波法和地震映像法等，会有的受人工活动干扰较大的问题；另一方面，在本方法最后一步中，用传统的地质雷达或瞬态瑞雷面波法对本方法结果进行修正和佐证，便可以准确获取地下管廊的位置和深度。
附图说明
图1为本发明的管廊探测方法流程图
图2为采集装置布置的主视示意图；
图3为采集装置布置的侧视示意图；
图4为采集到的地震波信号图；
图5为基模处理得到的剖面图。
图5是基模处理的结果，从图中可以看到高速异常区(圈出来的深色部分)，即是管廊的准确位置。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细说明：
一种城市地下管廊探测方法，如图1～5所示，它包括如下步骤：
步骤1：根据已有的城市地下管廊信息，确定待探测城市地下管廊的走向方向，在城市地下管廊的待探测区域上方的地面上沿管廊走向均匀布置多个地震波采集装置组，每个地震波采集装置组的设置位置作为一个采集区域，每个地震波采集装置组中的多个地震波采集装置在地面上沿管廊走向的垂直方向均匀布置，地震波采集装置组中的每个地震波采集装置为采集区域一个采样点；
步骤2：开启地震波采集装置，采集每个采集区域各采集点的地震波采集装置组所接收的地震波信号；
步骤3：利用每个采集区域各采集点的地震波信号进行互相关计算，计算每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数(相当于在最侧面采集装置处布置了一个震源，在其它采集装置出布置了接收器，提高后续计算的准确性)，具体公式如下：
Figure GDA0002188445380000051
其中，A1(t)表示一个最侧面采集点的地震波信号，An(t)表示第n个采集点的地震波信号，n由一个最侧面采集点外的第二个采集点算起，N为A1(t)和An(t)两个地震波信号的离散数据点的点数个数之和，m表示变数，
Figure GDA0002188445380000052
表示将式中m从0取到N-1，然后求和，t表示时刻，地震波信号为离散信号，A1(m)表示一个最侧面采集点地震波信号的第m个离散点的数值，An(m+t)表示第n个采集点地震波信号的第m+t个离散点的数值；
步骤4：利用步骤3中得到的地震波互相关函数，计算对应的经验格林函数(真实格林函数很难直接获取，在这里可以用经验格林函数来代替真实格林函数进行运算和处理)，具体操作如下：
每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数A1n(t)、经验格林函数
Figure GDA0002188445380000061
和真实格林函数G1n(t)之间的关系表示为：
Figure GDA0002188445380000062
其中，其中，
Figure GDA0002188445380000063
表示最侧面采集点的检波器作为震源，第n个采集点的检波器接收信号时的经验格林函数，G1n(t)表示最侧面采集点的检波器作为震源，第n个采集点的检波器接收信号时的真实格林函数，
Figure GDA0002188445380000064
表示第n个采集点的检波器作为震源，最侧面采集点的检波器接收信号时的经验格林函数，Gn1(t)表示第n个采集点的检波器作为震源，最侧面采集点的检波器接收信号时的真实格林函数；
由地震波互相关函数对时间微分即可得到经验格林函数；
Figure GDA0002188445380000065
An1(-t)表示每个采集区域中第n个采集点与一个最侧面采集点之间的地震波互相关函数；
Figure GDA0002188445380000066
An(m)表示第n个采集点地震波信号的第m个离散点的数值，A1(m+t)表示一个最侧面采集点地震波信号的第m+t个离散点的数值；
步骤5：利用步骤4获得的经验格林函数，提取面波频散曲线，其操作如下：
将步骤4得到的
Figure GDA0002188445380000067
组合成起来：
Figure GDA0002188445380000068
并对
Figure GDA0002188445380000069
进行二维傅里叶变换，得到频率波数谱F(f，k)，其中f表示频率，k表示波数；
根据波数的定义
Figure GDA0002188445380000071
其中，vs为横波波速；
将频率波数域记录变换到频率波速域：对频率波数域记录利用：
Figure GDA0002188445380000072
重新采样(即通过公式
Figure GDA0002188445380000073
把频率波数域转化到频率波速域)，得到频率波速谱，在频率波速域上的任意频率，从频率波速谱中取得的极值点，即可得到该频率对应相速度值，该频率指代上述频率波速域上的任意频率，从而得到该面波记录的频散曲线，该面波记录即为步骤2得到的地震波信号；
步骤6：由于待探测地下存在管廊，所以步骤5获得的频散曲线具有多模性，对频散曲线的基模和高阶模进行分解，可以得到频散曲线的基模部分和各高阶模部分；
步骤7：由于面波的特性，各个不同波长波的速度对应着地下相应深度位置的波速，于是根据步骤6所得的若干组频散曲线，可以得到相应的待探测区域的速度剖面图；
步骤8：对城市地下管廊的待探测区域进行探测，并与步骤7中得到的几幅剖面图进行对比分析，舍去不一致的剖面图，留下一致的剖面图，从剖面图中找出地下波速异常区，该地下波速异常区即对应地下管廊位置。
所述步骤8中，通过地质雷达或瞬态瑞雷面波法对城市地下管廊的待探测区域进行探测。
步骤8中，地下波速异常区表示波速高于周围波速的区域。
所述最侧面采集点为最左侧采集点或最右侧采集点。
相邻两个地震波采集装置组之间的距离范围为50～200米(在直线行进区域，在这个距离内有两个准确的端点，就可以确定管廊的具体位置)。
每个地震波采集装置组中相邻两个地震波采集装置之间的距离范围为0.15～2米(距离越小，测量精度越高，但需要考虑使用成本和效率问题)。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (6)

1.一种城市地下管廊探测方法，其特征在于，它包括如下步骤：
步骤1：根据已有的城市地下管廊信息，确定待探测城市地下管廊的走向方向，在城市地下管廊的待探测区域上方的地面上沿管廊走向均匀布置多个地震波采集装置组，每个地震波采集装置组的设置位置作为一个采集区域，每个地震波采集装置组中的多个地震波采集装置在地面上沿管廊走向的垂直方向均匀布置，地震波采集装置组中的每个地震波采集装置为采集区域一个采样点；
步骤2：开启地震波采集装置，采集每个采集区域各采集点的地震波采集装置组所接收的地震波信号；
步骤3：利用每个采集区域各采集点的地震波信号进行互相关计算，计算每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数，具体公式如下：
Figure FDA0002290648990000011
其中，A1(t)表示一个最侧面采集点的地震波信号，An(t)表示第n个采集点的地震波信号，n由一个最侧面采集点外的第二个采集点算起，N为A1(t)和An(t)两个地震波信号的离散数据点的点数个数之和，m表示变数，
Figure FDA0002290648990000012
表示将式中m从0取到N-1，然后求和，地震波信号为离散信号，A1(m)表示一个最侧面采集点地震波信号的第m个离散点的数值，An(m+t)表示第n个采集点地震波信号的第m+t个离散点的数值；
步骤4：利用步骤3中得到的地震波互相关函数，计算对应的经验格林函数，具体操作如下：
每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数A1n(t)、经验格林函数
Figure FDA0002290648990000013
和真实格林函数G1n(t)之间的关系表示为：
Figure FDA0002290648990000014
其中，
Figure FDA0002290648990000015
表示最侧面采集点的检波器作为震源，第n个采集点的检波器接收信号时的经验格林函数，G1n(t)表示最侧面采集点的检波器作为震源，第n个采集点的检波器接收信号时的真实格林函数，
Figure FDA0002290648990000021
表示第n个采集点的检波器作为震源，最侧面采集点的检波器接收信号时的经验格林函数，Gn1(t)表示第n个采集点的检波器作为震源，最侧面采集点的检波器接收信号时的真实格林函数；
由地震波互相关函数对时间微分即可得到经验格林函数；
Figure FDA0002290648990000022
An1(-t)表示每个采集区域中第n个采集点与一个最侧面采集点之间的地震波互相关函数；
Figure FDA0002290648990000023
An(m)表示第n个采集点地震波信号的第m个离散点的数值，A1(m+t)表示一个最侧面采集点地震波信号的第m+t个离散点的数值；
步骤5：利用步骤4获得的经验格林函数，提取面波频散曲线，其操作如下：
将步骤4得到的
Figure FDA0002290648990000024
组合成起来：
Figure FDA0002290648990000025
并对
Figure FDA0002290648990000026
进行二维傅里叶变换，得到频率波数谱F(f,k)，其中f表示频率，k表示波数；
根据波数的定义
Figure FDA0002290648990000027
其中，vs为横波波速；
将频率波数域记录变换到频率波速域，对频率波数域记录利用：
Figure FDA0002290648990000028
重新采样，得到频率波速谱，在频率波速域上的任意频率，从频率波速谱中取得的极值点，即可得到该频率对应相速度值，该频率指代上述频率波速域上的任意频率，从而得到该面波记录的频散曲线，该面波记录即为步骤2得到的地震波信号；
步骤6：由于待探测地下存在管廊，所以步骤5获得的频散曲线具有多模性，对频散曲线的基模和高阶模进行分解，得到频散曲线的基模部分和各高阶模部分；
步骤7：由于面波的特性，各个不同波长波的速度对应着地下相应深度位置的波速，根据步骤6所得的若干组频散曲线，得到相应的待探测区域的速度剖面图；
步骤8：对城市地下管廊的待探测区域进行探测，并与步骤7中得到的几幅剖面图进行对比分析，舍去不一致的剖面图，留下一致的剖面图，从剖面图中找出地下波速异常区，该地下波速异常区即对应地下管廊位置。
2.根据权利要求1所述的城市地下管廊探测方法，其特征在于：所述步骤8中，通过地质雷达或瞬态瑞雷面波法对城市地下管廊的待探测区域进行探测。
3.根据权利要求1所述的城市地下管廊探测方法，其特征在于：所述步骤8中，地下波速异常区表示波速高于周围波速的区域。
4.根据权利要求1所述的城市地下管廊探测方法，其特征在于：所述最侧面采集点为最左侧采集点或最右侧采集点。
5.根据权利要求1所述的城市地下管廊探测方法，其特征在于：相邻两个地震波采集装置组之间的距离范围为50～200米。
6.根据权利要求1所述的城市地下管廊探测方法，其特征在于：每个地震波采集装置组中相邻两个地震波采集装置之间的距离范围为0.25～2米。
CN201810898904.9A 2018-08-08 2018-08-08 一种地下管廊探测方法 Active CN109188504B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810898904.9A CN109188504B (zh) 2018-08-08 2018-08-08 一种地下管廊探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810898904.9A CN109188504B (zh) 2018-08-08 2018-08-08 一种地下管廊探测方法

Publications (2)

Publication Number Publication Date
CN109188504A CN109188504A (zh) 2019-01-11
CN109188504B true CN109188504B (zh) 2020-04-17

Family

ID=64920759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810898904.9A Active CN109188504B (zh) 2018-08-08 2018-08-08 一种地下管廊探测方法

Country Status (1)

Country Link
CN (1) CN109188504B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110673200A (zh) * 2019-09-09 2020-01-10 武汉市工程科学技术研究院 基于编码信号的智能管道定位装置及方法
CN112526597A (zh) * 2020-11-09 2021-03-19 中铁第四勘察设计院集团有限公司 一种微动技术的勘探方法、装置和存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102662191A (zh) * 2012-06-05 2012-09-12 上海通振建设工程有限公司 面向地下空间开发的地下管线层析成像系统
CN202794563U (zh) * 2012-08-13 2013-03-13 港华辉信工程塑料(中山)有限公司 用于地下塑料管网探测的定位系统
CN103424778A (zh) * 2013-07-22 2013-12-04 浙江省工程地震研究所 一种对敲法瑞雷波勘探方法
CN104459763A (zh) * 2014-12-05 2015-03-25 招商局重庆交通科研设计院有限公司 紧支集小波探测地下空洞位置的方法及系统
CN104678435A (zh) * 2014-10-27 2015-06-03 李欣欣 一种提取Rayleigh面波频散曲线的方法
CN106324687A (zh) * 2016-10-29 2017-01-11 西南石油大学 一种埋地铁质管线探测与精确定位方法及装置
CN107728195A (zh) * 2017-08-25 2018-02-23 武汉市工程科学技术研究院 基于无线数据传输的浅层地震勘探系统及勘探方法
CN108318918A (zh) * 2017-12-29 2018-07-24 福建省建筑设计研究院有限公司 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6549012B2 (en) * 2001-06-07 2003-04-15 Larry G. Stolarczyk Radio system for characterizing and outlining underground industrial developments and facilities
GB0724847D0 (en) * 2007-12-20 2008-01-30 Statoilhydro Method of and apparatus for exploring a region below a surface of the earth

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102662191A (zh) * 2012-06-05 2012-09-12 上海通振建设工程有限公司 面向地下空间开发的地下管线层析成像系统
CN202794563U (zh) * 2012-08-13 2013-03-13 港华辉信工程塑料(中山)有限公司 用于地下塑料管网探测的定位系统
CN103424778A (zh) * 2013-07-22 2013-12-04 浙江省工程地震研究所 一种对敲法瑞雷波勘探方法
CN104678435A (zh) * 2014-10-27 2015-06-03 李欣欣 一种提取Rayleigh面波频散曲线的方法
CN104459763A (zh) * 2014-12-05 2015-03-25 招商局重庆交通科研设计院有限公司 紧支集小波探测地下空洞位置的方法及系统
CN106324687A (zh) * 2016-10-29 2017-01-11 西南石油大学 一种埋地铁质管线探测与精确定位方法及装置
CN107728195A (zh) * 2017-08-25 2018-02-23 武汉市工程科学技术研究院 基于无线数据传输的浅层地震勘探系统及勘探方法
CN108318918A (zh) * 2017-12-29 2018-07-24 福建省建筑设计研究院有限公司 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Broadband ambient noise surface wave tomography across the United States;G. D. Bensen et al.;《JOURNAL OF GEOPHYSICAL RESEARCH》;20080515;第113卷;第1-21页 *
Electrical resistivity and seismic refraction tomography to detect buried cavities;Ettore Cardarelli et al.;《Geophysical Prospecting》;20101231;第685-695页 *
地震映像法在超深管线探测中的应用;肖顺等;《城市勘测》;20140228(第1期);第170-172页 *
地震绕射波法在管线探测上的数值模拟研究;刘文伍等;《城市勘测》;20170228(第1期);第158-163、168页 *
青藏东南缘背景噪声的瑞利波相速度层析成像及强震活动;王琼等;《中国科学: 地球科学》;20141231;第44卷(第11期);第2440-2450页 *

Also Published As

Publication number Publication date
CN109188504A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109343022B (zh) 估测层间土壤含水量的方法
WO2020248817A1 (zh) 基于探地雷达三维图像属性分析的供水管道漏损检测方法
CN103728659B (zh) 一种提高地下岩溶探测精度的方法
CN104678435A (zh) 一种提取Rayleigh面波频散曲线的方法
CN104459763A (zh) 紧支集小波探测地下空洞位置的方法及系统
CN105445800A (zh) 一种厚层砂体顶部分异岩性油藏的识别方法
CN109188504B (zh) 一种地下管廊探测方法
CN111045091B (zh) 一种玄武岩覆盖下古河道的识别定位方法
CN109884694B (zh) 一种基于挤压加窗傅里叶变换的高铁震源地震信号时频分析方法
CN102590874A (zh) 一种山地采煤沉陷水田地表裂缝探测方法
CN111812707A (zh) 一种基于微动和面波勘探的地铁沿线场地内孤石探测方法
CN111551927A (zh) 一种基于三维探地雷达的地下管道管径测量方法
Ortiz et al. Using infrasound waves to monitor tropospheric weather and crater morphology changes at Volcán Tungurahua, Ecuador
CN105467461A (zh) 一种利用二维视电阻率数据识别地质异常体的方法
CN104459789A (zh) 一种地震波初至拾取的方法
CN116774276A (zh) 一种基于量子测量的地震波采集方法与系统
CN113900145B (zh) 一种基于微动信号的地下横波速度分布的推测方法和装置
CN103376443A (zh) 探地雷达探测地面干扰快速消除方法
Sun et al. Research on detection and visualization of underground pipelines
CN113125559B (zh) 一种基于微动信息识别管线位置的方法
KR101672845B1 (ko) 소규모 배열식 지진계를 이용한 근거리 지하발파 감지 방법
Yin et al. Directional wells anti-collision technology based on detecting the drill bit vibration signal and its application in field
CN102890288B (zh) 一种地震波层速度反演方法
CN113536516A (zh) 基于偏移距向量片属性的观测系统评价方法
LU504248B1 (en) Method for detecting leaching and groundwater pollution range after acid in-situ leaching uranium mining

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