In this article, we establish the existence of nonnegative solutions to the following quasilinear and singular elliptic problems with supercritical nonlinearity: \(\begin{aligned} \left\{ \begin{aligned} {} -\Delta _p z-\Delta _q z&{}= \lambda \frac{h(x)}{z^\gamma }+z^\theta , \ z>0&\quad \text{ in } \, \Omega , \\ z&{}= 0&\quad \text{ on } \partial \Omega , \end{aligned} \right. \end{aligned}\) where \(\Omega \) is an open, bounded subset of \(\mathbb {R}^N (N\ge 3)\) with \(C^2\) boundary, h is a positive real-valued function, \(1<p<q<\infty \) and \(\lambda , \theta , \gamma \) are positive parameters. The main contribution of this paper is the treatment of singular and supercritical nonlinearities on the right-hand side in the presence of the nonhomogeneous \((p+q)\) -Laplace operator. To demonstrate the existence of a weak solution, we utilise the method of sub and supersolution.