Nonparametric Change Point Diagnosis Method of Concrete Dam Crack Behavior Abnormality

The study on diagnosis method of concrete crack behavior abnormality has always been a hot spot and difficulty in the safety monitoring field of hydraulic structure. Based on the performance of concrete dam crack behavior abnormality in parametric statistical model and nonparametric statistical model, the internal relation between concrete dam crack behavior abnormality and statistical change point theory is deeply analyzed from the model structure instability of parametric statistical model and change of sequence distribution law of nonparametric statistical model. On this basis, through the reduction of change point problem, the establishment of basic nonparametric change point model, and asymptotic analysis on test method of basic change point problem, the nonparametric change point diagnosis method of concrete dam crack behavior abnormality is created in consideration of the situation that in practice concrete dam crack behavior may have more abnormality points. And the nonparametric change point diagnosis method of concrete dam crack behavior abnormality is used in the actual project, demonstrating the effectiveness and scientific reasonableness of the method established. Meanwhile, the nonparametric change point diagnosis method of concrete dam crack behavior abnormality has a complete theoretical basis and strong practicality with a broad application prospect in actual project.


Introduction
China has the most water-conservancy and hydropower projects currently in the world.At the beginning of the 21st century, China entered the peak period of dam construction.There is a group of 320 m-level ultrahigh concrete arc dams completed or under construction in Southwest China.Concrete dam crack and its influence on the structure safety have always been a hot spot and difficulty in the discipline of hydraulic structure.Currently, the study is focused on whether the macrocrack of concrete structure will be expanded under load, which is an important issue for the evaluation of safety performance of concrete structure working with crack [1].In dam safety monitoring field, as the masses have higher and higher requirements to the safety of water conversancy projects, the problem of concrete dam crack has been more and more highlighted.It has also gradually become a hot spot and difficulty in the study on safety of hydraulic structure to study the transformation of crack behavior in concrete dam due to crack expansion and even its abnormality to cause influence on the safety of concrete structure.Such problem is defined as concrete dam crack behavior abnormality in this paper, and the emphasis is put on the study of diagnosis method of crack behavior abnormality, so as to diagnose the crack abnormality timely and rapidly.
For the hydraulic concrete structure working with crack, Li et al. [2], in combination with mutation theory created the grey cusp mutation model of hydraulic concrete structure crack abnormality.Bao and Yu [3], based on grey differential fitting modeling, created the grey state equation of concrete dam crack, analyzed the relation between the stability of grey state equation and crack subcritical expansion, and proposed the grey system method to distinguish and judge the subcritical expansion of crack.Given that concrete dam crack system is an open, dissipated, and complex nonlinear dynamical system, Gu et al. [4] proposed dynamical fuzzy crosscorrelation factor index method to preliminarily explore the dynamic abnormality diagnosis of concrete dam crack behavior.The current study plays an important promotion role in expanding the study thought and promoting discipline crossing, and development and forecasts that applied mathematics theory will have a strong development potential and broad development space in this field.
Deep analysis shows that all abnormality diagnosis methods above are based on the sequence of concrete dam crack mouth opening displacement (CMOD), but, according to numerous researches and investigations, the author finds that whether the CMOD can be taken as the basis for diagnosis of crack behavior abnormality is not deeply studied in the existing literatures, and most of the literatures are based on the fact that CMOD generally is the only monitoring item of concrete dam crack, causing that the existing concrete dam crack abnormality diagnosis method is lack of solid theoretical foundation.Based on this, the author, according to the fracture mechanics fictitious crack model and linear asymptotic superposition assumptions, creates the double crack tip opening displacement criterion of concrete crack behavior abnormality, providing a theoretical basis for the use of CMOD for crack behavior abnormality diagnosis in actual project [5].Meanwhile, the basis for the existing concrete dam crack behavior abnormality diagnosis method is the parametric statistical model of safety monitoring on concrete dam crack [6], which is established on the basis of this statistical model and in combination with certain or some applied mathematics theories and emphasizes the study from the aspect of information analysis, and it is still required to have a deep analysis on basic issues such as scientificalness and reasonableness of the abnormality diagnosis method.
The parametric statistical model of safety monitoring on concrete dam crack takes hydraulic pressure, temperature, and timeliness as the main influential factors, in which the theoretical basis of crack timeliness component expression is the creep theory of concrete and the component of cyclic load of temperature and so forth is fitted with periodic function, and it expresses the hydraulic component with cubic or quartic polynomial continuous function of water head based on the theory of engineering mechanics [7].So the parametric statistical model of safety monitoring on concrete dam crack is widely applied in actual project as it can reflect the evolution law of crack behavior truthfully in a certain extent.However, when the abnormality diagnosis of crack behavior is based on the parametric statistical model of safety monitoring on concrete dam crack, as the functional relationship between the influential factors of crack and the crack effect size is set in advance and the stochastic disturbance of model needs to meet the normal consumption, the effect of crack behavior abnormality diagnosis is not reasonable and comprehensive enough.On this basis, closely combining the hot issue in the latest statistics, change point theory, based on the analysis of change point characteristics of concrete dam crack behavior and through the reduction of change point issue, establishment of basic nonparametric change point model, and progressive analysis of basic change point test method as well as the study on multichange point diagnosis method of crack abnormality, the author creates the nonparametric change point diagnosis method of concrete dam crack behavior abnormality in this paper and applies it in the actual project, so as to inspect the scientific reasonableness of nonparametric change point diagnosis method of concrete dam crack behavior abnormality.

Analysis on Change Point Characteristics of Concrete Dam Crack Behavior Abnormality
Crack behavior abnormality is representation that the crack evolves from stable state to instable state.The abnormity of concrete dam crack behavior is presented as the instability of model structure in crack monitoring model, that is, the dynamic relation between environment size and crack effect size changes suddenly in several time points for some reason or changes continuously and slowly in a certain time interval.For example, according to the existing crack behavior abnormality diagnosis method, the hydraulic component and temperature component in the CMOD are recoverable, while the timeliness component is unrecoverable.Thus, most of the abnormality diagnosis methods of crack behavior are against the timeliness component of CMOD.When the timeliness component meets some conditions, it is thought that the abnormality occurs to the crack behavior, for example, the grey cusp mutation model method [8].The general form of such methods can be described as first, the general expression of timeliness component of CMOD is defined as where {  } is the influence factor related with timeliness and  = [ 1 ,  2 , . . .,   ]  is the -dimensional regression parameter of timeliness component; for example, the expression of timeliness component of CMOD is generally (  ) =  1   +  2 ln   .
Then the abnormality of crack behavior in the mode is presented as the regression parameter  changes in some time points that is, in time point , the regression coefficient  1∼ determined by samples 1 ∼  is significantly different from the regression coefficient  +1∼ determined by samples  + 1 ∼  statistically.
The above crack behavior abnormality is diagnosed based on the timeliness component.However, generally the hydraulic pressure and temperature components can also reflect the abnormality of crack behavior.For example, under the same water level and temperature conditions, when the concrete dam crack depth increases, the annual variation amplitude of hydraulic pressure component and temperature component in CMOD will be changed.In the statistical model of CMOD, it is reflected that the regression coefficients of hydraulic pressure component and temperature component are changed; when this change is significant statistically, the crack behavior is abnormal.
For the two given CMOD sample sequences, in order to test whether the crack is abnormal when generating the two sample sequences, the statistical two-sample analysis method can be adopted.In statistics, two-sample analysis method is mainly to analyze whether the probability distribution characteristics of two samples are significantly different statistically, for example, the Chow test [9].For concrete dam crack, assume that, in its monitoring model, the endogenous variable is , exogenous variables are  1 ,  2 ,  3 , . . .,   , their  sequence observation values   ,  1 ,  2 ,  3 , . . .,   ( = 1, 2, . . ., ) are from two times with different environments, in which the first  0 ( 0 < ) observation values are from a certain historical period and the rest  −  0 ones are from another period.At this moment, if the data of the pre- 0 period, and those of the post- −  0 period are, respectively, used to establish two linear models between  and  1 ,  2 ,  3 , . . .,   :  1 =  (1)  0 +  (1)  1  1 +  (1)  2  2 + ⋅ ⋅ ⋅ +  (1)     +  1  , 1 =  (2)  0 +  (2)  1  1 +  (2)  2  2 + ⋅ ⋅ ⋅ +  (2)     +  2  , Obviously, if the intercepts of  (1)  0 and  (2)  0 are different or the coefficients of a certain exogenous variable are different, the structural relation between  and  1 ,  2 ,  3 , . . .,   is instable.On the contrary, if the intercepts and coefficients of all exogenous variable are the same, the two formulas are the same model in nature.It is called that there exists a stable structural relation between  and  1 ,  2 ,  3 , . . .,   .On this basis, it is only needed to establish one model between all the  observation values between  and  1 ,  2 ,  3 , . . .,   : It can be seen that the study on structural stability of the model can be converted into the statistical test issue of model parameter in the following form: In order to test whether the different hypothesis on the two models is true, start from the original hypothesis that the two regression equations are the same, and then test whether the hypothesis can be refused.If the ordinary least squares method is used to estimate the models (2) and (3), respectively, as there is no limitation condition to the model coefficients, the sum of residual sum of squares of the two equations is used to represent the unconditional residual sum of squares, that is, ESS UR = ESS 1 + ESS 2 , in which ESS 1 and ESS 2 , respectively, represent the residual sums of squares of models (2) and (3), they are respectively subject to the  2 distribution of degrees of freedom of  0 − ( + 1) and  −  0 − ( + 1), and ESS UR is subject to the  2 distribution of degree of freedom of  − 2( + 1).Then the ordinary least squares method is used to evaluate the model (4) and calculate the residual sum of squares ESS  with condition limit, which is subject to the  2 distribution of degree of freedom  − ( + 1).If we assume that  0 is true, then ESS  will not exceed ESS UR too much.On this basis, the hypothesis  0 can be judged whether it is true according to whether the residual sum of squares ESS = ESS  − ESS UR is significant; the test statistics constructed is Obviously, the above test statistics is subject to the  distribution of first degree of freedom +1 and second degree of freedom  − 2( + 1).If the structure between  and  1 ,  2 ,  3 , . . .,   is stable, the parameters estimated in the models (2), (3), and (4) should be approximate.Thus the difference between ESS UR and ESS  is tiny; at this moment,  tends to be the smaller value; otherwise,  will be larger.Based on this, the structural stability of the model can be judged accordingly for the significance level  given; if  <   ( + 1,  − 2( + 1)),  0 is accepted; at this moment, it can be thought that the model structure is stable; that is, no abnormality occurs to the crack behavior at this moment; on the contrary, if  >   ( + 1,  − 2( + 1)),  0 is refused, and thus the structure between  and  1 ,  2 ,  3 , . . .,   is instable; that is, abnormality occurs to the crack behavior at this moment.
It can be seen that when the abnormality point of concrete dam crack behavior is known, the test analysis method of two-sample issue can be adopted to diagnose whether abnormality occurs to the crack behavior.But in practice, it is difficult to know the moment when abnormality occurs to the crack behavior.When the abnormality point is unknown, the issue on crack behavior abnormality diagnosis is converted into change point issue.
The main study of change point issue is to test and analyze whether the statistical characteristics of a given time sequence are significantly changed at an unknown moment.Statistically, change point  0 means that some statistical characteristics of the sequence or process change at an unknown moment  0 .Assuming a sample sequence  1 ,  2 , . . .,   , the corresponding distribution functions are  1 ,  2 , . . .,   .If there exists  0 , which makes   0 +1 that differs greatly from   0 in some characteristics (such as mean value and position parameter etc.), then  0 is called a change point of the sequence.More generally, if  1 ,  2 , . . .,   are divided into  groups, { 1 , . . .,   1 }, {  1 +1 , . . .,   2 }, . . ., {   +1 , . . .,   }, with sample distribution in each group relatively stable, while there are mutations at  1 ,  2 , . . .,   , then  1 ,  2 , . . .,   are the  change points of the sequence.
Change point is ubiquitous in nature and society, which reflects the process from quantitative change to qualitative change of things.The issue of change point, whether in theoretical exploration or actual data analysis, is very important, that is, to analyze whether the potential mechanism is changed, when it is changed, and how it is changed based on a group of data; if these problems are ignored, the accuracy of estimation and forecast will be reduced.
For the monitoring sequence of an effect size of concrete dam crack, for example, the CMOD monitoring sequence  1 ,  2 , . . .,   , whether abnormality occurs to the crack behavior can be expressed as the original hypothesis is true; that is, whether the distribution of CMOD sequence is changed somewhat, while the distribution change forms are diversified in theoretical exploration and practice.The distribution of observation value is changed according to a certain law and is changed to another law at an unknown moment, and this moment is the change point.Or the moment of specific change cannot be found out at all, and only the distribution law of observation is changed progressively just at the beginning.It is also possible that the sample distribution relies on a spatial parameter, while this mutation or progressive change relies on the position or interface in the space.Another possibility is that the distribution law of samples will be changed with change of system condition or an element in the system, and then the mutation or progressive change will rely on a specific nontime and space variable.Generally, the "structural change" in the parametric statistical model is often represented as the change of specific parameter .Another "structural change" is not completely decided by the change of some parameters, such issues are of diversified forms and their models are also relatively complex.The change point study on concrete dam crack behavior abnormality in this paper is mainly to detect whether the statistical behavior of the crack is significantly fluctuated, especially to detect whether the crack fluctuation exceeds its safety control range.When the crack is changed qualitatively (mainly exceeding the warning line of crack control), it is hoped that early warning can be released in time to avoid safety accident of the dam.

Construction of Concrete Dam Crack Behavior Change Point Diagnosis Method
It can be known from the above analysis on characteristics of concrete dam crack behavior abnormality that whether abnormality occurs to the crack behavior can be diagnosed with change point theory.According to the prior information depended, the existing change point theory diagnosis methods can be mainly classified into [10] parametric, semiparametric, and nonparametric diagnosis methods.In the study of change point issues, due to the uniqueness of nonparametric statistics; that is, no background of sample distribution is required, in recent years, the nonparametric methods in study on change point issue have been very active and their contents are also very rich.This paper is focused on the nonparametric change point study of concrete dam CMOD sequence, so as to establish the nonparametric diagnosis method of concrete dam crack behavior abnormality.

Reduction of Change Point
Issue.Any characteristic of the random sequence is determined by its probability distribution, that is, by the probability measure in a functional space.Therefore, the change of any characteristic of the random sequence is the change of its probability distribution generally.Based on this, it can be thought that the random sequence with the change point is "bonded" by the sequence fragments with different probability distributions [10], while the any-probability distribution of random sequence is decided by a limited-dimension distribution set meeting the conditions of consistency [11].Therefore, any change of the probability characteristics is the product of some change of any limited-dimension distribution.Consider the following single change point model: assume  = {  } ∞ =1 and  = {  } ∞ =1 are two strictly stationary random sequences, command 1 ≡  0 <  1 <  2 < ⋅ ⋅ ⋅ <   ,  < ∞, and define Assume   (⋅) ̸ =   (⋅) and consider the sequence  = {  } ∞ =1 , in which where {  } is any sequence.Then it is thought that  is the sequence with a change point; that is, its There are varieties of change point issues, and the different probability characteristics of stationary random sequence can be defined as the corresponding change point issues.Thus, only with the help of all finite-dimension distributions can the random sequence be fully described.Therefore, for the purpose of theoretical analysis and practical application of the change point issue, it is necessary to refine a basis change point type from so complex and various change point issues, so that other types of change point issues can be reduced to such basic type conveniently.It is thought that sequence  with a change point is a basic type of change point if the mathematical expectation of sequence  changes, that is, where   ̸ =   and {  } is any sequence.The method to diagnose the change point of the formula ( 9) is a basic change point diagnosis rule.Generally, in this basic type of change point, other characteristic parameters of the random sequence also may be changed; here other parameters are called redundant parameters.Meanwhile, study shows that [10], as long as enough diagnosis sequence {   }, is introduced; the test of any distribution function parameter change in sequence  = {  } ∞ =1 can be reduced to a basic type of change point test; that is, the reduction of change point issue makes it possible to adopt the same change point test method to test any information involved in the characteristics of random sequence change point as long as the basic change point test rule is analyzed.
→ .Assuming { (1) ()} and { (2) ()} are two strict and stationary independence random sequences, when certain conditions are satisfied, sequence where where  1 and  2 are two standard and independent Wiener processes.
When  (1) ≡  (2) ,  1 ≡  2 ≜ , and  1 ≡  2 ≜ , () can be simplified as The above method is constructed in consideration of the basic type of change point.Meanwhile, the change point test method of basic change point problem mentioned above can be spread to the change point test problem of any probability characteristics.

Asymptotic Analysis of Basic Change Point Test Method.
Through the asymptotic analysis of first and second types of probability errors and estimated probability error, it is able to choose the most optimal adjustment parameter in the family to meet the requirement of change point test quality in practice.Assuming  = {()}  =1 is an observation of a random sequence, the following form of statistics family should be considered: where 0 ≤  ≤ 1 and  = 1, 2, . . .,  − 1.
When  = 0,   () is the statistics.When  = 1/2,   () is the maximum likelihood statistics constructed by Hinkly [12] for Gaussian sequence.When  = 1,   () is the nonparametric statistics adopted by Deshayes and Picard [13,14].Assume the observed sequence  is from the following generating scheme: where 0 <  ≤  ≤  < 1,  > 0,  ∈ R 1 , ℎ ∈ R 1 , and ℎ ̸ = 0, Ẇ() is Gaussian white noise, and () is the standard Wiener process.For continuous time sequence, the statistics   () can be expressed as where    The first type of error probability can be expressed as The second type of error probability can be expressed as And estimated error probability can be expressed as For the first type of error probability, the following equation exists: where is of monotone increasing on the interval of 0 ≤  ≤ 1.So the asymptotic optimal method to reduce the first type of error probability is to assume  = 1 and then lim  → ∞ (| ln ()|/) = 2(/) 2 .
For the second type of error probability, under certain conditions, the following formula exists: It can be seen that the second type of error probability satisfies Therefore, the asymptotic optimal method to reduce the second type of error probability is to assume  = 0, so () ∼ exp(−(ℎ 2 /2 2 )()(1 − ) 2 ). where and Namely, when  → ∞, where (, ) satisfies max 0≤≤1 min ≤≤ (, ) = ℎ 2 /8 2 ; namely, when  = 1/2, the estimated error probability is the minimum.Therefore, selection appropriate parameter  can improve the quality of change point test.
The nonparametric diagnosis method of crack behavior abnormality of concrete dam is constructed theoretically by the reduction of any change point problems, construction of basic nonparametric change point model, and asymptotic analysis of basic change point test method.And its basic ideas are first, to construct corresponding diagnosis sequence based on the initial time sequence of CMOD according to the reduction guidelines of change point problem; then, to construct statistics sequence according to the formula (11); finally, to diagnose whether the change point exists or not according to the formula ( 13) under the given level of significance test.
However, in practice, for the given CMOD sequence of concrete dam, it is hard to ensure that there is only one abnormality point upon crack abnormality, so the number of crack abnormality point is unknown.Therefore, the multiple change point diagnosis method of crack behavior abnormality of concrete dam is built combing with the asymptotic property of the statistics family below.

Multiple Change Point Diagnosis Method of Crack Behavior Abnormality.
For many existing change point test methods, if assuming the diagnosis sequence can be described by a mathematical model of limited parameter set, the methods belong to the scope of parametric change point test method [15][16][17].When the phenomenological model of study object is known, the parametric method is effective.However, when the phenomenological model has model error during setting, the result obtained by parametric change point test method is doubtful.Therefore, this paper constructs corresponding nonparametric change point test method to diagnose the crack behavior abnormality of concrete dam which has less demand of prior information of diagnosis sequence; that is to say, the nonparametric method does not require the prior information of probability distribution of diagnosis sequence [18].Furthermore, for the given level of significance test, this method can almost test the abnormality characteristics of random crack behavior.Assuming a random sequence observation ( = 1, 2, . . .,  + 1) , where The asymptotic analysis of basic change point test method shows that when  = 1, the first type of error probability is minimum; namely, the false alarm probability is minimum; when  = 0, the second type of error probability is minimum; namely, the false dismissal probability is minimum; when  = 1/2, the estimated error probability of change point is minimum; namely, the probability of change point estimated value far away from the actual change point position is minimum.Therefore, selecting different  value in different change point test phase can improve the precision and credibility of the estimated result of change point [19].The basic steps of multiple abnormality point diagnosis method of crack behavior abnormality of concrete dam are as follows.
(1) Construction of the diagnosis sequence from the original CMOD sequence: construct diagnosis sequence from original signal to make the change of various characteristics of original signal convert into the change of diagnosis sequence mathematical expectation.
( In the formula (29),  is the threshold value of crack behavior abnormality point test which relies on the specified false alarm probability, sample size, sample variance, and correlation coefficient of sample where the false alarm probability is specified in advance, the smaller this parameter value, the more abnormality point will be of missing test.In this phase, to not miss testing any possible abnormality point, the false alarm probability should be bigger.
(3) Initial estimation of crack behavior abnormality point: the point  1 that makes the statistics |  (, 1)| global maximum can be used as the estimation of first abnormality point.This point divides the diagnosis sequence into two subsample set including  1 : where  is a smaller value relying on the sample size and the maximum value in statistics.Then, test the homogeneity of each subsample  1 and  2 by the method in step (2) and estimate the abnormality point by the method in step (3).Repeat above steps until each sample in the subsample set is homogeneous.A sequence of initial estimation  1 ,  2 , . . .,   of abnormality points is obtained by steps (2) and (3) where  is the initial estimation number of abnormality points.
(4) Eliminating suspicious crack behavior abnormality point: redivide the diagnosis sequence to make each subsample set only contain an abnormality point   of initial estimation, that is to say : : Conduct homogeneity test of each subsample set according to step (2).However, the false alarm probability should be smaller.If a certain subsample set is tested to be homogeneous, the corresponding abnormality points in the subsample set should be eliminated.
(5) Final estimation of crack behavior abnormality point: after step (4), for the remaining subsample set   with a size of   , calculate the statistics   (, 1/2), and the corresponding point of maximum value of the statistics is the final estimation of the  abnormality point.In the above abnormality point test method, the determination of threshold  is the key.To calculate the threshold , the homogeneity assumption of sequence should be workable.Then, according to the functional limit theorem of random sequence [20], the following formula exists under certain conditions lim where  makes the sequence When false alarm probability  f.a. is given, threshold  can be calculated by the following formula: It can be seen from the above analysis that the threshold  relies on  and  is unknown.To solve this problem, Piryatinska et al. [21] adopts the following method to estimate this threshold.First, when the false alarm probability  f.a. is given, assuming  = 1, the initial estimation of threshold   can be obtained according to the formulae (34) and (35).Then, divide the diagnosis sequence into two subsample sets including  1 : 1 ≤  ≤  1 − [] and  2 :  1 + [] ≤  ≤  by steps (2) and (3), calculate the mean value of each subsample set, and deduct the corresponding mean value from the subsample set to form a new sample set by merging.Finally, calculate the empiric variance of the sample set as the estimation σ2 of  2 .So the estimation of threshold  is It can be seen from the basic principle of the method that when there is no or only one abnormality point in the diagnosis sequence, the estimation σ2 is more precise and the estimation Ĉ is much reliable.However, when there are multiple abnormality points in the diagnosis sequence, the error of estimation σ2 obtained by this method is larger and the estimation Ĉ is suspicious.Meanwhile, it is assumed initially that  = 1 can affect the threshold estimation Ĉ to a certain degree.
In consideration of the problems existed, this paper improves this method and its basic ideas are under an assumed , estimate the abnormality points in the diagnosis sequence, then use empiric variance estimation  =  .according to similar method, and obtain corresponding threshold estimation C according to the formula (36).Then analyze the points of diagnosis sequence again by this threshold, in this cycle, until the abnormality point estimation for two times is the same,  estimation satisfies σ =  ., and threshold estimation is Ĉ = C.This improvement can eliminate the influence of multiple abnormality points and assumed  on threshold to a great extent.1), so as to diagnose the abnormality conditions of crack behavior.Upon the analysis of given CMOD sequence of concrete dam by adopting nonparametric change point diagnosis criteria, construct diagnosis sequence according to the original sequence first.And the CMOD means and amplitudes are the two essential characteristics where the mean value presents the base point of change of CMOD, namely, the average level of CMOD, and the amplitudes represent the amplitude of variation of crack under load in one cycle.Therefore, when there is change point in CMOD means or amplitudes, the crack of corresponding moment is deemed as abnormality.Based on this, this paper analyzes the crack behavior abnormality according to the diagnosis sequence of CMOD means, amplitudes, and correlation functions.

Project Case
The diagnosis sequence based on CMOD means can be expressed as The diagnosis sequence based on CMOD amplitudes can be expressed as where  is time delay.The diagnosis sequence based on CMOD correlation functions can be expressed as   = CMOD  CMOD + ,  = 1, 2, . . .,  − ;  = 0, 1, . . .,   , where   is time delay.The CMOD sequence and the diagnosis sequence of its corresponding means, amplitudes, and correlation function in this case should refer to Figures 2, 3, and 4.
This paper analyzes the diagnosis sequence according to the nonparametric change point diagnosis method of crack behavior abnormality built where Figure 5 is the diagnosis process of crack behavior abnormality based on the diagnosis sequence of CMOD means which is similar with the abnormality diagnosis process of diagnosis sequence based on CMOD amplitudes and correlation functions and not repeated here.Tables 1, 2, and 3 are the abnormality diagnosis results of diagnosis sequence based on CMOD means, amplitudes, and correlation functions.
For the diagnosis sequence of CMOD means, it can be seen from Figure 2  As analyzed above, the nonparametric change point diagnosis method of crack abnormality built in this paper can diagnose the abnormality condition of crack behavior.The diagnosis result not only contains obvious abnormality point of crack behavior, but also is of good diagnosis effect on the abnormality point of crack behavior that is not easy to find.

Conclusion
(1) Whether the concrete dam crack expands or not is an important subject of safety performance evaluation of concrete dam, and the study on the abnormality diagnosis of crack behavior of concrete dam has become the hot spot and difficult problems in hydraulic structure safety.However, as the functional relationship between influencing factor and effect size of crack is set in advance and the random noise item of model should satisfy the normality assumption, the abnormality diagnosis result of crack behavior based on the safety monitoring parametric statistics model of concrete dam crack is not reasonable enough and comprehensive.
(   (4) The nonparametric diagnosis method of crack behavior abnormality of concrete dam can be applied in practical engineering.The diagnosis result contains obvious abnormality point of crack behavior and has good diagnosis effect of crack behavior abnormality point that is not easy to find, demonstrating the effectiveness and scientific rationality of built method.Meanwhile, the nonparametric change point diagnosis method of crack behavior abnormality of concrete dam is of complete theoretic basis, strong practicability, and wide application prospect in practical engineering.
(5) In accordance with the concept of dam safety monitoring, the research in this paper and many other researches on the abnormality diagnosis of crack behavior in concrete dam belong to the a posteriori category; that is, the crack monitoring data is regarded as given and the sample size is known and statistical central limit theory is used to diagnose crack behavior abnormality.However, in order to achieve the real-time monitoring and early warning of dam safety in practice, the concern about crack behavior abnormality of concrete dam is whether the current crack behavior is abnormal and whether the former monitoring model of crack safety is suitable for the current status, which corresponds to online diagnosis of crack behavior abnormality.The crack monitoring data is regarded as changing and the sample size is unknown in the online diagnosis.This will be the main research content of the authors and the research group in the future.

10 Figure 1 :
Figure 1: Time sequence curve of the measured CMOD.

)
The internal relationship between the crack behavior abnormality of concrete dam and statistics change point theory is deeply analyzed to make the abnormality diagnosis method of crack behavior of concrete dam based on change point theory have complete theoretical basis from the perspective of structural instability of parametric statistics model and change in the distribution law of sequence of nonparametric statistics model based on the representation of crack behavior abnormality of concrete dam in parametric statistics model and nonparametric statistics model.(3)The test of parameters change of any distribution functions is reduced as the basic type and method of
2) Homogeneity test of diagnosis sequence: calculate the statistics   (, 1), 1 ≤  ≤  − 1 from diagnosis sequence.If Conduct abnormality diagnosis of the crack near 105 m elevation at the downstream face of a gravity arch dam based on the nonparametric change point diagnosis method of crack behavior abnormality of concrete dam.This dam is a concrete gravity arch dam of concentric variable radius with a crest elevation of 126.3 m and maximum dam height of 76.3 m and has 28 dam sections from left to right.It takes 12 years from dam pouring to substantial completion and is constructed by three phases.During Phase II concrete pouring, as the level rising speed is fast and interval of concrete lift is short, the contraction distortion of Phase II concrete is strongly constrained by Phase I concrete to make the top of Phase I concrete have cracks from 5 # dam block to 28 # dam block with a length of more than 300 m.So, the dam blocks of 14 # ∼ 20 # were of cross-seam joint bar grouting in 1973 and of epoxy resin grouting in 1987.This paper analyzes the 1680 CMOD monitoring data of joint meter inbuilt in 18 # dam block from October 10, 1974 to December 18, 2006 (as shown in Figure that the CMOD means have two obvious change points during the period ofJanuary 15, 1976to November 25, 1976 and September 26, 1985 to August 28, 1986, that is to say, the CMOD means have abnormality.Meanwhile, the CMOD means have increased since 1987; that is to say, the CMOD means increase gradually.However, it can be seen from Table1that upon the abnormality analysis of the diagnosis sequence of whole means, the abnormality point of CMOD means diagnosed first occurred on August

Table 1 :
Abnormality diagnosis results of diagnostic time sequence of the CMOD means.

Table 2 :
Abnormality diagnosis results of diagnostic time sequence of the CMOD amplitudes.

Table 3 :
Abnormality diagnosis results of diagnostic time sequence of the CMOD correlation functions.