A waveray method has been used in past studies in order to evaluate tsunami hazards. Though this method is useful for it's convenience, it is difficult to predict quantitative tsunami heights because of the problems of calculating method, such as divergences of refraction coefficient in the cross-section areas of waverays. In this study, to be able to use waveray method for tsunami height prediction, (1) we modeled the shape of wave sources and investigated the relationship between fault parameters, (2) we improved the method of calculating the refraction coefficient with waveray method.