The method of filling grid dummy values of potential field based on heat conduction model is proposed in this paper, the process of filling grid dummy in blank areas is regarded as a process of heat conduction to achieve grid data transferring from the known effective information to missing data areas, thus gradually completing data interpolation in blank areas. The numerical experiments show the method of filling grid dummy values by heat conduction model has high precision and is obviously effective. The contours are continuous and smooth on the edges, and also can truly reflect the features and trend of abnormal gradients. The application to the real airborne gravity data indicates that the method is feasible and practical, the effect of filling grid dummy can well conform with the variation law and reflect more geological information of regional gravity field, thus improving the accuracy of data processing and transformation in potential field. The method has practical significance.
[1] Briggs I. Machine contouring using minimum curvature[J].Geophysics,1974,39(1):39-48.[2] Swain C J. A fortran IV program for interpolating irregularly spaced data using the difference equations for minimum curvature[J].Computer & Geoscience,1976,1:231-240.[3] Smith W H F, Wessel P. Gridding with continuous curvature splines in tension[J].Geophysics, 1990,55(3):293-305.[4] Xiong L. Comparison of some gridding methods[J].The Leading Edge,1999:898-900.[5] Michael D O, Richard S S, Marc A V. Gridding aeromagnetic data using longitudinal and transverse horizontal gradiens with the minimum curvature operator[J]. The Leading Edge,2005,142-145.[6] 段本春,徐世浙.磁(重力)异常局部场与区域场分离处理中的扩边方法研究[J].物化探计算技术,1997, 19(4):298-304.[7] 王万银,邱之云,刘金兰,等.位场数据处理中的最小曲率扩边和补空方法研究[J].地球物理学进展, 2009,24(4):1327-1338.[8] 张志厚,徐世浙,余海龙,等.位场向下延拓的迭代法的扩边方法[J].浙江大学学报:工学版,2013,47 (5):918-924.[9] 马国庆,孟令顺,杜小娟,等.磁法数据处理中的扩边和优化中值滤波方法的研究[J].物探化探计算技术, 2010,32(2):194-199.[10] 熊元新,刘涤尘.傅里叶级数的收敛性与吉布斯现象[J].武汉大学学报:工学版,2001,34(1):69-71.[11] Xu S Z, Yang J Y, Yang C F, et al. The iteration method for downward continuation of a potential field from a horizontal plane[J]. Geophysical Prospecting,2007,55:883-890.[12] 徐世浙. 迭代法与FFT法位场向下延拓效果的比较[J].地球物理学报,2007,50(1):285-289.[13] 徐世浙.位场大深度向下延拓[J].物探化探计算技术,2006,28(s1):285-289.[14] 张辉,陈龙伟,任治新,等.位场向下延拓迭代法收敛性分析及稳健向下延拓方法研究[J].地球物理学报,2009,52(4):1107-1113.[15] 刘东甲,洪天求,贾志海,等.位场向下延拓的波数域迭代法及其收敛性[J].地球物理学报,2009,52(6): 1599-1605.[16] Richard S S, Michael D O. Interpolation and gridding of aliased geophysical data using constrained anisotropic diffusion to enhance trends[J]. Geophysics,2005,70(5):121-127.[17] 胡汉平.热传导理论[M].合肥:中国科学技术大学出版社,2010.[18] 李荣华.偏微分方程数值解法[M].北京:高等教育出版社,2010.[19] 陆金甫,关治.偏微分方程数值解法[M].北京:清华大学出版社,2006.[20] 王朔中,克达尔,秦川,等.应用热传导模型的微分方程图像修复[J].上海大学学报:自然科学版,2007, 13(4):331-336.[21] 吕桂霞,马福明.二维热传导方程有限差分区域分解算法[J].数值计算与计算机应用,2006,27(2):96-105.[22] 史策.热传导方程有限差分的MATLAB实现[J].咸阳师范学院学报,2009,24(4):27-29.[23] Dawson C N, Du Q, Dupont T F. A finite difference domain decomposition algorithm for numerical solution of the heat equation[J].Mathematics of Computation,1991,57(195):63-71.[24] Wan Z S, Zhang B L, Chen G N. Design and analysis of finite difference domain decomposition algorithms for the two-dimensional heat equation[C].Washington: IEEE computer Society,2002.[25] 李宛珊,王文洽. 二维热传导方程的有限差分区域分解算法[J].山东大学学报:理学版,2011,46(12): 1-12.[26] 邢家省,李争辉.热传导方程初值问题解的若干性质[J].聊城大学学报:自然科学版,2010,23(3): 6-8.