CC BY
83
УДК 622.276.03
АЛ Ковалев
Фильтрационно-прочностной расчет окрестности ствола вертикальной скважины с использованием критерия пластичности Друкера-Прагера
Ключевые слова:
модель
материала пласта, фильтрационнопрочностной расчет, напряженнодеформированное состояние, напряжения, деформации.
Keywords:
model of stratum material,
filtration-strength
calculation,
stress-deformed
condition,
stresses,
deformations.
При выборе конструкции забоя скважины в слабосцементированном коллекторе необходимо учитывать возможность разрушения пород, примыкающих к забойному оборудованию, и вынос частиц разрушенного материала с добываемой продукцией. На практике при планировании бурения на новом объекте в первую очередь руководствуются накопленным опытом эксплуатации скважин в схожих условиях. Параллельно ведутся экспериментальные и теоретические исследования по поиску математических моделей неупругого поведения горной породы и совершенствованию методов фильтрационно-прочностных расчетов на основе испытаний керна и данных геофизических исследований (ГИС) о составе пород продуктивного разреза. Важной частью этих исследований являются вычислительные эксперименты, для проведения которых создают специальные исследовательские модели и программы.
В настоящей статье рассматриваются модели горной породы в окрестности необсаженного ствола вертикальной скважины. Модели созданы в Центре исследований нефтегазовых пластовых систем и технологического моделирования ООО «Г азпром ВНИИГАЗ» и являются частью расчетного инструментария, используемого при фильтрационно-прочностной оценке конструкций скважин на конкретных объектах. Пористая порода предполагается изотропной, принимается цилиндрическая система координат, вертикальная ось которой совпадает с осью скважины. Рассматриваются установившиеся течения флюида. Предполагается, что граничные условия не зависят от угловых координат. В силу перечисленных допущений все искомые функции зависят от одной переменной - радиуса г.
Математическая модель
Рассматривается область пласта, ограниченная двумя соосными цилиндрическими поверхностями и двумя плоскостями, ортогональными оси цилиндров. Внутренняя цилиндрическая поверхность соответствует поверхности необсаженного ствола скважины, внешняя - выбирается таким образом, чтобы на ней можно было считать затухшим геомеханическое возмущение от вынутой из ствола породы и фильтрационного потока к нему:
Rw < г < R, (1)
где Rw, Rf - радиусы внутренней и внешней границ области.
Законы сохранения массы и количества движения в случае установившегося притока флюида к цилиндрической полости принимаются в виде
pw2nrh = pQ = const,
dP u oil — = ——w -ppww. dr k ' 1
(2)
(3)
где P - пластовое (поровое) давление; w - скорость фильтрации; h - толщина пласта; p, р - плотность и вязкость; к, р - коэффициенты проницаемости и турбулентности (коэффициент Форхгеймера).
№ 2 (18) / 2014
Актуальные вопросы исследований пластовых систем месторождений углеводородов
37
Закон (3) соответствует однофазной фильтрации, но излагаемая расчетная схема позволяет при необходимости учесть и многофазную фильтрацию, для чего следует заменить (3) соответствующей моделью.
В принятой системе координат действуют три главных эффективных (по Терцаги [1]) напряжения - радиальное (стД тангенциальное (стф) и осевое (или вертикальное) (ст2). Суммарный баланс количества движения для флюида и пористой среды записывается в виде (положительными при этом приняты растягивающие напряжения):
dQr _ dp + Qr -Оф _ о dr dr r
(4)
Задача решается для условий плоской деформации. Тензор полных деформаций коак-сиален тензору напряжений и полностью определяется перемещением вдоль радиуса - u(r). Радиальная, тангенциальная, осевая и объемная деформации определяются соотношениями
du
£ dr ’ ф
е _ £r +£ф+£
z 0
du +u dr r
+ £z
(5)
где e - деформация вдоль соответствующих осей; ez0 - начальная деформация вдоль вертикальной оси (в классических задачах плоской деформации принимается равной нулю); e - объемная деформация.
В области упругого состояния напряжения связаны с деформациями законом Гука-Терцаги-Био:
or = Xe + 2Ger + psKP, (6)
стф = Xe + 2GS(p + р^КР, (7)
CTz = Xe + 2Gezo + PSKP, (8)
где X - коэффициент Ламе; G - модуль сдвига; К - модуль объемного сжатия; ps - сжимаемость материала скелета породы.
Граница упругого и пластического состояний (поверхность течения) описывается критерием Кулона-Мора.
fc-M = М + tg9CT„ - sc = ° (9)
где х и сти - соответственно, касательное и нормальное напряжения, приложенные к площадке
скольжения; Sc - когезия (модуль сцепления); Ф - угол внутреннего трения. Круг Мора для напряжений, расположенных на поверхности течения, должен касаться прямой (9).
В зоне пластического состояния породы полная деформация складывается из двух частей - упругой, подчиняющейся закону Гука (6)-(8), и остаточной пластической.
£ _ £e +£p £ _ £e + £p £ _ £e +£p
cr cr^cr’ °ф сф^сф’ cz cz^cz-
(10)
По условию плоской деформации
£z _ £Z + £ Z _£ z
(10a)
Для определения пластических деформаций принят ассоциированный закон, согласно которому эти деформации связаны с поверхностью течения соотношениями
£Р _ Xf, £ф _ X fp _ X,f даф
do,
do.
(11)
где X - коэффициент пропорциональности.
Попутно отметим, что в описываемой расчетной схеме легко может быть использован и неассоциированный закон пластического течения, для чего в (11) функция пластичности (9) заменяется пластическим потенциалом.
Если взять производную по радиусу от тан-гециальной деформации (5) и исключить перемещение, то в результате получим условие совместимости
d £ф dr
£r _£ф
r
(12)
Использование напрямую критерия Кулона-Мора в радиальной задаче сопряжено с определенными трудностями. Дело в том, что касательное и нормальное напряжения в (9) необходимо выразить через два главных напряжения - минимальное и максимальное. Когда главных напряжений всего два, как в задаче со сферической симметрией, проблем с использованием критерия Кулона-Мора не возникает [2]. В радиальной задаче главных напряжений три, минимальным из них (по абсолютной величине) для условий прискважинной области всегда будет радиальное [3]. Остается выбрать из оставшихся двух (тангенциального и осевого) максимальное по абсолютной величине напряжение. Этот выбор может быть серьезно осложнен тем, что в пределах зоны пластических деформаций возможна смена в качестве
№ 2 (18) / 2014
38
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
максимального одного напряжения на другое. В работах [3, 4] эта проблема решается допущением существования в зоне пластических деформаций некой внутренней зоны, в которой тангенциальное напряжение равно осевому. В настоящей статье была использована аппроксимация критерия Кулона-Мора, известная как соотношение Друкера-Прагера [1], в котором задействованы все три главных напряжения. Для рассматриваемого случая соотношение будет иметь вид:
щуюся жидкость считать несжимаемой (слабосжимаемой) и пренебречь сжимаемостью породы, то указанную систему с учетом (5) можно преобразовать к виду
JP = JH_1 + PP||^_ (14)
dr 2nkh r 4n h r
du
dr
о -X\ -+£
X + 2G
(15)
fD - P + \/J2 C,
I = O + O+Oz .
J2 = 6 [(O -0Ф )2 +(0Ф -°z )2 + (Oz -O )2 ] >
a =
2sin ф
6 SC cos ф V3(3 ± sin ф) V3(3 ± sin ф)
C =
(13)
где I1 и J2 - первый инвариант тензора напряжений и второй инвариант девиатора тензора напряжений. Знак («+» или «-») в знаменателе выражений для коэффициентов показывает, какая именно поверхность Друкера-Прагера используется для аппроксимации - описанная («+») или вписанная («-»).
Применение соотношения Друкера-Пра-гера к осесимметричной задаче также вызывает больше трудностей, чем для сферы [5], поскольку связь между напряжениями на поверхности пластического течения перестает быть линейной. Поэтому для решения настоящей задачи был использован численный метод Рунге-Кутта.
Для замыкания математической модели считалось, что заданы поровые давления, радиальное напряжение на внешней границе области (1) и расход флюида. В общем случае существует граница r = R разделяющая область упругого состояния материала Rp < r < Rf и область с остаточными деформациями Rw < r < Rf. Возможны также случаи, когда одна из зон отсутствует и вся область находится в упругом или пластическом состоянии. Для того чтобы выяснить, какой из случаев имеет место, необходимо попытаться «сшить» на границе решения для упругой и пластической зон. При этом для зон решались следующие системы уравнений.
Уравнения для упругой зоны
Напряженно-деформированное состояние в упругой зоне описывается системой дифференциальных уравнений (2)-(4). Если фильтрую-
dо dP ог -оФ
—- =-------r----. (16)
dr dr r
Чтобы замкнуть систему (14)-(16), необходимо добавить выражение для тангенциального напряжения, которое может быть получено из (7) с учетом (5):
°ф=^ [dr+Его)+a+2G) и • (17)
Для разрешимости системы (14)-(17) к названным выше граничным условиям необходимо добавить радиальное перемещение на внешней границе. Поскольку такое перемещение неизвестно, то задавалось радиальное напряжение на границе упругой и пластической зон, а перемещение подбиралось таким, чтобы решение соответствовало этому и остальным граничным условиям. В частном случае (при отсутствии зоны пластических деформаций) граничным будет нулевое радиальное напряжение на стенке скважины.
Уравнения для пластической зоны
В систему уравнений для пластической зоны также входят уравнения (14) и (16), а кроме того уравнение для тангенциальной деформации, полученное из (12) подстановкой в него выражения для радиальной упругой деформации:
dе itO- -v(O+Oz)] + £p -£ф
—ф = E-----------------------, (18)
dr r
где E - модуль Юнга; v - коэффициент Пуассона.
Чтобы замкнуть систему (14), (16), (18), необходимо добавить выражения для оф, oz, и е£. В качестве таковых используется следующая система уравнений:
оД + л]J2 - C = 0, (19)
№ 2 (18) / 2014
Актуальные вопросы исследований пластовых систем месторождений углеводородов
39
{E[cz -v(ar +оф)]-ezо j(6аJ + 2аг -Оф -o) +
+ -p(6a J-ar -Оф + 2oz) = 0, (20)
{£ф - ^ I[°ф - V(0r + 0z )Ц + 2Ог -Оф - О z ) -
-ep (6aJ-or + 2Оф - oz) = 0. (21)
Уравнение (19) получено из (13) для условия пластического течения. Уравнения (20) и (21) - из (10а) и (10) с учетом (11) и (12). Система уравнений (19)-(21) является нелинейной и решается итерационным методом Ньютона.
Система уравнений (14), (16), (18)-(21) решается методом Рунге-Кутта от внутренней границы (стенки скважины). Для ее разрешимости на этой границе достаточно задать еф (радиальное напряжение здесь равно нулю, пластовое давление известно, если фильтрационные свойства пласта не зависят от напряжений). Для лучшей сходимости численного метода значения остальных напряжений, входящих в (16) и (18), на внутренней границе можно определить, составив и решив относительно сф и cz систему из двух уравнений - уравнения (19) и следующего:
E[оz -УОф ] (6а J + 2Оф - оz) + Еф (6а J -Оф + 2оz) -
- E [оф - уо z ] (6 aJ - Оф + 2о z ) + e z о (6 aJ + 2Оф - о z ) = 0. (22)
Данная расчетная схема справедлива для условий идеальной пластичности, но легко может быть адаптирована и для модели пластичности с упрочнением (билинейной модели). В этом случае когезия (сцепление) не постоянна, а зависит от накопленной в данной точке пластической деформации
Sc (х) = SC + (srk - SC) -f-,
—
max
X = V (-P )2 + (-ф )2 + (-P )2, (23)
где SC - модуль сцепления в состоянии перехода от упругих к пластическим деформациям; SCeak - модуль сцепления в состоянии максимального напряжения; х - параметр, характеризующий суммарную пластическую деформацию; -Pdx - максимальная пластическая деформация (деформация разрушения).
Совместное решение для упругой и пластической зон
Как указано выше, параметры пластической зоны можно определить полностью, задав на внутренней границе (стенке скважины) значение единственной переменной -полной тангенциальной деформации. Чтобы «сшить» решение для пластической зоны с решением для упругой зоны, необходимо определить еще один параметр - радиус границы между зонами. Условиями «сшивки» будут являться попарное равенство на границе всех трех напряжений (со стороны каждой из зон) и равенство нулю на границе пластических деформаций.
Указанная задача решалась процедурой подбора названных переменных - полной тангенциальной деформации на стенке скважины и радиуса границы между упругой и пластической зонами - до достижения с заданной точностью равенства напряжений на границе зон и равенства нулю пластической тангенциальной деформации на той же границе.
№ 2 (18) / 2014
40
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
Примеры расчета
В завершение приведем две серии расчетов. В первой серии сопоставляются три расчета (I, II и III) с идеальной пластичностью, отличающиеся уровнем вертикального (осевого) напряжения на внешней границе. Во второй серии один из вариантов первой серии сопоставляется с соответствующим по осевому напряжению вариантом с упрочнением (билинейная модель). В аппроксимации критерия пластичности использована описанная поверхность Друкера-Прагера. Для простоты в примерах рассматривается нелинейная фильтрация несжимаемой жидкости. Исходные данные для расчетов приведены в таблице. Результаты расчетов представлены на рис. 1-6.
На рис. 1 сопоставляются напряжения в трех расчетах с идеальной пластичностью. Наглядно отображено, что с повышением осевого напряжения на внешней границе увеличивается радиус зоны пластических деформаций. Для осевого напряжения -330,6 бар он составляет 1,00 м, для -430,6 бар - 1,12 м, для -490,6 бар - 1,40 м. В первом случае максимальным главным напряжением на протяжении всей пластической зоны остается тангенциальное напряжение, в двух других на границе упругой и пластической зон максимальным главным является осевое напряжение, в дальнейшем (по мере приближения к стволу скважины) его в этом качестве меняет тангенциальное напряжение.
На рис. 2 и 3 для указанных расчетов сопоставляются пластические деформации и перемещения. С увеличением осевого напряжения наблюдается рост этих параметров, и для напряжения -490,6 бар они достигают, соответственно, значений 141 (тангенциальная деформация) и 5,46 м (перемещение). Столь высокие значения деформаций и перемещений указывают на необходимость при использовании модели идеальной пластичности введения связанных с этими параметрами критериев разрушения материала пласта.
На рис. 4-6 сопоставляются напряжения, деформации и перемещения в варианте пластичности с упрочнением (билинейная модель) с осевым напряжением на внешней границе -430,6 бар и соответствующем по этому параметру варианте идеальной пластичности. Отображено, что в модели с упрочнением по сравнению с моделью идеальной пластичности сокращается радиус зоны пластичности (до 0,4 против 1,12 м), тангенциальное напряжение остается на том же уровне, многократно уменьшаются деформации и перемещения (соответственно, до 2,84 и 0,11 м против 112,84 и 4,37 м). В отличие от модели идеальной пластичности в модели с упрочнением уже есть критерий разрушения материала - деформация разрушения, по которому можно определить радиус зоны разрушения.
Исходные данные для примеров
Наименование параметра Единица измерения Значение
ИП I ИП II ИП III БМ
Радиус скважины м 0,1
Радиус внешней цилиндрической границы м 10
Толшцна пласта м 10
Проницаемо сть Дарси 5
Коэффициент Форхгеймера 1/см 639504,3
Дебит жидкости м3/сут 676
Вязкость жидкости сП 2
Плотность жидкости кг/м3 800
Пластовое давление на внешней границе бар 200
Горное давление на внешней границе бар -460
Вертикальное (осевое) напряжение на внешней границе бар -330,6 -430,6 -490,6 -430,6
Модуль Юнга бар 20000
Коэффициент Пуассона б/р 0,25
Модуль сцепления при переходе в пластическое состояние бар 2
Модуль сцепления при разрушении бар - - - 4
Угол внутреннего трения О 25
Деформация разрушения б/р - - - 0,25
Примечание'. ИП - идеальная пластичность, БМ - билинейная модель.
№ 2 (18) / 2014
Напряжение, бар Напряжение, бар Напряжение, бар
ii
iii
Рис. 1. Распределение ог, и ог в расчетах с идеальной пластичностью (i-iii)
i
№ 2 (18) / 2014
Напряжение, бар Перемещение по радиусу, м Деформация, б/р
42
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
Рис. 2. Распределение пластических деформаций £rPL, £фРЬ и £zPL в расчетах с идеальной пластичностью (I-III)
Рис. 3. Распределение и в расчетах с идеальной пластичностью (I—III)
Радиус, М
Рис. 4. Сопоставление распределения or, и oz в расчетах для идеальной пластичности и пластичности с упрочнением
№ 2 (18) / 2014
Актуальные вопросы исследований пластовых систем месторождений углеводородов
43
Рис. 5. Сопоставление распределения пластических деформаций zrPL, гфРЬ и £zPL в расчетах для идеальной пластичности и пластичности с упрочнением
Рис. 6. Сопоставление и в расчетах для идеальной пластичности и пластичности
с упрочнением
В заключение на основании рассмотренного материала можно сделать следующие выводы:
• автором предложен комплекс методик для фильтрационно-прочностного расчета напряжений и деформаций, возникающих в горной породе, окружающей необсаженный ствол вертикальной скважины, под действием установившегося течения пластового или закачиваемого флюида;
• указанные методики позволяют учесть линейный и двучленный законы фильтрации, а также многофазность потока;
• расчеты напряженно-деформированного состояния включают определение пластических деформаций, что позволяет использовать различные критерии разрушения, в том числе деформационные.
Список литературы
1. Fjaer E. Petroleum related rock mechanics /
E. Fjaer, R.M. Holt, P Horsrud et. al. -Amsterdam, Oxford: Elsevier, 2008. - 491 p.
2. Bartli R.K. Stability and Failure of Sand Arches / R.K. Bartli, R. Risnes // SPEJ. - April 1981. -
P. 236-248.
3. Risnes R. Sand Stresses Around a Wellbore /
R. Risnes, R.K. Bartli, P. Horsrud // SPEJ. -December 1982. - P. 883-898.
4. Пятахин М.В. Геомеханические проблемы при эксплуатации скважин / М.В. Пятахин. - М.:
Г азпром ВНИИГАЗ, 2011. - 266 с.
5. Ковалев А.Л. Фильтрационно-прочностной расчет окрестности перфорационной каверны / А. Л. Ковалев, Е.В. Шеберстов // Актуальные вопросы исследований пластовых систем месторождений углеводородов: сб. науч. ст. -М.: Газпром ВНИИГАЗ, 2012. - С. 117-129. -(Серия «Вести газовой науки»).
№ 2 (18) / 2014
