СИСТЕМА И СПОСОБ ВЫПОЛНЕНИЯ МОДЕЛИРОВАНИЯ ПЛАСТОВОГО РЕЗЕРВУАРА С ИСПОЛЬЗОВАНИЕМ ПРЕДОБУСЛОВЛИВАНИЯ

1. Способ решения линейных систем уравнений, возникающих при моделировании пластового резервуара, при этом способ содержит этапы, на которых:
(a) предоставляют линейную систему уравнений потока текучей среды, которая представляет свойства потока текучей среды в геологических породах подземного пластового резервуара, причем линейная система уравнений потока текучей среды содержит множество неизвестных переменных;
(b) применяют блочно-диагональное масштабирование к линейной системе уравнений потока текучей среды, чтобы получить блочно-диагональную матрицу;
(c) применяют неполную нижнюю-верхнюю приближенную факторизацию с пороговым значением к блочно-диагональной матрице, чтобы получить матрицу предобусловливания; и
(d) итеративно решают линейную систему уравнений потока текучей среды для неизвестных переменных с использованием матрицы предобусловливания.
2. Система для решения линейных систем уравнений, возникающих при моделировании пластового резервуара, при этом система содержит:
- базу данных, выполненную с возможностью сохранять данные, содержащие линейную систему уравнений потока текучей среды, которая представляет свойства потока текучей среды в геологических породах подземного пластового резервуара, причем линейная система уравнений потока текучей среды содержит множество неизвестных переменных;
- процессор компьютера, выполненный с возможностью принимать сохраненные данные из базы данных и выполнять программные инструкции с использованием сохраненных данных; и
- программу, выполняемую на процессоре компьютера, причем программа содержит машиночитаемые программные инструкции, содержащие:
- модуль блочно-диагонального масштабирования, выполненный с возможностью применять блочно-диагональное масштабирование к линейной системе уравнений потока текучей среды, чтобы получить блочно-диагональную матрицу;
- модуль факторизации, выполненный с возможностью применять неполную нижнюю-верхнюю приближенную факторизацию с пороговым значением к блочно-диагональной матрице, чтобы получить матрицу предобусловливания; и
- модуль итеративного решающего устройства, выполненный с возможностью итеративно решать линейную систему уравнений потока текучей среды для неизвестных переменных с использованием матрицы предобусловливания.
3. Постоянный машиночитаемый носитель, содержащий машиночитаемые программные инструкции, используемые для решения линейных систем уравнений, возникающих при моделировании пластового резервуара, причем программные инструкции содержат:
- модуль блочно-диагонального масштабирования, выполненный с возможностью применять блочно-диагональное масштабирование к линейной системе уравнений потока текучей среды, чтобы получить блочно-диагональную матрицу;
- модуль факторизации, выполненный с возможностью применять
неполную нижнюю-верхнюю приближенную факторизацию с пороговым значением к блочно-диагональной матрице, чтобы получить матрицу предобусловливания; и
- модуль итеративного решающего устройства, выполненный с возможностью итеративно решать линейную систему уравнений потока текучей среды для неизвестных переменных с использованием матрицы предобусловливания.
4. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором члены в блочно-диагональной матрице обнуляются, если члены ниже предварительно определенного порогового значения.
5. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором уравновешивание по столбцам и строкам Синкорна-Кноппа или по строкам и столбцам Синкорна-Кноппа применяется до применения неполной нижней-верхней приближенной факторизации с пороговым значением к блочно-диагональной матрице.
6. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором уравновешивание по строкам и столбцам Руиза применяется до применения неполной нижней-верхней приближенной факторизации с пороговым значением к блочно-диагональной матрице.
7. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый
носитель по любому из пп. 1-3, в котором упорядочение по переменным на основе диагонального доминирования по строкам и столбцам применяется до применения неполной нижней-верхней приближенной факторизации с пороговым значением к блочно-диагональной матрице.
8. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором упорядочение по переменным на основе полного заполнения применяется до применения неполной нижней-верхней приближенной факторизации с пороговым значением к блочно-диагональной матрице.
9. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором упорядочение по ячейкам на основе обратного алгоритма Катхилла-Макки или "красно-черное" упорядочение по ячейкам применяется до применения неполной нижней-верхней приближенной факторизации с пороговым значением к блочно-диагональной матрице.
10. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором неполная нижняя-верхняя приближенная факторизация с пороговым значением применяется без выбора ведущего элемента матрицы.
11. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором матрица предобусловливания содержит строго-нестрогую неполную
факторизацию, по меньшей мере, с одним пороговым значением.
12. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором матрица предобусловливания используется непосредственно в способе на основе предобусловленного подпространства Крылова или используется для предобусловливания второй стадии в способе предобусловливания на основе ограниченного остатка давления.
13. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором:
- линейная система уравнений потока текучей среды представляется на этапе (d) следующим образом:
<mrow><mover accent="true"><mi>A</mi><mo>˜</mo></mover><mi>x</mi><mo>=</mo><mi>r</mi><mtext>  c  </mtext><mi>r</mi><mo>=</mo><mi>b</mi><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><msup><mrow><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mi>b</mi><mo>;</mo></mrow>
и
- матрица предобусловливания представляется следующим образом:
<mrow><msubsup><mi>M</mi><mrow><mi>C</mi><mi>P</mi><mi>R</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msub><mi>D</mi><mi>C</mi></msub><mi>P</mi><mtext> </mtext><msup><mi>U</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>L</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>P</mi><mi>T</mi></msup><msub><mi>D</mi><mi>R</mi></msub><mo stretchy="false">[</mo><mi>I</mi><mo>−</mo><mover accent="true"><mi>A</mi><mo>˜</mo></mover><mi>C</mi><mtext> </mtext><msubsup><mi>A</mi><mrow><mi>p</mi><mi>p</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mtext> </mtext><msup><mover accent="true"><mi>W</mi><mo>˜</mo></mover><mi>T</mi></msup><mo stretchy="false">]</mo><mo>+</mo><mi>C</mi><mtext> </mtext><msubsup><mi>A</mi><mrow><mi>p</mi><mi>p</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mtext> </mtext><msup><mover accent="true"><mi>W</mi><mo>˜</mo></mover><mi>T</mi></msup></mrow>
,
<mrow><msup><mover accent="true"><mi>W</mi><mo>˜</mo></mover><mi>T</mi></msup><mo>=</mo><msup><mi>W</mi><mi>T</mi></msup><mi>b</mi><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow>
где
.
14. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором:
- линейная система уравнений потока текучей среды представляется на этапе (d) следующим образом:
<mrow><mover accent="true"><mi>A</mi><mo>^</mo></mover><mtext> </mtext><mi>y</mi><mo>=</mo><mtext> </mtext><mi>r</mi><mtext>  c  </mtext><mi>r</mi><mo>=</mo><msub><mi>D</mi><mi>R</mi></msub><mtext> </mtext><mi>b</mi><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><msup><mrow><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mi>b</mi><mo>;</mo></mrow>
и
- матрица предобусловливания представляется следующим образом:
<mrow><msubsup><mi>M</mi><mrow><mi>C</mi><mi>P</mi><mi>R</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mi>P</mi><msup><mi>U</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>L</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>P</mi><mi>T</mi></msup><mo stretchy="false">[</mo><mi>I</mi><mo>−</mo><mover accent="true"><mi>A</mi><mo>^</mo></mover><mtext> </mtext><mover accent="true"><mi>C</mi><mo>^</mo></mover><mtext> </mtext><msubsup><mi>A</mi><mrow><mi>p</mi><mi>p</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mtext> </mtext><msup><mover accent="true"><mi>W</mi><mo>^</mo></mover><mi>T</mi></msup><mo stretchy="false">]</mo><mo>+</mo><mover accent="true"><mi>C</mi><mo>^</mo></mover><mtext> </mtext><msubsup><mi>A</mi><mrow><mi>p</mi><mi>p</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mtext> </mtext><msup><mover accent="true"><mi>W</mi><mo>^</mo></mover><mi>T</mi></msup></mrow>
,
<mrow><msup><mover accent="true"><mi>W</mi><mo>^</mo></mover><mi>T</mi></msup><mo>=</mo><msup><mi>W</mi><mi>T</mi></msup><mi>b</mi><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo><msubsup><mi>D</mi><mi>R</mi><mrow><mo>−</mo><mn>1</mn></mrow></msubsup></mrow>
где
.
15. Способ решения линейных систем уравнений, система для решения линейных систем уравнений, постоянный машиночитаемый носитель по любому из пп. 1-3, в котором решение линейной системы уравнений потока текучей среды используется для того, чтобы решить одну из задачи оптимизации добычи, задачи адаптации исторических данных, либо используется для того, чтобы регулировать параметр скважины в пластовом месторождении.