Array Propagation Resistivity: Definition, LWD, and Geosteering
Array propagation resistivity is a logging-while-drilling (LWD) measurement technique that determines formation resistivity by transmitting high-frequency electromagnetic waves (typically 400 kHz to 2 MHz) from one or more transmitter coils on a drill collar and measuring both the phase difference and the amplitude attenuation of the wave at multiple receiver coil pairs spaced at fixed distances — commonly 16, 22, 28, 34, and 40 inches — along the same collar, yielding two distinct resistivity curves (phase-shift resistivity and attenuation resistivity) per spacing for a total of up to ten simultaneous resistivity channels with different radial depths of investigation. Because the tool operates continuously while drilling, it delivers a real-time multi-depth radial resistivity profile through mud-pulse or electromagnetic telemetry to the surface, enabling the driller and geosteering geologist to make immediate trajectory decisions to keep a horizontal well within a target pay zone — a capability that no post-drill wireline tool can replicate. The Schlumberger Array Resistivity Compensated (ARC) tool, the Halliburton Electromagnetic Wave Resistivity (EWR-Phase 4) tool, and the Baker Hughes Multiple Propagation Resistivity (MPR) tool are the principal commercial platforms; all use the same fundamental physics but differ in transmitter-receiver geometry, frequency set, azimuthal sector capability, and telemetry bandwidth. In the Western Canada Sedimentary Basin, array propagation resistivity has become the standard real-time formation evaluation and geosteering sensor in Montney, Cardium, Viking, and Duvernay horizontal completions, where the ability to detect bed boundaries ahead of or beside the bit is the difference between a productive 3,500-metre lateral and a poorly-placed wellbore that penetrates water-bearing zones.
Key Takeaways
- Phase shift and attenuation measurements — two curves per spacing: A propagation resistivity tool transmits an electromagnetic wave at a fixed frequency (e.g. 2 MHz for the short spacings, 400 kHz for the long spacings) from a transmitter coil T1. The wave propagates outward through the formation and is detected at two receiver coils R1 and R2 spaced a fixed distance apart on the same collar. Two independent measurements are extracted from the R1/R2 signal pair: the phase difference (delta-phi, in degrees) between the wave at R1 and the wave at R2, and the attenuation ratio (in dB) between the amplitudes at R1 and R2. Each measurement is converted to an apparent resistivity through a lookup table or analytic inversion derived from Maxwell's equations for the measured frequency and spacing. The phase-shift resistivity (Rps) has shallower radial investigation (typically 20 to 40 per cent less depth than attenuation) because the phase is more sensitive to near-borehole conductivity changes. The attenuation resistivity (Ratt) penetrates deeper into the formation and is therefore less affected by mud invasion in permeable zones. Together, the pair behaves analogously to the shallow and deep curves on a wireline array induction tool and is interpreted using the same invasion logic: positive separation (Ratt greater than Rps) indicates resistive invasion into a water-wet pay zone; negative separation indicates conductive invasion.
- Depth of investigation and the role of frequency and spacing: Propagation resistivity depth of investigation (the D50 radius at which the integrated geometric factor equals 50 per cent of the total) scales approximately with the transmitter-receiver spacing and inversely with the square root of frequency. At 2 MHz with a 16-inch spacing, D50 is approximately 8 to 12 inches from the borehole wall — effectively a near-borehole invaded-zone measurement. At 400 kHz with a 40-inch spacing, D50 extends to approximately 30 to 40 inches, approaching the undisturbed formation in moderately invaded formations. Modern multi-frequency tools (ARC5 from Schlumberger uses both 400 kHz and 2 MHz at each spacing) acquire both frequency channels simultaneously, doubling the number of radial investigation depths from five to ten curves and providing significantly more constraint for the 1D radial inversion to determine Rt, Rxo, and di from real-time LWD data without any post-drill wireline run. In tight formations where invasion is minimal (permeability below 0.1 mD in Montney or Duvernay), the phase and attenuation curves at all spacings read essentially the same value, and the deepest attenuation resistivity at 40-inch spacing is used directly as Rt for on-site Archie water saturation calculations delivered to the operator within 10 to 30 minutes of drilling through the zone.
- Azimuthal propagation resistivity for geosteering and formation imaging: A standard propagation resistivity tool has coils wound symmetrically around the drill collar and responds to formation resistivity averaged around the full 360-degree azimuth. An azimuthal array propagation resistivity tool — such as the Schlumberger ARC with the azimuthal deep reading feature, or the Halliburton EcoScope — adds a rotational encoding capability that records individual measurements in 16, 32, or 64 azimuthal sectors as the drill string rotates at 60 to 120 rpm. When the tool approaches a formation boundary (e.g. the bottom of a Cardium sand as the horizontal well begins to exit into the underlying marine shale), the sector facing the approaching boundary reads lower resistivity than the sector facing away, creating a measurable up/down or left/right resistivity asymmetry that the geosteering geologist uses to detect the boundary 1 to 4 metres before the drill bit crosses it. This "look-around" capability allows the driller to nudge the trajectory up or down to stay within the pay interval, improving reservoir contact by 15 to 30 per cent compared to non-azimuthal tools that only detect the boundary when the bit has already exited the target zone. In the WCSB Cardium play at Pembina and Willesden Green, azimuthal propagation resistivity is standard practice for all horizontal completions, where the Cardium A and B sand horizons may be only 2 to 4 metres thick and require constant trajectory corrections to maintain zone penetration along a 4,000-metre lateral.
- Anisotropy response and dip effects in horizontal wells: Formation resistivity is intrinsically anisotropic in laminated reservoirs: horizontal resistivity Rh (measured when current flows parallel to bedding) is typically 2 to 10 times lower than vertical resistivity Rv (measured when current flows perpendicular to bedding) due to the alternating high-porosity/high-conductivity and low-porosity/low-conductivity laminae. A horizontal propagation resistivity tool drills parallel to bedding and therefore predominantly measures Rh at shallow depths of investigation and an average of Rh and Rv at deeper investigation depths — causing the deep attenuation curve to appear higher than the shallow phase curve even in uninvaded formations, a pattern that superficially mimics resistive invasion but is actually an anisotropy artefact. In Montney horizontal wells where Rv/Rh ratios of 3 to 8 are common in the inter-laminated siltstone/carbonate facies, the Ratt/Rps separation at a 40-inch spacing can reach 3 to 5 ohm-m even in zones with no mud invasion, and the geosteering geologist must account for the anisotropy contribution before interpreting the curve separation as invasion-related hydrocarbon. Triaxial induction tools and deep azimuthal resistivity tools (see deep azimuthal resistivity) are used when quantifying Rv and Rh separately is required for thin-bed log analysis or fracture characterisation.
- Integration with real-time geosteering workflow and telemetry constraints: Array propagation resistivity data are transmitted to surface via mud-pulse telemetry at 3 to 12 bits per second (bps) for conventional wells, or via wired drill pipe at 57,600 bps for premium data systems. The telemetry bandwidth constraint means that typically only four to six resistivity channels (one phase and one attenuation per selected spacing) plus azimuthal sector data are transmitted in real time; the full 10-channel raw dataset is stored in the tool's 512 MB to 2 GB onboard memory and is downloaded during bit trips or at total depth. The geosteering workflow uses the transmitted channels — typically the 16-inch phase (near-borehole, invasion diagnostic), the 28-inch attenuation (medium, pay indicator), and the 40-inch attenuation (deep, true resistivity proxy) — alongside the gamma ray and inclination survey to update the geological model and propose the next trajectory waypoint. In the Montney at Dawson Creek, operators running 5,500-metre multilateral pad wells update the geosteering model every 90 to 150 metres of drilled footage, and the array propagation resistivity provides the primary formation boundary detection horizon through which all trajectory decisions are made.
Physics of High-Frequency Propagation: Wave Attenuation and Phase Shift in Resistive Formations
When a time-varying electromagnetic wave at frequency f is launched from a transmitter coil embedded in a conductive medium (the formation), the wave propagates radially outward as an exponentially attenuating wave of the form E(r) = E0 exp(-r/delta) exp(i x omega x r/v), where delta is the skin depth (delta = sqrt(2/(omega x mu x sigma)) in SI units), omega = 2*pi*f, mu is magnetic permeability, sigma is formation conductivity, and v is phase velocity. The skin depth delta sets the characteristic decay length: for sigma = 0.1 S/m (10 ohm-m formation, common in a Viking sand pay zone) and f = 2 MHz, delta is approximately 1.1 metres (43 inches), comparable to the tool spacing. For sigma = 1 S/m (1 ohm-m, shale), delta drops to 0.36 metres (14 inches), meaning the wave is heavily attenuated within the shale layer and the 40-inch spacing receiver detects a much weaker signal — directly encoding the formation conductivity in the amplitude ratio (attenuation measurement). For sigma = 0.005 S/m (200 ohm-m, tight carbonate or anhydrite), delta is approximately 5 metres and the wave propagates with minimal attenuation, and the attenuation measurement is dominated by geometric spreading rather than formation conductivity, reducing sensitivity in very high-resistivity formations above roughly 500 ohm-m.
The phase shift between receiver R1 at distance d1 from the transmitter and receiver R2 at distance d2 is delta-phi = (omega/v) x (d2 - d1) + phi0, where v = 1/sqrt(mu x epsilon_eff) is the phase velocity in the medium and epsilon_eff is the effective permittivity. In resistive media (typical reservoir rock at 2 MHz), the phase velocity is controlled primarily by resistivity and the phase shift is a nearly linear function of formation conductivity for mid-range resistivities of 1 to 200 ohm-m. At very high resistivity (above 300 ohm-m), the formation begins behaving as a dielectric rather than a conductor, and the permittivity term dominates the phase shift, causing the phase-shift resistivity curve to read anomalously low in tight carbonates and cemented sandstones — a phenomenon called "dielectric effect" that can cause Rps to underestimate true resistivity by a factor of 2 to 5 in formations above 300 ohm-m at 2 MHz. Modern tools apply a dielectric correction using the simultaneous 400 kHz channel (less dielectric-affected) and a joint inversion to separate resistive and dielectric contributions.