A computer-implemented method for generating a fractured reservoir model executed by one or more processors, the method comprising: receiving a seismic dataset of a surveyed subsurface, the seismic dataset comprising a set of parameters representing physical features of the surveyed subsurface; identifying a sensitivity of each parameter of the set of parameters on a dynamic response of the fractured reservoir model, wherein each parameter of the set is associated with a range of values; selecting a subset of parameters from the set of parameters based on the sensitivity of each parameter; sampling an outer boundary of a parameter uncertainty domain, the parameter uncertainty domain based on the range of values for each parameter of the subset of parameters; adjusting the range of values associated with each parameter of the subset based on the sampling; generating a geo-model based on the adjusted range of values associated with each parameter of the subset; generating a discrete fracture network model based on the geo-model; generating a scenario of a simulation model based on the discrete fracture network model and the geo-model; performing a forward simulation based on the scenario of the simulation model; determining that a misfit of the forward simulation is below a threshold by evaluating an objective function; in response to determining that the misfit of the forward simulation is not below the threshold, redefining the objective function by an assisted history matching process in which a new parameter uncertainty domain is sampled for each of the geo-model, the discrete fracture network model, and the scenario of a simulation model, and generating a history matched model based on the redefined objective function; and producing a 3D model based on the forward simulation. receiving a seismic dataset of a surveyed subsurface, the seismic dataset comprising a set of parameters representing physical features of the surveyed subsurface; identifying a sensitivity of each parameter of the set of parameters on a dynamic response of the fractured reservoir model, wherein each parameter of the set is associated with a range of values; selecting a subset of parameters from the set of parameters based on the sensitivity of each parameter; sampling an outer boundary of a parameter uncertainty domain, the parameter uncertainty domain based on the range of values for each parameter of the subset of parameters; adjusting the range of values associated with each parameter of the subset based on the sampling; generating a geo-model based on the adjusted range of values associated with each parameter of the subset; generating a discrete fracture network model based on the geo-model; generating a scenario of a simulation model based on the discrete fracture network model and the geo-model; performing a forward simulation based on the scenario of the simulation model; determining that a misfit of the forward simulation is below a threshold by evaluating an objective function; in response to determining that the misfit of the forward simulation is not below the threshold, redefining the objective function by an assisted history matching process in which a new parameter uncertainty domain is sampled for each of the geo-model, the discrete fracture network model, and the scenario of a simulation model, and generating a history matched model based on the redefined objective function; and redefining the objective function by an assisted history matching process in which a new parameter uncertainty domain is sampled for each of the geo-model, the discrete fracture network model, and the scenario of a simulation model, and generating a history matched model based on the redefined objective function; and producing a 3D model based on the forward simulation. <EOS>