Pre-stack inversion model fusion construction method and system based on synchronous extrusion wavelet transform

Abstract

The invention discloses a pre-stack inversion model fusion construction method and system based on synchronous extrusion wavelet transform. The method comprises the steps of obtaining a low-frequency speed model within a seismic data time range and a high-frequency speed model at the top and above of a fracture-cavity body; performing time-frequency analysis on the high-frequency model and the low-frequency model by utilizing synchronous extrusion wavelet transform to obtain a synchronous extrusion wavelet transform coefficient; analyzing an overlapping condition based on a synchronous extrusion wavelet transform coefficient to obtain a frequency spectrum at a certain moment; based on the frequency spectrum of a certain moment, performing frequency fusion processing on all moments of a certain channel of data to obtain a relatively wide-band time-frequency spectrum of the channel of data; and obtaining a wide-band time-frequency spectrum for all seismic traces to obtain relative broadband model data. According to the method, time-frequency domain fusion based on synchronous extrusion wavelet transform is carried out through velocity fields of two backgrounds, so that an initial inversion model for revealing deep information is constructed, and the problems of insufficient fracture-cavity body and lower drilling data and incomplete initial model information are effectively solved.

Description

Technical Field

[0001] The invention belongs to the technical field of pre-stack inversion initial model construction for ultra-deep carbonate reservoirs, and particularly relates to a method and system for fusion construction of a pre-stack inversion model based on synchronous squeezing wavelet transform. Background Art

[0002] In the process of advancing into deep and ultra-deep carbonate rocks, more and more oil and gas have been discovered, making it the most active region for ultra-deep oil and gas exploration and development in the world. The Ordovician fracture-cave reservoir in the Tarim Basin is the most special large-scale ultra-deep carbonate reservoir discovered in China. The fault-controlled karst fracture-cave reservoir is extremely developed along the fault trend, and the fault-controlled karst elements are very representative. Most wells in the Tarim Basin encounter karst when entering the Ordovician strata, and a large number of emptying and leakage occur. According to the drilling design principles, drilling is completed to target point B in principle; if the target layer section has complex situations such as emptying and leakage in advance, plugging measures should be considered to drill through the fracture zone as much as possible; if the drilling is complicated during the drilling process, resulting in the possibility of uncontrollable well control risks, the drilling can be completed in advance after the integrated discussion and decision-making of engineering geology. The logging instrument cannot enter the cave for measurement along the wall, and deeper information cannot be obtained, resulting in scarce drilling data for deep carbonate rocks. However, for the detailed characterization of oil and gas content based on prestack inversion parameters, the lack of internal information of caves makes the construction of the initial inversion model incomplete, which is not conducive to the detailed characterization of ultra-deep carbonate reservoirs. Summary of the invention

[0003] The purpose of the present invention is to provide a method and system for fusion construction of prestack inversion model based on synchronous squeezing wavelet transform, so as to solve the problem that the initial inversion model construction is incomplete due to the lack of internal information of the cave.

[0004] To achieve the above object, the present invention adopts the following technical solutions:

[0005] In a first aspect, the present invention provides a method for constructing a prestack inversion model fusion based on synchronous squeezing wavelet transform, comprising:

[0006] Collect the pre-stack migration velocity field seismic data and the information of the top of the fracture-cavity body, and obtain the low-frequency velocity model within the time range of the seismic data and the high-frequency velocity model at and above the top of the fracture-cavity body;

[0007] The synchronous squeezing wavelet transform is used to perform time-frequency analysis on the high-frequency model and the low-frequency model, and the synchronous squeezing wavelet transform coefficients are obtained;

[0008] Based on the analysis of overlapping conditions of synchronously squeezed wavelet transform coefficients, the spectrum at a certain moment is obtained;

[0009] Based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain a relatively wide-segment time-frequency spectrum of the data channel;

[0010] The wide-band time-frequency spectra are acquired for all seismic traces to obtain relatively broadband model data.

[0011] Optionally, collect pre-stack migration velocity field seismic data and information on the top of the fracture-cavity body to obtain a low-frequency velocity model within the time range of the seismic data and a high-frequency velocity model at and above the top of the fracture-cavity body:

[0012] The prestack migration velocity field is constructed based on seismic data and is used to establish a low-frequency velocity model within the time range of seismic data. The information on the top of the fracture-cavity body is obtained based on the logging data, and the inverted high-frequency velocity model is constructed through the information on the top of the cave body.

[0013] Optionally, the synchronous squeezing wavelet transform is used to perform time-frequency analysis on the high-frequency model and the low-frequency model to obtain the synchronous squeezing wavelet transform coefficients:

[0014] Using Synchrosqueezing Wavelet Transform to Analyze High Frequency Model V high (x, t) and the low-frequency model V low (x, t) to perform time-frequency analysis. high (x i , t) and V low (x i , t) are obtained high (x i , t, f) and F low (x i , t, f); for data V high (x i , t), firstly, we use formula (1) to perform wavelet transform to obtain WF high (x i , a, b):

[0015]

[0016] in, It is the complex conjugate of the mother wavelet. a and b are the scaling factor and translation factor, respectively. The scaling factor is related to the frequency, and the translation factor is related to the time. On this basis, the instantaneous frequency ω is calculated. s (x i , a, b):

[0017]

[0018] After obtaining the instantaneous frequency, the time-frequency spectrum of the wavelet transform is squeezed in the frequency direction to obtain the synchronous squeezed wavelet transform coefficient F high (x i , t, f), that is:

[0019]

[0020] For data V low (x, t) can also be obtained according to formulas (1)-(3) low (x i , t, f).

[0021] Optionally, the overlap is analyzed based on the synchronously squeezed wavelet transform coefficients to obtain the spectrum at a certain moment:

[0022] For the i-th data at a certain moment F high (x i , t i , f) and F low (x i , t i , f) perform time-frequency analysis; F high (x i , t i , f) and F low (x i , t i , f) There are two overlapping situations:

[0023] Where two spectra intersect, an inverted Hamming window is added at the intersection to smooth the value;

[0024] Where the two spectra do not intersect, the two spectra are directly added.

[0025] Optionally, in the first case, where the two spectra intersect, an inverted Hamming window is added to the intersection for smoothing and obtaining F all (x i , t i , f) is the t of the i-th data i The spectrum at the moment is:

[0026]

[0027] Optionally, where the two spectra do not cross, the two spectra are directly added as F all (x i , t i , f), that is, F all (x i , t i , f)=F 1ow (x i , t i , f)+F high (x i , t i , f).

[0028] Optionally, based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain a relatively wide-segment time-frequency spectrum of the data channel:

[0029] Frequency fusion processing is performed on all moments of the i-th channel data, and the relatively wide-band time-frequency spectrum F of the i-th channel data can be obtained. all (x i , t, f), and then use the synchronous squeezing wavelet inverse transform to transform F all (x i , t, f) is transformed into V all (x i , t), that is:

[0030]

[0031] That is, relatively wide-band seismic trace data.

[0032] In a second aspect, the present invention provides a prestack inversion model fusion construction system based on synchronous squeezing wavelet transform, comprising:

[0033] A data acquisition module is used to collect pre-stack migration velocity field seismic data and information on the top of the fracture-cavity body, and obtain a low-frequency velocity model within the time range of the seismic data and a high-frequency velocity model at and above the top of the fracture-cavity body;

[0034] A synchronous squeezing wavelet transform coefficient acquisition module is used to perform time-frequency analysis on the high-frequency model and the low-frequency model using synchronous squeezing wavelet transform to obtain synchronous squeezing wavelet transform coefficients;

[0035] The frequency fusion processing module is used to analyze the overlapping situation based on the synchronous squeezing wavelet transform coefficients to obtain the spectrum at a certain moment; based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain the relatively wide-segment time-frequency spectrum of the data channel; wide-segment time-frequency spectrum is obtained for all seismic channels to obtain relatively wide-band model data.

[0036] In a third aspect, the present invention provides a computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the steps of a method for constructing a prestack inversion model fusion based on synchronous squeezing wavelet transform are implemented.

[0037] In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of a method for constructing a prestack inversion model fusion based on synchronous squeezing wavelet transform.

[0038] Compared with the prior art, the present invention has the following technical effects:

[0039] The prestack migration velocity field is constructed based on seismic data, and can establish a low-frequency background velocity field within the time range of seismic data. However, when drilling into the Ordovician fracture-cavity body in the Tarim Basin, well logging only obtains information on the top of the fracture-cavity body, and this information can be used to construct an inverted high-frequency velocity field. The present invention uses the time-frequency domain fusion of the two background velocity fields based on synchronous extrusion wavelet transform to construct an inversion initial model that reveals deep information, effectively solving the problem of insufficient fracture-cavity body and lower drilling data and incomplete initial model information.

[0040] The present invention can supplement the fracture-cavity information when the well logging curve lacks fracture-cavity information and the inversion model requires relevant fracture-cavity information, thereby constructing a relatively complete initial inversion model.

[0041] The frequency fusion based on time-frequency analysis is realized by using synchronous squeezing wavelet transform. Compared with the conventional frequency fusion based on Fourier transform, its accuracy is higher and the precision is also improved.

[0042] The inverted Henning window is used to smoothly fuse the frequency overlapping parts to avoid frequency mutation while maintaining the effective frequency information. Compared with the conventional rectangular window, its spectrum is smoother, and the result obtained by mapping it to the time domain is more stable, which is also conducive to the stability of subsequent inversion results. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 It is a flow chart of the present invention.

[0044] Figure 2 This is a schematic diagram of frequency domain merging when two models overlap in the frequency domain.

[0045] Figure 3 Schematic diagram of frequency domain merging when two models do not overlap in the frequency domain.

[0046] Figure 4 This is a schematic diagram of the technical principle.

[0047] Figure 5 It is a low-frequency model for the entire time period.

[0048] Figure 6 It is a high-frequency model at and above the top of the fracture-cavity body.

[0049] Figure 7 It is a relatively broadband model.

[0050] Figure 8 It is the high-frequency model inversion result at and above the top of the fracture-cavity body in the embodiment.

[0051] Fig. 9 It is the inversion result of the relatively broadband model in the embodiment. DETAILED DESCRIPTION

[0052] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0053] In the description of the present invention, it should be understood that the terms “include” and “comprises” indicate the presence of described features, wholes, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or collections thereof.

[0054] It should also be understood that the terms used in the present specification are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the present specification and the appended claims, the singular forms "a", "an" and "the" are intended to include plural forms unless the context clearly indicates otherwise.

[0055] It should be further understood that the term "and / or" used in the present specification and the appended claims refers to any combination of one or more of the associated listed items and all possible combinations, and includes these combinations. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in the present invention generally indicates that the associated objects are in an "or" relationship.

[0056] It should be understood that, although the terms first, second, third, etc. may be used to describe preset ranges, etc. in the embodiments of the present invention, these preset ranges should not be limited to these terms. These terms are only used to distinguish preset ranges from each other. For example, without departing from the scope of the embodiments of the present invention, the first preset range may also be referred to as the second preset range, and similarly, the second preset range may also be referred to as the first preset range.

[0057] The word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrases "if it is determined" or "if (stated condition or event) is detected" may be interpreted as "when it is determined" or "in response to determining" or "when detecting (stated condition or event)" or "in response to detecting (stated condition or event)", depending on the context.

[0058] Various structural schematic diagrams of the embodiments disclosed in the present invention are shown in the accompanying drawings. These figures are not drawn to scale, and some details are magnified and some details may be omitted for the purpose of clear expression. The shapes of various regions and layers shown in the figures and the relative sizes and positional relationships therebetween are only exemplary, and may deviate in practice due to manufacturing tolerances or technical limitations, and those skilled in the art may additionally design regions / layers with different shapes, sizes, and relative positions according to actual needs.

[0059] The present invention is further described below in conjunction with the accompanying drawings:

[0060] See also Figures 1 to 9 , based on synchronous squeezing wavelet transform prestack inversion model fusion construction method, including:

[0061] Collect the pre-stack migration velocity field seismic data and the information of the top of the fracture-cavity body, and obtain the low-frequency velocity model within the time range of the seismic data and the high-frequency velocity model at and above the top of the fracture-cavity body;

[0062] The synchronous squeezing wavelet transform is used to perform time-frequency analysis on the high-frequency model and the low-frequency model, and the synchronous squeezing wavelet transform coefficients are obtained;

[0063] Based on the analysis of overlapping conditions of synchronously squeezed wavelet transform coefficients, the spectrum at a certain moment is obtained;

[0064] Based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain a relatively wide-segment time-frequency spectrum of the data channel;

[0065] The wide-band time-frequency spectra are acquired for all seismic traces to obtain relatively broadband model data.

[0066] The prestack migration velocity field is constructed based on seismic data, and can establish a low-frequency background velocity field within the time range of seismic data. However, when drilling into the Ordovician fracture-cavity body in the Tarim Basin, well logging only obtains information on the top of the fracture-cavity body, and this information can be used to construct an inverted high-frequency velocity field. The present invention uses the time-frequency domain fusion of the two background velocity fields based on synchronous extrusion wavelet transform to construct an inversion initial model that reveals deep information, effectively solving the problem of insufficient fracture-cavity body and lower drilling data and incomplete initial model information.

[0067] Specific:

[0068] A frequency fusion technology based on time-frequency analysis is proposed. The low-frequency model and the high-frequency model are analyzed by synchronous squeezing wavelet transform to obtain different time-frequency spectra of the same channel. Since there is some overlap between the frequencies, it is proposed to use the inverted Hamming window for numerical fusion.

[0069] The details are as follows:

[0070] Using Synchrosqueezing Wavelet Transform to Analyze High Frequency Model V high (x, t) and the low-frequency model V low (x, t) to perform time-frequency analysis. high (x i , t) and V low (x i , t) can get F high (x i , t, f) and F low (x i , t, f). For data V high (x i , t), firstly, we use formula (1) to perform wavelet transform to obtain WF high (x i , a, b):

[0071]

[0072] in, It is the complex conjugate of the mother wavelet. a and b are the scaling factor and translation factor, respectively. The scaling factor is related to frequency, and the translation factor is related to time. On this basis, the instantaneous frequency ω is calculated. s (x i , a, b):

[0073]

[0074] After obtaining the instantaneous frequency, the time-frequency spectrum of the wavelet transform is squeezed in the frequency direction to obtain the synchronous squeezed wavelet transform coefficient F high (x i , t, f), that is:

[0075]

[0076] For data V low (x, t) can also be obtained according to formulas (1)-(3) low (x i , t, f).

[0077] For the i-th data at a certain moment F high (x i , t i , f) and F low (x i , t i , f) for detailed analysis. high (x i , t i , f) and F low (x i , t i, f) may have two overlapping situations. The first one is as follows Figure 1 The two blue lines represent the shapes of the two spectra. Where the two spectra intersect, an inverted Hamming window is added to smooth the values ​​at the intersection to avoid frequency mutations, and the result is F all (x i , t i , f) is the t of the i-th data i The spectrum at the moment is:

[0078]

[0079] The second case is that the two spectra are slightly different, such as Figure 2 As shown, the two spectra can be directly added together to obtain F all (x i , t i , f), that is, F all (x i , t i , f)=F low (x i , t i , f)+F high (x i , t i , f)

[0080] Frequency fusion processing is performed on all moments of the i-th channel data, and the relatively wide-band time-frequency spectrum F of the i-th channel data can be obtained. all (x i , t, f). Then use the synchronous squeezing wavelet inverse transform to transform F all (x i , t, f) into V all (x i , t), that is:

[0081]

[0082] That is, relatively wide-band seismic trace data.

[0083] By processing all seismic traces as above, relatively broadband model data can be obtained.

[0084] The present invention can supplement the fracture-cavity information when the well logging curve lacks fracture-cavity information and the inversion model requires relevant fracture-cavity information, thereby constructing a relatively complete initial inversion model.

[0085] The frequency fusion based on time-frequency analysis is realized by using synchronous squeezing wavelet transform. Compared with the conventional frequency fusion based on Fourier transform, its accuracy is higher and the precision is also improved.

[0086] The inverted Henning window is used to smoothly fuse the frequency overlapping parts to avoid frequency mutation while maintaining the effective frequency information. Compared with the conventional rectangular window, its spectrum is smoother, and the result obtained by mapping it to the time domain is more stable, which is also conducive to the stability of subsequent inversion results.

[0087] Example:

[0088] Taking Yueman 22 well in Yueman West Work Area as an example, firstly, the low-frequency velocity model of the whole time period is obtained by migrating the velocity. Figure 5 As shown, its frequency range is 1-12 Hz. Then the high-frequency velocity model is obtained by interpolation through logging data. Figure 6 As shown in , its frequency range is 1-70Hz. Since the two models have overlapping parts in the frequency domain, the following Figure 2 The frequency domain merging strategy is shown.

[0089] The synchronous squeezing wavelet transform is performed on the low-frequency velocity model and the high-frequency velocity model of the entire time period respectively, and the frequency is merged by formula 4, where f1 is selected as 8Hz and f2 is selected as 16Hz. The synchronous squeezing wavelet transform is inversely transformed on the frequency-merged data body to obtain the relatively broadband model of the patented method as shown in FIG. Figure 7 shown.

[0090] Figure 8 In order to use the relatively broadband model obtained by fusion to perform prestack inversion of P-wave impedance,

[0091] Fig. 9 To use the high-frequency model at and above the top of the fracture-cavity body for pre-stack inversion of the longitudinal wave impedance.

[0092] Compared with the inversion results obtained by the high-frequency model, the inversion results of the broadband model have more complete fracture-cavity morphology and richer deep information.

[0093] In yet another embodiment of the present invention, a system for constructing a prestack inversion model based on synchronous squeezing wavelet transform fusion is provided, which can be used to implement the above-mentioned method for constructing a prestack inversion model based on synchronous squeezing wavelet transform fusion. Specifically, the system comprises:

[0094] A data acquisition module is used to collect pre-stack migration velocity field seismic data and information on the top of the fracture-cavity body, and obtain a low-frequency velocity model within the time range of the seismic data and a high-frequency velocity model at and above the top of the fracture-cavity body;

[0095] A synchronous squeezing wavelet transform coefficient acquisition module is used to perform time-frequency analysis on the high-frequency model and the low-frequency model using synchronous squeezing wavelet transform to obtain synchronous squeezing wavelet transform coefficients;

[0096] The frequency fusion processing module is used to analyze the overlapping situation based on the synchronous squeezing wavelet transform coefficients to obtain the spectrum at a certain moment; based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain the relatively wide-segment time-frequency spectrum of the data channel; wide-segment time-frequency spectrum is obtained for all seismic channels to obtain relatively wide-band model data.

[0097] The present invention performs time-frequency domain fusion based on synchronous squeezing wavelet transform through two background velocity fields, thereby constructing an inversion initial model that reveals deep information, effectively solving the problems of insufficient fracture-cavity and lower drilling data and incomplete initial model information.

[0098] The division of modules in the embodiments of the present invention is schematic and is only a logical function division. There may be other division methods in actual implementation. In addition, each functional module in each embodiment of the present invention may be integrated into one processor, or may exist physically separately, or two or more modules may be integrated into one module. The above-mentioned integrated modules may be implemented in the form of hardware or in the form of software functional modules.

[0099] In another embodiment of the present invention, a computer device is provided, the computer device including a processor and a memory, the memory is used to store a computer program, the computer program includes program instructions, and the processor is used to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc., which are the computing core and control core of the terminal, which are suitable for implementing one or more instructions, and are specifically suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function; the processor described in the embodiment of the present invention can be used for the operation of the pre-stack inversion model fusion construction method based on synchronous squeezed wavelet transform.

[0100] In another embodiment of the present invention, the present invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device for storing programs and data. It is understandable that the computer-readable storage medium here can include both built-in storage media in the computer device and, of course, extended storage media supported by the computer device. The computer-readable storage medium provides a storage space, which stores the operating system of the terminal. In addition, one or more instructions suitable for being loaded and executed by the processor are also stored in the storage space, and these instructions can be one or more computer programs (including program codes). It should be noted that the computer-readable storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk memory. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the corresponding steps of the method for fusion construction of the prestack inversion model based on synchronous squeezed wavelet transform in the above embodiment.

[0101] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0102] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0103] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0104] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0105] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. Based on synchronous squeezing wavelet transform prestack inversion model fusion construction method, It is characterized in that include: Collect the pre-stack migration velocity field seismic data and the information of the top of the fracture-cavity body, and obtain the low-frequency velocity model within the time range of the seismic data and the high-frequency velocity model at and above the top of the fracture-cavity body; The synchronous squeezing wavelet transform is used to perform time-frequency analysis on the high-frequency model and the low-frequency model, and the synchronous squeezing wavelet transform coefficients are obtained; Based on the analysis of overlapping conditions of synchronously squeezed wavelet transform coefficients, the spectrum at a certain moment is obtained; Based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain a relatively wide-segment time-frequency spectrum of the data channel; The wide-band time-frequency spectra are acquired for all seismic traces to obtain relatively broadband model data.

2. The method for constructing a prestack inversion model based on synchronous squeezing wavelet transform according to claim 1, It is characterized in that Collect the pre-stack migration velocity field seismic data and the information of the top of the fracture-cavity body, and obtain the low-frequency velocity model within the time range of the seismic data and the high-frequency velocity model at and above the top of the fracture-cavity body: The prestack migration velocity field is constructed based on seismic data and is used to establish a low-frequency velocity model within the time range of seismic data. The information on the top of the fracture-cavity body is obtained based on the logging data, and the inverted high-frequency velocity model is constructed through the information on the top of the cave body.

3. The method for constructing a prestack inversion model based on synchronous squeezing wavelet transform according to claim 1, It is characterized in that The synchronous squeezing wavelet transform is used to perform time-frequency analysis on the high-frequency model and the low-frequency model, and the synchronous squeezing wavelet transform coefficients are obtained: Using Synchrosqueezing Wavelet Transform to Analyze High Frequency Model V high (x, t) and the low-frequency model V low (x, t) to perform time-frequency analysis. high (x i ,t) and V low (x i ,t) are obtained high (x i ,t,f) and F low (x i ,t,f); for data V high (x i , t), firstly, we use formula (1) to perform wavelet transform to obtain WF high (x i ,a,b): in, It is the complex conjugate of the mother wavelet. a and b are the scaling factor and translation factor, respectively. The scaling factor is related to the frequency, and the translation factor is related to the time. On this basis, the instantaneous frequency ω is calculated. s (x i ,a,b): After obtaining the instantaneous frequency, the time-frequency spectrum of the wavelet transform is squeezed in the frequency direction to obtain the synchronous squeezed wavelet transform coefficient F high (x i , t, f), that is: For data V low (x, t) can also be obtained according to formulas (1)-(3) low (x i , t, f).

4. The method for constructing a prestack inversion model based on synchronous squeezing wavelet transform according to claim 1, It is characterized in that Based on the analysis of the overlap of the synchronously squeezed wavelet transform coefficients, the spectrum at a certain moment is obtained: For the i-th data at a certain moment F high (x i , t i , f) and F low (x i , t i , f) perform time-frequency analysis; F high (x i , t i , f) and F low (x i , t i , f) There are two overlapping situations: Where two spectra intersect, an inverted Hamming window is added at the intersection to smooth the value; Where the two spectra do not intersect, the two spectra are directly added.

5. The method for constructing a prestack inversion model based on synchronous squeezing wavelet transform according to claim 4, It is characterized in that For the first case where the two spectra cross, an inverted Hamming window is added at the crossing point for smooth value extraction to obtain F all (x i , t i , f) is the spectrum of the i-th trace at time t i , that is:

6. The method for constructing a prestack inversion model based on synchronous squeezing wavelet transform according to claim 4, It is characterized in that Where the two spectra do not cross, the two spectra are directly added as F all (x i , t i , f), that is, F all (x i , t i , f)=F low (x i , t i , f)+F high (x i , t i , f).

7. The method for constructing a prestack inversion model based on synchronous squeezing wavelet transform according to claim 1, It is characterized in that Based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain a relatively wide-segment time-frequency spectrum of the data channel: Frequency fusion processing is performed on all moments of the i-th channel data, and the relatively wide-band time-frequency spectrum F of the i-th channel data can be obtained. all (x i , t, f), and then use the synchronous squeezing wavelet inverse transform to transform F all (x i , t, f) into V all (x i , t), that is: That is, relatively wide-band seismic trace data.

8. Based on synchronous squeezing wavelet transform prestack inversion model fusion construction system, It is characterized in that include: A data acquisition module is used to collect pre-stack migration velocity field seismic data and information on the top of the fracture-cavity body, and obtain a low-frequency velocity model within the time range of the seismic data and a high-frequency velocity model at and above the top of the fracture-cavity body; A synchronous squeezing wavelet transform coefficient acquisition module is used to perform time-frequency analysis on the high-frequency model and the low-frequency model using synchronous squeezing wavelet transform to obtain synchronous squeezing wavelet transform coefficients; The frequency fusion processing module is used to analyze the overlapping situation based on the synchronous squeezing wavelet transform coefficients to obtain the spectrum at a certain moment; based on the spectrum at a certain moment, frequency fusion processing is performed on all moments of a certain data channel to obtain the relatively wide-segment time-frequency spectrum of the data channel; wide-segment time-frequency spectrum is obtained for all seismic channels to obtain relatively wide-band model data.

9. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, It is characterized in that When the processor executes the computer program, the steps of the method for fusion construction of prestack inversion model based on synchronous squeezing wavelet transform are implemented as described in any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program. It is characterized in that When the computer program is executed by a processor, the steps of the method for constructing a prestack inversion model based on synchronous squeezing wavelet transform fusion are implemented as described in any one of claims 1 to 7.

Images