Reachable set estimation for discrete-time bilinear systems with time-varying delays