松原壅水坝工程库区左岸浸没预测三维数值模拟.pdfVIP

第 28卷 第 3期 农 业 工 程 学 报 、b1．28 No-3 2012芷 2月 TransactionsoftheCSAE Feb．2012 129 松原壅水坝工程库区左岸浸没预测三维数值模拟 骆祖江 ，杨 林 ，李 占军2，王俊杰2 (1．河海大学地球科学与工程学院，南京 2100981 2 中水东北勘测设计研究有限责任公司，长春 130021) 摘 要：针对平原区水库蓄水后引起的浸没问题，以松原壅水坝工程库区左岸浸没为例，在概化出库区左岸水文地质概 念模型的基础上，建立了松原壅水坝工程库区左岸浸没预测的地下水三维非稳定流数值模拟模型，并采用有限差分模拟 方法，按无防渗墙和长为4．05、8．05及 12．05km 3种不同防渗墙蓄水工况方案进行预测，获得了无防渗墙蓄水工况4— 10月底逐月的和3种不同防渗墙蓄水工况8月底的浸没范围。计算结果表明：无防渗墙蓄水工况方案4月份浸没范围最 小，总面积为 127．62x10m。，随后逐渐增大；8月份浸没范围达到最大，总面积为 1473．94x10m ，而后又逐渐减小； 随着在潜水含水层沿库区左岸布置防渗墙的加长，浸没范围逐渐减小，但减幅不大。该研究可为今后水库蓄水引起的库 区浸没模拟预测提供参考。 关键词：数值模拟，三维，非稳定流，防渗，浸没，预测 doi：10．3969／j．issn．1002—6819．2012．03．023 中图分类号：P641 文献标志码：A 文章编号：1002—6819(2012)一03—0129—06 骆祖江，杨 林，李 占军，等．松原壅水坝工程库区左岸浸没预测三维数值模拟 [J]．农业工程学报，2012，28(3)：129 — 134． LuoZujiang，YangLin，LiZhanjun，eta1．Three—dimensionalnumericalsimulationforimmersionpredictionofleftbankof Songyuanreservoirproject[J]．TransactionsoftheCSAE，2012，28(3)：129—134．(inChinesewihtEnglishabstract) 地下水三维非稳定流数值模拟模型，并结合工程规划， 0 引 言 按无防渗墙和 3种不同防渗墙蓄水工况方案进行模拟预 水库蓄水后引起的壅水浸没是常见的环境地质问题fl。 测，计算了各方案的浸没范围，并通过对比分析，确定 以往对库区浸没的计算中，常采用解析解、地下水二维 了本次工程给浸没带来的影响，使计算结果更加真实、 流或准三维流模型I3】进行研究，且将地下水流态概化为稳 可信，为工程设计提供了可靠的科学依据。 定流，考虑的因素简单，在大多数情况下不能较好地刻 1 工程概况 画库区整个水文地质体及实际工况特征。如地下水二维 流模型切断了地下水系统各含水层之 间的垂向水力联 松原壅水坝位于第二松花江下游松原市城区段，前 系，根本无法模拟计算复杂的由许多不连续砂体、淤泥、 扶大桥下游约4．5km处，上距哈达山水利枢纽工程坝址 粘土组成的地下水系统。准三维模型假定地下水在含水 约32km。壅水建筑物正常蓄水位 l31．501TI，相应的回水 层中为平面二维流，在弱含水层中为垂向一维流，这种 长度约为 14km，库容为0．38~10rn，工程等别为III等， 方法虽然较二维模型前进了一步，但忽略了弱含水层 自 工程规模为中型，主要建筑物级别为 3级。水库区两岸 身的释水机理，简化了各含水层之间的水文地质机制， 分属不同的地貌单元，右岸为上更新统冲积平原，地面 也不能整体模拟地下水系统的三维结构和渗流特征，常 高程 145～155m，多为陡岸，高出江水面约 10～13