弹塑性变形油藏中多相渗流的数值模拟.pdfVIP

第 16 卷第 1 期 计 算 力 学 学 报 V o l. 16 N o. 1 1999 年 2 月 CH IN E SE JOU RN AL O F COM PU TA T ION AL M ECHAN IC S F eb ru ary 1999 弹塑性变形油藏中多相渗流的数值模拟 　　　冉启全 顾小芸 ( ) ( ) 中国石油天然气集团公司石油勘探开发科学研究院, 北京, 100083 　　 中国科学院力学研究所, 北京, 100080 摘　要　基于流固耦合力学理论, 建立了弹塑性变形油藏中多相渗流的数学模型, 该模 型考虑了渗流与变形的耦合作用, 以及注采交变载荷作用下油藏多孔介质的弹塑性变 形特征; 给出了耦合数值模拟方法和算例。 关键词　弹塑性变形; 多相渗流; 耦合作用; 数值模拟; 油藏 分类号　T E 3 19; O 357 3 1　前　　言 ( 计算渗流力学在石油工程中得到了相当重视和广泛应用。油藏渗流的计算机数值计算 即 ) 油藏数值模拟 已成为油藏工程师进行油藏工程分析必不可少的有力工具。但长期以来, 在研 究油藏开采过程中的渗流问题时, 均采用纯渗流力学理论进行研究, 不考虑油藏渗流与岩土变 形的耦合作用及其影响。近年来, 随着流固耦合力学理论和数值计算技术的发展, 人们对油藏 开采过程中的力学机制有了更深入的了解和认识: 岩土在多相流体渗流作用下会发生变形, 而 岩土变形反过来又影响多相流体的渗流。因此, 油藏开采过程中的多相渗流问题是流固耦合问 题, 而应用计算机对油藏开采过程中的流固耦合渗流问题进行数值计算和模拟, 称为流固耦合 油藏数值模拟。 国外在流固耦合油藏数值模拟方面已作过研究工作, 不少学者作了有益的探索[ 1～ 3 ] 。但这 [4, 5 ] 些研究均将岩土介质视为弹性介质, 对岩土受力变形过程作简化处理。大量的实验研究 表 明, 高温高压下的油藏岩土, 特别是在注采交变载荷作用下其受力变形特性极为复杂, 不仅有 可恢复的弹性变形, 更多的表现为不可恢复的塑性变形。为此, 本文将渗流力学与岩土力学相 结合, 充分考虑注采交变载荷作用下油藏岩土的弹塑性变形特性, 建立了体现多相渗流与弹塑 性变形的流固耦合模型, 用有限差分与有限元相结合的耦合数值计算方法对模型进行求解。 2　基本理论 流固耦合问题由流固耦合方程来描述。这组方程包括流固耦合渗流方程和流固耦合变形 方程, 而未知量含有描述渗流现象的变量和描述变形现象的变量, 流固耦合效应在两类方程中 均要得到体现。 2 1　渗流方程 油藏可以看成是势量所形成的场。当油藏内各处势量存在差异时就会发生渗流。首先引 入流动势的概念, 定义下式: = P - D ( 1) 收稿 日期: 　冉启全: 男, 1965 年生, 博士后, 高级工程师 1 期 　　　　　　　　 冉启全等: 弹塑性变形油藏中多相渗流的数值模拟 　　　　　　　25 式中= o 和 分别代表油相和水相, 为流体相的流动势, P 为流体相的