Gravity isostaty of Ordos block and its surrounding regions
-
摘要: 为更加精细地了解鄂尔多斯及其邻区的地壳均衡状态，本文在考虑地壳密度横向变化的基础上计算了该区域的重力均衡异常。首先，利用研究区内重力异常和真实莫霍面深度进行反演得到壳幔密度差分布，以艾里均衡理论为基础计算了莫霍面的理论深度；其次，通过莫霍深度理论值与真实值对比来获得研究区的重力均衡异常；最后，结合区内地震分布，讨论了重力均衡异常与地震活动的关系。结果表明：研究区内盆地以负均衡异常为主，山地则表现为正均衡异常，鄂尔多斯地块内部基本均衡；研究区内中强地震主要分布在重力均衡异常梯度带上。Abstract: In order to better understand the isostasy of the crust beneath Ordos block and its surrounding regions, we take into account laterally variation of crustal density in the calculation of isostaty anomalies. Firstly, obtained by the inversion for regional gravity anomaly and Moho depth when calculating the isostasy anomaly; Secondly, we compared theoretical Moho depth and real one which is inverted from the receiver function so as to get the crustal isostasy in the studied regions. Finally, we discussed the relationship between the isostatic gravity anomaly and distribution of earthquakes in this region. The results indicate that the negative isostatic anomaly mainly corresponds to basins, and positive anomaly to mountainous regions, therefore the interior of Ordos block is basically in isostatic state. Moreover, the most major earthquakes mainly occurred on the isostatic anomaly gradient belt in the studied region.
-
Keywords:
- Ordos block /
- gravity isostasy anomaly /
- crust-mantle density difference /
- Moho depth
-
引 言
工程地震动参数是指在工程应用领域使用较多并具有较高参考价值的地震动参数，例如地震动的幅值、持时、反应谱、阿里亚斯（Arias）强度和卓越周期等。地震动参数是工程抗震设计的依据，通常普通工业与民用建筑需要重点参考地震动的峰值加速度（peak ground acceleration，简写为PGA），而重大基础设施则需要额外参考频谱和持时，边坡稳定以及滑坡体的安全性等需要重点参考阿里亚斯强度等参数。本文采用Irikura （1986）等改进了的经验格林函数方法，其可靠性已经过日本多次地震动数值模拟所验证（Miyake et al，2003 ；Irikura，Miyake，2011；Irikura et al，2017 ），且被广泛认可。由于小震记录受震源的破裂过程及传播介质复杂性的影响（李启程，2010），所以由小震记录合成的大地震记录自然也需考虑震源和传播介质的复杂性（李宗超等，2016）。
Kanamori （1979）尝试利用对大小地震地震矩的比值进行取整的方法来确定所需小震数目，并针对所用小震与大震记录并非取自同一地震和同一台站的情况提出了修正公式。Irikura等（1986）考虑了大、小地震之间的相似性关系，将断层大小、平均位错量和上升时间等震源参数联系起来，通过理论推导，提出了大小地震满足一定相似条件下的经验格林函数方法，将其称为改进经验格林函数方法。Irikura （1986）从理论上证明，并结合实际震例验证了经验格林函数的有效性。Dan等（1989）通过考虑大、小地震之间的应力降与平均位错的不同，采用布龙模型的小震远场傅里叶谱对大震子源进行了修正。
Irikura和Kamae （1994）为了扩展经验格林函数的应用范围，在大小地震长宽比相等但不等于位错比的假设下，提出合成地震动需考虑断层面上凹凸体的影响，第一次将凹凸体模型与经验格林函数方法结合起来。Somerville （1999）等通过回归分析得到了确定凹凸体各参数的公式，并基于断层中凹凸体特征的统计分析，提出了确定断层破裂面上凹凸体的标准，即凹凸体为滑动量大于断层上平均滑动量1.5倍或以上的区域，而确定凹凸体的方法是去掉断层中滑动量小于0.3倍断层平均滑动量的断层边界后，得到新的矩形断层，重复此过程直到得到所有断层的边界均大于或等于0.3倍的断层平均滑动量，再将断层上大于等于平均滑动1.5倍的区域视为凹凸体。鉴于凹凸体模型具有与经验格林函数法相结合的优点且模拟结果更加精确，本文拟采用Irikura （1986）改进的经验格林函数法对2016年日本熊本MJ7.3 地震的地震动进行数值模拟，进行地震动模拟时同时将采用包含凹凸体的震源模型，以期获取熊本地震影响区域内诸如PGA，阿里亚斯强度，卓越周期等的地震动参数的空间分布特征。
1. 基本地震信息及震源参数估计
2016年4月16日，日本熊本县发生MJ7.3地震，震中位置为（32.74°N，130.76°E），位于熊本县城区，震源深度为13.1 km，发震断层尺度为40 km×16 km。该区域地形相对比较平缓（图1），东部为九州山地，西部为日本海。而MJ5.4余震的震中位置为（32.75°N，130.76°E），距离主震位置较近，震源深度为15 km，应力降为0.9 MPa。主震和余震详细信息列于表1。本文利用MJ5.4余震记录作为格林函数，选取K-net地震台网中可同时记录到这两次地震的47个强震台站的地震记录来模拟地震动场（National Research Institute for Earth Science and Disaster Resilience，2016），台站分布如图1所示。
![]()
图 1 震源位置、近场区域地形及台站分布Figure 1. Source location and distribution of stations and terrains in the near field of 2016 Kumamoto earthquake表 1 熊本MJ7.3地震和MJ5.3余震相关参数（Disaster Prevention Research Institute，Kyoto University，2016）Table 1. The relevant parameters of Kumamoto MJ7.3 earthquake and its aftershock MJ5.3（Disaster Prevention Research Institute，Kyoto University，2016）震级 震中位置 发震断层 震源深度/km 地震矩/（1017 N·m） 北纬/° 东经/° 走向/° 倾角/° 滑动角/° MJ7.3 32.74 130.76 235 88 −14 13.1 446 MJ5.4 32.75 130.76 234 74 15 15.0 2 2. 基本频谱分析
重大工程的抗震设计中，需要地震动时程作为地震输入来进行结构地震反应分析计算，而所输入地震动时程的质量和普适性与抗震分析结果的准确性直接相关（Dai et al，2014 ）。地震动时程的模拟结果则是衡量地震反演精准性的基本标准。本文所用程序为Miyake等（2003）修改后的egfm v1.2e程序。在此基础上，Li 等（2017，2018）对该程序的输出模式和内容进行调整，使其可直接输出模拟的波形数据和频谱数据，包括加速度时程、傅里叶谱、反应谱等数据。本文在地震动加速度时程数据的基础上计算得到熊本地震的地震动参数，包括47个模拟台站的PGA、阿里亚斯强度、卓越周期等参数，并将47个台站的数据进行拟合得到各个地震动参数的空间分布特征。将地震动基本频谱和参数的观测值与模拟值进行对比（图2）可以看出，各台站的模拟时程与观测时程拟合较好，从地震动持时及加速度的幅值上看，模拟的数据与目标台站观测数据总体上较为匹配。换言之，模拟时程与观测时程作用在建筑物上可产生相似甚至相同的破坏程度。此外，合成的地震动时程包络线与观测值的包络线具有较高的吻合度，模拟的时程可以较好地再现真实的地震动时程特征。
![]()
图 2 研究区内部分台站所记录的熊本地震观测值与模拟值的加速度时程图（a） FKO011台站；（b） FKO012台站；（c） KGS003台站；（d） KMM003台站；（e） MYZ005 台站Figure 2. The acceleration time history of the observed and synthetic of the Kumamoto earthquake recorded at several stations in the research area（a） Station FKO011；（b） Station FKO012；（c） Station KGS003；（d） Station KMM003；（e） Station MYZ005傅里叶谱的振幅描述了地震动过程的频谱特征，包括各频率分量的相位及幅值信息。本文中的傅里叶谱振幅的观测值和与模拟值整体上的相似性较好（图3），尤其当频率范围介于1—15 Hz之间时相似度极高；当频率低于0.5 Hz时，模拟的傅里叶谱的振幅一般较观测值的傅里叶谱的振幅低。这种现象也反映出经验格林函数法在模拟时具有高频段结果较好、而在低频段的模拟结果较差的特点。
![]()
图 3 熊本地震4个台站点位置的傅里叶谱Figure 3. The Fourier spectrum of four stations in Kumamoto earthquake3. 峰值加速度分布特征分析
峰值加速度是反映地震动特征的一个非常重要的参数，能够为地震烈度评估、地震危险区确定、城市未来规划以及建筑物的抗震设计提供科学的参考信息。本文通过对47个强震台站的PGA进行插值，得到了熊本地震影响区域内的PGA的分布情况，如图4所示。
通过图4对比可以看出，观测值与模拟值相似度较高；近场水平向PGA整体呈椭圆形向四周衰减扩散，长轴方向大致呈EW走向；EW方向分量的PGA最大值出现在KMM008台站，而NS方向分量PGA的最大值出现在KMM006台站。PGA呈椭圆形向四周衰减扩散的现象可能与断层的方向性效应和近地表断层上下盘断层效应有关。而在FKO015台站位置处，存在一个PGA值异常的区域，其模拟值比观测值要大，可能是由于台站与震源之间复杂地形的影响，但由于数值模拟时很难清晰地了解某些特殊地形对地震动的影响，故其中存在较大的不确定性。
![]()
图 4 熊本地震的PGA观测值与模拟值的空间分布特征（a） EW向观测值；（b） NS向观测值；（c） EW向模拟值；（d） NS向模拟值Figure 4. The distribution characteristics of observed and synthetic PGA record of Kumamoto earthquake（a） Observed EW component；（b） Observed NS component；（c） Synthetized EW component；（d） Synthetized NS component4. 阿里亚斯强度分布特征分析
阿里亚斯强度是表示地震能量的一个重要参量，描述了地震动能量的释放特征，包括振幅变化、频率时程等信息（Arias，1970），可以用于估计地震引起的滑坡、砂土液化以及结构破坏状态等（Arias，1970；Del Gaudio，Wasowski，2004），也可以应用于地震灾害预防以及震后地震快速救援等。
阿里亚斯强度可以量化地震发生过程中弹性系统的能量消耗，其表现形式为
$ {I_{\rm{a}}} {\text{＝}} \frac{\rm{\pi }}{{2{\rm{g}}}}\int_0^{{T_{\rm{d}}}} {a{{{(}t{)}}^2}} {\rm d}t {\text{，}} $
(1) 式中，a（t）为强震仪记录到的地震动加速度；Td为记录的持时；g为重力加速度（王秀英等，2010；李宗超，2017）。由于阿里亚斯强度比受干扰因素影响较多的峰值加速度参数能更准确地反映真正的地震动水平，因而其与地表破坏程度的关系也更为密切（王秀英等，2010）。通过对47个强震台站的阿里亚斯强度值进行插值拟合，得到整个研究区域的阿里亚斯强度分布情况，如图5所示。
![]()
图 5 熊本地震阿里亚斯强度观测值与模拟值空间分布特征（a） EW向观测值；（b） NS向观测值；（c） EW向模拟值；（d） NS向模拟值Figure 5. The distribution characteristics of Arias intensity of observed record and synthetic records for Kumamoto earthquake（a） Observed EW component；（b） Observed NS component；（c） Synthetized EW component；（d） Synthetized NS component同PGA的拟合结果相似，阿里亚斯强度的观测值与模拟值同样具备较高的吻合度，近场水平向整体呈椭圆形向四周衰减扩散，长轴方向大致呈EW走向。EW方向分量的阿里亚斯强度的最大值出现在KMM008台站，而NS方向分量的阿里亚斯强度的最大值出现在KMM006台站以及KMM009台站附近。阿里亚斯强度呈椭圆形向四周衰减扩散的现象可能与断层的方向性效应和近地表断层上下盘断层效应有关。
熊本地震的阿里亚斯强度的模拟值整体上要大于其观测值。EW方向分量差异较大，最大观测值为6.9 m/s，而最大模拟值为16.28 m/s；NS方向分量同样差异较大，最大观测值为6.87 m/s，而最大模拟值为10.54 m/s。阿里亚斯强度是与地震能量关联性较大的地震动参量，模拟值与观测值出现较大差异，可能是传播路径以及地形复杂因素所致，虽然经验格林函数法具有不考虑传播路径和场地条件影响的优势，但其只是个大致的假定，真实情况下，大、小地震的传播路径并非完全一致，存在较大的不确定性。
5. 卓越周期分布特征分析
卓越周期的本质是波形的共振，也就是说当地震波振动周期与岩石或土体的自然振动周期一致时，共振效应所引起的表面振动会加强（李宗超，2017）。当然如果某一周期（与场地卓越周期不一致）的地震动成分很高，即使在复杂传播路径及场地环境的影响下，这一成分的地震动反应仍会比较强烈，但是在当前工程应用领域中，对建筑物抗震设防影响较大的地震动周期成分和场地的卓越周期是一致的。
卓越周期是通过地震记录统计得到的，地基土硬度不同则所对应的卓越周期不同。根据卓越周期的不同，地基土可以分为4个级别，如表2所示。土越软越深，频谱中的长周期成分就越多。根据熊本地震的地震动数据可以得出该地震中S波的卓越周期，同时与模拟值进行对比，结果如图6所示。
表 2 根据卓越周期进行的地基土分类Table 2. Foundation soil classification based on predominant period序号 类别 卓越周期/s 平均周期/s 1 稳定基岩层 0.1—0.2 0.15 2 普通土层 0.21—0.4 0.27 3 软土层 介于第2类土与第4类土之间 介于第2类土与第4类土之间 4 异常土层 0.3—0.7 0.5 ![]()
图 6 熊本地震的卓越周期观测值与模拟值的空间分布特征（a） EW向观测值；（b） NS向观测值；（c） EW向模拟值；（d） NS向模拟值Figure 6. The distribution characteristics of predominant period for observed record and synthetic record（a） Observed EW component；（b） Observed NS component；（c） Synthetized EW component；（d） Synthetized NS component卓越周期整体上拟合较好，但不同方向分量上存在差异。在EW方向分量上，地震动观测值的卓越周期最大值为0.78 s，出现在KMM013台站，模拟值卓越周期最大值为0.54 s，出现在KGS004台站。卓越周期在近场、震源西北部、东部等区域拟合较好，而在震源南部区域拟合欠佳，地震动模拟值的卓越周期明显比观测值的卓越周期大，且在台站MYZ010和台站MYZ012的模拟结果明显比周围台站的大，与观测值的分布特征存在明显差异。
对于NS方向分量，观测值卓越周期最大值为0.72 s，出现在KMM010台站，模拟值卓越周期最大值为0.66 s，出现在KMM010台站，最大观测值与模拟值在同一位置出现。但在震源东部，例如MYZ001，MYZ019和OIT015台站，模拟值的卓越周期明显小于观测值的卓越周期，观测值在这些台站处的卓越周期异常值并未体现出来。在震源南部区域，模拟值的卓越周期与观测值的卓越周期也存在较大差异。
卓越周期是与场地条件相关性很大的参数，与场地土层类型以及场地周边地形特征等都有关系。经验格林函数方法重点考虑了震源环节的参数选取，对于场地则假定大小地震的场地具有相似性，此假设存在不确定性，这是导致卓越周期在某些区域的拟合结果差异较大，因此，后期在用经验格林函数法模拟卓越周期时，尤其是在地形以及场地条件复杂的区域，需要补充考虑场地条件的影响。
6. 讨论与结论
本文比较了地震动时程的观测值与模拟值的基本频谱，结果显示：模拟的地震动持时以及加速度幅值与目标台站总体上较为匹配；模拟的地震动时程包络线与观测值的包络线有较高的吻合度，模拟的时程可以较好地再现真实的地震动时程特征；振幅谱的观测值与模拟值整体上的相似性非常好，尤其是频率范围介于1—15 Hz之间的频段，相似度极高，但是，在频率低于0.5 Hz时，模拟的振幅谱一般比观测值的振幅谱低，这也反映出经验格林函数法对于高频段的模拟较好，在低频段的模拟结果不是很好。
通过对熊本地震地震动工程应用参数的观测值与模拟值的比较和分析得到如下结论：
1） PGA与阿里亚斯强度的观测值与模拟值在震源距小于50 km范围内的拟合结果较好，阿里亚斯强度的模拟值整体大于观测值，且二者均在近场区域以椭圆形向周围衰减扩散；
2） 卓越周期整体拟合较好，但在局部区域存在较大差异，模拟结果难以再现场地环境的复杂性。
本文也再次验证了使用经验格林函数模拟地震动是可行的。但也不否认，一些参数在一些地区的观测值与模拟值存在差异，甚至是较大的差异，并非所有参数都拟合得很好。造成这种差异的原因很多，例如传播路径和传播介质的不均匀特性，复杂地形引起来的地震波的反射和折射，模拟的过程还包括其它很多不确定性。考虑不确定性因素进行地震动模拟将是我们下一步需要重点探索的问题。
-
![]()
图 3 研究区布格重力异常（a）和重力异常（b）分布图
Figure 3. Bouguer gravity anomaly （a） and regional gravity anomaly （b） around Ordos block
![]()
图 5 鄂尔多斯地块及邻区的区域重力异常功率谱
Figure 5. Power spectrum of regional gravity anomaly in Ordos block and its surrounding regions
![]()
图 6 研究区的区域重力异常拟合情况对比图
（a） 理论区域重力异常；（b） 理论重力异常与实际重力异常间的残差（RES）分布
Figure 6. Comparison of simulative gravity anomalies with real ones
（a） Simulated gravity anomaly based on the theory model；（b） Residuals between simulative gravity anomalies based on the theory model and the real ones
![]()
图 7 研究区内壳幔密度差分布图
Figure 7. Density difference between crust and mantle in the studied area
![]()
图 8 研究区变密度差（a）和常密度差（b）对应的重力均衡状态
Figure 8. Gravity isostatic status from variable crust-mantle density difference （a） and from constant density difference （b） in the studied region
-
陈石,王谦身,徐伟民,蒋长胜,卢红艳. 2011. 华北地区热均衡、重力均衡与深部构造[J]. 地球物理学报,54(11):2864–2875. doi: 10.3969/j.issn.0001-5733.2011.11.016 Chen S,Wang Q S,Xu W M,Jiang C S,Lu H Y. 2011. Thermal isostasy of North China and its gravity isostasy and deep structure[J]. Chinese Journal of Geophysics,54(11):2864–2875 (in Chinese).
陈兆辉,王椿镛,楼海. 2018. 鄂尔多斯地块地壳上地幔速度结构及构造意义[J]. 科学通报,63(3):327–339. Chen Z H,Wang C Y,Lou H. 2018. Crust and upper mantle velocity structures beneath the Ordos blocks and its tectonic significance[J]. Chinese Science Bulletin,63(3):327–339 (in Chinese). doi: 10.1360/N972017-00741
方剑,马宗晋,许厚泽. 2006. 地形-均衡补偿重力、大地水准面异常频谱分析[J]. 地球物理学进展,21(1):25–30. doi: 10.3969/j.issn.1004-2903.2006.01.005 Fang J,Ma Z J,Xu H Z. 2006. Frequency analyze of topography-isostasy compensative gravity and geoid anomaly[J]. Progress in Geophysics,21(1):25–30 (in Chinese).
方盛明,冯锐,田长征,孙桂香,王志理,李长法. 1997. 亚洲中部地区均衡重力异常特征及地震活动性[J]. 地震学报,19(6):650–654. Fang S M,Feng R,Tian C Z,Sun G X,Wang Z L,Li C F. 1997. Equilibrium compensation depth and inversion analyses for elastic slab model of the lithosphere in the central Asian region[J]. Acta Seismologica Sinica,19(6):650–654 (in Chinese).
冯娟,孟小红,陈召曦,石磊,武粤,樊振军. 2014. 三维密度界面的正反演研究和应用[J]. 地球物理学报,57(1):287–294. Feng J,Meng X H,Chen Z X,Shi L,Wu Y,Fan Z J. 2014. The investigation and application of three-dimensional density interface[J]. Chinese Journal of Geophysics,57(1):287–294 (in Chinese).
付广裕. 2014. 川西盆地重力异常与地壳均衡研究[C]//中国地球科学联合学术年会. 北京: 中国地球物理学会: 1729-1732. Fu G Y. 2014. Research on gravity anomaly and Moho isostasy of western Sichuan[C]//Annual Meeting of Chinese Geoscience Union. Beijing: Chinese Geophysical Society: 1729-1732 （in Chinese）.
郭良辉,孟小红,石磊,陈召曦. 2012. 优化滤波方法及其在中国大陆布格重力异常数据处理中的应用[J]. 地球物理学报,55(12):4078–4088. doi: 10.6038/j.issn.0001-5733.2012.12.020 Guo L H,Meng X H,Shi L,Chen Z X. 2012. Preferential filtering method and its application to Bouguer gravity anomaly of Chinese continent[J]. Chinese Journal of Geophysics,55(12):4078–4088 (in Chinese).
胡敏章,李建成. 2010. 顾及地壳密度模型的Airy-Heiskanen均衡重力异常的计算[J]. 大地测量与地球动力学,30(5):48–52. Hu M Z,Li J C. 2010. Computation of Airy-Heiskanen isosatic gravity anomaly with considering crust density model[J]. Journal of Geodesy and Geodynamics,30(5):48–52 (in Chinese).
黎哲君,徐如刚,张玮晶,张毅,李辉,杨光亮,龙剑锋. 2015. 华北断块区均衡重力异常、构造特征及地震活动[J]. 地震地质,37(2):496–509. doi: 10.3969/j.issn.0253-4967.2015.02.012 Li Z J,Xu R G,Zhang W J,Zhang Y,Li H,Yang G L,Long J F. 2015. Isostatic gravity anomaly,structural features and seismicity of North China craton[J]. Seismology and Geology,37(2):496–509 (in Chinese).
楼海,王椿镛,刘琼林,姚志祥. 2012. 青藏高原东北缘和鄂尔多斯地区的均衡重力异常[J]. 国际地震动态,(6):72. Lou H,Wang C Y,Liu Q L,Yao Z X. 2012. Isostatic gravity anomalies of Ordos and northeast margin of the Tibetan Plateau[J]. Recent Developments in World Seismology,(6):72 (in Chinese).
任枭,徐志国,杨辉,陈宏峰,邹立晔. 2012. 鄂尔多斯地块东南缘地带Moho深度变化特征研究[J]. 地球物理学报,55(12):4089–4096. doi: 10.6038/j.issn.0001-5733.2012.12.021 Ren X,Xu Z G,Yang H,Chen H F,Zou L Y. 2012. Moho depth distribution character beneath the Ordos block’s southeastern margin areas[J]. Chinese Journal of Geophysics,55(12):4089–4096 (in Chinese).
唐新功. 2006. 华北地区重力学研究[D]. 北京: 北京大学: 1-34. Tang X G. 2006. Gravity Study of North China Region[D]. Beijing: Peking University: 1-34 （in Chinese）.
王谦身,滕吉文,张永谦,杨辉. 2008. 龙门山断裂系及邻区地壳重力均衡效应与汶川地震[J]. 地球物理学进展,23(6):1664–1670. Wang Q S,Teng J W,Zhang Y Q,Yang H. 2008. The effect of crustal gravity isostasy and Wenchuan earthquake in Longmenshan faults and adjacent area[J]. Progress in Geophysics,23(6):1664–1670 (in Chinese).
王谦身,滕吉文,张永谦,皮娇龙. 2015. 鄂尔多斯—中秦岭—四川东部的重力异常场与深部地壳结构[J]. 地球物理学报,58(2):532–541. Wang QS,Teng J W,Zhang Y Q,Pi J L. 2015. Gravity anomalies and deep crustal structure of the Ordos basin-middle Qinling orogeny-eastern Sichuan basin[J]. Chinese Journal of Geophysics,58(2):532–541 (in Chinese).
王谦身,滕吉文,张永谦,杨辉,赵彬彬,胡国泽. 2013. 中秦岭地带重力异常特征及地壳结构的探榷[J]. 地球物理学报,56(12):3999–4008. doi: 10.6038/cjg20131206 Wang Q S,Teng J W,Zhang Y Q,Yang H,Zhao B B,Hu G Z. 2013. Discussion on gravity anomalies and crustal structure of the middle Qinling Mountains[J]. Chinese Journal of Geophysics,56(12):3999–4008 (in Chinese).
许德树,曾华霖. 2000. 优选延拓技术及其在中国布格重力异常图处理上的应用[J]. 现代地质,14(2):215–222. doi: 10.3969/j.issn.1000-8527.2000.02.016 Xu D S,Zeng H L. 2000. Preferential continuation and its application to Bouguer gravity anomaly in China[J]. Geoscience,14(2):215–222 (in Chinese).
杨彦明,张国清,陈婧,黄瑞滨. 2016. 鄂尔多斯块体北缘及邻区Moho面深度特征[J]. 地震地磁观测与研究,37(6):1–8. Yang Y M,Zhang G Q,Chen J,Huang R B. 2016. Moho depth distribution character beneath the Ordos block’s northern margin areas and surrounding regions[J]. Seismological and Geomagnetic Observation and Research,37(6):1–8 (in Chinese).
殷秀华,刘铁胜,刘占坡. 1993. 均衡重力异常和地壳表、浅层地质结构[J]. 地震地质,15(2):149–156. Yin X H,Liu T S,Liu Z P. 1993. Isostatic anomalies of gravity and geologic structures of the surface and shallow crust[J]. Seismology and Geology,15(2):149–156 (in Chinese).
詹艳. 2008. 青藏高原东北缘地区深部电性结构及构造涵义[D]. 北京: 中国地震局地质研究所: 1–59. Zhan Y. 2008. Deep Electric Structures Beneath the Northeastern Margin of the Tibetan Plateau And Its Tectonic Impli-cations[D]. Beijing: Institute of Geology, China Earthquake Administration: 1–59 （in Chinese）.
张恩会. 2015. 基于人工地震与重力约束的川滇地区地壳结构研究[D]. 北京: 中国地震局地球物理研究所: 7–48. Zhang E H. 2015. Crustal Structure Across Sichuan and Yunnan Region From Wide-Angle Seismic And Gravity Data[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 7–48 （in Chinese）.
张永谦,王谦身,滕吉文. 2010. 川西藏东地区的地壳均衡异常及其与地震分布的关系[J]. 地球物理学报,53(11):2631–2638. Zhang Y Q,Wang Q S,Teng J W. 2010. The crustal isostatic anomaly beneath eastern Tibet and western Sichuan and its relationship with the distribution of earthquakes[J]. Chinese Journal of Geophysics,53(11):2631–2638 (in Chinese).
赵国泽,詹艳,王立凤,王继军,汤吉,陈小斌,肖骑彬. 2010. 鄂尔多斯断块地壳电性结构[J]. 地震地质,32(3):345–359. doi: 10.3969/j.issn.0253-4967.2010.03.001 Zhao G Z,Zhan Y,Wang L F,Wang J J,Tang J,Chen X B,Xiao Q B. 2010. The active fault belts in eastern Tibet margin inferred using magneto-tellurics[J]. Seismology and Geology,32(3):345–359 (in Chinese).
Airy G B. 1855. On the computation of the effect of the attraction of mountain-masses,as disturbing the apparent astronomical latitude of stations in geodetic surveys[J]. Philos Trans R Soc London,16(2):42–43.
Bagherbandi M. 2011. An Isostatic Earth Crustal Model and Its Applications[D]. Sweden: KTH Royal Institute of Technology.
Dorman L M,Lewis B T R. 1970. Experimental isostasy:Theory of the determination of the Earth’s isostatic response to a concentrated load[J]. J Geophys Res,75(17):3357–3365. doi: 10.1029/JB075i017p03357
Dutton C E. 1889. On some of the greater problems of physical geology[J]. Bull Phil Soc Wash,11:51–64.
Guo L H, Meng X H, Chen Z X,Li S L,Zheng Y M. 2013. Preferential filtering for gravity anomaly separation[J]. Comput Geosci,51(2):247–254.
Jarvis A, Reuter H I, Nelson A, Guevara E. 2008. Hole-filled seamless SRTM data V4, International Centre for Tropical Agriculture （CIAT）[EB/OL]. [2018–01–03]. http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp .
Li Y H,Gao M T,Wu Q J. 2014. Crustal thickness map of the Chinese mainland from teleseismic receiver functions[J]. Tectonophysics,611:51–60. doi: 10.1016/j.tecto.2013.11.019
Parker R L. 1973. The rapid calculation of potential anomalies[J]. Geophys J Int,31(4):447–455. doi: 10.1111/j.1365-246X.1973.tb06513.x
Pratt J H. 1855. On the attraction of the Himalaya mountains,and of the elevated regions beyond them,upon the plumb-line in India[J]. Philos Trans R Soc London,145:53–100. doi: 10.1098/rstl.1855.0002
Sandwell T,Smith W H F. 2009. Global marine gravity from retracked Geosat and ERS-1 altimetry:Ridge segmentation versus spreading rate[J]. J Geophys Res,114(B1):B01411.
Steffen R,Steffen H,Jentzsch G. 2011. A three-dimensional Moho depth model for the Tien Shan from EGM2008 gravity data[J]. Tectonics,30(5):TC5019. doi: 10.1029/2011TC002886
-
期刊类型引用(5)
1. 李宗超，纪志伟，孙吉泽，陈鲲，谢俊举. 基于安徽合肥M_S4.7地震快速估计本地区M_S6.0设定地震加速度时程及其破坏能力. 地球与行星物理论评(中英文). 2025(03): 338-351 . 
百度学术
2. 黄亭，李宗超，唐方头，孙吉泽，付长华，王晓辉，纪志伟，高孟潭，叶坤. 泸定地震大渡河区域灾害调查及格林函数库应用可行性探讨. 地球与行星物理论评(中英文). 2024(04): 440-452 . 百度学术
3. 李宗超，孙吉泽，高孟潭，陈学良，赵延娜，吴清. 青海玛多M7.4地震中野马滩大桥地震动特征初判. 地球与行星物理论评. 2022(01): 101-106 . 百度学术
4. 李宗超，高孟潭，孙吉泽，司洁戈，吴清，李奇. 海域俯冲带高频地震动模拟——以2021年2月13日日本福岛M_S7.1地震为例. 震灾防御技术. 2022(03): 516-528 . 百度学术
5. 罗全波，陈学良，高孟潭，李铁飞. 台湾双冬断层近场脉冲型地震动的数值模拟. 地震学报. 2019(03): 377-390+411 . 本站查看
其他类型引用(2)
下载:
计量
- 文章访问数: 1285
- HTML全文浏览量: 748
- PDF下载量: 71
- 被引次数: 7
