摘要(Abstract)：
为了揭示不同大孔隙参数取值条件下边坡水分场的变化规律，利用2个耦合的Richards方程描述大孔隙流，并联合运动波方程，建立坡面径流与大孔隙流耦合模型；基于有限元软件COMSOL的偏微分方程接口，实现所建立耦合模型的数值求解，并自行设计室内模型试验，验证数值结果，分析不同大孔隙参数对体积含水率和坡面积水深度的影响。结果表明：COMSOL软件的偏微分方程接口可以实现所建立耦合模型的数值求解，相比于无大孔隙坡面，考虑大孔隙时水分入渗深度明显更大；边坡饱和区深度与湿润锋深度均随大孔隙域占比的增大而增大，均随大孔隙域与基质域导水系数之比与经验参数的增大而减小；当降雨历时为30 min时，坡面积水深度区分度最大，表现在随着大孔隙域占比的增大而减小，均随大孔隙域与基质域导水系数之比和经验参数的增大而增大；3种大孔隙参数按照对边坡水分场影响程度由大到小的顺序为大孔隙域占比、经验参数、大孔隙域与基质域导水系数之比。
关键词(KeyWords)： 坡面径流;水分场;有限元法;大孔隙流;湿润锋
基金项目(Foundation): 国家自然科学基金项目(41772297);; 福建省中青年教师教育科研项目(JAT210542)
作者(Author): 李尚辉,刘代文,阙云,陈宪,丁辉
DOI: 10.13349/j.cnki.jdxbn.20230309.003
参考文献(References)：
[1] 杨寅，包红军，彭涛.台风“鲇鱼”强降水引发的地质灾害气象风险预警检验与分析[J].暴雨灾害，2019,38(3):221.
[2] 靳红华，王林峰，任青阳，等.降雨循环条件下高切坡稳定性演化过程及预测方法[J].土木与环境工程学报(中英文),2021,43(4):12.
[3] 许旭堂.非饱和残积土坡对降雨入渗的响应及其失稳演变过程研究[D].福州：福州大学，2015.
[4] 徐则民，黄润秋.山区流域高盖度斜坡对极端降雨事件的地下水响应[J].地球科学进展，2011,26(6):598.
[5] 徐宗恒，徐则民，曹军尉，等.土壤优先流研究现状与发展趋势[J].土壤，2012,44(6):905.
[6] BOGAARD T A,GRECO R.Landslide hydrology:from hydrology to pore pressure[J].Wiley Interdisciplinary Reviews:Water,2016,3(3):439.
[7] JARVIS N,KOESTEL J,LARSBO M.Understanding preferential flow in the vadose zone:recent advances and future prospects[J].Vadose Zone Journal,2016,15(12):1.
[8] 柯云斌，符元帅，阙云，等.降雨诱发残积土陡坡坡面冲刷破坏颗粒流模拟[J].济南大学学报(自然科学版),2019,33(5):453.
[9] 王乐，田东方.边坡渗流与坡面径流联合求解三维有限元模型[J].人民珠江，2019,40(5):24.
[10] TIAN D F,ZHENG H,LIU D F.A 2D integrated FEM model for surface water-groundwater flow of slopes under rainfall condition[J].Landslides,2017,14(2):577.
[11] 韩同春，苏钰钦，张宇.双层结构边坡降雨入渗与坡面径流耦合分析[J].工程科学与技术，2020,52(6):145.
[12] 童富果，田斌，刘德富.改进的斜坡降雨入渗与坡面径流耦合算法研究[J].岩土力学，2008,29(4):1035.
[13] SHIGORINA E,RüDIGER F,TARTAKOVSKY A M,et al.Multiscale smoothed particle hydrodynamics model development for simulating preferential flow dynamics in fractured porous media[J].Water Resources Research,2020,57(3):1.
[14] SUN X,LUO H,KENICHI S.A coupled thermal-hydraulic-mechanical-chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics[J].Journal of Zhejiang University:Science A,2018,19(8):600.
[15] 孟庆娟，石志飞.基于周期理论和COMSOL PDE的排桩减振特性研究[J].岩土力学，2018,39(11):4252.
[16] 刘俊新，刘育田，胡启军.非饱和地表径流-渗流和流固体耦合条件下降雨入渗对路堤边坡稳定性研究[J].岩土力学，2010,31(3):904.
[17] GERKE H H,VAN GENUCHTEN M T.A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media[J].Water Resources Research,1993,29(2):305.
[18] 雷志栋，杨诗秀，谢森传.土壤水动力学[M].北京：清华大学出版社，1988:4.
[19] 张培文，刘德富，郑宏，等.降雨条件下坡面径流和入渗耦合的数值模拟[J].岩土力学，2004,25(1):111.
[20] KHNE J M,MOHANTY B P.Water flow processes in a soil column with a cylindrical macropore:experiment and hierarchical modeling[J].Water Resources Research,2005,41(3):1.
[21] NI S M,ZHANG D Q,FENG S Y.Quantitative relationship between hydraulics parameters and soil erosion rate on remolded soil slopes with different textures[J].Acta Pedologica Sinica,2019,56(6):1336.
[22] 交通运输部公路科学研究院.公路土工试验规程：JTG 3430—2020[S].北京：人民交通出版社，2020:8.
