Multi-wave imaging method, system and equipment for suppressing stretching distortion and medium

By suppressing the large-angle outgoing signals of the forward and reverse propagation wave fields during the initial propagation and final imaging stages, the problem of stretching distortion in multiple wave imaging is solved, and efficient multiple wave imaging effect is achieved.

CN121784823APending Publication Date: 2026-04-03CHINA NAT OFFSHORE OIL CORP +1
View PDF 0 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2025-12-08
Publication Date
2026-04-03

AI Technical Summary

Technical Problem

In marine seismic data processing, stretching distortion occurs during multiple wave imaging, leading to reduced imaging resolution. Existing methods are computationally intensive and difficult to apply in industrial production.

Method used

In the initial propagation and final imaging stages, large-angle outgoing signals from the forward and reverse propagation wave fields are suppressed respectively. Through Fourier transform and outgoing angle calculation, large-angle signals are suppressed using a threshold to generate multiple wave imaging results.

Benefits of technology

With a small increase in computation, stretching distortion is effectively suppressed, improving the resolution and imaging quality of multiple wave imaging.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN121784823A_ABST
    Figure CN121784823A_ABST
Patent Text Reader

Abstract

The invention relates to a multiple wave imaging method, system and device for suppressing stretching distortion and a medium. The method comprises the following steps: acquiring a speed model and an initial wave field from seismic data; performing large-angle emergent signal suppression on the initial forward transmission wave field and the back transmission wave field; continuously propagating the suppressed wave field according to a given speed model to obtain a global forward propagation wave field and a reverse propagation wave field; performing large-angle emergent signal suppression in the imaging stage on the global forward transmission wave field and the back transmission wave field respectively; and based on the suppressed wave field, according to a preset imaging condition, generating a multiple imaging result. The method can be widely applied to the field of seismic data processing.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention belongs to the field of marine seismic data processing, specifically relating to a method, system, device, and medium for suppressing stretching distortion in multiple wave imaging. Background Technology

[0002] Multiples are one of the most common and critical types of noise in marine seismic data processing. Improper suppression can easily damage the effective signal or leave false axes in the imaging profile. In recent years, imaging techniques that utilize multiples as effective signals have been continuously developed and have achieved practical application results.

[0003] Multiple wave imaging treats all seismic records from receivers as virtual sources for secondary propagation, significantly improving illumination. Compared to primary wave imaging, multiple waves have longer propagation paths and smaller emission angles at the same offset, resulting in higher vertical resolution. However, during migration imaging, seismic waves suffer from stretching distortion; signals with large emission angles typically appear as low-frequency signals in the imaging gather, reducing imaging resolution. For primary wave imaging, the source location of each shot gather is fixed, and under most geological conditions, the emission angle is positively correlated with the offset, so stretching distortion can be suppressed by removing large offset signals from the offset gather. However, for multiple wave imaging, each receiver is a virtual source, making it impossible to calculate the offset gather.

[0004] Suppression of stretching distortion in multiple-wave imaging based on gathers requires calculating the imaging angle gather. The calculation of the angle gather involves decomposing the forward propagating wavefield according to the incident angle and the reverse propagating wavefield according to the reflection angle, based on conventional wavefield extrapolation. The difference between the incident and reflection angles is then divided by 2 to obtain the true reflection angle of the subsurface medium. The imaging results of the forward and reverse propagating wavefields with the same true reflection angle are calculated and summed to obtain the imaging angle gather. It is evident that the calculation of the angle gather involves multiple loops and is computationally intensive, rarely used in practical industrial production. Summary of the Invention

[0005] To address the aforementioned problems, the present invention aims to provide a method, system, device, and medium for efficient suppression of stretching distortion in multiple-wave imaging. By suppressing large-angle signals in both the initial propagation and final imaging stages, stretching distortion can be effectively suppressed with only a small increase in computational load.

[0006] To achieve the above objectives, the present invention adopts the following technical solution: In a first aspect, the present invention provides a multiple-wave imaging method for suppressing stretching distortion, comprising: Obtain velocity models and initial wavefields from seismic data; Suppress the large-angle outgoing signals of the initial forward propagation wavefield and the reverse propagation wavefield respectively; Based on the given velocity model, the initial forward propagation wave field and reverse propagation wave field after the large-angle outgoing signal is suppressed are continued to propagate to obtain the global forward propagation wave field and reverse propagation wave field. Suppress large-angle outgoing signals during the imaging stage of both the global forward propagation wavefield and the reverse propagation wavefield. Based on the large-angle outgoing signal suppression of the wave field during the imaging stage, multiple wave imaging results are generated according to preset imaging conditions.

[0007] Furthermore, the suppression of the large-angle emitted signals in the initial forward propagation wavefield and the reverse propagation wavefield respectively includes the following steps: ① For the initial forward propagation wave field Perform a time-domain Fourier transform to obtain the initial forward propagating wavefield in the frequency domain. ,in, Let be the coordinate point. For frequency; ② Initial forward propagating wave field in the frequency domain Perform a two-dimensional Fourier transform to obtain ,in, , They are respectively Wave number in direction; ③ Based on the surface velocity of the medium Calculate each frequency at each point on the sea surface. Corresponding wave number and the angle of departure ; ④ Based on a given threshold and the angle of departure The surface-emitted signal is suppressed to obtain the initial forward propagation wavefield after suppressing the large-angle-emitted signal. ; ⑤ The initial forward propagation wave field after suppressing the large-angle emitted signal Performing an inverse Fourier transform yields the initial forward propagation wavefield in the spatial domain that suppresses large-angle outgoing signals. ; ⑥ Using the same method as steps ① to ⑤, the initial reverse propagation wave field is... Suppressing large-angle outgoing signals yields the initial backpropagation wavefield in the spatial domain for suppressing large-angle outgoing signals. .

[0008] Furthermore, the wavenumber The calculation formula is:

[0009] Angle of departure The calculation formula is: .

[0010] Furthermore, the threshold-based and the angle of departure Suppressing the emitted signal from the surface layer specifically refers to suppressing signals exceeding a given threshold. The surface-emitted signal is suppressed to obtain the initial forward propagation wavefield after suppressing the large-angle emitted signal. The calculation formula is: .

[0011] Furthermore, the large-angle outgoing signal suppression during the imaging stage of the global forward and reverse propagation wavefields includes: ① Global forward propagation wave field in the spatial domain Perform a two-dimensional Fourier transform to obtain ; ②Based on the velocity of the medium Calculate the angle of departure at each point on the sea surface. ; ③ For emission angles greater than a given threshold The signal is suppressed to obtain the global forward propagation wavefield after suppressing the large-angle outgoing signal. ; ④ Global forward propagation wave field after suppressing large-angle emitted signals Performing an inverse Fourier transform yields the global forward propagation wavefield in the spatial domain after suppressing large-angle outgoing signals. ; ⑤ Using the same method as steps ① to ④, suppress the large-angle outgoing signal of the global backpropagation wavefield in the space domain to obtain the global backpropagation wavefield in the space domain after suppressing the large-angle outgoing signal. .

[0012] Furthermore, the emission angle The calculation formula is:

[0013] The global forward propagation wave field after suppressing energy emitted at large angles The calculation formula is: .

[0014] Furthermore, the calculation formula for the multiple wave imaging result is as follows:

[0015] in, The stability parameter is calculated using the following formula:

[0016] in, for . conjugate.

[0017] In a second aspect, the present invention provides a multiple-wave imaging system for suppressing stretching distortion, comprising: The data acquisition module is configured to acquire velocity models and initial wavefields from seismic data; The initial wavefield suppression module is configured to suppress the large-angle outgoing signals of the initial forward propagation wavefield and the reverse propagation wavefield, respectively. The propagation module is configured to continue propagating the initial forward propagation wavefield and reverse propagation wavefield after the large-angle outgoing signal is suppressed according to a given velocity model, so as to obtain the global forward propagation wavefield and reverse propagation wavefield. The imaging stage suppression module is configured to suppress large-angle outgoing signals in the global forward propagation wavefield and reverse propagation wavefield during the imaging stage, respectively. The imaging module is configured to generate multiple wave imaging results based on the suppressed wavefield of the large-angle outgoing signal during the imaging stage, according to preset imaging conditions.

[0018] Thirdly, the present invention provides a computer-readable storage medium for storing one or more programs, said one or more programs including instructions that, when executed by a computing device, cause the computing device to perform any method.

[0019] Fourthly, the present invention provides a computing device comprising: one or more processors and a memory, wherein the memory stores one or more programs and is configured to be executed by the one or more processors, the one or more programs including instructions for performing any method.

[0020] The present invention has the following advantages due to the adoption of the above technical solutions: This invention applies a suppression stretching distortion algorithm to the forward propagation wavefield and the reverse propagation wavefield in the initial propagation and final imaging stages, respectively. It can effectively suppress stretching distortion and achieve high imaging quality with only a small increase in computation.

[0021] Therefore, this invention can be widely applied in the field of seismic data processing. Attached Figure Description

[0022] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Throughout the drawings, the same reference numerals denote the same parts. In the drawings: Figure 1 This is a flowchart of the multiple-wave imaging method for suppressing stretching distortion provided in the embodiments of the present invention; Figure 2 This is the Marmousi velocity model provided in the embodiments of the present invention; Figure 3 (a) ~ Figure 3 (d) is the imaging result of a single shot set provided in an embodiment of the present invention, wherein, Figure 3 (a) indicates that a large-angle resection will not be performed; Figure 3 (b) indicates the removal of energy with an initial positive and negative wave field emission angle greater than 40°; Figure 3 (c) indicates the removal of energy with a global forward and reverse propagation wave field exit angle greater than 40°; Figure 3 (d) indicates that energy with an initial forward and reverse propagation wave field emission angle of 40° or higher and global forward and reverse propagation wave field emission angles of 40° or higher are simultaneously removed. Figure 4 (a) ~ Figure 4 (d) is the image overlay result of all the shot sets (400 shots) provided in the embodiment of the present invention, wherein, Figure 4 (a) indicates that a large-angle resection will not be performed; Figure 4 (b) indicates the removal of energy with an initial positive and negative wave field emission angle greater than 40°; Figure 4 (c) indicates the removal of energy with a global forward and reverse propagation wave field exit angle greater than 40°; Figure 4 (d) indicates that the energy of the initial positive and negative wave field emission angles above 40° and the global positive and negative propagation wave field emission angles above 40° are simultaneously removed. Detailed Implementation

[0023] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the described embodiments of the present invention are within the scope of protection of the present invention.

[0024] It should be noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the scope of exemplary embodiments according to the invention. As used herein, the singular form is intended to include the plural form as well, unless the context clearly indicates otherwise. Furthermore, it should be understood that when the terms "comprising" and / or "including" are used in this specification, they indicate the presence of features, steps, operations, devices, components, and / or combinations thereof.

[0025] In some embodiments of the present invention, a multiple wave imaging method for suppressing stretching distortion is provided. This method suppresses stretching distortion in both the initial propagation and final imaging stages of the seismic wave, effectively suppressing stretching distortion with only a small increase in computational load.

[0026] Correspondingly, in other embodiments of the present invention, a multiple-wave imaging system, device, and medium for suppressing stretching distortion are provided.

[0027] Example 1 like Figure 1 As shown, the present invention provides a method for suppressing stretching distortion in multiple-wave imaging, which includes the following steps: 1) Obtain the velocity model and initial wavefield from seismic data; 2) Suppress the initial forward propagation wavefield and reverse propagation wavefield at large angles respectively; 3) Based on the given velocity model, the initial forward propagation wavefield and reverse propagation wavefield after the large-angle outgoing signal is suppressed continue to propagate to obtain the global forward propagation wavefield and reverse propagation wavefield. 4) Suppress large-angle outgoing signals during the imaging stage of both the global forward and reverse propagation wavefields. 5) Based on the suppressed wavefield in step 4), generate multiple wave imaging results according to the preset imaging conditions.

[0028] Furthermore, in step 2) above, the suppression of the initial forward and reverse propagation wavefields by large-angle outgoing signals includes the following steps: 2.1) For the initial forward propagation wave field Perform a time-domain Fourier transform to obtain the initial forward propagating wavefield in the frequency domain. ,in, Let be the coordinate point. For frequency; 2.2) Initial forward propagating wave field in the frequency domain Perform a two-dimensional Fourier transform to obtain ,in, , They are respectively Wave number in direction.

[0029] 2.3) Based on the surface velocity of the medium Calculate each frequency at each point on the sea surface. Corresponding wave number and the angle of departure ; Among them, wave number The calculation formula is:

[0030] Angle of departure The calculation formula is:

[0031] 2.4) Based on a given threshold and the angle of departure The surface-emitted signal is suppressed to obtain the initial forward propagation wavefield after suppressing the large-angle-emitted signal. ; In this embodiment, based on a given threshold Suppressing the emitted signal from the surface layer specifically refers to suppressing signals exceeding a given threshold. The surface-emitted energy is suppressed to obtain the initial forward propagation wavefield after suppressing the large-angle emitted signal. The calculation formula is:

[0032] 2.5) Initial forward propagation wavefield after suppressing large-angle outgoing signals Performing an inverse Fourier transform yields the initial forward propagation wavefield in the spatial domain that suppresses large-angle outgoing signals. .

[0033] 2.6) Using the same method as steps 2.1) to 2.5), the initial backpropagation wavefield is processed. Suppressing large-angle outgoing signals yields the initial backpropagation wavefield in the spatial domain for suppressing large-angle outgoing signals. .

[0034] Furthermore, in step 3) above, such as Figure 2 As shown, based on the given velocity model, the initial forward propagation wavefield in the spatial domain for suppressing large-angle outgoing signals is... and reverse propagation wave field By continuing the propagation, the corresponding global forward propagation wave field can be obtained. and reverse propagation wave field .

[0035] Furthermore, in step 4) above, the large-angle outgoing signal suppression during the imaging stage of the global forward and reverse propagation wavefield specifically includes the following steps: 4.1) Global forward propagation wave field in the spatial domain Perform a two-dimensional Fourier transform to obtain ; 4.2) Based on the velocity of the medium Calculate the angle of departure at each point on the sea surface. :

[0036] 4.3) For emission angles greater than a given threshold The signal is suppressed to obtain the global forward propagation wavefield after suppressing the large-angle outgoing signal. :

[0037] 4.4) Global forward propagation wave field after suppressing large-angle emitted signals Performing an inverse Fourier transform yields the global forward propagation wavefield in the spatial domain after suppressing large-angle outgoing signals. ; 4.5) Using the same method as steps 4.1) to 4.4), suppress the large-angle outgoing signal of the global backpropagation wavefield in the space domain to obtain the global backpropagation wavefield in the space domain after suppressing the large-angle outgoing signal. .

[0038] Furthermore, in step 5) above, appropriate imaging conditions are selected, such as deconvolution imaging conditions. The calculation formula for the multiple wave imaging results is as follows:

[0039] in, The stability parameter is calculated using the following formula:

[0040] in, for . conjugate.

[0041] like Figure 3 (a) ~ Figure 3 As shown in (d), this is the imaging result of a single shot collection provided in an embodiment of the present invention; as Figure 4 (a) ~ Figure 4 As shown in (d), this is the image superposition result of all the shot sets provided in the implementation of the present invention. It can be seen that the method of the present invention can effectively suppress stretching distortion and greatly improve the image quality.

[0042] Example 2 The above-described embodiment 1 provides a multiple-wave imaging method for suppressing stretching distortion. Correspondingly, this embodiment provides a multiple-wave imaging system for suppressing stretching distortion. The system provided in this embodiment can implement the multiple-wave imaging method for suppressing stretching distortion of embodiment 1. This system can be implemented through software, hardware, or a combination of both. For example, the system may include integrated or separate functional modules or units to perform the corresponding steps in the methods of embodiment 1. Since the system in this embodiment is basically similar to the method embodiment, the description process in this embodiment is relatively simple. Relevant details can be found in the description of embodiment 1. The system embodiment provided in this embodiment is merely illustrative.

[0043] The multiple-wave imaging system for suppressing stretching distortion provided in this embodiment includes: The data acquisition module is configured to acquire velocity models and initial wavefields from seismic data; The initial wavefield suppression module is configured to suppress the large-angle outgoing signals of the initial forward propagation wavefield and the reverse propagation wavefield, respectively. The propagation module is configured to continue propagating the initial forward propagation wavefield and reverse propagation wavefield after the large-angle outgoing signal is suppressed according to a given velocity model, so as to obtain the global forward propagation wavefield and reverse propagation wavefield. The imaging stage suppression module is configured to suppress large-angle outgoing signals in the global forward propagation wavefield and reverse propagation wavefield during the imaging stage, respectively. The imaging module is configured to generate multiple wave imaging results based on the suppressed wavefield of the large-angle outgoing signal during the imaging stage, according to preset imaging conditions.

[0044] Example 3 This embodiment provides a processing device corresponding to the multiple-wave imaging method for compression-stretch distortion provided in Embodiment 1. The processing device can be a client-side processing device, such as a mobile phone, laptop, tablet computer, desktop computer, etc., to execute the method of Embodiment 1.

[0045] The processing device includes a processor, a memory, a communication interface, and a bus. The processor, memory, and communication interface are connected via the bus to enable communication between them. The memory stores a computer program that can run on the processor. When the processor runs the computer program, it executes the multi-wave imaging method for compression-stretch distortion provided in Embodiment 1.

[0046] Preferably, the memory may be high-speed random access memory (RAM), and may also include non-volatile memory, such as at least one disk storage device.

[0047] Preferably, the processor can be any type of general-purpose processor such as a central processing unit (CPU) or a digital signal processor (DSP), and there is no limitation herein.

[0048] Example 4 The multiple wave imaging method for suppressing and stretching distortion in Embodiment 1 can be specifically implemented as a computer program product. The computer program product may include a computer-readable storage medium on which computer-readable program instructions for executing the multiple wave imaging method for suppressing and stretching distortion described in Embodiment 1 are loaded.

[0049] A computer-readable storage medium can be a tangible device that holds and stores instructions for use by an instruction execution device. A computer-readable storage medium can be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any combination thereof.

[0050] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes. These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes. These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0051] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of the present invention. Any modifications or equivalent substitutions that do not depart from the spirit and scope of the present invention should be covered within the protection scope of the claims of the present invention.

Claims

1. A method for suppressing stretching distortion in multiple-wave imaging, characterized in that, include: Obtain velocity models and initial wavefields from seismic data; Suppress the large-angle outgoing signals of the initial forward propagation wavefield and the reverse propagation wavefield respectively; Based on the given velocity model, the initial forward propagation wavefield and reverse propagation wavefield after the large-angle outgoing signal is suppressed continue to propagate to obtain the global forward propagation wavefield and reverse propagation wavefield. Suppress large-angle outgoing signals during the imaging stage of both the global forward propagation wavefield and the reverse propagation wavefield. Based on the suppression of the wave field by the large-angle outgoing signal during the imaging stage, multiple wave imaging results are generated according to preset imaging conditions.

2. The multiple-wave imaging method for suppressing stretching distortion as described in claim 1, characterized in that, The method of suppressing the large-angle emitted signals of the initial forward propagation wavefield and the reverse propagation wavefield respectively includes the following steps: ① For the initial forward propagation wave field Perform a time-domain Fourier transform to obtain the initial forward propagation wavefield in the frequency domain. ,in, Let be the coordinate point. For frequency; ② Initial forward propagating wave field in the frequency domain Perform a two-dimensional Fourier transform to obtain ,in, , They are respectively Wave number in direction; ③ Based on the surface velocity of the medium Calculate each frequency at each point on the sea surface. Corresponding wave number and the angle of departure ; ④ Based on a given threshold and the angle of departure The surface-emitted signal is suppressed to obtain the initial forward propagation wavefield after suppressing the large-angle-emitted signal. ; ⑤ The initial forward propagation wave field after suppressing the large-angle emitted signal Performing an inverse Fourier transform yields the initial forward propagation wavefield in the spatial domain that suppresses large-angle outgoing signals. ; ⑥ Using the same method as steps ① to ⑤, the initial reverse propagation wavefield is processed. Suppressing large-angle outgoing signals yields the initial backpropagation wavefield in the spatial domain for suppressing large-angle outgoing signals. .

3. The multiple-wave imaging method for suppressing stretching distortion as described in claim 2, characterized in that, The wave number The calculation formula is: Angle of departure The calculation formula is: 。 4. The multiple-wave imaging method for suppressing stretching distortion as described in claim 2, characterized in that, The based on a given threshold and the angle of departure Suppressing the surface emitted signal specifically refers to: suppressing signals exceeding a given threshold. The surface-emitted signal is suppressed to obtain the initial forward propagation wavefield after suppressing the large-angle emitted signal. The calculation formula is: 。 5. The multiple-wave imaging method for suppressing stretching distortion as described in claim 1, characterized in that, The large-angle outgoing signal suppression during the imaging stage of the global forward and reverse propagation wavefields includes: ① Global forward propagation wave field in the spatial domain Perform a two-dimensional Fourier transform to obtain ; ②Based on the velocity of the medium Calculate the angle of incidence at each point on the sea surface. ; ③ For emission angles greater than a given threshold The signal is suppressed to obtain the global forward propagation wavefield after suppressing the large-angle outgoing signal. ; ④ Global forward propagation wave field after suppressing large-angle emitted signals Performing an inverse Fourier transform yields the global forward propagation wavefield in the spatial domain after suppressing large-angle outgoing signals. ; ⑤ Using the same method as steps ① to ④, suppress the large-angle outgoing signal of the global backpropagation wavefield in the space domain to obtain the global backpropagation wavefield in the space domain after suppressing the large-angle outgoing signal. .

6. The multiple-wave imaging method for suppressing stretching distortion as described in claim 5, characterized in that, The exit angle The calculation formula is: The global forward propagation wave field after suppressing energy emitted at large angles The calculation formula is: 。 7. The multiple-wave imaging method for suppressing stretching distortion as described in claim 5, characterized in that, The formula for calculating the multiple wave imaging results is as follows: in, The stability parameter is calculated using the following formula: in, for . conjugate.

8. A multiple-wave imaging system for suppressing stretching distortion, characterized in that, include: The data acquisition module is configured to acquire velocity models and initial wavefields from seismic data; The initial wavefield suppression module is configured to suppress the large-angle outgoing signals of the initial forward propagation wavefield and the reverse propagation wavefield, respectively. The propagation module is configured to continue propagating the initial forward propagation wavefield and reverse propagation wavefield after the large-angle outgoing signal is suppressed according to a given velocity model, so as to obtain the global forward propagation wavefield and reverse propagation wavefield. The imaging stage suppression module is configured to suppress large-angle outgoing signals in the global forward propagation wavefield and reverse propagation wavefield during the imaging stage, respectively. The imaging module is configured to generate multiple wave imaging results based on the suppressed wavefield of the large-angle outgoing signal during the imaging stage, according to preset imaging conditions.

9. A computer-readable storage medium for storing one or more programs, characterized in that, The one or more programs include instructions that, when executed by a computing device, cause the computing device to perform any of the methods described in claims 1 to 7.

10. A computing device, characterized in that, include: One or more processors and a memory, wherein the memory stores one or more programs and is configured to be executed by the one or more processors, the one or more programs including instructions for performing any of the methods described in claims 1 to 7.