Development of Information Triangulation Method for Prediction of Earthquake Center Zone Using Distributed Measurements of Emitted Electromagnetic Radiation

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 researches, 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.

In line with [8], speed of occurence of seismological signals before earthquakes depending on volume of cracked rock at depth is a primary local source of precursor signals. As it is noted in [9], researchers on possible earthquake prediction must be carried out on a deterministic basis. That implies the necessity to carry out combined ground and space observations of earthquake precursor. Earthquake precursory phenomes include not only emission of ULF electromagnetic waves rad, but also visible and near infrared light. According to [10], during seismic event, amplitude of electromagnetic radiation increases as long as the crack continues to grow because new atomic bonds are severed.
As it is noted in [11], seismo-electromagnetic precursory based prediction method appeared in the scientific literature for many decades but motivation is given recently in [12] [13]. Well-known triangulation method based on propagation of electromagnetic waves was considered in [14] [15]. The principle of this method is based on speed of movement of electromagnetic waves ant fixation of their arrival. Graphical realization of the method based on registration of arrived wave at three stations is illustrated in Figure 1 [11].
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.

Mathematical Basics
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 Figure 1, where the point A characterize any single seismosensor, which located at the distance A L from center of earthquake, receives electromagnetic radiation of frequency A F , and duration of receipt is A T (Figure 2).  Taking into account of above, we can find the signal/noise ratio ψ of the received signal as where ( ) 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 , 1, , where value of i T increases depending on i. In this case, the amount of information, which is gathered from all n sensors can be assessed as follows 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 where: During the whole period of information retrieval, total amount of gained information can be found as ( ) where i T -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.
( ) Taking into account above, the functional of effectiveness can be formed as where λ is L'Agrange multiplier.
According to the principle of optimal lowering of dimensionality [2], we should find such type of optimal function ( ) L T ϕ = which would lead the functional of effectiveness to its maximal value.
Solution of the formulated maximization task using Euler's formula gives us following type of said function 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.

Development of New Method
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 (Figure 3).  be in line with formula (3), which guarantee maximum value of total information gathered from seismosensors. Selected three seismosensors 1 2 3 , , S S S will be commutated during time period i T 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 i L , 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 lowfrequency electromagnetic radiation from the source of these signals as far as receiver in the first approximation can be found as where F F ψ ψ ∂ ′ = ∂ ; F-frequency of electromagnetic radiation.
During the whole period of information retrieval total amount of gathered information can be found as ( ) where i T -duration of information retrieval from the sensor numbered as i.
Then we use a limitation condition (8)

Computer
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: where T-time of information retrieval; λ -multiplier of L'Agrang.
According to the principle of optimal lowering of dimensionality [16], we should find such type of optimal function ( ) which would lead the functional of effectiveness (7) to its maximal value. Solution of above optimization task using Euler's formula gives us following type of said function 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 ( Figure 5), which consists of following steps: Positioning  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.

Conclusion
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 Yes No end 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.