华北地块及邻区重力异常特征及地质结构

第４５卷㊀第５期２０２３年９月
地㊀震㊀工㊀程㊀学㊀报
C H I N A E A R T H Q U A K EE N G I N E E R I N GJ O U R N A L
V o l ．４５㊀N o ．５
S e p
t e m b e r ,２０２３㊀㊀收稿日期:２０２３Ｇ０３Ｇ１３
㊀㊀基金项目:河北省地震灾害防御与风险评价重点实验室开放基金资助项目(F Z ２２３２０５
)㊀㊀第一作者简介:郭雨帆(１９９２－),硕士,工程师,主要从事地震监测预报方面的工作.E Ｇm a i l :g u o y u f a n _o u c ＠１２６．c o m .㊀㊀通信作者:肖㊀健(１９７５－),高级工程师,主要从事地震监测预报方面的工作.E Ｇm a i l :５９４６７１９９０＠q q
．c o m .郭雨帆,肖健,付光明,等．华北地块及邻区重力异常特征及地质结构[J ]．地震工程学报,２０２３,４５(５):１０７５Ｇ１０８４．D O I :１０．
２００００/j
．１０００Ｇ０８４４．２０２３０３１３００５G U O Y u f a n ,X I A OJ i a n ,F U G u a n g m i n g ,e t a l ．G r a v i t y a n o m a l y c h a r a c t e r i s t i c s a n d g e o l o g
i c a l s t r u c t u r e o f t h eN o r t hC h i n a b l o c k a n d i t s s u r r o u n d i n g a r e a s [J ]．C h i n a E a r t h q u a k e E n g i n e e r i n g J o u r n a l ,２０２３,４５(５):１０７５Ｇ１０８４．D O I :１０．２００００/j ．１０００Ｇ０８４４．２０２３０３１３００５
华北地块及邻区重力异常特征及地质结构
郭雨帆１,
２
,肖㊀健１,付光明１,许志山１,李清武１,黄仁桂１,赵㊀影１,徐长文１,肖建华１,袁㊀俊１
(１．江西省地震局,江西南昌３３００９６;
２．河北省地震灾害防御与风险评价重点实验室,河北三河０６５２０１
)摘要:研究华北地块及邻区重力异常特征㊁
地质结构和断裂构造,对该区地震风险研判和矿产资源远景规划具有重要意义.文章基于小波多尺度分解法㊁P a r k e r ＧO l d e n b u r g 迭代反演法和重力异常导数法,利用华北地块及邻区布格重力异常数据,获得布格重力异常小波分解逼近场和细节场㊁莫霍面深度和重力异常导数.结果表明:(１)逼近场和细节场揭示出华北地块可分为西部的鄂尔多斯克拉通块体和东部的华北似环状裂谷盆地两个子单元.(２)莫霍面深度范围在２８~５２k m 之间,
呈N N E 走向,由西向东逐渐变浅,具有显著的 东西分带 特征.尤为显著的是鄂尔多斯克拉通块体和华北似环状裂谷盆地,块体内部莫霍面变化较平缓,块体边界莫霍面变化剧烈.(３)根据重力异常导数结果识别并校正出研究区１３条块体边界及主要断裂带.
关键词:华北地块;重力异常;莫霍面;断裂中图分类号:P ３１５㊀㊀㊀㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:１０００Ｇ０８４４(２０２３)０５－１０７５－１０
D O I :１０．２００００/j
．１０００Ｇ０８４４．２０２３０３１３００５G r a v i t y a n o m a l y c h a r a c t e r i s t i c s a n d g e o l o g
i c a l s t r u c t u r e o f t h eN o r t hC h i n ab l o c ka n d i t s s u r r o u n d i n g a
r e a s G U O Y u f a n １,
２,X I A OJ i a n １,F U G u a n g m i n g １,X UZ h i s h a n １,L IQ i n g w u １,HU A N G R e n g
u i １
,Z H A O Y i n g １,X U C h a n g
w e n １,X I A OJ i a n h u a １,Y U A NJ u n １
(１．J i a n g x iE a r t h q u a k eA g e n c y ,N a n c h a n g ３３００９６,J i a n g
x i ,C h i n a ;２．H e b e iK e y L a b o r a t o r y o f E a r t h q
u a k eD i s a s t e rP r e v e n t i o na n dR i s kA s s e s s m e n t ,S a n h e ０６５２０１,H e b e i ,C h i n a )A b s t r a c t :E x a m i n i n g t h e g r a v i t y a n o m a l y c h a r a c t e r i s t i c s ,g e o l o g
i c a l s t r u c t u r e ,a n d f a u l t s t r u c t u r e o f t h eN o r t hC h i n ab l o c ka n ds u r r o u n d i n g a r e a s i ss i g n i f i c a n t f o ra s s e s s i n g e a r t h q
u a k er i s ka n d l o n g Ｇt e r m m i n e r a l r e s o u r c e f o r m a t i o n i n t h i s r e g i o n ．B a s e do n t h eB o u g u e r g r a v i t y a n o m a l y d
a t a f r o mt h eN o r t hC h i n a
b l o
c k a n
d s u r r o u n d i n g a r
e a s ,t h ew a v e l e t a p p
r o x i m a t i o n s a n d d e t a i l s o f t h e B o u g u e r g r a v i t y a n o m a l y ,M o h od e p t h ,a n d g r a v i t y a n o m a l y d
e r i v a t i v e s a r eo b t a i n e dv i am u l t i Ｇ
s c a l ew a v e l e t d e c o m p o s i t i o n,P a r k e rＧO l d e n b u r g i t e r a t i v e i n v e r s i o n,a n d g r a v i t y a n o m a l y d e r i v aＧt i v em e t h o d s．T h ek e y r e s u l t sa r ea s f o l l o w s:(１)t h ew a v e l e t a p p r o x i m a t i o n sa n dd e t a i l ss h o w t h a t t h eN o r t hC h i n ab l o c k c a nb e d i v i d e d i n t o t w o s u b u n i t s,n a m e l y,t h eO r d o sC r a t o nb l o c k i n t h ew e s t a n d t h eN o r t hC h i n a r i n gＧl i k e r i f t b a s i n i n t h e e a s t．(２)T h eM o h od e p t h i s i n t h e r a n g e o f２８－５２k m;i t s t r i k e sN N Ea n d g r a d u a l l y d e c r e a s e s f r o m w e s t t o e a s t,w i t h o b v i o u s e a s tＧw e s t z o n i n g． I n t h eO r d o sC r a t o nb l o c ka n d t h eN o r t hC h i n a r i n gＧl i k e r i f t b a s i n,t h eM o h od e p t h i nＧs i d e t h eb l o c k r e m a i n s s t a b l e,w h e r e a s t h a tb e l o wt h eb l o c kb o u n d a r i e s c h a n g e sd r a s t i c a l l y．(３) T h i r t e e nb l o c kb o u n d a r i e sa n d m a j o r f a u l tz o n e s i nt h es t u d y a r e aa r e i d e n t i f i e da n dc o r r e c t e d b a s e do n t h e d e r i v a t i v e r e s u l t s o f g r a v i t y a n o m a l y．
K e y w o r d s:N o r t hC h i n ab l o c k;g r a v i t y a n o m a l y;M o h o;f a u l t
０㊀引言
中国大陆及邻区可划分为东北地块㊁华北地块㊁华南地块㊁西域地块㊁青藏地块和滇缅地块[１Ｇ２],其中,华北地块位于六大地块中心位置,区域内岩浆和断裂带活跃㊁强震频发㊁矿产资源丰富[３].华北地块由西部的鄂尔多斯克拉通块体和东部的华北似环状裂谷盆地组成[４](图１).华北克拉通于古元古代(２５００~１６００M a)逐渐形成.中元古代 三叠纪(１６００~２０５M a)是华北地台发育期,形成了典型的地台型沉积.在三叠纪末期(约２０５M a),华北地块
虚线为块体边界及主要断裂:１．华北地块边界;２．鄂尔多斯克拉通块体与华北似环状裂谷盆地分界;３．河套断陷盆地断裂系;４．包头 张家口 平泉深断裂;５．济源 黄口深断裂;６．山西断陷盆地断裂系;７．渭河断陷盆地断裂系;８．太行山断裂带;９．邢台 唐山新断裂带;１０．沧东断裂带;１１．麻城 崇阳 新宁断裂;１２．郯城 庐江断裂带;１３．连云港 千里岩 开城断裂;红线为莫霍面测线,白圆点为莫霍面控制点(列于表１)
图１㊀研究区水深地形和断裂展布图
F i g．１㊀R e g i o n a l t o p o g r a p h y a n d f r a c t u r e d i s t r i b u t i o nm a p 和扬子地块碰撞形成秦岭 大别山构造带和郯庐断裂带,华北地台与蒙古微大陆通过兴蒙 吉黑造山带焊接.侏罗纪 白垩纪(２０５~６５M a),燕山运动造成华北地台分化,表现为华北地台东部大陆 活化 ,强烈的岩浆活动和构造运动导致地壳破坏改造,而西部鄂尔多斯地区保持稳定,基本不受影响.新生代以来(６５M a至今),华北地区进入大陆伸展构造发育的环境,东部裂谷盆地岩石圈拆沉㊁地壳减薄,而西部鄂尔多斯块体较稳定,并整体隆升[４].
目前华北地块及邻区重力研究成果相对较少. L i等[５]使用E GM２００８重力模型反演了华北克拉通岩石圈密度,发现克拉通活化和减薄尚未影响到华北克拉通中㊁西两个部分.玄松柏等[６]基于E GM２００８重力模型数据,利用P a r k e rＧO l d e n b u r g 方法反演了中国大陆莫霍面深度.石岚[７]采用N G D C７２０地磁模型和WGM２０１２重力模型分别反演了华北克拉通及邻区居里等温面和莫霍面,并探讨了这两个界面的构造意义.许田等[８]基于E GM２００８重力模型及E T O P O１地形模型,研究了华北地区重力和地形及其对克拉通破坏深部过程的约束作用.侍文等[９]使用WGM２０１２数据计算中国大陆岩石圈有效弹性厚度和初始载荷比,以探究强震构造区力学特征.上述研究没有针对华北地块及邻区地质结构进行系统而全面地分析,缺乏使用重力资料反演该区莫霍面并识别校正断裂的研究成果.基于上述原因,有必要使用重力数据,对华北地块及邻区地质结构与构造特征开展针对性研究.本文采用小波多尺度分解法对布格重力异常数据进行处理,得到华北地块及邻区布格重力异常小波分解逼近场和细节场.在此基础上,使用P a r k e rＧO l d e nＧb u r g法反演莫霍面深度,并结合断裂资料,利用重力异常导数法识别和校正块体边界和主要断裂带.
６７０１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
目的是为深入分析华北地块及邻区重力异常特征㊁
地质结构㊁断裂构造提供基础地球物理证据和相关
研究成果,为该区地震风险研判及矿产资源远景规划提供技术支撑和参考依据.
１㊀数据和方法
１．１㊀数据
布格重力异常数据来源于全球重力场数据库(B u r e a uG r a v i m ét r i q
u e I n t e r n a t i o n a l ,简称B G I )使用E GM ２００８模型计算得到的重力异常数据[１
０]
.该数据空间分辨率约为２．５ᶄˑ２．５ᶄ,地形改正使用了E T O P O １模型,布格改正应用了F u l l e a 等[１
１]
开发的F A ２B O U G H T 代码.E GM ２００８模型是由美国
国家地理空间情报局(U．S ．N a t i o n a lG e o s p a t i a l I n t e l l i g e n c eA g e n c y ,简称N G A )E GM 开发团队发布的,参考椭球面为W G S ８４椭球,它采用了G R A C E
卫星跟踪数据㊁卫星测高数据和地面重力数据.
本文主要研究区域为华北地块及邻区(１０２ʎ~
１２２ʎE ,２９ʎ~４５ʎN ),研究区布格重力异常如图２.重力异常总体呈N N E 向条带状分布,从西部到东部异常值逐渐增大.在西南部形成的重力异常极值圈闭对应青藏地块,它是整个研究区的重力异常低值区,最低值可达－５４３．４９m g a l .华北地块西部鄂尔多斯克拉通块体的重力异常值低于东部华北似环状裂谷盆地的布格重力异常,
表明这两个地区地壳
黑虚线为块体边界及主要断裂
图２㊀华北地块及邻区布格重力异常场
F i g ．２㊀B o u g u e r g r a v i t y a n o m a l y f
i e l d i n t h eN o r t hC h i n a b l o c ka n d s u r r o u n d i n g a
r e a s 结构和介质密度存在差异.从图２中还可以识别出秦岭 大别山构造带,其重力异常值呈NWW 向
延伸.
１．２㊀方法
小波多尺度分解法凭借其明显的技术优势已被应用到众多研究领域,在地球科学领域如重力异常场分解研究中也得到广泛应用,它可以将重力异常
场分解成不同尺度成分,实现重力异常场分离[
１２Ｇ１３]
.本次研究采用杨文采等[１２]和候遵泽等[１４]的二维多尺度小波分解法对布格重力异常数据进行处理,设重力异常为Δg (x ,y )
,其分解公式为:Δg (x ,y )
＝A i ＋D i ＋D i －１＋ ＋D １㊀(１)式中:A i 为重力异常的i 阶(i 为不小于２的整数)近似,即重力异常的低频成分;D i 为经i 次分解后得到的各阶小波细节,即重力异常的高频成分.在重力异常的小波多尺度分解中,确定最优小波基是一个关键的问题,使用不同的小波基函数会产生不
同的分解结果[１５].牟力等[１６
]经过理论层面分析㊁
理论模型对比试验和华北地区实测重力数据对比试验,认为d b １１小波基为重力数据多尺度分析的最优小波基,因此本次研究由d b １１小波构造二维尺度函
数和小波函数.刁博等[１７]
通过讨论信号长度㊁小波
母函数的支撑长度与分解阶次的关系,发现对重力异常进行小波分解时,恰当的小波分解阶次为５阶,故本次研究将小波分解尺度确定为５阶.
P a r k e r 给出频率域内计算连续密度界面在地
表平面内产生的重力异常公式[１８]
后,O l d e n b u r g 根据该公式的傅里叶级数展开,提取出第一项,得到密
度界面反演的迭代公式[
１９]
:F [h (x )]＝－F [Δg (x )]e (－k
z ０)２πG ρ
－
ðɕn ＝２(k )
n －１
n
!F [h n (x )](２)式中:F [h (x )]㊁F [Δg (
x )]分别是界面深度和重力异常的傅里叶变换;z ０为参考深度;k 为波数;G 为
万有引力常数;ρ为壳Ｇ幔密度差.本次研究中,设置截断频率WH ＝０．０１,S H ＝０．０１２.据表１的莫霍
面控制点[２０
],它们的平均深度为３８．８k m ,
结合反演的实际情况,设置z ０＝３８k m .
重力异常导数法是重力勘探中识别和分析断裂
构造最常用的方法.０ʎ水平方向导数㊁９０ʎ
水平方向导数和垂向二次导数(罗森巴赫公式)对识别不同产状断裂效果较好且各有优势.本文对剩余布格重力
异常小波分解３阶细节场[２１]
求０ʎ水平方向导数㊁
７７０１第４５卷第５期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀郭雨帆,
等:华北地块及邻区重力异常特征及地质结构㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀
９０ʎ
水平方向导数和垂向二次导数,并根据导数图的线性特征㊁极值带等标识特征,识别和校正基底及以上断裂[
２２Ｇ２３]
.２㊀结果与讨论
２．１㊀布格重力异常场小波多尺度分解结果
以d b １１为小波母函数进行小波多尺度分解,得到研究区布格重力异常场小波１~５阶逼近场(图
３
)和细节场(图４).各阶小波逼近场图像中反映的重力场源深度不同,但均可看出研究区不同大地构造单元的基本格局,即西南部的低值圈闭对应青藏地块,南部的华南地块与北部的华北地块通过秦岭 大别山构造带隔开,华北地块内部可以识别出重力异常较低的鄂尔多斯克拉通块体和重力异常较高的华北似环状裂谷盆地.低阶逼近场与原始布格重力异常形态接近,随着小波逼近阶数增大,不同构造单元及深大断裂的重力异常的浅源㊁局部信息得到了较好的抑制,区域异常信息显得更为突出.前人研究表明４阶逼近
场可能是莫霍面引起的异常[
１２,２１,２４Ｇ２７
].通过对比不同阶次的逼近场(图３),发现４阶逼近场既能一定程度地消除浅部信息影响,又对深部信息有较好反映,初步考虑使用４阶逼近场反演莫霍面深度.将
小波分析法和径向对数功率谱法[２８]
两种方法结合有利于对地质结构进行定量分析[１２,２４,２９Ｇ３１]
,本次研
究使用径向对数功率谱法估计４阶逼近异常场源埋深,计算结果为３７．３k m ,
与表１的莫霍面控制点[２０]
平均深度３８．８k m 接近.可以推断４阶逼近场可能
反映了研究区莫霍面的起伏,故选取４阶逼近场反演莫霍面深度
.
黑虚线为块体边界及主要断裂
图３㊀华北地块及邻区布格重力异常小波１~５阶逼近场
F i g ．３㊀T h e １s t t o ５t hw a v e l e t a p p r o x i m a t i o n s o fB o u g u e r g r a v i t y a n o m a l y f o r t h eN o r t hC h i n ab l o c ka n d s u r r o u n d i n g a
r e a s ㊀㊀１阶和２阶小波细节[图４(a )㊁(b
)]场源信息丰富而复杂,主要反映的是浅部高频重力异常信息.其分布特征显示出西南部青藏地块的地壳横向变化剧烈且比较破碎,而华北地块的重力异常变化则相对平缓,大部分数据在－５~５m g
a l 之间变化,尤其是华北似环状裂谷盆地的重力异常变化更为平缓.研究区重力异常的走向㊁规模㊁分布位置与区内断裂一致性较好,例如太行山断裂带东西两侧的重力异常差异明显,能清晰地勾勒出太行山断裂带轮廓.随着小波细节阶数的增加,重力异常变得宽缓平滑,
８７０１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
主要体现了深部场源所引起的较大规模的低频重力异常信息.同３阶小波细节[图４(c )]相比,４阶和
５阶小波细节[图４(d )㊁(４e
)]中小规模的条带状及串珠状异常逐渐消失,重力异常圈闭现象更加明显,
圈闭空间逐渐扩大,呈带状和块状分布.４阶和５阶小波细节显示块体边界处重力异常变化剧烈,表明块体边界深大断裂对重力异常的分布特征具有较
好的控制作用
.
黑虚线为块体边界及主要断裂
图４㊀华北地块及邻区布格重力异常小波１~５阶细节场
F i g ．４㊀T h e １s t t o ５t hw a v e l e t d e t a i l s o fB o u g u e r g r a v i t y a n o m a l y f o r t h eN o r t hC h i n ab l o c ka n d s u r r o u n d i n g a
r e a s ㊀㊀综上所述,逼近场和细节场共同显示华北地块内部出现西㊁东两个子单元,分别对应鄂尔多斯克拉通块体和华北似环状裂谷盆地.子单元内部变化较为平缓,而在子单元边界及块体边界处变化较剧烈.
２．２㊀莫霍面特征
莫霍面是地壳与地幔的分界面,本文通过P a r Ｇk e r ＧO l d e n b u r g 方法迭代反演得到莫霍面深度.图５(a )是华北地块及邻区莫霍面埋深平面图,将其与C R U S T １．０地壳厚度模型图[图５(b )]对比,发现两者变化趋势一致,但本次研究得到的莫霍面深度变化细节更加丰富.同时,与前人深地震探测得到的
莫霍面控制点[２０]
进行比较(表１),发现莫霍面深度
偏差较小,说明本次研究得到的莫霍面深度结果可信度较高.
由图５(a )可知,研究区莫霍面埋深范围在２８~
５２k m 之间,横向上起伏较大,等值线呈N N E 走
向,从西到东莫霍面逐渐变浅,呈现 东西分带 的特点.第一条带主要集中于研究区西南的青藏地块,该地区莫霍面深度可达５２k m ;第二条带包括四川盆地和鄂尔多斯克拉通块体,后者莫霍面深度略大于前者,约为４０k m ;
第三条带主要为华北似环状裂谷盆地(莫霍面平均深度为３４k m )
,及华南地块和东北地块的部分区域.石岚[７]使用WGM ２０１２重力模型反演了华北克拉通及邻区莫霍面深度,分析
认为华北克拉通西部 中部和东部莫霍面埋深分别为３７~４５k m 和３３~３６k m ,与本次研究 鄂尔多斯克拉通块体莫霍面深度约为４０k m ㊁华北似环状裂谷盆地莫霍面深度约为３４k m
的结果基本一致.值得注意的是,郯庐断裂带走向与其莫霍面等值线走向近乎平行,且莫霍面深度较同纬度其他区域更浅,可能指示了郯庐断裂带为地幔热物质上涌提供
９７０１第４５卷第５期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀郭雨帆,
等:华北地块及邻区重力异常特征及地质结构㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀
了通道[３２].除此之外,－３５k m莫霍面埋深等值线在(１１６ʎE,３０ʎN)附近由北部的N E向转为西部的近E W向,支持了乔计花等[３３]关于郯庐断裂带并未跨过长江继续向南延伸的观点.
黑线为莫霍面测线,黑圆点为莫霍面控制点(列于表１),黑虚线为块体边界及主要断裂
图５㊀华北地块及邻区莫霍面埋深平面图及C R U S T１．０地壳厚度模型图
F i g．５㊀M o h od e p t h p l a na n dC R U S T１．０c r u s t a l t h i c k n e s sm o d e l f o r t h e
N o r t hC h i n ab l o c ka n d s u r r o u n d i n g a r e a s
表１㊀莫霍面深度反演结果与数据控制点结果[２０]对比
T a b l e１㊀C o m p a r i s o nb e t w e e n i n v e r s i o n r e s u l t s a n d c o n t r o l
p o i n t r e s u l t s o fM o h od e p t h[２０]
编号
经度
/(ʎE)纬度
/(ʎN)控制点值
/k m反演结果
/k m偏差/k m
(１)１０６．８２３７．２５４２．００４２．２３－０．２３(２)１０７．５９３７．４７４１．００４２．５４－１．５４(３)１０８．８３３９．６８４２．００４０．６５１．３５(４)１１２．７３３５．９４４０．００３８．９２１．０８(５)１１３．２８３６．２８３８．００３９．０１－１．０１(６)１１３．４７４１．４８４１．００４０．４６０．５４(７)１１５．７６３６．１４３５．００３５．４５－０．４５(８)１１７．１３３６．３８３５．００３５．８３－０．８３(９)１２０．４６３７．０９３５．００３５．０００
结合研究区水深地形及构造单元特征,选取A Aᶄ㊁B Bᶄ㊁C Cᶄ三条测线绘制莫霍面埋深剖面图(图６).在A Aᶄ测线剖面图中,莫霍面埋深与水深地形呈现良好的 镜像对称 关系.从西北到东南,莫霍面逐渐变浅,形成四级上升台阶,分别对应①阿拉善地块㊁②鄂尔多斯克拉通块体㊁③华北似环状裂谷盆地陆区和④黄海海域.由B Bᶄ和C Cᶄ测线剖面图可以发现,莫霍面埋深在块体内部变化较平缓,而在块体边界呈 V 字形剧烈变化.C Cᶄ测线剖面的莫霍面比B Bᶄ更浅,可能与华北地块东部岩石圈拆沉㊁地壳减薄作用较西部更显著有关[５].
２．３㊀断裂体系的重力异常特征
收集和整理研究区断裂资料[３Ｇ５,３４Ｇ４３],根据断裂构造在重力异常导数图上识别出六种标志:线性重力高与重力低之间的过渡带㊁异常轴线明显错动的部位㊁串珠状异常的两侧或轴部所在位置㊁两侧异常特征明显不同的分界线㊁封闭异常等值线突变的部位㊁等值线同行扭曲部位[４４],在０ʎ水平导数图㊁９０ʎ水平导数图和垂向二次导数图上识别和校正了华北地块１３条边界及内部主要断裂带(图７).邓起东等[４５Ｇ４６]依据主要活动带上的几何学和运动学定量数据资料,编制了«１ʒ４００万中国活动构造图»,将其与这１３条边界及内部主要断裂带的平面分布位置进行对比分析,发现两者基本吻合,表明本次结果可靠性较高.由于篇幅所限,本文不再展示两者对比图件.下面描述这１３条边界及内部主要断裂带的重力异常场导数图特征.
一级块体边界:(１)华北地块边界.其西北边界为银川 吉兰泰断陷盆地断裂系,产状为４０ʎ/S E ø６０ʎ,正右旋性质,它在０ʎ水平导数图上沿着串珠状正负异常间分界线发育;北边界由西 中部的阴山北缘断裂带和东部的赤峰 开源深断裂组成,近
０８０１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
E W 走向,它们在９０ʎ
水平导数图上沿着串珠状正负异常间分界线分布,局部可见异常轴线错动;西南边界包括海原断裂和六盘山断裂,海原断裂NWW 走向,逆左旋性质,它在９０ʎ水平导数图上是正负高值与正负低值异常的分界线,识别度较高,而六盘山断
裂NW 走向,它在垂向二次导数图上沿串珠状正负
异常间分界线展布;南边界为秦岭 大别山构造带,NWW 走向,华南地块朝华北地块N N E 向逆冲,
它在９０ʎ水平导数图上大体沿着串珠状正负异常间分界线发育,局部存在异常轴线错动
.
①㊁②㊁③㊁④分别对应阿拉善地块㊁
鄂尔多斯克拉通块体㊁华北似环状裂谷盆地陆区㊁黄海海域图６㊀华北地块及邻区莫霍面埋深剖面图
F i g ．６㊀M o h od e p t h p r o f i l e s o f t h eN o r t hC h i n ab l o c ka n d s u r r o u n d i n g a
r e a
s 蓝色㊁绿色㊁黑色虚线分别代表一㊁二块体边界及三级断裂,名称及编号同图１
图７㊀华北地块及邻区重力异常导数及其断裂解释图
F i g ．７㊀
G r a v i t y a n o m a l y d e r i v a t i v e a n d i t s c o r r e l a t i o n t o f r a c t u r e s i n t h eN o r t hC h i n ab l o c ka n d s u r r o u n d i n g a
r e a s ㊀㊀二级块体边界:(２
)鄂尔多斯克拉通块体与华北似环状裂谷盆地分界,位于山西断陷盆地断裂系与太行山断裂带之间,N N E 走向,它在垂向二次导数图上沿异常轴线明显错动的部位和串珠状正负异常间分界线发育.
三级断裂:(３
)河套断陷盆地断裂系,近E W 走向,N 或S 倾向,陡倾角,正断性质,它在９０ʎ水平导数图上沿串珠状正负异常间分界线发育;(４)包头 张家口 平泉深断裂,近E W 向展布,正左旋性质,它在９０ʎ
水平导数图上沿着串珠状正负异常间分界线分布,局部可见异常轴线错动;(５)济源 黄口深断裂,NWW E W 走向,它在９０ʎ水平导数图上沿１８０１第４５卷第５期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀郭雨帆,
等:华北地块及邻区重力异常特征及地质结构㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀
着串珠状正负异常间分界线发育;(６)山西断陷盆地断裂系,产状为N N E N E/NW或S Eø(４０ʎ~８０ʎ),正右旋性质,它在０ʎ水平导数图上沿着串珠状正负异常间分界线和等值线同形扭曲部位分布;(７)渭河断陷盆地断裂系,产状为E W N E E/N或S ø(５０ʎ~８０ʎ),正左旋性质,它在９０ʎ水平导数图上沿着串珠状异常发育;(８)太行山断裂带,N N E N E走向,NWW NW或S E E S E倾向,正断性质,它在垂向二次导数图上沿异常轴部错动的部位展布;(９)邢台 唐山新断裂带,产状为５０ʎ/NW ø(７０ʎ~８５ʎ),正右旋性质,隐伏断裂,它在垂向二次导数图上沿异常轴线错动的部位和串珠状正负异常间分界线分布;(１０)沧东断裂带,走向约为３０ʎ,倾向为S E E,是新生代渤海湾盆地黄骅坳陷和沧县隆起共有的隐伏边界非活动断裂,它在垂向二次导数图上沿异常轴线错动的部位发育;(１１)麻城 崇阳 新宁断裂N E N N E走向,为鄂南㊁湘北湖盆区的东部边界断裂,它在垂向二次导数图上沿异常轴线错动的部位展布;(１２)郯城 庐江断裂带,包括北部的沂沐断裂和南部的宿迁 广济断裂,沂沐断裂产状为(５ʎ~２５ʎ)/S E或NWø(６０ʎ~８０ʎ),逆右旋性质,宿迁 广济断裂为N N E N E走向,逆右旋性质,它们在垂向二次导数图上沿串珠状正负异常间分界线和异常轴线错动的部位延伸;(１３)连云港 千里岩 开城断裂,产状为５０ʎ/S,正断性质,它在垂向二次导数图上沿串珠状正负异常间分界线分布.
３㊀结语
本文基于小波多尺度分解法㊁P a r k e rＧO l d e nＧb u r g迭代反演法和重力异常导数法对华北地块及邻区布格重力异常数据进行分析处理,获得了布格重力异常小波分解逼近场和细节场㊁莫霍面深度和重力异常导数,主要得出如下结论:
(１)布格重力异常小波分解逼近场图和细节场图揭示了华北地块可以分为西部的鄂尔多斯克拉通块体和东部的华北似环状裂谷盆地两个子单元,子单元内部变化较为平缓,而块体边界处变化较剧烈.
(２)莫霍面深度范围在２８~５２k m之间,呈现N N E走向,从西到东莫霍面逐渐变浅,具有显著的 东西分带 特征.尤为显著的是鄂尔多斯克拉通块体和华北似环状裂谷盆地,两者莫霍面平均深度分别为４０k m和３４k m,块体内部莫霍面变化较平缓,块体边界莫霍面呈 V 字形变化剧烈.
(３)根据重力异常导数结果识别并校正了研究区１３条块体边界及主要断裂带,包括１条一级块体边界,１条二级块体边界,１１条三级断裂带.
参考文献(R e f e r e n c e s)
[１]㊀张培震,邓起东,张国民,等．中国大陆的强震活动与活动地块[J]．中国科学(D辑:地球科学),２００３,３３(增刊１):１２Ｇ２０．
Z H A N G P e i z h e n,D E N G Q i d o n g,Z HA N G G u o m i n,e t a l．S t r o n g e a r t h q u a k e a c t i v i t y a n da c t i v eb l o c k s i nC h i n e s e M a i nＧl a n d[J]．S c i e n c e i nC h i n a(S e rD),２００３,３３(S u p p l０１):１２Ｇ２０．[２]㊀张国民,马宏生,王辉,等．中国大陆活动地块与强震活动关系[J]．中国科学(D辑:地球科学),２００４,３４(７):５９１Ｇ５９９．
Z H A N G G u o m i n,MA H o n g s h e n g,WA N G H u i,e ta l．R e l aＧt i o n s h i p b e t w e e n a c t i v e b l o c k a n d s t r o n g e a r t h q u a k e a c t i v i t y i n
C h i n e s eM a i n l a n d[J]．S c i e n c e i nC h i n a(S e rD),２００４,３４(７):
５９１Ｇ５９９．
[３]㊀李清武．华北中部地区N E和NW向断裂体系对现今构造应力场及地震孕育和触发控制的数值模拟研究[D]．青岛:中国海洋
大学,２０１５．
L IQ i n g w u．N u m e r i c a l s i m u l a t i o n s t u d y o n t h e p r e s e n t t e c t o n i c s t r e s s f i e l da n de a r t h q u a k e p r e p a r a t i o na n dt r i g g e r c o n t r o l b y N E a n d NW f a u l t s y s t e m si n c e n t r a l N o r t h C h i n a[D]．Q i n g d a o:O c e a nU n i v e r s i t y o fC h i n a,２０１５．
[４]㊀邓晋福,魏文博,邱瑞照．中国华北地区岩石圈三维结构及演化[M]．北京:地质出版社,２００７．
D E N G J i n f u,W E I W e n b o,Q I U R u i z h a o．T h r e eＧd i m e n s i o n a l
s t r u c t u r e o f l i t h o s p h e r e a n d i t s e v o l u t i o n i nN o r t hC h i n a[M]．
B e i j i n g:G e o l o g i c a l P u b l i s h i n g H o u s e,２００７．
[５]㊀L IY Y,Y A N G YS．G r a v i t y d a t a i n v e r s i o n f o r t h e l i t h o s p h e r i c
d e n s i t y s t r u c t u r e b e n e a t hN o r t hC h i n aC r a t o n f r o mE GM２００８
m o d e l[J]．P h y s i c so f t h eE a r t ha n dP l a n e t a r y I n t e r i o r s,２０１１,１８９(１/２):９Ｇ２６．
[６]㊀玄松柏,申重阳,谈洪波,等．利用E GM２００８重力数据反演中国大陆M o h o深度[J]．大地测量与地球动力学,２０１５,３５(２):３０９Ｇ３１１,３１７．
X U A NS o n g b a i,S H E NC h o n g y a n g,T A N H o n g b o,e t a l．I n v e rＧs i o no fM o h od e p t h i nC h i n aM a i n l a n d f r o m E GM２００８g r a v i t y
d a t a[J]．J o u r n a lo fG
e o d e s y a n d G e o d y n a m i c s,２０１５,３５(２):
３０９Ｇ３１１,３１７．
[７]㊀石岚．华北克拉通及邻区重磁异常特征及其构造意义[D]．昆明:云南大学,２０１９．
S H IL a n．C h a r a c t e r i s t i c so f g r a v i t y a n d m a g n e t i ca n o m a l i e s i n N o r t hC h i n aC r a t o na n di t sa d j a c e n ta r e a sa n dt h e i rt e c t o n i c s i g n i f i c a n c e[D]．K u n m i n g:Y u n n a nU n i v e r s i t y,２０１９．[８]㊀许田,黄金水．华北地区重力和地形及其对克拉通破坏深部过程的约束[J]．大地测量与地球动力学,２０２０,４０(５):５１７Ｇ５２１．
X U T i a n,HU A N GJ i n s h u i．C o n s t r a i n t s o f g r a v i t y a n d t o p o g r aＧp h y o n t h ed e s t r u c t i o no fN o r t hC h i n aC r a t o n[J]．J o u r n a lo f
G e o d e s y a n dG e o d y n a m i c s,２０２０,４０(５):５１７Ｇ５２１．[９]㊀侍文,陈石,韩建成,等．中国大陆岩石圈有效弹性厚度与强震构造区力学特征研究[J]．地震,２０２１,４１(１):１Ｇ１２．
２８０１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年
S H IW e n,C H E NS h i,H A NJ i a n c h e n g,e t a l．S t u d y o nt h ee fＧ
f e c t i v e e l a s t i c t h i c k n e s so f l i t h o s p h e r ea n d m e c h a n i c a l c h a r a cＧ
t e r i s t i c s o f g r e a t e a r t h q u a k e t e c t o n i c a r e a s i n t h eC h i n e s eM aＧ
i n l a n d[J]．E a r t h q u a k e,２０２１,４１(１):１Ｇ１２．
[１０]㊀P A V L I SN K,H O L M E SSA,K E N Y O NSC,e t a l．A ne a r t h
g r a v i t a t i o n a lm o d e l t od e g r e e２１６０:E GM０８[C]//t h e２００８
G e n e r a lA s s e m b l y o f t h eE u r o p e a n G e o s c i e n c e U n i o n．V i e nＧ
n a,A u s t r i a,２００８:１３Ｇ１８．
[１１]㊀F U L L E A J,F E R NÀN D E Z M,Z E Y E N H．F A２B O U G a
F O R T R A N９０c o d et oc o m p u t eB o u g u e r g r a v i t y a n o m a l i e s
f r o m
g r i d d e d f r e eＧa i r a n o m a l i e s:a p p l i c a t i o n t o t
h eA t l a n t
i c－
M e d i t e r r a n e a n t r a n s i t i o n z o n e[J]．C o m p u t e r s&G e o s c i e n c e s,
２００８,３４(１２):１６６５Ｇ１６８１．
[１２]㊀杨文采,施志群,侯遵泽,等．离散小波变换与重力异常多重分解[J]．地球物理学报,２００１,４４(４):５３４Ｇ５４１,５８２．
Y A N G W e n c a i,S H I Z h i q u n,HO U Z u n z e,e t a l．D i s c r e t e
w a v e l e tt r a n s f o r m f o r m u l t i p l e d e c o m p o s i t i o n o f g r a v i t y
a n o m a l i e s[J]．C h i n e s eJ o u r n a lo fG e o p h y s i c s,２００１,４４(４):
５３４Ｇ５４１,５８２．
[１３]㊀方盛明,张先康,嘉世旭,等．华北地区布格重力异常的多尺度分解特征与地震活动性[J]．大地测量与地球动力学,２００２,２２
(１):３４Ｇ３９．７．
F A N
G S h e n g m i n g,Z
H A N G X i a n k a n g,J
I A S h i x u,e t a l．
M u l t iＧs c a l ed e c o m p o s i t i o no fB o u g u e r g r a v i t y a n o m a l y a n d
s e i s m i c a c t i v i t y i n N o r t hC h i n a[J]．C r u s t a lD e f o r m a t i o na n d
E a r t h q u a k e,２００２,２２(１):３４Ｇ３９．
[１４]㊀侯遵泽,杨文采．小波多尺度分析应用[M]．北京:科学出版社,２０１２．
H O U Z u n z e,Y A N G W e n c a i．A p p l i c a t i o n o f w a v e l e t m u l t iＧ
s c a l e a n a l y s i s[M]．B e i j i n g:S c i e n c eP r e s s,２０１２．
[１５]㊀X U C,L I UZ W,L U OZC,e t a l．M o h o t o p o g r a p h y o f t h eT iＧ
b e t a nP l a t e a uu s i n g m u l t iＧs
c a l e g r a v i t y a n a l y s i sa n di t st e cＧ
t o n i c i m p l i c a t i o n s[J]．J o u r n a l o fA s i a nE a r t hS c i e n c e s,２０１７,
１３８:３７８Ｇ３８６．
[１６]㊀牟力,陈召曦．重力资料多尺度分析最优小波基的选择[J]．物探与化探,２０１５,３９(５):１０１３Ｇ１０１９．
MO U L i,C H E NZ h a o x i．T h eo p t i m a l c h o i c eo fw a v e l e t b a s e s
i n g r a v i t y d a t am u l t iＧs c a l ea n a l y s i s[J]．G e o p h y s i c a l a n dG e oＧ
c h e m i c a l E x p l o r a t i o n,２０１５,３９(５):１０１３Ｇ１０１９．
[１７]㊀刁博,王家林,程顺有．重力异常小波多分辨分析分解阶次的确定[J]．地球科学,２００７,３２(４):５６４Ｇ５６８．
D I A OB o,WA N GJ i a l i n,C H
E N GS h u n y o u．T h e c o n f i r m a t i o n
o f d e c o m p o s i t i o n l e v e l i nw a v e l e tm u l t iＧr e s o l u t i o n a n a l y s i s f o r
g r a v i t y a n o m a l i e s[J]．E a r t hS c i e n c e,２００７,３２(４):５６４Ｇ５６８．[１８]㊀P A R K E R R L．T h er a p i dc a l c u l a t i o no f p o t e n t i a l a n o m a l i e s [J]．G e o p h y s i c a l J o u r n a l o f t h eR o y a lA s t r o n o m i c a l S o c i e t y,
１９７３,３１(４):４４７Ｇ４５５．
[１９]㊀O L D E N B U R GD W．T h e i n v e r s i o n a n d i n t e r p r e t a t i o n o f g r a vＧ
i t y a n o m a l i e s[J]．G e o p h y s i c s,１９７４,３９(４):５２６Ｇ５３６．[２０]㊀熊小松,高锐,张兴洲,等．深地震探测揭示的华北及东北地区莫霍面深度[J]．地球学报,２０１１,３２(１):４６Ｇ５６．
X I O N G X i a o s o n g,G A O R u i,Z H A N G X i n g z h o u,e ta l．T h e
M o h od e p t ho fN o r t hC h i n aa n dn o r t h e a s tC h i n a r e v e a l e db y
s e i s m i c d e t e c t i o n[J]．A c t aG e o s c i e n t i c aS i n i c a,２０１１,３２(１):
４６Ｇ５６．
[２１]㊀雷建平．南海南部重力异常特征及深部构造研究[D]．青岛:中国海洋大学,２０１６．
L E I J i a n p i n g．C h a r a c t e r i s t i c s o f g r a v i t y a n o m a l y a n dd e e p t e cＧ
t o n i c s i nt h es o u t h e r nS o u t hC h i n aS e a[D]．Q i n g d a o:O c e a n
U n i v e r s i t y o fC h i n a,２０１６．
[２２]㊀陈晓慧．台湾海峡布格重力异常及构造解释[D]．青岛:中国海洋大学,２０１５．
C H E N X i a o h u i．B o u g u e r g r a v i t y a n o m a l y i nT a i w a nP r o v i n c e
S t r a i t a n d i t s t e c t o n i c i n t e r p r e t a t i o n[D]．Q i n g d a o:O c e a nU n iＧ
v e r s i t y o fC h i n a,２０１５．
[２３]㊀郭雨帆,雷建平,胡珊,等．华南地块及邻区重力异常特征及结构研究[J]．大地测量与地球动力学,２０１７,３７(增刊２):１６１Ｇ
１６７．
G U O Y u f a n,L E I J i a n p i n g,HUS h a n,e t a l．R e s e a r c ho n g r a v iＧ
t y a n o m a l y c h a r a c t e r i s t i c s a n d s t r u c t u r e o f S o u t hC h i n a b l o c k
[J]．J o u r n a l o f G e o d e s y a n d G e o d y n a m i c s,２０１７,３７
(S u p p l０２):１６１Ｇ１６７．
[２４]㊀周稳生．四川盆地重磁异常特征与深部结构[D]．南京:南京大学,２０１６．
Z H O U W e n s h e n g．C h a r a c t e r i s t i c s a n d d e e p s t r u c t u r e o f g r a v iＧ
t y a n d m a g n e t i ca n o m a l i e si n S i c h u a n B a s i n[D]．N a n j i n g:
N a n j i n g U n i v e r s i t y,２０１６．
[２５]㊀G A O DZ,H O U ZZ,T A N GJ．M u l t i s c a l e a n a l y s i so f g r a v i t y
a n o m a l i e s o nE a s tC h i n aS e a a n d a d j a c e n t r e g i o n s[J]．C h i n e s e
J o u r n a l o fG e o p h y s i c s,２０００,４３(６):８８９Ｇ８９６．
[２６]㊀谢天峰,付永涛．南黄海布格重力异常的小波多尺度分析[J]．海洋科学,２００８,３２(１１):６０Ｇ６３．
X I E T i a n f e n g,F U Y o n g t a o．W a v e l e tt r a n s f o r m a n d m u l t iＧ
s c a l ed e c o m p o s i t i o n o f B o u g u e r g r a v i t y d a t a o ft h e S o u t h
Y e l l o wS e a[J]．M a r i n eS c i e n c e s,２００８,３２(１１):６０Ｇ６３．[２７]㊀L E I JP,J I A N G S H,L ISZ,e ta l．G r a v i t y a n o m a l y i nt h e s o u t h e r nS o u t hC h i n aS e a:a c o n n e c t i o no fM o h o d e p t h t o t h e
n a t u r eo f t h es e d i m e n t a r y b a s i n s＇c r u s t[J]．G e o l o g i c a l J o u rＧ
n a l,２０１６,５１:２４４Ｇ２６２．
[２８]㊀S P E C T O R A,G R A N TFS．S t a t i s t i c a lm o d e l s f o r i n t e r p r e t i n g
a e r o m a g n e t i c d a t a[J]．G e o p h y s i c s,１９７０,３５(２):２９３Ｇ３０２．[２９]㊀X U C,L U O ZC,S U N R,e ta l．M u l t i l a y e rd e n s i t i e su s i n g a w a v e l e tＧ
b a s e d g r a v i t y m e t h o da n dt h e i r t e
c t o n i c i m p l i c a t i o n s
b e n e a t h t h eT i b e t a nP l a t e a u[J]．G e o p h y s i
c a l J o u r n a l I n t e r n aＧ
t i o n a l,２０１８,２１３(３):２０８５Ｇ２０９５．
[３０]㊀黎晋博,许闯,简光煜,等．中国南海重力梯度多尺度分析[J]．测绘科学,２０２０,４５(１１):１Ｇ７．
L I J i n b o,X UC h u a n g,J I A N G u a n g y u,e t a l．M u l t iＧs c a l e a n a l yＧ
s i s o f g r a v i t yg r a d i e n t s i nS o u t hC h i n aS e a[J]．S c i e n c e o f S u rＧ
v e y i n g a n d M a p p i n g,２０２０,４５(１１):１Ｇ７．
[３１]㊀张永奇,韩美涛,张恩会,等．渭河盆地及邻区重力异常小波多尺度分解与解释[J]．地震研究,２０２２,４５(１):７５Ｇ８７．
Z H A N G Y o n g q i,H A N M e i t a o,Z H A N G E n h u i,e ta l．M u l t iＧ
s c a l e w a v e l e t d e c o m p o s i t i o n a n di n t e r p r e t a t i o n o f g r a v i t y
３８０１
第４５卷第５期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀郭雨帆,等:华北地块及邻区重力异常特征及地质结构㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀。