Correction method, device, equipment and medium for time-shifted seafloor node seismic data

Abstract

The present invention discloses a method, device, equipment, and medium for correcting time-shifted seafloor node seismic data. The method includes: determining a correction range for the time-shifted seafloor node seismic data based on the trace integral response data of the time-shifted seafloor node seismic data, and determining a reference correction within the correction range; if the formation thickness is determined to be less than the tuning thickness, determining a tuning error correction based on the trace integral response data and well data; determining a full-area wellpoint error correction based on the wellpoint layering data and seismic horizon data; determining a total correction based on the reference correction, the tuning error correction, and the full-area wellpoint error correction, and correcting the time-shifted seafloor node seismic data based on the total correction. The technical solution of this method can improve the depth accuracy of time-shifted seafloor node seismic data.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas geophysical exploration, and in particular to a correction method, device, equipment and medium for time-shifted seabed node seismic data. Background Art

[0002] Ocean Bottom Node (OBN) is a key technology in marine seismic exploration. It enables effective detection and imaging of blind areas beyond the reach of conventional seismic exploration. It provides wide-azimuth, multi-azimuth, and multi-offset seismic datasets essential for effective seismic imaging of complex and high-velocity shield structures (such as subsalt shields). It also provides high-quality multi-wave data for oil and gas reservoir identification. Time-lapse seismic technology analyzes changes in seismic data response caused by changes in reservoir parameters over a specific time period, thereby studying reservoir dynamics and identifying remaining oil and gas. Time-lapse ocean bottom node seismic exploration is becoming a common exploration method for offshore oil and gas field development.

[0003] With the continuous deepening of offshore oil and gas field development, the potential of remaining oil and the deployment of infill wells have placed increasingly higher requirements on the accuracy of time-lapse seabed node seismic data. There is an urgent need for a correction method for time-lapse seabed node seismic data to improve the success rate of remaining oil prediction and drilling. Summary of the Invention

[0004] The present invention provides a correction method, device, equipment and medium for time-shifted seafloor node seismic data, so as to improve the depth accuracy of time-shifted seafloor node seismic data.

[0005] In a first aspect, an embodiment of the present invention provides a method for correcting time-shifted seafloor node seismic data, the method comprising:

[0006] determining a correction range of the time-lapsed seafloor node seismic data based on the trace integrated response data of the time-lapsed seafloor node seismic data, and determining a reference plane correction amount within the correction range;

[0007] If it is determined that the formation thickness is less than the tuning thickness, determining a tuning error correction amount based on the trace integrated response data and the well data;

[0008] Determine the error correction amount of well points in the entire area based on well point layer data and seismic layer data;

[0009] A total correction amount is determined based on the reference plane correction amount, the tuning error correction amount, and the whole-area well point error correction amount, and the time-lapse seabed node seismic data is corrected based on the total correction amount.

[0010] In a second aspect, an embodiment of the present invention further provides a device for correcting time-shifted seafloor node seismic data, the device comprising:

[0011] a reference plane correction amount determination module, configured to determine a correction range of the time-lapsed seafloor node seismic data based on the trace integral response data of the time-lapsed seafloor node seismic data, and determine a reference plane correction amount within the correction range;

[0012] a tuning error correction amount determination module, configured to determine a tuning error correction amount based on the trace integrated response data and the well data if it is determined that the formation thickness is less than the tuning thickness;

[0013] The module for determining the error correction amount of well points in the whole area is used to determine the error correction amount of well points in the whole area based on the well point layer data and seismic layer data;

[0014] The time-lapse seabed node seismic data correction module is used to determine the total correction amount based on the reference plane correction amount, the tuning error correction amount and the full-area well point error correction amount, and correct the time-lapse seabed node seismic data based on the total correction amount.

[0015] In a third aspect, an embodiment of the present invention further provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, a method for correcting time-shifted seafloor node seismic data as described in any one of the embodiments of the present invention is implemented.

[0016] In a fourth aspect, an embodiment of the present invention further provides a storage medium storing computer-executable instructions, which, when executed by a computer processor, are used to perform a method for correcting time-shifted seafloor node seismic data as described in any one of the embodiments of the present invention.

[0017] The technical solution of the embodiment of the present invention calculates the trace integral response of time-lapsed seafloor node seismic data, determines the correction range of the time-lapsed seafloor node seismic data, and calculates the reference level correction within the correction range. When the formation thickness is less than the tuning thickness, the tuning error correction is determined based on the trace integral response data and well data. The well point error correction for the entire area is determined based on the well point layer data and seismic horizon data. The total correction is determined based on the reference level correction, the tuning error correction, and the well point error correction for the entire area, and the time-lapsed seafloor node seismic data is corrected based on the total correction. This improves the depth accuracy of the time-lapsed seafloor node seismic data, thereby improving the success rate of remaining oil prediction and drilling.

[0018] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it intended to limit the scope of the present invention. Other features of the present invention will become readily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. 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.

[0020] Figure 1 This is a flow chart of a method for correcting time-shifted seafloor node seismic data provided by the first embodiment of the present invention;

[0021] Figure 2 This is a schematic diagram of the maximum amplitude attribute value of time-lapsed seafloor node seismic data provided by the first embodiment of the present invention;

[0022] Figure 3 is a schematic diagram of a correction range provided in Example 1 of the present invention;

[0023] Figure 4 This is a schematic diagram of a zero-phase layer provided by the first embodiment of the present invention;

[0024] Figure 5 is a schematic diagram of a reference plane correction amount provided by the first embodiment of the present invention;

[0025] Figure 6 1 is a schematic diagram of a zero-phase layer after tuning error correction provided by the first embodiment of the present invention;

[0026] Figure 7 is a schematic diagram of a tuning error correction amount provided by the first embodiment of the present invention;

[0027] Figure 8 This is a schematic diagram of a full-area well point error correction amount provided by the first embodiment of the present invention;

[0028] Figure 9 is a schematic diagram of a total error correction amount provided by the first embodiment of the present invention;

[0029] Figure 10 is a schematic diagram of a depth correction model provided in Example 1 of the present invention;

[0030] Figure 11 1 is a schematic diagram of a depth domain profile of time-lapsed seafloor node seismic data after correction provided by the first embodiment of the present invention;

[0031] Figure 12 This is a schematic structural diagram of a correction device for time-shifted seafloor node seismic data provided by a second embodiment of the present invention;

[0032] Figure 13 This is a structural diagram of an electronic device provided in Example 3 of the present invention. DETAILED DESCRIPTION

[0033] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions 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 embodiments described are only 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 making creative efforts should fall within the scope of protection of the present invention.

[0034] It should be noted that the terms "first", "second", etc. in the description and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that the numbers used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0035] Example 1

[0036] Figure 1 A flowchart of a method for correcting time-shifted seabed node seismic data is provided for embodiment 1 of the present invention. This embodiment is applicable to the situation where time-shifted seabed node seismic data is depth-corrected. The method can be executed by a device for correcting time-shifted seabed node seismic data. The device for correcting time-shifted seabed node seismic data can be implemented in the form of hardware and / or software, and the device for correcting time-shifted seabed node seismic data can be configured in an electronic device.

[0037] like Figure 1 As shown, the method includes:

[0038] S110 , determining a correction range of the time-shifted seafloor node seismic data according to the trace integrated response data of the time-shifted seafloor node seismic data, and determining a reference plane correction amount within the correction range.

[0039] Seismic trace integration is an unconstrained inversion technique using seismic data, that is, a technique for performing stratum analysis by obtaining the relative velocity profile of the stratum through integrating the seismic traces.

[0040] In this embodiment, trace integrated response data is calculated for the time-shifted seafloor node seismic data, and a correction range is determined based on the trace integrated response data. Depth correction is performed on the time-shifted seafloor node seismic data within the correction range.

[0041] Furthermore, based on the track-integrated response data of the time-shifted seafloor node seismic data, a correction range of the time-shifted seafloor node seismic data is determined, including: determining the track-integrated response data of the time-shifted seafloor node seismic data, and extracting the maximum amplitude attribute value from the track-integrated response data; determining the maximum amplitude attribute value range that matches the water flooding signal based on the spatial position of the production well and the spatial position of the injection well; and determining the correction range of the time-shifted seafloor node seismic data based on the maximum amplitude attribute value range that matches the water flooding signal.

[0042] In this embodiment, calculating the channel integral response value and extracting the maximum amplitude attribute value from the channel integral response value are both public technologies in the art, and this embodiment does not limit this. Figure 2 A schematic diagram of the maximum amplitude attribute value of the time-lapsed seafloor node seismic data is provided. Figure 2 It can be seen that water can be injected from injection well 1 and pushed to production wells 1, 2, 3, and 5. According to the spatial location of the production wells and injection wells, and the maximum amplitude attribute value corresponding to the node generating the flooding signal, the maximum amplitude attribute value range that needs to be corrected is determined. Figure 2 For example, the maximum amplitude attribute value range corresponding to the node generating the flooding signal is [0.04, 0.12]. Figure 3 A schematic diagram of the correction range is provided, such as Figure 3 As shown, the range within the dotted box is the range of the time-lapse seafloor node seismic data that requires depth correction.

[0043] Furthermore, determining the reference plane correction amount within the correction range includes: determining the zero phase level of the event axis of the channel integral response data within the correction range; and using the difference between the zero phase level and the oil-water interface as the reference plane correction amount.

[0044] Specifically, the zero phase level of the event axis of the channel integrated response data can be determined manually. Figure 2 Take the trace integral lines A-A' and B-B' in the example, Figure 4 A schematic diagram of the zero-phase level is provided, such as Figure 4 As shown, the zero phase plane corresponding to the integrated measurement line A-A' is Figure 4 The gray solid line corresponding to A-A' and the zero phase plane corresponding to the integrated measurement line B-B' are Figure 4 The gray solid line corresponding to B-B' in the figure. The oil-water interface can be determined in advance. This embodiment does not limit the specific value of the oil-water interface. For example, the oil-water interface can be -7244 ft. Figure 4The white dotted line in the figure is the oil-water interface. The difference between the zero-phase plane and the oil-water interface is the reference plane correction value. Figure 4 The difference between the medium grey solid line and the white dotted line represents the reference plane correction amount. The reference plane correction amount within the correction range is as follows: Figure 5 shown.

[0045] S120: If it is determined that the formation thickness is less than the tuning thickness, a tuning error correction amount is determined based on the trace integrated response data and the well data.

[0046] The tuning thickness can be represented by λ / 4. When the layer thickness is less than the tuning thickness, the thinning of the layer no longer manifests as a closer line connecting the peaks and troughs of the reflection event wave, but instead as a weakening of the amplitude. At this point, the zero-phase plane determined from the time-shifted channel integrated response data contains a function error with amplitude as the variable. In this embodiment, a tuning error correction is calculated to eliminate this error.

[0047] Furthermore, a tuning error correction amount is determined based on the trace integrated response data and the well data, including: determining the maximum value of the maximum amplitude attribute value within the correction range; determining the wavelength of the statistical wavelet of the wellside trace based on the well data; and determining the tuning error correction amount based on the maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, the wavelength, and the maximum value of the maximum amplitude attribute value within the correction range.

[0048] Specifically, the maximum value Max of the maximum amplitude attribute value within the statistical correction range is calculated as Figure 2 For example, the Max value is 0.12. Using the well data, the statistical wavelet of the well bypass is extracted. Taking a period of 10ms and an average velocity of 3000m / s as an example, the wavelength λ is calculated to be 100ft.

[0049] The tuning error correction amount is determined according to the following formula:

[0050] Wherein, T represents the tuning error correction amount, Max represents the maximum value of the maximum amplitude attribute value within the correction range, A represents the maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, and λ represents the wavelength of the statistical wavelet of the wellside trace.

[0051] exist Figure 4 On the basis of Figure 6 A schematic diagram of the zero phase plane after tuning error correction is provided, as shown in Figure 6 As shown in the figure, the zero phase layer corresponding to the original gray solid line is represented by the white thick dashed line after the tuning error correction. The schematic diagram of the tuning error correction amount within the correction range is shown in the figure. Figure 7 shown.

[0052] S130. Determine the error correction amount of the well points in the entire area based on the well point stratification data and the seismic layer data.

[0053] Based on the pre-acquired well point layer data and seismic layer data, the error correction amount of the well point in the whole area is calculated. The calculation method of the error correction amount of the well point in the whole area is a public technology in this field, and this embodiment does not limit it. The schematic diagram of the error correction amount of the well point in the whole area within the correction range is shown in FIG. Figure 8 shown.

[0054] S140, determining a total correction amount according to the reference plane correction amount, the tuning error correction amount, and the area-wide well point error correction amount, and correcting the time-lapse seafloor node seismic data according to the total correction amount.

[0055] Furthermore, the reference surface correction amount, the tuning error correction amount, and the error correction amount of the well point in the entire area are added together to obtain the total correction amount.

[0056] In this embodiment, the reference plane correction is used to eliminate the error based on the reference plane, the tuning error correction is used to eliminate the function error when the formation thickness is less than the tuning thickness, and the full-area well point error correction is used to eliminate the error based on the well point. By adding these three corrections, the depth correction of the time-shifted seafloor node seismic data can be performed more accurately. The schematic diagram of the total correction amount within the correction range is shown in Figure 1. Figure 9 shown.

[0057] Furthermore, the depth correction is performed on the time-shifted seafloor node seismic data according to the total correction amount through a depth correction model to obtain correction data that matches the time-shifted seafloor node seismic data.

[0058] In this embodiment, a depth correction model is pre-established, illustratively, Figure 10 A schematic diagram of a depth correction model is provided. According to the total correction amount, depth correction is performed on the time-lapsed seafloor node seismic data within the correction range to obtain corrected data. Figure 11 : is a schematic diagram of a depth domain profile of a corrected time-lapse seafloor node seismic data provided by the first embodiment of the present invention, such as Figure 11 As shown in the figure, the seismic layer data and well point layer data of the corrected time-lapse seafloor node seismic data are consistent with each other and have small errors.

[0059] The technical solution of the embodiment of the present invention calculates the trace integral response of time-lapsed seafloor node seismic data, determines the correction range of the time-lapsed seafloor node seismic data, and calculates the reference level correction within the correction range. When the formation thickness is less than the tuning thickness, the tuning error correction is determined based on the trace integral response data and well data. The well point error correction for the entire area is determined based on the well point layer data and seismic horizon data. The total correction is determined based on the reference level correction, the tuning error correction, and the well point error correction for the entire area, and the time-lapsed seafloor node seismic data is corrected based on the total correction. This improves the depth accuracy of the time-lapsed seafloor node seismic data, thereby improving the success rate of remaining oil prediction and drilling.

[0060] Example 2

[0061] Figure 12 This is a schematic diagram of a correction device for time-shifted seafloor node seismic data provided by the second embodiment of the present invention. Figure 12 As shown, the device includes: a reference correction amount determination module 210, a tuning error correction amount determination module 220, a full-area well point error correction amount determination module 230, and a time-shifted seafloor node seismic data correction module 240. Among them:

[0062] A reference plane correction amount determination module 210 is configured to determine a correction range of the time-shifted seafloor node seismic data based on the trace integrated response data of the time-shifted seafloor node seismic data, and determine a reference plane correction amount within the correction range;

[0063] a tuning error correction amount determination module 220 for determining a tuning error correction amount based on the trace integrated response data and the well data if it is determined that the formation thickness is less than the tuning thickness;

[0064] The module 230 for determining the error correction amount of well points in the entire region is used to determine the error correction amount of well points in the entire region based on the well point layer data and the seismic layer data;

[0065] The time-lapse seafloor node seismic data correction module 240 is used to determine a total correction amount based on the reference plane correction amount, the tuning error correction amount and the full-area well point error correction amount, and correct the time-lapse seafloor node seismic data based on the total correction amount.

[0066] The technical solution of the embodiment of the present invention calculates the trace integral response of time-lapsed seafloor node seismic data, determines the correction range of the time-lapsed seafloor node seismic data, and calculates the reference level correction within the correction range. When the formation thickness is less than the tuning thickness, the tuning error correction is determined based on the trace integral response data and well data. The well point error correction for the entire area is determined based on the well point layer data and seismic horizon data. The total correction is determined based on the reference level correction, the tuning error correction, and the well point error correction for the entire area, and the time-lapsed seafloor node seismic data is corrected based on the total correction. This improves the depth accuracy of the time-lapsed seafloor node seismic data, thereby improving the success rate of remaining oil prediction and drilling.

[0067] Based on the above embodiment, the reference plane correction amount determination module 210 includes:

[0068] A maximum amplitude attribute value extraction unit is used to determine the trace integral response data of the time-shifted seafloor node seismic data and extract the maximum amplitude attribute value from the trace integral response data;

[0069] a maximum amplitude attribute value range determining unit, configured to determine a maximum amplitude attribute value range matching the flooding signal based on the spatial location of the production well and the spatial location of the injection well;

[0070] The correction range determination unit is used to determine the correction range of the time-shifted seafloor node seismic data according to the maximum amplitude attribute value range that matches the flooding signal.

[0071] Based on the above embodiment, the reference plane correction amount determination module 210 includes:

[0072] a zero phase level determination unit, configured to determine the zero phase level of the event axis of the channel integrated response data within the correction range;

[0073] The reference plane correction amount determining unit is used to take the difference between the zero-phase plane and the oil-water interface as the reference plane correction amount.

[0074] Based on the above embodiment, the tuning error correction amount determination module 220 includes:

[0075] a maximum value determining unit, configured to determine a maximum value of a maximum amplitude attribute value within the correction range;

[0076] A wavelength determination unit, for determining the wavelength of the statistical wavelet of the well bypass channel according to the well data;

[0077] The tuning error correction amount determining unit is configured to determine the tuning error correction amount according to the maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, the wavelength, and the maximum value of the maximum amplitude attribute value within the correction range.

[0078] Based on the above embodiment, the tuning error correction amount determination unit is specifically configured to:

[0079] The tuning error correction amount is determined according to the following formula:

[0080]

[0081] Wherein, T represents the tuning error correction amount, Max represents the maximum value of the maximum amplitude attribute value within the correction range, A represents the maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, and λ represents the wavelength of the statistical wavelet of the wellside trace.

[0082] Based on the above embodiment, the time-shifted seafloor node seismic data correction module 240 includes:

[0083] The total correction amount determination unit is used to add the reference surface correction amount, the tuning error correction amount and the whole area well point error correction amount to obtain the total correction amount.

[0084] Based on the above embodiment, the time-shift seabed node seismic data correction module 240 includes:

[0085] A depth correction unit is used to perform depth correction on the time-lapsed seabed node seismic data according to the total correction amount through a depth correction model to obtain correction data that matches the time-lapsed seabed node seismic data.

[0086] The correction device for time-shifted seafloor node seismic data provided by an embodiment of the present invention can execute the correction method for time-shifted seafloor node seismic data provided by any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0087] Example 3

[0088] Figure 13 A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0089] like Figure 13As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., which is communicatively connected to the at least one processor 11. The memory stores a computer program that can be executed by the at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. Various programs and data required for the operation of the electronic device 10 can also be stored in the RAM 13. The processor 11, ROM 12, and RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0090] Multiple components in the electronic device 10 are connected to the I / O interface 15, including an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0091] The processor 11 may be any general-purpose and / or specialized processing component with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the correction method for time-shifted seafloor node seismic data.

[0092] In some embodiments, the correction method for time-shifted seafloor node seismic data can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the correction method for time-shifted seafloor node seismic data described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to perform the correction method for time-shifted seafloor node seismic data by any other appropriate means (e.g., by means of firmware).

[0093] Various embodiments of the systems and techniques described herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system that includes at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0094] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program may be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0095] In the context of the present invention, computer-readable storage media can be tangible media that can contain or store a computer program for use with an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. Computer-readable storage media can include but are not limited to electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, computer-readable storage media can be machine-readable signal media. More specific examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0096] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0097] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0098] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0099] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in the present invention can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved. This is not limited herein.

[0100] The above specific embodiments do not limit the scope of protection of the present invention. Those skilled in the art will appreciate that various modifications, combinations, sub-combinations, and substitutions may be made based on design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention are intended to be included within the scope of protection of the present invention.

Claims

1. A correction method for time-shifted seafloor node seismic data, characterized in that: include: determining a correction range of the time-lapsed seafloor node seismic data based on the trace integrated response data of the time-lapsed seafloor node seismic data, and determining a reference plane correction amount within the correction range; If it is determined that the formation thickness is less than the tuning thickness, determining a tuning error correction amount based on the trace integrated response data and the well data; Determine the error correction amount of well points in the entire area based on well point layer data and seismic layer data; A total correction amount is determined based on the reference plane correction amount, the tuning error correction amount, and the whole-area well point error correction amount, and the time-lapse seabed node seismic data is corrected based on the total correction amount.

2. The method according to claim 1, characterized in that According to the trace integral response data of the time-lapsed seafloor node seismic data, the correction range of the time-lapsed seafloor node seismic data is determined, including: Determine the trace integral response data of the time-lapse seafloor node seismic data, and extract the maximum amplitude attribute value from the trace integral response data; According to the spatial location of the production well and the spatial location of the water injection well, the maximum amplitude attribute value range matching the flooding signal is determined; The correction range of the time-shifted seafloor node seismic data is determined based on the maximum amplitude attribute value range that matches the flooding signal.

3. The method according to claim 2, characterized in that Determining a reference surface correction amount within the correction range includes: Determining the zero phase plane of the event axis of the channel integrated response data within the correction range; The difference between the zero-phase plane and the oil-water interface is used as the reference plane correction amount.

4. The method according to claim 2, characterized in that Determining a tuning error correction amount based on the trace integrated response data and the well data includes: determining a maximum value of the maximum amplitude attribute value within the correction range; According to the well data, the wavelength of the statistical wavelet of the well bypass is determined; The tuning error correction amount is determined according to the maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, the wavelength, and the maximum value of the maximum amplitude attribute value within the correction range.

5. The method according to claim 4, characterized in that Determining a tuning error correction amount according to a maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, the wavelength, and a maximum value of the maximum amplitude attribute value within the correction range, including: The tuning error correction amount is determined according to the following formula: Wherein, T represents the tuning error correction amount, Max represents the maximum value of the maximum amplitude attribute value within the correction range, A represents the maximum amplitude attribute value of the trace integrated response data of each well point within the correction range, and λ represents the wavelength of the statistical wavelet of the wellside trace.

6. The method according to claim 1, characterized in that Based on the datum correction, tuning error correction, and all-area well point error correction, the total correction is determined, including: The total correction is obtained by adding the reference correction, tuning error correction and all-area well point error correction.

7. The method according to claim 1, characterized in that Correcting the time-shifted seafloor node seismic data according to the total correction amount includes: The depth correction model is used to perform depth correction on the time-shifted seafloor node seismic data according to the total correction amount, thereby obtaining correction data that matches the time-shifted seafloor node seismic data.

8. A correction device for time-shifted seafloor node seismic data, characterized in that: include: a reference plane correction amount determination module, configured to determine a correction range of the time-lapsed seafloor node seismic data based on the trace integral response data of the time-lapsed seafloor node seismic data, and determine a reference plane correction amount within the correction range; a tuning error correction amount determination module, configured to determine a tuning error correction amount based on the trace integrated response data and the well data if it is determined that the formation thickness is less than the tuning thickness; The module for determining the error correction amount of well points in the whole area is used to determine the error correction amount of well points in the whole area based on the well point layer data and seismic layer data; The time-lapse seabed node seismic data correction module is used to determine the total correction amount based on the reference plane correction amount, the tuning error correction amount and the full-area well point error correction amount, and correct the time-lapse seabed node seismic data based on the total correction amount.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the method for correcting time-shifted seafloor node seismic data as described in any one of claims 1 to 7 is implemented.

10. A storage medium storing computer executable instructions, characterized in that: When the computer executable instructions are executed by a computer processor, they are used to perform the method for correcting time-shifted seafloor node seismic data according to any one of claims 1 to 7.

Images