发明名称 一种时间域航空电磁高度校正方法
摘要 本发明公开了一种时间域航空电磁高度校正方法,该方法利用x、z两个分量的感生电动势数据联合计算一种等效均匀半空间电阻率序列,再按等效均匀半空间电阻率为均匀半空间模型的电阻率参数值,在校正高度下正演计算均匀半空间模型的各瞬时感生电动势响应数据,即为最终的校正结果。其中,利用等效均匀半空间电阻率值实现观测感生电动势数据校正是本发明的创新部分。本发明计算速度快,校正结果的相对误差较小,符合航空电磁资料处理要求。
申请公布号 CN101710187B 申请公布日期 2013.01.09
申请号 CN200910263449.6 申请日期 2009.12.17
申请人 成都理工大学;中国地质科学院地球物理地球化学勘查研究所 发明人 毛立峰;胡平;孟庆敏;李文杰;王绪本
分类号 G01V3/38(2006.01)I 主分类号 G01V3/38(2006.01)I
代理机构 代理人
主权项 1.一种时间域航空电磁高度校正方法,其特征在于该方法包括以下步骤:(1)对实测的两感生电动势分量数据进行转换在实际飞行高度下,由两个感生电动势分量实测数据Vx、Vz计算对应各校正时刻T<sub>k</sub>=(t<sub>k</sub>+t<sub>k+L</sub>)/2的幅值比值数据<img file="FSB00000873076300011.GIF" wi="562" he="184" />所述k=1,2,…,n-L,n为采样时间个数,L为时间间隔个数;(2)计算幅度比值数据查询表给定m个离散化查询电阻率数据<img file="FSB00000873076300012.GIF" wi="51" he="58" />(i=1,2,...,m),分别作为等效均匀半空间的电阻率值,先计算它的两分量航空瞬变电磁响应<img file="FSB00000873076300013.GIF" wi="183" he="58" />再计算理论幅度比值查询数据体<maths num="0001"><![CDATA[<math><mrow><msup><mi>A</mi><mn>0</mn></msup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mn>0</mn></msubsup><mo>,</mo><msub><mi>T</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mrow><mo>[</mo><msubsup><mi>V</mi><mi>x</mi><mn>0</mn></msubsup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mn>0</mn></msubsup><mo>,</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>[</mo><msubsup><mi>V</mi><mi>x</mi><mn>0</mn></msubsup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mn>0</mn></msubsup><mo>,</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mi>L</mi></mrow></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></mrow><mrow><msup><mrow><mo>[</mo><msubsup><mi>V</mi><mi>z</mi><mn>0</mn></msubsup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mn>0</mn></msubsup><mo>,</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>[</mo><msubsup><mi>V</mi><mi>z</mi><mn>0</mn></msubsup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mn>0</mn></msubsup><mo>,</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mi>L</mi></mrow></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>并将其作为计算等效均匀半空间电阻率的查询数据表保存在一个二维数组中,所述i=1,2,…,m,k=1,2,…,n-L,n为采样时间个数,L为时间间隔个数;(3)计算等效均匀半空间电阻率值以各校正时刻T<sub>k</sub>对应的两分量幅度比值数据<img file="FSB00000873076300015.GIF" wi="203" he="59" />作为电阻率的单值函数关系,用局部样条插值的方法计算实测两分量数据的幅度比值数据A(T<sub>k</sub>)对应的电阻率值,即为所求的第k个校正时刻的等效均匀半空间电阻率值,逐一计算n-L个校正时刻的等效均匀半空间电阻率值,即完成了等效均匀半空间电阻率的计算,所述i=1,2,…,m;k=1,2,…,n-L,n为采样时间个数,L为时间间隔个数;(4)计算各校正时刻的感生电动势数据根据各校正时刻的等效均匀半空间电阻率值,在给定的校正高度h1下,逐一计算全部n-L个校正时刻对应的各分量感生电动势数据,即为最终的校正结果。
地址 610000 四川省成都市成华区二仙桥东三路1号