岩溶管道水系统物理模拟_以岩滩水电站板文地下河系为例

《岩溶管道水系统物理模拟_以岩滩水电站板文地下河系为例》由会员分享，可在线阅读，更多相关《岩溶管道水系统物理模拟_以岩滩水电站板文地下河系为例（7页珍藏版）》请在金锄头文库上搜索。

1、第15卷 第4期 中 国 岩 溶 Vo1. 15 No. 41996年12月 CARSOLOGICA SI N ICADec. 1996岩溶管道水系统物理模拟以岩滩水电站板文地下河系为例 李文兴 郭纯青(地质矿产部岩溶地质研究所)摘 要 从物理模拟的角度出发,论述应用相似原理概化岩溶管道水 系统和建立符合能量转化、 质量守恒的等效水箱,并采用变径管束组合 等效岩溶管道。还介绍了改变模型结构和不改变模型结构的两种模型 预报。为解决岩溶地区岩溶管道系统水流运动问题提供很好的方法。关键词 岩溶 管道水系统 物理模拟0 引 言为研究岩溶管道水运动而进行的物理模拟,从模拟内容上可以分成系统模拟与功能模拟

2、。 本文以岩滩水电站库区板文地下河系为例讨论系统模拟。 从能量转换、 质量守恒的基本物理法则出发,对系统概化出具有介质特征和水流特征的实体模型,根据相似理论对实体模型进行缩 小,组装成能在室内模拟的物理模型。物理模拟就是要组装并测量模型各物理量的过程。实体 模型从结构上可分为储水介质与导水介质。而物理模型从结构上可分为水能存储单元与水能 输送单元。水能存储单元与储水介质相对应;水能输送单元与导水介质相对应。本文举例的水 能存储单元包括地上、 地下储水库容。1 物理模拟模型设计以相似理论为基础,在模拟的原型(实体模型)与模型(物理模型)之间建立一种相似关系。 具体的做法是将原型的物理量换算成模型

3、的物理量,对模型的物理量进行模拟试验;然后将模 拟试验的结果转换到原型中去,就可以得到需要求解的物理量。 具体要求为:几何相似、 运动相 似、 动力相似、 边界条件相似。一般情况下,水流受到一定外部边界的约束,边界的变化使水流 运动状态也发生变化。 从水力学的观点上看,上述过程也就是作用于水的各种力共同作用的结果。要使原型与模型两者水流间形成相似关系,除了边界条件相似以外,更重要的是考虑到受第一作者简介:李文兴,男, 1957年出生,高级工程师, 1982年毕业于中山大学数学力学系力学专业。(541004)广西桂林市七星40号。 收稿日期: 1996- 01- 16;改回日期: 1996- 0

4、9- 04。力条件的相似。 根据相似原理的重力相似准则:如果物理模型与实体模型两个水流在重力作用 下相似,则它们的佛汝德数相等。即v2g l= 1(1)式中v为水流流速,g为重力加速度,l为长度。l为长度比:l原型?l模型;v为原型与模型水流速度之比:v原型?v模型;g为原型与模型重力加速度之比:g原型?g模型。因为原型与模型的液体都是水,因此有g原型=g模型,即:g= 1。所以公式(1)成为:v=gl015(2)从而,流量的比例尺 Q为:Q=Q原型 Q模型=v原型A原型 v模型A模型=vl2=l215(3)式中Q原型、Q模型和A原型、A模型分别为原型与模型的流量和管道及裂隙过水断面面积。2

5、岩滩电站库区板文地下河系物理模型本物理模拟模型的水能输送单元为岩滩库区板文地下河。板文地下河位于广西东兰县南部、 红水河岩滩水电站上游西侧,分布范围南北长约24km ,东西宽约2125km。该地下河系是桂西地区一条大型地下河,流域面积352km2,干流全长2613km。地下河出口高程为177m。 地下河上游源于江平一带峰丛洼地,中游(即淹没区)为峰丛(林)谷地,下游又为高峰洼地。当 中、 上游来水集中于中游谷地时,主要靠纳亮、 拉垌、 拉平三大消水洞群(共130个)消水至地下 河道,经过14km的下游峰丛区,由板文出口排入红水河(见图1)。 岩滩水库蓄水前板文地下河图1 板文地下河系水箱概化示

6、意图Fig. 1 Block diagram illustrating the tank structure in the Banwen subterranean river system1.流域边界; 2.地表河; 3.内涝区汇水界线; 4.地下河; 5.雨量站; 6.水位站253 中国岩溶 1996年出口的红水河水位为177m ,水库蓄水后水位抬升到222m ,使得本来就排水不畅的巴纳、 拉平 两谷地的水流受到的阻力更大,淹没的时间更长。本模型以巴纳、 拉平两谷地及地下岩溶含水 介质为水能储存单元,岩滩水库作为水能输送单元下游的蓄能单元。 两谷地到板文地下河出口 为水能输送单元。其中,水能

7、输送单元用大管道(胶管或玻璃管)加阻力元件来模拟,以阻力元 件等效岩溶管道阻力。 水能储存单元以物理水箱来实现,水箱水位模拟谷地水位及岩滩水库水位。为了模拟的方便,模型的起点高度(模型零点)定在实际标高的某一高度(本模型的起点标 高为177m)。 根据岩滩库区板文地下河的地质及水文地质条件,取模型的长度比为l= 100,应用(2)、(3)式可求得 v= 10,Q= 100000,具体见表1。表1 模型数据表Tab. 1 M odel parameters名 称原型模型 名 称原型模型长度比100m1m基点(建库前红水河水位)177100m0100cm面积比10000m21m2岩滩水库水位222

8、100m45100cm体积比1000000m31m3模拟巴纳93年初始水位301199m124199cm速度比10m?s1m?s模拟拉平93年初始水位289199m112199cm流量比100000m3?s1m3?s模拟巴纳93年最高水位307138m130138cm巴纳最低点297m120cm模拟拉平93年最高水位300184m123184cm拉平最低点290m113cm模拟巴纳93年最终水位302158m125158cm模拟拉平93年最终水位296125m119125cm图2 物理模拟试验过程图Fig. 2 Flow chart illustrating physicalsi mulati

9、on test3 物理模拟试验过程物理模拟试验过程如图2所示。3. 1 水箱水位及水量调试 水箱水位及水量的调试是根据地质、 地貌、 地形的特征概化出在不同水位标高下的水面面积及水体体积,这就是可控水 箱的基本数据,将其输入计算机,作为模拟 试验中物理水箱的输入水量。拉平和巴纳 两谷地在不同水位标高下的水量见表2及 表3。3. 2 管道流量及阻力调试 从图2可知,模拟试验的过程,就是根 据已知条件,主要求解管道的阻力参数的 过程。 而管道流量及阻力的模拟又是一个极其复杂的过程,要求的技术很高,由于篇幅有限,这 里只给出本模型的试验结果。根据能量方程有:Q(t) =v(t) =(h(t) -h0

10、(t)1?2(4)353第15卷 第4期 李文兴等:岩溶管道水系统物理模拟 表2 巴纳片不同水位高程的蓄水面积及库容数据表Tab. 2 Storage area and storage capacity at different water level elevation,Bana水位(m)面积(万m2)库容(万m3)水位(m)面积(万m2)库容(万m3)水位(m)面积(万m2)库容(万m3)297055103022904853077703415298206210303450861308820418529940821030458013903098605020300751311530567020

11、253108905900301180255103067302700表3 拉平片不同水位高程的蓄水面积及库容数据表Tab. 3 Storage area and storage capacity at different water level elevation,L aping水位(m)面积(万m2)库容(万m3)水位(m)面积(万m2)库容(万m3)水位(m)面积(万m2)库容(万m3)29008102971604681030426620051029120251029818062510305267227010292365410299200800103062862551102935495103

12、0021810081030729428411029475150103012301240103083033139152951002261530224514901030931034461029613533510303255174510310318376010式中:Q(t)为管道水流流量(m3?s);为管道过水断面面积(m2);v(t)为管道水流流速(m?s);为管道流量参数,即模型参数(m5?2?s);h(t)为谷地水位标高(m);h0(t)为地下河出口水位标高(m)。对模型进行试验后,得到巴纳、 拉平谷地的流量参数如下:= 11581(5)=01577 当h(t)2950135355 当h(t)

13、 295(6)参数 为岩溶管道系统的综合流量参数(在这里为了使叙述方便,将管道阻力参数的倒数叫流 量参数),反映了两谷地在不同的水位标高下,地下河对水流的作用。从模拟结果看,巴纳谷地 比拉平谷地的流量参数大,也即其阻力参数小,说明巴纳地下河比拉平地下河畅通。而从拉平 谷地的地下河流量参数看,在不同的水位标高下也是不同的。当水位高于等于295m时,流量 参数小,反之则流量参数大。它还反映了管道系统的分叉情况。3. 3 模型水位及流量的测试 在模拟试验中,物理水箱水位的测量采用测压管和压力传感器加微机监测采集来实现;管 道流量的测量采用转子流量计和压差传感器加微机监测采集来完成。其中传感器加微机监

14、测 采集是用来对模拟试验的结果进行收集并存盘保存,而测压管和转子流量计则用来校正取数 的准确性。另外,流量计还起到管道阻力元件的作用。在水箱不同的水位高程上,蓄水面积及体积的变化,用计算机通过电磁阀来控制其水箱进水量来实现。 微机自动控制采集如框图3所 示。453 中国岩溶 1996年图3 物理模型水位及流量自动监控示意图Fig. 3 Sketch show ing the waterlevel and discharge of physicalmodel on automatic control4 模型预报根据模拟的目的,模型预报可分为改变模型 结构的预报和不改变模型结构的预报。4. 1 改

15、变模型结构的预报 在水坝设计及矿区开采设计时,往往需要在现在模型的基础上,预测未来在系统结构发生变 化导致其流场的变化,或者在水坝建成及矿区开 采以后,要恢复原来的水流流场,都需要采用改变 模型结构的预报方法。对本物理模型采用改变模 型结构的模型预报,就是在目前的降雨状态下,假设岩滩电站不蓄水,即地下河出口的水位标高不 是目前水位,而是水库蓄水前的177m。在模型预 报时,就要改变模型的结构。 这样模拟出的两谷地 水位过程线见图4。图4 1993年巴纳片、 拉平片水位模拟曲线图Fig. 4 Si mulated water level fluctuation at Bana and L aping in 19934. 2 不改变模型结构的预报如果要了解模型在不同的降雨量及降雨强度下,将会产生的水位效应及流量效应,就要采用不改变模型结构的预报方法。具体做法为:保证在岩溶管道水系统水动力状态(即模型介质和参数)没有发生变化的情况下,根据当地历史降雨系列,求出不同降雨量发生的频率,设计出不同的降雨量,该降雨量即是模型的输入函数,对模型进行模拟,求解水位及流量等物理量。 对岩滩水电站库区板文地下河系,采用频率(P)分别为10%、20%、50% ,时段分别为1、3、7、15d的降雨量,作为模型的输入函数,对模型进行预报,结果见表4。553第15卷 第4期 李文兴等:岩溶管道水系统物理模拟