Over the last few decades, shale gas has become an increasingly important global source of natural gas, especially in the United States. According to Polczer (2009) and Krauss (2009), shale gas will greatly expand worldwide and is expected to supply as much as half the natural gas production in North America by 2020. Due to extremely low shale matrix permeability, shale is considered an unconventional source of gas and requires fractures to provide a flow path to the wellbore. Due to uncertainties in quantifying the gas-in-place and identifying flow behavior; estimating the ultimate recoveries in shale gas reservoirs requires new techniques. In this paper, we used four approaches to estimate the ultimate recovery in shale gas wells: two empirical methods (conventional and modified decline curve analysis), analytical modeling and numerical modeling. All four approaches were applied on wells from four different shale plays (Barnett, Haynesville, Marcellus and Woodford).