A high order method for orbital conjunctions analysis: Monte Carlo collision probability computation