油藏数值模拟实验报告记录

1.油藏数值模拟:就是应用计算机研究油气藏中多相流体渗流规律得数值计算方法,它能够解决油气藏开发过程中难以解析求解得极为复杂得渗流及工程问题,就是评价与优化油气藏开发方案得有力工具、2.油田开发得任务:就是从油田得客观实际出发,以最少得投资、采用最合理得速度,获得最高得最终采收率,也就就是获得最大得经济效益,油藏数值模拟就是达到以上目得较有效得方式。

3.油藏描述:油田开发后对油藏动态得认识,特别强调对剩余油饱与度得认识,便于提高开发效果。

不同于地质上得油藏静态描述。

直接观察法4５６/模拟法７4.钻观察井:用在勘探初期或油田开发过程中,它可以直接取芯分析油层岩性与物性以及流体性质与在油层中得分布;5.直接测试:测井、井间地震、试井、井间示踪剂测试等从微观孔隙结构到宏观井间连通6.开辟开发实验区:油田开发初期为了达到某种目得(如提高采收率),要在油田内部选择一个有代表性得地区进行试验。

如大庆得小井距进行单层注水与各种提高采收率得方法得试验等;优点:1、直观(瞧得见,摸得着) ; 2、准确。

缺点:1、有一定局限性,范围小;2、成本高,周期长; 3、不能重复进行。

7.模型基本上有两种,一类就是您摸得着得物理模拟,另一类就是您不摸着数学模拟。

8.物理模拟:相似模型、单元模型9.数学模拟:水电相似模型、解析模型与数值模型、10.水电相似:多孔介质中得渗流过程与导电介质中电得流动过程, 相似得原理来进行模拟研究,应用越来越少。

11.数值模型:就是一种离散化得近似方法,常用得方法就是有限差分方法、数值方法求得得解不就是一个数学函数关系式,而就是分布在足够多得点上得一系列函数值、逼近近似解,使复杂得偏微分方程得求解成为可能、目前已经建立了功能强大得软件并在油气田开发中得到广泛得应用——油藏数值模拟、12.优点:能重复开发,可以进行所谓得“多次开发”;可以在短期内进行“开发”;成本相对比较低;可以模拟各种非均质条件与开发要求,避免了直接观察法得缺点、结果可以直接用于油田开发。

第一章油藏数值模拟方法分析1.1油藏数值模拟1.1.1油藏数值模拟简述油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。

其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。

其基础理论是基于达西渗流定律。

油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。

基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值(生产动态)与实际吻合。

其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。

充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。

这组流动方程组由运动方程、状态方程和连续方程所组成。

油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。

具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。

油藏数值模拟技术从50 年代的提出到90 年代间历经40 年的发展，日益成熟。

现在进入另外一个发展周期。

近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。

在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用非常广泛。

油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。

图1 油藏数值模拟流程图1.1.2油藏数值模拟的类型油藏数值模拟类型的划分方法有多种，划分时最常用的标准是油藏类型、需要模拟的油藏流体类型和目标油藏中发生的开采过程，也可以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定，油气藏特性和油气性质不同，选择的模型也不同，还可以根据油藏数值模拟模型所使用的坐标系、空间维数和相态数来划分。

数值模拟上机实验1：三对角系数矩阵的解法 隐式差分格式出发点在（i,n ＋1)，取关于t 的一阶向后差商和关于x 的二阶差商。

取对于一维渗流问题的隐式差分方程组的系数矩阵为三对角矩阵，追赶法（THOMAS ）就是用来求解三对角矩阵方程组的一种比较简单、 应用也极为广泛的解法。

它的基本思路是将三对角矩阵A 分解成两个特定形式的三对角矩阵的乘积。

in i i n ii n i i d pb pa pc =++++++-11111tp pxp p pn in in i n i n i ∆-=∆+-++-+++121111122tx δ∆=∆()n in i n in ipppp-=++-++++-1111121δδδ追赶法程序如下:Dim m As Integer, n As Integer, i As IntegerDim P(1 To 10) As Single, x(1 To 10) As Single, y(1 To 10) As Single, a(1 To 10) As Single, b(1 To 10) As Single, c(1 To 10) As Single, d(1 To 10) As SingleDim l(1 To 10) As Single, u(1 To 10) As Single n = InputBox("请输入方程个数") For i = 1 To na(i) = InputBox("a(" & i & ")=?") Print "a(" & i & ")="; a(i); Next i PrintFor i = 1 To n - 1b(i) = InputBox("b(" & i & ")=?")Print "b(" & i & ")="; b(i);1d 2d xN d id acb1p 2p xN p ipPrintFor i = 2 To nc(i) = InputBox("c(" & i & ")=?")Print "c(" & i & ")="; c(i);Next iPrintFor i = 1 To nd(i) = InputBox("d(" & i & ")=?")Print "d(" & i & ")="; d(i);Next iPrintl(1) = a(1)For i = 2 To nu(i - 1) = b(i - 1) / l(i - 1)l(i) = a(i) - c(i) * u(i - 1)Next iy(1) = d(1) / l(1)Print "y(1)="; y(1);For i = 2 To ny(i) = (d(i) - c(i) * y(i - 1)) / l(i) Print "y(" & i & ")="; y(i);Printx(n) = y(n)For i = n - 1 To 1 Step -1x(i) = y(i) - u(i) * x(i + 1)Next iFor i = 1 To nPrint "x(" & i & ")="; x(i);Next i运行结果:实验2:单相流数值模拟方法已知一维均质等厚单相渗流，已知边界条件定压、初始压力分布，求任意时刻的压力分布。