WO2009108432A1 - Rock physics model for simulating seismic response in layered fractured rocks - Google Patents
Rock physics model for simulating seismic response in layered fractured rocks Download PDFInfo
- Publication number
- WO2009108432A1 WO2009108432A1 PCT/US2009/032007 US2009032007W WO2009108432A1 WO 2009108432 A1 WO2009108432 A1 WO 2009108432A1 US 2009032007 W US2009032007 W US 2009032007W WO 2009108432 A1 WO2009108432 A1 WO 2009108432A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- wave
- anisotropy
- model
- subsurface region
- anisotropic
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
Definitions
- This invention relates generally to the field of petroleum reservoir management or other fields where movement of pore fluids through fractured rock is important, and more particularly to seismic fracture characterization methods.
- the invention is a method for developing a rock physics model for layered fractured rocks to use in simulating seismic response.
- Fractures may serve as major conduits for movement of pore fluids (hydrocarbons or water) and dissolved chemicals through low porosity, low permeability reservoir or non- reservoir rocks. Understanding fluid flow and mass transport in fractured rocks is essential for optimal reservoir management as well as other applications such as assessing the ground- water resources of hard-rock aquifers, investigating the suitability of underground sites for hazardous waste disposal, and predicting the movement of hazardous chemicals if contamination occurs.
- a key strategy for fractured reservoir management is a quantitative description of the geology, geophysics and petrophysical attributes.
- aligned fractures will cause the anisotropic behavior of rock properties (such as wave velocities, attenuation, resistivity, thermal conductivity and permeability).
- rock properties such as wave velocities, attenuation, resistivity, thermal conductivity and permeability.
- Different fracture configurations/alignments will have different types of anisotropy.
- a resulting effective medium of an isotropic rock matrix permeated with a single set of aligned fractures will be transversely isotropic (TI).
- TI transversely isotropic
- two or three orthogonally-intersected fracture sets will give an orthorhombic system. In more general case where two or more fracture sets intercept at arbitrary angles, one has a monoclinic system.
- Quantifying fracture anisotropy with surface seismic data should provide an optimal strategy for fractured reservoir management by integrating the geophysical data from all scales with the engineering data.
- P-wave AVO/ A VA and azimuthal velocity anisotropy analyses can be important for inferring fracture properties.
- S- wave splitting (birefringence) of either pure S waves or converted PS waves can be analyzed to characterize special fracture distributions.
- a fractured rock layer often has limited thickness and coexists with other rock types, such as shale ( Figure 1).
- the wavelength can be as large as 100 meters. If the overall thickness of the mixed rock column is less than, for example, half of the wavelength, it is not possible to separate the fracture anisotropy from shale anisotropy and layering anisotropy. In practical applications, therefore, it is critical to have a rock physics model that can consistently consider various types of geological factors including, but not limited to, porosity, shale volume, fluid content, lithology, pore type and fractures when modeling seismic response.
- the invention is a method for modeling anisotropic elastic properties of a subsurface region comprising mixed fractured rocks and other geological bodies, in order to predict fluid flow or mass transport or for seismic interpretation in the subsurface region, comprising: (a) developing an anisotropic rock physics model for the subsurface region, wherein the model simulates seismic anisotropy via preferred orientations of clay-related pores, fracture anisotropy via alignment of fractures, and stress-induced anisotropy via preferred orientations of microcracks; (b) calculating an elastic stiffness/compliance tensor, comprising at least three anisotropy parameters, at selected discrete depth intervals using the anisotropic rock physics model;
- the upscaled anisotropic ID model is calibrated before generating the synthetic seismic response.
- the calibration is preferably performed using a calibration method capable of handling general anisotropy.
- the inventive method further comprises:
- the anisotropy parameters are calculated by steps comprising: (a) obtaining P-wave and cross-dipole (fast and slow) 5-wave well log data from a single deviated well in the subsurface region;
- This physics-based method for calculating the anisotropy parameters may be used for any purpose in which anisotropy parameters are needed and deviated well data are available, and the method is particularly advantageous where data are available from only a single deviated well.
- Fig. 1 illustrates a multi- layer system in which fracture layers coexist with non- fracture layers, but the two types of layers have very different anisotropy symmetry;
- FIG. 2 is a flow chart showing basic steps in one embodiment of the present inventive method;
- Fig. 3 illustrates data preparation and parameter estimation for anisotropic rock physics modeling;
- FIG. 4 is a schematic diagram illustrating construction of a rock physics model capable of treating fracture anisotropy, stress-induced anisotropy and shale anisotropy;
- Fig. 5 is a flow chart showing basic steps in an embodiment for constructing a rock physics model corresponding to the schematic diagram of Fig. 4;
- Fig. 6 is a workflow for anisotropy parameter estimation from limited well logs at deviated wells and seismic modeling;
- Fig. 7 illustrates a single fracture set, aligned perpendicular to the x-axis, in a mixed carbonate-clastic environment
- Fig. 8 illustrates two fracture sets, one aligned perpendicular to the x-axis and the other perpendicular to the y-axis, in a mixed carbonate-clastic environment
- Fig.9 illustrates P-, SV- and SH-wave velocities as a function of incident (deviation) angles
- Fig. 10 compares well log shear-wave splitting (the difference between the fast and slow S-wave logs) and that calculated using our rock physics model at two North Sea wells;
- Fig. 11 illustrates the upscaling results in VTI layered media using Backus Averaging
- Fig. 12 illustrates the effect of sonic log correction and anisotropy on AVO response using data from Fig.11.
- the present invention is a method for seismic modeling in a mixed geological environment where fractured layers are inter-bedded with non- fractured layers, an example of which is illustrated in Fig. 1.
- Layers 1 and 3 are shale
- 2 and 4 are fracture layers.
- the anisotropy symmetry in the fractured layers (typically HTI for a single fracture set) can be significantly different from that in the non-fractured layers (typically VTI).
- the thinness of the layers is indicated by comparison to the approximate seismic wavelength as shown by wavelet 5.
- Figure 2 is a flowchart showing basic steps in one embodiment of the present inventive method for accurate well/seismic tie and AVO (amplitude versus offset) modeling.
- step 21 parameters and data, which are necessary for anisotropy rock physics modeling, are collected and analyzed. This involves the following sub-steps, outlined in Fig.
- step 31 deviation survey data and bedding dipping data, if available, are obtained for a subsurface region, and typically are loaded into a computer system. The data are needed to define the relative angle between the well trajectory and the normal to bedding.
- Bedding dipping data may be interpreted from seismic data or well log data. If deviation survey data and bedding dipping data are not available, assumptions may be made on the relative angle between the well trajectory and bedding normal; for example, a vertical well with horizontal bedding may be assumed.
- step 32 temperature and pressure profiles are constructed to consider their effects on fluid properties.
- the stress data are also useful for the analysis of stress-induced anisotropy, especially in young unconsolidated sands.
- a plurality of rock properties e.g. porosity, clay content, mineralogy, permeability, and water saturation are estimated from a plurality of well logs such as gamma ray, density, neutron porosity, sonic, and various resistivity (deep, medium, shallow and micro) logs.
- fracture parameters such as fracture density and fracture orientation distribution are estimated from, for example, FMI (Formation Micro-Imager) and/or cross- dipole logs.
- the orientation distribution for clay-related pores may be estimated from high- resolution images, such as SEM (Scanning Electron Microscopy). Crack orientation distribution parameters may be calibrated using laboratory or/and well log data (Xu 2002).
- step 22 the elastic stiffness/compliance tensor at each depth is calculated using a rock physics model that is able to consistently handle different types of seismic anisotropy in geo-layers with or without fractures.
- a rock physics model will be described below.
- Results from steps 31-34 are used as input quantities in developing the rock physics model.
- P-wave logs, and fast and slow S-wave logs traveling in the direction along the wellbore are then calculated from the elastic stiffness/compliance tensor.
- step 23 the calculated P- and 5-wave logs may then be compared with measured data, if available. If their agreement is not satisfactory, return to step 21 or 22 to reevaluate the input parameters to the rock physics model and repeat the calculation again. Otherwise, continue to step 24 below.
- anisotropic upscaling capable of handling different types of anisotropy (e.g. VTI, HTI and orthorhombic) is applied to the calculated well logs from step 22 or 23, sampled at original fine intervals, e.g. 0.5 ft., for efficient seismic modeling.
- Step 24 involves two sub-steps: (a) generate a ID layered earth model, typically by detecting boundaries of geological layers from one or more well logs (e.g. P-, 5-wave and/or density logs); and
- checkshot data may be applied to calibrate the anisotropic earth model (the calibration step is recommended but not essential to the invention) and generate synthetic seismic response.
- (b) Generate synthetic CMP (common middle point) gathers using a program that is able to handle general anisotropy (e.g. orthorhombic) and perform AVO analyses.
- Tying seismic data to well logs is one of the main hydrocarbon reservoir characterization work processes to evaluate the quality of the seismic data, the well log data, and the validity of the rock physics model. Seismic data are in some respects meaningless until tied to well logs, which link to the geology. Data interpreters, for example, need seismic/well tie(s) to define the top or base reservoir on seismic data. Quality of seismic data is evaluated using goodness in the match between synthetic seismic data and the actual (field) seismic data.
- AVO analysis is typically used for fluid prediction (e.g., oil or water) and lithology discrimination. Synthetic seismic gathers are typically compared with real seismic gathers to quantify the fluid effect or lithology effect. However, AVO is affected by other factors, anisotropy in particular. These factors need to be considered in order to interpret the observed AVO anomalies correctly.
- FIG. 4 is a schematic diagram showing how the model works in practice.
- the sand grains are mixed with clay particles using a mixing law (e.g. the Voigt-Reuss-Hill average).
- the various pores are added using an anisotropic effective medium theory.
- the clay pores are added to the system using the differential effective medium (“DEM”) scheme or the anisotropic dry rock approximation (“ADRA”) scheme, which will be described later, to consider the mechanical interaction between the pores.
- DEM differential effective medium
- ADRA anisotropic dry rock approximation
- step 43 is elastic constants for the region of interest that are generated with consideration of not only stress-induced anisotropy and shale anisotropy but also fracture anisotropy.
- P-, and S-wave velocities and Thomsen parameters can be calculated from the effective elastic constants obtained using the above procedures. Development of the model in one embodiment of the invention is described in more detail next, with reference to the flow chart of Fig. 5 (Fig. 4 reference numbers shown in parentheses where applicable).
- Step 51 define pore types and partition pore space based on porosity, shale volume, stress and fracture data evaluated in steps 33 and 34 of Fig. 3.
- the pore system has five pore types: inter-particle, equant, clay pores, microcracks, and fractures. Treating more than five pore types is straightforward. Then, different minerals are mixed using a mixing law based on the mineralogy data evaluated in step 33 of Fig. 3.
- Step 52 (41) to consider the mineralogy effect, a mixing law, such as the Hill Average, is used to average the bulk and shear moduli of the constituent minerals, e.g. quartz, calcite, dolomite and clay particles.
- Step 53 (42) define which pores are communicating and which pores are effectively isolated according their pore sizes (local permeability) to consider the effect of fluid movement on P- and S-wave velocities at the pore scale. In general, clay pores and micro cracks which are smaller than, say, 1 ⁇ are practically isolated and always occupied by bound water.
- Step 54 (42) add water-wet micro pores (microcracks and clay pores) into the system using a theory (e.g. DEM or ADRA). This will give the effective compliance tensor of the composite with water- wet micro pores only. The effective compliance tensor is needed for anisotropic Gassmann fluid substitution later.
- a theory e.g. DEM or ADRA.
- Step 55 (42) add all pores (micro cracks, clay pores, inter-particle, equant pores and fractures) into the system using a theory.
- micro pores e.g. clay pores and microcracks
- These micro pores with fluid added at this point will give a high-frequency seismic response since the fluid contained in these pores is not allowed to communicate with fluids in other pores.
- Macro pores inter- particle, equant pores and fractures are added into the system without fluid at this point. This will give the effective compliance tensor of the "dry" rock.
- Step 56 (43) determine and mix the rest of the fluids using a fluid mixing law.
- Step 57 (43) put the mixed fluid into the pores without fluid using an anisotropic Gassmann fluid substitution method (Brown and Korringa 1977). This will give the final effective stiffness/compliance tensor of the saturated rock.
- a deviated well means a non-vertical well, but the method described next also works when the well is vertical but the layers in the earth are dipping.
- the term deviated well will be understood herein, including the claims, to include both situations.
- the following method is advantageous either as a stand-alone method for any application in which anisotropy parameters are needed or as part of the process for updating the rock physics model in step 23 of Fig. 2 when deviated well data are available.
- the method is physics-based rather than empirical.
- the result is an estimate of the anisotropy parameters ( ⁇ and ⁇ ) defined by Thomsen (1986) and the vertical P-wave velocity.
- the authors also proposed a procedure to estimate the anisotropy-corrected vertical sonic logs from sonic data recorded in a deviated well.
- the inputs are well deviation, P-wave sonic log, volume of shale or gamma ray data, and anisotropy parameters for rock with 100% shale volume.
- the compressional-wave group velocity surface is computed and the equivalent vertical P-wave sonic log is output.
- the equivalent vertical sonic log was used for standard seismic applications using isotropic velocity assumptions.
- a model-based (rather than empirical) method is disclosed herein that is highly constrained by rock physics principles.
- the estimated anisotropy parameters may then be validated using cross-dipole logs, by comparing the observed fast and slow shear-wave logs (hence the shear-wave splitting) at the well and that calculated from the model.
- the method is very flexible and can be applied to a single or multi-deviated well as long as cross-dipole logs are available.
- the Thomsen parameters ( ⁇ , ⁇ and ⁇ ) calculated using this method are internally consistent because they are calculated using the same method and same parameters (i.e. input data).
- step 61 collect and analyze parameters and data for an anisotropic rock physics model. This step is very similar to step 31 of Fig. 3, except this method is applicable to rocks with or without fractures. Therefore, fracture data may or may not be required depending on the technical issue to be addressed.
- step 62 determine working interval and select a starting depth.
- step 63 calculate full elastic constants from which Thomsen's (1986) anisotropy parameters ( ⁇ , ⁇ and ⁇ 5) are then derived, the elastic constants being calculated from the input parameters using an anisotropic rock physics model (e.g. Xu et al. 2006).
- P- wave, cross-dipole (fast and slow) 5-wave and density logs along the wellbore are calculated from the elastic tensors.
- step 64 compare measured cross-dipole (fast and slow) S-wave logs with the calculated to quantify the accuracy of the calculated anisotropy parameters.
- Step 64 is a significant step for the anisotropy parameter calibration using measured P-wave, fast and slow S-wave logs. This is because shear-wave splitting (the difference between SH-wave and SV-wave logs) is determined by the anisotropy parameters. Therefore, the goodness of fit between the calculated shear-wave splitting and that derived from the fast and slow S-wave logs indicates the accuracy of the calculated anisotropy parameters.
- step 65 apply an appropriate upscaling method to block the P-wave, S-wave, density and the calculated anisotropy parameter logs consistently.
- This step is very similar to step 24 except that fractures are optional here.
- the Backus Averaging Backus
- upscaling may be used if all the layers are transversely isotropic with a vertical symmetry axis (VTI) or isotropic (no fractures).
- VTI vertical symmetry axis
- isotropic no fractures.
- Schoenberg and Muir, (1989) will give good results.
- checkshot data may be used to calibrate the anisotropic earth model and generate synthetic seismic response. Again, this step is very similar to step 25, except the earth model can be VTI only, or mixed VTI with HTI or orthorhombic.
- the rock physics model of the present invention can predict seismic anisotropy of such a mixed system.
- the model is flexible enough to handle more general types of anisotropy (e.g. orthorhombic or monoclinic).
- the model is based on, but is developed beyond, the anisotropic dry rock approximation proposed by Xu, Saltzer and Keys (2006). It is, therefore, accurate and efficient.
- the present invention provides (1) a five- component rock physics model which is able to consistently predict elastic properties of different rock types with or without fractures, (2) a workflow to consistently evaluate parameters for such a model and correctly handle the layering effect (upscaling), and (3) a consistent method for checkshot calibration for rocks with arbitrary anisotropy.
- Pore partitioning [0063] A mixed geological environment with fractured rock layers coexisting with non- fractured rock types such as shale is assumed.
- One critical part of the present invention is the definition of the pore system.
- the invention uses effective medium theory and idealizes pores as ellipsoids.
- the pore shape is therefore characterized by the aspect ratio (short axis/long axis) of the ellipsoid.
- effective medium theory there is no limitation on the number of pore types (shapes) in our method.
- the present invention is not limited to the five pore types mentioned above.
- the rock physics model is able to handle a pore system with more than five pore types. How many pore types are used in the model primarily depends on the objectives of the work and the type of data. Also, not all the pore types may exist in a single rock type. For example, equant pores are frequently referred to as the modic or vuggy pores in carbonate rocks. In a clastic shaly-sand system, ⁇ Eqmnt may be set to zero. Each pore type can be characterized using some statistical parameters, such as pore orientation, pore size and pore shape (aspect ratio) and/or their distributions. [0065] Hornby et al.
- ⁇ ciay V S h ⁇ ( 2 ) where V sh is the shale volume normalized by grain matrix.
- V sh is the shale volume normalized by grain matrix.
- the anisotropy in shale is simulated via a preferred orientation distribution of the clay-related pores, which can be described by, for example, a Gaussian distribution.
- the amount of microcracks is often related to the in- situ stress conditions. This is how stress-induced anisotropy may be simulated in clastic rocks, such as unconsolidated sandstones. In tectonically passive regions, stress conditions are dominated by compaction.
- the model of the present invention can handle more than one fracture set.
- the fracture system has two orthogonally-intersected fracture sets.
- the coordinate system may be chosen so that one fracture set is aligned perpendicular to x-axis and the other to y-axis. Fractures are seldom perfectly aligned in practice and assuming perfectly aligned fractures often causes an overestimate of fracture anisotropy.
- a Gaussian distribution is again used to describe fracture orientation distribution for each fracture set.
- fracture density estimated from FMI logs preferably are calibrated using cross-dipole logs or multi-component VSP (vertical seismic profiling) data, or multi-azimuth seismic data if available.
- Equant pores may exist in some special rocks e.g. carbonates, volcanics. They typically represent modic or vuggy pores in carbonate rocks. Photographical data, such as, thin sections, 3D X-ray microtomography images, or SEM pictures, are needed to determine amount of equant pores in the system. The rest of porosity will be inter-granular pores, typically the most commonly seen pore type in sedimentary rocks. In this specific embodiment of the invention, both equant pores and inter-granular pores are assumed to be randomly orientation, giving isotropic seismic response.
- the model simulates seismic anisotropy via the preferred orientations of clay-related pores, fracture anisotropy via the alignment of fractures, and stress-induced anisotropy via the preferred orientation of microcracks, which is controlled by the differential stress.
- V 1 denotes the volume fraction for the f mineral
- C 1 is its stiffness tensor
- C M is the stiffness tensor of the mixed minerals.
- Nishizawa (1982) showed how to calculate the effective elastic stiffness tensor for a solid embedded with parallel pores.
- Nishizawa' s theory is extended to the case where the pores may have a preferred orientation distribution, and the differential equations are formulated according to the approach suggested by Hornby et al. (1994). Following Xu et al. (2006), one gets
- a is the aspect ratio of a pore type, with a porosity of ⁇ a .
- ⁇ and ⁇ are polar angle and azimuth, respectively.
- A is the compliance tensor of the effective medium and ⁇ is porosity.
- F a ( ⁇ , ⁇ ) is the probability density function for the orientation distribution for each pore type.
- Tensor Q is calculated in the local coordinate system, where its z-axis is chosen parallel to the symmetry axis of the inclusion.
- C and C 1 are the elastic stiffness tensors of the solid and fluid phases, respectively.
- 5 is the Eshelby (1957) tensor, which is a function of the pore shape (aspect ratio). As shown later, S can be calculated using equations given by Eshelby (1957) when the background solid is isotropic and using equations given by Mura (1987) when the background solid is transversely isotropic.
- equation (5) has to be solved iteratively (numerically). The number of iterations depends on porosity and the pore aspect ratios.
- / is the identity tensor.
- S tensor is approximately independent of porosity in a dry rock and, therefore, remains constant.
- a 0 is the compliance tensor of the solid matrix.
- P is a fourth-order tensor rather than a scalar. The anisotropic dry rock approximation dramatically increases the calculation speed.
- Equation (11) The compliance tensor in Equation (11) is the anisotropic dry rock approximation.
- the tensor (l - ⁇ ) P in Equation (11) can be defined by the power series expansion:
- the tensor P in Equation (12) is isomorphic to a 9 x 9 matrix. That is, there is a one-to-one transformation from the space of tensors of rank four onto the space of 9 x 9 matrices that preserves multiplication and addition. Consequently, the tensor P in Equation (12) can be replaced by a 9 x 9 matrix P, so that the power series in (12) is defined by ordinary matrix operations.
- the differential effective medium (DEM) theory or the Anisotropic Dry Rock Approximation (ADRA) described above may then be used to add the pore space into the system. Special attention is preferably paid to the fluid communication issue.
- the present inventive method is able to simulate elastic properties of the composite with relaxed (all the pore fluid are communicating), unrelaxed (no fluid communication, i.e. isolated pores), or partially relaxed pore fluids. Which pore fluid is communicating and which is not depends on many factors such as pore size, pore connectivity (local permeability), fluid viscosity and frequency.
- One embodiment of the present invention employs a four-step procedure to implement certain concepts (see Fig. 4).
- the effective compliance tensor (A 0 ) is calculated for the solid with water-wet micro pores (clay pores and microcracks) only (step
- step 45 the effective compliance tensor (A D ) is calculated for the solid with all the pores (water-wet for clay pores and microcracks and dry for inter-particle, equant pores and fractures).
- step 46 the effective water saturation and the effective fluid properties are calculated for the dry pores using the following equations (step 46):
- a mixing law e.g. the Wood suspension model, may be used to calculate the compressibility (c f ) of the pore fluid mixture: where C Brme and C Hydrocarbon are compressibilities of brine and hydrocarbon (oil or gas), respectively.
- sonic logs are overly sampled; sampling intervals of the sonic logs are usually much smaller (e.g., 0.5 ft) than the seismic wavelength (e.g., 200 ft).
- Synthetic seismic modeling performed on the original sampled sonic logs is computationally expensive and unnecessary.
- Log blocking/upscaling is commonly carried out to "up-sample” the sonic logs to allow for realistic computation of synthetic seismogram generation.
- Special treatment of sonic log blocking/upscaling is a feature of the present invention.
- Backus Averaging Backus 1962] is used to block sonic logs.
- each layer can be isotropic, vertically transversely isotropic (VTI), horizontally transversely isotropic (HTI) or orthorhombic.
- C 1 are the elastic constants of each individual layer and are the effective elastic constants of the equivalent medium.
- ⁇ denotes thickness-weighted arithmetic averaging.
- the preceding approach is one embodiment of the invention; it handles one vertically aligned fracture set or two orthogonally intersected fracture sets.
- a method such as that given by Schoenberg and Muir, (1989) may be used.
- Checkshot calibration in an anisotropic media is tricky. Traditionally, checkshot data are applied to modify a velocity (P-wave or S-wave) log only. This traditional approach does not maintain the internal consistency of the elastic constants in the cases of anisotropic media. In those embodiments of the invention using checkshot calibration, i.e. calibration with checkshot data, a four-step approach is preferred to maintain internal consistency:
- V MoM is the P- or S-wave velocity of the layered model in the direction of calibration
- V CheckShotl is the P- or 5-wave velocity derived from the checkshot data for the same layer.
- the Eshelby tensor S l]U describes the pore geometry effect on the overall effective elastic properties of the porous rock. This is an important parameter which controls the porosity-velocity relationship and, especially, anisotropy. Adding a particular pore into the system means modifying the effective properties of the porous composite in terms of the Eshelby tensor (step 42 in Fig. 4).
- a key step for the differential effective medium theory or the anisotropic dry rock approximation is the calculation of the S tensor.
- Mura (1987) has shown how to calculate the Eshelby tensor S a for a pore embedded in a solid with hexagonal symmetry. where C is the stiffness tensor of the solid.
- G is a fourth-rank geometry tensor which depends on the pore shape and the elastic properties of the solid.
- the 9 missing non-zero elements are: G 2121 , G 1221 , G 2112 , G 3131 , G 1331 , G 3113 , G 3232 , G 2332 , G 2332 with
- Equations 29 to 46 are valid for calculating the geometrical properties of an inclusion in a transversely isotropic (TI) medium. Since the integrals have to be done numerically, it can be extremely time-consuming to calculate S llU , especially when aspect ratios are small, e.g. less than 10 . Since the anisotropic dry rock approximation has analytical solution, S M does not have to be updated and, hence, can be calculated using Eshelby's equations for the isotropic matrix (Eshelby 1957),
- v Poisson's ratio which can be calculated using the following equation
- ⁇ and ⁇ are lame constants for an isotropic solid. Calculate P- andS-wave slowness at an arbitrary angle.
- the invention tries to minimize the difference between the calculated P-wave slowness along the wellbore and the measured P-wave log to invert porosity or other parameters.
- Full elastic constants C ljkl are calculated using the rock physics model for the given porosity, shale volume and water saturation, in addition to other parameters. Then the P-wave slowness in the direction of the wellbore will be calculated. In the case of a transversely isotropic (TI) medium, for example, there are five independent elastic constants.
- TI transversely isotropic
- P-, SV- and SH-wave phase velocities can be calculated using the following equations (Thomsen 1986).
- ⁇ is the angle between the symmetry axis of the TI medium and the wave propagation direction. In the case of a VTI medium, ⁇ is the polar angle.
- the group velocity (v Gmup ) for a particular wave (P, SV or SH) can be calculated from its phase velocity (V Phase ) using the following equation (Berryman 1979).
- ⁇ is the group angle which is determined by the following equation.
- a standard inversion method is used to invert porosity by minimizing the error between the calculated P-wave slowness (DT), as shown above, and the measured P-wave log at each depth.
- the anisotropy is HTI in clean limestone and VTI in pure shale, and orthorhombic in shaly (muddy) carbonate with vertical fractures.
- Azimual anisotropy is indicated by the difference between Cu and C 22 in track 73.
- Very strong azimuthal anisotropy can be seen in a clean limestone layer just below 4000 ft as seen from the large differences between Cn and C 22 (track 73) and that between Cn and C33 (track 74). Note that the nearly identical C22 and C33 (track 75) in the clean limestone zones, indicating a HTI symmetry.
- a second example (Fig. 8) has two fracture sets, a first set aligned perpendicular to the x-axis (the same fracture set as in Fig. 7) and the other perpendicular to y-axis, in the same mixed carbonate-clastic environment (Track 81 in Fig. 8 is the same as track 71 in Fig.
- the fracture porosity for the first fracture set is twice as great as that for the second fracture set, denoted by different shadings of gray in track 82.
- the anisotropy is VTI in pure shale and orthorhombic symmetry in muddy or clean carbonate.
- the azimuthal anisotropy (tracks 83 and 84) in the clean limestone layer just below 4000 ft is much weaker than that for the single fracture case.
- Curve P represents the P-wave velocity
- curve SH represents the horizontal 5-wave velocity
- curve SV represents the vertical 5-wave velocity.
- rock properties can be modeled using a VTI (transverse isotropy with a vertical symmetry axis) system (Xu, Saltzer and Keys 2006).
- shear-wave splitting (the separation between SH- and SF-wave logs) is controlled by anisotropy parameters. If the anisotropic rock physics correctly predicts the anisotropy parameters, the calculated fast and slow 5-wave logs should match the measured reasonably well.
- Fig. 10 shows the results (left three tracks for Well “A”, and right three tracks for Well “B”).
- the dark and light-shaded curves in track 101 are the measured Well “A” cross- dipole slow (typically, but not always, the SV wave) and fast (typically, but not always, the SH wave) shear-wave logs, respectively.
- Their separation is an indicator of shear-wave splitting, which is controlled by anisotropy parameters and the well deviation angle. At 35 degrees (Well "A"), their separation is small. This does not mean weak anisotropy at the well, but rather the deviation angle is not high enough to see the anisotropy.
- Vp-Vs relationship is also affected by anisotropy parameters, in addition to lithology and fluid content.
- VpIVs varies with deviation angles. Since S- wave logs are predicted from P- wave logs using the rock physics model, the goodness of fit between measured (light shade) and calculated (dark shade) fast S-wave logs (Tracks 103 and 106) can also be used to validate the anisotropic rock physics model.
- Fig. 11 shows vertical P- (Track 111), vertical S- wave (Track 112), density (Track 113) logs and the calculated ⁇ (Track 114), ⁇ (Track 115) and ⁇ (Track 116) before (light shade) and after (dark shade) upscaling.
- Fig. 12 shows synthetic gathers generated from well logs illustrated in Fig. 11 to demonstrate the effect of anisotropy on AVO response.
- the synthetic gather on the left was generated from an isotropic earth model, which was constructed using measured P-, slow S- wave and density logs.
- the synthetic gather on the right was generated from the anisotropic earth model which was constructed from full elastic tensors (honoring both P- and S-wave logs). The two synthetic gathers were aligned to the top of the logs.
- the anisotropy effect on travel time can be observed; see, for example, reflection event 121 as depicted in each gather.
Abstract
Description
Claims
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/863,867 US8184502B2 (en) | 2008-02-28 | 2009-01-26 | Rock physics model for simulating seismic response in layered fractured rocks |
CA2710607A CA2710607A1 (en) | 2008-02-28 | 2009-01-26 | Rock physics model for simulating seismic response in layered fractured rocks |
EP09714967A EP2248007A1 (en) | 2008-02-28 | 2009-01-26 | Rock physics model for simulating seismic response in layered fractured rocks |
AU2009217648A AU2009217648A1 (en) | 2008-02-28 | 2009-01-26 | Rock physics model for simulating seismic response in layered fractured rocks |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US6753408P | 2008-02-28 | 2008-02-28 | |
US61/067,534 | 2008-02-28 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2009108432A1 true WO2009108432A1 (en) | 2009-09-03 |
Family
ID=41016426
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2009/032007 WO2009108432A1 (en) | 2008-02-28 | 2009-01-26 | Rock physics model for simulating seismic response in layered fractured rocks |
Country Status (5)
Country | Link |
---|---|
US (1) | US8184502B2 (en) |
EP (1) | EP2248007A1 (en) |
AU (1) | AU2009217648A1 (en) |
CA (1) | CA2710607A1 (en) |
WO (1) | WO2009108432A1 (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011086443A3 (en) * | 2010-01-12 | 2012-01-05 | Schlumberger Technology B.V. | Walkaway vsp calibrated sonic logs |
US20120092958A1 (en) * | 2010-10-18 | 2012-04-19 | Geobiz Inc. | Estimation of anisotropy from compressional waves from array sonic waveforms in well logging |
US8494827B2 (en) | 2009-09-25 | 2013-07-23 | Exxonmobil Upstream Research Company | Method of predicting natural fractures and damage in a subsurface region |
RU2492447C1 (en) * | 2012-03-15 | 2013-09-10 | Федеральное государственное бюджетное учреждение науки Институт проблем нефти и газа РАН | Method to define anisotropy of pore space and position of main axes of permeability tensor of rocks on core |
US8725478B2 (en) | 2010-08-09 | 2014-05-13 | Conocophillips Company | Reservoir upscaling method with preserved transmissibility |
CN103827696A (en) * | 2011-07-27 | 2014-05-28 | 普拉德研究及开发股份有限公司 | Multi-well anisotropy inversion |
US8908473B2 (en) | 2008-12-23 | 2014-12-09 | Schlumberger Technology Corporation | Method of subsurface imaging using microseismic data |
EP2847429A4 (en) * | 2012-06-14 | 2016-01-27 | Halliburton Energy Services Inc | System, method,&computer program product to determine placement of fracture stimulation points using mineralogy |
CN106154320A (en) * | 2015-04-22 | 2016-11-23 | 中国石油化工股份有限公司 | The detection method of fracture fluid |
US9696442B2 (en) | 2013-10-29 | 2017-07-04 | Exxonmobil Upstream Research Company | Method for estimating subsurface properties from geophysical survey data using physics-based inversion |
EP2513674A4 (en) * | 2009-12-18 | 2017-11-08 | Chevron U.S.A., Inc. | Workflow for petrophysical and geophysical formation evaluation of wireline and lwd log data |
US10619481B2 (en) | 2015-07-22 | 2020-04-14 | Conocophillips Company | Micromechanical elastic properties solver |
CN111255435A (en) * | 2020-01-17 | 2020-06-09 | 西安石油大学 | Method for calculating shale content of complex reservoir |
US20220342113A1 (en) * | 2021-04-21 | 2022-10-27 | Baker Hughes Oilfield Operations Llc | Estimation of properties of a subterranean region using a synthetic physical model |
Families Citing this family (68)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8243549B2 (en) * | 2008-03-10 | 2012-08-14 | Schlumberger Technology Corporation | Estimating seismic anisotropy of shales |
WO2010056427A1 (en) * | 2008-11-14 | 2010-05-20 | Exxonmobil Upstream Research Company | Forming a model of a subsurface region |
US20100135115A1 (en) * | 2008-12-03 | 2010-06-03 | Chevron U.S.A. Inc. | Multiple anisotropic parameter inversion for a tti earth model |
US20110141851A1 (en) * | 2009-12-16 | 2011-06-16 | Chevron U.S.A. Inc. | System and method for integrated reservoir and seal quality prediction |
US8515720B2 (en) * | 2010-04-06 | 2013-08-20 | Schlumberger Technology Corporation | Determine field fractures using geomechanical forward modeling |
US8547788B2 (en) * | 2010-05-17 | 2013-10-01 | Schlumberger Technology Corporation | Methods for making acoustic anisotropy logging while drilling measurements |
US8792302B2 (en) * | 2010-05-26 | 2014-07-29 | Westerngeco L.L.C. | Estimating anisotropic parameters |
US8457940B2 (en) | 2010-07-29 | 2013-06-04 | Schlumberger Technology Corporation | Model-consistent structural restoration for geomechanical and petroleum systems modeling |
US9213122B2 (en) * | 2010-09-24 | 2015-12-15 | Schlumberger Technology Corporation | Single well anisotropy inversion using velocity measurements |
US9182913B2 (en) * | 2011-10-18 | 2015-11-10 | Ubiterra Corporation | Apparatus, system and method for the efficient storage and retrieval of 3-dimensionally organized data in cloud-based computing architectures |
US20130124162A1 (en) * | 2011-11-16 | 2013-05-16 | Conocophillips Company | Method of calculating a shape factor of a dual media fractured reservoir model from intensities and orientations of fracture sets for enhancing the recovery of hydrocarbins |
CA2866597C (en) * | 2012-03-06 | 2018-05-22 | Ion Geophysical Corporation | Model predicting fracturing of shale |
WO2013148623A1 (en) * | 2012-03-26 | 2013-10-03 | Schlumberger Canada Limited | Dynamic to static elastic property transforms for tiv rock material |
EP2856373A4 (en) * | 2012-05-24 | 2016-06-29 | Exxonmobil Upstream Res Co | System and method for predicting rock strength |
US9465140B2 (en) * | 2012-06-22 | 2016-10-11 | Exxonmobil Upstream Research Company | Petrophysical method for predicting shear strength anisotropy in fine-grained rock formations |
WO2014008931A1 (en) * | 2012-07-10 | 2014-01-16 | Statoil Petroleum As | Anisotropy parameter estimation |
CN102788994B (en) * | 2012-07-12 | 2015-01-21 | 恒泰艾普石油天然气技术服务股份有限公司 | Reservoir fracture determining method |
EP2859387A1 (en) | 2012-07-13 | 2015-04-15 | Halliburton Energy Services, Inc. | Method of estimating anisotropic formation resistivity profile using a multi-component induction tool |
EP2915064B1 (en) | 2012-11-05 | 2022-10-19 | Fnv Ip B.V. | Method of indicating the presence of gas hydrate and shallow gas in deepwater environment |
CA2920008C (en) * | 2013-07-29 | 2017-06-13 | Cgg Services Sa | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media |
AU2013402201B2 (en) | 2013-10-01 | 2017-07-13 | Landmark Graphics Corporation | In-situ wellbore, core and cuttings information system |
US9664816B2 (en) | 2013-12-06 | 2017-05-30 | Halliburton Energy Services, Inc. | Fracture detection and characterization using resistivity images |
JP6198242B2 (en) * | 2014-04-07 | 2017-09-20 | 株式会社神戸製鋼所 | Boundary detection method |
MX2016014193A (en) * | 2014-04-30 | 2017-05-03 | Halliburton Energy Services Inc | Characterizing a downhole environment using stiffness coefficients. |
US10221659B2 (en) * | 2014-10-08 | 2019-03-05 | Chevron U.S.A. Inc. | Automated well placement for reservoir evaluation |
US10400550B2 (en) | 2014-10-24 | 2019-09-03 | Halliburton Energy Services, Inc. | Shale fracturing characterization and optimization using three-dimensional fracture modeling and neural network |
WO2016083893A1 (en) * | 2014-11-25 | 2016-06-02 | Cgg Services Sa | Modeling an elastic stiffness tensor in a transverse isotropic subsurface medium |
CN105759310B (en) * | 2014-12-16 | 2018-05-08 | 中国石油化工股份有限公司 | The attenuation of seismic wave and velocity dispersion Forecasting Methodology in complicated heterogeneous reservoir medium |
WO2016099488A1 (en) * | 2014-12-17 | 2016-06-23 | Halliburton Energy Services, Inc. | Geomechanical model of stresses on an orthorhombic media |
CN104484573B (en) * | 2014-12-30 | 2017-09-15 | 中国石油天然气股份有限公司 | A kind of method of determination stratum stiffness coefficient |
WO2016118360A1 (en) * | 2015-01-23 | 2016-07-28 | Halliburton Energy Services, Inc. | A combination model for predicting stiffness coefficients absent stoneley wave velocity data |
US10724346B2 (en) * | 2015-03-31 | 2020-07-28 | Halliburton Energy Services, Inc. | Simplified geomechanical model of stresses on an orthorhombic media |
US10191167B2 (en) * | 2015-06-30 | 2019-01-29 | Halliburton Energy Services, Inc. | Correcting the effects of deviation and dispersion on sonic log measurements of deviated wells in laminated formations |
EP3326012B1 (en) * | 2015-07-22 | 2020-10-28 | ConocoPhillips Company | Micromechanical elastic properties solver |
MX2018002023A (en) * | 2015-09-18 | 2018-04-13 | Halliburton Energy Services Inc | Global inversion based estimation of anisotropy parameters for orthorhombic media. |
CN105425287B (en) * | 2015-11-09 | 2017-10-27 | 中国地质大学(北京) | The prestack separation method of seismic wave |
CN105629303B (en) * | 2015-12-28 | 2017-12-08 | 中国石油大学(北京) | Based on petrophysical pre-stack fracture quantitative forecasting technique and system |
EP3203270B1 (en) | 2016-02-03 | 2024-04-10 | CGG Services SAS | Complex pore geometry modeling by continuously varying inclusions (ci) method for elastic parameter prediction using inclusion models |
WO2017184123A1 (en) * | 2016-04-19 | 2017-10-26 | Halliburton Energy Services, Inc. | Fracture and stress characterization using layer thickness variation |
CN105866003B (en) * | 2016-05-16 | 2019-05-17 | 中国地质大学(武汉) | A kind of shale gas reservoir total porosity acquisition methods and device |
US10274625B2 (en) * | 2016-11-09 | 2019-04-30 | Chevron U.S.A. Inc. | System and method for porosity estimation in low-porosity subsurface reservoirs |
GB2562285B (en) * | 2017-05-11 | 2022-03-02 | Equinor Energy As | Fluid substitution |
US10990882B2 (en) | 2017-07-28 | 2021-04-27 | International Business Machines Corporation | Stratigraphic layer identification from seismic and well data with stratigraphic knowledge base |
US10991078B2 (en) | 2017-09-15 | 2021-04-27 | Saudi Arabian Oil Company | Inferring petrophysical properties of hydrocarbon reservoirs using a neural network |
US10983237B2 (en) | 2018-04-13 | 2021-04-20 | Saudi Arabian Oil Company | Enhancing seismic images |
US10891462B2 (en) | 2018-06-29 | 2021-01-12 | Saudi Arabian Oil Company | Identifying geometrical properties of rock structure through digital imaging |
CN110187008A (en) * | 2019-04-29 | 2019-08-30 | 中国科学院地理科学与资源研究所 | A kind of full elastic parameter detection device of lossless sound wave and method |
CN111983670B (en) * | 2019-05-21 | 2022-10-14 | 中国石油化工股份有限公司 | Rock physical characteristic distribution rule analysis method and computer system |
CN110107282B (en) * | 2019-05-30 | 2022-05-20 | 中国石油天然气股份有限公司 | Method for controlling while-drilling track of thin-layer shale oil horizontal well |
CN110362857B (en) * | 2019-06-04 | 2023-05-16 | 绍兴文理学院 | Theoretical model analysis method for anisotropic distribution characteristics of rock structural surface |
CN112213770B (en) * | 2019-07-09 | 2023-11-28 | 中国石油天然气股份有限公司 | Method and device for identifying hydrocarbon-containing sandstone in reservoir based on logarithmic domain difference distribution characteristics |
CN112558154B (en) * | 2019-09-25 | 2023-10-13 | 中国石油化工股份有限公司 | Orthogonal anisotropic medium speed modeling method and system |
US11372123B2 (en) | 2019-10-07 | 2022-06-28 | Exxonmobil Upstream Research Company | Method for determining convergence in full wavefield inversion of 4D seismic data |
CN112824940A (en) * | 2019-11-20 | 2021-05-21 | 中国石油天然气集团有限公司 | Amplitude matching method and device for reflected longitudinal wave and reflected converted transverse wave |
US11493658B2 (en) * | 2019-12-03 | 2022-11-08 | China Petroleum & Chemical Corporation | Computer-implemented method and system employing nonlinear direct prestack seismic inversion for poisson impedance |
CN111188612B (en) * | 2020-01-13 | 2022-12-13 | 中国石油天然气股份有限公司大港油田分公司 | Method for quickly identifying shale oil dessert with well logging multi-parameter fusion |
US20230098683A1 (en) * | 2020-03-13 | 2023-03-30 | Reveal Energy Services, Inc. | Determining a dimension associated with a wellbore |
CN111427082A (en) * | 2020-04-08 | 2020-07-17 | 中国石油天然气股份有限公司 | Method and device for determining seismic response template and storage medium |
CN113866824A (en) * | 2020-06-30 | 2021-12-31 | 中国石油化工股份有限公司 | Method for determining elastic parameters of stratum between salts, storage medium and computer equipment |
US11624277B2 (en) | 2020-07-20 | 2023-04-11 | Reveal Energy Services, Inc. | Determining fracture driven interactions between wellbores |
CN111911146B (en) * | 2020-08-17 | 2023-04-25 | 中国石油大学(华东) | Reservoir fracturing effect evaluation method based on fracture intersection plate |
CN112305589A (en) * | 2020-09-22 | 2021-02-02 | 中国石油天然气集团有限公司 | Method and device for imaging depth domain of anisotropic medium |
US11767751B2 (en) | 2020-11-04 | 2023-09-26 | Reveal Energy Services, Inc. | Determining a characteristic associated with a reservoir |
US11668847B2 (en) | 2021-01-04 | 2023-06-06 | Saudi Arabian Oil Company | Generating synthetic geological formation images based on rock fragment images |
CN113740896B (en) * | 2021-07-30 | 2022-06-03 | 中国石油大学(华东) | Method for predicting orthogonal anisotropic reservoir fractures and storage medium |
CN114002744B (en) * | 2021-10-29 | 2023-07-14 | 重庆科技学院 | Dense sandstone reservoir fluid identification method based on deep learning |
CN114280272A (en) * | 2021-12-13 | 2022-04-05 | 华能煤炭技术研究有限公司 | Coal rock physical parameter analysis method, system and electronic equipment |
CN116559959A (en) * | 2023-05-24 | 2023-08-08 | 中国石油大学(华东) | Deep fracture type reservoir reflection coefficient equation establishment method and related equipment |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020120429A1 (en) * | 2000-12-08 | 2002-08-29 | Peter Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
US6901333B2 (en) * | 2003-10-27 | 2005-05-31 | Fugro N.V. | Method and device for the generation and application of anisotropic elastic parameters |
US20060074561A1 (en) * | 2004-09-16 | 2006-04-06 | Schlumberger Technology Corporation | Methods for visualizing distances between wellbore and formation boundaries |
US20060153005A1 (en) * | 2005-01-07 | 2006-07-13 | Herwanger Jorg V | Determination of anisotropic physical characteristics in and around reservoirs |
US20070097787A1 (en) * | 2005-10-28 | 2007-05-03 | Geomechanics International, Inc. | Hydrocarbon saturation determination using acoustic velocities obtained through casing |
US20080020229A1 (en) * | 2004-09-08 | 2008-01-24 | University Of Washington | Energy absorbent material |
Family Cites Families (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4375090A (en) * | 1977-12-23 | 1983-02-22 | Chevron Research Company | Method for interpreting seismic records to yield indications of gas/oil in an earth formation such as a sandstone, limestone, or dolostone |
FR2733073B1 (en) * | 1995-04-12 | 1997-06-06 | Inst Francais Du Petrole | METHOD FOR MODELING A LAMINATED AND FRACTURED GEOLOGICAL ENVIRONMENT |
US5828981A (en) * | 1995-05-11 | 1998-10-27 | Texaco Inc. | Generating pore types and synthetic capillary pressure curves from wireline logs using neural networks |
US5675147A (en) * | 1996-01-22 | 1997-10-07 | Schlumberger Technology Corporation | System and method of petrophysical formation evaluation in heterogeneous formations |
US5869755A (en) * | 1997-12-31 | 1999-02-09 | Schlumberger Technology Corporation | Porosity estimation method in carbonate rock |
US6088656A (en) * | 1998-11-10 | 2000-07-11 | Schlumberger Technology Corporation | Method for interpreting carbonate reservoirs |
US6529833B2 (en) * | 1998-12-30 | 2003-03-04 | Baker Hughes Incorporated | Reservoir monitoring in a laminated reservoir using 4-D time lapse data and multicomponent induction data |
CA2383711C (en) * | 2000-06-29 | 2014-03-18 | Stephen R. Kennon | Method and system for solving finite element models using multi-phase physics |
GB0018480D0 (en) * | 2000-07-27 | 2000-09-13 | Geco Prakla Uk Ltd | A method of processing surface seismic data |
FR2823877B1 (en) * | 2001-04-19 | 2004-12-24 | Inst Francais Du Petrole | METHOD FOR CONSTRAINING BY DYNAMIC PRODUCTION DATA A FINE MODEL REPRESENTATIVE OF THE DISTRIBUTION IN THE DEPOSIT OF A CHARACTERISTIC PHYSICAL SIZE OF THE BASEMENT STRUCTURE |
US6795773B2 (en) * | 2001-09-07 | 2004-09-21 | Halliburton Energy Services, Inc. | Well completion method, including integrated approach for fracture optimization |
US6925031B2 (en) * | 2001-12-13 | 2005-08-02 | Baker Hughes Incorporated | Method of using electrical and acoustic anisotropy measurements for fracture identification |
WO2003067025A2 (en) * | 2002-02-01 | 2003-08-14 | Regents Of The University Of Minnesota | Interpretation and design of hydraulic fracturing treatments |
US6714480B2 (en) * | 2002-03-06 | 2004-03-30 | Schlumberger Technology Corporation | Determination of anisotropic moduli of earth formations |
US7035165B2 (en) * | 2003-01-29 | 2006-04-25 | Baker Hughes Incorporated | Imaging near-borehole structure using directional acoustic-wave measurement |
US6904365B2 (en) * | 2003-03-06 | 2005-06-07 | Schlumberger Technology Corporation | Methods and systems for determining formation properties and in-situ stresses |
WO2004109334A2 (en) * | 2003-05-30 | 2004-12-16 | Halliburton Energy Services, Inc. | Systems and methods for analyzing carbonate formations while drilling |
US6970397B2 (en) * | 2003-07-09 | 2005-11-29 | Gas Technology Institute | Determination of fluid properties of earth formations using stochastic inversion |
US7042802B2 (en) * | 2003-09-18 | 2006-05-09 | Schlumberger Technology Corporation | Determination of stress characteristics of earth formations |
US7286939B2 (en) * | 2003-10-28 | 2007-10-23 | Westerngeco, L.L.C. | Method for estimating porosity and saturation in a subsurface reservoir |
US6959246B2 (en) * | 2003-12-29 | 2005-10-25 | Schlumberger Technology Corporation | Carbonate permeability |
CA2580570C (en) | 2004-12-06 | 2014-03-11 | Exxonmobil Upstream Research Company | Integrated anisotropic rock physics model |
US20060219402A1 (en) * | 2005-02-16 | 2006-10-05 | Commonwealth Scientific And Industrial Research Organisation | Hydraulic fracturing |
-
2009
- 2009-01-26 WO PCT/US2009/032007 patent/WO2009108432A1/en active Application Filing
- 2009-01-26 EP EP09714967A patent/EP2248007A1/en not_active Withdrawn
- 2009-01-26 CA CA2710607A patent/CA2710607A1/en not_active Abandoned
- 2009-01-26 AU AU2009217648A patent/AU2009217648A1/en not_active Abandoned
- 2009-01-26 US US12/863,867 patent/US8184502B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020120429A1 (en) * | 2000-12-08 | 2002-08-29 | Peter Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
US6901333B2 (en) * | 2003-10-27 | 2005-05-31 | Fugro N.V. | Method and device for the generation and application of anisotropic elastic parameters |
US20080020229A1 (en) * | 2004-09-08 | 2008-01-24 | University Of Washington | Energy absorbent material |
US20060074561A1 (en) * | 2004-09-16 | 2006-04-06 | Schlumberger Technology Corporation | Methods for visualizing distances between wellbore and formation boundaries |
US20060153005A1 (en) * | 2005-01-07 | 2006-07-13 | Herwanger Jorg V | Determination of anisotropic physical characteristics in and around reservoirs |
US20070097787A1 (en) * | 2005-10-28 | 2007-05-03 | Geomechanics International, Inc. | Hydrocarbon saturation determination using acoustic velocities obtained through casing |
Non-Patent Citations (2)
Title |
---|
"Simula+ documentation", FREE SOFTWARE FOUNDATION, INC, 1991, BOSTON (MA), XP008169022, Retrieved from the Internet <URL:http://www.lpmm.fr/collard/pub/simula+/doc/html/classmateriol_1_t> [retrieved on 20090315] * |
BAYUK I.O. ET AL.: "Upscaling of elastic properties of anisotropic sedimentary rocks", GEOPHYSICAL JOURNAL INTERNATIONAL, vol. 172, no. ISS.2, 21 November 2007 (2007-11-21), pages 842 - 860, XP008140451, Retrieved from the Internet <URL:http://www3.interscience.wiley.com/joumal/119410236/abstract?CRETRY=18SRETRY=0> [retrieved on 20090312] * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8908473B2 (en) | 2008-12-23 | 2014-12-09 | Schlumberger Technology Corporation | Method of subsurface imaging using microseismic data |
US8494827B2 (en) | 2009-09-25 | 2013-07-23 | Exxonmobil Upstream Research Company | Method of predicting natural fractures and damage in a subsurface region |
EP2513674A4 (en) * | 2009-12-18 | 2017-11-08 | Chevron U.S.A., Inc. | Workflow for petrophysical and geophysical formation evaluation of wireline and lwd log data |
GB2489156B (en) * | 2010-01-12 | 2015-12-02 | Schlumberger Holdings | Method and system for processing subsurface data |
US8411529B2 (en) | 2010-01-12 | 2013-04-02 | Schlumberger Technology Corporation | Walkaway VSP calibrated sonic logs |
GB2489156A (en) * | 2010-01-12 | 2012-09-19 | Schlumberger Holdings | Walkaway vsp calibrated sonic logs |
WO2011086443A3 (en) * | 2010-01-12 | 2012-01-05 | Schlumberger Technology B.V. | Walkaway vsp calibrated sonic logs |
NO344426B1 (en) * | 2010-01-12 | 2019-12-02 | Schlumberger Technology Bv | Method and system for combining sonic and seismic log data from boreholes |
US8725478B2 (en) | 2010-08-09 | 2014-05-13 | Conocophillips Company | Reservoir upscaling method with preserved transmissibility |
US20120092958A1 (en) * | 2010-10-18 | 2012-04-19 | Geobiz Inc. | Estimation of anisotropy from compressional waves from array sonic waveforms in well logging |
CN103827696A (en) * | 2011-07-27 | 2014-05-28 | 普拉德研究及开发股份有限公司 | Multi-well anisotropy inversion |
US10436921B2 (en) | 2011-07-27 | 2019-10-08 | Schlumberger Technology Corporation | Multi-well anisotropy inversion |
RU2492447C1 (en) * | 2012-03-15 | 2013-09-10 | Федеральное государственное бюджетное учреждение науки Институт проблем нефти и газа РАН | Method to define anisotropy of pore space and position of main axes of permeability tensor of rocks on core |
EP2847429A4 (en) * | 2012-06-14 | 2016-01-27 | Halliburton Energy Services Inc | System, method,&computer program product to determine placement of fracture stimulation points using mineralogy |
US9835027B2 (en) | 2012-06-14 | 2017-12-05 | Halliburton Energy Services, Inc. | System, method, and computer program product to determine placement of fracture stimulation points using minerology |
US10598009B2 (en) | 2012-06-14 | 2020-03-24 | Halliburton Energy Services, Inc. | System, method, and computer program product to determine placement of fracture stimulation points using minerology |
US9696442B2 (en) | 2013-10-29 | 2017-07-04 | Exxonmobil Upstream Research Company | Method for estimating subsurface properties from geophysical survey data using physics-based inversion |
CN106154320B (en) * | 2015-04-22 | 2018-07-20 | 中国石油化工股份有限公司 | The detection method of fracture fluid |
CN106154320A (en) * | 2015-04-22 | 2016-11-23 | 中国石油化工股份有限公司 | The detection method of fracture fluid |
US10619481B2 (en) | 2015-07-22 | 2020-04-14 | Conocophillips Company | Micromechanical elastic properties solver |
CN111255435A (en) * | 2020-01-17 | 2020-06-09 | 西安石油大学 | Method for calculating shale content of complex reservoir |
US20220342113A1 (en) * | 2021-04-21 | 2022-10-27 | Baker Hughes Oilfield Operations Llc | Estimation of properties of a subterranean region using a synthetic physical model |
Also Published As
Publication number | Publication date |
---|---|
US8184502B2 (en) | 2012-05-22 |
CA2710607A1 (en) | 2009-09-03 |
US20100312534A1 (en) | 2010-12-09 |
EP2248007A1 (en) | 2010-11-10 |
AU2009217648A1 (en) | 2009-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8184502B2 (en) | Rock physics model for simulating seismic response in layered fractured rocks | |
US8473213B2 (en) | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model | |
EP1820137B1 (en) | Integrated anisotropic rock physics model | |
Kuster et al. | Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations | |
Pride et al. | Permeability dependence of seismic amplitudes | |
Wang | Seismic anisotropy in sedimentary rocks, part 2: Laboratory data | |
Müller et al. | Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks—A review | |
Dou et al. | Rock-physics-based carbonate pore type characterization and reservoir permeability heterogeneity evaluation, Upper San Andres reservoir, Permian Basin, west Texas | |
Han et al. | Effects of porosity and clay content on wave velocities in sandstones | |
Quan et al. | Seismic attenuation tomography using the frequency shift method | |
Vanorio et al. | The effect of chemical and physical processes on the acoustic properties of carbonate rocks | |
Prioul et al. | Forward modeling of fracture-induced sonic anisotropy using a combination of borehole image and sonic logs | |
Massaro et al. | Analyzing a suitable elastic geomechanical model for Vaca Muerta Formation | |
Teixeira et al. | Elastic properties of rock salt in the Santos Basin: Relations and spatial predictions | |
Cui et al. | Seismic forward modeling of fractures and fractured medium inversion | |
Helms et al. | Rock physics based velocity-porosity correlations developed for estimation of the elastic properties of the Bakken formations of the Williston Basin, North Dakota | |
Xu et al. | Carbonate rock physics: Analytical models and validations using computational approaches and lab/log measurements | |
Aziz et al. | Development a statistical relationship between compressional wave velocity and petrophysical properties from logs data for JERIBE Formation ASMARI Reservoir in FAUQI Oil Field | |
ZHAO et al. | Seismic rock physics characteristics and modeling of ultra-deep carbonate reservoirs | |
Saberi | A closer look at rock physics models and their assisted interpretation in seismic exploration | |
Collet et al. | Estimating azimuthal stress‐induced P‐wave anisotropy from S‐wave anisotropy using sonic log or vertical seismic profile data | |
Sun et al. | Quantifying the effects of pore structure and fluid saturation on acoustic wave velocity in carbonates | |
Malki et al. | Novel Approach for Elastic Anisotropy Assessment in Organic-Rich Shales: A Case Study from Early-Mature Region of the Bakken Formation | |
Gordin et al. | Relating Acoustic Anisotropy to Kerogen Content in Unconventional Formations-A Case Study in A Kerogen-Rich Unconventional Carbonate | |
Benson et al. | A modeling solution for predicting (a) dry rock bulk modulus, rigidity modulus and (b) seismic velocities and reflection coefficients in porous, fluid-filled rocks with applications to laboratory rock samples and well logs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09714967 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2710607 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009217648 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: 12863867 Country of ref document: US |
|
ENP | Entry into the national phase |
Ref document number: 2009217648 Country of ref document: AU Date of ref document: 20090126 Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009714967 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |