反褶积方法在试井解释中的应用.docxVIP

反褶积方法在试井解释中的应用

反褶积的应用

抗弯环比（也称为抗弯环比）是一种广泛用于积分变换的数学方法。用褶积解决试井解释中的问题,早就取得了很好成果;而反褶积,直到最近,Schroeter、Hollaender和Gringarten等人解决了其计算方法上的稳定性问题,使反褶积方法很快引起了试井界的广泛注意。

有专家认为,反褶积的应用是试井解释方法发展史上的又一次重大飞跃。他们预言,随着测试新工具和新技术的增加和应用,以及与其它专业研究成果的更紧密结合,试井在油气藏描述中的作用和重要性必将不断增大。

测量方法的确定

由实测资料构造出理想的、等效的、对应于相同时间段内以恒定产量生产条件下的压力变化,由此得到测试全部历程的压力响应,其探测范围要大得多,而且不存在叠加计算所带来的影响,也不存在产量史的不完整而造成的误差。寻找到这些信息,便可获得比常规试井解释更多、更可靠的测试解释结果——这就是反褶积的大意和意义。

1.复的褶积方程

从测试录取得到的资料是产量变化史q(t)和压力变化史Δp(t),这就是试井解释的对象,而试井解释的任务就是从这些资料识别测试层和测试井的类型并算出它们的各种参数。实际上,由实测的压力变化和产量资料计算dΔpudlnt=t?Δpu′dΔpudlnt=t?Δpu′,是解释过程中至关要紧的一个步骤。

根据杜哈美(Duhamel)原理,得到研究压力降落和压力恢复的褶积方程为

p(t)=pi-∫t0q(τ)?Δpu(t-τ)?(t-τ)dτ=pi-∫t0q(τ)Δp′u(t-τ)dτ(1)Δp(t)=pi-p(t)=∫t0q(τ)?Δpu(t-τ)?(t-τ)dτ=∫t0q(τ)Δp′u(t-τ)dτ(2)p(t)=pi?∫t0q(τ)?Δpu(t?τ)?(t?τ)dτ=pi?∫t0q(τ)Δp′u(t?τ)dτ(1)Δp(t)=pi?p(t)=∫t0q(τ)?Δpu(t?τ)?(t?τ)dτ=∫t0q(τ)Δp′u(t?τ)dτ(2)

其中Δpu=Κh0.9210qμBΔpΔp=pi-pwf(t)其中Δpu=Kh0.9210qμBΔpΔp=pi?pwf(t)

式中:pwf(t)——开井生产t时刻井底压力,MPa;

t——测试井开井生产时间,h;

q——测试井产量,m3/d;

pi——原始地层压力,MPa;

Δpu——重整压力或单位产量下重整压力响应(以恒定1单位产量生产所造成压力变化)。

(2)式中,测试录取的两项资料,其中Δp(t)是褶积,q(t)则是构成此褶积的一个函数。要通过解释得到构成此褶积的另一个函数dΔpudt=Δpu′dΔpudt=Δpu′,进而得到dΔpudlnt=t?Δpu′dΔpudlnt=t?Δpu′,此过程就是反褶积。

2.反褶积公式及其求解

Schroeter方法即是Schroeter、Hollaender和Gringarten反褶积计算方法的简称。这里只介绍此方法梗概,其计算方法可参阅参考文献。

通常是用试井解释诊断曲线——压力导数曲线(lgdΔp(t)dlnt=lg[tdΔp(t)dt]—lgt(lgdΔp(t)dlnt=lg[tdΔp(t)dt]—lgt曲线)来识别油藏类型。因而,想要得到的是dΔpudlntdΔpudlnt。为此,Schroeter等人参照压力导数曲线,定义了两个新变量,即

σ=lnΔtz(σ)=ln[dpu(Δt)dlnΔt]=ln[dpu(σ)dσ]

事实上它们就是压力导数曲线的两个坐标,则

Δt=eσdΔt=eσdσdpu(Δt)dlnΔt=dpu(σ)dσ=ez(σ)dΔp(Δt)d(Δt)=1ΔtdΔpu(Δt)dln(Δt)=1eσez(σ)

根据褶积的性质,有

f(t)?g(t)=g(t)?f(t)

褶积方程(1)和(2)可分别写成

p(t)=pi-∫lnt-∞q(t-eσ)ez(σ)dσ(3)Δp=∫lnt-∞q(t-eσ)ez(σ)dσ(4)

这里,构成褶积的两个函数是q(t-eσ)和ez(σ),它们都是复合函数。

运用上述褶积公式和最优化过程(此处Schroeter方法是非线性回归方法),找出一个分段线性函数z(σ),也就是一条z(σ)—σ曲线,要求它能很好地拟合实测和理论的压力导数曲线,而拟合的时间范围是自生产开始直至最后一个数据录取点,也可以只包括若干个实测压力可靠的时间段。这就是试井解释中的反褶积。

在这一反褶积中,最主要的“未知元”是z(σ)—σ曲线。还有两个附加可选“未知元”,一个是原始地层压力pi,另一个是产量史容许误差限(在最优化过程中,需要它控制其收敛)。最优化过程就是通过不断调整这三个未知