Study on isoparametric finite-element integral algorithm of gravity and magnetic anomaly for body with complex shape
Kim Kangsop1,2, Hu Xiangyun1, Cho Gyonglae2, Nam Myonggun2, Kang Jongnam2, Kim Gyongchil2, Liu Hui1
1. Institute of Geophysics and Geomatics,China University of Geosciences,Wuhan City,Hubei Province,430074,China;
2. Engineering Science Center,Kimchaek University of Technology,Pyongyang,DPRK.
Abstract：Intrinsically,the forward issue of gravity and magnetic exploration is all relative to integral computation.Currently,a great part of the forward computational methods of gravity and magnetic exploration can belong to the analytic integral method,but analytic formula is difficult to be deduced for some complex bodies.The paper presented a new algorithm computing gravity and magnetic anomaly,the principle of which is to use Gaussian numeric integral formula based on isoparametric finite element to directly compute the anomaly on area unit and volume unit and thus to realize integral.The theoretical study and computational result of model showed the effectiveness of the method,which guarantees higher precision without computing logarithmic function and antitrigonometric function and easily implements the program.Especially,the method can use the curve and curved surface to describe geologic body with any complex shapes,providing a new way for computing the gravity and magnetic anomaly of body with any complex shapes.
金钢燮, 胡祥云, 超敬来, 南明根, 康正男, 金敬七, 刘慧. 复杂形体重磁异常的等参数有限元积分算法研究[J]. 石油地球物理勘探, 2009, 44(2): 231-239.
Kim Kangsop, Hu Xiangyun, Cho Gyonglae, Nam Myonggun, Kang Jongnam, Kim Gyongchil, Liu Hui. Study on isoparametric finite-element integral algorithm of gravity and magnetic anomaly for body with complex shape. OGP, 2009, 44(2): 231-239.
