Seismic data processing method, device, computer equipment, storage medium and product

Abstract

The present application provides a seismic data processing method, device, computer equipment, storage medium, and product, belonging to the field of seismic exploration technology. The method includes: for each seismic trace, based on the time history information of the seismic trace, the sampling time point corresponding to the time history information and the sampling time point after it are used as the first sampling point of the seismic trace, and the sampling time point before the sampling time point corresponding to the time history information is used as the second sampling point of the seismic trace; based on the signal value of the first sampling point of the seismic trace, the signal value of the second sampling point of the seismic trace is transformed to obtain the first signal value of the second sampling point of the seismic trace; the first signal values ​​of the second sampling points of multiple seismic traces are respectively subjected to gain processing to obtain the second signal values ​​of the second sampling points of multiple seismic traces; for each seismic trace, based on the second signal value of the second sampling point of the seismic trace, the signal value of the first sampling point of the seismic trace is inversely transformed. This method improves the display effect of the vertical seismic profile.

Description

Technical Field

[0001] The present application relates to the field of seismic exploration technology, and in particular to a seismic data processing method, device, computer equipment, storage medium and product. Background Art

[0002] With the advancement of seismic exploration technology, a growing number of seismic observation methods have emerged, such as surface seismic profiling and vertical seismic profiling (VSP). The surface seismic profiling method collects seismic waves on the surface and generates horizontal seismic profiles based on adjacent seismic traces in the horizontal direction. The vertical seismic profiling method collects seismic waves in a well and generates vertical seismic profiles based on adjacent seismic traces in the depth direction. In the horizontal direction, adjacent traces have similar surface conditions and consistent excitation conditions, so there is no significant energy difference between adjacent traces in a horizontal seismic profile. However, in a vertical seismic profile, the distance between adjacent traces in the depth direction is large. When a thin energy shielding layer (such as igneous rock) exists in the subsurface, the energy difference between adjacent traces becomes significant, reducing the visual quality of the vertical seismic profile. Therefore, a seismic data processing method is needed to process the seismic profile data of a vertical seismic profile to reduce the energy difference between traces and thereby improve the visual quality of the vertical seismic profile. Summary of the Invention

[0003] The present invention provides a method, apparatus, computer device, storage medium, and product for processing seismic data, which can improve the display effect of vertical seismic profiles. The technical solution is as follows:

[0004] In one aspect, a method for processing seismic data is provided, the method comprising:

[0005] Determining, from seismic profile data of an observation well, time history information of a plurality of seismic traces corresponding to a plurality of depth points in the observation well, the seismic profile data including signal values ​​of the plurality of seismic traces at a plurality of sampling time points, the time history information of the plurality of seismic traces being one-way time sampling time points or two-way time sampling time points when seismic waves arrive at the seismic traces;

[0006] For each seismic trace, based on the time history information of the seismic trace, a sampling time point corresponding to the time history information and a sampling time point after the sampling time point corresponding to the time history information are used as the first sampling point of the seismic trace, and a sampling time point before the sampling time point corresponding to the time history information is used as the second sampling point of the seismic trace;

[0007] For each seismic trace, transforming a signal value at a second sampling point of the seismic trace based on a signal value at a first sampling point of the seismic trace to obtain a first signal value at the second sampling point of the seismic trace;

[0008] performing gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces respectively to obtain second signal values ​​of the second sampling points of the plurality of seismic traces;

[0009] For each seismic trace, the signal value of the first sampling point of the seismic trace is inversely transformed based on the second signal value of the second sampling point of the seismic trace, and the signal values ​​of the first sampling points of multiple seismic traces after inverse transformation are used to determine the vertical seismic profile.

[0010] In some embodiments, there are multiple first sampling points and multiple second sampling points, and for each seismic trace, transforming the signal value of the second sampling point of the seismic trace based on the signal value of the first sampling point of the seismic trace to obtain the first signal value of the second sampling point of the seismic trace includes:

[0011] For each seismic trace, determining a plurality of target sampling points among a plurality of second sampling points based on the number of the plurality of first sampling points, wherein the number of the plurality of target sampling points is the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points is the starting sampling point among the plurality of second sampling points;

[0012] The signal values ​​of the plurality of target sampling points are transformed into the signal values ​​of the plurality of first sampling points to obtain the first signal values ​​of the plurality of target sampling points.

[0013] In some embodiments, for each seismic trace, performing an inverse transformation on a signal value of a first sampling point of the seismic trace based on a second signal value of a second sampling point of the seismic trace comprises:

[0014] The signal values ​​of the plurality of first sampling points are transformed into second signal values ​​of the plurality of target sampling points.

[0015] In some embodiments, performing gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces to obtain the second signal values ​​of the second sampling points of the plurality of seismic traces includes:

[0016] For each sampling time point, the signal values ​​of the multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point;

[0017] Dividing the signal set based on a reference depth window to obtain a plurality of first subsets;

[0018] For each first subset, performing gain processing on the first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset;

[0019] Second signal values ​​of second sampling points of the plurality of seismic traces are determined from a plurality of second subsets respectively corresponding to the plurality of sampling time points.

[0020] In some embodiments, for each first subset, performing gain processing on the first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset includes:

[0021] For each first subset, determining a gain parameter based on a plurality of signal values ​​in the first subset and a number of the plurality of signal values;

[0022] Products of a plurality of signal values ​​in the first subset and the gain parameter are respectively determined to obtain a second subset corresponding to the first subset.

[0023] In some embodiments, performing gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces to obtain the second signal values ​​of the second sampling points of the plurality of seismic traces includes:

[0024] For each sampling time point, the signal values ​​of the multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point;

[0025] For each signal set, gain processing is performed on the signal set based on multiple signal values ​​included in the signal set, and second signal values ​​of the multiple second sampling points are determined from the multiple gain-processed signal sets corresponding to the multiple sampling time points respectively.

[0026] In another aspect, a seismic data processing device is provided, comprising:

[0027] an information determination module, configured to determine, from seismic profile data of an observation well, time history information of a plurality of seismic traces corresponding to a plurality of depth points in the observation well, the seismic profile data comprising signal values ​​of the plurality of seismic traces at a plurality of sampling time points, the time history information of the plurality of seismic traces being one-way time sampling time points or two-way time sampling time points when seismic waves arrive at the seismic traces;

[0028] a sampling point determination module configured to, for each seismic trace, based on time history information of the seismic trace, determine a sampling time point corresponding to the time history information and a sampling time point after the sampling time point corresponding to the time history information as a first sampling point of the seismic trace, and determine a sampling time point before the sampling time point corresponding to the time history information as a second sampling point of the seismic trace;

[0029] a transformation module, configured to transform, for each seismic trace, a signal value of a second sampling point of the seismic trace based on a signal value of a first sampling point of the seismic trace, to obtain a first signal value of the second sampling point of the seismic trace;

[0030] a gain module, configured to perform gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces respectively to obtain second signal values ​​of the second sampling points of the plurality of seismic traces;

[0031] An inverse transformation module is used to perform an inverse transformation on the signal value of the first sampling point of each seismic trace based on the second signal value of the second sampling point of the seismic trace, and the signal values ​​of the first sampling points of multiple seismic traces after inverse transformation are used to determine a vertical seismic profile.

[0032] In some embodiments, there are multiple first sampling points and multiple second sampling points, and the transformation module is used to:

[0033] For each seismic trace, determining a plurality of target sampling points among a plurality of second sampling points based on the number of the plurality of first sampling points, wherein the number of the plurality of target sampling points is the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points is the starting sampling point among the plurality of second sampling points;

[0034] The signal values ​​of the plurality of target sampling points are transformed into the signal values ​​of the plurality of first sampling points to obtain the first signal values ​​of the plurality of target sampling points.

[0035] In some embodiments, the inverse transform module is configured to:

[0036] The signal values ​​of the plurality of first sampling points are transformed into second signal values ​​of the plurality of target sampling points.

[0037] In some embodiments, the gain module is configured to:

[0038] For each sampling time point, the signal values ​​of the multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point;

[0039] Dividing the signal set based on a reference depth window to obtain a plurality of first subsets;

[0040] For each first subset, performing gain processing on the first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset;

[0041] Second signal values ​​of second sampling points of the plurality of seismic traces are determined from a plurality of second subsets respectively corresponding to the plurality of sampling time points.

[0042] In some embodiments, the gain module is configured to:

[0043] For each first subset, determining a gain parameter based on a plurality of signal values ​​in the first subset and a number of the plurality of signal values;

[0044] Products of a plurality of signal values ​​in the first subset and the gain parameter are respectively determined to obtain a second subset corresponding to the first subset.

[0045] In some embodiments, the gain module is configured to:

[0046] For each sampling time point, the signal values ​​of the multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point;

[0047] For each signal set, gain processing is performed on the signal set based on multiple signal values ​​included in the signal set, and second signal values ​​of the multiple second sampling points are determined from the multiple gain-processed signal sets corresponding to the multiple sampling time points respectively.

[0048] On the other hand, a computer device is provided, which includes one or more processors and one or more memories, wherein at least one program code is stored in the one or more memories, and the at least one program code is loaded and executed by the one or more processors to implement the seismic data processing method described in any of the above implementation methods.

[0049] On the other hand, a computer-readable storage medium is provided, in which at least one program code is stored. The at least one program code is loaded and executed by a processor to implement the seismic data processing method described in any of the above implementations.

[0050] On the other hand, a computer program product is provided, which includes a computer program code, the computer program code is stored in a computer-readable storage medium, a processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code, so that the computer device executes the seismic data processing method described in any of the above-mentioned implementation methods.

[0051] In the embodiment of the present application, since signals are only acquired at the sampling time point corresponding to the time-series information and the sampling time points thereafter, the first sampling point is the valid sampling time point for acquiring signals, and the signal value at the second sampling point is zero. Furthermore, since the second sampling point occurs before the first sampling point, the signal value of the second sampling point is transformed based on the signal value of the first sampling point, thereby shifting the signal value of the first sampling point forward. Furthermore, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, thereby aligning the signals of the multiple seismic traces at the sampling time points. Gain processing is then performed on this basis. Since the signals of the multiple seismic traces are aligned based on the sampling time points, the regularity and convenience of the gain processing can be improved, thereby enhancing the gain processing effect. The signal value of the first sampling point is then inversely transformed, and the gain-processed signal is restored to the first sampling point, thus achieving effective gain processing of the seismic profile data. This solves the energy step problem in the generated vertical seismic profile, thereby improving the display effect of the vertical seismic profile. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate the technical solutions in the embodiments of the present application, 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 application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0053] Figure 1 This is a schematic diagram of an implementation environment provided by an embodiment of the present application;

[0054] Figure 2 is a flow chart of a seismic data processing method provided in an embodiment of the present application;

[0055] Figure 3 is a flow chart of another seismic data processing method provided in an embodiment of the present application;

[0056] Figure 4 is a schematic diagram of a vertical seismic profile provided in an embodiment of the present application;

[0057] Figure 5 This is a schematic diagram of the effect of data changes provided by an embodiment of the present application;

[0058] Figure 6 is a schematic diagram of another vertical seismic profile provided in an embodiment of the present application;

[0059] Figure 7 is a schematic diagram of another vertical seismic profile provided in an embodiment of the present application;

[0060] Figure 8 is a schematic diagram of another vertical seismic profile provided in an embodiment of the present application;

[0061] Figure 9 is a schematic diagram of another vertical seismic profile provided in an embodiment of the present application;

[0062] Figure 10 is a block diagram of a seismic data processing device provided in an embodiment of the present application;

[0063] Figure 11 This is a block diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0064] In order to make the objectives, technical solutions and advantages of this application clearer, the implementation methods of this application will be further described in detail below with reference to the accompanying drawings.

[0065] The terms "first," "second," "third," and "fourth," etc. in the specification and claims of this application and the accompanying drawings are used to distinguish different objects, not to describe a specific order. 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 apparatus comprising a series of steps or elements is not limited to the listed steps or elements, but may optionally include steps or elements not listed, or may optionally include other steps or elements inherent to the process, method, product, or apparatus.

[0066] The above are only optional embodiments of the present application and are not intended to limit the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application shall be included in the scope of protection of the present application. It should be noted that the information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data used for analysis, stored data, displayed data, etc.) and signals involved in this application are all authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with the relevant laws, regulations and standards of the relevant countries and regions. For example, the seismic profile data involved in this application were all obtained with full authorization.

[0067] Figure 1This is a schematic diagram of the implementation environment of a seismic data processing method provided by this application. Figure 1 The implementation environment includes: a computer device; the computer device can be provided as a terminal 101, or can be provided as a terminal 101 and a server 102, which is not specifically limited.

[0068] If the computer device is provided as terminal 101 and a target application is installed on terminal 101, the user can log in to the target application and then process the seismic profile data using the method provided in this application to generate a vertical seismic profile diagram.

[0069] If the computer device is provided as a terminal 101 and a server 102, the terminal 101 and the server 102 can be connected via a wireless or wired network. Accordingly, the terminal 101 has a target application installed, and the server 102 is the server 102 corresponding to the target application. Users process seismic profile data using the target application, and the server 102 provides background services. In the embodiments of this application, the computer device is provided as a terminal 101 as an example for description.

[0070] The terminal 101 is at least one of a mobile phone, a tablet computer, a PC (Personal Computer), an intelligent voice interaction device, and an in-vehicle terminal. The server 102 is at least one of a single server, a server cluster consisting of multiple servers, a cloud server, a cloud computing platform, and a virtualization center.

[0071] The present invention provides a method for processing seismic data. Figure 2 , methods include:

[0072] 201. The computer equipment determines the time history information of multiple seismic channels corresponding to multiple depth points in the observation well from the seismic profile data of the observation well. The seismic profile data includes signal values ​​of the multiple seismic channels at multiple sampling time points. The time history information of the multiple seismic channels is the one-way sampling time point or the two-way sampling time point when the seismic wave reaches the seismic channel.

[0073] In the embodiments of the present application, the signal value is an amplitude value. A computer device performs first-arrival wave picking and dynamic motion correction processing on seismic data to obtain seismic profile data. The one-way time sampling point for each seismic trace refers to the sampling time point when the seismic wave first arrives at the seismic trace before dynamic motion correction. The two-way time sampling point for each seismic trace refers to the sampling time point when the seismic wave first arrives at the seismic trace after dynamic motion correction.

[0074] In the embodiment of the present application, if the time history information is one-way time information, that is, the one-way time sampling time point when the seismic wave arrives at the seismic trace, then the seismic profile data is the seismic profile data before dynamic dynamic correction, and accordingly, the computer device processes the seismic profile data before dynamic dynamic correction. If the time history information is two-way time information, that is, the two-way time sampling time point when the seismic wave arrives at the seismic trace, then the seismic profile data is the seismic profile data after dynamic dynamic correction, and accordingly, the computer device processes the seismic profile data after dynamic dynamic correction.

[0075] In the embodiment of the present application, the one-way sampling time point and the two-way sampling time point are respectively sampling time points among a plurality of sampling time points. The plurality of sampling time points can be distributed in sequence according to a fixed time interval, or the time interval can be set as needed. The plurality of seismic traces can be distributed in sequence according to a fixed depth interval, or the depth interval can be set as needed, which is not specifically limited here. Accordingly, when determining the time history information from the seismic profile data, if the time history information is one-way time information, the sampling time point at which the signal is first collected for each seismic trace in the seismic profile data before dynamic correction is used as the time history information of each seismic trace. If the time history information is two-way time information, the sampling time point at which the signal is first collected for each seismic trace in the seismic profile data after dynamic correction is used as the time history information of each seismic trace.

[0076] 202. For each seismic trace, the computer device uses, based on the time information of the seismic trace, a sampling time point corresponding to the time information and a sampling time point after the sampling time point corresponding to the time information as the first sampling point of the seismic trace, and uses a sampling time point before the sampling time point corresponding to the time information as the second sampling point of the seismic trace.

[0077] In this embodiment of the present application, for each seismic trace, no valid signal is collected at sampling time points before the time point at which the seismic wave is received, and thus the signal values ​​corresponding to these sampling time points are zero. Accordingly, the first sampling point is the sampling time point at which the seismic trace receives the seismic wave, and the second sampling point is the sampling time point at which the seismic trace does not receive the seismic wave.

[0078] 203. For each seismic trace, the computer device transforms the signal value of the second sampling point of the seismic trace based on the signal value of the first sampling point of the seismic trace to obtain the first signal value of the second sampling point of the seismic trace.

[0079] In the embodiment of the present application, the signal value of the first sampling point is assigned to the second sampling point. Since the time of the second sampling point is before the first sampling point, the sampling time point of the signal value of the first sampling point is shifted forward.

[0080] 204. The computer device performs gain processing on the first signal values ​​of the second sampling points of the multiple seismic channels to obtain the second signal values ​​of the second sampling points of the multiple seismic channels.

[0081] In the embodiment of the present application, the gain processing is used to equalize the first signal values ​​of the multiple second sampling points, that is, to reduce the difference between the first signal values ​​of the multiple second sampling points. Optionally, performing gain processing on the multiple signal values ​​means uniformly reducing the multiple signal values.

[0082] 205. The computer device performs an inverse transformation on the signal value of the first sampling point of each seismic trace based on the second signal value of the second sampling point of the seismic trace, and the signal values ​​of the first sampling points after the inverse transformation of the multiple seismic traces are used to determine the vertical seismic profile.

[0083] In the embodiment of the present application, the second signal value of the second sampling point obtained by gain processing is assigned to the signal value of the first sampling point, thereby shifting the sampling time point of the second signal value of the second sampling point backward.

[0084] In the embodiment of the present application, the signal value of the first sampling point is shifted to the second sampling point for gain, and then the gained signal value is restored to the first sampling point, thereby achieving gain on the signal value of the first sampling point.

[0085] In the embodiment of the present application, since signals are only acquired at the sampling time point corresponding to the time-series information and the sampling time points thereafter, the first sampling point is the valid sampling time point for acquiring signals, and the signal value at the second sampling point is zero. Furthermore, since the second sampling point occurs before the first sampling point, the signal value of the second sampling point is transformed based on the signal value of the first sampling point, thereby shifting the signal value of the first sampling point forward. Furthermore, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, thereby aligning the signals of the multiple seismic traces at the sampling time points. Gain processing is then performed on this basis. Since the signals of the multiple seismic traces are aligned based on the sampling time points, the regularity and convenience of the gain processing can be improved, thereby enhancing the gain processing effect. The signal value of the first sampling point is then inversely transformed, and the gain-processed signal is restored to the first sampling point, thus achieving effective gain processing of the seismic profile data. This solves the energy step problem in the generated vertical seismic profile, thereby improving the display effect of the vertical seismic profile.

[0086] above Figure 2 The embodiment is the basic process of seismic data processing. Figure 3 The embodiment of the present invention further illustrates the process of seismic data processing. Figure 3 , Figure 3 This is a flowchart of a seismic profile data processing method provided in an embodiment of the present application, which includes the following steps.

[0087] 301. The computer device acquires the seismic data collected in the observation well and obtains seismic profile data based on the seismic data.

[0088] In this embodiment, when collecting seismic data, a location on the surface at a reference distance from the wellhead of an observation well is determined as the source excitation point. A surface source is used to generate seismic waves, which are received by a three-component geophone in the observation well to obtain seismic data. The observation well can be a vertical well, a newly drilled well, an old well, or a renovated well. The reference distance between the source excitation point and the observation well meets the zero-well source excitation distance requirement specified in the industry standard SY / T 7450, ensuring both wellsite operational safety and drilling safety.

[0089] In an embodiment of the present application, the process of obtaining seismic profile data based on seismic data by a computer device includes the following steps: the computer device performs first arrival wave picking processing on the seismic data to obtain first profile data, where the first profile data includes one-way time sampling time points when the seismic waves respectively arrive at multiple seismic traces. The computer device performs dynamic correction processing on the first profile data to obtain second profile data, where the second profile data includes two-way time sampling time points when the seismic waves respectively arrive at multiple seismic traces. The first arrival wave picking processing and the dynamic correction processing process meet the industry standard SY / T 7450-2019.

[0090] 302. The computer device determines the time history information of multiple seismic channels corresponding to multiple depth points in the observation well from the seismic profile data of the observation well. The seismic profile data includes signal values ​​of the multiple seismic channels at multiple sampling time points. The time history information of the multiple seismic channels is the one-way sampling time point or the two-way sampling time point when the seismic wave reaches the seismic channel.

[0091] In the embodiment of the present application, the time history information is the one-way time information before dynamic correction or the two-way time information after dynamic correction. The one-way time information and the two-way time information refer to the one-way time sampling time point and the two-way time sampling time point when the seismic wave arrives at the seismic trace, respectively.

[0092] See also Figure 4 , Figure 4 : is a schematic diagram of a vertical seismic profile provided in an embodiment of the present application. Figure 4 The intersections of the thick oblique lines in (a) and multiple seismic traces are the one-way time information of the multiple seismic traces. Figure 4 The intersections of the thick oblique lines and multiple seismic traces in (b) are the two-way time information of the multiple seismic traces.

[0093] 303. For each seismic trace, the computer device uses, based on the time information of the seismic trace, a sampling time point corresponding to the time information and a sampling time point after the sampling time point corresponding to the time information as the first sampling point of the seismic trace, and uses a sampling time point before the sampling time point corresponding to the time information as the second sampling point of the seismic trace.

[0094] For example, see Figure 4 , Figure 4 (a) and Figure 4 The area where the wavy line in the lower left corner of (b) is located is the distribution area of ​​the first sampling point. Figure 4 (a) and Figure 4 The area where the vertical line in the upper right corner of (b) is located is the distribution area of ​​the second sampling point.

[0095] 304. For each seismic trace, the computer device transforms the signal value of the second sampling point of the seismic trace based on the signal value of the first sampling point of the seismic trace to obtain the first signal value of the second sampling point of the seismic trace.

[0096] In an embodiment of the present application, for each seismic trace, there are multiple first sampling points and multiple second sampling points. The number of first sampling points may be less than, greater than, or equal to the number of second sampling points. Accordingly, the computer device transforms the signal value of the second sampling point of each seismic trace based on the signal value of the first sampling point of the seismic trace to obtain the first signal value of the second sampling point of the seismic trace, including the following implementation methods.

[0097] A first implementation method: For each seismic trace, if the number of first sampling points is less than the number of second sampling points, the computer device determines a plurality of target sampling points from the plurality of second sampling points based on the number of the plurality of first sampling points. The number of the plurality of target sampling points is the same as the number of the plurality of first sampling points, and a starting sampling point of the plurality of target sampling points is the starting sampling point of the plurality of second sampling points. The computer device converts signal values ​​of the plurality of target sampling points into signal values ​​of the plurality of first sampling points to obtain first signal values ​​of the plurality of target sampling points.

[0098] The starting sampling point among the second sampling points is also the starting sampling time point among the multiple sampling time points, and the starting sampling time point is the earliest in time. Optionally, the starting sampling time point is zero. In this implementation, starting from the smallest sampling time point among the multiple second sampling points, a corresponding number of second sampling points are intercepted backward to obtain a plurality of target sampling points equal in number to the first sampling points. The multiple first sampling points are arranged sequentially in chronological order, and the multiple target sampling points are also arranged sequentially in chronological order, and then the signal value of each target sampling point is converted to the signal value of the first sampling point corresponding thereto. For example, the signal value of the first target sampling point is converted to the signal value of the first first sampling point, the signal value of the second target sampling point is converted to the signal value of the second first sampling point, and so on. Optionally, the signal values ​​of the second sampling points other than the target sampling points among the multiple second sampling points remain unchanged.

[0099] In this implementation, the signal values ​​of multiple target sampling points that are ranked first in the second sampling points and have the same number as the first sampling points are converted into signal values ​​of multiple first sampling points, so that the signals of the first sampling points are flattened to the starting sampling point, thereby achieving alignment of the signals of multiple seismic channels at the starting sampling point, thereby facilitating subsequent simultaneous gain processing of the signals of multiple seismic channels.

[0100] A second implementation method: For each seismic trace, if the number of first sampling points is greater than the number of second sampling points, the computer device determines a plurality of third sampling points based on the number of the plurality of first sampling points, where the plurality of third sampling points includes the plurality of second sampling points and a first number of first sampling points, the sum of the first number and the number of the plurality of second sampling points being equal to the number of the plurality of first sampling points, and the starting sampling point of the first number of first sampling points being the starting sampling point of the plurality of first sampling points. Signal values ​​of the plurality of third sampling points are converted into signal values ​​of the plurality of first sampling points to obtain first signal values ​​of the plurality of third sampling points.

[0101] In this implementation, starting from the smallest sampling time point among the multiple first sampling points, a corresponding number of first sampling points are truncated backward to obtain a plurality of third sampling points equal in number to the first sampling points. The multiple first sampling points are arranged sequentially in chronological order, and the multiple third sampling points are also arranged sequentially in chronological order. The signal value of each third sampling point is then transformed into the signal value of the first sampling point corresponding thereto. For example, the signal value of the first third sampling point is transformed into the signal value of the first first sampling point, the signal value of the second third sampling point is transformed into the signal value of the second first sampling point, and so on. In this implementation, when the number of first sampling points is greater than the number of second sampling points, the signal values ​​of some first sampling points are also transformed based on the signal values ​​of the first sampling points. This ensures the integrity of the signal values ​​of the first sampling points during the forward shift, thereby improving the precision and accuracy of the data and facilitating subsequent gain processing.

[0102] A third implementation method: For each seismic trace, if the number of first sampling points is equal to the number of second sampling points, the computer device converts the signal values ​​of the multiple second sampling points into the signal values ​​of the multiple first sampling points, thereby obtaining the first signal values ​​of the multiple second sampling points. The multiple first sampling points are arranged sequentially in chronological order, and the multiple second sampling points are also arranged sequentially in chronological order. The signal value of each second sampling point is then converted into the signal value of the first sampling point corresponding thereto. For example, the signal value of the first second sampling point is converted into the signal value of the first first sampling point, the signal value of the second second sampling point is converted into the signal value of the second first sampling point, and so on. In this implementation method, the conversion is performed directly, which improves conversion efficiency.

[0103] It should be noted that since the signal value at the first sampling point is subsequently inversely transformed based on the second signal value at the second sampling point, after the signal value at the second sampling point is transformed, the signal value at the first sampling point may remain unchanged or become zero. In the embodiments of the present application, the example in which the signal value at the first sampling point becomes zero is used for illustration, thereby facilitating subsequent gain processing.

[0104] For example, the number of the multiple sampling time points is M, the number of the multiple seismic traces is N, and M and N are integers greater than 0, then the seismic profile data before transformation can be expressed as A(M, N). Accordingly, the computer device can perform the transformation according to the following formula (1).

[0105] B(1:p j , j) = A(k j :k j +p j -1,j),j=1,2……N (1);

[0106] Where A represents the signal value of the first sampling point; B represents the first signal value of the second sampling point; k j represents the kth sampling time point corresponding to the time history information of the jth seismic trace; p j The number of sampling time points p from the kth sampling time point to the last sampling time point corresponding to the time history information of the jth seismic trace is k j +p j -1=M;1:p j represents the time from the first sampling point to the pth sampling point of the jth seismic trace; k j :k j +p j It represents the time from the kth sampling point to the k+p-1th sampling point of the jth seismic trace; that is, B(1:p j, j) represents the signal value from the first sampling time point to the pth sampling time point of the jth seismic trace; A(k j :k j +p j -1,j) represents the signal value from the kth sampling time point to the k+p-1th sampling time point of the jth seismic trace.

[0107] For example, see Figure 5 , Figure 5 This is a schematic diagram of the effect of data transformation provided by an embodiment of the present application. Figure 5 The upgoing wave field in (a) is transformed into Figure 5 In the upgoing wavefield in (b), the signals of multiple seismic traces are flattened at the initial sampling time point.

[0108] For example, see Figure 6 , Figure 6 Schematic diagram of a vertical seismic profile provided in an embodiment of the present application, wherein the vertical seismic profile is obtained based on signal values ​​transformed based on double-path time information, and the signals of multiple seismic traces are flattened to the starting sampling time point.

[0109] 305. The computer device performs gain processing on the first signal values ​​of the second sampling points of the multiple seismic channels to obtain the second signal values ​​of the second sampling points of the multiple seismic channels.

[0110] In some embodiments, the above-mentioned computer device performs gain processing on the first signal values ​​of the second sampling points of multiple seismic channels respectively to obtain the second signal values ​​of the second sampling points of the multiple seismic channels, which includes the following steps: for each sampling time point, the computer device sorts the signal values ​​of the multiple seismic channels at the sampling time point in sequence to form a signal set; wherein, for each seismic channel, if the sampling time point is the second sampling point of the seismic channel, the signal value of the seismic channel at the sampling time point is the first signal value of the second sampling point; based on the reference depth window, the signal set is divided to obtain multiple first subsets; for each first subset, based on the multiple signal values ​​included in the first subset, the first subset is gain processed to obtain the second subset corresponding to the first subset; and the second signal values ​​of the second sampling points of the multiple seismic channels are determined from the multiple second subsets corresponding to the multiple sampling time points.

[0111] It should be noted that for each seismic trace, if the number of first sampling points is less than the number of second sampling points, then for each sampling time point, if the sampling time point is the target sampling point of the seismic trace, the signal value of the seismic trace at that sampling time point is the first signal value of the target sampling point. For each seismic trace, if the number of first sampling points is greater than the number of second sampling points, then for each sampling time point, if the sampling time point is the third sampling point of the seismic trace, the signal value of the seismic trace at that sampling time point is the first signal value of the third sampling point. For each seismic trace, if the number of first sampling points is equal to the number of second sampling points, then for each sampling time point, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at that sampling time point is the first signal value of the second sampling point.

[0112] In the embodiment of the present application, a signal set X can be obtained by extracting a set of signal values ​​along the depth direction according to the following formula (2): i (N), then corresponding to the depth direction of M sampling time points, M groups of signal sets can be obtained.

[0113] X i (1:N)=B(i,1:N),i=1,2...M (2);

[0114] Among them, X i represents the signal set at the i-th sampling time point; 1:N represents the signal set from the 1st seismic trace to the Nth seismic trace. That is, the signal set at the i-th sampling time point consists of the signal values ​​of the 1st to Nth seismic traces at that sampling time point.

[0115] Among them, the computer equipment will signal set X i (1:N) is divided according to the reference depth window to obtain W first subsets X iw (d). Wherein, W represents the number of the first subset, and d represents the number of signal values ​​included in the first subset. W and d can be obtained by the following formula (3).

[0116]

[0117] Where Depth represents the depth of the measured observation well; ΔDep represents the depth interval, i.e., the reference depth window; Int represents a floor function, for example, Int(5 / 3)=1; and Cdep represents the signal acquisition depth interval. The reference depth window is a sliding depth window. The computer device can determine the reference depth window based on the law of formation change. For example, the reference depth window can be determined based on the formation lithology. If the change in formation lithology is not greater than the reference change degree, the reference depth window is greater than the preset depth window; if the change in formation lithology is greater than the reference change degree, the reference depth window is not greater than the preset depth window. For example, if the depth of the measured observation well is 8,000 meters and the reference depth window is 400 meters, 20 first subsets are obtained. The signal acquisition depth interval within each first subset is 20 meters, so a first subset includes 20 signal values.

[0118] In some embodiments, the above-mentioned computer device performs gain processing on each first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset, including the following steps: the computer device determines a gain parameter for each first subset based on multiple signal values ​​and the number of multiple signal values ​​in the first subset; and determines the product of multiple signal values ​​in the first subset and the gain parameter respectively to obtain the second subset corresponding to the first subset.

[0119] In some embodiments, the process of determining a gain parameter by a computer device based on multiple signal values ​​and the number of multiple signal values ​​in a first subset includes the following steps: the computer device determines the cumulative value of the squares of the multiple signal values, and uses the quotient of the number of multiple signal values ​​and the cumulative value as the gain parameter. In this embodiment, gain processing is performed on the signal values ​​based on the signal values ​​in the subset and the number of signal values, which can improve the pertinence and accuracy of the gain processing, thereby improving the effect of the gain processing. In the embodiments of the present application, the gain parameter can also be determined by other means, which are not specifically limited here. For example, the computer device uses the quotient of the number of multiple signal values ​​and the sum of the multiple signal values ​​as the gain parameter.

[0120] For example, referring to the following formula (4), the computer device performs gain processing on the first subset through the following formula (4).

[0121]

[0122] Among them, X iw represents the first subset of the i-th element, Y iw represents the second subset of i, X iw (l) represents the lth signal value in the first subset.

[0123] The multiple second subsets of the multiple sampling time points include the second signal values ​​of the second sampling points of the multiple seismic traces. Accordingly, for each sampling time point, the computer device sequentially arranges the multiple second subsets corresponding to the sampling time point in depth order to obtain the multiple signal values ​​after gain at the sampling time point. Step 305 is repeatedly executed for each sampling time point to obtain the gain data Y(M, W×d) after the automatic gain of the M groups of signal values. Then, for each seismic trace, based on the multiple second sampling points of the seismic trace, the second signal values ​​of the multiple second sampling points of the seismic trace are found from the gain data.

[0124] It should be noted that if the number of first sampling points is less than the number of second sampling points, then for each seismic trace, based on the multiple target sampling points of the seismic trace, the second signal values ​​of the multiple target sampling points of the seismic trace are found from the gain data. If the number of first sampling points is greater than the number of second sampling points, then for each seismic trace, based on the multiple third sampling points of the seismic trace, the second signal values ​​of the multiple third sampling points of the seismic trace are found from the gain data. If the number of first sampling points is equal to the number of second sampling points, then for each seismic trace, based on the multiple second sampling points of the seismic trace, the second signal values ​​of the multiple second sampling points of the seismic trace are found from the gain data.

[0125] In some cases, W×d≥N, then the first N columns of the seismic profile data Y can be selected, that is, Y(M,1:N) is used as the gain result of the seismic profile data.

[0126] For example, see Figure 7 , Figure 7 : is a schematic diagram of a vertical seismic profile provided in an embodiment of the present application. The vertical seismic profile is obtained based on the signal value after gain. Figure 6 and Figure 7 It can be seen that Figure 7 The energy difference between seismic traces in the vertical seismic profile is significantly reduced.

[0127] In the above embodiment, the gain processing is performed by dividing each signal set into multiple subsets as an example. In this way, the gain processing is performed on the signal values ​​in each subset based on the signal values ​​in the subset, thereby improving the accuracy of the gain processing.

[0128] In other embodiments, the computer device directly performs gain processing on each signal set. Accordingly, the process of performing gain processing on the first signal values ​​of the second sampling points of multiple seismic channels to obtain the second signal values ​​of the second sampling points of the multiple seismic channels further includes the following implementation method: for each sampling time point, the computer device sequentially sorts the signal values ​​of the multiple seismic channels at the sampling time point to form a signal set; wherein, for each seismic channel, if the sampling time point is the second sampling point of the seismic channel, the signal value of the seismic channel at the sampling time point is the first signal value of the second sampling point; for each signal set, based on the multiple signal values ​​included in the signal set, the signal set is gain processed, and the second signal values ​​of the multiple second sampling points are determined from the multiple gain-processed signal sets corresponding to the multiple sampling time points. In this embodiment, based on the multiple signal values ​​in the signal set, the signal values ​​in the signal set are directly gain processed, thereby improving the efficiency of the gain processing.

[0129] Among them, the process of the computer device performing gain processing on the signal set based on the multiple signal values ​​included in the signal set is the same as the process of the above-mentioned computer device performing gain processing on the first subset based on the multiple signal values ​​included in the first subset, and will not be repeated here.

[0130] 306. The computer device performs an inverse transformation on the signal value of the first sampling point of each seismic trace based on the second signal value of the second sampling point of the seismic trace, and the signal values ​​of the first sampling points after the inverse transformation of the multiple seismic traces are used to determine the vertical seismic profile.

[0131] In this embodiment, an inverse transformation is performed based on the above formula (1), and the signal value of the first sampling point is transformed into the second signal value of the second sampling point. Similarly, if the number of first sampling points is less than the number of second sampling points, the signal values ​​of the plurality of first sampling points are transformed into the second signal values ​​of the plurality of target sampling points. If the number of first sampling points is greater than the number of second sampling points, the signal values ​​of the plurality of first sampling points are transformed into the second signal values ​​of the plurality of third sampling points. If the number of first sampling points is equal to the number of second sampling points, the signal values ​​of the plurality of first sampling points are transformed into the second signal values ​​of the plurality of second sampling points.

[0132] In an embodiment of the present application, after the signal value of the first sampling point is inversely transformed, the signal value of the second sampling point is restored to the signal value before the transformation, that is, zero value, thereby ensuring that the distribution of the signal value at multiple sampling time points is the same as the distribution before the seismic profile data is processed.

[0133] In an embodiment of the present application, the seismic profile data may be data in a vertical seismic profile, and the vertical seismic profile changes with changes in the seismic profile data. Accordingly, after performing an inverse transformation on the signal value at the first sampling point, a changed vertical seismic profile is obtained. The seismic profile data may also be data independent of the vertical seismic profile, and the vertical seismic profile does not change with changes in the seismic profile data. Accordingly, after performing an inverse transformation on the signal value at the first sampling point, the computer device generates the vertical seismic profile based on the inverse-transformed signal value at the first sampling point, based on an instruction for generating a vertical seismic profile.

[0134] For example, see Figure 8 and Figure 9 , Figure 8 and Figure 9 They are schematic diagrams of a vertical seismic profile provided by the embodiments of the present application. Figure 8 The vertical seismic profile in is obtained based on the method provided in the embodiment of the present application, Figure 9 The vertical seismic profile in the figure is not obtained by the method provided in the embodiment of the present application. As can be seen from the figure, Figure 9 In the area pointed by the arrows in the figure, there are obvious energy steps between adjacent seismic traces. Figure 8 In the corresponding area, the energy step problem has been significantly improved.

[0135] The method provided in the embodiment of the present application can perform automatic gain in multiple directions such as one-way time and two-way time, solving the problem of obvious energy steps displayed in vertical seismic profiles due to strong shielding layers.

[0136] In the embodiment of the present application, since signals are only acquired at the sampling time point corresponding to the time-series information and the sampling time points thereafter, the first sampling point is the valid sampling time point for acquiring signals, and the signal value at the second sampling point is zero. Furthermore, since the second sampling point occurs before the first sampling point, the signal value of the second sampling point is transformed based on the signal value of the first sampling point, thereby shifting the signal value of the first sampling point forward. Furthermore, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, thereby aligning the signals of the multiple seismic traces at the sampling time points. Gain processing is then performed on this basis. Since the signals of the multiple seismic traces are aligned based on the sampling time points, the regularity and convenience of the gain processing can be improved, thereby enhancing the gain processing effect. The signal value of the first sampling point is then inversely transformed, and the gain-processed signal is restored to the first sampling point, thus achieving effective gain processing of the seismic profile data. This solves the energy step problem in the generated vertical seismic profile, thereby improving the display effect of the vertical seismic profile.

[0137] The present application also provides a seismic data processing device. Figure 10, the device comprises:

[0138] An information determination module 1001 is configured to determine, from seismic profile data of an observation well, time history information of multiple seismic traces corresponding to multiple depth points in the observation well. The seismic profile data includes signal values ​​of the multiple seismic traces at multiple sampling time points. The time history information of the multiple seismic traces is the one-way time sampling time point or the two-way time sampling time point when the seismic wave reaches the seismic trace.

[0139] The sampling point determination module 1002 is configured to, for each seismic trace, determine, based on the time history information of the seismic trace, a sampling time point corresponding to the time history information and a sampling time point after the sampling time point corresponding to the time history information as a first sampling point of the seismic trace, and a sampling time point before the sampling time point corresponding to the time history information as a second sampling point of the seismic trace;

[0140] A transformation module 1003 is configured to transform, for each seismic trace, a signal value at a second sampling point of the seismic trace based on a signal value at the first sampling point of the seismic trace, to obtain a first signal value at the second sampling point of the seismic trace;

[0141] A gain module 1004 is configured to perform gain processing on the first signal values ​​of the second sampling points of the multiple seismic traces to obtain second signal values ​​of the second sampling points of the multiple seismic traces;

[0142] The inverse transformation module 1005 is used to perform an inverse transformation on the signal value of the first sampling point of each seismic trace based on the second signal value of the second sampling point of the seismic trace. The signal values ​​of the first sampling points of multiple seismic traces after inverse transformation are used to determine the vertical seismic profile.

[0143] In some embodiments, there are multiple first sampling points and multiple second sampling points, and the transformation module 1003 is configured to:

[0144] For each seismic trace, determining a plurality of target sampling points among the plurality of second sampling points based on the number of the plurality of first sampling points, wherein the number of the plurality of target sampling points is the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points is the starting sampling point among the plurality of second sampling points;

[0145] The signal values ​​of the plurality of target sampling points are transformed into the signal values ​​of the plurality of first sampling points to obtain the first signal values ​​of the plurality of target sampling points.

[0146] In some embodiments, the inverse transform module 1005 is configured to:

[0147] The signal values ​​of the plurality of first sampling points are transformed into the second signal values ​​of the plurality of target sampling points.

[0148] In some embodiments, the gain module 1004 is configured to:

[0149] For each sampling time point, the signal values ​​of multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point;

[0150] Dividing the signal set based on the reference depth window to obtain multiple first subsets;

[0151] For each first subset, performing gain processing on the first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset;

[0152] Second signal values ​​of second sampling points of the plurality of seismic traces are determined from the plurality of second subsets respectively corresponding to the plurality of sampling time points.

[0153] In some embodiments, the gain module 1004 is configured to:

[0154] For each first subset, determining a gain parameter based on the plurality of signal values ​​and the number of the plurality of signal values ​​in the first subset;

[0155] Products of a plurality of signal values ​​in the first subset and the gain parameter are respectively determined to obtain a second subset corresponding to the first subset.

[0156] In some embodiments, the gain module 1004 is configured to:

[0157] For each sampling time point, the signal values ​​of multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point;

[0158] For each signal set, gain processing is performed on the signal set based on multiple signal values ​​included in the signal set, and second signal values ​​of multiple second sampling points are determined from multiple gain-processed signal sets corresponding to multiple sampling time points respectively.

[0159] In the embodiment of the present application, since signals are only acquired at the sampling time point corresponding to the time-series information and the sampling time points thereafter, the first sampling point is the valid sampling time point for acquiring signals, and the signal value at the second sampling point is zero. Furthermore, since the second sampling point occurs before the first sampling point, the signal value of the second sampling point is transformed based on the signal value of the first sampling point, thereby shifting the signal value of the first sampling point forward. Furthermore, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, thereby aligning the signals of the multiple seismic traces at the sampling time points. Gain processing is then performed on this basis. Since the signals of the multiple seismic traces are aligned based on the sampling time points, the regularity and convenience of the gain processing can be improved, thereby enhancing the gain processing effect. The signal value of the first sampling point is then inversely transformed, and the gain-processed signal is restored to the first sampling point, thus achieving effective gain processing of the seismic profile data. This solves the energy step problem in the generated vertical seismic profile, thereby improving the display effect of the vertical seismic profile.

[0160] Figure 11 A structural block diagram of a computer device 1100 provided by an exemplary embodiment of the present application is shown.

[0161] Typically, the computer device 1100 includes a processor 1101 and a memory 1102 .

[0162] The processor 1101 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 1101 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), or PLA (Programmable Logic Array). The processor 1101 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 1101 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 1101 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0163] Memory 1102 may include one or more computer-readable storage media, which may be non-transitory. Memory 1102 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices and flash memory storage devices. In some embodiments, the non-transitory computer-readable storage medium in memory 1102 is used to store at least one computer program, which is executed by processor 1101 to implement the seismic data processing method provided in the method embodiment of the present application.

[0164] In some embodiments, computer device 1100 may optionally include a peripheral device interface 1103 and at least one peripheral device. Processor 1101, memory 1102, and peripheral device interface 1103 may be connected via a bus or signal lines. Each peripheral device may be connected to peripheral device interface 1103 via a bus, signal lines, or circuit boards. Specifically, the peripheral device may include at least one of a radio frequency circuit 1104, a display screen 1105, a camera assembly 1106, an audio circuit 1107, and a power supply 1108.

[0165] The peripheral device interface 1103 can be used to connect at least one I / O (Input / Output)-related peripheral device to the processor 1101 and the memory 1102. In some embodiments, the processor 1101, the memory 1102, and the peripheral device interface 1103 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 1101, the memory 1102, and the peripheral device interface 1103 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.

[0166] The RF circuit 1104 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 1104 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 1104 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals into electrical signals. Optionally, the RF circuit 1104 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, and the like. The RF circuit 1104 can communicate with other computer devices via at least one wireless communication protocol. Such wireless communication protocols include, but are not limited to, the World Wide Web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 1104 may also include circuitry related to Near Field Communication (NFC), which is not limited in this application.

[0167] The display screen 1105 is used to display a UI (User Interface). The UI may include graphics, text, icons, multimedia resources, and any combination thereof. When the display screen 1105 is a touch screen display, the display screen 1105 also has the ability to collect touch signals on the surface or above the surface of the display screen 1105. The touch signal can be input as a control signal to the processor 1101 for processing. In this case, the display screen 1105 can also be used to provide virtual buttons and / or virtual keyboards, also known as soft buttons and / or soft keyboards. In some embodiments, there can be one display screen 1105, which is set on the front panel of the computer device 1100; in other embodiments, there can be at least two display screens 1105, which are respectively set on different surfaces of the computer device 1100 or in a folding design; in other embodiments, the display screen 1105 can be a flexible display screen, which is set on the curved surface or folding surface of the computer device 1100. Even more, the display screen 1105 can be set to a non-rectangular irregular shape, that is, a special-shaped screen. The display screen 1105 can be made of materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).

[0168] The camera component 1106 is used to capture images or multimedia resources. Optionally, the camera component 1106 includes a front camera and a rear camera. Typically, the front camera is set on the front panel of the computer device, and the rear camera is set on the back of the computer device. In some embodiments, there are at least two rear cameras, which are any one of a main camera, a depth of field camera, a wide-angle camera, and a telephoto camera, so as to realize the fusion of the main camera and the depth of field camera to realize the background blur function, the fusion of the main camera and the wide-angle camera to realize panoramic shooting and VR (Virtual Reality) shooting function or other fusion shooting functions. In some embodiments, the camera component 1106 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash refers to a combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0169] The audio circuit 1107 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, and convert the sound waves into electrical signals that are input into the processor 1101 for processing, or input into the radio frequency circuit 1104 to achieve voice communication. For the purpose of stereo sound collection or noise reduction, there may be multiple microphones, each located in different parts of the computer device 1100. The microphone may also be an array microphone or an omnidirectional collection microphone. The speaker is used to convert electrical signals from the processor 1101 or the radio frequency circuit 1104 into sound waves. The speaker may be a traditional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can not only convert electrical signals into sound waves audible to humans, but also convert electrical signals into sound waves inaudible to humans for purposes such as distance measurement. In some embodiments, the audio circuit 1107 may also include a headphone jack.

[0170] Power supply 1108 is used to power various components in computer device 1100. Power supply 1108 can be AC ​​power, DC power, a disposable battery, or a rechargeable battery. When power supply 1108 includes a rechargeable battery, the rechargeable battery can be a wired rechargeable battery or a wireless rechargeable battery. A wired rechargeable battery is charged via a wired line, while a wireless rechargeable battery is charged via a wireless coil. The rechargeable battery can also support fast charging technology.

[0171] In some embodiments, the computer device 1100 further includes one or more sensors 1109 , including but not limited to an acceleration sensor 1110 , a gyroscope sensor 1111 , a pressure sensor 1112 , an optical sensor 1113 , and a proximity sensor 1114 .

[0172] The accelerometer 1110 can detect the magnitude of acceleration along the three coordinate axes of the coordinate system established by the computer device 1100. For example, the accelerometer 1110 can be used to detect the components of gravity acceleration along the three coordinate axes. The processor 1101 can control the display screen 1105 to display the user interface in a landscape or portrait view based on the gravity acceleration signal collected by the accelerometer 1110. The accelerometer 1110 can also be used to collect game or user motion data.

[0173] The gyroscope sensor 1111 can detect the orientation and rotation angle of the computer device 1100. It can also work with the accelerometer 1110 to collect 3D motions of the user on the computer device 1100. Based on the data collected by the gyroscope sensor 1111, the processor 1101 can implement the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.

[0174] The pressure sensor 1112 can be installed on the side frame of the computer device 1100 and / or below the display screen 1105. When the pressure sensor 1112 is installed on the side frame of the computer device 1100, it can detect the user's grip signal of the computer device 1100. The processor 1101 can perform left and right hand recognition or shortcut operations based on the grip signal collected by the pressure sensor 1112. When the pressure sensor 1112 is installed below the display screen 1105, the processor 1101 controls the operational controls on the UI interface based on the user's pressure operation on the display screen 1105. The operational controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0175] Optical sensor 1113 is used to detect ambient light intensity. In one embodiment, processor 1101 can control the display brightness of display screen 1105 based on the ambient light intensity detected by optical sensor 1113. Specifically, when the ambient light intensity is high, the display brightness of display screen 1105 is increased; when the ambient light intensity is low, the display brightness of display screen 1105 is decreased. In another embodiment, processor 1101 can also dynamically adjust the shooting parameters of camera assembly 1106 based on the ambient light intensity detected by optical sensor 1113.

[0176] Proximity sensor 1114, also known as a distance sensor, is typically located on the front panel of computer device 1100. Proximity sensor 1114 is used to detect the distance between the user and the front of computer device 1100. In one embodiment, when proximity sensor 1114 detects that the distance between the user and the front of computer device 1100 is gradually decreasing, processor 1101 controls display screen 1105 to switch from a screen-on state to a screen-off state. When proximity sensor 1114 detects that the distance between the user and the front of computer device 1100 is gradually increasing, processor 1101 controls display screen 1105 to switch from a screen-off state to a screen-on state.

[0177] Those skilled in the art will understand that Figure 11 The structure shown in the figure does not constitute a limitation on the computer device 1100, and the computer device 1100 may include more or fewer components than shown in the figure, or combine some components, or adopt a different component arrangement.

[0178] An embodiment of the present application also provides a computer-readable storage medium, in which at least one program code is stored. The at least one program code is loaded and executed by a processor to implement the seismic data processing method of any of the above implementation methods.

[0179] An embodiment of the present application also provides a computer program product, which includes computer program code, the computer program code is stored in a computer-readable storage medium, the processor of the computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code, so that the computer device executes the seismic data processing method of any of the above-mentioned implementation methods.

[0180] In some embodiments, the computer program product involved in the embodiments of the present application can be deployed and executed on a computer device, or on multiple computer devices located at one location, or on multiple computer devices distributed at multiple locations and interconnected through a communication network. Multiple computer devices distributed at multiple locations and interconnected through a communication network can constitute a blockchain system.

Claims

1. A seismic data processing method, characterized in that: The method comprises: Determining, from seismic profile data of an observation well, time history information of a plurality of seismic traces corresponding to a plurality of depth points in the observation well, the seismic profile data including signal values ​​of the plurality of seismic traces at a plurality of sampling time points, the time history information of the plurality of seismic traces being one-way time sampling time points or two-way time sampling time points when seismic waves arrive at the seismic traces; For each seismic trace, based on the time history information of the seismic trace, a sampling time point corresponding to the time history information and a sampling time point after the sampling time point corresponding to the time history information are used as the first sampling point of the seismic trace, and a sampling time point before the sampling time point corresponding to the time history information is used as the second sampling point of the seismic trace; For each seismic trace, transforming a signal value at a second sampling point of the seismic trace based on a signal value at a first sampling point of the seismic trace to obtain a first signal value at the second sampling point of the seismic trace; performing gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces respectively to obtain second signal values ​​of the second sampling points of the plurality of seismic traces; For each seismic trace, the signal value of the first sampling point of the seismic trace is inversely transformed based on the second signal value of the second sampling point of the seismic trace, and the signal values ​​of the first sampling points of multiple seismic traces after inverse transformation are used to determine the vertical seismic profile.

2. The method according to claim 1, characterized in that There are a plurality of first sampling points and a plurality of second sampling points, and for each seismic trace, transforming a signal value of a second sampling point of the seismic trace based on a signal value of the first sampling point of the seismic trace to obtain a first signal value of the second sampling point of the seismic trace includes: For each seismic trace, determining a plurality of target sampling points among a plurality of second sampling points based on the number of the plurality of first sampling points, wherein the number of the plurality of target sampling points is the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points is the starting sampling point among the plurality of second sampling points; The signal values ​​of the plurality of target sampling points are transformed into the signal values ​​of the plurality of first sampling points to obtain the first signal values ​​of the plurality of target sampling points.

3. The method according to claim 2, characterized in that For each seismic trace, performing an inverse transformation on a signal value of a first sampling point of the seismic trace based on a second signal value of a second sampling point of the seismic trace comprises: The signal values ​​of the plurality of first sampling points are transformed into second signal values ​​of the plurality of target sampling points.

4. The method according to claim 1, wherein The performing gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces to obtain the second signal values ​​of the second sampling points of the plurality of seismic traces includes: For each sampling time point, the signal values ​​of the multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point; Dividing the signal set based on a reference depth window to obtain a plurality of first subsets; For each first subset, performing gain processing on the first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset; Second signal values ​​of second sampling points of the plurality of seismic traces are determined from a plurality of second subsets respectively corresponding to the plurality of sampling time points.

5. The method according to claim 4, characterized in that The step of performing gain processing on each first subset based on multiple signal values ​​included in the first subset to obtain a second subset corresponding to the first subset includes: For each first subset, determining a gain parameter based on a plurality of signal values ​​in the first subset and a number of the plurality of signal values; Products of a plurality of signal values ​​in the first subset and the gain parameter are respectively determined to obtain a second subset corresponding to the first subset.

6. The method according to claim 1, characterized in that The performing gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces to obtain the second signal values ​​of the second sampling points of the plurality of seismic traces includes: For each sampling time point, the signal values ​​of the multiple seismic traces at the sampling time point are sequentially sorted to form a signal set; wherein, for each seismic trace, if the sampling time point is the second sampling point of the seismic trace, the signal value of the seismic trace at the sampling time point is the first signal value of the second sampling point; For each signal set, gain processing is performed on the signal set based on multiple signal values ​​included in the signal set, and second signal values ​​of the multiple second sampling points are determined from the multiple gain-processed signal sets corresponding to the multiple sampling time points respectively.

7. A seismic data processing device, characterized in that: The device comprises: an information determination module, configured to determine, from seismic profile data of an observation well, time history information of a plurality of seismic traces corresponding to a plurality of depth points in the observation well, the seismic profile data comprising signal values ​​of the plurality of seismic traces at a plurality of sampling time points, the time history information of the plurality of seismic traces being one-way time sampling time points or two-way time sampling time points when seismic waves arrive at the seismic traces; a sampling point determination module configured to, for each seismic trace, based on time history information of the seismic trace, determine a sampling time point corresponding to the time history information and a sampling time point after the sampling time point corresponding to the time history information as a first sampling point of the seismic trace, and determine a sampling time point before the sampling time point corresponding to the time history information as a second sampling point of the seismic trace; a transformation module, configured to transform, for each seismic trace, a signal value of a second sampling point of the seismic trace based on a signal value of a first sampling point of the seismic trace, to obtain a first signal value of the second sampling point of the seismic trace; a gain module, configured to perform gain processing on the first signal values ​​of the second sampling points of the plurality of seismic traces respectively to obtain second signal values ​​of the second sampling points of the plurality of seismic traces; An inverse transformation module is used to perform an inverse transformation on the signal value of the first sampling point of each seismic trace based on the second signal value of the second sampling point of the seismic trace, and the signal values ​​of the first sampling points of multiple seismic traces after inverse transformation are used to determine a vertical seismic profile.

8. A computer device, characterized in that: The computer device includes one or more processors and one or more memories, wherein at least one program code is stored in the one or more memories, and the at least one program code is loaded and executed by the one or more processors to implement the seismic data processing method according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that At least one program code is stored in the storage medium, and the at least one program code is loaded and executed by the processor to implement the seismic data processing method according to any one of claims 1 to 6.

10. A computer program product, characterized in that The computer program product includes a computer program code, which is stored in a computer-readable storage medium. The processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code, so that the computer device executes the seismic data processing method according to any one of claims 1 to 6.

Images