›› 2018, Vol. 39 ›› Issue (3): 1063-1070.doi: 10.16285/j.rsm.2016.0599

• 数值分析 • 上一篇    下一篇

三维区域地面沉降数值模拟

罗 跃1, 2,叶淑君2,吴吉春2   

  1. 1. 东华理工大学 江西省数字国土重点实验室,江西 南昌 330013;2. 南京大学 地球科学与工程学院,江苏 南京 210046
  • 收稿日期:2016-03-28 出版日期:2018-03-12 发布日期:2018-06-06
  • 通讯作者: 叶淑君,女,1974年生,博士,教授,主要从事区域地面沉降的数值模拟研究工作。E-mail: sjye@nju.edu.cn E-mail:kuaikuaikaikai@126.com
  • 作者简介:罗跃,男,1985年生,博士,讲师,主要从事区域地面沉降的数值模拟研究工作。
  • 基金资助:

    东华理工大学江西省数字国土重点实验室开放研究基金资助项目(No.DLLJ201608);国家自然科学基金资助项目(No.41602258);江苏省自然科学基金资助项目(No.BK2012730);国土资源部地面沉降监测与防治重点实验室开放基金资助项目(No.20140012)

Numerical model for simulating 3D regional land subsidence

LUO Yue1, 2, YE Shu-jun2, WU Ji-chun2   

  1. 1. Key Laboratory for Digital Land Resources of Jiangxi Province, East China Institute of Technology, Nanchang, Jiangxi 330013, China; 2. School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu 210046, China
  • Received:2016-03-28 Online:2018-03-12 Published:2018-06-06
  • Supported by:

    This work was supported by the Key Laboratory for Digital Land and Resources of Jiangxi Province (DLLJ201608), the National Natural Science Foundation of China (41602258), the Natural Science Foundation of Jiangsu Province (BK2012730) and the Open Research Program of Key Laboratory of Land Subsidence Monitoring and Prevention (20140012).

摘要: 由于地裂缝研究及地表变形监测技术(例如GPS,InSAR等)的快速发展,抽取地下水引起的地层水平位移受到关注。传统区域地面沉降模型虽然求解快速但不能模拟水平位移;比奥模型虽然能够模拟土体的三维变形,但模型求解的计算量较大,较少应用于大尺度的区域地面沉降数值模拟。为解决以上问题,推导了解耦三维地面沉降数学模型,模型推导过程显示:比奥模型假设法向总应力和不变,则可简化为解耦三维地面沉降模型;解耦三维地面沉降模型假设土体仅有垂向一维变形,则可简化为传统区域地面沉降模型。同时通过数值试验验证了解耦三维地面沉降模型可以作为比奥模型的替代模型和传统区域地面沉降模型的改进模型,用来模拟抽取地下水引起的三维区域地面沉降。

关键词: 地面沉降, 地裂缝, 解耦三维地面沉降模型, 贮水率

Abstract: Horizontal deformation resulted from excessive groundwater exploitation has received attention in recent years because of earth fissures and technological development of deformation monitoring, such as GPS and InSAR. However, traditional uncoupling 1D model of land subsidence cannot simulate horizontal deformation. Though Biot model is able to simulate horizontal deformation, using this model is not appropriate for numerically solving the problem of regional land subsidence because the size of discrete model and time step is limited due to long computation time and ill-conditioned linear system matrix. To overcome inherent weaknesses above, uncoupling 3D mathematical model of land subsidence is developed by combining the benefits of uncoupling 1D model and Biot model. Uncoupling 3D mathematical model consists of groundwater flow and deformation equation. Both equations are coupled with parameters (i.e., Young modulus, Poisson ratio and specific water storage). Groundwater flow equation initially computes 3D flow field (equivalent fluid pore pressures) and then the deformation equation calculates 3D deformation based on the known flow field. The derivation of uncoupling 3D mathematical model shows that Biot model can be simplified to uncoupling 3D model based on the assumption that the total normal stress does not change in the process of groundwater flow; uncoupling 3D model can be simplified to uncoupling 1D model based on the assumption that the radial displacements vanish. Meanwhile, numerical experiment shows that uncoupling 3D model can provide deformation results similar to Biot model with much less running time compared with Biot model. Therefore, uncoupling 3D model can be considered as an alternative model to Biot model and an improved model for uncoupling 1D land subsidence model to simulate 3D regional land subsidence.

Key words: land subsidence, earth fissure, uncoupling 3D model of land subsidence, specific water storage

中图分类号: 

  • P 641.6

[1] 徐进, 王少伟, 杨伟涛. 水位变化下可压缩土层的黏弹性耦合变形分析[J]. 岩土力学, 2020, 41(3): 1065-1073.
[2] 熊仲明, 张 朝, 陈 轩. 地震作用下地裂缝场地地震动参数试验研究[J]. 岩土力学, 2019, 40(2): 421-428.
[3] 杨建民,霍王文,. 线状和面状群井抽水导致场外地面沉降呈s-lnr线性关系[J]. , 2018, 39(10): 3565-3572.
[4] 杨吉龙,袁海帆,胡云壮,胥勤勉,施佩歆,陈永胜,. 天津滨海地区深部黏土层弹塑性变形特征与地面沉降关系研究[J]. , 2018, 39(10): 3763-3772.
[5] 蔡奇鹏,吴宏伟,胡 平,陈星欣,郭力群,. 隐伏裂缝对互层胶结土中正断层扩展影响研究[J]. , 2017, 38(7): 2015-2021.
[6] 庄 超,周志芳,李兆峰,郭巧娜. 一种确定超固结弱透水层水力参数的方法[J]. , 2017, 38(1): 61-66.
[7] 王 瑞,胡志平,夏香波,苑韦虹,陈 悦,. 西安地裂缝界面处地震波传播规律[J]. , 2016, 37(S2): 642-648.
[8] 崔庆龙 ,沈水龙 ,吴怀娜 ,许烨霜 , . 广州岩溶地区深基坑开挖对周围环境影响的研究[J]. , 2015, 36(S1): 553-557.
[9] 邓亚虹 ,李 丽 ,慕焕东 ,王 鹏 ,李飞霞,. 西安地区地裂缝带Q3原状黄土流变特性试验研究[J]. , 2015, 36(7): 1847-1855.
[10] 杨吉龙 ,曹国亮 ,李 红 ,李 静 ,胡云壮 , 胥勤勉 ,秦雅飞 ,杜 东 ,方 成,. 天津滨海地区晚新生代地层自然固结与地面沉降研究[J]. , 2014, 35(9): 2579-2586.
[11] 王启耀 ,彭建兵 ,蒋臻蔚 ,滕宏泉,. 西安典型段地面沉降分层标观测及数值模拟[J]. , 2014, 35(11): 3298-3302.
[12] 徐玉龙 ,杨春和 ,李银平 ,施锡林 ,孔君凤 ,李浩然,. 基于改进的Mogi模型对樟树盐矿采空区地面沉降的分析[J]. , 2014, 35(10): 2894-2900.
[13] 林存刚 吴世明 张忠苗 刘俊伟 李宗梁 . 盾构掘进速度及非正常停机对地面沉降的影响[J]. , 2012, 33(8): 2472-2482.
[14] 许烨霜,马 磊,沈水龙. 上海市城市化进程引起的地面沉降因素分析[J]. , 2011, 32(S1): 578-0582.
[15] 李凯玲 ,门玉明 ,刘 洋 ,姜 容. 地裂缝剖面形态对地铁隧道变形影响模型试验研究[J]. , 2011, 32(6): 1690-1696.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!