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

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

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


Pages:     | 1 | 2 || 4 | 5 |   ...   | 6 |

Научные основы методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений нефти и газа

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

Полученные зависимости в дальнейшем учитывались при математическом моделировании НДС коллекторов.

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

= nп + тв,

где n – открытая пористость.

Рис.6. Примеры компрессионных кривых ряда месторождений.

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

Компрессионная зависимость дает полное представление о деформируемости порового пространства, т.к. с ее помощью можно найти коэффициент сжимаемости пор в заданном интервале всестороннего давления:

п = n/(n) = e/(e),

где n - пористость, e - коэффициент пористости.

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

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

4. Прогноз напряженно-деформированного состояния горных пород

при разработке месторождений углеводородов

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

где tij - полные напряжения; K, G – модуль объемного сжатия и модуль сдвига; р – поровое давление; - коэффициент разгрузки.

Для полной характеристики упругой пористой среды уравнения состояния (12) должны быть дополнены уравнением фильтрации находящегося в порах флюида:

где – модуль Био; k, - проницаемость породы и вязкость флюида; v – объемная деформация породы.

В наиболее общем варианте для расчета напряженно-деформированного состояния пористой среды необходимо совместное решение уравнений теории упругости (12) и уравнения фильтрации (13). Однако в ряде случаев, при определенных видах напряженного состояния, можно упростить вид уравнения фильтрации и решать его независимо от уравнений теории упругости. Например, если деформирование пород происходит с сохранением постоянного объема (v=0), то уравнение фильтрации преобразуется к виду

где f – сжимаемость флюида; 1 - сжимаемость породы при v=0.

Другие характерные виды напряженного состояния возникают при одномерном уплотнении (x=y=0, tz = const), при плоской деформации (z= 0, tx = ty = const) и при постоянстве среднего давления: (tx + ty + tz)/3 = const. В этих случаях в уравнении фильтрации будут фигурировать коэффициенты сжимаемости породы 2, 3 и 4. Т.е., уравнение фильтрации в общем виде записывается

где r – некоторый постоянный коэффициент сжимаемости породы. При этом можно показать, что при произвольном напряженном состоянии выполняется соотношение 1 r 4. Уравнение (14) решается независимо от уравнений теории упругости, т.к. не содержит в явном виде v. Коэффициент сжимаемости флюида обычно значительно больше сжимаемости породы, поэтому общая сжимаемость системы флюид-порода общ=(f+r) определяется в основном коэффициентом f. В качестве примера в табл.2 представлены коэффициенты сжимаемости, рассчитанные для условий, характерных для нефтяных месторождений севера Пермского края.

Таблица 2 - Коэффициенты сжимаемости системы флюид-порода

Сжимаемость породы r, МПа-1 1 2 3 4
3,2510-5 1,21410-4 1,65810-4 1,92510-4
Общая сжимаемость общ, МПа-1 1,03210-3 1,12110-3 1,16610-3 1,19210-3

Можно видеть, при изменении коэффициента сжимаемости породы от минимального значения 1 до максимального 4 общая сжимаемость системы общ изменяется довольно незначительно – на 15%, а для наиболее вероятных значений (от 2 до 3) изменяется еще меньше – на 5%. Подобное соотношение между коэффициентами сжимаемости характерно для большинства практических случаев, поэтому гидродинамическое моделирование разработки нефтяных и газовых месторождений производится на основе уравнений вида (14), т.е. деформационные свойства горных пород учитываются применением соответствующих коэффициентов сжимаемости r.

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

Наиболее простой вид напряженного состояния при добыче нефти и газа возникает при одномерном уплотнении, т.е. при отсутствии горизонтальных деформаций (x=y =0) и постоянстве полных вертикальных напряжений (tz = const). Для данных условий из уравнений закона Гука (12) можно найти вертикальную деформацию уплотнения:

где см - коэффициент одномерного уплотнения; р – падение давления.

Уравнение (15) можно применять для глубоко залегающих плотных низкопористых пород, деформации которых при больших давлениях с достаточной степенью точности соответствуют теории упругости. В более общем случае необходимо учитывать специфические законы деформации пористых сред. Работы К. Терцаги, М. Био, Ф. Гасмана, В.М.Добрынина показывают, что объемные деформации пористой среды складываются из деформаций твердой фазы и порового пространства при изменении внешнего всестороннего давления t и давления жидкости в порах:

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

Слабой стороной методики использования эффективных напряжений является неопределенность, связанная с учетом изменения коэффициента разгрузки при росте всестороннего давления. Однако в большинстве практических случаев зависимость коэффициента разгрузки от давления можно не учитывать, поэтому в расчетах использовались эффективные напряжения. Сжимаемость твердой фазы коллекторов принималась по табличным данным (для песчано-глинистых коллекторов тв = 0,030 ГПа-1, для карбонатных – тв = 0,025 ГПа-1), сжимаемость пор определялась по результатам компрессионных испытаний на основе зависимостей (1):

С учетом этих соотношений можно найти объемную деформацию коллектора при изменении эффективного давления от 1 до 2:

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

где h – уплотнение; h – начальная толщина коллектора; z,1, z,2 – начальное и конечное значение эффективных вертикальных напряжений.

Специальные сопоставительные расчеты показали, что использование табличных значений для сжимаемости твердой фазы коллекторов обеспечивает достаточную для практических целей точность. Даже при неучете сжимаемости твердой фазы погрешность расчета вертикальной деформации обычно не превышает 10%. Это подтверждает правомерность использования механических моделей «шатрового» типа. Поскольку величина одномерного уплотнения представляет собой верхний предел возможных оседаний поверхности, то полученные аналитические зависимости позволяют сделать предварительные выводы об общем уровне напряженного состояния массива. Это дает возможность в дальнейшем более обоснованно подходить к построению общих расчетных моделей изучаемых объектов.

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

{}={el}+{vp}.

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

где - параметр вязкости, играющий роль регулятора сходимости итерационного процесса. Непосредственно деформации находятся путем численного интегрирования (22) по времени:

Реализация моделей была выполнена в конечно-элементных программах “GEOTECH” (разработка ПермГТУ) и “ANSYS”.

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

{t} = {} + {m}p, (24)

где { t } = { t x, t y, t z, t xy, t yz, t zx}т – полные напряжения;

{} = {x, y, z, xy, yz, zx} т – эффективные напряжения;

{m} = {1, 1, 1, 0, 0, 0} т.

С учетом (24) можно получить общее соотношение МКЕ в виде:

[Ke ] {e} = {Fe} - {Fep}, (25)

где {Fep} – вектор дополнительных узловых сил, обусловленных давлением насыщающей жидкости.

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

Для анализа полного характера уплотнения коллекторов была выполнена серия расчетов осесимметричной модели, в которых варьировалось отношение G/Gk (Gk - модуль сдвига пород коллектора, G - модуль сдвига вмещающих пород) и отношение R/H (R – радиус коллектора, H – глубина залегания). По результатам расчетов строились графики относительного уплотнения (h/hmax) в зависимости от текущего радиуса r (рис.7) для анализа характера уплотнения на всем протяжении продуктивных пластов (максимальное значение hmax соответствует условиям одномерного уплотнения).

 Влияние упругих свойств вмещающих пород и геометрических размеров -23

Рис. 7. Влияние упругих свойств вмещающих пород и геометрических размеров

коллекторов на характер их уплотнения.

Было установлено, что уплотнение коллекторов в целом выше при более слабых вмещающих породах, особенно на границах области снижения давления. В то же время в центральной части уплотнение близко к максимальному и слабо зависит от жесткости вмещающих пород. Уплотнение также увеличивается с ростом радиуса коллектора, однако уже при R/H 2 влияние этого фактора становится несущественным. Наконец, относительное уплотнение коллектора зависит от его начальной толщины. При прочих равных условиях относительное уплотнение выше для более тонкого коллектора, что проявляется на флангах залежи; в центральной части влияние этого фактора практически отсутствует. В целом результаты данных расчетов показали, что даже при небольшом отношении R/H деформации коллекторов в центральной части нефтегазовых месторождений близки к условиям одномерного уплотнения. Однако в краевых частях области снижения исходного пластового давления условия одномерного уплотнения не выполняются и напряженное состояние имеет более сложный вид. Оседания земной поверхности, кроме вышеназванных параметров, зависят также от соотношения жесткости покрывающих и подстилающих пород. Наибольшие оседания земной поверхности имеют место при более жестких подстилающих породах. В целом степень влияния упругих свойств вмещающих пород на уплотнение коллекторов и оседания земной поверхности уменьшается с ростом размеров коллектора, т.е. деформации горного массива определяются прежде всего величиной уплотнения продуктивного слоя.

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

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

где RD, rD, рD и tD –безразмерные величины: RD = Rk/Rc; rD = r/Rc; tD = kt / (nобщR2c); pD = 2kh [p0-p] / (Q).

В этих формулах: r, t – текущий радиус точки и текущее время работы скважины; Rk- радиус контура питания; Rc – радиус скважины; k – проницаемость пласта; h – толщина пласта; - вязкость флюида; n – пористость; общ – коэффициент сжимаемости системы флюид-порода; p0, р – начальное и текущее пластовое давление.

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

 Пластовое давление и уплотнение коллектора в зоне влияния скважины Следующий-25

Рис. 8. Пластовое давление и уплотнение коллектора в зоне влияния скважины

Следующий рассматриваемый фактор связан с тем, что коллектора нефти и газа имеют куполообразное строение, т.к. приурочены к сводовым частям геологических структур-ловушек углеводородов. Поскольку углы наклона крыльев складок составляют первые градусы или даже минуты, расчетные модели обычно строят в виде прямолинейных пластов с постоянной глубиной залегания. Влияние подобной идеализации геометрии залежи было рассмотрено на численном примере коллектора толщиной 50 м, радиусом R=1500 м, залегающего на глубине Н=1500 м (рис.9). Предполагалось, что водонефтяной контакт расположен на 100м ниже купольной части, что при данных размерах залежи является довольно большой величиной. Рассчитывались оседания земной поверхности, которые сравнивались с аналогичными величинами для коллектора плоской формы с постоянной глубиной залегания Н=1500 м.

 Фрагмент конечно-элементной модели коллектора куполообразной формы В базовом-26

Рис. 9. Фрагмент конечно-элементной модели коллектора куполообразной формы



Pages:     | 1 | 2 || 4 | 5 |   ...   | 6 |
 





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

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