CN109188504A - 一种地下管廊探测方法 - Google Patents
一种地下管廊探测方法 Download PDFInfo
- Publication number
- CN109188504A CN109188504A CN201810898904.9A CN201810898904A CN109188504A CN 109188504 A CN109188504 A CN 109188504A CN 201810898904 A CN201810898904 A CN 201810898904A CN 109188504 A CN109188504 A CN 109188504A
- Authority
- CN
- China
- Prior art keywords
- wave
- collection point
- seismic
- piping lane
- acquisition device
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 239000006185 dispersion Substances 0.000 claims abstract description 20
- 238000004613 tight binding model Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 239000000284 extract Substances 0.000 claims abstract description 4
- 238000005314 correlation function Methods 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000001052 transient effect Effects 0.000 claims description 6
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000010835 comparative analysis Methods 0.000 claims 1
- 238000000034 method Methods 0.000 description 16
- 230000000694 effects Effects 0.000 description 3
- 230000005674 electromagnetic induction Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 206010038743 Restlessness Diseases 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052755 nonmetal Inorganic materials 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
-
- 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个采集点之间的地震波互相关函数,具体公式如下:
其中,A1(t)表示一个最侧面采集点的地震波信号,An(t)表示第n个采集点的地震波信号,n由一个最侧面采集点外的第二个采集点算起,N为A1(t)和An(t)两个地震波信号的离散数据点的点数个数之和,m表示变数,表示将式中m从1取到N-1,然后求和,t表示时刻,地震波信号为离散信号,A1(m)表示一个最侧面采集点地震波信号的第m个离散点的数值,An(m+t)表示第n个采集点地震波信号的第m+t个离散点的数值;
步骤4:利用步骤3中得到的地震波互相关函数,计算对应的经验格林函数,具体操作如下:
每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数A1n(t)、经验格林函数和真实格林函数G1n(t)之间的关系表示为:
其中,其中,表示最侧面采集点的检波器作为震源,第n个采集点的检波器接收信号时的经验格林函数,G1n(t)表示最侧面采集点的检波器作为震源,第n个采集点的检波器接收信号时的真实格林函数,表示第n个采集点的检波器作为震源,最侧面采集点的检波器接收信号时的经验格林函数,Gn1(t)表示第n个采集点的检波器作为震源,最侧面采集点的检波器接收信号时的真实格林函数;
由地震波互相关函数对时间微分即可得到经验格林函数;
An1(-t)表示每个采集区域中第n个采集点与一个最侧面采集点之间的地震波互相关函数;
An(m)表示第n个采集点地震波信号的第m个离散点的数值,A1(m+t)表示一个最侧面采集点地震波信号的第m+t个离散点的数值;
步骤5:利用步骤4获得的经验格林函数,提取面波频散曲线,其操作如下:
将步骤4得到的组合成起来:
并对进行二维傅里叶变换,得到频率波数谱F(f,k),其中f表示频率,k表示波数;
根据波数的定义
其中,vs为横波波速;
将频率波数域记录变换到频率波速域:对频率波数域记录利用:
重新采样,得到频率波速谱,在频率波速域上的任意频率,从频率波速谱中取得的极值点,即可得到该频率对应相速度值,从而得到该面波记录的频散曲线;
步骤6:由于待探测地下存在管廊,所以步骤5获得的频散曲线具有多模性,对频散曲线的基模和高阶模进行分解,可以得到频散曲线的基模部分和各高阶模部分;
步骤7:由于面波的特性,各个不同波长波的速度对应着地下相应深度位置的波速,于是根据步骤6所得的若干组频散曲线,可以得到相应的待探测区域的速度剖面图;
步骤8:对城市地下管廊的待探测区域进行探测,并与步骤7中得到的几幅剖面图进行对比分析,舍去不一致的剖面图,留下一致的剖面图,从剖面图中找出地下波速异常区,该地下波速异常区即对应地下管廊位置。
本发明的有益效果:
本发明采用的地震波方法,所以能分辨出电法,如电磁感应法和地质雷达法等难以分辨的PE管;同时,采用被动源方法,便不会有主动源法,如瞬态瑞雷面波法和地震映像法等,会有的受人工活动干扰较大的问题;另一方面,在本方法最后一步中,用传统的地质雷达或瞬态瑞雷面波法对本方法结果进行修正和佐证,便可以准确获取地下管廊的位置和深度。
附图说明
图1为本发明的管廊探测方法流程图
图2为采集装置布置的主视示意图;
图3为采集装置布置的侧视示意图;
图4为采集到的地震波信号图;
图5为基模处理得到的剖面图。
图5是基模处理的结果,从图中可以看到高速异常区(圈出来的深色部分),即是管廊的准确位置。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细说明:
一种城市地下管廊探测方法,如图1~5所示,它包括如下步骤:
步骤1:根据已有的城市地下管廊信息,确定待探测城市地下管廊的走向方向,在城市地下管廊的待探测区域上方的地面上沿管廊走向均匀布置多个地震波采集装置组,每个地震波采集装置组的设置位置作为一个采集区域,每个地震波采集装置组中的多个地震波采集装置在地面上沿管廊走向的垂直方向均匀布置,地震波采集装置组中的每个地震波采集装置为采集区域一个采样点;
步骤2:开启地震波采集装置,采集每个采集区域各采集点的地震波采集装置组所接收的地震波信号;
步骤3:利用每个采集区域各采集点的地震波信号进行互相关计算,计算每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数(相当于在最侧面采集装置处布置了一个震源,在其它采集装置出布置了接收器,提高后续计算的准确性),具体公式如下:
其中,A1(t)表示一个最侧面采集点的地震波信号,An(t)表示第n个采集点的地震波信号,n由一个最侧面采集点外的第二个采集点算起,N为A1(t)和An(t)两个地震波信号的离散数据点的点数个数之和,m表示变数,表示将式中m从1取到N-1,然后求和,t表示时刻,地震波信号为离散信号,A1(m)表示一个最侧面采集点地震波信号的第m个离散点的数值,An(m+t)表示第n个采集点地震波信号的第m+t个离散点的数值;
步骤4:利用步骤3中得到的地震波互相关函数,计算对应的经验格林函数(真实格林函数很难直接获取,在这里可以用经验格林函数来代替真实格林函数进行运算和处理),具体操作如下:
每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数A1n(t)、经验格林函数和真实格林函数G1n(t)之间的关系表示为:
其中,其中,表示最侧面采集点的检波器作为震源,第n个采集点的检波器接收信号时的经验格林函数,G1n(t)表示最侧面采集点的检波器作为震源,第n个采集点的检波器接收信号时的真实格林函数,表示第n个采集点的检波器作为震源,最侧面采集点的检波器接收信号时的经验格林函数,Gn1(t)表示第n个采集点的检波器作为震源,最侧面采集点的检波器接收信号时的真实格林函数;
由地震波互相关函数对时间微分即可得到经验格林函数;
An1(-t)表示每个采集区域中第n个采集点与一个最侧面采集点之间的地震波互相关函数;
An(m)表示第n个采集点地震波信号的第m个离散点的数值,A1(m+t)表示一个最侧面采集点地震波信号的第m+t个离散点的数值;
步骤5:利用步骤4获得的经验格林函数,提取面波频散曲线,其操作如下:
将步骤4得到的组合成起来:
并对进行二维傅里叶变换,得到频率波数谱F(f,k),其中f表示频率,k表示波数;
根据波数的定义
其中,vs为横波波速;
将频率波数域记录变换到频率波速域:对频率波数域记录利用:
重新采样(即通过公式把频率波数域转化到频率波速域),得到频率波速谱,在频率波速域上的任意频率,从频率波速谱中取得的极值点,即可得到该频率对应相速度值,从而得到该面波记录的频散曲线;
步骤6:由于待探测地下存在管廊,所以步骤5获得的频散曲线具有多模性,对频散曲线的基模和高阶模进行分解,可以得到频散曲线的基模部分和各高阶模部分;
步骤7:由于面波的特性,各个不同波长波的速度对应着地下相应深度位置的波速,于是根据步骤6所得的若干组频散曲线,可以得到相应的待探测区域的速度剖面图;
步骤8:对城市地下管廊的待探测区域进行探测,并与步骤7中得到的几幅剖面图进行对比分析,舍去不一致的剖面图,留下一致的剖面图,从剖面图中找出地下波速异常区,该地下波速异常区即对应地下管廊位置。
所述步骤8中,通过地质雷达或瞬态瑞雷面波法对城市地下管廊的待探测区域进行探测。
步骤8中,地下波速异常区表示波速高于周围波速的区域。
所述最侧面采集点为最左侧采集点或最右侧采集点。
相邻两个地震波采集装置组之间的距离范围为50~200米(在直线行进区域,在这个距离内有两个准确的端点,就可以确定管廊的具体位置)。
每个地震波采集装置组中相邻两个地震波采集装置之间的距离范围为0.15~2米(距离越小,测量精度越高,但需要考虑使用成本和效率问题)。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (6)
1.一种城市地下管廊探测方法,其特征在于,它包括如下步骤:
步骤1:根据已有的城市地下管廊信息,确定待探测城市地下管廊的走向方向,在城市地下管廊的待探测区域上方的地面上沿管廊走向均匀布置多个地震波采集装置组,每个地震波采集装置组的设置位置作为一个采集区域,每个地震波采集装置组中的多个地震波采集装置在地面上沿管廊走向的垂直方向均匀布置,地震波采集装置组中的每个地震波采集装置为采集区域一个采样点;
步骤2:开启地震波采集装置,采集每个采集区域各采集点的地震波采集装置组所接收的地震波信号;
步骤3:利用每个采集区域各采集点的地震波信号进行互相关计算,计算每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数,具体公式如下:
其中,A1(t)表示一个最侧面采集点的地震波信号,An(t)表示第n个采集点的地震波信号,n由一个最侧面采集点外的第二个采集点算起,N为A1(t)和An(t)两个地震波信号的离散数据点的点数个数之和,m表示变数,表示将式中m从1取到N-1,然后求和,地震波信号为离散信号,A1(m)表示一个最侧面采集点地震波信号的第m个离散点的数值,An(m+t)表示第n个采集点地震波信号的第m+t个离散点的数值;
步骤4:利用步骤3中得到的地震波互相关函数,计算对应的经验格林函数,具体操作如下:
每个采集区域中一个最侧面采集点与第n个采集点之间的地震波互相关函数A1n(t)、经验格林函数和真实格林函数G1n(t)之间的关系表示为:
其中,其中,表示最侧面采集点的检波器作为震源,第n个采集点的检波器接收信号时的经验格林函数,G1n(t)表示最侧面采集点的检波器作为震源,第n个采集点的检波器接收信号时的真实格林函数,表示第n个采集点的检波器作为震源,最侧面采集点的检波器接收信号时的经验格林函数,Gn1(t)表示第n个采集点的检波器作为震源,最侧面采集点的检波器接收信号时的真实格林函数;
由地震波互相关函数对时间微分即可得到经验格林函数;
An1(-t)表示每个采集区域中第n个采集点与一个最侧面采集点之间的地震波互相关函数;
An(m)表示第n个采集点地震波信号的第m个离散点的数值,A1(m+t)表示一个最侧面采集点地震波信号的第m+t个离散点的数值;
步骤5:利用步骤4获得的经验格林函数,提取面波频散曲线,其操作如下:
将步骤4得到的组合成起来:
并对进行二维傅里叶变换,得到频率波数谱F(f,k),其中f表示频率,k表示波数;
根据波数的定义
其中,vs为横波波速;
将频率波数域记录变换到频率波速域,对频率波数域记录利用:
重新采样,得到频率波速谱,在频率波速域上的任意频率,从频率波速谱中取得的极值点,即可得到该频率对应相速度值,从而得到该面波记录的频散曲线;
步骤6:由于待探测地下存在管廊,所以步骤5获得的频散曲线具有多模性,对频散曲线的基模和高阶模进行分解,得到频散曲线的基模部分和各高阶模部分;
步骤7:由于面波的特性,各个不同波长波的速度对应着地下相应深度位置的波速,根据步骤6所得的若干组频散曲线,得到相应的待探测区域的速度剖面图;
步骤8:对城市地下管廊的待探测区域进行探测,并与步骤7中得到的几幅剖面图进行对比分析,舍去不一致的剖面图,留下一致的剖面图,从剖面图中找出地下波速异常区,该地下波速异常区即对应地下管廊位置。
2.根据权利要求1所述的城市地下管廊探测方法,其特征在于:所述步骤8中,通过地质雷达或瞬态瑞雷面波法对城市地下管廊的待探测区域进行探测。
3.根据权利要求1所述的城市地下管廊探测方法,其特征在于:所述步骤8中,地下波速异常区表示波速高于周围波速的区域。
4.根据权利要求1所述的城市地下管廊探测方法,其特征在于:所述最侧面采集点为最左侧采集点或最右侧采集点。
5.根据权利要求1所述的城市地下管廊探测方法,其特征在于:相邻两个地震波采集装置组之间的距离范围为50~200米。
6.根据权利要求1所述的城市地下管廊探测方法,其特征在于:每个地震波采集装置组中相邻两个地震波采集装置之间的距离范围为0.25~2米。
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 true CN109188504A (zh) | 2019-01-11 |
CN109188504B 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) |
Cited By (2)
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 (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030048104A1 (en) * | 2001-06-07 | 2003-03-13 | Stolarczyk Larry G. | Radio system for characterizing and outlining underground industrial developments and facilities |
US20110046885A1 (en) * | 2007-12-20 | 2011-02-24 | Statoil Asa | Method of and apparatus for exploring a region below a surface of the earth |
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曲线的地下不良地质体无损探测方法及应用 |
-
2018
- 2018-08-08 CN CN201810898904.9A patent/CN109188504B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030048104A1 (en) * | 2001-06-07 | 2003-03-13 | Stolarczyk Larry G. | Radio system for characterizing and outlining underground industrial developments and facilities |
US20110046885A1 (en) * | 2007-12-20 | 2011-02-24 | Statoil Asa | Method of and apparatus for exploring a region below a surface of the earth |
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 (6)
Title |
---|
ETTORE CARDARELLI ET AL.: "Electrical resistivity and seismic refraction tomography to detect buried cavities", 《GEOPHYSICAL PROSPECTING》 * |
G. D. BENSEN ET AL.: "Broadband ambient noise surface wave tomography across the United States", 《JOURNAL OF GEOPHYSICAL RESEARCH》 * |
刘文伍等: "地震绕射波法在管线探测上的数值模拟研究", 《城市勘测》 * |
张鸿升等: "《地下管线探测原理、方法及技术》", 31 August 1998, 中国矿业大学出版社 * |
王琼等: "青藏东南缘背景噪声的瑞利波相速度层析成像及强震活动", 《中国科学: 地球科学》 * |
肖顺等: "地震映像法在超深管线探测中的应用", 《城市勘测》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110673200A (zh) * | 2019-09-09 | 2020-01-10 | 武汉市工程科学技术研究院 | 基于编码信号的智能管道定位装置及方法 |
CN112526597A (zh) * | 2020-11-09 | 2021-03-19 | 中铁第四勘察设计院集团有限公司 | 一种微动技术的勘探方法、装置和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109188504B (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106324687B (zh) | 一种埋地铁质管线探测与精确定位方法 | |
CN104459763A (zh) | 紧支集小波探测地下空洞位置的方法及系统 | |
CN109782360A (zh) | 一种低阻覆盖区深部砂体探测方法 | |
CN103969678A (zh) | 煤矿采空区无源地震探测方法 | |
CN109188504A (zh) | 一种地下管廊探测方法 | |
CN106371137B (zh) | 用于浅层异常体地震探测装置及其三维观测方法 | |
Zeng et al. | Turning a telecom fiber‐optic cable into an ultradense seismic array for rapid postearthquake response in an urban area | |
CN103499844A (zh) | 全程全空间瞬变电磁装置及测量方法 | |
CN102662191A (zh) | 面向地下空间开发的地下管线层析成像系统 | |
CN107271539A (zh) | 一种海底油气管道漏磁/甚低频接收检测定位方法及装置 | |
Ortiz et al. | Using infrasound waves to monitor tropospheric weather and crater morphology changes at Volcán Tungurahua, Ecuador | |
CN103984006A (zh) | 一种全断面观测系统的隧道超前地质探测方法 | |
CN109188542A (zh) | 一种波区相关性检测的远参考大地电磁阻抗计算方法 | |
CN102937725B (zh) | 一种基于过渡区与相叠合的位场异常边缘增强方法 | |
CN103376443A (zh) | 探地雷达探测地面干扰快速消除方法 | |
CN112180452A (zh) | 基于探地雷达和三维速度谱的地下管线埋深估计方法 | |
CN103615962A (zh) | 一种滑坡体地表位移测量方法 | |
CN113900145B (zh) | 一种基于微动信号的地下横波速度分布的推测方法和装置 | |
Sun et al. | Research on detection and visualization of underground pipelines | |
CN202649476U (zh) | 面向地下空间开发的地下管线层析成像系统 | |
Yang et al. | Long-term ambient seismic interferometry for constraining seasonal subsurface velocity variations in urban settings: a distributed acoustic sensing (DAS) case study | |
Huang et al. | Earthquake detection in the Jiangsu region, China using graphics-processing-unit-based Match & Locate and rapid earthquake association and location | |
Lu et al. | Tunnel concealed karst cave joint detection by tunnel seismic and transient electromagnetic | |
CN106569278B (zh) | 一种多道相似相干速度谱计算方法 | |
Kitov et al. | Detection, estimation of magnitude, and relative location of weak aftershocks using waveform cross-correlation: The earthquake of August 7, 2016, in the town of Mariupol |
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 |