[发明专利]基于波恩迭代法的航空瞬变电磁电导率三维反演方法

权利要求书

1.基于波恩迭代法的航空瞬变电磁电导率三维反演方法，其特征在于包括以下步骤：

1)数据预处理；

2)三维初始模型设定；

3)模型参数的更新；

4)二次场值的计算；

5)误差计算；

6)模型更新量的计算；

7)收敛条件判断。

2.如权利要求1所述基于波恩迭代法的航空瞬变电磁电导率三维反演方法，其特征在于在步骤1)中，所述数据预处理的具体方法为：在航空电磁法中，接收机在空中接收的电压信号为总场响应信号，并且接收过程中会受到一定的噪声干扰；通过对接收机接收数据进行平均，滤波和去噪处理，获取较为理想的总场响应信号；再将计算的理论均匀半空间响应从总响应信号中扣除，得到反演所需要的二次场响应信号。

3.如权利要求1所述基于波恩迭代法的航空瞬变电磁电导率三维反演方法，其特征在于在步骤2)中，所述三维初始模型设定的具体方法为：根据具体实验飞行的信息建立三维分层初始模型，初始模型的参数包括大地的分层信息，每层的分界面位置，每层介质的介电常数、电导率、磁导率信息以及异常体的位置和电导率、介电常数。

4.如权利要求1所述基于波恩迭代法的航空瞬变电磁电导率三维反演方法，其特征在于在步骤3)中，所述模型参数的更新的具体方法为：在电导率的三维反演算法中，主要是对计算得到的二次场值进行拟合，其中二次场的数值是通过地下介质和异常体的对比度函数体现的，对比度函数χ(r)的定义如下：

χ(r)=ϵ~(r)ϵ~q-1]]>

其中表示复介电常数，定义为r表示空间位置，j表示虚数单位，ω表示角频率；在每一次迭代计算中，通过求解成本函数得到模型参数的更新量δχ_n+1(r)，χ_n+1(r)即为对比度函数在空间的分布；计算得到第n+1次的更新量之后，通过上一次的对比度函数的分布，计算得到下一次迭代时需要带入的对比度函数，即χ_n+1(r)＝χ_n(r)+δχ_n+1(r)；χ_n(r)为第n次迭代时对比度函数的分布，χ_n+1(r)为第n+1迭代时对比度函数的分布；在一般反演情况下，通常采用均匀半空间为初始模型，即第一次迭代时对比度函数为0。