Mature shale anisotropy analysis method, electronic equipment, storage medium and device

Abstract

The invention discloses a mature shale anisotropy analysis method, electronic equipment, a storage medium and a device. The method comprises the following steps: constructing a mature shale rock physical model, and obtaining rock equivalent elastic parameters; establishing a mature shale anisotropy interpretation quantity plate based on the rock equivalent elastic parameters; carrying out intersection analysis on the parallel bedding stiffness and the vertical bedding stiffness of the mature shale based on a mature shale anisotropy interpretation measuring plate; mature shale variable confining pressure elastic parameter data measured in a laboratory are input into the mature shale anisotropy interpretation measurement plate, and anisotropy interpretation analysis of the mature shale is completed. According to the method, the mature shale anisotropy interpretation measuring plate is established, changes, caused by stress changes, of rock elasticity moduli in different directions are analyzed according to mineralogy and acoustic data, the source of shale stratum stress-related rigidity changes and anisotropy of media are evaluated, a proper rock skeleton model can be selected easily, and the rock stratum stress-related rigidity change source evaluation method can be used for evaluating the rock stratum stress-related rigidity change source evaluation method. And the hydraulic fracturing design is optimized.

Description

Technical Field

[0001] The present invention belongs to the field of geophysical technologies, and more specifically, relates to a method for analyzing anisotropy of mature shale, an electronic device, a storage medium, and a device. Background Art

[0002] According to the different degrees of kerogen maturity, shale can be divided into immature, mature, and overmature shale. Among them, kerogen in immature shale plays a role in supporting the load; kerogen in mature shale also supports the load, but there are oil-bearing pores developed in the kerogen; while in overmature shale, kerogen mainly exists in the form of inclusions in the rock skeleton. The change in the occurrence characteristics of kerogen makes shales with different maturities exhibit different elastic and anisotropic characteristics. For tight unconventional shale formations, their seismic anisotropy is generated by shale intrinsic anisotropy and external anisotropy. Shale intrinsic anisotropy includes thin layers, lenticular clay layers deposited nearly parallel to the bedding plane, as well as oriented organic matter, pores, and fractures, etc. The external anisotropy of shale comes from the oriented fractures generated by drilling or coring operations. The anisotropic characteristics of shale rocks will undergo mechanical changes (such as the compression of soft components arranged in an oriented manner and the closure of fractures) due to pressure loss during the development process, resulting in artificially induced changes in elastic and fluid transport characteristics. Rock physics models can be used to estimate rock stress and direction-dependent elastic moduli, as well as horizontal stress. Currently, the interpretation work on the elastic anisotropy of mature shale still needs to be continuously improved. Although multiple directions of elastic stiffness are required to simulate the seismic response of shale, usually only the stiffness in one direction is available. When the selection of rock physics model parameters is not restricted, the prediction of in-situ seismic parameters has great uncertainty. Therefore, it is necessary to develop an anisotropy interpretation template for mature shale, which can analyze the change in rock elastic modulus in different directions due to stress changes based on mineralogical and acoustic data, and evaluate the source of stress-related stiffness changes in the shale formation and the anisotropy of the medium. By deepening the understanding of the shale rock structure or symmetry, it helps to select an appropriate rock skeleton model and optimize the hydraulic fracturing design.

[0003] The information disclosed in the background art part of the present invention is only intended to deepen the understanding of the general background art of the present invention, and should not be regarded as an admission or any form of implication that this information constitutes the prior art known to those skilled in the art. Summary of the Invention

[0004] The objective of the present invention is to propose a method for analyzing anisotropy of mature shale, an electronic device, a storage medium, and a device, so as to be able to analyze the change in rock elastic modulus in different directions due to stress changes based on mineralogical and acoustic data, and evaluate the source of stress-related stiffness changes in the shale formation and the anisotropy of the medium.

[0005] To achieve the above object, the present invention provides a method for analyzing the anisotropy of mature shale, an electronic device, a storage medium, and a device.

[0006] According to a first aspect of the present invention, there is provided a method for analyzing the anisotropy of mature shale, comprising:

[0007] Constructing a petrophysical model of mature shale to obtain equivalent elastic parameters of the rock;

[0008] Based on the equivalent elastic parameters of the rock, establishing an anisotropy interpretation template for mature shale;

[0009] Based on the anisotropy interpretation template for mature shale, performing an intersection analysis on the stiffness parallel to the bedding and the stiffness perpendicular to the bedding of mature shale;

[0010] Inputting the data of the elastic parameters of mature shale under variable confining pressure measured in the laboratory into the anisotropy interpretation template for mature shale to complete the anisotropy interpretation analysis of mature shale.

[0011] Optionally, the constructing of the petrophysical model of mature shale includes:

[0012] Obtaining the equivalent elastic stiffness of the inorganic rock skeleton through anisotropic SCA;

[0013] Obtaining the equivalent elastic stiffness of the pore-containing kerogen through anisotropic DEM;

[0014] Based on the equivalent elastic stiffness of the inorganic rock skeleton and the equivalent elastic stiffness of the pore-containing kerogen, adding the pore-containing kerogen to the inorganic shale skeleton through anisotropic Backus averaging to construct the petrophysical model of mature shale.

[0015] Optionally, the expression of the anisotropic SCA is:

[0016]

[0017] where v r is the volume fraction of the r-th component (r = 1, 2,..., N) of the inorganic part of the shale, C r represents the stiffness tensor of the r-th component, C SCA represents the self-consistent equivalent elastic tensor of the inorganic skeleton, is a fourth-order tensor given by the stress Green's function integrated over the inclusion shape.

[0018] Optionally, the expression of the anisotropic DEM is:

[0019]

[0020] where C fluidRepresents the stiffness of hydrocarbons filled in kerogen-related pores, C DEM (v) represents the stiffness of kerogen with pores, where v is the content of kerogen-related pores, C DEM (0) represents the stiffness of solid kerogen (v = 0, the content of kerogen-related pores is 0).

[0021] Optionally, the expression of the anisotropic Backus average is:

[0022]

[0023]

[0024]

[0025]

[0026]

[0027] Where and are the equivalent stiffnesses calculated by Backus averaging for kerogen with pores and the inorganic shale skeleton, C 11 、C 13 、C 33 、C 44 and C 66 are the stiffnesses of kerogen with pores or the inorganic shale skeleton, <.> represents taking the weighted average of the properties in the brackets according to their volume ratios.

[0028] Optionally, the expressions of the parallel bedding stiffness and the perpendicular bedding stiffness are respectively:

[0029] C || =ρ×V || 2 ;

[0030]

[0031] Where C || is the parallel bedding stiffness of the rock, C ⊥ is the perpendicular bedding stiffness, ρ is the rock density, V || is the rock velocity parallel to the bedding, V ⊥ is the rock velocity perpendicular to the bedding.

[0032] Optionally, in an isotropic medium, DC || ≈DC ⊥ ;

[0033] In a VTI medium, DC || <DC ⊥ ;

[0034] In the HTI medium, DC || >DC ⊥ ;

[0035] wherein, DC || represents the change of elastic stiffness with stress in the direction parallel to the bedding plane, and DC ⊥ represents the change of elastic stiffness with stress in the direction perpendicular to the bedding plane.

[0036] According to the second aspect of the present invention, an anisotropy analysis device for mature shale is provided, including:

[0037] A construction module for constructing a petrophysical model of mature shale and obtaining equivalent elastic parameters of the rock;

[0038] An establishment module for establishing an anisotropy interpretation template for mature shale based on the equivalent elastic parameters of the rock;

[0039] An intersection analysis module for performing intersection analysis on the stiffness parallel to the bedding plane and the stiffness perpendicular to the bedding plane of mature shale based on the anisotropy interpretation template for mature shale;

[0040] An interpretation analysis module for inputting the data of the elastic parameters of mature shale under variable confining pressure measured in the laboratory into the anisotropy interpretation template for mature shale to complete the anisotropy interpretation analysis of mature shale.

[0041] According to the third aspect of the present invention, an electronic device is provided, the electronic device includes:

[0042] At least one processor; and,

[0043] A memory communicatively connected to the at least one processor; wherein,

[0044] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute any one of the mature shale anisotropy analysis methods in the first aspect.

[0045] According to the fourth aspect of the present invention, a non-transitory computer-readable storage medium is provided, and the non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute any one of the mature shale anisotropy analysis methods in the first aspect.

[0046] The beneficial effects of the present invention are as follows: By constructing a mature shale rock physics model, the present invention establishes an anisotropic interpretation template for mature shale, which can analyze the change in rock elastic modulus in different directions due to stress changes based on mineralogy and acoustic data, evaluate the source of stress-related stiffness changes in the shale formation, and the anisotropy of the medium.

[0047] The system of the present invention has other characteristics and advantages, which will be obvious from the accompanying drawings incorporated herein and the subsequent detailed description, or will be described in detail in the accompanying drawings incorporated herein and the subsequent detailed description. These accompanying drawings and detailed description are used together to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] By describing the exemplary embodiments of the present invention in more detail in conjunction with the accompanying drawings, the above and other objects, features, and advantages of the present invention will become more obvious. In the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

[0049] Figure 1 FIG. shows a flowchart of the steps of a method for analyzing the anisotropy of mature shale according to the present invention.

[0050] Figure 2 FIG. shows a flowchart of the steps of a method for analyzing the anisotropy of mature shale according to Embodiment 2 of the present invention.

[0051] Figure 3 FIG. shows a schematic diagram of the longitudinal wave anisotropy interpretation template of mature shale according to Embodiment 2 of the present invention.

[0052] Figure 4 FIG. shows a schematic diagram of the shear wave anisotropy interpretation template of mature shale according to Embodiment 2 of the present invention.

[0053] Figure 5 FIG. shows a schematic diagram of an apparatus for analyzing the anisotropy of mature shale according to Embodiment 3 of the present invention. DETAILED DESCRIPTION

[0054] The present invention will be described in more detail below with reference to the accompanying drawings. Although the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to make the present invention more thorough and complete, and to fully convey the scope of the present invention to those skilled in the art.

[0055] As Figure 1 shown, a method for analyzing the anisotropy of mature shale according to the present invention includes:

[0056] Construct a mature shale rock physics model to obtain the equivalent elastic parameters of the rock;

[0057] Establish an anisotropic interpretation template for mature shale based on the equivalent elastic parameters of the rock;

[0058] Conduct an intersection analysis of the stiffness parallel to the bedding and the stiffness perpendicular to the bedding of mature shale based on the anisotropic interpretation template for mature shale;

[0059] Input the data of the elastic parameters of mature shale under variable confining pressure measured in the laboratory into the anisotropic interpretation template for mature shale to complete the anisotropic interpretation analysis of mature shale.

[0060] Specifically, the present invention uses anisotropic SCA to mix quartz, calcite-like rock minerals, and oriented clay, etc. to obtain the equivalent elastic stiffness of the inorganic rock skeleton, uses anisotropic DEM to obtain the equivalent elastic stiffness of pore-containing kerogen, uses anisotropic Backus averaging to add the pore-containing kerogen to the shale skeleton to construct a mature shale rock physics model, and then uses the equivalent elastic parameters of the rock obtained from the mature shale rock physics model to establish an anisotropic interpretation template for mature shale; through the anisotropic interpretation template for mature shale, the stiffness parallel to the bedding (C || ) and the stiffness perpendicular to the bedding (C ⊥ ) of the shale are subjected to intersection analysis. Finally, the data of the elastic parameters of mature shale under variable confining pressure measured in the laboratory are projected onto the anisotropic interpretation template to complete the anisotropic interpretation analysis of mature shale. In an isotropic medium, the deformation of the rock in the direction parallel to the bedding || and the direction perpendicular to the bedding ⊥ is approximately equal, that is, DC || ≈DC ⊥ , where DC || and DC ⊥ respectively represent the change of the elastic stiffness with stress in the || and ⊥ directions. In a VTI medium, the compression of the || component (such as clay and kerogen layers) dominates, so increasing the confining stress will result in DC || <DC ⊥ . In an HTI medium, the deformation of the soft components (such as fractures) arranged perpendicular to the bedding is more significant, so the dominant deformation will be in the || direction, resulting in DC || >DC ⊥ . The present invention evaluates the source of the change in the stress-related stiffness of the shale formation and the anisotropy of the medium by establishing an anisotropic interpretation template for mature shale and analyzing the change in the elastic modulus of the rock in different directions due to stress changes based on mineralogical and acoustic data. By deepening the understanding of the shale rock structure or symmetry, it helps to select an appropriate rock skeleton model and optimize the hydraulic fracturing design.

[0061] In one example, constructing a mature shale rock physics model includes:

[0062] Obtain the equivalent elastic stiffness of the inorganic rock skeleton through anisotropic SCA;

[0063] Obtain the equivalent elastic stiffness of the pore-containing kerogen through anisotropic DEM;

[0064] Based on the equivalent elastic stiffness of the inorganic rock skeleton and the equivalent elastic stiffness of the pore-containing kerogen, add the pore-containing kerogen to the inorganic shale skeleton through anisotropic Backus averaging to construct a mature shale rock physics model.

[0065] In one example, the expression of anisotropic SCA is:

[0066]

[0067] where v r is the volume fraction of the r-th component (r = 1, 2,..., N) of the inorganic part of the shale, C r represents the stiffness tensor of the r-th component, C SCA represents the self-consistent equivalent elastic tensor of the inorganic skeleton, is a fourth-order tensor given by the stress Green's function integrated over the inclusion shape.

[0068] In one example, the expression of anisotropic DEM is:

[0069]

[0070] where C fluid represents the stiffness of the hydrocarbon-filled pores in the kerogen-related pores, C DEM (v) represents the stiffness of the pore-containing kerogen, v is the content of the kerogen-related pores, C DEM (0) represents the stiffness of the solid kerogen (v = 0, the content of the kerogen-related pores is 0).

[0071] In one example, the expression of anisotropic Backus averaging is:

[0072]

[0073]

[0074]

[0075]

[0076]

[0077] where, and The equivalent stiffness C obtained by Backus averaging for the porous kerogen and inorganic shale matrix 11 and C 13 and C 33 and C 44 and C 66 are the stiffnesses of the porous kerogen or inorganic shale matrix, and <.> represents the weighted average of the properties in the brackets according to their volume ratios.

[0078] Specifically, the transversely isotropic equivalent medium described by the anisotropic Backus averaging can be represented by the following five equivalent stiffnesses:

[0079]

[0080]

[0081]

[0082]

[0083]

[0084] The above five equivalent stiffnesses are the equivalent elastic parameters of the rock, such as where V kerogen is the volume fraction of kerogen, is the stiffness of kerogen, is the stiffness of the inorganic rock matrix.

[0085] In an example, the expressions for the stiffness parallel to the bedding and the stiffness perpendicular to the bedding are respectively:

[0086] C || = ρ × V || 2 ;

[0087]

[0088] where C || is the stiffness of the rock parallel to the bedding, C ⊥ is the stiffness perpendicular to the bedding, ρ is the rock density, V || is the rock velocity parallel to the bedding, and V ⊥ is the rock velocity perpendicular to the bedding.

[0089] Specifically, when studying the longitudinal wave velocity of the rock, When studying the shear wave velocity of the rock,

[0090] In an example, in an isotropic medium, DC|| ≈DC ⊥ ;

[0091] In VTI media, DC || <DC ⊥ ;

[0092] In HTI media, DC || >DC ⊥ ;

[0093] wherein, DC || represents the change of elastic stiffness with stress in the direction parallel to the bedding, and DC ⊥ represents the change of elastic stiffness with stress in the direction perpendicular to the bedding.

[0094] Specifically, the present invention completes the anisotropic interpretation and analysis of mature shale according to the following conditions. In isotropic media, the deformation of the rock in the direction parallel to the bedding and the direction perpendicular to the bedding is approximately equal, that is, DC || ≈DC ⊥ , wherein DC || and DC ⊥ respectively represent the change of elastic stiffness with stress in the direction parallel to the bedding and the direction perpendicular to the bedding. In VTI media, the compression of components parallel to the bedding (such as clay and kerogen layers) dominates, so increasing the confining stress will result in DC || <DC ⊥ . In HTI media, the deformation of soft components arranged perpendicular to the bedding (such as fractures) is more significant, so the dominant deformation will be in the direction parallel to the bedding, resulting in DC || >DC ⊥ .

[0095] The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments, but it is not a limitation of the present invention. It should be noted that, without conflict, the embodiments and features in the embodiments of the present invention can be combined with each other.

[0096] Embodiment 1

[0097] This embodiment provides a method for analyzing the anisotropy of mature shale, including:

[0098] Obtaining the equivalent elastic stiffness of the inorganic rock skeleton through anisotropic SCA; obtaining the equivalent elastic stiffness of the pore-containing kerogen through anisotropic DEM; adding the pore-containing kerogen to the inorganic shale skeleton based on the equivalent elastic stiffness of the inorganic rock skeleton and the equivalent elastic stiffness of the pore-containing kerogen through anisotropic Backus averaging to construct a rock physics model of mature shale and obtain the equivalent elastic parameters of the rock;

[0099] The expression of anisotropic SCA is:

[0100] where \(v\) r is the volume fraction of the \(r\)-th component (\(r = 1, 2, \ldots, N\)) of the inorganic part of the shale, \(C\) r represents the stiffness tensor of the \(r\)-th component, and \(C\) SCA represents the self-consistent equivalent elastic tensor of the inorganic skeleton, is a fourth-order tensor given by the stress Green's function integrated over the inclusion shape.

[0101] The expression for the anisotropic DEM is:

[0102] where \(C\) DEM (\(v\)) represents the stiffness of the kerogen with pores, \(v\) is the content of kerogen-related pores, and \(C\) DEM (0) represents the stiffness of the solid kerogen (\(v = 0\), the content of kerogen-related pores is 0).

[0103] The expression for the anisotropic Backus average is:

[0104]

[0105]

[0106]

[0107]

[0108]

[0109] where and are the equivalent stiffnesses obtained by using the Backus average for the kerogen with pores and the inorganic shale skeleton, and \(C\) 11 , \(C\) 13 , \(C\) 33 , \(C\) 44 , and \(C\) 66 are the stiffnesses of the kerogen with pores or the inorganic shale skeleton, and \(\langle.\rangle\) represents the weighted average of the properties in the brackets according to their volume ratios.

[0110] Based on the rock equivalent elastic parameters, an anisotropic interpretation template for mature shale is established; based on the anisotropic interpretation template for mature shale, an intersection analysis of the stiffness parallel to the bedding and the stiffness perpendicular to the bedding of mature shale is carried out; the expressions for the stiffness parallel to the bedding and the stiffness perpendicular to the bedding are: \(C\) || =\(\rho\times V\) || 2 ; where \(C\) ||is the stiffness of the rock parallel bedding, C ⊥ is the stiffness of the vertical bedding, ρ is the rock density, V || is the rock velocity parallel to the bedding, V ⊥ is the rock velocity perpendicular to the bedding. Input the data of the confining pressure elastic parameters of the mature shale measured in the laboratory into the anisotropic interpretation template of the mature shale to complete the anisotropic interpretation analysis of the mature shale. In an isotropic medium, DC || ≈DC ⊥ ; in a VTI medium, DC || <DC ⊥ ; in an HTI medium, DC || >DC ⊥ ; where DC || represents the change of the elastic stiffness with stress in the direction parallel to the bedding, and DC ⊥ represents the change of the elastic stiffness with stress in the direction perpendicular to the bedding.

[0111] Example 2

[0112] As Figure 2 shown, this embodiment provides a method for analyzing the anisotropy of mature shale, including:

[0113] Step 1: Use anisotropic SCA for quartz, calcite rock minerals, and oriented clays to obtain the equivalent elastic stiffness of the inorganic rock skeleton; Anisotropic SCA can be expressed by the following formula: where V r is the volume fraction of the r-th component (r = 1, 2,..., N) of the inorganic part of the shale. C r represents the stiffness tensor of the r-th component, and C SCA represents the self-consistent equivalent elastic tensor of the inorganic skeleton, is a fourth-order tensor given by the stress Green's function integrated over the inclusion shape;

[0114] Step 2: Use anisotropic DEM to obtain the equivalent elastic stiffness of the kerogen with pores; Anisotropic DEM can be expressed by the following formula:

[0115]

[0116] where C fluid represents the stiffness of the hydrocarbon-filled pores in the kerogen-related pores, C DEM (v) represents the stiffness of the kerogen with pores, v is the content of the kerogen-related pores, and C DEM (0) represents the stiffness of the solid kerogen (v = 0, the content of the kerogen-related pores is 0);

[0117] Step 3: Incorporate the porous kerogen into the shale framework using the anisotropic Backus average to construct a mature shale petrophysical model. The transversely isotropic equivalent medium described by the anisotropic Backus average can be represented by the following five equivalent stiffnesses:

[0118] where and are the equivalent stiffnesses calculated using the Backus average for the porous kerogen and the inorganic shale framework. C 11 , C 13 , C 33 , C 44 and C 66 are the stiffnesses of the porous kerogen or the inorganic shale framework. <.> represents the weighted average of the properties in the brackets according to their volume ratios.

[0119] Step 4: Establish an anisotropic interpretation template for mature shale using the rock equivalent elastic parameters obtained from the mature shale petrophysical model. Mainly, perform crossplot analysis on the stiffness parallel to the bedding (C || ) and the stiffness perpendicular to the bedding (C ⊥ ) of the shale. The stiffness parallel to the bedding and the stiffness perpendicular to the bedding of the rock can be expressed as: C || = ρ × V || 2 , where ρ is the density, and V || and V ⊥ are the rock velocities parallel and perpendicular to the bedding, respectively.

[0120] Step 5: Project the laboratory-measured data of the elastic parameters of mature shale under varying confining pressures onto the anisotropic interpretation template to complete the anisotropic interpretation analysis of mature shale. In an isotropic medium, the deformations of the rock in the directions parallel (||) and perpendicular (⊥) to the bedding are approximately equal, i.e., DC || ≈ DC ⊥ , where DC || and DC ⊥ represent the changes in elastic stiffness with stress in the || and ⊥ directions, respectively. In a VTI medium, the compression of the || components (such as clay and kerogen layers) dominates, so increasing the confining stress will result in DC || < DC ⊥ . In an HTI medium, the deformation of the soft components arranged perpendicular to the bedding (such as fractures) is more significant, so the dominant deformation will be in the || direction, resulting in DC || > DC ⊥ .

[0121] The anisotropic interpretation template for mature shale was applied to the elastic data of continental mature shale measured in the laboratory in the Sichuan Basin. The mature shale in the study area generally shows the characteristics of VTI, and some shale samples are isotropic.

[0122] Figure 3 This is the longitudinal wave anisotropic interpretation template for mature shale. This interpretation template can analyze the anisotropic characteristics of shale under different clay contents (1 - QF) and organic matter contents (OM); it can be seen that as the confining pressure increases, the longitudinal wave anisotropic parameters of most samples gradually decrease; this may be because the increase in pressure compresses the organic matter, making the rock show an effect similar to reducing the volume content of organic matter, thus weakening the anisotropy; for most samples, as the confining pressure increases, the change in elastic stiffness in the direction approximately along the isolith line (same quartz / feldspar content QF) is less than the change in elastic stiffness in the ⊥ direction, that is, DC || <DC ⊥ , the mature shale in the study area is generally a VTI medium as a whole; for some samples, as the confining pressure increases, the change in elastic stiffness in the || direction is equal to the change in elastic stiffness in the ⊥ direction, that is, DC || <=DC ⊥ , showing the characteristics of an isotropic medium. In addition, it can be seen from the template that the organic matter content of these samples is relatively low, which may be the main reason for their weak anisotropy.

[0123] Figure 4 This is the transverse wave anisotropic interpretation template for mature shale. The test samples show similar characteristics on the transverse wave anisotropic interpretation template to those on the longitudinal wave anisotropic interpretation template; it can also be known that the mature shale in the study area is generally a VTI medium as a whole; it should be noted that the elastic moduli of some samples do not appear at the same position on the transverse wave anisotropic interpretation template as on the longitudinal wave anisotropic interpretation template, which may be due to the fact that the template does not consider fractures or does not fully consider the influence of fluids.

[0124] Example 3

[0125] As Figure 5 shown, this embodiment provides an anisotropic analysis device for mature shale, including:

[0126] A construction module for constructing a petrophysical model of mature shale and obtaining rock equivalent elastic parameters;

[0127] An establishment module for establishing an anisotropic interpretation template for mature shale based on the rock equivalent elastic parameters;

[0128] An intersection analysis module for performing intersection analysis on the parallel bedding stiffness and perpendicular bedding stiffness of mature shale based on the anisotropic interpretation template for mature shale;

[0129] An interpretation and analysis module, configured to input the data of the confining pressure elastic parameters of mature shale measured in the laboratory into the anisotropic interpretation template of mature shale, and complete the anisotropic interpretation and analysis of mature shale.

[0130] Embodiment 4

[0131] This embodiment provides an electronic device, which includes:

[0132] At least one processor; and,

[0133] A memory communicatively connected to the at least one processor; wherein,

[0134] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the mature shale anisotropy analysis method in Embodiment 1.

[0135] The electronic device according to an embodiment of the present disclosure includes a memory and a processor, and the memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0136] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In an embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0137] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience effect, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included in the protection scope of the present disclosure.

[0138] For the detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, and details are not repeated herein.

[0139] Embodiment 5

[0140] This embodiment provides a non-transitory computer-readable storage medium, which stores computer instructions for causing a computer to execute the mature shale anisotropy analysis method in Embodiment 1.

[0141] A computer-readable storage medium according to an embodiment of the present disclosure stores non-transitory computer-readable instructions. When the non-transitory computer-readable instructions are run by a processor, all or part of the steps of the methods of the various embodiments of the present disclosure described above are executed.

[0142] The above computer-readable storage medium includes, but is not limited to: optical storage media (e.g., CD-ROMs and DVDs), magneto-optical storage media (e.g., MOs), magnetic storage media (e.g., magnetic tapes or removable hard disks), media with built-in rewritable non-volatile memories (e.g., memory cards), and media with built-in ROMs (e.g., ROM cartridges).

[0143] The various embodiments of the present invention have been described above. The above description is exemplary and not exhaustive, and is also not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for analyzing the anisotropy of mature shale, characterized in that, it includes: Construct a petrophysical model of mature shale to obtain the equivalent elastic parameters of the rock; Based on the equivalent elastic parameters of the rock, establish an anisotropy interpretation template for mature shale; Based on the anisotropy interpretation template of mature shale, conduct cross-analysis on the parallel bedding stiffness and perpendicular bedding stiffness of mature shale; Input the data of the variable confining pressure elastic parameters of mature shale measured in the laboratory into the anisotropy interpretation template of mature shale to complete the anisotropy interpretation analysis of mature shale.

2. The method for analyzing the anisotropy of mature shale according to claim 1, characterized in that, The construction of the petrophysical model of mature shale includes: Obtain the equivalent elastic stiffness of the inorganic rock skeleton through anisotropic SCA; Obtain the equivalent elastic stiffness of the pore-containing kerogen through anisotropic DEM; Based on the equivalent elastic stiffness of the inorganic rock skeleton and the equivalent elastic stiffness of the pore-containing kerogen, add the pore-containing kerogen to the inorganic shale skeleton through anisotropic Backus averaging to construct the petrophysical model of mature shale.

3. The method for analyzing the anisotropy of mature shale according to claim 2, characterized in that, The expression of the anisotropic SCA is: where, v r is the volume fraction of the r-th component (r = 1, 2, …, N) of the inorganic part of the shale, C r denotes the stiffness tensor of the r-th component, C SCC denotes the self-consistent equivalent elastic tensor of the inorganic skeleton, is a fourth-order tensor given by the stress Green's function integrated over the inclusion shape.

4. The method for analyzing the anisotropy of mature shale according to claim 2, characterized in that, The expression of the anisotropic DEM is: Among them, C fluid represents the stiffness of the hydrocarbon-filled pores in the kerogen-related pores, C DEM (v) represents the stiffness of the kerogen with pores, where v is the content of the kerogen-related pores, C DEM (0) represents the stiffness of the solid kerogen (v = 0, the content of the kerogen-related pores is 0).

5. The method for analyzing the anisotropy of mature shale according to claim 2, characterized in that, The expression of the anisotropic Backus averaging is: Among them, and are the equivalent stiffnesses calculated using the Backus average for the porous kerogen and the inorganic shale matrix. C 11 , C 13 , C 33 , C 44 and C 66 are the stiffnesses of the porous kerogen or the inorganic shale matrix. <.> represents the weighted average of the properties in the parentheses according to their volume ratios.

6. The method for analyzing the anisotropy of mature shale according to claim 1, characterized in that, The expressions of the parallel bedding stiffness and the perpendicular bedding stiffness are respectively: Among them, C || is the stiffness of the rock parallel to the bedding, C ⊥ is the stiffness perpendicular to the bedding, ρ is the rock density, V || is the rock velocity parallel to the bedding, V ⊥ is the rock velocity perpendicular to the bedding.

7. The method for analyzing the anisotropy of mature shale according to claim 6, characterized in that, In an isotropic medium, DC || ≈ DC ⊥ ; In VTI media, DC || <DC ⊥ ; In the HTI medium, DC || >DC ⊥ ; Among them, DC || represents the change of elastic stiffness with stress in the direction parallel to the bedding plane, and DC ⊥ represents the change of elastic stiffness with stress in the direction perpendicular to the bedding plane.

8. An electronic device, characterized in that, The electronic device includes: At least one processor; and, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the method for analyzing the anisotropy of mature shale according to any one of claims 1-7.

9. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute the method for analyzing the anisotropy of mature shale according to any one of claims 1-7.

10. An apparatus for analyzing the anisotropy of mature shale, characterized in that, it includes: A construction module for constructing a petrophysical model of mature shale to obtain the equivalent elastic parameters of the rock; An establishment module for establishing an anisotropy interpretation template for mature shale based on the equivalent elastic parameters of the rock; A cross-analysis module for conducting cross-analysis on the parallel bedding stiffness and perpendicular bedding stiffness of mature shale based on the anisotropy interpretation template of mature shale; An interpretation and analysis module, configured to input the data of the confining pressure-dependent elastic parameters of mature shale measured in the laboratory into the anisotropic interpretation template of the mature shale, and complete the anisotropic interpretation and analysis of the mature shale.

Images