-This paper presents production data analysis for two production wells located in the shale gas field, Canada, with the proper analysis method according to each production performance characteristics. In the case A production well, the analysis was performed by applying both time and superposition time because the production history has high variation. Firstly, the flow regimes were classified with a log-log plot, and as a result, only the transient flow was appeared. Then the area of simulated reservoir volume (SRV) analyzed based on flowing material balance plot was calculated to 180 acres of time, and 240 acres of superposition time. And the original gas in place (OGIP) also was estimated to 15, 20 Bscf, respectively. However, as the area of SRV was not analyzed with the boundary dominated flow data, it was regarded as the minimum one. Therefore, the production forecasting was conducted according to variation of b exponent and the area of SRV. As a result, estimated ultimate recovery (EUR) increased 1.2 and 1.4 times respectively depending on b exponent, which was 0.5 and 1. In addition, as the area of SRV increased from 240 to 360 acres, EUR increased 1.3 times.