[发明专利]一种探地雷达数据比对算法

摘要

本发明属于地球物理数据解释领域，针对探地雷达数据特点，对两组探地雷达数据，先利用道间配准算法去除丢道和道间距不均匀造成的探地雷达数据形变，再利用道内配准算法去除天气、季节、时间窗等差异造成的探地雷达数据形变，最后利用加窗相关系数比对算法实现探地雷达数据比对，找出两组探地雷达数据的差异范围。

主权项

一种探地雷达数据比对算法，该算法建立在概率论和信息论的理论基础上，特征在于，依次采用道间配准算法和道内配准算法对两组探地雷达数据进行配准，采用相关系数作为道间配准算法的评价参数，采用条件熵作为道内配准算法的评价参数，对配准后的两组探地雷达数据，采用加窗相关系数比较两组探地雷达数据的差异；其中，所述的道间配准算法包括以下步骤：步骤(A1)：待配准探地雷达数据的配准起始道Yc(c1),c1∈Z+，和探地雷达参照数据的配准起始道Yr(r1),r1∈Z+，满足|c1‑r1|＜K,K∈Z+，且c1＜2K,r1＜2K，在待配准探地雷达数据和探地雷达参照数据的Yc(1)～Yc(2K),Yr(1)～Yr(2K)范围内搜索配准起始道；步骤(A2)：将2K分为M′等份，求记作Yc(X1),Yc(X2),...,Yc(XM′)和记作Yr(X1),Yr(X2),...,Yr(XM′)之间的最大相关系数其中p,q∈{1,2,...,M′}，配准起始道Yc(c1),Yr(r1)在Yc(Xp),Yr(Xq)范围内；步骤(A3)：重复步骤(A2)，逐步缩小范围，找到配准起始道Yc(c1),Yr(r1)；步骤(A4)：搜索道数为α,α∈Z+，若某道待配准探地雷达数据Yc(ca),ca＝1,2,...,N′与探地雷达参照数据Yr(rb+1),Yr(rb+2),...,Yr(rb+α),rb+1,rb+2,...,rb+α∈Z+的相关系数最大值满足是阈值，则第rl道是探地雷达参照数据的最相关道，在待配准探地雷达数据的第ca‑1道和第ca道之间线性插值l‑a道探地雷达数据，若将Yc(ca)删除；步骤(A5)：通过找到最不相关待配准探地雷达数据和探地雷达参照数据范围是Yc(Xp′),Yr(Xq′)，阈值是最不相关探地雷达数据相关系数的上确界