WO2009126881A2 - Procédé pour produire des pseudo-carottes numériques utilisant des images de fond de puits, des échantillons de roche numériques, et des statistiques multi-emplacements - Google Patents
Procédé pour produire des pseudo-carottes numériques utilisant des images de fond de puits, des échantillons de roche numériques, et des statistiques multi-emplacements Download PDFInfo
- Publication number
- WO2009126881A2 WO2009126881A2 PCT/US2009/040198 US2009040198W WO2009126881A2 WO 2009126881 A2 WO2009126881 A2 WO 2009126881A2 US 2009040198 W US2009040198 W US 2009040198W WO 2009126881 A2 WO2009126881 A2 WO 2009126881A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- data
- reservoir
- image data
- borehole
- core
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000011435 rock Substances 0.000 title claims description 109
- 238000012545 processing Methods 0.000 claims abstract description 17
- 238000012549 training Methods 0.000 claims description 90
- 208000035126 Facies Diseases 0.000 claims description 50
- 239000011148 porous material Substances 0.000 claims description 45
- 238000004088 simulation Methods 0.000 claims description 41
- 238000004422 calculation algorithm Methods 0.000 claims description 33
- 230000035699 permeability Effects 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000005553 drilling Methods 0.000 claims description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 238000011084 recovery Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 7
- 238000003491 array Methods 0.000 claims description 6
- 230000010339 dilation Effects 0.000 claims description 5
- 238000011835 investigation Methods 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 230000001143 conditioned effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 238000010795 Steam Flooding Methods 0.000 claims 1
- 230000005251 gamma ray Effects 0.000 claims 1
- 238000002347 injection Methods 0.000 claims 1
- 239000007924 injection Substances 0.000 claims 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims 1
- 229910052753 mercury Inorganic materials 0.000 claims 1
- 238000010206 sensitivity analysis Methods 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 description 23
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 17
- 238000003384 imaging method Methods 0.000 description 13
- 239000003086 colorant Substances 0.000 description 8
- 238000013459 approach Methods 0.000 description 7
- 238000002591 computed tomography Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 238000010603 microCT Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000010354 integration Effects 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 235000015076 Shorea robusta Nutrition 0.000 description 3
- 244000166071 Shorea robusta Species 0.000 description 3
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 3
- 238000004624 confocal microscopy Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- CJDRUOGAGYHKKD-XMTJACRCSA-N (+)-Ajmaline Natural products O[C@H]1[C@@H](CC)[C@@H]2[C@@H]3[C@H](O)[C@@]45[C@@H](N(C)c6c4cccc6)[C@@H](N1[C@H]3C5)C2 CJDRUOGAGYHKKD-XMTJACRCSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 230000003750 conditioning effect Effects 0.000 description 2
- 238000013079 data visualisation Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000010884 ion-beam technique Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000001000 micrograph Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 235000021419 vinegar Nutrition 0.000 description 2
- 239000000052 vinegar Substances 0.000 description 2
- 101000965313 Legionella pneumophila subsp. pneumophila (strain Philadelphia 1 / ATCC 33152 / DSM 7513) Aconitate hydratase A Proteins 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000013050 geological sediment Substances 0.000 description 1
- 125000001475 halogen functional group Chemical group 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000001878 scanning electron micrograph Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000009718 spray deposition Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6163—Electromagnetic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6167—Nuclear
Definitions
- the invention is generally related to a method using a seminal Multi-point statistics (MPS) algorithm to generate numerical pseudocores from digital rock or core samples and borehole-imaging logs. More particularly, this patent specification relates to creating 3D numerical cores from computed X-ray tomography (CT scans) and formation micro-image (FMI) logs, and performing flow modeling in these numerical cores to understand fluid-flow paths and recovery factors in selected reservoir.
- MPS seminal Multi-point statistics
- Electrical and acoustic borehole-imaging tools are widely used to log subsurface boreholes to locate and map the boundaries between rock layers, e.g., bed boundaries, and to visualize and orient fractures and faults. Because electrical logging tools are pad-type devices with fixed arrays of electrodes, it is common to have gaps with missing information between the pads. Electrical and acoustic logs commonly have intervals with poor data quality due to non-functioning electrodes, insufficient pad pressure, borehole irregularities, rock debris, decentralized tools, or poor acoustic reflections.
- CTscans are the most widely used approach.
- CTscans are 2-dimensional (2D) cross sections generated by an X-ray source that rotates around the sample. Density is computed from X-ray attenuation coefficients. Scans of serial cross sections are used to construct 3D images of the sample. Because the density contrast is high between rock and fluid-filled pores, CT images can be used to visualize the rock-pore system. Resolutions are on the sub-millimeter to micron scale, depending on the device being used.
- Multi-point statistics are used to create simulations of spatial geological and reservoir property fields for subsurface reservoir modeling. These methods are conditional simulations that use known results, such as those measured in wellbores, as fixed or "hard” data that are absolutely honored during the simulations.
- MPS uses ID, 2D, or 3D "training images” as quantitative templates to model subsurface property fields.
- MPS modeling captures geological structures from training images and anchors them to data locations. These structures can be either a priori geological interpretations or conceptual models.
- Multipoint geostatistics is a new advanced geostatistics approach. It allows reservoir modelers to incorporate their prior knowledge, interpretations, or conceptual models into the reservoir modeling process through training images. These training images are numerical representations of the structures/features that are believed to exist in the reservoir under study. Once we have the training images, MPS can extract curvilinear structures or complex features from the training images and anchor them to the reservoir locations where the samples/observations are collected, leading to more realistic reservoir models. Introducing training images into reservoir modeling is a milestone. Note that there are two ingredients in the use of MPS: training images (conceptual models) and the real data (observations). These two pieces are typically separated. However, in realistic applications, generating representative training images, in particular in 3D, has proved to be a bottleneck in MPS applications. Generating a continuous variable training image is even more difficult than the creation of categorical training image.
- electrical borehole images may run in water-based (conductive) mud, such as Schlumberger's FMI (Formation Microlmager) log, which is based on dipmeter technology that has been commercially available since the 1950's.
- Electrical borehole- imaging tools are, in essence, sophisticated dipmeters.
- the imaging tools have microresistivity electrodes arranged around the wellbore on pads that are pressed against the borehole wall.
- the evolutionary trend from dipmeters to borehole images has been from a few electrodes to a complex array of electrodes on multiple pads. See Hurley, N. F., 2004, Borehole Images, in Asquith, G.
- Tools are first run into the hole with the pads closed. At the start of the log run, either four, six, or eight pads are pressed against the borehole wall. The number of pads depends on the logging device. Electrical current is forced into the rock through the electrodes, and sensors measure the current after it interacts with the formation. Raw data include multiple electrode readings, caliper readings from individual pads or pairs of pads, and x-, y-, and z-axis accelerometer and magnetometer readings. Borehole deviation and pad 1 (tool) orientation are determined from the magnetometers. The sample rate for electrode and accelerometer data is very high, normally 120 samples/ft (400 samples/m).
- Areal coverage of the borehole face is a function of width of the electrode arrays, number of pads, and borehole diameter. In general, 40 to 80% of the borehole face is imaged in typical boreholes. Non-imaged parts of the borehole appear as blank strips between the pads.
- Borehole images are created by assigning color maps to different bins or ranges of resistivity values. Colored pixels are then arranged in their proper geometric position around the wellbore. By convention, low-resistivity features, such as shales or fluid-filled fractures, are displayed as dark colors. High-resistivity features, such as sandstones and limestones, are displayed as shades of brown, yellow, and white.
- Static images are those which have had one contrast setting applied to the entire well. They provide useful views of relative changes in rock resistivity throughout the borehole.
- Dynamic images which have had variable contrast applied in a moving window, provide enhanced views of features such as vugs, fractures, and bed boundaries. Dynamic images bring out subtle features in rocks that have very low resistivities, such as shales, and very high resistivities, such as carbonates and crystalline rocks.
- Oil-Based (non- conductive) Mud in particular high mud resistivities (greater than 50 ohm-m), typical of oil- based muds, are unsuitable for most electrical borehole images.
- OBMI Oil-Base Microlmager
- This tool generates borehole images by passing electrical current into the formation from two large electrodes on each pad, which is at a high voltage (about 300V).
- GVR GeoVision Resistivity
- ADN Azimuthal Density Neutron
- a aspect of Borehole images can be Acoustic borehole images, also known as borehole televiewers, which are based on technology first developed in the 1960's. Zemanek, J., Glenn, E. E., Norton, L. J., and Caldwell, R. L., 1970, Formation evaluation by inspection with the borehole televiewer: Geophysics, v. 35, p. 254-269.
- the Ultrasonic Borehole Imager (UBI) is Schlumberger' s primary acoustic tool for open-hole applications.
- the UBI tool which is centralized in the well, has a rotating transducer that emits and records sound waves that bounce off of the borehole wall. Both acoustic amplitude and travel time are recorded and processed into images. Normally, borehole coverage is 100%, with no gaps in the images. However, poor-quality images may result when the tool is decentralized, or the borehole wall is irregular.
- Petrophysical Facies may be considered, among other things, as characteristic signatures on borehole-image logs, such as vugs, and resistive and conductive patches.
- a particular view by Dehghani et al. in 1999 suggested that zones of enhanced porosity and permeability exist in the vicinity of vugs. Dehghani, K., Harris, P. M., Edwards, K. A., and Dees, W. T., 1999, Modeling a vuggy carbonate reservoir: AAPG Bulletin, v. 83, p. 19-42.
- Dehghani et al. (1999) used thin sections, SEM images, and mini-permeability measurements to confirm their concept.
- Schindler (2005) and Tanprasat (2005) used image analysis of fluorescent-inked core photos to show that swarms of small vugs preferentially exist in the vicinity of large vugs. See Schindler, J., 2005, Quantification of vuggy porosity, Indian Basin field, New Mexico: Unpublished M.S. thesis, Colorado School of Mines, Golden, CO.; and Tanprasat, S., 2005, Petrophysical analysis of vuggy porosity in the Shu'aiba Formation of the United Arab Emirates: Unpublished M.S. thesis, Colorado School of Mines, Golden, CO. Such small vugs are below the resolution of the borehole-imaging tool, so they appear as dark regions, rather than as discrete vugs in the image logs.
- the published literature has many examples of numerical rocks built using techniques (or digital rock models of rocks and pores) that include reconstructions made from 2D thin sections or scanning-electron microscope (SEM) images, electrofacies interpreted from logs, computer-generated sphere packs, laser-scanning confocal microscopy, and various types of CTscans (conventional, microCT, and synchrotron-computed micro tomography) .
- techniques or digital rock models of rocks and pores
- SEM scanning-electron microscope
- Bakke and Oren (1997), Oren et al. (1998), and Oren and Bakke (2002) developed a technique that constructs 3D pore networks from 2D thin sections.
- Numerical Rocks (http://www.numericalrocks.com/) computes 3D pore models from 2D thin sections. See Bakke, S., and Oren, P.-E., 1997, 3-D pore-scale modeling of sandstones and flow simulations in the pore networks: SPE preprint 35,479, European 3-D Reservoir Modeling Conference, Stavanger, Norway, April 16-17, p. 136-149; Oren, P.-E., Bakke, S., and Arntzen, O.
- Dvorkin et al. (2003) described Digital Rock Physics technology, which consists of pore-scale numerical simulations derived from: (a) 2D thin sections and statistical indicator simulation, or (b) CTscans. See Dvorkin, J., Kameda, A., Nur, A., Mese, A., and Tutuncu, A. N., 2003, Real time monitoring of permeability, elastic moduli and strength in sands and shales using Digital Rock Physics: SPE preprint 82246, presented at the SPE European Formation Damage Conference, The Hague, Netherlands, May 13-14, 7 p. They built 3D models of virtual rock, and did flow simulations using the lattice-Boltzmann method. U.S. Patent 6,516,080 (below) is related to this work.
- Mini-models are populated using "principle rock types” (PRT), which "cover and categorize the full range of pore types, sizes, pore- throat size distributions, capillary entry pressures, relative permeabilities, etc.”
- PRT' s are organized into “rock type associations” (RTA), which are based on "sedimentary fabric” determined from borehole-image logs. RTA' s are distributed in the reservoir using borehole- image logs, and observed layering, facies models, and seismic data.
- Bosl et al. (1998) and Holt (2001) have generated similar digital rock models for flow experiments. See Bosl, W. J, Dvorkin, J., and Nur, A., 1998, A study of porosity and permeability using a lattice Boltzmann simulation: Geophysical Research Letters, v. 25, p. 1475-1478; and see Holt, R. M., 2001, Particle vs. laboratory modelling in in situ compaction: Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy, v. 26, Issue 1-2, p. 89-93.
- Multi-point (or multiple-point) statistical methods are a new family of spatial statistical interpolation algorithms proposed in the 1990's that are used to generate conditional simulations of discrete variable fields, such as geological facies. See Guardiano, F. and Srivastava, R. M., 1993, Multivariate geostatistics: beyond bivariate moments: Geostatistics-Troia, A. Soares. Dordrecht, Netherlands, Kluwer Academic Publications, v. 1, p. 133- 144.
- a training image is a numerical prior geological model that contains the facies structures and relationships believed to exist in realistic reservoirs. Training images are conceptual in nature and can be as simple as a hand-drawn map, or they can be created by computer tools.
- Strebelle's (2002) seminal MPS algorithm, called SNESIM has been used in many applications for reservoir modeling, and has become the reference tool for modeling fluvial channel deposits when combined with rotation and affinity transformations (Zhang, 2002; Caers and Zhang, 2004; Strebelle and Zhang, 2004). See Zhang, T., 2002, Multiple-point simulation of multiple reservoir facies: Unpublished M.S. thesis, Stanford University, California, 163 p.; see Caers, J. and Zhang, T., 2004, Multiple- point geostatistics: A quantitative vehicle for integration of geologic analogs into multiple reservoir models, in M. Grammer, P. M. Harris and G. P.
- an MPS algorithm such as FILTERSIM (Zhang 2006a). See Zhang, T., 2006a, Filter-based training image pattern classification for spatial pattern simulation: Unpublished Ph.D. thesis, Stanford University, California, 153 p.
- the FILTERSIM algorithm applies a set of local filters to the training image, which can be either categorical or continuous, to group local patterns into pattern classes. It then proceeds to simulate patterns on the basis of that classification.
- a filter is a local template (window) with a set of weights associated to each pixel location of the template.
- the RAM limitation also prevents us from considering too many categories or classes jointly, thus limiting seminal MPS algorithm to the simulation of categorical variables.
- the seminal MPS algorithm searches for exact replicates of the conditioning data event, builds the reservoir model one pixel at a time, conditioned to a multiple-point data event, and does not allow any filtering or averaging of the patterns found in the training image.
- FILTERSIM a new MPS algorithm named FILTERSIM was proposed by Zhang (2006a).
- the HLTERSIM algorithm applies a set of local filters to the training image, which can be either categorical or continuous, to group local patterns into pattern classes. It then proceeds to simulate patterns on the basis of that classification.
- a filter is a local template (window) with a set of weights associated to each pixel location of the template. Applying a filter to a local pattern results in a filter score, the score is viewed as a numerical summary of that local pattern.
- a set of default or use- defined filters is designed such that each filter can record different aspects of the training pattern seen within the template. These filters are used to transform training patterns into a filter score space. This pattern scoring provides a dimension reduction of patterns. By partitioning that score space of limited dimension, similar training patterns are classified based on their filter scores.
- the FILTERSIM algorithm starts with a classification of local training patterns in a filter score space of reduced dimension. Simulation proceeds along a sequential path through the simulation space, by determining which pattern class is most similar to the local conditioning data event, sampling a specific pattern from the pattern class, and then patching the sampled pattern onto the image at the simulation sites.
- the simulation random path and the sampling of patterns from pattern classes allow for different simulated realizations, yet all are conditional to the same original data. Because of the dimension reduction brought by the filter summaries of any pattern, and because patterns are grouped into classes, the algorithm is fast and reasonable in terms of RAM demand.
- the seminal MPS algorithm and FILTERSIM algorithm are able to honor absolute or so-called "hard” constraints from data acquired in wells or outcrops, and conditional or “soft” constraints from seismic data, facies probability fields, and rotation and affinity (or scale) constraint grids. All of these data are used in the stochastic modeling process to generate ID, 2D, or 3D maps of geological facies or rock properties. Because there is a random component involved in MPS simulations, individual realizations of property fields created by MPS algorithms differ, but the ensemble of realizations provides geoscientists and reservoir engineers with improved quantitative estimates of the spatial distribution and uncertainty of geological facies in a modeled reservoir volume. Moreover, these algorithms honor both hard and soft input data constraints (Zhang, 2006a). See Zhang, T., Switzer P., and Journel A., 2006b, Filter-based classification of training image patterns for spatial pattern simulation: Mathematical Geology, v. 38, p. 63-80.
- rescaled filters are applied over the respective grid (Zhang, 2006a).
- training images There are two types of training images: one with a very limited number of categories and another for continuous variables such as reservoir petrophysical properties.
- Multipoint geostatistical methods require ID, 2D, or 3D grids of training images as prior conceptual geological models that contain patterns of the spatial attributes under study.
- the shapes of different features appearing on the images are supposed to represent a model of real geological features, with each category typically representing a different geological facies or different kind of geological body.
- Training images are typically required to contain "stationary" patterns, i.e., the patterns must be independent of their location in space (invariant according to any translation) and must be repetitive over the training image area.
- this stationarity can consist, but is not limited to, geological object orientation stationarity (where directional objects/features do not rotate across the image) and geological scale stationarity (where the size of objects/features on the image does not change across the image) (Caers and Zhang, 2004). See Caers, J. and Zhang, T., 2004, Multiple-point geostatistics: A quantitative vehicle for integration of geologic analogs into multiple reservoir models, in M. Grammer, P. M. Harris and G. P. Eberli, eds.: Integration of Outcrop and Modern Analogs in Reservoir Modeling, AAPG. Memoir 80, p. 383-394.
- Training images are supposed to model or reproduce real geological features and should as much as possible be derived from existing geologically meaningful images. This requires research on statistical and image-processing methods that will allow use of images from any source, e.g., hand-drawn sketches, aerial photographs, satellite images, seismic volumes, geological object models, physical scale models, or forward geological process models. Compared to the creation of continuously variable training images, generating categorically variable training images is easier.
- An object-based approach is commonly used to generate training images with categorical variables.
- a region-based approach, combined with adding desired constraints, can be used to generate continuously variable training mages (Zhang et al., 2005).
- the present invention relates to a method to generate 3-dimensional (3D) computer models of rocks and pores, known as numerical pseudocores.
- the technique uses fullbore images, digital rock samples, and multi-point statistics (MPS) to reconstruct a 3D pseudocore for a logging interval where there is no real core collected, but there are logged borehole images.
- the digital core samples are used to guide the 3D feature reconstruction of the pseudocores by multiple-point geostatistics and the final pseudocores are constrained by the fullbore images.
- FIG. Ia shows a training image defined as 3 foot ( or 1 meter) interval of a borehole-image log, in particular, of a single-pass formation micro-imager (FMI) in a vuggy carbonate formation, according to at least one embodiment of the invention
- FIG. Ib illustrates the result of a Multi-Point Simulation (MPS) using FILTERISM of the same training image of FIG. Ia, according to at least one embodiment of the invention
- FIG. 2a shows the same training image of FIG. Ia, according to at least one embodiment of the invention
- FIG. 2b illustrates the result of a first realization of the training image of FIG. Ia and matches the original, measured data of the vuggy formation of Figure 1, such that the measured date are honored in the first (FIG. 2a), second (FIG. 2b) and third (FIG. 2c) realizations which provide subtle variations in the modeled areas, according to at least one embodiment of the invention;
- FIG. 2c illustrates the result of a second realization of the training image of FIG. Ia and matches the original, measured data of the vuggy formation of FIG. Ia, according to at least one embodiment of the invention
- FIG. 2d illustrates the result of a third realization of the training image of FIG. Ia and matches the original, measured data of the vuggy formation of FIG. Ia, according to at least one embodiment of the invention
- FIG. 3a shows a static image of vugs appearing as dark features in the logging- while drilling (LWD) images of a vuggy carbonate, such that conductive patches are illustrated in brown and non-conductive patches illustrated in white, according to at least one embodiment of the invention
- FIG. 3b shows a dynamic image of vugs appearing as dark features in the logging- while drilling (LWD) images of a vuggy carbonate, such that conductive patches are illustrated in brown and non-conductive patches illustrated in white, according to at least one embodiment of the invention;
- FIG. 4 shows a Fullbore image of vuggy porosity of FIGIa that shows contours (green lines) that outline the less-resistive areas of the electrical image, according to at least one embodiment of the invention;
- FIG. 5a and 5b show a block diagrammatical view of the method to generate numerical pseudocores using borehole images, digital rock samples and multipoint statistics, according to at least one embodiment of the invention
- FIG. 6 shows the training image of FIGIa in a translucent, 3D view of a CT scan of slabbed vuggy carbonate sample (approximate 6 inches in height and 4 inches in diameter), such that the vugs are illustrated in red to gray and the rock matrix illustrated in blue, according to an embodiment of the invention
- FIG. 7 shows borehole images warped to cylindrical shape, matching a borehole diameter (approximate borehole diameter 8.5 inches and approximate length in interval 3 feet or 1 meter), at a particular depth.
- CT scan training image of FIG 6 is correctly scaled and poisoned at its correct depth in the center of the borehole images; wherein the vugs are illustrated in red and the rock matrix illustrated in light blue, according to embodiments of the invention;
- FIG. 8 illustrates the result of a numerical two-facies pseudocore generated from the training image and fullbore image that corresponds to FIG 7, such that the vugs are illustrated in red and the rock matrix illustrated in dark blue, according to embodiments of the invention
- FIG. 9 illustrates the result of a numerical pseudocore generated from the training image and fullbore image that corresponds to FIG 7, such that the vugs are invisible and the rock matrix illustrated in gray.
- the arbitrary suspended slice shows that this is a 3D model, with an abundant of pores in every slice (approximate length of interval is 1 foot or 0.3 meters and an approximate borehole diameter of 8.5 inches), according to embodiments of the invention;
- FIG. 10 illustrates the result of a numerical pseudocore generated from the training image and fullbore image that corresponds to FIG 6.
- Conductive patches (illustrated in red) resemble contours around conductive patches as shown in FIG. 4. It is noted three petrophysical facies are shown: pores or vugs (illustrated in green), conductive patches (illustrated in red), and rock matrix (illustrated in blue), according to embodiments of the invention;
- FIG. 11 shows the numerical pseudocore as of FIG. 10, with an invisible rock matrix, which highlights the 3D interconnected nature of the pores or vugs (illustrated in green) and conductive patches (illustrated in red).
- FIG. 12 shows numerical pseudocore as Figure 10, and the complex 3D shape of the pores or vugs (illustrated green), and both the conductive patches and rock matrix are invisible, according to embodiments of the invention;
- FIG. 13 shows numerical pseudocores that can be resampled or regrided to radial grids wherein the radial grids can be layered, based upon layers observed in borehole images or other well logs, according to embodiments of the invention
- FIG. 14a shows relative permeability curve that is used for conductive patch in the flow simulation; wherein there are three different relative permeability curves for each of the three rock types in the pseudocore model: such as, the vugs, conductive patch and tight rock matrix, according to embodiments of the invention.
- FIG. 14b shows flow simulation result through a numerical pseudocore, such that a line of micro-injectors of water surrounds outer the diameter, a line of micro- producers surrounds the inner diameter of the pseudocore: the colors represent oil saturation (So), where heterogeneous fingering and breakthrough is shown in this flow model, according to embodiments of the invention.
- So oil saturation
- the present invention relates to a method to generate 3-dimensional (3D) computer models of rocks and pores, known as numerical pseudocores.
- the technique uses fullbore images, digital rock samples, and multi-point statistics (MPS) to reconstruct a 3D pseudocore for a logged interval where there is no real core collected, but there are borehole images.
- the digital rock samples are used to guide the 3D feature reconstruction of the pseudocores by multiple-point geostatistics and the final pseudocores are constrained by the fullbore images.
- the seminal idea lies in the use of training image: we directly use the data themselves [in fullbore creation, it is the original 2D incomplete image (continuous variable training image) that has >60% coverage of the entire region, while in the pseudocore reconstruction, the Catscan digital core is directly used as a 3D training image that is combined later with the fullbore image data]. Hence, the entire process of applying MPS becomes data-driven. This advantage should be stressed in our patent memo and provisions
- Borehole images produce oriented electrical and acoustic maps of the rocks and fluids encountered by a borehole.
- Fullbore images are complete, 360-degree views of the borehole wall. By design, most acoustic and logging-while-drilling tools generate fullbore images. Most resistivity image logs need to have fullbore images modeled because gaps exist between pads, and there may be damage to certain electrodes or pads. The modeling process uses continuous-variable algorithms developed within MPS.
- Fullbore images provide the physical location of features, such as bed boundaries, pores, and conductive and resistive patches within the rock on the borehole wall. These outline complex 3D volumes, known as petrophysical facies.
- Digital rocks or core samples can be generated from 2-dimensional thin sections, scanning-electron-microscope images, confocal-microscope images, or computer- generated sphere packs. Most digital rocks or core samples, however, are generated from computed- tomographic scans (CTscans) of rocks and fluids. Tomographic images are created when the scanner transmits X-rays at different angles through the rock to a receiver. X-ray attenuation is converted to density, and serial sections are built into 3D views of the rocks and pores.
- CTscans computed- tomographic scans
- each digital core sample for example, a CTscan, is directly taken as a training image. These are discrete variable training images with the attribute being the rock (white) or pore (black) at each pixel of the image.
- the training image can have any shape of boundaries or contain any number of irregular holes.
- Numerical pseudocores are created using discrete-variable algorithms within MPS. Integer values are assigned to each petrophysical facies, such as rock (0), pores (1), and conductive patches (2). Digital rock or core samples are used as training images, i.e., the quantitative templates used to model property fields. Fullbore images surround the numerical pseudocore with cylindrical envelopes to condition the models. Each numerical pseudocore absolutely honors the digital rock or core samples and fullbore images. Numerical pseudocores can be gridded into models suitable for fluid-flow simulations. Capillary pressure and relative permeability curves are provided by conceptual models, special core analysis, or established techniques of fine-scale pore-network modeling.
- FIG. Ia is a training image defined as a 3 foot (1 meter) interval of a borehole- image log, of a single-pass formation micro-imager (FMI) in a vuggy carbonate formation.
- Vugs are pores that are filled with water-based drilling mud, and they appear as dark, low- resistivity spots. Note the gaps between the pads in FIG. Ia.
- the pixel-based, user-defined 3x3 template (bottom) is moved through the training image, detecting patterns and giving filter scores to the neighborhoods around each measured pixel. This provides the basis for MPS simulation using FILTERSIM.
- the Fullbore image in FIG. Ib shows a FILTERSIM realization that uses the entire image of FIG. Ia as a training image.
- FIG. Ib illustrates the result of a Multi-Point statistics (MPS) using FILTERISM of the same training image of FIG. Ia.
- MPS Multi-Point statistics
- FIG. 2b, 2c and 2d show three realizations of a fullbore image in a vuggy carbonate formation, wherein FIG. 2a shows the same training image of FIG. Ia. Note that the measured data are honored in each realization, and there are subtle variations in the modeled areas.
- FIG. 2b illustrates the result of a first realization of the training image of FIG. 2a and matches the original, measured data of the vuggy formation of FIG. Ia.
- FIG. 2c illustrates the result of a second realization of the training image of FIG. 2a and matches the original, measured data of the vuggy formation of FIG. Ia.
- FIG. 3d illustrates the result of a third realization of the training image of FIG. 2a and matches the original, measured data of the vuggy formation of FIG. Ia.
- FIG. 3a shows a static image of vugs appearing as dark features in the logging-while drilling (LWD) images of a vuggy carbonate (see Xiao, L., Jun, C, Duo, Y. S., Han, S. Y., Xia, W. H., and Xi, W. Y., 2007: Fully integrated solution for LWD resistivity image application a case study from Beibu Gulf, China: 1 st SPWLA India Regional Conference, Formation Evaluation in Horizontal Wells, Mumbai, March 19-20, 10 p.), such that conductive patches are illustrated in brown and non-conductive patches illustrated in white.
- 3b shows a dynamic image of vugs appearing as dark features in the logging-while drilling (LWD) images of a vuggy carbonate, such that conductive patches are illustrated in brown and non-conductive patches illustrated in white.
- LWD logging-while drilling
- FIG. 5 shows a block diagrammatical view of the method to generate numerical pseudocores using borehole images, digital rock samples and multi-point statistics.
- Step 100 includes collecting depth-defined intervals of borehole-imaging logs from a reservoir that can be used as training images.
- the training images are oriented; 2-dimentional (2D) scalar arrays of continuously variable numerical values, with gaps between the pads and/or areas that need repair, such that the pads represent measured values and the gaps are non-image parts of the borehole.
- 2D 2-dimentional
- the original Fullbore images are generated by certain logging tools, such as acoustic devices and logging-while-drilling tools. In other cases, fullbore images must be generated. This is especially true for resistivity logs, which commonly have gaps between the pads.
- the interval could be 1, 3, or 10 ft (0.3, 1, or 3 m) of measured depth.
- the user may want to choose a thick or thin interval, depending on the observed amount of layering, fracturing, and other heterogeneous patterns.
- Fullbore images can be generated using the FILTERSIM algorithm within MPS as noted in Step 100 of FIG. 5.
- FILTERSIM uses filter scores to group and then simulate patterns in the gaps between the pads, where no measured data exist.
- Step 200 of FIG. 5 includes collecting at least one core sample from the reservoir and then generating a digital core sample from the collected core sample.
- the collected core sample similarly portrays the representative features and/or structure of one or more depth- defined interval of the borehole and/or of the reservoir. For example, digital rocks or core samples, such as CTscans, can be used as training images.
- Training images are 3D arrays of discrete numerical values.
- the rock has a numerical value of 0, and the pores have a numerical value of 1.
- the pores have a numerical value of 1
- the conductive patches have a numerical value of 2.
- Outlines of individual facies bodies (volumes) can have any shape or size.
- Step 300 of FIG. 5 includes pre-modeling of the collected borehole-imaging log Data such as generating fullbore images.
- Step 300a of FIG. 5 discloses making a selection of a depth-defined interval of the borehole-image log.
- Fullbore images are generated using a FILTERSIM algorithm within Multi-Point Statistics (MPS), wherein MPS modeling captures geological structures from training images and anchors them to data locations.
- MPS Multi-Point Statistics
- FILTERSIM uses filter scores to group and then simulate patterns in the gaps between the pads, where no measured data exists.
- Step 300b of FIG. 5 discloses truncating into petrophysical facies that identify, for example: pores or vugs; conductive and resistive patches of the rock on the borehole wall; and bed boundaries or the rock matrix.
- Step 300c of FIG. 5 discloses warping the Fullbore Images into Scaled Cylindrical Shapes. In particular, for routine interpretations, it is difficult to examine borehole images in 3D. Therefore, it is common to split the borehole along true north, and then unroll the cylinder until it becomes a 2D view.
- Step 400 of FIG. 5 discloses pre-modeling of the collected digital core sample.
- Step 400a converts from 2 facies to 3 facies so that the pores or vugs, conductive patches and rock matrix in the digital rock have been identified.
- the conductive patches have not been identified in the digital rock (training image)
- they can be simulated by dilation of the porous facies by a fixed number of voxels.
- large-scale heterogeneity can be captured in the numerical pseudocore.
- the voxel resolution of the digital rock or core training image is finer than the fullbore image resolution.
- the training image is coarsely sampled according to the ratio of the resolutions of the digital core to the fullbore image.
- Step 500 of FIG. 5 discloses generating realizations of numerical pseudocores by merging the structures borrowed from the collected digital core data and collected borehole- imaging log data.
- Step 500a discloses using the MPS SNESIN algorithm to generate realizations of numerical pseudocores for 2 or more facies, so as to condition the realizations to match the facies sizes and shapes observed in training images of the digital rock and the fullbore images. More than 2 facies can be modeled in cases where the conductive patches have been mapped ( Figures 3 and 4) in the fullbore images and the digital rocks. If the conductive patches have not been identified in the digital rock (training image), they can be simulated by dilation of the porous facies by a fixed number of voxels. In this way, large- scale heterogeneity can be captured in the numerical pseudocore.
- the radial size and height of the numerical pseudocore is limited only by the amount of computer memory that is available.
- Step 600 of FIG. 5 discloses resampling numerical pseudocores to a radial grid.
- the Cartisian numerical pseudocore model generated using the previous steps into a radial grid in the form of cylindrical coordinates.
- matrix, vug, or conductive patch For Cartesian each cell of the radial grid, because it consists of many Cartesian voxels, averaged porosities and permeabilities are generated.
- the averaged porosity is obtained by arithmetically averaging all porosity of the Cartesian voxels within the cell; the permeability is obtained by performing geometric average.
- Step 700 of FIG. 5 discloses perform flow simulations of the near-wellbore region.
- Numerical simulations of fluid flow e.g. water flooding
- a look-up table of capillary pressure and relative permeability for different facies in the numerical pseudocore provides values that are fed into a flow simulator, for example Ecllipse. This is the key step to quantify the impact of carbonate rock heterogeneity on fluid flow based on the pseudocore model.
- the capillary pressure and relative permeability could be obtained from SCAL or MICP data (if available) of core samples with the same rock type.
- FIG. 6 shows the training image of FIGIa in a translucent, 3D view of a CTscan of slabbed vuggy carbonate sample (approximate 6 inches in height and 4 inches in diameter), such that the vugs are illustrated in red to gray and the rock matrix illustrated in blue.
- FIG. 6 also shows the chosen training image, a CTscan of a vuggy carbonate (Gowelly, S., 2003, 3-D analysis of vug connectivity, Indian Basin field, New Mexico: Unpublished M.S.
- Figure Ia shows a 3-ft (1-m) interval of a borehole-image log, from the same depth and in the same formation, and that vugs are large, irregular pores, visible to the naked eye. Further still, in a borehole-image log, vugs can appear as dark spots because they are filled with water-based drilling mud, and they conduct electricity.
- FIG. 7 shows borehole images warped to cylindrical shape, matching a borehole diameter (approximate borehole diameter 8.5 inches and approximate length in interval 3 feet or 1 meter), at a particular depth.
- CT scan training image of FIG 6 is correctly scaled and poisoned at its correct depth in the center of the borehole images, wherein the vugs are illustrated in red and the rock matrix illustrated in light blue.
- the borehole images are warped to the 3D cylindrical shape, such that the fullbore images are not shown, because it allows a view at the relative scales of the training image (CTscan) (in the center of the cylinder) and the log image.
- CTscan training image
- FIG.8 illustrates the result of a numerical pseudocore generated from the training image and fullbore image that corresponds to FIG 7, such that the vugs are illustrated in red and the rock matrix illustrated in dark blue. It is noted that the length of the interval can be 3 feet (1 meter), and the borehole diameter can be 8.5 inches (22 cm).
- FIG. 9 illustrates the result of a numerical pseudocore generated from the training image and fullbore image that corresponds to FIG 7, such that the vugs are invisible and the rock matrix illustrated in gray. Further, the arbitrary suspended slice (see top of FIG. 9) shows that this is a 3D model, with an abundant of pores in every slice (approximate length of interval is 1 foot or 0.3 meters and an approximate borehole diameter of 8.5 inches).
- FIG. 10 illustrates the result of a numerical pseudocore generated from the training image and fullbore image that corresponds to FIG 6.
- the numerical pseudocore shows using a circular dilation with an 8-pixel radius around each pore. Conductive patches (illustrated in red) resemble contours around conductive patches as shown in FIG.
- This model therefore, has 3 petrophysical facies: pores or vugs (illustrated in green), conductive patches (illustrated in red), and rock matrix (illustrated in blue).
- the conductive patches provide 3D connectivity between the pores or vugs, and allow capturing the heterogeneity that is inherent in most carbonate rocks.
- FIG. 11 shows the numerical pseudocore as of FIG. 10, with an invisible rock matrix, which highlights the 3D interconnected nature of the pores or vugs (illustrated in green) and conductive patches (illustrated in red). It is noted that the length of the interval can be 3 feet (1 meter), and the borehole diameter can be 8.5 inches (22 cm).
- FIG. 12 shows the numerical pseudocore as Figure 10, and the complex 3D shape of the pores or vugs (illustrated green), and both the conductive patches and rock matrix are invisible. It is noted that the length of the interval can be 3 feet (1 meter), and the borehole diameter can be 8.5 inches (22 cm).
- FIG. 13 shows numerical pseudocores that can be resampled or regrided to radial grids wherein the radial grids can be layered, based upon layers observed in borehole images or other well logs.
- FIG. 13 also shows a radial grid that is generated to investigate the flow behavior of the numerical pseudocore model. Note that an inner part (approximately 4 inches in diameter) of the original numerical pseudocore has been drilled out to allow the arrangement of micro-producers. The micro-injectors can be located around the outer boundary of the pseudocore.
- FIG. 14a shows relative permeability curve that is used for conductive patch in the flow simulation; wherein there are three different relative permeability curves for each of the three rock types in the pseudocore model: such as, the vugs, conductive patches and tight rock matrix.
- FIG. 14a shows the relative permeability curve that was used for the conductive patch in the flow simulation.
- the pseudocore for both FIG. 14a and FIG. 14b can be approximately 1 foot (0.3 meters) high, the outer diameter is approximately 8.5 inches (22 cm) (or the width of the core is 8.5 in (22 cm) and the inner diameter is approximately 4 inches.
- FIG. 14b shows a flow simulation result through a numerical pseudocore, such that a line of micro-injectors of water surrounds outer the diameter, a line of micro-producers surrounds the inner diameter of the pseudocore: the colors represent oil saturation (So), where heterogeneous fingering and breakthrough are shown in this flow model, according to embodiments of the invention. It is noted that FIG. 14b shows the oil saturation profile of a waterflooded numerical pseudocore.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
BRPI0902889A BRPI0902889A2 (pt) | 2008-04-10 | 2009-04-10 | Método para criar um modelo de pseudonúcleo numérico, sistema para criação de um modelo de pseudonúcleo numérico, e sistema para criar um modelo de pseudonúcleo numérico. |
CN2009800001849A CN101802649B (zh) | 2008-04-10 | 2009-04-10 | 利用井眼图像、数字岩石样品以及多点统计算法生成数值假岩心的方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US4403108P | 2008-04-10 | 2008-04-10 | |
US61/044,031 | 2008-04-10 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2009126881A2 true WO2009126881A2 (fr) | 2009-10-15 |
WO2009126881A3 WO2009126881A3 (fr) | 2010-03-11 |
Family
ID=41162644
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2009/040198 WO2009126881A2 (fr) | 2008-04-10 | 2009-04-10 | Procédé pour produire des pseudo-carottes numériques utilisant des images de fond de puits, des échantillons de roche numériques, et des statistiques multi-emplacements |
Country Status (4)
Country | Link |
---|---|
CN (1) | CN101802649B (fr) |
BR (1) | BRPI0902889A2 (fr) |
RU (1) | RU2444031C2 (fr) |
WO (1) | WO2009126881A2 (fr) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2503010A (en) * | 2012-06-14 | 2013-12-18 | Reeves Wireline Tech Ltd | A method of processing geological log data |
CN104297791A (zh) * | 2014-09-25 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种基于地震优势频率的反演方法及系统 |
CN104570092A (zh) * | 2014-12-24 | 2015-04-29 | 中国石油天然气股份有限公司 | 一种识别碳酸盐岩大型缝洞型储层顶包络面的方法及装置 |
WO2014200996A3 (fr) * | 2013-06-10 | 2015-06-18 | Services Petroliers Schlumberger | Système pour compléter une lacune d'image de puits de forage |
US9563963B2 (en) | 2013-03-11 | 2017-02-07 | Reeves Wireline Technologies Limited | Methods of and apparatuses for identifying geological characteristics in boreholes |
EP3033484A4 (fr) * | 2013-08-13 | 2017-03-08 | Services Pétroliers Schlumberger | Analyse numérique de la sensibilité de carottes |
US20180225868A1 (en) * | 2015-07-28 | 2018-08-09 | Schlumberger Technology Corporation | Method and System for Generating A Virtual Core |
CN109684778A (zh) * | 2019-02-01 | 2019-04-26 | 长江地球物理探测(武汉)有限公司 | 一种数字岩心智能管理系统 |
US10295691B2 (en) | 2015-10-15 | 2019-05-21 | Reeves Wireline Technologies Limited | Method and apparatus for processing log data |
EP3394394A4 (fr) * | 2015-12-22 | 2019-08-14 | Landmark Graphics Corporation | Visualisation de tenseur de propriété de roche à base d'images de grille géocellulaire dans un environnement 3d dynamique |
US10753202B2 (en) | 2012-06-14 | 2020-08-25 | Reeves Wireline Technologies Limited | Geological log data processing methods and apparatuses |
CN111754627A (zh) * | 2020-05-30 | 2020-10-09 | 上海电力大学 | 一种结合卷积神经网络和页岩软数据的页岩重构方法 |
WO2020219035A1 (fr) * | 2019-04-24 | 2020-10-29 | Halliburton Energy Services, Inc. | Reconstructions de carottage de trou de forage à l'aide de balayages de trou de forage |
WO2020222757A1 (fr) * | 2019-04-29 | 2020-11-05 | Halliburton Energy Services, Inc. | Alignement de balayage de carotte automatisé |
CN112668254A (zh) * | 2020-12-28 | 2021-04-16 | 大连理工大学 | 一种考虑颗粒大小分布和孔隙度的致密数字岩心生成方法 |
US20210182597A1 (en) * | 2019-09-27 | 2021-06-17 | Carl Zeiss X-ray Microscopy, Inc. | Process parameter prediction using multivariant structural regression |
CN113052963A (zh) * | 2019-12-27 | 2021-06-29 | 中国石油天然气股份有限公司 | 三维岩性属性特征分析方法及装置 |
CN113553546A (zh) * | 2021-07-14 | 2021-10-26 | 中国海洋石油集团有限公司 | 一种电成像数据提取岩石连续胶结指数的方法、系统及计算机可读存储介质 |
WO2021258170A1 (fr) * | 2020-06-24 | 2021-12-30 | Petróleo Brasileiro S.A. - Petrobras | Script de programmation pour le traitement d'images et son application dans un procédé de détermination d'images faciologiques |
WO2022005993A1 (fr) * | 2020-07-01 | 2022-01-06 | Saudi Arabian Oil Company | Modélisation de réservoirs d'hydrocarbures à l'aide d'une classification de structure rocheuse dans les conditions du réservoir |
CN114166892A (zh) * | 2021-12-09 | 2022-03-11 | 中国矿业大学 | 一种基于网络并行电法的受载岩样损伤自电位成像方法 |
RU2773492C1 (ru) * | 2022-02-14 | 2022-06-06 | Федеральное государственное автономное образовательное учреждение высшего образования "Пермский государственный национальный исследовательский университет" | Способ оценки воздействия раствора на породу с целью повышения нефтеотдачи пласта |
US11353608B2 (en) | 2017-10-13 | 2022-06-07 | Schlumberger Technology Corporation | Method for determining a property of a material behind a borehole casing |
CN117251739A (zh) * | 2023-11-17 | 2023-12-19 | 核工业北京地质研究院 | 用于砂岩型铀资源定量预测的样本集的构建方法 |
WO2024018277A3 (fr) * | 2022-07-22 | 2024-04-04 | Cgg Services Sas | Analyse de données de trou de forage basée sur l'apprentissage automatique |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2263107A4 (fr) * | 2008-04-10 | 2016-12-28 | Services Petroliers Schlumberger | Procédé pour caractériser une formation géologique traversée par un forage |
CN102812203B (zh) * | 2011-04-01 | 2016-04-13 | Qri集团有限责任公司 | 用于通过性能量度的非对称分析来动态地评估石油储集层能力并提高产量和采收率的方法 |
US9097821B2 (en) * | 2012-01-10 | 2015-08-04 | Chevron U.S.A. Inc. | Integrated workflow or method for petrophysical rock typing in carbonates |
US20130262028A1 (en) * | 2012-03-30 | 2013-10-03 | Ingrain, Inc. | Efficient Method For Selecting Representative Elementary Volume In Digital Representations Of Porous Media |
RU2501046C1 (ru) * | 2012-05-10 | 2013-12-10 | Ильнар Ильясович Хасанов | Способ исследования керна горных пород |
US10241218B2 (en) * | 2012-05-30 | 2019-03-26 | Pgs Geophysical As | Methods and systems for computing notional source signatures from near-field measurements and modeled notional signatures |
RU2580174C1 (ru) * | 2012-06-09 | 2016-04-10 | Шлюмберже Текнолоджи Б.В. | Способ определения пористости образца породы |
US8666149B2 (en) * | 2012-08-01 | 2014-03-04 | Chevron U.S.A. Inc. | Method for editing a multi-point facies simulation |
WO2014104909A1 (fr) * | 2012-12-25 | 2014-07-03 | Общество С Ограниченной Ответственностью "Интровижн Ресерч Энд Девелопмент" | Procédé de production d'un modèle tridimensionnel caractéristique d'un échantillon de matériau poreux afin d'étudier les propriétés de perméabilité |
CN103114848B (zh) * | 2013-01-18 | 2015-09-30 | 西南石油大学 | 一种基于岩心测量的地层裂缝空间重构方法 |
WO2014168506A1 (fr) * | 2013-04-12 | 2014-10-16 | Schlumberger Canada Limited | Recuperation assistee de petrole a l'aide d'echantillon de forage numerique |
CN103306671B (zh) * | 2013-05-17 | 2015-04-08 | 中国石油天然气股份有限公司 | 一种四象限储层类型识别方法及系统 |
WO2015021182A1 (fr) * | 2013-08-06 | 2015-02-12 | Bp Corporation North America Inc. | Simulation numérique directe à base d'images des propriétés pétrophysiques sous des conditions simulées de contrainte et d'élongation |
US10422924B2 (en) | 2013-11-08 | 2019-09-24 | Schlumberger Technology Corporation | Stratigraphic and structural interpretation of deviated and horizontal wellbores |
FR3021796B1 (fr) * | 2014-06-02 | 2016-06-24 | Safran | Procede et dispositif de reconstruction numerique d'un volume elementaire representatif d'une microstructure de materiau composite |
CN104360415B (zh) * | 2014-10-31 | 2018-05-29 | 中国石油化工股份有限公司 | 一种致密砂岩储层裂缝识别的方法 |
RU2705629C1 (ru) * | 2015-11-09 | 2019-11-11 | Халлибертон Энерджи Сервисез, Инк. | Определение параметров ствола скважины с использованием ультразвуковых и микрорезистивных каверномеров |
CN106837315B (zh) * | 2015-12-03 | 2020-05-12 | 中国石油化工股份有限公司 | 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法 |
US10049172B2 (en) * | 2015-12-10 | 2018-08-14 | Saudi Arabian Oil Company | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution |
RU2651679C1 (ru) * | 2016-12-28 | 2018-04-23 | Акционерное общество "Всероссийский нефтегазовый научно-исследовательский институт имени академика А.П. Крылова" (АО "ВНИИнефть") | Способ создания синтетического образца керна с использованием трехмерной печати и компьютерной рентгеновской томографии |
BR112020009026A2 (pt) * | 2017-11-06 | 2020-10-27 | Khalifa University of Science and Technology | método e sistema para a determinação da permeabilidade de um meio poroso |
US20210165937A1 (en) * | 2017-12-14 | 2021-06-03 | Schlumberger Technology Corporation | System and Method for Simulating Reservoir Models |
CN108492331B (zh) * | 2018-03-20 | 2019-04-23 | 中国科学院地质与地球物理研究所 | 基于模拟压汞法获取岩石的孔径分布的方法和装置 |
CN109389128B (zh) | 2018-08-24 | 2021-08-27 | 中国石油天然气股份有限公司 | 电成像测井图像特征自动提取方法及装置 |
CN110134918B (zh) * | 2019-04-04 | 2021-05-28 | 中国石油大学(北京) | 基于滑动窗口法的岩心自动归位方法和系统 |
CN110987985A (zh) * | 2019-12-27 | 2020-04-10 | 西南石油大学 | 射孔损害室内评价数字岩心方法 |
CN112801205B (zh) * | 2021-02-20 | 2022-04-15 | 中国科学院地质与地球物理研究所 | 岩石纹层结构的获取方法及装置 |
CN112859197B (zh) * | 2021-03-31 | 2024-03-22 | 中国石油天然气集团有限公司 | 一种基于均质化地层电磁场理论的数字井筒电阻率模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4821164A (en) * | 1986-07-25 | 1989-04-11 | Stratamodel, Inc. | Process for three-dimensional mathematical modeling of underground geologic volumes |
US5838634A (en) * | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
US20070239359A1 (en) * | 2006-04-07 | 2007-10-11 | Chevron U.S.A. Inc. | Method for creating a stratigraphic model using pseudocores created from borehole images |
-
2009
- 2009-04-10 CN CN2009800001849A patent/CN101802649B/zh not_active Expired - Fee Related
- 2009-04-10 RU RU2009135608/28A patent/RU2444031C2/ru not_active IP Right Cessation
- 2009-04-10 WO PCT/US2009/040198 patent/WO2009126881A2/fr active Application Filing
- 2009-04-10 BR BRPI0902889A patent/BRPI0902889A2/pt not_active IP Right Cessation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4821164A (en) * | 1986-07-25 | 1989-04-11 | Stratamodel, Inc. | Process for three-dimensional mathematical modeling of underground geologic volumes |
US5838634A (en) * | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
US20070239359A1 (en) * | 2006-04-07 | 2007-10-11 | Chevron U.S.A. Inc. | Method for creating a stratigraphic model using pseudocores created from borehole images |
Cited By (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10753202B2 (en) | 2012-06-14 | 2020-08-25 | Reeves Wireline Technologies Limited | Geological log data processing methods and apparatuses |
US9715615B2 (en) | 2012-06-14 | 2017-07-25 | Reeves Wireline Technologies Limited | Geological log data processing methods and apparatuses |
GB2503010B (en) * | 2012-06-14 | 2018-04-18 | Reeves Wireline Tech Ltd | A method of processing geological log data |
US10255473B2 (en) | 2012-06-14 | 2019-04-09 | Reeves Wireline Technologies Limited | Geological log data processing methods and apparatuses |
GB2503010A (en) * | 2012-06-14 | 2013-12-18 | Reeves Wireline Tech Ltd | A method of processing geological log data |
GB2503061B (en) * | 2012-06-14 | 2018-04-25 | Reeves Wireline Tech Ltd | Geological log data processing methods and apparatuses |
US9563963B2 (en) | 2013-03-11 | 2017-02-07 | Reeves Wireline Technologies Limited | Methods of and apparatuses for identifying geological characteristics in boreholes |
US9576374B2 (en) | 2013-03-11 | 2017-02-21 | Reeves Wireline Technologies Limited | Methods of and apparatuses for identifying geological characteristics in boreholes |
US9704263B2 (en) | 2013-03-11 | 2017-07-11 | Reeves Wireline Technologies Limited | Methods of and apparatuses for identifying geological characteristics in boreholes |
WO2014200996A3 (fr) * | 2013-06-10 | 2015-06-18 | Services Petroliers Schlumberger | Système pour compléter une lacune d'image de puits de forage |
US10113411B2 (en) | 2013-06-10 | 2018-10-30 | Schlumberger Technology Corporation | Borehole image gap filling |
EP3033484A4 (fr) * | 2013-08-13 | 2017-03-08 | Services Pétroliers Schlumberger | Analyse numérique de la sensibilité de carottes |
AU2014307046B2 (en) * | 2013-08-13 | 2018-05-24 | Schlumberger Technology B.V. | Digital core sensitivity analysis |
US10465483B2 (en) | 2013-08-13 | 2019-11-05 | Schlumberger Technology Corporation | Digital core sensitivity analysis |
CN104297791A (zh) * | 2014-09-25 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种基于地震优势频率的反演方法及系统 |
CN104570092A (zh) * | 2014-12-24 | 2015-04-29 | 中国石油天然气股份有限公司 | 一种识别碳酸盐岩大型缝洞型储层顶包络面的方法及装置 |
CN104570092B (zh) * | 2014-12-24 | 2017-07-07 | 中国石油天然气股份有限公司 | 一种识别碳酸盐岩大型缝洞型储层顶包络面的方法及装置 |
EP3329091A4 (fr) * | 2015-07-28 | 2019-04-17 | Services Petroliers Schlumberger | Procédé et système de génération de carotte virtuelle |
US11982144B2 (en) | 2015-07-28 | 2024-05-14 | Schlumberger Technology Corporation | Method and system for generating a virtual core |
US20180225868A1 (en) * | 2015-07-28 | 2018-08-09 | Schlumberger Technology Corporation | Method and System for Generating A Virtual Core |
GB2543491B (en) * | 2015-10-15 | 2019-08-14 | Reeves Wireline Tech Ltd | Method of and apparatus for processing log data |
US10295691B2 (en) | 2015-10-15 | 2019-05-21 | Reeves Wireline Technologies Limited | Method and apparatus for processing log data |
US11060391B2 (en) | 2015-12-22 | 2021-07-13 | Landmark Graphics Corporation | Image based rock property tensor visualization of a geocellular grid in a dynamic 3D environment |
EP3394394A4 (fr) * | 2015-12-22 | 2019-08-14 | Landmark Graphics Corporation | Visualisation de tenseur de propriété de roche à base d'images de grille géocellulaire dans un environnement 3d dynamique |
US11353608B2 (en) | 2017-10-13 | 2022-06-07 | Schlumberger Technology Corporation | Method for determining a property of a material behind a borehole casing |
CN109684778A (zh) * | 2019-02-01 | 2019-04-26 | 长江地球物理探测(武汉)有限公司 | 一种数字岩心智能管理系统 |
WO2020219035A1 (fr) * | 2019-04-24 | 2020-10-29 | Halliburton Energy Services, Inc. | Reconstructions de carottage de trou de forage à l'aide de balayages de trou de forage |
US11836927B2 (en) | 2019-04-24 | 2023-12-05 | Halliburton Energy Services, Inc. | Borehole coring reconstructions using borehole scans |
WO2020222757A1 (fr) * | 2019-04-29 | 2020-11-05 | Halliburton Energy Services, Inc. | Alignement de balayage de carotte automatisé |
US20210182597A1 (en) * | 2019-09-27 | 2021-06-17 | Carl Zeiss X-ray Microscopy, Inc. | Process parameter prediction using multivariant structural regression |
CN113052963A (zh) * | 2019-12-27 | 2021-06-29 | 中国石油天然气股份有限公司 | 三维岩性属性特征分析方法及装置 |
CN111754627A (zh) * | 2020-05-30 | 2020-10-09 | 上海电力大学 | 一种结合卷积神经网络和页岩软数据的页岩重构方法 |
WO2021258170A1 (fr) * | 2020-06-24 | 2021-12-30 | Petróleo Brasileiro S.A. - Petrobras | Script de programmation pour le traitement d'images et son application dans un procédé de détermination d'images faciologiques |
WO2022005993A1 (fr) * | 2020-07-01 | 2022-01-06 | Saudi Arabian Oil Company | Modélisation de réservoirs d'hydrocarbures à l'aide d'une classification de structure rocheuse dans les conditions du réservoir |
US11592593B2 (en) | 2020-07-01 | 2023-02-28 | Saudi Arabian Oil Company | Modeling hydrocarbon reservoirs using rock fabric classification at reservoir conditions |
CN112668254A (zh) * | 2020-12-28 | 2021-04-16 | 大连理工大学 | 一种考虑颗粒大小分布和孔隙度的致密数字岩心生成方法 |
CN113553546B (zh) * | 2021-07-14 | 2023-08-18 | 中国海洋石油集团有限公司 | 一种电成像数据提取岩石连续胶结指数的方法、系统及计算机可读存储介质 |
CN113553546A (zh) * | 2021-07-14 | 2021-10-26 | 中国海洋石油集团有限公司 | 一种电成像数据提取岩石连续胶结指数的方法、系统及计算机可读存储介质 |
CN114166892A (zh) * | 2021-12-09 | 2022-03-11 | 中国矿业大学 | 一种基于网络并行电法的受载岩样损伤自电位成像方法 |
RU2773492C1 (ru) * | 2022-02-14 | 2022-06-06 | Федеральное государственное автономное образовательное учреждение высшего образования "Пермский государственный национальный исследовательский университет" | Способ оценки воздействия раствора на породу с целью повышения нефтеотдачи пласта |
WO2024018277A3 (fr) * | 2022-07-22 | 2024-04-04 | Cgg Services Sas | Analyse de données de trou de forage basée sur l'apprentissage automatique |
RU2792960C1 (ru) * | 2022-12-15 | 2023-03-28 | федеральное государственное автономное образовательное учреждение высшего образования "Пермский национальный исследовательский политехнический университет" | Способ оценки воздействия технологических жидкостей на горную породу |
CN117251739A (zh) * | 2023-11-17 | 2023-12-19 | 核工业北京地质研究院 | 用于砂岩型铀资源定量预测的样本集的构建方法 |
CN117251739B (zh) * | 2023-11-17 | 2024-01-30 | 核工业北京地质研究院 | 用于砂岩型铀资源定量预测的样本集的构建方法 |
Also Published As
Publication number | Publication date |
---|---|
RU2444031C2 (ru) | 2012-02-27 |
BRPI0902889A2 (pt) | 2017-08-29 |
CN101802649B (zh) | 2013-01-23 |
CN101802649A (zh) | 2010-08-11 |
WO2009126881A3 (fr) | 2010-03-11 |
RU2009135608A (ru) | 2011-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8725477B2 (en) | Method to generate numerical pseudocores using borehole images, digital rock samples, and multi-point statistics | |
RU2444031C2 (ru) | Способ генерирования численных псевдокернов с использованием изображений скважины, цифровых образов породы и многоточечной статистики | |
US9581723B2 (en) | Method for characterizing a geological formation traversed by a borehole | |
US9134457B2 (en) | Multiscale digital rock modeling for reservoir simulation | |
Skalinski et al. | Carbonate petrophysical rock typing: integrating geological attributes and petrophysical properties while linking with dynamic behaviour | |
Dell'Arciprete et al. | Comparison of three geostatistical methods for hydrofacies simulation: a test on alluvial sediments | |
US9939548B2 (en) | Systems, methods, and computer medium to produce efficient, consistent, and high-confidence image-based electrofacies analysis in stratigraphic interpretations across multiple wells | |
Liu et al. | Multiple-point simulation integrating wells, three-dimensional seismic data, and geology | |
Slatt et al. | Pore-to-regional-scale integrated characterization workflow for unconventional gas shales | |
Hurley et al. | Method to generate full-bore images using borehole images and multipoint statistics | |
Williams-Stroud et al. | Microseismicity-constrained discrete fracture network models for stimulated reservoir simulation | |
Zhang | MPS-driven digital rock modeling and upscaling | |
Shahinpour | Borehole image log analysis for sedimentary environment and clay volume interpretation | |
Catinat et al. | Estimating permeability in a limestone geothermal reservoir from NMR laboratory experiments | |
Holden et al. | Integration of production logs helps to understand heterogeneity of Mishrif reservoir in Rumaila | |
Hyndman et al. | Hydrogeophysical case studies at the local scale: The saturated zone | |
Poduri et al. | On the performance of pilot‐point based hydraulic tomography with a geophysical a priori model | |
Worthington | Reservoir characterization at the mesoscopic scale | |
Gao et al. | Seismic structure and texture analyses for fractured reservoir characterization: An integrated workflow | |
Hurley et al. | Multiscale workflow for reservoir simulation | |
Verma et al. | Imaging and interpretation: Seismic, rock physics and image log analysis workflows for deepwater systems | |
Martinsen et al. | Outcrops revitalized: Tools, techniques and applications | |
Jia et al. | Advances and challenges of reservoir characterization: A review of the current state-of-the-art | |
Williams-Rojas et al. | Geologic controls on reservoir performance in Muspac and Catedral gas fields, southeastern Mexico | |
Zhang et al. | High-Resolution Mineralogy Modeling-A Case Study in the Vaca Muerta Formation, Neuquén Basin, Argentina |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200980000184.9 Country of ref document: CN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009135608 Country of ref document: RU |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09730368 Country of ref document: EP Kind code of ref document: A2 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 09730368 Country of ref document: EP Kind code of ref document: A2 |
|
ENP | Entry into the national phase |
Ref document number: PI0902889 Country of ref document: BR Kind code of ref document: A2 Effective date: 20090921 |