Naturally fractured gas reservoirs have contributed significantly to global gas reserves and production. The classical gas-well decline analysis relies largely on Arps’ empirical decline models, or modern production decline analysis associating with pseudo-variables. The explicit original gas in place determination methodology is extended from homogeneous reservoir to naturally fractured reservoir under constant or variable bottom-hole pressure conditions in gas-well rate decline analysis. Then, the relationship between gas flow rate and average reservoir pseudo-pressure in the boundary-dominated flow period is re-derived. This formula is in the same format with the equation for homogeneous reservoir by due to the introduction of a new productivity index parameter that captures the inter-porosity flow between fracture and matrix in the natural fractured reservoir. The proposed step-by-step procedures are applied here, which enable the estimation of decline exponent and the explicit and straightforward determination of the original gas in place without any iterative calculations. Four simulated cases prove that our methodology can be successfully used in heterogeneous naturally fractured reservoirs with irregular boundary under constant or variable bottom-hole pressure conditions.