CN-121995493-A - Seismic data closure difference correction method, device, equipment and storage medium

CN121995493ACN 121995493 ACN121995493 ACN 121995493ACN-121995493-A

Abstract

The application discloses a method, a device, equipment and a storage medium for correcting a closing difference of seismic data, and belongs to the technical field of seismic exploration. The method comprises the steps of firstly determining a correlation coefficient between first seismic sub data and second seismic sub data, wherein the correlation coefficient can indicate the correlation degree between the seismic data of two work areas, so that the position of a splicing line determined according to the correlation coefficient can indicate the optimal splicing position of the seismic data of the two work areas. And then, based on the position of the stitching line, sequentially carrying out space integral closing difference correction and sample-point-by-sample dynamic closing difference correction on the first seismic data of the corrected working area, wherein the space integral closing difference correction carries out closing difference correction on the seismic channel layer, and the sample-point-by-sample dynamic closing difference correction carries out closing difference correction on the sampling point layer, so that the closing difference between the seismic data of the two working areas can be reduced to the greatest extent through the two correction processes, and the accuracy of the closing difference correction is improved.

Inventors

XIONG JINGXUAN
LIU HONG
CAO ZHONGLIN
HE GUANGMING
CHEN DAN
ZHOU QI

Assignees

中国石油天然气集团有限公司
中国石油集团东方地球物理勘探有限责任公司
中油油气勘探软件国家工程研究中心有限公司

Dates

Publication Date: 20260508
Application Date: 20241108

Claims (10)

1. A method of seismic data closure difference correction, the method comprising: acquiring first seismic data corresponding to each of a plurality of work areas in a target area; For a first working area and a second working area, determining correlation coefficients between first seismic sub-data and second seismic sub-data, wherein the first working area and the second working area are any two working areas with overlapping areas in the target area, the first seismic sub-data are data corresponding to a first seismic channel, the second seismic sub-data are data corresponding to a second seismic channel, the first seismic channel and the second seismic channel are seismic channels which are positioned at the same position in the overlapping area but belong to different working areas, the first seismic channel belongs to the first working area, and the second seismic channel belongs to the second working area; determining the positions of splicing lines corresponding to the first work area and the second work area based on the correlation coefficient; determining the first work area as a standard work area, and determining the second work area as a correction work area; based on the position of the splicing line, performing space whole channel closing difference correction on the first seismic data corresponding to the correction work area to obtain second seismic data; based on the first seismic sub data, carrying out sample-by-sample point dynamic closing difference correction on the second seismic data to obtain third seismic data; And based on the position of the splicing line, splicing the first seismic data corresponding to the standard work area and the third seismic data corresponding to the correction work area to obtain the integral seismic data corresponding to the target area.
2. The method of claim 1, wherein determining correlation coefficients between the first seismic sub-data and the second seismic sub-data comprises: determining a first average value based on the data of each first sampling point in the first seismic trace and the number of sampling points, wherein the first seismic sub-data comprises the data of each first sampling point in the first seismic trace; determining a second average value based on the data of each second sampling point in the second seismic trace and the sampling point number, wherein the second seismic sub-data comprises the data of each second sampling point in the second seismic trace; The correlation coefficient is determined based on the data of each first sampling point, the first average value, the data of each second sampling point, and the second average value.
3. The method of claim 1, wherein performing spatial integral closure difference correction on the first seismic data corresponding to the corrected work area based on the stitching line position to obtain second seismic data comprises: determining a correction time window corresponding to the correction work area based on the position of the splicing line; For any third seismic trace in the correction time window, constructing a first function, wherein the first function is used for expressing the relation among x coordinates, y coordinates and z coordinates of third sampling points, the third sampling points are any sampling points in the third seismic trace, the first seismic data corresponding to the correction work area comprises third seismic sub-data corresponding to each third seismic trace in the correction time window, and the third seismic sub-data corresponding to the third seismic trace comprises x coordinates, y coordinates and z coordinates of each third sampling point; Constructing a second function, wherein the second function is used for representing the relation among an x coordinate, a y coordinate and a z 'coordinate of the third sampling point, and the z' coordinate is obtained by smoothing the z coordinate of the third sampling point; Constructing a third function based on the z coordinate and the z 'coordinate of the third sampling point, wherein the third function is used for representing the relation between the displacement sum of a plurality of third sampling points corresponding to the third seismic channel and the z' coordinate of each third sampling point; Solving the third function, and performing space whole channel closing difference correction on third seismic sub-data corresponding to the third seismic channel based on a solving result to obtain fourth seismic sub-data corresponding to the third seismic channel; and determining the second seismic data based on fourth seismic sub-data corresponding to each third seismic channel in the correction time window and first seismic data which are positioned in the correction work area and are uncorrected outside the correction time window.
4. The method of claim 3, wherein the solving the third function, based on the solving result, performs spatial integral trace closure difference correction on the third seismic sub-data corresponding to the third seismic trace, to obtain the fourth seismic sub-data corresponding to the third seismic trace, including: Solving the third function by taking the minimum sum of the displacements as a target to obtain a z' coordinate of the third sampling point; And moving the z coordinate of each third sampling point based on the z' coordinate of each third sampling point of the third seismic trace to obtain fourth seismic sub-data corresponding to the third seismic trace.
5. The method of claim 1, wherein the second seismic data comprises fifth seismic sub-data corresponding to each second seismic trace within the overlap region and sixth seismic sub-data corresponding to each fourth seismic trace outside the overlap region; The performing sample-by-sample dynamic closing difference correction on the second seismic data based on the first seismic sub-data to obtain third seismic data, including: For each first sampling point in the first seismic trace, determining a target sampling point with the shortest distance from the first sampling point from a plurality of second sampling points in the second seismic trace, wherein the first seismic sub-data comprises x-coordinate, y-coordinate and z-coordinate of each first sampling point, and the fifth seismic sub-data comprises x-coordinate, y-coordinate and z-coordinate of each second sampling point; Moving the position of the target sampling point to the same position as the first sampling point to obtain seventh seismic sub-data; And determining the third seismic data based on the seventh seismic sub-data corresponding to each second seismic trace and the sixth seismic sub-data corresponding to each fourth seismic trace.
6. The method of claim 1, wherein acquiring the first seismic data corresponding to each of the plurality of work areas within the target area comprises: Acquiring post-stack seismic data corresponding to each of a plurality of work areas in the target area; according to an absolute coordinate system, performing gridding treatment on the post-stack seismic data corresponding to each of the plurality of work areas; And if the meshed seismic data corresponding to each of the plurality of work areas has the seismic data with inconsistent sampling rates, resampling the meshed seismic data corresponding to each of the plurality of work areas to obtain first seismic data corresponding to each of the plurality of work areas.
7. A seismic data closure difference correction apparatus, the apparatus comprising: The acquisition module is used for acquiring first seismic data corresponding to each of the plurality of work areas in the target area; The first determining module is used for determining correlation coefficients between first and second seismic sub-data for a first work area and a second work area, wherein the first work area and the second work area are any two work areas with overlapping areas in the target area, the first seismic sub-data are data corresponding to first seismic channels, the second seismic sub-data are data corresponding to second seismic channels, the first seismic channels and the second seismic channels are seismic channels which are positioned at the same position in the overlapping areas but belong to different work areas, the first seismic channels belong to the first work area, and the second seismic channels belong to the second work area; The second determining module is used for determining the positions of the splicing lines corresponding to the first work area and the second work area based on the correlation coefficient; a third determining module, configured to determine the first work area as a standard work area, and determine the second work area as a corrected work area; The first correction module is used for correcting the space whole channel closing difference of the first seismic data corresponding to the correction work area based on the position of the splicing line to obtain second seismic data; The second correction module is used for carrying out sample-by-sample dynamic closing difference correction on the second seismic data based on the first seismic sub data to obtain third seismic data; And the splicing module is used for splicing the first seismic data corresponding to the standard work area and the third seismic data corresponding to the correction work area based on the splicing line position to obtain the integral seismic data corresponding to the target area.
8. An electronic device comprising a processor and a memory, wherein the memory has stored therein at least one program code that is loaded and executed by the processor to implement the seismic data closure difference correction method of any of claims 1 to 6.
9. A computer readable storage medium having stored therein at least one program code loaded and executed by a processor to implement the seismic data closure difference correction method of any of claims 1 to 6.
10. A computer program product, characterized in that it has stored therein at least one program code, which is loaded and executed by a processor to implement the seismic data closure difference correction method according to any of claims 1 to 6.

Description

Seismic data closure difference correction method, device, equipment and storage medium Technical Field The present application relates to the field of seismic exploration technologies, and in particular, to a method, an apparatus, a device, and a storage medium for correcting a closing difference of seismic data. Background Along with the development of seismic exploration, the exploration force is continuously increased, the exploration target is deeper and finer, and the overall rule of in-band reservoir development, crack development characteristics, hydrocarbon source rock development conditions and the like in the centibasin area are more important to the overall knowledge of the area. The area generally comprises a plurality of blocks, the seismic data of the blocks are spliced together in a shortcut mode to form unified overall data, the interpretation precision and efficiency of the area can be improved, the fine patch processing of the local area can be carried out by the vector on the basis of overall knowledge, and the problem that the full basin or the full zone patch processing can not be completed due to the conventional hardware and processing time is solved. Because of the difference of the acquisition mode and the processing mode, the closing difference exists among the seismic data of a plurality of blocks, and therefore, the closing difference correction is required to be carried out on the seismic data of the plurality of blocks before the splicing. Disclosure of Invention The embodiment of the application provides a method, a device, equipment and a storage medium for correcting a closing difference of seismic data, which can improve the accuracy of correcting the closing difference of the seismic data. The technical scheme is as follows: In one aspect, a method for seismic data closure difference correction is provided, the method comprising: acquiring first seismic data corresponding to each of a plurality of work areas in a target area; For a first working area and a second working area, determining correlation coefficients between first seismic sub-data and second seismic sub-data, wherein the first working area and the second working area are any two working areas with overlapping areas in the target area, the first seismic sub-data are data corresponding to a first seismic channel, the second seismic sub-data are data corresponding to a second seismic channel, the first seismic channel and the second seismic channel are seismic channels which are positioned at the same position in the overlapping area but belong to different working areas, the first seismic channel belongs to the first working area, and the second seismic channel belongs to the second working area; determining the positions of splicing lines corresponding to the first work area and the second work area based on the correlation coefficient; determining the first work area as a standard work area, and determining the second work area as a correction work area; based on the position of the splicing line, performing space whole channel closing difference correction on the first seismic data corresponding to the correction work area to obtain second seismic data; based on the first seismic sub data, carrying out sample-by-sample point dynamic closing difference correction on the second seismic data to obtain third seismic data; And based on the position of the splicing line, splicing the first seismic data corresponding to the standard work area and the third seismic data corresponding to the correction work area to obtain the integral seismic data corresponding to the target area. In another aspect, there is provided a seismic data closure difference correction apparatus, the apparatus comprising: The acquisition module is used for acquiring first seismic data corresponding to each of the plurality of work areas in the target area; The first determining module is used for determining correlation coefficients between first and second seismic sub-data for a first work area and a second work area, wherein the first work area and the second work area are any two work areas with overlapping areas in the target area, the first seismic sub-data are data corresponding to first seismic channels, the second seismic sub-data are data corresponding to second seismic channels, the first seismic channels and the second seismic channels are seismic channels which are positioned at the same position in the overlapping areas but belong to different work areas, the first seismic channels belong to the first work area, and the second seismic channels belong to the second work area; The second determining module is used for determining the positions of the splicing lines corresponding to the first work area and the second work area based on the correlation coefficient; a third determining module, configured to determine the first work area as a standard work area, and determine the second work area as a corrected work area; The first co