一种基于探地雷达的地下管线参数自识别方法技术

本发明专利技术提供一种基于探地雷达的地下管线参数自识别方法，该方法通过检测和处理由地下管线的回波信号，自动识别地下管线目标的位置、尺寸等信息，其实现步骤为：1.预处理，利用SVD对探地雷达接收的原始剖面图像进行变换分析，分离背景和目标信号；2.类双曲线边缘提取，对预处理图像进行边缘检测，提取管线反射形成的类双曲线同相轴的位置及形状；3.目标初始特征参数提取，对步骤2产生的类双曲线同相轴图像，进行管线空间位置、上覆地层介质平均速度和管径等参数进行提取，获得管径反演计算的初始参数；4.速度扫描及管线参数反演迭代；5.Hough变换提取管线参数。本发明专利技术是基于探地雷达探测技术，对地下管线目标体进行管径参数的自识别方法。

【技术实现步骤摘要】
【专利摘要】本专利技术提供，该方法通过检测和处理由地下管线的回波信号，自动识别地下管线目标的位置、尺寸等信息，其实现步骤为：1.预处理，利用SVD对探地雷达接收的原始剖面图像进行变换分析，分离背景和目标信号；2.类双曲线边缘提取，对预处理图像进行边缘检测，提取管线反射形成的类双曲线同相轴的位置及形状；3.目标初始特征参数提取，对步骤2产生的类双曲线同相轴图像，进行管线空间位置、上覆地层介质平均速度和管径等参数进行提取，获得管径反演计算的初始参数；4.速度扫描及管线参数反演迭代；5.Hough变换提取管线参数。本专利技术是基于探地雷达探测技术，对地下管线目标体进行管径参数的自识别方法。【专利说明】—种基于探地雷达的地下管线参数自识别方法
本专利技术属于探地雷达地下目标识别的
，涉及一种基于探地雷达的地下管线参数的自识别方法，该方法通过处理由介质不连续处散射回来的高频脉冲电磁波，达到检测地下管线目标，并对其位置、尺寸等信息进行自识别的目的。
技术介绍
地质雷达是通过发射高频脉冲电磁波进行地下目标探测的技术方法。电磁波的产生及其传播规律是进行地质雷达数据处理及解释的研究基础。电磁波是通过变化的电场与变化的磁场相互激发在介质中进行传播并形成电磁场。介质中(尤其是土壤中)电磁场、电磁波的传播——波速、衰减、反射与折射的理论是地质雷达探测检测的理论基础。地质雷达是一种快速、高效、无损探测的物探方法。常用的雷达探测方法主要以反射式探测和透射式探测这两种。反射式雷达探测原理:反射探测是雷达的发射天线和接收天线都放置于被测介质表面，通过向介质中发射高频宽带电磁波脉冲信号，并接收介质中有差异地方所反射回的电磁波信号实现探测的过程。本专利技术专利就是基于该探测原理开展的。
技术实现思路
本专利技术的目的是对浅层地下管线目标实现自识别，提供目标的定位及管径等信肩、O本专利技术提出了一种基于反射回波的浅层地下管线目标探测预处理技术。由于地表的强反射回波、设备天线间的耦合、媒质内电导率的不均匀及其中大尺寸颗粒的反射等杂波信号的存在，使得管径等目标反射信号不明显，因此对回波信号进行预处理，从中滤除背景噪声及各类杂波信号对于后续的目标检测和目标识别非常关键，预处理技术包括:零偏矫正、数字滤波、小波软阈值去噪、测线间去均值、SVD变换分离背景和目标信号等。本专利技术还给出了一种类双曲线的检测技术，该技术采用小波模极大值方法(但不局限于此，还可用方向算子、梯度算子、canny算子及Sobel等方法进行)对预处理后雷达剖面进行边缘提取，在得到类双曲线同时，最大的保留了有用信号，实现了对不连续的双曲线进行自适应拾取、检测，并采用线连接方法将较近的类双曲线相连，得到连续的双曲线。本专利技术还给出了一种基于波动反射理论的目标定位、管径尺寸及介质速度的计算方法，该方法利用已知的接收天线位置和到达时间等信息，计算出目标的位置、管径尺寸及介质速度等参数，作为初始参考值，计算方法参照下面的公式(2)。本专利技术还给出了一种空间散射点真实位置的反演方法，由于读取误差的存在，利用公式(2)提取的参数与实际值存在一定误差，只能作为反演计算的初始值，对速度在一定范围内进行扫描，反演出各扫描参数对应的散射点位置，当参数值与真实值一致时，反演出的散射点位置组合成真实的管线外壁。本专利技术还给出了一种基于Hough变换的圆检测技术，结合峰值检测技术，得到管线目标位置及管径等参数的最佳值:当参数值与真实值一致时，反演出的曲线形状最接近于真实的管线外壁，也最类似于圆，而其他情况则类似于椭圆或双曲线，以反演的结果作为Hough的输入时，当曲线越接近于圆时，参数域中的‘亮点’值越大，通过峰值检测技术，在速度扫描范围内，找到最大值，并将该值对应的速度值记录为介质速度的估计值，从而对应得到管线的目标位置及管径等参数。本专利技术的特征在于，同时对预处理后的剖面进行边缘提取，在消除杂波的基础上，最大限度的保留了原始剖面中的目标双曲线边缘信号。本专利技术的特征在于，给出了一种基于波动反射理论的目标定位、管径及介质速度的计算方法，结合多个已知接收信号的到达时间和接收天线位置进行运算，即可计算出目标的位置、管径的尺寸以及介质的速度等参数，无需已知介质的传输速度，计算方法见公式(2)，由于并对参数在一定的范围内进行扫描计算，反演出一定参数范围内的类双曲线上信号对应空间真实散射点的位置范围，由于读取误差的存在，利用公式(2)提取的参数与实际值存在一定误差，只能作为反演计算的初始值。本专利技术的特征在于，以基于波动反射理论的计算结果为参考值，给出了目标真实散射点的反演方法，计算方法见公式(3)，对速度在一定范围内进行扫描时，反演出各扫描参数对应的散射点位置，当参数值与真实值一致时，反演出的散射点位置可以组合成真实的管线外壁。本专利技术的特征在于，以速度扫描范围内的反演结果作为Hough变换圆检测的输入时，只有当参数值与真实值一致时，反演出的曲线形状才最接近于圆，而其他情况则类似于椭圆或双曲线，曲线越接近于圆时，参数域中的‘亮点’积累的值越大，通过峰值检测技术，可以得到管线目标位置及管径等参数的最佳估计值，检测结果误差小，可信度高。为了达到上述目的，本专利技术采取的步骤如下: 1)利用对接收的雷达剖面进行预处理； 2)对预处理后的剖面同时进行边缘检测，找出类双曲线信号； 3)结合目标的特征和波动反射理论，根据已知的接收天线位置及信号到达时间，求取目标的位置、管径以及介质的速度等参数的初始值； 4)进行速度扫描，并针对各扫描值，反演出各空间散射点位置信息； 5)基于Hough变换，对步骤4)中的反演结果做圆检测，得到变换后参数域的值，并对其进行峰值检测，记录峰值对应的介质速度，并得到目标位置以及管径等参数的最佳估计值。在上述技术方案中，步骤I)所述预处理是从中接收信号中滤除背景噪声及各类杂波信号。预处理技术包括:零偏矫正、数字滤波、小波软阈值去噪、测线相关去除直达波等。步骤2)利用小波模极大值方法对剖面边缘进行检测，也可以借助其他图像边缘检测方法，包括方向算子、梯度算子、canny算子及Sobel算子等。步骤3)所述对地下管径目标进行雷达探测时，在雷达剖面内，已知的是各个接收天线的坐标和它所接收到信号的时刻，不知道的是目标的位置和尺寸以及介质传播速度。如图2，设目标位置的空间坐标为(?,?),管径为R，天线第I个接收点坐标为(ζ,Ι),检测到类双曲线信号的时刻为|，介质传播的速度为,在目标散射点和第I个接收天线之间的走时方程见公式(I)。【权利要求】1.，其特征在于，可以根据探地雷达接收到管线的散射信号反演出管径目标的位置及管径，包括如下步骤: 1:预处理，首先采用零偏矫正、数字滤波、小波软阈值去噪、从接收信号中滤除背景噪声及各类杂波信号，然后利用SVD变换的优势，分离背景和目标信号； 2:双曲线边缘提取，利用小波模极大值方法对预处理后的剖面进行边缘检测，找出目标曲线位置并最大限度的保留类双曲线信号在原始剖面中的形状特征； 3:目标初始特征参数求取，结合目标的特征和波动反射理论，求取介质的速度、目标的位置以及管径等参数，由于读取误差的存在，本步骤提取参数与实际值存在一定误差，只能作为反演计算的初始参数； 4本文档来自技高网...

【技术保护点】
一种基于探地雷达的地下管线参数自识别方法，其特征在于，可以根据探地雷达接收到管线的散射信号反演出管径目标的位置及管径，包括如下步骤：1：预处理，首先采用零偏矫正、数字滤波、小波软阈值去噪、从接收信号中滤除背景噪声及各类杂波信号，然后利用SVD变换的优势，分离背景和目标信号；2：双曲线边缘提取，利用小波模极大值方法对预处理后的剖面进行边缘检测，找出目标曲线位置并最大限度的保留类双曲线信号在原始剖面中的形状特征；3：目标初始特征参数求取，结合目标的特征和波动反射理论，求取介质的速度、目标的位置以及管径等参数，由于读取误差的存在，本步骤提取参数与实际值存在一定误差，只能作为反演计算的初始参数；4：以3中计算的参数作为初始参数，在一定的正负范围内，对速度参数值进行等间隔扫描，并针对每次的扫描值，计算出相应的位置和管径参数，反演出散射点的空间位置；5：基于Hough变换，对反演结果做圆检测，记录每个反演信息在参数域中对应的峰值，进行峰值检测，作为识别管径参数的依据，得到介质速度、目标位置以及管径等参数的最佳估计值，从而达到管径检测的目的。

【技术特征摘要】

【专利技术属性】
技术研发人员：郑晶，杨峰，杜翠，乔旭，夏云海，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京;11