We study homogenization of iterated randomized singular integrals and homeomorphic solutions to the Beltrami differential equation with a random Beltrami coefficient. More precisely, let $(F_j)_{j \geq 1}$ be a sequence of normalized homeomorphic solutions to the planar Beltrami equation $\overline{\partial} F_j (z)=μ_j(z,ω) \partial F_j(z),$ where the random dilatation satisfies $|μ_j|\leq k<1$ and has locally periodic statistics, for example of the type $$μ_j (z,ω)=φ(z)\sum_{n\in \mathbf{Z}^2}g(2^j z-n,X_{n}(ω)), $$ where $g(z,ω)$ decays rapidly in $z$, the random variables $X_{n}$ are i.i.d., and $φ\in C^\infty_0$. We establish the almost sure and local uniform convergence as $j\to\infty$ of the maps $F_j$ to a deterministic quasiconformal limit $F_\infty$. This result is obtained as an application of our main theorem, which deals with homogenization of iterated randomized singular integrals. As a special case of our theorem, let $T_1,\ldots , T_{m}$ be translation and dilation invariant singular integrals on ${\bf R}^d, $ and consider a $d$-dimensional version of $μ_j$, e.g., as defined above or within a more general setting. We then prove that there is a deterministic function $f$ such that almost surely as $j\to\infty$, $$ μ_j T_{m}μ_j\ldots T_1μ_j\to f \quad \textrm{weakly in } L^p,\quad 1 < p < \infty\ . $$