RU2523776C2 - Способ количественного расчета насыщенности трещинного коллектора углеводородами - Google Patents
Способ количественного расчета насыщенности трещинного коллектора углеводородами Download PDFInfo
- Publication number
- RU2523776C2 RU2523776C2 RU2011143258/28A RU2011143258A RU2523776C2 RU 2523776 C2 RU2523776 C2 RU 2523776C2 RU 2011143258/28 A RU2011143258/28 A RU 2011143258/28A RU 2011143258 A RU2011143258 A RU 2011143258A RU 2523776 C2 RU2523776 C2 RU 2523776C2
- Authority
- RU
- Russia
- Prior art keywords
- porosity
- fracture
- data
- results
- saturation
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 23
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 21
- 238000004364 calculation method Methods 0.000 title claims description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 36
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 238000004088 simulation Methods 0.000 claims abstract description 26
- 239000011148 porous material Substances 0.000 claims abstract description 23
- 238000005325 percolation Methods 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 230000035699 permeability Effects 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 4
- 238000005481 NMR spectroscopy Methods 0.000 claims description 3
- XOJVVFBFDXDTEG-UHFFFAOYSA-N Norphytane Natural products CC(C)CCCC(C)CCCC(C)CCCC(C)C XOJVVFBFDXDTEG-UHFFFAOYSA-N 0.000 claims description 3
- 238000001493 electron microscopy Methods 0.000 claims description 3
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 3
- 229910052753 mercury Inorganic materials 0.000 claims description 3
- 238000002459 porosimetry Methods 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 10
- 238000002474 experimental method Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 239000011435 rock Substances 0.000 description 7
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- -1 hydrocarbon hydrocarbon Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относится к геофизике. Сущность: способ включает определение пористости трещин и расчет показателя удельного сопротивления на различных глубинах трещинного коллектора на основе данных, полученных при помощи керна полного диаметра, и отображения данных каротажного зондирования; создание модели перколяционной сетки, сочетающей матрицу и трещину, при известных особенностях структуры пор; калибровку результатов численного моделирования в соответствии с моделью перколяционной сетки на основе данных эксперимента с использованием керна и анализа результатов, полученных при использовании герметизированого керна, с последующим установлением зависимости между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости; расчет насыщенности трещинного коллектора углеводородами посредством подбора интерполяционной функции. Технический результат: повышение точности. 9 з.п. ф-лы, 3 ил.
Description
Область техники
Настоящее изобретение относится к разработке нефтяных месторождений, в частности к методу количественного расчета насыщенности трещинного коллектора жидкими (газообразными) углеводородами на основе калибровки результатов численного моделирования с использованием данных анализа керна и герметизированого керна в сочетании с двойным боковым каротажным зондированием (БКЗ) и отображением данных каротажа в трещинном коллекторе, в частности в лаве или карбонатной породе.
Предшествующий уровень техники
Оставшиеся месторождения углеводородов в земной коре распределены главным образом между четырьмя источниками, включая залежи литологического типа, прибрежные залежи, залежи среднего и глубинного уровня наложенных и внутренних впадин, из которых залежи литологического типа составляют более 60% от вновь разведанных запасов и таким образом представляет собой основной источник замещения законсервированных месторождений нефти и увеличения ее запасов. Значительная часть обнаруженных за последние годы крупных нефтяных карбонатных и вулканических коллекторов и месторождений газа представляет собой коллекторы литологического типа. Указанные коллекторы часто имеют комплексную пористость и содержат разнообразные поры, в частности каверны и трещины. В связи с этим количественный расчет насыщенности трещинного коллектора углеводородами имеет большое значение для повышения эффективности разработки литологического коллектора.
До настоящего времени электрическое зондирование, примером которого является двойное боковое каротажное зондирование, является основным способом оценки насыщенности углеводородами. Оценку насыщенности углеводородами с использованием электрического зондирования обычно проводят на основе формулы Арчи и расширенных формул. Расширенные формулы для расчета насыщенности углеводородами трещинного коллектора выведены из модели простой трещины. Исследования по расчету насыщенности трещинного коллектора описаны в монографии [Yandong Tan. Interpretation Model and Evaluation Method of Well Logging in Fractured Reservoir (1987, Petroleum Industrial Press)]. На основе модели простой трещины автором получены аналитические выражения для удельного сопротивления и показателя удельного сопротивления для породы с горизонтальными, вертикальными и сетчатыми трещинами. В монографии [Liangxiao Zhao. Evaluation Technique of Well Logging in Carbonate Rock Reservoir (1994, Petroleum Industrial Press)] предложен способ расчета насыщенности коллектора трещинами и порами. Все способы количественного расчета насыщенности трещинного коллектора, разработанные в более поздних исследованиях, основаны на описанных выше классических способах.
В классических способах расчета насыщенности трещинного коллектора существует три проблемы. Во-первых, модель простой трещины не отражает закономерностей реального распределения трещин в коллекторе. Во-вторых, способ на основе простого последовательного и/или параллельного расчета не полностью учитывает электрические свойства породы, находящейся под воздействием трещин. Влияние трещин на электрические свойства породы включает как перенос электричества самими трещинами, так и изменение электрических свойств, вызванное влиянием насыщенности матрицы. В-третьих, эти модели часто включают такие параметры, как показатель трещинной пористости (mf), показатель насыщенности трещины (nf), насыщенность трещины связанной водой (Sfb) и т.д., точное определение которых крайне затруднительно. На практике наблюдается большое расхождение между результатами расчета согласно классической модели насыщенности в трещинных коллекторах и результатами анализа герметизированого керна.
Кроме того, чисто численное моделирование электрических свойств трещинного коллектора учитывает только влияние их относительного изменения, которое сильно отличается от характеристик реального коллектора и, соответственно, не может использоваться на практике для оценки коллектора посредством зондирования. В связи с этим оценка насыщенности трещинного коллектора углеводородами представляет собой одну из подлежащих решению проблем.
Сущность изобретения
Задачей настоящего изобретения является разработка способа количественного расчета насыщенности трещинного коллектора углеводородами с целью повышения точности оценки каротажной диаграммы трещинного коллектора.
Для достижения вышеуказанной задачи настоящее изобретение включает следующие стадии С1-С7.
С1. На основе известных данных керна полного диаметра определяют зависимость удельного сопротивления (I) от водонасыщенности (Sw) матрицы;
на основе анализа данных герметизированого керна определяют реальную насыщенность водой (нефтью/газом) нетронутой (исходной) формации трещинного коллектора.
С2. Определяют трещинную пористость на различных глубинах трещинного коллектора с использованием графического отображения известных данных каротажного зондирования. Показатели удельного сопротивления на различных глубинах рассчитывают с использованием данных двойного БКЗ в сочетании с данными анализа керна.
Показатели удельного сопротивления на различных глубинах рассчитывают на стадии С2 следующим образом: удельное сопротивление нетронутой формации Rt определяют по данным глубинного БКЗ, а удельное сопротивление водонасыщенной формации R0 рассчитывают по формуле Арчи на основе диаграммы каротажа пористости в сочетании с данными анализа керна.
Отображение данных каротажного зондирования на стадии С2 представляет собой отображение данных полнопроходного сканирования микросопротивления или электроакустического каротажа.
С3. Создают модель перколяционной сетки, включающей матрицу и трещину, на основе известных данных о региональных особенностях структуры пор, пористости и проницаемости, и посредством дренажного моделирования определяют взаимосвязь между показателем удельного сопротивления (I) и водонасыщенностью (Sw) трещинного коллектора.
Модель перколяционной сетки, включающей матрицу и трещину, на стадии С3 создают посредством определения радиусов пор и каналов и их распределений, определения соединений между порами и установления трехмерной модели пористой структуры.
Региональные особенности структуры пор, описанные на стадии С3, могут быть определены посредством анализа методами ртутной порометрии, ядерного магнитного резонанса и электронной микроскопии.
При создании модели на стадии С3 параметры модели, включая радиусы и их распределение, пористость и проницаемость матрицы, соответствуют данным анализа керна.
С4. На основе данных экспериментов с использованием керна и герметизированого керна производят калибровку результатов численного моделирования в соответствии с моделью перколяционной сетки.
Результаты численного моделирования на стадии С4 калибруют следующим образом:
С4-1. Результаты численного моделирования при трещинной пористости, равной 0, на стадии С3 калибруют на основе результатов экспериментов по определению электрических свойств керна (зависимость I-Sw) в заданном интервале. По своей сути способ заключается в аппроксимации параметров при численном моделировании таким образом, чтобы результаты численного моделирования при трещинной пористости, равной 0, совпали с результатами экспериментов по измерению электрических свойств керна.
С4-2. Затем результаты моделирования, полученные на стадии С3, калибруют с использованием данных анализа водонасыщенности, полученных с использованием герметизированого керна, показателя удельного сопротивления на соответствующей глубине, рассчитанного на стадии С2, и трещинной пористости на соответствующей глубине, рассчитанной на стадии С1. По своей сути способ заключается в аппроксимации параметров при численном моделировании таким образом, чтобы результаты численного моделирования при трещинной пористости на соответствующей глубине совпали с результатами анализа герметизированого керна.
Рассмотренные выше параметры моделирования включают главным образом структуру пор, свойства флюида и т.д.
С5. После калибровки получают зависимость между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости посредством численного моделирования; трещинная пористость, описанная на стадии С5, находится в пределах 0,1%-0,5%.
С6. Посредством интерполяции получают зависимость удельного сопротивления от водонасыщенности при различной трещинной пористости; оптимальная интерполяционная функция, полученная на стадии С6, имеет вид:
где I - показатель удельного сопротивления; Sw - водонасыщенность; а, и и с - константы.
При матричной пористости 3% и трещинной пористости 0,0%: a=0,03, b=3,44, c=0,52;
при матричной пористости 3% и трещинной пористости 0,1%: a=0,06, b=2,70, c=0,52;
при матричной пористости 3% и трещинной пористости 0,2%: a=0,20, b=1,47, c=0,67;
при матричной пористости 3% и трещинной пористости 0,3%: a=0,38, b=0,85, c=0,77;
при матричной пористости 3% и трещинной пористости 0,4%: a=0,52, b=0,53, c=0,91.
С7. Насыщенность трещинного коллектора углеводородами рассчитывают в соответствии с трещинной пористостью, определенной на стадии С2, и с использованием интерполяционной функции, выбранной на стадии С6.
Интерполяционную функцию, используемую на стадии С7, выбирают с целью определения соответствующей зависимости показателя удельного сопротивления от водонасыщенности. Водонасыщенность и насыщенность углеводородами на определенной глубине рассчитывают на основе показателя удельного сопротивления, рассчитанного на стадии С2.
Краткое описание графических материалов
На фиг.1 представлена модель пористой структуры, включающей матрицу и трещину; средний слой соответствует трещине, а верхний и нижний слои соответствуют матрице. В матрице все пространство пор изображено в виде шариков, а каналы - в виде линий, при этом размер шарика соответствует объему поры, а толщина линии соответствует размеру канала.
На фиг.2 представлены данные для матричной пористости 3%, а трещинная пористость (FV РА) имеет значения от 0% до 0,4%. На оси абсцисс изображена водонасыщенность, а на оси ординат - показатель удельного сопротивления. Крайняя правая кривая на фиг.2 соответствует экспериментальной зависимости показателя удельного сопротивления от водонасыщенности, полученной из анализа керна при трещинной пористости, равной 0. Крайняя левая кривая на
фиг.2 соответствует результатам моделирования на основе калибровки данных анализа герметизированого керна. Остальные кривые на фиг.2 соответствуют зависимостям показателя удельного сопротивления от водонасыщенности, полученным посредством численного моделирования с использованием в качестве ограничений данных анализа керна и анализа герметизированого керна при трещинной пористости соответственно 0,1%, 0,2% и 0,3%.
На фиг.3 приведен пример результатов обработки результатов каротажа с использованием предлагаемого в настоящем изобретении способа количественного расчета насыщенности трещинного коллектора водой (углеводородами).
Подробное описание изобретения
Ниже приводится подробное описание конкретного осуществления настоящего изобретения со ссылкой на прилагаемые фигуры.
1. Выполняют эксперимент с керном полного диаметра и анализ герметизированного керна трещинного коллектора: на основе данных эксперимента с керном полного диаметра определяют зависимость показателя удельного сопротивления (I) от насыщенности матрицы водой (Sw); реальную насыщенность водой (углеводородами) нетронутой формации трещинного коллектора определяют посредством анализа герметизированного керна.
2. Трещинную пористость на различных глубинах определяют посредством обработки отображения данных зондирования (FMI, STAR и т.д.). Показатели удельного сопротивления на различных глубинах рассчитывают с использованием данных двойного БКЗ в сочетании с данными анализа керна. Удельное сопротивление нетронутой формации Rt определяют по данным глубинного БКЗ, а удельное сопротивление водонасыщенной формации RO рассчитывают по формуле Арчи на основе пористости по данным каротажа в сочетании с данными анализа керна.
3. Модель перколяционной сетки, включающей матрицу и трещину, в исследуемой области устанавливают на основе данных анализа структурных особенностей пор (методами ртутной порометрии, ядерного магнитного резонанса, электронной микроскопии и т.д.), анализируемой пористости и данных о проницаемости. Зависимость между показателем удельного сопротивления (I) и насыщенностью трещинного коллектора водой (Sw) изучают посредством дренажного моделирования. При создании численной модели ее параметры, такие как радиусы и распределения пор и каналов, пористость, проницаемость и т.д., должны с достаточно высокой точностью соответствовать результатам анализа керна. На фиг.1 представлена полученная модель пористой структуры, включающей матрицу и трещину; средний слой соответствует трещине, а верхний и нижний слои соответствуют матрице. В матрице все пространство пор изображено в виде шариков, а каналы - в виде линий, при этом размер шарика соответствует объему поры, а толщина линии соответствует размеру канала.
4. Результаты численного моделирования калибруют на основе данных эксперимента с использованием керна и анализа герметизированого керна. Результаты моделирования при трещинной пористости, равной 0, калибруют в соответствии с результатами экспериментов по определению электрических свойств керна (зависимость I-Sw) в заданном интервале. Последовательность осуществления способа включает: приведение результатов численного моделирования при трещинной пористости, равной нулю, в соответствие с результатами эксперимента по определению электрических свойств керна (крайняя правая кривая на фиг.2) за счет подгонки параметров модели, полученной посредством численного моделирования (к указанным параметрам относятся структура пор, свойства флюида и т.д.).
Результаты моделирования калибруют с использованием водонасыщенности по данным анализа герметизированого керна, показателя удельного сопротивления на соответствующей глубине, рассчитанного на рассмотренной выше стадии 2, и трещинной пористости на соответствующей глубине, рассчитанной на рассмотренной выше стадии 1. Подробная последовательность расчета в соответствии с данным способом включает:
приведение результатов численного моделирования при трещинной пористости, соответствующей заданной глубине (указанное значение определено по данным отображения результатов каротажа), в соответствии с результатами анализа герметизированного керна (крайняя левая кривая на фиг.2) за счет подгонки параметров модели, полученной посредством численного моделирования (к указанным параметрам относятся структура пор, свойства флюида и т.д.).
5. После калибровки с использованием данных анализа керна и герметизированного керна моделирование осуществляют повторно с целью установления зависимости между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости. На фиг.2 представлена зависимость между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при трещинной пористости 0,1%, 0,2% и 0,3%.
6. Устанавливают интерполяционную функцию между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости. На основе результатов моделирования, представленных на фиг.2, по данным интерполяционного анализа определены следующие оптимальные зависимости между показателем удельного сопротивления (I) и водонасыщенностью (Sw):
при трещинной пористости, равной 0,0%:
при трещинной пористости, равной 0.1%:
при трещинной пористости, равной 0,2%:
при трещинной пористости, равной 0,3%:
при трещинной пористости, равной 0,4%:
7. Рассчитывают насыщенность трещинного коллектора углеводородами. Количественный расчет насыщенности трещинного коллектора водой (углеводородами) производят в соответствии с трещинной пористостью, определенной по данным отображения результатов каротажа, выбранной функцией зависимости между установленными ранее показателем удельного сопротивления (I) и водонасыщенностью (Sw), и данными глубинного БКЗ. Например, в соответствии с результатами анализа отображения данных каротажа при выборе для расчета насыщенности углеводородами приведенной выше формулы (2) трещинная пористость на определенной глубине составляет 0,1%. На фиг.3 приведена практическая обработка результатов в определенном интервале для определенной скважины в нефтяном месторождении Daqing в Китае, где профиль, ближайший к крайнему правому, соответствует результатам анализа насыщенности газом в соответствии с традиционной моделью трещины, а крайний правый профиль соответствует результату расчета насыщенности газом в соответствии со способом, предлагаемым в настоящем изобретении. Сопоставление представленных данных показывает, что результат расчета насыщенности газом в соответствии со способом, представленным в настоящем изобретении, хорошо согласуется с результатом анализа данных, полученных с использованием герметизированного керна.
Промышленное применение
В традиционных способах расчета насыщенности трещинного коллектора существует три проблемы. Во-первых, модель простой трещины не отражает закономерностей реального распределения трещин в коллекторе. Во-вторых, способ на основе простого последовательного и/или параллельного расчета не полностью учитывает электрические свойства породы, находящейся под воздействием трещин. Влияние трещин на электрические свойства породы включает как перенос электричества самими трещинами, так и изменение электрических свойств, вызванное влиянием насыщенности матрицы. В-третьих, эти модели часто включают такие параметры, как показатель трещинной пористости (mf), показатель насыщенности трещины (nf), насыщенность трещины связанной водой (Sfb) и т.д., точное определение которых крайне затруднительно. На практике наблюдается большое расхождение между результатами расчета согласно классической модели насыщенности в трещинных коллекторах и результатами анализа с использованием герметизированного керна.
Кроме того, чисто численное моделирование электрических свойств трещинного коллектора учитывает только влияние их относительного изменения, которое сильно отличается от характеристик реального коллектора и, соответственно, не может использоваться на практике для оценки коллектора путем каротажного зондирования. В связи с этим оценка насыщенности трещинного коллектора углеводородами представляет собой одну из подлежащих решению проблем.
В сравнении с традиционными моделями для расчета насыщенности трещинного коллектора настоящее изобретение способно значительно повысить точность расчета насыщенности трещинного коллектора углеводородами. При матричной пористости 3%, трещинной пористости 0,2% и показателе удельного сопротивления 4 насыщенность углеводородами в соответствии с традиционным способом расчета составляет 0,49, в то время как в соответствии со способом, предлагаемым в настоящем изобретении, она составляет 0,67, что превышает результат, полученный при использовании традиционного способа расчета, на 0,18.
Claims (10)
1. Способ количественного расчета насыщенности трещинного коллектора углеводородами, отличающийся тем, что включает следующие стадии:
С1) на основе данных керна полного диаметра определяют зависимость удельного сопротивления (I) от водонасыщенности (Sw) матрицы; на основе анализа данных герметизированного керна определяют реальную насыщенность водой (нефтью/газом) исходной формации трещинного коллектора;
С2) определяют трещинную пористость на различных глубинах трещинного коллектора с использованием графического отображения известных данных каротажного зондирования; рассчитывают показатели удельного сопротивления на различных глубинах с использованием данных двойного бокового каротажного зондирования в сочетании с данными анализа керна;
С3) создают модель перколяционной сетки, включающей матрицу и трещину, в соответствии с известными региональными особенностями структуры пор, пористостью и проницаемостью исследованной области, и определяют зависимость между показателем удельного сопротивления (I) и водонасыщенностью (Sw) трещинного коллектора посредством дренажного моделирования;
С4) осуществляют калибровку результатов численного моделирования в соответствии с моделью перколяционной сетки на основе данных анализа керна и результатов анализа герметизированого керна;
С5) после калибровки определяют вариации зависимости между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости посредством численного моделирования;
С6) определяют функцию зависимости между показателем удельного сопротивления и водонасыщенностью при различной трещинной пористости посредством интерполяции;
С7) рассчитывают насыщенность трещинного коллектора углеводородами в соответствии с пористостью трещин, определенной на стадии С2, посредством выбора интерполяционной модели на основе функций зависимости между показателем удельного сопротивления и водонасыщенностью, определенных на стадии С6.
С1) на основе данных керна полного диаметра определяют зависимость удельного сопротивления (I) от водонасыщенности (Sw) матрицы; на основе анализа данных герметизированного керна определяют реальную насыщенность водой (нефтью/газом) исходной формации трещинного коллектора;
С2) определяют трещинную пористость на различных глубинах трещинного коллектора с использованием графического отображения известных данных каротажного зондирования; рассчитывают показатели удельного сопротивления на различных глубинах с использованием данных двойного бокового каротажного зондирования в сочетании с данными анализа керна;
С3) создают модель перколяционной сетки, включающей матрицу и трещину, в соответствии с известными региональными особенностями структуры пор, пористостью и проницаемостью исследованной области, и определяют зависимость между показателем удельного сопротивления (I) и водонасыщенностью (Sw) трещинного коллектора посредством дренажного моделирования;
С4) осуществляют калибровку результатов численного моделирования в соответствии с моделью перколяционной сетки на основе данных анализа керна и результатов анализа герметизированого керна;
С5) после калибровки определяют вариации зависимости между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости посредством численного моделирования;
С6) определяют функцию зависимости между показателем удельного сопротивления и водонасыщенностью при различной трещинной пористости посредством интерполяции;
С7) рассчитывают насыщенность трещинного коллектора углеводородами в соответствии с пористостью трещин, определенной на стадии С2, посредством выбора интерполяционной модели на основе функций зависимости между показателем удельного сопротивления и водонасыщенностью, определенных на стадии С6.
2. Способ по п.1, отличающийся тем, что показатели удельного сопротивления на различных глубинах рассчитывают следующим образом: удельное сопротивление нетронутой формации Rt определяют по данным глубинного бокового каротажного зондирования, а удельное сопротивление водонасыщенной формации R0 рассчитывают по формуле Арчи на основе данных каротажного зондирования пористости в сочетании с данными анализа керна.
3. Способ по п.1, отличающийся тем, что в качестве используемых на стадии С2 данных каротажного зондирования используют отображения данных полнопроходного сканирования микросопротивления или электроакустического каротажа.
4. Способ по п.1, отличающийся тем, что модель перколяционной сетки, включающей матрицу и трещину, на стадии С3 создают посредством определения радиусов и распределения пор и определения радиусов и распределения каналов, определения соединений между порами и создания трехмерной модели структуры пор.
5. Способ по п.1, отличающийся тем, что региональные особенности структуры пор на стадии С3 могут быть определены методами ртутной порометрии, ядерного магнитного резонанса и электронной микроскопии.
6. Способ по п.1, отличающийся тем, что при создании модели на стадии С3 используемые в модели параметры, главным образом радиусы и распределение пор в матрице, пористость, проницаемость и аналогичные параметры, соответствуют результатам анализа керна.
7. Способ по п.1, отличающийся тем, что результаты численного моделирования на стадии С4 калибруют следующим образом:
результаты численного моделирования при трещинной пористости, равной 0, полученные на стадии С3, калибруют на основе результатов измерения электрических свойств керна (зависимость I-Sw) в заданном интервале, в частности, результаты численного моделирования при трещинной пористости, равной 0, приводят в соответствие с результатами измерения электрических свойств керна посредством подгонки параметров численного моделирования;
затем результаты моделирования, полученные на стадии С3, калибруют с использованием значений водонасыщенности, определенных по результатам анализа герметизированного керна, показателя удельного сопротивления на соответствующей глубине, рассчитанного на стадии С2, и трещинной пористости на соответствующей глубине, рассчитанной на стадии С1; в частности, результаты численного моделирования при трещинной пористости, соответствующей указанной глубине, приводят в соответствие с результатами анализа герметизированого керна посредством подгонки параметров численного моделирования.
результаты численного моделирования при трещинной пористости, равной 0, полученные на стадии С3, калибруют на основе результатов измерения электрических свойств керна (зависимость I-Sw) в заданном интервале, в частности, результаты численного моделирования при трещинной пористости, равной 0, приводят в соответствие с результатами измерения электрических свойств керна посредством подгонки параметров численного моделирования;
затем результаты моделирования, полученные на стадии С3, калибруют с использованием значений водонасыщенности, определенных по результатам анализа герметизированного керна, показателя удельного сопротивления на соответствующей глубине, рассчитанного на стадии С2, и трещинной пористости на соответствующей глубине, рассчитанной на стадии С1; в частности, результаты численного моделирования при трещинной пористости, соответствующей указанной глубине, приводят в соответствие с результатами анализа герметизированого керна посредством подгонки параметров численного моделирования.
8. Способ по любому из пп.1-7, отличающийся тем, что параметры численного моделирования включают параметры структуры пор и параметры свойств флюида.
9. Способ по п.1, отличающийся тем, что значения пористости на стадии С5 находятся в пределах от 0,1% до 0,5%.
10. Способ по п.1, отличающийся тем, что оптимальная интерполяционная функция, определяемая на стадии С6, имеет вид:
где I - показатель удельного сопротивления; Sw - водонасыщенность; a, b и c - константы,
при матричной пористости 3% и трещинной пористости 0,0%: a=0,03, b=3,44, c=0,52;
при матричной пористости 3% и трещинной пористости 0,1%: a=0,06, b=2,70, c=0,52;
при матричной пористости 3% и трещинной пористости 0,2%: a=0,20, b=1,47, c=0,67;
при матричной пористости 3% и трещинной пористости 0,3%: a=0,38, b=0,85, c=0,77;
при матричной пористости 3% и трещинной пористости 0,4%: a=0,52, b=0,53, c=0,91.
где I - показатель удельного сопротивления; Sw - водонасыщенность; a, b и c - константы,
при матричной пористости 3% и трещинной пористости 0,0%: a=0,03, b=3,44, c=0,52;
при матричной пористости 3% и трещинной пористости 0,1%: a=0,06, b=2,70, c=0,52;
при матричной пористости 3% и трещинной пористости 0,2%: a=0,20, b=1,47, c=0,67;
при матричной пористости 3% и трещинной пористости 0,3%: a=0,38, b=0,85, c=0,77;
при матричной пористости 3% и трещинной пористости 0,4%: a=0,52, b=0,53, c=0,91.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200910087474.3 | 2009-06-22 | ||
| CN200910087474A CN101929973B (zh) | 2009-06-22 | 2009-06-22 | 裂缝储层含油气饱和度定量计算方法 |
| PCT/CN2010/000420 WO2010148628A1 (zh) | 2009-06-22 | 2010-03-31 | 裂缝储层含油气饱和度定量计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2011143258A RU2011143258A (ru) | 2013-07-27 |
| RU2523776C2 true RU2523776C2 (ru) | 2014-07-20 |
Family
ID=43369269
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2011143258/28A RU2523776C2 (ru) | 2009-06-22 | 2010-03-31 | Способ количественного расчета насыщенности трещинного коллектора углеводородами |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US8805661B2 (ru) |
| EP (1) | EP2447469B1 (ru) |
| CN (1) | CN101929973B (ru) |
| AU (1) | AU2010265746B2 (ru) |
| RU (1) | RU2523776C2 (ru) |
| WO (1) | WO2010148628A1 (ru) |
Families Citing this family (102)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102175832B (zh) * | 2011-01-10 | 2013-10-16 | 中国石油天然气股份有限公司 | 一种确定典型储层最佳饱和度计算模型的方法 |
| WO2012114146A2 (en) * | 2011-02-23 | 2012-08-30 | Total Sa | Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region |
| CN102262041B (zh) * | 2011-04-20 | 2013-03-13 | 中国石油天然气股份有限公司 | 一种基于多谱孔隙结构分析的饱和度确定方法 |
| CN102243196B (zh) * | 2011-04-26 | 2012-07-04 | 中国石油大学(华东) | 一种实验室快速测定饱和度指数n的方法 |
| EP2769243A4 (en) * | 2011-10-18 | 2017-03-15 | Saudi Arabian Oil Company | 4d saturation modeling |
| CN103064114B (zh) * | 2011-10-18 | 2016-05-04 | 中国石油化工股份有限公司 | 一种裂缝储层的表征方法及装置 |
| CN102540265B (zh) * | 2011-12-21 | 2014-06-25 | 西南石油大学 | 一种基于网络模拟的多孔介质含水饱和度计算方法 |
| CN103376469B (zh) * | 2012-04-26 | 2017-09-26 | 中国石油集团长城钻探工程有限公司 | 一种基于超声成像测井的裂缝定量评价方法 |
| US9835027B2 (en) | 2012-06-14 | 2017-12-05 | Halliburton Energy Services, Inc. | System, method, and computer program product to determine placement of fracture stimulation points using minerology |
| CN103543474A (zh) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | 一种基于非导电孔隙模型的含油气饱和度评价方法 |
| WO2014154242A1 (en) * | 2013-03-25 | 2014-10-02 | Statoil Petroleum As | Determining a condition of a subsurface reservoir |
| CN103499679B (zh) * | 2013-09-11 | 2015-06-03 | 中国石油天然气股份有限公司 | 一种测定孔洞型储层的饱和度的方法及设备 |
| CN103616731B (zh) * | 2013-11-19 | 2016-04-06 | 中国石油天然气股份有限公司 | 一种油气勘探中蚀变火山岩有效储层确定方法及装置 |
| ITMI20131986A1 (it) * | 2013-11-28 | 2015-05-29 | Eni Spa | Metodo per la misurazione della saturazione in gas intrappolato in un campione di roccia |
| CN103675945B (zh) * | 2013-12-17 | 2017-03-08 | 中国石油天然气股份有限公司 | 一种测定孔洞型储层的饱和度的方法及设备 |
| CN103698494B (zh) * | 2013-12-30 | 2015-04-29 | 中国石油大学(北京) | 一种确定岩性圈闭含油气饱和度的方法及装置 |
| CN103728671B (zh) * | 2014-01-15 | 2016-11-16 | 中国石油大学(北京) | 确定碳酸盐岩潜山含油气性的方法和装置 |
| CN104358565B (zh) * | 2014-10-31 | 2019-03-22 | 中国石油化工股份有限公司 | 一种气藏含气饱和度的确定方法 |
| CN104809283B (zh) * | 2015-04-22 | 2018-04-06 | 中国石油天然气股份有限公司 | 一种稠油油藏内的夹层识别方法 |
| CN104948176B (zh) * | 2015-05-08 | 2017-10-24 | 西南石油大学 | 一种基于渗透增大率识别碳酸盐岩储层裂缝的方法 |
| CN104865610B (zh) * | 2015-05-26 | 2017-06-13 | 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 | 一种深部储层岩心的数值表征方法 |
| CN106285656A (zh) * | 2015-06-09 | 2017-01-04 | 中国石油化工股份有限公司 | 用于确定页岩气储层的游离气饱和度的方法及系统 |
| FR3038408B1 (fr) * | 2015-06-30 | 2017-08-25 | Services Petroliers Schlumrberger | Modelisation de la saturation et permeabilite de reservoir de champ petrolifere |
| CN105021506A (zh) * | 2015-07-09 | 2015-11-04 | 中国石油大学(华东) | 一种基于孔隙网络模型的三相相对渗透率的计算方法 |
| CN105158803B (zh) * | 2015-07-27 | 2017-10-17 | 中国石油天然气集团公司 | 一种选取优势储层的方法 |
| CN105093340A (zh) * | 2015-08-05 | 2015-11-25 | 中国石油天然气股份有限公司 | 煤层可动水体积的获取方法及获取装置 |
| CN105240001B (zh) * | 2015-09-23 | 2018-03-09 | 中国石油大学(北京) | 核磁共振测井孔隙度校正方法及装置 |
| CN105608628B (zh) * | 2015-12-21 | 2020-12-18 | 中国石油大学(华东) | 储层中深部异常高孔隙带类型划分及定量厘定方法 |
| CN105528732B (zh) * | 2015-12-22 | 2020-01-10 | 中国石油天然气股份有限公司 | 一种试气井产能预测方法 |
| CN105629316B (zh) * | 2016-01-13 | 2018-10-12 | 中国石油大学(北京) | 获得原状地层的流体径向油水饱和度变化的方法及装置 |
| CN105844011B (zh) * | 2016-03-22 | 2018-10-16 | 中国石油天然气集团公司 | 一种基于毛管模型的渗透率计算方法 |
| CN105930932B (zh) * | 2016-04-25 | 2019-11-12 | 中石化石油工程技术服务有限公司 | 基于含气指数的页岩气层标准化无阻流量的获取方法 |
| CN107355219B (zh) * | 2016-05-09 | 2020-09-04 | 中国石油化工股份有限公司 | 裂缝性地层模型和其使用方法 |
| CN106093083B (zh) * | 2016-07-01 | 2018-08-07 | 清能艾科（深圳）能源技术有限公司 | 一种采用数字岩心模拟建立岩电关系的方法 |
| CN106204302B (zh) * | 2016-07-05 | 2019-09-10 | 中国石油天然气股份有限公司 | 一种计算原始含水饱和度的方法与应用 |
| CN106443802B (zh) * | 2016-09-20 | 2018-11-13 | 中石化石油工程技术服务有限公司 | 一种基于形态学滤波的电成像储层缝洞体的定量表征方法及系统 |
| CN107885893B (zh) * | 2016-09-29 | 2021-06-01 | 中国石油化工股份有限公司 | 描述二氧化碳驱储层非均质性的构建方法 |
| CN106437695B (zh) * | 2016-10-17 | 2023-07-21 | 中国石油大学(华东) | 一种适用于稠油渗流规律的新型电模拟装置 |
| CN106370582B (zh) * | 2016-10-28 | 2023-03-17 | 陕西延长石油(集团)有限责任公司研究院 | 一种模拟裂缝性特低渗油藏动态渗吸的实验装置及其应用 |
| CN107143331A (zh) * | 2017-06-26 | 2017-09-08 | 中国石油大学(华东) | 一种非均质储集层含油饱和度测量实验装置 |
| CN107807396B (zh) * | 2017-09-15 | 2019-07-09 | 中国石油天然气股份有限公司 | 地层基质电阻率的确定方法和装置 |
| CN107808032B (zh) * | 2017-09-28 | 2021-07-16 | 中国石油化工股份有限公司 | 页岩气储层含水饱和度计算方法及系统 |
| CA3089719A1 (en) * | 2018-02-06 | 2019-08-15 | Conocophillips Company | 4d seismic as a method for characterizing fracture network and fluid distribution in unconventional reservoir |
| CN108562944A (zh) * | 2018-03-12 | 2018-09-21 | 中国石油天然气集团有限公司 | 一种基于电成像孔隙度谱的孔隙结构参数计算方法 |
| CN110716024A (zh) * | 2018-07-11 | 2020-01-21 | 中国石油天然气股份有限公司 | 物性曲线的确定方法及装置 |
| CN110807235B (zh) * | 2018-08-01 | 2023-07-25 | 中国石油天然气股份有限公司 | 一种模拟k气藏断缝孔渗流的方法 |
| CN109376481B (zh) * | 2018-08-16 | 2022-12-09 | 清能艾科（深圳）能源技术有限公司 | 基于多gpu的数字岩心相渗曲线的计算方法、装置及计算机设备 |
| CN109281662B (zh) * | 2018-08-16 | 2022-03-29 | 中国石油天然气股份有限公司 | 水平井不同节段位置储层相对可压性指数的确定方法 |
| CN109903387B (zh) * | 2018-08-16 | 2022-12-09 | 清能艾科（深圳）能源技术有限公司 | 孔隙尺度的聚合物驱替模型建立方法及装置、电子设备 |
| KR102203949B1 (ko) * | 2018-08-30 | 2021-01-18 | 한양대학교 산학협력단 | 오일 생산량 예측 방법 및 그 장치 |
| CN110939428B (zh) * | 2018-09-21 | 2023-08-01 | 中国石油化工股份有限公司 | 一种致密砂岩油气藏储层裂缝的识别方法 |
| CN111563306B (zh) * | 2019-01-25 | 2022-11-04 | 中国石油天然气股份有限公司 | 用于岩石孔隙结构仿真的逾渗网络模型构建方法及设备 |
| CN111485874B (zh) * | 2019-01-28 | 2023-09-26 | 中国石油天然气股份有限公司 | 确定水平井储层含水饱和度的方法及装置 |
| CN111521532B (zh) * | 2019-02-01 | 2023-02-10 | 中国石油天然气股份有限公司 | 一种致密储层含气饱和度的测定方法 |
| CN109944587B (zh) * | 2019-03-12 | 2020-04-14 | 东北石油大学 | 定量表征变形带对流体注采效率影响的方法 |
| CN111693427B (zh) * | 2019-03-15 | 2023-02-10 | 中国石油天然气股份有限公司 | 油气藏流体可动性的分析方法 |
| CN110806435A (zh) * | 2019-03-26 | 2020-02-18 | 中国石油大学(北京) | 一种co2-水-岩相互作用对致密岩石孔喉结构影响的交流阻抗谱评价方法 |
| CN110056347A (zh) * | 2019-04-28 | 2019-07-26 | 中国石油天然气集团有限公司 | 一种基于包络幅度的油气含量计算方法 |
| CN110163533B (zh) * | 2019-06-03 | 2020-03-03 | 西南石油大学 | 一种页岩气缝网综合可压性评价方法 |
| CN112145164B (zh) * | 2019-06-28 | 2024-02-23 | 中国石油化工股份有限公司 | 确定地层水体积压缩系数和天然气体积压缩系数的方法 |
| CN110334868B (zh) * | 2019-07-08 | 2020-12-08 | 西南石油大学 | 一种耦合流体流动与地质应力预测最优焖井时间的方法 |
| CN112922587B (zh) * | 2019-11-20 | 2023-10-27 | 中国石油天然气股份有限公司 | 盐霜位置的识别方法及装置 |
| CN110853138B (zh) * | 2019-11-21 | 2023-08-18 | 科吉思石油技术咨询(北京)有限公司 | 双重介质碳酸盐岩孔隙-裂缝双重网络模型构建方法 |
| CN111008482B (zh) * | 2019-12-18 | 2023-04-07 | 长江大学 | 一种变质岩储层含水饱和度计算方法及装置 |
| CN110955982B (zh) * | 2019-12-18 | 2023-04-07 | 长江大学 | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 |
| CN112145165B (zh) * | 2019-12-26 | 2024-04-23 | 中国海洋石油集团有限公司 | 一种微裂缝-孔隙型储层动静态渗透率转换方法 |
| CN113049321B (zh) * | 2019-12-26 | 2023-01-10 | 中国石油天然气股份有限公司 | 储层模型的制备方法及储层模型 |
| CN111353218B (zh) * | 2020-02-20 | 2023-03-24 | 西安石油大学 | 一种煤层气-致密气储层合压性的测井定量评价方法 |
| CN111274528B (zh) * | 2020-03-02 | 2021-09-17 | 中国石油大学(北京) | 储层裂缝渗吸质量预测方法及系统 |
| CN111398117B (zh) * | 2020-03-15 | 2023-05-26 | 长江大学 | 一种岩心电频散分析方法及装置 |
| CN111322065A (zh) * | 2020-04-09 | 2020-06-23 | 中国石油天然气集团有限公司 | 一种基于等效电阻率迭代的水淹层含水饱和度计算方法 |
| CN113625360B (zh) * | 2020-05-08 | 2024-02-23 | 中国石油化工股份有限公司 | 微裂缝地层产量预测方法、系统、电子设备及介质 |
| CN113685174B (zh) * | 2020-05-19 | 2023-09-26 | 中国石油天然气股份有限公司 | 相态变化对致密油井产能影响计算方法及装置 |
| CN111622750B (zh) * | 2020-05-25 | 2023-07-25 | 中国石油天然气集团有限公司 | 一种常规测井资料评价砂砾岩储层有效性的计算方法 |
| CN111814347B (zh) * | 2020-07-20 | 2024-06-11 | 中国石油大学(华东) | 预测油藏内气窜通道的方法及系统 |
| CN112098293B (zh) * | 2020-08-03 | 2021-06-01 | 西南石油大学 | 基于孔隙裂缝双重介质气藏非稳态气水两相渗流模拟方法 |
| CN111894568A (zh) * | 2020-08-04 | 2020-11-06 | 中国地质大学（北京） | 裂缝性碳酸盐岩油藏饱和度模型数字岩心分析方法 |
| CN114086945B (zh) * | 2020-08-06 | 2024-04-30 | 中国石油天然气集团有限公司 | 一种低矿化度成因低阻油层的录井识别方法 |
| CN111963162B (zh) * | 2020-09-09 | 2023-07-25 | 中国石油天然气股份有限公司 | 基于岩性与水性的双参数流体性质识别方法 |
| CN114185108B (zh) * | 2020-09-14 | 2024-04-12 | 中国石油化工股份有限公司 | 裂缝孔隙度确定方法、存储介质和计算机设备 |
| CN112379416B (zh) * | 2020-10-13 | 2024-02-06 | 北京恒泰兴科信息技术有限公司 | 煤岩岩石物理建模预测横波方法、装置及电子设备 |
| CN112431588B (zh) * | 2020-10-15 | 2023-07-25 | 中国石油天然气股份有限公司 | 一种侧钻选井优选方法 |
| CN112392476B (zh) * | 2020-12-02 | 2022-02-15 | 西南石油大学 | 利用常规测井资料确定低渗裂缝孔渗参数的方法 |
| CN112784404B (zh) * | 2021-01-05 | 2023-12-22 | 中国石油天然气集团有限公司 | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 |
| CN113655540A (zh) * | 2021-01-11 | 2021-11-16 | 中国石油天然气股份有限公司 | 基于常规测井的非均质储层物性参数计算方法及装置 |
| CN112862261B (zh) * | 2021-01-18 | 2024-04-30 | 长江大学 | 考虑裂缝中游离气的页岩气井控制储量计算方法 |
| CN112986091B (zh) * | 2021-02-02 | 2021-09-07 | 东北石油大学 | 基于电阻率测试的层理缝发育页岩的含油饱和度评价方法 |
| CN112967147B (zh) * | 2021-02-04 | 2024-04-02 | 中海石油(中国)有限公司海南分公司 | 一种考虑多尺度裂缝的基岩产量贡献率计算方法 |
| CN112966380B (zh) * | 2021-03-10 | 2021-12-21 | 东北石油大学 | 一种确定陡坡砂体型汇聚脊油气储量的方法 |
| CN112814656B (zh) * | 2021-03-17 | 2022-07-22 | 成都理工大学 | 底水砂岩油藏开发大尺度高温高压模拟装置及模拟方法 |
| CN113553546B (zh) * | 2021-07-14 | 2023-08-18 | 中国海洋石油集团有限公司 | 一种电成像数据提取岩石连续胶结指数的方法、系统及计算机可读存储介质 |
| CN113719280A (zh) * | 2021-09-08 | 2021-11-30 | 陕西能源职业技术学院 | 一种低阻气层测井识别方法 |
| CN114037805B (zh) * | 2021-10-26 | 2022-07-15 | 重庆科技学院 | 一种非常规储层饱和度模型构建方法 |
| CN114417593B (zh) * | 2022-01-13 | 2023-06-06 | 东北石油大学 | 一种评价孔隙连通性对页岩采收率影响的实验方法 |
| CN114495432B (zh) * | 2022-01-19 | 2022-10-14 | 中国矿业大学 | 一种煤层顶底板含氢流体灾害的监测预警方法 |
| US20240019599A1 (en) * | 2022-07-15 | 2024-01-18 | Halliburton Energy Services, Inc. | Reservoir properties derived using ultra-deep resistivity inversion data |
| CN116502557B (zh) * | 2023-05-04 | 2024-01-30 | 西南石油大学 | 基于网格分析法的储层压裂后可动流体饱和度计算方法 |
| CN116882576B (zh) * | 2023-07-17 | 2024-03-08 | 北京福瑞宝能源科技有限公司 | 一种致密气砂岩气储层的产水概率预测方法 |
| CN117332301A (zh) * | 2023-10-17 | 2024-01-02 | 大庆油田有限责任公司 | 用于储层分类评价的水淹层解释方法 |
| CN117495085B (zh) * | 2023-10-31 | 2024-06-04 | 大庆油田有限责任公司 | 一种井位实施风险定量评价方法 |
| CN117386359B (zh) * | 2023-12-13 | 2024-03-08 | 中国石油集团川庆钻探工程有限公司 | 一种缝洞型碳酸盐岩储层孔隙度随钻定量评价方法 |
| CN118011509A (zh) * | 2024-04-09 | 2024-05-10 | 中国石油大学(华东) | 地层裂缝有效性评价方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4751646A (en) * | 1986-05-15 | 1988-06-14 | Restech, Inc. | Method for determining original saturations in a produced field |
| SU1795095A1 (en) * | 1990-05-28 | 1993-02-15 | Yurochko Aleksandr | Method to determine coefficients of open porosity and oil and gas saturation of pay bed |
| RU2018887C1 (ru) * | 1991-09-02 | 1994-08-30 | Западно-Сибирский научно-исследовательский институт геофизических методов разведки | Способ определения характера насыщения пластов-коллекторов |
| CN101413388A (zh) * | 2008-12-02 | 2009-04-22 | 大庆油田有限责任公司 | 获得油水同层原始含油饱和度的方法及估算未试油油水同层原始含油饱和度的方法 |
| RU2360233C1 (ru) * | 2007-12-19 | 2009-06-27 | Открытое акционерное общество "Томский научно-исследовательский и проектный институт нефти и газа Восточной нефтяной компании" ОАО "ТомскНИПИнефть ВНК" | Способ определения нефтенасыщенности породы |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU397636A1 (ru) * | 1971-08-18 | 1973-09-17 | Г. С. Степанова , А. М. Мамед заде Институт геологии , разработки горючих испокаемых | Способ определения давления насыщения пластовых жидкостей в пористой среде |
| US4517836A (en) * | 1983-10-25 | 1985-05-21 | Mobil Oil Corporation | Method for determining oil saturation in a subsurface formation |
| US4752882A (en) * | 1986-05-05 | 1988-06-21 | Mobil Oil Corporation | Method for determining the effective water saturation in a low-resistivity hydrocarbon-bearing rock formation based upon rock matrix conductance |
| CN1032663C (zh) | 1988-07-25 | 1996-08-28 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
| FR2844360B1 (fr) * | 2002-09-11 | 2004-10-15 | Inst Francais Du Petrole | Methode pour determiner l'indice de resistivite en fonction de la saturation en eau, de certaines roches de porosite complexe |
| ATE447192T1 (de) | 2004-12-17 | 2009-11-15 | Schlumberger Technology Bv | Verfahren zur bestimmung der wassersättigung einer untergrundformation |
| US7363161B2 (en) * | 2005-06-03 | 2008-04-22 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
| US7825659B2 (en) * | 2005-06-03 | 2010-11-02 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
| US8311789B2 (en) * | 2006-02-24 | 2012-11-13 | Saudi Arabian Oil Company | Monte Carlo simulation of well logging data |
| CN100588990C (zh) * | 2006-08-31 | 2010-02-10 | 中国石油大学(北京) | 地层流体饱和度评价方法 |
| US7565278B2 (en) * | 2006-12-04 | 2009-07-21 | Chevron U.S.A. Inc. | Method, system and apparatus for simulating fluid flow in a fractured reservoir utilizing a combination of discrete fracture networks and homogenization of small fractures |
| US8230918B2 (en) * | 2007-05-24 | 2012-07-31 | Saudi Arabian Oil Company | Method of characterizing hydrocarbon reservoir fractures in situ with artificially enhanced magnetic anisotropy |
-
2009
- 2009-06-22 CN CN200910087474A patent/CN101929973B/zh active Active
-
2010
- 2010-03-31 WO PCT/CN2010/000420 patent/WO2010148628A1/zh active Application Filing
- 2010-03-31 RU RU2011143258/28A patent/RU2523776C2/ru active
- 2010-03-31 EP EP10791135.6A patent/EP2447469B1/en active Active
- 2010-03-31 AU AU2010265746A patent/AU2010265746B2/en active Active
- 2010-03-31 US US13/379,706 patent/US8805661B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4751646A (en) * | 1986-05-15 | 1988-06-14 | Restech, Inc. | Method for determining original saturations in a produced field |
| SU1795095A1 (en) * | 1990-05-28 | 1993-02-15 | Yurochko Aleksandr | Method to determine coefficients of open porosity and oil and gas saturation of pay bed |
| RU2018887C1 (ru) * | 1991-09-02 | 1994-08-30 | Западно-Сибирский научно-исследовательский институт геофизических методов разведки | Способ определения характера насыщения пластов-коллекторов |
| RU2360233C1 (ru) * | 2007-12-19 | 2009-06-27 | Открытое акционерное общество "Томский научно-исследовательский и проектный институт нефти и газа Восточной нефтяной компании" ОАО "ТомскНИПИнефть ВНК" | Способ определения нефтенасыщенности породы |
| CN101413388A (zh) * | 2008-12-02 | 2009-04-22 | 大庆油田有限责任公司 | 获得油水同层原始含油饱和度的方法及估算未试油油水同层原始含油饱和度的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101929973A (zh) | 2010-12-29 |
| WO2010148628A1 (zh) | 2010-12-29 |
| US20120109603A1 (en) | 2012-05-03 |
| EP2447469A4 (en) | 2017-02-22 |
| AU2010265746B2 (en) | 2014-09-04 |
| CN101929973B (zh) | 2012-10-17 |
| US8805661B2 (en) | 2014-08-12 |
| AU2010265746A1 (en) | 2011-08-18 |
| EP2447469A1 (en) | 2012-05-02 |
| AU2010265746A2 (en) | 2011-08-25 |
| RU2011143258A (ru) | 2013-07-27 |
| EP2447469B1 (en) | 2018-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2523776C2 (ru) | Способ количественного расчета насыщенности трещинного коллектора углеводородами | |
| Dexter | Soil physical quality: Part III: Unsaturated hydraulic conductivity and general conclusions about S-theory | |
| CN109753755B (zh) | 一种确定储层含水饱和度的方法 | |
| US8598882B2 (en) | Method of monitoring a hydrocarbon reservoir | |
| CN108763711A (zh) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 | |
| Kuanzhi et al. | A reserve calculation method for fracture-cavity carbonate reservoirs in Tarim Basin, NW China | |
| CN104833789A (zh) | 利用致密砂岩微观孔隙结构确定气水关系的方法 | |
| CN106126936B (zh) | 一种致密低渗透储层裂缝有效性的综合评价方法 | |
| RU2455483C2 (ru) | Способ оценки степени трещиноватости карбонатных пород через параметр диффузионно-адсорбционной активности | |
| Du | Prediction of permeability and its anisotropy of tight oil reservoir via precise pore-throat tortuosity characterization and “umbrella deconstruction” method | |
| CN106285652B (zh) | 确定页岩游离气体饱和度的方法 | |
| CN110593857B (zh) | 油润湿致密砂岩核磁共振测井流体识别方法及其图版 | |
| CN109901238B (zh) | 一种基于应力差电阻率实验的高应力地层电阻率校正方法 | |
| CN105844011A (zh) | 一种基于毛管模型的渗透率计算方法 | |
| CN108255990A (zh) | 一种定量表征不同河型河道砂体的几何参数关系的方法 | |
| CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
| CN114035236A (zh) | 基于新三水模型的核磁和电阻率联合反演含水饱和度方法 | |
| CN116146170A (zh) | 一种三维感应各向异性泥页岩储层饱和度定量计算方法和装置 | |
| Alreshedan et al. | Sensitivity analysis of the effect of pore structure and geometry on petrophysical and electrical properties of tight media: random network modeling | |
| CN114441405B (zh) | 基于压实和胶结减孔趋势的次生增孔幅度定量评价方法 | |
| Adeboye et al. | PERMEABILITY ESTIMATION AND HYDRAULIC ZONE PORE STRUCTURES IDENTIFICATION USING CORE AND WELL LOGS DATA. | |
| Bertotti et al. | Predicting multiscale fracture patterns in buried reservoirs: The importance of outcrop data in a coherent workflow | |
| CN111428907B (zh) | 老油田开发甜点优选方法 | |
| RU2778354C1 (ru) | Способ и система оценивания запасов углеводородов в неоднородном пласте | |
| CN107976390B (zh) | 一种基于覆压孔渗实验确定平均孔喉比的计算方法 |