一维变参数溶质运移模型的数值研究

　　摘 要：研究地下水中溶质的迁移规律可以为地下水环境质量预报、控制和管理服务。首次采用数值方法探讨与位置坐标相关的溶质运移参数模型，采用即时时间域显式差分法格式双精度Fortran语言编程计算了变溶质运移参数模型的一维溶质运移问题，并进行了算法可靠性验证;然后对本文首次提出的3種与位置坐标相关的变溶质运移参数模型进行了参数分析，获得了这3种模型的初步认识。还将时间域显式差分方程计算结果与解析解和实验结果进行了验证对比分析，说明算法可靠，显式法列式简单，更适于复杂模型计算。本文提出的算法和处理技巧可以应用到许多与地质相关的溶质运移问题的研究，通过建立不同模式的溶质运移参数模型，从而实现复杂条件下溶质运移过程的预测和模拟。

　　关键词：一维溶质迁移;与深度相关;变溶质运移参数;建模与模拟;即时显式差分法

　　由于工农业生产、雨水径流冲刷和人为因素等作用，造成水流中溶质瞬间注入或连续注入从而引起河流或地下水的水质改变。研究溶质瞬时注入或连续注入情况下不同时刻溶质浓度在河流或地下水流中的分布状况、及不同位置时溶质浓度随时间的变化情况，可通过数学模型来模拟或预报水质在时间或空间的变化，达到为水环境质量预报、控制和水资源管理服务的目的[1-22]。

　　数值模拟水流溶质运移的主要方法有有限元法、边界单元法、有限差分法和其他数值法等。由于进行水质模拟过程中涉及的参数在不同的时间和空间上变化较大，在水质模拟过程中，会造成数值弥散、数值波动和计算量巨大等问题。国内主要的研究有徐玉佩[2-3]用Galerkin法推导了二维地下水运动的有限元方程和溶质运移有限元方程，并编制了计算程序应用于明沟排水和田面保持淹灌水层情况下的冲洗种稻改良盐碱地的分析，还与现有解析解和实验资料进行了对照。徐玉佩还用水动力弥散理论进行了地下水中水盐动态的研究，可直接用于盐碱地的预防和改良。黄康乐[4]则提出了一种求解二维饱和-非饱和溶质运移问题的交替方向特征有限单元法(ADCG)，在克服数值弥散、数值波动和提高计算速度方面具有显著的优越性。

　　张效先等[5]提出了一种计算地下水中溶质运移的数值法，将多维问题分解成几个一维问题， 减少了计算工作量， 并消除了数值波动现象，具有较高的精度。冯绍元[6]提出了特征有限元法的基本理论和实施过程，其计算实例表明能消除大部分的数值弥散而又不增加解的振荡性。姚磊华[7]首先基于分步求解的思想， 用广义迎风对偶单元均衡法求解对流定解问题， 对扩散定解问题采用一般的Galerkin有限元法求解， 不仅避免了数值弥散和过量问题， 而且避免了求节点速度这一步， 简化了运算步骤。邱克俭等[8]结合拉格朗日观点的运动坐标与欧拉观点的固定坐标，分别求解对流-弥散方程的对流效应和弥散效应， 并通过与理论解和试验结果做比较，验证了模型的可靠性。

　　任理等[9]采用混合拉普拉斯变换有限单元法求解对流占优的地下水溶质运移问题，能有效地消除数值扩散和过量现象，特别适于大区域地下水污染的长期预报。李焕荣等[10-11]利用广义差分法建立了一维非饱和水流问题的守恒形式的数值模型， 具有计算量小和稳定性好的特点， 其提出的非粘性土壤水中溶质运移问题的守恒混合元格式，精度高且数值稳定。段德宏等[14]以黄河兰州段为例，采用一维稳态河流水质模型，对岸边19个污水排放口及共约32.9 km长的全河段污染状况成功地进行了模拟。

　　徐文彬等[15]根据数值随时间步长变化的特点，提出了针对数值波动、弥散和过量现象的新改进方案，并用实例验证了方法的可行性。何丕文[16]利用一维水质模型，采用有限差分法求解河流中污染物浓度，以淮河淮南段为例进行了验证。马东豪等[17]研究了土壤溶质迁移的两流区与两区模型，结果表明与两区模型相比，两流区模型可以更好地描述优先流情况下的溶质穿透曲线。对于无优先流情况， 连续流输入情况下两个模型的参数一致性较好;而在脉冲输入时一致性较差。国外在溶质运移模型及其数值计算方法上的研究也很多，如JAISWAL等[18]运用Chebyshev配点法求解一维溶质运移问题的数值解。SILAVWE等[19]讨论了某些数值方法用于一维溶质运移模型的参数估计，通过将数值解与试验结果比较，反演溶质运移模型参数。JAISWAL等[20]采用算子矩阵法对多孔介质的非线性微分方程进行了数值分析，该法应用切比雪夫多项式并结合切比雪夫导数运算矩阵和谱配置法对解进行逼近。

　　其优点是可以将这些问题转化为易于求解的代数方程组。CARR[21]利用拉普拉斯变换导出了层状多孔介质溶质运移的新半解析解，该法在相邻层的界面引入表示运移量的未知函数，将多层问题置于拉普拉斯域上分层解决，然后再将其数值反演回时域。ALLWRIGHT等[22] 对地下水溶质运移方程在对流明显强于扩散时出现的数值波动提出了两种新的数值近似计算方法，即只适用于对流项的迎风Crank-Nicolson格式和加权迎风-下风格式，并对这些新提出的格式进行了数值稳定性分析和比较。发现如果Crank-Nicolson格式只用于对流项时很合宜。此外，显式加权迎风-下风有限差分格式是对传统显式一阶迎风格式的改进，而隐式加权一阶迎风-下风有限差分格式在给定适当的权重因子时无条件稳定。

　　本文主要贡献在于对一维溶质对流弥散基本方程，首先，建立了一维溶质运移3种变溶质运移参数模型，即线性变化模型/多项式变化模型、指数变化模型和幂变化模型;其次，采用显式时域差分法成功地进行了溶质浓度时空分布的数值计算，有效克服了很多研究人员遭遇的数值波动、数值弥散和计算量大等难题，并与解析解和实验结果做了对比，验证了该数值方法的正确性;最后，对3种溶质运移参数模型进行了参数分析，获得了3种模型的初步认识。

　　1 基本原理

　　研究溶质在水流中的迁移规律，试验测定可采用电解质脉冲示踪法等来进行。电解质脉冲在水流中的运移规律符合溶质对流弥散基本方程，根据夏卫生等[12-13]的研究，其溶质运移参数随水流位置而变化，在大规模的地下水流问题中，溶质运移参数亦会随空间位置而变化。对其溶质运移规律的模拟可采用近似解析解和数值解。

　　推荐阅读：自然科学研究人员发表论文期刊

期刊VIP网，您身边的高端学术顾问

文章名称： 一维变参数溶质运移模型的数值研究

文章地址： http://www.qikanvip.com/lunwen/zonghelunwen/2021/0707/58330.html