留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

2021年  第40卷  第3期

目次
目次
2021, 40(3).
摘要:
深部地热开采专栏
干热岩勘查评价指标与形成条件
刘德民, 张昌生, 孙明行, 韦梅华, 关俊朋, 康志强, 熊爱民, 周天禹
2021, 40(3): 1-11. doi: 10.19509/j.cnki.dzkq.2021.0316
摘要:
干热岩具有利用率高、无污染、储量巨大、分布较广、持续稳定、安全性好等特点,被全球公认为在21世纪能够取代化石能源的一种最优质的、可再生的新型清洁能源。目前对于干热岩的成因机理还没有定论,而对于干热岩的勘探寻找、远景区的圈定以及资源评价也有不同的观点,为了有利于干热岩勘探开发,本文在综述前人研究工作基础上,总结归纳了干热岩勘查评价指标,几个重要的指标及其特征如下。第一个指标是岩石圈厚度和莫霍面埋深,岩石圈厚度较小并且莫霍面埋深比较浅,是评价干热岩远景区的一个重要指标。莫霍面埋深较浅指示深部热源(幔源热)更加接近地表,故埋深较浅(较薄的厚度)并具有上隆的特征有利于深部的热量向上传导,为干热岩的孕热环境提供良好的条件。第二个指标是居里等温面,埋深较浅的居里等温面是干热岩远景区一个重要的评价指标。居里等温面是地球内部一个非常特殊的温度(热物质)界面,它不仅能指示地下温度场的分布特征,还可指示地壳深部热能分布特征,对干热岩及地震的成因研究具有十分重要的意义;如果居里面埋深较浅,则热量传导到地表的距离比较短,深部的热流活动更容易向地表传送,不仅是有利的高温干热岩孕热环境,也有利于储存的热能快速向上传导。第三个指标为地温梯度,地温梯度较大是寻找干热岩远景区的一个重要指标;如果一个地区具有较高的地温梯度(≥ 35℃/km),则随着深度的增加深部的地温增加较快,在相对较浅的地方就可以获得温度较高的岩石体。第四个指标是大地热流,大地热流值较大(≥ 75 mW/m2)就指示地球深部有存在高温岩石体的可能;大地热流是地温场的综合性热参数,能够准确地反映区域内的地温场特征。第五个指标是新构造运动,这也是人们寻找干热岩时容易忽略的一个重要指标。新构造运动包括火山、地震及活动断裂构造等。地震和火山是极具破坏性的自然灾害,两者的发生都表明了地球内部的热能汇聚到一定程度,从而打破了地球内部平衡状态而以地震或火山的形式把热量进行释放的地球系统行为。该指标中,如地震震级大(>3级)、震源深度浅(10~15 km)、频度大,火山活动时间新(活火山、休眠火山、中新世以来的死火山)规模大都说明地球深部存在不稳定的高热状态,易形成干热岩;如果能提前找到该区域的干热岩,可以先取出其中的热,那么地震和火山就有可能不会发生,这样可以达到取热减灾减排的作用。活动的深大断裂即能产生一部分热,也能将深部的热传输到浅部,尤其是活动性强的走滑拉张断裂,其深部具有韧性剪切特征,直接指示了深部的高温体的存在。第六个指标是高温温泉与气田等。温泉、气田的形成通常与深部的热储关系密切,一般认为地下水沿某个通道向下渗透接触到深部高温热储被加热后再沿某一通道流出地表而形成温泉;所以,温泉的出露指示了深部存在高温的岩石体(干热岩);区域地温场异常明显,地表热泉等高温水热型地热田较密集的区域有望在深部寻找到干热岩,这也是一种就热(水热)找热(干热)的常规方法。作为固体矿产资源的干热岩,其形成具备四个必要条件:源、通、储、盖。第一个条件是要有丰富的动态热源如来源于深部地幔(幔源热),来源于晚新生代活动的控热构造系统-活动的韧性剪切带,来源于地壳内的低速低阻体(中下地壳热)以及来源于高放射性中新生代花岗岩体(壳源热)。第二个条件是要有优良的导热通道,如壳内15~25 km低速层不仅是热源,同时具有将深部地幔热能向上传导的作用;软流圈地幔上隆时具活动性的深大断裂(深部具韧性特征、浅部具脆性特征)常常具有很好的导热功能;地壳浅表层次的脆性断裂系统往往不是干热岩的热通道,而是水热型地热能的导水、释热构造。第三个条件是要有巨大的储热岩石体,除埋深要适中(3~6 km)并具有较高温度(≥150℃)外,其规模要大(蕴含丰富的热能),热导率大(>2 W/mK),裂隙少(不含水或含少量不流动的水);当然热储层可以是变质岩、岩浆岩,也可以为沉积岩。第四个条件是要有良好的保热盖层,盖层(被子)导热率低(< 2 W/mK)、厚度适中(>1 km)(具有良好的保温效果),地温梯度高(≥40℃/km)、大地热流值高(≥70 mW/m2)(指示深部存在高温特征)是深部赋存有高温地热资源的必要条件。
雄安新区高阳地热田东北部深部古潜山聚热机制
王思琪, 张保建, 李燕燕, 邢一飞, 袁文真, 李郡, 高俊, 赵甜
2021, 40(3): 12-21. doi: 10.19509/j.cnki.dzkq.2021.0319
摘要:
为查明雄安新区深部碳酸盐岩热储的赋存特征与聚热机制,中国地质调查局近两年组织实施的雄安新区地热清洁能源调查评价工作,综合地温场背景、深部地质结构、凹凸相间构造格局、深浅断裂系统控热导热、地球物理场特征、区域水动力场、水化学场与温度场的多场耦合指示等综合探测手段,分析总结了高阳地热田东北部深部古潜山地热资源的最新勘查成果,阐述了雄安新区深部高温古潜山热田的聚热机制:①受太平洋板块向西俯冲影响,华北克拉通东部破坏明显,岩石圈、地壳厚度大嵋减薄,有利于幔源热向浅部的传导;②区域马西断层等深大走滑断裂切穿了岩石圈,有利于深部的幔源热物质和深部岩浆上侵入地壳;③热流由低热导率的凹陷区向高热导率的凸起区聚集,地下水沿断裂的深循环对流加热,区域岩溶地下水的汇流区等也有一定的聚热效应。在雄安新区深部碳酸盐岩热储聚热机制理论指导下,在高阳地热田东北部古潜山布设施工了两眼高温、高产地热井,井口水温分别达到109.2℃,123.4℃,揭示验证了雄安新区深部存在地热开发的第二空间——蓟县系碳酸盐岩热储,深部碳酸盐岩地热水温度高、流量大、开发潜力巨大。研究成果对雄安新区乃至华北盆地高温、高产地热井位的布设与深部碳酸盐岩热储开发利用的科学部署具有很好的示范与指导意义。
涿鹿盆地三维多裂隙地质模型地温场数值模拟
陈金龙, 罗文行, 窦斌, 周洋, 宁文涛
2021, 40(3): 22-33. doi: 10.19509/j.cnki.dzkq.2021.0317
摘要:
地温场预测是地热资源储量评价的重要前提,深部裂隙岩体渗流-传热耦合温度场研究正受到越来越多的关注。涿鹿盆地是新生代地堑式断陷盆地,地表温泉出露点多,地热资源丰富。依据实际地质调查结果和3条地球物理勘探测线解译的地质构造剖面,建立了矿场尺度的涿鹿盆地三维不平行多裂隙地质模型和地热成因模型,赋值调整参数,得到考虑地下水流动情况下多裂隙岩体渗流-传热数学模型,模型的合理性得到验证后,模拟计算了盆地内部温度场,指出盆地地热异常区集中在盆地中心蓟县系(70℃左右)和盆地北东长城系部位(90~98℃),最后研究分析了水力梯度的变化对盆地内部温度场的影响规律。研究结果为涿鹿盆地或同类型地区地热资源勘探、评价和开发利用提供一定的依据。
西藏错那地热田水化学特征与物源机制
余浩文, 刘昭, 荣峰, 陈康, 男达瓦, 刘仕娟, 刘绍赟
2021, 40(3): 34-44. doi: 10.19509/j.cnki.dzkq.2021.0318
摘要:
错那县位于西藏自治区山南地区南部、喜马拉雅山脉中段,是我国重要的边境口岸。为了研究错那地区地热田的水化学特征与物源机制,通过对研究区的温泉点水样的水化学数据进行化验分析,结合研究区的地热地质条件,使用piper三线图和Gibbs图等分析了水化学特征和物质来源。在阐明几种主要的地热温度计原理和适用条件的基础上,利用各种地球化学温度计和多矿物平衡图解法对研究区地下热水进行了温度估算。基于综合分析,最终得出热储温度为117℃左右,分析出地下热水的水化学类型以HCO3-Na型为主,并进行了相应的地热资源评价。这有助于查明地下热水物质来源,探明该地区地热资源,以便在今后为该地区的地热资源开发利用提供技术支持,同时也推动了生态环境的保护和能源结构的逐步改变。
CO2基增强型地热系统中流体-花岗岩相互作用研究进展及展望
薛卉, 舒彪, 陈科平, 路伟, 张森, 胡永鹏
2021, 40(3): 45-53. doi: 10.19509/j.cnki.dzkq.2021.0021
摘要:
CO2与水相比,膨胀性大、黏度低、与岩石反应程度低,并且在作为增强型地热系统(EGS)渗流传热流体时,比水具有更高的换热效率。对CO2-EGS生产过程中储层岩石物性变化的研究具有重要意义,从理论研究、实验研究、数值模拟3个方面,对CO2基增强型地热系统CO2-EGS中流体-岩石相互作用的研究现状进行了总结,并且从矿物成分、微观孔结构和力学性质3个方面对储层岩石性质的变化进行了评价。结果表明,CO2-水-岩石相互作用参与反应的矿物主要为石英、长石类;而沉淀的矿物为蒙脱石、伊利石及方解石等。CO2-水-岩石相互作用会导致储层岩石的力学性质下降,孔隙结构特征改变。最后,讨论了CO2作为EGS渗流传热流体仍需攻克的难点问题,包括:CO2的热动力学特征、换热效率,CO2-水-岩石的相互作用及其对储层性质的改变,影响CO2-EGS经济性的因素,以及CO2-EGS数值模型的研究等。针对这些方面的研究可为今后CO2-EGS的开发奠定基础。
资源地质与工程
综合沉积正演与多点地质统计模拟辫状河三角洲: 以塔河T区为例
贺婷婷, 谈心, 段太忠, 赵磊, 张文彪, 乔勇, 刘彦锋
2021, 40(3): 54-66. doi: 10.19509/j.cnki.dzkq.2021.0301
摘要:
为有效提高三角洲沉积相模拟的精度,以塔河T区为例,应用地震、岩心及测井资料,通过定量评估可容纳空间、物源供给和沉积搬运之间复杂的关系,建立该地区辫状河三角洲砂体沉积正演模型,并将此转化为多点地质统计模拟的三维训练图像,进行研究区辫状河三角洲多点地质统计模拟。研究表明:辫状河三角洲砂体分布受沉积正演模拟控制参数影响,其中砂质供应体积分数、沉积物流入比例、水流载荷量及洪水期间隔均属于强敏感性参数,对研究区三角洲砂体分布影响较大;砂岩含量沉积正演模型体现了辫状河三角洲退积式沉积过程,符合研究区沉积特征,将其转化为三维训练图像,多点地质统计模拟结果在三维空间具有受训练图像约束的特征,体现了训练图像所反映的相带接触关系,并且与钻井认识一致。所提出的这种新的获取三角洲沉积体三维训练图像的方法,综合了沉积正演模拟与多点地质统计模拟的优势,取得了一定的应用效果,对类似沉积体三维地质建模具有一定的借鉴作用。
龙马溪组层状页岩微观非均质性及力学各向异性特征
解经宇, 陆洪智, 陈磊, 金显鹏, 王丹, 付国强
2021, 40(3): 67-77. doi: 10.19509/j.cnki.dzkq.2021.0302
摘要:
层状页岩的微观非均质性及力学各向异性对研究井壁稳定以及水力裂缝扩展形态具有重要意义。为了向页岩优化钻井、压裂工艺参数提供一定的理论和试验依据,对沿不同角度取心的页岩试样开展单轴压缩实验,配合场发射扫描电镜、原子力显微镜观测实验和波速测试等,研究龙马溪组层状页岩微观非均质性及力学各向异性特征,并讨论这些物理力学特征对水力裂缝形态的影响规律。结果表明:受层理面的影响,龙马溪组页岩呈现出较强的微观非均质性和宏观力学各向异性特征。具体的,微观孔隙结构特征方面,随着观测方向与层理方向之间夹角β的增大,微观孔隙结构的发育程度逐渐增加,说明气体的储集和空间呈增加趋势;宏观力学特征方面,单轴压缩条件下,随着加载方向与层理方向间夹角θ的增加,页岩试样的破坏模式从贯穿层理面的张拉破坏,先转变为剪切破坏,再变为劈裂-剪切混合破坏;龙马溪组层状页岩的单轴抗压强度、泊松比随着θ的增加呈现出先减小后增大的"U"形各向异性模式,弹性模量、横纵波速则逐渐减小,胶结程度较弱的页岩层理面会先于基质体发生破坏,进而显著影响岩石整体的力学性质;页岩微观非均质性及力学各向异性特征在一定程度上影响压裂过程中水力裂缝的扩展行为,以及停泵后压裂液的渗流路径。研究结果可为页岩压裂工艺参数优选提供一定依据。
缝洞型碳酸盐岩靶向酸压目标体分类与建模
宋志峰, 张建光
2021, 40(3): 78-84. doi: 10.19509/j.cnki.dzkq.2021.0303
摘要:
塔河油田碳酸盐岩井周缝洞型储集体靶向酸压改造工艺针对性强,对于靶向目标精准定位与表征要求很高。围绕井周酸压增储提效,精准定量表征井筒附近储集体主应力方位、距离、规模、空间分布成为亟待解决的关键问题。从井孔与井周储层识别角度切入,采用一维资料-三维资料,多方法结合定量识别与刻画不同储层类型,在分类表征基础上,实现了模型的分类整合,构建了单元级次的多类型缝洞融合储层模型。然后,在单元级次模型基础上,创建了井周不同范围的可视化三维定量模型构建技术方法,建立了实际单井定量模型,并基于酸压靶向目标体三级八因素分类划分指标体系,确定了7种类型靶向目标体及适应的酸压技术方法。靶向酸压工艺矿场应用显示:靶向酸压井措施效果显著,同样区块范围井,比常规酸压方法平均增产倍比大,显著提高了酸压工艺的针对性和靶向性,有较好的推广应用前景。
珠江口盆地白云-荔湾深水区幔源CO2充注的黏土矿物成岩响应
刘妍鷨, 陈红汉, 王艳飞, 韩晋阳, 李倩
2021, 40(3): 85-95. doi: 10.19509/j.cnki.dzkq.2021.0020
摘要:
珠江口盆地白云-荔湾深水区油气勘探中钻遇大量高含量CO2气层,使得如何规避高CO2风险成为当前面临的挑战。通过工区充注CO2的储层砂岩黏土矿物X射线衍射(XRD)、流体包裹体显微测温和稀有气体同位素组成等分析,明确了幔源无机CO2的成因,并运用黏土矿物特征讨论高温幔源无机CO2充注的黏土矿物成岩响应:有序带倒转和大量自生高岭石的生成。依据CO2气井的XRD全定量分析数据剖面图将幔源无机CO2气层段的黏土矿物组合分为3类:Ⅰ类为I/S混层黏土矿物中S体积分数介于10%~25%之间,有序度处于R2~R3带,高岭石体积分数高;Ⅱ类为I/S混层黏土矿物中S体积分数介于10%~15%,有序度为R2~R3带,自生高岭石体积分数低;Ⅲ类为I/S混层黏土矿物中S体积分数介于50%~60%,有序度为R0带,高岭石体积分数低。总结出幔源CO2在珠江组和珠海组中存在2种运聚模式:一是幔源CO2通过与基性岩浆连通的深大断裂垂直向上运聚;二是垂直向上运移后通过连通砂体侧向长距离运聚。研究结果对珠江口盆地白云-荔湾深水区下一步油气勘探规避高含量CO2风险具有指导作用。
地震沉积学在河流-浅水三角洲沉积相研究中的应用: 以渤海海域蓬莱A构造区馆陶组为例
付鑫, 杜晓峰, 官大勇, 李瑾, 李晓辉
2021, 40(3): 96-108. doi: 10.19509/j.cnki.dzkq.2021.0304
摘要:
庙西南凸起蓬莱A构造区馆陶组有良好油气发现,但储层厚度多小于5 m,难以被精确识别刻画。以薄层砂体为研究对象,综合地震、测井、钻井岩心、古生物等基础资料,明确研究区馆陶组沉积时期受北东向物源供给沉积,发育河流-浅水三角洲相。基于90°相位化处理确定砂体地震响应特征;分频解释处理判断相对等时地震同相轴,构建馆陶组高频层序地层格架(三级层序SQ1-3),优选60Hz地震数据体进行地层切片划分;地层切片高精度沉积学解释揭示研究区馆陶组SQ1层序至SQ3层序依次发育河流-席状化浅水三角洲-朵叶状浅水三角洲-水下河道型浅水三角洲。其中SQ1层序中河道砂体以及SQ2-3层序中浅水三角洲前缘砂体皆为有利储层。此次研究丰富发展了渤海海域河流-浅水三角洲相中地震沉积学在砂体预测中的研究成果。
九瑞矿集区鸡笼山铜金矿床含矿岩体成因和地球动力学背景
詹小飞, 魏俊浩
2021, 40(3): 109-123. doi: 10.19509/j.cnki.dzkq.2021.0023
摘要:
为了研究九瑞矿集区鸡笼山铜金矿床含矿岩体成因和地球动力学背景,对矿集区内成矿关系十分密切的花岗闪长斑岩体进行了主微量元素、锆石U-Pb年代学和Sr-Nd同位素测试。鸡笼山花岗闪长斑岩具有相对较高的SiO2含量和碱质,相对富集Al2O3、CaO,贫MgO、TiO2和P2O5。LREE相对富集、HREE较为亏损的右倾特征,无明显的负Eu异常,相对富集Rb、Ba、K、Sr等大离子亲石元素(LILE)和亏损Ta、Nb等高场强元素(HFSE),获得样品的形成年龄为(140±1)Ma(MSWD=1.02),表明其均为高钾钙碱性准铝质系列,属于壳幔混合成因的I型花岗岩,岩体可能为富集岩石圈地幔在部分熔融过程中混染了一定程度的古老下地壳物质的产物。岩体侵入年代与长江中下游地区中生代大规模岩浆活动时限一致,可能与古太平洋板块的俯冲作用和区域岩石圈伸展减薄条件下深部壳幔物质相互作用有关。
鸡西盆地主力煤层水可动性及其孔渗控制
周文宇, 王小明, 曾凡桂, 党正, 朱冲, 陈文文, 王志壮, 涂明恺
2021, 40(3): 124-131. doi: 10.19509/j.cnki.dzkq.2021.0305
摘要:
水可动性是影响煤层气产出的重要因素,分析孔渗对水可动性的作用对鸡西盆地煤层气开发具有重要意义。以鸡西盆地不同矿区主力煤层为研究对象,开展了低场核磁共振(NMR)及渗透率实验,同时,结合称重测定煤样含水饱和度的方法,分析了煤层孔渗特征对煤中水可动性的影响。结果表明:①研究区煤样中吸附孔较为发育,平均占比为58.36%;渗流孔和裂隙发育程度相当,平均占比分别为21.23%和20.41%,且两者之间连通性较好;②确定了煤样达到束缚水状态的离心力为1.38 MPa,此压力下煤中可排出半径约0.1 μm半开放或开放孔隙中的可动水。煤样可动水饱和度为27.84%~60.87%,平均值37.86%;束缚水饱和度为39.13%~72.16%,平均值62.14%,可动水含量相对束缚水含量较低。随着离心压力的增加,可动水首先沿着大裂隙流动,随后经过渗流孔流出,束缚水主要存在于吸附孔中,需要克服较大的毛细管阻力,难以流动;③尽管煤样变质程度相近,但其内部水的可动性及赋存状态存在明显差异,这可能是造成同一区块不同煤层气井气水产出差异的原因之一;④研究区煤中半径>1 000 nm段孔喉越发育、煤岩渗透率越大,煤的可动水饱和度越高,水在其中越容易流出;煤中半径0~100 nm段孔喉越发育,煤的可动水饱和度越低,水在其中被束缚而难以流动。
温压条件对非饱和低渗透砂岩中CH4突破压力影响的实验研究
赵世宇, 李铱, 明亮, 张勇, 李旭峰, 郑长远
2021, 40(3): 132-139. doi: 10.19509/j.cnki.dzkq.2021.0306
摘要:
突破压力在气藏开采和盖层封闭性评价中起着重要的作用。为模拟地层温度压力变化对低渗透性砂岩中CH4突破压力的影响,采用逐步法对取自鄂尔多斯盆地早二叠世非饱和低渗透砂岩,进行了不同温压组合条件下的CH4突破压力实验。结果表明:随着温度或压力的升高,CH4突破压力和突破时间均呈下降趋势,且压力变化对CH4突破过程的影响更加显著。分析发现,温度压力对CH4突破压力的影响,是由于温压变化改变了界面张力和接触角的大小,进而控制突破压力的大小;另外,突破压力随着两相黏度比的增大而减小,且温度压力越高,黏度比对突破压力的影响越小。因此,实验温压范围内,对于非饱和低渗透砂岩储层,温度和压力越高越有利于气藏开采;对于非饱和低渗透砂岩盖层,温度和压力越低,盖层封闭性越好且越安全。
松辽盆地大庆长垣南缘四方台组砂岩型铀矿碳氧同位素特征及铀成矿意义
蔡宁宁, 邹妞妞, 付勇, 魏帅超, 李建国, 汤超, 魏佳林
2021, 40(3): 140-150. doi: 10.19509/j.cnki.dzkq.2021.0307
摘要:
为进一步探讨松辽盆地大庆长垣南缘四方台组砂岩型铀矿的成岩环境、铀沉淀机制及成矿机理,通过钻井岩心观察,岩石薄片,含铀砂岩碳酸盐胶结物碳、氧同位素分析,对研究区砂岩岩石类型、碳酸盐胶结物类型、成岩阶段以及碳、氧同位素特征进行了系统研究。结果表明:含铀砂岩岩石类型为长石岩屑砂岩、岩屑长石砂岩,碳酸盐胶结物可分为3个期次,以第一和第二期次为主,胶结物包括方解石、铁方解石,见极少量含铁白云石。δ13CPDB值为-22.45‰~-13.65‰,平均值为-18.33‰,δ18OPDB值为-17.56‰~-9.46‰,平均值为-13.52‰,揭示有一定有机质参与;古盐度Z值介于75.98~90.97之间,为淡水沉积;古温度介于67.31~94.92℃之间,结合以上岩石特征、胶结物类型及成岩作用判定成岩环境处于早成岩—中成岩阶段A期。综合分析认为:早成岩阶段砂岩的压实、胶结作用使含铀砂岩孔隙度减小,抑制含氧富铀流体运移;随着成岩作用进行,埋藏深度增加,有机质成熟度升高发生热脱羧基作用,与压实作用产生的吸附水形成酸性溶液对长石、碳酸盐胶结物等进行溶蚀形成次生孔隙,增大砂岩孔隙度,促进含氧富铀流体下渗运移,下伏油气由断裂向上逸散,为地层提供还原剂,使铀在四方台组底部大规模沉淀和富集。
南岭及邻区第四纪地层分区与综合地层格架
史元润, 林晓, 徐亚东, 余舒超, 孙清泉
2021, 40(3): 151-162. doi: 10.19509/j.cnki.dzkq.2021.0309
摘要:
南岭及邻区的第四纪沉积以冲积物、残积物和洞穴堆积物为主,长期以来缺乏统一的地层格架和时间框架。基于剖面资料和实地调查,将南岭第四纪地层划分为永州-郴州、桂林-贺州、韶关-清远、赣州、道县-阳山等5个地层小区;通过厘定每个小区的岩石地层序列、生物地层和气候地层及其测年成果,建立了综合地层的对比格架。结果表明,南岭第四系与中国更新统泥河湾阶、周口店阶、萨拉乌苏阶和尚未建阶的全新统可一一对比;其生物地层以早更新世巨猿动物群和中更新世晚期以来的马坝人、道县人和柳江人等智人演化为特点;气候地层以洞穴石笋和高山泥炭重建的古气候记录为代表,主要反映东亚季风背景下的南岭局地气候。但由于测年数据较少,冲积层和残积层的对比仍然存在困难;生物演化阶段与气候变化之间也未能详细对比。因此,今后南岭的第四纪地层研究应围绕提高测年精度与建立统一的年代标尺开展,促进综合地层对比的精细化、延伸气候地层的时间尺度并拓展与南海海相地层的对比;这一基础性工作对水文、环境、工程等实践应用亦具有重要意义。
工程地质与地质灾害
单轴压缩条件下不同含水率黑云母二长花岗岩破坏特征与机制
曹洋兵, 陈玉华, 张朋, 黄真萍, 张向向, 陈杨涛
2021, 40(3): 163-172. doi: 10.19509/j.cnki.dzkq.2021.0308
摘要:
花岗岩在不同含水率条件下的变形破坏特征和机制对此类工程岩体稳定性评价具有重要的意义。开展不同含水率黑云母二长花岗岩单轴压缩试验,分析破坏特征和应力-应变曲线特征,开展断口扫描电镜试验,分析微观形貌特征,研究破坏机制。试验结果表明:黑云母二长花岗岩具有明显的应变软化特征;随含水率增大,曲线上微裂隙压密阶段长度逐渐增加,稳定破裂阶段及非稳定破裂阶段长度均逐渐缩短,但所占比例增大,曲线上峰前阶段涨落交替现象加剧;饱和时单轴抗压强度和弹性模量相比干燥时分别降低了40.68%,20.3%;变形破坏过程可大致分为以下5个阶段:平静期、裂纹萌生期、裂纹扩展伴随颗粒弹射期、片状碎片剥落伴随颗粒弹射期及崩落式破坏期;随含水率增大,花岗岩破坏时的剧烈程度、发出的声响及脆性程度均逐渐降低;花岗岩破坏机制为拉-剪复合破坏,低含水率时以压致拉张破坏为主,随含水率增大呈现拉张破坏减少而剪切破坏增多的趋势,饱和时以剪破坏为主。研究结果可为黑云二长花岗岩与水之间的耦合模型构建提供理论支撑,对水-岩耦合环境下工程岩体稳定性分析具有重要科学意义。
基于线激光扫描的岩石激光钻孔的三维重建和可视化
官东林, 文国军, 王玉丹, 童志伟, 吴玲玲, 罗耀坤
2021, 40(3): 173-183. doi: 10.19509/j.cnki.dzkq.2021.0310
摘要:
激光钻进岩石形成的钻孔的孔形较为复杂,具有较小的孔直径和较高的孔壁粗糙度,使得利用传统方法进行钻孔尺寸的测量较为困难。为了精确钻孔测量和方便孔形研究,提出了一种基于线激光扫描及逆向建模的钻孔建模方法。首先,搭建了线激光扫描平台,建立了空间坐标系,以获取钻孔的三维坐标,构建了钻孔的初始点云数据。其次,在MATLAB中对获取的点云数据进行无效点移除及多视角点云配准,其中,无效点移除利用顺序查找法实现,多视角点云配准则基于迭代最近点(ICP)算法,包括初始配准和精确配准两个阶段。最后,基于Delaunay三角网格划分及曲面重建算法,实现了钻孔模型的重建和可视化。此外,还采用滴液法和切割法进行实际钻孔容积值测量及钻孔轮廓线获取,并与由点云重建的钻孔模型上获取的测算结果进行对比分析,以验证所述方法建立的钻孔模型的精度。结果表明:重建的钻孔模型与实际钻孔之间的误差小于4%,重建的模型能够满足激光岩石钻进钻孔的测量要求,证实了所述方法的可行性。与传统测量方法相比,所述方法属于非接触、非破坏性方法,可重复性测量。
原状及重塑黄土双轴试验微观力学特征离散元模拟
井旭, 谢婉丽, 单帅
2021, 40(3): 184-193. doi: 10.19509/j.cnki.dzkq.2021.0311
摘要:
黄土作为一种特殊的颗粒材料,微观上颗粒组成的结构决定了其力学特性。原状及重塑黄土因结构的差异而具有不同的力学特性。针对黄土结构性如何影响其力学特征这一基本问题,开展基于电镜扫描获取细观颗粒信息,同时考虑颗粒形状、颗粒破碎可能性进行建模的离散单元法进行原状黄土和饱和重塑土在恒定应变速率双轴试验下的宏观力学和细观力学性能研究。研究结果显示:试样微观结构的差异对变形破坏过程产生显著影响。当轴向应力较低时原状黄土及重塑黄土力链多分布于大型骨架颗粒附近,随着轴向应力增加,原状黄土力链形成网状图案但仍具备主要传导区域,重塑黄土无明显主要传导区,呈现均匀网状。原状土及重塑土骨架颗粒簇周围多形成张拉裂隙,剪切裂隙多数形成于骨架颗粒簇内部,又以颗粒簇相互挤压接触时最为明显。使用该建模方法,可以有效反映原状及重塑黄土由于内部结构组成差异导致相同应力条件下产生的不同内部应力状态。基于以上研究结论,给出了黄土结构性对宏观强度影响的微观解释。研究成果可为黄土地区地质灾害防治提供一定依据。使用该建模方法,可以有效反映原状及重塑黄土由于内部结构组成差异导致相同应力条件下产生的不同内部应力状态。基于以上研究结论,给出了黄土结构性对宏观强度影响的微观解释。
环境地质与水文地质
新疆阿克苏典型山前洪积扇内高氟地下水的化学特征及氟富集机制
潘欢迎, 邹常健, 毕俊擘, 刘运德, 黄丽文
2021, 40(3): 194-203. doi: 10.19509/j.cnki.dzkq.2021.0312
摘要:
在内陆干旱区,作为重要饮用水源的地下水常面临氟含量超标问题。查明内陆干旱区高氟地下水的分布规律,了解氟在地下水中的富集过程及其影响因素,既可丰富高氟地下水的研究体系,也是保证内陆干旱区饮水安全的重要基础。以新疆阿克苏地区典型山前洪积扇——依格齐艾肯河-喀拉玉尔滚河河间地带为研究区,基于水文地球化学调查结果,刻画了高氟地下水的分布区;结合氟离子含量与特征性水化学指标间的关系,揭示了高氟地下水的成因机制。结果表明:①地下水中氟含量的变化范围为0.8~6.1 mg/L,83%的水样氟含量超过《生活饮用水卫生标准》(GB 5749-2006)规定的上限(1.0 mg/L);②总体上,氟含量沿地下水流动路径逐渐增大,低氟地下水(ρ(F-)≤1.0 mg/L)分布在国道314以北的补给区,高氟地下水(ρ(F-)>1.0 mg/L)分布在国道314以南的径流区和排泄区;③高氟地下水的水化学类型以Cl·HCO3-Na型为主,而低氟地下水则以Cl·SO4-Na型为主,高氟地下水相比于低氟地下水优势阴离子偏向于HCO3-;④地下水的pH值范围为7.9~8.9(均值为8.4),表明其处于弱碱环境中。地下水中ρ(F-)与pH值呈正相关,此外构成浅层含水层的上更新统沉积物中含有黑云母、氟磷灰石等矿物,其表面存在一定数量的可交换F-,这表明水中OH-与矿物表面F-间的阴离子交换可能对氟的富集有一定贡献;⑤地下水的F-含量与Ca2+含量呈负相关,即高氟地下水中ρ(Ca2+)小于低氟地下水。考虑到氟化钙(CaF2)是自然界中的主要含氟矿物,也是地下水中氟的主要来源,ρ(F-)与ρ(Ca2+)间的这种负相关指示着高氟地下水中可能存在去Ca2+、Mg2+作用,如阳离子交替吸附或碳酸盐岩沉淀等。研究区地下水样中ρ(F-)与ρ(Mg2+)间也呈负相关关系,且和ρ(F-)与ρ(Ca2+)间的关系高度相似,也佐证了高氟地下水中去Ca2+、Mg2+作用的存在;⑥绝大部分地下水样品都位于氯碱性指数图的负值区域,且ρ(F-)与CAI-1和CAI-2均呈较好负相关,CAI-1和CAI-2都随ρ(F-)的增大而减小,这表明高氟地下水中存在Ca2+、Mg2+与Na+间更强的交换作用,对氟富集起着重要作用。地下水中ρ(F-)与SAR间呈较好正相关关系,且高氟地下水样的SAR均值(5.71)远大于低氟地下水SAR均值(1.67),这也进一步证明高氟地下水中的Ca2+、Mg2+与含水介质的Na+间存在强烈的交替作用,对氟的富集起着重要作用;⑦所有地下水样中的萤石均处于未饱和状态,且萤石的饱和指数(SI)与F-含量间呈现较好的正相关,这表明地下水对含氟矿物(主要是萤石)的持续溶解应是导致研究区地下水中氟富集的主要原因。与之相反,研究区所有地下水样中的方解石均处于过饱和状态(SI>0)。这表明CaCO3的沉淀可能促进了CaF2的溶解,导致地下水中氟离子质量浓度增高;⑧研究区低氟地下水的δ18O值介于-11.20‰~-10.67‰间,平均值为-10.94‰,而高氟地下水的δ18O值介于-11.65‰~-11.21‰间,平均值为-11.49‰,即低氟地下水较高氟地下水富集δ18O。此外,F-质量浓度较低(ρ(F-)≤3.0 mg/L)的地下水样中δ18O值与F-质量浓度呈负相关,即低氟地下水具有更正的δ18O值;F-质量浓度较高(ρ(F-)≥4.8 mg/L)的地下水样中δ18O值与F-质量浓度的相关性不显著,随F-质量浓度的增高,δ18O值基本维持不变。以上表明蒸发浓缩作用对地下水中氟的富集贡献较小;⑨研究区地下水中ρ(F-)/ρ(Cl-)比值与ρ(F-)间呈现正相关,即ρ(F-)/ρ(Cl-)比值随ρ(F-)增高呈增大趋势,这也说明地下水中氟富集的主要原因是含氟矿物的溶解,而不是蒸发浓缩作用。此外,Gibbs图也提供了证据:研究区地下水样基本处于水岩作用主导区域,表明地下水化学特征(包括氟的富集)主要受水岩作用控制,蒸发浓缩影响很小。总之,地下水中氟的富集主要由溶解作用引起,OH-与矿物表面F-间的交换也有贡献,但蒸发浓缩作用影响微弱。含氟矿物持续溶解的驱动机制是阳离子交替吸附(地下水中Ca2+与岩土颗粒表面Na+之间)及方解石沉淀所引起的地下水中Ca2+的衰减。
武汉市东西湖区主要湖泊表层沉积物重金属污染特征与生态风险评价
龚杰, 孙紫童, 冯璐, 胡胜华, 敖鸿毅, 吴辰熙
2021, 40(3): 204-210. doi: 10.19509/j.cnki.dzkq.2021.0313
摘要:
人类活动对水环境影响巨大,同时水环境影响着人类的生产生活。采用地质累积指数法及潜在生态危害指数法对武汉市东西湖区主要湖泊表层沉积物重金属情况展开评价发现东西湖主要湖泊表层沉积物平均污染程度处于相对较低水平。地累积指数法评价结果表明,东大湖和金湖中Cd元素处于中度污染水平累积,墨水湖和上金湖中Cd元素为偏中度污染累积水平,其他湖泊沉积物中重金属均处于轻度或无污染水平。潜在生态危害指数法评价结果表明金湖和东大湖沉积物中重金属的潜在风险处于重度生态风险水平,墨水湖、上金湖、下银湖处于中度风险水平,其余湖泊沉积物中重金属的潜在风险均处于低风险水平。主成分分析表明,Hg、Cr、Ni、Cu、Zn、Cd、Pb可能来自于工业和交通源的排放,而As可能与农药的使用有关。后期金湖、东大湖、墨水湖、上金湖和下银湖的内源污染治理过程中需考虑重金属的生态风险问题。
新技术新方法
基于关联规则与随机森林的地震多属性砂体厚度预测
曲志鹏, 王芳芳, 张云银, 李晓晨
2021, 40(3): 211-218. doi: 10.19509/j.cnki.dzkq.2021.0314
摘要:
地震属性技术是砂体厚度预测的重要手段,由于目前可从地震数据中提取的地震属性种类较多,在利用地震属性技术前,必须优化出对砂体厚度最敏感的地震属性组合,以减少地震属性信息的重复与冗余。为此提出了一种联合关联规则与随机森林回归算法的地震多属性砂体厚度预测方法。随机森林回归算法能够建立地震多属性与砂体厚度之间的非线性关系,并能进行属性选择,但是该算法无法识别地震多种属性中的冗余特征。关联规则能够发现地震属性之间的非线性关联,并能借助卡方检验消除地震属性间的冗余性。分别采用了随机森林回归算法(RFR)、联合关联规则与随机森林回归(AR-RFR)及BP神经网络回归的算法(AR-BP)对滩坝砂岩合成模型和某实际工区进行了砂体厚度预测。对比结果表明,基于关联规则的属性优选得到的属性间相关性低,关联规则与随机森林算法的结合提高了砂体厚度的预测精度。数值实验证明了该方法的有效性。
CryoSat DEM与多种南极DEM的对比分析
岳文丽, 董玉森, 阿布都拉·阿不都卡地尔, 崔静月, 张帮政
2021, 40(3): 219-227. doi: 10.19509/j.cnki.dzkq.2021.0315
摘要:
南极洲被巨厚冰雪覆盖,地质构造以南极横断山脉为界,总体分为东南极地盾和西南极活动带。数字高程模型(DEM)是研究南极冰盖变化的基础数据之一。通过多期次数字高程模型相比较获得高程的变化信息,是分析南极冰盖厚度变化和物质平衡的重要手段。然而不同类型DEM之间存的平面误差和垂直误差影响分析结果的精度。首先利用配准消除DEM间的水平误差,然后计算并按坡度提取CryoSat DEM与其他DEM的平均高程差和标准差,最后分析高程差的时空变化特征。通过分析发现,DEM之间存在不同的平面误差。其中TanDEM_X DEM与CryoSat DEM的高程平面偏差最小,而ICESat DEM与CryoSat DEM的高程平面偏差最大。在垂直方向上,0°~1°的坡度范围内,CryoSat DEM与TanDEM_X DEM的平均高程差在3.5~5.5 m之间,标准差小于18.0 m;CryoSat DEM和Bamber 1km DEM的平均高程差在-2.5~+1.0 m之间,标准差小于24.2 m;CryoSat DEM与ICESat DEM的平均高程差在-25.0~-1.0 m之间,标准差小于47.2 m;CryoSat DEM与RAMPv2 DEM的平均高程差在1.3~3.2 m之间,标准差小于45.6 m。通过研究发现南极冰盖内部高程增加,但西南极冰盖和东南极冰盖高程均在降低,且西南极降低明显,同时南极边缘地区高程降低明显。本研究为全球变化研究和南极物质平衡研究提供了重要参考。