Method and device for determining harmonic influence range of vibroseis, equipment and medium

Abstract

The invention relates to a method, a device and equipment for determining the harmonic influence range of a vibroseis and a medium, and relates to the technical field of seismic acquisition data analysis, and the method comprises the steps: obtaining a scanning signal of the vibroseis, and carrying out the calculation to obtain the instantaneous frequency of the scanning signal; according to the instantaneous frequency, harmonic instantaneous frequency curves corresponding to harmonic signals of different sliding times and different orders are generated; generating a scanning signal instantaneous frequency curve corresponding to different dislocation time according to the instantaneous frequency; and determining the harmonic influence range of the vibroseis according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve. According to the scheme, the harmonic influence range can be conveniently and efficiently determined, and the method is suitable for the situation that a nonlinear scanning signal or a scanning signal cannot be accurately expressed through an analytical expression.

Description

Technical Field

[0001] The present disclosure relates to the technical field of seismic acquisition data analysis, and in particular to a method, device, equipment and medium for determining the influence range of controllable vibroseis harmonics. Background Art

[0002] Vibroseis has become the predominant excitation source for onshore seismic exploration due to its numerous advantages over explosive sources, including controllable frequency and energy, safety, environmental friendliness, low cost, and high efficiency. The traditional scanning method for vibroseis sources requires each source to wait for the previous source's scan to complete (including the vibration time and monitoring time, i.e., after receiving the echo) before starting, resulting in low efficiency. To improve data acquisition efficiency, a sliding scanning technique has been proposed during vibroseis data acquisition. This technique uses overlapping excitation (i.e., starting the next source scan before the previous one has completed) to improve acquisition efficiency. However, while sliding scanning improves efficiency, it inevitably introduces some interference, and harmonic interference is the most common type of interference in sliding scanning. Therefore, determining the impact range of harmonic interference is crucial to assist in optimizing construction parameters to suppress harmonic interference.

[0003] In the process of realizing the concept of the present disclosure, the inventors found that there are at least the following technical problems in the related technology: there are currently some empirical formulas to calculate the harmonic influence range of linear scanning, but the above empirical formulas are only applicable to calculating the harmonic influence range of linear scanning. When the scanning signal is not a linear scanning signal or the scanning signal cannot be accurately expressed by an analytical expression, the above empirical formula cannot be used to determine the influence range of harmonic interference. Summary of the Invention

[0004] In order to solve the above technical problems or at least partially solve the above technical problems, the embodiments of the present disclosure provide a method, apparatus, device and medium for determining the harmonic influence range of a controllable vibroseis source.

[0005] In a first aspect, embodiments of the present disclosure provide a method for determining the harmonic influence range of a vibrator. The method comprises: obtaining a scanning signal from the vibrator and calculating the instantaneous frequency of the scanning signal; generating, based on the instantaneous frequency, harmonic instantaneous frequency curves corresponding to different sliding times and harmonic signals of different orders; generating, based on the instantaneous frequency, scanning signal instantaneous frequency curves corresponding to different shift times; and determining the harmonic influence range of the vibrator based on the intersection of the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve.

[0006] In some embodiments, the harmonic influence range of the controllable vibrator is determined based on the intersection relationship between the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve, including: in response to the above-mentioned intersection relationship indicating that the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve have an intersection position, determining the intersection start time and the intersection end time according to the intersection position; determining the interference start time and the interference end time of the corresponding order harmonic signal according to the above-mentioned intersection start time, the intersection end time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve; and determining the time period corresponding to the above-mentioned interference start time and the above-mentioned interference end time as the harmonic influence range of the above-mentioned controllable vibrator at the corresponding order and the corresponding sliding time.

[0007] In some embodiments, the interference start time and interference end time of the corresponding order harmonic signal are determined based on the above-mentioned intersection start time, intersection end time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve, including: adding the above-mentioned intersection start time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve to obtain the interference start time of the corresponding order harmonic signal; adding the above-mentioned intersection end time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve to obtain the interference end time of the corresponding order harmonic signal.

[0008] In some embodiments, the harmonic influence range of the controllable seismic source is determined based on the intersection relationship between the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve, and also includes: in response to the above-mentioned intersection relationship indicating that the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve do not have an intersection position, determining that the above-mentioned controllable seismic source is not affected by harmonic interference within the corresponding offset time.

[0009] In some embodiments, calculating the instantaneous frequency of the scanning signal includes: obtaining the instantaneous frequency corresponding to the baseband signal from the scanning signal based on an analytical signal method or a time-frequency analysis method as the instantaneous frequency of the scanning signal.

[0010] In some embodiments, the sweep signal includes one of the following types: a nonlinear sweep signal and a sweep signal without an analytical expression. Based on the instantaneous frequency, generating harmonic instantaneous frequency curves corresponding to harmonic signals of different sliding times and orders includes: modifying the instantaneous frequency expression by multiplying the instantaneous frequency by the order coefficient, and introducing a sliding time term into the time parameter of the instantaneous frequency to represent the time shift present in the harmonic propagation process, thereby obtaining a harmonic instantaneous frequency expression related to the sliding time and order coefficient; and plotting the harmonic instantaneous frequency expression to obtain a harmonic instantaneous frequency curve.

[0011] In some embodiments, based on the above-mentioned instantaneous frequency, a harmonic instantaneous frequency curve corresponding to different sliding times and different-order harmonic signals is generated, and it also includes: receiving adjustment information for at least one of the above-mentioned order coefficient or sliding time item; adjusting the above-mentioned harmonic instantaneous frequency expression according to the above-mentioned adjustment information, and drawing based on the adjusted expression to obtain an updated harmonic instantaneous frequency curve.

[0012] In the second aspect, an embodiment of the present disclosure provides a device for determining the harmonic influence range of a controllable vibrator. The above-mentioned determination device includes: a signal acquisition module, an instantaneous frequency calculation module, a first curve generation module, a second curve generation module and an influence range determination module. The above-mentioned signal acquisition module is used to obtain the scanning signal of the controllable vibrator. The above-mentioned instantaneous frequency calculation module is used to calculate the instantaneous frequency of the above-mentioned scanning signal. The above-mentioned first curve generation module is used to generate harmonic instantaneous frequency curves corresponding to different sliding times and different orders of harmonic signals according to the above-mentioned instantaneous frequency. The above-mentioned second curve generation module is used to generate scanning signal instantaneous frequency curves corresponding to different offset times according to the above-mentioned instantaneous frequency. The above-mentioned influence range determination module is used to determine the harmonic influence range of the controllable vibrator according to the intersection relationship between the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve.

[0013] In a third aspect, embodiments of the present disclosure provide an electronic device. The electronic device includes a processor, a communication interface, a memory, and a communication bus. The processor, communication interface, and memory communicate with each other via the communication bus. The memory is configured to store a computer program. The processor is configured to implement the aforementioned method for determining the harmonic influence range of a vibrator when executing the program stored in the memory.

[0014] In a fourth aspect, embodiments of the present disclosure provide a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for determining the harmonic influence range of a vibroseis as described above.

[0015] The above technical solutions provided by the embodiments of the present disclosure have at least some or all of the following advantages: In nonlinear variable frequency scanning (such as pseudo-random scanning or custom variable frequency scanning), the instantaneous frequency of the scanning signal of the controllable vibrator is calculated based on the nonlinear modulation characteristics of the signal and the instantaneous frequency propagation law; according to the above instantaneous frequency, the harmonic instantaneous frequency curves corresponding to different sliding times and different orders of harmonic signals are generated; according to the intersection relationship between the above harmonic instantaneous frequency curves and the instantaneous frequency curves of the scanning signal corresponding to different slip times, the harmonic influence range of the controllable vibrator is determined; the harmonic influence range can be determined conveniently and efficiently and is applicable to nonlinear scanning signals or situations where the scanning signal cannot be accurately expressed by analytical expressions. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and, together with the description, serve to explain the principles of the present disclosure.

[0017] In order to more clearly illustrate the embodiments of the present disclosure or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or related technical descriptions. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0018] Figure 1 The flowchart of the method for determining the harmonic influence range of a controllable source according to an embodiment of the present disclosure is schematically shown.

[0019] Figure 2 The figure schematically shows a scanning signal diagram of a controllable vibrator according to an embodiment of the present disclosure.

[0020] Figure 3 The instantaneous phase diagram of the scanning signal of the controllable vibrator according to an embodiment of the present disclosure is schematically shown.

[0021] Figure 4 The figure schematically shows an instantaneous frequency diagram of a scanning signal of a controllable vibrator according to an embodiment of the present disclosure.

[0022] Figure 5 FIG. 1 schematically shows an instantaneous frequency diagram of a second-order harmonic according to an embodiment of the present disclosure.

[0023] Figure 6 FIG. 1 schematically shows an instantaneous frequency diagram of a third-order harmonic according to an embodiment of the present disclosure.

[0024] Figure 7 The figure schematically shows a curve diagram of the instantaneous frequency of the scanning signal corresponding to a certain shift time according to an embodiment of the present disclosure.

[0025] Figure 8 The figure schematically shows an intersection diagram of a second-order harmonic instantaneous frequency curve and a family of scan signal instantaneous frequency curves according to an embodiment of the present disclosure.

[0026] Figure 9 The figure schematically shows an intersection diagram of a third-order harmonic instantaneous frequency curve and a family of scan signal instantaneous frequency curves according to an embodiment of the present disclosure.

[0027] Figure 10 The structural block diagram of the electronic device provided by the embodiment of the present disclosure is schematically shown. DETAILED DESCRIPTION

[0028] To make the objectives, technical solutions, and advantages of the embodiments of the present disclosure more clear, the technical solutions in the embodiments of the present disclosure will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present disclosure, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present disclosure without making any creative efforts shall fall within the scope of protection of the present disclosure.

[0029] A first exemplary embodiment of the present disclosure provides a method for determining a vibroseis harmonic influence range, which can be applied to electronic devices with computing capabilities.

[0030] Figure 1 The flowchart of the method for determining the harmonic influence range of a controllable source according to an embodiment of the present disclosure is schematically shown.

[0031] Reference Figure 1 As shown, the method for determining the harmonic influence range of a vibroseis provided in an embodiment of the present disclosure includes the following steps: S110, S120, S130 and S140.

[0032] In step S110, a scanning signal of the vibrator is acquired, and the instantaneous frequency of the scanning signal is calculated.

[0033] In the embodiments of the present disclosure, the sweep signal of the controllable vibrator includes one of the following types: a nonlinear sweep signal, or a sweep signal for which no analytical expression is available. The method of this embodiment can be applied to these types of sweep signals to address the issue of the inability to adapt empirical formulas in related technologies. It should be noted that the method provided in the embodiments of the present disclosure is a general method and can also be applied to scenarios involving linear sweep signals or other types of sweep signals.

[0034] In some embodiments, in step S110 , calculating the instantaneous frequency of the scanning signal includes: obtaining the instantaneous frequency corresponding to the baseband signal from the scanning signal based on an analytical signal method or a time-frequency analysis method as the instantaneous frequency of the scanning signal.

[0035] For example, based on the analytical signal method, for the scanning signal , based on the Hilbert transform, the analytical signal is obtained The analytical signal satisfies the following expression: , (1) in, () represents the Hilbert transform function, j is the imaginary part.

[0036] The analytical signal The phase derivative is the instantaneous frequency corresponding to the fundamental frequency signal. The expression of instantaneous frequency is as follows: , (2) in, Indicates the instantaneous frequency of the sweep signal; Represents the analytic signal of the complex form The corresponding main value of the argument (corresponding to the instantaneous phase of the scanning signal , the unit is rad (radian)); express About time t The derivative of .

[0037] Video analysis methods, for example, can use energy ridge tracking in time-frequency distributions. By analyzing the signal's time-frequency representation (such as a time-frequency plot generated by a short-time Fourier transform or wavelet transform), the path along which characteristic frequencies change over time (i.e., the instantaneous frequency trajectory) can be determined. This method is primarily applied to non-stationary signals (such as frequency-modulated signals and speech signals) and can reveal the signal's dynamic characteristics in both time and frequency.

[0038] Figure 2 The figure schematically shows a scanning signal diagram of a controllable vibrator according to an embodiment of the present disclosure. Figure 3 The instantaneous phase diagram of the scanning signal of the controllable vibrator according to an embodiment of the present disclosure is schematically shown. Figure 4 The figure schematically shows an instantaneous frequency diagram of a scanning signal of a controllable vibrator according to an embodiment of the present disclosure.

[0039] For example, refer to Figure 2 As shown in one embodiment, the selected scanning signal frequency is a nonlinear scanning of 3Hz~84Hz, wherein the scanning length is 16s and the sampling interval is 2ms. Based on the analytical signal method, the instantaneous phase is solved for the scanning signal. , the differential of the instantaneous phase is the instantaneous frequency. Figure 3 As shown, the curve of instantaneous phase changing with time is shown; Figure 4 As shown in FIG, by differentiating the instantaneous frequency with respect to time (or describing it as a derivative), a curve of the instantaneous frequency changing with time is obtained.

[0040] In step S120, based on the instantaneous frequency, harmonic instantaneous frequency curves corresponding to different sliding times and harmonic signals of different orders are generated.

[0041] In some embodiments, in the above step S120, based on the above instantaneous frequency, a harmonic instantaneous frequency curve corresponding to different sliding times and different order harmonic signals is generated, including: modifying the expression of the above instantaneous frequency, multiplying the order coefficient before the instantaneous frequency, and introducing a sliding time term in the time parameter of the instantaneous frequency to represent the time shift existing in the harmonic propagation process, thereby obtaining a harmonic instantaneous frequency expression related to the sliding time and order coefficient; and drawing based on the above harmonic instantaneous frequency expression to obtain a harmonic instantaneous frequency curve.

[0042] In some embodiments, in the above step S120, based on the above instantaneous frequency, a harmonic instantaneous frequency curve corresponding to different sliding times and different order harmonic signals is generated, and the method further includes: receiving adjustment information for at least one of the above order coefficient or sliding time term; adjusting the above harmonic instantaneous frequency expression according to the above adjustment information, and drawing based on the adjusted expression to obtain an updated harmonic instantaneous frequency curve. That is, after obtaining the general formula of the harmonic instantaneous frequency curve, it supports adjusting the values ​​of the order coefficient, sliding time term, etc. In the case of adjustment, the subsequent steps S130 and S140 are also adapted to make adjustments and determine the harmonic influence range based on the intersection relationship of the adjusted curves.

[0043] Based on the principle that harmonics are integer multiple frequency distortions of the fundamental frequency signal and that their instantaneous frequency changes with time and scales synchronously with the fundamental frequency, the instantaneous frequency of the harmonics is preliminarily constructed by multiplying the instantaneous frequency by the order. At the same time, considering that there may be time shifts in harmonic propagation in actual seismic records (such as stratum filtering effects), it is necessary to introduce a sliding time window for correction. After correction, the instantaneous frequency expression of the harmonics is obtained.

[0044] In the specific implementation process, the above transformation process can be divided into two stages or completed in one go. Here is an example of the staged transformation process. First, the instantaneous frequency is multiplied by the order coefficient to meet the scaling principle. It is expressed by the following expression: , (3-1) in, represents the instantaneous frequency of the harmonics obtained by preliminary construction; k Indicates the order corresponding to the harmonic, which is used here as the order coefficient.

[0045] Afterwards, a sliding time window is introduced for correction, and the instantaneous frequency expression of the harmonic related to the sliding time and order coefficient is obtained: , (3-2) in, It represents the instantaneous frequency expression of harmonics obtained after correction based on the sliding time window; Represents sliding time (can also be described as a sliding time window).

[0046] Figure 5 FIG. 1 schematically shows an instantaneous frequency diagram of a second-order harmonic according to an embodiment of the present disclosure. Figure 6 FIG. 1 schematically shows an instantaneous frequency diagram of a third-order harmonic according to an embodiment of the present disclosure.

[0047] Reference Figure 5 and Figure 6 It can be seen that by respectively substituting the specific order coefficients 2 and 3 into the harmonic instantaneous frequency expressions related to the sliding time and order coefficient, the instantaneous frequency diagrams of the second-order harmonic and the third-order harmonic corresponding to the above scanning signal can be obtained. Figure 4 、 Figure 5 and Figure 6 From the above, we can see that the instantaneous frequency range of the scanning signal is 3Hz~84Hz, the instantaneous frequency range of the second-order harmonic is 6Hz~168Hz, and the instantaneous frequency range of the third-order harmonic is 9Hz~252Hz.

[0048] In step S130 , instantaneous frequency curves of the scanning signal corresponding to different shift times are generated according to the instantaneous frequency.

[0049] Those skilled in the art will appreciate that the shift time represents the difference between the arrival time of a harmonic signal (such as the second-order and third-order harmonics) and the theoretical distortion-free propagation time.

[0050] Since the instantaneous frequency corresponding to the scanning signal is obtained in the above step S110, the expression corresponding to the instantaneous frequency curve of the scanning signal corresponding to different shift times can be generated according to the instantaneous frequency: , (4) in, Indicates the time of slip The corresponding instantaneous frequency curve of the scanning signal; Indicates the shift time, with a specific range of -T~T, where T is the scan length of the scan signal.

[0051] Figure 7 The figure schematically shows a curve diagram of the instantaneous frequency of the scanning signal corresponding to a certain shift time according to an embodiment of the present disclosure.

[0052] Reference Figure 7 As shown, the instantaneous frequency curve of the scanning signal corresponding to the shift time range of -16s~16s is illustrated.

[0053] In step S140, the harmonic influence range of the vibrator is determined based on the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve.

[0054] In some embodiments, in the above step S140, determining the harmonic influence range of the vibrator according to the intersection relationship between the above harmonic instantaneous frequency curve and the above scanning signal instantaneous frequency curve includes: In response to the intersection relationship indicating that the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve have an intersection position, the intersection start time is determined according to the intersection position. and the intersection end time ; According to the above intersection starting time , intersection end time The sliding time corresponding to the above harmonic instantaneous frequency curve , determine the corresponding order (here shown as k The interference starting time of the harmonic signal (order) and the end time of the interference ; For example, the above intersection starting time Sliding time corresponding to the above harmonic instantaneous frequency curve Perform sum calculation to obtain the interference start time of the corresponding order harmonic signal interference ; The above intersection end time Sliding time corresponding to the above harmonic instantaneous frequency curve Perform sum calculation to obtain the end time of interference of the corresponding order harmonic signal ; The above interference starting time and the end time of the above interference The corresponding time period is determined as the above-mentioned vibrator at the corresponding order k and the corresponding sliding time The harmonic influence range is as follows.

[0055] That satisfies the following expression: , (5) . (6) Figure 8 The figure schematically shows an intersection diagram of a second-order harmonic instantaneous frequency curve and a family of scan signal instantaneous frequency curves according to an embodiment of the present disclosure. Figure 9 The diagram schematically shows the intersection of the instantaneous frequency curve of the third-order harmonic and the instantaneous frequency curve family of the scanning signal according to an embodiment of the present disclosure. Figure 8 、 Figure 9 In the figure, red lines are used to indicate the instantaneous frequency curves of the second-order harmonic and the third-order harmonic respectively.

[0056] Reference Figure 8As shown, from the intersection diagram of the second-order harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve family, it can be concluded that the intersection starting time of the second-order harmonic is -13.04s, the intersection end time is -0.12s. If the sliding time If it is 10s, the corresponding interference start time can be calculated by formula (5) and formula (6) respectively -3.04s, corresponding to the end time of interference If the value is 9.88s, it is determined that the corresponding period of harmonic interference of the controllable source under the second-order harmonic and the corresponding sliding time of 10s is: -3.04s~9.88s.

[0057] Reference Figure 9 As shown, from the intersection of the third-order harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve family, it can be concluded that the intersection starting time of the second-order harmonic is -14.4s, the intersection end time =-0.26s, and the interference start time corresponding to the sliding time of 10s is calculated by formulas (5) and (6) respectively -4.4s, corresponding to the end time of interference For example, if the third-order harmonic and a corresponding sliding time of 10 seconds are used, the harmonic interference period for the vibrator is -4.4 seconds to 9.74 seconds. The interference range for higher-order harmonics can be adjusted using the above method. For other sliding times, the same approach can be used to adjust the corresponding sliding time.

[0058] In some embodiments, the step S140 further includes a negative branch: in response to the intersection relationship indicating that the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve do not intersect, determining that the vibrator is not affected by harmonic interference within the corresponding shift time.

[0059] In an embodiment including steps S110 to S140, in a nonlinear variable frequency scan (such as a pseudo-random scan or a custom variable frequency scan), the instantaneous frequency of the scanning signal of the controllable vibrator is calculated based on the nonlinear modulation characteristics of the signal and the instantaneous frequency propagation law; according to the above instantaneous frequency, a harmonic instantaneous frequency curve corresponding to different sliding times and different-order harmonic signals is generated; according to the intersection relationship between the above harmonic instantaneous frequency curve and the instantaneous frequency curve of the scanning signal corresponding to different offset times, the harmonic influence range of the controllable vibrator is determined; the harmonic influence range can be determined conveniently and efficiently and is applicable to situations where nonlinear scanning signals or scanning signals cannot be accurately expressed by analytical expressions.

[0060] A second exemplary embodiment of the present disclosure provides a device for determining a vibroseis harmonic influence range.

[0061] The above-mentioned determination device includes: a signal acquisition module, an instantaneous frequency calculation module, a first curve generation module, a second curve generation module and an influence range determination module.

[0062] The signal acquisition module is used to acquire a scanning signal of a controllable vibrator. In some embodiments, the scanning signal includes one of the following types: a nonlinear scanning signal, and a scanning signal without an analytical expression.

[0063] The instantaneous frequency calculation module is used to calculate the instantaneous frequency of the scanning signal. In some embodiments, calculating the instantaneous frequency of the scanning signal includes: obtaining the instantaneous frequency corresponding to the baseband signal from the scanning signal based on an analytical signal method or a time-frequency analysis method, and using the instantaneous frequency as the instantaneous frequency of the scanning signal.

[0064] The first curve generating module is used to generate harmonic instantaneous frequency curves corresponding to different sliding times and different orders of harmonic signals according to the instantaneous frequency.

[0065] In some embodiments, the first curve generation module includes: an expression transformation submodule and a drawing submodule.

[0066] The above expression transformation submodule is used to transform the above instantaneous frequency expression by multiplying the order coefficient before the instantaneous frequency and introducing a sliding time term into the time parameter of the instantaneous frequency to represent the time shift in the harmonic propagation process, thereby obtaining the harmonic instantaneous frequency expression related to the sliding time and order coefficient.

[0067] The drawing submodule is used to draw a graph based on the harmonic instantaneous frequency expression to obtain a harmonic instantaneous frequency curve.

[0068] In some embodiments, the first curve generating module further includes: an updating submodule.

[0069] The updating submodule is configured to receive adjustment information for at least one of the order coefficient or the sliding time term; and adjust the harmonic instantaneous frequency expression according to the adjustment information.

[0070] The drawing submodule is further used to draw a graph based on the adjusted expression to obtain an updated harmonic instantaneous frequency curve.

[0071] The second curve generating module is used to generate instantaneous frequency curves of the scanning signal corresponding to different shift times according to the instantaneous frequency.

[0072] The above-mentioned influence range determination module is used to determine the harmonic influence range of the controllable vibrator according to the intersection relationship between the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve.

[0073] In some embodiments, the above-mentioned influence range determination module includes: an intersection start and end time determination submodule, an interference start and end time determination submodule, and an influence range generation submodule.

[0074] The intersection start and end time determination submodule is used to determine the intersection start time and the intersection end time according to the intersection position when the intersection relationship indicates that the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve have an intersection position.

[0075] The above-mentioned interference start and end time determination submodule is used to determine the interference start time and interference end time of the corresponding order harmonic signal based on the above-mentioned intersection start time, intersection end time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve. In some embodiments, the above-mentioned interference start and end time determination submodule is specifically used to add the above-mentioned intersection start time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve to obtain the interference start time of the corresponding order harmonic signal; add the above-mentioned intersection end time and the sliding time corresponding to the above-mentioned harmonic instantaneous frequency curve to obtain the interference end time of the corresponding order harmonic signal.

[0076] The above-mentioned influence range generating submodule is used to determine the time period corresponding to the above-mentioned interference start time and the above-mentioned interference end time as the harmonic influence range of the above-mentioned controllable source under the corresponding order and corresponding sliding time.

[0077] In some embodiments, the above-mentioned influence range generation submodule is also used to determine that the above-mentioned controllable seismic source has no harmonic interference influence within the corresponding offset time when the above-mentioned intersection relationship indicates that the above-mentioned harmonic instantaneous frequency curve and the above-mentioned scanning signal instantaneous frequency curve do not have an intersection position.

[0078] The device for determining the harmonic influence range of a controllable vibrator of this embodiment can determine the harmonic influence range conveniently and efficiently and is applicable to situations where nonlinear scanning signals or scanning signals cannot be accurately expressed by analytical expressions.

[0079] For more details of this embodiment, please refer to the relevant description of the first embodiment, which will not be repeated here.

[0080] Any number of the functional modules included in the determination device of this embodiment can be combined into a single module, or any one of these modules can be split into multiple modules. Alternatively, at least part of the functionality of one or more of these modules can be combined with at least part of the functionality of other modules and implemented in a single module. At least one of the functional modules included in the determination device of this embodiment can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on a chip, a system on a substrate, a system on a package, an application-specific integrated circuit (ASIC), or can be implemented in hardware or firmware through any other reasonable method of circuit integration or packaging, or can be implemented in any one of the three implementation methods of software, hardware, and firmware, or any appropriate combination of any of these. Alternatively, at least one of the functional modules included in the determination device of this embodiment can be at least partially implemented as a computer program module, which, when executed, can perform the corresponding function.

[0081] A third exemplary embodiment of the present disclosure provides an electronic device.

[0082] Figure 10 The structural block diagram of the electronic device provided by the embodiment of the present disclosure is schematically shown.

[0083] Reference Figure 10 As shown, the electronic device 1000 provided by the embodiment of the present disclosure includes a processor 1001, a communication interface 1002, a memory 1003 and a communication bus 1004, wherein the processor 1001, the communication interface 1002 and the memory 1003 communicate with each other through the communication bus 1004; the memory 1003 is used to store computer programs; the processor 1001 is used to implement the method for determining the harmonic influence range of the controllable source as described above when executing the program stored in the memory.

[0084] A fourth exemplary embodiment of the present disclosure further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for determining the harmonic influence range of a vibroseis as described above.

[0085] The computer-readable storage medium may be included in the device or apparatus described in the above embodiments, or may exist independently without being incorporated into the device or apparatus. The computer-readable storage medium carries one or more programs, and when the one or more programs are executed, the method according to the embodiments of the present disclosure is implemented.

[0086] According to embodiments of the present disclosure, a computer-readable storage medium may be a non-volatile computer-readable storage medium, such as, but not limited to, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present disclosure, a computer-readable storage medium may be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0087] In some embodiments, a computer program product is further provided, wherein a computer program is stored in the product. When the computer program is executed by a processor, the method for determining the harmonic influence range of a controllable vibrator as described above is implemented.

[0088] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or device comprising the element.

[0089] The foregoing description is intended only to provide specific embodiments of the present disclosure, intended to enable those skilled in the art to understand and implement the present disclosure. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present disclosure. Therefore, the present disclosure is not to be limited to the embodiments shown herein, but is to be construed in the broadest manner consistent with the principles and novel features claimed herein.

Claims

1. A method for determining the influence range of vibroseis harmonics, characterized in that: include: Acquiring a scanning signal of a controllable vibrator and calculating an instantaneous frequency of the scanning signal; Generating harmonic instantaneous frequency curves corresponding to harmonic signals of different sliding times and different orders according to the instantaneous frequency; Generating instantaneous frequency curves of the scanning signal corresponding to different shift times according to the instantaneous frequency; The harmonic influence range of the controllable vibrator is determined according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve.

2. The determination method according to claim 1, characterized in that Determining the harmonic influence range of the vibrator according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve includes: In response to the intersection relationship indicating that the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve have an intersection position, determining an intersection start time and an intersection end time according to the intersection position; Determine the interference start time and interference end time of the corresponding order harmonic signal according to the intersection start time, the intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve; The time period corresponding to the interference start time and the interference end time is determined as the harmonic influence range of the controllable source under the corresponding order and corresponding sliding time.

3. The determination method according to claim 2, characterized in that: Determining the interference start time and the interference end time of the corresponding order harmonic signal according to the intersection start time, the intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve, including: The intersection start time and the sliding time corresponding to the harmonic instantaneous frequency curve are added and calculated to obtain the interference start time when the interference of the corresponding order harmonic signal occurs; The intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve are added together to calculate and obtain the interference end time when the interference occurs in the harmonic signal of the corresponding order.

4. The determination method according to claim 2, characterized in that: Determining the harmonic influence range of the vibrator according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, further comprising: In response to the intersection relationship indicating that there is no intersection between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, it is determined that the vibrator is not affected by harmonic interference within the corresponding shift time.

5. The determination method according to claim 1, characterized in that: Calculating the instantaneous frequency of the scanning signal includes: Based on an analytical signal method or a time-frequency analysis method, an instantaneous frequency corresponding to the baseband signal is obtained from the scanning signal as the instantaneous frequency of the scanning signal.

6. The determination method according to any one of claims 1 to 5, characterized in that: The scanning signal includes one of the following types: a nonlinear scanning signal, a scanning signal without an analytical expression; According to the instantaneous frequency, generating harmonic instantaneous frequency curves corresponding to different sliding times and different order harmonic signals, including: The instantaneous frequency expression is modified by multiplying the instantaneous frequency by the order coefficient and introducing a sliding time term into the time parameter of the instantaneous frequency to represent the time shift in the harmonic propagation process, thereby obtaining an expression of the harmonic instantaneous frequency related to the sliding time and the order coefficient; Drawing is performed based on the harmonic instantaneous frequency expression to obtain a harmonic instantaneous frequency curve.

7. The determination method according to claim 6, characterized in that: Generating harmonic instantaneous frequency curves corresponding to harmonic signals of different sliding times and different orders according to the instantaneous frequency, further comprising: receiving adjustment information for at least one of the order coefficient or the sliding time term; The harmonic instantaneous frequency expression is adjusted according to the adjustment information, and a plot is performed based on the adjusted expression to obtain an updated harmonic instantaneous frequency curve.

8. A device for determining the influence range of vibroseis harmonics, characterized in that: include: A signal acquisition module, used for acquiring a scanning signal of a controllable vibrator; An instantaneous frequency calculation module, used to calculate the instantaneous frequency of the scanning signal; A first curve generating module is used to generate harmonic instantaneous frequency curves corresponding to harmonic signals of different sliding times and different orders according to the instantaneous frequency; A second curve generating module is used to generate instantaneous frequency curves of the scanning signal corresponding to different stagger times according to the instantaneous frequency; The influence range determination module is used to determine the harmonic influence range of the controllable vibrator according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve.

9. An electronic device, characterized in that: The processor, the communication interface, the memory and the communication bus are connected to each other via the communication bus. Memory for storing computer programs; A processor, configured to implement the method according to any one of claims 1 to 7 when executing a program stored in a memory.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images