Методика физико-геологического моделирования объектов с переменной плотностью и намагниченностью

Виноградов, Владислав Борисович

Бесплатный автореферат и диссертация по наукам о земле на тему
Методика физико-геологического моделирования объектов с переменной плотностью и намагниченностью
ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

Автореферат диссертации по теме "Методика физико-геологического моделирования объектов с переменной плотностью и намагниченностью"

На правах рукописи

Г)

ВИНОГРАДОВ Владислав Борисович

МЕТОДИКА ФИЗИКО-ГЕОЛОГИЧЕСКОГО МОДЕЛИРОВАНИЯ ОБЪЕКТОВ С ПЕРЕМЕННОЙ ПЛОТНОСТЬЮ И НАМАГНИЧЕННОСТЬЮ

Специальность: 25.00.10-«Геофизика, геофизические методы поисков полезных ископаемых»

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата геолого-минералогических наук

2 6 НОЯ 2009

Екатеринбург - 2009

003484966

Работа выполнена в ГОУ ВПО «Уральский государственный горный университет»

Научный руководитель - доктор геолого-минералогических наук,

профессор, заслуженный геолог Российской Федерации Филатов Владимир Викторович

Официальные оппоненты: доктор геолого-минералогических наук,

доцент Писецкий Владимир Борисович

кандидат геолого-минералогических наук Глухих Игорь Иванович

Ведущая организация - Горный институт УрО РАН

Защита диссертации состоится «3» декабря 2009 г. в 1430час. на заседании диссертационного совета Д 212.280.01 при ГОУ ВПО «Уральский государственный горный университет» по адресу:

620144, г. Екатеринбург, ул. Куйбышева, 30, 3-й учебный корпус, ауд. 3326.

С диссертацией можно ознакомиться в научной библиотеке ГОУ ВПО «Уральский государственный горный университет»

Автореферат разослан «2» ноября 2009 г.

Ученый секретарь диссертационного совета --

А.Б. Макаров

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Для достижения цели были решены следующие задачи:

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

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

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

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

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

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

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

Защищаемые положения:

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

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

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

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

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

Апробация работы. Результаты проведенных исследований обсуждались на сессиях международного семинара им. Д.Г. Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (Днепропетровск, 1991; Ухта, 1998, 2008; Екатеринбург, 1999, 2002, 2006; Пермь, 2005; Казань 2009), на конференции «Геофизические методы при разведке недр и экологических исследованиях» (Томск, 2002), на региональных научно-практических на конференциях «Геология и полезные ископаемые Западного Урала» (Пермь, 2008, 2009), на «Пятых научных чтениях памяти Ю.П. Булашевича» (Екатеринбург, 2009).

Практическая значимость работы. Разработанные петрофизические модели могут применяться при интерпретации геофизических данных для решения задач как нефтегазовой, так и рудной геологии разных масштабов. Результаты исследований внедрены в ПГО «Севвостгеология» и ПО «Сильвинит». Методики оценки НДС могут применяться при истолкова-

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

Достоверность получаемых результатов основана на результатах расчетов, выполненных для теоретических моделей.

Фактический материал и личный вклад. Представленные в работе методики моделирования сложных геологических сред, программы расчета полей и параметров НДС созданы автором. Практические примеры их применения основаны на материалах предоставленных для выполнения хоздоговоров, а некоторые из них основаны на опубликованных в печати материалах, что оговорено в тексте.

Публикации. По теме диссертации опубликовано 14 работ, в том числе две статьи в журналах из перечня ВАК.

Объем работы и структура. Работа состоит из введения, четырех глав, заключения, общим объемом 105 страниц, 55 иллюстраций, 3 таблиц, списка литературы, включающего 110 наименований.

Работа выполнена на факультете геологии и геофизики УГГУ под руководством доктора геолого-минералогических наук В.В. Филатова, которому автор выражает благодарность за всестороннюю помощь и поддержку при выполнении работы.

Автор благодарит своих коллег В.И. Бондарева, С.М. Крылаткова, A.C. Долгаля, О.В. Зотеева, A.B. Давыдова, Ю.Н. Субботина, В.М. Сапожнико-ва, В.В. Бабенко, В.А. Кочнева, чьими советами и помощью он неоднократно пользовался.

Автор выражает благодарность участникам семинара им. Д.Г. Успенского за обсуждение работы, сотрудникам объединения «Уралкалий», ПГО «Сев-востгеология», ФГУП «Челябгеолсъемка», ОАО «Баженовская экспедиция», оказавшим помощь в процессе работы и предоставившим материалы для работы.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

В четвертой главе описана методика моделирования геологической среды с целью вычисления параметров НДС по гравитационному и магнитному полям с учетом формы геологических объектов.

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

Защищаемые положения

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

В геологической практике встречается как скачкообразное, так и непрерывное изменение физических свойств (рис. 1). В прикладной геофизике при проведении интерпретации используются преимущественно модели с кусочно-постоянными свойствами. Подробно для эндогенных месторождений это явление рассмотрено Г.А. Соловьевым. Физико-геологические модели для отдельных типов месторождений с непрерывным латеральным изменением свойств предложены А.Н. Тимофеевым, Ю.П. Булашевичем, В.М. Новоселицким, Л.Я. Ерофеевым, З.М. Слепаком, В,И. Гольдшмидтом, В.Н. Страховым, А.И. Кобруновым и другими. Влияние формы объекта на распределение магнитных свойств рассмотрено Ю.И. Блохом. Закономерности изменения физических свойств под влиянием механических напряжений подробно рассмотрены В.В. Филатовым.

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

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

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

Рис. 1. Вертикальные карты изолиний магнитной восприимчивости по четырем профилям на одном из эксплуатационных блоков Гусевогорского месторождения (расстояние между профилями 5 м)

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

ления полей является главным критерием практической применимости методики моделирования на основе принятых математических конструкций.

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

Автором предложены и детально разработаны способы описания геологических поверхностей в зависимости от их вида с помощью сплайн-функций различного типа: кубических, тригонометрических, параметрических [2*,1,2]. Один из них - описание поверхности кубическим сплайном:

г(х)=а1+ь,{х-х!)+с1(х-х1у , хе[х„ х,н], (1)

где а\, Ь\, с„ - коэффициенты сплайна; х, - узлы сплайна; г(х) - аппликата поверхности. Как частный случай выражение (1) включает в себя описание негладкой поверхности, в этом случае с,- и А в выражении (1) равны нулю. Совместное применение кусочно-постоянных функций и сплайн-функций для описания поверхностей обеспечивает адекватное описание любой геологической поверхности.

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

0(лс)=а(х,)+Ь1(х-х,)+с,(х-х,У+с!,(х-х,У , хе[х„ (2)

Описание сплайнами закона изменения намагниченности в пространстве позволяет строить эффективные алгоритмы расчета магнитных полей вблизи аномалиеобразующих тел [2].

Физико-геологические сплайн-модели позволяют заменить петрофизи-ческие модели А.Н. Тимофеева, Ю.П. Булашевича, Г,А. Соловьева и других авторов, что подтверждается проведенными расчетами и теоретическими преобразованиями (рис. 2). Для этого необходимо правильно выбрать расстояние между узлами сплайна, которое может быть вычислено теоретически. Например, для модели Г.А. Соловьева расстояние между узлами

сплайна (Ах) определяется неравенством:

х , (3)

к4 К\ 2

где Х\-&Х - магнитная восприимчивость петрофизической 1-й зоны и ее избыток; к, - ширина 1-й переходной области между петрофизическими зонами.

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

определено кубическим сплайном. В точках оси X (х,0,0) аномалия определяется формулой [6]:

4Г(*)=£

,^23. 4

л\

+и

4+^(382+«2)

1п/~

3 4

/ Л2 2

-их -е

1 — $ , 1и е

—+е\ГГб

Ш "I „

яагар — +—Ши + & б

02 | Ь2 | с2 о/ М с/

(4)

где х,- - абсциссы узлов сплайна; х - абсцисса точки расчета; ,

и=у-у, хо, у о', ъ> - координаты источника; R=4t*+lF+si ;

Лй=о1+Ъ1(х-х1) + с1{х-х1)2+с11{х-х^, | Й | " I * " характерные размеры прямоугольного параллелепипеда по осям х, у, г] ст» Ъ„ с,, -коэффициенты сплайна плотности.

А/м

Модель Булашевича

Модель Тимофеева

х, м

Рис. 2. Аппроксимация моделей латеральной изменчивости намагниченности А.Н. Тимофеева и Ю. П. Булашевича сплайн-моделями - точками на графике показаны узлы интерполяции

б. г/см3

Рис. 3. Описание кубическими сплайнами латеральной изменчивости физических свойств: а - магнитных свойств пород в зоне влияния гранитного интрузивного массива (по В.В. Бродовому и др.); б - плотности в отложениях турнейского яруса Ямашинской структуры; в - плотности в отложениях ассельского яруса Улеминского поднятия (по З.М. Слепаку)

Магнитное поле модели, состоящей из N параллелепипедов, вектор намагниченности которых направлен вертикально, а его величина относительно мала и определена кубическим сплайном:

2П (х—хЛу,.

¿ «ОД5-^+Ых~х0)+Ю+

У о Я

+ /),,„,„( Д+^^А)) | 5 | £ | (5)

Уо "о

где | | " | % - характерные размеры прямоугольного параллелепипеда по осям х, у, 2 .

Сеточная модель содержит параллелепипеды разных размеров, незакономерно расположенных в пространстве. Формула для вычисления интенсивности магнитного поля такой сеточной модели, состоящей из К параллелепипедов с латеральным изменением намагниченности, примет вид [10, 12]:

и-*о)2о

4=1 ¡=1

Д С,{У-У0)

У-Уо

' (у~Уо)Л , -+

У-У о Я

~--С,[у-у0)

апИ-

л лл

агЛ- I „2 \ Ь2 I с2

+-[х-х0)-А] + 2С.1П(2-Л) + 20,Л1 I I " I , (6)

У-Уд 1 1 1

где Я=М'{х-х0)2Цу-у0)2+(1-г0)г ; С1=чгЪЛ,г1 ; ;

; ¡>¡=<1^ ^/с - число звеньев сплайна в Л-м параллелепипеде,

| п/ | Й | с; - характерные размеры прямоугольного параллелепипеда по осям х, у. 2.

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

х[х)=ах1+Ьх2+сх+<1у3+еу2+1у+£23+кг2 + \>г+™ . (7)

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

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

Страхова В.Н. и Бродского М.А. установлена единственность решения обратной задачи в вышеприведенной постановке, если проведена точная непрерывная регистрация поля на изучаемом участке. Однако наблюдаемые поля всегда осложнены ошибками, наблюдения проводятся в узлах дискретной сети пунктов на ограниченном участке, поэтому решение всегда будет представлено семейством е-эквивалентных по полю сплайнов намагниченности (плотности при интерпретации поля Дg). В тексте диссертации приведен пример построения семейства эквивалентных по полю решений обратной задачи магниторазведки.

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

Рис. 4. Обобщенная физико-геологическая модель зоны разлома: графики изменения плотности (а), поля силы тяжести (б), магнитной восприимчивости (в), магнитного поля (г): 1 - зона упругих деформаций, 2 - зона пликативных деформаций, 3 - зона дробления

На рис. 5 представлена физико-геологическая модель Байкальского разлома, построенная по опубликованным данным о непрерывным изменении степени катаклаза и плотности.

Расчеты полей на основе описания физических свойств кубическими сплайнами проводились на Гусевогорском железорудном месторождении, молибдено-вольфрамовом месторождении, Карамкенском золоторудном месторождении [1*, 6, 10].

Второе защищаемое положение: установлено три типа зависимости параметров напряженно-деформированного состояния среды н силы

тяжести для тел простой геометрической формы: взаимно однозначная, корреляционная и многозначная.

•). ДдмГал

-2,65

-2,60

■---- 3 ~ ---------- ■ 6 —1 9 12 х, км

-н--Н- -Р -ь + 44- 4- 4--т. 44- ~ + 4- ~ ~ 4" 4- 4-4- 4- -¥- 4=-

и г+ча паз игл*

Рис. 5. Петроплотностная модель деформированной зоны, построенная на основе аппроксимации плотности сплайном: 1 - породы со слабым катаклазом; 2 - бластомилониты; 3 - породы с умеренным катаклазом; 4 - породы с сильным катаклазом; 5 - сплайн плотности; 6 - теоретическая кривая Д§ кусочно-постоянной модели; 7 - теоретическая кривая Дg сплайн-модели

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

К настоящему времени имеется решение прямой задачи тектонограви-метрии во второй постановке для материальной точки. Впервые решение получено Р. Миндлиным и Д. Ченем методом зеркальных изображений для компонентов вектора смещения и тензора деформации, обусловленных действием сосредоточенной силы в упругом полупространстве. Позднее решения этой задачи другими методами получены А. И. Лурье и К. Бреб-бия.

В. В. Филатовым также на основе решения Миндлина получены интегральные соотношения для компонентов вектора смещения плотностной неоднородности произвольной формы, находящейся в упругом полупространстве. Компоненты вектора смещения получены им в виде произведения вторых производных пзавитационного потенциала на некотопые гЬгнк--

ции координат. В.В. Филатовым и Л.А. Болотновой проведена оценка НДС геологической среды г. Екатеринбурга и его окрестностей и других территорий.

Рассмотрим НДС, обусловленное телами правильной геометрической формы, для которых известны формулы силы тяжести. Для этого сначала получим выражения для параметров напряженного состояния среды для объектов произвольной формы. С помощью соотношений Коши получены компоненты тензора деформации объектов произвольной формы на поверхности полупространства (г=0). Одна из них имеет вид:

_ 1+у) 8лЕ(1 -V)

Щ-Цх-х{

--->2 п2

5 ¿+(3-4УД

1 К2

. -Ч х—Хп) } 2(г+г0 2 аУ—

Зр^г(1 + у)с

1—2У)( 1 + У) с К\(И+г+г0)-{х-ха)1(2Ъг + 1+20) '

I Л^+г+Г,) ' ^

где g без потери точности вычисляемых параметров можно считать постоянным и равным среднему нормальному значению по Земле; р - относительная плотность; V - коэффициент Пуассона; Е - модуль Юнга.

В гравиразведке самым распространенным элементарным объектом, применяемым при решении прямых и обратных задач, является прямоугольный параллелепипед. Набором прямоугольных параллелепипедов удобно аппроксимировать геологические объекты, ограниченные поверхностями с субвертикальным падением, и крутопадающие разломы. Методы построения плотностной модели на основе аппроксимации среды прямоугольными параллелепипедами хорошо освоены, подбор ведется в диалоговом или автоматизированном режиме (В.А. Кочнев, Е.Г. Булах и др.) [7].

На поверхности упругого полупространства (г=0) одна из формул для

вычисления КТД прямоугольного параллелепипеда примет вид

_¿

- | « | ^ (9)

где | °с! | " | Й - характерные размеры прямоугольного параллелепипеда по осям х, у, 2.

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

Автором получены зависимости горизонтальной и вертикальной компонент вектора смещения и компонент тензора чистой деформации единичной массы (М=1) от величины [5]:

Ех„

~ё0

Яо-

1 —2У

2,39(1-V2)

Ех„

ё1

2(1-

1 + у

—~ё о

Ех,

0.5

ёК2-3^) | 2-2

,(1+УЫ(ЗЕо-2У) £(1-у)4,5

(10)

где ё0 - относительная величина силы тяжести. Зависимости для других КТД приведены в тексте диссертации.

Из формул следует, что связь КТД, КВС и для материальной точки функциональная, однозначная. То же можно утверждать относительно других параметров НДС. По известным значениям гравитационного поля можно вычислить компоненты тензора деформации. Зависимости и(А£), е!/(4?) имеют вид дробно-рациональной функции. Амплитуда вертикальных смещений м> монотонно возрастает с увеличением поля силы тяжести, зависимость м> и Дg взаимно однозначная. Зависимость горизонтальных смещений и от величины поля силы тяжести выражается однозначной функцией, но не является взаимно однозначной. Графики этих зависимостей приведены на рис. 6.

Для изучения связи КТД и Ag горизонтального материального стержня воспользуемся известными выражениями Д^, а КТД найдены автором. КТД материального стержня не зависят от линейной плотности стержня:

Графики зависимостей ен(Л§), построенные по результа-

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

Для вертикального материального стержня уравнения в замкнутом виде исследуемой зависимости получить не удается, несмотря на то, что аналитические выражения имеются и для КТД, и для Д§. Тем не менее зависимость и в этом случае функциональная, однозначная, что подтверждено численными расчетами. Графики исследуемых зависимостей, построенные по результатам вычислений н> и Д^ для вертикального и горизонтального стержня, представлены на рис. 7. Для горизонтального и вертикального стержней по величине силы тяжести можно вычислить характеристики НДС.

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

чены на основе численных расчетов. Они являются однозначными функци-

ями своего аргумента. Для вертикального уступа, флексуры, наклонной плоскости по известным величинам Ag можно вычислить параметры НДС. и.\Л/ ^

\Л/

а б

Рис. 6. Зависимость КВС (а) и КТД (б) от величины поля силы тяжести на поверхности упругого полупространства, напряженное состояние которого обусловлено материальной точкой

В результате проведенных расчетов установлено, что для наклонных пластов, горизонтальной пластины и некоторых объектов правильной геометрической формы зависимость характеристик НДС от ^ корреляционная [5]. Примеры корреляционных зависимостей силы тяжести и КТД для Днепровского массива приведены на рис. 8, а,б.

Для прямоугольного параллелепипеда и горизонтальной пластины выражения КТД и очень сложные. Прямоугольный параллелепипед с квадратным сечением в плоскости ХО У имеет узкую полосу значений е,,, соответствующих одному значению Ag. Степень многозначности зависимости параметров НДС от величины поля силы тяжести меняется при изменении сечения параллелепипеда.

а б

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

Рис. 8. Аппроксимация зависимостей КТД е,у от величины силы тяжести кубическим сплайном (Днепровский массив). Кружками показаны точки, в которых известны Дя и е,>

На основе зависимостей между величиной КВС и КТД и интенсивностью силы тяжести для тел правильной геометрической формы автором предложено проводить оценку НДС по следующей схеме [5]:

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

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

по которьм отыскиваются главные значения и главные направления тензора деформации. Для этой группы также возможен прямой пересчет д^ в КТД. Связь между параметрами НДС и величиной силы тяжести корреляционная.

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

Форма геологических объектов, их распределение в пространстве, плотность уникальны пля конкоетной геологической ситуации. Рассмотренные

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

Если результаты интерпретации поля силы тяжести позволяют выбрать объект, для которого зависимость параметров НДС от величины поля силы тяжести является однозначной функцией, то возможен пересчет Д^ в характеристики НДС. Необходимо выбрать 50-70 точек наблюдения поля силы тяжести, равномерно покрывающих диапазон изменения для выбранных точек по плотностной модели вычислить КТД и определить коэффициенты аппроксимирующей конструкции. Далее, используя полученную зависимость Еу(Дя), для всех точек вычисляют КТД, а по ним, в свою очередь, главные значения и главные направления тензора деформации, инварианты тензора деформации и иные характеристики НДС.

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

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

По результатам измерений силы тяжести была проведена оценка напряженного состояния для южной части Верхнекамского месторождения калийных солей. Результаты вычислений приведены на рис. 9 [3]. Главные горизонтальные деформации вблизи земной поверхности представляют собой растяжения, за исключением юго-западной области, где среда испытывает сжатие. Третья главная деформация направлена вертикально и всюду, кроме южной части, представляет собой сжатие.

Наиболее общим подходом при оценке НДС по полю силы тяжести является решение прямой геомеханической задачи для плотностной модели исследуемой территории. Методика создания и верификации плотностной модели описана в работе [8]. Для оценки НДС существенную роль играет адекватность интерпретационной модели реальной ситуации. Оценка НДС

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

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

Рис. 9. Результаты расчетов компонентов тензора деформаций вблизи земной поверхности в южной части Верхнекамского месторождения: горизонтальные компоненты (а) и план изолиний вертикальной компоненты тензора деформации южной части Верхнекамского месторождения (б) (по В.В. Филатову и В.Б. Виноградову):

1 - величина и направление деформации сжатия; 2 — величина и направление деформации растяжения; 3 - населенные пункты; 4, 5, 6 - изолинии егг положительные, нулевые, отрицательные соответственно

ПОЕЗСШз

При проведении работ нередко оказывается, что сеть магнитных наблюдений гораздо детальнее, чем гравиметрических. В таких случаях необходимо на основе аналитической аппроксимации магнитного поля вычислить псевдогравитационное поле и по нему выполнить оценку НДС [9]. После построения карт изолиний параметров НДС выделяют зоны максимальных и минимальных деформаций, области сжатия и растяжения. Проводят их анализ с привлечением всех имеющихся геолого-геофизических данных. Горные породы сопротивляются растяжению намного хуже, чем сжатию, поэтому зонам растяжения уделяют особое внимание. Пример вычисления Д^рэ по данным магнитной съемки масштаба 1:10000 для небольшого участка Карамкенского месторождения размерами 2x2 км приведен на рис. 10 [11].

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

х/

-и» \

/ ^ / ^

Ж <

■^МЙ

-яр

.У

(а

/ /

I г-I

7 7 —

У ' ~

У / /

" /( / I

\ I I

1 > ;

К../- ; С, ^

а б в

Рис. 10. Магнитное поле (в нТл) (а), псевдогравитационное поле (в мГал) (б) и план изолиний вертикальной компоненты тензора деформации (в) одного участка Карамкенского месторождения

К участкам оруденения приурочены мощные зоны разуплотнения, их размеры по простиранию достигают 2 км, ширина до 300 м, глубина превышает 500 м, плотность пород в зоне разуплотнения уменьшается до 2,40 г/см3, а дефект плотности по отношению к вмещающим породам достигает -0,2 г/см3. Магнитное поле рудных зон спокойное, с отрицательными значениями до 100 нТл.

Была создана объемная физико-геологическая модель зоны разуплотнения, составленная из прямоугольных параллелепипедов с разной глубиной до нижней кромки (рис. 11) [11]. Для модели были вычислены значения дилатации и коэффициентов Лодэ - Надаи. Результаты расчетов в виде планов изолиний приведены на рис. 11. Даже для простой модели план изолиний коэффициента Лодэ - Надаи имеет очень сложный вид, на нем присутствует 4 экстремума. По локальным аномалиям поля силы тяжести были вычислены главные компоненты ТД, построен план изолиний дилатации (рис. 12). По совокупности геолого-геофизических данных автором построена схема основных сжимающих и растягивающих усилий.

Мозаично-блоковое геологическое строение участка обусловило сложный характер наблюдаемых физических полей и поля дилатации. Палео-кальдера характеризуется знакопеременными значениями дилатации от -20-10 до 15-10, т. е. амплитуда изменений невелика. Хорошо в поле дилатации отразились разрывные нарушения широтного и северо-западного простирания. Четко на плане проявлен наиболее древний и долгоживущий разлом субмеридионального простирания.

Рис. 11. План изолиний коэффициента Лодэ - Надаи (1), план изолиний дилатации (2) и контур плотностной модели зоны разуплотнения (3)

На плане изолиний дилатации Карамкенского участка в околоразлом-ной зоне меридионального разлома наблюдается аномалия, аналогичная

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

Рис. 12. План изолиний дилатации Карамкенского участка -оцифровка изолиний дилатации в 10"9: 1 - изолинии отрицательных значений; 2 - изолинии нулевых значений; 3 - изолинии положительных значений; 4 - разрывные нарушения; 5 - контур палеокальдеры; оцифровка изолиний в 10"9

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

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

В работе представлена технология описания геологических поверхностей сплайнами различного вида, в зависимости от характера поверхности.

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

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

Заключение

Получены выражения расчета гравитационного и магнитного полей для моделей с произвольным законом изменения плотности и намагниченности в вертикальном и горизонтальном направлениях. Предлагаемая технология позволяет вычислить погрешность, с которой проводится аппроксимация среды. Для расчетов составлены программы для ЭВМ.

Разработаны алгоритмы решения обратных задач гравиразведки и магниторазведки для геологической среды с непрерывным пространственным изменением физических свойств.

Установлена взаимосвязь между параметрами напряженно-деформированного состояния геологической среды и величиной силы тяжести для объектов простой геометрической формы.

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

Создана методика расчета параметров НДС, учитывающая геометрическую форму изучаемых геологических тел. Методика реализована в программном обеспечении ЭВМ.

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

Основные положения диссертация опубликованы в работах: Статьи, опубликованные в ведущих рецензируемых журналах, входящих в перечень ВАК

1". Виноградов В.Б. О методике моделирования железорудных тел //Изв. вузов. Горный журнал. - 200В. - № 8. - С. 85 -88.

2'. Виноградов В.Б. Способ описания поверхности геологических объектов сплайнами// Литосфера. - 2009,- № 3, - С. 91 - 93.

Работы, опубликованные в других изданиях

1. Виноградов В.Б., Филатов В.В. Применение сплайн-функций для описания поверхности геологических объектов // Геофизические методы поисков и разведки рудных и нерудных месторождений. Межвуз. науч. темат. сб. Вып. 12 - Свердловск, 1986. - С. 51 -55.

2. Филатов В.В., Виноградов, В.Б. Прямая задача магниторазведки на основе представления поверхности геологических объектов сплайнами // Геофизический журнал. -1988. - Т. 10, №2. - С. 67-73.

3. Виноградов В.Б. Определение напряжений и деформаций в земной коре по полю силы тяжести // Геофизические аспекты изучения геологического строения месторождений калийных солей: сб. науч. трудов ВНИИГ. - Л., 1989. - С. 162 - 167.

4. Виноградов В.Б., Кашкаров А.А. Оценка напряженно-деформированного состояния геологической среды по геофизическим данным // Инф. лист. № 287-90. - Свердловск, 1991.-4 с.

5. Виноградов В.Б. Методика расчета деформаций среды по гравитационному полю // Известия УГГТА. - Вып. 5. - Сер.: Геол. и геофизика. - Екатеринбург, 1996. - С. 110111.

6. Виноградов В.Б. Магнитные и плотностные интерпретационные сплайн-модели/ В.Б. Виноградов // Известия УГГГА. Вып. 8. Сер.: Геология и геофизика. — Екатеринбург, 1998. —С. 142-144.

7. Иванов Д.Б., Виноградов В.Б., Чурсин A.B. Моделирование угленосных структур по данным магниторазведки // Геофизические методы при разведке недр и экологических исследованиях. - Вып. 2. - Томск; Изд. ТГУ, 1999. - С. 212 - 214.

8. Виноградов, В.Б. О методах интерпретации гравитационных и магнитных полей восточноуральских угольных месторождений // Известия УГТТА. - Вып. 15. Сер.: Геология и геофизика. - Екатеринбург, 2002. —С. 209 - 214.

9. Виноградов В.Б. Выбор вида и параметров аналитической аппроксимации для расчета трансформаций // Вопросы теории и практики геологической интерпретации гравитационных магнитных и электрических полей: материалы 32-й сессии семинара им. Д.Г. Успенского. - Пермь: Горный институт УрО РАН, 2005. - С. 44.

10. Виноградов В.Б. Совершенствование методики обработки данных КМВ на железорудных месторождениях // Геология и полезные ископаемые Западного Урала: материалы региональной научно-практической конференции. - Пермь: Пермский гос. унт, 2008. - С. 224-227.

11. Виноградов В.Б. Модели золоторудного месторождения// Вопросы теории и практики геологической интерпретации гравитационных магнитных и электрических полей: материалы Зб-й сессии семинара им. Д.Г. Успенского. - Казань: Изд-во Казан, гос. ун.-та, 2009. - С. 86 -89.

12. Виноградов В.Б. Моделирование объектов с переменной намагниченностью// Пятые научные чтения памяти Ю.П. Булашевича: материалы.-Екатсркнбург: ИГф УрО РАН, 2009. - С. 74 - 77.

Подписано в печать 23.10.09 г. Формат 60x84 1/16. Бумага офсетная. Печать на ризографе. Печ. л. 1,0. Тираж 100. Заказ 1%>'2-

Отпечатано с оригинал-макета в лаборатории множительной техники издательства УГТУ

620144, г. Екатеринбург, ул. Куйбышева, 30 Уральский государственный горный ун-т

Содержание диссертации, кандидата геолого-минералогических наук, Виноградов, Владислав Борисович

Введение.

Глава 1. Физико-геологическое моделирование в гравиразведке и магниторазведке.

1.1 . Моделирование геологических объектов.

1.2. Напряженное состояние земных недр.

1.3. Задачи тектоногравиметрии и тектономагнитометрии.

Глава 2. Интерпретационные модели с произвольным законом изменения плотности и намагниченности.

2.1. Описание поверхности геологических тел кубическими сплайнами.26 сплайн-функци

2.2. Аппроксимация изменчивости плотности и магнитных свойств.

2.3. Магнитные модели с латеральной изменчивостью намагниченности.

2.4. Модели с изменяющейся по латерали плотностью.

2.5. Решение обратных задач для сред с непрерывным латеральным изменением физических свойств

Глава 3. Взаимосвязь гравитационного поля и напряженного состояния верхней части земной коры.

3.1 Напряженное состояние упругого полупространства, обусловленное сосредоточенной силой.

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

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

3.4 Зависимость параметров напряженного состояния от величины силы тяжести для сложных моделей.

Глава 4. Методика тектонофизического анализа состояния геологической среды.

4.1. Методика оценки напряженного состояния геологической среды по полю силы тяжести.

4.2. Методика оценки напряженного состояния геологической среды по магнитному полю.

4.3. Интерпретационные модели разломных зон.

4.4. Оценка напряженного состояния геологической среды.

4.4.1. Карамкенский участок.

4.4.2. Полетаевская площадь.

4.4.3. Балахонцевский участок.

Введение Диссертация по наукам о земле, на тему "Методика физико-геологического моделирования объектов с переменной плотностью и намагниченностью"

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

Изучение напряженного состояния геологической среды позволяет прогнозировать динамические явления на ранних этапах геологических исследований, а также влечет за собой расширение круга решаемых задач в условиях сокращения рынка геофизических услуг. Всестороннее изучение геологической среды позволяет полнее использовать имеющиеся ресурсы. Актуальность темы подтверждается публикациями на эту тему в ведущих геофизических журналах страны (Булах Е. Г. 2007, 2008; Мартышко П. С. 2005 и другие).

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

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

Для достижения цели были решены следующие задачи:

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

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

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

Защищаемые положения:

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

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

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

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

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

Апробация работы. Результаты проведенных исследований обсуждались на 5 Уральской конференции "Применение математических методов и ЭВМ при обработке информации на геологоразведочных работах" (1986), на научно-технической конференции Свердловского горного института'(1990), на сессиях международного семинара им. Д.Г. Успенского «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (Ухта, 1998, 2008; Екатеринбург, 1999, 2002, 2006; Пермь, 2005, Казань, 2009), на конференции «Геофизические методы при разведке недр и экологических исследованиях» (Томск, 2002), на региональных научно-практических конференциях «Геология и полезные ископаемые Западного Урала» (Пермь, 2008, 2009), на «Пятых научных чтениях памяти Ю.П., Булашевича (Екатеринбург, 2009).

Практическая значимость работы; Разработанные петрофизические модели могут применяться при интерпретации геофизических данных для решения задач как нефтегазовой, так и рудной геологии разных масштабов. Результаты исследований внедрены в ПГО «Севвостгеологйя» и ПО «Сильвинит».

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

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

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

Публикации. По теме диссертации опубликовано 14 работ, в том числе две статьи в ведущих рецензируемых журналах, входящих в перечень ВАК.

Объем работы и структура. Работа состоит из введения, четырех глав, заключения, общим объемом 105 страниц, 55 иллюстраций, 3 таблицы, списка литературы, включающего 110 наименований.

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

Автор благодарит своих коллег проф. В. И. Бондарева, доц. С. М. Крылатко-ва, д.г.-м. н. А. С. Долгаля, проф. А. В. Давыдова, чл. корр. РАН Ю. Н. Субботина, проф. В. М. Сапожникова, проф. В. В. Бабенко, д.ф.-м. н. В. А. Кочнева, проф. В. Б. Писецкого, д. т. н. О. В. Зотеева, доц. Ж. Н. Александрову, нач. партии Л. Д. Нояксову, доц. В. Е. Петряева советами и помощью которых он неоднократно пользовался. В ходе работы автор получал помощь многих сотрудников и студентов ИГИГ. Всем указанным лицам автор выражает глубокую благодарность.

Автор выражает благодарность участникам семинара им. Д.Г. Успенского за обсуждение работы, сотрудникам объединения «Уралкалий», ПГО «Севвост-геология», ФГУП «Челябгеолсъемка», ОАО «Баженовская экспедиция», оказавшим помощь в процессе работы и предоставившим материалы для работы.

При решении вышеперечисленных задач проводились теоретические исследования, физическое и математическое моделирование с помощью программ, созданных автором, в том числе с применением математических пакетов "Mathematica" и "Mathcad".

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Виноградов, Владислав Борисович

ЗАКЛЮЧЕНИЕ

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

В работе представлена технология описания геологических поверхностей сплайнами различного вида, в зависимости от характера поверхности. Технология разработана для двухмерного и трехмерного вариантов.

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

Получены выражения для вычисления гравитационного и магнитного полей прямоугольного параллелепипеда, свойства которого описаны полиномом третьей степени.

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

Разработана методика описания изменчивости плотности и намагниченности в трехмерном пространстве и способы решения прямых задач для таких физико-геологических моделей

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

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

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

УСЛОВНЫЕ ОБОЗНАЧЕНИЯ И СОКРАЩЕНИЯ

НДС - напряженно-деформированное состояние геологической среды;

ТД - тензор деформации;

КТД - компоненты тензора деформации;

КТН - компоненты тензора напряжений;

КВС - компоненты вектора смещения;

Ag - поле силы тяжести; u,v,w - компоненты вектора смещения;

Z - вертикальная составляющая вектора магнитной индукции; сту - компоненты тензора напряжений; sy - компоненты тензора деформации;

Е - модуль Юнга; v - коэффициент Пуассона; к - гравитационная постоянная; а - плотность;

X - магнитная восприимчивость; J— намагниченность.

Библиография Диссертация по наукам о земле, кандидата геолого-минералогических наук, Виноградов, Владислав Борисович, Екатеринбург

1. Алексидзе М.А. Решение некоторых основных задач гравиметрии/ М.А. Алексидзе. Тбилиси: Мицниереба, 1985. - 412 с.

2. Аронов В.И. Обработка на ЭВМ значений силы тяжести при произвольном рельефе поверхности наблюдений/ В.И. Аронов. М.: Недра, 1976. - 129 с.

3. Аронов В.И. О некоторых вопросах применения математических методов и ЭВМ в геологии и геофизике/ В.И. Аронов//Геофизический журнал, 1990. - Т. 12,№2. -С. 62-72.

4. Артюшков Е.В. Физическая тектоника/ Е.В. Артющков. М.: Наука, 1993, -454 с.

5. Дилатансионные процессы в консолидированной коре. 1/ Беличенко П.В., П.В., Гинтов О.В., Исай В.М. и др.//Геофизический журнал, 1990, т. 12, №2, -С. 32 44. 2. т.13, №2. - С. 48 - 54.

6. Болотнова Л.А. Методика изучения деформационного состояния геологической среды района Екатеринбурга по гравиметрическим данным: Автореф. дис. . канд. геол.-минерал. наук: 25.00.10/- Екатеринбург, 2007. 23 с.

7. Бор Карл де. Практическое руководство по сплайнам/Карл де Бор. М.: Радио и связь, 1985. - 303 с.

8. Бреббия К. Методы граничных элементов: пер. с англ./ К. Бреббия, Ж. Тел-лес, Вроубел М.: Мир, 1987, - 534 с.

9. Булах Е.Г. Применение метода минимизации для решения задач структурной геологии по данным гравиразведки/ Булах Е.Г., Ржаницын В.А., Маркова М.Н. Киев.: Наукова думка, 1976. - 220 с

10. Булашевич Ю.П. Магнитное поле горизонтального пласта с неоднородным распределением магнитных минералов. УФАН // Геофизический сборник. №2: Труды горно-геологического института. Свердловск. Вып. 30. 1957. - С. 100 — 110.

11. Викторов Г.Г. Некоторые результаты изучения магнитной анизотропии в геодинамических зонах / Викторов Г.Г., Свирина И.Н //Изв. Вузов. Геология и разведка. М., 1982. Деп. №6314-82.

12. Виноградов В.Б. Прямая задача магниторазведки на основе представления поверхности геологических объектов сплайнами / В.В. Филатов, В.Б. Виноградов // Геофизический журнал. 1988. - Т. 10, № 2. - С. 67 - 73.

13. Виноградов В.Б. Определение напряжений и деформаций в земной коре по полю силы тяжести // Геофизические аспекты изучения геологического строения месторождений калийных солей: сб. науч. трудов ВНИИГ. Л., 1989. - С. 162- 167.

14. Виноградов В.Б. Оценка напряженно-деформированного состояния геологической среды по геофизическим данным / В.Б. Виноградов, А.А. Кашкаров // Инф. лист. № 287-90. Свердловск, 1991. - 4 с.

15. Виноградов, В.Б. Методика расчета деформаций среды по гравитационному полю // Известия УГГГА, Вып. 5. Сер.: Геол. и геофизика. Екатеринбург. 1996, С. 110-111.

16. Виноградов В.Б. Магнитные и плотностные интерпретационные сплайн-модели/ В.Б. Виноградов // Известия УГГГА. Вып. 8. Сер.: Геология и геофизика,-1998,—С. 142- 144.

17. Виноградов В.Б. Моделирование угленосных структур по данным магниторазведки / Д.Б. Иванов, В.Б. Виноградов, А.В. Чурсин // Геофизические методы при разведке недр и экологических исследованиях. Вып. 2. Томск: Изд. ТГУ, -1999. С. 212-214.

18. Виноградов В.Б. О методах интерпретации гравитационных и магнитных полей восточноуральских угольных месторождений // Известия УГГГА. Вып. 15. Сер.: Геология и геофизика. 2002. — С. 209 - 214.

19. Виноградов, В.Б. Совершенствование методики обработки данных КМВ на железорудных месторождениях // Геология и полезные ископаемые Западного Урала: материалы региональной научно-практической конференции. Пермь: Пермский гос. ун-т, 2008. С. 224-227.

20. Виноградов В.Б. О методике моделирования железорудных тел //Изв. вузов. Горный журнал. 2008. - № 8. - С. 85-88.

21. Виноградов В.Б. Способ описания поверхности геологических объектов сплайнами// Литосфера. 2009. № 3, - С. 91 - 93.

22. Виноградов В.Б. Моделирование объектов с переменной намагниченностью// Пятые научные чтения памяти Ю.П. Булашевича: материалы. Екатеринбург: ИГф УрО РАН, 2009. С. 74 - 77.

23. Виттке В. Механика скальных пород: пер. с нем. М.: Недра, 1990. 439 с.

24. Гзовский М.В. Основы тектонофизики. М.: Наука, 1975. - 536 с.

25. Гинтов О.В., Исай В.М., Кобылянский В.Б. Механизм разломообразования в коре Украинского щита.// Геофизический журнал, 1991, ч. 1, № 1. - С. 3 - 16, ч. 2, № 2. - С. 35 - 41, ч. 3, №3. - С. 3 - 11, ч. 4, №4. - С. 3 - 16.

26. Гореванов Д.Е., Конорев В.И. Структурно-вещественные преобразования горных пород в зоне Приморского глубинного разлома//Геология и геофизика, 1991.-№5.-С. 54-61.

27. Глазнев В.Н. Оценка напряженного состояния земной коры северо-востока Балтийского щита на основе плотностной модели/ Глазнев В.Н., Маслов Л.А., Комова О.С. // Физика Земли, 1988. №10. - С. 62 - 67.

28. Гудман Р. Механика скальных пород: пер. с англ./ М.: Стройиздат, 1987. -232 с.

29. Ерофеев Л.Я. Магнитное поле и природа аномалий на месторождениях золота. Томск: Изд. Томского у-та. - 1989. - 160 с.

30. Динник А.Н. Статьи по горному делу. М.: Углеиздат. - 1957. - 195 с.

31. Звягинцев Л.И. Деформации горных пород и эндогенное рудообразование. -М.: Наука. 1984.-164 с.

32. Иванов Б.М. Исследования напряженного состояния призабойной части пласта посредством определения плотности угля / Научные сообщения ИГД АН СССР. Т. 2. 1959. - С. 44 - 50.

33. Иванов Н.А., Шапиро В.А. Характеристики динамической намагниченности, создаваемой в образцах естественных ферромагнетиков / Н.А. Иванов, В.А. Шапиро //Изв. АН СССР. Сер. Физика Земли. 1972. - №3. - С. 80 - 86.

34. Калашников А.Г. Магнитные свойства неоднородно намагниченных призм//Известия АН СССР. Сер: геофиз. 1956. - №12. - С. 1369 - 1383.

35. Каплан В.Е. К обобщению и практическому приложению магнитной задачи Булашевича // Известия АН СССР. Физика Земли. 1973. - №5. - С. 102 - 106.

36. Каратаев Г.И. Вопросы теории временных возмущений гравитационного и магнитного полей и движений земной поверхности в связи с современными тектоно-физическими процессами в Земле/ Г.И. Каратаев, В.К. Панкрушин,

37. Комова О.С. Решение прямой трехмерной задачи геомеханики для аномальных плотностных тел/О.С. Комова, Л.А. Маслов, И.К. Туезов //Тихоокеанская геология. 1984. - №3. - С. 103 -105.

38. Комплексные инженерно-геологические исследования при строительстве гидротехнических сооружений, /под ред. Савича А.И. и Куюнджича В.Д. М.: Недра, 1990.-462. с.

39. Кондратьев Б Л. Теория потенциала. Новые методы и задачи с решениями. -М.: Мир, 2007.-512 с;

40. Ластовин И.В. Решение прямой задачи гравиразведки палеточным способом на основе сплайн-аппроксимации источников поля. Днепропетровск: ДГИ. -1983. 13 с. Деп. 482Ук-Д83.

41. Лурье А.И. Теория упругости. М.:. Наука, 1970. - 940 с.

42. Лучицкий А.И. Решение прямой задачи гравиметрии и магнитометрии для многогранника с линейной и квадратичной плотностью/А.И. Лучицкий, Л.В. Гричук // Теория и практика интерпретации гравитационных аномалий. — М., Изд. ИФЗ. 1982. - С. 138 - 203.

43. Магницкий В.А. Основы физики Земли. М.: Геодезиздат. - 1953. - 290 с.

44. Малинин Н.Н. Прикладная теория пластичности и ползучести. М.: Машиностроение, 1975. - 318 с.

45. Марков Г.А. Напряженное состояние пород и горное давление в структурах гористого рельефа / Г.А. Марков, С.И. Савченко JL: Наука. 1984. - 124 с.

46. Мартынова Т.А., Вадковский В.Н. Исследование распределения напряжений в зонах разломов./Сб. Магнетизм горных пород и палеомагнетизм. М., 1969, С. 45-47.

47. Маслов JI.A. Расчет скоростей движений земной коры по аномалиям высот геоида.//Физика Земли. -№11.- 1989. С. 96 - 102.

48. Маслов Л.А. О связи вертикальных перемещений поверхности Земли с гравитационным полем, вызванным внутренними источниками//Тихоокеанская геология. 1982. - №4, С. 58 -64.

49. Маслов Л.А. Модель напряженно-деформированного состояния среды в области проявления гравитационной аномалии: Препринт/ Л.А. Маслов, А.Е. Молчанов // М., ИФЗ АН СССР, 1980, №2, 10 с.

50. Маслов Л.А., Комова О.С. Численное моделирование глубинных геодинамических процессов в активных окраинах. /Физика Земли, №3, 1990 С. 53

51. Матусевич А.В. Объемное моделирование геологических объектов на ЭВМ. -М.: Недра, 1988.- 184 с.

52. Мегеря В.М. Применение детальной гравиразведки в комплексе с другими методами при изучении нефтегазоносности структур на склонах Сургутского и Нижневартовского сводов. Автореф. дис. . канд. геол.-минерал. наук. 04.00.12 Пермь, 1982. - 24 с

53. Миндлин Р., Чень Д. Сосредоточенная сила в упругом полу пространствен/Механика. 1952. - №4(14). - С. 118 - 133.

54. Мохова Е.Н. Прямоугольная призма постоянной восприимчивости в однородном магнитном поле.// Изв. АН СССР. Сер. геоф., -1958. №3. - С. 387 -390.

55. Нагата Т. Магнетизм горных пород: пер.с яп./ М.: Мир, 1965. 346

56. Напряженно-деформированное состояние и устойчивость скальных склонов и бортов карьеров. (Материалы 6 Всесоюзной конференции по механике горных пород.) Фрунзе, изд. Илим, 1979. 402 с.

57. Новоселицкий В.М., Гордин В.М. Построение горизонтальной слоисто-зональной модели плотностного разреза осадочной толщи по данным гравираз-ведки./Сб. "Геофизические изыскания". Вып. 1, 1975. С. 115 128.

58. Нульман А.А. Магнитная восприимчивость сильномагнитных горных пород в условиях гидростатического сжатия: Автореф. дис . канд. ф.-м. наук. Свердловск, 1980.-23 с.

59. Осокина Д.Н. Пластичные и упругие низкомодульные оптически-активные материалы для исследования напряжений в земной коре методом моделировании //. М.: Изд. АН. СССР. 1963, 196 с.

60. Опарин В.Н. К задаче определения ориентации главных напряжений методом скважинной электрометрии. //Геофизические методы контроля напряжений в горных породах. Новосибирск, 1983, С. 3 8.

61. Паркинсон У. Введение в геомагнетизм. М.: Мир, 1986. 528 с.

62. Пахомов М.И., Пахомов В.И. Петрофизический метод выделения и оценки метасоматитов. — М.: Недра. 1988. 152 с.

63. Петрова Г.Н., Юхновец Н.И. Изменение магнитных свойств горных пород в зоне разлома/ Г.Н. Петрова, Н.И. Юхновец //Изв. АН СССР. Сер. геоф., 1953.- №2. -С. 115-123.

64. Петрофизика: Справочник. В 3 кн. Кн. 1. Горные породы и полезные ископаемые/Под ред. Дортман Н.Б. М.: Недра, 1992. - 391 с.

65. То же. Кн. 2. Земная кора и мантия. 286 с.

66. Поля напряжений и деформаций в земной коре: Сб. ст./АН СССР ИФЗ. Отв. Ред. Буланже Ю. М.: Наука, 1987. 183 с.

67. Поля напряжений и деформаций в литосфере. М.: Наука, 1979, - 256 с.

68. Пономарев B.C. Особенности напряженного состояния неравновесной геофизической среды.//Известия Ан СССР, Сер. Физика Земли, 1987, №4, С. 94-97.

69. Ржевский В.В., Новик Г.Я. Основы физики горных пород. 3-е изд, пере-раб. и доп. М.: недра, 1978. - 390 с.

70. Ревякин П.С. Высокоточная магниторазведка/ Ревякин П.С., Бродовой В.В., Ревякина Э.А // М.: Недра, 1986, 272 с.

71. Рыльникова М.В. Геомеханика / М.В. Рыльникова, О.В.Зотеев. М.: Руда и металлы, 2003. - 240 с.

72. Савинский И.Д. Программные системы обработки и интерпретации гравитационных и магнитных данных// Геофизика. 1995. - №1. - С. 24 - 28.

73. Сила тяжести и тектоника. М.: Мир, 1975. - 504 с.

74. Слепак З.М. Применение гравиразведки при поисках нефтеперспектив-ных структур. М.: Недра, 1989. 200 с.

75. Слепак З.М. Применение гравиразведки для изучения нефтегазоносных структур. М.: Недра. 1980. - 152 с.

76. Сначев В.И., Шулькин Е.П., Муркин В.П., Кузнецов Н.С. Магматизм Восточно-Уральского пояса Южного Урала./ВНЦ УРО АН СССР. Уфа, 1990. 179 с.80.* Соловьев F.А. Петрофизическая характеристика эндогенных месторождений. М.: Недра, 1984. - 224 с.

77. Старостин В.И. Геодинамика и петрофизика рудных полей и месторождений.-М:: 1984.-204 с

78. Старостин В.И. Палеотектоническиё режимы и механизмы формирования структур рудных месторождений. -М.: Недра, 1988, 256 с.

79. Стечкин С.Б., Субботин Ю. Н. Сплайны в вычислительной математике. -М.: Наука, 1976.-248 с.

80. Стовас М.В. Избранные труды. Ч. 1 -М.: Недра, 1975. 156 с.

81. Страхов В.Н. Основные направления теории и методологии интерпретации геофизических данных на рубеже XXI столетия // Геофизика. -1995. -№3.-С.9- 18.

82. Страхов В.Н. Об общих решениях обратных задач гравиметрии и магни-тометрии//Известия вузов. Геология и разведка. №4. - 1978. - С. 105-117.

83. Страхов В.Н. Некоторые вопросы интерпретации геомагнитных измерений в океане./ Магнитные аномалии океанов и новая глобальная тектоника. М.: Наука, 1981, С. 20-59.

84. Теркот Д., Шуберт Дж. Геодинамика: Геологические приложения физики сплошных сред. Ч. 1: Пер. с англ. М.: Мир, 1985. - 376 с.

85. Тимофеев А.Н. Об интерпретации магнитных аномалий в случае изменяющейся магнитной восприимчивости горных пород. /Вопросы разведочной геофизики. Вып. 1. 1960. - С. 137 -146.

86. Туезова Н.А. Физические свойства горных пород Западно-Сибирской нефтегазоносной провинции/ Н.А. Туезова,Л.М. Дорогницкая, Р.Г. Демина, Н.И. Брюзгина //М.: Недра, 1975. 184 с.

87. Турчанинов И.А. Основы механики горных пород./ И.А.Турчанинов, М.А. Иофис, Э.В. Каспарян -Л.: Недра, 1977. 503 с.

88. Троян В.Н. Методы аппроксимации геофизических данных на ЭВМ/ Тро-ян В.Н., Соколов Ю.М. Л.: Изд. ЛГУ, 1989. - 304 с.

89. Трубицын А.П. Упругие напряжения связанные с неровностями плотно-стных границ раздела в Земле // Физика Земли. 1979. - №12. - С. 15-22.

90. Трухин В.И., Жиляева В.А., Зинчук Н.Н. и др. Магнетизм кимберлитов и траппов./ В.И. Трухин, В.А. Жиляева, Н.Н. Зинчук и др. М.: Изд. МГУ, 1989

91. Тяпкин К.Ф. Изучение разломных структур докембрия геолого-геофизическими методами. Киев: Наукова думка, 1986. - 166 с.

92. Ферхуген Дж., Тернер Ф., Вейс Л. и др. Земля Введение в общую геологию: Пер. с англ. М.: Мир,, т. 1, 1974. 392 с.

93. Физические свойства горных пород и полезных ископаемых (петрофизи-ка. Справочник геофизика. М.: Недра, 1984. - 455 с.

94. Филатов В.В., Кузнецов Н.С. Применение гравиметрии для тектоно-физического анализа / В.В. Филатов, Н.С. //Известия вузов. Геология и разведка. 1989. - №4. - С. 97 - 102.

95. Филатов В.В. Магниторазведка: магнитная восприимчивость, индуцированная и остаточная намагниченности: научное издание. Екатеринбург: Изд-во УГГУ, 2008. 236 с.

96. Филатов В.В. Теоретические основы метода оценки деформаций по гравитационному полю/Сб. научн. Трудов ВНИИТ, Геофизические аспектыгеологического строения месторождений калийных солей., JI. 1989, С. 155 -161.

97. Филатов В.В. Теория и практика геодинамического анализа гравитационного поля (на примере рудных районов Урала): Дис. . д-ра геол.-минерал, наук: 04.00.12. Свердловск, 1990. - 377 с.

98. Фотади Э.Э. К теории временных возмущений гравитационного и магнитного полей в связи с современными тектоно-физическими процессами в земле/ Э.Э. Фотади, Г.И. Каратаев, В.И. Щеглов //ДАН СССР, 1966. т. 171, №3. - С. 590-592.

99. Цок Н.О. Решение обратной задачи гравиразведки для контактной поверхности аппроксимированной с помощью сплайнов// Докл. АН УССР, 1982.- №8.-С. 31-33.

100. Шолпо JI.E. Использование магнетизма горных пород для решения геологических задач. Д.: Недра, 1977. - 182 с.

101. Fajklewicz Z. Gebirgslage und grawimetrische Mikroanjmalien./Freib/ Forsch. H. С 368/ Geowissenschaften/ Geophysik. 1981. S. 163 166/

102. Pedersen L.B., Rasmussen T.M. The gradient tensor of potential fielid anomalies: Some implications on data collection and data processing of maps. 55. 1558-1556. 1990.1. Фондовая литература

103. Филатов B.B., Виноградов В.Б. Совершенствование методики оценки напряженно-деформированного состояния среды по геофизическим данным. Отчет по теме 62-204-89./ Свердловский горный институт; Руководитель В.В. Филатов; Свердловск. 1990.

104. Ольховников Ю.П. Физико-механичесике характеристики пестроцвет-ной известково-песчанистой толщи / ПО «Уралкалий». Пермь. 1977.

105. Гайнанов Ш.Х. Определение физико-механических свойств пород на шахтных полях ПО "Уралкалий" в процессе геологоразведочных работ / УФ ВНИИГа. Пермь. 1986.

106. СОГЛАСОВАНО ' РекюрСпроректорЭвузаj„---'----ядатель1. Ь №ф.Багаутинов Г.А.

Информация о работе

Виноградов, Владислав Борисович
кандидата геолого-минералогических наук
Екатеринбург, 2009
ВАК 25.00.10

Диссертация

Методика физико-геологического моделирования объектов с переменной плотностью и намагниченностью - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы