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...
Saved in:
| Main Author: | |
|---|---|
| Other Authors: | , |
| Published: |
2025
|
| Subjects: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1864513521259118592 |
|---|---|
| author | Ali Asif (23717367) |
| author2 | Abdul Salam Abd (14151648) Ahmad Abushaikha (17148349) |
| author2_role | author author |
| author_facet | Ali Asif (23717367) Abdul Salam Abd (14151648) Ahmad Abushaikha (17148349) |
| author_role | author |
| dc.creator.none.fl_str_mv | Ali Asif (23717367) Abdul Salam Abd (14151648) Ahmad Abushaikha (17148349) |
| dc.date.none.fl_str_mv | 2025-12-27T09:00:00Z |
| dc.identifier.none.fl_str_mv | 10.1007/s10596-025-10391-5 |
| dc.relation.none.fl_str_mv | https://figshare.com/articles/journal_contribution/A_residual-accelerated_Jacobian_method_for_rapid_convergence_in_reservoir_simulation/32075541 |
| dc.rights.none.fl_str_mv | CC BY 4.0 info:eu-repo/semantics/openAccess |
| dc.subject.none.fl_str_mv | Engineering Resources engineering and extractive metallurgy Mathematical sciences Applied mathematics Linearization Finite difference Numerical methods Reservoir simulation |
| dc.title.none.fl_str_mv | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| dc.type.none.fl_str_mv | Text Journal contribution info:eu-repo/semantics/publishedVersion text contribution to journal |
| description | <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> |
| eu_rights_str_mv | openAccess |
| id | Manara2_15946c4158a00765c4c03f72074b54bb |
| identifier_str_mv | 10.1007/s10596-025-10391-5 |
| network_acronym_str | Manara2 |
| network_name_str | Manara2 |
| oai_identifier_str | oai:figshare.com:article/32075541 |
| publishDate | 2025 |
| repository.mail.fl_str_mv | |
| repository.name.fl_str_mv | |
| repository_id_str | |
| rights_invalid_str_mv | CC BY 4.0 |
| spelling | A residual-accelerated Jacobian method for rapid convergence in reservoir simulationAli Asif (23717367)Abdul Salam Abd (14151648)Ahmad Abushaikha (17148349)EngineeringResources engineering and extractive metallurgyMathematical sciencesApplied mathematicsLinearizationFinite differenceNumerical methodsReservoir 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 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>2025-12-27T09:00:00ZTextJournal contributioninfo:eu-repo/semantics/publishedVersiontextcontribution to journal10.1007/s10596-025-10391-5https://figshare.com/articles/journal_contribution/A_residual-accelerated_Jacobian_method_for_rapid_convergence_in_reservoir_simulation/32075541CC BY 4.0info:eu-repo/semantics/openAccessoai:figshare.com:article/320755412025-12-27T09:00:00Z |
| spellingShingle | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation Ali Asif (23717367) Engineering Resources engineering and extractive metallurgy Mathematical sciences Applied mathematics Linearization Finite difference Numerical methods Reservoir simulation |
| status_str | publishedVersion |
| title | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| title_full | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| title_fullStr | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| title_full_unstemmed | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| title_short | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| title_sort | A residual-accelerated Jacobian method for rapid convergence in reservoir simulation |
| topic | Engineering Resources engineering and extractive metallurgy Mathematical sciences Applied mathematics Linearization Finite difference Numerical methods Reservoir simulation |