双水平集逼近油藏模型特征的数值模拟方法.pdfVIP

中国科学: 数学 2012 年 第42 卷 第5 期: 455 472 双水平集逼近油藏模型特征的数值模拟方法 ∗ 姚昌辉 贾尚晖 郑州大学数学系, 郑州 530001; 中央财经大学应用数学学院, 北京 100875 E-mail: chyao@lsec.cc.ac.cn, shjia@lsec.cc.ac.cn 收稿日期: 2011-12-22; 接受日期: 2012-03-07; * 通信作者 国家自然科学基金(批准号:和 资助项目 摘要 本文使用双水平集函数逼近油藏模型特征, 构造出 Uzawas 算法进行数值模拟. 对于两相流渗 透率的数值求解问题, 可以通过测量油井数据和地震波数据来实现. 将构造出来的带限制的最优化问 题使用变异的Lagrange 方法求解. 如果使用双水平集函数逼近渗透率函数, 则需要对Lagrange 函数 进行修正, 从而将带限制的最优化问题转化成无限制的最优化问题. 由于双水平集函数的优越性, 进一 步构造出最速梯度下降Uzawas 算法和算子分裂格式Uzawas 算法进行求解对应的最优化子问题. 数 值算例表明设计的算法是高效的、稳定的. 关键词 油藏模型 双水平集 历史匹配问题 渗透率 主题分类 49Q10, 35R30, 65J20, 74G75 引言 本文考虑带有各向同性渗透率的不可压、不混溶、水平方向的多孔介质油水二相模型 ∂S (κ()κ (S ) ) Φ() o − ∇ · ro o ∇p = f (), (1.1) ∂t µo o o Φ() ∂Sw − ∇ · (κ()κrw (Sw ) ∇pw ) = fw (), (1.2) ∂t µw 其中 (, t) ∈ Ω × [0, T], Ω ⊂ 2 是一个有界地质油藏, i 表示o 和w , 分别代表油相和水相, S , µ 和p i i i 分别为石油和水的密度、粘度和压力, κri 表示相对渗透率, κ() 表示绝对渗透率, Φ() 表示多孔介质 孔隙度. 假设油和水充满岩石的孔隙空间, 且 So + Sw = 1, (1.3) 并且假设毛细管压力函数Pc 已知, 且 p − p = P , (1.4) o w c 其中 Φ, κ, κri 和Pc 都依赖于多孔介质而不能直接测量. 假定Φ, κri 已知, 且设Pc 为0, 本文的目的是估计绝对渗透率κ(). 为了方便恢复渗透率, 本文使 用从油井和地震波中测量到的信息, 然而却不能直接得到渗透率的任何信息. 从流体方程 (1.1)–(1.4), 英文引用格式 Yao C H, Jia S H. Simulation of reservoir model characterization approximation of binary level set method (in Chinese). Sci Sin Math, 2012, 42(5): 455–472, doi: 10.1360/012011-1019 姚昌辉等: 双水平