We are trying to set up GWB simulation in which we simulate an injection through a series of injection wells within a 50x50 domain where the chemistry of injected water changes periodically. Since GWB only allows for two injection fluids in one run, we are intending to run this problem using multiple simulations in which the end result of one simulation will be the initial conditions of the next simulation and so on. We are compiling the new input file using FORTRAN to create the node by node input. We have successfully created and run subsequent injection periods, but we have run into some challenges.
The problem is at the end of the first simulation, we have different water types (pH, and equilibrium minerals) within different parts of the domain. Thus, “minerals in the system” (as listed in the output file at the end of the simulation) in one area of the domain around injection wells are much different than those in other areas of the domain. In other words, there is NO one mineral that we can set that constrains species concentrations. We can set species initial concentrations on node-by-node basis (from the output file of the previous simulation) but we don’t know how, if possible, to also define total moles of the existing “minerals in the system” on node-by-node basis in addition to species concentrations.
GWB 9.0 64-bit
Thanks!
Silver_lobo
