Seismic data fusion method, device, equipment and medium

Abstract

The invention provides a seismic data fusion method, apparatus and device, and a readable medium. The method comprises the steps of obtaining original seismic data of different frequency bands; respectively converting the original seismic data of different frequency bands into a frequency domain; calculating an energy correction factor between the seismic data of different frequency bands based on the original seismic data of different frequency bands; based on the energy correction factor obtained through calculation, broadband seismic data in a frequency domain are calculated through an inversion method; and converting the broadband seismic data in the frequency domain into broadband seismic data in the time domain. By using the scheme of the invention, seismic data of different frequency bands can be effectively fused together so as to obtain broadband seismic data, and the precision of seismic data fusion can be improved.

Description

Technical Field

[0001] The present invention relates to the field of geophysical exploration, and more particularly to a method, apparatus, device and readable medium for seismic data fusion. Background Art

[0002] In early seismic exploration, explosives were the main excitation sources. In recent years, with the increasingly strict requirements for safety and environmental protection, green, safe and environmentally friendly excitation sources have received more and more attention. The vibrator has many advantages over the explosive source, such as controllable frequency energy, safety and environmental protection, low cost, high construction efficiency, etc., and is used more and more widely. The vibrator excitation is to excite a series of signals with changing frequencies by the source, which are reflected back to the surface through the underground and recorded by the geophones, and then the recorded signals are correlated with the theoretical reference signals to obtain seismic data. With the refinement of exploration targets, the seismic exploration has higher and higher requirements for the frequency band, and wide-band exploration has become the mainstream technology. The high-efficiency acquisition technology of the vibrator is an important means to achieve the "two-width and one-high" technology. Due to the limitation of the vibrator manufacturing process, it is difficult to simultaneously take into account low frequencies and high frequencies in the manufacture of the source. To solve this problem, the frequency division acquisition technology has been proposed. By splitting the conventional wide-band signal into multiple narrow-band signals such as low frequency, medium frequency, and high frequency, and finally obtaining the wide-band signal through subsequent processing. The traditional frequency division acquisition data is fused by direct superposition. Since there are overlapping frequencies in different frequency bands, direct superposition will reduce the accuracy of the data. Summary of the Invention

[0003] In view of this, the purpose of the embodiments of the present invention is to provide a method, apparatus, device and readable medium for seismic data fusion. By using the technical solution of the present invention, seismic data in different frequency bands can be effectively fused together to obtain wide-band seismic data, and the accuracy of seismic data fusion can be improved.

[0004] Based on the above purpose, an aspect of the embodiments of the present invention provides a method for seismic data fusion, including the following steps:

[0005] Obtain the original seismic data in different frequency bands;

[0006] Convert the original seismic data in different frequency bands into the frequency domain respectively;

[0007] Calculate the energy correction factor between the seismic data in different frequency bands based on the original seismic data in different frequency bands;

[0008] Calculate the wide-band seismic data in the frequency domain based on the calculated energy correction factor and by the method of inversion;

[0009] Convert the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain.

[0010] According to an embodiment of the present invention, obtaining the original seismic data of different frequency bands includes:

[0011] Setting the highest frequency, the lowest frequency of the first frequency band, and the first scanning time of the seismic data, and scanning the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the first frequency band;

[0012] Setting the highest frequency, the lowest frequency of the second frequency band, and the second scanning time of the seismic data, and scanning the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the second frequency band, wherein the frequency parts in the first frequency band and the second frequency band overlap.

[0013] According to an embodiment of the present invention, respectively converting the original seismic data of different frequency bands into the frequency domain includes:

[0014] Performing one-dimensional Fourier transform on the original seismic data of different frequency bands respectively to convert the original seismic data into the frequency domain.

[0015] According to an embodiment of the present invention, calculating the energy correction factor between the seismic data of different frequency bands based on the original seismic data of different frequency bands includes:

[0016] Using the formula: Calculating the energy correction factor between the seismic data of different frequency bands, where c is the energy correction factor of the seismic data of two different frequency bands, F i 1 and F i 2 respectively represent the amplitude spectra of the two different frequency band data at frequency i, and i1 and i2 represent the starting frequency and the ending frequency of the overlapping frequency of the two different frequency band data.

[0017] According to an embodiment of the present invention, it further includes:

[0018] Constructing low-frequency and high-frequency filters, and the filter is expressed as:

[0019] where the diagonal elements in the matrix E are non-zero, and other elements are all zero, and the number and positions of the non-zero elements are determined by the frequency bands of the frequency division scanning signal.

[0020] According to an embodiment of the present invention, calculating the broadband seismic data in the frequency domain based on the calculated energy correction factor and by the inversion method includes:

[0021] Using the formula: and M = (f h 2 + c 2 f l 2 )-1 (f h T R h +cf l T R l ) Calculate the broadband seismic data in the frequency domain, where f l represents the filtering factor of the low-frequency band data, and f h represents the filtering factor of the high-frequency band data, and f l and f h can be represented by the matrix E, M represents the broadband seismic data in the frequency domain, and R l represents the narrowband seismic data of the low-frequency band in the frequency domain, and R h represents the narrowband seismic data of the high-frequency band in the frequency domain, c is the energy correction factor for the data of two different frequency bands, and f h T represents the transpose of the filtering factor of the high-frequency band data, and f l T represents the transpose of the filtering factor of the low-frequency band data.

[0022] According to an embodiment of the present invention, converting the broadband seismic data in the frequency domain into the broadband seismic data in the time domain includes:

[0023] Performing a one-dimensional inverse Fourier transform on the broadband seismic data in the frequency domain to convert the broadband seismic data in the frequency domain into the broadband seismic data in the time domain.

[0024] Another aspect of the embodiments of the present invention further provides a device for seismic data fusion, and the device includes:

[0025] An acquisition module, configured to acquire the original seismic data of different frequency bands;

[0026] A first conversion module, configured to respectively convert the original seismic data of different frequency bands into the frequency domain;

[0027] A calculation module, configured to calculate the energy correction factor between the seismic data of different frequency bands based on the original seismic data of different frequency bands;

[0028] An obtaining module, configured to calculate the broadband seismic data in the frequency domain based on the calculated energy correction factor and by means of inversion;

[0029] A second conversion module, configured to convert the broadband seismic data in the frequency domain into the broadband seismic data in the time domain.

[0030] Another aspect of the embodiments of the present invention further provides a computer device, and the computer device includes:

[0031] At least one processor; and

[0032] A memory that stores computer instructions executable on the processor, and when the instructions are executed by the processor, the steps of any one of the above methods are implemented.

[0033] Another aspect of an embodiment of the present invention further provides a computer-readable storage medium that stores a computer program, and when the computer program is executed by a processor, the steps of any one of the above methods are implemented.

[0034] The present invention has the following beneficial technical effects: The method for seismic data fusion provided by the embodiments of the present invention, through the technical solutions of acquiring original seismic data in different frequency bands, respectively converting the original seismic data in different frequency bands into the frequency domain, calculating the energy correction factor between the seismic data in different frequency bands based on the original seismic data in different frequency bands, calculating the broadband seismic data in the frequency domain based on the calculated energy correction factor and through an inversion method, and converting the broadband seismic data in the frequency domain into the broadband seismic data in the time domain, can effectively fuse the seismic data in different frequency bands to obtain broadband seismic data, and can improve the accuracy of seismic data fusion. Description of the Drawings

[0035] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention, and those of ordinary skill in the art can also obtain other embodiments based on these drawings without creative efforts.

[0036] Figure 1 It is a schematic flowchart of a method for seismic data fusion according to an embodiment of the present invention;

[0037] Figure 2 It is a schematic diagram of a low-frequency sweep signal and its spectrum in a frequency band according to an embodiment of the present invention;

[0038] Figure 3 It is a schematic diagram of a high-frequency sweep signal and its spectrum in a frequency band according to an embodiment of the present invention;

[0039] Figure 4 It is a schematic diagram of a single-shot cross-well data acquisition and data fusion result according to an embodiment of the present invention;

[0040] Figure 5 It is a schematic diagram of a cross-well data acquisition profile and data fusion result according to an embodiment of the present invention;

[0041] Figure 6Schematic diagram of the wavelet and spectrum after inversion data fusion according to an embodiment of the present invention;

[0042] Figure 7 Schematic diagram of the wavelet and spectrum of direct superposition fusion in the time domain according to an embodiment of the present invention;

[0043] Figure 8 Schematic diagram of the device for seismic data fusion according to an embodiment of the present invention;

[0044] Figure 9 Schematic diagram of the computer device according to an embodiment of the present invention;

[0045] Figure 10 Schematic diagram of the computer-readable storage medium according to an embodiment of the present invention. Detailed implementation manners

[0046] To make the objectives, technical solutions and advantages of the present invention clearer and more understandable, the following further elaborates on the embodiments of the present invention in detail with reference to specific embodiments and the accompanying drawings.

[0047] Based on the above objectives, the first aspect of the embodiments of the present invention proposes an embodiment of a method for seismic data fusion. Figure 1 The schematic flowchart of the method is shown.

[0048] As Figure 1 shown, the method may include the following steps:

[0049] S1 Obtain the original seismic data of different frequency bands. Set the highest frequency, the lowest frequency of the first frequency band, and the first scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the first frequency band. Set the highest frequency, the lowest frequency of the second frequency band, and the second scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the second frequency band, where the frequency parts in the first frequency band and the second frequency band overlap. Of course, seismic data of more frequency bands can also be obtained, and the method is the same as the above.

[0050] S2 Convert the original seismic data of different frequency bands into the frequency domain respectively. Perform one-dimensional Fourier transform on the original seismic data of different frequency bands respectively to convert the original seismic data into the frequency domain.

[0051] S3 Calculate the energy correction factor between the seismic data of different frequency bands based on the original seismic data of different frequency bands. The formula: can be used to calculate the energy correction factor between the seismic data of different frequency bands, where c is the energy correction factor of the seismic data of two different frequency bands, F i 1 and Fi 2 The distribution represents the amplitude spectrum with frequency i of data in two different frequency bands. i1 and i2 represent the starting frequency and the ending frequency of the overlapping frequency of data in the two different frequency bands.

[0052] S4 calculates the broadband seismic data in the frequency domain based on the calculated energy correction factor and by means of inversion. Before the inversion calculation, filters for low frequency and high frequency also need to be constructed. The filters can be expressed as:

[0053] Where the diagonal elements in matrix E are non-zero and other elements are all zero. The number and positions of the non-zero elements are determined by the frequency bands of the frequency division scanning signal. Then the formula: and M = (f h 2 + c 2 f l 2 ) -1 (f h T R h + cf l T R l is used to calculate the broadband seismic data in the frequency domain, where f l represents the filtering factor of the low-frequency band data, f h represents the filtering factor of the high-frequency band data, f l and f h can be represented by matrix E. M represents the broadband seismic data in the frequency domain. R l represents the narrowband seismic data of the low-frequency band in the frequency domain. R h represents the narrowband seismic data of the high-frequency band in the frequency domain. c is the energy correction factor of data in two different frequency bands. f h T represents the transpose of the filtering factor of the high-frequency band data, and f l T represents the transpose of the filtering factor of the low-frequency band data.

[0054] S5 converts the broadband seismic data in the frequency domain into broadband seismic data in the time domain. The broadband seismic data in the frequency domain is subjected to one-dimensional inverse Fourier transform to convert the broadband seismic data in the frequency domain into broadband seismic data in the time domain.

[0055] By using the technical solution of the present invention, seismic data in different frequency bands can be effectively fused together to obtain broadband seismic data, and the accuracy of seismic data fusion can be improved.

[0056] In a preferred embodiment of the present invention, obtaining the original seismic data in different frequency bands includes:

[0057] Set the highest frequency, lowest frequency of the first frequency band, and the first scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the first frequency band;

[0058] Set the highest frequency, lowest frequency of the second frequency band, and the second scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the second frequency band, where the frequency parts in the first frequency band and the second frequency band overlap. The original seismic data of different frequency bands obtained can be prestack data or post-stack seismic data. For example, to obtain seismic data of two different frequency bands, data collected separately by two different frequency bands excited by a vibrator can be used. The frequency range collected in the low frequency band is 1.5 - 12 Hz, and the scanning length is 12 s, as Figure 2 shown, the frequency range collected in the medium and high frequency band is 6 - 96 Hz, and the scanning length is 16 s, as Figure 3 shown, the overlapping frequency of the scanning signals of the two different frequency bands is 6 - 12 Hz, and the original single-shot data collected is Figure 4 shown, after performing stacking processing on the original single-shot data collected, the stacked profile of the original frequency division acquisition is obtained, as Figure 5 shown.

[0059] In a preferred embodiment of the present invention, respectively converting the original seismic data of different frequency bands into the frequency domain includes:

[0060] Respectively perform one-dimensional Fourier transform on the original seismic data of different frequency bands to convert the original seismic data into the frequency domain. Convert the original narrow-band seismic data into the frequency domain. For each trace of the prestack or post-stack data, it is transformed into the frequency domain through one-dimensional Fourier transform. Generally, for the original data of different frequency bands, it is necessary to limit its highest frequency and lowest frequency, and the highest and lowest frequencies of different frequency bands should be consistent with the highest frequency and lowest frequency of their scanning signals.

[0061] In a preferred embodiment of the present invention, calculating the energy correction factor between the seismic data of different frequency bands based on the original seismic data of different frequency bands includes:

[0062] Use the formula: Calculate the energy correction factor between the seismic data of two different frequency bands, where c is the energy correction factor of the seismic data of two different frequency bands, F i 1 and F i 2The distribution represents the amplitude spectrum with frequency i of data in two different frequency bands. i1 and i2 represent the starting frequency and the ending frequency of the overlapping frequency of data in the two different frequency bands. Calculate the energy correction factors for data in different frequency bands according to the above formula. For the above example, the theoretical energy correction factors for low-frequency and medium-high-frequency data are 0.1555. For the actually collected data, due to the filtering effect of the earth, there are differences between the calculated energy correction factors and the theoretical factors. An energy correction factor can be calculated for each trace, or an average energy correction factor can be applied to the entire gather data. For the above example, an energy correction factor of 0.182 is applied to the entire gather.

[0063] In a preferred embodiment of the present invention, it further includes:

[0064] Construct filters for low frequency and high frequency. The filters are represented as:

[0065] where the diagonal elements in matrix E are non-zero, and other elements are all zero. The number and positions of the non-zero elements are determined by the frequency bands of the frequency division scanning signals.

[0066] In a preferred embodiment of the present invention, calculating the wide-band seismic data in the frequency domain based on the calculated energy correction factors and through an inversion method includes:

[0067] Using the formula: and M = (f h 2 + c 2 f l 2 ) -1 (f h T R h + cf l T R l ) to calculate the wide-band seismic data in the frequency domain, where f l represents the filtering factor of the low-frequency band data, f h represents the filtering factor of the high-frequency band data, f l and f h can be represented by matrix E, M represents the wide-band seismic data in the frequency domain, R l represents the narrow-band seismic data of the low-frequency band in the frequency domain, R h represents the narrow-band seismic data of the high-frequency band in the frequency domain, c is the energy correction factor for data in two different frequency bands, f h T represents the transpose of the filtering factor of the high-frequency band data, f l TThe transpose of the filtering factor representing the low-frequency band data. Different narrow-band seismic data can be approximately regarded as filtering the wide-band seismic data, and this process can be represented by the formula which can be obtained by calculating the formula M = (f h 2 +c 2 f l 2 ) -1 (f h T R h +cf l T R l ) to obtain the wide-band seismic data in the frequency domain. To obtain the wide-band seismic data in the frequency domain by the inversion method, it is first necessary to construct low-frequency and high-frequency filters, which can be represented by the matrix E shown in the above formula. Except for the non-zero elements on the diagonal part, other elements are all 0. The number and position of the non-zero elements are determined by the frequency band of the frequency division scanning signal. After constructing the filter, the wide-band seismic data can be inverted by the above formula.

[0068] In a preferred embodiment of the present invention, converting the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain includes:

[0069] Performing one-dimensional inverse Fourier transform on the wide-band seismic data in the frequency domain to convert the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain. The result inverted above is in the frequency domain. By performing one-dimensional inverse Fourier transform on it, the wide-band seismic data in the time domain can be obtained. Figure 6 are the wide-band seismic wavelets and spectra obtained by the inversion method, Figure 7 are the fused wide-band wavelets and spectra obtained by direct superposition. By comparing the two, it can be seen that the influence of overlapping frequencies is eliminated by the inversion method, and the accuracy and resolution of the wide-band wavelets are improved. Figure 4 c and Figure 5 c are respectively the fused wide-band single-shot data and the wide-band stacked section obtained after inversion.

[0070] By using the technical solution of the present invention, it is possible to effectively fuse seismic data of different frequency bands to obtain wide-band seismic data, and improve the accuracy of seismic data fusion.

[0071] It should be noted that those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The above program can be stored in a computer-readable storage medium. When the program is executed, it can include the processes of the above method embodiments. The storage medium can be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM), etc. The embodiments of the above computer program can achieve the same or similar effects as the corresponding foregoing method embodiments.

[0072] In addition, the method disclosed according to the embodiments of the present invention can also be implemented as a computer program executed by a CPU, and the computer program can be stored in a computer-readable storage medium. When the computer program is executed by the CPU, the above functions defined in the method disclosed in the embodiments of the present invention are executed.

[0073] Based on the above object, in the second aspect of the embodiments of the present invention, a device for seismic data fusion is proposed, as Figure 8 shown, the device 200 includes:

[0074] An acquisition module, configured to acquire original seismic data of different frequency bands;

[0075] A first conversion module, configured to respectively convert the original seismic data of different frequency bands into the frequency domain;

[0076] A calculation module, configured to calculate the energy correction factor between the seismic data of different frequency bands based on the original seismic data of different frequency bands;

[0077] An obtaining module, configured to calculate the broadband seismic data in the frequency domain based on the calculated energy correction factor by an inversion method;

[0078] A second conversion module, configured to convert the broadband seismic data in the frequency domain into broadband seismic data in the time domain.

[0079] Based on the above object, in the third aspect of the embodiments of the present invention, a computer device is proposed. Figure 9 The following shows a schematic diagram of an embodiment of the computer device provided by the present invention. As Figure 9 shown, the embodiments of the present invention include the following devices: at least one processor 21; and a memory 22, the memory 22 stores computer instructions 23 that can run on the processor, and when the instructions are executed by the processor, the following methods are implemented:

[0080] Acquire original seismic data of different frequency bands;

[0081] Convert the original seismic data in different frequency bands into the frequency domain respectively;

[0082] Calculate the energy correction factor between the seismic data in different frequency bands based on the original seismic data in different frequency bands;

[0083] Calculate the broadband seismic data in the frequency domain based on the calculated energy correction factor and by means of inversion;

[0084] Convert the broadband seismic data in the frequency domain into the broadband seismic data in the time domain.

[0085] In a preferred embodiment of the present invention, obtaining the original seismic data in different frequency bands includes:

[0086] Set the highest frequency, the lowest frequency of the first frequency band, and the first scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the first frequency band;

[0087] Set the highest frequency, the lowest frequency of the second frequency band, and the second scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the second frequency band, wherein the frequencies in the first frequency band and the second frequency band partially overlap.

[0088] In a preferred embodiment of the present invention, converting the original seismic data in different frequency bands into the frequency domain respectively includes:

[0089] Perform one-dimensional Fourier transform on the original seismic data in different frequency bands respectively to convert the original seismic data into the frequency domain.

[0090] In a preferred embodiment of the present invention, calculating the energy correction factor between the seismic data in different frequency bands based on the original seismic data in different frequency bands includes:

[0091] Use the formula: Calculate the energy correction factor between the seismic data in different frequency bands, where c is the energy correction factor of the seismic data in two different frequency bands, F i 1 and F i 2 respectively represent the amplitude spectra with frequency i of the data in two different frequency bands, and i1 and i2 represent the starting frequency and the ending frequency of the overlapping frequencies of the data in two different frequency bands.

[0092] In a preferred embodiment of the present invention, it further includes:

[0093] Construct low-frequency and high-frequency filters, and the filter is expressed as:

[0094] Among them, the diagonal elements in matrix E are non-zero, and the other elements are all zero. The number and positions of the non-zero elements are determined by the frequency band of the frequency division scanning signal.

[0095] In a preferred embodiment of the present invention, calculating the wide-band seismic data in the frequency domain based on the calculated energy correction factor and by means of inversion includes:

[0096] Using the formula: and M = (f h 2 + c 2 f l 2 ) -1 (f h T R h + cf l T R l ) to calculate the wide-band seismic data in the frequency domain, where f l represents the filtering factor of the low-frequency band data, f h represents the filtering factor of the high-frequency band data, f l and f h can be represented by matrix E, M represents the wide-band seismic data in the frequency domain, R l represents the narrow-band seismic data of the low-frequency band in the frequency domain, R h represents the narrow-band seismic data of the high-frequency band in the frequency domain, c is the energy correction factor of the two different frequency band data, f h T represents the transpose of the filtering factor of the high-frequency band data, f l T represents the transpose of the filtering factor of the low-frequency band data.

[0097] In a preferred embodiment of the present invention, converting the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain includes:

[0098] Performing a one-dimensional inverse Fourier transform on the wide-band seismic data in the frequency domain to convert the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain.

[0099] For the above purpose, a fourth aspect of the embodiments of the present invention proposes a computer-readable storage medium. Figure 10 Shown is a schematic diagram of an embodiment of the computer-readable storage medium provided by the present invention. As Figure 10 shown, the computer-readable storage medium 31 stores a computer program 32 that, when executed by a processor, executes the following method:

[0100] Obtain the original seismic data of different frequency bands;

[0101] Convert the original seismic data in different frequency bands into the frequency domain respectively;

[0102] Calculate the energy correction factor between the seismic data in different frequency bands based on the original seismic data in different frequency bands;

[0103] Calculate the broadband seismic data in the frequency domain based on the calculated energy correction factor and by means of inversion;

[0104] Convert the broadband seismic data in the frequency domain into the broadband seismic data in the time domain.

[0105] In a preferred embodiment of the present invention, obtaining the original seismic data in different frequency bands includes:

[0106] Set the highest frequency, the lowest frequency of the first frequency band and the first scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the first frequency band;

[0107] Set the highest frequency, the lowest frequency of the second frequency band and the second scanning time of the seismic data, and scan the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the second frequency band, wherein the frequencies in the first frequency band and the second frequency band partially overlap.

[0108] In a preferred embodiment of the present invention, converting the original seismic data in different frequency bands into the frequency domain respectively includes:

[0109] Perform one-dimensional Fourier transform on the original seismic data in different frequency bands respectively to convert the original seismic data into the frequency domain.

[0110] In a preferred embodiment of the present invention, calculating the energy correction factor between the seismic data in different frequency bands based on the original seismic data in different frequency bands includes:

[0111] Use the formula: Calculate the energy correction factor between the seismic data in different frequency bands, where c is the energy correction factor of the seismic data in two different frequency bands, F i 1 and F i 2 respectively represent the amplitude spectra with frequency i of the data in two different frequency bands, and i1 and i2 represent the starting frequency and the ending frequency of the overlapping frequencies of the data in two different frequency bands.

[0112] In a preferred embodiment of the present invention, it further includes:

[0113] Construct filters for low frequency and high frequency, and the filters are expressed as:

[0114] Among them, the diagonal elements in matrix E are non-zero, and other elements are all zero. The number and positions of the non-zero elements are determined by the frequency band of the frequency division scanning signal.

[0115] In a preferred embodiment of the present invention, calculating the wide-band seismic data in the frequency domain based on the calculated energy correction factor and by means of inversion includes:

[0116] Using the formula: and M=(f h 2 +c 2 f l 2 ) -1 (f h T R h +cf l T R l ) to calculate the wide-band seismic data in the frequency domain, where f l represents the filtering factor of the low-frequency band data, f h represents the filtering factor of the high-frequency band data, f l and f h can be represented by matrix E, M represents the wide-band seismic data in the frequency domain, R l represents the narrow-band seismic data of the low-frequency band in the frequency domain, R h represents the narrow-band seismic data of the high-frequency band in the frequency domain, c is the energy correction factor of the two different frequency band data, f h T represents the transpose of the filtering factor of the high-frequency band data, f l T represents the transpose of the filtering factor of the low-frequency band data.

[0117] In a preferred embodiment of the present invention, converting the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain includes:

[0118] Performing a one-dimensional inverse Fourier transform on the wide-band seismic data in the frequency domain to convert the wide-band seismic data in the frequency domain into the wide-band seismic data in the time domain.

[0119] In addition, the method disclosed according to the embodiments of the present invention can also be implemented as a computer program executed by a processor, and this computer program can be stored in a computer-readable storage medium. When this computer program is executed by the processor, the above functions defined in the method disclosed according to the embodiments of the present invention are executed.

[0120] In addition, the above method steps and system units can also be implemented by using a controller and a computer-readable storage medium for storing a computer program that enables the controller to implement the above step or unit functions.

[0121] Those skilled in the art will also understand that the various exemplary logical blocks, modules, circuits, and algorithmic steps described in connection with the disclosure herein can be implemented as electronic hardware, computer software, or a combination of both. To clearly illustrate this interchangeability of hardware and software, a general description has been given of the functionality of the various illustrative components, blocks, modules, circuits, and steps. Whether this functionality is implemented as software or hardware depends upon the particular application and the design constraints imposed on the overall system. The functionality that can be implemented in various ways by those skilled in the art for each particular application, but such implementation decisions should not be construed as causing a departure from the scope of the disclosure of the embodiments of the present invention.

[0122] In one or more exemplary designs, the functionality can be implemented in hardware, software, firmware, or any combination thereof. If implemented in software, the functionality can be stored on or transmitted via a computer-readable medium as one or more instructions or code. Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one location to another. The storage media can be any available media that can be accessed by a general or special purpose computer. By way of example and not limitation, the computer-readable media can include RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and that can be accessed by a general or special purpose computer or a general or special purpose processor. Additionally, any connection is properly termed a computer-readable medium. For example, if software is transmitted from a website, server, or other remote source using a coaxial cable, fiber optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave, then the coaxial cable, fiber optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave are included in the definition of medium. As used herein, disk and disc include compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk, and Blu-ray disc where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above should also be included within the scope of computer-readable media.

[0123] The above are exemplary embodiments disclosed in the present invention, but it should be noted that various changes and modifications can be made without departing from the scope of the disclosure of the embodiments of the present invention as defined by the claims. The functions, steps, and / or actions of the method claims according to the disclosed embodiments herein need not be performed in any particular order. Additionally, although the elements of the embodiments disclosed in the present invention may be described or claimed in individual form, they can also be understood as plural unless explicitly limited to the singular.

[0124] It should be understood that, as used herein, unless the context clearly supports the exception, the singular form "a" is intended to also include the plural form. It should also be understood that the "and / or" used herein refers to any and all possible combinations of one or more of the associated listed items.

[0125] The serial numbers of the disclosed embodiments of the present invention above are only for description and do not represent the advantages or disadvantages of the embodiments.

[0126] Those of ordinary skill in the art can understand that all or part of the steps to implement the above embodiments can be completed by hardware, or can be completed by instructing relevant hardware through a program, and the program can be stored in a computer-readable storage medium. The above-mentioned storage medium can be a read-only memory, a magnetic disk or an optical disc, etc.

[0127] Those of ordinary skill in the art should understand that: the discussion of any of the above embodiments is only exemplary and is not intended to imply that the scope of the disclosure of the embodiments of the present invention (including the claims) is limited to these examples; under the concept of the embodiments of the present invention, the technical features between the above embodiments or different embodiments can also be combined, and there are many other variations in different aspects of the embodiments of the present invention as above, which are not provided in detail for the sake of brevity. Therefore, any omission, modification, equivalent replacement, improvement, etc. made within the spirit and principle of the embodiments of the present invention shall be included in the protection scope of the embodiments of the present invention.

Claims

1. A method for seismic data fusion, characterized in that, The following steps are involved: Obtain raw seismic data in different frequency bands; Convert the original seismic data of different frequency bands into the frequency domain respectively; Calculating energy correction factors between seismic data of different frequency bands based on original seismic data of different frequency bands; Based on the calculated energy correction factor, broadband seismic data in the frequency domain are calculated by inversion method; Convert broadband seismic data in the frequency domain into broadband seismic data in the time domain.

2. The method according to claim 1, characterized in that, Obtaining raw seismic data in different frequency bands includes: Setting the highest frequency, the lowest frequency and the first scanning time of the seismic data of the first frequency band, and scanning the seismic data between the highest frequency and the lowest frequency to obtain the original seismic data of the first frequency band; The highest frequency, the lowest frequency and the second scanning time of the seismic data of the second frequency band are set, and the seismic data between the highest frequency and the lowest frequency are scanned to obtain the original seismic data of the second frequency band, wherein the frequencies in the first frequency band and the second frequency band partially overlap.

3. The method according to claim 1, wherein Converting the original seismic data of different frequency bands into the frequency domain includes: One-dimensional Fourier transform is performed on the original seismic data of different frequency bands to convert the original seismic data into the frequency domain.

4. The method according to claim 1, characterized in that, Calculating the energy correction factor between seismic data of different frequency bands based on the original seismic data of different frequency bands includes: Use the formula: Calculate the energy correction factor between seismic data of different frequency bands, where c is the energy correction factor of seismic data of two different frequency bands, F i 1 and F i 2 distributions represent the amplitude spectra of frequency i of data of two different frequency bands, and i1 and i2 represent the starting frequency and ending frequency of the overlapping frequency of data of two different frequency bands.

5. The method according to claim 1, characterized in that, Also includes: Construct low-frequency and high-frequency filters. The filters are expressed as: Among them, the diagonal elements in matrix E are non-zero, and other elements are all zero. The number and positions of the non-zero elements are determined by the frequency band of the frequency division scanning signal.

6. The method according to claim 5, wherein Based on the calculated energy correction factor and by inversion method, broadband seismic data in the frequency domain are calculated including: Use the formula: and to calculate the broadband seismic data in the frequency domain, where f l represents the filtering factor of the low-frequency band data, f h represents the filtering factor of the high-frequency band data, f l and f h can be represented by the matrix E, M represents the broadband seismic data in the frequency domain, R l represents the narrowband seismic data of the low-frequency band in the frequency domain, R h represents the narrowband seismic data of the high-frequency band in the frequency domain, c is the energy correction factor of the two different frequency band data, f h T represents the transpose of the filtering factor of the high-frequency band data, f l T represents the transpose of the filtering factor of the low-frequency band data.

7. The method according to claim 1, characterized in that, Converting broadband seismic data in the frequency domain into broadband seismic data in the time domain includes: The broadband seismic data in the frequency domain is subjected to a one-dimensional inverse Fourier transform to convert the broadband seismic data in the frequency domain into broadband seismic data in the time domain.

8. An apparatus for seismic data fusion, characterized in that, The device comprises: An acquisition module, wherein the acquisition module is configured to acquire raw seismic data of different frequency bands; A first conversion module, wherein the first conversion module is configured to convert raw seismic data of different frequency bands into a frequency domain respectively; A calculation module, wherein the calculation module is configured to calculate an energy correction factor between seismic data of different frequency bands based on original seismic data of different frequency bands; An obtaining module, the obtaining module being configured to calculate broadband seismic data in the frequency domain by an inversion method based on the calculated energy correction factor; A second conversion module, wherein the second conversion module is configured to convert broadband seismic data in the frequency domain into broadband seismic data in the time domain.

9. A computer device, characterized in that, include: at least one processor; as well as A memory storing computer instructions executable on the processor, wherein the instructions, when executed by the processor, implement the steps of the method according to any one of claims 1 to 7.

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

Images