Construction method, device and equipment of rock physical model, medium and product

Abstract

The invention discloses a rock physical model construction method, device and equipment, a medium and a product. The method comprises the following steps: determining a compaction factor according to porosity data and stratum thickness data; determining the elasticity modulus of the dry rock skeleton according to the mineral component content and the compaction factor; determining a fluid bulk modulus and a fluid density according to the water saturation data and the gas saturation data; determining a saturated rock elastic modulus according to the dry rock skeleton elastic modulus and the fluid volume modulus; determining the saturated rock density according to the fluid density and the rock matrix equivalent density, and further determining the longitudinal wave velocity and the transverse wave velocity of the saturated rock; and determining a rock physical model according to the longitudinal wave velocity and the transverse wave velocity. According to the technical scheme, the problem of poor reliability of rock physical modeling is solved, the compaction factor is introduced in the rock matrix elastic modulus determination process, the rock matrix elastic modulus considering the mechanical compaction effect is calculated, and the accuracy and effectiveness of the rock physical model can be improved.

Description

Technical Field

[0001] The present invention relates to the field of geophysical technology, and in particular to a method, device, equipment, medium and product for constructing a rock physics model. Background Art

[0002] Currently, the demand for oil and gas resources for production and daily life is gradually increasing. With the continuous progress of oil and gas exploration, the output of conventional oil fields is clearly insufficient to meet the demand. Therefore, the focus of oil and gas reservoir development is gradually shifting to the development of unconventional oil and gas. Among them, low-porosity and low-permeability gas reservoirs are a typical unconventional reservoir with abundant reserves and good exploration prospects.

[0003] Low-porosity, low-permeability reservoirs have complex underground environments and are challenging to explore. Reservoir characteristics, and thus reservoir location, are typically determined by analyzing the relationship between reservoir physical properties and elasticity. Seismic data inversion can be used to obtain rock elastic parameters, namely compressional wave velocity, shear wave velocity, and density. Elastic parameters such as bulk modulus, shear modulus, and Poisson's ratio can be calculated using these elastic parameters. Well logging data can also provide a wealth of physical properties of reservoir rocks and fluids. Seismic and well logging data are not matched in terms of volume, necessitating a rock physics model as an effective bridge connecting reservoir and seismic characteristics. This model takes on the task of connecting seismic and well logging data to explain the relationship between elastic parameters and physical properties.

[0004] However, existing technologies for constructing rock physics models typically only consider the effect of mineral composition on elastic modulus, that is, only the effect of chemical compaction. In reality, compaction, a major factor contributing to low porosity and low permeability in reservoirs, includes not only chemical compaction but also mechanical compaction. Therefore, how to incorporate the effects of mechanical compaction into rock physics modeling is an urgent issue. Summary of the Invention

[0005] The present invention provides a method, device, equipment, medium and product for constructing a rock physics model to address the problem of poor reliability of rock physics modeling. By introducing a compaction factor in the process of determining the rock matrix elastic modulus and calculating the rock matrix elastic modulus taking into account the mechanical compaction effect, the accuracy and effectiveness of the rock physics model can be improved.

[0006] According to one aspect of the present invention, a method for constructing a rock physics model is provided, the method comprising:

[0007] Determine the compaction factor based on the porosity data and formation thickness data of the target layer;

[0008] Determine the elastic modulus of the rock matrix based on the content of each mineral component in the target layer and the compaction factor, and determine the elastic modulus of the dry rock skeleton based on the elastic modulus of the rock matrix;

[0009] Determining the fluid bulk modulus and fluid density based on the water saturation data and gas saturation data of the target layer, and determining the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus;

[0010] Determining the saturated rock density based on the fluid density and a predetermined rock matrix equivalent density determined based on the content and density of each mineral component in the target interval;

[0011] The longitudinal wave velocity and the shear wave velocity of the saturated rock are determined according to the elastic modulus and the density of the saturated rock, and the rock physics model of the target layer is determined according to the longitudinal wave velocity and the shear wave velocity of the saturated rock.

[0012] According to another aspect of the present invention, there is provided a device for constructing a rock physics model, the device comprising:

[0013] A compaction factor determination module is used to determine the compaction factor based on the porosity data and formation thickness data of the target layer;

[0014] a dry rock skeleton modulus determination module, configured to determine the rock matrix elastic modulus according to the content of each mineral component in the target layer and the compaction factor, and determine the dry rock skeleton elastic modulus according to the rock matrix elastic modulus;

[0015] a saturated rock modulus determination module, configured to determine the fluid bulk modulus and fluid density based on the water saturation data and the gas saturation data of the target layer, and to determine the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus;

[0016] A saturated rock density determination module is used to determine the saturated rock density based on the fluid density and a predetermined rock matrix equivalent density; the rock matrix equivalent density is determined based on the content and density of each mineral component in the target layer;

[0017] The rock physics model determination module is used to determine the longitudinal wave velocity and the shear wave velocity of the saturated rock according to the elastic modulus and the density of the saturated rock, and determine the rock physics model of the target layer according to the longitudinal wave velocity and the shear wave velocity of the saturated rock.

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

[0019] 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 execute the method for constructing a rock physics model described in any embodiment of the present invention.

[0020] 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 constructing a rock physical model according to any embodiment of the present invention when executed.

[0021] According to another aspect of the present invention, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the method for constructing a rock physics model according to any embodiment of the present invention.

[0022] The technical solution of the embodiment of the present invention determines a compaction factor based on porosity data and formation thickness data of the target layer; determines a rock matrix elastic modulus based on the content of each mineral component in the target layer and the compaction factor; and determines a dry rock skeleton elastic modulus based on the rock matrix elastic modulus; determines a fluid bulk modulus and fluid density based on water saturation data and gas saturation data of the target layer; and determines a saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus; determines a saturated rock density based on the fluid density and a predetermined rock matrix equivalent density; the rock matrix equivalent density is determined based on the content of each mineral component and the mineral density in the target layer; determines the longitudinal and shear wave velocities of the saturated rock based on the saturated rock elastic modulus and the saturated rock density; and determines a rock physics model of the target layer based on the longitudinal and shear wave velocities of the saturated rock. This technical solution solves the problem of poor reliability of rock physics modeling. By introducing the compaction factor into the determination of the rock matrix elastic modulus and calculating the rock matrix elastic modulus that takes into account mechanical compaction, the accuracy and effectiveness of the rock physics model can be improved.

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

[0024] 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.

[0025] Figure 1 is a flow chart of a method for constructing a rock physics model according to the first embodiment of the present invention;

[0026] Figure 2 This is a rock physics modeling flow chart provided according to a specific application scenario 1 of the present invention;

[0027] Figure 3 1 is a schematic structural diagram of a rock physics model construction device provided in accordance with a second embodiment of the present invention;

[0028] Figure 4 It is a structural schematic diagram of an electronic device for implementing the method for constructing a rock physics model according to an embodiment of the present invention. DETAILED DESCRIPTION

[0029] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions 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 embodiments described 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 making creative efforts should fall within the scope of protection of the present invention.

[0030] 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 interchangeable 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 clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices. The acquisition, storage, use, processing, etc. of data in the technical solution of this application comply with the relevant provisions of national laws and regulations.

[0031] Example 1

[0032] Figure 1A flowchart of a method for constructing a rock physics model is provided for the first embodiment of the present invention. This embodiment is applicable to oil and gas exploration scenarios, especially the construction of rock physics models. The method can be executed by a rock physics model construction device, which can be implemented in the form of hardware and / or software, and the device can be configured in an electronic device. Figure 1 As shown, the method includes:

[0033] S110 , determining a compaction factor based on porosity data and formation thickness data of the target layer.

[0034] This solution can be executed by electronic devices such as computers and servers, which can pre-acquire well logging data and core data from wells drilled in the target area. The well logging data can include physical property curves such as porosity, shale content, water saturation, and gas saturation for the wells drilled in the target area. The core data can be mineral composition data of the target layer obtained by X-ray diffraction experiments on core samples. For example, this data may include information such as the mineral composition types, the content and ratio of each type of mineral composition, etc.

[0035] Electronic equipment can identify lithologic interfaces based on well logging data from drilling within the target area and determine the formation thickness data for each layer. Based on the porosity data and formation thickness data for the target layer, the electronic equipment can calculate a compaction factor. The compaction factor can be used to characterize the mechanical compaction effect on the target layer. Specifically, the calculation formula for the compaction factor can be expressed as:

[0036]

[0037] Where c represents the compaction factor, h represents the formation thickness of the target layer, and Δh represents the formation thickness changed by mechanical compaction. Indicates the porosity of the target layer, represents the initial porosity determined based on the rock lithology of the target layer. Indicates the corrected porosity of the target layer.

[0038] S120. Determine the rock matrix elastic modulus based on the content of each mineral component in the target layer and the compaction factor, and determine the dry rock skeleton elastic modulus based on the rock matrix elastic modulus.

[0039] It is understood that the rock matrix elastic modulus includes the bulk modulus and shear modulus of the rock matrix. The electronic device can extract the content of each mineral component in the target layer from the core data. Based on the content of each mineral component in the target layer, the electronic device calculates the first bulk modulus and the first shear modulus using the Voigt upper limit formula, and calculates the second bulk modulus and the second shear modulus using the Ruess lower limit formula.

[0040] The electronic device can determine weight coefficients for the first bulk modulus and the second bulk modulus based on the compaction factor, and use the weighted sum of the first bulk modulus and the second bulk modulus as the first equivalent bulk modulus, i.e., the bulk modulus of the rock matrix. Similarly, the electronic device can determine weight coefficients for the first shear modulus and the second shear modulus based on the compaction factor, and use the weighted sum of the first shear modulus and the second shear modulus as the first equivalent shear modulus, i.e., the shear modulus of the rock matrix.

[0041] In this solution, the rock matrix elastic modulus includes a first equivalent bulk modulus and a first equivalent shear modulus;

[0042] The calculation formula of the first equivalent bulk modulus is:

[0043] K m =cK V +(1-c)K R ;

[0044] The calculation formula of the first equivalent shear modulus is:

[0045] μ m =cμ V +(1-c)μ R ;

[0046] Among them, K V represents the first bulk modulus, K R represents the second bulk modulus, μ V represents the first shear modulus, μ R represents the second shear modulus, i represents the mineral component identifier, N represents the number of mineral components, f i represents the content of mineral component i, μ i represents the shear modulus of mineral component i, K i represents the bulk modulus of mineral component i, c represents the compaction factor, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus.

[0047] After obtaining the rock matrix elastic modulus, the electronic device can calculate the dry rock skeleton elastic modulus based on the Kuster-Toksoz model. The dry rock skeleton elastic modulus includes a second equivalent bulk modulus and a second equivalent shear modulus, where the second equivalent bulk modulus is the bulk modulus of the dry rock skeleton, and the second equivalent shear modulus is the shear modulus of the dry rock skeleton.

[0048] The second equivalent bulk modulus satisfies:

[0049]

[0050] The second equivalent shear modulus satisfies:

[0051]

[0052] Among them, K KT represents the second equivalent bulk modulus, μ KT represents the second equivalent shear modulus, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus, P m and Q m are coefficients determined based on the pore shape. Indicates the porosity of the target layer.

[0053] In one feasible solution, the pore shape can be needle-shaped due to compaction, P m and Q m The expressions can be expressed as:

[0054]

[0055] Among them, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus, γ m represents the elastic coupling parameter corresponding to the pore shape,

[0056] S130. Determine the fluid bulk modulus and fluid density based on the water saturation data and gas saturation data of the target layer, and determine the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus.

[0057] Based on the Wood formula, the electronic device can determine the fluid bulk modulus and fluid density based on the water saturation data and gas saturation data of the target layer. The fluid bulk modulus calculation formula is:

[0058]

[0059] The fluid density calculation formula is:

[0060] ρ fl =S W ρ W +S G ρ G ;

[0061] Among them, K fl represents the bulk modulus of the fluid, ρ fl represents the fluid density, SW Indicates water saturation, S G Indicates gas saturation, ρ W represents the density of mineral water, ρ G Represents the density of natural gas.

[0062] After obtaining the fluid bulk modulus and fluid density, the electronic device can determine the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus based on the Gassmann equation. The saturated rock elastic modulus includes the saturated rock bulk modulus and shear modulus. Specifically, the saturated rock bulk modulus satisfies:

[0063]

[0064] The calculation formula of the shear modulus of the saturated rock is:

[0065] μ sat =μ KT ;

[0066] Among them, K sat represents the bulk modulus of saturated rock, μ sat represents the shear modulus of saturated rock, Indicates the porosity of the target layer, K KT represents the second equivalent bulk modulus, μ KT represents the second equivalent shear modulus, K m represents the first equivalent bulk modulus, K fl represents the bulk modulus of the fluid.

[0067] S140. Determine the saturated rock density based on the fluid density and a predetermined rock matrix equivalent density; the rock matrix equivalent density is determined based on the content and density of each mineral component in the target layer.

[0068] The electronic device can pre-calculate the rock matrix equivalent density based on the content and density of each mineral component in the target layer. Based on the rock matrix equivalent density and the fluid density, the electronic device can calculate the saturated rock density. Specifically, the saturated rock density calculation formula is:

[0069]

[0070] Among them, ρ sat represents the saturated rock density, represents the porosity of the target layer, ρ m represents the rock matrix equivalent density, i represents the mineral component identifier, N represents the number of mineral components, f i represents the content of mineral component i, ρ irepresents the density of mineral component i, ρ fl Represents the fluid density.

[0071] S150. Determine the longitudinal wave velocity and the transverse wave velocity of the saturated rock according to the elastic modulus and the density of the saturated rock, and determine a rock physics model of the target layer according to the longitudinal wave velocity and the transverse wave velocity of the saturated rock.

[0072] The electronic device can calculate the longitudinal wave velocity and the transverse wave velocity of the saturated rock according to the elastic modulus and the density of the saturated rock. Specifically, the longitudinal wave velocity calculation formula of the saturated rock is:

[0073]

[0074] The calculation formula for the shear wave velocity of the saturated rock is:

[0075]

[0076] Among them, K sat represents the bulk modulus of saturated rock, μ sat represents the shear modulus of saturated rock, ρ sat represents the saturated rock density.

[0077] After obtaining the P-wave velocity and S-wave velocity of the saturated rock, the electronic equipment can use the P-wave velocity data and S-wave velocity data in the logging data to verify the calculated P-wave velocity and S-wave velocity of the saturated rock and determine the accuracy of the rock physics model.

[0078] Specifically, if the difference between the P-wave velocity data obtained from well logging and the calculated P-wave velocity of the saturated rock is within a first preset error range, and the difference between the S-wave velocity data obtained from well logging and the calculated S-wave velocity of the saturated rock is within a second preset error range, the electronic device may determine that the rock physics model meets the accuracy requirements and output the rock physics model. If the difference between the P-wave velocity data obtained from well logging and the calculated P-wave velocity of the saturated rock is not within the first preset error range, or the difference between the S-wave velocity data obtained from well logging and the calculated S-wave velocity of the saturated rock is not within the second preset error range, the electronic device may adjust the rock physics parameters in the rock physics model that are reasonably assumed due to inability to measure, reconstruct the rock physics model, and verify the rock physics model until a rock physics model that meets the accuracy requirements is obtained.

[0079] The technical solution of the embodiment of the present invention determines a compaction factor based on porosity data and formation thickness data of the target layer; determines a rock matrix elastic modulus based on the content of each mineral component in the target layer and the compaction factor; and determines a dry rock skeleton elastic modulus based on the rock matrix elastic modulus; determines a fluid bulk modulus and fluid density based on water saturation data and gas saturation data of the target layer; and determines a saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus; determines a saturated rock density based on the fluid density and a predetermined rock matrix equivalent density; the rock matrix equivalent density is determined based on the content of each mineral component and the mineral density in the target layer; determines the longitudinal and shear wave velocities of the saturated rock based on the saturated rock elastic modulus and the saturated rock density; and determines a rock physics model of the target layer based on the longitudinal and shear wave velocities of the saturated rock. This technical solution solves the problem of poor reliability of rock physics modeling. By introducing the compaction factor into the determination of the rock matrix elastic modulus and calculating the rock matrix elastic modulus that takes into account mechanical compaction, the accuracy and effectiveness of the rock physics model can be improved.

[0080] Specific application scenario 1

[0081] Figure 2 This is a rock physics modeling flow chart provided according to a specific application scenario 1 of the present invention. Figure 2 As shown, the method includes:

[0082] Step 1: Laboratory core experiments

[0083] Obtain the composition of different minerals in the core sample and the content of each mineral.

[0084] Step 2: Obtain well logging data

[0085] Obtain physical property curves such as porosity, mud content, water saturation and gas saturation.

[0086] Step 3: Calculate the compaction factor:

[0087] The compaction factor is calculated based on the well logging porosity, initial porosity, and formation thickness; the compaction factor is used to characterize the degree of rock matrix compaction.

[0088] Step 4: Constructing the rock matrix in the rock physics model

[0089] Based on the compaction factor, the improved Voigt-Reuss-Hill model is used to calculate the rock matrix elastic modulus considering the mechanical compaction effect.

[0090] Step 5: Constructing the dry rock skeleton of the rock physics model

[0091] Based on the Kuster-Toksoz model, the elastic modulus of the rock matrix is used to calculate the elastic modulus of the dry rock skeleton.

[0092] Step 6: Construct the fluid in the rock physics model

[0093] The fluid bulk modulus is calculated based on the water saturation and gas saturation using the Wood formula, and the fluid density is calculated based on the measured mineral water density, natural gas density, water saturation, and gas saturation.

[0094] Step 7: Construct the saturated rock in the rock physics model

[0095] Based on the Gassmann equation, the elastic modulus of the dry rock skeleton is combined with the bulk modulus of the fluid to obtain the elastic modulus of the saturated rock, and the density of the saturated rock is obtained by combining the fluid density with the equivalent density of the rock matrix.

[0096] Step 8: Calculate the sound wave velocity in saturated rock

[0097] The longitudinal wave velocity and shear wave velocity of the saturated rock are calculated based on the elastic modulus and density of the saturated rock.

[0098] Compared with the existing technology, the advantages of this solution are:

[0099] 1. By improving the Voigt-Reuss-Hill model using the compaction factor and considering the impact of mechanical compaction on rock physics modeling, it is possible to effectively simulate the geological environment of oil and gas reservoirs and improve the adaptability and accuracy of the rock physics model.

[0100] 2. Due to the difficulty in measuring shear wave velocity, the amount of shear wave velocity data in actual work usually cannot meet research needs. The constructed rock physics model can be used to accurately calculate shear wave velocity data and provide reasonable and effective logging data for work area research.

[0101] Example 2

[0102] Figure 3 This is a schematic diagram of a rock physics model construction device provided in the second embodiment of the present invention. Figure 3 As shown, the device includes:

[0103] A compaction factor determination module 310 is used to determine a compaction factor based on porosity data and formation thickness data of a target layer;

[0104] a dry rock skeleton modulus determination module 320 for determining a rock matrix elastic modulus based on the content of each mineral component in the target layer and the compaction factor, and determining a dry rock skeleton elastic modulus based on the rock matrix elastic modulus;

[0105] a saturated rock modulus determination module 330 for determining the fluid bulk modulus and fluid density based on the water saturation data and the gas saturation data of the target layer, and determining the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus;

[0106] The saturated rock density determination module 340 is used to determine the saturated rock density based on the fluid density and a predetermined rock matrix equivalent density; the rock matrix equivalent density is determined based on the content and density of each mineral component in the target layer;

[0107] The rock physics model determination module 350 is used to determine the longitudinal wave velocity and the shear wave velocity of the saturated rock according to the elastic modulus and the density of the saturated rock, and determine the rock physics model of the target layer according to the longitudinal wave velocity and the shear wave velocity of the saturated rock.

[0108] In this solution, the calculation formula of the compaction factor is:

[0109]

[0110] Where c represents the compaction factor, h represents the formation thickness of the target layer, and Δh represents the formation thickness changed by mechanical compaction. Indicates the porosity of the target layer, represents the initial porosity determined based on the rock lithology of the target layer. Indicates the corrected porosity of the target layer.

[0111] On the basis of the above solution, the rock matrix elastic modulus includes a first equivalent bulk modulus and a first equivalent shear modulus;

[0112] The calculation formula of the first equivalent bulk modulus is:

[0113] K m =cK V +(1-c)K R ;

[0114] The calculation formula of the first equivalent shear modulus is:

[0115] μ m =cμ V +(1-c)μ R ;

[0116] Among them, K V represents the first bulk modulus, K R represents the second bulk modulus, μ V represents the first shear modulus, μ R represents the second shear modulus, i represents the mineral component identifier, N represents the number of mineral components, f i represents the content of mineral component i, μ i represents the shear modulus of mineral component i, K i represents the bulk modulus of mineral component i, c represents the compaction factor, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus.

[0117] In this embodiment, the dry rock skeleton elastic modulus includes a second equivalent bulk modulus and a second equivalent shear modulus;

[0118] The second equivalent bulk modulus satisfies:

[0119]

[0120] The second equivalent shear modulus satisfies:

[0121]

[0122] Among them, K KT represents the second equivalent bulk modulus, μ KT represents the second equivalent shear modulus, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus, P m and Q m are coefficients determined based on the pore shape. Indicates the porosity of the target layer.

[0123] In a feasible solution, the fluid bulk modulus calculation formula is:

[0124]

[0125] The fluid density calculation formula is:

[0126] ρ fl =S W ρ W +S G ρ G ;

[0127] Among them, K fl represents the bulk modulus of the fluid, ρ fl represents the fluid density, S W Indicates water saturation, S G Indicates gas saturation, ρ W represents the density of mineral water, ρ GRepresents the density of natural gas.

[0128] On the basis of the above scheme, the elastic modulus of the saturated rock includes the bulk modulus and shear modulus of the saturated rock;

[0129] The bulk modulus of the saturated rock satisfies:

[0130]

[0131] The calculation formula of the shear modulus of the saturated rock is:

[0132] μ sat =μ KT ;

[0133] The saturated rock density calculation formula is:

[0134]

[0135] Among them, K sat represents the bulk modulus of saturated rock, μ sat represents the shear modulus of saturated rock, ρ sat represents the saturated rock density, represents the porosity of the target layer, ρ m represents the rock matrix equivalent density, i represents the mineral component identifier, N represents the number of mineral components, f i represents the content of mineral component i, ρ i represents the density of mineral component i, K KT represents the second equivalent bulk modulus, μ KT represents the second equivalent shear modulus, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus, K fl represents the bulk modulus of the fluid, ρ fl Represents the fluid density.

[0136] The rock physics model construction device provided in the embodiment of the present invention can execute the rock physics model construction method provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0137] Example 3

[0138] Figure 4A schematic diagram of the structure of an electronic device 410 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can 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 required herein.

[0139] like Figure 4 As shown, the electronic device 410 includes at least one processor 411, and a memory connected to the at least one processor 411, such as a read-only memory (ROM) 412, a random access memory (RAM) 413, etc., wherein the memory stores a computer program that can be executed by the at least one processor, and the processor 411 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 412 or the computer program loaded from the storage unit 418 to the random access memory (RAM) 413. Various programs and data required for the operation of the electronic device 410 can also be stored in the RAM 413. The processor 411, ROM 412 and RAM 413 are connected to each other via a bus 414. An input / output (I / O) interface 415 is also connected to the bus 414.

[0140] Multiple components in electronic device 410 are connected to I / O interface 415, including an input unit 416, such as a keyboard, mouse, etc.; an output unit 417, such as various types of displays, speakers, etc.; a storage unit 418, such as a magnetic disk, optical disk, etc.; and a communication unit 419, such as a network card, modem, wireless communication transceiver, etc. The communication unit 419 allows electronic device 410 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0141] Processor 411 can be any general-purpose and / or specialized processing component with processing and computing capabilities. Some examples of processor 411 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, digital signal processors (DSPs), and any other suitable processor, controller, microcontroller, etc. Processor 411 executes the various methods and processes described above, such as the method for constructing a rock physics model.

[0142] In some embodiments, the rock physics model construction method can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as storage unit 418. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 410 via ROM 412 and / or communication unit 419. When the computer program is loaded into RAM 413 and executed by processor 411, one or more steps of the rock physics model construction method described above can be performed. Alternatively, in other embodiments, processor 411 can be configured to execute the rock physics model construction method in any other suitable manner (e.g., by means of firmware).

[0143] Various embodiments of the systems and techniques described 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), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system that includes 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.

[0144] Computer programs for implementing the methods of the present invention can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable rock physics model construction device, such that when executed by the processor, the computer programs implement the functions / operations specified in the flowcharts and / or block diagrams. The computer programs can be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0145] In the context of the present invention, computer-readable storage media can be tangible media that can contain or store a computer program for use with an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. Computer-readable storage media can include but are not limited to electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, computer-readable storage media can be machine-readable signal media. More specific examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0146] To provide interaction with a user, the systems and techniques described herein can 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 pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0147] The systems and techniques described herein can be implemented in a computing system that includes back-end 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 front-end components (e.g., a user computer with a graphical user interface or 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 back-end components, middleware components, or front-end components. The components of the system can 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.

[0148] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0149] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in the present invention can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved. This is not limited herein.

[0150] The above specific embodiments do not limit the scope of protection of the present invention. Those skilled in the art will appreciate that various modifications, combinations, sub-combinations, and substitutions may be made based on design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention are intended to be included within the scope of protection of the present invention.

Claims

1. A method for constructing a rock physics model, characterized in that: The method comprises: Determine the compaction factor based on the porosity data and formation thickness data of the target layer; Determine the elastic modulus of the rock matrix based on the content of each mineral component in the target layer and the compaction factor, and determine the elastic modulus of the dry rock skeleton based on the elastic modulus of the rock matrix; Determining the fluid bulk modulus and fluid density based on the water saturation data and gas saturation data of the target layer, and determining the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus; Determining the saturated rock density based on the fluid density and a predetermined rock matrix equivalent density determined based on the content and density of each mineral component in the target interval; The longitudinal wave velocity and the shear wave velocity of the saturated rock are determined according to the elastic modulus and the density of the saturated rock, and the rock physics model of the target layer is determined according to the longitudinal wave velocity and the shear wave velocity of the saturated rock.

2. The method according to claim 1, characterized in that The calculation formula of the compaction factor is: Where c represents the compaction factor, h represents the formation thickness of the target layer, and Δh represents the formation thickness changed by mechanical compaction. Indicates the porosity of the target layer, represents the initial porosity determined based on the rock lithology of the target layer. Indicates the corrected porosity of the target layer.

3. The method according to claim 2, characterized in that The rock matrix elastic modulus includes a first equivalent bulk modulus and a first equivalent shear modulus; The calculation formula of the first equivalent bulk modulus is: K m =cK V +(1-c)K R ; The calculation formula of the first equivalent shear modulus is: μ m =cμ V +(1-c)μ R ; Among them, K V represents the first bulk modulus, K R represents the second bulk modulus, μ V represents the first shear modulus, μ R represents the second shear modulus, i represents the mineral component identifier, N represents the number of mineral components, f i represents the content of mineral component i, μ i represents the shear modulus of mineral component i, K i represents the bulk modulus of mineral component i, c represents the compaction factor, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus.

4. The method according to claim 3, characterized in that The dry rock skeleton elastic modulus includes a second equivalent bulk modulus and a second equivalent shear modulus; The second equivalent bulk modulus satisfies: The second equivalent shear modulus satisfies: Among them, K KT represents the second equivalent bulk modulus, μ KT represents the second equivalent shear modulus, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus, P m and Q m are coefficients determined based on the pore shape. Indicates the porosity of the target layer.

5. The method according to claim 3, characterized in that The calculation formula of the fluid bulk modulus is: The fluid density calculation formula is: r fl =S W r W +S G r G ; Among them, K fl represents the bulk modulus of the fluid, ρ fl represents the fluid density, S W Indicates water saturation, S G Indicates gas saturation, ρ W represents the density of mineral water, ρ G Represents the density of natural gas.

6. The method according to claim 5, characterized in that The elastic modulus of the saturated rock includes the bulk modulus and shear modulus of the saturated rock; The bulk modulus of the saturated rock satisfies: The calculation formula of the shear modulus of the saturated rock is: m sat =μ KT ; The saturated rock density calculation formula is: Among them, K sat represents the bulk modulus of saturated rock, μ sat represents the shear modulus of saturated rock, ρ sat represents the saturated rock density, represents the porosity of the target layer, ρ m represents the rock matrix equivalent density, i represents the mineral component identifier, N represents the number of mineral components, f i represents the content of mineral component i, ρ i represents the density of mineral component i, K KT represents the second equivalent bulk modulus, μ KT represents the second equivalent shear modulus, K m represents the first equivalent bulk modulus, μ m represents the first equivalent shear modulus, K fl represents the bulk modulus of the fluid, ρ fl Represents the fluid density.

7. A device for constructing a rock physics model, characterized in that: The device comprises: A compaction factor determination module is used to determine the compaction factor based on the porosity data and formation thickness data of the target layer; a dry rock skeleton modulus determination module, configured to determine the rock matrix elastic modulus according to the content of each mineral component in the target layer and the compaction factor, and determine the dry rock skeleton elastic modulus according to the rock matrix elastic modulus; a saturated rock modulus determination module, configured to determine the fluid bulk modulus and fluid density based on the water saturation data and the gas saturation data of the target layer, and to determine the saturated rock elastic modulus based on the dry rock skeleton elastic modulus and the fluid bulk modulus; A saturated rock density determination module is used to determine the saturated rock density based on the fluid density and a predetermined rock matrix equivalent density; the rock matrix equivalent density is determined based on the content and density of each mineral component in the target layer; The rock physics model determination module is used to determine the longitudinal wave velocity and the shear wave velocity of the saturated rock according to the elastic modulus and the density of the saturated rock, and determine the rock physics model of the target layer according to the longitudinal wave velocity and the shear wave velocity of the saturated rock.

8. 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 to enable the at least one processor to execute the method for constructing a rock physical model according to any one of claims 1 to 6.

9. 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 constructing a rock physical model according to any one of claims 1 to 6 when executed.

10. A computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the computer program implements the method for constructing a rock physical model according to any one of claims 1 to 6.

Images