To improve the accuracy of microseismic event locations, we developed a new inversion method with double-difference constraints for determining the hypocenters and the anisotropic velocity model for unconventional reservoirs. We applied this method to a microseismic data set monitoring a Middle Bakken completion in the Beaver Lodge area of North Dakota. Geophone arrays in four observation wells improved the ray coverage for the velocity inversion. Using an accurate anisotropic velocity model is important to correctly assess the height growth of the hydraulically induced fractures in the Middle Bakken. Our results showed that (1) moderate-to-strong anisotropy exists in all studied sedimentary layers, especially in the Upper and Lower Bakken shale formations, where the Thomsen parameters ( and ) can be greater than 0.4, (2) all the events selected for high signal-to-noise ratio and used for the joint velocity inversion are located in the Bakken and overlying Lodgepole formations, i.e., no events are detected in the Three Forks formation below the Bakken, and (3) more than half of the strong events are in two clusters at approximately 100 and 150 m above the Middle Bakken. Reoccurrence of strong, closely clustered events suggested activation of natural fractures or faults in the Lodgepole formation. The sensitivity analysis for the inversion results showed that the relative uncertainty in parameter is larger than other anisotropy parameters. The microseismic event locations and the anisotropic velocity model are validated by comparing synthetic and observed seismic waveforms and by S-wave splitting.