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

Abstract

The present disclosure relates to a method, device, equipment and medium for determining the harmonic influence range of a controllable source, and relates to the technical field of seismic acquisition data analysis. The method comprises: obtaining a sweep signal of the controllable source, and calculating the instantaneous frequency of the sweep signal; generating a harmonic instantaneous frequency curve corresponding to different slip times and different order harmonic signals according to the instantaneous frequency; generating a sweep signal instantaneous frequency curve corresponding to different slip times according to the instantaneous frequency; and determining the harmonic influence range of the controllable source according to the intersection relationship between the harmonic instantaneous frequency curve and the sweep signal instantaneous frequency curve. The above scheme can conveniently and efficiently determine the harmonic influence range and is suitable for the case that the sweep signal is nonlinear or cannot be accurately expressed by an analytical expression.

Description

TECHNICAL FIELD

[0001] The present disclosure relates to the technical field of seismic acquisition data analysis, and particularly relates to a method and device for determining a harmonic influence range of a vibrator, equipment and a medium. BACKGROUND

[0002] Due to the controllable frequency energy, safety and environmental protection, low cost, high construction efficiency and many other advantages of the vibrator compared with the explosive source, the vibrator has become the main excitation source for land seismic exploration. The traditional scanning method of the vibrator is that each vibrator must wait for the scanning of the previous vibrator to be completely finished (including vibration time + monitoring time, i.e. after receiving the echo) before starting, which is low in efficiency. In order to improve the data acquisition efficiency, a sliding scanning technology is proposed during the data acquisition of the vibrator, which improves the acquisition efficiency by overlapping excitation (i.e. starting the scanning of the next vibrator when the scanning of the previous vibrator has not been completed). However, the sliding scanning inevitably brings some interference while improving the efficiency, and the harmonic interference is one of the most common interferences in the sliding scanning. Therefore, it is crucial to determine the influence range of the harmonic interference, which helps to optimize the construction parameters for harmonic interference suppression.

[0003] In the process of implementing the present disclosure, the inventors have found that there are at least the following technical problems in the related art: There are some empirical formulas to calculate the harmonic influence range of linear scanning, but the above-mentioned empirical formulas are only applicable to calculating the harmonic influence range of linear scanning, and when the scanning signal is not a linear scanning signal or the scanning signal cannot be accurately expressed by an analytical expression, the influence range of the harmonic interference cannot be determined by the above-mentioned empirical formulas. SUMMARY

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

[0005] In a first aspect, embodiments of the present disclosure provide a method for determining a harmonic influence range of a vibrator. The above-mentioned determination method comprises: obtaining a scanning signal of the vibrator, and calculating an instantaneous frequency of the scanning signal; generating a harmonic instantaneous frequency curve corresponding to different harmonic signals of different slip times and different orders according to the instantaneous frequency; generating a scanning signal instantaneous frequency curve corresponding to different slip times according to the instantaneous frequency; and 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.

[0006] In some embodiments, the harmonic influence range of the controllable seismic source is determined according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, including: in response to the intersection relationship indicating that there is an intersection position between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, determining an intersection start time and an intersection end time according to the intersection position; determining an interference start time and an interference end time at which the corresponding order harmonic signal appears to be interfered according to the intersection start time, the intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve; and determining a corresponding time period of the interference start time and the interference end time as the harmonic influence range of the controllable seismic source under the corresponding order and the corresponding sliding time.

[0007] In some embodiments, the interference start time and the interference end time at which the corresponding order harmonic signal appears to be interfered are determined according to the intersection start time, the intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve, including: adding the intersection start time and the sliding time corresponding to the harmonic instantaneous frequency curve to obtain the interference start time at which the corresponding order harmonic signal appears to be interfered; and adding the intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve to obtain the interference end time at which the corresponding order harmonic signal appears to be interfered.

[0008] In some embodiments, the harmonic influence range of the controllable seismic source is determined according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, further including: in response to the intersection relationship indicating that there is no intersection position between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, determining that there is no harmonic interference influence of the controllable seismic source within the corresponding slip time.

[0009] In some embodiments, the instantaneous frequency of the scanning signal is calculated, including: based on the analytical signal method or the time-frequency analysis method, obtaining the instantaneous frequency corresponding to the fundamental frequency signal from the scanning signal as the instantaneous frequency of the scanning signal.

[0010] In some embodiments, 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, the harmonic instantaneous frequency curve corresponding to the harmonic signal of different sliding time and different order is generated, including: transforming the expression of the instantaneous frequency, multiplying the order coefficient before the instantaneous frequency, and introducing the sliding time item in the time parameter of the instantaneous frequency to represent the time shift existing in the harmonic propagation process, to obtain the harmonic instantaneous frequency expression related to the sliding time and the order coefficient; and plotting based on the harmonic instantaneous frequency expression to obtain the harmonic instantaneous frequency curve.

[0011] In some embodiments, the harmonic instantaneous frequency curve corresponding to different sliding time and different order harmonic signals is generated according to the instantaneous frequency, and the method further comprises: receiving adjustment information for at least one of the order coefficient or the sliding time term; adjusting the harmonic instantaneous frequency expression according to the adjustment information, and drawing the updated harmonic instantaneous frequency curve based on the adjusted expression.

[0012] In a second aspect, the embodiments of the present disclosure provide a device for determining a harmonic influence range of a vibroseis. The device comprises 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 signal acquisition module is configured to acquire a scanning signal of the vibroseis. The instantaneous frequency calculation module is configured to calculate an instantaneous frequency of the scanning signal. The first curve generation module is configured to generate a harmonic instantaneous frequency curve corresponding to different sliding time and different order harmonic signals according to the instantaneous frequency. The second curve generation module is configured to generate a scanning signal instantaneous frequency curve corresponding to different slip times according to the instantaneous frequency. The influence range determination module is configured to determine 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.

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

[0014] In a fourth aspect, the embodiments of the present disclosure provide a computer readable storage medium. The computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the method for determining the harmonic influence range of the vibroseis.

[0015] The above technical solutions provided by the embodiments of the present disclosure have at least some or all of the following advantages:

[0016] In the nonlinear frequency scanning (such as pseudo-random scanning or custom frequency scanning), based on the nonlinear modulation characteristics of the signal and the instantaneous frequency propagation law, the instantaneous frequency of the scanning signal of the controllable source is calculated; according to the above instantaneous frequency, the harmonic instantaneous frequency curves corresponding to different slip times and different order harmonic signals are generated; according to the intersection relationship of the above harmonic instantaneous frequency curves and the scanning signal instantaneous frequency curves corresponding to different slip times, the harmonic influence range of the controllable source is determined; the harmonic influence range can be determined conveniently and efficiently and is suitable for nonlinear scanning signals or scanning signals that cannot be accurately expressed by analytical expressions. BRIEF DESCRIPTION OF DRAWINGS

[0017] The accompanying drawings, which are incorporated herein and form part of the specification, illustrate embodiments consistent with the present disclosure and, together with the description, further serve to explain the principles of the disclosure.

[0018] In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure or the prior art, the accompanying drawings needed to be used in the embodiments or related description will be briefly introduced. Obviously, those skilled in the art can obtain other drawings according to these drawings without any creative effort.

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

[0020] Figure 2 A scanning signal schematic diagram of a controllable source according to an embodiment of the present disclosure is schematically shown.

[0021] Figure 3 An instantaneous phase diagram of a scanning signal of a controllable source according to an embodiment of the present disclosure is schematically shown.

[0022] Figure 4 An instantaneous frequency diagram of a scanning signal of a controllable source according to an embodiment of the present disclosure is schematically shown.

[0023] Figure 5 An instantaneous frequency diagram of a second-order harmonic according to an embodiment of the present disclosure is schematically shown.

[0024] Figure 6 An instantaneous frequency diagram of a third-order harmonic according to an embodiment of the present disclosure is schematically shown.

[0025] Figure 7 A scanning signal instantaneous frequency curve diagram corresponding to a certain slip time according to an embodiment of the present disclosure is schematically shown.

[0026] Figure 8A plot of a second-order harmonic instantaneous frequency curve intersecting with a family of scan signal instantaneous frequency curves is schematically shown according to an embodiment of the present disclosure.

[0027] Figure 9 A plot of a third-order harmonic instantaneous frequency curve intersecting with a family of scan signal instantaneous frequency curves is schematically shown according to an embodiment of the present disclosure.

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

[0029] To make the objectives, technical solutions, and superiorities of the embodiments of the present disclosure clearer, the technical solutions in the embodiments of the present disclosure will be described below in connection with the drawings in the embodiments of the present disclosure. Obviously, the described embodiments are only some of the embodiments of the present disclosure, rather than all the embodiments. Based on the embodiments in the present disclosure, any other embodiments obtained by those of ordinary skill in the art without creative efforts are within the scope of the present disclosure.

[0030] A first exemplary embodiment of the present disclosure provides a method for determining a harmonic influence range of a controllable seismic source. The determining method can be applied in an electronic device with computing capability.

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

[0032] REFERENCE Figure 1 As shown in the figure, the method for determining the harmonic influence range of the controllable seismic source provided by an embodiment of the present disclosure includes the following steps: S110, S120, S130, and S140.

[0033] In step S110, a scan signal of a controllable seismic source is acquired, and an instantaneous frequency of the scan signal is calculated.

[0034] In an embodiment of the present disclosure, the scan signal of the controllable seismic source includes one of the following types: a nonlinear scan signal, a scan signal without an analytical expression. The method of the present embodiment can be applied to the above-mentioned types of scan signals to solve the problem that the empirical formula in the related art cannot be adapted; it should be noted that the method provided by an embodiment of the present disclosure is a general method, which can also be applied to the scenario of a linear scan signal or other various types of scan signals.

[0035] In some embodiments, in step S110, the instantaneous frequency of the scan signal is calculated, including: based on an analytical signal method or a time-frequency analysis method, an instantaneous frequency corresponding to a fundamental frequency signal of the scan signal is acquired from the scan signal as the instantaneous frequency of the scan signal.

[0036] For example, based on the analytical signal method, for scanning signals... Analytical signals are obtained based on Hilbert transform. An analytic signal satisfies the following expression:

[0037] (1)

[0038] in, () represents the Hilbert transform function. j This is the imaginary part.

[0039] This analytical signal The phase derivative is the instantaneous frequency corresponding to the fundamental frequency signal. The expression for the instantaneous frequency is as follows:

[0040] (2)

[0041] in, Indicates the instantaneous frequency of the scanning signal; Analytic signals representing complex numbers The corresponding principal argument value (corresponding to the instantaneous phase of the scanning signal) (Unit: rad (radians)) express Regarding time t The derivative of .

[0042] Video analysis methods, such as energy ridge tracing based on time-frequency distributions, determine the path of characteristic frequencies over time (i.e., instantaneous frequency trajectories) by analyzing the time-frequency representation of the signal (e.g., time-frequency graphs generated by short-time Fourier transform or wavelet transform). This method is primarily applied to non-stationary signals (such as FM signals and speech signals) and can reveal the dynamic characteristics of the signal in time and frequency.

[0043] Figure 2 A schematic diagram of a scanned signal of a controllable seismic source according to an embodiment of the present disclosure is shown. Figure 3 An instantaneous phase diagram of a scanned signal from a controllable seismic source according to an embodiment of the present disclosure is schematically shown. Figure 4 An instantaneous frequency diagram of a scanned signal from a controllable seismic source according to an embodiment of the present disclosure is schematically shown.

[0044] For example, refer to Figure 2 As shown, in one embodiment, the selected scanning signal frequency is a nonlinear scan of 3Hz to 84Hz, where the scan length is 16s and the sampling interval is 2ms. Based on the analytical signal method, the instantaneous phase of the scanning signal is solved. The differential of the instantaneous phase is the instantaneous frequency. (Refer to...) Figure 3 The figure illustrates the curve of instantaneous phase change over time; refer toFigure 4 As shown, the curve of the instantaneous frequency over time is obtained by taking the derivative of the instantaneous frequency with respect to time (or described as a derivative).

[0045] In step S120, according to the above-mentioned instantaneous frequency, the harmonic instantaneous frequency curve corresponding to different sliding time and different order harmonic signals is generated.

[0046] In some embodiments, in step S120, according to the above-mentioned instantaneous frequency, the harmonic instantaneous frequency curve corresponding to different sliding time and different order harmonic signals is generated, including: modifying the expression of the above-mentioned instantaneous frequency, multiplying the order coefficient before the instantaneous frequency, and introducing the sliding time item in the time parameter of the instantaneous frequency to represent the time shift existing in the harmonic propagation process, to obtain the harmonic instantaneous frequency expression related to the sliding time and the order coefficient; based on the above-mentioned harmonic instantaneous frequency expression, the harmonic instantaneous frequency curve is obtained.

[0047] In some embodiments, in step S120, according to the above-mentioned instantaneous frequency, the harmonic instantaneous frequency curve corresponding to different sliding time and different order harmonic signals is generated, and further includes: receiving adjustment information for at least one of the order coefficient or the sliding time item; according to the adjustment information, the harmonic instantaneous frequency expression is adjusted, and based on the adjusted expression, the updated harmonic instantaneous frequency curve is obtained. That is, after obtaining the general formula of the harmonic instantaneous frequency curve, the values of the order coefficient, the sliding time item, etc. are adjusted. In the case of adjustment, the subsequent steps S130 and S140 are also adapted to adjust and determine the harmonic influence range based on the intersection relationship of the adjusted curve.

[0048] Based on the principle that the harmonic is an integer multiple frequency distortion of the fundamental signal, and the instantaneous frequency changes with time synchronously scaled with the fundamental frequency, by multiplying the order number before the instantaneous frequency, the harmonic instantaneous frequency is preliminarily constructed; at the same time, considering that in actual seismic records, there may be time shift in harmonic propagation (such as stratum filtering effect), it is necessary to correct it by introducing a sliding time window, and after correction, the harmonic instantaneous frequency expression is obtained.

[0049] In the specific implementation process, the above-mentioned modification process is divided into two stages or completed at one time. Here, the process of stage modification is shown, first multiplying the order coefficient before the instantaneous frequency, satisfying the scaling principle, and using the following expression:

[0050] , (3-1)

[0051] Wherein, represents the harmonic instantaneous frequency preliminarily constructed; k represents the order corresponding to the harmonic, which is used as the order coefficient here.

[0052] Then, the harmonic instantaneous frequency expression related to the sliding time and order coefficient is obtained by introducing the sliding time window for correction:

[0053] , (3-2)

[0054] wherein, represents the harmonic instantaneous frequency expression obtained after correction based on the sliding time window; represents the sliding time (which can also be described as the sliding time window).

[0055] Figure 5 The instantaneous frequency diagram of the second-order harmonic according to an embodiment of the present disclosure is schematically shown. Figure 6 The instantaneous frequency diagram of the third-order harmonic according to an embodiment of the present disclosure is schematically shown.

[0056] Referring to Figure 5 and Figure 6 , it can be seen that the specific order coefficients 2, 3 are respectively brought into the harmonic instantaneous frequency expression related to the sliding time and order coefficient, and the instantaneous frequency diagram of the second-order harmonic corresponding to the above scanning signal and the instantaneous frequency diagram of the third-order harmonic can be obtained. Comparing Figure 4 , Figure 5 and Figure 6 , it can be seen 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.

[0057] In step S130, according to the above instantaneous frequency, the scanning signal instantaneous frequency curve corresponding to different jitter times is generated.

[0058] As known by those skilled in the art, the jitter time represents the difference between the arrival time of the harmonic signal (such as the second-order and third-order harmonic) and the theoretical undistorted propagation time.

[0059] Since the instantaneous frequency corresponding to the scanning signal is obtained in the above step S110, the expression of the scanning signal instantaneous frequency curve corresponding to different jitter times can be generated according to the instantaneous frequency:

[0060] , (4)

[0061] wherein, represents the jitter time corresponding to the scanning signal instantaneous frequency curve; represents the jitter time, and the specific range is -T~T, T being the scanning length of the scanning signal.

[0062] Figure 7The diagram schematically illustrates the instantaneous frequency curve of the scan signal at a certain time of misalignment according to an embodiment of the present disclosure.

[0063] Reference Figure 7 As shown, the instantaneous frequency curve of the scanning signal is illustrated within the range of -16s to 16s.

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

[0065] In some embodiments, in step S140 above, determining the harmonic influence range of the controllable source based on the intersection relationship between the instantaneous harmonic frequency curve and the instantaneous frequency curve of the scanning signal includes:

[0066] In response to the aforementioned intersection relationship indicating that the instantaneous frequency curve of the harmonics and the instantaneous frequency curve of the scanning signal intersect at a certain point, the intersection start time is determined based on the intersection point. and the end of the intersection ;

[0067] Based on the above intersection start time Intersection End Time The sliding time corresponding to the instantaneous frequency curve of the above harmonics Determine the corresponding order (here illustrated as) k The onset time of interference from (order) harmonic signals. and the end of the interference For example, the above intersection start time The sliding time corresponding to the above harmonic instantaneous frequency curve By performing summation calculations, the interference onset time of the corresponding order harmonic signal can be obtained. The intersection time mentioned above ends The sliding time corresponding to the above harmonic instantaneous frequency curve By performing summation calculations, the time when the interference of the corresponding order harmonic signal ends can be obtained. ;

[0068] The above interference start time and the time when the above interference ends The corresponding time period is determined to be the controllable seismic source at the corresponding order k and the corresponding sliding time. The range of harmonic influence.

[0069] That is, satisfying the following expression:

[0070] (5)

[0071] (6)

[0072] Figure 8 Fig. 6 schematically shows the intersection of the second-order harmonic instantaneous frequency curve and the family of scanning signal instantaneous frequency curves according to an embodiment of the present disclosure. Figure 9 Fig. 7 schematically shows the intersection of the third-order harmonic instantaneous frequency curve and the family of scanning signal instantaneous frequency curves according to an embodiment of the present disclosure. Figure 8 Figure 9 In Fig. 6, the second-order harmonic instantaneous frequency curve and the third-order harmonic instantaneous frequency curve are respectively shown by red lines.

[0073] Referring to Fig. 6, it can be seen from the intersection of the second-order harmonic instantaneous frequency curve and the family of scanning signal instantaneous frequency curves that the intersection start time of the second-order harmonic is -13.04 s, and the intersection end time is -0.12 s. Figure 8 If the sliding time T is 10 s, the corresponding interference start time T0 is -3.04 s, and the corresponding interference end time T1 is 9.88 s, which can be calculated by formula (5) and formula (6) respectively.

[0074] Referring to Fig. 7, it can be seen from the intersection of the third-order harmonic instantaneous frequency curve and the family of scanning signal instantaneous frequency curves that the intersection start time of the second-order harmonic is -14.4 s, and the intersection end time is -0.26 s. Figure 9 If the sliding time T is 10 s, the corresponding interference start time T0 is -4.4 s, and the corresponding interference end time T1 is 9.74 s, which can be calculated by formula (5) and formula (6) respectively.

[0075] In some embodiments, the step S140 further includes a negative branch: in response to the above intersection relationship indicating that there is no intersection position between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, it is determined that there is no harmonic interference effect on the controllable source within the corresponding slip time.

[0076] ​​​​​​​​​​In the embodiment comprising steps S110-S140, in a nonlinear frequency sweep (such as a pseudo-random sweep or a self-defined frequency sweep), based on the nonlinear modulation characteristics of the signal and the instantaneous frequency propagation law, the instantaneous frequency of the sweep signal of the controllable source is calculated; based on the above instantaneous frequency, the harmonic instantaneous frequency curves corresponding to different slip times and different order harmonic signals are generated; based on the intersection relationship between the above harmonic instantaneous frequency curves and the sweep signal instantaneous frequency curves corresponding to different slip times, the harmonic influence range of the controllable source is determined; the harmonic influence range can be conveniently and efficiently determined and is suitable for nonlinear sweep signals or sweep signals that cannot be accurately expressed by analytical expressions.

[0077] A second exemplary embodiment of the present disclosure provides a device for determining the harmonic influence range of a controllable source.

[0078] The above-mentioned determination device comprises 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.

[0079] The above-mentioned signal acquisition module is used to acquire the sweep signal of the controllable source. In some embodiments, the sweep signal comprises one of the following types: a nonlinear sweep signal, a sweep signal without an analytical expression.

[0080] The above-mentioned instantaneous frequency calculation module is used to calculate the instantaneous frequency of the above-mentioned sweep signal. In some embodiments, the instantaneous frequency of the above-mentioned sweep signal is calculated, comprising: based on the analytical signal method or the time-frequency analysis method, the instantaneous frequency corresponding to the fundamental frequency signal is acquired from the above-mentioned sweep signal as the instantaneous frequency of the above-mentioned sweep signal.

[0081] The above-mentioned first curve generation module is used to generate the harmonic instantaneous frequency curves corresponding to different slip times and different order harmonic signals based on the above-mentioned instantaneous frequency.

[0082] In some embodiments, the above-mentioned first curve generation module comprises an expression modification submodule and a drawing submodule.

[0083] The above-mentioned expression modification submodule is used to modify the expression of the above-mentioned instantaneous frequency, multiply the order coefficient before the instantaneous frequency, and introduce the slip time item in the time parameter of the instantaneous frequency to represent the time shift existing in the harmonic propagation process, to obtain the harmonic instantaneous frequency expression related to the slip time and the order coefficient.

[0084] The above-mentioned drawing submodule is used to draw based on the above-mentioned harmonic instantaneous frequency expression to obtain the harmonic instantaneous frequency curve.

[0085] In some embodiments, the above-mentioned first curve generation module further comprises an updating submodule.

[0086] 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.

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

[0088] The second curve generating module is configured to generate a scanning signal instantaneous frequency curve corresponding to different slip times according to the instantaneous frequency.

[0089] The influence range determining module is configured to determine a harmonic influence range of the controllable source according to an intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve.

[0090] In some embodiments, the influence range determining module includes an intersection start and end time determining submodule, an interference start and end time determining submodule, and an influence range generating submodule.

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

[0092] The interference start and end time determining submodule is configured to determine an interference start time and an interference end time of interference of a harmonic signal of a corresponding order according to the intersection start time, the intersection end time, and a sliding time corresponding to the harmonic instantaneous frequency curve. In some embodiments, the interference start and end time determining submodule is specifically configured to add the intersection start time and the sliding time corresponding to the harmonic instantaneous frequency curve to obtain the interference start time of the interference of the harmonic signal of the corresponding order, and add the intersection end time and the sliding time corresponding to the harmonic instantaneous frequency curve to obtain the interference end time of the interference of the harmonic signal of the corresponding order.

[0093] The influence range generating submodule is configured to determine a period corresponding to the interference start time and the interference end time as a harmonic influence range of the controllable source under the corresponding order and the corresponding sliding time.

[0094] In some embodiments, the influence range generating submodule is further configured to determine that there is no harmonic interference influence of the controllable source under the corresponding slip time when the intersection relationship indicates that the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve do not have the intersection position.

[0095] The device for determining the harmonic influence range of the controllable source in the embodiment can conveniently and efficiently determine the harmonic influence range and is applicable to the case that the nonlinear scanning signal or the scanning signal cannot be accurately expressed by an analytical expression.

[0096] More details of the embodiment can be referred to the related description of the first embodiment, which will not be repeated here.

[0097] Any of the functional modules included in the determining device of the embodiment can be combined in one module, or any of the modules can be split into multiple modules. Alternatively, at least part of the function of one or more of the modules can be combined with at least part of the function of other modules, and implemented in one module. At least one of the functional modules included in the determining device of the 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 chip, a system on board, a system on package, an application specific integrated circuit (ASIC), or any other reasonable way of integrating or packaging a circuit, etc. hardware or firmware, or any one of software, hardware and firmware or any appropriate combination of several of them. Alternatively, at least one of the functional modules included in the determining device of the embodiment can be at least partially implemented as a computer program module which can perform the corresponding function when executed.

[0098] The third exemplary embodiment of the present disclosure provides an electronic device.

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

[0100] Referring to Figure 10 As shown in the figure, 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 complete the communication among each other through the communication bus 1004; the memory 1003 is used for storing a computer program; the processor 1001 is used for executing the program stored on the memory to realize the method for determining the harmonic influence range of the controllable source as described above.

[0101] The fourth exemplary embodiment of the present disclosure also provides a computer readable storage medium. The computer readable storage medium stores a computer program, and the computer program is executed by a processor to realize the method for determining the harmonic influence range of the controllable source as described above.

[0102] The computer readable storage medium can be included in the devices or apparatuses described in the above embodiments; it can also exist independently, and is not assembled into the devices or apparatuses. The computer readable storage medium described above bears one or more programs, which, when executed, implement the method according to the embodiments of the present disclosure.

[0103] According to the embodiments of the present disclosure, the computer readable storage medium can be a non-volatile computer readable storage medium, which can include, but is not limited to, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any appropriate combination thereof. In the present disclosure, the computer readable storage medium can be any tangible medium that bears, or stores, a program that can be used by, or in connection with, an instruction execution system, apparatus, or device.

[0104] In some embodiments, a computer program product is also provided, which stores a computer program, and the computer program, when executed by a processor, implements the method for determining the harmonic influence range of a controllable source as described above.

[0105] It should be noted that, in the present document, the relationship terms such as “first” and “second” are merely used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply that there is any such actual relationship or sequence between these entities or operations. Moreover, the terms “include”, “contain” or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device that includes a series of elements not only includes those elements, but also includes other elements not explicitly listed or inherent to such a process, method, article or device. Without more limitations, the element defined by the statement “including a…” does not exclude the presence of other identical elements in the process, method, article or device that includes the element.

[0106] The above description is merely one specific implementation of the present disclosure, which enables those skilled in the art to understand or implement the present disclosure. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present disclosure. Therefore, the present disclosure will not be limited to these embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features presented herein.

Claims

1. A method of determining a harmonic influence range of a controllable source, characterized in that, The method comprises the following steps: acquiring a sweep signal of a controllable source and calculating an instantaneous frequency of the sweep signal; generating a harmonic instantaneous frequency curve corresponding to different slip times and different order harmonic signals according to the instantaneous frequency; generating a sweep signal instantaneous frequency curve corresponding to different slip times according to the instantaneous frequency; determining a harmonic influence range of the controllable source according to an intersection relationship between the harmonic instantaneous frequency curve and the sweep signal instantaneous frequency curve; determining a harmonic influence range of the controllable source according to an intersection relationship between the harmonic instantaneous frequency curve and the sweep signal instantaneous frequency curve, comprising: in response to the intersection relationship indicating that there is an intersection position between the harmonic instantaneous frequency curve and the sweep signal instantaneous frequency curve, determining an intersection start time and an intersection end time according to the intersection position; determining an interference start time and an interference end time of the interference of the corresponding order harmonic signal according to the intersection start time, the intersection end time and the slip time corresponding to the harmonic instantaneous frequency curve; determining a corresponding period of the interference start time and the interference end time as the harmonic influence range of the controllable source under the corresponding order and the corresponding slip time; in response to the intersection relationship indicating that there is no intersection position between the harmonic instantaneous frequency curve and the sweep signal instantaneous frequency curve, determining that there is no harmonic interference influence of the controllable source within the corresponding slip time; determining an interference start time and an interference end time of the interference of the corresponding order harmonic signal according to the intersection start time, the intersection end time and the slip time corresponding to the harmonic instantaneous frequency curve, comprising: performing addition calculation on the intersection start time and the slip time corresponding to the harmonic instantaneous frequency curve to obtain the interference start time of the interference of the corresponding order harmonic signal; performing addition calculation on the intersection end time and the slip time corresponding to the harmonic instantaneous frequency curve to obtain the interference end time of the interference of the corresponding order harmonic signal.

2. The determination method according to claim 1, characterized in that, The method for calculating the instantaneous frequency of the sweep signal comprises the following steps: based on the analytical signal method or the time-frequency analysis method, acquiring the instantaneous frequency corresponding to the fundamental frequency signal from the sweep signal as the instantaneous frequency of the sweep signal.

3. The determination method according to claim 1 or 2, characterized in that, The sweep signal comprises one of the following types: a nonlinear sweep signal, a sweep signal without an analytical expression; The method for generating a harmonic instantaneous frequency curve corresponding to different slip times and different order harmonic signals according to the instantaneous frequency comprises the following steps: transforming the expression of the instantaneous frequency, multiplying the order coefficient before the instantaneous frequency, and introducing a slip time item in the time parameter of the instantaneous frequency to represent the time shift existing in the harmonic propagation process, to obtain a harmonic instantaneous frequency expression related to the slip time and the order coefficient; based on the harmonic instantaneous frequency expression, drawing to obtain the harmonic instantaneous frequency curve.

4. The determination method according to claim 3, characterized in that, The method for generating a harmonic instantaneous frequency curve corresponding to different slip times and different order harmonic signals according to the instantaneous frequency further comprises the following steps: receiving adjustment information for at least one of the order coefficient or the slip time item; According to the adjustment information, the harmonic instantaneous frequency expression is adjusted, and a drawing is made based on the adjusted expression to obtain an updated harmonic instantaneous frequency curve.

5. A device for determining the influence range of harmonics from a controllable seismic source, characterized in that, Comprise: The signal acquisition module is used for acquiring the scanning signal of the controllable source; The instantaneous frequency calculation module is used for calculating the instantaneous frequency of the scanning signal; The first curve generation module is used for generating the harmonic instantaneous frequency curve corresponding to the harmonic signal of different slip times and different orders according to the instantaneous frequency; The second curve generation module is used for generating the scanning signal instantaneous frequency curve corresponding to different slip times according to the instantaneous frequency; The influence range determination module is used for determining the harmonic influence range of the controllable source according to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve; According to the intersection relationship between the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve, the harmonic influence range of the controllable source is determined, comprising: 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 and the intersection end time are determined according to the intersection position; According to the intersection start time, the intersection end time and the slip time corresponding to the harmonic instantaneous frequency curve, the interference start time and the interference end time of the interference of the corresponding order harmonic signal are determined; The interference start time and the interference end time corresponding to the period are determined as the harmonic influence range of the controllable source under the corresponding order and the corresponding slip time; In response to the intersection relationship indicating that the harmonic instantaneous frequency curve and the scanning signal instantaneous frequency curve do not have an intersection position, it is determined that the controllable source does not have harmonic interference influence within the corresponding slip time; According to the intersection start time, the intersection end time and the slip time corresponding to the harmonic instantaneous frequency curve, the interference start time and the interference end time of the interference of the corresponding order harmonic signal are determined, comprising: The intersection start time and the slip time corresponding to the harmonic instantaneous frequency curve are added to obtain the interference start time of the interference of the corresponding order harmonic signal; The intersection end time and the slip time corresponding to the harmonic instantaneous frequency curve are added to obtain the interference end time of the interference of the corresponding order harmonic signal.

6. An electronic device, comprising: Comprise processor, communication interface, memory and communication bus, wherein the processor, the communication interface and the memory complete the communication among each other through the communication bus; The memory is used for storing the computer program; The processor is used for executing the program stored on the memory, and the method of any one of claims 1-4 is realized.

7. A computer-readable storage medium having stored thereon a computer program, characterized in that The computer program is executed by the processor to realize the method of any one of claims 1-4. The computer program is executed by the processor to realize the method of any one of claims 1-4.

Images