江苏沿江开发带地下水开采三维渗流场一地面沉降耦合数值模拟
|
摘要
随着国民经济的发展,地下水过量开采引发的地面沉降地质灾害问题日益严重。严峻的形势对地下水资源评价与管理提出了新的要求—地下水资源评价必须与地面沉降预测相结合,以达到在获得最大地下水开采量的同时,将地面沉降控制到最小的目的。
本文以江苏沿江开发带第四纪松散孔隙含水层为研究对象,建立了研究区的地下水开采三维渗流场与地面沉降耦合模型,运用有限差分法对该耦合模型进行数值求解,对该区域地下水位和地面沉降的动态变化趋势进行了预测并提出了优化开采方案。其研究成果为该区域地下水资源的合理开发、地面沉降的有效控制提供了科学依据。主要研究内容如下:
1、系统分析沿江开发带的地质、水文地质和工程地质条件,概化该区域水文地质特点和地面沉降产生机理,建立地下水开采与地面沉降地质概念模型。
2、建立沿江开发带地下水开采三维渗流场与地面沉降耦合模型,运用有限差分法对模型进行求解并编制计算机计算程序。
3、选取2003年1月1日至2003年12月31日的实测水位观测资料,对模型的水文地质参数进行识别、校正,采用类比法确定地面沉降参数,运用编制的程序建立了可视化模型。
4、预测该区域在地下水现状开采模式下今后10年的地下水位和地面沉降量变化趋势。针对现状开采模式下地面沉降预测量过大的问题,在维持地下水现有开采量和尽可能控制地面沉降的前提下,通过调整开采布局,提出地下水优化开采方案。
With the development of the national economy, the land subsidence of geological disaster becomes worse because of excessive pumping groundwater. Under the pressure of the urgent circumstances, the new demand from evaluation and management of the groundwater has been put forward that the evaluation of the groundwater must be combined with subsidence for controlling the subsidence to the greatest extent while getting the maximum quantity of groundwater resources.The quaternary-loose-pore aquifer in developing belt along Yangtze River at Jiangsu Province is studied. A coupling model between 3D groundwater seepage and subsidence is built. By using finite-difference approach, the coupling model is simulated numerically. After modifying and calibrating the numerical model, dynamic changes of groundwater head and subsidence are assessed. All conclusions we sum up in the paper provide scientific reference for the reasonable evaluation of groundwater resources and efficient control of subsidence development for the belt. The main contents are as follows:Firstly, the geologic, hydro-geologic and engineering geologic conditions of the belt are analyzed systemically. After summarizing the hydro-geologic feature and the mechanism of subsidence development, a conceptual model between groundwater and subsidence is built.Secondly, a coupling model between 3D seepage field of groundwater exploitation and subsidence is set up. With the finite-difference approach, the coupling model is simulated numerically and a computer program is written.Thirdly, by using observation borehole datum from January 1, 2003 to December 31, 2003, the hydro-geologic parameters are identified and calibrated. With the method of analogy, the parameters of subsidence are confirmed. Then a visualized computer model is built by calculation program according to the geological condition.Finally, the changes of the hydraulic head and land subsidence for the next ten years are predicted with the existing exploitation's scheme. Aiming at the problem of the excessive subsidence with the existing exploitation scheme, an exploitation project is put forward on the condition of controlling land subsidence possibly which is exploited with the existing exploitation's scheme by adjusting exploitation layout.
本文以江苏沿江开发带第四纪松散孔隙含水层为研究对象,建立了研究区的地下水开采三维渗流场与地面沉降耦合模型,运用有限差分法对该耦合模型进行数值求解,对该区域地下水位和地面沉降的动态变化趋势进行了预测并提出了优化开采方案。其研究成果为该区域地下水资源的合理开发、地面沉降的有效控制提供了科学依据。主要研究内容如下:
1、系统分析沿江开发带的地质、水文地质和工程地质条件,概化该区域水文地质特点和地面沉降产生机理,建立地下水开采与地面沉降地质概念模型。
2、建立沿江开发带地下水开采三维渗流场与地面沉降耦合模型,运用有限差分法对模型进行求解并编制计算机计算程序。
3、选取2003年1月1日至2003年12月31日的实测水位观测资料,对模型的水文地质参数进行识别、校正,采用类比法确定地面沉降参数,运用编制的程序建立了可视化模型。
4、预测该区域在地下水现状开采模式下今后10年的地下水位和地面沉降量变化趋势。针对现状开采模式下地面沉降预测量过大的问题,在维持地下水现有开采量和尽可能控制地面沉降的前提下,通过调整开采布局,提出地下水优化开采方案。
With the development of the national economy, the land subsidence of geological disaster becomes worse because of excessive pumping groundwater. Under the pressure of the urgent circumstances, the new demand from evaluation and management of the groundwater has been put forward that the evaluation of the groundwater must be combined with subsidence for controlling the subsidence to the greatest extent while getting the maximum quantity of groundwater resources.The quaternary-loose-pore aquifer in developing belt along Yangtze River at Jiangsu Province is studied. A coupling model between 3D groundwater seepage and subsidence is built. By using finite-difference approach, the coupling model is simulated numerically. After modifying and calibrating the numerical model, dynamic changes of groundwater head and subsidence are assessed. All conclusions we sum up in the paper provide scientific reference for the reasonable evaluation of groundwater resources and efficient control of subsidence development for the belt. The main contents are as follows:Firstly, the geologic, hydro-geologic and engineering geologic conditions of the belt are analyzed systemically. After summarizing the hydro-geologic feature and the mechanism of subsidence development, a conceptual model between groundwater and subsidence is built.Secondly, a coupling model between 3D seepage field of groundwater exploitation and subsidence is set up. With the finite-difference approach, the coupling model is simulated numerically and a computer program is written.Thirdly, by using observation borehole datum from January 1, 2003 to December 31, 2003, the hydro-geologic parameters are identified and calibrated. With the method of analogy, the parameters of subsidence are confirmed. Then a visualized computer model is built by calculation program according to the geological condition.Finally, the changes of the hydraulic head and land subsidence for the next ten years are predicted with the existing exploitation's scheme. Aiming at the problem of the excessive subsidence with the existing exploitation scheme, an exploitation project is put forward on the condition of controlling land subsidence possibly which is exploited with the existing exploitation's scheme by adjusting exploitation layout.
引文
[1] 谭荣初.吴江市地面沉降与开采地下水关系的研究[J].地质灾害与环境保护,2001,12(3):24-28
[2] 缪晓图.苏州城市规划区Ⅱ承压水开采与地面沉降预防控制研究[J].中国地质灾害与防治学报,2004,15(3):55-59
[3] 朱兴贤,朱锦旗.苏锡常地区地面沉降灾害与经济损失分析[J].水文地质工程地质,1997,(3)24-25.
[4] 薛禹群.地下水资源与江苏地面沉降研究[J].江苏地质,2001,25(4):193-195
[5] 盛海洋,孟秋立,朱殿华等.我国地下水开发利用中的水环境问题及其对策[J].水土保持研究,2006,13(1):51-53
[6] 张阿根,魏子新.中国地面沉降研究[M].上海,上海科学技术出版社2005:4-9
[7] 张云.苏南地区地面沉降概述[J].地质灾害与环境保护,1999,10(3):66-71
[8] 许乃政,姜月华.我国东南沿海地面沉降类型及其特点[J].灾害学,2005,20(4):67-72
[9] 吴世明,周健,徐建平.(岩土工程新进展丛书)环境岩土工程[M].上海:同济大学出版社,2000.
[10] 赵常洲,龚固培,王晖.地面沉降成因与危害[J].西部探矿工程,2006(1):261-263
[11] 郑铣鑫,武强.关于城市地面沉降研究的几个前沿问题[J].地球学报,2002,23(3):2:79-282
[12] 凌荣华等.地面沉降模拟分析中几个基本问题的讨论[J].中国地质灾害与防治学报,1995,6(4):53-561
[13] 长江三角洲上海地区第四纪地质水文地质、工程地质及地面沉降数学模型研究.中国-比利时联合研究报告[R].1989
[14] 张云,薛禹群.抽水地面沉降数学模型的研究现状与展望[J].中国地质灾害与防治学报,2002.13(2):1-6
[15] 武胜忠,方鹏飞.地面沉降的计算理论和方法[J].太原理工大学学报,2000,21(3):162-166
[16] 阎世骏,刘长礼.城市地面沉降研究现状与展望[J].地学前缘,1996.3(1-2):96-98
[17] RolandoBravo,ARivera,等.地质矿产部地质环境管理司组织翻译.地面沉降—第四届地面沉降国际讨论会译文选集[C]..北京:地震出版社,1994.
[18] 宫辉力.郑州市地下水渗流场的数值模拟和优化管理[J].北京大学学报(自然科学版),1998,34(6):827-835
[19] 骆祖江,李朗等.复合含水层地区深基坑降水三维渗流场数值模拟-以上海环球金融中心为例[J].工程地质学报,2006,14(1)72-77
[20] 曹文炳.应用结合水渗流机制说明粘性土释水机制的初步探讨[J].全国第二届地面沉降 学术讨论会,1980
[21] 曹文炳,李克文.水位升降引起的粘性土层释水、吸水与越流发生过程的室内研究方法[J].勘察科学技术.No.4,1986
[22] 朱百里,沈珠江.计算土力学[M].上海:上海科学技术出版社,1990.
[23] A C Waltham. Ground Subsidence [M].Blackie&Son Limited, London, 1989
[24] Burbey.T.J, D.C.Helm.Modeling three-dimensional deformation in response to pumping of unconsolidated aquifers [J].Environmental and environmental 5, 1999, no.2:199-212
[25] Jiang Li.A nonlinear elastic solution for subsidence due to ASR applications to multi-aquifer systems [A].In: edited by Laura Carbognin, Giuseppe Gambolati A. Ivan Johnson. Proceedings of the Sixth International Symposium on Land Subsidence[C].Ravenna, Italy, 2000
[26] 吴林高,缪俊发.抽灌水作用下土层变形及应力应变本构律的研究[J].地球科学,1995(5):581-588
[27] JiangLi, Donald C Helm.A nonlinear viscous model for aquifer compression associated with ASR applications [A].In: Edited by Laura Carbognin, Giuseppe Gambolati, and A. Ivan Johnson. Proceedings of the Sixth Intemational Symposium on Land Subsidence [C], Ravenna, Italy, 2000
[28] 冉启全,顾小芸.考虑流变特性的流固耦合地面沉降计算模型[J].中国地质灾害与防治学报,1998.9(2).
[29] X Y Gu, Q Q Ran. 3-D coupled with consideration of theological properties [A]. In: edited by Laura Carbognin, Giuseppe Gambolati, A. Ivan Johnson. Proceedings of the Sixth International Symposium on Land Subsidence [C], Ravenna, Italy, 2000
[30] Terzaghi K. Theoretical soil mechanics [M].Wiley, New York, 1943
[31] G Gambolati, R Allan Frezze. Mathematical simulation of the subsidence of Venice [J]. Theory Water Resource Research,1973,9(3):721-733
[32] 崔小东.Modflow和IDP在天津地面沉降数值计算中的应用与开发[J].中国地质灾害与防治学报,1998,9(2):122-128
[33] T R Shearer. A numerical model to calculate land subsidence, applied at Han gu in China [J]. Engineering Geology, 1998, 49(2):85-93
[34] R W Lewis, B.Schrefler.A fully coupled consolidation model of the subsidence of Venice[J].Water Resource Research, 1978, 14(2).
[35] J P His, J P Carter, J C Small.Surface subsidence and drawdown of the watertable due to pumping [J]. Geotech—nique44.1994(3):381-396
[36] Miau B Su, C L Su, C J Chang, et al.A numerical model of ground deformation induced by single well pumping[J]. Computer and Geotechnics, 1998(2):39-60.
[37] 张克绪等.开采地下水引起地面变形的分析[J].自然灾害学报,1996,5(4):50-591
[38] Biot M A. General theory of three dimensional consolidation [J].J.Appl.Phys, 1941, 12:155～164.
[39] Biot M A. General solution of the equation of elasticity and consolidation for a porous material [J]. J.Appl.Phys. 1956, 78:91～96
[40] Chai J C, Shen S L, Zhu H H, Zhang X L. Land subsidence due to groundwater drawdown in Shanghai [J]. Geotechnique, 2004, 54(2): 143～148.
[41] Shen S L, Tohno I, Nishigaki M, Miura N. 3D-analysis of land subsidence due to withdrawal of groundwater [A]. Geotechnical Engineering in Urban Construction [C].Beijing: Tsinghua University Press, 2003.431-436.
[42] Gambolati G, Monaco A, Galeati G, et al. New approaches and applications in subsurface flow modeling: 3D finite element analysis of dewateringfor an electro-nuclear plan [A]. In: Groundwater Flow and Quality Modeling[C].Boston: D. Reidel Publishing Co., 1998.11～15.
[43] Cassiani G, Palozzo W, Zoccatelli C. Simulation of Fig.8 Steady flow system with an impervious wall subsidence caused by gas production from an offshore field and comparison with field data [J]. Journal of Petroleum Science and Engineering, 2002, 45(1): 123～134.
[44] 上海市环境地质站,比利时国家科技部.上海市地面沉降教学模型研究[R].1989
[45] 朱锡冰,奚建国等.上海浦东新区“地下水水量—水位—沉降联合数学模型应用”[R].1995
[46] 方正,李勤奋等.上海市区域地下水调查[R].1999
[47] 陈崇希,裴顺平.地下水开采—地面沉降模型研究[J].水文地质工程地质,2001,2:5-8
[48] 方鹏飞,朱向荣.太原市地面沉降的计算与预测[J].煤田地质与勘探,2002,30(4):44-46
[49] 崔亚莉,邵景力,谢振华等.基于MODFLOW的地面沉降模型研究—以北京市区为例[J].工程勘察,2003,5:19-22
[50] 张利生,朱国荣.鹏山水源地开采条件下的含水层固结模型[J].岩土力学,2005,26(7):1141-1147
[51] 吴振样.地下热水三维有限元数值模拟与地面沉降研究[D],2004
[52] 卢立志,隋兆显等.江苏省沿江工业走廊水文地质、工程地质、环境地质综合评价报告[R].1989.
[53] Jorn Hoffmann, S.A.Leake, etc. MODFLOW-2000 Ground-Water Model—User Guide to the Subsidence and Aquifer-System CompactionPackage[R].2000
[54] 钱家欢,殷宗泽等.土工计算与原理[M]..北京:中国水利水电出版社,2003
[55] 朱学愚.地下水资源评价[M].南京:南京大学出版社,1987
[56] 骆祖江,武永霞.盐城市地下水资源规划评价三维数值模型[J].水资源保护,2005,21(5):37-41
[57] Terzaghi K. Theoretical soil mechanics [M].Wiley, New York, 1943
[58] 罗国煜,李生林.工程地质学基础[M].南京:南京大学出版社,1990.
[59] 殷宗泽.土体的沉降与固结[M].北京:中国电力出版社,1998.
[60] Poland, J.F., and Davis, G.H. Land subsidence due to withdrawals of fluids, in Varnes, D.J., and Kiersch, G., eds.: Geological Society of America Reviews in Engineering Geology, v. 2, p. 187-269.
[61] 薛禹群,朱学愚等.地下水动力学[M].北京:地质出版社,1997:5.
[62] Wen Hsing Chiang, Wolfgang Kinzelbach.3D Groundwater modeling with PMWIN—A Simulation System for Modeling Groundwater Flow and Pollution[M]. Springer, 2001.
[63] Jorgensen, D.G., 1980, Relationships between basic soils-engineering equations and basic ground-waterflow equations:[R].U.S. Geological Survey Water-Supply Paper 2064, 40 p.
[64] Johnson, A.I., Moston, R.P., and Morris, D.A., 1968, Physical and hydrologic properties of water-bearing deposits in subsiding areas in California[R].U.S. Geological Survey Professional Paper 497-A, 71 p.
[65] Riley, F.S., 1969, Analysis of borehole extensometer data from central California: [A].International Association of Scientific Hydrology Publication 89, p. 423-431.
[66] Poland, J.F., and Ireland, R.L., 1988, Land subsidence in the Santa Clara Valley, California, as of 1982: [A] U.S. Geological Survey Professional Paper 497-F, 61 p.
[67] Heywood, C.E., 1997, Piezometric-extensometric estimations of specific storage in the Albuquerque Basin, New Mexico, in Prince, K.R. and Leake, S.A., eds., U.S. Geological Survey Subsidence lnterest Group Conference: Proceedings of the technical meeting, Las Vegas, Nevada, February 14-16, [A]. 1995: U.S. Geological Survey Open-File Report 97-47, p. 21-26
[68] Joseph F. Poland. Mechanics of aquifer systems—The scientific legacy of, Land subsidence case studies and current research: [A].Proceedings of the Dr. Joseph F. Poland Symposium on Land Subsidence, Association of Engineering Geologists Special Publication 1998,no. 8, p. 13-27.
[69] A C Waltham. Ground Subsidence[M].Blackie&Son Limited, London, 1989.
[70] T.J.Burbey. Stress-strain analyses for aquifer-system characterization [J].Ground Water, 2001, 39(1):128-136.
[71] T J.Berbey.Effects of horizontal strain in estimating specific storage and compaction in confined and leaky aquifer systems [J]. Hydrogeology Journal, 1999, (7)
[72] 陈崇希,唐仲华.地下水流动问题数值方法[M].武汉:中国地质大学出版社,1994.
[73] 孙讷正.地下水流的数值模型和数值方法[M].北京:地质出版社,1981.
[74] 孙峰根,王心田.水文地质计算的数值方法[M].徐州:中国矿业大学出版社,1995.
[75] 张蔚榛.地下水非稳定流计算和地下水资源评价[M].北京:科学出版社,1983.
[76] Leake, S.A., and Prudic, D.E., Documentation of a computer program to simulate aquifer-systemcompactionusingthemodular finite-differenceground-water flow model[J].U.S. Geological Survey Techniques of Water-Resources Investigations, book 6, 1991,chap.A2,68 p.
[77] 陈玉田.偏微分方程数值解法[M].南京:河海大学出版社,1999
[78] 顾尔祚.流体力学中的有限差分法基础[M].上海:上海交通大学出版社,1988
[79] 蒋亚萍,陈余道.MODFLOW— 一套水文地质学实用计算软件[J].广西地质,1999,12(3):75-78
[80] User's manual for Visual MODFLOW. [R].Waterloo Hydrogeologic Inc. 1999
[81] 谭浩强,田淑清.FORTRAN语言程序设计[M].北京:高等教育出版社,1997
[82] 马瑞民,衣治安.FORTRAN90程序设计[M].哈尔滨工程大学出版社.1998
[83] 范啸涛,季光明.预优矩阵及其构造技术[J].成都理工大学学报(自然科学版),2003.30(4):P432-435
[84] 李裕伟,任效颖等..二维地质统计学[M].北京:地质出版社,1994.
[85] 李国敏,陈崇希.克立格法在地下水数值模型中的应用[J].地质科技情报,1995,14(3):95-99
[86] MITASOVA H, MITAS L, BROWN W M, et al. Modeling spatially and temporally distributed phenomena: new method and tools for GRASS GIS[J].IntJGIS, 1995(4):433-446.
[87] 楚明敏,王彩会等.苏州城市规划区浅层地下水开发利用研究报告[R].2003
[88] 郭坤一,骆祖江等.长江三角洲地区地下水资源与地质灾害调查评价[R].2003
[89] 于军,宗开红等.基于GMS平台的苏锡常地区地下水流模型研究报告[R].2003
[90] 祝晓彬,吴吉春等.长江三角洲(长江以南)地区深层地下水三维数值模拟[J].地理科学,2005,25(1):68-73
[91] 姜蓓蕾,吴吉春等.常州市滆湖典型地块浅层地下水数值模拟[J].水文地质工程地质,2004(6):29-32
[92] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第一季度)[R].2003
[93] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第二季度)[R].2003
[94] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第三季度)[R].2003
[95] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第四季度)[R].2004
[96] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第一季度)[R].2003
[97] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第二季度)[R].2003
[98] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第三季度)[R].2003
[99] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第四季度)[R].2004
[100] 江苏省扬州市水利局.扬州市地下水动态监测汇总表(2003年度)[R].2004
[101] 国家技术监督局.中华人民共和国国家标准(GBT14497-93)地下水资源管理模型工作技术要求[S].1994,3:1～15
[2] 缪晓图.苏州城市规划区Ⅱ承压水开采与地面沉降预防控制研究[J].中国地质灾害与防治学报,2004,15(3):55-59
[3] 朱兴贤,朱锦旗.苏锡常地区地面沉降灾害与经济损失分析[J].水文地质工程地质,1997,(3)24-25.
[4] 薛禹群.地下水资源与江苏地面沉降研究[J].江苏地质,2001,25(4):193-195
[5] 盛海洋,孟秋立,朱殿华等.我国地下水开发利用中的水环境问题及其对策[J].水土保持研究,2006,13(1):51-53
[6] 张阿根,魏子新.中国地面沉降研究[M].上海,上海科学技术出版社2005:4-9
[7] 张云.苏南地区地面沉降概述[J].地质灾害与环境保护,1999,10(3):66-71
[8] 许乃政,姜月华.我国东南沿海地面沉降类型及其特点[J].灾害学,2005,20(4):67-72
[9] 吴世明,周健,徐建平.(岩土工程新进展丛书)环境岩土工程[M].上海:同济大学出版社,2000.
[10] 赵常洲,龚固培,王晖.地面沉降成因与危害[J].西部探矿工程,2006(1):261-263
[11] 郑铣鑫,武强.关于城市地面沉降研究的几个前沿问题[J].地球学报,2002,23(3):2:79-282
[12] 凌荣华等.地面沉降模拟分析中几个基本问题的讨论[J].中国地质灾害与防治学报,1995,6(4):53-561
[13] 长江三角洲上海地区第四纪地质水文地质、工程地质及地面沉降数学模型研究.中国-比利时联合研究报告[R].1989
[14] 张云,薛禹群.抽水地面沉降数学模型的研究现状与展望[J].中国地质灾害与防治学报,2002.13(2):1-6
[15] 武胜忠,方鹏飞.地面沉降的计算理论和方法[J].太原理工大学学报,2000,21(3):162-166
[16] 阎世骏,刘长礼.城市地面沉降研究现状与展望[J].地学前缘,1996.3(1-2):96-98
[17] RolandoBravo,ARivera,等.地质矿产部地质环境管理司组织翻译.地面沉降—第四届地面沉降国际讨论会译文选集[C]..北京:地震出版社,1994.
[18] 宫辉力.郑州市地下水渗流场的数值模拟和优化管理[J].北京大学学报(自然科学版),1998,34(6):827-835
[19] 骆祖江,李朗等.复合含水层地区深基坑降水三维渗流场数值模拟-以上海环球金融中心为例[J].工程地质学报,2006,14(1)72-77
[20] 曹文炳.应用结合水渗流机制说明粘性土释水机制的初步探讨[J].全国第二届地面沉降 学术讨论会,1980
[21] 曹文炳,李克文.水位升降引起的粘性土层释水、吸水与越流发生过程的室内研究方法[J].勘察科学技术.No.4,1986
[22] 朱百里,沈珠江.计算土力学[M].上海:上海科学技术出版社,1990.
[23] A C Waltham. Ground Subsidence [M].Blackie&Son Limited, London, 1989
[24] Burbey.T.J, D.C.Helm.Modeling three-dimensional deformation in response to pumping of unconsolidated aquifers [J].Environmental and environmental 5, 1999, no.2:199-212
[25] Jiang Li.A nonlinear elastic solution for subsidence due to ASR applications to multi-aquifer systems [A].In: edited by Laura Carbognin, Giuseppe Gambolati A. Ivan Johnson. Proceedings of the Sixth International Symposium on Land Subsidence[C].Ravenna, Italy, 2000
[26] 吴林高,缪俊发.抽灌水作用下土层变形及应力应变本构律的研究[J].地球科学,1995(5):581-588
[27] JiangLi, Donald C Helm.A nonlinear viscous model for aquifer compression associated with ASR applications [A].In: Edited by Laura Carbognin, Giuseppe Gambolati, and A. Ivan Johnson. Proceedings of the Sixth Intemational Symposium on Land Subsidence [C], Ravenna, Italy, 2000
[28] 冉启全,顾小芸.考虑流变特性的流固耦合地面沉降计算模型[J].中国地质灾害与防治学报,1998.9(2).
[29] X Y Gu, Q Q Ran. 3-D coupled with consideration of theological properties [A]. In: edited by Laura Carbognin, Giuseppe Gambolati, A. Ivan Johnson. Proceedings of the Sixth International Symposium on Land Subsidence [C], Ravenna, Italy, 2000
[30] Terzaghi K. Theoretical soil mechanics [M].Wiley, New York, 1943
[31] G Gambolati, R Allan Frezze. Mathematical simulation of the subsidence of Venice [J]. Theory Water Resource Research,1973,9(3):721-733
[32] 崔小东.Modflow和IDP在天津地面沉降数值计算中的应用与开发[J].中国地质灾害与防治学报,1998,9(2):122-128
[33] T R Shearer. A numerical model to calculate land subsidence, applied at Han gu in China [J]. Engineering Geology, 1998, 49(2):85-93
[34] R W Lewis, B.Schrefler.A fully coupled consolidation model of the subsidence of Venice[J].Water Resource Research, 1978, 14(2).
[35] J P His, J P Carter, J C Small.Surface subsidence and drawdown of the watertable due to pumping [J]. Geotech—nique44.1994(3):381-396
[36] Miau B Su, C L Su, C J Chang, et al.A numerical model of ground deformation induced by single well pumping[J]. Computer and Geotechnics, 1998(2):39-60.
[37] 张克绪等.开采地下水引起地面变形的分析[J].自然灾害学报,1996,5(4):50-591
[38] Biot M A. General theory of three dimensional consolidation [J].J.Appl.Phys, 1941, 12:155～164.
[39] Biot M A. General solution of the equation of elasticity and consolidation for a porous material [J]. J.Appl.Phys. 1956, 78:91～96
[40] Chai J C, Shen S L, Zhu H H, Zhang X L. Land subsidence due to groundwater drawdown in Shanghai [J]. Geotechnique, 2004, 54(2): 143～148.
[41] Shen S L, Tohno I, Nishigaki M, Miura N. 3D-analysis of land subsidence due to withdrawal of groundwater [A]. Geotechnical Engineering in Urban Construction [C].Beijing: Tsinghua University Press, 2003.431-436.
[42] Gambolati G, Monaco A, Galeati G, et al. New approaches and applications in subsurface flow modeling: 3D finite element analysis of dewateringfor an electro-nuclear plan [A]. In: Groundwater Flow and Quality Modeling[C].Boston: D. Reidel Publishing Co., 1998.11～15.
[43] Cassiani G, Palozzo W, Zoccatelli C. Simulation of Fig.8 Steady flow system with an impervious wall subsidence caused by gas production from an offshore field and comparison with field data [J]. Journal of Petroleum Science and Engineering, 2002, 45(1): 123～134.
[44] 上海市环境地质站,比利时国家科技部.上海市地面沉降教学模型研究[R].1989
[45] 朱锡冰,奚建国等.上海浦东新区“地下水水量—水位—沉降联合数学模型应用”[R].1995
[46] 方正,李勤奋等.上海市区域地下水调查[R].1999
[47] 陈崇希,裴顺平.地下水开采—地面沉降模型研究[J].水文地质工程地质,2001,2:5-8
[48] 方鹏飞,朱向荣.太原市地面沉降的计算与预测[J].煤田地质与勘探,2002,30(4):44-46
[49] 崔亚莉,邵景力,谢振华等.基于MODFLOW的地面沉降模型研究—以北京市区为例[J].工程勘察,2003,5:19-22
[50] 张利生,朱国荣.鹏山水源地开采条件下的含水层固结模型[J].岩土力学,2005,26(7):1141-1147
[51] 吴振样.地下热水三维有限元数值模拟与地面沉降研究[D],2004
[52] 卢立志,隋兆显等.江苏省沿江工业走廊水文地质、工程地质、环境地质综合评价报告[R].1989.
[53] Jorn Hoffmann, S.A.Leake, etc. MODFLOW-2000 Ground-Water Model—User Guide to the Subsidence and Aquifer-System CompactionPackage[R].2000
[54] 钱家欢,殷宗泽等.土工计算与原理[M]..北京:中国水利水电出版社,2003
[55] 朱学愚.地下水资源评价[M].南京:南京大学出版社,1987
[56] 骆祖江,武永霞.盐城市地下水资源规划评价三维数值模型[J].水资源保护,2005,21(5):37-41
[57] Terzaghi K. Theoretical soil mechanics [M].Wiley, New York, 1943
[58] 罗国煜,李生林.工程地质学基础[M].南京:南京大学出版社,1990.
[59] 殷宗泽.土体的沉降与固结[M].北京:中国电力出版社,1998.
[60] Poland, J.F., and Davis, G.H. Land subsidence due to withdrawals of fluids, in Varnes, D.J., and Kiersch, G., eds.: Geological Society of America Reviews in Engineering Geology, v. 2, p. 187-269.
[61] 薛禹群,朱学愚等.地下水动力学[M].北京:地质出版社,1997:5.
[62] Wen Hsing Chiang, Wolfgang Kinzelbach.3D Groundwater modeling with PMWIN—A Simulation System for Modeling Groundwater Flow and Pollution[M]. Springer, 2001.
[63] Jorgensen, D.G., 1980, Relationships between basic soils-engineering equations and basic ground-waterflow equations:[R].U.S. Geological Survey Water-Supply Paper 2064, 40 p.
[64] Johnson, A.I., Moston, R.P., and Morris, D.A., 1968, Physical and hydrologic properties of water-bearing deposits in subsiding areas in California[R].U.S. Geological Survey Professional Paper 497-A, 71 p.
[65] Riley, F.S., 1969, Analysis of borehole extensometer data from central California: [A].International Association of Scientific Hydrology Publication 89, p. 423-431.
[66] Poland, J.F., and Ireland, R.L., 1988, Land subsidence in the Santa Clara Valley, California, as of 1982: [A] U.S. Geological Survey Professional Paper 497-F, 61 p.
[67] Heywood, C.E., 1997, Piezometric-extensometric estimations of specific storage in the Albuquerque Basin, New Mexico, in Prince, K.R. and Leake, S.A., eds., U.S. Geological Survey Subsidence lnterest Group Conference: Proceedings of the technical meeting, Las Vegas, Nevada, February 14-16, [A]. 1995: U.S. Geological Survey Open-File Report 97-47, p. 21-26
[68] Joseph F. Poland. Mechanics of aquifer systems—The scientific legacy of, Land subsidence case studies and current research: [A].Proceedings of the Dr. Joseph F. Poland Symposium on Land Subsidence, Association of Engineering Geologists Special Publication 1998,no. 8, p. 13-27.
[69] A C Waltham. Ground Subsidence[M].Blackie&Son Limited, London, 1989.
[70] T.J.Burbey. Stress-strain analyses for aquifer-system characterization [J].Ground Water, 2001, 39(1):128-136.
[71] T J.Berbey.Effects of horizontal strain in estimating specific storage and compaction in confined and leaky aquifer systems [J]. Hydrogeology Journal, 1999, (7)
[72] 陈崇希,唐仲华.地下水流动问题数值方法[M].武汉:中国地质大学出版社,1994.
[73] 孙讷正.地下水流的数值模型和数值方法[M].北京:地质出版社,1981.
[74] 孙峰根,王心田.水文地质计算的数值方法[M].徐州:中国矿业大学出版社,1995.
[75] 张蔚榛.地下水非稳定流计算和地下水资源评价[M].北京:科学出版社,1983.
[76] Leake, S.A., and Prudic, D.E., Documentation of a computer program to simulate aquifer-systemcompactionusingthemodular finite-differenceground-water flow model[J].U.S. Geological Survey Techniques of Water-Resources Investigations, book 6, 1991,chap.A2,68 p.
[77] 陈玉田.偏微分方程数值解法[M].南京:河海大学出版社,1999
[78] 顾尔祚.流体力学中的有限差分法基础[M].上海:上海交通大学出版社,1988
[79] 蒋亚萍,陈余道.MODFLOW— 一套水文地质学实用计算软件[J].广西地质,1999,12(3):75-78
[80] User's manual for Visual MODFLOW. [R].Waterloo Hydrogeologic Inc. 1999
[81] 谭浩强,田淑清.FORTRAN语言程序设计[M].北京:高等教育出版社,1997
[82] 马瑞民,衣治安.FORTRAN90程序设计[M].哈尔滨工程大学出版社.1998
[83] 范啸涛,季光明.预优矩阵及其构造技术[J].成都理工大学学报(自然科学版),2003.30(4):P432-435
[84] 李裕伟,任效颖等..二维地质统计学[M].北京:地质出版社,1994.
[85] 李国敏,陈崇希.克立格法在地下水数值模型中的应用[J].地质科技情报,1995,14(3):95-99
[86] MITASOVA H, MITAS L, BROWN W M, et al. Modeling spatially and temporally distributed phenomena: new method and tools for GRASS GIS[J].IntJGIS, 1995(4):433-446.
[87] 楚明敏,王彩会等.苏州城市规划区浅层地下水开发利用研究报告[R].2003
[88] 郭坤一,骆祖江等.长江三角洲地区地下水资源与地质灾害调查评价[R].2003
[89] 于军,宗开红等.基于GMS平台的苏锡常地区地下水流模型研究报告[R].2003
[90] 祝晓彬,吴吉春等.长江三角洲(长江以南)地区深层地下水三维数值模拟[J].地理科学,2005,25(1):68-73
[91] 姜蓓蕾,吴吉春等.常州市滆湖典型地块浅层地下水数值模拟[J].水文地质工程地质,2004(6):29-32
[92] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第一季度)[R].2003
[93] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第二季度)[R].2003
[94] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第三季度)[R].2003
[95] 江苏省常州市水文水资源勘测局.常州市(城区、金坛、溧阳)地下水水情报告(2003年第四季度)[R].2004
[96] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第一季度)[R].2003
[97] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第二季度)[R].2003
[98] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第三季度)[R].2003
[99] 江苏省南通市水利局.南通市深层地下水监测报告(2003年第四季度)[R].2004
[100] 江苏省扬州市水利局.扬州市地下水动态监测汇总表(2003年度)[R].2004
[101] 国家技术监督局.中华人民共和国国家标准(GBT14497-93)地下水资源管理模型工作技术要求[S].1994,3:1～15