裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

裂缝性油藏流固耦合渗流

本文给出了考虑介质变形的双重孔隙介质流固耦合渗流模型,并考虑渗流参数随有效应力而变化的非线性双重孔隙介质流固耦合渗流,在此基础上,本文还推导了双重孔隙介质非线性系数非线性等流固耦合流流计算,并给出了算例。     

裂缝性低渗透油藏流固耦合渗流分析

在低渗透油田的开发过程中，油藏流体渗流和储层岩土之间存在明显的耦合作用。本文首先研究给出了低渗裂缝性储层孔渗参数的等效方法，然后将渗流力学和岩土力学相结合，给出了低渗透裂缝性储层流固耦合渗流的数学模型，该模型不仅可以反映基质孔渗参数在开发中的变化，而且更能反映裂缝开度变化所引起的渗透率变化，而这对于低渗透裂缝性油田而言十分重要。最后对一实际井网进行了流固耦合油藏数值模拟，给出了开发过程中孔渗参数的变化及其耦合效应对油田开发的影响．     

基于XFEM-MBEM的嵌入式离散裂缝模型流固耦合数值模拟方法

离散缝网的表征与模拟是目前国内外研究的热点. 在非常规油气开发过程中, 由于地应力场的存在会对裂缝的流动属性产生显著影响, 若将裂缝视为静态对象, 与矿场数据会出现极大偏差, 因此要基于动态裂缝做更深入的研究. 本文针对致密油藏应力场?渗流场耦合力学问题, 提出了一种高效的混合数值离散化方法, 其中采用扩展有限元法 (XFEM) 求解岩石的弹性形变, 采用了混合边界元法 (MBEM) 精确计算基岩与裂缝间的非稳态窜流, 这两种数值格式是完全耦合的, 并对整体计算格式的时间项进行了全隐式求解, 可准确表征致密油藏开采过程中的裂缝变形及流体流动机理. 此外, 本文采用了嵌入式离散裂缝前处理算法显式表征大尺度水力压裂缝, 并考虑了支撑剂的作用; 采用了双孔有效应力原理和双重介质隐式裂缝表征方法, 可捕捉基质与小尺度天然裂缝的动态信息; 由此, 本文所提出的混合模型综合表征了基质?天然裂缝?水力压裂缝共同组成的致密油藏复杂渗流环境, 并通过几个实例论证了模型的准确性, 研究表明: 对致密油藏压裂水平井进行产能评价时, 应力场所引起渗流参数的改变及裂缝开度降低的影响不可忽略. 本文研究可为非常规油气资源的开发提供理论指导.      

薄膜结构流固耦合的CFD数值模拟研究

基于弱耦合分区求解策略,在CompaqVisualFortran6．5环境下搭建了薄膜结构三维流固耦合效应的CFD数值模拟平台。程序采用模块化编程思想,主要包含几何建模、流体分析、结构分析和数据交换四个模块。其中几何建模模块采用自行编制的膜结构找形分析程序,流体分析模块采用经过二次开发的计算流体力学软件FLUENT6．0,结构分析模块采用自行编制的膜结构动力分析程序MDLFX;在数据交换模块中,编制了基于薄板样条法的插值计算程序,以实现流固交界面上不同区域网格间的数据传递问题,编制了基于代数法和迭代法的动网格变形程序,以实现流固耦合运算中的动网格更新。基于该软件平台,对单向柔性屋盖和鞍形膜结构屋盖进行了流固耦合数值模拟,验证了方法的有效性。     

人内耳前庭系统膜迷路流固耦合数值模拟

基于健康志愿者右侧颞骨完整的组织切片数据,将人内耳前庭系统膜迷路内内淋巴液描述成不可压缩牛顿流体,将位于膜迷路壶腹端的平衡感受器壶腹嵴视为大变形线弹性材料,数值模拟人在立位转身、回顾时膜迷路内复杂的流固耦合作用,建立了内耳前庭系统膜迷路的三维有限元生物力学模型,分别得到不同旋转速度激励内淋巴液的压强场和壶腹嵴的位移场. 研究结果有助于从生物力学的角度定量分析人在旋转时前庭膜迷路内内淋巴液和壶腹嵴相互作用的力学特性.      

基于连接元法的多尺度裂缝性油藏数值模拟

针对油藏不同尺度复杂几何特征描述和动态连通性识别等难题,近年来发展了一种基于非欧物理连通网络具有无网格特征的油藏数值模拟连接元方法.文章将连接元法推广到裂缝性油藏,从流体流动的角度,利用连接单元将油藏离散为物理连通网络.根据节点物性参数、影响域半径和加权最小二乘法给出了压力扩散项的广义差分近似.结合物质守恒方程计算节点控制体积、基质节点间传导率、裂缝节点间传导率以及基质节点与裂缝节点间传导率,从而构建渗流控制方程组的全隐式离散格式,求解压力、饱和度以及含水率等生产动态参数.引入图论深度优先搜索算法,基于每个时间步求解的节点间压力梯度,计算各时间步注入井的劈分系数,定量表征井节点间的流动关系和连通性.算例验证表明,相较基于网格体系的传统方法,该方法能够自由灵活地刻画包括裂缝复杂分布、不规则油藏边界在内的复杂油藏几何,在粗化模型情况下能够保留更丰富的流动拓扑结构,实现计算精度和计算效率的更优平衡,能更好满足实际大规模裂缝性油藏的生产动态模拟预测需求,同时为具有多尺度几何特征的裂缝性油藏及复杂边界油藏的数值模拟提供了新思路.     

带有柔性薄板三维方柱流固耦合的数值模拟研究

采用双向流固耦合方法,对带有柔性薄板三维方柱的流场变化特性进行了研究。通过对比单方柱,分析了带有柔性薄板三维方柱阻力系数、升力系数以及斯特劳哈尔数的变化规律。研究表明,在方柱尾流区域附加一柔性薄板可以使其阻力系数降低34.6%,同时其变化幅值大大减小;其升力系数的均方根减小84.8%,流场脉动大幅度减小;斯特劳哈尔数降低79.5%。研究结果表明,在三维方柱后设置柔性薄板可有效抑制涡脱落,从而改善三维方柱的尾流特性。     

基于子循环-预测校正的三维SRM点火瞬态流固耦合数值模拟

为了分析固体火箭发动机点火瞬态工作过程,通过采用基于子循环-预测校正的紧祸合算法,研究高效的耦合器耦合计算流体动力学代码及计算结构动力学代码,对某翼柱型装药固体火箭发动机点火瞬态进行数值模拟,结果表明:翼柱型发动机点火过程中,尾部药柱壁面最开始点燃,同时火焰快速向翼槽和头部传播,翼段是最容易发生破坏的地方;紧耦合算法能...     

基于LS-DYNA的高速破片水中运动特性流固耦合数值模拟

基于大型有限元分析软件ANSYS/LS-DYNA,建立三维长方体高速破片在水介质中运动的有限元动力分析模型,采用ALE方法对破片在水下运动过程进行流固耦合数值模拟,获得了破片的速度衰减曲线。研究了速度衰减规律、破片墩粗变形规律以及冲击波传播过程。得到高速破片的侵彻能力随速度的变化规律:当初速度大于910~1 115m/s时破片头部将产生显著变形,并大大影响其侵彻阻力;当破片速度较小时,水中侵彻距离随破片初速的增大而增大,当破片速度达到某临界值以后,侵彻距离将随初始速度的增大而逐渐减小。     

悬臂止推箔片轴承三维流-固-热耦合数值分析

工作温度是影响轴承性能关键因素之一,为了深入分析轴承气膜温度变化对箔片轴承承载特性的影响,本文中建立了考虑冷却气实际流动特性的悬臂止推箔片轴承三维流-固-热耦合模型,研究了轴承间隙、转速以及冷却气流速变化时悬臂箔片轴承气膜压力场和温度场等的变化规律.研究结果表明：随着轴承间隙减小以及轴承转速增大,气膜表面温度逐渐升高,其中气膜高压区温度上升更快,气膜温度与轴承间隙呈现线性变化趋势,气膜温度与转速呈现非线性增长的趋势,相比于轴承间隙,轴承转速升高对润滑气膜温度场影响更大,气膜温度过高可能会造成轴承热失效,因此需要深入研究工作温度对轴承性能的影响;随着冷却气流速增大,轴承承载力和气膜温度逐渐减小.冷却气流速增大到一定程度时,降温效果不再发生明显变化.     

近场水下爆炸瞬态强非线性流固耦合无网格数值模拟研究

近场水下爆炸涉及多相流体的掺杂耦合以及结构的大变形、损伤和断裂等瞬态强非线性现象, 传统的网格算法在模拟近场水下爆炸时面临结构网格畸变、多相界面捕捉精度不足等难题, 鉴于此, 本文建立了完全无网格的近场水下爆炸冲击波和气泡全物理过程瞬态强非线性流固耦合动力学模型. 流体采用基于黎曼求解器的光滑粒子流体动力学(SPH)方法求解, 结构采用重构核粒子法(RKPM)求解, 并基于法向通量边界条件实现流固耦合. 为提高SPH对流场间断的求解精度, 引入黎曼问题思想并结合MUSCL重构算法, 为解决流场粒子体积变化剧烈导致的精度下降问题, 应用了自适应粒子分割与合并方法. 为模拟水下爆炸对结构造成的损伤断裂, 基于退化实体几何表述, 采用Lemaitre损伤算法, 建立了RKPM壳结构断裂损伤模型. 依据所建立的SPH-RKPM流固耦合模型, 对近场水下爆炸冲击波传播、气泡脉动与射流以及结构毁伤进行了模拟, 将得到的冲击波载荷、气泡演化以及结构响应与实验值和其他数值解对比, 验证了当前建立的SPH-RKPM流固耦合模型的有效性和精度, 并给出了水下爆炸载荷特性及其对结构的流固耦合毁伤机制与规律, 旨在为近场水下爆炸载荷预报提供理论和基础性技术支撑, 为毁伤威力评估和舰船防护结构设计提供参考.      

弹塑性多孔介质流固耦合新理论:混合耦合理论

在全球气候变化和双碳政策的大背景下,多孔介质中固体的变形和流体的输运问题变得尤为重要。然而,在多孔介质中建立流固耦合模型仍面临的挑战之一是需要考虑跨越宏观尺度到纳米尺度的耦合作用。本文利用基于非平衡热力学的混合耦合理论,提出了一个弹塑性多孔介质流固耦合新模型,在同一个理论框架内研究了弹性变形、塑性变形和液体渗流之间跨尺度的耦合,考虑了耗散过程中的熵产,并利用Helmholtz自由能来连接宏观尺度上的力学变形和纳米尺度上的液体输运之间的相互作用。在应力-应变关系中采用了弹塑性刚度系数以反映塑性的影响。同时,经典的达西定律扩展为考虑固体的塑性变形。通过与文献中模型的比较,验证了该模型的有效性。最后,数值分析表明在多孔介质的流固耦合中塑性变形具有比较显著的影响。     

MPS-FEM数值分析带自由面的流固耦合问题

随着计算科学的发展,研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中,耦合的移动粒子半隐式方法 MPS(Moving Particle Semi-Implicit method)和有限单元法FEM(Finite Element method)为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上,对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外,在该求解器框架内,基于FEM法拓展了针对结构场进行求解的功能。首先,对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后,采用该自研MPS-FEM耦合求解器,数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比,验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。     

化学-力学耦合理论与数值方法

该文研究了化学场中的质量扩散与力学耦合问题,构造了化-力耦合情况下的力学本构关系与质量扩散的本构关系,并由这些本构关系和化学场、力学场的控制方程,得到化-力场耦合的有限元方程.通过数值算例,详细分析了由应力场引起的质量重分布和由化学场引起的结构变形.研究表明,力学与化学之间存在明显的相互作用,并采用有限元数值方法进行了分析.     

非线性流固耦合问题的ALE分步有限元数值方法

推导了一种基于ALE（Arbitrary Lagrangian-Eulerian,即任意的拉格朗日-欧拉）描述下的带有自由液面不可压液体与运动结构非线性耦合问题的一种有限元数值计算方法,在时间域上采用分步格式计算,这种方法与混合插值方法比较其优点是速度和压力可使用同阶线性插值,给数值计算带来很大方便,将所得到的方法应用到结构与TLD（调谐式（液体阻尼器）装置之间的耦合问题,数值计算结果验证了本文所推导的方法。     

线性与非线性流固耦合动力学数值方法的进展及应用

本文综述了线性与非线性流固耦合问题数值方法的进展及工程应用。讨论了四种数值分析方法：(1)混合有限元–子结构–子区域数值模型,以求解有限域线性流固耦合问题,如流体晃动,声腔–结构耦合,流体中的压力波,化工容器的地震响应,坝水耦合等;(2)混合有限元–边界元数值模型,以求解涉及无限域的线性流固耦合问题,如大型浮体承受飞机降落冲击,船舰的炮击回应等;(3)混合有限元–有限差分(体积)数值模型,以求解不涉及破浪和两相分离的非线性流固耦合问题;(4)混合有限元–光滑粒子数值模型,以求解涉及破浪和两相分离的非线性流固耦合问题。文中推荐分区迭代求解过程,以便应用现有的固体及流体求解器,于毎一时间步长分别求解固体及流体的方程,通过耦合迭代收敛,向前推进以达问题求解。文中选用的工程应用例子包含气–液–壳三相耦合,液化天然气船水晃动,人体步行冲击引起的声腔–建筑结构耦合,大型浮体承受飞机降落冲击的瞬态动力回应,涉及破浪和两相分离的气–翼耦合及结构于水上降落的冲击。数值分析结果与可用的实验或计算结果作了比较,以说明所述方法的精度及工程应用价值。文中列出了基于流固耦合的波能采积装置模型,以应用线性系统的共振及非线性系统的周期解原理,有效地采积波能。本文列出了231篇参考文献,以便读者进一步研讨所感兴趣方法。     

固-液两相流黑水管道冲蚀磨损的数值模拟研究

煤气化黑水处理系统管道由于其流体介质高含固体颗粒和腐蚀性介质,且工作在高温、高压差环境中,极易受到冲蚀磨损和腐蚀的耦合作用而失效,影响其服役寿命.采用计算流体力学(CFD)方法数值模拟研究了煤气化黑水处理系统固-液两相流管道的冲蚀磨损行为和机理,以及流体介质速度和固体颗粒粒径对管道冲蚀磨损的影响规律,并分析了盲通管和涡室结构对弯管冲蚀磨损行为的优化改善效果.研究结果显示,煤气化黑水处理系统管线的冲蚀高危区主要分布在弯管外拱和变径管等结构突变区域;管道冲蚀磨损行为与其内部流体的运动和颗粒冲击特性有关;管道的冲蚀率均随着流体速度的增加而加剧,而粒径对弯管和变径管冲蚀率的影响并非单调关系,这与颗粒受力作用有关;弯管优化分析显示,涡室结构可以降低弯管的最大冲蚀率,减缓弯管的冲蚀磨损.     

单裂隙流-固耦合渗流的试验研究

通过对较大尺寸的单裂隙岩体试块进行不同侧面加载的渗流试验,在实验室里开展了单裂隙流 固耦合渗流研究,模拟核废料贮藏库的围岩自由面的最危险部位的渗流量 应力耦合状态。分析了裂隙岩体渗流与应力的耦合机理,获得了几种典型情况下的试验数据,并拟合出不同应力条件下单裂隙岩体渗流量与应力间数学经验公式。从而说明并非任一方向的应力增加都能使渗流量减小,而是裂隙岩体的渗流量随着双向压应力的增加而减少,随着平行于裂隙面方向的单向压应力的增加而增加。缝隙开度虽然随着法向应力的增加而逐渐减小,但最终不可能完全闭合,所以,此时流量不可能为零。同时,在试验过程中还通过闭环控制来实现被加载面的均匀受力,这为大尺寸岩体试验提供了一种很好的加载方法。     

油藏数值模拟发展现状

油藏数值模拟是以多相流体在多孔介质中流动的理论为基础,目前已广泛应用于油藏工程.本文就流体流动方程、数值计算方法、软件的开发和应用等方面,阐述了所谓黑油模型发展现状,同时介绍了组分模型和热采模型等其他类型模型.还就国内如何开展油藏模拟工作提出了建议.