专题四 FEFOW应用实例.doc

专题四 FEFLOW应用实例 第五章 FEFLOW在宿迁市的应用 5.1 水文地质条件概化 本研究区周长为758km,面积为7489km2。模拟区边界除西北部骆马湖和东南部洪泽湖概化为一类水头边界(1st head)外，其余均概化为三类水头边界(3rd head)。补给主要为大气降雨入渗。排泄为人工开采地下水和垂向蒸发。从地下水动态资料来看，模拟区的孔隙水水位呈非稳定状态。模拟区内部结构视为非均质各向同性，根据不同区域的水文地质条件对不同参数进行相应的参数分区，同一参数分区内可视为均质。水流服从达西定律，可视为平面二维流。 5.2 数学模型的建立 根据上述水文地质概念模型，可建立起相应的数学模型： ， ， , t0 ， 式中： H ：承压含水层水位(m)； h ：潜水含水层水位(m)； ：弹性储水（释水）系数； ：越流系数； T：导水系数(m2/d)； QV：承压水开采强度(m/d)； H0：初始水位(m)； H1：计算区一类边界水位(m)； Γ1：一类水头边界。 D：计算区范围。 5.3 含水层系统的识别 地下水流系统的模拟，要解决两大问题，即参数识别问题和预报问题。解决这两个问题的前提是要建立一个尽可能反映实际含水层地质体与边界条件的数学模型。本次评价中，因无最新资料可用来进行模型识别，所以直接参考1999年《江苏省宿迁市地下水资源评价报告》，仅利用模型进行预报。 5.3.1 空间离散 根据研究区含水层系统的垂向分布特征，将研究区含水层系统分为三层（layer）四片（slice），上层为潜水和微承压含水层，下层为Ⅱ、Ⅲ承压水含水层，中间为弱透水层。平面上将其剖分成32991个有限元网格，共23680个节点（图5—1）。 图5—1 平面网格剖分图 5.3.2 时间离散 FEFLOW中的时间步长形式主要有三种（见图5—2）：定步长（constant time steps）、变步长（varying time steps）、自动时间步长（automatic time steps）。为了数据的输入和输出的方便，此处选择变步长。预报期选定在2007年1月1日到2020年12月31日，共5114天。采用变步长，按照月为单位计，共168个时间步长。 图5—2 时间步长控制图 5.3.3 初始流场 初始水位根据2000年第一季度《宿迁市深层地下水监测季报》（表5—1、5—2）而确定。研究区深层的初始流场如图5—3所示。 表5—1 宿迁市2000年1—3月份深层地下水监测井水位埋深统计表 县（区） 井号 地理位置 1月1日埋深（m） 宿城区 1302002 市供电局 5.94 泗 洪 县 1318001 龙集水利站 27.37 1318003 天岗湖乡政府 12.69 1318004 金锁镇自来水厂 4.50 1318005 陈圩乡水利站 10.00 1318006 魏营乡王拐村王元庄 24.82 1312008 成河自来水厂 1.82 1312009 孙园乡井村 16.06 1311010 双沟中心卫生院 11.98 1312011 峰山乡黄岗村翟岗 33.78 1312012 重岗乡水利站 8.97 1313013 界集自来水二厂 20.06 1313014 大楼乡 11.01 1318016 县人民银行 9.11 泗 阳 县 1328002 王集自来水厂 7.92 1328003 新源自来水厂 19.53 1328004 八集粮管所 7.33 宿豫区 1338001 侍岭缫丝厂（3号井） 16.76 沭 阳 县 1342001 高墟乡政府 2.73 1342002 吴集乡柳树村 4.49 1342004 吴集乡柳树村 2.80 1348005 悦来乡双蔡村 16.74 1348006 耿圩乡沂南村 10.02 1348007 沂涛乡敬老院 3.24 1348009 胡集乡盆荡村 6.06 1348010 扎下乡 11.44 1342011 刘集乡吉湖村 10.08 表5—2 宿迁市2000年1—3月份逐日深层地下水监测井水位埋深统计表 县（区） 井号 地理位置 1月1日埋深（m） 沭阳县 1348012 华冲镇淀粉厂 1.63 1348013 庙头村乡政府 2.44 宿城区 1301001 井头乡彩塑厂 2.7 泗 洪 县 1318002 太平乡粮管所 14.89 1311007 归仁镇政府 11.77 1313015 半城镇大新庄 14.91 泗 阳 县 1321001 县绢纺厂 14.14