A new method for recovering paleoporosity of sandstone: case study of middle Es3 member of Paleogene formation in Niuzhuang Sag, Dongying Depression, Bohai Bay Basin in China