You must log in to edit PetroWiki. Help with editing
Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information
Water saturation determination
DO NOT APPROVE THIS CHANGE: Changes are not visible unless the change has been approved by a moderator or administrator. These comments are not visible to the general public; they are only visible to the moderators and administrators. And also as long as you don't mark the change as a "minor edit", because minor edits bypass the moderators and go straight to me (Denise).
Water saturation (S_{w}) determination is the most challenging of petrophysical calculations and is used to quantify its more important complement, the hydrocarbon saturation (1 – S_{w}). Complexities arise because there are a number of independent approaches that can be used to calculate S_{w}. The complication is that often, if not typically, these different approaches lead to somewhat different S_{w} values that may equate to considerable differences in the original oil in place (OOIP) or original gas in place (OGIP) volumes. The challenge to the technical team is to resolve and to understand the differences among the S_{w} values obtained using the different procedures, and to arrive at the best calculation of S_{w} and its distribution throughout the reservoir vertically and areally. In OOIP and OGIP calculations, it is important to remember the relative importance of porosity and S_{w}. A 10% pore volume (PV) change in S_{w} has the same impact as a 2% bulk volume (BV) change in porosity (in a 20% BV porosity reservoir).
Techniques for calculating water saturation
S_{w} in wellbores can be determined by the following primary methods:
- S_{w} calculations from resistivity well logs by application of a model relating S_{w} to porosity, connate-water resistivity, and various rock electrical properties.
- S_{w} calculations from laboratory capillary pressure/saturation (P_{c}/S_{w}) measurements by application of a model relating S_{w} to various rock and fluid properties and height above the free-water level.
- S_{w} calculations using oil-based mud (OBM)-core-plug Dean-Stark water-volume determinations.
- Combinations of these methods.
This listing is the chronological order in which data are likely to become available, not in a ranked order based on the accuracy of the various methods. The choice of which S_{w}-calculation approach to use is often controlled by the availability of the various types of data. If no OBM cores have been cut, then this technique cannot be used unless funds are spent to acquire such data from one or more newly drilled wells. This is not a high incremental cost when OBM use is planned for other purposes. Resistivity logs are run in all wells, so these data are available for making standard-log-analysis S_{w} calculations. A key consideration when making calibrated S_{w} calculations is the availability of special-core-analysis (SCAL) data on core samples from the particular reservoir; that is, the number of laboratory electrical-property and P_{c}/S_{w} core-plug measurements that have been made.
The technique chosen to calculate S_{w} is often a hybrid that combines the use of two of these basic data sources. For example, the OBM-core S_{w} data can be used in combination with the resistivity logs to expand the data set used to include all wells and the whole of the hydrocarbon column. Alternatively, the OBM-core S_{w} data can be used in combination with the P_{c}/S_{w} data. In this way, the OBM-core S_{w} data define the S w values for the majority of the reservoir, whereas the P_{c}/S_{w} data define the S_{w} values in the interval just above the fluid contact and perhaps in areas of the field where P_{c} data are available but OBM-core data are not.
Data availability and data quality
This section discussed the input-data availability and data quality issues for each S_{w} technique. These considerations often control the initial choice of methodology to calculate S_{w} and need to be addressed at the start of the project to determine whether it is practically possible to fill gaps in the database in order to use a more accurate S_{w}-calculation approach. This discussion assumes that accurate porosity values are available from the routine-core-analysis database and that porosity is calculated point by point from the well logs. The discussion focuses on particular aspects that affect the choice of S_{w} methodology. Many of the database considerations are discussed on the petrophyscial database page.
Resistivity logs
Wells generally have one variety or another of laterolog or induction resistivity log because they are broadly useful and because government regulations typically demand that they be recorded. This generally provides point-by-point data from the top of the hydrocarbon column down through any aquifer intervals that are present. However, in many fields, the early wells are spread thinly over the reservoir area, but the later development wells are drilled only in areas chosen to maximize rate and recovery while minimizing costs. This means that, often, few wells are drilled downdip where the hydrocarbon column thins because of an underlying aquifer, or in the potentially thin updip limits of the reservoir. In such areas, there may, therefore, be few resistivity logs.
Laterologs are preferred to induction logs when the drilling mud has moderate to high salinity. This limitation of induction tools arises because of the excessive conductivity signal from the borehole and the mud-filtrate-invaded zone. Deep laterolog tools read too high when measuring immediately beneath anhydrite and salt, ^{[1]} and alternative resistivity curves should be selected. When formation resistivity, R_{t}, is very high, previous generations of induction tools had limited accuracy, but current tools are much improved. Although the deep induction measurement is a running average over many vertical feet, modern tools include systems to deconvolve the raw log and provide a final log with a good vertical resolution.
Deep invasion of water-based mud (WBM) filtrate affects all resistivity logs, and, in the extreme, the available resistivity log may be used only qualitatively. At the opposite extreme, when oil-based mud (OBM) filtrate invades a hydrocarbon reservoir, the invading OBM filtrate generally displaces only the reservoir oil and gas, leaving the S_{w} unchanged. Here, invasion of OBM does not usually change the deep-formation or the invaded-zone resistivity. For moderate invasion depths, the logging company charts are sometimes used to correct the deep-reading log to provide a better estimate of R_{t}.
Pad-mounted shallow-reading microresistivity logs measure R_{xo}, the resistivity of the mud-filtrate-invaded zone. When used together with the deeper-reading tools, these logs provide valuable information about the mobility of the reservoir fluids, including the presence of tar. In WBM wells, they also provide an estimate of the residual-hydrocarbon saturation, S_{orw}.
Connate-brine resistivity data
An accurate value of connate-brine resistivity, R_{w}, or its values and distribution throughout the reservoir, are required for accurate S_{w} calculations using resistivity logs. Temperature estimates are also required.
A first check on the R_{w} of the aquifer is to back-calculate the apparent R_{w} with the Archie equation using the invasion-corrected resistivity logs and the best estimates of a and m parameters. Because S_{w} is typically 100% PV in the aquifer interval, the n value is not relevant here.
The spontaneous-potential (SP) log provides a second method to calculate R_{w} in wells drilled with WBM. Information on the mud-filtrate composition and temperature is used with the SP deflection to calculate R_{w}.^{[1]}^{[2]} The moderately accurate calculation process is valid in the aquifer but is also valid in the hydrocarbon column if high resistivity does not suppress the SP response. When OBM-core salinity measurements are not available, the SP log provides the only evidence of possible R_{w} variations in the hydrocarbon column.
A third estimate of aquifer-water composition and R_{w} is often taken from samples recovered during flow tests of the aquifer interval; however, the R_{w} of the oil and/or gas column is not always the same as that of the aquifer interval. ^{[3]}^{[4]}^{[5]} Aquifer-interval flow tests must be validated and checked for contamination from mud-filtrate invasion.
For the oil or gas column, the determination of the R_{w} value or values is far more of a challenge because the reservoir water will not flow. The typical, but not necessarily correct, first assumption is that the hydrocarbon-column R_{w} is the same as that of the underlying aquifer. If wells have been cored with OBM, core plugs from the hydrocarbon- and water-bearing intervals can be analyzed for both their water volume and their salt content, particularly the chloride ion that in almost all cases dominates the anion side of the salinity determination. ^{[3]}
Fig. 1 shows variation of chloride concentration with depth for a reservoir in Ecuador. ^{[4]} The chloride value can generally be used to quantify the reservoir-water salinity, from which the R_{w} at reservoir conditions can be calculated using standard water-resistivity vs. chloride charts or algorithms. For reservoirs in which there is a considerable CO_{2} content (3+ mol%), the ion distribution at surface conditions will differ from that at reservoir temperature and pressure. Equilibrium ion-distribution calculations need to be made when adjusting the surface-salinity measurement to reservoir conditions.
Fig. 1 – Formation-water salinity variation within the hydrocarbon column (Villano oil field, Ecuador).^{[4]} Hollin is the name of the oil-reservoir formation, and WOC is the water/oil contact. Chloride concentration in the reservoir brine varies from approximately 2,000 to 35,000 ppm. Lowest invasion refers to limited OBM-filtrate invasion into the cores.
Formation temperature affects the S_{w} estimates because, for constant formation-water composition, R_{w} varies with temperature. ^{[1]} Maximum downhole temperature is measured with most log runs and drillstem-tool (DST) tests, and these are widely used to estimate a temperature vs. depth profile. It can be argued that the temperature required for resistivity-derived S_{w} estimates is the prevailing temperature in the rock volume seen by the tool at the time of logging. At this time, the relevant rock is likely to be cooler than the original formation temperature. The error induced by the usual maximum temperature simplification is not large, and the cooling issue is generally ignored.
Electrical-property SCAL data
The third aspect of making these S_{w} calculations is the choice of the model for the "electrical network" within the rock. These models relate S_{w} to several formation variables including the bulk-formation resistivity and the formation-water resistivity. A number of models have been published, for example^{[6]}^{[7]} :
- Archie
- Waxman-Smits-Thomas (WST)
- Dual-water (DW)
- Indonesia
Laboratory measurements of two or more types of electrical properties are taken. All of these models assume a homogeneous rock sample.
Archie exponents
First, a set of cleaned core plugs with a range of porosities are fully saturated with brine of known resistivity, and the bulk resistivity of each core plug is measured. For this simplest model, the slope of a line fitted to a log-log plot of the data set gives the cementation exponent, m, and the intercept is the cementation constant, ^{[8]} a (see Fig. 2, where a = 1 and m = 1.77). These parameters are used to predict point-by-point F from porosity; leading to predictions of R_{0} and S_{w}.
where F = formation factor, R_{w} = brine-water resistivity, and R_{0} = rock resistivity with zero oil and gas saturation (100% PV S_{w}). The plotted logarithmic data (log_{10}F and log_{10}ϕ) are fitted with a linear model of the form,
where ϕ = porosity, a = cementation constant, and m = cementation exponent. Therefore, m = − change in log 10 F /change in log 10 ϕ (the slope of the line-fit) and a = F at 100% BV porosity (the line-fit intercept).
This model was developed by Archie, ^{[8]} who proposed a = 1.0 and m = 1.8 to 2.0 for his data set. Subsequent work by Exxon researchers for several sandstone rocks recommended a = 0.61 and m = 2.15 (the Humble formula). ^{[10]} Carbonates also have been studied and yielded a recommendation to use m = 1.87 + 0.019/ϕ below 9% BV (the Shell formula). ^{[1]} However, carbonate pore and fracture networks vary greatly, and m values from 1.0 to 3.0 may be required. Clearly, m is not a constant, but varies with rock type.
When plotting these formation-factor data, it is typically assumed that the rock samples have similar pore geometry, but with differing levels of porosity and diagenesis. Reservoir-specific exponent values are likely to provide more-accurate S_{w} results than worldwide correlations. However, before reservoir-specific values are determined, descriptive and experimental data need to be studied to determine whether they need to be subdivided into various groupings that relate to distinct differences in lithological properties like:
- Grain size
- Sorting
- Clay-mineral content
In partially brine-saturated rocks, a related experimental study involves measuring electrical properties as a function of water saturation. In these experiments, the resistivity index (I_{R}), the ratio of the desaturated-rock resistivity to the 100% PV brine-saturated rock resistivity (R_{t}/R_{0}), is measured as a function of brine saturation. For example, in a porous-plate apparatus, S_{w} is changed by increasing gas pressure, and therefore capillary pressure, at the gas/water interface in the pores. Brine flows from the base of the plug via a porous-plate. From the measurements on each core plug, a log-log plot of I_{R} vs. S_{w} is made (see Fig. 3, where n = 1.64). The slope of the line (almost always forced through I_{R} = 1.0 at S_{w} = 100% PV) is the Archie saturation exponent n (see Eqs 3 and 4). On the basis of experimental data, Archie^{[8]} recommended that n = 2.0, and this value is still widely used when no experimental data are available. Although cementation exponents can be determined from log analysis, saturation exponents cannot and, therefore, require external information from core data.
where n = saturation exponent, the slope from the origin of a line-fit of several data points; I_{R} = resistivity index; and S_{w} = fractional brine water saturation.
A straight line-fit is usually used, but curved line-fits can be considered where necessary. Curvature is often the result of the clay-mineral content but may also result from an inhomogeneous water distribution at the pore scale (e.g., when microporous rock grains are present). When significant amounts of clay minerals are present in the rocks, other models are required to extend the Archie relationships. The WST model, discussed next, is based on laboratory SCAL measurements including cation-exchange capacity (CEC).
Waxman-Smits-Thomas exponents and cation-exchange capacity
WST cementation and saturation exponents (m* and n*) are required to apply the WST shaly-sand-model equation discussed below. The quantity of cation-exchange sites per gram of rock sample (CEC) may be measured in the laboratory by several methods and, after converting to CEC per unit PV, is used as the model parameter Q_{V}.^{[11]}^{[12]} The most reliable measurement of Q_{V} involves carrying out bulk-rock resistivity, R_{0}, tests at several brine resistivities and, therefore, is time consuming. The rock conductivity values (1/R_{0}) are plotted vs. brine conductivity (1/R_{w}) to identify the excess conductivity resulting from the shales and clay minerals. The slope of the fitted line is the reciprocal of F*, the WST formation factor. The excess conductivity is modeled as being equal to BQ_{v}/F*, and B is presumed in this model to be always positive. The parameter B is the equivalent counter-ion conductance, ^{[11]}^{[12]} which is a function of temperature and the free-water resistivity. Q_{v} is estimated from the values of F* and B. Core resistivities are also measured when S_{w} is less than 100% PV and both the WST exponents m* and n* are derived (see Fig. 4 and Fig. 5, respectively). It should be noted that m* > m and n* > n, except in "clean" sands.
Other CEC methods require the breaking up, disaggregation, and consequent partial loss of the real geometry of the rock’s electrical network. These simpler methods, such as the ammonia method, use analytical-chemistry methods to measure CEC. After measuring porosity and grain density, this practical laboratory unit is converted to the required Q_{v} parameter. ^{[14]} These simpler CEC measurements are often made on sidewall cores and are used together with exponent values measured on cores from neighboring wells.
Numerous other shaly-sand models have been developed, and, unlike WST, many are calculated from effective porosity. These types of models are generally applied using Archie exponents. When using SCAL electrical-property data, there must be consistency between the electrical-network model used to derive the laboratory parameters and the model used in the final S_{w} calculations from the porosity and the resistivity logs (e.g., if the laboratory provides standard Archie n values, these are not appropriate input to the WST equation).
Capillary pressure SCAL data
P_{c} data are a different type of SCAL data that can be taken experimentally in several ways. All P_{c} saturation tests respond to the pore-size distribution of the rock and the interfacial properties of the various solid/fluid systems. These data are obtained by desaturating core plugs, either using a centrifuge or a porous-plate apparatus. Initially, cleaned and dry plugs are saturated with either water or oil. The liquid is then displaced by air or nitrogen. Because air is very nonwetting compared with either water or oil, using these fluid pairs (air/water or air/oil) means that, as the P_{c} increases, the air will first occupy the largest pores. As the P_{c} and air saturation increase, the air will occupy smaller and smaller pores. The core plug begins the experiment saturated with the wetting phase, so the desaturation process provides data for the drainage P_{c} curve. After completing the drainage process, the core plug can be spun under the liquid in a centrifuge experiment, the liquid saturation will increase, and the imbibition P_{c} curve will be generated. Usually, only drainage P_{c}/S_{w} data are taken, and for most reservoir situations, these are the relevant data because they correspond to the original oil (or gas) trap-filling process.
Mercury-injection capillary pressure (MICP) data are taken on cleaned and dried irregular core pieces. The core pieces are evacuated to a low vacuum, and mercury is injected with increasing pressure, up to 20,000 psi and sometimes higher. Corrections for clay-mineral adsorbed-water removed during drying can be made with the Hill-Shirley-Klein method. ^{[15]} The MICP experiment has the advantage of being run rapidly but is not a true wetting/nonwetting system. The sample cannot be used for subsequent SCAL tests because some mercury is retained within the core pieces at the end of the testing sequence. MICP data are widely used to measure pore-size distribution, but, when considering whether they should be used for accurate S_{w} calculations, MICP should be compared with air/water or air/oil P_{c}/S_{w} data.
The P_{c}/S_{w} data are usually compared first on a Leverett "J-function" basis. ^{[16]} The P_{c} data are converted to the J-function basis by multiplying each P_{c} value by the square root of its permeability divided by porosity and then dividing by the fluid-pair IFT multiplied by the contact angle (see Eq. 5). J-function values differ depending on whether they are calculated in oilfield or metric units. The J-function approach assumes similar pore-size distribution in all cores tested. In this way, the various P_{c}/S_{w} data tend to converge when the underlying assumptions are met; however, there may still be enough scatter to suggest that the data need to be divided into two or more groupings (see Fig. 6).
From a J-function vs. S_{w} plot, the technical team can determine whether enough data have been obtained, whether new data need to be gathered to fill in portions of the data ranges, and whether the data indicate that subgroupings are appropriate and needed. Also, this plot indicates whether there are significant outliers that should be excluded or examined in more detail. A drawback to this averaging method is the introduction into S_{w} determination of four measured parameters and their associated errors (i.e., porosity, permeability, IFT, and contact angle).
Fig. 6 – Example plots of centrifuge air/brine capillary pressure (P_{c}/S_{w}) data and its conversion to a J-function basis; data from an Asian gas field. The air/brine P_{c} value of 100 psi is equivalent to a height above free-water level of 200 to 350 ft, depending on the reservoir fluid’s properties and the temperature and pressure conditions.
Capillary pressure data may also be averaged by various models. ^{[17]}^{[18]} The relationship of S_{w} with permeability, and then porosity, is examined and is followed by examination of the height dependency.
OBM core water saturation data
The last type of S_{w} data discussed here is that obtained from routine core analysis of core plugs cut from OBM cores, either preserved as whole cores or else with core plugs cut at the wellsite and preserved individually. These data are taken foot-by-foot and are direct measurements of reservoir S_{w} values. ^{[4]}^{[5]}^{[19]}^{[20]}^{[21]}^{[22]} Many fields may never have had any wells cored with OBM; others may have only one or two OBM-cored wells. Even a single OBM coring of the full reservoir interval offers significant data that may impact the technical team’s methodology for making the S_{w} calculations. It is better to have at least two wells fully cored with OBM from different areas of the reservoir.
To evaluate the OBM-core S_{w} data, they should be plotted as S_{w} vs. log_{10} (permeability) or vs. porosity to identify outliers and trends in the data. Particularly, the low-porosity/low-permeability data range should be examined for potential measurement problems. Sometimes, the raw laboratory measurements of water volume and PV data need to be reviewed for problem points and recalculations made where appropriate. Finally, the data should be divided into various possible interval groupings so that any needed zonation can be identified.
If reservoir connate water has flowed out of the core plug at any stage before the laboratory measurement, the OBM-core S_{w} data are clearly not representative of the in-situ reservoir S_{w}. This certainly occurs in water-bearing formations and can also occur in the lowest intervals of the oil/water or gas/water transition zones. These lowest intervals, which may be a few feet to approximately 30-ft thick, are precisely the same intervals in which a water cut is expected with the initial oil production. The mobile-water intervals can be identified in OBM wells where the shallow-reading induction-log resistivity is higher than the deep-reading induction-log resistivity. This pattern indicates higher oil saturations in the invaded zone compared with the original oil saturations. Where mobile water is observed, the OBM-core S_{w} measurements do not represent in-situ S_{w} and are too low.
Application of each water saturation technique
Methodologies for quantifying S_{w} at the wellbore are discussed here. The main features of each approach are described; however, in some cases, there are variations that are not addressed. For each technique, its strengths and weaknesses are discussed.
Resistivity well logs and relational model
The most common technique for calculating S_{w} is the use of resistivity logs with a model (empirical or theoretical) that relates S_{w} to R_{t}, R_{w}, and porosity. As mentioned previously, a large number of R_{t}/S_{w} models have been published. The models are applied at every data point in the reservoir where deep resistivity, porosity, and shale-volume estimates, if required, are available. The evaluation of all other necessary parameters (constant or variable R_{w} values, a, m, n, Q_{V}, V_{sh}, R_{0} = F • R_{w}, etc.) has been discussed previously. Several commercial software packages are available that perform these S w calculations for a variety of log models.
Clean sand (Archie) model
and, alternatively,
This model^{[8]} is used for field studies in the many sandstone and carbonate reservoirs in which the clay-mineral content is low. This decision is strengthened after SCAL data have demonstrated that the simplest solution is satisfactory. When a significant fraction of smectite (montmorillonite) is present and where finely laminated sand and shale sequences occur, one of the shaly-sand models is very likely to be required. Low-resistivity pay is an issue in several oil-producing areas, such as the US Gulf Coast, Egypt, and Indonesia, and hydrocarbon reserves can be missed and left undiscovered as a result of the resistivity suppression by clay minerals and shales.
Shaly-sand model
In the clean-sand model, the formation water is the only electrically conductive medium. In shaly rocks, R_{t} is suppressed and Archie S_{w} calculations are too high. As clay-mineral-rich rocks were studied and experimentally tested, more-complicated electrical models were developed to account for the effects of the geometries of conductive clay minerals and shale on rock resistivity. The primary goal of the shaly-sand models is to determine a working relationship between S_{w} using parameters similar to the Archie model, but also incorporating the quantity and specific electrical properties of the clay-mineral/shale. All of the shaly-sand models reduce to the Archie equation when the shale component is zero. For simplicity, in all of the shaly-sand models, the cementation constant, a, is taken to be 1.0 but, if required, can be easily associated again with the R_{w} term.
Laminated sand/shale model
A parallel resistor model might be used for laminated sands, with multiple thin parallel layers of 100% shale interbedded with clean-sand layers. Thin, in this context, means that there are several beds within the vertical resolution of the resistivity-logging tool.
where the clean-sand resistivity . For this laminated shale/sand model, effective porosity depends simply on the sand fraction of the bulk volume:
The value of ϕ_{sd} may be assumed from neighboring thick sands, and all of the parameters, except the S_{w} of the sand, S_{wsd}, can be estimated.
Poupon-Leveaux (Indonesia) model
The Indonesia model was developed by field observation in Indonesia rather than by laboratory experimental measurement support. ^{[9]} It remains useful because it is based on readily available standard log-analysis parameters and gives reasonably reliable results. The formula was empirically modeled with field data in water-bearing shaly sands, but the detailed functionality for hydrocarbon-bearing sands is unsupported, except by common sense and long-standing use. S_{w} results from the formula are comparatively easy to calculate and, because it is not a quadratic equation, it gives results that are always greater than zero. Several of the other quadratic and iterative-solution models can calculate unreasonable negative S_{w} results.
The Indonesia model, ^{[9]} and other similar models, are often used when field-specific SCAL rock electrical-properties data are unavailable but are also sometimes used where the SCAL exponents do not measure the full range of shale volumes. Although it was initially modeled on the basis of Indonesian data, the Indonesia model can be applied everywhere. The inputs are the effective porosity, ϕ_{e}, shale volume and resistivity (V_{sh} and R_{sh}), and water and deep resistivities (R_{w} and R_{t}). The S_{w} output is usually taken to be the water saturation of the effective porosity, but it has been recently suggested that the output is likely to estimate S_{wt}.^{[7]} Many other log-based shaly-sand models have been proposed^{[23]} but, for brevity, are not discussed here.
Waxman-Smits-Thomas and dual-water models
S_{wt}, the water saturation of the total porosity, is calculated at each reservoir data point by iterative solution of the complex multiparameter Waxman-Smits-Thomas (WST) and dual-water (DW) equations (Eqs. 7 and 8). For brevity, the details^{[1]}^{[11]}^{[12]}^{[13]}^{[14]}^{[24]}^{[25]} of the solution methods are not presented here. The WST and DW models are total-porosity/S_{w} system models.
The WST model is based on laboratory measurements of resistivity, porosity, and saturation of real rocks. ^{[11]}^{[12]}^{[13]}^{[14]} Q_{v} is the cation-exchange capacity (CEC) per unit PV.
where S_{wt} = water saturation of the total porosity as shown schematically in Fig. 1, B = specific cation conductance in (1/ohm•m)/(meq/mL), and Q_{V} = CEC in meq/mL of total PV. The exponents m* and n* apply to the total PV.
The DW model^{[7]}^{[24]}^{[25]} is also based on the WST data. It uses clay-bound-water conductivity instead of WST’s BQ_{v} factor (see Eqs. 7 and 8) and an alternative shale-volume descriptor, S_{wb}, the saturation of physically bound water in the total PV (see Fig. 1). ^{[1]}^{[14]} When V_{sh} is zero, S_{wb} is zero; and when V_{sh} is 100% BV, S_{wb} and S_{wt} are also 100% PV.
where R_{wb} = resistivity of clay-bound water in the shales, and R_{wf} = resistivity of free formation water in the shale-free water zones. Because of the different model assumptions, DW exponents m_{o} and n_{o} must always be smaller than the WST exponents^{[24]} and may be values similar to "clean" sand exponents. Where the WST and DW models have been properly applied, the hydrocarbon pore volume (V_{HCP}) results should be equal. All S_{wt} calculations from the WST and DW methods must be checked to ensure that they are greater than S_{wb}. After this check, they are used with ϕ_{t} to obtain the V_{HCP}. For the DW model, when the outputs require conversion to effective porosity, ϕ_{e}, and effective water saturation, S_{we}, the properties are converted with Eqs. 9 and 10, respectively.
Strengths and weaknesses of resistivity log-based calculations
The greatest strength of S_{w} calculations from the R_{t} logs is that these calculations can be made at each net-pay depth with valid data for all wells within the log database. The calculations can account for any subsets of input parameters related to the individual zones.
The weaknesses of the R_{t}-based S_{w} calculations are that one has to select a model to describe the relationship of S_{w} to R_{t}, R_{w}, and a variety of other input parameters. Any model is an approximation to the real nature of the reservoir pore system and, typically, has limitations such as how the clay-mineral conductivity is modeled. Log-analysis estimates of V_{sh} are rather uncertain, so sands that are substantially free of clay minerals can easily, and incorrectly, be assigned significant clay volumes. In these circumstances, complex shaly-sand models may have been applied when it is more appropriate to model the sand as clean sand. Effective porosity is also impacted by the uncertain V_{sh} estimates. The R_{w} is often assumed to be constant within the hydrocarbon column, and usually there is little data regarding R_{w} other than from aquifer samples. In several cases in which the R_{w} distribution has been studied in depth, it was found to vary in systematic ways within the hydrocarbon column and not necessarily be the same as in the underlying aquifer. ^{[3]}^{[4]}^{[5]}
CEC can be measured in the laboratory, but in the reservoir it must be estimated by correlations with porosity or V_{sh}. For the laboratory CEC measurements, there are fundamental uncertainties such as the degree to which the clay-mineral geometry is altered by the disaggregation of the core. The total surface area and CEC may be enhanced by comminution (i.e., grinding to grain-size particles). ^{[26]}
The other input parameters for the S_{w}/R_{t} models are either based on "worldwide experience" (such as default exponent parameters in commercial software packages) or developed from SCAL rock-electrical-property measurements on a relatively small number of core plugs from the reservoir interval. Hence, there are relatively few data determining the parameters that are used for the log point-by-point S_{w} calculations. It has to be assumed that the manner in which the water saturation is distributed in the core plugs during these laboratory experiments is like that of the real reservoir. Because water is present during the laboratory measurements, clay minerals are rehydrated at the time of the tests.
Laboratory capillary-pressure/saturation measurements
A second S_{w} method that is totally independent of the resistivity logs uses laboratory-measured P_{c}/S_{w} data. The underlying concept of the use of capillary pressure data is that the reservoir has come to capillary equilibrium over geologic time (the millions of years since hydrocarbons have entered and filled the reservoir trap). This equilibrium is reproduced in laboratory experiments using the centrifuge, porous-plate, and mercury-injection capillary pressure (MICP) methods. The P_{c}/S_{w} data are measured on a selected set of reservoir core plugs representing a range of porosity and permeability values (and possibly also lithologies).
Centrifuge experiments are typically made on 1-in. core plugs over a period of several days in the intense gravitational field (up to 1000 G) of the centrifuge and are assumed to be equivalent to what occurs in a hydrocarbon reservoir over millions of years in a 1-G gravitational field and over lengths of 10 to hundreds of feet. These assumptions are broadly accepted as being reasonable, provided that the samples are not damaged during testing in the centrifuge. The reported P_{c}/S_{w} values are not the raw laboratory data. In the laboratory, the average saturation is determined at each centrifuge speed, and those raw data are input to a mathematical model to convert them to a tabulation of endface saturations and P_{c} values.
Porous-plate P_{c} tests are made on core plugs at several different gas pressures and are generally carried out at the same time as the resistivity experiments. After reaching equilibrium with no further brine flow at each pressure, the S_{w} is constant along each plug and is calculated from its weight loss.
MICP tests are made on dried core pieces and the volume of injected mercury, the nonwetting phase, is converted to an S_{w} value. This is considered to be total S_{w} if, at high enough pressures, mercury enters both the microporosity and dry clay-mineral porosity. Conversely, for centrifuge or porous-plate tests, where brine is present as the wetting phase, clay minerals probably hydrate, and their physically bound water is unlikely to be displaced during the test. Brine-related P_{c}/S_{w} measurements may give total or effective S_{w}, depending on the specific porosity measurement method used (i.e., whether the porosity occupied by the clay-mineral physically-bound water is included or excluded from the porosity calculation). Effective S_{w} values are always lower than total S_{w} values and should be very low at high capillary pressures if there is little nonclay-mineral-related microporosity. ^{[27]}
The conversion of the laboratory P_{c}/S_{w} data to reservoir conditions requires knowledge of the IFT and contact angle of the fluid pair used in the laboratory and properties of the brine and hydrocarbon fluids at reservoir conditions. These are needed to calculate the density of each phase and to estimate the interfacial-tension (IFT) between the fluid pair at reservoir conditions. The P_{c} values (in psi) are converted to vertical height above the hydrocarbon/water contact, H_{hwc} (in feet), with the following formula:
where fluid densities (ρ) are in g/cm^{3}, and the subscripts are r = reservoir, s = surface, h = hydrocarbon, and w = water. Table 1 lists some typical values^{[28]} for IFT, σ, and contact angle, θ, used in 11. ^{[29]} and provides approximate ranges for the factors for converting P_{c}-laboratory data to height above a reservoir free-water level. Height-P_{c} conversion factors are similar for many oil and gas reservoirs; the footnotes in Table 1 describe the values that were assumed to calculate these ranges. More details of correlations for brine/hydrocarbon IFT as a function of oil or gas gravity have been published. ^{[30]} The reservoir-condition contact angle, θ, is usually taken as 0 for gas reservoirs and 0 or 30° for oil reservoirs because, generally, data are not available at reservoir conditions.
The suite of P_{c}/S_{w} data is typically converted to a mathematical relationship between S_{w} as the dependent variable and the independent variables—porosity, permeability, and H_{owc} or H_{gwc}.^{[18]}^{[17]}^{[31]} Because permeability is usually determined as a function of porosity, it is often not included as an independent variable. Two of the mathematical forms that have been used are
where A, B, C, D, and E are curve-fit constants. In Eq. 13, B permits the removal of singularities at zero height.
In developing the coefficients for these relationships, any zonation of the reservoir intervals needs to be applied, and then separate sets of coefficients developed for each zone. The zonation can be based on geological interpretation of the reservoir depositional and diagenetic history and/or variation in the P_{c}/S_{w} curves for different parts of the reservoir interval.
The depth of the reservoir’s oil/water contact (OWC) or gas/water contact (GWC) must be known in order to make S_{w} calculations using the P_{c}/S_{w} methodology. The calculations of S w are made only above this depth. In reality, the H_{owc} or H_{gwc} is referenced to the free water level (FWL) (i.e., the depth at which P_{c} =0 and which is deeper than the observed OWC or GWC). For a gas reservoir consisting of good-quality rocks, the difference between the FWL and the GWC is typically 1 ft or less. However, for an oil reservoir containing a heavier oil, this difference can be 10 feet or more, and, given four-way closure on an anticlinal structure, the impact on the OOIP volume between using the FWL vs. the observed OWC as the H_{owc} = 0 depth can equal a few percent of OOIP.
Once the various sets of coefficients have been developed and the P_{c} to H_{owc} (or H_{gwc}) conversion made, an S_{w} value can be calculated at each data point within the log database that has a valid porosity value and is above the OWC or GWC. Hence, there will be the same number of, or more, S_{w} values available from this S_{w} methodology as when using the R_{t} logs.
Strengths and weaknesses of capillary pressure-based calculations
The strength of S_{w} calculations from P_{c}/S_{w} data is that, after making a correlation with porosity and height, a unique S_{w} value is available for all wells at all net-pay depths with valid porosity values in the log database. This also applies to the whole hydrocarbon column anywhere in the reservoir once the wellbore porosity values have been propagated into the full geocellular model grid. These calculations can account for any zonation and subsets of input parameters related to the individual zones.
A potential weakness in the P_{c} approach to S_{w} calculations is whether the laboratory measurements have been allowed sufficient time to reach equilibrium. If not, the S_{w} values, particularly at high P c values, will be too high. Another potential weakness is the accuracy of the IFT value used in converting from surface to reservoir conditions; fortunately, these values vary over a limited range for most hydrocarbon/brine pairs. A third potential weakness is the definition of the FWL depth compared with the observed OWC or GWC. A fourth potential weakness is whether enough data have been taken to be representative, both vertically and areally, of the zones in the reservoir. ^{[32]}
The fifth potential weakness concerns the complexity of the reservoir’s hydrocarbon-filling and structural history. In simple oil-reservoir situations and most gas-reservoir situations, this is not an issue. However, for oil reservoirs with tar mats and heavy-oil zones there is a complication because of the varying oil density near the OWC, including the possibility that the tar mat has a hydrocarbon density very close to that of the connate brine. Another aspect may be whether all or portions of the hydrocarbon column are on the imbibition cycle where imbibition P_{c}/S_{w} data are needed for the S_{w} calculations, not the typical drainage P_{c}/S_{w} data. ^{[32]}
OBM-core-plug Dean-Stark water-volume determinations
The third method for determining the S_{w} in a reservoir’s hydrocarbon column is to cut OBM cores and perform Dean-Stark water-volume determinations on the routine core plugs. Foot-by-foot S_{w} values can be calculated from these water volumes and the associated core-plug PVs. OBM cores are typically cut only in a few wells in a particular field. These S_{w} data can be applied to other, uncored wells in the reservoir if strong correlations between these values and porosity and/or permeability are identified. These data are not valid in the oil/water or gas/water transition zone or in the aquifer, intervals in which the connate brine is mobile. OBM-core S_{w} values may be found to be either higher or lower than those from the other two methods described previously.
Strengths and weaknesses of OBM-core values
The strength of S_{w} values from routine-core-analysis Dean-Stark S_{w} data is that these data are the most direct measure of reservoir connate S_{w} values above the oil/water or gas/water transition zone. Relative to the two methods discussed previously and the variations of these methods, the OBM-core S_{w} approach is a direct S_{w} determination and the other methods are indirect S_{w}-calculation approaches that require many more assumptions and inferences.
The weaknesses of the OBM S_{w} method are that it does not apply to the lowest parts of the oil/water or gas/oil transition zone where the brine phase has mobility and that, generally, the amount of OBM core S_{w} data is limited because the operator cuts cores with OBM only in a limited number of wells because of the expense. The first of these weaknesses can be overcome if the OBM S_{w} data are used in combination with either the resistivity logs or with P_{c}/S_{w} data.
Another consideration is that the whole project, from the mud formulation to the core-handling and -preservation procedures on through the routine-core-analysis measurements, needs to be monitored and reviewed in detail to ensure all steps were executed properly. This demands that considerable time and effort be spent by the technical team to ensure success; however, to some extent the same comment applies to the P_{c}/S_{w} and resistivity-log/S_{w} calculation approaches discussed previously.
Integration of water saturation data from different methods
Depending on the data availability in a particular reservoir situation, a combination of the various S_{w} approaches may prove superior to the use of a single type of data. The first step in going to a combination approach is to review the reservoir’s database to identify any significant gap in vertical, or areal, coverage. The most obvious gap often occurs near the fluid contact, because there is little reason to drill wells in downdip locations, particularly during a reservoir’s development phase. Three examples of combination approaches are described below.
Resistivity-log and capillary pressure data
Resistivity-log-derived S_{w} results may not be available throughout the hydrocarbon column of a reservoir. To fill gaps and average the point-by-point data set, it is common practice to plot S_{w} as a function of height, to omit nonpay points, and to identify various porosity ranges by coding the data points. Resistivity-log-derived S_{w} data frequently shows V- or U-shaped patterns on these plots because of the shoulder/bed effects near nonpay sections (shales). The most accurate S_{w} values in such patterns are usually at the lowest S_{w} values where the thin-bed correction is minimized. In a manner similar to that described in the previous P_{c}/S_{w} section, height/saturation curves are often fitted to these resistivity-log-derived S_{w} data to enable reservoir hydrocarbons-in-place volumes to be calculated. The function forms are similar to or are the same as those described above for P_{c}/S_{w}. ^{[31]}
Routine OBM core with capillary pressure data
Because there is a need to define the S_{w} characteristics of the oil/water or gas/water transition zone and because the OBM-core S_{w} data can be incorrect and too low in this interval, one approach is to use P_{c}/S_{w} data in combination with the routine OBM-core S_{w} data. This can be done by first correlating the OBM-core S_{w} data to porosity and assuming that this relationship is valid above the oil/water or gas/water transition zone. The functional form of this first relationship might be
The second step is to create a tabular data set in which the S_{w}/porosity correlation is used to calculate an array of S_{w} values for large H_{owc} or H_{gwc} values and a range of porosity values. For this part of the data set, S_{w} is assumed to be independent of the H_{owc} or H_{gwc} values. The P_{c}/S_{w} data converted to reservoir conditions is used to provide data points for low H_{owc} or H_{gwc} values and various porosity values. Statistical calculations are applied to the whole of this data set. The functional form of this second relationship might be
With this functional form, the boundary conditions of the first step are automatically met in the second step.
Routine OBM-core with resistivity-log data
To address the lack of valid OBM-core S_{w} data in the oil/water or gas/water transition zone discussed previously, it is also possible to combine OBM-core S_{w} data with the resistivity-log data to develop an overall S_{w} methodology. This approach assumes that a number of wells have been drilled through the OWC or GWC so that there are log resistivity values through the oil/water or gas/water transition zone. In this approach, the OBM-core S_{w} data are used to back-calculate the saturation-exponent, n, values over each zone so that the core-based V_{HCP} value equals that calculated from the resistivity logs (see Eq. 2). Then the core-based saturation-exponent, n, values are applied to the noncored well’s resistivity logs to calculate S_{w} point-by-point throughout the reservoir interval in all wells. ^{[29]}^{[33]} This approach assumes that the R_{w}, a, and m values have been determined from other experimental and fluid-sample data so that R_{0} can be calculated.
where R_{0} is the bulk resistivity at S_{w} = 100% PV and is calculated with Eqs. 4 and 5. R_{t} is the deep-reading resistivity-log reading, and S_{wc} is the OBM-core S_{w} above the mobile-water transition zone. The resulting back-calculated n values at the core-plug depths are averaged for the zone. In some instances, n may be found to have an areal variation within a zone that should be taken into account in subsequent calculations.
Adjustments to water saturation data from different methods
We have described three methodologies water saturation determination. They are basically independent methods; hence, they can be used together to determine the accuracy of the S_{w} calculations throughout the hydrocarbon column. Because the methods are based on very different technical approaches and assumptions, if the different methods give essentially the same S_{w} answer, then it is highly likely that this is the correct S_{w}.
However, the challenge comes when, as is often the case, the different methods result in different S_{w} values and distributions. The OBM-core S_{w} values might be either higher or lower than those from the other two methods. The common misunderstanding that OBM-core S_{w} is likely to be too low is unsubstantiated. In a very large reservoir, it could go both ways depending on where one is in the reservoir. ^{[29]}^{[33]} If the values are quite different, two aspects of the calculations need to be reviewed in depth. First, the quality of the input laboratory data needs to be checked and how it was converted from raw data into the input values to the S_{w} calculations needs to be reviewed. Second, the assumptions and models used for the S_{w} calculations need to be checked. For example, with the P_{c}/S_{w} data, the assumed oil/water density difference may be considerably in error, or the shaly-sand S_{w} model may be inappropriate for the particular reservoir. As well as the S_{w} averages, the zone-average V_{HCP} values from the various methods should be compared, which includes porosity in the comparison calculations.
Core, total, and effective systems compatibility
The Archie R_{t}-based S_{w} equation models "clean" sands. Various other shaly-sand models use either the effective or the total-porosity systems. It is well known that these basic models, if applied properly to the same formation, must produce the same final V_{HCP} from their different calculation procedures (see Fig. 1 and Eqs. 6 through 8). ^{[1]} ϕ_{t} is greater than or equal to ϕ_{e}; however, at the same time, S_{wt} is greater than or equal to S_{we} and, when used together, the appropriate combinations must give the same V_{HCP} result. For the total-porosity system, V_{HCP} = ϕ_{t} × (1-S_{wt}), whereas for the effective-porosity system, V_{HCP} = ϕ_{e} × (1-S_{we}).
The V_{HCP} can also be estimated from a combination of core porosity and Dean-Stark S_{w} measured on preserved OBM cores. The several systems—core, total, and effective—must all give the same fundamental results, and the most accurate of them (the OBM-core method) can be used to calibrate and test the less accurate methods. When properly adjusted and applied (e.g., by improving the V_{sh} estimates or IFT values), all three methods give the same final V_{HCP}. If they do not agree, the likely sources of uncertainty and error must be examined.
It is clearly inconsistent and incorrect to mix the systems by, for example, reporting an effective porosity with a total S_{w}, a total porosity with an effective S_{w}, or a standard-core porosity with an effective S_{w}. System compatibility must also be maintained by correct use of the SCAL measurements and log-analysis formulae, when these are used to calibrate the resistivity logs and P_{c}/S_{w} methods. The differences should be resolved as much as possible. To the extent that they are not, the differences can be considered to be a measure of the uncertainty in the S_{w} calculations.
Uncertainties
It is the uncertainty of the hydrocarbon saturation (1 − S_{w}) that is economically important, not the absolute uncertainty in S_{w}. When uncertainties in S_{w} are evaluated, their importance in terms of S_{o} and S_{g} should be accounted for. The uncertainties of the several S_{w}-evaluation methods vary widely.
OBM-core water saturation data
The water volume extracted from a single core plug may have a random and known systematic uncertainty of ± 0.05 cm^{3}, where each uncertainty refers to one SD. The PV of a typical 1-in. core plug is 4.0 cm^{3} if the porosity is 20% BV. The water-volume uncertainty alone equates to an S_{w} uncertainty of ± 1% PV (0.05/4.0). The uncertainties in porosity have a further effect on this calculation. ^{[19]} An OBM-core S_{w} of 20% PV, therefore, has a combined 1-SD range from approximately 18 to 22% PV. At lower porosities and higher S_{w} values, the water-volume uncertainty may be ± 0.1 cm^{3}, leading to an S_{w} uncertainty of ± 3% PV, when the porosity is 15% BV. As porosity decreases, the uncertainty grows. Before the measurements are made, any water in the toluene and the Dean-Stark apparatus must be removed, or the S_{w} values will be overstated. The extraction time required to recover the water adsorbed on the clay minerals adds to the uncertainty.
The uncertainty of the average core S_{w} will be improved when plugs are selected at one or two per foot with equal spacing and without regard for the rock quality. However, as discussed previously, plug samples are not always selected at random, so care must be exercised, especially regarding the S_{w} values predicted at depths where core is not available. From a broader perspective, it must also be remembered that 1-in. core plugs only sample approximately 2% of the full-core volume. Because of these many factors, the authors estimate that uncertainties similar to those given concerning porosity also propagate to the zone-average OBM-core S_{w} values. Measurements in which larger core plugs are analyzed will reduce several of the uncertainties.
Resistivity-log-derived water saturation values
The log readings, typical SCAL-derived Archie exponents, and all of the other associated parameters are uncertain. For example, the resistivity-log uncertainty may be ± 50% when R t is 500 ohm•m. The most important uncertainty contributors at low S_{w} values are likely to be R_{t} and n. S_{w} uncertainty in this circumstance is estimated at ± 5% PV (i.e., if S_{w} is calculated as 10% PV, the 1-SD range is 5 to 15% PV). ^{[19]} At lower porosity values and higher water saturations, similar methods led to uncertainty estimates of ± 9% PV. Given that further uncertainty in the final calculated S_{w} may arise from shaly-sand effects and many other sources, the authors believe that the ranges given apply equally to the overall systematic uncertainty of the S_{w} zone-average values. These estimates are all 1 SD; therefore, in 32% of cases, zone-average uncertainties are considered likely to be greater than the ranges given.
Capillary pressure-derived water saturation values
The uncertainty estimates are the sum of several factors. Most of these factors have their greatest impact on the S_{w} calculations in the first 100 to 200 ft of the hydrocarbon column above the fluid contact. Therefore, because the transition zone is considerably longer in many oil reservoirs than in a gas reservoir, their impacts will be greater in most oil reservoirs. Above 200 ft, the S_{w} values are usually only changing slowly; hence, the primary consideration above the transition zone is whether the laboratory measurements are taken at equilibrium conditions.
The first factor in the uncertainty analysis is the fundamental assumption as to whether the drainage or imbibition P_{c}/S_{w} data should be used. In most cases, the drainage curves should be used, but, in a few situations, the reservoir may be on the imbibition cycle. In these situations, the improper choice of using the drainage P_{c} curve can lead to a +5 to 20% PV S_{w} error in the first 100 to 200 ft above the OWC. ^{[34]}^{[35]}
The second factor concerns the laboratory P_{c}/S_{w} measurements. If the measurements are not taken to equilibrium, then the S_{w} values at a particular P_{c} value will be too high. This can be +1 to 10% PV effect for the large H_{owc} or H_{gwc} range. The other key aspects for reported centrifuge laboratory results are how the raw laboratory measurements of water volumes were determined and how these data have been converted to the reported endface saturations. The water-volume measurements have the same-size potential error as discussed for the OBM Dean-Stark S_{w} measurements (± 1 to 3% PV). Differences in the laboratory calculation procedures can result in further variations of ± 1 to 3% PV in reported P_{c}/S_{w} results when using the same raw laboratory data. For porous-plate tests and others, the repeated handling of poorly cemented or uncemented core plugs can cause grain loss, which, after the final calculations, translates into small errors in S_{w}.
The third factor is how the suite of raw laboratory data for a particular reservoir interval are curve-fitted and presented in the final laboratory report as tabulations of P_{c}/S_{w} values for each core plug. Uncertainty in the application arises from how these reported values are averaged for use in the S_{w} calculations over the full range of reservoir porosity and permeability values. This uncertainty includes how the data are weighted and whether some potential outlier data from one or two core plugs distorts the averaged P_{c}/S_{w} curves. These uncertainties primarily affect the first 100 to 200 ft above the H_{owc} or H_{gwc} so that their impact depends on how thick the hydrocarbon column is and its distribution as a function of H_{owc} or H_{gwc}.
The final factor is the conversion of the averaged P_{c}/S_{w} curves (or equation) from surface to reservoir conditions, all of which affect the conversion of P_{c} values to H_{owc} or H_{gwc} values. This includes a number of subfactors, each with its own uncertainty level: IFT at surface and reservoir conditions, fluid-pair density difference at reservoir conditions, contact angles, and depth of the actual in-situ FWL compared with the OWC or GWC. The contact angles at surface and reservoir conditions are generally taken to be the same because no data are available to proceed otherwise. For these other factors, the uncertainty is considerably greater for an oil reservoir than for a gas reservoir; because the IFT values can be low and compared with those for a gas reservoir, the density differences are significantly less particularly if there is a vertical oil-gravity variation that results in a heavy-oil interval just above the OWC. All of these factors affect the H_{owc} or H_{gwc} values; therefore, their impact on the S w calculations is predominantly in the first 100 to 200 ft above the fluid contact.
In summary, the use of P_{c}/S_{w} data can result in S_{w} uncertainty of ± 5 to 15% PV in the oil/water or gas/water transition zone. Above that transition zone, the uncertainty is related to whether the laboratory data were taken at equilibrium conditions and how the various P_{c}/S_{w} curves have been averaged together. In this range, the uncertainty is likely to be 3 to 10% PV.
Nomenclature
a | = | Archie cementation constant |
a* | = | Waxman-Smits cementation constant |
A | = | Coefficient in various equations of this chapter |
B | = | Specific cation conductance, [(1/ohm•m) / (meq/mL)] |
C | = | Coefficient in various equations |
D | = | Coefficient in various equations |
E | = | Coefficient in various equations |
F | = | Archie formation factor |
F* | = | Waxman-Smits-Thomas formation factor |
F_{HCP} | = | hydrocarbon pore feet, L, ft [m] |
H_{gwc} | = | height above the gas/water contact, L, ft [m] |
H_{hwc} | = | height above the hydrocarbon/water contact, L, ft [m] |
H_{owc} | = | height above the oil/water contact, L, ft [m] |
I_{R} | = | resistivity index |
J(S_{w}) | = | Leverett J-function |
k | = | permeability, L^{2}, md [μm^{2}] |
m | = | Archie cementation exponent |
m* | = | Waxman-Smits-Thomas cementation exponent |
m_{o} | = | dual-water cementation exponent |
n | = | Archie saturation exponent |
n* | = | Waxman-Smits-Thomas saturation exponent |
n_{o} | = | dual-water saturation exponent |
P_{c} | = | capillary pressure, m/Lt^{2}, psi |
P_{ce} | = | entry capillary pressure, m/Lt^{2}, psi |
Q_{v} | = | cation-exchange capacity of total PV, meq/mL |
r | = | correlation coefficient |
R_{0} | = | rock resistivity with 100% PV water saturation, ohm•m |
R_{sd} | = | clean-sand resistivity, ohm•m |
R_{sh} | = | shale resistivity, ohm•m |
R_{t} | = | true resistivity of uninvaded, deep formation, ohm•m |
R_{w} | = | connate-brine resistivity, ohm•m |
R_{wb} | = | clay-bound water resistivity, ohm•m |
R_{wf} | = | free-formation-water resistivity, ohm•m |
R_{xo} | = | shallow-reading invaded-zone microresistivity, ohm•m |
S_{g} | = | gas saturation, %PV |
S_{o} | = | oil saturation, %PV |
S_{orw} | = | residual-oil saturation to water displacement, %PV |
S_{w} | = | water saturation, %PV |
S_{wb} | = | saturation of clay-bound water in the total porosity, %PV |
S_{wc} | = | connate water saturation, %PV |
S_{wc} | = | core water saturation, %PV |
S_{we} | = | water saturation of the effective porosity, %PV |
S_{wsd} | = | sand water saturation, %PV |
S_{wt} | = | water saturation of the total porosity, %PV |
V_{cl} | = | clay content, %BV |
V_{HCP} | = | hydrocarbon pore volume, L^{3}, ft^{3} [m^{3}] |
V_{sh} | = | shale content, %BV |
θ | = | contact angle, degrees |
ρ_{b} | = | formation bulk density, m/L^{3}, g/cm^{3} |
ρ_{fl} | = | fluid density, m/L^{3}, g/cm^{3} |
ρ_{h} | = | hydrocarbon density, m/L^{3}, g/cm^{3} |
ρ_{ma} | = | matrix or grain density, m/L^{3}, g/cm^{3} |
ρ_{w} | = | water density, m/L^{3}, g/cm^{3} |
σ | = | interfacial tension, m/t^{2}, dynes/cm |
ϕ | = | porosity, %BV |
ϕ_{c} | = | core porosity, %BV |
ϕ_{cl} | = | clay porosity, %BV |
ϕ_{e} | = | effective porosity, %BV |
ϕ_{sd} | = | sand porosity, %BV |
ϕ_{sh} | = | shale porosity, %BV |
ϕ_{t} | = | total porosity, %BV |
References
- ↑ ^{1.0} ^{1.1} ^{1.2} ^{1.3} ^{1.4} ^{1.5} ^{1.6} Log Interpretation Principles/Applications. 1989. Houston, Texas: Schlumberger.
- ↑ Log Interpretation Charts. 2000. Sugar Land, Texas: Schlumberger.
- ↑ ^{3.0} ^{3.1} ^{3.2} McCoy, D.D., Jr., H.R.W., and Fisher, T.E. 1997. Water-Salinity Variations in the Ivishak and Sag River Reservoirs at Prudhoe Bay. SPE Res Eng 12 (1): 37-44. SPE-28577-PA. http://dx.doi.org/10.2118/28577-PA.
- ↑ ^{4.0} ^{4.1} ^{4.2} ^{4.3} ^{4.4} Rathmell, J., Atkins, L.K., and Kralik, J.G. 1999. Application of Low Invasion Coring and Outcrop Studies to Reservoir Development Planning for the Villano Field. Presented at the Latin American and Caribbean Petroleum Engineering Conference, Caracas, Venezuela, 21-23 April 1999. SPE-53718-MS. http://dx.doi.org/10.2118/53718-MS.
- ↑ ^{5.0} ^{5.1} ^{5.2} Rathmell, J.J., Bloys, J.B., Bulling, T.P. et al. 1995. Low Invasion, Synthetic Oil-Base Mud Coring in the Yacheng 13-1 Gas Reservoir for Gas-in-Place Calculation. Presented at the International Meeting on Petroleum Engineering, Beijing, China, 14-17 November 1995. SPE-29985-MS. http://dx.doi.org/10.2118/29985-MS.
- ↑ Barber, T.D. 1985. Introduction to the Phasor Dual Induction Tool. J Pet Technol 37 (9): 1699-1706. SPE-12049-PA. http://dx.doi.org/10.2118/12049-PA.
- ↑ ^{7.0} ^{7.1} ^{7.2} ^{7.3} ^{7.4} Woodhouse, R. and Warner, H.R. 2005. Sw and Hydrocarbon Pore Volume Estimates in Shaly Sands - Routine Oil-Based-Mud Core Measurements Compared With Several Log Analysis Models. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 9-12 October 2005. SPE-96618-MS. http://dx.doi.org/10.2118/96618-MS.
- ↑ ^{8.0} ^{8.1} ^{8.2} ^{8.3} Archie, G.E. 1942. The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics. Trans. of AIME 146 (1): 54-62. http://dx.doi.org/10.2118/942054-G.
- ↑ ^{9.0} ^{9.1} ^{9.2} ^{9.3} Poupon, A. and Leveaux, J. 1971. Evaluation of Water Saturations in Shaly Formations. The Log Analyst 12 (4).
- ↑ Winsauer, W.O., Shearin H.M., Masson P.H., and Williams M. 1952. Resistivity of Brine Saturated Sands in Relation to Pore Geometry. AAPG Bull. 36 (2): 253-277.
- ↑ ^{11.0} ^{11.1} ^{11.2} ^{11.3} Waxman, M.H. and Smits, L.J.M. 1968. Electrical Conductivities in Oil-Bearing Shaly Sands. SPE J. 8 (2): 107–122. SPE-1863-PA. http://dx.doi.org/10.2118/1863-PA.
- ↑ ^{12.0} ^{12.1} ^{12.2} ^{12.3} Waxman, M.H. and Thomas, E.C. 1974. Electrical Conductivities in Shaly Sands-I. The Relation Between Hydrocarbon Saturation and Resistivity Index; II. The Temperature Coefficient Of Electrical Conductivity. J Pet Technol 26 (2): 213-225. SPE-4094-PA. http://dx.doi.org/10.2118/4094-PA.
- ↑ ^{13.0} ^{13.1} ^{13.2} ^{13.3} Keelan, D.K. and McGinley, D.C. 1979. Application of Cation Exchange Capacity in a Study of the Shannon Sand of Wyoming. Paper KK presented at the 1979 SPWLA Annual Symposium, June.
- ↑ ^{14.0} ^{14.1} ^{14.2} ^{14.3} Juhasz, I. 1979. The Central Role of Q v and Formation Water Salinity in the Evaluation of Shaly Formations. The Log Analyst 20 (4).
- ↑ Hill, H.J., Shirley, O.J., and Klein, G.E. 1979. Bound Water in Shaly Sands—Its Relation to Q v and Other Formation Properties. The Log Analyst 20 (3): 3.
- ↑ Leverett, M.C. 1941. Capillary Behavior in Porous Media. Trans., AIME 142: 152.
- ↑ ^{17.0} ^{17.1} Heseldin, G.M. 1974. A Method of Averaging Capillary Pressure Curves. The Log Analyst 4 (3).
- ↑ ^{18.0} ^{18.1} Johnson, A. 1987. Permeability Averaged Capillary Data: A Supplement to Log Analysis in Field Studies. Paper EE presented at the 1987 SPWLA Annual Symposium, London, June.
- ↑ ^{19.0} ^{19.1} ^{19.2} Woodhouse, R. 1998. Accurate Reservoir Water Saturations from Oil-Mud Cores: Questions and Answers from Prudhoe Bay and Beyond. The Log Analyst 39 (3): 23.
- ↑ Richardson, J.G., Holstein, E.D., Rathmell, J.J. et al. 1997. Validation of As-Received Oil-Based-Core Water Saturations From Prudhoe Bay. SPE Res Eng 12 (1): 31-36. SPE-28592-PA. http://dx.doi.org/10.2118/28592-PA.
- ↑ Dawe, B.A. and Murdock, D.M. 1990. Laminated Sands: An Assessment of Log Interpretation Accuracy by an Oil-Base Mud Coring Program. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 23-26 September 1990. SPE-20542-MS. http://dx.doi.org/10.2118/20542-MS.
- ↑ Egbogah, E.O. and Amar, Z.H.B.T. 1997. Accurate Initial / Residual Saturation Determination reduces Uncertainty in Further Development and Reservoir Management of the Dulang Field, Offshore Peninsular Malaysia. Presented at the SPE Asia Pacific Oil and Gas Conference and Exhibition, Kuala Lumpur, Malaysia, 14-16 April 1997. SPE-38024-MS. http://dx.doi.org/10.2118/38024-MS.
- ↑ Worthington, P.F. 1985. The Evolution of Shaly-Sand Concepts in Reservoir Evaluation. The Log Analyst 23 (1).
- ↑ ^{24.0} ^{24.1} ^{24.2} Clavier, C., Coates, G., and Dumanoir, J. 1984. Theoretical and Experimental Bases for the Dual-Water Model for Interpretation of Shaly Sands. SPE J. 24 (2): 153-168. SPE-6859-PA. http://dx.doi.org/10.2118/6859-PA.
- ↑ ^{25.0} ^{25.1} Best, D.L., Gardner, J.S., and Dumanoir, J.L. 1979. A Computer-Processed Wellsite Log Computation. Paper Z presented at the 1979 SPWLA Annual Symposium.
- ↑ Huff, G.F. 1987. A Correction for the Effect of Comminution on the Cation Exchange Capacity of Clay-Poor Sandstones. SPE Form Eval 2 (3): 338-344. SPE-14877-PA. http://dx.doi.org/10.2118/14877-PA.
- ↑ Bryant, W.T. and Robert B. Truman, I. 2002. Proper Core-Based Petrophysical Analysis Doubles Size of Ha'py Field. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 29 September-2 October 2002. SPE-77638-MS. http://dx.doi.org/10.2118/77638-MS.
- ↑ Fundamentals of Rock Properties. 2002. Aberdeen: Core Laboratories UK Ltd.
- ↑ ^{29.0} ^{29.1} ^{29.2} Holstein, E.D. and Warner, J., H. R. 1994. Overview of Water Saturation Determination For the Ivishak (Sadlerochit) Reservoir, Prudhoe Bay Field. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 25-28 September 1994. SPE-28573-MS. http://dx.doi.org/10.2118/28573-MS.
- ↑ Katz, D.L. and Firoozabadi, A. 1978. Predicting Phase Behavior of Condensate/Crude-Oil Systems Using Methane Interaction Coefficients. J Pet Technol 30 (11): 1649–1655. SPE-6721-PA. http://dx.doi.org/10.2118/6721-PA.
- ↑ ^{31.0} ^{31.1} Harrison, B. and Jing, X.D. 2001. Saturation Height Methods and Their Impact on Volumetric Hydrocarbon in Place Estimates. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 30 September-3 October 2001. SPE-71326-MS. http://dx.doi.org/10.2118/71326-MS.
- ↑ ^{32.0} ^{32.1} Richardson, J.G. and Holstein, E.D. 1994. Comparison of Water Saturations from Capillary Pressure Measurements with Oil-Based-Mud Core Data, Ivishak (Sadlerochit) Reservoir, Prudhoe Bay Field. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 25-28 September 1994. SPE-28593-MS. http://dx.doi.org/10.2118/28593-MS.
- ↑ ^{33.0} ^{33.1} McCoy, D.D. and Grieves, W.A. 1997. Use of Resistivity Logs To Calculate Water Saturation at Prudhoe Bay. SPE Res Eng 12 (1): 45-51. SPE-28578-PA. http://dx.doi.org/10.2118/28578-PA.
- ↑ Lucia, F.J. 2000. San Andres and Grayburg Imbibition Reservoirs. Presented at the SPE Permian Basin Oil and Gas Recovery Conference, Midland, Texas, 21-23 March 2000. SPE-59691-MS. http://dx.doi.org/10.2118/59691-MS.
- ↑ Thai, B.N., Hsu, C.F., Bergersen, B.M. et al. 2000. Denver Unit Infill Drilling and Pattern Reconfiguration Program. Presented at the SPE Permian Basin Oil and Gas Recovery Conference, Midland, Texas, 21-23 March 2000. SPE-59548-MS. http://dx.doi.org/10.2118/59548-MS.
Noteworthy papers in OnePetro
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read
External links
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro
See also
Log analysis in shaly formations