A mapping amplitude preserving method, apparatus and device

Abstract

Embodiments of the present specification provide a mapping amplitude preserving method, device and equipment, wherein the method comprises: obtaining a first seismic signal; wherein the first seismic signal is an original seismic signal; obtaining a second seismic signal; wherein the second seismic signal is a seismic signal after high resolution processing on the original seismic signal; performing time-frequency analysis on the first seismic signal to obtain an effective frequency band; establishing a mapping relationship in amplitude between the first seismic signal and the second seismic signal within the effective frequency band range; and performing amplitude preserving processing on a full frequency band signal of the second seismic signal according to the mapping relationship to obtain a seismic signal after amplitude preserving processing. In the embodiments of the present specification, the mapping relationship in amplitude between the seismic signals before and after high resolution processing within the effective frequency band range can be used to accurately restore the relative amplitude of the seismic signal after high resolution processing.

Description

TECHNICAL FIELD

[0001] Embodiments of the present specification relate to the technical field of seismic data processing, and in particular to a mapping amplitude preservation method, device and equipment. BACKGROUND

[0002] In actual seismic exploration, there are many factors affecting the amplitude of seismic waves, including excitation, reception, absorption and processing, etc. If one aspect goes wrong, it will be difficult to obtain high-fidelity seismic imaging results. In view of these aspects, the relevant amplitude processing methods are divided into two categories: true amplitude recovery technology and relative amplitude processing technology. Among them, the relative amplitude processing technology is to eliminate the amplitude information irrelevant to the geological information in the seismic data processing process, and to maintain the relative change trend of the amplitude in time and space to reflect the relative change relationship of the energy of the seismic wave propagating in the stratum.

[0003] The relative amplitude processing technology in the prior art usually preserves the wideband signal energy of the seismic data after high-resolution processing to the narrowband energy before processing in the time domain. This way of amplitude preservation processing in the time domain will damage the true amplitude of the seismic data to some extent, and will damage the amplitude information and phase information including the original frequency band. Therefore, using the technical solution in the prior art cannot accurately restore the full-band amplitude information of the seismic data after high-resolution processing, thereby providing a good foundation for quantitative analysis of reservoir lithology and interpretation of geological phenomena.

[0004] At present, there is no effective solution to the above problems. SUMMARY

[0005] The embodiments of the present specification provide a mapping amplitude preservation method, device and equipment to solve the problem that the prior art cannot accurately restore the full-band amplitude information of the seismic data after high-resolution processing, thereby providing a good foundation for quantitative analysis of reservoir lithology and interpretation of geological phenomena.

[0006] The embodiments of the present specification provide a mapping amplitude preservation method, comprising: obtaining a first seismic signal; wherein the first seismic signal is an original seismic signal; obtaining a second seismic signal; wherein the second seismic signal is a seismic signal after high-resolution processing on the original seismic signal; performing time-frequency analysis on the first seismic signal to obtain an effective frequency band; establishing a mapping relationship between the first seismic signal and the second seismic signal in amplitude within the effective frequency band range; and performing amplitude preservation processing on the full-band signal of the second seismic signal according to the mapping relationship to obtain a seismic signal after amplitude preservation processing.

[0007] The embodiment of the present specification further provides a mapping amplitude preserving device, comprising: a first acquisition module configured to acquire a first seismic signal; wherein the first seismic signal is an original seismic signal; a second acquisition module configured to acquire a second seismic signal; wherein the second seismic signal is a processed seismic signal; a time-frequency analysis module configured to perform time-frequency analysis on the first seismic signal to obtain an effective frequency band; a building module configured to build a mapping relationship in amplitude between the first seismic signal and the second seismic signal within the effective frequency band; and an amplitude preserving processing module configured to perform amplitude preserving processing on a full frequency band signal of the second seismic signal according to the mapping relationship to obtain an amplitude preserved seismic signal.

[0008] The embodiment of the present specification further provides a mapping amplitude preserving device, comprising a processor and a memory for storing processor executable instructions, wherein the processor implements the steps of the mapping amplitude preserving method when executing the instructions.

[0009] The embodiment of the present specification further provides a computer readable storage medium having computer instructions stored thereon, wherein the instructions implement the steps of the mapping amplitude preserving method when executed.

[0010] The embodiment of the present specification provides a mapping amplitude preserving method, which can obtain an original seismic signal and a seismic signal processed by high resolution, and perform processing on the original seismic signal to obtain an effective frequency band of the original seismic signal. Further, a mapping relationship in amplitude can be built between the seismic signals before and after high resolution processing within the effective frequency band, and the mapping relationship can be applied to a full frequency band seismic signal of the seismic signal processed by high resolution to perform amplitude preserving processing on the full frequency band seismic signal, thereby obtaining an amplitude preserved seismic signal. Thus, the relative amplitude of the seismic signal processed by high resolution can be accurately restored, the high and low frequency energy can be effectively recovered, and more reservoir structure information can be obtained by using the amplitude preserved seismic signal. Moreover, the mapping amplitude preserving method has small amount of calculation data, is widely applicable, is not limited by geological environmental factors, and has good universality. BRIEF DESCRIPTION OF DRAWINGS

[0011] The accompanying drawings, which are included to provide a further understanding of the embodiments of the present specification and constitute a part of the present specification, do not limit the present specification. In the drawings:

[0012] Figure 1 is a step schematic diagram of the mapping amplitude preserving method provided by the embodiment of the present specification;

[0013] Figure 2 is a wavelet transform time-frequency decomposition schematic diagram provided by the embodiment of the present specification;

[0014] Figure 3 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0015] Figure 4 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0016] Figure 5 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0017] Figure 6 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0018] Figure 7 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0019] Figure 8 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0020] Figure 9 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0021] Figure 10 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0022] Figure 11 is a wavelet transform time-frequency decomposition schematic diagram provided according to an embodiment of the present specification;

[0023] Figure 12 is a single-shot original data seismic profile schematic diagram provided according to an embodiment of the present specification;

[0024] Figure 13 is a single-shot original data seismic profile schematic diagram provided according to an embodiment of the present specification;

[0025] Figure 14 is a single-shot seismic profile schematic diagram after high-resolution processing provided according to an embodiment of the present specification;

[0026] Figure 15 is a single-shot seismic profile schematic diagram after band-pass filtering of processed single-shot data provided according to an embodiment of the present specification;

[0027] Figure 16 is an amplitude-preserved processed single-channel seismic profile schematic diagram provided according to an embodiment of the present specification;

[0028] Figure 17 is a single-shot original data seismic profile schematic diagram provided according to an embodiment of the present specification;

[0029] Figure 18 is a schematic diagram of a single-shot seismic profile after band-pass filtering of raw data according to an embodiment of the present specification;

[0030] Figure 19 is a schematic diagram of a single-shot seismic profile after high-resolution processing according to an embodiment of the present specification;

[0031] Figure 20 is a schematic diagram of a single-shot seismic profile after band-pass filtering of processed single-shot data according to an embodiment of the present specification;

[0032] Figure 21 is a schematic diagram of a single-channel seismic profile after amplitude-preserved processing according to an embodiment of the present specification;

[0033] Figure 22 is a schematic diagram of the structure of an amplitude-preserved mapping device according to an embodiment of the present specification;

[0034] Figure 23 is a schematic diagram of the structure of an amplitude-preserved mapping device according to an embodiment of the present specification. DETAILED DESCRIPTION

[0035] The principles and spirits of the embodiments of the present specification will be described below with reference to a number of exemplary embodiments. It should be understood that these embodiments are given only to enable those skilled in the art to better understand and implement the embodiments of the present specification, and do not limit the scope of the embodiments of the present specification in any way. On the contrary, these embodiments are provided to make the disclosure of the embodiments of the present specification more thorough and complete, and to fully convey the scope of the disclosure to those skilled in the art.

[0036] Those skilled in the art know that the embodiments of the present specification can be implemented as a system, a device, a method or a computer program product. Therefore, the embodiments of the present specification can be specifically implemented as follows: complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0037] Although the processes described below include a number of operations that appear in a particular order, it should be clear that these processes can include more or fewer operations, and that the operations can be performed in sequence or in parallel (for example, using parallel processors or multi-threaded environments).

[0038] Referring to Figure 1 , the present embodiment can provide an amplitude-preserved mapping method. The amplitude-preserved mapping method described above can include the following steps.

[0039] S101: Obtain a first seismic signal; wherein the first seismic signal is a raw seismic signal.

[0040] In the embodiment, the first seismic signal can be obtained, wherein the first seismic signal can be an original seismic signal. The original seismic signal can be an original single-shot multi-channel seismic record collected by a seismic exploration instrument (or a seismic recording instrument), and it can be understood that the original seismic signal can also be other forms of seismic signals, and the specific form can be determined according to actual conditions, and the specification does not limit this.

[0041] In the embodiment, the first seismic signal can be obtained in the following manner: pulled from a preset database or received by a user input. It can be understood that the first seismic signal can also be obtained in other possible manners, for example, searched in a webpage according to certain search conditions, and the specific manner can be determined according to actual conditions, and the specification does not limit this.

[0042] S102: Obtain a second seismic signal; wherein the second seismic signal is a seismic signal processed by high-resolution processing on an original seismic signal.

[0043] In the embodiment, the second seismic signal can be obtained, wherein the second seismic signal can be a seismic signal processed by high-resolution processing on an original seismic signal. High-resolution processing is a key to further improve the reservoir resolution capability, and therefore, the original seismic signal can be processed by high-resolution processing. The high-resolution processing method can include HHT point spectrum whitening (Hilbert-Huang Transform point spectrum whitening), pre-stack seismic data noise ratio improvement processing, main frequency improvement of reflection signals and effective signal band widening, frequency division processing method, optimized iterative stacking method, signal direction constraint prediction denoising method, and the like. The high-resolution processing method is not limited to the above examples, and other changes can be made by those skilled in the art under the inspiration of the technical essence of the specification, as long as the functions and effects achieved are the same or similar to those of the specification, and they should be covered by the protection scope of the specification.

[0044] In the embodiment, the second seismic signal can be obtained in the following manner: pulled from a preset database or received by a user input. It can be understood that the second seismic signal can also be obtained in other possible manners, for example, searched in a webpage according to certain search conditions, and the specific manner can be determined according to actual conditions, and the specification does not limit this.

[0045] S103: Perform time-frequency analysis on the first seismic signal to obtain an effective frequency band.

[0046] In the embodiment, since the first seismic signal contains multiple frequency bands, in order to improve the effectiveness of the mapping relationship and reduce the data processing amount, time-frequency analysis can be performed on the first seismic signal to obtain an effective frequency band. The effective frequency band can be a frequency band range with consistent amplitudes from top to bottom, that is, a frequency band with uniform energy distribution.

[0047] In the embodiment, the time-frequency analysis includes a technology of studying a signal from time domain and frequency domain at the same time by using various time-frequency representations. The time-frequency analysis is neither observation of a one-dimensional signal (a real or complex function with one dimension as a one-dimensional real line) nor some transform (another function with one dimension as a one-dimensional real line is obtained from the original signal through the operation), but a two-dimensional signal is studied, that is, a function with two dimensions as a two-dimensional real plane is obtained from the signal through time-frequency transform.

[0048] In the embodiment, the form of time-frequency analysis can include Fourier transform, short-time Fourier transform, wavelet transform, empirical mode decomposition, Hilbert-Huang transform, and the like. The wavelet transform is preferred. The wavelet transform refers to representing a signal by using a finite-length or fast-decaying oscillation waveform called a “mother wavelet”. The waveform is scaled and shifted to match the input signal. The wavelet transform has better time resolution at high frequencies and better frequency resolution at low frequencies, which meets the resolution requirements of signal analysis at high and low frequencies. Of course, the form of time-frequency analysis is not limited to the above examples. Those skilled in the art can also make other changes under the inspiration of the technical essence of the embodiments of the present application. As long as the functions and effects achieved are the same or similar to those of the embodiments of the present application, they should be covered within the protection scope of the embodiments of the present application.

[0049] S104: Establish a mapping relationship between the first seismic signal and the second seismic signal in amplitude within the effective frequency band range.

[0050] In the embodiment, since the local frequency band component of the seismic signal has certain similarity with the characteristics of the full-band signal, the mapping relationship between the first seismic signal and the second seismic signal in amplitude can be established within the effective frequency band range.

[0051] In the embodiment, the mapping relationship can be obtained by comparing and analyzing the amplitudes of the first seismic signal within the effective frequency band range and the amplitudes of the second seismic signal within the effective frequency band range based on the mapping theory. That is, the mapping relationship can be used to represent the corresponding relationship in amplitude between the local original seismic signal within the effective frequency band and the whole seismic signal after high-resolution processing. The mapping relationship can be recorded in the form of data pairs or charts, and of course, it can be understood that the mapping relationship can also be recorded in any other possible form. The specific form can be determined according to actual conditions, and the embodiments of the present application do not limit this.

[0052] S105: Perform amplitude-preserved processing on the full-band signal of the second seismic signal according to the mapping relationship, to obtain an amplitude-preserved processed seismic signal.

[0053] In the embodiment, since high-resolution processing generally modifies the seismic waveform and amplitude to meet the fidelity requirement of the processed seismic signal for oil and gas reservoir exploration, the mapping relationship between the first seismic signal and the second seismic signal in amplitude in the effective frequency range can be applied to the full-band signal of the second seismic signal, so as to realize amplitude-preserved processing on other frequency components, and obtain an amplitude-preserved processed seismic signal.

[0054] In the embodiment, relative amplitude-preserved processing refers to that the amplitude attribute of a seismic signal remains unchanged or is proportional after one or more processing processes. For a forward model, the theoretical reflectivity of a reflection interface in the model is equal to or proportional to the reflectivity of the same interface after processing, that is, the incident wavelet of seismic data is basically consistent with the outgoing wavelet.

[0055] In the embodiment, it is assumed that the actually obtained original seismic signal is X0, the true amplitude seismic signal is X, the seismic signal obtained after high-resolution processing on the original seismic signal is Y, and the operator of the high-resolution processing is D. Then, Y=D(X0). In actual situations, the spectral bandwidth of Y is approximately equal to the bandwidth of X, but the amplitudes of the processed seismic signal Y and the true amplitude signal X are greatly different, and therefore, Y needs to be amplitude-preserved to X.

[0056] In the embodiment, on a seismic record, the amplitude value of a reflected wave is not only determined by the reflection coefficient of an interface, but also affected by the gain control of a seismic amplifier and the attenuation caused by divergence and absorption of the wave in a medium, and therefore, the original seismic record can not reflect the true amplitude value. When a parameter of a dynamic characteristic is needed, especially when oil is directly sought, the true amplitude value of a seismic signal must be obtained. True amplitude recovery can include two steps: the first is gain recovery, and the second is compensation of the amplitude value lost due to attenuation. Gain recovery of a digital record is to multiply the recorded signal value by the corresponding gain value, and amplitude compensation is to remove the amplitude by using an attenuation coefficient related to divergence and absorption. In some embodiments, the first seismic signal can also be a seismic signal obtained after true amplitude recovery on an original seismic signal, and the specific implementation can be determined according to actual situations, which is not limited in the embodiments of the present disclosure.

[0057] From the above description, it can be seen that the embodiments of the present specification achieve the following technical effects: the effective frequency band of the original seismic signal can be obtained by obtaining the original seismic signal and the seismic signal after high-resolution processing of the original seismic signal, and processing the original seismic signal. Further, the mapping relationship between the seismic signals before and after high-resolution processing can be established in the effective frequency band range, and the mapping relationship is applied to the full frequency band seismic signal of the seismic signal after high-resolution processing, to realize the amplitude-preserving processing of the full frequency band seismic signal, and obtain the seismic signal after amplitude-preserving processing. Thus, the relative amplitude of the seismic signal after high-resolution processing can be accurately restored, the high and low frequency energy can be effectively restored, and more reservoir structure information can be obtained by using the seismic signal after amplitude-preserving processing. And the mapping amplitude method has small calculation data, is widely used in various fields, is not limited by geological environmental factors, and has good universality.

[0058] In one embodiment, before the mapping relationship between the first seismic signal and the second seismic signal in the amplitude in the effective frequency band range is established, it can also include: respectively performing band-pass filtering on the first seismic signal and the second seismic signal in the effective frequency band range to obtain the filtered first seismic signal and the filtered second seismic signal.

[0059] In the present embodiment, the band-pass filter is a device that allows waves of a specific frequency band to pass while shielding other frequency bands. An ideal band-pass filter should have a completely flat passband, without amplification or attenuation in the passband, and all frequencies outside the passband are completely attenuated. The range of band-pass filtering can be the above-mentioned effective frequency band range, so that the signals of other frequency bands in the first seismic signal and the second seismic signal except the effective frequency band can be filtered out, and the data obtained after filtering has good amplitude preservation.

[0060] In one embodiment, corresponding to the establishment of the mapping relationship between the first seismic signal and the second seismic signal in the amplitude in the effective frequency band range, it can include: establishing a mapping relationship between the filtered first seismic signal and the filtered second seismic signal in the amplitude in the effective frequency band range.

[0061] In the present embodiment, it is assumed that the actually obtained original seismic signal is X0, the true amplitude seismic signal is X, the seismic signal obtained after high-resolution processing of the original seismic signal is Y, and the high-resolution processing operator is D. Then Y=D(X0). Assuming that the band-pass filtering operator is B, the mapping relationship established is P, and the envelope operator is H, the following can be obtained: The above formula satisfies: P(B(Y))=B(X0).

[0062] In the embodiment, the mapping relationship P can be further analyzed from the perspective of an operator, and the following operator relationship can be obtained:

[0063] Y = D(X0)

[0064] X = P(Y)

[0065] The seismic signal in the effective frequency band has the following relationship:

[0066] B(X) = B(X0)

[0067] Therefore, the following can be obtained:

[0068] B(P(Y)) = B(X0) = B(D -1 (Y))

[0069] When the mapping relationship operator P is an approximate linear operator, the order of the band-pass filtering operator B can be exchanged, and the following can be obtained:

[0070] B(X0) = D -1 (B(Y)) = P(B(Y))

[0071] It can be seen that the mapping relationship in the full frequency band can be approximately regarded as the mapping relationship in the effective frequency band range, and the mapping relationship in the effective frequency band range can also be used for the amplitude-preserving processing of the full frequency band signal. Therefore, the mapping relationship P can be applied to the full frequency band signal Y to realize the amplitude-preserving processing of Y.

[0072] In one embodiment, after the amplitude-preserving processing of the full frequency band signal of the second seismic signal according to the mapping relationship, the true amplitude seismic signal of the first seismic signal can be obtained, and the error between the amplitude-preserving processed seismic signal and the true amplitude seismic signal can be further determined. If the error is less than a preset threshold, the amplitude-preserving processing process can be ended.

[0073] In the embodiment, in order to ensure the accuracy of the amplitude-preserving processing result, if the error is greater than or equal to the preset threshold, the above effective frequency band can be adjusted, the mapping relationship can be re-determined based on the adjusted effective frequency band, and the amplitude-preserving processing can be performed again until the error between the amplitude-preserving processed seismic signal and the true amplitude seismic signal is less than the above preset threshold.

[0074] In the embodiment, on the seismic record, the amplitude value of the reflected wave is not only determined by the reflection coefficient of the interface, but also affected by the gain control of the seismic amplifier and the attenuation of the wave propagation in the medium, so the original seismic record may not reflect the true amplitude value. When the parameters of the dynamic characteristics are needed, especially when the oil is directly found, the true amplitude value of the seismic signal must be obtained. True amplitude recovery can include two steps: the first is gain recovery; the second is to compensate the amplitude value lost due to attenuation. The gain recovery of the digital record is to multiply the recorded signal value by the corresponding gain value, and the amplitude compensation is to remove the amplitude with the attenuation coefficient related to divergence and absorption.

[0075] In the embodiment, the true amplitude seismic signal X can represent the true amplitude value of the first seismic signal, and the amplitude preserving processing result should satisfy: ||P(Y)-X||2<ε, where ε is a preset threshold, the mapping relationship P can realize amplitude preserving in the error accuracy allowed range for the full frequency band signal Y. If the actual amplitude preserving processing result is greatly different from the expected processing result and does not meet the accuracy requirement, the effective frequency band range can be selected again to adjust the mapping relationship P between the data.

[0076] In one embodiment, the first seismic signal, the second seismic signal, the amplitude preserving processed seismic signal and the like can be displayed in trace normalization and global display, so as to better observe the amplitude of the whole data and the relative size relationship between the data of the traces.

[0077] In one scene example, the time-frequency analysis on the original seismic signal can obtain the wavelet transform time-frequency decomposition diagram as shown in Figures 2 to 11 , wherein Figures 2 to 11 The left half part is the original single shot record, and the right graph is the frequency division profile, Figures 2 to 11 The corresponding main frequency and frequency band range are: 6.36hz, 5hz-7hz; 9.00hz, 7hz-10hz; 12.73hz, 10hz-14hz; 18.00hz, 14hz-20hz; 25.46hz, 20hz-30hz; 36hz, 30hz-42hz; 50.92hz, 42hz-60hz; 72.00hz, 60hz-84hz; 101.92hz, 84hz-120hz; 144.00hz, 120hz-170hz. Figures 2 to 11 The horizontal coordinate is the trace number (TRACENO), and the vertical coordinate is the travel time (Time / ms). The effective frequency band of the original seismic signal can be selected according to Figures 2 to 11 The amplitude of the frequency division data with the frequency band of 20hz (hertz) to 30hz is stable from top to bottom, and the amplitude preserving is good, so the effective frequency band can be selected as 20hz to 30hz.

[0078] In one scenario example, Figure 12 is a single-shot raw data seismic profile, Figure 13 is a single-shot raw data seismic profile after band-pass filtering, according to Figure 12 and Figure 13 It can be seen that the amplitudes of the filtered data are consistent from top to bottom, and the effective frequency band is accurately selected. Figure 14 is a single-shot seismic profile after high-resolution processing, Figure 15 is a single-shot seismic profile after band-pass filtering on the processed single-shot data, which can establish a mapping relationship between the amplitudes of the data in Figure 13 and Figure 15 and apply the mapping relationship to the full-band data, thereby realizing amplitude-preserving processing of the full-band signal. Figure 16 is a single-channel seismic profile after amplitude-preserving processing, compared with Figure 12 the original data seismic profile, it can be seen that the energy of different frequency components of the signal is well recovered and maintained, and the true amplitude is recovered to the maximum extent. Figures 12 to 16 The seismic profiles are all in the normalized display mode, in which the normalized display mode can be to divide each channel of the seismic data by its maximum value, and the maximum value of each channel is 1. In the normalized display mode, Figures 12 to 16 the horizontal coordinate is the channel number (TRACENO), and the vertical coordinate is the travel time (Time / ms).

[0079] In one scenario example, Figure 17 is a single-shot raw data seismic profile, Figure 18 is a single-shot raw data seismic profile after band-pass filtering, according to Figure 17 and Figure 18 It can be seen that the amplitudes of the filtered data are consistent from top to bottom, and the effective frequency band is accurately selected. Figure 19 is a single-shot seismic profile after high-resolution processing, Figure 20 is a single-shot seismic profile after band-pass filtering on the processed single-shot data, which can establish a mapping relationship between the amplitudes of the data in Figure 18 and Figure 20 and apply the mapping relationship to the full-band data, thereby realizing amplitude-preserving processing of the full-band signal. Figure 21 is a single-channel seismic profile after amplitude-preserving processing, compared with Figure 17 the original data seismic profile, it can be seen that the high and low frequency energy is obviously recovered, the amplitude spatial consistency is improved, and the relative amplitude relationship is well maintained. Figures 17 to 21 The seismic profiles are all in the globalized display mode, in which the globalized display mode is to find the maximum value in all channels of the seismic data, and the displayed data is the single-shot data divided by the maximum value, and the relative size relationship between the channels is maintained. In the globalized display mode, Figures 17 to 21The horizontal coordinate is trace number (TRACENO) and the vertical coordinate is travel time (Time / ms).

[0080] Based on the same inventive concept, the embodiments of the present specification also provide a mapping amplitude preserving device, as follows. Since the mapping amplitude preserving device solves problems by similar principles to the mapping amplitude preserving method, the implementation of the mapping amplitude preserving device can be referred to the implementation of the mapping amplitude preserving method, and the repeated parts will not be described herein. The term "unit" or "module" used below can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware or a combination of software and hardware is also possible and is contemplated. Figure 22 is a structural block diagram of the mapping amplitude preserving device of the embodiments of the present specification, as shown in Figure 22 The structure can include a first obtaining module 201, a second obtaining module 202, a time-frequency analysis module 203, a establishing module 204, an amplitude preserving processing module 205, and the structure will be described below.

[0081] The first obtaining module 201 can be used to obtain a first seismic signal; wherein the first seismic signal is an original seismic signal;

[0082] The second obtaining module 202 can be used to obtain a second seismic signal; wherein the second seismic signal is a processed seismic signal;

[0083] The time-frequency analysis module 203 can be used to perform time-frequency analysis on the first seismic signal to obtain an effective frequency band;

[0084] The establishing module 204 can be used to establish a mapping relationship between the first seismic signal and the second seismic signal in amplitude within the effective frequency band;

[0085] The amplitude preserving processing module 205 can be used to perform amplitude preserving processing on the full frequency band signal of the second seismic signal according to the mapping relationship to obtain an amplitude preserved seismic signal.

[0086] The embodiments of the present specification also provide an electronic device, which can be specifically referred to Figure 23The electronic device shown in the mapping amplitude preserving method provided by the embodiments of the present application includes an input device 31, a processor 32 and a memory 33. The input device 31 is configured to input a first seismic signal and a second seismic signal. The first seismic signal is an original seismic signal, and the second seismic signal is a seismic signal processed by high resolution on the original seismic signal. The processor 32 is configured to perform time-frequency analysis on the first seismic signal to obtain an effective frequency band, establish a mapping relationship between the first seismic signal and the second seismic signal in amplitude within the effective frequency band, and perform amplitude preserving processing on a full frequency band signal of the second seismic signal according to the mapping relationship to obtain a seismic signal processed by amplitude preserving. The memory 33 is configured to store the seismic signal processed by amplitude preserving and the like.

[0087] In the embodiment, the input device can be one of the main devices for information exchange between the user and the computer system. The input device can include a keyboard, a mouse, a camera, a scanner, a light pen, a handwriting input board, a voice input device, etc. The input device is used to input raw data and programs for processing the data into the computer. The input device can also obtain data transmitted by other modules, units and devices. The processor can be implemented in any appropriate manner. For example, the processor can take the form of a microprocessor or a processor and a computer readable medium storing computer readable program code (such as software or firmware) executable by the (micro) processor, logic gates, switches, application specific integrated circuits (ASIC), programmable logic controllers and embedded microcontrollers, etc. The memory can be a memory device used to save information in modern information technology. The memory can include multiple levels, and in a digital system, as long as it can save binary data, it can be a memory; in an integrated circuit, a circuit without a physical form and with a storage function is also called a memory, such as RAM, FIFO, etc.; in a system, a storage device with a physical form is also called a memory, such as a memory stick, a TF card, etc.

[0088] In the embodiment, the functions and effects of the electronic device can be explained in comparison with other embodiments, and will not be described here.

[0089] The embodiment of the present specification also provides a computer storage medium based on a mapping amplitude-preserving method, which stores computer program instructions, and when the computer program instructions are executed, the following can be achieved: obtaining a first seismic signal; wherein the first seismic signal is an original seismic signal; obtaining a second seismic signal; wherein the second seismic signal is a seismic signal after high-resolution processing of the original seismic signal; performing time-frequency analysis on the first seismic signal to obtain an effective frequency band; establishing a mapping relationship between the first seismic signal and the second seismic signal in amplitude within the effective frequency band; and performing amplitude-preserving processing on the full-band signal of the second seismic signal according to the mapping relationship to obtain an amplitude-preserving processed seismic signal.

[0090] In the embodiment, the storage medium includes but is not limited to a random access memory (RAM), a read-only memory (ROM), a cache, a hard disk drive (HDD), or a memory card. The storage medium can be used to store computer program instructions. The network communication unit can be an interface set according to a standard specified by a communication protocol, and is used for network connection communication.

[0091] In the embodiment, the functions and effects of the program instructions stored in the computer storage medium can be explained in comparison with other embodiments, and will not be described here.

[0092] Obviously, those skilled in the art should understand that each module or each step of the above-mentioned embodiment of the present specification can be realized by a general computing device, which can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, and optionally, they can be realized by program codes executable by a computing device, so that they can be stored in a storage device and executed by a computing device, and in some cases, the steps shown or described can be executed in an order different from here, or they can be manufactured into each integrated circuit module, or multiple modules or steps among them can be manufactured into a single integrated circuit module. Therefore, the embodiment of the present specification is not limited to any specific combination of hardware and software.

[0093] Although the method operations are described in a particular, sequential order, one of ordinary skill in the art will understand that the method operations can be performed in other sequences other than the described, or that the method operations can be performed in parallel, or that the method operations can be performed at different times, or in different systems. The method operations can be implemented using any combination of hardware, firmware, or software methods or entities designed for hardware implementation, including one or more application specific integrated circuits (ASICs), digital signal processors (DSPs), general purpose processors, or processors communicating via a system bus.

[0094] It is to be understood that the above description is intended to be illustrative, and not restrictive. Many other embodiments will be apparent to those of skill in the art upon reading the above description. The scope of embodiments should, therefore, be determined not with reference to the above description, but should instead be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. The disclosures of patents and patent documents herein are incorporated by reference in their entirety. In the event of any inconsistency between the description provided herein and the disclosure in the incorporated references, the description provided herein will control.

[0095] The above description is intended to be illustrative, and not restrictive. Many other embodiments will be apparent to those of skill in the art upon reading the above description. The scope of embodiments should, therefore, be determined not with reference to the above description, but should instead be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. The disclosures of patents and patent documents herein are incorporated by reference in their entirety. In the event of any inconsistency between the description provided herein and the disclosure in the incorporated references, the description provided herein will control.

Claims

1. A mapping amplitude preserving method characterized by, The method comprises the following steps: obtaining a first seismic signal; wherein the first seismic signal is an original seismic signal; obtaining a second seismic signal; wherein the second seismic signal is a seismic signal processed by high-resolution processing on the original seismic signal; performing time-frequency analysis on the first seismic signal to obtain an effective frequency band; wherein the time-frequency analysis on the first seismic signal is performed by using wavelet transform to obtain the effective frequency band; the first seismic signal comprises a plurality of frequency bands; the effective frequency band is a frequency band with the most consistency in amplitude from top to bottom and the most uniform energy distribution among the plurality of frequency bands; establishing a mapping relationship in amplitude between the first seismic signal and the second seismic signal within the effective frequency band range; performing amplitude-preserving processing on a full-band signal of the second seismic signal according to the mapping relationship to obtain an amplitude-preserved seismic signal; obtaining a true-amplitude seismic signal of the first seismic signal; determining an error between the amplitude-preserved seismic signal and the true-amplitude seismic signal; in the case that the error is less than a preset threshold, ending the amplitude-preserving processing procedure; in the case that the error is greater than or equal to the preset threshold, adjusting the effective frequency band until the error between the amplitude-preserved seismic signal and the true-amplitude seismic signal is less than the preset threshold.

2. The method of claim 1, wherein, Before the step of establishing the mapping relationship in amplitude between the first seismic signal and the second seismic signal within the effective frequency band range, the method further comprises the following steps: respectively performing band-pass filtering on the first seismic signal and the second seismic signal to obtain a filtered first seismic signal and a filtered second seismic signal; wherein the band-pass filtering range is the effective frequency band range.

3. The method of claim 2, wherein, The step of establishing the mapping relationship in amplitude between the first seismic signal and the second seismic signal within the effective frequency band range comprises the step of establishing the mapping relationship in amplitude between the filtered first seismic signal and the filtered second seismic signal within the effective frequency band range.

4. The method of claim 1, wherein, After the step of performing amplitude-preserving processing on the full-band signal of the second seismic signal according to the mapping relationship, the method further comprises the following step:

5. An amplitude preserving mapping device, characterized by performing trace normalization display and global display on the amplitude-preserved seismic signal. The method comprises the following steps: a first obtaining module is configured to obtain a first seismic signal; wherein the first seismic signal is an original seismic signal; a second obtaining module is configured to obtain a second seismic signal; wherein the second seismic signal is a processed seismic signal; a time-frequency analysis module is configured to perform time-frequency analysis on the first seismic signal to obtain an effective frequency band; wherein the time-frequency analysis on the first seismic signal is performed by using wavelet transform to obtain the effective frequency band; the first seismic signal comprises a plurality of frequency bands; the effective frequency band is a frequency band with the most consistency in amplitude from top to bottom and the most uniform energy distribution among the plurality of frequency bands; an establishing module is configured to establish a mapping relationship in amplitude between the first seismic signal and the second seismic signal within the effective frequency band range; an amplitude-preserving processing module is configured to perform amplitude-preserving processing on a full-band signal of the second seismic signal according to the mapping relationship to obtain an amplitude-preserved seismic signal; the mapping amplitude-preserving device is further configured to: taking a true amplitude seismic signal of the first seismic signal; determining an error between the amplitude-preserved seismic signal and the true amplitude seismic signal; ending the amplitude-preservation processing procedure if the error is less than a preset threshold; adjusting the effective frequency band until the error between the amplitude-preserved seismic signal and the true amplitude seismic signal is less than the preset threshold if the error is greater than or equal to the preset threshold.

6. A mapping preserving device, characterized by A computer program product comprising a computer readable medium having stored thereon computer instructions, the instructions, when executed by a computer, implementing the steps of any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, A computer program product comprising a computer readable medium having stored thereon computer instructions, the instructions, when executed by a computer, implementing the steps of any one of claims 1 to 4.

Images