This submission includes synthetic seismic modeling data for the Push-Pull project at Brady Hot Springs, NV. The synthetic seismic is all generated by finite-difference method regarding different fracture and rock properties.

In numerical reservoir simulation naturally fractured reservoirs are commonly divided into matrix and fracture systems. The high permeability fractures are usually entirely responsible for flow between blocks and flow to the wells. The flow in these fractures is modeled using Darcy`s law and its extension to multiphase flow by means of relative permeabilities. The influence and measurement of fracture relative permeability for two-phase flow in fractured porous media have not been studied extensively, and the few works presented in the literature are contradictory. Experimental and numerical work on two-phase flow in fractured porous media has been initiated. An apparatus for monitoring this type of flow was designed and constructed. It consists of an artificially fractured core inside an epoxy core holder, detailed pressure and effluent monitoring, saturation measurements by means of a CT-scanner and a computerized data acquisition system. The complete apparatus was assembled and tested at conditions similar to the conditions expected for the two-phase flow experiments. Fine grid simulations of the experimental setup-were performed in order to establish experimental conditions and to study the effects of several key variables. These variables include fracture relative permeability and fracture capillary pressure. The numerical computations show that the flow is dominated by capillary imbibition, and that fracture relative permeabilities have only a minor influence. High oil recoveries without water production are achieved due to effective water imbibition from the fracture to the matrix. When imbibition is absent, fracture relative permeabilities affect the flow behavior at early production times.

The PoroTomo team has completed inverse modeling of the three data sets (seismology, geodesy, and hydrology) individually, as described previously. The estimated values of the material properties are registered on a three-dimensional grid with a spacing of 25 meters between nodes. The material properties are listed an Excel file. Figures show planar slices in three sets: horizontal slices in a planes normal to the vertical Z axis (Z normal), vertical slices in planes perpendicular to the dominant strike of the fault system (X normal), and vertical slices in planes parallel to the dominant strike of the fault system (Y normal). The results agree on the following points. The material is unconsolidated and/or fractured, especially in the shallow layers. The structural trends follow the fault system in strike and dip. The geodetic measurements favor the hypothesis of thermal contraction. Temporal changes in pressure, subsidence rate, and seismic amplitude are associated with changes in pumping rates during the four stages of the deployment in 2016. The modeled hydraulic conductivity is high in fault damage zones. All the observations are consistent with the conceptual model: highly permeable conduits along faults channel fluids from shallow aquifers to the deep geothermal reservoir tapped by the production wells.

TOUGH (Transport Of Unsaturated Groundwater and Heat) is a continuum-scale simulator developed by Lawrence Berkeley National Laboratory. TOUGH2 is used to simulate the coupled transport of water, vapor, non-condensible gas, and heat in porous and fractured media. TOUGH consists of a suite of tools including TOUGH3, iTOUGH2, TOUGHREACT, TOUGH2, and TOUGH+