溶质运移理论水动力弥散系数计算方法.pptVIP

* 二、二维水动力弥散-瞬时投放示踪剂 4.弥散晕面积求参法 表明t时刻的一系列浓度与对应的椭圆面积在半对数坐标上为一直线，斜率为 求出横轴截距d，有（5-61）得 可求有效孔隙度 由（5-60）得 求出浓度C等值线围成的椭圆面积随时间变化值，则S/t与t在半对数坐标中呈直线。同理可求上述参数 联立，可得DL、DT (5-62) * 二、二维水动力弥散-瞬时投放示踪剂 算例 在示踪剂浓度溶液中氯化物浓度为580.7mg/L，体积共0.7m3，分5井瞬时注入。由于含水层非均质，示踪剂分为流速快、慢移动。慢带实际平均速度 m/s，121d后，如图所示，试求慢带参数 * 二、二维水动力弥散-瞬时投放示踪剂 计算步骤： 1.量取不同浓度的椭圆长、短轴，按比例换算成实际长度，求出相应面积及两轴比值。 2.舍去异常值，求出a/b平均值，此处为1.7614。 3.作lnC-S关系图，求出斜率，此处K=-18.453 * 二、二维水动力弥散-瞬时投放示踪剂 计算步骤： 得 则 则 弥散晕面积法优点： 1.不一定要求观测孔在主流线上（但要求较多观测孔）； 2.求弥散系数不需知道流速且可求孔隙度 3.数据在时空上均有一定延伸，代表性较好 * 二、二维水动力弥散-瞬时投放示踪剂 5.lnC-t曲线拐点法 对式 两边取自然对数并求时间的二阶导再取二阶导为0，得 ti表示lnC-t曲线上出现拐点时相应的时间。 (5-65) * 二、二维水动力弥散-瞬时投放示踪剂 5.lnC-t曲线拐点法 若（x’,y’）处另一观测孔中出现拐点的时间为ti’有 联立，可得DL、DT 优缺点： 1.不需求流速且不要求观测点在主流线上； 2.拐点的选择有一定任意性。 (5-66) * 二、二维水动力弥散-瞬时投放示踪剂 6.优选拟合法 将上述方法所计算出的参数视作初值，再据弥散试验中某孔的观测资料与（5-65）式的计算相拟合，从而通过计算机进行参数优选，采用目标函数 来反映计算浓度值与实测浓度值之间的拟合程度 * 二、二维水动力弥散-瞬时投放示踪剂 算例 在一组观测孔中所测得的电导率数据为基础，采用优选拟合法计算水动力弥散参数。因电导率与浓度近似线性，直接采用电导率ρ进行拟合计算。结果如下： * 二、二维水动力弥散-连续投放示踪剂 连续投放示踪剂 对 求t的一阶导，得到浓度随时间的变化率后，可借用瞬时投放的解法 1.lgWc（A, β）-lgA标准曲线法 记 对（4-70）（5-68’）两边取对数，得 （4-70） （ 5-68’ ） * 二、二维水动力弥散-连续投放示踪剂 1.lgWc（A, β）-lgA标准曲线法 (1)作对应不同β值的标准曲线族，做出实测曲线与标准曲线进行拟合； (2)找出拟合最好的β0 ； (3)找出与A=1，WC=1相对应的C0、t0 ； 实测曲线与理论曲线形状相同，只发生坐标平移 * 二、二维水动力弥散-连续投放示踪剂 将u、DT、DL视作待求参数，联立上述3式可得结果 * 二、二维水动力弥散-连续投放示踪剂 算例 取 * 二、二维水动力弥散-连续投放示踪剂 * 二、二维水动力弥散-连续投放示踪剂 2.lg[eBK0（B）]-lgB标准曲线法 当 时，有 观测孔位于x轴上，有 其中， 实测曲线与理论曲线形状相同 只发生坐标平移 * 二、二维水动力弥散-连续投放示踪剂 2.lg[eBK0（B）]-lgB标准曲线法 类似地，找出B=1,eBK0（B）=1的点x0、C0，可得 地下水溶质运移理论及模型 第五章 水动力弥散系数的计算方法 中国地质大学环境学院 2014春 * 第五章 水动力弥散系数的计算方法 假设： （1）地下水流速较大，可忽略分子扩散系数； （2）弥散系数与孔隙平均流速呈线性关系，可先求弥散系数再除以孔隙平均流速而获取弥散度或反之。 计算方法： （1）利用解析公式直接或间接求解； （2）采用标准曲线对比法 不讨论三维弥散的情况 * 一、一维水动力弥散 利用正态分布函数性质求参 （1）瞬时注入示踪剂 其解为 * 一、一维水动力弥散 DL求解步骤： 1.确定Cmax的值及其对应的数学期望， 2.取0.607 Cmax，从浓度曲线上找出拐点位