A method based on polynomial regression algorithm (PRA) is proposed in this paper to compensate the nonlinear phase noise in optical frequency domain reflection (OFDR) systems. In this method, the nonlinear phase of OFDR systems is represented by the polynomial phase function, and then the coefficients of the polynomial phase function are estimated by PRA. Finally, the nonlinearity is compensated by match Fourier transform (MFT). Simulation results demonstrate that the proposed algorithm has good performance in compensating both weak and strong nonlinear phase noises of OFDR systems.