An equivalent fracture length-based numerical method for modeling nonlinear flow in 2D fracture networks