[go: up one dir, main page]

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

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

Info

Publication number
RU2014104788A
RU2014104788A RU2014104788/03A RU2014104788A RU2014104788A RU 2014104788 A RU2014104788 A RU 2014104788A RU 2014104788/03 A RU2014104788/03 A RU 2014104788/03A RU 2014104788 A RU2014104788 A RU 2014104788A RU 2014104788 A RU2014104788 A RU 2014104788A
Authority
RU
Russia
Prior art keywords
wetting
fluids
sample
wetting fluids
fluid
Prior art date
Application number
RU2014104788/03A
Other languages
English (en)
Other versions
RU2593853C2 (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 RU2014104788A publication Critical patent/RU2014104788A/ru
Application granted granted Critical
Publication of RU2593853C2 publication Critical patent/RU2593853C2/ru

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials
    • G01N33/241—Earth materials for hydrocarbon content
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Preparation Of Compounds By Using Micro-Organisms (AREA)
  • Sampling And Sample Adjustment (AREA)

Abstract

1. Способ моделирования движения отдельных фаз смачивающих текучих сред и несмачивающих текучих сред через пористую среду, содержащий следующие этапы:а) создание трехмерного цифрового представления пористой среды (Образец), содержащего весь объем текучих сред, представляющих собой смачивающие текучие среды и несмачивающие текучие среды,б) определение первой доли всего объема текучих сред, которая содержит смачивающие текучие среды, и второй доли всего объема текучих сред, которая содержит несмачивающие текучие среды,в) определение значения скорости потока всего объема текучих сред, проходящих через Образец,г) оценка свойств смачивающих текучих сред и несмачивающих текучих сред,д) определение начальных условий для насыщения (Sw) смачивающих текучих сред, насыщения (Sn) несмачивающих текучих сред, впускного давления (Pw) смачивающих текучих сред и впускного давления (Pn) несмачивающих текучих сред,е) установление условий на впускной поверхности Образца, где несмачивающие текучие среды и смачивающие текучие среды входят в поры образца в отдельных и отличающихся площадях, иж) вычисление давлений, насыщений и векторов скорости внутри Образца,з) вычисление скоростей (Qn) потока несмачивающих текучих сред, проходящих через Образец, скоростей (Qw) потока смачивающих текучих сред, проходящих через Образец, и давлений на выпуске Образца,и) повторение этапов а) - з) для заданного числа приращений времени t, и периодическое регулирование впускных давлений Pn и Pw с использованием алгоритма регулирования с обратной связью, при этом получают значения Qn и Qw квазистационарного состояния.2. Способ по п. 1, в котором пористая среда я�

Claims (48)

1. Способ моделирования движения отдельных фаз смачивающих текучих сред и несмачивающих текучих сред через пористую среду, содержащий следующие этапы:
а) создание трехмерного цифрового представления пористой среды (Образец), содержащего весь объем текучих сред, представляющих собой смачивающие текучие среды и несмачивающие текучие среды,
б) определение первой доли всего объема текучих сред, которая содержит смачивающие текучие среды, и второй доли всего объема текучих сред, которая содержит несмачивающие текучие среды,
в) определение значения скорости потока всего объема текучих сред, проходящих через Образец,
г) оценка свойств смачивающих текучих сред и несмачивающих текучих сред,
д) определение начальных условий для насыщения (Sw) смачивающих текучих сред, насыщения (Sn) несмачивающих текучих сред, впускного давления (Pw) смачивающих текучих сред и впускного давления (Pn) несмачивающих текучих сред,
е) установление условий на впускной поверхности Образца, где несмачивающие текучие среды и смачивающие текучие среды входят в поры образца в отдельных и отличающихся площадях, и
ж) вычисление давлений, насыщений и векторов скорости внутри Образца,
з) вычисление скоростей (Qn) потока несмачивающих текучих сред, проходящих через Образец, скоростей (Qw) потока смачивающих текучих сред, проходящих через Образец, и давлений на выпуске Образца,
и) повторение этапов а) - з) для заданного числа приращений времени t, и периодическое регулирование впускных давлений Pn и Pw с использованием алгоритма регулирования с обратной связью, при этом получают значения Qn и Qw квазистационарного состояния.
2. Способ по п. 1, в котором пористая среда является горной породой, грунтом, цеолитом, биологической тканью, деревом, пробкой, цементом, керамикой, песком, глиной, неорганическим составом, органическим составом или металлом.
3. Способ по п. 1, в котором Образец содержит несколько упорядоченных плоскостей вокселов, при этом каждый из вокселов представляет пору (воксел поры) или твердое вещество (воксел зерна).
4. Способ по п. 1, в котором свойства смачивающих текучих сред содержат вязкость, краевой угол, межфазное натяжение, другие физические или химические свойства или любые их комбинации.
5. Способ по п. 1, в котором свойства несмачивающих текучих сред содержат вязкость, краевой угол, межфазное натяжение, другие физические или химические свойства или любые их комбинации.
6. Способ по п. 1, в котором Образец содержит (a) впускную поверхность и выпускную поверхность, при этом впускная поверхность и выпускная поверхность параллельны друг другу, и (б) три или больше поверхностей ортогональных впускной поверхности и выпускной поверхности, при этом три или больше ортогональных поверхностей являются непроницаемыми для потока смачивающих текучих сред и несмачивающих текучих сред.
7. Способ по п. 6, в котором впускная поверхность дополнительно содержит буферную зону, параллельную впускной поверхности, содержащую по меньшей мере одну плоскость вокселов.
8. Способ по п. 7, в котором межфазное натяжение между смачивающей текучей средой и несмачивающей текучей средой устанавливается на ноль для всех вычислений в буферной зоне.
9. Способ по п. 8, в котором межфазное натяжение между смачивающей текучей средой и несмачивающей текучей средой устанавливается на ноль, и вязкости смачивающей текучей среды и несмачивающей текучей среды увеличиваются по меньшей мере около в 10 раз для всех вычислений в буферной зоне.
10. Способ по п. 1, в котором несмачивающие текучие среды и смачивающие текучие среды входят в Образец через поры на впускной поверхности образца.
11. Способ по п. 10, в котором поры на впускной поверхности содержат отдельные и отличающиеся площади, образованные с помощью размещения вокселов пор непосредственно примыкающими к вокселам зерен для потока смачивающих текучих сред, (Aw), и остальные вокселы пор размещаются для потока несмачивающих текучих сред, (An).
12. Способ по п. 11, в котором Aw увеличивается с помощью дополнительного размещения вокселов пор, примыкающих к Aw для потока смачивающих текучих сред, (Aw) и остальные вокселы пор размещаются для потока несмачивающих текучих сред, (An).
13. Способ по п. 11, в котором (сумма вокселов в An)/((сумма вокселов в An) + (сумма вокселов в Aw)) составляет приблизительно 0,5 или меньше.
14. Способ по п. 1, в котором вычисление содержит расчетную гидродинамику.
15. Способ по п. 14, в котором расчетная гидродинамика содержит решеточный метод Больцмана.
16. Способ по п. 1, в котором приращение времени t выражается в секундах, миллисекундах или других единицах времени.
17. Способ по п. 1, в котором число временных приращений составляет 10000 или больше.
18. Способ по п. 1, в котором алгоритм регулирования с обратной связью содержит отдельный алгоритм регулирования с обратной связью для установления впускного давления для смачивающей текучей среды и отдельный алгоритм регулирования с обратной связью для установления впускного давления для несмачивающей текучей среды, при этом впускное давление для смачивающей текучей среды и впускное давление для несмачивающей текучей среды устанавливаются независимо.
19. Способ по п. 1, в котором алгоритм регулирования с обратной связью содержит один алгоритм регулирования с обратной связью для установления впускного давления как для смачивающей, так и для несмачивающей текучей среды, при этом впускное давление для смачивающей текучей среды и впускное давление для несмачивающей текучей среды являются равными.
20. Способ по п. 1, в котором алгоритм регулирования с обратной связью содержит контур пропорционально-интегрально-дифференциального управления, адаптивное управление, иерархическое управление, интеллектуальное управление, оптимальное управление, робастное управление, нейронно-сетевое управление, управление с нечеткой логикой или стохастическое управление.
21. Способ по п. 1, в котором алгоритм регулирования с обратной связью является алгоритмом регулирования с отрицательной обратной связью.
22. Способ по п. 1, в котором алгоритм регулирования с обратной связью содержит контур пропорционально-интегрально-дифференциального (ПИД) управления.
23. Способ по п. 22, в котором контур ПИД управления содержит входное отклонение Ew и выдает новое впускное давление Pw, где
Pw=Pi+Pi* π
Figure 00000001
 w
Pi - начальное давление, установленное в начале моделирования
π
Figure 00000002
 w - f(Ew) так что K p ⋅ E w + K I ∫ E w d t + K D d E w d t
Figure 00000003
Kp - константа пропорционального регулирования,
KI - константа интегрального регулирования,
KD - константа дифференциального регулирования,
Ew=QwT-QW, и
QwT=QT∙Fw,
QT - заданная суммарная скорость потока, проходящего через Образец,
Fw - доля смачивающей текучей среды, которая входит во впускную поверхность образца,
QW - среднее значение Qwt, и
Qwt - вычисленная скорость потока смачивающей текучей среды и интервал времени t.
24. Способ по п. 23, в котором QW содержит арифметическое средневзвешенное значение, геометрическое средневзвешенное значение, гармоническое средневзвешенное значение, простое скользящее среднее, экспоненциально взвешенное скользящее среднее или полученное другим способом среднего ряда чисел.
25. Способ по п. 22, в котором контур ПИД управления содержит входное отклонение En и выдает новое впускное давление Pn, где
Pn=Pi+Pi* π
Figure 00000002
 n
Pi - начальное давление, установленное в начале моделирования
π
Figure 00000002
 n=f(En) так что K p ⋅ E n + K I ∫ E n d t + K D d E n d t
Figure 00000004
Kp - константа пропорционального регулирования,
KI - константа интегрального регулирования,
KD - константа дифференциального регулирования,
En=QnT-QN, и
QnT=QT∙Fn,
QT - заданная суммарная скорость потока, проходящего через Образец,
Fn - доля несмачивающей текучей среды, которая входит во впускную поверхность Образца,
QN - среднее значение Qnt, и
Qnt - вычисленная скорость потока несмачивающей текучей среды и интервал времени t.
26. Способ по п. 25, в котором QN содержит арифметическое средневзвешенное значение, геометрическое средневзвешенное значение, гармоническое средневзвешенное значение, простое скользящее среднее, экспоненциально взвешенное скользящее среднее или полученное другим способом среднего ряда чисел.
27. Способ по п. 22, в котором контур ПИД управления содержит входное отклонение Ec и выдает новое впускное давление Pc, при этом
Pc=Pi+Pi* π
Figure 00000002
 c
Pi = начальное давление, установленное в начале моделирования
π
Figure 00000002
 с = f(Ec) так что K p ⋅ E c + K I ∫ E c d t + K D d E c d t
Figure 00000005
Kp = константа пропорционального регулирования,
KI = константа интегрального регулирования,
KD = константа дифференциального регулирования,
Ec=QT-QC, и
QT = заданная суммарная скорость потока, проходящего через Образец,
QC = среднее значение Qwt+Qnt,
Qnt = вычисленная скорость потока несмачивающей текучей среды и интервал времени t, и
Qwt = вычисленная скорость потока смачивающей текучей среды и интервал времени t.
28. Способ по п. 27, в котором QT содержит арифметическое средневзвешенное значение, геометрическое средневзвешенное значение, гармоническое средневзвешенное значение, простое скользящее среднее, экспоненциально взвешенное скользящее среднее или полученное другим способом среднее ряда чисел.
29. Способ по п. 1, в котором периодическое регулирование впускных давлений происходит через каждые около 10 или больше приращений времени.
30. Способ по п. 29, в котором число приращений времени для последовательных периодических регулировок впускного давления является различным.
31. Способ по п. 1, в котором периодическое регулирование впускных давлений проводится чаще в первой половине общего времени моделирования, чем во второй половине моделирования.
32. Способ по п. 31, в котором периодическое регулирование впускных давлений в первой половине общего времени моделирования проводится по меньшей мере в 10 раз чаще, чем во второй половине.
33. Способ по п. 1, в котором квазистационарное состояние является состоянием, где вычисленные значения Qn, Qw, Pn, Pw и/или насыщения изменяются не больше заданного значения.
34. Способ по п. 1, в котором предопределенное число приращений времени t, устанавливается достаточно большим для получения квазистационарного состояния.
35. Способ определения относительной проницаемости смачивающих и несмачивающих текучих сред, проходящих через пористую среду, содержащий
a) определение ряда пар несмачивающих текучих сред и смачивающих текучих сред, причем каждая пара подлежит продавливанию через Образец,
б) установление начального насыщения Образца,
в) продавливание каждой пары несмачивающих и смачивающих текучих сред через Образец,
г) запись вычисленных значений QN, QW и насыщения Sw смачивающей текучей среды для каждой пары смачивающих и несмачивающих текучих сред,
д) вычисление значений относительной проницаемости kw смачивающей текучей среды; вычисление значений относительной проницаемости kn несмачивающей текучей среды; и вычисление значений водонасыщения Sw, и
е) построение графика значений kw и kn в зависимости от Sw.
36. Способ по п. 35, в котором пористая среда является горной породой, грунтом, цеолитом, биологической тканью, деревом, пробкой, цементом, керамикой, песком, глиной, горной породой, неорганическим составом, органическим составом или металлом.
37. Способ по п. 35, в котором пара несмачивающей и смачивающей текучей среды содержит состав с долями (Fn, Fw) несмачивающей текучей среды и дробной величиной смачивающей текучей среды.
38. Способ по п. 37, в котором пара несмачивающей и смачивающей текучей среды содержит множество пар долей (Fn, Fw), при этом любая комбинация Fn и Fw в сумме равна 1.
39. Способ по п. 37, в котором пара несмачивающей и смачивающей текучей среды содержит множество пар долей (Fn′, Fw′), при этом Fn′ вычисляется, как величина (Fn×R), и Fw′ вычисляется, как величина (1-(Fn×R)), при этом R является соотношением вязкостей фазы низкой вязкости и фазы высокой вязкости смачивающихся и несмачивающихся текучих сред, и любая комбинация долей Fn′ и Fw′ в сумме равна 1.
40. Способ по п. 37, в котором пары несмачивающей и смачивающей текучей среды содержат упорядоченный ряд значений, в котором Fn уменьшается ступенчато до нуля и затем увеличивается до 1,0.
41. Способ по п. 35, в котором начальное насыщение Образца является полным насыщением смачивающей текучей средой, Sw=1,0, и полным насыщением несмачивающей текучей средой, Sn=1,0.
42. Способ по п. 35, в котором начальное насыщение Образца является конечным насыщением условий предыдущего моделирования.
43. Способ по п. 35, в котором доля площади An пор на впускной поверхности, отведенная для нагнетания несмачивающей текучей среды, уменьшается, когда Fn меньше около 0,2.
44. Способ по п. 43, в котором An уменьшена до около 0,4 или меньше.
45. Система определения скорости отдельных фаз многофазного многокомпонентного потока, проходящего через пористую среду, содержащая:
a) сканер, выполненный с возможностью получения трехмерного цифрового изображения пористой среды,
б) компьютер, содержащий по меньшей мере один процессор с функциональной возможностью исполнения компьютерной программы для группирования элементов в трехмерное цифровое изображение твердого вещества (зерен) и пор (полостей),
в) компьютер (аналогичный компьютеру б) или отличающийся)), содержащий по меньшей мере один процессор с функциональной возможностью исполнения компьютерной программы для выполнения вычислений, при этом вычисления содержат (I) создание трехмерного цифрового представления пористой среды (Образец), содержащего весь объем текучих сред, представляющих собой смачивающие текучие среды и несмачивающие текучие среды, (II) определение первой доли всего объема текучих сред, которая содержит смачивающие текучие среды, и определение второй доли всего объема текучих сред, которая содержит несмачивающие текучие среды, (III) определение значения скорости потока всего объема текучих сред, проходящих через Образец, (IV) оценка свойств смачивающих текучих сред и несмачивающих текучих сред, (V) определение начальных условий для насыщения (Sw) смачивающих текучих сред, насыщения (Sn) несмачивающих текучих сред, впускного давления (Pw) смачивающих текучих сред и впускного давления (Pn) несмачивающих текучих сред, (VI) установление условий на впускной поверхности Образца, где несмачивающие текучие среды и смачивающие текучие среды входят в поры образца в отдельных и отличающихся площадях, и (VII) вычисление давлений, насыщений и векторов скорости внутри Образца пористой среды, (VIII) вычисление скоростей (Qn) потока несмачивающих текучих сред, проходящих через Образец, скоростей (Qw) потока смачивающих текучих сред, проходящих через Образец и давлений на выпуске Образца, (IX) повторение этапов (I)-(VIII) для заданного числа приращений времени t, и (Х) периодическое регулирование впускных давлений Pn и Pw с использованием алгоритма регулирования с обратной связью, где достигаются значения Qn и Qw квазистационарного состояния, и
г) по меньшей мере одно устройство для отображения, печати или сохранения результатов вычислений.
46. Система по п. 45, в которой сканер представляет собой компьютерный томографический (КТ) сканер, сканирующий электронный микроскоп, сканирующий электронный микроскоп со сфокусированным ионным лучом или аналогичное устройство, выполненное с возможностью получения трехмерного цифрового изображения пористой среды.
47. Система по п. 45, в которой устройство содержит запоминающее устройство для сохранения с возможностью осуществления выборки результатов вычислений.
48. Компьютерный программный продукт на машиночитаемом носителе который, когда выполняется на контроллере в компьютеризованном устройстве побуждает компьютер выполнять способ по п. 45.
RU2014104788/03A 2011-07-12 2012-07-02 Способ моделирования движения отдельных фаз многофазного/многокомпонентного потока, проходящего через пористую среду RU2593853C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161506680P 2011-07-12 2011-07-12
US61/506,680 2011-07-12
PCT/US2012/045220 WO2013009512A2 (en) 2011-07-12 2012-07-02 Method for simulating fractional multi-phase/multi-component flow through porous media

Publications (2)

Publication Number Publication Date
RU2014104788A true RU2014104788A (ru) 2015-08-20
RU2593853C2 RU2593853C2 (ru) 2016-08-10

Family

ID=46598929

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2014104788/03A RU2593853C2 (ru) 2011-07-12 2012-07-02 Способ моделирования движения отдельных фаз многофазного/многокомпонентного потока, проходящего через пористую среду

Country Status (10)

Country Link
US (1) US9183326B2 (ru)
EP (1) EP2732135B1 (ru)
CN (1) CN103906893A (ru)
AU (1) AU2012283030B2 (ru)
BR (1) BR112014000758B1 (ru)
CA (1) CA2840942C (ru)
CO (1) CO6930314A2 (ru)
MX (1) MX350511B (ru)
RU (1) RU2593853C2 (ru)
WO (1) WO2013009512A2 (ru)

Families Citing this family (67)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9140117B2 (en) * 2012-07-13 2015-09-22 Ingrain, Inc. Method for evaluating relative permeability for fractional multi-phase, multi-component fluid flow through porous media
US9047513B2 (en) 2012-08-10 2015-06-02 Ingrain, Inc. Method for improving the accuracy of rock property values derived from digital images
CN103176055A (zh) * 2013-02-07 2013-06-26 吉林大学 层状结构岩石垂直纹理方向介电性质测量样品的制备方法
CA2912674C (en) * 2013-05-16 2021-05-18 Exa Corporation Mass exchange model for relative permeability simulation
CN103344540B (zh) * 2013-06-28 2016-02-03 武汉理工大学 一种连续表征水泥石密实度的测试方法
CN103806905B (zh) * 2013-07-31 2016-08-24 中国石油大学(华东) 制备具有双重孔隙的碳酸盐岩储层物理模型的装置及方法
AU2014307046B2 (en) * 2013-08-13 2018-05-24 Schlumberger Technology B.V. Digital core sensitivity analysis
CN105683494A (zh) * 2013-09-16 2016-06-15 兰德马克绘图国际公司 伪相生产模拟：在储层流动模拟中经由连续相似的阶跃函数相对渗透率受控模型来评估准多相流生产以便对多个岩石物理学实现排序的信号处理方法
WO2015102508A1 (en) * 2013-12-30 2015-07-09 Schlumberger Holdings Limited Method for estimating petrophysical properties of a hydrocarbon reservoir
US10705250B2 (en) * 2014-01-27 2020-07-07 Schlumberger Technology Corporation Generation of isotherm datasets for reservoir volumetric estimation
CA2942729A1 (en) 2014-03-18 2015-09-24 Exa Corporation Fluid blob tracking for evaluation of multiphase flow simulations
CN104535472B (zh) * 2014-12-22 2017-04-12 中国石油大学(北京) 一种煤岩动态渗透率检测方法及装置
GB2554013B (en) * 2015-05-01 2021-03-24 Geoquest Systems Bv Multiphase flow in porous media
CN105043935A (zh) * 2015-05-18 2015-11-11 华南理工大学 基于数字图像处理多孔金属材料浸润性能测量装置及方法
CN105156091B (zh) * 2015-07-13 2018-03-27 中国海洋石油总公司 基于自适应井下钻具工具面动态控制的钻机系统及钻井方法
CN105156021B (zh) * 2015-07-13 2017-07-14 中国海洋石油总公司 基于自适应井下钻具工具面动态控制的钻机系统及钻井方法
WO2017035433A1 (en) * 2015-08-26 2017-03-02 Board Of Regents, The University Of Texas System Systems and methods for measuring relative permeability from unsteady state saturation profiles
CN105158039B (zh) * 2015-09-08 2018-05-08 河海大学 一种用于裂隙岩体渗流试验仿真裂隙制作方法
CN105203439B (zh) * 2015-09-11 2018-01-05 珠海格力电器股份有限公司 空气净化控制方法及装置
US11017289B2 (en) * 2016-05-09 2021-05-25 1Qb Information Technologies Inc. Method and system for improving a policy for a stochastic control problem
WO2018022956A1 (en) * 2016-07-28 2018-02-01 Board Of Regents, The University Of Texas System Systems and methods for determining unsteady-state two-phase relative permeability
CN106353233B (zh) * 2016-08-05 2019-01-11 山东大学 一种岩体粗糙裂隙渗流仿真模型及制作方法
CN108459034A (zh) * 2016-11-18 2018-08-28 中国石油化工股份有限公司 一种砂岩酸岩反应效果可视化定量评价方法
CN106780724A (zh) * 2016-12-21 2017-05-31 中国石油天然气股份有限公司 获取漏失通道模型的方法与装置
CN110573842B (zh) * 2017-01-26 2021-05-11 达索系统西姆利亚公司 基于流体占据时间的多相流可视化图
CN108986167B (zh) * 2017-06-05 2022-01-11 梭特科技股份有限公司 置晶设备的校正方法及使用该方法的置晶设备
JP6869621B2 (ja) * 2017-08-30 2021-05-12 住友重機械工業株式会社 シミュレーション方法及びシミュレーション装置
US11081212B2 (en) 2017-10-06 2021-08-03 DigiM Solution LLC System and method for computing drug controlled release performance using images
US10830713B2 (en) 2017-11-20 2020-11-10 DigiM Solution LLC System and methods for computing physical properties of materials using imaging data
US10426424B2 (en) 2017-11-21 2019-10-01 General Electric Company System and method for generating and performing imaging protocol simulations
US11714040B2 (en) 2018-01-10 2023-08-01 Dassault Systemes Simulia Corp. Determining fluid flow characteristics of porous mediums
US11009497B2 (en) 2018-06-22 2021-05-18 Bp Corporation North America Inc. Systems and methods for estimating mechanical properties of rocks using grain contact models
US11530598B2 (en) 2018-08-21 2022-12-20 Dassault Systemes Simulia Corp. Determination of oil removed by gas via miscible displacement in reservoir rock
WO2020070571A1 (en) * 2018-10-01 2020-04-09 King Abdullah University Of Science And Technology Physics-preserving impes scheme and system
US11314909B2 (en) 2018-10-05 2022-04-26 Qatar Foundation For Education, Science And Community Development Methods and systems for simulating multiphase flow through porous media
US11119025B2 (en) * 2018-12-18 2021-09-14 Shell Oil Company Method for characterizing the porosity of rock
AU2019406627B2 (en) 2018-12-18 2022-05-19 Shell Internationale Research Maatschappij B.V. Method for digitally characterizing the permeability of rock
RU2692369C1 (ru) * 2018-12-26 2019-06-24 Публичное акционерное общество "Газпром нефть" Способ выбора системы разработки месторождения
CN109736790B (zh) * 2018-12-28 2022-05-24 新疆工程学院 基于深度信念网络的测井岩性识别方法
CN109520894B (zh) * 2018-12-29 2020-04-24 西南石油大学 一种页岩有机质强制自吸量预测方法
CN109556996B (zh) * 2019-01-04 2021-06-25 中国石油大学(华东) 油水两相干扰压力梯度的测量方法
US11613984B2 (en) * 2019-09-04 2023-03-28 Dassault Systemes Simulia Corp. Determination of hydrocarbon mobilization potential for enhanced oil recovery
US11875091B2 (en) 2019-09-05 2024-01-16 Toyota Motor Engineering & Manufacturing North America, Inc. Method for data-driven comparison of aerodynamic simulations
WO2021086388A1 (en) 2019-10-31 2021-05-06 Halliburton Energy Services, Inc. Scale-coupled multiscale model simulation
US11079581B2 (en) * 2019-11-25 2021-08-03 Saudi Arabian Oil Company Resolution preserving methodology to generate continuous log scale reservoir permeability profile from petrographic thin section images
CN111027155A (zh) * 2019-12-13 2020-04-17 上海市计量测试技术研究院 基于cfd技术的飞机空气循环制冷系统的模拟分析方法
EP4097447A4 (en) * 2020-03-26 2023-10-11 Hewlett-Packard Development Company, L.P. Fluid dynamics modeling to determine a pore property of a screen device
CN115398216B (zh) 2020-04-22 2025-05-13 国际壳牌研究有限公司 用于估计岩石的流体饱和度的方法
US11847391B2 (en) 2020-06-29 2023-12-19 Dassault Systemes Simulia Corp. Computer system for simulating physical processes using surface algorithm
EP4172661B1 (en) 2020-06-30 2025-04-09 Shell Internationale Research Maatschappij B.V. Method for estimating hydrocarbon saturation of a rock
US12392090B2 (en) * 2020-07-14 2025-08-19 Peridot Print Llc Pore removal from screen devices to increase flow uniformity
CN111879674B (zh) * 2020-07-15 2022-02-01 西南石油大学 基于页岩渗吸渗透率确定合理闷井时间的测试装置及方法
CN112179815B (zh) * 2020-09-21 2022-04-05 西南石油大学 一种基于孔隙网络模型的单相非稳态渗流模型建立方法
US11907625B2 (en) 2020-12-29 2024-02-20 Dassault Systemes Americas Corp. Computer simulation of multi-phase and multi-component fluid flows including physics of under-resolved porous structures
EP4323904A4 (en) * 2021-04-13 2025-01-01 Hewlett-Packard Development Company, L.P. Modification of digital pores in screen digital models
US20240200279A1 (en) * 2021-04-13 2024-06-20 Hewlett-Packard Development Company, L.P. Decreasing fiber density in wet parts formed from slurries
CN115248965B (zh) * 2021-04-25 2025-11-21 广州汽车集团股份有限公司 一种离合器充油时间的计算方法、装置及计算机可读存储介质
CN114235860A (zh) * 2021-11-29 2022-03-25 北京师范大学 Ct扫描技术和离子示踪技术量化土壤中溶质运移的方法
US11746623B2 (en) 2022-01-27 2023-09-05 Halliburton Energy Services, Inc. System and method to calibrate digital rock wettability
CN117548156A (zh) * 2022-08-05 2024-02-13 北京计算科学研究中心 用于多孔介质的流场调节装置及调节方法
US20240053245A1 (en) * 2022-08-12 2024-02-15 Halliburton Energy Services, Inc. Determining relative permeability in a rock sample
AU2023333457A1 (en) * 2022-08-29 2025-03-06 University Of Wyoming Methods and devices for dynamic pore network modeling of two-phase flow in mixed-wet porous media
CN115753543B (zh) * 2022-11-05 2024-01-23 西南石油大学 考虑概率分布的页岩支撑裂缝相对渗透率测定装置及方法
CN115711120B (zh) * 2022-12-02 2024-04-12 西南石油大学 一种油气井水泥环完整性及损伤演化定量评价方法
CN115983161B (zh) * 2023-03-17 2023-06-27 西南石油大学 一种多孔介质应力作用下相对渗透率计算方法
CN116879122B (zh) * 2023-05-31 2024-06-21 中国石油大学(北京) 一种无限可扩展的模块化可视渗流模拟装置
CN116663370B (zh) * 2023-07-05 2023-10-24 成都理工大学 一种深层缝洞型碳酸盐岩油藏流固耦合数值模拟方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2345935A (en) 1942-04-06 1944-04-04 Shell Dev Method and apparatus for permeability measurements
US4506542A (en) 1983-04-22 1985-03-26 Chandler Engineering Company Apparatus and procedure for relative permeability measurements
SU1183868A1 (ru) * 1983-06-09 1985-10-07 Всесоюзный нефтегазовый научно-исследовательский институт Устройство дл измерени водонасыщенности пористых материалов
SU1695176A1 (ru) * 1989-06-07 1991-11-30 Московский Геологоразведочный Институт Им.Серго Орджоникидзе Способ определени относительных фазовых проницаемостей пористой среды
US5497321A (en) * 1994-01-11 1996-03-05 Schlumberger Technology Corporation Well logging method for determining fractional flow characteristics of earth formations
FR2772483B1 (fr) * 1997-12-15 2000-01-14 Inst Francais Du Petrole Methode pour modeliser des deplacements de fluides dans un milieu poreux
US6516080B1 (en) 2000-04-05 2003-02-04 The Board Of Trustees Of The Leland Stanford Junior University Numerical method of estimating physical properties of three-dimensional porous media
EP1260920A1 (en) * 2001-05-23 2002-11-27 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Method and apparatus for computing an interface of a fluid in a space
FR2836228B1 (fr) * 2002-02-21 2005-08-19 Inst Francais Du Petrole Methode et dispositif pour evaluer des parametres physiques d'un gisement souterrain a partir de debris de roche qui y sont preleves
WO2004010234A2 (en) * 2002-07-19 2004-01-29 Celerity Group, Inc. Methods and apparatus for pressure compensation in a mass flow controller
FR2864244B1 (fr) * 2003-12-19 2006-08-04 Inst Francais Du Petrole Methode et dispositif pour mesurer des caracteristiques physiques d'un echantillon solide poreux
US7716028B2 (en) * 2006-05-24 2010-05-11 Schlumberger Technology Corporation Method for modeling a reservoir using a 3D wettability map generated from a wettability logging tool
RU2010152658A (ru) * 2008-05-23 2012-06-27 Те Острейлиан Нэшнл Юниверсити (Au) Обработка данных изображения
US8155377B2 (en) 2008-11-24 2012-04-10 Ingrain, Inc. Method for determining rock physics relationships using computer tomographic images thereof
US8081802B2 (en) 2008-11-29 2011-12-20 Ingrain, Inc. Method for determining permeability of rock formation using computer tomograpic images thereof
CN101556703B (zh) * 2009-05-16 2011-01-26 中国石油大学(华东) 基于连续切片图像的网络模型建立方法
BRPI0924623A2 (pt) * 2009-06-30 2016-03-01 Prad Res & Dev Ltd método numérico para cálculo de transporte de calor de massa, químico e elétrico para um meio tridimensional poroso.
CN201575965U (zh) * 2009-12-25 2010-09-08 北京首瑞测控技术有限公司 混凝土渗透性测试仪
CN102087195B (zh) * 2010-12-03 2012-07-04 宁波大学 全自动微滤膜孔径分布测定仪及其自动测定方法和应用

Also Published As

Publication number Publication date
EP2732135A2 (en) 2014-05-21
CA2840942C (en) 2017-05-02
RU2593853C2 (ru) 2016-08-10
BR112014000758A2 (pt) 2019-10-01
BR112014000758B1 (pt) 2021-01-19
EP2732135B1 (en) 2016-08-24
WO2013009512A3 (en) 2013-11-07
CN103906893A (zh) 2014-07-02
AU2012283030A1 (en) 2014-01-16
US9183326B2 (en) 2015-11-10
MX2014000429A (es) 2014-09-01
CA2840942A1 (en) 2013-01-17
CO6930314A2 (es) 2014-04-28
AU2012283030B2 (en) 2016-01-28
MX350511B (es) 2017-09-04
WO2013009512A2 (en) 2013-01-17
US20130018641A1 (en) 2013-01-17

Similar Documents

Publication Publication Date Title
RU2014104788A (ru) Способ моделирования движения отдельных фаз многофазного/многокомпонентного потока, проходящего через пористую среду
Carrel et al. Biofilms in 3D porous media: Delineating the influence of the pore network geometry, flow and mass transfer on biofilm development
Zhang Pore-scale modelling of relative permeability of cementitious materials using X-ray computed microtomography images
Shmulevich et al. Interaction between soil and a wide cutting blade using the discrete element method
Wang et al. Effects of soil water repellency on infiltration rate and flow instability
Or Scaling of capillary, gravity and viscous forces affecting flow morphology in unsaturated porous media
CN111982740B (zh) 基于长膜泥沙的水体元素迁移扩散的分析方法
Juez et al. Experimental and numerical simulation of bed load transport over steep slopes
Schnoor Modeling chemical transport in lakes, rivers, and estuarine systems
CN114021399B (zh) 一种考虑多个化学反应过程的水盐模拟方法
Pearson et al. Predictive modelling of supercritical CO2 dewatering of capillary tubes
Lu et al. Crossflow microfiltration of mono-dispersed deformable particle suspension
Boyce et al. Intrusion of a Liquid Droplet into a Powder under Gravity
Li et al. A machine learning based-method to generate random circle-packed porous media with the desired porosity and permeability
Chashechkin et al. Unsteady vortex dynamics past a uniformly moving tilted plate
Dushin et al. Two phase flows in porous media under microgravity conditions
Cejas et al. Kinetics of gravity-driven water channels under steady rainfall
Petrolo et al. Experimental study on radial gravity currents flowing in a vegetated channel
Andreasen et al. Relating water and air flow characteristics in coarse granular materials
Werber Experimental Evaluation of Three Backward Transit Time Distributions (bTTD) for Solute Storage and Release During Hyporheic Exchange
Obodovich et al. Mathematical modeling of mass transfer in baromembrane processes
Wandel et al. Spreading rate and dispersion behavior of evaporation-suppressant monolayer on open water surfaces: Part 2–Under wind stress
Van der Meer et al. Numerical simulations of upstream and downstream overdeepening
Iribarne Transport properties of internal gravity waves
Zeidman et al. Coarse-grained Monte Carlo simulation of the distribution and capillarity of multiple fluids in porous media

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20190703