Well, you could always do it via a coordinate transformation to a hypercube, at which point you could call HCubature.jl or similar (which has the advantage of being an adaptive scheme rather than fixed-order, though there are published works on adaptive quadrature for simplices).
See also Integral of a function inside non-rectangular domain - #9 by stevengj