  We consider the massive graviton phenomenological model based on the
graviton's dispersion terms included into phase of gravitational wave's
waveform. Such model was already considered in many works but it was based on a
single leading-order dispersion term only. Here we derive a relation between
relativistic gravitons emission and absorption time intervals computed up to
${\cal O}(\gamma^{-6})$, where $\gamma$ is the Lorentz factor. Including the
dispersion terms into the phase of gravitational wave's waveform results in two
non-GR parameters of the $1^{st}$ and the $-2^{nd}$ post-Newtonian orders whose
posteriors are used to put a constraint on the graviton's rest mass. We use the
TaylorF2 waveform model to analyse the event GW170817 and report the following
$95\%$-confidence upper bounds on the graviton's rest mass:
$m^{Low\,Spin}_{g}\leq1.305\times10^{-54}$g and
$m^{High\,Spin}_{g}\leq2.996\times10^{-54}$g for the high and low spin priors.
