Закономерности образования осадков при нейтрализации фторсодержащих кислотных составов после их взаимодействия с породой терригенного коллектора Текст научной статьи по специальности «Химические науки»

УДК 622.276.63

Л.А. Магадова1, e-mail: [email protected]; З.Р. Давлетов1, e-mail: [email protected]; М.М. Мухин1, e-maii: [email protected]

1 Российский государственный университет нефти и газа (Национальный исследовательский университет) им. И.М. Губкина (Москва, Россия).

Закономерности образования осадков при нейтрализации фторсодержащих кислотных составов после их взаимодействия с породой терригенного коллектора

Согласно проведенному литературному обзору сложность химизма протекающих реакций и отсутствие надежных методик лабораторной оценки обусловливают недостаточность сведений в области осадкообразования при кислотных обработках терригенных пластов, в связи с чем в настоящей статье представлены результаты изучения процессов осаждения продуктов реакций при нейтрализации растворов кислот после растворения породы терригенного коллектора. В ходе экспериментальной работы установлено, что количество осадка значительно зависит от концентрации фторсодержащего реагента и степени нейтрализации растворов кислот. Наряду с этим время обработки породы влияет на процесс осадкообразования в несколько меньшей степени, что обусловлено низкими скоростями реакций алюмосиликатов с фторсодержащими кислотными составами. Методом рентгенофазового анализа в структуре осадков идентифицированы малладрит (гексафторсиликат натрия) и криолит (гексафторалюминат натрия). На основании расчета отношений ин-тенсивностей некоторых характеристических рефлексов указанных фаз показано значительное снижение концентрации гексафторсиликата натрия и одновременное увеличение содержания гексафторалюмината натрия при повышении рН растворов кислот. При нейтрализации растворов кислот до рН = 7 в структуре осадка полностью исчезает фаза гексафторсиликата натрия. Установлено, что по мере нейтрализации степень кристалличности полученных осадков постепенно уменьшается, что, вероятно, вызвано протеканием гидролиза гексафторсиликата с образованием аморфной фазы гидратированных форм кремнезема. Результаты инфракрасной (ИК) спектроскопии нарушенного полного внутреннего отражения подтверждают присутствие в структуре осадка кремнекислородных соединений, при этом расчет отношений максимумов поглощения для связей Si-O и Si-F указывает на повышение доли гидратированных форм кремнезема при повышении значений рН растворов.

Ключевые слова: кислотная обработка, кислотный состав, терригенный коллектор, осадкообразование, нейтрализация, рентгенофазовый анализ, ИК-спектроскопия, нарушенное полное внутреннее отражение.

L.A. Magadova1, e-mail: [email protected]; Z.R. Davletov1, e-mail: [email protected]; M.M. Mukhin1, e-mail: [email protected]

1 Gubkin Russian State Oil and Gas University (National Research University) (Moscow, Russia).

Laws of deposits formation in the neutralization of fluorine-containing acid compounds after their interaction with the reservoir rock sandstone

According to a literature review, the complexity of the chemistry of the reactions and the absence of reliable methods of laboratory evaluation cause information lack for sedimentation area in acidizing of terrigenous reservoirs, therefore this article presents the results of studying the processes of reaction products deposition in the neutralization of acids solutions after reservoir rock sandstone dissolving. During the experimental work the significant dependence of sediment amount on the concentration of fluorine-containing reagent solution and pH of the neutralized acids was established. At the same time, the rock processing time affects sedimentation process to a lesser extent, due to the slow speeds of aluminosilicates reactions with fluorinated acid compositions. X-ray diffraction in the structure of deposits identified malladrit (sodium hexafluorosilicate) and cryolite (sodium hexafluoroaluminate). On the basis of

OIL AND GAS PRODUCTION

the calculation of the intensity ratio for certain characteristic reflexes of specified phases a significant decrease of sodium hexafluorosilicate concentration was showed with the simultaneous increase of odium hexafluoroaluminate content at acid solutions pH increase. When neutralizing acid solutions to a pH = 7 the sodium hexafluorosilicate phase completely disappears from the sediment structure. It is found that while neutralization the degree of obtained by neutralization deposits crystallinity gradually decreases, probably due to hexafluorosilicate hydrolysis with a formation of amorphous phase of hydrated form of silica. The results of the infrared (IR) spectroscopy of attenuated total internal reflection confirm the presence of silicon-oxygen compounds in deposit structure, and the calculation of the relationship of the maximum absorption for Si-O and Si-F bonds indicates the increase of hydrated silica forms rate at higher neutralization pH.

Keywords: acidizing, acid composition, sandstone reservoirs, sedimentation, neutralization, X-ray diffraction, infrared spectroscopy, attenuated total internal reflection.

Взаимодействие фторсодержащих кислотных составов с минералами породы терригенных коллекторов сопряжено с протеканием интенсивных процессов осадкообразования, что способствует закупориванию пор,значительно ухудшая фильтрационно-емкостные свойства пласта [1]. При рассмотрении взаимодействия плавиковой кислоты с алюмосиликатами зачастую выделяют первичные, вторичные и третичные реакции, представленные схемами (1-3), соответственно [2-5].

HF + M-AL-Si + H+ -> H2SiF6 + ALF (3-x) + M+,

2 6 x '

H2SiF6 + M-AL-Si + H+ + H2O ->

2 6 2

ALFx(3-x) + M+ + nSiO2 mH2O,

(1)

(2)

А^х<3-х> + М-А1^1 + Н+ + Н20

_> А^у<3-У> + М+ + nSiO2 тн20, х > у, (3)

где М-А1^' - алюмосиликат, М+ - катион металла.

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

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

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

оо 7 ю

CVJ

б

Ol CVÎ CM

X з

-------------1-------------1------------- 1 1 1 1 1 1 ------------- \

у = 0,543x2 - 3,765x + 6,775 R2 = 0,967

1 1

1 1 У 1 1 / k____________

_______________________________-yi-A 1 1 ^r

1 1 1 JÈf ___________4_______/^Г____________ _____________

T 1 —-4— 1 -

5

pH

Рис. 1. Зависимость отношения интенсивностей рефлекса при 20 = 22,9° для гексафторалюмината натрия и рефлекса при 20 = 26,8° для гексафторсиликата натрия от рН растворов

Fig. 1. Dependence of the reflex intensity ratios at 20 = 22.9° for sodium hexafluoroaluminate and reflex at 20 = 26.8° for sodium hexafluorosilicate on neutralization pH of acid solutions

Ссылка для цитирования (for citation):

Магадова Л.А., Давлетов З.Р., Мухин М.М. Закономерности образования осадков при нейтрализации фторсодержащих кислотных составов после их взаимодействия с породой терригенного коллектора // Территория «НЕФТЕГАЗ». 2016. № 10. C. 54-59.

Magadova L.A., DavLetov Z.R., Mukhin М.М. Laws of deposits formation in the neutralization of fluorine-containing acid compounds after their interaction with the reservoir rock sandstone (In Russ.). Territorija «NEFTEGAZ» = OiL and Gas Territory, 2016, No. 10, P. 54-59.

TERRITORIJA NEFTEGAS - OIL AND GAS TERRITORY No. 10 october 2016

55

Массовые доли осадков, полученных при повышении рН кислотных составов Mass fractions of deposits obtained by increasing the pH of acidic compounds

№ Item Содержание, % масс. Content, wt. % Время обработки, ч Treatment time, h Массовая доля осадка, % масс. Mass fraction of deposit , wt. %

nh4f HCl рН = 3 рН = 4 рН = 5 рН = 6 рН = 7

1 1,85 8,0 1,0 2,7 5,4 8,3 13,9 18,7

2 3,70 10,0 1,0 7,1 11,9 15,5 22,7 27,6

3 5,55 12,0 1,0 12,2 19,6 25,4 31,9 36,2

4 1,85 8,0 4,0 4,5 8,2 11,4 16,7 21,8

5 3,70 10,0 4,0 10,4 15,8 19,7 27,2 31,3

6 5,55 12,0 4,0 16,2 23,5 29,2 36,1 39,9

7 1,85 8,0 8,0 6,9 11,8 14,1 18,4 23,2

8 3,70 10,0 8,0 14,1 19,9 23,3 27,5 32,6

9 5,55 12,0 8,0 20,7 27,6 33,1 39,3 42,8

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

Развитие данного метода оценки осадкообразования представлено в работе [7] в виде комплексной программы иссле-

дований, включающей три последовательных раздела:

• оценка вторичных осадков при взаимодействии раствора кислотного состава с образцом дезинтегрированного керна в свободном объеме;

• физическое моделирование кислотной обработки на моделях фильтрации в условиях, приближенных к пластовым;

• литолого-минералогическая характеристика образцов кернов до и после проведения фильтрационных экспериментов.

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

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

80

70

Ъ 60

50

40

------------- ------------- ------------- -------------

1— ------------- ------------- -------------

у = -0,679х2 R2 f 0,236х + 77,043 = 0,995

5

pH

Рис. 2. Зависимость степени кристалличности осадка от рН кислотного состава 5,6 % NH4F/12 % HCL

Fig. 2. Dependence of the deposit degree of crystaLLinity on neutralization pH of the acid composition of 5.6 % NH4F/12 % HCL

OIL AND GAS PRODUCTION

экспериментальная часть

В работе использовался измельченный керновый материал терригенного коллектора следующего минералогического состава (% масс.): кварц - 60, полевые шпаты - 32, глинистые минералы -5, слюды - 3, пирит - следы. В качестве кислотных составов исследованы системы NH4F/HCL и HF/HCL, содержащие эквивалентные количества фтора. Для определения минералогического состава породы и идентификации структуры осадков использовали порошковый рентгеновский дифрактометр Thermo Scientific ARL X'TRA (источник излучения - рентгеновская трубка с медным анодом). Дифрактограммы регистрировали в интервале углов от 4 до 60° 20. Качественный и количественный рентгенофазовый анализ проводили с помощью программ CrystaLLographica Search-Match и Siroquant, соответственно. Степень кристалличности рассчитывали по полученным дифрактограммам с использованием программного пакета CrystaLLinity.

ИК-спектры осадков снимали с помощью инфракрасного Фурье-спектрометра Thermo Scientific NicoLet iS10 с приставкой однократного нарушенного полного внутреннего отражения Smart iTR (кристалл ZnSe). Спектральный диапазон съемки составлял 4000...650 см-1. Кислотную обработку кернового материала проводили при 80 °С в течение 1; 4 и 8 часов. По истечении времени реакции кислоту отфильтровывали. Далее порцию отработанной кислоты нейтрализовывали до значений рН в диапазоне 3-7. Повышение рН осуществляли путем добавления к прореагировавшим кислотным составам раствора щелочи в соответствии с методическими подходами, описанными в работах [6, 7]. Полученный осадок переносили на бумажный фильтр, промывали дистиллированной водой, сушили до постоянной массы. Массовую долю осадка определяли путем расчета отношения массы полученного осадка к исходной навеске породы с учетом объема порции кислоты, взятой для анализа.

обсуждение результатов

Массовые доли полученных осадков в зависимости от концентрации реагентов, продолжительности обработки и рН растворов на примере составов NH4F/ HCl представлены в таблице. Было обнаружено, что интенсивность осадкообразования в значительной степени зависит от концентрации фтор-содержащего реагента и степени нейтрализации растворов кислот. Наряду с этим время обработки породы влияет на данный процесс в несколько меньшей степени, что обусловлено низкими скоростями реакций алюмосиликатов с фторсодержащими кислотными составами.

В структуре осадка методом рентге-нофазового анализа идентифицированы малладрит (гексафторсиликат натрия) Na2[SiF6] и криолит (гекса-фторалюминат натрия) Na3[AlF6]. Для оценки количественного содержания указанных фаз были рассчитаны отношения интенсивностей некоторых рефлексов гексафторалюмината и

ООО «ЛИTAC» ООО «ЛИТАС РЕНТГЕН» ООО «ЛИТАС СЕРВИС»

ООО «АСК-РЕНТГЕН-ИДЕЛЬ» web: www.litas-x-ray.ru~,www.litas.ru

г. КАЗАНЬ _

Я + 7(843) 564-65-51,

564-48-37, 557-34-07 ПРОИЗВОДСТВ О, e-mail: [email protected], [email protected], П D О Д 3 Ж 3

[email protected] f « *

Сервис.

на правах рекламы

ДОБЫЧА НЕФТИ И ГАЗА

0,6

0,5

0,4

I I" 0,3

F

0,2

0,1

0

1280

-------- -------- ________ -------- ________ -------- -------- ________ -------- _____ / \

Si-O Si-F

1180

1080

980

880

780

680

Волновое число, см1 Wave number, cm-1

Нейтрализация до рН = 3 ™ Нейтрализация до рН = 5 ™ Нейтрализация до рН = 7 Neutralization to рН = 3 Neutralization to рН = 5 Neutralization to рН = 7

Нейтрализация до рН = 4 ~ Нейтрализация до рН = б Neutralization to рН = 4 Neutralization to рН = 6

Рис. 3. ИК-спектры поглощения осадков (НПВО), полученных при различных значениях рН растворов

Fig. 3. IR absorption spectra of deposits (ATR) obtained at different values of neutralization pH

Рис. 4. Зависимость отношения максимумов поглощения на ИК-спектре для связей Si-O и Si-F в осадке от рН растворов

Fig. 4. Dependence of the absorption maximum ratio in the IR spectrum for the Si-O and Si-F bonds in the deposit on neutralization pH

гексафторсиликата натрия. В случае Na3[ALF6] был выбран рефлекс при двойном угле рассеяния 20, равном 22,9°, а для Ма2[Б1Т6] - при 20 = 26,8°. Данные рефлексы характеризуются высокой интенсивностью и не перекрываются рефлексами второй фазы. Расчет указанного отношения (рис. 1) свидетельствует о значительном снижении содержания Иа2[Б1Т6] и одновременном увеличении концентрации Ыа3[Д1Р6] в структуре осадка. Присутствие кристаллической фазы Иа2[Б1Т6], по всей видимости, обусловлено осаждением гексафторсиликат-анионов (4), образовавшихся при растворении алюмосиликатов, под действием гидроксида натрия [8]:

Н2[Б1Т6] + 2Иа0Н

- Иа^и + 2Н20. (4)

Последующее уменьшение содержания гексафторсиликата натрия в структуре осадка по мере повышения рН раствора (5) объясняется его гидролитической неустойчивостью [9]:

Иа2[Б1Т6] + 4Иа0Н ->

6Иа16 + БЮ2.2Н20. (5)

По данным рентгенофазового анализа, при рН = 7 фаза Иа2[Б1Т6] в составе осадка полностью исчезает. Согласно полученным результатам при повышении рН кислотных составов происходит существенное снижение степени кристалличности осадка (рис. 2), что, вероятно, обусловлено интенсивным осаждением гидратированных форм кремнезема, имеющих аморфную структуру.

Для подтверждения результатов рентгенофазового анализа и идентификации соединений аморфной фазы образцы осадков были изучены с помощью ИК-спектроскопии НПВО. Характеристические максимумы на ИК-спектрах поглощения осадков (рис. 3) были соотнесены с поглощением связей Б1-0 и Б1-Р [10-12].

Расчет отношений интенсивностей поглощения связей Б1-0 и Б1^ при повышении степени нейтрализации кислотных составов свидетельствует об увеличении содержания гидратированных форм

58

№ 10 октябрь 2016 ТЕРРИТОРИЯ НЕФТЕГАЗ

OIL AND GAS PRODUCTION

кремнезема и одновременном снижении количества гексафторсиликата натрия (рис. 4), подтверждая протекание разложения Na2[SiF6] с образованием SiO2.nH2O.

В то же время результаты ИК-спек-троскопии указывают на присутствие анионов Р1Т6]2- в структуре всех исследованных осадков. Отсутствие Ма2Р1Т6] при рН = 7, по данным рентгенофазо-вого анализа, может быть объяснено образованием гидратированных форм кремнезема на поверхности кристаллов гексафторсиликата натрия в ходе

реакции гидролиза [13], что приводит к «рентгеноаморфности» фазы Na2[SiF6] и затрудняет ее идентификацию с помощью рентгеновской дифракции. В случае системы HF/HCl были получены аналогичные результаты и закономерности, что и для кислотного состава NH4F/HCl.

выводы

1. Количество образующегося осадка при нейтрализации кислотных составов после обработки породы терригенно-го коллектора значительно зависит от

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

2. По мере нейтрализации кислотных составов наблюдается осаждение гексафторсиликата и гексафторалюмината натрия, а также гидролиз гексафторсиликата с образованием аморфной фазы гидратированных форм кремнезема.

Литература:

1. Магадова Л.А. и др. Осадкообразование при взаимодействии кислотных составов с минералами терригенного коллектора // Нефтепромысловое дело. 2015. № 9. С. 31-36.

2. Gdanski R.D. Kinetics of the Primary Reaction of HF on ALumino-SiLicates. SPE Production & Facilities, November, 2000. P. 279-287.

3. Gdanski R.D. Kinetics of the Secondary Reaction of HF on ALumino-SiLicates. SPE Production & FaciLities, November, 1999. P. 260-268.

4. Shuchart C.E., Buster D.C. Determination of the Chemistry of HF Acidizing with the Use of 19F NMR Spectroscopy. SPE Conference Paper 28975-MS, 1995.

5. Gdanski R.D. Kinetics of the Tertiary Reactions of HydrofLuoric Acid on ALuminosiLicates. SPE Production & FaciLities, May, 1998. P. 75-80.

6. Мазаев В.В., Томчук Н.Н., Лавренова Н.А. Адаптация кислотных составов для обработки призабойной зоны пласта скважин месторождений Западной Сибири // Пути реализации нефтегазового потенциала ХМАО - Югры. Ханты-Мансийск: ИздатНаукаСервис, 2006. Т. 1. С. 401-411. (Цит. по: Глущенко В.Н., Силин М.А. Нефтепромысловая химия: Изд. в 5 т. Т. 4. Кислотная обработка скважин / Под ред. И.Т. Мищенко. М.: Интерконтакт Наука, 2010. 703 с.)

7. Макиенко В.В. и др. Разработка методики по оценке вторичного осадкообразования при кислотных обработках терригенных коллекторов // Мат-лы I Международной (IX Всероссийской) науч.-практ. конф. «Нефтепромысловая химия». М.: Изд. центр РГУ нефти и газа им. И.М. Губкина, 2014. С. 31-33.

8. Лидин Р.А., Молочко В.А., Андреева Л.Л. Химические свойства неорганических веществ: учеб. пособие для вузов / Под ред. Р.А. Лидина. 3-е изд., испр. М.: Химия, 2000. 480 с.

9. Мышляева Л.В., Краснощеков В.В. Аналитическая химия кремния. М.: Наука, 1972. 212 с.

10. MiLLer F.A., WiLkins C.H. Infrared Spectra and Characteristic Frequencies of Inorganic Ions. AnaLyticaL Chemistry, 1952, VoL. 24, No. 8. P. 1253-1294.

11. Nakamoto K. Infrared and Raman Spectra of Inorganic and Coordination Compounds. Part A: Theory and AppLications in Inorganic Chemistry, 6th Edition, John WiLey & Sons, Inc., 2009. 432 pp.

12. Плюснина И.И. Инфракрасные спектры минералов. М.: Изд-во Моск. ун-та, 1976. 175 с.

13. Ахметов Т.Г. и др. Химическая технология неорганических веществ: учеб. пособие. В 2 кн. Кн. 2. М.: Высш. шк., 2002. 533 с.

References:

1. Magadova L.A. et. al. Sedimentation in the interaction of acidic compounds with minerals sandstone reservoirs. Neftepromyslovoe delo = Petroleum Engineering, 2015, No. 9, P. 31-36. (In Russian)

2. Gdanski R.D. Kinetics of the Primary Reaction of HF on Alumino-Silicates. SPE Production & Facilities, November, 2000. P. 279-287.

3. Gdanski R.D. Kinetics of the Secondary Reaction of HF on Alumino-Silicates. SPE Production & Facilities, November, 1999. P. 260-268.

4. Shuchart C.E., Buster D.C. Determination of the Chemistry of HF Acidizing with the Use of 19F NMR Spectroscopy. SPE Conference Paper 28975-MS, 1995.

5. Gdanski R.D. Kinetics of the Tertiary Reactions of Hydrofluoric Acid on Aluminosilicates. SPE Production & Facilities, May, 1998. P. 75-80.

6. Mazaev V.V., Tomchuk N.N., Lavrenova N.A. Adaptation of acid compositions for the treatment of wells bottomhole zone in Western Siberia fields formation. Puti realizatsii neftegazovogo potentsiala KhMAO - Yugry = Ways of realization of oil and gas potential KhMAD - Ugra, Khanty-Mansiysk, IzdatNaukaServis, 2006, Vol. 1, P. 401-411 (Cit. Glushenko V.N., Silin M.A. Oilfield chemistry, Ed. 5 vols, Vol. 4, Wells acidizing; ed. by I.T. Mischenko, Moscow, Intercontact Nauka, 2010, 703 pp.) (In Russian)

7. Makienko V.V. et. al. Development of techniques for the evaluation of secondary sedimentation in terrigenous reservoirs acidizing. Proceedings of the I International (IX All-Russia) Scientific-practical conference "Oilfield Chemistry". Moscow, Publishing Center of Gubkin Russian State University of Oil and Gas, 2014, P. 31-33. (In Russian)

8. Lidin R.A., Molochko V.A., Andreeva L.L. Chemical properties of inorganic substances: Manual for higher schools. Ed. by R.A. Lidina; 3rd rev., revised, Moscow, Chemistry, 2000, 480 pp. (In Russian)

9. Myshlyaeva L.V., Krasnoshchekov V.V. Analytical chemistry of silicon. Moscow, Nauka, 1972, 212 pp. (In Russian)

10. Miller F.A., Wilkins C.H. Infrared Spectra and Characteristic Frequencies of Inorganic Ions. Analytical Chemistry, 1952, Vol. 24, No. 8. P. 1253-1294.

11. Nakamoto K. Infrared and Raman Spectra of Inorganic and Coordination Compounds. Part A: Theory and Applications in Inorganic Chemistry, 6th Edition, John Wiley & Sons, Inc., 2009. 432 pp.

12. Plyusnina I.I. Infrared spectra of minerals. Moscow University Publishing House, 1976, 175 pp. (In Russian)

13. Akhmetov T.G. et. al. Chemical technology of inorganic substances, Manual in 2 vols, Vol. 2, Moscow, Vysshaya shkola, 2002, 533 pp. (In Russian)

TERRITORIJA NEFTEGAS - OIL AND GAS TERRITORY No. 10 october 2016

59