Time-lapse seismic matching inversion method and device based on temperature-sensitive rock physics experiment

Abstract

The present application relates to a kind of time-lapse seismic matching inversion method and device based on temperature-sensitive rock physics experiment, and the temperature-sensitive rock physics model is constructed by improving Gassmann equation, time-lapse well curve reconstruction, seismic identifiable scale 3D isostatic reservoir structure framework building and high-precision time-lapse seismic matching inversion and temperature field data conversion, obtain SAGD development different period oil sand reservoir temperature profile and temperature plane, fine delineation with temperature increase in the development process Reservoir internal space changes, improve oil and gas reservoir dynamic monitoring prediction accuracy, and combined with reservoir static attribute can be collectively for remaining oil distribution prediction direction.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas reservoir prediction, and in particular to a time-lapse seismic matching inversion method and device based on temperature-sensitive rock physics experiments. Background Art

[0002] With the development of geophysical technology, seismic inversion methods have been widely used in reservoir prediction, reserve calculation, and dynamic monitoring of oil reservoirs. They are one of the important methods for characterizing underground reservoir characteristics. Reservoir time-lapse seismic monitoring technology monitors the changes in seismic response characteristics generated during oil and gas reservoir development, and obtains changes in elastic parameters such as P-wave impedance and density. It further converts these changes into changes in reservoir parameters based on rock physics models to obtain changes in reservoir properties caused by oil and gas reservoir development, providing guidance for predicting remaining oil distribution and adjusting well deployment. In particular, when the study area has abundant seismic data but lacks well logging data and the sedimentation is complex, direct reservoir prediction using conventional seismic inversion methods will cause uncertainty in the inversion results, making it impossible to accurately characterize the internal differences in reservoir characteristics caused by development.

[0003] Oil sands reservoirs are often developed using steam-assisted gravity drainage (SAGD) using horizontal wells. Time-lapse seismic matching and inversion are effective methods for dynamic monitoring of oil sands SAGD development. However, the mechanisms for geophysical monitoring of reservoirs during SAGD development remain unclear. Furthermore, the lack of well log data makes secondary acquisition difficult and expensive, and the mismatch between the measured data and the seismic data acquisition years makes well log data inaccurately characterizing reservoir variations between years. Furthermore, the complex sedimentary conditions in the study area make conventional low-frequency model construction methods incapable of characterizing reservoir internal structure. Consequently, conventional inversion techniques cannot accurately characterize the differential characteristics of oil sands reservoirs caused by SAGD development. Summary of the Invention

[0004] The purpose of the present invention is to provide a time-lapse seismic matching inversion method and device based on temperature-sensitive rock physics experiments to solve the problem that conventional inversion technology in the prior art cannot accurately characterize the differential characteristics of oil sand reservoirs caused by SAGD development.

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

[0006] The present invention provides a time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments, comprising:

[0007] Based on temperature-sensitive rock physics experiments, a temperature-dependent rock physics model was established by improving the Gassmann equation, and velocity-temperature, density-temperature, and wave impedance-temperature versions of oil sands physical quantities were constructed.

[0008] Based on the constructed velocity-temperature and density-temperature plates, combined with the well temperature curves obtained through well logging, well logging curves that match the multi-period time-lapse seismic data are reconstructed;

[0009] Establish a 3D isochronal reservoir structure framework at a seismically recognizable scale and construct a low-frequency model that matches multi-period time-lapse seismic data;

[0010] Using the reconstructed logging curves and the constructed low-frequency model, high-precision seismic inversion results were obtained through the time-lapse seismic inversion method. Based on the constructed wave impedance-temperature plate, the high-precision seismic inversion results were converted into temperature field data to obtain the characteristic characterization of the temperature profile and temperature plane of the oil sand reservoir at different periods of SAGD development.

[0011] Furthermore, it also includes a method for determining the physical mechanism of oil sands: the determination method includes a laboratory ultrasonic testing system based on temperature-sensitive rock physics experiments, measuring the change pattern of elastic parameters of oil sand samples under different temperatures, pressures, and oil saturations, determining the rock physical mechanism of the change of elastic parameters of oil sand samples with temperature based on the test results, and improving the Gassmann equation after the physical mechanism of oil sands is determined.

[0012] Furthermore, the invention also includes an improved method of the Gassmann equation: the improved method includes introducing a rock skeleton weakening correction parameter and a pore volume change correction parameter to improve the traditional Gassmann equation;

[0013] The traditional rock physics modeling Gassmann equation is expressed as follows:

[0014]

[0015] μ sat =μ dry ;

[0016] Where K sat and K dry are the bulk moduli of rock under saturated and dry conditions, respectively; K g is the bulk modulus of mineral particles; K f is the pore fluid modulus; φ is the porosity; μ sat and μ dry are the shear moduli of rock under saturated and dry conditions, respectively;

[0017] The following parameters are introduced to correct the changes in rock skeleton caused by temperature changes, namely:

[0018] K dry =K dry (T);

[0019] μ dry =μ dry (T);

[0020] Where K dry is the rock skeleton weakening correction parameter, μ dry Correction parameter for pore volume changes.

[0021] Furthermore, the rock skeleton weakening correction parameter and pore volume change correction parameter are based on the heavy oil liquid point temperature (T lp ) is taken as the benchmark, and its parameter expression is further as follows:

[0022] K dry (T) = K dry (T lp )+ΔK×(T lp -T), T<T lp ;

[0023] K dry (T) = K dry (T lp )-ΔK×(T lp -T), T≥T lp ;

[0024] μ dry (T) = μ dry (T lp )+Δμ×(T lp -T), T<T lp ;

[0025] μ dry (T) = μ dry (T lp )-Δμ×(T lp -T), T≥T lp ;

[0026] Where ΔK and Δμ are the increments of bulk modulus and shear modulus with temperature, respectively.

[0027] Furthermore, the invention also includes a method for establishing a 3D isochronous reservoir structure framework at a seismically identifiable scale: the method includes establishing a 3D isochronous reservoir structure framework at a seismically identifiable scale by using a dip scanning calculation method based on geological body-guided coherence estimation and a dip-guided isochronous stratigraphic framework construction method;

[0028] The calculation method of the dip scanning for the geological body steering coherence estimation is as follows:

[0029]

[0030] Where p and q represent the x and y-direction components of the apparent tilt angle, τ represents time, and x j and y j is the local coordinate of the jth data based on the analysis point, Δxj and Δy j That is, the distance between the jth channel and the analysis point in the x and y directions, J is the total number of channels in the analysis window, K is the number of upper and lower sampling points centered on the analysis point in the time window, u and u H are seismic data and their Hilbert transform with respect to time;

[0031] The expression of the dip-oriented isochronous stratigraphic framework construction method is as follows:

[0032] Z i，j，k =(Z i-1，j，k +ΔX×Xdip i-1，j，k +Z i，j-1，k +ΔY×Ydip i，j-1，k ) / 2k=(0, 1, 2, ..., n);

[0033] Where Z i，j，k is the temporal depth of the kth layer at position (i, j), ΔX and ΔY are the sampling intervals in the X and Y directions, respectively, and Xdip and Ydip are the apparent dip angles in the X and Y directions, respectively, in ms / m.

[0034] Furthermore, it also includes a method for obtaining high-precision seismic inversion results, which includes obtaining deterministic inversion results including wave impedance and density parameters through a time-lapse seismic inversion method, using the wave impedance and density parameters as input parameters and participating in high-precision geostatistical inversion to obtain a high-precision wave impedance data body.

[0035] The present invention also provides an analysis device for a time-lapse seismic matching inversion method, comprising a processor configured to execute the above-mentioned time-lapse seismic matching inversion method.

[0036] The present invention adopts the above technical solution, which has the following beneficial effects:

[0037] 1. This invention determines the rock physical mechanism that is mainly affected by temperature during oil sand SAGD development based on temperature-sensitive rock physics experiments. By improving the traditional Gassmann equation to establish a rock physics model, the relationship between oil sand velocity, density, wave impedance and temperature is obtained.

[0038] 2. Based on the relationship between velocity, density and temperature, combined with the well temperature curve, reconstruct the time-lapse logging curve that matches the multi-period seismic data to resolve the uncertainty of the inversion results caused by the missing logging curve;

[0039] 3. Establish a seismically identifiable scale isochronous stratigraphic framework to address the challenge of characterizing the internal structure of reservoirs in complex sedimentary environments, improve the accuracy of identifying the internal structure of reservoirs, and overcome the multi-solution problem of inversion results caused by the low accuracy of low-frequency models in seismic inversion.

[0040] 4. By reconstructing logging curves and building a low-frequency model to match multi-period time-lapse seismic data, the time-lapse seismic inversion method is used to obtain high-precision seismic inversion results. Furthermore, the wave impedance-temperature relationship is used to convert the inversion results into temperature field data, and the temperature profile and temperature plane of the oil sand reservoir at different stages of SAGD development are obtained. The changes in the internal space of the reservoir as the temperature rises during the development process are precisely portrayed, and the accuracy of dynamic monitoring and prediction of oil and gas reservoirs is improved. Combined with the static properties of the reservoir, it can provide a direction for the prediction of remaining oil distribution. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Various other advantages and benefits will become apparent to those skilled in the art upon reading the detailed description of the preferred embodiment below. The accompanying drawings are for illustration purposes only and are not to be considered as limiting the present invention. Throughout the drawings, the same reference numerals are used to denote the same components. In the drawings:

[0042] Figure 1 This is a flow chart of a time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments provided by an embodiment of the present invention;

[0043] Figure 2 Figure (a) is a velocity-temperature relationship diagram of an oil sand sample obtained based on the improved Gassmann equation rock physics modeling of the present invention, Figure (b) is a density-temperature relationship diagram of an oil sand sample of the present invention, and Figure (c) is a wave impedance-temperature relationship diagram of an oil sand sample of the present invention;

[0044] Figure 3 This is a schematic diagram of time-lapse logging curves reconstructed based on temperature-sensitive rock physical quantities;

[0045] Figure 4 This is a schematic diagram of the technical process for constructing a 3D isochronous reservoir structure framework at a seismically recognizable scale;

[0046] Figure 5 It is a high-precision isochronous stratigraphic framework built in the research target area;

[0047] Figure 6 It is a high-precision wave impedance profile of different years obtained by time-lapse seismic matching inversion;

[0048] Figure 7 It is a temperature profile inversion process driven by rock physics;

[0049] Figure 8 Figure (a) is a schematic cross-sectional view of the temperature field in different years obtained by the method of the present invention, and Figure (b) is a schematic plan view of the temperature field in different years obtained by the method of the present invention. DETAILED DESCRIPTION

[0050] Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments described herein. Rather, these embodiments are provided to enable a more thorough understanding of the present invention and to fully convey the scope of the present invention to those skilled in the art.

[0051] Conventional inversion techniques cannot accurately depict the differential characteristics of oil sand reservoirs caused by SAGD development. This invention provides a time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments. By improving the Gassmann equation to construct a temperature-sensitive rock physics model, reconstructing time-lapse logging curves, building a 3D isochronous reservoir structure framework at a seismically recognizable scale, and performing high-precision time-lapse seismic matching inversion and temperature field data conversion, the method obtains temperature profiles and temperature planes of oil sand reservoirs at different stages of SAGD development. This method accurately depicts the spatial changes within the reservoir as temperature rises during development, improves the accuracy of dynamic monitoring and prediction of oil and gas reservoirs, and, combined with the static properties of the reservoir, provides guidance for the prediction of remaining oil distribution.

[0052] The scheme of the present invention is described in detail below through examples.

[0053] Example

[0054] like Figure 1 The present invention provides a time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments, comprising:

[0055] Based on temperature-sensitive rock physics experiments, a temperature-dependent rock physics model was established by improving the Gassmann equation, and velocity-temperature, density-temperature, and wave impedance-temperature versions of oil sands physical quantities were constructed.

[0056] Based on the constructed velocity-temperature and density-temperature plates, combined with the well temperature curves obtained through well logging, well logging curves that match the multi-period time-lapse seismic data are reconstructed;

[0057] Establish a 3D isochronal reservoir structure framework at a seismically recognizable scale and construct a low-frequency model that matches multi-period time-lapse seismic data;

[0058] Using the reconstructed logging curves and the constructed low-frequency model, high-precision seismic inversion results were obtained through the time-lapse seismic inversion method. Based on the constructed wave impedance-temperature plate, the high-precision seismic inversion results were converted into temperature field data to obtain the characteristic characterization of the temperature profile and temperature plane of the oil sand reservoir at different periods of SAGD development.

[0059] Furthermore, it also includes a method for determining the physical mechanism of oil sands. The determination method includes a laboratory ultrasonic testing system based on temperature-sensitive rock physics experiments, measuring the change pattern of elastic parameters of oil sand samples under different temperatures, pressures, and oil saturations, determining the rock physical mechanism of the change of elastic parameters of oil sand samples with temperature based on the test results, and improving the Gassmann equation after the physical mechanism of oil sands is determined.

[0060] Furthermore, the improved method of Gassmann equation includes introducing the correction parameter of rock skeleton weakening and the correction parameter of pore volume change to improve the traditional Gassmann equation;

[0061] The traditional rock physics modeling Gassmann equation is expressed as follows:

[0062]

[0063] μ sat =μ dry ;

[0064] Where K sat and K dry are the bulk moduli of rock under saturated and dry conditions, respectively; K g is the bulk modulus of mineral particles; K f is the pore fluid modulus; φ is the porosity; μ sat and μ dry are the shear moduli of rock under saturated and dry conditions, respectively;

[0065] The following parameters are introduced to correct the changes in rock skeleton caused by temperature changes, namely:

[0066] K dry =K dry (T);

[0067] μ dry =μ dry (T);

[0068] Where K dry is the rock skeleton weakening correction parameter, μ dry Correction parameter for pore volume changes.

[0069] Among them, if the rock skeleton weakening correction parameter and pore volume change correction parameter are based on the heavy oil liquid point temperature (T lp ) is taken as the benchmark, and its parameter expression is further as follows:

[0070] K dry (T) = K dry (T lp )+ΔK×(T lp -T), T<T lp ;

[0071] K dry (T) = K dry (T lp )-ΔK×(T lp -T), T≥T lp ;

[0072] μ dry (T) = μ dry (T lp )+Δμ×(T lp -T), T<T lp ;

[0073] μ dry (T) = μ dry (T lp )-Δμ×(T lp -T), T≥T lp ;

[0074] Where ΔK and Δμ are the increments of bulk modulus and shear modulus with temperature, respectively.

[0075] Figure 2 This is a quantitative version of the relationship between the elastic parameters of oil sand samples and temperature, obtained based on temperature-sensitive rock physics experiments and improved Gassmann equation rock physics modeling. The temperature-sensitive rock physics experiments show that the rock physical properties during oil sand SAGD development are mainly dominated by temperature, and as the temperature increases, the oil sands will undergo a phase change from solid to liquid. Because temperature increases can cause oil sand phase changes (49°C), correction parameters for rock skeleton weakening and pore volume change are introduced, and the traditional Gassmann equation is improved. A rock physics model that takes temperature changes into account is established to obtain a quantitative version of the relationship between oil sand sample parameters such as velocity and density and temperature.

[0076] Figure 3 This time-lapse well logging curve is reconstructed based on the relationship between elastic parameters and temperature in oil sand samples and temperature curves from temperature monitoring wells. Given the velocity and density curves for 2003 and temperature curves for different years, forward modeling of velocity and density curves for 2003, 2013, and 2015 is performed. The figure shows that the reconstructed 2003 velocity and density curves are consistent with the actual velocity and density curves, demonstrating the accuracy of the reconstructed well logging curves. The reconstructed 2013 and 2015 velocity and density curves show a decreasing trend with increasing temperature, consistent with experimental research. Reconstructing time-lapse well logging curves provides reliable basic data support for seismic inversion, overcoming the problem of missing well logging curves.

[0077] Furthermore, the invention also includes a method for establishing a 3D isochronous reservoir structure framework at a seismically recognizable scale, which includes integrating seismic, geological, and well logging information to perform fine well-seismic calibration to extract a comprehensive wavelet, and establishing a 3D isochronous reservoir structure framework at a seismically recognizable scale through a dip scanning calculation method based on geological body-guided coherence estimation and a dip-guided isochronous stratigraphic framework construction method.

[0078] The calculation method of the dip scanning for the geological body steering coherence estimation is as follows:

[0079]

[0080] Where p and q represent the x and y-direction components of the apparent tilt angle, τ represents time, and x j and y j is the local coordinate of the jth data based on the analysis point, Δx j and Δy j That is, the distance between the jth channel and the analysis point in the x and y directions, J is the total number of channels in the analysis window, K is the number of upper and lower sampling points centered on the analysis point in the time window, u and u H are seismic data and their Hilbert transform with respect to time;

[0081] The expression of the dip-oriented isochronous stratigraphic framework construction method is as follows:

[0082] Z i，j，k =(Z i-1，j，k +ΔX×Xdip i-1，j，k +Z i，j-1，k +ΔY×Ydip i，j-1，k ) / 2k=(0, 1, 2, ..., n);

[0083] Where Z i，j，k is the temporal depth of the kth layer at position (i, j), ΔX and ΔY are the sampling intervals in the X and Y directions, respectively, and Xdip and Ydip are the apparent dip angles in the X and Y directions, respectively, in ms / m.

[0084] Figure 4 This is a flow chart for constructing a high-precision 3D isochronous seismic stratigraphic framework in complex sedimentary areas. The dip angle calculated by the dip scanning method can accurately reflect the inclination and lateral fluctuation of the actual seismic phase axis, and is used to guide the tracking of small layers, thereby establishing an isochronous stratigraphic framework that is consistent with the seismic phase axis, solving the problem of characterizing the internal structure of the reservoir and improving the recognition accuracy of the internal spatial structure of the reservoir. The isochronous seismic stratigraphic framework obtained by this method is as follows: Figure 5 As shown, the accuracy of the low-frequency model is improved.

[0085] Figure 6This high-precision seismic inversion result is obtained through seismic inversion methods based on time-shift matching of well logs and time-shift matching of low-frequency models. As can be seen from the figure, while maintaining the original seismic structure, the inversion results improve the accuracy of reservoir spatial structure, lateral boundary identification, and vertical phase division, thereby enhancing the accuracy of differentiated characterization of target reservoir intervals.

[0086] Furthermore, the method for obtaining high-precision seismic inversion results includes obtaining deterministic inversion results including wave impedance and density parameters through time-lapse seismic inversion, using the wave impedance and density parameters as input parameters for high-precision geostatistical inversion to obtain a high-precision wave impedance data volume. High-precision geostatistical inversion is a conventional inversion method that includes inversion analysis based on lithologic spatial distribution trends, different lithologic ranges, and probability density functions of different lithologic types.

[0087] Figure 7 is based on Figure 2 The temperature field data conversion flow chart of the temperature-sensitive rock physical quantity version was established. The converted temperature field plane and profile results are as follows: Figure 8 shown.

[0088] The method of the present invention is applied to the actual work area to obtain the seismic high-resolution wave impedance inversion results (such as Figure 6 The rock physics version is further used to convert the inversion results into temperature field data (as shown in Figure 7 As shown in the figure), the temperature field monitoring planes and profiles of different years are obtained by conversion (as shown in the figure). Figure 8 As shown in Figure 2 ), the inversion results can accurately depict the dynamic characteristics of steam chamber development during SAGD development of oil sand reservoirs. With the continuous injection of high-temperature steam during reservoir development, the steam chamber exhibits dynamic vertical diffusion. This method provides a reliable basis for guiding the prediction of remaining oil distribution.

[0089] The present invention also provides an analysis device for a time-lapse seismic matching inversion method, comprising a processor configured to execute the above-mentioned time-lapse seismic matching inversion method.

[0090] The present invention is described in terms of flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to specific embodiments. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as a 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 generate 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 1A device that provides the functions specified in a block or multiple blocks.

[0091] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work 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 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

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

[0093] 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 aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments, characterized in that: The time-lapse seismic matching inversion method comprises: Based on temperature-sensitive rock physics experiments, a temperature-dependent rock physics model was established by improving the Gassmann equation, and velocity-temperature, density-temperature, and wave impedance-temperature versions of oil sands physical quantities were constructed. Based on the constructed velocity-temperature and density-temperature plates, combined with the well temperature curves obtained through well logging, well logging curves that match the multi-period time-lapse seismic data are reconstructed; Establish a 3D isochronal reservoir structure framework at a seismically recognizable scale and construct a low-frequency model that matches multi-period time-lapse seismic data; Using reconstructed well logging curves and a constructed low-frequency model, high-precision seismic inversion results are obtained through time-lapse seismic inversion. These high-precision seismic inversion results are then converted into temperature field data based on the constructed wave impedance-temperature plate to characterize the temperature profile and temperature plane of the oil sands reservoir during different periods of SAGD development. Also included is an improved method for the Gassmann equation: the improved method includes introducing a rock skeleton weakening correction parameter and a pore volume change correction parameter to improve the traditional Gassmann equation; The traditional rock physics modeling Gassmann equation is expressed as follows: ； ； Where, and are the bulk moduli of rock under saturated and dry conditions, respectively; is the bulk modulus of mineral particles; is the pore fluid modulus; is the porosity; and are the shear moduli of rock under saturated and dry conditions, respectively; The following parameters are introduced to correct the changes in rock skeleton caused by temperature changes, namely: ； ； Where, is the rock skeleton weakening correction parameter, Correction parameters for pore volume changes; The correction parameters for rock skeleton weakening and pore volume change are based on the heavy oil liquid point temperature. As a benchmark, its parameter expression is further: ； ； ； ； Where, and are the increments of bulk modulus and shear modulus with temperature, respectively.

2. The time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments according to claim 1 is characterized in that: It also includes a method for determining the physical mechanism of oil sands: the determination method includes a laboratory ultrasonic testing system based on temperature-sensitive rock physics experiments, measuring the change pattern of elastic parameters of oil sand samples under different temperatures, pressures, and oil saturations, determining the rock physical mechanism of the change of elastic parameters of oil sand samples with temperature based on the test results, and improving the Gassmann equation after the oil sand physical mechanism is determined.

3. The time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments according to claim 1 is characterized in that: The invention also includes a method for establishing a 3D isochronous reservoir structure framework at a seismically recognizable scale. The method includes integrating seismic, geological, and well logging information to perform fine well-seismic calibration to extract a comprehensive wavelet, and establishing a 3D isochronous reservoir structure framework at a seismically recognizable scale through a dip scanning calculation method using geological body-guided coherence estimation and a dip-guided isochronous stratigraphic framework construction method. The calculation method of the dip scanning for the geological body steering coherence estimation is as follows: Where, and Respectively and The apparent inclination component of the direction, Indicates time, and The first is based on the analysis point The local coordinates of the track data, and That is The path is relative to the analysis point and The spacing of the direction, is the total number of channels in the analysis window, is the number of up and down sampling points centered on the analysis point in the time window, and are seismic data and their Hilbert transform with respect to time; Indicates time intervals; The expression of the dip-oriented isochronous stratigraphic framework construction method is as follows: ； Where, yes Position The temporal depth of the layer, and They are and The sampling interval in the direction, and They are Direction and The apparent inclination angle of the direction, in units of .

4. The time-lapse seismic matching inversion method based on temperature-sensitive rock physics experiments according to claim 1 is characterized in that: It also includes a method for obtaining seismic inversion results, which includes obtaining deterministic inversion results including wave impedance and density parameters through a time-lapse seismic inversion method, using the wave impedance and density parameters as input parameters and participating in high-precision geostatistical inversion to obtain a high-precision wave impedance data body.

5. An analysis device for a time-lapse seismic matching inversion method, comprising a processor, wherein the processor is configured to execute the method according to any one of claims 1 to 4.

Images