A variable-coefficient forced Korteweg-de Vries equation with spacial inhomogeneity is investigated in this paper. Under constraints, this equation is transformed into its bilinear form, and multi-soliton solutions are derived. Effects of spacial inhomogeneity for soliton velocity, width and background are discussed. Nonlinear tunneling for this equation is presented, where the soliton amplitude can be amplified or compressed. Our results might be useful for the relevant problems in fluids and plasmas.