最新全球重力数据库V23中陆域重力资料质量评估
张春灌1, 袁炳强1, 张国利2
1.西安石油大学 地球科学与工程学院, 陕西 西安 710065
2.中国地质调查局 天津地质调查中心, 天津 300170

作者简介:张春灌(1981-),男,江西于都人,副教授,主要从事综合地球物理勘探及构造地球物理研究.E-mail:chunguan-zhang@163.com

摘要

为了评估最新全球重力数据库V23中陆域重力资料的质量问题,选择了美国、秘鲁、爱尔兰、南非、肯尼亚、澳大利亚、中国等7个国家的8个区块约133万km2的重力数据进行对比分析。基于V23中自由空间重力异常数据,利用纯球坐标系内重力校正方法获得了这些区块的计算布格重力异常。利用相关分析的方法对这些区块的实测布格重力异常和计算布格重力异常进行相关分析,得到了相应的相关系数。通过对实测布格重力异常和计算布格重力异常之间相关系数及差值特征的综合分析,评估了全球重力数据库V23中选取的8个区块重力数据的质量。研究结果表明,斯克里普斯海洋研究所发布的最新全球重力数据库V23中陆域重力资料融合了大量测量年代较早的地面重力数据或者航空重力数据,而近十几年新完成的地面重力资料或者航空重力资料可能未融合进该数据库。具有大量年代较早的地面或者航空重力测量资料的区域,该数据库中陆域重力资料的质量较高,而重力工作程度较低的区域,其资料的质量则较低。

关键词: 质量评估; 陆域重力资料; 全球重力数据库; 斯克里普斯海洋研究所
中图分类号:P312.9 文献标志码:A 文章编号:1001-8166(2017)01-0075-08
Quality Evaluation of Land Gravity Data in the Latest Global Gravity Database V23
Zhang Chunguan1, Yuan Bingqiang1, Zhang Guoli2
1.School of Earth Science and Engineering, Xi’an Shiyou University, Xi’an 710065, China;
2.Tianjin Geological Survey Center, China Geological Survey, Tianjin 300170, China

First author:Zhang Chunguan (1981-), male, Yudu County, Jiangxi Province, Associate professor. Research areas include integrated geophysical exploration and tectonophysics. E-mail:chunguan-zhang@163.com

Abstract

In order to evaluate the quality of the land gravity data in the latest global gravity database V23, the authors chose the gravity data of eight blocks with a total area of 1 330 000 km2 to compare and analyze, and these blocks involved seven countries, including the United States, Peru, Ireland, South Africa, Kenya, Australia, and China. Based on the free-air gravity anomaly data of the latest global gravity database V23, the authors obtained the calculated Bouguer gravity anomaly using the gravity correction methods within the pure spherical coordinate system in these eight blocks. Then, the correlation coefficients between the measured Bouguer gravity anomaly and the calculated Bouguer gravity anomaly were calculated by the correlation analysis method in these eight blocks. Finally, through comprehensive analysis of the features of these correlation coefficients and differences between the measured Bouguer gravity anomaly and the calculated Bouguer gravity anomaly in these eight blocks, the quality of the gravity data of these eight blocks was evaluated in the latest global gravity database V23. The results showed that the latest global gravity database V23, released by Scripps Institution of Oceanography, integrated a large number of the ground or airborne gravity data measuring in an earlier era, and the newly surveyed ground or airborne gravity data may not be integrated into the database. The quality of land gravity data is relatively high in the areas with a large number of older ground or airborne gravity data, otherwise it is low in the zones with a lower gravity working degrees.

Keyword: Quality evaluation; Land gravity data; Global gravity database; Scripps Institution of Oceanography.
1 引言

2014年10月, 加利福尼亚大学圣地亚哥分校斯克里普斯海洋研究所(Scripps Institution of Oceanography, University of California San Diego)发布了最新全球海域自由空间重力数据库V23, 其数据覆盖了80.738° S~80.738° N的区域。虽然该重力数据库被称为全球海洋重力(Global marine gravity)数据库, 但是上述范围内的陆域亦有相应的重力数据。海域重力数据基于卫星测高数据计算而来[1~3], 数据精度较高。现代卫星测高重力异常精度达到4× 10-5m/s2左右、空间分辨率接近8 km, 可以满足1∶ 100 万比例尺面积性海洋重力工作的需要[4~8]。陆域重力数据没有明确的来源, 可能融合了卫星重力数据、航空重力数据、地面重力数据等多种数据。不同地区重力工作程度差异较大, 使得陆域重力数据精度在不同地区差异较大, 完全基于卫星重力的区域其异常精度较低[9~11], 而主要基于航空重力数据及地面重力数据的区域其异常精度显然较高。重力数据在资源勘探、国家安全、军事等领域具有重要的应用价值[12~14], 在人们使用最新全球重力数据库V23中陆域重力资料时, 显然涉及资料质量的高低问题。

本文拟基于全球重力数据库V23中自由空间重力异常数据以及高精度地形高程数据, 计算陆域不同地区的布格重力异常, 并对计算的布格重力异常与实测布格重力异常进行对比分析, 进而对这些地区重力资料质量进行评估。

2 方法与技术
2.1 纯球坐标系内重力校正

利用全球重力数据库V23中自由空间重力异常资料计算布格重力异常, 本文采用纯球坐标系内重力校正方法进行地形校正以及中间层校正[15]

球坐标系内重力地形校正和中间层校正, 就是要消除水准面到观测面之间的质量在测点处产生的重力值。按照中华人民共和国国土资源部发布的《区域重力调查规范》(DZ/T 0082-2006)规定, 上述2项校正的计算范围为0~166.7 km。

定义地壳中“ 球壳型扇形地块” 为:具有球壳的曲率, 上下表面之间的厚度为(R2-R1), 东西侧面与z轴之间夹角之差为(θ 21)。该扇形地块南北侧面的纬度分别为φ 2(φ 2=90° -θ 2)和φ 1(φ 1=90° -θ 1)。则该球壳型扇形地块在由地心出发、向上过北极的z轴上产生的Δ g(zi)表达式为:

Δ g(zi)= (λ2-λ1)360°× 2πGρzi3× rzi2+rzicosθ+3cos2θ-2× rzi2-2rzicosθ+1-3cosθ sin2θ

× ln rzi-cosθ+rzi2-2rzicosθ+1|r=R2, θ =θ 2r=R1, θ =θ 1=Δ g(zi, r=R2, θ 2, θ 1)-Δ g(zi, r=R1, θ 2, θ 1) (1)

为了使每一个重力测点P均成为“ 北极点” , 需要计算出每一个地形点Q相对于该“ 北极点” P的经度和纬度。为此要借助于球面三角公式:

α=arccos[sinφ·sinφ0+cosφ·cosφ0·cos(λ-λ0)β=arctancosφ·sin(λ-λ0)sinφ·cosφ0-cosφ·sinφ0·cos(λ-λ0)(2)

式中:(λ 0, φ 0)和(λ , φ )分别是重力测点P和地形点Q的地理经纬度, (λ -λ 0)是Q点和P点的经度差。

公式(1)~(2)即可实现重力地形校正值的计算。在极限情况下, 即zi=Ri=R2, R1=R0, θ 1=0° , θ 2, λ 21时, 按洛必达极限法则可推导出中间层校正计算公式。

2.2 相关分析

为了详细分析实测布格重力异常与计算布格重力异常之间的关系, 本文利用了相关分析方法[16]

现象之间的相互联系, 常表现为一定的关系, 其中一个或若干个起着影响作用的变量用X表示, 与之相关的另一现象变量用Y表示。通常用相关系数r来表征这2个量的相互关联程度。

r= (x-x̅)(y-y̅)n(x-x̅)n2(y-y̅)n2=(x-x̅)(y-y̅)(x-x̅)2(y-y̅)2(3)

式中:n是变量X的个数, x̅是变量X的平均值, y̅是变量Y的平均值。

相关系数r的取值范围:-1≤ r≤ 1。r> 0为正相关, r< 0为负相关; |r|=0表示不存在相关关系; |r|=1表示完全线性相关; 0< |r|< 1表示存在不同程度线性相关; |r|≤ 0.3为不存在线性相关; 0.3< |r|≤ 0.5为低度线性相关; 0.5< |r|≤ 0.8为显著线性相关; |r|> 0.8为高度线性相关。

3 对比与分析

鉴于海域卫星测高重力异常精度较高, 已经用于中等比例尺的海洋地质与资源调查工作[4], 因此本文对于全球重力数据库V23中海域数据不进行对比与分析, 主要是对于陆域的重力数据进行对比与分析。作者选择了爱尔兰、美国、中国等北半球3个国家的4个区块以及澳大利亚、南非、秘鲁、肯尼亚等南半球4个国家的4个区块进行了数据的对比与分析, 区块分布见图1。所选的8个区块分布于六大洲的不同纬度、不同海拔、不同重力工作程度的地区, 总面积约133万km2。地形高程数据利用斯克里普斯海洋研究所于2014年11月发布的最新全球地形高程数据库SRTM30_PLUS V10, 数据覆盖了80.738° S~80.738° N的区域[17, 18]。作者首先基于全球重力数据库V23中提取的自由空间重力异常数据和地形高程数据库SRTM30_PLUS V10, 利用纯球坐标系内重力校正方法对所选取的8个区块进行地形校正以及中间层校正, 获得这些区块的计算布格重力异常, 然后利用相关分析的方法对所选取8个区块的实测布格重力异常和计算布格重力异常分别进行相关分析, 获得相应的相关系数, 结果见图2和图3。

图1 数据对比区块分布图Fig.1 The block distribution map for data comparison

图2为澳大利亚、南非、秘鲁、肯尼亚等南半球4个国家的4个区块的数据对比图。其中, 澳大利亚区块的实测布格重力异常(图2(a1))数据来源于澳大利亚矿产资源局于1979年出版的1∶ 1 000万全国布格重力异常图[19]; 南非区块的实测布格重力异常(图2(b1))数据来源于Webb等[20]于2004年发表的布格重力异常图, 该图数据整合了1939— 1986年的地面重力资料; 秘鲁区块的实测布格重力异常(图2(c1))数据来源于早期的地面重力路线测量资料, 重力测点分布见图2(c4); 肯尼亚区块的实测布格重力异常(图2(d1))数据来源于2008年实测的1∶ 25万地面重力资料[21], 重力测点分布见图2(d4)中的黑色十字点, 而早期该地区有地面重力路线测量资料(图2(d4))。

澳大利亚及南非区块的实测布格重力异常与计算布格重力异常之间存在明显的相关关系, 这种相关关系在相关系数图(图2(a4)和图2 (b4))上有非常明显的显示。相关系数大于0.8的区域占绝大多数, 仅有很小面积的相关系数为0.5~0.8, 表明计算布格重力异常与实测布格重力异常的变化特征一致性非常好。澳大利亚区块实测布格重力异常与计算布格重力异常的差值绝对值小于10× 10-5m/s2的区域占总面积的97.99%, 而南非区块则为87.79%, 表明这2种数据的差异总体较小(图2(a3)和图2 (b3))。实测布格重力异常与计算布格重力异常之间的相关系数及差值特征表明, 全球重力数据库V23中澳大利亚及南非区块的重力数据质量较高, 可能融合了地面重力资料或者航空重力资料。

秘鲁区块实测布格重力异常与计算布格重力异常之间宏观上呈一定的相关关系(图2c), 而肯尼亚区块实测布格重力异常与计算布格重力异常之间宏观上并无明显的相关关系(图2d)。秘鲁区块相关系数大于0.8的区域约占一半面积, 和重力测点分布基本一致; 相关系数介于0.5~0.8的面积较小, 而介于-0.9~0.5的面积较大, 基本处于无重力测点分布的区域(图2(c4))。肯尼亚区块相关系数大于0.8的区域仅占很小的面积, 和早期的重力路线测点分布基本一致; 相关系数介于0.5~0.8的面积也较小, 而介于-1~0.5的面积则较大(图2(d4))。秘鲁区块实测布格重力异常与计算布格重力异常的差值绝对值小于10× 10-5m/s2的区域占总面积的58.01%, 而肯尼亚区块则为82.17%, 表明这2种数据的差异在肯尼亚区块总体较小, 秘鲁区块则较大(图2(c3)和图2 (d3))。实测布格重力异常与计算布格重力异常之间的相关系数及差值特征表明, 全球重力数据库V23中秘鲁及肯尼亚区块可能融合了早期的重力路线测量资料, 这些重力测点较少, 以致重力数据总体质量较低, 而在重力测点较密集的局部区域的数据质量则较高。

图2 南半球区块数据对比图
(a)澳大利亚区块; (b)南非区块; (c)秘鲁区块; (d)肯尼亚区块; 1.实测布格重力异常; 2.计算布格重力异常; 3.实测布格重力异常与计算布格重力异常的差值; 4.实测布格重力异常与计算布格重力异常的相关系数; c4中黑色十字点为重力路线测点; d4中黑色十字点为2008年重力测点、红色圆点为早期重力路线测点
Fig.2 The data comparison map in the Southern Hemisphere
(a)Australia block; (b)South Africa block; (c)Peru block; (d)Kenya block; 1.Measured Bouguer gravity anomaly; 2.Calculated Bouguer gravity anomaly; 3.The difference between the measured Bouguer gravity anomaly and the calculated Bouguer gravity anomaly; 4.The correlation coefficient between the measured Bouguer gravity anomaly and the calculated Bouguer gravity anomaly; The black cross points are the points of gravity lines in c4; The black cross points are the gravity points in 2008 and the red dots are the points of early gravity lines in d4

图3为爱尔兰、美国、中国等北半球3个国家的4个区块的数据对比图。其中, 爱尔兰区块的实测布格重力异常(图3(a1))数据来源于Readman等[22]于1997年发表的布格重力异常图, 该图数据整合了1987年以前的地面重力资料, 原始重力测点密度为0.42点/km2, 相当于点距为1.55 km; 美国区块的实测布格重力异常(图3(b1))数据来源于Seeley等[23]于2003年发表的布格重力异常图, 该图数据整合了1997年以前的地面重力资料; 渭河区块的实测布格重力异常(图3(c1))数据来源于1991年编制的1∶ 50万布格重力异常图; 青藏区块的实测布格重力异常(图3(d1))数据来源于2006年完成的1∶ 100万地面重力资料[24]

图3 北半球区块数据对比图
(a)爱尔兰区块; (b)美国区块; (c)渭河区块; (d)青藏区块; 1.实测布格重力异常; 2.计算布格重力异常; 3.实测布格重力异常与计算布格重力异常的差值; 4.实测布格重力异常与计算布格重力异常的相关系数
Fig.3 The data comparison map in the Northern Hemisphere
(a)Ireland block; (b)America block; (c)Weihe block; (d)Qinghai Tibet block; 1.Measured Bouguer gravity anomaly; 2.Calculated Bouguer gravity anomaly; 3.The difference between the measured Bouguer gravity anomaly and the calculated Bouguer gravity anomaly; 4.The correlation coefficient between the measured Bouguer gravity anomaly and the calculated Bouguer gravity anomaly

爱尔兰区块实测布格重力异常与计算布格重力异常之间存在明显的相关关系(图3(a4)), 相关系数大于0.8的区域占绝大多数, 仅有很小面积的相关系数介于0.5~0.8以及中北部极小面积的相关系数介于0.3~0.5。该区块实测布格重力异常与计算布格重力异常的差值绝对值小于10× 10-5m/s2的区域占总面积的99.93%, 表明这2种数据的差异总体很小(图3(a3))。实测布格重力异常与计算布格重力异常之间的相关系数及差值特征表明, 全球重力数据库V23中该区块可能融合了地面重力资料, 重力数据总体质量较高。

美国及渭河区块实测布格重力异常与计算布格重力异常之间宏观上呈明显的相关关系, 相关系数大于0.8的区域占大多数(图3(b4)和图3 (c4))。美国区块中北部有大面积相关系数介于0.5~0.8以及北部边缘处有小面积相关系数介于0.2~0.5, 表明计算布格重力异常与实测布格重力异常在该区中北部一致性较差; 渭河区块北部和南部有大面积相关系数介于0.5~0.8以及小面积相关系数介于0.1~0.5, 表明计算布格重力异常与实测布格重力异常在该区北部和南部一致性较差。美国区块实测布格重力异常与计算布格重力异常的差值绝对值小于10× 10-5m/s2的区域占总面积的29.99%, 而渭河区块则为85.31%, 表明这2种数据的差异在渭河区块总体较小, 美国区块则较大(图3(b3)和图3(c3))。实测布格重力异常与计算布格重力异常之间的相关系数及差值特征表明, 全球重力数据库V23中美国及渭河区块的重力数据可能融合了地面重力资料, 但是美国区块的重力异常数值可能做过适当的处理, 所以渭河区块重力数据总体质量较高而美国区块重力数据总体质量较低。

青藏区块实测布格重力异常与计算布格重力异常之间无明显的相关关系(图3(d4)), 相关系数皆小于0.8, 介于0.5~0.8的面积亦非常小, 表明实测布格重力异常与计算布格重力异常的变化特征一致性很差。该区块实测布格重力异常与计算布格重力异常的差值绝对值小于10× 10-5m/s2的区域占总面积的43.52%, 表明这2种数据的差异较大(图3(d3))。实测布格重力异常与计算布格重力异常之间的相关系数及差值特征表明, 全球重力数据库V23中青藏区块的重力数据总体质量很低。

4 结 语

上述8个区块的数据对比与分析表明, 斯克里普斯海洋研究所发布的最新全球重力数据库V23中陆域重力资料融合了大量测量年代较早的地面重力数据或者航空重力数据, 而近十几年新完成的地面重力资料或者航空重力资料可能未融合进该数据库。也就是说, 对于该数据库中陆域重力资料, 具有大量年代较早的地面或者航空重力测量资料的区域, 其数据质量较高, 而那些重力工作程度较低的区域, 其数据质量则较低。

鉴于不同地区重力资料质量不同, 在利用全球重力数据库V23中陆域重力资料时, 可以先查询所涉及工作区的重力资料分布情况或者重力工作程度情况(比如GETECH地球物理公司可以查询), 然后再分析该区重力资料质量的高低。

The authors have declared that no competing interests exist.

参考文献
[1] Sand well D T, Muller R D, Smith W H F, et al. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure[J]. Science, 2014, 346(6 205): 65-67. [本文引用:1]
[2] Sand well 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, 2009, 114(1): 1-18. [本文引用:1]
[3] Sand well D T, Garcia E, Soofi K, et al. Towards 1-mGal Global Marine Gravity from CryoSat-2, Envisat, and Jason-1[J]. The Leading Edge, 2013, 32(8): 892-899. [本文引用:1]
[4] Zhang Minghua, Zhang Jiaqiang. Resolution of modern satellite altimetric gravity anomaly and its application to marine geological survey[J]. Geophysical & Geochemical Exploration, 2005, 29(4): 295-298.
[张明华, 张家强. 现代卫星测高重力异常分辨能力分析及在海洋资源调查中应用[J]. 物探与化探, 2005, 29(4): 295-298. ] [本文引用:2]
[5] Liang Ziliang, Chen Lu, Xie Kun, et al. Inversion of marine gravity anomalies and geoid using multi-satellite altimeter data[J]. Journal of Geodesy and Geodynamics, 2015, 35(1): 40-44.
[梁子亮, 陈路, 解琨, . 利用多代卫星测高数据反演海洋重力异常及大地水准面[J]. 大地测量与地球动力学, 2015, 35(1): 40-44. ] [本文引用:1]
[6] Liu Shanwei, Li Jiajun, Wan Jianhua, et al. Calculation of gravity anomalies over China Sea and its vicinity based on multi-generation satellite altimetry data[J]. Marine Sciences, 2015, 39(12): 130-134.
[刘善伟, 李家军, 万剑华, . 利用多代卫星测高数据计算中国近海及邻域重力异常[J]. 海洋科学, 2015, 39(12): 130-134. ] [本文引用:1]
[7] Guan Yihe, Sheng Hui, Liu Shanwei, et al. Inversion of the gravity anomalies by using multi-generation satellite altimeter data in the South China Sea[J]. Hydrographic Surveying and Charting, 2016, 36(1): 11-14.
[管一鹤, 盛辉, 刘善伟, . 联合多代卫星测高资料反演中国南海重力异常[J]. 海洋测绘, 2016, 36(1): 11-14. ] [本文引用:1]
[8] Hu yi, Wang Liming, Zhong Guicai, et al. Gravity and magnetic characteristics of the Weddell Sea and its tectonic significance[J]. Advances in Earth Science, 2015, 30(11): 1 231-1 238.
[胡毅, 王立明, 钟贵才, . 威德尔海的重磁场特征及其构造意义[J]. 地球科学进展, 2015, 30(11): 1 231-1 238. ] [本文引用:1]
[9] Xu Haijun, Zhang Yongzhi, Duan Hurong, et al. Gravity anomaly detected by GOCE satellite in China[J]. Progress in Geophysics, 2012, 27(2): 404-408.
[徐海军, 张永志, 段虎荣, . GOCE卫星监测的中国区域重力异常[J]. 地球物理学进展, 2012, 27(2): 404-408. ] [本文引用:1]
[10] Xie Rong, Liu Yawen, Li Xiangxiang. Key technologies of Earth observation satellite data integration system under big data environment[J]. Advances in Earth Science, 2015, 30(8): 855-862.
[谢榕, 刘亚文, 李翔翔. 大数据环境下卫星对地观测数据集成系统的关键技术[J]. 地球科学进展, 2015, 30(8): 855-862. ] [本文引用:1]
[11] Miao Chunsheng, Cheng Yuan, Wang Jianhong, et al. Data fusion of offshore SST from China FY and HY2 satellites and its application[J]. Advances in Earth Science, 2015, 30(10): 1 127-1 143.
[苗春生, 程远, 王坚红, . 中国风云卫星与海洋卫星近海SST资料融合技术及应用研究[J]. 地球科学进展, 2015, 30(10): 1 127-1 143. ] [本文引用:1]
[12] Liu Daizhi. A number of research directions in geophysics for national security[J]. Progress in Geophysics, 2007, 22(4): 1 327-1 331.
[刘代志. 国家安全地球物理学的若干研究方向[J]. 地球物理学进展, 2007, 22(4): 1 327-1 331. ] [本文引用:1]
[13] Xu Zunyi, Yan Lei, Ning Shunian, et al. Situation and development of marine gravity aided navigation system[J]. Progress in Geophysics, 2007, 22(1): 104-111.
[徐遵义, 晏磊, 宁书年, . 海洋重力辅助导航的研究现状与发展[J]. 地球物理学进展, 2007, 22(1): 104-111. ] [本文引用:1]
[14] Wang Liupeng, Guo Yanping, Feng Wei. GOCE satellite in the military applications[J]. Geospatial Information, 2011, 9(1): 11-15.
[王留朋, 郭燕平, 冯炜. GOCE重力卫星在军事上的应用前景分析[J]. 地理空间信息, 2011, 9(1): 11-15. ] [本文引用:1]
[15] An Yulin, Zhang Minghua, Huang Jinming, et al. The computation scheme and computation process for gravity correction values within the pure spherical coordinate system[J]. Geophysical & Geochemical Exploration, 2010, 34(6): 1-9.
[安玉林, 张明华, 黄金明, . 纯球坐标系内各项重力校正值计算方案和过程[J]. 物探与化探, 2010, 34(6): 1-9. ] [本文引用:1]
[16] Zhang C, Dong Y, Yuan B, et al. A genesis analysis of the regional gravity and magnetic anomalies in the northern part of eastern Xinjiang, Northwest China[J]. Petroleum Science and Technology, 2014, 32(17): 2 075-2 085. [本文引用:1]
[17] Smith W H F, Sand well D T. Global seafloor topography from satellite altimetry and ship depth soundings[J]. Science, 1997, 277(5 334): 1 957-1 962. [本文引用:1]
[18] Becker J J, Sand well D T, Smith W H F, et al. Global bathymetry and elevation data at 30 Arc seconds resolution: SRTM30_PLUS[J]. Marine Geodesy, 2009, 32(4): 355-371. [本文引用:1]
[19] Wellman P, Murray A S. Bouguer Gravity Anomalies[R]. Australia: Bureau of Mineral Resources, Geology and Geophysics, Department of National Development, 1979. [本文引用:1]
[20] Webb S J, Cawthorn R G, Nguuri T, et al. Gravity modeling of Bushveld Complex connectivity supported by Southern African Seismic Experiment results[J]. South African Journal of Geology, 2004, 107: 207-218, doi: DOI:10.2113/107.1-2.207. [本文引用:1]
[21] Yuan B, Xie W, Liu G, et al. Gravity field and tectonic features of Block L2 in the Lamu Basin[J]. Geophysical Prospecting, 2012, 60(1): 161-178. [本文引用:1]
[22] Readman P W, O’Reilly B M, Murphy T. Gravity gradients and upper-crustal tectonic fabrics, Ireland [J]. Journal of the Geological Society, 1997, 154: 817-828. [本文引用:1]
[23] Seeley J M, Rand y Keller G. Delineation of subsurface Proterozoic Unkar and Chuar Group sedimentary basins in northernArizona using gravity and magnetics: Implications for hydrocarbon source potential[J]. AAPG Bulletin, 2003, 87: 1 299-1 321, doi: DOI:10.1306/0319.03200198. [本文引用:1]
[24] Zhang Minghua, He Hao, Qiao Jihua, et al. Integration of Regional Gravity Survey in Qinghai-Tibet Plateau[M]. Beijing: Geological Publishing House, 2015.
[张明华, 贺颢, 乔计花, . 青藏高原区域重力调查成果综合研究[M]. 北京: 地质出版社, 2015. ] [本文引用:1]