华北地区岩石圈有效弹性厚度及其地球动力学意义.pdf

第4 3卷 第9期2 0 2 3年9月大 地 测 量 与 地 球 动 力 学J o u r n a l o fG e o d e s ya n dG e o d y n a m i c sV o l.4 3N o.9S e p t.,2 0 2 3收稿日期:2 0 2 2-1 1-0 1项目来源:国家自然科学基金(4 2 2 7 4 0 8 7,4 1 8 7 4 1 1 9,4 2 1 7 4 0 8 3)。第一作者简介:孙文魁,硕士生,主要从事重力学和地球动力学研究,E-m a i l:2 0 2 0 7 1 3 1 4y a n g t z e u.e d u.c n。D O I:1 0.1 4 0 7 5/j.j g g.2 0 2 3.0 9.0 1 1文章编号:1 6 7 1-5 9 4 2(2 0 2 3)0 9-0 9 3 6-0 9华北地区岩石圈有效弹性厚度及其地球动力学意义孙文魁1,2 唐新功1,21 长江大学油气资源与勘探技术教育部重点实验室,武汉市大学路1 1 1号,4 3 0 1 0 02 长江大学地球物理与石油资源学院,武汉市大学路1 1 1号,4 3 0 1 0 0摘 要:以E I G E N-6 C 4模型布格重力异常数据、V 1 9.1地形数据和C R U S T 1.0模型地壳数据为基础,采用导纳法和相关函数法联合反演,计算华北地区岩石圈有效弹性厚度(Te)。结果表明,华北地区Te值的分布存在明显的横向不均匀性,西部鄂尔多斯地块为高值区,Te厚度约为3 06 0k m;中部华北平原为低值区,Te厚度约为52 5k m;鲁东黄海地块为中、高值区,Te厚度约为2 05 0k m。华北克拉通高Te值区与稳定的鄂尔多斯地块相对应,中、低Te值区则对应构造相对活跃的太行山构造带和郯庐断裂带。从Te值与地表热流的关系来看,高Te值区一般对应低热流值区,而低Te值区一般热流值也较高,这与岩石圈厚度和地表热流的一般分布规律一致。从Te与地震的空间分布来看,华北地区的地震活动主要集中在中、低Te值区,表明Te较低的地区稳定性较差,可能更有利于地震的孕育和发生。关键词:岩石圈;有效弹性厚度;华北克拉通;地震活动;地表热流中图分类号:P 3 1 3 文献标识码:A 岩石圈的力学强度反映了岩石圈长期受构造应力作用的响应,是影响岩石圈变形演化的主要因素,表征岩石圈力学强度的最佳方式之一是通过计算岩石圈有效弹性厚度(e f f e c t i v ee l a s t i ct h i c k n e s so f t h e l i t h o s p h e r e,Te)来实现。Te定义为与岩石圈板块中实际应力分布所产生的弯矩相等的弯曲弹性薄板的厚度1。Te概念源于区域均衡补偿模型2或挠曲模型、弹性板模型。不同于局部补偿模型3-4(如A i r y模型和P r a t t模型),区域均衡模型或挠曲模型其载荷部分是由覆盖在软流圈5之上的岩石圈板内的弹性应力支撑并在一定区域范围内发生补偿。在挠曲模型中,岩石圈弹性板的响应往往使用其挠曲刚度(D)来表征。挠曲模型与A i r y模型类似,均是通过岩石圈向下(或向上)发生挠曲以及地壳增厚(或减薄)来实现对地表地形负载的补偿。不同的是,作为局部补偿模式的A i r y模型对应的挠曲刚度为0,即代表弱地壳;而挠曲模型则代表具有一定横向强度的地壳6。基于地形和重力异常数据求取Te的导纳方法最早由D o r m a n等7提出,但在早期利用导纳法(a d m i t t a n c e)估计Te时均是仅考虑地表的负载,而忽略地下负载加载的情况。实际上岩石圈不仅会受到地表负载的施压发生形变而产生山根,其同样可能受到地下负载的加载而产生上隆,例如洋中脊与地幔热物质的上涌会造成地壳和地幔物质的密度变化,从而形成地下负载8。F o r-s y t h1首次提出基于地形和布格重力异常的相关函数法(c o h e r e n c e),引入地表和地下负载的概念,并在导纳法中率先加入地下负载导纳的计算方法。J i n等9在F o r s y t h基础上又提出准相关理论,该理论的关键是上下负载的相关角,当=9 0 时,地表和地下负载正交,相关函数为传统的相关函数;当=0 时,地表和地下负载完全相关,相关函数为1;当0 9 0 时,地表和地下负载具有不同的相关度。由于采用传统傅里叶变换求取地形和重力异常功率谱方法(称为周期图法)会存在频率泄露问题,文献1 0-1 1 又分别将多窗谱法引入导纳法和相关函数法以改善此问题。随后,K i r b y等1 2基于叠加的M o r l e t小波,提出一种F a n小波相关函数法来反演岩石圈的有效弹性厚度与初始负载比分布,该方法可克服采用固定尺度窗口而造成的 第4 3卷第9期孙文魁等:华北地区岩石圈有效弹性厚度及其地球动力学意义数据频域信息不完整的缺陷1 3。华北地区位于欧亚大陆东南部,自中生代以来,相对稳定的华北克拉通遭受强烈的活化改造,发生大规模的构造变形、岩石圈减薄、岩浆活动和盆地再造,致使其原有的克拉通结构遭到明显的破坏1 4,沉积层结构也相当复杂1 5(图1),因此华北克拉通被认为是世界上研究古老岩石圈减薄、改造、置换的最佳场所和天然实验室1 6。近年来,国内外学者围绕华北克拉通破坏先后提出过不同的动力学模式,主要有热模式、化学侵蚀模式和拆沉模式等。然而无论哪种模式,都存在一些难以回答的问题,只能够解释部分地质现象1 7。华北地区活动断裂丰富,地震频发,历史上曾发生过多次7级以上强震,因此认识华北地区的深部构造特征及岩石圈动力演化对于认识华北岩石圈力学强度和克拉通破坏具有重要的科学意义。部分学者利用不同的数据和方法对华北地区的岩石圈有效弹性厚度进行估算,如郑勇等1 8采用F a n小波分析法计算华北克拉通Te值及其各向异性;陈波等6采用多窗谱法和F a n小波分析法计算中国及邻区Te值,并分析其动力学意义;L u等1 9和侍文等2 0分别基于导纳法和相关函数法,以E T O P O 1地形数据、WGM 2 0 1 2布格重力异常数据和C r u s t 1.0地壳结构数据为基础计算中国大陆地区Te值分布,分别对中国大陆强震构造区力学特征和影响岩石圈强度的因素进行探讨。由于华北地区的地下结构非常复杂,前人使用的数据和方法存在差异,导致对华北地区Te的计算结果不尽相同,同时部分学者计算时也未考虑沉积层的影响。因此,本文使用精度更高的重力异常和地形数据,在完成对重力异常数据的沉积层校正后,基于导纳法和相关函数法联合反演并结合贝叶斯最优参数估计方法,得到华北地区的岩石圈有效弹性厚度以及初始负载比,讨论岩石圈有效弹性厚度与地震活动分布以及地表热流分布的关系,以期为认识华北岩石圈力学特征和克拉通破坏机制提供重要的重力学参考。1 方法及数据1.1 导纳法与相关函数法联合反演方法导纳法的优点是物理意义明确、易于理解和实现;而相关函数法对不同均衡响应模式的转换波长更为敏感,计算精度更高2 0。A u d e t2 1提出结合导纳法和相关函数法的联合反演方法,该方法被国内外学者广泛采用并应用于Te值的估算2 0-2 2。本文采用基于连续平面小波变换计算得到的图1 华北地区地质构造F i g.1 G e o l o g i c a l s t r u c t u r em a po fn o r t hC h i n a导纳和相关函数进行联合反演,估计华北地区Te值,同时为判断结果的准确性,结合贝叶斯推断法来评价特定情况下模型参数的统计分布。导纳和相关函数的表达式分别为1,2 1:Z(k)=g(k)H(k)(1)r2(k)=|2(2)式中,Z(k)为导纳,g(k)为重力异常谱,H(k)为地形谱,k为波数,表示平均谱,*表示共轭。为更准确地求取岩石圈有效弹性厚度,首先需要评价和判断负载来源,即首先考虑地表负载和地下负载同时加载的情况(图2)。图中上部浅色区域代表地壳,深色区域代表地幔,Hi和Wi分别为初始地表和地下负载,HT和WT表示在地表负载作用下产生的地表和莫霍面变形,HB和WB表示在地下负载作用时产生的地表和莫霍面变形,H和W表示在共同作用下最终的地表和莫霍面变形。图2 弹性板挠曲响应示意图1F i g.2 F l e x u r a l r e s p o n s ed i a g r a mo f e l a s t i cp l a t e1739大 地 测 量 与 地 球 动 力 学2 0 2 3年9月采用发生挠曲变形前地表负载与地下负载的初始负载比F来判断负载的主要来源2 3-2 4,初始负载比定义为:F=f1+f(3)式中,f表示地下负载产生的高程谱HB与地表地形产生的高程谱HT之间的负载系数(l o a d i n gf a c t o r)。当F=0时,f为无穷小,代表负载完全来自于地表;F=1时,f为无穷大,代表负载完全来自于地下。使用连续平面小波变换分别计算重力异常谱和地形谱,得到每个网格节点的导纳和相干曲线。该方法的优点是可以避免将网格划分为一系列重叠窗口,从而大大降低频谱泄漏的概率。设定好一系列初始模型参数并计算其理论导纳和相关性之后,再使用非线性最小二乘法搜索与实测导纳和相关性拟合效果最佳的Te和F值来实现对正演数据的反演。理论与实测导纳或相关函数误差使用如下卡方函数(C h i-s q u a r e df u n c t i o n)进行计算2 1:2(F,Te)=1(2N-2)2j=1Ni=1di j-si j(F,Te)i j 2(4)式中,j表示实际导纳(j=1)或实测相关性(j=2),d表示实测导纳或相关性,为理论与实测导纳或相关性函数之间的方差,s为给定模型的正演结果,N为采样网格数量。1.2 贝叶斯推断使用贝叶斯推断方法2 5在每个采样网格对导纳法或相关函数法获得的有效弹性厚度(Te)与载荷比(F)数据模型M()进行评估,以确定该网格 最 优 的 由 Te,F组 成 的 模 型 参 数 组合1 9-2 2。贝叶斯推断方法表达式为:P(|x)=P(x|)P()P(x)=P(x|)P()P(x|)P()d(5)式中,x代表采样网格数据点构成的集合x1,x2,xn;为采样网格数据点所对应的模型参数,由Te,F 组成;P(|x)为给定一组数据x后得到的模型参数的概率(即后验);P(x|)为给定参数时生成数据x的概率(即似然);P()为独立于数据x的的分布(称为先验);P(x)为归一化因子,以保证在整个模型参数空间上评估的后验概率为1。从贝叶斯推断方法中随机选择3处区域来显示模型参数统计分布特征和曲线拟合效果,包括华北克拉通西部(1 0 9.2 E,3 7.8 N)(图3和图4)、中部(1 1 4.1 E,3 8.4 N)(图5和图6)以及东部地区(1 1 8.6 E,3 7.5 N)(图7和图8),使用贝叶斯推断法得到的模型参数统计分布以及采用不同反演方法得到的导纳和相关性拟合曲线如图38所示,其中每处均使用4个独立的马尔科夫链(MCMC)从后验序列中抽取20 0 0个样本,并忽略前5 0 0个样本。通过对比相关函数法和联合反演法结果,以展示联合反演方法的优势。3个区域分别使用相关函数法和联合反演法得到的平均值、拟合误差和最大后验估计值见表1。图中蓝色符号为观测导纳与相关性曲线图3 华北克拉通西部不同方法反演的导纳与相关性曲线拟合F i g.3 A d m i t t a n c ea n dc o h e r e n c e f i t t i n go fd i f f e r e n ti n v e r s i o nm e t h o d s i nt h ew e s t e r nn o r t hC h i n aC r a t o n图中为平均值,为标准差,9 5%C I对应9 5%置信区 间,MA P为最大后验估计图4 华北克拉通西部贝叶斯推断法模型参数统计分布F i g.4 S t a t i s t i c a l d i s t r i b u t i o no fm o d e l p a r a m e t e r sf r o mt h eB a y e s i a n i n f e r e n c em e t h o d i nt h ew e s t e r nn o r t hC h i n aC r a t o n从表1以及图4、图6和图8这3组计算结果可以看出,仅使用相关函数法进行反演时,得到的839 第4 3卷第9期孙文魁等:华北地区岩石圈有效弹性厚度及其地球动力学意义图5 华北克拉通中部不同方法反演的导纳和相关性曲线拟合F i g.5 A d m i t t a n c ea n dc o h e r e n c e f i t t i n go fd i f f e r e n ti n v e r s i o nm e t h o d s i nt h ec e n t r a l n o r t hC h i n aC r a t o n图6 华北克拉通中部贝叶斯推断法模型参数统计分布F i g.6 S t a t i s t i c a l d i s t r i b u t i o no fm o d e l p a r a m e t e r sf r o mt h eB a y e s i a n i n f e r e n c em e t h o d i nt h ec e n t r a l n o r t hC h i n aC r a t o n图7 华北克拉通东部不同方法反演的导纳和相关性曲线拟合F i g.7 A d m i t t a n c ea n dc o h e r e n c e f i t t i n go fd i f f e r e n ti n v e r s i o nm e t h o d s i nt h ee a s t e r nn o r t hC h i n aC r a t o nTe和F值不确定性更大,最大后验估计MA P远大于平均值,Te值高于联合反演法,标准差也明显大于联合反演法。从曲线拟合效果来看(图图8 华北克拉通东部贝叶斯推断法模型参数统计分布F i g.8 S t a t i s t i c a l d i s t r i b u t i o no fm o d e l p a r a m e t e r sf r o mt h eB a y e s i a n i n f e r e n c em e t h o d i nt h ee a s t e r nn o r t hC h i n aC r a t o n3、图5和图7),采用联合反演法得到的导纳曲线与相关性曲线拟合效果更佳,高于单独使用相关函数法的反演效果。因此,不论Te值高低,相比于单独使用相关函数法,联合反演法均能获得更加准确的Te和F结果。表1 不同方法的模型参数统计T a b.1 M o d e lp a r a m e t e r s t a t i s t i c so fd i f f e r e n tm e t h o d s位置方法参数平均值标准差最大后验估计MA PTe/k mFTe/k mFTe/k mF西部相关函数法4 5 0.4 62 0 0.3 51 2 6 0.9 9联合反演法4 0 0.7 71 0.0 24 0 0.7 7中部相关函数法1 8 0.0 75 0.0 71 2 7 1.0 0联合反演法1 1 0.3 51 0.0 71 1 0.3 4东部相关函数法3 4 0.0 67 0.0 31 2 6 0.9 8联合反演法2 4 0.1 41 0.0 22 4 0.1 41.3 重力数据本文以华北地区(1 0 7 1 2 3 E,3 0 4 2 N)为研究区,以最新的基于卫星、地面和航空等重力数据融合的高精度全球重力模型数据E I G E N-6 C 42 6为基础,采用导纳法和相关函数法进行联合反演,求取华北地区岩石圈有效弹性厚度Te。E I G E N-6 C 4模型融合了G O C E卫星重力数据,在长波长部分较E GM 2 0 0 8模型数据具有更高的精度2 7-2 8。地形数据来自斯克里普斯海洋研究所(S I O)的V 1 9.1模型,空间分辨率为1 1。莫霍面深度数据和地壳密度数据来自C R U S T 1.0模型2 9。使用939大 地 测 量 与 地 球 动 力 学2 0 2 3年9月谱方法进行计算时,需要地形和重力数据在长波长部分具备良好的相关性,但短波长的沉积层会降低地形的长波长成分,导致地形数据与重力异常之间的相关性降低,在计算Te过程中会引入较大噪声,影响Te值的准确求取。因此本文根据C R U S T 1.0模型中沉积层厚度数据,采用P a r k e r3 0正演法计算沉积层产生的重力异常,然后从布格重力异常中减去其结果,从而实现重力异常数据沉积层校正。2 反演结果图9为反演得到的华北地区有效弹性厚度Te与初始负载比F分布图。由图9(a)可以看出,华北地区Te值分布存在明显的横向差异性,整体呈现东西高、中间低,北部高、南部低特征。鄂尔多斯盆地Te值较高,大约为3 06 0k m;华北平原Te值总体较低,约为52 5k m;鲁东-黄海地块Te值较高,约为2 0 5 0k m。初始负载比F(图9(b)结果显示,华北地区东南部和太行山构造带北部F值较高(大于0.8),鄂尔多斯地区F介于0.4和0.6之间,其他地区特别是华北平原地区整体F值较低(小于0.4),表明在华北地区东南部和太行山构造带负载主要来自于地表,而在华北盆地地区负载主要来自于地下,这也与华北盆地岩石圈拉张减薄具有较好的对应关系。这种明显的不均匀性证实了华北克拉通东部发生克拉通破坏和岩石圈减薄,而西部仍保持相对稳定的特征。其原因可能是自中生代以来太平洋板块俯冲和停滞,地幔过渡带俯冲板片的停滞和脱水造成含水量增加,显著降低岩石圈强度,引起上地幔减压熔融或地幔物质向上流动。这些热物质上升后,不断通过热作用和化学作用从底部侵蚀岩石圈,从而逐渐导致华北克拉通岩石圈的减薄和破坏。为考察图9结果的准确性,计算有效弹性厚度Te和初始负载比F的误差,结果如图1 0所示。可以看出,两者误差分别介于05k m和00.1 2,且大部 分 区 域 的 误 差 分 别 在03k m和00.0 6之内,表明本文得到的Te和F结果均能较好地反映华北地区岩石圈力学强度及影响岩石圈力学强度的动力学来源。侍文等2 0和L u等1 9采用联合反演计算得到华北地区Te值在55 0k m之间,与本文计算结果基本一致,且高低值分布区也较为吻合,说明本文计算结果的可靠性较高。图中灰色圆点表示1 9 0 02 0 2 2年发生在华北地区5级及以上地震图9 华北地区Te值和F分布F i g.9 D i s t r i b u t i o no fTea n dFi nn o r t hC h i n a3 结果分析3.1 Te分布与华北克拉通破坏从图9可以看出,华北克拉通西部鄂尔多斯地区表现为Te高值区(3 06 0k m),表明其岩石圈强度较大,构造相对稳定;F值介于0.40.6之间,表明鄂尔多斯地区的地表负载和地下负载大致相当。前人研究结果表明,鄂尔多斯地块上地幔呈现为明显的高速异常3 1,说明鄂尔多斯地块仍然是一个未遭破坏、相对稳定的克拉通单元。南部秦岭-大别造山带Te值较低,岩石圈活动性较强,秦岭-大别造山带东段F值较高,表明该地区岩石圈强度主要受地下负载影响。推测是扬子板块对华北地块的深向俯冲,强烈挤压隆起和碰撞造山等复杂构造演化,导致秦岭-大别造山带东段内部出现地壳垂向错断,并且秦岭-大别造山带东段与周围地块之间存在明显的深部构造转换带。在发生克拉通破坏时,东西向延伸的秦岭-大别-苏鲁局部增厚岩石圈将来自于地幔深处的热量转流到北部受扬子板块深俯冲作用影响较小的049 第4 3卷第9期孙文魁等:华北地区岩石圈有效弹性厚度及其地球动力学意义图1 0 华北地区Te误差和F误差分布F i g.1 0 D i s t r i b u t i o no fTee r r o ra n dFe r r o r i nn o r t hC h i n a冀中、渤海湾以及松辽地区的岩石圈底部,使得这些区域的岩石圈成为减薄中心1 6。华北克拉通中部呈现为沿太行山构造带的Te低值区(52 5k m),通常Te值较低表明岩石圈的力学强度较低,构造活动性相对较强。重力均衡研究结果表明,太行山构造带处于欠均衡状态3 2。当岩石圈底部受到软流圈侵蚀、拆沉或机械拉张等破坏机制作用时,力学强度较低的地区将优先被破坏或破坏程度较为显著。北部燕山构造带地壳内普遍存在高、低速互层现象,表现出大规模岩浆活动特征3 3;地震成像研究表明,燕山地区壳幔边界强烈的速度差异可能来自下地壳和岩石圈拆沉诱发的软流圈抬升,从而形成速度显著不同的物质分界面3 4。古太平洋的西向俯冲作用导致华北克拉通下方产生不稳定的地幔流动体系,并引起上地幔减压熔融或地幔物质向上流动3 5,强烈的对流、橄榄岩-熔体相互作用及伸展作用也可能导致岩石圈软化和克拉通破坏3 3。这与该地区低Te值(高活动性)和高F值(地表与地下负载共同作用)特征相吻合。华北克拉通东部为Te中、低值区(5 3 5k m)。郯庐断裂带起源于华北与华南板块的碰撞造山过程,郯庐断裂带较低的岩石圈强度及有利的岩浆活动通道,使其在克拉通破坏中成为伸展活动带和岩浆活动带,是华北克拉通破坏最严重的区域3 4,3 6。郯庐断裂带全线及其次级断层成为软流圈物质上涌的首选通道,这些热物质上升后,不断通过热作用和化学作用从底部侵蚀岩石圈,从而逐渐导致华北克拉通岩石圈的减薄和破坏。3.2 Te分布与地震活动性关系部分学者认为,岩石圈的强度能够控制岩石圈在构造力作用下的局部变形,因此与地震分布具有一定关联性3 7-3 8。从图9中Te值与地震空间分布来看,华北地区地震主要沿Te值较低(52 0k m)或高低值过渡带区域分布,如岩石圈强度较低的山西断裂带、太行山构造带和渤海湾盆地北部等地区。而在Te值较高区,如华北克拉通西部的鄂尔多斯地块内部,地震则较少发生。这表明在Te值较高区其岩石圈强度也较高,而Te值较低区岩石圈强度也较低,在构造应力作用下更易于释放积累的应变,从而更易诱发地震。由于岩石圈有效弹性厚度并不完全代表岩石圈的真实厚度,地震活动性与Te的关系也较为复杂,Te低值区并不一定能够代表地震频发区,地震活动性主要与其所处的构造环境有关。如太行山造山带南部和秦岭-大别造山带Te值虽然较小,但很少发生地震。秦岭-大别造山带形成于扬子板块和华北克拉通碰撞造山运动,区域内发育大规模的岩浆和构造变形活动3 9-4 0,表明该区域岩石圈强度的降低主要是由于深部热物质上升和流动过程中发生强烈的壳幔改造作用,发生地震所需的能量被释放,导致该区域地震活动较少。在高Te值的黄海地区,地震活动却很频繁,这是因为江苏陆区和南黄海海域的基底同属于华南下扬子地块4 1,嘉山-响水断裂和苏州-湖州断裂等大型断裂带从陆地延伸至海域,作为主要的边界断裂控制着该区域的地震活动特征。此外,地震分布与F值大小之间并未发现明显的对应关系,因此初始负载比F对孕震环境的影响作用尚不明确。3.3 Te与地表热流关系地球内部的热作用是地质构造运动和地震活动的主要动力来源之一。大地热流包含丰富的岩石圈热结构和地球动力学信息4 2。Te值大小还受岩石圈热结构影响,因此Te与地表热流分布之149大 地 测 量 与 地 球 动 力 学2 0 2 3年9月间存在特定联系6。图1 1为根据中国大陆地区地表热流数据库4 3采用克里金插值方法绘制的华北地区地表热流分布图,图中圆点代表地表热流观测点。可以看出,华北地区地表热流数据观测点较少,且分布不均匀。华北地区热流值整体较高,一般来说,高地表热流反映岩石圈构造较为活跃,代表岩石圈强度较弱,即Te与地表热流值呈现负相关关系。例如低地表热流值和高Te值的鄂尔多斯地区,构造活动较为稳定,代表稳定的克拉通区域;而高地表热流值和低Te值的山西断裂带与华北平原,构造活动相对活跃。但是在苏北盆地和渤海湾地区却表现出高热流值和高Te值,这可能是因为受到华北克拉通破坏的影响。在太平洋板块俯冲以及后撤作用下,中国东部构造体制由早期的挤压转换成伸展4 4,在该背景之下,深部地幔上隆、岩石圈伸展、大量热物质上涌带来巨大的热量,使整个地区产生活化作用,导致苏北盆地和渤海湾地区地表热流值较高。图1 1 华北地区地表热流分布F i g.1 1 D i s t r i b u t i o no f t e r r e s t r i a lh e a t f l o wi nn o r t hC h i n a4 结 语本文使用E I G E N-6 C 4模型布格重力异常、V 1 9.1地形和C RU S T 1.0模型地壳结构数据,采用导纳和相关函数联合反演方法,得到华北地区岩石圈有效弹性厚度Te和初始负载比F分布,讨论Te与华北克拉通破坏、地震活动和地表热流的关系,得到以下结论:1)华北克拉通Te值分布整体呈现横向不均匀性,这种明显的不均匀性证实了东部发生克拉通破坏和岩石圈减薄,而西部仍保持相对稳定的特征。2)从地震空间分布来看,地震活动与Te值关系复杂。华北地区地震活动主要沿Te相对低值区发生,表明在岩石圈强度薄弱区容易释放积累的应力与能量,更有利于地震的孕育和发生。而Te高值区构造则相对稳定,整体具有较强的抵抗岩石圈变形的能力。地震是构造应力长期积累和突然释放的结果,地震活动性主要受其所处的深部构造环境影响。3)Te值大小还受岩石圈热结构的影响,Te值与地表热流值整体对应关系较好,但在苏北盆地和渤海湾地区,Te值与地表热流值呈现负相关关系。其原因可能是由于热物质上升和流动,浅部地壳被破坏,产生高地表热流值,而深部地壳仍具有克拉通属性。参考文献1 F o r s y t hDW.S u b s u r f a c eL o a d i n g a n dE s t i m a t e s o f t h eF l e x u r-a lR i g i d i t yo fC o n t i n e n t a lL i t h o s p h e r eJ.J o u r n a l o fG e o p h y s i-c a lR e s e a r c h:S o l i dE a r t h,1 9 8 5,9 0(B 1 4):1 26 2 3-1 26 3 22 V e n i n g-M e i n e s zFA.C o mm e n t so nI s o s t a s yM.W a s h-i n g t o n:N a t i o n a lR e s e a r c hC o u n c i l,1 9 3 23 A i r yGB.O nt h eC o m p u t a t i o no f t h eE f f e c to f t h eA t t r a c-t i o no fM o u n t a i n-M a s s e s,a sD i s t u r b i n gt h eA p p a r e n tA s-t r o n o m i c a lL a t i t u d eo fS t a t i o n si n G e o d e t i cS u r v e y sJ.P h i l o s o p h i c a lT r a n s a c t i o n so f t h eR o y a l S o c i e t yo fL o n d o n,1 8 5 5,1 4 54 P r a t t VJH.O nt h eA t t r a c t i o no f t h eH i m a l a y aM o u n t a i n sa n do f t h eE l e v a t e dR e g i o n sb e y o n d t h e m,u p o n t h eP l u m b-L i n e i nI n d i aJ.M o n t h l yN o t i c e so f t h eR o y a lA s t r o n o m i-c a lS o c i e t y,1 8 5 5,1 6(2):3 6-4 15 B a r r e l l J.T h eS t r e n g t ho f t h eE a r t hsC r u s tJ.T h e J o u r-n a l o fG e o l o g y,1 9 1 4,2 2(7):6 5 5-6 8 36 陈波.中国及邻区岩石圈有效弹性厚度及其动力学意义D.武汉:中国地质大学,2 0 1 3(C h e nB o.T h eE f f e c t i v eE l a s t i cT h i c k n e s so v e r C h i n aa n d S u r r o u n d i n g sa n dI t sL i t h o s p h e r eD y n a m i c I m p l i c a t i o nD.W u h a n:C h i n aU n i-v e r s i t yo fG e o s c i e n c e s,2 0 1 3)7 D o r m a nLM,L e w i sBTR.E x p e r i m e n t a l I s o s t a s y:T h e o-r yo f t h eD e t e r m i n a t i o no f t h eE a r t hs I s o s t a t i cR e s p o n s e t oaC o n c e n t r a t e dL o a dJ.J o u r n a l o fG e o p h y s i c a lR e s e a r c h,1 9 7 0,7 5(1 7):33 5 7-33 6 58 金煜,姜效典.岩石圈动力学M.北京:科学出版社,2 0 0 2(J i nY u,J i a n gX i a o d i a n.L i t h o s p h e r eD y n a m i c sM.B e i j i n g:S c i e n c eP r e s s,2 0 0 2)9 J i nY,M c N u t tM K,Z h uYS.E v i d e n c e f r o mG r a v i t ya n dT o p o g r a p h yD a t af o rF o l d i n go fT i b e tJ.N a t u r e,1 9 9 4,3 7 1(64 9 9):6 6 9-6 7 41 0M c K e n z i eD,F a i r h e a dD.E s t i m a t e s o f t h eE f f e c t i v eE l a s t i cT h i c k n e s so f t h eC o n t i n e n t a lL i t h o s p h e r e f r o mB o u g u e r a n dF r e eA i rG r a v i t yA n o m a l i e sJ.J o u r n a l o fG e o p h y s i c a lR e-s e a r c h:S o l i dE a r t h,1 9 9 7,1 0 2(B 1 2):2 75 2 3-2 75 5 21 1P r e z-G u s s i n y M,L o w r yAR,W a t t sAB,e ta l.O nt h eR e c o v e r yo f E f f e c t i v e E l a s t i c T h i c k n e s s U s i n g S p e c t r a lM e t h o d s:E x a m p l e s f r o mS y n t h e t i cD a t aa n d f r o mt h eF e n-249 第4 3卷第9期孙文魁等:华北地区岩石圈有效弹性厚度及其地球动力学意义n o s c a n d i a nS h i e l dJ.J o u r n a lo fG e o p h y s i c a lR e s e a r c h:S o l i dE a r t h,2 0 0 4,1 0 9(B 1 0)1 2K i r b yJF,S w a i nCJ.G l o b a l a n dL o c a l I s o s t a t i cC o h e r e n c ef r o mt h e W a v e l e tT r a n s f o r mJ.G e o p h y s i c a lR e s e a r c hL e t t e r s,2 0 0 4,3 1(2 4)1 3佘雅文,付广裕,高原,等.华北地区中东部岩石圈挠曲与均衡特性以及地震活动性分析J.地球物理学报,2 0 1 8,6 1(1 1):44 4 8-44 5 8(S h eY a w e n,F uG u a n g y u,G a oY u a n,e ta l.F l e x u r eo ft h eL i t h o s p h e r e,I s o s t a t i cC h a r a c t e r i s t i ca n dS e i s m i c i t y i nM i d d l eE a s tA r e ao fN o r t hC h i n aJ.C h i-n e s eJ o u r n a l o fG e o p h y s i c s,2 0 1 8,6 1(1 1):44 4 8-44 5 8)1 4高山,章军锋,许文良,等.拆沉作用与华北克拉通破坏J.科学通报,2 0 0 9,5 4(1 4):19 6 2-19 7 3(G a oS h a n,Z h a n gJ u n f e n g,X uW e n l i a n g,e t a l.D i s m a n t l i n ga n dS u b-s i d e n c ea n dt h eF a i l u r eo fN o r t hC h i n aC r a t o nJ.C h i n e s eS c i e n c eB u l l e t i n,2 0 0 9,5 4(1 4):19 6 2-19 7 3)1 5唐新功,陈永顺,严良俊,等.应用布格重力异常研究太行山地区地壳密度结构J.西北地震学报,2 0 0 8,3 0(4):3 0 5-3 0 9(T a n gX i n g o n g,C h e nY o n g s h u n,Y a nL i a n g j u n,e t a l.R e s e a r c ho nC r u s t a lD e n s i t yS t r u c t u r ei nt h eP i e d-m o n tF a u l tZ o n eo fT a i h a n gM o u n t a i nA r e aU s i n g t h eB o u-g u e rG r a v i t yD a t aJ.N o r t h w e s t e r nS e i s m o l o g i c a l J o u r n a l,2 0 0 8,3 0(4):3 0 5-3 0 9)1 6嵇少丞