Effect of fulvic acid concentration on modeling electrostatic and heterogeneity effects in proton binding reactions