A method and system for fusion processing of repeatedly collected seismic data

Abstract

The present invention provides a method for fusing and processing repeatedly acquired seismic data, belonging to the field of oil and gas geophysical exploration. The method for fusing and processing repeatedly acquired seismic data provided by the present invention comprises the following steps: S1, acquiring at least two acquisition data of the same position, and inverting velocity models with different physical points and calculating static corrections accordingly; S2, performing CMP stacking on the corresponding acquisition data using the static corrections to obtain corresponding CMP stacking data; S3, calculating the relative time difference between at least two CMP stacking data using cross-correlation; S4, applying the relative time difference decomposition to the static corrections described in step S1 to obtain the corrected static corrections; S5, using the corrected static corrections to stack and fuse the at least two acquisition data. The present invention has good application effects in the fields of time-shift seismic processing, repeated acquisition fusion processing, multi-block concatenation processing, etc.

Description

Technical Field

[0001] The invention belongs to the field of oil and gas geophysical exploration, and in particular relates to a method and system for fusing and processing repeatedly acquired seismic data. Background Art

[0002] In current earth exploration, fusion processing is to fuse two or more three-dimensional seismic data (overlapping acquisition at different times) of the same exploration area together for processing. Time-delay seismic data processing is to synchronously process three-dimensional seismic data (overlapping acquisition at different times) acquired successively, and analyze and monitor the dynamic information of reservoir exploitation by comparing the subtle attribute changes of the processing results of three-dimensional seismic data acquired successively. Joint processing is to process multiple three-dimensional data (acquired adjacently at different times) together. The common point of these three is that they all need to eliminate the time difference problem between data acquired at different times. In the process of time-delay seismic processing, fusion processing, and joint processing, the time difference between the two acquisitions will cause non-co-phase superposition, and the slight time difference will also reduce the signal-to-noise ratio and resolution of the processing results, so it must be eliminated.

[0003] The source type, excitation well depth, observation system, etc. of the seismic data collected successively often differ, which will cause the time difference between the seismic data collected successively. In addition, due to seasonal changes, groundwater exploitation and other reasons, the water table changes, and the travel time of seismic waves in the weathering layer is different, which will also cause the time difference between the data collected successively.

[0004] There are currently two main methods for correcting the time difference between repeated acquisitions. The first is to time-shift one of the blocks as a whole. This method can correct the overall time difference recorded by the instrument between the two acquisitions, but it cannot correct the spatially varying time difference caused by changes in near-surface velocity. The second method is non-rigid matching technology. This technology is based on graphic image processing technology. By comparing the data of the two acquisitions, the data of one of the acquisitions is used as the objective function to correct the seismic data of the other acquisition. This technology is usually used for consistency processing of post-stack data of time-shift seismic, and the method process has no clear geophysical significance. Summary of the invention

[0005] In order to solve the above technical problems, the present invention provides a method for fusion processing of repeatedly acquired seismic data, comprising the following steps:

[0006] S1, obtain at least two acquisition data at the same location, and use different physical points to invert the velocity model and calculate the static correction amount based on it;

[0007] S2, using the static correction amount to perform CMP stacking on the corresponding collected data to obtain corresponding CMP stacking data;

[0008] S3, calculating the relative time difference between at least two CMP stacking data by cross-correlation;

[0009] S4, decomposing the relative time difference and applying it to the static correction amount described in step S1 to obtain a corrected static correction amount;

[0010] S5. Using the corrected static correction amount, superimpose and fuse the at least two acquisition data.

[0011] In step S1, at least two acquisition data of the same position are distinguished by different gun station numbers so that they have different physical points.

[0012] Wherein, after step S2 and before step S3, the method further includes: comparing at least two CMP overlay data, and performing correction if there is a time difference.

[0013] Among them, the step S3 specifically includes: calculating the cross-correlation function between at least two of the CMP stacking data; picking up the time corresponding to the maximum amplitude value of the cross-correlation function to obtain the relative time difference between two acquisitions of each CMP; calculating the spatial weighting coefficient of the time difference, and calculating the total time difference of each CMP stacking data in combination with the relative time difference, which is recorded as the correction amount of each CMP stacking data.

[0014] Wherein, the step S4 specifically includes: comparing the inversion velocity models corresponding to the at least two acquisition data, determining whether the cause of the acquisition time difference is the shot point well depth correction or the near-surface velocity; if it is caused by the shot point well depth correction, then directly decomposing the correction amount of the CMP stacking data to the corresponding shot point; if it is caused by the near-surface velocity, then distributing the correction amount of the CMP stacking data according to the travel time of the weathering layer of the shot inspection point, and combining the static correction amount to obtain the corrected static correction amount.

[0015] After step S5, the method further includes: determining whether the effect after superposition and fusion is better than the CMP superposition data, and if so, terminating the process; if not, returning to step S3.

[0016] The present invention also provides a system for repeatedly acquiring seismic data and fusing and processing the seismic data, comprising:

[0017] Static correction calculation unit: used to obtain at least two acquisition data at the same location, and invert the velocity model with different physical points and calculate the static correction accordingly;

[0018] A CMP stacking data acquisition unit: performs CMP stacking on the corresponding collected data using the static correction amount to acquire corresponding CMP stacking data;

[0019] A cross-correlation calculation unit: using cross-correlation to calculate the relative time difference between at least two CMP stacking data;

[0020] A static correction value correction unit: applying the relative time difference decomposition to the static correction value in the static correction value calculation unit to obtain a corrected static correction value;

[0021] A superposition and fusion unit is used to superimpose and fuse the at least two acquisition data using the corrected static correction amount.

[0022] Wherein, in the static correction calculation unit, at least two acquisition data of the same position are distinguished by different gun station numbers so that they have different physical points.

[0023] Among them, the following steps are run in the cross-correlation calculation unit: calculate the cross-correlation function between at least two of the CMP stacking data; pick up the time corresponding to the maximum amplitude value of the cross-correlation function to obtain the relative time difference between two acquisitions of each CMP; calculate the spatial weighting coefficient of the time difference, and calculate the total time difference of each CMP stacking data in combination with the relative time difference, which is recorded as the correction amount of the CMP stacking data.

[0024] The present invention also provides a computer-readable storage medium, characterized in that: the computer-readable storage medium stores at least one program executable by a computer, and when the at least one program is executed by the computer, the computer executes the steps in the method for fusion processing of repeatedly acquired seismic data as described in any one of claims 1-6.

[0025] Compared with the prior art, the present invention has the following beneficial effects:

[0026] The present invention aims at the time difference problem between repeatedly collected seismic data, based on the near-surface velocity model obtained by tomographic inversion, according to the method and principle of static correction calculation, starting from the static correction calculation of the shot check point, and correcting the time difference problem between repeatedly collected data from the source, thereby improving the in-phase superposition of seismic data and effectively improving the seismic imaging effect. The present invention can be applied to the fusion processing of repeatedly collected data, enabling each collected data to be superimposed in phase, thereby improving the data quality; when processing multiple blocks in a piecemeal manner, it can improve the consistency of the piecemeal data, effectively improve the splicing boundary problem, and reduce the structural illusion of the splicing area, and has good application effects in the data processing fields of both.

[0027] The processing method adopted by the present invention for the time difference problem between repeatedly collected seismic data can obtain more reasonable structural imaging in time-shift seismic processing than traditional non-rigid matching technology and also has good application effect. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1Flow chart of the method for fusion processing of repeatedly acquired seismic data provided by the present invention;

[0029] Figure 2 The present invention provides a near-surface velocity model comparison of two acquisition inversions of repeated acquisition seismic data;

[0030] Figure 3 The time difference correction method for repeated acquisition of seismic data provided by the present invention adopts space-varying time difference calculated by cross-correlation;

[0031] Figure 4 The time difference correction method for repeated acquisition of seismic data provided by the present invention uses a stacking comparison diagram before and after the time difference correction;

[0032] Figure 5 A structural diagram of a system for fusion processing of repeatedly acquired seismic data provided by the present invention. DETAILED DESCRIPTION

[0033] The present invention is further described in detail below with reference to specific embodiments.

[0034] [Example 1]

[0035] The purpose of this patent is to address the problem of time difference or weak structural difference between the currently repeatedly collected seismic data. Based on the near-surface velocity model, the space-varying time difference caused by changes in near-surface conditions between at least two acquisitions is calculated, and the surface consistency is decomposed into static correction quantities for shot points and receiver points.

[0036] Taking the fusion processing of a work area repeatedly collected once as an example, the processing flow implemented by the present invention is as follows: Figure 1 As shown, the principle and implementation steps of the method are as follows:

[0037] S1. Acquire the data collected twice at the same location, and use different physical points to invert the velocity model and calculate the static correction amount accordingly; distinguish the two data collected at the same location by different gun station numbers so that they have different physical points.

[0038] Unify the data of the two acquisitions. When the excitation and reception of the two acquisitions are at the same coordinate position, the positions of the shot point and the detection point are the same during the two acquisitions. In order to distinguish the two acquisitions, we add a number of acquisitions before the shot detection station number. For example, for the first acquisition, the original 30015011 shot station number is changed to 130015011, and the second acquisition station number at the same position is changed to 230015011 (you can also modify only one of the acquisition records). After such arrangement, the records of the points at the same physical position acquired twice can be used to invert the velocity model and calculate the static correction using two different physical points (shot points or detection points at the same coordinate position). The two acquisitions are recorded separately, and the same shot station is recorded as a different shot station number in the two records. That is, the two acquisitions at the same position are distinguished by different shot station numbers, so that they have different physical points.

[0039] The near-surface velocity inversion is performed for each data collection. Pick the first arrival, use the first arrival wave tomography technology, use the same tomography method and parameters for each data collection, and perform near-surface velocity inversion for each data collection. The inverted near-surface velocity model is recorded as M 1 、M 2 .like Figure 2 The figure shows the velocity inversion results of near-surface tomography from two acquisitions at the same location. It can be seen that the velocities inverted from the two acquisitions are different, such as Figure 2 As shown, there is a difference in speed between 0.2s and 0.3s.

[0040] S2. According to the result of step 2, the static correction is calculated respectively. The static correction of the two collected shot points is recorded as S 1s , S 2s The static correction of the two acquisition points is denoted as S 1r , S 2r Each time data is collected, a seismic slope surface can be obtained. Each seismic slope surface is composed of several cmps, and each cmp is composed of many shot checks (offset distances). The static correction amount is used to perform CMP stacking on each collected data. The CMP stacking data is recorded as D 1 and D 2 , compare D 1 and D 2 , if D 1 and D 2 If there is a time difference, correction is required.

[0041] S3. Calculate the time difference using cross-correlation. (e.g. Figure 3 As shown in the figure, the part of the data or graph in the figure is the time difference) calculates the CMP superposition data twice (D 1 With D 2), pick up the time corresponding to the maximum amplitude value of the cross-correlation function, and obtain the relative time difference T between the two acquisitions of each CMP (common center point gather) xy (x, y are coordinates, CMP position, unit is meter). When calculating the cross-correlation function, if there is a difference in seismic wavelets between the two acquisitions, it is necessary to apply shaping filtering before cross-correlation to eliminate the influence of the wavelets in advance. In order to improve the accuracy of cross-correlation, the main frequency filter is used to improve D 1 With D 2 A high signal-to-noise ratio is also necessary (preprocessing).

[0042] Calculate the spatial weighting coefficient Q of the time difference xy (x, y are coordinates), and the time difference is reasonably distributed to the two acquisitions. The weight coefficient is expressed as a percentage. If D 1 The weight coefficient is Q xy , corresponding to D 2 1-Q xy , then for each CMP, D 1 With D 2 The total time difference allocated is T xy *Q xy and T xy *(1-Q xy ), denoted as T 1 and T 2 , where T xy =T 1 +T 2 The weight coefficient should be determined according to the data quality. 1 With M 2 And the corresponding superposition data D1 and D 2 , analyze and compare the rationality of structural morphology in the spatial plane. Generally, the later collected data uses more advanced equipment and technology, the data quality is higher, and the structure is more reasonable. As the objective function of correction, the weight coefficient is 0. For the boundary of a certain collection of data, the coverage area is not satisfactory, the structural reliability is low, and a complete correction is required. The weight coefficient is given 1, and other positions are determined and interpolated according to the data quality. That is, the spatial weighting coefficient of the time difference is calculated, and the total time difference of each CMP stacking data is calculated in combination with the relative time difference, which is recorded as the correction amount of the CMP stacking data.

[0043] S4, CMP correction amount T 1 , T 2 Decompose to the shot check point. Check the original acquisition parameters and compare and analyze the inversion velocity model (M) corresponding to the two acquisition data. 1 and M 2 ), find out the main reason for the time difference between the two acquisitions. If it is caused by the correction of the shot point depth, then the correction amount T of the CMP stacking data is1 and T 2 Applied to the shot point, if the cause is the near-surface velocity, the correction amount of the CMP stacking data is allocated according to the weathering layer travel time of the shot check point, and the corrected static correction amount is obtained by combining the static correction amount. For example, the weathering layer travel time of the shot check point is 12ms and 24ms respectively, and the CMP correction amount is 6ms. Then the correction amount allocated to the shot point is 6*(12 / (12+24))=2ms, and the detection point is 4ms. After the decomposition is completed, the correction time difference of the two acquisition shot points is T 1s and T 1r The time differences of the two acquisition detection points are T 2s and T 2r Static correction S for inversion calculation 1s , S 1r , S 2s and S 2r Apply time difference T 1s , T 1r , T 2s and T 2r , and obtain the corrected static correction S 1s '=S 1s +T 1s , S 1r '=S 1r +T 1r , S 2s '=S 2s +T 2s and S 2r '=S 2r +T 2r .

[0044] S5. Apply static correction, compare and stack images, and perform quality control. Figure 4 As shown in the figure, a is the first acquisition data stacking, b is the second acquisition data stacking, and the static correction value S after correction is applied 1s '、S 1r '、S 2s '、S 2r ', the two collected data are fused and superimposed to form data D 12 , D 12 For Figure 4 c shown, due to D 12 The information of the two data collections is integrated, and the time difference between the two data collections is eliminated. In theory, D 12 The signal-to-noise ratio is higher. If the actual stacking effect is better than the single acquisition stacking D 1 , D 2 ,Right now Figure 4If the signal-to-noise ratio and structural rationality of c are higher than those of a and b, it can be determined that the time difference correction is reasonable and the correction process can be terminated. It is determined whether the effect after superposition and fusion is better than the CMP superposition data. If so, the process is terminated; if not, the process returns to step S3. That is, if D 12 The signal-to-noise ratio is lower than D 1 , D 2 If any one of the above is found, it means that the previous decomposition process is unreasonable, and it is necessary to return to step S3 and adjust the gun inspection weighting coefficient in step S3.

[0045] [Example 2]

[0046] like Figure 5 As shown, this embodiment provides a time difference correction system for repeatedly acquiring seismic data, comprising:

[0047] Static correction calculation unit: used to obtain two acquisition data at the same location, and use different physical points to invert the velocity model and calculate the static correction value based on it; the two acquisition data at the same location are distinguished by different gun station numbers so that they have different physical points.

[0048] Unify the data of the two acquisitions. When the excitation and reception of the two acquisitions are at the same coordinate position, the positions of the shot point and the detection point are the same during the two acquisitions. In order to distinguish the two acquisitions, we add a number of acquisitions before the shot detection station number. For example, for the first acquisition, the original 30015011 shot station number is changed to 130015011, and the second acquisition station number at the same position is changed to 230015011 (you can also modify only one of the acquisition records). After such arrangement, the records of the points at the same physical position acquired twice can be used to invert the velocity model and calculate the static correction using two different physical points (shot points or detection points at the same coordinate position). The two acquisitions are recorded separately, and the same shot station is recorded as a different shot station number in the two records. That is, the two acquisitions at the same position are distinguished by different shot station numbers, so that they have different physical points.

[0049] The near-surface velocity inversion is performed for each data collection. Pick the first arrival, use the first arrival wave tomography technology, use the same tomography method and parameters for each data collection, and perform near-surface velocity inversion for each data collection. The inverted near-surface velocity model is recorded as M 1 、M 2 .like Figure 2 The figure shows the velocity inversion results of near-surface tomography from two acquisitions at the same location. It can be seen that the velocities inverted from the two acquisitions are different, such as Figure 2 As shown, there is a difference in speed between 0.2s and 0.3s.

[0050] A CMP stacking data acquisition unit: performs CMP stacking on the corresponding collected data using the static correction amount to acquire corresponding CMP stacking data;

[0051] According to the result of step 2, the static correction is calculated respectively. The static correction of the two collected shot points is recorded as S 1s , S 2s The static correction of the two acquisition points is denoted as S 1r , S 2r Each time data is collected, a seismic slope surface can be obtained. Each seismic slope surface is composed of several cmps, and each cmp is composed of many shot checks (offset distances). The static correction amount is used to perform CMP stacking on each collected data. The CMP stacking data is recorded as D 1 and D 2 , compare D 1 and D 2 , if D 1 and D 2 If there is a time difference, correction is required.

[0052] A cross-correlation calculation unit: using cross-correlation to calculate the relative time difference between at least two CMP stacking data;

[0053] Use cross-correlation to calculate the time difference. (eg Figure 3 As shown in the figure, the part of the data or graph in the figure is the time difference) calculates the CMP superposition data twice (D 1 With D 2 ), pick up the time corresponding to the maximum amplitude value of the cross-correlation function, and obtain the relative time difference T between the two acquisitions of each CMP (common center point gather) xy (x, y are coordinates, CMP position, unit is meter). When calculating the cross-correlation function, if there is a difference in seismic wavelets between the two acquisitions, it is necessary to apply shaping filtering before cross-correlation to eliminate the influence of the wavelets in advance. In order to improve the accuracy of cross-correlation, the main frequency filter is used to improve D 1 With D 2 A high signal-to-noise ratio is also necessary (preprocessing).

[0054] Calculate the spatial weighting coefficient Q of the time difference xy (x, y are coordinates), and the time difference is reasonably distributed to the two acquisitions. The weight coefficient is expressed as a percentage. If D 1 The weight coefficient is Q xy , corresponding to D 2 1-Q xy , then for each CMP, D 1 With D 2 The total time difference allocated is T xy *Q xy and T xy *(1-Qxy ), denoted as T 1 and T 2 , where T xy =T 1 +T 2 The weight coefficient should be determined according to the data quality. 1 With M 2 And the corresponding superposition data D1 and D 2 , analyze and compare the rationality of structural morphology in the spatial plane. Generally, the later collected data uses more advanced equipment and technology, the data quality is higher, and the structure is more reasonable. As the objective function of correction, the weight coefficient is 0. For the boundary of a certain collection of data, the coverage area is not satisfactory, the structural reliability is low, and a complete correction is required. The weight coefficient is given 1, and other positions are determined and interpolated according to the data quality. That is, the spatial weighting coefficient of the time difference is calculated, and the total time difference of each CMP stacking data is calculated in combination with the relative time difference, which is recorded as the correction amount of the CMP stacking data.

[0055] Static correction unit: applying the relative time difference decomposition to the static correction amount in step S1 to obtain a corrected static correction amount;

[0056] Calculate the cross-correlation function between at least two CMP stacking data; pick up the time corresponding to the maximum amplitude value of the cross-correlation function to obtain each spatial weighting coefficient, and calculate the relative time difference between the two acquisitions of each CMP in combination with the relative time difference; calculate the total time difference of the time difference P stacking data, and record it as the correction amount of the CMP stacking data.

[0057] A superposition and fusion unit is used to superimpose and fuse the at least two acquisition data using the corrected static correction amount.

[0058] Apply static correction, compare and stack images, and perform quality control. Figure 4 As shown in the figure, a is the first acquisition data stacking, b is the second acquisition data stacking, and the static correction value S after correction is applied 1s '、S 1r '、S 2s '、S 2r ', the two collected data are fused and superimposed to form data D 12 , D 12 For Figure 4 c shown, due to D 12 The information of the two data collections is integrated, and the time difference between the two data collections is eliminated. In theory, D 12 The signal-to-noise ratio is higher. If the actual stacking effect is better than the single acquisition stacking D 1 , D 2 ,Right now Figure 4If the signal-to-noise ratio and structural rationality of c are higher than those of a and b, it can be determined that the time difference correction is reasonable and the correction process can be terminated. It is determined whether the effect after superposition and fusion is better than the CMP superposition data. If so, the process is terminated; if not, the process returns to step S3. That is, if D 12 The signal-to-noise ratio is lower than D 1 , D 2 If any one of the above is found, it means that the previous decomposition process is unreasonable, and it is necessary to return to step S3 and adjust the gun inspection weighting coefficient in step S3.

[0059] The present invention adopts the method and principle of static correction calculation based on the near-surface velocity model obtained by tomographic inversion, starts from the static correction calculation of the shot check point, and corrects the time difference problem of repeated data collection from the source.

[0060] This method first re-arranges the data collected twice in a unified manner, and uses each collected data to perform near-surface velocity inversion and static correction calculation, and superimposes the superimposed data body collected each time. The cross-correlation of the superimposed data collected twice is calculated, and the maximum amplitude time of the cross-correlation function is picked up to obtain the space-varying time difference. Finally, according to the near-surface velocity distribution, the time difference is reasonably distributed and applied to the collected shot inspection static correction. Based on the near-surface velocity model obtained by tomographic inversion, the present invention corrects the time difference problem between the repeatedly collected data from the source of the shot inspection point according to the method and principle of static correction calculation. It is suitable for the pre-processing of pre-stack migration imaging, which can improve the consistency of the successively collected data, and eliminate the problems such as false structures caused by the time difference problem, thereby improving the in-phase superposition of seismic data, improving the effect of data processing results, and effectively improving the seismic imaging effect. The present invention has good application effects in the fields of time-shift seismic processing, repeated acquisition fusion processing, and multi-block contiguous processing.

Claims

1. A method for fusion processing of repeatedly acquired seismic data, characterized in that: The steps include: S1, obtain at least two acquisition data at the same location, and use different physical points to invert the velocity model and calculate the static correction amount based on it; S2, using the static correction amount to perform CMP stacking on the corresponding collected data to obtain corresponding CMP stacking data; S3, calculating the relative time difference between at least two CMP stacking data by cross-correlation; S4, decomposing the relative time difference into two acquisitions, and further decomposing it into shot inspection points, and applying it to the static correction amount described in step S1 to obtain a corrected static correction amount; S5, using the corrected static correction amount to perform superposition and fusion on the at least two acquisition data; The step S3 specifically includes: calculating a cross-correlation function between at least two of the CMP stack data; Picking up the time corresponding to the maximum amplitude value of the cross-correlation function to obtain the relative time difference between two acquisitions of each CMP; Calculating a spatial weighting coefficient of the time difference, and calculating a total time difference of each CMP stacking data in combination with the relative time difference, which is recorded as a correction amount of each CMP stacking data; The step S4 specifically includes: Comparing the inversion velocity models corresponding to the at least two acquisition data, determining whether the cause of the acquisition time difference is the shot point well depth correction or the near-surface velocity; If the reason is the correction of the well depth of the shot point, the correction amount of the CMP stacking data is directly decomposed to the corresponding shot point; If it is caused by the near-surface velocity, the correction amount of the CMP stacking data is distributed according to the weathering layer travel time of the shot inspection point, and combined with the static correction amount to obtain the corrected static correction amount.

2. The method for fusion processing of repeatedly acquired seismic data according to claim 1, characterized in that In step S1, at least two acquisition data of the same position are distinguished by different gun station numbers so that they have different physical points.

3. The method for fusion processing of repeatedly acquired seismic data according to claim 1, characterized in that: After step S2 and before step S3, the method further comprises: Compare at least two CMP stacking data and make corrections if there is a time difference.

4. The method for fusion processing of repeatedly acquired seismic data according to claim 1, characterized in that: After step S5, the method further comprises: Determine whether the effect after superposition and fusion is better than the CMP superposition data, if so, end the process; If not, return to step S3.

5. A system for fusion processing of repeatedly acquired seismic data, used to implement the method for fusion processing of repeatedly acquired seismic data as claimed in any one of claims 1 to 4, comprising: Static correction calculation unit: used to obtain at least two acquisition data at the same location, and invert the velocity model with different physical points and calculate the static correction accordingly; A CMP stacking data acquisition unit: performs CMP stacking on the corresponding collected data using the static correction amount to acquire corresponding CMP stacking data; A cross-correlation calculation unit: using cross-correlation to calculate the relative time difference between at least two CMP stacking data; A static correction value correction unit: applying the relative time difference decomposition to the static correction value in the static correction value calculation unit to obtain a corrected static correction value; A superposition and fusion unit is used to superimpose and fuse the at least two acquisition data using the corrected static correction amount.

6. The system for repeated acquisition of seismic data fusion processing according to claim 5, characterized in that: In the static correction calculation unit, at least two acquisition data at the same position are distinguished by different gun station numbers so that they have different physical points.

7. The system for repeated acquisition of seismic data fusion processing according to claim 5, characterized in that: The cross-correlation calculation unit is used for: calculating a cross-correlation function between at least two of the CMP stack data; Picking up the time corresponding to the maximum amplitude value of the cross-correlation function to obtain the relative time difference between two acquisitions of each CMP; The spatial weighting coefficient of the time difference is calculated, and the total time difference of each CMP stacking data is calculated in combination with the relative time difference, which is recorded as the correction amount of the CMP stacking data.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores at least one computer-executable program, and when the at least one program is executed by the computer, the computer executes the steps in the method for fusion processing of repeatedly acquired seismic data as described in any one of claims 1-4.

Images