Numerical solution to binary quasi-linear partial differential equations for drilling surge pressures based on the characteristic-lines method