ALGORITHMS FOR CONTROLLING THE ONSET OF
ANTHROPOGENIC PROCESSES IN URBAN DEVELOPMENT AREAS IN SEISMICALLY ACTIVE
REGIONS
Telman Aliev, Naila Musaeva
Institute of Control Systems, Azerbaijan National Academy of Sciences Azerbaijan University of Architecture and Construction telmancyber@gmail.com , musanaila@gmail.com
Abstract
It is noted that during changes in the geodetic state of urban development area, as well as the technical condition of the objects located there, the noisy signals coming from the corresponding sensors contain noise correlated with the useful signal. It is shown that the characteristics of the relationship between the useful signal and the noise are informative attributes of the beginning of anomalous natural phenomena, as well as changes in the stress-strain state of objects of the urban development complex. We have developed algorithms for calculating the estimates of the relay cross-correlation functions, normalized cross-correlation functions, as well as the correlation coefficient between the useful signal and the noise. It is noted that the use of the algorithms for calculating these characteristics in control systems allows detecting nascent anomalous processes, as well as defects of building structures at the initial latent stage and to prevent the possibility of accident situations.
Keywords: control, seismically active regions, noise, noisy signal, cross-correlation function, correlation coefficient
I. Introduction
It is known that the main objective of urban development is the design and development of cities, districts, complexes, residential and industrial zones, etc. Therefore, architectural and urban development objects are the territories where cities, villages, administrative districts, etc. are located. And buildings, structures, residential and industrial complexes, etc. are considered to be architectural and urban development objects.
At the same time, urban development areas are subjected to various kinds of natural and man-made impacts. For instance, man-made impacts are a result of the activities of industrial and military facilities, the movement of rail, road, sea, river transport, communications network, etc. Depending on the duration, cyclidty and degree of anthropogenic impacts are determined by the impact on civil and military complexes, residential buildings, national economy facilities, etc. Other dangerous effects on urban development areas are such natural processes as earthquakes, landslides, hurricanes, floods, etc. And one of the most dangerous natural phenomena are landslides, which lead to displacement of rocks as a result of seismic shocks, overwetting, washout of the slope and other processes. Such phenomena can occur on rocky terrain as well as on sand, clay, etc. Powerful landslides can cause an earthquake and vice versa, landslides can occur as a result of an earthquake.
In order to avoid the negative impact of natural and anthropogenic impacts on the state of objects, urban development monitoring is carried out. For this purpose, geodetic monitoring is
first carried out, i.e., soil processes are studied, subsidence, displacement and subsidence of massifs are evaluated, groundwater is analyzed, possible seismically hazardous and landslide zones and areas are identified. At the same time at the stage of construction as well as during operation such defects of building structures and constructions as misalignments, cracks, shifts, displacements are identified and their stress-strained state is assessed. Thus, the current state and dynamics of changes in the technical condition of buildings and structures are assessed, possible accidents are predicted and appropriate recommendations are given.
In order to carry out such monitoring, appropriate sensors, such as displacement sensors, pressure sensors, reinforcement and anchor loads, crack aperture sensors, inclinometers, accelerometers, etc., are installed in the information locations of the analyzed area and building structures. The data on the geodetic state of the urban development complex as well as of the building structures and constructions are transmitted in the form of signals to the central monitoring station. The obtained digital information is recorded, processed, analyzed and conclusions are made about the geodetic condition of the urban development area and the stressstrain state of the building structures, the presence of defects and the degree of damage. In addition, conclusions are made about the causes that led to this or that defect.
However, existing systems of urban monitoring do not provide control of the latent period of change of the geodetic state of the territory under study, as well as the technical condition of the objects located there [1-2]. This is especially important for territories affected by seismic vibrations, landslides, dangerous vibrations and other natural processes, because after each weak but frequent earthquake or a small but repeated landslide, invisible and undetectable microscopic cracks, deformations, bends, etc., appear which subsequently can lead to serious damage that requires large costs [1-2].
In [1, 3-4] it is shown that changes in the geodetic state, as well as formation of even the smallest damages of objects of the urban planning complex is accompanied by the appearance of additive noise, which are correlated with the useful signals of noisy signals coming from the above sensors. The conducted studies have shown that the use of the technology for processing and analyzing noise and characteristics of the interrelation of noisy signals in the systems of noise monitoring of the urban development complex allows identifying the early latent period of defect and damage initiation, determining the dynamics of its development, and reducing the risks of accidents of urban development objects in seismically active and landslide-prone regions.
II. Problem statement
It is known that in practice, real signals are the sum of useful signals X(t) and noise E(t), i.e., G(t) = X(t) + E(t). Because of the contamination of the useful signals X(t) with the noise E(t) there are tangible errors in determining the estimates of their correlation functions Rxx(p). In this case, the total noise E(t) is the sum of noise Ei(t) from the influence of external factors and the noise E2(t), correlated with the useful signal, which arises when the geodetic state of the investigated area changes and a defect is generated during the operation of objects, i.e., E(t) = E1(t)+E2(t).
Suppose that G(t) is a sampled stationary random signal with normal distribution, consisting of the useful signal X(t) and the noise E(t) with mathematical expectation mE=0. In this case the formula for calculating the estimate of the correlation function Rcc (p) can be represented as [1, 34]:
Rgg(p) = ^ZiLi G(iAt)G((i + ju)At) = Rxx(m) + Rxe(m) + Rex(m) + Ree(j) , (1) where RXX(j) = ~Zf=i X(iAt)x((i + j)At) is the autocorrelation function of the useful signal
X(t); RXE([i) iX(iAt)E((i + p)At), REX(fi) = ±lf= 1E(iAt) x((i + p)At) are the cross-
correlation functions between X(t) and E(t); REE(p) = DE = iE(iAt)E((i + pi)At) is the variance of the noise ( ).
Given that R EE (p) =0 at the total error will be
f 2 Rxe( 0 ) + R EE( 0 ) wh e n p = 0 g(p) = {R xeO-i) + Rex(p) when p±0 . (2)
Because of this, the inequality Rxx (p) ^ RGG (p) takes place. It is clear from this that the basic information needed for monitoring is contained in the estimates of the noise characteristics ( ), Rex(jp)> Ree(0) . Since the noise E(t) cannot be isolated from the noisy signal G(t), it is impossible to draw adequate conclusions about the geodetic state of construction objects on the estimate ( ) in urban development areas it is not possible to make adequate conclusions about the geodetic state and the stress-strain state of construction objects. In this regard, there is an obvious need to create algorithms and technologies to find the estimates of the noise variance the and cross-correlation functions RX E(p), REX(p) between the useful signal and the noise.
III. Algorithms for calculating the characteristics of the relationship between the useful signal and the noise
As shown above, in the normal geodetic state of the urban development area and objects located there, noise ( ) ( ) arises from random external factors that have no correlation with the useful signal ( ). However, at the beginning of the latent period of change in the geodetic state of the area and the technical condition of objects, the noise ( ) emerges which is correlated with the useful signal. Therefore, from this point on, the correlation between the useful signal ( ) and the total noise ( ) ( ) ( ) is non-zero. In this case, the initiation and development of malfunctions is essentially reflected in the estimates of the cross-correlation functions and the correlation coefficient between X(t) and E(t) [1, 3-4]. Therefore, for controlling the latent period of the origin and the dynamics of the development of natural anomalies and technical malfunctions, it is quite often more appropriate to use in monitoring systems the estimates ( ), ( ), as well as the correlation coefficient between the useful signal and the noise as a carrier of diagnostic information. In the following paragraphs, we show the possibility of using estimates of the relay cross-correlation functions, normalized cross-correlation functions, as well as the correlation coefficient between the useful signal and the noise as informative attributes.
It is known that the estimates of the relay cross-correlation functions can be calculated from the expression [1]:
rxe(p) = i sgnX(iAt) E((i + p)At), (3)
(+1 w hen X(iAt) > 0 where sgnX(iAt) = { 0 when X(iAt) = 0 .
( - 1 w he n X( iA t) < 0
Obviously, to use this formula it is necessary to determine the samples of the noise E( iA t) and the useful signal X(iAt) that cannot be measured directly or extracted from the noisy signal G(t) [1]. Let us consider one of possible variants of the approximate estimation of the relay cross-correlation function ( ) between the useful signal ( ) and noise ( ) as a result of calculation of the relay correlation function ( ) of the noisy signal ( ) [1]:
RrGg(p) = 111 1 sgnG(iAt) G((i + p)At), (4)
(+1 w he n G( iA t) > 0 where sgn G(iAt) = { 0 when G(iAt) = 0 .
( - 1 w he n G( iA t) < 0
It is shown in [1, 3-4] that the estimate of the relay cross-correlation function RrX*E(jAt) at different time shifts i between X(iAt) and E(iAt) can be determined from the formula:
Rtxe(j) = 1Zi=1 sgnG(iAt) ( g((i + j)At) - 2 g((i + j + 1)At) + g((i + j + 2 )At)). (5)
It is clear that with the normal geodetic state of the area and the technical condition of the objects due to the lack of correlation between X (t) and E (t) the estimate of the relay cross-correlation function ( ) between the useful signal and the noise will be close to zero. When various anomalies and defects occur, the value of the estimate of the relay cross-correlation function will vary depending on the degree of correlation between ( ) ( ). And the distinctive feature of this algorithm is that if there is a correlation between ( ) ( ), differences in the values of the estimates RxE(j = 0), Rx*E(j = At), RxE(j = 2At), ... at different time instants unambiguously reflect the dynamics of malfunction development, which makes it possible to obtain reliable information about the risk of accidents.
In addition, to assess the latent period of the initiation of anomalous seismic processes and landslides, as well as the technical condition of objects in the urban development area, it is advisable to use estimates of the normalized cross-correlation function, as well as the correlation coefficient between the useful signal ( ) and the noise ( ).
In this case, using formula (5), the calculation of the normalized cross-correlation function can be reduced to the form [1]:
pUJ) = (6)
^71 E
where the mean square deviation o^ of the noise can be calculated from the expressions [1, 3-4]:
JRcc(0 ) — 2RcC(At) + Rcg(2A t) for the general case
G g — a __. (7)
WR C C( 0 ) — R C C(A t) fo r th e s p e ci al c as e
It is known that the value of the normalized cross-correlation function at ^=0 is the correlation coefficient. Therefore, the value of the correlation coefficient between the useful signal X(t) and the noise E(t) can be calculated from the expression:
rx*E — pU 0) — (8)
¡1 ..
^71 E
IV. Technologies for determining the latent period of changes in the geodetic state and technical condition of objects of the urban development area
To assess the geodetic state and technical condition of objects in the urban development area, the estimates of the relay and normalized cross-correlation functions, as well as the correlation coefficients between the useful signal X(t) and the noise E(t) for each of the monitored parameters should be calculated at different time instants . The obtained values of the
estimates are entered into the database of informative attributes, which can be represented in the form of matrices:
51 =
52 =
Rxl El(ß) t 1 Rrx[ El(ß) t2 ■ ■ ■ Rxi El(ß)tk Rx2 E2(ß) 11 Rxl E2(ß)
RXnEn(ß)t 1
)t2
Rrx*nEnUt)t2
Prxl El(ß) t1 PX\ El(ß)t2 Pxl E2(ß) 11 Px*2 E2(ß)t2
XlEl\ r* /
x2 e2(
r x:*№tk
XnEn
( )
(9)
PX\ El(ß) tk Px2 E2(P)tk
.PxnEn(ß)t1 PrxnEnU*)t2 ■ ■ ■ Px"nEnUi)tk.
(10)
53 = r vr* ' xlEltl ' X2E2tl yr* 'xlElt2 ■ yr* <x2B2t2 ■ yr* lxiEnk yr* f*2*2tk
yr* - XnEnt1 yr* ' xnEnt2 ' XnEntk
(11)
After appropriate training, the geodetic state and technical condition of the objects in the urban development area is identified. For instance, if the geodetic indicators of the noise characteristics have changed in some limited area, and the estimates of the noise characteristics of the technical condition have changed only for nearby objects, it can be regarded as the beginning of a hidden period of landslide occurrence. If the geodetic indices of noise characteristics have changed for quite an extensive territory, and the estimates of noise characteristics of technical conditions have changed for all objects simultaneously on this territory, this indicates the beginning of a latent period of preparation for an earthquake. It is possible to judge about the dynamics of geodetic processes by the intensity of changes in the noise characteristics. And the noise-characteristics will take different values depending on the type of terrain, e.g., rock, sand, clay, etc. At the same time it is necessary to carry out monitoring of technical condition of buildings, constructions, structures and other objects in the urban development already at the initial stage of occurrence of anomalous natural phenomena.
Depending on the values of the relay and normalized cross-correlation functions, as well as the correlation coefficients between the useful signal X(t) and noise E(t) at different time instants, conclusions are made about the geodetic state of the urban development area, for instance: 1 -ground waters are leaking; 2 - displacements and shifts in the soil are observed; 3 - subsidence of soil masses are observed; 4 - a source of vibrations is detected; 5 - a source of seismic vibrations is detected; 6 - a landslide zone is detected, etc. For each of these states the intensity of the process is defined: a - not developing; b - developing very slowly; c - developing slowly; d - developing with average intensity; e - developing intensively; f - developing very intensively; g - developing rapidly, etc. Then the matrices of the geodetic state of the urban development area take the form:
G51 =
rxi El(ß)t 1 rxi eM) t2 ■ ■ ■ Rxn El GO
XXlEl
DT*
*X2B2
R ZMt 1
RX*2 E2(ß)t2
)t k
R xn^(ß)tk
XlEl\ (
1X*nEnUt)t 1
Rl
M)
t2
R
x-ne-,
( )
1 — a
2 — b
(12)
GS 2 =
PXI El(ß) 11 Px*2 eM) 11
Pxl El(ß) t2
Px2 E2(p) t2
.PxnEn(ß)t 1 PX*nEn(ß)t2
Pxl El(ß)tk
Px2 E2(P)tk pX*nEn(ß)tk.
1 -a
2 — b
(13)
GS3 = ryr* 'xlEltl yr* 'x2B2tl yr* xiEit2 ' yr* <X2E2t2 ■ yr* " tk yr* " tk
y-r* - XnEnt1 yr* 'XnEnt2 ' y-r* xnEntk
1 -a
2 — b
Similar conclusions are made about the presence or absence of latent damage and the dynamics of their development for structures and objects of the urban development area, for instance: 0 - technical condition is serviceable; 1 - operational without damage; 2 - limited operational with minor damage; 3 - limited operational with minor damage developing intensively; 4 - non- operational; 5 - pre-emergency; 6 - emergency, etc., depending on the values of the noise characteristic estimates at a given time instant. Then the matrices of the technical condition of structures and objects take the form:
751 =
E,(P)t 1 Rx\E,(p) t
RxUM t2
*X2E2
RX* E, (P)tk
R ZMt 1
?r*
R ry*^(M)tk
R];
M)t
»
t2
R
X-nEi
(p)
tk
ts2 =
PXl E,(P)t 1 Px\ E,(P)t2 Px2 E2(P)t 1 Pll E2(P)t2
p£ E„(p) 11 Prx*,E„(P)t2
Px1 E1(P) tk
Px2 E2Ut)tk p£ EvU^tk-
(15)
(16)
753 = rrr* 'xlEltl yV* 'x2E2tl yT* xiEit2 ' yV* <X2E2t2 ■ yV* " tk yV* " '*2*2tfc
. Xn En 11 yV* 'XnEnt2 ' xnEntk
(17)
IV. Discussion
The conducted analysis shows that using the matrices of noise characteristics (9)-(11), (12)-(14), (15)-(17), it is possible to identify an early stage of occurrence of changes in the geodetic state of an urban development area, as well as the technical condition of the objects located there. At the same time the objective of geodetic monitoring is to identify the latent period of the emergence of hazardous natural phenomena. The main objective of technical monitoring is to detect latent defects in the roof, facade, microcracks in the walls, the possibility of deformations, etc. This can be done by calculating the characteristics of the relationship between the useful signal and the noise. The combination of geodetic and technical control makes it possible to objectively assess the condition of the urban development complex; to prevent in a timely manner the risk of the risk of destruction of buildings, structures and constructions in seismically active regions; to identify the causes and sources of risks both for individual buildings and for the entire complex; to reduce the cost of maintenance and major repairs; to alert personnel to the possibility of destruction of objects with catastrophic consequences.
References
[1] Aliev, T. A. Noise Control of the Beginning and Development Dynamics of Accidents, Springer, 2019, 201 p. DOI: 10.1007/978-3-030-12512-7
[2] Jooste, A.F. & de Kock, I.H. & Musango, J.K. (2019). A systematic literature review of sustainable urban planning challenges associated with developing countries. South African Journal of Industrial Engineering. Vol. 30(3), Special Edition, pp. 253-261 253. DOI:10.7166/30-3-2247
[3] Aliev, T.A. & Musaeva, N.F. & Suleymanova, M.T. (2018). Algorithms for Indicating the Beginning of Accidents Based on the Estimate of the Density Distribution Function of the Noise of Technological Parameters. Automatic Control and Computer Science. Vol. 52. Issue 3. pp. 231-242. DOI: 10.3103/S0146411618030021
[4] Aliev, T.A. & Musaeva, N.F. (2019). Technologies for Early Monitoring of Technical Objects Using the Estimates of Noise Distribution Density. Journal of Automation and Information Sciences. Vol. 51, no. 9, pp. 12-23.