2、土壤污染模拟
土壤是一个多相的疏松的多孔介质，固相中有大量的有机和无机胶体。石油是一种天然的粘油状液体，主要成分为烃类化合物(占80%一90%) 。烃类化合物是非极性有机物，其偶极矩
式中:Cs 为平衡时固相中的浓度(mg/kg)；Ce 为平衡时液相中的浓度(mg/l)根据土壤中溶质运移模型和石油类污染物质在土壤中的迁移转化过程，考虑吸附作用而忽略石油的挥发，建立石油类污染物质在土壤中迁移转化二维综合模型。它包括水运动方程和石油运动方程。
土壤中水运动方程:
土壤中石油类运动方程:
式中：C(h) 为比水容量（cm -1）；K x 、K z 分别为横向纵向水力传导系数（cm/d）；Dxx 、Dzz 分别为横向纵向弥散系数（cm 2/d）；Rd 为滞留因子；c 为液相中石油的浓度(mg/l)；qx 、qz 分别为x 和z 方向的达西流速（cm/d）；θ为含水量(%) ；λ为降解系数（d -1）；h 为土壤中压力水头(cm) 。
初始条件和边界条件
根据监测的结果和落地油的分布特征，预测石油类在土壤中迁移过程及石油是否会对地下水造成污染，选择预测范围为：长80m ，深6m 剖面区域。并对部分问题可进行简化处理, 作一些基本假设。假设土壤水最初不含石油，即未受到污染，但土壤中存在一定的本底值，经取样测定取平均值为40.3mg/kg。在土壤的预测范围内，土壤被认为是均质的。
对水运动方程上边界确定为Cauchy 边界，下边界为Neumann 边界。
初始条件：
边界条件：
式中，ε为降雨-蒸发系数（mm/d）
对于溶质运动方程的边界条件和初始条件确定为：初始条件：
边界条件：
式中，n x 、n z 分别为土壤中x 、z 的外法向向量。
2、土壤污染模拟
土壤是一个多相的疏松的多孔介质，固相中有大量的有机和无机胶体。石油是一种天然的粘油状液体，主要成分为烃类化合物(占80%一90%) 。烃类化合物是非极性有机物，其偶极矩
式中:Cs 为平衡时固相中的浓度(mg/kg)；Ce 为平衡时液相中的浓度(mg/l)根据土壤中溶质运移模型和石油类污染物质在土壤中的迁移转化过程，考虑吸附作用而忽略石油的挥发，建立石油类污染物质在土壤中迁移转化二维综合模型。它包括水运动方程和石油运动方程。
土壤中水运动方程:
土壤中石油类运动方程:
式中：C(h) 为比水容量（cm -1）；K x 、K z 分别为横向纵向水力传导系数（cm/d）；Dxx 、Dzz 分别为横向纵向弥散系数（cm 2/d）；Rd 为滞留因子；c 为液相中石油的浓度(mg/l)；qx 、qz 分别为x 和z 方向的达西流速（cm/d）；θ为含水量(%) ；λ为降解系数（d -1）；h 为土壤中压力水头(cm) 。
初始条件和边界条件
根据监测的结果和落地油的分布特征，预测石油类在土壤中迁移过程及石油是否会对地下水造成污染，选择预测范围为：长80m ，深6m 剖面区域。并对部分问题可进行简化处理, 作一些基本假设。假设土壤水最初不含石油，即未受到污染，但土壤中存在一定的本底值，经取样测定取平均值为40.3mg/kg。在土壤的预测范围内，土壤被认为是均质的。
对水运动方程上边界确定为Cauchy 边界，下边界为Neumann 边界。
初始条件：
边界条件：
式中，ε为降雨-蒸发系数（mm/d）
对于溶质运动方程的边界条件和初始条件确定为：初始条件：
边界条件：
式中，n x 、n z 分别为土壤中x 、z 的外法向向量。