[发明专利]一种基于球坐标系的大地电磁三维正演方法

摘要

本发明公开了一种基于球坐标系的大地电磁三维正演方法，属于地球物理勘探技术领域，包括以下步骤：a、建立大地电磁控制方程；b、在球坐标系中，划分成若干个小的倒立四棱柱网格单元；c、设置球坐标模型参数，构建磁场离散的球坐标交错网格单元；d、对单个频率进行循环，将大地电磁控制方程在球坐标交错网格单元进行数值离散；e、求解线性方程组；f、由球坐标张量阻抗公式计算阻抗，再带入卡尼亚视电阻率计算公式求取测点处的视电阻率和相位。本发明能够有效克服地球曲率对大地电磁三维深部探测的干扰，避免因地球曲率所带来的计算误差，适用于大地电磁三维正演数值模拟，能够与地球模型更好的匹配，精确度高。

主权项

1.一种基于球坐标系的大地电磁三维正演方法，其特征在于，包括以下步骤：a、建立大地电磁控制方程；b、在球坐标系中，沿r、θ、三个坐标轴方向，分别用若干平行的球面以不同的间距划分成若干个小的倒立四棱柱网格单元；c、设置球坐标模型参数，包括网格节点坐标、单元、节点编号和单元电阻率，构建磁场离散的球坐标交错网格单元；d、对单个频率进行循环，将步骤a中的大地电磁控制方程在步骤c中的球坐标交错网格单元进行数值离散，得到系数矩阵和右端项，组装到线性方程组中；e、求解线性方程组，迭代过程中开展散度校正，计算得到各个节点的场值；f、根据大地电磁场源，由球坐标张量阻抗公式计算阻抗，再带入卡尼亚视电阻率计算公式求取测点处的视电阻率和相位。