A C program based on the diffusion equation that simulates three-dimensiona
l (3-D) basin infilling is described. The model solves the sets of non-line
ar differential equations by a fully implicit finite-difference method usin
g a Newton's scheme. DIBAFILL (DIffusion BAsin inFILL) can be employed for
the prediction of 3-D morphology, sediment accumulation maps, and fraction
distributions, as well as for studying process-response systems in sediment
ary basins. The sensitivity of the diffusion coefficient to the lithology,
depositional environment and grain size, makes it possible to adapt the pro
gram for a wide range of shallow-water marine depositional systems, includi
ng siliciclastic, mixed siliciclastic/carbonate and carbonate environments.
The program includes an option to run the model with a depth-dependent car
bonate production law. Applications to idealized basins are presented to sh
ow the capabilities of the program at various spatial scales (reservoir and
basin scales), and for different depositional environments (carbonate and
siliciclastic). (C) 2000 Elsevier Science Ltd. All rights reserved.