Результат интеллектуальной деятельности: СПОСОБ ОПРЕДЕЛЕНИЯ РАВНОВЕСНОЙ СМАЧИВАЕМОСТИ ПОВЕРХНОСТИ РАЗДЕЛА ПУСТОТНОГО ПРОСТРАНСТВА И ТВЕРДОЙ ФАЗЫ ОБРАЗЦА ГОРНОЙ ПОРОДЫ

ОБЛАСТЬ ТЕХНИКИ

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

УРОВЕНЬ ТЕХНИКИ

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

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

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

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

В большинстве отраслей промышленности основным методом, используемым для оценки смачиваемости, является метод измерения контактного угла между поверхностью твердого тела и межфазной поверхностью двух смачивающих жидкостей (см., например, W. Abdallah et al., Fundamentals of Wettability, Oilfield review, vol. 19, no. 2, pp. 44-61, 2007). Недостатком этого метода является длительное время, необходимое для установления равновесного угла смачивания (до 1000 ч.), гистерезис контактного угла, обусловленный множеством причин, таких как, но не ограничиваясь ими: поверхностная структура, поверхностные неоднородности, и др. Другим недостатком метода является то, что метод применим для оценки смачиваемости гладких однородных поверхностей и его довольно трудно реализовать для характеризации поверхностей пористых сред и гетерогенных поверхностей.

В нефтяной и газовой промышленности смачиваемость сильно смачиваемого водой или сильно смачиваемого нефтью резервуара в некоторых случаях может быть определена по результатам скважных геофизических исследований. В большинстве же случаев, для исследования смачиваемости используются лабораторные исследования керна горной породы. Методы USBM и метод Амота-Харви или их комбинации считаются отраслевыми стандартами для оценки смачиваемости керна (см., например, Е. Donaldson and W. Alam, Wettability, Houston: Gulf Publising Company, 2008, pp. 2-27). Оба перечисленных метода основаны на лабораторном моделировании процессов вытеснения нефти водой и воды нефтью, происходящих в пласте. В результате экспериментов по определению смачиваемости данными методами, как правило для каждого образца керна находят лишь одно число, характеризующее смачиваемость: между -1 и +1 для метода Амота и между -∞ и +∞ для USBM. Это число не характеризует смачиваемость напрямую, скорее оно характеризует эффективность вытеснения жидкостей в указанных процессах и в связи с тем, что эффективность вытеснения зависит от смачиваемости этот параметр имеет отношение к состоянию смачиваемости. В то же время полученное число не позволяет определить ни контактный угол смачивания жидкости на поверхности и не дает никакой информации о распределении смачиваемости на поверхности керна. Сложная внутренняя структура пор оказывает дополнительное влияние на результаты тестов вне зависимости от состояния смачиваемости поверхностей.

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

РАСКРЫТИЕ ИЗОБРЕТЕНИЯ

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

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

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

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

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

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

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

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

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

КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ

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

ПОДРОБНОЕ ОПИСАНИЕ ИЗОБРЕТЕНИЯ

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

Так как образовавшиеся адсорбированные пленки нефтяных компонент прочно связаны с поверхностью и, с большой вероятностью, не десорбируются при разработке месторождения (часто не десорбируются, даже с помощью упомянутых выше методов очистки керна), для характеризации смачиваемости горной породы необходимо определить смачиваемость породы покрытой пленками данных отложений, а не непосредственно смачиваемость твердой фазы породы W⁰(A). Именно это состояние поверхности горной породы (толщиной пленок нефтяных отложений мы пренебрегаем) мы будем называть равновесным состоянием. Под определением равновесной смачиваемости мы понимаем нахождение параметра смачиваемости W(A), для каждой точки поверхности А породы, соответствующего равновесному состоянию в пласте.

Как показано на фиг. 1 и 2, предлагаемый способ включает получение изображения (сканирование) образца 1 горной породы и получение 2 трехмерного изображения 3 образца 1. Сканирование может быть осуществлено, например, методом компьютерной рентгеновской микротомографии. Метод рентгеновской томографии позволяет определять пространственное распределение линейных коэффициентов ослабления (ЛКО) рентгеновского излучения внутри образца. В ходе эксперимента исследуемый образец просвечивают рентгеновским пучком с разных сторон. Прошедшее излучение (теневая проекция) при этом каждый раз регистрируется двумерным рентгеновским детектором. Далее, полученный массив теневых проекций восстанавливают каким-либо из известных алгоритмов реконструкции (например, обратное проецирование, алгебраические методы). Восстановленные данные являются послойным представлением внутренней структуры образца.

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

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

Затем проводят процедуру обработки 4 полученного изображения 3 образца 1, которая состоит в выделении пустотного пространства и твердой фазы, а также определении поверхности 5 раздела пустотного пространства и твердой фазы (пора-порода) образца.

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

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

Пусть А - ориентированная поверхность раздела пустотного пространства и твердой фазы исследуемого образца горной породы с нормалью, направленной от породы к поре. В общем случае для описания поверхности А может быть использована любая из известных систем координат трехмерного пространства (например, прямоугольная декартова (x, y, z)). Для определения А по полученному трехмерному изображению внутренней структуры пористого образца горной породы, выраженному в значениях коэффициента ослабления, необходимо выполнить следующие операции: сегментация трехмерного изображения (бинаризация); выделение точек трехмерного изображения (вокселей), соответствующих границам раздела пора-порода; определение нормалей к поверхности в каждом вокселе.

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

Сегментация трехмерного изображения может быть осуществлена любым из известных методов, наиболее подходящим в каждом конкретном случае в зависимости от качества данных и структуры объекта исследования. Один из возможных вариантов представлен в заявке WO/2014/003596. В результате должен быть получен трехмерный массив (набор двумерных изображений), в котором в случае бинаризации значение «0» соответствует поровому пространству, а значение «1» соответствует породе. Вместо «1» может быть взято любое другое заранее установленное число. Например, для 8-битных изображений часто используют значение «255». Не ограничивая общности, далее всюду используется значение «1».

Выделение граничных вокселей необходимо для создания дискретного представления искомой поверхности 5. Эта операция может быть выполнена несколькими способами. Например, в ходе выполнения сегментации воксели, которые имели близкие вероятности принадлежности к обоим классам (пора «0» и порода «1»), могут быть классифицированы как граничные. Иная возможность состоит в определении границы известным морфологическим методом «эрозия». Суть операции состоит в том, что всем вокселям со значением «1», отстоящим от какого-либо из вокселей со значением «0» на расстояние меньшее наперед заданного d присваивают значение «0». Далее, полученный массив может быть поэлементно вычтен из исходного, сегментированного. В результате будет получен массив со значением «1» на вокселях, принадлежащих границе пора-порода толщиной d. Представляется оптимальным использовать d=1. В описанном варианте исходные несегментированные изображения не требуются.

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

Затем проводят минеральный анализ 7, т.е. определяют распределение 8 минералов на поверхности 5 раздела пустотного пространства и твердой фазы образца. Различные минералы обладают разными свойствами смачиваемости. По этой причине определение свойств смачиваемости поверхности пустотного пространства требует учета ее минерального состава. Обозначим через М(А) соответствие между точками на поверхности А и их минеральным составом М в малой окрестности каждой точки. Таким образом, значения М должны однозначно указывать на тип минерала соответствующей точки. Технически это удобно реализовать в виде таблицы целых чисел, каждое из которых кодирует тот или иной тип минерала. Далее эти числа присваиваются всем точкам поверхности А в зависимости от их минеральной принадлежности. В простейшем случае мономинеральной поверхности М(А)=М, где М соответствует, например, кварцу, доломиту и т.д. В более сложных случаях возникает необходимость определения М(А), что может быть осуществлено различными способами. Ясно, что в силу ограниченности пространственного разрешения полученного трехмерного изображения, в одной точке (вокселе) поверхности А может находиться несколько минералов (смесь минералов). В таком случае представляется целесообразным выделять подобные смеси в отдельные условные типы минералов и все дальнейшие процедуры проводить для них, как для независимых минералов. Если количество таких точек на поверхности А незначительно, то введение дополнительных условных минералов, соответствующих смесям реальных минералов, не имеет смысла и лишь усложняет осуществление изобретения.

Распределение М(А) можно определить напрямую из данных рентгеновской микротомографии, а именно реконструированных значений ЛКО и их сравнения с теоретическими. Величина ЛКО для каждого вещества (минерала, в частности) зависит от его химического состава, массовой плотности и энергии зондирующего рентгеновского пучка. В случае известного энергетического спектра пучка тип минерала в точке может быть приблизительно определен путем корреляции реконструированного и теоретически рассчитанного значений ЛКО (см., например, J.Н. Hubbel and S.S.M, "NIST X-ray Attenuation Databases," The National Institute of Standards and Technology (NIST), [Online]. Available: http://www.nist.gov/pml/data/xraycoef/index.cfm.) Некоторые минералы имеют слишком близкие значения ЛКО и, как правило, неразличимы.

Возможно определение М(A) из данных рентгеновской микротомографии, дополненных результатами энерго-дисперсионного рентгеновского микроанализа какой-либо области (сечения) исследуемого образца (Energy Dispersive X-Ray microanalysis, EDX), выполнимого в системах сканирующей электронной микроскопии (Scanning Electron Microscopy, SEM) (см., например. G. Ying, "Automated Scanning Electron Microscope Based Mineral Liberation Analysis," Journal of Minerals & Materials Characterization & Engineering, vol. 2, no. 1, pp. 33-41, 2003). В методе EDX-SEM происходит регистрация характеристического рентгеновского спектра (возбуждаемого электронным пучком микроскопа) каждой точки сканируемой поверхности образца. Для наилучшего результата рекомендуется подготовить плоскую поверхность на специализированных станках в соответствии с требованиями метода EDX-SEM (заливка эпоксидной смолой исходного образца, разрезывание в желаемом месте, шлифовка, полировка и напыление токопроводящего материала). Полученные спектры могут быть проинтерпретированы и представлены в виде количественного химического состава. Для дальнейшего определения распределения минералов на сканированной поверхности необходимо использовать библиотеку спектров минералов, которую потребуется создать самостоятельно. Иная возможность состоит в использовании коммерчески доступного программного обеспечения для проведения минеральной идентификации. Хотя это всего лишь информация о поверхности (или некоторого сечения), полученное распределение минералов может быть представительным, для того чтобы распространить распределение на трехмерное изображение микротомографии и/или напрямую на А (например, на основании реконструированных значений ЛКО с предварительным пространственным совмещением распределения минералов и 3D изображения микротомографии).

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

Существуют и иные методы, которые могут быть использованы для определения объемного содержания минералов и их распределения 8 на поверхности М(А), в частности. К таким методам могут быть отнесены, например: синхротронная микротомография с использованием монохроматического рентгеновского излучения, рентгеновская флуоресцентная микротомография, рентгеновская топо-томография.

Помимо распределения минералов, для более корректного определения свойств смачиваемости на поверхности образца могут быть приняты во внимание топографические особенности поверхности. Эти особенности могут быть учтены через процедуру 9 определения параметров локальной шероховатости поверхности. В общем случае, шероховатость поверхности должна описываться целым набором параметров. Среди них наиболее важным для определения смачиваемости является фактор шероховатости 10, который вычисляется на поверхности А:

где S_True(A) - действительная площадь поверхности раздела пора-порода в некоторой малой окрестности рассматриваемой точки, a S_CT - площадь поверхности А, построенной по данным микротомографии, вычисленная в той же окрестности той же точки. Фактор шероховатости (1) отражает разницу в площадях поверхности, вычисленных на А и с помощью других методов более высокого разрешения.

Процедура 9 оценки значений 10 параметра F(A) в каждой точке поверхности А может быть выполнена разными способами. Например, возможны прямые измерения любым доступным методом топографии поверхности для каждого минерала, входящего в состав образца. Далее, производится расчет фактора шероховатости для каждого из них F(M). Наконец, с учетом ранее определенного распределения М(А) несложно рассчитать F(A).

Фактор шероховатости для конкретного минерала может быть приближенно оценен на основании литературных данных. Для определенной ранее поверхности А раздела пустотного пространства и твердой фазы исследуемого пористого образца несложно вычислить полную площадь (площадь поверхности образца, включая пустотное пространство). Вычисленную таким образом площадь можно сравнить с аналогичной величиной, измеренной лабораторно на этом же (желательно) или подобном образце. Площадь поверхности образца горной породы может быть измерена, например, широко известным методом BET (см., например, S. Brunauer, P. Emmett and Е. Teller, "Adsorption of Gases in Multimolecular Layers," J. Am. Chem. Soc, vol. 60, p. 309, 1938) или любым другим. Далее, фактор шероховатости выбирается таким, чтобы найденные площади удовлетворяли (1).

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

Затем определяют 11 смачиваемость твердой фазы в каждой точке образца с учетом полученного поверхностного распределения минералов М(А) и фактора шероховатости F(A). Определение 11 смачиваемости 12 твердой фазы горной породы W⁰(A) обозначает, что для каждой точки поверхности А определена величина:

где индексы S (solid), О (oil), W (water) обозначают соответственно поверхность горной породы, нефть и воду, γ^XY - поверхностное натяжение на границе между фазами Х и Y.

Для определения W⁰(A) могут быть использованы несколько способов.

Возможно измерение смачиваемости W(M) поверхностей минералов М, составляющих исследуемую поверхность, нефтью и пластовой водой. Смачиваемость каждого минерала может быть измерена отдельно с помощью методов измерения контактного угла, калориметрии, газовой адсорбции или любого другого известного метода измерения смачиваемости. Далее, зная распределение минералов М(А) по поверхности А, несложно определить W⁰(A).

В случае, когда выделение каждого минерала для исследований его свойств смачиваемости данными нефтью и водой не представляется возможным или связано со значительными сложностями возможен следующий вариант определения W⁰(A). Смачиваемость минералов, составляющих поверхность, определяется с помощью компонентной теории поверхностного натяжения (surface tension component theory). При этом для определения компонент поверхностного натяжения твердой поверхности каждого минерала могут быть использованы литературные данные, например, K. Mittal, Acid-Base Interactions: Relevance to Adhesion Science and Technology, Mittal, Ed., 2000, pp. 317-347, а для определения компонентных составляющих поверхностного натяжения исследуемых нефти и воды могут быть использованы экспериментальные измерения компонентных составляющих с помощью подхода трех жидкостей и измерений межфазных поверхностных натяжений по методу сидящей, висящей или вращающейся капли, либо любой другой известный метод измерения поверхностных натяжений жидкость-жидкость. Далее, как и ранее, зная распределение минералов М(А) по поверхности А, несложно определить W⁰(A).

В случае, когда подобные исследования межфазных натяжений жидкостей не представляются возможными, оценка смачиваемости каждого минерала производится на основе доступных литературных данных о смачиваемости минералов наиболее близкими по свойствам минеральным раствором и нефтью. Далее, зная распределение минералов М(А) по поверхности А, можно определить W⁰(A).

Как упоминалось выше, под определением равновесной смачиваемости 16 мы понимаем нахождение параметра смачиваемости W(A), для каждой точки поверхности А породы, соответствующему равновесному состоянию в пласте, после формирования нефтегазового месторождения. По этой причине далее предлагается провести гидродинамическое моделирование 15, а именно, моделировать процесс миграции нефти в пустотное пространство породы, заполненное пластовой водой на начальной стадии формирования месторождения, в соответствии с нашими знаниями и/или предположениями об истории этого процесса в данном пласте. В ходе моделирования система приходит в равновесное состояние. На основании получившегося равновесного распределения флюидов (нефть, пластовая вода) определяется итоговое распределение 16 параметра смачиваемости W(A). Как будет показано ниже, проведение данного моделирования осуществляется на построенной модели 6 с использованием найденных распределений смачиваемости 12 и фактора шероховатости 10.

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

Моделирование обоих сценариев миграции нефти может быть осуществлено разными способами. Например, методом расчета распределения жидкостей в пространстве пористой среды на основе минимизации полной свободной энергии Гельмгольца. Подобные расчеты могут быть выполнены с помощью метода функционала плотности в гидродинамике (A.Y. Demianov, О.Y. Dinariev and N.V. Evseev, Introduction to the density functional method in hydrodynamics, Moscow: Fizmatlit, 2014, pp. 37-70). Данный метод позволяет моделировать гидродинамические и диффузионные процессы в многокомпонентных жидкостях, протекающие в сложных трехмерных геометриях.

В случае, если нефть мигрирует в нефтяной пласт в качестве отдельных молекул, равновесное состояние смачиваемости горной породы W(A) может быть определено, если заданы насыщенность исследуемого образца нефтью S_O (отношение объема пористого пространства, занятого нефтью к суммарному объему пористого пространства образца), смачиваемость 12 твердой фазы горной породы W⁰(A), фактор шероховатости 10 поверхности F(A), и модель пустотного пространства 6.

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

Смачиваемость твердой фазы горной породы совместно с фактором шероховатости поверхности F(A) определяют эффективную начальную смачиваемость 14 (W_эфф⁰(А)) в точках поверхности А, которая затем в явном виде используются при численном моделировании 15 гидродинамики системы. Процедура 13 определения W_эфф⁰ (А) через W⁰(A) и F(A) может быть осуществлена разными способами, например, согласно (3).

В результате расчетов по данному методу, система приходит в равновесное состояние, соответствующее минимуму полной свободной энергии Гельмгольца. Например, в случае однородного распределения смачиваемости твердой фазы по поверхности W_эфф⁰(А)=W_эфф⁰=const, в зависимости от значения W_эфф⁰ возможны следующие четыре варианта равновесного распределения жидкостей (A.Y. Demianov, О.Y. Dinariev and N.V. Evseev, Introduction to the density functional method in hydrodynamics, Moscow: Fizmatlit, 2014, pp. 57-70, 76-86):

- W_эфф⁰ < -1: Вся поверхность породы покрыта слоем воды, нефть занимает центральную часть пор, не прикасаясь к поверхности породы.

- -1 < W_эфф⁰ < 0: Вода занимает все поры с большим отношением площадь/объем (ПО - отношение площади поверхности данной поры к ее объему). В то время как нефть займет поры с минимальным отношением ПО. Равновесные распределения определятся насыщенностями флюидов.

- 0 < W_эфф⁰ < 1: Наоборот, нефть займет все поры с большим отношением ПО, а вода с меньшим.

- W_эфф⁰ > 1: Вся поверхность породы будет покрыта слоем нефти. Вода займет центральную часть пор, не прикасаясь к поверхности.

Для определения равновесного распределения смачиваемости всей поверхности А далее делается предположение, что после прихода системы в равновесное состояние, те поверхности твердой фазы, которые оказались в равновесном состоянии в контакте с нефтью (А_О), изменили смачиваемость, в то время как смачиваемость оставшейся части поверхности (A_W, покрытой водой) - осталась неизменной. Для задания измененной смачиваемости поверхности А_О могут использоваться различные методы, например, может отдельно исследоваться смачиваемость образовавшихся нефтяных отложений в системе вода-нефть-пленка. В простейшем же случае можно принять W(A_O)=-1. Таким образом, в данном случае равновесная смачиваемость породы будет определяться смачиваемостью твердой фазы W(A_W)=W_эфф⁰(А_W) той части поверхности, которая соприкасалась только с водой (по результатам моделирования равновесного распределения), смачиваемость остальной части поверхности образца принимается W(A_O)=-1.

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

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

Для сложной геометрии поверхности горной породы при рассмотрении баланса сил необходимо также учитывать дополнительное давление Лапласа, возникающее за счет кривизны поверхности пленки вблизи поверхности породы:

где R₁ и R₂ главные радиусы кривизны межфазной поверхности нефть-вода, а R=R₁R₂/(R₁+R₂) - эффективная локальная кривизна поверхности. В дальнейшем, несмотря на то, что локальная кривизна этой поверхности может несколько отличаться от локальной кривизны твердой поверхности, вблизи которой находится эта пленка, мы будем предполагать, что это отличие несущественно, принимая их равными друг другу. С учетом того, что толщина тонких пленок составляет единицы - сотни нанометров, это предположение кажется оправданным при описании данного метода.

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

Для определения локальной эффективной кривизны поверхности в каждой точки поверхности А может использоваться несколько методов. Например, возможно определение прямым способом из дискретного представления поверхности в шаге 2 с последующим расчетом необходимых производных величин координа (V.A. Toponogov, Differential Geometry of Curves and Surfaces. Birkhauser, Boston, 2006, pp. 83-86).

Вычисление кривизны также может быть выполнено по исходному объемному изображению (см., например, В. Rieger, F.J. Timmermans, L.J. van Vliet and P.W. Verbeek, "On Curvature Estimation of ISO Surfaces in 3D Gray-Value Images and the Computation of Shape Descriptors," IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 26, no. 8, pp. 1088-1094, 2004). В этом случае в каждой точке трехмерного изображения рассчитываются структурные градиентные тензоры интенсивности, для которых затем решается задача поиска собственных векторов (СВ) и собственных значений (СЗ). Для вокселей, лежащих на поверхности раздела пора-порода, СВ с большим СЗ соответствует нормали к этой поверхности, в то время как два других - двум главным направлениям кривизн поверхности в данной точке трехмерного пространства. Нормы векторов производных от векторного поля нормалей по этим двум направлениям в каждой точке поверхности пора-порода (поверхность А) определяют пару значений R₁ и R₂ главных кривизн поверхности в рассматриваемой точке.

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

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

После того, как найдена эффективная кривизна поверхности можно рассчитать давление Лапласа (4). Далее, если известны изотермы расклинивающего давления (из литературных данных или измерены) тонкой пленки воды, образованной между поверхностью данного минерала и нефтью, может быть рассчитана равновесная толщина пленки воды (см., например, L. Boinovich and A. Emelyanenko, "The prediction of wettability of curved surfaces on the basis of the isotherms of the disjoining pressure," Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 383, pp. 10-16, 2011). Таким образом можно определить, в каких областях поверхности пленка воды отсутствует (толщина пленки равна нулю) и, следовательно, полярные компоненты нефти могли изменить смачиваемость данной области поверхности. Смачиваемость данной части поверхности меняется на W(A_O)=-1, а смачиваемость оставшейся части поверхности остается неизменной.

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

Для выбора критической величина R' могут быть использованы несколько вариантов. Если известно соотношение площадей поверхности покрытых адсорбированными нефтяными компонентами к общей площади поверхности - коэффициент гидрофобности, то значение R' можно подобрать итеративно, рассчитав площадь поверхности, соответствующей поверхности с кривизной меньше R' и больше R'. Если соотношение площадей этих поверхностей неизвестно, то критическое значение R' может быть подобрано итеративно с помощью следующей процедуры. В качестве начального значения R' может быть взято любое значение, например, R'=0. Для всей поверхности, обладающей меньшей кривизной условия смачиваемости не меняются, а для поверхности обладающей большей кривизной условия смачивания меняются на W=-1. Далее с помощью компьютерного моделирования, например, методом функционала плотности, могут быть рассчитаны параметры, зависящие от смачиваемости, такие как, например, относительные фазовые проницаемости, электрическая проводимость, и т.д. В результате моделирования выбирается параметр R', позволяющий приблизить результаты моделирования к аналогичным экспериментальным результатам наилучшим образом.

После того, как критическое значение кривизны поверхности определено, равновесная смачиваемость поверхности определяется следующим образом: смачиваемость поверхности, обладающей локальной кривизной, меньше критического значения остается неизменной (равной смачиваемости твердой фазы породы W⁰), а смачиваемость остальной части поверхности меняется на W(A_O)=-1.

Основные процедуры описанного способа (получение изображения 2, выполнение сегментации 4, определение распределения минералов 7, учет фактора шероховатости 9 и, безусловно, само моделирование процесса первичной миграции нефти 15) реализуются с использованием электронных вычислительных систем. Эта система может быть единственной для всех упомянутых процедур. Однако, как правило, используются различные вычислительные машины для каждой из задач. Общие требования ко всем этим вычислительным системам состоят в наличии:

- процессора для выполнения необходимых команд и вычислений;

- памяти для хранения временной и постоянной информации (объемом достаточным для работы с получаемыми данными);

- устройств ввода и вывода для управления и контроля выполняемых процедур;

- а также возможности передачи информации через сеть или съемные носители.

В простейшем случае реализации изобретения такой вычислительной системой может быть персональный компьютер.