First I would add a 5th microphone, so that the maths are simpler. 
(edit: https://bioacoustics.stackexchange.com/questions/74/what-are-the-minimum-number-of-micro-hydro-phones-needed-to-perform-3d-tracking)

I would measure (as precise as you can; e.g. with a laser) the separation between each pair of microphones and the distance from multiple reference points.

Then we need a little bit of maths to solve the multiple measurements together with the dime delays into a source location.

Traditionally, if the locations are known as x,y,z then maths is straight forward and discussed also in this SE. The interesting aspect is to use measured microphone distances to estimate also microphone locations. I have not done that yet, but may try.

Edit: it all translate to triangular geometry

Edit2: here comes a basic outline of the procedure

[![enter image description here][1]][1]


  [1]: https://i.sstatic.net/tfEb4.png