In this study, we examine a particular type of fractional \((\phi _1,\phi _2)\) -Laplacian problem, which can be represented by the following equation: \(\begin{aligned} \displaystyle {\left\{ \begin{array}{ll} (-\Delta )^{\textrm{s}}_{\phi _1(\cdot )}\hbox {u}+(-\Delta )^{\textrm{s}}_{\phi _{2}(\cdot )}\hbox {u}+\hbox {V}(\epsilon \hbox {x})(\phi _1(|\hbox {u}|)\hbox {u}+\phi _2(|\hbox {u}|)\hbox {u})=\mathrm{f(x,u)} &{} \text{ in } \mathbb {R}^{\textrm{n}},\\ u\in W_\epsilon ,\quad u>0\quad \text{ in } \mathbb {R}^{\textrm{n}}. \end{array}\right. } \end{aligned}\) Here, \(s\in (0, 1)\) , \(\epsilon > 0\) , and \((-\Delta )^s_{\phi _i(.)}\) (where \(i=1,2\) ) is a fractional \(\phi _i\) -Laplacian operator. The potential function \(V:\mathbb {R}^n\rightarrow \mathbb {R}\) is a continuous, possibly unbounded, function and \(f:\mathbb {R}\rightarrow \mathbb {R}\) is a continuous nonlinearity with Orlicz subcritical growth. Our goal is to investigate the multiplicity and concentration properties of the solutions to this problem as \(\epsilon \) approaches zero, using the Nehari manifold and the penalization method as our primary tools.