Naphthenic acids (NAs) naturally occur in crude oil and its associated produced water, presenting significant challenges, such as corrosion, in refinery apparatus and ecotoxicity in aquatic habitats. This study delineates a multivariate method to quantify NAs in produced water via electrospray ionization coupled with high-resolution Orbitrap mass spectrometry (ESI-Orbitrap MS). By employing liquid-liquid extraction, followed by direct infusion ESI(-)-Orbitrap MS, we characterized and quantified NAs employing a partial least-squares regression (PLS) model enhanced by the ordered predictor selection (OPS) algorithm. Thirty-six produced water samples were utilized, with 24 allocated for calibration and 12 designated for validation. The PLS-OPS model demonstrated notable accuracy in predicting NA concentrations in simulated and actual produced water samples ranging from ∼30 to 300 mg·L-1. This methodology offers a rapid yet robust alternative for quantifying NAs using mass spectrometry augmented by PLS and the OPS. Its significance is underscored by its potential to equip the petroleum industry with a swift and reliable monitoring mechanism for NAs in produced water, thereby aiding in mitigating environmental and operational risks.