###### Abstract  
In order to clarify the water invasion law and residual gas distribution characteristics in edge and bottom water gas reservoirs with multi-well development, a large-scale three-dimensional physical simulation model was developed and a physical simulation experiment method for the water invasion law of multi-well development in sandstone gas reservoirs with strong water drives was established. Water invasion physical simulation experiments of multi-well development under the conditions of different water body multiples and production systems were conducted. The results show the following: (1) Gas wells near fractures and high-permeability zones experience the earliest water breakthrough. The larger the water body multiple, the faster the rate of water invasion, the earlier the water breakthrough time of gas wells, the more severe the degree of water invasion in gas reservoirs, and the lower the ultimate recovery rate. (2) Shutting in low-position gas wells immediately after water breakthrough reduces the overall water production of the gas reservoir and extends the overall water-free gas production period. However, the ultimate recovery rate is lower than when the wells are not shut in. (3) The residual gas in the fracture model is mainly distributed around the fracture and the edge of the gas reservoir, with the ultimate recovery rate ranging from 38.5% to 58.2%. The residual gas in the fracture-high-permeability zone model is mainly distributed around the fracture-high-permeability zone and the edge of the gas reservoir, with the ultimate recovery rate ranging from 28.32% to 41.8%. The experimental results have important guiding significance for the economical and efficient development of similar gas reservoirs. multi-well development; large-scale three-dimensional physical simulation model; water invasion law; residual gas distribution characteristics 26 July 202428 August 202428 August 20243 September 20249 September 20  
declines rapidly [9]. Meanwhile, much gas is trapped during water invasion by the cut-off, bypass flow, blind end, and water lock, reducing the ultimate recovery rate of gas reservoirs [10; 11; 12; 13; 14]. Moreover, with the aggravation of water invasion, the wellbore fluid increases, the abandonment pressure rises, and the gas-water simultaneous production, and even the shut-in due to waterflooding, production cost, and development difficulty significantly increase [15; 16; 17; 18; 19]. In recent years, many scholars have conducted numerous physical simulation experiments to study the law of water invasion in gas reservoirs. [PERSON] et al. [20] considered the influences of different water body multiples, gas reservoir pressure, and single-well production allocation, carrying out physical simulation experiments utilizing three kinds of full-diameter long-core combination models designed to study the water invasion mechanism of fractured bottom water gas reservoirs. [PERSON] et al. [21] established a physical simulation experimental model for multi-layer commingled production by combining typical reservoir cores in parallel according to the similarity criterion. Then, the model was used to study the influence of inter-layer heterogeneity, drawdown pressure, water saturation, and water invasion on the gas supply capacity. [PERSON] et al. [22] conducted experiments using a visualization sand filling tube to study the effects of the dissolved gas on the law of the gas-water contact changes in the water-soluble gas release process. [PERSON] et al. [23] developed physical simulation experiments of water invasion in fractured gas reservoirs using the fracture core samples generated by artificial fracturing to analyze the effects of water body multiples, gas production rates, matrix permeability, and fracture apertures on gas production. [PERSON] et al.