Oleg Nechaev

Due to the ongoing global warming on the Earth, permafrost degradation has been extensively taking place, which poses a substantial threat to civil and industrial facilities and infrastructure elements, as well as to the utilization of natural resources in the Arctic and high-latitude regions. In order to prevent the negative consequences of permafrost thawing under the foundations of constructions, various geophysical techniques for monitoring permafrost have been proposed and applied so far: temperature, electrical, seismic and many others. At the same time, non-stationary electromagnetic methods seem to have significant potential and a number of important practical advantages. We propose a cross-borehole exploration system for the highest localization of target objects in the cryolithozone. A novel mathematical apparatus for three-dimensional modeling of transient electromagnetic (TEM) signals by the vector finite element method has been developed. The original combination of the latter, the Sumudu integral transform and artificial neural networks makes it possible to examine spatially heterogeneous objects of the cryolithozone with a high contrast of geoelectric parameters, significantly reducing computational costs. We have studied the sensitivity of the TEM signals to the boundary between thawed and frozen rocks depending on the inter-borehole distance and spatial orientations of the transmitters and receivers. Numerical simulation results of the TEM monitoring of industrial facilities located on permafrost are considered. The formation of a talik has been shown to significantly manifest itself in the measured electromagnetic responses, which enables timely prevention of industrial disasters and environmental catastrophes.

An algorithm for the numerical simulation of spontaneous potential (SP) signals based on the finite element method is proposed. The SP signals are passively recorded electrical signals resulting from the physical and chemical processes of various nature occurring inside the well and adjacent formations. The proposed algorithm allows one to simulate the SP signals for the case of axial symmetry of the studied area and an electrochemical SP source. Using sigmoidal functions to describe the distribution of petrophysical parameters in the study area allows one to obtain an analytical expression for the density of the SP source. The correctness of the proposed algorithm is verified on a simplified three-layer model. The dependence of the obtained SP signals on the thickness of the permeable formation and on its salinity is demonstrated.

(1) This article is devoted to the development of a theoretical and algorithmic basis for numerical modeling of the spontaneous potential method (SP) as applied to the study of sandy-argillaceous reservoirs. (2) In terms of coupled flows, we consider a physical–mathematical model of SP signals from an electrochemical source, with regards to the case of fluid-saturated shaly sandstone. (3) An algorithm for 2D finite-element modeling of SP signals was developed and implemented in software, along with its internal and external testing with analytical solutions. The numerical SP modeling was carried out, determining the dependences on the reservoir thickness and porosity, the amount of argillaceous material and the type of minerals. We performed a comparative analysis of the simulated and field SP data, using the results of laboratory core examinations taken from wells in a number of fields in the Latitudinal Ob Region of Western Siberia. (4) The results of the study may be used either for the development of the existing SP techniques, by providing them with a consistent computational model, or for the design of new experimental approaches.

This paper is dedicated to the topical problem of examining permafrost’s state and the processes of its geocryological changes by means of geophysical methods. To monitor the cryolithozone, we proposed and scientifically substantiated a new technique of pulsed electromagnetic cross-well sounding. Based on the vector finite-element method, we created a mathematical model of the cross-well sounding process with a pulsed source in a three-dimensional spatially heterogeneous medium. A high-performance parallel computing algorithm was developed and verified. Through realistic geoelectric models of permafrost with a talik under a highway, constructed following the results of electrotomography field data interpretation, we numerically simulated the pulsed sounding on the computing resources of the Siberian Supercomputer Center of SB RAS. The simulation results suggest the proposed system of pulsed electromagnetic cross-well monitoring to be characterized by a high sensitivity to the presence and dimensions of the talik. The devised approach can be oriented to addressing a wide range of issues related to monitoring permafrost rocks under civil and industrial facilities, buildings, and constructions.

The paper is dedicated to the topical problem of examining permafrost state and the processes of its geocryological changes by means of geophysical methods. To monitor the cryolithozone, we propose and scientifically substantiate a new technique of pulsed electromagnetic cross-well sounding. Based on the vector finite-element method, we created a mathematical model of the cross-well sounding process with a pulsed source in a three-dimensional spatially heterogeneous medium. A high-performance parallel computing algorithm was developed and verified. Through realistic geoelectric models of permafrost with a talik under a highway, constructed following the results of electrotomography field data interpretation, we numerically simulated the pulsed sounding on the computing resources of the Siberian Supercomputer Center of SB RAS. The simulation results suggest the proposed system of pulsed electromagnetic cross-well monitoring to be characterized by a high sensitivity to the presence and dimensions of the talik. The devised approach can be oriented to addressing a wide range of issues related to monitoring permafrost rocks under civil and industrial facilities, buildings and constructions.

In this article, we are the first to formulate the direct and inverse problems of resistivity logging on determining the components of the electrical resistivity tensor of rocks from a set of high-frequency induction and lateral logging sounding measurements. Using a finite element approximation, high-order hierarchical basis functions, computationally efficient multilevel methods and a multistart algorithm with the DFO-LS local optimization method, we investigate the capability of reconstructing the horizontal and vertical resistivity components, as well as the tilt of the resistivity tensor principal axes with regard to the study of complex geological objects. A separate consideration is given to a realistic generalized geoelectric model of the unique hydrocarbon source with hard-to-recover reserves, the Bazhenov Formation.