多孔介质热流固耦合有限元分析.docxVIP

PAGE 1 PAGE 1 多孔介质热流固耦合有限元分析  探讨了应用有限元分析软件ADINA进行多孔介质热流固耦合的分析方法。以石油生产中的热力采油为例，对多孔介质热流固耦合问题进行数值模拟。分析了在热力采油过程中由于注蒸汽以及流体流淌等因素对储油层的影响。对开采选址与提高石油产量有一定的参考意义。  0引言  在热流固耦合问题研究领域，最早是对地热资源的开发研究，探讨地热区域内地应力、温度变化和岩石渗透率变化的影响；近年来，国内学者研究了非等温条件下煤层中瓦斯渗流的规律，并建立了数学模型，考虑了温度变化对瓦斯渗流的影响。在多孔介质热流固耦合研究方面，特殊是在石油工程的研究过程中。Lewis等进行了非等温条件下油藏渗流的规律以及因石油开采而引起的地面下沉等问题的研究，考虑了温度场的变化对油藏渗流和地面沉降的影响以及流体渗流对温度场的影响等；但是，这些研究都没有考虑温度场的变化对多孔介质变形的影响(即没有考虑温度与固体之间、温度与流体之间的耦合作用)，因此没有实现真正意义上的热-流-固三场耦合。  由于热流固耦合理论开展研究的时间比较晚以及自身的复杂性等原因造成了其还处于初期的探索阶段。所以现阶段对热流固耦合的研究主要是借鉴流固耦合理论的思想和方程，例如双向相互作用的思想对于研究热流固耦合有很大启发，流固耦合理论中的岩石变形方程可以通过转变本构关系推广到热流耦合理论中的岩石变形方程等。  1热流固耦合的两种模型  对于热流固耦合理论目前还处于探索阶段，所以对该理论的深入研究主要以流固耦合理论为平台，借鉴其中的一些思想和方程。现在对该问题的研究主要有两个方向，并提出了两类模型：  第一类称为非完全耦合模型，主要是以流固耦合理论为基础进行研究。在此模型中，认为骨架温度的变化会引起岩体的体积热应变，热应变与其他因素引起的岩体应变共同组成了固体骨架的总应变，骨架的总应变引起了岩体的变形和流体渗流，从而成为影响固体变形方程和流体渗流方程的重要因素。将这种模型称为非完全耦合模型是因为在模型中将温度当作假想的与其他力学载荷等效的载荷(即所谓的“热载荷”)加载到岩体的固相骨架上，这种“热载荷”只是温度函数，并没有体现多孔介质变形、流体渗流的耦合项，也就是说在该模型中只有流体和固体两者之间的完全耦合，而热应变方程只是温度的函数，并没有实现完全耦合。此目的是为了简化热流固耦合模型，且使模型求解简洁，便于解决实际问题。  第二类称为完全耦合模型，该模型认为温度场的变化会对多孔介质的变形以及流体的渗流造成影响，并且还认为二者反过来也会对温度场造成影响。温度场和多孔介质变形、流体渗流之间的关系与流体渗流和多孔介质变形之间的双向相互作用类似，温度场的变化与岩石变形之间也存在双向相互作用，即热一固耦合效应，一方面温度场的变化要产生热应力、热应变从而导致多孔介质的变形，而另一方面，多孔介质的变形反过来又要导致温度场的转变。这样，流体渗流、固体变形以及温度场的变化三者同为基本变量，从而形成热流固三者的完全耦合。  2多孔介质热流固耦合模型建立  2．1建模的基本假设  多孔介质热流固耦合模型的建立基于以下基本假设：  1)多孔介质材料的性质和压力是匀称各向同性的，并且在传热传质过程中不发生变化。即孔隙率和肯定渗透率等参数维持不变；  2)比热、导热系数(固、液、气)均为常数；  3)多孔介质内流体流淌的速度小，其动能和惯性的影响可不予考虑；  4)固体相是一个固定的、不变形的骨架；  5)流体始终布满多孔介质的孔隙结构。  2．2建模的基本步骤  应用ADINA对多孔介质热流固耦合模型进行计算分析主要有：建立模型、加载求解、分析计算结果三部分组成。详细的步骤如图1所示： 图1ADlNA分析流程图 2．3建立模型  模型模拟的是石油工程热力采油过程，几何模型为6m×4m×6m，分三个地层，中间层为储油层，为多孔介质属性。设有三口井，其中一口作为热力采油主人蒸汽的人口，其余两口井为出油口，同时中间层的两侧也有流体压力存在，由于压差的存在，将导致中间层发生变形。模型采用流固耦合算法，有结构和流体两个模型，温度作为“热荷载”施加到模型上。  2．3．1建立结构模型  在进行流固耦合计算时，不管是流体模型还是结构模型都必需施加流固耦合边界，模型中有三个流固耦合边界条件，分别为中间地层的两侧和井的人口。模型中，因为要考虑温度的影响，即施加温度荷载，所以就必需指定模型在施加温度荷载之前的初始荷载。并且在划分网格之前需要定义网格密度，由于模型中各边的尺寸相差很大，所以应采用不同的方法进