s_w \, V_{mw} \cdot \frac{(\epsilon -1)(2\epsilon+1)}{9\epsilon} + \left( 1-s_ow \, V_{mo} + s_g \, V_{mg} \right) \cdot \frac{\epsilon-1}{\epsilon+2} = P = s_w P_w + s_o P_g + s_g P_g

P_w = V_{mw} \cdot \frac{(\epsilon_w-1)(2\epsilon_w+1)}{9 \, \epsilon_w}

P_o(T) = V_{mw} \cdot \frac{\epsilon_o-1}{\epsilon_o+2}

P_g(T) = V_{mw} \cdot\frac{\epsilon_g-1}{\epsilon_g+2}

\epsilon_o(T) = 16 ÷ 320

\epsilon_og(T) = 1 ÷ 2