What is the best way to get the phase between noisy signals
The obvious answer would be to take the FFT of each signal. The phase angle of the highest (magnitude) bin gives you the phase of the fundamental sinewave buried in each signal. Compare and contrast as needed.

Use an appropriate window function to minimize edge effects.