吉林大学学报(地球科学版) ›› 2025, Vol. 55 ›› Issue (3): 893-905.doi: 10.13278/j.cnki.jjuese.20230299

• 地质工程与环境工程 • 上一篇    下一篇

考虑边界层流体黏度变化的黏性土低速非达西渗流数学模型

王福刚1,2,管小桐1,2,何庆成3,程辉1,2,杨国华1,2,程中乐1,2,王耀辉1,2   

  1. 1.地下水资源与环境教育部重点实验室(吉林大学),长春130021
    2.吉林省水资源与水环境重点实验室,长春130021
    3.中国地质科学院,北京100037

  • 出版日期:2025-05-26 发布日期:2025-06-06
  • 作者简介:王福刚(1975-),男,教授,博士生导师,主要从事地下水资源形成演化理论和数值模拟方面的研究,E-mail:wangfugang@jlu.edu.cn
  • 基金资助:
    国家自然科学基金项目(42072280,U2244215,41172205);国家重点研发计划项目(2019YFC0409103) 

 A Mathematical Model of Low-Velocity Non-Darcy Flow in Clayey Soil Considering Change of Viscosity in Boundary Layer Fluid

Wang Fugang1,2, Guan Xiaotong1,2, He Qingcheng3, Cheng Hui1,2, Yang Guohua1,2, Cheng Zhongle1,2, Wang Yaohui1,2#br#

#br#
  

  1. 1. Key Laboratory of Groundwater Resources and Environment (Jilin University), Ministry of Education, Changchun 130021, China
    2. Jilin Provincial Key Laboratory of Water Resources and Water Environment, Changchun 130021, China
    3. Chinese Academy of Geological Sciences, Beijing 100037, China
  • Online:2025-05-26 Published:2025-06-06
  • Supported by:
    Supported by the National Natural Science Foundation of China (42072280,U2244215,41172250) and the National Key Research and Development Program of China (2019YFC0409103)

摘要: 边界层流体及其黏度变化对低渗介质中流体的流动具有重要影响。现有低速非达西流模型通常忽略流体黏度的变化,或简单地将流体黏度处理为一个与孔隙半径有关的函数;而边界层流体黏度通常与接触角有关。本文基于毛细管束模型和分形理论,考虑接触角对边界层流体黏度的影响,构建了一种新的低速非达西渗流模型,通过实验对模型进行了可靠性验证,并对模型中的关键参数进行了敏感性分析。结果表明:接触角对边界层流体黏度有显著影响,进而影响着流体的流动行为。考虑边界层流体黏度变化后,渗流速度下降,本文实验中,考虑流体黏度变化与不考虑流体黏度变化的情况相比,渗流速度的下降率可达57.1%。渗流速度随孔隙分布分形维数Df的增大而下降,随接触角的增大而增大。本文中,在同一压力梯度下,当Df从1.2增至1.8,每增加0.2,渗流速度分别减小18.4%、23.3%、29.1%;当接触角从0°增至80°,每增加20°,渗流速度分别增长10.9%、12.3%、14.0%、16.3%。随压力梯度增大,边界层流体参与流动,有效渗流断面和流体有效黏度增大,流体黏度增大对其流动产生阻滞效应,但作用程度小于渗流断面增大对流速的增益效应,因此,总体上呈现出渗流速度随压力梯度增大而增大。


关键词: 地下水, 黏性土, 低速非达西流, 边界层, 流体黏度, 分形理论

Abstract:  The boundary layer fluid and its viscosity changeshave a significant impact on the flow of fluids in low permeability media. Existing low-velocity non-Darcy flow models usually ignore the change of fluid viscosity, or simply treat the fluid viscosity as a function of pore radius. The viscosity of boundary layer fluids is usually related to contact angle. In this paper, based on capillary model and fractal theory, a new low-velocity non-Darcy flow model is proposed considering the influence of contact angle on the viscosity of boundary layer fluid. The reliability of the model was verified through experiments and the sensitivity analysis of key parameters was conducted. The results show that the contact angle has a significant impact on the viscosity of the boundary layer fluid, which in turn affects the flow behavior of the fluid. After considering the change of fluid viscosity in the boundary layer, the flow velocity decreases. In this experiment, compared with ignoring the change of fluid viscosity, the decrease in flow velocity can reach 57.1%. The flow velocity decreases with the increase of fractal dimension of pore distribution (Df) and increases with the increase of contact angle. In this paper, under the same pressure gradient, when Df increases from 1.2 to 1.8, for every 0.2 increase, the flow velocity decreases by 18.4%, 23.3%, and 29.1%, respectively; When the contact angle increases from 0° to 80°, for every 20° increase, the flow velocity increases by 10.9%, 12.3%, 14.0%, and 16.3%, respectively. As the pressure gradient increases, the boundary layer fluid participates in the flow, and the effective seepage section and fluid effective viscosity increase. The increase of fluid viscosity has retarding effect on the flow, but the effect is smaller than the gain effect of the increase in flow section on the flow velocity. Therefore, in general, the flow velocity increases with the increase of the pressure gradient.


Key words:  , groundwater, clayey soil, low-velocity non-Darcy flow, boundary layer, fluid viscosity, fractal theory

中图分类号: 

  • TU43
[1] 胡启智, 王秀娟, 丘锦荣, 康迪, 曾经文, 蔡倩怡, 刘人涛, 刘娜. 潮汕地区地下水饮用水源地酸性浅层地下水特征及其成因探讨[J]. 吉林大学学报(地球科学版), 2025, 55(3): 906-918.
[2] 张晟瑀, 申文超, 苏小四. 基于随机森林法的区域地下水硝酸盐污染风险评价[J]. 吉林大学学报(地球科学版), 2025, 55(3): 943-956.
[3] 倪嘉楠, 洪勇, 姜奕辰, 于超, 贺可强, 李亮. 基于离散元方法的黏性土剪切特性细观机理[J]. 吉林大学学报(地球科学版), 2025, 55(1): 223-232.
[4] 段磊, 曾经文, 赵显林, 丘锦荣, 刘娜, 陶钧实, 周建利. 珠江三角洲典型水产养殖区浅层地下水碘化物分布、来源及健康风险评估[J]. 吉林大学学报(地球科学版), 2024, 54(5): 1657-1674.
[5] 韩红闪, 朱琳, 郭高轩, 李炳华, 卢灿. 新水情背景下北京平原朝阳—通州沉降区地面沉降特征及其影响[J]. 吉林大学学报(地球科学版), 2024, 54(4): 1326-1338.
[6] 杨鹏, 张寿川, 郭本力, 杨崇敬, 冯启原, 吕有成, 解露茜. 山东日照地区不同类型含水层地下水位多年动态变化特征及其影响因素[J]. 吉林大学学报(地球科学版), 2024, 54(3): 967-979.
[7] 束龙仓, 张彤豪, 澈丽木格, 刘圣昱, 刘波. 降水和地下水开采对通辽市科尔沁区潜水位的影响[J]. 吉林大学学报(地球科学版), 2024, 54(3): 954-966.
[8] 刘国东, 杜成鸿, 侯杰, 杨梦溪, 陈宇, 谢杨.

基于用高斯分布生成管道的岩溶地下水流数值模拟 [J]. 吉林大学学报(地球科学版), 2024, 54(2): 592-603.

[9] 闫佰忠, 盖俊百, 王昕洲, 占新凯, 马苗苗. 石家庄滹沱河山前冲洪积扇地下水位动态演变特征及影响机制[J]. 吉林大学学报(地球科学版), 2023, 53(6): 1880-1891.
[10] 王福刚 杨国华, 程辉, 管小桐, 袁益龙. 热变形对饱和黏性土渗透系数的影响[J]. 吉林大学学报(地球科学版), 2023, 53(6): 1835-1844.
[11] 赵勇胜, 戴贞洧.

海泡石改性土-膨润土泥浆阻截墙阻截地下水重金属阳离子污染 [J]. 吉林大学学报(地球科学版), 2023, 53(5): 1549-1559.

[12] 李世杰, 骆祖江, 徐明钻, 张斌. 沿海地区地下水水质特征及其对地下水源热泵系统潜在性危害[J]. 吉林大学学报(地球科学版), 2023, 53(4): 1204-1215.
[13] 许天福, 陈敬宜, 冯波, 姜振蛟. 地热资源开发过程中潜在地下水环境问题[J]. 吉林大学学报(地球科学版), 2023, 53(4): 1149-1162.
[14] 王福刚, 盘惠林, 李胜伟, 王东辉, 吴铭杰, 平世飞, 曹玉清. 含膏岩地层地下水环境对混凝土建筑侵蚀实验[J]. 吉林大学学报(地球科学版), 2023, 53(4): 1185-1194.
[15] 王清, 吕作俊, 姚萌, 董佳祺, 夏玮彤, 杨天亮, 牛岑岑. 崇明东滩吹填区黏性土层抗剪强度随时间变化特征及机理[J]. 吉林大学学报(地球科学版), 2023, 53(4): 1163-1163.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 何辉, 孔垂显, 蒋庆平, 邓西里, 肖芳伟, 李顺明. 准噶尔盆地西北缘二叠系火山岩储层裂缝发育特征及分布预测—以金龙2井区佳木河组为例[J]. 吉林大学学报(地球科学版), 2015, 45(5): 1278 -1288 .
[2] 李三忠, 杨朝, 赵淑娟, 李玺瑶, 索艳慧, 郭玲莉, 余珊, 戴黎明, 李少俊, 牟墩玲. 全球早古生代造山带(Ⅱ):俯冲-增生型造山[J]. 吉林大学学报(地球科学版), 2016, 46(4): 968 -1004 .
[3] 袁红旗, 王蕾, 于英华, 张冬杰, 许凤鸣, 刘海涛. 沉积学粒度分析方法综述[J]. 吉林大学学报(地球科学版), 2019, 49(2): 380 -393 .
[4] 王世成, 杨仲杰, 杨菊, 张璟, 孙守亮, 刘长纯. 辽东石庙沟岩体岩石地球化学特征、锆石U-Pb年龄、Hf同位素及其地质意义[J]. 吉林大学学报(地球科学版), 2021, 51(2): 429 -441 .
[5] 王惠初, 陆松年, 初航, 相振群, 张长捷, 刘欢. 辽阳河栏地区辽河群中变质基性熔岩的锆石U-Pb年龄与形成构造背景[J]. J4, 2011, 41(5): 1322 -1334 .
[6] 田梦宇, 狄永军, 王帅, 贾一龙. 广西云开地区那蓬岩体黑云母二长花岗岩年代学、地球化学特征及成因[J]. 吉林大学学报(地球科学版), 2021, 51(3): 749 -766 .
[7] 谭晓淼, 高锐, 王海燕, 侯贺晟, 李洪强, 匡朝阳. 中亚造山带东段深地震反射剖面大炮揭露下地壳与Moho结构——数据处理与初步解释[J]. 吉林大学学报(地球科学版), 2021, 51(3): 898 -908 .
[8] 刘建民, 赵国春, 徐刚, 邱海成, 李建锋, 肖昌浩, 沙德铭, 刘福兴, 毕广源, 房兴, 张家奇, 郭祺, 于婳. 辽东半岛金矿成矿作用与深部资源勘查[J]. 吉林大学学报(地球科学版), 2021, 51(6): 1613 -1635 .
[9] 张兵强, 赵富远, 杨清毫, 黄毅, 李俊海, 刘松. 贵州省盘县架底金矿床成矿地质条件及找矿方向[J]. 吉林大学学报(地球科学版), 2022, 52(1): 94 .
[10] 王清, 汪洲, 韩梦霞, 刘经. 疏水材料处理吉林西部乾安地区盐渍土的力学强度试验[J]. 吉林大学学报(地球科学版), 2021, 51(3): 804 -814 .