Hi,
I want to simulate a pumping well above seawater intrusion in SEAWAT using flopy and want to cease pumping automaticly when the salinity concentration in the well's cell reaches to a certain level, for example, 5% of relative salinity. In other words, I want the well only extracts freshwater and when the salinity starts increasing in the well (due to saltwater up-coning), stop pumping. I really appreciate it if someone helps me to do this task.
Pumping well can be simulated by utilizing optimization models are based on the sharp interface approximation and the Ghyben Herzberg relation.
Hi Amir,
As far as I'm aware of, the packages in SEAWAT do not allow you to do this internally.
A workaround with Flopy would be:
If you need this in code: if you can send me a script to create a simple test model, I can help you out with converting the above workflow to a script. It is then probably the most convenient to exchange scripts on Github instead of on Researchgate.
FYI: The USGS and Deltares are currently working on including the Basic Model Interface (BMI) https://csdms.colorado.edu/wiki/BMI_Description
for Modflow6, which allows users to directly change states in MODFLOW6 at runtime and allow them to create these kind of decision schemes in a pythonscript, to circumvent having to continuously write new input with Flopy.
This is still under development though, so I would not wait for it at the moment to fix this specific problem you are solving, but it could be interesting for you in the future.
Good luck!
Joeri
Hi Joeri,
Thanks so much for your answer and suggested method. I was looking for a way to modify the Seawat's Fortran subrutines and compile modified version of Seawat.
Cheers,
Amir
- Badges
- Science topic
- Similar topics
- Phytochemistry
- Extracts