The fractional Feynman-Kac equations describe the distributions of functionals of non-Brownian motion,or anomalous diffusion,including two types called the forward and backward fractional Feynman-Kac equations,where the nonlocal time-space coupled fractional substantial derivative is involved.This paper focuses on the more widely used backward version.Based on the newly proposed approximation operators for fractional substantial derivative,we establish compact finite difference schemes for the backward fractional Feynman-Kac equation.The proposed difference schemes have the q-th (q =1,2,3,4) order accuracy in temporal direction and fourth order accuracy in spatial direction,respectively.The numerical stability and convergence in the maximum norm are proved for the first order time discretization scheme by the discrete energy method,where an inner product in complex space is introduced.Finally,extensive numerical experiments are carried out to verify the availability and superiority of the algorithms.Also,simulations of the backward fractional Feynman-Kac equation with Dirac delta function as the initial condition are performed to further confirm the effectiveness of the proposed methods.