[發明專利]一種基于大地電磁數據反演電阻率和磁化率的方法及系統

有效

• 申請號: 201410258344.2
• 申請日: 2014-06-11
• 公開（公告）號: CN104102814B
• 公開（公告）日: 2017-07-11
• 發明（設計）人: 于昌明;邢寶山
• 申請（專利權）人: 中國科學院地質與地球物理研究所
• 主分類號: G06F19/00
• 分類號: G06F19/00;G01V3/38
• 代理公司: 北京輕創知識產權代理有限公司11212
• 代理人: 楊立

1.一種基于大地電磁數據反演電阻率和磁化率的方法，其特征在于，包括：

步驟1，在TM極化模式下測量地面視電阻率實測值ρ_ai，并基于該視電阻率實測值ρ_ai建立初始地質模型，獲得該初始地質模型各層的電阻率ρ_i和磁化率x_i；

步驟2，在水平層狀介質條件下，基于初始地質模型的電阻率ρ_i和磁化率x_i，采用一維正演方法計算地面的視電阻率理論值ρ_ati；

步驟3，逐層修正初始地質模型的電阻率ρ_i，修正公式為：

${\mathrm{\ρ}}_i^k={\mathrm{\ρ}}_i^{k-1}+\mathrm{\α}\frac{\mathrm{\ξ}|{\mathrm{\ρ}}_{ai}-{\mathrm{\ρ}}_{ati}^{k-1}{|}^2}{1+\mathrm{\ξ}|{\mathrm{\ρ}}_{ai}-{\mathrm{\ρ}}_{ati}^{k-1}|};$

逐層修正初始地質模型的磁化率x_i，修正公式為：

$x_i^k=x_i^{k-1}+\mathrm{\β}\frac{|{\mathrm{\ρ}}_{ai}-{\mathrm{\ρ}}_{ati}^{k-1}|}{{\mathrm{\ρ}}_{ati}^{k-1}(1+\mathrm{\ξ}|{\mathrm{\ρ}}_{ai}-{\mathrm{\ρ}}_{ati}^{k-1}\left|\right)};$

其中，i為地質模型的層數，k為模型修正次數；α和β為設定的步長；ξ為控制電阻率與磁化率所占比例的權參數，其取值范圍為0至∞；

步驟4，基于修正后的各層電阻率和磁化率重新計算地面視電阻率理論值若地面視電阻率實測值ρ_ai和重新計算的地面視電阻率理論值的擬合差小于預先設定的誤差期望值，則停止步驟3中對電阻率和磁化率的修正。