Оценка эффективности гидроразрыва пласта с двумя трещинами в окрестности одиночной скважины

оригинальная статья DOI: https://doi.org/10.18599/grs.2018.2.108-114

Оценка эффективности гидроразрыва пласта с двумя трещинами в окрестности одиночной скважины

Ю.И. Яхина

ООО «Термосим», Казань, Россия E-mail: vib-210@mail.ru

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

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

Для цитирования: Яхина Ю.И. (2018). оценка эффективности гидроразрыва пласта с двумя трещинами в окрестности одиночной скважины. Георесурсы, 20(2), ^ 108-114. DOI: https://doi.org/10.18599/grs.2018.2.108-114

уДК 532.546:622.276

Введение

с целью вовлечения в разработку недренируемых и слабодренируемых запасов полезных ископаемых применяется гидравлический разрыв пласта (грп). Данная технология способствует созданию высокопроводимой трещины, что позволяет увеличить продуктивность добывающих или приёмистость нагнетательных скважин. По экспертным оценкам гидродинамических исследований (ГИс) гидравлический разрыв продуктивного пласта является наиболее результативным геолого-техническим мероприятием, гарантирующим повышение эффективности разработки плохо проницаемых пластов. одно из первых исследований стационарного однофазного притока к скважине с одиночной трещиной грп представлено в работе (Prats, 1961). Теоретические расчеты показывают, что ГрП позволяет в несколько раз увеличить дебит скважин за счет вовлечения в эксплуатацию слабодренируе-мых зон, причем объем добычи или закачки зависит от проводимости и длины трещины (Cinco-Ley et. al., 1978; Meehan et. al., 1989).

обзор ряда других публикаций и анализ эффективности ГрП в скважине с одной вертикальной трещиной в круговом контуре питания выполнен в работе (Морозов, 2016) на основе различных оценок продуктивности и обобщенного понятия псевдоскин-фактора (Economides et. al., 2002).

При существенном снижении темпов добываемой продукции относительно первоначальных режимов проводят

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

Повторный гидравлический разрыв пласта является одним из наиболее распространенных методов интенсификации добычи запасов полезных ископаемых и повышения продуктивности пласта, что делает актуальной задачу анализа и прогноза его эффективности. В частном случае, задача о стационарном притоке к скважине с несколькими вертикальными трещинами гидроразрыва при равномерно распределенном притоке вдоль трещин в области, ограниченной круговом контуром питания, рассматривалась в работе (Raghavan, Joshi, 1993). Аналогичная задача о стационарном равномерно распределенном притоке к двум трещинам с различными геометрическими характеристиками при заданном угле пересечения решалась численно в (Лихтарев, Пестриков, 2010).

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

© 2018 Ю.И. Яхина

SCIENT1FK AND TECHNICAL JOURNAL

Ш GEORESOURCES www.geors.ru

Общее интегральное представление возмущения поля пластового давления

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

Основной характеристикой фильтрационных процессов является поле пластового давления, которое в случае неустановившейся упругой фильтрации в пористой среде с распределенными источниками описывается уравнением пъезопроводности (Чарный, 1963): дР а*

где х = к/(цвт0) - коэффициент пъезопроводности, определяющий скорость распространения возмущений давления Р(х,у,2, ¿) в пористой среде с проницаемостью к, пористостью т0, коэффициентом упругоемкости в, вязкостью жидкости ц; щ*(х,у,х, ¿) - плотность распределенных источников, которая представляет собой количество жидкости, поступающей в единицу объема среды за единицу времени ^ Д - дифференциальный оператор Лапласа в системе декартовых координат (х,у^).

Пусть в начальный момент времени возмущение давления в пласте отсутствует, тогда решение уравнения пъезопроводности при нулевых начальных условиях можно записать в общей интегральной форме с помощью метода мгновенных точечных источников (Карслоу, Егер, 1964; Тихонов, Самарский, 1999):

£ +оо

р{х,у,г,€) = I т

т0р

— тdr, dr,

где С*(^,п,С,х,у,г,1 - т) - функция влияния мгновенного точечного источника, которая представляет собой возмущение давления в точке (х,у^) в момент времени вызванное единичным мгновенным точечным источником в точке в момент времени / = т:

V. С. У. г, г-т) = ( 1

Общий (суммарный) приток или отток жидкости 0 в скважине с трещинами за единицу времени по определению равен:

g 4x

<2(0 = fu q*О df dr, d<;.

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

Рис. 1. Нагнетательная скважина с пресекающимися плоски -ми трещинами

С целью упрощения дальнейших рассуждений предположим, что трещины протяженностью 2xf и шириной д f расположены симметрично относительно оси скважины, имеют одинаковую геометрию и идентичные гидродинамические свойства. Расход жидкости Q(t), закачиваемый в скважину, распределен равномерно по мощности пласта (по высоте трещин) с поверхностной плотностью q(x,t). Здесь под q понимается количество жидкости, закачиваемой в единицу мощности пласта с единицы длины трещины за единицу времени. При этом плотность распределенного источника вдоль трещины с учетом малости раскрытия разрыва д определяется соотношением: q*(x,y,z,t) ~ q(x,t)/д f (при q* = 0 вне трещины, |y| > df /2).

В силу линейности рассматриваемой задачи и идентичности трещин создаваемое ими возмущение поля давления можно рассматривать как сумму возмущений, вызванных одиночными трещинами, одна из которых расположена, как показано на рисунке 1, вдоль оси x, а другая, повторная трещина, направлена вдоль оси с координатами: feos а — sin а\ /х\ > а ) XV)'

ей

sin а

cosi

При допущении об однородности пласта и равномерности нагнетания по вертикали, общее интегральное представление поля давления (1) после интегрирования по п преобразуется к виду:

-*/

POcyt)=JLf [МЛ nx'y'tJ 4nkJ J

(*-fi2+y2 (*'-?) +У'2

g 4X

d(dr.

Суммарный расход жидкости 0 примет вид:

(¿(Ь) = 2^1 dx.

~хг

Гидродинамика одиночной трещины

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

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

kf8f d2Pf

q; 0 < x < Xf, t > 0.

II дх2

По предположению однородная трещина состоит из двух разнонаправленных симметричных относительно скважины «крыльев», что позволяет рассматривать область фильтрации жидкости только в одном направлении разрыва (например, в интервале 0 < х < х), как показано на рисунке 2.

Рис. 2. Геометрическая модель крыла трещины

На забое скважины задано возмущение давления жидкости Р^ относительно пластового, а приток с торца считается пренебрежимо малым. Таким образом, уравнение (3) дополняется краевыми условиями:

_ дРГ Pf{x, t)\x=0 — Pw,

= 0.

X=Xf

решение сформулированной выше задачи, описывающий процесс фильтрации в трещине, может быть представлено в виде:

* */

Pfix.t)

Pw kf8f

f fq(f,t)df + x J q(f,t)df

0 < x < xf.

Соотношение (4) связывает распределение давления вдоль трещины с плотностью расхода жидкости q(x,t), и должно рассматриваться совместно с уравнением (2) при у = 0:

с хг

т)

Ш

0 -Xf

fr'-fl+y'2

е 4*

(x-fí2

dr,

где 0 < x < x, t > 0, x' = x cosa и y' = x sina.

Для определения распределенного притока вдоль трещины, объединяя уравнения (4) и (5), получим интегральное уравнение: t хг

q

р -JLf [Ш w~4nkj J (t-

0 -Xf

T)

(X-O2

+

kfSf

xf

J fq((,t)d( + x j q(Z,t)dt

d%dT+

где 0 < x < x, и t > 0.

Квазистационарный режим скважины с пересекающимися трещинами

Как показывают предварительные оценки при неизменных условиях эксплуатации скважины (постоянное давление на забое Р или постоянный расход О) квазистационарный режим фильтрации в пласте, характеризующийся неизменным во времени профилем относительного распределения притока вдоль трещины (?(*) = q(x,t)/< q >, где < q >= (¿Ю/фИ-х^), устанавливается при к 1/х2е >100.

При этом интегральное уравнение (6) принимает вид:

м<2СОГ 1 (.* + 02\ .

+ Ei

(-V + V + y^ + aLULjrУ

V 4>rt J \ 4*t J

+ Ei

4xt

( (*'-02+У'2М 1 Г ч

хг

+

* J q(0

df

Здесь Ei

i»

Г e~l

J и

du - интегральная показатель-

ная функция, для малых значений аргумента (Абрамовиц, Стиган, 1979):

Е^-х) « 1п(1.781х)-

Соответственно в квазистационарном приближении вместо выражения (7) для больших значений времени t будем иметь:

хг

AhXf j 4nk.

+

v о

1П\22458xt) + Ы { 2.2458 xt ) +

1

+

+ ln[

0

xf

+

Ж / чСО

где X = x cosa и y' = x sina.

Численное решение задачи о распределенном расходе жидкости вдоль трещины при наличии фильтрационного сопротивления

Здесь и в дальнейшем в расчетах принимаются следующие типичные для промысловых условий нагнетания жидкости (воды) значения параметров: давление на забое - Рк = 4.053-106 Па (40 атм), полудлина трещины - х{ = 100 м, ширина трещины - де = 0.2 м, коэффициент пъезопроводности - н = 2.78 м2/с (104 м2/час),

SCIENT1FK AND TECHNICAL JOURNAL

GEORESOURCES

проницаемость пласта - k = 10-13 м2 (100 мД), вязкость флюида - ^ = 0.001 Па-с (1 спз), время эксплуатации скважины - t = 360000 с (100 часов).

результаты вычислительных экспериментов, представленные на рисунке 3 а, показали, что распределение притока вдоль трещины существенно зависит от ее относительного фильтрационного сопротивления kx f/ (kfdf). На конце трещины наблюдается резкое увеличение оттока, что связано с изменением характера течения жидкости с линейного на квазирадиальный.

В свою очередь, при известной плотности закачки q(x,t), учитывающей фильтрационное сопротивление в трещинах, квазистационарное распределение давления в пласте может быть рассчитано на основе асимптотического представления общего уравнения (2):

где x = x cosa - y sina и y' = x sina + y cosa.

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

Рассмотрим распределение притока в пласт при равномерном распределении давления вдоль трещин Р е = Р ^ при отсутствии фильтрационного сопротивления (к е ^ да). Анализ решения обратной задачи о распределении оттока при заданном давлении показывает, что при увеличении угла раскрытия между трещинами градиенты давления в околоскважинной зоне пласта снижаются и уменьшается локальный отток жидкости, как представлено на рисунке 5. А на удаленной концевой части трещины наблюдается обратный эффект - плотность потока возрастает по причине снижения взаимовлияния между трещинами.

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

Продуктивность скважины и псевдоскин-фактор

Обобщая традиционное понятие скин-фактора (Economides et а1., 2002; Морозов, 2016), определим псевдоскин-фактор скважины с пересекающимися трещинами

и

2 30

г

^20

10

а)

10 20 30 4-0 50 60 70 S0 90

X, м

<¿?

3S-353430252624-

б)

10 20 39 40 50 60 70 S0 90

X, м

Рис. 3. Распределение притока в трещине (а) и давления вдоль трещины (б) в окрестности скважины при взаимно перпендикулярном расположении трещин и различных проницаемостях: 1 - к= 5104мД и 2 - к= 104мД

40-,

Сц 2Ь-

а)

''":,';<ЛХХ --г"" -1

у, м loo^^1^^-]^ х, М 200

-200

-200

б)

у, м

200

X, М

Рис. 4. Возмущение пластового давления в окрестности скважины при взаимно перпендикулярном расположении трещин с проницаемостью:

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

ГЕйРЕСУРСЫ

100

50

80

О 70 СЙ №

й

-200

X, м

200

а) 1 0 2 О 3 0 40 0 0 6 0 7 0 8 0 9 О X, м

б)

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

конечной проницаемости при заданной депрессии. Из соотношения (8) при х =0 непосредственно следует:

Рш = -

гадц «

7 пк Ч2.2458>Л/

2кШхг

о

*>

1п\

2.2458XI

где гш - радиус скважины. В результате, продуктивность скважины будет задаваться формулой: <2 2 к1гк 1

5 = - 1п(х//ё). Кроме того, согласно оценкам, полученным в (Чарный, 1963), на основе метода последовательной смены стационарных состояний, выражение л/2.2458я£ можно рассматривать как радиус подвижного текущего кругового контура питания г . В результате приведенное выше соотношение (9) может быть переписано в эквивалентной форме:

1п

\хг/е)'

где псевдоскин-фактор равен: Х1

о

В частном случае одиночной трещины бесконечной проводимости (ц = 1) псевдоскин-фактор в (9) равен

которое совпадает с решением (Raghavan, Joshi, 1993).

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

Повторный гидравлический разрыв пласта позволяет увеличить псевдоскин-фактор при высокой проницаемости трещин до 5-10 %, как показано на рисунке 6а. Эффективность повторного гидроразрыва пласта согласно

Рис. 6. Псевдоскин-фактор при заданном забойном давлении Р^ с проницаемостью (а) - к г = 106 мД (кх г/ (к^З^) ~0,05) и (б) - к г = 104

а)

. П

мД

8С1ЕИТ1 (ТС АЫй ТЕСНМЮЛ.

б)

рисунку 6б возрастает до 10-20 % в случае небольшой проницаемости трещин, что согласуется с результатами численного решения, полученного в (Лихтарев, Пестриков, 2010).

Заключение

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

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

Повторный гидравлический разрыв продуктивного пласта позволяет увеличить псевдоскин-фактор при высокой проницаемости (к/ ~ 106 мД, кх/ (//) ~ 0,05) трещин до 5-10 %. Эффективность повторного гидроразрыва пласта возрастает до 10-20 % в случае небольшой проницаемости (к/~104 мД, кх//) ~5).

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

Благодарности

Автор выражает благодарность научному руководителю профессору А.Н. Саламатину за ценные советы при планировании и реализации исследования и рекомендации по оформлению работы. Также автор благодарит профессора Г.Т. Булгакову за выбор тематики и направления исследований и выражает признательность рецензенту за ценные замечания и предложения, которые способствовали улучшению данной статьи.

литература

Абрамовиц М., Стиган И. (1979). Справочник по специальным функциям. Москва: Наука, с. 58-59.

Карслоу Г., Егер Д. (1964). Теплопроводность твердых тел. Москва: Наука, с. 251-259.

Лихтарев А.В., Пестриков А.В. (2010). Определение коэффициента продуктивности скважины, пересеченной двумя трещинами гидроразрыва пласта разного азимута, на основе математического моделирования. Научно-технический вестник ОАО «НК»Роснефть», с. 12-14.

Морозов П.Е. (2016). Псевдоскин-фактор и оптимальная проводимость вертикальной трещины гидавлического разрыва пласта. Межд. научно-практ. конф.: Инновации в разведке и разработке нефтяных и газовых месторождений, Казань, с. 53-56.

Тихонов А.Н., Самарский А.А. (1999). Уравнения математической физики. Москва: МГу, с. 477-485.

Чарный И.А. (1963). Подземная гидрогазодинамика. Москва, Гостоптехиздат, с. 250-313.

Cinco-Ley H., Samaniego V.F., Domínguez A.N. (1978). Transient pressure behavior for a well with a finite-conductivity vertical fracture. SPE Journal, 18(4), pp. 253-264.

Economides M., Oligney R., Valkó P. (2002). Unified fracture design: bridging the gap between theory and practice. Alvin, Texas: Orsa Press, 25 p.

Meehan D. N., Horne R.N., Ramey H.J. (1989). Interference testing of finite conductivity hydraulically fractured wells. SPE Annual Technical Conference and Exhibition. Society of Petroleum Engineers, pp. 137-152.

Prats M. (1961). Effect of vertical fractures on reservoir behavior-incompressible fluid case. SPE Journal, 1(02), pp. 105-118.

Raghavan R., Joshi S.D. (1993). Productivity of multiple drainholes or fractured horizontal wells. SPE Form. Eval., 8(01), pp. 11-16.

сведения об авторе

Юлия Ильшатовна Яхина - младший инженер-программист в инженерно-аналитической группе, отдел научных исследований, ООО «Термосим»

Россия, 420108, Казань, ул. Магистральная, 59/1

E-mail: vib-210@mail.ru

Статья поступила в редакцию 29.01.2018;

Принята к публикации 28.04.2018;

Опубликована 30.06.2018

^ш in English

Hydraulic fracturing efficiency evaluation in the vicinity of a single well for a reservoir with two fractures

Yu.I. Yakhina

Termosim LLC, Kazan, Russian Federation

Abstract. The solution of fluid flow problem in an unbounded homogeneous reservoir penetrated by a single well with two intersecting fractures with taking into account their hydrodynamic resistance is constructed and investigated. A general integral representation of the perturbed pressure field is obtained using the instantaneous point-source. As a particular case, the quasi-stationary operating mode of

the well is considered. The accuracy and limitations of the obtained asymptotic solutions are estimated. A comparative analysis obtained results is done. Well productivity and the pseudoskin factor for the fractures system are determined, explicit analytical expressions for these characteristics are constructed. In the course of computational experiments, the interaction of intersecting fractures at different opening

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

angles are investigated. An estimate of the efficiency of repeated hydraulic fracturing of the productive formation is obtained. It is shown that the maximum flow is achieved for a perpendicular arrangement of the fractures, and the distribution of the outflow (inflow) along the flat vertical fracture essentially depends on its relative filtration resistance.

Keywords: hydraulic fracturing, crack azimuth of repeated fracturing, efficiency of repeated fracturing, quasi-stationary operation of well, pseudoskin factor of repeated fracturing, filter resistance in cracks, instantaneous point-source, fracture hydrodynamics

Recommended citation: Yakhina Yu.I. (2018). Hydraulic fracturing efficiency evaluation in the vicinity of a single well for a reservoir with two fractures. Georesursy = Georesources, 20(2), pp. 108-114. DOI: https://doi.org/10.18599/grs.2018.2.108-114

Acknowledgements

The author expresses gratitude to the scientific adviser, Professor A.N. Salamatin for valuable advice in the planning and implementation of research and recommendations on the formatting of work. The author thanks Professor G.T. Bulgakova for choosing the topic and area of research and expresses gratitude to the reviewer for valuable comments and suggestions that contributed to the improvement of this article.

References

Abramowitz M., Stegun A.

Carslaw H.S, Jaeger J.C.

Charnyy I. A.

Cinco-Ley H., Samaniego V.F., Domínguez A.N. (1978). Transient pressure behavior for a well with a finite-conductivity vertical fracture. SPE Journal, 18(4), pp. 253-264.

Economides M., Oligney R., Valkó P. (2002). Unified fracture design: bridging the gap between theory and practice. Alvin, Texas: Orsa Press, 25 p.

Lihtarev A. V., Pestrikov A. V.

Meehan D. N., Horne R.N., Ramey H.J. (1989). Interference testing of finite conductivity hydraulically fractured wells. SPE Annual Technical Conference and Exhibition. Society of Petroleum Engineers, pp. 137-152.

Morozov P.E.

Prats M. (1961). Effect of vertical fractures on reservoir behavior-incompressible fluid case. SPE Journal, 1(02), pp. 105-118.

Raghavan R., Joshi S.D. (1993). Productivity of multiple drainholes or fractured horizontal wells. SPE Form. Eval., 8(01), pp. 11-16.

Tikhonov A. N., Samarskii A. A.

About the Author

Yulia I. Yakhina - Junior Software Engineer, Engineering-Analytical Group, Department of Scientific Research, Termosim LLC

59/1 Magistralnaya str., Kazan, 420108, Russian Federation E-mail: vib-210@mail.ru

Manuscript received 29 January 2018;

Accepted 28 April 2018;

Published 30 June 2018

SCIENT1FK AND TECHNICAL JOURNAL

GEDRESDURCES