Аннотация:The presentation will provide an overview of the problems of mathematical modeling ofgeomechanical processes occurring in rocks during the exploration and development of reservoirs [3,4, 16]. These problems were solved using CAE Fidesys software [8, 11, 24], which allows solvingstatic and dynamic problems of geomechanics and geophysics using finite (FEM) and spectral (SEM)element methods [9, 10, 23] of a variable approximation order in space at non-conformal unstructuredmeshes [17-21]. The mathematical formulation is based on the theory of repeated superposition oflarge deformations [13-15]. A numerical discretization of the posed boundary problems of interactingsolids [12, 22, 23] is performed using a discontinuous spectral element method and multi-pointcontstraints [1, 2, 5-7] at non-matching mesh interfaces between interacting solid rock structures.Several industrial applications of the developed approach are considered. Seismic wavepropagation in the heterogeneous media with initial geomechanical stresses is considered. A modellingof an induced anisotropy is performed by the superposition of dynamic deformations onto initialgenerally finite strains. Use of variable order spectral elements at non-conformal meshes allows one tosimplify the process of unstructured mesh generation for the discretization of complex geologicalmodels and to set the local spatial order of the SEM discretization depending on the speed of seismicwaves in geological structures, which significantly reduces the computational costs when conductingnumerical modeling and lowers the requirements to the model preprocessing and mesh quality.Geomechanical simulation of the multistep drilling of the wellbore taking into account a historydependence of the strains due to the sequential origination of stress concentrators as well asmechanical properties change is presented. Using SEM, we simulated stresses and strains of the nearwell bore poroelastoplastic zone taking into account the dynamic change in the porepressure field ofthe saturating fluid. Assessment of deformation of the wellbore depending on the pressure of thedrilling fluid inside it expands the possibilities of applying geomechanical modeling and allows one toclarify recommendations for the drilling parameters. The effect of the development of the plasticregion is the cause of the non-uniform and, in real conditions, non-linear distribution of themechanical and filtration-capacitive properties of the rock. Physical processes resulting from thiseffect, in turn, affect well productivity. A geomechanical model of the near-wellbore zone, built for aspecific task, can reduce risks during drilling and trajectory, assess potentially promising areas fordrilling and stimulation, and solve the problems of wellbore stability.The considered approach allows predicting in more detail the behavior of the rock duringreservoir development, taking into account different stages of the field deformations. In particular, theredistribution of accumulated deformations during multistep loading and / or changes in the structure(topology) of the loaded body, as well as contact conditions of adhesion / sliding at the interlayerboundaries and bonded contacts are taken into account.The reported study was partially funded by Russian Science Foundation project № 19-77-10062and by RFBR and Moscow city Government according to the research project № 19-38-70001.