**6. Conclusions**

The iterative virtual events method, when used to suppress internal multiples, involves first using the traditional virtual events method to build the initial multiple models, and then iterating the initial multiple models to complete the suppression of internal multiples. This approach not only simplifies the process of multiple suppression and improves computational efficiency and accuracy, but also obtains more abundant multiple information. From the model test, field data processing, and profile and spectrum analysis before and after processing, we can conclude the following:


These results indicate that this method can effectively suppress the internal multiples in seismic data and can be widely used in the internal multiples suppression of field data.

**Author Contributions:** Conceptualization, M.X. and W.W.; methodology, J.X. and W.L.; software, J.S.; validation, W.L., J.S. and B.J.; formal analysis, Y.Z. and W.L.; investigation, T.Z. and Y.Z.; resources, Y.W. and B.J.; data curation, J.X., Y.W. and T.Z.; writing—original draft preparation, M.X. and Y.W.; writing—review and editing, W.W. and J.X.; visualization, B.J. and M.X.; supervision, W.W. and W.L.; project administration, Y.Z.; funding acquisition, T.Z. and W.W. All authors have read and agreed to the published version of the manuscript.

**Funding:** We thank the prospective basic projects of CNPC under grant no. 2021DJ-0304 and the National Natural Science Foundation of China under grant nos. 41974116 and 4227040468 for supporting this work.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.
