band-limited function

Band-limited function in digital well log data processing refers to the spectral constraint inherent in wireline and logging-while-drilling formation evaluation measurements, where the physical response of each sensor tool and the digital sampling interval of the acquisition system together determine the maximum spatial frequency at which genuine formation features can be resolved in the recorded log curve, and where the Nyquist criterion establishes that the highest detectable frequency in a log digitized at a depth sampling interval of dz is 1/(2dz) cycles per metre, placing a hard upper limit on the bed thickness resolution achievable from any digitally recorded log regardless of the intrinsic vertical resolution of the sensor itself. In Western Canada Sedimentary Basin formation evaluation practice, wireline logs are acquired at a standard sampling interval of 0.1524 m (0.5 feet) in conventional wells and at 0.0381 m (0.125 feet) in high-resolution acquisition mode for thin-bed analysis in laminated Cretaceous reservoirs such as the Mannville channel sands and the McMurray Formation; at the standard 0.1524 m sampling interval, the Nyquist frequency is 3.28 cycles per metre, meaning that a bed thinner than 0.152 m (approximately 6 inches) cannot be resolved as a distinct formation feature in the log data even if the tool's vertical resolution (typically 0.3 to 0.6 m for standard gamma ray and resistivity tools) would theoretically have detected it with continuous recording. The band-limited nature of digital log data has three practical consequences in WCSB formation evaluation: first, thin-bed analysis of the McMurray oil sands, Cardium tight sand laminations, and Montney siltstone sequences requires resampling the log data at the acquisition sampling rate or finer before applying net-pay cutoffs, because downsampling (decimating to coarser depth intervals for display or correlation) applies a further low-pass filter that aliases the thin-bed character into a smooth curve that obscures hydrocarbon-bearing laminations; second, depth matching between logs acquired at different sampling intervals (e.g., high-resolution image logs at 0.01 m versus standard gamma ray at 0.1524 m) requires interpolation that introduces no new frequency content above the Nyquist limit of the coarser log, so the matched pair retains the spectral bandwidth of the lower-resolution measurement; third, spectral analysis of log sequences (Fourier transform of GR or resistivity over a depth interval) identifies cyclical depositional patterns (parasequences, metre-scale bedding cycles, lamination frequency) only up to the Nyquist frequency, so the sampling interval must be chosen to satisfy the Nyquist criterion for the smallest cycle the geologist intends to detect before the acquisition program is designed. Understanding how the band-limited nature of digital well logs constrains bed resolution in WCSB thin-bed reservoirs, how anti-aliasing filters are applied before downsampling, how Nyquist frequency limits the spatial frequency content of log cyclicity analysis, and how tool vertical resolution and digital sampling interval interact to set the effective band limit of the final log curve gives WCSB petrophysicists, formation evaluation engineers, and log analysts the spectral framework to correctly interpret thin-bed log signatures and to design acquisition programs with sufficient sampling rates for the geological targets of interest.

• Nyquist frequency and bed thickness resolution limits in WCSB standard wireline acquisition: A wireline log sampled at 0.1524 m depth intervals has a Nyquist frequency of 3.28 cycles/m, corresponding to a minimum resolvable bed thickness of 0.152 m (6 inches) under ideal conditions of perfect sensor response. The standard gamma ray tool in WCSB Mannville channel sand evaluation has a vertical resolution of 0.46 m (1.5 feet) due to detector geometry and logging speed, which sets an effective band limit of approximately 1.09 cycles/m regardless of the 3.28 cycles/m Nyquist limit imposed by sampling; the effective band limit of the log data is therefore the minimum of the Nyquist limit and the tool resolution limit. For WCSB McMurray oil sands thin-bed analysis, high-resolution GR tools with vertical resolution of 0.15 to 0.20 m are logged at 0.0381 m sampling intervals (Nyquist 13.1 cycles/m) to ensure the sampling rate does not limit the tool's intrinsic resolution; the effective band limit of 5 cycles/m is then set by the tool resolution of 0.15 to 0.20 m rather than the sampling interval.
• Anti-aliasing filtering before log data decimation for WCSB database loading: WCSB formation evaluation databases typically store log data at 0.1524 m sampling for standard well comparisons, requiring that logs acquired at finer intervals (0.0381 m high-resolution, 0.01 m image log depth tracks) be decimated before loading. Decimation without anti-aliasing filtering aliases the high-frequency content above the Nyquist frequency of the target sampling rate into the low-frequency portion of the decimated log, creating artificial cyclicity that misrepresents the formation. WCSB petrophysicists apply a Butterworth or Gaussian low-pass filter with a cutoff frequency of 2.5 cycles/m (80% of the 3.28 cycles/m Nyquist) before decimating from 0.0381 m to 0.1524 m, attenuating frequency content above the Nyquist limit before it can alias into the coarser sampled log. The phase response of the anti-aliasing filter must be linear (zero-phase or symmetric FIR) to preserve depth alignment between the filtered log and the original formation boundaries.
• Spectral analysis of WCSB GR logs for parasequence identification in Cretaceous stratigraphy: Fourier transform analysis of gamma ray logs over stratigraphic intervals in WCSB Cretaceous formations (Mannville, Colorado, Viking) identifies metre-scale depositional cycles by converting the depth-domain GR signal to a spatial frequency spectrum and identifying peaks that correspond to repeating fining-upward or coarsening-upward packages. A WCSB Mannville log over a 200 m interval sampled at 0.1524 m (1,312 depth samples) resolved at a fundamental spatial frequency of 0.005 cycles/m (200 m cycle) up to the Nyquist limit of 3.28 cycles/m; dominant spectral peaks at 0.2 to 0.4 cycles/m (2.5 to 5 m cycles) in the Clearwater interval correlated to ichnofacies-bounded parasequences confirmed by core description, while a secondary peak at 1.1 cycles/m (0.9 m cycles) corresponded to tidal lamination bundles identified in the core. The band-limited nature of the log data meant that lamination bundles thinner than 0.15 m were not represented in the spectral output, limiting the parasequence analysis to features detectable at the 0.1524 m sampling interval.
• Depth matching band-limited logs in WCSB LWD and wireline data integration: WCSB horizontal Montney and Duvernay wells typically acquire formation evaluation data from both LWD tools (logged while drilling at 0.1 to 0.3 m depth increments depending on ROP and acquisition rate) and post-drill wireline logs (acquired at 0.1524 m); integrating these two datasets requires depth matching that must respect the different band limits of each dataset. The LWD log acquired at 0.3 m effective depth resolution has a Nyquist frequency of 1.67 cycles/m, while the wireline log at 0.1524 m has a Nyquist of 3.28 cycles/m; depth matching by crosscorrelation of gamma ray curves must use only the frequency content common to both logs (0 to 1.67 cycles/m) to avoid the wireline GR's higher-frequency content from distorting the correlation result. WCSB formation evaluation software (Techlog, IP) applies automatic band-matching before LWD-wireline depth correlation, low-pass filtering the wireline GR to the LWD Nyquist frequency before computing the depth shift, preserving the spatial resolution of both datasets after registration.
• Upscaling band-limited log data for WCSB reservoir simulation models: Reservoir simulation models for WCSB Montney and Duvernay tight reservoirs use grid cells of 0.5 to 5 m vertical thickness, requiring that petrophysical properties derived from 0.1524 m wireline logs be upscaled (averaged) to the simulation grid cell scale. Upscaling is mathematically equivalent to applying a low-pass filter with a cutoff frequency equal to 0.5/(cell thickness) cycles/m; for a 1 m simulation cell, the cutoff is 0.5 cycles/m, removing all formation variation at scales finer than 1 m from the simulation model. WCSB petrophysicists use power-law averaging (arithmetic for horizontal permeability, harmonic for vertical) rather than simple arithmetic mean when upscaling from the 0.1524 m log scale to the simulation cell scale, because the heterogeneity within each cell is statistically predictable from the variance of the band-limited log data and must be preserved through the upscaling transformation to avoid underestimating effective permeability in laminated reservoirs.

Nyquist Aliasing Causing False Cyclicity in a WCSB Mannville Log Analysis

A WCSB formation evaluation team performing cyclicity analysis on Mannville GR logs from 15 wells exported the logs from their petrophysical software at a 0.3048 m (1-foot) sampling interval for loading into a cycle analysis spreadsheet, inadvertently applying a 2:1 decimation from the original 0.1524 m acquisition data without anti-aliasing filtering. The Nyquist frequency of the 0.3048 m sampled output was 1.64 cycles/m; GR frequency content between 1.64 and 3.28 cycles/m from the original data aliased into the 0.4 to 1.64 cycles/m band of the decimated log. Fourier analysis of the decimated logs showed a dominant spectral peak at 0.8 cycles/m (1.25 m cycles) across all 15 wells, which the team initially interpreted as a regionally persistent parasequence boundary. Verification by re-running the analysis on the original 0.1524 m data with proper anti-aliasing filtering revealed that the 0.8 cycles/m peak was an aliasing artifact of frequencies near 2.46 cycles/m (0.41 m cycles) in the original data, not a real depositional cycle. The false parasequence correlation was discarded, preventing a well placement decision that would have targeted the erroneous cycle boundary rather than the actual Mannville pay interval.

Related Terms

Band-limited function is the primary entry covering the seismic processing application where wavelet spectral bandwidth constrains seismic resolution; this companion entry covers the well log data processing application, where the Nyquist frequency of the digital sampling interval and tool vertical resolution together set the effective band limit that constrains bed thickness detection in WCSB formation evaluation. Wireline log is the formation evaluation measurement whose digital acquisition creates the band-limited dataset analyzed in this entry; standard WCSB acquisition at 0.1524 m sampling sets the Nyquist limit for gamma ray, resistivity, density, neutron, and sonic logs used in petrophysical interpretation and stratigraphic correlation. Gamma ray log is the most commonly Fourier-analyzed wireline log in WCSB Cretaceous stratigraphic studies; its band-limited spatial frequency content at the acquisition sampling rate determines which depositional cycle thicknesses can be resolved in parasequence and sequence stratigraphic analysis. Logging while drilling (LWD) produces formation evaluation logs with effective sampling intervals controlled by rate of penetration and acquisition rate rather than a fixed depth wheel, creating variable Nyquist limits that must be accounted for when depth-matching LWD and wireline datasets in WCSB horizontal well formation evaluation. Petrophysics is the formation evaluation discipline that applies band-limited log data to derive reservoir properties; upscaling petrophysical properties from the 0.1524 m wireline log scale to simulation grid cells is a low-pass filtering operation whose cutoff frequency is determined by the cell thickness of the WCSB reservoir model.