Terrestrial shale petrophysical modeling method, device, equipment and medium

Abstract

The present disclosure provides a method, device, equipment and medium for modeling the lithology of continental shale, and relates to the field of reservoir lithology in geological exploration. The method comprises the following steps: dividing the minerals into a laminated part equivalent mineral proportion and a non-laminated part equivalent mineral proportion; obtaining the laminated part equivalent mineral proportion according to the laminated density, the laminated thickness and the unit thickness, wherein the laminated thickness of each laminated layer is the same; adding the laminated part equivalent mineral proportion to the lithology modeling by using a Backus average model; and adding the non-laminated part equivalent mineral proportion to the lithology modeling. According to the present disclosure, the model has better applicability and higher prediction accuracy.

Description

Technical Field

[0001] This disclosure relates to the field of reservoir rock physics in geological exploration, and in particular to a method, apparatus, equipment and medium for rock physics modeling of continental shale. Background Technology

[0002] Shale gas has large recoverable reserves and is an important energy resource. However, the unique geological conditions of shale reservoirs result in complex mineral composition and pore fluid characteristics. Horizontally arranged clay minerals and laminae give shale strong anisotropy, and the coexistence of multiple pore types and organic kerogen demonstrate its heterogeneity, making gas layer identification difficult and the elastic characteristics of the reservoir unclear. Seismic rock physics models are an important technique for guiding the macroscopic elastic analysis of shale reservoirs. With the increasing demands of seismic exploration, more accurate shale reservoir physical property information requires the development of shale rock physics models that better reflect actual underground conditions.

[0003] It should be noted that the information disclosed in the background section above is only used to enhance the understanding of the background of this disclosure, and therefore may include information that does not constitute prior art known to those skilled in the art. Summary of the Invention

[0004] This disclosure provides a method, apparatus, equipment, and medium for rock physics modeling of continental shale, which at least partially solves the problem of insufficient prediction accuracy of conventional rock physics models.

[0005] Other features and advantages of this disclosure will become apparent from the following detailed description, or may be learned in part from practice of this disclosure.

[0006] According to one aspect of this disclosure, a method for petrological modeling of continental shale is provided, comprising:

[0007] The minerals are divided into the equivalent mineral percentage of the lamellar portion and the equivalent mineral percentage of the non-laminated portion.

[0008] Based on the lamellar density, lamellar thickness, and unit thickness, the equivalent mineral proportion of the lamellar portion is obtained, wherein the lamellar thickness of each lamellar is the same;

[0009] The Backus average model was used to incorporate the equivalent mineral proportion of the lamellar portion into the rock physics modeling;

[0010] The equivalent mineral percentage of the non-laminated portion is added to the rock physics model.

[0011] In one embodiment of this disclosure, the method further includes:

[0012] The SCA theory is applied to incorporate mineral matrix and intergranular pores into rock physics modeling.

[0013] In one embodiment of this disclosure, the method further includes:

[0014] Add kerogen to rock physics modeling.

[0015] In one embodiment of this disclosure, the equivalent mineral percentage of the laminar portion is obtained based on laminar density, laminar thickness, and unit thickness, including:

[0016] The total number of laminae is obtained based on the laminae density and unit thickness;

[0017] The equivalent mineral percentage of the lamellar portion is obtained based on the total number of lamellar lines, lamellar thickness, and unit thickness.

[0018] In one embodiment of this disclosure, the method further includes:

[0019] Based on the total mineral content and the equivalent mineral percentage of the lamellar portion, the equivalent mineral percentage of the non-laminated portion is obtained.

[0020] In one embodiment of this disclosure, the method further includes:

[0021] Using the Eshelby-Cheng model of anisotropic media, fluid-containing ellipsoidal fractures along the three-axis perpendicular direction were added to isotropic rocks.

[0022] According to another aspect of this disclosure, a petrological modeling apparatus for continental shale is provided, comprising:

[0023] The mineral classification module is used to divide minerals into equivalent mineral percentages in lamellar regions and equivalent mineral percentages in non-laminated regions.

[0024] The calculation module is used to obtain the equivalent mineral proportion of the laminar portion based on the laminar density, laminar thickness, and unit thickness, wherein the laminar thickness of each laminar portion is the same;

[0025] The first modeling module is used to incorporate the equivalent mineral proportion of the laminar part into the rock physics model using the Backus average model;

[0026] The second modeling module is used to incorporate the equivalent mineral proportions of non-laminated parts into the rock physics modeling.

[0027] According to another aspect of this disclosure, an electronic device is provided, comprising: a memory for storing instructions; and a processor for calling the instructions stored in the memory to implement the above-described method for petrological modeling of continental shale.

[0028] According to another aspect of this disclosure, a computer-readable storage medium is provided having computer instructions stored thereon, which, when executed by a processor, implement the above-described method for petrophysical modeling of continental shale.

[0029] According to another aspect of this disclosure, a computer program product is provided, which stores instructions that, when executed by a computer, cause the computer to perform the above-described method for modeling continental shale rock physics.

[0030] According to another aspect of this disclosure, a chip is provided, including at least one processor and an interface;

[0031] An interface is used to provide program instructions or data to at least one processor;

[0032] At least one processor is used to execute program instructions to implement the above-described petrological modeling method for continental shale.

[0033] The petrophysical modeling method for continental shale provided in this embodiment takes into account the characteristics of continental shale such as complex lithology, strong anisotropy, and well-developed laminae. It qualitatively considers the laminae, and since the laminae thickness is the same, the overall influence of all laminae on the model can be obtained based on the laminae density. The Backus model is used to incorporate the laminae into the petrophysical model, which has better applicability and higher prediction accuracy.

[0034] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and are not intended to limit this disclosure. Attached Figure Description

[0035] The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments consistent with this disclosure and, together with the description, serve to explain the principles of this disclosure.

[0036] Obviously, the accompanying drawings described below are merely some embodiments of this disclosure. Those skilled in the art can obtain other drawings based on these drawings without any creative effort.

[0037] Figure 1 This diagram illustrates a flowchart of a petrophysical modeling method for terrestrial shale according to an embodiment of the present disclosure;

[0038] Figure 2 This diagram illustrates a flowchart of another method for rock physics modeling of continental shale in an embodiment of this disclosure;

[0039] Figure 3 This diagram illustrates the influence of lamination on a rock physical model in an embodiment of this disclosure.

[0040] Figure 4 This illustrates the P-wave prediction results for well A in an embodiment of this disclosure;

[0041] Figure 5 shows a comparison between Young's modulus and Poisson's ratio of the well A model in the embodiment of this disclosure and the measured values;

[0042] Figure 6 This diagram illustrates a physical modeling apparatus for terrestrial shale in an embodiment of the present disclosure.

[0043] Figure 7 A structural block diagram of an electronic device according to an embodiment of the present disclosure is shown. Detailed Implementation

[0044] The exemplary implementation will now be described more fully with reference to the accompanying drawings.

[0045] It should be noted that the example implementation can be implemented in many forms and should not be construed as being limited to the examples set forth herein.

[0046] Earthquake rock physics theory shows that the microscopic properties of underground rocks, such as lithology, mineral composition, cementation, pressure, temperature, pore fluid properties (formation water salinity, crude oil density, gas-oil ratio, etc.), and fluid distribution within pores (uniform or patchy distribution, etc.), are important factors affecting the macroscopic elastic characteristics of rocks.

[0047] The inventors discovered that shale reservoirs are classified into continental shale, marine shale, and transitional marine-continental shale, with continental shale also known as lacustrine shale. Continental shale oil and gas reservoirs are characterized by abundant clastic particles, poor sorting, weak compaction, numerous narrow intergranular pores, and well-developed horizontal bedding. Another prominent feature of continental shale is the distinctive laminar structure formed by interbedded argillaceous and silty materials, typically 1 mm or a few millimeters thick. This laminar structure significantly influences anisotropy parameters. However, conventional shale petrophysical modeling does not consider the laminar structure's influence. Incorporating laminar density into the petrophysical model can improve prediction accuracy. Therefore, this study focuses on continental shale, investigating the impact of laminar structure on reservoir seismic petrophysical analysis, and proposing an effective modeling process. This research can provide strong support and guidance for subsequent seismic inversion and interpretation.

[0048] The following detailed description of this exemplary implementation method is provided in conjunction with the accompanying drawings and embodiments.

[0049] Figure 1 This diagram illustrates a flowchart of a petrophysical modeling method for continental shale according to an embodiment of the present disclosure, as follows: Figure 1 As shown, the terrestrial shale rock physics modeling method provided in this embodiment includes steps S110-S140.

[0050] In S110, minerals are divided into equivalent mineral proportions in the lamellar portion and equivalent mineral proportions in the non-laminated portion.

[0051] In S120, the equivalent mineral percentage of the lamellar portion is obtained based on the lamellar density, lamellar thickness, and unit thickness, wherein the lamellar thickness of each lamellar portion is the same.

[0052] In S130, the Backus average model is used to incorporate the equivalent mineral proportion of the lamellar portion into the rock physics modeling.

[0053] In S140, the equivalent mineral percentage of the non-laminated portion is added to the rock physics model.

[0054] In this embodiment, minerals are divided into lamellar and non-laminated portions with equivalent mineral proportions. The Backus average model and the SCA model are respectively incorporated into the rock physics modeling to consider the influence of lamellarity on the rock physics modeling.

[0055] In some embodiments, this disclosure also applies SCA theory to incorporate mineral matrix and intergranular pores into rock physical modeling.

[0056] In some embodiments, this disclosure also incorporates kerogen into rock physical modeling.

[0057] In some embodiments, the equivalent mineral percentage of the laminar portion is obtained based on laminar density, laminar thickness, and unit thickness, including: obtaining the total number of laminar lines based on laminar density and unit thickness; and obtaining the equivalent mineral percentage of the laminar portion based on the total number of laminar lines, laminar thickness, and unit thickness.

[0058] In some embodiments, this disclosure also obtains the equivalent mineral percentage of the non-laminated portion based on the total mineral content and the equivalent mineral percentage of the lamination portion.

[0059] In some embodiments, this disclosure also utilizes the Eshelby-Cheng model of anisotropic media to add fluid-containing ellipsoidal fractures along the three-axis perpendicular direction to isotropic rocks.

[0060] Figure 2 This diagram illustrates a flowchart of a petrophysical modeling method for continental shale according to an embodiment of the present disclosure, as follows: Figure 2 As shown, in this embodiment, the mineral matrix and intergranular pores are incorporated using the SCA theory. Berryman (1980, 1995) gave the general form of the self-compatible approximation for N-phase mixtures as follows:

[0061]

[0062] in,

[0063] Adding kerogen: The self-composite modulus of the two-phase mixture is estimated using the following equation:

[0064]

[0065] Among them, K mixture ,μ mixtureThese are the bulk modulus and shear modulus of the mixed matrix after the addition of organic matter, K. kerogen ,μ kerogen These are the bulk modulus and shear modulus of organic kerogen, V. kerogen P represents the volume content of kerogen. *k Q *k It is a geometric factor related to the shape of the added kerogen.

[0066] like Figure 3 Laminar density ρ w This refers to the number of laminae per unit thickness. Minerals are divided into laminate and non-laminated portions with equivalent mineral proportions. The Backus average model and the SCA model are respectively incorporated into the rock physics modeling to consider the impact of laminae on rock physics modeling. Assuming the laminate thickness is ΔL, the number of laminae in the model is calculated, thus yielding the equivalent mineral proportion of the laminate portion.

[0067] D = ρ w *H (3)

[0068] X=D*ΔL / H (4)

[0069] Y = TX (5)

[0070] Where D is the total number of laminae, ρ w H is the laminar density, H is the unit thickness, X is the equivalent mineral percentage of the laminar portion, T is the total mineral content, and Y is the equivalent mineral percentage of the non-laminated portion.

[0071] The laminar content is averaged using Backus's method and incorporated into the model:

[0072]

[0073]

[0074]

[0075]

[0076]

[0077] Where <·> represents the average value of the closed property weighted by volume.

[0078] Using the Eshelby-Cheng model of anisotropic media, fluid-containing ellipsoidal fractures along the three-axis perpendicular direction were added to isotropic rocks.

[0079]

[0080] The modulus and density of mixtures of different types of fluids are calculated using the Wood formula (Wood, 1955). The formula is as follows:

[0081]

[0082]

[0083] Among them, K R f is the Reuss average of the mixture, and ρ is the average density. i K i and ρ i These are the volume content, bulk modulus, and density of each component.

[0084] Calculate the longitudinal wave velocity and physical parameters of saturated fractured rocks.

[0085]

[0086] In this embodiment, considering the complex lithology, strong anisotropy, and well-developed laminae of continental shale, we qualitatively consider the laminae, assuming that the laminae thickness is the same and ΔL. Based on the laminae density, we can obtain the overall influence of all laminae on the model. The Backus model is used to incorporate the laminae into the rock physics model. This model was verified in the Jurassic continental shale of the Sichuan Basin, and the P-wave velocity, S-wave modulus, Poisson's ratio, and measured data showed good agreement, verifying the applicability of the model. Figure 4 The P-wave prediction results for well A are shown, where the dashed line represents the prediction results and the solid line represents the measured results. Figure 5a and Figure 5b The graph shows a comparison between Young's modulus and Poisson's ratio of the A-well model and the measured values. The dashed line represents the predicted results, and the solid line represents the measured results.

[0087] Furthermore, although the steps of the method in this disclosure are described in a specific order in the accompanying drawings, this does not require or imply that the steps must be performed in that specific order, or that all the steps shown must be performed to achieve the desired result.

[0088] In some embodiments, certain steps may be omitted, multiple steps may be combined into one step for execution, and / or one step may be broken down into multiple steps for execution.

[0089] Based on the same inventive concept, this disclosure also provides a porosity prediction device, as described in the following embodiments. Since the principle by which this device solves the problem is similar to that of the method embodiments described above, the implementation of this device embodiment can refer to the implementation of the method embodiments described above, and repeated details will not be elaborated further.

[0090] Figure 6This diagram illustrates a petrological modeling apparatus for terrestrial shale, as shown in an embodiment of the present disclosure. Figure 6 As shown, the terrestrial shale rock physics modeling device 600 includes:

[0091] The mineral classification module 602 is used to classify minerals into equivalent mineral proportions of lamellar portions and equivalent mineral proportions of non-laminated portions.

[0092] The calculation module 604 is used to obtain the equivalent mineral proportion of the laminar portion based on the laminar density, laminar thickness and unit thickness, wherein the laminar thickness of each laminar portion is the same;

[0093] The first modeling module 606 is used to incorporate the equivalent mineral proportion of the lamellar portion into the rock physics model using the Backus average model;

[0094] The second modeling module 608 is used to add the equivalent mineral proportion of the non-laminated part to the rock physics modeling.

[0095] In some embodiments, the continental shale rock physics modeling apparatus 600 further includes a third modeling module. The third modeling module is used to apply SCA theory to incorporate mineral matrix and intergranular pores into the rock physics modeling.

[0096] In some embodiments, the third modeling module is also used to incorporate kerogen into rock physical modeling.

[0097] In some embodiments, the calculation module 604 is used to obtain the total number of laminae based on the laminae density and unit thickness; and to obtain the equivalent mineral percentage of the laminae based on the total number of laminae, laminae thickness, and unit thickness.

[0098] In some embodiments, the continental shale petrophysical modeling apparatus 600 further includes a fourth modeling module. The fourth modeling module is used to obtain the equivalent mineral percentage of the non-laminated portion based on the total mineral content and the equivalent mineral percentage of the laminated portion.

[0099] In some embodiments, the fourth modeling module is further configured to add fluid-containing ellipsoidal fractures along the three-axis perpendicular direction to isotropic rocks using the Eshelby-Cheng model of anisotropic media.

[0100] In some embodiments, the fourth modeling module is also used to calculate the modulus and density of different types of fluids after mixing using formulas (12) and (13).

[0101] The concepts of "first" and "second" mentioned in this disclosure are used only to distinguish different devices, modules or units, and are not used to define the order of functions performed by these devices, modules or units or their interdependencies.

[0102] Regarding the continental shale rock physics modeling device in the above embodiments, the specific methods by which each module performs its operation have been described in detail in the embodiments related to the continental shale rock physics modeling method, and will not be elaborated here.

[0103] It should be noted that although several modules or units of the device used for action execution are mentioned in the detailed description above, this division is not mandatory.

[0104] In fact, according to embodiments of this disclosure, the features and functions of two or more modules or units described above can be embodied in one module or unit. Conversely, the features and functions of one module or unit described above can be further divided and embodied by multiple modules or units.

[0105] Some of the block diagrams shown in the accompanying drawings are functional entities and do not necessarily correspond to physically or logically independent entities. These functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0106] The following reference Figure 7 This describes the electronic device provided in the embodiments of this disclosure. Figure 7 The electronic device 700 shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments disclosed herein.

[0107] Figure 7 This diagram illustrates the architecture of an electronic device 700 provided in an embodiment of the present invention. Figure 7 As shown, the electronic device 700 includes, but is not limited to, at least one processor 710 and at least one memory 720.

[0108] Memory 720 is used to store instructions.

[0109] In some embodiments, memory 720 may include a readable medium in the form of volatile memory cells, such as random access memory (RAM) 7201 and / or cache memory 7202, and may further include read-only memory (ROM) 7203.

[0110] In some embodiments, the memory 720 may also include a program / utility 7204 having a set (at least one) program module 7205, such program module 7205 including but not limited to: an operating system, one or more application programs, other program modules, and program data, each or some combination of these examples may include an implementation of a network environment.

[0111] In some embodiments, the memory 720 may store an operating system. This operating system may be a real-time operating system (RTX), such as Linux, UNIX, Windows, or OS X.

[0112] In some embodiments, the memory 720 may also store data.

[0113] As an example, processor 710 can read data stored in memory 720, which may be stored at the same memory address as the instruction, or the data may be stored at a different memory address than the instruction.

[0114] Processor 710 is configured to invoke instructions stored in memory 720 to implement the steps described in the "Exemplary Methods" section above, according to various exemplary embodiments of this disclosure. For example, processor 710 may execute the steps of the above method embodiments.

[0115] It should be noted that the processor 710 described above can be a general-purpose processor or a special-purpose processor. The processor 710 may include one or more processing cores, and the processor 710 executes various functional applications and data processing by running instructions.

[0116] In some embodiments, processor 710 may include a central processing unit (CPU) and / or a baseband processor.

[0117] In some embodiments, the processor 710 may determine an instruction based on the priority identifier and / or function category information carried in each control instruction.

[0118] In this disclosure, the processor 710 and the memory 720 can be configured separately or integrated together.

[0119] As an example, the processor 710 and memory 720 can be integrated on a single board or a system-on-a-chip (SOC).

[0120] like Figure 7 As shown, the electronic device 700 is embodied in the form of a general-purpose computing device. The electronic device 700 may also include a bus 730.

[0121] Bus 730 can represent one or more of several types of bus structures, including a memory bus or memory controller, peripheral bus, graphics acceleration port, processor, or a local bus using any of the various bus structures.

[0122] Electronic device 700 can also communicate with one or more external devices 740 (e.g., keyboard, pointing device, Bluetooth device, etc.), and with one or more devices that enable a user to interact with electronic device 700, and / or with any device that enables electronic device 700 to communicate with one or more other computing devices (e.g., router, modem, etc.). Such communication can be performed through input / output (I / O) interface 750.

[0123] Furthermore, the electronic device 700 can also communicate with one or more networks (such as local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) via the network adapter 760.

[0124] like Figure 7 As shown, the network adapter 760 communicates with other modules of the electronic device 700 via the bus 730.

[0125] It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in conjunction with the electronic device 700, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0126] It is understood that the structure illustrated in the embodiments of this disclosure does not constitute a specific limitation on the electronic device 700. In other embodiments of this disclosure, the electronic device 700 may include more than Figure 7 This may involve more or fewer components, or combining certain components, or splitting certain components, or different component arrangements. Figure 7 The components shown can be implemented in hardware, software, or a combination of both.

[0127] This disclosure also provides a computer-readable storage medium storing computer instructions thereon, which, when executed by a processor, implement the terrestrial shale rock physics modeling method described in the above method embodiments.

[0128] In this embodiment of the disclosure, the computer-readable storage medium is a computer instruction that can be sent, propagated, or transmitted for use by or in conjunction with an instruction execution system, apparatus, or device.

[0129] As an example, a computer-readable storage medium is a non-volatile storage medium.

[0130] In some embodiments, more specific examples of computer-readable storage media in this disclosure may include, but are not limited to: electrical connections having one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, USB flash drives, portable hard drives, or any suitable combination of the foregoing.

[0131] In this embodiment of the disclosure, the computer-readable storage medium may include data signals propagated in baseband or as part of a carrier wave, wherein computer instructions (readable program code) are carried.

[0132] The transmitted data signal can take many forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof.

[0133] In some examples, computational instructions contained on a computer-readable storage medium may be transmitted using any suitable medium, including but not limited to wireless, wired, optical fiber, RF, etc., or any suitable combination thereof.

[0134] This disclosure also provides a computer program product that stores instructions that, when executed by a computer, cause the computer to implement the terrestrial shale rock physics modeling method described in the above method embodiments.

[0135] The aforementioned instructions can be program code. In practice, the program code can be written using any combination of one or more programming languages.

[0136] Programming languages ​​include object-oriented programming languages—such as Java and C++—as well as conventional procedural programming languages—such as the "C" language or similar programming languages.

[0137] The program code can be executed entirely on the user's computing device, partially on the user's computing device, as a standalone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0138] In cases involving remote computing devices, the remote computing devices can be connected to user computing devices via any type of network, including local area networks (LANs) or wide area networks (WANs), or they can be connected to external computing devices (e.g., via the Internet using an Internet service provider).

[0139] This disclosure also provides a chip, including at least one processor and an interface;

[0140] An interface is used to provide program instructions or data to at least one processor;

[0141] At least one processor is used to execute program instructions to implement the terrestrial shale rock physics modeling method described in the above method embodiments.

[0142] In some embodiments, the chip may further include a memory for storing program instructions and data, the memory being located within or outside the processor.

[0143] Those skilled in the art will understand that all or part of the steps of the above embodiments can be specifically implemented in the following forms: a completely hardware implementation, a completely software implementation (including firmware, microcode, etc.), or a combination of hardware and software implementations, which can be collectively referred to as "circuit", "module" or "system".

[0144] Other embodiments of this disclosure will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein.

[0145] This disclosure is intended to cover any variations, uses, or adaptations of this disclosure that follow the general principles of this disclosure and include common knowledge or customary techniques in the art not disclosed herein. The description and examples are to be considered exemplary only, and the true scope and spirit of this disclosure are indicated by the appended claims.

Claims

1. A method for rock physical modeling of continental shale, characterized in that, include: The minerals are divided into the equivalent mineral percentage of the lamellar portion and the equivalent mineral percentage of the non-laminated portion. Based on the lamellar density, lamellar thickness, and unit thickness, the equivalent mineral proportion of the lamellar portion is obtained, wherein the lamellar thickness of each lamellar is the same; The equivalent mineral proportion of the lamellar portion was incorporated into the rock physics model using the Backus average model. The equivalent mineral proportion of the non-laminated portion is added to the rock physics modeling; Based on laminar density, laminar thickness, and unit thickness, the equivalent mineral percentage of the laminar portion is obtained, including: The total number of laminae is obtained based on the laminae density and unit thickness; Based on the total number of laminae, laminae thickness, and unit thickness, the equivalent mineral percentage of the laminae portion is obtained; Based on the total mineral content and the equivalent mineral percentage of the lamellar portion, the equivalent mineral percentage of the non-laminated portion is obtained.

2. The method according to claim 1, characterized in that, The method further includes: The SCA theory is applied to incorporate mineral matrix and intergranular pores into rock physics modeling.

3. The method according to claim 1, characterized in that, The method further includes: Add kerogen to rock physics modeling.

4. The method according to claim 1, characterized in that, The method further includes: Using the Eshelby-Cheng model of anisotropic media, fluid-containing ellipsoidal fractures along the three axes perpendicular to the rock were added to isotropic rocks.

5. The method according to claim 1, characterized in that, The method further includes: The modulus and density of a mixture of different types of fluids can be calculated using the following formula: ； ； in, It is the Reuss average of the mixture. It is the average density; This indicates the volume content of each component. This represents the bulk modulus of each component. This indicates the density of each component.

6. A physical modeling device for terrestrial shale, characterized in that, include: The mineral classification module is used to divide minerals into equivalent mineral percentages in lamellar regions and equivalent mineral percentages in non-laminated regions. The calculation module is used to obtain the equivalent mineral proportion of the laminar portion based on the laminar density, laminar thickness, and unit thickness, wherein the laminar thickness of each laminar portion is the same; The first modeling module is used to incorporate the equivalent mineral proportion of the lamellar portion into the rock physics modeling using the Backus average model; The second modeling module is used to add the equivalent mineral proportion of the non-laminated part to the rock physical modeling; Based on laminar density, laminar thickness, and unit thickness, the equivalent mineral percentage of the laminar portion is obtained, including: The total number of laminae is obtained based on the laminae density and unit thickness; Based on the total number of laminae, laminae thickness, and unit thickness, the equivalent mineral percentage of the laminae portion is obtained; Based on the total mineral content and the equivalent mineral percentage of the lamellar portion, the equivalent mineral percentage of the non-laminated portion is obtained.

7. An electronic device, characterized in that, include: Memory, used to store instructions; A processor is configured to invoke instructions stored in the memory to implement the petrophysical modeling method for continental shale as described in any one of claims 1-5.

8. A computer-readable storage medium storing computer instructions thereon, characterized in that, When the computer instructions are executed by the processor, they implement the petrological modeling method for terrestrial shale as described in any one of claims 1-7.

Images