Measurements of strong gravitational lensing jointly with type Ia supernovae (SNe Ia) observations have been used to test the validity of the cosmic distance duality relation (CDDR), $D_L(z)/[(1+z)^2D_A(z)]=\eta=1$, where $D_L(z)$ and $D_A(z)$ are the luminosity and the angular diameter distances to a given redshift $z$, respectively. However, several lensing systems lie in the interval $1.4 \leq z \leq 3.6$ i.e., beyond the redshift range of current SNe Ia compilations ($z \approx 1.50$), which prevents this kind of test to be fully explored. In this paper, we circumvent this problem by testing the CDDR considering observations of strong gravitational lensing along with SNe Ia and { a subsample from} the latest gamma-ray burst distance modulus data, whose redshift range is $0.033 \leq z \leq 9.3$. { We parameterize their luminosity distances with a second degree polynomial function and search for possible deviations from the CDDR validity by using four different $\eta(z)$ functions: $\eta(z)=1+\eta_0z$, $\eta(z)=1+\eta_0z/(1+z)$, $\eta(z)=(1+z)^{\eta_0}$ and $\eta(z)=1+\eta_0\ln(1+z)$. Unlike previous tests done at redshifts lower than $1.50$, the likelihood for $\eta_0$ depends strongly on the $\eta(z)$ function considered, but we find no significant deviation from the CDDR validity ($\eta_0=0$). However, our analyses also point to the fact that caution is needed when one fits data in higher redshifts to test the CDDR as well as a better understanding of the mass distribution of lenses also is required for more accurate results.