СИСТЕМЫ И СПОСОБЫ ДЛЯ ОЦЕНКИ МОМЕНТОВ ПРОРЫВА ФЛЮИДА В МЕСТОНАХОЖДЕНИЯХ ДОБЫВАЮЩИХ СКВАЖИН

1. Способ оценки момента прорыва флюида в добывающей скважине на основании данных моделирования распространения флюида, содержащий:
идентификацию данных отслеживания линии тока;
вычисление среднего времени пробега линии тока в каждой ячейке сети, на основании данных отслеживания линии тока;
идентификацию кратчайшей или быстрейшей линии тока для добывающей скважины, используя среднее время пробега линии тока в каждой ячейке сети;
вычисление среднего времени пролета для кратчайшей или быстрейшей линии тока через каждую пересекаемую ячейку сети, используя процессор вычислительной машины; и
оценку момента прорыва флюида в добывающей скважине, используя данные моделирования распространения флюида и среднее время пролета для кратчайшей или быстрейшей линии тока.
2. Способ по п. 1, в котором данные моделирования распространения флюида содержат момент вторжения флюида, представленный числом повторений моделирования, необходимых для достижения флюидом добывающей скважины от нагнетательной скважина через одну или более ячеек сети, представляющих собой модель свойств пласта.
3. Способ по п. 1, в котором данные отслеживания линии тока содержат число сегментов линии тока, пересекающих каждую ячейку сети, время пробега для каждого сегмента линии тока в каждой ячейке сети, индексы для каждой ячейки сети и общее число ячеек сети, пересекаемых всеми линиями тока, соединяющими нагнетательную скважину с добывающей скважиной.
4. Способ по п. 3, в котором среднее время пробега линии тока в каждой ячейке сети вычисляют по формуле:
<mrow><mo>∂</mo><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow></msub></mrow></munderover><mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow></mrow></mstyle></mrow>
(2)
SLN
<mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow></mrow>
где (N) является числом сегментов линий тока, пересекающих каждую ячейку сети, а
является временем пробега для каждого сегмента линий тока в каждой ячейке сети.
5. Способ по п. 1, в котором кратчайшая или быстрейшая линия тока для добывающей скважины представляет собой линию тока с наименьшей суммой средних времен пробега линии тока в ячейках сети, пересекаемых линией тока между нагнетательной скважиной и добывающей скважиной.
6. Способ по п. 5, в котором среднее время пролета для кратчайшей или быстрейшей линии тока вычисляют по каждой пересекаемой ячейке сети, используя наименьшую сумму средних времен пробега линии тока для кратчайшей или быстрейшей линии тока, общему числу ячеек сети, пересекаемых кратчайшей или быстрейшей линией тока.
7. Способ по п. 6, в котором среднее время пролета для кратчайшей или быстрейшей линии тока вычисляют по формуле:
<mrow><msup><mrow><mrow><mo>〈</mo><mrow><mi>T</mi><mi>O</mi><mi>F</mi></mrow><mo>〉</mo></mrow></mrow><mrow><mi>min</mi></mrow></msup><mo>=</mo><mfrac><mn>1</mn><mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow></munderover><mrow><mo>∂</mo><msubsup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mi>u</mi><mrow><mi>min</mi></mrow></msubsup></mrow></mstyle></mrow>
<mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow>
<mrow><mo>∂</mo><msubsup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mi>u</mi><mrow><mi>min</mi></mrow></msubsup></mrow>
где (
) представляет собой общее число всех ячеек сети, пересекаемых кратчайшей или быстрейшей линией тока, (
) представляет собой наименьшую сумму средних времен пробега для кратчайшей или быстрейшей линии тока, а (u) представляет собой число пробегов по всем индексам ячеек сети, пересекаемых кратчайшей или быстрейшей линей тока.
8. Способ по п. 2, в котором момент прорыва флюида в добывающей скважине оценивают по формуле:
<mrow><msup><mi>T</mi><mrow><mi>B</mi><mi>T</mi></mrow></msup><mo>=</mo><msup><mrow><mrow><mo>〈</mo><mrow><mi>T</mi><mi>O</mi><mi>F</mi></mrow><mo>〉</mo></mrow></mrow><mrow><mi>min</mi></mrow></msup><mo>×</mo><mfrac><mrow><msubsup><mi>t</mi><mrow><mi>I</mi><mi>N</mi><mi>V</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msubsup></mrow><mrow><msub><mi>N</mi><mi>p</mi></msub></mrow></mfrac><mo>×</mo><mfrac><mrow><msubsup><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup></mrow><mrow><msub><mi>N</mi><mrow><mi>x</mi><mi>y</mi><mi>z</mi></mrow></msub></mrow></mfrac></mrow>
xyz
p
min
<mrow><msubsup><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup></mrow>
<mrow><msubsup><mi>t</mi><mrow><mi>I</mi><mi>N</mi><mi>V</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msubsup></mrow>
где (N) и (N) представляют собой общий размер модели свойств пласта и общее число добывающих скважин, соответственно, (<TOF>) представляет собой среднее время пролета для кратчайшей или быстрейшей линии тока, (
) представляют собой общее число ячеек сети, пересекаемых всеми линиями тока, соединяющими нагнетательную скважину с добывающей скважиной, а (
)представляют собой время вторжения флюида.
9. Способ по п. 1, дополнительно содержащий повторение этапов по п. 1 для каждой добывающей скважины.
10. Способ по п. 1, в котором модель свойств пласта является моделью проницаемости.
11. Устройство постоянного хранения, материально несущее исполняемые компьютером команды для оценки момента прорыва флюида в добывающей скважине на основании данных моделирования распространения флюида, команды исполняют для реализации:
идентификации данных отслеживания линии тока;
вычисления среднего времени пробега линии тока в каждой ячейке сети, на основании данных отслеживания линии тока;
идентификации кратчайшей или быстрейшей линии тока для добывающей скважины, используя среднее время пробега линии тока в каждой ячейке сети;
вычисления среднего времени пролета для кратчайшей или быстрейшей линии тока через каждую пересекаемую ячейку сети; и
оценки момента прорыва флюида в добывающей скважине, используя данные моделирования распространения флюида и среднее время пролета для кратчайшей или быстрейшей линии тока.
12. Устройство носителя программы по п. 11, в котором данные моделирования распространения флюида содержат момент вторжения флюида, представленный числом повторений моделирования, необходимых для достижения флюидом добывающей скважины от нагнетательной скважина через одну или более ячеек сети, представляющих собой модель свойств пласта.
13. Устройство носителя программы по п. 11, в котором данные отслеживания линии тока содержат число сегментов линии тока, пересекающих каждую ячейку сети, время пробега для каждого сегмента линии тока в каждой ячейке сети, индексы для каждой ячейки сети и общее число ячеек сети, пересекаемых всеми линиями тока, соединяющими нагнетательную скважину с добывающей скважиной.
14. Устройство носителя программы по п. 13, в котором среднее время пробега линии тока в каждой ячейке сети вычисляют по формуле:
<mrow><mo>∂</mo><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow></msub></mrow></munderover><mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow></mrow></mstyle></mrow>
(2)
SLN
<mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow></mrow>
где (N) является числом сегментов линий тока, пересекающих каждую ячейку сети, а
является временем пробега для каждого сегмента линий тока в каждой ячейке сети.
15. Устройство носителя программы по п. 11, в котором кратчайшая или быстрейшая линия тока для добывающей скважины представляет собой линию тока с наименьшей суммой средних времен пробега линии тока в ячейках сети, пересекаемых линией тока между нагнетательной скважиной и добывающей скважиной.
16. Устройство носителя программы по п. 15, в котором среднее время пролета для кратчайшей или быстрейшей линии тока вычисляют по каждой пересекаемой ячейке сети, используя наименьшую сумму средних времен пробега линии тока для кратчайшей или быстрейшей линии тока и общему числу ячеек сети, пересекаемых кратчайшей или быстрейшей линией тока.
17. Устройство носителя программы по п. 16, в котором среднее время пролета для кратчайшей или быстрейшей линии тока вычисляют по формуле:
<mrow><msup><mrow><mrow><mo>〈</mo><mrow><mi>T</mi><mi>O</mi><mi>F</mi></mrow><mo>〉</mo></mrow></mrow><mrow><mi>min</mi></mrow></msup><mo>=</mo><mfrac><mn>1</mn><mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow></munderover><mrow><mo>∂</mo><msubsup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mi>u</mi><mrow><mi>min</mi></mrow></msubsup></mrow></mstyle></mrow>
<mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow>
<mrow><mo>∂</mo><msubsup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mi>u</mi><mrow><mi>min</mi></mrow></msubsup></mrow>
где (
) представляет собой общее число всех ячеек сети, пересекаемых кратчайшей или быстрейшей линией тока, (
) представляет собой наименьшую сумму средних времен пробега для кратчайшей или быстрейшей линии тока, а (u) представляет собой число пробегов по всем индексам ячеек сети, пересекаемых кратчайшей или быстрейшей линией тока.
18. Устройство носителя программы по п. 12, в котором момент прорыва флюида в добывающей скважине оценивают по формуле:
<mrow><msup><mi>T</mi><mrow><mi>B</mi><mi>T</mi></mrow></msup><mo>=</mo><msup><mrow><mrow><mo>〈</mo><mrow><mi>T</mi><mi>O</mi><mi>F</mi></mrow><mo>〉</mo></mrow></mrow><mrow><mi>min</mi></mrow></msup><mo>×</mo><mfrac><mrow><msubsup><mi>t</mi><mrow><mi>I</mi><mi>N</mi><mi>V</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msubsup></mrow><mrow><msub><mi>N</mi><mi>p</mi></msub></mrow></mfrac><mo>×</mo><mfrac><mrow><msubsup><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup></mrow><mrow><msub><mi>N</mi><mrow><mi>x</mi><mi>y</mi><mi>z</mi></mrow></msub></mrow></mfrac></mrow>
xyz
p
min
<mrow><msubsup><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup></mrow>
<mrow><msubsup><mi>t</mi><mrow><mi>I</mi><mi>N</mi><mi>V</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msubsup></mrow>
где (N) и (N) представляют собой общий размер модели свойств пласта и общее число добывающих скважин, соответственно, (<TOF>) представляет собой среднее время пролета для кратчайшей или быстрейшей линии тока, (
) представляют собой общее число ячеек сети, пересекаемых всеми линиями тока, соединяющими нагнетательную скважину с добывающей скважиной, а (
)представляют собой время вторжения флюида.
19. Устройство носителя программы по п. 11, дополнительно содержащий повторение этапов по п. 1 для каждой добывающей скважины.
20. Устройство носителя программы по п. 11, в котором модель свойств пласта является моделью проницаемости.