以鄂尔多斯南缘地区布格重力异常数据及2014—2017年相对重力观测数据为基础，采用欧拉反褶积方法，对引起重力变化的场源深度进行反演，并对空间分布规律予以解释。通过构建理论模型，反演得到最优模型参数，并对实际数据进行计算和分析。为了减弱和消除欧拉解的发散性，利用水平梯度滤波法优化反演结果。结果表明：①构造指数为1时，适合对相对重力数据进行反演，当滑动窗口选择4—8倍测点间距时，可以获得较为可靠的场源参数；②在主要断裂附近，不同时间尺度的场源呈较好的一致性特征；③引起重力变化的场源深度集中在10—30 km，与鄂尔多斯南缘地区的震源深度基本一致；④不同时间尺度的场源位置相对分散，表明引起重力变化的物质流动具有随机性和波动性。
In this paper, based on the complete Bouguer gravity anomaly data and the relative gravity data from 2014 to 2017 in the southern margin of Ordos block, we adopt the Euler deconvolution method to invert for the field source depths that caused the gravity variation and interpret the characteristics of their spatial distribution. By constructing the theoretical model, the optimal model parameters are obtained, and the actual data are calculated and analyzed based on the optimal model parameters. In order to reduce and eliminate the divergence of Euler's solution, we use the horizontal gradient filtering method to optimize the inversion results. The results show that:when the structural index is 1, it is suitable for the inversion of relative gravity data, and when the sliding window is 4-8 times the distance between observation points, more reliable field source parameters can be obtained. The horizontal gradient filtering method can effectively weaken and eliminate the divergence of inversion results. The inversion results show that the field source characteristics of different time scales are in good agreement near the main faults, and the depth of the field source that causes gravity variation mainly concentrates within 10-30 km, which is basically the same as the focal depth in the southern margin of Ordos block. In addition, the location of field sources at different time scales is relatively scattered, indicating that the material flow that causes gravity changes is random and volatile.
2021,42(3): 84-96 收稿日期：2020-11-17
DOI：10.3969/j.issn.1003-3246.2021.03.010
基金项目：陕西省地震局启航与创新基金课题（项目编号：QC202013，QC202023）；陕西省自然科学基础研究计划青年项目（项目编号：2021JQ-962，2021JQ-980）
作者简介：张永奇(1985-),男,工程师,从事基于大地测量资料的地壳形变正反演研究工作。E-mail:zyq19851115@163.com
参考文献：
曹建平, 祝意青, 韩美涛, 等. 陕西重力网优化改造及监测能力评估[J]. 地震地磁观测与研究, 2014, 35(Z1):165-170.
陈石, 王谦身, 祝意青, 等. 青藏高原东缘重力导纳模型均衡异常时空特征[J]. 地球物理学报, 2011a, 54(1):22-34.
陈石, 王谦身, 徐伟民, 等. 华北地区热均衡、重力均衡与深部构造[J]. 地球物理学报, 2011b, 54(11):2 864-2 875.
陈石, 王谦身, 祝意青, 等. 汶川MS 8.0震前区域重力场变化与震质中研究[J]. 地球物理学进展, 2011c, 26(4):1 147-1 156.
陈石, 王青华, 王谦身, 等. 云南鲁甸MS 6.5地震震源区和周边三维密度结构及重力场变化[J]. 地球物理学报, 2014, 57(9):3 080-3 090.
陈石, 徐伟民, 蒋长胜. 中国大陆西部重力场变化与强震危险性关系[J]. 地震学报, 2015, 37(4):562-574.
陈运泰, 顾浩鼎, 卢造勋. 1975年海城地震和1976年唐山地震前后的重力变化[J]. 地震学报, 1980, 2(1):21-31.
杜建军. 渭河盆地东南缘主要断裂晚更新世以来的活动性及灾害效应[D]. 北京:中国地质科学院, 2016.
范美宁. 欧拉反褶积方法的研究与应用[D]. 吉林大学, 2006.
高倩, 陈石. 基于欧拉反褶积方法计算川滇交界重力变化场源特征[J]. 地球物理学进展, 2015, 30(2):503-509.
管志宁. 地磁场与磁力勘探[M]. 北京:地质出版社, 2005.
刘芳, 祝意青, 梁伟锋, 等. 基于欧拉反褶积方法估算华北地区重力变化的等效源参数[J]. 地震, 2016, 36(4):163-170.
鲁宝亮, 范美宁, 张原庆. 欧拉反褶积中构造指数的计算与优化选取[J]. 地球物理学进展, 2009, 24(3):1 027-1 031.
马宗晋, 高祥林, 宋正范. 中国布格重力异常水平梯度图的判读和构造解释[J]. 地球物理学报, 2006, 49(1):106-114.
申重阳, 李辉, 孙少安, 等. 2008年于田MS 7.3地震前重力场动态变化特征分析[J]. 大地测量与地球动力学, 2010,
30(4):1-7.
申重阳, 谈洪波, 郝洪涛, 等. 2009年姚安MS 6.0地震重力场前兆变化机理[J]. 大地测量与地球动力学, 2011, 31(2):17-23.
石磊, 郭良辉, 孟小红. 估计磁场源位置的倾斜角梯度算法改进[J]. 物探与化探, 2012, 36(2):287-292.
王卫东. 西安地裂缝形成的区域稳定性动力学背景研究[D]. 长安大学, 2009.
姚长利, 管志宁, 吴其斌, 等. 欧拉反演方法分析及实用技术改进[J]. 物探与化探, 2004, 28(2):150-155.
翟丽娜, 贾晓东, 焦明若, 等. 基于欧拉反褶积方法计算辽宁地区重力变化场源特征[J]. 地震研究, 2020, 43(2):388-393.
张永奇, 曹建平, 郝洪涛, 等. 陕西地区流动重力历史资料清理分析[J]. 地震地磁观测与研究, 2014, 35(Z1):171-176.
郑秋月, 王青华, 刘东, 等. 基于时变重力数据的川滇地区场源模型反演与解释[J]. 地震, 2021, 41(1):205-218.
祝意青, 徐云马, 吕弋培, 等. 龙门山断裂带重力变化与汶川8.0级地震关系研究[J]. 地球物理学报, 2009, 52(10):2 538-2 546.
祝意青, 梁伟锋, 赵云峰, 等. 2017年四川九寨沟MS 7.0地震前区域重力场变化[J]. 地球物理学报, 2017, 60(10):
4 124-4 136.
祝意青, 张勇, 张国庆, 等. 21世纪以来青藏高原大地震前重力变化[J]. 科学通报, 2020, 65(10):622-632.
Chen Shi, Liu M, Xing L, et al. Gravity increase before the 2015 MW 7.8 Nepal earthquake[J]. Geophysical Research Letter, 2016a, 43. doi:10.1002/2015GL066595.
Chen Shi, Jiang C S, Zhuang J C. Statistical evaluation of efficiency and possibility of earthquake predictions with gravity field variation and its analytic signal in western China[J]. Pure and Applied Geophysics, 2016b, 173(1):305-319.
Pavlis N K, Holmes S A, Kenyon S C, et al. The Development and Evolution of the Earth Gravitational Model 2008 (EGM2008)[J]. Journal of Geophysical Research, 2012, 117:B04406.
Reid A B, Allsop J M, Granser H, et al. Magnetic interpretation in three dimensions using Euler deconvolution[J]. Geophysics, 1990, 55(1):80-91.
Sandwell D T, Smith W H F. Global Marine Gravity from Retracked Geosat and ERS-1 Altimetry:Ridge Segmentation Versus Spreading Rate[J]. Journal of Geophysical Research:Solid Earth, 2009, 114(B1).
Thompson D T. EULDPH:A new technique for making computer-assisted depth estimates from magnetic data[J]. Geophysics, 1982, 47(1):31-37.
