湖泊水库水温预测.pptVIP

湖泊水库水温预测 经验公式法 指数函数型垂向水温分布经验公式 幂函数型垂向水温分布经验公式 模型模拟法 垂向一维对流扩散模型 风力混合层模型 指数函数型垂向水温分布经验公式 水利水电工程水文计算规范中推荐的： 计算公式 幂函数型垂向水温分布经验公式 李怀恩1990年提出 幂函数型垂向水温分布经验公式 有关参数的确定 当有实测数据时：利用实际数据进行参数拟合，可采用最小二乘法。分月率定Tc,hc,A,B。一般来说B值比较稳定，全年可取固定值，一般取3 缺乏实测资料 1.直接利用有实测水温资料的相似水库的优选参数 2. 综合分析参数变化规律确定 B值比较固定，全年可取常数 Tc,i与上月气温Ta,i密切相关，可根据类似水库建立好的这种关系，由本水库的气温资料确定 水库表层温度、底部温度容易根据水库所处的纬度、气温和水深H确定。可以带入前式，推求hc和A 数学模型法 垂向一维对流扩散模型 模型的建立 温度在水平方向基本均匀一致，垂直方向则表现出很大的变化。 微层体积元的水热交换 数学模型法 模型的建立 对于微时段dt该微元的各项热交换为： 纵向移流的净热输送 垂向移流的净热输送 垂向分子扩散作用的净热输送 数学模型法 模型的建立 垂向太阳辐射的净热输送 长波辐射、蒸发和传导的热交换都只发生在水的表层，只有太阳辐射会进入到下面： 表层以下的净热输送 数学模型法 模型的建立：水温模型基本方程： 数学模型法 模型的建立 ?按下式计算 ， ?= ?I (1-β） exp[-η(zS - z)] 其中 zS ， z— 水面及水面下某一界面的高程，m； ?I —水面的太阳净辐射，等于总入射量I减去反射量RI， ?I = I – RI β—?I 在表面被吸收的分数，约0.4~0.5； η—消光系数，和透明度、SS含量及矿化度有关，一般在0.1~2.0m-1之间 数学模型法 模型建立 数学模型法 模型建立- 对垂直湍流混合作用的处理 本方程中未包含此项。因为只是在湖面温水区域，当温度分布不稳定而引起密度分布不稳定时，才发生对流混作引起的热迁移。 对实际计算来水，每当计算中出现这种不稳定时，即上层的密度大于下层密度时候，就要对一定深度的水温进行平均。 这种平均一直进行到由温度引起的不稳定密度分布消失为止的深度。 数学模型法 MIT水温模型的数值方法 温度模型的差分方程 差分体系：显示体系 水库分层与坐标, 底层为原点， 向上为正 数学模型 MIT水温模型的数值方法 ΔT在表层以下的计算公式如下： 数学模型 MIT数值方法 对表层水体还要计算水表面对太阳辐射的吸收项ΔT5和表面热损失项ΔT6，此时： MIT数值方法 模型的运算流程 1）计算大气热交换 由时断末的气象资料计算j+1时断的φI = I – RI、 φG=G-RG、S、 φE、 φC。 2）各层入、出流计算 由j时段的水文资料和初始、边界条件 经调节计算，求得j+1时段的库水位zsj+1、分层数n 计算入流水温及流速高程分布TI，ij+1、 uI，ij+1 出流的流速沿高程分布uO，ij+1； 自下而上分层作水量平衡计算，求得各层界面处的垂向流速vij+1 MIT数值方法 模型的运算流程 3）计算各层水温 按上述ΔT的差分计算公式，由（1）、（2）结果计算各层的ΔTij+1，进而求得Tij+1。 4）考虑热力不稳定作用对水温分布进行调整 当上步计算的水温分布使上层水的密度大于下层时，即发生对流混合。4度以上时，密度和温度负相关。所以上步算完后，立即自上而下检查水温分布，如表层水温Tj+1s低于下层的Tj+1s-1就使这2层混合，如果混合后的温度又比s-2层的低，就再混合，直到混合层的温度不比下面那一层低为止。照此调整，直到整个水温分布不存在不稳定结构为止，此时得到的才是最后要求的（j+1）?t时刻的水温分布。 MIT数值方法 MIT模拟过程 5）计算下一时段的水温分布 以第（4）步计算的水温作为该时段的初始条件，及可采用（1）~（4）步的方法推求该时段的水温分布。 模型应用中几个问题的处理 计算步长问题 本模型是一个显式差分体系，计算步长的选取必须保证计算的稳定性。稳定条件为： 模型中几个问题的处理 水表面高程及表层厚度的确定 时段末的水位根据时段内的入、出流量和起始水位，通过水量平衡求得。表层的厚度Δz受水位涨落变化控制，不是一个固定值。通过下式计算 模型中几个问题的处理 水库的入流水温及入、出流垂向分布问题 模型中几个问题的处理 水库的入流水温及入、出流垂向分布问题 用类似方法假定出流流速为在泄