基于变差函数的插值方法计算华北地区重力场变化

摘要/Abstract

摘要： 定点重复重力测量是获取区域重力场变化的主要手段之一, 重力场特征与地形起伏、构造走向等因素相关。以华北地区为例, 考虑区域内地形和构造的北东向分布规律, 从EGM2008重力模型中拟合各向异性变差函数参数, 利用变差函数网格化插值, 对华北地区2009—2013年期间重力场观测数据进行网格化重建, 获取华北地区重力场时空变化结果。研究结果表明基于区域地形、构造特征的各向异性变差函数插值方法, 获得的空间重力场变化在重力异常梯级带上更加明显, 重力变化与活动构造分布具有更好的一致性。本文研究方法对于恢复区域时空重力场异常具有重要意义, 有助于提高应用重力资料划分潜在地震危险区空间位置的精度, 为华北地区震情研判和构建地震预报定量指标体系提供可靠的地球物理场数据。

关键词: 重力场变化, 区域重力异常, 流动重力, 变差函数, Kriging插值, 华北地区

Abstract: Fixed-point gravity repeat measurement is one of the main approaches to get regional gravity changes. Regional gravity field is related to topographic relief, structural faults direction and other factors. We took north China as an example, considering regularity of topographic relief and structural faults with the north-east direction (NE), then fitted the anisotropic parameters of variogram based on EGM2008 gravity model. After that,we rebuilt the gravity variation field with the observations from 2009 to 2013 in north China by anisotropic variogram in Kriging gridding, for further studying temporal and spatial gravity variations in northern China. The results show that there are obvious temporal and spatial gravity variations, especially in gravity anomaly gradient belt through the variogram approach, which is based on regional topography and structure characteristics. Besides, the gravity anomaly field had a better consistency with the distribution of tectonic faults, which could provide the geophysical basis for seismic risk analysis in north China.

Key words: Gravity variation, Regional gravity anomaly, Mobile gravity measurement, North China

中图分类号:

P315.7

徐伟民, 孙少波, 李晓一, 卢红艳, 郑秋月. 基于变差函数的插值方法计算华北地区重力场变化[J]. 地震, 2016, 36(4): 171-185.

XU Wei-min, SUN Shao-bo, LI Xiao-yi, LU Hong-yan, ZHENG Qiu-yue. Calculation of Gravity Variations in North China based on Kriging Gridding with Variogram Approach[J]. EARTHQUAKE, 2016, 36(4): 171-185.

参考文献

[1] 孙和平. 重力场的时间变化与地球动力学[J]. 中国科学院院刊, 2004, 19(3): 189-193.
[2] 王谦身, 等. 重力学[M]. 北京: 地震出版社, 2003.
[3] 顾功叙, Kuo J T, 刘克人, 等. 中国京津唐张地区时间上连续的重力变化与地震的孕育和发生[J]. 科学通报, 1997, 42(18): 1920-1930.
[4] 牛之俊, 马宗晋, 陈鑫连, 等. 中国地壳运动观测网络[J]. 大地测量与地球动力学, 2002, 22( 3): 88-93.
[5] 李强, 游新兆, 杨少敏, 等. 中国大陆构造变形高精度大密度GPS监测现今速度场[J]. 中国科学: 地球科学, 2012, 42(5): 629-632.
[6] 李辉, 申重阳, 孙少安, 等. 中国大陆近期重力场动态变化图像[J]. 大地测量与地球动力学, 2009, 29(3): 1-10.
[7] Chen Y T, G. H. D, Lu Z X. Variations of gravity before and after the Haichengearthquake, 1975, and the Tangshan earthquake, 1976[J]. Phys.Earth.Planet.Int, 18: 330-338.
[8] Kuo J T, Zheng J H, Song S H, et al. Determination of earthquake epicentroids by inversion of gravity variation data in the BTTZ region, China[J]. Tectonophysics, 1999, 312: 267-281.
[9] Zhu Y Q, Z. F. B., Zhou J C, et al. Gravity Measurements and Their Variations before the 2008 WenchuanEarthquake[J]. BSSA, 2010, 100(5B): 2815-2824.
[10] 陈石, 王谦身, 祝意青, 等. 汶川M_S8.0震前区域重力场变化与震质中研究[J]. 地球物理学进展, 2011, 26(4): 1147-1156.
[11] 陈石, 王青华, 王谦身, 等. 云南鲁甸M_S6.5地震震源区和周边三维密度结构及重力场变化[J]. 地球物理学报, 2014, 57(9): 3080-3090.
[12] Chen S, Liu M, Xing L, et al. Gravity increase before the 2015 M_W7.8 Nepal earthquake [J]. Geophysical Research Letters, 2016, 43, doi:10.1002/2015GL066595.
[13] Chen S, Jiang C S, Zhuang J C. Statistical evaluation of efficiency and possibility of earthquake predictions with gravity field variations and its analytic signal in western china [J]. Pure and Applied Geophysics, 2016, 173(1): 305-319.
[14] 陈石, 徐伟民, 蒋长胜. 中国大陆西部重力场变化与强震危险性关系[J]. 地震学报, 2015, 37(4): 562-574.
[15] Cordell L. A scattered equivalentsource method for interpolationand gridding of potentialfield data in three dimension[J]. Geophysics, 1992, 57: 629-636.
[16] 刘泽容, 杜庆龙, 蔡忠. 应用变差函数定量研究储层非均匀性[J]. 地质评论, 1993, 39(4): 297-301.
[17] Hansen R. Interpretive gridding by anisotropic kriging[J]. Geophysics, 1993, 58(10): 1491-1497.
[18] Deutsch C V, Journel A G. GSLIB: Geostatistical Software Library and User’s Guide[M]. NewYork: OxfordUniversity Press, 1998, 1-369.
[19] Journel A G, Huijbregts C. Mining geostatistics[M]. London: Academic Press, 1978: 148-218.
[20] 张培震, 邓起东, 张竹琪, 等. 中国大陆的活动断裂、地震灾害及其动力过程[J]. 中国科学: 地球科学, 2013, 43: 1607-1620.
[21] 殷秀华, 史志宏, 刘占坡, 等. 中国大陆区域重力场的基本特征[J]. 地震地质, 1980, 2(4): 69-75.
[22] 殷秀华, 史志宏, 刘占坡, 等. 华北北部均衡重力异常的初步研究[J]. 地震地质, 1982, 4(4): 27-34.
[23] 李辉, 徐如刚, 申重阳, 等.大华北地震动态重力监测网分形特征研究[J].大地测量与地球动力学, 2010, 30(5): 15-18.
[24] 祝意青, 闻学泽, 张晶, 等.华北中部重力场的动态变化及其强震危险含义[J].地球物理学报, 2013, 56(2): 531-541.
[25] Davis J. Statistics and Data Analysis in Geology[M]. Wiley, 1973.
[26] Randal B. Variogram Tutorial[M], Golden Software, Inc.
[27] 靳松, 朱筱敏, 钟大康.变差函数在沉积微相自动识别中的应用[J].石油学报, 2006, 27(3): 57-60.
[28] 陈永清, 黄静宁, Zhai Xiaoming, 等.多尺度地球化学勘查聚焦找矿靶区以滇东Pt地球化学省勘查研究为例[J]. 中国科学D辑: 地球科学, 2009, 52(5): 627-637.
[29] 李晓军, 王长虹, 朱合华. Kriging插值方法在地层模型生成中的应用[J]. 岩土力学, 2009, 30(1): 157-162.
[30] 王云专, 付检刚, 曲忠勇, 等. 变差函数在储层地质建模中的应用[J]. 科学技术与工程, 2010, 10(29): 7147-7150.
[31] 孟小红, 侯建全, 梁宏英, 等. 离散光滑插值方法在地球物理位场中的快速实现[J]. 物探与化探, 2002, 26(4): 302-306.
[32] 郭良辉, 孟小红, 郭志宏, 等. 地球物理不规则分布数据的空间网格化法[J]. 物探与化探, 2005, 29(5): 438-442.
[33] 孙树海, 刘雪莹.用地震资料预测孔隙度的广义克里格法[J]. 地球物理学报, 1993, 36(6): 798-804.
[34] 孟健, 马小明. Kriging空间分析法及其在城市大气污染中的应用[J]. 数学实践与认识, 2002, 32(2): 309-312.
[35] 常文渊, 戴新刚, 陈洪武.地质统计学在气象要素场插值的实例研究[J].地球物理学报, 2004, 47(6): 982-990.
[36] 程勖, 杨毅恒. 变差函数中加权系数拟合方法的改善[J]. 地球物理学进展, 2009, 24(1) : 362-366.
[37] Sandwell D T, Smith W H F. Marine gravity anomaly from Geosat and ERS 1 satellite altimetry[J]. Journal of Geophysical Research: Solid Earth, 1997, 102: 10039-10054.
[38] Sandwell D T, Smith W H F. Global marine gravity from retrackedGeosat and ERS-1 altimetry: Ridge Segmentation versus spreading rate[J]. J Geophys Res, 2009, 114: B01411, doi:10.1029/2008JB006008.
[39] Sandwell D, Garcia E, Soofi K, et al. Toward 1-mGal accuracy in global marine gravity from CryoSat-2, Envisat, and Jason-1[J]. The Leading Edge, 2013, 32: 892-899.
[40] 章传银, 郭春喜, 陈俊勇, 等. EGM2008地球重力场模型在中国大陆适用性分析[J]. 测绘学报, 2009, 38(4): 283-289.
[41] 杨金玉, 张训华, 张菲菲, 等. EGM2008地球重力模型数据在中国大陆地区的精度分析[J]. 地球物理学进展, 2012, 27(4): 1298-1306.
[42] 付广裕, 祝意青, 高尚华, 等. 川西地区实测自由空气重力异常与EGM2008模型结果的差异[J]. 地球物理学报, 2013, 56(11): 3761-3769.
[43] 张国民, 马宏生, 王辉, 等.中国大陆活动地块边界带与强震活动[J]. 地球物理学报, 2005, 48(3): 602-610.