基于探地雷达的运营期地下管道管径测定方法技术

技术编号:9841225 阅读:142 留言:0更新日期:2014-04-02 04:20
本发明专利技术是一种基于探地雷达的运营期地下管道管径测定方法,其原理利用基于探地雷达电磁波反射机理以及管道特殊圆形形状,实现在地下管道运营阶段管道内充满有耗介质时还能够在非开挖、无扰动下探测与解译地下管道管径的新方法。实现步骤为:①在确定地下管道走向、埋深、对应地表位置与合适雷达频谱参数的基础上,开展探地雷达探测并测得管道上方雷达测线方向3个测点的雷达图像;②在对雷达图像去噪处理基础上,提取3个测点坐标对应电磁波单道波形图;③从单道波形曲线峰谷值确定从3个测点到管壁的电磁波传播时间;④确定场地电磁波波速基础上,得到3个测点到管壁的传播距离;⑤利用3个测点坐标与其到管壁的电磁波传播距离,计算管道管径。

【技术实现步骤摘要】
【专利摘要】本专利技术是一种,其原理利用基于探地雷达电磁波反射机理以及管道特殊圆形形状,实现在地下管道运营阶段管道内充满有耗介质时还能够在非开挖、无扰动下探测与解译地下管道管径的新方法。实现步骤为:①在确定地下管道走向、埋深、对应地表位置与合适雷达频谱参数的基础上,开展探地雷达探测并测得管道上方雷达测线方向3个测点的雷达图像;②在对雷达图像去噪处理基础上,提取3个测点坐标对应电磁波单道波形图;③从单道波形曲线峰谷值确定从3个测点到管壁的电磁波传播时间;④确定场地电磁波波速基础上,得到3个测点到管壁的传播距离;⑤利用3个测点坐标与其到管壁的电磁波传播距离,计算管道管径。【专利说明】
本专利技术属于地球物理测试的领域,涉及一种。在地下管道运营期时管道内充满有耗介质(水、油、气)使得很难测到管底土管交界处电磁反射波的情况下,基于探地雷达的电磁波反射机理与地下管道特殊的圆形形状,提出对地下管道管径进行非开挖、无扰动的现场探测与解译的方法。
技术介绍
探地雷达作为一种新兴的地球物理测试方法。它由天线、发射机、接收机、信号处理器和终端设备等几部分组成。探地雷达在探测地下管道的工作机理如下:根据电磁波在有耗介质中的传播特性,探地雷达通过发射天线向地下发射高频脉冲电磁波,当电磁波在地下介质中传播时,其传播路径、电磁场强度和波形将随着所通过介质的电性和几何形态的变化而变化,雷达接收机将对此部分的反射波进行适时接收获得测试数据和雷达图像,利用专用数据和图像处理软件对其处理,可确定地下界面、地质体空间位置和结构特性等。现场最为直接的探地雷达图像是以脉冲反射波的单道波形形式记录的,单道波形的正负峰值分别以黑白或者是灰阶或是彩色表示,同相轴(等灰度线、等色线)的差异即可形象地反映出地下反射面或是异常物体。随着雷达天线不断的探测推进,不同位置时的单道波形图按顺序显示在一起就形成了最终的雷达图像。基于探地雷达测定地下管道管径的现有方法仅适用于在非运行期地下管道为空管的情况。即:利用雷达反射波的单道波形图,在单道波形图上寻找电磁波在管道上、下壁面反射点的对应峰值点,两点的传播距离差即可间接估算管道管径。然而这一方法主要基于理想电磁波模型的理论计算,实际探测环境中雷达接收到的回波信号是复杂的并含有各类干扰信号的时域波形,电磁波在地层介质与探测目标中经过一系列反射与折射后,很难清楚准确获得管道下壁面的电磁波反射位置。尤其在管道处于运行期时内部充满着吸收电磁波的有耗介质(水、油、气),更会发生根本测不到管底部位的电磁波反射波信号的情况。
技术实现思路
技术问题:本专利技术针对在地下管道在运行期时管道内充满着吸收电磁波的有耗介质(水、油、气),使得雷达的接收天线无法接收到电磁波在管道底部的反射波信号。提出一种只利用电磁波在管道顶部的反射波信号,以及地下管道一股为圆形的几何形状特征,通过精确获得电磁波在管道顶部三个反射点的空间坐标信息,实现一种。从而解决地下管线运营期管道充满有耗介质时非开挖、无扰动探测与解译地下管道管径。技术方案:本专利技术的一种具体步骤如下:步骤1:对拟探测的区域平整场地,保证含有地下管道的场地地表水平平整;步骤2:采用低频率雷达对地下管道进行预探测,并得到拟探测地下管道的走向、埋深及其相对地表位置,从而确定对于地下管道管径解译的探地雷达精测的测线方向与雷达频谱参数,应保证探地雷达测线必须垂直于管道的走向,最后建立雷达探测与解译示意图,其中XOY坐标系以测线方向为X轴与以管道埋深为Y轴;步骤3:按照步骤2获得的测线方向与雷达频谱参数,开展针对地下管道管径解译的探地雷达精细探测,探测时雷达天线应沿着测线方向匀速推进,在地下管道相对地表位置点的正上方及其附近沿测线选择三个探测点,要保证这三个点在一条直线上,并且竖向高度一致,探测的路线上这三个测点的平面坐标分别记为(Xl,yi),(x2,y2),(x3,y3);步骤4:根据步骤3,探地雷达探测后得到含有电磁波对地下管道反射的雷达回波图像,首先对雷达图像进行去噪化频谱处理,然后进行旨在使得图像清晰化的图像边缘化处理;步骤5:确定雷达电磁波的波速:采用矩阵模板法确定雷达电磁波波速,利用雷达试验的电磁波波速生成曲线模板,然后用模板匹配目标的雷达图像,匹配最好的对应波速就是雷达波速,雷达电磁波波速模板的矩阵模板为:【权利要求】1.一种,其特征在于,具体步骤如下:步骤1:对拟探测的区域平整场地,保证含有地下管道的场地地表水平平整;步骤2:采用低频率雷达对地下管道进行预探测,并得到拟探测地下管道的走向、埋深及其相对地表位置,从而确定对于地下管道管径解译的探地雷达精测的测线方向与雷达频谱参数,应保证探地雷达测线必须垂直于管道的走向,最后建立雷达探测与解译示意图,其中XOY坐标系以测线方向为X轴与以管道埋深为Y轴;步骤3:按照步骤2获得的测线方向与雷达频谱参数,开展针对地下管道管径解译的探地雷达精细探测,探测时雷达天线应沿着测线方向匀速推进,在地下管道相对地表位置点的正上方及其附近沿测线选择三个探测点,要保证这三个点在一条直线上,并且竖向高度一致,探测的路线上这三个测点的平面坐标分别记为(Xl,yi),(x2,y2), (x3,y3);步骤4:根据步骤3,探地雷达探测后得到含有电磁波对地下管道反射的雷达回波图像,首先对雷达图像进行去噪化频谱处理,然后进行旨在使得图像清晰化的图像边缘化处理;步骤5:确定雷达电磁波的波速:采用矩阵模板法确定雷达电磁波波速,利用雷达试验的电磁波波速生成曲线模板,然后用模板匹配目标的雷达图像,匹配最好的对应波速就是雷达波速,雷达电磁波波速模板的矩阵模板为: 2.根据权利要求1所述的,其特征在于,所述的几何关系如下,假设三个测点坐标分别为(Xl,yi)> (x2, y2)> (χ3,y3),电磁波从3个测点传播到地下管道上壁面的圆半径分别为A、r2、r3,设所求地下管道的圆心点坐标为(X,y),管径的半径为r,故根据其几何关系可得到方程组: 【文档编号】G01V3/12GK103675922SQ201310684886【公开日】2014年3月26日 申请日期:2013年12月13日 优先权日:2013年12月13日 【专利技术者】张鹏, 马彬, 董韬, 王旭东 申请人:南京工业大学本文档来自技高网
...

【技术保护点】
一种基于探地雷达的运营期地下管道管径测定方法,其特征在于,具体步骤如下:步骤1:对拟探测的区域平整场地,保证含有地下管道的场地地表水平平整;步骤2:采用低频率雷达对地下管道进行预探测,并得到拟探测地下管道的走向、埋深及其相对地表位置,从而确定对于地下管道管径解译的探地雷达精测的测线方向与雷达频谱参数,应保证探地雷达测线必须垂直于管道的走向,最后建立雷达探测与解译示意图,其中XOY坐标系以测线方向为X轴与以管道埋深为Y轴;步骤3:按照步骤2获得的测线方向与雷达频谱参数,开展针对地下管道管径解译的探地雷达精细探测,探测时雷达天线应沿着测线方向匀速推进,在地下管道相对地表位置点的正上方及其附近沿测线选择三个探测点,要保证这三个点在一条直线上,并且竖向高度一致,探测的路线上这三个测点的平面坐标分别记为(x1,y1),(x2,y2),(x3,y3);步骤4:根据步骤3,探地雷达探测后得到含有电磁波对地下管道反射的雷达回波图像,首先对雷达图像进行去噪化频谱处理,然后进行旨在使得图像清晰化的图像边缘化处理;步骤5:确定雷达电磁波的波速:采用矩阵模板法确定雷达电磁波波速,利用雷达试验的电磁波波速生成曲线模板,然后用模板匹配目标的雷达图像,匹配最好的对应波速就是雷达波速,雷达电磁波波速模板的矩阵模板为: H [ i , j ] = r ( X , t ) = 1 v z l 3 u ′ ( t - 2 l v ) const 式中,r(X,t)为由费涅耳和基尔霍夫的衍射理论及球面波在介质中的传播特征得到雷达在不同位置接收到目标信号;t=dt·j,dt是采样时间间隔;X=dx·i,dx是雷达水平采样间隔;v是雷达试验电磁波波速;z是目标深度;u′为散射信号与发射信号的延时的导数,const为任意常数,不影响生成的模板中的双曲线形状,每个雷达试验波速对应一个模板;1是雷达天线到对应位置的距离;假定雷达中的目标回波曲线的图像为T[M,N],均值和方差分别为μT和σT,雷达的目标回波图像与生成的对应模板之间的匹配度定义为两个图像的相关系数: γ ( T , H ) = Σ i = 1 M Σ j = 1 N ( T [ i , j ] - μ T ) ( H [ i , j ] - μ H ) MN σ T σ H 式中,γ(T,H)为每个速度模板对应的匹配度;T[i,j]为雷达图像上任意一个像素点;H[i,j]为矩阵模板上任意一个像素点;μT和σT为T[M,N]的均值和方差;μH和σH为H[M,N]的均值和方差;M,N为图像的尺寸大小;这样每个速度模板对应的匹配度γ(T,H)可以得到,最大匹配度...

【技术特征摘要】

【专利技术属性】
技术研发人员:张鹏马彬董韬王旭东
申请(专利权)人:南京工业大学
类型:发明
国别省市:江苏;32

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1