The microtremor exploration technique has been widely used in the geological exploration of a complex urban environment due to its strong spatial adaptability and anti-mechanical interference ability.The forward modeling of microtremor characteristic curves is all based on a horizontally layered medium model,while actual strata are typical inhomogeneous geological models.The inconsistency leads to the unsatisfactory fitting degree of characteristic curves.Therefore,it is necessary to conduct forward modeling research based on the condition of inhomogeneous media.Since the geological model based on the assumption of horizontally layered media cannot meet the requirements for high-precision forward modeling of microtremor characteristic curves,forward modeling needs to be conducted using an inhomogeneous model.Based on the dynamic optimization method,this study proposed the forward modeling technique of microtremor characteristic curves to realize the forward modeling of high-precision characteristic curves of inhomogeneous media.In this technique,the aperture range and the optimization coefficients varying with frequency were calculated based on the theoretical characteristic curves of equivalent horizontally layered media.Taking spectral ratio curves as an example,the forward modeling technique of microtremor characteristic curves using the dynamic optimization method can obtain the shear wave distribution of inhomogeneous media through inversion,achieving higher lateral resolution and exploration precision.
Aki K. Space and time spectra of stationary stochastic waves, with special reference to microtremors[J]. Bulletin of Earthquake Research Institution, 1957, 35: 415-456.
[2]
Toksöz M N, Lacoss R T. Microseisms: Mode structure and sources[J]. Science, 1968:159.
Ran W Y, Wang Z D. The long-wave microtremors method and its advances[J]. Geophysical and Geochemical Exploration, 1994, 18(1): 28-34.
[4]
Haghshenas E, Bard P Y, Theodulidis N, et al. Empirical evaluation of microtremor H/V spectral ratio[J]. Bulletin of Earthquake Engineering, 2008, 6(1): 75-108.
doi: 10.1007/s10518-007-9058-x
[5]
Arai H, Tokimatsu K. S-wave velocity profiling by inversion of microtremor H/V spectrum[J]. Bulletin of the Seismological Society of America, 2004, 94(1): 53-63.
doi: 10.1785/0120030028
Liu T H, Liu Tie, Cheng G H, et al. Research progress of geophysical exploration technology in complex urban environment[J]. Chinese Journal of Engineering Geophysics, 2020, 17(6) :711-720.
[7]
Rayleigh J W S. On waves propagated along the plane surface of an elastic solid[J]. Proceedings of the London Mathematical Society, 1885, 17: 4-11.
[8]
Haskell N A. The dispersion of surface waves on multilayered[J]. Bulletin of the seismological Society America, 1953, 43(1): 17-34.
doi: 10.1785/BSSA0430010017
[9]
Aki K, Richards P G. Quantitative Seismology: Theory and methods[M]. San Francisco: CA, 1980.
Xu P F, Li C J, Ling S Q, et al. Mapping collapsed columns in coal mines utilizing microtremor survey methods[J]. Chinese Journal of Geophysics, 2009, 52(7):1923-1930.
Wang W J, Liu L B, Chen Q F, et al. Applications of microtremor H/V spectral ratio and array techniques in assessing the site effect and near surface velocity structure[J]. Chinese Journal of Geophysics, 2009, 52(6):1515-1525.
Liu T H. Research and application of geological exploration methods in urban drilling blind area[J]. Chinese Journal of Engineering Geophysics, 2019, 16(5):572-579.
Zhang M H, Wu Z B, Ma L X, et al. Research progress of passive source detection technology based on short-period dense seismic array[J]. Progress in Geophysics, 2020, 35(2):495-511.
Liu T H. Analysis and application of geological exploration methods in Urban drilling blind area[J]. Journal of Railway Engineering Society, 2019, 36(10):88-93.
Liu W, Huang T, Wang T Y, et al. The application of integrated geophysical prospecting methods to the exploration of urban buried fault[J]. Geophysical and Geochemical Exploration, 2021, 45(4): 1077-1087.
Dong Y, Li G H, Gao P J, et al. The application of fretting exploration technology in the exploration of middle and deep clean energy[J]. Geophysical and Geochemical Exploration, 2020, 44(6):1345-1351.
[17]
Knopoff L. A matrix method for elastic wave problems[J]. Bulletin of the Seismological Society of America, 1964, 54(1):431-438.
doi: 10.1785/BSSA0540010431
[18]
Abo-Zena A. Dispersion function computations for unlimited frequency values[J]. Geophysical Journal International, 1979, 58(1):91-105.
doi: 10.1111/j.1365-246X.1979.tb01011.x
[19]
Chen X. A systematic and efficient method of computing normal modes for multilayered half-space[J]. Geophysical Journal International, 1993, 115(2):391-409.
doi: 10.1111/j.1365-246X.1993.tb01194.x
Fan Y H, Liu J Q, Xiao B X. Fast vector-transfer algorithm for computation of Rayleigh wave dispersion curves[J]. Journal of Hunan University:Natural Sciences Edition, 2002(5):25-30.
[21]
Arai H, Tokimatsu K. S-wave velocity profiling by inversion of microtremor H/V spectrum[J]. Bulletin of the Seismological Society of America, 2004, 94(1): 53-63.
doi: 10.1785/0120030028
[22]
张立. 层状介质中瑞利面波波场特征分析和反演方法研究[D]. 成都: 西南交通大学, 2009.
[22]
Zhang L. An Approach to the wave-field characteristic analysis and inversion algorithm of Rayleigh surface wave for layered medium[D]. Chengdu: Southwest Jiaotong University, 2009.