完整可压缩渗流驱动问题的有限元--差分方法.pdfVIP

完全可压缩渗流驱动问题的有限元一差分方法 完全可压缩渗流驱动问题的有限元一一差分方法 夏炳松 (山东大学数学与系统科学学院。济南) 中文摘要 多孔介质中渗流驱动问题的科学研究对于工程技术领域的某些方面如地下 石油开采有着十分重要的意义，我们可以用现代计算方法对流体模型进行 数值模拟，为选择井位、i圭水量、生产量等提供依据。本文建立丁二相完 全可压缩渗流驱动模型，构造了有限元一差分格式，f其中压力方程采用有 限元方法求解．饱和度方程采用沿特征线修正的差分方法求冁)．得到了最 优误差估计。 ，本文主要分为三节。 第一节：描述完全可压缩渗流驱动问题的数学模型。 m)害十V·“+w(c)f，即一gD甲c-印=叮． 矿赛球)害坤如·Vic坷(DVc)州咖印呤啪 初值条件为： r1．xEQ： c(x，0)=co(r)．p(r．0)=P。I 完全可压缩渗流驱动问题的有限元一一差分方法 边界条件为： zf．疗=DVc．n=0r∈aQ(Q的边界)，n：aQ上单位外法向 矢量 第二节：数值计算格式 对压力方程，其有限元离散格式基于如下弱形式 讲 p∈H1(Q) 己知P”∈N。，以及C”，求尸”‘∈N。，采如下格式： 一(gDVqC”)·VP”，V)=((，，V)，V∈N^ 对饱和度方程： vc可以写成沿特征方向 注意饱和度方程中输运项≯宴+P(c)u (V 导数妒-c3c，这里 fr 妒=(≯二+P：(c)Il，|二)” 嘉=p弋妒妄w(c)孙V)0t Cl 因此我们在饱和度方程的离敞格式中采用对输运项直接沿特i芷线方向导数 逼近，于是构造出下面的格式： 丸学一h(DVhC1)。 完全可压缩溘流驱动问题的有限元一差分方法 第三节：误差估计 定理设初边值问题(11～(I4)解适当光滑．且是Q-周期的；如果剖分参 数满足 At=o(hp)，h。=o(hP)，I≥l， 则有限元．特征线修正差分格式(22)．(23)的解与真解之间有如下误差估计 max妒一c”忆+妙一∥U。}M(Al+h。+^；) 0≤Hs，df 其中肘依赖c、-及c关于特征方向的二阶导数窘。、) 关键词：可压缩渗漏洋动问题：特J压线修正差分；有限元一差分方法；误 差估计 致谢：首先感谢在三年研究生学习过程中为我授业解惑的老师们一袁益让 教授、羊丹平教授、粱栋教授、王文洽教授等．他们在学业上对我的启迪 帮助很大，特别要感谢我的导师程爱杰老师．他在我的论文写作过程中遇 到的许多问题给予了详细、耐心的解答和指导，程老师对学问研究一丝不 苟，这一点特别值得我学习。另外，对崔明博士的热心帮助表示谢意。、 完全可压缩渗流驱动问履的有限元一差分方法 te举多孔介质中渗流驱动问愿的科学研究对于工程技术领域的某些方面 如地下石油开采有着十分重要的意义，我们可以用现代计算方法对流体模 型进行数值模拟。为选择井位、注水量、生产量等提供依据．本文建立了 二相完全可压缩渗流驱动模型，构造了有限元—差分格式，(其中压力方程 采用有限元方法求解，饱和度方程采用沿特征线修正的差分方法求解)，得 到了最优误差估计。 关■诃：可压缩渗流驱动问题；特征线修正差分；有限元一差分方法；误 差估计 第—节数学文曩 在多孔介质中，我们考虑一种可压缩流体被另一种可压缩流体(可混溶) 驱动：取平面ZR2上单位厚度有界区域Q．忽略重力项，压力P和流体Darcy 速度u满足Deu∞?定律： 、 “：一掣vp：一口(c}Vp