This study presents an approximate analytical solution for predicting drawdown temperature transient behaviors of a fully penetrating vertical well in a two-zone radial composite reservoir system. The inner zone may represent a damaged (skin) zone, and the outer (non-skin) zone represents an infinitely extended reservoir. The analytical solution is obtained by solving the decoupled isothermal pressure diffusivity equation and temperature equation for the inner and outer zones with the Boltzmann transformation. The convection, transient adiabatic expansion and Joule-Thomson heating effects are all accounted for in the solution. The developed solution compares well with the results of a thermal numerical simulator. The analytical solution is used as a forward model for estimating the parameters of interest by nonlinear regression built on a gradient-based maximum likelihood estimation (MLE) method. A methodology, based on semilog analyses of pressure and temperature data as well as log-log diagnostic plots of pressure- and temperature-derivative data, is proposed to obtain good initial guesses of parameters, which derive the MLE objective function to have reliable optimized estimates. The statistical measures such as estimated standard deviation of noise in pressure and temperature data, confidence intervals for parameters and correlation coefficients between parameter pairs are used to evaluate the goodness of fit and reliability of the estimated parameters from history matching pressure and/or temperature data sets. The results show that the rock, fluid and thermal properties of the skin zone and non-skin zone can be reliably estimated by regressing on temperature transient data jointly with pressure transient data in the presence of noise.