RU2549216C2 - Измерение параметров, связанных с прохождением текучих сред в пористом материале - Google Patents

Измерение параметров, связанных с прохождением текучих сред в пористом материале Download PDF

Info

Publication number
RU2549216C2
RU2549216C2 RU2012136121/28A RU2012136121A RU2549216C2 RU 2549216 C2 RU2549216 C2 RU 2549216C2 RU 2012136121/28 A RU2012136121/28 A RU 2012136121/28A RU 2012136121 A RU2012136121 A RU 2012136121A RU 2549216 C2 RU2549216 C2 RU 2549216C2
Authority
RU
Russia
Prior art keywords
volume
pressure
coefficient
sample
porosity
Prior art date
Application number
RU2012136121/28A
Other languages
English (en)
Other versions
RU2012136121A (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 RU2012136121A publication Critical patent/RU2012136121A/ru
Application granted granted Critical
Publication of RU2549216C2 publication Critical patent/RU2549216C2/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/082Investigating permeability by forcing a fluid through a sample
    • G01N15/0826Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Fluid Mechanics (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

Изобретение относится к измерению физических свойств, связанных с прохождением текучей фазы в пористом материале. Способ оценки физических параметров пористого материала, находящегося в потоке текучих сред, содержит этапы, на которых образец (2) материала помещают в герметичную камеру (1) таким образом, чтобы входная сторона (3) образца сообщалась с первым объемом (V0) и чтобы его выходная сторона (4) сообщалась со вторым объемом. В первом объеме осуществляют модуляцию давления и в течение времени измеряют изменения соответствующих давлений в первом объеме и во втором объеме. При помощи дифференциального уравнения, параметрами которого являются собственная проницаемость материала, его пористость и его коэффициент Клинкенберга, производят цифровой анализ изменений измеряемых давлений для оценки по меньшей мере собственной проницаемости и коэффициента Клинкенберга, а также предпочтительно его пористости в ходе одного эксперимента. Техническим результатом является повышение оценки проницаемости kI и коэффициента Клинкенберга b, а также возможность одновременно производить оценку пористости ϕ в ходе одного эксперимента. 12 з.п. ф-лы, 28 ил., 4 табл.

Description

Изобретение относится к измерению физических свойств, связанных с прохождением текучей фазы в пористом материале.
В частности, изобретение находит свое применение для материалов с очень малыми диаметрами проходных каналов в масштабе поры, то есть для материалов с очень большим сопротивлением прохождению текучей среды (в отличие от собственной проницаемости). В качестве неограничивающих примеров этих материалов можно указать сверхплотные породы газоносных пласт-коллекторов (″tight gas reservoirs″), покрывающие породы в потенциальных местах хранения, материалы, используемые в уплотнительных устройствах, композиционные материалы и т.д.
На уровне характерного блока материала прохождение текучей среды через пористое вещество зависит от трех собственных физических характеристик вещества, которыми являются:
- его жидкостная или собственная проницаемость kI, выраженная в м2 или, что более распространенно, в Д (дарси: 1 Д ≈ 0,987 × 10-12 м2);
- его коэффициент Клинкенберга b, выраженный в Па, если речь идет о малопроницаемом веществе и при прохождении газового потока с низким уровнем давления, или его коэффициент Форхаймера β, выраженный в м-1, называемый также коэффициентом инерционного сопротивления, если речь идет о высоких уровнях напора, приводящих к инерционным эффектам;
- его пористость ϕ, равная отношению объема пустот материала к его общему объему.
Ни один современный метод не позволяет определять эти три параметра одновременно в ходе одного эксперимента. В частности, пористость часто измеряют отдельно от двух других параметров при помощи пикнометрического метода (гелиевая пикнометрия, ртутная пикнометрия,…) или посредством взвешивания.
Для измерения газопроницаемости материала различают методы в стационарном режиме и методы в нестационарном режиме. См., например, J.A. Rushing et al., Klinkenberg-corrected permeability measurements in tight gas sands: Steady-state versus unsteady-state techniques, SPE 89867 1-11, 2004.
Недостатком способа в стационарном режиме является то, что он требует довольно длительного времени для установления стационарного режима потока, чтобы получить точку измерения. Время достижения этого стационарного режима меняется обратно пропорционально значению kI и прямо пропорционально квадрату толщины образца. Оно может легко дойти до нескольких часов при очень низкой проницаемости. Раздельное определение собственной проницаемости kI и коэффициента Клинкенберга b требует нескольких точек измерения и, следовательно, получения такого же числа стационарных состояний. Это может занять много времени, поэтому данный метод не годится для измерения в диапазоне низких значений проницаемости. Кроме того, эта технология требует измерения расхода текучей среды, которое может создать проблемы при очень низкой проницаемости.
Для устранения этих недостатков предпочтительно производить измерение в переходных условиях. Обычно эксперимент при нестационарном состоянии состоит в регистрации изменения дифференциального давления ΔP(t) между концами образца. Каждый конец образца соединяют с соответствующей емкостью и на один из них воздействуют импульсом давления. Этот метод называют ″Pulse decay″. Вариант, в котором выходной резервуар представляет собой бесконечный объем (атмосфера), называется ″Draw down″.
Анализ ΔP(t) дает возможность идентифицировать проницаемость вещества. Часто эта технология не учитывает эффекты Клинкенберга.
В патенте US №2 867 116 был предложен приближенный метод экспериментального определения пористости, эффективной проницаемости (то есть включая эффекты Клинкенберга) и собственной проницаемости. В этих работах значения kI, b и ϕ определяют приближенно, повторяя эксперимент три раза при постоянном соотношении между значением первоначального импульса давления и первоначальным давлением образца. Первый эксперимент проводят, замеряя время уменьшения импульса давления до заданной доли (например, 55%) от его первоначального значения. Второй эксперимент является идентичным первому, но его осуществляют, просто меняя уровень давления импульса и первоначальное давление образца таким образом, чтобы разность между ними была такой же, что и в первом эксперименте. При этом опять замеряют время, чтобы импульс давления уменьшился до такой же доли (55%) от первоначального значения. Третий эксперимент является идентичным двум первым, но с изменением объема камеры, используемой для создания импульса давления. Определение kI, b и ϕ на основании этих трех экспериментов осуществляют приближенно, применяя график и используя эмпирическое линейное поведение. На самом деле в общем случае трудно рассчитать реальное влияние этих приближений. Кроме того, следует отметить сложность эксперимента, связанную с устройством и с временем выполнения, необходимым для доведения образца до кондиции при разных давлениях.
В публикации ″A method for the simultaneous determination of permeability and porosity in a low permeability cores″, SPE 15379, 1-11, 1988, S.E. Haskett et al. был предложен способ определения проницаемости kI, причем без учета эффектов Клинкенберга, и пористости ϕ. Способ требует проведения эксперимента до установления конечного равновесия давления во входном и выходном резервуарах. Он основан на измерении в течение времени разности давления между входным и выходным резервуарами. Эта конфигурация не отличается точностью и не является оптимальной для определения параметров.
В публикации ″A detailed analysis of permeability and Klinkenberg coefficient estimation from unsteady-state pulse decay or draw-down experiments″, Symp. Soc. Core Analysts, Calgary, 10-13 September, 5CA2007-08, 2007, Y. Jannot et al. было предложено пересмотреть метод ″Pulse decay″ без выдвижения каких-либо особых упрощающих гипотез. Просто было указано, что образец представляет собой твердую матрицу, не деформирующуюся под действием потока измерительного газа, и что поток газа является слабосжимаемым, изотермическим и ползучим. В этом контексте физическую задачу, описывающую опыт ″Pulse decay″ в общем случае, можно выразить следующим образом:
x [ ( P + b ) P x ] = φ μ k I P t п р и 0 < x < e и t > 0 ( 1 )
Figure 00000001
при начальных условиях:
P ( 0.0 ) = P 0 i ( 2 )
Figure 00000002
P ( x ,0 ) = P 1 i п р и x > 0 ( 3 )
Figure 00000003
и при граничных условиях:
k I S μ V 0 [ P ( 0, t ) + b ] P x ( 0, t ) = P t ( 0, t ) ( 4 )
Figure 00000004
k I S μ V 1 [ P ( e , t ) + b ] P x ( e , t ) = P t ( e , t ) ( 5 )
Figure 00000005
где: Р является давлением в момент t и в положении x вдоль образца, при этом x = 0 соответствует входной стороне образца, x = е соответствует его выходной стороне, причем импульс давления прикладывают в t = 0;
S является площадью сечения образца;
е является длиной образца;
V0 и V1 соответственно являются объемами входного резервуара (высокое давление) и выходного резервуара (низкое давление), сообщающихся с образцом, в которых первоначально (в t = 0) присутствует давление с соответствующими значениями Р0i и Р1i;
µ является динамической вязкостью газа, предположительно постоянной.
В конфигурации ″Draw down″ второе граничное условие заменено классическим условием Дирихле: P(e,t) = P1 = Р1i. В данном случае предполагают, что первоначально образец находится под действием атмосферного давления, с которым он в нормальных условиях находится в состоянии равновесия.
На входе образца обязательно присутствует мертвый объем, то есть между вентилем, который изолирует образец от входного резервуара, и входной стороной образца. Желательно, чтобы объем V0 был очень малым (в идеале порядка объема поры в образце), чтобы повысить чувствительность измерений к пористости ϕ, но в этом случае ее точное определение с целью учета в условии (4) становится очень деликатным, так как предполагает точное знание мертвого объема. Таким образом, существование этого мертвого объема оказывает большое влияние на оценочные значения kI и b. Кроме того, открывание вентиля в момент начала эксперимента ″Pulse decay″ приводит к расширению текучей среды в мертвом объеме, что влечет за собой термические и гидродинамические помехи, которые, хотя и являются наблюдаемыми, очень трудно учитывать точно при моделировании. Вышеуказанные уравнения (1)-(5) не учитывают эти термические и гидродинамические эффекты.
Ошибка в определении пористости ϕ существенно влияет на оценочные значения проницаемости kI и коэффициента Клинкенберга b. Качественная оценка этих двух параметров требует таким образом точного знания ϕ, если эту величину используют в качестве входного параметра. Применяемые для этого пикнометрические методы занимают много времени и приводят к оценке собственной пористости, а не пористости под напряжением (коэффициент загрузки), обычно применяемой для анализа реального материала.
Существует потребность в экспериментальном методе, позволяющем улучшить оценку проницаемости kI и коэффициента Клинкенберга b (для значений низкой проницаемости, вместо которого определяют коэффициент Форхаймера для значений высокой проницаемости). Кроме того, желательно одновременно производить оценку пористости ϕ в ходе одного эксперимента.
В связи с этим изобретением предложен способ оценки физических параметров материала, содержащий следующие этапы:
- образец материала помещают в герметичную камеру таким образом, чтобы входная сторона образца сообщалась с первым объемом и чтобы выходная сторона образца сообщалась со вторым объемом;
- в первом объеме осуществляют модуляцию давления;
- в течение времени измеряют изменения соответствующих давлений в первом объеме и во втором объеме; и
- при помощи дифференциального уравнения, параметрами которого являются собственная проницаемость материала, пористость материала и, по меньшей мере, один другой коэффициент, характерный для материала, и граничным условием которого является изменение измеряемого давления в первом объеме, производят цифровой анализ изменения измеряемого давления во втором объеме для оценки, по меньшей мере, собственной проницаемости и указанного другого коэффициента.
Чтобы устранить проблемы, связанные с мертвым объемом на входе образца, в качестве исходной данной уже берут не только исходную данную, которую единственно рассматривают как значение импульса давления P0i, которое служит для моделирования изменения P(0,t), для осуществления инверсии. Наоборот, предусматривают резервуар конечного объема V1 с выходной стороны и рассматривают две отдельные данные, причем обе измеряемые: сигнал давления на входе P(0,t) = P0(t) и сигнал давления на выходе P(1,t) = P1(t). Сигнал P0(t) может служить входным сигналом на этапе анализа, который состоит в цифровой инверсии дифференциального уравнения, производимой на выходном сигнале P1(t). Поскольку P0(t) получен в результате измерения, а не моделирования, он может содержать погрешности, связанные с термическими явлениями, с наличием мертвого объема и т.д., хотя это и не является причиной искажения по отношению к модели, используемой в процедуре инверсии.
Другим коэффициентом, характерным для материала и получаемым посредством оценки одновременно с его собственной проницаемостью kI, обычно является коэффициент Клинкенберга b, если известно, что исследуемый материал находится в интервале низких значений проницаемости (примерно менее 10-16 м2). Если проницаемость находится в более высоком интервале, другим коэффициентом может быть коэффициент Форхаймера β. Может существовать диапазон проницаемости, для которого в модели можно учитывать одновременно коэффициент Клинкенберга b и коэффициент Форхаймера β.
В случае, когда вместе с собственной проницаемостью kI определяют коэффициент Клинкенберга b, на этапе анализа производят цифровую инверсию (1) на выходном сигнале P1(t). Граничное условие (4) заменяют условием Дирихле по давлению P(0,t) = P0(t), где P0(t) измеряют при помощи манометра в первом объеме V0. Физическая задача больше не зависит ни от V0, ни от мертвого объема, который таким образом можно игнорировать.
Модуляцию давления в первом объеме применяют не просто моментально, а в шкале времени, большей, чем шкала времени импульса давления. Обычно ее производят в шкале времени, зависящей от диапазона проницаемости материала, но, как правило, в течение более одной минуты. Эту модуляцию давления в первом объеме можно, в частности, осуществлять при помощи последовательности импульсов давления.
В варианте выполнения цифровой анализ изменений измеряемых давлений содержит мониторинг изменения во времени пониженной чувствительности давления P1(t), измеряемого во втором объеме, к собственной проницаемости и изменения во времени пониженной чувствительности P1(t) к коэффициенту Клинкенберга или Форхаймера. Это позволяет убедиться, что модуляция давления была применена в первом объеме таким образом, чтобы избежать стабилизации соотношения между этими двумя чувствительностями, что не позволило бы правильно оценить проницаемость и соответствующий коэффициент.
В предпочтительном варианте выполнения цифровой анализ изменений измеряемых давлений P0(t), P1(t) осуществляют таким образом, чтобы оценить пористость ϕ материала в дополнение к оценке его собственной проницаемости kI и его коэффициента Клинкенберга b (или Форхаймера β).
В классическом эксперименте типа ″Pulse decay″ чувствительность P1(t) быстро становится постоянной, причем по истечении достаточно короткого времени, чтобы этот параметр можно было оценить правильно. Чтобы повысить эту чувствительность, следует увеличить число краткосрочных эффектов, чтобы накапливание текучей среды в порах материала происходило многократно в течение эксперимента. Поскольку способ содержит измерение P0(t), которое становится данной для осуществления инверсии на P1(t), возможно любое задаваемое изменение P0(t). Таким образом, на входе образца генерируют последовательность импульсов давления, чтобы активировать емкостное поведение системы с целью облегчения оценки пористости.
Цифровой анализ изменений измеряемых давлений может содержать мониторинг изменения во времени пониженной чувствительности измеряемого давления P1(t) к пористости. Это позволяет убедиться, что модуляция давления была применена в первом объеме таким образом, чтобы избежать стабилизации этой пониженной чувствительности к пористости, что не позволило бы правильно оценить пористость ϕ.
Для облегчения схождения оценки параметров в некоторых случаях можно производить предварительную оценку собственной проницаемости kI и коэффициента Клинкенберга b при помощи давлений, измеренных в интервалах времени, в которых давление во втором объеме меняется по существу линейно.
В этом случае предпочтительный вариант выполнения содержит исследование изменения во времени давления во втором объеме. Если это исследование показывает, что давление во втором объеме меняется по существу линейно в зависимости от времени, это давление оставляют изменяться по существу линейно, чтобы получить значения для предварительной оценки собственной проницаемости и коэффициента, после чего в первом объеме применяют новый импульс давления.
Другие отличительные признаки и преимущества настоящего изобретения будут более очевидны из нижеследующего описания неограничивающего примера осуществления со ссылками на прилагаемые чертежи, на которых:
на фиг.1 показана схема установки, используемой для применения заявленного способа оценки физических параметров;
на фиг.2 представлен график, показывающий значения пониженной чувствительности к проницаемости, к коэффициенту Клинкенберга и к пористости в варианте осуществления способа;
на фиг.3 представлен график, показывающий изменение моделируемого давления на выходе образца в примере применения способа;
на фиг.4 представлен график, показывающий изменение значений пониженной чувствительности к проницаемости, к коэффициенту Клинкенберга и к пористости в примере, представленном на фиг.3;
на фиг.5 представлен график, показывающий изменение соотношения между значениями пониженной чувствительности к проницаемости и к коэффициенту Клинкенберга в примере, представленном на фиг.3;
на фиг.6 представлен график, показывающий изменение соотношения между значениями пониженной чувствительности к проницаемости и к пористости в примере, представленном на фиг.3;
на фиг.7-10 представлены графики, аналогичные графикам, показанным на фиг.3-6, в другом примере применения способа;
на фиг.11-14 представлены графики, аналогичные графикам, показанным на фиг.3-6, в еще одном примере применения способа;
на фиг.15 и 16 представлены графики, показывающие изменение моделируемых давлений на входе и на выходе образца в случае теста, согласно способу;
на фиг.17 и 18 представлены графики, показывающие изменение моделируемых давлений на входе и на выходе образца в испытании на сосновом образце;
на фиг.19 представлен график, показывающий остаточное давление на выходе образца в испытании, изображенном на фиг.17 и 18, при этом остаток представляет собой разность между давлением, вычисленным при помощи модели, описывающей физику испытания, и давлением, измеренным во время испытания;
на фиг.20-22 представлены графики, аналогичные графикам, показанным на фиг.17-19, во время первого испытания на образце породы;
на фиг.23-25 представлены графики, аналогичные графикам, показанным на фиг.17-19, во время второго испытания на этом же образце породы;
на фиг.26-28 представлены графики, аналогичные графикам, показанным на фиг.17-19, во время третьего испытания на этом же образце породы.
Установка, показанная на фиг.1, содержит камеру Хасслера, в которую помещают образец 2 материала, для которого необходимо определить физические параметры в присутствии потока текучей среды. В качестве неограничивающего примера текучей среды можно указать газ, такой как азот или гелий.
Как известно, камера Хасслера представляет собой полый цилиндр, в котором образец 2 цилиндрической формы и с сечением S закрывают с уплотнением, чтобы заставить газ проходить через пористую структуру материала. Образец 2 имеет входную сторону 3 и выходную сторону 4, которые сообщаются с двумя резервуарами 5, 6, объемы которых обозначены соответственно V0 и V1.
Манометры 7, 8 позволяют измерять давления в резервуарах 5, 6. Газ, который пропускают через образец, поступает из баллона 10, соединенного с входным объемом V0 через вентиль 11 и редуктор 12. Со стороны выхода объем V1 соединен с баллоном-коллектором 15 через вентиль 16 и редуктор 17. Между редуктором 12 и входным объемом V0 и между редуктором 17 и выходным объемом V1 установлены другие вентили 18, 19, чтобы выборочно устанавливать сообщение между редукторами и резервуарами 5, 6.
Другой вентиль 20 установлен между входным резервуаром 5 и камерой Хасслера 1 для создания импульсов давления на уровне входной стороны 3 образца. Для приложения первого импульса давления к образцу 2 вентиль 19 позиционируют таким образом, чтобы установить выходной резервуар 6 на исходное давление Р1i (например, атмосферное давление), тогда как вентиль 20 закрыт. После достижения равновесия давления вентиль 19 закрывают. Открывают вентили 11, 18 и регулируют редуктор 12 по значению давления, необходимому для импульса. Таким образом, газ заполняет входной объем V0 при необходимом давлении. После этого закрывают вентиль 18 и открывают вентиль 20 для приложения импульса давления к образцу 2. При помощи манометров 7 и 8 отслеживают снижение давления во входном объеме V0 и повышение давления в выходном объеме V1. Изменение измеряемых давлений P0(t) и P1(t) записывают для последующего цифрового анализа. Чтобы в дальнейшем приложить еще один импульс давления к образцу 2, редуктор 12 регулируют по новому необходимому значению давления, открывают вентиль 18 для заполнения объема Vo при необходимом давлении, затем его закрывают.
Перед приложением первого импульса давления вентиль 20 закрыт, образец 2 находится в состоянии равновесия с выходным объемом V1 и можно применить первоначальное условие (3). В случае когда эффектами Форхаймера можно пренебречь, физической задачей, решаемой в цифровом виде для оценки параметров, является следующая задача (1)-(3)-(4′)-(5):
x [ ( P + b ) P x ] = φ μ k I P t п р и 0 < x < e и t > 0 ( 1 )
Figure 00000001
при первоначальном условии:
P ( x ,0 ) = P 1 i п р и x > 0 ( 3 )
Figure 00000003
и при граничных условиях:
P ( 0, t ) = P 0 ( t ) п р и t 0 ( 4 ' )
Figure 00000006
k I S μ V 1 [ P ( e , t ) + b ] P x ( e , t ) = P t ( e , t ) ( 5 )
Figure 00000005
В выражении этой задачи давление P0(t) на входе образца 2 является данной. Физическими параметрами материала образца 2, участвующими в системе, являются его пористость ϕ, его собственная проницаемость kI и его коэффициент Клинкенберга b.
Возможность оценки параметров на основании сигнала f(t) можно определить посредством исследования чувствительности. В нашем случае можно, например, рассмотреть сигнал f(t) - P1(t). Чувствительность f(t) к параметру ψ, который необходимо оценить, определяется отношением ( t ) ψ
Figure 00000007
. Из практических соображений чаще всего применяют Σ ψ = ψ ( t ) ψ
Figure 00000008
, что позволяет получить эти величины в единицах давления. Анализ изменения этих величин в зависимости от времени позволяет спрогнозировать возможность оценки параметра ψ на основании сигнала f(t). Эта оценка оказывается возможной, если:
- изменения Σ ψ = ψ ( t ) ψ
Figure 00000008
являются значительными в достаточно большом интервале времени перед шагом временной дискретизации сигнала. Понятие «значительные» означает, что Σψ превышает точность измерительного инструмента (датчики давления 7, 8), используемого для измерения f(t);
- если искомыми являются несколько параметров (например, kI, b и даже ϕ), значения пониженной чувствительности ко всем этим параметрам необходимо вывести из взаимосвязи, то есть они не должны быть пропорциональны между собой. В противном случае это значит, что изменения, наблюдаемые на f(t), нельзя по отдельности отнести к тому или иному из параметров, что делает их одновременную оценку невозможной на основании только одного сигнала f(t).
На фиг.2 показано изменение во времени значений пониженной чувствительности к проницаемости kI, к коэффициенту Клинкенберга b и к пористости ϕ, вычисленных в случае только одного импульса давления (типа ″Pulse decay″) в следующих условиях: kI = 10-19 м2, b = 13,08 бар, ϕ = 0,02, е = 5 см, диаметр образца d = 5 см, V0 = 10-3 м3, V1 = 2,5.10-3 м3 при исходном давлении 15 бар во входном объеме V0 и 1 бар в выходном объеме V1.
Эти значения чувствительности вычисляют на основании сигналов P1(t), смоделированных при помощи физической модели (1)-(3)-(4′)-(5). При этом по истечении нескольких десятков минут пониженная чувствительность Σϕ к пористости стабилизируется, поэтому после этого времени измерение давления больше не дает представления о пористости ϕ. Таким образом, измерения, произведенные в условиях, представленных на фиг.2, могут оказаться недостаточными для определения пористости ϕ. Тем не менее, они могут быть использованы для определения проницаемости kI и коэффициента Клинкенберга b, если значение пористости ϕ известно. Эти оценки kI и b можно получить, не прибегая к точному определению объема V0 и соответствующего мертвого объема с входной стороны образца 2 и не обращая внимания на возможные погрешности P0(t), которые после измерения уже не надо вычислять.
Для повышения чувствительности к пористости ϕ и для обеспечения ее оценки следует увеличить число краткосрочных эффектов, чтобы накапливание газа в порах материала происходило многократно в течение всего времени эксперимента. Это проиллюстрировано ниже при помощи трех примеров.
Пример 1 (фиг.3-6)
В этом примере анализ чувствительности был произведен посредством моделирования на материале, имеющем собственную проницаемость kI = 10-17 м2, коэффициент Клинкенберга b = 2,49 бар, пористость ϕ = 0,02 в течение времени эксперимента tf = 500 с. Были применены три последовательных импульса давления, первый - 5 бар при t = 0, второй - 10 бар при t = tf/3 и третий - 15 бар при t = 2tf/3.
Объем входного резервуара 5 был равен V0 = 10-3 м3, а объем выходного резервуара 6 был равен V1 = 2,5.10-3 м3.
На фиг.3 показано изменение во времени давления P1(t) на выходе образца. На фиг.4 показано изменение во времени значений пониженной чувствительности ΣkI, Σb и Σϕ давления P1(t) к собственной проницаемости kI, к коэффициенту Клинкенберга b и к пористости ϕ. На фиг.5 показано изменение во времени соотношения между значениями пониженной чувствительности ΣkI, Σb к собственной проницаемости kI и к коэффициенту Клинкенберга b. На фиг.6 показано изменение во времени соотношения между значениями пониженной чувствительности ΣkI, Σϕ к собственной проницаемости kI и к пористости ϕ.
Пример 2 (фиг.7-10)
В этом примере анализ чувствительности был произведен в тех же условиях, что и в примере 1, на материале, имеющем собственную проницаемость kI = 10-17 м2, коэффициент Клинкенберга b = 2,49 бар, пористость ϕ = 0,1 в течение времени эксперимента tf = 200 с. На фиг.7-10 представлены графики для примера 2, аналогичные графикам на фиг.3-6.
Пример 3 (фиг.11-14)
В этом примере анализ чувствительности был произведен в тех же условиях, что и в примерах 1 и 2, на материале, имеющем собственную проницаемость kI = 10-19 м2, коэффициент Клинкенберга b = 13,08 бар, пористость ϕ = 0,02 в течение времени эксперимента tf = 13000 с. На фиг.11-14 представлены графики для примера 3, аналогичные графикам на фиг.3-6.
Все три примера показывают на трех разных материалах, что повышения давления P1(t) на выходе образца являются измеряемыми величинами, даже если выбрать относительно большой объем V1 (2,5 литра). Этот объем был специально выбран большим, чтобы подчеркнуть, как можно минимизировать относительную ошибку при измерении. Выбор меньшего объема приводит к более значительным повышениям и можно убедиться, что это не влияет на значения чувствительности. Выбор большого объема для V1 не мотивирован, как в методе ″Draw down″, получением большого изменения среднего давления в течение времени. Действительно, очень значительное изменение среднего давления в данном случае вызвано последовательными импульсами давления.
Во всех случаях значения чувствительности являются очень значительными и хорошо декоррелированы по парам. Это позволяет произвести одновременную оценку трех параметров kI, b и ϕ. Если сравнить фиг.4, 8 или 12 с фиг.2, то можно заметить, что модуляция входного давления при помощи нескольких последовательных импульсов существенно улучшает чувствительность к пористости ϕ, что облегчает ее оценку.
Чтобы проиллюстрировать возможность одновременной оценки трех параметров kI, b и ϕ при помощи заявленного способа, был проведен ряд тестов на основании сигналов, сформированных в цифровом виде при помощи физической модели, при P0i=1бар. Чтобы лучше представить реальное измерение, на два смоделированных сигнала P0(t) и P1(t) был наложен произвольный шум посредством δР0 = 0,01 × dP × s × Р0max/3 и δP1 = 0,01 × dP × s × P1max/3. Эти шумы таковы, что s является случайным числом со стандартным отклонением от 1, а dP является ошибкой от P0(t) и P1(t) (взятые в % от измерений). Коэффициент 3 был взят таким образом, чтобы интервалы P0(t) ± δР0 и P1(t) ± δP1 включали в себя 99,7% значений, если бы они были действительно измерены. В этих моделях взяли dP = 0,1%, что является типовым значением для датчика давления, исключая случаи 14, 15 и 16, для которых взяли dP = 1%.
Таблица I
Случай № kI, (м2) b (бар) ϕ (%) V03) V13) P0i (бар) dP (%) N tf (с) М
1 10-17 2,49 0,02 10-5 10-2 3 0,1 1000 1000 200
2 10-17 2,49 0,02 10-2 10-2 3 0,1 1000 1000 200
3 10-17 2,49 0,02 10-5 10-2 3 0,1 100 1000 200
4 10-17 2,49 0,02 10-5 10-2 3 0,1 1000 300 200
5 10-17 2,49 0,02 10-5 10-2 3 0,1 1000 1000 50
6 10-17 2,49 0,1 10-5 10-2 3 0,1 1000 1000 200
7 10-17 13,08 0,02 10-5 10-2 3 0,1 1000 100000 200
8 10-17 2,49 0,02 10-2 10-4 3 0,1 1000 1000 200
9 10-17 2,49 0,1 10-2 10-4 3 0,1 1000 1000 200
10 10-17 13,08 0,02 10-2 10-4 3 0,1 1000 10000 200
11 10-17 13,08 0,02 10-2 10-4 5 0,1 1000 10000 200
12 10-17 2,49 0,02 10-2 10-4 3 1 1000 1000 200
13 10-17 2,49 0,1 10-2 10-4 3 1 1000 1000 200
14 10-17 13,08 0,02 10-2 10-4 5 1 1000 10000 200
Параметрами, применяемыми в последовательности тестов, являются параметры, указанные в таблице I, в том числе число N точек измерения давлений P0(t) и P1(t), продолжительность tf эксперимента и число М шагов дискретизации в пространстве по толщине е образца, используемых для инверсии задачи (1)-(3)-(4′)-(5). В каждом случае в моменты 0, tf/3 и 2tf/3 были применены три импульса давления, доводящие давление во входном резервуаре до P0i, 2P0i и 3P0i. Модуляция давления, примененная в этой последовательности тестов, позволяет в этом случае квалифицировать метод измерения как ″Step Decay″. На фиг.15 и 16 показано изменение давлений P0(t) и P1(t) в барах, синтезированных на входе и на выходе образца в случае №1.
Результаты инверсий представлены в таблице II с отклонениями dkI, db, dϕ в % между исходными значениями kI, b и ϕ и значениями, полученными посредством инверсионных вычислений.
Таблица II
Случай № kI2) dkI (%) b (бар) db (%) ϕ (%) dϕ (%)
1 1,007.10-17 0,10 2,456 0,26 1,96 0,28
2 1,003.10-17 0,11 2,469 0,34 1,97 0,50
3 1,067.10-17 0,61 2,159 1,70 1,60 2,00
4 1.005.10-17 0,09 2,460- 0,31 1,98 0,21
5 1,006.10-17 0,10 2,460 0,26 1,97 0,28
6 1,007.10-17 0,09 2,454 0,28 9,93 0,18
7 1,030.10-1У 0,31 12,55 0,44 1,96 0,25
8 1,004.10-17 0,11 2,460 0,35 1,95 0,53
9 1,002.10-17 0,08 2,472 0,38 9,95 0,38
10 1,027.10-17 0,21 12,58 0,47 1,96 0,30
11 1,007.10-17 0,21 13,08 0,37 2,00 0,30
12 1,001.10-17 1,07 2,49 3,47 2,06 5,11
13 1,003.10-17 0,81 2,42 3,85 9,67 3,81
14 1,088.10-17 1,81 11,24 3,37 1,78 2,99
Эти результаты позволяют сделать следующие выводы:
- точность оценок является высокой (часто лучше чем 1%) для трех оценочных параметров и остается вполне приемлемой с измерительным шумом в 1% от максимального значения давления (случаи №№12-14);
- точность мало зависит от объема Vo. Вполне приемлемым является объем от 0,1 до 10 литров, предпочтительно от 0,1 до 1 литра;
- объем V1 необходимо выбирать таким образом, чтобы повышение давления было достаточным, чтобы его можно было измерить с высокой точностью. Значение 0,1 литра приводит к удовлетворительным результатам в случае рассматриваемых материалов и является достаточно высоким, чтобы ограничить погрешность, связанную с мертвым объемом. В общем случае можно выбрать объем V1 от 0,05 до 10 литров;
- работа с высокими уровнями давления (5, 10 и 15 бар) приводит к более высокой точности в случае, когда b = 13,08 бар;
- число 1000 экспериментальных точек представляется адекватным, при этом точность в определенной степени понижается, если это число уменьшить до 100 (случай №3);
- увеличение продолжительности эксперимента сверх определенного предела существенно не повышает точность (сравнение случаев №7 и 10);
- приемлемой продолжительностью для оценки трех параметров является 20 минут для kI около 10-17 м и 3 часа для kI около 10-19 м2. В целом следует применять модуляцию давления в первом объеме в шкале времени от нескольких десятков минут до нескольких часов, но в любом случае более одной минуты. Следует отметить, что модуляция давления при помощи последовательности импульсов является частным случаем и что можно применять другие варианты модуляции в достаточной шкале времени для применения изобретения, учитывая, что измеряемый профиль давления на входе P0(t) может иметь любую форму.
Преимуществом относительно большого объема V0 является то, что между двумя импульсами давление P0(t) меняется лишь незначительно. Кроме того, если объем V1 выбирают достаточно большим, чтобы повышение давления было мало значительным по сравнению с Р0, эксперимент проходит в условиях слабонестационарного (квазистационарного) состояния. В этих условиях можно осуществить предварительную оценку хорошего качества на участках P1(t), соответствующих каждому импульсу давления. Эта предварительная качественная оценка позволяет обеспечить более легкое схождение оценки с полной моделью по всему сигналу.
Можно применять автоматизацию проведения эксперимента. Действительно, появление линейного или квазилинейного режима на P1(t) при каждом импульсе давления (фиг.16) соответствует квазистационарному режиму с потерей чувствительности P1(t) к пористости ϕ (по определению, в квазистационарном режиме эффект накапливания в порах образца исчезает). Следовательно, эксперимент можно проводить таким образом, чтобы каждый импульс давления имел продолжительность, позволяющую P1(t) достичь квазилинейного поведения в зависимости от времени. Этот квазилинейный режим оставляют на короткий период, позволяющий произвести хорошую предварительную оценку kI и b. Таким образом, использование достаточно больших объемов V0 и V1 позволяет управлять экспериментом напрямую, чтобы оптимизировать его продолжительность и получать при этом совпадающие результаты.
Были проведены лабораторные испытания при V0 = 1,02.10-3 м и V1 = 2,26.10-3 м3 согласно следующему экспериментальному протоколу:
- образец 2 поместили в камеру Хасслера 1;
- в камере, внешней по отношению к камере Хасслера, создали давление;
- закрыли вентиль 20 и дождались равновесия;
- открыли вентиль 18 и отрегулировали редуктор 12 для получения Р0 = Р0i1;
- начали регистрацию давлений P0(t) и P1(t);
- закрыли вентиль 18 и открыли вентиль 20;
- отрегулировали редуктор 12 для получения Р0 = P0i2;
- после времени t1 на несколько секунд открыли вентиль 18;
- отрегулировали редуктор 12 для получения Р0 = P0i3;
- после времени t2 на несколько секунд открыли вентиль 18;
- после времени t3 остановили регистрацию, открыли вентиль 19 и извлекли образец 2.
Для некоторых из этих испытаний была произведена предварительная оценка проницаемости kI и коэффициента Клинкенберга b следующим образом:
- произвели оценку крутизны трех частей кривых P1(t), соотнесенных с прямыми, соответствующими последовательным импульсам давления;
- на основании этого вывели три значения эффективной проницаемости kg. Для этих вычислений использовали значения Р0 и P1, равные половине суммы значений на концах каждого интервала;
- построили график значений kg в зависимости от 1/Рсред = 2/(Р0+P1) и посредством линейной регрессии классически получили предварительную оценку kI и b с учетом того, что kg = kI. (1+b/Рсред).
Эти предварительные оценочные значения kI и b в дальнейшем использовали в качестве исходных данных для конечной оценки kI, b и ϕ путем инверсии на полном сигнале P1(t), при этом сигнал P0(t) является входной данной, при этом использовали модель (1)-(3)-(4′)-(5).
Пример 4 (Фиг.17-19)
Согласно вышеуказанному протоколу эксперимента провели два испытания на деревянном образце из высушенной сосны, имеющем следующие размерные характеристики: d-38,5 мм и е = 60 мм. Пористость (без напряжений) образца измерили посредством пикнометрии и она составила ϕ = 0,27.
Во втором испытании проницаемость kI и коэффициент Клинкенберга b предварительно оценили в значениях 1,76·10-16 м и 0,099 бар. Конечные результаты оценки приведены в таблице III с типовыми относительными отклонениями σkI, σb и σϕ на трех параметрах, определенных путем одновременной оценки.
Таблица III
Испытание № kI, (м2) σkI (%) b (бар) σb (%) ϕ σϕ (%)
1 1,64.10'1b 0,045 0,230 0,98 0,257 0,38
2 1,70.10"1b 0,038 0,187 1Д4 0,252 0,39
Изменения измеренных давлений P0(t) в барах и ΔP1(t) = P1(t) - P1(0) в миллибарах представлены на фиг.17-18. На фиг.19 показан остаток P1(t) в миллибарах после оценки. Оценки имеют отличное качество, о чем свидетельствуют кривые измеренного и оценочного P1(t) и особенно кривая остатка, что подтверждают незначительные относительные типовые отклонения σkI, σb, σϕ.
Пример 5 (фиг.20-28)
Согласно вышеуказанному протоколу эксперимента провели три испытания на керне породы, имеющем следующие размерные характеристики: d=38 мм и е = 60,3 мм. Пористость (без напряжений) керна измерили посредством пикнометрии и она составила ϕ = 0,06.
Предварительная оценка проницаемости kI и коэффициента Клинкенберга b дала значения 3,34·10-17 м2 и 1,47 бар во втором испытании и 3,86·10-17 м2 и 0,97 бар в третьем испытании. Конечные результаты оценки приведены в таблице IV.
Таблица IV
Испытание № kI2) σkI (%) b (бар) σb (%) ϕ σϕ (%)
1 3,41.10-17 0,058 1,57 0,18 0,055 0,23
2 3,51.10-17 0,023 1,41 0,11 0,055 0,28
3 3,77.10-17 0,052 1,16 0,41 0,048 1,13
Изменения измеренных давлений P0(t) в барах и ΔP1(t) = P1(t) - P1(0) в миллибарах представлены на фиг.20-21 для первого испытания, на фиг.23-24 для второго испытания и на фиг.26-27 для третьего испытания. Остатки по P1(t) в миллибарах после оценки показаны на фиг.22 для первого испытания, на фиг.25 для второго испытания и на фиг.28 для третьего испытания. В данном случае оценки тоже имеют отличное качество, о чем свидетельствуют кривые измеренных и оценочных P1(t) и кривые остатка и что подтверждают незначительные относительные типовые отклонения σkI, σb, σϕ.

Claims (13)

1. Способ оценки физических параметров пористого материала, находящегося в потоке текучей среды, содержащий этапы, на которых:
- помещают образец (2) материала в герметичную камеру (1) так, чтобы входная сторона (3) образца сообщалась с первым объемом (V0), а выходная сторона (4) образца сообщалась со вторым объемом (V1);
- осуществляют модуляцию давления в первом объеме;
- измеряют изменения по времени соответствующих давлений (P0(t), P1(t)) в первом объеме и во втором объеме; и
- при помощи дифференциального уравнения, параметрами которого являются собственная проницаемость (kI) материала, пористость (ϕ) материала и по меньшей мере один другой коэффициент (b, β), характерный для материала, и граничным условием которого является измеряемое изменение давления (P0(t)) в первом объеме, производят цифровой анализ измеряемого изменения давления (P1(t)) во втором объеме для оценки по меньшей мере собственной проницаемости и указанного другого коэффициента.
2. Способ по п.1, в котором модуляцию давления в первом объеме (V0) применяют в шкале времени, большей, чем шкала времени импульса давления.
3. Способ по п.1, в котором модуляцию давления в первом объеме (V0) производят в шкале времени, превышающей одну минуту.
4. Способ по п.1, в котором модуляцию давления в первом объеме (V0) вызывают последовательностью импульсов давления.
5. Способ по п.1, в котором цифровой анализ изменений измеряемых давлений (P0(t), P1(t)) содержит отслеживание изменения во времени пониженной чувствительности (ΣkI) давления P1(t), измеряемого во втором объеме (V1), к собственной проницаемости (kI) и изменения во времени пониженной чувствительности (Σb) давления P1(t), измеряемого во втором объеме, к указанному другому коэффициенту (b, β).
6. Способ по п.1, в котором цифровой анализ изменений измеряемых давлений (P0(t), P1(t)) осуществляют так, чтобы дополнительно оценить пористость (ϕ) материала.
7. Способ по п.6, в котором цифровой анализ изменений измеряемых давлений (P0(t), P1(t)) содержит отслеживание изменения во времени пониженной чувствительности (Σϕ) давления (P1(t)), измеряемого во втором объеме (V1), к пористости (ϕ).
8. Способ по п.1, в котором цифровой анализ измеряемых давлений (P0(t), P1(t)) содержит в интервалах времени, в которых давление (P1(t)) во втором объеме меняется по существу линейно, предварительную оценку собственной проницаемости (kI) и указанного коэффициента (b) для облегчения схождения оценки.
9. Способ по п.8, в котором исследуют изменения во времени давления (P1(t)) во втором объеме (V1) и, если это исследование показывает, что давление во втором объеме меняется по существу линейно в зависимости от времени, это давление оставляют изменяться по существу линейно, чтобы получить значения для предварительной оценки собственной проницаемости (kI) и указанного коэффициента (b), затем в первом объеме (V0) применяют новый импульс давления.
10. Способ по п.1, в котором первый объем (V0) составляет от 0,1 до 10 литров.
11. Способ по п.1, в котором второй объем (V1) составляет от 0,05 до 10 литров.
12. Способ по п.1, в котором указанный другой коэффициент является коэффициентом Клинкенберга (b).
13. Способ по п.1, в котором указанный другой коэффициент является коэффициентом Форхаймера (β).
RU2012136121/28A 2010-01-22 2011-01-21 Измерение параметров, связанных с прохождением текучих сред в пористом материале RU2549216C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1050437 2010-01-22
FR1050437A FR2955662B1 (fr) 2010-01-22 2010-01-22 Mesure de parametres lies a l'ecoulement de fluides dans un materiau poreux
PCT/FR2011/050123 WO2011089367A1 (fr) 2010-01-22 2011-01-21 Mesure de parametres lies a l'ecoulement de fluides dans un materiau poreux

Publications (2)

Publication Number Publication Date
RU2012136121A RU2012136121A (ru) 2014-02-27
RU2549216C2 true RU2549216C2 (ru) 2015-04-20

Family

ID=42671805

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2012136121/28A RU2549216C2 (ru) 2010-01-22 2011-01-21 Измерение параметров, связанных с прохождением текучих сред в пористом материале

Country Status (9)

Country Link
US (1) US20130054157A1 (ru)
EP (1) EP2526400A1 (ru)
CN (1) CN102906556A (ru)
AU (1) AU2011208574C1 (ru)
BR (1) BR112012018093A2 (ru)
CA (1) CA2787771A1 (ru)
FR (1) FR2955662B1 (ru)
RU (1) RU2549216C2 (ru)
WO (1) WO2011089367A1 (ru)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2982949A1 (fr) * 2011-11-23 2013-05-24 Diam Bouchage Dispositif pour la mesure de la permeabilite de bouchons de bouteilles et methode correspondante
FR3002632B1 (fr) * 2013-02-27 2020-08-07 Brgm Dispositif d'analyse en milieu percolant
PT107408B (pt) * 2014-01-17 2021-01-21 Amorim Cork Research & Services, Lda. Processo e dispositivo para verificação de estanquidade de rolhas de cortiça
WO2015160691A1 (en) * 2014-04-14 2015-10-22 Schlumberger Canada Limited Methods for measurement of ultra-low permeability and porosity
US10288517B2 (en) 2014-04-14 2019-05-14 Schlumberger Technology Corporation Apparatus and calibration method for measurement of ultra-low permeability and porosity
US10274411B2 (en) 2014-04-14 2019-04-30 Schlumberger Technology Corporation Methods for measurement of ultra-low permeability and porosity
US10108762B2 (en) * 2014-10-03 2018-10-23 International Business Machines Corporation Tunable miniaturized physical subsurface model for simulation and inversion
WO2016179593A1 (en) * 2015-05-07 2016-11-10 The Uab Research Foundation Full immersion pressure-pulse decay
US10365202B2 (en) * 2015-05-11 2019-07-30 Schlumberger Technology Corporation Method for measurement of ultra-low permeability and porosity by accounting for adsorption
CN104990857A (zh) * 2015-07-23 2015-10-21 重庆大学 真三轴环境下煤岩渗透率的检测方法及装置
CN107907661A (zh) * 2017-12-15 2018-04-13 中国科学院地质与地球物理研究所兰州油气资源研究中心 盆地深部储层岩石与流体相互作用模拟装置及使用方法
CN108801872B (zh) * 2018-04-18 2021-02-12 中国矿业大学 一种岩土材料渗流系数相关偏度确定方法
WO2019236466A1 (en) 2018-06-05 2019-12-12 Saudi Arabian Oil Company Systems and methods for analyzing natural gas flow in subterranean reservoirs
CN109975140B (zh) * 2019-04-16 2022-02-22 重庆地质矿产研究院 超临界二氧化碳脉冲致裂与渗透率测试一体化的实验装置及方法
US11079313B2 (en) * 2019-05-17 2021-08-03 Saudi Arabian Oil Company Methods and systems for determining core permeability pulse decay experiments
CN110320149B (zh) * 2019-08-02 2021-08-03 西南石油大学 一种流向可调式不规则岩样高压渗透装置及测试方法
FR3111706B1 (fr) 2020-06-19 2022-06-03 Ifp Energies Now Procédé pour déterminer le volume poreux d'un échantillon de milieu poreux
CN112362556B (zh) * 2020-11-13 2024-03-29 重庆大学 获得煤矿采动稳定区渗透系数连续函数的方法
CN112557276B (zh) * 2020-11-26 2021-09-10 清华大学 一种同时测量多孔介质渗透率和孔隙率的方法
CN112924357B (zh) * 2021-01-29 2022-02-01 西南石油大学 一种地层压力下致密岩石孔渗联测装置及方法
CN115096784B (zh) * 2022-06-01 2024-04-19 湖南大学 基于Forchheimer定理对高渗混凝土非线性低流速时渗透系数的测定

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2166747C1 (ru) * 2000-04-13 2001-05-10 Уфимский государственный нефтяной технический университет Устройство для определения распределения пор по размерам
FR2863052A1 (fr) * 2003-12-02 2005-06-03 Inst Francais Du Petrole Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2867116A (en) 1954-12-27 1959-01-06 Socony Mobil Oil Co Inc Method of measuring characteristics of porous material
US5417104A (en) * 1993-05-28 1995-05-23 Gas Research Institute Determination of permeability of porous media by streaming potential and electro-osmotic coefficients
US5442950A (en) * 1993-10-18 1995-08-22 Saudi Arabian Oil Company Method and apparatus for determining properties of reservoir rock
US5513515A (en) * 1995-05-15 1996-05-07 Modern Controls, Inc. Method for measuring permeability of a material
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
FR2853071B1 (fr) * 2003-03-26 2005-05-06 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

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2166747C1 (ru) * 2000-04-13 2001-05-10 Уфимский государственный нефтяной технический университет Устройство для определения распределения пор по размерам
FR2863052A1 (fr) * 2003-12-02 2005-06-03 Inst Francais Du Petrole Methode pour determiner les composantes d'un tenseur de permeabilite effectif d'une roche poreuse

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Stefan Finsterle et al, Determining permeability of tight rock samples using inverse modeling, Water Resources Research, Volume 33, стр. 1803"1811, August 1997. *

Also Published As

Publication number Publication date
BR112012018093A2 (pt) 2016-05-03
AU2011208574B2 (en) 2014-03-27
EP2526400A1 (fr) 2012-11-28
FR2955662B1 (fr) 2014-08-22
RU2012136121A (ru) 2014-02-27
CA2787771A1 (fr) 2011-07-28
CN102906556A (zh) 2013-01-30
WO2011089367A1 (fr) 2011-07-28
AU2011208574A1 (en) 2012-08-30
AU2011208574C1 (en) 2014-07-31
FR2955662A1 (fr) 2011-07-29
US20130054157A1 (en) 2013-02-28

Similar Documents

Publication Publication Date Title
RU2549216C2 (ru) Измерение параметров, связанных с прохождением текучих сред в пористом материале
US10571384B2 (en) Methods and systems for determining gas permeability of a subsurface formation
US10416064B2 (en) Methods and systems for determining gas permeability of a subsurface formation
CN103528934B (zh) 测量超低渗岩石渗透率应力敏感性的互相关方法
CN104237099B (zh) 测定致密岩心径向渗透率的装置及方法
US11371905B2 (en) Methods for detecting leakage in permeability measurement system
CN103335928B (zh) 一种测量孔隙岩石渗透率的方法和装置
CN107831103A (zh) 一种压力脉冲衰减气测渗透率测试装置的精度评估方法
CN110927035A (zh) 一种低渗致密砂岩束缚水饱和度计算方法
CN107462936A (zh) 利用压力监测资料反演低渗透储层非达西渗流规律的方法
US20160334347A1 (en) Prediction of gas production rates from time-dependent nmr measurements
Lasseux et al. The “Step Decay”: a new transient method for the simultaneous determination of intrinsic permeability, Klinkenberg coefficient and porosity on very tight rocks
MX2007000778A (es) Metodo y aparato para calcular una distribucion de permeabilidad durante una prueba de pozo.
US5036193A (en) Earthen core analyzing means and method
Hinderdael et al. An analytical amplitude model for negative pressure waves in gaseous media
JP5286427B2 (ja) コンクリート構造物の空洞検査方法
CN115200977A (zh) 一种高温高压条件下岩心应力敏感评价装置及方法
Jannot et al. A simultaneous determination of permeability and Klinkenberg coefficient from an unsteady-state pulse-decay experiment
WO2019209121A1 (en) Method of testing an integrity of a structure separating a chamber from an adjacent environment, and related apparatus
Prince et al. Permeability estimation in tight gas sands and shales using NMR–a new interpretive methodology
CN211553699U (zh) 一种高温高压降内压孔隙度应力敏感性的测试装置
Abahri et al. Total pressure gradient incidence on hygrothermal transfer in highly porous building materials
Lasseux et al. Experimental measurement of the permeability of die-formed exfoliated graphite compression packings
Qin et al. Measurement of porosity for gas shale and estimation of its permeability
CN116341412A (zh) 一种基于方差分解的岩石渗透率模型参数敏感性评价方法

Legal Events

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

Effective date: 20210122