Construction of bedrock geological model in shallow overburden area: Application of a joint inversion of gravity anomaly and borehole data in the Tanlu Fault Zone, Sihong area
-
摘要:
利用重力异常数据建立覆盖区的基岩地质模型往往需要大量的地质与岩石物性资料，同时受地质先验信息、人为因素等干扰严重。运用重力异常钻孔约束反演的方法尝试构建浅覆盖区基岩地质模型，探讨在人为介入较少的情况下快速高效构建包含基岩面起伏特征、断裂展布、密度分布等要素的地质模型的可行性。在郯庐断裂泗洪地区的地质调查实际应用中，对区内布格重力异常采用小波多尺度分解的方法进行基岩面异常场提取，随后联合运用小波断裂分析、线性回归反演、视密度填图等方法构建了该区域的隐伏基岩地质模型。结果显示，在不需要大量地质先验信息的情况下，该系列方法能够快速高效地构建覆盖区地质模型。结果模型能够较好地反映基岩面构造及物性特征，并且细化郯庐断裂带泗洪地区的区域构造特征，对于郯庐断裂地区的构造演化研究具有一定的参考意义，同时为探索较少人为干预的地球物理数据反演覆盖区隐伏基岩的方法研究提供一定的实践价值。
Abstract:The establishment of the bedrock geological model with gravity anomaly data in the covered area usually requires a lot of geological and rock-physical data.The results are often disturbed by geological prior information and human factors.A joint inversion method of gravity anomaly and borehole data for constructing the geological bedrock model in shallow overburden area was introduced to explore the feasibility of rapidly and efficiently constructing the geological model of bedrock surface undulation, fault distribution and density distribution in the case of less human intervention.In the practical application of the Tanlu Fault Zone, Sihong area, the wavelet multi-scale decomposition method was used to extract the gravity anomaly field of bedrock surface; furthermore, the geological model of bedrock in this region was constructed by means of wavelet fault analysis, linear regression inversion and apparent density mapping.The results show that this method can construct the geological model of overburden area quickly and efficiently without requiring a lot of geological prior information.This model can truthfully reflect the structural and physical characteristics of the bedrock surface, and refine the regional tectonic characteristics of the Tanlu Fault, which is meaningful for the better understanding of the regional tectonic evolution and geophysical data inversion of buried bedrock.
-
浅覆盖区覆盖层的地质调查和研究与深部隐伏基岩的物性、构造特征研究密不可分。一方面，深部基岩面构造特征制约了上覆盖层的沉积演化，沉积层中的同沉积构造变形可进一步细化区域构造演化史，另一方面，沉积物也可以通过侵蚀、搬运、负载等方式对基岩的构造特征进行改造，二者是一对耦合关系。常规的地表地质调查很难获得隐伏基岩及断裂的特征信息，故本次研究采用地质钻探与重力异常相结合的方法进行重力异常钻孔约束反演，重力位场数据横向分辨能力较高，地质钻孔数据可提供岩石物性特征并进行深度约束，因此，可有效减少重力异常数据地质解释的多解性，构建较准确的基岩地质模型。研究区郯庐断裂带泗洪地区距离曾发生8.5级地震的郯城仅130余千米，由于该区域被第四纪—新近纪沉积物广为覆盖，前人对其深部构造特征的研究较薄弱，故本次研究结合区域第四纪—新近纪覆盖层沉积结构，探讨断裂的空间展布特征，具有一定科学意义。
1. 区域地质及基岩物性特征
1.1 区域地质背景
郯庐断裂带是贯穿中国东部大陆边缘的一条NNE向走滑剪切断裂带。研究表明，这条断裂在晚新生代仍然活动，是一条重要的地震活动带，其分段特征明显，从南到北可分为苏皖段、山东段及东北段3段[1-2], 其中苏皖段途经地区人口稠密、经济发达，受到国内外学者的广泛关注。中生代以来，郯庐断裂带经历了华南-华北板块碰撞、太平洋板块俯冲，以及白垩纪、古近纪强烈的地幔热作用等一系列重大地质事件，构造演化复杂[3-5]。第四纪以来，郯庐断裂带依然活动频繁，发生过多次中强级以上地震，如1668年莒县-郯城8.5级地震，1888年渤海湾7.5级地震，1969年渤海7.4级地震，1975年海城7.3级地震等，其中1668年莒县-郯城8.5级地震是中国东部有记录的最强烈的地震。郯庐断裂带泗洪地区位于郯城以南约130 km，归属郯庐断裂带苏皖段的中北部，位于华北板块、苏鲁造山带及扬子板块的结合部位(图 1-a)，该地区被第四纪—新近纪沉积物大面积覆盖(图 1-b)，基岩及断裂以隐伏断裂为主，仅在重岗山、朱山附近有基岩出露。前人认为，郯庐断裂带在泗洪地区主要由山左口-泗洪断裂(F1)、新沂-新店断裂(F2)、墨河-凌城断裂(F3)、王集-朱山断裂(F4)和桥北镇-宿迁断裂(F5)5条分支断裂延伸组成，宽20~30 km，其中F2、F3和F4分别由山东段的安丘-莒县断裂、沂水-汤头断裂和鄌郚-葛沟断裂延续而来[1, 4, 6-7]。该区域性断裂带西靠华北板块，东临苏鲁造山带，同时受到邻近扬子板块挤压、东侧太平洋板块俯冲及远端喜马拉雅抬升的共同作用，构造演化极其复杂，前人将郯庐断裂带的构造演化总结为4个阶段，分别为早中生代左旋转换走滑，晚侏罗世—早白垩世早期左行平移，早白垩世—古近纪大规模伸展正断及新近纪以来挤压逆冲兼具右行走滑[8-18]。不同区段内部受物性结构及动力学背景差异影响，构造演化特征存在偏差。根据采取的研究手段精度及研究时间尺度的不同，对郯庐断裂带构造演化仍存在争议，如万天丰等[19]将其演化总结为：1次左行走滑(250~208 Ma)；2次正断层(135~52 Ma与23.3~0.73 Ma)；3次挤压逆冲断层(208~135 Ma、52~23.3 Ma、0.73 Ma以来)。而张岳桥等[20]将郯庐断裂带中生代的运动学演化概括为“两大运动时期、五个发展阶段”。Li等[21]的研究显示，郯庐断裂带泗洪地区缺失古近纪地层，指示泗洪地区在晚白垩世末期—古近纪经历挤压抬升，同时泗洪县南部的反射地震剖面显示其下伏断裂为第四纪正断层，说明郯庐断裂带部分地区在第四纪仍处于伸展状态。
![]()
图 1 郯庐断裂带构造分区(a)和泗洪地区地质概况(b)NCB—华北板块；SLO—苏鲁造山带；YZB—扬子板块；HXB—华夏板块；DBO—大别造山带；TLFZ—郯庐断裂带；JXF—嘉山-响水断裂；ZTF—镇江-铜陵断裂；JSF—江绍断裂Figure 1. Tectonic division of the Tanlu Fault Zone (a)and geological map of the Sihong area (b)1.2 基岩地层及物性特征
泗洪地区仅沿郯庐断裂带零星分布的孤丘可见基岩露头，为中—新元古代地层组成的丘陵岗地，其余大部分地区被第四纪松散层覆盖，其沉积厚度差异较大、地层发育不全。
基岩地层单元以郯庐断裂带为界，划分为3个大区，西侧属华北地层区，在朱山附近出露新元古代八公山群兰陵组(Pt3l)及九顶山组上段(Pt3jd2)，分别为石英岩及微晶灰岩；断裂带内为晚白垩世红色砂岩、砂砾岩，在重岗山地区出露，岩石普遍风化较强，砂砾岩中砾石主要为片麻岩巨砾、漂砾，且沿重岗山东、西两侧分别可见其逆冲于新近系、第四系之上，断层走向近正北；东侧为苏鲁造山带锦屏岩群，地表未见露头，钻孔揭露显示基底为花岗片麻岩。据2016—2018年研究区完成的29个地质钻孔岩心特征，研究区普遍缺失古近系。基岩及覆盖层时代特征及填图单位分布如表 1所示。
表 1 基岩填图单位Table 1. The mapping units of the bedrock中新生代盖层 华北板块 苏鲁造山带 时代 古近纪(缺失) 白垩纪 新元古代 新元古代—中元古代 中元古代—古元古代 地层及代号 三垛组(E2-3s) 戴南组(E2d) 阜宁组(E1f) 泰州组(E1t) 王氏组(K2w) 青山群(K1qs) 淮河群 八公山群 云台岩群(Pt2-3Y) 锦屏岩群 九顶山组(Pt3jd) 张渠组(Pt3zh) 魏集组(Pt3wj) 兰陵组(Pt3l) 石英岩片麻岩组(Pt1-2Jqg) 磷灰岩片麻岩组(Pt1-2Jpg) 岩石物性资料来源于中国地质调查局南京地质调查中心的实测物性成果资料。在重岗山、柳山、朱山、大庄村基岩出露点用取样机进行取心，通过对取心样品和钻孔内基岩样品进行物性测量，得到样品密度参数。本次物性测量共采集密度标本191块，野外采集中新生代、新元古代露头143块，采集钻孔标本48块，经加工为圆柱体后对标本的密度进行测定(表 2)。
表 2 泗洪地区部分岩石密度参数Table 2. The densities of the bedrock in the Sihong area年代 岩性 标本数/块 密度/(103 kg·m-3) 最小值 最大值 平均值 标准差 白垩纪 粉砂岩 13 2.033 2.346 2.191 0.131 白垩纪 砂岩 43 2.243 2.535 2.330 0.048 白垩纪 火山岩 23 2.254 2.368 2.330 0.028 元古宙 花岗片麻岩 19 2.426 2.669 2.585 0.065 元古宙 片麻岩 24 2.545 2.806 2.658 0.062 元古宙 灰岩 32 2.616 2.857 2.800 0.050 元古宙 石英岩 34 2.598 2.676 2.643 0.014 元古宙 斜长角闪岩 3 2.456 2.616 2.552 0.069 从岩石密度统计表(表 2)可以看出，元古宙地层密度高于中新生代盖层密度，元古宙地层中灰岩密度最高，均值为2.800×103 kg/m3，其次为片麻岩和石英岩，密度均值分别为2.658×103 kg/m3、2.643×103 kg/m3，斜长角山岩密度均值为2.552×103 kg/m3，花岗片麻岩密度均值为2.585×103 kg/m3。结果显示，中新生代地层密度值在2.191×103~2.330×103 kg/m3之间。
2. 数据及处理
郯庐断裂带泗洪地区及其邻区重力异常数据来源于中国地质调查局南京地质调查中心，数据网格化间隔为1500 m，参与处理与反演的重力异常均为布格重力异常。
本次重力异常所选区域略大于研究区，泗洪地区及其邻区布格重力异常特征如图 2所示，从NWW至SEE体现了较明显的分区特征：郯庐断裂带在异常图中显示出明显的NNE向线性梯度带、串珠状异常特征，异常规模较大，宽约30 km且纵贯整个区域中部，断裂带内部以中等正异常为主。以此条中等正异常串珠状条带为界，西侧为低重力起伏区，东侧为高重力起伏区，东南角则为高重力平缓区，分别对应华北板块、苏鲁造山带及扬子板块的基底特征。
![]()
图 2 郯庐断裂带泗洪地区及邻区布格重力异常特征Figure 2. Bouguer gravity anomaly in the Sihong area of the Tanlu Fault Zone区域布格异常显示该区的主导构造应力方向与构造单元分布情况，郯庐断裂带泗洪地区布格异常指示构造特征方向为NNE向，内部主要包括华北地块、郯庐断裂带及苏鲁造山带。郯庐断裂带在泗洪地区各分支断裂的展布形态在布格异常图中不明显，仅重岗山西侧、朱山东侧发育的梯级带较好地指示了F2和F4断裂的发育特征。
由于重力异常数据是不同深度的异常场源在某一观测面的叠加，因此区域异常场和局部异常场的分离是异常反演和解释的关键性问题之一[22]。随着近年场源分离技术方法及相关反演理论的发展，重力异常数据解释也越来越多地被应用于覆盖区隐伏基岩的调查与研究中[23-25]，常用的分离区域异常场方法有多项式曲面拟合法[26]、最小曲率法[27]、有限元分析[28]、小波多尺度分解[29-32]等，本次选用小波多尺度分解与对数功率谱相结合方法进行重力异常的位场分离及深度估算。
小波分析可将位场信号分解到不同的尺度空间，并且通过伸缩、平移聚焦到信号的任一细节进行分析[29]。根据小波多尺度分解方法原理[30-32]，选取db11小波基[33-34]，采用Mallat重构算法[35]，重力异常场Δf(x, y)的一阶分解为：
Δf(x,y)=A1Δf(x,y)+Dh1Δf(x,y)+Dv1Δf(x,y)+Dd1Δf(x,y) (1) 其中，A1Δf(x, y)代表分解中的低频部分，Dh1Δf(x,y),Dv1Δf(x,y),Dd1Δf(x,y)为第一阶分解中的EW(水平)、SN(垂直)和NE-SW(对角线)方向上的高频细节部分，可用来进行断裂识别，在这里将公式(1)简写为:
Δf(x,y)=A1Δf(x,y)+D1Δf(x,y) (2) 再对A1Δf(x, y)进行分解得到：
A1Δf(x,y)=A2Δf(x,y)+D2Δf(x,y) (3) 以此类推，可将Δf(x, y)分解为:
Δf(x,y)=AnΔf(x,y)+D1Δf(x,y)+D2Δf(x,y)+D3Δf(x,y)+⋯+DnΔf(x,y) (4) 其中，AnΔf(x, y)为第n阶小波变换的逼近部分，DnΔf(x, y)为第n阶小波变换高频细节部分。
通过对泗洪地区布格异常进行2阶小波分解，得到1、2阶小波细节异常(D1G和D2G)及2阶小波逼近异常(A2G)(图 3-a~c)。随后运用对数功率谱法[36]对研究区的各阶重力异常细节进行场源深度估算(图 3-d、e)，得到D1G和D2G反映的平均场源深度约为1 km和2 km。
![]()
图 3 郯庐断裂带泗洪地区及其邻区布格重力异常二阶小波分解结果(a~c)及一、二阶小波细节异常功率谱测深结果(d, e)Figure 3. The second order wavelet multi-scale decomposition of the bouguer gravity anomaly in the Sihong area of Tanlu Fault Zone (a~c)and the logarithmic power spectrum of the first and second order gravity anomaly details (d, e)根据对重力小波细节异常场源深度的估算，D1G和D2G可以反映泗洪地区地壳顶部不同密度体引起的重力异常，如隐伏基岩不同年代地层的分布和呈线性特征的断裂等，其中D1G的不均一性较明显，且场源深度较浅，可以作为基岩面顶部物性分布及起伏特征的反映，而D2G的线性特征更突出，场源深度更深，更适合研究基岩断裂的平面展布特征；A2G则较好地反映了区域构造分区特征，其特征显示出华北板块与苏鲁造山带的物性差异。结合泗洪地区钻孔揭露的基岩面平均深度(约200 m)，本次选取D1G进行隐伏基岩埋深和密度特征反演，同时对一、二阶小波细节异常的水平、垂直和对角线方向上的高频细节特征(据公式4)进行分析，推测区内基岩断裂的展布特征。
由于研究区构造线走向主要为NNE—近SN向，郯庐断裂带也以此走向为主，故本区的断裂划分以D1G和D2G在SN向上的分量特征作为主要依据，结合对角线方向与东西向的细节特征，勾画研究区断裂展布特征。结果显示，郯庐断裂带在泗洪地区的5条分支断裂在SN向特征分量上(图 4-b、e)呈明显的近SN向连续的重力梯级带，同时在近NW—SE向上(图 4-a、d)也有许多串珠状异常呈线性排布，指示一系列走滑剪切断裂发育。该走滑断裂系在研究区中部与SN向和NNE向的断裂相互交切，形成了较复杂的断裂体系，且近EW向与NW向断裂多切割NNE—NE向的郯庐断裂带各分支断裂。据此推断，本区存在多期次的构造叠加，郯庐断裂NNE向构造体系被较新的NW向构造运动改造，形成现今的断裂体系。
![]()
图 4 郯庐断裂带泗洪地区小波断裂分析结果a~c—重力一阶小波细节分别在NE-SW向、SN向和EW向的特征；d~f—重力二阶小波细节分别在NE—SW向、SN向和EW向的特征; TLFZ—郯庐断裂带；F1—山左口-泗洪断裂；F2—新沂-新店断裂；F3—墨河-淩城断裂；F4—王集-朱山断裂；F5—桥北镇-宿迁断裂；F6—太平-孙园断裂；F7—界集-太平断裂；F8~F10—局部走滑断裂Figure 4. Wavelet analysis of the fault distribution in the Sihong area of the Tanlu Fault Zone3. 结果与讨论
3.1 基岩面埋深及密度分布
根据线性回归近似反演理论[37-38]，如果密度界面起伏平缓，可以认为重力变化与密度界面的起伏近似呈线性关系，即：
h=a+b⋅Δg (5) 式中，h为界面深度，Δg为界面起伏引起的重力异常，a，b为2个系数，它们与重力异常起算点处界面深度和界面上下层物质层的密度差有关。
如果已有n个已知基岩面深度，设其深度为hi(i=1, 2, ..., n)，根据最小二乘原理，为了确定a，b，应使得各点深度hi与计算出的深度hi偏差的平方和最小，即：
ϕ(a,b)=n∑i=1(h−ˉhi)=min (6) 令\partial \varphi /\partial a、\partial \varphi /\partial b分别等于零，可得：
\sum {{{\rm{h}}_{\rm{i}}}} - {\rm{na}} - \sum {\rm{b}} \Delta {g_i} = 0 (7) \sum {{h_i}} \Delta {g_i} - \sum a \Delta {g_i} - \sum b \Delta g_i^2 = 0 (8) 两式联立，解之得：
a = \frac{{\sum \Delta {g_i}^2\sum {{h_i}} - \sum \Delta {g_i}\sum \Delta {g_i}{h_i}}}{{n\sum \Delta {g_i}^2 - {{\left( {\sum \Delta {g_i}} \right)}^2}}} (9) b = \frac{{n\sum {\Delta {g_i}} {h_i} - \sum {{h_i}} \sum \Delta {g_i}{\rm{ }}}}{{n\sum {\Delta {g_i}^2} - {{\left( {\sum \Delta {g_i}} \right)}^2}}} (10) a，b求出后，则可根据公式(5)计算出各测点下方密度界面的深度。
本次研究应用前期完成的26个地质钻孔的基岩面埋深数据及朱山基岩露头高程数据作为已知数据，钻孔位置见图 1-b，其对应的重力一阶小波细节异常值如表 3所示。
表 3 用于线性回归反演的地质钻孔信息Table 3. The geological boreholes for linear regression inversion序号 钻孔编号及基岩出露点 钻孔高程/m 基岩面高程/m Δg值/(10-6 m·s-2) 1 SHJ01 18.88 -41.4 0.082 2 SHJ03 18.47 -234.11 -0.097 3 SHJ05 18.12 -232.78 -0.061 4 SHJ07 22.47 -146.31 -0.034 5 SHJ08 19.59 -195.38 -0.062 6 SHJ09 18.71 -206.85 0.037 7 SHJ10 18.69 -220 0.086 8 SHJ11 18.13 -223.4 -0.086 9 SHJ14 16.37 -154 -0.207 10 SHJ15 16.71 -113.6 -0.188 11 SHJ17 16.09 -83.65 0.005 12 SHJ18 18.82 -171.8 0.026 13 SHJ19 46.04 -162.36 0.013 14 SHJ20 22.77 -158.72 -0.021 15 SHJ21 16.47 -77.86 0.087 16 SHJ22 15.46 -117 0.009 17 SHJ23 13.02 -194.1 -0.065 18 SHJ24 20.42 -254.77 -0.209 19 SHJ25 20.86 -165.05 -0.087 20 SHJ26 18.39 -125.67 0.057 21 SHJ28 15.32 -184.8 0.051 22 SHJ31 16.73 -103.49 -0.096 23 SHJ34 15.03 -106.3 -0.036 24 SHJ35 19.63 -155.18 0.088 25 SHJ37 16.38 -139.25 0.017 26 SHJ38 26.07 -108 -0.072 27 朱山 / 157.2 0.458 从钻孔基岩面高程数据(图 5-a)与对应重力小波一阶细节异常数据(图 5-b)的分布特征图中，可以看出二者具有一定的线性相关性。进一步依据相关系数的计算公式：
{R_{j, z}} = \frac{{{\rm{cov}}\left( {j, z} \right)}}{{{S_j} \times {S_z}}} (11) ![]()
图 5 钻孔揭露的基岩面高程(a)、钻孔处重力小波一阶细节异常数据(b)及二者相关性分析(c)Figure 5. Elevations of bedrocks exposed by boreholes (a), the first order detail data of the gravity anomaly in the location of boreholes (b)and the correlation analysis between them (c)计算得出二者的相关系数为0.684，属于中度相关(图 5-c)，可以进行线性回归分析[39]，其中cov(j, z)为基岩面高程数据与重力异常数据的协方差，Sj和Sz分别为基岩面高程数据和重力异常数据的标准差。
据公式(6)~(10)得到基岩面高程与重力一阶小波细节异常线性方程为：
h = 433.49 \times \Delta g - 140.24 (12) 根据线性回归方程(公式(12))对重力一阶小波细节异常未知基岩面高程数据点进行计算，得到研究区基岩面高程分布特征(图 6-a)：研究区的基岩波动起伏较明显，起伏范围为-220~20 m，西北角及中部偏西位置为2个较明显的隆起，走向NNE向，与朱山、重岗山位置较吻合，2处基岩出露岩性分别为新元古代白云质灰岩及晚白垩世泥质粉砂岩，西侧由F1~F5组成的郯庐断裂带内基岩面埋深起伏较大，形态上呈明显的坳隆相间分布，构造上为地堑地垒交替出现样式，朱山东侧和重岗山附近均发育NNE向展布的地堑，朱山附近、上塘—黑塔东一线及重岗山东侧为3处地垒构造。据基岩起伏特征及小波断裂分析结果可勾画出区内基岩断裂的发育情况。
![]()
图 6 郯庐断裂带泗洪地区隐伏基岩地质模型a—隐伏基岩起伏特征及基底断裂分布；b—基岩面附近岩石视密度分布Figure 6. Geological model of buried bedrocks in the Sihong area, Tanlu Fault本次采用频率域视密度填图方法[40]对基底密度分布特征进行反演，用于反演的重力异常为重力一阶小波细节异常，选取本次线性回归反演得到的基岩面深度为密度层上界面深度，下界面深度选为对数功率谱标定的D1G平均深度1 km，平均密度设为2.3×103 kg/m3，经6次迭代得到近基岩面深度视密度分布特征(图 6-b)。郯庐断裂带泗洪地区地堑构造内为明显的低密度特征，密度值在2.14×103~2.26×103 kg/m3之间，与晚白垩世粉砂岩、砂岩的密度特征吻合；西北角朱山地垒为明显的高密度特征，密度值在2.42×103~2.6×103 kg/m3之间，与新元古代石英岩物性特征相似；而上塘-黑塔地垒与重岗山东地垒体现出的高密度体可能与断裂带内深部岩体侵入或古老变质结晶基底残留有关；东部苏鲁造山带内则总体为高密度分布，局部凹陷内表现为低密度特征，对比低密度区钻孔岩心特征，其基岩岩性为晚白垩世粉砂岩，结合区内断裂分布解释结果，说明苏鲁造山带内的晚白垩世地层为局部断陷沉积。
3.2 讨论
本文依据钻孔数据，对目标重力异常场界面起伏特征进行了反演，得到泗洪地区隐伏基岩面起伏特征与Jiang等[41]总结的区域构造地质特征较吻合，同时结合重力异常小波断裂分析、视密度填图等方法对断裂展布和岩性特征分布进行了解释，得到郯庐断裂带泗洪地区较精细的隐伏基岩地质模型。该模型具有一定的可信度，但受限于重力异常数据处理与解释过程中的多解性及反演方法基本假定的不精确性，该模型仍存在一些不足。首先，受位场分离方法精度限制，提取的小波一阶细节异常可近似反演隐伏基岩异常特征，精度不足；其次，线性回归反演理论假设重力变化与密度界面的起伏近似呈线性关系，忽略了重力异常变化一般与密度正相关这一因素，故在基岩岩石密度较高地区，界面深度反演结果偏深(或矮)，而在基岩密度较低区，深度偏浅(或高)。
4. 结论
(1) 已知钻孔约束下的重力异常反演方法可以在较少人为干预下，较高效地构建浅覆盖区隐伏基岩的地质模型，能够在覆盖区的地质填图中起到深部基岩构造的辅助解释作用，但其精度受制于位场分离方法与线性回归算法的粗糙性，仍有待进一步提升。
(2) 郯庐断裂带泗洪地区基岩地质模型显示：①泗洪地区郯庐断裂带的5条分支断裂在研究区中部呈NNE向平行排列，且发现区内多条NW向左行走滑断裂发育，并切割郯庐断裂带，说明该区域受到新的构造应力叠加；②郯庐断裂带在泗洪地区的隐伏基岩面呈“凹-隆-凹-隆”起伏特征; ③基岩视密度填图结果显示, 泗洪地区郯庐断裂带内的基岩物性特征存在较大的不均一性，高密度体物性特征与苏鲁造山带锦屏岩群相似，但断裂带内岩浆活动频繁，也不排除其为岩浆侵入成因。
致谢: 感谢两位审稿专家提出的宝贵意见及中国地质调查局地质力学研究所胡健民研究员在本文写作过程中的帮助。 -
![]()
图 1 郯庐断裂带构造分区(a)和泗洪地区地质概况(b)
NCB—华北板块；SLO—苏鲁造山带；YZB—扬子板块；HXB—华夏板块；DBO—大别造山带；TLFZ—郯庐断裂带；JXF—嘉山-响水断裂；ZTF—镇江-铜陵断裂；JSF—江绍断裂
Figure 1. Tectonic division of the Tanlu Fault Zone (a)and geological map of the Sihong area (b)
![]()
图 2 郯庐断裂带泗洪地区及邻区布格重力异常特征
Figure 2. Bouguer gravity anomaly in the Sihong area of the Tanlu Fault Zone
![]()
图 3 郯庐断裂带泗洪地区及其邻区布格重力异常二阶小波分解结果(a~c)及一、二阶小波细节异常功率谱测深结果(d, e)
Figure 3. The second order wavelet multi-scale decomposition of the bouguer gravity anomaly in the Sihong area of Tanlu Fault Zone (a~c)and the logarithmic power spectrum of the first and second order gravity anomaly details (d, e)
![]()
图 4 郯庐断裂带泗洪地区小波断裂分析结果
a~c—重力一阶小波细节分别在NE-SW向、SN向和EW向的特征；d~f—重力二阶小波细节分别在NE—SW向、SN向和EW向的特征; TLFZ—郯庐断裂带；F1—山左口-泗洪断裂；F2—新沂-新店断裂；F3—墨河-淩城断裂；F4—王集-朱山断裂；F5—桥北镇-宿迁断裂；F6—太平-孙园断裂；F7—界集-太平断裂；F8~F10—局部走滑断裂
Figure 4. Wavelet analysis of the fault distribution in the Sihong area of the Tanlu Fault Zone
![]()
图 5 钻孔揭露的基岩面高程(a)、钻孔处重力小波一阶细节异常数据(b)及二者相关性分析(c)
Figure 5. Elevations of bedrocks exposed by boreholes (a), the first order detail data of the gravity anomaly in the location of boreholes (b)and the correlation analysis between them (c)
![]()
图 6 郯庐断裂带泗洪地区隐伏基岩地质模型
a—隐伏基岩起伏特征及基底断裂分布；b—基岩面附近岩石视密度分布
Figure 6. Geological model of buried bedrocks in the Sihong area, Tanlu Fault
表 1 基岩填图单位
Table 1 The mapping units of the bedrock
中新生代盖层 华北板块 苏鲁造山带 时代 古近纪(缺失) 白垩纪 新元古代 新元古代—中元古代 中元古代—古元古代 地层及代号 三垛组(E2-3s) 戴南组(E2d) 阜宁组(E1f) 泰州组(E1t) 王氏组(K2w) 青山群(K1qs) 淮河群 八公山群 云台岩群(Pt2-3Y) 锦屏岩群 九顶山组(Pt3jd) 张渠组(Pt3zh) 魏集组(Pt3wj) 兰陵组(Pt3l) 石英岩片麻岩组(Pt1-2Jqg) 磷灰岩片麻岩组(Pt1-2Jpg) 
下载: 导出CSV
表 2 泗洪地区部分岩石密度参数
Table 2 The densities of the bedrock in the Sihong area
年代 岩性 标本数/块 密度/(103 kg·m-3) 最小值 最大值 平均值 标准差 白垩纪 粉砂岩 13 2.033 2.346 2.191 0.131 白垩纪 砂岩 43 2.243 2.535 2.330 0.048 白垩纪 火山岩 23 2.254 2.368 2.330 0.028 元古宙 花岗片麻岩 19 2.426 2.669 2.585 0.065 元古宙 片麻岩 24 2.545 2.806 2.658 0.062 元古宙 灰岩 32 2.616 2.857 2.800 0.050 元古宙 石英岩 34 2.598 2.676 2.643 0.014 元古宙 斜长角闪岩 3 2.456 2.616 2.552 0.069
下载: 导出CSV
表 3 用于线性回归反演的地质钻孔信息
Table 3 The geological boreholes for linear regression inversion
序号 钻孔编号及基岩出露点 钻孔高程/m 基岩面高程/m Δg值/(10-6 m·s-2) 1 SHJ01 18.88 -41.4 0.082 2 SHJ03 18.47 -234.11 -0.097 3 SHJ05 18.12 -232.78 -0.061 4 SHJ07 22.47 -146.31 -0.034 5 SHJ08 19.59 -195.38 -0.062 6 SHJ09 18.71 -206.85 0.037 7 SHJ10 18.69 -220 0.086 8 SHJ11 18.13 -223.4 -0.086 9 SHJ14 16.37 -154 -0.207 10 SHJ15 16.71 -113.6 -0.188 11 SHJ17 16.09 -83.65 0.005 12 SHJ18 18.82 -171.8 0.026 13 SHJ19 46.04 -162.36 0.013 14 SHJ20 22.77 -158.72 -0.021 15 SHJ21 16.47 -77.86 0.087 16 SHJ22 15.46 -117 0.009 17 SHJ23 13.02 -194.1 -0.065 18 SHJ24 20.42 -254.77 -0.209 19 SHJ25 20.86 -165.05 -0.087 20 SHJ26 18.39 -125.67 0.057 21 SHJ28 15.32 -184.8 0.051 22 SHJ31 16.73 -103.49 -0.096 23 SHJ34 15.03 -106.3 -0.036 24 SHJ35 19.63 -155.18 0.088 25 SHJ37 16.38 -139.25 0.017 26 SHJ38 26.07 -108 -0.072 27 朱山 / 157.2 0.458
下载: 导出CSV
-
刘保金, 酆少英, 姬计法, 等. 郯庐断裂带中南段的岩石圈精细结构[J]. 地球物理学报, 2015, 58(5): 1610-1621. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201505013.htm 张鹏, 王良书, 钟锴, 等. 郯庐断裂带的分段性研究[J]. 地质论评, 2007, 53(5): 721-722. https://www.cnki.com.cn/Article/CJFDTOTAL-DZLP200705001.htm 徐嘉炜. 郯城-庐江平移断裂系统[J]. 构造地质论丛, 1984, 3: 19-22. 王小凤, 李中坚, 陈柏林, 等. 郯庐断裂带[M]. 北京: 地质出版社, 2000. 朱光, 王道轩, 刘国生, 等. 郯庐断裂带的演化及其对西太平洋板块运动的响应[J]. 地质科学, 2004, 40(1): 36-49. doi: 10.3321/j.issn:0563-5020.2004.01.005 张鹏, 李丽梅, 张景发, 等. 郯庐断裂带江苏段第四纪活动特征及其动力学背景探讨[J]. 防灾减灾工程学报, 2011, 31(4): 389-396. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201104008.htm 刘备, 朱光, 胡红雷, 等. 郯庐断裂带江苏段新构造活动规律分析[J]. 地质学报, 2015, 89(8): 1352-1366. doi: 10.3969/j.issn.0001-5717.2015.08.002 朱光, 刘国生, 牛漫兰, 等. 郯庐断裂带的平移运动与成因[J]. 地质通报, 2003, 22(3): 200-207. doi: 10.3969/j.issn.1671-2552.2003.03.009 朱光, 牛漫兰, 刘国生, 等. 郯庐断裂带早白垩世走滑运动中的构造、岩浆、沉积事件[J]. 地质学报, 2002, 76(3): 325-334. doi: 10.3321/j.issn:0001-5717.2002.03.005 朱光, 王道轩, 刘国生, 等. 郯庐断裂带的伸展活动及其动力学背景[J]. 地质科学, 2001, 36(3): 269-278. doi: 10.3321/j.issn:0563-5020.2001.03.002 朱光, 张力, 谢成龙, 等. 郯庐断裂带构造演化的同位素年代学制约[J]. 地质科学, 2009, 44(4): 1327-1342. doi: 10.3321/j.issn:0563-5020.2009.04.019 刘国生, 朱光, 牛漫兰, 等. 合肥盆地对郯庐断裂带同造山走滑活动的沉积响应[J]. 合肥工业大学学报(自然科学版), 2005, 28(10): 1233-1237. doi: 10.3969/j.issn.1003-5060.2005.10.001 刘国生, 朱光, 宋传中, 等. 郯庐断裂带新近纪以来的挤压构造与合肥盆地的反转[J]. 安徽地质, 2002, 12(2): 81-85. doi: 10.3969/j.issn.1005-6157.2002.02.001 侯明金, 朱光, Mercier J. 郯庐断裂带(安徽段)及邻区的动力学分析与区域构造演化[J]. 地质科学, 2007, 42(2): 362-381. doi: 10.3321/j.issn:0563-5020.2007.02.010 Grimmer J C, Jonckheere R, Enkelmann E, et al. Cretaceous-Cenozoic history of the southern Tan-Lu Fault Zone: apatite fission-track and structural constraints from the Dabie Shan(eastern China)[J]. Tectonophysics, 2002, 359(3): 225-253.
牛漫兰, 朱光, 刘国生. 郯庐断裂带中-南段中生代岩浆活动的构造背景与深部过程[J]. 地质科学, 2002, 37(4): 393-404. doi: 10.3321/j.issn:0563-5020.2002.04.002 孙卫东, 凌明星, 汪方跃, 等. 太平洋板块俯冲与中国东部中生代地质事件[J]. 矿物岩石地球化学通报, 2007, 27(s1): 218-225. https://www.cnki.com.cn/Article/CJFDTOTAL-KYDH200803003.htm 王勇生, 朱光, 宋传中, 等. 大别山东端郯庐断裂带由走滑向伸展运动转换的40Ar-39Ar年代学记录[J]. 地质科学, 2006, 41(2): 242-255. doi: 10.3321/j.issn:0563-5020.2006.02.007 万天丰, 朱鸿, 赵磊, 等. 郯庐断裂带的形成与演化: 综述[J]. 现代地质, 1996, 10(2): 159-168. https://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ602.002.htm 张岳桥, 董树文. 郯庐断裂带中生代构造演化史: 进展与新认识[J]. 地质通报, 2008, 27(9): 1371-1390. doi: 10.3969/j.issn.1671-2552.2008.09.002 Li C B, Jiang R, Zeng J W, et al. Deep structures underneath the Sihong Segment of the Tan-Lu Fault Zone, Eastern China: Interpretations of gravity anomaly and seismic profiles[J]. Journal of Asian Earth Sciences, 2019, 176: 229-243. doi: 10.1016/j.jseaes.2019.02.014
Sun B, Wang L, Dong P, et al. Integrated analysis on gravity and magnetic fields of the Hailar Basin, NE China: Implications for basement structure and deep tectonics[J]. Pure and Applied Geophysics, 2012, 169(11): 2011-2029. doi: 10.1007/s00024-012-0452-1
喻劲松, 荆磊, 王乔林, 等. 特殊地质地貌区填图物化探技术应用[J]. 地质力学学报, 2016, 22(4): 893-906. doi: 10.3969/j.issn.1006-6616.2016.04.008 刘彦, 严加永, 吴明安, 等. 基于重力异常分离方法寻找深部隐伏铁矿——以安徽泥河铁矿为例[J]. 地球物理学报, 2012, 55(12): 4181-4193. doi: 10.6038/j.issn.0001-5733.2012.12.030 邓震, 孟贵祥, 汤贺军, 等. 浅覆盖区1:5万基岩地质填图实践探索——以准噶尔北缘克什克涅绍喀尔(L45E009020)图幅为例[J]. 地球学报, 2019(5): 651-660. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201905002.htm Beltrão J F, Silva J B C, Costa J C. Robust polynomial fitting method for regional gravity estimation[J]. Geophysics, 1991, 56(1): 80-89. doi: 10.1190/1.1442960
Mickus K L, Aiken C L V, Kennedy W D. Regional-residual gravity anomaly separation using the minimum-curvature technique[J]. Geophysics, 1991, 56(2): 279-283. doi: 10.1190/1.1443041
Mallick K, Sharma K K. A finite element method for computation of the regional gravity anomaly[J]. Geophysics, 1999, 64(2): 461-469. doi: 10.1190/1.1444551
侯遵泽, 杨文采. 中国重力异常的小波变换与多尺度分析[J]. 地球物理学报, 1997, 40(1): 85-95. doi: 10.3321/j.issn:0001-5733.1997.01.010 刘天佑, 吴招才, 詹应林, 等. 磁异常小波多尺度分解及危机矿山的深部找矿: 以大冶铁矿为例[J]. 地球科学——中国地质大学学报, 2007, (1): 135-140. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200701020.htm Xu Y, Hao T, Li Z, et al. Regional gravity anomaly separation using wavelet transform and spectrum analysis[J]. Journal of Geophysics & Engineering, 2009, (3): 279-287.
杨文采, 施志群, 侯遵泽, 等. 离散小波变换与重力异常多重分解[J]. 地球物理学报, 2001, 44(4): 534-541. doi: 10.3321/j.issn:0001-5733.2001.04.012 牟力, 陈召曦. 重力资料多尺度分析最优小波基的选择[J]. 物探与化探, 2015, 39(5): 1013-1019. https://www.cnki.com.cn/Article/CJFDTOTAL-WTYH201505022.htm Daubechies I. Orthonormal base of compactly supported wavelets[J]. Communications on Pure and Applied Mathematics, 1998, XLI: 909-996.
Mallat S, Hwang W L. Singularity detection and processing with wavelets[J]. IEEE Transactions on Information Theory, 2002, 38(2): 617-643.
Spector A, Grant F S. Statistical models for interpreting aeromagnetic data[J]. Geophysics, 1970, 35(2): 293-302. doi: 10.1190/1.1440092
Cook R D. Detection of influential observation in linear regression[J]. Technometrics, 1977, 42(1): 65-68.
Fan D, Li S, Meng S, et al. Predicting submarine topography by linear regression analysis[J]. Journal of Chinese Inertial Technology, 2018, 26(1): 2018.
汪冬华. 多元统计分析与SPSS应用[M]. 上海: 华东理工大学出版社, 2010. 王谦身, 安玉林, 张赤军, 等. 重力学[M]. 北京: 地震出版社, 2003. Jiang R, Cao K, Zeng J, et al. Late Cenozoic tectonic evolution of the southern segment of the Tan-Lu fault zone, Eastern China[J]. Journal of Asian Earth Sciences, 2019, 182: 103932. doi: 10.1016/j.jseaes.2019.103932
-
期刊类型引用(8)
1. 孙旭鹏，李超楠. 瞬变电磁法在煤矿老窑采空区积水探测中的应用. 煤炭技术. 2024(03): 208-212 . 
百度学术
2. 胡斌. 瞬变电磁法在煤矿超前地质勘查中的应用. 能源与环保. 2024(06): 70-75 . 百度学术
3. WEI Laonao，LIU Yunhe，ZHANG Bo. UAV-based transient electromagnetic 3D forward modeling and inversion and analysis of exploration capability. Global Geology. 2024(03): 154-166 . 必应学术
4. 都兵，张旭晴，孙涛，陈峰，王琳炜. 基于时序InSAR的吉林省松原市宁江区地下水活动研究. 世界地质. 2024(03): 424-433 . 百度学术
5. 翟皓，李桐林，康鑫泽，吴宇豪. 半航空瞬变电磁数据处理与解释在贺斯格乌拉露天煤矿水灾害勘探中的应用. 世界地质. 2024(03): 434-443 . 百度学术
6. 秦莉旻. 物探技术在地下水资源污染防治中的应用研究. 环境科学与管理. 2023(03): 72-76 . 百度学术
7. 范莹琳，潘树仁，杜松，李萌，赵岳，张玉峰，丁晏，宋思彤，车巧慧，王锋利. 基于半航空瞬变电磁法识别复杂地形废弃煤矿富水空间的应用研究. 煤炭科学技术. 2023(12): 79-89 . 百度学术
8. 林康利. 基于航空物探的铁路地质选线研究. 工程地球物理学报. 2022(05): 603-609 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 2612
- HTML全文浏览量: 253
- PDF下载量: 1585
- 被引次数: 9
