Method, device and equipment for analyzing elasticity modulus of rock and medium

Abstract

The invention discloses a method, a device, equipment and a medium for analyzing the elasticity modulus of rock. The method comprises the following steps: acquiring a rock sample to be analyzed, acting fillers of at least two pore types on the rock sample aiming at a first experiment, and determining an experiment result corresponding to the rock sample based on an initial elastic modulus and the fillers; aiming at the experiment after the first experiment, taking the rock sample after the action of all the fillers in the previous experiment as the rock sample of the current experiment, taking the second elastic modulus in the previous experiment as the initial elastic modulus of the current experiment, and repeatedly executing the operation of the previous experiment; and determining the influence of different fillers on the elastic modulus of the rock sample based on experimental results of multiple experiments. The problem that the elastic parameter analysis effect is poor due to the fact that a single pore type is added every time is solved, the calculation precision of the elastic modulus of the rock is improved, and the effect of improving the accuracy of elastic parameter analysis is achieved.

Description

Technical Field

[0001] The invention relates to the field of rock oil reservoir exploration, and in particular to a method, device, equipment and medium for analyzing rock elastic modulus. Background Art

[0002] With the development of oil and gas exploration technology, the study of rocks with complex pore structures is usually involved in the process of studying reservoirs rich in oil and gas resources and predicting fluid reservoirs. For example, when seismic waves propagate in rocks, different porosities and pore shapes, as well as different pore fluids, will affect the elastic parameters of the rock, such as the modulus. This pore fluid with different porosities and different pore structure types can be used to perform rock physics modeling and analyze rock elastic parameters to invert the physical properties, elastic parameters and fluid characteristics of underground reservoirs, and then provide reliable technical support for reservoir oil and gas exploration and development based on the inversion analysis results.

[0003] The existing rock elastic parameter analysis method usually adds different types of pore fluids to the same rock matrix, and analyzes the influence of different pore fluids on the elastic parameters of the rock matrix based on the property information of the rock matrix after addition. This method can only add a single pore type at a time, and the calculation process is cumbersome. Introducing multiple pore structures requires repeated iterations, resulting in poor elastic parameter analysis results and high costs. Summary of the invention

[0004] The present invention provides a method, device, equipment and medium for analyzing the elastic modulus of rock, so as to improve the calculation accuracy of the elastic modulus of rock, and at the same time achieve the technical effect of improving the accuracy of elastic parameter analysis and reducing analysis costs.

[0005] According to one aspect of the present invention, a method for analyzing rock elastic modulus is provided, the method comprising:

[0006] Obtain a rock sample to be analyzed, and for a first experiment, apply fillers of at least two pore types to the rock sample, and determine an experimental result corresponding to the rock sample based on an initial elastic modulus and the fillers, wherein the experimental result includes a first elastic modulus and rock porosity of the rock sample under the action of the fillers of each pore type, and a second elastic modulus of the rock sample under the action of all the fillers;

[0007] For the experiment after the first experiment, the rock sample after all the fillers in the previous experiment is used as the rock sample of the current experiment, the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment, and the operation of applying fillers of at least two pore types to the rock sample is repeatedly performed, and the experimental result corresponding to the rock sample is determined based on the initial elastic modulus and the fillers;

[0008] Based on the experimental results corresponding to the multiple experiments, the influence properties of the different fillers on the elastic modulus of the rock sample are determined.

[0009] According to another aspect of the present invention, there is provided a device for analyzing rock elastic modulus, the device comprising:

[0010] A first experiment module, used for obtaining a rock sample to be analyzed, for the first experiment, applying fillers of at least two pore types to the rock sample, and determining an experimental result corresponding to the rock sample based on an initial elastic modulus and the fillers, wherein the experimental result includes a first elastic modulus and rock porosity of the rock sample under the action of the fillers of each pore type, and a second elastic modulus of the rock sample under the action of all the fillers;

[0011] a subsequent experiment module, for an experiment after the first experiment, using the rock sample after all the fillers in the previous experiment as the rock sample of the current experiment, using the second elastic modulus in the previous experiment as the initial elastic modulus of the current experiment, repeatedly performing the operation of applying fillers of at least two pore types to the rock sample, and determining an experimental result corresponding to the rock sample based on the initial elastic modulus and the fillers;

[0012] The analysis module is used to determine the influence properties of different fillers on the elastic modulus of the rock sample based on the experimental results corresponding to multiple experiments.

[0013] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0014] at least one processor; and

[0015] a memory communicatively connected to the at least one processor; wherein,

[0016] The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for analyzing rock elastic modulus described in any embodiment of the present invention.

[0017] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for analyzing rock elastic modulus described in any embodiment of the present invention when executed.

[0018] The technical solution of the embodiment of the present invention is to obtain a rock sample to be analyzed, and for the first experiment, at least two pore types of fillers are applied to the rock sample, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the fillers, wherein the experimental results include the first elastic modulus and rock porosity of the rock sample under the action of each pore type of filler, and the second elastic modulus of the rock sample under the action of all fillers; for the experiment after the first experiment, the rock sample after all fillers in the previous experiment is used as the rock sample of the current experiment, and the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment, and the at least two pore types of fillers are applied to the rock sample. The invention relates to an operation of applying a filler of a certain type to a rock sample, and determining the experimental results corresponding to the rock sample based on the initial elastic modulus and the filler; determining the influence properties of different fillers on the elastic modulus of the rock sample based on the experimental results corresponding to multiple experiments, solving the problem of poor elastic parameter analysis results caused by the experimental method of adding a single pore type each time in the prior art, and realizing the calculation accuracy of the elastic modulus of the rock by adding fillers of different pore types to the rock sample at the same time and in multiple times, and taking the rock sample after all the fillers in the previous experiment as the rock sample of the current experiment, and taking the second elastic modulus in the previous experiment as the initial elastic modulus of the current experiment. Then, based on the first elastic modulus and rock porosity of the rock sample under the action of each pore type filler obtained in each experiment, the influence properties of different fillers on the elastic modulus of the rock sample are determined, achieving the technical effect of improving the accuracy of elastic parameter analysis and reducing analysis costs.

[0019] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0021] Figure 1 This is a flow chart of a method for analyzing rock elastic modulus provided according to Embodiment 1 of the present invention;

[0022] Figure 2 is a schematic diagram of characterizing the change of rock bulk modulus with porosity according to the first embodiment of the present invention;

[0023] Figure 3 is a schematic diagram of characterizing the change of rock shear modulus with porosity provided according to the first embodiment of the present invention;

[0024] Figure 4 is a schematic diagram of a method for analyzing rock elastic modulus provided according to the second embodiment of the present invention;

[0025] Figure 5 is a schematic structural diagram of a device for analyzing rock elastic modulus provided according to Embodiment 3 of the present invention;

[0026] Figure 6 It is a schematic diagram of the structure of an electronic device for implementing the method for analyzing rock elastic modulus according to an embodiment of the present invention. DETAILED DESCRIPTION

[0027] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0028] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0029] Embodiment 1

[0030] Figure 1This is a flow chart of a method for analyzing rock elastic modulus according to the first embodiment of the present invention. This embodiment is applicable to rock physics modeling. The method can be performed by a device for analyzing rock elastic modulus. The device for analyzing rock elastic modulus can be implemented in the form of hardware and / or software. The device for analyzing rock elastic modulus can be configured in a computing device. Figure 1 As shown, the method includes:

[0031] S110, obtaining a rock sample to be analyzed, and for a first experiment, applying fillers of at least two pore types to the rock sample, and determining an experimental result corresponding to the rock sample based on an initial elastic modulus and the fillers.

[0032] Among them, the rock sample can be a rock matrix, for example, it can be a carbonate rock composed of quartz, calcite, dolomite, clay minerals, or other types of rocks. The pore type can characterize the pore structure, such as coin-shaped, disc-shaped, needle-shaped, sphere-shaped, etc. The filler can be air, water, oil, etc. The initial elastic modulus during the first experiment can be the elastic modulus of the untreated rock sample, including shear modulus and bulk modulus. The experimental results include the first elastic modulus and rock porosity of the rock sample under the action of each pore type filler, and the second elastic modulus of the rock sample under the action of all fillers. The first elastic modulus and the second elastic modulus both include two types of elastic property information: shear modulus and bulk modulus.

[0033] In practical applications, after selecting the rock sample to be analyzed, the initial elastic modulus of the rock sample can be determined based on the rock attribute information of the rock sample, or an external force is applied to the rock sample, and the initial elastic modulus is determined based on the external force and the deformation of the rock sample. During the first experiment, at least two types of pore fillers can be added to the rock sample at the same time. It can be seen that the rock sample will produce different changes under the action of fillers of different pore types, and will eventually change under the action of fillers of all different pore types. Based on the initial elastic modulus and the information of the added fillers of various pore types, the elastic modulus (i.e., the first elastic modulus) and rock porosity of the rock sample under the action of each pore type of filler can be determined, as well as the elastic modulus (i.e., the second elastic modulus) of the rock sample under the action of all fillers. The filler information can characterize the attribute characteristics of the filler, such as mass, elastic modulus, volume ratio with the rock sample, mass ratio, etc.

[0034] In this embodiment, in the process of applying at least two pore types of fillers to a rock sample and determining experimental results corresponding to the rock sample based on the initial elastic modulus and the fillers, at least two pore types of fillers can be applied to the rock sample to determine the initial elastic modulus corresponding to the rock sample and filling characteristic information of each pore type of filler; and the experimental results are obtained based on the initial elastic modulus and the filling characteristic information of each pore type of filler.

[0035] Among them, the filling characteristic information includes filling elastic modulus, filling porosity and pore aspect ratio. The filling elastic modulus includes the shear modulus and bulk modulus of the filler. The pore aspect ratio can be used to characterize the morphological characteristics of rock pores, and its value can be between 0 and 1. In practical applications, after selecting the filler to be added to the rock sample, the filling characteristic information such as the filling elastic modulus, filling porosity and pore aspect ratio corresponding to the filler can be determined based on the filler's own information and the rock sample's own information. For example, the filling porosity corresponding to the filler can be the quotient of the volume of the filler divided by the volume of the rock sample. Further, when the first elastic modulus of the rock sample under the action of each pore type of filler is obtained based on the initial elastic modulus and the filling characteristic information of the filler of each pore type, the implementation method can be: for each pore type of filler, based on the initial elastic modulus, the filling elastic modulus corresponding to the filler, the filling porosity and the pore aspect ratio, the first elastic modulus of the rock sample under the action of the filler is determined. It should be noted that the method for determining the first elastic modulus of the rock sample under the action of each pore type filler is the same, and the method for determining the first elastic modulus of the rock sample under the action of any pore type filler can be used as an example for introduction.

[0036] In this embodiment, when determining the first elastic modulus of the rock sample under the action of any filler, the elastic modulus formula for calculating the dry rock skeleton can be used to determine the first elastic modulus of the rock sample under the action of the filler based on the initial elastic modulus, the filling elastic modulus corresponding to the filler, the filling porosity and the pore aspect ratio. The geometric shape factor of the filler can also be determined based on the filling elastic modulus, the filling porosity and the pore aspect ratio of the filler; the first elastic modulus of the rock sample under the action of the filler can be determined based on the initial elastic modulus, the filling elastic modulus, the geometric shape factor and the filling porosity. Among them, the geometric shape factor can be used to characterize the ability of the rock to resist crack propagation under the action of external force.

[0037] In practical applications, the geometric shape factor corresponding to the filler under the shear modulus can be determined based on the filler's shear modulus, filling porosity, and pore aspect ratio; the geometric shape factor corresponding to the filler under the bulk modulus can also be determined based on the filler's bulk modulus, filling porosity, and pore aspect ratio. For example, based on a geometric shape factor calculation function corresponding to the filler, the geometric shape factor corresponding to the filler can be solved according to the elastic modulus, filling porosity, and pore aspect ratio of the filler.

[0038] For example, the geometric shape factor R corresponding to the bulk modulus of the coin-shaped gap type filler is mi The geometric shape factor S corresponding to the shear modulus mi The calculation function is as follows:

[0039]

[0040]

[0041] The geometric shape factor R corresponding to the bulk modulus of the filler of the disc-like pore type mi The geometric shape factor S corresponding to the shear modulus mi The calculation function is as follows:

[0042]

[0043]

[0044] The geometric shape factor R corresponding to the bulk modulus of the needle-like pore type filler mi The geometric shape factor S corresponding to the shear modulus mi The calculation function is as follows:

[0045]

[0046]

[0047] The geometric shape factor R corresponding to the bulk modulus of the spherical pore type filling mi The geometric shape factor S corresponding to the shear modulus mi The calculation function is as follows:

[0048]

[0049]

[0050] Among them, A m Expressed as the bulk modulus of the rock matrix m, B m Expressed as the shear modulus of the rock matrix, A i represents the bulk modulus of the i-th filler, Bi represents the shear modulus of the i-th filler, and α is the aspect ratio of the fracture pore. In this item, the A and B superscripts are consistent with the β superscript. In this item, the A and B superscripts are consistent with the β superscript. In this item, the A and B superscripts are consistent with the ξ superscript.

[0051] In this embodiment, the initial elastic modulus, the filling elastic modulus of the filler, the geometric shape factor and the filling porosity can be input into the elastic modulus calculation function to determine the first elastic modulus of the rock sample under the action of the filler.

[0052] For example, the calculation function of the bulk modulus in the elastic modulus formula may be:

[0053]

[0054] Among them, φ i is the porosity of the i-th filler, It is expressed as the bulk modulus of the rock sample after the action of the i-th filler, and N represents the total number of fillers.

[0055] The calculation function of the shear modulus in the elastic modulus formula can be:

[0056]

[0057] in, It is expressed as the shear modulus of the rock sample after the action of the i-th filling material.

[0058] In this embodiment, the implementation method of determining the first elastic modulus of the rock sample under the action of the filler based on the initial elastic modulus, the filling elastic modulus, the geometric shape factor and the filling porosity can also be: based on the filling porosity, determining the filling ratio corresponding to the filler; based on the filling ratio, the initial elastic modulus, the filling elastic modulus and the geometric shape factor, determining the first elastic modulus of the rock sample under the action of the filler.

[0059] Among them, the filling ratio can be used to characterize the volume ratio of the filling material to the rock sample.

[0060] In this embodiment, when a filler is added to a rock sample, the filling ratio can be calculated based on the ratio of the content of the added filler to the total content of the rock sample; the filling porosity of the filler can also be used as the filling ratio. Further, the first elastic modulus of the rock sample under the action of the filler can be determined based on the calculation function corresponding to the elastic modulus, the filling ratio, the initial elastic modulus, the filling elastic modulus and the geometric shape factor.

[0061] Exemplarily, the calculation function of the bulk modulus is shown in formula (1), and the calculation function of the shear modulus is shown in formula (2),

[0062] Formula (1):

[0063] Formula (2):

[0064] Wherein, x represents the filling ratio and h represents the number of times the filler is added.

[0065] It should be noted that when a filler with a volume content ratio of x is added to the rock sample in each experiment, if the filler is air, that is, pores are added to the rock sample, in experiments other than the first experiment, the newly added filler replaces not only the rock matrix, but also part of the previously added filler. Therefore, the filler added to the rock matrix each time is not the actual filler volume content dφ, and only (1-φ)x contributes to the actual filler volume content component dφ during the iterative experiment. φ is the porosity of the rock matrix. The change in the actual filler volume content during each iteration is:

[0066] dφ=(1-φ)x; The relationship between the actual filling volume content in the rock and x is: φ=1-e -∑x ;∑x is the total content of filler. If the number of iterative experiments is H, then φ=1-e Hx .

[0067] Based on this, based on the filling porosity, the implementation method of determining the filling ratio of the filler can be: based on the filling porosity and a preset first numerical value, determine a first intermediate value; perform logarithmic processing on the first intermediate value to obtain a second intermediate value; determine a preset number of experiments for the experiment; based on the second intermediate value and the preset number of experiments, obtain the filling ratio of the filler.

[0068] The preset first value may be 1. The preset number of experiments may be a preset number of experiments, such as 10 times, which requires ten iterations of experiments. The specific value may be determined by technicians according to actual working conditions and is not limited here.

[0069] In practical applications, the preset first value and the filling porosity can be subtracted to obtain a difference, and then the preset first value and the difference can be processed as a quotient, and the obtained quotient can be used as a first intermediate value. The first intermediate value is processed logarithmically, and the obtained value can be used as a second intermediate value. The second intermediate value and the preset number of experiments are processed as a quotient, and the obtained quotient is used as the filling ratio of the filler.

[0070] For example, the filling ratio x may be determined based on formula (3):

[0071] Formula (3) can be expressed as:

[0072] Where H represents the preset number of experiments. φ represents the filling porosity, which represents the volume content of the actual filling material in the rock sample.

[0073] In this embodiment, based on the initial elastic modulus and the filling characteristic information of the fillers of each pore type, the second elastic modulus of the rock sample under the action of all fillers is obtained by: for each filler, based on the filling characteristic information and the initial elastic modulus of the filler, determining the third intermediate value corresponding to the filler; determining the fourth intermediate value based on the third intermediate values ​​corresponding to all fillers and the preset second value; determining the second elastic modulus based on the initial elastic modulus and the fourth intermediate value.

[0074] The preset second value is 1. The second elastic modulus includes a bulk modulus and a shear modulus. The implementation method of determining the third intermediate value corresponding to each filler is the same, and the method of determining the third intermediate value corresponding to any filler is described as an example.

[0075] In this embodiment, the third intermediate value of the filler under the bulk modulus can be determined according to the first function based on the initial elastic modulus and the bulk modulus of the filler, the geometric shape factor under the bulk modulus and the filling porosity. The first function can be: The third intermediate value of the filler under shear modulus can be determined according to the second function based on the initial elastic modulus and the shear modulus of the filler, the geometric shape factor under shear modulus and the filler porosity. The second function can be: The third intermediate values ​​corresponding to all fillers are summed, the sum is added to the preset second value, and the negative first power of the obtained value is used as the fourth intermediate value. The product of the initial elastic modulus and the fourth intermediate value is used as the second elastic modulus.

[0076] For example, the bulk modulus of the rock sample under the action of the filler can be calculated based on formula (1), and the shear modulus of the rock sample under the action of the filler can be calculated based on formula (2). The bulk modulus of the rock sample under the action of all fillers can be determined based on the calculation function (4) of the bulk modulus, and the shear modulus of the rock sample under the action of all fillers can be determined based on the calculation function (5) of the shear modulus.

[0077] Among them, the calculation function (4) can be expressed as

[0078] The calculation function (5) can be expressed as

[0079] S120. For the experiment after the first experiment, the rock sample after all the fillers in the previous experiment is used as the rock sample of the current experiment, and the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment. The operation of applying fillers of at least two pore types to the rock sample is repeated, and the experimental result corresponding to the rock sample is determined based on the initial elastic modulus and the fillers.

[0080] In this embodiment, when conducting an experiment after the first experiment, the rock sample after all the fillers in the previous experiment can be used as the rock sample of the current experiment, and the second elastic modulus obtained in the previous experiment can be used as the initial elastic modulus of the rock sample of the current experiment. In this way, during the current experiment, the operation of applying fillers of at least two pore types to the rock sample and determining the experimental results corresponding to the rock sample based on the initial elastic modulus and the fillers is repeated. Accordingly, the experimental results obtained after multiple experiments can be obtained. For example, in the second round of experiments, the rock sample A after all the fillers in the first experiment is used as the rock sample B of the second experiment, and the second elastic modulus C obtained in the first experiment is used as the initial elastic modulus D of the rock sample in the second experiment. At least two pore types of fillers can be applied to the rock sample B, and the experimental results corresponding to the rock sample B can be determined based on the initial elastic modulus D and the fillers.

[0081] Exemplarily, the bulk modulus of the rock sample under the action of all fillers can be determined based on the calculation function (6) of the bulk modulus, and the shear modulus of the rock sample under the action of all fillers can be determined based on the calculation function (7) of the shear modulus:

[0082]

[0083]

[0084] S130. Based on the experimental results corresponding to the multiple experiments, determine the influence properties of different fillers on the elastic modulus of the rock sample.

[0085] In this embodiment, after multiple experiments, the first elastic modulus and rock porosity of the rock sample under the action of each pore type filler in each experiment can be obtained. The influence of different pore shapes on the elastic properties of the rock can be analyzed based on the change of the first elastic modulus with the change of rock porosity.

[0086] For example, see Figure 2 , which can characterize the change of rock bulk modulus with porosity when different pore types of fillers are added. The pore types are: fracture pores, spherical pores, convex triangular pores, and concave quadrilateral pores. Figure 3 , which can characterize the change of rock shear modulus with porosity when different pore types of fillers are added. It can be seen that different pore types have an impact on the rock elastic modulus, especially cracks, whose content changes have the greatest impact on the rock elastic properties.

[0087] The technical solution of this embodiment is to obtain a rock sample to be analyzed, and for the first experiment, at least two pore types of fillers are applied to the rock sample, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the fillers, and the experimental results include the first elastic modulus and rock porosity of the rock sample under the action of each pore type of filler, and the second elastic modulus of the rock sample under the action of all fillers; for the experiment after the first experiment, the rock sample after all fillers in the previous experiment is used as the rock sample of the current experiment, and the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment, and the at least two pore types are repeatedly applied to the rock sample. The filler acts on the rock sample, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the filler; based on the experimental results corresponding to multiple experiments, the influence properties of different fillers on the elastic modulus of the rock sample are determined, which solves the problem of poor elastic parameter analysis results caused by the experimental method of adding a single pore type each time in the prior art, and realizes that the calculation accuracy of the elastic modulus of the rock is improved by adding fillers of different pore types to the rock sample at the same time and in multiple times, and taking the rock sample after all the fillers in the previous experiment as the rock sample of the current experiment, and taking the second elastic modulus in the previous experiment as the initial elastic modulus of the current experiment. Then, based on the first elastic modulus and rock porosity of the rock sample under the action of each pore type of filler obtained in each experiment, the influence properties of different fillers on the elastic modulus of the rock sample are determined, so as to achieve the technical effect of improving the accuracy of elastic parameter analysis and reducing analysis costs.

[0088] Embodiment 2

[0089] As an alternative embodiment of the above embodiment, Figure 4 Schematic diagram of a method for analyzing rock elastic modulus provided according to Embodiment 2 of the present invention. For details, please refer to the following specific contents.

[0090] See also Figure 4, the upper and lower limits of the elastic modulus of the rock matrix (i.e., rock sample) can be obtained by the Voigt and Reuss method, and then the elastic modulus of the final rock matrix can be obtained by averaging as the initial elastic modulus. Furthermore, the constructed rock physics model can be used to simulate the addition of fillers (such as fracture pores) to the rock matrix, and the newly obtained porous rock can be used as the new rock matrix. The elastic modulus of the filler fluid can be calculated using the Brie model, and then the Gassmann equation can be used to introduce the filler fluid into the dry rock skeleton to simulate the internal structure of the rock under the actual formation environment, and the elastic modulus of the rock matrix after the action can be calculated. This method not only considers the influence of fillers on the rock matrix, but also conforms to the internal pore structure and distribution of rocks in actual reservoirs.

[0091] In this embodiment, the rock physics model can be constructed in the following manner: considering the influence of various pore type fillers on rock elastic parameters, the rock physics model can be constructed by adding fillers of different pore types to the rock matrix multiple times and simultaneously. Based on the rock physics model, an analysis experiment of rock elastic parameters is conducted to obtain the first elastic modulus and rock porosity of the rock sample under the action of each pore type filler in each experiment, and the second elastic modulus of the rock sample under the action of all fillers, so as to determine the influence of different fillers on the elastic modulus of the rock sample based on the first elastic modulus and rock porosity obtained from multiple experiments.

[0092] Exemplarily, the calculation function of the physical model for determining the bulk modulus of a rock sample under the action of all fillers may be:

[0093]

[0094] The calculation function of the physical model for determining the shear modulus of a rock sample under the action of all fillings can be:

[0095]

[0096] Among them, A m is the bulk modulus of the rock sample, B m is the shear modulus of the rock sample, A i Expressed as the bulk modulus of the i-th filler, B i Expressed as the shear modulus of the i-th filler, R mi is the geometric shape factor of the bulk modulus of the i-th filler in the rock sample, S mi is the geometric shape factor of the shear modulus of the i-th filling in the rock sample, φ i is the porosity of the i-th filler, is the calculated bulk modulus of the rock sample after the action, It is the shear modulus of the rock sample calculated after the action.

[0097] By introducing the idea of ​​differential equivalence, a small amount of pores (i.e. fillers) can be added to the rock sample each time to make it meet the low porosity condition (the ratio of the porosity of the rock sample to the pore aspect ratio is much less than 1). The new equivalent medium after the pore effect is regarded as a new rock sample, and the process is iterated multiple times to calculate the equivalent elastic modulus of the rock at high porosity.

[0098] During multiple experiments, the calculation function of the physical model for determining the bulk modulus of rock samples under the action of all fillers can be:

[0099]

[0100] During multiple experiments, the calculation function of the physical model for determining the shear modulus of the rock sample under the action of all filling materials can be:

[0101]

[0102] A h is the bulk modulus of the rock sample obtained after h iterations of the experiment, B h is the shear modulus of the rock sample obtained after h iterations of the experiment, is the bulk modulus of the rock sample after h+1 iterations of experiments, is the shear modulus of the rock sample after h+1 iterations of the experiment.

[0103] It should be noted that if the filler content (i.e., filling ratio) added each time is x, when the filler is air, it is added to the pores. In the experiments except the first experiment, the newly added filler replaces not only the rock sample, but also part of the filler. Therefore, the filler added to the rock matrix each time is not the actual filler volume content dφ, only (1-φ)x of the rock sample is replaced, so dφ=(1-φ)x, and the relationship between the actual filler volume content in the rock and x is: φ=1-e -∑x ;∑x is the total content of filler. If the number of iterative experiments is H, then φ=1-e Hx . x can be written as: Substitute x for φ in the calculation function i , the calculation function of the physical model of the bulk modulus of the rock sample under the action of the filler after adding a single pore type filler is obtained:

[0104]

[0105] The calculation function of the physical model for determining the shear modulus of the rock sample under the action of the filler is:

[0106]

[0107] The technical solution of this embodiment constructs a rock physics model by simulating the addition of porous rock porosity and pore structure to rock samples, and then analyzes the influence of complex pore and fracture morphology on the elastic modulus of elastic waves of fluid-saturated reservoirs through the rock physics model. While truly simulating the internal structure of the rock formation, it also improves the accuracy of complex pore and fracture reservoir evaluation, provides an important theoretical basis for the design of oil and gas reservoir development plans and successful trial production, and provides reliable technical support for reservoir oil and gas exploration and development.

[0108] Embodiment 3

[0109] Figure 5 FIG. 1 is a schematic diagram of a structure of a device for analyzing rock elastic modulus according to Embodiment 3 of the present invention. Figure 5 As shown, the device includes: a first experiment module 210, a subsequent experiment module 220 and an analysis module 230.

[0110] Among them, the first experiment module 210 is used to obtain a rock sample to be analyzed, and for the first experiment, at least two pore types of fillers are applied to the rock sample, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the fillers, wherein the experimental results include the first elastic modulus and rock porosity of the rock sample under the action of the fillers of each pore type, and the second elastic modulus of the rock sample under the action of all the fillers; the subsequent experiment module 220 is used for the experiment after the first experiment, the rock sample after all the fillers in the previous experiment is used as the rock sample of the current experiment, the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment, and the operation of applying at least two pore types of fillers to the rock sample is repeatedly performed, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the fillers; the analysis module 230 is used to determine the influence attributes of different fillers on the elastic modulus of the rock sample based on the experimental results corresponding to multiple experiments.

[0111] The technical solution of this embodiment is to obtain a rock sample to be analyzed, and for the first experiment, at least two pore types of fillers are applied to the rock sample, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the fillers, and the experimental results include the first elastic modulus and rock porosity of the rock sample under the action of each pore type of filler, and the second elastic modulus of the rock sample under the action of all fillers; for the experiment after the first experiment, the rock sample after all fillers in the previous experiment is used as the rock sample of the current experiment, and the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment, and the at least two pore types are repeatedly applied to the rock sample. The filler acts on the rock sample, and the experimental results corresponding to the rock sample are determined based on the initial elastic modulus and the filler; based on the experimental results corresponding to multiple experiments, the influence properties of different fillers on the elastic modulus of the rock sample are determined, which solves the problem of poor elastic parameter analysis results caused by the experimental method of adding a single pore type each time in the prior art, and realizes that the calculation accuracy of the elastic modulus of the rock is improved by adding fillers of different pore types to the rock sample at the same time and in multiple times, and taking the rock sample after all the fillers in the previous experiment as the rock sample of the current experiment, and taking the second elastic modulus in the previous experiment as the initial elastic modulus of the current experiment. Then, based on the first elastic modulus and rock porosity of the rock sample under the action of each pore type of filler obtained in each experiment, the influence properties of different fillers on the elastic modulus of the rock sample are determined, so as to achieve the technical effect of improving the accuracy of elastic parameter analysis and reducing analysis costs.

[0112] On the basis of the above-mentioned device, optionally, the first experiment module 210 includes a filling feature information determination unit and an experiment result determination unit.

[0113] A filling characteristic information determining unit, used for applying fillers of at least two pore types to the rock sample to determine an initial elastic modulus corresponding to the rock sample and filling characteristic information of the fillers of each pore type;

[0114] The experimental result determination unit is used to obtain the experimental result based on the initial elastic modulus and the filling characteristic information of the filler of each pore type.

[0115] Based on the above device, optionally, the filling characteristic information includes filling elastic modulus, filling porosity and pore aspect ratio; and the experimental result determination unit includes a first result determination unit.

[0116] The result determines a first unit, which is used to determine the first elastic modulus of the rock sample under the action of the filler for each pore type of the filler based on the initial elastic modulus, the filling elastic modulus corresponding to the filler, the filling porosity and the pore aspect ratio.

[0117] Based on the above device, optionally, the result determining first unit includes a geometric shape factor determining unit and a first elastic modulus determining unit.

[0118] a geometric shape factor determination unit, configured to determine a geometric shape factor of the filler based on the fill elastic modulus, the fill porosity and the pore aspect ratio;

[0119] The first elastic modulus determination unit is used to determine the first elastic modulus of the rock sample under the action of the filler based on the initial elastic modulus, the filling elastic modulus, the geometric shape factor and the filling porosity.

[0120] On the basis of the above device, optionally, the first elastic modulus determining unit includes a filling ratio determining subunit and a first elastic modulus determining subunit.

[0121] A filling ratio determination subunit, used to determine the filling ratio corresponding to the filler based on the filling porosity;

[0122] The first elastic modulus determination subunit is used to determine the first elastic modulus of the rock sample under the action of the filler based on the filling ratio, the initial elastic modulus, the filling elastic modulus and the geometric shape factor.

[0123] On the basis of the above device, optionally, the filling ratio determination subunit includes a first intermediate value determination subunit, a second intermediate value determination subunit, a preset number of experiments determination subunit and a filling ratio determination subunit.

[0124] A first intermediate value determination subunit is used to determine a first intermediate value based on the filling porosity and a preset first value;

[0125] A second intermediate value determination subunit is used to perform logarithmic processing on the first intermediate value to obtain a second intermediate value;

[0126] The preset number of experiments determines the small unit, which is used to determine the preset number of experiments to be conducted;

[0127] The filling ratio determination subunit is used to obtain the filling ratio of the filler based on the second intermediate value and the preset number of experiments.

[0128] Based on the above device, optionally, the experimental result determination unit includes a second result determination unit. The second result determination unit includes a third intermediate value determination unit, a fourth intermediate value determination unit and a second elastic modulus determination unit.

[0129] A third intermediate value determining unit, configured to determine, for each of the fillers, a third intermediate value corresponding to the filler based on the filling characteristic information of the filler and the initial elastic modulus;

[0130] a fourth intermediate value determining unit, configured to determine a fourth intermediate value based on the third intermediate values ​​corresponding to all the fillers and a preset second value;

[0131] The second elastic modulus determining unit is configured to determine the second elastic modulus based on the initial elastic modulus and the fourth intermediate value.

[0132] The device for analyzing the elastic modulus of rock provided in the embodiment of the present invention can execute the method for analyzing the elastic modulus of rock provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0133] Example 6

[0134] Figure 6 1 is a schematic diagram of the structure of an electronic device for implementing the method for analyzing the elastic modulus of rock according to an embodiment of the present invention. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workbenches, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0135] like Figure 6As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0136] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0137] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a method for analyzing the elastic modulus of rock.

[0138] In some embodiments, the method for analyzing the elastic modulus of rock can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the method for analyzing the elastic modulus of rock described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to perform the method for analyzing the elastic modulus of rock in any other appropriate manner (for example, by means of firmware).

[0139] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0140] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0141] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, device, or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices, or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0142] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0143] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0144] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0145] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0146] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for analyzing rock elastic modulus, characterized in that: include: Obtain a rock sample to be analyzed, and for a first experiment, apply fillers of at least two pore types to the rock sample, and determine an experimental result corresponding to the rock sample based on an initial elastic modulus and the fillers, wherein the experimental result includes a first elastic modulus and rock porosity of the rock sample under the action of the fillers of each pore type, and a second elastic modulus of the rock sample under the action of all the fillers; For the experiment after the first experiment, the rock sample after all the fillers in the previous experiment is used as the rock sample of the current experiment, the second elastic modulus in the previous experiment is used as the initial elastic modulus of the current experiment, and the operation of applying fillers of at least two pore types to the rock sample is repeatedly performed, and the experimental result corresponding to the rock sample is determined based on the initial elastic modulus and the fillers; Based on the experimental results corresponding to the multiple experiments, the influence properties of the different fillers on the elastic modulus of the rock sample are determined.

2. The method according to claim 1, characterized in that The step of applying fillers of at least two pore types to the rock sample and determining an experimental result corresponding to the rock sample based on an initial elastic modulus and the fillers comprises: Apply fillers of at least two pore types to the rock sample to determine an initial elastic modulus corresponding to the rock sample and filling characteristic information of the fillers of each pore type; The experimental results are obtained based on the initial elastic modulus and the filling characteristic information of the filler of each pore type.

3. The method according to claim 2, characterized in that The filling characteristic information includes filling elastic modulus, filling porosity and pore aspect ratio; the first elastic modulus of the rock sample under the action of the filler of each pore type is obtained based on the initial elastic modulus and the filling characteristic information of the filler of each pore type, including: For the filler of each pore type, the first elastic modulus of the rock sample under the action of the filler is determined based on the initial elastic modulus, the filling elastic modulus corresponding to the filler, the filling porosity and the pore aspect ratio.

4. The method according to claim 3, characterized in that The determining of a first elastic modulus of the rock sample under the action of the filler based on the initial elastic modulus, the filling elastic modulus corresponding to the filler, the filling porosity and the pore aspect ratio comprises: Determining a geometric shape factor of the filler based on the fill elastic modulus, fill porosity, and pore aspect ratio; A first elastic modulus of the rock sample under the action of the filler is determined based on the initial elastic modulus, the filling elastic modulus, the geometric shape factor and the filling porosity.

5. The method according to claim 4, characterized in that The determining, based on the initial elastic modulus, the filling elastic modulus, the geometric shape factor and the filling porosity, of a first elastic modulus of the rock sample under the action of the filler comprises: Based on the filling porosity, determining a filling ratio corresponding to the filler; Based on the filling ratio, the initial elastic modulus, the filling elastic modulus and the geometric shape factor, a first elastic modulus of the rock sample under the action of the filler is determined.

6. The method according to claim 5, characterized in that The determining of the filling ratio corresponding to the filler based on the filling porosity includes: Determining a first intermediate value based on the filling porosity and a preset first value; Performing logarithmic processing on the first intermediate value to obtain a second intermediate value; Determine the preset number of trials for conducting the experiment; Based on the second intermediate value and the preset number of experiments, a filling ratio of the filler is obtained.

7. The method according to claim 2, characterized in that The step of obtaining a second elastic modulus of the rock sample under the action of all the fillers based on the initial elastic modulus and the filling characteristic information of the fillers of each pore type comprises: For each of the fillers, based on the filling characteristic information of the filler and the initial elastic modulus, determining a third intermediate value corresponding to the filler; Determining a fourth intermediate value based on the third intermediate values ​​corresponding to all the fillers and the preset second value; Based on the initial elastic modulus and the fourth intermediate value, the second elastic modulus is determined.

8. A device for analyzing rock elastic modulus, characterized in that: include: A first experiment module, used for obtaining a rock sample to be analyzed, for the first experiment, applying fillers of at least two pore types to the rock sample, and determining an experimental result corresponding to the rock sample based on an initial elastic modulus and the fillers, wherein the experimental result includes a first elastic modulus and rock porosity of the rock sample under the action of the fillers of each pore type, and a second elastic modulus of the rock sample under the action of all the fillers; a subsequent experiment module, for an experiment after the first experiment, using the rock sample after all the fillers in the previous experiment as the rock sample of the current experiment, using the second elastic modulus in the previous experiment as the initial elastic modulus of the current experiment, repeatedly performing the operation of applying fillers of at least two pore types to the rock sample, and determining an experimental result corresponding to the rock sample based on the initial elastic modulus and the fillers; The analysis module is used to determine the influence properties of different fillers on the elastic modulus of the rock sample based on the experimental results corresponding to multiple experiments.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can perform the method for analyzing rock elastic modulus as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for analyzing rock elastic modulus according to any one of claims 1 to 7 when executed.

Images