在这里插入图片描述

饱和地下水流和溶质运移常见的模拟软件

常用的求解地下水水流和溶质运移方程的数学方法有两种：有限差分法和有限元法。两者主要的差别在于离散模型区的方法不同。基于不同的数学方法，当前市场上有一些地下水模拟图形用户界面，它们在基本功能和软件包上基本相似，但是又各有特点。建模软件一般根据模拟目标和软件功能需求进行选择。饱和地下水流和溶质运移常见的模拟软件有Visual MODFLOW、GMS和Feflow。

Visual MODFLOW：

Visual MODFLOW是行业标准的用于地下水流和污染物运移模拟的三维地下水模拟软件。该软件基于矩形网格的有限差分法，无缝集成MODFLOW-2000、MODFLOW-2005、MODPATH、MT3D99、MT3D、WinPEST等软件的基础上，建立了合理的Windows 菜单界面与可视化功能，增强了模型数值的可操作性。界面设计包括三大彼此联系但又相对独立的模块，即前处理模块、计算模块和后处理模块。现有Visual MODFLOW 4.1和Visual MODFLOW Flex 2014.1中文版，以及Visual MODFLOW Flex 2015.1英文版，以满足客户的专业需求。

GMS

由Brigham Young大学环境模拟研究实验室（Environmental Modeling Research8Laboratory）开发的基于概念模型的地下水环境模拟软件。GMS在环境保护、水资源利用与管理、采矿、建筑等许多行业和部门得到了广泛的应用，成为最为普及的地下水运移数值模拟的计算软件。GMS基于矩形网格的有限差分法，无缝集成MODFLOW-2000、MODFLOW-2005、MODPATH、MT3D、WinPEST、PHT3D等软件。GMS 全面包容了模拟地下水流每一个阶段所需的工具，如边界概化、建模、后处理、调参、可视化。GMS是惟一支持TINs、Solids、钻孔数据、2D或者3D地质统计学的系统。现有的最近版本为GMS10.0。

FEFLOW

FEFLOW是由德国WASY公司开发的，该软件包具有图形人机对话、地理信息系统数据接口、自动产生空间各种有限单元网、空间参数区域化及快速精确的数值算法和先进的图形视觉化技术等特点。在FEFLOW系统中，用户可以很方便迅速地产生空间有限三角网格单元，设置模型参数和定义边界条件，运行数值模拟以及实时图形显示结果与成图。

常用地下水与溶质运移模拟软件说明表

序号	模型名称	维数	功能	模型解类型	学习模型需要的时间	获取途径	研发单位（个人）	备注
1	AQTESOLV EPA-WHPA	二维	半解析解模型，用于确定含有多个注水井和抽水井的均质、二维地下水稳态流的捕获区。它包括Monte-Carlo模拟参数影响模块。	解析解水流模型	所用天数适中		国际地下水模型中心
3	WALTON		WALTON35 是由35个简单的水流、溶质或热迁移及淡水-盐水界面解析解与数值解模型组成。	解析解水流溶质运移模型
4	AT123D	三维	点、线、面或体积源在恒定、水平均质流条件下发生溶质迁移，能模拟放射性核素、化学物质或热废料等。	解析解溶质运移模型	所用天数适中		国际地下水模型中心；GEMS	被广泛应用，3D版本存在一些问题。
5	EPA-WHPA	二维	WHPA是一个地下水流模拟程序，用于勾勒井口保护区捕获区。它可应用于二维稳态水流、均质的含水层。	解析解溶质运移模型
6	MAS3DB	三维	二元高斯在y、z平面源无限、均匀的流场	解析解溶质运移模型	仅需几个小时		滑铁卢水文地质公司
7	NODK2D, 3D DK2D,3D	二维、三维	反演解析解模型（衰减可选），能解决弥撒度、源维数、浓度及持续时间，筛选因子等问题。	解析解溶质运移模型	所用天数适中		P.Domenoco. Texas A&M	不能不加区别的使用
8	PLUME3D	三维	用于点源溶质迁移，均匀流场	解析解溶质运移模型	仅需几个小时		国际地下水模型中心
9	PRINCE	一维、二维及三维	PRINCE是一个软件包，包括10个解析解模型。其中，7个一维、二维与三维质量迁移模型，3个二维水流模型。该软件包界面友好，制图功能强大，能导入数字或DXF位图文件。一维、二维与三维对流弥散质量迁移模型包含一级衰减和线性延迟；二维条带和高斯源质量迁移包含对流弥散、一级衰减和线性延迟；多溶质注入井采用Wilson-Miller方法；多注水-抽水井采用二维各向异性泰斯方法；渗漏、各向异性承压含水层采用二维Hantush-Jacob方法。	解析解溶质运移模型	所用小时数适中	商业软件	滑铁卢水文地质公司；科学软件集团	用户界面友好，软件包驱动
10	SLUG	二维	均匀流场中微水溶质运移	解析解溶质运移模型	仅需几个小时		国际地下水模型中心
11	SLUG3D	三维	均匀流场中微水溶质运移	解析解溶质运移模型	仅需几个小时		国际地下水模型中心
12	WMPLUME	二维	点源、均匀流场中溶质迁移		仅需几个小时		国际地下水模型中心
13	CFEST	三维	CFEST是一个用于各向异性、非均质、承压/半承压、潜水含水层的三维有限元技术的地下水流和迁移程序，包括对流、弥散、扩散、吸附、衰减、盐溶解及渗滤等。	数值解饱和流模型
14	HST3D	三维	用于模拟三维地下水流及热与溶质运移的计算机代码	数值解饱和流模型
15	MODFLOW	准三维	MODFLOW是一个有限差分模型，用于模拟承压含水层、潜水含水层及潜水与承压结合含水层的三维水流。也能模拟与水流相关的外部压力，如井区补给、蒸散发、排泄及溪流作用等。	数值解饱和流模型	需要很多天		国际地下水模型中心，美国地质调查局或科学软件集团	被广泛应用，能与MOC,MODPATH,MT3D,PATH3D等结合进行运移计算
16	FEMWATER	二维、三维	使用同化能力标准划分井口保护区。模拟地下水流考虑一下几点：（1）考虑到各种地质构造的均质、各向异性介质；（2）考虑分布点的、随时间变化的源汇项；（3）兼容4中类型的边界条件。	数值解饱和流模型	需要几个月时间		美国橡树岭国家实验室	能与FEMWASTE结合进行运移计算
17	AQUA	二维	瞬态非等温、非均质、各向异性含水层的水流和运移；运移过程包括对流、平流、吸附和衰减过程。	数据解饱和水流与溶质迁移模型	需要几个星期		科学软件集团	包括交互式图形，数据生成选项
18	BIOPLUME Ⅱ	二维	模拟在有限氧的生物降解、复氧和厌氧生物降解（作为一级衰减）的影响下溶解性碳氢化合物的运移。模型建立在USGS 2-D溶质迁移程序基础上，模拟对流、弥散、混合机生物降解。可以计算每个时间步长碳氢化合物与氧气分布羽。氧气与碳氢化合物的瞬时反应是假设用来模拟天然生物降解过程，延迟羽，及现场生物修复计划。另外，注水井可认为是氧源。	数据解饱和水流与溶质迁移模型	需要几个小时		NTIS
19	CHEMFLO	一维	模拟非饱和土壤中水流与化学物质运移。水流模拟使用Richard’s方程。污染物迁移使用对流弥散方程，模拟结果可以通过水量导水率、密度通量、化学物浓度随时间和空间变化的图形表示。土壤剖面假设均质。	数据解饱和水流与溶质迁移模型	需要几个小时		美国环保局
20	FLONET: FLOWTRANS	二维	用于二维、各向异性、垂直剖面非均质、承压或潜水含水层中稳态流与瞬时流运移。另外，可以模拟生物降解和吸附。	数据解饱和水流与溶质迁移模型	需要几天时间		滑铁卢水文地质公司；科学软件集团
21	FLOWPATH	二维	用于各向同性、非均质、承压、渗漏或潜水含水层中稳态流与迹线；计算迹线、迁移时间、流速与捕获区。	数据解饱和水流与溶质迁移模型	需要几天时间		滑铁卢水文地质公司；科学软件集团	应用较广泛
22	FTWORK	三维	用于各向异性、非均质、承压、潜水含水层中稳态流或瞬时流与运移。运移过程有对流、弥散、吸附与衰减。对于稳态流可进行参数评估。	数据解饱和水流与溶质迁移模型	需要几个星期时间		Geo Trans	已经应用于低放射性废物评估区治理
23	Golder Groundwater Package	二维	用于分层含水层系统中瞬态流和化学物运移	数据解饱和水流与溶质迁移模型	需要几天时间		Golder Associates	在澳大利亚被广泛应用
24	GWTHERM	二维	用于各向异性、非均质、潜水含水层的瞬态、非等温流。包括对流、弥散、扩散、吸附、衰减及颗粒示踪等选项。	数据解饱和水流与溶质迁移模型	需要几天时间		Dames & Moore	用于核废料库
25	HST3D	三维	模拟三维地下水流与热和溶质相关的溶质运移。通过裂隙孔隙速度，流体密度与温度、压力及溶质质量关系、及溶质粘度对温度、溶质量关系等控制方程耦合在一起。溶质运移方程是单一溶质组分对线性平衡吸附和线性衰减。	数据解饱和水流与溶质迁移模型	需要几天时间		国际地下水模型中心	有时计算时间较长
26	MOCDENSE	二维	USGS2D-MOC的密度耦合版	数据解饱和水流与溶质迁移模型	需要几天时间		国际地下水模型中心
27	MT3D	"二维
三维"	用于各向异性、非均质、承压、半承压或潜水含水层中稳态流或瞬时流。运移过程包括对流、弥散、非线性吸附与衰减。需要由中心差分水流模型获得结果，比如MODFLOW 作为输入。	数据解饱和水流与溶质迁移模型	需要几个星期		国际地下水模型中心；Papadopoulos&Associates
28	PORFLOW	三维	用于各向异性、非均质裂隙的承压半承压含水层的瞬态非等温水流。运移过程包括对流、弥散及扩散、吸附、衰减、多相运移及溶解。	数据解饱和水流与溶质迁移模型	需要几个星期		国际地下水模型中心
29	SWENT	三维	用于各向异性、非均质、承压含水层的瞬态、非等温及密度耦合水流与运移。模拟运移的过程包括对流、弥散、吸附、衰减、盐溶解及淋滤等。	数据解饱和水流与溶质迁移模型	需要几个星期		Intera Environmental Consultants
30	SWICHA	二维、三维	用于各向异性、非均质、承压含水层中瞬态,密度耦合的水流与运移。运移过程包括对流、弥散及扩散。	数据解饱和水流与溶质迁移模型	需要几个星期		Geotrans
31	SWIFT Ⅱ	三维	用于各向异性、非均质裂隙、承压、半承压或潜水含水层中瞬态、非等温、密度耦合的水流与溶质运移。运移过程包括对流、弥散、扩散、吸附、衰减、盐溶解、淋滤及双孔隙度。	数据解饱和水流与溶质迁移模型	需要几个月时间		SNL Geotrans	应用于热储存、地下水污染、核废物隔离及地下深注工程
32	TARGET	二维、三维	用于各向异性、非均质裂隙、承压、半承压或潜水含水层中瞬态、、密度耦合的水流与溶质运移。运移过程包括对流、弥散、扩散、吸附及衰减。包括湖泊水位预测模块。	数据解饱和水流与溶质迁移模型	二维：需要几天时间；三维：需要几个星期		国际地下水模型中心	在Peclet数较大时不稳定
33	USGS2D-MOC	二维	模拟饱和地下水系统中非保守溶质运移。它计算由于对流迁移、水动力学弥散、混合、或补给引起的溶解及化学反应等造成的不同时空浓度分布。化学反应包括一级不可逆速率反应、不可逆平衡控制的线性、Freundlich 或Langmuir 等温吸附，离子交换控制的可逆平衡。含水层可以是非均质的，或者各向异性。	数据解饱和水流与溶质迁移模型	需要几天时间		国际地下水模型中心
34	VLEACH	一维	模拟包气带中挥发性、吸附的污染物的渗滤，其模拟4个主要过程：液相对流、固相对流、固相吸附及气相平衡。该程序基于一些假设，包括稳态非饱和带流，无降解等。	数据解饱和水流与溶质迁移模型			美国环保局
35	CHAMP	三维	用于各向异性、非均质、裂隙、变形多孔介质的瞬态的水流与运移过程。运移过程包括对流、弥散、扩散、吸附及衰减等。	数值解可变饱和水流与溶质运移模型	需要几个月时间		LBL	用于废物处置中的化学物运移
36	FEMWASTE、3DLEWASTE	"一维
二维
三维"	用于各向异性、非均质多孔介质的瞬态、密度/含水量耦合的运移模型。运移过程包括对流、弥散、扩散、吸附及衰减（含3个吸附模型）。	数值解可变饱和水流与溶质运移模型	需要几个月时间		美国橡树岭国家实验室
37	GS3	三维	用于各向异性、非均质、承压或潜水含水层中含有自由地表蒸发的瞬态流和溶质运移模型。运移过程包括对流、弥散及扩散。	数值解可变饱和水流与溶质运移模型	需要几个星期时间		Bechtel	应用于饱和与非饱和带化学物质迁移
38	MOTIF	三维	用于各向异性、非均质、承压或潜水、裂隙或多孔介质含水层中瞬态流和密度耦合的溶质运移。运移过程包括对流、弥散、扩散、吸附及衰减。	数值解可变饱和水流与溶质运移模型	需要几个月时间		加拿大原子能公司
39	NAMMU	三维	用于各向异性、非均质、承压或潜水、裂隙或多孔介质含水层中瞬态流和密度耦合的溶质运移。运移过程包括对流、弥散、扩散、吸附、衰减及放射性核素链衰变。	数值解可变饱和水流与溶质运移模型	需要几个月时间		英国原子能管理局Harwell 实验室
40	SUTRA	二维	用于各向异性、非均质、多孔介质含水层中瞬态流和密度耦合、非等温水流和溶质运移。运移过程包括对流、弥散、扩散、吸附及衰减（3个吸附模型）	数值解可变饱和水流与溶质运移模型	需要几个星期时间		国际地下水模型中心
41	SWANFLOW	三维	用于各向异性、非均质、承压或潜水含水层中瞬态、三相（水、气、非水相液体）的水流与溶质迁移密度耦合、非等温水流和溶质运移。运移过程包括对流、弥散、扩散、吸附及衰减（3个吸附模型）	数值解可变饱和水流与溶质运移模型
42	TOUGH2	三维	用于各向异性、非均质、裂隙或多孔介质中瞬态、非等温的水流与溶质运移。运移过程包括对流、扩散、溶解及多项流。	数值解可变饱和水流与溶质运移模型	需要几个月时间		LBL
43	TRACR3D		用于各向异性、非均质、裂隙介质中瞬态水流与溶质运移。运移过程包括对流、弥散、扩散、吸附、衰减、多组分、二相流及随机参数输入等。	数值解可变饱和水流与溶质运移模型	需要几个月时间		LNAL
44	3DADE	三维	用于评价一系列的三维对流弥散溶质迁移方程的解析解。迁移方程包括项：对流弥散、溶质延迟、一级衰减、零级沉积。该程序能直接解决下列问题：对于具体的模型参数，浓度计算结果为时间和空间的函数；对具体的实验室数据进行解析法拟合，以评估所选择的参数。程序采用Fortran计算机语言。	解析解饱和水流与溶质运移模型		免费	美国盐渍土实验室
45	AGU-10		一个可用于筛选均质、各向同性流场的水流与溶质迁移程序的程序集，其基于美国地球物理协会的水资源专著10。其由5个Fortran语言编写的模拟程序和2个Microsoft Basic的预/后处理器组成。	解析解饱和水流与溶质运移模型		免费	美国地球物理协会
46	AT123D	一维、二维、三维	AT123D是一个瞬态的一维、二维、三维的解析解地下水迁移模型。其计算含水层系统中污染物浓度的时空分布，并地下水含水层中污染羽的瞬态扩散。溶质迁移过程包括对流、扩散、吸附及衰减。AT123D能解决两类污染源：瞬态源与连续源；三类污染源种类：放射性核素、化学污染物及热；四类污染源分布结构：点源、线源、面源及体积源；四种类型的含水层尺寸：有限深度与有限宽度、有限深度与无限宽度、无限深度与有限宽度及无限深度与无限宽度。	解析解饱和水流与溶质运移模型
47	CXTFIT		在可利用准确迁移参数基础上能成功预测地下水溶质迁移。其基于Levenberg-Marquardt算法对均质土壤一维稳定流动中溶质易混合置换实验所得到的穿透曲线(Breakthrough curve,BTC)进行非线性最小二乘拟合。	解析解饱和水流与溶质运移模型		免费	CXTFIT是由美国盐渍土实验室研制的计算机软件。
48	EPA-VHS		EPA-VHS筛选模型，用于预测连续污染源下游指定距离的最大污染物浓度。模型基于保守成分的解析解方法。	解析解饱和水流与溶质运移模型			美国环保局
49	HPS	三维	HPS是一个解析解模型，用于在模拟一个均质区域地下水中水平面源三维污染物迁移的模型。该模型嵌入对流弥散迁移、污染物延迟及一级衰减。	解析解饱和水流与溶质运移模型
50	MAP		MAP监测分析软件包是一个监测网络的设计工具。其由污染羽生成模型（PLUME）,监测效率模型（MEMO）及污染概率模型（COPRO）组成。	解析解饱和水流与溶质运移模型		商业软件
51	ONE-D		ONE-D是一个含有5个解析解模型的软件包，包括一维对流弥散方程、线性吸附、零级沉淀、一级衰减。控制方程项包括对流，弥散及扩散。	解析解饱和水流与溶质运移模型		商业软件
52	PLUME		PLUME（也称PLUME-PH）是一个解析解模型，用于计算连续的、随时间变化释放的水平线源中单个放射性或非放射性组分在均质含水层中点浓度。	解析解饱和水流与溶质运移模型
53	PLUME 2D		PLUME 2D是一个解析解模型，基于一个或多个示踪剂的瞬时或连续释放的非保守溶质迁移方程的封闭式解法。	解析解饱和水流与溶质运移模型
54	Visual MODFLOW	三维	Visual MODFLOW是综合已有的MODFLOW、MODPATH、MT3D、RT3D和WinPEST等地下水模型开发出来的可视化模拟系统，可进行三维水流模拟、溶质运移模拟（溶质化学反应、水岩相互作用等）。	数值解水流与溶质运移模型		商业软件	加拿大滑铁卢公司
55	RWH		RWH是一个广发使用RANDOM WALK 水流和溶质运移模型的分析版本，专门用于设计筛选和教育目的。其用于模拟均质、各向同性、等厚度、有限范围的承压含水层介质。	解析解饱和水流与溶质运移模型
56	RANDOM WALK		RANDOM WALK 用于研究溶质迁移问题，该模型比建立有限差分或有限元水流模型简单。	解析解饱和水流与溶质运移模型		商业软件
57	SOLUTE		SOLUTE软件包是5个程序的集合，该软件基于非保守示踪溶质的对流弥散运移解析解模型。	解析解饱和水流与溶质运移模型		商业软件
58	USGS-SOL	一维、二维、三维	USGS-SOL显示了一系列偏微分方程的解析解，其描述各种类型边界和一维、二维、三维的非保守溶质对流弥散运移。	解析解饱和水流与溶质运移模型
59	VERTPAK-1		VERTPAK-1包含流体流动、岩石变形、裂隙中溶质迁移、多孔介质中溶质运移等8个解析解模型。	解析解饱和水流与溶质运移模型		商业软件
60	WALTON35		WALTON35 是由35个简单的水流、溶质或热迁移及淡水-盐水界面解析解与数值解模型组成。	解析解饱和水流与溶质运移模型
61	ASM	二维	ASM(Aquifer Simulation Model)是解决地下水流与溶质运移问题的二维水平或垂直模拟的数值模型。其通过中心节点有限差分法求解水流方程。	数值解解饱和水流与溶质运移模型
62	BIOPLUMEⅢ		有限差分模型，用于模拟地下水中由于对流、弥散、吸附及生物降解等作用的有机污染物自然降解过程。生物转化过程在有机污染的含水层修复中非常重要。因此，这些过程对烃类污染物修复方案进行评价。	数值解解饱和水流与溶质运移模型		免费	美国环保局
63	FEFLOW		FEFLOW是一款用于模拟地下水流与溶解性成分或热的溶质迁移过程模型的专业软件包。	数值解解饱和水流与溶质运移模型		商业软件	DHI GROUP
64	FLOWPATH Ⅱ	二维	FLOWPATH 是一个二维有限差分、稳定流地下水流模型，用于计算井口保护区、水头、地下水流速、时间迹线、捕获区、水平衡及稳定状态降深分布。	数值解解饱和水流与溶质运移模型
65	FTWORK	一维、二维、三维	FTWORK 程序是一个以块为中心的有限差分模型，用于模拟承压含水层与潜水含水层中饱和介质的一维、二维与三维的稳定流和瞬态流的溶质运移。	数值解解饱和水流与溶质运移模型
66	GMS	二维、三维	GMS是一个用于地下水模拟各阶段的大软件工具包，包括评估区刻画、模型建立、后处理、律定及可视化等。GMS是唯一一个支持三角不规则网、立方体、钻孔数据、二维与三维地质统计以及二维与三维中有限元和有限差分模型的系统。目前，其支持的模型包括MODFLOW,MODPATH,MT3D,RT3D,FEMWATER与SEEPSD.	数值解解饱和水流与溶质运移模型			Brigham young university
67	Groundwater Vistas	三维	Groundwater Vistas是一个复杂的地下水三维水流与溶质运移的Windows图形用户界面，其将全面的图形分析工具和模型设计体系耦合在一起。	数值解解饱和水流与溶质运移模型
68	HST2D/3D ENHANCED		HST 2D/3D是一个用于模拟饱和含水层多孔介质中热与溶质迁移的有限差分模型。	数值解解饱和水流与溶质运移模型			美国地质调查局
69	HST3D	三维	用于模拟三维地下水流及热与溶质运移的计算机代码。	数值解解饱和水流与溶质运移模型		免费
70	HYDRUS 2D/3D	二维、三维	基于微软Windows环境，用于模拟二维和三维的变饱和介质中水流和溶质迁移的软件。图形用户界面颜色丰富，精度高。	数值解解饱和水流与溶质运移模型		商业软件
71	MIGRATE	二维	该软件很快且准确地模拟二维地表或地下的多个污染源溶质运移。模型考虑对流弥散迁移、吸附、放射核素与生物降解及裂隙迁移等。	数值解解饱和水流与溶质运移模型		商业软件
72	MOC		一个二维有限元程序，用于模拟裂隙或颗粒含水层中地下水流和溶质运移，并能处理渗漏承压与水位系统。	数值解解饱和水流与溶质运移模型		免费	美国地质调查局
73	MOC3D	三维	MOC3D是一个三维水流与溶质运移模型。其将Modflow与运移模型刻画集成，并考虑到对流、弥散、流体源混合、线性吸附及放射性核素衰减等。	数值解解饱和水流与溶质运移模型		免费	美国地质调查局
74	MOCDENSE	二维	MOCDENSE是一个二维横断面模型，用于分析盐水入侵。其模拟地下水系统中密度流中保守的一个或两个组分的溶质迁移扩散。	数值解解饱和水流与溶质运移模型		免费	美国地质调查局
75	MODFLOWT	三维	MODFLOWT是MODFLOW的一个新版本，其包括模拟三维溶质运移模块。	数值解解饱和水流与溶质运移模型			GeoTrans，Inc
76	MT3DMS	三维	MT3DMS模型是一个三维溶质迁移模型模块。其能够模拟地下水系统中的对流、扩散、衰减、溶质化学反应、线性与非线性吸附作用等，能够对承压含水层、不承压含水层、承压与不承压交替含水层以及倾斜的和单元厚度变化的含水层进行空间离散。	数值解解饱和水流与溶质运移模型			Authored by Zheng at S. S. Papadopulos & Associates, Inc.
77	MULAT	三维	MULAT是一个分析三维稳态地下水流与颗粒对流运移的程序，其包括吸附与解吸。模型假设含水层中水流主要是水平流，垂直流由一个含水层通过阻水层到另一个含水层。	数值解解饱和水流与溶质运移模型
78	PHREEQM-2D		PHREEQM-2D是一个多组分溶质迁移模型，其由两个耦合的模拟程序组成：HST2D与PHREEQE。PHREEQE是一个地球化学反应模型（USGS编写），PHREEQM-2D包括溶解、沉淀及离子交换，支持矿物动力学计算。	数值解解饱和水流与溶质运移模型			美国地质调查局
79	PMWIN	三维	PMWIN是一个模拟三维水流和溶质迁移的图形界面。它支持文件输入的准备工作和分析MODFLOW相关程序（包括MODFLOW本身，离子示踪程序MODPATH和PMPATH,运移程序MT3D及参数估计程序PEST）的结果。	数值解解饱和水流与溶质运移模型