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

Abstract

Актуальность исследования обусловлена необходимостью расчета напряженно-деформированного состояния околоскважинной зоны при разработке нефтяных и газовых месторождений для прогноза устойчивости обсадной колонны и сохранности цементного камня. Цель: на основе сопоставительного анализа методов аналитического и численного моделирования расчета напряжений вблизи вертикальной скважины с применением упругой модели определить распределение радиальных и угловых напряжений, сопоставить точность их вычисления разными методами и выявить достоинства и недостатки каждого из них. Объекты: околоскважинная зона терригенного коллектора ачимовских отложений одного из месторождений Ханты-Мансийского автономного округа. Методы: аналитический и численный конечно-элементный методы расчета напряженно-деформированного состояния околоскважинной зоны с учетом основных конструктивных элементов скважины и с использованием линейно-упругой модели. Результаты. Рассмотрены аналитические соотношения, используемые для расчета радиальных и угловых напряжений в колонне, цементном камне и породе-коллекторе, а также уравнения, применяемые при численном конечно-элементом моделировании напряжений вблизи вертикальной скважины. Разработана конечно-элементная схема околоскважинной зоны, включающая ее основные конструктивные элементы. Представлены результаты расчета основных компонент тензора напряжений в конструктивных элементах скважины в зависимости от радиальной координаты для величин забойного давления 20, 40 и 60 МПа. Проведен сопоставительный анализ результатов расчетов напряжений указанными методами. Показано, что наибольшее расхождение между аналитическим и численным методом составило 2 %, что соответствует радиальным напряжениям для варианта расчета с забойным давлением 20 МПа. В среднем расхождения составили: для радиальных напряжений - 0,04 %, для угловых - 0,72 %. Сделан вывод о том, что при использовании модели линейно-упругой среды и при граничных условиях в виде закрепления модели в верхней и нижней части по нормали к поверхности, а также без учета распределения давления в депрессионной воронке модели достаточно использовать аналитический метод расчета. Если же предполагается использовать комбинированные граничные условия, пороупругую модель, учет вязкопластических деформаций, то наиболее предпочтительно использовать метод численного моделирования.

The relevance of the research is caused by the scientific interest in calculating the near-wellbore zone stress-strain state to improve the development of oil and gas fields, to predict the casing stability and the safety of the cement stone. The main aim: based on a comparative analysis of the methods of analytical and numerical simulation of stress calculation near a vertical well using an elastic model, determine the distribution of radial and tangential stresses, compare the accuracy of their calculation by different methods and identify the advantages and disadvantages of each of them. Objects: near-wellbore zone of the terrigenous reservoir of the Achimov deposits of one of the fields of the Khanty-Mansiysk autonomous region. Methods: analytical and numerical finite element methods for stress-strain state calculating of the near-wellbore zone, taking into account the main structural elements of the well and using a linear elastic model. Results. The paper considers the analytical relationships used to calculate the radial and tangential stresses in the casing, cement stone and reservoir rock, as well as the equations used in numerical finite element modeling of stresses near a vertical well. The authors have developed a finite element scheme of the near-wellbore zone, including its main structural elements. The paper introduces the results of calculation of the main components of the stress tensor in the structural elements of the well depending on the radial coordinate for the bottom hole pressure of 20, 40 and 60 MPa. The authors carried out the comparative analysis of the stress calculations results by the methods used. It is shown that the largest discrepancy between the analytical and numerical methods was 2 % that corresponds to radial stresses for the calculation option with a bottom hole pressure of 20 MPa. On average, the discrepancies were: for radial stresses - 0,04 %, for tangential stresses - 0,72 %. It is concluded that when using the model of a linearly elastic medium and under boundary conditions in the form of fixing the model in the upper and lower parts along the normal to the surface, and also without taking into account the pressure distribution in the depression funnel of the model, it is sufficient to use the analytical method of calculation. If it is supposed to use combined boundary conditions, a poroelastic model, taking into account viscoplastic deformations, then it is most preferable to use the numerical simulation method.