Designing effective coverage routes for unmanned surface vehicles (USVs) is crucial to improve the efficiency of offshore bathymetric surveys. However, existing coverage planning methods for practical use are limited, primarily due to the large-scale surveying areas and intricate region geometries caused by coastal features. This study aims to address these challenges by introducing a coverage path planning framework for USV-assisted bathymetric mapping, specifically aimed at the joint optimization of paths to cover numerous complex regions. Initially, we conceptualize the large-scale bathymetric survey mission as an integer programming model. The model uses four distinct decision variables to meticulously formulate length calculations, inter-regional connections, entry and exit point selections, and line sweep direction. Then, a novel hierarchical algorithm is devised to solve the problem. The method first incorporates a bisection-based convex decomposition method to achieve optimal partitioning of complex regions. Additionally, a hierarchical heuristic optimization algorithm that seamlessly integrates the optimization of all influencing factors is designed, which includes order generation, candidate pattern finding, tour finding, and final optimization. The reliability of the framework is validated through semi-physical simulations and lake trials using a real USV. Through comparative studies, our model demonstrates clear advantages in computational efficiency and optimization capability compared to state-of-the-arts, with its superiority becoming more pronounced as the problem scale increases. The results from lake trials further affirm the efficient and reliable performance of our model in practical bathymetric survey tasks.