Even when geologically based methods are used to determine fault rock permeabilities and thicknesses for input into flow simulators, a wide range of simplifying assumptions regarding fault structure and content are still present. Many of these assumptions are addressed by defining quantitative and flexible methods for realistic parameterization of fault-related uncertainties, and by defining automated methods for including these effects routinely in full-field flow simulation modelling. The fault effects considered include: the two-phase properties of fault rocks; the spatial distributions of naturally variable or uncertain single-phase fault rock properties and fault throws; and the frequencies and properties of sub-resolution fault system or fault zone complexities, including sub-seismic faults, normal drag and damage zones, paired slip surfaces and fault relay zones. Innovative two-phase or geometrical upscaling approaches implemented in a reservoir simulator preprocessor provide transmissibility solutions incorporating the effect, but represented within the geometrical framework of the full-field flow simulation model. The solutions and flexible workflows are applied and discussed within the context of a sensitivity study carried out on two faulted versions of the same full-field flow simulation model. Significant influence of some of these generally neglected fault-related assumptions and uncertainties is revealed.

Faults are represented conventionally in production flow simulation models using transmissibility multipliers which capture the single phase, but not the two phase, fault-rock properties. Available data indicate that fault-rocks have similar two phase properties to sediments of the same permeability, hence existing methods can be applied to estimate two phase fault-rock properties from their intrinsic permeabilities. Two methods of representing the two phase fault-rock properties implicitly in the flow simulator are compared, using one-dimensional numerical flow models containing water-wet faults with imbibition capillary pressure curves. The method which is the closer two phase analogue of the single phase transmissibility multiplier is inappropriate, as the implementation is unreasonably unwieldy. A simpler implementation is to derive pseudo-relative permeability functions including the fault-rock properties in the upstream grid block; these properties are then incorporated directly in the simulator. Relative transmissibility multiplier functions can be back-calculated from the pseudo-relative permeability functions, and indicate how closely the single phase multiplier approximates two phase flow through the fault. Implementation in a 3D model with complex fault juxtapositions validates the approach, and a practical workflow for the routine inclusion of two phase fault-rock properties in conventional faulted flow simulation models is outlined.

No abstract

Fault representation and scaling in flow models are examined with respect to fault zone properties, the accuracy with which they can be determined, and how these variables and fault geometries can be incorporated realistically in to flow models. Outcrop data show that fault displacement/thickness ratios and permeability vary widely. For simple single fault models, results for numerical models are compared with analytical and statistical methods. Representation of a fault as a transmissibility surface conflates the effects of four variables -fault zone thickness and permeability, grid-block size and matrix (host-rock) cell permeability. Random spatial variation of transmissibility factor values is well represented by a uniform transmissibility factor which is the arithmetic mean of the values representing log-normally distributed permeability and thicknesses. Realistic ranges of fault zone thicknesses can be represented without grid-block refinement by an upscaling method based on simple transformation of transmissibility factor curves derived from a range of coarse grid-block models. Sub-seismic faults have significant effects on effective permeability of model volumes at kilometre scales only when the faults are assigned a permeability less than c. 0.001 of the matrix permeability.

S U M M A R YIt has become increasingly important to develop fast and accurate automatic procedures to process and fully exploit increasing large seismic data sets. Traditionally these data sets are processed manually, which requires significant amounts of both manpower and time with sometimes-variable results. We have developed a cost minimization approach to train three automatic pickers: an Sta/Lta, T pd and the PAI-K picker at each station within a dense temporary network located in northern Chile and southern Bolivia. The optimum picking parameters for each station show regional variability and need to be adjusted individually to achieve the best performance. We developed a weighting scheme that uses four independent predictors of weight calibrated using a handpicked data subset, which mimics the picking by an expert seismologist. We use the fact that each of the three pickers highlights different properties of the observed seismic trace to combine two pickers that work in tandem. The first makes an initial pick before the second picker refines and improves the accuracy of the automatic pick. We find the tandem pickers improve the accuracy of the automatic picks when compared to the single automatic pickers. We demonstrate that following the cost minimization procedure described here the automatic picks have sufficient accuracy that they would be suitable for high-precision earthquake location, focal mechanism determination or high-resolution seismic tomography.