An Identification Method of Geological Boundary Using Elastic Waves

Authors: Masamitsu Chikaraishi, Mutsuto Kawahara

Abstract:

This paper focuses on a technique for identifying the geological boundary of the ground strata in front of a tunnel excavation site using the first order adjoint method based on the optimal control theory. The geological boundary is defined as the boundary which is different layers of elastic modulus. At tunnel excavations, it is important to presume the ground situation ahead of the cutting face beforehand. Excavating into weak strata or fault fracture zones may cause extension of the construction work and human suffering. A theory for determining the geological boundary of the ground in a numerical manner is investigated, employing excavating blasts and its vibration waves as the observation references. According to the optimal control theory, the performance function described by the square sum of the residuals between computed and observed velocities is minimized. The boundary layer is determined by minimizing the performance function. The elastic analysis governed by the Navier equation is carried out, assuming the ground as an elastic body with linear viscous damping. To identify the boundary, the gradient of the performance function with respect to the geological boundary can be calculated using the adjoint equation. The weighed gradient method is effectively applied to the minimization algorithm. To solve the governing and adjoint equations, the Galerkin finite element method and the average acceleration method are employed for the spatial and temporal discretizations, respectively. Based on the method presented in this paper, the different boundary of three strata can be identified. For the numerical studies, the Suemune tunnel excavation site is employed. At first, the blasting force is identified in order to perform the accuracy improvement of analysis. We identify the geological boundary after the estimation of blasting force. With this identification procedure, the numerical analysis results which almost correspond with the observation data were provided.

Keywords: Parameter identification, finite element method, average acceleration method, first order adjoint equation method, weighted gradient method, geological boundary, navier equation, optimal control theory.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1332318

References:

[1] B.M.Chaparro, S.Thuillier, L.F.Menezes, P.Y.Manach and J.V.Fernandes : "Material parameters identification: Gradient-based, genetic and hybrid optimization algorithms", J. Comput. Mat. Sci., Vol.44, No.2, pp.339-346, 2008.
[2] F.X.Le Dimet, I.M.Navon and D.N.Daescu : "Second Order Information in Data Assimilation", Month. Weath. Rev., Vol.130:3, pp.629-648, 2002.
[3] F.Yoshida, M.Urabe, V.V.Toropov : "Identification of Material Parameters in Constitutive Model for Sheet Metals from Cyclic Bending Tests", Int. J. Mech. Sci., Vol.40, pp.237-249, 1998.
[4] G.Swoboda, Y.Ichikawa, Q.Dong and M.Zaki : "Back Analysis of Large Geotechnical Models", Int. J. Numer. Analy. Meth. Geomech., John Wiley and Sons Ltd., Vol.23, No.13, pp.1455-1472, 1999.
[5] G.Wu : "Dynamic response analysis of saturated granular soils to blast loads using a single phase model", AGRA Earth and Environmental Ltd., Burnaby, 1995.
[6] K.Kojima, T.Kodama and N.Kaneko : "Parameter Identification and Control of Ground Temperature", Int. J. Comput. Fl. Dyn., Vol.7, pp.193-200, 1996.
[7] M.Kawahara, K.I.Sasaki and Y.Sano : "Parameter Identification and Optimal Control of Ground Temperature", Int. J. Numer. Meth. Fl. Vol.20, pp.789-801, 1995.
[8] M.Piasecki and N.D.Katopodes : "Dispersion Parameter Estimation in Free Surface Flow Region", J. Hyd. Engrg., ASCE, Vol.125, No.7, pp.714-724, 1999.
[9] M.Sezaki, ¨O.Aydan, Y.Ichikawa and T.Kawamoto : "Properties for Initial Design of Tunnels by NATM Using Rock Mass Data-Base System", Japan Society of Civil Engineers Magazine, No.421/VI-13, pp.125-134, 1990.
[10] N.Koizumi and M.Kawahara : "Parameter Identification Method for Determination of Elastic Modulus of Rock Based on Adjoint Equation and Blasting Wave Measurements", Int. J. Numer. Analy. Meth. Geomech., 2008. (In press)
[11] R.Mahnken, E.Stein : "Parameter Identification for Viscoplastic Models Based on Analytical Derivatives of a Least-Squares Functional and Stability Investigations", Int. J. Plasticity, Vol.12, No.4, pp.451-479, 1996.
[12] R.Mahnken, E.Kuhl : "Parameter Identification of Gradient Enhanced Damage Models with the Finite Element Method", Eur. J. Mech., A/Solids, Vol.18, pp.819-835, 1999.
[13] T.Asai and M.Kawahara : "Parameter Identification of Soil Seepage Interraction Problem Using Finite Element Method", Bulletin of the Faculty of Science and Engineering, Chou Univ. Vol.41, pp.15-31, 1998.
[14] T.Ohkami and G.Swoboda : "Parameter Identification of Viscoelastic Materials", J. Comput. Geotech., Elsevier Science Ltd., Vol.24, No.4, pp.279-295, 1999.
[15] U.S.Army Corps of Engineers : "Systematic Drilling and Blasting for Surface Excavations", Engineering Manuals, EM 1110-2-3800, Office of the chief, U. S. Army Corps of Engineers, Washington, D.C., 1972.
[16] Y.Huang and Z.Liu : "Time Domain Parameters Identification of Foundation-structure Interaction System", J. Appl. Math. Mech., Vol.26, No.7, pp.855-864, 2005.
[17] Y.Takahashi and M.Kawahara : "Identification of Reynolds Number Using Automatic Differentiation", Inverse Problems, Design and Optimization Symposium, Rio de Janeiro, Brazil, 2004.
[18] Z.Xiang, G.Swoboda and Z.Cen : "Optimal Layout of Displacement Measurements for Parameter Identification Process in Geomechanics", Int. J. Geomech., Vol.3, No.2, pp.205-216, 2003.