多含水层系统的程序设计及其意义.pdfVIP

第29卷第4期 中国矿业大学学报 Vcl．29N。．4 of 2000年7月 China of Journal UniversityMining＆TechnologyJul．2000 文章■号，1000-1964(2000)04—0419—04 多含水层系统的程序设计及其意义 许光泉1，李慈君2，张铭。 (1．淮南工业学院资源与环境工程系，安徽淮中国地质大学，北京100083) 摘要：针对我国平原地区含水层系统的特点，建立了多层含水层系统数学模型，给出了程序设计 的总体思路．对于多节点形成的高度稀疏方程组，运用共轭梯度或sSoR一共轭梯度对此进行 求解，程序的设计对于多含水层系统的水资源评价具有理论和现实意义． 关键词：多舍水层系统}模型；程序设计 中图分类号：P64 文献标识码；A 我国许多冲积型山前平原如华北平原、松嫩平 原其地下水一般为多层含水层．随着对该类型平原 无关，弱隔水层可有储水释放； 地区地下水的进一步开采，深层地下水开采资料及 深层观测孔的长期系列资料越来越丰富，使得准确 小时段； 计算多层含水层系统的地下水模拟模型成为可能． 6)不能混合抽水和混合观测． 两层以上的含水层系统称为多层的含水层系统，它 1．2模型建立 们之间存在越流、‘天窗’的联系，系统边界存在补 根据以上模型要求可用如下微分方程的定解 排关系．含水层系统的第一层多为潜水层，它存在 问题来描述： ． 大气降雨、地表水体(河流、湖泊、沼泽、水库及回灌 入渗)的入渗和潜水通过饱气带的蒸发、蒸腾作用 的交换，从理论上说，这种系统可分为M层．在人 为因素影响下，对各含水层之间进行开采时存在着 相互干扰的情况．从含水层形成的古沉积环境来 cM(^㈣一¨一口。=晶警 说，不外分以下3种情况[1]：1)以河流相沉积为主 (m=2，3，…，肘一1)； (2) 的多含水层系统，其含水层富水性强、储存量较大， 第^f层div(TM·gradh_iIf)+ 地下水资源丰富}2)以湖相沉积为主的多含水层 系统，其富水性差，储存量小，地下水资源贫乏；3) 介于以上二者之间河湖交替相的多含水层系统；本 文以此为地质概念模型进行程序设计． 1数学模型 (m=1，2，…，M)， (5) 1．1建模要求 对松散的多含水层系统来说，模型要求如下： 1)具有相互联系的多含水层系统其水流近似 水平运动，通过弱隔水层时水流是垂直的； 2)含水层的水位降深与含水层厚度之比很 小； 3)含水层为非均质、各向异性； 收藕日期l1999一ll一02 万方数据 作誊倚介。许光泉(1967一)，男．安徽省稚南市人．淮南工业学院讲师，工学硬士，从事水资薄与环境研完 中国矿业大学学报 第29卷 生成时段文件输出，在其过程中对各分区(按行政 水强度，(L／T)，q。=∑‰(奎，f)d(；～赫。)，其中 酉 分区或地质单元分区)的大气降雨人渗量、灌溉入 ％o，f)为第m含水层在；点的抽水强度；z。，为第渗、河渠人渗、及农业开采及各层的工业开采经处 m含水层的第i个抽水点的坐标；Ⅳ。为第m含水层 理后，生成一个后处理文件，为