非稳定初始面对自由渗流影响规律的研究
|
摘要
河渠渗流、土坝渗流、边坡渗流等皆属于有自由面的无压渗流问题。这类问题的自由面位置是待求的，事先未知，如何简捷、有效地确定自由面的位置历来都是无压渗流场分析的关键。自由面渗流分析方法很多，归纳起来主要有解析法、数值法和模拟法等三大类；在自由面渗流规律研究过程中，往往依赖解析法。关于自由面渗流规律的研究，现有文献给出的模型，其初始自由面都是稳态的；以河渠一半无限含水层中的潜水渗流模型为例，初始自由面不仅是稳态的、而且往往是水平的，这既不符合大多数实际情况，也不具有普遍性。
本文以河渠控制下的半无限潜水含水层中的潜水渗流模型为例，构造初始面非稳定条件下的一维潜水非稳定渗流模型；利用Boussinesq第一线性化方法，通过Laplace变换，再结合迭加原理，导出模型的解析解；并以其为基础，讨论模型中参数的计算方法、研究非稳定初始面对自由渗流的影响规律。
研究表明，文中建立的潜水非稳定渗流模型合理，所给出解析解考虑了非稳定初始面因素，计算的自由面位置比原潜水渗流模型的解更接近观测值；规律研究显示非稳定初始面对自由渗流的影响与导压系数(α)、时间因子(T)和时间(t)成反比，与河渠原水位和最初初始面水位之差(△H_1)及距河渠距离(X)成正比。
The seepage of waterway, soil-dam and side-slope are belongs to not press groundwater flow problems with free surfaces. The free surfaces' position is unknown and requested. How to calculate the free surfaces' position briefly and effectively is a sticking point. The analysis methods of flow with free surface can divide into three major type: analytical method, mathematical method and simulation method. The analytical method used currently in free seepage regulation researching. Presently research of free surface flow regulation always use steady beginning free surface. Take the potential flow model of waterway- semi-infinity porous medium for example, all now study are all use steady beginning free surface, and is usually horizontal. This is obviously not agreed with actual circumstances, also not have the theories meaning of catholicity.
Take the potential flow model of waterway- semi- infinity porous medium for example, the one-dimensional seepage model of a infinite aquifer have not perpendicularity seepage with unsteady free surface is described. Based on the first linearized Boussnesq equation, the analytical solution of the model mentioned above can be obtained by Laplace transform combined with the principle of superposition. Base on the analytical solution, the account method of parameters in the model and the affect disciplinary of unsteady free surface are analyzed.
The result shows that the parameter-computation method proposed by the paper is reasonable. Comparing with the original solution, the new analytical solution subject to unsteady beginning-surface is more accurately approaches to the observation value. Regulation research show that the affection of unsteady beginning-surface in inverse proportion to press transmit coefficient (α)、time gene(T) and the time(t),but in direct proportion to the water level difference(△H_1) between the waterway and prime beginning-surface and the distance to waterway(x).
本文以河渠控制下的半无限潜水含水层中的潜水渗流模型为例，构造初始面非稳定条件下的一维潜水非稳定渗流模型；利用Boussinesq第一线性化方法，通过Laplace变换，再结合迭加原理，导出模型的解析解；并以其为基础，讨论模型中参数的计算方法、研究非稳定初始面对自由渗流的影响规律。
研究表明，文中建立的潜水非稳定渗流模型合理，所给出解析解考虑了非稳定初始面因素，计算的自由面位置比原潜水渗流模型的解更接近观测值；规律研究显示非稳定初始面对自由渗流的影响与导压系数(α)、时间因子(T)和时间(t)成反比，与河渠原水位和最初初始面水位之差(△H_1)及距河渠距离(X)成正比。
The seepage of waterway, soil-dam and side-slope are belongs to not press groundwater flow problems with free surfaces. The free surfaces' position is unknown and requested. How to calculate the free surfaces' position briefly and effectively is a sticking point. The analysis methods of flow with free surface can divide into three major type: analytical method, mathematical method and simulation method. The analytical method used currently in free seepage regulation researching. Presently research of free surface flow regulation always use steady beginning free surface. Take the potential flow model of waterway- semi-infinity porous medium for example, all now study are all use steady beginning free surface, and is usually horizontal. This is obviously not agreed with actual circumstances, also not have the theories meaning of catholicity.
Take the potential flow model of waterway- semi- infinity porous medium for example, the one-dimensional seepage model of a infinite aquifer have not perpendicularity seepage with unsteady free surface is described. Based on the first linearized Boussnesq equation, the analytical solution of the model mentioned above can be obtained by Laplace transform combined with the principle of superposition. Base on the analytical solution, the account method of parameters in the model and the affect disciplinary of unsteady free surface are analyzed.
The result shows that the parameter-computation method proposed by the paper is reasonable. Comparing with the original solution, the new analytical solution subject to unsteady beginning-surface is more accurately approaches to the observation value. Regulation research show that the affection of unsteady beginning-surface in inverse proportion to press transmit coefficient (α)、time gene(T) and the time(t),but in direct proportion to the water level difference(△H_1) between the waterway and prime beginning-surface and the distance to waterway(x).
引文
[1]孔祥言．高等渗流力学[M]．合肥：中国科学技术大学出版社，1999
[2]雅．贝尔(Bear J．)著，许涓铭等译．地下水水力学[M]．北京：地质出版社，1985
[3]薛禹群主编．地下水动力学原理[M]．北京：地质出版社，1986
[4]Darcy H．Les fontaines publiques de la ville de Dijon[M]．Paris：Victor Dalmont，1856
[5]陶月赞，姚梅．地下水渗流力学的发展进程与动向[J]．吉林大学学报(地球科学版)，2007(2)：221-230
[6]A.J.E.J.Dupuit.Etudes theoriques et practiques sur le mouvement des aux dans les canaux decouverts et a travers les terrains permeables[M].2nd ed.:Dunod,Paris,1863
[7]J.Boussinesq.Recherches theoriques sur l'ecoulement des.nappes d'eau infiltrees dans le sol et sur le debit des sources[J].Journal de Math.Pures Appl.1904,5(10):5-78
[8]A.Theim.,Verfahress for Naturlicher Groundwater geschwindegkiten[J].Polyt.Notizblatt,1887,42:229
[9]P.Forchheimer,Uber die Ergebigkeit von Brunnen Anlagen und Sickerschlitzen[J].Zeitschrift des Architekten-und Ingenieur Vereins zu Hannover.1886,32:539-564
[10]Thomas E.R.A Brief History of Contribution to Ground Water Hydrogeology by U.S.Geological Survey[J].Ground Water,2004,42(2):625-631
[11]Oscar Edward Meinzer.Compressibility and elasticity of artesian aquifers[J].Economic Geology,1928,23(3):263-291
[12]张宏仁等．地下水非稳定流理论的发展和应用[M]．地质出版社，1975
[13]Theis,C.V.The relation between the lowering of the piezometric surface and the rate and discharge of a well using ground water storage[J].Transactions of the American Geophysical Union,1935,16:519-524
[14]Slichter,Charles S.Theoretical investigation of the motion of ground water[R].U.S.Geological Survey Nineteenth Annual Report,1899,Part 2:295-384
[15]李佩成．地下水动力学[M]．北京：农业出版社，1993
[16]Hunt R．J．Ground Water Modeling Applications Using the Analytic Element Method[J]Ground Water，2006，44(1)：5-15
[17]郭尚平，刘慈群．渗流力学(中国大百科全书在线)，http：／／www．baicle．com：8080／cp／resource／articles／451
[18]杜延龄，许国安．渗流分析的有限元法和电网络法[M]．北京：水利电力出 版社，1992
[19]吴林高，缪俊发，张瑞等．渗流力学[M]．上海：上海科学技术文献出版社，1996
[20]丁留谦，许国安．三维非达西渗流的有限元分析[J]．水利学报，1990，(10)：49-54
[21]朱伯芳．有限单元法原理与应用(第二版)[M]．北京：中国水利水电出版社，1998：247-269
[22]Hunt R．J．Ground Water Modeling Applications Using the Analytic Element Method[J]．Ground Water，2006，44(1)：5-15
[23]陶月赞．河渠控制下潜水含水层中潜水蒸发规律与计算[J]．自然科学进展，2006，16(5)：590-595
[24]陶月赞，曹彭强，席道瑛．垂向入渗与河渠边界影响下潜水非稳定流参数的求解[J]水利学报，2006，37(8)：913-917
[25]Tao Y．Z．，Xi D．Y．Rule of transient phreatic flow subjected to vertical and horizontal seepage[J]．Applied mathematics and mechanics，2006，27(1)：53-59
[26]顾慰慈．渗流计算原理及应用[M]．北京：中国建材工业出版社，2000
[27]关明芳，陈洪凯．渗流自由面求解方法综述[J]．重庆交通学院学报，2005，24(5)：68-73
[28]薛禹群，谢春红．水文地质学的数值法[M]．北京：煤炭工业出版社，1980
[29]陈崇希，唐仲华．地下水流动问题数值方法[M]．武汉：中国地质大学出版社，1990
[30]苑莲菊，李振栓，武胜忠等．工程渗流力学及应用[M]．北京：中国建材工业出版社，2001
[31]介玉新，揭冠周，李广信．用适体坐标变换方法求解渗流[J]．岩土工程学报，2004，26(1)：52-56
[32]Min-Ho Koo,Darrell I.Leap.Modeling Three-Dimensional Groundwater Flows[J].Transport in Porous Media,1998,30:345-362
[33]Zienkiewicz O.C.,Mayer P.Solution of anisotropic seepage by Finite Elements.engMech Div,ASCE,1966,92,EMI:111-120
[34]France,P.W.etal,Numerical Analysis of Free Surface Seepage Problems [J]Irrig.Drain Div,ASCE,1971,97-112
[35]RAbier Sylvain,Medale Mare.Computation of free surface flow with a projection FEM in a moving mesh frame work.[J]Computer Methods in Applied Mechanics and Engineering,2003,10(192):4703-4721
[36]魏泽光，黄俊，许国安等．三维稳定渗流的有限元计算[J]．水利学报，1982， (1)：39-47
[37]李春华．稳定渗流的有限元计算和决定自由面表面位置的一种方法[M]．水利水电科学院论文集第8集(岩土工程)，北京：水利电力出版社，1982：148-156
[38]张有天．用边界元求解有排水孔的渗流场[J]．水利学报，1982，(7)：36-40
[39]杜延龄，许国安．渗流分析的有限元法和电网络法[M]．北京：水利电力出版社，1992
[40]陈惠远．土石坝有限元分析[M]．南京：河海大学出版社，1988：1-119．
[41]崔俊芝，梁俊．现代有限元软件方法[M]．北京：国防工业出版社，1995
[42]王勖成，邵敏．有限单元法基本原理和数值方法(第二版)[M]．北京：清华大学出版社，1997．
[43]黄蔚，刘迎曦，周承芳．三维无压渗流场的有限元算法研究[J]．水利学报，2001，(6)：33-36
[44]陈士俊等．工程渗流有限元计算方法研究[J]．人民黄河，2000，22(9)：34-36
[45]杜延龄，许国安．复杂岩基三维渗流有限元分析研究[J]．水利学报，1991，(7)：19-26
[46]Neuman S.P.Saturate-unsaturated seepage by finite element[J].Hydraulic Div,ACSE,1973,99(HY12):2233-2250
[47]Oden J.T.,Kikuchi N.Rencent advances:Thoory of veriational inequalities with application to problema of flow through porous media[J].International Journal of Engineering Science,1980(18):173-184.
[48]Bruch J.C.A survey of free boundary value problems in the theory of fluid flow through porous media:Variational inequality approach[J].Variational inequality approach.ADV.Water res,1980(3):65-80
[49]Desai C.S.Finite element residual schemes for uneonfitmed flow [J]Int Num Method Eng.1976,10(6):1415-1418
[50]Desai C.S.,Li G.C.A residual flow procedure and application for free surface flow in porous media[J]Advances in Water Resources,1983(6):27-35
[51]Bathe K.J.,Khoshgoftan M R.Finite element free surface seepage analysis without mesh itemfion[J].International Joumal for Numerical and Analytic Methods in Geomechanic,1979(3):13-22
[52]李春华．稳定渗流有限元计算时采用固定网格法的初步研究[D]．第三届全国渗流力学学术讨论会论文汇编(3)，长江科学院，1986
[53]张有天，陈平，王镭．有自由面渗流分析的初流量法[J]．水利学报，1988，(8)：18-26
[54]王媛．求解有自由面渗流问题的初流量法的改进[J]．水利学报，1998，(3)： 68-73
[55]吴梦喜，张学勤．有自由面渗流数值分析的虚单元法[J]．水利学报，1994，(8)：67-71
[56]彭华，罗谷华．有自由面渗流分析的子单元及其应用[J]．水电站设计，1996，12(4)：24-27
[57]彭华，曹定胜，朱劲大等．有自由面渗流分析的虚单元法[J]．中国农村水利水电(农田水利与小水电)，1997，(3)：26-28
[58]彭华，曹定胜等．有自由面渗流分析的弃单元法及其应用[J]．人民长江，1997，28(8)：38-40
[59]凌道盛．有自由面渗流分析的虚节点法[J]．浙江大学学报(工学版)，2002，36(3)：243-246
[60]王均星，吴雅峰，白呈富．有自由面渗流分析的流形单元法[J]．水电能源科学，2003，21(4)：23-26
[61]王贤能，黄润秋．有自由面渗流分析的高斯点法[J]．水文地质工程地质，1997，(6)：1-4
[62]邓建霞．有自由面渗流分析的加密高斯点单元传导矩阵调整法[D]四川大学，2004
[63]周创兵，熊文林，梁业国．求解无压渗流场的一种新方法[J]．水动力学研究与进展，1996，11(5)：528-534
[64]速宝玉，沈振中，赵坚．用变分不等式理论求解渗流问题的截止负压法[J]．水利学报，1996，(3)：22-29
[65]李广信，葛镜宏，介玉新．有自由面渗流的无单元法[J]．清华大学学报(自然科学版)，2002，42(4)：1552-1553
[66]葛锦宏，李广信，介玉新．无单元法在有自由面渗流计算中的应用[J]．计算力学学报，2003，20(2)：241-243
[67]王均星，吴雅峰，百呈富．有自由面渗流分析的流形单元法[J]．水电能源科学，2003，21(4)：23-25
[68]朱军，刘光廷．改进的单元渗透矩阵调整法求解无压渗流场[J]．水利学报，2001，(8)：49-52
[69]徐书平，喻国安，扬剑飞．渗流分析边界条件处理的算法实现[J]．土工基础，2002，16(2)：51-54
[70]朱岳明，龚道勇．三维饱和—非饱和渗流场求解及其逸出面边界条件处理[J]．水科学进展，2003，14(1)：67-71
[71]金生，耿艳芬，王志力．利用饱和—非饱和渗流模型计算坝体自由面渗流[J]．大连理工大学学报，2004，44(1)：110-113
[72]陈崇希，林敏．地下水动力学[M]．武汉：中国地质大学出版社，1999
[2]雅．贝尔(Bear J．)著，许涓铭等译．地下水水力学[M]．北京：地质出版社，1985
[3]薛禹群主编．地下水动力学原理[M]．北京：地质出版社，1986
[4]Darcy H．Les fontaines publiques de la ville de Dijon[M]．Paris：Victor Dalmont，1856
[5]陶月赞，姚梅．地下水渗流力学的发展进程与动向[J]．吉林大学学报(地球科学版)，2007(2)：221-230
[6]A.J.E.J.Dupuit.Etudes theoriques et practiques sur le mouvement des aux dans les canaux decouverts et a travers les terrains permeables[M].2nd ed.:Dunod,Paris,1863
[7]J.Boussinesq.Recherches theoriques sur l'ecoulement des.nappes d'eau infiltrees dans le sol et sur le debit des sources[J].Journal de Math.Pures Appl.1904,5(10):5-78
[8]A.Theim.,Verfahress for Naturlicher Groundwater geschwindegkiten[J].Polyt.Notizblatt,1887,42:229
[9]P.Forchheimer,Uber die Ergebigkeit von Brunnen Anlagen und Sickerschlitzen[J].Zeitschrift des Architekten-und Ingenieur Vereins zu Hannover.1886,32:539-564
[10]Thomas E.R.A Brief History of Contribution to Ground Water Hydrogeology by U.S.Geological Survey[J].Ground Water,2004,42(2):625-631
[11]Oscar Edward Meinzer.Compressibility and elasticity of artesian aquifers[J].Economic Geology,1928,23(3):263-291
[12]张宏仁等．地下水非稳定流理论的发展和应用[M]．地质出版社，1975
[13]Theis,C.V.The relation between the lowering of the piezometric surface and the rate and discharge of a well using ground water storage[J].Transactions of the American Geophysical Union,1935,16:519-524
[14]Slichter,Charles S.Theoretical investigation of the motion of ground water[R].U.S.Geological Survey Nineteenth Annual Report,1899,Part 2:295-384
[15]李佩成．地下水动力学[M]．北京：农业出版社，1993
[16]Hunt R．J．Ground Water Modeling Applications Using the Analytic Element Method[J]Ground Water，2006，44(1)：5-15
[17]郭尚平，刘慈群．渗流力学(中国大百科全书在线)，http：／／www．baicle．com：8080／cp／resource／articles／451
[18]杜延龄，许国安．渗流分析的有限元法和电网络法[M]．北京：水利电力出 版社，1992
[19]吴林高，缪俊发，张瑞等．渗流力学[M]．上海：上海科学技术文献出版社，1996
[20]丁留谦，许国安．三维非达西渗流的有限元分析[J]．水利学报，1990，(10)：49-54
[21]朱伯芳．有限单元法原理与应用(第二版)[M]．北京：中国水利水电出版社，1998：247-269
[22]Hunt R．J．Ground Water Modeling Applications Using the Analytic Element Method[J]．Ground Water，2006，44(1)：5-15
[23]陶月赞．河渠控制下潜水含水层中潜水蒸发规律与计算[J]．自然科学进展，2006，16(5)：590-595
[24]陶月赞，曹彭强，席道瑛．垂向入渗与河渠边界影响下潜水非稳定流参数的求解[J]水利学报，2006，37(8)：913-917
[25]Tao Y．Z．，Xi D．Y．Rule of transient phreatic flow subjected to vertical and horizontal seepage[J]．Applied mathematics and mechanics，2006，27(1)：53-59
[26]顾慰慈．渗流计算原理及应用[M]．北京：中国建材工业出版社，2000
[27]关明芳，陈洪凯．渗流自由面求解方法综述[J]．重庆交通学院学报，2005，24(5)：68-73
[28]薛禹群，谢春红．水文地质学的数值法[M]．北京：煤炭工业出版社，1980
[29]陈崇希，唐仲华．地下水流动问题数值方法[M]．武汉：中国地质大学出版社，1990
[30]苑莲菊，李振栓，武胜忠等．工程渗流力学及应用[M]．北京：中国建材工业出版社，2001
[31]介玉新，揭冠周，李广信．用适体坐标变换方法求解渗流[J]．岩土工程学报，2004，26(1)：52-56
[32]Min-Ho Koo,Darrell I.Leap.Modeling Three-Dimensional Groundwater Flows[J].Transport in Porous Media,1998,30:345-362
[33]Zienkiewicz O.C.,Mayer P.Solution of anisotropic seepage by Finite Elements.engMech Div,ASCE,1966,92,EMI:111-120
[34]France,P.W.etal,Numerical Analysis of Free Surface Seepage Problems [J]Irrig.Drain Div,ASCE,1971,97-112
[35]RAbier Sylvain,Medale Mare.Computation of free surface flow with a projection FEM in a moving mesh frame work.[J]Computer Methods in Applied Mechanics and Engineering,2003,10(192):4703-4721
[36]魏泽光，黄俊，许国安等．三维稳定渗流的有限元计算[J]．水利学报，1982， (1)：39-47
[37]李春华．稳定渗流的有限元计算和决定自由面表面位置的一种方法[M]．水利水电科学院论文集第8集(岩土工程)，北京：水利电力出版社，1982：148-156
[38]张有天．用边界元求解有排水孔的渗流场[J]．水利学报，1982，(7)：36-40
[39]杜延龄，许国安．渗流分析的有限元法和电网络法[M]．北京：水利电力出版社，1992
[40]陈惠远．土石坝有限元分析[M]．南京：河海大学出版社，1988：1-119．
[41]崔俊芝，梁俊．现代有限元软件方法[M]．北京：国防工业出版社，1995
[42]王勖成，邵敏．有限单元法基本原理和数值方法(第二版)[M]．北京：清华大学出版社，1997．
[43]黄蔚，刘迎曦，周承芳．三维无压渗流场的有限元算法研究[J]．水利学报，2001，(6)：33-36
[44]陈士俊等．工程渗流有限元计算方法研究[J]．人民黄河，2000，22(9)：34-36
[45]杜延龄，许国安．复杂岩基三维渗流有限元分析研究[J]．水利学报，1991，(7)：19-26
[46]Neuman S.P.Saturate-unsaturated seepage by finite element[J].Hydraulic Div,ACSE,1973,99(HY12):2233-2250
[47]Oden J.T.,Kikuchi N.Rencent advances:Thoory of veriational inequalities with application to problema of flow through porous media[J].International Journal of Engineering Science,1980(18):173-184.
[48]Bruch J.C.A survey of free boundary value problems in the theory of fluid flow through porous media:Variational inequality approach[J].Variational inequality approach.ADV.Water res,1980(3):65-80
[49]Desai C.S.Finite element residual schemes for uneonfitmed flow [J]Int Num Method Eng.1976,10(6):1415-1418
[50]Desai C.S.,Li G.C.A residual flow procedure and application for free surface flow in porous media[J]Advances in Water Resources,1983(6):27-35
[51]Bathe K.J.,Khoshgoftan M R.Finite element free surface seepage analysis without mesh itemfion[J].International Joumal for Numerical and Analytic Methods in Geomechanic,1979(3):13-22
[52]李春华．稳定渗流有限元计算时采用固定网格法的初步研究[D]．第三届全国渗流力学学术讨论会论文汇编(3)，长江科学院，1986
[53]张有天，陈平，王镭．有自由面渗流分析的初流量法[J]．水利学报，1988，(8)：18-26
[54]王媛．求解有自由面渗流问题的初流量法的改进[J]．水利学报，1998，(3)： 68-73
[55]吴梦喜，张学勤．有自由面渗流数值分析的虚单元法[J]．水利学报，1994，(8)：67-71
[56]彭华，罗谷华．有自由面渗流分析的子单元及其应用[J]．水电站设计，1996，12(4)：24-27
[57]彭华，曹定胜，朱劲大等．有自由面渗流分析的虚单元法[J]．中国农村水利水电(农田水利与小水电)，1997，(3)：26-28
[58]彭华，曹定胜等．有自由面渗流分析的弃单元法及其应用[J]．人民长江，1997，28(8)：38-40
[59]凌道盛．有自由面渗流分析的虚节点法[J]．浙江大学学报(工学版)，2002，36(3)：243-246
[60]王均星，吴雅峰，白呈富．有自由面渗流分析的流形单元法[J]．水电能源科学，2003，21(4)：23-26
[61]王贤能，黄润秋．有自由面渗流分析的高斯点法[J]．水文地质工程地质，1997，(6)：1-4
[62]邓建霞．有自由面渗流分析的加密高斯点单元传导矩阵调整法[D]四川大学，2004
[63]周创兵，熊文林，梁业国．求解无压渗流场的一种新方法[J]．水动力学研究与进展，1996，11(5)：528-534
[64]速宝玉，沈振中，赵坚．用变分不等式理论求解渗流问题的截止负压法[J]．水利学报，1996，(3)：22-29
[65]李广信，葛镜宏，介玉新．有自由面渗流的无单元法[J]．清华大学学报(自然科学版)，2002，42(4)：1552-1553
[66]葛锦宏，李广信，介玉新．无单元法在有自由面渗流计算中的应用[J]．计算力学学报，2003，20(2)：241-243
[67]王均星，吴雅峰，百呈富．有自由面渗流分析的流形单元法[J]．水电能源科学，2003，21(4)：23-25
[68]朱军，刘光廷．改进的单元渗透矩阵调整法求解无压渗流场[J]．水利学报，2001，(8)：49-52
[69]徐书平，喻国安，扬剑飞．渗流分析边界条件处理的算法实现[J]．土工基础，2002，16(2)：51-54
[70]朱岳明，龚道勇．三维饱和—非饱和渗流场求解及其逸出面边界条件处理[J]．水科学进展，2003，14(1)：67-71
[71]金生，耿艳芬，王志力．利用饱和—非饱和渗流模型计算坝体自由面渗流[J]．大连理工大学学报，2004，44(1)：110-113
[72]陈崇希，林敏．地下水动力学[M]．武汉：中国地质大学出版社，1999