Analytical solutions for two-phase subsurface flow to a leaky fault considering vertical flow effects and fault properties