К оглавлению журнала

А.Г. Арье (ВНИГНИ)

Проблема определения фильтрационно-емкостных свойств в межскважинном пространстве подземной гидросферы методами сейсморазведки является кардинальной для поисков и разведки нефтяных и газовых залежей. Ей посвящено большое число публикаций как в отечественной, так и в мировой специальной литературе. К их числу можно отнести работы Н.Н. Пузырева (1959),Е.А. Копилевича, В.С. Славкина, Е.С. Шараповой (1988), Д.Н. Крылова (1992), А. Нура (1990), Д.В. Неффа (1993), Е.А. Копилевича (1995) и многие другие.

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

Впервые в своих построениях это обстоятельство попытался учесть Е.А.Копилевич [2]. Так, он считает, что в слабопроницаемых отложениях сейсмическая волна ведет себя как в непористом твердом теле, и на этом основании выделяет емкостную составляющую, относящуюся исключительно к пластам-коллекторам. Тем самым он косвенно показывает, что скорость прохождения сейсмической волны в принципе зависит от проницаемости пористой среды. В качестве основания для такого утверждения использован широко известный факт удержания слабопроницаемыми породами поровых жидкостей в связанном состоянии, по консистенции близком к твердому. Понятно, что предшественники, не учитывавшие этого обстоятельства, имели все основания считать, что скорость прохождения сейсмического возбуждения в геологической пористой среде определяется исключительно пористостью среды и не зависит от ее проницаемости [3].

Тем не менее задача оценки проницаемости в межскважинном геологическом пространстве не теряла своей актуальности и вследствие этого работа над ее решением продолжалась. Как правило, главным направлением исследований служила попытка прямого коррелирования параметров пористости и проницаемости. При этом исследователи почему-то упускали из виду, что рассматриваемая корреляция должна учитывать и величину удельной поверхности, как это следует из известного соотношения Козени — Кармана. Иногда такая корреляция реализуется, но вследствие приведенных соображений это, скорее, частный случай, особенно в тех ситуациях, когда коррелятивная связь оказывается прямолинейной. На основании изложенного это направление исследований нельзя признать теоретически обоснованным и, следовательно, перспективным для достижения поставленной цели.

Как упоминалось, скорость сейсмической волны в заполненной жидкостью пористой среде в традиционной постановке задачи, т.е. в неявном предположении о том, что движение поровой жидкости подчиняется закону Дарси, определяется пористостью среды и не зависит от ее проницаемости. Такое утверждение справедливо только для сред, характеризующихся относительно высокими значениями проницаемости [1]. В общем же случае такая постановка задачи оказывается недостаточно корректной, так как в ней не в полной мере учитывается влияние взаимодействия жидкой фазы с поверхностью скелета минеральной составляющей среды. Если же такое взаимодействие принять в расчет, то вывод о независимости скорости распространения упругой волны в водонасыщенной пористой среде от ее проницаемости оказывается недостаточно убедительным. Рассмотрим этот процесс, используя традиционный алгоритм анализа, применяемый обычно для решения такого рода задач.

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

где n — пористость образца; V_c, V_ж — скорость прохождения волны через материал скелета и жидкости соответственно.

Если обе части формулы (1) поделить на l, то в соответствии с определением понятия "интервальное время" получим

где t — интервальное время, характерное для всей системы в целом; t_c, t_ж - то же для материала скелета и жидкости соответственно.

Из теоретической физики (Китайгородский А.И., 1959) известно, что

Второе слагаемое целесообразно рассмотреть более детально. По определению пористость среды — константа, а вот величина βζ, по существу, являющаяся коэффициентом упругоемкости пласта (Щелкачев В.Н., 1959), может быть выражена через фильтрационные характеристики пласта, такие как коэффициент проницаемости, коэффициент пьезопроводности среды и вязкость фильтрующейся жидкости:

Такой вид выражение (6) имеет, если считать.что процесс фильтрации подчиняется закону Дарси. Если иметь в виду, что этот закон является частным случаем обобщенного закона фильтрации (Арье А.Г., 1981)

где V— скорость фильтрации; k₀ — коэффициент проницаемости по Дарси, т.е. при J>>J₀; J- градиент давлений, развиваемый источником колебаний; J₀ - начальный градиент фильтрации пористой среды.

Очевидно, что при J = J₀ вязкость жидкости возрастает и приближается по значению к вязкости твердого тела, так как закон справедлив только в интервале J₀<J<=J_кр

Этот, казалось бы, физически парадоксальный вывод означает, ЧТО при УСЛОВИИ J = j₀ поровая жидкость целиком переходит в связанное состояние и сжимаемость ее становится "неотличимой" от сжимаемости твердого тела *.

Подставляя выражение (10) в формулу (5), будем иметь:

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

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

где Сз — некоторая постоянная; R — усредненный радиус поровых каналов, определяющий значение проницаемости. Следовательно, окончательно можно записать, что

причем C₁, C₂, C₃ — постоянные величины.

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

В заключение представляется важным проанализировать физический смысл формулы (11) с учетом выражения для начального градиента, представленного уравнением (15). Из нее видно, что при достаточно низком начальном градиенте фильтрации, характерном для пород с относительно высокой проницаемостью, как это следует из (15), формула (11) вырождается в зависимость (5) или тождественную ей (1). Другими словами, при достаточно высокой проницаемости традиционные представления о зависимости интервального времени только от пористости породы достаточно справедливы; вернее, зависимостью этого показателя от проницаемости 'можно пренебречь без заметного ущерба для точности оценок.

Напротив, при весьма низкой проницаемости в соответствии с формулой (5) начальный градиент фильтрации становится соразмерным градиенту напора, вызванному сейсмическим возбуждением, и при jq = J второй член двучлена (11) обращается в нуль. Физически это означает, что поровый флюид жестко связан стенками пор и, следовательно, по деформационным характеристикам и агрегатному состоянию приближается к свойствам твердого тела, т.е. становится квазижестким. В такой ситуации пористость, фигурирующая в качестве слагаемого в первом члене выражения (11), становится "равной" нулю или мало от него отличной. Тогда формула (11) вырождается в соотношение

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

* На самом деле сжимаемость связанной жидкости не равна сжимаемости твердого тела; здесь имеется в виду, что βж << βж⁰