扇形区域内有井储和表皮的井底压力.docVIP

扇形区域内有井储和表皮的井底压力 摘要 利用边界元方法求解扇形区域内有井存储和表皮效应的井底瞬时压力，给出了边界封闭或定压时，无量纲井底瞬时压力及其导数的双对数图；分析了扇形角度，井位对井底压力曲线的影响。比较了Laplace空间上的边界元解与实空间上的边界元解。 主题词 试井 边界元偏移（扇形区域） 井底压力 导数 方程 一、前 言 在现代试井分析中，无论是采用计算机手动拟合方法[1]、还是采用非线性回归方法[2]来反求参数（如渗透率、表皮因子、井筒存储等），都要涉及到理论曲线（即样板曲线）的计算。目前计算样板曲线的方法大都是解析方法或半解析法，这只能适用边界极简单的方程求解，因而要想解决任意外边界形状的试井分析问题只能借助数值方法。最常用的数值方法就是差分法，但差分法对复杂边界形状处理困难。有限元法可以解决复杂形状的外边界问题，但有限元方法的计算量和存储量太大。总之，无论是差分法，还是有限元方法，都要涉及到时间及空间网格的划分，都要储存大量的中间变量，因而就必须用较大内存的计算机。如果再考虑井筒存储和表皮因子，问题就变得更加复杂，其求解速度和解的精确度有可能不能满足试井分析的需要。 除上述两种数值方法外，边界元方法也可以用来求解复杂外边界形状地层的瞬时井底压力，KlKanl[3]等人用边界元方法求解了任意外边界形状点源井的解。KlKanl等人认为边界