We develop a new method to locate microseismic events induced by hydraulic fracturing with simultaneous anisotropic tomography, using differential arrival times and differential backazimuths. Compared to the existing double-difference method, our method incorporates backazimuth information to better constrain microseismic locations in the case of downhole linear seismic arrays used for monitoring induced seismicity. The tomography is constrained to a 1-D layered VTI (transversely isotropic structure with a vertical symmetry axis) structure to improve inversion stability given the limited passive seismic data. We derive analytical sensitivities for the elastic moduli (C ij) and layer thickness L, and verify the analytical results with numerical calculations. The forward modelled traveltimes and sensitivities are all calculated analytically without weak anisotropy assumption. By incorporating the relative information among events, the extended double-difference method can provide better relative locations for events and, therefore, can characterize the fractures with higher accuracy. In the two tests with synthetic data, our method provides more accurate relative locations than the traditional methods, which only use absolute information. With fast speed and high accuracy, our inversion scheme is suitable for real-time microseismic monitoring of hydraulic fracturing.