[go: up one dir, main page]

US20100313633A1 - Estimating effective permeabilities - Google Patents

Estimating effective permeabilities Download PDF

Info

Publication number
US20100313633A1
US20100313633A1 US12/637,958 US63795809A US2010313633A1 US 20100313633 A1 US20100313633 A1 US 20100313633A1 US 63795809 A US63795809 A US 63795809A US 2010313633 A1 US2010313633 A1 US 2010313633A1
Authority
US
United States
Prior art keywords
rocks
measurements
predicted
permeabilities
effective
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US12/637,958
Other languages
English (en)
Inventor
Vivek Anand
Robert Freedman
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Schlumberger Technology Corp
Original Assignee
Schlumberger Technology Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Schlumberger Technology Corp filed Critical Schlumberger Technology Corp
Priority to US12/637,958 priority Critical patent/US20100313633A1/en
Assigned to SCHLUMBERGER TECHNOLOGY CORPORATION reassignment SCHLUMBERGER TECHNOLOGY CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: ANAND, VIVEK, FREEDMAN, ROBERT
Priority to BRPI1001536-1A priority patent/BRPI1001536A2/pt
Priority to MX2010005336A priority patent/MX2010005336A/es
Publication of US20100313633A1 publication Critical patent/US20100313633A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/084—Detection of potentially hazardous samples, e.g. toxic samples, explosives, drugs, firearms, weapons
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/448—Relaxometry, i.e. quantification of relaxation times or spin density

Definitions

  • Implementations of various technologies described herein generally relate to techniques for determining effective permeabilities of earth formations and, more particularly, to techniques for determining such effective permeabilities using measurements in the earth formations.
  • the interpretation of well logging and geophysical measurements generally involves formulating and solving a mathematical inverse problem. That is, typically, one would like to predict the physical properties of some underlying physical system, such as effective permeabilities, from a suite of measurements.
  • the suite of measurements could be from a borehole logging tool or a suite of tools for which the underlying physical system is the porous fluid-filled rock formations surrounding the borehole.
  • the physical properties predicted from the measurements might include porosities, fluid types and saturations, and bed thicknesses.
  • the measurements could be surface measurements of reflected seismic wave energy as a function of wavelength made at different receiver locations.
  • the underlying physical system is the subsurface consisting of layers of porous sediments. The physical properties of most interest are those of the hydrocarbon-bearing layers.
  • Effective permeability is one of the physical properties of the underlying physical system that may be predicted by solving the inverse problems. Effective permeability is the ability to preferentially flow or transmit a particular fluid through a rock in the presence of other immiscible fluids in the reservoir. The estimation of effective permeability assists in reservoir development and management. For example, permeability is used for determining production rates and optimal drainage points, optimizing completion and perforation design, and devising enhanced oil recovery (EOR) patterns and injection conditions.
  • EOR enhanced oil recovery
  • the first method uses pressure transient tests, such as formation testers (e.g., Schlumberger Modular Dynamics Tester (MDT)), to measure the transient build up in the pore pressure following an extraction of a fixed volume of formation fluid.
  • MDT Schlumberger Modular Dynamics Tester
  • the effective permeability (k e ) of the formation can be related to the pressure build up.
  • transmissivity (k e h/ ⁇ ) during radial flow
  • h reservoir thickness
  • viscosity of the fluid ⁇
  • pressure measurements are influenced by the presence of skin in the near-probe region.
  • the estimates of permeability from pressure measurements may be incorrect.
  • the pressure transient tests usually yield the effective permeability of the mud filtrate in the invaded zone rather than the effective permeability of formation fluids measured.
  • the estimation of permeability from transient tests requires matching the transient to type curves and formation models. Because of these factors, the permeability estimates from pressure transient tests remain qualitative.
  • the second method for in-situ estimation of permeability uses continuous log data. These data provide a continuous survey of formation properties such as porosity, irreducible water saturation and Nuclear Magnetic Resonance (NMR) parameters. Empirical and semi-empirical correlations have been developed that relate the absolute permeability of the formation to NMR parameters. The following two correlations, called the Schlumberger Doll Research (SDR) model and the Timur-Coates model respectively, are commonly employed for estimation of permeability from NMR log data:
  • SDR Schlumberger Doll Research
  • Timur-Coates model are commonly employed for estimation of permeability from NMR log data:
  • T 2LM is the logarithmic mean of the water T 2 distribution
  • FFI and BVI are the free fluid index and bound volume irreducible.
  • a lithology specific T 2,cutoff is employed to partition the T 2 distribution into bound and free fluid components.
  • a major limitation of determining permeabilities using the equations above is that the parameters a SDR , a coates and T 2cutoff are not universal and need to be calibrated for each reservoir area. Additionally, the correlations provide estimates of the absolute permeability (permeability at 100% water saturation) of the formation and not the effective permeability, which is the more useful parameter.
  • the third method for in-situ estimation of permeability includes using production tests and production history.
  • An estimate of in-situ permeability can be obtained from flow rate and pressure data during steady-state production, preferably from specific tests at different flow rates.
  • Another estimation method involves adjusting the permeability to match a history of production data. Both methods suffer from non-uniqueness of the solution of highly non-linear inverse problems. Furthermore, only an average value of permeability can be obtained.
  • a method for determining effective permeabilities of earth formations may include receiving a database having one or more measurements made on a collection of fluid filled rocks. The method may then include dividing the measurements into input measurements and output measurements. The input measurements may describe one or more measured properties of the fluid filled rocks and the output measurements may describe effective permeabilities of the fluid filled rocks. The method may then include constructing a mapping function using the input measurements and the output measurements. After the mapping function is constructed, the method may include receiving one or more input measurements made on one or more rocks that are not part of the collection of fluid filled rocks. The method may then include predicting effective permeabilities of the rocks using the mapping function and the input measurements made on the rocks.
  • FIG. 1 illustrates a schematic diagram of a logging apparatus in accordance with implementations of various techniques described herein.
  • FIG. 2 illustrates a graph indicating a predicted absolute permeability for carbonate cores estimated by a SDR and a Timur-Coates model versus measured effective permeability to oil in accordance with implementations of various techniques described herein.
  • FIG. 3 illustrates a graph indicating a predicted absolute permeability for sandstone cores estimated by a SDR and a Timur-Coates model versus measured effective permeability to oil in accordance with implementations of various techniques described herein.
  • FIG. 4 illustrates flow diagram of a method for estimating effective permeabilities in accordance with implementations of various techniques described herein.
  • FIG. 5 illustrates a graph indicating a predicted effective permeability to oil for carbonate cores estimated by a radial basis function interpolation technique in accordance with implementations of various techniques described herein.
  • FIG. 6 illustrates a graph indicating a predicted effective permeability to oil for sandstone cores estimated by a radial basis function interpolation technique in accordance with implementations of various techniques described herein.
  • FIG. 7 illustrates a computer network into which implementations of various technologies described herein may be implemented.
  • a computer application may receive a database that includes measurements made on a collection of fluid filled rocks. The measurements may have been made using a well logging device or in a laboratory. In either case, the computer application may divide the measurements in the database into input measurements and output measurements. The input measurements may include one or more measured properties of the fluid filled rocks, and the output measurements may include the corresponding effective permeabilities of the fluid filled rocks. The computer application may then generate the mapping function by correlating the input measurements to their corresponding output measurements (i.e., effective permeabilities).
  • the computer application may receive one or more input measurements pertaining to one or more rocks that were not part of the collection of fluid filled rocks that was used to create the mapping function. The computer application may then predict the output measurements or the effective permeabilities of the rocks that were not part of the collection of fluid filled rocks using the mapping function and the input measurements of the rocks.
  • FIGS. 1-6 illustrate one or more implementations of various techniques described herein in more detail.
  • FIG. 1 shows a borehole 32 that has been drilled in formations 31 with drilling equipment, and typically, using drilling fluid or mud that results in a mudcake represented at 35 .
  • a logging device 100 is shown, and can be used in connection with various implementations described herein.
  • the logging device 100 may be suspended in the borehole 32 on an armored multiconductor cable 33 .
  • Known depth gauge apparatus (not shown) is provided to measure cable displacement over a sheave wheel (not shown), and thus the depth of the logging device 100 in the borehole 32 .
  • Circuitry 51 represents control and communication circuitry for the investigating apparatus. Although circuitry 51 is shown at the surface, portions thereof may typically be downhole. Also shown at the surface are processor 50 and recorder 90 .
  • the logging device 100 is shown to be a wireline logging tool, it should be noted that other tools, such as a logging while drilling tool, may be used in connection with various implementations described herein.
  • the logging device 100 may represent any type of logging device that takes measurements from which formation characteristics can be determined, for example, by solving complex inverse problems.
  • the logging device 100 may be an electrical type of logging device (including devices such as resistivity, induction, and electromagnetic propagation devices), a nuclear logging device, a sonic logging device, or a fluid sampling logging device, or combinations thereof.
  • Various devices may be combined in a tool string and/or used during separate logging runs. Also, measurements may be taken during drilling and/or tripping and/or sliding.
  • Examples of the types of formation characteristics that can be determined using these types of devices include: determination, from deep three-dimensional electromagnetic measurements, of distance and direction to faults or deposits, such as salt domes or hydrocarbons; determination, from acoustic shear and/or compressional wave speeds and/or wave attenuations, of formation porosity, permeability, and/or lithology; determination of formation anisotropy from electromagnetic and/or acoustic measurements; determination, from attenuation and frequency of a rod or plate vibrating in a fluid, of formation fluid viscosity and/or density; determination, from resistivity and/or nuclear magnetic resonance (NMR) measurements, of formation water saturation and/or permeability; determination, from count rates of gamma rays and/or neutrons at spaced detectors, of formation porosity and/or density; and determination, from electromagnetic, acoustic and/or nuclear measurements, of formation bed thickness.
  • determination from deep three-dimensional electromagnetic measurements, of distance and direction to faults or deposits, such as salt
  • FIG. 2 illustrates a graph 200 indicating a predicted absolute permeability for carbonate cores estimated by a SDR and a Timur-Coates model versus measured effective permeability to oil in accordance with implementations of various techniques described herein.
  • Graph 200 compares the absolute permeabilities of 37 carbonate cores estimated by the SDR and the Timur model with the effective oil permeabilities measured in a laboratory.
  • the solid black line is the best-fit line and the dashed lines are located at a deviation factor of 3.
  • the SDR and Timur-Coates model estimates provide a poor correlation between the estimated absolute permeabilities and the measured effective permeabilities for carbonate cores.
  • FIG. 3 illustrates a graph 300 indicating a predicted absolute permeability for sandstone cores estimated by a SDR and a Timur-Coates model versus measured effective permeability to oil in accordance with implementations of various techniques described herein.
  • Graph 300 compares the absolute permeabilities of 80 sandstone cores estimated by the SDR and the Timur model with the effective oil permeabilities measured in a laboratory.
  • the solid black line is the best-fit line and the dashed lines are located at a deviation factor of 3.
  • the SDR and Timur-Coates model estimates provide a poor correlation between the estimated absolute permeabilities and the measured effective permeabilities for sandstone cores.
  • FIGS. 2 and 3 illustrates a traditional approach (e.g., using simple empirically derived equations like the SDR and Timur-Coates models) for solving mathematical inverse problems that may be used to interpret well logging measurements obtained from logging device 100 or to interpret geophysical measurements obtained from a laboratory.
  • the traditional approach includes fitting a theoretical or empirically derived forward model (e.g., SDR model, Timur-Coates model) to measurement data (see e.g., the book by A. Tarantola, “Inverse Problem Theory: Methods For Data Fitting And Model Parameter Estimation”, published by Elsevier, Amsterdam, The Netherlands, 1987).
  • the forward model is a function of a set of model parameters that are either identical to or related to the physical properties of the underlying physical system. Selecting the values of the model parameters that minimize the difference between the actual measurements and those predicted by the forward model is assumed to solve the inverse problem.
  • This basic assumption is itself fraught with difficulties and can lead to erroneous solutions because most well logging and geophysical inverse problems are ill posed, i.e., the solutions are not unique.
  • This traditional approach has other inherent limitations and caveats that render it unsuitable or too computationally expensive for providing accurate solutions to many problems of interest.
  • FIG. 4 illustrates flow diagram of a method 400 for estimating effective permeabilities in accordance with implementations of various techniques described herein. It should be understood that while method 400 indicates a particular order of execution of the operations, in some implementations, certain portions of the operations might be executed in a different order. In one implementation, method 400 may be performed by a system computer which will be described in more detail with reference to FIG. 7 .
  • the system computer may receive a database containing measurements that have been made on a collection of fluid filled rocks.
  • the measurements may be obtained from laboratory measurements with core plugs.
  • the laboratory measurements may be obtained by performing one or more tests on the collection of fluid filled rocks to determine certain characteristics of the fluid filled rocks, such as fluid porosity, saturation, Nuclear Magnetic Resonance (NMR) T2 response, NMR T1 response, viscosity, effective permeabilities and the like.
  • the data may be obtained from well logging devices, such as the logging device 100 illustrated in FIG. 1 .
  • the fluid filled rocks may be rocks or earth formations of any type that may be found at or near wells, such as carbonates, sandstones and the like.
  • the input measurements may include well log measurements that may be made routinely by logging service companies.
  • Porosity may be the most basic well logging measurement. Porosity may be determined from neutron and density logs, NMR derived porosities or combinations thereof. Porosity may also be derived from acoustic or dielectric well logging measurements.
  • Water saturation may also be derived from resistivity and dielectric logs using porosity and other log inputs. Water saturation may also be derived from NMR tool diffusion measurements. The derivation of porosity and water saturation from well log data is well-known to anyone skilled in the art of well-logging formation evaluation. Nuclear Magnetic Resonance (NMR) T2 response may be derived from NMR logging tool measurements.
  • NMR Nuclear Magnetic Resonance
  • the system computer may divide the measurements in the database into input and output measurements.
  • dividing the measurements into input and output measurements may include designating a portion of the measurements made on the collection of fluid filled rocks as input measurements and the remaining portion as output measurements.
  • Both the input and output measurements may include various formation properties of the collection of fluid filled rocks, but the output measurements may include information that is being sought.
  • the output measurements may include the effective permeabilities of the collection of fluid filled rocks because one may be seeking to predict the effective permeabilities of one or more rocks that are not part of the collection of fluid filled rocks.
  • the method described herein is directed at estimating effective permeabilities, it should be noted that the method described herein may also be used to estimate various other properties of rocks.
  • the system computer may generate a mapping function based on the input measurements and the output measurements identified at step 420 .
  • the mapping function may approximate the underlying physical relationship between the input measurements and the output measurements.
  • the mapping function may approximate the relationship between the characteristics of the fluid filled rocks, such as fluid porosity, saturation, NMR T2 response, NMR T1 response (i.e., input measurements) and the effective permeabilities of the fluid filled rocks (i.e., output measurements).
  • the mapping function may be visualized as a multivariate interpolation between the input measurements in the database and the effective permeabilities.
  • the mapping function may be generated or constructed using a linear combination of one or more non-linear functions or using a weighted sum of one or more non-linear functions.
  • the mapping function may be generated using radial basis functions.
  • ⁇ right arrow over ( ⁇ ) ⁇ ( ⁇ right arrow over (x) ⁇ ), ⁇ right arrow over (x) ⁇ R′′ and ⁇ right arrow over ( ⁇ ) ⁇ R′′′ be a real valued vector function of n variables, and let values of ⁇ right arrow over ( ⁇ ) ⁇ ( ⁇ right arrow over (x) ⁇ i ) ⁇ right arrow over (y) ⁇ i be given at N distinct points, ⁇ right arrow over (x) ⁇ i .
  • the interpolation problem is to construct function ⁇ right arrow over (F) ⁇ ( ⁇ right arrow over (x) ⁇ ), that approximates ⁇ right arrow over (f) ⁇ ( ⁇ right arrow over (x) ⁇ ) and satisfies the interpolation equations,
  • RBF interpolation chooses an approximating or mapping function of the form
  • the non-linear functions ⁇ ( ⁇ right arrow over (x) ⁇ right arrow over (x i ) ⁇ ) are called “radial” because the argument of the function depends only on the distance between ⁇ right arrow over (x i ) ⁇ and an arbitrary input vector ⁇ right arrow over (x) ⁇ .
  • the argument is given by the Euclidean norm in the n-dimensional hyper space, i.e.,
  • the weights or coefficients, ⁇ right arrow over (c) ⁇ i in Equation (2) are determined by requiring that the interpolation equations (1) be satisfied exactly.
  • the system computer may calibrate the coefficients of the mapping function such that the interpolation of the input measurements to the output measurements is exact. As such, the coefficients are a linear combination of the given function,
  • Generating the mapping function may include solving inverse problems that involve predicting the physical properties of an underlying system (i.e., output measurements), given the set of input measurements.
  • the mapping function may include solving inverse problems that involve predicting the physical properties of an underlying system (i.e., output measurements), given the set of input measurements.
  • the output measurements ⁇ right arrow over (y) ⁇ i represent samples of the function that the system computer may want to approximate and ⁇ right arrow over (x) ⁇ i , are the distinct points at which the function is given.
  • the input measurements, ⁇ right arrow over (x) ⁇ i represent the measurements from which the output measurements, ⁇ right arrow over (y) ⁇ i , of the underlying system are to be predicted.
  • the output measurements, ⁇ right arrow over (y) ⁇ i may include the physical properties of the underlying system, such as effective permeabilities.
  • the mapping function may be configured such that given input measurements ⁇ right arrow over (x) ⁇ that are not in the database received at step 410 , the system computer may predict the output measurements, ⁇ right arrow over (y) ⁇ ( ⁇ right arrow over (x) ⁇ ), (i.e., permeabilities) of the physical system consistent with the input measurements. As such, the mapping function solves the inverse problem by predicting the physical properties of the system from the input measurements.
  • the radial basis functions used in the implementations described herein may be normalized Gaussian radial basis functions defined by the equation,
  • radial basis functions such as exponential, multiquadrics, or inverse multiquadrics, may also be used. These functions may be normalized in the sense that the summation over the input measurements, ⁇ right arrow over (x) ⁇ i , is equal to unity for all ⁇ right arrow over (x) ⁇ , i.e.,
  • the width, S i of the Gaussian radial basis function centered at ⁇ right arrow over (x) ⁇ i is representative of the range or spread of the function in the input space.
  • the optimal widths, for accurate approximations, should be of the order of the nearest neighbor distances in the input space. The idea is to pave the input space with basis functions that have some overlap with nearest neighbors but negligible overlap for more distant neighbors.
  • Equation 8 An intuitive understanding of how the mapping function in Equation 8 predicts an output vector for an input vector not in the data base can be gained by considering the Nadaraya-Watson Regression Estimator (NWRE).
  • NWRE Nadaraya-Watson Regression Estimator
  • the NWRE is based on a simple approximation for the weight vectors (S. Haykin, Neural Networks: A Comprehensive Foundation , Prentice Hall, Hamilton Ontario, Canada, 1999).
  • the interpolation equations for the mapping function in Equation 8 can be written in the form,
  • Equation (9) can be neglected if one neglects the overlap of the data base radial basis functions.
  • NWRE approximation assumes that the interpolation matrix in Equation 4 is diagonal and leads to a simple approximation for the coefficient vectors,
  • ⁇ right arrow over (F) ⁇ ( ⁇ right arrow over (x) ⁇ ) approaches the sample mean of the data base output vectors.
  • ⁇ right arrow over (F) ⁇ ( ⁇ right arrow over (x) ⁇ ) approaches the output vector ⁇ right arrow over (y) ⁇ j corresponding to the database input vector ⁇ right arrow over (x) ⁇ j that is closest ⁇ right arrow over (x) ⁇ .
  • ⁇ right arrow over (F) ⁇ ( ⁇ right arrow over (x) ⁇ ) is a weighted average of the database output vectors with weighting factors determined by the closeness of ⁇ right arrow over (x) ⁇ the database input vectors. It can be observed that the NWRE approximation in Equation 11 does not satisfy the interpolation conditions in Equation 1.
  • the NWRE approximation can be improved upon by determining optimal coefficient vectors such that the interpolation equations are satisfied.
  • the problem is linear if the widths of the Gaussian radial basis functions are fixed.
  • the interpolation conditions lead to a set of linear equations for the coefficient vectors whose solution can be written in matrix form, i.e.,
  • Equation 12 the matrix ⁇ whose inverse appears in Equation 12 is the, N ⁇ N , positive definite matrix of Gaussian radial basis functions, i.e.,
  • the N ⁇ m matrix, Yin Equation 12 contains the database output vectors, e.g.,
  • Equation 12-16 improves on the NWRE approximation by determining optimal coefficient vectors with the caveat of having fixed widths for the Gaussian radial basis functions. It can be proved mathematically that the matrix ⁇ is non-singular for certain functional forms of RBFs, including Gaussian, multiquadric, and inverse quadrics. This property ensures that the mapping function of Equation (2) is unique.
  • a mapping, i.e., interpolating, function that is consistent with the measurements can be uniquely defined from Equation (13).
  • the desired output may then be obtained by evaluating the mapping function at the corresponding input i.e.,
  • the system computer may receive one or more input measurements pertaining to one or more rocks that are not part of the collection of fluid filled rocks in the database received at step 410 .
  • the rocks that are not part of the collection of fluid filled rocks include rocks that are related to hydrocarbons.
  • the input measurements pertaining to the rocks that are not part of the collection of fluid filled rocks may include characteristics, such as fluid porosity, saturation, NMR T2 response, NMR T1 response and the like.
  • the system computer may then predict the effective permeabilities (i.e., output measurements) from the input measurements received at step 440 using the mapping function created at step 430 .
  • the mapping function generated at step 430 may be used to predict the effective permeabilities of rocks using the input measurements on those rocks.
  • effective permeabilities are usually denoted by the symbol, k o (S w ).
  • many oil reservoirs encountered in practice are at irreducible water saturation, k o (S wi ), and the mapping functions generated at step 430 may be directly applied to borehole well logging measurements in such reservoirs to predict the corresponding effective permeabilities of the rocks within the borehole.
  • the mapping functions generated at step 430 may be directly applied to borehole well logging measurements in reservoirs that are at other saturations to predict the corresponding effective permeabilities of the rocks within the borehole.
  • the effective permeability to water at irreducible water saturation is defined as k w (S wi ) ⁇ 0.
  • oil reservoirs at irreducible water saturation should flow oil and no water (i.e., the water cut should be zero).
  • a continuous depth log of effective permeabilities to oil, k o (S wi ) would be a useful new reservoir quality characterization parameter.
  • the continuous depth log of effective permeabilities may be derived by the predicted effective permeabilities.
  • the continuous depth log of effective permeabilities to oil may be useful in making well to well comparisons of reservoir quality in a development field.
  • the continuous depth log of effective permeabilities to oil may also be useful for selecting zones to complete a single well and in choosing perforation depths for optimal flow rates in a single well.
  • the continuous depth logs of the earth formations in a reservoir may be predicted using the predicted effective permeabilities, predicted flow rates, predicted mobilities, a fractional water flow of rocks, or combinations thereof.
  • One or more implementations for predicting flow rates, mobilities and fractional water flow of rocks are described below.
  • the predicted effective permeabilities may be used to predict the flow rates of one or more fluids in the rocks that are not part of the collection of fluid filled rocks.
  • the predicted effective permeabilities may be used to construct a “producibility index” from the oil mobility (M) which is defined as the effective permeability to oil divided by oil viscosity ( ⁇ ), e.g.,
  • the flow rate is proportional to mobility and a depth log of mobility in a borehole may be used as a parameter for choosing completion zone depths and location of perforations in a well.
  • the effective permeabilities may also be used to predict the mobilities of fluids in the rocks that are not part of the collection of fluid filled rocks.
  • the mobilities of fluids may then be used to select casing perforation depths in a borehole in order to optimize production rates.
  • Effective permeabilities may also be used as inputs in reservoir simulation models.
  • a more quantitative prediction of flow rates i.e., production rate
  • a more quantitative prediction of flow rates may be obtained by solving Navier-Stokes equations or Darcy's equations for multi-phase flow. Both of these equations are employed in reservoir engineering simulations.
  • One input to the simulations may include the effective permeabilities of the reservoir fluids as a function of the wetting phase fluid saturation.
  • These equations and their solutions are well-known to those skilled in reservoir engineering and flow in porous media.
  • Darcy's equation in differential form for the flow rate of oil in the presence of water can be written in the form, (e.g., see R. E. Collins, Flow of Fluids Through Porous Media, pp. 60-62),
  • Equation 19 ⁇ right arrow over (v) ⁇ o is the mean flow rate per unit area, k o is the effective permeability to oil and is a function of the wetting phase saturation, p o is the oil density, ⁇ o is the oil viscosity in the reservoir, and ⁇ o is the flow potential for the oil phase. It is understood that a similar equation for water flow can be written.
  • the mobilities of oil and water may then be used to predict the fraction of the total flow that will be water.
  • the fraction of the total flow that will be water is the water-cut and is given approximately by (see equation 6-28 in Collins),
  • Equation (20) is the water cut when water is displacing oil and may be used for a secondary recovery water flood or for primary production by a water drive.
  • the simple form shown in Equation (20) may neglect capillary pressure and gravity effects.
  • method 400 may also be used to predict relative permeabilities.
  • Relative permeabilities may be predicted by calculating the ratio of the effective permeability to an absolute permeability.
  • the absolute permeability is a measurement of the permeability of a rock filled with a single fluid.
  • the database containing measurements received at step 410 may include relative permeability data.
  • the prediction of accurate effective fluid permeabilities may be an ingredient in the prediction of reservoir performance.
  • method 400 detailed above will lead to improved predictions of key reservoir performance factors including production rates, water cut, recoverable reserves, residual oil saturation, ultimate recovery and the like.
  • FIG. 5 illustrates a graph indicating a predicted effective permeability to oil for carbonate cores estimated by a radial basis function interpolation technique in accordance with implementations of various techniques described herein.
  • the radial basis function interpolation technique may describe the mapping function generated at step 430 .
  • FIG. 5 uses a world-wide rock database consisting of petrophysical input measurements on carbonate cores to predict effective permeabilities of the carbonate cores.
  • the carbonate cores may have been obtained from formations around the world and may have included carbonate cores of different geologic ages ranging from the Miocene period to the Ordovician period.
  • Carbonate cores with a wide range of petrophysical and geological properties such as lithology, texture and porosity types may also have been included in the database.
  • the database used in FIG. 5 consisted of carbonate rocks of two important lithologies namely, limestone and dolostone.
  • the carbonate cores comprised of grainstone, packstone, wackestone, mudstone, boundstone, and crystalline texture.
  • the porosity types of the cores included interparticle, intraparticle, intercrytalline, intracrystalline, moldic, vugular (touching and isolated) and fenestral.
  • the pore filling minerals included calcite, dolomite, chert, anhydrite, clay and solid hydrocarbon.
  • the input measurement database consisted of core porosity, irreducible water saturation, effective permeability to oil at irreducible water saturation, and NMR response at irreducible water saturation.
  • the porosity of the cores measured using helium expansion method varied from 5% to 35%.
  • the effective permeability of oil at irreducible water saturation and 5000 psig confining pressure varied from 0.1 md to 1000 md.
  • NMR response of oil saturated plugs at irreducible water saturation was measured at 2 MHz and 0.2 ms echo spacing.
  • the effective permeabilities of the carbonate cores may be expressed as a linear combination of RBFs as shown below.
  • ⁇ right arrow over (A) ⁇ r is the input vector which includes porosity, irreducible water saturation and normalized amplitudes of the T 2 distribution.
  • ⁇ right arrow over (A) ⁇ T is defined as:
  • the amplitudes of the T 2 distribution for each sample are normalized with the largest respective values to eliminate the dependence on hardware and software settings.
  • the widths of the Gaussian RBFs are proportional to the Euclidean nearest neighbor distance in the input space.
  • Equation (21) the effective permeabilities of the carbonate cores were estimated from porosity, irreducible water saturation and amplitudes of T 2 distribution.
  • the effective permeabilities of the samples were calculated from Equation (21) using the using the leave-one out method.
  • the widths were heuristically determined to be half the nearest neighbor distances in the input space.
  • the comparison of the estimated permeabilities with those measured in the laboratory is plotted in the graph of FIG. 5 . As seen in FIG.
  • the effective oil permeabilities may be accurately predicted for suites of carbonate rocks using measurements of total rock porosity ( ⁇ ), irreducible water saturation (S wi ), and a T2 distribution.
  • ⁇ total rock porosity
  • S wi irreducible water saturation
  • T2 total rock porosity
  • the permeability is estimated within a factor of 3 which is a significant improvement compared to the estimates of physical models as shown in FIG. 2 .
  • FIG. 6 illustrates a graph indicating a predicted effective permeability to oil for sandstone cores estimated by a radial basis function interpolation technique in accordance with implementations of various techniques described herein.
  • a world-wide rock database consisting of petrophysical measurements on sandstone cores was used to determine the effective permeabilities to oil.
  • the plugs were obtained from formations from around the world and were of different geologic ages ranging from the Pliocene period to the Devonian period.
  • the database incorporated core plugs from formations with a wide range of petrophysical and geological properties such as grain size, sorting, degree of consolidation and cement types.
  • the grain size of the core plugs ranged from silt size ( ⁇ 0.06 mm) to pebble size (>2 mm).
  • the grain sorting varied from very well to very poor sorting.
  • the database also included unconsolidated sands as well as consolidated sands with varying degree of consolidation.
  • the input measurements with the core plugs included porosity, irreducible water saturation, effective permeability to oil at irreducible water saturation, absolute permeability (100% water saturated), and T 2 distribution at irreducible water saturation.
  • the porosity of the cores measured using helium expansion method varied from 5% to 35%.
  • the effective permeability of oil measured at irreducible water saturation at 5000 psig confining pressure varied from 0.1 md to 1000 md.
  • the T 2 distributions of the oil saturated core plugs at irreducible water saturation were measured at proton Larmor frequency of 2 MHz and 0.2 ms echo spacing.
  • the NMR response and effective permeability were measured with different cores that were derived from the same larger-sized ( ⁇ 1 foot) plug.
  • the interpolation between T 2 distribution and effective permeability incorporates an error due to the heterogeneity of the formation over the length scale of the plug. This is not necessarily a disadvantage because log data, in general, may have similar or worse vertical resolution.
  • the effective permeabilities of the cores can be expressed as a linear combination of RBFs as shown in Equations 21-22.
  • FIG. 6 shows the comparison of the effective permeabilities estimated using Equation (21) with those measured in the lab.
  • the effective oil permeabilities may be accurately predicted for suites of sandstones using measurements of total rock porosity ( ⁇ ), irreducible water saturation (S wi ), and a T2 distribution.
  • the effective permeability is estimated within a factor of 3, which is a significant improvement compared to the estimates of the physical models as shown in FIG. 3 .
  • FIG. 7 illustrates a computing system 700 , into which implementations of various techniques described herein may be implemented.
  • the computing system 700 (system computer) may include one or more system computers 730 , which may be implemented as any conventional personal computer or server.
  • system computers 730 may be implemented as any conventional personal computer or server.
  • HTTP hypertext transfer protocol
  • FIG. 7 illustrates a computing system 700 , into which implementations of various techniques described herein may be implemented.
  • the computing system 700 may include one or more system computers 730 , which may be implemented as any conventional personal computer or server.
  • HTTP hypertext transfer protocol
  • hand-held devices hand-held devices
  • multiprocessor systems multiprocessor systems
  • microprocessor-based or programmable consumer electronics network PCs
  • minicomputers minicomputers
  • mainframe computers and the like.
  • the system computer 730 may be in communication with disk storage devices 729 , 731 , and 733 , which may be external hard disk storage devices. It is contemplated that disk storage devices 729 , 731 , and 733 are conventional hard disk drives, and as such, will be implemented by way of a local area network or by remote access. Of course, while disk storage devices 729 , 731 , and 733 are illustrated as separate devices, a single disk storage device may be used to store any and all of the program instructions, measurement data, and results as desired.
  • measurements received at step 410 in method 400 may be stored in disk storage device 731 .
  • the system computer 730 may retrieve the appropriate data from the disk storage device 731 to predict effective permeabilities according to program instructions that correspond to implementations of various techniques described herein.
  • the program instructions may be written in a computer programming language, such as C++, Java and the like.
  • the program instructions may be stored in a computer-readable medium, such as program disk storage device 733 .
  • Such computer-readable media may include computer storage media and communication media.
  • Computer storage media may include volatile and non-volatile, and removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules or other data.
  • Computer storage media may further include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid state memory technology, CD-ROM, digital versatile disks (DVD), or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the system computer 730 .
  • Communication media may embody computer readable instructions, data structures or other program modules.
  • communication media may include wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above may also be included within the scope of computer readable media.
  • system computer 730 may present output primarily onto graphics display 727 , or alternatively via printer 728 .
  • the system computer 730 may store the results of the methods described above on disk storage 1029 , for later use and further analysis.
  • the keyboard 726 and the pointing device (e.g., a mouse, trackball, or the like) 725 may be provided with the system computer 730 to enable interactive operation.
  • the system computer 730 may be located at a data center remote from the region were the earth formations were obtained from.
  • the system computer 730 may be in communication with the logging device described in FIG. 1 (either directly or via a recording unit, not shown), to receive signals indicating the measurements on the earth formations.
  • These signals after conventional formatting and other initial processing, may be stored by the system computer 730 as digital data in the disk storage 731 for subsequent retrieval and processing in the manner described above.
  • these signals and data may be sent to the system computer 730 directly from sensors, such as well logs and the like.
  • the system computer 730 may be described as part of an in-field data processing system.
  • system computer 730 may process seismic data already stored in the disk storage 731 .
  • system computer 730 may be described as part of a remote data processing center, separate from data acquisition.
  • the system computer 730 may be configured to process data as part of the in-field data processing system, the remote data processing system or a combination thereof. While FIG. 7 illustrates the disk storage 731 as directly connected to the system computer 730 , it is also contemplated that the disk storage device 731 may be accessible through a local area network or by remote access.
  • disk storage devices 729 , 731 are illustrated as separate devices for storing input seismic data and analysis results, the disk storage devices 729 , 731 may be implemented within a single disk drive (either together with or separately from program disk storage device 733 ), or in any other conventional manner as will be fully understood by one of skill in the art having reference to this specification.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Molecular Biology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Toxicology (AREA)
  • Remote Sensing (AREA)
  • Geophysics And Detection Of Objects (AREA)
US12/637,958 2009-06-11 2009-12-15 Estimating effective permeabilities Abandoned US20100313633A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US12/637,958 US20100313633A1 (en) 2009-06-11 2009-12-15 Estimating effective permeabilities
BRPI1001536-1A BRPI1001536A2 (pt) 2009-06-11 2010-05-07 mÉtodo para determinar permeabilidades eficazes de formaÇÕes terrestres
MX2010005336A MX2010005336A (es) 2009-06-11 2010-05-14 Estimacion de permeabilidades efectivas.

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US18621009P 2009-06-11 2009-06-11
US12/637,958 US20100313633A1 (en) 2009-06-11 2009-12-15 Estimating effective permeabilities

Publications (1)

Publication Number Publication Date
US20100313633A1 true US20100313633A1 (en) 2010-12-16

Family

ID=43305209

Family Applications (1)

Application Number Title Priority Date Filing Date
US12/637,958 Abandoned US20100313633A1 (en) 2009-06-11 2009-12-15 Estimating effective permeabilities

Country Status (3)

Country Link
US (1) US20100313633A1 (pt)
BR (1) BRPI1001536A2 (pt)
MX (1) MX2010005336A (pt)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140005945A1 (en) * 2010-12-10 2014-01-02 Schlumberger Technology Corporation Method For Estimation Of Borehole And Formation Properties From Nuclear Logging Measurements
CN103674803A (zh) * 2013-11-18 2014-03-26 中国石油天然气集团公司 一种计算煤岩割理渗透率的方法及系统
US20140214324A1 (en) * 2013-01-25 2014-07-31 Schlumberger Technology Corporation Predicting Mineralogy Properties From Elemental Compositions
US20140257702A1 (en) * 2013-03-08 2014-09-11 Saudi Arabian Oil Company Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data
CN104634804A (zh) * 2013-11-08 2015-05-20 中国石油天然气股份有限公司 一种利用核磁共振t2谱确定储层相对渗透率的方法
US20160201447A1 (en) * 2015-01-14 2016-07-14 Chevron U.S.A. Inc. Hydraulic fracture permeability characterization from electromagnetic measurements
WO2018075019A1 (en) * 2016-10-18 2018-04-26 Halliburton Energy Services, Inc. Relative permeability estimation methods and systems employing downhole pressure transient analysis, saturation analysis, and porosity analysis
US10337315B2 (en) 2015-11-25 2019-07-02 International Business Machines Corporation Methods and apparatus for computing zonal flow rates in reservoir wells
CN110348736A (zh) * 2019-07-10 2019-10-18 中国海洋石油集团有限公司 一种地下水封洞库施工巷道集中渗漏段识别方法
CN110662962A (zh) * 2017-03-14 2020-01-07 沙特阿拉伯石油公司 烃源岩属性的协同感测与预测
US10585083B2 (en) * 2016-06-13 2020-03-10 Schlumberger Technology Corporation System and method for predicting viscosity of heavy oil formations based on nuclear magnetic resonance (NMR) measurements
CN112855089A (zh) * 2021-02-01 2021-05-28 重庆科技学院 一种计算二维填砂模型有效渗透率的应用方法
CN112946006A (zh) * 2021-02-04 2021-06-11 武汉大学 裂隙岩心渗透注浆充填效果检测方法及系统
US11073463B2 (en) * 2018-05-15 2021-07-27 W.D. Von Gonten Laboratories, LLC Methods for modeling permeability in layered rock formations
US11280931B2 (en) 2017-05-09 2022-03-22 Halliburton Energy Services, Inc. Reservoir formation characterization from NMR T1/T2 ratio
US20230259670A1 (en) * 2022-02-16 2023-08-17 Saudi Arabian Oil Comoany Permeability modeling in a reservoir simulation model using dynamic pressure transient analysis
CN116625898A (zh) * 2023-04-03 2023-08-22 内蒙古自治区交通运输科学发展研究院 一种冻融条件下温度与渗流特性确定方法及相关设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111827995B (zh) * 2020-07-31 2023-07-25 中国石油天然气集团有限公司 一种基于核磁共振的超压砾岩储层渗透率的计算方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040093937A1 (en) * 2001-01-18 2004-05-20 Hashem Mohamed Naguib Determining the in situ effective mobility and the effective permeabilty of a formation
US20050257611A1 (en) * 2004-05-21 2005-11-24 Halliburton Energy Services, Inc. Methods and apparatus for measuring formation properties
US20060055403A1 (en) * 2004-04-30 2006-03-16 Schlumberger Technology Corporation Method for determining characteristics of earth formations
US20060155475A1 (en) * 2003-03-27 2006-07-13 Hezhu Yin Method to convert seismic traces into petrophysical property logs
US20080066534A1 (en) * 2006-09-18 2008-03-20 Lennox Reid Obtaining and evaluating downhole samples with a coring tool
US20080133186A1 (en) * 2006-12-04 2008-06-05 Chevron U.S.A. Inc. Method, System and Apparatus for Simulating Fluid Flow in a Fractured Reservoir Utilizing A Combination of Discrete Fracture Networks and Homogenization of Small Fractures
US20090114385A1 (en) * 2007-09-26 2009-05-07 Peter Lumbye Method of stimulating a well
US20090175125A1 (en) * 2008-01-03 2009-07-09 Nonlinear Seismic Imaging, Inc. Direct Mapping of Oil-Saturated Subsurface Formations
US20090255669A1 (en) * 2008-04-15 2009-10-15 Schlumberger Technology Corporation Tool and method for determining formation parameter
US20100198570A1 (en) * 2009-01-30 2010-08-05 Pallav Sarma System and method for predicting fluid flow in subterranean reservoirs
US20100254217A1 (en) * 2009-04-03 2010-10-07 Dezhi Chu Method for Determining the Fluid/Pressure Distribution of Hydrocarbon Reservoirs from 4D Seismic Data

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040093937A1 (en) * 2001-01-18 2004-05-20 Hashem Mohamed Naguib Determining the in situ effective mobility and the effective permeabilty of a formation
US20060155475A1 (en) * 2003-03-27 2006-07-13 Hezhu Yin Method to convert seismic traces into petrophysical property logs
US20060055403A1 (en) * 2004-04-30 2006-03-16 Schlumberger Technology Corporation Method for determining characteristics of earth formations
US20050257611A1 (en) * 2004-05-21 2005-11-24 Halliburton Energy Services, Inc. Methods and apparatus for measuring formation properties
US20080066534A1 (en) * 2006-09-18 2008-03-20 Lennox Reid Obtaining and evaluating downhole samples with a coring tool
US20080133186A1 (en) * 2006-12-04 2008-06-05 Chevron U.S.A. Inc. Method, System and Apparatus for Simulating Fluid Flow in a Fractured Reservoir Utilizing A Combination of Discrete Fracture Networks and Homogenization of Small Fractures
US20090114385A1 (en) * 2007-09-26 2009-05-07 Peter Lumbye Method of stimulating a well
US20090175125A1 (en) * 2008-01-03 2009-07-09 Nonlinear Seismic Imaging, Inc. Direct Mapping of Oil-Saturated Subsurface Formations
US20090255669A1 (en) * 2008-04-15 2009-10-15 Schlumberger Technology Corporation Tool and method for determining formation parameter
US20100198570A1 (en) * 2009-01-30 2010-08-05 Pallav Sarma System and method for predicting fluid flow in subterranean reservoirs
US20100254217A1 (en) * 2009-04-03 2010-10-07 Dezhi Chu Method for Determining the Fluid/Pressure Distribution of Hydrocarbon Reservoirs from 4D Seismic Data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Definition of "permeability", Free Online Dictionary, 2012. *

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140005945A1 (en) * 2010-12-10 2014-01-02 Schlumberger Technology Corporation Method For Estimation Of Borehole And Formation Properties From Nuclear Logging Measurements
US9696453B2 (en) * 2013-01-25 2017-07-04 Schlumberger Technology Corporation Predicting mineralogy properties from elemental compositions
US20140214324A1 (en) * 2013-01-25 2014-07-31 Schlumberger Technology Corporation Predicting Mineralogy Properties From Elemental Compositions
US9671525B2 (en) * 2013-03-08 2017-06-06 Saudi Arabian Oil Company Determining continuous capillary pressure curves for subsurface earth formations using saturation and NMR log data
US20140257702A1 (en) * 2013-03-08 2014-09-11 Saudi Arabian Oil Company Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data
WO2014137863A3 (en) * 2013-03-08 2015-01-08 Saudi Arabian Oil Company Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data
CN105190366A (zh) * 2013-03-08 2015-12-23 沙特阿拉伯石油公司 使用饱和度及核磁共振测井数据确定地表下地层的连续毛细管压力曲线
CN104634804A (zh) * 2013-11-08 2015-05-20 中国石油天然气股份有限公司 一种利用核磁共振t2谱确定储层相对渗透率的方法
CN103674803A (zh) * 2013-11-18 2014-03-26 中国石油天然气集团公司 一种计算煤岩割理渗透率的方法及系统
US20160201447A1 (en) * 2015-01-14 2016-07-14 Chevron U.S.A. Inc. Hydraulic fracture permeability characterization from electromagnetic measurements
US9938818B2 (en) * 2015-01-14 2018-04-10 Chevron U.S.A. Inc. Hydraulic fracture permeability characterization from electromagnetic measurements
WO2016114896A1 (en) * 2015-01-14 2016-07-21 Chevron U.S.A. Inc. Hydraulic fracture permeability characterization from electromagnetic measurements
US10337315B2 (en) 2015-11-25 2019-07-02 International Business Machines Corporation Methods and apparatus for computing zonal flow rates in reservoir wells
US10585083B2 (en) * 2016-06-13 2020-03-10 Schlumberger Technology Corporation System and method for predicting viscosity of heavy oil formations based on nuclear magnetic resonance (NMR) measurements
WO2018075019A1 (en) * 2016-10-18 2018-04-26 Halliburton Energy Services, Inc. Relative permeability estimation methods and systems employing downhole pressure transient analysis, saturation analysis, and porosity analysis
CN110662962A (zh) * 2017-03-14 2020-01-07 沙特阿拉伯石油公司 烃源岩属性的协同感测与预测
US11280931B2 (en) 2017-05-09 2022-03-22 Halliburton Energy Services, Inc. Reservoir formation characterization from NMR T1/T2 ratio
US11073463B2 (en) * 2018-05-15 2021-07-27 W.D. Von Gonten Laboratories, LLC Methods for modeling permeability in layered rock formations
CN110348736A (zh) * 2019-07-10 2019-10-18 中国海洋石油集团有限公司 一种地下水封洞库施工巷道集中渗漏段识别方法
CN112855089A (zh) * 2021-02-01 2021-05-28 重庆科技学院 一种计算二维填砂模型有效渗透率的应用方法
CN112946006A (zh) * 2021-02-04 2021-06-11 武汉大学 裂隙岩心渗透注浆充填效果检测方法及系统
US20230259670A1 (en) * 2022-02-16 2023-08-17 Saudi Arabian Oil Comoany Permeability modeling in a reservoir simulation model using dynamic pressure transient analysis
CN116625898A (zh) * 2023-04-03 2023-08-22 内蒙古自治区交通运输科学发展研究院 一种冻融条件下温度与渗流特性确定方法及相关设备

Also Published As

Publication number Publication date
BRPI1001536A2 (pt) 2011-07-26
MX2010005336A (es) 2010-12-13

Similar Documents

Publication Publication Date Title
US20100313633A1 (en) Estimating effective permeabilities
US9696453B2 (en) Predicting mineralogy properties from elemental compositions
US10445669B2 (en) System and method for mapping reservoir properties away from the wellbore
Heidari et al. Improved estimation of mineral and fluid volumetric concentrations from well logs in thinly bedded and invaded formations
Chang et al. Evaluation and comparison of residual oil saturation determination techniques
US7555390B2 (en) Petrophysical interpretation of multipass array resistivity data obtained while drilling
Mondol Well logging: Principles, applications and uncertainties
US10495780B2 (en) Correcting shale volume and measuring anisotropy in invaded zone
CA2546139A1 (en) Method for determining characteristics of earth formations
US20210255358A1 (en) Method for estimating the subsurface total organic carbon (TOC) from well-log data
US10151197B2 (en) Hydrocarbon density determination method
Yao et al. Reservoir permeability estimation from time-lapse log data
US11788402B2 (en) Methods for determining a volume fraction of water and water salinity in organic shale reservoirs using multi-frequency dielectric well logging measurements
Bennett et al. Use of geostatistical modeling and automatic history matching to estimate production forecast uncertainty–a case study
US20230288604A1 (en) Hydrocarbon Reservoir Saturation Logging
Davies Permeability Modelling of a Sandstone Reservoir in Parts of the Niger Delta
Biswal et al. Petrophysical modeling in the absence of core data in a small, marginal field in the Mumbai Offshore Basin, India
Ghosh et al. Petrophysical Characterization of Mineralogically Complex Panna Formation: An Alternate Approach
Wiki Standard interpretation
Cavalleri et al. Rock Permeability Conundrum-Workflow for Continuous Quantitative Assessment as Input to Real-Time Interpretation for Producibility and Well Deliverability at the Early Stage
Nguyen Pham et al. Unlocking Potential in a Thinly Laminated Reservoir by Integrating Electric Logs, Images and Core Analysis: First Case Study in Vietnam
Podberezhny et al. Determination of Invaded Zone Properties from Time-Lapse Induction Logging Measurements Based on Invasion Simulation and Electromagnetic Modeling
Frenkel et al. Integration of LWD and wireline array technologies to improve estimation of hydrocarbon reserves
Raymer et al. The role of well logs in reservoir modeling
Cheng et al. Improvement in Heavy-Oil Reservoir Evaluation Using Nuclear Magnetic Resonance: Long Lake and Kinosis SAGD Projects, Alberta, Canada

Legal Events

Date Code Title Description
AS Assignment

Owner name: SCHLUMBERGER TECHNOLOGY CORPORATION, TEXAS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ANAND, VIVEK;FREEDMAN, ROBERT;REEL/FRAME:023686/0093

Effective date: 20091218

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION