“…WhereY i 0 is measured value of surrounding rock deformation, Y i is calculated value of surrounding rock deformation,m is the number of measured value, x i is i-th parameter, n is number of the parameter, x i a and x i b is Upper and lower limits of x i. Y i is a complicated nonlinear implicit function which selects parameter of rock x i (x 1 ,x 2 ,….,x n )as variables,introduced global optimization method-Difference Evolution (DE) Algorithm to global optimize parameter x i .Calculation process of back analysis of DE -FEM Coupling which have been formed by imbeding the FEM in DE [6] is as follows:…”