天然气水合物多场耦合数值模拟方法研究取得新进展

近日，自然资源部中国地质调查局青岛海洋地质研究所在水合物多场耦合数值模拟方法研究方面取得了一项新进展。该研究提出了一种基于控制体有限元法和有限元混合的流-固耦合数值模拟方法，为水合物开采的产能和力学稳定性评价提供了一种新的数值模拟研究手段。成果发表在岩土工程领域顶级期刊Computers and Geotechnics(中科院一区Top，IF=4.956)。

天然气水合物的开采是一个涉及渗流、传热、水合物相变及沉积物力学变形的THMC多物理场耦合过程。从研究对象来看，其中渗流、传热和水合物相变可视作流体力学问题，而沉积物变形则是固体力学问题。由于研究对象的不同，流体和固体系统应当使用不同的数值方法来处理。目前，在水合物THMC耦合数值模拟方法有两种研究思路：第一种是使用有限体积法计算流体系统，而使用有限元法处理固体系统；第二种思路则是全部使用有限元方法同时处理流体系统和固体系统。然而，第一种方法在流固耦合过程中存在需要插值以及流体系统的网格不适用于固体系统的问题；第二种方法则存在着流体系统局部不守恒从而造成不收敛和非物理意义解的问题。

本研究提出了一种控制体有限元和有限元混合的新思路，在全非结构网格的基础上构建一套对偶的控制体网格，在控制体网格上用有限体积法结合有限元的插值函数进行流体系统的求解，而固体系统则在原网格上使用有限元法求解。这样一套新的研究思路有几个方面的优势：①构建的控制体使得流体系统和固体系统的求解变量位于相同的节点上，在保证流体系统求解局部守恒性的同时避免了流-固耦合的插值和网格不适应的问题；②引入有限元插值函数进行多点流量近似(Multiple-point flux approximation, MPFA)，避免了传统有限体积法对网格正交性的要求，可以使用完全非结构网格构建复杂的几何模型。基于上述算法，该团队开发了一套具有完全自主知识产权的数值模拟器QIMGHyd-THMC，利用室内模拟实验的结果和国际主流数值模拟软件对该模拟器进行了验证，为水合物多场耦合数值模拟方法提供了一种新的思路和框架。 

该研究得到了国家自然科学基金(41906187、41976205)、山东省自然科学基金(ZR2019BD058)、青岛海洋科学与技术国家实验室海洋矿产资源评价与探测技术功能实验室自主课题(MMRZZ201807)和泰山学者特聘专家计划(ts201712079)的共同资助。

开发的QIMGHyd-THMC模拟器与国际主流模拟器计算结果的对比