In the bond percolation model on a lattice, we colour vertices with colours independently at random according to Bernoulli distributions. A vertex can receive multiple colours and each of these colours is individually observable. The colours colour the entire component into which they fall. Our goal is to estimate the parameters of the model: the probabilities of colouring of single vertices and the probability with which an edge is open. The input data is the configuration of colours once the complete components have been coloured, without the information which vertices were originally coloured or which edges are open. We use a Monte Carlo method, the method of simulated moments, to achieve this goal. Under the unproven assumption of identifiability, we show that this method is a strongly consistent estimator by proving a uniform strong law of large numbers for the vertices' weakly dependent colour values. Our proof method quantifies dependence among the spatially arranged random variables using percolation theory: the FKG and BK inequalities, and the exponential decay of the cluster size distribution. We evaluate the method in computer tests. We have made our software publicly available. The motivating application is cross-contamination rate estimation for digital PCR in lab-on-a-chip microfluidic devices. [ABSTRACT FROM AUTHOR]