Equal distribution of load among fibrils in contact with a substrate is an important characteristic of fibrillar structures used by many small animals and insects for contact and adhesion. This is in contrast with continuum systems where stress concentration dominates interfacial failure. In this work, we study how adhesion strength of a fibrillar system depends on substrate roughness and variability of the fibril structure, which are modelled using probability distributions for fibril length and fibril attachment strength. Monte Carlo simulations are carried out to determine the adhesion strength statistics where fibril length follows normal or uniform distribution and attachment strength has a power-law form. Our results indicate that the strength distribution is Gaussian (normal) for both the uniform and the normal distributions for length. However, the fibrillar structure having normally distributed lengths has higher strength and lower toughness than one having uniformly distributed lengths. Our simulations also show that an increase in the compliance of the fibrils can compensate for both the substrate roughness and the attachment strength variation. We also show that, as the number of fibrils n increases, the load-carrying efficiency of each fibril goes down. For large n, this effect is found to be small. Furthermore, this effect is compensated by the fact that the standard deviation of the adhesive strength decreases as 1/ square root n.