Second-order adjoint equation method for parameter identification of rock based on blast waves in tunnel excavation

Abstract

The objective of this study is to present a method for identifying the elastic moduli of ground rock via the first- and second-order adjoint methods using blast vibration measurements during tunnel excavation. For identifying these parameters, the magnitudes of the blast force should be identified beforehand. Parameter identification is a minimization problem of the square sum of the discrepancy between the computed and observed velocities. The magnitudes of the three components of borehole pressure are assumed to be independent in each direction. The propagation of an elastic wave is assumed because the amplitude of such a wave is infinitesimal. The three-dimensional finite element method and the linear acceleration method are used effectively. The extended performance function can be expanded into a series of small constants to derive the necessary condition of minimization. The adjoint equation and the dynamic equation of motion can be used to obtain the gradient and the Hessian product of the extended performance function with respect to the parameters. The weighted gradient method and Broyden–Flecher–Goldfarb–Shanno method are successfully employed for the minimization. By applying the present identification technique at the Ohyorogi tunnel site, the fact that the computed and observed velocities are well in agreement is verified. The present method can be shown to be useful for tunnel excavation.

Keywords:

• parameter identification
• finite element method
• blasting wave
• first-order adjoint equation
• second-order adjoint equation

Acknowledgements

The computation was carried out using Dospara Prime 557848.