[发明专利]一种地质雷达三维数据偏移归位方法

权利要求书

1.一种地质雷达三维数据偏移归位方法，该方法以地质雷达二维剖面为基础，首先通过边界延拓和插值获得稠密剖面，然后通过预估波速的频率-波数域偏移方法对稠密剖面二维成像，最后通过像间插值构建三维数据，具体步骤如下：

步骤1：对地质雷达二维剖面D(x,y,t)，其中x＝1,2,...,m，y＝0，t＝1,2,...,n，m为二维剖面天线移动方向的道数，n为二维剖面电磁波传播方向的时间域采样点数，对二维剖面进行边界延拓得到

其中x′＝1,2,...,m,m+1,...,m+N，N为边界延拓道数，再进行道间三阶插值

其中，x′为插值后稠密剖面的道序号，d_l为与的距离，l＝1,2,...,N，N为边界延拓道数，通常N＝6；

步骤2：稠密剖面其中x′＝1,2,...,m,m+1,...,2m，y＝0，t＝1,2,...,n采用预估波速的频率-波数域偏移方法，首先经过傅里叶变换，得到

其中ω为角频率，k_x为天线移动方向的波数域分量，再经过傅里叶逆变换，得到

其中，v_x′y为剖面位置(x′,y)处的局部波速，k_y为电磁波传播方向的波数域分量，x′＝1,2,...m,m,+1,..m.，,2y＝1,2,...,n，t＝0；

步骤3：对成像结果进行三维加权插值，由于t＝0，可以简化为由多组二维成像结果构成的三维数据可以表示为其中z＝1,2,...,s，s为二维成像结果的总数，此时三维数据的任意一点可以表示为

其中，Δx∈(x′-1,x′+1)为距离x′小于1的某个值，Δy∈(y-1,y+1)为距离y小于1的某个值，Δz∈(z-1,z+1)为距离z小于1的某个值，α_i为第i个权系数值，β_i为第i个权系数值，γ_i为第i个权系数值，其中i＝0,1,2，x′＝1,2,...,m,m+1,...,2m，y＝1,2,...,n，z＝1,2,...,s。