Method and device for determining quality parameters of a controllable source and computer equipment

Abstract

The application provides a method and device for determining a quality parameter of a controllable source and a computer device, and belongs to the technical field of geophysical exploration. The method comprises the following steps: determining a first force signal and a second force signal, the first force signal being a time-domain force signal excited by the controllable source, and the second force signal being a frequency-domain force signal corresponding to the first force signal; segmenting the second force signal to obtain a plurality of sub-signals, the scanning bandwidth of each sub-signal being less than one time of the starting scanning frequency of the sub-signal; determining a target frequency band in the second force signal based on the effective scanning point number and the plurality of sub-signals, the target frequency band being a frequency band corresponding to a fundamental wave signal in the first force signal; performing time-domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal; and obtaining the quality parameter of the controllable source based on the first force signal and the fundamental wave signal, the quality parameter being used for determining the quality of the first force signal, and the accuracy of determining the quality of the force signal of the controllable source is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geophysical exploration, and particularly relates to a method and device for determining a quality parameter of a controllable source and a computer device. BACKGROUND

[0002] Using a controllable source to excite seismic waves is an important method for seismic exploration. When the controllable source excites seismic waves, the vibrator system and the hydraulic system of the controllable source apply seismic waves according to a preset scanning signal, and the output force signal is the seismic signal to be obtained. Because the coupling relationship between the vibrator system, the hydraulic system, etc. of the controllable source and the stratum is irregular, the force signal not only includes the fundamental wave signal in the scanning signal, but also includes a harmonic signal that interferes with the fundamental wave signal. Therefore, some parameters are needed to determine the interference degree of the harmonic signal on the fundamental wave signal to determine the quality of the force signal.

[0003] In related technologies, a distortion parameter is generally used to determine the quality of the force signal. The distortion parameter is used to represent the difference between the force signal and the scanning signal, so that the distortion parameter can only qualitatively determine the interference degree of the harmonic signal on the fundamental wave signal, thereby reducing the accuracy of determining the quality of the force signal. SUMMARY

[0004] Embodiments of the present application provide a method and device for determining a quality parameter of a controllable source and a computer device, which can improve the accuracy of determining the quality of the force signal of the controllable source. The technical solution is as follows:

[0005] In one aspect, a method for determining a quality parameter of a controllable source is provided, and the method comprises:

[0006] determining a first force signal and a second force signal, the first force signal being a time-domain force signal excited by the controllable source, and the second force signal being a frequency-domain force signal corresponding to the first force signal;

[0007] segmenting the second force signal to obtain a plurality of sub-signals, and the scanning bandwidth of each sub-signal being less than one time the starting scanning frequency of the sub-signal;

[0008] determining the number of effective scanning frequency points of the second force signal;

[0009] based on the number of effective scanning frequency points and the plurality of sub-signals, determining a target frequency band in the second force signal, the target frequency band being a frequency band corresponding to the fundamental wave signal in the first force signal;

[0010] performing time-domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal;

[0011] obtaining a quality parameter of the vibroseis based on the first force signal and the fundamental signal, the quality parameter being used to determine the quality of the first force signal.

[0012] In a possible implementation, the segmenting the second force signal to obtain the plurality of sub-signals comprises:

[0013] obtaining a scanning signal corresponding to the first force signal;

[0014] In a case where the scanning signal is a linear signal, segmenting the second force signal based on a frequency sampling interval of a complete force signal to obtain the plurality of sub-signals, the complete force signal being a force signal of a full frequency range corresponding to the first force signal;

[0015] In a case where the scanning signal is a nonlinear signal, determining a plurality of frequency multiplication time points, and segmenting the second force signal based on the plurality of frequency multiplication time points to obtain the plurality of sub-signals.

[0016] In a possible implementation, the segmenting the second force signal based on the frequency sampling interval of the complete force signal to obtain the plurality of sub-signals comprises:

[0017] obtaining a starting scanning frequency of the first force signal;

[0018] In a case where the sub-signal is a first sub-signal, determining the starting scanning frequency of the first force signal as a starting scanning frequency of the sub-signal, and determining a difference between twice the starting scanning frequency of the first force signal and the frequency sampling interval as an ending scanning frequency of the sub-signal;

[0019] In a case where the sub-signal is not the first sub-signal, determining a sum of an ending scanning frequency of a previous sub-signal and the frequency sampling interval as a starting scanning frequency of the sub-signal, and determining a difference between twice the starting scanning frequency of the sub-signal and the frequency sampling interval as an ending scanning frequency of the sub-signal.

[0020] In a possible implementation, the segmenting the second force signal based on the plurality of frequency multiplication time points to obtain the plurality of sub-signals comprises:

[0021] obtaining a starting scanning frequency of the first force signal;

[0022] In a case where the sub-signal is a first sub-signal, determining the starting scanning frequency of the first force signal as a starting scanning frequency of the sub-signal, and taking a scanning frequency corresponding to a first frequency multiplication time point as an ending scanning frequency of the sub-signal;

[0023] In the case that the sub-signal is not the first sub-signal, a time sampling interval of a complete force signal is obtained, a sum of the time sampling interval and a previous frequency multiplication time point is determined as a starting time point of the sub-signal, a scanning frequency corresponding to the starting time point is determined as a starting scanning frequency of the sub-signal, a frequency multiplication time point corresponding to the starting scanning frequency is determined, and a scanning frequency corresponding to the frequency multiplication time point is determined as a terminal scanning frequency of the sub-signal, wherein the complete force signal is a force signal of a full frequency band corresponding to the first force signal.

[0024] In a possible implementation, the process of determining the plurality of frequency multiplication time points comprises:

[0025] For each sub-signal, a starting scanning frequency and a target scanning frequency of the sub-signal are determined, the target scanning frequency is less than twice the starting scanning frequency, and the target scanning frequency is a scanning frequency corresponding to a preset energy value;

[0026] An output parameter corresponding to a scanning frequency of each frequency sampling point between the starting scanning frequency and the target scanning frequency of the sub-signal is determined;

[0027] Based on the output parameters corresponding to the scanning frequencies of the plurality of frequency sampling points, a target time period required to reach the preset energy value is determined through energy relationship data, wherein the energy relationship data is relationship data between the output parameters corresponding to the scanning frequencies of the plurality of frequency sampling points and the target time period;

[0028] A sum of a starting time point of the starting scanning frequency and the target time period is determined as a frequency multiplication time point of the sub-signal.

[0029] In a possible implementation, the determining of the target frequency band in the second force signal based on the effective scanning frequency point number and the plurality of sub-signals comprises:

[0030] A plurality of effective signals in the plurality of sub-signals are determined, the plurality of effective signals are a plurality of sub-signals corresponding to a number of the effective scanning frequency point number after a starting scanning frequency of the second force signal;

[0031] A frequency band between the starting scanning frequency of the second force signal and a terminal scanning frequency of a last effective signal in the plurality of effective signals is determined as the target frequency band.

[0032] In a possible implementation, the determining of the effective scanning frequency point number of the second force signal comprises:

[0033] A difference between a terminal scanning frequency and a starting scanning frequency of the second force signal is determined as a scanning frequency width of the second force signal.

[0034] determine a quotient of the scanning bandwidth and a frequency sampling interval of a complete force signal as the effective scanning frequency point number, the complete force signal being a force signal of a full frequency band corresponding to the first force signal.

[0035] In a possible implementation, the quality parameters include fundamental wave energy, harmonic wave energy, signal-to-noise ratio, and lost energy; and the obtaining of the quality parameters of the controllable vibrator based on the first force signal and the fundamental wave signal includes:

[0036] determining a difference between the first force signal and the fundamental wave signal as a harmonic wave signal in the first force signal;

[0037] determining fundamental wave energy and harmonic wave energy corresponding to the fundamental wave signal and the harmonic wave signal, respectively;

[0038] determining a quotient of the fundamental wave energy and the harmonic wave energy as the signal-to-noise ratio;

[0039] obtaining a preset energy of a scanning signal corresponding to the first force signal;

[0040] determining a quotient of the preset energy and the fundamental wave energy as the lost energy.

[0041] In another aspect, a device for determining quality parameters of a controllable vibrator is provided, and the device includes:

[0042] a first determining module configured to determine a first force signal and a second force signal, the first force signal being a time-domain force signal excited by the controllable vibrator, and the second force signal being a frequency-domain force signal corresponding to the first force signal;

[0043] a segmenting module configured to segment the second force signal to obtain a plurality of sub-signals, a scanning bandwidth of each sub-signal being less than one time of a starting scanning frequency of the sub-signal;

[0044] a second determining module configured to determine an effective scanning frequency point number of the second force signal;

[0045] a third determining module configured to determine a target frequency band in the second force signal based on the effective scanning frequency point number and the plurality of sub-signals, the target frequency band being a frequency band corresponding to a fundamental wave signal in the first force signal;

[0046] a converting module configured to perform time-domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal;

[0047] a fourth determining module configured to obtain quality parameters of the controllable vibrator based on the first force signal and the fundamental wave signal, the quality parameters being used to determine a quality of the first force signal.

[0048] In a possible implementation, the segmenting module comprises:

[0049] an acquisition unit, configured to acquire a scanning signal corresponding to the first force signal;

[0050] a first segmenting unit, configured to, in a case where the scanning signal is a linear signal, segment the second force signal based on a frequency sampling interval of a complete force signal to obtain the plurality of sub-signals, the complete force signal being a full-band force signal corresponding to the first force signal;

[0051] a second segmenting unit, configured to, in a case where the scanning signal is a nonlinear signal, determine a plurality of frequency multiplication time points, and segment the second force signal based on the plurality of frequency multiplication time points to obtain the plurality of sub-signals.

[0052] In a possible implementation, the first segmenting unit is configured to:

[0053] acquire a starting scanning frequency of the first force signal;

[0054] in a case where the sub-signal is a first sub-signal, determine that the starting scanning frequency of the first force signal is a starting scanning frequency of the sub-signal, and determine a difference between twice the starting scanning frequency of the first force signal and the frequency sampling interval as a terminal scanning frequency of the sub-signal;

[0055] in a case where the sub-signal is not the first sub-signal, determine a sum of a terminal scanning frequency of a previous sub-signal and the frequency sampling interval as a starting scanning frequency of the sub-signal, and determine a difference between twice the starting scanning frequency of the sub-signal and the frequency sampling interval as a terminal scanning frequency of the sub-signal.

[0056] In a possible implementation, the second segmenting unit is configured to:

[0057] acquire a starting scanning frequency of the first force signal;

[0058] in a case where the sub-signal is a first sub-signal, determine that the starting scanning frequency of the first force signal is a starting scanning frequency of the sub-signal, and determine a scanning frequency corresponding to a first frequency multiplication time point as a terminal scanning frequency of the sub-signal;

[0059] in a case where the sub-signal is a first sub-signal, determine the starting scanning frequency of the first force signal, and determine a scanning frequency corresponding to a first frequency multiplication time point as a terminal scanning frequency of the sub-signal;

[0060] In a case that the sub-signal is not the first sub-signal, a time sampling interval of a complete force signal is obtained, a sum of the time sampling interval and a last frequency multiplication time point is determined as a starting time point of the sub-signal, a scanning frequency corresponding to the starting time point is determined as a starting scanning frequency of the sub-signal, a frequency multiplication time point corresponding to the starting scanning frequency is determined, and a scanning frequency corresponding to the frequency multiplication time point is taken as a terminal scanning frequency of the sub-signal, wherein the complete force signal is a force signal of a full frequency band corresponding to the first force signal.

[0061] In a possible implementation, the second segmentation unit is configured to:

[0062] For each sub-signal, a starting scanning frequency and a target scanning frequency of the sub-signal are determined, the target scanning frequency is less than twice the starting scanning frequency, and the target scanning frequency is a scanning frequency corresponding to a preset energy value;

[0063] An output parameter corresponding to a scanning frequency of each frequency sampling point between the starting scanning frequency and the target scanning frequency of the sub-signal is determined.

[0064] Based on the output parameters corresponding to the scanning frequencies of the plurality of frequency sampling points, a target time period required to reach the preset energy value is determined through energy relationship data, wherein the energy relationship data is relationship data between the output parameters corresponding to the scanning frequencies of the plurality of frequency sampling points and the target time period.

[0065] A sum of a starting time point of the starting scanning frequency and the target time period is determined as a frequency multiplication time point of the sub-signal.

[0066] In a possible implementation, the third determination module is configured to:

[0067] A plurality of effective signals in the plurality of sub-signals are determined, the plurality of effective signals are a plurality of sub-signals corresponding to a number of the effective scanning frequency points after the starting scanning frequency of the second force signal.

[0068] A frequency band between the starting scanning frequency of the second force signal and a terminal scanning frequency of a last effective signal in the plurality of effective signals is taken as the target frequency band.

[0069] In a possible implementation, the second determination module is configured to:

[0070] A difference between the terminal scanning frequency and the starting scanning frequency of the second force signal is determined as a scanning frequency width of the second force signal.

[0071] Determine a quotient of the scanning bandwidth and a frequency sampling interval of a complete force signal as the effective scanning frequency point number, the complete force signal being a force signal of a full frequency band corresponding to the first force signal.

[0072] In a possible implementation, the quality parameters include fundamental wave energy, harmonic wave energy, signal-to-noise ratio, and lost energy; and the fourth determining module is configured to:

[0073] Determine a difference between the first force signal and the fundamental wave signal as a harmonic wave signal in the first force signal.

[0074] Determine a fundamental wave energy and a harmonic wave energy corresponding to the fundamental wave signal and the harmonic wave signal, respectively.

[0075] Determine a quotient of the fundamental wave energy and the harmonic wave energy as the signal-to-noise ratio.

[0076] Obtain a preset energy of a scanning signal corresponding to the first force signal.

[0077] Determine a quotient of the preset energy and the fundamental wave energy as the lost energy.

[0078] In another aspect, a computer device is provided, which includes one or more processors and one or more memories, the one or more memories storing at least one instruction, the at least one instruction being loaded and executed by the one or more processors to implement operations performed by the method for determining quality parameters of a controllable source described in any of the implementations.

[0079] In another aspect, a computer readable storage medium is provided, which stores at least one instruction, the at least one instruction being loaded and executed by a processor to implement operations performed by the method for determining quality parameters of a controllable source described in any of the implementations.

[0080] In another aspect, a computer program product or a computer program is provided, which includes computer program code stored in a computer readable storage medium. A processor of a computer device reads the computer program code from the computer readable storage medium, and the processor executes the computer program code, so that the computer device performs operations performed by the method for determining quality parameters of a controllable source described above.

[0081] The technical scheme provided by the embodiments of the present application has at least the following beneficial effects:

[0082] The embodiment of the present application provides a method for determining a quality parameter of a controllable source, since the method segments the second force signal in the frequency domain, the frequency band corresponding to the fundamental wave signal in the first force signal in the time domain can be obtained, and then the separation of the fundamental wave signal and the harmonic wave signal in the first force signal can be realized; in this way, the quality parameter of the controllable source is obtained through the fundamental wave signal and the harmonic wave signal in the first force signal, the quantitative determination of the quality of the first force signal can be realized, and therefore the accuracy of determining the quality of the force signal of the controllable source is improved. BRIEF DESCRIPTION OF DRAWINGS

[0083] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiment description will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative effort on the basis of these drawings.

[0084] Figure 1 is a flow chart of a method for determining a quality parameter of a controllable source provided by the embodiment of the present application;

[0085] Figure 2 is a schematic diagram of segmenting a force signal in the frequency domain provided by the embodiment of the present application;

[0086] Figure 3 is a schematic diagram of a scanning signal provided by the embodiment of the present application;

[0087] Figure 4 is a schematic diagram of a first force signal provided by the embodiment of the present application;

[0088] Figure 5 is a schematic diagram of a second force signal provided by the embodiment of the present application;

[0089] Figure 6 is a schematic diagram of a fundamental wave signal provided by the embodiment of the present application;

[0090] Figure 7 is a schematic diagram of a harmonic wave signal provided by the embodiment of the present application;

[0091] Figure 8 is a schematic diagram of a quality parameter provided by the embodiment of the present application;

[0092] Figure 9 is a block diagram of a device for determining a quality parameter of a controllable source provided by the embodiment of the present application;

[0093] Figure 10 is a block diagram of a computer device provided by the embodiment of the present application. DETAILED DESCRIPTION

[0094] For the purpose of making the purpose, technical solutions and advantages of the present application clearer, the embodiments of the present application will be described in further detail below with reference to the drawings.

[0095] The terms "first", "second", "third", and "fourth" and the like in the description and claims of the present application and the drawings are used to distinguish different objects, and are not used to describe a particular order. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device including a series of steps or units is not limited to the listed steps or units, but can optionally further include steps or units not listed, or can optionally further include other steps or units inherent to the process, method, product or device.

[0096] The embodiments of the present application provide a method for determining a quality parameter of a controllable source, referring to Figure 1 , the method comprises:

[0097] Step 101: The computer device determines a first force signal and a second force signal.

[0098] The first force signal is a time-domain force signal excited by the controllable source, and the second force signal is a frequency-domain force signal corresponding to the first force signal. The first force signal is a force signal of a complete force signal excited by the controllable source in a preset time period; the preset time period can be set and changed as needed.

[0099] The abscissa of the first force signal is time, and the ordinate is amplitude. The abscissa of the second force signal is frequency, and the ordinate is amplitude. Optionally, the first force signal is Fourier transformed to obtain the second force signal, and the second force signal is inverse Fourier transformed to obtain the first force signal.

[0100] Step 102: The computer device segments the second force signal to obtain a plurality of sub-signals.

[0101] The scanning width of each sub-signal is less than one time of the starting scanning frequency of the sub-signal. It should be noted that the harmonic signal of a general controllable source is a multiple frequency harmonic signal, that is, the frequency of the harmonic signal is an integer multiple of the frequency of the fundamental signal. The occurrence time of the multiple frequency harmonic signal in the time-domain force signal is obtained by the following formula one:

[0102] Formula one:

[0103] Wherein, F1 is the starting scanning frequency, T1 is the time when the multiple frequency harmonic of the starting scanning frequency appears, F2 is the terminal scanning frequency, T2 is the time when the multiple frequency harmonic of the terminal scanning frequency appears, T is the scanning length, that is, the time length of the force signal, W is the scanning width, and the superscript 2 indicates the second harmonic.

[0104] From the above formula one, no matter how many sampling points are taken in the time domain, the harmonic signal will always appear at the time point of the starting scanning frequency / scanning frequency width*scanning length. It is known that it is difficult to separate the fundamental signal and the harmonic signal in the time domain. However, for the frequency domain, since the harmonic frequency is the multiple of the fundamental frequency, the ending scanning frequency in the force signal is set to be less than the frequency at which the second harmonic appears, that is, the harmonic signal and the fundamental signal are separated in the frequency domain. Then, a suitable frequency band range can be selected to separate the fundamental signal and the harmonic signal in the force signal. See Figure 2 , Figure 2 The figure is a schematic diagram for segmenting the force signal in the frequency domain. The starting scanning frequency is 1.5 Hz, the scanning interval is 1.5-2.8 Hz, and then the frequency interval of the corresponding harmonic signal is 3-5.6 Hz.

[0105] In this step, the computer device obtains a scanning signal corresponding to the first force signal. The scanning signal is one of a linear signal and a nonlinear signal. Step 102 includes the following two implementation cases:

[0106] The first case: in the case where the scanning signal is a linear signal, the computer device segments the second force signal based on the frequency sampling interval of the complete force signal to obtain a plurality of sub-signals. The complete force signal is the force signal of the full frequency band corresponding to the first force signal.

[0107] The first implementation case includes the following steps (1)-(3):

[0108] (1) The computer device obtains the starting scanning frequency of the first force signal.

[0109] The starting scanning frequency is obtained from the seismic data of the controllable source.

[0110] (2) In the case where the sub-signal is the first sub-signal, the computer device determines that the starting scanning frequency of the first force signal is the starting scanning frequency of the sub-signal, and determines that the difference between twice the starting scanning frequency of the first force signal and the frequency sampling interval is the ending scanning frequency of the sub-signal.

[0111] The computer device obtains the frequency sampling interval by the following way:

[0112] The computer device obtains the frequency width and the number of sampling points of the complete force signal. The computer device determines the quotient of the frequency width and the number of sampling points of the complete force signal as the frequency sampling interval, see formula two:

[0113] Formula two:

[0114] Wherein, F is the frequency width of the complete force signal, SN is the number of sampling points of the complete force signal, and df is the frequency sampling interval.

[0115] It should be noted that the starting scanning frequency of the first sub-signal is the starting scanning frequency of the first force signal. Optionally, f 11 represents the starting scanning frequency of the first force signal, i.e. the starting scanning frequency of the first sub-signal, f 12 represents the ending scanning frequency of the first sub-signal, then f 12 satisfies f 12 <2×f 11 .

[0116] (3) In the case where the sub-signal is not the first sub-signal, the computer device determines the sum of the ending scanning frequency of the last sub-signal and the frequency sampling interval as the starting scanning frequency of the sub-signal, and determines the difference between twice the starting scanning frequency of the sub-signal and the frequency sampling interval as the ending scanning frequency of the sub-signal.

[0117] For example, in the case where the sub-signal is the next sub-signal of the first sub-signal, i.e. in the case where the sub-signal is the second sub-signal, the starting scanning frequency of the sub-signal is the sum of the ending scanning frequency of the first sub-signal and the frequency sampling interval, and the ending scanning frequency of the sub-signal is the difference between twice the starting scanning frequency of the sub-signal and the frequency sampling interval. Optionally, f 21 represents the starting scanning frequency of the second sub-signal, f 22 represents the ending scanning frequency of the second sub-signal, then f 21 =f 12 +df, f 22 <2×f 21 . In this way, the segmentation of the entire frequency band of the second force signal is completed.

[0118] The second case: In the case where the scanned signal is a nonlinear signal, the computer device determines a plurality of frequency multiplication time points, segments the second force signal based on the plurality of frequency multiplication time points, and obtains a plurality of sub-signals.

[0119] Among them, the computer device segments the second force signal based on the plurality of frequency multiplication time points to obtain a plurality of sub-signals, which is realized by the following way:

[0120] The computer device acquires the starting scan frequency of the first force signal. If the sub-signal is the first sub-signal, the computer device determines the starting scan frequency of the first force signal as the starting scan frequency of that sub-signal, and uses the scan frequency corresponding to the first harmonic time point as the ending scan frequency of that sub-signal. If the sub-signal is not the first sub-signal, the computer device acquires the time sampling interval of the complete force signal, determines the sum of the time sampling interval and the previous harmonic time point as the starting time point of that sub-signal, determines the scan frequency corresponding to the starting time point as the starting scan frequency of that sub-signal, determines the harmonic time point corresponding to the starting scan frequency, and uses the scan frequency corresponding to the harmonic time point as the ending scan frequency of that sub-signal.

[0121] The process of determining multiple frequency doubling time points using computer equipment includes steps A1-A4:

[0122] A1: For each sub-signal, the computer device determines the starting scan frequency and the target scan frequency of that sub-signal. The target scan frequency is less than twice the starting scan frequency and is the scan frequency corresponding to the preset energy value.

[0123] The preset energy value for each sub-signal is the energy value reached by that sub-signal in the force signal excited by the controllable seismic source. For example, the energy value reached by the sub-signal is 1000 kg.

[0124] A2: The computer equipment determines the output parameters corresponding to the scanning frequency of each frequency sampling point between the starting scanning frequency and the target scanning frequency of the sub-signal.

[0125] The output parameter is used to represent the energy value corresponding to any scanning frequency; optionally, for scanning frequencies of 1Hz and 2Hz, the corresponding output parameters are 100 kg and 150 kg, respectively.

[0126] A3: The computer equipment determines the target time period required to reach the preset energy value based on the output parameters corresponding to the scanning frequencies of multiple frequency sampling points and the energy relationship data. The energy relationship data is the relationship data between the output parameters corresponding to the scanning frequencies of multiple frequency sampling points and the target time period.

[0127] The energy relationship data are as follows:

[0128]

[0129] Among them, f min f is the starting scan frequency. i Let η be the target scanning frequency, and η be the frequency. Here, PSD represents the output parameter, and DF represents the energy density. 2 For the output value, t i The target time period.

[0130] It should be noted that since the scanning signals corresponding to the first force signal and the second force signal are nonlinear signals, the output force parameters corresponding to each scanning frequency are irregular. Therefore, it is necessary to accumulate and sum the output force parameters of each scanning frequency to obtain the frequency doubling time point.

[0131] A4: The sum of the starting time point and the target time period for determining the starting scan frequency by the computer equipment is used as the frequency doubling time point of the sub-signal.

[0132] In this implementation, the target time period and frequency doubling time point of each sub-signal are obtained by summing the output parameters of the scanning frequency. Then, based on the frequency doubling time point, the start scanning frequency and end scanning frequency of each sub-signal are determined.

[0133] Step 103: The computer device determines the number of effective scanning frequency points of the second force signal.

[0134] Step 103 includes the following steps (1)-(2):

[0135] (1) The computer equipment determines the difference between the termination scanning frequency and the start scanning frequency of the second force signal as the scanning bandwidth of the second force signal.

[0136] It should be noted that the scanning bandwidth of the second force signal is the same as that of the first force signal.

[0137] (2) The computer equipment determines the ratio of the scanning bandwidth to the frequency sampling interval of the complete force signal as the effective number of scanning frequency points.

[0138] The effective number of scanning frequency points is obtained using Formula 3:

[0139] Formula 3:

[0140] Among them, f 12 f is the termination scan frequency of the second force signal. 11 f is the starting scan frequency of the second force signal. 12 with f 11 The difference is the scanning bandwidth, df is the frequency sampling interval, and n eff The effective number of scanning frequency points.

[0141] Step 104: The computer device determines the target frequency band in the second force signal based on the number of effective scan frequency points and multiple sub-signals.

[0142] The target frequency band is the frequency band corresponding to the fundamental wave signal in the first force signal. Step 104 includes the following steps (1)-(2):

[0143] (1) The computer device determines a plurality of valid signals in the plurality of sub-signals, the plurality of valid signals being a plurality of sub-signals corresponding to a number of valid scanning frequency points after a starting scanning frequency of the second force signal.

[0144] For example, if the number of valid scanning frequency points is 10, the plurality of valid signals are the first 10 sub-signals after the starting scanning frequency of the second force signal.

[0145] (2) The computer device takes a frequency band between the starting scanning frequency of the second force signal and the ending scanning frequency of the last valid signal in the plurality of valid signals as a target frequency band.

[0146] In the embodiment of the present application, since the plurality of valid signals are sub-signals corresponding to the number of valid frequency points, that is, the plurality of valid signals correspond to the fundamental wave signal; in this way, the frequency band corresponding to the plurality of valid signals in the second force signal is taken as the target frequency band, that is, the separation of the fundamental wave signal and the harmonic signal in the second force signal is realized, and the fundamental wave signal is obtained.

[0147] Step 105: The computer device performs time domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal.

[0148] It should be noted that the computer device uses the Fourier inverse transform method to perform time domain conversion on the target frequency band, and the force signal in the time domain obtained is the fundamental wave signal in the first force signal.

[0149] In the embodiment of the present application, since the second force signal in the frequency domain is segmented at the frequency multiplication point, the separation of the fundamental wave signal and the harmonic signal is realized, and then the fundamental wave signal in the frequency domain is converted to the time domain, that is, the fundamental wave signal in the time domain is obtained, and the effective separation of the fundamental wave signal and the harmonic signal in the time domain is realized.

[0150] Step 106: The computer device obtains a quality parameter of the controllable source based on the first force signal and the fundamental wave signal,

[0151] The quality parameter is used to determine the quality of the first force signal. The quality parameter includes the fundamental wave energy, the harmonic energy, the signal-to-noise ratio and the lost energy, and step 106 includes the following steps (1)-(5):

[0152] (1) The computer device determines the difference between the first force signal and the fundamental wave signal as the harmonic signal in the first force signal.

[0153] Referring to Figure 3 , Figure 3 is a schematic diagram of a certain segment of the scanning signal obtained by segmenting the complete scanning signal. Referring to Figure 4 , Figure 4 is the first force signal corresponding to the scanning signal, and comparing Figure 3 and Figure 4It can be seen that there is a difference between the scanning signal and the first force signal. See Figure 5 , Figure 5 is a second force signal corresponding to the first force signal. See Figure 6 , Figure 6 is a fundamental wave signal in the first force signal obtained by the method provided by the embodiment of the present application. See Figure 7 , Figure 7 is a harmonic wave signal in the first force signal obtained based on the difference between the first force signal and the fundamental wave signal.

[0154] (2) The computer device determines the fundamental wave energy and the harmonic wave energy corresponding to the fundamental wave signal and the harmonic wave signal respectively.

[0155] The fundamental wave energy refers to the energy of the fundamental wave component in the first force signal, and is the effective energy of the quantitative display output. The harmonic wave energy refers to the energy of the harmonic wave component in the first force signal, and is the interference energy of the quantitative display output. The computer device extracts data sequences of the fundamental wave signal and the harmonic wave signal respectively, and performs norm calculation on the data sequences of the fundamental wave signal and the harmonic wave signal respectively to obtain the fundamental wave energy and the harmonic wave energy.

[0156] The fundamental wave energy is E f , and the harmonic wave energy is E d . s is the extracted fundamental wave data sequence, d is the extracted harmonic wave data sequence, E f is the fundamental wave energy, and E d is the harmonic wave energy.

[0157] (3) The computer device determines the quotient of the fundamental wave energy and the harmonic wave energy as the signal-to-noise ratio.

[0158] The signal-to-noise ratio refers to the ratio of the fundamental wave component to the harmonic wave component interference in the first force signal, that is, the ratio of the fundamental wave energy to the harmonic wave energy, which is used to quantitatively reflect the response accuracy of the controllable seismic source system to a given scanning signal.

[0159] The signal-to-noise ratio is obtained by the following formula four:

[0160] Formula four:

[0161] The signal-to-noise ratio is R f , the fundamental wave energy is E f , and the harmonic wave energy is E d .

[0162] (4) The computer device obtains a preset energy of the scanning signal corresponding to the first force signal.

[0163] The computer device extracts a data sequence of the scanning signal, and performs norm calculation on the data sequence to obtain the preset energy of the scanning signal.

[0164] Wherein, E=||S||2, S is a data sequence of the scanning signal corresponding to the first force signal, and E is a preset energy.

[0165] (5) The computer device determines a quotient of the preset energy and the fundamental wave energy as the loss energy.

[0166] Wherein, the loss energy refers to a difference between the fundamental wave component in the first force signal and the corresponding scanning signal component energy, that is, a quotient of the preset energy and the fundamental wave energy, and is used for quantitatively reflecting an efficiency of the controllable seismic source system in transforming effective energy.

[0167] The loss energy is obtained through the following Formula Five.

[0168] Formula Five:

[0169] Wherein, D i / o is the loss energy, E is the preset energy, and E f is the fundamental wave energy.

[0170] It should be noted that the loss energy is low, indicating that the energy conversion rate of the preset scanning signal is high. The force signal quality is good with a high signal-to-noise ratio, a large output fundamental wave energy value, and a small harmonic energy value. See Figure 8 , Figure 8 The quality parameter is obtained through the method provided in the embodiments of the present application, and the fundamental wave energy, the harmonic energy, the loss energy, and the signal-to-noise ratio are displayed from left to right. Wherein, the quality of the force signal of different types of seismic sources and different exploration areas can be evaluated through the quality parameter obtained through the method provided in the embodiments of the present application; and the evaluation standard of the force signal itself can be improved, which is conducive to the optimization of the excitation method and the landing of high-precision seismic exploration, thereby providing an evaluation basis for the improvement of the controllable seismic source itself and the development of high-precision scanning signals.

[0171] The embodiments of the present application provide a method for determining a controllable seismic source quality parameter. Since the method segments the second force signal in the frequency domain, the frequency band corresponding to the fundamental wave signal in the first force signal in the time domain can be obtained, and the separation of the fundamental wave signal and the harmonic signal in the first force signal can be realized. In this way, the quality parameter of the controllable seismic source is obtained through the fundamental wave signal and the harmonic signal in the first force signal, and the quantitative determination of the quality of the first force signal can be realized, thereby improving the accuracy of determining the quality of the force signal of the controllable seismic source.

[0172] The embodiments of the present application also provide a device for determining a controllable seismic source quality parameter. See Figure 9 , the device comprises:

[0173] The first determining module 901 is configured to determine a first force signal and a second force signal, the first force signal being a time-domain force signal excited by a controllable vibrator, and the second force signal being a frequency-domain force signal corresponding to the first force signal;

[0174] The segmenting module 902 is configured to segment the second force signal to obtain a plurality of sub-signals, a scanning frequency width of each sub-signal being less than one time of a starting scanning frequency of the sub-signal;

[0175] The second determining module 903 is configured to determine an effective scanning frequency point number of the second force signal;

[0176] The third determining module 904 is configured to determine a target frequency band in the second force signal based on the effective scanning frequency point number and the plurality of sub-signals, the target frequency band being a frequency band corresponding to a fundamental wave signal in the first force signal;

[0177] The converting module 905 is configured to perform time-domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal;

[0178] The fourth determining module 906 is configured to obtain a quality parameter of the controllable vibrator based on the first force signal and the fundamental wave signal, the quality parameter being used to determine a quality of the first force signal.

[0179] In a possible implementation, the segmenting module 902 includes:

[0180] The acquisition unit is configured to acquire a scanning signal corresponding to the first force signal;

[0181] The first segmenting unit is configured to, in a case where the scanning signal is a linear signal, segment the second force signal based on a frequency sampling interval of a complete force signal to obtain the plurality of sub-signals, the complete force signal being a force signal of a full frequency band corresponding to the first force signal;

[0182] The second segmenting unit is configured to, in a case where the scanning signal is a nonlinear signal, determine a plurality of frequency multiplication time points, and segment the second force signal based on the plurality of frequency multiplication time points to obtain the plurality of sub-signals.

[0183] In a possible implementation, the first segmenting unit is configured to:

[0184] acquire a starting scanning frequency of the first force signal;

[0185] in a case where the sub-signal is a first sub-signal, determine that the starting scanning frequency of the first force signal is a starting scanning frequency of the sub-signal, and determine a difference between twice the starting scanning frequency of the first force signal and the frequency sampling interval as a terminal scanning frequency of the sub-signal;

[0186] In a case where the sub-signal is not the first sub-signal, a sum of a terminal scanning frequency of the previous sub-signal and the frequency sampling interval is determined as a starting scanning frequency of the sub-signal, and a difference between twice the starting scanning frequency of the sub-signal and the frequency sampling interval is determined as a terminal scanning frequency of the sub-signal.

[0187] In a possible implementation, the second segmentation unit is configured to:

[0188] obtain a starting scanning frequency of the first force signal;

[0189] In a case where the sub-signal is the first sub-signal, the starting scanning frequency of the first force signal is determined as the starting scanning frequency of the sub-signal, and a scanning frequency corresponding to the first frequency multiplication time point is determined as the terminal scanning frequency of the sub-signal.

[0190] In a case where the sub-signal is not the first sub-signal, a time sampling interval of a complete force signal is obtained, a sum of the time sampling interval and a previous frequency multiplication time point is determined as a starting time point of the sub-signal, a scanning frequency corresponding to the starting time point is determined as the starting scanning frequency of the sub-signal, a frequency multiplication time point corresponding to the starting scanning frequency is determined, and a scanning frequency corresponding to the frequency multiplication time point is determined as the terminal scanning frequency of the sub-signal, the complete force signal being a force signal of a full frequency range corresponding to the first force signal.

[0191] In a possible implementation, the second segmentation unit is configured to:

[0192] For each sub-signal, a starting scanning frequency of the sub-signal and a target scanning frequency are determined, the target scanning frequency being less than twice the starting scanning frequency, and the target scanning frequency being a scanning frequency corresponding to a preset energy value;

[0193] determining, for each frequency sampling point between the starting scanning frequency and the target scanning frequency of the sub-signal, a scanning frequency corresponding to an output parameter;

[0194] based on the output parameters corresponding to the scanning frequencies of the plurality of frequency sampling points, determining, through energy relationship data, a target time period required to reach the preset energy value, the energy relationship data being relationship data between the output parameters corresponding to the scanning frequencies of the plurality of frequency sampling points and the target time period;

[0195] determining, as a frequency multiplication time point of the sub-signal, a sum of a starting time point of the starting scanning frequency and the target time period.

[0196] In a possible implementation, the third determination module 904 is configured to:

[0197] determining a plurality of effective signals in the plurality of sub-signals, the plurality of effective signals being a plurality of sub-signals corresponding to a number of effective scanning frequency points after the starting scanning frequency of the second force signal;

[0198] A frequency band between the start scanning frequency of the second force signal and the end scanning frequency of the last one of the plurality of valid signals is taken as a target frequency band.

[0199] In a possible implementation, the second determining module 903 is configured to:

[0200] determine a difference between the end scanning frequency and the start scanning frequency of the second force signal as a scanning frequency width of the second force signal.

[0201] determine a quotient of the scanning frequency width and a frequency sampling interval of a complete force signal as a valid scanning frequency point number, the complete force signal being a force signal corresponding to a full frequency band of the first force signal.

[0202] In a possible implementation, the quality parameter includes fundamental wave energy, harmonic wave energy, signal-to-noise ratio, and lost energy; and the fourth determining module 906 is configured to:

[0203] determine a difference between the first force signal and the fundamental wave signal as a harmonic wave signal in the first force signal.

[0204] respectively determine fundamental wave energy and harmonic wave energy corresponding to the fundamental wave signal and the harmonic wave signal.

[0205] determine a quotient of the fundamental wave energy and the harmonic wave energy as a signal-to-noise ratio.

[0206] obtain a preset energy of a scanning signal corresponding to the first force signal.

[0207] determine a quotient of the preset energy and the fundamental wave energy as a lost energy.

[0208] Figure 10 A structural block diagram of a computer device 1000 provided by an example embodiment of the present application is shown. The computer device 1000 can be a portable mobile computer device, such as a smart phone, a tablet computer, an MP3 player, an MP4 player, a notebook computer, or a desktop computer. The computer device 1000 can also be referred to as a user device, a portable computer device, a laptop computer device, a desktop computer device, or other names.

[0209] Generally, the computer device 1000 includes a processor 1001 and a memory 1002.

[0210] The processor 1001 can include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 1001 can be implemented in at least one of a hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), a PLA (Programmable Logic Array). The processor 1001 can also include a main processor and a coprocessor, the main processor being a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit), and the coprocessor being a low-power processor for processing data in a standby state. In some embodiments, the processor 1001 can be integrated with a GPU (Graphics Processing Unit) for rendering and drawing content required to be displayed by the display screen. In some embodiments, the processor 1001 can further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

[0211] The memory 1002 can include one or more computer-readable storage media that can be non-transitory. The memory 1002 can also include high-speed random access memory and nonvolatile, computer-readable storage media such as one or more magnetic disk storage devices, flash memory devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 1002 is used to store at least one instruction for being executed by the processor 1001 to implement the method for determining the quality parameter of the seismic source provided by the method embodiments in the present application.

[0212] In some embodiments, the computer device 1000 can also optionally include a peripheral device interface 1003 and at least one peripheral device. The processor 1001, the memory 1002, and the peripheral device interface 1003 can be connected through a bus or a signal line. Each peripheral device can be connected to the peripheral device interface 1003 through a bus, a signal line, or a circuit board. Specifically, the peripheral device includes at least one of a radio frequency circuit 1004, a display screen 1005, a camera assembly 1006, an audio circuit 1007, a positioning assembly 1008, and a power supply 1009.

[0213] The peripheral interface 1003 can be used to connect at least one I / O (Input / Output) related peripheral device to the processor 1001 and the memory 1002. In some embodiments, the processor 1001, the memory 1002 and the peripheral interface 1003 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 1001, the memory 1002 and the peripheral interface 1003 can be implemented on a separate chip or circuit board, and the present embodiments are not limited in this regard.

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

[0215] The display screen 1005 is configured to display a UI (User Interface). The UI can include graphics, text, icons, video, and any combination thereof. When the display screen 1005 is a touch display screen, the display screen 1005 is further configured to capture touch signals on or above the surface of the display screen 1005. The touch signals can be input to the processor 1001 as control signals for processing. In this case, the display screen 1005 can also be configured to provide virtual buttons and / or virtual keyboard, also known as soft buttons and / or soft keyboard. In some embodiments, the display screen 1005 can be one, disposed on the front panel of the computer device 1000; in other embodiments, the display screen 1005 can be at least two, respectively disposed on different surfaces of the computer device 1000 or in a folding design; in other embodiments, the display screen 1005 can be a flexible display screen, disposed on a curved surface or a folding surface of the computer device 1000. Even, the display screen 1005 can also be disposed in an irregular shape other than a rectangle, i.e., a special-shaped screen. The display screen 1005 can be made of materials such as LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), etc.

[0216] The camera assembly 1006 is configured to capture images or videos. Optionally, the camera assembly 1006 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is disposed on the front panel of the computer device, and the rear-facing camera is disposed on the back of the computer device. In some embodiments, the rear-facing camera is at least two, which are any one of a main camera, a depth-of-field camera, a wide-angle camera, and a telephoto camera, to realize the background blur function by fusing the main camera and the depth-of-field camera, the panoramic shooting and VR (Virtual Reality) shooting function by fusing the main camera and the wide-angle camera, or other fusion shooting functions. In some embodiments, the camera assembly 1006 can further include a flash. The flash can be a single-color-temperature flash or a dual-color-temperature flash. The dual-color-temperature flash refers to a combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0217] The audio circuit 1007 can include a microphone and a speaker. The microphone is used to collect sound waves of a user and an environment, and convert the sound waves into an electrical signal input to the processor 1001 for processing, or input to the radio frequency circuit 1004 to realize voice communication. The microphone can be multiple for the purpose of stereo sound collection or noise reduction, and arranged at different parts of the computer device 1000. The microphone can also be an array microphone or an omnidirectional collection microphone. The speaker is used to convert an electrical signal from the processor 1001 or the radio frequency circuit 1004 into sound waves. The speaker can be a traditional diaphragm speaker, or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only can it convert an electrical signal into a sound wave audible to humans, but also can convert an electrical signal into an inaudible sound wave to humans for ranging purposes. In some embodiments, the audio circuit 1007 can also include a headphone jack.

[0218] The positioning component 1008 is used to position the current geographic location of the computer device 1000 to realize navigation or LBS (Location Based Service). The positioning component 1008 can be a positioning component based on the GPS (Global Positioning System) of the United States, the Beidou system of China, or the Galileo system of Russia.

[0219] The power supply 1009 is used to supply power to various components in the computer device 1000. The power supply 1009 can be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the power supply 1009 includes a rechargeable battery, the rechargeable battery can be a wired charging battery or a wireless charging battery. The wired charging battery is a battery charged through a wired line, and the wireless charging battery is a battery charged through a wireless coil. The rechargeable battery can also be used to support fast charging technology.

[0220] In some embodiments, the computer device 1000 further includes one or more sensors 1010. The one or more sensors 1010 include, but are not limited to, an acceleration sensor 1011, a gyroscope sensor 1012, a pressure sensor 1013, a fingerprint sensor 1014, an optical sensor 1015, and a proximity sensor 1016.

[0221] The acceleration sensor 1011 can detect the acceleration magnitude in three coordinate axes of a coordinate system established by the computer device 1000. For example, the acceleration sensor 1011 can be used to detect the components of gravitational acceleration in three coordinate axes. The processor 1001 can control the display screen 1005 to display a user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 1011. The acceleration sensor 1011 can also be used for gaming or collection of user motion data.

[0222] The gyroscope sensor 1012 can detect the body direction and rotation angle of the computer device 1000, and the gyroscope sensor 1012 can cooperate with the acceleration sensor 1011 to collect the 3D action of the user on the computer device 1000. According to the data collected by the gyroscope sensor 1012, the processor 1001 can realize the following functions: action sensing (such as changing the UI according to the user's tilt operation), image stabilization when shooting, game control, and inertial navigation.

[0223] The pressure sensor 1013 can be arranged on the side frame of the computer device 1000 and / or the lower layer of the display screen 1005. When the pressure sensor 1013 is arranged on the side frame of the computer device 1000, the user's holding signal on the computer device 1000 can be detected, and the left and right hand recognition or shortcut operation can be performed by the processor 1001 according to the holding signal collected by the pressure sensor 1013. When the pressure sensor 1013 is arranged on the lower layer of the display screen 1005, the controllable control on the UI interface can be controlled by the processor 1001 according to the pressure operation of the user on the display screen 1005. The controllable control includes at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0224] The fingerprint sensor 1014 is used to collect the fingerprint of the user, and the identity of the user can be recognized by the processor 1001 according to the fingerprint collected by the fingerprint sensor 1014, or by the fingerprint sensor 1014 according to the collected fingerprint. When the identity of the user is recognized as a trusted identity, the processor 1001 authorizes the user to perform related sensitive operations, including unlocking the screen, viewing encrypted information, downloading software, payment, and changing settings. The fingerprint sensor 1014 can be arranged on the front, back or side of the computer device 1000. When the computer device 1000 is provided with a physical button or a manufacturer's logo, the fingerprint sensor 1014 can be integrated with the physical button or the manufacturer's logo.

[0225] The optical sensor 1015 is used to collect the ambient light intensity. In one embodiment, the processor 1001 can control the display brightness of the display screen 1005 according to the ambient light intensity collected by the optical sensor 1015. Specifically, when the ambient light intensity is high, the display brightness of the display screen 1005 is increased; when the ambient light intensity is low, the display brightness of the display screen 1005 is decreased. In another embodiment, the processor 1001 can also dynamically adjust the shooting parameters of the camera assembly 1006 according to the ambient light intensity collected by the optical sensor 1015.

[0226] The proximity sensor 1016, also referred to as a distance sensor, is usually arranged on the front panel of the computer device 1000. The proximity sensor 1016 is used to collect the distance between the user and the front face of the computer device 1000. In an embodiment, when the proximity sensor 1016 detects that the distance between the user and the front face of the computer device 1000 gradually decreases, the display screen 1005 is switched from the bright screen state to the screen-off state under the control of the processor 1001; when the proximity sensor 1016 detects that the distance between the user and the front face of the computer device 1000 gradually increases, the display screen 1005 is switched from the screen-off state to the bright screen state under the control of the processor 1001.

[0227] Those skilled in the art can understand that, Figure 10 The structure shown in the figure is not a limitation on the computer device 1000, and can include more or fewer components than shown, or combine certain components, or use a different arrangement of components.

[0228] The embodiment of the present application further provides a computer readable storage medium, at least one instruction is stored in the computer readable storage medium, and the at least one instruction is loaded and executed by a processor to implement the operations performed by the method for determining the quality parameter of the force signal of the controllable source.

[0229] The embodiment of the present application further provides a computer program product or a computer program, the computer program product or the computer program includes computer program code, and the computer program code is stored in a computer readable storage medium. The processor of the computer device reads the computer program code from the computer readable storage medium, and the processor executes the computer program code, so that the computer device executes the operations performed by the method for determining the quality parameter of the force signal of the controllable source.

[0230] In some embodiments, the computer program related to the embodiment of the present application can be deployed to execute on one computer device, or execute on multiple computer devices located in one place, or execute on multiple computer devices distributed in multiple places and interconnected through a communication network, and the multiple computer devices distributed in multiple places and interconnected through a communication network can constitute a blockchain system.

[0231] The embodiment of the present application provides a method for determining the quality parameter of the controllable source. Since the method segments the second force signal in the frequency domain, the frequency band corresponding to the fundamental wave signal in the first force signal in the time domain can be obtained, and the separation of the fundamental wave signal and the harmonic signal in the first force signal can be realized. In this way, the quality parameter of the controllable source can be obtained through the fundamental wave signal and the harmonic signal in the first force signal, and the quantitative determination of the quality of the first force signal can be realized, thereby improving the accuracy of determining the quality of the force signal of the controllable source.

[0232] The above merely preferred embodiments of the present application and are not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method of determining a quality parameter of a vibroseis, characterized in that, The method comprises: determining a first force signal and a second force signal, the first force signal being a time domain force signal of a vibroseis excitation, and the second force signal being a frequency domain force signal corresponding to the first force signal; obtaining a scanning signal corresponding to the first force signal; in a case where the scanning signal is a linear signal, segmenting the second force signal based on a frequency sampling interval of a complete force signal to obtain a plurality of sub-signals, the complete force signal being a full-band force signal corresponding to the first force signal; in a case where the scanning signal is a nonlinear signal, determining a plurality of frequency multiplication time points, segmenting the second force signal based on the plurality of frequency multiplication time points to obtain a plurality of sub-signals, and a scanning bandwidth of each sub-signal being less than one time of a starting scanning frequency of the sub-signal; determining an effective scanning frequency point number of the second force signal; determining a plurality of effective signals in the plurality of sub-signals, the plurality of effective signals being a plurality of sub-signals corresponding to the effective scanning frequency point number after the starting scanning frequency of the second force signal; and taking a frequency band between the starting scanning frequency of the second force signal and an ending scanning frequency of the last effective signal in the plurality of effective signals as a target frequency band in the second force signal, the target frequency band being a frequency band corresponding to a fundamental wave signal in the first force signal; performing time domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal; obtaining a quality parameter of the vibroseis based on the first force signal and the fundamental wave signal, the quality parameter being used to determine a quality of the first force signal; wherein the quality parameter comprises a fundamental wave energy, a harmonic wave energy, a signal-to-noise ratio, and a loss energy; and the obtaining of the quality parameter of the vibroseis based on the first force signal and the fundamental wave signal comprises: determining a difference between the first force signal and the fundamental wave signal as a harmonic wave signal in the first force signal; determining the fundamental wave energy and the harmonic wave energy corresponding to the fundamental wave signal and the harmonic wave signal respectively; determining a quotient of the fundamental wave energy and the harmonic wave energy as the signal-to-noise ratio; obtaining a preset energy of a scanning signal corresponding to the first force signal; and determining a quotient of the preset energy and the fundamental wave energy as the loss energy.

2. The method of claim 1, wherein, The segmenting of the second force signal based on the frequency sampling interval of the complete force signal to obtain a plurality of sub-signals comprises: obtaining a starting scanning frequency of the first force signal; in a case where the sub-signal is a first sub-signal, determining the starting scanning frequency of the first force signal as a starting scanning frequency of the sub-signal, and determining a difference between two times of the starting scanning frequency of the first force signal and the frequency sampling interval as an ending scanning frequency of the sub-signal; in a case where the sub-signal is not the first sub-signal, determining a sum of the ending scanning frequency of a previous sub-signal and the frequency sampling interval as a starting scanning frequency of the sub-signal, and determining a difference between two times of the starting scanning frequency of the sub-signal and the frequency sampling interval as an ending scanning frequency of the sub-signal.

3. The method of claim 1, wherein, The second force signal is segmented based on the multiple frequency multiplication time points to obtain multiple sub-signals, including: Obtaining a starting scanning frequency of the first force signal; In the case that the sub-signal is the first sub-signal, determining the starting scanning frequency of the first force signal as the starting scanning frequency of the sub-signal, and taking the scanning frequency corresponding to the first frequency multiplication time point as the terminal scanning frequency of the sub-signal; In the case that the sub-signal is not the first sub-signal, obtaining a time sampling interval of a complete force signal, determining the sum of the time sampling interval and the previous frequency multiplication time point as the starting time point of the sub-signal, determining the scanning frequency corresponding to the starting time point as the starting scanning frequency of the sub-signal, determining the frequency multiplication time point corresponding to the starting scanning frequency, and taking the scanning frequency corresponding to the frequency multiplication time point as the terminal scanning frequency of the sub-signal, the complete force signal being a force signal of a full frequency band corresponding to the first force signal.

4. The method of claim 1, wherein, The process of determining the multiple frequency multiplication time points includes: For each sub-signal, determining the starting scanning frequency and a target scanning frequency of the sub-signal, the target scanning frequency being less than twice the starting scanning frequency, the target scanning frequency being a scanning frequency corresponding to a preset energy value; Determining the output parameter corresponding to the scanning frequency of each frequency sampling point between the starting scanning frequency and the target scanning frequency of the sub-signal; Based on the output parameter corresponding to the scanning frequency of the multiple frequency sampling points, determining the target time period required to reach the preset energy value through energy relationship data, the energy relationship data being relationship data between the output parameter corresponding to the scanning frequency of the multiple frequency sampling points and the target time period; Determining the sum of the starting time point of the starting scanning frequency and the target time period as the frequency multiplication time point of the sub-signal.

5. The method of claim 1, wherein, The determination of the effective scanning frequency point number of the second force signal includes: Determining the difference between the terminal scanning frequency and the starting scanning frequency of the second force signal as the scanning frequency width of the second force signal; Determining the quotient of the scanning frequency width and the frequency sampling interval of a complete force signal as the effective scanning frequency point number, the complete force signal being a force signal of a full frequency band corresponding to the first force signal.

6. An apparatus for determining a quality parameter of a vibroseis, the apparatus comprising: The device includes: A first determination module for determining a first force signal and a second force signal, the first force signal being a time-domain force signal excited by a controllable source, and the second force signal being a frequency-domain force signal corresponding to the first force signal; A segmentation module for obtaining a scanning signal corresponding to the first force signal; in the case that the scanning signal is a linear signal, segmenting the second force signal based on the frequency sampling interval of a complete force signal to obtain multiple sub-signals, the complete force signal being a force signal of a full frequency band corresponding to the first force signal; in the case that the scanning signal is a nonlinear signal, determining multiple frequency multiplication time points, and segmenting the second force signal based on the multiple frequency multiplication time points to obtain multiple sub-signals, the scanning frequency width of each sub-signal being less than one time of the starting scanning frequency of the sub-signal. a second determining module, configured to determine an effective scanning frequency point number of the second force signal; a third determining module, configured to determine a plurality of effective signals in the plurality of sub-signals, the plurality of effective signals being a plurality of sub-signals corresponding to the effective scanning frequency point number after a starting scanning frequency of the second force signal; taking a frequency band between the starting scanning frequency of the second force signal and a terminal scanning frequency of a last effective signal in the plurality of effective signals as a target frequency band in the second force signal, the target frequency band being a frequency band corresponding to a fundamental wave signal in the first force signal; a converting module, configured to perform time domain conversion on the target frequency band in the second force signal to obtain the fundamental wave signal in the first force signal; a fourth determining module, configured to obtain a quality parameter of the controllable seismic source based on the first force signal and the fundamental wave signal, the quality parameter being used to determine a quality of the first force signal; wherein the quality parameter includes a fundamental wave energy, a harmonic wave energy, a signal-to-noise ratio and a loss energy; the fourth determining module is configured to: determine a difference between the first force signal and the fundamental wave signal as a harmonic wave signal in the first force signal; respectively determine a fundamental wave energy and a harmonic wave energy corresponding to the fundamental wave signal and the harmonic wave signal; determine a quotient of the fundamental wave energy and the harmonic wave energy as the signal-to-noise ratio; obtain a preset energy of a scanning signal corresponding to the first force signal; and determine a quotient of the preset energy and the fundamental wave energy as the loss energy.

7. A computer device, comprising: The computer device includes one or more processors and one or more memories, the one or more memories storing at least one instruction, the at least one instruction being loaded and executed by the one or more processors to implement the operations performed by the method for determining a quality parameter of a controllable seismic source according to any one of claims 1 to 5.

Images