Numerical modeling of solid-phase microextraction: binding matrix effect on equilibrium time