【技术实现步骤摘要】
本专利技术涉及一种基于孔隙网络模型的数值重构毛细管压力计算方法，属于环境岩土工程中的毛细管压力计算。
技术介绍
1、目前，我国能源、化工等行业污染尤为严重，产生大量有机污染物，集中分布在京津冀、长江经济带、粤港澳大湾区等沿海、沿江重点区块。此类有机污染物因在水中溶解度非常小被统称为非水相流体(nonaqueous phase liquids，napls)，根据其密度与水密度的对比，又分为轻质非水相污染物(lnapls)和重质非水相污染物(dnapls)。这些dnapls包括但不限于氯化乙烯、四氯化碳、三氯乙烯等有机溶剂，它们具有挥发性和毒性，对地下水系统和生态环境造成潜在威胁。它们的运移行为复杂多变，受到多孔介质的物理和化学特性以及地下水流场的影响。准确描述和预测dnapls在地下水和土壤中的行为对于环境保护和污染治理至关重要。
2、毛细压力-饱和度-相对渗透率本构模型(k-s-p曲线)在研究dnapls在多孔介质中的迁移规律方面起着至关重要的作用。k-s-p曲线能够考虑多相流动的毛细压力作用、不同相之间相对渗透率变化，从而更准确地模拟和预测dnapls的迁移行为，对于制定有效的环境污染治理策略至关重要，以减少或阻止dnapls的进一步扩散和影响。
3、传统试验手段获取毛细压力曲线通常涉及繁琐的实验流程和较长的时间成本。此外，这些试验还可能受到实验条件的限制，难以完全模拟实际场景。另外，数值模拟方法虽然也可以用来估计毛细压力曲线，但常常受到土体微观结构的复杂性导致模拟计算非常复杂需要大量计算资源。
...【技术保护点】
1.一种基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤1中，采用PFC离散元软件构建土体微观颗粒堆积结构三维图形，并将三维图形导出为VRML格式文件。
3.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤2中，利用AVIZO图像处理软件将三维图形转换成若干二维切片图片，并将二维切片图片保存为RAW格式文件。
4.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤3中，球棒孔隙网络模型的结构信息包括孔隙、喉道、截面形状因子、孔隙半径、喉道半径、喉道长度、孔隙与喉道的位置、孔隙与喉道的连接关系。
5.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤4中，使用准静态侵入渗流算法对流动前方的孔喉进行排序来识别流动阻力最小的路径，然后非湿润相侵入阻力最小的喉道，即毛细管压力最小的喉道，继续搜索下一次驱替过程中毛细管压
6.根据权利要求5所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述毛细管压力由Young-Laplace方程计算得到，具体如下：
7.一种计算机设备，包括存储器、处理器，以及存储在所述存储器中并能够在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至6任一项所述的基于孔隙网络模型的数值重构毛细管压力计算方法的步骤。
8.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述的基于孔隙网络模型的数值重构毛细管压力计算方法的步骤。
...【技术特征摘要】
1.一种基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤1中，采用pfc离散元软件构建土体微观颗粒堆积结构三维图形，并将三维图形导出为vrml格式文件。
3.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤2中，利用avizo图像处理软件将三维图形转换成若干二维切片图片，并将二维切片图片保存为raw格式文件。
4.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤3中，球棒孔隙网络模型的结构信息包括孔隙、喉道、截面形状因子、孔隙半径、喉道半径、喉道长度、孔隙与喉道的位置、孔隙与喉道的连接关系。
5.根据权利要求1所述的基于孔隙网络模型的数值重构毛细管压力计算方法，其特征在于，所述步骤4中，使用准静...
【专利技术属性】
技术研发人员：刘昂，何家兴，林文丽，丁志涛，张尔康，
申请(专利权)人：南京工业大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。