Low-lying coastal cities are vulnerable to compound floods caused by many factors including river flows, tides and local rainfall. Many previous studies focus on the impacts of rainfall and tidal levels (two driving factors) on estuaries or regions near the main single river, while research about the three influencing factors on the floods for complex urban river networks remain relatively scarce. This paper proposes a framework to simulate the coastal compound flooding process, with a series of scenarios designed to explore the main flooding drivers and their worst possible combination. The approach is used to quantify the induced flooding consequences, in which a two-dimensional model is employed to conduct a case study for a city with complex river networks. The Qianshan River basin in Zhuhai City, located in the Guangdong-Hong Kong-Macao Greater Bay Area (GBA) of China is used as the case study. Application results show that: (i) the occurrence and properties of flooding in this coastal area are jointly affected by upstream river flows, downstream tidal levels and local rainfall intensity with different contribution levels, (ii) compared to tidal levels, local rainfall intensity has a larger impact on the flooding levels of the study region, and (iii) the worst scenario where a 50-year return period upstream river flow co-occurring with high tide and a 100-year return period rainfall can produce a flooding area up to 50 km² for the study region, with an average of 1.3 m flooded depth. The proposed framework with its findings offers an in-depth understanding and quantification of coastal urban flood characteristics within the GBA of China.