Ground-hole-tunnel combined fine detection method and system

A ground-borehole-tunnel combined fine detection system, comprising:

• a pilot receiving station array for rock-breaking seismic source, a tunnel receiving station array, a ground surface receiving station array, a far-borehole condition geological drilling receiving station array, a near-borehole condition geological drilling receiving station array, a time synchronization system, and a seismic wave data processing instrument system;
• wherein,
• the pilot receiving station array for rock-breaking seismic source is positioned behind a cutterhead of a tunnel boring machine (TBM), the tunnel receiving station array is installed on a middle portion of a TBM body or surrounding rock of tunnel sidewalls, the ground surface receiving station array is positioned on the ground surface in front of the tunneling face, and the far-borehole and near-borehole condition geological drilling receiving station arrays are positioned inside geological drilling boreholes with different distances in front of the tunneling face;
• the time synchronization system is configured to perform precise synchronization on each receiving station array;
• and
• the seismic wave data processing instrument system is configured to receive and store the observation data of each receiving station array and perform fast processing;
• wherein:
• a seismic record integrated processing method comprises:
• preprocessing received signals;
• performing near-borehole condition fine extraction;
• performing far-borehole condition wavefield optimization;
• performing face and body wave separation on data acquired by the ground surface receiving station array;
• performing cross-correlation and deconvolution processing on seismic source signals and processed received signals;
• importing a coordinate of an observation system, and performing automatic first break picking;
• performing spectral analysis and band-pass filtering;
• performing intra-trace equalization and inter-trace equalization;
• suppressing ineffective reflection waves, remaining effective reflection waves, and performing P-wave and S-wave separation;
• performing extraction of surface wave frequency dispersion curves, and frequency dispersion of an energy diagram;
• performing transmitted wave-surface wave joint inversion on data acquired in an optimized far-borehole condition geological drilling and surface wave data frequency dispersion curves, and acquiring a velocity model in front of the working face of the TBM using a joint inversion method;
• and
• acquiring the seismic section in front of the working face of the TBM using the velocity model acquired through joint inversion and adopting reverse time migration imaging;
• wherein, an objective function of the transmitted wave-surface wave joint inversion is as follows:
• Φ
• =
• ∑
• 
• W
• d
• [
• d
• o
• ⁢
• b
• ⁢
• s
• -
• G
• ⁡
• (
• m
• )
• ]
• 
• 2
• +
• λ
• ⁢
• 
• W
• c
• (
• m
• -
• m
• 0
• )
• 
• 2
• +
• β
• ⁢
• 
• τ
• ⁡
• (
• m
• )
• 
• ;
• wherein, a first item of a right formula is a data fitting item, dobs is an acquired transmitted wave data and extracted surface wave frequency dispersion curve, G(m) is a transmitted wave data and forward frequency dispersion curve obtained through inversion model forward, a second item of the right formula is a model fitting item, m is a model obtained through current inversion, m0 is an initial model, a third item of the right formula is a cross-gradient item, τ(m)=∇mp(x,y,z)×∇ms(x,y,z), my represents a P-wave velocity, ms represents a S-wave velocity, and Wd, Wc, λ, and β are respectively weights.