CC BY
60
УДК 622.244 ББК И Р. А. Гасумов, В. А. Толпаев, О. С. Кондренко
г. Ставрополь, Россия
Вывод системы уравнений для приближенного расчета температуры бурового раствора при проводке скважины в многолетнемерзлой породе
Предложен приближенный метод расчета количества тепла, которое буровой раствор отдает в многолетнемерзлую породу, окружающую скважину. Выведены дифференциальные уравнения теплового баланса для скважины с теплоизолированной стенкой. Для расчета фактического распределения температуры применяется поправочный коэффициент. Результаты статьи проиллюстрированы тестовым расчетом, представленным в виде графика.
Ключевые слова: скважина, многолетнемерзлая порода, теплопроводность.
R. A. Gasumov, V. A. Tolpaev, O. S. Kondrenko
Stavropol, Russia
Derivation of Equations Set for the Approximate Calculation of the Drilling Fluid Temperature in the Permafrost Well
An approximate method of calculating the amount of heat that drill mud delivers into the permafrost, surrounding the hole. The differential equations of thermal balance for a well with a thermally insulated wall are derived. To calculate the actual temperature distribution a correction factor is used. Our results are illustrated by the test calculations presented in a kind of graphic.
Keywords: the well, permafrost, thermal conductivity.
Проектирование скважин, предназначенных для разработки газовых месторождений в районах Крайнего Севера, осуществляется с обязательным учетом оценки радиусов протаивания многолетнемерзлых пород (ММП) вокруг строящейся скважины. Расчет радиусов протаивания ММП необходим для правильного выбора бурового раствора при бурении скважины под кондуктор, для нахождения его температуры на входе, а также для определения оптимальной скорости проводки.
Расчетная схема теплового взаимодействия нисходящего и восходящего потоков бурового раствора приведена на рис. 1:
О
а3 $ X
2 х + Ах
Ф)
Д Л3 $
/ X
Рис. 1.
где 1 — сегмент буровой колонны, 2 — сегмент ствола скважины,
3 — нисходящий поток бурового раствора, 4 — восходящий поток бурового раствора, До — внешний диаметр бурильной колонны, Дс — диаметр скважины
60 © Гасумов Р. А., Толпаев В. А., Кондренко О. С., 2011
В процессе бурения скважины по буровой колонне (1) вниз к долоту направляется буровой раствор (3), который затем выходит из отверстий долота и направляется обратно вверх (4), как схематически показано на рис. 1. Восходящий поток бурового раствора выносит на дневную поверхность шлам разбуренной породы и воду, получающуюся при таянии льда ММП. На рис. 1 температуры нисходящего (3) и восходящего (4) потоков бурового раствора обозначены соответственно как То(х) и Ті(х) (в градусах Кельвина). Через х и I (в м) соответственно обозначены декартовые координаты текущей точки наблюдения и местоположения долота. Начало отсчета на направленной вертикально вниз координатной оси идет от дневной поверхности.
Буровой раствор, имея на входе начальную температуру Твх и двигаясь вниз по бурильной колонне, доходит до долота и получает от него дополнительное количество тепла, выделяемого долотом при работе. Поэтому температура То(х) в нисходящем потоке бурового раствора растет по пути к долоту, как это показано в виде графика, выходящего из точки (1,0) на рис. 2. Скорость роста температуры нисходящего потока различна на разных участках. Точка излома графика зеленого цвета соответствует началу участка интенсивного тепловыделения, исходящего от работающего долота.
Схема распределения температур в нисходящем и восходящем потоках бурового раствора для случаев, когда стенка скважины теплоизолирована от окружающей породы и когда есть тепловой контакт с окружающей ММП, приведена на рис. 2:
-----------------------1----------------------1----------------------1----------------------1----------------------1----------------------1-----------------------1----------------------1----------------------1----------------------1
О 0,2 0,4 0,6 0,8 1
х
I
Рис.2.
Здесь линии, выходящие из точек (0), (1,0) и (1,1), представляют распределение температур соответственно в восходящем потоке с учетом теплообмена, в нисходящем потоке и в восходящем потоке без учета теплообмена с окружающей породой.
Дополнительно нагретый работающим долотом буровой раствор выходит из отверстий долота и далее движется вверх, вынося шлам от разбуренной породы и подогревая его до своей температуры, а также растапливая лед в разбуренной ММП. Поэтому при обратном движении вверх температура бурового раствора в восходящем потоке падает, вследствие теплопотерь на подогрев шлама ММП и растопление льда. Кроме того, температура восходящего потока по пути вверх падает вследствие теплообмена через стенку бурильной трубы с холодным нисходящим потоком. Если не учитывать теплового взаимодействия восходящего потока бурового раствора с окружающей ММП (например, в предположении, что стенка разбуренной скважины теплоизолирована), то падение температуры Т\(х) по пути вверх в восходящем потоке на рис. 2 будет определяться верхним графиком. Без учета теплообмена с окружающей ММП на выходе буровой раствор, как показывают расчеты, будет иметь температуру ТВых более высокую, чем на входе Твх. Поэтому график температуры Тх{х) на рис. 2 располагается выше графика температуры То(х). Однако в действительности стенка разбуренной скважины находится в тепловом контакте с окружающей ММП, и фактическая температура Твых бурового раствора на выходе бывает ниже температуры на входе Твх. Поэтому график (красного цвета) фактической температуры Тц(х) (из точки (0)) восходящего потока бурового раствора в точке выхода оказывается ниже графика (зеленого цвета) температуры То(х) (из точки (1,0)).
Разница температур Т\(х) и Тц, представленных на рис. 2 соответственно из точек (1,1) и (0), обусловлена оттоком тепла из двигающегося вверх бурового раствора в окружающую скважину ММП. Поэтому для расчета радиусов протаивания ММП вокруг скважины нужно знать функци-
ональные зависимости температуры То(х) в нисходящем потоке, а также теоретической Ті(х) (без учета теплообмена с ММП) и фактической Тц(х) (учитывающей теплообмен с ММП) температур в восходящем потоке.
Целью статьи является вывод уравнений для расчета названных температур. Перечислим исходные данные, которые потребуются для расчета температур Т0(х), Т1(х) и Т11(х):
1) объемные скорости (м3/с) бурового раствора до и ді в нисходящем и восходящем потоках;
2) объемные теплоемкости (в ) бурового раствора сор и с\р в нисходящем и восходящем
потоках соответственно, спор - разбуриваемой породы, св - воды;
3) скрытая объемная теплота плавления льда Ь - ^г’,
4) льдистость разбуриваемой многолетнемерзлой породы си(х) - ^§;
5) коэффициент теплопередачи о. (в м2^™к ) через стенку бурильной колонны от нисходящего потока бурового раствора к восходящему;
6) внешний диаметр бурильной колонны До и разбуриваемой скважины Вс (в м);
7) то(х) > 0 и «ті (х) > 0 - абсолютные величины линейных плотностей распределения сторонних тепловых источников (в нисходящем потоке бурового раствора) и стоков (в восходящем потоке) соответственно - в .
м-с
С помощью заданной плотности то(х) тепловых источников учитывается генерация тепловой энергии работой долота на забое (силами трения бурильной трубы о буровой раствор) и работой сил давления по циркуляции раствора. С помощью заданной плотности ті(х) тепловых стоков учитывается тепловая энергия, которая расходуется на подогрев разбуренной породы до температуры бурового раствора, на подогрев льда до температуры плавления, на плавление льда и на подогрев талой воды до температуры бурового раствора.
Дифференциальные уравнения относительно температур То(х) и Ті(х) бурового раствора соответственно в нисходящем и восходящем потоках получим путем составления уравнений теплового баланса. При этом считаем, что стенка пробуренной скважины теплоизолирована от окружающей ММП. Поэтому вывод дифференциальных уравнений делается при условии, что теплопередача от восходящего потока бурового раствора в окружающую породу отсутствует. Рассмотрим на рис. 1 слой £(х, х + Дх). В нисходящем потоке сечение х будет входным, а сечение х + Дх - выходным из слоя сечением. Для восходящего потока наоборот, сечение х + Дх - входное, а сечение х - выходное из слоя 5.
Количество теплоты фі (Дж), вносимой нисходящим потоком бурового раствора в слой 5(х,х + Дх) за зафиксированный отрезок времени Ді равно фі = То(х)доДі сор. Количество теплоты ^2, вырабатываемой сторонними тепловыми источниками в слое £(х, х + Дх) буровой колонны за зафиксированный отрезок времени Ді, равно ^2 = то(х)ДіДх. Количество тепла, передаваемого от нисходящего потока бурового раствора к восходящему потоку, осуществляется путем теплового контакта потоков через стенку бурильной трубы. Уравнение теплопередачи при контакте тел с разными температурами определяется на основании известного закона Ньютона - Рихмана в виде фз = апДоДхДі[То(х) — Ті(х)].
Разность Дф = (фі + ^2) — фз между поступающим за время Ді в слой 5 теплом фі + ^2 и выходящим из него теплом фз идет на изменение температуры раствора со значения То(х) до значения То(х + Дх). Таким образом, из слоя £(х, х + Дх) буровой раствор выносит тепло, равное ф4 = То(х + Дх) доДі сор.
Уравнение теплового баланса фі + ф2 — ф3 = ф4 после алгебраических преобразований и перехода к пределу при Дх ^ 0, примет вид дифференциального уравнения
^ = ао(х) - атгВ0[Т0(х) - Ті(х)]. (1)
Уравнение теплового баланса для восходящего потока бурового раствора составляется аналогично: <31 + <32 — <3з = <?4, где ф1 - тепло на входе в слой 5, <32 - тепло, поступающее в слой 5 от
нисходящего потока, фз - тепло, отбираемое из слоя 5 тепловыми потоками, ф4 - тепло на выхо-
де из слоя 5. Отсюда после алгебраических преобразований и предельного перехода при Дх ^ 0 получаем дифференциальное уравнение
^ = <7і (х) - атгВ0[Т0(х) - Ті (ж)]. (2)
Уравнения (1) и (2) определяют систему дифференциальных уравнений, описывающую тепловое взаимодействие нисходящего и восходящего потоков бурового раствора при условии, что стенка разбуренной скважины не передает тепла от восходящего потока бурового раствора в окружающую ММП.
Вычитая из уравнения (1) уравнение (2), получим равенство
^То(ж) ^Ті(ж)
Ч0С0р~сїх-------Ч1°1р <1х = °"0^ ~~ °нж)’
интегрируя которое найдем первый интеграл системы уравнений (1), (2):
X
9оСорТо(х) - дісірїі(ж) = J[оо(*0 - аі^)]^ + (7Ь (3)
о
где С71 - произвольная постоянная. Первый интеграл (3) позволяет исключить в системе (1), (2) одну из функций и свести задачу к интегрированию обыкновенного дифференциального уравнения относительно второй неизвестной функции. Например, исключая функцию Т.(ж), относительно функции То (ж) получим дифференциальное уравнение
X
Яосор^у^- + смгД)(1 - /З)То(ж) = сго(ж) - а7Г^° ([а0(г) - ах(г)\<1г + С\, (4)
«Ж д.С.р .]
о
где
Я1С1р Я1С1р
Интегрирование дифференциального уравнения (4) приводит к следующему выражению для распределения температуры нисходящего потока:
Т0(ж) = (1 - е~кх) + С2е~кх+
«осорк
{X XX
[ екгао(г)(1г--------- [ екгс1г [\ао(т) — аит)Ыт
.) Я1С1р ] ]
о 0 0
где С. и С2 - произвольные постоянные,
*=^(1-/3).
«осор
Температура Т. (ж) восходящего потока бурового раствора в воображаемой скважине с теплоизолированной от ММП боковой стенкой по найденной температуре То (ж) вычисляется при помощи интеграла (3). Произвольные постоянные С., С2 в формулах распределения температур То (ж) и Т. (ж) находятся из следующих граничных условий. В точке ж = 0 задается температура Твх входящего в скважину потока бурового раствора, а в точке расположения долота ж = I задается равенство температур нисходящего и восходящего потоков: То|х=о = То(0) = Твх и То(1) = Т.(I).
После того, как температура Т. (ж) восходящего потока бурового раствора в воображаемой скважине с теплоизолированной боковой стенкой будет найдена, приступаем к расчету методом аппроксимации температуры Тц(ж) в восходящем потоке бурового раствора для реальной скважины. Для этого вначале по найденной температуре Т. (ж) вычисляем температуру Т*ых бурового раствора на выходе из скважины, которая была бы в условиях термостатического теплообмена Т*ых = Т.(0). Затем по фактическому значению температуры Твых бурового раствора на выходе и теоретическому значению Т*ых вычисляем переменный поправочный коэффициент
сн = ^ + (1-^)г. (5)
Т вых \ Т вых/ 1
Значение коэффициента С (ж) на выходе из скважины равно, согласно (5), отношению температур Твых/Тв“ых, а в точке ж = I расположения бурового долота равно единице. Во всех же промежуточных точках значения поправочного коэффициента С(ж) по крайним значениям вычисляются методом линейной аппроксимации. Приближенное к фактическому распределение температуры Тц(ж) по температуре Т. (ж) и поправочному коэффициенту £ (ж) вычисляем по формуле
Тц(ж) = С (ж)Т.(ж). (6)
Качественные графики распределения температур То(ж), Т.(ж) и Тц(ж), полученных из системы дифференциальных уравнений (1), (2) и формулы (6) в конкретном тестовом примере, выше были приведены на рис. 2.
Расчет радиусов протаивания ММП вокруг строящихся скважин по рассчитанным температурам Т. (ж) и Тц(ж) осуществляется по методике, ранее опубликованной авторами в работе [1].
Список литературы
1. Гасумов Р. А., Кондренко О. С., Толпаев В. А. Расчет изотермического фронта протаивания многолетнемерзлой породы по данным, замеряемым при проводке скважины // НТЖ «Строительство нефтяных и газовых скважин на суше и на море». №2, 2011. С. 37-41.
Рукопись поступила в редакцию 18 апреля 2011 г.
