Numerical analysis of sedimentary compaction: Implications for porosity and layer thickness variation

To understand the formation and evolution of a sedimentary basin in basin analysis and modelling studies, it is important to analyze the thickness and age range of sedimentary layers infilling a basin. Because the compaction effect reduces the thickness of sedimentary layers during burial, basin modelling studies typically restore the reduced thickness using the relation of porosity and depth (compaction trend). Based on the compilation plots of published compaction trends of representative sedimentary rocks (sandstone, shale and carbonate), this study estimates the compaction trend ranges with exponential curves and equations. Numerical analysis of sedimentary compaction is performed to evaluate the variation of porosity and layer thickness with depth at key curves within the compaction trend ranges. In sandstone, initial porosity lies in a narrow range and decreases steadily with increasing depth, which results in relatively constant thickness variations. For shale, the porosity variation shows two phases which are fast reduction until ~2,000 m in depth and slow reduction at deeper burial, which corresponds to the thickness variation pattern of shale layers. Carbonate compaction is characterized by widely distributed porosity values, which results in highly varying layer thickness with depth. This numerical compaction analysis presents quantitatively the characteristics of porosity and layer thickness variation of each lithology, which influence on layer thickness reconstruction, subsidence and thermal effect analyses to understand the basin formation and evolution. This work demonstrates that the compaction trend is an important factor in basin modelling and underlines the need for appropriate application of porosity data to produce accurate analysis outcomes.