格林样条插值算法及其应用
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
新近发展起来了一种基于双调和算子格林函数计算的数据插值方法———格林样条插值法,这是用中心点位于各观测数据点的多个格林函数进行加权叠加而解析地计算出插值曲面(曲线)的全局插值方法.本文介绍了该方法的基本原理、插值方法和实用程序的发展概况及使用中可能出现的问题.利用格林样条插值法和另两种常用插值法处理了某地布格重力异常的实测数据,对比插值结果说明,格林样条法在抑制虚假异常、稳定显示局部异常和消除奇异边缘效应等方面,具有一定的优越性.
A recently developed interpolator,the greenspline,is based on the application of green's functions of the biharmonic operator.The interpolated surface is constructed by globally summing the weighted Green's functions centered at every data point.The basic mathematic principle of the interpolation,it's developing state and possible problems in application of this technique are reviewed in this paper.The greenspline interpolator and the surface interpolating program in GMT tool were used to test a 2D interpolation of a set of GTOPO30 data.The result indicates that the greenspline acquired a more accurate interpolated topography pattern than the surface did.The greenspline and other two traditional interpolators were also applied to processing the Bouguer gravity anomaly observations obtained somewhere.Comparison of the interpolated results indicates that the greenspline interpolator is superior to the other two in suppressing false anomalies,stably displaying local anomalies and removing singular boundary effect.
引文
[1]陆艳洪,陆文凯,翟正军.一种边缘保持的地震数据插值方法[J].地球物理学报,2012,55(3):991-997.Lu Y H,Lu W K,Zhai Z J.An edge-preserving seismic datainterpolation method[J].Chinese J.Geophys.(in Chinese),2012,55(3):991-997.
    [2]孙章庆,孙建国,韩复兴.复杂地表条件下基于线性插值和窄带技术的地震波走时计算[J].地球物理学报,2009,52(11):2846-2853.Sun Z Q,Sun J G,Han F X.Travel times computation usinglinear interpolation and narrow band technique under complextopographical conditions[J].Chinese J.Geophys.(inChinese),2009,52(11):2846-2853.
    [3]李信富,李小凡.分形插值地震数据重建方法研究[J].地球物理学报,2008,51(4):1196-1201.Li X F,Li X F.Seismic data reconstruction with fractalinterpolation[J].Chinese J.Geophys.(in Chinese),2008,51(4):1196-1201.
    [4]郑伟,许厚泽,钟敏,等.插值公式、相关系数和采样间隔对GRACE Follow-On星间加速度精度的影响[J].地球物理学报,2012,55(3):822-832.Zheng W,Hsu H T,Zhong M,et al.Impacts of interpolationformula,correlation coefficient and sampling interval on theaccuracy of GRACE Follow-On intersatellite range-acceleration[J].Chinese J.Geophys.(in Chinese),2012,55(3):822-832.
    [5]郭良辉,孟小红,郭志宏.位场数据网格化的反插值法[J].地球物理学报,2006,49(4):1183-1190.Guo L H,Meng X H,Guo Z H.Inverse interpolation forgridding of potential field data[J].Chinese J.Geophys.(inChinese),2006,49(4):1183-1190.
    [6]徐世浙,余海龙.位场曲化平的插值-迭代法[J].地球物理学报,2007,50(6):1811-1815.Xu S Z,Yu H L.The interpolation-iteration method forpotential field continuation from undulating surface to plane[J].Chinese J.Geophys.(in Chinese),2007,50(6):1811-1815.
    [7]李晓军,王长虹,朱合华.Kriging插值方法在地层模型生成中的应用[J].岩土力学,2009,30(1):157-162.Li X J,Wang C H,Zhu H H.Kriging interpolation and itsapplication to generating stratum model[J].Rock and SoilMechanics(in Chinese),2009,30(1):157-162.
    [8]夏吉庄.双调和样条内插方法在测井和地震资料整合中的应用[J].石油与天然气地质,2010,31(2):244-249.Xia J Z.An application of biharmonic spline interpolation tointegration of logging and seismic data[J].Oil&Gas Geology(in Chinese),2010,31(2):244-249.
    [9]王家华.克里金地质绘图技术,计算机的模型和算法[M].北京:石油工业出版社,1999.Wang J H.Kriking Geological Mapping technology,ComputerModel and Algorithms(in Chinese)[M].Beijing:PetroliaIndustry Press,1999.
    [10]Schultz C A,Myers S C,Hipp J,et al.Nonstationary Bayesiankriging:apredictive technique to generate spatial correctionsfor seismic detection,location,and identification[J].Bull.Seism.Soc.Am.,1998,88(5):1275-1288.
    [11]潘常周,靳平,肖卫国.利用克里金技术标定新疆及附近地区P/S震相幅值比及其在地震事件识别中的应用[J].地震学报,2007,29(6):625-634.Pan C Z,Jin P,Xiao W G.Calibration of P/S amplituderatios for seismic events in Xinjiang and adjacent areas basedon a Bayesian Kriging method[J].Acta Seismologica Sinica(in Chinese),2007,29(6):625-634.
    [12]涂俐兰,黄丹.插值法在数据修正中的应用[J].数学理论与应用,2012,32(3):112-118.Tu L L,Huang D.Application of interpolation in datacorrection[J].Mathematical Theory and Applications(inChinese),2012,32(3):112-118.
    [13]赵宇,蔺宏伟,鲍虎军.细分曲面拟合的局部渐进插值方法[J].计算机研究与发展,2012,49(8):1699-1707.Zhao Y,Lin H W,Bao H J.Local progressive interpolationfor subdivision surface fitting[J].J.Computer Res.Develop(in Chinese),2012,49(8):1699-1707.
    [14]Sandwell D T.Biharmonic spline interpolation of GEOS-3andSEASAT altimeter data[J].Geophysical Research Letters,1987,14(2):139-142.
    [15]Wessel P.A general-purpose Green’s function-based interpolator[J].Computer&Geosciences,2009,35(6):1247-1254.
    [16]蔺小林.计算方法(高等学校本科规划教材)[M].西安:西安电子科技大学出版社,2009.Lin X L.Method of Calculation,A College Textbook(inChinese)[M].Xi’an:Xi’an Electronic TechnologyUniversity Press,2009.
    [17]清华大学、北京大学《计算方法》编写组.数值计算方法,上册[M].北京:科学出版社,1974.Editorial Group of Method of Calculation.Method ofCalculation,Vol.I(in Chinese)[M].Beijing:SciencePress,1974.
    [18]Roach G F.Green’s Functions[M].Cambridge:CambridgeUniversity Press,1982:344.
    [19]Smith W H F,Wessel P.Gridding with continuous curvaturesplines in tension[J].Geophysics,1990,55(3):293-305.
    [20]MitásováH,Mitás L.Interpolation by regularized splinewith tension:I.Theory and implementation[J].Mathematical Geology,1993,25(6):641-655.
    [21]MitásováH,Hofierka J.Interpolation by regularized splinewith tension:II.Application to terrain modeling and surfacegeometry analysis[J].Mathematical Geology,1993,25(6):657-669.
    [22]Wessel P,Bercovici D.Interpolation with splines in tension:a Green’s function approach[J].Mathematical Geology,1998,30(1):77-93.
    [23]Wessel P,Becker J M.Interpolation using ageneralized Green’sfunction for a spherical surface spline in tension[J].Geophysical Journal International,2008,174(1):21-28.
    [24]Renka R J.Interpolation on the surface of a sphere[J].ACMTransaction of Mathematical Software,1984,10(4):437-439.
    [25]Parker R L.Geophysical Inverse Theory[M].Princeton,NJ:Princeton University Press,1994:386.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心