ISSN 1000-3665 CN 11-2202/P
  • 中文核心期刊
  • CSCD核心期刊
  • 中科双效期刊
  • 中国科技核心期刊
  • Caj-cd规范获奖期刊
欢迎扫码关注“i环境微平台”

留言板

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

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

优先发表

优先发表栏目展示本刊经同行评议确定正式录用的文章,这些文章目前处在编校过程,尚未确定卷期及页码,但可以根据DOI进行引用。
显示方式:
安阳市地下水源热泵系统建设水资源管理区划研究
朱文举, 平建华, 侯俊山, 宁艺武, 耿文斌
, doi: 10.16030/j.cnki.issn.1000-3665.202011003
摘要:
为从赋存条件和水资源管理两个层面进行地下水源热泵系统建设的区划研究,实施技术与管理的结合,通过分析水文地质条件、水动力条件和水化学条件,结合水资源管理分区,建立了安阳市地下水源热泵系统水资源管理区划评价体系。采用云模型改进的层次分析法进行了一级评价,在此基础上结合水资源管理分区利用GIS空间分析功能完成二级评价,将研究区地下水源热泵系统划分了3个等级。结果表明:研究区范围内地下水源热泵系统适宜发展区面积为117.45 km2,主要分布在安阳河冲洪积扇扇体中心强富水区,部分分布在扇缘的外围区域;限制发展区面积为459.26 km2,分布在扇缘的西南和北部丘陵弱富水区以及扇体中心的地下水降落漏斗区;禁止发展区面积为24.02 km2,分布在水源地和南水北调保护区以及铁路和高速公路两旁,在研究区交错分布。在适宜性分区的基础上结合水资源管理的区划研究更全面合理,可为地下水源热泵系统科学布局及合理的开发利用提供参考。 为从赋存条件和水资源管理两个层面进行地下水源热泵系统建设的区划研究,实施技术与管理的结合,通过分析水文地质条件、水动力条件和水化学条件,结合水资源管理分区,建立了安阳市地下水源热泵系统水资源管理区划评价体系。采用云模型改进的层次分析法进行了一级评价,在此基础上结合水资源管理分区利用GIS空间分析功能完成二级评价,将研究区地下水源热泵系统划分了3个等级。结果表明:研究区范围内地下水源热泵系统适宜发展区面积为117.45 km2,主要分布在安阳河冲洪积扇扇体中心强富水区,部分分布在扇缘的外围区域;限制发展区面积为459.26 km2,分布在扇缘的西南和北部丘陵弱富水区以及扇体中心的地下水降落漏斗区;禁止发展区面积为24.02 km2,分布在水源地和南水北调保护区以及铁路和高速公路两旁,在研究区交错分布。在适宜性分区的基础上结合水资源管理的区划研究更全面合理,可为地下水源热泵系统科学布局及合理的开发利用提供参考。
基于接力排水的强夯法在滨海回填区地基处理中的试验研究
张军舰, 李鹏, 殷坤宇, 罗玉磊, 郭幔
, doi: 10.16030/j.cnki.issn.1000-3665.202104076
摘要:
本文对山东半岛海岸带滨海杂填土、饱和粉细砂、淤泥质土等特殊复杂地层地基处理方法进行了研究。以经济高效的强夯法为基础,提出复杂地层整体排水概念,设计了浅层、深层竖向排水和水平排水的接力排水系统,并进行了现场试验研究。监测数据表明,强夯荷载作用下,接力排水系统整体协同排水,可快速排出各个地层中地下水、消散超孔隙水压力。7 h左右可基本消除强夯引起的地下水上升及孔隙水压力消散。持续降水,地表沉降为上部土体厚度的0.7%~2.0%。强夯动力荷载作用下,表层土体压缩为上部土体厚度的8.7%~10.9%。埋深3~7 m土体沉降约为土体厚度的5‰、3‰,埋深7~10 m土体沉降为土体厚度的2‰。检测数据表明,在强夯有效影响深度内地基处理效果明显,土体工程性状改善明显。表层承载力及变形模量满足设计要求,4 m以下淤泥承载力平均值略低于设计要求,下部淤泥质土计算平均固结度为77%。夯后1个月监测数据表明,地表沉降量在25 mm以内,已逐步趋于稳定,分层沉降、孔隙水压力数值整体稳定略有下降。 本文对山东半岛海岸带滨海杂填土、饱和粉细砂、淤泥质土等特殊复杂地层地基处理方法进行了研究。以经济高效的强夯法为基础,提出复杂地层整体排水概念,设计了浅层、深层竖向排水和水平排水的接力排水系统,并进行了现场试验研究。监测数据表明,强夯荷载作用下,接力排水系统整体协同排水,可快速排出各个地层中地下水、消散超孔隙水压力。7 h左右可基本消除强夯引起的地下水上升及孔隙水压力消散。持续降水,地表沉降为上部土体厚度的0.7%~2.0%。强夯动力荷载作用下,表层土体压缩为上部土体厚度的8.7%~10.9%。埋深3~7 m土体沉降约为土体厚度的5‰、3‰,埋深7~10 m土体沉降为土体厚度的2‰。检测数据表明,在强夯有效影响深度内地基处理效果明显,土体工程性状改善明显。表层承载力及变形模量满足设计要求,4 m以下淤泥承载力平均值略低于设计要求,下部淤泥质土计算平均固结度为77%。夯后1个月监测数据表明,地表沉降量在25 mm以内,已逐步趋于稳定,分层沉降、孔隙水压力数值整体稳定略有下降。
基于物理模型试验的碎屑流流态化运动特征分析
龙艳梅, 宋章, 王玉峰, 程谦恭, 李坤, 吴越
, doi: 10.16030/j.cnki.issn.1000-3665.202011035
摘要:
流态化运动是高速远程滑坡的主要运动形式,是揭示高速远程滑坡运动机理的重要基础。本文基于粒子图像测速(PIV)分析方法,采用物理模型试验对不同粒径组成条件下的颗粒流内部的速度分布特征、剪切变形特征及流态特征进行了研究,并对高速远程滑坡流态化运动特征进行了讨论分析。结果表明,碎屑流流态化运动特征与颗粒粒径呈显著的相关性,随着粒径的减小或细颗粒含量的增加,颗粒流底部相对于边界的滑动速度以及整体的运动速度均呈现逐渐减小的趋势,颗粒流内部剪切变形程度增加,颗粒的运动形式由“滑动”向“流动”转变。当颗粒粒径较小或细颗粒含量较高时,颗粒流内部剪切速率增大的趋势在颗粒流底部更加显著,反映了粒径减小有助于促进颗粒流内部剪切向底部的集中。在同一颗粒流的不同运动阶段及不同纵向深度,其流态特征具有显著差别,颗粒流前缘及尾部主要呈惯性态,颗粒间以碰撞作用为主,而主体部分则主要呈密集态,颗粒间以摩擦接触作用为主;在颗粒流表面及底部,颗粒间相互作用方式主要是碰撞作用,中间部分则以摩擦作用为主。对于不同粒径的颗粒流,随着粒径的增大或粗颗粒含量的增加,颗粒流内部颗粒的碰撞作用加强,颗粒流整体趋于向惯性态转变。 流态化运动是高速远程滑坡的主要运动形式,是揭示高速远程滑坡运动机理的重要基础。本文基于粒子图像测速(PIV)分析方法,采用物理模型试验对不同粒径组成条件下的颗粒流内部的速度分布特征、剪切变形特征及流态特征进行了研究,并对高速远程滑坡流态化运动特征进行了讨论分析。结果表明,碎屑流流态化运动特征与颗粒粒径呈显著的相关性,随着粒径的减小或细颗粒含量的增加,颗粒流底部相对于边界的滑动速度以及整体的运动速度均呈现逐渐减小的趋势,颗粒流内部剪切变形程度增加,颗粒的运动形式由“滑动”向“流动”转变。当颗粒粒径较小或细颗粒含量较高时,颗粒流内部剪切速率增大的趋势在颗粒流底部更加显著,反映了粒径减小有助于促进颗粒流内部剪切向底部的集中。在同一颗粒流的不同运动阶段及不同纵向深度,其流态特征具有显著差别,颗粒流前缘及尾部主要呈惯性态,颗粒间以碰撞作用为主,而主体部分则主要呈密集态,颗粒间以摩擦接触作用为主;在颗粒流表面及底部,颗粒间相互作用方式主要是碰撞作用,中间部分则以摩擦作用为主。对于不同粒径的颗粒流,随着粒径的增大或粗颗粒含量的增加,颗粒流内部颗粒的碰撞作用加强,颗粒流整体趋于向惯性态转变。
软土地区双线区间盾构隧道施工对周边地表以及建筑物沉降的影响
潘涛
, doi: 10.16030/j.cnki.issn.1000-3665.202106014
摘要:
研究盾构隧道施工对周围地面以及建筑物沉降造成的影响,是软土地区盾构隧道安全施工和正常运营的基础课题。为了分析宁波轨道交通5号线同德路站—石碶站区间双线盾构隧道施工对周边地表和建筑物的影响,本文在建立盾构隧道动态施工过程三维有限元模型的基础上,基于地表以及建筑物沉降数值模拟结果与现场监测值的对比,分析了隧道开挖对隧道周围地表沉降与建筑物沉降的影响。结果表明,掘进完成时,开挖方向沉降槽往上行线隧道方向偏移、呈现倒梯形形态,横断面影响区域为距离双线隧道轴线中心小于3倍隧道直径;上行线在下行线开挖后并不会增加地表沉降,但增大了沉降槽宽度;下行线到达前产生的沉降占最终累计沉降的67%;当盾构掘进面刚到达建筑物时、建筑物的倾斜方向与盾构掘进方向一致,当盾构掘进面离开建筑物时、建筑物将沿着盾构掘进的反方向倾斜;建筑物两侧沉降值较中部沉降值降低了83%;双线贯通后建筑物沉降呈“U”形分布,最大沉降量发生在远离隧道一侧距建筑物中心0.5 m处。 研究盾构隧道施工对周围地面以及建筑物沉降造成的影响,是软土地区盾构隧道安全施工和正常运营的基础课题。为了分析宁波轨道交通5号线同德路站—石碶站区间双线盾构隧道施工对周边地表和建筑物的影响,本文在建立盾构隧道动态施工过程三维有限元模型的基础上,基于地表以及建筑物沉降数值模拟结果与现场监测值的对比,分析了隧道开挖对隧道周围地表沉降与建筑物沉降的影响。结果表明,掘进完成时,开挖方向沉降槽往上行线隧道方向偏移、呈现倒梯形形态,横断面影响区域为距离双线隧道轴线中心小于3倍隧道直径;上行线在下行线开挖后并不会增加地表沉降,但增大了沉降槽宽度;下行线到达前产生的沉降占最终累计沉降的67%;当盾构掘进面刚到达建筑物时、建筑物的倾斜方向与盾构掘进方向一致,当盾构掘进面离开建筑物时、建筑物将沿着盾构掘进的反方向倾斜;建筑物两侧沉降值较中部沉降值降低了83%;双线贯通后建筑物沉降呈“U”形分布,最大沉降量发生在远离隧道一侧距建筑物中心0.5 m处。
考虑基质吸力的非饱和土邓肯-张统计损伤修正模型
谭维佳, 魏云杰, 王俊豪, 高敬轩
, doi: 10.16030/j.cnki.issn.1000-3665.202105022
摘要:
为反映非饱和地基土变形全过程,以某基础工程非饱和土为研究对象,开展固结排水三轴压缩试验。研究发现,非饱和土的偏应力-应变曲线形态近似双曲线,基质吸力对土体力学行为影响十分明显,基质吸力越大,土体偏应力越高。根据非饱和土变形特性及工程特点,选取邓肯-张双曲线模型作为基础模型,引入统计损伤理论,假设非饱和土微元强度服从Weibull概率密度分布,建立邓肯-张统计损伤模型。通过搭建初始切线模量与基质吸力的联系,建立一种新的考虑基质吸力的非饱和土邓肯-张统计损伤模型。给出参数解析方法,得到Weibull分布参数经验表达式,从而修正模型。分析不同基质吸力条件下非饱和土损伤累积规律,采用所建模型和传统邓肯-张模型对比验证非饱和土偏应力-应变试验曲线,证明所建模型的可行性和合理性。研究成果为非饱和土的力学特性研究及辨识模拟提供一定参考。 为反映非饱和地基土变形全过程,以某基础工程非饱和土为研究对象,开展固结排水三轴压缩试验。研究发现,非饱和土的偏应力-应变曲线形态近似双曲线,基质吸力对土体力学行为影响十分明显,基质吸力越大,土体偏应力越高。根据非饱和土变形特性及工程特点,选取邓肯-张双曲线模型作为基础模型,引入统计损伤理论,假设非饱和土微元强度服从Weibull概率密度分布,建立邓肯-张统计损伤模型。通过搭建初始切线模量与基质吸力的联系,建立一种新的考虑基质吸力的非饱和土邓肯-张统计损伤模型。给出参数解析方法,得到Weibull分布参数经验表达式,从而修正模型。分析不同基质吸力条件下非饱和土损伤累积规律,采用所建模型和传统邓肯-张模型对比验证非饱和土偏应力-应变试验曲线,证明所建模型的可行性和合理性。研究成果为非饱和土的力学特性研究及辨识模拟提供一定参考。
冻融作用下污泥固化轻质土动力特性及结构演化
杨爱武, 王斌彬, 钟晓凯
, doi: 10.16030/j.cnki.issn.1000-3665.202103013
摘要:
冻融作用会影响污泥固化轻质土动力特性和结构特性。为研究不同冻融循环次数和冻结温度下的污泥固化轻质土动力特性及结构演化,对冻融循环作用后的污泥固化轻质土进行动三轴试验和固结试验。试验结果表明:在冻融循环作用下,污泥固化轻质土动应力-应变曲线呈弱应变硬化型。污泥固化轻质土动强度和变形随着冻融次数(n)的增加和冻结温度(T)的降低而减小。经历前4次冻融作用后,对固化土变形和动强度影响较大,经历8次冻融循环后基本趋于稳定。冻融循环对土体动力特性的影响本质上是对土体结构性的影响。固化土冻融结构势(mdσn)的变化规律与动强度类似,随冻融次数的增加和冻结温度的降低而下降。冻融循环次数是影响污泥固化轻质土动力特性和结构性的主要因素,冻结温度为次要因素。 冻融作用会影响污泥固化轻质土动力特性和结构特性。为研究不同冻融循环次数和冻结温度下的污泥固化轻质土动力特性及结构演化,对冻融循环作用后的污泥固化轻质土进行动三轴试验和固结试验。试验结果表明:在冻融循环作用下,污泥固化轻质土动应力-应变曲线呈弱应变硬化型。污泥固化轻质土动强度和变形随着冻融次数(n)的增加和冻结温度(T)的降低而减小。经历前4次冻融作用后,对固化土变形和动强度影响较大,经历8次冻融循环后基本趋于稳定。冻融循环对土体动力特性的影响本质上是对土体结构性的影响。固化土冻融结构势(mdσn)的变化规律与动强度类似,随冻融次数的增加和冻结温度的降低而下降。冻融循环次数是影响污泥固化轻质土动力特性和结构性的主要因素,冻结温度为次要因素。
基于数字表面模型的岩体结构面产状获取
宣程强, 章杨松, 许文涛
, doi: 10.16030/j.cnki.issn.1000-3665.202104029
摘要:
针对传统现场接触式测量获取岩体结构面参数效率低、工作量大、结果精确性受人为因素影响等问题,本文结合数字摄影测量技术与运动法(structure from motion,SFM)进行岩体三维数字表面模型重建,并在此基础上建立了岩体结构面自动识别方法。岩体数字表面模型重建步骤主要为岩体影像资料采集,基于尺度不变特征变换(Scale-Invariant Feature Transform,SIFT)算法进行图像特征匹配、稀疏点云构建、点云稠密化以及岩体曲面模型重构。结构面识别方法流程主要为:首先平滑岩体数字表面模型;通过改变搜索半径和角度阈值实现模型平面分割;基于区域生长原理进行结构面搜索;最后基于随机采样一致性拟合结构面得到结构面产状。将该方法应用于甘肃北山地下实验巷道,实现了巷道三维数字表面模型的重建与结构面产状数据获取,最后将识别到的结构面分组表征在模型表面。与人工实地测量方法以及现有的结构面识别软件相比,本文提出的方法具有良好的准确性,可为工程应用提供一定的参考。 针对传统现场接触式测量获取岩体结构面参数效率低、工作量大、结果精确性受人为因素影响等问题,本文结合数字摄影测量技术与运动法(structure from motion,SFM)进行岩体三维数字表面模型重建,并在此基础上建立了岩体结构面自动识别方法。岩体数字表面模型重建步骤主要为岩体影像资料采集,基于尺度不变特征变换(Scale-Invariant Feature Transform,SIFT)算法进行图像特征匹配、稀疏点云构建、点云稠密化以及岩体曲面模型重构。结构面识别方法流程主要为:首先平滑岩体数字表面模型;通过改变搜索半径和角度阈值实现模型平面分割;基于区域生长原理进行结构面搜索;最后基于随机采样一致性拟合结构面得到结构面产状。将该方法应用于甘肃北山地下实验巷道,实现了巷道三维数字表面模型的重建与结构面产状数据获取,最后将识别到的结构面分组表征在模型表面。与人工实地测量方法以及现有的结构面识别软件相比,本文提出的方法具有良好的准确性,可为工程应用提供一定的参考。
基于固结-固化复合技术对温州淤泥加固的试验研究
袁波, 邵吉成, 骆嘉成, 叶宏峰
, doi: 10.16030/j.cnki.issn.1000-3665.202103035
摘要:
为了对软弱淤泥土进行加固,使之满足工程建设所需的一定承载力,在试验基地采用固结-固化复合技术对淤泥进行加固研究。试验时,将一定厚度淤泥分为浅层固化加固层和深层固结加固层;浅层(≤1 m)淤泥采用固化技术进行加固,使之形成高强度硬壳层;对于深层(>1 m)淤泥,采用真空预压技术进行加固,以提高深层淤泥的承载力并控制加固土体的后期沉降量。试验结果表明,当固化剂掺量为0.6%~5.0%时,浅层固化淤泥承载力特征值在109~330 kPa;固结-固化复合技术对土体加固效果突出,分层加固土体整体的承载力特征值在89~230 kPa;浅层淤泥经过固化处理后,土体强度较高,对地表荷载起到了明显的扩散作用,有效地减小了地表荷载在下卧层土体中产生的附加应力;多数试验单元浅层固化土的应力扩散角在19.474°~26.303°之间。 为了对软弱淤泥土进行加固,使之满足工程建设所需的一定承载力,在试验基地采用固结-固化复合技术对淤泥进行加固研究。试验时,将一定厚度淤泥分为浅层固化加固层和深层固结加固层;浅层(≤1 m)淤泥采用固化技术进行加固,使之形成高强度硬壳层;对于深层(>1 m)淤泥,采用真空预压技术进行加固,以提高深层淤泥的承载力并控制加固土体的后期沉降量。试验结果表明,当固化剂掺量为0.6%~5.0%时,浅层固化淤泥承载力特征值在109~330 kPa;固结-固化复合技术对土体加固效果突出,分层加固土体整体的承载力特征值在89~230 kPa;浅层淤泥经过固化处理后,土体强度较高,对地表荷载起到了明显的扩散作用,有效地减小了地表荷载在下卧层土体中产生的附加应力;多数试验单元浅层固化土的应力扩散角在19.474°~26.303°之间。
双排桩双梁组合支护刚度计算的改进与位移分析
宫凤梧, 刘晨, 郭文娟, 李晨光, 封占英, 张银铂
, doi: 10.16030/j.cnki.issn.1000-3665.202103030
摘要:
双排桩支护组合体系作为一种新型悬臂类支护结构,其整体刚度的提升有利于保持基坑边侧的安全稳定。本文依托于张家口万全区某双排桩基坑支护工程案例,以现有双排桩冠梁刚度系数计算方法为基础,引入冠梁与连梁作用效应系数优化改进考虑连梁和冠梁作用的基坑矩形双排桩支护结构横向支撑刚度的计算方法,并对双梁组合支护体系下不同土性对双排桩前后排桩桩身最大横向位移的影响进行探讨。结果显示:(1)在双排桩结构计算中需考虑冠梁与连梁对双排支护桩的共同横向约束作用,并将冠梁与连梁的刚性连接作为一个整体以提高矩形双排桩双梁横向支撑刚度系数。(2)双梁组合支护体系组合刚度对桩顶位移有较大影响,组合刚度为40~50 MN/m下的位移与观测值较为贴近;冠梁计算长度与引入的冠梁与连梁作用效应系数对双梁组合支护体系组合刚度影响较大,计算长度对组合刚度呈负相关,效应系数对组合刚度呈正相关。(3)双梁组合支护体系下双排桩横向支撑刚度受前后排桩竖向与横向位移差影响,前后排桩桩身最大横向位移受土层内摩擦角、黏聚力和土体水平抗力比例系数影响;改变抗拉强度不会影响双排桩桩体位移。在基坑埋深以下及桩底范围内桩身存在位移拐点,拐点处各不同内摩擦角、不同黏聚力条件下位移相等。 双排桩支护组合体系作为一种新型悬臂类支护结构,其整体刚度的提升有利于保持基坑边侧的安全稳定。本文依托于张家口万全区某双排桩基坑支护工程案例,以现有双排桩冠梁刚度系数计算方法为基础,引入冠梁与连梁作用效应系数优化改进考虑连梁和冠梁作用的基坑矩形双排桩支护结构横向支撑刚度的计算方法,并对双梁组合支护体系下不同土性对双排桩前后排桩桩身最大横向位移的影响进行探讨。结果显示:(1)在双排桩结构计算中需考虑冠梁与连梁对双排支护桩的共同横向约束作用,并将冠梁与连梁的刚性连接作为一个整体以提高矩形双排桩双梁横向支撑刚度系数。(2)双梁组合支护体系组合刚度对桩顶位移有较大影响,组合刚度为40~50 MN/m下的位移与观测值较为贴近;冠梁计算长度与引入的冠梁与连梁作用效应系数对双梁组合支护体系组合刚度影响较大,计算长度对组合刚度呈负相关,效应系数对组合刚度呈正相关。(3)双梁组合支护体系下双排桩横向支撑刚度受前后排桩竖向与横向位移差影响,前后排桩桩身最大横向位移受土层内摩擦角、黏聚力和土体水平抗力比例系数影响;改变抗拉强度不会影响双排桩桩体位移。在基坑埋深以下及桩底范围内桩身存在位移拐点,拐点处各不同内摩擦角、不同黏聚力条件下位移相等。
系统维度对变密度溶质运移的影响研究
叶逾, 蔡芳敏, 谢一凡, 井淼, 鲁春辉
, doi: 10.16030/j.cnki.issn.1000-3665.202103015
摘要:
关于地下水变密度流与溶质运移的研究通常局限于二维系统中,三维系统中的研究非常有限。然而,仍有零星研究表明,三维系统中的对流过程与二维系统不同。本文通过数值模拟,系统研究了二维和三维系统中的变密度溶质运移过程,利用舍伍德数、空间矩和稀释指数量化了系统的不稳定性和溶质扩散、稀释程度。结果表明,二维系统产生分散指流,而三维系统中因扩散程度的增强使得分散指流的产生受到抑制,但是,三维系统的不稳定性比二维系统更强,对流入渗更快,与指流的产生与否并无直接关联。另外,三维系统中溶质的稀释程度大于二维系统,但溶质在二维系统中将更快接近稀释的最大值,传统二阶中心矩可能会造成自由对流和不稳定性溶质运移过程中溶质扩散和稀释程度的错误估计。 关于地下水变密度流与溶质运移的研究通常局限于二维系统中,三维系统中的研究非常有限。然而,仍有零星研究表明,三维系统中的对流过程与二维系统不同。本文通过数值模拟,系统研究了二维和三维系统中的变密度溶质运移过程,利用舍伍德数、空间矩和稀释指数量化了系统的不稳定性和溶质扩散、稀释程度。结果表明,二维系统产生分散指流,而三维系统中因扩散程度的增强使得分散指流的产生受到抑制,但是,三维系统的不稳定性比二维系统更强,对流入渗更快,与指流的产生与否并无直接关联。另外,三维系统中溶质的稀释程度大于二维系统,但溶质在二维系统中将更快接近稀释的最大值,传统二阶中心矩可能会造成自由对流和不稳定性溶质运移过程中溶质扩散和稀释程度的错误估计。
基于时空维度耦合的地质灾害发育程度评价研究
曲雪妍, 李媛, 房浩, 杨旭东, 谢振桦, 尹春荣, 张艳玲, 佟彬
, doi: 10.16030/j.cnki.issn.1000-3665.201912061
摘要:
在现有地质灾害综合评价指标体系中,地质灾害发育程度评价因子只包括反映空间维度的点数量、点密度、体积数、体积密度等,不包括反映时间维度的评价因子。文章引入反映时间维度的受灾年份因子充实之。通过对地质灾害体积取中位数对数法、对受灾年份采用概率密度法进行处理和采用熵权法确定评价因子权重来建立时空耦合评价模型;采用斜率法将评价结果划分为高、较高、中、低四级,以反映已发生灾害的发育状况。选取中国(不包括港、澳、台)2011—2020年已发生的地质灾害点,按此方法对地质灾害发育程度进行评价与区划,结果表明:地质灾害高发育区共323个县,涉及19个省份,主要分布在湖南大部、湖北西部和东部、江西东南部、广东东部、福建中西部、四川东部、重庆北部、甘肃南部等地。较高发育区共566个县,涉及25个省份,主要分布在四川大部、云南西部和南部、湖南东南部、湖北北部、江西大部、广东大部、广西大部、福建大部、浙江大部、贵州西南部、陕西南部以及吉林东部。中发育区623个县,涉及30个省份,主要包括新疆大部、西藏大部、甘肃大部、河北大部、山西大部、河南大部。低发育区共1336个县,涉及30个省份,主要包括青海、西藏大部、甘肃北部、内蒙、黑龙江、吉林西部、辽宁西部、河北东部及南部、山东、河南东部、安徽中北部、江苏等地。本次评价结果与国家防灾部署、灾害发育分布具有较高的吻合度。野外实地调研证明,引入受灾年份后的评价结果较引入前更科学、合理。 在现有地质灾害综合评价指标体系中,地质灾害发育程度评价因子只包括反映空间维度的点数量、点密度、体积数、体积密度等,不包括反映时间维度的评价因子。文章引入反映时间维度的受灾年份因子充实之。通过对地质灾害体积取中位数对数法、对受灾年份采用概率密度法进行处理和采用熵权法确定评价因子权重来建立时空耦合评价模型;采用斜率法将评价结果划分为高、较高、中、低四级,以反映已发生灾害的发育状况。选取中国(不包括港、澳、台)2011—2020年已发生的地质灾害点,按此方法对地质灾害发育程度进行评价与区划,结果表明:地质灾害高发育区共323个县,涉及19个省份,主要分布在湖南大部、湖北西部和东部、江西东南部、广东东部、福建中西部、四川东部、重庆北部、甘肃南部等地。较高发育区共566个县,涉及25个省份,主要分布在四川大部、云南西部和南部、湖南东南部、湖北北部、江西大部、广东大部、广西大部、福建大部、浙江大部、贵州西南部、陕西南部以及吉林东部。中发育区623个县,涉及30个省份,主要包括新疆大部、西藏大部、甘肃大部、河北大部、山西大部、河南大部。低发育区共1336个县,涉及30个省份,主要包括青海、西藏大部、甘肃北部、内蒙、黑龙江、吉林西部、辽宁西部、河北东部及南部、山东、河南东部、安徽中北部、江苏等地。本次评价结果与国家防灾部署、灾害发育分布具有较高的吻合度。野外实地调研证明,引入受灾年份后的评价结果较引入前更科学、合理。
近水平层状坝基岩体渗透结构及其工程意义
王俊智, 李清波, 王贵军, 张一
, doi: 10.16030/j.cnki.issn.1000-3665.202105024
摘要:
查明坝基岩体的渗透结构对于大坝工程防渗方案的制定具有重要意义。以黄河古贤水利枢纽工程近水平分布的红层坝基岩体为研究对象,使用“连续、高分辨率”压水试验数据处理与统计方法,分析了岩体透水率与高程、岩性及顺层剪切带的关系,阐明了坝基岩体的渗透结构,提出了工程防渗建议。研究发现,坝基岩体渗透结构表现出结构性和随机性双重特征。结构性表现为岩体透水率随高程的降低而减小;随机性表现为岩体透水率会因岩性和剪切构造等随机变化。对于河床及滩地坝基岩体,风化卸荷带底高程在450 m附近,1 Lu关键透水率指标界线的底高程位于340 m附近。对于岸坡坝基岩体,3 Lu关键透水率指标界线的底高程位于560~580 m。在坝基勘察深度范围内,350~360 m高程段连续分布具有一定厚度的软岩类岩体,剪切构造不发育且岩体透水微弱,可视为相对隔水层,对坝基防渗帷幕优化设计具有指示意义。文章提出的统计、分析方法具有一定普适性,可为类似工程所借鉴。 查明坝基岩体的渗透结构对于大坝工程防渗方案的制定具有重要意义。以黄河古贤水利枢纽工程近水平分布的红层坝基岩体为研究对象,使用“连续、高分辨率”压水试验数据处理与统计方法,分析了岩体透水率与高程、岩性及顺层剪切带的关系,阐明了坝基岩体的渗透结构,提出了工程防渗建议。研究发现,坝基岩体渗透结构表现出结构性和随机性双重特征。结构性表现为岩体透水率随高程的降低而减小;随机性表现为岩体透水率会因岩性和剪切构造等随机变化。对于河床及滩地坝基岩体,风化卸荷带底高程在450 m附近,1 Lu关键透水率指标界线的底高程位于340 m附近。对于岸坡坝基岩体,3 Lu关键透水率指标界线的底高程位于560~580 m。在坝基勘察深度范围内,350~360 m高程段连续分布具有一定厚度的软岩类岩体,剪切构造不发育且岩体透水微弱,可视为相对隔水层,对坝基防渗帷幕优化设计具有指示意义。文章提出的统计、分析方法具有一定普适性,可为类似工程所借鉴。
循环荷载下水泥土桩复合体动力参数试验研究
叶观宝, 秦粮凯, 张振, 郑文强, 陈勇
, doi: 10.16030/j.cnki.issn.1000-3665.202103037
摘要:
水泥土桩被广泛应用于软土路基加固工程中。然而,人们对水泥桩与桩间土形成的加固体的动力特性尚缺乏认识,无法合理评价水泥土桩复合地基的长期性能。基于此,本文开展水泥土桩复合体大型动三轴试验,研究围压、静偏应力、置换率及分级加卸载路径对其动力参数的影响,并分析了动力参数的波动性。试验结果表明:随着静偏应力增加,复合单元体的动弹模量减小,阻尼比增大,临界动应力比减小。随着置换率增加,动弹模量略有增加,阻尼比略有减小。逐级卸载造成复合单元体的动力参数劣化。阻尼比具有较强的波动性,复合单元体阻尼比的变异系数是动弹模量的2.8~7.0倍。相比于软土,复合体动弹模量提高了2~6倍,静偏应力越大,提高系数越大。 水泥土桩被广泛应用于软土路基加固工程中。然而,人们对水泥桩与桩间土形成的加固体的动力特性尚缺乏认识,无法合理评价水泥土桩复合地基的长期性能。基于此,本文开展水泥土桩复合体大型动三轴试验,研究围压、静偏应力、置换率及分级加卸载路径对其动力参数的影响,并分析了动力参数的波动性。试验结果表明:随着静偏应力增加,复合单元体的动弹模量减小,阻尼比增大,临界动应力比减小。随着置换率增加,动弹模量略有增加,阻尼比略有减小。逐级卸载造成复合单元体的动力参数劣化。阻尼比具有较强的波动性,复合单元体阻尼比的变异系数是动弹模量的2.8~7.0倍。相比于软土,复合体动弹模量提高了2~6倍,静偏应力越大,提高系数越大。
缺水基岩山区蓄水构造类型划分及其属性分析
张福存, 王新峰, 李伟, 程旭学, 韩双宝, 潘晓东, 龚磊
, doi: 10.16030/j.cnki.issn.1000-3665.202107042
摘要:
近十多年来,服务脱贫攻坚的基岩山区找水打井工作充分运用并丰富了蓄水构造理论,但缺乏对其类型的统一划分和属性的深入分析。以蓄水构造为理论指导,考虑地貌汇水、岩性赋水、构造控水等控制因素,以构造为主控因素,综合地形地貌、地层岩性和植被发育及开发利用实际价值等因素,系统总结划分出缺水基岩山区水平岩层、单斜、褶皱、断裂、接触、风化壳、复合、洞藏等8类20型35式蓄水构造,并借鉴、结合新构造控水理论和地下水系统理论分析总结不同尺度、不同类型蓄水构造的共性特征,从蓄水构造空间组合特征、边界水文地质性质、水动力条件、尺度效应、控水共生特征、开发利用实际价值等方面,提出了岩性与地质构造的相关性、透水与隔水的相对性、汇水与蓄水的平衡性、尺度与系统的统一性、形成与改造的继承性和规模与目的的匹配性等蓄水构造属性。研究成果细化完善了缺水基岩山区地下水蓄水构造类型划分,总结提炼了蓄水构造属性特征,深化了对蓄水构造概念内涵的理解与认识。既进一步推动了基岩山区水文地质学的研究,也可为缺水区找水定井提供参考,并指导缺水区水资源可持续利用。 近十多年来,服务脱贫攻坚的基岩山区找水打井工作充分运用并丰富了蓄水构造理论,但缺乏对其类型的统一划分和属性的深入分析。以蓄水构造为理论指导,考虑地貌汇水、岩性赋水、构造控水等控制因素,以构造为主控因素,综合地形地貌、地层岩性和植被发育及开发利用实际价值等因素,系统总结划分出缺水基岩山区水平岩层、单斜、褶皱、断裂、接触、风化壳、复合、洞藏等8类20型35式蓄水构造,并借鉴、结合新构造控水理论和地下水系统理论分析总结不同尺度、不同类型蓄水构造的共性特征,从蓄水构造空间组合特征、边界水文地质性质、水动力条件、尺度效应、控水共生特征、开发利用实际价值等方面,提出了岩性与地质构造的相关性、透水与隔水的相对性、汇水与蓄水的平衡性、尺度与系统的统一性、形成与改造的继承性和规模与目的的匹配性等蓄水构造属性。研究成果细化完善了缺水基岩山区地下水蓄水构造类型划分,总结提炼了蓄水构造属性特征,深化了对蓄水构造概念内涵的理解与认识。既进一步推动了基岩山区水文地质学的研究,也可为缺水区找水定井提供参考,并指导缺水区水资源可持续利用。
随机网络裂隙中DNAPLs运移规律研究
常兴, 骆乾坤, 邓亚平, 马雷, 钱家忠
, doi: 10.16030/j.cnki.issn.1000-3665.202103055
摘要:
重非水相流体(DNAPLs)在裂隙介质中运移和分布受多种因素影响,包括DNAPLs本身的物理化学性质、裂隙的几何性质以及泄漏条件等。以往研究多基于单个裂隙,随机网络裂隙中DNAPLs运移规律的研究较少。文章基于蒙特卡罗方法生成随机网络裂隙,运用像素扫描识别并输出裂隙的坐标和宽度,然后采用PetraSim模拟四氯乙烯(PCE)在随机网络裂隙中的运移,主要探讨裂隙宽度空间变异性和泄漏条件(包括泄漏速率和泄漏位置)对DNAPLs运移规律的影响。数值模拟结果表明裂隙空间变异性的不同,同样影响DNAPLs的运移路径和空间展布。随着网络裂隙中裂隙宽度空间变异性增大,出现优势通道,DNAPLs运移的速率加快,DNAPLs的质心位置和饱和度空间分布发生明显变化;泄露速率影响DNAPLs的运移范围和空间展布,泄漏速率越高,DNAPLs运移速率越快,模型底部蓄积的DNAPLs的饱和度越大,DNAPLs的空间展布也越大;同一网络裂隙中,泄漏位置的不同,同样会导致DNAPLs的运移路径及分布范围不同,不同的泄露位置重力方向裂隙空间变异性不同,导致DNAPLs运移路径和空间展布各不相同。结果可丰富裂隙介质中DNAPLs运移机理研究,为裂隙介质中DNAPLs污染修复提供模型参考。 重非水相流体(DNAPLs)在裂隙介质中运移和分布受多种因素影响,包括DNAPLs本身的物理化学性质、裂隙的几何性质以及泄漏条件等。以往研究多基于单个裂隙,随机网络裂隙中DNAPLs运移规律的研究较少。文章基于蒙特卡罗方法生成随机网络裂隙,运用像素扫描识别并输出裂隙的坐标和宽度,然后采用PetraSim模拟四氯乙烯(PCE)在随机网络裂隙中的运移,主要探讨裂隙宽度空间变异性和泄漏条件(包括泄漏速率和泄漏位置)对DNAPLs运移规律的影响。数值模拟结果表明裂隙空间变异性的不同,同样影响DNAPLs的运移路径和空间展布。随着网络裂隙中裂隙宽度空间变异性增大,出现优势通道,DNAPLs运移的速率加快,DNAPLs的质心位置和饱和度空间分布发生明显变化;泄露速率影响DNAPLs的运移范围和空间展布,泄漏速率越高,DNAPLs运移速率越快,模型底部蓄积的DNAPLs的饱和度越大,DNAPLs的空间展布也越大;同一网络裂隙中,泄漏位置的不同,同样会导致DNAPLs的运移路径及分布范围不同,不同的泄露位置重力方向裂隙空间变异性不同,导致DNAPLs运移路径和空间展布各不相同。结果可丰富裂隙介质中DNAPLs运移机理研究,为裂隙介质中DNAPLs污染修复提供模型参考。
地下水机器学习方法研究—水位监测数据驱动的区域补排边界识别
齐永强, 李文鹏, 郑跃军, 李慧, 王成见
, doi: 10.16030/j.cnki.issn.1000-3665.202104030
摘要:
随着我国地下水监测工作的高速发展,高频率高密度水位监测数据的出现催生了对其进行深入信息挖掘的需求。在传统地下水模型研究中,地下水位监测值常位于模型构建过程的下游,当水位监测的时空密度逐渐增大时,新增信息无法有效传导至模型的规划阶段并指导概念模型的修订。文章提出了一种地下水系统补排边界的识别方法,在不建立地下水数值模型的前提下,以监测井空间位置为节点,按照德劳内原则建立三角网格。在此网格系统中,首先定义一个水力梯度变换函数gradF,以求取网格中任意位置的水力梯度;借鉴机器学习领域的优化算法,使用水力梯度场驱动含水层中随机分布质点的运行轨迹,并以此推断和识别区域内地下水补给和排泄边界。在环境地学计算平台EnviFusion-CGS中实现,并构建了详细工作流程。以山东省青岛市大沽河中下游含水层作为示范区,对含水系统的补给区和排泄区的空间分布及其动态变化进行了分析,取得了良好效果。本研究为构建和修订已有的含水层概念模型提供了新思路。 随着我国地下水监测工作的高速发展,高频率高密度水位监测数据的出现催生了对其进行深入信息挖掘的需求。在传统地下水模型研究中,地下水位监测值常位于模型构建过程的下游,当水位监测的时空密度逐渐增大时,新增信息无法有效传导至模型的规划阶段并指导概念模型的修订。文章提出了一种地下水系统补排边界的识别方法,在不建立地下水数值模型的前提下,以监测井空间位置为节点,按照德劳内原则建立三角网格。在此网格系统中,首先定义一个水力梯度变换函数gradF,以求取网格中任意位置的水力梯度;借鉴机器学习领域的优化算法,使用水力梯度场驱动含水层中随机分布质点的运行轨迹,并以此推断和识别区域内地下水补给和排泄边界。在环境地学计算平台EnviFusion-CGS中实现,并构建了详细工作流程。以山东省青岛市大沽河中下游含水层作为示范区,对含水系统的补给区和排泄区的空间分布及其动态变化进行了分析,取得了良好效果。本研究为构建和修订已有的含水层概念模型提供了新思路。
利用示踪试验时间-浓度曲线分析岩溶管道结构特征
陈亚洲, 董维红
, doi: 10.16030/j.cnki.issn.1000-3665201909067
摘要:
岩溶管道结构特征对岩溶水资源的合理取用、保护及地下工程安全施工具有重要影响。目前在利用示踪试验曲线分析岩溶管道结构特征时,存在曲线叠加形态、钝锋形态、曲线不规则上升和下降与多条岩溶管道连接关系、地下湖所处位置、地下水状态难以准确分析的问题。文章通过地下水溶质运移理论,推导出岩溶管道流溶质运移模型,根据模型绘制三维溶质运移形态,根据三维溶质运移形态推导出理论时间-浓度曲线形态,利用理论时间-浓度曲线形态结合水力学相关知识,对岩溶管道连接关系、地下湖位置及岩溶管道流形态变化对应的时间-浓度曲线进行解释。得出以下结论:(1)曲线出峰个数对应岩溶管道条数,由于管道径流长度及流速存在差异,双管道并联曲线存在三种模型;(2)单管道曲线下降梯度个数对应地下湖个数,多管道需结合管道个数和地下湖位置具体分析下降梯度和地下湖个数关系,根据地下湖位置将双管道并联岩溶管道划分四种类型;(3)曲线形态极速变化标志着管道流发生表流和承压流的相互转化。 岩溶管道结构特征对岩溶水资源的合理取用、保护及地下工程安全施工具有重要影响。目前在利用示踪试验曲线分析岩溶管道结构特征时,存在曲线叠加形态、钝锋形态、曲线不规则上升和下降与多条岩溶管道连接关系、地下湖所处位置、地下水状态难以准确分析的问题。文章通过地下水溶质运移理论,推导出岩溶管道流溶质运移模型,根据模型绘制三维溶质运移形态,根据三维溶质运移形态推导出理论时间-浓度曲线形态,利用理论时间-浓度曲线形态结合水力学相关知识,对岩溶管道连接关系、地下湖位置及岩溶管道流形态变化对应的时间-浓度曲线进行解释。得出以下结论:(1)曲线出峰个数对应岩溶管道条数,由于管道径流长度及流速存在差异,双管道并联曲线存在三种模型;(2)单管道曲线下降梯度个数对应地下湖个数,多管道需结合管道个数和地下湖位置具体分析下降梯度和地下湖个数关系,根据地下湖位置将双管道并联岩溶管道划分四种类型;(3)曲线形态极速变化标志着管道流发生表流和承压流的相互转化。
共和盆地干热岩体人工裂隙带结构的控热机理与产能优化
陈炫沂, 姜振蛟, 徐含英, 冯波
, doi: 10.16030/j.cnki.issn.1000-3665.202104058
摘要:
人工压裂是获取干热岩型地热资源的关键环节,压裂后的人工裂隙带结构对开采条件下水热传递过程具有重要控制作用。结合我国共和盆地干热岩储层地质条件,采用数值模拟方法着重分析干热岩不同产状人工裂隙带的渗透率与宽度对热储中水热传递过程的影响机理,明确不同人工裂隙结构条件下水热产出能力,进而优化井间距。结果表明:当人工裂隙带渗透率较小时(小于5 D),裂隙带规模越大,开采井温度越高;当渗透率较大时(大于10 D),在水平裂隙带中,随着裂隙带规模的增加,由于注入冷水的快速扩散导致整体低温区域增加,开采井温度反而降低。在水平裂隙带中注入冷水主要为水平向流动,随着渗透率的增加,开采井温度更易受注入冷水的影响而降低;但在垂直裂隙带及倾斜裂隙带中,随着渗透率的增加,垂向自由对流增强,注入冷水更易于向储层底部高温区域流动,经加热后到达开采井,使得开采温度提升。综合比较,同一井间距条件下,低渗水平裂隙带以及高渗垂直裂隙带的产热能力较其他裂隙带更强。 人工压裂是获取干热岩型地热资源的关键环节,压裂后的人工裂隙带结构对开采条件下水热传递过程具有重要控制作用。结合我国共和盆地干热岩储层地质条件,采用数值模拟方法着重分析干热岩不同产状人工裂隙带的渗透率与宽度对热储中水热传递过程的影响机理,明确不同人工裂隙结构条件下水热产出能力,进而优化井间距。结果表明:当人工裂隙带渗透率较小时(小于5 D),裂隙带规模越大,开采井温度越高;当渗透率较大时(大于10 D),在水平裂隙带中,随着裂隙带规模的增加,由于注入冷水的快速扩散导致整体低温区域增加,开采井温度反而降低。在水平裂隙带中注入冷水主要为水平向流动,随着渗透率的增加,开采井温度更易受注入冷水的影响而降低;但在垂直裂隙带及倾斜裂隙带中,随着渗透率的增加,垂向自由对流增强,注入冷水更易于向储层底部高温区域流动,经加热后到达开采井,使得开采温度提升。综合比较,同一井间距条件下,低渗水平裂隙带以及高渗垂直裂隙带的产热能力较其他裂隙带更强。
通过永定河补水渗漏补给地下水的适宜性评估研究
曹欣怡, 翟远征, 李木子, 潘成忠, 郑富新, 芦红, 夏雪莲, 滕彦国, 王金生
, doi: 10.16030/j.cnki.issn.1000-3665.202107069
摘要:
受气候变化和人类活动等因素的综合影响,地下水超采和含水层水量亏空已成为备受关注的全球性问题。为了弥补含水层水量亏空和促进地下水资源涵养,已探索出多种人为干预地下水补给的措施,其中通过河道开展地下水人工补给由于具有明显优势而受到重视。已有实践表明,受渗漏补给潜力和包气带调蓄能力等的限制,并不是所有河流或河段都适宜开展地下水人工补给工作。目前尚缺乏比较为大家广泛接受的适宜性评估方法。本研究以永定河生态补水为契机,以北京平原段河道和下伏含水层为研究区开展了案例研究,并用实测数据对研究结果进行了检验。应用指标体系法建立的适宜性评估模型(LMBGITSC模型)包括河床土地利用类型、河床介质类型、河床宽度、河床地形坡度、包气带介质类型、包气带厚度、包气带给水度、包气带水平渗透性等8个指标。案例研究结果表明,沿河流流向,通过河道渗漏补给地下水的适宜性由好变差(防渗河段除外)。该规律性主要受渗漏补给潜力和包气带调蓄能力控制,因此适宜性也呈现出“阶梯式”演变规律。检验结果表明所建的评估方法具有较好的适用性。该方法具有较好的可移植性,因此也可为其他同类地区提供参考。 受气候变化和人类活动等因素的综合影响,地下水超采和含水层水量亏空已成为备受关注的全球性问题。为了弥补含水层水量亏空和促进地下水资源涵养,已探索出多种人为干预地下水补给的措施,其中通过河道开展地下水人工补给由于具有明显优势而受到重视。已有实践表明,受渗漏补给潜力和包气带调蓄能力等的限制,并不是所有河流或河段都适宜开展地下水人工补给工作。目前尚缺乏比较为大家广泛接受的适宜性评估方法。本研究以永定河生态补水为契机,以北京平原段河道和下伏含水层为研究区开展了案例研究,并用实测数据对研究结果进行了检验。应用指标体系法建立的适宜性评估模型(LMBGITSC模型)包括河床土地利用类型、河床介质类型、河床宽度、河床地形坡度、包气带介质类型、包气带厚度、包气带给水度、包气带水平渗透性等8个指标。案例研究结果表明,沿河流流向,通过河道渗漏补给地下水的适宜性由好变差(防渗河段除外)。该规律性主要受渗漏补给潜力和包气带调蓄能力控制,因此适宜性也呈现出“阶梯式”演变规律。检验结果表明所建的评估方法具有较好的适用性。该方法具有较好的可移植性,因此也可为其他同类地区提供参考。
基于可靠度方法的全基质吸力段土-水特征模型研究
刘青灵, 简文彬, 许旭堂, 聂闻
, doi: 10.16030/j.cnki.issn.1000-3665.202104033
摘要:
土-水特征曲线是非饱和土入渗和流-固耦合分析的基础,目前该模型的研究成果由于试验方法和设备的局限,仅能够得到局部基质吸力段(低吸力或高吸力)的土-水特征曲线;而如何根据局部试验数据集推断与试验结果相吻合的全基质吸力段范围的土-水特征曲线是一个难题。对东南沿海广泛分布的原状残积土进行了不同基质吸力段的土-水特征试验。根据残积土低吸力段和高吸力段的试验数据,运用可靠度分析方法,从土壤孔隙微元破裂失稳的角度,建立了土壤全基质吸力段的土-水特征模型。根据UNSODA2.0数据库提取砂土、粉土、黏土等3个类型土样本(52个)的试验数据,讨论了模型参数意义和不同土类的模型适用性,进行了不同土-水特征模型的比较分析。对比验证研究表明,形状参数σ和中值比例参数\begin{document}$ \mu $\end{document}使提出的模型具有很强的适应性,并且具有依据局部试验数据(低或高基质吸力段)快速获得全基质吸力段土壤土-水特征曲线的能力。与其它模型相比,模型计算过程适应性强,模型决定系数R2>0.98。研究结果对快速获取不同类型土壤的全基质吸力段土-水特征曲线具有参考价值。 土-水特征曲线是非饱和土入渗和流-固耦合分析的基础,目前该模型的研究成果由于试验方法和设备的局限,仅能够得到局部基质吸力段(低吸力或高吸力)的土-水特征曲线;而如何根据局部试验数据集推断与试验结果相吻合的全基质吸力段范围的土-水特征曲线是一个难题。对东南沿海广泛分布的原状残积土进行了不同基质吸力段的土-水特征试验。根据残积土低吸力段和高吸力段的试验数据,运用可靠度分析方法,从土壤孔隙微元破裂失稳的角度,建立了土壤全基质吸力段的土-水特征模型。根据UNSODA2.0数据库提取砂土、粉土、黏土等3个类型土样本(52个)的试验数据,讨论了模型参数意义和不同土类的模型适用性,进行了不同土-水特征模型的比较分析。对比验证研究表明,形状参数σ和中值比例参数$ \mu $使提出的模型具有很强的适应性,并且具有依据局部试验数据(低或高基质吸力段)快速获得全基质吸力段土壤土-水特征曲线的能力。与其它模型相比,模型计算过程适应性强,模型决定系数R2>0.98。研究结果对快速获取不同类型土壤的全基质吸力段土-水特征曲线具有参考价值。
长江干流水化学成因与风化过程CO2消耗通量解析
张鸿, 周权平, 姜月华, 金阳, 杨国强, 顾轩, 梅世嘉, 王晓龙
, doi: 10.16030/j.cnki.issn.1000-3665.202106027
摘要:
长江流域面积巨大,岩性多变,加之三峡大坝等重大水利工程的影响,干流河水的水化学成因存在较大争议。此外,以往研究中流域矿物风化过程的碳汇通量估算一般基于阳离子来源分析,但该算法通常涉及多种矿物端元的参数选取,结果具有不确定性。本次研究对长江干流水化学的时空演变进行了整体分析,并基于上游河水样品\begin{document}$\rm{HCO}_3^{-}$\end{document}含量的校正与计算,提出了一种计算矿物风化过程碳汇通量的新方法。研究结果表明,蒸发盐溶解、循环盐作用、矿物风化及硫酸盐溶解是控制长江干流河水离子组成的主要水文地球化学作用,而人类活动主要影响了离海距离3 000 km以内河水\begin{document}$\rm{NO}_3^{-}$\end{document}含量;长江上游干流硅酸盐风化消耗CO2速率为1.16×105 mol/(km2·a),碳酸盐风化消耗CO2速率为4.75×105 mol/(km2·a)。本研究有助于加深对长江干流主要水文地球化学作用的认识,丰富和完善碳循环研究理论。 长江流域面积巨大,岩性多变,加之三峡大坝等重大水利工程的影响,干流河水的水化学成因存在较大争议。此外,以往研究中流域矿物风化过程的碳汇通量估算一般基于阳离子来源分析,但该算法通常涉及多种矿物端元的参数选取,结果具有不确定性。本次研究对长江干流水化学的时空演变进行了整体分析,并基于上游河水样品$\rm{HCO}_3^{-}$含量的校正与计算,提出了一种计算矿物风化过程碳汇通量的新方法。研究结果表明,蒸发盐溶解、循环盐作用、矿物风化及硫酸盐溶解是控制长江干流河水离子组成的主要水文地球化学作用,而人类活动主要影响了离海距离3 000 km以内河水$\rm{NO}_3^{-}$含量;长江上游干流硅酸盐风化消耗CO2速率为1.16×105 mol/(km2·a),碳酸盐风化消耗CO2速率为4.75×105 mol/(km2·a)。本研究有助于加深对长江干流主要水文地球化学作用的认识,丰富和完善碳循环研究理论。
吉木萨尔超排液氢氧同位素特征及地下水含量计算
贾俊, 王亮, 孟云涛
, doi: 10.16030/j.cnki.issn.1000-3665.202103045
摘要:
压裂返排液对环境安全和人类健康构成巨大挑战,也为地质和工程应用研究提供了宝贵资料。新疆吉木萨尔芦草沟组压裂返排液存在超排现象,为了揭示超排液来源与组成,获取原始地层卤水地球化学信息,在缺乏地下水采样情况下,系统采集研究区内两口水平井压裂液、压裂返排液和地表水样品,测定并分析其氢、氧同位素组成及随返排时间动态变异特征。在对水岩作用可能引起的同位素交换分析厘定基础上,分别测算地下水氧同位素组成及返排液中地下水含量。结果表明:相对于富D贫18O的地表水和压裂液,返排液样点明显偏离该地区大气降水线,18O同位素发生了正漂移,D同位素呈现负漂移,其水体来源发生了变化。依据返排液δ18O实测数据随返排时间的收敛模型,计算地下水δ18O为−6.902‰;返排第60天,返排液中的地下水含量分别为84%、81%,尚有43 283.4,39 150.3 m3压裂液未返出,即随着压裂改造缝不断延伸,返排时间的延长,前期注入的压裂液绝大部分仍滞留于地层中,超排液来源为持续混入的地下水。 压裂返排液对环境安全和人类健康构成巨大挑战,也为地质和工程应用研究提供了宝贵资料。新疆吉木萨尔芦草沟组压裂返排液存在超排现象,为了揭示超排液来源与组成,获取原始地层卤水地球化学信息,在缺乏地下水采样情况下,系统采集研究区内两口水平井压裂液、压裂返排液和地表水样品,测定并分析其氢、氧同位素组成及随返排时间动态变异特征。在对水岩作用可能引起的同位素交换分析厘定基础上,分别测算地下水氧同位素组成及返排液中地下水含量。结果表明:相对于富D贫18O的地表水和压裂液,返排液样点明显偏离该地区大气降水线,18O同位素发生了正漂移,D同位素呈现负漂移,其水体来源发生了变化。依据返排液δ18O实测数据随返排时间的收敛模型,计算地下水δ18O为−6.902‰;返排第60天,返排液中的地下水含量分别为84%、81%,尚有43 283.4,39 150.3 m3压裂液未返出,即随着压裂改造缝不断延伸,返排时间的延长,前期注入的压裂液绝大部分仍滞留于地层中,超排液来源为持续混入的地下水。