We develop a linearized boundary control method for the inverse boundary value problem of determining a potential in the acoustic wave equation from the Neumann-to-Dirichlet map. When the linearization is at the zero potential, we derive a reconstruction formula based on the boundary control method and prove that it is of Lipschitz-type stability. When the linearization is at a nonzero potential, we prove that the problem is of H"{o}lder-type stability in two and higher dimensions. The proposed reconstruction formula is implemented and evaluated using several numerical experiments to validate its feasibility.