Method and device for processing fault interruption of seismic data imaging

Abstract

The present invention discloses a method and device for processing fault interruption of seismic data imaging, the method comprising: calculating the imaging value and illumination times of each seismic trace at each imaging point underground according to the seismic recording data, the shot point position and the receiving point position of each seismic trace in the seismic data; sequentially superimposing the imaging value and illumination times of each seismic trace at each imaging point underground in the order of the number of seismic traces from small to large; outputting the stacking data of the seismic traces currently stacked when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces; continuously executing the above steps until each seismic trace in the seismic data participates in the stacking to obtain multiple seismic trace stacking data; performing stacking calculation on the multiple seismic trace stacking data to obtain a seismic imaging data body. The present invention can solve the problem of resource consumption caused by the need to re-image when the seismic data imaging fault is interrupted, and improve the timeliness of imaging of massive seismic data.

Description

Technical Field

[0001] The invention relates to the technical field of geophysical oil and gas exploration, and in particular to a method and device for processing fault interruption of seismic data imaging. Background Art

[0002] This section is intended to provide a background or context to the embodiments of the invention recited in the claims. No description herein is admitted to be prior art by inclusion in this section.

[0003] Seismic exploration is a method that uses the differences in elasticity and density of underground media to observe and analyze the response of the earth to artificially excited seismic waves to infer the underground structural morphology and rock properties. It has become the most important means of oil and natural gas resource exploration and is also widely used in coalfield and engineering geological surveys, regional geological research, and crustal research. At present, oil and gas resource exploration mainly uses reflection wave seismic exploration, which consists of three parts: seismic data acquisition, processing, and interpretation. Among them, seismic data migration imaging is the key core link of seismic data processing.

[0004] With the development of seismic data acquisition technology, wide-azimuth, high-density 3D seismic acquisition technology with small bins, small sampling time intervals, and large arrays has become the mainstream. Each shot records tens of thousands of traces, and the amount of data in a work area is often as high as several T (1T = 10 12 ). Faced with such a huge amount of big data, 3D seismic imaging requires a lot of calculations and takes a long time. For a long time, geophysicists have been committed to developing seismic data migration imaging, mainly Kirchhoff depth migration, wave equation depth migration and reverse time migration. Currently, the most commonly used in industry is Kirchhoff depth migration, which can realize reflection wave imaging, scattered wave imaging and anisotropic imaging. Kirchhoff 3D depth migration imaging consists of two parts: travel time calculation and integral imaging.

[0005] The calculation of the three-dimensional integral imaging method of seismic big data can be roughly divided into two modes: "multiple input and single output" and "single input and multiple output". Among them, the "single input and multiple output" mode is more widely used, avoiding the time-consuming problem of multiple readings of massive seismic big data in the "multiple input and single output" mode. The "single input and multiple output" mode reads one seismic trace at a time and distributes it to each point in the three-dimensional imaging space according to the Kirchhoff migration principle. In the specific implementation process, the "single input and multiple output" mode inputs a seismic record each time and saves the imaging value in the three-dimensional array until all seismic record traces have completed the migration imaging, and then outputs the imaging value in the memory to the disk file. The purpose of this is to improve the calculation efficiency, because memory operations are much faster than disk operations, especially for angle domain imaging data volumes, whose data volume is dozens of times that of conventional imaging data, avoiding the disk storage problem in the in-and-out process.

[0006] In order to improve computing efficiency, people have developed a variety of parallel processing solutions, such as internal parallel imaging calculation and parallel gather input. However, these methods all have a risk of interruption due to faults, such as power outages, machine hardware failures, and system failures (commonly known as "crashes"). Since the three-dimensional imaging of massive seismic data takes a long time, and the imaging results are retained in the memory during the migration imaging process, once the machine is restarted, the previous calculation results will be lost. Recalculation will greatly waste machine time and cause a huge resource consumption problem. Summary of the invention

[0007] The embodiment of the present invention provides a method for processing seismic data imaging fault interruption, which is used to solve the problem of resource consumption caused by re-imaging when seismic data imaging fault interruption occurs. The method includes:

[0008] According to the seismic record data, shot point position and receiving point position of each seismic channel in the seismic data, the imaging value and illumination times of each imaging point of each seismic channel in the underground are calculated;

[0009] The imaging value and illumination times of each seismic trace at each imaging point underground are sequentially stacked in the order of the number of seismic traces from small to large. If the stacking process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and illumination times of each seismic trace at each imaging point underground are sequentially stacked again in the order of the number of seismic traces from small to large.

[0010] When the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the stacking data of the currently stacked seismic traces is output; the output imaging seismic trace number is used to represent the number of seismic traces that can be stacked within a preset time length; the seismic trace stacking data is used to represent the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at each imaging point underground;

[0011] Continue to perform the above steps until each seismic trace in the seismic data is involved in the stacking, and obtain multiple seismic trace stacking data;

[0012] The stacking calculation is performed on multiple seismic trace stacking data to obtain a seismic imaging data volume.

[0013] The embodiment of the present invention further provides a processing device for seismic data imaging fault interruption, which is used to solve the problem of resource consumption caused by re-imaging when seismic data imaging fault interruption occurs. The device includes:

[0014] The imaging value and illumination times calculation module is used to calculate the imaging value and illumination times of each imaging point of each seismic channel in the underground according to the seismic record data, the shot point position and the receiving point position of each seismic channel in the seismic data;

[0015] The stacking module is used to stack the imaging value and illumination times of each seismic trace at each imaging point underground in order of the number of seismic traces from small to large. If the stacking process is interrupted due to a fault, the maximum number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and illumination times of each seismic trace at each imaging point underground are stacked again in order of the number of seismic traces from small to large.

[0016] The seismic trace stacking data output module is used to: output the stacking data of the seismic traces currently stacked when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces; the output imaging seismic trace number is used to represent the number of seismic traces that can be stacked within a preset time length; the seismic trace stacking data is used to represent the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at each imaging point underground;

[0017] A continuous execution module is used to: continuously execute the above steps until each seismic trace in the seismic data is stacked to obtain multiple seismic trace stacking data;

[0018] The seismic imaging data volume calculation module is used to: perform superposition calculation on multiple seismic trace superposition data to obtain a seismic imaging data volume.

[0019] An embodiment of the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for processing seismic data imaging fault interruption when executing the computer program.

[0020] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program for executing the above-mentioned method for processing seismic data imaging fault interruption.

[0021] In an embodiment of the present invention, the imaging value and illumination times of each seismic trace at each underground imaging point are calculated according to the seismic recording data, the shot point position and the receiving point position of each seismic trace in the seismic data; the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the number of seismic traces from small to large; wherein, if the superposition process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace superposition data is obtained, and starting from the number of seismic traces after the maximum value, the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed again in the order of the number of seismic traces from small to large; when the number of seismic traces participating in the superposition reaches an integer multiple of the preset output imaging seismic trace number, the currently superimposed seismic trace superposition data is output; the output imaging seismic trace number is used to characterize the number of seismic traces that can be superimposed within a preset time length The seismic trace stacking data is used to characterize the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at various imaging points underground; the above steps are continuously performed until each seismic trace in the seismic data participates in the stacking to obtain multiple seismic trace stacking data; the multiple seismic trace stacking data are stacked and calculated to obtain a seismic imaging data body, so that when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the currently stacked seismic trace stacking data is output, and the output seismic trace stacking data can be controlled in stages. When the seismic data imaging is interrupted by a fault, the stacked seismic trace stacking data can be protected. Compared with the prior art, the situation in which the completed seismic trace stacking data is completely lost due to the interruption of the seismic data imaging fault is avoided, the problem of resource consumption caused by the need to re-image is solved, and the timeliness of the imaging of massive seismic data is improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0023] Figure 1 A schematic flow chart of a method for processing a seismic data imaging fault interruption according to an embodiment of the present invention;

[0024] Figure 2 It is a flow chart of an example of a method for processing a seismic data imaging fault interruption in an embodiment of the present invention;

[0025] Figure 3 A specific flow chart of a method for processing a seismic data imaging fault interruption in an embodiment of the present invention;

[0026] Figure 4 A common center point gather map input in an embodiment of the present invention;

[0027] Figure 5 is a longitudinal cross-sectional view of a reflection wave imaging body in an embodiment of the present invention;

[0028] Figure 6 is a transverse cross-sectional view of a reflection wave imaging body in an embodiment of the present invention;

[0029] Figure 7 It is a longitudinal cross-sectional view of an imaging body before being interrupted in an embodiment of the present invention;

[0030] Figure 8 It is a transverse cross-sectional view of an imaging body before the improvement in an embodiment of the present invention;

[0031] Fig. 9 It is a longitudinal cross-sectional view of an improved interrupted imaging body in an embodiment of the present invention;

[0032] Fig.10 is a transverse cross-sectional view of an improved interrupted imaging body in an embodiment of the present invention;

[0033] Fig.11 It is a structural schematic diagram of a device for processing fault interruption of seismic data imaging in an embodiment of the present invention;

[0034] Fig.12 This is a specific example diagram of a device for processing seismic data imaging fault interruption in an embodiment of the present invention. DETAILED DESCRIPTION

[0035] To make the purpose, technical solution and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below in conjunction with the accompanying drawings. Here, the exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0036] At present, wide-azimuth, high-density three-dimensional seismic data acquisition has acquired a massive amount of seismic big data, which has brought great challenges to big data imaging processing. In order to improve computing efficiency, people have developed a variety of parallel processing schemes, such as internal parallel imaging calculations and parallel track set input. However, these methods all have a risk of fault interruption, such as power outages, machine hardware failures, and system failures (commonly known as "freezes"). Since the three-dimensional imaging of massive large-scale seismic data takes a long time, and the imaging results are retained in the memory during the offset imaging process, once the machine is restarted, the previous calculation results will disappear. Recalculation will greatly waste machine time and cause great resource consumption problems. The embodiment of the present invention provides a method for processing seismic data imaging fault interruptions, which provides a set of solutions for seismic big data integral imaging interruptions, can avoid the risks of computer software and hardware accidents, and ensure the timeliness of seismic massive data imaging while ensuring uniform imaging energy distribution. The method for processing seismic data imaging fault interruptions, such as Figure 1 As shown, the following steps may be included:

[0037] Step 101: Calculate the imaging value and illumination times of each imaging point in the underground of each seismic channel according to the seismic record data, shot point position and receiving point position of each seismic channel in the seismic data;

[0038] Step 102: In the order of the number of seismic traces from small to large, the imaging value and the illumination times of each seismic trace at each imaging point underground are sequentially stacked respectively; wherein, if the stacking process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and the illumination times of each seismic trace at each imaging point underground are sequentially stacked again in the order of the number of seismic traces from small to large;

[0039] Step 103: when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, output the stacking data of the seismic traces currently stacked; the output imaging seismic trace number is used to represent the number of seismic traces that can be stacked within a preset time length; the seismic trace stacking data is used to represent the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at each imaging point underground;

[0040] Step 104: Continue to perform the above steps until each seismic trace in the seismic data is stacked to obtain multiple seismic trace stacking data;

[0041] Step 105: Perform stacking calculation on multiple seismic trace stacking data to obtain a seismic imaging data volume.

[0042] In an embodiment of the present invention, the imaging value and illumination times of each seismic trace at each underground imaging point are calculated according to the seismic recording data, the shot point position and the receiving point position of each seismic trace in the seismic data; the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the trace number of the seismic trace from small to large; wherein, if the superposition process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace superposition data is obtained, and starting from the number of the seismic trace after the maximum value, the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the trace number of the seismic trace from small to large again; when the number of each seismic trace involved in the superposition reaches an integer multiple of the preset output imaging seismic trace number, the currently superimposed seismic trace superposition data is output; the output imaging seismic trace number is used to represent the number of seismic traces that can be superimposed within a preset time length ; The seismic trace stacking data is used to characterize the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at various imaging points underground; the above steps are continuously performed until each seismic trace in the seismic data participates in the stacking to obtain multiple seismic trace stacking data; the multiple seismic trace stacking data are stacked and calculated to obtain a seismic imaging data body, so that when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the currently stacked seismic trace stacking data is output, and the output seismic trace stacking data can be controlled in stages. When the seismic data imaging is interrupted by a fault, the stacked seismic trace stacking data can be protected. Compared with the prior art, the situation in which the completed seismic trace stacking data is completely lost due to the interruption of the seismic data imaging fault is avoided, the problem of resource consumption caused by the need to re-image is solved, and the timeliness of the imaging of massive seismic data is improved.

[0043] In specific implementation, firstly, the imaging value and illumination times of each underground imaging point of each seismic channel are calculated according to the seismic record data, shot point position and receiving point position of each seismic channel in the seismic data.

[0044] In an embodiment, the method for processing seismic data imaging fault interruption provided by the embodiment of the present invention may also include: obtaining the travel time, ray azimuth, ray polar angle and Jacobian determinant element value from the shot point and receiving point of each seismic channel in the seismic data to each underground imaging point; calculating the imaging value and illumination times of each seismic channel at each underground imaging point according to the seismic recording data, shot point position and receiving point position of each seismic channel in the seismic data, which may include: calculating the imaging value and illumination times of each seismic channel at each underground imaging point according to the seismic recording data, shot point position, receiving point position of each seismic channel in the seismic data, the travel time, ray azimuth, ray polar angle and Jacobian determinant element value from the shot point and receiving point to each underground imaging point.

[0045] In the above embodiment, by calculating the imaging value and illumination times of each seismic trace at each imaging point underground, it can be ensured that the imaging energy distribution of the seismic imaging data volume obtained in the subsequent steps is uniform.

[0046] In specific implementation, after calculating the imaging value and illumination times of each seismic channel at each underground imaging point according to the seismic record data, shot point position and receiving point position of each seismic channel in the seismic data, the imaging value and illumination times of each seismic channel at each underground imaging point are sequentially superimposed in the order of the seismic channel number from small to large; wherein, if the superposition process is interrupted due to a fault, the maximum value of the seismic channel number in the output seismic channel superposition data is obtained, and starting from the seismic channel number after the maximum value, the imaging value and illumination times of each seismic channel at each underground imaging point are sequentially superimposed in the order of the seismic channel number from small to large again.

[0047] In the above embodiment, by sequentially superimposing the imaging value and illumination times of each seismic trace at each imaging point underground, it is helpful to perform superimposition calculation on multiple seismic trace superposition data in subsequent steps to obtain a seismic imaging data body. If the superposition process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace superposition data is obtained, and the number of seismic traces after the maximum value is started, and the imaging value and illumination times of each seismic trace at each imaging point underground are sequentially superimposed in the order of the number of seismic traces from small to large. It can be ensured that when the superposition process is interrupted due to a fault, the currently superimposed seismic trace superposition data can be output, and by starting with the number of seismic traces after the maximum value, the number of seismic traces is restarted in the order of the number of seismic traces from small to large, the breakpoint continuation of seismic trace superposition can be achieved, and the output seismic trace superposition data can be controlled in stages. When the seismic data imaging is interrupted due to a fault, the superimposed seismic trace superposition data can be protected, and the failure risk caused by power failure or machine failure in the middle of seismic big data imaging can be avoided, and the resource consumption problem caused by the need to re-image can be solved.

[0048] In a specific implementation, after the imaging value and the number of illuminations of each seismic channel at each underground imaging point are sequentially superimposed in the order of the number of seismic channels from small to large, when the number of seismic channels participating in the superposition reaches an integer multiple of the preset number of output imaging seismic channels, the currently superimposed seismic channel superposition data is output; the output imaging seismic channel number is used to characterize the number of seismic channels that can be superimposed within a preset time length; the seismic channel superposition data is used to characterize the sum of the imaging values ​​and the sum of the illumination times of multiple seismic channels at each underground imaging point; the above steps are continuously performed until each seismic channel in the seismic data participates in the superposition, and multiple seismic channel superposition data are obtained.

[0049] In the embodiment, when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the stacking data of the currently stacked seismic traces is output. In order to solve the problem of machine time loss caused by power outage or machine failure in the three-dimensional integral imaging of massive seismic big data, the stacking data of the currently stacked seismic traces can be output in stages, so that it is not necessary to restart the stacking of seismic trace data from the beginning when the stacking process is interrupted due to a failure, thus avoiding the situation where all the completed seismic trace stacking data are lost due to the interruption of seismic data imaging failure, solving the problem of resource consumption caused by the need to re-image, and improving the timeliness of massive seismic data imaging.

[0050] In the above embodiment, under the premise of ensuring the imaging quality, the three-dimensional integral imaging data volume and the illumination frequency volume can be output in stages and synthesized later. This can not only ensure the uniform distribution of imaging energy, but also avoid the risk of failure caused by power outage or machine failure in the middle of seismic big data imaging.

[0051] In a specific implementation, after obtaining a plurality of seismic trace superposition data, superposition calculation is performed on the plurality of seismic trace superposition data to obtain a seismic imaging data volume.

[0052] In an embodiment, performing superposition calculation on multiple seismic trace superposition data to obtain a seismic imaging data body may include: performing superposition calculation on multiple seismic trace superposition data to respectively obtain the sum of imaging values ​​and the sum of illumination times of all seismic traces in the seismic data at each imaging point underground; and calculating the seismic imaging data body based on the sum of imaging values ​​and the sum of illumination times of all seismic traces in the seismic data at each imaging point underground.

[0053] In the above embodiment, the stacking calculation of multiple seismic trace stacking data is performed according to the following formula to obtain the sum of imaging values ​​and the sum of illumination times of all seismic traces in the seismic data at each imaging point underground:

[0054]

[0055]

[0056] Among them, IMG (ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI (ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; img m (ix, iy, iz) is the imaging value of each underground imaging point of the seismic trace stacking data output for the mth time; ni m(ix, iy, iz) is the illumination times of each underground imaging point of the m-th output seismic trace stacking data; m is the number of output seismic trace stacking data; x, y, z are the spatial coordinates of the underground imaging points respectively; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively.

[0057] In the above embodiment, the seismic imaging data volume can be calculated according to the sum of the imaging values ​​of all seismic traces in the seismic data obtained respectively at each imaging point underground and the sum of the illumination times according to the following formula:

[0058]

[0059] (ix=0,1,…,nx-1; iy=0,1,…,ny-1; iz=0,1,…,nz-1)

[0060] Among them, (ix, iy, iz) is the seismic imaging data volume; IMG (ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI (ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively; nx, ny and nz are the total number of discrete points of the underground imaging points on the x, y and z coordinate axes respectively.

[0061] In the above embodiment, a seismic imaging data volume can be obtained by performing superposition calculation on multiple seismic trace superposition data, thereby achieving the purpose of generating a seismic imaging data volume.

[0062] A specific example is given below to illustrate the specific application of the method of the present invention. Figure 2 As shown, Figure 2 This is a flow chart of an example of a method for processing a seismic data imaging fault interruption in an embodiment of the present invention. This specific example may include the following three steps:

[0063] Step 201, online parallel imaging: according to the seismic record data, shot point position and receiving point position of each seismic channel in the seismic data, the imaging value and illumination times of each seismic channel at each imaging point underground are calculated; the imaging value and illumination times of each seismic channel at each imaging point underground are sequentially stacked in the order of the number of seismic channels from small to large; when the number of seismic channels participating in the stacking reaches an integer multiple of the preset number of output imaging seismic channels, the stacking data of the seismic channels currently stacked are output; the above steps are continuously performed until each seismic channel in the seismic data participates in the stacking, and a plurality of seismic channel stacking data are obtained;

[0064] Step 202, interruption of subsequent processing: if the stacking process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and illumination times of each seismic trace at each imaging point underground are sequentially stacked starting from the number of the seismic trace after the maximum value, and the number of the seismic traces is sequentially increased from small to large;

[0065] Step 203, post-imaging synthesis: performing stacking calculation on multiple seismic trace stacking data to obtain a seismic imaging data volume.

[0066] The following describes the steps in this example in detail in combination with specific application scenarios:

[0067] 1. Online parallel imaging

[0068] (a) Input a seismic record data u(it), (it = 0, 1, ..., nt-1), where nt is the number of sampling points of the record. At the same time, input the shot point position S of the record k and receiving point location G k , where k is the seismic channel number;

[0069] (b) Read the shot point position S from the pre-calculated travel time table k and receiving point location G k Information such as travel time, ray azimuth, ray polar angle, Jacobian determinant elements, etc. to each underground imaging point position M(ix, iy, iz) (ix = 0, 1, ..., nx-1; iy = 0, 1, ..., ny-1; iz = 0, 1, ..., nz-1), where nx, ny, nz are the total number of discrete points in the x, y, and z directions of the three-dimensional imaging area, respectively;

[0070] (c) Based on the theory of reflection wave imaging, scattered wave imaging or anisotropic imaging, the imaging value im91 of each underground imaging point is calculated in parallel k (ix, iy, iz) and lighting times ni1 k (ix, iy, iz), where k is the seismic trace number;

[0071] (d) Online superposition imaging value img1 k (ix, iy, iz) and lighting times ni1 k (ix, iy, iz), and obtain the imaging value of each underground imaging point img(ix, iy, iz) = ∑ k img1 k (ix, iy, iz) and the number of illumination times ni(ix, iy, iz) = ∑ k ni1 k (ix, iy, iz). This completes the imaging process of the current input channel;

[0072] (e) Output intermediate results for interrupt processing:

[0073] (e.1) Determine whether the current input channel number is an integer multiple of the preset parameter OIT (i.e., the preset output imaging seismic channel number), where OIT can be determined by the user based on computer performance and the amount of seismic data imaging calculations, and is generally set to 10,000 channels. The principle of selecting this parameter is the number of channels that can be calculated in one day;

[0074] (e.2) If the number of input channels is not an integer multiple of OIT, skip the intermediate result output shown in the following steps (e.3)-(e.5) and continue with the following step (f);

[0075] (e.3) If the number of input channels is an integer multiple of OIT, the imaging value img (ix, iy, iz) and the number of illumination times ni (ix, iy, iz) are output from the memory array to the disk file. The file name can be specified in advance. In order to save disk space, the same file name is used for each output. This output is completed using a separate output module outputimg(), in which the output imaging value can be processed, such as applying imaging factors, channel header processing, etc. However, the value of the array img (ix, iy, iz) must not be changed because it will continue to be used for the imaging of the next input channel;

[0076] (e.4) After outputting the imaging value and the number of illuminations, the current number of input channels needs to be recorded in a log file, where the log file name can be specified by the user in advance;

[0077] (e.5) Use fflush(logfp) to clear the data in the buffer, where logfp is the stream pointing to the log file;

[0078] (f) Determine whether the current input channel is the last channel. If not, read the next seismic record and continue with the above steps (a)-(e); if it is the last channel, proceed to the next step;

[0079] (g) Call the output module outputimg() to output the imaging value img(ix, iy, iz) and the number of illumination times ni(ix, iy, iz), and finally close the imaging output file and the log file;

[0080] 2. Interrupt follow-up processing

[0081] If there is no interruption during the operation of the online parallel imaging job, skip this step and continue to step 3. When the online parallel imaging job suddenly fails due to a power outage, the imaging job has been running for many days. The usual practice in the prior art is to resend the imaging job, so that all previous calculations will be invalidated. In order to solve this problem, the present invention example takes the following steps:

[0082] (h) Read the last recorded input channel number from the log file to know the number of seismic records completed in the previous round of imaging operation and its imaging value img (ix, iy, iz) and the number of illumination times ni (ix, iy, iz). For the convenience of explanation, they are respectively denoted as img m (ix, iy, iz) and lighting times ni m (ix, iy, iz), where m is the number of times the imaging job is sent;

[0083] (i) Eliminate the input trace data that has been imaged from the massive seismic big data, use the remaining seismic record data as the input for the next round of imaging operations, and change the file name of the imaging output;

[0084] (j) sending the next round of imaging jobs and repeating the calculation in step (1);

[0085] (k) If there is another interruption in the calculation process of this round of operation, repeat the above steps (a)-(c) until all seismic traces are imaged.

[0086] 3. Post-production imaging synthesis

[0087] (a) Read the imaging value img from the file output by multiple rounds of imaging m (ix, iy, iz) and lighting times ni m (ix, iy, iz), calculate the total imaging value IMG (ix, iy, iz) = ∑ m img m (ix, iy, iz) and total lighting times NI(ix, iy, iz) = ∑ m ni m (ix, iy, iz), where m is the number of times the imaging job is sent.

[0088] (b) Calculate the final seismic imaging data volume:

[0089]

[0090] (ix=0,1,…,nx-1; iy=0,1,…,ny-1; iz=0,1,…,nz-1)

[0091] The embodiment of the present invention also provides a specific implementation process of the above example. Figure 2 Step 201 may include the following steps 301-307, such as Figure 3 As shown:

[0092] Step 301: Input a seismic record data u(it), (it=0, 1, ..., nt-1), where nt is the number of sampling points of the record. At the same time, input the shot point position S of the record. k and receiving point location G k , where k is the seismic channel number. In the embodiment of the present invention, the number of sampling points nt per channel is 3500 and the sampling interval is 2 milliseconds. Figure 3 A common center point gather is shown, which has 430 channels. This embodiment consists of 13 beam lines, each beam line has 30 longitudinal lines, and each longitudinal line has 1000 Figure 3 The common center point gather shown in the figure has a total of 167.7 million seismic records, with a total amount of 2.34T (2.34×10 12 ). This embodiment adopts parallel computing, calculates one beam of lines at a time, implements external parallelism for the longitudinal survey line number, calculates one longitudinal survey line for each operation, and adopts 2 nodes to implement online parallel imaging.

[0093] Step 302: Read the shot point positions S and S from the pre-calculated travel time table. k and receiving point location G k The information of travel time, ray azimuth, ray polar angle, Jacobian determinant elements, etc. to each underground imaging point position M(ix, iy, iz) (ix=0, 1, ..., nx-1; iy=0, 1, ..., ny-1; iz=0, 1, ..., nz-1), wherein nx, ny, nz are the total number of discrete points in the x, y, and z directions of the three-dimensional imaging area, respectively. According to the computer disk space and computer computing power, the travel time table grid in this embodiment is larger than the imaging area grid, so when the travel time and other information are applied, bilinear interpolation is used to obtain the information on the specific imaging point position. The grid size used for the travel time table and other information in the embodiment of the present invention is 50 meters × 100 meters × 20 meters, and the number of grid points is nxv=241, nyv=91, nzv=401. The imaging grid size is 12.5 meters × 25 meters × 5 meters, and the number of imaging grid points is nx=500, ny=60, nz=1401.

[0094] Step 303: According to the theory of reflection wave imaging, scattered wave imaging or anisotropic imaging, the imaging value img1 of each underground imaging point is calculated in parallel. k (ix, iy, iz) and lighting times ni1 k(ix, iy, iz), wherein the imaging range is controlled by parameters such as imaging aperture, anti-aliasing operator, scattering angle range and illumination inclination angle range. In the embodiment of the present invention, the reflection wave imaging aperture parameters are as follows: longitudinal aperture xoffsetmax=6000 meters, yoffsetmax=4500 meters; the anti-aliasing operator adopts the depth domain anti-aliasing operator (patent number ZL 2018 1 0521842.X); the scattering angle range is 40°; and the illumination inclination angle range is ±30°.

[0095] Step 304: Online superposition of imaging value img1 k (ix, iy, iz) and lighting times ni1 k (ix, iy, iz), and obtain the imaging value of each underground imaging point img(ix, iy, iz) = ∑ k img1 k (ix, iy, iz) and the number of illumination times ni(ix, iy, iz) = ∑ k ni1 k (ix, iy, iz). This completes the imaging process of the current input channel.

[0096] Step 305: Output intermediate results for interrupt processing

[0097] (1) Determine whether the current input trace number is an integer multiple of the preset parameter OIT, where OIT is determined by the user based on computer performance and the amount of seismic data imaging calculations. In the embodiment of the present invention, OIT is set to 10,000 traces. The principle of selecting this parameter is the number of traces that can be calculated in one day.

[0098] (2) If the number of input channels is not an integer multiple of OIT, skip the intermediate result output shown in the following steps (3)-(5) and continue to the following step S206.

[0099] (3) If the number of input channels is an integer multiple of OIT, the imaging value img (ix, iy, iz) and the number of illumination times ni (ix, iy, iz) are output from the memory array to a disk file. The file name can be specified in advance. In order to save disk space, the same file name is used for each output. This output is completed using a separate output module outputimg, in which the output imaging value can be processed, such as applying imaging factors, channel header processing, etc. However, the value of the array img (ix, iy, iz) cannot be changed because it will continue to be used for the imaging of the next input channel.

[0100] (4) After outputting the imaging value and the number of illumination times, the current input channel number needs to be recorded in a log file, where the log file name can be specified in advance by the user.

[0101] (5) Use fflush(logfp) to clear the data in the buffer, where logfp is the stream pointing to the log file.

[0102] Step 306, determine whether the current input channel is the last channel, if not, read the next earthquake record, and continue with the above steps S201-S205; if it is the last channel, proceed to the next step.

[0103] Step 307, calling the output module outputimg(), outputting the imaging value img(ix, iy, iz) and the number of illumination times ni(ix, iy, iz), and finally closing the imaging output file and the log file.

[0104] In the above Figure 2 In step 202, if there is no interruption during the operation of the online parallel imaging job, this step 202 is skipped and step 203 is continued. In the embodiment of the present invention, one interruption occurs. When the online parallel imaging job suddenly occurs due to a power outage or other failure, the imaging job has been running for many days. The usual practice is to resend the imaging job, so that all previous calculations will be invalidated. In order to solve this problem, the following steps are taken:

[0105] (a) Read the last recorded input channel number from the log file to know the number of seismic records completed in the previous round of imaging operation and its imaging value img(ix, iy, iz) and the number of illumination times ni(ix, iy, iz). For the convenience of explanation, they are respectively denoted as img m (ix, iy, iz) and lighting times ni m (ix, iy, iz), where m is the number of times the imaging job is sent.

[0106] (b) Eliminate the input trace data that has been imaged from the massive seismic big data, use the remaining seismic record data as the input for the next round of imaging operations, and change the file name of the imaging output.

[0107] (c) Send the next round of online parallel imaging jobs and repeat the calculation of S101.

[0108] (d) If there is another interruption in the calculation process of this round of operation, repeat the above steps (a)-(c) until all seismic traces are imaged.

[0109] Figure 2 The post-imaging synthesis in step 303 may include:

[0110] (a) Read the imaging value img from the file output by multiple rounds of imaging m (ix, iy, iz) and lighting times ni m(ix, iy, iz), calculate the total imaging value IMG (ix, iy, iz) = ∑ m img m (ix, iy, iz) and total lighting times NI(ix, iy, iz) = ∑ m ni m (ix, iy, iz), where m is the number of times the imaging job is sent.

[0111] (b) Calculate the final seismic imaging data volume:

[0112]

[0113] Wherein (ix=0, 1, ..., nx-1; iy=0, 1, ...; ny-1; iz=0, 1, ..., nz-1), in the embodiment of the present invention, nx=500, ny=60, nz=1401.

[0114] Compared with the prior art, in order to solve the problem of machine time loss caused by power outages or machine failures encountered in the three-dimensional integral imaging process of massive seismic big data, while ensuring the imaging quality, the embodiments of the invention can output the currently stacked seismic trace stacking data in stages, and output the three-dimensional integral imaging data volume and the illumination frequency volume and post-synthesis in stages. This can not only ensure uniform imaging energy distribution, but also avoid the risk of failure caused by power outages or machine failures in the middle of seismic big data imaging.

[0115] Figure 4 is a common center point gather map input in the embodiment of the present invention. The gather has 430 channels, each channel has 3500 sampling points, and the time sampling interval is 2 milliseconds. In this embodiment, a beam line consists of 30 longitudinal survey lines, each longitudinal survey line has 1000 sampling points. Figure 4 The common center point gather is shown in the figure. There are 13 such beamlines in the entire work area, with a total of 167.7 million seismic records and a total amount of 2.34T (2.34×10 12 ).

[0116] Figure 5 is a longitudinal cross-sectional view (line number 1450) of a reflection wave imaging body in an embodiment of the present invention, Figure 6 It is a transverse cross-sectional view (CDP1666) of the reflected wave imaging body in an embodiment of the present invention. Figure 5 and Figure 6 The image data volume of a longitudinal survey line number 1480 in this massive seismic big data is displayed. The imaging area ranges from line number 1372 to line number 1490, CDP1502-CDP2500, and depth 3000 meters to 7000 meters. The computer model used is Dell Precision Tower7910, and the CPU model is CPU E5-2660V4@2.00GHz, which has 2 physical CPUs, each with 14 cores, total logical CPU number is 56. 15 parallel operations can be performed at the same time, and the operation of inputting each longitudinal survey line lasts 2 months.

[0117] For example, when calculating for 45 days, the power is cut off due to line maintenance. Before the method for processing seismic data imaging fault interruption provided by the embodiment of the present invention is adopted, the calculation of these 45 days is wasted. If the method provided by the embodiment of the present invention is adopted, Figure 5 and Figure 6 The imaging results shown thus resolve the losses due to power outage.

[0118] Figure 7 is a longitudinal cross-sectional view of an imaging body before being interrupted in an embodiment of the present invention (line number 1450), Figure 8 is a transverse cross-sectional view of an imaging body before improvement in an embodiment of the present invention (CDP1666), Fig. 9 is a longitudinal cross-sectional view of an improved interrupted imaging body in an embodiment of the present invention (line number 1450), Fig.10 It is a transverse cross-sectional view of the imaging body after improvement (CDP1666) in an embodiment of the present invention.

[0119] During the imaging process, parameters such as imaging aperture, anti-aliasing operator, scattering angle range, and illumination inclination angle range control the imaging range. Figure 8 and Fig.10 A longer imaging track appears on the middle right, which is mainly determined by the ±30° illumination inclination parameter range.

[0120] contrast Figure 7 and Fig. 9 , and contrast Figure 8 and Fig.10 It can be seen that the lateral energy distribution on the longitudinal and lateral sections of the imaging body before the improvement of the embodiment of the present invention is uneven, while the lateral energy distribution is uniform after the improvement. Because the imaging factor is applied to the output imaging body before the improvement, while the imaging factor is applied to the output imaging body after the improvement while ensuring that the imaging body in the memory remains unchanged, it shows that this interrupt processing solution can provide guarantee for the industrial application of seismic big data processing.

[0121] In an embodiment of the present invention, the imaging value and illumination times of each seismic trace at each underground imaging point are calculated according to the seismic recording data, the shot point position and the receiving point position of each seismic trace in the seismic data; the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the trace number of the seismic trace from small to large; wherein, if the superposition process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace superposition data is obtained, and starting from the number of the seismic trace after the maximum value, the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the trace number of the seismic trace from small to large again; when the number of each seismic trace involved in the superposition reaches an integer multiple of the preset output imaging seismic trace number, the currently superimposed seismic trace superposition data is output; the output imaging seismic trace number is used to represent the number of seismic traces that can be superimposed within a preset time length ; Seismic trace stacking data is used to characterize the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at each imaging point underground; continue to perform the above steps until each seismic trace in the seismic data participates in the stacking, and obtains multiple seismic trace stacking data; perform stacking calculations on multiple seismic trace stacking data to obtain a seismic imaging data body, so that when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the currently stacked seismic trace stacking data is output, and the output seismic trace stacking data can be controlled in stages. When the seismic data imaging is interrupted by a fault, the stacked seismic trace stacking data can be protected. Compared with the prior art, the situation in which the completed seismic trace stacking data disappears due to the interruption of the seismic data imaging fault is avoided, and the resource consumption problem caused by the need to re-image is solved, while the timeliness of the imaging of massive seismic data is improved. In summary, the embodiments of the present invention can be applied to reflection wave imaging, scattered wave imaging and anisotropic imaging of massive seismic big data.

[0122] The present invention also provides a device for processing seismic data imaging fault interruption, such as the following embodiment. Since the principle of solving the problem by the device is similar to the method for processing seismic data imaging fault interruption, the implementation of the device can refer to the implementation of the method for processing seismic data imaging fault interruption, and the repeated parts will not be repeated.

[0123] An embodiment of the present invention provides a device for processing seismic data imaging fault interruption, such as Fig.11 As shown, it may include:

[0124] The imaging value and illumination times calculation module 01 is used to calculate the imaging value and illumination times of each imaging point of each seismic channel in the underground according to the seismic record data, the shot point position and the receiving point position of each seismic channel in the seismic data;

[0125] The stacking module 02 is used to stack the imaging value and illumination times of each seismic trace at each imaging point underground in order of the number of seismic traces from small to large. If the stacking process is interrupted due to a fault, the maximum number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and illumination times of each seismic trace at each imaging point underground are stacked again in order of the number of seismic traces from small to large.

[0126] The seismic trace stacking data output module 03 is used to: output the stacking data of the seismic traces currently stacked when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces; the output imaging seismic trace number is used to represent the number of seismic traces that can be stacked within a preset time length; the seismic trace stacking data is used to represent the sum of the imaging values ​​and the sum of the illumination times of multiple seismic traces at each imaging point underground;

[0127] The continuous execution module 04 is used to: continuously execute the above steps until each seismic trace in the seismic data is stacked to obtain a plurality of seismic trace stacking data;

[0128] The seismic imaging data volume calculation module 05 is used to perform stacking calculation on multiple seismic trace stacking data to obtain a seismic imaging data volume.

[0129] In one embodiment, Fig.12 As shown, it also includes: a seismic parameter acquisition module 06, which is used to: obtain the travel time, ray azimuth, ray polar angle and Jacobian determinant element value from the shot point and receiving point of each seismic channel in the seismic data to each underground imaging point;

[0130] The imaging value and lighting times calculation module is specifically used for:

[0131] Based on the seismic record data, shot point location, receiving point location, travel time from the shot point and receiving point to each underground imaging point, ray azimuth, ray polar angle and Jacobian determinant element values ​​of each seismic channel in the seismic data, the imaging value and illumination times of each seismic channel at each underground imaging point are calculated.

[0132] In one embodiment, the seismic imaging data volume calculation module is specifically used to:

[0133] Performing superposition calculation on the superposition data of multiple seismic traces, respectively obtaining the sum of imaging values ​​and the sum of illumination times of all seismic traces in the seismic data at each imaging point underground;

[0134] The seismic imaging data volume is calculated according to the sum of the imaging values ​​of all seismic traces in the seismic data obtained respectively at each imaging point underground and the sum of the illumination times.

[0135] In one embodiment, the seismic imaging data volume calculation module is specifically used to:

[0136] In one embodiment, the seismic imaging data volume calculation module is specifically used to:

[0137] Seismic imaging data volume calculation module, specifically used for:

[0138] The stacking calculation of multiple seismic trace stacking data is performed according to the following formula to obtain the sum of the imaging values ​​and the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground:

[0139]

[0140]

[0141] Among them, IMG (ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI (ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; img m (ix, iy, iz) is the imaging value of each underground imaging point of the seismic trace stacking data output for the mth time; ni m (ix, iy, iz) is the illumination times of each underground imaging point of the m-th output seismic trace stacking data; m is the number of output seismic trace stacking data; x, y, z are the spatial coordinates of the underground imaging points respectively; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively.

[0142] In one embodiment, the seismic imaging data volume calculation module is specifically used to:

[0143] The stacking calculation of multiple seismic trace stacking data is performed according to the following formula to obtain the sum of the imaging values ​​and the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground:

[0144]

[0145]

[0146] Among them, IMG (ix, iy, iz) is the sum of the imaging values ​​of each underground imaging point obtained by multiple imaging operations in the seismic trace stacking data; NI (ix, iy, iz) is the sum of the illumination times of each underground imaging point obtained by multiple imaging operations in the seismic trace stacking data; img m (ix, iy, iz) is the imaging value of each underground imaging point obtained by the mth imaging operation; ni m(ix, iy, iz) is the number of illuminations of each underground imaging point obtained in the mth imaging operation; m is the number of times the imaging operation is sent; x, y, z are the spatial coordinates of the underground imaging points respectively; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively.

[0147] In one embodiment, the seismic imaging data volume calculation module is specifically used to:

[0148] The seismic imaging data volume is calculated according to the sum of the imaging values ​​of all seismic traces at each imaging point underground and the sum of the illumination times in the seismic data obtained separately, using the following formula:

[0149]

[0150] (ix=0,1,…,nx-1; iy=0,1,…,ny-1; iz=0,1,…,nz-1)

[0151] Wherein, I(ix, iy, iz) is the seismic imaging data volume; IMG(ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI(ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively; nx, ny and nz are the total number of discrete points of the underground imaging points on the x, y and z coordinate axes respectively.

[0152] An embodiment of the present invention further provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for processing seismic data imaging fault interruption when executing the computer program.

[0153] An embodiment of the present invention further provides a computer-readable storage medium, which stores a computer program for executing the above-mentioned method for processing seismic data imaging fault interruption.

[0154] In an embodiment of the present invention, the imaging value and illumination times of each seismic trace at each underground imaging point are calculated according to the seismic recording data, the shot point position and the receiving point position of each seismic trace in the seismic data; the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the trace number of the seismic trace from small to large; wherein, if the superposition process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace superposition data is obtained, and starting from the number of the seismic trace after the maximum value, the imaging value and illumination times of each seismic trace at each underground imaging point are sequentially superimposed in the order of the trace number of the seismic trace from small to large again; when the number of each seismic trace involved in the superposition reaches an integer multiple of the preset output imaging seismic trace number, the currently superimposed seismic trace superposition data is output; the output imaging seismic trace number is used to represent the number of seismic traces that can be superimposed within a preset time length ; The seismic trace stacking data is used to characterize the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at various imaging points underground; the above steps are continuously performed until each seismic trace in the seismic data participates in the stacking to obtain multiple seismic trace stacking data; the multiple seismic trace stacking data are stacked and calculated to obtain a seismic imaging data body, so that when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the currently stacked seismic trace stacking data is output, and the output seismic trace stacking data can be controlled in stages. When the seismic data imaging is interrupted by a fault, the stacked seismic trace stacking data can be protected. Compared with the prior art, the situation in which the completed seismic trace stacking data is completely lost due to the interruption of the seismic data imaging fault is avoided, the problem of resource consumption caused by the need to re-image is solved, and the timeliness of the imaging of massive seismic data is improved.

[0155] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented 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.

[0156] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0157] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0158] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0159] The above specific embodiments further illustrate the purpose, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above are only specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for processing fault interruption of seismic data imaging, characterized in that: include: According to the seismic record data, shot point position and receiving point position of each seismic channel in the seismic data, the imaging value and illumination times of each imaging point of each seismic channel in the underground are calculated; The imaging value and illumination times of each seismic trace at each imaging point underground are sequentially stacked in the order of the number of seismic traces from small to large. If the stacking process is interrupted due to a fault, the maximum value of the number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and illumination times of each seismic trace at each imaging point underground are sequentially stacked again in the order of the number of seismic traces from small to large. When the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces, the stacking data of the currently stacked seismic traces is output; the output imaging seismic trace number is used to represent the number of seismic traces that can be stacked within a preset time length; the seismic trace stacking data is used to represent the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at each imaging point underground; Continue to perform the above steps until each seismic trace in the seismic data is involved in the stacking, and obtain multiple seismic trace stacking data; The stacking calculation is performed on multiple seismic trace stacking data to obtain a seismic imaging data volume.

2. The method according to claim 1, characterized in that Also includes: Obtain the travel time, ray azimuth, ray polar angle and Jacobian determinant element values ​​from the shot point and receiving point of each seismic channel in the seismic data to each underground imaging point; According to the seismic record data, shot point location and receiving point location of each seismic channel in the seismic data, the imaging value and illumination times of each imaging point of each seismic channel in the underground are calculated, including: Based on the seismic record data, shot point location, receiving point location, travel time from the shot point and receiving point to each underground imaging point, ray azimuth, ray polar angle and Jacobian determinant element values ​​of each seismic channel in the seismic data, the imaging value and illumination times of each seismic channel at each underground imaging point are calculated.

3. The method according to claim 1, characterized in that Perform stack calculation on multiple seismic trace stack data to obtain seismic imaging data volume, including: Performing superposition calculation on the superposition data of multiple seismic traces, respectively obtaining the sum of imaging values ​​and the sum of illumination times of all seismic traces in the seismic data at each imaging point underground; The seismic imaging data volume is calculated according to the sum of the imaging values ​​of all seismic traces in the seismic data obtained respectively at each imaging point underground and the sum of the illumination times.

4. The method according to claim 3, characterized in that The stacking calculation of multiple seismic trace stacking data is performed according to the following formula to obtain the sum of the imaging values ​​and the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground: Among them, IMG (ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI (ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; img m (ix, iy, iz) is the imaging value of each underground imaging point of the seismic trace stacking data output for the mth time; ni m (ix, iy, iz) is the illumination times of each underground imaging point of the m-th output seismic trace stacking data; m is the number of output seismic trace stacking data; x, y, z are the spatial coordinates of the underground imaging points respectively; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively.

5. The method according to claim 3, characterized in that The seismic imaging data volume is calculated according to the sum of the imaging values ​​of all seismic traces at each imaging point underground and the sum of the illumination times in the seismic data obtained separately, using the following formula: Wherein, I(ix, iy, iz) is the seismic imaging data volume; IMG(ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI(ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively; nx, ny and nz are the total number of discrete points of the underground imaging points on the x, y and z coordinate axes respectively.

6. A processing device for seismic data imaging fault interruption, characterized in that: include: The imaging value and illumination times calculation module is used to calculate the imaging value and illumination times of each imaging point of each seismic channel in the underground according to the seismic record data, the shot point position and the receiving point position of each seismic channel in the seismic data; The stacking module is used to stack the imaging value and illumination times of each seismic trace at each imaging point underground in order of the number of seismic traces from small to large. If the stacking process is interrupted due to a fault, the maximum number of seismic traces in the output seismic trace stacking data is obtained, and the imaging value and illumination times of each seismic trace at each imaging point underground are stacked again in order of the number of seismic traces from small to large. The seismic trace stacking data output module is used to: output the stacking data of the seismic traces currently stacked when the number of seismic traces participating in the stacking reaches an integer multiple of the preset number of output imaging seismic traces; the output imaging seismic trace number is used to represent the number of seismic traces that can be stacked within a preset time length; the seismic trace stacking data is used to represent the sum of imaging values ​​and the sum of illumination times of multiple seismic traces at each imaging point underground; A continuous execution module is used to: continuously execute the above steps until each seismic trace in the seismic data is stacked to obtain multiple seismic trace stacking data; The seismic imaging data volume calculation module is used to: perform superposition calculation on multiple seismic trace superposition data to obtain a seismic imaging data volume.

7. The device according to claim 6, characterized in that Also includes: The seismic parameter acquisition module is used to obtain the travel time, ray azimuth, ray polar angle and Jacobian determinant element value from the shot point and receiving point of each seismic channel in the seismic data to each underground imaging point; The imaging value and lighting times calculation module is specifically used for: Based on the seismic record data, shot point location, receiving point location, travel time from the shot point and receiving point to each underground imaging point, ray azimuth, ray polar angle and Jacobian determinant element values ​​of each seismic channel in the seismic data, the imaging value and illumination times of each seismic channel at each underground imaging point are calculated.

8. The device according to claim 6, characterized in that Seismic imaging data volume calculation module, specifically used for: Performing superposition calculation on the superposition data of multiple seismic traces, respectively obtaining the sum of imaging values ​​and the sum of illumination times of all seismic traces in the seismic data at each imaging point underground; The seismic imaging data volume is calculated according to the sum of the imaging values ​​of all seismic traces in the seismic data obtained respectively at each imaging point underground and the sum of the illumination times.

9. The device according to claim 8, characterized in that Seismic imaging data volume calculation module, specifically used for: The stacking calculation of multiple seismic trace stacking data is performed according to the following formula to obtain the sum of the imaging values ​​and the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground: Among them, IMG (ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI (ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; img m (ix, iy, iz) is the imaging value of each underground imaging point of the seismic trace stacking data output for the mth time; ni m (ix, iy, iz) is the illumination times of each underground imaging point of the m-th output seismic trace stacking data; m is the number of output seismic trace stacking data; x, y, z are the spatial coordinates of the underground imaging points respectively; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively.

10. The device according to claim 8, characterized in that Seismic imaging data volume calculation module, specifically used for: The seismic imaging data volume is calculated according to the sum of the imaging values ​​of all seismic traces at each imaging point underground and the sum of the illumination times in the seismic data obtained separately, using the following formula: Wherein, I(ix, iy, iz) is the seismic imaging data volume; IMG(ix, iy, iz) is the sum of the imaging values ​​of all seismic traces in the seismic data at each imaging point underground; NI(ix, iy, iz) is the sum of the illumination times of all seismic traces in the seismic data at each imaging point underground; ix, iy and iz are the numbers of the underground imaging points on the x, y and z coordinate axes respectively; nx, ny and nz are the total number of discrete points of the underground imaging points on the x, y and z coordinate axes respectively.

11. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 5 is implemented.

12. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program for executing the method according to any one of claims 1 to 5.

Images