The Tanlu fault zone is the most active fault zone in eastern China, exhibiting multi-phase, complex, and segmented activity characteristics. Based on its geometric structure and activity, the Tanlu fault zone can be divided into four sections from south to north: the Jiashan-Guangji, Weifang-Jiashan, Xialiaohe-Laizhouwan, and Hegang-Tieling sections. Among these, the Weifang-Jiashan section displays the strongest tectonic activity, within which the F₅ fault （Anqiu-Juxian fault） is identified as the most recent Holocene active fault. Studying its shallow tectonic characteristics is critical for understanding fault evolution and assessing seismic hazards.

Domestic and international scholars have reached a consensus on the developmental scale, strike, and activity characteristics of the F₅ fault. However, rapid and effective methods for detecting fault zone width and internal medium structure remain lacking. While core drilling offers high precision, its high cost and lengthy duration limit widespread application. Reflection wave exploration can locate fault points by analyzing co-phase axis deformation and diffracted wave development, but it fails to resolve detailed internal medium structures. Advances in observation technology and data processing now make high-resolution surface wave exploration a promising solution.

Rayleigh wave exploration, an emerging geophysical method, enables rapid and cost-effective shear wave velocity measurements, making it a research hotspot. Rayleigh wave methods are categorized into active-source and passive-source techniques based on seismic excitation modes. The active-source method yields high-frequency signals but has limited depth resolution, whereas the passive-source method provides low-frequency signals for deeper exploration. Combining both approaches enhances effective exploration depth while maintaining shallow-structure accuracy.

This study employs a joint detection method for active and passive surface waves, acquiring data using identical array configurations for both sources. Dispersion spectra are superimposed to improve dispersion curve quality and expand the observable frequency range. A genetic algorithm inverts the 2D S-wave velocity structure.

Results reveal significant lateral heterogeneity in the crustal medium, with a segmented S-wave velocity structure. Fault points are identified at distances 110 m （F_5-1-1）, 160 m （F_5-1）, and 640 m （F_5-2）. The F_5-1 and F_5-2 faults form opposing boundary faults （dip angles 70° and 65°, respectively）, creating a graben structure with a 330 m-wide fracture zone. The western boundary is not a single fault but a composite zone comprising F_5-1 and its secondary fault F_5-1-1. The fracture zone exhibits overall low S-wave velocities but contains localized high-velocity anomalies, likely reflecting F₅ fault activity.

Comparisons with shallow seismic exploration and borehole profiles confirm that the joint method accurately locates main fault surfaces, while also providing clearer delineation of fracture zone width and internal material structure. This approach is particularly effective for shallow faults or high-interference environments like urban areas.