книги / Инструкция по созданию наблюдательных станций и производству инструментальных наблюдений за процессами сдвижения земной поверхности при разработке нефтяных месторождений в регионе
..pdf
- 21 -
fдоп(план) = 2 (m1,план)2 + (m2,план)2 + (m3,план)2 + .... + (mn,план)2 , fдоп(h) = 2 (m1,h )2 + (m2,h )2 + (m3,h )2 + .... + (mn,h )2 ,
где mi,план и mi,h – погрешности измерения сторон (векторов) замкнутого полигона в
плане и высоте, мм; которые вычисляются по формулам, если иного не предусмотрено проектом:
где D – длина вектора, км.
Контроль качества полевых измерений (тип
решения векторов, невязки в 
полигонах, контрольные
измерения)
Уточненное решение векторов спутниковой сети
для каждого периода наблюдений.
Оценка внутренней сходимости, выявление
грубых измерений, получение уравненных геометрических связей между пунктами сети.
Совокупность устойчивых опорных
пунктов.
Каталог координат пунктов спутниковой сети данного периода и их оценки.
Смещения пунктов и их 
оценки.
mi,план = (5 + 0, 5D) , mi,h = (10 + 0, 5D) ,
1. Предварительная (полевая) обработка
спутниковых измерений.
2. Заключительная (камеральная) обработка спутниковых измерений.
3. Предварительное (свободное) уравнивание спутниковой сети каждого
периода наблюдений.
4. Анализ устойчивости опорных пунктов.
5. Заключительное уравнивание спутниковых измерений относительно
координат устойчивых опорных пунктов.
6. Вычисление рассогласований в координатах рабочих пунктов между циклами
наблюдений.
Ежедневные файлы спутниковых наблюдений, данные полевых журналов.
Координаты исходных пунктов, точные эфемериды, параметры атмосферы и др. данные для улучшения точности спутниковых измерений.
Выбор стратегии назначения весов спутниковых измерений и других статистических параметров.
Свободно-уравненные спутниковые сети сравниваемых периодов.
Координаты опорных пунктов.
Каталоги координат пунктов спутниковой сети разных периодов.
Графики, карты-схемы |
|
7. Графическое |
|
|
|
изображение процесса смещений |
|
Картографический |
|
вертикальных и |
|
|
||
|
земной поверхности на |
|
материал, геологические |
|
горизонтальных смещений |
|
|
||
|
наблюдаемом объекте за |
|
разрезы, параметры |
|
земной поверхности на |
|
|
||
|
данный период и его |
|
разработки объектов и т.п. |
|
наблюдаемом объекте. |
|
|
||
|
интерпретация. |
|
|
|
|
|
|
|
Рис.6.2. Блок-схема обработки циклов спутниковых измерений
- 22 -
Расхождения между результатами определения векторов из разных сеансов устанавливаются, исходя из допуска на разность двойных измерений по правилу:
–доп = 2m 2 – для парных измерений,
–доп = 2m – для уклонения от средних значений при числе определений
базисной линии больше 2,
где m – погрешность измерения вектора в плане или высоте, вычисляемая по формулам, приведенным выше, если иного не предусмотрено проектом.
Расхождения, превышающие допустимые, не принимаются. В этом случае необходимо выполнить повторные наблюдения.
6.3.5.Заключительная обработка спутниковых наблюдений может выполняться с использованием стандартного программного обеспечения фирмыизготовителя спутниковых приемников, сертифицированного для применения на территории Российской Федерации и обеспечивающего необходимую точность мониторинга. Рекомендуется выполнять обработку в специальных (профессиональных) программах от независимых участников, например: Bernese, GIPSY, GAMIT.
6.3.6.Обработка спутниковых измерений выполняется методом отдельных независимых линий и заключается в вычислении приращений координат каждого определяемого пункта по отношению к исходному в соответствии с геометрическим построением сети.
6.3.7.В качестве исходного используется пункт, имеющий полные значения координат в геоцентрической системе с точностью не более 1 м. Большие ошибки
вкоординатах местоположения исходного пункта приводят к ошибкам вычислений базисных линий, ошибкам масштаба, элементов разворота и могут сказываться на неудачном разрешении неоднозначности.
6.3.8.Циклы спутниковых наблюдений должны обрабатываться в одной системе геоцентрических координат.
6.3.9.При обработке спутниковых измерений необходимо придерживаться следующих принципов:
−алгоритм разрешения неоднозначности должен быть основан на методе использования фазовых измерений несущих частот или их комбинации;
−угол маски должен быть установлен на значение 15°;
−применять модели стандартной атмосферы для минимизации тропосферной и стратосферной рефракций;
−использовать двухчастотные наблюдения для исключения влияния ионосферы;
−использовать параметры калибровок приемных антенн;
−при вычислении векторов, длиной свыше 10 км, использовать уточненные (апостериорные) эфемериды спутников;
−рекомендуется проводить обработку сетей разных циклов с одними и теми же параметрами обработки векторов.
6.3.10. Результаты заключительной обработки спутниковых измерений должны удовлетворять критериям п. 6.3.4.
-23 -
6.3.11.Уравнивание спутниковых измерений выполняется строгим способом по принципу наименьших квадратов, т.е. минимизацией выражения pv2 , где p –
веса измерений, v – соответствующие к ним поправки.
6.3.12.Предварительное уравнивание выполняется как свободное, т.е. либо без закрепления координат пунктов, либо с фиксацией положения исходного пункта. Целью свободного уравнивания является устранение неоднозначности и увеличения точности спутниковых определений, а также получение уравненных геометрических связей между пунктами сети. Итогом свободного уравнивания является каталог координат пунктов для дальнейшего анализа. Рекомендуется получать каталог координат пунктов сети относительно координат исходного пункта. Для первой серии наблюдений свободное уравнивание является заключительным, а полученный каталог координат пунктов окончательным.
6.3.13.Оценка устойчивости опорных пунктов выполняется путем анализа изменения их взаимного положения между текущей и начальной сериями наблюдений по результатам свободного уравнивания спутниковых измерений соответствующих периодов. Наиболее устойчивым из опорных является пункт, имеющий наименьшую сумму квадратов равновесных смещений в рассматриваемый период. Изменения положений других опорных пунктов относительно устойчивого не должны превышать утроенного значения погрешностей определения данных изменений.
6.3.14.Рекомендуется по результатам первых двух серий наблюдений получить наиболее вероятные приращения координат между опорными пунктами, т.е. уточнить их положение при условии устойчивости последних. В этом случае оценку устойчивости в каждой последующей серии наблюдений необходимо выполнять относительно вероятнейшего положения опорных пунктов.
6.3.15.Заключительное уравнивание спутниковых измерений сравниваемых периодов выполняется строгим способом по принципу наименьших квадратов относительно взаимных положений устойчивых опорных пунктов. В качестве фиксируемых взаимных положений устойчивых опорных пунктов используют приращения координат, определенные по данным свободного уравнивания первой серии наблюдений, или их вероятнейшие значения, определенные по результатам первых двух серий наблюдений (см. п. 6.3.14). В ходе заключительного уравнивания получают: пространственные прямоугольные (X,Y,Z) и эллипсоидальные (B,L,H) координаты рабочих пунктов; средние квадратические погрешности координат рабочих пунктов, определенные через полную ковариационную матрицу, получаемую в процессе уравнивания; размеры плановых эллипсов ошибок и их ориентацию. Точностные характеристики результатов заключительного уравнивания должны соответствовать требованиям проекта мониторинга и данной инструкции.
6.3.16.Расчет смещений рабочих пунктов рекомендуется выполнять с использованием пространственной прямоугольной системы координат по следующим формулам:
- 24 -
ξxi0− k |
|
− sin B cos L |
− sin B sin L |
cos B |
||||||
|
|
|
|
|
i |
i |
i |
i |
|
i |
|
0−k |
− sin Li |
|
cos Li |
|
0 |
||||
ξyi |
|
= |
|
|
||||||
|
η |
0−k |
|
|
cos B cos L |
cos B sin L |
sin B |
|||
|
i |
|
|
i |
i |
i |
i |
|
i |
|
|
|
|
|
|
||||||
|
|
|
|
|
ξ0−k = ξx0− k 2 |
+ ξy0−k 2 |
, |
|||
|
|
|
|
|
i |
|
i |
|
i |
|
XikYi kkZi
−Xi0
−Yi0 ,
−Zi0
где X ik ,Yik , Zik и X i0 ,Yi0 , Zi0 – пространственные прямоугольные координаты i-го
пункта, определенные по результатам строгого уравнивания относительно взаимных положений устойчивых опорных пунктов, соответственно в k-м и
начальном цикле наблюдений; Bi , Li – эллипсоидальные координаты i-го пункта (широта, долгота), определенные в начальном цикле наблюдений; ξxi0−k , ξyi0− k – горизонтальные смещения i-го пункта между начальным и k-м циклами наблюдений по соответствующим осям; ξi0− k – значение полного вектора горизонтального смещения i-го пункта между начальным и k-м циклами наблюдений; ηi0−k – вертикальное смещение i-го пункта между начальным и k-м
циклами наблюдений.
Направление полного вектора горизонтального смещения пункта находят из решения обратной геодезической задачи.
Для выполнения корректных вычислений пространственные прямоугольные и эллипсоидальные координаты пунктов должны быть отнесены к одному эллипсоиду.
6.3.17. Средние квадратические погрешности горизонтальных и вертикальных смещений пунктов вычисляются по формулам:
0−k |
) = |
|
0−k |
) |
2 |
|
|
0−k |
) |
|
2 |
, |
|
m(ξi |
m(ξxi |
|
|
+ m(ξyi |
|
|
|||||||
m(ξx0− k ) = (mx0 )2 |
+ (mxk )2 |
, |
|
|
|
||||||||
|
i |
|
|
|
i |
|
|
|
i |
|
|
|
|
m(ξy0− k ) = (my0 )2 |
+ (myk )2 |
, |
|
|
|
||||||||
|
i |
|
|
|
i |
|
|
|
i |
|
|
|
|
m(ηi0− k ) = (mHi0 )2 + (mHik )2 , |
|
|
|
||||||||||
где mxi0 , myi0 , mH i0 и mxik , myik , mH ik |
– |
средние |
квадратические погрешности |
||||||||||
определения координат и высот i-го рабочего пункта, определенные по результатам строгого уравнивания относительно взаимных положений устойчивых опорных пунктов, соответственно в начальном и k-м цикле наблюдений.
Полученные погрешности горизонтальных и вертикальных смещений пунктов должны удовлетворять требованиям проекта мониторинга и данной инструкции.
6.3.18. Результаты вычислений смещений пунктов должны быть представлены в графическом виде за весь срок мониторинга и при необходимости за отдельные периоды. Горизонтальные смещения отображаются на плане наблюдательной станции в виде векторных траекторий с указанием значения полного вектора плановой подвижки. Вертикальные смещения отображаются на плане наблюдательной станции в виде изолиний или столбчатых диаграмм для каждого пункта. При необходимости вертикальные подвижки представляют в виде графиков смещений отдельных пунктов или по заданным направлениям.
- 25 -
Масштаб(ы) отображения смещений выбирают, исходя из информативности и наглядности картины исследуемого процесса.
6.3.19. Схемы (графики) смещений сопровождают изображением распределения их предельных средних квадратических погрешностей, рассчитанных по нижеследующим формулам, если иного не предусмотрено проектом мониторинга:
mξi Пред. = k × m(ξi0−k ) , mηi Пред. = k × m(ηi0−k ) ,
где k – коэффициент кратности интеграла вероятностей нормального закона распределения, устанавливающий ширину доверительного интервала значений измеряемой величины в соответствии с заданной доверительной вероятностью. При доверительной вероятности равной 0.997, k = 3.
7.ТРЕБОВАНИЯ, ПРЕДЪЯВЛЯЕМЫЕ К ПРИБОРАМ, ИНСТРУМЕНТАМ
ИОРГАНИЗАЦИИ-ИСПОЛНИТЕЛЮ РАБОТ
7.1. При выборе аппаратуры для проведения наблюдений следует воспользоваться следующими рекомендациями:
•Для проведения наблюдений рекомендуется применять приборы и инструменты известных фирм-товаропроизводителей.
•Рекомендуется использовать мультисистемные спутниковые приемники, принимающие сигналы, как минимум, двух радионавигационных систем: GPS Navstar и ГЛОНАСС.
•Спутниковые приёмники должны обеспечивать измерение фазы несущего радиосигнала на принимаемых диапазонах частот, ее интеграции и сохранения на внутреннем или внешнем накопителе для последующей обработки.
•Приборы и инструменты должны иметь сертификат Комитета Российской Федерации по стандартизации метрологии и сертификации и быть зарегистрированы в государственном реестре средств измерений.
•Приборы и инструменты должны иметь свидетельство о метрологической поверке, действительной на период выполнения измерений.
7.2. К выполнению работ по составлению проекта маркшейдерскогеодезического мониторинга, выполнению и обработке инструментальных наблюдений допускаются организации, имеющие опыт выполнения подобных работ на территории ВКМКС и действующую лицензию на производство маркшейдерских работ с соответствующим выполняемым работам видом деятельности.
-26 -
8.ТРЕБОВАНИЯ К ТОЧНОСТИ И ПЕРИОДИЧНОСТИ ИНСТРУМЕНТАЛЬНЫХ
НАБЛЮДЕНИЙ
8.1. Для изучения современных движений земной поверхности и техногенных процессов на территориях эксплуатируемых нефтепромыслов, расположенных на территории ВКМКС, точность определения элементов наблюдательной сети должна быть не менее:
•средняя квадратическая погрешность (СКП) определения превышений в нивелирных ходах по профильным линиям должна быть не более ±1 мм/км и не более ±0,3 мм на станции;
•СКП измерения расстояний при изучении оползневых, карстовых процессов и зон активных тектонических разломов должна быть не более
±3 мм;
•СКП измерения векторов в сети с помощью спутниковых радионавигационных систем должна быть не более величин, устанавливаемых п.6.3.4;
•СКП определения планового положения пункта спутниковой сети относительно опорных пунктов должно быть не более ±7 мм;
•СКП определения высотного положения пункта спутниковой сети относительно опорных пунктов должна быть не более ±11 мм.
8.2. СКП определения векторов вертикальных смещений относительно опорных пунктов должна быть не более:
• методом высокоточного нивелирования ±10 мм;
• методом спутниковых измерений ±15 мм;
• при наблюдении за оползневыми и карстовыми процессами ±2 мм;
8.3. СКП определения векторов горизонтальных смещений относительно опорных пунктов должна быть:
•методом спутниковых измерений ±10 мм;
•при наблюдении за оползневыми и карстовыми процессами ±3 мм.
8.4.Начальная серия наблюдений на полигоне должна быть выполнена не менее чем за 6 месяцев до начала пробной эксплуатации месторождения, но не ранее чем за 3 месяца после окончания работ по закладке реперов наблюдательной станции.
8.5.Периодичность наблюдений согласно требованиям нормативного документа «Правила промышленной безопасности при освоении месторождений нефти на площадях залегания калийных солей» (Постановление Госгортехнадзора России от 4.02.2002 г. за № 8) устанавливается один раз в год.
Вслучае возникновения сейсмических явлений или при интенсивном развитии процессов оседания (максимальное оседание – более 100 мм в год) необходимо выполнять наблюдения два раза в год.
8.6. При обнаружении мест концентрированных деформаций, связанных с наличием тектонических нарушений, при наличии в данных местах инженерных сооружений (нефтепроводы, водоводы, кусты скважин, ДНС, КНС, зоны промышленной и гражданской застройки и т.д.) рекомендуется увеличить количество выполняемых циклов наблюдений в таких местах до двух ежегодно и рассмотреть целесообразность организации наблюдений непосредственно за объектами.
-27 -
8.7.При возникновении в районе месторождения техногенного сейсмического явления рекомендуется, по возможности, в кратчайшие сроки выполнить дополнительный цикл наблюдений.
9.ОБРАБОТКА РЕЗУЛЬТАТОВ НАБЛЮДЕНИЙ
9.1. По результатам нулевого и последующих циклов нивелирования значения оседаний, наклонов и кривизны между циклами наблюдений следует получать по формулам:
|
|
η |
= H n − H n + m ; |
i = |
ηi + 1 − ηi |
; |
K |
= |
ik + 1 − ik |
|
; |
|
|
|
i |
i |
i |
k |
lk |
|
|
(lk + 1 + lk ) / 2 |
|
||
|
|
|
|
|
|
|
|
|
||||
где η |
i |
– оседание i -го репера за период m − n ; |
H n – высотная отметка i -го |
|||||||||
|
|
|
|
|
|
|
|
|
i |
|
|
|
репера из предыдущей |
n -й серии нивелирования; |
H n + m – |
высотная отметка |
|||||||||
|
|
|
|
|
|
|
|
|
|
i |
|
|
i -го репера из последующей серии нивелирования; |
ik – наклон k -го интервала; |
|||||||||||
ηi + 1 – оседание переднего репера; ηi – оседание заднего репера; lk – длина k -го интервала; K – кривизна в общей точке двух смежных интервалов k + 1 и k ;
ik |
– наклон предыдущего интервала; ik + 1 – наклон последующего интервала; |
lk |
– длина предыдущего интервала; lk + 1 – длина последующего интервала. |
9.2.В процессе камеральной обработки оформляются ведомости вычислений средних превышений, оценки стабильности реперов исходной основы, уравнивания и распределения поправок в превышения, вероятнейших отметок и вертикальных смещений реперов с оценкой их надежности. Составляется схема уравнивания нивелирной сети. По значениям вертикальных смещений реперов и их доверительным интервалам строятся соответствующие графики, а при достаточно сгущенной сети – изолинии оседаний, наложенные на план местности.
9.3.Значения и направления векторов горизонтальных сдвижений пунктов спутниковой сети следует определять, руководствуясь п. 6.3.16 данной инструкции.
9.4.Значение горизонтального сдвижения ξ рабочего пункта локальной наблюдательной станций для изучения оползневых процессов и концентрированных горизонтальных деформаций в местах контактов тектонических блоков, заложенной прямолинейно, следует определять как разницу горизонтальных проложений от исходного до данного пункта из предыдущего и последующего циклов измерений. При этом горизонтальное смещение исходного пункта должно быть учтено в расчете смещения рабочего пункта.
9.5.Значения горизонтальных деформаций ε определяются как отношение
разницы горизонтальных сдвижений ξ к горизонтальной проекции длины интервала между пунктами наблюдений.
-28 -
9.6.При построении графиков (схем) сдвижений и деформаций следует руководствоваться следующими правилами:
– графики вертикальных сдвижений и деформаций строят по каждой профильной линии нивелирования, а в случае необходимости – дополнительно по их отдельным участкам;
– по оси абсцисс откладывают горизонтальные расстояния по профильной линии от опорного (узлового) репера, а по оси ординат – рассчитанные величины вертикальных сдвижений и деформаций;
– вертикальный и горизонтальный масштабы отображения подбирают, исходя из информативности и наглядности картины исследуемого процесса;
– на графиках показывают распределения предельных средних квадратических погрешностей рассчитанных величин сдвижений и деформаций,
атакже границы пересечения профильных линий с контуром нефтеносности;
– при выявлении зон опасного сдвижения земной поверхности (см. приложение 1) графики совмещают с геологическими разрезами, промысловыми (скважины, ДНС, трубопроводы и т.п.) и гражданскими объектами;
– построение схем горизонтальных сдвижений и деформаций по данным спутниковых наблюдений регламентировано п. 6.3.18 данной инструкции.
10.ПЕРЕЧЕНЬ МАТЕРИАЛОВ, ПОДЛЕЖАЩИХ СДАЧЕ
10.1.При завершении очередного цикла наблюдений Заказчику подлежат сдаче следующие материалы:
–план-схема наблюдательной станции с изображением:
•опорных и рабочих реперов нивелирной сети, профильных линий нивелирования;
•опорных и рабочих пунктов спутниковой сети, базисных линий;
•контура ВНК объекта мониторинга и внешней границы зоны опасного влияния его разработки;
•промышленной инфраструктуры объекта мониторинга (пром. объекты, пром. площадки, инженерные коммуникации);
•элементов топографии (населенные пункты, гражданские объекты, дорожная сеть, гидрография, растительность и пр.);
•границы тектонических блоков, разломы, крупные линеаменты, выявленные по результатам аэрокосмогеологических исследований, структурного бурения и др. методов;
•иных источников деформаций земной поверхности (мест возможного возникновения оползневых и карстовых процессов, горных отводов действующих и проектируемых рудников, соседних нефтяных месторождений и т.п.), совмещенных с объектом мониторинга или находящихся в непосредственной близости от него;
•другой графической информации, необходимой для выявления основных закономерностей процесса сдвижения.
(Масштаб план-схемы выбирают, исходя из достаточной информативности и наглядности изображения объекта мониторинга и наблюдательной станции. Рекомендуемым масштабом является 1:25 000. Заказчик имеет право в техническом задании на выполнение работ установить требуемые значения масштабов выходной графической документации);
-29 -
–каталоги уравненных координат и высот пунктов наблюдательной станции
исредние квадратические ошибки их определения относительно опорных пунктов;
–ведомости вертикальных и горизонтальных смещений рабочих пунктов наблюдательной станции, величины горизонтальных деформаций на локальных наблюдательных станциях и средние квадратические ошибки их определения;
–графики оседаний реперов по профильным линиям нивелирования и предельных погрешностей их определения;
–схемы горизонтальных смещений пунктов спутниковой сети и предельных погрешностей их определения.
10.2.Дополнительные материалы, касающиеся используемых инструментов
иоценки точности выполненных работ:
–результаты первичной (периодической) метрологической поверки используемых приборов (нивелира, реек, спутниковых приемников, антенн);
–ведомости невязок в замкнутых полигонах нивелирной (спутниковой) сети и их допустимые значения;
–определение случайных средних квадратических ошибок нивелирования по разностям двойных измерений и результатам строгого уравнивания;
–исследование остаточных систематических ошибок нивелирования;
–анализ устойчивости опорных реперов (пунктов);
–расчет точности определения оседаний;
–расчет точности определения горизонтальных сдвижений пунктов спутниковой сети.
10.3.По требованию заказчику передаются полевые журналы, данные измерений электронных приборов (нивелира, спутникового приемника) в одном из обменных форматов (RINEX, GSI, M5).
10.4.По завершении очередного цикла наблюдений организацияисполнитель выполняет анализ и обобщение результатов работ с написанием технического отчета (приложение 7). Полный анализ и обобщение результатов работ выполняют не реже одного раза в пять лет специализированные организации, привлекаемые нефтедобывающими предприятиями по согласованию с территориальными органами МПР. Результаты анализа
оформляются в виде отчетов и хранятся в фондах данных предприятий. О выполненных работах информируются контролирующие и заинтересованные организации. По запросам оформленные отчеты передаются калийным предприятиям, Пермгеолкому, Западно-Уральскому управлению Ростехнадзора.
- 30 -
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ
1.Инструкция по созданию наблюдательных станций и производству инструментальных наблюдений за процессами сдвижения земной поверхности при разработке нефтяных месторождений в регионе Верхнекамского месторождения калийно-магниевых солей (ВКМКС). – Пермь: Изд-во Перм. гос.
техн. ун-та, 2003. – 34 с.
2.Инструкция по производству маркшейдерских работ (утв. Постановлением Госгортехнадзора России № 73 06.06.2003 г.). – М., 2003.
3.Правила промышленной безопасности при освоении месторождений нефти на площадях залегания калийных солей (Постановл. Госгортехнадзора России от 4.02.2002 г. за № 8). – М., 2002.
4.Исследование влияния добычи нефти на напряженно-деформированное состояние горного массива в районе ВКМКС: отчет о НИР / руков. работы Ю.А. Кашников // Фонды ЗАО «ЛУКОЙЛ-ПЕРМЬ». – Пермь, 1998. – 190 с.
5.Производство инструментальных наблюдений на наблюдательной станции Сибирского месторождения: отчет о НИР / руков. работы Ю.А. Кашников // Фонды ЗАО «ЛУКОЙЛ-ПЕРМЬ». – Пермь, 1999 г. – 87 с.
6.Правила закладки центров и реперов на пунктах геодезической и нивелирной сетей. – М.: Картгеоцентр: Геодезиздат, 1993.
7.Инструкция по нивелированию I, II, III, IV классов. ГКИНП (ГНТА) -03-010- 03 / Федерал. служба геодезии и картографии России. – М.: Изд-во ЦНИИГАиК, 2003.
8.Измерение вертикальных смещений сооружений и анализ устойчивости реперов / В.Н. Ганьшин, А.Ф. Стороженко, А.Г. Ильин [и др.]. – М.: Недра, 1981. – 215 с.
9.Лазарев Г.Е., Самошкин Е.М. Основы высшей геодезии. – М.: Недра, 1980. – 424 c.
10.Кашников Ю.А., Ашихмин С.Г. Механика горных пород при разработке месторождений углеводородного сырья. – М.: Недра, 2007. – 467 с.
11.Механика горных пород применительно к проблемам разведки и добычи нефти / под. ред. В. Морн, Д. Фурментро. – М.: МИР: Эльф-Акитен, 1994.
12.Кузьмин Ю.О. Современная геодинамика и оценка геодинамического риска при недропользовании. – М.: Агентство экономических новостей, 1999. – 220 с.
13.Иофис М.А., Шмелев А.И. Инженерная геомеханика при подземных разработках. – М.: Недра, 1985. – 248 с.