航空瞬变电磁数据和航空大地电磁数据联合反演方法技术

本发明专利技术公开了一种航空瞬变电磁数据和航空大地电磁数据联合反演方法，属于地球物理勘探技术领域，包括以下步骤：a、采用非线性共轭梯度反演法对航空大地电磁数据进行二维反演处理；b、将电阻率取倒数转换为电导率，并进行网格化形成规则的电导率网格数据；c、构建航空瞬变电磁数据一维反演的参考模型；d、获得地下介质的视电导率和视深度；e、构建航空瞬变电磁数据一维反演的初始模型；f、进行约束，完成航空瞬变电磁数据的一维组合正则化约束反演；g、输出航空瞬变电磁和航空大地电磁联合解释的反演结果。本发明专利技术能够为大面积区域的航空地球物理勘探提供精度更高的反演解释结果，以提高综合勘探效率。

Joint inversion method of Aeronautical transient electromagnetic data and aeronautical magnetotelluric data

【技术实现步骤摘要】
航空瞬变电磁数据和航空大地电磁数据联合反演方法
本专利技术涉及到地球物理勘探
，尤其涉及一种航空瞬变电磁数据和航空大地电磁数据联合反演方法。
技术介绍
航空地球物理是利用飞机平台搭载地球物理勘探设备，在空中进行观测的地球物理勘探技术，具有探测面积大、速度快、相对成本低等优点。在戈壁、海域、高山和森林覆盖等地面人员很难进行勘查作业的复杂区域，航空地球物理勘探能取得较好的勘探效果。航空瞬变电磁法(Transientelectromagneticmethod,TEM)是航空地球物理中使用较多，也是发展较快的一种勘探方法，进入21世纪的十几年，航空TEM在观测仪器和数据处理与反演解释方法等方面都取得了巨大的进步，已经广泛应用于矿产资源勘查、地质填图、油气勘查、水资源勘查和环境监测等领域。然而，航空TEM的勘探深度较浅，有效勘探深度在300米左右。航空大地电磁法(Magnetotelluric,MT)是以天然电磁场为场源，通过飞机搭载接收线圈观测电磁场的垂直分量，并以地面基准点的水平分量为参考的一种航空地球物理方法。航空MT具有勘探深度大的优点，能很好的反映出地下深部异常体的电性结构，但是航空MT在浅部的探测分辨率较低，很难反映出地下浅部介质的分布情况。现有的航空TEM勘探深度较低和航空MT浅部探测分辨率较低的问题限制了航空TEM和航空MT数据的有效解释。公开号为CN106338774A，公开日为2017年01月18日的中国专利文献公开了一种基于电导率-深度成像的时间域航空电磁数据反演方法，其特征在于包括如下步骤：(1)根据均匀半空间地电模型的电导率和各时间道的电磁响应值建立均匀半空间地电模型电磁响应与电导率之间的函数关系表；(2)将待解释的地电模型的时间域航空电磁正演模拟数据代入步骤(1)中所述的函数关系表中，获得待解释地电模型各时间道电磁响应对应的视电导率；(3)将视电导率代入深度转换公式中计算各时间道对应的成像视深度；(4)由视电导率和视深度构建初始模型，计算初始模型的理论电磁响应值；(5)将初始模型理论电磁响应值进行反演，并计算初始模型的理论电磁响应值与反演数据之间的误差；(6)判断误差是否满足收敛条件，如果满足则反演结束，否则修正模型参数，判断修正后模型的理论电磁响应值与其反演数据之间的误差是否满足收敛条件，如果满足则反演结束，否则继续修正模型参数，直至达到预期设定的收敛条件；(7)输出反演结果。该专利文献公开的基于电导率-深度成像的时间域航空电磁数据反演方法，勘探深度较浅，不能获得深部地下介质的电性分布结构，无法为航空地球物理勘探提供精度更高的反演解释结果。
技术实现思路
本专利技术为了克服上述现有技术的缺陷，提供一种航空瞬变电磁数据和航空大地电磁数据联合反演方法，本专利技术不仅能够解决航空TEM勘探深度较浅，不能获得地下介质深部电性分布的问题，而且能够解决航空MT在浅部的分辨能力较低的问题，从而为大面积区域的航空地球物理勘探提供精度更高的反演解释结果，以提高综合勘探效率。本专利技术通过下述技术方案实现：航空瞬变电磁数据和航空大地电磁数据联合反演方法，其特征在于，包括以下步骤：a、采用非线性共轭梯度反演法对航空大地电磁数据进行二维反演处理，获得一个二维的电导率剖面结果，反演的目标函数通过式1获得；Φ(m)＝(d-F(m))TV-1(d-F(m))+λmTLTLm式1其中，λ为正则化因子；V为与误差向量e有关的协方差矩阵；L为与模型参数m相关的二维微分矩阵；d为观测数据；F(m)为正演数值模拟值；m为模型参数向量；T为矩阵的转置；b、将二维的电导率剖面结果中的电阻率取倒数转换为电导率，并进行网格化形成规则的电导率网格数据；c、根据电导率网格数据构建航空瞬变电磁数据一维反演的参考模型；d、对航空瞬变电磁数据进行电导率-深度成像处理，获得地下介质的视电导率和视深度；e、根据地下介质的视电导率和视深度构建航空瞬变电磁数据一维反演的初始模型；f、采用步骤c中构建的参考模型和步骤e中构建的初始模型进行约束，完成航空瞬变电磁数据的一维组合正则化约束反演，反演的目标函数通过式2获得；其中，Φd为数据拟合项；为参考模型约束项；为垂直模型粗糙度约束项；λ为总体正则化因子；αr为控制参考模型约束项对应权重的权重因子；αs为垂直模型粗糙度约束项对应权重的权重因子；g、输出航空瞬变电磁和航空大地电磁联合解释的反演结果。所述步骤f中，参考模型约束项是指采用航空大地电磁数据反演结果通过插值构建的参考模型对航空瞬变电磁数据反演中深部电性结构进行约束。所述参考模型约束项通过式3计算获得；Φrm＝[mref-m]TWr[mref-m]式3其中，为参考模型约束项；mref为参考模型对应的模型参数；m为模型参数；T为各层的厚度；Wr为一个Np×Np的对角矩阵；Np为未知模型参数的个数；Np＝2Nl-1；Nl为层状模型的层数。本专利技术的有益效果主要表现在以下方面：1、本专利技术，“a、采用非线性共轭梯度反演法对航空大地电磁数据进行二维反演处理，获得一个二维的电导率剖面结果，反演的目标函数通过式1获得；b、将二维的电导率剖面结果中的电阻率取倒数转换为电导率，并进行网格化形成规则的电导率网格数据；c、根据电导率网格数据构建航空瞬变电磁数据一维反演的参考模型；d、对航空瞬变电磁数据进行电导率-深度成像处理，获得地下介质的视电导率和视深度；e、根据地下介质的视电导率和视深度构建航空瞬变电磁数据一维反演的初始模型；f、采用步骤c中构建的参考模型和步骤e中构建的初始模型进行约束，完成航空瞬变电磁数据的一维组合正则化约束反演，反演的目标函数通过式2获得；g、输出航空瞬变电磁和航空大地电磁联合解释的反演结果”，作为一个完整的技术方案，较现有技术而言，通过将航空瞬变电磁数据和航空大地电磁数据进行联合反演，不仅能够解决航空TEM勘探深度较浅，不能获得地下介质深部电性分布的问题，而且能够解决航空MT在浅部的分辨能力较低的问题，从而为大面积区域的航空地球物理勘探提供精度更高的反演解释结果，以提高综合勘探效率。2、本专利技术，较公开号为CN106338774A，公开日为2017年01月18日的中国专利文献而言，联合了航空大地电磁数据反演，将其反演结果作为航空瞬变电磁数据反演的参考模型，因此能克服航空瞬变电磁勘探深度低的问题，从而获得地下深部介质的电性结构。3、本专利技术，步骤f中，参考模型约束项是指采用航空大地电磁数据反演结果通过插值构建的参考模型对航空瞬变电磁数据反演中深部电性结构进行约束，使得反演结果中不仅包含航空瞬变电磁反演得到的浅部电性信息，同时还保留了航空大地电磁反演得到深部电性信息，利于提高反演解释结果的精度。附图说明下面将结合说明书附图和具体实施方式对本专利技术作进一步的具体说明：图1为本专利技术的流程框图；图2为具体实例中的二维地电模型示意图；图3为单独的航空大地电磁二维非线性共轭梯度反演电阻率剖面结果示意图；图4为单独的航空瞬变电磁一维反演电阻率剖面结果示意图；图5为航空瞬变电磁和航空大地电磁联合反演的电阻率剖面结果示意图。具体实施方式实施例1航空瞬变电磁数据和航空大地电磁数据联合反演方法，包括以下步骤：a、采用非线性共轭梯度反演法对航空大地电磁数据进行二维反演处理，获本文档来自技高网...

【技术保护点】
1.航空瞬变电磁数据和航空大地电磁数据联合反演方法，其特征在于，包括以下步骤：a、采用非线性共轭梯度反演法对航空大地电磁数据进行二维反演处理，获得一个二维的电导率剖面结果，反演的目标函数通过式1获得；Φ(m)＝(d‑F(m))

【技术特征摘要】
1.航空瞬变电磁数据和航空大地电磁数据联合反演方法，其特征在于，包括以下步骤：a、采用非线性共轭梯度反演法对航空大地电磁数据进行二维反演处理，获得一个二维的电导率剖面结果，反演的目标函数通过式1获得；Φ(m)＝(d-F(m))TV-1(d-F(m))+λmTLTLm式1其中，λ为正则化因子；V为与误差向量e有关的协方差矩阵；L为与模型参数m相关的二维微分矩阵；d为观测数据；F(m)为正演数值模拟值；m为模型参数向量；T为矩阵的转置；b、将二维的电导率剖面结果中的电阻率取倒数转换为电导率，并进行网格化形成规则的电导率网格数据；c、根据电导率网格数据构建航空瞬变电磁数据一维反演的参考模型；d、对航空瞬变电磁数据进行电导率-深度成像处理，获得地下介质的视电导率和视深度；e、根据地下介质的视电导率和视深度构建航空瞬变电磁数据一维反演的初始模型；f、采用步骤c中构建的参考模型和步骤e中构建的初始模型进行约束，完成航...

【专利技术属性】
技术研发人员：王绪本，赵宁，余小东，杨聪，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川,51