IP Library › Granted Patent US 12,723,964
Granted Patent B2
US 12,723,964 · App. 18/256,007 · Granted Sep 1, 2026

Process for identifying a sub-sample and a method for determining the petrophysical properties of a rock sample

Inventor: Marco Miarelli (San Donato Milanese, IT)
Assignee: ENI S.P.A.
G01N15/0806G01N33/24G06V10/443G06V10/462G06V10/54G06V20/64G01N2015/0846
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,723,964
App. No.
18/256,007
Granted
Sep 1, 2026
Kind
B2
Abstract

A process to identify a representative sub-sample of a rock sample includes the following steps: digitally acquiring at low resolution first 2D or 3D images of the rock sample, and by a recomposition of the first images along a direction representing a digital sample; subdividing with a volumetric subdivision the digital sample defining 3D digital blocks having a volume of substantially similar value; analysing each 3D digital block on the basis of key points or texture descriptors obtained elaborating said plurality of the first images on the basis of the properties or structures of the rock, identifying two or more homogeneous classes of 3D digital blocks, each homogeneous class defining a portion of digital sample having 3D digital blocks equivalent to each other as to the technical characteristics of the rock; and selecting at least one representative block for each homogeneous class identified, locating and extracting a sub-sample from the rock sample.

Claims (35)

1 . A process for identifying a sub-sample representative of

a rock sample, which comprises:

digitally acquiring at low resolution a plurality of first 2D images (N slice , N el ) or 3D images of said rock sample, said plurality of first 2D images (N slice , N el ) or 3D images being adapted to represent a digital sample of said rock sample;

characterised in that it comprises:

a recomposition of the plurality of first 2D images (N slice , N el ) or 3D images along a direction (T) obtaining said digital sample;

subdividing with a volumetric subdivision said digital sample defining a plurality of three-dimensional digital blocks having a volume (V) of substantially similar value (V 1 );

analysing each three-dimensional digital block of the plurality of three-dimensional digital blocks on the basis of a plurality of key points or texture descriptors obtained elaborating said plurality of first 2D images (N slice , N el ) or 3D images on the basis of properties or structures of the rock sample, and identifying two or more homogeneous classes of three-dimensional digital blocks, each homogeneous class of the identified two or more homogeneous classes defining a portion of digital sample having three-dimensional digital blocks equivalent to each other as to technical characteristics of the rock;

selecting at least one representative block for each homogeneous class identified on said digital sample, localizing and extracting a sub-sample from said rock sample, said sub-sample comprising said at least one representative blocks for each homogeneous class.

2 . The process according to claim 1 , further comprising:

performing a rock analysis of said plurality of first 2D images (N slice , N el ) and defining a plurality of first processed 2D images (N key ) each comprising the plurality of key points or texture descriptors;

digitally processing said plurality of first processed 2D images (N key ) on the basis of a digital similarity in proximity to said key points or texture descriptors and defining a plurality of first clustered 2D images (N Clust ) wherein said key points or texture descriptors are subdivided into at least two groups of similar points;

representing said rock sample with a second digital sample superimposing, according to the direction T, said plurality of first clustered 2D images (N Clust );

performing the volumetric subdivision of said second digital sample defining said plurality of three dimensional digital blocks.

3 . The process according to claim 1 , further including:

statistically processing and subdividing said plurality of three dimensional digital blocks on the basis of the density of said key points or texture descriptors to define said two or more homogeneous classes of said plurality of three dimensional digital blocks;

extracting by a selection said at least one representative blocks from each of said two or more homogeneous classes using statistical analysis and/or filtering methods for screening said 3D digital blocks.

4 . The process according to claim 3 , further including providing for a preliminary processing which digitally processes said plurality of first 2D images (N slice ) or 3D images and which determines for each first 2D image (N slice ) or 3D image a corresponding first processed 2D image (N el ) or processed 3D image, said rock analysis being performed on said plurality of first processed 2D images (N el ) or processed 3D images so as to define first processed 2D images (N key ) or processed 3D images comprising said key points or texture descriptors and corresponding regions of interest or volumes of interest surrounding said key points or texture descriptors, a subsequent grouping of said key points or texture descriptors, on the basis of the technical characteristics in said regions of interest or volumes of interest, identifies first clustered 2D images (N Clust ) or clustered 3D images.

5 . The process according to claim 1 , wherein the localizing comprises mapping said at least one representative block on the digital sample subdivided into said two or more homogeneous classes and identifying by means of an affine transformation said at least one representative block in said rock sample in order to identify said sub-sample.

6 . The process according to claim 1 , wherein said value (V 1 ) is determined with a statistical analysis through suitable processing and/or statistical algorithms of the technical characteristics identified by said similar key points or texture descriptors included in said 3D digital blocks and/or in that said value (V 1 ) has a minimum volume value (V 1 min ) referred to a predefined computational accuracy.

7 . A method for determining the physical or petrophysical properties of a rock sample, the method including the following steps:

identifying and extracting from said rock sample at least one sub-sample using the process of claim 1 , said at least one sub-sample comprising at least one representative block for each of the two or more homogeneous classes of three-dimensional digital blocks,

digitally acquiring, at high resolution, a plurality of third digital images of said at least one sub-sample at said at least one representative block, and

performing a physical or petrophysical analysis of said plurality of third digital images acquired to determine the physical or petrophysical properties of said rock sample starting from properties of said at least one representative block for each of said two or more homogeneous classes.

8 . The process according to claim 1 , further comprising the step of extracting the sub-sample of the rock sample having at least two representative blocks, said at least two representative blocks having the volume (V) of substantially similar value.

9 . A data processing system comprising:

a first tomograph configured to perform a first digital scan at low resolution on a rock sample defining a plurality of first 2D images (N slice ) or 3D images, and

a first processor configured to acquire said plurality of first 2D images (N slice ) or 3D images, and to execute the process for identifying the sub-sample representative of said rock sample according to claim 1 , and

comprising:

a second tomograph configured to perform a second digital scan at high resolution of the sub-sample defining a plurality of third 2D images or 3D images, and

a second processor configured to execute a method for determining the physical or petrophysical properties of said rock sample, the method including the following steps:

identifying and extracting from said rock sample at least one sub-sample, said at least one sub-sample comprising at least one representative block for each of the two or more homogeneous classes of three-dimensional digital blocks,

digitally acquiring, at high resolution, a plurality of third digital images of said at least one sub-sample at said at least one representative block, and

performing a physical or petrophysical analysis of said plurality of third digital images acquired to determine the physical or petrophysical properties of said rock sample starting from properties of said at least one representative block for each of said two or more homogeneous classes.

10 . A non-transitory memory storing a computer program, the computer program comprising instructions which when the program is executed by a computer, the computer executes the process according to claim 1 .

11 . A non-transitory memory storing a computer program, the computer program comprising instructions which when the program is executed by a computer, the computer executes the method according to claim 7 .

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jun 20, 2023
From: MIARELLI, MARCO
To: ENI S.P.A.
Reel/Frame 063999/0257 →
Priority Claims (1)
IT 102020000029744 · Dec 3, 2020 · national
Continuity (1)
Related Publication 20240053246A1 · Feb 15, 2024
References Cited (38)
US 8515721B2 · Morton · 2013 [cited by examiner]
US 9507047B1 · Dvorkin · 2016 [cited by examiner]
US 10198804B2 · Sungkorn · 2019 [cited by examiner]
US 10223782B2 · Sungkorn · 2019 [cited by examiner]
US 10514372B2 · De Prisco · 2019 [cited by examiner]
US 10830713B2 · Zhang · 2020 [cited by examiner]
US 20100128932A1 · Dvorkin et al. · 2010 [cited by applicant]
US 20100228486A1 · Wu · 2010 [cited by examiner]
US 20120221306A1 · Hurley · 2012 [cited by examiner]
US 20130179080A1 · Skalinski · 2013 [cited by examiner]
US 20130262028A1 · De Prisco · 2013 [cited by examiner]
US 20150331145A1 · Grachev · 2015 [cited by examiner]
US 20160307312A1 · Sungkorn · 2016 [cited by examiner]
US 20170018096A1 · Sungkorn · 2017 [cited by examiner]
AU 2017239499A1 · 2017 [cited by examiner]
WO WO2015153505A1 · 2015 [cited by examiner]
WO WO2015153506A1 · 2015 [cited by examiner]
Estimating permeability from thin sections without reconstruction: Digital rock study of 3D properties from 2D images, Nishank Saxena et al., Elsevier, 2017, pp. 79-99 (Year: 2017). [cited by examiner]
Development of a Digital Rock Physics workflow for the analysis of sandstones and tight rocks, I. Verri et al., Elsevier, 2017, pp. 790-800 (Year: 2017). [cited by examiner]
Estimation of Petrophysical Parameters from 3D images of Carbonate Core, R. M. Sok et al., Apr. 2007, pp. 1-15 (Year: 2007). [cited by examiner]
Digital Rock Physics, Hasan Al-Marzouqi, IEEE, 2018, pp. 121-131 (Year: 2018). [cited by examiner]
Identifying the Representative Elementry Volume for Permeability in Heterolithic Deposits Using Numerical Rock Models, Kjetil Nordahl et al., Springer, 2008, pp. 753-771 (Year: 2008). [cited by examiner]
Determination of petrophysical properties of sedimentary rocks by optical methods, D. Korte et al., Elsevier, 2017, pp. 72-79 (Year: 2017). [cited by examiner]
Petrophysical properties of porous medium from Petrographic Image Analysis data, Adrian Cerepi et al., Elsevier, 2001, pp. 233-256 (Year: 2001). [cited by examiner]
United Arab Emirates Office Action for Application No. P6001290/2023, dated Dec. 13, 2024, 9 pages. [cited by applicant]
Al-Marzouqi, “Digital Rock Physics: Using CT scans to compute rock properties”, IEEE Signal Processing Magazine vol. 35, No. 2, Mar. 1, 2018, pp. 121-131. [cited by applicant]
International Search Report for International Application No. PCT/IB2021/061223, dated Mar. 3, 2022, 4 pages. [cited by applicant]
Written Opinion for International Application No. PCT/IB2021/061223, dated Mar. 3, 2022, 12 pages. [cited by applicant]
Al-Raoush R. and Papadopoulos A., 2010-Representative elementary volume analysis of porous media using X-ray computed tomography. Powder Technology 200:69-77, DOI 10.1016/j.powtec.2010.02.011. [cited by applicant]
Arthur D, Vissilvitskii S (2007) K-means++: the advantages of careful seeding. Proc of the Annual ACM-SIAM Symposium on Discrete Algorithms 8:1027-1035, DOI 10.1145/1283383.1283494. [cited by applicant]
Blunt M. et al 2013—Pore-scale imaging and modelling. Advances in Water Resources 51:197-216, DOI 10.1016/j.advwatres.2012.03.003. [cited by applicant]
Clewlow L. and Strickland C., 2000-Energy derivatives: Princing and risk management. [cited by applicant]
Mostaghimi et al., 2012—Computations of absolute permeability on micro-ct images; Mathematical Geosciences 45, 24 pages, DOI 10.1007/s1 1004-012-9431-4. [cited by applicant]
Nordahl K. and Ringrose P, 2008—Identifying the representative elementary volume for permeability in heterolithic deposits using numerical rock models. Mathematical Geosciences 40:753-771. [cited by applicant]
Parrinello CM et al, 2016—Interactive Outlier Removal: A method for Identifying Outliers in Laboratory Recalibration Studies. Clinical Chemistry 62(7):966-972 DOI 10.1373/clinchem.2016.255216. [cited by applicant]
Rousseeuw PJ, 1987—Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. Journal of Computational and Applied Mathematics 20:53-65, DOI https://doi.org/10.1016/03777-0427(87)90125-7. [cited by applicant]
Verri et al. 2017—Development of a digital rock physics workflow for the analysis of sandstones and tight rocks. Journal of Petroleum Science and Engineering 156:790-800, DOI https://doi.org/10.1016/j.petrol.2017.06.053. [cited by applicant]
Saudi Arabia Office Action for Application No. 523440988, dated Jan. 14, 2025, 17 pages with translation. [cited by applicant]