Method to generate numerical pseudocores using borehole images, digital rock samples, and multi-point statistics

• US 20090259446A1
• Filed: 04/08/2009
• Published: 10/15/2009
• Est. Priority Date: 04/10/2008
• Status: Active Grant

• Current Assignee: Schlumberger Technology Corporation (Schlumberger NV)
• Original Assignee: Schlumberger Technology Corporation (Schlumberger NV)
• Inventors: Zhang, Tunfeng, Hurley, Neil Francis, Zhao, Weishu
• Granted Patent: US 8,725,477 B2
• Time in Patent Office: Days
• US Class Current: 703/2
• CPC Class Codes: E21B 49/00 Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
• CPC Class Codes: G01V 1/44 using generators and receivers in the same well G01V1/52 takes precedence
• CPC Class Codes: G06F 2111/08 Probabilistic or stochastic CAD
• CPC Class Codes: G06F 2111/10 Numerical modelling
• CPC Class Codes: G06F 30/20 Design optimisation, verification or simulation optimisation, verification or simulation of circuit designs G06F30/30

Abstract

Methods and systems for creating a numerical pseudocore model, comprising: a) obtaining logging data from a reservoir having depth-defined intervals of the reservoir, and processing the logging data into interpretable borehole image data having unidentified borehole image data; b) examining one of the interpretable borehole image data, other processed logging data or both to generate the unidentified borehole image data, processing the generated unidentified borehole image data into the interpretable borehole image data to generate warped fullbore image data; c) collecting one of a core from the reservoir, the logging data or both and generating a digital core data from one of the collected core, the logging data or both such that generated digital core data represents features of one or more depth-defined interval of the reservoir; and d) processing generated digital core data, interpretable borehole image data or the logging data to generate realizations of the numerical pseudocore model.

Claims

1. A method for creating a numerical pseudocore model comprising:

a) obtaining logging data from a reservoir that includes depth-defined intervals of the reservoir, and processing the logging data into a portion of at least one interpretable borehole image data having unidentified borehole image data;

b) examining one of the portion of the at least one interpretable borehole image data, other processed logging data or both to generate the unidentified borehole image data, then processing the generated unidentified borehole image data into the portion of the at least one interpretable borehole image data so as to generate at least one warped fullbore image data;

c) collecting one of at least one core from the reservoir, a portion of the logging data or both and generating a digital core data from one of the collected at least one core, the portion of logging data or both such that the generated digital core data represents one of features or structures of one or more depth-defined interval of the reservoir; and

d) processing a portion of the generated digital core data, the portion of at least one interpretable borehole image data or a portion of the logging data to generate realizations of the numerical pseudocore model.

2. The method of claim 1, wherein the portion of the at least one interpretable borehole image data is a training image that is two-dimensional (2D) scalar arrays of continuously variable numerical values.

3. The method of claim 1, wherein the portion of the at least one interpretable borehole image data consists of processed raw data that consists of measured values and non-measured values.

4. The method of claim 1, wherein the unidentified borehole image data represents one of non-measured values of the portion of the at least one interpretable borehole image data or data gaps in the portion of the at least one interpretable borehole image data.

5. The method of claim 4, wherein the data gaps is from the group consisting of one of at least one damaged pad in the reservoir, at least one damaged area in the reservoir, at least one pad with inadequate pad pressure against a borehole wall in the reservoir, at least one pad obstructed from contacting the borehole wall in the reservoir or at least one inoperable pad in the reservoir.

6. The method of claim 1, wherein the other processed logging data is from the group consisting of one of logging-while-drilling data or wireline logging data, or some combination thereof.

7. The method of claim 6, wherein the wireline logging data includes at least one of gamma ray, density, sonic, neutron, caliper or resistivity logs.

8. The method of claim 1, wherein the at least one warped fullbore image data is one of two-dimensional (2D), three-dimensional (3D) or both.

9. The method of claim 1, wherein the generating at least one fullbore image consists of:

1) selecting a depth-defined interval of the borehole image log from the depth-defined intervals of the reservoir, processing the depth-defined interval of the borehole image log utilizing a multi-point statistics (MPS) algorithm, wherein the MPS algorithm provides for a modeling capturing unidentified geological structures from the portion of the at least one interpretable borehole image data so as to identify data locations in the portion of the at least one interpretable borehole image data;

2) processing the identified data locations using filter scores to group and then simulate patterns in the unidentified borehole image data whereby the unidentified borehole image data represents one of non-measured values of the portion of the at least one interpretable borehole image or data gaps in the portion of the at least one interpretable borehole image;

3) identifying petrophysical facies of the reservoir, wherein the petrophysical facies is from the group consisting of one of fractures, pores, vugs, conductive patches of rocks in a borehole wall of the reservoir, resistive patches of rocks in the borehole wall of the reservoir, bed boundaries or a rock matrix; and

4) processing the unidentified borehole image data into the portion of the at least one interpretable borehole image so as to generate the at least one fullbore image, then with a known borehole diameter warping the at least one fullbore image into an original three-dimensional (3D) shape whereby the numerical pseudocore is conditioned so as to match one of the warped portion of the fullbore image or the portion of the at least one interpretable borehole image.

10. The method of claim 1, wherein the logging data is from the group consisting of one of logging data having multiple depths of investigation, analog models of core data from one or more reservoirs or some combination thereof.

11. The method of claim 1, wherein the generating of the digital core data consists of identifying petrophysical facies from one of the collected at least one core or the portion of logging data, wherein the petrophysical facies is from the group consisting of one of fractures, pores, vugs, conductive patches of rocks in a borehole wall of the reservoir, resistive patches of rocks in the borehole wall of the reservoir, bed boundaries or a rock matrix.

12. The method of claim 11, wherein identifying the petrophysical facies includes simulating via dilation of other petrophysical facies of the collected at least one core by a fixed number of voxels.

13. The method of claim 1, wherein the at least one core is collected from at least one other reservoir.

14. The method of claim 1, wherein a second core from the at least one core is obtained from one or more reservoirs.

15. The method of claim 1, wherein step (d) includes plotting a digital file of the generated realizations of the numerical pseudocore model onto one of a digital media or hard copy media.

16. The method of claim 1, wherein the generating of realization of the numerical pseudocore model of step (d) includes using a multi-point statistics (MPS) algorithm to generate realizations of numerical pseudocores for two or more petrophysical facies, so as to condition the realizations to match facies sizes and shapes observed in the portion of the at least one interpretable borehole image of the portion of the generated digital core data and the generated at least one warped fullbore image data.

17. The method of claim 16, wherein the two or more petrophysical facies are from the group consisting of one of the collected at least one core, the portion of logging data or both, wherein the petrophysical facies is from the group consisting of one of fractures, pores, vugs, conductive patches of rocks in a borehole wall of the reservoir, resistive patches of rocks in the borehole wall of the reservoir, bed boundaries or a rock matrix.

18. The method of claim 15, further comprising:

e) resampling the numerical pseudocores of step (d) to a radial grid.

19. The method of claim 18, wherein the radial grid provides for a flow investigation of the numerical pseudocore model and includes:

1) building a cylindrical grid having one of one or more shape or one or more layer;

2) resampling the numerical pseudocore model from Cartesian coordinates, each voxel of which has a constant porosity or a constant permeability according to an associated rock type of the voxel, to cylindrical coordinates, wherein each cylindrical cell includes multiple Cartesian voxels of the numerical pseudocore model;

3) assigning each cylindrical cell an average porosity and an average permeability based on the included multiple Cartesian voxels; and

4) assigning each cylindrical cell relative permeability and capillary pressure curves based on the dominant rock type of the included multiple Cartesian voxels, using a reference table of one of capillary pressure, relative permeability or both for different facies in the numerical pseudocore that is fed into a flow simulator, whereby to quantify an impact of rock heterogeneity on fluid flow based on the pseudocore model referring to one of special core analysis (SCAL), mercury injection capillary pressure (MICP) data, pore network modeling or some combination thereof, of one or more core samples having the same rock type.

20. The method of claim 19, wherein the associated rock type is from the group consisting of one of fractures, pores, vugs, conductive patches of rocks in a borehole wall of the reservoir, resistive patches of rocks in the borehole wall of the reservoir, bed boundaries or a rock matrix.

21. The method of claim 18, further comprising:

f) performing flow simulations of one or more geographically associated reservoir such as a near-wellbore region.

22. The method of claim 21, wherein performing flow simulations includes one of a single phase or multiphase flow, such that the flow simulation is performed on the calculated numerical pseudocore model of step (e) so as to estimate one or more parameter.

23. The method of claim 22, wherein the one or more parameter is from the group consisting of porosity, permeability, capillary pressure, relative permeability, water cut, oil recovery factor or recovery efficiency.

24. The method of claim 21, further comprising:

g) conducting a sensitivity analysis on flow-related parameters, so as to determine a process for one of water flooding, gas flooding, steam flooding or enhanced oil recovery (EOR).

25. The method of claim 24, wherein the flow-related parameters is from the group consisting of porosity, permeability, capillary pressure, relative permeability, water cut, oil recovery factor or recovery efficiency.

26. The method of claim 21, further comprising:

g) determining a process for evaluating wettability effects in estimated parameters, wherein the estimated parameters includes one of capillary pressure, relative permeability, water cut, oil recovery factor or recovery efficiency.

27. A method for creating a numerical pseudocore model comprising:

a) obtaining logging data from a reservoir that includes depth-defined intervals of the reservoir, and processing the logging data into a portion of at least one interpretable borehole image data;

b) examining one of the portion of the at least one interpretable borehole image data, other processed logging data or both to generate at least one warped fullbore image data;

c) collecting one of at least one core from the reservoir, a portion of the logging data or both and generating a digital core data from one of the collected at least one core, the portion of logging data or both such that the generated digital core data represents one of features or structures of one or more depth-defined interval of the reservoir; and

d) processing a portion of the generated digital core data, the portion of at least one interpretable borehole image data, the portion of the logging data or some combination thereof, to generate realizations of the numerical pseudocore model.

28. The method of claim 27, wherein the portion of the at least one interpretable borehole image data is a training image that is two-dimensional (2D) scalar arrays of continuously variable numerical values.

29. The method of claim 27, wherein the other processed logging data is from the group consisting of one of logging-while-drilling data or wireline logging data, or some combination thereof.

30. The method of claim 27, wherein the logging data is from the group consisting of one of logging data having multiple depths of investigation or analog models of core data from one or more reservoirs or some combination thereof.

31. The method of claim 27, wherein the generating of the digital core data consists of identifying petrophysical facies from one of the collected at least one core or the portion of logging data, wherein the petrophysical facies is from the group consisting of one of fractures, pores, vugs, conductive patches of rocks in a borehole wall of the reservoir, resistive patches of rocks in the borehole wall of the reservoir, bed boundaries or a rock matrix.

32. The method of claim 31, wherein identifying the petrophysical facies includes simulating via dilation of other petrophysical facies of the collected at least one core by a fixed number of voxels.

33. A system for creating a numerical pseudocore model comprising:

a) obtaining logging data from a reservoir that includes depth-defined intervals of the reservoir, and processing the logging data into a portion of at least one interpretable borehole image data having unidentified borehole image data;

b) examining one of the portion of the at least one interpretable borehole image data, other processed logging data or both to generate the unidentified borehole image data, then processing the generated unidentified borehole image data into the portion of the at least one interpretable borehole image data so as to generate at least one warped fullbore image data;

c) collecting one of at least one core from the reservoir, a portion of the logging data or both and generating a digital core data from one of the collected at least one core, the portion of logging data or both such that the generated digital core data represents one of features or structures of one or more depth-defined interval of the reservoir; and

d) processing a portion of the generated digital core data, the portion of at least one interpretable borehole image data, portion of the logging data or some combination thereof, to generate realizations of the numerical pseudocore model;

wherein the generating of realization of the numerical pseudocore model includes using a multi-point statistics (MPS) algorithm to generate realizations of numerical pseudocores for two or more facies, so as to condition the realizations to match facies sizes and shapes observed in the portion of the at least one interpretable borehole image of the portion of the generated digital core data and the generated at least one warped fullbore image data.

34. A system for creating a numerical pseudocore model comprising:

a) obtaining logging data from a reservoir that includes depth-defined intervals of the reservoir, and processing the logging data into a portion of at least one interpretable borehole image data;

b) examining one of the portion of the at least one interpretable borehole image data, other processed logging data or both to generate at least one warped fullbore image data;

c) collecting one of at least one core from the reservoir, a portion of the logging data or both and generating a digital core data from one of the collected at least one core, the portion of logging data or both such that the generated digital core data represents one of features or structures of one or more depth-defined interval of the reservoir; and

d) processing a portion of the generated digital core data, the portion of at least one interpretable borehole image data, the portion of the logging data or some combination thereof, to generate realizations of the numerical pseudocore model;

wherein the generating of the digital core data consists of identifying facies from one of the collected at least one core or the portion of logging data, wherein the facies is from the group consisting of one of fractures, pores, vugs, conductive patches of rocks in a borehole wall of the reservoir, resistive patches of rocks in the borehole wall of the reservoir, bed boundaries or a rock matrix.