^{1}

^{*}

^{1}

The new method for prediction of earthquake center zone is suggested. The method is based on feature of amount of registered information to reach its maximum upon some condition regulating interrelation of major parameters of used distributed measuring system. The mathematical basis of suggested is based on known integrated Shannon formula of amount of information and integral limitation condition, expressing fixed position of used sensors. As a result of held research es, new method of information trangulation method for determination of earthquake center zone is suggested. The mathematical grounding and the operational algorithm of the method are given.

Development of methods for prediction of earthquake attracted attention of many scientists [

According to [

In line with [

As it is noted in [

At the same time, numerous sums of effects of rocks cracking and random character of exposures of seismical activity lead to statistical nature of forerunners of earthquakes. Particularly, as regard a normal electromagnetic irradiation, it is practically not possible to select the most informative frequencies, which should be observed. The optimal azimuth angle of receiving of such radiation usually is also unknown, which creates additional problems for optimal solution of earthquake center prediction problem. Due to the above, it can be concluded, that the system destined for gathering and processing of seismological information of electromagnetic radiation type should guarantee receipt and processing maximal amount of such information. This condition leads us to construct the single informational model of electromagnetic radiation before earthquakes, optimization of which makes it possible to predict the center of earthquakes using informational criteria. Thus the aim of this research is development of new information.

First of all, we should stress out that the intensity of electromagnetic radiation depends on two factors:

1) Distance between the receiver and center zone of earthquake—L.

2) Frequency of electromagnetic radiation—F.

Above dependences in both case expressed by fading of intensity of electromagnetic radiation.

Taking into account above, we propose, that each seismosensor can be characterized in the three-sized space as it is shown in

Taking into account of above, we can find the signal/noise ratio ψ of the received signal as

ψ = ( L , F ) = ψ 0 + ∂ ψ ∂ L ⋅ d L + ∂ ψ ∂ F ⋅ d F , (1)

where ψ 0 = ψ ( L = L 0 ; F = F 0 ) .

Now considering the whole net of seismosensors consisting of n ones, and taking into account, that optimal system of processing is designed in such a manner, that it receives signals of duration T i , i = 1 , n ¯ , where value of T i increases depending on i. In this case, the amount of information, which is gathered from all n sensors can be assessed as follows

Φ = ∑ i = 1 n T i Δ T log 2 [ ψ 0 + ∂ ψ ∂ L d L + ∂ ψ ∂ F d F + 1 ] . (2)

Using formula (2), we can solve two type of optimization task, which allow us to form two methods for prediction of earthquakes center using static distributed set of receivers and maximal informational criteria.

1) Receive of electromagnetic radiation at the same frequency.

2) Receive of electromagnetic radiation at the same distance used for confirmation of found center of earthquake on first method.

In the first case, the ratio of signal/noise ψ in the channel of propagation of Low-frequency electromagnetic radiation from the source of these signals as far as receiver in the first approximation can be found as

ψ = ψ 0 + ψ ′ L , (3)

where: ψ 0 = ψ ( L = 0 ) ; ψ ′ = d ψ d L .

During the whole period of information retrieval, total amount of gained information can be found as

M 0 = ∑ i = 1 n M i = ∑ i = 1 n T i Δ t log 2 ( ψ 0 + ω ′ L i + 1 ) , (4)

where T i —duration of information retrieval from the sensor numbered as i.

Next, we use a limitation condition, formed as a result of fixed construction of the set of sensors, i.e.

∑ i = 1 n ( ψ 0 + ψ ′ L i + 1 ) = c o n s t .

Taking into account above, the functional of effectiveness can be formed as

Φ = ∫ 0 T max [ T Δ t log 2 ( ψ 0 + ψ ′ L + 1 ) + λ ( ψ 0 + ψ ′ L + 1 ) ] d T ,

where λ is L’Agrange multiplier.

According to the principle of optimal lowering of dimensionality [

Solution of the formulated maximization task using Euler’s formula gives us following type of said function

L = ψ 0 + 1 | ψ ′ | − T ψ 0 T max | ψ ′ | . (5)

As a result, we have the possibility to carry out an adaptive control of seismosensors, i.e. all system of gathering and processing of seismic information.

On the basis of above result, we can propose a new method of informational for casting of epicenter of expected earthquakes. We assume that seismosensors are placed on the territory with high seismic risk (

In order to forecast the epicenter of earthquake we should designate set of values { L i } , and set of values { T i } , where dependence between L i and T i should

be in line with formula (3), which guarantee maximum value of total information gathered from seismosensors. Selected three seismosensors S 1 , S 2 , S 3 will be commutated during time period T i determined by formula (3). Under above conditions computer should detect maximum amount of information, if the point A is actually the center of expected earthquake. Selected set of sensors should be moved along the territory in order to detect maximum of gathered information. More strictly no sensors, but selector contour of three sensors should be moved over the territory via fixed set of sensors and movement of this contour should be stopped if maximal amount of information is reached. Then using known set of L i , we can find needed placecenter of earthquake.

Receive of electromagnetic radiation at the same distance used for confirmation of found center of the earthquake on first method.

In this case, the ratio signal/noise ψ in the channel of propagation of low-frequency electromagnetic radiation from the source of these signals as far as receiver in the first approximation can be found as

ψ = ψ 0 + ψ ′ F ⋅ F , (6)

where ψ ′ F = ∂ ψ ∂ F ; F—frequency of electromagnetic radiation.

During the whole period of information retrieval total amount of gathered information can be found as

M 0 = ∑ i = 1 n M i = ∑ i = 1 n T i Δ t log 2 ( ψ 01 + ψ ′ F ⋅ F + 1 ) , (7)

where T i —duration of information retrieval from the sensor numbered as i.

Then we use a limitation condition

∑ i = 1 n ( ψ 01 + ψ ′ F ⋅ F i + 1 ) = c o n s t

or

∑ i = 1 n F i = c o n s t . (8)

A limitation condition (6) means, that band of reasonably received frequencies is limited.

Taking into account of (5) and (6) we can form the functional of effectiveness as follows:

Φ = ∫ 0 T max [ T Δ T log 2 ( ψ 01 + ψ ′ F ⋅ F + 1 ) + λ ( ψ 01 + ψ ′ F ⋅ F + 1 ) ] d T (9)

where

T—time of information retrieval; λ —multiplier of L’Agrang.

According to the principle of optimal lowering of dimensionality [

Solution of above optimization task using Euler’s formula gives us following type of said function

F = φ ( T ) = ψ 0 + 1 | ψ ′ F | − T ⋅ ψ 0 T max | ψ ′ F | (10)

As a result, we obtain the possibility to carry out an adaptive control of seismosensors, i.e. whole system of gathering and processing of seismic information. This does mean, that sensor with the lowest frequency of received signal should be examined during uppermost time period in order to reach maximal efficiency of the system.

On the basis of above result, we can propose a new second method for informational forecasting of center of expected earthquakes. We assume that seismosensors are placed on the territory with high seismic risc, forming a rectangular Net (

In order to forecast the place of earthquake we should designate a set of values { F i } ; i = 1 , n ¯ , and set of values { T i } ; i = 1 , n ¯ , where dependence between F i and T i should be in line with formula (8), which guarantees reaching of maximum value of total information, gathered from seismosensors. Here we should note, that each sensor consists of n receiver with fixed frequency { F i } ; i = 1 , n ¯ .

Hence, each sensor is compound on n number of receiver and makes it possible to receive in the frequency band { F 1 , F 2 , ⋯ , F n } .

Selected contour of four seismosensors (in

It should be noted that the above second method is used for confirmation of place of earthquake center, determined by the first method.

On the basis of two above methods we can formulate single general algorithm for prediction of earthquake’s center zone (

1) Using of fixed rectangular net of seismosensors for determination of center zone on first method.

2) Construction of non-rectangular mobile micro-net consisting of three or four sensors located at the same distance from found center zone.

3). Application of the second method.

4) Confirmation, that two results on points 1 and 3 are the same.

Taking into consideration of above we conclude, that a proposed new informational method of forecasting of earthquake center zone may increase effectiveness of measures taken for global forecasting of seismic events.

Thus, the new method for prediction of earthquake center zone is suggested. The method is based on feature of amount of registered information to reach its maximum upon some condition regulating interrelation of major parameters of used distributed measuring system. The mathematical basis of suggested method is based on known integrated Shannon formula of amount of information and integral limitation condition, expression fixed position of used sensors.

Now we can note major result of held research.

The new method of information triangulation method for determination of earthquake center zone is suggested.

The mathematical grounding and the operational algorithm of the method are given.

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

Asadov, H.H. and Nuriyeva, L.I. (2020) Development of Information Triangulation Method for Prediction of Earthquake Center Zone Using Distributed Measurements of Emitted Electromagnetic Radiation. Positioning, 11, 1-9. https://doi.org/10.4236/pos.2020.111001