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

Abstract

The invention provides a seismic data processing method and device, computer equipment, a storage medium and a product, and belongs to the technical field of seismic exploration. The method comprises the following steps: for each seismic channel, based on time history information of the seismic channel, taking a sampling time point corresponding to the time history information and a sampling time point after the time history information as a first sampling point of the seismic channel, and taking a sampling time point before the sampling time point corresponding to the time history information as a second sampling point of the seismic channel; 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 a first signal value of 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 channels to obtain second signal values of the second sampling points of the plurality of seismic channels; and for each seismic trace, performing inverse transformation on the signal value of the first sampling point of the seismic trace based on the second signal value of the second sampling point of the seismic trace. The method improves the display effect of the vertical seismic profile map.

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 development of seismic exploration technology, there are more and more seismic observation methods, such as the ground seismic profile method and the vertical seismic profile (Vertical Seismic Profile, VSP) method. Among them, the ground seismic profile method is to collect seismic waves on the ground and generate a horizontal seismic profile based on adjacent seismic traces in the horizontal direction; the vertical seismic profile method is to collect seismic waves in the well and generate a vertical seismic profile based on adjacent seismic traces in the depth direction. Among them, in the horizontal direction, since the surface conditions between adjacent traces are close and the excitation conditions are consistent, there is no obvious difference in energy between adjacent traces in the horizontal seismic profile. In the vertical seismic profile, since the distance between adjacent traces in the depth direction is large, when there is a thin energy shielding layer (such as igneous rock) in the underground stratum, the energy difference between adjacent traces will be obvious, which reduces the display effect of the vertical seismic profile. Therefore, a seismic data processing method is needed to process the seismic profile data of the vertical seismic profile to reduce the energy difference between seismic traces, thereby improving the display effect of the vertical seismic profile. Summary of the invention

[0003] The embodiments of the present application provide a seismic data processing method, device, computer equipment, storage medium and product, 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] Determine, from the seismic profile data of the observation well, time history information of multiple seismic traces corresponding to multiple depth points in the observation well, wherein the seismic profile data includes signal values ​​of the multiple seismic traces at multiple sampling time points, and the time history information of the multiple seismic traces is a one-way time sampling time point or a two-way time sampling time point when the seismic wave reaches the seismic trace;

[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, based on the signal value of the first sampling point of the seismic trace, transform the signal value of the second sampling point of the seismic trace to obtain the first signal value of 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, 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, including:

[0011] For each seismic trace, based on the number of the plurality of first sampling points, determine a plurality of target sampling points among the plurality of second sampling points, the number of the plurality of target sampling points being the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points being 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 multiple seismic traces to obtain the second signal values ​​of the second sampling points of the multiple 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] Based on the reference depth window, the signal set is divided to obtain a plurality of first subsets;

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

[0019] The second signal values ​​of the 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.

[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 multiple seismic traces to obtain the second signal values ​​of the second sampling points of the multiple 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 respectively corresponding to the multiple sampling time points.

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

[0027] An information determination module, for determining, from the seismic profile data of the 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 a one-way time sampling time point or a two-way time sampling time point when the seismic wave reaches the seismic trace;

[0028] A sampling point determination module is used for, for each seismic trace, based on the time history information of the seismic trace, taking the sampling time point corresponding to the time history information and the sampling time point after the sampling time point corresponding to the time history information as the first sampling point of the seismic trace, and taking the sampling time point before the sampling time point corresponding to the time history information as the 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, used for 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;

[0031] The 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 the 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, based on the number of the plurality of first sampling points, determine a plurality of target sampling points among the plurality of second sampling points, the number of the plurality of target sampling points being the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points being 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 used 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 used 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] Based on the reference depth window, the signal set is divided to obtain a plurality of first subsets;

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

[0041] The second signal values ​​of the 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.

[0042] In some embodiments, the gain module is used 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 used 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 respectively corresponding to the multiple sampling time points.

[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-mentioned 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, wherein 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 implementations.

[0051] In the embodiment of the present application, since the signal is collected only at the sampling time point corresponding to the time history information and the sampling time point thereafter, that is, the first sampling point is the effective sampling time point for collecting the signal, and the signal value of the second sampling point is zero. And since the time of the second sampling point is 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, so that the signal value of the first sampling point is shifted forward. Then, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, so that the alignment of the signals of multiple seismic traces at the sampling time point is achieved. On this basis, gain processing is performed again. Since the signals of multiple seismic traces are aligned based on the sampling time point, the regularity and convenience of the gain processing can be improved, thereby improving the gain processing effect. Then, the signal value of the first sampling point is inversely transformed, and the signal after gain processing is restored to the first sampling point, that is, effective gain processing of the seismic profile data is achieved, and the energy step problem in the generated vertical seismic profile is solved, 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 drawings required for use in the description of the embodiments will be briefly introduced below. 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 creative work.

[0053] Figure 1 It 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] Fig. 9 is a schematic diagram of another vertical seismic profile provided in an embodiment of the present application;

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

[0063] Fig.11 It 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 the present application clearer, the implementation methods of the present application will be further described in detail below in conjunction with the accompanying drawings.

[0065] The terms "first", "second", "third" and "fourth" etc. in the specification and claims of the present application and the drawings are used to distinguish different objects, rather than 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 device that includes a series of steps or units is not limited to the listed steps or units, but optionally includes steps or units that are not listed, or optionally includes other steps or units inherent to these processes, methods, products or devices.

[0066] The above are only optional embodiments of the present application and are not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application shall be included in the protection scope 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 obtained with full authorization.

[0067] Figure 1is a schematic diagram of an implementation environment of a seismic data processing method provided by the present 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, a user can log in to the target application, and then process the seismic profile data through 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 is installed with a target application, and the server 102 is a server 102 corresponding to the target application. The user processes the seismic profile data through the target application, and the server 102 provides background services. In the embodiment of the present application, only the computer device provided as a terminal 101 is used 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, a vehicle-mounted terminal, etc. The server 102 is at least one of a server, a server cluster consisting of multiple servers, a cloud server, a cloud computing platform, and a virtualization center.

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

[0072] 201. The computer equipment determines the time course 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 course information of the multiple seismic channels is the one-way sampling time point or the two-way sampling time point when the seismic wave arrives at the seismic channel.

[0073] In the embodiment of the present application, the signal value is an amplitude value. The computer device performs first arrival wave picking processing and dynamic correction processing on the seismic data to obtain seismic profile data. The one-way time sampling time point of each seismic trace refers to the sampling time point when the seismic wave first arrives at the seismic trace before dynamic correction. The two-way time sampling time point of each seismic trace refers to the sampling time point when the seismic wave first arrives at the seismic trace after dynamic 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, the seismic profile data is the seismic profile data before the dynamic correction, and the computer device processes the seismic profile data before the dynamic correction accordingly. 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, the seismic profile data is the seismic profile data after the dynamic correction, and the computer device processes the seismic profile data after the dynamic correction accordingly.

[0075] In the embodiment of the present application, the one-way time sampling point and the two-way time sampling point are respectively sampling time points among a plurality of sampling time points. A plurality of sampling time points may be distributed in sequence at fixed time intervals, or the time intervals may be set as needed. A plurality of seismic traces may be distributed in sequence at fixed depth intervals, or the depth intervals may 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 from each seismic trace in the seismic profile data before the 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 from each seismic trace in the seismic profile data after the dynamic correction is used as the time history information of each seismic trace.

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

[0077] In the embodiment of the present application, for each seismic trace, no valid signal is collected at the sampling time point 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 moved forward.

[0080] 204. The computer device performs gain processing on the first signal values ​​of the second sampling points of the multiple seismic channels respectively 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 of the multiple seismic traces after the inverse transformation 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 realizing the backward shift of the sampling time point of the second signal value of the second sampling point.

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

[0085] In the embodiment of the present application, since the signal is collected only at the sampling time point corresponding to the time history information and the sampling time point thereafter, that is, the first sampling point is the effective sampling time point for collecting the signal, and the signal value of the second sampling point is zero. And since the time of the second sampling point is 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, so that the signal value of the first sampling point is shifted forward. Then, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, so that the alignment of the signals of multiple seismic traces at the sampling time point is achieved. On this basis, gain processing is performed again. Since the signals of multiple seismic traces are aligned based on the sampling time point, the regularity and convenience of the gain processing can be improved, thereby improving the gain processing effect. Then, the signal value of the first sampling point is inversely transformed, and the signal after gain processing is restored to the first sampling point, that is, effective gain processing of the seismic profile data is achieved, and the energy step problem in the generated vertical seismic profile is solved, 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 describes the process of seismic data processing. Figure 3 , Figure 3 It is a flowchart of a seismic profile data processing method provided in an embodiment of the present application, and the method includes the following steps.

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

[0088] In the embodiment of the present application, when collecting seismic data, the position of the reference distance from the wellhead of the observation well on the surface is determined as the source excitation point, and the surface source is used to trigger seismic waves, and the seismic waves are received by the three-component detector in the observation well to obtain seismic data. The observation well is a vertical well or a newly drilled well, an old well, a reconstructed well, etc. that is approximately a vertical well. The reference distance between the source excitation point and the observation well meets the zero well source distance excitation distance requirement in the industry standard SY / T 7450, so as to take into account both the well site operation safety and the drilling safety.

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

[0090] 302. The computer equipment determines the time course 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 course information of the multiple seismic channels is the one-way sampling time point or the two-way sampling time point when the seismic wave arrives at 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 respectively 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.

[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 and multiple seismic traces in (a) 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, based on the time information of the seismic trace, uses the sampling time point corresponding to the time information and the sampling time point after the sampling time point corresponding to the time information as the first sampling point of the seismic trace, and uses the 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 the 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 the number of second sampling points, may be greater than the number of second sampling points, or may be equal to the number of second sampling points. Accordingly, the process in which the above-mentioned 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 for each seismic trace to obtain the first signal value of the second sampling point of the seismic trace includes the following implementation methods.

[0097] The 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 multiple target sampling points among the multiple second sampling points based on the number of multiple first sampling points. The number of multiple target sampling points is the same as the number of multiple first sampling points, and the starting sampling point of the multiple target sampling points is the starting sampling point among the multiple second sampling points; the computer device transforms the signal values ​​of the multiple target sampling points into the signal values ​​of the multiple first sampling points to obtain the first signal values ​​of the multiple target sampling points.

[0098] Among them, the starting sampling point in the second sampling point is also the starting sampling time point in the multiple sampling time points, and the time of the starting sampling time point is the earliest. Optionally, the starting sampling time point is zero. In this implementation, the smallest sampling time point among the multiple second sampling points is used as the starting point, and the corresponding number of second sampling points are intercepted backward to obtain multiple target sampling points with the same number as the first sampling points. The multiple first sampling points are arranged in time order, and the multiple target sampling points are also arranged in time order, and then the signal value of each target sampling point is transformed into the signal value of the first sampling point corresponding to the order. For example, the signal value of the first target sampling point is transformed into the signal value of the first first sampling point, and the signal value of the second target sampling point is transformed into the signal value of the second first sampling point, etc. Optionally, the signal values ​​of the second sampling points other than the target sampling points in 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 transformed into signal values ​​of multiple first sampling points, so that the signals of the first sampling points are flattened to the starting sampling point, and then the signals of multiple seismic channels can be aligned at the starting sampling point, which facilitates the subsequent gain processing of the signals of multiple seismic channels at the same time.

[0100] The 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, the plurality of third sampling points including the plurality of second sampling points and the first number of first sampling points, the sum of the first number and the number of the plurality of second sampling points is equal to the number of the plurality of first sampling points, and the starting sampling point of the first number of first sampling points is the starting sampling point of the plurality of first sampling points. The signal values ​​of the plurality of third 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 third sampling points.

[0101] In this implementation, the minimum sampling time point among the plurality of first sampling points is taken as the starting point, and the corresponding number of first sampling points are intercepted backward to obtain a plurality of third sampling points with the same number as the first sampling points. The plurality of first sampling points are arranged in time sequence, and the plurality of third sampling points are also arranged in time sequence, and then the signal value of each third sampling point is 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, and the signal value of the second third sampling point is transformed into the signal value of the second first sampling point, etc. In this implementation, when the number of the first sampling points is greater than the number of the second sampling points, the signal values ​​of some of the first sampling points are also transformed based on the signal values ​​of the first sampling points, so as to ensure the integrity of the forward shift of the signal values ​​of the first sampling points, that is, to improve the precision and accuracy of the data, and facilitate the subsequent gain processing.

[0102] The 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 transforms the signal values ​​of multiple second sampling points into the signal values ​​of multiple first sampling points to obtain the first signal values ​​of multiple second sampling points. Multiple first sampling points are arranged in time sequence, and multiple second sampling points are also arranged in time sequence, and then the signal value of each second sampling point is transformed into the signal value of the first sampling point corresponding to the sequence. For example, the signal value of the first second sampling point is transformed into the signal value of the first first sampling point, and the signal value of the second second sampling point is transformed into the signal value of the second first sampling point, etc. In this implementation method, the transformation is performed directly, which improves the efficiency of the transformation.

[0103] It should be noted that, since the signal value of the first sampling point is subsequently inversely transformed based on the second signal value of the second sampling point, after the signal value of the second sampling point is transformed, the signal value of the first sampling point may remain unchanged or become zero. In the embodiment of the present application, the signal value of the first sampling point becomes zero as an example for explanation, so as to facilitate subsequent gain processing.

[0104] For example, the number of the multiple sampling time points is M, the number of the multiple seismic channels is N, and M and N are integers greater than 0, respectively. Then the seismic profile data before transformation can be expressed as A(M, N). Accordingly, the computer device can perform 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 represents the kth sampling time point to the k+p-1th sampling time 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 : 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 the signal value after double-path time information transformation, 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 respectively 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 multiple seismic channels, including 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; 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 the 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 the 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 the 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 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 the Nth seismic traces at that sampling time point.

[0115] The computer device sends the 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] Wherein, Depth represents the depth of the measured observation well; ΔDep represents the depth interval, i.e., the reference depth window; Int represents the rounding function, for example, Int(5 / 3)=1; Cdep represents the signal acquisition depth interval. Wherein, the reference depth window is a sliding depth window. The computer device can determine the reference depth window according to the law of formation change. For example, the reference depth window is determined based on the formation lithology. Wherein, 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, the measured depth of the observation well is 8,000 meters, and the reference depth window is 400 meters, then 20 first subsets are obtained. The signal acquisition depth interval in each first subset is 20 meters, then a first subset includes 20 signal values.

[0118] In some embodiments, the above-mentioned computer device performs 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, 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 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 embodiment 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 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 the ith element, X iw (l) represents the lth signal value in the first subset.

[0123] Among them, 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 arranges the multiple second subsets corresponding to the sampling time point in sequence according to the depth order to obtain the multiple signal values ​​after the gain of 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 the first sampling points is less than the number of the 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 the first sampling points is greater than the number of the 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 the first sampling points is equal to the number of the 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, so that 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 in which the computer device performs 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 multiple seismic channels also includes the following implementation method: for each sampling time point, the computer device sequentially sorts the signal values ​​of multiple seismic channels at the sampling time points 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 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 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 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 of the multiple seismic traces after the inverse transformation 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 the first sampling points is less than the number of the 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 the first sampling points is greater than the number of the 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 the first sampling points is equal to the number of the 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 the signal value of the first sampling point is inversely transformed, the 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 the signal value of the first sampling point is inversely transformed, the computer device generates the vertical seismic profile based on the vertical seismic profile generation instruction and the signal value of the first sampling point after the inverse transformation.

[0134] For example, see Figure 8 and Fig. 9 , Figure 8 and Fig. 9 They are schematic diagrams of a vertical seismic profile provided in 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, Fig. 9 The vertical seismic profile in is not obtained by the method provided in the embodiment of the present application. It can be seen from the figure that Fig. 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, thereby 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 the signal is collected only at the sampling time point corresponding to the time history information and the sampling time point thereafter, that is, the first sampling point is the effective sampling time point for collecting the signal, and the signal value of the second sampling point is zero. And since the time of the second sampling point is 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, so that the signal value of the first sampling point is shifted forward. Then, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, so that the alignment of the signals of multiple seismic traces at the sampling time point is achieved. On this basis, gain processing is performed again. Since the signals of multiple seismic traces are aligned based on the sampling time point, the regularity and convenience of the gain processing can be improved, thereby improving the gain processing effect. Then, the signal value of the first sampling point is inversely transformed, and the signal after gain processing is restored to the first sampling point, that is, effective gain processing of the seismic profile data is achieved, and the energy step problem in the generated vertical seismic profile is solved, thereby improving the display effect of the vertical seismic profile.

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

[0138] The information determination module 1001 is used to determine the time history information of multiple seismic traces corresponding to multiple depth points in the observation well from the seismic profile data of the observation well, wherein the seismic profile data includes signal values ​​of the multiple seismic traces at multiple sampling time points, and 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 used for, for each seismic trace, based on the time history information of the seismic trace, taking the sampling time point corresponding to the time history information and the sampling time point after the sampling time point corresponding to the time history information as the first sampling point of the seismic trace, and taking the sampling time point before the sampling time point corresponding to the time history information as the second sampling point of the seismic trace;

[0140] A transformation module 1003 is used for transforming 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;

[0141] A gain module 1004, 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;

[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 used to:

[0144] For each seismic trace, based on the number of the plurality of first sampling points, determining a plurality of target sampling points among the plurality of second sampling points, the number of the plurality of target sampling points being the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points being 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 used 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 used 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] Based on the reference depth window, the signal set is divided to obtain a plurality of first subsets;

[0151] For each first subset, based on a plurality of signal values ​​included in the first subset, perform gain processing on 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 used 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] The products of the 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 used 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 the signal is collected only at the sampling time point corresponding to the time history information and the sampling time point thereafter, that is, the first sampling point is the effective sampling time point for collecting the signal, and the signal value of the second sampling point is zero. And since the time of the second sampling point is 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, so that the signal value of the first sampling point is shifted forward. Then, the signal values ​​of the first sampling points of multiple seismic traces are uniformly shifted forward, so that the alignment of the signals of multiple seismic traces at the sampling time point is achieved. On this basis, gain processing is performed again. Since the signals of multiple seismic traces are aligned based on the sampling time point, the regularity and convenience of the gain processing can be improved, thereby improving the gain processing effect. Then, the signal value of the first sampling point is inversely transformed, and the signal after gain processing is restored to the first sampling point, that is, effective gain processing of the seismic profile data is achieved, and the energy step problem in the generated vertical seismic profile is solved, thereby improving the display effect of the vertical seismic profile.

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

[0164] In some embodiments, the computer device 1100 may also optionally include: a peripheral device interface 1103 and at least one peripheral device. The processor 1101, the memory 1102 and the peripheral device interface 1103 may be connected via a bus or a signal line. Each peripheral device may be connected to the peripheral device interface 1103 via a bus, a signal line or a circuit board. Specifically, the peripheral device includes: 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 may be used to connect at least one peripheral device related to I / O (Input / Output) 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 may be implemented on a separate chip or circuit board, which is not limited in this embodiment.

[0166] The radio frequency circuit 1104 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 1104 communicates with the communication network and other communication devices through electromagnetic signals. The radio frequency circuit 1104 converts the electrical signal into an electromagnetic signal for transmission, or converts the received electromagnetic signal into an electrical signal. Optionally, the radio frequency 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 radio frequency circuit 1104 can communicate with other computer devices through at least one wireless communication protocol. The wireless communication protocol includes, but is not limited to: the World Wide Web, a metropolitan area network, an intranet, various generations of mobile communication networks (2G, 3G, 4G and 5G), a wireless local area network and / or a WiFi (Wireless Fidelity) network. In some embodiments, the radio frequency circuit 1104 may also include circuits related to NFC (Near Field Communication), 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 display screen, 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 to the processor 1101 as a control signal for processing. At this time, 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, the display screen 1105 can be one, arranged on the front panel of the computer device 1100; in other embodiments, the display screen 1105 can be at least two, respectively arranged 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, arranged on a curved surface or a folding surface of the computer device 1100. Even, the display screen 1105 can also be arranged as a non-rectangular irregular figure, 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 assembly 1106 is used to collect images or multimedia resources. Optionally, the camera assembly 1106 includes a front camera and a rear camera. Usually, 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 assembly 1106 may also include a flash. The flash can be a monochrome 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 at 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 and input them into the processor 1101 for processing, or input them into the radio frequency circuit 1104 to achieve voice communication. For the purpose of stereo acquisition or noise reduction, there may be multiple microphones, which are respectively arranged at different parts of the computer device 1100. The microphone may also be an array microphone or an omnidirectional acquisition microphone. The speaker is used to convert the electrical signal from the processor 1101 or the radio frequency circuit 1104 into sound waves. The speaker may be a traditional film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can not only convert the electrical signal into sound waves audible to humans, but also convert the electrical signal into sound waves inaudible to humans for purposes such as ranging. In some embodiments, the audio circuit 1107 may also include a headphone jack.

[0170] The power supply 1108 is used to power various components in the computer device 1100. The power supply 1108 can be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the 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 a battery that is charged through a wired line, and a wireless rechargeable battery is a battery that is charged through a wireless coil. The rechargeable battery can also be used to 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 acceleration sensor 1110 can detect the magnitude of acceleration on the three coordinate axes of the coordinate system established by the computer device 1100. For example, the acceleration sensor 1110 can be used to detect the components of gravity acceleration on the three coordinate axes. The processor 1101 can control the display screen 1105 to display the user interface in a horizontal view or a vertical view according to the gravity acceleration signal collected by the acceleration sensor 1110. The acceleration sensor 1110 can also be used to collect motion data of games or users.

[0173] The gyro sensor 1111 can detect the body direction and rotation angle of the computer device 1100, and the gyro sensor 1111 can cooperate with the acceleration sensor 1110 to collect the user's 3D actions on the computer device 1100. The processor 1101 can implement the following functions based on the data collected by the gyro sensor 1111: motion sensing (such as changing the UI according to the user's tilt operation), image stabilization during shooting, game control, and inertial navigation.

[0174] The pressure sensor 1112 may be provided at the side frame of the computer device 1100 and / or the lower layer of the display screen 1105. When the pressure sensor 1112 is provided at the side frame of the computer device 1100, it may detect the user's gripping signal of the computer device 1100, and the processor 1101 may perform left and right hand recognition or shortcut operation according to the gripping signal collected by the pressure sensor 1112. When the pressure sensor 1112 is provided at the lower layer of the display screen 1105, the processor 1101 may control the operability controls on the UI interface according to the user's pressure operation on the display screen 1105. The operability controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0175] The optical sensor 1113 is used to collect the ambient light intensity. In one embodiment, the processor 1101 can control the display brightness of the display screen 1105 according to the ambient light intensity collected by the optical sensor 1113. Specifically, when the ambient light intensity is high, the display brightness of the display screen 1105 is increased; when the ambient light intensity is low, the display brightness of the display screen 1105 is reduced. In another embodiment, the processor 1101 can also dynamically adjust the shooting parameters of the camera assembly 1106 according to the ambient light intensity collected by the optical sensor 1113.

[0176] The proximity sensor 1114, also called a distance sensor, is usually disposed on the front panel of the computer device 1100. The proximity sensor 1114 is used to collect the distance between the user and the front of the computer device 1100. In one embodiment, when the proximity sensor 1114 detects that the distance between the user and the front of the computer device 1100 is gradually decreasing, the processor 1101 controls the display screen 1105 to switch from the screen-on state to the screen-off state; when the proximity sensor 1114 detects that the distance between the user and the front of the computer device 1100 is gradually increasing, the processor 1101 controls the display screen 1105 to switch from the screen-off state to the screen-on state.

[0177] Those skilled in the art will understand that Fig.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 less 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 modes.

[0179] An embodiment of the present application also provides a computer program product, 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 of any of the above-mentioned implementation methods.

[0180] In some embodiments, the computer program product involved in the embodiments of the present application may 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 by a communication network. Multiple computer devices distributed at multiple locations and interconnected by a communication network may constitute a blockchain system.

Claims

1. A method for processing seismic data, characterized in that: The method comprises: Determine, from the seismic profile data of the observation well, time history information of multiple seismic traces corresponding to multiple depth points in the observation well, wherein the seismic profile data includes signal values ​​of the multiple seismic traces at multiple sampling time points, and the time history information of the multiple seismic traces is a one-way time sampling time point or a two-way time sampling time point when the seismic wave reaches the seismic trace; 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, based on the signal value of the first sampling point of the seismic trace, transform the signal value of the second sampling point of the seismic trace to obtain the first signal value of the second sampling point of the seismic trace; Performing gain processing on the first signal values ​​of the second sampling points of the multiple seismic traces respectively to obtain the second signal values ​​of the second sampling points of the multiple 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, based on the signal value of the first sampling point of the seismic trace, transforming the signal value of the second sampling point of the seismic trace to obtain the first signal value of the second sampling point of the seismic trace comprises: For each seismic trace, based on the number of the plurality of first sampling points, determine a plurality of target sampling points among the plurality of second sampling points, the number of the plurality of target sampling points being the same as the number of the plurality of first sampling points, and a starting sampling point among the plurality of target sampling points being 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, based on a second signal value of a second sampling point of the seismic trace, inversely transforming a signal value of a first 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, characterized in that: The step of 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 comprises: 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; Based on the reference depth window, the signal set is divided to obtain a plurality of first subsets; For each first subset, based on a plurality of signal values ​​included in the first subset, perform gain processing on the first subset to obtain a second subset corresponding to the first subset; The second signal values ​​of the 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.

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 step of 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 comprises: 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 respectively corresponding to the multiple sampling time points.

7. A seismic data processing device, characterized in that: The device comprises: An information determination module, for determining, from the seismic profile data of the 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 a one-way time sampling time point or a two-way time sampling time point when the seismic wave reaches the seismic trace; A sampling point determination module is used for, for each seismic trace, based on the time history information of the seismic trace, taking the sampling time point corresponding to the time history information and the sampling time point after the sampling time point corresponding to the time history information as the first sampling point of the seismic trace, and taking the sampling time point before the sampling time point corresponding to the time history information as the 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, used for 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; The 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 the 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 as described in 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. 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 performs the seismic data processing method as described in any one of claims 1 to 6.

Images