简介： 本发明提供了一种层状地质条件下任意位置偶极源电磁场响应的数值计算方法，包括以下步骤：步骤S100：根据偶极源在层状介质中的位置定义各地层；步骤S200：正交分解偶极源；步骤S300：推导单位偶极源在各个地层产生的洛伦兹矢量势表达式；步骤S400：根据洛伦兹矢量势表达式，推导各地层内电磁场表达式；步骤S500：使用高斯外推法计算电磁场表达式内贝塞尔函数积分，得电磁场的数值解；步骤S600使用数值解乘以偶极源极矩得到偶极源电磁场响应。本发明运用相邻地层之间洛伦兹矢量势所满足的关系式，利用偶极源所在地层特解循环递推得到所有地层洛伦兹矢量势表达式，再通过洛伦兹矢量势得到电磁场表达式。最后使用高斯外推法计算电磁场表达式中贝塞尔函数积分，得到偶极源的电磁场响应。具有理论清晰，计算精度高的特点。