基于N-S方程的储层渗流与孔眼管流耦合分析.docVIP

基于N-S方程的储层渗流与孔眼管流 (辽宁工程技术大学辽宁阜新123000）摘 要：用多孔介质渗流力学、流固耦合理论，建立了应力作用下油藏流固耦合数学模型，依据流体动力学管流理论用Navier-Stokes方程控制井筒和孔眼中流体流动；采用多重物理量耦合分析软件COMSOL Multiphysic对所建的模型进行数值模拟分析，实现了对四个偏微分方程组的全耦合求解，对模拟结果研究表明孔眼末端压降和应变最大，该处岩石容易失稳模拟结果证明所建模型和求解方法是有效的。 关键词：，耦合，数值模拟 中图分类号：Analysis of Reservoir Seepage and Perforation Pipeline Coupling Under Stress Based on N-S Equation LIU Beizhan LIANG Bin (Science College of Liaoning Technical University, Fuxin 123000） Abstract：eory. On the basis of the hydrokinetics pipe flow theory, the flow in well and well bore is controlled by the Navier-Stokes equation. The established model is analyzed numerically with COMSOL Multiphysic, a multi-physical quantity coupling analyzing software. And four systems of partial differential equations are resolved with unity coupling. The study on simulation results shows that the pressure drop and strain are the biggest at the perforation. The rocks at the spot are easy to collapse. The simulation result proves that the established model and solution are effective. Key Words: N-S equation, seepage, coupling, numerical simulation 1 引 言 完井的数值模拟研究是国内外油田较常使用的一种方法, 长期以来不少学者在这方面做了大量的工作, 采用的方法主要用有限元或有限差分法，取得了一些成果，这无疑对油田优化设计提供了理论依据[1,2]。但射孔完井是一个流、固、多相流多场耦合过程，求解多物理耦合问题通常首选全耦合算法，因为在理论上它能给出最真实的结果。然而由于同时求解多物理耦合过程异常复杂，到目前为止还没有这方面令人信服的研究成果[3]。因此各学者不得不对模型进行简化和假设，使得这些研究不能很好地反映流体汇流入井眼和井筒内流体的流动关系[1]。基于此，本文以COMSOL Multiphysics多场物理量耦合软件为平台对应力作用下油气藏渗流、孔眼和井筒流体流动进行数值模拟。由于油气藏的渗流过程是在流固耦合作用下进行的，因此建立了应力场作用下的渗流控制方程[4]；孔眼和井筒为大孔道的流体流动，流动特性与流体动力学中管流相似，因此本文用Navier-Stokes方程来模拟孔眼和井筒中流体流动过程；并对建立的应力场方程、渗流场方程和流体流动方程这一偏微分方程组用COMSOL Multiphysics软件成功地进行了求解，从而实现同时求解多物理耦合过程，而且用数值模拟结果证明了本文建立的数学模型和求解方法是有效的。 2 基本理论 油气在储层中运移和产出是一个多场耦合的复杂过程，油储层是一种多孔介质，石油在储层中的运移属于渗流问题，常用arcy运动方程来描述，而石油在孔眼和井筒中的流动属于管流，用流体动力学中的Navier-Stokes方程来控制其运动，同时流体的运动受储层应力的影响，因此，油气的运移和产出是一个复杂的流-固耦合问题。在建立方程前首先引入几个假设：①油储层是各向同性的线弹性材料；②不考虑温度的影响；③储层固体骨架变形遵从Terzaghi有效应力原理： （1） 其中，为有效应力，为总应力，为孔隙流体压力，为Biot常数，为Kroneker张量。 2.1 油藏岩体应力场方程 应力平衡方程： （2） 几何方程：