Способ определения коэффициента эффективной пористости горных пород в скважинах Российский патент 2023 года по МПК G01V11/00 G01V1/28 G01V1/30 

Изобретение относится к области нефтегазовой промышленности и предназначено для определения коэффициента эффективной пористости горных пород продуктивных и водонасыщенных интервалов и может найти применение при подсчете запасов и в гидродинамической модели пласта по глубине исследуемого разреза в скважине.

Коэффициент эффективной пористости Кпэф, характеризующий долю объема пор в котором происходит течение пластовых флюидов, является важной петрофизической характеристикой коллектора и используется при выделении пластов-коллекторов, при построении постоянно действующих геолого-промысловых моделей месторождений нефти и газа, при подсчете запасов углеводородов, при определении эффективных и фазовых проницаемостей, при гидродинамическом моделировании процессов разработки нефтяного или газового месторождения.

Известны способы определения трещинной пористости пород, основанные на предварительном формировании набора образцов исследуемой породы, экспериментальном определении общей пористости каждого из упомянутых образцов в атмосферных условиях, определении скорости распространения продольной волны и общей пористости в образцах исследуемой породы в условиях, моделирующих пластовые условия, расчете величины трещинной пористости для каждого из образцов исследуемой породы и определении поровой пористости, как разницы между общей пористостью и трещинной пористостью (RU 2516392, 2012 г., RU 2615051, 2015 г.).

Однако известные способы направлены на определение трещинной пористости и не обеспечивают получение информации об эффективной пористости пласта.

Известен способ определения коэффициента эффективной пористости на образцах керна, включающий высушивание образца до постоянной массы, насыщение под вакуумом моделью пластовой воды, создание в образце остаточной водонасыщенности путем вытеснения воды вытесняющим флюидом до прекращения выхода воды за счет центрифугирования при максимальных оборотах или за счет создания максимального давления в капилляриметре с полупроницаемой мембраной, при этом используют для измерений образец после экстракции растворителями или после экстракции неагрессивным агентом или промытый неэкстрагированный образец, полностью насыщенный моделью пластовой воды образец взвешивают в модели пластовой воды и в вытесняющем флюиде, после создания остаточной водонасыщенности образец взвешивают в вытесняющем флюиде и определяют коэффициент эффективной пористости по соотношению разницы веса образца, полностью насыщенного моделью пластовой воды в вытесняющем флюиде, и веса образца, насыщенного вытесняющим флюидом и остаточной водой, к разнице веса полностью насыщенного образца в вытесняющем флюиде и веса полностью насыщенного образца в пластовой воде (RU 2483291, 2011 г.).

Предлагаемый способ определения коэффициента эффективной пористости горных пород применим для случаев коллекторов с высоким содержанием разбухающей глинистой компоненты, а также образцов керна неэкстрагированных или экстрагированных с применением неагрессивных агентов.

Известен способ определения коэффициента эффективной пористости по данным метода ядерно-магнитного резонанса путем определения разности между пористостью и объемным содержанием воды, связанной с глинистым цементом в зависимости от содержания и типа цемента (US 5557200, 1996 г.). Указанный способ предусматривает проведение в скважине исследований методом ядерно-магнитного резонанса, а также специальной процедуры обоснования разделения сигнала от подвижной и связанной воды, что усложняет его реализацию.

Также известен способ определения коэффициента эффективной пористости по разности пористости и содержания глинистого цемента, при этом пористость определяют по зависимости ее изменения по глубине, которая формируется по результатам определения пористости в данном регионе. Глинистость определяется по данным гамма-метода или гамма-спектрометрии с использованием опорных пластов с максимальным и нулевым содержанием глинистого цемента и его плотности (US 4346592, 1982 г.).

Недостатком способа является отсутствие учета связанной воды матрицы (скелета) пласта, отсутствие учета неоднородности глинистого цемента, необходимость обоснования зависимости изменения пористости с глубиной, следствием чего является недостаточно высокая достоверность определения эффективной пористости.

Из известных технических решений наиболее близким к предлагаемому по технической сущности и достигаемому результату является способ определения коэффициента эффективной пористости пласта по данным геофизических исследований скважин (ГИС) и рассчитанных опорных характеристических значений физических параметров ГИС с использованием теоретических петрофизических моделей физических свойств коллектора для нескольких методов ГИС. Для расчета эффективной пористости проводится объединение парных сопоставлений методов ГИС с использованием оптимизационного алгоритма, основанного на ортонормализация Грама-Шмидта, требующего применения специализированного программного обеспечения (Коваленко К.В., Система петрофизического обеспечения моделирования залежей нефти и газа на основе эффективной пористости гранулярных коллекторов: специальность 25.00.10 "Геофизика, геофизические методы поисков полезных ископаемых": автореферат диссертации на соискание ученой степени доктора геолого-минералогических наук. - Москва, 2015. EDN ZPUYWN, стр. 182-202).

В указанном способе при расчете опорных характеристических значений физических параметров требуется наличие априорной информации о составе и свойствах обломочной и тонкодисперсной фракций пласта-коллектора, их объемном содержании, использование справочных значений физических свойств компонентов породы и не регламентируются процедуры специальных петрофизических исследований на образцах керна, что ограничивает контроль входных данных и приводит к снижению достоверности определения эффективной пористости.

Технической проблемой, на решение которой направлено предлагаемое изобретение, является повышение достоверности определения коэффициента эффективной пористости горных пород продуктивных и водонасыщенных интервалов информирование диаграммы объективных значений коэффициента эффективной пористости по глубине исследуемого разреза в скважине, пригодных для использования при подсчете запасов и в гидродинамической модели пласта и обеспечивающих учет диапазона изменения коэффициента граничной пористости.

Указанная проблема решается описываемым способом определения коэффициента эффективной пористости в скважинах, заключающимся в том, что предварительно на промытых неэкстрагированных контрастных по фильтрационно-емкостным свойствам образцах керна определяют коэффициенты пористости и остаточной водонасыщенности, значения объемной плотности и интервального времени пробега продольной акустической волны в условиях, моделирующих пластовые, проводят гранулометрический анализ образцов и определяют объемные содержания песчаной и алевритовой фракций в объеме породы, затем путем регрессионного анализа определяют параметр линейной связи между объемными содержаниями песчаной и алевритовой фракции образцов, по которому находят коэффициент максимальной пористости пласта и объемное содержание глинистой фракции образцов, и параметр линейной связи между объемным содержанием глинистой фракции и коэффициентом остаточной водонасыщенности образцов, по которому определяют значение коэффициента минимальной остаточной водонасыщенности, соответствующего коэффициенту максимальной пористости пласта, после чего по значению коэффициента максимальной пористости пласта и значению коэффициента минимальной остаточной водонасыщенности определяют коэффициент максимальной эффективной пористости пласта и коэффициент граничной пористости для каждого образца по следующей формуле:

где Кп, Кво - измеренные на керне коэффициенты пористости и остаточной водонасыщенности, Кпmax, К п эф max - коэффициенты максимальной пористости и максимальной эффективной пористости пласта, все коэффициенты в долях единицы,

путем анализа данных по полученной выборке устанавливают коэффициенты минимальной Кпгрmin и максимальной Кпгрmax граничной пористости пласта и путем регрессионного анализа определяют параметры C1, С2, С3, С4 линейных связей значений объемной плотности σ в г/см³ и интервального времени ΔT вмкс/м с коэффициентом пористости Кпв % по формулам:

где параметры С1 и С3- отношения изменений значений объемной плотности и интервального времени к изменению коэффициента пористости; С2 и С4 - значения объемной плотности в г/см³ и интервального времени в мкс/м при коэффициенте пористости, равном нулю, n - число образцов,

с использованием полученных параметров С1, С2, С3, С4 рассчитывают значения интервального времени ΔT₁ и объемной плотности σ₁, соответствующие коэффициенту максимальной пористости Кпmax пласта, значения интервального времени ΔТ₂ и объемной плотности σ₂, соответствующие коэффициенту минимальной граничной пористости пласта Кпгрmin, и значения интервального времени ΔТ₃ и объемной плотности σ₃, соответствующие коэффициенту максимальной граничной пласта Кпгрmin, далее измеряют в скважине по глубине со стандартным шагом квантовании значения объемной плотности σ плотностным гамма-гамма методом (ГГМ-П) и интервального времени пробега продольной акустической волны ΔT акустическим методом(АМ), и рассчитывают коэффициент эффективной пористости в каждой точке по глубине исследования пласта в скважине по формуле:

где σ - значения объемной плотности по данным ГГМ-П, ΔT - значения интервального времени пробега продольной акустической волны по данным AM, σ₁ и ΔT₁ - значения объемной плотности и интервального времени, соответствующие коэффициенту максимальной пористости Кпmax, σ₂ и ΔТ₂ - значения объемной плотности и интервального времени, соответствующие коэффициенту минимальной граничной пористости Кпгрmin, σ₃ и ΔТ₃ - значения объемной плотности и интервального времени, соответствующие коэффициенту максимальной граничной Кпгрmax, Кп эф max - коэффициент максимальной эффективной пористости пласта в %, значения объемной плотности выражаются в г/см³, значения интервального времени пробега продольной акустической волны в мкс/м.

Достигаемый технический результат заключается в оптимизации алгоритма определения коэффициента эффективной пористости пластана основе данных двух независимых методов исследования путем комплексирования результатов лабораторных исследований кернового материала по модели реликтовой остаточной воды и геофизических исследований скважин, вскрывающих межзерновые (гранулярные) породы-коллекторы и, как следствие, обеспечения повышения достоверности определения коэффициента эффективной пористости пласта за счет повышения объема контролируемой информации.

Сущность способа заключается в определении коэффициента эффективной пористости пласта по данным геофизических исследований скважин (ГИС) и рассчитанных опорных характеристических значений объемной плотности и интервального времени с использованием вычисленных параметров связей объемной плотности а и интервального времени с коэффициентом пористости по модели реликтовой остаточной воды и по установленным опорным характеристическим значениям коэффициентов пористости на образцах керна.

Способ осуществляют следующим образом.

Предварительно по керновому материалу определяют опорные характеристические коэффициенты пористости в исследуемом интервале.

Под характеристическими коэффициентами пористости понимаются коэффициент максимальной пористости Кпmax неглинистого пласта-коллектора, соответствующий ему коэффициент максимальной эффективной пористости Кпэфmax и коэффициенты граничной пористости пласта Кпгрmin и Кпгрmax при значениях Кво=100% или нулевых значениях коэффициента эффективной пористости.

Исследования проводят на промытых неэкстрагированных контрастных по фильтрационно-емкостным свойствам образцах.

При этом используют образцы керна с использованием модели реликтовой остаточной воды, физические свойства которой определяют в интервалах коллекторов с наихудшими фильтрационно-емкостными свойствами, в которых реликтовая остаточная вода преобладает.

Количество образцов выбирают согласно действующему РД 39-0147716-505-85 «Порядок отбора, привязки, хранения, движения, и комплексного исследования керна и грунтов нефтегазовых скважин», в соответствие с которым частота отбора образцов керна для комплексных лабораторных исследований определяется литологическим составом, изменчивостью физических свойств и характером насыщения изучаемых пород.

На отобранных образцах керна измеряют коэффициенты пористости Кп и остаточной водонасыщенности Кво, объемную плотность σ и интервальное время пробега продольной акустической волны ΔT в условиях, моделирующих пластовые.

Проводят гранулометрический анализ отобранных образцов и определяют объемные содержания песчаной и алевритовой фракций в объеме породы.

Затем путем регрессионного анализа определяют параметр линейной связи S1 в % между объемными содержаниями песчаной и алевритовой фракции образцов по формуле:

где Кпесч, Калев - объемные содержаний песчаной и алевритовой фракций в объеме породы в %, n - число образцов.

По рассчитанному параметру S1 находят коэффициент максимальной пористости пласта по формуле:

Кпmax=100-S1 (2)

и объемное содержание глинистой фракции образцов Кглв % по формуле:

Кгл=Кпmax-Кп (3)

Путем регрессионного анализа определяют параметр линейной связи S2 в % между объемным содержанием глинистой фракции и коэффициентом остаточной водонасыщенности образцов Кво по формуле:

где Кгл - объемное содержание глинистой фракции в %, Кво- коэффициент остаточной водонасыщенности %, n - число образцов, по которому определяют значение минимального коэффициента остаточной водонасыщенности Квоmin=S2 в %, соответствующего коэффициенту максимальной пористости пласта.

После чего рассчитывают коэффициент максимальной эффективной пористости Кп эф max в % по формуле:

Кп эф max=Кпmax(1-Квоmin/100) (5)

Далее для каждого образца определяют коэффициент граничной пористости по формуле:

где Кп, Кво - измеренные на керне коэффициенты пористости и остаточной водонасыщенности, Кпmax, Кпэфmax - коэффициенты максимальной пористости и максимальной эффективной пористости пласта, все коэффициенты в долях единицы.

Проводят анализ полученных данных по всем образцам и выбирают коэффициенты минимальной Кпгрmin и максимальной Кпгрmax граничной пористости пласта.

Затем проводят регрессионный анализ полученных в ходе лабораторных измерений на образцах керна значений объемной плотности и интервального времени и определяют параметры C1, С2, С3, С4 линейных связей этих параметров с коэффициентом пористости по формулам:

где параметры С1 и С3 - отношения изменений значений объемной плотности и интервального времени к изменению коэффициента пористости; С2 и С4 - значения объемной плотности в г/см³ и интервального времени в мкс/м при коэффициенте пористости, равном нулю, n - число образцов,

Используя полученные параметры C1, С2, С3, С4 определяют:

- значения интервального времени ΔT₁ и объемной плотности σ₁, соответствующие коэффициенту максимальной пористости пласта Кпmax:

σ₁=С1 ⋅ Кпmax+С2 (11)

ΔT₁=С3 ⋅ Кпmax+С4, (12)

(при необходимости значения корректируют за влияние фильтрата промывочной жидкости и остаточной нефте- или газонасыщенности с использованием петрофизической модели неглинистого коллектора и стандартных методик интерпретации);

- значения интервального времени ΔТ₂ и объемной плотности ⋅σ₂, соответствующие коэффициенту минимальной граничной пористости пласта Кпгрmin:

σ₂=C1 ⋅ Кпгрmin+C2 (13)

ΔT₂=С3 ⋅ Кпгрmin+C4; (14

- значения интервального времени ΔT₃ и объемной плотности ⋅σ₃, соответствующие коэффициенту максимальной граничной пористости пласта Кпгрmax:

σ₃=С1 ⋅ Кпгрmax+С2 (15)

ΔТ₃=С3 ⋅ Кпгрmax+С4 (16)

Далее проводят геофизические измерения в скважине (ГИС)в исследуемом интервале вскрытого скважиной коллектора и измеряют по глубине со стандартным шагом квантования значения объемной плотности а плотностным гамма-гамма методом (ГГМ-П) и интервального времени пробега продольной акустической волны ΔT акустическим методом(АМ).

Значения коэффициента эффективной пористости в каждой точке по глубине исследования в скважине находят по формуле:

Приведенная формула позволяет рассчитать коэффициент эффективной пористости с использованием значений объемной плотности ⋅σ и интервального времени ΔT по данным ГИС и параметров ΔТ_1-3, σ_1-2, Кпэф max, полученных по результатам лабораторных исследований керна, т.е. рассчитать среднее арифметическое значение эффективной пористости по данным двух независимых методов исследования, линейно связанных с пористостью, с использованием средних значений диапазона изменения физических параметров при нулевой эффективной пористости.

Ниже приведен пример осуществления способа.

Сформировали набор из 27 образцов песчаника, измерили коэффициенты пористости Кп, остаточной водонасыщенности Кво и значения объемной плотности σ. Сформировали набор из 6 образцов и измерили значения интервального времени пробега продольной акустической волны ΔТв условиях, моделирующих пластовые, результаты измерений представлены в таблицах 1 и 2.

Сформировали набор из 19 образцов песчаника, провели гранулометрический анализ с определением объемных содержаний песчаной и алевритовой фракций в объеме породы. Результаты анализа представлены в таблице 3.

Рассчитали путем регрессионного анализа параметр линейной связи S1 по формуле (1):

Рассчитали коэффициент максимальной пористости пласта по формуле (2):

Кпmax=100-S1=100-73.49=26.51%.

Рассчитали для каждого образца объемное содержание глинистой фракции в объеме породы по формуле (3), результат показан в таблице 1.

Рассчитали путем регрессионного анализа соответствующий коэффициенту Кпmax коэффициент остаточной водонасыщенности Квоmin=S2 (параметр линейной связи S2) по формуле (4):

- Рассчитали коэффициент максимальной эффективной пористости пласта по формуле (5):

Кп эф max=Кпmax(1-Квоmin/100)=26.51(1-0.0212)=25.95%.

- С использованием формулы (6) рассчитали для каждого образца коэффициент граничной пористости, результат показан в таблице 1.

Установили путем анализа данных по полученной выборке коэффициенты минимальной Кпгрmin=1.49% и максимальной Кпгрmax=8.08% граничной пористости пласта.

- Определили путем регрессионного анализа по формулам (7-10) параметры C1, С2, С3, С4 линейных связей значений объемной плотности σ₁ и интервального времени ΔT₁ с коэффициентом пористости:

- Рассчитали по формулам (11-12) с использованием полученных параметров C1, С2, С3, С4 значения интервального времени ΔТ₁ и объемной плотности σ₁, соответствующие коэффициенту максимальной пористости Кпmax пласта:

σ₁=С1 ⋅ Кпmax+С2=-0.016 ⋅ 26.51+2.615=2.19 г/см³

ΔT₁=С3 ⋅ Кпmax+С4=5.3786 ⋅ 26.51+148.18=290.78 мкс/м,

По формулам (13-14) определили значения интервального времени ΔТ₂ и объемной плотности σ₂, соответствующие коэффициенту минимальной граничной пористости Кпгрmin пласта:

σ₂=С1 ⋅ Кпгрmin+С2=-0.016 ⋅ 1.49+2.615=2.59 г/см³

ΔТ₂=С3 ⋅ Кпгрmin+С4=5.3786 ⋅ 1.49+148.18=156.20 мкс/м

и по формулам (15-16) значения интервального времени ΔТ₃ и объемной плотности σ₃, соответствующие коэффициенту максимальной граничной Кпгрmax пласта:

σ₃=C1 ⋅ Кпгрmax+C2=-0.016 ⋅ 8.08+2.615=2.49 г/см³

ΔT₃=С3 ⋅ Кпгрmax+С4=5.3786 ⋅ 8.08+148.18=191.64 мкс/м.

Измерили в скважине по глубине со стандартным шагом квантования значения объемной плотности σ плотностным гамма-гамма методом (ГГМ-П) и интервального времени пробега продольной акустической волны ΔT акустическим методом(АМ).

- Рассчитали коэффициент эффективной пористости в каждой точке по глубине исследования пластав скважине по формуле (17).

Например, измеренные в скважине на глубине ХХ48.0 м аппаратурой ГГМ-П значения объемной плотности составляют 2.34 г/см³ и измеренные аппаратурой AM значения интервального времени составляют 224.73 мкс/м. Тогда, по формуле (17):

Аналогично провели расчет коэффициента эффективной пористости на каждом кванте глубины исследования.

В таблице 4 приведены результаты измерений значений объемной плотности σ и интервального времени ΔT в скважине и результат расчета коэффициента эффективной пористости Кп эф.

Для сравнения рассчитали коэффициент эффективной пористости по результатам лабораторных исследований керна по формуле Кпэф=Кп (1-Кво/100).

На фиг. 1 приведены результаты измерений значений объемной плотности σ и интервального времени ΔT в скважине методами ГГМ-П и АМ и сравнение расчетов коэффициента эффективной пористости по ГИС (линия) с расчетами коэффициента эффективной пористости по результатам лабораторных исследований керна по формуле Кпэф=Кп (1-Кво/100)(маркеры).

Таким образом, предлагаемый способ определения коэффициента эффективной пористости горных пород с использованием данных лабораторных измерений на образцах керна по модели реликтовой остаточной воды, обоснованием опорных характеристических параметров и использованием геофизических измерений в скважине (методы ГГМ-П и AM) позволяет повысить достоверность определения коэффициента эффективной пористости за счет комплексирования данных двух независимых методов исследования скважин, вскрывающих межзерновые (гранулярные) породы-коллекторы, включая коллекторы с высоким содержанием разбухающей глинистой или полиминеральной компоненты и неоднородных по величине остаточной водонасыщенности.

Предлагаемый способ позволяет рассчитать достоверные значения коэффициента эффективной пористости неоднородного пласта-коллектора для любого соотношения между коэффициентами пористости и остаточной водонасыщенности, в том числе, когда корректные статистические связи Кп(Кво) и Кпэф(Кп) отсутствуют.

Изобретение обеспечивает формирование диаграммы объективных значений коэффициента эффективной пористости по глубине исследуемого разреза в скважине, пригодных для использования при подсчете запасов и в гидродинамической модели пласта.

Изобретение относится к области нефтегазовой промышленности, предназначено для определения коэффициента эффективной пористости горных пород продуктивных и водонасыщенных интервалов и может найти применение при подсчете запасов и в гидродинамической модели пласта по глубине исследуемого разреза в скважине. Способ определения коэффициента эффективной пористости горных пород в скважинах заключается в том, что предварительно на промытых неэкстрагированных контрастных по фильтрационно-емкостным свойствам образцах керна определяют коэффициенты пористости и остаточной водонасыщенности, значения объемной плотности и интервального времени пробега продольной акустической волны в условиях, моделирующих пластовые. Проводят гранулометрический анализ образцов и определяют объемные содержания песчаной и алевритовой фракций в объеме породы. Затем путем регрессионного анализа определяют параметр линейной связи между объемными содержаниями песчаной и алевритовой фракции образцов, по которому находят коэффициент максимальной пористости пласта и объемное содержание глинистой фракции образцов, и параметр линейной связи между объемным содержанием глинистой фракции и коэффициентом остаточной водонасыщенности образцов, по которому определяют значение коэффициента минимальной остаточной водонасыщенности, соответствующего коэффициенту максимальной пористости пласта. После чего по значению коэффициента максимальной пористости пласта и значению коэффициента минимальной остаточной водонасыщенности определяют коэффициент максимальной эффективной пористости пласта и коэффициент граничной пористости для каждого образца. Путем анализа данных по полученной выборке устанавливают коэффициенты минимальной К_{п гр min} и максимальной К_{п гр max} граничной пористости пласта и путем регрессионного анализа определяют параметры C1, С2, С3, С4 линейных связей значений объемной плотности σ в г/см³ и интервального времени ΔT в мкс/м с коэффициентом пористости К_пв, %. Далее измеряют в скважине по глубине со стандартным шагом квантования значения объемной плотности σ плотностным гамма-гамма методом (ГГМ-П) и интервального времени пробега продольной акустической волны ΔT акустическим методом (AM) и рассчитывают коэффициент эффективной пористости в каждой точке по глубине исследования пласта в скважине. Технический результат заключается в оптимизации алгоритма определения коэффициента эффективной пористости пласта на основе данных двух независимых методов исследования путем комплексирования результатов лабораторных исследований кернового материала по модели реликтовой остаточной воды и геофизических исследований скважин, вскрывающих межзерновые (гранулярные) породы-коллекторы, и, как следствие, обеспечения повышения достоверности определения коэффициента эффективной пористости пласта за счет повышения объема контролируемой информации. 1 ил.

Способ определения коэффициента эффективной пористости горных пород в скважинах, заключающийся в том, что предварительно на промытых неэкстрагированных контрастных по фильтрационно-емкостным свойствам образцов керна определяют коэффициенты пористости и остаточной водонасыщенности, значения объемной плотности и интервального времени пробега продольной акустической волны в условиях, моделирующих пластовые, проводят гранулометрический анализ образцов и определяют объемные содержания песчаной и алевритовой фракций в объеме породы, затем путем регрессионного анализа определяют параметр линейной связи между объемными содержаниями песчаной и алевритовой фракции образцов, по которому находят коэффициент максимальной пористости пласта и объемное содержание глинистой фракции образцов, и параметр линейной связи между объемным содержанием глинистой фракции и коэффициентом остаточной водонасыщенности образцов, по которому определяют значение коэффициента минимальной остаточной водонасыщенности, соответствующего коэффициенту максимальной пористости пласта, после чего по значению коэффициента максимальной пористости пласта и значению коэффициента минимальной остаточной водонасыщенности определяют коэффициент максимальной эффективной пористости пласта и коэффициент граничной пористости для каждого образца по следующей формуле:

где К_п, К_во - измеренные на керне коэффициенты пористости и остаточной водонасыщенности, К_{п max}, К_{п эф max} - коэффициенты максимальной пористости и максимальной эффективной пористости пласта, все коэффициенты в долях единицы, путем анализа данных по полученной выборке устанавливают коэффициенты минимальной К_{п гр min} и максимальной К_{п гр max} граничной пористости пласта и путем регрессионного анализа определяют параметры C1, С2, С3, С4 линейных связей значений объемной плотности σ в г/см³ и интервального времени ΔT в мкс/м с коэффициентом пористости К_пв % по формулам:

,

где параметры С1 и С3 - отношения изменений значений объемной плотности и интервального времени к изменению коэффициента пористости; С2 и С4 - значения объемной плотности в г/см³ и интервального времени в мкс/м при коэффициенте пористости, равном нулю, n - число образцов, с использованием полученных параметров C1, С2, С3, С4 рассчитывают значения интервального времени ΔT₁ и объемной плотности σ₁, соответствующие коэффициенту максимальной пористости К_{п max} пласта, значения интервального времени ΔТ₂ и объемной плотности σ₂, соответствующие коэффициенту минимальной граничной пористости пласта К_{п. гр min}, и значения интервального времени ΔТ₃ и объемной плотности σ₃, соответствующие коэффициенту максимальной граничной пористости пласта К_{п гр min}, далее измеряют в скважине по глубине со стандартным шагом квантования значения объемной плотности σ плотностным гамма-гамма методом (ГГМ-П) и интервального времени пробега продольной акустической волны ΔT акустическим методом (AM) и рассчитывают коэффициент эффективной пористости в каждой точке по глубине исследования пласта в скважине по формуле:

,

где σ - значения объемной плотности по данным ГГМ-П, ΔT - значения интервального времени пробега продольной акустической волны по данным AM, σ₁ и ΔT₁ - значения объемной плотности и интервального времени, соответствующие коэффициенту максимальной пористости К_{п max}, σ₂ и ΔТ₂ - значения объемной плотности и интервального времени, соответствующие коэффициенту минимальной граничной пористости К_{п гр min}, σ₃ и ΔТ₃ - значения объемной плотности и интервального времени, соответствующие коэффициенту максимальной граничной пористости К_{п гр max}, К_{п эф max} - коэффициент максимальной эффективной пористости пласта в %, значения объемной плотности выражаются в г/см³, значения интервального времени пробега продольной акустической волны в мкс/м.