Conventional back-projection (BP) method uses 1-D earth model to image rupture process of earthquakes. This may result in discrepancies while using data from arrays with different azimuths. In our study, we integrate a 3-D earth model into BP method to obtain more consistent results when data from different arrays are used for a given earthquake.