Математическое моделирование нового метода обработки сейсмосигнала на примере разреза Саратовского Правобережья. Часть 2

Abstract

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

This paper presents a new method of seismic exploration with two buried receivers on the same vertical, which helps in comparing the experimental seismogram cleared of non-vertical signals with the calculated theoretical seismogram, which depends on the parameters of the layered structure, such as the thickness of the layers, their density, velocity of propagation of longitudinal acoustic waves in rock, and seismic quality factors of media. The method of least squares aided in the search for the global minimum of the discrepancy function of the calculated and experimental seismic pulse waveforms. The implementation of the algorithm for solving the inverse seismic problem was carried out in the MatLab software package. The construction of a theoretical seismogram includes the calculation of the spectrum of the input impulse, calculation of the frequency dependences of the reflection coefficients and the passage of waves at the interfaces of rocks with different physical parameters, the calculation of the shape of the seismic impulses received by the receivers. A comparison of the calculated and real parameters of the geological structure obtained in the COMSOL Multiphysics software package shows that the developed algorithm results in obtaining a seismogram that practically coincides with the experimental seismogram corresponding to the real rock section, with appropriate geological parameters of media, such as thickness, density and propagation velocity of acoustic waves, as well as quality factors.