不同方法求解疏排水引起的地面沉降对比研究
|
摘要
地下水开发利用和疏排水过程引起的地面沉降是一个渗流场与应力场耦合的过程。基于算例模型,分别采用目前常用的基于渗流-沉降分步计算的土力学经验公式、渗流-沉降部分耦合的GMS软件中SUB模型以及目前少用的基于渗流-沉降全耦合的COMSOL Multiphysics模型对疏排水引起的地面沉降量进行了计算。对比不同方法计算结果表明,COMSOL Multiphysics计算疏排水引起的地面沉降可行且计算结果较经验公式法、GMS中SUB模型的结果更为合理、更符合实际沉降特征。在此基础上,采用COMSOL Multiphysics和GMS中SUB模型对比研究了不同渗透系数对地面沉降计算的影响,结果进一步验证了COMSOL Multiphysics求解疏排水引起地面沉降的可靠性。
Land subsidence caused by groundwater exploitation and drainage is a process of coupling between the seepage field and the stress field. An example based on the model,using the current flow-settlement formula of soil mechanics calculation based on the commonly used step by step,seepage and settlement of partially coupled GMS software SUB model and the current use of drainage caused by COMSOL multiphysics model based on coupling seepage and settlement of ground settlement are calculated. Comparing the results of different methods,it is feasible to calculate land subsidence caused by drainage by COMSOL Multiphysics.The calculated results are more reasonable than the empirical formula method and the SUB model in GMS,which are more consistent with the actual settlement characteristics. The influence of different coefficients of permeability on the calculation of land subsidence is compared with the SUB model in COMSOL multiphysics and GMS. The results further verify the reliability of COMSOL multiphysics for solving the land subsidence caused by drainage.
Land subsidence caused by groundwater exploitation and drainage is a process of coupling between the seepage field and the stress field. An example based on the model,using the current flow-settlement formula of soil mechanics calculation based on the commonly used step by step,seepage and settlement of partially coupled GMS software SUB model and the current use of drainage caused by COMSOL multiphysics model based on coupling seepage and settlement of ground settlement are calculated. Comparing the results of different methods,it is feasible to calculate land subsidence caused by drainage by COMSOL Multiphysics.The calculated results are more reasonable than the empirical formula method and the SUB model in GMS,which are more consistent with the actual settlement characteristics. The influence of different coefficients of permeability on the calculation of land subsidence is compared with the SUB model in COMSOL multiphysics and GMS. The results further verify the reliability of COMSOL multiphysics for solving the land subsidence caused by drainage.
引文
[1]薛禹群,谢春红.地下水数值模拟[M].北京:科学出版社,2007:43-49.[XUE Y Q,XIE C H.Numerical simulation of groundwater[M].Beijing:Science Press,2007:43-49.(in chinese)]
[2]叶淑君,薛禹群,张云,等.上海区域地面沉降模型中土层变形特征研究[J].岩土工程学报,2005,27(2):140-147.[YE S J,XUE Y Q,ZHANG Y,et al.Study on the deformation characteristics of soil layer in the Shanghai regional ground subsidence model[J].Journal of Geotechnical Engineering,2005,27(2):140-147.(in Chinese)]
[3]SHI X Q,XUE Y Q,YE S J,et al.Characterization of land subsidence induced by groundwater withdrawals in Su-Xi-Chang area[J].China Environmental Geology,2007,52(1):27-40.
[4]ZHANG Y,XUE Y Q,WU J C,et al.Stress-strain measurements of deforming aquifer systems that underlie Shanghai,China[J].Environmental&Engineering Geoscience,2007,13(3):217-228.
[5]熊小锋,施小清,吴剑峰,等.弹塑性变形条件下抽水引起的地面沉降三维数值模拟[J].水文地质工程地质,2017,44(2):151-159.[XIONG X F,SHI X Q,WU J F,et al.3D numerical Simulation of elasto-plastic Land subsidence induced by groundwater pumping[J].Hydrogeology&Engineering Geology,2017,44(2):151-159.(in Chinese)]
[6]Lewis R W,B Schrefler.A fully coupled consolidation model of the subsidence of Venice[J].Water Resources Research,1978,14(2):223-230.
[7]金玮泽,骆祖江,陈兴贤,等.地下水渗流与地面沉降耦合模拟[J].地球科学,2014,39(5):611-619.[JIN Z W,LUO Z J,CHEN X X,et al.Coupled simulation of groundwater seepage and ground subsidence[J].Geoscience,2014,39(5):611-619.(in Chinese)]
[8]孙培德,杨东全,陈奕柏.多物理场耦合模型及数值模拟导论[M].北京:中国科学技术出版社,2007:166-172.[SUN P D,YANG D Q,CHEN Y B.Introduction to the multi physical field coupling model and numerical simulation[M].Beijing:China Science and Technology Press,2007:166-172.(in Chinese)]
[9]纪佑军,刘建军,程林松,等.考虑流-固耦合的隧道开挖数值模拟[J].岩土力学,2011,32(4):1229-1233.[JI Y G,LIU J J,CHENG L S,et al.Considering the numerical simulation of tunnel excavation with fluid solid coupling[J].Geodynamics,2011,32(4):1229-1233.(in Chinese)]
[10]刘波,张功,江永华,等.基于变渗透系数的深基坑单井抽水沉降研究[J].工程地质学报,2014,22(6):1123-1127.[LIU B,ZHANG G,JIANG Y H,et al.Study on single well pumping settlement of deep foundation pit based on variable permeability coefficient[J].Engineering Geology Journal,2014,22(6):1123-1127.(in Chinese)]
[11]耿冬青,贾振宇,韩建聪.关于降水引起地面沉降问题的探讨[J].施工技术,2005(增刊1):55-56.[GENG D Q,JIA Z Y,HAN J C.On the ground subsidence caused by precipitation[J].Construction Technology,2005(Sup1):55-56.(in chinese)]
[12]王艺伟.苏北沿海地区区域地面沉降模拟研究[D].南京:南京大学,2016.[WANG Y W.Simulation of regional ground subsidence in the coastal areas of North Jiangsu[D].Nanjing:Nanjing University,2016.(in Chinese)]
[13]Carslaw H S,J C Jaeger,H Feshbach.Conduction of heat in solids[M].Carolina:Clarendon Press,1959:1054-1057.
[14]李培超,孔祥言,卢德唐.饱和多孔介质流固耦合渗流的数学模型[J].水动力学研究与进展,2003,18(4):419-426.[LI P C,KONG X Y,LU T D.Mathematical model of fluid solid coupling seepage in saturated porous media[J].Hydrodynamics,2003,18(4):419-426.(in Chinese)]
[15]冉启全,李士伦.流固耦合油藏数值模拟中物性参数动态模型研究[J].石油勘探与开发,1997,24(3):61-65.[RAN Q Q,LI S L.Study on dynamic model of physical parameters in the numerical simulation of fluid solid coupling reservoir.[J].Petroleum Exploration and Development,1997,24(3):61-65.(in Chinese)]
[16]贺小桐,叶淑君,于军,等.基于固体颗粒速度场的三维地面沉降模拟[J].南京大学学报(自然科学),2015,51(6):1268-1278.[HE X T,YE S J,YU J,et al.Simulation of three-dimensional ground subsidence based on solid particle velocity field[J].Journal of Nanjing University(Natural Science),2015,51(6):1268-1278.(in chinese)]
[17]万伟锋.西安市地下水开采—地面沉降数值模拟及防治方案研究[D].西安:长安大学,2008.[WAN W F.Xi'an groundwater mining ground subsidence numerical simulation and prevention and control scheme[D].Xi'an:Chang'an University,2008.(in chinese)]
[18]苏晨,崔亚莉,邵景力.基于理论法的地面沉降数值模拟进展[J].水文地质工程地质,2014,41(6):147-152.[SU C,CUI Y L,SHAO J L.Anvances in numerical models of land subsidence based on comprehensive theories[J].Hydrogeology&Engineering Geology,2014,41(6):147-152.(in Chinese)]
[19]董成志,江思珉,夏学敏,等.基于MODFLOW-SUB模块的地层压缩与地面沉降模拟研究[J].工程勘察,2017(11):42-48.[DONG C Z,JIANG S M,XIA X M,et al.Simulation of stratigraphic compression and ground subsidence based on MODFLOW-SUB module[J].Geotechnical Investigation&Surveying,2017(11):42-48(in chinese)]
[20]熊小锋,罗跃,施小清,等.基于TOUGH2-FLAC3D耦合的三维地面沉降数值模拟及控制策略研究[J].高校地质学报,2017,23(1):172-180.[XIONG X F,LUO Y,SHI X Q,et al.Numerical simulation of three-dimensional ground subsidence based on TOUGH2-FLAC3D coupling and control strategy research[J].University Geo Journal,2017,23(1):172-180.(in chinese)]
[21]SHI X Q,XUE Y Q.Modlification of SUB package to simulation aquifer-system compaction and land subsidence and land subsidence[C]//Proceedings of the 7th international symposium on ground settlement.Shanghai:Shanghai Science and Technology Press,2005:24-29.
[2]叶淑君,薛禹群,张云,等.上海区域地面沉降模型中土层变形特征研究[J].岩土工程学报,2005,27(2):140-147.[YE S J,XUE Y Q,ZHANG Y,et al.Study on the deformation characteristics of soil layer in the Shanghai regional ground subsidence model[J].Journal of Geotechnical Engineering,2005,27(2):140-147.(in Chinese)]
[3]SHI X Q,XUE Y Q,YE S J,et al.Characterization of land subsidence induced by groundwater withdrawals in Su-Xi-Chang area[J].China Environmental Geology,2007,52(1):27-40.
[4]ZHANG Y,XUE Y Q,WU J C,et al.Stress-strain measurements of deforming aquifer systems that underlie Shanghai,China[J].Environmental&Engineering Geoscience,2007,13(3):217-228.
[5]熊小锋,施小清,吴剑峰,等.弹塑性变形条件下抽水引起的地面沉降三维数值模拟[J].水文地质工程地质,2017,44(2):151-159.[XIONG X F,SHI X Q,WU J F,et al.3D numerical Simulation of elasto-plastic Land subsidence induced by groundwater pumping[J].Hydrogeology&Engineering Geology,2017,44(2):151-159.(in Chinese)]
[6]Lewis R W,B Schrefler.A fully coupled consolidation model of the subsidence of Venice[J].Water Resources Research,1978,14(2):223-230.
[7]金玮泽,骆祖江,陈兴贤,等.地下水渗流与地面沉降耦合模拟[J].地球科学,2014,39(5):611-619.[JIN Z W,LUO Z J,CHEN X X,et al.Coupled simulation of groundwater seepage and ground subsidence[J].Geoscience,2014,39(5):611-619.(in Chinese)]
[8]孙培德,杨东全,陈奕柏.多物理场耦合模型及数值模拟导论[M].北京:中国科学技术出版社,2007:166-172.[SUN P D,YANG D Q,CHEN Y B.Introduction to the multi physical field coupling model and numerical simulation[M].Beijing:China Science and Technology Press,2007:166-172.(in Chinese)]
[9]纪佑军,刘建军,程林松,等.考虑流-固耦合的隧道开挖数值模拟[J].岩土力学,2011,32(4):1229-1233.[JI Y G,LIU J J,CHENG L S,et al.Considering the numerical simulation of tunnel excavation with fluid solid coupling[J].Geodynamics,2011,32(4):1229-1233.(in Chinese)]
[10]刘波,张功,江永华,等.基于变渗透系数的深基坑单井抽水沉降研究[J].工程地质学报,2014,22(6):1123-1127.[LIU B,ZHANG G,JIANG Y H,et al.Study on single well pumping settlement of deep foundation pit based on variable permeability coefficient[J].Engineering Geology Journal,2014,22(6):1123-1127.(in Chinese)]
[11]耿冬青,贾振宇,韩建聪.关于降水引起地面沉降问题的探讨[J].施工技术,2005(增刊1):55-56.[GENG D Q,JIA Z Y,HAN J C.On the ground subsidence caused by precipitation[J].Construction Technology,2005(Sup1):55-56.(in chinese)]
[12]王艺伟.苏北沿海地区区域地面沉降模拟研究[D].南京:南京大学,2016.[WANG Y W.Simulation of regional ground subsidence in the coastal areas of North Jiangsu[D].Nanjing:Nanjing University,2016.(in Chinese)]
[13]Carslaw H S,J C Jaeger,H Feshbach.Conduction of heat in solids[M].Carolina:Clarendon Press,1959:1054-1057.
[14]李培超,孔祥言,卢德唐.饱和多孔介质流固耦合渗流的数学模型[J].水动力学研究与进展,2003,18(4):419-426.[LI P C,KONG X Y,LU T D.Mathematical model of fluid solid coupling seepage in saturated porous media[J].Hydrodynamics,2003,18(4):419-426.(in Chinese)]
[15]冉启全,李士伦.流固耦合油藏数值模拟中物性参数动态模型研究[J].石油勘探与开发,1997,24(3):61-65.[RAN Q Q,LI S L.Study on dynamic model of physical parameters in the numerical simulation of fluid solid coupling reservoir.[J].Petroleum Exploration and Development,1997,24(3):61-65.(in Chinese)]
[16]贺小桐,叶淑君,于军,等.基于固体颗粒速度场的三维地面沉降模拟[J].南京大学学报(自然科学),2015,51(6):1268-1278.[HE X T,YE S J,YU J,et al.Simulation of three-dimensional ground subsidence based on solid particle velocity field[J].Journal of Nanjing University(Natural Science),2015,51(6):1268-1278.(in chinese)]
[17]万伟锋.西安市地下水开采—地面沉降数值模拟及防治方案研究[D].西安:长安大学,2008.[WAN W F.Xi'an groundwater mining ground subsidence numerical simulation and prevention and control scheme[D].Xi'an:Chang'an University,2008.(in chinese)]
[18]苏晨,崔亚莉,邵景力.基于理论法的地面沉降数值模拟进展[J].水文地质工程地质,2014,41(6):147-152.[SU C,CUI Y L,SHAO J L.Anvances in numerical models of land subsidence based on comprehensive theories[J].Hydrogeology&Engineering Geology,2014,41(6):147-152.(in Chinese)]
[19]董成志,江思珉,夏学敏,等.基于MODFLOW-SUB模块的地层压缩与地面沉降模拟研究[J].工程勘察,2017(11):42-48.[DONG C Z,JIANG S M,XIA X M,et al.Simulation of stratigraphic compression and ground subsidence based on MODFLOW-SUB module[J].Geotechnical Investigation&Surveying,2017(11):42-48(in chinese)]
[20]熊小锋,罗跃,施小清,等.基于TOUGH2-FLAC3D耦合的三维地面沉降数值模拟及控制策略研究[J].高校地质学报,2017,23(1):172-180.[XIONG X F,LUO Y,SHI X Q,et al.Numerical simulation of three-dimensional ground subsidence based on TOUGH2-FLAC3D coupling and control strategy research[J].University Geo Journal,2017,23(1):172-180.(in chinese)]
[21]SHI X Q,XUE Y Q.Modlification of SUB package to simulation aquifer-system compaction and land subsidence and land subsidence[C]//Proceedings of the 7th international symposium on ground settlement.Shanghai:Shanghai Science and Technology Press,2005:24-29.