New treatment of interfaces for 3D seismic wave propagation problem using generalized finite difference method

Local surface sedimentary structures are made up of stacking strata. Layers of rock or sediment characterized by certain lithological properties that distinguish them from adjacent layers, separated by bedding planes here called interfaces. These strata are typically parallel, but their geometry can be complex, with inclined layers of variable thicknesses produced by tectonic movements or erosion processes. In this work, a formulation is proposed for the treatment of interfaces in seismic wave propagation problems with 3D domains using generalized finite difference method. This capacity is of great relevance, since it is the complexity of the soil structures what makes this meshless method interesting in comparison to the computationally efficient finite difference method. To validate the proposal, a set of examples are solved and their results are compared to those derived from analytical expressions or through finite element method models.