北方岩溶区断裂带水文地质性质及结构模型

刘元晴, 周乐, 王新峰, 吕琳, 路小慧, 于开宁, 张伟峰. 北方岩溶区断裂带水文地质性质及结构模型[J]. 中国岩溶, 2022, 41(6): 975-985. doi: 10.11932/karst20220609
引用本文: 刘元晴, 周乐, 王新峰, 吕琳, 路小慧, 于开宁, 张伟峰. 北方岩溶区断裂带水文地质性质及结构模型[J]. 中国岩溶, 2022, 41(6): 975-985. doi: 10.11932/karst20220609
LIU Yuanqing, ZHOU Le, WANG Xinfeng, LV Lin, LU Xiaohui, YU Kaining, ZHANG Weifeng. Hydrogeological structure model of the fault zone in the karst area of north China[J]. Carsologica Sinica, 2022, 41(6): 975-985. doi: 10.11932/karst20220609
Citation: LIU Yuanqing, ZHOU Le, WANG Xinfeng, LV Lin, LU Xiaohui, YU Kaining, ZHANG Weifeng. Hydrogeological structure model of the fault zone in the karst area of north China[J]. Carsologica Sinica, 2022, 41(6): 975-985. doi: 10.11932/karst20220609

北方岩溶区断裂带水文地质性质及结构模型

  • 基金项目: 山东省中央引导地方科技发展资金项目(YDZX20203700002937);河北省高校生态环境地质应用技术研发中心开放课题(JSYF-202101);中国地质调查局地质调查项目( DD2022175403)
详细信息
    作者简介: 刘元晴(1988-),高级工程师,主要从事基岩山区水文地质调查工作。E-mail:lyq198896@126.com
    通讯作者: 于开宁(1965-),教授,主要从事水文地质、环境地质方面的研究。E-mail:1211931193@qq.com
  • 中图分类号: P641.1

Hydrogeological structure model of the fault zone in the karst area of north China

More Information
  • 断层带结构和内部流体流动特性是水文地质研究领域的难点问题。石油地质领域,在油气运移与成藏方向已形成较成熟的断层封闭性定量评价技术手段。相比较,断裂的水文地质性质研究尚停留在断裂的力学性质对断层导水、阻水特性的定性评价阶段,尚未详细开展断裂带结构、渗透性各向异性等方面的研究工作。文章梳理总结国外断裂带水文地质性质研究中关于结构组成、断裂带演化、渗透率影响因素等方面的研究成果,引入断裂带渗透率结构模型,并以中国北方岩溶区碳酸盐岩与碎屑岩互层含水岩组为例,构建断裂带水文地质结构模型。断裂带研究尺度和精度不同、断裂带发育部位不同,导致其结构及水文地质性质亦不相同,如何建立起精确、典型的断裂带水文地质结构模型,需要各领域数据共享及多学科融合共同开展研究工作。

  • 氦气是国防军工和高科技产业发展不可或缺的稀有战略性物资之一(李玉宏等,2022),在航天、国防和高端能源系统,半导体和光纤制造等工业领域,医学成像与深潜水等民生领域应用广泛(贾凌霄,2022)。地球上氦气资源极为有限且分布极不均匀。美国地质调查局2021年报告显示,全球氦气总资源量约为484×108 m3,主要分布在美国、卡塔尔、阿尔及利亚、俄罗斯等4国(张宇轩等,2022),其资源量总和占全球总量的87%。自20世纪80年代以来,全球氦气资源长期短缺,中国氦气年需求量近3000万m3,但其产量不足100万m3,长期依赖进口。目前,国内已发现的氦气资源主要集中在中−西部盆地和东部郯庐断裂带两侧含油气盆地。中−西部主要发现于四川盆地、塔里木盆地、柴达木盆地、鄂尔多斯盆地及南部渭河盆地(李玉宏等,2016陶小晚等,2019Wang et al.,2020晁海德等,2022贺政阳等,2022赵安坤等,2022周俊林等,2022;)。此外,东北地区的松辽、海拉尔、辽河、铁法等盆地发现了具有工业价值的氦气资源,晋中的临汾–运城盆地、江西省的温泉中也发现了氦气资源(邹勇军等,2022司庆红等,2023张健等,2023)。

    氦在自然界中有3He和4He两种稳定的同位素(徐永昌,1997) ,其来源主要有3个:大气源、壳源(放射性来源)和幔源,其中大气中的氦含量极少,且提取难度极大,故可忽略不计。中国东部含氦盆地中的氦气来源以壳幔混合为主,而中西部地区的氦气为壳源成因,在20世纪70年代,渭河盆地的天然气井中发现有良好的富氦天然气显示,氦气为壳源成因,氦源岩为盆地内分布的多期富Th、U的花岗岩(李玉宏等,2011)。目前,中国对氦气资源的研究程度较低,资源量认识也不够清晰,仅在四川威远气田开展过工业制氦,因其开发时间较长,天然气资源几乎枯竭(刘凯旋等,2022),寻找新的富He天然气藏刻不容缓。近年来,随着对鄂尔多斯盆地油气勘探开发的进一步深入,在盆地东缘石西区块中的煤层气与砂岩气中均发现有氦气显示(刘超等,2021),多为0.05%~0.15%,为含He天然气;盆地北部东胜气田天然气中氦气含量为0.045%~0.487%,达到含He−富He气田标准,且氦气来源于太古宇—元古宇变质岩−花岗岩系衰变产生(何发岐等,2022),且普遍认为鄂尔多斯盆地氦气为典型的壳源氦特征(戴金星等,2005孙晓等,2021)。盆地边部石西区块与东胜气田含He天然气的发现,表明该盆地在一定的地质条件下会形成富He天然气,为该盆地氦气资源的进一步调查提供依据。但盆地内所开展的天然气中氦气的相关研究较少,对于氦气分布规律认识并不充分。

    盆地断裂体系及基底岩性研究对于认识氦气分布规律具有关键作用。关于鄂尔多斯盆地及周缘的基底结构与断裂特征前人已做了一定的研究(周正,2009李明等,2010何紫娟,2013许文强等,2015包洪平等,2019李冰等,2019;),但上述研究大多是针对整个盆地或是盆地南缘展开的,缺乏对盆地东南部构造特征系统的研究。此外,目前对于鄂尔多斯盆地东南缘基底岩性也鲜有研究(李明等,2012)。区域性重磁资料具有经济、覆盖面积广、横向分辨率高等优势,已广泛应用于盆地构造研究,以重磁技术为主的地球物理方法在氦气资源分布规律研究方面也取得了一定的成果(路利春等,2017张春灌等,2017Feng et al.,2022)。笔者以重磁资料为主,推断鄂尔多斯盆地东南部的断裂构造及潜在的氦源岩,为研究氦气分布规律提供重要的地球物理资料支撑。

    鄂尔多斯盆地位于华北地块西部,是一个古生代稳定沉降、中生代坳陷东移、新生代多断陷的多旋回克拉通盆地(杨华等,2014黄志刚等,2016吴浩等,2017魏柳斌等,2021),盆地演化发展受到华北板块演变的直接影响。盆地基底主要由新太古界(Ar3)—古元古界(Pt1)结晶变质岩系组成,岩石组成极为复杂,大多经历了较强的区域变质作用,属于变质程度较深的区域变质岩系,主要是各种片岩、片麻岩、变粒岩、混合岩、大理岩及花岗片麻岩等(包洪平等,2019),盆地基底形态为东高西低,北高南低,呈不对称状(密文天等,2016)。盆地盖层发育的地层包括中生界的三叠系(T)、侏罗系(J)、下白垩统(K1)和新生界的古近系(E)、新近系(N)及第四系(Q)等(何紫娟,2013汤超,2014),各个地层在不同区域横向与纵向上分布差别比较大,主要地层为三叠系、侏罗系和白垩系。

    研究区位于鄂尔多斯盆地的东南部(图1),处于伊陕斜坡和渭北隆起的过渡部位,其东北区域为晋西挠褶带的一部分,盆地的基底岩性在磁性上存在着较大差异,是由不同岩性的岩石所组成,在其分界处可能存在有大型基底断裂,其方向主要以NE向为主,研究区内主要断裂的方向与盆地内部断裂的方向大致相同。伊陕斜坡东南部构造简单,断裂发育也较少(孙建博等,2018);而渭北隆起的构造相对较为复杂,为中生代以来形成的断块隆起(任战利等,2014),断裂则较为发育,由于鄂尔多斯盆地与秦岭造山带长期相互作用,导致中新生代沉积盖层发生了较为强烈的构造变形,从北到南构造规模由小到大,构造变形由弱到强。鄂尔多斯地块周缘有广泛的岩浆岩出露,且从太古宇到白垩纪均有分布。鄂尔多斯盆地北缘有古元古代2期以及早石炭世末期—二叠纪、早—中侏罗世、早白垩世共5期岩浆活动;南缘有新元古代、古生代、早中生代和晚中生代4期岩浆活动;东缘则以发育中生代以来的岩浆活动为主。此外,盆地内部还发育多个已被证实的隐伏岩浆岩体,包括龙门岩体、吴忠岩体和乌海东岩体等(何登发等,2021)。

    图 1.  研究区位置示意图
    Figure 1.  Location of study area

    综合分析和总结研究区岩石、地层密度特征是进行重力异常处理与解释的前提和重要依据(孟军海等,2021)。笔者通过收集前人在研究区及邻区的物性研究成果,分析总结研究区岩石、地层密度特征(李冰等,2019宁媛丽等,2020)(表1)。研究区内除缺失泥盆系、志留系以外,前寒武系至第四系均有不同程度的出露,区内多发育三叠系、侏罗系、白垩系,只有在研究区南部出露少量二叠系—石炭系、奥陶系—寒武系。区内地层密度随时代变老而逐渐增大,新生界岩性以红色、黄色的黏土、黄土为主,表现为低密度特征,密度值为1.87~2.38 g/cm3;中生界岩性主要以砂岩、泥岩、页岩为主,表现出中等密度特征,密度值为2.40~2.55 g/cm3;古生界以白云岩、石英岩、灰岩为主,表现出中高密度特征,密度值为2.60~2.70 g/cm3;前寒武系密度值大于2.70 g/cm3。盆地内部暂未发现岩浆岩露头,盆地周缘的闪长岩、花岗岩、玄武岩、二长岩等侵入岩的密度值为2.56~2.69 g/cm3,表现为中高密度特征。

    表 1.  鄂尔多斯盆地岩石物性特征表(据李冰等,2019宁媛丽等,2020
    Table 1.  The petrophysical property in Ordos basin
    地层/岩体 岩石密度ρ(g/cm3磁化率(10−5 SI)
    密度值 平均值 密度值 平均值
    新生界第四系1.63~2.101.8736~7956
    古近系
    新近系
    2.30~2.392.385~139.6
    中生界白垩系2.4010.5~58.918.6
    侏罗系2.41~2.442.436.5~17.310.5
    三叠系2.558.8~18.615.3
    古生界二叠系2.602.5~22.016.3
    石炭系3.2~20.89.6
    奥陶系2.705.6~14.68.8
    寒武系10.5~16.812.7
    元古界2.705.4~9.27.8
    太古界变粒岩2.761300~106005600
    片麻岩8.6~1400780
    混合岩10.5~16.813.6
    大理岩11.6~24518.4
    侵入岩闪长岩2.5~3.340~95455332
    花岗闪长岩2.562.561374.18
    花岗岩2.580~60483206
    玄武岩2.8~3.32.6910~27131250
    二长石2.47~2.712.59771.16
     | Show Table
    DownLoad: CSV

    鄂尔多斯盆地沉积岩几乎无磁性,磁异常主要来自结晶基底的变质岩系和侵入岩。统计数据显示(表1),研究区内地层分为弱–无磁性层、磁性层。新生界—元古界的磁化率为7.8×10−5~56×10−5 SI,呈弱磁性或无磁性特征。太古界片麻岩、变粒岩等变质岩的磁化率值一般为780×10−5~5600×10−5 SI,分布比较广泛,为区域磁性层。除此之外,盆地周缘分布有多种侵入岩,如闪长岩、花岗岩、玄武岩和二长石等,其磁化率一般为700×10−5~5400×10−5 SI,是引起盆地周缘地区高磁异常的主要原因。

    研究区布格重力异常(图2)资料的比例尺为1∶50万,异常整体呈现北西低、南东高的特征,与研究区构造特征相吻合。研究区东南部大荔及其北东区域,重力异常呈NEE向高低相间的条带,为渭河盆地的一部分,其重力异常梯级带可能是渭河盆地北缘断裂的反映,反映了该盆地呈隆坳相间的特征。研究区中南部铜川–合阳–黄龙一带,重力异常主要呈NE向高值条带,该区为渭北隆起区,其为中生代以来形成的断块隆起,寒武系等老地层出露,中生界等地层较薄,是引起高重力异常的主要原因。研究区北西部安塞至延安一带,呈区域性重力低,局部异常多呈NE向、近NS向展布。该区为伊陕斜坡的一部分,地层以较厚的中生界为主,是重力异常呈现平缓低重力特征的主要原因。研究区中部偏北西的富县–宜川–大宁一带,以平缓的中低重力异常为主,富县–延川存在局部NNE向中高重力异常,可能反映了伊陕斜坡内基底存在局部的凸起,在低重异常与高重力异常过渡区域可能存在有断裂。

    图 2.  布格重力异常图
    Figure 2.  Bouguer gravity anomaly

    研究区化极磁力异常资料(图3)由1∶5万~1∶20万不等比例尺的磁测资料拼合而成,整体呈现为“两高夹一低”的形态。研究区东南部合阳一带表现为高磁异常,整体磁力异常值为200~600 nT,呈NE向带状展布,推测该高磁异常主要由基底强磁性变质岩如片麻岩、变粒岩引起,在合阳东南侧存在着NE向的磁力异常梯级带,可能为基底断裂的反映。盆地中部以铜川–黄龙–宜川–大宁和富县–延川为界,中部区域表现为整体的磁力低,推测主要由基底弱磁性变质岩如大理岩、混合岩等引起。在黄龙、宜川等地存在局部高磁异常,可能由盆地基底之中强磁性变质岩引起的,延安至延川西北区域表现为中高磁力异常,磁力异常值为0~200 nT。鄂尔多斯盆地的沉积盖层磁性一般较弱,只有结晶基底中片麻岩、变粒岩等属于强磁性岩体,因此大范围的高磁异常是强磁性结晶基底的反映。基底岩性的不同造成了区域性磁力异常的不同,化极磁力异常规模较大的梯级带可能为基底大断裂的反映,其也为基底岩性的分界线。

    图 3.  化极磁力异常图
    Figure 3.  Magnetic anomaly reduced to the pole

    断裂活动破坏原有地质体的连续性,造成物性(密度、磁性)上的横向差异,使得断裂两侧呈现明显的重力和磁力异常。归一化总水平导数垂向导数(NVDR–THDR)边缘识别技术(Wang et al. ,2009)是一种有效的识别断裂的方法,通过极大值连线位置或极大值错断位置确定断裂构造特征线,已广泛应用于断裂识别之中(王万银等,2014纪晓琳等,2019马涛等,2020王学发等,2020)。研究区存在3个级别的断裂,其中一级断裂为盆地基底岩性的区域性分界线,二级断裂主要为基底断裂,控制了基底的隆坳格局,三级断裂则为沉积盖层内部断裂。由表1可知,鄂尔多斯盆地基底岩性在磁性上存在较大差异,密度差异很小,因此以化极磁力异常为主进行一级断裂的划分,重磁异常相结合对二级断裂进行划分,三级断裂则利用重力异常进行划分。

    为了进一步说明断裂识别方法的可靠性,在研究区的中南部沿一条地震剖面(魏国齐等,2019)A–A’(位置如图3)提取化极磁力异常、化极磁力异常NVDR–THDR及化极磁力垂向一阶导数进行对比分析,结果如图4所示。结果表明,化极磁力异常在此剖面上自NW向SE呈减小趋势,且在断裂位置处的明显减小(图4b);化极磁力异常NVDR–THDR剖面的极大值处与地震剖面断裂相叠合(图4c),这与该技术的断裂识别标志相吻合;在化极磁力异常垂向一阶导数剖面中,断裂与高低异常值的分界处及0值线处相叠合(图4d)是断裂在化极磁力异常垂向一阶导数中较好的反映。综上所述,利用该方法结合地质资料可以有效识别出研究区的断裂构造特征。

    图 4.  A–A’剖面与磁力异常对比分析图
    a. 地震资料断裂剖面特征(魏国齐等,2019);b. 化极磁力异常剖面;c. 化极磁力异常NVDR–THDR剖面;d. 化极磁力异常垂向一阶导数剖面
    Figure 4.  Comparative analysis of A–A’ section and magnetic anomalies

    笔者以重磁异常NVDR–THDR为主,结合其他异常,最终在研究区推断出了一级断裂3条、二级断裂8条以及三级断裂51条(图5a)。断裂的主要走向为NE向和近NEE向,其次是NW向和近EW向。

    图 5.  推断断裂分布与磁力异常图
    a. 推断断裂平面分布图;b. 断裂与化极磁力异常图;c. 断裂与化极磁力异常NVDR–THDR图;d. 断裂与化极磁力异常垂向一阶导数图
    Figure 5.  The distribution of inferred faults and magnetic anomalies

    研究区内共推断一级断裂3条,断裂长度均在200 km以上,该类断裂主要为基底不同岩性的分界线。断裂两侧重磁异常特征明显,表现为大型线性异常带、梯级带或不同重磁异常的分界线,异常连续性好且延伸较长。

    (1)F1-1(延安–延川断裂):该断裂位于延安–延川一带,断裂总长为218.8 km,呈NE走向,断裂两侧多出露新近系—第四系、侏罗系,少量出露三叠系、白垩系。该断裂两侧化极磁力异常特征明显不同(图5b),断裂北西侧为区域性高磁异常,而其南东侧以大范围平缓低磁异常为主,表明该断裂可能为不同岩性基底的分界线,可能反映了鄂尔多斯盆地基底早期的拼合特征。在化极磁力异常NVDR–THDR图上(图5c)表现为一系列NE、NEE向的极大值条带,可能反映了断裂附近构造活动较为强烈。在化极磁力异常垂向一阶导数图上(图5d),断裂处于高低异常值的分界线上,在断裂的上半部分与下半部分表现尤为明显。该断裂在重力异常图上也有较明显的反映,在布格重力异常(图6a)上,断裂两侧异常特征不同,北西侧以低重力异常为主,而南东侧以中高重力异常为主。在布格重力异常NVDR–THDR图上(图6b),由几个NE向极大值条带组成,异常特征较为明显。由此可见,该断裂也控制了基底之上沉积层的展布,使得断裂两侧呈现出了不同的重力异常特征。

    图 6.  推断断裂分布与重力异常图
    a. 断裂与布格重力异常图;b. 断裂与布格重力异常NVDR–THDR图
    Figure 6.  The distribution of inferred faults and gravity anomalies

    (2)F1-2(渭北隆起北缘断裂):该断裂为渭北隆起北界,长约为342.7 km,呈NNE向展布。断裂两侧主要出露三叠纪、侏罗纪,并见少量新近系。断裂两侧表现为截然不同的化极磁力异常特征(图5b),其北西侧为大范围平缓低磁异常为主,而断裂的南东侧则是以NE向的大型高磁异常带为主,表明断裂可能为不同岩性基底的分界线,其两侧表现为不同的区域构造格局。在化极磁力异常NVDR–THDR图上(图5c)表现为一系列NE、NEE向的极大值条带,该断裂为逆断层,规模较大,构造活动性强,是渭北隆起与陕北斜坡的边界断裂(许文强等,2015)。在化极磁力异常垂向一阶导数图中(图5d),断裂处于高低异常值的分界处与0值线上,有较好的断裂构造显示。该断重力异常图上也有较为明显的反映,在布格重力异常(图6a)上,断裂的南段处于高重力异常区上,断裂北段两侧则表现为不同的重力异常特征。在布格重力异常NVDR–THDR图上(图6b),断裂的南段没有明显的异常特征,断裂北段由几个NNE向的极大值条带组成。

    (3)F1-3(渭河盆地北缘断裂):断裂位于渭北隆起的南缘,为渭北隆起与渭河盆的边界断裂,总长约为224 km,呈NE向分布。断裂两侧主要出露的地层为侏罗系和二叠系,整体处于高磁异常带上(图5b),其两侧均表现为区域性高磁异常特征,大荔及其研究区东南部则为局部低磁异常,该断裂为北升南降的高角度正断层,构造活动复杂,规模大,活动期长,近期仍有活动(许文强等,2015)。在化极磁力异常NVDR–THDR图中(图5c),断裂处于极大值连线上,与极大值对应较好。在化极磁力异常垂向一阶导数图中没有明显的表现。该断裂在重力异常图上也有一定的反映,在布格重力异常(图6a)上,断裂的北段处于高重力异常上,中部则处于局部低重力异常区,南段处于重力异常梯级带上,在布格重力异常NVDR–THDR图中(图6b),由几个NE向极大值条带组成,异常特征较为明显,可能反映了断裂附近构造活动较为强烈,小规模断裂较为发育。

    研究区内共推断二级断裂8条,断裂长度为100~275 km,该类断裂重磁异常特征较为明显,具有一定规模的线性异常带、梯级带、串珠异常带等,在化极磁力异常NVDR–THDR图上表现为极大值连续性较好。

    (1)NW向断裂系(F2-1~F2-4):在整个研究区内部,除了发育北东向的大型盆地基底断裂,还发育了一些NW、NWW向的二级断裂,长度为160~230 km。处于具有一定规模的化极磁力异常线性带、梯级带上(图5b),断裂系的北西部分处于高磁异常区,中部则处于区域性低磁异常区,南东部分处于高磁异常带上,表明该断裂的形成时间明显晚于NE向构造,对后者起到切割和错断作用,使得研究区中部,整体分段向左滑动。在化极磁力异常NVDR–THDR图上(图5c),断裂系处于极大值条带上或极大值的错断处。在化极磁力异常垂向一阶导数图及重力异常图上则没有明显的反映。

    (2)NE向断裂系(F2-5~F2-8):研究区发育了一些NE向的二级断裂,断裂长度为100~275 km,大致呈等间距排列,断裂方向与一级断裂方向大致相同,形成时间早于NW向断裂系,被切割错断,导致NW向断裂系向左滑动。在化极磁力异常图上(图5b),F2-5、F2-8处于磁力异常高值区,F2-6则处于磁力异常梯级带上,F2-7处于局部高磁异常上。在化极磁力异常NVDR–THDR图上(图5c)也有一定的反映,断裂处于极大值或极大值错断处。在化极磁力异常垂向一阶导数图中(图5d),F2-6处于0值线处,是断裂构造的标志。在布格重力异常(图6a)上,F2-5、F2-7、F2-8处于重力异常梯级带上,F2-6则处于重力异常高值带上;在布格重力异常NVDR–THDR图中(图6b),极大值与断裂对应较好,表明该断裂系在一定程度上控制着研究区沉积层的分布,使得断裂两侧表现为不同的重力异常特征。

    除一级、二级断裂外,盆地发育了大量的三级断裂,共计51条,其断裂方向以NE、NW向为主,分布于整个研究区。三级断裂主要为沉积层内部断裂,与磁力异常的关系不明显,主要利用布格重力异常与布格重力异常NVDR–THDR进行划分,表现为布格重力异常局部错断位置、重力异常梯级带或NVDR–THDR的极大值或局部错断位置(图6a图6b)。该类断裂规模小,数量多,可能为氦气运移的通道。

    国内研究较多的含He盆地(如渭河盆地)的勘探实践表明,地球物理是研究盆地氦气资源分布规律的有效方法,其中磁力资料能在很大程度上反映氦源岩的分布特征。例如,渭河盆地内部高氦气含量井在空间分布上往往与高磁异常有关(李玉宏等,2011张春灌等,2017)。鄂尔多斯盆地氦气相关研究成果表明盆地氦气为壳源成因,其成因为壳内岩体放射性衰变。花岗岩被认为是渭河盆地和四川威远气藏的主要氦源,花岗岩中U、Th元素放射性衰变形成氦气;变质岩中U、Th元素含量也较高(陈道公等,2004),也可作为氦源岩,其在东胜气田得到了证实(何发岐等,2022)。因此,如何识别盆地内部强磁性变质岩是氦气资源远景调查的关键之一。笔者根据化极磁力异常垂向一阶导数、剩余化极磁力异常与划分的一级、二级断裂,并结合地质资料进行分析,将研究区划分为西北部、中部及东南部3个大区,并识别出有效的氦源岩。

    图3显示,3个大区的化极磁力异常存在明显的差异,引起该异常变化的原因可能有2种:①磁性结晶基底的起伏所引起的差异。②不同岩性的磁性结晶基底引起磁力异常变化。为了明确引起异常变化的原因,在研究区沿一条近SN向的剖面B–B’(徐兴雨,2020)(图3),进行磁力剖面拟合,通过修改模型参数,当拟合的正演曲线与实测曲线差异最小时,认为拟合效果达到了最优。

    图7为B–B’剖面拟合结果,剖面位于研究区的西部,长约为184 km。不同磁性基底所构成的基底模型共分为9个,拟合时假设沉积层无磁性。由拟合结果可知(图7b),结晶基底由北到南逐渐升高,其两侧的磁化强度远高于中部,拟合曲线与实测曲线差异较小,故比较符合地质构造特征。同一磁性的结晶基底的磁化率磁化强度不变,磁化强度设置为0.04 A/m,模型不变,正演出的曲线与实测曲线相差较大(图7c),不符合该地质特征。综上所述,结晶基底磁性差异是引起区域性磁力异常不同的主要原因。

    图 7.  磁力剖面拟合图(图中数值为磁化强度)
    a. 南北向地震剖面(徐兴雨,2020);b. 基底不同磁性的拟合结果;c. 基底同一磁性的拟合结果
    Figure 7.  Magnetic profile fitting(The value in the figure is the magnetization intensity)

    西北与东南区域分别以F1-1、F1-2断裂为边界,西北部以大范围的中高磁力异常为主,东南区域也分布有大型NE向高磁异常带(图8a)。盆地基底岩石中大量的片麻岩、变粒岩等强磁性岩体,是引起区域性高磁异常的主要原因,这些区域内可能存在大规模的强磁性变质岩等氦源岩,具有较好的气源条件。研究区中部F1-1、F1-2断裂之间的区域,为区域性低磁异常区(图8a),反映了基底岩性整体以混合岩、大理岩等弱磁性岩体为主,局部等轴状、椭圆状高磁异常可能为局部分布的片麻岩等强磁性变质岩的反映(如G1~G21)(图8b)。剩余化极磁力异常(图8b)中,局部高磁异常G1~G21主要集中在富县–宜川–黄龙一带,在西南缘及东北缘也有分布,其可能为基底强磁性变质岩的反映,为潜在的氦源岩。

    图 8.  化极磁力异常垂向一阶导数图与剩余化极磁力异常图
    a. 化极磁力异常垂向一阶导数图;b. 剩余化极磁力异常图
    Figure 8.  First vertical derivative of RTP magnetic anomalies and residual polarized magnetic anomalies

    氦源岩与断裂分布也存在一定的联系,如G1~G3岩体分布于研究区的东北缘,岩体面积为60~140 km2,处于F1-2、F2-1、F2-6断裂的交汇区域;G4分布F2-6断裂的右侧,岩体面积为310 km2,岩体规模较大,富县–宜川–黄龙一带的G5~G16岩体亦处于断裂交汇区域,大面积分布的氦源岩使得该地区有着较好的氦气资源潜力。

    (1)以NVDR–THDR位场边缘识别技术为主对该地区的断裂进行划分,推断出一级断裂3条、二级断裂8条和三级断裂51条,其中一级断裂为盆地基底岩性的区域性分界线,呈NE向分布,横穿整个研究区,可能反映鄂尔多斯盆地基底早期的拼合特征。二级断裂为基底断裂,呈NE向、NW向展布,其中NW向断裂形成时间晚于NE向断裂,控制基底的隆坳格局。三级断裂多为沉积层内部断裂,规模小,数量多,可能为氦气运移的通道。

    (2)利用化极磁力异常垂向一阶导数、剩余化极异常并结合区域地质资料对氦源岩的分布进行研究,研究区西北、东南部大面积的高磁异常反映该区基底岩性可能以大规模的强磁性变质岩为主,具有较为充足的氦源条件;中部延川–富县–黄龙–宜川一带基底岩性整体以弱磁性大理岩、混合岩等为主,局部高磁异常为强磁性变质岩的反映,其也为潜在的氦源岩。

    (3)研究区断裂为氦气运移提供了通道,氦源岩与断裂的关系表明,大规模的氦源岩多分布于断裂附近及断裂的交汇部位,形成了研究区氦气资源分布及富集的地质–地球物理背景。

  • 图 1  断裂带结构组成图[28]

    Figure 1. 

    图 2  不同尺度断裂带演化对比模型图[13]

    Figure 2. 

    图 3  各领域采用的断裂带水文地质研究方法对比[11]

    Figure 3. 

    图 4  断裂带渗透率模型[28]

    Figure 4. 

    图 5  碳酸盐岩断裂带水文地质结构模型及应用

    Figure 5. 

    图 6  不同研究精度与尺度条件下断裂带及其渗透率结构模型

    Figure 6. 

    图 7  距离断层尖端的不同位置处断裂带结构模型[47]

    Figure 7. 

    表 1  断裂带组成样式及渗透率结构[28]

    Table 1.  Fault zone architectural styles and permeability structures

    渗透率结构断裂带结构断层核破碎带实例应用的流动模型
    局部导水沿单一曲面或长的离散
    平面发生局部滑动
    缺失或发育较窄缺失或发育较少Shawangunk Mountains
    断层
    具有平行外壁的离散导管
    分散导水沿分布表面和裂缝
    分布发生滑移
    缺乏或发育较窄,狭窄、
    离散和不连续带
    发育的离散滑动面和
    相关的断裂网络
    Hill断裂带等效多孔介质
    局部阻水破碎带内局部滑动发育较好的断层
    核碎裂岩
    缺失或发育较少San Gabriel
    碎裂岩带
    高渗透性含水层(原岩)中的
    弱透水层(断层核)
    复合导水—
    阻水
    变形适应于局部破碎带
    和次生构造分布区
    发育较好的断层
    核碎裂岩
    发育的离散滑动面和
    相关的断裂网络
    Stillwater断裂带夹在两个含水层之间的
    弱透水层(断层核)
    下载: 导出CSV
  • [1]

    Wibberley C A J, Shimamoto T. Internal structure and permeability of major strike-slip fault zones: the Median Tectonic Line in Mie Prefecture, Southwest Japan[J]. Journal of Structural Geology, 2003, 25(1):59-78. doi: 10.1016/S0191-8141(02)00014-7

    [2]

    Singhal B B S, Gupta R P. Applied Hydrogeology of Fractured Rocks (second edition)[M]. Springer, 2010: 1-408.

    [3]

    Childs C, Watterson J, Walsh J J. A model for the structure and development of fault zones[J]. Journal of the Geological Society, 1996, 153:337-340. doi: 10.1144/gsjgs.153.3.0337

    [4]

    Tondi E, Cilona A, Agosta F, Aydin A, Rustichelli A, Renda P, Giunta G. Growth processes, dimensional parameters and scaling relationships of two conjugate sets of compactive shear bands in porous carbonate grainstones, Favignana Island, Italy[J]. Journal of Structural Geology, 2012, 37:53-64. doi: 10.1016/j.jsg.2012.02.003

    [5]

    Faulkner D R, Jackson C A L, Lunn R J, Schlische R W, Wibberley C A J, Withjack M O. A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones[J]. Journal of Structural Geology, 2010, 32:1557-1575. doi: 10.1016/j.jsg.2010.06.009

    [6]

    Aydin A. Fractures, faults, and hydrocarbon entrapment, migration and flow[J]. Marine & Petroleum Geology, 2000, 17(7):797-814.

    [7]

    付晓飞, 许鹏, 魏长柱, 吕延防. 张性断裂带内部结构特征及油气运移和保存研究[J]. 地学前缘, 2012, 19(6):200-212.

    FU Xiaofei, XU Peng, WEI Changzhu, LV Yanfang. Internal structure of normal fault zone and hydrocarbon migration[J]. Earth Science Frontiers, 2012, 19(6):200-212.

    [8]

    Douglas M, Clark I D, Raven K, Bottomley D. Groundwater mixing dynamics at a Canadian Shield mine[J]. Journal of Hydrology, 2000, 235(1-2):88-103. doi: 10.1016/S0022-1694(00)00265-1

    [9]

    Dockrill B, Shipton Z K. Structural controls on leakage from a natural CO2 geologic storage site: Central Utah, U. S. A.[J]. Journal of Structural Geology, 2010, 32(11):1768-1782. doi: 10.1016/j.jsg.2010.01.007

    [10]

    Person M, Banerjee A, Hofstra A, Sweetkind D, Gao Y. Hydrologic models of modern and fossil geothermal systems in the Great Basin: Genetic implications for epithermal Au-Ag and Carlin-type gold deposits[J]. Geosphere, 2008, 4(5):888-917. doi: 10.1130/GES00150.1

    [11]

    Bense V F, Gleeson T, Loveless S E, Bour O, Scibek J. Fault zone hydrogeology[J]. Earth-Science Reviews, 2013, 127:171-192. doi: 10.1016/j.earscirev.2013.09.008

    [12]

    宋佳佳, 孙建孟, 王敏, 傅爱兵, 高建申. 断层内部结构研究进展[J]. 地球物理学进展, 2018, 33(5):1956-1966.

    SONG Jiajia, SUN Jianmeng, WANG Min, FU Aibing, GAO Jianshen. Research progress in the internal structure of the fault[J]. Progress in Geophysics, 2018, 33(5):1956-1966.

    [13]

    Wibberley C A J, Kurt W, Imber J. The internal structure of fault zones: Implications for mechanical and fluid-flow properties[M]. Geological Society, 2008: 1-367.

    [14]

    Bense V F, Shipton Z K, Kremer Y, Kampman N. Fault zone hydrogeology: Introduction to the special issue[J]. Geofluids, 2016, 16(4):655-657. doi: 10.1111/gfl.12205

    [15]

    Giwelli A, Delle Piane C, Esteban L, Clennell M B, Dautriat J, Raimon J, Kager S, Kiewiet L. Laboratory observations of fault transmissibility alteration in carbonate rock during direct shearing[J]. Geofluids, 2016, 16(4):658-672. doi: 10.1111/gfl.12183

    [16]

    Scibek J, Gleeson T, Mckenzie J M. The biases and trends in fault zone hydrogeology conceptual models: global compilation and categorical data analysis[J]. Geofluids, 2016(4):782-798.

    [17]

    付广, 李世朝, 杨德相. 断裂输导油气运移形式分布区预测方法及其应用[J]. 沉积学报, 2017, 35(3):592-599.

    FU Guang, LI Shichao, YANG Dexiang. A method forecasting distribution areas of fault transporting oil-gas migration and its application[J]. Acta Sedimentologica Sinica, 2017, 35(3):592-599.

    [18]

    付晓飞, 方德庆, 吕延防, 付广, 孙永河. 从断裂带内部结构出发评价断裂垂向封闭性的方法[J]. 地球科学, 2005, 30(3):328-336.

    FU Xiaofei, FANG Deqing, LV Yanfang, FU Guang, SUN Yonghe. Method of evaluating vertical sealing of faults in terms of the internal structure of fault zones[J]. Earth Science, 2005, 30(3):328-336.

    [19]

    吴智平, 陈伟, 薛雁, 宋国奇, 刘惠民. 断裂带的结构特征及其对油气的输导和封堵性[J]. 地质学报, 2010, 84(4):570-578.

    WU Zhiping, CHEN Wei, XUE Yan, SONG Guoqi, LIU Huimin. Structural characteristics of faulting zone and its ability in transporting and sealing oil and gas[J]. Acta Geologica Sinica, 2010, 84(4):570-578.

    [20]

    潘晓东, 曾洁, 任坤, 焦友军, 彭聪, 兰干江. 贵州毕节岩溶斜坡地带地下水赋存规律与钻探成井模式[J]. 地球学报, 2018, 39(5):606-612.

    PANG Xiaodong, ZENG Jie, REN Kun, JIAO Youjun, PENG Cong, LAN Ganjiang. Groundwater occurrence characteristics and drilling well models in karst slope zone, Bijie, Guizhou Province[J]. Acta Geoscientica Sinica, 2018, 39(5):606-612.

    [21]

    Yamashita T, Tsutsumi A. Involvement of fluids in earthquake ruptures[M]. Springer Japan, 2018: 1-185.

    [22]

    Butler C A, Holdsworth R E, Strachan R A. Evidence for Caledonian sinistral strikeslip motion and associated fault zone weakening, Outer Hebrides Fault Zone, Scotland[J]. Journal of the Geological Society, 1995, 152(5):743-746. doi: 10.1144/gsjgs.152.5.0743

    [23]

    Schulz S E, Evans J P. Spatial variability in microscopic deformation and composition of the Punchbowl fault, Southern California: Implications for mechanisms, fluid–rock interaction, and fault morphology[J]. Tectonophysics, 1998, 295(1-2):223-244. doi: 10.1016/S0040-1951(98)00122-X

    [24]

    Chester F M, Friedman M, Logan J M. Foliated cataclasites[J]. Tectonophysics, 1985, 111(1-2):139-146. doi: 10.1016/0040-1951(85)90071-X

    [25]

    Jefferies S P, Holdsworth R E, Wibberley C A J, Shimamoto T, Spiers C J, Niemeijer A R, Lloyd G E. The nature and importance of phyllonite development in crustal-scale fault cores: An example from the Median Tectonic Line, Japan[J]. Journal of Structural Geology, 2006, 28(2):220-235. doi: 10.1016/j.jsg.2005.10.008

    [26]

    Faulkner D R, Lewis A C, Rutter E H. On the internal structure and mechanics of large strike-slip fault zones: Field observations of the Carboneras fault in Southeastern Spain[J]. Tectonophysics, 2003, 367:235-251. doi: 10.1016/S0040-1951(03)00134-3

    [27]

    Chester F M, Logan J M. Implications for mechanical properties of brittle faults from observations of the Punchbowl fault zone, California[J]. Pure Applied Geophysics, 1986, 124(1-2):79-106. doi: 10.1007/BF00875720

    [28]

    Caine J S, James P E, Craig B F. Fault zone architecture and permeability structure[J]. Geology, 1996, 24(11):1025-1028. doi: 10.1130/0091-7613(1996)024<1025:FZAAPS>2.3.CO;2

    [29]

    Gudmundsson A, Berg S S, Lyslo K B, Skurtveit E. Fracture networks and fluid transport in active fault zones[J]. Journal of Structural Geology, 2001, 23(2):343-353.

    [30]

    Caine J S, Bruhn R L, Craig B F. Internal structure, fault rocks, and inferences regarding deformation, fluid flow, and mineralization in the seismogenic Stillwater normal fault, Dixie Valley, Nevada[J]. Journal of Structural Geology, 2010, 32:1576-1589. doi: 10.1016/j.jsg.2010.03.004

    [31]

    James P E, Craig B F, Jame V G. Permeability of fault-related rocks, and implications for hydraulic structure of fault zones[J]. Journal of Structural Geology, 1997, 19(11):1393-1404. doi: 10.1016/S0191-8141(97)00057-6

    [32]

    罗胜元, 何生, 王浩. 断层内部结构及其对封闭性的影响[J]. 地球科学进展, 2012, 27(2):154-164.

    LUO Shengyuan, HE Sheng, WANG Hao. Review on fault internal structure and the influence on fault sealing ability[J]. Advances in Earth Science, 2012, 27(2):154-164.

    [33]

    Bruhn R L, Parry W T, Yonkee W A, Thompson T. Fracturing and hydrothermal alteration in normal fault zones[J]. Pure& Applied Geophysics, 1994, 142(3):609-644.

    [34]

    Scholz C H. Wear and gouge formation in brittle faulting[J]. Geology, 1987, 15(6):493-495. doi: 10.1130/0091-7613(1987)15<493:WAGFIB>2.0.CO;2

    [35]

    Hull J. Thickness-displacement relationships for deformation zones[J]. Journal of Structural Geology, 1988, 10(4):431-435. doi: 10.1016/0191-8141(88)90020-X

    [36]

    Wibberley C A J, Petit J P, Rives T. Micromechanics of shear rupture and the control of normal stress[J]. Journal of Structural Geology, 2000, 22(4):411-427. doi: 10.1016/S0191-8141(99)00158-3

    [37]

    Chester F M, Evans J P, Biegel R L. Internal structure and weakening mechanisms of the San-Andreas fault[J]. Journal of Geophysical Research Atmospheres, 1993, 98(B1):771-786. doi: 10.1029/92JB01866

    [38]

    Caine J S, Forster C B. Fault zone architecture and fluid flow: Insights from field data and numerical modeling[C]// Haneberg W C, Mozley P S, Moore J C. Faults and Subsurface Fluid Flow in the Shallow Crust. Washington DC, AGU, 1999: 101–127.

    [39]

    Darnault C J G. Overexploitation and contamination of shared groundwater resources [M]. Netherland: Nato Security Through Science, 2011: 203–226.

    [40]

    王焰新. 我国北方岩溶泉域生态修复策略研究: 以晋祠泉为例[J]. 中国岩溶, 2022, 41(3):331-344.

    WANG Yanxin. Study on ecological restoration strategy of karst spring region in North China: Taking Jinci spring as an example[J]. Carsologica Sinica, 2022, 41(3):331-344.

    [41]

    梁永平, 王维泰, 赵春红, 王玮, 唐春雷. 中国北方岩溶水变化特征及其环境问题[J]. 中国岩溶, 2013, 32(1):34-42.

    LIANG Yongping, WANG Weitai, ZHAO Chunhong, WANG Wei, TANG Chunlei. Variations of karst water and environmental problems in North China[J]. Carsologica Sinica, 2013, 32(1):34-42.

    [42]

    梁永平, 申豪勇, 赵春红, 王志恒, 唐春雷, 赵一, 谢浩, 石维芝. 对中国北方岩溶水研究方向的思考与实践[J]. 中国岩溶, 2021, 40(3):363-380.

    LIANG Yongping, SHEN Haoyong, ZHAO Chunhong, WANG Zhiheng, TANG Chunlei, ZHAO Yi, XIE Hao, SHI Weizhi. Thinking and practice on the research direction of karst water in northern China[J]. Carsologica Sinica, 2021, 40(3):363-380.

    [43]

    高旭波, 王万洲, 侯保俊, 高列波, 张建友, 张松涛, 李成城, 姜春芳. 中国北方岩溶地下水污染分析[J]. 中国岩溶, 2020, 39(3):287-298.

    GAO Xubo, WANG Wanzhou, HOU Baojun, GAO Liebo, ZHANG Jianyou, ZHANG Songtao, LI Chengcheng, JIANG Chunfang. Analysis of karst groundwater pollution in Northern China[J]. Carsologica Sinica, 2020, 39(3):287-298.

    [44]

    卢海平, 张发旺, 赵春红, 夏日元, 梁永平, 陈宏峰. 我国南北方岩溶差异[J]. 中国矿业, 2018, 27(S2):317-319. doi: 10.12075/j.issn.1004-4051.2018.S2.077

    LU Haiping, ZHANG Fawang, ZHAO Chunhong, XIA Riyuan, LIANG Yongping, CHEN Hongfeng. Differences between southern karst and northern karst besides scientific issues that need attention[J]. China Mining Magazine, 2018, 27(S2):317-319. doi: 10.12075/j.issn.1004-4051.2018.S2.077

    [45]

    梁永平, 申豪勇, 高旭波. 中国北方岩溶地下水的研究进展[J]. 中国岩溶, 2022, 41(5):199-219.

    LIANG Yongping, SHEN Haoyong, GAO Xubo. Review of research progress of karst groundwater in Northern China[J]. Carsologica Sinica, 2022, 41(5):199-219.

    [46]

    Matonti C, Lamarche J, Guglielmi Y, Marie L. Structural and petrophysical characterization of mixed conduit/seal fault zones in carbonates: Example from the Castellas fault (SE France)[J]. Journal of Structural Geology, 2012, 39:103-121. doi: 10.1016/j.jsg.2012.03.003

    [47]

    Agosta F. Rock physical properties of carbonate fault rocks, Fucino Basin (central Italy): Implications for fault seal in platform carbonates[J]. Geofluids, 2007, 7(1):19-32. doi: 10.1111/j.1468-8123.2006.00158.x

    [48]

    Ferrill D A, Morris A P. Dilational normal faults[J]. Journal of Structural Geology, 2003, 25(2):183-196. doi: 10.1016/S0191-8141(02)00029-9

    [49]

    Childs C, Manzocchi T, Walsh J J, Bonson C G, Nicol A, Schopfer M P J. A geometric model of fault zone and fault rock thickness variations[J]. Journal of Structural Geology, 2009, 31(2):117-127. doi: 10.1016/j.jsg.2008.08.009

    [50]

    Micarelli L, Benedicto A. Normal fault terminations in limestones from the SE-Basin (France): Implications for fluid flow[C]//Wibberley C A J, Kurz W, Imber J. The internal structure of fault zones: Implications for mechanical and fluid-flow properties. Geological Society, 2016: 123-138.

    [51]

    Peacock D C P, Xing Z. Field examples and numerical modelling of oversteps and bends along normal faults in cross-section[J]. Tectonophysics, 1994, 234(1-2):147-167. doi: 10.1016/0040-1951(94)90209-7

  • 期刊类型引用(2)

    1.  李贵仁,刘大金,赵珍,折书群. 矿山地下水三维空间流场模型构建及应用. 矿业研究与开发. 2024(03): 152-158 . 百度学术
    2.  金文正. 云南省洱源县断裂特征及其对地热的控制作用. 中国岩溶. 2024(01): 57-71 . 本站查看

    其他类型引用(1)

  • 加载中

(7)

(1)

计量
  • 文章访问数:  1951
  • PDF下载数:  408
  • 施引文献:  3
出版历程
收稿日期:  2021-09-23
刊出日期:  2022-12-25

目录

  • 表 1.  鄂尔多斯盆地岩石物性特征表(据李冰等,2019宁媛丽等,2020
    Table 1.  The petrophysical property in Ordos basin
    地层/岩体 岩石密度ρ(g/cm3磁化率(10−5 SI)
    密度值 平均值 密度值 平均值
    新生界第四系1.63~2.101.8736~7956
    古近系
    新近系
    2.30~2.392.385~139.6
    中生界白垩系2.4010.5~58.918.6
    侏罗系2.41~2.442.436.5~17.310.5
    三叠系2.558.8~18.615.3
    古生界二叠系2.602.5~22.016.3
    石炭系3.2~20.89.6
    奥陶系2.705.6~14.68.8
    寒武系10.5~16.812.7
    元古界2.705.4~9.27.8
    太古界变粒岩2.761300~106005600
    片麻岩8.6~1400780
    混合岩10.5~16.813.6
    大理岩11.6~24518.4
    侵入岩闪长岩2.5~3.340~95455332
    花岗闪长岩2.562.561374.18
    花岗岩2.580~60483206
    玄武岩2.8~3.32.6910~27131250
    二长石2.47~2.712.59771.16
     | Show Table
    DownLoad: CSV