本發(fā)明提供了一種層狀地質(zhì)條件下任意位置偶極源電磁場(chǎng)響應(yīng)的數(shù)值計(jì)算方法，包括以下步驟：步驟S100：根據(jù)偶極源在層狀介質(zhì)中的位置定義各地層；步驟S200：正交分解偶極源；步驟S300：推導(dǎo)單位偶極源在各個(gè)地層產(chǎn)生的洛倫茲矢量勢(shì)表達(dá)式；步驟S400：根據(jù)洛倫茲矢量勢(shì)表達(dá)式，推導(dǎo)各地層內(nèi)電磁場(chǎng)表達(dá)式；步驟S500：使用高斯外推法計(jì)算電磁場(chǎng)表達(dá)式內(nèi)貝塞爾函數(shù)積分，得電磁場(chǎng)的數(shù)值解；步驟S600使用數(shù)值解乘以偶極源極矩得到偶極源電磁場(chǎng)響應(yīng)。本發(fā)明運(yùn)用相鄰地層之間洛倫茲矢量勢(shì)所滿足的關(guān)系式，利用偶極源所在地層特解循環(huán)遞推得到所有地層洛倫茲矢量勢(shì)表達(dá)式，再通過(guò)洛倫茲矢量勢(shì)得到電磁場(chǎng)表達(dá)式。最后使用高斯外推法計(jì)算電磁場(chǎng)表達(dá)式中貝塞爾函數(shù)積分，得到偶極源的電磁場(chǎng)響應(yīng)。具有理論清晰，計(jì)算精度高的特點(diǎn)。