The ultradeep Ordovician limestone fracture-vug reservoirs in the Fuman oilfield on the south bank of the Tahe River in the Tarim Basin is the core target area for the production capacity construction of the oilfield.These reservoirs in the study area contain massive dissolution vugs formed by the formation fracturing due to the strike-slip faulting and thus are highly heterogeneous.The low-frequency models based on conventional wave impedance inversion are built using data on horizons,faults,and logs,and thus they cannot characterize the fault-controlled heterogeneity of the reservoirs.For this reason,this study proposed an iterative inversion method constrained by fault fractured zone facies to predict reservoirs.The technical process of this method is as follows.First,the original seismic data were interpreted to obtain the attribute volumes that can reflect the characteristics of the fault fractured zones of the reservoirs.Then,the attribute volumes were proportionally fused with the initial low-frequency model and the spatial profile data of high-quality reservoirs obtained from the conventional inversion based on wave impedance.As a result,a new heterogeneous low-frequency model was formed.Using this model,multiple rounds of iterative inversions were conducted.The inversion results can truly describe the characteristics of the fault fractured zones in the heterogeneous reservoirs.The prediction results of 10 wells had coincidence rates of up to 92.86%.As indicated by the application,the method proposed in this study can improve the reservoir prediction performance and the prediction precision of fault-controlled reservoirs,thus effectively supporting the work in the study area.