Abstract: To further develop structural reliability analysis methods for high-pile wharf and shore slope systems, this paper, based on the KL-Copula cross-correlated random field simulation method, comprehensively considers the influence of spatial variability of slope soil parameters, its related structures, and dimensions. It proposes a numerical implementation procedure for shore slope-high-pile wharf systems cross-correlated random field simulation and reliability solution, formulating the high-pile wharf-shore slope system cross-correlation random field simulation and reliability analysis plug-in solver KLCPRP based on the PLAXIS computing program platform. Building on this, non-intrusive stochastic analysis methods are further employed for related example analysis to explore the impact of different slope strength parameter Copula function structures on cross-correlated random field simulation and high-pile wharf structural reliability solutions. Computational analysis indicates: there exist certain differences in the linear correlation coefficients and rank correlation coefficients of simulation samples corresponded to different Copula function structures; for limited experimental data scenarios, the Gaussian Copula function structure tends to yield conservative random field computation results. The KL-Copula cross-correlated random field numerical implementation procedure and the corresponding plug-in solver KLCPRP presented in this paper possess favorable accuracy and effectiveness, offering a theoretical reference for the related engineering design and safety evaluation of high-pile wharf and shore slope interaction systems.