CN107807356B - 一种gpr绕射波速度分析方法 - Google Patents

一种gpr绕射波速度分析方法 Download PDF

Info

Publication number
CN107807356B
CN107807356B CN201711069363.0A CN201711069363A CN107807356B CN 107807356 B CN107807356 B CN 107807356B CN 201711069363 A CN201711069363 A CN 201711069363A CN 107807356 B CN107807356 B CN 107807356B
Authority
CN
China
Prior art keywords
gpr
diffraction
distance
diffracted wave
velocity
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.)
Expired - Fee Related
Application number
CN201711069363.0A
Other languages
English (en)
Other versions
CN107807356A (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 Shiyou University
Original Assignee
Xian Shiyou University
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 Shiyou University filed Critical Xian Shiyou University
Priority to CN201711069363.0A priority Critical patent/CN107807356B/zh
Publication of CN107807356A publication Critical patent/CN107807356A/zh
Application granted granted Critical
Publication of CN107807356B publication Critical patent/CN107807356B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Traffic Control Systems (AREA)

Abstract

一种GPR绕射波速度分析方法，步骤为：1）依据绕射波时距双曲线特征，在含有m道的GPR剖面上拾取出每一记录道上绕射波旅行时t j ；2）确定发射天线与接收天线之间的距离a，在GPR剖面上确定绕射点在地面投影位置M，计算出发射天线与绕射点在地面投影位置M之间的距离x j 和接收天线与绕射点在地面投影位置M之间的距离(x j ‑a)；3）基于绕射波时距方程，将t j 、x j 、(x j ‑a)带入GPR绕射波速度分析方程中，给定一个绕射点埋藏深度值，计算每一个GPR道的绕射波速度v j ；4）绘制绕射波速度分析曲线图；5）重复第三步、第四步，直到求取的所有GPR道对应的速度值相等，该速度就是真实的GPR绕射波速度，据此速度可直接对GPR数据进行常速度偏移成像处理及地质解释。

Description

一种GPR绕射波速度分析方法
技术领域
本发明属于探地雷达（GPR）数据处理技术领域，特别涉及一种GPR绕射波速度分析方法。
背景技术
当地下存在孔洞、暗穴、采空区等地质异常体时，在GPR剖面上会出现一类与地质构造有关的绕射波，因此GPR绕射波是识别这类地质异常体的有效波。在GPR剖面法勘探过程中，要想利用绕射波对地质异常体的特征做出准确判断，还需要对GPR数据做偏移成像处理，而偏移成像处理的一个重要环节是需要获得准确的GPR成像速度模型。因此速度参数获取的准确程度将直接影响到GPR绕射波偏移成像的质量，也将直接影响到GPR地质解释的准确程度。然而GPR剖面法获得的数据不包含多次覆盖的信息，因此难以利用传统的反射波速度分析技术求取速度参数。
发明内容
为了克服上述现有技术的不足，本发明的目的在于提出了一种GPR绕射波速度分析方法，该方法不但能有效、准确获取GPR绕射波速度，而且获得的速度参数可直接作为GPR数据常速度偏移成像处理的速度模型，还可用于GPR资料地质解释。
为了实现上述目的，本发明采用的技术方案是：一种GPR绕射波速度分析方法，其特征在于，包括如下步骤：
第一步，依据绕射波时距双曲线特征，在含有m道的GPR剖面上拾取出每一记录道上绕射波旅行时t j ；
第二步，通过野外数据采集参数确定发射天线T与接收天线R之间的距离a，在GPR剖面上确定绕射点在地面投影位置M，并计算出发射天线T与绕射点在地面投影位置M之间的距离x j 和接收天线R与绕射点在地面投影位置M之间的距离(x j -a)；
第三步，基于绕射波时距方程，将t j 、x j 、(x j -a)带入到GPR绕射波速度分析方程中，在给定一个绕射点埋藏深度值h的情况下，计算每一个GPR道的绕射波速度v j ；
依据下列公式计算GPR绕射波速度v j ：
GPR绕射波时距方程：
Figure DEST_PATH_IMAGE001
GPR绕射波速度分析方程：
Figure 951779DEST_PATH_IMAGE002
式中，t j 为绕射波旅行时，v j 为绕射波速度，x j 为发射天线T与绕射点在地面投影位置M之间的距离，j为GPR道序号，j=1，2，…，m，m为GPR记录总道数，a为发射天线T与接收天线R之间的距离，h为绕射点埋藏深度；
第四步，绘制速度与GPR记录道坐标相关联的绕射波速度分析曲线图，给定一个绕射点埋藏深度值h，可绘出一条GPR绕射波速度曲线图；
第五步，重复第三步、第四步，直到求取的所有GPR道对应的速度值相等，则理论计算的GPR绕射波时距曲线与实际拾取的GPR绕射波时距曲线完全重合，此时认为给定的绕射点埋藏深度值h与实际的绕射点埋藏深度相符合，该求取的速度就是真实的GPR绕射波速度，依据此速度可直接对GPR数据进行常速度偏移成像处理及地质解释。
本发明的有益效果是：
本发明的方法基于绕射波时距曲线满足双曲线规律的特征，通过绕射波时距曲线拟合的思想求取GPR绕射波速度，有效解决了剖面法获得的GPR资料中不存在多次覆盖信息，从而导致难以利用传统的GPR反射波速度分析技术获取速度的问题。通过该方法获得的速度参数可直接作为GPR数据常速度偏移成像处理的速度模型，还可用于GPR资料地质解释。
附图说明
图1为本发明实测GPR剖面。
图2为本发明点绕射地质模型及GPR绕射波观测示意图。
图3为本发明取4个埋藏深度值对应的GPR绕射波速度分析曲线图。
图4为本发明取v=0.23 m/ns时，理论计算获得的GPR绕射波时距曲线与实际拾取的GPR绕射波时距曲线对比图。
图5为本发明取v=0.23 m/ns时，采用Kirchhoff常速度偏移技术处理获得的GPR成像剖面图。
具体实施方法
下面结合附图和实施例对本发明进一步详细说明。
参见图1、2、3、4、5，一种GPR绕射波速度分析方法，包括如下步骤：
第一步，依据绕射波时距双曲线特征，在含有m道的GPR剖面上拾取出每一记录道上绕射波旅行时t j ，参见图1和图4；
第二步，通过野外数据采集参数确定发射天线T与接收天线R之间的距离a，在GPR剖面上确定绕射点在地面投影位置M，并计算出发射天线T与绕射点在地面投影位置M之间的距离x j 和接收天线R与绕射点在地面投影位置M之间的距离(x j -a)，参见图2；
第三步，基于绕射波时距方程，将t j 、x j 、(x j -a)带入到GPR绕射波速度分析方程中，在给定一个绕射点埋藏深度值h的情况下，计算每一个GPR道的绕射波速度v j ，参见图2和图3；
依据下列公式计算GPR绕射波速度v j ：
GPR绕射波时距方程：
Figure DEST_PATH_IMAGE003
GPR绕射波速度分析方程：
Figure 557335DEST_PATH_IMAGE004
式中，t j 为绕射波旅行时，v j 为绕射波速度，x j 为发射天线T与绕射点在地面投影位置M之间的距离，j为GPR道序号，j=1，2，…，m，m为GPR记录总道数，a为发射天线T与接收天线R之间的距离，h为绕射点埋藏深度；
第四步，绘制速度与GPR记录道坐标相关联的绕射波速度分析曲线图，给定一个绕射点埋藏深度值h，可绘出一条GPR绕射波速度曲线图，参见图3；
第五步，重复第三步、第四步，直到求取的所有GPR道对应的速度值相等，则理论计算的GPR绕射波时距曲线与实际拾取的GPR绕射波时距曲线完全重合，此时认为给定的绕射点埋藏深度值h与实际的绕射点埋藏深度相符合，该求取的速度就是真实的GPR绕射波速度，依据此速度可直接对GPR数据进行常速度偏移成像处理及地质解释，参见图3、图4和图5。
实施例1
将一套含有273道、道距为0.3 m，发射天线与接收天线之间的距离为1 m的实测GPR数据为例说明实施步骤：
第一步，依据绕射波时距双曲线特征，在含有273道的GPR剖面上拾取出每一记录道上绕射波旅行时t j ，参见图1和图4；
第二步，通过野外数据采集参数确定发射天线T与接收天线R之间的距离a=1m，在GPR剖面上确定绕射点在地面投影位置M=57.0 m，并计算出发射天线T与绕射点在地面投影位置M之间的距离x j 和接收天线R与绕射点在地面投影位置M之间的距离(x j -a)，参见图2；
第三步，基于绕射波时距方程，将t j 、x j 、(x j -a)带入到GPR绕射波速度分析方程中，在给定一个绕射点埋藏深度值h的情况下，计算每一个GPR道的绕射波速度v j ，参见图2和图3；
依据下列公式计算GPR绕射波速度v j ：
GPR绕射波时距方程：
Figure 100002_DEST_PATH_IMAGE005
GPR绕射波速度分析方程：
Figure 100002_DEST_PATH_IMAGE006
式中，t j 为绕射波旅行时，v j 为绕射波速度，x j 为发射天线与绕射点在地面投影位置之间的距离，j为GPR道序号，j=1，2，…，273，a=1 m为发射天线与接收天线之间的距离，h为绕射点埋藏深度；
第四步，绘制速度与GPR记录道坐标相关联的绕射波速度分析曲线图，给定一个绕射点埋藏深度值h，可绘出一条GPR绕射波速度曲线图，参见图3；
第五步，重复第三步、第四步，直到求取的所有GPR道对应的速度值相等，则理论计算的GPR绕射波时距曲线与实际拾取的GPR绕射波时距曲线完全重合，此时认为给定的绕射点埋藏深度值h与实际的绕射点埋藏深度相符合，该求取的速度就是真实的GPR绕射波速度，依据此速度可直接对GPR数据进行常速度偏移成像处理及地质解释，参见图3、图4和图5。
实例效果说明：
图1是一套实测GPR数据，其中横坐标为距离（由道号乘以道距换算而来，单位：m），纵坐标为雷达波旅行时（单位：ns），该GPR数据中存在一组明显的且成双曲线特征分布的绕射波。
图2是点绕射地质模型及GPR绕射波观测示意图，T为发射天线，R为接收天线，M为绕射点在地面的投影位置，a为发射天线与接收天线之间的距离，x为发射天线与绕射点在地面投影位置之间的距离，h为绕射点埋藏深度，v为绕射波速度，t为旅行时。
图3是取埋藏深度h=25 m、30 m、33 m和35 m分别对应的GPR绕射波速度分析曲线图，其中横坐标为距离（单位：m），纵坐标为速度（单位：m/ns）。当h=33 m时，速度曲线接近于一条直线，其平均值为0.23 m/s。受实际GPR资料信噪比、旅行时拾取的精度、地层的复杂性、异常体的空间尺寸大小等因素的影响，即使给定的埋藏深度很精确，实际计算的速度曲线也难以达到一条理想直线的状态，但当速度曲线接近于一条直线时，认为其平均值对应的速度值就是要求取的绕射波速度。
图4是取v=0.23 m/ns时，理论计算获得的绕射波时距曲线与实际拾取的绕射波时距曲线对比图，两条绕射波时距曲线几乎重叠，说明通过GPR绕射波速度分析提取的速度值是准确的。
图5是取v=0.23 m/ns时，采用Kirchhoff常速度偏移技术处理获得的GPR成像剖面，GPR绕射波已完全收敛，异常体偏移成像质量较高，也说明通过GPR绕射波速度分析提取的速度是准确的。

Claims (2)

1.一种GPR绕射波速度分析方法，其特征在于，包括如下步骤：
第一步，依据绕射波时距双曲线特征，在含有m道的GPR剖面上拾取出每一记录道上绕射波旅行时t j ；
第二步，通过野外数据采集参数确定发射天线T与接收天线R之间的距离a，在GPR剖面上确定绕射点在地面投影位置M，并计算出发射天线T与绕射点在地面投影位置M之间的距离x j 和接收天线R与绕射点在地面投影位置M之间的距离(x j -a)；
第三步，基于绕射波时距方程，将t j 、x j 、(x j -a)带入到GPR绕射波速度分析方程中，在给定一个绕射点埋藏深度值h的情况下，计算每一个GPR道的绕射波速度v j ；
依据下列公式计算GPR绕射波速度v j ：
GPR绕射波时距方程：
Figure 807675DEST_PATH_IMAGE002
GPR绕射波速度分析方程：
Figure 193657DEST_PATH_IMAGE004
式中，t j 为绕射波旅行时，v j 为绕射波速度，x j 为发射天线T与绕射点在地面投影位置M之间的距离，j为GPR道序号，j=1，2，…，m，m为GPR记录总道数，a为发射天线T与接收天线R之间的距离，h为绕射点埋藏深度；
第四步，绘制速度与GPR记录道坐标相关联的绕射波速度分析曲线图，给定一个绕射点埋藏深度值h，可绘出一条GPR绕射波速度曲线图；
第五步，改变绕射点埋藏深度值h，然后重复第三步、第四步，直到求取的所有GPR道对应的速度值相等，则理论计算的GPR绕射波时距曲线与实际拾取的GPR绕射波时距曲线完全重合，此时认为给定的绕射点埋藏深度值h与实际的绕射点埋藏深度相符合，该求取的速度就是真实的GPR绕射波速度，依据此速度可直接对GPR数据进行常速度偏移成像处理及地质解释。
2.根据权利要求1所述的一种GPR绕射波速度分析方法，其特征在于，包括以下步骤：
第一步，依据绕射波时距双曲线特征，在含有273道的GPR剖面上拾取出每一记录道上绕射波旅行时t j ；
第二步，通过野外数据采集参数确定发射天线T与接收天线R之间的距离a=1 m，在GPR剖面上确定绕射点在地面投影位置M=57.0 m，并计算出发射天线T与绕射点在地面投影位置M之间的距离x j 和接收天线R与绕射点在地面投影位置M之间的距离(x j -a)；
第三步，基于绕射波时距方程，将t j 、x j 、(x j -a)带入到GPR绕射波速度分析方程中，在给定一个绕射点埋藏深度值h的情况下，计算每一个GPR道的绕射波速度v j ；
依据下列公式计算GPR绕射波速度v j ：
GPR绕射波时距方程：
Figure DEST_PATH_IMAGE005
GPR绕射波速度分析方程：
Figure DEST_PATH_IMAGE006
式中，t j 为绕射波旅行时，v j 为绕射波速度，x j 为发射天线与绕射点在地面投影位置之间的距离，j为GPR道序号，j=1，2，…，273，a=1 m为发射天线与接收天线之间的距离，h为绕射点埋藏深度；
第四步，绘制速度与GPR记录道坐标相关联的绕射波速度分析曲线图，给定一个绕射点埋藏深度值h，可绘出一条GPR绕射波速度曲线图；
第五步，改变绕射点埋藏深度值h，然后重复第三步、第四步，直到求取的所有GPR道对应的速度值相等，则理论计算的GPR绕射波时距曲线与实际拾取的GPR绕射波时距曲线完全重合，此时认为给定的绕射点埋藏深度值h与实际的绕射点埋藏深度相符合，该求取的速度就是真实的GPR绕射波速度，依据此速度可直接对GPR数据进行常速度偏移成像处理及地质解释。
CN201711069363.0A 2017-11-03 2017-11-03 一种gpr绕射波速度分析方法 Expired - Fee Related CN107807356B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711069363.0A CN107807356B (zh) 2017-11-03 2017-11-03 一种gpr绕射波速度分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711069363.0A CN107807356B (zh) 2017-11-03 2017-11-03 一种gpr绕射波速度分析方法

Publications (2)

Publication Number Publication Date
CN107807356A CN107807356A (zh) 2018-03-16
CN107807356B true CN107807356B (zh) 2020-10-13

Family

ID=61591656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711069363.0A Expired - Fee Related CN107807356B (zh) 2017-11-03 2017-11-03 一种gpr绕射波速度分析方法

Country Status (1)

Country Link
CN (1) CN107807356B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108627875A (zh) * 2018-05-14 2018-10-09 中国神华能源股份有限公司 雷达偏移速度确定和偏移处理方法、存储介质、电子设备
CN109343110B (zh) * 2018-09-26 2020-04-07 西安石油大学 一种均值阻滤波分离地震绕射波的方法
CN113945987B (zh) * 2021-10-15 2022-04-26 中国矿业大学(北京) 病害地质体检测的方法、装置及电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004257990A (ja) * 2003-02-27 2004-09-16 Yamatake Corp 位置検出方法および装置
CN102818852A (zh) * 2011-06-10 2012-12-12 杭州海儿科技有限公司 一种桥梁预应力孔道灌浆密实度测试方法及系统
CN105372718A (zh) * 2015-12-01 2016-03-02 西安石油大学 一种2d位场构建3d位场的方法
CN106443674A (zh) * 2016-09-23 2017-02-22 中国地质大学（武汉） 一种基于衍射和成像与最小熵技术的探地雷达波速估计方法
CN106707277A (zh) * 2017-01-16 2017-05-24 中国地质大学（北京） 基于起伏剧烈地表探地雷达数据高精度逆时偏移成像方法
CN106772593A (zh) * 2016-11-10 2017-05-31 中国矿业大学(北京) 绕射波的成像方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004257990A (ja) * 2003-02-27 2004-09-16 Yamatake Corp 位置検出方法および装置
CN102818852A (zh) * 2011-06-10 2012-12-12 杭州海儿科技有限公司 一种桥梁预应力孔道灌浆密实度测试方法及系统
CN105372718A (zh) * 2015-12-01 2016-03-02 西安石油大学 一种2d位场构建3d位场的方法
CN106443674A (zh) * 2016-09-23 2017-02-22 中国地质大学（武汉） 一种基于衍射和成像与最小熵技术的探地雷达波速估计方法
CN106772593A (zh) * 2016-11-10 2017-05-31 中国矿业大学(北京) 绕射波的成像方法及装置
CN106707277A (zh) * 2017-01-16 2017-05-24 中国地质大学（北京） 基于起伏剧烈地表探地雷达数据高精度逆时偏移成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
enhance GPR diffracted waves and its applications;hongyan shen;《international geophysical conference，qingdao，china》;20170420;1014-1016 *
偏移及速度估计在探地雷达信号处理中的研究;修志杰;《中国优秀硕士学位论文全文数据库 信息科技辑》;20070215;I136-382 *
地震散射波成像技术研究;沈鸿雁;《中国博士学位论文全文数据库 基础科学辑》;20101115;A011-21 *
探地雷达数值模拟与成像方法研究;张虎;《中国优秀硕士学位论文全文数据库 基础科学辑》;20140315;A011-175 *

Also Published As

Publication number Publication date
CN107807356A (zh) 2018-03-16

Similar Documents

Publication Publication Date Title
CN107807356B (zh) 一种gpr绕射波速度分析方法
CN101630016B (zh) 一种提高垂直地震剖面成像质量的方法
CN101609160B (zh) 一种抗背景噪音强的地震数据初至波自动拾取方法
CN102012521A (zh) 一种地震储层预测中叠前裂缝的检测方法
CN114460649B (zh) 一种深海近底拖曳式多道地震接收阵列形态重建方法
CN113625337B (zh) 一种极浅水高精度地震资料快速成像方法
CN109343118B (zh) 一种异常初至时间修正方法
CN102313901A (zh) 一种初至波迭代拾取的方法
CN105259570A (zh) 一种海上时移地震接收点位移校正方法
CN107843919B (zh) 一种微地震监测数据振幅均衡方法
CN111103621A (zh) 一种主动源共成像点叠加多道面波分析方法
CN104570122A (zh) 基于反射波的地震数据静校正方法及装置
CN113156427A (zh) 探地雷达数据的反演方法和装置
CN110389377B (zh) 基于波形互相关系数相乘的微震偏移成像定位方法
CN112180452B (zh) 基于探地雷达和三维速度谱的地下管线埋深估计方法
CN111538082B (zh) 一种地震波时频域初至自动拾取方法
CN108375789A (zh) 联合采集地震数据的同步匹配方法
CN106443674B (zh) 一种基于衍射和成像与最小熵技术的探地雷达波速估计方法
CN113075732B (zh) 一种消除高分辨率小多道地震地层异常起伏的方法
CN102759749B (zh) 一种瑞雷面波速度分析方法
CN111650645B (zh) 一种变偏移距vsp弯线校正处理方法和装置
CN107942373A (zh) 基于裂缝性油气储层断裂系统检测的相干算法
CN109490964B (zh) 一种改进的高精度avo弹性参数快速反演方法
CN112014816A (zh) 基于改进层位追踪算法的双程走时计算方法
CN118131229B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201013

Termination date: 20211103

CF01 Termination of patent right due to non-payment of annual fee