Precise estimation of the effective petrophysical characteristics for the oil or gas-bearing reservoir plays a vital role in production control problems. It is preferable to use nondestructive measurement methods for porosity, conductivity, geomechanical modulus, and some other parameter fields. For the estimation of anisotropic permeability tensor flows in different directions need to be simulated in the core plug digital image void space, which is very difficult (if possible) to conduct with the real rock sample. The creation of the digital core image includes three stages: construction of its internal structure based on the computed tomography (CT) sinogram, filtration, and segmentation. Routine practice is the filtration of 2D slices pack of the core plug CT image because of a lack of computational power and memory limitations. Unfortunately, this can generate a directional orientation error orthogonal to the pack of slices. The total 3D edge-preserving filtration with the usage of two approaches: modified implicit anisotropic diffusion and discrete orthogonal transforms can reduce this error. Universal code based on the MPI+OpenAcc programming paradigm was tested on different high-performance computing systems incorporating various accelerators like GPGPU and heterogeneous processors such as Sunway 26010.

No abstract