Seismic record string sensing identification method and device, electronic equipment and storage medium

Abstract

The application provides a seismic record string induction identification method and device, electronic equipment and storage medium, and belongs to the technical field of geological exploration. The method comprises the following steps: determining a current calculation time window based on near-surface structure data of a seismic exploration area; determining a signal-to-noise ratio of seismic waves of the seismic exploration area based on seismic wave data in a historical calculation time window of the seismic exploration area; calculating an identification threshold of string induction based on the signal-to-noise ratio of the seismic waves; determining a similarity coefficient of adjacent seismic traces in the current calculation time window based on seismic wave data of the adjacent seismic traces; and identifying whether the adjacent seismic traces in the current calculation time window contain string induction according to the similarity coefficient and the identification threshold. By using the seismic record string induction identification method provided by the application, whether the adjacent seismic traces contain string induction can be quickly identified by calculating the similarity coefficient of the adjacent seismic traces; the problem that string induction is not easy to identify and the identification efficiency is low, thereby affecting the quality of seismic data acquisition, is solved.

Description

Technical Field

[0001] This application relates to the field of geological exploration technology, and more specifically, to a method, apparatus, electronic device, and storage medium for seismic record string induction identification. Background Technology

[0002] Crosstalk is a common type of interference encountered during seismic data acquisition. It arises from electromagnetic induction interference from external sources and has been a long-standing problem affecting the quality of seismic data acquisition for many years. In raw seismic records, crosstalk typically appears simultaneously on some or all of a continuous seismic trace, and its waveform resembles a series of pulses, hence the name "crosstalk." Crosstalk is characterized by its wide bandwidth and the fact that it affects multiple seismic data traces, severely damaging seismic data.

[0003] In existing technologies, the identification of cross-current induction relies on the experience of researchers, involving visual identification of seismic wave data from each shot. This method not only demands a high level of experience from the personnel but is also prone to misjudgment and has low efficiency. Therefore, the difficulty in identifying cross-current induction and its low efficiency, which affects the quality of seismic data acquisition, is a problem that urgently needs to be solved. Summary of the Invention

[0004] This application provides a method, apparatus, electronic device, and storage medium for identifying seismic record string induction, aiming to solve the problems of difficult identification and low identification efficiency of string induction, which affect the quality of seismic data acquisition.

[0005] The first aspect of this application provides a method for seismic record string induction identification, including:

[0006] Based on near-surface structure data of the seismic exploration area, determine the current calculation window;

[0007] Based on the seismic wave data within the historical calculation window of the seismic exploration area, the signal-to-noise ratio of the seismic waves in the seismic exploration area is determined;

[0008] Based on the signal-to-noise ratio of seismic waves in the seismic exploration area, the identification threshold of cross-induction in the seismic exploration area is calculated.

[0009] Based on the seismic wave data of adjacent seismic traces within the current calculation window, determine the similarity coefficient of adjacent seismic traces;

[0010] Based on the similarity coefficient of adjacent seismic traces and the identification threshold of cross-induction within the seismic exploration area, it is determined whether there is cross-induction in adjacent seismic traces within the current calculation window.

[0011] Optionally, based on the similarity coefficient of adjacent seismic traces and the identification threshold of cross-induction within the seismic exploration area, it is determined whether adjacent seismic traces within the current calculation window contain cross-induction, including:

[0012] When the similarity coefficient of adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area, it is determined that there is cross-induction in the adjacent seismic traces in the current calculation window.

[0013] When the similarity coefficient of adjacent seismic traces is less than the identification threshold of cross-induction within the seismic exploration area, it is determined that there is no cross-induction in the adjacent seismic traces within the current calculation window.

[0014] Optionally, based on near-surface structure data of the seismic exploration area, the current calculation window is determined, including:

[0015] Based on near-surface structure data, the first time window for the first arrival of seismic waves and the second time window for the appearance of reflected waves in the seismic exploration area were determined.

[0016] Within the time window where seismic waves occur, any time window is determined as the current calculation time window. The time window where seismic waves occur does not include the first time window and the second time window.

[0017] Optionally, the near-surface structure data may include at least: the distance between the excitation point and the receiver point, the thickness and velocity of the weathering layer, and the velocity and thickness of the deceleration layer.

[0018] Optionally, the seismic wave data within the historical calculation window of the seismic exploration area includes the amplitude values ​​of the series induction samples and the amplitude values ​​of the environmental noise samples within the historical calculation window; based on the seismic wave data within the historical calculation window of the seismic exploration area, the signal-to-noise ratio of the seismic waves in the seismic exploration area is determined, including:

[0019] Based on the amplitude values ​​of the series induction samples and the ambient noise samples, the energy values ​​of the series induction samples and the ambient noise samples are determined.

[0020] The signal-to-noise ratio of seismic waves is determined based on the ratio of the energy values ​​of series induction samples to the energy values ​​of environmental noise samples.

[0021] Optionally, the seismic wave data of adjacent seismic traces within the current calculation window includes the amplitude value of each sample point in each adjacent seismic trace; based on the seismic wave data of adjacent seismic traces within the current calculation window, the similarity coefficient of adjacent seismic traces is determined, including:

[0022] The similarity coefficient between adjacent seismic traces is determined based on the amplitude value of each sample point in the adjacent seismic trace.

[0023] A second aspect of this application provides a seismic record string induction identification device, comprising:

[0024] The first determination module is used to determine the current calculation window based on near-surface structure data of the seismic exploration area.

[0025] The second determination module is used to determine the signal-to-noise ratio of seismic waves in the seismic exploration area based on the seismic wave data within the historical calculation time window of the seismic exploration area.

[0026] First calculation module: used to calculate the identification threshold of cross-induction in the seismic exploration area based on the signal-to-noise ratio of seismic waves in the seismic exploration area;

[0027] The third determination module is used to determine the similarity coefficient of adjacent seismic traces based on the seismic wave data of adjacent seismic traces within the current calculation window.

[0028] Identification module: Used to identify whether there is cross-induction in adjacent seismic traces within the current calculation window, based on the similarity coefficient of adjacent seismic traces and the identification threshold of cross-induction in the seismic exploration area.

[0029] Optionally, the recognition module includes:

[0030] First determining unit: used to determine that there is cross-induction in adjacent seismic traces within the current calculation window when the similarity coefficient of adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area.

[0031] The second determining unit is used to determine that there is no cross-induction in adjacent seismic traces within the current calculation window when the similarity coefficient of adjacent seismic traces is less than the identification threshold of cross-induction in the seismic exploration area.

[0032] A third aspect of this application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the seismic record string induction identification method provided in the first aspect of this application.

[0033] The fourth aspect of this application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the seismic record string induction identification method provided in the first aspect of this application.

[0034] Beneficial effects:

[0035] The seismic record crosstalk identification method, device, electronic device, and storage medium provided in this application are based on the high consistency of crosstalk between adjacent seismic traces. By calculating the similarity coefficient between adjacent seismic traces within the current calculation window, and obtaining the intensity of environmental noise and crosstalk within the seismic exploration area based on historical data of crosstalk-inducing seismic traces acquired within historical calculation windows, the signal-to-noise ratio (SNR) of seismic waves in the seismic exploration area is determined. Using the SNR of seismic waves in the seismic exploration area, the identification threshold for crosstalk within the seismic exploration area is calculated. By comparing the similarity coefficient of adjacent seismic traces with the identification threshold for crosstalk within the seismic exploration area, it is possible to quickly identify whether adjacent seismic traces within the current calculation window contain crosstalk, thereby achieving rapid monitoring of seismic acquisition data quality. This solves the problems of difficult and inefficient identification of crosstalk, which affects the quality of seismic data acquisition. Attached Figure Description

[0036] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0037] Figure 1 This is a flowchart of a seismic record string induction identification method proposed in an embodiment of this application;

[0038] Figure 2A This is a schematic diagram of the related post-seismic data traces containing cross-induction within a historical calculation window according to an embodiment of this application;

[0039] Figure 2B This is a schematic diagram of the time spectrum of the related post-seismic data trace containing cross-induction within a historical calculation time window proposed in an embodiment of this application;

[0040] Figure 3A This is a schematic diagram of one set of correlated post-seismic records containing serial induction, as proposed in an embodiment of this application;

[0041] Figure 3B This is a schematic diagram of another set of correlated post-seismic records containing serial induction, as proposed in one embodiment of this application;

[0042] Figure 3C This is a schematic diagram of another set of correlated post-seismic records containing serial induction, as proposed in an embodiment of this application;

[0043] Figure 3D This is a schematic diagram of the similarity coefficient curves of adjacent seismic traces containing series induction, as proposed in an embodiment of this application.

[0044] Figure 4 This is a structural block diagram of a seismic record string induction identification device proposed in an embodiment of this application. Detailed Implementation

[0045] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0046] In related technologies, the identification of cross-current induction relies on the experience of researchers, involving visual identification of seismic wave data from each shot. This method not only demands a high level of experience from the personnel but is also prone to misjudgment and has low identification efficiency. Therefore, the difficulty in identifying cross-current induction and its low identification efficiency, which affects the quality of seismic data acquisition, is a problem that urgently needs to be solved.

[0047] In light of this, we conducted extensive experiments using previously acquired historical seismic data. The experiments revealed that crosstalk manifests as a series of pulses with varying energy. By performing correlation operations on this series of pulses—compressing them into a single pulse signal—the crosstalk in the correlated seismic record appears as continuous vibrations with frequency varying over time. Adjacent seismic traces exhibit high consistency, especially before the first arrival, in terms of energy and waveform, while showing significant differences from normal seismic traces. (Refer to...) Figure 2A This shows a schematic diagram of the relevant post-seismic data traces containing series induction within the historical calculation window of this application, with reference to... Figure 2B This illustrates a time-spectrum diagram of the correlated post-seismic data traces containing cross-induction within the historical calculation window of this application, as shown below. Figure 2A and Figure 2B As shown, there is a high degree of consistency in both energy and waveform between adjacent seismic traces, and the frequency varies over time. Based on the high consistency of cross-induction between adjacent seismic traces, this application can quickly identify cross-induction in seismic records by calculating and analyzing the similarity between adjacent seismic data traces.

[0048] Example 1

[0049] Reference Figure 1 The flowchart of a seismic record string induction identification method according to this application is shown, as follows: Figure 1 As shown, the seismic record string induction identification method of this application includes the following steps:

[0050] Step S1: Determine the current calculation window based on near-surface structure data of the seismic exploration area.

[0051] In this embodiment, the current calculation window is the window for calculating the similarity coefficient between current seismic traces. Based on the near-surface structure data of the seismic exploration area, the time window of the seismic wave that will affect the similarity coefficient calculation is calculated. The selection of the current calculation window avoids the time window of the seismic wave that will affect the similarity coefficient calculation, thereby avoiding the problem of low accuracy of the similarity coefficient calculation.

[0052] Step S2: Determine the signal-to-noise ratio of seismic waves in the seismic exploration area based on the seismic wave data within the historical calculation window of the seismic exploration area.

[0053] The seismic wave data within the historical calculation window of the seismic exploration area are seismic trace data containing cross-induction collected within the historical calculation window. The signal-to-noise ratio of the seismic waves in the seismic exploration area is determined by the intensity of environmental noise and the intensity of cross-induction within the seismic exploration area.

[0054] Step S3: Calculate the identification threshold of cross-induction within the seismic exploration area based on the signal-to-noise ratio of seismic waves in the seismic exploration area.

[0055] The threshold for identifying cross-induction within a seismic exploration area is used to compare with the similarity coefficient of adjacent seismic traces to determine whether there is cross-induction within adjacent seismic traces.

[0056] Step S4: Determine the similarity coefficient of adjacent seismic traces based on the seismic wave data of adjacent seismic traces within the current calculation window.

[0057] Within the current calculation window, seismic wave data for each adjacent seismic trace is acquired. The similarity coefficient between adjacent seismic traces is then calculated using this data. The number of adjacent seismic traces can be two, three, or any number greater than two; no limitation is imposed here.

[0058] Step S5: Based on the similarity coefficient of adjacent seismic traces and the identification threshold of cross-induction in the seismic exploration area, identify whether there is cross-induction in adjacent seismic traces within the current calculation window.

[0059] Reference Figure 3A This illustrates a schematic diagram of one set of correlated post-seismic records containing series induction, as shown in this application. Figure 3B This illustrates another set of correlated post-seismic records containing series induction, as shown in the schematic diagram of this application. Figure 3C This illustrates a schematic diagram of another set of correlated post-seismic records containing series induction, as shown in the present application. Figure 3DThis paper illustrates a schematic diagram of the similarity coefficient curves of adjacent seismic traces containing series induction, as shown in the example. Figure 3A , 3B As shown in Figures 3C and 3D, the calculation time window is 0-200 milliseconds. Cross-induction traces can be quickly identified from the similarity coefficient curves. This seismic record contains three sets of cross-induction traces of different intensities. Figure 3D The boxed sections, from left to right, are as follows: Figure 3A , 3B Similarity coefficients corresponding to three different intensities of cross-induction in 3C. By comparing the similarity coefficients of adjacent seismic traces with the identification threshold of cross-induction within the seismic exploration area, it is possible to identify whether there is cross-induction in adjacent seismic traces within the current calculation window.

[0060] The seismic record cross-induction identification method provided in this application is based on the high consistency of cross-induction between adjacent seismic traces. It calculates the similarity coefficient between adjacent seismic traces within the current calculation window, and obtains the intensity of environmental noise and cross-induction within the seismic exploration area based on historical data of seismic traces containing cross-induction acquired within historical calculation windows. This determines the signal-to-noise ratio (SNR) of the seismic waves in the exploration area. The cross-induction identification threshold within the exploration area is calculated using the SNR. By comparing the similarity coefficient of adjacent seismic traces with the cross-induction identification threshold, it is possible to quickly identify whether adjacent seismic traces within the current calculation window contain cross-induction, thus achieving rapid monitoring of seismic acquisition data quality. This solves the problems of difficult and inefficient identification of cross-induction, which affects the quality of seismic data acquisition.

[0061] Based on the above-mentioned seismic record string induction identification method, this application provides the following specific implementation examples. Under the premise of not contradicting each other, the various examples can be arbitrarily combined to form another seismic record string induction identification method. It should be understood that any seismic record string induction identification method formed by combining any examples should fall within the protection scope of this application.

[0062] In one feasible implementation, the method for identifying whether adjacent seismic traces within the current calculation window contain cross-induction is based on the similarity coefficient of adjacent seismic traces and the identification threshold of cross-induction within the seismic exploration area. This includes the following steps:

[0063] Step S51: When the similarity coefficient of adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area, it is determined that the adjacent seismic traces in the current calculation window contain cross-induction.

[0064] When the similarity coefficient of adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area, cross-induction can be quickly identified in adjacent seismic traces within the current calculation window through numerical comparison. This allows for the rapid output of seismic trace information containing cross-induction, enabling rapid monitoring of the quality of seismic acquisition data.

[0065] Step S52: When the similarity coefficient of adjacent seismic traces is less than the identification threshold of cross-induction in the seismic exploration area, it is determined that there is no cross-induction in the adjacent seismic traces in the current calculation window.

[0066] When the similarity coefficient of adjacent seismic traces is less than the identification threshold of cross-induction in the seismic exploration area, that is, by comparing the values, it can be quickly identified that there is no cross-induction in the adjacent seismic traces in the current calculation window. Then the data acquired by that seismic trace can be used normally, realizing rapid monitoring of the quality of seismic acquisition data.

[0067] In one feasible implementation, the current calculation window is determined based on near-surface structure data of the seismic exploration area, including:

[0068] Step S11: Based on near-surface structure data, determine the first time window for the first arrival of seismic waves and the second time window for the appearance of reflected waves in the seismic exploration area.

[0069] Based on near-surface structure data of the seismic exploration area, the first time window for the first arrival wave and the second time window for the reflected wave, which will affect the similarity coefficient calculation, are calculated.

[0070] Step S12: Within the time window of seismic wave occurrence, determine any time window as the current calculation time window. The time window of seismic wave occurrence does not include the first time window and the second time window.

[0071] The presence of first arrival waves and reflected waves can interfere with the calculation of similarity coefficients of adjacent seismic traces. Therefore, the current calculation time window needs to avoid the time window when first arrival waves and reflected waves appear, thereby improving the accuracy of similarity coefficient calculation.

[0072] In another feasible implementation, the near-surface structure data includes at least: the distance between the excitation point and the receiver point, the thickness and velocity of the weathering layer, and the velocity and thickness of the deceleration layer, which are calculated for the first and second time windows using the following formula:

[0073]

[0074]

[0075]

[0076]

[0077] Wherein, x is the distance between the shot point and the receiving point, v1 is the velocity of the weathered layer, t1 is the time for the direct wave to reach the receiving point, h1 is the thickness of the weathered layer, v2 is the velocity of the low-velocity layer, t2 is the time for the refracted wave to reach the receiving point, when t1<t2, the first-break wave is a direct wave, when t1>t2, the first-break wave is a refracted wave, h is the buried depth of the top surface of the reflection layer, v is the stacking velocity, and t is the time for the reflected wave to reach the receiving point.

[0078] In a feasible implementation, the seismic wave data within the historical calculation time window of the seismic exploration area includes the amplitude values of string induction sample points within the historical calculation time window and the amplitude values of environmental noise sample points within the historical calculation time window; determining the signal-to-noise ratio of seismic waves in the seismic exploration area based on the seismic wave data within the historical calculation time window of the seismic exploration area comprises the following steps:

[0079] Step S21: determining the energy value of the string induction sample points and the energy value of the environmental noise sample points based on the amplitude values of the string induction sample points and the amplitude values of the environmental noise sample points.

[0080] A plurality of sample points are arranged in the seismic exploration area, the amplitude value of string induction and the amplitude value of environmental noise of each sample point within the historical calculation time window of the seismic exploration area are obtained, the square of the amplitude value is the energy value, so as to determine the energy value of string induction and the energy value of environmental noise of each sample point within the historical calculation time window.

[0081] Step S22: determining the signal-to-noise ratio of the seismic wave based on the ratio of the energy value of the string induction sample points to the energy value of the environmental noise sample points.

[0082]

[0083] Wherein, SNR is the signal-to-noise ratio, i is the sequence number of the sample point in the current calculation time window, S i is the amplitude value of the string induction sample point in the current calculation time window, n i is the amplitude value of the environmental noise sample point in the current calculation time window, S i 2 is the energy value of the string induction sample point in the current calculation time window, n i 2 is the energy value of the environmental noise sample point in the current calculation time window.

[0084] In a feasible implementation, the seismic wave data of adjacent seismic traces within the current calculation time window includes the amplitude value of each seismic trace sample point in the adjacent seismic traces; determining the similarity coefficient of adjacent seismic traces based on the seismic wave data of adjacent seismic traces within the current calculation time window comprises the following steps:

[0085] Step S41: determining the similarity coefficient of the adjacent seismic traces based on the amplitude value of each seismic trace sample point in the adjacent seismic traces.

[0086]

[0087] Where r represents the similarity coefficient between adjacent seismic traces, i is the sample point sequence number within the calculation window, and x and y represent two adjacent seismic traces, respectively. i and y i This represents the amplitude value of each sample point in the adjacent seismic trace with index i.

[0088] In one feasible implementation, the identification threshold for cross-induction within the seismic exploration area is calculated based on the signal-to-noise ratio of seismic waves in the seismic exploration area, including the following steps:

[0089] When the calculation window contains only random environmental noise and cross-induction, the seismic trace data within the calculation window can be represented as:

[0090] x(t)=s(t)+n(t) (3)

[0091] Where x(t) is the output data of the seismic trace within the calculation window, s(t) is the series induction, n(t) is the environmental noise, and t is the time.

[0092] Substituting formula (3) into formula (2) yields formula (4):

[0093]

[0094] Where j is the seismic data trace sequence number, i is the sample point sequence number within the calculation window, and S j,i Let n be the amplitude value of the series induction sample point corresponding to sample point i in seismic data channel j. j,i Let j be the amplitude value of the environmental noise sample point corresponding to the seismic data channel j and sample point i.

[0095] Based on the assumption that environmental noise is random, equation (4) is simplified to:

[0096]

[0097] Based on the characteristics of cross-induction in the relevant post-earthquake records, formula (5) can be further simplified to:

[0098]

[0099] Substituting formula (1) into formula (6) yields the recognition threshold R for string sensing:

[0100]

[0101] Example 2

[0102] Based on the same inventive concept, another embodiment of this application provides a seismic record string induction identification device, which is used to perform the seismic record string induction identification method as provided in Embodiment 1 of this application; see reference Figure 4 The diagram shows the structural block diagram of the seismic record string induction identification device, such as... Figure 4 As shown, the seismic record string induction identification device includes:

[0103] First Determination Module 11: Used to determine the current calculation window based on near-surface structure data of the seismic exploration area;

[0104] Second determination module 12: used to determine the signal-to-noise ratio of seismic waves in the seismic exploration area based on the seismic wave data within the historical calculation time window of the seismic exploration area;

[0105] First calculation module 13: used to calculate the identification threshold of cross-induction in the seismic exploration area based on the signal-to-noise ratio of seismic waves in the seismic exploration area;

[0106] The third determination module 14 is used to determine the similarity coefficient of adjacent seismic traces based on the seismic wave data of adjacent seismic traces within the current calculation window.

[0107] Identification module 15: Used to identify whether there is cross-induction in adjacent seismic traces within the current calculation window, based on the similarity coefficient of adjacent seismic traces and the identification threshold of cross-induction in the seismic exploration area.

[0108] The seismic record crosstalk identification device provided in this application is based on the high consistency of crosstalk between adjacent seismic traces. It calculates the similarity coefficient between adjacent seismic traces within the current calculation window, and obtains the intensity of environmental noise and crosstalk within the seismic exploration area based on historical data of seismic traces containing crosstalk within the calculation window. This determines the signal-to-noise ratio (SNR) of the seismic waves in the seismic exploration area. Using the SNR of the seismic waves in the seismic exploration area, it calculates the identification threshold for crosstalk within the seismic exploration area. By comparing the similarity coefficient of adjacent seismic traces with the identification threshold for crosstalk within the seismic exploration area, it can quickly identify whether adjacent seismic traces within the current calculation window contain crosstalk, thereby achieving rapid monitoring of the quality of seismic acquisition data. This solves the problems of difficult and inefficient identification of crosstalk, which affects the quality of seismic data acquisition.

[0109] In one feasible implementation, the identification module 15 includes:

[0110] First determining unit 151: used to determine that there is cross-induction in adjacent seismic traces within the current calculation window when the similarity coefficient of adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area.

[0111] The second determining unit 152 is used to determine that there is no cross-induction in the adjacent seismic traces within the current calculation window when the similarity coefficient of the adjacent seismic traces is less than the identification threshold of cross-induction in the seismic exploration area.

[0112] In one feasible implementation, the first determining module 11 includes:

[0113] The third determining unit 111 is used to determine the first time window for the first arrival of seismic waves and the second time window for the appearance of reflected waves in the seismic exploration area based on near-surface structure data.

[0114] Fourth determining unit 112: used to determine any time window as the current calculation time window within the time window of the occurrence of seismic waves, wherein the time window of the occurrence of seismic waves does not include the first time window and the second time window.

[0115] In one feasible implementation, the near-surface structure data includes at least: the distance between the excitation point and the receiver point, the thickness and velocity of the weathering layer, and the velocity and thickness of the deceleration layer.

[0116] In another feasible implementation, the seismic wave data within the historical calculation window of the seismic exploration area includes the amplitude values ​​of the series induction samples within the historical calculation window and the amplitude values ​​of the environmental noise samples within the historical calculation window. The second determining module 12 includes:

[0117] Fifth determining unit 121: used to determine the energy value of the series induction sample and the energy value of the environmental noise sample based on the amplitude value of the series induction sample and the amplitude value of the environmental noise sample;

[0118] The sixth determining unit 122 is used to determine the signal-to-noise ratio of seismic waves based on the ratio of the energy value of the series induction sample points to the energy value of the environmental noise sample points.

[0119] In one feasible implementation, the seismic wave data of adjacent seismic traces within the current calculation window includes the amplitude value of each sample point in the adjacent seismic traces; the third determination module 14 includes:

[0120] The seventh determining unit 141 is used to determine the similarity coefficient of adjacent seismic traces based on the amplitude value of each sample point in the adjacent seismic trace.

[0121] Example 3

[0122] This application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the computer program is executed by the processor, it implements the seismic record string induction identification method provided in the first aspect of this application.

[0123] The electronic device provided in this application, based on the high consistency of cross-induction between adjacent seismic traces, calculates the similarity coefficient between adjacent seismic traces within the current calculation window. It also obtains the intensity of environmental noise and cross-induction within the seismic exploration area based on historical data of cross-induced seismic traces acquired within historical calculation windows, thereby determining the signal-to-noise ratio (SNR) of seismic waves in the exploration area. Using the SNR of the seismic waves in the exploration area, it calculates the identification threshold for cross-induction within the exploration area. By comparing the similarity coefficient of adjacent seismic traces with the identification threshold for cross-induction within the exploration area, it can quickly identify whether adjacent seismic traces within the current calculation window contain cross-induction, thus achieving rapid monitoring of seismic acquisition data quality. This solves the problem of difficult and inefficient identification of cross-induction, which affects the quality of seismic data acquisition.

[0124] Example 4

[0125] This application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the seismic record string induction identification method provided in the first aspect of this application.

[0126] The computer-readable storage medium provided in this application, based on the high consistency of cross-induction between adjacent seismic traces, calculates the similarity coefficient between adjacent seismic traces within the current calculation window. It then obtains the intensity of environmental noise and cross-induction within the seismic exploration area based on historical data of cross-induced seismic traces acquired within historical calculation windows, determining the signal-to-noise ratio (SNR) of seismic waves in the exploration area. Using the SNR of the seismic waves in the exploration area, it calculates the identification threshold for cross-induction within the exploration area. By comparing the similarity coefficient of adjacent seismic traces with the identification threshold for cross-induction within the exploration area, it can quickly identify whether adjacent seismic traces within the current calculation window contain cross-induction, thus achieving rapid monitoring of seismic acquisition data quality. This solves the problems of difficult and inefficient identification of cross-induction, which affects the quality of seismic data acquisition.

[0127] As the device embodiment is basically similar to the method embodiment, the description is relatively simple, and relevant parts can be found in the description of the method embodiment.

[0128] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0129] It should be understood that although preferred embodiments of the present application have been described in this specification, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of this application.

[0130] The seismic record string induction identification method, device, electronic device, and storage medium provided in this application have been described in detail above. Specific examples have been used to illustrate the principles and implementation methods of this application. The description of the above embodiments is only for the purpose of helping to understand the method and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for inductive identification of seismic record strings, characterized in that, The method includes: Based on near-surface structure data of the seismic exploration area, determine the current calculation window; Based on the seismic wave data within the historical calculation time window of the seismic exploration area, the signal-to-noise ratio of the seismic waves in the seismic exploration area is determined; Based on the signal-to-noise ratio of seismic waves in the seismic exploration area, calculate the identification threshold of cross-induction in the seismic exploration area. Based on the seismic wave data of adjacent seismic traces within the current calculation window, the similarity coefficient of the adjacent seismic traces is determined; Based on the similarity coefficient of the adjacent seismic traces and the identification threshold of cross-induction within the seismic exploration area, it is determined whether the adjacent seismic traces within the current calculation window contain cross-induction. Among them, based on near-surface structure data of the seismic exploration area, the current calculation window is determined, including: Based on the near-surface structure data, the first time window for the first arrival of seismic waves and the second time window for the occurrence of reflected waves in the seismic exploration area are determined. Within the time window in which the seismic wave occurs, any time window is determined as the current calculation time window, wherein the time window in which the seismic wave occurs does not include the first time window and the second time window.

2. The method according to claim 1, characterized in that, Based on the similarity coefficient of the adjacent seismic traces and the identification threshold of cross-induction within the seismic exploration area, the method identifies whether the adjacent seismic traces within the current calculation window contain cross-induction, including: When the similarity coefficient of the adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area, it is determined that the adjacent seismic traces in the current calculation window contain cross-induction. When the similarity coefficient of the adjacent seismic traces is less than the identification threshold of cross-induction in the seismic exploration area, it is determined that the adjacent seismic traces in the current calculation window do not contain cross-induction.

3. The method according to claim 1, characterized in that, The near-surface structure data includes at least: the distance between the excitation point and the receiver point, the thickness and velocity of the weathering layer, and the velocity and thickness of the deceleration layer.

4. The method according to claim 1, characterized in that, The seismic wave data within the historical calculation window of the seismic exploration area includes the amplitude values ​​of the series induction samples and the amplitude values ​​of the environmental noise samples within the historical calculation window; Based on the seismic wave data within the historical calculation window of the seismic exploration area, the signal-to-noise ratio of the seismic waves in the seismic exploration area is determined, including: Based on the amplitude values ​​of the series induction samples and the amplitude values ​​of the environmental noise samples, the energy values ​​of the series induction samples and the environmental noise samples are determined. The signal-to-noise ratio of the seismic wave is determined based on the ratio of the energy values ​​of the series induction samples to the energy values ​​of the environmental noise samples.

5. The method according to claim 1, characterized in that, The seismic wave data of adjacent seismic traces within the current calculation window includes the amplitude value of each sample point in the adjacent seismic traces; Based on the seismic wave data of adjacent seismic traces within the current calculation window, the similarity coefficient of the adjacent seismic traces is determined, including: The similarity coefficient of the adjacent seismic traces is determined based on the amplitude value of each sample point in the adjacent seismic traces.

6. A seismic record string induction identification device, characterized in that, The device includes: The first determination module is used to determine the current calculation window based on near-surface structure data of the seismic exploration area. The second determining module is used to determine the signal-to-noise ratio of seismic waves in the seismic exploration area based on the seismic wave data within the historical calculation time window of the seismic exploration area. First calculation module: used to calculate the identification threshold of cross-induction in the seismic exploration area based on the signal-to-noise ratio of seismic waves in the seismic exploration area; The third determining module is used to determine the similarity coefficient of adjacent seismic traces based on the seismic wave data of adjacent seismic traces within the current calculation window. Identification module: used to identify whether there is cross-induction in the adjacent seismic traces within the current calculation window, based on the similarity coefficient of the adjacent seismic traces and the identification threshold of cross-induction in the seismic exploration area; The first determining module includes: The third determining unit is used to determine the first time window for the first arrival of seismic waves and the second time window for the appearance of reflected waves in the seismic exploration area based on near-surface structure data. The fourth determining unit is used to determine any time window as the current calculation time window within the time window in which the seismic wave occurs. The time window in which the seismic wave occurs does not include the first time window and the second time window.

7. The apparatus according to claim 6, characterized in that, The recognition module includes: First determining unit: used to determine that the adjacent seismic traces in the current calculation window contain cross-induction when the similarity coefficient of the adjacent seismic traces is greater than or equal to the identification threshold of cross-induction in the seismic exploration area; The second determining unit is used to determine that there is no cross-induction in the adjacent seismic traces within the current calculation window when the similarity coefficient of the adjacent seismic traces is less than the identification threshold of cross-induction in the seismic exploration area.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the seismic record string induction identification method as described in any one of claims 1 to 5.

9. A computer-readable storage medium storing a computer program that, when executed by a processor, implements the seismic record string induction identification method as described in any one of claims 1 to 5.

Images