In this paper, a novel method based on block pulse functions is proposed to identify continuous-time fractional-order systems with time delays. First, the operational matrices of block pulse functions for fractional integral operator and time delay operator are derived. Then, these operational matrices are applied to convert the continuous-time fractional-order systems with time delays to an algebraic equation. Finally, the system’s parameters along with the differentiation orders and the time delays are all simultaneously estimated through minimizing a quadric error function. The proposed method reduces the computation complexity of the identification process, and also it does not require the system’s differentiation orders to be commensurate. The effectiveness of the proposed method are demonstrated by several numerical examples.