A receiver function (RF) is the response of the Earth's structure below a seismometer to an incident teleseismic wave and consists of a series of P-to-S (Ps) or S-to-P (Sp) converted waves generated at structural interfaces (Langston, 1979), mainly those generated by the velocity discontinuities in the crust and upper mantle below the seismometer. Crustal thickness (H) and P-wave and S-wave velocity ratio (κ) are important parameters reflecting the crustal structure and internal material composition and can provide an important basis for regional tectonics and dynamics. Zhu and Kanamori (2000) noted that when the Moho converted phase Ps and crustal multiple waves are considered at the same time, H and κ can be estimated under the assumption of the average crustal P wave velocity. This is currently the most commonly used H-κ method. The H-κ algorithm stacks the amplitudes of RFs at predicted arrival times for Ps and crustal multiples (PpPs and PpSs + PsPs; here referred to as M1 and M2, respectively, for convenience, as in J. Li et al., 2019) for different values of H and κ in a grid search.The assumptions of the H-κ method are very simple, which leads to its limitations. It assumes that the crustal medium is isotropic, the Moho surface is horizontal, and the P-wave velocity of the crust is fixed. However, the above assumptions cannot be satisfied under complex geological settings, so the results of the H-κ method may be inaccurate. The Chinese continent is known to have highly heterogeneous crustal structure and extremely variable Moho depth, as revealed for example, by seismic tomography (