авторефераты диссертаций БЕСПЛАТНАЯ РОССИЙСКАЯ БИБЛИОТЕКА - WWW.DISLIB.RU

АВТОРЕФЕРАТЫ, ДИССЕРТАЦИИ, МОНОГРАФИИ, НАУЧНЫЕ СТАТЬИ, КНИГИ

 
<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ


Pages:     | 1 || 3 |

Численное моделирование термогидродинамических процессов в подземной гидросфере

-- [ Страница 2 ] --

Здесь h – толщина пласта, -постоянна Эйлера, rc - приведенный радиус скважины.

В третьем разделе рассмотрены вопросы тепломассопереноса в пористых средах и других породах литосферы. Здесь актуально исследование температурных изменений в породах и давления фильтрующихся жидкостей. Соответствующая система уравнений для пласта записывается в виде [3]

 (4) Здесь T - температура, m– пористость, w- скорость фильтрации, - коэффициент-9 (4)

Здесь T - температура, m– пористость, w- скорость фильтрации, - коэффициент Джоуля-Томсона, - коэффициент адиабатического сжатия (расширения) флюида, ср, сп- коэффициенты теплоемкости флюида и пласта, п- коэффициент теплопроводности пласта.

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

(5)

здесь Т0 – температура в скважине в момент прекращения теплового воздействия (остановка скважины); Тп – температура пород на данном срезе до начала возмущения; Т – текущая температура в момент времени t после остановки; t0 - длительность теплового возмущения (закачки или добычи), - коэффициент температуропроводности пород.

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

(6)

Здесь – E и Еп - модули Юнга горных пород и пласта, H- глубина залегания пласта, - коэффициент Пуассона пород, - вспомогательная функция, связанная с распределением давления в пласте, определяемая условиями конкретной задачи.

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

(7)

здесь Т- температура, с, - коэффициенты теплоемкости и теплопроводности, V- скорость вертикального движения, z -вертикальная координата, H- нижняя граница литосферы.

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

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

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

 Кривая восстановления давления скважины Архангельского месторождения. -20

Рис. 1. Кривая восстановления давления скважины Архангельского месторождения.

Рис.2. Кривая восстановления давления (сплошная) и динамическая гидропроводность (пунктирная) в пласте с ухудшенной призабойной зоной (2=1000 с, ’=300 с).

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

В первом разделе рассматриваются решения уравнения (2) для случая 10, 20, p=0, w=0 (трещиновато-пористая среда без учета неравновесного характера фильтрации). Кривые восстановления давления, полученные из решения системы (1) имеют два линейных участка на больших и малых временах. С уменьшением ’ (увеличением разницы между упругоемкостями трещинного и блокового пространств), величина первого линейного участка сокращается. Для реальных кривых восстановления давления, на временах наблюдения от десятков секунд, характерно отсутствие первого линейного участка, что, возможно, свидетельствует о малой упругоемкости трещин. Первоначальные участки реальных кривых представляют собой период относительно слабого изменения давления и имеют длительность 103-104 секунд (рис.1). Увеличение продолжительности первоначального участка наблюдается при росте приведенного радиуса скважины и уменьшении пьезопроводности пласта, увеличении гидропроводности призабойной зоны и наличии послепритока q(t) из пласта. В случае наличия притока q(t) обработка кривых восстановления давления весьма удовлетворительна в координатах [lg t, (p-p0)/(Q0-q(t))].

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

Во втором разделе рассматриваются решения (2) для случая 1=0, 2=0, p0, w0 (пористая среда с релаксационными эффектами). Выявлено, что при представлении кривой восстановления давления, состоящей условно из трех участков, значение феноменологического параметра p может быть оценено как время, соответствующее середине второго участка. Релаксационный эффект, определяемый параметром p, отражается в смещении линии КВД в сторону больших времен и запаздывании проявления неоднородностей на графике гидропроводности. Он также не позволяет определять значения гидродинамических параметров ближней к скважине области с достаточной точностью. В разделе подробно рассматриваются и анализируются различные ситуации, реализуемые на практике, связанные с наличием послепритока, неоднородностью пласта, зависимостью проницаемости от давления.

В третьем разделе рассматривается общий случай: 10, 20, p0, w0 : трещиновато-пористая среда с релаксационными свойствами. Показано, что учет неравновесности фильтрации приближает вид модельных кривых восстановления давления в трещиновато-пористых пластах к реальным (рис.3). При этом величина p сказывается на малых, а соотношение упругоемкостей трещин и блоков (’,2) – на больших временах наблюдения. Для исследованных ситуаций из сопоставления с экспериментальными кривыми восстановления давления порядок значений временных параметров p,2 может быть оценен как 103 -104 с.

Рис. 3. Кривая восстановления давления скважины Архангельского месторождения (точки - экспериментальные данные, линия – модельная кривая).

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

В третьей главе рассмотрены вопросы массообмена флюидов в системе блоки-трещины.

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

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

, (8)

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

В разделе 3 рассматривается обмен жидкостью между блоками и трещинами при периодическом зондировании пласта.

 Безразмерная величина массообмена между блоками и трещинами в зависимости от-23

Рис.4. Безразмерная величина массообмена между блоками и трещинами в зависимости от периода циклического воздействия; 1- 2=1800 с, 2- 2=3600 с, 3- 2=7200 с, 4- 2=14400 с.

Для анализа массобмена между блоками и трещинами при периодическом гидродинамическом воздействии на пласт [1,2] можно использовать соотношение

. (9)

Здесь

, Pc – амплитуда колебаний давления на скважине, - сдвиг фаз между дебитом и давлением на ней.

Как показывают проведенные расчеты, для каждой точки в пласте массообмен между блоками и трещинами в зависимости от периода воздействия имеет максимум. Значение периода, соответствующего максимуму, его можно называть оптимальным периодом, находится в прямой зависимости, прежде всего, от параметра 2. Причем с увеличением 2 увеличивается и ширина интервала периодов, обеспечивающих максимальный обмен (рис.4). С удалением от скважины длительность оптимальных периодов возрастает.

В четвертом разделе приводятся результаты использования амплитудно-частотных характеристик (АЧХ) - и – фазо-частотных характеристик - (ФЧХ) системы для решения задачи идентификации линейной фильтрационной модели, а также методика использования диаграмм Найквиста для анализа исследуемого периодического процесса (здесь AQ, Ap– амплитуды, Q, p- фазы сигналов по расходу и давлению в скважине). Действительно, в рамках концепции линейного отклика передаточную функцию любой сложной системы можно представить в виде произведения передаточных функций простых подсистем, если их соединить последовательно. При задании периодических колебаний давления на скважине комплексная передаточная функция для (2) будет выглядеть для установившихся режимов как

(10) Пример амплитудно-частотной и фазо-частотной характеристик пласта (точки --28(10)

 Пример амплитудно-частотной и фазо-частотной характеристик пласта (точки --30 Пример амплитудно-частотной и фазо-частотной характеристик пласта (точки --31

Рис. 5. Пример амплитудно-частотной и фазо-частотной характеристик пласта (точки - экспериментальные данные, линии - расчетные).

Модельные АЧХ и ФЧХ аналогичны характеристикам карбонатных коллекторов, полученным экспериментально [2] (рис.5), что свидетельствует о правомерности описания трещиновато-пористых коллекторов используемой моделью. Частотные характеристики могут быть использованы для оценки временных и гидродинамических параметров. Приведенный подход является не только радиофизическим взглядом на рассматриваемую проблему, но и используется для подбора адекватных значений временных параметров для реальных систем.

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

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

В разделе 1 рассматривается распределение температуры на квазистационарном режиме, вводится понятие эффективного коэффициента Джоуля-Томсона (ЭКДТ) как проявление связи гидродинамических и теплофизических процессов в пористых средах [3], где Tсм-температура смеси газа из пласта, TГ - геотермическая температура на уровне кровли пласта, pк - давление на контуре питания пласта, pз - давление на забое, возможности его определения и практического применения для определения давления в пласте, оценки проницаемости пласта.

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

 Динамическая гидропроводность: пунктирная определенная из изменения давления,-33

Рис. 6. Динамическая гидропроводность: пунктирная определенная из изменения давления, сплошная - из изменения температуры (задано скачкообразное изменение проницаемости пласта на r*=5 м).

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

Вторжение в коллектор флюида с температурой, отличной от окружающих пород, искажает невозмущенное температурное поле. Однако по скважине, в которой можно было бы произвести измерение температуры, идет поток, вносящий дополнительные возмущения в температурное поле. Для регистрации температуры используются кратковременные остановки скважины. В реальных условиях для определения t0 (см. (5)) используется результаты восстановления температуры по слою пород, заведомо не содержащему коллекторы (соли, глины), невозмущенная температура которого известна. Для реализации предлагаемого метода необходимо зарегистрировать не менее 3 термограмм по глубине на нестационарных режимах. На основе моделирования процесса восстановления температуры определяется время записи каждой из термограмм с целью обеспечения наибольшей точности расчета температуры пород при условии сокращения времени остановки скважины. Отклонение расчетной температуры от невозмущенной геотермы сверх указанной точности определяет места поглощения. Метод был реализован на Совхозном ПХГ, нефтяных месторождениях Татарии и Пермской области. На способ получен патент РФ № 2013533. Приводятся примеры реализации метода (рис.7).

 Расчетная температура пород (линия с точками) и невозмущенная геотерма-34

Рис.7. Расчетная температура пород (линия с точками) и невозмущенная геотерма (сплошная) скважины Ромашкинского месторождения (интервал скопления в районе 125 м).

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

 Термограммы после остановки скважины, цифры соответствуют времени после-35

Рис.8. Термограммы после остановки скважины, цифры соответствуют времени после остановки скважины в часах (поглощающий интервал 57.5-60 м).

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

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



Pages:     | 1 || 3 |
 





 
© 2013 www.dislib.ru - «Авторефераты диссертаций - бесплатно»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.