Disclosure of Invention
The embodiment of the invention provides a seismic profile display method, which is used for effectively solving the problems that the same phase axis is discontinuous and the real seismic trend can not be expressed when a thin reservoir profile with a larger inclination angle is displayed, and comprises the following steps:
performing dip angle analysis on seismic data of a preset work area to obtain dip angle data;
performing reverse time shifting on the front and rear seismic channel data of a preset seismic channel in a preset work area according to the inclination data, and flattening the seismic event axis;
interpolation processing is carried out between two adjacent seismic channels along the horizontal direction, so that each pixel point between two adjacent sampling points on the two adjacent seismic channels obtains a corresponding value;
setting colors for the values corresponding to each pixel point by using a preset color table;
and restoring the data of the front and rear seismic channels of the preset seismic channel according to the dip angle data to finish the seismic section display.
Optionally, the restoring the seismic channel data before and after the preset seismic channel according to the dip angle data to complete the seismic section display includes:
restoring the front and rear seismic channel data of the preset seismic channel according to the dip angle data;
and carrying out interpolation processing on the uncovered pixel points along the longitudinal direction according to a preset color table to finish the seismic section display.
Optionally, the method for performing dip analysis on the seismic data of the preset work area is a plane wave decomposition algorithm.
Optionally, the dip angle data is a time difference between a sampling point on a preset seismic channel and a corresponding sampling point on the same reflection interface on the next seismic channel of the preset seismic channel.
The embodiment of the invention also provides a seismic profile display device, which is used for effectively solving the problems that the same phase axis is discontinuous and the real seismic trend can not be expressed when the thin reservoir profile with a larger inclination angle is displayed, and comprises the following components:
the data analysis module is used for carrying out dip angle analysis on the seismic data of the preset work area to obtain dip angle data;
the time shifting module is used for carrying out reverse time shifting on the front and rear seismic channel data of a preset seismic channel in a preset work area according to the inclination data and flattening the seismic event-phase axis;
the interpolation processing module is used for carrying out interpolation processing between two adjacent seismic channels along the horizontal direction so as to enable each pixel point between two adjacent sampling points on the two adjacent seismic channels to obtain a corresponding value;
the color setting module is used for setting colors for the values corresponding to the pixel points by utilizing a preset color table;
and the seismic section display module is used for restoring the data of the front and rear seismic channels of the preset seismic channel according to the dip angle data to finish the seismic section display.
Optionally, the seismic profile display module is further configured to:
restoring the front and rear seismic channel data of the preset seismic channel according to the dip angle data;
and carrying out interpolation processing on the uncovered pixel points along the longitudinal direction according to a preset color table to finish the seismic section display.
Optionally, the method for performing dip analysis on the seismic data of the preset work area is a plane wave decomposition algorithm.
Optionally, the dip angle data is a time difference between a sampling point on a preset seismic channel and a corresponding sampling point on the same reflection interface on the next seismic channel of the preset seismic channel.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above method is stored.
In the embodiment of the invention, dip angle analysis is carried out on seismic data of a preset work area to obtain dip angle data, backward time shifting is carried out on front and back seismic channel data of a preset seismic channel in the preset work area according to the dip angle data, the seismic event is leveled, interpolation processing is carried out between two adjacent seismic channels along the horizontal direction to enable each pixel point between two adjacent sampling points on the two adjacent seismic channels to obtain a corresponding value, a preset color table is used for setting the color of the corresponding value of each pixel point, and finally, restoration processing is carried out on the front and back seismic channel data of the preset seismic channel according to the dip angle data, so that seismic section display can be completed. The whole process is subjected to data time shifting, horizontal difference value and time shifting restoration, the problems that the same phase axis is discontinuous when a thin reservoir section with a larger inclination angle is displayed and the real earthquake trend cannot be expressed can be effectively solved, and the earthquake interpretation precision is higher.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
Fig. 1 is a flowchart of a seismic profile display method according to an embodiment of the present invention, and as shown in fig. 1, the method includes:
step 101, performing dip angle analysis on seismic data of a preset work area to obtain dip angle data.
In this embodiment, the method for performing dip analysis on seismic data of a preset work area is a plane wave decomposition algorithm. The dip angle data is the time difference between a sampling point on a preset seismic channel and a corresponding sampling point on the same reflection interface on the next seismic channel of the preset seismic channel (see the attached figure 3), and the dip angle data is consistent with the seismic profile data range.
And 102, performing reverse time shifting on the front and rear seismic channel data of the preset seismic channel in the preset work area according to the inclination data, and flattening the seismic event.
In specific implementation, for example, when displaying a seismic section between the nth trace and the (n +1) th trace, the seismic data needs to be preprocessed by the dip data: and (3) the nth data is fixed, and the (n +1) th, n +2 th and n-1 th seismic channel data are subjected to reverse time shifting according to the dip angle, namely, the dip angle is reset to zero, and the seismic event axis is leveled, and specifically, the method can refer to the attached figures 4 and 5, wherein the attached figure 4 is a schematic diagram of the pre-processed seismic channel of the work area, and the attached figure 5 is a schematic diagram of the pre-processed seismic channel of the work area.
Wherein, the n +1 th seismic data time shift Tshift1Comprises the following steps:
wherein I is the sampling point, IsampleIs the sampling interval, dip [ i ]]The inclination angle of the ith sampling point of the nth channel is calculated by using a plane wave decomposition algorithmThe time difference.
The result data11 after the time shift of the n +1 th seismic data is:
data11[i]=data1[i+Tshift1] (2)
the result data12 after the time shift of the n-1 st seismic data is:
data12[i]=data1[i-Tshift1] (3)
after the time shift, the n +2 th seismic data2 is:
wherein dip1[ i ] is the tilt moveout of the ith sampling point of the (n +1) th channel calculated by using a plane wave decomposition algorithm.
And 103, performing interpolation processing between two adjacent seismic channels along the horizontal direction to enable each pixel point between two adjacent sampling points on the two adjacent seismic channels to obtain a corresponding value.
Instep 102, in the implementation, interpolation processing is performed in the horizontal direction using formulas (1), (2), (3), and (4), and three data passes are interpolated between every two passes, as shown in fig. 6. And then, continuously interpolating between the nth and the n +1 th tracks by using the three interpolated data, so that each pixel point between the ith and the i +1 th sampling points between the nth and the n +1 th tracks on the screen obtains a corresponding value.
And 104, setting colors for the corresponding values of each pixel point by using a preset color table.
In the embodiment, for the variable density color display of the seismic data, the reason is to set each pixel of the screen area of the seismic display to be the corresponding color. When setting the color, the order from left to right and from top to bottom is adopted, and the color is set for the value corresponding to each pixel point by interpolation between every two channels and between every two sampling points and by using a preset color table, as shown in fig. 7.
And 105, restoring the front and rear seismic channel data of the preset seismic channel according to the dip angle data to finish the seismic section display.
In specific implementation, based on the above steps, after the pixel between the i-th sampling point and the i + 1-th sampling point between the n-th track and the n + 1-th track is fully inserted, the data position is restored, that is, the data between the sampling points (i and i +1) is shifted to the actual position according to the inclination direction (see fig. 8). When moving, the horizontal x value and the vertical y value of the pixel coordinate are not changedpixiThe movement is as follows:
wherein, yoripixiFor the value y interpolated instep 103, dip [ i ]]Is the inclination angle time difference, x, of the ith sampling point of the nth channel calculated by a plane wave decomposition algorithm1Is the x-direction pixel coordinate of the n-th data, x2X-direction pixel coordinate, y, of the n +1 th datastartIs the y-direction pixel coordinate of the ith sampling point, yendIs the y-direction pixel coordinate of the (i +1) th sampling point.
In this embodiment,step 105 includes:
restoring the front and rear seismic channel data of the preset seismic channel according to the dip angle data;
and (4) carrying out interpolation processing on the uncovered pixel points along the longitudinal direction according to a preset color table to finish seismic section display, and referring to fig. 9.
The same method is adopted for other seismic channels, and the steps are repeated until the seismic display area of the screen is fully paved.
Further, for boundary lanes, for example, for the left boundary: virtualizing n-1 tracks to be consistent with n tracks of data, and for a right boundary: virtualizing n +2 tracks to be consistent with n +1 tracks, and then processing the above steps.
As can be seen from fig. 1, in the seismic section display method provided in the embodiment of the present invention, dip angle analysis is performed on seismic data of a preset work area to obtain dip angle data, reverse time shift is performed on front and rear seismic channel data of a preset seismic channel in the preset work area according to the dip angle data, a seismic event is leveled, interpolation processing is performed between two adjacent seismic channels along a horizontal direction to obtain a corresponding value for each pixel point between two adjacent sampling points on the two adjacent seismic channels, a preset color table is used to set a color for the value corresponding to each pixel point, and finally, the front and rear seismic channel data of the preset seismic channel are restored according to the dip angle data, so that seismic section display can be completed. The whole process is subjected to data time shifting, horizontal difference value and time shifting restoration, the problems that the same phase axis is discontinuous when a thin reservoir section with a larger inclination angle is displayed and the real earthquake trend cannot be expressed can be effectively solved, and the earthquake interpretation precision is higher.
Based on the same inventive concept, embodiments of the present invention further provide a seismic profile display apparatus, as described in the following embodiments. The principle of the seismic profile display device for solving the problems is similar to that of the seismic profile display method, so the implementation of the seismic profile display device can be referred to the implementation of the seismic profile display method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
Fig. 2 is a schematic structural diagram of a seismic profile display device according to an embodiment of the present invention, and as shown in fig. 2, the device includes:
thedata analysis module 201 is configured to perform dip angle analysis on the seismic data of the preset work area to obtain dip angle data.
And themigration module 202 is configured to perform reverse time shifting on the front and rear seismic channel data of the preset seismic channel in the preset work area according to the inclination data, and flatten the seismic event.
And theinterpolation processing module 203 is configured to perform interpolation processing between two adjacent seismic channels along the horizontal direction, so that each pixel point between two adjacent sampling points on the two adjacent seismic channels obtains a corresponding value.
And acolor setting module 204, configured to set a color for the corresponding value of each pixel point by using a preset color table.
And the seismicsection display module 205 is configured to perform reduction processing on the seismic channel data before and after the preset seismic channel according to the dip angle data, so as to complete seismic section display.
In an embodiment of the present invention, the seismicprofile display module 205 is further configured to:
restoring the front and rear seismic channel data of the preset seismic channel according to the dip angle data;
and carrying out interpolation processing on the uncovered pixel points along the longitudinal direction according to a preset color table to finish the seismic section display.
In the embodiment of the invention, the method for carrying out dip analysis on the seismic data of the preset work area is a plane wave decomposition algorithm.
In the embodiment of the invention, the dip angle data is the time difference between a sampling point on a preset seismic channel and a sampling point on the next seismic channel of the preset seismic channel.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the above method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above method is stored.
In summary, the invention performs dip analysis on seismic data of a preset work area to obtain dip data, performs reverse time shift on front and rear seismic channel data of a preset seismic channel in the preset work area according to the dip data, flattens a seismic event, performs interpolation processing between two adjacent seismic channels along a horizontal direction to enable each pixel point between two adjacent sampling points on the two adjacent seismic channels to obtain a corresponding value, sets a color for the corresponding value of each pixel point by using a preset color table, and finally performs reduction processing on the front and rear seismic channel data of the preset seismic channel according to the dip data, thereby completing seismic section display. The whole process is subjected to data time shifting, horizontal difference value and time shifting reduction, the problems that the same phase axis is discontinuous when a thin reservoir section with a larger inclination angle is displayed and the real earthquake trend cannot be expressed can be effectively solved, the consistency with the actual geological condition is higher, the earthquake interpretation precision is higher, and a foundation is provided for follow-up work.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may 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, and the like) having computer-usable program code embodied therein.
The present 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 flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may 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, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.