An improved TSVD-GCV inversion algorithm of pore size distribution in time-domain induced polarization using migration Hankel matrix

Reservoir porosity obtained from time-domain induced polarization (TDIP) well logging plays a vital role in estimating the hydraulic properties and obtain the reservoir parameters in a water-flood oil-field. Improving the inversion accuracy of the reservoir porosity can enhance the oil recovery in t...

Full description

Saved in:
Bibliographic Details
Published in:Journal of petroleum science & engineering Vol. 183; p. 106368
Main Authors: Zhou, Kaibo, Zhang, Shuo, Huang, Zhen, Zhang, Jianyu
Format: Journal Article
Language:English
Published: Elsevier B.V 01.12.2019
Subjects:
Backus-Gilbert theory
Migration Hankel matrix
Pore size distribution
Time-domain induced polarization
TSVD-GCV
Migration Hankel matrix
Backus-Gilbert theory
Pore size distribution
TSVD-GCV
Time-domain induced polarization
ISSN:0920-4105, 1873-4715
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Reservoir porosity obtained from time-domain induced polarization (TDIP) well logging plays a vital role in estimating the hydraulic properties and obtain the reservoir parameters in a water-flood oil-field. Improving the inversion accuracy of the reservoir porosity can enhance the oil recovery in the water-flood oil-filed. Evaluating reservoir pore size distribution through induced polarization decay curve is confronted with the problems of poor applicability of data pre-processing, low accuracy and lacking of evaluation criteria for inversion results of pore size distribution. The basic principles of TDIP are introduced and the relationship between pore relaxation time and pore diameter is given. Combining the mathematical characteristics of polarization decay curve data, the performance and the limitations of existing pre-processing algorithms are analyzed and pointed out, respectively. An improved data pre-processing algorithm using the spatial characteristics of linear transformation based on migration Hankel matrix is proposed, and this method improves the inversion accuracy of pore size distribution greatly. In the engineering application, 2-logarithmic sampling method is proposed to sample the polarization decay data for more efficient petroleum exploration with less sample points. The different regularization methods, regularization matrix and regularization parameter determination methods are compared and analyzed for the inversion of the pore size distribution. The numerical simulation experimental results show that the stability and accuracy of Truncated Singular Value Decomposition - Generalized Cross Validation (TSVD-GCV), Truncated Singular Value Decomposition - L Curve (TSVD-L) and Tikhonov-I-L are appropriate for the inversion of pore size distribution. Because of the truth that the pore size distribution of rock is unknown, Backus-Gilbert (BG) theory is introduced to evaluate the inversion results of rock polarization decay curve data of a mining area in Jilin Province. The rock sample experiment shows that the TSVD-GCV inversion algorithm has the best performance. •Separating noise from sampling signal by introducing migration Hankel matrix to migrate feature space.•Obtaining high inversion accuracy of pore size distribution within less sampling points using 2-logarithmic sampling method.•The TSVD-GCV algorithm has the best pore size distribution inversion performance by introducing the BG theory.
ISSN:0920-4105
1873-4715
DOI:10.1016/j.petrol.2019.106368