Fault displacement modelling using 3D vector fields

In history matching and sensitivity analysis, flexibility in the structural modelling is of great importance. The ability to easily generate multiple realizations of the model will have impact both on the updating workflow in history matching and uncertainty studies based on Monte Carlo simulations. The main contribution to fault modelling by the work presented in this paper is a new algorithm for calculating a 3D displacement field applicable to a wide range of faults due to a flexible representation. This gives the possibility to apply this field to change the displacement and thereby moving horizons and fault lines. The fault is modelled by a parametric format where the fault has a reference plane defined by a centre point, dip and strike angles. The fault itself is represented as a surface defined by a function z = f(x,y), where x, y and z are coordinates local to the reference plane, with the z-axis being normal to the plane. The displacement associated with the fault outside the fault surface is described by a 3D vector field. The displacement on the fault surface can be found by identifying the intersection lines between horizons and the fault surface (fault lines), and using kriging techniques to fill in values at all points on the surface. Away from the fault surface the displacement field is defined by a monotonic decreasing function which ensures zero displacement at a specified distance from the fault. An algorithm is developed where the displacement can be increased or decreased according to user-defined parameters. This means that the whole displacement field is changed and points on horizons around the fault can be moved accordingly by applying the modified displacement field on them. The interaction between several faults influencing the same points is taken care of by truncation rules and the ordering of the faults. The method is demonstrated on a realistic synthetic case based on a real reservoir.