Pre-stack gather event leveling method and device for mobile mapping

Abstract

The embodiment of the invention relates to a pre-stack gather event leveling method and device for mobile mapping. The method comprises the following steps: acquiring an original pre-stack gather to be subjected to event leveling; calculating a first leveling amount of adjacent traces in the original pre-stack gather; generating a local leveling gather based on the first leveling amount and the original pre-stack gather; calculating a second leveling amount of adjacent traces in the local leveling trace set; on the basis of a preset rule, moving mapping leveling is carried out on the second leveling amount; and calculating a time difference correction value after leveling and preset well-seismic matching calibration to complete pre-stack gather event leveling. In this way, the gather data is subjected to leveling processing and event leveling, second-class AVO features are reserved, and meanwhile it is ensured that no stretching distortion exists; the method has an accurate and good leveling effect for the gather with the second-class AVO characteristics, and ensures that the energy of the leveled gather is focused, and the correction value is reasonable and distortionless.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of geophysical processing technology, and in particular to a method and device for leveling event axes of pre-stack gathers of mobile mapping. Background Art

[0002] In seismic data processing, stacking is an effective denoising method. However, in actual seismic data processing, the dynamic correction velocity is based on the NMO algorithm with horizontal layer assumptions, and it is the interpolation velocity of discrete points, which is difficult to calculate accurately; when there are non-uniform geological bodies in the seismic ray path, even the time migration algorithm cannot obtain the leveled gathers. Therefore, there will be a certain amount of residual dynamic correction on the data after dynamic correction, which will reduce the profile resolution when doing in-phase stacking.

[0003] Whether it is the common reflection point stacking in the seismic processing stage or the AVA analysis in the seismic data interpretation stage, they are all based on the common reflection point gathers after event leveling. Uneven events will affect the quality of subsequent processing and even get wrong results. In processing, event leveling (dynamic correction) is often performed based on the time-distance curve equation of the common reflection point gathers. Although dynamic correction can ensure the basic leveling of the events of the common reflection point gathers, dynamic correction often fails to achieve the leveling of the common reflection point gathers due to various factors (such as inaccurate static correction due to the drastic fluctuation of the surface elevation, and the conventional time-distance curve equation error due to the assumption of horizontal layered isotropic media), which reduces the signal-to-noise ratio of the common reflection point stacking and leads to inaccurate or even wrong results of the elastic parameters obtained by AVA analysis.

[0004] At present, the premise of performing gather leveling is to obtain accurate time shift of pre-stack gathers. Many residual dynamic correction methods have been proposed, such as the earlier coherent and high-precision stacking, radial filtering technology, and the technology of treating residual dynamic correction as local random static time difference. Some researchers have also proposed to solve the problem of residual dynamic correction through residual velocity analysis, such as the time shift method of extracting the time difference of sampling points at different times. By dividing the model trace and the common reflection point gather into regions and opening time windows for cross-correlation calculation, the time shift of the significant extreme point of the seismic trace can be determined more accurately, and other data can be stretched or contracted and then moved and leveled as a whole.

[0005] The main channel gather leveling algorithms currently have the following problems: the leveling algorithm needs to predetermine a reference channel, and then compare the other channels in the channel gather with the reference channel. When the phase axis displacement of the two compared seismic channels is large (more than one wavelength), it will lead to erroneous leveling results; most leveling algorithms cannot solve the Class II AVO problem of reverse polarity of the phase axis; the current channel gather leveling algorithm has poor lateral stability and may be over-leveled. Summary of the invention

[0006] In view of this, in order to solve the above technical problems or part of the technical problems, an embodiment of the present invention provides a method and device for leveling event axes of pre-stack gathers with mobile mapping.

[0007] In a first aspect, an embodiment of the present invention provides a method for leveling events of pre-stack gathers of mobile mapping, comprising:

[0008] Obtaining the original pre-stack gathers to be subjected to event leveling;

[0009] Calculate the first leveling amount of adjacent traces in the original pre-stack trace gather;

[0010] Generate a local leveling gather based on the first leveling amount and the original pre-stack gather;

[0011] Calculate the second leveling amount of the adjacent track in the local leveling track set;

[0012] Based on a preset rule, performing mobile mapping leveling on the second leveling amount;

[0013] The pre-stack gather event leveling is completed by calculating the time difference correction after leveling and the preset well-seismic matching calibration.

[0014] In a possible implementation, the method further includes:

[0015] Acquire an original pre-stack gather of a time-domain common imaging point to be subjected to event leveling, wherein the original pre-stack gather has a residual dynamic correction amount;

[0016] Based on the residual normal error, a range of the residual normal error of the original pre-stack gather is calculated from the data header.

[0017] In a possible implementation, the method further includes:

[0018] By comparing and analyzing the similarity of the waveforms of adjacent seismic traces in the original pre-stack gather, a first leveling amount of adjacent traces in the original pre-stack gather of the time domain common imaging point to be leveled for the event is calculated;

[0019] The length of the sliding time window is controlled by the number of points where the amplitude value of the seismic trace waveform in the original pre-stack gather passes through zero.

[0020] In a possible implementation, the method further includes:

[0021] Applying the first leveling amount to the original pre-stack gathers through a first formula to generate a local leveling gather, wherein the first formula is: I1(t, θ)=d(t+dt1, θ);

[0022] Among them, I1(t, θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling value.

[0023] In a possible implementation, the method further includes:

[0024] Based on the similarity of the waveforms of adjacent seismic traces in the local leveling gather, the second leveling amount of the adjacent traces in the local leveling gather is calculated.

[0025] In a possible implementation, the method further includes:

[0026] By using the second formula, a near-track superposition track is selected to perform mobile mapping leveling on the second leveling amount, and the second formula is: I2(t, θ)=d(t+dt2, θ);

[0027] Wherein, I2(t, θ) represents the gather after mobile mapping leveling, and θ represents the incident angle.

[0028] In a possible implementation, the method further includes:

[0029] Based on the gathers after mobile mapping leveling, calculate the time difference correction after leveling;

[0030] The well-seismic matching consistency is calibrated by the well-seismic matching calibration method;

[0031] The pre-stack gather event leveling is completed based on the post-leveling time difference correction and the preset well-seismic matching calibration.

[0032] In a second aspect, an embodiment of the present invention provides a pre-stack gather event leveling device for mobile mapping, comprising:

[0033] An acquisition module, used for acquiring original pre-stack gathers to be subjected to event leveling;

[0034] A calculation module, used for calculating the first leveling amount of adjacent traces in the original pre-stack trace gather;

[0035] A generating module, configured to generate a local leveling gather based on the first leveling amount and the original pre-stack gather;

[0036] The calculation module is further used to calculate the second leveling amount of the adjacent track in the local leveling track set;

[0037] A leveling module, configured to perform mobile mapping leveling on the second leveling amount based on a preset rule;

[0038] The calculation module is also used to calculate the time difference correction amount after leveling and the preset well-seismic matching calibration to complete the pre-stack gather event leveling.

[0039] In a third aspect, an embodiment of the present invention provides an electronic device, comprising: a processor and a memory, wherein the processor is used to execute a mobile mapped pre-stack trace gather event leveling program stored in the memory to implement the mobile mapped pre-stack trace gather event leveling method described in the first aspect above.

[0040] In a fourth aspect, an embodiment of the present invention provides a storage medium, comprising: the storage medium stores one or more programs, and the one or more programs can be executed by one or more processors to implement the pre-stack trace gather in-phase axis leveling method of mobile mapping described in the first aspect above.

[0041] The mobile mapping pre-stack gather event leveling scheme provided in the embodiment of the present invention obtains the original pre-stack gather to be leveled; calculates the first leveling amount of adjacent tracks in the original pre-stack gather; generates a local leveling gather based on the first leveling amount and the original pre-stack gather; calculates the second leveling amount of adjacent tracks in the local leveling gather; performs mobile mapping leveling on the second leveling amount based on preset rules; calculates the time difference correction amount after leveling and the preset well-seismic matching calibration to complete the event leveling of the pre-stack gather. Compared with the current mainstream gather leveling algorithm, the wrong leveling results appear, the problem of the second type of AVO with reverse polarity of the event cannot be solved well, and there is the problem of over-leveling. According to this scheme, the gather data is leveled, the event is leveled, the second type of AVO features are retained, and at the same time, no stretching distortion is ensured; for the gathers with the second type of AVO features, accurate and good leveling effect is achieved, and the energy of the gathers after leveling is ensured to be focused, and the correction amount is reasonable and distortion-free. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] Figure 1 A schematic flow chart of a method for leveling event axes of pre-stack gathers of mobile mapping provided in an embodiment of the present invention;

[0043] Figure 2 A schematic flow chart of another method for leveling events of pre-stack gathers using mobile mapping provided by an embodiment of the present invention;

[0044] Figure 3 A schematic diagram of a pre-stack gather with uneven events provided by an embodiment of the present invention;

[0045] Figure 4 A schematic diagram of a pre-stack gather after mobile mapping event leveling provided by an embodiment of the present invention;

[0046] Figure 5 A schematic diagram of residual dynamic correction with uneven event axis provided by an embodiment of the present invention;

[0047] Figure 6A schematic diagram of the residual dynamic correction amount after the moving mapping event axis is leveled provided in an embodiment of the present invention;

[0048] Figure 7 A schematic structural diagram of a mobile mapping pre-stack gather event leveling device provided in an embodiment of the present invention;

[0049] Figure 8 A schematic diagram of the structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0050] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0051] To facilitate understanding of the embodiments of the present invention, specific embodiments will be further explained below in conjunction with the accompanying drawings. The embodiments do not constitute a limitation on the embodiments of the present invention.

[0052] Embodiment 1:

[0053] Figure 1 A flow chart of a method for leveling event axes of pre-stack gathers of mobile mapping provided by an embodiment of the present invention is shown in FIG. Figure 1 As shown, the method specifically includes:

[0054] S11. Obtain original pre-stack gathers to be subjected to event leveling.

[0055] In the embodiment of the present invention, the original pre-stack gathers of the time domain common imaging points that need to be leveled for the event axis can be input by the user end, such as Figure 3 As shown. The input original pre-stack gathers to be subjected to event leveling have a certain amount of residual dynamic correction. The approximate range of the residual dynamic correction of the original pre-stack gathers is calculated from the data header, and the range is used to judge the effectiveness of the method.

[0056] S12, calculating the first leveling amount of adjacent traces in the original pre-stack trace gather.

[0057] By comparing and analyzing the similarity of the waveforms of adjacent seismic traces in the original pre-stack trace gather, the first leveling amount dt1 of adjacent traces in the original pre-stack trace gather of the time domain common imaging point to be leveled is calculated. For example, the leveling amount of the first and third traces can be determined by adding the leveling amount of the first and second traces and the leveling amount of the second and third traces.

[0058] Furthermore, the length of the sliding time window is controlled by the number of points where the amplitude of the seismic trace waveform in the original pre-stack gather passes through zero.

[0059] S13. Generate a local leveling gather based on the first leveling amount and the original pre-stack gather.

[0060] The first leveling value dt1 is applied to the original pre-stack gather of the time domain common imaging point through the first formula to generate a local leveling gather. The first formula is as follows:

[0061] I1(t,θ)=d(t+dt1,θ)

[0062] Among them, I1(t, θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling value.

[0063] S14, calculating the second leveling amount of the adjacent tracks in the local leveling track set.

[0064] Based on the similarity of the waveforms of adjacent seismic traces in the locally leveled gather obtained above, the second leveling amount dt2 of the adjacent traces in the locally leveled gather is calculated.

[0065] S15. Perform motion mapping leveling on the second leveling amount based on preset rules.

[0066] Through the second formula, select the short-cut "superposition track" and implement mobile mapping leveling for the second leveling amount. The second formula is as follows:

[0067] I2(t1θ)=d(t+dt2,θ)

[0068] Wherein, I2(t, θ) represents the gather after mobile mapping leveling, and θ represents the angle.

[0069] S16, calculating the time difference correction amount after leveling and the preset well-seismic matching calibration to complete the pre-stack gather event leveling.

[0070] Calculate the time difference correction of the pre-stack gather after leveling to ensure that the energy of the gather after leveling is focused and the correction is reasonable without distortion. Use the well-seismic matching calibration method to calibrate the well-seismic matching to ensure the leveling effect of the pre-stack gather. For the original pre-stack gather of the common imaging point in the time domain, use the moving mapping leveling method to obtain the pre-stack gather of the common imaging point in the time domain after leveling, such as Figure 4 As shown. Figure 4 It can be seen that the event axis is effectively leveled, retaining the second type AVO characteristics and maintaining the relative amplitude relationship.

[0071] Figure 5 is the residual NMO due to unevenness of the event axis; Figure 6 is the residual NMO after applying the pre-stack gather event leveling method using motion mapping. Figure 6 It can be seen that the residual dynamic correction amount ±20ms in the pre-stack gather is eliminated, and the energy is more focused. A pre-stack gather event leveling method with mobile mapping provided by an embodiment of the present invention can effectively eliminate the residual dynamic correction amount, level the energy focus of the gather after leveling, and the correction amount is reasonable and without distortion.

[0072] The mobile mapping pre-stack gather event leveling method provided in the embodiment of the present invention obtains the original pre-stack gather to be leveled; calculates the first leveling amount of adjacent tracks in the original pre-stack gather; generates a local leveling gather based on the first leveling amount and the original pre-stack gather; calculates the second leveling amount of adjacent tracks in the local leveling gather; performs mobile mapping leveling on the second leveling amount based on preset rules; calculates the time difference correction amount after leveling and the preset well-seismic matching calibration to complete the pre-stack gather event leveling. Compared with the current mainstream gather leveling algorithm, the wrong leveling results appear, the problem of the second type of AVO with reverse polarity of the event cannot be solved well, and there is the problem of over-leveling. By this method, the gather data is leveled, the event is leveled, the second type of AVO features are retained, and at the same time, no stretching distortion is ensured; this method has an accurate and good leveling effect for gathers with the second type of AVO features, ensuring that the energy of the gather is focused after leveling, and the correction amount is reasonable and distortion-free.

[0073] Embodiment 2:

[0074] Figure 2 A schematic flow chart of another method for leveling the event axis of a pre-stack gather using mobile mapping provided by an embodiment of the present invention is shown in FIG. Figure 2 As shown, the method specifically includes:

[0075] S21, obtaining original pre-stack gathers of time-domain common imaging points to be subjected to event leveling, wherein the original pre-stack gathers have residual dynamic correction amounts.

[0076] In the embodiment of the present invention, the original pre-stack gathers of the time domain common imaging points that need to be leveled for the event axis can be input by the user end, such as Figure 3 As shown in Figure 2, the input original pre-stack gathers to be leveled for event axis have a certain amount of residual dynamic correction.

[0077] S22. Based on the residual dynamic correction amount, calculate the range of the residual dynamic correction amount of the original pre-stack gather from the data header.

[0078] The approximate range of the residual dynamic correction amount of the original pre-stack gather is calculated from the data header, and the range is used to judge the effectiveness of the method.

[0079] S23, calculating the first leveling amount of adjacent traces in the original pre-stack gather of the time domain common imaging point to be leveled for the event by comparing and analyzing the similarity of the waveforms of adjacent seismic traces in the original pre-stack gather.

[0080] By comparing and analyzing the similarity of the waveforms of adjacent seismic traces in the original pre-stack trace set, the first leveling amount dt1 of adjacent traces in the original pre-stack trace set of the time domain common imaging point to be leveled is calculated. For example, the leveling amount of the first trace and the third trace can be obtained by adding the leveling amount of the first trace and the second trace and the leveling amount of the second trace and the third trace.

[0081] S24, controlling the length of the sliding time window by the number of points where the amplitude value of the seismic trace waveform in the original pre-stack gather passes through zero.

[0082] The length of the sliding window is controlled by the number of points where the amplitude of the seismic trace waveform in the original gather passes through zero, which can ensure that the shallow, medium and deep layers of the seismic traces in the gather have a suitable sliding window length.

[0083] S25. Apply the first leveling amount to the original pre-stack gather using a first formula to generate a local leveling gather.

[0084] The first leveling amount is applied to the original pre-stack gathers through the first formula to generate local leveling gathers. The first formula is as follows:

[0085] I1(t,θ)=d(t+dt1,θ)

[0086] Among them, I1(t, θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling value.

[0087] Specifically, any sliding time window is selected in the lateral direction of the original pre-stack gather; the middle track of the sliding time window is used as the reference track, and the first leveling amount is applied to the slide time window to level the gather; the gather leveling result is placed on the middle track of the sliding time window to obtain a locally leveled gather.

[0088] S26. Based on the similarity of the waveforms of adjacent seismic traces in the local leveling gather, calculate the second leveling amount of adjacent traces in the local leveling gather.

[0089] The second leveling amount dt2 of the adjacent track in the local leveling gather is calculated again using the local leveling gather obtained above.

[0090] S27. Select a short-distance superposition track by using the second formula, and perform mobile mapping leveling on the second leveling amount.

[0091] Through the second formula, select the short-cut "superposition track" and implement mobile mapping leveling for the second leveling amount. The second formula is as follows:

[0092] I2(t,θ)=d(t+dt2,θ)

[0093] Wherein, I2(t, θ) represents the gather after mobile mapping leveling, and θ represents the incident angle.

[0094] S28. Calculate the time difference correction amount after leveling based on the track gather after mobile mapping leveling.

[0095] S29. Calibrate the well-seismic matching consistency by using the well-seismic matching calibration method.

[0096] S210, completing pre-stack gather event leveling based on the post-leveling time difference correction amount and the preset well-seismic matching calibration.

[0097] The following is a unified description of S28-S210:

[0098] Calculate the time difference correction of the pre-stack gather after leveling to ensure that the energy of the gather after leveling is focused and the correction is reasonable without distortion. Use the well-seismic matching calibration method to calibrate the well-seismic matching to ensure the leveling effect of the pre-stack gather. For the original pre-stack gather of the common imaging point in the time domain, use the moving mapping leveling method to obtain the pre-stack gather of the common imaging point in the time domain after leveling, such as Figure 4 As shown. Figure 4 It can be seen that the event axis is effectively leveled, retaining the second type AVO characteristics and maintaining the relative amplitude relationship.

[0099] Figure 5 is the residual NMO due to unevenness of the event axis; Figure 6 is the residual NMO after applying the pre-stack gather event leveling method using motion mapping. Figure 6 It can be seen that the residual dynamic correction amount ±20ms in the pre-stack gather is eliminated, and the energy is more focused. A pre-stack gather event leveling method with mobile mapping provided by an embodiment of the present invention can effectively eliminate the residual dynamic correction amount, level the energy focus of the gather after leveling, and the correction amount is reasonable and without distortion.

[0100] The mobile mapping pre-stack trace gather event leveling method provided in an embodiment of the present invention is as follows: obtaining the original pre-stack trace gather to be leveled; calculating the first leveling amount of adjacent traces in the original pre-stack trace gather; generating a local leveling trace gather based on the first leveling amount and the original pre-stack trace gather; calculating the second leveling amount of adjacent traces in the local leveling trace gather; performing mobile mapping leveling on the second leveling amount based on a preset rule; calculating the time difference correction amount after leveling and the preset well-seismic matching calibration to complete the pre-stack trace gather event leveling. By this method, the trace gather data is leveled, the event is flattened, the second type of AVO features are retained, and at the same time, no stretching distortion is ensured; this method has an accurate and good leveling effect for the trace gather with the second type of AVO features, ensuring that the energy of the trace gather is focused after leveling, and the correction amount is reasonable and distortion-free.

[0101] Embodiment three:

[0102] Figure 7 A schematic structural diagram of a mobile mapping pre-stack gather event leveling device provided in an embodiment of the present invention specifically includes:

[0103] An acquisition module 701 is used to acquire original pre-stack gathers to be subjected to event leveling;

[0104] A calculation module 702, configured to calculate a first leveling amount of adjacent traces in the original pre-stack trace gather;

[0105] A generating module 703, configured to generate a local leveling gather based on the first leveling amount and the original pre-stack gather;

[0106] The calculation module 702 is further used to calculate the second leveling amount of the adjacent track in the local leveling track set;

[0107] A leveling module 704, configured to perform movement mapping leveling on the second leveling amount based on a preset rule;

[0108] The calculation module 702 is also used to calculate the post-leveling time difference correction amount and the preset well-seismic matching calibration to complete the pre-stack gather event leveling.

[0109] In a possible implementation, the acquisition module 701 is also used to acquire the original pre-stack trace set of the time domain common imaging point to be subjected to in-phase axis leveling, wherein the original pre-stack trace set has residual dynamic correction amount; based on the residual dynamic correction amount, the range of the residual dynamic correction amount of the original pre-stack trace set is calculated from the data header.

[0110] In one possible implementation, the calculation module 702 is also used to calculate the first leveling amount of adjacent traces in the original pre-stack trace set of the time domain common imaging point to be leveled for the in-phase axis by comparing and analyzing the similarity of the waveforms of adjacent seismic traces in the original pre-stack trace set; and to control the length of the sliding time window by the number of points where the amplitude value of the seismic trace waveform in the original pre-stack trace set passes through zero.

[0111] In a possible implementation, the generation module 703 is also used to apply the first leveling amount to the original pre-stack gather through a first formula to generate a local leveling gather, and the first formula is: I1(t, θ) = d(t+dt1, θ); wherein I1(t, θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling amount.

[0112] In a possible implementation, the calculation module 702 is further configured to calculate a second leveling amount of adjacent traces in the local leveling gather based on similarity of waveforms of adjacent seismic traces in the local leveling gather.

[0113] In a possible implementation, the leveling module 704 is further used to select a near-track superposition track through a second formula to perform mobile mapping leveling on the second leveling amount, and the second formula is: I2(t, θ) = d(t+dt2, θ); wherein I2(t, θ) represents the track set after mobile mapping leveling, and θ represents the incident angle.

[0114] In a possible implementation, the calculation module 702 is also used to calculate the time difference correction after leveling based on the track set after the mobile mapping leveling; calibrate the well-seismic matching consistency through the well-seismic matching calibration method; and complete the pre-stack track set event axis leveling based on the time difference correction after leveling and the preset well-seismic matching calibration.

[0115] The pre-stack gather event leveling device for mobile mapping provided in this embodiment can be as follows: Figure 7 The pre-stack gather event leveling device for mobile mapping shown in FIG. can perform the following operations: Figure 1-2 All steps of the pre-stack gather event leveling method based on mobile mapping are achieved Figure 1-2 The technical effect of the pre-stack gather event leveling method for mobile mapping is shown in Figure 2. Figure 1-2 For the sake of brevity, the relevant description is not repeated here.

[0116] Figure 8 A schematic diagram of the structure of an electronic device provided by an embodiment of the present invention, Figure 8The electronic device 800 shown includes: at least one processor 801, a memory 802, at least one network interface 804 and other user interfaces 803. The various components in the electronic device 800 are coupled together via a bus system 805. It is understood that the bus system 805 is used to achieve connection and communication between these components. In addition to the data bus, the bus system 805 also includes a power bus, a control bus and a status signal bus. However, for the sake of clarity, the bus system 805 is not described in detail. Figure 8 Various buses are labeled as bus system 805 .

[0117] The user interface 803 may include a display, a keyboard, or a pointing device (eg, a mouse, a trackball, a touch pad, or a touch screen).

[0118] It can be understood that the memory 802 in the embodiment of the present invention can be a volatile memory or a non-volatile memory, or can include both volatile and non-volatile memories. Among them, the non-volatile memory can be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory can be a random access memory (RAM), which is used as an external cache. By way of example but not limitation, many forms of RAM are available, such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDRSDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link dynamic random access memory (SLDRAM), and direct RAM bus random access memory (DRRAM). The memory 802 described herein is intended to include, but is not limited to, these and any other suitable types of memory.

[0119] In some implementations, the memory 802 stores the following elements, executable units or data structures, or a subset thereof, or an extended set thereof: an operating system 8021 and an application program 8022 .

[0120] The operating system 8021 includes various system programs, such as a framework layer, a core library layer, a driver layer, etc., which are used to implement various basic services and process hardware-based tasks. The application program 8022 includes various application programs, such as a media player (Media Player), a browser (Browser), etc., which are used to implement various application services. The program for implementing the method of the embodiment of the present invention can be included in the application program 8022.

[0121] In the embodiment of the present invention, by calling the program or instruction stored in the memory 802, specifically, the program or instruction stored in the application 8022, the processor 801 is used to execute the method steps provided by each method embodiment, for example, including:

[0122] The invention relates to obtaining an original pre-stack trace set for in-phase axis leveling; calculating a first leveling amount of adjacent traces in the original pre-stack trace set; generating a local leveling trace set based on the first leveling amount and the original pre-stack trace set; calculating a second leveling amount of adjacent traces in the local leveling trace set; performing mobile mapping leveling on the second leveling amount based on a preset rule; and calculating a time difference correction amount after leveling and a preset well-seismic matching calibration to complete in-phase axis leveling of the pre-stack trace set.

[0123] In one possible implementation, an original pre-stack trace gather of a time-domain common imaging point to be subjected to event leveling is obtained, wherein the original pre-stack trace gather has a residual dynamic correction amount; and based on the residual dynamic correction amount, a range of the residual dynamic correction amount of the original pre-stack trace gather is calculated from a data trace header.

[0124] In one possible implementation, by comparing and analyzing the similarities of the waveforms of adjacent seismic traces in the original pre-stack trace set, the first leveling amount of adjacent traces in the original pre-stack trace set of the time domain common imaging point to be leveled with the event axis is calculated; and the length of the sliding time window is controlled by the number of points where the amplitude values ​​of the seismic trace waveforms in the original pre-stack trace set pass through zero.

[0125] In one possible implementation, the first leveling amount is applied to the original pre-stack gather through a first formula to generate a local leveling gather, and the first formula is: I1(t, θ)d(t+dt1, θ); wherein I1(t, θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling amount.

[0126] In a possible implementation, based on the similarity of waveforms of adjacent seismic traces in the local leveled trace gather, the second leveling amount of adjacent traces in the local leveled trace gather is calculated.

[0127] In one possible implementation, a near-track superposition track is selected through a second formula, and mobile mapping leveling is performed on the second leveling amount. The second formula is: I2(t, θ) = d(t+dt2, θ); wherein I2(t, θ) represents the track set after mobile mapping leveling, and θ represents the incident angle.

[0128] In one possible implementation, based on the track gather after mobile mapping leveling, the time difference correction after leveling is calculated; the well-seismic matching consistency is calibrated through the well-seismic matching calibration method; and the pre-stack track gather event axis leveling is completed based on the time difference correction after leveling and the preset well-seismic matching calibration.

[0129] The method disclosed in the above embodiment of the present invention can be applied to the processor 801, or implemented by the processor 801. The processor 801 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the hardware integrated logic circuit in the processor 801 or the instruction in the form of software. The above processor 801 can be a general processor, a digital signal processor (Digital Signal Processor, DSP), an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), a field programmable gate array (Field Programmable Gate Array, FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components. The methods, steps and logic block diagrams disclosed in the embodiments of the present invention can be implemented or executed. The general processor can be a microprocessor or the processor can also be any conventional processor. The steps of the method disclosed in the embodiment of the present invention can be directly embodied as a hardware decoding processor to execute, or the hardware and software units in the decoding processor can be executed. The software unit can be located in a mature storage medium in the field such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory 802, and the processor 801 reads the information in the memory 802 and completes the steps of the above method in combination with its hardware.

[0130] It is understood that the embodiments described herein may be implemented in hardware, software, firmware, middleware, microcode, or a combination thereof. For hardware implementation, the processing unit may be implemented in one or more application specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPDevice, DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), general purpose processors, controllers, microcontrollers, microprocessors, other electronic units for performing the functions described in the present application, or a combination thereof.

[0131] For software implementation, the technology described herein can be implemented by a unit that performs the functions described herein. The software code can be stored in a memory and executed by a processor. The memory can be implemented in the processor or outside the processor.

[0132] The electronic device provided in this embodiment may be Figure 8 The electronic device shown in FIG. 1 may perform the following steps: Figure 1-2 All steps of the pre-stack gather event leveling method based on mobile mapping are realized Figure 1-2 The technical effect of the pre-stack gather event leveling method for mobile mapping is shown in Figure 2. Figure 1-2 For the sake of brevity, the relevant description is not repeated here.

[0133] The embodiment of the present invention also provides a storage medium (computer-readable storage medium). The storage medium here stores one or more programs. The storage medium may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a read-only memory, a flash memory, a hard disk or a solid-state drive; the memory may also include a combination of the above-mentioned types of memory.

[0134] When one or more programs in the storage medium can be executed by one or more processors, the above-mentioned method for leveling the event axis of the pre-stack gather of mobile mapping executed on the electronic device side is implemented.

[0135] The processor is used to execute the mobile mapped pre-stack gather event leveling program stored in the memory to implement the following steps of the mobile mapped pre-stack gather event leveling method executed on the electronic device side:

[0136] The invention relates to obtaining an original pre-stack trace set for in-phase axis leveling; calculating a first leveling amount of adjacent traces in the original pre-stack trace set; generating a local leveling trace set based on the first leveling amount and the original pre-stack trace set; calculating a second leveling amount of adjacent traces in the local leveling trace set; performing mobile mapping leveling on the second leveling amount based on a preset rule; and calculating a time difference correction amount after leveling and a preset well-seismic matching calibration to complete in-phase axis leveling of the pre-stack trace set.

[0137] In one possible implementation, an original pre-stack trace gather of a time-domain common imaging point to be subjected to event leveling is obtained, wherein the original pre-stack trace gather has a residual dynamic correction amount; and based on the residual dynamic correction amount, a range of the residual dynamic correction amount of the original pre-stack trace gather is calculated from a data trace header.

[0138] In one possible implementation, by comparing and analyzing the similarities of the waveforms of adjacent seismic traces in the original pre-stack trace set, the first leveling amount of adjacent traces in the original pre-stack trace set of the time domain common imaging point to be leveled with the event axis is calculated; and the length of the sliding time window is controlled by the number of points where the amplitude values ​​of the seismic trace waveforms in the original pre-stack trace set pass through zero.

[0139] In one possible implementation, the first leveling amount is applied to the original pre-stack gather through a first formula to generate a local leveling gather, and the first formula is: I1(t, θ)d(t+dt1, θ); wherein I1(t, θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling amount.

[0140] In a possible implementation, based on the similarity of waveforms of adjacent seismic traces in the local leveled trace gather, the second leveling amount of adjacent traces in the local leveled trace gather is calculated.

[0141] In one possible implementation, a near-track superposition track is selected through a second formula, and mobile mapping leveling is performed on the second leveling amount. The second formula is: I2(t, θ) = d(t+dt2, θ); wherein I2(t, θ) represents the track set after mobile mapping leveling, and θ represents the incident angle.

[0142] In one possible implementation, based on the track gather after mobile mapping leveling, the time difference correction after leveling is calculated; the well-seismic matching consistency is calibrated through the well-seismic matching calibration method; and the pre-stack track gather event axis leveling is completed based on the time difference correction after leveling and the preset well-seismic matching calibration.

[0143] The professionals should further realize that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented by electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in the above description according to the function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0144] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0145] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is 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 leveling event axes of pre-stack gathers by mobile mapping, characterized in that: include: Obtaining the original pre-stack gathers to be subjected to event leveling; Calculate the first leveling amount of adjacent traces in the original pre-stack trace gather; Generate a local leveling gather based on the first leveling amount and the original pre-stack gather; Calculate the second leveling amount of the adjacent track in the local leveling track set; Based on a preset rule, performing mobile mapping leveling on the second leveling amount; The pre-stack gather event leveling is completed by calculating the time difference correction after leveling and the preset well-seismic matching calibration.

2. The method according to claim 1, characterized in that The method of obtaining the original pre-stack gathers to be subjected to event leveling comprises: Acquire an original pre-stack gather of a time-domain common imaging point to be subjected to event leveling, wherein the original pre-stack gather has a residual dynamic correction amount; Based on the residual normal error, a range of the residual normal error of the original pre-stack gather is calculated from the data header.

3. The method according to claim 1, characterized in that The calculating of the first leveling amount of adjacent traces in the original pre-stack trace gather comprises: By comparing and analyzing the similarity of the waveforms of adjacent seismic traces in the original pre-stack gather, a first leveling amount of adjacent traces in the original pre-stack gather of the time domain common imaging point to be leveled for the event is calculated; The length of the sliding time window is controlled by the number of points where the amplitude value of the seismic trace waveform in the original pre-stack gather passes through zero.

4. The method according to claim 3, characterized in that The generating of a local leveling gather based on the first leveling amount and the original pre-stack gather comprises: Applying the first leveling amount to the original pre-stack gathers through a first formula to generate a local leveling gather, wherein the first formula is: I1(t,θ)=d(t+dt1,θ); Among them, I1(t,θ) represents the local leveling gather, θ represents the incident angle, and dt1 is the first leveling value.

5. The method according to claim 4, characterized in that The calculating of the second leveling amount of the adjacent track in the local leveling track set comprises: Based on the similarity of the waveforms of adjacent seismic traces in the local leveling gather, the second leveling amount of the adjacent traces in the local leveling gather is calculated.

6. The method according to claim 5, characterized in that The performing mobile mapping leveling on the second leveling amount based on a preset rule includes: By using the second formula, a short-distance superposition track is selected to perform mobile mapping leveling on the second leveling amount, and the second formula is: I2(t,θ)=d(t+dt2,θ); Among them, I2(t,θ) represents the gather after mobile mapping leveling, and θ represents the incident angle.

7. The method according to claim 6, characterized in that The calculation of the post-leveling time difference correction amount and the preset well-seismic matching calibration completes the pre-stack gather event leveling, including: Based on the gathers after mobile mapping leveling, calculate the time difference correction after leveling; The well-seismic matching consistency is calibrated by the well-seismic matching calibration method; The pre-stack gather event leveling is completed based on the post-leveling time difference correction and the preset well-seismic matching calibration.

8. A mobile mapping pre-stack gather event leveling device, characterized in that: include: An acquisition module, used for acquiring original pre-stack gathers to be subjected to event leveling; A calculation module, used for calculating the first leveling amount of adjacent traces in the original pre-stack trace gather; A generating module, configured to generate a local leveling gather based on the first leveling amount and the original pre-stack gather; The calculation module is further used to calculate the second leveling amount of the adjacent track in the local leveling track set; A leveling module, configured to perform mobile mapping leveling on the second leveling amount based on a preset rule; The calculation module is also used to calculate the time difference correction amount after leveling and the preset well-seismic matching calibration to complete the pre-stack gather event leveling.

9. An electronic device, characterized in that: include: A processor and a memory, wherein the processor is used to execute a mobile mapped pre-stack gather event leveling program stored in the memory to implement the mobile mapped pre-stack gather event leveling method according to any one of claims 1 to 7.

10. A storage medium, characterized in that: The storage medium stores one or more programs, and the one or more programs can be executed by one or more processors to implement the pre-stack gather event leveling method for mobile mapping according to any one of claims 1 to 7.

Images