CURVELET‐BASED 3D RECONSTRUCTION OF DIGITAL CORES USING THE POCS METHOD

With the development of shale-gas exploration and exploitation, exploitation is necessary to study the 3D spatial distribution of shale fractures for research on shale rock physics. Because of the limitation of instruments, accurate shale slice is discontinuous in depth, and the minimum interval between adjacent slices is inconsistent with horizontal resolution of digital cores. Those are the main factors which can prevent accuracy improvement of fracture characterization and physical modeling for digital cores. In order to study the 3D spatial distribution of fractures, this paper achieves reconstruction of 3D digital cores using the curvelet transform and Project Onto Convex Sets (POCS) method. Firstly, sample one per two slices in depth regularly for 3D dataset of sandstone obtained by X ray scanner, then reconstruct the 3D dataset using the proposed method. The reconstructed result is consistent with the original one and superior over the current method (spgl1), which approves the validity and superiority of the proposed method. Secondly, double the accurate 2D shale slices in depth that are obtained using focused ion beam scanning electron microscopy (FIB-SEM). It generates the 3D digital cores in which the minimum interval between adjacent slices in depth is basically consistent with the horizontal resolution, weakens the discontinuity of shale slices in depth caused by the instrument limitations and makes the fracture distribution more clear. Tests on sandstone and accurate shale rocks demonstrate the validity and superiority of the proposed method.

[1]  H Stark,et al.  High-resolution image recovery from image-plane arrays, using convex projections. , 1989, Journal of the Optical Society of America. A, Optics and image science.

[2]  D. Youla,et al.  Image Restoration by the Method of Convex Projections: Part 1ߞTheory , 1982, IEEE Transactions on Medical Imaging.

[3]  F. Herrmann,et al.  Non-linear primary-multiple separation with directional curvelet frames , 2007 .

[4]  Mauricio D. Sacchi,et al.  Beyond alias hierarchical scale curvelet interpolation of regularly and irregularly sampled seismic data , 2010 .

[5]  Michael P. Friedlander,et al.  Probing the Pareto Frontier for Basis Pursuit Solutions , 2008, SIAM J. Sci. Comput..

[6]  Jianmeng Sun,et al.  Reconstruction of 3-D digital cores using a hybrid method , 2009 .

[7]  R. Abma,et al.  3D interpolation of irregular data with a POCS algorithm , 2006 .

[8]  Mauricio D. Sacchi,et al.  Convergence improvement and noise attenuation considerations for beyond alias projection onto convex sets reconstruction , 2013 .

[9]  Jinghuai Gao,et al.  Curvelet-based POCS interpolation of nonuniformly sampled seismic records , 2012 .

[10]  Emmanuel J. Candès,et al.  The curvelet transform for image denoising , 2002, IEEE Trans. Image Process..

[11]  L. Bregman The relaxation method of finding the common point of convex sets and its application to the solution of problems in convex programming , 1967 .

[12]  S. Bakke,et al.  Process Based Reconstruction of Sandstones and Prediction of Transport Properties , 2002 .

[13]  Guochang Liu,et al.  Irregular seismic data reconstruction based on exponential threshold model of POCS method , 2010 .

[14]  Mike E. Davies,et al.  Iterative Hard Thresholding for Compressed Sensing , 2008, ArXiv.

[15]  R. Neelamani,et al.  Coherent and random noise attenuation using the curvelet transform , 2008 .

[16]  Laurent Demanet,et al.  Fast Discrete Curvelet Transforms , 2006, Multiscale Model. Simul..

[17]  Felix J. Herrmann,et al.  Non-parametric seismic data recovery with curvelet frames , 2008 .

[18]  T. Blumensath,et al.  Iterative Thresholding for Sparse Approximations , 2008 .

[19]  F. Herrmann,et al.  Simply denoise: Wavefield reconstruction via jittered undersampling , 2008 .

[20]  M. Blunt,et al.  Prediction of permeability for porous media reconstructed using multiple-point statistics. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Jean-Luc Starck,et al.  Very high quality image restoration by combining wavelets and curvelets , 2001, SPIE Optics + Photonics.

[22]  Jianmeng Sun,et al.  Numerical simulation of rock electrical properties based on digital cores , 2009 .

[23]  Ignace Loris,et al.  Nonlinear regularization techniques for seismic tomography , 2008, J. Comput. Phys..