The application of geostatistical inversion in shale lithofacies prediction: a case study of the Lower Silurian Longmaxi marine shale in Fuling area in the southeast Sichuan Basin, China

Based on cores, well logs and seismic data, we established the isochronous sequence stratigraphic framework of the Lower Silurian Longmaxi Formation and predicted the shale lithofacies distribution within the sequence stratigraphic framework using geostatistical inversion. The results of our study show that the Lower Member of the Longmaxi Formation is a third order sequence that includes a transgressive systems tract (TST), an early highstand systems tract (EHST) and a late highstand systems tract (LHST). Four lithofacies units have been recognized, specifically siliceous shale, argillaceous shale, calcareous shale and mixed shale. The results of geostatistical inversion reveal that the TST is characterized by flaky siliceous shale and some sparsely distributed calcareous shale. The EHST is dominated by mixed shale with minor amounts of siliceous shale, which occurs in only a small area. Moreover, in the LHST, argillaceous shale occupies almost the entire study region. Comparing to traditional geological research with geophysical research, the vertical resolution of the predictive results of geostatistical inversion could reach 1–2 m. Geostatistical inversion effectively solves the problem of precisely identifying the lithofacies in the Fuling shale gas field and predicting their spatial distribution. This successful study showcases the potential of this method for carrying out marine shale lithofacies prediction in China and other locations with similar geological backgrounds.