同位网格界面速度插值方法研究

讨论了同位网格下,离散的连续性方程、动量方程及标量方程中控制容积界面上速度的计算方法．分别采用动量插值技术和线性插值技术计算了动量方程和标量方程的离散系数中的界面速度,并将两种方法得到的计算结果进行了比较．指出当采用线性插值技术去计算离散方程系数中的界面速度时,离散系数中的质量残余必须等于0,这样才能保证数值解的准确性和计算的收敛速度．     

高温下钙蒙脱石膨胀特性的分子动力学模拟

高温下蒙脱石的膨胀特性在核废料深部封存、二氧化碳封存及页岩气开发等应用中有着重要影响,但相关机理尚不明确.本工作使用分子动力学模拟为技术手段计算5 MPa和298—500 K等条件下,1.40—4.00 nm晶面间距(d)的一系列饱和钙蒙脱石的膨胀压力.以模拟所得的数值结果为依据,基于水化效应、双电层效应和离子关联效应等模型推演膨胀压力随温度与d的变化规律,并与相应的实验数据进行对比.模拟结果表明,当d较小时,因为高温会弱化水化力的强度,钙蒙脱石膨胀压力震荡的幅度降低,同时水化力作用的d的范围减小.当d较大时,因为高温强化离子关联效应,膨胀压力降低,同时双电层力的作用的d的范围增加.在较高温度和较大d时,膨胀压力为收缩力,阻碍膨胀.这些膨胀压力的变化规律与前期钠蒙脱石体系的研究类似.然而,通过对比两种蒙脱石体系的模拟结果,发现两种体系存在显著的差异—钙蒙脱石比钠蒙脱石更难膨胀到较大的d.此模拟结果与前人实验观测的结果相符.我们进一步将此差异归于钙蒙脱石的离子关联效应要远大于钠蒙脱石.有别于分子模拟中对于离子关联效应的精确描述,连续化的Poisson-Boltzmann方程因为忽略了离子关联效应,从而无法表达出与两种体系模拟结果都相吻合的膨胀压力变化规律.     

高温下钙蒙脱石膨胀特性的分子动力学模拟

高温下蒙脱石的膨胀特性在核废料深部封存、二氧化碳封存及页岩气开发等应用中有着重要影响,但相关机理尚不明确.本工作使用分子动力学模拟为技术手段计算5 MPa和298—500 K等条件下,1.40—4.00 nm晶面间距(d)的一系列饱和钙蒙脱石的膨胀压力.以模拟所得的数值结果为依据,基于水化效应、双电层效应和离子关联效应等模型推演膨胀压力随温度与d的变化规律,并与相应的实验数据进行对比.模拟结果表明,当d较小时,因为高温会弱化水化力的强度,钙蒙脱石膨胀压力震荡的幅度降低,同时水化力作用的d的范围减小.当d较大时,因为高温强化离子关联效应,膨胀压力降低,同时双电层力的作用的d的范围增加.在较高温度和较大d时,膨胀压力为收缩力,阻碍膨胀.这些膨胀压力的变化规律与前期钠蒙脱石体系的研究类似.然而,通过对比两种蒙脱石体系的模拟结果,发现两种体系存在显著的差异—钙蒙脱石比钠蒙脱石更难膨胀到较大的d.此模拟结果与前人实验观测的结果相符.我们进一步将此差异归于钙蒙脱石的离子关联效应要远大于钠蒙脱石.有别于分子模拟中对于离子关联效应的精确描述,连续化的Poisson-Boltzmann方程因为忽略了离子关联效应,从而无法表达出与两种体系模拟结果都相吻合的膨胀压力变化规律.     

页岩凝析气藏相平衡的快速准确计算方法

对页岩油气藏中复杂流体的相平衡计算需要建立考虑毛细作用效应的先进的数值模型, 并设计出快速可靠的算法以应对实际工况中储层流体包含多达数十种组分的复杂情况. 本文将基于适合页岩油气藏常见组分的真实流体状态方程, 即Peng?Robinson状态方程构建具有热力学一致性的VT型孔观相平衡计算体系. 通过引入描述毛细压力做功的数学模型实现对页岩流体热力学性质更准确的刻画. 结合扩散界面模型建立动力学演化格式, 采用成熟的凸分裂方法求解摩尔数和体积分数的演变, 从而描述相平衡的动态过程. 在此基础上, 本文开发了一套具有自适应性的深度学习算法, 设计了独特的双网络结构以实现对不同流体中不同组分的广泛适用性. 该神经网络的输入和输出参数均在热力学分析的基础上选取关键的热力学性质参数, 并进行了全面的超参调试以确定最合适的网络架构和最后形成的预测模型的基本结构, 且通过多种深度学习技术解决了过拟合问题, 在显著加速了传统的基于迭代方法的闪蒸计算的同时保证了相平衡状态预测的准确性, 得到了较好的预测效果. 相分离判定自动整合在预测结果中, 且从最终预测结果可以显著地捕捉到毛细作用的影响. 这一套快速、准确、可靠地基于深度学习算法的页岩油气孔观相平衡计算体系可以为后续的多相流动模拟提供具有物理意义的相分布初场, 确定系统内各个阶段的相数, 并可以作为构建具有物理守恒性的多相数值模型的热力学基础.      

Calculation of cell face velocity of non-staggered grid system

In this paper,the cell face velocities in the discretization of the continuity equation,the momentum equation,and the scalar equation of a non-staggered grid system are calculated and discussed.Both the momentum interpolation and the linear interpolation are adopted to evaluate the coefficients in the discretized momentum and scalar equations.Their performances are compared.When the linear interpolation is used to calculate the coefficients,the mass residual term in the coefficients must be dropped to maintain the accuracy and convergence rate of the solution.     

双孔双渗介质中气体流动的POD模型研究

将特征正交分解(POD)方法引入到页岩气流动双孔双渗模型的数值计算,建立了POD投影模型。在建模过程中发现将气体在基质和缝网中的流动控制方程分别进行POD投影的方式,虽然在数学上正确且满足基质或缝网的质量守恒,但会人为地在基质与缝网间附加一个质量传递投影项,无法保证整个流动系统的质量守恒,从而影响了POD模型的健壮性和精度。针对此问题提出了满足基质-缝网总体质量守恒的POD建模方式,在保证高精度(平均相对偏差约0.8%)和快速计算(计算速度提高720倍)的同时,大大提高了POD模型的健壮性。