Blood oxygen level-dependent functional magnetic resonance imaging (BOLD fMRI) arises from a physiological and physical cascade of events taking place at the level of the cortical microvasculature which constitutes a medium with complex geometry. Several analytical models of the BOLD contrast have been developed, but these have not been compared directly against detailed bottom-up modeling methods. Using a 3D modeling method based on experimentally measured images of mice microvasculature and Monte Carlo simulations, we quantified the accuracy of two analytical models to predict the amplitude of the BOLD response from 1.5 to 7T, for different echo time (TE) and for both gradient echo and spin echo acquisition protocols. We also showed that accounting for the tridimensional structure of the microvasculature results in more accurate prediction of the BOLD amplitude, even if the values for SO2 were averaged across individual vascular compartments. A secondary finding is that modeling the venous compartment as two individual compartments results in more accurate prediction of the BOLD amplitude compared with standard homogenous venous modeling, arising from the bimodal distribution of venous SO2 across the microvasculature in our data.