试井分析中的离散数据拉普拉斯变换方法.docxVIP

试井分析中的离散数据拉普拉斯变换方法 试井分析中的离散数据拉普拉斯变换方法 开发试采 天 然 气 工 业 2021年3月 试井分析中的离散数据拉普拉斯变换方法 张杰 李翠楠 张向阳 李治平 (1. 西南石油学院 2. 中国石油西南油气田分公司重庆气矿 3. 中国地质大学 北京) 张杰等. 试井分析中的离散数据拉普拉斯变换方法. 天然气工业, 2021; 24(3) :100~102 摘 要 随着计算机技术和试井分析理论的日益发展, 人们将拉普拉斯数值变换方法引入到试井分析中, 并得到越来越广泛的应用, 拉普拉斯变换已经成为处理试井测试数据的一种重要手段。以前人们研究较多的是基于连续函数的拉普拉斯变换方法, 然而, 试井测试所获得的数据通常都是某一时间区间上受到一定测试噪声影响的有序离散数据点列, 所以, 以前的那些方法不能直接用于处理试井测试数据。为此, 文章在考虑试井测试数据特点和拉普拉斯变换对原函数要求的基础上, 通过将整个半无穷界时间域分成三段, 对每个时间段分别采取合适的离散数据积分算法, 推导建立了适用于试井分析的离散数据拉普拉斯变换方法。该方法具有方便快捷、计算过程简单、计算稳定可靠的特点, 特别适合于在计算机上进行处理。 主题词 油气田数据 拉普拉斯变换 测试 噪声 试井解释 分析 1 拉普拉斯变换的定义 设函数f (t ) 当t 0时有定义, 而且积分 个测试点之前的时间区间; 第二段[t 1, t N ], 是测试时间区间, 该部分包含较多的有效信息; 第三段(t N , ) , 是最后一个测试时间点之后的时间区间: = f (x ) exp(-st ) d t (t 为复参变量) 在复平面上某一 区域内收敛, 则此积分所确定的函数: F (s ) = f (t ) exp(-1 f (t ) exp(-st ) d t + f (t ) exp(-st ) d t st ) d t (2) f (x ) exp(-st ) d t (1) 这里需要指出的是, 如果t 1对应的时间点就是0时刻, 那么, 整个半无穷时间区间将只被划分为[0, t N ]和(t N , ) 两个时间区间。 为了实现测试数据的拉普拉斯变换, 只需给出三个积分所对应的算法即可。每个积分对应的算法都有很多, 不同的算法组合将推导得到不同的拉普拉斯变换方法。 第一个积分就是对第一个测试时间点t 1进行向前外推。因为对应于时间0点的数据通常是知道的, 例如, 压力恢复测试t =0时 p =0, 所以端点t 1的外推变成了插值。如果应用简单的线性插值方法来处理第一个积分部分, 那么 I 1= 这称为函数的拉普拉斯变换, 用f 表示。 从该定义可以很明显的看出, 拉普拉斯变换是连续函数f (t ) 在半无穷界(0, ) 上的函数变换。 2 离散数据的拉普拉斯变换 在压力不稳定测试中, 人们获得的数据通常是在一个有限时间区域[t 1, t N ]内的有序离散时间点t i (i =1, 2, , N ) 上的压力值(p wf i ) , 或压力变化值( p i ) , 有时可能还有流量数据(q i ) 。这样的数据表现为基本函数f (t ) 的受噪音影响的离散数据点, 这里的f (t ) 既可以是测试的压力序列, 也可以是测试的流量序列。因此, 为了将测试所得的离散数据点变换到整个半无界拉普拉斯域上, 首先必须对拉普拉斯变换进行处理, 即按照测试起止时间点将整个半无穷界时间区间分成三段, 第一段(0, t 1) , 是第一 f (x ) e d t f 0f 1e =S -S d 0(1-e +S 作者简介:张杰, 1976年生, 西南石油学院在读博士。地址:(610500) 四川省成都市新都区西南石油学院博2002级。电话:(028)E -mail:s wpu@tom. com 第24卷第3期 天 然 气 工 业 开发试采 式中:d 0为早期段直线的斜率。 第二个积分部分包含的信息较多, 测试所得的数据绝大部分都在这个部分, 即 f (t ) e d t = f (t ) e d t 如果在测试数据任意两点之间插值, 将得到: f i e S f i +1e - 式中:d i 为相邻两点间直线段的斜率。 第三个积分包含的信息较少, 处理起来比较复杂, 需要对测试末期的数据进行外推。为了避免处理无穷界的不便, 首先进行如