Next Article in Journal
Nondestructive Evaluation of Concrete Bridge Decks with Automated Acoustic Scanning System and Ground Penetrating Radar
Next Article in Special Issue
Square-Root Unscented Information Filter and Its Application in SINS/DVL Integrated Navigation
Previous Article in Journal
Polystyrene Oxygen Optodes Doped with Ir(III) and Pd(II) meso-Tetrakis(pentafluorophenyl)porphyrin Using an LED-Based High-Sensitivity Phosphorimeter
Previous Article in Special Issue
Data Fusion Architectures for Orthogonal Redundant Inertial Measurement Units
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fire Source Range Localization Based on the Dynamic Optimization Method for Large-Space Buildings

1
Luoyang Institute of Science and Technology, Luoyang 471000, China
2
College of Electrical Engineering, Henan University of Technology, Zhengzhou 450001, China
3
Bureau of Geophysical Prospecting, China National Petroleum Corporation, Zhuozhou 072751, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(6), 1954; https://doi.org/10.3390/s18061954
Submission received: 27 April 2018 / Revised: 9 June 2018 / Accepted: 11 June 2018 / Published: 15 June 2018
(This article belongs to the Collection Multi-Sensor Information Fusion)

Abstract

:
This paper is concerned to the fire localization problem for large-space buildings. Two kinds of circular fire source arrangement localization methods are proposed on the basis of the dynamic optimization technology. In the Range-Point-Range frame, a dynamic optimization localization is proposed to globally estimate the circle center of the circular arrangement to be determined based on all the point estimates of the fire source. In the Range-Range-Range frame, a dynamic optimization localization method is developed by solving a non-convex optimization problem. In this way, the circle center and the radius are obtained simultaneously. Additionally, the dynamic angle bisector method is evaluated. Finally, a simulation with three simulation scenes is provided to illustrate the effectiveness and availability of the proposed methods.

1. Introduction

In recent decades, more and more large-space buildings have been utilized as product storage and manufacture places, which are usually designed with complicated structures and are crowded with many kinds of materials and products. For these places, fire safety has become one of the most important and difficult problems, and increasing attention has been by governments, factories, firemen, researchers, and engineers [1,2,3,4,5,6]. For fighting against fire disasters, a critical important issue is how to determine the location of a fire source, in order to provide accurate information and necessary guidance for firefighting. Therefore, it is necessary to develop fire location methods to determine the fire source location as quickly and accurately as possible, for intelligent firefighting, especially in large-space buildings.
In fact, the fire source location method is one kind of wireless indoor positioning technique [7,8,9,10,11]. The existing fire source location methods are mainly developed on the basis of image processing technology, fiber-optic temperature measurement, temperature sensor arrays, and smoke sensor arrays. The fire location systems based on image information are commonly used for the large open spaces, and are easy impacted by barriers [12,13,14]. The fiber-based temperature measurement methods utilize the impacts of temperature on the anti-Stokes spectral lines in the Raman scattering process in an optical fiber to determine the fire sources position, which is suitable for tunnels and other such scenes, however, the cost of installation and maintenance of optical fiber is too expensive [15,16,17,18]. Moreover, several temperature field-based positioning methods aiming at environments, such as mines and forests, are given in [19,20,21,22]. These methods usually require that the fire releases more heat. Recently, the fire monitoring and location methods based on temperature/smoke sensor arrays have been paid increasing attention, in which the measurements are sometimes non-line-of-sight propagated. The fire source burning process produces the hot gases or smoke, which rise up to the ceiling and spread in a nearly circular shape [23,24,25,26,27,28]. With two temperature sensor arrays, a prototype system for determining the fire location was proposed in [18], in which three temperature sensors of a four-sensor array were utilized to estimate the angle of the fire source to the sensor array, and the fire source location can be obtained with the help of two estimated angles and the distance between two sensor arrays. This method was identified as the far-field mode in [24], due to the assumption that the distance between the fire source and the sensor array is much large than the distance intervals in the sensor array. In [24], a far-field fire location approach and a near-field method were presented and compared. In the near-field fire location method, the distance and the angle of the fire source to the sensor array were estimated simultaneously, with signal source localization of the planar wave fronts, but not the planar wave fronts, which was utilized in the far-field approach. The position estimation of the smolder source was realized in [28], using the smoke sensor arrays arranged in the planar circular mode. In [25], a fire location position algorithm was also developed with the signal source localization of a bilinear sensor array. As shown in [27,29], the output results of the above fire source location methods were a number of dispersedly estimated point coordinates of the fire source. The fire source point coordinate estimates were solved locally and independently with the sensor data sampled in the current time, but not the global sensor data sampled in the time series. Therefore, the clustering algorithm was introduced in the fire source localization method in the Range-Point-Range (RPR) frame in [27]. Additionally, two effective fire range estimation approaches were proposed in the Range-Range-Range (RRR) frame in [29]. The circle range estimates were obtained globally in the form of the inscribed circle and the circum-circle of a quadrangle, which was solved recursively using a dynamic angular bisector method. However, it did not discuss the optimality of the circle ranges solved by the dynamic angular bisector method in [29].
Inspired by the above discussion, in this paper, the method of determining the circle range which covers the fire source is developed based on the dynamic optimization theory. Firstly, in the RPR frame, all the point estimates of the fire source obtained by the fire source point location method are considering to obtain a circle center of the circular fire source range to be determined, the sum distance from which each fire point estimate is minimizes. The radius is the standard deviation of the distances from the circle center to every fire point estimate. Actually, compared with to the clustering algorithm in [27], this method determines the circle center and radius based on all the fire point estimates. Therefore, the circular fire source range is determined globally, according to [29]. However, the dynamic optimization problem solved in this algorithm is not recursive. When several new fire point estimates are obtained by using the current measurements of sensor arrays, all the obtained fire point estimates in the time series should be considered to determine the circle center and radius of the circular fire source range to be determined. Although the location accuracy can be improved along with the fire source range determining process, the computing complexity is increasing poor. In order to reduce the computing complexity of the dynamic optimization problem, a circular fire source range location method is also developed in the RRR frame. In [29], the quadrangle is solved by using the dynamic angle bisector method with global information, which, in this paper, is utilized to determine a new circum-circle, the circle center and radius of which is solved by minimizing the distances of the circle center to the four vertices of the given quadrangle. An interesting discovery is that the circle center is similar to the one obtained by the dynamic angle bisector method. It seems that the circle center obtained by using the dynamic angle bisector method is optimal in the sense of minimizing the distances of the circle center to the four vertices of the quadrangle solved with the global information.
The major contributions of this paper are three-fold. First, with all the fire point estimates, a circular fire source range estimation method is proposed in the RPR frame. As shown in [27], the clustering algorithms used to obtain the circular fire source range inevitably discard some fire point estimates. However, in this paper, the circular fire source range is solved globally by using a dynamic optimization algorithm based on all the point estimates of the fire source. Secondly, the above circular fire source range location method is developed in the RRR frame to reduce the computational complexity. The global angle information in time series are recursively utilized to obtain a quadrangle by using the dynamic angle bisector method. What is different with the method in [29] is that the circle center and the radius are determined by solving a non-convex optimization problem. In fact, the non-convex optimal problem can be treated as a standard to evaluate the dynamic angle bisector method. Additionally, the optimality of the dynamic angle bisector method is analyzed by the numeric examples in the simulation section.
It should be noted that, for the case where the angle information are not uniform on the two sides of the bisectors, the proposed circular fire source arrangement methods have better location results than the dynamic angle bisector method in [29], as illustrated in the simulation.
The rest of this paper is organized as follows: In Section 2, the fire location problem is presented, and a fire source point localization principle is briefly introduced. On this basis, two fire source range localization methods based on the dynamic optimization algorithm are proposed in Section 3. The first one is designed with the help of the fire point estimates in the RPR frame in Section 3.1, and in Section 3.2, the second method is developed with the quadrangle solved by the dynamic angle bisector method, which is studied in the RRR frame. In Section 4, a simulation with three different simulation scenes are provided to prove the effectiveness of the proposed methods. Finally, the conclusion of this paper is given in Section 5.

2. Problem Formulation

The structured fire scene studied in this paper is similar to the ones in [6,12,29]. An ignition fire source is considered in a large-space building with static wind. The fire source burning process produces the hot gases, which rise up to the ceiling and spread in a nearly circular shape at a constant current velocity from, and around, the center and form a hot upper layer [6,12,29], as shown in Figure 1. In Figure 1, two temperature sensor arrays are placed in different locations on the ceiling to monitor the air temperature. In each sensor array, four homogeneous temperature sensors are arranged in a square shape with the side length d . The distance between the two sensor arrays is L . The distances from the fire source to the two sensor arrays are D 1 and D 2 , respectively. Then the expectation of the temperature field around the fire source can be described as [12]:
E { T ( x , y , t ) } | ( x x 0 ) 2 + ( y y 0 ) 2 = c 2 = f ( t )
where ( x 0 , y 0 ) is the point coordinate of fire source and f ( t ) is the temperature function at location ( x , y ) , the distance from which to the fire source is c .
For the structured fire scene described above, several effective fire source point location methods are proposed in [6]. Taking the famous far-field algorithm as an example, the main steps of the fire source point location scheme can be sketched as follows:
(1)
Delay Estimation. As shown in Figure 1, there are different distances from the fire source to different temperature sensors. According to the temperature field expectation in Equation (1), the delay time of the same temperature time series spread to different sensors can be estimated by the correlation function method [12], or the gray relation analysis method [29].
(2)
Angle Estimation. Denote τ i , j , ( i , j = 1 , 2 , 3 , 4 ) as the delay time from S i to S j , f A as the sample rate of the temperature sensors, α ( k ) is an intersection angle crossed by the horizontal line and the line from the sensor array to the fire source, as shown in Figure 1. Based on the planar waves assumption of the far-field algorithm, for the sensor set { S 1 , S 2 , S 3 } , one can obtain:
| v | τ 12 f A = d cos ( α ( k ) )
| v | τ 13 f A = d sin ( α ( k ) )
then:
α ^ ( k ) = arctan ( τ 12 τ 13 ) , v ^ = f A d cos ( α ^ ( k ) ) τ 12
Similarly, for different sensor sets in the same sensor array, such as { S 1 , S 2 , S 4 } and { S 1 , S 3 , S 4 } , different estimates of α ( k ) can also be obtained. In the sensor array, the sensor at S 1 is the reference node. For the sake of convenience, the angle estimate is denoted as α ^ i , j ( k ) , in which k indicates the sampled times, i { A , B } represents the different sensor arrays, and j = 1 , 2 , 3 signifies the different sensor sets.
(3)
Fire Source Point Estimation. For every combination { α ^ A , j ( k ) , α ^ B , l ( k ) } , ( j , l = 1 , 2 , 3 ) , the fire source point can be estimated as follows [12]:
x ^ 0 , j , l ( k ) = L tan ( α ^ B , l ( k ) ) tan ( α ^ A , j ( k ) ) + tan ( α ^ B , l ( k ) ) + x 1
y ^ 0 , j , l ( k ) = L tan ( α ^ B , l ( k ) ) tan ( α ^ A , j ( k ) ) + tan ( α ^ B , l ( k ) ) + y 1
where ( x 1 , y 1 ) is the coordinate of S 1 in the sensor array A.
According to the fire source point estimation method given in Equations (5) and (6), nine estimates of the fire source can be solved at one time. Afterward, one can obtain a great deal of discretely ruleless estimate points. This kind of result limits the ability to guide firefighters. Therefore, some modified fire location methods are developed based on the fusion method [12,29]: the clustering algorithms [27]. In the next section, two kinds of location methods are designed to determine more accurate and compact ranges of the fire source, based on the above fire source point estimation method.

3. Main Results

3.1. Dynamic Optimization Localization Method in the RPR Frame

According to the fire source point estimation method mentioned in the above section, several ruleless discrete location points can be solved at one time. In order to fuse this location information into one estimate point, the mean of the three angle estimates for each sensor array is taken as the final angle estimate at one time, then the fusion fire source points are solved by Equations (5) and (6) in [12]. As noted in [29], this fusion result is solved only based on the sensor data at the current time, but not all of the sensor data in the time series. Therefore, in [27], a fire source range location method in the RPR frame is proposed based on an improved clustering algorithm. The fire source point estimates in the time series are utilized to determine a circular range which contains the real fire source with a large probability. However, this fire source circle range does not contain all the fire source point estimates in the clustering algorithm [29].
In this section, a fire source range location method is studied in the RPR frame, based on the dynamic optimization algorithm. All the estimates ( x ^ 0 , j , l ( t ) , y ^ 0 , j , l ( t ) ) , j , l = 1 , 2 , 3 ; t = 1 , 2 , , k obtained by Equations (5) and (6) contain some information of the fire source. Thus, all of them should be considered to deal with the globally determined circular range. Therefore, in this section, the circle center of the circular sharp arrangement to be determined is considered as the solution of the following optimization problem:
min s u m ( ( x ^ 0 k x 0 ) 2 + ( y ^ 0 k y 0 ) 2 )
where ( x ^ 0 k , y ^ 0 k ) is the matrix including all the fire source point estimates in the time series obtained by the fire source point estimation method mentioned in the above section.
It is indicated in Equation (7) that the sum of the distances between the circle center to be solved and every fire source point estimates in the time series should be the minimum.
Denote the solution of the optimization problem of Equation (1) as ( x ^ 0 ( k ) , y ^ 0 ( k ) ) . Take it as the circle center of the circular range to be determined, and take the standard deviation of the distances from ( x ^ 0 ( k ) , y ^ 0 ( k ) ) to every fire source point estimate ( x ^ 0 , j , l ( t ) , y ^ 0 , j , l ( t ) ) , j , l = 1 , 2 , 3 ; t = 1 , 2 , , k as the radius of the circular range to be determined. Namely:
r ( k ) = s t d { ( x ^ 0 , j , l ( t ) x ^ 0 ( k ) ) 2 + ( y ^ 0 , j , l ( t ) y ^ 0 ( k ) ) 2 } , j , l = 1 , 2 , 3 ; t = 1 , 2 , , k
where, ( x ^ 0 , j , l ( t ) , y ^ 0 , j , l ( t ) ) , j , l = 1 , 2 , 3 ; t = 1 , 2 , , k are just the constituent elements of the matrix ( x ^ 0 k , y ^ 0 k ) .
In this way, a circular fire source range is determined by solving the optimization problem in Equation (7).
Remark 1.
All the fire source point estimates in the time series are utilized to solve the dynamic optimization problem in Equation (7). Therefore, the circle center is solved based on the global information.
However, the dynamic optimization problem solved in this algorithm is not recursive. When several new fire point estimates are obtained by using the current measurements of sensor arrays, all the obtained fire point estimates in the time series should be considered to determine the circle center and radius of the circular fire source range to be determined. Although the location accuracy can be improved along with the fire source range determining process, the computing complexity is increasing poor. In order to reduce the computing complexity of the dynamic optimization problem in the RPR frame, a circular fire source range location method is also developed in the RRR frame in the next subsection.

3.2. Dynamic Optimization Localization Method in the RRR Frame

In the fire source localization scene shown in Figure 1, for each sensor array i , i { A , B } , three angle estimates can be obtained at each time k , according to Equations (2) and (3). In the RRR frame [29], the up and down bounds of the angle estimate range can be obtained by using the statistical mean and variance of the three angle estimates:
α ^ i ( k ) = 1 3 j = 1 3 α ^ i , j ( k ) , δ i 2 = 1 2 j = 1 3 ( α ^ i , j ( k ) α ^ i ( k ) ) 2
α ^ i H ( k ) = α ^ i ( k ) + δ i , α ^ i L ( k ) = α ^ i ( k ) δ i
For the global angle estimates in the time series, the corresponding statistical mean and variance can be recursively calculated by:
α ^ i k = α ^ i ( k ) + ( k 1 ) α ^ i k 1 k , ( δ i k ) 2 = ( α ^ i ( k ) α ^ i k ) 2 + ( k 1 ) ( δ i k 1 ) 2 k 1
Then, the up and down bounds of global angle estimate range can be given by:
α ^ i k , H ( k ) = α ^ i k ( k ) + δ i k , α ^ i k , L ( k ) = α ^ i k ( k ) δ i k
The overlapped area of the global angle estimate ranges of the two sensor arrays contains the fire source with large probability, as shown in Figure 2. In [29], a dynamic angle bisector method was proposed to determine a circular fire source arrange, the circle center and radius of which were calculated on the basis of this overlapped area, which is a quadrangle. Actually, using the statistical features of the global data to obtain the quadrangle which contains the fire source with large probability is an excellent data compression method, which can also be used to reduce the computing complexity of the dynamic optimization problem in the proposed fire source range location method in Section 3.1.
As show in Figure 2, set the reference vertices’ coordinates of two sensor arrays as ( x A , y A ) , ( x B , y B ) , then the four vertices C j ( c j k , d j k ) , ( j = 1 , 2 , 3 , 4 ) of the quadrangle overlapped by the two global angle estimate ranges can be solved by [29]:
{ c 1 k = tan ( α ^ A k , L ) x 1 tan ( α ^ B k , H ) x 2 y 1 + y 2 tan ( α ^ A k , L ) tan ( α ^ B k , H ) ; d 1 k = y 2 / tan ( α ^ B k , H ) y 1 / tan ( α ^ A k , L ) x 2 + x 1 1 / tan ( α ^ B k , H ) 1 / tan ( α ^ A k , L ) c 2 k = tan ( α ^ A k , H ) x 1 tan ( α ^ B k , H ) x 2 y 1 + y 2 tan ( α ^ A k , H ) tan ( α ^ B k , H ) ; d 2 k = y 2 / tan ( α ^ B k , H ) y 1 / tan ( α ^ A k , H ) x 2 + x 1 1 / tan ( α ^ B k , H ) 1 / tan ( α ^ A k , H ) c 3 k = tan ( α ^ A k , L ) x 1 tan ( α ^ B k , L ) x 2 y 1 + y 2 tan ( α ^ A k , L ) tan ( α ^ B k , L ) ; d 3 k = y 2 / tan ( α ^ B k , L ) y 1 / tan ( α ^ A k , L ) x 2 + x 1 1 / tan ( α ^ B k , L ) 1 / tan ( α ^ A k , L ) c 4 k = tan ( α ^ A k , H ) x 1 tan ( α ^ B k , L ) x 2 y 1 + y 2 tan ( α ^ A k , H ) tan ( α ^ B k , L ) ; d 4 k = y 2 / tan ( α ^ B k , L ) y 1 / tan ( α ^ A k , H ) x 2 + x 1 1 / tan ( α ^ B k , H ) 1 / tan ( α ^ A k , L )
On this basis, the fire localization problem can be formulated to the following dynamic optimization problem, in which the circle center and radius of the circular fire region to be determined can be solved simultaneously:
min r 2 s . t . ( x 0 c j k ) 2 + ( y 0 d j k ) 2 r 2 ( j = 1 , 2 , 3 , 4 )
where ( x 0 , y 0 ) and r are the center’s coordinate to be solved, respectively.
It should be noted that Equation (14) is only one of many optimization formulations and not the one and only. It is noted that Equation (14) is a non-convex optimization problem and it is difficult to solve. Thereby, in order to successfully solve this non-convex optimization, there is an available method to transform it to a set of convex optimization problems. In fact, the non-convex optimization problem given by Equation (14) can be expressed equivalently as the following set of four convex optimization problems, namely:
min ( x 0 c 1 k ) 2 + ( y 0 d 1 k ) 2 s . t . ( x 0 c j k ) 2 + ( y 0 d j k ) 2 ( x 0 c 1 k ) 2 + ( y 0 d 1 k ) 2 ( j = 2 , 3 , 4 )
min ( x 0 c 2 k ) 2 + ( y 0 d 2 k ) 2 s . t . ( x 0 c j k ) 2 + ( y 0 d j k ) 2 ( x 0 c 2 k ) 2 + ( y 0 d 2 k ) 2 ( j = 1 , 3 , 4 )
min ( x 0 c 3 k ) 2 + ( y 0 d 3 k ) 2 s . t . ( x 0 c j k ) 2 + ( y 0 d j k ) 2 ( x 0 c 3 k ) 2 + ( y 0 d 3 k ) 2 ( j = 1 , 2 , 4 )
min ( x 0 c 4 k ) 2 + ( y 0 d 4 k ) 2 s . t . ( x 0 c j k ) 2 + ( y 0 d j k ) 2 ( x 0 c 4 k ) 2 + ( y 0 d 4 k ) 2 ( j = 1 , 2 , 3 )
Denote the optimal solutions of Equations (15)–(18) as ( x 0 , j k , y 0 , j k ) , ( j = 1 , 2 , 3 , 4 ) , and set:
r 0 , j k = ( x 0 , j k c 1 k ) 2 + ( y 0 , j k d 1 k ) 2 , j = 1 , 2 , 3 , 4
Then, the optimal global solution on the radius of the circle to be determined is:
r 0 k = min { r 0 , 1 k , r 0 , 2 k , r 0 , 3 k , r 0 , 4 k } = r 0 , p k
where p { 1 , 2 , 3 , 4 } . The corresponding estimate of center’s coordinate of the expected circle is:
x ^ 0 k = x 0 , p k , y ^ 0 k = y 0 , p k
Therefore, in the RRR frame, a fire source localization arrangement covering the fire source is presented with the dynamic optimization theory and the data compression method based on the statistical features of the global data.
Remark 2.
Although the set of four convex optimization problems shown in Equations (15)–(18) can equivalently transform the non-convex optimization problem of Equation (14), but it cannot be directly solved with the CVX software package, which requires that constraints must be in the most simplified form. Thus, the optimization problems in (15)–(18) should be equivalently expressed as:
min ( x 0 c 1 k ) 2 + ( y 0 d 1 k ) 2 s . t . 2 x 0 ( c 1 k c j k ) + 2 y 0 ( d 1 k d j k ) 2 ( c 1 k ) 2 + ( d 1 k ) 2 ( c j k ) 2 ( d j k ) 2 ( j = 2 , 3 , 4 )
min ( x 0 c 2 k ) 2 + ( y 0 d 2 k ) 2 s . t . 2 x 0 ( c 2 k c j k ) + 2 y 0 ( d 2 k d j k ) 2 ( c 2 k ) 2 + ( d 2 k ) 2 ( c j k ) 2 ( d j k ) 2 ( j = 1 , 3 , 4 )
min ( x 0 c 3 k ) 2 + ( y 0 d 3 k ) 2 s . t . 2 x 0 ( c 3 k c j k ) + 2 y 0 ( d 3 k d j k ) 2 ( c 3 k ) 2 + ( d 3 k ) 2 ( c j k ) 2 + ( d j k ) 2 ( j = 1 , 2 , 4 )
min ( x 0 c 4 k ) 2 + ( y 0 d 4 k ) 2 s . t . 2 x 0 ( c 4 k c j k ) + 2 y 0 ( d 4 k d j k ) 2 ( c 4 k ) 2 + ( d 4 k ) 2 ( c j k ) 2 + ( d j k ) 2 ( j = 1 , 2 , 3 )
which can be directly solved with the CVX software package.
Remark 3.
Theoretically speaking, the angle bisector method in [29] subjectively fixes the circle center and the radius of the circular fire source arrangement to be determined, while the dynamic optimization method looks more reasonable because it can obtain a dynamic optimal solution of the estimate of circle center and the radius which, in the meantime, is also global. Clearly, there is a question that the subjective circle center cannot theoretically be guaranteed to be reasonable and perfect. The optimization problem in this section can be taken as a standard to evaluate the effectiveness of the angle bisector method. It should be noted that the circle center and the radius are estimated, respectively, although they have a very close connection. In fact, the circle center can be considered as a solution of the optimization problem that the minimum distance from a point to be determined in the quadrangle to cover all the four vertices, and the circum-circle radius estimated by the angle bisector method is the shortest one of the four distances from the vertices to the determined circle center, but not the distance solved by the last optimization problem. Therefore, there are some areas in the quadrangle that cannot be covered by the circum-circle estimated by the angle bisector method. For the dynamic optimal method in the RRR frame, the circle center and radius are estimated simultaneously in a consistent optimization standard. Thus, the circum-circle determined by the estimated circle center and radius cover the whole quadrangle. Therefore, the output of the dynamic optimization localization method in the RRR frame is another kind of circum-circle, which covers the whole quadrangle. It differs from the output of the angle bisector method with the circum-circle.
Remark 4.
Similar to the angle bisector method in [29], the dynamic optimal localization method in the RRR frame is sensitive to the statistical character of the measurement data on the fire position. For cases where the mean of the fire source point estimates are not zero, the location performances will be reduced to some extent. Nevertheless, the dynamic optimization localization method in the RPR frame has stronger robustness for the change of statistical characters of the measurement data.

4. Simulation

In this section, a numerical simulation with three different simulation scenes are provided to prove the effectiveness and availability of the two proposed dynamic optimal localization methods. The simulation settings are similar to those set in [12,29], as shown in Table 1 and Table 2.
No matter the localization methods designed in the RPR frame or the RRR frame, the processes of delay time estimation and angle estimation are coincident, even for the principle of fire source point estimation. Therefore, the simulation of angle estimates do not influence the validity of the simulation results. The angle estimate simulation method in [29] is adopted to simplify the simulation process. Namely, several fire source point estimates are simulated firstly, and then the angle estimates can be computed according to the triangle principle.
In order to show the effectiveness and robustness of the proposed fire localization methods, they are compared with three other kinds of fire location methods in three different fire source localization simulation scenes.
Simulation Scene 1. The fire source point estimates are simulated based on the real fire coordinate by adding Gaussian noise. The estimation results of Algorithm A and Algorithm B1 are compared in Figure 3.
Although Algorithm A and Algorithm B1 can both be implemented with global data to obtain the circular fire source arrangement estimation, the estimated fire source arrangements are different. This is because the standard to guide the algorithm is different, which has been analyzed with detail in Remark 3. In addition, in Algorithm B1, it is subjective to take the standard deviation of the distances from the circle center to every fire source point estimate as the radius of the circular range to be determined. Thus, the circular fire source arrangement determined by Algorithm B1 is greater than Algorithm A, and the circle centers of the real fire point and its estimates by Algorithms A and B1 are shown in Table 3.
Simulation Scene 2. The fire source point estimates are simulated based on the real fire coordinate by adding Gaussian noise. The estimation results of Algorithm A and Algorithm B2 are compared in Figure 4.
As illustrated in Figure 4, the circle centers of the circular fire source arrangement estimated by Algorithm A and Algorithm B2 are very similar, as shown in Table 4, the relative error is O ( 10 2 ) . Although the circles estimated by Algorithm A and Algorithm B2 are both the circum-circle of the quadrangle, the circum-circle obtained by Algorithm B2 can cover the whole quadrangle, while the one solved by Algorithm B2 cannot. The reason for this result has been analyzed in Remark 3.
Simulation Scene 3. The fire source point estimates are simulated based on the real fire coordinate by adding random mixed noise combined with Gaussian noise, Rayleigh noise, and sine wave signals. As shown in Remark 4, the dynamic optimal localization method in the RRR frame and the angle bisector method are sensitive to the statistical characters of the measurement data. Thus, the localization results are only compared by three robust localization methods: with the dynamic optimization localization method in the RPR frame, the localization method in the RPR frame with clustering technology, and the VB-ASCK method.
The simulation results given in Figure 5 are the average results of 100 Monte Carlo simulations. It is inevitable that the data from several measurements cannot be utilized to estimate the circle center and radius of the circular fire source arrangement. Therefore, the other two localization methods, with global measurement data, have better accuracy than Algorithm D. In Algorithm C, the mixed noise is taken as Gaussian noise with unknown covariance. Although the specific variational Bayesian technology is used to adaptively estimate the unknown covariance, the estimation accuracy of Algorithm C is worse than that of the dynamic optimization localization method in the RPR frame. As shown in Table 5, the mean estimation error of Algorithm D is 1.7, while that of Algorithm C is 1.2 and that of Algorithm B1 is 0.9.

5. Conclusions

In this paper, two circular fire source arrangement localization methods are proposed on the basis of the dynamic optimization technology. In the RPR frame, the dynamic optimization localization method is global, but not recursive, so the computing complexity of the dynamic optimization process is increasingly poor. Inspired by the data compression technology in the angle bisector method, a dynamic optimization localization method is developed by solving a non-convex optimization problem. The non-convex optimization problem can be treated as a standard to evaluate the dynamic angle bisector method. It should be noted that the circle radius is subjectively obtained by the dynamic optimization localization method. How to build an optimal index to solve the circle radius and the circle center optimally and synchronously is still an open topic for future research.

Author Contributions

G.W. and X.F. designed the fire source location method; G.W. and Z.Z. implemented the experiments and analyzed the experimental data; and X.F. and Z.Z. wrote and modified the paper.

Funding

This work was supported by the National Natural Science Foundation of China (grant nos. 61503174, 61333005, and 61304258), the Natural Science Foundation of Henan Province of China (grant nos.162102210196 and 15A413011), and the Scientific Research Foundation of Henan University of Technology (grant no. 2015RCJH14).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cram, D.; Hatch, C.; Tyler, S.; Ochoa, C. Use of Distributed Temperature Sensing Technology to Characterize Fire Behavior. Sensors 2016, 16, 1712. [Google Scholar] [CrossRef] [PubMed]
  2. Lu, Y.; Wang, H.; Wang, K. Fire Location for High and Large-Span Space Buildings based on Binocular Stereo Vision. Int. J. Smart Home 2016, 10, 271–282. [Google Scholar]
  3. Sun, Y.; Huang, H.; Li, F.; Yu, D.; Hao, C. Research on Zoning Control of Fire in Large Space Building. J. Wuhan Univ. Technol. 2017, 2, 144–147. [Google Scholar]
  4. Mao, Y.; Wang, H.; Lu, Y.; Qin, L. Precise Localisation Method for Fire Sources in Large-space Buildings. Comput. Appl. Softw. 2016, 2, 169–172. [Google Scholar]
  5. Ying, L.; Wang, H.; Qin, L. Accurate fire location method in high and large-span space buildings. J. Jilin Univ. 2016, 6, 2067–2073. [Google Scholar]
  6. Zhang, Q.; Yang, X. Numerical Simulation and Detection Response Analysis of Fire in a Large Space Building. In Proceedings of the International Conference on Intelligent Computation Technology and Automation, Changsha, China, 25–26 October 2014; pp. 388–391. [Google Scholar]
  7. Song, T.; Tang, B.; Zhao, M.; Deng, L. An accurate 3-D fire location method based on sub-pixel edge detection and non-parametric stereo matching. Measurement 2014, 4, 160–171. [Google Scholar] [CrossRef]
  8. Hackner, A.; Oberpriller, H.; Ohnesorge, A.; Hechtenberg, V.; Müller, G. Heterogeneous sensor arrays: Merging cameras and gas sensors into innovative fire detection systems. Sens. Actuators B Chem. 2016, 231, 497–505. [Google Scholar] [CrossRef]
  9. Song, H.; Lee, Y. A Development of Fire Detective Control System using the Fusion of a CCD Image and an Infrared Image. J. Korean Inst. Illum. Elec. Install. Eng. 2017, 5, 60–66. [Google Scholar] [CrossRef]
  10. Jiang, M.; Yao, B.; Cheng, Y. Response characteristics of distributed optical fiber temperature testing system in highway tunnel fire detection. Fire Saf. Sci. 2013, 4, 194–200. [Google Scholar]
  11. Feng, R. Optical Fiber Grating Temperature Sensing Tunnel Fire Alarm System. Trans. Sci. Technol. 2009, 1, 89–91. [Google Scholar]
  12. Sun, M.; Tang, Y.; Yang, S.; Dong, F. Fiber optic distributed temperature sensing for fire source localization. Meas. Sci. Technol. 2017, 8, 085102. [Google Scholar] [CrossRef]
  13. Sun, M.; Tang, Y.; Yang, S.; Li, J.; Sigrist, M.W.; Dong, F. Fire Source Localization Based on Distributed Temperature Sensing by a Dual-Line Optical Fiber System. Sensors 2016, 16, 829. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, X.; Li, K.; Li, J.; Hong, Y. Security technology research on high-temperature fire area blasting in surface mine. Opencast Min. Technol. 2010, 4, 73–75. [Google Scholar]
  15. Martínez-de Dios, J.R.; Merino, L.; Caballero, F.; Ollero, A. Automatic Forest-Fire Measuring Using Ground Stations and Unmanned Aerial Systems. Sensors 2011, 11, 6328–6353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Shan, Y.; Ma, Y.; Fu, H.; Li, W.; Wang, C. Application of distributed optical fiber temperature measurement system in coal mine fire monitoring system. Chin. J. Sens. Actuators 2014, 5, 22–25. [Google Scholar]
  17. Noordin, N.; Ney, H. Localization in wireless sensor network for forest fire detection. In Proceedings of the 3rd IEEE International Symposium on Telecommunication Technologies, Kuala Lumpur, Malaysia, 28–30 November 2016; pp. 87–90. [Google Scholar]
  18. Berentsen, M.; Kaiser, T. Fire Location Estimation using Temperature Sensor Arrays. In Proceedings of the 12th International Conference on Automatic Fire Detection, Gaithersburg, MD, USA, 25 March 2001; pp. 432–443. [Google Scholar]
  19. Wang, S.; Berentsen, M.; Kaiser, T. Signal processing algorithms for fire localization using temperature sensor arrays. Fire Saf. J. 2005, 8, 689–697. [Google Scholar] [CrossRef]
  20. Stern-Gottfried, J.; Rein, G.; Bisby, L.; Torero, J.L. Experimental review of the homogeneous temperature assumption in post-flashover compartment fires. Fire Saf. J. 2010, 4, 249–261. [Google Scholar] [CrossRef]
  21. Duan, S.; Cheng, K.; Ge, Q.; Wen, C. Wireless Intelligent Fire Fighting Systems Software Platform R&D. In Proceedings of the 2011 International Conference on Internet of Things and International Conference on Cyber, Physical and Social Computing, Dalian, China, 19–22 October 2011; pp. 94–99. [Google Scholar]
  22. Duan, S.; Ge, Q.; Wen, C. New method for fire source clustering localization based on sensor arrays. J. Nanjing Univ. Sci. Technol. 2011, 35, 59–65. [Google Scholar]
  23. Zhuang, Z.; Zhang, X.; Kalin, L. Method Research on Early Fire Source Localization Based on Planar Circular Sensor Arrays. Chin. J. Sens. Actuators 2009, 8, 1208–1212. [Google Scholar]
  24. Ge, Q.; Wen, C.; Duan, S. Fire Localization Based on Range-Range-Range Model for Limited Interior Space. IEEE Trans. Instrum. Meas. 2014, 9, 2223–2237. [Google Scholar] [CrossRef]
  25. Liu, H.; Darabi, H.; Banerjee, P.; Liu, J. Survey of Wireless Indoor Positioning Techniques and Systems. IEEE Trans. Syst. Man Cybern. Part C 2007, 37, 1067–1080. [Google Scholar] [CrossRef]
  26. He, S.; Hu, T.; Chan, S.H.G. Contour-based Trilateration for Indoor Fingerprinting Localization. In Proceedings of the 2015 ACM Conference on Embedded Networked Sensor Systems, Seoul, Korea, 1–4 November 2015; pp. 225–238. [Google Scholar]
  27. Sun, G.; Chen, J.; Guo, W.; Liu, K.R. Signal processing techniques in network-aided positioning: A survey of state-of-the-art positioning designs. IEEE Signal Process. Mag. 2005, 22, 12–23. [Google Scholar]
  28. Fonollosa, J.; Solórzano, A.; Marco, S. Chemical Sensor Systems and Associated Algorithms for Fire Detection: A Review. Sensors 2018, 18, 553. [Google Scholar] [CrossRef] [PubMed]
  29. Sawada, A.; Higashino, T.; Oyabu, T.; Takei, Y.; Nanto, H.; Toko, K. Gas sensor characteristics for smoldering fire caused by a cigarette smoke. Sens. Actuators B Chem. 2008, 130, 88–93. [Google Scholar] [CrossRef]
Figure 1. The fire source localization scene.
Figure 1. The fire source localization scene.
Sensors 18 01954 g001
Figure 2. The fire source arrangement localization principle in the RRR frame.
Figure 2. The fire source arrangement localization principle in the RRR frame.
Sensors 18 01954 g002
Figure 3. The circular fire source arrangement solved by Algorithm A and Algorithm B1.
Figure 3. The circular fire source arrangement solved by Algorithm A and Algorithm B1.
Sensors 18 01954 g003
Figure 4. The circular fire source arrangement solved by Algorithm A and Algorithm B2.
Figure 4. The circular fire source arrangement solved by Algorithm A and Algorithm B2.
Sensors 18 01954 g004
Figure 5. The estimation error curves of the fire source by the Algorithms B1, C, and D.
Figure 5. The estimation error curves of the fire source by the Algorithms B1, C, and D.
Sensors 18 01954 g005
Table 1. The simulation settings.
Table 1. The simulation settings.
Simulation BackgroundSimulation Setting
The length of the buildings11 m (meters)
The width of the buildings11 m
the real point of fire source(5 m, 5 m)
The distance between the two temperature sensor arrays L = 9   m
The distance between the two sensors in a array d 0 = 0.1   m
The coordinate of the reference node of sensor array A(0.5 m, 1 m)
The coordinate of the reference node of sensor array B(1 m, 10 m)
The sampling frequency500 Hz
Table 2. The abbreviations of the algorithms used in this section.
Table 2. The abbreviations of the algorithms used in this section.
The AlgorithmsAbbreviations
The angle bisector method with the circum-circleAlgorithm A
The dynamic optimization localization method in the RPR frameAlgorithm B1
The dynamic optimization localization method in the RRR frameAlgorithm B2
The localization method based on VB-ASCKF in the RRR frameAlgorithm C
The localization method in the RRR frame with clustering technologyAlgorithm D
Table 3. The circle centers of real fire point and its estimates by Algorithms A and B1.
Table 3. The circle centers of real fire point and its estimates by Algorithms A and B1.
The Real Fire PointAlgorithm AAlgorithm B1
(5, 5)(4.849, 4.935)(4.873, 4.98)
Table 4. The circle centers of real fire points and its estimates by Algorithms A and B2.
Table 4. The circle centers of real fire points and its estimates by Algorithms A and B2.
The Real Fire PointAlgorithm AAlgorithm B2
(5, 5)(5.027, 5.144)(5.027, 5.133)
Table 5. Mean estimation error of the fire source by the Algorithms B1, C, and D.
Table 5. Mean estimation error of the fire source by the Algorithms B1, C, and D.
The AlgorithmsAlgorithm B1Algorithm CAlgorithm D
Mean estimation error0.91.21.7

Share and Cite

MDPI and ACS Style

Wang, G.; Feng, X.; Zhang, Z. Fire Source Range Localization Based on the Dynamic Optimization Method for Large-Space Buildings. Sensors 2018, 18, 1954. https://doi.org/10.3390/s18061954

AMA Style

Wang G, Feng X, Zhang Z. Fire Source Range Localization Based on the Dynamic Optimization Method for Large-Space Buildings. Sensors. 2018; 18(6):1954. https://doi.org/10.3390/s18061954

Chicago/Turabian Style

Wang, Guoyong, Xiaoliang Feng, and Zhenzhong Zhang. 2018. "Fire Source Range Localization Based on the Dynamic Optimization Method for Large-Space Buildings" Sensors 18, no. 6: 1954. https://doi.org/10.3390/s18061954

APA Style

Wang, G., Feng, X., & Zhang, Z. (2018). Fire Source Range Localization Based on the Dynamic Optimization Method for Large-Space Buildings. Sensors, 18(6), 1954. https://doi.org/10.3390/s18061954

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop