分析二维水动力弥散试验数据的线性回归法.pdfVIP

维普资讯 西北水资源与水工程 第 l0卷 l999年 第 l期 Water ~esourees ＆ Water Engineering 分析二维水动力弥散试验数据的 ，f c 线性回归法 — 郭 建 青 马 致 远 ， (西安工程学院水工系 西安 710054) 摘 要 文中将描述瞬时投放示踪剂，一维流动二维地下水弗散过程的解析解进行 了适 当数学变换 ．得到一个二元线性 回归方程。该方程中的自变量和因变量 中仅包古着试验数据 ， 回归系数中包含着待求参数。通过对两个观测孔中的观测数据进行分析．可以得出五个包古 待求参数的方程。由这些方程，推导出了计算 、t、“、 和 五个参数值的公式。线性回归 法的主要优点为在示踪剂不论是否为惰性的情况下均可以应用。 关键词 弥散试验 参数确定 线性回归 确定主要是通过分析实际弥散试验数据而实现的。目前已有诸如试算法、直线图解法、元量 纲标准曲线法[、逐点求参法、弥散晕面积法和优选拟台法[2j_等分析一维流动二维弥散、瞬 时投放示踪剂条件下的弥散试验数据，确定含水层介质纵向与横向弥散度的方法。在这些方 法中，除优选拟台法外，其余方法都不能在示踪剂为可降解或可衰变情况下应用。也就是说 它们仅能在示踪荆为惰性的情况下应用。另一方面，从理论上来讲 ，优选拟合法可以在示踪 剂为非惰性的情况下应用 然而，在欲求参数个数较多时，适当地选择参数初始值并非容易。 因此，作者在文中提出一种在示踪剂不论是否为惰性情况下均可以应用的数据分析方法 ，其 能够同时直接计算含水层介质的纵向和横向弥散度 和卵、地下水流速 “、含水层的有效孔 隙率 和示踪剂的一级反应速率系数 ^等五个参数值 1 弥散试验 所谓一维流动二维水动力弥散试验是在如图1所示的一维流场 (天然或人工形成)内的 投放井中，瞬时投放一定质量的示踪剂，取投放时刻为流场 中示踪剂浓度测定起始时刻，按 试验设计要求的时问间隔测定观测并中示踪剂浓度值 ，由此可以得到各个观测井中示踪剂 浓度随时间的变化过程 收稿 刚朝：1998—0914 — 14 一 维普资讯 在含水层为水平等厚 ，均质各 向同性与 无限延伸的条件下，若取投放点为坐标原 yI 点，水流方向为x轴方向，垂直水流方向为 I ．(．y) 轴方向，则描述在投入示踪剂后，示踪剂浓 。 l 度在时空上分布的解析表达式为。： 彝。t_——————— —_：I c：— 塑 一 … p[一 l 4矶～ 乱“ I 一 — L - ]’ p(一 ) (1) 4aLut 式中：c一 在时间f、流场 中点 、)处的示 围1 一维流动二雏弥散示意图 踪剂浓度值 (mL )；△ 一 瞬时投入示踪剂的质量 (m)；n一 含水层的有效孔隙率 (无量 纲 )；m一含水层厚度 (L)；“一 地下水流速 (LT )；dL一 含水层介质的纵向弥散度 (L)； 钾一 含水层介质的横向弥散度 (L)； 一时问坐标 (T)；、一 空间坐标 (L)；——示踪剂 的一级反应速率常数 (T )，在示踪剂为可降解情况下，其代表降解系数；在示踪剂为可衰 变情况下，其代表衰变系数。 2 计算原理与方法 为了能够以式 (1)为基础，建立分析弥散试验数据的方法，对式 (1)将进行如下数学变 换 在式 (1)两端同乘以f，然后取对数可以得出： 1n( ) (— )一 一 一 4 H “V r 然后在式 (2)两端再同乘以f，并经化简与整理后可得： (fc)一一