买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网获悉中国地质大学(武汉)申请的专利一种数值模拟器的开发方法、装置、系统及存储介质获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115481553B 。

龙图腾网通过国家知识产权局官网在2025-11-14发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211142428.0，技术领域涉及：G06F30/23；该发明授权一种数值模拟器的开发方法、装置、系统及存储介质是由李晓东;宁伏龙;雷刚;孙嘉鑫;程万设计研发完成，并于2022-09-20向国家知识产权局提交的专利申请。

本一种数值模拟器的开发方法、装置、系统及存储介质在说明书摘要公布了：本发明公开了一种数值模拟器的开发方法、装置、系统及存储介质，所述开发方法包括：采用网格剖分软件，针对目标研究区域进行网格剖分，并输出网格信息文件，转化为目标研究区域的节点坐标信息矩阵和节点编号信息矩阵；根据推导的增量解析公式，计算开采过程中天然气水合物的压力、温度、三相饱和度分布，同时采用伽辽金有限元方法计算区域内的位移分布；运用增量解析公式时围绕网格顶点构造控制体积；不同物理场之间的解耦将以显式递归的形式进行。该方法在用增量方案求解渗流传热系统时，无需求解线性方程组，避免了大型系数矩阵的生成，能节省一定的内存。该方案在考虑天然气水合物沉积物的弹塑性特点时将有良好的适应性。

本发明授权一种数值模拟器的开发方法、装置、系统及存储介质在权利要求书中公布了：1.一种数值模拟器的开发方法，用于水合物沉积物的多物理场耦合计算，其特征在于，包括： S100：采用网格剖分软件，针对目标研究区域进行网格剖分，并输出网格信息文件； S200：导入网格信息文件，并且转化为所述目标研究区域的每个节点坐标信息矩阵和每个网格所包含节点编号的信息矩阵以及边界区域的网格编号和节点编号信息矩阵； S300：根据水合物的温度场、渗流场、力学场和相变场的初始条件和边界条件，将温度、孔隙压力、位移、应力、水饱和度、气体饱和度、水合物饱和度和孔隙度进行初始化； S400：以网格顶点为中心构造控制体，计算所述控制体的各个边界的通量，并进行累加，以得到每个控制体的通量差； 计算所述控制体的各个边界的通量具体包括水通量和气体通量，其中： 水关于某一个截面的通量 1 气体关于某一个截面的通量 2 其中：是水相的渗透率，是气相的渗透率，和分别是根据有限元基函数和节点函数值组装而成的孔隙水压力和孔隙气压力，，，是有限元基函数，是网格顶点的数量，是该截面的单位法向量，代表当前截面的微分； S500：将每个控制体的通量差代入增量公式中,根据当前时刻的各场参数，计算下一时刻每一个网格顶点的温度、孔隙气压力、孔隙水压力、气相饱和度、水相饱和度和水合物饱和度的增量； 所述孔隙气压力增量的具体表达式如下： 3 所述孔隙水压力的增量的具体表达式如下： 4 所述气相饱和度增量的具体表达式如下： 5 所述水相饱和度增量的具体表达式如下： 6 所述水合物饱和度增量的具体表达式如下： 7 式中：，和， 其中：代表水通量的散度，代表气体通量的散度，代表气相饱和度，代表水相饱和度，代表水合物饱和度，代表有效水饱和度，代表气相的刚度，代表水相的刚度，代表水合物的刚度，是气相孔隙压力，是水相孔隙压力，是孔隙度，是体积应变，是时间步长，是气相的摩尔质量，是水相的摩尔质量，是水合物的摩尔质量，是水和数，是气相密度，是水相密度，是水合物密度，是气相热膨胀系数，是水相热膨胀系数，是天然气水合物热膨胀系数，是沉积物颗粒的热膨胀系数，是水合物分解的摩尔数，是温度，是毛细管压力，是孔隙气压力，是孔隙水压力，是气体饱和度，是水饱和度，是水合物饱和度的增量； 所述水通量的散度由式8计算得到： 8 所述气体通量的散度由式9计算得到： 9 式中：代表控制体积在每个周围网格的那一组成部分的编号，代表每个网格属于该控制体积的那部分区域，代表水的密度，代表气体密度，代表水相的渗流速度，代表气相的渗流速度，是当前控制体的面积，代表控制体，代表控制体的边界，是控制体边界的外法线单位向量； 所述计算下一时刻每一个网格顶点的温度的增量的计算表达式为： 10 其中：是温度的增量,是各相比热的加权平均值，是气相比热，是水相比热，是水合物比热，是热传导系数张量，是水合物分解导致的焓变，是有效应力，是塑性应变，代表热通量关于控制体积的散度； 600：针对含水合物沉积物多物理场耦合数学模型的力学控制方程，使用伽辽金有限元法求解纯弹性问题，使用增量变刚度法求解弹塑性问题； S700：根据计算得到的体积应变、温度和孔隙压力，更新储层的孔隙度，并根据所述孔隙度更新绝对渗透率； S800：判定计算结果是否满足稳定性条件，如果满足则继续循环计算，如果不满足，则降低时间步长，重新计算，直到满足稳定性条件为止，并以场图形式输出当前时刻的各物理场计算结果。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人中国地质大学(武汉)，其通讯地址为：430000 湖北省武汉市洪山区鲁磨路388号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。

以上内容由龙图腾AI智能生成。