A method and device for quantitatively evaluating the influence of dissolution on the elasticity of sedimentary rock

Abstract

The application relates to the field of oil and gas exploration and development and rock physics, and discloses a quantitative evaluation of the influence of dissolution on the elasticity of sedimentary rock. The method comprises the following steps: determining the mineral components of a carbonate rock sample of a target layer, the content of each mineral component, the first pore characteristics and the first elastic parameters; determining the dissolution environment parameters according to the sedimentary environment, the underground fluid condition and the burial history of the area where the carbonate rock sample is located, and performing dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameters; re-determining the second pore characteristics and the second elastic parameters of the carbonate rock sample; establishing a quantitative evaluation model according to the mineral components, the content of each mineral component, the first pore characteristics, the first elastic parameters, the second pore characteristics and the second elastic parameters, and evaluating the influence of dissolution on the elasticity of the carbonate rock according to the quantitative evaluation model. The application helps to improve the accuracy of seismic reservoir prediction.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil and gas exploration and development and petrophysics, and in particular to a method and device for quantitatively evaluating the influence of dissolution on the elasticity of sedimentary rock, a storage medium and an electronic device. BACKGROUND

[0002] The purpose of the background description provided herein is to generally present the context of the application. The statements in this section are not to be taken as admissions of prior art.

[0003] In recent years, with the progress of exploration technology, carbonate rock oil and gas exploration has gradually been pushed to deep layers. For example, several deep large gas fields such as Tahe, Puguang, Yuanba, Anyue and Shanzhong have been discovered in Tarim, Sichuan and Ordos basins.

[0004] However, due to deep burial and old age, carbonate rocks in China have been subjected to the influence of high temperature and high pressure, diagenesis and cementation for a long time, and there are few primary pores left. Secondary pores, fractures and caves are important oil and gas reservoir spaces, and dissolution is the key to the formation of secondary pores.

[0005] At present, in the research of carbonate reservoirs, scholars have carried out a lot of research work on petrology and reservoir performance evaluation. For example, some scholars believe that tectonic-pressure coupling controls fractures and dissolution, and fluid-rock interaction controls deep dissolution and pore preservation; some scholars point out that meteoric water leaching and dissolution during penecontemporaneous period is a necessary condition for the formation of large-scale high-quality reservoirs. For the research literature on dissolution, high-temperature and high-pressure dissolution kinetics simulation experimental devices are used to carry out experimental research on the burial dissolution temperature window and pore evolution pattern of carbonate rocks, mainly focusing on the formation and mechanism of secondary pores in the dissolution process.

[0006] After dissolution diagenesis, carbonate rocks develop secondary pores, and the pore network structure is complex, which seriously affects the elastic properties of carbonate reservoirs, and the rock velocity-porosity relationship shows strong dispersion. In the existing research on the influence of dissolution on the elastic properties of rocks, a large number of rock samples are collected to study the velocity variation law during the dissolution process. For example, some scholars collected multiple core samples of different ages from different regions, measured the P-wave and S-wave velocities in the laboratory, and analyzed the influence of dissolution on rock velocity; some scholars selected multiple dense limestone samples for ultrasonic velocity measurement to analyze the influence of pore structure on the elastic characteristics of rock skeleton and fluid sensitivity.

[0007] The prior art is to analyze the pore structure of natural rock samples and combine the elastic parameter test data to analyze and study the influence law of different pore structures on the elastic parameters. Since the carbonate rock pore and fracture reservoir is very complex, the deep drilling core is less, and it is very difficult to obtain natural samples with different rock structures and pore structures. Because the sample characteristics of different regions and different time strata are different, it is difficult to obtain systematic and quantitative evaluation of the influence of different dissolution and dissolution degree of pore structure on rock elastic properties, and the post-dissolution-dissolution and deep burial-hydrothermal reconstruction are the key geological conditions for the development of large carbonate rock reservoirs.

[0008] Therefore, it is one of the problems to be solved in the field to develop a quantitative evaluation method for simulating the influence of dissolution on the elastic properties of carbonate rocks using natural rocks. SUMMARY

[0009] To solve the above problems, the present application provides a quantitative evaluation method and device for the influence of dissolution on the elasticity of sedimentary rocks, a storage medium and an electronic device. The quantitative evaluation method for the influence of dissolution on the elastic properties of carbonate rocks disclosed in the present application controls the reconstruction degree of dissolution on the pore-fracture of carbonate rock reservoir by simulating different solution media, temperature, pressure and dissolution time, quantitatively analyzes the influence mechanism of dissolution pore and fracture on rock matrix and its elastic parameters, adopts different pore-fracture reservoir rock physical models based on different response mechanisms of the reservoir, and evaluates the seismic attribute response of pore-fracture carbonate rock reservoir, which helps to improve the accuracy of seismic reservoir prediction.

[0010] In a first aspect, the present application provides a quantitative evaluation method for the influence of dissolution on the elasticity of sedimentary rocks, which comprises:

[0011] determining the mineral components of a carbonate rock sample of a target layer, the content of each mineral component, the first pore characteristics and the first elastic parameters; wherein the pore characteristics include porosity and pore aspect ratio, and the elastic parameters include P-wave velocity and S-wave velocity;

[0012] determining the dissolution environment parameters according to the sedimentary environment, underground fluid condition and burial history of the area where the carbonate rock sample is located, and performing dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameters;

[0013] redetermining the second pore characteristics and the second elastic parameters of the carbonate rock sample;

[0014] establishing a quantitative evaluation model according to the mineral components, the content of each mineral component, the first pore characteristics, the first elastic parameters, the second pore characteristics and the second elastic parameters, and evaluating the influence of dissolution on the elasticity of carbonate rocks according to the quantitative evaluation model.

[0015] Further, the re-determining the second porosity characteristics and the second elastic parameters of the carbonate rock sample comprises:

[0016] According to the dissolution environment parameters, the carbonate rock sample is subjected to multiple dissolution simulation processes in sequence to obtain multiple second porosity characteristics and second elastic parameter data of different dissolution degrees.

[0017] Further, the dissolution environment parameters comprise:

[0018] dissolution temperature, dissolution pressure, dissolution reaction solution, solubility of the dissolution reaction solution, and dissolution treatment time; wherein,

[0019] In the case that the carbonate rock sample is subjected to deep formation organic acid dissolution, the dissolution reaction solution comprises acetic acid dissolution solution;

[0020] In the case that the carbonate rock sample is subjected to near-surface weathering leaching dissolution, the dissolution reaction solution comprises carbonic acid dissolution solution.

[0021] Further, the establishing a quantitative evaluation model according to the mineral components, the contents of the mineral components, the first porosity characteristics, the first elastic parameters, the second porosity characteristics, and the second elastic parameters comprises:

[0022] establishing a relationship between porosity and pressure according to the first porosity characteristics, the first elastic parameters, the second porosity characteristics, and the second elastic parameters, and determining a fitting coefficient of the relationship between porosity and pressure;

[0023] establishing a quantitative evaluation model according to the mineral components, the contents of the mineral components, and the relationship between porosity and pressure.

[0024] Further, the relationship between porosity and pressure comprises:

[0025]

[0026] wherein, is the porosity, P is the pressure, is the normal pressure porosity, and b is the fitting coefficient.

[0027] Further, the quantitative evaluation model comprises:

[0028]

[0029]

[0030] wherein, K is the bulk modulus, μ is the shear modulus, P is the formation pressure, and γ is the pore aspect ratio. Where C is porosity, C1 is a constant, and a is a constant coefficient;

[0031] The critical porosity of the rock;

[0032] f is the equivalent bulk modulus of the mineral components. i K represents the volume percentage of the i-th mineral component in the rock. i Let be the bulk modulus of the i-th mineral component in the rock;

[0033]

[0034] μ i Let be the shear modulus of the i-th mineral component in the rock.

[0035] Furthermore, the quantitative evaluation model also includes:

[0036]

[0037]

[0038] Among them, V p V represents the longitudinal wave velocity of the formation. s ρ is the shear wave velocity of the formation, ρ is the density of the carbonate rock sample, and C2 is a constant.

[0039] A second aspect of this application provides a quantitative evaluation device for the impact of dissolution on the elasticity of sedimentary rocks, the device comprising:

[0040] The first determining module is used to determine the mineral composition, the content of each mineral component, the first pore characteristics, and the first elastic parameters of the carbonate rock sample of the target layer; wherein, the pore characteristics include porosity and pore aspect ratio, and the elastic parameters include longitudinal wave velocity and transverse wave velocity;

[0041] The dissolution simulation processing module is used to determine the dissolution environment parameters based on the sedimentary environment, underground fluid conditions and burial history of the carbonate rock sample area, and to perform dissolution simulation processing on the carbonate rock sample based on the dissolution environment parameters.

[0042] The second determining module is used to redetermine the second pore characteristics and the second elastic parameter of the carbonate rock sample;

[0043] The quantitative evaluation model determining module is configured to establish a quantitative evaluation model according to the mineral components, the contents of the mineral components, the first pore characteristics, the first elastic parameters, the second pore characteristics and the second elastic parameters, and evaluate the influence of the dissolution on the elasticity of the carbonate rock according to the quantitative evaluation model.

[0044] In a third aspect, a computer readable storage medium storing a computer program is provided, the computer program being executable by one or more processors to implement the steps of the method as described above.

[0045] In a fourth aspect, an electronic device is provided, comprising a memory and one or more processors, the memory storing a computer program, the memory and the one or more processors being communicatively connected, the computer program being executable by the one or more processors to implement the steps of the method as described above.

[0046] Compared with the prior art, the technical scheme of the present application has the following advantages or beneficial effects:

[0047] The quantitative evaluation of the influence of the dissolution on the elasticity of the sedimentary rock disclosed in the present application controls the degree of reconstruction of the carbonate rock reservoir by simulating different solution media, temperatures, pressures and dissolution times, obtains the relationship between the porosity and the pressure, quantitatively analyzes the influence mechanism of the dissolution pores and fractures on the rock skeleton and the elastic parameters thereof, establishes a quantitative evaluation model of the elastic parameters of the carbonate rock and the characteristic parameters of the dissolution pores and fractures, evaluates the seismic attribute response of the porous and fractured carbonate rock reservoir, and is helpful to improve the accuracy of the seismic reservoir prediction. BRIEF DESCRIPTION OF DRAWINGS

[0048] In order to more clearly illustrate the technical scheme in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only embodiments of the present application, and those skilled in the art can obtain other drawings according to the provided drawings without any creative effort.

[0049] In addition, it should be further pointed out that, for the convenience of description, only the parts related to the present disclosure are shown in the drawings. The drawings accompanying the specification form part of the present application and are used to provide a further understanding of the present application. The schematic embodiments and their descriptions in the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0050] Figure 1 A flowchart of a quantitative evaluation method of the influence of the dissolution on the elasticity of the sedimentary rock is provided for the embodiments of the present application.

[0051] Figure 2A flow chart of another quantitative evaluation method of the influence of dissolution on the elasticity of sedimentary rock provided by the embodiments of the present application is shown in FIG. 6.

[0052] Figure 3 A CT scan diagram of the inside of a sample before and after dissolution provided by the embodiments of the present application is shown in FIG. 7.

[0053] Figure 4 A statistical diagram of the aspect ratio of pores of a sample before dissolution provided by the embodiments of the present application is shown in FIG. 8.

[0054] Figure 5 A statistical diagram of the aspect ratio of pores of a sample after dissolution provided by the embodiments of the present application is shown in FIG. 9.

[0055] Figure 6 A diagram of the variation of the P-wave velocity and the S-wave velocity of a sample with pressure before and after dissolution provided by the embodiments of the present application is shown in FIG. 10.

[0056] Figure 7 Another diagram of the variation of the P-wave velocity and the S-wave velocity of a sample with pressure before and after dissolution provided by the embodiments of the present application is shown in FIG. 11.

[0057] Figure 8 A structural diagram of a device provided by the embodiments of the present application is shown in FIG. 12.

[0058] Figure 9 A connection block diagram of an electronic device provided by the embodiments of the present application is shown in FIG. 13. DETAILED DESCRIPTION

[0059] The embodiments of the present application will be described in detail below with reference to the accompanying drawings and embodiments, so that how the present application applies technical means to solve technical problems and achieves corresponding technical effects can be fully understood and implemented. The embodiments of the present application and each feature in the embodiments can be combined with each other on the premise of no conflict, and the formed technical solutions are all within the protection scope of the present application.

[0060] It should be clear that the embodiments described below are only some of the embodiments of the present application, rather than all the embodiments. All other embodiments obtained by those skilled in the art based on the embodiments in the present application without creative labor are within the protection scope of the present application.

[0061] As known from the background, in recent years, with the increasing progress of exploration technology, carbonate rock oil and gas exploration has gradually been pushed to deep layers, for example, many deep large gas fields such as Tahe, Puguang, Yuanba, Anyue and Shanzhong have been discovered in Tarim, Sichuan and Ordos basins.

[0062] Because of deep burial, old age, long-term high temperature and high pressure, diagenesis and cementation, there are few primary pores in carbonate rocks in China. The secondary pores, fractures and caves are the important oil and gas reservoir spaces, and the dissolution is the key to the formation of secondary pores.

[0063] At present, in the research of carbonate reservoirs, scholars have carried out a lot of research work on petrology and reservoir performance evaluation. For example, some scholars believe that tectonic-pressure coupling controls fractures and dissolution, and fluid-rock interaction controls deep dissolution and pore preservation; some scholars point out that the quasi-synthetic atmospheric carbon water leaching dissolution is a necessary condition for the formation of large-scale high-quality reservoirs. For the research literature on dissolution, high temperature and high pressure dissolution kinetics simulation experiment device is used to carry out experimental research on carbonate burial dissolution temperature window and pore evolution pattern, mainly focusing on the formation and mechanism of secondary pores in the dissolution process.

[0064] After the dissolution diagenesis, the carbonate secondary pores develop, and the pore network structure is complex, which seriously affects the elastic properties of carbonate reservoirs, and the rock velocity-porosity relationship shows strong dispersion. In the existing research on the influence of dissolution on rock elastic properties, a large number of rock samples are collected to analyze the velocity change rule in the dissolution process. For example, some scholars collect multiple core samples of different ages in different regions, measure the P-wave and S-wave velocity in the laboratory, and analyze the influence of dissolution on rock velocity; some scholars select multiple dense limestone samples to measure ultrasonic velocity, and analyze the influence of pore structure on rock skeleton elastic characteristics and fluid sensitivity.

[0065] The existing technology is to analyze the pore structure of natural rock samples, and combine the elastic parameter test data to analyze and study the influence rule of different pore structures on elastic parameters. Because the carbonate pore-fracture-cave reservoir is very complex, it is difficult to drill core in deep layer, and it is difficult to obtain natural samples with different rock structures and pore structures. Because the samples of different regions and different time periods have large differences in rock skeleton characteristics, it is difficult to obtain systematic and quantitative evaluation of the influence of different dissolution and dissolution degree on rock elastic properties. Deep oil and gas exploration, especially deep carbonate oil and gas reservoirs, post-dissolution-dissolution and deep burial-hydrothermal reconstruction are the key geological conditions for the development of large carbonate reservoirs.

[0066] In view of this, the application provides a quantitative evaluation of the influence of dissolution on the elasticity of sedimentary rocks.

[0067] Embodiment one

[0068] The embodiment provides a quantitative evaluation method of the influence of dissolution on the elasticity of sedimentary rocks, Figure 1 A flowchart of the quantitative evaluation method of the influence of dissolution on the elasticity of sedimentary rocks provided by the embodiment is shown in Figure 1 The method disclosed by the embodiment includes the following steps:

[0069] Step 110, determining the mineral components of a carbonate rock sample of a target layer, the content of each mineral component, the first pore characteristics and the first elastic parameters; wherein the pore characteristics include porosity and pore aspect ratio, and the elastic parameters include P-wave velocity and S-wave velocity;

[0070] Step 120, determining the dissolution environment parameters according to the sedimentary environment, underground fluid condition and burial history of the area where the carbonate rock sample is located, and performing dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameters;

[0071] Step 130, re-determining the second pore characteristics and the second elastic parameters of the carbonate rock sample;

[0072] Step 140, establishing a quantitative evaluation model according to the mineral components, the content of each mineral component, the first pore characteristics, the first elastic parameters, the second pore characteristics and the second elastic parameters, and evaluating the influence of dissolution on the elasticity of carbonate rocks according to the quantitative evaluation model.

[0073] In some embodiments, the dissolution environment parameters include:

[0074] dissolution temperature, dissolution pressure, dissolution reaction solution, solubility of the dissolution reaction solution and dissolution treatment time; wherein,

[0075] In the case that the carbonate rock sample is subjected to deep formation organic acid dissolution, the dissolution reaction solution includes acetic acid dissolution solution;

[0076] In the case that the carbonate rock sample is subjected to near-surface weathering leaching dissolution, the dissolution reaction solution includes carbonic acid dissolution solution.

[0077] As an example, in determining the dissolution treatment condition, the dissolution environment parameters can be determined according to the deposition environment, underground fluid condition, burial history of the area where the sample is located, the dissolution environment parameters including temperature, pressure, reaction solution (medium) for dissolution and solubility, and experimental time.

[0078] If the deep formation is in contact with the source rock formation, the dissolution is mainly caused by organic acid, and acetic acid dissolution liquid is used.

[0079] If the dissolution is mainly caused by near-surface weathering and leaching, carbonic acid dissolution liquid is used.

[0080] In some embodiments, the re-determining the second porosity characteristics and the second elastic parameters of the carbonate rock sample comprises:

[0081] According to the dissolution environment parameters, the carbonate rock sample is subjected to multiple dissolution simulation treatments in sequence to obtain multiple second porosity characteristics and second elastic parameter data of different dissolution degrees.

[0082] As an example, the carbonate rock sample described above can be subjected to high-temperature and high-pressure dissolution simulation experiment by means of a dissolution simulation device, and after the experiment is completed, the sample is taken out and subjected to drying treatment. Multiple dissolution can also be performed to obtain reservoir characteristics and elastic characteristics of different dissolution degrees.

[0083] In some embodiments, the establishing a quantitative evaluation model according to the mineral components, the content of each mineral component, the first porosity characteristics, the first elastic parameters, the second porosity characteristics and the second elastic parameters comprises:

[0084] According to the first porosity characteristics, the first elastic parameters, the second porosity characteristics and the second elastic parameters, a relationship between porosity and pressure is established, and a fitting coefficient of the relationship between porosity and pressure is determined.

[0085] According to the mineral components, the content of each mineral component and the relationship between porosity and pressure, a quantitative evaluation model is established.

[0086] In some embodiments, the relationship between porosity and pressure comprises:

[0087]

[0088] wherein, porosity; P is pressure, unit MPa; is the normal pressure porosity; b is the fitting coefficient, which is obtained from experimental data.

[0089] In some embodiments, the quantitative evaluation model comprises:

[0090]

[0091]

[0092] wherein K is the bulk modulus, unit GPa; μ is the shear modulus, unit GPa; P is the formation pressure, unit MPa; γ is the pore aspect ratio; is the porosity; C1 is a constant; a is a coefficient related to the fluid, wherein η is the fluid viscosity, unit Pa.s;

[0093] is the critical porosity of the rock;

[0094] is the equivalent bulk modulus of the mineral components (unit GPa); f i is the volume proportion of the i-th mineral component in the rock; K i is the bulk modulus of the i-th mineral component in the rock (unit GPa);

[0095]

[0096] μ i is the shear modulus of the i-th mineral component in the rock (unit GPa).

[0097] In some embodiments, the quantitative evaluation model further comprises:

[0098]

[0099]

[0100] wherein V p is the formation P-wave velocity, unit m / s; V s is the formation S-wave velocity, unit m / s; ρ is the density of the carbonate rock sample, unit g / cm 3 ; C2 is a constant.

[0101] The quantitative evaluation method for the influence of dissolution on the elastic properties of carbonate rocks disclosed in the present application controls the reconstruction degree of pores and fractures of carbonate rock reservoirs by simulating different solution media, temperatures, pressures and dissolution times, quantitatively analyzes the influence mechanism of dissolved pores and fractures on the rock framework and its elastic parameters, adopts different pore and fracture reservoir rock physical models based on different response mechanisms of the reservoirs, evaluates the seismic attribute response of the pore and fracture carbonate rock reservoirs, and is helpful to improve the seismic reservoir prediction accuracy.

[0102] Embodiment Two

[0103] The embodiment further explains and illustrates the quantitative evaluation method of the influence of the dissolution on the elasticity of the sedimentary rock disclosed in the first embodiment, as shown in the following. Figure 1 The method disclosed in the embodiment includes the following steps:

[0104] In step 110, the mineral components of the carbonate rock sample of the target layer, the content of each mineral component, the first pore characteristics, and the first elastic parameters are determined. The pore characteristics include porosity and pore aspect ratio, and the elastic parameters include P-wave velocity and S-wave velocity.

[0105] It should be noted that the carbonate rock in the embodiment includes limestone and dolomite.

[0106] In some embodiments, before the determination of the mineral components of the carbonate rock sample of the target layer, the content of each mineral component, the first pore characteristics, and the first elastic parameters, the method further includes:

[0107] The carbonate rock sample of the target layer is selected and prepared.

[0108] As an example, first, the carbonate rock sample of the target layer is selected according to the experimental purpose; then, the sample is prepared into a standard sample according to the physical property and elastic analysis requirements. The porosity and elastic test requirements standard plunger sample: diameter 25 mm, length about 2 times the diameter, the plunger sample is polished at both ends, the end surface is perpendicular to the cylindrical surface and has no damage, and the sample is subjected to drying treatment. Thin section analysis, the sample is prepared into a thin section.

[0109] As an example, in the determination of the mineral components of the carbonate rock sample of the target layer, the content of each mineral component, the first pore characteristics, and the first elastic parameters, the following methods can be adopted:

[0110] (1) The mineral components and content of the sample are determined by XRD or other methods; the phase identification and quantitative calculation of the content of each mineral component are performed by X-ray powder diffraction;

[0111] (2) The pore aspect ratio characteristics inside the sample before and after dissolution are obtained by CT scanning or scanning electron microscopy;

[0112] (3) The porosity of the sample before dissolution is obtained by the overpressure porosity tester;

[0113] (4) The elastic parameters such as P-wave velocity and S-wave velocity under the formation temperature and pressure conditions are obtained by the ultrasonic testing instrument.

[0114] It should be noted that the thin section identification analysis, XRD diffraction mineral analysis, CT scanning, overpressure porosity test, and ultrasonic elastic parameter test are all conventional test and analysis methods in the field, for example, the thin section identification analysis refers to SY / T 5368-2016 Rock Thin Section Identification, the XRD diffraction mineral analysis refers to SYT 5163-2018 Clay Mineral and Common Non-Clay Mineral X-Ray Diffraction Analysis Method in Sedimentary Rock, the overpressure porosity and permeability test refers to SY / T 6385-2016 Method for Determining Rock Porosity and Permeability under Overpressure, and the ultrasonic elastic parameter test refers to SY / T 6351-2012 Laboratory Measurement Specification for Rock Sample Acoustic Wave Characteristics.

[0115] In step 120, a dissolution environment parameter is determined according to the sedimentary environment, underground fluid condition, and burial history of the area where the carbonate rock sample is located, and the carbonate rock sample is subjected to dissolution simulation processing according to the dissolution environment parameter.

[0116] In some embodiments, the dissolution environment parameter includes:

[0117] dissolution temperature, dissolution pressure, dissolution reaction solution, dissolution reaction solution solubility, and dissolution processing time; wherein,

[0118] In the case of deep formation organic acid dissolution of the carbonate rock sample, the dissolution reaction solution includes acetic acid dissolution liquid;

[0119] In the case of near-surface weathering leaching dissolution of the carbonate rock sample, the dissolution reaction solution includes carbonic acid dissolution liquid.

[0120] As an example, in determining the dissolution processing conditions, the dissolution environment parameter including temperature, pressure, dissolution reaction solution (medium), and solubility can be determined according to the sedimentary environment, underground fluid condition, and burial history of the area where the sample is located.

[0121] If the deep formation is in contact with the source rock layer and is mainly subjected to organic acid dissolution, acetic acid dissolution liquid is used;

[0122] If it is mainly subjected to near-surface weathering leaching dissolution, carbonic acid dissolution liquid is used.

[0123] As an example, the carbonate rock sample can be subjected to high-temperature and high-pressure dissolution simulation experiment by means of a dissolution simulation device, and after the experiment is completed, the sample is taken out and subjected to drying treatment; and the sample after dissolution is subjected to CT scanning to obtain the pore characteristics of the sample, test the porosity, and test the longitudinal wave velocity, transverse wave velocity, and other elastic parameters under the formation temperature and pressure conditions.

[0124] Step 130, re-determining the second porosity characteristics and the second elastic parameters of the carbonate rock sample.

[0125] As an example, the CT scanning of the sample after dissolution is performed to obtain the porosity characteristics of the sample, and the porosity is tested, and the elastic parameters such as the longitudinal wave velocity and the transverse wave velocity under the formation temperature and pressure conditions are tested.

[0126] In some embodiments, the re-determining the second porosity characteristics and the second elastic parameters of the carbonate rock sample comprises:

[0127] According to the dissolution environment parameters, the carbonate rock sample is subjected to multiple dissolution simulation processes in sequence to obtain multiple sets of second porosity characteristics and second elastic parameter data of different dissolution degrees.

[0128] Further, the step 120 can be repeatedly performed to perform multiple dissolution processes on the above-mentioned carbonate rock sample, thereby obtaining reservoir characteristics and elastic characteristics of different dissolution degrees.

[0129] It should be noted that the formation temperature and pressure conditions refer to the formation temperature and the overburden formation static pressure. The simulation experiment process can determine the formation temperature and the overburden formation static pressure according to the formation depth, the geothermal gradient of the research area, and the formation static pressure gradient.

[0130] Step 140, establishing a quantitative evaluation model according to the mineral components, the contents of the mineral components, the first porosity characteristics, the first elastic parameters, the second porosity characteristics, and the second elastic parameters, and evaluating the influence of dissolution on the elasticity of the carbonate rock according to the quantitative evaluation model.

[0131] In some embodiments, the establishing a quantitative evaluation model according to the mineral components, the contents of the mineral components, the first porosity characteristics, the first elastic parameters, the second porosity characteristics, and the second elastic parameters comprises:

[0132] According to the first porosity characteristics, the first elastic parameters, the second porosity characteristics, and the second elastic parameters, a relationship between the porosity and the pressure is established, and a fitting coefficient of the relationship between the porosity and the pressure is determined.

[0133] A quantitative evaluation model is established according to the mineral components, the contents of the mineral components, and the relationship between the porosity and the pressure.

[0134] As an example, the changes of the reservoir characteristics and the elastic characteristics before and after the dissolution are quantitatively evaluated: the changes of the mineral components, the porosity, the porosity characteristics, the longitudinal wave velocity, the transverse wave velocity, the ratio of the longitudinal wave velocity to the transverse wave velocity, and other parameters in steps 110 and 120 are compared, the relationship between the porosity and the pressure is established, and the fitting coefficient is determined.

[0135] In some embodiments, the relationship between the porosity and the pressure comprises:

[0136]

[0137] wherein, is the porosity; P is the pressure, in MPa; is the normal porosity; b is a fitting coefficient, obtained from experimental data.

[0138] In some embodiments, the quantitative evaluation model comprises:

[0139]

[0140]

[0141] wherein, K is the bulk modulus, in GPa; μ is the shear modulus, in GPa; P is the formation pressure, in MPa; γ is the pore aspect ratio; is the porosity; C1 is a constant; a is a fluid-related coefficient, wherein, η is the fluid viscosity, in Pa.s;

[0142] is the critical porosity of the rock;

[0143] is the equivalent bulk modulus of the mineral components, in GPa; f i is the volume fraction of the i-th mineral component in the rock; K i is the bulk modulus of the i-th mineral component in the rock, in GPa;

[0144]

[0145] μ i is the shear modulus of the i-th mineral component in the rock, in GPa.

[0146] Preferably, the value of C1 is 4 / 15.

[0147] As understood by those skilled in the art, the fluid viscosity η is a physical parameter of the fluid contained in the rock, and is a known quantity.

[0148] Further, the quantitative evaluation model can also be converted into the formation longitudinal wave velocity (V p ) and the transverse wave velocity (V s ) as follows:

[0149]

[0150]

[0151] wherein, V p is the longitudinal wave velocity of the formation, in m / s; V s is the transverse wave velocity of the formation, in m / s; p is the density of the carbonate sample, in g / cm 3 ; C2 is a constant.

[0152] Preferably, C2 has a value of 4 / 3.

[0153] This embodiment controls the degree of reconstruction of the pore-fracture of the carbonate reservoir by simulating different solution media, temperatures, pressures and dissolution times, quantitatively analyzes the influence mechanism of the dissolution pore and fracture on the rock skeleton and its elastic parameters, adopts different dissolution pore-fracture reservoir rock physics models based on different response mechanisms of the reservoir, evaluates the seismic attribute response of the pore-fracture carbonate reservoir, and is helpful to improve the accuracy of seismic reservoir prediction.

[0154] Embodiment Three

[0155] This embodiment is a specific example of the present application.

[0156] At present, large oil and gas reservoirs have been found in deep carbonate rocks of the Ordovician system in the Tarim Basin, the Changxing Formation and the Dengying Formation in the Sichuan Basin, etc. Although the primary pores are few due to the long-term effects of high temperature and high pressure, diagenesis and cementation, the developed secondary pores, fractures and caves become important oil and gas reservoir spaces, and dissolution is the key to the formation of the secondary pores.

[0157] In this embodiment, lithological samples in the Sichuan Basin are selected to carry out dissolution simulation experiments to compare and analyze the changes in the reservoir performance and elastic property characteristics of different types of rocks before and after dissolution.

[0158] Specifically, as shown in Figure 2 the method disclosed in this embodiment includes the following steps:

[0159] Step 1. First, according to the core observation description and combined with the reservoir geological research, different types of carbonate rock samples in the target interval are selected.

[0160] The different types of carbonate rock samples include one or more of oolitic dolarenite and micritic dolomite with oolitic pore development, dolarenite and powdery dolomite with dissolution pore and cave development, fractured sparry grainstone and micritic-powdery dolomite, etc., which are mainly determined according to thin section identification. In this embodiment, two groups of dolomite samples are selected.

[0161] Then, according to the physical properties and elastic analysis requirements, the sample is prepared into a standard sample. Among them, the porosity and elastic test requires a standard plunger sample: diameter 25 mm, length about 2 times of the diameter, generally length 40-50 mm, the plunger sample is polished at both ends, the end surface is perpendicular to the cylindrical surface, without damage, and is treated by drying. For thin section analysis, the sample is prepared into a thin section.

[0162] Step 2, storage and elastic property test of the sample before dissolution, in particular:

[0163] (1) The mineral composition and content of the sample are determined by XRD or other methods; the phase identification and quantitative calculation of the content of each mineral component are carried out by X-ray powder diffraction. The X-ray diffraction analysis data of the carbonate rock sample are shown in Table 1. From the structure, sample 1 and sample 2 are mainly composed of dolomite, accounting for 85.8% and 88.4% respectively, containing 11.1% and 4.0% of quartz, containing 1.6% and 6.0% of clay, and containing 1.5% of plagioclase and other minerals. Combined with thin section identification analysis, it is determined that the lithology is siliceous powder crystal dolomite and micritic dolomite.

[0164] (2) The pore characteristics of the sample before dissolution are obtained by CT scanning or scanning electron microscopy. Among them, sample 1 has low porosity, mainly intergranular pores, and poor connectivity of pore cracks; sample 2 has developed dissolution pores, which are dissolution holes, and has large pores, medium porosity and good pore connectivity. As shown in part a of Figure 3 , it is a CT scanning image of the inside of a dissolution experiment sample. The aspect ratio distribution characteristics of the sample pores are shown in Figure 4 , and the main pore aspect ratio is 0.01.

[0165] (3) The porosity and permeability of the sample before dissolution are obtained by using the overpressure porosity and permeability tester (see Table 1).

[0166] (4) The ultrasonic testing instrument is used to obtain the elastic parameters such as longitudinal wave velocity (V p ), transverse wave velocity (V s ), longitudinal to transverse wave velocity ratio (V p / V s ) under the conditions of formation temperature and pressure (see Tables 3, Figure 6 and Figure 7 , wherein, Figure 6 and Figure 7 are the Vp and Vs data of the first elastic parameters before dissolution).

[0167] Step 3, determine the dissolution experiment conditions, in particular:

[0168] According to the sedimentary environment, underground fluid situation and burial history of the sample, the dissolution experiment environmental parameters are determined, including temperature, pressure, reaction solution (medium) and solubility, and experiment time.

[0169] In this study, the deep formation is in contact with the source rock layer, mainly affected by organic acid dissolution, and 1ml / L acetic acid solution (pH = 4.002) is used. The temperature used in this experiment is 60℃, the pressure is in the range of 10-52Mpa, and the experimental time is 6h.

[0170] Step 4, using the dissolution simulation device, the carbonate rock sample is subjected to high temperature and high pressure dissolution simulation experiment. After the experiment, the sample is taken out and the sample after dissolution is placed for half a day to dry naturally, and then dried at 110℃ in a constant temperature environment.

[0171] Step 5, the sample after dissolution is subjected to CT scanning to obtain the pore characteristics of the sample. From the CT image of the sample after dissolution, it can be seen that the intergranular pore of sample 1 is generally dissolved and enlarged, and some cracks are generated, and the aspect ratio of the sample as a whole is slightly smaller. The pore of sample 2 is dissolved and enlarged, and the local particle is dissolved into a solution hole (hole), which evolves into a hole type (b part in Figure 3 Figure 5 The aspect ratio of the main pore is 0.3. The aspect ratio of the sample as a whole increases from 0.01 to about 0.3 before and after dissolution (reference Figure 4 、 Figure 5 ). The porosity and permeability are tested under the same temperature and pressure conditions as step 2 (see table 2). The longitudinal wave velocity, transverse wave velocity and other elastic parameters of the sample are tested under the same formation temperature and pressure conditions as step 2 (reference Figure 6 、 Figure 7 , wherein, Figure 6 and Figure 7 the hollow square in the figure is the second elastic parameter data after dissolution).

[0172] Step 6, quantitative evaluation of the change of reservoir characteristics and elastic characteristics before and after dissolution experiment: compare the change of mineral composition, porosity, permeability and other parameters of step 2 and step 5 (table 2). The dissolution of mud crystal dolomite reservoir porosity is increased by 1.76%, which is about 117% of the original, and the permeability is increased to 220%. The porosity of the dense pore type siliceous powder crystal dolomite is increased by 0.71%, which is 120% of the original. The permeability is greatly improved from 0.026 to 0.26, which is 10 times of the original. The elastic parameters of the sample before and after dissolution experiment are shown in table 3.

[0173] Table 1 mineral content, porosity and permeability test data before and after dissolution experiment

[0174]

[0175] Table 2 change of porosity and permeability of dissolution experiment ​

[0176]

[0177] Table 3 Elastic parameters of samples before and after dissolution experiment

[0178]

[0179] Step 7, quantitatively analyze the influence of dissolution pores and cracks on the elastic parameters of the rock, and establish a quantitative evaluation model of the elastic parameters of carbonate rocks and the characteristic parameters of dissolution pores and cracks.

[0180] According to the analysis of the relationship between the longitudinal and transverse wave velocities and the pressure before and after dissolution, for sample 2: after the dissolution experiment, the velocity gradually decreases with the increase of the pressure, the pores and cracks increase, and the velocity ratio approaches the value of 2, which reflects the characteristics of the dolomite mineral of the rock skeleton; for sample 1: after dissolution, the velocity decreases, the pores and cracks increase, especially the small aspect ratio cracks increase, and the saturated fluid velocity changes obviously with the pressure. In describing the elastic properties of the reservoir rock, the relationship between the elastic parameters of carbonate rocks and the characteristic parameters of dissolution pores and cracks can be evaluated by using the characteristics of the formation velocity changing with the pressure of different pore structures, and the propagation law of the real deep underground seismic wave can be analyzed, which can improve the accuracy of reservoir prediction and fluid identification in oil and gas seismic exploration and development.

[0181] In this embodiment, by simulating different solution media, temperature, pressure and dissolution time, the degree of reconstruction of carbonate reservoirs by dissolution can be controlled, and the influence law and mechanism analysis of the reconstruction of the surface water, deep hydrothermal dissolution and other dissolution effects on the shallow and deep carbonate reservoir skeleton, pore type, pore structure and reservoir elastic properties can be studied. By using the porosity of the test sample before and after dissolution, the relationship between porosity and pressure is established, the influence of dissolution pores and cracks on the rock skeleton and its elastic parameters is quantitatively analyzed, the rock physics model of the dissolution pore and crack reservoir is more truly described, the seismic attribute response of the pore and crack carbonate reservoir is evaluated, which helps to improve the accuracy of seismic reservoir prediction.

[0182] Embodiment four

[0183] The embodiment provides a quantitative evaluation device for the influence of dissolution on the elasticity of sedimentary rocks. The device embodiment can be used to execute the method embodiments of the application. For details not disclosed in the device embodiment, please refer to the method embodiments of the application. Figure 8 The structure diagram of a device provided by the embodiment of the application is shown in Figure 8 The device 800 disclosed in the embodiment includes:

[0184] The first determining module 801 is configured to determine mineral components of a carbonate rock sample of a target layer, contents of each mineral component, first pore characteristics, and first elastic parameters; wherein the pore characteristics include porosity and pore aspect ratio, and the elastic parameters include longitudinal wave velocity and transverse wave velocity.

[0185] The dissolution simulation processing module 802 is configured to determine dissolution environment parameters according to a deposition environment, underground fluid condition, and burial history of a region where the carbonate rock sample is located, and perform dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameters.

[0186] The second determining module 803 is configured to redetermine second pore characteristics and second elastic parameters of the carbonate rock sample.

[0187] The quantitative evaluation model determining module 804 is configured to establish a quantitative evaluation model according to the mineral components, the contents of each mineral component, the first pore characteristics, the first elastic parameters, the second pore characteristics, and the second elastic parameters, and evaluate the influence of dissolution on the elasticity of carbonate rock according to the quantitative evaluation model.

[0188] In some embodiments, the second determining module 803 is configured to perform multiple dissolution simulation processes on the carbonate rock sample according to the dissolution environment parameters in sequence to obtain multiple pieces of second pore characteristic and second elastic parameter data of different dissolution degrees.

[0189] In some embodiments, the dissolution environment parameters include:

[0190] dissolution temperature, dissolution pressure, dissolution reaction solution, solubility of the dissolution reaction solution, and dissolution processing time; wherein,

[0191] In the case that the carbonate rock sample is subjected to deep formation organic acid dissolution, the dissolution reaction solution includes acetic acid dissolution liquid.

[0192] In the case that the carbonate rock sample is subjected to near-surface weathering leaching dissolution, the dissolution reaction solution includes carbonic acid dissolution liquid.

[0193] In some embodiments, the quantitative evaluation model determining module 804 includes a determining unit and an establishing unit; wherein,

[0194] The determining unit is configured to establish a relationship between porosity and pressure according to the first pore characteristics, the first elastic parameters, the second pore characteristics, and the second elastic parameters, and determine a fitting coefficient of the relationship between the porosity and the pressure.

[0195] The establishing unit is configured to establish a quantitative evaluation model according to the mineral components, the contents of each mineral component, and the relationship between the porosity and the pressure.

[0196] In some embodiments, the relationship between the porosity and the pressure comprises:

[0197]

[0198] wherein, is the porosity; P is the pressure, in MPa; is the normal porosity; b is a fitting coefficient, obtained from experimental data.

[0199] In some embodiments, the quantitative evaluation model comprises:

[0200]

[0201]

[0202] wherein, K is the bulk modulus, in GPa; μ is the shear modulus, in GPa; P is the formation pressure, in MPa; γ is the pore aspect ratio; is the porosity; C1 is a constant; a is a fluid-related coefficient, wherein, η is the fluid viscosity, in Pa.s;

[0203] is the critical porosity of the rock;

[0204] is the equivalent bulk modulus of the mineral components, in GPa; f i is the volume fraction of the i-th mineral component in the rock; K i is the bulk modulus of the i-th mineral component in the rock, in GPa;

[0205]

[0206] μ i is the shear modulus of the i-th mineral component in the rock, in GPa.

[0207] In some embodiments, the quantitative evaluation model further comprises:

[0208]

[0209]

[0210] wherein, V p is the formation P-wave velocity, in m / s; V s is the formation S-wave velocity, in m / s; ρ is the density of the carbonate rock sample, in g / cm 3 ; C2 is a constant.

[0211] Those skilled in the art can understand that the structure shown in the above embodiments is not a limitation to the device of the present application, and can include more or less modules / units than the drawings, or combine certain modules / units, or arrange different modules / units. Figure 8 The structure shown in the above embodiments is not a limitation to the device of the present application, and can include more or less modules / units than the drawings, or combine certain modules / units, or arrange different modules / units.

[0212] Those skilled in the art should understand that the modules or steps of the present application described above can be realized by a general computing device, which can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices. Alternatively, they can be realized by program codes executable by a computing device, so that they can be stored in a storage device and executed by a computing device, and in some cases, the steps shown or described can be executed in an order different from that shown here, or they can be manufactured into individual integrated circuit modules, or multiple modules or steps among them can be manufactured into a single integrated circuit module.

[0213] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of each module of the device for quantitatively evaluating the effect of dissolution on the elasticity of sedimentary rocks can refer to the corresponding process in the foregoing method embodiments, and the present embodiment will not be repeated here.

[0214] Embodiment Five

[0215] The present embodiment provides a computer readable storage medium. The computer readable storage medium stores a computer program, and the computer program can implement the method in the foregoing method embodiments when executed by a processor.

[0216] The computer readable storage medium can also individually include a computer program, a data file, a data structure, etc., or a combination thereof. The computer readable storage medium or the computer program can be specifically designed and understood by those skilled in the computer software field, or can be known and available to those skilled in the computer software field. Examples of the computer readable storage medium include: magnetic media, such as a hard disk, a floppy disk, and a magnetic tape; optical media, such as a CD ROM disk and a DVD; a magneto-optical medium, such as an optical disk; and a hardware device specifically configured to store and execute a computer program, such as a read-only memory (ROM), a random access memory (RAM), a flash memory; or a server, an app application store, etc. Examples of the computer program include machine code (e.g., code generated by a compiler) and a file containing high-level code that can be executed by a computer by using an interpreter. The described hardware device can be configured to function as one or more software modules to perform the above-described operations and methods, and vice versa. In addition, the computer readable storage medium can be distributed in a networked computer system, and the program code or computer program can be stored and executed in a distributed manner.

[0217] Embodiment six

[0218] The embodiment provides a computer program product. The computer program product includes a computer program or instructions, which, when executed by a processor, implement all or part of the steps of the method in the foregoing method embodiments, which will not be repeated here.

[0219] Further, the computer program product can include one or more computer executable components configured to perform embodiments when the program is run; the computer program product can also include a computer program tangibly embodied on a computer readable medium, the computer program containing program code for performing any of the methods in the embodiments of the present disclosure. In such embodiments, the computer program can be downloaded and installed from a network by a communication part, and / or installed from a removable medium.

[0220] Embodiment seven

[0221] The embodiment provides an electronic device. Figure 9 A connection block diagram of an electronic device provided by the embodiment of the present application is shown in Figure 9 The electronic device 900 can include one or more processors 901, a memory 902, a multimedia component 903, an input / output (I / O) interface 904, and a communication component 905, as shown in the figure.

[0222] The one or more processors 901 are configured to perform, wholly or partly, the steps of the methods described above. The memory 902 is configured to store various types of data, which can include, for example, instructions of any application programs or methods in the electronic device, and application-related data.

[0223] The memory 902 can be implemented by any type of volatile or nonvolatile memory or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disc, or optical disc.

[0224] The one or more processors 901 can be implemented by an application specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic elements, to perform, wholly or partly, the steps of the methods described above.

[0225] The multimedia component 903 can include a screen, which can be a touch screen, and an audio component for outputting and / or inputting audio signals. For example, the audio component can include a microphone for receiving external audio signals. The received audio signals can be further stored in the memory or transmitted through the communication component. The audio component also includes at least one speaker for outputting audio signals.

[0226] The I / O interface 904 provides an interface between the one or more processors 901 and other interface modules, which can be a keyboard, a mouse, a button, etc. These buttons can be virtual buttons or physical buttons.

[0227] The communication component 905 is configured to perform wired or wireless communication between the electronic device 900 and other devices. The wired communication includes communication through a network port, a serial port, etc., and the wireless communication includes Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, 4G, 5G, or a combination of one or more of them.

[0228] In summary, the present application provides a quantitative evaluation method and device for the influence of dissolution on the elasticity of sedimentary rock, a computer readable storage medium and an electronic device. The quantitative evaluation method for the influence of dissolution on the elastic properties of carbonate rock disclosed in the present application controls the reconstruction degree of pores and fractures of carbonate rock reservoirs by simulating different solution media, temperatures, pressures and dissolution times, quantitatively analyzes the influence mechanism of dissolution pores and fractures on the rock skeleton and its elastic parameters, adopts different pore and fracture reservoir rock physical models based on different response mechanisms of the reservoir, evaluates the seismic attribute response of the pore and fracture carbonate rock reservoir, and helps to improve the seismic reservoir prediction accuracy.

[0229] It should also be understood that the methods or apparatus disclosed in the embodiments of the present application can also be implemented in other ways. The above-described method or apparatus embodiments are only illustrative, for example, the flowchart and block diagram in the drawings show the possible implementation architecture, function and operation of the method and apparatus according to the embodiments of the present application. In this regard, each block in the flowchart or block diagram can represent a module, a computer program segment or a part of a computer program, which includes one or more computer programs for implementing the specified logic function. It should also be noted that in some alternative implementations, the functions noted in the block can occur in different orders from that noted in the drawings, and in fact, they can be executed substantially in parallel, and sometimes in reverse order, depending on the function involved. It should also be noted that each block in the block diagram and / or flowchart, and the combination of blocks in the block diagram and / or flowchart, can be implemented by a dedicated hardware-based system for performing the specified function or action, or can be implemented by a combination of dedicated hardware and computer programs.

[0230] In this application, the terms "comprise", "contain", or any other variant thereof are intended to cover a non-exclusive inclusion, so that processes, methods, articles, or devices that comprise a list of elements do not only include those elements, but also include other elements that are not expressly listed, or further include elements inherent in such processes, methods, articles, or devices. Without more limitations, the element defined by the statement "comprises a" does not exclude the presence of additional identical elements in the process, method, device or equipment comprising the element; if there is a description of "first", "second", etc., it is only for the purpose of description, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of indicated technical features or the sequence of indicated technical features; in the description of the present application, unless otherwise specified, the meaning of the term "a plurality of" or "a plurality" is at least two; if there is a description of a server, it should be noted that the server can be a stand-alone physical server or terminal, or a server cluster composed of multiple physical servers, or a cloud server capable of providing cloud server, cloud database, cloud storage and CDN and other basic cloud computing services; if there is a description of a smart terminal or a mobile device in the present application, it should be noted that the smart terminal or mobile device can be a mobile phone, a tablet computer, a smart watch, a netbook, a wearable electronic device, a personal digital assistant (PDA), an augmented reality technology device (AR), a virtual reality device (VR), a smart television, a smart sound, a personal computer (PC), etc., but is not limited thereto, and the specific form of the smart terminal or mobile device is not specially limited in the present application.

[0231] Finally, it should be noted that in the description of the present application, the description of the terms "one embodiment", "some embodiments", "example", "one example" or "some examples" means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present application, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0232] Although the embodiments of the present application have been shown and described above, it is understood that all the above-described embodiments are exemplary only, the contents described are merely adopted for the purpose of facilitating the understanding of the present application, and are not intended to limit the present application. Any person skilled in the art to which the present application belongs can make any modification and change in the form and details without departing from the spirit and scope of the present application, but the protection scope of the present application shall be subject to the scope defined by the appended claims.

Claims

1. A method for quantitatively evaluating the effect of dissolution on the elasticity of sedimentary rocks, characterized by, The method comprises: determining mineral components, content of each mineral component, first pore characteristics and first elastic parameters of a carbonate rock sample of a target layer; wherein the pore characteristics comprise porosity and pore aspect ratio, and the elastic parameters comprise longitudinal wave velocity and transverse wave velocity; determining a dissolution environment parameter according to a sedimentary environment, underground fluid condition and burial history of a region where the carbonate rock sample is located, and performing dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameter; redetermining second pore characteristics and second elastic parameters of the carbonate rock sample; establishing a quantitative evaluation model according to the mineral components, content of each mineral component, first pore characteristics, first elastic parameters, second pore characteristics and second elastic parameters, and evaluating the influence of dissolution on the elasticity of carbonate rock according to the quantitative evaluation model; wherein the quantitative evaluation model comprises: wherein, is the bulk modulus, is the shear modulus, is the pressure, is the pore aspect ratio, is the porosity, is a constant, is a constant coefficient; , is the critical porosity of the rock; K i is the bulk modulus of the i-th mineral component, Vi is the volume fraction of the i-th mineral component in the rock, Ki is the bulk modulus of the i-th mineral component in the rock. ； , G, is the shear modulus of the ith mineral component in the rock; The quantitative evaluation model further comprises: wherein, is the formation P-wave velocity, is the formation S-wave velocity, is the density of the carbonate sample, is a constant.

2. The method for quantitative evaluation of the effect of dissolution on the elasticity of sedimentary rocks according to claim 1, characterized in that, The redetermination of the second pore characteristics and the second elastic parameters of the carbonate rock sample comprises: performing multiple dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameter in sequence to obtain multiple second pore characteristics and second elastic parameter data of different dissolution degrees.

3. The method for quantitatively evaluating the elastic effect of dissolution on sedimentary rock according to claim 1, characterized in that, The dissolution environment parameter comprises: dissolution temperature, dissolution pressure, dissolution reaction solution, solubility of the dissolution reaction solution and dissolution processing time; wherein, in the case that the carbonate rock sample is subjected to deep formation organic acid dissolution, the dissolution reaction solution comprises acetic acid dissolution liquid; in the case that the carbonate rock sample is subjected to near-surface weathering leaching dissolution, the dissolution reaction solution comprises carbonic acid dissolution liquid.

4. The method for quantitatively evaluating the elastic effect of dissolution on sedimentary rock according to claim 1, characterized in that, The establishment of the quantitative evaluation model according to the mineral components, content of each mineral component, first pore characteristics, first elastic parameters, second pore characteristics and second elastic parameters comprises: establishing a relationship between porosity and pressure according to the first pore characteristics, first elastic parameters, second pore characteristics and second elastic parameters, and determining a fitting coefficient of the relationship between porosity and pressure; establishing a quantitative evaluation model according to the mineral components, content of each mineral component and the relationship between porosity and pressure.

5. The method for quantitatively evaluating the elastic effect of dissolution on sedimentary rocks according to claim 4, characterized in that, The relationship between porosity and pressure comprises: wherein, is the porosity, is the pressure, is the normal pressure porosity, is the fitting coefficient.

6. A device for quantitatively evaluating the effect of dissolution on the elasticity of sedimentary rocks, characterized by, comprises: a first determination module configured to determine mineral components, content of each mineral component, first pore characteristics and first elastic parameters of a carbonate rock sample of a target layer; wherein the pore characteristics comprise porosity and pore aspect ratio, and the elastic parameters comprise longitudinal wave velocity and transverse wave velocity; a dissolution simulation processing module configured to determine a dissolution environment parameter according to a sedimentary environment, underground fluid condition and burial history of a region where the carbonate rock sample is located, and perform dissolution simulation processing on the carbonate rock sample according to the dissolution environment parameter; a second determination module configured to redetermine second pore characteristics and second elastic parameters of the carbonate rock sample; The quantitative evaluation model determining module is configured to establish a quantitative evaluation model according to the mineral components, the content of each mineral component, the first pore characteristics, the first elastic parameters, the second pore characteristics and the second elastic parameters, and evaluate the influence of the dissolution on the elasticity of the carbonate rock according to the quantitative evaluation model. The quantitative evaluation model includes: wherein, is the bulk modulus, is the shear modulus, is the pressure, is the pore aspect ratio, is the porosity, is a constant, is a constant coefficient; , is the critical porosity of the rock; is the equivalent bulk modulus of the mineral components, is the volume fraction of the i-th mineral component in the rock, is the bulk modulus of the i-th mineral component in the rock; ； , G, is the shear modulus of the ith mineral component in the rock; The quantitative evaluation model further includes: wherein, is the formation P-wave velocity, is the formation S-wave velocity, is the density of the carbonate sample, is a constant.

7. A computer-readable storage medium, characterized in that, The computer program stored in the computer readable storage medium, when executed by one or more processors, implements the quantitative evaluation method of the influence of the dissolution on the elasticity of the sedimentary rock according to any one of claims 1 to 5.

8. An electronic device, comprising: The computer program stored in the computer readable storage medium, when executed by one or more processors, implements the quantitative evaluation method of the influence of the dissolution on the elasticity of the sedimentary rock according to any one of claims 1 to 5.

Images