A residual-accelerated Jacobian method for rapid convergence in reservoir simulation
<p dir="ltr">Reservoir simulation requires efficient algorithms for complex, nonlinear systems. We introduce a Residual–Accelerated Jacobian (RAJ) for fully implicit reservoir simulation. RAJ assembles the Jacobian by finite differences with a residual–adaptive step that enlarges far...
محفوظ في:
| المؤلف الرئيسي: | |
|---|---|
| مؤلفون آخرون: | , |
| منشور في: |
2025
|
| الموضوعات: | |
| الوسوم: |
إضافة وسم
لا توجد وسوم, كن أول من يضع وسما على هذه التسجيلة!
|
| الملخص: | <p dir="ltr">Reservoir simulation requires efficient algorithms for complex, nonlinear systems. We introduce a Residual–Accelerated Jacobian (RAJ) for fully implicit reservoir simulation. RAJ assembles the Jacobian by finite differences with a residual–adaptive step that enlarges far from convergence and shrinks near tolerance; saturation probes are projected to physical bounds. The governing residual is unchanged, so accuracy is preserved while Jacobian columns remain stable across nonlinear episodes. This adaptive mechanism ensures that the RAJ method remains responsive and can effectively adjust to the changing dynamics of the simulation, enabling faster convergence without compromising the accuracy of the solution, a crucial aspect in handling complex, non-linear systems in reservoir simulations. We evaluate RAJ on SPE10 (top five layers) and the Norne field (corner-point with NNCs), reporting CPU time, Newton/linear iterations, and robustness indicators (wasted steps/iterations), alongside accuracy parity. On SPE10 water injection, RAJ runs in <b>226.73</b> s vs FD <b>254.42</b> s (–<b>10.9%</b>); all methods use 43 Newton iterations and 2215 linear iterations, with no wasted steps. On PE10 gas injection, RAJ completes in <b>1600.55</b> s vs FD <b>1772.12</b> s (<b>–9.7%</b>) and lowers wasted work (wasted time-step fraction 31.98% vs 36.38%; wasted linear iterations 5385 vs 6644). On Norne, RAJ takes<b> 276.81</b> s vs FD <b>313.85 </b>s (<b>–11.8%</b>) with 51 vs 53 Newton iterations and zero wasted steps. As expected, analytical derivatives are fastest where available (water <b>154.816</b> s, gas <b>1177.69</b> s, Norne <b>223.935</b> s). Overall, RAJ delivers comparable or better CPU times than fixed-step FD while preserving accuracy and reducing wasted work, offering a practical, drop-in alternative when analytical Jacobians are unavailable.</p><h2 dir="ltr">Other Information</h2><p dir="ltr">Published in: Computational Geosciences<br>License: <a href="https://creativecommons.org/licenses/by/4.0" target="_blank">https://creativecommons.org/licenses/by/4.0</a><br>See article on publisher's website: <a href="https://dx.doi.org/10.1007/s10596-025-10391-5" target="_blank">https://dx.doi.org/10.1007/s10596-025-10391-5</a></p> |
|---|