CN112147689A

CN112147689A - Carbonate reservoir interpretation method and interpretation quantity version establishment method

Info

Publication number: CN112147689A
Application number: CN202010875557.5A
Authority: CN
Inventors: 李闯; 潘建国; 张虎权; 王宏斌; 赵建国; 李慧珍; 陈军; 丰超; 王振卿; 周俊峰
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2020-08-27
Filing date: 2020-08-27
Publication date: 2020-12-29
Anticipated expiration: 2040-08-27
Also published as: CN112147689B

Abstract

The invention provides a carbonate reservoir interpretation method and an interpretation quantity version establishment method. The method for establishing the interpretation quantity version comprises the following steps: establishing an interpretation version of a carbonate reservoir at low frequency and/or medium frequency and/or high frequency; acquiring elastic parameter test results of rock cores with different porosities of a target layer of a work area under low-frequency and/or medium-frequency and/or high-frequency conditions at different oil saturation degrees; the low frequency is consistent with the earthquake main frequency, the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; and (3) casting the elastic parameter test results of the cores with different porosities under the low-frequency and/or medium-frequency and/or high-frequency conditions under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the low-frequency and/or medium-frequency and/or high-frequency conditions, namely an interpretation quantity chart under the low-frequency and/or medium-frequency and/or high-frequency conditions of the carbonate reservoir.

Description

Carbonate reservoir interpretation method and interpretation quantity version establishment method

Technical Field

The invention relates to the geophysical technology in the petroleum industry, in particular to a carbonate reservoir interpretation method and an interpretation quantity edition building method.

Background

The quantitative explanation of earthquake is realized by means of rock physics research and establishing the accurate corresponding relation of reservoir, fluid and earthquake elastic parameters. Seismic petrophysical research involves two aspects: 1. carrying out rock physical theory modeling research; 2. research on rock physics laboratory testing technology.

The seismic petrophysical method which is widely applied in the petroleum industry is a rock physical theory modeling method, and works such as well logging curve prediction, reservoir fluid quantitative interpretation and the like are carried out by establishing a rock physical model of a research target. For the process of rock physical modeling, a quality control means is usually provided, and the process can be verified by using logging data and calibrated by referring to the result of a high-frequency test of a laboratory core.

This procedure has the problem of considerable poor accuracy, even if there is still a considerable deviation in interpretation by these calibrations. The effect is particularly bad when the established interpretation volume is used for interpreting seismic data, because the general rock physical modeling does not consider dispersion factors, the data for calibrating the interpretation volume are high frequency (ultrasonic measurement) and medium frequency (logging data), and the inversion result to be interpreted is low frequency (seismic frequency band).

Therefore, a method for establishing quantitative interpretation quantity versions of various frequency bands needs to be found again, quantitative interpretation can be carried out on data of different frequency dispersions, and interpretation results are more accurate and reasonable than those of the prior art in which a rock physical model is uniformly used for interpretation.

Disclosure of Invention

Aiming at the defects of the prior art, the invention aims to provide a method for establishing an interpretation quantity version of a carbonate reservoir. The quantitative version established by the method can meet the requirements of reservoir logging data quantitative interpretation, reservoir earthquake quantitative interpretation and reservoir high-frequency ultrasonic test data quantitative interpretation. Can be well used for carbonate reservoir interpretation.

The invention also aims to provide a carbonate reservoir interpretation method, which realizes quantitative interpretation of seismic inversion data, quantitative interpretation of well logging data and quantitative interpretation of high-frequency ultrasonic test data, predicts the porosity and saturation of underground carbonate and solves the problem of high-efficiency exploration of oil and gas.

In order to achieve the above object, the present invention provides a method for creating an interpretation quantity version of a carbonate reservoir, wherein the interpretation quantity version of the carbonate reservoir comprises interpretation quantity versions of the carbonate reservoir at low frequency and/or medium frequency and/or high frequency, the method comprising:

acquiring elastic parameter test results of rock cores with different porosities of a target layer of a work area under low-frequency and/or medium-frequency and/or high-frequency conditions at different oil saturation degrees; the elastic parameter test result is put into a coordinate system with coordinate axes of longitudinal and transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal and transverse wave parameter ratio and longitudinal wave impedance values of cores with different porosities and different oil saturation under low-frequency and/or medium-frequency and/or high-frequency conditions is established, namely a carbonate reservoir low-frequency and/or medium-frequency and/or high-frequency interpretation plate;

the low frequency used by the low-frequency conditional elastic parameter test result is consistent with the earthquake dominant frequency, the intermediate frequency used by the intermediate-frequency conditional elastic parameter test result is consistent with the logging frequency, and the high frequency used by the high-frequency conditional elastic parameter test result is consistent with the formation ultrasonic testing frequency;

the elastic parameters comprise a longitudinal wave parameter ratio and a longitudinal wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio.

In the method for establishing the carbonate reservoir interpretation version, the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees refer to the following steps: and (3) testing the elastic parameter of the rock cores with different porosities under the condition of a low-frequency seismic source under different oil saturation states. The medium-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees refer to that: and (3) testing the elastic parameter of the rock cores with different porosities under the condition of the medium-frequency seismic source under different oil saturation states. The high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees refer to that: and (3) testing the elastic parameter of the rock cores with different porosities under the condition of a high-frequency seismic source under different oil saturation states.

In the method for establishing the carbonate reservoir interpretation version, preferably, the elastic parameter test results of the cores with different porosities of the target zone in the low-frequency and/or medium-frequency and/or high-frequency conditions under different oil saturation degrees are obtained by the following steps:

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

and respectively carrying out elastic parameter tests on the cores with different porosities under different oil saturation states under the conditions of a low-frequency and/or medium-frequency and/or high-frequency seismic source to obtain the elastic parameter test results of the cores with different porosities under the conditions of low frequency and/or medium-frequency and/or high-frequency under different oil saturation.

In the method for establishing the carbonate reservoir interpretation version, preferably, the obtained work area target layer core is pretreated before the porosity is determined, and the pretreatment comprises sample preparation treatment and decontamination treatment; wherein, the decontamination treatment is preferably carried out by soaking by using an organic solution; the sample preparation treatment preferably comprises the steps of processing the core into a proper size and polishing the surface of the core.

In a specific embodiment, the different oil saturations include 100% oil saturation, 90% oil saturation, 80% oil saturation, 70% oil saturation, 60% oil saturation, 50% oil saturation, 40% oil saturation, 30% oil saturation, 20% oil saturation, 10% oil saturation.

In one embodiment, the low frequency is 20 Hz.

In one embodiment, the intermediate frequency is 10000 Hz.

In one embodiment, the high frequency is 100000 Hz.

In the above method for establishing the carbonate reservoir interpretation volume version, preferably, the establishing the carbonate reservoir interpretation volume version under the low frequency further includes: after the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal-transverse wave parameter ratio and longitudinal wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a first rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; the first rock physical theory model is calibrated by using the low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees.

In the above method for establishing the carbonate reservoir interpretation quantity version, preferably, establishing the frequency-down interpretation quantity version in the carbonate reservoir further includes: after the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a second rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; and the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

In the above method for establishing the carbonate reservoir interpretation volume version, preferably, the establishing the carbonate reservoir interpretation volume version at high frequency further includes: after the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal wave parameter ratio and a transverse wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a third rock physical theory model to determine the longitudinal wave parameter ratio and the transverse wave impedance value of the rock cores with different porosities under different oil saturation degrees; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

Due to the limited number of the cores and the randomness of the porosity, the actual porosity values are not uniformly distributed, and the mode of extrapolation and interpolation is adopted to supplement the throw values of more cores with different porosities under different oil saturation degrees, so that the established gauge plate has a wider interpretation range.

In the method for establishing the carbonate reservoir interpretation quantity version, preferably, the first petrophysical theoretical model is calibrated by using test data; the test data refers to the obtained low-frequency condition elastic parameter test result of the rock core with different porosities of the target layer under different oil saturation degrees. In one embodiment, the results of the calculations are calibrated with test data as the parameters are input to the theoretical model. The test data is used for calibrating the theoretical model, so that the progress of the theoretical model can be improved; the elastic parameters of the rock with different fluids can be calculated through a theoretical model, and therefore an explanation quantity version is established. Part of parameters in the theoretical model need to give initial values, in the prior art, the initial values are usually empirical values or rock elastic parameter values obtained by consulting a rock physics manual, but due to the complexity of rocks, the empirical values or the values recorded in the rock physics manual are not necessarily suitable for rocks in a research area, so that the parameter root can be attached to the actual rock condition by calibrating the theoretical model through test data, and the data obtained by model calculation is more accurate.

In the method for establishing the carbonate reservoir interpretation quantity version, preferably, the second rock physical theory model is calibrated by using test data; the test data refers to the obtained intermediate frequency condition elastic parameter test result of the rock core with different porosities of the target layer under different oil saturation degrees. In one embodiment, the results of the calculations are calibrated with test data as the parameters are input to the theoretical model. The test data is used for calibrating the theoretical model, so that the progress of the theoretical model can be improved; the elastic parameters of the rock with different fluids can be calculated through a theoretical model, and therefore an explanation quantity version is established. Part of parameters in the theoretical model need to give initial values, in the prior art, the initial values are usually empirical values or rock elastic parameter values obtained by consulting a rock physics manual, but due to the complexity of rocks, the empirical values or the values recorded in the rock physics manual are not necessarily suitable for rocks in a research area, so that the parameter root can be attached to the actual rock condition by calibrating the theoretical model through test data, and the data obtained by model calculation is more accurate.

In the method for establishing the carbonate reservoir interpretation quantity version, preferably, the third petrophysical theoretical model is calibrated by using test data; the test data refers to the obtained high-frequency condition elastic parameter test result of the rock core with different porosities of the target layer under different oil saturation degrees. In one embodiment, the results of the calculations are calibrated with test data as the parameters are input to the theoretical model. The test data is used for calibrating the theoretical model, so that the progress of the theoretical model can be improved; the elastic parameters of the rock with different fluids can be calculated through a theoretical model, and therefore an explanation quantity version is established. Part of parameters in the theoretical model need to give initial values, in the prior art, the initial values are usually empirical values or rock elastic parameter values obtained by consulting a rock physics manual, but due to the complexity of rocks, the empirical values or the values recorded in the rock physics manual are not necessarily suitable for rocks in a research area, so that the parameter root can be attached to the actual rock condition by calibrating the theoretical model through test data, and the data obtained by model calculation is more accurate.

In the above method for establishing the carbonate reservoir interpretation volume version, preferably, the establishing the carbonate reservoir interpretation volume version under the low frequency further includes: after the elastic parameter test results of the cores with different porosities under the low-frequency conditions with different oil saturation degrees are cast into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different oil saturation degrees under the low-frequency condition and/or fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the low-frequency condition.

In the above method for establishing the carbonate reservoir interpretation quantity version, preferably, establishing the frequency-down interpretation quantity version in the carbonate reservoir further includes: the method comprises the steps of putting test results of elastic parameters of cores with different porosities under the medium-frequency condition under different oil saturation into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and fitting a relation trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities under the medium-frequency condition and/or fitting a relation trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities under the medium-frequency condition.

In the above method for establishing the carbonate reservoir interpretation volume version, preferably, the establishing the carbonate reservoir interpretation volume version at high frequency further includes: the method comprises the steps of putting test results of elastic parameters of cores with different porosities under high-frequency conditions under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different oil saturation degrees under the high-frequency condition and/or fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the high-frequency condition.

The invention also provides a carbonate reservoir interpretation method, wherein the method comprises the following steps:

performing carbonate reservoir interpretation based on pre-stack seismic inversion results and/or performing carbonate reservoir interpretation based on logging data and/or performing carbonate reservoir interpretation based on formation ultrasonic test data;

the carbonate reservoir interpretation based on the pre-stack seismic inversion result comprises the following steps: establishing a low-frequency interpretation quantity version of the carbonate reservoir of the target layer in the work area by adopting the establishment method of the carbonate reservoir interpretation quantity version;

acquiring a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data volume and a longitudinal wave impedance data volume according to the pre-stack seismic inversion result;

intersecting longitudinal and transverse wave parameters with a data body and a longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation amount version under low frequency;

wherein the carbonate reservoir interpretation based on the well log data comprises:

establishing a carbonate reservoir under-frequency interpretation quantity version of a target layer of a work area by adopting the establishment method of the carbonate reservoir interpretation quantity version;

acquiring target logging data of a work area, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the logging data;

intersecting longitudinal and transverse wave parameters with a data body and a longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation quantity version under intermediate frequency;

wherein, the carbonate reservoir interpretation based on the formation ultrasonic test data comprises:

establishing a high-frequency interpretation version of the carbonate reservoir in the target layer of the work area by adopting the establishment method of the carbonate reservoir interpretation version;

acquiring stratum ultrasonic test data of a work area, and determining a longitudinal and transverse wave parameter ratio data body and a longitudinal wave impedance data body according to the stratum ultrasonic test data;

and intersecting the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation quantity version under high frequency.

The invention also provides a system for establishing the carbonate reservoir interpretation quantity version, wherein the system comprises:

an elastic parameter acquisition module: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees, and/or obtaining medium-frequency condition elastic parameter test results of the rock cores with different porosities of the target layer of the work area under different oil saturation degrees, and/or obtaining high-frequency condition elastic parameter test results of the rock cores with different porosities of the target layer of the work area under different oil saturation degrees; the low frequency is consistent with the earthquake main frequency, the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

an interpretation quantity version establishing module: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; and/or, the method is used for casting the test results of the elastic parameters of the rock cores with different porosities under the medium-frequency condition under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing query charts of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition, namely the interpretation quantity chart of the rock carbonate reservoir under the medium frequency; and/or, the method is used for casting the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis as the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition, namely the carbonate reservoir high-frequency interpretation quantity chart.

In the above system for establishing an interpretation version of a carbonate reservoir, preferably, the elastic parameter obtaining module includes:

a core obtaining submodule: the method comprises the steps of obtaining a work area target layer rock core;

a porosity determination submodule: the porosity of a target layer core of the work area is determined;

an elastic parameter acquisition submodule: the method is used for respectively testing the elastic parameters of the rock cores with different porosities under different oil saturation states under the condition of a low-frequency seismic source to obtain the low-frequency elastic parameter test results of the rock cores with different porosities under different oil saturations; and/or, the elastic parameter testing device is used for respectively testing the elastic parameters of the rock cores with different porosities under different oil saturation states under the condition of the medium-frequency seismic source to obtain the elastic parameter testing results of the rock cores with different porosities under the medium-frequency conditions under different oil saturations; and/or, the elastic parameter testing device is used for respectively testing the elastic parameters of the rock cores with different porosities under different oil saturation states under the condition of a high-frequency seismic source to obtain the high-frequency elastic parameter testing results of the rock cores with different porosities under different oil saturations;

more preferably, the elasticity parameter acquiring module further comprises:

a core processing submodule: and the pretreatment device is used for pretreating the obtained work area target layer rock core, and the pretreatment comprises sample preparation treatment and decontamination treatment.

In the above system for establishing an interpretation quantity version of a carbonate reservoir, preferably, the interpretation quantity version establishing module is configured to cast the low-frequency condition elastic parameter test results of cores with different porosities at different oil saturation degrees into a coordinate system whose coordinate axes are a longitudinal-transverse wave parameter ratio and a longitudinal wave impedance, and when an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the low-frequency condition is the low-frequency interpretation quantity version of the carbonate reservoir, the interpretation quantity version establishing module is further configured to cast the low-frequency condition elastic parameter test results of the cores with different porosities at different oil saturation degrees into the coordinate system whose coordinate axes are the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance, and then further supplement the cast point values of the cores with different porosities at different oil saturation degrees in the coordinate system by using an extrapolation and interpolation method, determining longitudinal and transverse wave parameter ratios and longitudinal wave impedance values of rock cores with different porosities under different oil saturation degrees by adopting a first rock physical theory model through an extrapolation and interpolation method; the first rock physical theory model is calibrated by using low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

the interpretation quantity version establishing module is used for casting the intermediate frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and when establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the intermediate frequency condition, namely the interpretation quantity version under the intermediate frequency of the carbonate reservoir, the interpretation quantity version establishing module is further used for casting the intermediate frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into the coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and then further supplementing the casting values of the rock cores with different porosities under different oil saturation degrees in the coordinate system by adopting an extrapolation and interpolation method, wherein the extrapolation and interpolation method adopts a second rock physical theory model to perform longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees Determining; the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

the interpretation quantity version establishing module is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of the longitudinal wave parameter ratio and the longitudinal wave impedance, establishing an inquiry chart of the longitudinal wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely the interpretation quantity version of the carbonate reservoir under the high frequency, the interpretation quantity version establishing module is further used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into the coordinate system with the coordinate axis of the longitudinal wave parameter ratio and the longitudinal wave impedance value, further supplementing the casting values of the rock cores with different porosities under different oil saturation degrees in the coordinate system by adopting an extrapolation and interpolation method, wherein the extrapolation and interpolation method adopts a third rock physical theory model to perform longitudinal wave parameter ratio and longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees Determining; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

In the above system for establishing an interpretation quantity version of a carbonate reservoir, preferably, the interpretation quantity version establishing module is configured to cast the low-frequency condition elastic parameter test results of cores with different porosities under different oil saturation degrees into a coordinate system whose coordinate axes are a longitudinal-transverse wave parameter ratio and a longitudinal wave impedance, and when an inquiry chart of longitudinal-transverse wave parameter ratio and longitudinal wave impedance values of cores with different porosities and different oil saturation degrees under a low-frequency condition is the low-frequency interpretation quantity version of the carbonate reservoir, the interpretation quantity version establishing module is further configured to cast the low-frequency condition elastic parameter test results of cores with different porosities under different oil saturation degrees into a coordinate system whose coordinate axes are a longitudinal-transverse wave parameter ratio and a longitudinal wave impedance, and then fit a relationship trend line between the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance of the cores with different oil saturation degrees under the low-frequency condition and/or fit a relationship trend line between the longitudinal-transverse wave parameters and longitudinal wave impedance of the cores with different porosity degrees under the low-frequency condition A trend line of the relationship between the parameter ratio and the longitudinal wave impedance value;

the interpretation quantity plate establishing module is used for casting the test results of the medium-frequency condition elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, when an inquiry chart of longitudinal and transverse wave parameter ratios and longitudinal wave impedance values of rock cores with different porosities and oil saturation degrees under the medium-frequency condition is established, namely an interpretation chart under the medium-frequency condition of a carbonate reservoir, the interpretation quantity plate establishing module is further used for casting the test results of the medium-frequency condition elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, fitting a relation trend line of longitudinal and transverse wave parameter ratios of different oil saturation degrees of each porosity rock core and longitudinal wave impedance values under the medium-frequency condition and/or fitting a relation trend line of longitudinal and transverse wave parameter ratios of different porosity rock cores and longitudinal wave impedance values under the medium-frequency condition;

the interpretation quantity plate establishing module is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal and transverse wave parameter ratio and longitudinal wave impedance, when an inquiry chart of longitudinal and transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation under a high-frequency condition is established, namely an interpretation chart of carbonate reservoir under high frequency, the interpretation quantity plate establishing module is further used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal and transverse wave parameter ratio and longitudinal wave impedance, and fitting a relation trend line of the longitudinal wave impedance value and the longitudinal wave parameter ratio of the core with different oil saturation degrees under the high-frequency condition and/or fitting a relation trend line of the longitudinal wave impedance value and the longitudinal wave parameter ratio of the core with different oil saturation degrees under the high-frequency condition.

The invention also provides a carbonate reservoir interpretation system, wherein the system comprises:

the system for establishing the carbonate reservoir interpretation quantity version is used for establishing a low-frequency interpretation quantity version of a target layer carbonate reservoir in a work area and/or a medium-frequency interpretation quantity version of the target layer carbonate reservoir and/or a high-frequency interpretation quantity version of the target layer carbonate reservoir;

work area target layer data acquisition module: the method comprises the steps of obtaining a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the pre-stack seismic inversion result; and/or acquiring target logging data of a work area, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the logging data; and/or acquiring stratum ultrasonic test data of a work area, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the stratum ultrasonic test data;

reservoir quantitative interpretation module: the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body determined according to a pre-stack seismic inversion result, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation version under low frequency; and/or, the device is used for utilizing the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body to be intersected, and carrying out porosity and/or saturation quantitative interpretation by referring to an interpretation quantity version under the intermediate frequency; and/or, the device is used for carrying out intersection by utilizing the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body, and carrying out porosity and/or saturation quantitative interpretation by referring to the interpretation amount version under high frequency.

In the carbonate reservoir interpretation system described above, preferably, the pre-stack seismic inversion results include a compressional impedance data volume and a shear impedance data volume.

In the carbonate reservoir interpretation system described above, the well log data preferably comprises shear wave well logs, compressional wave well logs and density well logs.

In the carbonate reservoir interpretation system described above, the formation ultrasonic test data preferably includes compressional wave data, shear wave data, and density data.

The invention also provides a device for establishing the carbonate reservoir interpretation volume version, which comprises a processor and a memory; wherein the content of the first and second substances,

a memory for storing a computer program;

and the processor is used for realizing the steps of the method for establishing the carbonate reservoir interpretation volume edition when executing the program stored in the memory.

The invention also provides a carbonate reservoir interpretation method and a device, which comprise a processor and a memory; wherein the content of the first and second substances,

a memory for storing a computer program;

and the processor is used for realizing the steps of the carbonate reservoir interpretation method when executing the program stored in the memory.

The present invention also provides a computer readable storage medium storing one or more programs, the one or more programs being executable by one or more processors to perform the steps of the method for establishing an interpretation volume version of a carbonate reservoir as described above.

The present invention also provides a computer readable storage medium having one or more programs stored thereon that are executable by one or more processors to perform the steps of the carbonate reservoir interpretation method described above.

In the last decade, the technology for testing core samples in a laboratory is greatly developed, particularly the technology for testing full frequency bands is provided, the technical scheme provided by the invention establishes quantitative interpretation quantity versions of low, medium and high frequency bands by testing the variable fluid saturation of the low, medium and high frequency bands of the reservoir core, and meets the multiple purposes of logging evaluation, reservoir earthquake quantitative interpretation and high-frequency ultrasonic testing analysis. Meanwhile, full-band calibration data are provided for researching the rock physical theoretical model with the frequency dispersion characteristic, and the development and the perfection of the frequency dispersion theoretical model are promoted.

The establishment of quantitative explanation quantity versions in the prior art is obtained based on a rock physical theory model. However, the rock physics theory models are various, the applicable conditions of each model are different, and the input control parameters are complex. According to the technical scheme provided by the invention, based on the elastic parameter test results of the target layer core in the research area under the conditions of low frequency, medium frequency and high frequency, a low frequency quantitative interpretation quantity version is established for interpreting seismic data, a medium frequency quantitative interpretation quantity version is established for interpreting logging data, and a high frequency quantitative interpretation quantity version is established for interpreting high frequency ultrasonic test data. The measurement version established by the method is obtained by testing the same high-frequency according to the actual core of the research area, the same low-frequency of earthquake, the same medium-frequency of logging and the same high-frequency of high-frequency ultrasonic testing, and the porosity, the saturation and the fluid-containing property of the known core, so that the method has higher reliability. The technical scheme provided by the invention solves the problem of the applicable condition of the model in the process of establishing the quantitative version of the theoretical model; compared with high-frequency rock physical experiment tests, the method overcomes the problem of frequency dispersion in the tests.

The carbonate reservoir interpretation method provided by the invention is characterized in that seismic inversion data are quantitatively interpreted by means of a low-frequency interpretation quantity version established based on a fluid-containing core elastic parameter test result under a low-frequency condition, logging data are quantitatively interpreted by means of a medium-frequency interpretation quantity version established based on a fluid-containing core elastic parameter test result under a medium-frequency condition, high-frequency ultrasonic test data are quantitatively interpreted by means of a high-frequency interpretation quantity version established based on a fluid-containing core elastic parameter test result under a high-frequency condition, and quantitative interpretation of the porosity and/or the saturation of the underground carbonate reservoir is realized. Compared with the prior art, the carbonate reservoir interpretation method provided by the invention obviously improves the prediction precision of the reservoir porosity and solves the problem of high-efficiency exploration of oil and gas; specifically, the improvement of the porosity prediction precision is beneficial to accurately estimating the oil and gas resource quantity, and the effective reservoir drilling rate, the drilling success rate and the recovery benefit are improved.

Drawings

Fig. 1 is a schematic flow chart of the low-frequency interpretation volume version in the method for establishing the carbonate reservoir interpretation volume version according to the embodiment of the present invention.

Fig. 2 is a schematic diagram of a process for establishing a medium-frequency interpretation volume version in the method for establishing a carbonate reservoir interpretation volume version according to an embodiment of the present invention.

Fig. 3 is a schematic diagram of a process for establishing a high-frequency interpretation volume edition in the method for establishing a carbonate reservoir interpretation volume edition according to an embodiment of the present invention.

Fig. 4 is a schematic flow chart of carbonate reservoir interpretation based on pre-stack seismic inversion results in the carbonate reservoir interpretation method according to an embodiment of the present invention.

Fig. 5 is a schematic flow chart of carbonate reservoir interpretation based on well log data in the carbonate reservoir interpretation method according to an embodiment of the present invention.

Fig. 6 is a schematic flow chart of carbonate reservoir interpretation based on formation ultrasonic test data in the carbonate reservoir interpretation method according to an embodiment of the present invention.

Fig. 7 is a low frequency interpretation version of a carbonate reservoir established in an embodiment of the invention.

Fig. 8 is a frequency interpreted version of a carbonate reservoir established in an embodiment of the present invention.

Fig. 9 is a high frequency interpretation of a carbonate reservoir established in an embodiment of the invention.

Fig. 10 is a low frequency, medium frequency, and high frequency interpretation magnitude comparison chart of a carbonate reservoir established in an embodiment of the present invention.

Fig. 11 is a schematic structural diagram of a system for establishing an explanatory quantity version of a carbonate reservoir according to an embodiment of the present invention.

Fig. 12 is a schematic structural diagram of an explanation system for a carbonate reservoir according to an embodiment of the present invention.

Fig. 13 is a schematic structural diagram of an apparatus for establishing an explanatory quantity version of a carbonate reservoir according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in detail and completely with reference to the drawings in the embodiments of the present invention. It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The principles and spirit of the present invention are explained in detail below with reference to several representative embodiments of the invention.

In order to achieve the above object, an embodiment of the present invention provides a method for creating an interpretation volume version of a carbonate reservoir, wherein the method includes creating a low-frequency interpretation volume version of the carbonate reservoir.

Another embodiment of the invention provides a method for establishing an interpretation quantity version of a carbonate reservoir, wherein the method comprises establishing a frequency interpretation quantity version of the carbonate reservoir.

Another embodiment of the invention provides a method for establishing an interpretation quantity version of a carbonate reservoir, wherein the method comprises establishing a high-frequency interpretation quantity version of the carbonate reservoir.

Another embodiment of the invention provides a method for establishing an interpretation quantity version of a carbonate reservoir, wherein the method comprises establishing a low-frequency interpretation quantity version of the carbonate reservoir and establishing a medium-frequency interpretation quantity version of the carbonate reservoir.

Another embodiment of the invention provides a method for establishing an interpretation quantity version of a carbonate reservoir, wherein the method comprises establishing a low-frequency interpretation quantity version of the carbonate reservoir and establishing a high-frequency interpretation quantity version of the carbonate reservoir.

Another embodiment of the invention provides a method for establishing an interpretation quantity version of a carbonate reservoir, wherein the method comprises establishing a medium frequency interpretation quantity version of the carbonate reservoir and establishing a high frequency interpretation quantity version of the carbonate reservoir.

Another embodiment of the invention provides a method for establishing an interpretation quantity version of a carbonate reservoir, wherein the method comprises establishing a low-frequency interpretation quantity version of the carbonate reservoir, establishing a medium-frequency interpretation quantity version of the carbonate reservoir and establishing a high-frequency interpretation quantity version of the carbonate reservoir.

Referring to fig. 1, establishing a low-frequency interpretation version of a carbonate reservoir includes:

step S11: obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the low frequency is consistent with the seismic dominant frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

step S12: and (3) casting the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the low-frequency condition, namely a carbonate reservoir low-frequency interpretation quantity chart.

With reference to fig. 2, establishing a lower frequency interpretation version of a carbonate reservoir comprises:

step S21: acquiring intermediate frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the intermediate frequency is consistent with the logging frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

step S22: and (3) casting the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal-transverse wave parameter ratios and longitudinal wave impedance, and establishing a query plate of the longitudinal-transverse wave parameter ratios and the longitudinal wave impedance values of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition, namely a carbonate reservoir medium-frequency interpretation plate.

With reference to fig. 3, establishing a high-frequency interpretation version of a carbonate reservoir includes:

step S31: obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

step S32: and (3) casting the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition, namely a carbonate reservoir high-frequency interpretation chart.

Wherein, the low-frequency condition elasticity parameter test results of the rock cores with different porosities under different oil saturation degrees refer to that: and (3) testing the elastic parameter of the rock cores with different porosities under the condition of a low-frequency seismic source under different oil saturation states.

Wherein, the medium frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees refer to: and (3) testing the elastic parameter of the rock cores with different porosities under the condition of the medium-frequency seismic source under different oil saturation states.

Wherein the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees refer to that: and (3) testing the elastic parameter of the rock cores with different porosities under the condition of a high-frequency seismic source under different oil saturation states.

In one embodiment, the low-frequency condition elasticity parameter test results of the cores with different porosities of the target zone of the work area under different oil saturation degrees are obtained by the following steps:

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

and respectively carrying out elastic parameter tests on the cores with different porosities under different oil saturation states under the condition of a low-frequency seismic source to obtain the elastic parameter test results of the cores with different porosities under the low-frequency conditions under different oil saturations.

In one embodiment, the medium-frequency condition elasticity parameter test results of the cores with different porosities of the target zone of the work area under different oil saturation degrees are obtained by the following steps:

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

and respectively carrying out elastic parameter tests on the rock cores with different porosities under different oil saturation states under the condition of a medium-frequency seismic source to obtain medium-frequency elastic parameter test results of the rock cores with different porosities under different oil saturations.

In one embodiment, the high-frequency condition elasticity parameter test results of the cores with different porosities of the target zone of the work area under different oil saturation degrees are obtained by the following steps:

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

and respectively carrying out elastic parameter tests on the rock cores with different porosities under different oil saturation states under the condition of a high-frequency seismic source to obtain high-frequency elastic parameter test results of the rock cores with different porosities under different oil saturations.

In one embodiment, the obtained work area target layer core is pretreated before the porosity is determined, wherein the pretreatment comprises sample preparation treatment and decontamination treatment; wherein, the decontamination treatment is preferably carried out by soaking by using an organic solution; the sample preparation treatment preferably comprises the steps of processing the core into a proper size and polishing the surface of the core. The obtained core usually contains a certain amount of residual crude oil and asphalt, and the core is soaked in an organic solution for decontamination treatment, so that the residual crude oil and asphalt in the core can be effectively removed.

In one embodiment, the core of the target zone is obtained by coring a borehole.

In one embodiment, the core was processed into a sample having a diameter of 3.8 centimeters and a length of greater than 5 centimeters.

In one embodiment, the different oil saturations include 100% oil saturation, 90% oil saturation, 80% oil saturation, 70% oil saturation, 60% oil saturation, 50% oil saturation, 40% oil saturation, 30% oil saturation, 20% oil saturation, 10% oil saturation.

In one embodiment, the low frequency is 20 Hz.

In one embodiment, the intermediate frequency is 10000 Hz.

In one embodiment, the high frequency is 100000 Hz.

In one embodiment, creating the carbonate reservoir low frequency interpretation volume further comprises: after the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal-transverse wave parameter ratio and longitudinal wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a first rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; wherein a first petrophysical theoretical model (e.g., a jet model of a dual-hole structure hypothesis derived by Gurevich et al (2010)) is calibrated using low-frequency conditional elastic parameter test results of cores of different porosities at different oil saturations.

In one embodiment, establishing the lower interpretation volume in the carbonate reservoir further comprises: after the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a second rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; wherein a second petrophysical theoretical model (e.g., a jet model of a dual-pore structure hypothesis derived by Gurevich et al (2010)) is calibrated using medium frequency conditional elastic parameter test results of cores of different porosities at different oil saturations.

In one embodiment, establishing the interpretation at high frequency version of the carbonate reservoir further comprises: after the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal wave parameter ratio and a transverse wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a third rock physical theory model to determine the longitudinal wave parameter ratio and the transverse wave impedance value of the rock cores with different porosities under different oil saturation degrees; wherein a third petrophysical theoretical model (e.g., a jet model of a dual-hole structure hypothesis derived by Gurevich et al (2010)) is calibrated using high-frequency conditional elastic parameter test results of cores of different porosities at different oil saturations.

Due to the limited number of the cores and the randomness of the porosity, the actual porosity values are not uniformly distributed, and the interpolation method supplements the throw-in values of more cores with different porosities under different oil saturation degrees, so that the established gauge plate has a wider interpretation range. And when the theoretical model input parameters are selected, the calculated result is calibrated by using the test data.

In one embodiment, creating the carbonate reservoir low frequency interpretation volume further comprises: after the elastic parameter test results of the cores with different porosities under the low-frequency conditions with different oil saturation degrees are cast into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different oil saturation degrees under the low-frequency condition and/or fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the low-frequency condition.

In one embodiment, establishing the lower interpretation volume in the carbonate reservoir further comprises: the method comprises the steps of putting test results of elastic parameters of cores with different porosities under the medium-frequency condition under different oil saturation into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and fitting a relation trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities under the medium-frequency condition and/or fitting a relation trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities under the medium-frequency condition.

In one embodiment, establishing the interpretation at high frequency version of the carbonate reservoir further comprises: the method comprises the steps of putting test results of elastic parameters of cores with different porosities under high-frequency conditions under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different oil saturation degrees under the high-frequency condition and/or fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the high-frequency condition.

An embodiment of the invention provides a carbonate reservoir interpretation method, wherein the method comprises the following steps: and (4) carrying out carbonate reservoir interpretation based on the pre-stack seismic inversion result.

Another embodiment of the present invention provides a carbonate reservoir interpretation method, wherein the method comprises: carbonate reservoir interpretation is performed based on the well log data.

Another embodiment of the present invention provides a carbonate reservoir interpretation method, wherein the method comprises: and carrying out carbonate reservoir interpretation based on the formation ultrasonic testing data.

Another embodiment of the present invention provides a carbonate reservoir interpretation method, wherein the method comprises: and performing carbonate reservoir interpretation based on the pre-stack seismic inversion result and performing carbonate reservoir interpretation based on the logging data.

Another embodiment of the present invention provides a carbonate reservoir interpretation method, wherein the method comprises: and performing carbonate reservoir interpretation based on the pre-stack seismic inversion result and performing carbonate reservoir interpretation based on the formation ultrasonic test data.

Another embodiment of the present invention provides a carbonate reservoir interpretation method, wherein the method comprises: and carrying out carbonate reservoir interpretation based on the logging data and carrying out carbonate reservoir interpretation based on the formation ultrasonic testing data.

Another embodiment of the present invention provides a carbonate reservoir interpretation method, wherein the method comprises: carbonate reservoir interpretation is performed based on pre-stack seismic inversion results, carbonate reservoir interpretation is performed based on logging data, and carbonate reservoir interpretation is performed based on formation ultrasonic testing data.

Referring to fig. 4, the carbonate reservoir interpretation based on the prestack seismic inversion results includes:

step S41: obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the low frequency is consistent with the seismic dominant frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

step S42: the method comprises the steps of putting test results of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation into a coordinate system with the coordinate axis being longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing query charts of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturations under the low-frequency conditions, namely a carbonate reservoir low-frequency interpretation chart;

step S43: acquiring a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data volume and a longitudinal wave impedance data volume according to the pre-stack seismic inversion result;

step S44: intersecting longitudinal and transverse wave parameters with a data body and a longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation amount version under low frequency;

therein, referring to fig. 5, carbonate reservoir interpretation based on well log data includes:

step S51: acquiring intermediate frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the intermediate frequency is consistent with the logging frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

step S52: the method comprises the steps that (1) a test result of medium-frequency condition elastic parameters of rock cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate;

step S53: acquiring target logging data of a work area, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the logging data;

step S54: intersecting longitudinal and transverse wave parameters with a data body and a longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation quantity version under intermediate frequency;

referring to fig. 6, the carbonate reservoir interpretation based on formation ultrasonic test data includes:

step S61: obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

step S62: the method comprises the steps of putting high-frequency condition elastic parameter test results of cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis being longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing query charts of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition, namely carbonate reservoir high-frequency interpretation charts;

step S63: acquiring stratum ultrasonic test data of a work area, and determining a longitudinal and transverse wave parameter ratio data body and a longitudinal wave impedance data body according to the stratum ultrasonic test data;

step S64: and intersecting the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation quantity version under high frequency.

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

In one embodiment, the low frequency is 20 Hz.

In one embodiment, the intermediate frequency is 10000 Hz.

In one embodiment, the high frequency is 100000 Hz.

In one embodiment, creating the carbonate reservoir low frequency interpretation volume further comprises: after the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal-transverse wave parameter ratio and longitudinal wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a first rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; the first rock physical theory model is calibrated by using the low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees.

In one embodiment, establishing the lower interpretation volume in the carbonate reservoir further comprises: after the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a second rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; and the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

In one embodiment, establishing the interpretation at high frequency version of the carbonate reservoir further comprises: after the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal wave parameter ratio and a transverse wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a third rock physical theory model to determine the longitudinal wave parameter ratio and the transverse wave impedance value of the rock cores with different porosities under different oil saturation degrees; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

In one embodiment, the pre-stack seismic inversion results include a compressional impedance data volume and a shear impedance data volume. The conventional pre-stack seismic inversion result generally comprises a longitudinal wave impedance data volume, a transverse wave impedance data volume and a density data volume, and the wave impedance data volume and the transverse wave impedance data volume are directly selected.

In one embodiment, the well log data includes shear wave well logs, compressional wave well logs, and density well logs. The transverse wave logging curve and the longitudinal wave logging curve may be impedance logging curves or velocity logging curves, and the mutual conversion between the velocity and the impedance can be completed according to an impedance and velocity conversion formula, specifically, the impedance is velocity × density.

In one embodiment, the formation ultrasonic test data includes compressional wave data, shear wave data, and density data.

In another embodiment of the present invention, a carbonate reservoir interpretation method is provided for interpreting the contents of the carbonate reservoir of the aoto family in a basin of western China.

The effective reservoir stratum of the carbonate rock in the area is mainly distributed in a high-energy phase zone area, and primary pores are well preserved and are only filled in a small amount. Lithology is mainly dolomite and grey matter dolomite, and lithology transverse distribution is relatively stable. The reservoir development area has reflection characteristics on a seismic section, but the waveform changes are complex. By utilizing an inversion method such as poststack inversion, the development and distribution area of the favorable reservoir can be qualitatively predicted. In order to further develop efficiently, a good-quality dessert is selected, and pre-stack inversion is further carried out on the region to obtain a three-dimensional data volume of longitudinal wave impedance and transverse wave impedance. And further putting the two data together for comprehensive interpretation by establishing a quantitative interpretation quantity plate to obtain quantitative interpretation results of different porosities.

Wherein, the method comprises the following steps:

1) acquiring low-frequency condition elastic parameter test results, medium-frequency condition elastic parameter test results and high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the low frequency is consistent with the seismic main frequency and is 20Hz, the medium frequency is consistent with the logging frequency and is 10000Hz, and the high frequency is consistent with the ultrasonic testing frequency and is 100000 Hz; the elastic parameters comprise the velocity ratio of longitudinal waves to transverse waves and the impedance of longitudinal waves; specifically, the method comprises the following steps:

1.1) obtaining a work area target layer rock core;

1.2) carrying out sample preparation treatment and decontamination treatment on the rock core; wherein the decontamination treatment is carried out by soaking with an organic solution; the sample preparation treatment comprises the steps of processing the core into a sample with the diameter of 3.8 cm and the length of 5 cm, and polishing the surface of the core;

1.3) determining the porosity of the rock core;

1.4) respectively carrying out elastic parameter tests on the cores with different porosities under different oil saturation states under the condition of a low-frequency seismic source to obtain low-frequency elastic parameter test results of the cores with different porosities under different oil saturations; specifically, the method comprises the following steps:

calculating the amount of fluid filled in the pores of the rock core according to the numerical value of the total porosity, and controlling the change of the saturation; keeping a low-frequency seismic source consistent with the seismic main frequency (20HZ), carrying out oil-containing fluid displacement on the rock core (firstly carrying out water saturation, and then carrying out fluid displacement by using white oil), carrying out laboratory low-frequency condition elastic parameter tests every 10%, namely testing the low-frequency condition elastic parameters of each rock core under 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90% and 100% oil saturation, recording the longitudinal and transverse wave velocity ratio and longitudinal wave impedance data, and finishing the test results;

1.5) changing the seismic source condition, and repeating the step 1.4) to obtain the intermediate frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees;

1.6) changing the seismic source condition, and repeating the step 1.4) to obtain the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

2) The method comprises the following steps of firstly, casting a low-frequency condition elastic parameter test result of rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal and transverse wave velocity ratio and longitudinal wave impedance; supplementing the throw-point values of the rock cores with different porosities in the coordinate system under different oil saturation degrees by adopting an extrapolation-interpolation method, wherein the extrapolation-interpolation method adopts a first rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; fitting a relation trend line of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios of different oil saturation degrees of each porosity core under a low-frequency condition, and fitting a relation trend line of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios of different porosity cores under different oil saturation degrees under a low-frequency condition; thereby establishing a query plate of longitudinal and transverse wave parameter ratio and longitudinal wave impedance values of cores with different porosities and different oil saturation under the low-frequency condition, namely a low-frequency interpretation plate of the carbonate reservoir, and the result is shown in fig. 7;

the method comprises the steps that (1) the medium-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of longitudinal and transverse wave velocity ratio and longitudinal wave impedance; adding the throw-point values of the rock cores with different porosities in the coordinate system under different oil saturation degrees by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a second rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; (ii) a Fitting a relation trend line of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios of different oil saturation degrees of each porosity core under the medium-frequency condition, and fitting a relation trend line of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios of different porosity cores under the medium-frequency condition; thereby establishing a query plate of longitudinal and transverse wave parameter ratios and longitudinal wave impedance values of cores with different porosities and different oil saturation under the medium-frequency condition, namely a carbonate reservoir medium-frequency interpretation plate, and the result is shown in fig. 8;

the method comprises the steps of (1) casting high-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave velocity ratio and longitudinal wave impedance; adding the throw-point values of the rock cores with different porosities in the coordinate system under different oil saturation degrees by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a third rock physical theory model to determine the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; fitting a relation trend line of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios of different oil saturation degrees of the cores with different porosities under a high-frequency condition, and fitting a relation trend line of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios of the cores with different oil saturation degrees under the high-frequency condition; thereby establishing a query plate of longitudinal and transverse wave parameter ratio and longitudinal wave impedance values of cores with different porosities and different oil saturation under a high-frequency condition, namely a high-frequency interpretation plate of the carbonate reservoir, and the result is shown in fig. 9;

projecting the longitudinal wave velocity ratio (equal to the longitudinal wave impedance ratio) value and the longitudinal wave impedance value of each rock core experimental test result into a coordinate system with the y-axis being the longitudinal wave velocity ratio (Vp/Vs) and the x-axis being the longitudinal wave impedance (rho. Vp); because the porosity of each core is different, a change curve of the saturation under different porosities can be established; five rock cores are selected, the porosity degrees are respectively 1%, 2.3%, 3.5%, 4.8% and 5.2%, the maximum porosity degree of the test is 5.2%, so that the establishment of a larger porosity measurement version is supplemented by rock physical modeling, and if more rock core samples exist, a laboratory test interpretation version with a wider porosity degree range can be obtained;

in fig. 7, 8, and 9, the setpoints at different saturations for porosity of 1%, 2.3%, 3.5%, 4.8%, and 5.2% were based on the longitudinal-to-transverse wave ratio and longitudinal wave impedance data obtained from 5 core laboratory tests; due to the limitation of the porosity distribution of the sample, the rock theory physical model is adopted to calculate the longitudinal and transverse wave ratio value and the longitudinal wave impedance data value under different saturation degrees corresponding to the porosity of 10%, 20% and 30%, and an extrapolation interpolation method is used to perform supplementary dotting on the coordinate axis as the supplement of the measurement version; the used first rock theoretical physical model is over-calibrated by using the longitudinal and transverse wave ratio and the longitudinal wave impedance data under the low-frequency condition obtained by 5 rock core laboratory tests, the used second rock theoretical physical model is over-calibrated by using the longitudinal and transverse wave ratio and the longitudinal wave impedance data under the seed-frequency condition obtained by 5 rock core laboratory tests, and the used third rock theoretical physical model is over-calibrated by using the longitudinal and transverse wave ratio and the longitudinal wave impedance data under the high-frequency condition obtained by 5 rock core laboratory tests. If the rock samples are more and the porosity range distribution is larger, the point is added without using an extrapolation and interpolation method;

FIG. 10 is a comparison graph of the overlapping of the quantitative interpretation quantity versions in different frequency bands, and it can be seen from the graph that the interpretation quantity versions have differences at different frequencies, and a certain data point has different readings in different quantity versions, which illustrates the necessity of different data interpretation by using different quantity versions; only if the number is matched, a more accurate interpretation result can be obtained.

3) And performing carbonate reservoir interpretation based on the pre-stack seismic inversion result:

3.1) obtaining a target prestack seismic inversion result of a work area, and determining a first longitudinal-transverse wave velocity ratio data volume and a first longitudinal wave impedance data volume according to the prestack seismic inversion result;

utilizing a longitudinal wave impedance data volume and a transverse wave impedance data volume obtained by seismic inversion before the target stack of a work area; calculating a three-dimensional data volume of a longitudinal-transverse wave velocity ratio Vp/Vs by using longitudinal wave impedance ÷ transverse wave impedance;

3.2) intersecting the first longitudinal wave impedance data body by using the first longitudinal-transverse wave velocity ratio data body, and quantitatively interpreting the porosity and the saturation by referring to a low-frequency quantitative interpretation quantity version to finish carbonate reservoir interpretation; specifically, the method comprises the following steps:

carrying out two-dimensional coordinate intersection on the first longitudinal wave impedance data volume and the first longitudinal wave velocity ratio data volume, and projecting the record of each sampling point in the space into a coordinate system; and projecting the established low-frequency quantitative interpretation version in a coordinate system, expressing different porosity distribution ranges by using different colors according to 'marks' of the interpretation version, selecting sample points of the data body of the effective oil-containing (gas-containing) reservoir region corresponding to the sample points, and displaying the sample points in a three-dimensional data body, so that the spatial distribution range of the effective porosity of the carbonate reservoir can be predicted.

4) And carrying out carbonate reservoir interpretation based on the logging data:

4.1) obtaining logging data of a target layer of a work area, and determining a second longitudinal wave velocity ratio data volume and a second longitudinal wave impedance data volume according to logging;

determining a three-dimensional data body of a longitudinal wave velocity ratio Vp/Vs and a longitudinal wave impedance data body by utilizing a transverse wave logging curve, a longitudinal wave logging curve and a density logging curve obtained by logging in a target zone of a work area and utilizing the transverse wave logging curve, the longitudinal wave logging curve and the density logging curve;

4.2) intersecting the second longitudinal wave impedance data body by using the second longitudinal-transverse wave velocity ratio data body, and quantitatively interpreting the porosity and the saturation by referring to a medium-frequency quantitative interpretation quantity version to finish carbonate reservoir interpretation; specifically, the method comprises the following steps:

carrying out two-dimensional coordinate intersection on a second longitudinal wave impedance data volume and a second longitudinal wave velocity ratio data volume, and projecting the record of each sampling point in the space into a coordinate system; the established intermediate frequency quantitative interpretation version is projected in a coordinate system, different porosity distribution ranges are represented by using different colors according to 'marks' of the interpretation version, effective oil-containing (gas-containing) reservoir region data body sampling points corresponding to the intermediate frequency quantitative interpretation version are selected, and the sampling points are displayed in a three-dimensional data body, so that the effective porosity spatial distribution range of the carbonate reservoir can be predicted.

5) Carrying out carbonate reservoir interpretation based on formation ultrasonic testing data:

5.1) acquiring ultrasonic test data of a target layer of a work area, and determining a third longitudinal-transverse wave velocity ratio data volume and a third longitudinal wave impedance data volume according to the ultrasonic test data;

calculating a three-dimensional data volume of a longitudinal wave velocity ratio Vp/Vs and a longitudinal wave impedance data volume by using data obtained by ultrasonic testing of a target layer of a work area;

5.2) intersecting the third longitudinal wave impedance data body by using a third longitudinal-transverse wave velocity ratio data body, and quantitatively interpreting porosity and saturation by referring to a high-frequency quantitative interpretation quantity version to finish carbonate reservoir interpretation; specifically, the method comprises the following steps:

carrying out two-dimensional coordinate intersection on a third longitudinal wave impedance data volume and a third longitudinal wave velocity ratio data volume, and projecting the record of each sampling point in the space to a coordinate system; the established high-frequency quantitative interpretation version is projected in a coordinate system, different porosity distribution ranges are represented by using different colors according to 'marks' of the interpretation version, effective oil-containing (gas-containing) reservoir region data body sampling points corresponding to the high-frequency quantitative interpretation version are selected, and the sampling points are displayed in a three-dimensional data body, so that the effective porosity spatial distribution range of the carbonate reservoir can be predicted.

The embodiment of the invention also provides a system for establishing the carbonate reservoir interpretation volume version, and the system is preferably used for realizing the method embodiment.

Fig. 11 is a block diagram of a system for establishing an interpretation measure of a carbonate reservoir according to an embodiment of the present invention, as shown in fig. 11, the system including: an elasticity parameter obtaining module 111 and an interpretation quantity version establishing module 112.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the low frequency is consistent with the seismic dominant frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method is used for casting the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the low-frequency condition, namely a carbonate reservoir low-frequency interpretation quantity chart.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the intermediate frequency is consistent with the logging frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method is used for casting the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance, and establishing a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition, namely a carbonate reservoir medium-frequency interpretation plate.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely a carbonate reservoir high-frequency interpretation quantity chart.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein, the low frequency is consistent with the earthquake main frequency, and the intermediate frequency is consistent with the logging frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; the method is used for casting the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance, and establishing a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition, namely a carbonate reservoir medium-frequency interpretation plate.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the low frequency is consistent with the seismic main frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; the method is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely a carbonate reservoir high-frequency interpretation quantity chart.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method comprises the steps that a test result of medium-frequency elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate; the method is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely a carbonate reservoir high-frequency interpretation quantity chart.

In one embodiment:

the elastic parameter obtaining module 111: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the low frequency is consistent with the earthquake main frequency, the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 112: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; the method comprises the steps that a test result of medium-frequency elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate; the method is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely a carbonate reservoir high-frequency interpretation quantity chart.

In an embodiment, the elasticity parameter obtaining module 111 may include:

in an embodiment, the elasticity parameter obtaining module 111 further includes:

In an embodiment, the interpretation quantity plate establishing module 112 is configured to throw the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedances, and when an inquiry plate of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities and different oil saturation degrees under the low-frequency condition is the low-frequency interpretation quantity plate of the carbonate reservoir, the interpretation quantity plate establishing module 112 is further configured to throw the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedances, and further supplement the throw values of the cores with different porosities under different oil saturation degrees in the coordinate system by using an extrapolation and interpolation method, where the extrapolation and interpolation method uses a first rock physical theory model to perform transverse wave parameter test on the cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; the first rock physical theory model is calibrated by using low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

in an embodiment, the interpretation quantity plate establishing module 112 is configured to throw the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal-transverse wave parameter ratios and longitudinal wave impedances under different oil saturation degrees, and when an inquiry plate of the longitudinal-transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities and different oil saturation degrees under the intermediate-frequency condition is the interpretation quantity plate under the intermediate-frequency condition of the carbonate reservoir, the interpretation quantity plate establishing module 112 is further configured to throw the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of longitudinal-transverse wave parameter ratios and longitudinal wave impedances under different oil saturation degrees by using an extrapolation and interpolation method, wherein the extrapolation and interpolation method uses a second rock physical theory model to perform transverse-wave parameter test on the cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

in an embodiment, the interpretation quantity plate establishing module 112 is configured to throw the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and when an inquiry plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is the interpretation quantity plate of the carbonate reservoir under the high frequency, the interpretation quantity plate establishing module 112 is further configured to throw the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and further supplement the throw values of the cores with different porosities under different oil saturation degrees in the coordinate system by using an extrapolation and interpolation method, wherein the extrapolation and interpolation method uses a third rock physical theory model to perform transverse wave parameter test on the cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

In an embodiment, the interpretation quantity plate establishing module 112 is configured to cast the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and when an inquiry plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the low-frequency condition is the interpretation quantity plate under the low-frequency of the carbonate reservoir, the interpretation quantity plate establishing module 112 is further configured to cast the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance, and then fit a relationship trend line of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance ratio of the cores with different oil saturation degrees under the low-frequency condition and/or fit a relationship of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance ratio of the cores with different porosity degrees under the low-oil saturation degrees under the low-frequency condition A potential line;

in an embodiment, the interpretation quantity plate establishing module 112 is configured to cast the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and when an inquiry plate of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities and different oil saturation degrees under the intermediate-frequency condition is the interpretation quantity plate under the intermediate frequency of the carbonate reservoir, the interpretation quantity plate establishing module 112 is further configured to cast the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance ratios, and then fit a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance ratios of the cores with different oil saturation degrees under the intermediate frequency condition and/or fit a relationship of the longitudinal and transverse wave parameter ratios of the cores with different porosity degrees under the intermediate frequency condition A potential line;

in an embodiment, the interpretation quantity plate establishing module 112 is configured to throw the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and when an inquiry plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is the interpretation quantity plate of the carbonate reservoir under the high frequency, the interpretation quantity plate establishing module 112 is further configured to throw the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and then fit a relationship trend line of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance ratio of the cores with different oil saturation degrees under the high frequency condition and/or fit a relationship of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance ratio of the cores with different porosity degrees under the high frequency condition A potential line.

The embodiment of the invention also provides a carbonate reservoir interpretation system which is preferably used for realizing the method embodiment.

Fig. 12 is a block diagram of a carbonate reservoir interpretation system according to an embodiment of the present invention, as shown in fig. 12, the system including: the system comprises an elastic parameter acquisition module 121, an interpretation quantity version establishing module 122, a work area target layer data acquisition module 123 and a reservoir quantitative interpretation module 124.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the low frequency is consistent with the seismic dominant frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate;

the work area target layer data acquisition module 123: the method comprises the steps of obtaining a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the pre-stack seismic inversion result;

reservoir quantitative interpretation module 124: the method is used for conducting intersection by utilizing the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body determined according to the pre-stack seismic inversion result, and conducting porosity and/or saturation quantitative interpretation by referring to the interpretation quantity version under the low frequency.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the intermediate frequency is consistent with the logging frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of medium-frequency elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate;

the work area target layer data acquisition module 123: the system comprises a data acquisition module, a data acquisition module and a data acquisition module, wherein the data acquisition module is used for acquiring target logging data of a work area and determining a longitudinal and transverse wave parameter ratio data volume and a longitudinal wave impedance data volume according to the logging data;

reservoir quantitative interpretation module 124: the method is used for performing intersection by utilizing the longitudinal and transverse wave parameters, the data body and the longitudinal wave impedance data body, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation quantity version under the medium frequency.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of high-frequency condition elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is established, namely a carbonate reservoir high-frequency interpretation quantity chart;

the work area target layer data acquisition module 123: the system comprises a data volume, a longitudinal wave impedance data volume and a transverse wave parameter ratio data volume, wherein the data volume is used for acquiring stratum ultrasonic test data of a work area and determining a longitudinal wave parameter ratio data volume and a longitudinal wave impedance data volume according to the stratum ultrasonic test data;

reservoir quantitative interpretation module 124: the method is used for carrying out intersection by utilizing the longitudinal and transverse wave parameters to compare with the data body and the longitudinal wave impedance data body, and carrying out porosity and/or saturation quantitative interpretation by referring to the interpretation amount version under high frequency.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; wherein, the low frequency is consistent with the earthquake main frequency, and the intermediate frequency is consistent with the logging frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; the method comprises the steps that a test result of medium-frequency elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate;

the work area target layer data acquisition module 123: the method comprises the steps of obtaining a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the pre-stack seismic inversion result; the system comprises a data acquisition module, a data acquisition module and a data acquisition module, wherein the data acquisition module is used for acquiring target logging data of a work area and determining a longitudinal and transverse wave parameter ratio data volume and a longitudinal wave impedance data volume according to the logging data;

reservoir quantitative interpretation module 124: the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body determined according to a pre-stack seismic inversion result, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation version under low frequency; the method is used for performing intersection by utilizing the longitudinal and transverse wave parameters, the data body and the longitudinal wave impedance data body, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation quantity version under the medium frequency.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the low frequency is consistent with the seismic main frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; the method comprises the steps that a test result of high-frequency condition elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is established, namely a carbonate reservoir high-frequency interpretation quantity chart;

the work area target layer data acquisition module 123: the method comprises the steps of obtaining a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the pre-stack seismic inversion result; the system comprises a data volume, a longitudinal wave impedance data volume and a transverse wave parameter ratio data volume, wherein the data volume is used for acquiring stratum ultrasonic test data of a work area and determining a longitudinal wave parameter ratio data volume and a longitudinal wave impedance data volume according to the stratum ultrasonic test data;

reservoir quantitative interpretation module 124: the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body determined according to a pre-stack seismic inversion result, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation version under low frequency; the method is used for carrying out intersection by utilizing the longitudinal and transverse wave parameters to compare with the data body and the longitudinal wave impedance data body, and carrying out porosity and/or saturation quantitative interpretation by referring to the interpretation amount version under high frequency.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of medium-frequency elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate; the method comprises the steps that a test result of high-frequency condition elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is established, namely a carbonate reservoir high-frequency interpretation quantity chart;

the work area target layer data acquisition module 123: the system comprises a data acquisition module, a data acquisition module and a data acquisition module, wherein the data acquisition module is used for acquiring target logging data of a work area and determining a longitudinal and transverse wave parameter ratio data volume and a longitudinal wave impedance data volume according to the logging data; the system comprises a data volume, a longitudinal wave impedance data volume and a transverse wave parameter ratio data volume, wherein the data volume is used for acquiring stratum ultrasonic test data of a work area and determining a longitudinal wave parameter ratio data volume and a longitudinal wave impedance data volume according to the stratum ultrasonic test data;

reservoir quantitative interpretation module 124: the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation quantity version under intermediate frequency; the method is used for carrying out intersection by utilizing the longitudinal and transverse wave parameters to compare with the data body and the longitudinal wave impedance data body, and carrying out porosity and/or saturation quantitative interpretation by referring to the interpretation amount version under high frequency.

In one embodiment:

the elastic parameter obtaining module 121: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining medium-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the method comprises the steps of obtaining high-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees; the low frequency is consistent with the earthquake main frequency, the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave parameter ratio and a transverse wave impedance, and the longitudinal wave parameter ratio is a longitudinal wave impedance ratio or a longitudinal wave velocity ratio;

interpretation volume creation module 122: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions is established, namely a carbonate reservoir low-frequency interpretation plate; the method comprises the steps that a test result of medium-frequency elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and a query plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely a carbonate reservoir medium-frequency interpretation quantity plate; the method comprises the steps that a test result of high-frequency condition elastic parameters of cores with different porosities under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is established, namely a carbonate reservoir high-frequency interpretation quantity chart;

work area target layer data acquisition module: the method comprises the steps of obtaining a pre-stack seismic inversion result of a work area target, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the pre-stack seismic inversion result; the system comprises a data acquisition module, a data acquisition module and a data acquisition module, wherein the data acquisition module is used for acquiring target logging data of a work area and determining a longitudinal and transverse wave parameter ratio data volume and a longitudinal wave impedance data volume according to the logging data; the system comprises a data volume, a longitudinal wave impedance data volume and a transverse wave parameter ratio data volume, wherein the data volume is used for acquiring stratum ultrasonic test data of a work area and determining a longitudinal wave parameter ratio data volume and a longitudinal wave impedance data volume according to the stratum ultrasonic test data;

reservoir quantitative interpretation module: the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body determined according to a pre-stack seismic inversion result, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation version under low frequency; the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation quantity version under intermediate frequency; the method is used for carrying out intersection by utilizing the longitudinal and transverse wave parameters to compare with the data body and the longitudinal wave impedance data body, and carrying out porosity and/or saturation quantitative interpretation by referring to the interpretation amount version under high frequency.

In an embodiment, the elasticity parameter obtaining module 121 may include:

in an embodiment, the elasticity parameter obtaining module 121 further includes:

In an embodiment, the interpretation quantity plate establishing module 122 is configured to throw the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedances, and when an inquiry plate of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities and different oil saturation degrees under the low-frequency condition is the low-frequency interpretation quantity plate of the carbonate reservoir, the interpretation quantity plate establishing module 122 is further configured to throw the low-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedances, and further supplement the throw values of the cores with different porosities under different oil saturation degrees in the coordinate system by using an extrapolation and interpolation method, where the extrapolation and interpolation method uses a first rock physical theory model to perform transverse wave parameter test on the cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; the first rock physical theory model is calibrated by using low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

in an embodiment, the interpretation quantity plate establishing module 122 is configured to throw the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of longitudinal-transverse wave parameter ratios and longitudinal wave impedances under different oil saturation degrees, and when an inquiry plate of the longitudinal-transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosities and different oil saturation degrees under the intermediate-frequency condition is the interpretation quantity plate under the intermediate-frequency condition of the carbonate reservoir, the interpretation quantity plate establishing module 122 is further configured to throw the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of longitudinal-transverse wave parameter ratios and longitudinal wave impedances under different oil saturation degrees by using an extrapolation interpolation method, wherein the extrapolation interpolation method uses a second rock physical theory model to perform transverse-wave parameter test on the cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

in an embodiment, the interpretation quantity plate establishing module 122 is configured to throw the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into a coordinate system with coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and when an inquiry plate of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is the interpretation quantity plate of the carbonate reservoir under the high frequency, the interpretation quantity plate establishing module 122 is further configured to throw the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees into the coordinate system with coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and further supplement the throw-point values of the cores with different porosities under different oil saturation degrees in the coordinate system by using an extrapolation interpolation method, wherein the extrapolation interpolation method uses a third rock physical theory model to perform the transverse wave parameter test on the longitudinal wave parameters of the cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

In an embodiment, the interpretation quantity plate establishing module 122 is configured to, when the test results of the low-frequency condition elastic parameters of the cores with different porosities under different oil saturation degrees are cast into the coordinate system with the coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and the query chart of the longitudinal wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the low-frequency condition is the interpretation quantity plate under the low-frequency condition of the carbonate reservoir, the interpretation quantity plate establishing module 122 is further configured to, after the test results of the low-frequency condition elastic parameters of the cores with different porosities under different oil saturation degrees are cast into the coordinate system with the coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, fit a trend line of the relationship between the longitudinal wave parameter ratio and the longitudinal wave impedance ratio of the cores with different oil saturation degrees under the low-frequency condition and/or fit a trend line of the relationship between the longitudinal wave parameter ratio and the longitudinal wave impedance ratio of the cores with different porosity degrees under the oil saturation degrees under the low-frequency condition A potential line;

in an embodiment, the interpretation quantity plate establishing module 122 is configured to, when the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees are cast into the coordinate system with the coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and the query chart of the longitudinal wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the intermediate-frequency condition is established as the interpretation quantity plate under the intermediate frequency of the carbonate reservoir, the interpretation quantity plate establishing module 122 is further configured to, after the intermediate-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees are cast into the coordinate system with the coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, fit a trend line of a relationship between the longitudinal wave parameter ratio and the longitudinal wave impedance ratio of the cores with different oil saturation degrees under the intermediate frequency condition and/or fit a trend of a relationship between the longitudinal wave parameter ratio and the longitudinal wave impedance ratio of the cores with different porosity degrees under the intermediate frequency condition A potential line;

in an embodiment, the interpretation quantity plate establishing module 122 is configured to, when the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees are cast into the coordinate system with the coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, and the query chart of the longitudinal wave parameter ratio and the longitudinal wave impedance value of the cores with different porosities and different oil saturation degrees under the high-frequency condition is established as the interpretation quantity plate of the carbonate reservoir under the high frequency, the interpretation quantity plate establishing module 122 is further configured to, after the high-frequency condition elastic parameter test results of the cores with different porosities under different oil saturation degrees are cast into the coordinate system with the coordinate axes of the longitudinal wave parameter ratio and the longitudinal wave impedance, fit a trend line of the relationship between the longitudinal wave parameter ratio and the longitudinal wave impedance ratio of the cores with different oil saturation degrees under the high frequency condition and/or fit a trend line of the relationship between the longitudinal wave parameter ratio and the longitudinal wave impedance ratio of the cores with different porosity degrees under the oil saturation degrees under the high frequency condition A potential line.

In one embodiment, the pre-stack seismic inversion results include a compressional impedance data volume and a shear impedance data volume.

In one embodiment, the well log data includes shear wave well logs, compressional wave well logs, and density well logs.

Fig. 13 is a schematic diagram of a device for establishing an explanatory volume version of a carbonate reservoir according to an embodiment of the invention. The device for creating the interpretation version of the carbonate reservoir shown in fig. 13 is a general-purpose data processing device, which includes a general-purpose computer hardware structure, which includes at least a processor 1000 and a memory 1111; the processor 1000 is configured to execute the carbonate reservoir interpretation quantity version establishing program stored in the memory, so as to implement the carbonate reservoir interpretation quantity version establishing method according to each method embodiment (for a specific method, refer to the description of the method embodiment, and no further description is given here).

One embodiment of the invention provides a schematic diagram of a carbonate reservoir interpretation device, which comprises a general computer hardware structure, at least comprising a processor and a memory; the processor is configured to execute the carbonate reservoir interpretation program stored in the memory to implement the carbonate reservoir interpretation method according to each method embodiment (for a specific method, refer to the description of the method embodiment, and no further description is given here).

The embodiment of the present invention further provides a computer-readable storage medium, where one or more programs are stored, and the one or more programs may be executed by one or more processors to implement the method for establishing the carbonate reservoir interpretation volume version in each method embodiment (for a specific method, refer to the description of the above method embodiment, and are not described herein again).

The embodiment of the present invention further provides a computer-readable storage medium, where one or more programs are stored, and the one or more programs may be executed by one or more processors to implement the carbonate reservoir interpretation method described in each method embodiment (for a specific method, refer to the description of the above method embodiment, and are not described herein again).

The preferred embodiments of the present invention have been described above with reference to the accompanying drawings. The many features and advantages of the embodiments are apparent from the detailed specification, and thus, it is intended by the appended claims to cover all such features and advantages of the embodiments which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the embodiments of the invention to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope thereof.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on at least one computer-usable storage medium (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims

1. A method for establishing carbonate reservoir interpretation volume edition, wherein the carbonate reservoir interpretation volume edition comprises carbonate reservoir interpretation volume editions under low frequency and/or medium frequency and/or high frequency, the establishing method comprises the following steps:

the low frequency is consistent with the seismic main frequency, the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency;

2. The setup method according to claim 1,

the elastic parameter test results of the rock cores with different porosities of the target layer of the work area under the low-frequency and/or medium-frequency and/or high-frequency conditions with different oil saturation degrees are obtained by the following method:

acquiring a work area target layer rock core;

determining the porosity of a target layer rock core in a work area;

3. The building method according to claim 2, wherein the obtained work area target layer core is subjected to pretreatment before the porosity is determined, and the pretreatment comprises sample preparation treatment and decontamination treatment.

4. The setup method according to claim 1,

when the interpretation version of the carbonate reservoir under the low frequency is established, the method further comprises the following steps: after the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal wave parameter ratio and a transverse wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation and interpolation method, wherein the extrapolation and interpolation method adopts a first rock physical theory model to determine the longitudinal and transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities under different oil saturation degrees; the first rock physical theory model is calibrated by using low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

when establishing the interpretation quantity version under the frequency of the carbonate reservoir, the method further comprises the following steps: after the test results of the medium-frequency elastic parameters of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axes of longitudinal-transverse wave parameter ratios and longitudinal wave impedance, further supplementing the cast values of the rock cores with different porosities under different oil saturation degrees in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a second rock physical theory model to determine the longitudinal-transverse wave parameter ratios and the longitudinal wave impedance values of the rock cores with different porosities under different oil saturation degrees; the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

when the interpretation version of the carbonate reservoir under high frequency is established, the method further comprises the following steps: after the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees are cast into a coordinate system with the coordinate axis of a longitudinal wave parameter ratio and a transverse wave impedance, the cast values of the rock cores with different porosities under different oil saturation degrees are further supplemented in the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a third rock physical theory model to determine the longitudinal wave parameter ratio and the transverse wave impedance value of the rock cores with different porosities under different oil saturation degrees; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

5. The establishment method of claim 1 or 4,

establishing a carbonate reservoir low-frequency interpretation version, the method further comprising: after the elastic parameter test results of the cores with different porosities under the low-frequency condition under different oil saturation degrees are cast into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the low-frequency condition and/or fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the low-frequency condition;

establishing a lower frequency interpretation version of the carbonate reservoir, the method further comprising: the method comprises the steps of putting test results of elastic parameters of rock cores with different porosities under medium-frequency conditions with different oil saturation into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and fitting a relation trend line of the longitudinal and transverse wave parameter ratios with different oil saturations of the rock cores with different porosities and the longitudinal wave impedance value under the medium-frequency condition and/or fitting a relation trend line of the longitudinal and transverse wave parameter ratios with the longitudinal wave impedance value of the rock cores with different porosities under the medium-frequency condition;

establishing a carbonate reservoir interpretation version at high frequency, the method further comprising: the method comprises the steps of putting test results of elastic parameters of cores with different porosities under high-frequency conditions under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, and fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different oil saturation degrees under the high-frequency condition and/or fitting a relationship trend line of the longitudinal and transverse wave parameter ratios and the longitudinal wave impedance values of the cores with different porosity degrees under the high-frequency condition.

6. A carbonate reservoir interpretation method, wherein the method comprises:

wherein the carbonate reservoir interpretation based on the prestack seismic inversion result comprises: establishing a low-frequency interpretation quantity version of a target carbonate reservoir in a work area by adopting the method for establishing the interpretation quantity version of the carbonate reservoir in any one of claims 1 to 5;

intersecting the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation quantity version under low frequency;

establishing a frequency lower interpretation quantity version of a work area target layer carbonate reservoir by adopting the establishment method of the carbonate reservoir interpretation quantity version of any one of claims 1-5;

intersecting the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body, and quantitatively interpreting the porosity and/or the saturation by referring to an interpretation quantity version under intermediate frequency;

wherein the carbonate reservoir interpretation based on formation ultrasonic test data comprises:

establishing a high-frequency interpretation quantity version of a target carbonate reservoir in a work area by adopting the method for establishing the interpretation quantity version of the carbonate reservoir in any one of claims 1 to 5;

acquiring stratum ultrasonic test data of a work area, and determining a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body according to the stratum ultrasonic test data;

7. The interpretation method according to claim 6, wherein,

the pre-stack seismic inversion result comprises a longitudinal wave impedance data volume and a transverse wave impedance data volume;

the logging data comprises a transverse wave logging curve, a longitudinal wave logging curve and a density logging curve;

the formation ultrasonic testing data comprises longitudinal wave data, transverse wave data and density data.

8. A system for establishing an interpretation version of a carbonate reservoir, wherein the system comprises:

an elastic parameter acquisition module: the method comprises the steps of obtaining low-frequency condition elastic parameter test results of rock cores with different porosities of a target layer of a work area under different oil saturation degrees, and/or obtaining medium-frequency condition elastic parameter test results of the rock cores with different porosities of the target layer of the work area under different oil saturation degrees, and/or obtaining high-frequency condition elastic parameter test results of the rock cores with different porosities of the target layer of the work area under different oil saturation degrees; the low frequency is consistent with the seismic main frequency, the medium frequency is consistent with the logging frequency, and the high frequency is consistent with the formation ultrasonic testing frequency; the elastic parameters comprise a longitudinal wave and transverse wave parameter ratio and longitudinal wave impedance, and the longitudinal wave and transverse wave parameter ratio is a longitudinal wave and transverse wave impedance ratio or a longitudinal wave and transverse wave velocity ratio;

an interpretation quantity version establishing module: the method comprises the steps that a test result of elastic parameters of cores with different porosities under low-frequency conditions under different oil saturation degrees is cast into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and query charts of longitudinal-transverse wave parameter ratio and longitudinal wave impedance value of cores with different porosities and different oil saturation degrees under the low-frequency conditions are established, namely a carbonate reservoir low-frequency interpretation quantity chart; and/or, the method is used for casting the test results of the elastic parameters of the rock cores with different porosities under the medium-frequency condition under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing query charts of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition, namely the interpretation quantity chart of the rock carbonate reservoir under the medium frequency; and/or, the method is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis as the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely the carbonate reservoir high-frequency interpretation quantity chart.

9. The setup system of claim 8, wherein the elasticity parameter acquisition module comprises:

an elastic parameter acquisition submodule: the method is used for respectively testing the elastic parameters of the rock cores with different porosities under different oil saturation states under the condition of a low-frequency seismic source to obtain the low-frequency elastic parameter test results of the rock cores with different porosities under different oil saturations; and/or, the elastic parameter testing device is used for respectively testing the elastic parameters of the rock cores with different porosities under different oil saturation states under the condition of the medium-frequency seismic source to obtain the elastic parameter testing results of the rock cores with different porosities under the medium-frequency conditions under different oil saturations; and/or the elastic parameter testing device is used for respectively testing the elastic parameters of the rock cores with different porosities under different oil saturation states under the condition of a high-frequency seismic source to obtain the high-frequency elastic parameter testing results of the rock cores with different porosities under different oil saturations.

10. The setup system of claim 9, wherein the elasticity parameter acquisition module further comprises:

11. The build system of claim 8,

the interpretation quantity version establishing module is used for casting the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the low-frequency condition, namely the carbonate reservoir low-frequency interpretation quantity version, and the interpretation quantity version establishing module is further used for casting the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into the coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and then further supplementing the casting values of the rock cores with different porosities under different oil saturation degrees into the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a first rock physical theory model to perform longitudinal-transverse wave parameter test on the rock cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; the first rock physical theory model is calibrated by using low-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

the interpretation quantity version establishing module is used for casting the intermediate frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the intermediate frequency condition, namely the interpretation quantity version of the carbonate reservoir under the intermediate frequency, and the interpretation quantity version establishing module is further used for further supplementing the casting values of the rock cores with different porosities under different oil saturation degrees into the coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance by adopting an extrapolation and interpolation method after casting the intermediate frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into the coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, wherein the extrapolation and interpolation method adopts a second rock physical theory model to perform longitudinal-transverse wave parameter test on the rock cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; the second rock physical theory model is calibrated by using the medium-frequency condition elastic parameter test results of rock cores with different porosities under different oil saturation degrees;

the interpretation quantity version establishing module is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and establishing an inquiry chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition, namely the carbonate reservoir high-frequency interpretation quantity version, and the interpretation quantity version establishing module is further used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into the coordinate system with the coordinate axis of longitudinal-transverse wave parameter ratio and longitudinal wave impedance, and then further supplementing the casting values of the rock cores with different porosities under different oil saturation degrees into the coordinate system by adopting an extrapolation interpolation method, wherein the extrapolation interpolation method adopts a third rock physical theory model to perform longitudinal-transverse wave parameter test on the rock cores with different porosities under different oil saturation degrees Determining a numerical ratio and a longitudinal wave impedance value; and the third rock physical theory model is calibrated by using the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees.

12. The building system according to claim 8 or 11,

the interpretation quantity plate establishing module is used for casting the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, when the query chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock core with different porosities and different oil saturation under the low-frequency condition is established as the interpretation quantity chart of the carbonate reservoir under the low-frequency condition, the interpretation quantity plate establishing module is further used for casting the low-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal and transverse wave parameter ratio and longitudinal wave impedance, fitting a relation trend line of longitudinal and transverse wave parameter ratios of different oil saturation degrees of each porosity rock core and longitudinal wave impedance values under a low-frequency condition and/or fitting a relation trend line of longitudinal and transverse wave parameter ratios of different porosity rock cores and longitudinal wave impedance values under different oil saturation degrees of each porosity rock core under a low-frequency condition;

the interpretation quantity plate establishing module is used for casting the test results of the medium-frequency condition elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, when the query chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the medium-frequency condition is established, namely the interpretation quantity chart of the carbonate reservoir under the medium-frequency condition, the interpretation quantity plate establishing module is further used for casting the test results of the medium-frequency condition elastic parameters of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, fitting a relation trend line of longitudinal and transverse wave parameter ratios of different oil saturation degrees of each porosity rock core and longitudinal wave impedance values under the medium-frequency condition and/or fitting a relation trend line of longitudinal and transverse wave parameter ratios of different porosity rock cores and longitudinal wave impedance values under the medium-frequency condition;

the interpretation quantity plate establishing module is used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axes of longitudinal and transverse wave parameter ratios and longitudinal wave impedance, when the query chart of the longitudinal-transverse wave parameter ratio and the longitudinal wave impedance value of the rock cores with different porosities and different oil saturation degrees under the high-frequency condition is established, namely the interpretation chart of the carbonate reservoir under the high frequency is established, the interpretation quantity plate establishing module is further used for casting the high-frequency condition elastic parameter test results of the rock cores with different porosities under different oil saturation degrees into a coordinate system with the coordinate axis of longitudinal and transverse wave parameter ratio and longitudinal wave impedance, and fitting a relation trend line of the longitudinal wave impedance value and the longitudinal wave parameter ratio of the core with different oil saturation degrees under the high-frequency condition and/or fitting a relation trend line of the longitudinal wave impedance value and the longitudinal wave parameter ratio of the core with different oil saturation degrees under the high-frequency condition.

13. A carbonate reservoir interpretation system, wherein the system comprises:

a carbonate reservoir interpretation volume creation system as defined in any one of claims 8 to 12: the method is used for establishing a low-frequency interpretation quantity version of a target layer carbonate reservoir in a work area and/or a medium-frequency interpretation quantity version of the target layer carbonate reservoir and/or a high-frequency interpretation quantity version of the target layer carbonate reservoir;

reservoir quantitative interpretation module: the device is used for performing intersection by utilizing a longitudinal wave parameter ratio data body and a longitudinal wave impedance data body determined according to the pre-stack seismic inversion result, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation version under low frequency; and/or, the device is used for utilizing the longitudinal and transverse wave parameters to intersect the data body and the longitudinal wave impedance data body, and performing porosity and/or saturation quantitative interpretation by referring to an interpretation quantity version under the medium frequency; and/or, the device is used for carrying out intersection by utilizing the longitudinal wave parameter ratio data body and the longitudinal wave impedance data body, and carrying out porosity and/or saturation quantitative interpretation by referring to an interpretation amount version under high frequency.

14. The device for establishing the carbonate reservoir interpretation volume edition comprises a processor and a memory; wherein the content of the first and second substances,

a memory for storing a computer program;

a processor for implementing the steps of the method for establishing a carbonate reservoir interpretation volume according to any of claims 1 to 5 when executing the program stored on the memory.

15. A carbonate reservoir interpretation device comprises a processor and a memory; wherein the content of the first and second substances,

a memory for storing a computer program;

a processor for implementing the steps of the carbonate reservoir interpretation method of claim 6 or 7 when executing the program stored on the memory.

16. A computer readable storage medium storing one or more programs, the one or more programs being executable by one or more processors to perform the steps of the method for establishing a carbonate reservoir interpretation volume as defined in any one of claims 1 to 5.

17. A computer readable storage medium storing one or more programs, the one or more programs being executable by one or more processors to perform the steps of the carbonate reservoir interpretation method of claim 6 or 7.