三水动力弥散方程的解析解法.ppt

* 半无限长多孔介质柱体，一端为定浓度边界 该式可用于地表水体。如一条均匀的长渠道，在x=0处定浓度C0,并以稳定速度u流动，只需讲Dl改成Dm 余补误差函数erfc（η）随着η的增大而减少，当x足够大或t足够长，右端第二项可忽略不计，即 * 瞬时注入示踪剂-平面瞬时点源 均质各项同性、等厚的承压含水层中存在一维稳定流动，孔隙平均流速为u，取x坐标轴平行地下水流向，产生，在x方向为纵向弥散系数DL，在y方向为横向弥散系数DT。 * 三、一维稳定流动二维水动力弥散问题 假定（1）示踪剂的注入不改变地下水渗流状态 或（2）取样点离注入井足够远 数学模型为 * 瞬时注入示踪剂-平面瞬时点源 通过动坐标以及变换x、y坐标尺度的方法，与基本解产生联系 令 则 同理 得 * 瞬时注入示踪剂-平面瞬时点源 引入动坐标 令 记 套用基本解，有 整理得 * 瞬时注入示踪剂-平面瞬时点源 因DLDT,a是x的长半轴 当t与C为定值时，上式为常数，记为-A，并设X=x-ut，上式变为 为中心坐标（ut，0）, 的椭圆 * 瞬时注入示踪剂-平面瞬时点源 某浓度等值线所围面积随t的变化，先是增大然后变小。低浓度所围面积随时间增大的持续时间长，而高浓度持续时间短 * 连续注入示踪剂-平面连续点源 连续点源的作为视为无数瞬时点源之和 设单位时间注入示踪剂的质量为ml(=C0Q)，dt’注入示踪剂质量为mldt’有 经两重换元并化简后，得 * 连续注入示踪剂-平面连续点源 当t较长，简化为 平面稳定连续注入点源的解 * 注入拟稳定条件下示踪剂的径向弥散 设在水平、等厚（B）、无限展布的均质各项同性承压含水层中有一口完整井，井径rw。通过井向其连续注入定流量Q且示踪剂浓度C0的水。忽略天然流速，井的附近形成拟稳定二维径向流。 以井为中心，任意半径为r的圆周通量 平均流速 * 三、一维稳定流动二维水动力弥散问题 转换成极坐标，令 代入式 有 对于均质各向同性介质，无天然流速，弥散是对称的 即 * 三、一维稳定流动二维水动力弥散问题 其中 有 若忽略分子扩散模型，则 其近似解 用于确定实测的纵向弥散度 空间瞬时点源 渗透系数K为均质各向同性， 若存在一维稳定流动，流速为 。某点o处瞬时注入示踪剂m。 取o点为坐标原点，x轴平行于 ，且方向相同。 * 三、一维稳定流动三维水动力弥散问题 对弥散系数D来说是各向异性的，它属于二度各向异性，即 ，弥散方程写成 * 三、一维稳定流动三维水动力弥散问题 通过坐标变换变成各向同性，令 * 三、一维稳定流动三维水动力弥散问题 则 有 得 采用动坐标，令 ，记 * 三、一维稳定流动三维水动力弥散问题 方程改成 套用基本解，得 进行坐标反变换，得 * 三、一维稳定流动三维水动力弥散问题 讨论： （1）随时弥散维数的增加，浓度C衰减速度也加快。 （2）对式进行变换，得 椭球方程。等浓度面为一个旋转椭球面，呈橄榄球桩，长轴沿x方向 * 三、一维稳定流动三维水动力弥散问题 空间稳定连续点源 假定：示踪剂注入并不改变渗流场的原始特征，即示踪剂是理想的，且保持原来的一维稳定流动。 设单位时间注入示踪剂的质量为ml 连续点源视为无数瞬时点源组成，注入时间记为t’。 t’时刻于坐标原点处注入示踪剂质量 空间瞬时点源的解： * 三、一维稳定流动三维水动力弥散问题 空间稳定连续点源 解得 （具体求解过程见P52-54） 式中 * 三、一维稳定流动三维水动力弥散问题 当 故 式子变成 地下水溶质运移理论及模型 * 一、基本解 基本解 将瞬时注入点源问题的解称为基本解。由基本解出发，利用叠加原理到处线源、面源、多点源及连续注入问题的解。 三维空间瞬时点源 （1）均质各向同性； （2）静止流场 ，弥散系数为常数，即 ，流体密度为常数； （3）t=0时，在原点处瞬时注入质量为m的溶质； （4）瞬时点源位置为坐标原点； * 一、基本解 对流弥散方程简化成 D表示多孔介质分子扩散系数 浓度C对称于原点分布 取半径为R和R+dR的两个球面所构成的单元体为均衡段，根据质量均衡得 * 一、基本解 问题写成 略去高阶变量 * 一、基本解 问题写成 略去高阶变量 * 一、基本解 根据因次分析中的π定理设 将m、n合并成新变量m/n，得 和 对该问题，有两个独立的π参数，依π定理有 π1、π2可有多种组合，但上述组合可得到最简单的常微分方