CN107153190A - 用于处理两层介质的多偏移距绕射叠加的成像方法 - Google Patents
用于处理两层介质的多偏移距绕射叠加的成像方法 Download PDFInfo
- Publication number
- CN107153190A CN107153190A CN201710340756.4A CN201710340756A CN107153190A CN 107153190 A CN107153190 A CN 107153190A CN 201710340756 A CN201710340756 A CN 201710340756A CN 107153190 A CN107153190 A CN 107153190A
- Authority
- CN
- China
- Prior art keywords
- wavelets
- source
- mrow
- msup
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Abstract
本发明提供了一种用于处理两层介质的多偏移距绕射叠加的成像方法，包括：将探测区域离散网格化，分成M*N个网格，每个网格作为一个子波源，其中M、N为正整数；对于第i个子波源，从预处理后的标准化探测仪数据集中确定该子波源被每个接收天线接收到的能量，获得该子波源在每帧数据中对应的能量，其中i＝1，2，……，M*N；对于第i个子波源，将该子波源在每帧数据中的能量按照帧数进行叠加，获得该子波源的总能量；以及遍历M*N个子波源，从而获得整个探测区域的地质情况进行成像。通过采用偏移成像中的绕射叠加方法，适用于具有一定高度、多偏移距、多输入多输出的雷达成像要求，且成像清晰，同时对于处理非均质、分层介质具有良好的优势。
Description
技术领域
本发明属于地表和地下探测与成像领域，涉及一种用于处理两层介质的多偏移距绕射叠加的成像方法。
背景技术
月壤结构探测仪LRPR(Lunar Regolith Penetrating Radar)是一种基于嫦娥五号着陆器平台的无载频皮秒脉冲信号体制的高分辨率月壤表面穿透探测雷达，是实现嫦娥五号(CE-5)任务科学探测目标的重要载荷之一。月壤结构探测仪的探测任务是月球次表层结构探测，主要用于如下两个方面：月壤厚度和结构探测；以及为钻取采样过程提供信息支持。
月壤结构探测仪安装在着陆器地板上，距离月面高度约为90cm，由于着陆器系统是静止不动的，故月壤结构探测仪也在静止的状态下工作，采用12个超宽带时域天线组成天线阵列，通过电扫描的方式形成多偏移距的阵列雷达，其天线布局示意图如图1所示，能够实现每一个天线既能发射信号又能接收信号；当其中一个天线发射信号时，其余所有天线均能进行接收。故月壤探测仪具有如下探测特点：
(1)在静止状态下工作；
(2)阵列雷达；
(3)多偏移距，多个天线发射、接收；
(4)距离月面有一个不可忽视的高度，约90cm。
而现有的商业探地雷达通常在移动模式下工作，且只有一个发射天线和一个接收天线，目前大部分的成像技术均是针对商业探地雷达的，对于月壤结构探测仪的成像不适用，因此需要研究一种新的成像算法，来满足在静止状态下、距离地面有一定的高度、多偏移距、多输入多输出的阵列雷达成像需求。
发明内容
(一)要解决的技术问题
本发明提供了用于处理两层介质的多偏移距绕射叠加的成像方法，以至少部分解决以上所提出的技术问题。
(二)技术方案
根据本发明的一个方面，提供了一种用于处理两层介质的多偏移距绕射叠加的成像方法，包括：将探测区域离散网格化，分成M*N个网格，每个网格作为一个子波源，其中M、N为正整数；对于第i个子波源，从预处理后的标准化探测仪数据集中确定该子波源被每个接收天线接收到的能量，获得该子波源在每帧数据中对应的能量，其中i＝1，2，……，M*N；对于第i个子波源，将该子波源在每帧数据中的能量按照帧数进行叠加，获得该子波源的总能量；以及遍历M*N个子波源，从而获得整个探测区域的地质情况进行成像。
在本发明的一个实施例中，确定第i个子波源被每个接收天线接收到的能量的方法，包括：根据天线、折射点和聚焦点的几何位置关系确定第i个子波源和每个接收天线之间的射线传播路径。
在本发明的一个实施例中，两层介质分别为第一介质和第二介质，探测仪发射的信号从第一介质进入第二介质，发生折射，确定第i个子波源和每个接收天线之间的射线传播路径，包括：由入射角、折射角的定义以及斯涅尔定律确定第一介质和第二介质分界面的折射点。
在本发明的一个实施例中，按照如下公式确定所述第一介质和第二介质分界面的折射点：
其中，n为反射系数；s为天线到折射点的矢量在x轴方向上的投影；xf为天线到聚焦点的矢量在x轴方向上的投影；h为天线到折射点的矢量在y轴方向上的投影；d为折射点到聚焦点的矢量在y轴方向上的投影；ε1为第一介质的相对介电常数；ε2为第二介质的相对介电常数。
在本发明的一个实施例中，网格为正方形网格，每个网格的边长满足如下条件：小于1/4波长；波长的大小等于c/f，c为光速，f为获取探测仪数据集的探测频率。
在本发明的一个实施例中，获取探测仪数据集的探测仪进行探测的天线阵列为包含j个天线的多输入、多输出天线阵列，其中每个天线能够接收其余天线发射的信号，每个天线接收的另外一个发射天线的能量对应一帧数据，按照如下公式将每个子波源在每帧数据中的能量按照帧数进行叠加：
其中，I(x，y)为子波源(x，y)按照帧数叠加后的振幅值；Am(*，x，y)表示子波源(x，y)在第m帧数据中的振幅值；Z为多输入、多输出天线阵列获取的数据帧数总量；Ra，t和Ra，r分别为分界面两个折射点到发射天线和接收天线的距离；Rg，t和Rg，r分别为子波源(x，y)到两个折射点的距离；c为光速；vg为电磁波在第二介质中的速度。
在本发明的一个实施例中，多输入、多输出天线阵列获取的数据帧数总量Z由下式确定：
Z＝j(j-1)
其中，j为多输入、多输出天线阵列包含的天线总数。
在本发明的一个实施例中，在将该子波源在每帧数据中的能量按照帧数进行叠加，获得该子波源的总能量的步骤之后还包括：根据叠加后的振幅值判断每个子波源是否为真实反射点或反射面，从而确定该子波源对应的探测区域的地质情况进行成像。
在本发明的一个实施例中，判断每个子波源是否为真实反射点或反射面的方法为：当子波源是一个真实的反射点或反射面时，叠加后的振幅值增大；反之，叠加后的振幅值减小。
在本发明的一个实施例中，探测区域为月球表面，其探测范围的宽*高为：1.6m*3m。
(三)有益效果
从上述技术方案可以看出，本发明提供的用于处理两层介质的多偏移距绕射叠加的成像方法，至少具有以下有益效果之一：
(1)通过采用偏移成像中的绕射叠加方法，将探测区域离散网格化，通过确定每个网格作为子波源辐射的能量在天线阵列获得的每帧数据中的位置以及将子波源在每帧数据的能量进行叠加，最终获取每个子波源的总能量，得到整个探测区域内所有子波源的能量分布，从而确定其地质情况进行成像，适用于具有一定高度、多偏移距、多输入多输出的雷达成像要求，且成像清晰；
(2)绕射叠加方法对于处理非均质、分层介质具有良好的优势，由于月壤结构探测仪距离月球表面具有约90cm高度的真空区域，和被探测的月壤的相对介电常数是不同的，构成了分层介质，本发明提供的方法在确定折射点以及射线传播路径方面具有便捷的优势。
附图说明
图1为现有技术月壤结构探测仪的天线布局示意图。
图2为根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法流程图。
图3A为根据本发明实施例基于绕射叠加确定第i个子波源被每个接收天线接收到能量的射线传播路径示意图；图3B为根据本发明实施例基于绕射叠加确定第i个子波源被每个接收天线接收到能量的光路几何示意图。
图4为根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法进行数值模拟实验得到的结果图，其中，(a)为数值模拟试验的模型；(b)为成像结果。
图5为根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法进行试验场地验证得到的结果图，其中(a)为试验场地及预埋目标设置；(b)为成像结果。
【符号说明】
矢量(加粗)：Ra，t，Ra，r，Rg，t，Rg，r；
距离：Ra，t，Ra，r，Rg，t，Rg，r；
位置点字母：R1，R2，R；
长度符号：xf，s，d，h；
角度符号：θi，θt；
坐标符号：(x，y)
反射系数：n；
相对介电常数：ε1，ε2。
具体实施方式
本发明提供了一种用于处理两层介质的多偏移距绕射叠加的成像方法，通过采用偏移成像中的绕射叠加方法，将探测区域离散网格化，通过确定每个网格作为子波源辐射的能量在天线阵列获得的每帧数据中的位置以及将子波源在每帧数据的能量进行叠加，最终获取每个子波源的总能量，得到整个探测区域内所有子波源的能量分布，从而确定其地质情况进行成像，适用于具有一定高度、多偏移距、多输入多输出的雷达成像要求；同时对于处理非均质、分层介质具有良好的优势。
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步详细说明。
图1为现有技术月壤结构探测仪的天线布局示意图。参照图1所示，月壤探测仪的成像过程如下：在雷达控制器的控制下，选通12个超宽带时域天线组成的天线阵列中的其中一个天线作为发射天线，其余所有天线均作为接收天线进行回波接收，当遍历所有天线均作为发射天线进行探测之后，月壤结构探测仪便能获得天线下方区域月壤的一次完整的探测数据，于是获得12*11帧数据，即132帧完整的探测数据，通过对探测数据进行处理和成像，便能获取一幅探测区域内月壤分层结构图像。
在本发明的一个示例性实施例中，提供了用于处理两层介质的多偏移距绕射叠加的成像方法。其中，两层介质分别为第一介质和第二介质，雷达信号从第一介质进入第二介质，发生折射，在本实施例中，第一介质为空气，第二介质为月壤。
图2为根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法流程图。由图2可知，用于处理两层介质的多偏移距绕射叠加的成像方法，包括：
步骤S202：将探测区域离散网格化，分成M*N个网格，每个网格作为一个子波源；
根据惠更斯原理：球形波面上的每一点或面源都是一个次级球面波的子波源，子波的波速和频率等于初级波的波速和频率，此后每一时刻的子波波面的包络就是该时刻总的波动的波面，据此，将M*N个网格中的每个网格视为一个子波源；
其中，M、N为正整数，数值越大，则精度越高，具体M、N值由探测的范围和离散的网格精度确定；
因为探测精度的需求，每个网格的边长至少应小于四分之一个波长；该波长大小等于c/f，c为光速，c＝3*108m/s，f为探测频率，在本实施例中等于月壤结构探测仪的中心频率2G，其值为：2G＝2*109Hz，则每个网格的边长应小于((3*108)/(2*109))/4＝0.375m；
本实施例针对月壤结构探测仪来说，其探测范围的宽*深为1.6m*3m，故M、N取值满足：M≥43，N≥80；
步骤S204：对于第i个子波源，从预处理后的标准化探测仪数据集中确定该子波源被每个接收天线接收到的能量，获得该子波源在每帧数据中对应的能量，其中i＝1，2，……，M*N；
由于月壤结构探测仪距离月球表面具有约90cm高度的真空区域，和被探测的月壤的相对介电常数是不同的，构成了分层介质，因此射线在两种分界面必然会产生折射，因此求出分界面的折射点是确定射线传播路径的关键所在，这对于确定第i个子波源被每个接收天线接收到的能量也是至关重要的；
图3A为根据本发明实施例基于绕射叠加确定第i个子波源被每个接收天线接收到能量的射线传播路径示意图；图3B为根据本发明实施例基于绕射叠加确定第i个子波源被每个接收天线接收到能量的光路几何示意图，参照图3A可知，从发射天线发出的射线在介电常数为ε1的空气中传播，沿着矢量Ra，t到达位置点R1，发生折射，然后在介电常数为ε2的月壤中传播，对应矢量Rg，t，到达子波源(x，y)，然后在该子波源处发生反射沿着矢量Rg，r在月壤中前进，到达月壤与空气的界面处，在R2处发生折射，进入空气，在空气中沿着矢量Ra，r前进，到达接收天线；由上可知：确定射线传播路径，即需要求出Ra，t、Ra，r、Rg，t和Rg，r，而求出这几个矢量的关键在于确定折射点R1，R2；那么由图3B可知，根据天线、折射点和聚焦点的几何位置关系，由入射角、折射角的定义以及斯涅尔定律可知：
其中，θi为入射角；θt为折射角；s为天线到折射点的矢量在x轴方向上的投影；h为天线到折射点的矢量在y轴方向上的投影；xf为天线到聚焦点的矢量在x轴方向上的投影；d为折射点到聚焦点的矢量在y轴方向上的投影；n为反射系数；
然后结合公式(1)～(3)，得到关于s的一元四次方程，表达式如下：
通过求解如公式(4)所示关于s的一元四次方程，便可以求得s，即获得折射点的位置，从而确定每个子波源在每帧数据中对应的能量；
步骤S206：对于第i个子波源，将该子波源在每帧数据中的能量按照帧数进行叠加，获得该子波源的总能量；以及遍历M*N个子波源，从而获得整个探测区域的地质情况进行成像。
每个子波源的位置坐标为(x，y)，将每个子波源在每帧数据中的能量按照帧数进行叠加，可以获得在132帧数据中的总能量，能量叠加的表达式如下：
其中，Am(*，x，y)表示子波源(x，y)在第m帧数据中的振幅值；Ra，t和Ra，r分别为分界面两个折射点到发射天线和接收天线的距离；Rg，t和Rg，r分别为子波源(x，y)到两个折射点的距离；c为光速；vg为电磁波在月壤中的速度；
当该子波源是一个真实的反射点或反射面时，由于各帧数据中的振幅接近同相叠加，则叠加后的振幅值增大；反之，不是真实的反射点或反射面时，由于各帧数据的振幅为非同相叠加，则叠加后的振幅值自然会减小；由此可以判断每个子波源是否为真实反射点或反射面；然后采用同样的方式对M*N个子波源都获取其总能量并进行判定，从而获得整个探测区域的地质情况进行成像。
需要特别说明的是，以上实施例中提供的用于处理两层介质的多偏移距绕射叠加的成像方法除了可以满足月壤结构探测仪的成像需求之外，还可以用于其他分层介质探测的雷达成像，只要满足：静止状态成像、距离地面有一定的高度、多偏移距、多输入多输出的阵列雷达成像特点即可。
为了验证本发明的可行性和有益效果，根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法分别进行了数值模拟验证和地面验证试验。
图4为根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法进行数值模拟实验得到的结果图，其中，(a)为数值模拟试验的模型；(b)为成像结果。
参照图4中(a)图，数值模拟区域的深*宽的范围为3m*1.6m，天线的布局情况与月壤结构探测仪完全一致，一共设置有五个目标体，包括：两个石块，对应图4(a)中的点，深度分别为0.8米和2.2米；三个石板，对应图4(a)中的线条，深度分别为1米、1.5米和2米；数值模拟后成像结果如图4中(b)图所示，可以清晰地看出五个目标体，且可以识别出上下目标体的上下界面。
图5为根据本发明实施例所示用于处理两层介质的多偏移距绕射叠加的成像方法进行试验场地验证得到的结果图，其中(a)为试验场地及预埋目标设置；(b)为成像结果。
参照图5中(a)图，试验场地为一个长*宽*深尺寸为7m*3m*2.5m的池子，在池子底部铺设了金属板，其中方框里面的部分即为月壤结构探测仪探测时所覆盖的范围，其中包括4个目标体：深度为0.5m的玄武岩块2，深度为0.8m的金属球，深度为1m的聚四氟乙烯板以及深度为2m的花岗岩石板；探测成像后的结果如图5中(b)图所示，其中，椭圆框为四个已知目标体成像的结果，从中可以看出，可以清晰的识别四个目标体。
综上所述，本发明提供了一种用于处理两层介质的多偏移距绕射叠加的成像方法，通过采用偏移成像中的绕射叠加方法，将探测区域离散网格化，确定每个网格子波源辐射的能量在每帧数据中的位置以及将子波源在每帧数据的能量进行叠加，最终获取每个子波源的总能量，得到整个探测区域内所有子波源的能量分布，从而确定其地质情况进行成像，适用于具有一定高度、多偏移距、多输入多输出的雷达成像要求，并且成像清晰；同时在确定分层介质的折射点以及射线传播路径方面也具有便捷的优势。
当然，根据实际需要，本发明提供的用于处理两层介质的多偏移距绕射叠加的成像方法，还包含其他的常用方法和步骤，由于同发明的创新之处无关，此处不再赘述。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种用于处理两层介质的多偏移距绕射叠加的成像方法，包括：
将探测区域离散网格化，分成M*N个网格，每个网格作为一个子波源，其中M、N为正整数；
对于第i个子波源，从预处理后的标准化探测仪数据集中确定该子波源被每个接收天线接收到的能量，获得该子波源在每帧数据中对应的能量，其中i＝1，2，……，M*N；
对于第i个子波源，将该子波源在每帧数据中的能量按照帧数进行叠加，获得该子波源的总能量；以及
遍历M*N个子波源，从而获得整个探测区域的地质情况进行成像。
2.根据权利要求1所述的成像方法，其中，确定所述第i个子波源被每个接收天线接收到的能量的方法，包括：
根据天线、折射点和聚焦点的几何位置关系确定第i个子波源和每个接收天线之间的射线传播路径。
3.根据权利要求2所述的成像方法，其中，所述两层介质分别为第一介质和第二介质，探测仪发射的信号从第一介质进入第二介质，发生折射，确定第i个子波源和每个接收天线之间的射线传播路径，包括：
由入射角、折射角的定义以及斯涅尔定律确定第一介质和第二介质分界面的折射点。
4.根据权利要求3所述的成像方法，按照如下公式确定所述第一介质和第二介质分界面的折射点：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<mi>s</mi>
<mn>4</mn>
</msup>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>x</mi>
<mi>f</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>x</mi>
<mi>f</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mo>+</mo>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>x</mi>
<mi>f</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mfrac>
<msqrt>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
</msqrt>
<msqrt>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
</msqrt>
</mfrac>
</mrow>
其中，n为反射系数；s为天线到折射点的矢量在x轴方向上的投影；xf为天线到聚焦点的矢量在x轴方向上的投影；h为天线到折射点的矢量在y轴方向上的投影；d为折射点到聚焦点的矢量在y轴方向上的投影；ε1为第一介质的相对介电常数；ε2为第二介质的相对介电常数。
5.根据权利要求1所述的成像方法，其中，所述网格为正方形网格，每个网格的边长满足如下条件：小于1/4波长；
所述波长的大小等于c/f，c为光速，f为获取探测仪数据集的探测频率。
6.根据权利要求1所述的成像方法，其中，获取探测仪数据集的探测仪进行探测的天线阵列为包含j个天线的多输入、多输出天线阵列，其中每个天线能够接收其余天线发射的信号，每个天线接收的另外一个发射天线的能量对应一帧数据，按照如下公式将每个子波源在每帧数据中的能量按照帧数进行叠加：
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Z</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>g</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>g</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>v</mi>
<mi>g</mi>
</msub>
</mfrac>
<mo>,</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
其中，I(x，y)为子波源(x，y)按照帧数叠加后的振幅值；Am(*，x，y)表示子波源(x，y)在第m帧数据中的振幅值；Z为多输入、多输出天线阵列获取的数据帧数总量；Ra，t和Ra，r分别为分界面两个折射点到发射天线和接收天线的距离；Rg，t和Rg，r分别为子波源(x，y)到两个折射点的距离；c为光速；vg为电磁波在第二介质中的速度。
7.根据权利要求6所述的方法，其中，多输入、多输出天线阵列获取的数据帧数总量Z由下式确定：
Z＝j(j-1)
其中，j为多输入、多输出天线阵列包含的天线总数。
8.根据权利要求7所述的成像方法，在将该子波源在每帧数据中的能量按照帧数进行叠加，获得该子波源的总能量的步骤之后还包括：根据叠加后的振幅值判断每个子波源是否为真实反射点或反射面，从而确定该子波源对应的探测区域的地质情况进行成像。
9.根据权利要求8所述的成像方法，其中，所述判断每个子波源是否为真实反射点或反射面的方法为：当子波源是一个真实的反射点或反射面时，叠加后的振幅值增大；反之，叠加后的振幅值减小。
10.根据权利要求1至9任一项所述的成像方法，所述探测区域为月球表面，其探测范围的宽*高为：1.6m*3m。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710340756.4A CN107153190B (zh) | 2017-05-12 | 2017-05-12 | 用于处理两层介质的多偏移距绕射叠加的成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710340756.4A CN107153190B (zh) | 2017-05-12 | 2017-05-12 | 用于处理两层介质的多偏移距绕射叠加的成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107153190A true CN107153190A (zh) | 2017-09-12 |
| CN107153190B CN107153190B (zh) | 2020-06-26 |
Family
ID=59793912
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710340756.4A Active CN107153190B (zh) | 2017-05-12 | 2017-05-12 | 用于处理两层介质的多偏移距绕射叠加的成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107153190B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706583A (zh) * | 2009-10-16 | 2010-05-12 | 西安交通大学 | 一种多偏移距vsp成像的局域化相空间方法 |
| CN101900833A (zh) * | 2010-06-02 | 2010-12-01 | 西安石油大学 | 一种地震散射p-p波成像速度分析方法 |
| CN102495426A (zh) * | 2011-12-02 | 2012-06-13 | 中国科学院地质与地球物理研究所 | 一种克希霍夫积分地震成像方法 |
| CN102520443A (zh) * | 2011-12-13 | 2012-06-27 | 中国科学院地质与地球物理研究所 | 一种绕射波叠前偏移的方法 |
| CN103439694A (zh) * | 2013-08-29 | 2013-12-11 | 中国科学院国家天文台 | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 |
| CN106019272A (zh) * | 2016-06-13 | 2016-10-12 | 中国科学院电子学研究所 | 一种采用超宽带雷达技术探测月壤分层结构的系统及方法 |
| US20170010351A1 (en) * | 2014-12-19 | 2017-01-12 | Institute Of Electronics, Chinese Academy Of Sciences | Method and system for detecting geological structure of an extraterrestrial solid planet by using a single-transmitter and multiple-receiver radar |
-
2017
- 2017-05-12 CN CN201710340756.4A patent/CN107153190B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706583A (zh) * | 2009-10-16 | 2010-05-12 | 西安交通大学 | 一种多偏移距vsp成像的局域化相空间方法 |
| CN101900833A (zh) * | 2010-06-02 | 2010-12-01 | 西安石油大学 | 一种地震散射p-p波成像速度分析方法 |
| CN102495426A (zh) * | 2011-12-02 | 2012-06-13 | 中国科学院地质与地球物理研究所 | 一种克希霍夫积分地震成像方法 |
| CN102520443A (zh) * | 2011-12-13 | 2012-06-27 | 中国科学院地质与地球物理研究所 | 一种绕射波叠前偏移的方法 |
| CN103439694A (zh) * | 2013-08-29 | 2013-12-11 | 中国科学院国家天文台 | 一种验证探地雷达在月球次表层探测深度和分辨率的方法 |
| US20170010351A1 (en) * | 2014-12-19 | 2017-01-12 | Institute Of Electronics, Chinese Academy Of Sciences | Method and system for detecting geological structure of an extraterrestrial solid planet by using a single-transmitter and multiple-receiver radar |
| CN106019272A (zh) * | 2016-06-13 | 2016-10-12 | 中国科学院电子学研究所 | 一种采用超宽带雷达技术探测月壤分层结构的系统及方法 |
Non-Patent Citations (2)
| Title |
|---|
| X. FENG ,ET AL: "Pre-stack migration applied to GPR for landmine detection", 《INVERSE PROBLEMS》 * |
| 王正成: "隧道质量无损检测的GPR新技术研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107153190B (zh) | 2020-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102313885B (zh) | 多维度微波成像系统及方法 | |
| CN100365434C (zh) | 激光雷达系统和方法 | |
| US6590519B2 (en) | Method and system for identification of subterranean objects | |
| Barnes et al. | Geological analysis of Martian rover‐derived digital outcrop models using the 3‐D visualization tool, Planetary Robotics 3‐D Viewer—Pro3D | |
| RU2265866C1 (ru) | Способ повышения радиолокационного разрешения, система для его осуществления и способ дистанционного выявления системой малоразмерных объектов | |
| CN101539627B (zh) | 一种电离层立体探测星载sar成像处理平台的构建方法 | |
| CN102419453A (zh) | 长导线源瞬变电磁地空探测方法 | |
| CN104237956A (zh) | 电性源瞬变电磁地空探测方法 | |
| CN102830401B (zh) | 一种探地雷达加窗加权后向投影成像方法 | |
| CN102053247A (zh) | 一种多基线合成孔径雷达三维成像的相位校正方法 | |
| WO2016095208A1 (zh) | 采用一发多收雷达探测地外固态星体地质结构的方法与系统 | |
| CN102445690B (zh) | 合成孔径雷达三维成像的qr分解方法 | |
| CN103389497B (zh) | 机载稀疏阵列天线下视三维成像雷达系统和成像方法 | |
| CN109283590B (zh) | 基于小波变换的多源重力数据融合方法 | |
| WO2016168412A1 (en) | Synthetic aperture radar mineral prospector | |
| CN106772284A (zh) | 考虑近海面大气环境影响的编队内雷达远场预测方法 | |
| CN105824019B (zh) | 大型分布式空间监测雷达的最优化波束对齐方法 | |
| CN107153190A (zh) | 用于处理两层介质的多偏移距绕射叠加的成像方法 | |
| Zhang et al. | Performance comparison of compact phased arrays and traditional seismic networks for microseismic monitoring at a CO2 sequestration test site | |
| CN110297237B (zh) | 考虑天线方向图的探地雷达绕射叠加成像方法及系统 | |
| CN108120955A (zh) | 一种用于长基线水声定位系统性能评估的方法 | |
| Zhang et al. | A detection performance analysis of Sanya incoherent scatter radar tristatic system | |
| CN102478663A (zh) | 一种三维地震观测系统偏移噪声获取方法及装置 | |
| Le Gall et al. | An imaging HF GPR using stationary antennas: Experimental validation over the Antarctic ice sheet | |
| Arvidson et al. | Rocky 7 prototype Mars rover field geology experiments: 1. Lavic Lake and Sunshine Volcanic Field, California |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |