浅谈线性方程组的数值解法及其应用.pdf

浅谈线性方程组的数值解法及其应用

1、相关定义

1.1、分形油藏基本概念

定义3.1[32]:维数为的分形渗透网嵌入到d(d=2,3)维岩块中,即整

个导流系统是一个分形体,称具有这种特性的油藏为分形油藏。df

3.1.1分形孔隙度θf和渗透率Kf(r)假设分形体内流体储集在体积

为的座点处(设每个座点体积相同),座点密度为。分形体中,座点孔

隙体积为常数。用描述某种相应对称性(如,VsN(r)VsBB=A2π

h和4π分别描述直线对称,圆柱面对称和球对称),a为位置-浓度参数

[33],为岩块的欧几里德维数,定义分形孔隙度dθf有:θf=aBVsrd

f?d(3-1)这说明分形网格的θf不再是常数,而是随波及半径r成幂

律关系。取,则有:wr=rdfdfwwrθθrθw=aVBsrwdf?d,得到分形

孔隙度=,θw为r=rw处的孔隙度。同样渗透率定义

为:Kf(r)=aVsBmrdf?d?θ(3-2)()r=rw处的渗透率Kw=aVsBmr

wdf?d?θ,得到渗透率Kfr=Kwrrwdf?d?θ。3.1.2分形参

数的物理意义(1)分形维数df分形维数df严格地是一个分形体的几

何特征,是分形体复杂程度的重要标志。一般认为,df值不同,复杂程

度也不一样。随复杂程度加剧,df值会愈高。15

1.2、分数阶的基本定义

从十七世纪分数阶微积分诞生之日起,数学家们就不断的探讨分数

阶算子的理论体系。后经多位数学家的努力,从不同的角度入手,建立

了多种不同形式的分数阶算子定义,现在主要通用的三种定义[1]形

式为:(1).Grümwald-Letnikov定义:对于任意的实数α,记α的整数部

分为[α]([α]为小于α的最大整数),假如函数f(t)在区间[α,t]上

有m+1阶连续的导数,α0时,m至少取[α],则定义分数阶α阶导

数为:()li0m0()nGaDtαftΔnhh→=t?ah?αi=∑????iα???ft?

ih(1-1)()(1)(2)(1)!iii其中,?α=?α?α+?α+L?α+?。故

上式可化为:()()()()()()(1)0()1()11GmkktmmatkaDftft

tfdkmαα==∑Γα?α+?α+?++Γ?α++∫?η?α++ηη(1-2)这个

定义是从寻找阶导数与n次积分的统一性出发,衍生出来的。可以

扩展为负整数。nn(2).Riemann-Liouville定义:Riemann-Liouville

定义的分数阶导数为:()(())1()1,01nnRatntnandfnNDf

tddtfdndtntααααηη?+ηα==∈Γ?∫?≤?

1.3、Bernsrein多项式的定义

定义2.2.2.1n次Bernsrein多项式定义[42]为如下形式:B,(x)n

xi(1)niinix(2-21)由(1x)ni的二项式展开得:nnii(1)ni(1)knnii

kixxkikx(2-22)0令(x)[B0,n(x),B1,n(x),L,Bn,n(x)]T,

则(x)An(x)(2-23)其中-12-(1)0n1nnn0nn00(1)001(1)0n0