一维层状模型大地电磁测深 (MT) 和可控源音频大地电磁测深 (CSAMT) 正演计算
TcheL, seistche@gmail.com, USTC
根据（朱仁学，2003），在一维情况下，$ N $ 层模型的地表阻抗可由下式计算：
其中 $ \omega $ 为角频率，$ \mu_0 = 4 \pi \times 10^{ - 7} H/m $ 为空气中的磁导率，$ k_j = \sqrt{ i \omega \mu_0 / \rho_j } $ 为第 j 层的传播常数（或称波数），$ \rho_j $ 为第 j 层的电阻率，$ h_j $ 为第 j 层的厚度。上式仅适用于 $ k_{N - 1} / k_N > 1 $ 的情况；当 $ k_{N - 1} / k_N < 1 $ 时，应采用双曲正切形式的计算公式。
则视电阻率可表示为：
阻抗相位可表示为：
波长定义为：
根据（付海涛 等，2019）及（张健，2011），$ N $ 层水平层状介质表面上水平电偶极子产生的电磁场各分量的表达式为：
其中，$ m_j = \sqrt{\lambda^2 - k_j^2}
上式中关于 Bessel 函数的积分均可由 (Guptasarma and Singh, 1997) 给出的方法进行数值计算得到。
特别地，当 $ N = 1 $ 时，可得到均匀半空间表面电磁场各分量的表达式为：
其中 $ I_j(\cdot) $ 和 $ K_j(\cdot) $ 分别为 j 阶第一和第二类虚宗量 Bessel 函数（即修正/变型 Bessel 函数）。
由柱坐标和直角坐标系之间的转换关系得：
采用不同的方法来定义均匀大地的电阻率或非均匀大地的视电阻率：
其中 $ P_E = I \text{dL} $ 为电偶极距。
- A 型地层：
- H 型地层：
- K 型地层：
- Q 型地层：
- 均匀半空间：
- D 型地层：
- G 型地层：
- H 型地层：
- K 型地层：
- KH 型地层：
- 朱仁学，2003. 大地电磁测深法[M].
- 付海涛，罗维斌，丁志军，余其林，张世宽，2019. 水平电偶极源层状模型垂直磁场全区视电阻率计算方法[J]. 物探与化探, 43(6), 1309-1319. http://doi.org/10.11720/wtyht.2019.0325
- 张健，2011. 可控源音频大地电磁法一维正反演研究[D]. 四川成都：成都理工大学.
- Guptasarma and Singh, 1997. New digital linear filters for Hankel J_0 and J_1 transforms[J]. Geophysical Prospecting, 45, 745-762.