基于FEFLOW的吉林西部地下水数值模拟研究
|
摘要
吉林省西部地区水资源缺乏,土地荒漠化、盐渍化现象严重,是典型的农牧交错带和生态环境脆弱区。本文以“吉林省西部地下水过度开发及其地表生态效应研究”为依托,简述了该区地质、水文地质条件,建立了水文地质概念模型及数学模型,应用FEFLOW软件对三维地下水流数值模型求解,并对未来15年地下水流场进行了预报。在水流模型的基础上进一步建立了吉林西部溶质运移模型,选取总氮为模拟对象。并应用该模型对未来10年地下水中总氮含量的演化趋势进行了预报。预测结果表明:区内地下水位逐年持续下降并有扩大趋势,潜水下降最明显处为洮儿河扇形地前缘和洮南、镇赉附近,其次为松拉河间地块;承压水只有在西部山前平原和低平原的过渡带水位下降值相对较大。研究区内氮浓度变化整体呈上升趋势,尤以洮南、洮北地段浓度升高幅度最大,前郭灌区次之。针对水量、水质预报结果分析原因并提出防治建议,为该区水资源的优化配置及地下水污染防治提供基本依据。
The paper which is based on the project“the research on groundwater overexploitation and its surface ecological effect in the West of Jilin Province”of Land Mineral Resources Management in Jilin Province has appraised and forecasted for groundwater resources in Western Jilin Province. The water resources are scarce in the west of Jilin Province, per-capita quantity of water resource is 809m3/a, which is 53.2% of the hole province (1520m3/a) and only 36.4% of China(2220m3/a). Land desertification and salinization are very serious in the research area, natural calamity such as the sandstorm takes place occasionally, so this area is a typical farming grazing transitional belt and vulnerable ecological environment area. In recent years, with the development of urban construction and national economy, water demand, the discharge of waste water and waste residue increased year by year, the contradiction of supply and requirements about water resources and the problems of environmental pollution become outstanding day by day, and have produced the inestimable negative effect on the ecological environment and human health. Considering these questions, this paper utilized the existing dynamic observation data of groundwater of the western region, spatial distribution materials of groundwater quality and some relevant data, set up the reasonable modeling of groundwater quantity and quality, and get clear understanding of exploitable resources amount of groundwater and groundwater environmental quality in the research area. The paper offers scientific basis for protecting the groundwater resources and realizing the sustainable development of regional water resource.
This text introduces numerical model, especially the current situation of domestic and international researches, development trends and existing main problems of the water quality model. On the basis of fully collecting and analyzing relative materials such as geology, hydrogeology and environment of the west of Jilin Province, the paper explaines systematically the cause, era and lithology of the main stratum in the research area, and analyzes the law of recharge, runoff, drainage and dynamic characteristics of the groundwater. After confirming the range of simulation district and object layer for calculation and generalizating hydraulic characteristics, vertical and lateral boundary, a groundwater mathematical model is set up, according to the hydrogeological conceptual model. And FEFLOW(Finite element subsurface flow system)which is based on finite element method is employed to solve the numerical equations.
Before data input, overlap parameter subarea and administrative region, then divide the source sink data of each administrative region into the recharge of sub area. This simulation used the water level observed value of selecting January of 2000 - January of 2002 to identify and inspect the model, the fitting effect is all right, and it explains that the determination of aquifer structure and hydrogeological parameter, the generalization boundary condition are reasonable, so the model can be used in forecast. Based on the model, the paper predicts the groundwater flow field of the following 15 years in the research area. The prediction results demonstrate that exploiting according to the predicted enlargement exploiting amount, the water tables of the phreatic water and the artesian water have the tendency to decline continuously year by year, comparing to the water table of the artesian water, the water table of the phreatic water declines in evidence, especially the water table of the advancing edge of Tao-er River Fan and Taonan and Zhenlai has the largest change, is 5-8m, the water table of phreatic water of plot among Songla River appears obvious declining too, the average decreasing amplitude in year is about 0.5m. The water table of the artesian water declines obviously in the center of the basin whose recharge condition is worse and western plain in front of mountain and transitional belt of low plain that the penetrability of aquifer is worse and the exploitation quantity is large, decreasing amplitude is 2-4m.
Based on the water quantity model, the groundwater solute transport model of the West Jilin is set up. Firstly, the Factor Analysis Method is used to analyze the main ion of groundwater chemical composition in the research area, the result shows that nitrate nitrogen is the main pollutant influenced by human activity. Ignoring the conversional mechanism between three nitrogen, we confirm the total nitrogen( NO3 ? ?N、NH 4+ ?N、NO 2? ?N) as simulation target of this research. Because the content of the total nitrogen in artesian water is low, and changes unconspicuously in many years, the simulation didn’t have actual meaning, so this paper gives the simulation and the forecast in total nitrogen of the artesian water. It chooses the water- chemical materials which keep the same time with flow model to identification and inspection, from the fitting picture of groundwater calculation and actual measurement concentration field, we can know that the fitting curve reflected the variation tendency, the precision of model agrees with the standard, so the modelcan be used to predict. The forecast time is ten years. The prediction results demonstrate that under the condition of area shape pollution because of using of nitrogenous fertilizer, the nitrogen concentration of the research area assume up tendency, especially in Baicheng and Taonan, because domestic wastewater, industrial sewage and solid waste that are piled up at will, the change of the nitrogen concentration is very obvious. Qianguo takes second place. In the other place, with the flow of the groundwater, nitrogen is discharge outside the area, and the concentration of the nitrogen does not change much on the whole.
The groundwater resources quantity of the research area is limited, and the intensity of development and utilization is high, so at the present time the water table drops obviously, which could not satisfy water need of long-term development of the whole area, so there are some suggestions that on the base of reasonable utilizing the groundwater resources, strengthen the building of hydraulic engineering of the river pass through, fully utilize surface water resource that pass through, and according to the allocation and management of surface water and groundwater, reach the long-range objective that can satisfy water demand and improve the fragile ecological environment of the research area. At the same time, reduce the artificial pollution of groundwater, control the discharge of waste water of industry and rural enterprise and the piling up of solid waste strictly, and reduce the pollution of the groundwater from industry and the continuation depravation of ecological environment. Improve the technology of digging well, especially the technology that agriculture well drilling, improve its close craft, reduce hydraulic interaction of the phreatic water and the artesian water, and control the pollution of the artesian water quality when the phreatic water is infiltrating. Then, satisfying growth in gear of crops as the precondition, improve the technology of the using of chemical fertilizer, and develop ecological agriculture and green agriculture to reduce the pollution to groundwater from pesticide and the using of chemical fertilizer.
The paper which is based on the project“the research on groundwater overexploitation and its surface ecological effect in the West of Jilin Province”of Land Mineral Resources Management in Jilin Province has appraised and forecasted for groundwater resources in Western Jilin Province. The water resources are scarce in the west of Jilin Province, per-capita quantity of water resource is 809m3/a, which is 53.2% of the hole province (1520m3/a) and only 36.4% of China(2220m3/a). Land desertification and salinization are very serious in the research area, natural calamity such as the sandstorm takes place occasionally, so this area is a typical farming grazing transitional belt and vulnerable ecological environment area. In recent years, with the development of urban construction and national economy, water demand, the discharge of waste water and waste residue increased year by year, the contradiction of supply and requirements about water resources and the problems of environmental pollution become outstanding day by day, and have produced the inestimable negative effect on the ecological environment and human health. Considering these questions, this paper utilized the existing dynamic observation data of groundwater of the western region, spatial distribution materials of groundwater quality and some relevant data, set up the reasonable modeling of groundwater quantity and quality, and get clear understanding of exploitable resources amount of groundwater and groundwater environmental quality in the research area. The paper offers scientific basis for protecting the groundwater resources and realizing the sustainable development of regional water resource.
This text introduces numerical model, especially the current situation of domestic and international researches, development trends and existing main problems of the water quality model. On the basis of fully collecting and analyzing relative materials such as geology, hydrogeology and environment of the west of Jilin Province, the paper explaines systematically the cause, era and lithology of the main stratum in the research area, and analyzes the law of recharge, runoff, drainage and dynamic characteristics of the groundwater. After confirming the range of simulation district and object layer for calculation and generalizating hydraulic characteristics, vertical and lateral boundary, a groundwater mathematical model is set up, according to the hydrogeological conceptual model. And FEFLOW(Finite element subsurface flow system)which is based on finite element method is employed to solve the numerical equations.
Before data input, overlap parameter subarea and administrative region, then divide the source sink data of each administrative region into the recharge of sub area. This simulation used the water level observed value of selecting January of 2000 - January of 2002 to identify and inspect the model, the fitting effect is all right, and it explains that the determination of aquifer structure and hydrogeological parameter, the generalization boundary condition are reasonable, so the model can be used in forecast. Based on the model, the paper predicts the groundwater flow field of the following 15 years in the research area. The prediction results demonstrate that exploiting according to the predicted enlargement exploiting amount, the water tables of the phreatic water and the artesian water have the tendency to decline continuously year by year, comparing to the water table of the artesian water, the water table of the phreatic water declines in evidence, especially the water table of the advancing edge of Tao-er River Fan and Taonan and Zhenlai has the largest change, is 5-8m, the water table of phreatic water of plot among Songla River appears obvious declining too, the average decreasing amplitude in year is about 0.5m. The water table of the artesian water declines obviously in the center of the basin whose recharge condition is worse and western plain in front of mountain and transitional belt of low plain that the penetrability of aquifer is worse and the exploitation quantity is large, decreasing amplitude is 2-4m.
Based on the water quantity model, the groundwater solute transport model of the West Jilin is set up. Firstly, the Factor Analysis Method is used to analyze the main ion of groundwater chemical composition in the research area, the result shows that nitrate nitrogen is the main pollutant influenced by human activity. Ignoring the conversional mechanism between three nitrogen, we confirm the total nitrogen( NO3 ? ?N、NH 4+ ?N、NO 2? ?N) as simulation target of this research. Because the content of the total nitrogen in artesian water is low, and changes unconspicuously in many years, the simulation didn’t have actual meaning, so this paper gives the simulation and the forecast in total nitrogen of the artesian water. It chooses the water- chemical materials which keep the same time with flow model to identification and inspection, from the fitting picture of groundwater calculation and actual measurement concentration field, we can know that the fitting curve reflected the variation tendency, the precision of model agrees with the standard, so the modelcan be used to predict. The forecast time is ten years. The prediction results demonstrate that under the condition of area shape pollution because of using of nitrogenous fertilizer, the nitrogen concentration of the research area assume up tendency, especially in Baicheng and Taonan, because domestic wastewater, industrial sewage and solid waste that are piled up at will, the change of the nitrogen concentration is very obvious. Qianguo takes second place. In the other place, with the flow of the groundwater, nitrogen is discharge outside the area, and the concentration of the nitrogen does not change much on the whole.
The groundwater resources quantity of the research area is limited, and the intensity of development and utilization is high, so at the present time the water table drops obviously, which could not satisfy water need of long-term development of the whole area, so there are some suggestions that on the base of reasonable utilizing the groundwater resources, strengthen the building of hydraulic engineering of the river pass through, fully utilize surface water resource that pass through, and according to the allocation and management of surface water and groundwater, reach the long-range objective that can satisfy water demand and improve the fragile ecological environment of the research area. At the same time, reduce the artificial pollution of groundwater, control the discharge of waste water of industry and rural enterprise and the piling up of solid waste strictly, and reduce the pollution of the groundwater from industry and the continuation depravation of ecological environment. Improve the technology of digging well, especially the technology that agriculture well drilling, improve its close craft, reduce hydraulic interaction of the phreatic water and the artesian water, and control the pollution of the artesian water quality when the phreatic water is infiltrating. Then, satisfying growth in gear of crops as the precondition, improve the technology of the using of chemical fertilizer, and develop ecological agriculture and green agriculture to reduce the pollution to groundwater from pesticide and the using of chemical fertilizer.
引文
[1] 盛连喜,刘长生,周道纬. 吉林生态与生态建设[M].长春:东北师范大学出版社,2001.
[2] 魏林宏,束龙仓,郝振纯. 地下水流数值模拟的研究现状和发展趋势[J].重庆大学学报(自然科学版),2000,(23)(增刊):50-52.
[3] Pinder G. F., J.D. Bredehoeft. Application of the digital computer for aquifer evaluation . Water Resource[J],1968,4(5):1069-1093.
[4] 薛禹群,谢春红. 水文地质学的数值法[M].北京:煤炭工业出版社,1979.
[5] 丁继红,周德亮,马生忠. 国内外地下水模拟软件的发展现状与趋势[J].勘察科学技术,2002,(1):37-42.
[6] 周念清,朱蓉,朱学愚. MODFLOW 在宿迁市地下水资源评价中的应用[J].水文地质与工程地质,2000,(6):9-13.
[7] 林学钰. 地下水资源的基本属性和我国水文地质科学的发展[J].地学前缘,2002,23(3)(增刊):50～52.
[8] 吴吉春,薛禹群,黄海等. 山西柳林局部区域运移二维数值模拟[J].水利学报.2001,(8):38-43.
[9] 杨天行. 求解对流占优地下水水质模型的有限单元高精度广义迎风单元均衡格式及应用[J].长春地质学院学报,1996,(2):196-201.
[10] Freg Pohll etal.Modeling ground water flow and radioactive transport in a fratured aquifer. Ground Water [J], 1999,37(5):770-784.
[11] Henk Haitjema eatl.Selecting MODFLOW cell sizes for accurate flow field. Ground Water [J], 2001,39(6):931-938.
[12] Hank E.Julian etal.Numerial simulation of natural gradient tracer experiment for the natural attenuation study flow and physical transport. Ground Water [J], 2001,39(4):534-545.
[13] Jonathan Levy and Eric E.Ludy.Uncertainty quantification for delineation of wellhead protection areas using the gauss hermite quadrature approach. Ground Water[J] , 2000,38(1):63-75.
[14] A.M.Ebraheem·S.Riad.P.Wycisk.A.M.Sefelnasr. A local-scale groundwater flowmodel for groundwater resources management in Dakhla Oasis, SW Egypt. Hydrogeology journal [J], 2004,(12):714-722.
[15] 蒋辉. 河南沁阳市地下水水质的数值模拟[J].水文地质工程地质,1996,(2):14-20.
[16] 段青梅. 西辽河平原三维地质地质建模及地下水数值模拟研究[D].中国地质大学,2006.5:14-15.
[17] 周志芳,李艳. 复杂岩体地下水运动问题的有限分析法[J].水科学进展,1997,(3):26-42.
[18] 岳梅. 淮南市浅层地下水水质污染预测模型探讨[J].地下水,2002,(2):77-79.
[19] 吴剑锋,朱学愚. 由 MODFLOW 浅谈地下水流数值模拟软件的发展趋势[J].工程勘察,2000,(2):12-15.
[20] 王惠清. 吉林省气候[M].北京:气象出版社,1997.
[21] 肖长来. 吉林西部地下水资源评价与水资源可持续开发利用研究[D].吉林大学,2001,6:6-13.
[22] 卢远. 吉林西部土地利用土地覆盖变化及其生态效应[D].吉林大学,2005,6:69-82.
[23] 赵凤琴. 吉林西部土地生态环境安全研究[D].吉林大学,2005,5:43-44.
[24] 吉林省地方志编纂委员会编纂[M].吉林省志.卷五.人口志.长春:吉林人民出版社,1992.
[25] 穆罕默德,赵莹莹,周道伟. 吉林省西部草原区地下水水质研究[J].生态学杂志,2002,21(2):14-16.
[26] 张勃夫. 吉林省地下水资源形成与分布的基本特点[J].吉林地质,1981,(2):24-30.
[27] 张博,洪梅,赵勇胜. GIS 支持下吉林省西部地下水特征组分分布特征研究[J].世界地质,2001,20(2):175-182.
[28] 刘春国,段长春,曹玉清. 含氟初始风化水化学组成的模型研究[J].长春地质学院学报,1995,25(4):423-427.
[29] 宋长春,邓伟. 吉林西部地下水特征及其与土壤盐渍化的关系[J].地理科学,2000,20(3):246-250.
[30] 曾昭华. 地下水中氟形成的控制因素及其分布规律[J].吉林地质,1997,16(4):26-30.
[31] 谷朝君,潘颖,潘明杰. 内梅罗指数法在地下水水质评价中的应用及存在问题[J].环境保护科学,2002,28(1):45-47.
[32] 吉林省地下水动态资料(白城、松原分册).吉林省水文水资源局.1998.
[33] 卞玉梅. 吉林省西部地下水动态研究[D].吉林大学,2005,6.
[34] 潘世兵,王忠静,邢卫国. 河流—含水层系统数值模拟方法探讨[J].水文,2002,22(4):19-21.
[35] 陈崇希,裴顺平,王逊. 非完整河的数值模拟方法及建模中的若干问题讨论[J].勘察科学技术,1999,(4):3-12.
[36] 边际,郭书会,司炳金. 数值模拟方法在评价地下水资源时区内河流的处理方法[J].工程勘察,1998,(5):32-35.
[37] 邵景力,崔亚莉,赵云章. 黄河下游影响带(河南段)三维地下水流数值模拟及应用[J].吉林大学学报(地球科学版),2003,33(1):51-55.
[38] 卢文喜. 地下水运动数值模拟过程中边界条件问题探讨[J].水利学报,2003,(3):33-36.
[39] 雅·贝尔. 地下水水力学[M].北京:地质出版社,1985:333-356.
[40] 孙纳正. 地下水流的数学模型和数值方法[M].北京:地质出版社,1981:42-56.
[41] 李同斌,邹立芝. 地下水动力学[M].长春:吉林大学出版社,1998:32-41.
[42] Wheeler M F. Roberson R R, Chilakapati A. Three-dimensional bioremediation modelling in heterogeneous porous media [J]. Water Resources. 1992,(2):299-315.
[43] 曹剑峰,冶雪艳,李升等. 河南境内黄河流域地下水系统划分与系统分析[J],吉林大学学报(地球科学版).2002,32(3):62-67.
[44] 薛禹群,谢春红,吴吉春. 水文地质数值法存在的问题及其对策[J].地球科学进展,1996,11(5):472-474.
[45] 陈崇希,裴顺平,王逊. 非完整河的数值模拟方法及建模中的若干问题讨论[J].勘察科学技术,1999,(4):3-12.
[46] 张力春. 吉林省西部水资源可持续利用的优化配置研究[D].吉林大学,2005,6.
[47] Bear, J. L.and Y. Bachmat.Introduction to modeling of transport phenomena in porous media [M]. Kluwer Academic Publishers,1991,(14):16～21.
[48] 吉林省水文地质图(1:500000)及其说明书[R].长春:吉林省地矿局水文队,1980.
[49] 吉林省地下水资源评价报告及附图(1:1000000)[R].长春:吉林省地质矿产局,1984.
[50] 吉林省地下水动态观测资料.长春:吉林省水文水资源局,1979~2000.
[51] 王大纯,张人权.水文地质学基础[M]. 北京:地质出版社,1995.
[52] 朱学愚,钱孝星,刘新仁. 地下水资源评价[M].南京:南京大学出版社,1987.
[53] 钱孝星. 水文地质计算[M].北京:水利水电出版社,1995.
[54] 周忠一. 区域水资源评价[M].北京:中国水利水电出版社,1996.
[55] 刘廷玺,朱仲元,马龙等. 通辽地区次降雨入渗补给系数的分析确定[J].内蒙古农业大学学报.2002,23(2):34-40.
[56] 陈洵洪. 用三维数学模型分析河流与含水层的相互作用,21 世纪中国水文地质学研究的新问题、新技术和新方法[M],北京:科学出版社,2001.
[57] Mark S. Nemeth1, Helena M.Solo-Gabriele Evaluation of the use of reach transmissivity to quantify exchange between groundwater and surface water [J]. Journal of Hydrology,2002,(5):201-208.
[58] 蒋放业,张新友. 河流与含水层耦合模型及其应用[J].地理学报,1999,54(6):48-55.
[59] 李宝元. 天堂河流区域浅层地下水运动数值模拟研究[D].中国农业大学,2002,3.
[60] 刘光亚. 水文地质数值模拟的拟合问题[J].水文地质工程地质,1992,19(4):2-4.
[61] 卢文喜. 地下水系统的模拟预测和优化管理[M].北京:科学出版社,1999.
[62] 左其亭,马军霞. 地下水系统中不确定性信息及其处理方法[J].水文地质工程地质.1994,21(5):41~43.
[63] 卢文喜. 地下水模拟预报工程中降水量的预报[J].勘察科学技术,1995,(4):8~10.
[64] 郭生练. 气温和降水序列的长期预测研究[J].水能源科学,1994,12(2):126~131.
[65] 卢文喜,邹立芝. 降水量的随机模拟[A],环境地学问题研究论文集[C].北京:石油工业出版社,1995:24-32.
[66] 卢晓燕,施鑫源,眭克仁等. 地下水资源系统时序管理模型[J].工程勘察,1999,(4):5-38.
[67] Hall, David;Fitzpatrick, Caroline S.B.Spectral analysis of pressure variations during combined air and water backwash of rapid gravity filters [J]. Water Research,1999, 33(17):3666-3672.
[68] 刘贯群,王秉忱,郑西来等. 内蒙古孪井地区降水量的时间序列分析预报模型[J].工程勘察,2003,(4):23-26.
[69] 杨位钦,顾岚. 时间序列分析与动态数据建模[M].北京:北京工业学院出版社,1986:122-130.
[70] Knotters M,PEV Van Walsum.Estimating fluctuation quantities from timeseries of water-table depths using models with a stochastic component [J]. Hydrol,1997,197(1):25-46.
[71] 叶青,亢锐等. 例说 Visual Basic 6.0 中文版程序设计[M].北京:北京科学技术出版社,2000.
[72] 徐振邦,娄元仁. 数学地质基础[M].北京:北京大学出版社,1994.:106-108.
[73] 卢文喜. 长春南湖水质成分的因子分析[J].世界地质,1999,18(3):61-65.
[74] 朱琳,苏小四. 吉林西部地区第四系潜水水质影响因素的 R 型因子分析[J].地球科学与环境学报,2006,28(1):51-56.
[75] 吕耀. 农业农业生态系统中造成的非点源污染[J].农业环境保护,1998,(1):36-40.
[76] 张思聪,吕贤弼,黄永刚. 灌溉施肥条件下氮素在土壤中迁移转化的研究[J].水利水电技术,1999,30(5):6-8.
[77] 赵勇胜,林学钰. 地下水污染模拟及污染的控制和处理[M].长春:吉林科学技术出版社,1994:50-53.
[78] 崔俊芳. 莱西地区地下水系统硝酸盐污染过程的研究[D].中国海洋大学,2003,6.
[79] 郑佳. 白城市地下水数值模拟与评价[D].吉林大学,2006,6:58-59.
[80] 郑西来、钱会、杨喜成. 地下水含水介质的弥散度测定[J].西安工程学院学报,1998,20(4):33~36.
[81] 成建梅. 考虑可信度的弥散度尺度效应分析[J].水利学报,2002,(2):90~94.
[82] 杨艳,王全九. 层状土溶质运移特征及其参数分析[J].灌溉排水学报,2005,24(6):19~21.
[83] www.people.com..cn 新农村.2007:中国肥价扑朔迷离.
[84] 王丽娟. 如何从环境学角度对地下水质量现状进行评价[J].西部探矿工程,2002,(5):73-74.
[85] 袁峰明. 土壤中的氮素淋失,氮素、产量、环境[M].北京:中国农业科技出版社,1996:196-216.
[86] 张玉珍. 九龙江上游五川流域农业非点源污染研究[D].厦门大学,2003,6:98-100.
[2] 魏林宏,束龙仓,郝振纯. 地下水流数值模拟的研究现状和发展趋势[J].重庆大学学报(自然科学版),2000,(23)(增刊):50-52.
[3] Pinder G. F., J.D. Bredehoeft. Application of the digital computer for aquifer evaluation . Water Resource[J],1968,4(5):1069-1093.
[4] 薛禹群,谢春红. 水文地质学的数值法[M].北京:煤炭工业出版社,1979.
[5] 丁继红,周德亮,马生忠. 国内外地下水模拟软件的发展现状与趋势[J].勘察科学技术,2002,(1):37-42.
[6] 周念清,朱蓉,朱学愚. MODFLOW 在宿迁市地下水资源评价中的应用[J].水文地质与工程地质,2000,(6):9-13.
[7] 林学钰. 地下水资源的基本属性和我国水文地质科学的发展[J].地学前缘,2002,23(3)(增刊):50～52.
[8] 吴吉春,薛禹群,黄海等. 山西柳林局部区域运移二维数值模拟[J].水利学报.2001,(8):38-43.
[9] 杨天行. 求解对流占优地下水水质模型的有限单元高精度广义迎风单元均衡格式及应用[J].长春地质学院学报,1996,(2):196-201.
[10] Freg Pohll etal.Modeling ground water flow and radioactive transport in a fratured aquifer. Ground Water [J], 1999,37(5):770-784.
[11] Henk Haitjema eatl.Selecting MODFLOW cell sizes for accurate flow field. Ground Water [J], 2001,39(6):931-938.
[12] Hank E.Julian etal.Numerial simulation of natural gradient tracer experiment for the natural attenuation study flow and physical transport. Ground Water [J], 2001,39(4):534-545.
[13] Jonathan Levy and Eric E.Ludy.Uncertainty quantification for delineation of wellhead protection areas using the gauss hermite quadrature approach. Ground Water[J] , 2000,38(1):63-75.
[14] A.M.Ebraheem·S.Riad.P.Wycisk.A.M.Sefelnasr. A local-scale groundwater flowmodel for groundwater resources management in Dakhla Oasis, SW Egypt. Hydrogeology journal [J], 2004,(12):714-722.
[15] 蒋辉. 河南沁阳市地下水水质的数值模拟[J].水文地质工程地质,1996,(2):14-20.
[16] 段青梅. 西辽河平原三维地质地质建模及地下水数值模拟研究[D].中国地质大学,2006.5:14-15.
[17] 周志芳,李艳. 复杂岩体地下水运动问题的有限分析法[J].水科学进展,1997,(3):26-42.
[18] 岳梅. 淮南市浅层地下水水质污染预测模型探讨[J].地下水,2002,(2):77-79.
[19] 吴剑锋,朱学愚. 由 MODFLOW 浅谈地下水流数值模拟软件的发展趋势[J].工程勘察,2000,(2):12-15.
[20] 王惠清. 吉林省气候[M].北京:气象出版社,1997.
[21] 肖长来. 吉林西部地下水资源评价与水资源可持续开发利用研究[D].吉林大学,2001,6:6-13.
[22] 卢远. 吉林西部土地利用土地覆盖变化及其生态效应[D].吉林大学,2005,6:69-82.
[23] 赵凤琴. 吉林西部土地生态环境安全研究[D].吉林大学,2005,5:43-44.
[24] 吉林省地方志编纂委员会编纂[M].吉林省志.卷五.人口志.长春:吉林人民出版社,1992.
[25] 穆罕默德,赵莹莹,周道伟. 吉林省西部草原区地下水水质研究[J].生态学杂志,2002,21(2):14-16.
[26] 张勃夫. 吉林省地下水资源形成与分布的基本特点[J].吉林地质,1981,(2):24-30.
[27] 张博,洪梅,赵勇胜. GIS 支持下吉林省西部地下水特征组分分布特征研究[J].世界地质,2001,20(2):175-182.
[28] 刘春国,段长春,曹玉清. 含氟初始风化水化学组成的模型研究[J].长春地质学院学报,1995,25(4):423-427.
[29] 宋长春,邓伟. 吉林西部地下水特征及其与土壤盐渍化的关系[J].地理科学,2000,20(3):246-250.
[30] 曾昭华. 地下水中氟形成的控制因素及其分布规律[J].吉林地质,1997,16(4):26-30.
[31] 谷朝君,潘颖,潘明杰. 内梅罗指数法在地下水水质评价中的应用及存在问题[J].环境保护科学,2002,28(1):45-47.
[32] 吉林省地下水动态资料(白城、松原分册).吉林省水文水资源局.1998.
[33] 卞玉梅. 吉林省西部地下水动态研究[D].吉林大学,2005,6.
[34] 潘世兵,王忠静,邢卫国. 河流—含水层系统数值模拟方法探讨[J].水文,2002,22(4):19-21.
[35] 陈崇希,裴顺平,王逊. 非完整河的数值模拟方法及建模中的若干问题讨论[J].勘察科学技术,1999,(4):3-12.
[36] 边际,郭书会,司炳金. 数值模拟方法在评价地下水资源时区内河流的处理方法[J].工程勘察,1998,(5):32-35.
[37] 邵景力,崔亚莉,赵云章. 黄河下游影响带(河南段)三维地下水流数值模拟及应用[J].吉林大学学报(地球科学版),2003,33(1):51-55.
[38] 卢文喜. 地下水运动数值模拟过程中边界条件问题探讨[J].水利学报,2003,(3):33-36.
[39] 雅·贝尔. 地下水水力学[M].北京:地质出版社,1985:333-356.
[40] 孙纳正. 地下水流的数学模型和数值方法[M].北京:地质出版社,1981:42-56.
[41] 李同斌,邹立芝. 地下水动力学[M].长春:吉林大学出版社,1998:32-41.
[42] Wheeler M F. Roberson R R, Chilakapati A. Three-dimensional bioremediation modelling in heterogeneous porous media [J]. Water Resources. 1992,(2):299-315.
[43] 曹剑峰,冶雪艳,李升等. 河南境内黄河流域地下水系统划分与系统分析[J],吉林大学学报(地球科学版).2002,32(3):62-67.
[44] 薛禹群,谢春红,吴吉春. 水文地质数值法存在的问题及其对策[J].地球科学进展,1996,11(5):472-474.
[45] 陈崇希,裴顺平,王逊. 非完整河的数值模拟方法及建模中的若干问题讨论[J].勘察科学技术,1999,(4):3-12.
[46] 张力春. 吉林省西部水资源可持续利用的优化配置研究[D].吉林大学,2005,6.
[47] Bear, J. L.and Y. Bachmat.Introduction to modeling of transport phenomena in porous media [M]. Kluwer Academic Publishers,1991,(14):16～21.
[48] 吉林省水文地质图(1:500000)及其说明书[R].长春:吉林省地矿局水文队,1980.
[49] 吉林省地下水资源评价报告及附图(1:1000000)[R].长春:吉林省地质矿产局,1984.
[50] 吉林省地下水动态观测资料.长春:吉林省水文水资源局,1979~2000.
[51] 王大纯,张人权.水文地质学基础[M]. 北京:地质出版社,1995.
[52] 朱学愚,钱孝星,刘新仁. 地下水资源评价[M].南京:南京大学出版社,1987.
[53] 钱孝星. 水文地质计算[M].北京:水利水电出版社,1995.
[54] 周忠一. 区域水资源评价[M].北京:中国水利水电出版社,1996.
[55] 刘廷玺,朱仲元,马龙等. 通辽地区次降雨入渗补给系数的分析确定[J].内蒙古农业大学学报.2002,23(2):34-40.
[56] 陈洵洪. 用三维数学模型分析河流与含水层的相互作用,21 世纪中国水文地质学研究的新问题、新技术和新方法[M],北京:科学出版社,2001.
[57] Mark S. Nemeth1, Helena M.Solo-Gabriele Evaluation of the use of reach transmissivity to quantify exchange between groundwater and surface water [J]. Journal of Hydrology,2002,(5):201-208.
[58] 蒋放业,张新友. 河流与含水层耦合模型及其应用[J].地理学报,1999,54(6):48-55.
[59] 李宝元. 天堂河流区域浅层地下水运动数值模拟研究[D].中国农业大学,2002,3.
[60] 刘光亚. 水文地质数值模拟的拟合问题[J].水文地质工程地质,1992,19(4):2-4.
[61] 卢文喜. 地下水系统的模拟预测和优化管理[M].北京:科学出版社,1999.
[62] 左其亭,马军霞. 地下水系统中不确定性信息及其处理方法[J].水文地质工程地质.1994,21(5):41~43.
[63] 卢文喜. 地下水模拟预报工程中降水量的预报[J].勘察科学技术,1995,(4):8~10.
[64] 郭生练. 气温和降水序列的长期预测研究[J].水能源科学,1994,12(2):126~131.
[65] 卢文喜,邹立芝. 降水量的随机模拟[A],环境地学问题研究论文集[C].北京:石油工业出版社,1995:24-32.
[66] 卢晓燕,施鑫源,眭克仁等. 地下水资源系统时序管理模型[J].工程勘察,1999,(4):5-38.
[67] Hall, David;Fitzpatrick, Caroline S.B.Spectral analysis of pressure variations during combined air and water backwash of rapid gravity filters [J]. Water Research,1999, 33(17):3666-3672.
[68] 刘贯群,王秉忱,郑西来等. 内蒙古孪井地区降水量的时间序列分析预报模型[J].工程勘察,2003,(4):23-26.
[69] 杨位钦,顾岚. 时间序列分析与动态数据建模[M].北京:北京工业学院出版社,1986:122-130.
[70] Knotters M,PEV Van Walsum.Estimating fluctuation quantities from timeseries of water-table depths using models with a stochastic component [J]. Hydrol,1997,197(1):25-46.
[71] 叶青,亢锐等. 例说 Visual Basic 6.0 中文版程序设计[M].北京:北京科学技术出版社,2000.
[72] 徐振邦,娄元仁. 数学地质基础[M].北京:北京大学出版社,1994.:106-108.
[73] 卢文喜. 长春南湖水质成分的因子分析[J].世界地质,1999,18(3):61-65.
[74] 朱琳,苏小四. 吉林西部地区第四系潜水水质影响因素的 R 型因子分析[J].地球科学与环境学报,2006,28(1):51-56.
[75] 吕耀. 农业农业生态系统中造成的非点源污染[J].农业环境保护,1998,(1):36-40.
[76] 张思聪,吕贤弼,黄永刚. 灌溉施肥条件下氮素在土壤中迁移转化的研究[J].水利水电技术,1999,30(5):6-8.
[77] 赵勇胜,林学钰. 地下水污染模拟及污染的控制和处理[M].长春:吉林科学技术出版社,1994:50-53.
[78] 崔俊芳. 莱西地区地下水系统硝酸盐污染过程的研究[D].中国海洋大学,2003,6.
[79] 郑佳. 白城市地下水数值模拟与评价[D].吉林大学,2006,6:58-59.
[80] 郑西来、钱会、杨喜成. 地下水含水介质的弥散度测定[J].西安工程学院学报,1998,20(4):33~36.
[81] 成建梅. 考虑可信度的弥散度尺度效应分析[J].水利学报,2002,(2):90~94.
[82] 杨艳,王全九. 层状土溶质运移特征及其参数分析[J].灌溉排水学报,2005,24(6):19~21.
[83] www.people.com..cn 新农村.2007:中国肥价扑朔迷离.
[84] 王丽娟. 如何从环境学角度对地下水质量现状进行评价[J].西部探矿工程,2002,(5):73-74.
[85] 袁峰明. 土壤中的氮素淋失,氮素、产量、环境[M].北京:中国农业科技出版社,1996:196-216.
[86] 张玉珍. 九龙江上游五川流域农业非点源污染研究[D].厦门大学,2003,6:98-100.