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

Abstract

Methods for creating digital grid models of geological surfaces based on approximation by bicubic B-splines are widely used in solving problems of mathematical geology. The variational-gridding method of geological mapping is a flexible and powerful tool that allows to use a large amount of source data for mapping, as well as a priori information about the spatial distribution of the mapped parameter. The smoothness of the basis functions does not allow the direct use of this effective method for mapping geological surfaces complicated by faults. This fact requires adaptation of the variational-gridding method to mapping with faults. The article discusses the technology of separate construction of the fault and plicative (smooth) components of the structural map. The fault component is represented in the form of an antiplane shear field of an elastic membrane described in a stationary two-dimensional Laplace equation. The fault network is modeled by narrow contours, at the boundaries of which the values of tectonic displacements are set. The Laplace equation is solved by the method of boundary integral equations, which allows one to calculate the displacement field at an arbitrary point of the mapped region, as well as to most accurately approximate the complex geometry of faults. Modeling of the plicative component of the structural surface takes place on the basis of a variational-gridding approach with adjustment of tectonic displacements. The approach combines the advantages of the spline approximation method and the accuracy of the semi-analytical solution for the fault component. It does not impose restrictions on the geometry of faults, and also allows for more efficient mathematical operations with structural surfaces, and their differential and integral characteristics.

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