Convergence of a Numerical Scheme for Stratigraphic Modeling

In this paper, we consider a multilithology diffusion model used in the field of stratigraphic basin simulations to simulate large scale depositional transport processes of sediments described as a mixture of L lithologies. This model is a simplified one for which the surficial fluxes are proportional to the slope of the topography and to a lithology fraction with unitary diffusion coefficients. The main variables of the system are the sediment thickness h, the L surface concentrations cis in lithology i of the sediments at the top of the basin, and the L concentrations ci in lithology i in the sediments inside the basin. For this simplified model, the sediment thickness decouples from the other unknowns and satisfies a linear parabolic equation. The remaining equations account for the mass conservation of the lithologies, and couple, for each lithology, a first order linear equation for cis with a linear advection equation for ci for which cis appears as an input boundary condition. For this coupled system, a weak formulation is introduced. The system is discretized by an implicit time integration and a cell centered finite volume method. This numerical scheme is shown to satisfy stability estimates and to converge, up to a subsequence, to a weak solution of the problem.

[1]  Roland Masson,et al.  Multi‐lithology stratigraphic model under maximum erosion rate constraint , 2004 .

[2]  O. Ladyženskaja.  Linear and Quasilinear Equations of Parabolic Type , 1968 .

[3]  John W. Harbaugh,et al.  Simulating Clastic Sedimentation , 1989 .

[4]  G. Tucker,et al.  Erosional dynamics, flexural isostasy, and long-lived escarpments: A numerical modeling study , 1994 .

[5]  D. L. Turcotte,et al.  Morphology of a delta prograding by bulk sediment transport , 1985 .

[6]  Jan C. Rivenæs.  Impact of sediment transport efficiency on large‐scale sequence architecture: results from stratigraphic computer simulation , 1997 .

[7]  D. W. Peaceman.  Fundamentals of numerical reservoir simulation , 1977 .

[8]  Didier Granjeon,et al.  Modélisation stratigraphique déterministe : conception et applications d'un modèle diffusif 3d multilithologique , 1996 .

[9]  Didier Granjeon,et al.  Concepts and Applications of a 3-D Multiple Lithology, Diffusive Model in Stratigraphic Modeling , 1999 .

[10]  Jan C. Rivenæs.  Application of a dual‐lithology, depth‐dependent diffusion equation in stratigraphic simulation , 1992 .

[11]  Marie-Hélène Vignal.  Convergence of a finite volume scheme for an elliptic-hyperbolic system , 1996 .

[12]  G. Chavent.  Mathematical models and finite elements for reservoir simulation , 1986 .

[13]  Marie Hélène Vignal,et al.  Numerical and theoretical study of a Dual Mesh Method using finite volume schemes for two phase flow problems in porous media , 1998, Numerische Mathematik.

[14]  Peter B. Flemings,et al.  A synthetic stratigraphic model of foreland basin development , 1989 .

[15]  Todd Arbogast,et al.  A Nonlinear Mixed Finite Eelement Method for a Degenerate Parabolic Equation Arising in Flow in Porous Media , 1996 .

[16]  R. Eymard,et al.  Finite Volume Methods , 2019, Computational Methods for Fluid Dynamics.

[17]  Thierry Gallouët,et al.  Convergence d'un schéma de type éléments finis-volumes finis pour un système formé d'une équation elliptique et d'une équation hyperbolique , 1993 .