Astrometric accuracy of complex modern VLBI arrays cannot be calculated analytically. We study the astrometric accuracy of phase-referenced VLBI observations for the VLBA, EVN and global VLBI array by simulating VLBI data for targets at declinations -25$^circ$, 0$^circ$, 25$^circ$, 50$^circ$, 75$^circ$ and 85$^circ$. The systematic error components considered in this study are calibrator position, station coordinate, Earth orientation and troposphere parameter uncertainties. We provide complete tables of the astrometric accuracies of these arrays for a source separation of 1$^circ$ either along the right ascension axis or along the declination axis. Astrometric accuracy is 50microas at mid declination and is 300microas at low (-25$^circ$) and high (85$^circ$) declinations for the VLBA and EVN. In extending our simulations to source separations of 0.5$^circ$ and 2$^circ$, we establish the formula for the astrometric accuracy of the VLBA: Delta = (Delta_1$^circ$-14)*d+ 14 (microas) where Delta_1$^circ$ is the astrometric accuracy for a separation d=1$^circ$ provided in our tables for various declinations and conditions of the wet troposphere. We argue that this formula is also valid for the astrometric accuracy of the EVN and global VLBI array.