Fast interpolation method for surfaces with faults by multi-scale second-derivative optimization
IFP Energies nouvelles, 1 et 4 avenue de Bois-Préau, 92852 Rueil-Malmaison Cedex, France
* Corresponding author: firstname.lastname@example.org
Accepted: 10 July 2020
We present a smooth surface interpolation method enabling to take discontinuities (e.g. faults) into account that can be applied to any dataset defined on a regular mesh. We use a second-derivative multi-scale minimization based on a conjugate gradient method. Our multi-scale approach allows the algorithm to process millions of points in a few seconds on a single-unit workstation. The interpolated surface is continuous, as well as its first derivative, except on some lines that have been specified as discontinuities. Application in geosciences are numerous, for instance when a structural model is to be built from points picked on seismic data. The resulting dip of interpolation extends the dip of the input data. The algorithm also works if faults are given by broken lines. We present results from a synthetic and real examples taking into account fault network.
© M. Léger & V. Clochard, published by IFP Energies nouvelles, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.