×

Implicit structural modeling by minimization of the bending energy with moving least squares functions. (English) Zbl 1421.86033

Summary: In this paper, an implicit structural modeling method using locally defined moving least squares shape functions is proposed. The continuous bending energy is minimized to interpolate between data points and approximate geological structures. This method solves a sparse problem without relying on a complex mesh. Discontinuities such as faults and unconformities are handled with minor modifications of the method using meshless optic principles. The method is illustrated on a two-dimensional model with folds, faults and an unconformity. This model is then modified to show the ability of the method to handle sparsity, noise and different reliabilities in the data. Key parameters of the shape functions and the pertinence of the bending energy for structural modeling applications are discussed. The predefined values deduced from these studies for each parameter of the method can also be used to construct other models.

MSC:

86A60 Geological problems
74G65 Energy minimization in equilibrium problems in solid mechanics
86A32 Geostatistics
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] ARANZ Geo (2018) Leapfrog geo. http://www.leapfrog3d.com/products/leapfrog-geo. Accessed 24 Oct 2018
[2] Aug C (2004) Modélisation géologique 3D et caractérisation des incertitudes par la méthode du champ de potentiel. PhD thesis
[3] Belytschko T, Lu YY, Gu L (1994) Element-free Galerkin methods. Int J Numer Methods Eng 37(2):229-256 · Zbl 0796.73077
[4] Belytschko T, Krongauz Y, Fleming M, Organ D, Liu WKS (1996) Smoothing and accelerated computations in the element free Galerkin method. J Comput Appl Math 74(1):111-126 · Zbl 0862.73058
[5] Brezzi F (1974) On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. Rev Française d’autom Inf Rech Opér Anal Numér 8(R2):129-151 · Zbl 0338.90047
[6] Calcagno P, Chilès JP, Courrioux G, Guillen A (2008) Geological modelling from field data and geological knowledge. Part I. Modelling method coupling 3D potential-field interpolation and geological rules. Phys Earth Planet Inter 171(1-4):147-157
[7] Carmichael T, Ailleres L (2016) Method and analysis for the upscaling of structural data. J Struct Geol 83:121-133
[8] Caumon G (2009) Vers une intégration des incertitudes et des processus en géologie numérique. PhD thesis
[9] Caumon G, Collon-Drouaillet P, De Veslud CLC, Viseur S, Sausse J (2009) Surface-based 3D modeling of geological structures. Math Geosci 41(8):927-945 · Zbl 1178.86002
[10] Caumon G, Gray G, Antoine C, Titeux MO (2013) Three-dimensional implicit stratigraphic model building from remote sensing data on tetrahedral meshes: theory and application to a regional model of la popa basin, ne mexico. IEEE Trans Geosci Remote Sens 51(3):1613-1621
[11] Chilès JP, Aug C, Guillen A, Lees T (2004) Modelling the geometry of geological units and its uncertainty in 3D from structural data: the potential-field method. In: orebody modeling and strategic mine planning—spectrum 14(July):22-24
[12] Collon P, Steckiewicz-Laurent W, Pellerin J, Laurent G, Caumon G, Reichart G, Vaute L (2015) 3D geomodelling combining implicit surfaces and Voronoi-based remeshing: a case study in the Lorraine coal basin (France). Comput Geosci 77:29-43
[13] Courant R (1950) Dirichlet’s principle, conformal mapping, and minimal surfaces. Interscience, New York · Zbl 0040.34603
[14] Cowan E, Beatson R, Ross H, Fright W, McLennan T, Evans T, Carr J, Lane R, Bright D, Gillman A, Oshust P, Titley M (2003) Practical implicit geological modelling. In: 5th international mining geology conference, vol 8, pp 89-99
[15] Cuomo S, Galletti a, Giunta G, Starace A (2013) Surface reconstruction from scattered point via RBF interpolation on GPU. In: 2013 federated conference on computer science and information systems (FedCSIS), pp 433-440. arXiv:1305.5179v1
[16] De la Varga M, Schaaf A, Wellmann F (2019) Gempy 1.0: open-source stochastic geological modeling and inversion. Geosci Model Dev (in review)
[17] Dubrule O (1984) Comparing splines and kriging. Comput Geosci 10(2-3):327-338
[18] Duchon J (1977) Splines minimizing rotation-invariant semi-norms in Sobolev spaces. In: Constructive theory of functions of several variables, Springer, pp 85-100 · Zbl 0342.41012
[19] Fossen H (2016) Structural geology. Cambridge University Press, Cambridge
[20] Frank T, Tertois AL, Mallet JL (2007) 3d-reconstruction of complex geological interfaces from irregularly distributed and noisy point data. Comput Geosci 33(7):932-943
[21] Fries TP, Matthias HG (2004) Classification and overview of meshfree methods. Department of Mathematics and Computer Science. Technical University of Braunschweig, p 64
[22] Godefroy G, Caumon G, Ford M, Laurent G, Jackson CAL (2018) A parametric fault displacement model to introduce kinematic control into modeling faults from sparse data. Interpretation 6(2):B1-B13
[23] Gonçalves ÍG, Kumaira S, Guadagnin F (2017) A machine learning approach to the potential-field method for implicit modeling of geological structures. Comput Geosci 103:173-182
[24] Grose L, Laurent G, Ailleres L, Armit R, Jessell M, Caumon G (2017) Structural data constraints for implicit modeling of folds. J Struct Geol 104:80-92
[25] Hillier MJ, Schetselaar EM, de Kemp EA, Perron G (2014) Three-dimensional modelling of geological surfaces using generalized interpolation with radial basis functions. Math Geosci 46(8):931-953 · Zbl 1323.86042
[26] Houlding SW (1994) 3D geoscience modelling. Computer techniques for geological characterisation, Springer, New York
[27] Intrepid-Geophysics (2018) Geomodeller. http://www.intrepid-geophysics.com/ig/index.php?page=geomodeller. Accessed 24 Oct 2018
[28] Jessell M, Aillères L, Kemp ED, Lindsay M, Wellmann F, Hillier M, Laurent G, Carmichael T, Martin R (2014) Next generation three-dimensional geologic modeling and inversion. SEG Spec Publ 18:261-272
[29] Karimi-Fard M, Durlofsky LJ (2016) A general gridding, discretization, and coarsening methodology for modeling flow in porous formations with discrete geological features. Adv Water Resour 96:354-372
[30] Lajaunie C, Courrioux G, Manuel L (1997) Foliation fields and 3D cartography in geology: principles of a method based on potential interpolation. Math Geol 29(4):571-584
[31] Lancaster P, Salkauskas K (1981) Surfaces generated by moving least squares methods. Math Comput 37(155):141-158 · Zbl 0469.41005
[32] Laurent G (2016) Iterative thickness regularization of stratigraphic layers in discrete implicit modeling. Math Geosci 48(7):811-833
[33] Laurent G, Ailleres L, Grose L, Caumon G, Jessell M, Armit R (2016) Implicit modeling of folds and overprinting deformation. Earth Planet Sci Lett 456:26-38
[34] Liu GR, Gu Y (2005) An introduction to meshfree method and their programming. Springer, New York
[35] Mallet JL (1988) Three-dimensional graphic display of disconnected bodies. Math Geol 20(8):977-990
[36] Mallet JL (1992) Discrete smooth interpolation. Comput Aided Des 24(4):178-191 · Zbl 0808.65005
[37] Mallet JL (1997) Discrete modeling for natural objects. Math Geol 29(2):199-219 · Zbl 0970.86505
[38] Mallet JL (2002) Geomodeling. Oxford University Press Inc, Oxford · Zbl 1009.05001
[39] Mallet JL (2014) Elements of mathematical sedimentary geology: The GeoChron model. EAGE publications, Houten
[40] Marechal, A.; Verly, G. (ed.); David, M. (ed.); Journel, AG (ed.); Marechal, A. (ed.), Kriging seismic data in presence of faults, 271-294 (1984), Dordrecht
[41] Martin R, Boisvert JB (2017) Iterative refinement of implicit boundary models for improved geological feature reproduction. Comput Geosci 109:1-15
[42] Matheron G (1981) Splines and kriging: their formal equivalence. In: Down-to-earth statistics: solutions looking for geological problems, vol 8. Syracuse University Geological Contributions, pp 77-95
[43] McLain DH (1976) Two dimensional interpolation from random data. Comput J 19(2):178-181 · Zbl 0321.65009
[44] Nguyen VP, Rabczuk T, Bordas S, Duflot M (2008) Meshless methods: a review and computer implementation aspects. Math Comput Simul 79(3):763-813 · Zbl 1152.74055
[45] Organ D, Fleming M, Terry T, Belytschko T (1996) Continuous meshless approximations for nonconvex bodies by diffraction and transparency. Comput Mech 18(3):225-235 · Zbl 0864.73076
[46] Paradigm (2018) Skua-gocad. http://www.pdgm.com/products/skua-gocad/. Accessed 24 Oct 2018
[47] Pellerin J, Lévy B, Caumon G, Botella A (2014) Automatic surface remeshing of 3D structural models at specified resolution: a method based on Voronoi diagrams. Comput Geosci 62:103-116
[48] Piegl L, Tiller W (1997) The NURBS book. Monographs in visual communication. Springer, New York
[49] Schlumberger (2018) Petrel. https://www.software.slb.com/products/petrel. Accessed 24 Oct 2018
[50] Shewchuk J (2002) What is a good linear finite element? Interpolation, conditioning, anisotropy, and quality measures (preprint). University of California at Berkeley, vol 73, p 137
[51] Souche L, Iskenova G, Lepage F, Desmarest D, et al (2014) Construction of structurally and stratigraphically consistent structural models using the volume-based modelling technology: applications to an Australian dataset. In: International petroleum technology conference
[52] Wahba G (1990) Spline models for observational data, vol 59. Society for Industrial and Applied Mathematics, Philaedelphia · Zbl 0813.62001
[53] Wendland H (1995) Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. Adv Comput Math 4(1):389-396 · Zbl 0838.41014
[54] Wu Z (1995) Compactly supported positive definite radial functions. Adv Comput Math 4(1):283 · Zbl 0837.41016
[55] Yokota R, Barba LA, Knepley MG (2010) Petrbf—a parallel o(n) algorithm for radial basis function interpolation with Gaussians. Comput Methods Appl Mech Eng 199(25):1793-1804 · Zbl 1231.65026
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.