Вы здесь

Об'ємне комп'ютерне моделювання складних рудних покладів за комплексом геолого-геофізичних даних

Автор: 
Нікулін Сергій Леонідович
Тип работы: 
Дис. канд. наук
Год: 
2002
Артикул:
3402U001611
99 грн
(320 руб)
Добавить в корзину

Содержимое

раздел 2).
В. Скважинная информация должна быть организована в массив K, как описано в
разделе 2 (см. рисунок 2.1а).
Г. В случае, если в какой-то точке скважины отсутствуют значения некоторых
дополнительных характеристик, они должны быть заменены определённым кодом
отсутствия данных (например, 9.9Е+9). Если в точке не проводились измерения
моделируемого параметра V, то эта точка должна быть исключена из дальнейшего
рассмотрения.
Д. В случае, если используется одно из гравитационных полей, в число t
дополнительных характеристик должна входить плотность; при использовании
магнитных полей – модуль и направление намагниченности.
E. Если в некоторой точке скважины отсутствуют замеры требуемых плотностных или
магнитных свойств, они должны быть оценены (спрогнозированы) любым известным
способом.
Простейший из них заключается в следующем. Если в точке Р известен тип породы
V(P), но отсутствует замер, например, плотности d, то для оценки последней
можно рассчитать среднюю плотность в точках, для которых V=V(P), и присвоить её
d(P).
Несколько более сложный вариант этого метода предполагает перемножение
плотностей на весовые коэффициенты, обратно пропорциональные евклидовым
расстояниям до точки P.
Если определения типов пород не производились, для оценки может быть
использована другая характеристика (или их набор), предположительно связанная
некоторой, максимально тесной, зависимостью с плотностью; кроме того, может
быть использован математический аппарат распознавания образов.
Ж.  Внешние потенциальные поля должны быть, по возможности, освобождены от
влияния региональной составляющей одним из известных способов [144]; также
желательно максимальное подавление помехи.
Далее для простоты будет рассматриваться работа только с гравитационным полем
(кроме случаев, указанных отдельно).
4.2. Построение блочных моделей
выполняется в 5 этапов.
1. Вначале для каждой из Г точек измерений, формирующих массив K, находятся M
«наилучших» ссылок с наименьшими значениями F. Последние рассчитываются по
формуле (2.1). Процесс нахождения ссылок детально описан в разделе 2. В
расчётах могут участвовать значения не только моделируемого параметра V, но и
любых других, измеренных или вычисленных в точках скважин.
В результате для каждой точки измерения Ki, i=1,…, Г, получаем набор из M
«наилучших» ссылок, соединяющих её c точками, которые с наибольшей вероятностью
принадлежат тому же самому геологическому телу, что и Ki (см. рисунок 4.1б);
2. Вводится трехмерная решётка (регулярный 3D грид) G размером q = I Ч J Ч L
ячеек, представляющая собой «каркас» строящейся блочной модели. Для центра
каждой ячейки (далее называемого «узлом») в результате будут определены
значение параметра V, а также степень надёжности этого значения. Размеры
единичной ячейки решетки определяются, исходя из необходимой детальности
изучения геологической среды (чем меньше ячейка, тем выше детальность).
Значения I, J и L выбираются таким образом, чтобы создаваемая 3D решетка
покрывала всю область исследования. При этом необходимо учитывать
технологические трудности, которые могут возникнуть в случае слишком большого
значения q  (существенные объёмы компьютерной памяти, необходимые для хранения
промежуточных и конечных результатов, значительное время для расчёта блочной
модели и манипуляций с нею).
Далее для каждого узла Р решетки G выполняются следующие действия:
- для каждой из M ранее найденных ссылок вычисляются значения Т по формуле
(2.2);
- из H ссылок (Hі1) с наименьшими значениями Т формируется массив ссылок Q,
характеризующий данный узел интерполяции. Концы этих ссылок представляют собой
точки измерений в скважинах, с наибольшей степенью вероятности находящиеся
внутри того же геологического тела, что и текущий узел решетки G. Число H
вводится интерпретатором; его увеличение приводит к получению более усреднённых
результатов.
Повторение этих двух операции для всех узлов G позволяет получить множество
массивов ссылок , где q = I Ч J Ч L.
В результате для каждого узла P решетки G получаем набор из H ссылок,
соединяющих узел с точками скважин, которые с наибольшей степенью вероятности
принадлежат к тому же самому геологическому телу, что и P (см. рисунок 4.2а;
при H=1).
3. Выполняется построение блочной плотностной модели путём вычисления плотности
в узлах G. Для этого в каждом узле Р из массива Q(P) выбирается ссылка с
наименьшим значением T, далее называемая «наилучшей» (рис. 4.1г), и вычисляется
значение плотности по формуле
, (4.1)
где А и B – концы ссылки (см. рис. 2.6); d(A), d(B) – значения плотности в
точках А и B. Найденные плотности приписываются окрестностям узлов решетки G.
В результате выполнения этого этапа создаётся блочная плотностная модель,
представляющая собой совокупность значений плотности, приписанных узлам
трёхмерной решётки G. Внешнее поле от построенной модели практически наверняка
не будет согласовано с измеренным. Увязка полей производится на следующем
этапе. Основные принципы создания плотностных моделей изложены в разделе 3.
Следует отметить, что вместо плотностной модели может строиться магнитная, или
же обе они одновременно; однако здесь, без потери общности, рассматривается
случай работы только с плотностной моделью.
4. Целью этапа является нахождение новых «наилучших» ссылок для узлов решетки
G, по которым в дальнейшем будет выполняться интерполяция значений
характеристики V. Ссылки выбираются из массивов ссылок Q таким образом, чтобы в
результате поле от построенной блочной плотностной модели было согласовано с
измеренным внешним гр