Research and Modeling of Nonlinear Acoustic Processes in a Layered Nonlinear Medium with a Porous Fluid-Saturated Inclusion of a Hierarchical Type

Olga Hachay¹, Veniamin Dryagin², Andrey Khachay³

¹Institute of Geophysics, Ural Branch of Russian Academy of Sciences, Yekaterinburg, Russia.
²LLC “NPF Intensonik”, Yekaterinburg, Russia.
³Ural Federal University, Yekaterinburg, Russia.

DOI: 10.4236/ojg.2019.99032                

Abstract

Problem statement: The results of the study of seism acoustic emission arising in a porous two-phase geological environment under acoustic influence are presented. Acoustic emission arising in reservoirs of oil fields using good observations is considered. The regularity of the emission processes of acoustic emission, which manifests itself in the form of discrete spectra of signals similar to oscillations of nonlinearly coupled oscillators, is shown. Spectra have special characteristics for each type of rock. Applied method and design: An algorithm for modeling the process of resonant acoustic response of a porous fluid-saturated reservoir with hierarchical structure and plastic properties on acoustic frequency excitation is developed. That algorithm is developed as an iterative process for the solution integral and integral-differential equations. The frequencies that are parameters of the direct problem are used from the spectra of observed data of acoustic emission in the oil wells. Typical results: For the first time, it had been found the relation between resonant frequencies of the acoustic emission and plastic properties, these values of frequencies had been used in the algorithm of modeling distribution of longitudinal waves in the fluid saturated nonlinear plastic environment. Concluding note (Practical value/implications): The analysis of these emission processes can serve as a source of information about the filtration-capacitive properties of productive reservoirs of a porous type with a hierarchical structure. It is used by practical data of oil fields of Western Siberia.

Keywords

Seismic Emission after Acoustic Impact, Reservoir of Hierarchical Structure with Plastic Properties, Algorithm of 2D Modeling Connected with Energy Spectrograms

1. Introduction

Using the method of Acoustic emission (AE) research, it is possible to obtain information about the processes occurring inside the substance of the reservoir in the well. It has an important advantage over other geophysical methods due to its high resolution in the frequency range. For understanding processes in a saturated porous medium and building adequate models, it is not enough to take into account only the linear relationship between the AE and dynamic processes in the medium, which is a nonlinear system. The principal feature of a nonlinear medium is the interaction of a limited number of waves associated with the conditions of frequency resonance and wave vectors [1] . The study of acoustic emission in samples of sandstones of terrigenous reservoirs under the action of mechanical stresses was carried out at the UIK-AE core research unit [2] . The energy of the acoustic emission signal, which is released in the process of acoustic impact, was determined by the method of calculating the spectral energy density in the whole recorded frequency range from units of hertz to 20 kHz by the Intel Graph program. In accordance with the theory of nonlinear systems, when nonlinearly coupled oscillators interact, processes are possible when the system generates oscillations with combinational frequencies. In this case, the condition of resonance of frequencies is fulfilled and it is indicated that energy can be exchanged between oscillators, in which the energy of the excited high-frequency oscillator is transferred to two low-frequency oscillators or the reverse process, the merging of low-frequency oscillations. Such an interaction has, for example, a place in the system of three nonlinearly coupled oscillators described by a system of equations [1] .

${\ddot{x}}_i+{\omega }_i^2{x}_i=\mu {\alpha }_i{x}_j{x}_k;i=1,2,3;j,k\ne i$

For a small value of µ, this is a weekly nonlinear system. Its behavior is close to the superposition of quasiharmonic oscillations of oscillators with slowly varying amplitudes. Due to the nonlinear coupling, oscillations of two oscillators with frequencies ω₁ and ω₂ generate oscillations with a combination of frequencies ω₁ and ω₂ in the system. The action of a small nonlinearity accumulates if the frequency resonance condition is fulfilled:

${\omega }_1+{\omega }_2={\omega }_3$

When sandstone is deformed, which is an inhomogeneous medium characteristic of the oil reservoir, and when an acoustic field is applied, the formation of wave processes takes place in it. The resulting acoustic emission is represented by both a discrete, pulsed, and continuous noise-like component [3] [4] .

2. Nonlinearity of Caused Acoustic Emission Investigation in the Oil Well

The field of acoustic emission and reception of it signals in a wide frequency range is produced by devices located in a single borehole geophysical device that can move along the wellbore during well testing with a given equipment operation algorithm [2] . The equipment operation algorithm consists of the sequential execution of operations for recording the micro-seismic background in the well, acoustic impact on the mountain massif and re-registration of the signals of emission immediately after the impact. The energy of the induced acoustic emission, which is released in the process of acoustic impact, was determined by the method of calculating the spectral energy density similarly to the method of core research using the Intel Graph 2 program. Figure 1 shows an example of an acoustic emission signal in the BS10 reservoir of a terrigenous type of the Tevlin-Russkinskoye field in Western Siberia during an acoustic impact. On average, the increase in the SAE signal after the acoustic impact in the reservoir saturated with oil is several tens of percent relative to the background level. At the same time, individual acoustic emission events in the form of single actions of emission sources occur randomly and have characteristic parameters of the signal pulses of a finite duration of a certain form. At the same time, the presence of discrete frequencies with a certain maximum energy value is clearly manifested. Spectral analysis of the acoustic emission signal with a duration of one second (Figure 2) is performed using a window Fourier transform, in a sliding window of a given length. The transformation parameters were selected in an optimal way to identify the details of the process in time, which are characteristic of this type of reservoir. The spectrogram shows the distinct discrete nature of the set of frequencies in the signal caused by acoustic emission. At the same time, there is an increase in the energy of acoustic emission at discrete frequencies by several times as compared with the micro-seismic background. In this spectrum, you can select some discrete frequencies, interconnected nonlinear relationship (Table 1). In a fractured carbonate reservoir, after acoustic impact,

Table 1. Discrete frequencies from the spectrum of the acoustic emission signal at a depth of 2840 meters.

Figure 1. Acoustic emission signal and its spectrogram after acoustic impact in an oil well at a depth of 2840 meters in a terrigenous-type reservoir. U is the amplitude of acoustic responce, axes X is time in sec.

Figure 2. The spectrum of the acoustic emission signal at the point of 2840 m of the terrigenous collector before (AT, blue) and after (AT, bown) the acoustic impact. S is the spectrum amplitude.

the AE signal also has a pronounced discrete character. Moreover, after the acoustic impact, there was a significant increase in the energy of the induced acoustic emission compared with the background and the formation of pronounced discrete frequencies of energy in the spectrum of the emission signal. The nonlinearity of the radiation process caused by acoustic emission is manifested in the ratio of the main frequencies of the spectrum. But at the same time, there is a significant difference in the composition of discrete frequencies of the terrigenous and carbonate type spectra. The fractured carbonate type of collector has a much smaller number of discrete frequencies and their location in a narrower band as compared with the emission spectrum of a terrigenous collector. This is due to the different structure of the pore space of these reservoirs and the nature of changes in the dynamic processes of the stress state in them and the filtration of fluids.

In the column ω (N) it is the frequency from the spectrum in Figure 2, in the column ω (∑) is the sum of the frequencies of the various vibration modes, and in the column ∆ω is the difference of the frequency values measured and calculated.

The main result of this part of our research consists on that fact that the resonant frequencies of the acoustic response and their sums are in a good agreement with the measured ones: so ω₁ + ω₂ = ω₃: (1259 + 2648) = 3907 and the measured frequency ω₃ is 4048, then the difference between the measured and calculated value of ω₃ is −3.5%, that value is less than the error of experimental data and that coincides with the result of connection two oscillations in a nonlinear medium [1] .

3. Sound Diffraction Simulation on a Two-Dimensional Porous Oil-Saturated Reservoir of a Hierarchical Type, Located in an N-Layer Elastic Medium after Acoustic Excitation

Here we shall construct an algorithm for modeling the resonant acoustic response process of a two-phase porous hierarchical medium according to the results written in Table 1. We shall assume that the number N in the table will coincide with the number of the rank of hierarchical inclusions in the elastic reservoir. These inclusions will be of different features: the may have elastic features, plastic features and can be oil saturated. They can be a lаrge inclusion, that consists of some number of fewer inclusions. And these features are linked with the sequence of 1 up to 9 of N in Table 1.

In the paper [5] algorithm for modeling sound diffraction on 2D elastic hierarchical inclusion located in J-th layer of the N-layered medium is written using such designations and functions.

$G_{Sp,j}\left(M,M^0\right)$.source function of seismic field, for which the boundary problem had been formulated in the paper [6] ; $k_{1ji}^2={\omega }^2\left({\sigma }_{ji}/{\lambda }_{ji}\right)$.wave number for the longitudinal wave; the index ji denotes the identity of the properties of the medium within the heterogeneity, ja—out of heterogeneity; λ—constant Lame; σ—density of the medium; ω—circular frequency; $u=grad\phi$.displacement vector; φ⁰—potential of a normal seismic field in a layered medium in the absence of heterogeneity: ${\phi }_{ji}^0={\phi }_{ja}^0$. $l=1,\cdots ,L$.rank hierarchy of inclusions. Let l = 1. ${\lambda }_{ji}={\lambda }_{ja}+{\omega }_1{\lambda }_{ji\left(l=1\right)}$.6] . ${\phi }_{\left(l=1\right)-1}^0$.orresponds to the exciting field at the reservoir level. The system of equations [5] will be rewritten as (1):

$\begin{array}{l}\frac{\left(k_{1ji\left(l=1\right)}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=1\right)}{\iint }{\phi }_{\left(l=1\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+\frac{{\sigma }_{ja}}{{\sigma }_{ji}}{\phi }_{\left(l=1\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=1\right)}\right)}{{\sigma }_{ji\left(l=1\right)}2\text{π}}{\displaystyle \underset{C\left(l=1\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=1\right)}}{\partial n}\text{d}c={\phi }_{\left(l=1\right)}\left(M^0\right),M^0\in S_{C\left(l=1\right)},\\ \frac{{\sigma }_{ji\left(l=1\right)}\left(k_{1ji\left(l=1\right)}^2-k_{1j}^2\right)}{\sigma (M^0)2\text{π}}{\displaystyle \underset{Sc\left(l=1\right)}{\iint }{\phi }_{\left(l=1\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=1\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=1\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=1\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=1\right)}}{\partial n}\text{d}c={\phi }_{\left(l=1\right)}\left(M^0\right),M^0\notin S_{C\left(l=1\right)}\end{array}$.(1)

l = l + 1. that is: l = 2, inclusions $S{c}_{\left(l=2\right)}$.nd ${С}_{\left(l=2\right)}$.re located inside sections $S{c}_{\left(l=1\right)}$.nd contours ${С}_{\left(l=1\right)}$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=1\right)}+{\omega }_2{\lambda }_{ji\left(l=2\right)}$. the inclusions are plastic, ω₂ is taken from Table 1. ${\phi }_{\left(l=2\right)-1}^0={\phi }_{\left(l=1\right)}$. The system of Equations (1) is rewritten as:

$\begin{array}{l}\frac{\left(k_{1ji\left(l=2\right)}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=2\right)}{\iint }{\phi }_{\left(l=2\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=2\right)}}{\phi }_{\left(l=2\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=2\right)}\right)}{{\sigma }_{ji\left(l=2\right)}2\text{π}}{\displaystyle \underset{C\left(l=2\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=2\right)}}{\partial n}\text{d}c={\phi }_{\left(l=2\right)}\left(M^0\right),M^0\in S_{C\left(l=2\right)},\\ \frac{{\sigma }_{ji\left(l=2\right)}\left(k_{1ji\left(l=2\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=2\right)}{\iint }{\phi }_{\left(l=2\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=2\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=2\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=2\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=2\right)}}{\partial n}\text{d}c={\phi }_{\left(l=2\right)}\left(M^0\right),M^0\notin S_{C\left(l=2\right)}.\end{array}$.(2)

l = l + 1. that is: l = 3, inclusions $S{c}_{\left(l=3\right)m}$.nd ${С}_{\left(l=3\right)m}$.re located inside $S{c}_{\left(l=2\right)}$.nd ${С}_{\left(l=2\right)}$. $m=\text{1},\cdots ,M$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=2\right)}+{\omega }_3{\lambda }_{ji\left(l=3\right)}$. the inclusions are plastic also, ω₃ is taken from Table 1. ${\phi }_{\left(l=3\right)-1}^0={\phi }_{\left(l=2\right)}$.

The system of Equations (1) will be rewritten as [7] :

$\begin{array}{l}{\displaystyle \underset{1}{\overset{M}{\sum }}\{\frac{\left(k_{1ji\left(l=3\right)m}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=3\right)m}{\iint }{\phi }_{\left(l=3\right)m}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M}+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=3\right)m}}{\phi }_{\left(l=3\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=3\right)m}\right)}{{\sigma }_{ji\left(l=3\right)m}2\text{π}}{\displaystyle \underset{C\left(l=3\right)m}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=3\right)m}}{\partial n}\text{d}c={\phi }_{\left(l=3\right)m}\left(M^0\right)+{\left(\alpha p_2\right)}_m\},M^0\in S_{C\left(l=3\right)m},\\ \frac{{\sigma }_{ji\left(l=3\right)}\left(k_{1ji\left(l=3\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=3\right)}{\iint }{\phi }_{\left(l=3\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=3\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=3\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=3\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=3\right)}}{\partial n}\text{d}c={\phi }_{\left(l=3\right)}\left(M^0\right),M^0\notin S_{C\left(l=3\right)m}m=1,\cdots ,M.\end{array}$.(3)

${\left(\alpha p_2\right)}_m={\left(\left(1-\chi -\frac{K}{K_0}\right)p_2\right)}_m$

K—compression module, $\chi$.porosity, $K_0$.true module of phase compressibility, p₂—pore hydrostatic pressure for the m-th fluid-saturated inclusion and ${\omega }_3={\omega }_1+{\omega }_2$. (these frequency values are taken from Table 1). As we can see there are some inclusions which contain the fluid. So for each inclusion the features of $\left(\alpha p_2\right)$.or each m are different and they are located inside the inclusions sections $S{c}_{\left(l=\text{2}\right)}$.nd contours ${С}_{\left(l=\text{2}\right)}$.

l = l + 1, that is l = 4, inclusions $S{c}_{\left(l=4\right)}$.nd ${С}_{\left(l=4\right)}$.re located inside sections $S{c}_{\left(l=2\right)}$.nd contours ${С}_{\left(l=2\right)}$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=2\right)}+{\omega }_4{\lambda }_{ji\left(l=4\right)}$. ${\omega }_4$.s taken from Table 1. ${\phi }_{\left(l=4\right)-1}^0={\phi }_{\left(l=2\right)}$.The system of equations [5] will be rewritten as:

$\begin{array}{l}\frac{\left(k_{1ji\left(l=4\right)}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=4\right)}{\iint }{\phi }_{\left(l=4\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=4\right)}}{\phi }_{\left(l=4\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=4\right)}\right)}{{\sigma }_{ji\left(l=4\right)}2\text{π}}{\displaystyle \underset{C\left(l=4\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=4\right)}}{\partial n}\text{d}c={\phi }_{\left(l=4\right)}\left(M^0\right),M^0\in S_{C\left(l=4\right)},\\ \frac{{\sigma }_{ji\left(l=4\right)}\left(k_{1ji\left(l=4\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=4\right)}{\iint }{\phi }_{\left(l=4\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=4\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=4\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=4\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=4\right)}}{\partial n}\text{d}c={\phi }_{\left(l=4\right)}\left(M^0\right),M^0\notin S_{C\left(l=4\right)}.\end{array}$ (4)

l = l + 1, that is l = 5, inclusions $S{c}_{\left(l=5\right)mm}$.nd ${С}_{\left(l=5\right)mm}$.re located inside the sections $S{c}_{\left(l=\text{2}\right)}$.nd contours ${С}_{\left(l=\text{2}\right)}$. $mm=1,\cdots ,MM$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=4\right)}+{\omega }_5{\lambda }_{ji\left(l=5\right)}$. ${\phi }_{\left(l=5\right)-1}^0={\phi }_{\left(l=2\right)}$. The system of equations [5] will be rewritten as [7] :

$\begin{array}{l}{\displaystyle \underset{1}{\overset{MM}{\sum }}\{\frac{\left(k_{1ji\left(l=5\right)mm}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=5\right)mm}{\iint }{\phi }_{\left(l=5\right)mm}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M}+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=5\right)mm}}{\phi }_{\left(l=5\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=5\right)mm}\right)}{{\sigma }_{ji\left(l=5\right)mm}2\text{π}}{\displaystyle \underset{C\left(l=5\right)mm}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=5\right)mm}}{\partial n}\text{d}c={\phi }_{\left(l=5\right)mm}\left(M^0\right)+{\left(\alpha p_2\right)}_{mm}\},M^0\in S_{C\left(l=5\right)mm},\\ \frac{{\sigma }_{ji\left(l=5\right)}\left(k_{1ji\left(l=5\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=5\right)}{\iint }{\phi }_{\left(l=5\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=5\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=5\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=5\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=5\right)}}{\partial n}\text{d}c={\phi }_{\left(l=5\right)}\left(M^0\right),M^0\notin S_{C\left(l=5\right)mm}mm=1,\cdots ,MM.\end{array}$.(5)

${\left(\alpha p_2\right)}_{mm}={\left(\left(1-\chi -\frac{K}{K_0}\right)p_2\right)}_{mm}$. ${\omega }_5={\omega }_2+{\omega }_3$. (these frequency values are taken from Table 1).

l = l + 1, that is l = 6, inclusions $S{c}_{\left(l=\text{6}\right)mn}$.nd ${С}_{\left(l=\text{6}\right)mn}$.re located inside the sections $S{c}_{\left(l=\text{2}\right)}$.nd contours ${С}_{\left(l=\text{2}\right)}$. $mn=\text{1},\cdots ,MN$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=2\right)}+{\omega }_6{\lambda }_{ji\left(l=6\right)}$. ${\phi }_{\left(l=6\right)-1}^0={\phi }_{\left(l=2\right)}$.The system of equations [5] will be rewritten as [7] :

$\begin{array}{l}{\displaystyle \underset{1}{\overset{MN}{\sum }}\{\frac{\left(k_{1ji\left(l=6\right)mn}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=6\right)mn}{\iint }{\phi }_{\left(l=6\right)mn}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M}+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=6\right)mn}}{\phi }_{\left(l=6\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=6\right)mn}\right)}{{\sigma }_{ji\left(l=6\right)mn}2\text{π}}{\displaystyle \underset{C\left(l=6\right)mn}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=6\right)mn}}{\partial n}\text{d}c={\phi }_{\left(l=6\right)mn}\left(M^0\right)+{\left(\alpha p_2\right)}_{mn}\},M^0\in S_{C\left(l=6\right)mn},\\ \frac{{\sigma }_{ji\left(l=6\right)}\left(k_{1ji\left(l=6\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=6\right)}{\iint }{\phi }_{\left(l=6\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=6\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=6\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=6\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=6\right)}}{\partial n}\text{d}c={\phi }_{\left(l=6\right)}\left(M^0\right),M^0\notin S_{C\left(l=6\right)mn}mn=1,\cdots ,MN.\end{array}$.(6)

${\left(\alpha p_2\right)}_{mn}={\left(\left(1-\chi -\frac{K}{K_0}\right)p_2\right)}_{mn}$. ${\omega }_6={\omega }_2+{\omega }_4$. (these frequency values are taken from Table 1).

l = l + 1. that is l = 7, inclusions $S{c}_{\left(l=7\right)nm}$. ${С}_{\left(l=7\right)nm}$.re located inside the section $S{c}_{\left(l=1\right)}$.nd contour ${С}_{\left(l=1\right)}$. $nm=\text{1},\cdots ,NM$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=1\right)}+{\omega }_7{\lambda }_{ji\left(l=7\right)}$. ${\phi }_{\left(l=7\right)-1}^0={\phi }_{\left(l=1\right)}$. The system of Equations (1) will be rewritten as [8] :

$\begin{array}{l}{\displaystyle \underset{1}{\overset{NM}{\sum }}\{\frac{\left(k_{1ji\left(l=7\right)nm}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=7\right)nm}{\iint }{\phi }_{\left(l=7\right)nm}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M}+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=7\right)nm}}{\phi }_{\left(l=7\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=7\right)nm}\right)}{{\sigma }_{ji\left(l=7\right)nm}2\text{π}}{\displaystyle \underset{C\left(l=7\right)nm}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=7\right)nm}}{\partial n}\text{d}c={\phi }_{\left(l=7\right)nm}\left(M^0\right)+{\left(\alpha p_2\right)}_{nm}\},M^0\in S_{C\left(l=7\right)nm},\\ \frac{{\sigma }_{ji\left(l=7\right)}\left(k_{1ji\left(l=7\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=7\right)}{\iint }{\phi }_{\left(l=7\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=7\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=7\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=7\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=7\right)}}{\partial n}\text{d}c={\phi }_{\left(l=7\right)}\left(M^0\right),M^0\notin S_{C\left(l=7\right)nm}nm=1,\cdots ,NM.\end{array}$.(7)

${\left(\alpha p_2\right)}_{nm}={\left(\left(1-\chi -\frac{K}{K_0}\right)p_2\right)}_{nm}$. ${\omega }_7={\omega }_1+{\omega }_6$. (these frequency values are taken from Table 1).

l = l + 1, that is l = 8, inclusions with sections $S{c}_{\left(l=8\right)m}$. and contours ${С}_{\left(l=8\right)m}$. are located inside the section $S{c}_{\left(l=4\right)}$.nd contour ${С}_{\left(l=4\right)}$. $m^{\prime }=1,\cdots ,M^{\prime }$.orrespondingly, ${\lambda }_{ji}={\lambda }_{ji\left(l=4\right)}+{\omega }_8{\lambda }_{ji\left(l=8\right)m^{\prime }}$. ${\phi }_{\left(l=8\right)-1}^0={\phi }_{\left(l=4\right)}$. The system of equations [5] will be rewritten as [7] :

$\begin{array}{l}{\displaystyle \underset{1}{\overset{M^{\prime }}{\sum }}\{\frac{\left(k_{1ji\left(l=8\right)m^{\prime }}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=8\right)m^{\prime }}{\iint }{\phi }_{\left(l=8\right)m^{\prime }}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M}+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=8\right)m^{\prime }}}{\phi }_{\left(l=8\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=8\right)m^{\prime }}\right)}{{\sigma }_{ji\left(l=8\right)m^{\prime }}2\text{π}}{\displaystyle \underset{C\left(l=8\right)m^{\prime }}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=8\right)m^{\prime }}}{\partial n}\text{d}c={\phi }_{\left(l=8\right)m^{\prime }}\left(M^0\right)+{\left(\alpha p_2\right)}_m\},M^0\in S_{C\left(l=8\right)m^{\prime }},\\ \frac{{\sigma }_{ji\left(l=8\right)}\left(k_{1ji\left(l=8\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=8\right)}{\iint }{\phi }_{\left(l=8\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=8\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=8\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=8\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=8\right)}}{\partial n}\text{d}c={\phi }_{\left(l=8\right)}\left(M^0\right),M^0\notin S_{C\left(l=8\right)m^{\prime }}{m}^{\prime }=1,\cdots ,M^{\prime }.\end{array}$ (8)

${\left(\alpha p_2\right)}_{m^{\prime }}={\left(\left(1-\chi -\frac{K}{K_0}\right)p_2\right)}_{m^{\prime }}$. ${\omega }_8={\omega }_4+{\omega }_3$. (these frequency values are taken from Table 1).

l = l + 1, that is l = 9, inclusions with sections $S{c}_{\left(l=9\right)mmm}$.nd contours ${С}_{\left(l=9\right)mmm}$.re located inside the sections $S{c}_{\left(l=3\right)}$.nd $S{c}_{\left(l=5\right)}$.nd contours ${С}_{\left(l=3\right)}$. and ${С}_{\left(l=5\right)}$ $m=\text{1},\cdots ,M$.nd $mm=\text{1},\cdots ,MM$.orrespondingly,

$\begin{array}{l}{\displaystyle \underset{1}{\overset{MMM}{\sum }}\{\frac{\left(k_{1ji\left(l=9\right)mm}^2-k_{1j}^2\right)}{2\text{π}}{\displaystyle \underset{Sc\left(l=9\right)mmm}{\iint }{\phi }_{\left(l=9\right)mmm}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M}+\frac{{\sigma }_{ja}}{{\sigma }_{ji\left(l=9\right)mmm}}{\phi }_{\left(l=9\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=9\right)mmm}\right)}{{\sigma }_{ji\left(l=9\right)mmm}2\text{π}}{\displaystyle \underset{C\left(l=9\right)mmm}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=9\right)mmm}}{\partial n}\text{d}c={\phi }_{\left(l=9\right)mmm}\left(M^0\right)+{\left(\alpha p_2\right)}_{mmm}\},M^0\in S_{C\left(l=5\right)mmm},\\ \frac{{\sigma }_{ji\left(l=9\right)}\left(k_{1ji\left(l=9\right)}^2-k_{1j}^2\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{Sc\left(l=9\right)}{\iint }{\phi }_{\left(l=9\right)}}\left(M\right)G_{Sp,j}\left(M,M^0\right)\text{d}{\tau }_M+{\phi }_{\left(l=9\right)-1}^0\left(M^0\right)\\ -\frac{\left({\sigma }_{ja}-{\sigma }_{ji\left(l=9\right)}\right)}{\sigma \left(M^0\right)2\text{π}}{\displaystyle \underset{C\left(l=9\right)}{\oint }{G}_{Sp,j}}\frac{\partial {\phi }_{\left(l=9\right)}}{\partial n}\text{d}c={\phi }_{\left(l=9\right)}\left(M^0\right),M^0\notin S_{C\left(l=9\right)mmm}mmm=1,\cdots ,MMM.\end{array}$.(9)

${\left(\alpha p_2\right)}_{mmm}={\left(\left(1-\chi -\frac{K}{K_0}\right)p_2\right)}_{mmm}$. ${\omega }_9={\omega }_5+{\omega }_3$. (these frequency values are taken from Table 1).

4. Conclusions

Thus, acoustic emission in porous media of reservoir rocks has the ability of permanently acting radiation. The main feature of acoustic emission is a reflection of stress relaxation and high sensitivity to external influences. In the case of fluid-filled pore space, the medium generates an induced energy of elastic energy that is adequately associated with its properties and state, and the radiation is non-linear.

Using data from Table 1, a variant of the system of structural adjustment algorithms for a nonlinear two-phase porous medium after its dynamic excitation is constructed. The values of the resonant frequencies are taken from the table. The rank of fluid saturated inclusions located in a porous medium with plastic properties is correlated with the number of the resonant frequencies. Such a system of algorithms was developed for the first time. That system can vary by the features (elastic, plastic and saturated) of the inclusions. By that analysis, we can take into account how many saturated inclusions can move inside the reservoir and can increase the potential oil recovery. After the first acoustic impact if it is needed, we can repeat the second impact, formulate the analogue of Table 1 and construct the system of algorithms linked with data of a new table.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Rabinovich, M. and Trubetskov, D. (2000) Introduction to the Theory of Oscillations and Waves of SIC “Regular and Chaotic Dynamics” Moscow.
[2]	Dryagin, V.V. (2018) The Use of Induced Acoustic Emission from Reservoirs for the Detection and Extraction of Hydrocarbons. Georesources, 20, 246-260. https://doi.org/10.18599/grs.2018.3.246-260
[3]	Chebotareva, I.Y., Volodin, I.A. and Dryagin, V.V. (2016) Generation of a Low-Frequency Branch of Acoustic Emission in Affected Rocks. Doklady Earth Sciences, 468, 205-208. https://doi.org/10.1134/S1028334X16050160
[4]	Chebotareva, I.Y., Volodin, I.A. and Dryagin, V.V. (2017) Acoustic Effects in the Deformation of Structually Inhomogeneous Media. Acoustical Physics, 63, 86-94. https://doi.org/10.1134/S106377101606004X
[5]	Hachay, O.A. and Khachay, A.Y. (2013) Modeling of Seismic and Electromagnetic Fields in Hierarchically Heterogeneous Environments. In: Bulletin of the South Ural State University. Series “Computational Mathematics and Computer Science”, 2, 48-56. https://doi.org/10.14529/cmse130204
[6]	Hachay, O. and Khachay, A. (2017) Acoustic Wave Monitoring of Fluid Dynamic in the Rock Massif with Anomaly Density, Stressed and Plastic Hierarchic Inclusion. Computational and Experimental Studies of Acoustic Waves, Intec Open, 2017, Chapter 4, pp. 63-80. https://doi.org/10.5772/intechopen.70590
[7]	Hachay, O. and Khachay, A. (2014) Reflection of the Processes of Non-Equilibrium Two-Phase Filtration in Oil-Saturated Hierarchical Media in the Data of Active Wave Geophysical Monitoring. GIAB, 4, 232-238.
[8]	Hachay, O. and Khachay, A. (2016) Modeling the Propagation of a Seismic Field in a Layered-Block Elastic Medium with Hierarchical Plastic Inclusions. GIAB, 12, 318-326.