Since you are using cplex, there is no need to linearize really, all you need is a second-order cone model.
Write as $$(p_i y_i - \sum_{j \in J}b_{j} x_{ji}) \le t,~t \leq \sqrt{\sum_{j \in J} x_{ji}^2 \sigma_j^2}$$ and square the nonlinear term (possible as the problematic term is non-negative, so there is no loss in generality to assume $t\geq 0$) $$t^2 \leq \sum_{j \in J} x_{ji}^2 \sigma_j^2$$ and use $x_{ji}^2 = x_{ji}$. From that, it follows that you have a convex quadratic constraint, i.e. second-order cone representable.