The B6 metamorphic buried hill condensate gas reservoir exhibits a highly compact matrix, leading to a rapid decline in bottom-hole pressure during initial production. The minimal difference between formation and saturation pressures results in severe retrograde condensation, with multiphase flow further increasing resistance. Conventional numerical simulations often overestimate reservoir energy supply due to their failure to account for this additional resistance, leading to inaccuracies in bottom-hole pressure predictions and gas–oil ratio during history matching. To address these challenges, this study conducted research on nonlinear numerical simulation for buried hill condensate gas reservoirs and established a method for calculating a multiphase pressure sweep range based on the well testing theory. By correcting and fitting the pressure propagation boundaries with numerical simulation, the nonlinear flow parameters applicable to the B6 gas field were inversed. This study revealed that conventional Darcy flow is inadequate for predicting pressure propagation boundaries and that it is possible to reasonably characterize the pressure sweep range through nonlinear flow. This approach resulted in an improvement in the accuracy of historical matching for bottom-hole pressure and gas–oil ratio, which improve the historical fitting accuracy to 85%, providing valuable insights for the development of similar reservoirs.