Earthquake first arrival automatic picking method, device, equipment and storage medium

Abstract

The application provides a seismic first arrival automatic picking method, device and equipment and a storage medium, and belongs to the field of oil and gas geophysical prospecting engineering. In the application, the double sliding window method is used to detect the first arrival take-off time in the seismic first arrival signal, wherein two adjacent windows are relatively static in the process of sliding relative to the received signal. In this method, the first arrival is basically reliable in the first picking process. Then, based on the picked first arrival, a linear smoothing intercept window is picked, the first arrival travel time in the window is fitted with a first arrival function, so that a relatively accurate first arrival travel time can be obtained to realize automatic picking, so that manual intervention is not required, the workload of the staff is reduced, the timeliness is improved, and the seismic wave first arrival picking effect has the characteristics of high precision and strong noise resistance. Even in an area with low signal-to-noise ratio, good results can be achieved.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of oil and gas geophysical prospecting engineering, in particular to a method and device for automatically picking up seismic first arrival, equipment and storage medium. BACKGROUND

[0002] Currently, the static correction problem in seismic exploration has been puzzling experts and scholars in the field of seismic exploration for accurate exploration and well placement in western, southern mountainous and loess plateau areas with dramatic surface relief and complex surface structure. Currently, the static correction processing of seismic data in the relief surface area is still mainly based on elevation static correction, refraction static correction and tomographic static correction methods. However, the effect of such static correction methods is mostly dependent on the picking accuracy of seismic first arrival travel time; in VSP (vertical seismic profile) seismic data processing, accurate seismic first arrival picking is the basis for all VSP data processing and interpretation. The accurate picking of first arrival travel time not only affects the effect of static correction, but also is crucial for the accuracy of shallow velocity modeling. Accurate shallow velocity modeling can be applied not only in the oil exploration industry, but also in natural earthquakes, water and soil resources, and engineering geophysical prospecting industries. Conventional near-surface first arrival travel time tomography can use first arrival time difference to invert shallow velocity, but it needs to input accurate seismic first arrival time. In addition, the early arrival waveform inversion proposed in recent years is also a relatively accurate shallow seismic velocity inversion method, but it still needs to input accurate first arrival travel time for early arrival waveform interception.

[0003] When facing complex surface structure, dramatic terrain relief, high-speed rock outcrop and gravel layer, how to accurately pick up the first arrival is an important problem. Currently, there are many automatic picking methods for seismic first arrival, such as algorithms based on energy, frequency, waveform similarity, vibration trajectory, polarization characteristics, etc. However, most of them need manual intervention, which leads to huge workload, long data processing period and inability to achieve timeliness.

[0004] Therefore, how to solve the above problems is an urgent problem to be solved at present. SUMMARY

[0005] The embodiments of the present application provide a method, device, equipment and storage medium for automatically picking up seismic first arrival, which can reduce workload and improve timeliness.

[0006] To achieve the above purpose, the present application adopts the following technical solutions:

[0007] In a first aspect, a method for automatically picking up seismic first arrival is provided. The method for automatically picking up seismic first arrival comprises:

[0008] The first arrival time in the seismic first arrival signal in the seismic data to be processed is acquired based on a double sliding window method;

[0009] The first arrival time is filtered based on a filter operator to obtain a filtered travel time;

[0010] The seismic data to be processed is time windowed based on the filtered travel time to obtain first arrival travel time information;

[0011] The first arrival travel time information is fitted with a first arrival function to obtain a first arrival time result;

[0012] The seismic data to be processed is cut based on the first arrival time result and the time window, and a target first arrival time is output.

[0013] According to the method of the first aspect, the first arrival time in the seismic first arrival signal in the seismic data to be processed is acquired based on a double sliding window method, and two adjacent windows are relatively static in the process of sliding relative to the received signal. In this method, the first arrival is basically reliable in the first picking process. Then, the picked first arrival is used to linearly smooth the cutting window, and the first arrival travel time in the window is fitted with a first arrival function. Therefore, the accurate first arrival travel time can be obtained, the picking is automated, manual intervention is not required, the workload of the staff is reduced, the timeliness is improved, and the first arrival picking effect of the seismic wave has the characteristics of high precision and strong noise resistance, and good results can be achieved even in an area with low signal-to-noise ratio.

[0014] In a possible design scheme, before the first arrival time in the seismic first arrival signal in the seismic data to be processed is acquired based on the double sliding window method, the method further includes: pixelating original seismic data to obtain the seismic data to be processed.

[0015] According to the above design scheme, by pixelating the original seismic data before the data is processed by the double sliding window method, the data processing efficiency of the double sliding window method can be effectively improved, and the accuracy can be improved.

[0016] In a possible design scheme, the first arrival time in the seismic first arrival signal in the seismic data to be processed is acquired based on the double sliding window method, including: detecting the first arrival takeoff time in the seismic first arrival signal in the seismic data to be processed by the double sliding window method to obtain the first arrival time, wherein the first arrival takeoff time is the time when the energy ratio starts to decrease when the seismic first arrival signal starts to enter the next time window.

[0017] According to the above design scheme, the first arrival takeoff time when the energy ratio starts to decrease is acquired by using the double sliding window method, so that the accurate first arrival time can be obtained, and the accuracy and precision of picking the first arrival of the earthquake are improved.

[0018] In a possible design, the filtering the first arrival time based on the filter operator to obtain a filtered travel time comprises: processing a seismic waveform picked up from the to-be-processed seismic data based on a one-dimensional Gaussian filter function to obtain the filter operator; and filtering the first arrival time based on the filter operator to obtain a filtered travel time.

[0019] According to the design, the one-dimensional Gaussian filter function is used to process the seismic waveform picked up from the to-be-processed seismic data, so that the filter operator is obtained, and the first arrival time is filtered based on the filter operator, so that the filtered travel time is more accurate, and the accuracy of picking up the first arrival time of the earthquake is improved.

[0020] In a possible design,

[0021] wherein T(itrace) represents a seismic waveform picked up from the to-be-processed seismic data of the ith trace, represents a one-dimensional Gaussian smoothing filter function.

[0022] In a possible design, the first arrival function fitting is performed on the first arrival time information to obtain a first arrival time result, comprising: performing first arrival function fitting on the first arrival time information based on a Biharmonic spline interpolation algorithm to obtain a first arrival time function; and outputting a first arrival time result according to the first arrival time function.

[0023] According to the design, the Biharmonic spline interpolation algorithm is used to perform linear smoothing windowing on the picked first arrival time information, and the first arrival time in the window is polynomially interpolated, so that the first arrival time function is accurately reconstructed, and more accurate first arrival time is obtained, thereby improving the accuracy of picking up the first arrival of the earthquake.

[0024] In a possible design, the first arrival time function satisfies

[0025]

[0026] wherein, is a Biharmonic operator, and δ is a Dirac function.

[0027] In a second aspect, a device for automatically picking up a first arrival of an earthquake is provided. The device for automatically picking up a first arrival of an earthquake comprises:

[0028] A first arrival time calculation module is configured to obtain a first arrival time in a first arrival signal of to-be-processed seismic data based on a double sliding window method.

[0029] a filter travel time calculation module, configured to filter the first arrival travel time based on a filter operator to obtain a filter travel time;

[0030] a second first arrival travel time calculation module, configured to perform time window interception on the to-be-processed seismic data according to the filter travel time to obtain first arrival travel time information;

[0031] a third first arrival travel time calculation module, configured to perform first arrival function fitting on the first arrival travel time information to obtain a first arrival travel time result;

[0032] an output module, configured to perform interception on the to-be-processed seismic data according to the first arrival travel time result and the time window, and output a target first arrival travel time.

[0033] In a possible design scheme, the first first arrival travel time calculation module is specifically configured to: obtain a first arrival travel time in a seismic first arrival signal in the to-be-processed seismic data by a double sliding window method, wherein the first arrival travel time is a time corresponding to when an energy ratio starts to drop when the seismic first arrival signal starts to enter a next time window.

[0034] In a possible design scheme, the seismic first arrival automatic picking device further includes a data processing module, configured to perform pixelization on original seismic data to obtain to-be-processed seismic data before obtaining a first arrival travel time in a seismic first arrival signal in the to-be-processed seismic data based on the double sliding window method.

[0035] In a possible design scheme, the filter travel time calculation module is specifically configured to: process a seismic waveform picked up from the to-be-processed seismic data based on a one-dimensional Gaussian filter function to obtain the filter operator; and filter the first arrival travel time based on the filter operator to obtain a filter travel time.

[0036] In a possible design scheme, the filter operator satisfies:

[0037]

[0038] wherein T(itrace) represents a seismic waveform picked up from the to-be-processed seismic data of an ithtrace channel, represents a one-dimensional Gaussian smoothing filter function.

[0039] In a possible design scheme, the third first arrival travel time calculation module is specifically configured to: perform first arrival function fitting on the first arrival travel time information based on a Biharmonic spline interpolation algorithm to obtain a first arrival travel time function; and output a first arrival travel time result according to the first arrival travel time function.

[0040] In a possible design scheme, the first arrival travel time function satisfies: wherein is the Biharmonic operator, and δ is the Dirac function.

[0041] In addition, the technical effects of the seismic first arrival automatic picking device according to the second aspect can refer to the technical effects of the seismic first arrival automatic picking method according to the first aspect, which will not be described here again.

[0042] In a third aspect, the present application provides an electronic device, which comprises a memory for storing executable instructions, and a processor for executing the executable instructions stored in the memory to implement the seismic first arrival automatic picking method according to any one of the first aspect.

[0043] In a fourth aspect, the present application provides a computer readable storage medium, which stores a computer program. When the computer program is run by a processing device, the steps of the seismic first arrival automatic picking method according to any one of the first aspect are executed. BRIEF DESCRIPTION OF DRAWINGS

[0044] Figure 1 A flowchart of a seismic first arrival automatic picking method according to an embodiment of the present application;

[0045] Figure 2 The first picking result obtained through step S201 and the data intercepted through step S203 according to the present application;

[0046] Figure 3 The picking result finally obtained through step S307 according to the present application;

[0047] Figure 4 A flowchart of a seismic first arrival automatic picking method according to an embodiment of the present application;

[0048] Figure 5 A functional module schematic diagram of a seismic first arrival automatic picking device according to an embodiment of the present application.

[0049] Figure 6 A structural schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0050] The technical solutions in the present application will be described below with reference to the accompanying drawings.

[0051] The present application will present various aspects, embodiments or features around a system which can comprise a plurality of devices, components, modules, etc. It should be understood and appreciated that each system can comprise additional devices, components, modules, etc., and / or can not comprise all the devices, components, modules, etc. discussed in connection with the accompanying drawings. In addition, combinations of these solutions can also be used.

[0052] In addition, in the embodiments of the present application, the words "example", "for example", and the like are used herein to mean serving as an instance, illustration, or demonstration. Any embodiment or design presented as an "example" in the present application should not be interpreted as being more preferred or advantageous than other embodiments or design schemes. In fact, the word "example" is used to present a concept in a specific manner. In addition, in the embodiments of the present application, the meaning expressed by "and / or" can be both, or can be one of the two.

[0053] Example One

[0054] The embodiments of the present application will be described below in conjunction with Figure 1 The method for automatically picking up the first arrival of an earthquake provided by the embodiments of the present application will be described in detail. Exemplarily, Figure 1 The flowchart of the method for automatically picking up the first arrival of an earthquake provided by the embodiments of the present application is shown in Figure One .

[0055] As Figure 1 shown, the method for automatically picking up the first arrival of an earthquake includes the following steps:

[0056] In step S201, the first arrival travel time in the first arrival signal in the seismic data to be processed is obtained based on the double sliding window method.

[0057] As an implementation manner, in step S201, the first arrival takeoff time in the first arrival signal in the seismic data to be processed is detected by the double sliding window method, and the first arrival travel time is obtained, wherein the first arrival takeoff time is the time corresponding to when the energy ratio starts to decrease when the first arrival signal starts to enter the next time window.

[0058] It should be understood that, when the first arrival takeoff time in the first arrival signal is detected by the double sliding window method, when a seismic data is processed, the two adjacent windows are relatively stationary, and when the two time windows are both in the noise section or both in the noise signal section, the signal energy ratio calculated by the two time windows is relatively flat. When the first arrival signal in a seismic data starts to enter the previous time window, the next time window is still in the noise stage, and until all the first arrival data enter the previous time window, the energy ratio can reach a peak value. When the first arrival signal starts to enter the next time window, the energy ratio starts to decrease, and the first arrival takeoff time (i.e., the first arrival travel time) can be extracted.

[0059] For example, as Figure 2 shown, the diagram is an experimental data diagram, the horizontal axis is a seismic trace (unit: meter), and the vertical axis is a sampling interval (unit: millisecond). As can be seen from the experimental data diagram, the first arrival takeoff time, i.e., the coordinate position of the point where the first arrival takes off, is: horizontal axis coordinate: 540 meters, and vertical axis coordinate: about 520.

[0060] In a possible implementation, before step S201, the seismic first arrival automatic picking method further includes: pixelizing the original seismic data to obtain the seismic data to be processed.

[0061] The original seismic data refers to a seismic data picture that has not been pixelized. For example, the original seismic data can be seismic data obtained by implementing acquisition, or historical seismic data.

[0062] That is, pixelization refers to a process of converting a seismic data picture into a pixelized effect.

[0063] In the implementation process, the original seismic data is pixelized before being processed by the double sliding window method, so that the seismic data picture is converted into a pixelized effect, thereby effectively improving the data processing efficiency and accuracy of the double sliding window method.

[0064] Step S202: filtering the first arrival travel time based on a filter operator to obtain a filtered travel time.

[0065] As an implementation, step S202 includes: processing the seismic waveform picked up from the seismic data to be processed based on a one-dimensional Gaussian filter function to obtain the filter operator; and filtering the first arrival travel time based on the filter operator to obtain the filtered travel time.

[0066] The filter operator satisfies:

[0067]

[0068] T(itrace) represents the seismic waveform picked up from the seismic data to be processed at the ithtrace, is a one-dimensional Gaussian smoothing filter function.

[0069] It can be understood that the filtered travel time refers to the first arrival travel time after filtering.

[0070] In the implementation process, the seismic waveform picked up from the seismic data to be processed is processed based on the one-dimensional Gaussian filter function to obtain the filter operator, so that the first arrival travel time is filtered based on the filter operator, a more accurate filtered travel time is obtained, and the accuracy of picking up the seismic first arrival travel time is improved.

[0071] Step S203: performing time window cutting on the seismic data to be processed based on the filtered travel time to obtain first arrival travel time information.

[0072] For example, as shown in FIG. 2, a small rectangular frame as shown in FIG. 2 is used for time window cutting. Figure 2 Figure 2 For example, as shown in FIG. 2, a small rectangular frame as shown in FIG. 2 is used for time window cutting.​

[0073] It is understandable that the time window (i.e., time window) adopted in this application can be set according to specific needs and is not specifically limited here. For example, the size of the time window can be in milliseconds or seconds, minutes, etc.

[0074] Step S204: performing first arrival function fitting on the first arrival wave travel time information to obtain a first arrival travel time result.

[0075] As an implementation method, step S204 includes: performing first arrival function fitting on the first arrival wave travel time information based on a Biharmonic spline interpolation algorithm to obtain a first arrival travel time function; and outputting a first arrival travel time result according to the first arrival travel time function.

[0076] Specifically, for a shot's seismic data, its first-arrival travel time can be expressed as a function. Through multi-point Green's function interpolation, the first-arrival travel time function can be accurately reconstructed. Specifically, for a known interpolation fitting problem, the picked first-arrival travel times can be considered as the N control points of the first-arrival function P. Thus, the interpolation problem can be converted into a system of equations (i.e., the first-arrival travel time function) using biharmonic splines:

[0077]

[0078] in, is the Biharmonic operator, and δ is the Dirac function.

[0079] Therefore, the general solution to the above equation is: is the m-dimensional Green's function, where m is 1.

[0080] Among them, the coefficient a j satisfy: Where w(p) represents the value of the initial arrival function P.

[0081] In the above implementation process, by using the Biharmonic spline interpolation algorithm to linearly smooth the picked first arrival wave travel time information to intercept the window, and performing polynomial interpolation on the first arrival wave travel time information within the window, the first arrival travel time function can be accurately reconstructed, thereby obtaining a more accurate first arrival travel time, so as to improve the accuracy of picking the first arrival of the earthquake.

[0082] Step S205 , intercepting the seismic data to be processed according to the first arrival travel time result and the time window, and outputting the target first arrival travel time.

[0083] For example, if Figure 3 , which is a schematic diagram of the final first arrival travel time (ie, the target first arrival travel time) obtained by executing step S205.

[0084] Example Two

[0085] Exemplarily, the embodiment of the present application also provides another method for automatically picking up the first arrival of earthquake, as shown in the flow chart, the method specifically comprises the following steps: Figure 4

[0086] Step S301, obtaining the first arrival time of the first arrival signal in the seismic data to be processed based on the double sliding window method.

[0087] Step S302, filtering the first arrival time based on a filter operator to obtain a filtered arrival time.

[0088] Step S303, performing time window cutting on the seismic data to be processed according to the filtered arrival time to obtain first arrival time information.

[0089] Step S304, performing first arrival function fitting on the first arrival time information to obtain a first arrival time result.

[0090] It should be noted that the specific implementation process of steps S301 to S304 is described above, and will not be repeated here.

[0091] Step S305, cutting the seismic data to be processed according to the first arrival time result and the time window to obtain an intermediate first arrival time.

[0092] Step S306, determining whether the intermediate first arrival time meets the preset requirement.

[0093] The preset requirement can be a pre-set rule, for example, when the intermediate first arrival time is located in a pre-set threshold interval, it is determined that the intermediate first arrival time meets the preset requirement, otherwise, it does not meet. For example, when the intermediate first arrival time is greater than or less than a certain threshold, it is determined that the intermediate first arrival time meets the preset requirement, otherwise, it does not meet.

[0094] It should be understood that the above is only an example, not a limitation.

[0095] Step S307, if yes, taking the intermediate first arrival time as a target first arrival time, and outputting the target first arrival time.

[0096] Step S308, if no, repeating step S302.

[0097] That is, when the intermediate first arrival time does not meet the preset requirement, step S302 is repeatedly executed to cyclically execute the above steps until the target first arrival time is output.

[0098] ​As another implementation, the method for automatically picking up the seismic first arrival provided by the present application can directly repeat the step S302 after the step S305 is performed until the number of times of repeating the step S302 reaches a certain threshold, and then output a result, which is the target first arrival travel time.

[0099] In summary, in combination with Figure 1 and Figure 4 the method for automatically picking up the seismic first arrival shown in FIG. 1, since the double sliding window method is used to obtain the first arrival travel time of the seismic first arrival signal in the seismic data to be processed, the two adjacent windows are relatively static in the process of sliding relative to the received signal. In this method, the first arrival is basically reliable in the initial picking process. Then, based on the picked first arrival travel time information, the linear smoothing intercept window is obtained, and the first arrival function fitting is performed on the first arrival travel time in the window. Therefore, the relatively accurate first arrival travel time can be obtained to realize the automatic picking, so that the manual intervention is required, the workload of the staff is reduced, and the timeliness is improved. In addition, the seismic first arrival picking effect provided by the present application has the characteristics of high precision and strong noise resistance, and good results can be obtained even in the area with low signal-to-noise ratio.

[0100] Example Three

[0101] Based on Figure 1 and Figure 4 The method for automatically picking up the seismic first arrival provided by the embodiment of the present application is described in detail below in combination with Figure 5 The seismic first arrival automatic picking device for performing the method for automatically picking up the seismic first arrival provided by the embodiment of the present application is described in detail.

[0102] Exemplarily, Figure 5 is a structural schematic diagram of the seismic first arrival automatic picking device provided by the embodiment of the present application. Figure One As Figure 5 shown, the seismic first arrival automatic picking device 500 includes a first first arrival travel time calculation module 510, a filtered travel time calculation module 520, a second first arrival travel time calculation module 530, a third first arrival travel time calculation module 540, and an output module 550. For the convenience of description, Figure 5 only the main components of the seismic first arrival automatic picking device are shown. The seismic first arrival automatic picking device 500 performs the functions of the method for automatically picking up the seismic first arrival shown in FIG. 1 or 2. Figure 1

[0103] The first first arrival travel time calculation module 510 is configured to obtain the first arrival travel time in the seismic first arrival signal in the seismic data to be processed based on the double sliding window method.

[0104] ​Optionally, the first first-arrival travel time calculation module 510 is specifically configured to: detect a first-arrival take-off time in a seismic first-arrival signal in the seismic data to be processed by using a double sliding window method to obtain a first-arrival travel time, wherein the first-arrival take-off time is a time corresponding to when an energy ratio starts to decrease when the seismic first-arrival signal starts to enter a later time window.

[0105] It should be understood that, by detecting the first-arrival take-off time in the seismic first-arrival signal by using the double sliding window method, when a seismic data is processed, two adjacent windows are relatively static, and the signal energy ratios calculated by the two time windows are relatively flat when both the two time windows are in a noise section or in a noisy signal section. When the first-arrival signal in a seismic data starts to enter an earlier time window, the later time window is still in a noise stage, and the energy ratio can reach a peak value until all the first-arrival data enter the earlier time window. When the first-arrival signal starts to enter a later time window, the energy ratio starts to decrease, and the first-arrival take-off time (i.e., the first-arrival travel time) can be extracted.

[0106] In a possible embodiment, the seismic first-arrival automatic picking device 500 further comprises a data processing module configured to pixelize original seismic data to obtain seismic data to be processed before obtaining a first-arrival travel time in a seismic first-arrival signal in the seismic data to be processed based on the double sliding window method.

[0107] The filtered travel time calculation module 520 is configured to filter the first-arrival travel time based on a filter operator to obtain a filtered travel time.

[0108] Optionally, the filtered travel time calculation module 520 is specifically configured to: process a seismic waveform picked up from the seismic data to be processed based on a one-dimensional Gaussian filter function to obtain the filter operator; and filter the first-arrival travel time based on the filter operator to obtain the filtered travel time.

[0109] Optionally, the filter operator satisfies:

[0110]

[0111] wherein T(itrace) represents a seismic waveform picked up from the seismic data to be processed in an ithtrace channel, represents a one-dimensional Gaussian smoothing filter function.

[0112] The second first-arrival travel time calculation module 530 is configured to perform time window interception on the seismic data to be processed based on the filtered travel time to obtain first-arrival travel time information.

[0113] The third first-arrival travel time calculation module 540 is configured to perform first-arrival function fitting on the first-arrival travel time information to obtain a first-arrival travel time result.

[0114] Optionally, the third first arrival travel time calculation module 540 is specifically configured to perform first arrival function fitting on the first arrival wave travel time information based on a Biharmonic spline interpolation algorithm to obtain a first arrival travel time function; and output a first arrival travel time result according to the first arrival travel time function.

[0115] Optionally, the first arrival travel time function satisfies wherein, is a Biharmonic operator, and δ is a Dirac function.

[0116] The output module 550 is configured to perform interception on the to-be-processed seismic data according to the first arrival travel time result and the time window, and output a target first arrival travel time.

[0117] In a possible embodiment, the output module 550 is further configured to perform interception on the to-be-processed seismic data according to the first arrival travel time result and the time window to obtain a plurality of intermediate first arrival travel times; determine whether the intermediate first arrival travel times meet a preset requirement; if not, repeatedly run the filter travel time calculation module 520; and if yes, take the intermediate first arrival travel times as the target first arrival travel time, and output the target first arrival travel time.

[0118] It should be noted that the seismic first arrival automatic picking device 500 can be a chip (system) or other components or assemblies provided in an apparatus, and the present application does not limit this.

[0119] In addition, the technical effects of the seismic first arrival automatic picking device 500 can refer to the technical effects of the seismic first arrival automatic picking method provided in any one of Figures 1-2

[0120] Example Four

[0121] Exemplarily, Figure 6 is a structural schematic diagram of an electronic device provided by an embodiment of the present application. The electronic device can be a terminal device or a server device, and can also be a chip (system) or other components or assemblies provided in a terminal device. As Figure 6 shown, the electronic device 2000 can include a processor 2001.

[0122] Optionally, the electronic device 2000 can further include a memory 2002 and / or a transceiver 2003. The processor 2001 is coupled with the memory 2002 and the transceiver 2003, and can be connected through a communication bus.

[0123] The various constituent components of the electronic device 2000 will be specifically introduced below: Figure 6

[0124] ​​The processor 2001 is the core in the electronic device 2000, which can be one processor or a combination of a plurality of processing elements. For example, the processor 2001 can be one or more central processing units (CPUs), application specific integrated circuits (ASICs), or one or more integrated circuits configured to perform the operations of the embodiments of the present application, such as one or more digital signal processors (DSPs), or one or more field programmable gate arrays (FPGAs).

[0125] Optionally, the processor 2001 can execute various functions of the electronic device 2000 by running or executing software programs stored in the memory 2002 and calling data stored in the memory 2002.

[0126] In a specific implementation, as an example, the processor 2001 can include one or more CPUs, such as the CPU0 and the CPU1 shown in FIG. 10. Figure 6

[0127] In a specific implementation, as an example, the electronic device 2000 can also include a plurality of processors, such as the processor 2001 and the processor 2004 shown in FIG. 10. Each of the processors can be a single-CPU or a multi-CPU. The processor herein can refer to one or more devices, circuits, and / or processing cores for processing data (e.g., computer program instructions). Figure 6

[0128] The memory 2002 is configured to store software programs for implementing the solutions of the present application, and the processor 2001 is configured to control the execution of the software programs. The specific implementation can refer to the above method embodiments, and will not be described here.

[0129] ​​Optionally, the memory 2002 can be a read-only memory (ROM) or other type of static storage device that can store static information and instructions, a random access memory (RAM) or other type of dynamic storage device that can store information and instructions, an electrically erasable programmable read-only memory (EEPROM), a compact disc read-only memory (CD-ROM) or other optical disk storage, a magnetic disk storage or other magnetic storage devices, or any other medium capable of storing desired program code in the form of instructions or data structures and that can be accessed by a computer, but is not limited to this. The memory 2002 can be integrated with the processor 2001 or exist independently and be coupled to the processor 2001 through an interface circuit (not shown in the figure) of the electronic device 2000, and embodiments of the present application do not make specific limitations hereon. Figure 6

[0130] The transceiver 2003 is configured to communicate with other terminal devices. For example, the electronic device 2000 is a terminal device, and the transceiver 2003 can be configured to communicate with a network device or another terminal device.

[0131] Optionally, the transceiver 2003 can include a receiver and a transmitter (not shown separately in the figure). The receiver is configured to implement the receiving function, and the transmitter is configured to implement the transmitting function. Figure 6

[0132] Optionally, the transceiver 2003 can be integrated with the processor 2001 or exist independently and be coupled to the processor 2001 through an interface circuit (not shown in the figure) of the electronic device 2000, and embodiments of the present application do not make specific limitations hereon. Figure 6

[0133] It should be noted that the structure of the electronic device 2000 shown in the figure does not constitute a limitation on the electronic device, and an actual electronic device can include more or fewer components than those shown in the figure, or combine certain components, or different component arrangements. Example Five

[0134] In addition, the technical effects of the electronic device 2000 can refer to the technical effects of the automatic first-arrival picking method described in the above method embodiments, which will not be described here again.

[0135] ​ ​​​​

[0136] The embodiment also provides a computer readable storage medium, and the computer readable storage medium stores a computer program. The computer program is run by a processing device to perform the steps of any one of the automatic first arrival picking methods provided by the above embodiments.

[0137] It should be understood that the processor in the embodiments of the present application can be a central processing unit (CPU), and can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor.

[0138] It should also be understood that the memory in the embodiments of the present application can be a volatile memory or a non-volatile memory, or can include both volatile and non-volatile memory. The non-volatile memory can be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically EPROM (EEPROM), or a flash memory. The volatile memory can be a random access memory (RAM) used as an external cache. By way of example, and not limitation, many forms of random access memory (RAM) are available, such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchlink DRAM (SLDRAM), and direct rambus RAM (DR RAM).

[0139] The above-described embodiments can be implemented in whole or in part by software, hardware (e.g., circuitry), firmware, or any combination thereof. When implemented in software, the above-described embodiments can be implemented in the form of a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded and executed on a computer, the processes or functions described in the embodiments of the present application are wholly or partially generated. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium or transferred from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions can be transferred from one website, computer, server, or data center to another website, computer, server, or data center through wired (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium accessible by a computer or a data storage device such as a server, data center, etc. containing one or more available medium collections. The available medium can be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium. The semiconductor medium can be a solid-state disk.

[0140] It should be understood that the term "and / or" herein merely describes an association relationship of associated objects, which means that there can be three relationships, for example, A and / or B can represent the following three cases: A exists alone, A and B exist together, and B exists alone, where A and B can be singular or plural. In addition, the character " / " herein generally represents an "or" relationship between the associated objects before and after it, but it can also represent an "and / or" relationship, which can be understood according to the context before and after it.

[0141] In this application, "at least one" means one or more, and "multiple" means two or more. "At least one of the following" or similar expressions means any combination of these items, including any combination of single or multiple items. For example, at least one of a, b, or c can represent a, b, c, a-b, a-c, b-c, or a-b-c, where a, b, and c can be single or multiple.

[0142] It should be understood that in various embodiments of the present application, the size of the sequence number of the above-described processes does not mean the order of execution, and the execution order of the processes should be determined according to their functions and inherent logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0143] Those skilled in the art can clearly understand that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be realized by electronic hardware or a combination of computer software and electronic hardware. Whether the functions are realized in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to realize the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0144] Those skilled in the art can clearly understand that, for the convenience and brevity of the description, the specific working processes of the above-described system, device and unit can refer to the corresponding processes in the foregoing method embodiments, which will not be repeated here.

[0145] In several embodiments provided in the present application, it should be understood that the disclosed system, device and method can be implemented in other ways. For example, the above-described device embodiments are only schematic, for example, the division of the units is only a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some interface, device or unit, and can be electrical, mechanical or other forms.

[0146] The units described as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, that is, they can be located in one place, or can be distributed on a plurality of network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the embodiment.

[0147] In addition, each functional unit in each embodiment of the present application can be integrated into a processing unit, or each unit can exist physically independently, or two or more units can be integrated into one unit.

[0148] If the functions are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application essentially or the parts that contribute to the prior art or parts of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various media that can store program codes.

[0149] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art can easily think of changes or replacements within the technical scope disclosed in the present application, which should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for automatically picking up earthquake first arrivals, characterized in that: The method comprises: The first arrival travel time of the earthquake first arrival signal in the earthquake data to be processed is obtained based on the double sliding window method; filtering the first arrival travel time based on a filtering operator to obtain a filtered travel time; Performing a second time window interception on the seismic data to be processed according to the filtered travel time to obtain first arrival wave travel time information; Performing first arrival function fitting on the first arrival wave travel time information to obtain a first arrival travel time result; intercepting the seismic data to be processed according to the first arrival travel time result and the second time window, and outputting a target first arrival travel time; The step of intercepting the seismic data to be processed according to the first arrival travel time result and the second time window and outputting the target first arrival travel time comprises: intercepting the seismic data to be processed according to the first arrival travel time result and the second time window to obtain an intermediate first arrival travel time, judging whether the intermediate first arrival travel time meets a preset requirement, and if so, taking the intermediate first arrival travel time as the target first arrival travel time and outputting the target first arrival travel time; if not, repeating the filtering, intercepting, and function fitting steps; The method of obtaining the first arrival travel time of the earthquake first arrival signal in the seismic data to be processed based on the double sliding window method includes: The first arrival take-off time of the earthquake first arrival signal in the earthquake data to be processed is detected by the double sliding window method to obtain the first arrival travel time, wherein the first arrival take-off time is: when the earthquake first arrival signal begins to enter the next time window, the time corresponding to the time when the signal energy ratio calculated by the two time windows begins to decrease, and the two adjacent time windows are relatively static.

2. The method according to claim 1, characterized in that Before obtaining the first arrival travel time of the earthquake first arrival signal in the seismic data to be processed based on the double sliding window method, the method further includes: The original seismic data is pixelated to obtain the seismic data to be processed.

3. The method according to claim 1, characterized in that The filtering the first arrival travel time based on the filtering operator to obtain the filtered travel time includes: Processing the seismic waveform obtained by picking up the seismic data to be processed based on a one-dimensional Gaussian filter function to obtain the filtering operator; The first arrival travel time is filtered according to the filtering operator to obtain a filtered travel time.

4. The method according to claim 3, characterized in that The filtering operator satisfies: ； in, It represents the seismic waveform obtained by picking up the seismic data to be processed in the itrace channel. Represents a one-dimensional Gaussian smoothing filter function.

5. The method according to claim 1, wherein The performing first arrival function fitting on the first arrival wave travel time information to obtain the first arrival travel time result includes: Performing first arrival function fitting on the first arrival wave travel time information based on a Biharmonic spline interpolation algorithm to obtain a first arrival travel time function; Output the first arrival travel time result according to the first arrival travel time function.

6. The method according to claim 5, characterized in that The first arrival travel time function satisfies: ； in, is the Biharmonic operator, 𝛿 is the Dirac function, is the coefficient, the coefficient satisfy: .

7. An automatic earthquake first arrival pickup device, characterized in that: The device comprises: A first first arrival travel time calculation module is used to obtain the first arrival travel time of the earthquake first arrival signal in the earthquake data to be processed based on a double sliding window method; A filtered travel time calculation module, configured to filter the first arrival travel time based on a filtering operator to obtain a filtered travel time; A second first arrival travel time calculation module is used to perform a second time window interception on the seismic data to be processed according to the filtered travel time to obtain first arrival wave travel time information; A third first arrival travel time calculation module is used to perform first arrival function fitting on the first arrival wave travel time information to obtain a first arrival travel time result; an output module, configured to intercept the seismic data to be processed according to the first arrival travel time result and the second time window, and output a target first arrival travel time; The output module is further configured to intercept the seismic data to be processed according to the first arrival travel time result and the second time window to obtain an intermediate first arrival travel time, determine whether the intermediate first arrival travel time meets a preset requirement, and if so, use the intermediate first arrival travel time as a target first arrival travel time and output the target first arrival travel time; if not, repeat the above filtering, intercepting, and function fitting steps; The first arrival travel time calculation module is also used to detect the first arrival takeoff time in the earthquake first arrival signal in the earthquake data to be processed through a double sliding window method to obtain the first arrival travel time, wherein the first arrival takeoff time is: when the earthquake first arrival signal begins to enter the next time window, the time corresponding to the time when the signal energy ratio calculated by the two time windows begins to decrease, and the two adjacent time windows are relatively stationary.

8. An electronic device, characterized in that: include: a memory for storing executable instructions; The processor is configured to implement the earthquake first arrival automatic picking method as described in any one of claims 1 to 6 when executing the executable instructions stored in the memory.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is run by a processing device, the computer program executes the steps of the earthquake first arrival automatic picking method according to any one of claims 1 to 6.

Images