[发明专利]一种基于大地电磁数据反演电阻率和磁化率的方法及系统

摘要

本发明涉及一种基于大地电磁数据反演电阻率和磁化率的方法，所述方法包括步骤1，在TM极化模式下测量地面的视电阻率实测值，并建立初始地质模型，获得初始地质模型各层的电阻率和磁化率；步骤2，在水平层状介质条件下，基于初始地质模型的电阻率和磁化率，采用一维正演方法计算地面的视电阻率理论值；步骤3，逐层修正初始地质模型的电阻率和磁化率；步骤4，基于修正后的电阻率和磁化率，重新计算地面的视电阻率理论值，若视电阻率实测值和重新计算的视电阻率理论值的拟合差小于设定的误差期望值，则停止对电阻率和磁化率的修正。本发明的方法不用进行繁复偏导数矩阵计算，迭代过程稳定收敛，可逐层判断地层的地电阻率参数和磁化率参数。

主权项

一种基于大地电磁数据反演电阻率和磁化率的方法，其特征在于，包括：步骤1，在TM极化模式下测量地面视电阻率实测值ρai，并基于该视电阻率实测值ρai建立初始地质模型，获得该初始地质模型各层的电阻率ρi和磁化率xi；步骤2，在水平层状介质条件下，基于初始地质模型的电阻率ρi和磁化率xi，采用一维正演方法计算地面的视电阻率理论值ρati；步骤3，逐层修正初始地质模型的电阻率ρi，修正公式为：ρik=ρik-1+αξ|ρai-ρatik-1|21+ξ|ρai-ρatik-1|;]]>逐层修正初始地质模型的磁化率xi，修正公式为：xik=xik-1+β|ρai-ρatik-1|ρatik-1(1+ξ|ρai-ρatik-1|);]]>其中，i为地质模型的层数，k为模型修正次数；α和β为设定的步长；ξ为控制电阻率与磁化率所占比例的权参数，其取值范围为0至∞；步骤4，基于修正后的各层电阻率和磁化率重新计算地面视电阻率理论值若地面视电阻率实测值ρai和重新计算的地面视电阻率理论值的拟合差小于预先设定的误差期望值，则停止步骤3中对电阻率和磁化率的修正。