In this work, we consider Richardson extrapolation of the Euler scheme for backward stochastic differential equations (BSDEs). First, applying the Adomian decomposition to the nonlinear generator of BSDEs, we introduce a new system of BSDEs. Then we theoretically prove that the solution of the Euler scheme for BSDEs admits an asymptotic expansion, in which the coefficients in the expansions are the solutions of the system. Based on the expansion, we propose Richardson extrapolation algorithms for solving BSDEs. Finally, some numerical tests are carried out to verify our theoretical conclusions and to show the stability, efficiency and high accuracy of the algorithms.