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

11 一、一维水动力弥散-浓度曲线线性化求参 则绘出其曲线图后，可根据斜率和截距求出弥散系数 运用全部试验资料 12 二、二维水动力弥散-瞬时投放示踪剂 1.逐点求参法 3个时刻t1、t2、t3对应浓度C1、C2、C3 ，由式 当y=0时，有 联立， 瞬时投放示踪剂 （5-16） （4-65） 13 二、二维水动力弥散-瞬时投放示踪剂 则 若只考虑（5-16）式 令y=0，由（4-65）导出横向弥散系数 理论上，每3组浓度-时间数据可求出u、DL、DT 14 二、二维水动力弥散-瞬时投放示踪剂 实际中，仅用3组数组求出的参数具有一定随机性甚至失去物理意义，故需要一些列C-t数据。可先去掉不符合物理意义的数值，再将其余参数算术平均值即可视作待求参数的近似值。如下表： 15 二、二维水动力弥散-瞬时投放示踪剂 2.直线图解法 （5-23） 16 一、二维水动力弥散-瞬时投放示踪剂 若观测孔位于x轴上，（5-23）可简化成 有 可通过求u来求DL 改写式子 令 则 （5-27） （5-30） 综合（5-27）（5-30） 17 二、二维水动力弥散-瞬时投放示踪剂 通过证明可得出下式（过程略）： （5-32） 取 （5-32）写成 参数计算的具体步骤 1.从实测的C-t数据序列中找出Cm、tm值，当观测数据较少时，可先作出C-t曲线后，从曲线上查出Cm、tm值 2.计算两组X、Y，绘在直角坐标系中，两轴比例一致，再量取R 18 二、二维水动力弥散-瞬时投放示踪剂 3.代入上述公式求DL、u 4.利用不在x轴上的观测孔资料，求出其tm’，根据 (5-23)求出DT 19 二、二维水动力弥散-瞬时投放示踪剂 算例 在某次瞬时投放示踪剂室内弥散试验中，测得一距投放井80cm且位于x轴上的观测点处的电导率和时间数据（见课本P62），相应的实测曲线见图5-5 20 二、二维水动力弥散-瞬时投放示踪剂 在表5-2中得出特征数据如下： 根据式 求出X、Y ρ-t曲线的上升段与下降段分别对应X-Y曲线的两条直线段。理想情况下，两直线重合，实际中其交角的大小反映了理论模型与实际模型的偏离程度，分别量取斜率R1和R2，提出DL和u后平均，再求取DT。 21 二、二维水动力弥散-瞬时投放示踪剂 直线法将高斯浓度转化为简单直线。 不足： 1.两个观测孔资料，其中一个位于x轴； 2.转化X、Y数据时，计算量大 22 二、二维水动力弥散-瞬时投放示踪剂 3.标准曲线法 若分子扩散可忽略，有 式 改成 若y=0（主流线上） 当 时， 因 ,故计算得u偏大。 23 二、二维水动力弥散-瞬时投放示踪剂 x是测得Cm的位置，此时 (5-42) (5-46) (5-43) (5-45) 24 二、二维水动力弥散-瞬时投放示踪剂 式中 CD是rD和tD函数，可作标准曲线求参 25 二、二维水动力弥散-瞬时投放示踪剂 对给定rD，K为定值。故CD-t与W-tD两曲线相似，可利用曲线拟合求参 由(5-42)得 26 二、二维水动力弥散-瞬时投放示踪剂 参数计算步骤： （1）若已知流向 在注入孔正下游设一取样孔（x1,0）,偏离正下游设孔（x2，y2） 27 二、二维水动力弥散-瞬时投放示踪剂 参数计算步骤： 1.作lgCD-lgt实测曲线； 2.与标准曲线拟合，得rD1和rD2； 3.根据（5-45）、（5-42）与（5-43），易知 若观测孔偏离流向，会出现弥散度偏低 28 二、二维水动力弥散-瞬时投放示踪剂 29 二、二维水动力弥散-瞬时投放示踪剂 参数计算步骤： （2）若流向不确定 计算方法改变，用（x1，y1）（x2,y2）两个观测孔 解得 30 二、二维水动力弥散-瞬时投放示踪剂 4.弥散晕面积求参法 前面已证 弥散晕为椭圆，圆心为（ut，0），以孔隙平均流速向前移动，长轴a和短轴b之比 以浓度C为等值线的椭圆面积为 (4-68) (5-61) 31 二、二维水动力弥散-瞬时投放示踪剂 4.弥散晕面积求参法 表明t时刻的一系列浓度与对应的椭圆面积在半对数坐标上为一直线，斜率为 求出横轴截距d，有（5-61）得 可求有效孔隙度 由（5-60）得 求出浓度C等值线围成的椭圆面积随时间变化值，则S/t与t在半对数坐标中呈直线。同理可求上述参数 联立，可得DL、DT (5-62) 32 二、二维水动力弥散-瞬时投放示踪剂 算例 在示踪剂浓度溶液中氯化物浓度为580.7mg/L，体积共0.7m3，分5井瞬时注入。由于含水层非均质，示踪剂分为流速快、慢移动。慢带