首页/文章/ 详情

石油、天然气、地下水等在多孔介质当中的渗流力学前沿研究

5月前浏览11339

文一:

 

基于嵌套局部网格细化的嵌入式离散裂缝模型的裂缝瞬态流动建模

摘要:

水力裂缝周围的瞬态流动是致密油藏研究的一个重要课题。在这项工作中,嵌入式离散裂缝模型(EDFM)与嵌套局部网格细化(LGR)相结合,以提高其模拟近裂缝瞬态流的准确性。开发了一种自动网格细化程序,用于在复杂裂缝周围迭代进行LGR。在每次迭代中,都会选择并细化与裂缝相交或靠近裂缝的单元。可以调整细化区域的大小和细化级别。通过这一过程,近断裂区域的细胞大小迅速减小,而总单元计数保持在相对较低的水平。还分析了正则笛卡尔网格的最优细化率。此外,EDFM是在具有嵌套LGR的网格中实现的,并且矩阵和裂缝之间的连接被适当地确定,以计算作为储层模拟器输入的传输系数。首先针对不同裂缝方向的对数LGR模型验证了所开发的方法,结果表明,可以可靠地模拟不同的流态。随后,用这种方法研究了非正交、非平行和非平面水力裂缝的流态,以说明复杂裂缝周围独特的瞬态流动模式。这项工作为复杂裂缝几何形状的储层进行瞬态分析提供了一种通用而方便的方法。

 

图:储层尺寸、井位和裂缝方向。

 

图:说明使用带嵌套局部网格细化的EDFM来提高求解精度的基本思想。(a) 原始矩阵网格和裂缝。(b) –(e)裂缝周围细化后的矩阵网格。从(b)到(e),细化级别从一级增加到四级。

 

图:显示裂缝周围自动细化步骤的流程图。

 

图:具有八条相交裂缝的嵌套LGR示例(如红线所示)。

 

图:具有嵌套LGR的三维笛卡尔网格中压裂段之间连接的两个示例。(a) 平行于x-y平面的裂缝位置。(b) 平行于z轴的两个裂缝的位置。(c) (a)中裂缝的裂缝段之间的连接(粉红色线)。(d) (b)中裂缝的裂缝段之间的连接(粉红色线)。

 

图:储层尺寸、井筒位置和裂缝位置。紫色的线代表裂缝。如图所示,裂缝周围的基质单元被细化。

 

图:水力裂缝周围不同时间的压降。

原文PDF文件

Modelling fracture transient flow using the Embeded Discrete Fracture Model with nested local grid refinement.pdf


文二:

 

水合物油藏双井采气动态数值模拟

摘要:

多井系统可以增加排水面积,从而提高天然气水合物(NGH)生产的生产力。从根本上讲,双井系统是理解多井系统和制定含水合物沉积物采气现场策略的关键和基础。然而,通过降压回收NGH的生产性能和多井生产策略尚未得到详细处理。在这里,我们开发了一个三维模型来描述低渗透水合物油藏的双井生产特征,并分析了井的部署策略。在相同条件下,这种双井系统的生产力比单井高1.5–2.5,这是由于扩大了排水面积。井距改变了渗流特征和相应的地球物理场,这是双井系统采气的关键变量。其范围可以根据生产率和水合物离解前沿的组合来确定。此外,由于双井系统提高产能的能力有限,采用储层改造技术和增产措施的多井系统对实现天然气的商业开发具有重要的实用价值。

 

图:GMGS3-W19现场示意图。

 

图:模拟中的对称三维模型。

 

图:数值结果与实验数据的比较: (a)产气; (b)压力演化。

 

图:天然气水合物生产15天、30天、90天和180天后的水合物饱和度分布。

 

图:天然气水合物生产15天、30天、90天和180天后的含气饱和度分布。

 

图:天然气水合物生产过程中水合物和含水饱和度随 AB 线的变化: (a)水合物饱和度; (b)天然气饱和度。

 

图:生产过程中的流线: (a)和(b)流场;(c)速度。

 

图:双井系统与单井系统的比较: (a)产能; (b)生产360天后水合物的分布。

 

图:单井和双井系统地层响应的比较。

 

图:天然气水合物生产多井系统的研制。

原文PDF文件:

Numerical simulation of gas extraction performance from hydrate reservoirs using double-well systems.pdf


文三:

 

含水合物沉积物物理特性建模:考虑发生模式的影响

摘要:

研究含水合物沉积物的物理性质对水合物赋存模式和饱和度的依赖性,对于天然气水合物资源的勘探至关重要。然而,由于当前实验技术和数值建模方法的局限性,无法准确控制水合物饱和度和分布模式,因此对这些影响的精确分析仍然具有挑战性。在这项工作中,为了解决这些问题,我们提出了一种新的混合建模方法,该方法集成了X射线CT成像技术、形态运算算法和四元结构生成集方法。然后,我们生成了75个含有孔隙漂浮、胶结和桥接水合物的样本,每个样本都具有预定义的饱和度水平,并全面研究了水合物分布模式和饱和度对HBS的孔隙半径、配位数、相关函数、渗透率、电导率和弹性模量的影响。研究结果表明,不同水合物类型的孔隙和喉部半径分布的不均匀性不同。孔隙漂浮水合物的增加导致孔喉半径、曲率和孔隙空间相关性下降最快,但它使平均配位数变大,而其他配位数则减少。此外,在水合物饱和度较低的情况下,胶结模式对渗透率和电导率的损害最小。当水合物饱和度大于约45%时,桥接模式对质量传输特性的影响最大。此外,孔隙漂浮和桥接模式分别导致弹性模量的最大和最小增加。水合物的出现模式对体积模量的影响大于对剪切模量的影响。

 

图:HBS中天然岩心规模的水合物(a),以及实验室观测中水合物赋存模式的概念示意图(b)

 

图:从原始CT图像中选择的子体积(a),选择的具有300个体素的3D RVE(b),以及孔隙率随子体积大小的变化(c)。

 

图:五组三维数字样品中具有不同水合物饱和度的水合物不同赋存模式的三维外部视图。每一行代表一个组。灰色、蓝色和黄色分别表示固相、孔隙空间和水合物。所有图像的分辨率和大小分别为18μm/体素和5400μm。

 

图:第一组3D数字样品的2D薄片中具有不同饱和度的水合物的不同赋存模式。

 

图:第一组数字样品的孔隙网络包含具有不同饱和度的水合物的不同赋存模式。

 

图:不同饱和度水合物赋存形态数字样品的配位数和孔隙度分布。

 

图:第一组数字样本在不同饱和度下的压力场。

原文PDF文件:

Modelling the physical properties of hydrate-bearing sediments-considering the effects of occurrence patterns.pdf


文四:

 

光滑粒子流体力学在非均匀多孔介质多相流模拟中的应用

摘要:

多孔介质中的多相流建模,特别是在石油和天然气行业中,涉及数值求解微分方程。求解这些方程必须考虑化学、物理、机械和地质等几种现象,而这种方法没有解析解。通过无法涵盖广泛参数和操作条件的实验室实验,可以对这些现象进行表征。在这项研究中,我们研究了光滑粒子流体动力学(SPH)方法在复杂和随机几何多孔介质中建模的实际应用。本研究的主要目的是调查在实验室研究的同时使用SPH作为补充方法的可能性。首先,对单相流动进行了表征,提取了绝对渗透率张量和弯曲度等基本宏观参数。然后,对多孔介质中的两相流模型进行了研究。这些研究表明,可以获得实际参数,如局部饱和度变化、流体驱替、流动滞后、穿透时间、平均流速、相对渗透率曲线、Buckley−Leverett理论、毛细管去饱和曲线和采收率。最后,对非均匀多孔介质中的多相流进行了敏感性分析。为此,研究了三个参数(多孔介质的几何形状、粘度比和润湿性)对通过多孔介质的多相流的影响。结果表明,润湿性是影响轻相回收率的主要参数。在下一个层面上,多孔介质的几何形状对生产率的影响最为显著,而最低的影响与粘度比有关。研究结果表明,SPH方法可以为油藏工程师预测多孔介质中的多相流动行为提供一个合适的视角。

 

图:格子玻尔兹曼数值模拟的结果。

 

图:(A)在横截面y=9处沿着多孔介质的压差分布,以及(B)在稳态条件下流过g=0.1125的类型1的多孔介质的流的压差颜色图。

 

图:(A) 通过(1)GT#1、(2)GT#2和(3)GT#3进行单相建模的速度剖面;(B) 通过(4)GT#1、(5)GT#2和(6)GT#3,不同g的单相建模中顶面和底面的交叉流速。

 

图:部分润湿条件下燃气轮机2号两相流色图

 

图:(A)不同湿润条件下两相流的相对渗透率与饱和度和(B)分数流动曲线,(1)部分非湿润,(2)中性湿润,和(3)部分湿润。

 

图:对于所有类型的几何形状和粘度比,两相流中的相对渗透率与饱和度,(A)在不同的润湿条件下,(1)强非润湿,(2)部分非湿润,(3)中性润湿,(4)部分润湿,和(5)强润湿,以及(B)在不同润湿条件下交叉点的相对位置(即Krw=Kro)。

原文PDF文件:

mohammadi-et-al-2023-application-of-smoothed-particle-hydrodynamics-for-modeling-of-multiphase-fluid-flow-in-non.pdf


文五:

 

页岩气吸附机理: 来自页岩吸附甲烷热力学和动力学研究的证据

摘要:

研究页岩气吸附机理对页岩气成藏和储量评价具有重要意义。为了从甲烷吸附热力学和动力学的角度研究页岩气的吸附机理,对四川盆地南部下志留统龙马溪页岩样品在40.6°C、60.6°C,75.6°C和95.6°C压力高达52MPa的条件下进行了高压甲烷吸附和吸附动力学实验。获得了甲烷的吸附等温线和动力学曲线,并进行了详细的分析。结果表明:(1)在0~52MPa条件下,甲烷在页岩上的绝对吸附等温线具有I型吸附等温线的特征。温度对甲烷的最大过量和绝对吸附有重要影响。在相同的温度下,甲烷在页岩上的绝对吸附量在更高的压力下增加得较慢,这表明甲烷吸附速率在更高压力下降低。(2) 甲烷在页岩上的平均等温吸附热为21.06kJ/mol,表明甲烷在页岩的主要吸附过程可能是物理吸附。吸附的同位热随着甲烷绝对吸附量的增加而增加,表明吸附热主要受吸附的甲烷分子之间相互作用的影响。(3) Bangham动力学模型可以用来描述甲烷在页岩上的动态吸附过程。较高的温度和压力导致Bangham吸附速率常数较低,这使得页岩更难吸附甲烷分子。这与甲烷在页岩上的绝对吸附等温线的热力学研究结论是一致的。

 

图:研究页岩样品的岩心照片(a)和高压吸附分析仪示意图(b)。

 

图:不同温度下还原部分的过量吸附量与压力之间的关系。

 

图:甲烷吸附的标准熵(ΔS °)与等吸附热(Qst)曲线。

原文PDF文件:

Mechanisms of shale gas adsorption-evidence from thermodynamics and kinetics study of methane adsorption on shale.pdf

   


来源:STEM与计算机方法
ACTSystem断裂多相流化学动网格多孔介质通用ADSUM理论自动驾驶数字孪生人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-05-12
最近编辑:5月前
江野
博士 等春风得意,等时间嘉许。
获赞 46粉丝 41文章 300课程 0
点赞
收藏
作者推荐

高分子材料(Polymer materials)的数值模拟前沿研究分享

文一: 聚合物复合材料断裂的梯度相间增强相场模拟方法摘要:与微米级夹杂物相比,纳米级填料颗粒在改善聚合物的弹性和断裂行为方面的优势已通过几项实验观察得到证实。然而,对这种非均质材料中裂纹扩展的精确建模涉及复杂裂纹拓扑结构的解析,同时考虑多个裂纹的聚结和分支。这种复杂性使得传统的尖锐裂纹建模方法,例如基于单位分割思想的方法,适用性有限,尤其是在3D中,因为这些方法涉及对演化裂纹表面的显式跟踪。因此,相场断裂方法已经成为传统尖锐裂纹模型的一种有吸引力的替代方法,尤其是当需要处理复杂的裂纹拓扑结构时。此外,标准的两阶段连续谱模型由于缺乏必要的长度尺度,无法捕捉到前面提到的较小即较强的尺寸效应。为了弥补这一缺陷,在有限应变超弹性的背景下,我们在先前的工作中引入了标准一阶连续体模型的基于梯度界面的增强。本贡献将梯度界面的思想扩展到相场断裂方法。在此,考虑了如实验所观察到的具有连续变化或分级材料性质的填料颗粒周围的界面区域。复合材料的断裂行为可以通过用于确定界面区域内的弹性和断裂特性的分级程度来控制。与标准相场断裂模型相比,分级相间参数的最佳组合产生了更硬的宏观响应,如在更硬相间的情况下所观察到的。通过广泛的数值实验,描述了所引入的技术对广泛的实验观察到的断裂行为进行建模的适当性,这取决于填料和基体相之间的粘附程度。 图:连续体的尖端裂纹示意图和弥散相场裂缝模型的示意图。 图:梯度相界概念。 图:建模示意图。 图:标准算例的断裂相场建模示意图。 图:在低相间断裂能的情况下,分级相间弹性对复合材料断裂响应的影响。文二: 模拟聚合物复合材料尺寸效应的增强计算均匀化技术摘要:几项实验研究证实,纳米夹杂物是比微米夹杂物更有效的增强聚合物的增强材料。标准的一阶计算均匀化方案由于缺乏必要的长度尺度来模拟这种尺度效应而存在不足,需要对标准方案进行改进。本文基于界面能量学的思想,对一个这样的扩展进行了深入的评估。系统的数值试验和分析表明了上述方法在复合材料力学行为建模中的局限性,即填充材料比基体更加刚性。提出了一种基于连续分级中间相思想的方法。通过具有代表性的数值算例对该方法进行了综合评价,表明该方法适用于不同填料-基体刚度组合的纳米复合材料的建模。 图:计算均匀化过程示意图,描绘了两个相关尺度下的BVP。 图:RVE 由基体材料(蓝色显示)和圆形夹杂物(红色显示)组成,具有(a)含能界面或(b)厚度为 t 的梯度界面(绿色显示)。 图:二氧化硅填料颗粒在环氧树脂基体中的相间性质随指数 n 值的变化。 图:IECH尺寸效应缩放:环氧树脂基体中橡胶颗粒的宏观Piola应力随填料颗粒尺寸的变化。 图:在单轴拉伸载荷作用下,橡胶颗粒在环氧树脂基体中的粒径和界面强度发生变化时,相应的柯西应力分量 σ11的节点值的云图。 图:在选定的颗粒半径下,RVE 在剪切载荷作用下的变形形状。为了显示填料和基体的相对位移,将变形形状叠加在未变形形状上。 图:相间厚度对容积加载尺寸效应标度的影响。 图:三维RVE的结构化六边形网格,a=100,r=25(和t=5),分别用于IECH和GICH。 图:三维强度尺度参数对两种方法的影响。本文研究了在体积载荷作用下,环氧树脂基体中二氧化硅颗粒在 a = 100,r = 25(和 t = 5)的固定微结构构型。文三: 扩展一个通用的快速粗粒度分子动力学模型来研究接枝聚合物纳米复合材料的力学行为摘要:聚合物纳米复合材料因其高通用性和良好的力学性能以及低密度而成为工程应用的一类重要材料。通过将聚合物链直接连接到纳米填料上,即所谓的接枝,实现了基质和填料之间更好的负载转移,此外,还获得了填料更好的分散性。两者都能增强力学性能。由于纳米尺度的实验研究极具挑战性,因此需要补充的数值研究来揭示聚合物纳米复合材料的力学行为。为此,分子动力学非常适合,因为它可以捕捉微观结构,但在数值上计算成本也很高。因此,这一贡献为研究接枝聚合物纳米复合材料的力学行为提供了一个快速粗粒分子动力学模型。为此,我们通过接枝键扩展了现有的模型,这使我们能够直接比较未处理和接枝填料的效果。特别是,我们研究了填料含量、接枝度和填料尺寸对聚合物(接枝)纳米复合材料刚度和强度的影响。我们得出的结论是,与未处理的填料相比,接枝键对刚度的影响很小,而强度显著提高,这与文献一致。所提出的聚合物接枝纳米复合材料的分子动力学模型为进一步研究,特别是关键的基质-填料界面提供了基础。此外,这一贡献将分子动力学的见解转化为力学性能,从而弥合了与工程规模的差距,从而代表着朝着开发聚合物(接枝)纳米复合材料的全部潜力迈出了一步。 图:嵌入热塑性基体的无机填料: a)未经处理填料的聚合物纳米复合材料(PNC) ; b)接枝填料的聚合物接枝纳米复合材料(PGN) ,即聚合物链与填料表面之间的共价键。 图:单轴拉伸模拟的刚度和强度。 图:填料含量的影响。 图:接枝键在聚合物链中的位置。文四: 非线性光弹性的数值模拟摘要:当分子光开关,如偶氮苯或降冰片二烯,嵌入到足够柔软的聚合物基体中时,所得化合物可以经历由特定波长的光引起的机械变形。这些光敏化合物具有作为软激励器应用的潜力,而不需要硬接线电子器件或单独的能源。这种特性在微型机器人的设计中 特别有吸引力,但其他应用,如高速数据传输或将光子能量转换为机械响应,也有很大的前景。尽管有这些几乎未来主义的可能性,光敏聚合物在工业应用中还没有得到足够的关注。提高这类软智能材料接受度的一个重要因素是结合数值模拟方法制定严格的本构建模方法。因此,在这篇文章中,我们提出了一种光机械建模方法,与之前发表的基本原理不同。我们简要介绍了必要的本构方程,这些方程随后与各自的平衡定律结合使用到有限元实现中。最后,通过一个简单的二维基准点示例说明了数值求解方法的能力,并随后将其扩展到更复杂的三维问题。 图:光致变形下未变形体和变形体的示意图。 图:(左)在可变强度光源照明下的材质样本草图。(右)施加10 mV电势差时有限元模型的变形。颜色映射和箭头指的是均匀化电子序参数的大小和方向。 图:(左) X 方向(虚线) ,Y 方向(虚线)和不同穿透深度的大小(实线)材料样品顶边中心点的位移。(右)变形材料样本。颜色是指在左边的图中选择的颜色。 图:(左)花瓣的一半与140个六面体单元相啮合的模型。(中)闭合的花,由八个完整的花瓣组成,以2 V的施加电势排列成花朵。(右)开放的花,以0.2 V的施加电位开放。彩色图是指位移大小。 图:(左)对于各种材质厚度值,花瓣尖端在Y方向(虚线)、Z方向(短划线)和大小(实线)上的位移。(右)变形的花瓣。着色是指变形。文五: 柔性管状金属-聚合物粘接接头的数值研究摘要:粘结管状接头已在各种工程和工业应用中得到应用,为连接管状基材提供了一种高效耐用的方法,同时确保了良好的密封性能。在这项工作中,提出了一种新的设计方法,用于在复杂载荷和被粘物大变形下承受大扭转旋转的粘结柔性管状金属-聚合物接头。为此,开发了两个不同的有限元模型。一种是针对较小的旋转值优化接头性能,另一种是使用内聚区建模并引入聚合物管状粘附物中的屈曲效应来预测整体接头强度。为了评估超弹性特性对接头性能和强度的影响,还考虑了不同的聚合物材料。三种不同接头布局的实验结果与数值结果吻合良好。结果表明,对于较小的旋转值,可以使用二维模型来准确预测关节行为。尽管如此,对于聚合物材料的高旋转值和大变形,复杂的3D模型(包括屈曲的影响)对于预测接头强度是必要的。 图:纯I型和II型的线性软化牵引分离定律。 图:管接头的布局、尺寸和胶粘剂厚度(以毫米为单位) : (a)-布局 X,(b)-布局 Y 和(c)-布局 Z。 图:用于制造管状接头的模具。 图:扭转试验装置。 图:不同关节布局的扭矩-扭转曲线实验研究。 图:通过解析解和本模型得到粘结层中的剪应力分布,采用粗网格和细网格。 图:变形形状,最大旋转角度沿轴线 = 180 ° 。 图:三维模型用于预测接头的最终旋转,网格表示。 图:对于布局 X、布局 Y 和布局 Z,前三个折叠模式形状。 图:不同布局的数值和实验扭矩-扭转曲线。 图:在旋转800 ° 时 X 排样的变形形状和粘接层中的 SDEG 损伤变量。 来源:STEM与计算机方法

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈