A 3D potential on a 3D topography is approximately regarded as a potential on an imaginary horizontal plane. A fast Fourier transform (FFT) is applied to calculate the outward normal derivative of the potential on the horizontal plane. An approximation can be made such that the calculated derivative is used as the outward normal derivative of the potential on the 3D topographic surface. Based on the potential and the approximated normal derivative on the topography, Green's formula is used to obtain the potential at an arbitrary point above the topography. When the potential at a flat level above the topography is obtained, an FFT is used again to determine the potential at other levels above the source of the potential.
A model test shows that the results from this method compare well with analytic solutions. The method has high computation speed and can be used for continuation of 3D potential fields for large data sets, e.g., aeromagnetic data.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址:北京市海淀区学院路29号 邮编:100083 电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700 |