基于径向基函数的1∶5万规则分布重力数据插值参数优选  

Optimization of interpolation parameters for 1∶50 000 regular distribution gravity data based on radial basis function

作　　者：许海红[1,2] 韩小锋 袁炳强[3] 张春灌 王宝文 赵飞 段瑞锋[4] XU Hai-Hong;HAN Xiao-Feng;YUAN Bing-Qiang;ZHANG Chun-Guan;WANG Bao-Wen;ZHAO Fei;DUAN Rui-Feng(Xi'an Center of Geological Survey,China Geological Survey,Xi'an 710054,China;Key Laboratory of Paleozoic Oil and Gas Geology in North China,Xi'an 710054,China;School of Earth Sciences and Engineering,Xi'an Shiyou University,Xi'an 710065,China;Shaanxi Geo-mining Geophysical and Geochemical Exploration Team Co.Ltd.,Xi'an 710043,China)

机构地区：[1]中国地质调查局西安地质调查中心,陕西西安710054 [2]北方古生界油气地质重点实验室中心,陕西西安710054 [3]西安石油大学地球科学与工程学院,陕西西安710065 [4]陕西地矿物化探队有限公司,陕西西安710043

出　　处：《物探与化探》2021年第6期1539-1552,共14页Geophysical and Geochemical Exploration

基　　金：中国地质调查局地质调查项目“银额盆地西部—北山盆地群油气地质调查”(DD20190092);“银额盆地及周缘油气基础地质调查”(DD20160172)。

摘　　要：为了求取1∶5万规则分布重力数据的最优插值参数,为数据网格化提供定量插值依据,以理论模型的重力数据为例,采用径向基函数法对插值核函数、搜索邻域等插值参数进行优选,根据标准偏差指标评价不同参数对应的插值结果。研究认为:自然三次样条核函数对应的插值精度最高;R^(2)参数处于第一区间(0~1)时插值稳定且精度高;搜索邻域为椭圆形时插值精度高,优选的插值参数分别为:搜索半径R_(1)=3 km、R_(2)=4.5 km,搜索扇区为4个,搜索角度为32°,各向异性比率为0.667,各向异性角度为32°,从所有扇区使用的最大的数据个数为80个,从每个扇区使用的最大的数据个数为20个,所有扇区的最小数据个数(更少则白化节点)为8个,如果空白扇区多于3个则白化节点。In order to select the optimized interpolation parameters of 1∶50 000 regular distribution gravity data to provide quantitative interpolation basis for data gridding.We take the gravity data of the theoretical model as an example,use the radial basis function method to optimize the interpolation parameters,such as the interpolation kernel function and the search neighborhood,and using the standard deviationindex to evaluate the interpolation results corresponding to different parameters.The results indicate that the natural cubic spline kernel function corresponds to the highest interpolation accuracy,when the R^(2) parameter is in the first interval(0~1),the interpolation is stable and accurate.The interpolation accuracy is highest when the search neighborhood is elliptical,and the preferred interpolation parameters are as follows:the search radius R_(1)=3 km,R_(2)=4.5 km,the number of sectors to search is 4,the search angle is 32°,the anisotropy ratio is 0.667,the anisotropy angle is 32°,the maximum number of data to use from all sectors is 80,the maximum number of data to use from each sector is 20,the minimum number of data in all sectors(node is blanked if fewer) is 8,the node is blanked if more than 3 sectors are empty.

关 键 词：径向基函数 SURFER软件 规则分布 1∶5万重力数据 插值参数 

分 类 号：P631[天文地球—地质矿产勘探]