Projection on convex sets with equality and inequality constraints

Question

I want to find the projection of a vector called "a" on a closed and convex set with linear constraints. The set is in the following form: \begin{array}{ll} & Ax = b \\ & Bx \le d \\ &x \ge 0. \end{array} I know I need to solve the following problem to get the projection: \begin{array}{ll} \text{minimize} & || x-a||^2 \\ & Ax = b \\ & Bx \le d \\ &x \ge 0. \end{array} However, I am looking for an efficient way (instead of solving with commercial software like CPLEX) using the structure, convexity, and linearity of the set. I was wondering if anyone knows how to find the projection on this type of feasible region and space properties efficiently. Any resource that has addressed this issue will also help a lot. Note that I found this "Projection of $z$ onto the affine set $\{x\mid Ax = b\}$" here, but the KKT condition requires other conditions in my problem that lead to non-linear constraints. Thanks in advance for your help.