Surface wave analysis for heterogeneous geological formations in geothermal fields: Effect of wave propagation direction

We applied advanced surface wave analysis for multichannel and multishot seismic data to estimate S-wave velocity structure with high lateral resolution. Although a horizontally layered structure is typically assumed in surface wave analysis, this assumption might be violated in environments such as geothermal fields because of their heterogeneous geological formation. The lateral variation of phase velocity can be effectively estimated with common midpoint (CMP) cross-correlation (CMPCC) analysis. In this study, we introduced two additional approaches into the CMPCC analysis workflow, to further improve lateral resolution of phase velocity estimates. One approach is window-controlled CMPCC analysis, which applies a spatial window for the CMP gathers while maintaining the accuracy of the phase velocity estimates. In this analysis, we found that it is difficult to improve phase velocity estimates at lower frequencies, as a wider spatial window must be kept to maintain the accuracy of the dispersion curve. Therefore, we introduced another approach in which we consider the direction of surface wave propagation from sources to receivers. Selection of cross-correlations based on source-CMP direction provides a significant improvement in dispersion curve resolution in the presence of lateral velocity variations, for a wide frequency range. We applied direction-controlled CMPCC analysis to seismic data acquired in the Yamagawa geothermal field, Kyushu Island, southwest Japan, and obtained dispersion curves in the heterogeneous geological setting. We then obtained S-wave velocity profiles by applying genetic algorithm inversion to the dispersion curves. The S-wave velocity profiles in the geothermal field resolve shallow and local heterogeneous structures (e.g. volcanic ash, pumice and igneous intrusions) that cannot be identified on reflection seismic profiles.