The current-voltage characteristics of an organic heterojunction (HJ), where random site trap density at the donor-acceptor (D-A) interface is included, is numerically investigated based on the polaronpair (PP) model (Phys. Rev. B. 82, 155305 (2010)). To derive the analytic equation, the electric field FI at the D-A interface, which is included in the field parameter b, is first approximated to estimate the polaronpair dissociation rate, k(ppd). It is then assumed that the quasi-Fermi level E-Fp/n at the donor/acceptor lies in the trap energy level (E-trap) at random sites without a direct dependence on applied bias (V-a), but on a HOMO/LUMO with V-a dependence. As a result, the low diode current (DC) at forward bias is attributed to the Etrap, which is lower than the LUMO at acceptors, and is higher than the HOMO at donors. This implies that the pattern of DC is at forward biases depending upon the position of E-trap