基于瞬时幅频参数的探地雷达信号预处理方法及装置制造方法
【专利摘要】本发明公开了一种基于瞬时幅频参数的探地雷达信号预处理方法及装置。所述方法包括：提取探地雷达信号的瞬时幅度和瞬时频率，分别得到瞬时幅度矩阵Amatrix和瞬时频率矩阵fmatrix；删除Amatrix中的第一行得到对瞬时频率矩阵fmatrix进行归一化处理，得到归一化瞬时频率矩阵计算探地雷达信号的联合瞬时参数矩阵采用了本发明的探地雷达信号的瞬时幅度信息能够精确地反映分界面的埋深和各个介质层对发射脉冲的吸收能力，同时探地雷达信号的瞬时频率信息可以对晚时信号起到增强作用，能够精确地反映相邻不同介质的相对介电常数的大小信息。
【专利说明】基于瞬时幅频参数的探地雷达信号预处理方法及装置
【技术领域】
[0001]本发明涉及探地雷达【技术领域】，具体涉及一种基于瞬时幅频参数的探地雷达信号预处理方法及装置。
【背景技术】
[0002]探地雷达是利用高频超宽带信号来探测地球浅层介质分布规律的地球物理勘探装置，探地雷达信号的反演是通过分析接收信号的波形、幅值和时间等特征来估计地下介质的几何机构、埋深和空间位置等信息。不同介质的相对介电常数差别较大，探地雷达发射的脉冲信号在不同介质中传播时，其传播速度在0.01?0.3m/ns之间，由于传播速度很快，为了提高反演结果的精度，需要对探地雷达信号进行预处理，以提高相邻不同介质分界面的定位精度。
[0003]探地雷达信号具有高频率、短波长、强介质吸收和强干扰等性质，以傅里叶变换为基础的经典信号处理方法不能达到较好的预处理效果。目前披露的采用时频分析进行探地雷达信号预处理的方法包括广义S变换方法、Hilbert谱分析方法等。
[0004]广义S变换方法在地下介质的层位识别应用中效果较好，但是广义S变换的系统参数需要靠经验调节，模型可扩展性较差，限制了该方法的应用推广。Hilbert谱分析方法对薄层位置具有一定的增强作用，但实测探地雷达信号中的杂波和干扰信号也得到了增强，限制了其对浅层地下介质的反演精度。
【发明内容】
[0005]本发明提供一种基于瞬时幅频参数的探地雷达信号预处理方法及装置。采用了本发明的探地雷达信号的瞬时幅度信息能够精确地反映分界面的埋深和各个介质层对发射脉冲的吸收能力，同时探地雷达信号的瞬时频率信息可以对晚时信号起到增强作用，能够精确地反映相邻不同介质的相对介电常数的大小信息。
[0006]下面阐述本发明的技术方案。
[0007]—种基于瞬时幅频参数的探地雷达信号预处理方法，所述方法包括:
[0008]提取探地雷达信号的瞬时幅度和瞬时频率，分别得到瞬时幅度矩阵Aniato和瞬时频率矩阵fmatrix ;
[0009]删除Amatrix中的第一行得到Jlatrix；
[0010]对瞬时频率矩阵fmatHx进行归一化处理，得到归一化瞬时频率矩阵/Ittix ;
[0011]计算探地雷达信号的联合瞬时参数矩阵JmaMx，K'?/二，其中符号?
表示两个矩阵处于相同位置的元素相乘。
[0012]所述对瞬时频率矩阵fmatHx进行归一化处理的公式为:
【权利要求】
1.一种基于瞬时幅频参数的探地雷达信号预处理方法，所述方法包括: 提取探地雷达信号的瞬时幅度和瞬时频率，分别得到瞬时幅度矩阵Amataix和瞬时频率矩阵 fmatrix ； 删除Amatox中的第一行得到^atrix; 对瞬时频率矩阵fmatHx进行归一化处理，得到归一化瞬时频率矩阵； 计算探地雷达信号的联合瞬时参数矩阵Jmafe J = A:^ ? ，其中符号@表示两个矩阵处于相同位置的元素相乘。
2.根据权利要求1所述的方法，其中所述对瞬时频率矩阵fmataix进行归一化处理的公式为:
3.一种基于瞬时幅频参数的探地雷达信号预处理装置，包括: 提取探地雷达信号的瞬时幅度和瞬时频率，分别得到瞬时幅度矩阵Amataix和瞬时频率矩阵fmatrix的模块； 删除Aniatox中的第一行得到疋atrix的模块； 对瞬时频率矩阵fmatHx进行归一化处理，得到归一化瞬时频率矩阵￡atrix的模块； 按照Λ*抓=?/I?，计算联合瞬时参数矩阵Jfflatrlx的模块。
【文档编号】G01V3/38GK104020498SQ201410205408
【公开日】2014年9月3日 申请日期:2014年5月16日 优先权日:2014年5月16日
【发明者】欧阳缮, 张海如, 王国富, 谢跃雷, 张法全, 叶金才 申请人:桂林电子科技大学