Research on early-warning index of the spatial temperature field in concrete dams

Warning indicators of the dam body’s temperature are required for the real-time monitoring of the service conditions of concrete dams to ensure safety and normal operations. Warnings theories are traditionally targeted at a single point which have limitations, and the scientific warning theories on global behavior of the temperature field are non-existent. In this paper, first, in 3D space, the behavior of temperature field has regional dissimilarity. Through the Ward spatial clustering method, the temperature field was divided into regions. Second, the degree of order and degree of disorder of the temperature monitoring points were defined by the probability method. Third, the weight values of monitoring points of each regions were explored via projection pursuit. Forth, a temperature entropy expression that can describe degree of order of the spatial temperature field in concrete dams was established. Fifth, the early-warning index of temperature entropy was set up according to the calculated sequential value of temperature entropy. Finally, project cases verified the feasibility of the proposed theories. The early-warning index of temperature entropy is conducive to the improvement of early-warning ability and safety management levels during the operation of high concrete dams.

are critical to perform safe operations and determine crucial topics in dam safety monitoring field research.
To date, the warning methods of dam body's temperature have obtained a series of research achievements. Temperature early-warning index (Su et al. 2011(Su et al. , 2012Zhou et al. 2012) is an important index to control cracks of concrete dam. When the current concrete temperature exceeds early-warning index of temperature, effective temperature control measures must be taken to control the concrete temperature. In 2004, Wang took some temperature measuring points as constraint conditions, using three-dimensional finite element analysis method to represent the characters of the dam's temperature field (Wang 2004). Wu applied sensitivity analysis method to determine temperature field's main impact factors (Wu and Song 2011). Combined with the situation of the practical concrete dam temperature control, Qu et al. applied expert evaluation and the theory of entropy weight to establish a multi-objective fuzzy mathematical model for the concrete dams (Qu et al. 2012). Although the above-mentioned theories and methods complement and improve traditional methods in solving difficult problems in dam safety, but traditional warnings theories are targeted at a single point which have limitations, and the scientific warning theories on the global behavior of the spatial temperature field are non-existent. Thus, a scientific and accurate early-warning index based on the spatial temperature field should be studied to improve the warning ability of the concrete dams.
In this paper, accordingly, the spatial temperature field at different elevations and positions have different behaviors, that is, the regional dissimilarity exist. Through Ward spatial clustering, the spatial temperature field was divided into regions. On this basis, an expression of temperature entropy was proposed based on the synergetics and entropy. This expression can comprehensively evaluate the overall variation of temperature field in concrete dams. With this expression, the sequential value and the early-warning index value of temperature entropy was analyzed and determined via the small probability method. Finally, project cases verified the feasibility of the proposed theories.

Research on the partitioning method of the spatial temperature field
The traditional partitioning method adopts the average value of the monitoring value time series. This method can only express the average variation condition of the monitoring values. Thus, information on the time dimension is considerably lost and an invisible hypothesis is generated, that is, the measuring points experience the same directional variation in the time dimension. Reflecting the real change rule of the structural behavior with time is difficult and unreasonable. The development law of the monitoring values of the measuring points in different periods is described in Fig. 1. If the traditional method is used, then measuring points 1 and 3 should be classified as one category. However, if the entire variation process is considered, then measuring points 2 and 3 should be classified as one category. If no hypothesis is formulated for the spatial temperature field, then partitioning is conducted according to the intimacy degree among all the observed values and based on Ward spatial clustering (Bertamini et al. 2016;Hu et al. 2016;Nascimento et al. 2011;Li et al. 2013). The partitioning principle attempts to make the similarity degree of the observed value change rule within the partitions considerably close while making the similarity degree of the observed value change rule between the partitions substantially low. The spatial temperature field partitioning needs to process two core problems: (1) selecting the type of statistical magnitude to characterize the similarity degree among the opening displacements of the measuring points and (2) selecting the criteria that should be used to determine the similarity between the partitions.

Measurement method of similarity degree among the monitoring values of the temperature measuring points
N is assumed as the number of all temperature monitoring points, T is the monitoring time series, and x it (i = 1, 2, . . . , N ; t = 1, 2, . . . , T ) represents the temperature data set. For set x it , S t is the standard deviation at time t, and the following three distances are defined to describe the similarity degree among the monitoring values of the temperature measuring points.
The "absolute distance" between measuring points m and n is simplified as d mn (AQED) and is expressed as follows: where d mn (AQED) characterizes the distance between measuring points m and n within the entire period.
The "incremental distance" between measuring points m and n is simplified as d mn (ISED) and is expressed as follows: where d mn (ISED) depicts the difference in variation trends of the monitoring values of the temperature measuring points m and n. If the monitoring values presents samedirectional change across time, then change is substantially coordinated, the two values will be similar, and d mn (ISED) will be considerably small. If the monitoring values presents reverse-directional change across time, then the similarity is poor and d mn (ISED) will be substantially high.
The "comprehensive distance" between measuring points m and n is simplified as d mn (CED) and is expressed as follows: The "comprehensive distance, " which is a weighted array of "absolute distance" and "incremental distance, " comprehensively describes the similarity of change of the measured values at the measuring points.

Measurement method of similarity degree among the different spatial temperature field partitions
All of the measuring points are assumed to be divided into k partitions. These partitions are G 1 , G 2 , . . . , G k , N l is the number of measuring points of the G l category, X l is the center of gravity of the measuring values of the G l category, and X il is the monitoring value of the i(i = 1, 2, …, N l ) th measuring point in the G l category. For the monitoring data of the N measuring points within T periods, the sum of the squares of deviations of the different measuring point series in the G l partition of the dam is expressed as follows: The sum of the squares of deviations of the k partition is expressed as follows: where W * l is the sum of the squares of deviations of the N l measuring values, X it is the monitoring value of the measuring point i in period t, Y it = X it X it−1 represents the relative increment of the monitoring values of the temperature measuring point i within the t period in G l , X it = X it − X it−1 represents the absolute quantity difference of the monitoring values of the temperature measuring point i in the t and t − 1 periods in G l , The assumption is that n combinations are conducted in the partitioning process based on the threshold value method, and the ratio of the distance between partitions in the lth partitioning to that in the last partitioning is S l = D l D n−1 . If the difference between S l and S l+1 is small and that between S l and S l−1 is large, then the corresponding distance D l between the partitions can be the threshold value of the partition.

Partitioning flow
This study employed Ward spatial clustering as basis to propose a complete partitioning flow. The procedure involved is as follows.
Step 2: The calculated coefficient values are substituted into expression (5), and the comprehensive distance d mn (CED) between every two measuring points among the N measuring points is calculated.
Step 3: Initially, all measuring points self-form a partition, the number of partitions is k = N, distance matrix D (1) between partitions is built, and the ith partition is Step 4: According to the principle of the minimum sum of the squares of deviations, two partitions with the minimum comprehensive distance are combined as a new partition, the comprehensive distance d ij (CED) between the new partition and other partitions is calculated, a new distance matrix is obtained, and Steps 4 and 5 are repeated until the partitioning ends.
Step 5: The optimal partition combination of the measuring points is obtained based on the threshold value method. Thereafter, the optimal number K of partitions is obtained.

Temperature entropy
Previous studies have reported that the temperature field changes gradually, whereas the concrete structure has an obvious nonlinear self-organization mechanism and a multiscale coupling effect during the operation of concrete dams (Andrie and Chen 1975;Xie 2004;Bai 2008). The macro-mechanical properties of the concrete structure have a collaborative self-organized phenomenological response of internal multiscale physical quantities to various effects. Currently, the multiscale synergic evolution of a system is difficult to elaborate because of the limited theoretical basis and poor calculation methods. Although the system has different evolution equations at different scales and levels, energy is a general physical quantity that could run across scales and levels. To date, some researchers have studied system evolution from the perspective of energy. Reference analyzed the macro physical significance of entropy based on the available energy of closed thermal systems (Yu 1995). Reference discussed the variation law of deformation energy and the mechanism of sudden energy changes during the gradual evolution of concrete dams (Wu and Guo 2010). Reference studied the analytical methods of multi-objective decision-making during concrete dam operation based on the entropy weight (Maken et al. 2014). The abovementioned studies imply that the orderly evolution of systems can be effectively studied from the perspective of energy.

Entropy features of the spatial temperature field
Entropy is a basic concept in thermodynamics that is used to express a state function of a system. This concept was introduced by Rudolf Clausius in the 1850s while discussing the Carnot cycle.
where dS is the differential of entropy of the state function, δQ represents the heat exchanged between the system and environment, rev indicates that the process is reversible, and T is the temperature of the heat source. Hence, entropy is the state function of a system. Currently, the physical significance of entropy is mostly based on the microscopic expression of entropy, which was proposed by Ludwing Boltzmann in the 1870s.
where k is the Boltzmann's constant and Ω is the amount of microscopic quantum states, that is, the probability of occurrence of macroscopic state. Boltzmann's definition states that the microscopic physical significance of entropy can be expressed as follows: entropy is system uncertainty and a measurement of the system disorder degree. According to entropy theory, when a spatial temperature field is in a more dangerous state, the system is in greater disorder and its entropy value is smaller.

Methods of characterizing contributions of single observation point
According to principle of synergetics, if the temperature monitoring points of high concrete dams are used as the characteristic points of system analysis, the whole temperature field can be expressed by evolution equations of all characteristic points. Because synergetics theory mainly studies that how to use system's own internal synergy to generate the ordered structure spontaneously in time, space and function. We can see that the value of single temperature monitoring points is out-of-order. However, if we take all of the temperature monitoring points as a whole, we can discover there similarity and relevance. This paper gave a quantitative description to the overall temperature of dam based on entropy structure. This is named as "temperature entropy". Temperature entropy is a comprehensive index that integrates the effects of various factors. This index can be used to describe degree of order of the whole temperature field. The hierarchical structure of temperature entropy of the spatial temperature field is shown in Fig. 2.

Constructing temperature entropy of single measured value
The degree of order for a single temperature monitoring point, a dimensionless standardized value, is introduced before calculating the single-point temperature entropy.
(10) S = k ln This value measures the degree of temperature order at single point. During the operation period, a concrete dam body is mainly influenced by the dead load and temperature load; consequently, the temperature variation can be viewed as a probability event (Wu 2003). For those complicated mechanical properties caused by materials, load, and environment and so on, they are reflected comprehensively in the measured temperature data from thermometers placed in the different part of dam. Therefore, when we proposed the temperature entropy expression that can describe degree of order of the temperature field in concrete dams, we didn't consider various mechanical and thermal properties. The small probability method for calculating the monitoring index discloses that the smaller the probability of temperature reoccurrence is, the more dangerous the dam will be. Therefore, the probability of temperature occurrence at monitoring points represents the safety degree of the temperature load and can be used as an important index to measure the degree of order of the temperature field.
Here, the degree of order of the measured temperature at each monitoring point can be defined. Generally speaking, the increasing dam temperature is viewed as a positive change. Therefore, the degree of order of the jth measured temperature at the ith monitoring point (μ ij ) can be defined as follows: 1. When temperature increases, u ij is 2. When temperature decreases, u ij is where f i (δ) is the probability density function of the measured temperature at the ith monitoring point and F(x) is the corresponding probability distribution of x.
Formulas (11) and (12) state that larger deviations of the measured temperature of a concrete dam from the initial temperature would lead to larger u ij values; otherwise, the u ij will be smaller.
According to the definition of information entropy, the temperature entropy of x ij can be defined as: where u 1 ij = u ij and u 2 ij = 1 − u ij are degree of order and degree of disorder of the monitoring point, respectively.
Obviously, Formulas (11) and (12) indicate that 0 < u 1 ij < 1 and 0 < u 2 ij < 1. Meanwhile, u 1 ij + u 2 ij = 1, which implies that the sum of the degree of order and the degree of disorder is 1. The degree of disorder decreases while the degree of order increases, and vice versa.

Weight entropy
The overall temperature entropy includes two levels: overall temperature entropy at the top and single-point temperature entropy at the bottom. The monitoring points interact with each other and influence the evolution of overall temperature entropy together. Let {ω i |i = 1, 2, · · · , n}(n is the number of total monitoring points) represent weight distribution of all temperature monitoring points. Apparently, ω i is non-negative (ω i ≥ 0 ). Meanwhile, the weight index meets the normalization conditions: n i=1 ω i = 1. Therefore, the weight distribution entropy of all monitoring points is: The temperature monitoring point with greater contributions to the overall temperature entropy possesses the higher weight. According to the principle of synergetics, the synergic evolution equation involves both stable and unstable modals. Structural evolution is mainly determined by unstable modals. Therefore, only unstable modes dominate the orderly variation of the dam temperature when the overall temperature field change is under critical phase change. The stable modals slightly influence the orderly variation of the dam temperature.

Weight optimization based on projection pursuit
Given that weight can be viewed as contributions of different models to the overall system evolution, weight calculation has a very important role in analysis of temperature entropy. In order to determine the weight distribution of monitoring points, this paper explored the implied information from the spatial temperature field data of a dam by projection pursuit analysis (PPA) (Friedman and Turkey 1974;Flick et al. 1990;Kumar et al. 2008). As a high-dimensional data analysis method, PPA could identify data structures or characteristics in the projection space by projecting high-dimensional monitoring data into low-dimensional spaces. PPA can analyze weights of evaluation index values according to the global data distribution and local data concentrations. Data dimension reduction is the main goal of PPA. The analysis integrates m-dimensional sequential values of measured dam temperature ( y ij i = 1 ∼ n; j = 1 ∼ m ) into one-dimensional comprehensive projection value (G(i)) toward the projection direction of P: where P is the unit length-width vector, P = {p 1 , p 2 , . . . , p m }.
Given the standard sample set of the measured dam temperature, P determines the form of the projection index function (Q(P)). Consequently, Q(P) only changes while P changes. Therefore, the optimum projection direction can be estimated by calculating the maximum Q(P): Meanwhile, P must satisfy the following constraint: 0.5 is the dispersity of projections; n j=1 R − r ij · f R − r ij is a projection index function, which represents the local density of one-dimensional data points along P; ḡ(i) is the mean value of sequence {G(i), i = 1, 2, …, n}; R is the window radius of local density, which is calculated as 0.1 S G in this paper; r ij is the distance between projections (r ij = |G(i) − G(j)|); f(t) is a unit step function. When t ≥ 0, f(t) = 1; when t < 0, f(t) = 0.
By substituting the optimum P* into Formula (15), the projection of sample point could be gained. Subsequently, the projection is normalized to obtain the weights of the deformation monitoring points: where G * (i) is the optimum projection of the ith evaluation index and ω i is the weight of the ith evaluation index.

Construction of the temperature entropy of the spatial temperature field
The degree of order sequence ({u 1 ij }), the degree of disorder sequence ({u 2 ij }), the deformation entropy sequence ({S j i }), and the weight distribution entropy of the spatial temperature field were calculated by using Formula (11)-(14). Based on these equations, the formula to calculate the multi-point temperature entropy was deduced as follows. The degree of order and the weights of the monitoring points were combined. For the ith monitoring point, the contribution of its degree of order to the overall temperature entropy is ω i u 1 ij , whereas the contribution of its degree of disorder to the overall deformation entropy is ω i u 2 ij . According to definition of generalized information entropy, the orderly entropy of the whole temperature field is: Formula (19) can be rewritten as: Two items in Formula (19) can be simplified into: By substituting Formulas (20) and (21) into Formula (19), and the overall temperature entropy can be simplified into: Formula (22) reveals that the overall temperature entropy involves the weight distribution entropy (S j ω ) and the weighted average of the temperature entropy of different monitoring points ( n i=1 ω i S j i ).

Establishment of an early-warning index based on the temperature entropy
The measured temperature data of concrete dams collected from each survey were transformed through Formula (22) to obtain the sequential value of temperature entropy ({S i }). Combined with practical situations of dams, the most unfavorable temperature entropy for the load combination (S mi ) was chosen. S mi is a random variable. A sample space with N samples can be gained from monitoring sequence: Subsequently, the distribution of this sample space was tested by small-sample statistical testing methods (e.g., the A-D method or the K-S method) to determine the distribution function F(X) of its probability density f(x). K-S (Kolmogorov-Smirnov) method bases on empirical distribution function, which is used to determine whether a sample is from a specific distribution. A-D (Anderson-Darling) method is a correction of A-D method and gives weight to the distribution of the tail. Besides, K-S inspection has nothing to do with the specific distribution, that's to say, its critical value doesn't depend to the tested specific distribution.
S m is defined as the extreme value of the temperature entropy (Jaynes 1963). When S ≤ S m , the dam will show abnormalities or danger. The failure probability of the dam is: After the distribution of deformation entropy (S) was calculated, most attention focused on the determination of the failure probability (P α , herein referred to as α) to estimate S m . Based on α and f(x), the early-warning index of temperature entropy can be determined: The spatial temperature entropy index (S m ) can be determined as soon as the optimum fitting distribution function was identified (Fig. 3).
S m is the early-warning index value under the failure probability of α. Dam security can be evaluated by comparing the temperature entropy and S m . If temperature entropy is smaller or equal to S m , attention should be paid to determine probable causes, to strengthen the monitoring system, and to analyze whether other monitoring items have abnormalities. An appropriate α should be chosen to calculate S m . The value of α is determined by various factors, such as engineering grade, engineering scale, operation progress and so on. Meanwhile, several confidence values of α shall be set for the dam risk management: α 1 , α 2 , …, α n (α 1 > α 2 > ···α n ).The corresponding S m (S m 1 , S m 2 , . . . , S m n ) should be calculated to get the multi-level early warning index of the overall temperature field in dams.

Project overview
The studied high concrete dam is mainly composed of the roller compacted concrete gravity dam, the exterior overflow dam section, the discharge chute and stilling pool, the underport at right bank for flood discharge (sand flushing), the underport at left bank for sand flushing, the permanent buildings (e.g. dam-rear power plant and access tunnel), and the temporary buildings (e.g. diversion tunnel and cofferdam). The elevation at the top is 1424 m. The biggest height is 160 m and the top length is 640 m. This paper mainly studied the 6th dam block. The thermometer distribution is shown in Fig. 4.
To study the temperature field changes during dam operation and the safe impounding period, temperature monitoring data from August 1, 2007 to January 20, 2014 were used as the sample data in this paper. Temperature data of all monitoring points within the 6th dam block were used to calculate temperature entropy and early-warning index value.

Partitioning calculation results
All measuring points were finally partitioned into three major categories, and the dam was divided into three regions according to the locations of the measuring points. Variation trends and rules of the measuring points within the categories obtained through clustering were identical. Each category of the measuring points could comprehensively describe the features of the spatial temperature field. A temperature field partition map of the 6th dam block is illustrated in Fig. 5. Table 1 provides the measuring point clustering partition table of the spatial temperature field.

Calculation results of the weight ω i
Figures 6 and 7 reveal that the environmental temperature obviously influenced the dam body's temperature. Temperature change was divided into the stage of temperature rise   Table 2.
Calculation of spatial temperature entropy Figures 8,9,10 show the temperature entropy curve of all partitions. Table 3 shows the annual extreme temperature entropies of all partitions.

Calculation of early-warning index of concrete dam
In this paper, two-stage warning indicators were set according to the practical running of this project and danger: α = 5% is the primary warning that is mainly used to discriminate and handle early dangerous case, whereas α = 1% is the secondary warning that is mainly used to determine grave danger and prevent urgent danger. Table 4 shows the K-S test results of all partitions. Table 5 shows the parameter values of the probability density function of all partitions. Table 6 shows the early-warning index values of temperature entropy of all partitions.
In this paper, the early-warning index values of all partitions was analyzed through the theoretical method. If the temperature entropy value of partition I reaches 0.569363, the dam is in the state of primary warning. If the temperature entropy value of partition I reaches 0.381156, the dam is in the state of secondary warning. If the temperature entropy value of partition II reaches 1.04325, the dam is in the state of primary warning. If the temperature entropy value of partition II reaches 1.002615, the dam is in the state of secondary warning. If the temperature entropy value of partition III reaches 0.971805, the dam is in the state of primary warning. If the temperature entropy value of partition III reaches 0.876503, the dam is in the state of secondary warning.  Table 2 Weight table of