An improved deconvolution algorithm using B-splines for well-test data analysis in petroleum engineering
文献类型:期刊论文
作者 | Liu WC; Liu YW(刘曰武)![]() ![]() |
刊名 | JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING
![]() |
出版日期 | 2017 |
卷号 | 149页码:306-314 |
通讯作者邮箱 | liuyuewulxs@126.com |
ISSN号 | 0920-4105 |
产权排序 | [Liu, Wenchao; Liu, Yuewu; Han, Guofeng; Wan, Yizhao] Chinese Acad Sci, Inst Mech, Beijing 100190, Peoples R China; [Zhang, Jianye] Petro China, Tarim Oilfield Co, Res Inst Explorat & Dev, Korla 841000, Xinjiang, Peoples R China |
通讯作者 | Liu, YW (reprint author), Chinese Acad Sci, Inst Mech, Beijing 100190, Peoples R China. |
中文摘要 | Ilk et al.'s deconvolution algorithm using B-splines involves the Laplace transformation of the convolution equation with respect to production rate and wellbore pressure based on Duhamel principle. However, for common cases, the production rate function has "discontinuity" with respect to production time; it does not satisfy the precondition that the function to be transformed by Laplace transformation should be continuous. This inherent defect may directly cause enormous amount of computational time or even the failure of the numerical Laplace inversion in the deconvolution process. Based on these concerns, a fundamentally improved deconvolution algorithm using B-splines is presented here. In the convolution equation, the wellbore pressure derivative corresponding to constant unit production rate as the target of deconvolution is still represented by weighted summation of second-order B-splines; however, the computation process of the deconvolution is kept in the level of integral in the real time space instead of the Laplace space, for the reason that there will be no continuity requirement for the production rate function in the application of Duhamel principle for the deconvolution computation problem. According to the real production rate history, a technique of piecewise analytical integration is adopted for obtaining the elements of sensitivity matrix of a linear system with respect to weight coefficients; the linear system is generated by substituting the measured wellbore pressure data and corresponding variable 'production rate data into the convolution equation containing B-splines. The proposed direct analytical solution method of the integration for calculating the elements of the sensitivity matrix can not only guarantee the success of the deconvolution computation, but also can largely enhance the deconvolution computation speed. Moreover, in order to further improve the computation speed, a binary search method is also applied to find which production segments (with constant production rate) the measured wellbore pressure data points locate at in the deconvolution computation process. Another linear system with respect to weight coefficients for the regularization from Ilk et al.'s deconvolution algorithm is appended in order to overcome the effect of data errors. The two linear systems are combined together as an over-determined linear system, which can be solved by the least square method. Eventually, the reconstructed wellbore pressure and its derivative by B-splines corresponding to the constant unit production rate can be obtained. Numerical experimental tests demonstrate that the improved deconvolution algorithm exhibits good accuracy, computation speed and stability of data error tolerance. And the statement on how to perform the regularization when data error exists is also made in order to deconvolve the correct wellbore pressure derivative. The improved deconvolution algorithm is also applied into an actual field example. It is found that the deconvolution results by the improved deconvolution algorithm have good agreement with the ones by Von Schroeter et al.'s deconvolution algorithm and by Levitan et al.'s deconvolution algorithm as a whole; and the feature of typical log-log curves of the wellbore pressure drop and the wellbore pressure derivative corresponding to the improved algorithm is very close to the one of typical log-log curves calculated directly from the wellbore pressure data in the well shut-in period. In addition, through many numerical experimental tests, it is also concluded that as the quantity of data largely increases, the improved Ilk et al.'s deconvolution algorithm exhibits the big advantage in fast computational speed over von Schroeter et al.'s algorithm and Levitan et al.'s algorithm. |
分类号 | 一类 |
类目[WOS] | Energy & Fuels ; Engineering, Petroleum |
研究领域[WOS] | Energy & Fuels ; Engineering |
关键词[WOS] | Deconvolution algorithm ; Well testing ; Duhamel principle ; Least square method |
收录类别 | SCI ; EI |
原文出处 | http://dx.doi.org/10.1016/j.petrol.2016.10.064 |
语种 | 英语 |
WOS记录号 | WOS:000393007900028 |
源URL | [http://dspace.imech.ac.cn/handle/311007/60102] ![]() |
专题 | 力学研究所_流固耦合系统力学重点实验室(2012-) |
推荐引用方式 GB/T 7714 | Liu WC,Liu YW,Han GF,et al. An improved deconvolution algorithm using B-splines for well-test data analysis in petroleum engineering[J]. JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING,2017,149:306-314. |
APA | Liu WC,刘曰武,Han GF,Zhang JY,Wan YZ,&Liu YW.(2017).An improved deconvolution algorithm using B-splines for well-test data analysis in petroleum engineering.JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING,149,306-314. |
MLA | Liu WC,et al."An improved deconvolution algorithm using B-splines for well-test data analysis in petroleum engineering".JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING 149(2017):306-314. |
入库方式: OAI收割
来源:力学研究所
浏览0
下载0
收藏0
其他版本
除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。