TY - JOUR
T1 - Effective seismic force retrieval from surface measurement for SH-wave reconstruction
AU - Guidio, Bruno
AU - Goh, Heedong
AU - Jeong, Chanseok
N1 - Funding Information:
This material is based upon work supported by the National Science Foundation, USA under Awards CMMI-2044887 and CMMI-2053694 . Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation. The authors are also grateful for the support by the Faculty Research and Creative Endeavors (FRCE) Research Grant- 48058 at Central Michigan University.
Publisher Copyright:
© 2022 Elsevier Ltd
PY - 2023/2
Y1 - 2023/2
N2 - We present a new method to obtain dynamic body force at virtual interfaces to reconstruct shear wave motions induced by a source outside a truncated computational domain. Specifically, a partial differential equation (PDE)-constrained optimization method is used to minimize the misfit between measured motions at a limited number of sensors on the ground surface and their counterparts reconstructed from optimized forces. Numerical results show that the optimized forces accurately reconstruct the targeted ground motions in the surface and the interior of the domain. The proposed optimization framework yields a particular force vector among other valid solutions allowed by the domain reduction method (DRM). Per this optimized or inverted force vector, the reconstructed wave field is identical to its reference counterpart in the domain of interest but may differ in the exterior domain from the reference one. However, we remark that the inverted solution is valid and introduce a simple post-process that can modify the solution to achieve an alternative force vector corresponding to the reference wave field. We also study the desired sensor spacing to accurately reconstruct the wave responses for a given dominant frequency of interest. We remark that the presented method is omnidirectionally applicable in terms of the incident angle of an incoming wave and is effective for any given material heterogeneity and geometry of layering of a reduced domain. The presented inversion method requires information on the wave speeds and dimensions of only a reduced domain. Namely, it does not need any information on the geophysical profile of an enlarged domain or a seismic source profile outside a reduced domain. Thus, the computational cost of the method is compact even though it leads to the high-fidelity reconstruction of wave response in the reduced domain, allowing for studying and predicting ground and structural responses using real seismic measurements.
AB - We present a new method to obtain dynamic body force at virtual interfaces to reconstruct shear wave motions induced by a source outside a truncated computational domain. Specifically, a partial differential equation (PDE)-constrained optimization method is used to minimize the misfit between measured motions at a limited number of sensors on the ground surface and their counterparts reconstructed from optimized forces. Numerical results show that the optimized forces accurately reconstruct the targeted ground motions in the surface and the interior of the domain. The proposed optimization framework yields a particular force vector among other valid solutions allowed by the domain reduction method (DRM). Per this optimized or inverted force vector, the reconstructed wave field is identical to its reference counterpart in the domain of interest but may differ in the exterior domain from the reference one. However, we remark that the inverted solution is valid and introduce a simple post-process that can modify the solution to achieve an alternative force vector corresponding to the reference wave field. We also study the desired sensor spacing to accurately reconstruct the wave responses for a given dominant frequency of interest. We remark that the presented method is omnidirectionally applicable in terms of the incident angle of an incoming wave and is effective for any given material heterogeneity and geometry of layering of a reduced domain. The presented inversion method requires information on the wave speeds and dimensions of only a reduced domain. Namely, it does not need any information on the geophysical profile of an enlarged domain or a seismic source profile outside a reduced domain. Thus, the computational cost of the method is compact even though it leads to the high-fidelity reconstruction of wave response in the reduced domain, allowing for studying and predicting ground and structural responses using real seismic measurements.
KW - Discretize-then-optimize (DTO) approach
KW - Domain reduction method (DRM)
KW - Effective seismic force vector
KW - Full-waveform inversion
KW - Passive-seismic inversion
KW - Reconstruction of seismic responses
UR - http://www.scopus.com/inward/record.url?scp=85145614672&partnerID=8YFLogxK
U2 - 10.1016/j.soildyn.2022.107682
DO - 10.1016/j.soildyn.2022.107682
M3 - Article
AN - SCOPUS:85145614672
VL - 165
JO - Soil Dynamics and Earthquake Engineering
JF - Soil Dynamics and Earthquake Engineering
SN - 0267-7261
M1 - 107682
ER -