In this work, we focus our attention on the existence of nontrivial solutions to the following supercritical Schrödinger-Poisson type system with $(p,q)$-Laplacian:

$$\{\begin{array}{ll}-{\mathrm{\Delta }}_{p}u-{\mathrm{\Delta }}_{q}u+\varphi |u{|}^{q-2}u=f\left(x,u\right)+\mu |u{|}^{s-2}u& \text{in}\mathrm{\Omega },\\ -\mathrm{\Delta }\varphi =|u{|}^q& \text{in}\mathrm{\Omega },\\ u=\varphi =0& \text{on}\mathrm{\partial }\mathrm{\Omega },\end{array}$$

where $\mathrm{\Omega }\subset {\mathbb{R}}^N$ is a bounded smooth domain, $\mu >0,N>1$, and $-{\mathrm{\Delta }}_{\mathrm{\wp }}\phi =div(|\mathrm{\nabla }\phi {|}^{\mathrm{\wp }-2}\mathrm{\nabla }\phi )$, with $\mathrm{\wp }\in \{p,q\}$, is the homogeneous $\mathrm{\wp }$-Laplacian. $1<p<q<\frac{q^{\ast }}{2}$, $q^{\ast }:=\frac{Nq}{N-q}<s$, and $q^{\ast }$ is the critical exponent to $q$. The proof is accomplished by the Moser iterative method, the mountain pass theorem, and the truncation technique. Furthermore, the $(p,q)$-Laplacian and the supercritical term appear simultaneously, which is the main innovation and difficulty of this paper.

In this article, we consider the following a class of Schrödinger-Poisson systems with $p$-Laplacian in ${\mathbb{R}}^3$ of the form:

$$\{\begin{array}{ll}-{\mathrm{\Delta }}_{p}u+(\lambda b(x)+1)|u{|}^{p-2}u+\varphi |u{|}^{s-2}u=g(u)& \text{in}{\mathbb{R}}^3,\\ -\mathrm{\Delta }\varphi =|u{|}^s& \text{in}{\mathbb{R}}^3,\end{array}$$

where $1<p<3$, $\frac{p}{2}<s<p$, ${\mathrm{\Delta }}_{p}u:=div(|\mathrm{\nabla }u{|}^{p-2}\mathrm{\nabla }u)$ is the $p$-Laplacian operator, $\lambda$ is a positive parameter. Assume that the nonnegative function $b$ possesses a potential well int $(b^{-1}(\{0\}))$, which is composed of $k$ disjoint components ${\mathrm{\Omega }}_1,{\mathrm{\Omega }}_2,\cdots ,{\mathrm{\Omega }}_k$ and consider the nonlinearity $g$ with subcritical growth. Using the variational methods and Morse iteration technique, the existence of positive multi-bump solutions are obtained.