Help request with MRST

[Jun Li]

Dear MRST developers,

 

I am Jun Li from China University of Geosciences, and works on the research of reservoir simulations. Recently I have been learning MRST and hope it can be used in the future research work.

 

I set up a 3D simulation model of CO2 injection under the back-oil framework modified from the example file with a geological model (in Eclipse format, see the attached files). The input m-file can also be found from the attached.

If use default linear solver, the simulation is very slow. So I use AMG linear solver in the setup. During the injection period, the simulation is OK. But after shut in, it will crash as shown in the following image, even though very small time step is used.

 

Do you have any good ideas or comments to help the problem?

 

Thank you!

Kind Regards,

 

 Jun Li

---

[Knut-Andreas Lie]

Hi,

I have had a look at your problem, and it is fortunately is simple to solve. In MRST, it is not sufficient to just set the rate of the well to zero to turn in off. You must also set it explicitly as inactive. That is, add the following after line 163 in your script

      schedule.control(2).W.status = false;

so that the declaration of the schedule now reads

      schedule.control    = struct('W', W, 'bc', bc);

      schedule.control(2) = struct('W', W, 'bc', bc);

      schedule.control(2).W.val = 0;

      schedule.control(2).W.status = false;

At least, this seemed to do the trick on my computer and will hopefully also work on yours. 😉

Best regards,

Knut-Andreas

--

Knut-Andreas Lie, Professor, PhD
Chief Scientist, SINTEF Digital, Mathematics & Cybernetics, Oslo
Phone: +47 930 58 721 (mobile)
http://folk.ntnu.no/andreas

---

From: Jun Li <lij...@cug.edu.cn>
Sent: Friday, December 23, 2022 1:35 PM
To: MRST-mailliste <mr...@sintef.no>
Cc: sintef-mrst <sinte...@googlegroups.com>
Subject: Help request with MRST

 

[Jun Li]

Thank you so much, Knut-Andreas!

I will try.

Best regards,

Jun

---

 

From: Knut-Andreas Lie

Date: 2023-01-15 01:08

To: Jun Li; MRST-mailliste

CC: sintef-mrst

Subject: Re: Help request with MRST

Hi,

I have had a look at your problem, and it is fortunately is simple to solve. In MRST, it is not sufficient to just set the rate of the well to zero to turn in off. You must also set it explicitly as inactive. That is, add the following after line 163 in your script

      schedule.control(2).W.status = false;

so that the declaration of the schedule now reads

      schedule.control    = struct('W', W, 'bc', bc);

      schedule.control(2) = struct('W', W, 'bc', bc);

      schedule.control(2).W.val = 0;

      schedule.control(2).W.status = false;

At least, this seemed to do the trick on my computer and will hopefully also work on yours. ">

Best regards,

Knut-Andreas

--

Knut-Andreas Lie, Professor, PhD
Chief Scientist, SINTEF Digital, Mathematics & Cybernetics, Oslo
Phone: +47 930 58 721 (mobile)
http://folk.ntnu.no/andreas

---

From: Jun Li <lij...@cug.edu.cn>
Sent: Friday, December 23, 2022 1:35 PM
To: MRST-mailliste <mr...@sintef.no>
Cc: sintef-mrst <sinte...@googlegroups.com>
Subject: Help request with MRST