A fast numerical method and optimization of 3D discrete fracture network considering fracture aperture heterogeneity. (April 2022)