一种温度场耦合驱动下土壤-地下水有机污染运移数值模拟方法及系统.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN115983143A(43)申请公布日2023.04.18(21)申请号202211045827.5G06F119/08(2020.01)(22)申请日2022.08.30G06F119/14(2020.01)(71)申请人河海大学地址210000江苏省南京市江宁区佛城西路8号申请人浙江省环境科技有限公司(72)发明人杨蕴崔孜铭张弛张蔚王锦国窦智陈舟庄超(74)专利代理机构南京经纬专利商标代理有限公司32200专利代理师罗运红(51)Int.Cl.G06F30/28(2020.01)G06F111/10(2020.01)G06F113/08(2020.01)权利要求书3页说明书7页附图5页(54)发明名称一种温度场耦合驱动下土壤-地下水有机污染运移数值模拟方法及系统(57)摘要本发明公开了一种温度场耦合驱动下土壤‑地下水有机污染运移数值模拟方法及系统，用于模拟温度场驱动下有机污染物的运移规律。该模型能够考虑温度改变有机污染物的理化性质及其驱动参数，进而影响其在土水介质中的迁移及空间分布，构建了耦合非等温多相流、溶质运移和物理化学驱动过程的模型来模拟温度场影响场地有机污染物迁移的过程。本发明所提供的方法包括：考虑温度改变有机污染物的理化性质及其驱动参数，构建耦合非等温多相流、溶质运移和物理化学驱动过程的模型，开展温度场驱动下有机污染物的迁移分布规律数值模拟。本发明所提供的模拟技术可为场地有机污染防治提供科学合理的理论依据。CN115983143ACN115983143A权利要求书1/3页1.一种温度场耦合驱动下土壤‑地下水有机污染运移数值模拟方法，其特征在于，该方法包括如下步骤：步骤SS1：确定温度场耦合驱动下有机污染物运移模型的基本参数，所述模型基本参数包括：地质结构模型的地层密度、孔隙度、渗透率，相对渗透率，毛细压力，非等温多相流模型的大气压强、水力梯度和温度梯度，非等温多相流溶质运移模型的有机污染物密度、粘度、溶解度、扩散率，以及有机物泄露位置及泄露速率；步骤SS2：根据步骤SS1收集的模型参数设置地层的岩性参数以建立地质结构模型，并进行网格剖分；步骤SS3：在步骤SS2中所建地质结构模型的基础上，根据质量和能量守恒的基本方程、多相流控制方程、温度场方程建立非等温多相流模型，步骤SS4：确定非等温多相流模型的初始条件、边界条件、进行求解，以计算模拟区域的渗流速度和温度的时间和空间分布；步骤SS5：在步骤SS4获得的非等温多相流模型结果的基础上，结合温度与理化性质及驱动参数耦合数学方程建立非等温多相流溶质运移模型；步骤SS6：设置非等温多相流溶质运移模型中有机污染物的理化性质、扩散率、以及泄露速率参数，并确定非等温多相流溶质运移模型的初始条件，源汇项；步骤SS7：求解非等温多相流溶质运移模型以得到模拟区内时间上和空间上有机物各相饱和度的分布，刻画出有机污染物泄露后的迁移规律，完成温度场耦合驱动下有机污染物土壤‑地下水系统中运移的数值模拟。2.根据权利要求1所述的温度场耦合驱动下土壤‑地下水有机污染运移数值模拟方法及系统，其特征在于，所述步骤SS3中的非等温多相流模型构建过程如下：(1)采用渗流场方程、质量和能量守恒方程、多相流控制方程计算模型中渗流速度的空间和时间分布、各组分质量、组分各相饱和度：a.渗流场方程为：式中，ρ为流体密度；vβx,vβy,vβz为渗流速度v在三个坐标轴上的分量；Δx、Δy、Δz为单元体体积，t为时间，n为孔隙度，S为相饱和度；b.质量和能量守恒方程：k式中：m为流动区单元体表面的外向单位矢量；M为组分κ在单位土壤介质中质量，Vm为k流动单元体的体积；Γm为流动单元体的表面积；F为进入到流体单元体的组分κ的总通量；qk为组分κ在单元体的源汇相；c.多相流控制方程：式中，β为指示相，为组分κ的β相通量，k为总渗透率，Kγβ为β相的相对渗透率，处于02CN115983143A权利要求书2/3页到1之间，ρβ为β相的密度，μβ是β相的动力黏度，Pβ是β相中的流体压力，g重力矢量；(2)使用温度场方程计算得到在考虑热量交换以及各种反应引起的热量变化导致温度变化下模型在x、y、z方向上的温度梯度分布情况：3式中：Cm、Cw分别为多孔介质和水的热容量，J/(m·℃)；T为模拟区的温度，℃；λxx、λyy、λzz为x、y、z方向的热动力弥散系数，J/(m·d·℃)，vx、vy、vz为x、y、z方向的渗流速度，f1为单位时间单位体积含水层内由其它的化学反应或者微生物活动引起的热量变化，J/(t·m3)。3.根据权利要求2所述的温度场耦合驱动下土壤‑地下水有机污染运移数值模拟方法及系统，其特征在于，所述步骤SS5中，结合温度与有机物理化性质及驱动参数耦合