Бесплатный автореферат и диссертация по наукам о земле на тему

Экспериментальное исследование динамики захвата частиц и изменения проницаемости при фильтрации суспензии через пористую среду

ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

Автореферат диссертации по теме "Экспериментальное исследование динамики захвата частиц и изменения проницаемости при фильтрации суспензии через пористую среду"

На правах рукописи

тг_

Рыжиков Никита Ильич

ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ДИНАМИКИ ЗАХВАТА ЧАСТИЦ И ИЗМЕНЕНИЯ ПРОНИЦАЕМОСТИ ПРИ ФИЛЬТРАЦИИ СУСПЕНЗИИ ЧЕРЕЗ ПОРИСТУЮ СРЕДУ

Специальность 25.00.10 - геофизика, геофизические методы поисков полезных ископаемых

Л"4-

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

005554840

Москва 2014

005554840

Работа выполнена на кафедре «теоретическая и экспериментальная физика геосистем» Московского Физико-Технического Института (Государственного университета)

Научный руководитель: кандидат физико-математических наук

Михайлов Дмитрий Николаевич

Официальные оппоненты: Доктор физико-математических наук

Гарагаш Игорь Александрович

Доктор физико-математических наук Филиппов Анатолий Николаевич

Ведущая организация: Башкирский государственный университет

Защита диссертации состоится «16» октября 2014 г. в 13 часов на заседании диссертационного совета Д 002.050.01 при Федеральном государственном бюджетном учреждении науки Институте динамики геосфер Российской академии наук по адресу: 119334, г. Москва, Ленинский проспект, д. 38, корп. 1.

С диссертацией можно ознакомиться в библиотеке и на сайте Федерального государственного бюджетного учреждения науки Института динамики геосфер Российской академии наук, http://idgl.chph.ras.ru/

Автореферат разослан «/у »(Зр+ОЧ^УЫ 2014 г.

Ученый секретарь Диссертационного совета,

кандидат физико-математических наук В.А. Рыбаков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность

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

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

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

На сегодняшний день получено несколько разновидностей феноменологических соотношений для описания течения различного рода дисперсных систем и полимеров. Тем не менее, вопрос о диапазоне применимости полученных зависимостей и их однозначности остается открытым.

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

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

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

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

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

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

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

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

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

Традиционные лабораторные фильтрационные эксперименты по исследованию воздействия суспензии на образец керна позволяют определить только изменение

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

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

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

Основные задачи работы:

• создание методики проведения фильтрационных экспериментов;

• проведение фильтрационных экспериментов с различными составами суспензий и различными образцами горных пород;

• создание экспериментальных методик для количественной оценки параметров зоны кольматации в образцы горных пород;

• интерпретация экспериментальных данных с использованием полученных количественных параметров зоны кольматации в образцах;

Научная новизна.

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

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

3. Метод окрашивания компонент суспензии с последующим анализом сколов образца для получения распределения проникших частиц суспензии.

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

5. Метод количественной интерпретация эффекта изменения скорости распространения акустических волн в образце песчаника из-за проникновения твердых компонент суспензии.

Защищаемые положения.

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

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

3. Метод окрашивания компонент суспензии с последующим анализом сколов образца для получения распределения проникших частиц суспензии.

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

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

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

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

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

Результаты данной работы легли в основу 3-х патентов и 3-х патентных заявок.

Личный вклад

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

Апробация работы.

По теме диссертации автором опубликовано 9 научных работ, в том числе 3 статьи в рецензируемых журналах (из них 3 в журналах из перечня ВАК), 6 статей в

7

научных сборниках и трудах конференций. Результаты исследований докладывались на семинарах в ИДГ РАН и в Московском Научном Центре Шлюмберже. На российских научных конференциях в Московском Физико-Техническом Институте (2010, 2011, 2013), в РГУ им. Губкина (2013, 2014), в БашГУ(2014), КГТУ (2012). Всего автор принял участие в 7 научных конференциях, из них 2 международных.

Структура и объем работы.

Диссертация состоит из введения, шести глав, заключения и списка использованной литературы из 139 наименований. Общий объем диссертации составляет 150 страниц, на которых размещено 84 рисунков и 4 таблицы.

Благодарности.

Автор выражает глубокую признательность своему научному руководителю к.ф.-м.н. Михайлову Д.Н., коллективу московского научно-исследовательского центра Шлюмберже (Шако В.В., Бурухину А.А., Жарниковой А.В., Надееву А.А., Якимчуку И.В.) за помощь и поддержку в ходе работ по теме диссертации. Автор также признателен Турунтаеву С.Б. и коллективу Института Динамики Геосфер за экспертную оценку работу и ценные советы.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

Первая глава посвящена обзору существующих теоретических и экспериментальных подходов к исследованию течения суспензии через пористую среду, механизмов захвата частиц, изменения проницаемости и структуры порового пространства. Приведен обзор литературы по теоретическому описанию исследуемых процессов, включая классическую ("deep-bed filtration") модель, стохастические и сеточные модели.

Подробно рассмотрена классическая модель (deep-bed), в рамках которой течение суспензии в пористой среде описывается системой уравнений, состоящей из закона Дарси, уравнения транспорта твердых компонент, уравнения кинетики захвата

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

Уравнение определяющее интенсивность захвата частиц кольматанта в поровом пространстве представляется в виде :

— =Ä0F(<T)WC (1)

dt

Здесь, w - фильтрационная скорость, С - концентрация частиц в суспензии, Х0 -коэффициент захвата; F(a) - корректирующая функция, учитывающая зависимость интенсивности захвата частиц от доли объема пористой среды, занимаемой захваченными компонентами.

Для связи между проницаемостью и объемной долей захваченных частиц наиболее распространены эмпирические зависимости вида (2)(см. Herzig [Herzig, 1970]).

к _ 1

к0~ l+ßo (2)

А также вида (3), основанная на обобщении уравнения Козени-Кармана [Михайлов, 1987; Шехтман, 1961; Воек, 2012; Herzig, 1970; Tien, 1979].

к ( а\м

— = (1-У—) (3)

к0 \ т0/

где /?, у, M — эмпирические параметры.

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

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

Во второй главе изложены условия и процедура проведения фильтрационных экспериментов. Приведены параметры исследуемых образцов горных пород и суспензий. Схема проведения фильтрационных экспериментов представлена на Рис. 1.. Описана используемая в работе установка ПИК-ИБР, произведенная компанией ООО «Геологика». Последовательность процедур во ходе фильтрационного эксперимента представлена на Рис. 2.

Схема эксперимента

1 ^^У

■ ...^..Х.----^

|........| Фильграг

□Гидравлическое масло

I 1 Суспензия

—* Направление

жидкости в режиме циркуляции бурового раствора

—* Направление

обратной прокачки

Рис. 1. Принципиальная схема экспериментов по прямой и обратной прокачке жидкостей через образцы кернов. Буквами обозначены различные механизмы повреждения горной породы: А - внешняя корка, Б — внутренняя корка, В - проникновение фильтрата.

Измерение проницаемости А

Фильтрация суспензии А

Обратная прокачка Л_

о 500 1000 1soo 2000 2500

Время, с

Рис. 2. Динамика перепада давления и расхода в ходе фильтрационного

эксперимента.

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

Томография выполнялась на приборе SkyScan 1172. Для идентификации проникших в пористую среду частиц требуется высокое разрешение 2-3 мкм на воксель (воксель — минимальный элемент трехмерного изображения кубической формы), что достигается малым диаметром сканируемой области (менее 10 мм). Сканирование области внутри образца стандартного керна приводит к высокому уровню сигнал/шум. Поэтому для сканирования с требуемым разрешением из стандартного образца керна (диаметр - 3 см, длина — 6 см) выпиливались образцы диаметром 8 мм и длиной от 1 до 1.5 см.

Для расчета профилей концентрации захваченных частиц разработан метод анализа гистограмм [Рыжиков, 2013]. Гистограмма сечения томографии - это график статистического распределения элементов цифрового изображения с различной яркостью, в котором по горизонтальной оси представлена яркость, а по вертикали — относительное число пикселов с конкретным значением яркости.

Гистограмма градаций яркости, нормализованная на общее количество

составляющих ее пикселей, является функцией плотности вероятности значения

серого цвета на выделенном сечении. Согласно [Gonzalez R. С., 1992], такая функция

11

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

где г — значение яркости, р(г) - аппроксимация функции плотности вероятности от всего сечения изображения, р,(г) - функция плотности вероятности для /-той подобласти, п — общее число подобластей, на которое разбивается изображение.

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

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

Таким образом, для характеризации функции плотности вероятности отдельного материала на сечении требуется три параметра: интенсивность (Л,), полуширина (а,) и среднее значения

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

п

(4)

¡=1

О-и)2 р,00 = ¿¡е 2ст1

(5)

о

X ь

к о о.

си

СП А I-и О т ь о с с

0.018 0.016 0.014 0.012 0.01 0.008 0.006 0.004 0.002 0

Яркость

- ФПВ сечения

- Скелет породы

- Поровый флюид

- Загрязнитель

- Переходная зона (Наполнитель - Скелет)

- Переходная зона (Скелет - Загрязнитель)

Рис. 3. Пример гистограммы сечения рентгеновской компьютерной микротомографии

Области на границе двух материалов называются «переходной зоной», пиксели которой в литературе называются микселями (пихек). Эти области (Рис. 4) имеют собственную функцию плотности вероятности рм, параметры которой зависят от параметров функций плотности вероятности граничащих материалов (¡'-го и к-го) и характеризуются только одним параметром - интенсивностью, которая пропорциональна протяжённости границы между материалами.

Рис. 4. Переходная зона (А) и профиль изменения серости в переходной

зоны (Б)

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

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

О (6)

У=1 ¡=1 ¿-1,к=1

Число свободных параметров в такой минимизации равно Зп+п(п-1)/2, где п —

количество материалов представленных на сечении. Относительные объемы

материалов в сечении ^ определяются выражением (7). где - площадь, занимаемая

чистым материалом с индексом /, 8р/ц - площадь переходной зоны между к-м и г'-м

материалами.

тг-1

1 V С (7)

кфЬ

Для расчета профиля объемной доли различных компонент вдоль анализируемого образца описанная выше процедура последовательно применяется ко всем сечениям рентгеновской компьютерной микротомографии, причем результаты минимизации для /-ого сечения используются как начальное приближение при минимизации для (1+1)-ого сечения. Кроме того, на результат минимизации для (1+1)-ото сечения накладывается ограничение в виде наперед заданного отклонения е от результатов /-ого, т.е. для любого параметра минимизации А должно выполняться условие \1-А1+1/А,\<£, где / - порядковый номер сечения, на котором А/ является решением обратной задачи. В практических расчётах £ варьировалась в диапазоне от 3% до 10 %.

Применение предложенного метода для определения параметров захвата частиц описано в шестой главе.

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

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

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

Рис. 5. Фотографии сколов образцов с различной проницаемостью через, которые фильтровалась суспензия бентонита. Двумя латинскими буквами зашифровано название породы из которой были сделаны образцы: IL -Indiana Limestone, GB - Grey Berea, BB - Buff Berea, CG - Castlegatc, BT-

Bentheimer

В среде MATLAB была создана программа анализа изменения цвета фотографии образцов по трем каналам: красный (R), зеленый (G) и синий (В). Фотография представляется в виде трех двумерных матриц размером MxN, соответствующих распределению интенсивности цветов формата RGB: R¡j, G¡j. Ви, где i — индекс элемента матрицы по направлению вдоль оси фильтрации, a j - поперек. В качестве характеристики количества проникших частиц окрашенного бентонита S¡j использовалась сумма абсолютных значений изменения цвета относительно фоновых значений (8):

% = К - Ко| + |Су - С0| + \Вц - В0|, <8)

где фоновые значения R0, G0, В0 полученные на чистом образце.

Для построения профиля распределения частиц S¡ вдоль направления фильтрации матрица S,; суммируется поу'-ому индексу.

Анализ корреляции между а(х) и S(x) выходит за рамки работы, поскольку для определения параметров захвата частиц суспензии необходима только форма

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

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

Для излучения и приема акустических волн использовались датчики РапатеЦтсБ V ЮЗ-ЯМ, апертура датчика равнялась 1.3 см, а основная частота - 1 МГц (Рис. 6). С помощью системы позиционирования (см. Рис. 7) осуществлялось перемещение датчиков и их диаметральный прижим к цилиндрическому образцу.

На Рис. 8 синим маркером представлен пример результатов профилирования для образца ВТ11 (песчаник ВеЩЬе1шег. пористость — 0.23, проницаемость 2900 мД), через который фильтровалась взвесь частиц карбида кремния в полимере (25 г/л частиц БЮ и 1 г/л полимер ХапЙап). Как видно, во входной части образца наблюдается повышение скорости ультразвуковых волн, которое однозначно можно связать с проникновением частиц карбида кремния в пористую среду.

Рис. 6. Датчики РапатеЫса VЮЗ-ЯМ и латунные держатели (акустические линзы).

Рис. 7. Система позиционирования датчиков

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

Для описания изменения скорости распространения акустических волн использовалась система уравнений теории Френкеля-Био-Николаевского для распространения упругих волн в пористой среде в одномерном случае в линейном приближении (9)-( 12).

сН>2 д \ ПП

ах к

а^=К„е1+гр,; кр=кьЛо\ е = рР = (1-^^, + тор2> (12)

где индексы 1, 2 - твердая фаза и насыщающая жидкость (газ) соответственно; р, - плотности фаз; V, - массовые скорости; е\ - продольная деформация матрицы; аец - эффективное напряжение Терцаги; р - давление в насыщающем флюиде; т0 -пористость; £ — проницаемость; ц - вязкость флюида; Кь - модуль объемной упругости пористой матрицы; в - модуль сдвига пористой матрицы; Р, — сжимаемости материала, составляющего пористую матрицу, и насыщающего флюида.

Для данной системы искалось решение в виде гармонического сигнала с

частотой со и волновым числом которое в результате подстановки давало

17

дисперсионное соотношение (13), являющееся биквадратным уравнением относительно Ç.

Два решения этого уравнения характеризуют две продольные волны распространяющихся в пористой среде. Интерес для нас представляет только решение соответствующее продольной волне с низким затуханием и высокой скоростью.

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

Параметр Е0 определялся из скорости распространения акустических волн в чистом образце, а параметр а - с помощью данных рентгеновской компьютерной микротомографии, согласно которой, на глубине 5 мм, объемная концентрация частиц а составляет 5%. Полученная зависимость была использована для восстановления профиля концентрации захваченных частиц вдоль образца, представленного на Рис. 8 (зеленый маркер).

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

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

В случае классического описания процесса кольматации с постоянной интенсивностью захвата F(cг) = 1 существует аналитическое решение для объемной доли захваченных частиц (16) и для концентрации взвешенных в суспензии частиц

(13)

771 = 771g — с

(14)

Е = Е0 + acr

(15)

" (0 _ ХТП ^ ^ 0,711 (_Л ог- пч \Л71 >1 г СN гш

0,

если < хт0

С0 ехр(-Л0;с), если > хта

С(х) =

(17)

0,

если < хт0

где х - координата вдоль оси керна; \\'1п] - объем закаченной в образец жидкости, отнесенный к площади сечения образца; Х.0 - коэффициент захвата.

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

Для анализа данных фильтрационных экспериментов удобно использовать гидравлическое (фильтрационное) сопротивление 3, соответствующее отношению текущего перепада давления Ар(1) к текущему расходу (¡(I). На Рис. 9 представлены гидродинамические данные для эксперимента по закачки суспензии взвешенных частиц карбида кремния через песчаник ВепЛетег с последующей обратной прокачкой. Предполагая, что при начале обратной прокачки внешняя фильтрационная корка разрушается практически мгновенно, был соотнесен начальный уровень гидравлического сопротивления при обратной прокачке (обозначен штриховой линией) с уровнем гидравлического сопротивления образца при наличии только внутренней фильтрационной корки. Таким образом в фильтрационных экспериментах определялся объем суспензии ()*, потраченный на формирование внутренней фильтрационной корки. В данном случае объем составил <2* = 21-23 мл.

0.08

£

Ч 0.03

0.01

0.02

0

О

100

Объем закачки, мл

Рис. 9. Динамика гидравлического сопротивления в ходе фильтрационного эксперимента. Образец - песчаник ВепШетег, загрязнитель - взвесь

частиц в воде. Цифрами на рисунке обозначены области: 1 -формирование внутренней фильтрационной корки, 2 - формирование внешней фильтрационной корки, 3 - обратная прокачка; <2*- объем суспензии, потраченный на формирование внутренней фильтрационной

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

где Ь - длина образца; ц - вязкость суспензии.

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

Для точного определения параметра Хо будем использовать профиль полученный с помощью анализа данных рентгеновской компьютерной микротомографии (Рис. 10).

корки

(18)

I 14

S

9 12

сч

—-■■-... Данные томографии — • - Аппроксимация 1 (Ь = 0, Хо = 1150 м"л) ---Аппроксимация 2 (Ь = 850, Хо = 120 м ')

\

\

0.5 1 1.5 2 2.5 В 3.5 Расстояние от входного торца, мм

Рис. 10. Профиль захваченных частиц карбида кремния в образце песчаника Bentheimer и аппроксимации профиля с помощью выражений

(16) и (20).

Проведенный анализ показал, что данный профиль не может быть аппроксимирован с достаточной точностью с помощью решения (16) при объеме прокачки Q* = 21-23 мл и заданной концентрации суспензии. Решение (16) позволяет аппроксимировать только часть профиля при Q* = 17 мл и не воспроизводит резкое изменение доли захваченных частиц вблизи входного торца образца.

Для адекватной аппроксимации всего профиля рассмотрим зависимость интенсивности захвата частиц от доли захваченных частиц в объеме пористой среды в виде (19) [Herzig, 1970;Boek, 2012]. да

— = A0Cw (1 + b&) (19)

öt

При таком выражении для захвата частиц существует аналитическое решение [Шехтман, 1961; Civan, 2007] (20):

!1 ехрС-ЛоСпЬ^тДО-этпо]) . .

bl-eXpCA0i)-exp(-A0C„b[win;(t)-xm0])' inj W - ХТП, 0 ^

0, Winj(t) < xm0

Аппроксимация профиля с помощью решения (27) методом наименьших квадратов приведена на Рис. 10 штрих-пунктирной линией. Полученные параметры: b = 850, Х0 — 120 мл.

Воспроизведен профиль захваченных частиц и тем самым зафиксированы параметры для уравнения захвата (19). После обоснования модели захвата частиц

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

Объем закачки, мл — Эксперимент

--Ро = 1000, \о = 125 м~\ b = 840, Аппроксимация (1)

— • - М = 20, \о = 125 м"', b = 840, Аппроксимация (2)

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

В случае анализа данных по фильтрации бентонитовой глины использовались профили распределения захваченных частиц, полученные с помощью анализа распределения интенсивности цветов на фотографий сколов образцов. На Рис. 12 представлены результаты анализа изменения цвета на сколе образца песчаника Buff Berea после закачки суспензии бентонитовой глины 10 г/л и аппроксимация полученного профиля с помощью аналитического решения (16). Коэффициент захвата для данной задачи составил \,= 152 м"1.

12 3 4

Расстояние от входного торца, см

- Результаты анализа

---- Аппроксимация

Рис. 12. Результат построения профиля цветовой характеристики проникновения частиц окрашенной бентонитовой глины в образец песчанника Buff Вегеа.

По анализу динамики гидравлического сопротивления образца, представленной на Рис. 13. был определен параметр снижения проницаемости (3 = 1850.

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

30

25

а

S 20

ч

S 15

t

-

UJ 10

я

5

0

- Эксперимент

Внутренняя корка; [5 = 1850 -Внешняя корка; kmc - 0.04 мД

10 15 20

Объем закачкн, мл

30

Рис. 13. Динамика гидравлического сопротивления образца песчаника Buff Вегеа в ходе закачки взвеси бентонитовой глины 10 г/л.

В Заключении изложены основные выводы и результаты.

По теме диссертации опубликованы следующие работы:

Рыжиков Н.И., Михайлов Д.Н., Шако В.В. Метод расчета профилей распределения пористости и объемных долей материалов в пористой среде с помощью анализа данных рентгеновской микротомографии //труды МФТИ. — 2013. -том 5.-№4(20).

Михайлов Д.Н., Рыжиков Н.И., Шако В.В., Комплексный экспериментальный подход к определению параметров зоны кольматации продуктивных пластов» //Вестник ЦКР РОСНЕДРА. -2014. -№1.

Михайлов Д.Н., Рыжиков Н.И., Шако В.В., «Комплексный экспериментальный подход к определению параметров проникновения и захвата компонентов бурового раствора и сопутствующего изменения проницаемости породы коллектора //труды РГУ им. Губкина. - 2013. - № 4.

Михайлов Д.Н., Рыжиков Н.И., Шако В.В. Комплексный экспериментальный подход к определению кинетики кольматации пористых сред. //Нефтяное Хозяйство, (принято в печать)

Михайлов Д.Н., Рыжиков Н.И., Шако В.В., Theuveny В. Комплексный экспериментальный подход к определению изменения свойств пород пласта под воздействием бурового раствора. -SPE-171291-RU-2014. (принято в печать)

Михайлов Д.Н., Рыжиков Н.И., Шако В.В. Комплексный экспериментальный подход для определения влияния проникших компонентов бурового раствора на изменения свойств призабойной зоны, //тезисы конференции «Геофизические исследования в процессе добычи и освоения скважин. Интеллектуальные скважины», 2014.

Рыжиков Н.И. Экспериментальное исследование структуры зоны кольматации и кинетики ее формирования, //тезисы 68-ой международной молодежной конференции «Нефть и Газ 2014». — 2014.

Рыжиков Н.И., Михайлов Д.Н. Экспериментальное исследование динамики захвата частиц и сопутствующего ухудшения проницаемости при фильтрации суспензии через пористую среду //тезисы XX конференции «Губкинские чтения». -2013.

Рыжиков H.H., Михайлов Д.Н. Комплекс экспериментальных методов для определения параметров захвата частиц и снижения проницаемости при фильтрации суспензии через пористую среду //тезисы 56-й научной конференции МФТИ. — 2013.

Рыжиков H.H., Михайлов Д.Н. Метод расчета профилей распределения пористости и объемных долей материалов в пористой среде с помощью анализа данных рентгеновской микротомографии //тезисы 54-й научной конференции МФТИ. -2011.

Рыжиков H.H., Михайлов Д.Н. Расчет профилей твердых компонент пористого образца с помощью рентгеновской томографии //тезисы 53-й научной конференции МФТИ.-2010.

Патенты и патентные заявки:

Патент № W02012165991 Al «A method for determination of spatial distribution and concentration of contrast components in a porous or / and heterogeneous sample», выдан 6.12.2012

Патент РФ № 2467315 «Способ определения пространственного распределения и концентрации глины в образце керна», выдан 23.06.2012

Патент РФ № 2013135670 «Способ прогнозирования изменения свойств призабойной зоны пласта под воздействием бурового раствора», выдан 30.05.2014

Заявка РФ № 2013156000 «Способ определения изменений параметров пористой среды под действием загрязнителя», подана 18.12.2012

Заявка РФ № 2013157418 «Способ определения изменения свойств околоскважинной зоны пласта-коллектора под воздействием бурового раствора», подана 25.12.2013

Заявка РФ № 2013157415 «Способ определения распределения и профиля загрязнителя в пористой среде», подана 25.12.2013

Подписано в печать:

11.09.2014

Заказ № 10198 Тираж - 100 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru