Method, device, electronic equipment and medium for vibroseis spectrum shaping

Abstract

This disclosure relates to a method, apparatus, electronic device, and medium for controlled-source seismic spectrum shaping. The method includes: determining a first frequency domain sequence of an explosive excitation spectrum and a second frequency domain sequence of a controlled-source excitation spectrum; the explosive excitation spectrum and the controlled-source excitation spectrum are acquired signal spectra obtained by excitation at the same exploration location; performing amplitude absolute value taking and amplitude normalization processing on the first frequency domain sequence and the second frequency domain sequence respectively to obtain corresponding first target frequency domain sequence and second target frequency domain sequence; determining an amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence; and shaping the amplitude of the second target frequency domain sequence based on the amplitude adjustment coefficient to obtain a second shaped frequency domain sequence. This method can improve the display signal-to-noise ratio of controlled-source data and reduce the display differences between controlled-source data and explosive data.

Description

Technical Field

[0001] This disclosure relates to the field of geophysical exploration technology, and in particular to a method, apparatus, electronic device and medium for controllable source spectrum shaping. Background Technology

[0002] Driven by cost, environmental, and safety concerns, as well as the continuous improvement of controlled-source technology, controlled-source excitation has gradually replaced explosive excitation as the primary source excitation method in seismic exploration. The replacement of traditional explosive-based collection methods with controlled-source collection has thus become a mainstream trend.

[0003] However, in realizing the concept disclosed herein, the inventors discovered at least the following technical problems in the related technology: when comparing and analyzing explosive data (data in this field can be understood as acquired data) and source data, the difference in the dominant frequency bands between controlled source data and explosive data leads to a visual difference in the signal-to-noise ratio of the two acquired data. For example, when comparing the post-stack cross-sections of controlled source data and explosive data, the video rate and wave group characteristics of explosive data and controlled source data are significantly different, resulting in a certain gap in the displayed imaging effect, causing some confusion for subsequent data interpreters. Summary of the Invention

[0004] To solve the above-mentioned technical problems, or at least partially solve them, embodiments of this disclosure provide a method, apparatus, electronic device, and medium for controllable source spectrum shaping.

[0005] In a first aspect, embodiments of this disclosure provide a method for controlled source spectrum shaping. The method includes: determining a first frequency domain sequence of an explosive excitation spectrum and a second frequency domain sequence of a controlled source excitation spectrum; the explosive excitation spectrum and the controlled source excitation spectrum are acquired signal spectra obtained by excitation at the same exploration location; performing absolute amplitude calculation and amplitude normalization on the first frequency domain sequence and the second frequency domain sequence respectively to obtain corresponding first target frequency domain sequence and second target frequency domain sequence; determining an amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence; and shaping the amplitude of the second target frequency domain sequence based on the amplitude adjustment coefficient to obtain a second shaped frequency domain sequence.

[0006] In some embodiments, determining the amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence includes: performing a quotient operation on the amplitudes at the same frequency in the first target frequency domain sequence and the second target frequency domain sequence to obtain the amplitude ratio at each frequency; using the amplitude ratio at each frequency as a proposed coefficient, performing a coefficient product operation on the amplitude at the corresponding frequency of the second target frequency domain sequence to obtain a proposed spectrum adjustment result; comparing and determining whether the proposed spectrum adjustment result is consistent with the bandwidth of the second target frequency domain sequence; if they are consistent throughout the entire frequency range, determining the proposed coefficient at each frequency in the entire frequency range as the amplitude adjustment coefficient at the corresponding frequency.

[0007] In some embodiments, determining the amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence further includes: adjusting the proposed coefficients within the target frequency domain when there is inconsistency in the target frequency segments, until the proposed spectrum adjustment result of the target frequency segment corresponding to the adjusted proposed coefficient is consistent with the bandwidth of the second target frequency domain sequence; and determining the proposed coefficients after the target frequency segment adjustment and the proposed coefficients in other frequency segments as the amplitude adjustment coefficients at the corresponding frequencies.

[0008] In some embodiments, if the bandwidth of the proposed spectrum adjustment result starting from the target frequency is greater than the bandwidth of the second target frequency domain sequence, then the frequency segment above the target frequency is determined as the target frequency segment.

[0009] In some embodiments, the amplitude of the second target frequency domain sequence is shaped according to the amplitude adjustment coefficient to obtain a second shaped frequency domain sequence, including: multiplying the amplitude adjustment coefficient and the amplitude of the second target frequency domain sequence at the corresponding frequency to obtain a second shaped frequency domain sequence.

[0010] In some embodiments, determining the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum includes: loading the first time domain sequence of the explosive excitation spectrum and the second time domain sequence of the controllable source excitation spectrum; performing Fourier transform processing on the first time domain sequence and the second time domain sequence to obtain the corresponding first frequency domain sequence and second frequency domain sequence.

[0011] In some embodiments, the method for shaping the spectrum of the controllable source further includes: performing an inverse Fourier transform on the second shaped frequency domain sequence to obtain the shaped time domain sequence corresponding to the controllable source.

[0012] Secondly, embodiments of this disclosure provide an apparatus for controlling the shaping of a seismic source spectrum. The apparatus includes: a frequency sequence determination module, a normalization processing module, a shaping coefficient determination module, and a shaping processing module. The frequency sequence determination module is used to determine a first frequency domain sequence of the explosive excitation spectrum and a second frequency domain sequence of the controlled seismic source excitation spectrum; the explosive excitation spectrum and the controlled seismic source excitation spectrum are acquired signal spectra obtained by excitation at the same exploration location. The normalization processing module is used to perform amplitude absolute value taking and amplitude normalization processing on the first frequency domain sequence and the second frequency domain sequence, respectively, to obtain corresponding first target frequency domain sequence and second target frequency domain sequence. The shaping coefficient determination module is used to determine an amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence. The shaping processing module is used to shape the amplitude of the second target frequency domain sequence based on the amplitude adjustment coefficient to obtain a second shaped frequency domain sequence.

[0013] Thirdly, embodiments of this disclosure provide an electronic device. The electronic device includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, communication interface, and memory communicate with each other via the communication bus; the memory stores computer programs; and the processor, when executing the program stored in the memory, implements the controllable source spectrum shaping method described above.

[0014] Fourthly, embodiments of this disclosure provide a computer-readable storage medium. The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method for controllable source spectrum shaping as described above.

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

[0016] By taking the absolute value and normalizing the amplitudes of the first and second frequency domain sequences of controlled seismic sources and explosives at the same exploration location, the amplitude energy levels between the two can be unified, facilitating comparison of relative amplitude differences. An amplitude adjustment coefficient is determined based on the amplitude ratio of the first and second target frequency domain sequences. This coefficient is used to achieve matched filtering of controlled seismic source data, adjusting the relative amplitudes at different frequencies to achieve shaping, improving the signal-to-noise ratio of the controlled seismic source data, and reducing the comparative display differences between controlled seismic source data and explosive data. This helps later structural interpreters to conduct geological structural analysis (e.g., stratigraphic tracing) by comparing and analyzing the shaped controlled seismic source data with the explosive data. Attached Figure Description

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

[0018] To more clearly illustrate the technical solutions in the embodiments of this disclosure or the prior art, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on these drawings without creative effort.

[0019] Figure 1 A flowchart illustrating a method for controllable source spectrum shaping according to some embodiments of the present disclosure is shown schematically.

[0020] Figure 2 A schematic comparison diagram of the excitation spectrum of explosives and the excitation spectrum of a controlled source is shown for some embodiments of this disclosure;

[0021] Figure 3A A schematic diagram of a time-domain sequence profile of the excitation spectrum of a concentrated explosive charge in a seismic channel, according to some embodiments of the present disclosure, is shown.

[0022] Figure 3B A schematic diagram of a time-domain sequence profile of the excitation spectrum of a controllable source in a seismic trace set according to some embodiments of the present disclosure is shown.

[0023] Figure 4 The diagram illustrates the results of taking the absolute value of the amplitude of the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum, respectively, according to some embodiments of the present disclosure.

[0024] Figure 5 The diagram illustrates, schematically, a first target frequency domain sequence after normalization of the explosive spectrum and a second target frequency domain sequence after normalization of the controllable source spectrum, according to some embodiments of the present disclosure.

[0025] Figure 6 A detailed implementation flowchart of step S130 of some embodiments of this disclosure is illustrated schematically;

[0026] Figure 7 The diagram schematically illustrates a comparison of the proposed coefficients within the target frequency band before and after adjustment, according to some embodiments of this disclosure.

[0027] Figure 8 An enlarged schematic diagram of the proposed coefficients in other frequency bands of some embodiments of this disclosure is shown schematically;

[0028] Figure 9The diagram schematically illustrates a comparison of the second shaped frequency domain sequence obtained by shaping according to some embodiments of the present disclosure, the second target frequency domain sequence before shaping, and the first target frequency domain sequence corresponding to the explosive excitation spectrum. The dashed line represents the first target frequency domain sequence corresponding to the explosive excitation spectrum, the solid line represents the second target frequency domain sequence before shaping, and the dotted line represents the second shaped frequency domain sequence obtained by shaping.

[0029] Figure 10 A structural block diagram of a controllable source spectrum shaping apparatus according to an embodiment of the present disclosure is schematically shown; and

[0030] Figure 11 A schematic block diagram of an electronic device provided in an embodiment of the present disclosure is shown. Detailed Implementation

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

[0032] The first exemplary embodiment of this disclosure provides a method for controllable source spectrum shaping.

[0033] Figure 1 A flowchart illustrating a method for controllable source spectrum shaping according to some embodiments of the present disclosure is shown schematically.

[0034] Reference Figure 1 As shown, the method for controllable source spectrum shaping provided in this embodiment includes the following steps: S110, S120, S130 and S140.

[0035] In step S110, the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum are determined; the above-mentioned explosive excitation spectrum and the above-mentioned controllable source excitation spectrum are the acquired signal spectra obtained by excitation at the same exploration location respectively.

[0036] Figure 2 The diagram schematically illustrates a comparison of the excitation spectrum of explosives and the excitation spectrum of a controlled source, representing some embodiments of this disclosure.

[0037] Reference Figure 2As shown, the frequency domain sequence of the explosive excitation spectrum (or described as a spectral curve) is represented by a solid line, and the frequency domain sequence of the controlled-source excitation spectrum is represented by a dashed line. Within the limited frequency band of the controlled-source, in band region A (7Hz–23Hz), the spectral curves of the controlled-source excitation spectrum and the explosive excitation spectrum highly overlap, and the corresponding visual imaging effects of the controlled-source excitation imaging profile and the explosive excitation imaging profile (e.g., bandpass 7Hz–23Hz imaging) are comparable. In band region B (23Hz–68Hz), the spectral curves of the controlled-source excitation spectrum and the explosive excitation spectrum differ significantly. This difference leads to differences in the visual imaging effects (e.g., bandpass 23Hz–68Hz imaging) between the controlled-source excitation imaging profile and the explosive excitation imaging profile.

[0038] Figure 3A A schematic diagram of a time-domain sequence profile of the excitation spectrum of a concentrated explosive charge in a seismic channel, according to some embodiments of the present disclosure, is shown. Figure 3B A schematic diagram of a time-domain sequence profile of the excitation spectrum of a controllable source in a seismic trace set according to some embodiments of the present disclosure is shown.

[0039] When comparing and analyzing explosive data and controlled-source data, the significant differences in the spectral curves of the controlled-source excitation spectrum and the explosive excitation spectrum in frequency band B affect the display of the controlled-source data. For example, comparing... Figure 3A and Figure 3B The time-domain sequences shown lack distinct characteristics of the controlled-source wavegroups, hindering geological structural analysis. In contrast, explosive data exhibits a higher signal-to-noise ratio, making it more suitable for geological structural analysis. Therefore, in this embodiment, the explosive excitation spectrum is used as an amplitude reference to shape the frequency domain data of the controlled-source excitation spectrum, achieving matched filtering of the controlled-source data.

[0040] In some embodiments, in step S110 above, determining the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum includes: loading the first time domain sequence of the explosive excitation spectrum and the second time domain sequence of the controllable source excitation spectrum; performing Fourier transform processing on the first time domain sequence and the second time domain sequence to obtain the corresponding first frequency domain sequence and second frequency domain sequence.

[0041] For example, a(t) is the first time-domain sequence of the explosive excitation spectrum, and b(t) is the second time-domain sequence of the controlled source excitation spectrum, where t represents time. By using Fourier transform (FFT) to transform the time-domain sequence into a frequency-domain sequence, we obtain the first frequency-domain sequence of the explosive excitation spectrum as a(f) and the second frequency-domain sequence of the controlled source excitation spectrum as b(f), where f represents frequency.

[0042] For example, in a specific instance, the acquisition gather of both the controlled-source and explosive-source sources contains 1814 data channels. Each data channel in these two types of source acquisition gathers is a time-domain sequence value. By reading one data channel from the controlled-source and the corresponding data channel from the explosive-source, and converting the time-domain sequence into a frequency-domain sequence using Fourier transform, the corresponding values ​​of frequency and amplitude for that channel are obtained. This yields the excitation spectrum of the explosive and the spectral curves of the controlled-source source, respectively.

[0043] In step S120, the first frequency domain sequence and the second frequency domain sequence are subjected to amplitude absolute value taking and amplitude normalization processing, respectively, to obtain the corresponding first target frequency domain sequence and second target frequency domain sequence.

[0044] Figure 4 The diagram illustrates the results of taking the absolute value of the amplitude of the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum, respectively, according to some embodiments of the present disclosure.

[0045] Reference Figure 4 As shown, in step S120, the absolute values ​​of the amplitude values ​​of the first frequency domain sequence a(f) and the second frequency domain sequence b(f) are taken to obtain the true amplitude spectrum curve A(f) of the explosive excitation data and the true amplitude spectrum curve B(f) of the controllable source excitation data.

[0046] Figure 5 The diagram illustrates a first target frequency domain sequence after normalization of the explosive spectrum and a second target frequency domain sequence after normalization of the controllable source spectrum, representing some embodiments of this disclosure.

[0047] Reference Figure 5 As shown, normalization operations are performed on the true amplitude spectrum curve A(f) of the explosive excitation data and the true amplitude spectrum curve B(f) of the controllable source excitation data to obtain the first target frequency domain sequence A after the explosive spectrum is normalized. normal (f) and the second target frequency domain sequence B after the controllable source spectrum normalization. normal (f). By taking the absolute value and normalizing the amplitudes of the first frequency domain sequence and the second frequency domain series of the controllable source and explosive at the same exploration location, the amplitude energy levels between the two can be unified and the relative amplitude differences can be easily compared.

[0048] For the i-th frequency f in the true amplitude spectrum curve A(f) of the explosive excitation data i The corresponding amplitude value, where i is the frequency index, and i takes the value of an integer; this frequency f i The corresponding amplitude value A(f) i Substitute the values ​​into the normalization formula for calculation:

[0049]

[0050] Among them, A min A is the minimum value in the true amplitude spectrum curve A(f) of the explosive excitation data. max The maximum value in the true amplitude spectrum curve A(f) of the explosive excitation data.

[0051] For the i-th frequency f in the true amplitude spectrum curve B(f) of the controlled source excitation data i The corresponding amplitude value, the frequency f i The corresponding amplitude value B(f) i Substitute the values ​​into the normalization formula for calculation:

[0052]

[0053] Among them, B min B is the minimum value in the true amplitude spectrum curve B(f) of the controlled source excitation data. max It represents the maximum value in the true amplitude spectrum curve B(f) of the controlled source excitation data.

[0054] In step S130, the amplitude adjustment coefficient is determined based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence.

[0055] In this embodiment of step S130, considering the different natural excitation widths of the two excitation sources—explosives and controllable seismic sources—explosives are broadband excitations, while seismic sources are finite-band excitations. Limited by the scanning frequency, superposition is a frequency reduction process, which can only be narrower than the designed bandwidth. Therefore, it is necessary to optimize the amplitude adjustment coefficient (as a shaping factor) in the high-frequency region to prevent aliasing at the boundaries of the spectrum during high-frequency scanning. Therefore, the embodiments of this disclosure further improve the process of determining the amplitude adjustment coefficient. The consistency verification between the proposed spectrum adjustment result and the bandwidth of the second target frequency domain sequence is used to determine whether to adjust the proposed coefficient within the frequency range. This ensures that after spectral shaping of the controllable seismic source data, the high signal-to-noise ratio frequency band signal is highlighted without affecting the data bandwidth, and avoids the possibility of aliasing in some high-frequency bands due to using the amplitude ratio as the amplitude adjustment coefficient for shaping across the entire frequency range.

[0056] Figure 6 A detailed implementation flowchart of step S130 of some embodiments of this disclosure is illustrated schematically.

[0057] In some embodiments, refer to Figure 6 As shown, in step S130 above, the amplitude adjustment coefficient is determined based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence, including the following steps: S610, S620, S630 and S640a.

[0058] In step S610, the amplitudes at the same frequency in the first target frequency domain sequence and the second target frequency domain sequence are divided to obtain the amplitude ratio at each frequency.

[0059] For example, at a frequency of 30Hz, A normal The amplitude value corresponding to (30Hz) is 0.4, B normal If the corresponding amplitude value for 30Hz is 0.8, then the amplitude ratio at 30Hz is 0.5. The processing for other frequencies is similar.

[0060] In step S620, the amplitude ratios at each of the above frequencies are used as the proposed coefficients, and the amplitudes at the corresponding frequencies of the second target frequency domain sequence are multiplied by the coefficients to obtain the proposed spectrum adjustment result.

[0061] For example, at a frequency of 30Hz, the second target frequency domain sequence B normal Multiply the amplitude of (30Hz) by 0.5 to obtain the amplitude adjustment result corresponding to 30Hz; multiply all the corresponding frequencies by the corresponding proposed coefficients to obtain the amplitude adjustment results for all frequencies, which is the proposed spectrum adjustment result.

[0062] In step S630, it is determined whether the proposed spectrum adjustment result is consistent with the bandwidth of the second target frequency domain sequence.

[0063] In step S640a, if the frequency is consistent throughout the entire frequency range, the proposed coefficients at each frequency within the entire frequency range are determined as the amplitude adjustment coefficients at the corresponding frequencies.

[0064] For example, the i-th frequency f i The proposed coefficient C(f) i The expression for ) is:

[0065] C(f i )= normal ( i ) / normal ( i (3)

[0066] The set coefficients for the entire frequency range constitute the sequence C(f) of amplitude adjustment coefficients.

[0067] In some embodiments, refer to Figure 6 As shown, the above step S130 also includes the following negation branches: S640b and S650b.

[0068] In step S640b, if there is an inconsistency in the target frequency segments, the proposed coefficients in the target frequency segments are adjusted until the proposed spectrum adjustment result in the target frequency segment corresponding to the adjusted proposed coefficients is consistent with the bandwidth of the second target frequency domain sequence.

[0069] Because the amplitude of the first target frequency domain sequence corresponding to the explosive is used as the target during the spectral amplitude adjustment process, and the natural excitation widths of the two excitation sources, explosive and controllable seismic source, are different. Explosive is a broadband excitation, while seismic source is a finite-band excitation, which is limited by the scanning frequency. Therefore, by verifying the consistency between the proposed spectral adjustment result and the bandwidth of the second target frequency domain sequence, it is possible to ensure that the spectral shaping of the controllable seismic source data can highlight the high signal-to-noise ratio frequency band signal without affecting the data bandwidth, and to avoid the possibility of false frequencies appearing in some high-frequency sections due to the use of amplitude ratio as the amplitude adjustment coefficient for shaping across the entire frequency range.

[0070] Figure 7 The diagram schematically illustrates a comparison of adjustments made to predetermined coefficients within a target frequency band according to some embodiments of this disclosure. Figure 7 To ensure font size and clear indication of the horizontal coordinates, the horizontal coordinate values ​​arranged in a row are indicated by vertical spacing. It can be understood that the numbers are still in ascending order to indicate frequency, and the misalignment is only for clarity.

[0071] In one embodiment, after the entire frequency range was shaped with a predetermined coefficient, a frequency scan of the pre-stack time-off profile revealed that a false frequency appeared at the high-pass 65Hz after the controllable source spectrum was shaped; at the same time, the bandwidth of the frequency band starting from the target frequency (65Hz) in the predetermined spectrum adjustment results was greater than the bandwidth corresponding to the second target frequency domain sequence mentioned above.

[0072] In some embodiments, if the bandwidth of the proposed spectrum adjustment result starting from the target frequency is greater than the bandwidth corresponding to the second target frequency domain sequence, then the frequency segment above the target frequency is determined as the target frequency segment. In this embodiment, for example, 65Hz to 123Hz is used as the target frequency segment. When the bandwidth of the proposed spectrum adjustment result within the target frequency segment is greater than the bandwidth of the second target frequency domain sequence, the proposed coefficients for the target frequency segment are reduced to obtain adjusted proposed coefficients. For example, refer to... Figure 7 As shown, the proposed coefficients within the target frequency band are adjusted, for example, in this embodiment, they are adjusted by reducing the frequency range.

[0073] In step S650b, the adjusted coefficients of the target frequency segment and the adjusted coefficients of other frequency segments are determined as the amplitude adjustment coefficients at the corresponding frequencies.

[0074] Figure 8 An enlarged schematic diagram of the proposed coefficients in other frequency bands of some embodiments of this disclosure is shown schematically.

[0075] Reference Figure 7 and Figure 8 As shown, the proposed coefficients in other frequency segments that do not require adjustment and the proposed coefficients after adjustment in the target frequency segment are determined as the amplitude adjustment coefficients for each corresponding frequency in the entire frequency range.

[0076] In step S140, the amplitude of the second target frequency domain sequence is shaped according to the amplitude adjustment coefficient to obtain the second shaped frequency domain sequence.

[0077] Figure 9 The diagram schematically illustrates a comparison of the second shaped frequency domain sequence obtained by shaping according to some embodiments of the present disclosure, the second target frequency domain sequence before shaping, and the first target frequency domain sequence corresponding to the explosive excitation spectrum. The dotted line represents the first target frequency domain sequence corresponding to the explosive excitation spectrum, the solid line represents the second target frequency domain sequence before shaping, and the dotted line represents the second shaped frequency domain sequence obtained by shaping.

[0078] For example, the expression for the second integer frequency domain sequence B'(f) is:

[0079] B'(f) = normal (f)×C(f), (4)

[0080] In formula (4), multiplying the sequences means multiplying the amplitudes at corresponding frequencies, and the results are arranged according to frequency to form a frequency domain curve.

[0081] In some embodiments, in step S140 above, shaping the amplitude of the second target frequency domain sequence according to the amplitude adjustment coefficient to obtain a second shaped frequency domain sequence includes: multiplying the amplitude adjustment coefficient and the amplitude of the second target frequency domain sequence at the corresponding frequency to obtain the second shaped frequency domain sequence. (Refer to...) Figure 9 As shown, matched filtering adjusts the relative amplitudes of different frequencies in controlled source data. The shaped amplitude curves are close to the spectrum of explosive excitation, which improves the signal-to-noise ratio of the controlled source data and reduces the differences in the comparison between controlled source data and explosive data. This helps structural interpreters to conduct geological structural analysis (such as stratigraphic tracing) by comparing and analyzing the shaped controlled source data and explosive data.

[0082] In some embodiments, the method for shaping the spectrum of the controllable source further includes: performing an inverse Fourier transform on the second shaped frequency domain sequence to obtain the shaped time domain sequence corresponding to the controllable source.

[0083] In one embodiment, the results of each channel are processed using the same method to obtain a second shaped frequency domain sequence of 1814 controllable sources. The result is then processed by inverse Fourier transform (IFFF) to obtain the shaped controllable source time domain sequence, which is then output to provide data support for geological structure and stratigraphic tracking analysis by structural analysts in the later stages.

[0084] A second exemplary embodiment of this disclosure provides an apparatus for controllable source spectrum shaping.

[0085] Figure 10 A structural block diagram of a controllable source spectrum shaping apparatus according to an embodiment of the present disclosure is shown schematically.

[0086] Reference Figure 10 As shown, the controllable source spectrum shaping device 1000 provided in this embodiment includes: a frequency sequence determination module 1001, a normalization processing module 1002, a shaping coefficient determination module 1003, and a shaping processing module 1004.

[0087] The frequency sequence determination module 1001 is used to determine the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum; the explosive excitation spectrum and the controllable source excitation spectrum are the acquired signal spectra obtained by excitation at the same exploration location.

[0088] The normalization processing module 1002 is used to perform amplitude absolute value taking and amplitude normalization processing on the first frequency domain sequence and the second frequency domain sequence respectively, to obtain the corresponding first target frequency domain sequence and second target frequency domain sequence.

[0089] The shaping coefficient determination module 1003 is used to determine the amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence.

[0090] The shaping module 1004 is used to shape the amplitude of the second target frequency domain sequence according to the amplitude adjustment coefficient to obtain the second shaped frequency domain sequence.

[0091] The aforementioned device 1000 also includes a time-frequency conversion module.

[0092] The aforementioned time-frequency conversion module is used to perform Fourier transform processing on the aforementioned first time-domain sequence and the aforementioned second time-domain sequence to obtain the corresponding first frequency-domain sequence and second frequency-domain sequence.

[0093] The aforementioned time-frequency conversion module is used to perform inverse Fourier transform processing on the aforementioned second shaped frequency domain sequence to obtain the shaped time domain sequence corresponding to the controllable source.

[0094] All relevant content from the first embodiment can be incorporated into this embodiment, and will not be repeated here.

[0095] Any plurality of the functional modules included in the aforementioned device 1000 may be combined into one module, or any one of the modules may be split into multiple modules. Alternatively, at least a portion of the functionality of one or more of these modules may be combined with at least a portion of the functionality of other modules and implemented in one module. At least one of the functional modules included in the aforementioned device 1000 may be at least partially implemented as hardware circuitry, such as a Field Programmable Gate Array (FPGA), a Programmable Logic Array (PLA), a System-on-Chip, a System-on-Substrate, a System-on-Package, an Application-Specific Integrated Circuit (ASIC), or any other reasonable means of integrating or packaging circuitry, or implemented in hardware or firmware, or in any one of software, hardware, and firmware implementations, or in a suitable combination of any of these. Alternatively, at least one of the functional modules included in the aforementioned device 1000 may be at least partially implemented as a computer program module, which, when run, can perform corresponding functions.

[0096] A third exemplary embodiment of this disclosure provides an electronic device.

[0097] Figure 11 The schematic diagram illustrates a structural block diagram of an electronic device provided in an embodiment of the present disclosure.

[0098] Reference Figure 11 As shown, the electronic device 1100 provided in this embodiment includes a processor 1101, a communication interface 1102, a memory 1103, and a communication bus 1104. The processor 1101, the communication interface 1102, and the memory 1103 communicate with each other through the communication bus 1104. The memory 1103 is used to store computer programs. When the processor 1101 executes the program stored in the memory, it implements the controllable source spectrum shaping method described above.

[0099] A fourth exemplary embodiment of this disclosure also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method for controllable source spectrum shaping as described above.

[0100] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments; or it may exist independently and not assembled into the device / apparatus. The computer-readable storage medium carries one or more programs that, when executed, implement the method according to the embodiments of this disclosure.

[0101] According to embodiments of this disclosure, the computer-readable storage medium can be a non-volatile computer-readable storage medium, such as including, but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this disclosure, the computer-readable storage medium can 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.

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

[0103] The above description is merely a specific embodiment of this disclosure, enabling those skilled in the art to understand or implement it. 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 this disclosure. Therefore, this disclosure is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A method for controllable source spectrum shaping, characterized in that, include: Determine the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum; The excitation spectrum of the explosive and the excitation spectrum of the controllable source are the acquired signal spectra obtained by excitation at the same exploration location; The first frequency domain sequence and the second frequency domain sequence are respectively subjected to amplitude absolute value taking and amplitude normalization processing to obtain the corresponding first target frequency domain sequence and second target frequency domain sequence; The amplitude adjustment coefficient is determined based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence. Based on the amplitude adjustment coefficient, the amplitude of the second target frequency domain sequence is shaped to obtain the second shaped frequency domain sequence.

2. The method according to claim 1, characterized in that, Based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence, the amplitude adjustment coefficient is determined, including: The amplitude ratio at each frequency is obtained by quotienting the amplitudes at the same frequency in the first target frequency domain sequence and the second target frequency domain sequence. Using the amplitude ratio at each frequency as a set coefficient, the amplitude at the corresponding frequency of the second target frequency domain sequence is multiplied by the coefficient to obtain the set spectrum adjustment result. The proposed spectrum adjustment result is compared to determine whether the bandwidth of the second target frequency domain sequence is consistent. When the frequency is consistent throughout the entire frequency range, the proposed coefficients at each frequency within the entire frequency range are determined as the amplitude adjustment coefficients at the corresponding frequencies.

3. The method according to claim 2, characterized in that, Determining the amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence further includes: In the event of inconsistencies in target frequency segments, the proposed coefficients within the target frequency segments are adjusted until the proposed spectrum adjustment result corresponding to the adjusted proposed coefficients in the target frequency segments is consistent with the bandwidth of the second target frequency domain sequence. The adjusted coefficients for the target frequency range and the adjusted coefficients for other frequency ranges are determined as the amplitude adjustment coefficients for the corresponding frequencies.

4. The method according to claim 3, characterized in that, If the proposed spectrum adjustment result shows a bandwidth greater than the bandwidth corresponding to the second target frequency domain sequence starting from the target frequency, then the frequency segment above the target frequency is determined as the target frequency segment.

5. The method according to claim 1, characterized in that, Based on the amplitude adjustment coefficient, the amplitude of the second target frequency domain sequence is shaped to obtain a second shaped frequency domain sequence, including: The amplitude adjustment coefficient and the amplitude of the second target frequency domain sequence at the corresponding frequency are multiplied to obtain the second shaped frequency domain sequence.

6. The method according to claim 1, characterized in that, The first frequency domain sequence for determining the excitation spectrum of the explosive and the second frequency domain sequence for the excitation spectrum of the controllable source include: The first time-domain sequence of the excitation spectrum of the loaded explosive and the second time-domain sequence of the excitation spectrum of the controllable source; Perform Fourier transform on the first time-domain sequence and the second time-domain sequence to obtain the first frequency-domain sequence and the second frequency-domain sequence respectively.

7. The method according to claim 1, characterized in that, Also includes: The second shaped frequency domain sequence is subjected to inverse Fourier transform processing to obtain the shaped time domain sequence corresponding to the controllable source.

8. A device for controllable source spectrum shaping, characterized in that, include: The frequency sequence determination module is used to determine the first frequency domain sequence of the explosive excitation spectrum and the second frequency domain sequence of the controllable source excitation spectrum. The excitation spectrum of the explosive and the excitation spectrum of the controllable source are the acquired signal spectra obtained by excitation at the same exploration location; The normalization processing module is used to perform amplitude absolute value taking and amplitude normalization processing on the first frequency domain sequence and the second frequency domain sequence respectively, to obtain the corresponding first target frequency domain sequence and second target frequency domain sequence. The shaping coefficient determination module is used to determine the amplitude adjustment coefficient based on the amplitude ratio of the first target frequency domain sequence and the second target frequency domain sequence; The shaping module is used to shape the amplitude of the second target frequency domain sequence according to the amplitude adjustment coefficient to obtain the second shaped frequency domain sequence.

9. An electronic device, characterized in that, It includes 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; Memory, used to store computer programs; A processor, when executing a program stored in memory, implements the method of any one of claims 1-7.

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

Images