We present detailed velocity and density models for the Montney unconventional resource play in northeast British Columbia, Canada. The new models are specifically essential for robust hypocenter determination in the areas undergoing multistage hydraulic‐fracturing operations and for detailed analysis of induced seismicity processes in the region. For the upper 4 km of the sedimentary structure, we review hundreds of well logs and select sonic and density logs from 19 locations to build the representative models. For depths below 4 km, we extend our models using data from the southern Alberta refraction experiment (Clowes et al., 2002). We provide one set of models for the entire Montney play along with two separated sets for the southern and northern areas. Specifically, the models for the southern and northern Montney play are based on logs located in and around the Kiskatinaw Seismic Monitoring and Mitigation Area and the North Peace Ground Motion Monitoring area, respectively. To demonstrate the usefulness of our detailed velocity model, we compare the hypocenter location of earthquakes that occurred within the Montney play as determined with our model and the simple two‐layered model (CN01) routinely used by Natural Resources Canada. Locations obtained by our velocity model cluster more tightly with the majority of events having root mean square residual of compared with that of when the CN01 model is used. Cross sections of seismicity versus depth across the area also show significant improvements in the determination of focal depths. Our model results in a reasonable median focal depth of for events in this area, which is consistent with the completion depths of hydraulic‐fracturing operations. In comparison, most solutions determined with the CN01 model have fixed focal depths (0 km) due to the lack of depth resolution.