Curved-grid kriging gridding method and its application to near-station terrain correction in gravity measurement
TANG Xiao-Ping1,2, LIU Kuan-Hou1, GENG Tao1, FENG Zhi-Han1
1. Xi'an Center of Geology Survey, China Geology Survey, Xi'an 710054, China;
2. Key laboratory for the Study of Focused Magmatism and Giant ore Deposits, MLR, Xi'an 710074, China
In this paper,a high-accuracy and strong stability near-zone terrain correction method is discussed.Firstly,combining Kriging gridding method with curved-grid computation technology,a Curved-Grid Kriging Gridding Method(short for CGKG-M)and its extending edge techniques are implemented to realize data-gridding in circle area;Next,this method is used to compute near-zone terrain correction value and realize the gravity measurement near-zone terrain correction value's computation which contains any angles and circles.Finally,some modeling tests,real data computation and error analysis of this method are carried out.The research result has shown that CGKG-M near-zone terrain correction value computation method has some characteristics,such as high computation accuracy,strong stability and flexibility.It is a simple,efficient and fast method for near-zone terrain correction value computation which matches the measurement data form of GTCS-1 gravity measurement near-station correction instrument.
[1] DZ0004-91重力调查技术规范(1:50000)[S].1991,中华人民共和国国土资源部.[2] DZ/T0171-1997大比列尺重力勘察规范[S].1997,中华人民共和国国土资源部.[3] Hammer S.Terain corrections for gravimeter stations[J].Geophysics,1939,4:184-194.[4] Nagy D.The gravitational attraction of a right rectangular prism[J].Geophysics,1966,31:362-371.[5] Forberg R.Gravity field terrain effect computations by FFT[J].Bull Geod,1985,59:342-360.[6] Li Y C,Sideris M.Improved gravimetric terrain correction[J].Geophysical Journal of International,1994,119:740-752.[7] 冯治汉.区域重力调查中的中区地形改正方法及精度[J].物探与化探,2007,31(5):455-458.[8] 张俊,张宝松,邸兵叶,等.高程数据网格间距对重力中区地形改正精度的影响[J].物探与化探,2014,38 (1):157-161.[9] Schiavone D,Capolongo D,Loddo M.High resolution dems for near station terrain correction in gravimetery[R].EGM2007 International workshop,2007.[10] Schiavone D,Capolongo D,Loddo M.Near-station topographic mases correction for high-accuracy gravimetric [J].Geophysics Prodpecting,2009,57(5):739-751.[11] 刘宽厚,耿涛,杨怀英,等.基于便携激光测距仪的重力测量近区地形改正系统[J].物探与化探,2012,36(3):403-408.[12] Thompson J F,Warsi Z U A,Mastin C W.Numerical Grid Generation: Foundations and Applications[P].North Holland,ISBN:044400985X,1985.[13] Fomberg B.The pseudo spectral method: accurate representation of interfaces in elastic wave calculations[J].Geophysics,1988,53:625-637.[14] Nielsen P,Flemming I,Berg P,et al. Using the pseudo-spectral method on curved grids for 2D elastic forward modeling[J].Geophysical Prospecting,1995,43:369-395.[15] 蒋丽丽,孙建国.基于Poisson方程的曲网格生成技术[J].世界地质,2008,27(3):298-305.[16] 孙建国,蒋丽丽.用于起伏地表条件下的地球物理场数值模拟的正交曲网格生成技术[J].石油地球物理勘探,2009,44(4):494-500.[17] 郭良辉,孟小红,郭志宏,等.地球物理不规则分布数据的空间网格化法[J].物探与化探,2005,29(5):438-442.[18] Matheron G.Principles of geostatistics[J].Economic Geology,1963,58:1246-1266.[19] Matheron G.The intrinsic random functions,and their applications[J].Adv.Appl.Prob.,1973,5:439-468.