The permeability of fractures, including natural and hydraulic, are essential parameters for the modeling of fluid flow in conventional and unconventional fractured reservoirs. Traditional analytical (Cubic Law-based) models used to estimate fracture permeability show unstable performance when dealing with different complexities of fracture cases. This work presents a data-driven, physics included model based on machine learning as an alternative to traditional methods. The workflow for the development of the data-driven model includes four steps. 1. Identify uncertain parameters and perform Latin Hypercube Sampling (LHS). We first identify the uncertain parameters that affect the fracture permeability. We then generate training samples using LHS. 2. Perform training simulations and collect inputs and outputs. In this step, high-resolution simulations with parallel computing for the Navier-Stokes (NS) equations are run for each of the training samples. We then collect the inputs and outputs from the simulations. 3. Construct an optimized data-driven surrogate model. A data-driven model based on machine learning is then built to model the nonlinear mapping between the inputs and outputs collected from Step 2. Herein, artificial neural network (ANN) coupling with the Bayesian optimization algorithm is implemented to obtain the optimized surrogate model. 4. Validate the proposed data-driven model. In this step, we conduct blind validation on the proposed model with high fidelity simulations. We further test the developed surrogate model with newly generated fracture cases with a broad range of roughness and tortuosity under different Reynolds numbers. We then compare its performance to the reference NS equation solutions. Results show that the developed data-driven model delivers a good accuracy, exceeding 90% for all training, validation, and test samples. This work introduces an integrated workflow for developing a data-driven, physics included model using machine learning to estimate fracture permeability under complex physics, e.g., inertial effect. To our knowledge, this technique is introduced for the first time for the upscaling of rock fractures. The proposed model offers an efficient and accurate alternative to the traditional upscaling methods that can be readily implemented in reservoir characterization and modeling workflows. © 2022 Aramco Services Company. All rights reserved.