水平井试井模型教案分析.doc

一 符号说明和无量纲量 字符说明 英文符号： B = 流体积系数； C = 井筒存储系数（m3/MPa）； ct = 总压缩系数（1/MPa）； h = 有效层厚（m）； j(v) = j 阶第类修正Bessel函数；kj = j方向渗透率（μm2），j = h、v； Kj() = j 阶第二类修正Bessel函数； L = 水平井筒半长（m）； P = 压力（MPa)； △P = 压差（MPa)； q = 产量（m3/d）； rw = 水平井筒半径（m）； Skin = 表皮因子； s = Laplace变换量； t = 时间（hour）； x、y = 坐标； xe、ye = 外边界距离（m）； z = 坐标； zw = 水平井垂向位置，m； 希腊符号： μ =流体粘度（mPa·s）； ??= 孔隙度； ?j?= j 方向扩散系数，定义为?j = kj /??ctj； 下标： D = 无量纲； w = 井壁； i = 初始状态； eD = 无量纲外边界无量纲量定义 采用中国SI单位制，定义如下无量纲量群： 无量纲压力： ；； 表皮压力降： ， 无量纲时间： ， 无量纲井筒存储系数： ， 无量纲Fair变井筒存储参数： ， 无量纲坐标： ，，，， 无量纲井筒半长： ；； 无量纲渗透边界补给系数： 无量纲双孔介质参数： ； 二 试井分析理论模型 由于试井分析模型较多，本章将分类介绍主要模型。 模型1 平面无限延伸平板地层 物理模型如图所示，垂向有上下平板型封闭、上定压下封闭、上封闭下定压及上下皆定压四种情形可选，水平方向无限延伸，均匀各向异性或双重孔隙介质任选其一，常井储和Fair变井储任选其一。物理模型如图所是2-1A~D。 无量纲参数：LD、zwD、HLCD、CDe2S、βD、C(D、ω、λ； 需要说明的是，由于定压边界的影响，在某些情况下上下定压和上封闭下定压、下封闭上定压情形区别不大。 压力分布解式：考虑垂向存在平板型封闭边界、水平方向无穷延伸的均匀或各向异性地层中一口均匀流量水平井的不定常渗流理论，它是水平井不定常渗流基本问题之一。关于这个专题以前的研究，较早而且比较完整的结果主要有Goode和Thambynayagam（1985）、Daviau（1985）、Clonts和Ramey（1986）、Ozkan和Raghavan（1987、1991）以及Kuchuk和Goode（1988）等人的论文，他们的结果各有特点。而针对双重介质的研究者有Carvalho（1988）、Aguilera（1989）以及刘慈群（1991）等人，以下分别简述之。 按SPE的文献发表时间排序，Goode和Thambynayagam的论文是占先的。他们将水平井看成垂直于x-y平面的条带源，采用Fourier余弦变换结合Laplace变换求解了三维Fourier方程，分别给出了压力降落和压力恢复问题的实时域解式。虽然他们的计算分析过于简单，也没有进一步给出快速计算曲线方面的研究结果，但他们给出了压力降落早期、中期、中晚期及晚期的近似表达式，这对后来的研究者有很大的启发。另外他们追溯了利用水平井开采的思想直至四十年代，并列举了包括东欧研究者在内的一些论文。同年，Daviau等人也研究了类似问题。他们利用已有的Gringarten和Ramey（1973）的瞬时点源函数式直接构造了实时域解式。Daviau等人的理论研究比较简洁，采用一种近似方法考虑了井筒存储的影响，但他们计算的主要参数影响的半对数曲线图版流动段很清楚，基本反应了这类问题的压力反应特征。Ozkan和Raghavan的论文更注重计算方法方面的研究，对解式的渐近变化表述的比较清楚。尤其是他们对压力导数曲线的分析和计算是前所未有的，所附录的数值计算结果也为后来的研究者提供了验证方便。相比来看，Kuchuk等人的研究结果集众家之长而富有代表性。他们在数学模型中吸取了Goode和Thambynayagam论文中的精华，又引入了一种新的混合型垂向边界条件，利用其极限形式研究了有底水或气顶的各向异性介质中水平井的瞬时压力。他们的解析解式几乎包含了在此以前的有关结果。其中值得注意的是他们采用沿着水平井段取积分平均的方法解决了测压点的选取问题。 物理模型：物理模型和Green-Newman解法思路如图所示： 线源有量纲解式：按Green-Newman解法，有量纲解式为： 式中，根据井和所在的地层条件，按前文的推导结果，可给出各维Green函数如下： （1）x方向为半长是L的直线源，取算术平均： （2）y方向要么为无限延伸： （3）z方向要么上下平板型封闭，要么有一个为定压边界而另一个为平板型封闭： 线源无量纲解式：有两种方法可以对褶积通式进行无量纲化（Kuchuk，1988；Ozka