A high-efficiency 3D boundary element method for estimating the stress/displacement field induced by complex fracture networks