The equations of the “Formation evaluation from well logs” lessons with their conditions and units, and what each symbol means in this discipline.
Sampling interval of LWD curves transmitted in real time
Δz=3600Rvb
v in m/h, b in bits, R in bit/s, Δz in m; 3600 is the number of seconds in an hour
Holds when: The whole channel carries the logging curves. In practice the directional survey and service data take part of it, and the interval is then larger. The rate of penetration and the telemetry rate are constant. We ignore data compression.
Reading in the middle of a thin bed (uniform-window model)
xc=Lhxb+(1−Lh)xs,h<L
any units of the measured quantity, the same for all; h and L in metres
Holds when: The tool averages the property evenly over a window of length L. The real response is smoother: the middle of the window counts for more than the edges. The bed lies between thick enclosing rocks that are identical. The quantity mixes linearly. This is reasonable for GR and density, but not for resistivity, where the mixing depends on the type of sonde. For h≥L the reading in the middle of the bed is xc=xb.
Statistical noise at a different logging speed and after averaging
σ2=σ1v1v2,σk=kσ1
noise in the units of the curve (g/cm³, API, or porosity units: one porosity unit is 0.01, or 1%, of porosity); speeds in the same units
Holds when: Counting statistics set the noise: the number of pulses per sample is proportional to Δz/v. We ignore other interference. Adjacent samples are independent and the sampling interval Δz is the same. The bed is uniform over the averaged samples.
Holds when: The tool and the filter are uniform windows, as in the thin-bed model (vertical resolution). We replace the two windows by one window with the same second moment (a measure of spread): when two averagings are applied one after the other, their variances (squared spreads) add. Compared with the exact calculation of two windows, the approximation slightly exaggerates the loss of contrast (in the example below, 0.149 against 0.152).
Holds when: Both recordings are depth-matched. We choose the depths in uniform parts of beds, not at boundaries. The errors of the two recordings are independent and equal: the difference then has a standard deviation of σ2. The standard error of the mean difference is sΔ/N. A mean that is several times larger than this error is a systematic offset.
Holds when: Between markers the shift changes smoothly (cable stretch, thermal elongation). Interpolation does not correct jerks of the tool at ledges in the hole. Above the uppermost marker and below the lowest one we take the shift as constant, equal to the shift at the nearest marker.
Depth of an LWD sensor from the time of the sample
zs=zb(t)−ds,δz=vδt
m, h
Holds when: We know the bit depth from the surface record. Its own errors (pipe stretch, heave at sea: the up-and-down motion of the rig) are not included here. The second formula holds as long as the rate of penetration is constant within the clock disagreement.
Holds when: The window contains distinct beds. On a uniform stretch r hardly depends on the shift, and the maximum cannot be trusted. The shift is constant within the window. The curves were recorded by tools of similar resolution.
Diameter of the invaded zone from the volume of filtrate
di=dh2+πhϕSxo4Vf
any consistent system (here SI)
Holds when: Piston-like displacement (the filtrate pushes the fluids ahead of it like a piston): beyond the filtrate front the formation is unchanged, and we neglect the transition zone. Invasion is uniform around the circumference and over the thickness of the bed, without gravity segregation. The filtrate displaces the formation water and the moveable hydrocarbons rather than mixing with them. We take the volume of filtrate as given. In reality the mud, the cake, the overbalance and time determine it.
Ratio of the resistivities of the flushed zone and the uninvaded formation
RtRxo=RwRmf(SxoSw)2
any, the same for all the resistivities (usually Ω·m)
Holds when: A clean formation with intergranular porosity (pores between the grains): only the water in the pores conducts current. Archie's law with a saturation exponent n=2 and the same formation resistivity factor in both zones. A step profile. We obtained Rxo and Rt after corrections for the borehole, the cake and invasion. Rmf and Rw have been brought to formation temperature.
Holds when: Filtration at rest, without circulation. The cake is incompressible, and its permeability and composition are constant. The resistance of the formation is small compared with the resistance of the cake. The cake is thin compared with the borehole radius (a planar problem). We neglect the instant initial filtration (spurt loss). In practice an initial volume Vsp is added to the formula.
Holds when: A step invasion profile. Corrections for the borehole and the enclosing rocks have already been applied. A conductive mud (the laterolog does not work in oil-based muds).
as indicated; Cg/L/ρ is grams of salt per kilogram of solution
Holds when: The density and the salinity refer to the same temperature. The density formula is an approximation for NaCl solutions at room temperature. For waters of complex composition we measure the density.
Holds when: The gradient is constant over the section. This is a good approximation for a uniform sedimentary sequence without thick salt layers and without active groundwater movement. TTD is the stabilised formation temperature, not the temperature measured immediately after circulation.
Recovery of the bottom-hole temperature after circulation stops
Tbh(Δt)=Tres−mlgΔttc+Δt
°C; tc and Δt in the same units; lg is the decimal logarithm
Holds when: The cooling of the bottom of the hole is like the removal of heat by a line source of constant power (a thin line that takes away heat) during tc. The readings are taken at one depth, and the time Δt is not too short (usually more than a few hours). We need at least three readings (two points cannot show whether the data lie on a line), recorded together with the time circulation stopped.
Resistivity of an NaCl solution at 24 °C from its concentration
R24≈0.0123+Cppm0.9553647.5
empirical: R in Ω·m, C in ppm
Holds when: An empirical approximation (a formula fitted to tables of measurements) for NaCl solutions, from hundreds of ppm to nearly saturated. The error is usually within a few per cent, and larger for nearly saturated brines. For a water of another composition, C is the equivalent NaCl concentration.
Converting the resistivity of a solution to another temperature (the Arps formula)
R2=R1T2+21.5T1+21.5
°C; the resistivities in any identical units
Holds when: An empirical approximation for NaCl solutions in the usual range of formation temperatures. In logging we apply it to all aqueous solutions: formation water, filtrate and drilling mud. For mudcake and for waters of complex composition the accuracy is lower.
How an error in saturation leads to an error in the hydrocarbon pore thickness of a bed
HδH=−1−SwδSw=−1−SwSw⋅SwδSw
dimensionless (relative changes)
Holds when: The thickness and porosity of the bed do not change; only Sw does. The error is small (linear approximation). The bed stays on the same side of the cut-offs. Otherwise its pore thickness changes in a jump, all at once.
Sonde reading as a mixture of the responses of the rock components
L=i∑ViLi,i∑Vi=1
units of the method: g/cm³ for density, barn/cm³ for the volumetric photoelectric absorption U, µs/m for interval transit time, fraction for the hydrogen index
Holds when: The sonde sees a uniform volume, the bed is thicker than the vertical resolution, and the borehole corrections have been applied. For bulk density the rule is exact, because it is a mass balance. For photoelectric absorption, it is not Pe that mixes but the product U=Peρe. For the neutron method the rule is approximate: the matrix mineral, neutron absorbers and gas all have an effect. For interval transit time it is the empirical time-average equation, valid for compacted, cemented rocks. For resistivity the rule does not hold: current flows only through the water in the pores, and the relationship is nonlinear (Archie's equation, Archie's equation and saturation of clean reservoirs).
Gamma ray reading and the content of radioactive minerals
GR=ρb∑iρiViAi
API; densities in any units, the same for all terms
Holds when: Corrections for the hole diameter, the mud and the casing have already been applied. The contributions of the minerals are independent and add up. The bed is uniform within the volume the tool sees and thicker than the tool's vertical resolution.
any consistent gamma ray units; the index is dimensionless
Holds when: The radioactivity of the bed is created only by shaly material. The shales in the bed are the same as in the reference shale packages, and the clean part is the same as in the reference clean bed. The baselines are chosen in the same hole, on the same record and close in depth.
Total gamma ray from the potassium, uranium and thorium contents (approximate)
GR≈16K+8U+4Th
GR in API; potassium in mass per cent, uranium and thorium in g/t
Holds when: The coefficients are approximate and depend on the tool and its calibration. Radioactive equilibrium holds in the uranium and thorium series. Borehole corrections have been applied.
Holds when: The count rate is constant within a sample (a uniform bed). The background radiation and the dead time of the detector (the short pause after each count, when it cannot count) are ignored. The samples are independent, with no filtering over depth.
mV; resistivities in any consistent units; lg is the base-10 logarithm
Holds when: A clean (shale-free) water-bearing bed between thick shales that behave as ideal cation membranes. The solutions are close in composition to an NaCl solution; with a large share of other salts the coefficient at 25 °C differs from 71 mV. The streaming potentials are small or cancel each other. Both resistivities and the coefficient K refer to formation temperature.
Holds when: This is a qualitative equivalent circuit. The currents in the bed and in the shales spread in a complicated way, and we do not measure rsd and rsh directly. The electrochemical EMF is concentrated at the shale–mud and filtrate–formation water boundaries.
Equivalent formation water resistivity from the static SP
Rwe=Rmfe⋅10SSP/K
Ω·m; SSP and K in mV
Holds when: The SSP was taken in a thick clean water-bearing bed or corrected for the thickness and resistivity of the bed. Rmfe and K refer to formation temperature. The water and the filtrate are close to NaCl solutions, and the streaming potentials are small.
Electrochemical coefficient in terms of the transport numbers
K=(1+t−−t+)F2.303RT=2t−F2.303RT
SI: K comes out in volts, multiply by 1000 for millivolts; 2.303=ln10 converts the natural logarithm to the base-10 one
Holds when: The shales are an ideal membrane: they pass only cations. The bed is clean and is not itself a membrane. Both solutions are solutions of the same salt, NaCl; the transport numbers are constant over the range of concentrations considered.
Holds when: The factor 106 converts seconds to microseconds. The velocity and the interval transit time refer to the same wave, the compressional one unless stated otherwise.
Velocities of the compressional and shear waves in an elastic medium
Vp=ρbK+34G,Vs=ρbG
SI: moduli in Pa, density in kg/m³, velocities in m/s
Holds when: A uniform, isotropic medium (the same in all directions) that behaves elastically (it springs back after it is deformed). The wavelength is much greater than the grains and pores, and the rock behaves as a continuous medium with averaged properties. For a liquid and a gas G=0: there is no shear wave in them, and Vp=K/ρ, where K and ρ are the bulk modulus and the density of the fluid itself.
Travel time of the head wave from the transmitter to the receiver
t=LΔt+(aT+aR)Δtmcosθc,sinθc=ΔtmΔt
times in µs, distances in m, interval transit times in µs/m
Holds when: The ray approximation: the wave travels along straight segments, from the transmitter to the wall, along the wall and from the wall to the receiver. The rock is faster than the mud (Δt<Δtm); otherwise there is no head wave. The wall in the path of the wave is smooth and parallel to the tool, and the rock along it is uniform.
Holds when: Both receivers trigger on the same oscillation of the same wave, the first arrival of the compressional head wave. The gaps between the tool and the wall at both receivers are equal. The rock is uniform within the span: the measured Δt is the average over an interval of length s.
Travel time through an interval and the average velocity
t=i∑Δtihi,Vˉ=t106∑ihi
Δt in µs/m, h in m, t in µs, Vˉ in m/s (the factor 106 converts microseconds to seconds)
Holds when: The wave travels vertically; in a deviated well, vertical thicknesses are used. The velocity at logging frequencies (kilohertz) is usually somewhat higher than at seismic frequencies (tens of hertz), and the tool sees only the rock at the wall. So we calibrate the integrated time against a check-shot survey or vertical seismic profiling (VSP). Both are direct measurements of the travel time from the surface to a receiver in the well.
Dynamic shear modulus and Poisson's ratio from the sonic log
G=Δts2ρb,ν=R2−121R2−1,R=ΔtpΔts
SI: density in kg/m³, interval transit time in s/m, modulus in Pa
Holds when: A uniform isotropic elastic rock. The moduli are dynamic: they describe the small, fast deformations in a sound wave. Static moduli, which are measured by loading core slowly, are usually lower. For engineering calculations we obtain them from relationships checked on core. Δts is measured, not estimated from Δtp.
Detector count and electron density (calibration relationship)
lnN=a−bρe
N in counts/s (natural logarithm), ρe in g/cm³
Holds when: An approximate relationship found by experiment, for a fixed tool design, energy window and pad contact. The distance to the detector is greater than a few photon mean free paths (the average distance a photon travels between collisions). Right next to the source the relationship is reversed: a denser rock scatters more photons to the detector. The rock is uniform in the volume of investigation, with no mudcake or standoff.
Holds when: Calibration on limestone saturated with fresh water: for calcite and for water ρa coincides with the true density. The constants 1.0704 and 0.1883 come from requiring ρa=ρb for these two substances. So for other substances the equality is approximate.
Holds when: The gas is methane (1.334≈1.0704⋅1.2467). For a heavier gas the coefficients are slightly different. The calibration formula is linear, and the fractions of the components add up to one. So the apparent densities of the components mix by volume in the same way as the true ones. The apparent density of gas is a calculated quantity for such mixing, not a physical density. It can even turn out negative.
Largest fraction of energy that a neutron loses in a single collision
EΔEmax=(A+1)24A
energies in the same units; the fraction is dimensionless
Holds when: An elastic collision (the colliding bodies keep their total energy of motion) with a nucleus at rest. In a glancing collision the loss is smaller, on average half the largest.
Holds when: Elastic scattering that is isotropic (equal in all directions) in the centre-of-mass frame. This is a good approximation for light nuclei at the energies of logging sources. For A=1 the second term vanishes and ξ=1. A medium made of nuclei of one kind. In a rock, the contributions of the nuclei add up according to their number and the probability of scattering on each.
Neutron porosity of a clean formation (simplified model)
ϕN≈ϕ[SxoHImf+(1−Sxo)HIh]
fraction
Holds when: A clean formation without shale. The curve is converted to the matrix of the formation and corrected for the borehole. The excavation effect is ignored: in a gas-bearing formation a real tool shows an even lower porosity (advanced section). Neutron absorbers have no influence (or the tool is epithermal).
Quick estimate of the porosity of a gas-bearing formation
ϕ≈2ϕN2+ϕD2
fraction
Holds when: An empirical rule (it comes from experience) for clean gas-bearing formations. Both curves are corrected for the borehole. For formations with liquid on limestone scales, we use the arithmetic mean (ϕN+ϕD)/2. It depends little on whether the formation is sandstone, limestone or dolomite. Charts and calculations give an accurate estimate, with the excavation effect and the gas saturation of the flushed zone (Gas, shale, carbonate and pressure effects on porosity logs).
Holds when: Hydrogen nuclei: other nuclei have a different gyromagnetic ratio and do not respond at the frequency of the tool. The design of the tool sets the field in the sensitive zone. In a non-uniform field, different distances from the tool correspond to different frequencies.
Holds when: One fluid with a single value of T1. In a rock, T1 has a spread of values, and we apply the formula to each part. There is no magnetisation before polarisation starts: the rock has just entered the field of the tool, or the previous cycle has completely “erased” it.
Holds when: The fluids are fully polarised and have a hydrogen index of 1. Otherwise the amplitudes are too low (see the section on hydrocarbons). Components with T2 comparable to TE or shorter die away before the first echoes, so they are poorly determined. The amplitudes of the late echoes are small compared with noise. So for a reliable decomposition, we stack the signal over several cycles (add the cycles together).
SI; it is convenient to take ρ2 in µm/s and S/V in µm⁻¹, and their product comes out in s⁻¹
Holds when: Fast diffusion: a molecule travels around the pore many times during the relaxation time, ρ2r/D≪1. The pore then relaxes as a whole. The surface term refers to the wetting phase. A non-wetting phase (oil in a water-wet rock, or gas) does not touch the walls and relaxes only by bulk and diffusion relaxation. The diffusion term is written for unrestricted diffusion in a uniform gradient. In small pores the walls limit the movement, and the term is smaller. Neighbouring pores do not exchange molecules during the measurement. Otherwise their times average out.
Holds when: The sum is taken over an agreed part of the distribution. In this lesson it is the effective porosity, T2>3 ms. Use the same convention in the calculation and in the calibration on core.
fraction of bulk volume (or per cent, that is, porosity units)
Holds when: The hydrogen index of all fluids is about 1, and polarisation is complete. Otherwise parts of the distribution are too low. The echo spacing is short enough to see clay-bound water. At a TE of about 1 ms the tool gives only effective porosity. The cutoffs are chosen for this rock (typical values are in the table below).
empirical: k in mD, ϕe in per cent; with another way of writing porosity, C changes
Holds when: This is an empirical relationship proposed for sandstones. We choose C (one often starts from a value of 10) and the exponents from the core of the formation. The cutoff separating BVI and FFI is calibrated for this rock. BVI is water only, and FFI is all the mobile fluids, with corrections for the hydrogen index and incomplete polarisation.
empirical: k in mD, ϕe as a fraction, T2LM in ms
Holds when: This is an empirical relationship. We choose a and the exponents from the core of the formation. Without core, one often starts from a of about 4 mD/ms² for sandstones. In carbonates the coefficients are different (a is much smaller). The pores are filled with water, and the distribution reflects pore sizes, not fluid properties. T2LM is calculated over the same part of the distribution as in the calibration.
NMR porosity with hydrocarbons in the flushed zone
ϕNMR=ϕ[Sxo+(1−Sxo)HIh(1−e−TW/T1h)]
porosity and saturation as fractions; TW and T1h in the same units
Holds when: Water-based filtrate with a hydrogen index of 1, fully polarised. The hydrocarbon signal is not lost because T2 is too short (for heavy oil this does not hold). The hydrogen index of water is taken as 1; for salt water it is slightly lower.
Holds when: A clean formation with a known matrix density; both tools see the same flushed zone. Water is fully polarised, and the hydrogen index of the filtrate is 1. The properties of the gas (ρg, HIg, T1g) are estimated for reservoir pressure and temperature.
Holds when: A tool with one working frequency. A multi-frequency tool, which measures several shells in turn, gathers data faster by about as many times as it has frequencies. The cycle time is polarisation plus the echo train. Service pauses are ignored. The distance the tool moves during stacking is simply added to the antenna length. This is a rough estimate of the vertical resolution.
Holds when: The noise is random and the same for all echoes. The time grid is set in advance and covers the whole possible range, from fractions of a millisecond to several seconds. A penalty in the form of a sum of squares is one common option. Penalties on the curvature of the distribution are also used.
SI. In logging, conductivity is usually given in mS/m: C=1000/R with R in Ω·m
Holds when: The sample is uniform, and the current flows along it evenly. Resistivity and conductivity refer to one temperature: in solutions and rocks they change a lot when the temperature rises.
Holds when: The electrodes are points, and the return electrode B is infinitely far away. The constant K is derived for a uniform medium. A borehole is not uniform, so here the formula only defines what we call apparent resistivity. By the reciprocity principle, the current and measuring electrodes can swap places, and the reading does not change.
Holds when: A uniform non-magnetic medium. It is an estimate for a plane wave. The field of a real sonde is more complicated, but the order of magnitude is the same. Displacement currents (the part of the current that comes from a changing electric field) are neglected; at frequencies of tens of kilohertz this is justified.
Readings of sondes with different depths of investigation in the step-invasion model
Rk=Jk(di)Rxo+[1−Jk(di)]Rt,k=1,2,3
any consistent units of resistivity, the same for every term
Holds when: The profile is a step: a flushed zone with Rxo out to the diameter di, and beyond it the formation with Rt. There is no transition zone and no annulus. This form is for the laterolog, where resistances add up. For induction logging the conductivities 1/R are added in the same way, with geometric factors. The bed is thicker than the vertical resolution of the sondes. Corrections for the borehole and the enclosing rocks have already been made.
Potential of a point current source in a uniform medium
U(r)=∫r∞4πr′2RIdr′=4πrRI
SI
Holds when: The medium is uniform, isotropic (the same in all directions) and infinite. The electrode is a point, and the return electrode is infinitely far away.
Reading of a focused microresistivity sonde in the step model
Rpad=JmcRmc+JxoRxo+(1−Jmc−Jxo)Rt
any identical units of resistivity
Holds when: A focused sonde: the measuring current crosses the zones one after the other, so the resistivities add up. The model does not suit the unfocused sondes of the microlog. A step profile: mudcake of thickness hmc, the flushed zone up to the invasion front, and beyond it the uninvaded formation. Jmc grows with the thickness of the mudcake, and the share of the formation grows as the front comes closer to the wall. The factors depend on the design of the tool and, to some extent, on the resistivity contrast, so real corrections come from the charts of the tool maker.
Holds when: A clean formation (one without shale). In a shaly one we use models of shaly reservoirs with the same substitutions (Water saturation in shaly sands). The water in the flushed zone is practically only filtrate. If the filtrate is mixed with formation water, the resistivity of that water is lower than Rmf, and a calculation with Rmf gives Sxo too high. A water-based mud: the filtrate conducts current. Rmf is brought to formation temperature.
Holds when:Sw and Sxo are calculated with the same porosity and the same coefficients. The filtrate displaced the formation fluids as in waterflooding (water pushing oil out of the rock). The hydrocarbons it moved are counted as moveable. Before the invasion, the flushed zone was the same as the uninvaded formation: the same bed and the same saturation.
Formation resistivity factor from the resistivity of the flushed zone
F=RmfRxo=ϕma
dimensionless
Holds when: A clean water-bearing bed completely flushed by the filtrate: Sxo=1, and the water in the flushed zone is filtrate. The invasion front is deeper than the depth of investigation of the microresistivity sonde, and the mudcake correction has been applied. In shaly beds, the conductivity of the shales makes F too low.
the interval transit times in the same units, µs/m or µs/ft; porosity as a fraction
Holds when: A clean rock without shale, with a known mineral composition. The rock is compacted and cemented. In loose sands the equation gives porosity that is too high. The pores are small and evenly spread, either intergranular or intercrystalline (between grains or between crystals). Vugs and fractures hardly enter the result. The pores in the zone of investigation are filled with liquid: the mud filtrate. The equation is empirical. It is not derived from the mechanics of elastic waves. It describes experimental data for such rocks.
Correction for undercompaction in the time-average equation
ϕS=Δtf−ΔtmaΔt−Δtma⋅Cp1,Cp≈328Δtsh
Δtsh in µs/m; if Δtsh is in µs/ft, the divisor is 100
Holds when: Estimating Cp from the neighbouring shales is a rough rule of thumb. Shales compact together with sands. A Δtsh above 328 µs/m (100 µs/ft) is taken as a sign of undercompaction. At smaller values we apply no correction (Cp=1). It is more reliable to find Cp in a reference bed where the porosity is known independently (see the calibration section below). The bed is clean, and its pores are filled with liquid.
Empirical transform of interval transit time to porosity
ϕS=CΔtΔt−Δtma
the interval transit times in the same units; porosity as a fraction
Holds when: A clean rock, with the pores filled with liquid. The formula approximates an experimental relationship in the range of ordinary porosity, up to about 35–37%. In looser rocks the grains no longer make a rigid framework, and the relationship is different. For gas-bearing beds we reduce the constant (see the section on oil and gas).
Effective porosity of a shaly bed from the sonic log
ϕe=ϕS−VshϕS,sh,ϕS,sh=Δtf−ΔtmaΔtsh−Δtma
the interval transit times in the same units; porosity and shale volume as fractions
Holds when: Laminated shale: the wave crosses the sand and shale layers one after the other, and the travel times add up. For thin layers this is an approximation. The shale layers are the same as the neighbouring shales from which Δtsh was taken. The shale volume Vsh was found independently (Determining shale volume). The rock is compacted, with liquid in the pores. Shale material dispersed in the pores of a compacted sandstone changes Δt little. For it the correction is too large.
Reducing correction of sonic porosity for oil and gas
ϕ=khcϕS
fraction
Holds when: A rough rule of thumb for high-porosity sands with shallow invasion. The multiplier depends on the saturation of the zone of investigation and on the stiffness of the rock, so it cannot be exact. We apply the correction only when the overestimate is obvious, for example by comparison with the water-bearing part of the same bed or with density porosity. In the empirical transform, for gas-bearing beds we take a smaller constant, C=0.6, instead of a multiplier.
Estimating secondary porosity by comparing methods
ϕ2=ϕt−ϕS
fraction
Holds when: A clean rock of known mineral composition, with liquid in the pores. Gas and shales distort the three methods in different ways. The parameters for calculating ϕS have been checked in an interval without vugs. There ϕS must agree with ϕt. Otherwise the error of the equation itself enters the difference. Vugs are not quite invisible to the wave, and they increase Δt a little. So the difference is usually somewhat smaller than their true volume. We exclude from the calculation the points with cycle skipping, which fractures often cause. This is an index, not a measurement: the errors of the two methods add up.
Holds when: An empirical relationship for clean water-bearing rocks. The porosity is less than about 37%. In looser rocks the grains no longer make a rigid framework, and the relationship is different.
Bulk density of a rock as the sum of the contributions of its components
ρb=i∑Viρi,i∑Vi=1
g/cm³ or kg/m³: any density units, the same for all terms
Holds when: All components are included: the volume fractions add up to one. The rock is uniform within the volume the tool sees: 10–15 cm from the borehole wall and 15–45 cm along the borehole. The reading is corrected for the mudcake (the layer of mud solids on the borehole wall) and the borehole. For ordinary minerals and liquids the tool reads practically the true density; for gas, salt and gypsum it does not (why).
Holds when: The fractions are volume fractions, as thin sections give them (thin slices of rock studied under a microscope). Mass fractions from X-ray analysis must first be divided by the mineral densities and then normalised (scaled so that they add up to one). The composition is constant throughout the formation for which the value is adopted.
Filtrate density at reservoir conditions (estimate)
ρmf≈(1+0.73P)kT,p
g/cm³
Holds when: The factor 1+0.73P is an approximation for NaCl solutions at 24 °C and atmospheric pressure. We assume that brine expands when heated and compresses under pressure roughly like fresh water. The uncertainty of the estimate is up to 0.01 g/cm³. For filtrates of complex composition and for oil-based muds, the density is measured or taken from the mud data.
Fluid density of the flushed zone with residual hydrocarbons
ρf=Sxoρmf+(1−Sxo)ρh
g/cm³
Holds when: The zone of investigation contains only filtrate and residual hydrocarbons; no formation water is left, or it cannot be told from the filtrate by density. For gas, use the apparent density (the value the density tool shows, which differs from the true density of gas). For methane, ρa,g≈1.334ρg−0.188 (derivation).
Holds when: The errors are small compared with the difference ρma−ρf (linear approximation). Independent random errors add in quadrature (as the square root of the sum of their squares). Systematic errors add with their signs.
Effective porosity of a shaly formation from density
ϕe=ϕD−VshϕD,sh,ϕD,sh=ρma−ρfρma−ρsh
densities in g/cm³; porosity and shale volume as fractions
Holds when: The shale in the formation is the same as in the reference shale unit; for dispersed shale this holds less well. The shale volume Vsh has been found independently (Determining shale volume). There is no gas in the zone of investigation.
Density of the sand laminae of a laminated formation
ρb,sd=1−Vshρb−Vshρsh
g/cm³
Holds when: The shale is laminated: the laminae are thinner than the vertical resolution of the tool (the thinnest bed it can separate), so the tool averages them over the volume. The shale laminae have the same density as the adjacent thick shales.
Holds when: The formation is clean; the flushed zone contains only filtrate and residual gas. Archie's equation is applied to the flushed zone with a=1 and m=n=2. With other coefficients the filtrate fraction itself depends on porosity, and the problem has to be solved by successive approximations. Rxo is measured with a microresistivity tool and corrected for mudcake. For light oil, use its density in place of ρa,g.
Neutron porosity as the sum of the contributions of fluid and shale
ϕN≈ϕHIfl+VshϕN,sh
fraction
Holds when: The contributions add up in proportion to the volumes. This is an approximation that suits fluid and shale. The matrix itself contributes nothing when the reading is in units of its own matrix. If the matrix is different, first convert the reading using the tool's conversion curve. Gas lowers the reading further because of the excavation effect (see the section on gas). There are no hydrogen-bearing minerals, such as gypsum or coal, in the formation.
Effective porosity of a shaly formation from the neutron log
ϕe≈ϕN−VshϕN,sh
fraction
Holds when: The pores are filled with water or oil (hydrogen index about 1); there is no gas. The shale material in the formation is the same as in the reference shale unit. The shale volume has been found independently, for example from the gamma ray log (GR) or SP (Determining shale volume).
Holds when: The formation is clean; the reading has been corrected for environmental conditions and converted to the formation matrix. The zone of investigation contains only filtrate and residual gas. The hydrogen index of methane is HIg=2.25ρg, with the gas density in g/cm³ (derivation).
Holds when: It approximates the results of neutron transport calculations. The numbers 2 and 0.04 are fitted to those results. Enter porosity and SwH as fractions. The formula describes the additional lowering. The lack of hydrogen itself is accounted for by the factor SwH.
Holds when: The first formula is derived in the topic on the neutron log, with atomic masses rounded to 12 and 1. If the composition is known, take x from it. For dry gas (almost pure methane), x=4 and HIh=2.25ρh. For oil, x≈2 and HIh≈1.29ρh. The estimate x≈4−2.5ρh is empirical (it comes from experience, not from theory). It reduces light hydrocarbons to methane, and oil with a density of about 0.8 g/cm³ to the formula CH2.
Apparent density of a hydrocarbon for the density tool
ρa,h=1.070412+1.01x2(6+x)ρh−0.1883
g/cm³
Holds when: The tool is calibrated on limestone saturated with fresh water: ρa=1.0704ρe−0.1883 (electron and apparent density). The fraction is the ratio 2∑Z/M for the group CHx. The group has 6+x electrons (six from carbon and x from hydrogen), and its molar mass is 12+1.01x. For methane the formula becomes ρa,g≈1.334ρg−0.188 (apparent gas density), and for oil with x=2 it becomes ρa,h≈1.222ρh−0.188.
Density and neutron porosity of a formation with residual hydrocarbons
ϕD=ϕ+kDϕShr,ϕN=ϕ−kNϕShr−Δϕex
fractions; densities in g/cm³
Holds when: The formation is clean and its matrix is known. Both curves have been corrected for environmental conditions and brought to the formation matrix. Both tools see the same zone, the first 10–30 cm behind the wall, with filtrate and residual hydrocarbons. The hydrogen index of the flushed-zone fluid SwH=1−kNShr is the same as in the excavation formula (more detail), with HImf=1.
Porosity and residual gas saturation from density and neutron porosity
ϕ=kD+kNkNϕD+kDϕN′,ϕShr=kD+kNϕD−ϕN′
fractions
Holds when: The formation is clean and its lithology is known. Both porosities are calculated for one, correct matrix and corrected for environmental conditions. In a shaly formation, they are used after the corrections for shale volume (next section). kD and kN are calculated from the gas density at reservoir conditions. The excavation correction depends on the unknowns ϕ and Shr, so repeat the calculation with the new correction until the result stops changing. Two or three repetitions are usually enough. Both tools see the same zone with the same saturation.
Density and neutron porosity corrected for shale volume
ϕD,c=ϕD−VshϕD,sh,ϕN,c=ϕN−VshϕN,sh
fraction
Holds when: The shale material in the formation is the same as in the reference shale unit. The contributions add in proportion to the volumes, as in the topics on each method. The gas shifts and the excavation effect apply to the pores outside the shale, so the solution for gas is applied to the corrected values.
Holds when: The same shales are compared; the borehole in the shales is smooth and density is reliable. The exponential trend follows from the compaction law ϕ=ϕ0e−cσ′ at normal pressure, when the effective stress grows in proportion to depth. We build the trend from the normally compacted part of the section and extend it downwards. The error of the extrapolation grows with distance.
Reading of a clean formation whose matrix is a mixture of two minerals
L=(1−f)L1(ϕ)+fL2(ϕ)
units of the method
Holds when: The formation is clean, the pores are filled with liquid of known density, and environmental corrections have been applied. For density the formula is an exact mass balance. For neutron porosity and interval transit time it is an approximation. The contributions of minerals do not add up in a strictly straight-line way, and on real charts the mixture lines are slightly curved.
Holds when: This is not a statement that the formation is limestone. It is only a scale common to all the lines. Choose the chart by the filtrate: for fresh and salty muds the lines are different.
Shift of the point of a gas-bearing formation on the neutron–density crossplot
ΔϕD,ls=kD,lsϕShr,ΔϕN=−kNϕShr−Δϕex
fractions
Holds when: The coefficients are the same as in the topic on gas, shale and carbonates, but density porosity is calculated for limestone. We estimate the gas volume ϕShr from the saturation of the flushed zone (found with microresistivity tools) and a first estimate of porosity, and then refine it. For oil, kD,ls is about 0.15 and kN is close to zero: the point shifts one or two units upwards.
Pe in barn/electron, U in barn/cm³, ρe in g/cm³ (as a measure of the number of electrons per unit volume)
Holds when: The components are mixed uniformly in the small volume behind the borehole wall that the tool sees. We take the electron density of each component from the calibration of the tool (limestone with fresh water), so the volumes mix the way the tool sees them. There is no barite or heavy minerals in the mudcake and the mud.
densities in g/cm³, U in barn/cm³, Pe in barn/electron
Holds when: The formation is clean, both minerals are known, and the pores are filled with liquid of known density and composition. The readings have been corrected for environmental conditions, and there is no barite in the mudcake or the mud. With gas in the flushed zone, porosity will come out too high (see the example with a gas-bearing limestone below).
Holds when: The formation is clean, and the pores are filled with liquid of known density. The apparent porosity must be correct: an error of 0.01 shifts ρmaa by about 0.02 g/cm³ (see below).
Holds when: The readings have been corrected for the borehole; there is no potassium chloride in the mud. The components have a constant composition. Natural minerals differ in contents several times over.
Organic carbon content from uranium (local calibration)
TOC≈aU+b
TOC in mass %, uranium in g/t
Holds when: Uranium is associated with organic matter, not with phosphates, fractures or secondary deposits. The calibration is valid for one unit and one depositional environment (the setting where the sediment settled). Do not transfer it to other units. Average the uranium curve over the interval, because individual readings are noisy.
Δt in µs/m, density in g/cm³; M is a conventional quantity (µs/ft per 1 g/cm³, divided by 100)
Holds when: A clean bed with liquid in the pores, and the readings corrected for environmental conditions. The interval transit time is linear in volumes (the time-average equation). With the empirical transform, M changes slightly with porosity. The multiplier is a conventional scale: 0.01 with Δt in µs/ft, that is, 0.01⋅0.3048=0.003048 with µs/m. The calculation and the chart must use the same multiplier.
Parameter N: the slope on the neutron–density plane
N=ρb−ρfϕNf−ϕN
neutron porosity as a fraction, density in g/cm³; N in cm³/g
Holds when: A clean bed with liquid. The neutron curve is in limestone units of the same tool for which the mineral points were built. N does not depend on porosity only if the neutron reading is linear in volumes. For sandstone and dolomite it is not linear, and their N changes by 0.02–0.04 between 0 and 30% pores.
Secondary porosity from the rise of the point on the M–N plot
ϕ2=(1−ϕt)(MmaM−1)
fraction
Holds when: A clean bed with liquid, made of one known mineral. The sonic log sees only the intergranular porosity ϕ1 and follows the time-average equation. The density and neutron logs see the total porosity ϕt=ϕ1+ϕ2.
Holds when: This form is for the time-average equation. The empirical transform has its own form (below). The mineral points on the chart were built with the same transform as the calculation.
wi and fi are fractions; M and N are in their own units (see above)
Holds when: The responses of all three methods are linear in volumes. For sandstone and dolomite, which have a non-linear neutron response, we take the mineral points at the porosity of the bed.
Holds when: A clean rock: only the water in the pores conducts current, and the surface conductivity of the grains and of the shale is negligible. Intergranular porosity, more or less uniform through the volume. a and m are fitted for the given rock, and they are used in the range of porosity where they were found. R0 and Rw refer to one temperature. For a bed, this is the formation temperature.
Holds when: A clean rock, and only water conducts. The water forms a connected system: films and filled small pores (a water-wet rock). Rt and R0 refer to one rock, one water and one temperature.
dimensionless; all resistivities in the same units
Holds when: A clean (shale-free) bed with intergranular porosity. Only the formation water conducts current. Rt is the resistivity of the uninvaded part of the formation, with corrections for the borehole, the enclosing rocks and invasion (Borehole environment: invaded zone and mudcake). Rw is at formation temperature. a, m and n suit the rock (from core or from a checked analogy). The rock is water-wet.
Holds when: The same as for Sw. The water in the flushed zone is practically only filtrate. A water-based mud (with an oil-based mud the filtrate does not conduct current).
Holds when: The errors are small compared with the quantities themselves (a linear approximation). The errors are independent. When we count them together, we add them with their signs, or, for random errors, in quadrature (as the square root of the sum of squares).
Estimating the filtrate and mudcake resistivities from the mud resistivity
Rmf≈0.75Rm,Rmc≈1.5Rm
Ω·m; all three quantities at the same temperature
Holds when: A statistical relationship for water-based muds in which NaCl is the main salt. It is unreliable for gypsum muds, for muds treated heavily with chemicals, and for heavily weighted muds (muds made dense with added solids). The scatter is large: tens of per cent. The heavier the mud, the smaller the ratio Rmf/Rm usually is. Correlations that take the mud density into account exist. If a sample has been measured, use the measurement.
Converting the equivalent resistivity of a brine to the true one (a formula fitted to the chart)
Rw,24≈146−337Rwe,2477Rwe,24+5
empirical: both resistivities in Ω·m at 24 °C (75 °F)
Holds when: The water is close to a NaCl solution. The formula applies when Rwe,24 is not greater than 0.12 Ω·m, that is, when the salinity is above about 50 g/kg. The procedure is: bring Rwe from formation temperature to 24 °C, convert it, and bring Rw back to formation temperature. The real chart for different temperatures is more complicated. At formation temperatures up to 150–200 °C the formula differs from it by up to about 10%, and at higher temperatures by more. The numbers 77, 5, 146 and 337 are fitting coefficients. They were chosen to make the curve fit the chart, and they have no physical meaning. When Rwe,24 is greater than 0.12 Ω·m, take Rw≈Rwe for NaCl-type waters, to within about 10%. For fresh waters of complex composition the relationship is different (see the advanced section).
Holds when: A clean bed with intergranular porosity (pores between the grains). Shale and conductive minerals distort the result. Rt has been corrected for the borehole, the enclosing rocks and invasion (Resistivity logs: electrical and induction tools). With shallow invasion, take the reading of the deep sonde. Use the same porosity and the same a and m that will later go into the saturation calculation. The rock type (lithology) and the water are the same in all the beds that we compare.
Water saturation from the apparent water resistivity
Sw=(RwaRw)1/n
fraction; the resistivities in the same units
Holds when: The same as for Rwa. We find the exponent n from core; without core data we usually assume 2. Rw and Rwa refer to the same temperature. If the depths differ a lot, convert Rw to the temperature of the bed.
Water saturation with a wrong cementation exponent, calibrated on a water-bearing bed
Sw,calc=Sw(ϕ2ϕ1)(m−m0)/n
dimensionless
Holds when: Both beds are clean, with the same water and the same relationship between the formation resistivity factor and porosity. The coefficient a is the same in both calculations, and we have taken the exponent n correctly. We take as Rw the Rwa of the reference bed, calculated with the same m.
Holds when: Clean beds of one rock type (lithology) with the same water. Rt and porosity are corrected in the same way as for the calculation of Rwa.
Converting the equivalent resistivity to the true one for fresh formation waters (a formula fitted to the chart)
Rw,24≈100.69Rwe,24−0.24−0.58
empirical: both resistivities in Ω·m at 24 °C (75 °F)
Holds when:Rwe,24 is greater than 0.12 Ω·m, and the water is an “average” fresh formation water with a noticeable share of divalent ions. The formula averages over waters of different composition, so for one particular water the error may be large. The numbers 0.69, 0.24 and 0.58 are fitting coefficients. The difference from the chart is about 10%.
Holds when: Clean beds of one lithology (a and m are the same) with the same formation water at a close temperature. Rt is corrected for the borehole, the enclosing rocks and invasion, and the porosity is corrected for environmental conditions and lithology. The exponent n is the same for all the beds of the interval.
Holds when: The exponent m is chosen in advance: the horizontal scale is built for it. The matrix, the fluid and the water are constant through the interval. The reading of the porosity method is linear in porosity: exactly for the density log, for the sonic log with the time-average equation, and approximately for the neutron log. The beds are clean, and Rt is corrected.
Ratio of the resistivities of the flushed zone and the uninvaded formation
RtRxo=RwRmf(SxoSw)n
any identical units of resistivity (usually Ω·m); Rmf and Rw at the same temperature
Holds when: A clean bed with intergranular porosity: only the water in the pores conducts current. Archie's equation with the same formation resistivity factor and the same exponent n in both zones. The water of the flushed zone is filtrate, and that of the uninvaded formation is formation water. The mud is water-based. Rxo and Rt are the resistivities of the zones after corrections for the borehole, the mudcake, the enclosing rocks and invasion, not the readings of the sondes.
Water saturation from the resistivity ratio, with the empirical relation between the saturations of the zones
Sw=Q5/8=(Rmf/RwRxo/Rt)5/8
dimensionless
Holds when: The empirical relation Sxo=Sw1/5: an oil-bearing bed, moderate invasion, an ordinary residual oil saturation. It has not been checked for gas, heavy oil, carbonates with vugs and oil-wet rocks. The saturation exponent is n=2. At another n the power is 5/(4n) (advanced section). All the assumptions of the ratio of the resistivities of the zones.
Water saturation from the resistivity ratio, with a known residual saturation of the flushed zone
Sw=(1−Shr)Q1/n
dimensionless
Holds when:Shr belongs to this bed and to the same process: displacement by filtrate during drilling. The value from waterflooding experiments can differ. All the assumptions of the ratio of the resistivities of the zones.
The curve from the resistivity ratio and its separation from the SP
ER=−KlgRtRxo,ΔE=ER−ESP=−nKlgSxoSw
mV; T in °C
Holds when: Clean beds with the same formation water. The SP is not distorted by the bed thickness and by streaming potentials (voltages caused by the flow of filtrate through the rock), or it is distorted in the same way in the beds that we compare. We shift the curve ER so that in a bed that is certainly water-bearing it coincides with the SP. The shift absorbs the difference between the equivalent and the true water resistivities, and constant errors of the ratio. The mud is fresh (Rmf>Rw), and Rxo/Rt was found with a correction for invasion. The second formula follows from the ratio of the resistivities of the zones, if the SP of the bed equals the SP of the water-bearing bed.
Resistivity of the water of the flushed zone with incomplete flushing
Rz1=Rwz+Rmf1−z
Ω·m; z is a fraction
Holds when: The conductivity of the mixture equals the volume-weighted conductivity of the waters. This is an approximation for solutions of one salt. Real mixtures conduct a little better, and the effect is slightly stronger. The share z is the same through the whole flushed zone. Its values in the example are assumed.
Porosity from the saturation found by the resistivity ratio
Fc=RwSwRnRt,ϕc=(Fca)1/m
dimensionless; Rt and Rw in the same units
Holds when: A clean bed, and the coefficients a, m and n are known. The porosity ϕt that we compare with is independent: from the density or neutron log, or from core.
An empirical estimate of water saturation with an annulus
Sw≈SwA(SwRSwA)1/4
fraction
Holds when: An empirical correction without a rigorous derivation. We use it only when the disagreement is really caused by an annulus. SwA and SwR were found from the same corrected Rt and Rxo.
Water saturation with the empirical relation at any saturation exponent
Sw=Q4n5
dimensionless
Holds when: The empirical relation Sxo=Sw1/5 and all the assumptions of the resistivity ratio. The relation was published for ordinary sandstones, and at an unusual n its applicability is also in doubt.
Ratio of the resistivities of the zones from the readings of two sondes
RtRxo≈ARdRs−(A−1),A=Js1
dimensionless
Holds when: A step profile. The deep sonde hardly sees the flushed zone, and Js is the same for all the beds. In practice A is chosen so that the formula agrees with the invasion correction chart in the chosen range of diameters. This also absorbs the contribution of the flushed zone to the deep sonde.
any GR units, the same for all three readings (API, the standard gamma ray unit, or µR/h); the index has no units
Holds when: Only the shale material makes the bed radioactive. The shale material in the bed is the same as in the reference shale bed, and the clean part is the same as in the reference clean bed. The bed is thicker than the vertical resolution of the tool (the thinnest bed the tool can read correctly). We cut off values outside the baselines: we use 0 if IGR<0, and 1 if IGR>1.
Holds when: It is an empirical relationship: the coefficient C reflects the rocks it was fitted to, and we refine it with core. The index uses correctly chosen clean-formation and shale baselines. The radioactivity comes from the shale material.
Holds when: It is an empirical relationship fitted to Tertiary sand–shale deposits. It is a special case of a family of curves with the curvature parameter A=3. The general form is in the advanced section.
Holds when: It is an empirical relationship. Geometrically, it is an arc of a circle with its centre at the point (−0.7;1.7), passing through (0;0) and (1;1).
Holds when: The bed is water-bearing. Hydrocarbons also reduce the SP amplitude, and then the estimate is too high. The bed is thick enough, or the amplitude is corrected for its thickness and resistivity. The resistivities of the mud filtrate (the liquid that seeps from the drilling mud into the rock) and of the formation water differ clearly. The water in the bed is the same as in the reference clean bed. The mud conducts electricity: the SP is not recorded in oil-based muds.
Shale volume from the gap between neutron and density porosity
Vsh≈ϕN,sh−ϕD,shϕN−ϕD
fractions
Holds when: The pores of the bed are filled with liquid. Gas brings the curves together or swaps them, and then the estimate is too low. The matrix is chosen correctly. In dolomite and in rocks with heavy minerals, the gap appears even without shale. The shale in the bed is the same as in the reference bed. The readings have been corrected for the effect of the borehole.
Mixing by mass, and the general form of the Stieber relationship
Vsh=A−(A−1)IGRIGR,A=ρclρsh
dimensionless; the densities are in the same units
Holds when: The index equals the mass fraction of shale: the shale in the bed is the same as in the reference bed, and only shale is radioactive. The densities of the components are constant.
Larionov coefficient from the shale volume at an index of 0.5
C=2log2(V0.51−1)
dimensionless
Holds when: The shale volume from core is defined in the same way as from logs: the volume of shale material, not only of clay minerals. The index in the core intervals is calculated with the same baselines as in the whole zone. One point gives only a first approximation. With all the data, we fit C by the least-squares method (the method that makes the sum of the squared errors as small as possible).
Qv is in meq per cm³ of pores. The 100 converts the CEC from 100 g to 1 g of rock.
Holds when: The CEC is measured on a sample of the same rock. Chemical methods can do this (for example, the uptake of methylene blue, a dye). So can electrical measurements at several salinities (salt contents of the water), described in the advanced section. All the exchangeable cations are in the pore water, that is, they belong to the whole pore volume ϕt.
Holds when: The saturation exponent is n=2. For another n, we solve the equation numerically (by repeated trial calculations). α and γ do not depend on water saturation. At γ=0 we get Archie's equation.
Holds when: The clean sandstone has the same porosity throughout the interval, and the shale material in the bed is the same as in the neighbouring shales. Dispersed shale fills the pores between the grains without pushing the grains apart. At Vsh=ϕsd, the pores are completely filled. Structural shale replaces grains without touching the pores between the grains. Shale volume and total porosity are found independently, for example from the gamma ray log and the density log (Gas, shale, carbonate and pressure effects on porosity logs).
Holds when: All the porosity of the shales is filled with bound water: the shales have no effective porosity. The shale material in the bed is the same as in the neighbouring shales. Vsh is the volume of shale material together with its water, as in the topic Determining shale volume.
Holds when: The layers are parallel to the bedding and thinner than the vertical resolution of the sonde (the measuring part of the logging tool), and the current flows along them. In a vertical borehole, an induction sonde (a resistivity tool that works with magnetic fields) measures this resistivity Rh (more). The sandy layers are clean and obey Archie's equation, and the shale layers are the same as the neighbouring shales. In a deviated borehole, or with dipping beds, the sonde sees a mixture of resistivities along and across the layers, and the model must be extended.
Holds when: It is an empirical equation: the shale term is proportional to shale volume and water saturation. For n=2, it is a quadratic equation of the general form with α=ϕem/(aRw) and γ=Vsh/Rsh. We take Rsh from the neighbouring shales. In practice people often adjust it so that water-bearing shaly beds give Sw≈1. Then it is a parameter of the model, not the resistivity of the shales.
Holds when: The same as for the Simandoux equation. The factor 1/(1−Vsh) comes from the laminated-sand model at m=2: the effective porosity is concentrated in the share 1−Vsh of the rock volume. In the version with m=n=2 it is also called the total shale equation: it does not assume a particular form of shale distribution.
Holds when: It is an empirical equation fitted to sections with fresh and slightly salty formation waters and high shale volume. The exponent 1−Vsh/2 changes from 1 at low shale volume to 0.5 in pure shale. It can be solved directly for Sw: Sw=[Rt−1/2/(…)]2/n.
Holds when: The exchangeable cations carry current in parallel with the pore water. When hydrocarbons push the water out, the cations stay in the water that remains, and their concentration grows as Qv/Swt. m∗ and n∗ are measured on core with a correction for shale (advanced section). They are larger than the “apparent” Archie exponents of the same core. B grows with temperature and falls in fresh water. At 25 °C in salty water it is about 4 (S/m)/(meq/cm³), and at 80 °C about three times larger. For the formation, we take it from published correlations for the formation temperature, not from the value at laboratory temperature. Qv is taken from the CEC of core, or from log estimates calibrated to core.
Holds when: There are two kinds of water in the pores: clay-bound water (volume Swb, conductivity Cwb) and free formation water. Hydrocarbons displace only the free water, so Swt≥Swb. From logs, we take Swb and Cwb from the neighbouring shales: Swb=Vshϕt,sh/ϕt, Rwb=1/Cwb=Rshϕt,shm/a. In very fresh water the layer of bound water becomes thicker, so we increase its share. In the examples of this lesson the water is salty enough to ignore this.
Conductivity of a water-saturated shaly sample against the brine conductivity
C0=F∗1Cw+F∗BQv
conductivities in S/m at the measurement temperature
Holds when: The sample is fully saturated with brine, Swt=1. B is constant. This holds for brines that are salty enough. In fresh brines B decreases, and the points fall below the line, so we find the slope from the salty brines.
Holds when: The shale is spread evenly in the pores and conducts current together with the water, like a mixture of electrolytes. The cementation and saturation exponents equal 2. We estimate ϕim from the sonic log, which records the speed of sound in the rock (ϕim≈ϕS). For the sonic log, the dispersed shale “looks like” pore water. We estimate q from the gap between density and sonic porosity. If the density of the shale is close to the density of the matrix, the density porosity is close to the effective porosity ϕim(1−q), and so q≈(ϕS−ϕD)/ϕS.
Holds when: The resistivity of the dispersed shale is several times larger than Rw, and the bed is not too shaly. The other assumptions are the same as for the full model.
Induction sonde reading with the mud in the borehole
Ra1=RmGm+Rt1−Gm
Ω·m; Gm has no units
Holds when: This is the geometric-factor theory: the skin effect (a weakening of the signal in conductive rock) is negligible or already corrected. The bed is thicker than the vertical resolution of the sonde (the thinnest bed it can read correctly), and there is no invasion. These are handled by the later corrections. Gm belongs to a particular tool, borehole diameter and tool position. The values in the example are assumed.
Holds when: A step profile. The bed is thicker than the vertical resolution of the sondes. The corrections for the borehole and the enclosing rocks have been applied. Laterologs: the zones add up in series, like resistances. For induction sondes the conductivities add up, and the chart is different. The microresistivity tool shows Rxo. The factors J do not depend on the resistivity contrast (for real sondes they depend on it only weakly).
Holds when: The scales are parallel, and heights are measured from a common horizontal line. f1, f2 and f3 are any monotonic functions (functions that only rise or only fall).
gaps between the tool and the wall at the transmitter and at the receiver, mSonic log
aw,amf
chemical activities of the formation water and the filtrate at formation temperature, the “effective” ion concentrations, —Spontaneous potential (SP) log
interval transit time of the fluid in the pores of the zone of investigation (the part of the rock that the tool reads), which is the mud filtrate, µs/mPorosity from the sonic log
interval transit time of the fluid in the zone of investigation (the part of the rock that the tool measures): 620 µs/m (189 µs/ft) for fresh filtrate (the liquid that enters the rock from the drilling mud), about 607 µs/m (185 µs/ft) for salty filtrate, µs/mM–N and MID plots for complex lithologies
formation resistivity factor (the ratio of the resistivity of the fully water-saturated rock to that of the water), cementation exponent and saturation exponent, found with allowance for shale and based on total porosity, —Water saturation in shaly sands
hydrogen index of the fluid in the zone of investigation (its hydrogen content relative to water): about 1 for water and oil, less for gas, —Porosity from the neutron log
pseudo-geometric factors of the mudcake and the flushed zone: their shares in the reading. The remainder 1−Jmc−Jxo belongs to the uninvaded part of the formation, —Microresistivity tools and flushed-zone resistivity
ratio of the density of water at reservoir temperature and pressure to its density at the surface: about 1.00 at 50 °C and 15 MPa, 0.97 at 100 °C and 30 MPa, 0.94 at 150 °C and 45 MPa, —Porosity from the density log
readings of formations made of pure mineral 1 and pure mineral 2 at the same porosity and with the same fluid, the sameLithology and porosity from crossplots
slope of the line: the drop in temperature per logarithmic cycle (one tenfold change of the ratio). It depends on the cooling power and on the thermal conductivity of the rocks, °CFormation water and formation temperature
interval transit time of the matrix, that is, of a tight rock with no pores (table below), µs/mPorosity from the sonic log
interval transit time of the matrix for this transform: 184 µs/m (56 µs/ft) for quartz sandstone, 161 (49) for limestone, 144 (44) for dolomite, µs/mPorosity from the sonic log
Δtma,Δtf
interval transit times of the matrix and of the fluid, the same as in the calculation of ϕS, µs/mPorosity from the sonic log
resistances of the sections of the current path: the mud in the borehole, the bed with its invaded zone, the enclosing shales (resistances of bodies, not resistivities), ΩSpontaneous potential (SP) log
equivalent filtrate resistivity at formation temperature; for ordinary fresh muds ≈0.85Rmf, Ω·mSpontaneous potential (SP) log
Rmfe,Rwe
equivalent resistivities of the filtrate and the formation water at formation temperature, which are inversely proportional to the activities, Ω·mSpontaneous potential (SP) log
resistivity of the uninvaded part of the bed: the reading of the deep sonde (a tool that looks far into the rock) after corrections, Ω·mDetermining formation water resistivity
the horizontal coordinate of the plot (the abscissa); at m=2 it is 1/Rt, that is, the square root of the conductivity, (Ω⋅m)−1/mResistivity–porosity crossplots
static amplitude: the deflection from the shale baseline opposite a clean water-bearing bed when no currents flow; negative when the formation water is more saline than the filtrate, mVSpontaneous potential (SP) log
transport numbers of the Cl⁻ anion and the Na⁺ cation in an NaCl solution, t−+t+=1; for a dilute solution at 25 °C t−≈0.60, —Spontaneous potential (SP) log
excavation-effect correction with the lithology coefficient K: 1 for sandstone, 1.046 for limestone, 1.173 for dolomite, fractionLithology and porosity from crossplots
θc
critical angle between the ray in the mud and the perpendicular to the wall, —Sonic log
dynamic Poisson's ratio (the ratio of the sideways shortening to the lengthening when the rock is stretched), —Sonic log
ξ
mean logarithmic energy loss per collision, —Neutron log
ρ
density of the solution; for NaCl solutions near 20 °C, approximately ρ≈1+0.73w, where w is the mass fraction of salt, g/cm³ (= kg/L)Formation water and formation temperature
effective porosity: the volume of the pores outside the shale material, as a fraction of the volume of the whole rock, fractionPorosity from the sonic log
effective porosity: the pore volume outside the shale, as a fraction of the volume of the whole rock, fractionPorosity from the density log
effective porosity: the pore volume outside the shale material, as a fraction of the rock volume, fractionPorosity from the neutron log
neutron porosity in limestone units (the tool is scaled to read the true porosity in limestone), after corrections for environmental conditions, fractionM–N and MID plots for complex lithologies
ϕN,ϕD
neutron and density porosity calculated for the same (correct) matrix, fractionNeutron log
neutron and density porosity of the bed, calculated for the same matrix, fractionDetermining shale volume