Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

STUDY OF CONVERGENCE OF ALGORITHMS AND COMPUTATIONAL SCHEMES OF METHOD OF GEOLOGICAL ENVIRONMENT INTERPOLATION IN A FRAGMENTED SET OF DATA

Kobrunov A.I. 1 Motryuk E.N. 1 Lominskiy D.O. 2
1 FGBOU VPO «Ukhta State Technical University»
2 FGBOU VPO «Russian State University of Oil and Gas named after Gubkin»
A study of the convergence of the algorithm-sided interpolation to restore the bulk medium system profiles, which allows you to take into account as much as possible a priori information about the geological and geophysical parameters of the medium. The convergence of algorithms and computational schemes for constructing the volume density of structural and geological environment model in model examples is investigated. The step of boundaries interpolation should be chosen so that it could take into account the characteristics of the geological environment with enough accuracy for solutions of the following tasks as well as close to the step to which the original data is set. The lower boundary of the step size depends on the type of work. Also, the choice of calculated lines within the area significantly influences the result of constructing. It is determined by the position of the anomaly under consideration. The estimated line has to pass through the central part of the anomaly in question transversely. Taking into account the comments made, the interpolation process is stable and can be used for modeling geological environments with complex structure.
gravity prospecting
complex structure environment
direct gravity problem
interpolation method of boundaries
geological environment model
convergence of the algorithm
1. Aminov L.Z., Kobrunov A.I., Moiseenkova S.V., Motrjuk E.N., Shilova S.V., Muzhikova A.V. Metodika integrirovannoj interpretacii gravimetricheskih dannyh v uslovijah slaboj izuchennosti s celju postroenija obemnyh regionalnyh plotnostnyh modelej sedimentacionnyh bassejnov / Geologija i mineralnye resursy Evropejskogo Severo-vostoka Rossii: materialy XIV geologicheskogo s#ezda Respubliki Komi. Tom IV. Syktyvkar: Geoprint, 2004. рр. 79–81.
2. Kobrunov A.I. Algoritm uchjota vlijanija vmeshhajushhej sredy pri inversii gravimetricheskih dannyh / A.I. Kobrunov, M.I. Barabanov // Izvestija Komi NC UrO RAN. 2014. no. 1 (17). рр. 84–88.
3. Kobrunov A.I. Matematicheskie metody modelirovanija v prikladnoj geofizike (izbrannye glavy) v 2-h chastjah. Ch.1 Funkcionalno-analiticheskie osnovy // Mezhdunarodnyj zhurnal jeksperimentalnogo obrazovanija. 2014. no. 11, ISSN 199-3947. 78 р. Ch. 2 Sistemnyj analiz i modelirovanie v uslovijah neopredelennosti. Uhta, UGTU, 2014. 154 р.
4. Kobrunov A.I., Petrovskij A.P., Danilenko A.N., Moiseenkova S.V., Aminov L.Z. Teorija i metody kolichestvennoj kompleksnoj interpretacii geofizicheskih dannyh // Aktualnye nauchno-tehnicheskie problemy razvitija geologo-geofizicheskih promyslovyh i poiskovo-razvedochnyh rabot v Respublike Komi. Kniga 1. Uhta, 2001. рр. 261–316.
5. Kobrunov A.I. Matematicheskie osnovy teorii interpretacii geofizicheskih dannyh: ucheb. posobie. Uhta: UGTU, 2007. 286 р.: il.
6. Motrjuk E.N., Muzhikova A.V. Ob#emnoe strukturno-plotnostnoe modelirovanie geologicheskoj sredy metodom odnostoronnej interpoljacii na primere tel prostejshej formy // Severgeojekoteh 2003: materialy IV Mezhregionalnaja molodezhnaja nauchnaja konferencija. Uhta, UGTU, 2003. рр. 554–556.
7. Motrjuk E.N., Muzhikova A.V., Shilova S.V. Metodika vybora approksimacii dlja reshenija prjamoj zadachi gravirazvedki // Severgeojekoteh 2002: tezisy dokladov Mezhregionalnaja molodezhnaja nauchnaja konferencija. Uhta, 2002. рр. 31–32.

Прогноз нефтегазоносности требует наличия объемной структурно-плотностной модели геологической среды, удовлетворяющей гравитационному полю. Однако для построения имеются только фрагментарно заданные данные. Построенные двумерные модели могут быть реконструированы решением структурной или плотностной обратных задач гравиразведки [1–5]. Для преобразования двумерных данных в трехмерные существует много прикладных программных средств (например, Surfer, Коскад 3D и др.), алгоритмов (послойной интерполяции, методы крайгинга и др.), выбор которых зависит от характера имеющихся данных. К тому же для реальных геологических сред погрешности, лежащие в основе этих продуктов, могут быть велики и сильно влиять на качество модели. В данной работе представлена апробация алгоритмов и вычислительных схем односторонней интерполяции профильных данных в объемные плотностные модели для случая сложнопостроенных сред. Их основные положения изложены в [6, 7]. Для установления точности вычислительных схем разработанной интерполяции геологических границ написана программа GRAST на языке Си++.

Модель среды

Модель среды представляет собой призму слоистой структуры, состоящую из трех пластов с плотностями σ0 = 2,6 г/cм3, σ1 = 3,0 г/cм3 и σ2 = σ0 = 2,6 г/cм3. Источником гравитационной аномалии является второй пласт, включающий в себя шар радиуса R = 2500 м, расположенный в центре призмы. Нижняя граница призмы korbunov01.wmf, верхняя korbunov02.wmf (рис. 1). Тестирование интерполяции проведено на трех вариантах расположения профилей, в статье приведен один вариант. В природе не встречаются такого вида аномалии, чаще всего они типа уступа, так как геологическая среда слоистая ввиду своего вида образования. Значит, в нашем случае все вычисления будут давать большее значение погрешности построения, чем для реального объекта. Шаг интерполяции призмы step_x = step_y принимает значения 100, 250, 300, 400, 400, 600, 800, 1000 м. В результате получается пространственная сетка K = {(x, y, z): x, y ∈ (10000, 30000), z ∈ (1, 10000)}. Шаг, с которым заданы значения границ по профилям step_r = 250 м.

Модель задачи

При расчете значений глубин использовались формулы, выведенные для односторонней интерполяции границ в [6, 7]. Для вычисления точных значений глубин на сфере используем уравнение

korbunov03.wmf

и уравнение прямой korbunov04.wmfв противном случае (вне сферы).

Для характеристики полученной при помощи односторонней интерполяции границы рассмотрим среднеквадратическое отклонение точных и рассчитанных значений при разном шаге интерполяции step_x:

korbunov05.wmf (1)

где Kx + 1, Ky + 1 – количество координат x и у сетки K.

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

Модель

Задано восемь тестируемых профилей, занумерованных 1–8 (рис. 2). Особенностью расположения профилей является то, что 7 и 8 профили имеют диагональную направленность, а остальные параллельны линиям сетки интерполяции.

Граница, представленная на рис. 3, не является симметричной относительно центра поверхностью, поскольку не все заданные координаты профилей совпадают с узлами сетки K. Граница имеет более округлый вид, чем в моделях 1 и 2, потому что в данном случае используется больше геолого-геофизической информации. В этом случае погрешность вычисленной границы (рис. 7–8) равна 2,3 м, что составляет 0,09 % от амплитуды рассматриваемой границы. Граница, представленная на рис. 4, не является симметричной относительно центра поверхностью, поскольку не все заданные координаты профилей совпадают с узлами сетки K. В этом случае погрешность вычисленной границы (рис. 3–4) равна 2,7 м и составляет 0,09 % от амплитуды рассматриваемой границы.

В случае, представленном на рис. 5, погрешность вычисленной границы (см. рис. 7–8) равна 6,4 м и составляет 0,26 % от амплитуды рассматриваемой границы. Граница, представленная на рис. 6, не является симметричной относительно центра поверхностью, поскольку не все заданные координаты профилей совпадают с узлами сетки K и к тому же очень велик шаг интерполяции.

Таким образом, для Модели 2 наилучшим шагом интерполяции для step_r = 250 м является шаг step_x до 800 м (что в 3,2 больше step_r), для step_r = 500 м – step_x до 800 м (что в 1,6 больше step_r). Далее отклонение от аппроксимирующей график невязки тренда-функции будет возрастать и иметь непредвиденные скачки.

pic_21.tif

Рис. 1. Тестируемая модель

pic_22.tif

Рис. 2. Сечение Модели 3 плоскостью z = 5000 м

pic_23.tif

Рис. 3. Результаты интерполяции второй границы Модели 3 при step_x = 250

pic_24.tif

Рис. 4. Результаты интерполяции второй границы Модели 3 при step_x = 300

pic_25.tif

Рис. 5. Результаты интерполяции второй геологической границы Модели 3 при step_x = 800

pic_26.tif

Рис. 6. Результаты интерполяции второй геологической границы Модели 3 при step_x = 500, step_r = 1000

pic_27.tif

Рис. 7. График зависимости среднеквадратической погрешности исходной и интерполируемой границ от step_x при step_r = 250 м

pic_28.tif

Рис. 8. График зависимости среднеквадратической погрешности исходной и интерполируемой границ от step_x при step_r = 500 м

Из графиков, представленных на рис. 7, 8, невязка в 5 мГал достигается при step_x = 670 м для step_r = 250 м, при step_x = 400 м для step_r = 500 м. При увеличении числа профилей в 1,3 раза фиксированное значение невязки достигается при шаге step_x, большем в 3,35 (при step_r = 250 м) и 2,66 (при step_r = 250 м) раза относительно исходного. Этот факт показывает, насколько важно для получения наиболее близкой к реальной среде модели использовать как можно больше геолого-геофизической информации.

Проведенные выше исследования позволяют сказать, что шаг интерполяции границ следует выбирать таким, чтобы он мог учитывать особенности геологической модели с достаточной для решения последующих задач точностью. Шаг, с которым происходит интерполяция исходных данных, лучше брать близким к шагу, с которым исходные данные задаются для интерполяции, до увеличения в 1,6 раз. Нижняя граница выбора step_x зависит от вида работ. Например, для детальных допустим самый маленький шаг step_x – 50 м, что для региональных исследований применять нецелесообразно. Существенно влияет на результат построения также выбор расчетных линий в пределах участка. Он определяется положением интересующей нас аномалии. Расчетная линия должна проходить через центральную часть наиболее интересующей нас аномалии, вкрест ее простирания. С учетом сделанных замечаний процесс интерполяции является устойчивым и может применяться для моделирования геологических сред сложного строения. Таким образом, разработанная интерполяционная процедура может использоваться для объемного моделирования осадочных бассейнов. Созданные способы восстановления объемной среды по системе профилей (односторонняя, «змейка», рекурсия) позволяют максимально учитывать априорную информацию о геолого-геофизических параметрах среды.

Рецензенты:

Бурмистрова О.Н., д.т.н., заведующая кафедрой технологии и машин лесозаготовок, ФГБОУ ВПО «Ухтинский государственный технический университет», г. Ухта;

Андронов И.Н., д.т.н., профессор кафедры сопротивления материалов и деталей машин, ФГБОУ ВПО «Ухтинский государственный технический университет», г. Ухта.