发明名称 基于单基地测量的地表介质参数反演方法
摘要 基于单基地测量的地表介质参数反演方法,涉及探地雷达。利用单天线垂直入射地面的测量结果对大地表层介质介电常数、电导率进行反演。本发明首先基于矢量网络分析仪测试获得单个喇叭天线垂直照射土壤表面(或金属板)的反射信号S<sub>11</sub>,根据测试结果引入Matrix Pencil(矩阵束)算法,进行不同界面反射信号的辨识分离,估计提取出由土壤表面(或金属板表面)造成的真正的反射信号,从而得到土壤表面的真正的反射系数,利用该方法能有效将各个界面反射分量加以区分,滤除非土壤表面反射信号及其他噪声信号,使得估计精度得到提高,对快速准确获得土壤表层电磁参数具有重要的意义。
申请公布号 CN105467222A 申请公布日期 2016.04.06
申请号 CN201510865939.9 申请日期 2015.12.01
申请人 厦门大学 发明人 刘颜回;俞子伟;张晓瑶;姜美伊;柳清伙;陈忠;刘海;张谅;朱春辉;杨晶
分类号 G01R27/26(2006.01)I;G01S13/88(2006.01)I 主分类号 G01R27/26(2006.01)I
代理机构 厦门南强之路专利事务所(普通合伙) 35200 代理人 马应森
主权项 基于单基地测量的地表介质参数反演方法,其特征在于包括如下步骤:(1)在外场环境中预先将喇叭端口朝上,对准无障碍物的天空,使用矢量网络分析仪测量喇叭端口的反射信号<img file="FDA0000864517300000011.GIF" wi="109" he="79" />若电缆为稳相电缆,则该反射信号可事先在暗室环境中测量;(2)将喇叭垂直对准待测土壤表面,测得反射信号<img file="FDA0000864517300000012.GIF" wi="102" he="79" />(3)在待测土壤表面放置平面金属反射板,用于截获喇叭天线波束辐射的主要能量,测得反射信号<img file="FDA0000864517300000013.GIF" wi="109" he="71" />(4)将<img file="FDA0000864517300000014.GIF" wi="77" he="79" />和<img file="FDA0000864517300000015.GIF" wi="83" he="79" />分别减去喇叭自身反射信号<img file="FDA0000864517300000016.GIF" wi="103" he="78" />得到土壤和金属板的反射信号Y<sub>1</sub>和Y<sub>2</sub>;(5)分别对反射信号Y<sub>1</sub>和Y<sub>2</sub>进行矩阵束处理,滤除其他反射信号及噪声,得到土壤表面和金属板表面的反射信号Y<sub>1</sub><sup>s</sup>和Y<sub>2</sub><sup>s</sup>;(6)根据平面波对理想介质分界的垂直入射、以及反射系数与介质波阻抗的关系,获得地面表层介质的平均波阻抗<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&eta;</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>1</mn><mo>+</mo><msup><mi>R</mi><mi>s</mi></msup></mrow><mrow><mn>1</mn><mo>-</mo><msup><mi>R</mi><mi>s</mi></msup></mrow></mfrac><msub><mi>&eta;</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000864517300000017.GIF" wi="674" he="135" /></maths>其中η<sub>0</sub>为空气波阻抗;另外,弱导电媒介的本征波阻抗可近似表述为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&eta;</mi><mn>1</mn></msub><mo>=</mo><msqrt><mfrac><mi>&mu;</mi><mrow><msub><mi>&epsiv;</mi><mn>0</mn></msub><msub><mi>&epsiv;</mi><mi>r</mi></msub></mrow></mfrac></msqrt><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mi>&omega;</mi><mrow><msub><mi>j&omega;&epsiv;</mi><mn>0</mn></msub><msub><mi>&epsiv;</mi><mi>r</mi></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>&ap;</mo><msqrt><mfrac><mi>&mu;</mi><mrow><msub><mi>&epsiv;</mi><mn>0</mn></msub><msub><mi>&epsiv;</mi><mi>r</mi></msub></mrow></mfrac></msqrt><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>j</mi><mfrac><mi>&omega;</mi><mrow><msub><mi>j&omega;&epsiv;</mi><mn>0</mn></msub><msub><mi>&epsiv;</mi><mi>r</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000864517300000018.GIF" wi="1621" he="183" /></maths>其中μ为自由空间中的磁导率,ε<sub>0</sub>为自由空间中的介电常数,ε<sub>r</sub>为表层土壤的相对介电常数,σ为表层土壤的电导率,ω为角频率,根据式(2)可得到表层土壤的电磁参数估算公式:相对介电常数<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>r</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>&eta;</mi><mn>0</mn></msub><msub><mi>&eta;</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000864517300000019.GIF" wi="826" he="174" /></maths>电导率<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&sigma;</mi><mo>=</mo><mfrac><mrow><mi>i</mi><mi>m</mi><mi>a</mi><mi>g</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mo>*</mo><mi>r</mi><mi>e</mi><mi>a</mi><mi>l</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>*</mo><msub><mi>&omega;&epsiv;</mi><mn>0</mn></msub><msub><mi>&epsiv;</mi><mi>r</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA00008645173000000110.GIF" wi="1023" he="167" /></maths>
地址 361005 福建省厦门市思明南路422号