本發(fā)明公開了，包括以下步驟：步驟1、將三維地質(zhì)體剖分成多個四面體單元；分別對觀測點(diǎn)和三維地質(zhì)體建立對應(yīng)的八叉樹結(jié)構(gòu)；步驟2、將三維地質(zhì)體劃分為近區(qū)與遠(yuǎn)區(qū)，步驟3、分別計(jì)算近區(qū)與遠(yuǎn)區(qū)重力參數(shù)；基于近區(qū)所包含的四面體單元的重力積分無奇異性解析表達(dá)式直接計(jì)算近區(qū)的積分，得到近區(qū)重力參數(shù)；并基于觀測點(diǎn)八叉樹和源八叉樹，采用自適應(yīng)快速多極展開算法計(jì)算遠(yuǎn)區(qū)的積分，得到遠(yuǎn)區(qū)重力參數(shù)；步驟4、將步驟3中計(jì)算得到的近區(qū)重力參數(shù)與遠(yuǎn)區(qū)重力參數(shù)求和，得到觀測點(diǎn)重力參數(shù)。本發(fā)明能夠基于非結(jié)構(gòu)化網(wǎng)格剖分，對任意分布的觀測點(diǎn)進(jìn)行重力正演計(jì)算，大幅度地提高了重力正演計(jì)算的效率。