基于L1L2混合范数的半航空瞬变电磁数据反演方法及装置制造方法及图纸

本发明专利技术公开了基于L1L2混合范数的半航空瞬变电磁数据反演方法，采用地面发射瞬变电磁信号，并利用无人机搭载接收装置获取电磁响应数据，其包括以下步骤：利用接收装置采集待测区域的电磁响应数据；构建待检测区域的地层模型，并预设迭代终止条件；求得电磁响应数据的L2范数约束项对应的横向约束加权矩阵；求得或更新电磁响应数据的L1范数约束项对应的纵向约束矩阵；将数据拟合梯度项、L1范数约束项、L2范数约束项构成目标函数，并求得对地层模型的增量求导等于零后的展开公式；利用展开公式反演计算求得地层模型增量并修正地层模型，达到终止迭代条件后获得反演成像图像。终止迭代条件后获得反演成像图像。终止迭代条件后获得反演成像图像。

【技术实现步骤摘要】
基于L1L2混合范数的半航空瞬变电磁数据反演方法及装置


[0001]本专利技术涉及地球物理航空电磁勘探
，尤其是基于L1L2混合范数的半航空瞬变电磁数据反演方法及装置。

技术介绍

[0002]目前，现有技术中的地球物理电磁勘探多采用航空电磁法，其采用直升机或固定翼飞机搭载发射和观测系统，其无法运用在城市地下空间探测、地面地质调查等场景。另外，现有技术中的航空电磁法采用一维单点反演算法，其在横向上没有约束，数据噪声较强时剖面电阻率连续性较差；在纵向对电阻率约束上通常采用正则化约束，这种正则化约束主要缺点在于使得模型过于光滑，不能有效刻画电性突变界面信息。
[0003]例如专利申请号为“202110820825.8”、名称为“一种半航空瞬变电磁电导率
‑
深度成像方法及设备”的中国专利技术专利，其根据实际情况预先模拟出符合实际情况的电磁响应查询数据库，再将探测到的电磁响应数据在“库”中进行搜索，无需像反演一样进行多次迭代计算，使得半航空瞬变电磁法可以进行快速成像，并快速获得初步成像结果及反演初始模型。
[000本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于L1L2混合范数的半航空瞬变电磁数据反演方法，采用地面发射瞬变电磁信号，并利用无人机搭载接收装置获取电磁响应数据，其特征在于，包括以下步骤：利用接收装置采集待测区域的电磁响应数据；构建待检测区域的地层模型，并预设迭代终止条件；求得电磁响应数据的L2范数约束项对应的横向约束加权矩阵；求得或更新电磁响应数据的L1范数约束项对应的纵向约束矩阵；将数据拟合梯度项、L1范数约束项、L2范数约束项构成目标函数，并求得对地层模型的增量求导等于零后的展开公式；利用展开公式反演计算求得地层模型增量并修正地层模型，达到终止迭代条件后获得反演成像图像。2.根据权利要求1所述的基于L1L2混合范数的半航空瞬变电磁数据反演方法，其特征在于，还包括：采用自适应正则化因子调节L1范数约束项和L2范数约束项在反演过程中的权重。3.根据权利要求1所述的基于L1L2混合范数的半航空瞬变电磁数据反演方法，其特征在于，还包括：获取所述接收装置的基础信息；所述基础信息包括接收线圈面积、飞行高度、抽道时间数和抽道时间；获取发射装置的基础性息；所发射装置的基础性息包括线源长度、电流峰值；获取待测区域内单条测线的测点数、总测线数、任一测点的坐标、偏移距、高程、时间道对应的相应值。4.根据权利要求1所述的基于L1L2混合范数的半航空瞬变电磁数据反演方法，其特征在于，所述迭代终止条件包括最大迭代次数和最小拟合误差。5.根据权利要求1或2或3或4所述的基于L1L2混合范数的半航空瞬变电磁数据反演方法，其特征在于，所述地层模型的参数包括总地层数、总地层厚度、起始地层厚度、任一层的电阻率。6.根据权利要求1或2或3或4所述的基于L1L2混合范数的半航空瞬变电磁数据反演方法，其特征在于，所述横向约束加权矩阵的表达式为：其中，r表示横向约束加权矩阵对角元素，n表示第n个测点p表示n测点的第p个临点。7.根据权利要求6所述的基于L1L2混合范数的半航空瞬变电磁数据反演方法，其特征在于，所述L1范数的正则项的表达式为：
其中，V...

【专利技术属性】
技术研发人员：王绪本，何可，王向鹏，郭明，高文龙，路俊涛，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：