INTERPRETIVE-GUIDED VELOCITY MODELING SEISMIC IMAGING METHOD AND SYSTEM, MEDIUM AND DEVICE

Abstract

The present disclosure belongs to the technical field of seismic exploration imaging, and relates to an interpretive-guided velocity modeling seismic imaging method and system, a medium and a device. The method comprises the following steps: S1. performing first imaging on a given initial velocity model to obtain a first imaging result; S2. performing relative wave impedance inversion on the first imaging result to obtain a relative wave impedance profile; S3. performing Curvelet filtering on the relative wave impedance profile to obtain a first interpretation scheme; S4. superposing the first interpretation scheme and the initial velocity model to obtain a new migration velocity field; S5. performing second imaging on a new migration velocity field to obtain a second imaging result; and S6. repeating steps S2-S4 for the obtained second imaging result until a final seismic imaging result is obtained.

Claims

1. An interpretive-guided velocity modeling seismic imaging method, comprising: performing first imaging on a given initial velocity model to obtain a first imaging result; performing relative wave impedance inversion on the first imaging result to obtain a relative wave impedance profile; performing Curvelet filtering on the relative wave impedance profile to obtain a first interpretation scheme; superposing the first interpretation scheme and the given initial velocity model to obtain a new migration velocity field; performing second imaging on the new migration velocity field to obtain a second imaging result; and repeating the steps of relative wave impedance inversion and Curvelet filtering for the obtained second imaging result until a final seismic imaging result is obtained, the final seismic imaging result is a clear image of a real geological structure.

2. The interpretive-guided velocity modeling seismic imaging method of claim 1, wherein the first imaging and the second imaging are obtained by respectively inputting the given initial velocity model, and the new migration velocity field into a least-square inverse-time migration algorithm.

3. The interpretive-guided velocity modeling seismic imaging method of claim 1, wherein the relative wave impedance inversion is directly performed based on a deconvolution method.

4. The interpretive-guided velocity modeling seismic imaging method of claim 3, wherein a relative wave impedance is calculated on the relative wave impedance profile by the relative wave impedance inversion method, and a standard impedance is obtained by performing migration velocity analysis, the standard impedance is normalized and calibrated to obtain a relative velocity profile.

5. The interpretive-guided velocity modeling seismic imaging method of claim 1, wherein a similarity, between the imaging result obtained each time and an imaging result of a real fault-karst model, is calculated to verify whether or not, each of the obtained imaging result is similar to an imaging result of each interpretive-guided velocity modeling.

6. The interpretive-guided velocity modeling seismic imaging method of claim 1, wherein the given initial velocity model is obtained by conventional velocity modeling, the sizes and number of horizontal and vertical grids of the initial velocity model are given, and shot gather records based on the initial velocity model are obtained by a finite difference method.

7. The interpretive-guided velocity modeling seismic imaging method of claim 6, wherein the initial velocity model comprises the following model parameters: sizes of the horizontal and vertical grids, spacing of the horizontal and vertical grids, wavelet duration and main frequency, time sampling interval, total time length, number of seismic sources, spacing between the seismic sources and starting positions of horizontal and vertical coordinates of the seismic sources.

8. An interpretive-guided velocity modeling seismic imaging system, comprising: a primary imaging module configured for performing first imaging on a given initial velocity model to obtain a first imaging result; a relative wave impedance inversion module configured for performing relative wave impedance inversion on the first imaging result to obtain a relative wave impedance profile; an interpretation module configured for performing Curvelet filtering on the relative wave impedance profile to obtain a first interpretation scheme; a superposing module configured for superposing the first interpretation scheme and the initial velocity model to obtain the new migration velocity field; a second imaging module configured for performing second imaging on the new migration velocity field to obtain a second imaging result; and a circulation module configured for inputting the obtained second imaging result into the relative wave impedance inversion module and circulating the obtained second imaging result in the interpretation module until the final seismic imaging result is obtained.

9. A computer-readable storage medium storing one or more programs, wherein the one or more programs comprise instructions, and when the instructions are executed by a computing device, the computing device executes the interpretive-guided velocity modeling seismic imaging method according to claim 1.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0018] FIG. 1 is a flowchart of an interpretive-guided velocity modeling seismic imaging system in an embodiment of the present disclosure;

[0019] FIG. 2 is a real fault-karst velocity model in an embodiment of the present disclosure;

[0020] FIG. 3 is a least-square reverse time migration imaging profile of the real fault-karst velocity model in an embodiment of the present disclosure;

[0021] FIG. 4 is a migration velocity field of the real fault-karst velocity model in an embodiment of the present disclosure;

[0022] FIG. 5 is a first imaging result diagram of the migration velocity field of the real fault-karst velocity model in an embodiment of the present disclosure;

[0023] FIG. 6 is a relative wave impedance profile of the first imaging result diagram of the real fault-karst velocity model in an embodiment of the present disclosure;

[0024] FIG. 7 is an interpretation scheme diagram based on the relative wave impedance profile of the real fault-karst velocity model in an embodiment of the present disclosure;

[0025] FIG. 8 is a superposition result of the interpretation scheme and the migration velocity field, and the updated image of the migration velocity field in an embodiment of the present disclosure; and

[0026] FIG. 9 is a second imaging result diagram in an embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

Embodiment I

[0027] The present embodiment discloses an interpretive-guided velocity modeling seismic imaging method and is illustrated by taking a fault-karst velocity model as an example; and a real model image of the fault-karst velocity model is as shown in FIG. 2. Synthetic seismic shot gather data is generated based on the fault-karst velocity model, and a real velocity model migration imaging profile as shown in FIG. 3 is obtained by adopting the least-square reverse time migration (RTM—Reverse Time Migration) imaging. As shown in FIG. 1, the method comprises the following steps:

[0028] S1. performing first imaging on a given initial velocity model to obtain a first imaging result;

[0029] The given initial velocity model is as shown in FIG. 4 a ad is obtained by conventional velocity modeling, the sizes and number of horizontal and vertical grids of the initial velocity model are given, and shot-gather records based on the initial velocity model are obtained by a finite difference method.

[0030] The initial velocity model comprises the following model parameters: sizes of the horizontal and vertical grids, spacing of the horizontal and vertical grids, wavelet duration and main frequency, time sampling interval, total time length, number of seismic sources, spacing between the seismic sources and starting positions of horizontal and vertical coordinates of the seismic sources.

[0031] S2. performing relative wave impedance inversion on the first imaging result to obtain a relative wave impedance profile.

[0032] The initial velocity model is imaged by adopting a least-square inverse time migration method based on the parameters of the migration velocity model and the shot-gather data to obtain the first imaging result, and the first imaging result in the present embodiment is as shown in FIG. 5, wherein the parameters of the model are the same as those of the initial velocity model.

[0033] Relative wave impedance inversion is directly performed based on a deconvolution method, which has the characteristics of simple calculation and being not restricted by geological data and well data. Specifically, the relative wave impedance is calculated on the relative wave impedance profile, the standard impedance is obtained by migration velocity analysis, and is normalized and calibrated to obtain a relative velocity profile, specifically as shown in FIG. 6, the relative wave impedance is obtained.

[0034] S3. performing Curvelet filtering on the relative wave impedance profile to highlight the geological structure and obtain a first interpretation scheme, the rest of which is as shown in FIG. 7.

[0035] S4. superposing the first interpretation scheme and the initial velocity model to obtain a new migration velocity field, as shown in FIG. 8.

[0036] S5. performing second imaging on the new migration velocity field to obtain the second imaging result, as shown in FIG. 9.

[0037] The new migration velocity field is imaged by adopting the least-square reverse time migration method based on the parameters of the migration velocity model and the shot-gather data to obtain the second imaging result.

[0038] S6. repeating S2-S4 for the obtained second imaging result until the final seismic imaging result is obtained. Through such iteration, the seismic imaging profile gradually approaches the real geological structure.

[0039] The similarity between the imaging result obtained each time and the imaging result of the real fault-karst model is calculated to verify the imaging result of each interpretive-guided velocity modeling and judge the effectiveness of the present method. The similarities between the first imaging result and the imaging result of the real fault-karst model and between the second imaging result and the imaging result of the real fault-karst model (that is, a real velocity imaging profile obtained based on the image in FIG. 2, that is, FIG. 3) are respectively calculated. The greater the similarity, the better the effect of the interpretive-guided velocity modeling seismic imaging method is. The similarity between the first imaging profile and the real velocity imaging profile is 71.46%, and the similarity between the second imaging profile and the real velocity imaging profile is 81.28%. After two iterations, the imaging effect is improved by nearly 10%.

Embodiment II

[0040] Based on the same inventive concept, the present disclosure provides an interpretive-guided velocity modeling seismic imaging system, which comprises:

[0041] a primary imaging module configured for performing first imaging on the given initial velocity model to obtain the first imaging result;

[0042] a relative wave impedance inversion module configured for performing relative wave impedance inversion on the first imaging result to obtain the relative wave impedance profile;

[0043] an interpretation module configured the performing Curvelet filtering on the relative wave impedance profile to obtain the first interpretation scheme;

[0044] a superposing module configured for superposing the first interpretation scheme and the initial velocity model to obtain the new migration velocity field;

[0045] a second imaging module configured for performing second imaging on the new migration velocity field to obtain a second imaging result; and

[0046] a circulation module configured for inputting the obtained second imaging result into the relative wave impedance inversion module and the interpretation module for circulation until the final seismic imagine result is obtained.

Embodiment III

[0047] Based on the same inventive concept, the present disclosure provides a computer-readable storage medium storing one or more programs, the one or more programs comprise instructions, and when the instructions are executed by a computing device, the computing device executes the interpretive-guided velocity modeling seismic imaging method of any one of the above.

Embodiment IV

[0048] Based on the same inventive concept, the present disclosure provides the computing device, comprising one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs are used for executing the interpretive-guided velocity modeling seismic imaging method according to any of the above.