Let \(\Omega\subset\mathbb{R}^{N}\) be a bounded domain with smooth boundary and \(0\in\Omega\). For \(0<s<1\), \(1\le r<q<p\), \(0\le\alpha<ps<N\) and a positive parameter \(\lambda\), we consider the fractional \((p,q)\)-Laplacian problems involving a critical Sobolev-Hardy exponent. This model comes from a nonlocal problem of Kirchhoff type $$\displaylines{ \big(a+b[u]_{s,p}^{(\theta-1)p}\big)(-\Delta)_{p}^{s}u+(-\Delta)_{q}^{s}u =\frac{|u|^{p_{s}^{*}(\alpha)-2}u}{|x|^{\alpha}}+\lambda f(x)\frac{|u|^{r-2}u}{|x|^{c}}\quad \hbox{in }\Omega,\cr u=0\quad\text{in }\mathbb{R}^{N}\setminus\Omega, }$$ where \(a,b>0\), \(c<sr+N(1-r/p)\), \(\theta\in(1,p_{s}^{*}(\alpha)/p)\) and \(p_{s}^{*}(\alpha)\) is critical Sobolev-Hardy exponent. For a given suitable \(f(x)\), we prove that there are least two nontrivial solutions for small \(\lambda\), by way of the mountain pass theorem and Ekeland's variational principle. Furthermore, we prove that these two solutions converge to two solutions of the limiting problem as \(a\to 0^{+}\). For the limiting problem, we show the existence of infinitely many solutions, and the sequence tends to zero when \(\lambda\) belongs to a suitable range. For more information see https://ejde.math.txstate.edu/Volumes/2021/66/abstr.html