龙慧，谢兴隆，李凤哲，任政委，王春辉，郭淑君
(中国地质调查局 水文地质环境地质调查中心，河北 保定 071051)
雄安新区是千年大计、国家大事，建设地上、地下两个雄安是新区规划建设的重要目标，地下雄安建设必须建立在对地下地质结构具有充分认识的基础上。雄安新区位于冀中平原中部，石油与国土部门在此开展过大量的油气勘探、地热勘探与构造解释研究[1-6]，积累了丰富的基础地质与地球物理资料。但是，目前的研究主要集中在地表和深部，针对工程建设层(200 m以浅)研究相对较少。自雄安新区设立以来，为满足城市地下空间规划、资源利用开发、地质灾害评价和可持续发展等方面的需求，雄安新区起步区浅层地质结构和隐伏断裂特征的探测和研究这项工作亟待解决。
地球物理探测技术在地下结构探测方面具有不破坏地表、探测效率高、成本低廉等明显优势，槽探、地震勘探、电磁法和探地雷达等技术在浅层地质结构探测中发挥了重要作用，从应用效果来看，这些方法的探测深度、分辨率、适用性以及解决地质问题的能力各有差异。槽探可称为最经典的浅层结构探测方法，但受制于观测区域小、经济成本昂贵、对地表环境破坏性大等，难以实现重复性探测[7-8]。对于电磁法探测技术，由于城市建筑和电场环境的影响，方法局限性较强[9-10]；
相对而言，高密度电阻率测深受带电输电线路、杆塔等影响较小，在工程勘察中的应用越来越广泛。浅层地震无疑是划分地层结构最好的方法之一，能够对第四系底界面及内部结构进行地球物理意义上的精细划分，识别隐伏活动断裂及其空间展布特征[11-12]，但城市噪声污染严重，对地表环境影响和后期数据处理要求较高，在城市浅层地下空间精细探测上，还需要配合其他方法进行综合分析。探地雷达是采用脉冲电磁波发射和接收，数据采集和处理类似于地震零偏移距自激自收方式，具备分辨率和效率“双高”的特征，但由于脉冲电磁波在地下介质中的强衰减特征，其探测深度极其有限[12]。综上所述，在雄安新区大规模建设阶段开展地球物理勘查工作，面临城市强干扰背景、探测精度、地表环境等因素的影响，无损、高效、抗干扰能力强的探测技术在提供城市地区地下地质信息方面可以发挥关键作用。
本文以浅层地震勘查为主，高密度电阻率测深为辅在雄安新区起步区开展综合地球物理探测，综合钻孔、录井和测井资料，对研究区浅层地质结构进行精细探测和描述，同时构建200 m以浅三维地质结构模型，支撑和服务“透明雄安”。
雄安新区辖区河北雄县、容城、安新3县及周边部分区域，构造位置属于渤海湾盆地冀中坳陷中部[13]，横跨徐水凹陷、容城凸起、牛背斜坡、牛驼镇凸起、霸州凹陷和文安斜坡等次级构造单元，南为保定凹陷、高阳低凸起、饶阳凹陷，北为廊固凹陷。冀中坳陷在中生代燕山运动挤压和新生代太平洋板块俯冲弧后拉张伸展的构造作用下形成的一系列正断层[14]，研究区内规模较大、对区域构造起控制作用的隐伏基岩断裂主要为保定—石家庄断裂、容西断裂、容东断裂、牛东断裂和徐水断裂[15-17]。新生代断裂活动主要集中在古近系沙河街组四段(Es4)和新近系馆陶组(N2g)、明化镇组(N2m)沉积时期，地表未见构造活动迹象。
沉积地层序列由下而上依次为太古宇变质岩，以角闪岩和片麻岩为主；
元古宇青白口系(Qb)、蓟县系(Jx)和长城系(Ch)，青白口系以灰岩为主，蓟县系和长城系以白云岩为主，还包括玄武岩和页岩；
下古生界寒武系(∈)由泥岩、灰岩和白云岩组成；
下古生界奥陶系(O)，上古生界石炭系(C)、二叠系(P)，中生界侏罗系(J)、白垩系(K)分布零星，在起步区缺失；
新生界古近系(E)和新近系(N)分别以湖相和河流相沉积为主，岩性为砂岩和泥岩，第四系(Q)以来以冲洪积、冲积和冲湖积物为主，岩性包括砾石、粗砂、中砂、细砂、粉砂、粉土、粉质黏土、黏土[18-20]。
前人研究表明，雄安新区浅层沉积物整体上呈现低密度、中低磁化率、低电阻率、低波速的物性响应特征，第四系内部地层物性特征差异较小[18,20]。本文通过浅层地震勘查和高密度电阻率测深两种地球物理方法探测雄安新区起步区工程建设层地质结构。
2.1 浅层地震勘查
浅层纵波反射地震是城市地下空间探测常用的勘探技术，对于浅地表50～300 m深度范围内的砂泥岩互层的变化特征和断层结构特征揭示效果较好，通过测井约束和振幅分析可获取高分辨率的砂泥岩的地层岩性界面[21]，还可获取S波速度和弹性模量等参数信息[12]。但是反射波地震数据采集和处理通常会存在“盲区”，使浅表(0～50 m)不能获得可靠的成像信息。
为获取高质量反射波信号，数据采集前开展了大量现场野外试验，在了解了研究区低速带情况，掌握了波场特征的基础上，确定了观测系统。此次采用美国Geometrics公司Strata Visor Nzxp24型地震仪及其扩展模块Geode，M18-612型可控震源(扫描频率5～250 Hz，最大出力18 t)及其配套电控系统，L40检波器，扫描方式线性升频、扫描频率20～160 Hz、扫描长度8 s、驱动幅度70%、叠加1～3次的激发参数，120道、道间距5 m、炮检距20 m、采样间隔0.5 ms，记录长度2 s、覆盖次数最高达40次的接收参数。数据处理采用自主开发模块、GeoEast V3.2、CGG、vista处理系统等软件，主要处理过程包括折射静校正、地表一致性处理、频谱分析、弯线共中心点道集编排、动校正、速度分析等。同时，为了获得更准确的地震速度资料，提高解释精度，结合了剖面上钻孔的录井，电阻率和声波时差测井，叠加速度谱等资料进行综合分析，获取研究区的速度模型用于地震资料的时深转换。
图1为雄安新区起步区R2-1测线地震时间剖面成果，200 m以浅第四系地层共有2个明显的地震反射波组(图2)。T0反射面双程反射时间为82～132 ms，地层推断解释为全新统—中更新统底界，埋深约20～80 m；
TQ反射面双程反射时间为250～292 ms，地层推断解释为第四系下更新统(Qp1)底界(即新近系明化镇组(N2m)顶界)，埋深约100～220 m。该剖面分别在距离505 m和5 135 m处200 ms以下下更新统和明化镇组地层中发生扰动，两侧波组特征发生明显变化，解释为F1和F7断层。F1断层倾向NW，倾角19.50°，上断点埋深168 m，断层上下盘地层变化较大，从断层的特征和规模来看，推测为容南断裂，F1-1、F1-3为其分支断层。F7断层倾向NE，倾角10.90°，上断点埋深349 m，断距随着地层深度增加逐渐增大。
图1 研究区R2-1测线地震时间剖面
二维地震探测结果揭示雄安新区起步区地层序列自上而下包括新生界第四系全新统—中更新统(Qh～Qp2)、下更新统(Qp1)、新近系明化镇组(N2m)、馆陶组(N2g)、古近系东营组(Ed)、沙河街组(Es)、孔店组(Ek)、古生界寒武系(∈)、中—新元古界青白口系(Qb)顶底板等15套地层(图2)，其中，明化镇组、馆陶组、东营组和沙河街组地层在容城凸起区域均有不同程度的缺失，孔店组和古生界地层仅在容城凸起西侧，即南张镇—小李镇—三台镇西北一带发育。
图2 研究区主要地层序列分布
依据同相轴发生错段，上下盘同相轴连续性差异明显等特征，推断存在F1～F9等断裂(图3)。区域性断裂主要有容南断裂(F1)，控制着容城凸起和保定凹陷的发育；
容东断裂(F2)，为容城凸起和牛北斜坡的构造分界线；
徐水断裂(F8)，是分割徐水凹陷和保定凹陷的边界断裂。F1、F2、F4、F5、F8、F9分别存在分支断层，F1、F4、F5、F8、F9与其分支断层均呈“Y”字形组合。断裂构造隐伏于地表之下，走向以NE向为主，控制了第四系、新近系、古近系、中生界和中—新元古界的空间分布；
除F1及其次级断裂F1-3对第四系下更新统地层有控制作用外，其他断裂早更新世以来活动不明显。
图3 研究区断裂构造分布
2.2 高密度电阻率测深
高密度电阻率测深主要用于探查浅部岩土层电性结构和隐伏断裂发育特征，被广泛应用于100 m以浅水文地质和工程地质领域，具有自动化程度高、数据信息采集量大、反映地下电性结构直观等优势，是一种现代化勘查技术方法[22-23]，最重要的是补充了地震勘查的浅部盲区。为获取高质量电阻率数据，同时满足一定探测深度，此次使用稳定性较好的骄鹏E60DN高密度电阻率测深仪，5～10 m采样点距，温纳α装置模式测量，参数设置为最大电压100～300 V、最大电流2 A、单点供电1 s、重复观测2次。数据处理采用RES2DINV软件进行反演，反演方法为最小二乘法。
从高密度电阻率测深结果来看，雄安新区起步区200 m以浅地层视电阻率分层特征明显，纵向上表现为“低—高—低”3个明显电性结构特征(图4)。结合钻孔信息可解释为：0～15 m，视电阻率值较低，范围在5～15 Ω·m之间，主要以回填土、粉土、黏土为主；
15～100 m，视电阻率值15～70 Ω·m，整体表现为高阻特征，岩性以粉土、粉质黏土、中砂、细砂为主，局部成团状分布，形成砂层透镜体；
100～200 m，视电阻率值10～50 Ω·m，为粉质黏土与中砂、细砂互层。
图4 雄安新区高密度电阻率测深反演结果
高密度电阻率测深不仅提供了起步区浅层地—电结构特征，结合工程地质孔和地下水矿化度资料，划分了含水层有利区段为15～30 m和40～80 m，还定量表示了地下水矿化度(TDS)与视电阻率之间的关系(图5)，为调查区咸淡水分界面的划分提供现实依据。
图5 地下水矿化度(TDS)与视电阻率关系
本次选用Schlumberger开发出品的Petrel三维可视化地质建模软件构建雄安新区起步区200 m以浅三维地质结构模型。该软件能够综合利用地震数据、钻孔数据、测井资料和地质信息等多源数据，实时对断层的平面组合与空间展布规律进行研究，实现地质资料与地震数据的有机结合，从而构建三维构造模型，进而更加直观准确地分析雄安新区浅地表各种复杂的地质现象[24]。此次建模数据包括 21 236 物理点二维地震数据、26眼地热钻孔分层数据和17眼工程地质钻孔测井数据。
3.1 建模方法
平原地区工程建设层地质结构多为第四系层状地层沉积，可采用自动建模的方法[20,25]。模型构建的过程大致归纳为以下三步：首先，收集、整理和筛选钻孔编录，测井，地形，二维地震解释层位和断层等数据，按照软件格式要求将标准化后的数据导入，通过插值和拟合算法由上至下生成层状地质界面，并用地热钻孔分层数据作为约束；
其次，利用边界和层面数据确定矢量化的建模范围，选择不同的插值方法剖分三维可视化网格，同时采用Structural Modeling模块进行整个系统的构造建模；
最后，经过多次拟合调整，最终选用多重网格逼近法进行插值计算，剖分网格243 360个，逐层构建地质实体图。
3.2 三维地质结构模型
雄安新区重点调查区0～200 m三维可视化地质模型如图6，发育新生界第四系全新统—中更新统、下更新统、新近系明化镇组地层，其中明化镇组地层仅在工作区东北部地区发育。地表高程模型采用13 689个RTK实测数据获得，RTK控制网采用Ⅱ级网，数据来源于国家测绘局资料与千寻网络位置。测量的投影方式采用UTM投影，投影的中央子午线为117°E。从图6可以看出，雄安新区重点调查区地表高程范围2.00～17.01 m，平均海拔15.01 m，高差较小，总体呈现北高南低、东高西低的趋势，小李镇—三台镇以西一带高程较两侧低。
图6 研究区0～200 m三维可视化地质模型
重点调查区新生界第四系松散层底界埋深100.60～211.24 m，平均110.64 m，厚度114.64～218.40 m，平均厚度103.76 m(图7)，总体变化趋势为北高南低，小李镇和安新县城北存在局部凸起，三台镇西南和大王镇西北存在局部凹槽。第四系全新统—中更新统底界埋深20.70～76.56 m，平均55.86 m，厚度26.94～84.46 m，平均厚度57.52 m。下更新统地层厚度59.30～147.77 m，平均厚度88.47 m，发育F1、F1-3断层。F1及其次级断裂位于三台镇地区，呈“Y”字形发育。F1为走向由NW92.92至NE55.97，倾向由NE至NW，倾角由92.17°～94.26°的正断层，上断点最小埋深103.44 m；
F1-3为走向NE81.46，倾向SE，倾角7.08°的正断层，上断点埋深119.75 m。
图7 研究区第四系底界埋深
根据雄安新区起步区实际探测效果，对比分析此次二维地震勘查和高密度电阻率测深解决地质问题的能力(表1)，总结物探技术在沉积盆地的综合探测效果。结合在雄安新区开展的其他地球物理方法，建立不同深度地下地质结构探测模式和城市地下空间探测方法技术体系。
表1 雄安新区地球物理方法适用性分析
1)通过在雄安新区核心区开展二维地震和高密度电阻率测深联合探测，将雄安新区200 m以浅工程建设层刻画为4套地层，地层序列自上而下包括新生界第四系全新统—中更新统、下更新统和新近系明化镇组，下更新统及以下地层发育容南断裂和其次级断裂。
2)以地球物理勘查数据为主要建模数据、钻孔和测井为约束条件的多源数据构建雄安新区200 m以浅三维地质模型，能够更加真实地反映地质体的三维空间结构特征。建模过程中，不同地球物理方法之间的数据融合难度较大，仍需开展进一步研究。
3)结合实际探测效果，对地球物理探测方法解决地质问题的能力进行分析和研究，能够指导沉积盆地城市地下空间地球物理探测。
猜你喜欢 第四系浅层高密度 晋西黄土区极端降雨后浅层滑坡调查及影响因素分析水土保持学报(2022年5期)2022-10-10基于三元混合模型的滨海金矿矿井涌水水源分析有色金属（矿山部分）(2022年2期)2022-07-13高密度养殖南美白对虾或者更容易成功当代水产(2022年3期)2022-04-26高密度电法在断裂构造探测中的应用资源信息与工程(2021年5期)2022-01-15渤海湾西第四纪以来差异性沉积特征及其构造指示城市地质(2021年3期)2021-11-06高密度存储服务器可靠性设计与实现①科学家(2021年24期)2021-04-25高密度存储服务器可靠性设计与实现科学家(2021年24期)2021-04-25浅层换填技术在深厚软土路基中的应用建材发展导向(2021年24期)2021-02-12基于无人机测绘-物探-钻探组合技术的某场地第四系地层特征研究科学与生活(2021年27期)2021-01-11江苏丹阳地区第四系物性特征研究赤峰学院学报·自然科学版(2019年9期)2019-09-10