One-dimension forward modeling of MT(python, numerical)

最新推荐文章于 2026-01-04 21:21:17 发布

原创

最新推荐文章于 2026-01-04 21:21:17 发布 · 813 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#python

本文介绍了一维地震电阻率法（MT）的正演模拟，假设地层为均匀各向同性。提供了用Python实现的计算地下介质电导率、相位和阻抗的函数。该函数接受频率、层深和电导率作为输入，并基于经典MT公式进行计算。附带了测试文件和结果展示。

One-dimension forward modeling of MT

Assumptions:

Layer medium
Uniform, homogeneous

Here is Python implementation:

import cmath as cm
import numpy as np
import math
import scipy.linalg as la

def mt1dte(freq, dz, sig):
    mu = 4.0E-7 * math.pi
    ii = cm.sqrt(-1)
    omega = 2.0 * math.pi * freq
    nz = len(sig)
    sig = np.append(sig, sig[nz - 1])
    # dz.append(math.sqrt(2.0 / (sig[nz] * omega * mu)))
    dz = np.append(dz, math.sqrt(2.0 / (sig[nz] * omega * mu)))
    diagA = []
    offdiagA = []

    for ki in range(nz):
        diagA.append(ii*omega*mu*(sig[ki]*dz[ki]+sig[ki+1]*dz[ki+1])-2/dz[ki]-2/dz[ki+1])

    for ki in range(nz-1):
        offdiagA.append(2/dz[ki+1])

    mtxA = np.diag(diagA) + np.diag(offdiagA, k=-1) + np.diag(offdiagA, k=1)
    rhs = np.zeros(nz)
    rhs[0] = - 2.0 / dz