渗流力学问题中的数值反演解ξπ∫π∫.pdfVIP

第 14 卷　第 1 期 应 用 力 学 学 报 . 14　 . 1 V o l N o 1997年3月 CH INESE JOURNAL OF APPL IED M ECHAN ICS M ar. 1997 渗流力学问题中的数值反演解 何光渝　　　　　郑书英 ( ) ( ) 西安石油学院　西安　710061 　 西北建筑工程学院　西安　710061 摘　　要 介绍了一种拉普拉斯变换数值反演的新方法——C rum p 方法以及它在渗流力 学问题中的应用。在数值解例中, 将 C rum p 方法和精确解及 Stehfest 方法进行了比 较, C rum p 方法的解与精确解吻合良好, 较之 Stehfet 方法有很大改进。在油气井现 代试井分析中, 用C rum p 方法绘制的典型曲线更精确。 关键词: 渗流力学; 拉普拉斯数值反演 1　前　言 在渗流力学中, 很多问题采用拉普拉斯变换来解决多孔介质中物质输运问题。简单的问题 可解析反演求得解析解; 而能应用于生产实践稍复杂的问题, 无法解析反演, 一般均采取数值 [ 1 ] 及演求解。目前国内外经常采用的数值反演方法是 Stehfest 方法 , 该方法简单易行, 计算时 间短; 但在一些有非单调解的问题中会产生振荡、不收敛以致于无法得到数值反演解。本文介 [ 2 ] 绍的C rum p 方法 既保持了 Stehfest 方法简单易行的优点, 又克服了上述缺点。采用 C rum p 方法应用于石油工程中的试井问题, 取得了良好的效果。 2　C rump 方法 众所周知, 拉普拉斯变换及其反演公式为: ∞ - s t 1 a+ i∞ s t F ( ) e f ( ) t, 　　　f ( ) e F ( ) s ( ) s = ∫ t d t = 2i ∫ s d 1 0 a- i∞ 式中 a ( ) ( ) ( ) 0 是任意的, 但必须比 F s 的所有奇点的实部要大。由 1 式可得 2ea t ∞ ( ) ( ( ) )