RU2011149638A - Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора - Google Patents

Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора Download PDF

Info

Publication number
RU2011149638A
RU2011149638A RU2011149638/02A RU2011149638A RU2011149638A RU 2011149638 A RU2011149638 A RU 2011149638A RU 2011149638/02 A RU2011149638/02 A RU 2011149638/02A RU 2011149638 A RU2011149638 A RU 2011149638A RU 2011149638 A RU2011149638 A RU 2011149638A
Authority
RU
Russia
Prior art keywords
drainage
grid
well
block
perforated
Prior art date
Application number
RU2011149638/02A
Other languages
English (en)
Other versions
RU2549147C2 (ru
Inventor
Юсуф МИДДИЯ
Али Х. ДОГРУ
Original Assignee
Сауди Арабиан Ойл Компани
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Сауди Арабиан Ойл Компани filed Critical Сауди Арабиан Ойл Компани
Publication of RU2011149638A publication Critical patent/RU2011149638A/ru
Application granted granted Critical
Publication of RU2549147C2 publication Critical patent/RU2549147C2/ru

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/06Measuring temperature or pressure
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Measuring Fluid Pressure (AREA)
  • Electrical Discharge Machining, Electrochemical Machining, And Combined Machining (AREA)

Abstract

1. Система для определения статических давлений в скважине для одной или нескольких скважин произвольной формы (119) в однородных и неоднородных коллекторах (101), система содержит компьютерный имитатор (30) коллектора, имеющий совокупность процессоров (32), устройство ввода (34), устройство вывода (34), память (40), содержащую первую, вторую, третью, четвертую и пятую базы данных, система дополнительно отличается тем, что содержит:первую базу данных (92), включающую в себя совокупности векторов потока флюида (107), вторую базу данных (94), включающую в себя совокупность границ дренирования (110), третью базу данных (96), включающую в себя совокупность эффективных объемов дренирования (111), четвертую базу данных (98), включающую в себя совокупность давлений блока сетки, и пятую базу данных (99), включающую в себя совокупность расчетных статических давлений в скважине,программный продукт (60) определения потока флюида для определения потока флюида между соседними блоками сетки (103), причем программный продукт (60) определения потока флюида хранится в памяти (40) компьютерного имитатора коллектора (30) и содержит инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютерному имитатору коллектора (30) осуществлять, для каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) каждой из одной или нескольких скважин (119), на каждой ньютоновой итерации имитатора коллектора на основе сетки конечных разностей, операции- вычисления совокупности векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103), и- сохра�

Claims (50)

1. Система для определения статических давлений в скважине для одной или нескольких скважин произвольной формы (119) в однородных и неоднородных коллекторах (101), система содержит компьютерный имитатор (30) коллектора, имеющий совокупность процессоров (32), устройство ввода (34), устройство вывода (34), память (40), содержащую первую, вторую, третью, четвертую и пятую базы данных, система дополнительно отличается тем, что содержит:
первую базу данных (92), включающую в себя совокупности векторов потока флюида (107), вторую базу данных (94), включающую в себя совокупность границ дренирования (110), третью базу данных (96), включающую в себя совокупность эффективных объемов дренирования (111), четвертую базу данных (98), включающую в себя совокупность давлений блока сетки, и пятую базу данных (99), включающую в себя совокупность расчетных статических давлений в скважине,
программный продукт (60) определения потока флюида для определения потока флюида между соседними блоками сетки (103), причем программный продукт (60) определения потока флюида хранится в памяти (40) компьютерного имитатора коллектора (30) и содержит инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютерному имитатору коллектора (30) осуществлять, для каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) каждой из одной или нескольких скважин (119), на каждой ньютоновой итерации имитатора коллектора на основе сетки конечных разностей, операции
- вычисления совокупности векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103), и
- сохранения совокупности векторов потока флюида (107) в первой базе данных (92),
программный продукт (70) определения границ дренирования для определения трехмерных границ эффективного объема дренирования (111) для каждой из одной или нескольких скважин (119), причем программный продукт (70) определения границ дренирования хранится в памяти (40) компьютерного имитатора коллектора (30) и содержит инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютерному имитатору коллектора (30) осуществлять, для каждой из одной или нескольких плоскостей дренирования (115) соответствующего одного или нескольких сегментов скважины (117), каждая из которых ориентирована перпендикулярно направлению одного или нескольких сегментов скважины (117) соответствующей одной или нескольких скважин (119) и включает в себя грань (105) перфорированного блока сетки (113) и грань (105) каждого из совокупности блоков сетки (103), соседствующего с ним, операции
- осуществления на каждом из совокупности блоков сетки (103), расположенных в каждом отдельном одном из совокупности направлений, задающих совокупность траекторий (120), выходящих наружу из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115)
- отслеживания знака каждого из совокупности расчетных векторов потока флюида (107) для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
- сравнения величины фазового потока вдоль двух перпендикулярных осей для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), с величиной соответствующего фазового потока на перфорированном блоке сетки (113) для соответствующей плоскости дренирования (115), и
- определения расстояния каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
- определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), для каждого из одного или нескольких сегментов скважины (117) для каждой из одной или нескольких скважин (119), для оценивания, таким образом, эффективного объема дренирования (111) каждой из одной или нескольких скважин (119), причем положение каждой границы дренирования (110) каждой отдельной плоскости дренирования (115) одной или нескольких плоскостей дренирования (115) является соответствующим одним из следующих положений, находящихся на меньшем расстоянии от соответствующего перфорированного блока сетки (113), чем все остальные из следующих положений:
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на соответствующем перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где расстояние соответствующего одного из совокупности блоков сетки (103) от соответствующего перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
- сохранения определенной границы дренирования (110) во второй базе данных (94), и
программный продукт (80) определения статического давления в скважине для определения статического давления в скважине для каждой из одной или нескольких скважин (119), причем программный продукт (80) определения статического давления в скважине хранится в памяти (40) компьютерного имитатора коллектора (30) и содержит инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютерному имитатору коллектора (30) осуществлять операции
- определения оценки эффективного объема дренирования (111) каждой из одной или нескольких скважин (119) в соответствии с соответствующими определенными границами дренирования (110) каждой из одной или нескольких плоскостей дренирования (115), связанных с соответствующей одной или несколькими скважинами (119),
- сохранения оценки эффективного объема дренирования (111) каждой из одной или нескольких скважин (119) в третьей базе данных (96),
- определения динамического давления блока сетки каждого из совокупности блоков сетки (103), содержащихся в соответствующем эффективном объеме дренирования (111) каждой из одной или нескольких скважин (119),
- сохранения динамических давлений блока сетки в четвертой базе данных (98),
- определения усредненного по объему перового пространства динамического давления блока сетки, по меньшей мере, существенного подмножества совокупности блоков сетки (103), содержащихся в каждом соответствующем эффективном объеме дренирования (111) каждой соответствующей скважины (119) для задания, таким образом, оценочного статического давления в скважине для каждой из одной или нескольких скважин (119), и
- сохранения оценочного статического давления в скважине для каждой отдельной одной из одной или нескольких скважин (119) в пятой базе данных (99).
2. Система по п.1, в которой операция определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115) для каждого из одного или нескольких сегментов скважины (117) для каждой из одной или нескольких скважин (119), для оценивания, таким образом, эффективного объема дренирования (111) каждой из одной или нескольких скважин (119), включает в себя:
аппроксимацию одной или нескольких из определенных границ дренирования (110) с использованием непрерывных прямых линий.
3. Система по п.1, в которой программный продукт (80) определения статического давления в скважине дополнительно включает в себя инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютерному имитатору коллектора (30) осуществлять операцию
упрощения формы эффективного объема дренирования (111) каждой из одной или нескольких скважин (119) на каждом временном шаге имитатора коллектора до формы прямоугольной трехмерной приближенной ячейки объема дренирования (111), содержащей совокупность слоев (118), причем, по меньшей мере, один из слоев (118) имеет другую толщину, чем, по меньшей мере, один другой из слоев (118).
4. Система по п.1, в которой модель одной или нескольких скважин (119), используемая программным продуктом (80) определения статического давления в скважине, содержит совокупность плоскостей дренирования (115), и эффективный объем дренирования (111) одной или нескольких скважин (119) оценивается на основании определенных границ дренирования (110) двух перпендикулярных плоскостей дренирования (115) из совокупности плоскостей дренирования (115) для снижения, таким образом, потребности в вычислительных ресурсах.
5. Система по любому из пп.1-4, в которой модель коллектора, используемая программным продуктом (70) определения границ дренирования, делится среди совокупности процессоров (32) для осуществления операции определения положения каждой отдельной границы дренирования (110) каждой соответствующей плоскости дренирования (115) для одной или нескольких скважин (119),
причем модель одной из одной или нескольких скважин (119), используемая программным продуктом (80) определения статического давления в скважине, включает в себя мультидоменную плоскость дренирования (115), имеющую границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119), и
программный продукт (70) определения границ дренирования дополнительно включает в себя инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютеру (30) осуществлять, когда внешний предел границы домена (203) первого процессора (32) достигается при прохождении мультидоменной плоскости дренирования (115) вдоль одной из совокупности траекторий (120), выходящих из перфорированного блока сетки (113) мультидоменной плоскости дренирования (115), операции
- сохранения данных расстояния для пройденного блока сетки (103), расположенного вдоль одной из совокупности траекторий (120) на внешнем пределе границы домена (203) первого процессора (32),
- сохранения данных направления, указывающих направление прохождения вдоль одной из совокупности траекторий (120), на которой достигается соответствующий внешний предел границы домена (203) первого процессора (32), и
- синхронизации данных расстояния и направления со вторым из совокупности процессоров (32), имеющим границу домена (203) рядом с внешним пределом границы домена (203) первого процессора (32).
6. Система по любому из пп.1-4, в которой модель коллектора, используемая программным продуктом (70) определения границ дренирования, делится среди совокупности процессоров (32) для осуществления операции определения положения каждой отдельной границы дренирования (110) каждой соответствующей плоскости дренирования (115) для одной или нескольких скважин (119),
причем модель одной из одной или нескольких скважин (119), используемая программным продуктом (80) определения статического давления в скважине, включает в себя мультидоменную плоскость дренирования (115), имеющую границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119), в, по меньшей мере, одном из совокупности направлений, и
программный продукт (70) определения границ дренирования дополнительно включает в себя инструкции, которые, при выполнении компьютерным имитатором коллектора (30), предписывают компьютеру (30) осуществлять операции формирования для каждого из совокупности направлений, в которых соответствующая граница дренирования (110) выходит за пределы границы домена (203) первого процессора (32), контура с длиной контура, равной количеству совокупности процессоров (32).
7. Система по любому из пп.1-4, в которой оценочное статическое давление в скважине для каждой из одной или нескольких скважин (119) равно оценочному среднему статическому давлению в скважине (p_static) соответствующего эффективного объема дренирования (111), определенному согласно следующему вычислению:
p_static=[Sum(Vp(k))·Р(k)]Sum(Vp(k)),
где k относится ко всем блокам сетки (103) в эффективном объеме дренирования (111),
Vp (k) - объем перового пространства блока сетки (103) k в эффективном объеме дренирования (111),
P(k) - динамическое давление блока сетки (103) k, и
Sum - оператор суммирования.
8. Система по любому из пп.1-4, в которой операция определения оценки эффективного объема дренирования (111) каждой из одной или нескольких скважин (119) содержит определение оценки эффективного объема дренирования (111) каждой из одной или нескольких скважин (119) на каждой из совокупности ньютоновых итераций в полностью неявном имитаторе коллектора.
9. Система по любому из пп.1-4, в которой операция определения оценки эффективного объема дренирования (111) каждой из одной или нескольких скважин (119) осуществляется для неоднородного коллектора (101).
10. Система по любому из пп.1-4, в которой каждое из совокупности направлений от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115) для каждой из одной или нескольких плоскостей дренирования (115) соответствующего одного или нескольких сегментов скважины (117), по меньшей мере, одной из одной или нескольких скважин (119) отличается от любого другого на 90°.
11. Программный продукт (50) моделирования коллектора для определения статических давлений в скважине для одной или нескольких скважин произвольной формы (119) в однородных и неоднородных коллекторах (101), причем программный продукт (50) содержит набор инструкций, хранящихся на вещественном носителе, считываемом компьютером (30), которые, при выполнении компьютером (30), предписывают компьютеру (30) осуществлять различные операции, включающие в себя, для каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), вычисление совокупности векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103), операции дополнительно отличаются тем, что содержат:
выбор одной или нескольких плоскостей дренирования (115) соответствующего одного или нескольких сегментов (117) скважины (119), причем каждая из одной или нескольких плоскостей дренирования (115) ориентирована перпендикулярно направлению соответствующего одного или нескольких сегментов (117) скважины (119) и включает в себя грань (105) перфорированного блока сетки (113) и грань (105) каждого из совокупности блоков сетки (103), соседствующего с ним,
осуществление на каждом из совокупности блоков сетки (103), расположенных в каждом отдельном одном из совокупности направлений, задающих совокупность траекторий (120), выходящих наружу из соответствующего перфорированного блока сетки (113) на соответствующей одной из одной или нескольких плоскостей дренирования (115), для определения, таким образом, каждой границы дренирования (110) эффективного объема дренирования (111) скважины (119)
- отслеживания знака каждого из совокупности расчетных векторов потока флюида (107) для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
- сравнения величины фазового потока для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), с величиной соответствующего фазового потока на перфорированном блоке сетки (113) для соответствующей плоскости дренирования (115), и
- определения расстояния каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
определение положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115) для каждого из одного или нескольких сегментов (117) скважины (119), для оценивания, таким образом, эффективного объема дренирования (111) скважины (119), причем положение каждой границы дренирования (110) каждой отдельной плоскости дренирования (115) одной или нескольких плоскостей дренирования (115) является соответствующим одним из следующих положений, находящихся на меньшем расстоянии от соответствующего перфорированного блока сетки (113), чем все остальные из следующих положений:
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на соответствующем перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где расстояние соответствующего одного из совокупности блоков сетки (103) от соответствующего перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
определение оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110), по меньшей мере, одной из одной или нескольких плоскостей дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111),
определение динамического давления блока сетки каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), и
определение усредненного по объему перового пространства динамического давления блока сетки, по меньшей мере, существенного подмножества совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), для задания, таким образом, оценочного статического давления в скважине для скважины (119).
12. Программный продукт (50) моделирования коллектора по п.11, в котором операция определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115) для каждого из одного или нескольких сегментов (117) скважины (119)) для определения, таким образом, эффективного объема дренирования (111) скважины (119), включает в себя:
аппроксимацию одной или нескольких из определенных границ дренирования (110) с использованием непрерывных прямых линий.
13. Программный продукт (50) моделирования коллектора по п.11, в котором операции дополнительно отличаются тем, что
определение оценки эффективного объема дренирования (111) скважины (119), включает в себя операцию
упрощения формы эффективного объема дренирования (111) скважины (119) на каждом временном шаге имитатора коллектора до формы прямоугольной трехмерной ячейки объема дренирования (111), содержащей совокупность слоев (118), причем, по меньшей мере, один из слоев (118) имеет другую толщину, чем, по меньшей мере, один другой из слоев (118).
14. Программный продукт (50) моделирования коллектора по п.11, в котором модель скважины (119) содержит совокупность плоскостей дренирования (115), и эффективный объем дренирования (111) скважины (119) оценивается на основании определенных границ дренирования (110) двух перпендикулярных плоскостей дренирования (115) из совокупности плоскостей дренирования (115) для снижения, таким образом, потребности в вычислительных ресурсах.
15. Программный продукт (50) моделирования коллектора по любому из пп.11-14, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления операции определения положения каждой отдельной границы дренирования (110) каждой соответствующей плоскости дренирования (115) для скважины (119),
причем модель скважины (119) включает в себя мультидоменную плоскость дренирования (115), имеющую границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119), и
когда внешний предел границы домена (203) первого процессора (32) достигается при прохождении мультидоменной плоскости дренирования (115) вдоль одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113), осуществляются операции
- сохранения данных расстояния для пройденного блока сетки (103), расположенного вдоль одной из совокупности траекторий (120) на внешнем пределе границы домена (203) первого процессора (32),
- сохранения данных направления, указывающих направление прохождения вдоль одной из совокупности траекторий (120), на которой достигается соответствующий внешний предел границы домена (203) первого процессора (32), и
- синхронизации данных расстояния и направления со вторым из совокупности процессоров (32), имеющим границу домена (203) рядом с внешним пределом границы домена (203) первого процессора (32).
16. Программный продукт (50) моделирования коллектора по любому из пп.11-14, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления операции определения положения каждой отдельной границы дренирования (110) каждой соответствующей плоскости дренирования (115) для скважины (119),
причем модель скважины (119) включает в себя мультидоменную плоскость дренирования (115), имеющую границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119), в, по меньшей мере, одном из совокупности направлений, и
операции дополнительно отличаются тем, что содержат формирование для каждого из совокупности направлений, в которых соответствующая граница дренирования (110) выходит за пределы границы домена (203) первого процессора (32), контура с длиной контура, равной количеству совокупности процессоров (32).
17. Программный продукт (50) моделирования коллектора по любому из пп.11-14, в котором оценочное статическое давление в скважине для скважины (119) равно оценочному среднему статическому давлению в скважине эффективного объема дренирования (111) (p_static), определенному согласно следующему вычислению:
p_static=[Sum(Vp(k))·P(k)]Sum(Vp(k)),
где k относится ко всем блокам сетки (103) в эффективном объеме дренирования (111),
Vp(k) - объем перового пространства блока сетки (103) k в эффективном объеме дренирования (111),
P(k) - динамическое давление блока сетки (103) k, и
Sum - оператор суммирования.
18. Программный продукт (50) моделирования коллектора по любому из пп.11-14, в котором операция определения оценки эффективного объема дренирования (111) скважины (119) содержит:
определение оценки эффективного объема дренирования (111) скважины (119) на каждой из совокупности ньютоновых итераций в полностью неявном имитаторе коллектора.
19. Программный продукт (50) моделирования коллектора по любому из пп.11-14, в котором операция определения оценки эффективного объема дренирования (111) осуществляется для неоднородного коллектора (101).
20. Программный продукт (50) моделирования коллектора по любому из пп.11-14, в котором каждое из совокупности направлений от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115) для каждой из одной или нескольких плоскостей дренирования (115) соответствующего одного или нескольких сегментов скважины (117), по меньшей мере, одной из одной или нескольких скважин (119) отличается от любого другого на 90°.
21. Программный продукт (50) моделирования коллектора для определения статических давлений в скважине для одной или нескольких скважин произвольной формы (119) в однородных и неоднородных коллекторах (101), причем программный продукт содержит набор инструкций, хранящихся на вещественном носителе, считываемом компьютером (30), которые, при выполнении компьютером (30), предписывают компьютеру (30) осуществлять различные операции, отличающийся тем, что содержит:
определение положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из перфорированного блока сетки (113) на плоскости дренирования (115) сегмента скважины (117) скважины, для оценивания, таким образом, эффективного объема дренирования (111) скважины (119), причем
- для имитаций коллектора, где множественные скважины (119) состязаются за отдачу флюида, положение соответствующей границы дренирования (110) определяется как
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) из перфорированного блока сетки (113) на плоскости дренирования (115), и
- для имитаций коллектора, где множественные скважины (119) незначительно состязаются за отдачу флюида, в результате чего возможность впервые встретить изменение потока флюида в силу конкуренции скважин не предполагается, положение соответствующей границы дренирования (110) определяется как соответствующее одно из следующих положений, находящихся на меньшем расстоянии от перфорированного блока сетки (113), чем другие из следующих положений:
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) из перфорированного блока сетки (113) на плоскости дренирования (115), и
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где расстояние соответствующего одного из совокупности блоков сетки (103) из перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) из перфорированного блока сетки (113) на плоскости дренирования (115),
определение динамического давления блока сетки каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), и
определение усредненного по объему перового пространства динамического давления блока сетки, по меньшей мере, существенного подмножества совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), для задания, таким образом, оценочного статического давления в скважине для скважины (119).
22. Программный продукт (50) моделирования коллектора по п.21, в котором плоскость дренирования (115) является выбранной плоскостью дренирования (115) соответствующего сегмента скважины (117) скважины (119), которая ориентирована перпендикулярно направлению сегмента скважины (117) скважины (119) и включает в себя грань (105) перфорированного блока сетки (113) и грань (105) каждого из совокупности блоков сетки (103), соседствующего с ним, причем операции дополнительно отличаются тем, что содержат:
осуществление для каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), и на каждой ньютоновой итерации вычисления каждого временного шага имитатора коллектора на основе сетки конечных разностей, операции вычисления совокупности векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103), и
осуществление на каждом из совокупности блоков сетки (103), расположенных в каждом отдельном одном из совокупности направлений, задающих совокупность траекторий (120), выходящих наружу из перфорированного блока сетки (113) на плоскости дренирования (115), для определения, таким образом, каждой границы дренирования (110) эффективного объема дренирования (111) скважины (119)
- отслеживания знака каждого из совокупности расчетных векторов потока флюида (107) для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113), на плоскости дренирования (115),
- сравнения величины фазового потока для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113), с величиной соответствующего фазового потока на перфорированном блоке сетки (113) для плоскости дренирования (115), и
- определения расстояния каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113), на плоскости дренирования (115).
23. Программный продукт (50) моделирования коллектора по п.21, в котором операция определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из перфорированного блока сетки (113) на плоскости дренирования (115), включает в себя определение положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113), как для имитаций коллектора, где множественные скважины (119) состязаются за отдачу флюида, так и для имитаций коллектора, где множественные скважины (119) не состязаются за отдачу флюида, как одного из следующих положении, находящихся на меньшем расстоянии от перфорированного блока сетки (113), чем все остальные из следующих положений:
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) из перфорированного блока сетки (113),
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на соответствующем перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) из перфорированного блока сетки (113), и
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где расстояние соответствующего одного из совокупности блоков сетки (103) из перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) из перфорированного блока сетки (113).
24. Программный продукт (50) моделирования коллектора по п.21, в котором операция определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из перфорированного блока сетки (113) на плоскости дренирования (115) сегмента скважины (117) скважины, для оценивания, таким образом, эффективного объема дренирования (111) скважины (119), включает в себя:
аппроксимацию одной или нескольких из определенных границ дренирования (110) с использованием непрерывных прямых линий.
25. Программный продукт (50) моделирования коллектора по любому из пп.21-24 дополнительно отличается тем, что операция
определения оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) плоскости дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111), включает в себя операцию
упрощения формы объема дренирования (111) скважины (119) на каждом временном шаге имитатора коллектора до формы прямоугольной трехмерной ячейки объема дренирования (111), содержащей совокупность слоев (118), причем, по меньшей мере, один из слоев (118) имеет другую толщину, чем, по меньшей мере, один другой из слоев (118).
26. Программный продукт (50) моделирования коллектора по любому из пп.21-24, в котором модель скважины (119) содержит совокупность плоскостей дренирования (115), и эффективный объем дренирования (111) скважины (119) оценивается на основании определенных границ дренирования (110) двух перпендикулярных плоскостей дренирования (115) из совокупности плоскостей дренирования (115).
27. Программный продукт (50) моделирования коллектора по любому из пп.21-24, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления операции определения положения каждой отдельной границы дренирования (110) плоскости дренирования (115) сегмента скважины (117) для скважины (119),
причем плоскость дренирования (115) является мультидоменной плоскостью дренирования (115), имеющей границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего скважину (119), и
когда внешний предел границы домена (203) первого процессора (32) достигается при прохождении мультидоменной плоскости дренирования (115) вдоль одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113), осуществляются операции
- сохранения данных расстояния для пройденного блока сетки (103), расположенного вдоль одной из совокупности траекторий (120) на внешнем пределе границы домена (203) первого процессора (32),
- сохранения данных направления, указывающих направление прохождения вдоль одной из совокупности траекторий (120), на которой достигается соответствующий внешний предел границы домена (203) первого процессора (32), и
- синхронизации данных расстояния и направления со вторым из совокупности процессоров (32), имеющим границу домена (203) рядом с внешним пределом границы домена (203) первого процессора (32).
28. Программный продукт (50) моделирования коллектора по любому из пп.21-24, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления операции определения положения каждой отдельной границы дренирования (110) плоскости дренирования (115) для скважины (119),
причем плоскость дренирования (115) является мультидоменной плоскостью дренирования (115), имеющей границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119) в, по меньшей мере, одном из совокупности направлений, и
операции дополнительно отличаются тем, что содержат формирование для каждого из совокупности направлений, в которых соответствующая граница дренирования (110) выходит за пределы границы домена (203) первого процессора (32), контура с длиной контура, равной количеству совокупности процессоров (32).
29. Программный продукт (50) моделирования коллектора по любому из пп.21-24, операции дополнительно отличаются тем, что содержат:
определение оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) плоскости дренирования (115) для идентификации, таким образом, блоков сетки (103), содержащихся в объеме дренирования (111),
причем операция определения эффективного объема дренирования (111) осуществляется для неоднородного коллектора (101) с применением численного решения, и
оценочное статическое давление в скважине для скважины (119) равно оценочному среднему статическому давлению в скважине (р static) эффективного объема дренирования (111), определенному согласно следующему вычислению:
p_static=[Sum(Vp(k))·P(k)]Sum(Vp(k)),
где k относится ко всем блокам сетки (103) в эффективном объеме дренирования (111),
Vp(k) - объем перового пространства блока сетки (103) k в эффективном объеме дренирования (111),
P(k) - динамическое давление блока сетки (103) k, и
Sum - оператор суммирования.
30. Программный продукт (50) моделирования коллектора по любому из пп.21-24, операции дополнительно отличаются тем, что
определение оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) плоскости дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111), включает в себя операцию
определения оценки эффективного объема дренирования (111) скважины (119) на каждой из совокупности ньютоновых итераций в полностью неявном имитаторе коллектора.
31. Компьютерно-реализуемый способ определения совокупности приближенных статических давлений в скважине для одной или нескольких скважин произвольной формы (119) в коллекторе, компьютерно-реализуемый способ содержит этап, на котором вычисляют совокупность векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины, причем способ дополнительно отличается тем, что содержит этапы, на которых:
осуществляют этап вычисления совокупности векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103) для каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины, и на каждой ньютоновой итерации вычисления каждого временного шага имитатора коллектора на основе сетки конечных разностей,
выбирают одну или несколько плоскостей дренирования (115) соответствующего одного или нескольких сегментов скважины (117) для скважины (119), причем каждая плоскость дренирования (115) включает в себя перфорированный блок сетки (113) и перпендикулярна к соответствующему сегменту скважины (117) скважины (119),
определяют положение каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на каждой из одной или нескольких плоскостей дренирования (115), для оценивания, таким образом, эффективного объема дренирования (111) каждой скважины (119), причем
- для имитаций коллектора, где множественные скважины (119) состязаются за отдачу флюида, в результате чего возможность впервые встретить изменение потока флюида не предполагается, определяют положение соответствующей границы дренирования (110) как
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
- для имитаций коллектора, где множественные скважины (119) незначительно состязаются за отдачу флюида, определяют положение соответствующей границы дренирования (110) как соответствующее одно из следующих положений, находящихся на меньшем расстоянии от соответствующего перфорированного блока сетки (113), чем другие из следующих положений:
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на соответствующем перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где расстояние соответствующего одного из совокупности блоков сетки (103) от соответствующего перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
определяют оценку эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) каждой из одной или нескольких плоскостей дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111),
определяют давление блока сетки каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), и
определяют усредненное по объему перового пространства динамическое давление блока сетки, по меньшей мере, существенного подмножества совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119) для задания, таким образом, оценочного статического давления в скважине для скважины (119).
32. Компьютерно-реализуемый способ по п.31, отличающийся тем, что на каждом из совокупности блоков сетки (103), расположенных на каждой отдельной одной из совокупности траекторий (120), выходящих наружу из соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), для определения, таким образом, каждой границы дренирования (110) эффективного объема дренирования (111) скважины (119), дополнительно осуществляют этапы, на которых:
отслеживают знак каждого из совокупности расчетных векторов потока флюида (107) для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
сравнивают величину фазового потока для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), с величиной соответствующего фазового потока на перфорированном блоке сетки (113) для плоскости дренирования (115), и
определяют расстояние каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113) на соответствующей плоскости дренирования (115).
33. Компьютерно-реализуемый способ по п.31, в котором на этапе определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на каждой из одной или нескольких плоскостей дренирования (115), определяют положение каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на каждой из одной или нескольких плоскостей дренирования (115), как для имитаций коллектора, где множественные скважины (119) состязаются за отдачу флюида, так и для имитаций коллектора, где множественные скважины (119) не состязаются за отдачу флюида, как одно из следующих положений, находящихся на меньшем расстоянии от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), чем все остальные из следующих положений:
положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115),
положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на соответствующем перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), и
положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115), где расстояние соответствующего одного из совокупности блоков сетки (103) от соответствующего перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) от соответствующего перфорированного блока сетки (113) на соответствующей плоскости дренирования (115).
34. Компьютерно-реализуемый способ по любому из пп.31-33, в котором на этапе определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113) на каждой из одной или нескольких плоскостей дренирования (115), для оценивания, таким образом, эффективного объема дренирования (111) скважины (119)
аппроксимируют одну или несколько определенных границ дренирования (110) с использованием непрерывных прямых линий.
35. Компьютерно-реализуемый способ по любому из пп.31-33, в котором на этапе определения оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) одной или нескольких плоскостей дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111),
упрощают форму эффективного объема дренирования (111) скважины (119) на каждом временном шаге полностью неявного имитатора коллектора до формы прямоугольной трехмерной ячейки объема дренирования (111), содержащей совокупность слоев (118), причем, по меньшей мере, один из слоев (118) имеет другую толщину, чем, по меньшей мере, один другой из слоев (118).
36. Компьютерно-реализуемый способ по любому из пп.31-33, в котором модель скважины (119) содержит совокупность плоскостей дренирования (115), и эффективный объем дренирования (111) скважины (119) оценивается на основании определенных границ дренирования (110) двух перпендикулярных плоскостей дренирования (115) из совокупности плоскостей дренирования (115) для снижения, таким образом, потребности в вычислительных ресурсах.
37. Компьютерно-реализуемый способ по любому из пп.31-33, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления этапа определения положения каждой отдельной границы дренирования (110) соответствующей плоскости дренирования (115) сегмента скважины (117) для скважины (119)
причем, по меньшей мере, одна плоскость дренирования (115) из одной или нескольких плоскостей дренирования (115) является мультидоменной плоскостью дренирования (115), имеющей границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего скважину (119), и
когда внешний предел границы домена (203) первого процессора (32) достигается при прохождении мультидоменной плоскости дренирования (115) вдоль одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113), осуществляют этапы, на которых:
сохраняют данные расстояния для пройденного блока сетки (103), расположенного вдоль одной из совокупности траекторий (120) на внешнем пределе границы домена (203) первого процессора (32),
сохраняют данные направления, указывающие направление прохождения вдоль одной из совокупности траекторий (120), на которой достигается соответствующий внешний предел границы домена (203) первого процессора (32), и
синхронизируют данные расстояния и направления со вторым из совокупности процессоров (32), имеющим границу домена (203) рядом с внешним пределом границы домена (203) первого процессора (32).
38. Компьютерно-реализуемый способ по любому из пп.31-33, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления этапа определения положения каждой отдельной границы дренирования (110) соответствующей плоскости дренирования (115) для скважины (119),
причем, по меньшей мере, одна плоскость дренирования (115) из одной или нескольких плоскостей дренирования (115) является мультидоменной плоскостью дренирования (115), имеющей границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119) в, по меньшей мере, одном из совокупности направлений,
причем способ дополнительно отличается тем, что содержит этап, на котором формируют для каждого из совокупности направлений, в которых соответствующая граница дренирования (110) выходит за пределы границы домена (203) первого процессора (32), контур с длиной контура, равной количеству совокупности процессоров (32).
39. Компьютерно-реализуемый способ по любому из пп.31-33, дополнительно отличается тем, что содержит этап, на котором:
определяют оценку эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) каждой из одной или нескольких плоскостей дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111),
причем этап оценивания эффективного объема дренирования (111) осуществляют для неоднородного коллектора (101) с применением численного решения, и
оценочное статическое давление в скважине для скважины (119) равно оценочному среднему статическому давлению в скважине (p_static) эффективного объема дренирования (111), определенному согласно следующему вычислению:
p_static=[Sum(Vp(k))·P(k)]Sum(Vp(k)),
где k относится ко всем блокам сетки (103) в эффективном объеме дренирования (111),
Vp(k) - объем перового пространства блока сетки (103) k в эффективном объеме дренирования (111),
Р(k) - динамическое давление блока сетки (103) k, и
Sum - оператор суммирования.
40. Компьютерно-реализуемый способ по любому из пп.31-33, в котором на этапе определения оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) каждой из одной или нескольких плоскостей дренирования (115) для идентификации, таким образом, блоков сетки (103), содержащихся в объеме дренирования (111),
определяют оценку эффективного объема дренирования (111) скважины (119) на каждой из совокупности ньютоновых итераций в полностью неявном имитаторе коллектора.
41. Компьютерно-реализуемый способ определения совокупности приближенных статических давлений в скважине для одной или нескольких скважин произвольной формы (119) в коллекторе (101), компьютерно-реализуемый способ отличается тем, что содержит этапы, на которых:
определяют положение каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из перфорированного блока сетки (113) на плоскости дренирования (115) сегмента скважины (117) скважины, для оценивания, таким образом, эффективного объема дренирования (111) скважины (119), причем
- для имитаций коллектора, где множественные скважины (119) состязаются за отдачу флюида, определяют положение соответствующей границы дренирования (110) как
- положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) из перфорированного блока сетки (113) на плоскости дренирования (115), и
- для имитаций коллектора, где множественные скважины (119) незначительно состязаются за отдачу флюида, в результате чего, возможность впервые встретить изменение потока флюида в силу конкуренции скважин не предполагается, определяют положение соответствующей границы дренирования (110) как соответствующее одно из следующих положений, находящихся на меньшем расстоянии от перфорированного блока сетки (113), чем другие из следующих положений:
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) из перфорированного блока сетки (113) на плоскости дренирования (115), и
- положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где расстояние соответствующего одного из совокупности блоков сетки (103) из перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) из перфорированного блока сетки (113) на плоскости дренирования (115),
определяют динамическое давление блока сетки каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), и
определяют усредненное по объему перового пространства динамическое давление блока сетки, по меньшей мере, существенного подмножества совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119) для задания, таким образом, оценочного статического давления в скважине для скважины (119).
42. Компьютерно-реализуемый способ по п.41, в котором плоскость дренирования (115) является выбранной плоскостью дренирования (115) соответствующего сегмента скважины (117) скважины (119), которая ориентирована перпендикулярно направлению сегмента скважины (117) скважины (119) и включает в себя грань (105) перфорированного блока сетки (113) и грань (105) каждого из совокупности блоков сетки (103), соседствующего с ним, причем способ дополнительно отличается тем, что содержит этапы, на которых:
осуществляют для каждого из совокупности блоков сетки (103), содержащихся в эффективном объеме дренирования (111) скважины (119), и на каждой ньютоновой итерации вычисления каждого временного шага имитатора коллектора на основе сетки конечных разностей, этап вычисления совокупности векторов потока флюида (107) на каждой из соответствующей совокупности границ раздела с совокупностью соседних блоков сетки (103), и
осуществляют на каждом из совокупности блоков сетки (103), расположенных в каждом отдельном одном из совокупности направлений, задающих совокупность траекторий (120), выходящих наружу из перфорированного блока сетки (113) на плоскости дренирования (115), для определения, таким образом, каждой границы дренирования (110) эффективного объема дренирования (111) скважины (119) этапы, на которых
- отслеживают знак каждого из совокупности расчетных векторов потока флюида (107) для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113) на плоскости дренирования (115),
- сравнивают величину фазового потока для каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113), с величиной соответствующего фазового потока на перфорированном блоке сетки (113) для плоскости дренирования (115), и
- определяют расстояние каждого из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120), выходящей из перфорированного блока сетки (113) на плоскости дренирования (115).
43. Компьютерно-реализуемый способ по п.41, в котором на этапе определения положения каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из перфорированного блока сетки (113) на плоскости дренирования (115) определяют положение каждой отдельной границы дренирования (110) вдоль каждой соответствующей отдельной одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113), как для имитаций коллектора, где множественные скважины (119) состязаются за отдачу флюида, так и для имитаций коллектора, где множественные скважины (119) не состязаются за отдачу флюида, как одно из следующих положений, находящихся на меньшем расстоянии от перфорированного блока сетки (113), чем все остальные из следующих положений:
положение одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где смена знака соответствующего вектора потока флюида (107) для одного из совокупности блоков сетки (103) относительно знака соседнего ранее пройденного одного из совокупности блоков сетки (103) вдоль соответствующей траектории (120) впервые встречается при прохождении соответствующей траектории (120) из перфорированного блока сетки (113),
положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где отношение фазовых потоков, а именно величины фазового потока соответствующего одного из совокупности блоков сетки (103) к величине соответствующего фазового потока на соответствующем перфорированном блоке сетки (113), впервые принимает значение ниже заранее заданного предельного отношения при прохождении соответствующей траектории (120) из перфорированного блока сетки (113), и
положение одного из совокупности блоков сетки (103), расположенных вдоль соответствующей траектории (120) из перфорированного блока сетки (113), где расстояние соответствующего одного из совокупности блоков сетки (103) из перфорированного блока сетки (113) впервые принимает значение, которое больше или равно заранее заданному предельному расстоянию, при прохождении соответствующей траектории (120) из перфорированного блока сетки (113).
44. Компьютерно-реализуемый способ по п.41, в котором на этапе определения каждая граница дренирования (110) эффективного объема дренирования (111) скважины (119)
аппроксимируют одну или несколько из определенных границ дренирования (110) с использованием непрерывных прямых линий.
45. Компьютерно-реализуемый способ по п.41 дополнительно отличается тем, что на этапе
определения оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) плоскости дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111)
упрощают форму объема дренирования (111) скважины (119) на каждом временном шаге имитатора коллектора до формы прямоугольной трехмерной ячейки объема дренирования (111), содержащей совокупность слоев (118), причем, по меньшей мере, один из слоев (118) имеет другую толщину, чем, по меньшей мере, один другой из слоев (118).
46. Компьютерно-реализуемый способ по любому из пп.41-45, в котором модель скважины (119) содержит совокупность плоскостей дренирования (115), и эффективный объем дренирования (111) скважины (119) оценивается на основании определенных границ дренирования (110) двух перпендикулярных плоскостей дренирования (115) из совокупности плоскостей дренирования (115).
47. Компьютерно-реализуемый способ по любому из пп.41-45, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления этапа определения положения каждой отдельной границы дренирования (110) плоскости дренирования (115) сегмента скважины (117) для скважины (119),
причем плоскость дренирования (115) является мультидоменной плоскостью дренирования (115), имеющей границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего скважину (119), и
когда внешний предел границы домена (203) первого процессора (32) достигается при прохождении мультидоменной плоскости дренирования (115) вдоль одной из совокупности траекторий (120), выходящих из соответствующего перфорированного блока сетки (113), осуществляют этапы, на которых
- сохраняют данные расстояния для пройденного блока сетки (103), расположенного вдоль одной из совокупности траекторий (120) на внешнем пределе границы домена (203) первого процессора (32),
- сохраняют данные направления, указывающие направление прохождения вдоль одной из совокупности траекторий (120), на которой достигается соответствующий внешний предел границы домена (203) первого процессора (32), и
- синхронизируют данные расстояния и направления со вторым из совокупности процессоров (32), имеющим границу домена (203) рядом с внешним пределом границы домена (203) первого процессора (32).
48. Компьютерно-реализуемый способ по любому из пп.41-45, в котором модель коллектора делится среди совокупности процессоров (32) для осуществления этапа определения положения каждой отдельной границы дренирования (110) плоскости дренирования (115) для скважины (119),
причем плоскость дренирования (115) является мультидоменной плоскостью дренирования (115), имеющей границу дренирования (110), которая выходит за пределы границы домена (203) первого процессора (32) из совокупности процессоров (32), содержащего соответствующую скважину (119) в, по меньшей мере, одном из совокупности направлений, и
дополнительно содержащий этапы, на которых формируют для каждого из совокупности направлений, в которых соответствующая граница дренирования (110) выходит за пределы границы домена (203) первого процессора (32), контур с длиной контура, равной количеству совокупности процессоров (32).
49. Компьютерно-реализуемый способ по любому из пп.41-45, дополнительно отличается тем, что содержит этап, на котором:
определяют оценку эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) плоскости дренирования (115) для идентификации, таким образом, блоков сетки (103), содержащихся в объеме дренирования (111),
причем этап определения эффективного объема дренирования (111) осуществляется для неоднородного коллектора (101) с применением численного решения, и
оценочное статическое давление в скважине для скважины (119) равно оценочному среднему статическому давлению в скважине (р_static) эффективного объема дренирования (111) определенному согласно следующему вычислению:
p_static=[Sum(Vp(k))·P(k)] Sum(Vp(k)),
где k относится ко всем блокам сетки (103) в эффективном объеме дренирования (111),
Vp(k) - объем перового пространства блока сетки (103) k в эффективном объеме дренирования (111),
Р(k) - динамическое давление блока сетки (103) k, и
Sum - оператор суммирования.
50. Компьютерно-реализуемый способ по любому из пп.41-45, дополнительно отличается тем, что на этапе
определения оценки эффективного объема дренирования (111) скважины (119) в соответствии с определенными границами дренирования (110) плоскости дренирования (115) для идентификации, таким образом, блоков сетки (103, 113), содержащихся в эффективном объеме дренирования (111)
определяют оценку эффективного объема дренирования (111) скважины (119) на каждой из совокупности ньютоновых итераций в полностью неявном имитаторе коллектора.
RU2011149638/08A 2009-05-07 2010-05-06 Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора RU2549147C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US17626109P 2009-05-07 2009-05-07
US61/176,261 2009-05-07
PCT/US2010/033848 WO2010129754A2 (en) 2009-05-07 2010-05-06 Systems, computer implemented methods, and computer readable program products to compute approximate well drainage pressure for a reservoir simulator

Publications (2)

Publication Number Publication Date
RU2011149638A true RU2011149638A (ru) 2013-06-20
RU2549147C2 RU2549147C2 (ru) 2015-04-20

Family

ID=42542983

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2011149638/08A RU2549147C2 (ru) 2009-05-07 2010-05-06 Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора

Country Status (5)

Country Link
US (1) US8589135B2 (ru)
CN (1) CN102640155B (ru)
NO (1) NO337774B1 (ru)
RU (1) RU2549147C2 (ru)
WO (1) WO2010129754A2 (ru)

Families Citing this family (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9488752B2 (en) * 2013-06-04 2016-11-08 Saudi Arabian Oil Company System for computing the radius of investigation in a radial, composite reservoir system
CA2702965C (en) 2007-12-13 2014-04-01 Exxonmobil Upstream Research Company Parallel adaptive data partitioning on a reservoir simulation using an unstructured grid
CN102640163B (zh) 2009-11-30 2016-01-20 埃克森美孚上游研究公司 用于储层模拟的适应性牛顿法
US9134454B2 (en) 2010-04-30 2015-09-15 Exxonmobil Upstream Research Company Method and system for finite volume simulation of flow
WO2012003007A1 (en) 2010-06-29 2012-01-05 Exxonmobil Upstream Research Company Method and system for parallel simulation models
US9058445B2 (en) 2010-07-29 2015-06-16 Exxonmobil Upstream Research Company Method and system for reservoir modeling
CA2805446C (en) 2010-07-29 2016-08-16 Exxonmobil Upstream Research Company Methods and systems for machine-learning based simulation of flow
CA2803066A1 (en) 2010-07-29 2012-02-02 Exxonmobil Upstream Research Company Methods and systems for machine-learning based simulation of flow
WO2012039811A1 (en) 2010-09-20 2012-03-29 Exxonmobil Upstream Research Company Flexible and adaptive formulations for complex reservoir simulations
EA201390853A1 (ru) * 2010-12-16 2014-01-30 Лэндмарк Графикс Корпорейшн Система и способ декомпозиции двухмерной области при параллельном моделировании природного резервуара углеводородов
AU2011360213B2 (en) * 2011-02-21 2016-09-29 Exxonmobil Upstream Research Company Method and system for field planning
CN103827441A (zh) * 2011-07-28 2014-05-28 普拉德研究及开发股份有限公司 用于执行井眼压裂作业的系统和方法
US8731892B2 (en) 2011-08-02 2014-05-20 Saudi Arabian Oil Company Systems and program product for performing a fully automated workflow for well performance model creation and calibration
US8688426B2 (en) 2011-08-02 2014-04-01 Saudi Arabian Oil Company Methods for performing a fully automated workflow for well performance model creation and calibration
CA2843929C (en) 2011-09-15 2018-03-27 Exxonmobil Upstream Research Company Optimized matrix and vector operations in instruction limited algorithms that perform eos calculations
FR2987149B1 (fr) * 2012-02-16 2014-10-31 IFP Energies Nouvelles Procede d'exploitation d'un gisement a partir d'une technique de selection des positions de puits a forer
US10036829B2 (en) 2012-09-28 2018-07-31 Exxonmobil Upstream Research Company Fault removal in geological models
KR101589798B1 (ko) * 2013-12-30 2016-01-28 연세대학교 산학협력단 해외 가스전 개발 프로젝트의 지속가능성 평가 시스템 및 방법
WO2016018723A1 (en) 2014-07-30 2016-02-04 Exxonmobil Upstream Research Company Method for volumetric grid generation in a domain with heterogeneous material properties
AU2015339883B2 (en) 2014-10-31 2018-03-29 Exxonmobil Upstream Research Company Methods to handle discontinuity in constructing design space for faulted subsurface model using moving least squares
CA2963416A1 (en) 2014-10-31 2016-05-06 Exxonmobil Upstream Research Company Handling domain discontinuity in a subsurface grid model with the help of grid optimization techniques
US10242136B2 (en) * 2015-05-20 2019-03-26 Saudi Arabian Oil Company Parallel solution for fully-coupled fully-implicit wellbore modeling in reservoir simulation
EP3423672A1 (en) * 2016-03-04 2019-01-09 Saudi Arabian Oil Company Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation
US10233749B2 (en) 2017-05-03 2019-03-19 Saudi Arabian Oil Company Multi-layer reservoir well drainage region
US10570706B2 (en) 2017-06-23 2020-02-25 Saudi Arabian Oil Company Parallel-processing of invasion percolation for large-scale, high-resolution simulation of secondary hydrocarbon migration
WO2019048599A1 (en) * 2017-09-08 2019-03-14 Roxar Software Solutions As WELL FRACTURE MODELING
US10519768B2 (en) 2018-02-21 2019-12-31 Saudi Arabian Oil Company Systems and methods for operating hydrocarbon wells to inhibit breakthrough based on reservoir saturation
CN110021220B (zh) * 2019-04-28 2024-03-26 中国华能集团清洁能源技术研究院有限公司 一种积木式地热尾水回灌分析系统及使用方法
US11009623B2 (en) 2019-07-16 2021-05-18 Saudi Arabian Oil Company Calculating shut-in bottom-hole pressure in numerical reservoir simulations
US11560785B2 (en) * 2020-01-28 2023-01-24 Enverus, Inc. Determining spacing between wellbores
US11193370B1 (en) 2020-06-05 2021-12-07 Saudi Arabian Oil Company Systems and methods for transient testing of hydrocarbon wells
CN111594161B (zh) * 2020-07-08 2022-03-11 西南石油大学 一种巨厚或多层油气藏的编织井网提高采收率方法
CN114371116B (zh) * 2021-12-30 2024-06-18 北京红山信息科技研究院有限公司 一种基于渗透率的路测质量评估方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU364948A1 (ru) * 1971-04-20 1972-12-28 В. П. Старинский , В. В. Сержантов Научно исследовательский институт водных проблем УСТРОЙСТВО дл МОДЕЛИРОВАНИЯ ВОДОЗАБОРОВ ПОДЗЕМНЫХ ВОД
AU2001293809A1 (en) * 2000-09-12 2002-03-26 Sofitech N.V. Evaluation of multilayer reservoirs
CA2538577C (en) * 2003-09-30 2010-08-03 Exxonmobil Upstream Research Company Characterizing connectivity in reservoir models using paths of least resistance
CN2884196Y (zh) * 2005-11-21 2007-03-28 西南石油学院 高温高压钻井液密度测量装置
FR2925726B1 (fr) * 2007-12-20 2010-04-23 Inst Francais Du Petrole Methode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'echange geologique et transitoire entre blocs matriciels et fractures
US8190414B2 (en) * 2008-03-26 2012-05-29 Exxonmobil Upstream Research Company Modeling of hydrocarbon reservoirs containing subsurface features

Also Published As

Publication number Publication date
US20100286971A1 (en) 2010-11-11
NO20111482A1 (no) 2012-02-07
RU2549147C2 (ru) 2015-04-20
WO2010129754A2 (en) 2010-11-11
CN102640155B (zh) 2015-01-28
CN102640155A (zh) 2012-08-15
US8589135B2 (en) 2013-11-19
NO337774B1 (no) 2016-06-20
WO2010129754A3 (en) 2012-01-05

Similar Documents

Publication Publication Date Title
RU2011149638A (ru) Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора
Lacasta et al. An optimized GPU implementation of a 2D free surface simulation model on unstructured meshes
Tsubota et al. Simulation study on effects of hematocrit on blood flow properties using particle method
Brodtkorb et al. Efficient shallow water simulations on GPUs: Implementation, visualization, verification, and validation
Adelsberger et al. 3D incompressible two-phase flow benchmark computations for rising droplets
RU2510969C2 (ru) Компьютеризованный способ для оптимизированного по затратам расчета аэродинамических сил на летательном аппарате
CN101661514A (zh) 一种油藏黑油模型数值模拟系统
NO339744B1 (no) Hybrid, lokal ikke-avstemmende metode for flerfasestrømningssimuleringer i heterogent sprukne medier
AU2011289788A1 (en) Reservoir upscaling method with preserved transmissibility
CN102968529A (zh) 一种供水管网模型计算结果不确定性区间的量化方法
CN103838852A (zh) 一种快速查找多块结构化网格对接关系的方法
CN106250102A (zh) 交错网格有限差分正演模拟优化的方法
CN104820826A (zh) 一种基于数字高程模型的坡形提取与识别方法
US8346584B2 (en) Method and apparatus for determining design modification effects on a computerized engineering model
Li et al. Modeling subgrid-scale topographic effects on shallow marsh hydrodynamics and salinity transport
Gao Pipe roughness estimation in water distribution networks using head loss adjustment
He A new adaptive Cartesian-grid CIP method for computation of violent free-surface flows
CN114065593A (zh) 一种多尺度裂隙基质系统中地下水耦合模拟的方法和装置
CN102063577A (zh) 一种基于等势面的地下水渗流量计算方法
KR101106548B1 (ko) 텐서형 와점성계수를 가진 2차원 하천흐름모형을 이용하여 천수흐름을 해석하는 방법
CN107832482A (zh) 致密储层多尺度裂缝网络建模及模拟方法
Xiong et al. A novel stencil selection method for the gradient reconstruction on unstructured grid based on OpenFOAM
Owen et al. Recent ship hydrodynamics developments in the parallel two-fluid flow solver Alya
CN111310349A (zh) 适用于离散元计算信息连续化展示的数据处理分析方法
CN102930087B (zh) 一种模拟仿真技术中的相邻粒子搜索方法