Jump to content
Geochemist's Workbench Support Forum

ASchu2010

Members
  • Posts

    4
  • Joined

  • Last visited

ASchu2010's Achievements

Newbie

Newbie (1/14)

0

Reputation

  1. Tom, Great, yes I am familiar with those two sources and they have indeed proved quite helpful! Also, thank you for looking at my script again. The system is running great now even with all of the minerals and water data included but I won't be able to make any quantitative results statements without solving this time issue and appreciate very much your expertise and suggestions. Thanks!! Anne
  2. Tom- Your comments have been immensely helpful! After revisiting the rate constants for slower rates and entering them, the system ran fine with all of the minerals. One of the critical questions I am trying to answer with respect to this project is how long the system takes to buffer the CO2 after injection has stopped. After running the system successfully several times and looking at plots in GtPlot I noticed that despite changing the start/end times the pH, for example, remained in the same pattern. So, running a system starting at 0 and ending at 2 days gave a pH of 5.2 at 0.5 days. Changing the start/end times to 0-30 days and looking at the pH vs. time chart again you would assume that it would still show 5.2 at 0.5 days, but instead the chart showed the pH decrease stretched out over the whole 30 days. I looked in the reaction modeling guide and input the kinetic dissolution of albite example on page 36/37. I then plotted the time vs. Q/K value just like the example showed on page 38. The plot was for 1500 years and so I wondered if I shortened the time frame as I did with my own system if GWB would shrink the profile to fit or if it would show a cut off version of what was shown for the 1500 year simulation. GWB didn't shrink the profile, it showed a cut off version which is what I would like to happen with my own system. Am I typing something in wrong or maybe missing something? Also, I was looking in the Short Course Notebook from Denver last year and on page 96 there is a table with rate constants. Do you know what the reference is for the values in the table? I have attached my updated script which now runs successfully using slower rate constants . Thank you so much for your help. -Anne Reynolds8_Avg_KineticMins4.rea
  3. Hi Tom, Thanks for such a prompt response! 1. I spent some time reading over the Kinetic Reactions section in the modeling guide today. I did notice the mention of the step size which corresponds to the "delxi" variable. I tried increasing this value by small incriments as well as by orders of magnitude to see if I could get the "Xi step size too small" error to disappear but my efforts were unsuccessful. Am I missing something? 2. I'm assuming by "initial speciation model" you mean the original values I used in SpecE8? If so, those values are the formation water data and converge fine in SpecE8. When I do a batch run with the water and minerals in REACT the system converges successfully as long as there are no clays present (if there are clays the system error reads: Didn't wake up, abandoning path). I tried your suggestion of titrating the minerals in one at a time but found none of them to allow the system to run. 3. I fixed my rate constants because, as you mentioned, they were negative and I realized I had log k rate constants. I converted them to the appropriate values which all turned out to be around +0.99. Despite this change, the system still showed the "Xi step size too small" error. With this change to non log values am I accounting for the (1-Q/K) portion of the rate constant equation that you refered to? 4. Additionally, I have heard from another GWB user that starting with a water/rock ratio of 1 can sometimes help, so I recalculated my rock masses to reflect this W/R=1. Even if I can get these kinetic reactions to converge and be successfull will me lowering the W/R ratio back down to more realistic reservoir ratios affect the reactions again? I have attached my edited script with the new rate constants. Reynolds8_Avg_KineticMins3.rea
  4. I am working on a carbon sequestration project and have both formation water and reservoir mineralogy data to use as inputs. I have calculated the surface area and looked up rate constants to input the minerals as kinetic minerals, but when I try to get my model to run I am getting an error message that says: --Can't converge, abandoning path --Xi step is too small I tried looking some on the forum to see if anyone else mentioned this Xi step but couldn't find anything. I have attached my script so you don't have to enter the data manually. Any help on why this error might be occuring would be helpful! Thanks! -Anne Reynolds8_Avg_KineticMins2.rea
×
×
  • Create New...