Earthquake rock physical modeling method and device for turbidite sandstone reservoir

Abstract

The invention relates to the technical field of rock physical modeling, and provides a turbidite sandstone reservoir earthquake rock physical modeling method and device. The method comprises the following steps: acquiring target turbidite sandstone reservoir information; determining a target rock matrix according to the target turbidite sandstone reservoir information; establishing a dry rock skeleton model by using the average coordination number and the critical porosity of the sandy mineral particles corresponding to the target turbidite sandstone reservoir; and establishing a fluid saturated rock model according to the dry rock skeleton model and a target rock matrix. According to the embodiment of the invention, the average coordination number and the critical porosity of the sandy mineral particles are comprehensively considered, so that the seismic rock physical modeling precision of the turbidite sandstone reservoir can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of rock physics modeling, and in particular to a method and device for seismic rock physics modeling of turbidite sandstone reservoirs. Background Art

[0002] Compared with conventional sandstone reservoirs, turbidite sandstone reservoirs are more susceptible to sedimentary effects such as compaction cementation, and the types of cementation are diverse, such as contact cementation, conventional cementation, and looseness. On the one hand, industry researchers have conducted the following analysis on cementation types: based on the types, causes, and properties of sandstone cementation, the effects of different types of cementation on rock physical properties are discussed; quartz crystallization in sandstone under different types of sedimentary and diagenetic backgrounds is studied, and the relationship between cementation type and the porosity and permeability of sandstone is emphasized; the controlling factors of quartz cementation in sandstone under deep burial conditions are explored, and the types of cementation, sources of cementing materials, and their effects on reservoir properties are analyzed. On the other hand, the existing turbidite sandstone reservoir modeling mainly focuses on the following aspects: based on the concept of rock physics quantitative plate, the relationship between seismic velocity and rock properties of turbidite sandstone and carbonate reservoirs is discussed, and the importance of using quantitative plate to interpret seismic data is emphasized; rock physics modeling of turbidite sandstone reservoirs is carried out based on the correlation between seismic data and rock properties, and the modeling results are used to predict reservoir properties; a comprehensive seismic reservoir characterization method is proposed, and the physical property modeling of turbidite sandstone reservoirs is discussed by combining seismic interpretation and modeling; multi-scale rock physics modeling and uncertainty quantification methods in complex turbidite sandstone reservoirs are discussed, and the relationship between pore structure and geophysical properties is analyzed. However, the existing technology lacks a rock physics modeling method for turbidite sandstone reservoirs that comprehensively considers the average coordination number and critical porosity of sandy mineral particles, and there is a problem of low modeling accuracy. Summary of the invention

[0003] In view of the low accuracy of rock physics modeling of turbidite sandstone reservoirs at present, which does not take into account the average coordination number and critical porosity of sandy mineral particles, the present solution is proposed to overcome the above problems or at least partially solve the above problems.

[0004] On the one hand, an object of some embodiments of this specification is to provide a method for seismic rock physics modeling of turbidite sandstone reservoirs, the method comprising:

[0005] Obtain target turbidite sandstone reservoir information;

[0006] Determining a target rock matrix according to the target turbidite sandstone reservoir information;

[0007] A dry rock skeleton model is established using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir;

[0008] A fluid-saturated rock model is established based on the dry rock skeleton model and the target rock matrix.

[0009] Furthermore, the target turbidite sandstone reservoir information at least includes well logging curves and rock matrix mud and sand component modulus information.

[0010] Further, determining the target rock matrix according to the target turbidite sandstone reservoir information includes:

[0011] Determining the shale content of the target turbidite sandstone reservoir by using the well logging curve;

[0012] According to the rock matrix mud and sand component modulus information and the mud content, the rock matrix modulus interval information is determined; wherein the rock matrix mud and sand component modulus information includes the bulk modulus of the sandy mineral component constituting the rock matrix, the bulk modulus of the muddy mineral component, the shear modulus of the sandy mineral component and the shear modulus of the muddy mineral component; the rock matrix modulus interval information includes the upper limit and lower limit of the rock matrix bulk modulus, and the upper limit and lower limit of the rock matrix shear modulus;

[0013] The bulk modulus and shear modulus of the target rock matrix are determined according to the rock matrix modulus interval information.

[0014] Furthermore, a dry rock skeleton model is established using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir, including:

[0015] Determine the bulk modulus and shear modulus of the cemented rock using the average coordination number and critical porosity of the sandy mineral particles;

[0016] The dry rock skeleton model is established according to the bulk modulus and shear modulus of the cemented rock.

[0017] Furthermore, the bulk modulus and shear modulus of the cemented rock are determined according to the following formula using the average coordination number and critical porosity of the sandy mineral particles:

[0018] in,

[0019] in,

[0020] Among them, K b is the bulk modulus of cemented rock, C is the average coordination number of sandy mineral particles, φ c is the critical porosity of the rock in its initial sedimentary state, M shale is the longitudinal wave modulus of mud minerals in rocks, ρ shale is the density of argillaceous minerals, V Pshaleis the longitudinal wave velocity of argillaceous minerals, S n is the normal flexibility of the cemented sandy mineral particles, G b is the shear modulus of cemented rock, G shale is the shear modulus of argillaceous minerals, V Sshale is the shear wave velocity of argillaceous minerals, S τ It is the longitudinal flexibility of cemented sandy mineral particles.

[0021] Furthermore, the normal flexibility and longitudinal flexibility of the cemented sandy mineral particles are determined using the following formula:

[0022] S n =A n (Λ n )α 2 +B n (Λ n )α+C n (Λ n )

[0023] in,

[0024]

[0025]

[0026]

[0027] S τ =A τ (Λ τ ,v)α 2 +B τ (Λ τ ,v)α+C τ (Λ τ ,v)

[0028] in,

[0029]

[0030]

[0031]

[0032] Among them, A n ,Λ n , B n , C n , A τ ,Λ τ , B τ , C τare all intermediate variables, α represents the cementation mode of sandy mineral particles, π is the circumference, G sand is the shear modulus of the matrix mineral, v sand represents the Poisson's ratio of sandy minerals, v shale represents the Poisson's ratio of mud minerals, and v represents the Poisson's ratio taking into account the bulk modulus and shear modulus of cemented rock.

[0033] Furthermore, the cementation mode of the sandy mineral particles is determined according to the average coordination number of the sandy mineral particles and the effective reservoir pressure using the following formula:

[0034]

[0035]

[0036] Among them, P c is the effective reservoir pressure, S is the percentage of pore space occupied by cement, and φ is the porosity of fluid-saturated rock.

[0037] Furthermore, according to the bulk modulus and shear modulus of the cemented rock, the dry rock skeleton model is established using the following formula:

[0038]

[0039]

[0040]

[0041] Among them, K dry is the bulk modulus of the dry rock skeleton, φ b is the cementation porosity, K ma is the bulk modulus of the target rock matrix, G ma is the shear modulus of the target rock matrix, G dry is the shear modulus of the dry rock skeleton, and z is an intermediate quantity.

[0042] Further, a fluid-saturated rock model is established based on the dry rock skeleton model and the target rock matrix, including:

[0043] Determining a mixed fluid corresponding to the fluid-saturated rock model according to the effects of pressure and temperature changes on the fluid;

[0044] Determine the bulk modulus and shear modulus of the fluid-saturated rock based on the bulk modulus of the mixed fluid, the dry rock skeleton model, and the target rock matrix;

[0045] According to the bulk modulus and shear modulus of fluid-saturated rock and the density of the mixed fluid, the longitudinal wave velocity and shear wave velocity of the equivalent fluid-saturated rock after the rock physical fluid replacement are determined to establish a fluid-saturated rock model.

[0046] Further, according to the influence of pressure and temperature changes on the fluid, determining the mixed fluid corresponding to the fluid saturated rock model includes:

[0047] Determine the flow rate and density of oil and the flow rate and density of brine based on the effects of pressure and temperature changes on the fluid;

[0048] Based on the flow rate and density of the oil, and the flow rate and density of the brine, the bulk modulus and density of the mixed fluid consisting of the oil and the brine are determined.

[0049] Furthermore, according to the influence of pressure and temperature changes on the fluid, the flow rate and density of the oil, as well as the flow rate and density of the brine, are determined using the following formula:

[0050]

[0051] ρ Oil =[ρ0+(0.00277P-1.71×10 -7 P 3 )(ρ0-1.15) 2 +3.49×10 -4 P] / [0.972+3.81×10 -4 (T+17.78) 1.175 ]

[0052] Among them, V Oil is the oil flow rate considering the changes in temperature and pressure, ρ0 is the oil density at normal temperature and pressure, T is the temperature, P is the pressure, ρ Oil Oil density to account for temperature and pressure changes;

[0053] V brine =V water +S water 1.5 (780-10P+0.16P 2 )-820S water 2 +S water (1170-9.6T+0.055T 2 -8.5×10 -5 T 3 +2.6P-0.0029TP-0.0476P 2 )

[0054] ρ water =1+1×10 -6 (-80T-3.3T 2 +0.00175T 3 +489P-2TP+0.016T2 P-1.3×10 -5 T 3 P-0.333P 2 -0.002TP 2 )

[0055] ρ brine =ρ water +0.668S water +0.44S water 2 +10 -6 S water [300P-2400PS water +T(80+3T-3300S water -13P+47PS water )]

[0056] Among them, V brine V is the brine flow rate considering the changes in temperature and pressure. water is the pure water speed, S water is the salt concentration in the brine, ρ water is the density of pure water considering temperature and pressure changes, ρ brine is the density of salt water taking into account changes in temperature and pressure.

[0057] Furthermore, based on the flow rate and density of the oil, and the flow rate and density of the brine, the bulk modulus and density of the mixed fluid consisting of the oil and brine are calculated using the following formula:

[0058]

[0059] ρ fl =S brine ρ brine +S Oil ρ Oil

[0060] Among them, K fl is the bulk modulus of the mixed fluid in the pore, S brine is the brine saturation, ρ brine V is the density of salt water considering the changes in temperature and pressure. brine is the brine flow rate considering the changes in temperature and pressure, S Oil is the oil saturation, ρ Oil V is the oil density considering the changes in temperature and pressure. Oil is the oil flow rate considering temperature and pressure changes, ρ fl is the density of the mixed fluid in the pores.

[0061] Furthermore, based on the bulk modulus of the mixed fluid, the dry rock skeleton model and the target rock matrix, the bulk modulus and shear modulus of the fluid-saturated rock are determined using the following formula:

[0062]

[0063] G sat =G dry

[0064] Among them, K sat is the bulk modulus of fluid-saturated rock, K dry is the bulk modulus of the dry rock skeleton, K ma is the bulk modulus of the target rock matrix, φ is the porosity of the fluid-saturated rock, K fl is the bulk modulus of the mixed fluid in the pore, G dry is the shear modulus of the dry rock skeleton, G sat is the shear modulus of fluid-saturated rock.

[0065] On the other hand, some embodiments of the present specification further provide a turbidite sandstone reservoir seismic rock physics modeling device, the device comprising:

[0066] A receiving module, used for acquiring target turbidite sandstone reservoir information;

[0067] A rock matrix determination module, used to determine a target rock matrix according to the target turbidite sandstone reservoir information;

[0068] A dry rock skeleton determination module is used to establish a dry rock skeleton model using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir;

[0069] A modeling module is used to establish a fluid-saturated rock model based on the dry rock skeleton model and the target rock matrix.

[0070] On the other hand, some embodiments of the present specification further provide a computer device, including a memory, a processor, and a computer program stored in the memory, wherein when the computer program is executed by the processor, the instructions of the above method are executed.

[0071] On the other hand, some embodiments of the present specification further provide a computer storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor of a computer device, the instructions of the above method are executed.

[0072] On the other hand, some embodiments of the present specification further provide a computer program product, wherein the computer program product comprises a computer program, and when the computer program is executed by a processor of a computer device, the instructions of the above method are executed.

[0073] One or more technical solutions provided by some embodiments of this specification have at least the following technical effects:

[0074] In the embodiment of the present specification, when performing seismic rock physics modeling on a turbidite sandstone reservoir, target turbidite sandstone reservoir information is obtained to determine a target rock matrix, and a dry rock skeleton model is established that comprehensively considers the average coordination number and critical porosity of sandy mineral particles of cemented rock based on the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir information. Thereafter, fluid-saturated rock is determined based on the dry rock skeleton model and the target rock matrix, thereby improving the rock physics modeling accuracy of the turbidite sandstone reservoir.

[0075] The above description is only an overview of the technical solutions of some embodiments of this specification. In order to more clearly understand the technical means of some embodiments of this specification, they can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of some embodiments of this specification more obvious and easy to understand, the specific implementation methods of some embodiments of this specification are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0076] In order to more clearly illustrate some embodiments of this specification or technical solutions in the prior art, the following briefly introduces the drawings required for the embodiments or prior art descriptions. Obviously, the drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative labor. In the drawings:

[0077] Figure 1 A schematic diagram of an implementation system of a turbidite sandstone reservoir seismic rock physics modeling method in some embodiments of this specification is shown;

[0078] Figure 2 A flowchart of a method for seismic rock physics modeling of turbidite sandstone reservoirs in some embodiments of this specification is shown;

[0079] Figure 3 A schematic diagram of the process of seismic rock physics modeling of turbidite sandstone reservoirs in some embodiments of this specification;

[0080] Figure 4 A schematic diagram of the steps for determining a target rock matrix in some embodiments of this specification;

[0081] Figure 5 A schematic diagram of the steps of establishing a dry rock skeleton model in some embodiments of this specification;

[0082] Figure 6 A comparative schematic diagram of bonding types in some embodiments of this specification;

[0083] Figure 7 A schematic diagram of the steps of establishing a fluid-saturated rock model in some embodiments of this specification;

[0084] Figure 8 A schematic diagram of the steps of determining a mixed fluid in some embodiments of this specification;

[0085] Figure 9a This is a first schematic diagram of comparison between actual well logging data and rock physics modeling results in some embodiments of this specification;

[0086] Figure 9b A second schematic diagram of comparison between actual well logging data and rock physics modeling results in some embodiments of this specification;

[0087] Fig.10 This is a schematic structural diagram of a turbidite sandstone reservoir seismic rock physics modeling device in some embodiments of this specification;

[0088] Fig.11 This is a schematic diagram of the computer device structure provided in some embodiments of this specification.

[0089] [Description of Reference Numerals]

[0090] 101. Terminal;

[0091] 102. Server;

[0092] 1001, receiving module;

[0093] 1002. Rock matrix determination module;

[0094] 1003. Dry rock skeleton determination module;

[0095] 1004, modeling module;

[0096] 1102. Computer equipment;

[0097] 1104. Processor;

[0098] 1106. Memory;

[0099] 1108. Driving mechanism;

[0100] 1110, input / output interface;

[0101] 1112. Input device;

[0102] 1114. Output device;

[0103] 1116. Presentation equipment;

[0104] 1118. Graphical user interface;

[0105] 1120, network interface;

[0106] 1122. Communication link;

[0107] 1124. Communication bus. DETAILED DESCRIPTION

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

[0109] It should be noted that the terms "first", "second", etc. in the specification and claims of this article and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of this article described here can be implemented in an order other than those illustrated or described here. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, device, product or equipment 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 equipment.

[0110] It should be noted that the acquisition, storage, use, and processing of data in the technical solution of this application are in compliance with the relevant provisions of relevant laws and regulations.

[0111] like Figure 1The schematic diagram of the implementation system of a turbidite sandstone reservoir seismic rock physics modeling method according to an embodiment of the present invention is shown, which may include: a terminal 101 and a server 102, wherein the terminal 101 and the server 102 communicate with each other through a network, and the network may include a local area network (LAN), a wide area network (WAN), the Internet or a combination thereof, and is connected to a website, a user device (such as a computing device) and a back-end system. A staff member may send a turbidite sandstone reservoir seismic rock physics modeling request to the server 102 through the terminal 101, and after receiving the turbidite sandstone reservoir seismic rock physics modeling request, the server 102 calls the target turbidite sandstone reservoir information in the database for calculation and processing, obtains a fluid saturated rock modeling result, and sends the fluid saturated rock modeling result to the terminal 101, so that the staff member processes the business according to the fluid saturated rock modeling result.

[0112] In the embodiments of this specification, the server 102 can be an independent physical server, or a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, content delivery networks (CDN, Content Delivery Network), and big data and artificial intelligence platforms.

[0113] In an optional embodiment, the terminal 101 may include but is not limited to electronic devices such as self-service terminal devices, desktop computers, tablet computers, laptop computers, smart wearable devices, etc. Optionally, the operating system running on the electronic device may include but is not limited to Android system, IOS system, Linux, Windows, etc. Of course, the terminal 101 is not limited to the above-mentioned electronic devices with a certain entity, and it can also be software running in the above-mentioned electronic devices.

[0114] In addition, it should be noted that Figure 1 What is shown is only an application environment provided by the present disclosure. In actual application, multiple terminals 101 may be included, and this specification does not limit this.

[0115] Figure 2It is a flowchart of a method for seismic rock physics modeling of turbidite sandstone reservoirs provided by an embodiment of the present invention. This specification provides method operation steps as described in the embodiment or flowchart, but may include more or fewer operation steps based on conventional or non-creative labor. The order of steps listed in the embodiment is only one way of executing the steps among many orders, and does not represent the only order of execution. When the system or device product is executed in practice, it can be executed in the order or in parallel according to the method shown in the embodiment or the accompanying drawings. Specifically, Figure 2 As shown, applied to the above-mentioned server side, the method may include:

[0116] S201: Acquire target turbidite sandstone reservoir information;

[0117] S202: determining a target rock matrix according to the target turbidite sandstone reservoir information;

[0118] S203: establishing a dry rock skeleton model using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir;

[0119] S204: Establishing a fluid-saturated rock model according to the dry rock skeleton model and the target rock matrix.

[0120] In the embodiment of the present specification, when performing seismic rock physics modeling on a turbidite sandstone reservoir, target turbidite sandstone reservoir information is obtained to determine a target rock matrix, and a dry rock skeleton model is established that comprehensively considers the average coordination number and critical porosity of sandy mineral particles of cemented rock based on the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir information. Thereafter, fluid-saturated rock is determined based on the dry rock skeleton model and the target rock matrix, thereby improving the rock physics modeling accuracy of the turbidite sandstone reservoir.

[0121] It can be understood that, in some embodiments, the target turbidite sandstone reservoir information at least includes well logging curves and rock matrix mud and sand component modulus information. Specifically, the well logging curve is a curve drawn from the logging information obtained by various logging methods to reflect the different lithology and stratigraphic characteristics in the formation. The target rock matrix can be determined according to the well logging curve in the target turbidite sandstone reservoir information. Then, the dry rock skeleton model can be determined according to the average coordination number and critical porosity of the sandy mineral particles corresponding to the target turbidite sandstone reservoir. The target rock matrix and the dry rock skeleton model are used to construct a turbidite sandstone seismic rock physics model that considers the particle coordination number and critical porosity. It is necessary to establish a quantitative relationship between the physical properties of turbidite sandstone (including porosity and saturation), pore structure, and formation pressure and elastic parameters (including P-wave velocity and S-wave velocity), so as to provide a theoretical model for the seismic response characteristic analysis of turbidite sandstone, reservoir prediction, and fluid identification, and to provide reliable technology and methods for the exploration and development of turbidite sandstone.

[0122] Further, see Attachment Figure 3 In some embodiments, rock matrix is ​​a term in petrology and seepage mechanics, which means the rock part in the fractured oil reservoir. A variety of minerals are mixed to form the rock matrix. Intuitively, the pores cut and reshape the rock matrix to form a dry rock skeleton. Specifically, when establishing a dry rock skeleton model, in addition to considering the rock matrix, it is also necessary to consider the influence of two types of factors, the average coordination number of sandy mineral particles and the critical porosity, on the dry rock skeleton. After obtaining the dry rock skeleton model, when the dry rock skeleton model is replaced with rock physical fluid, in addition to analyzing the bulk modulus of the mixed fluid, etc., since the basis for establishing saturated rock and dry rock skeleton (which can be understood as a "foundation") is the target rock matrix, it is also necessary to determine the final fluid-saturated rock model in combination with the bulk modulus of the target rock matrix to improve the modeling accuracy of the final fluid-saturated rock model.

[0123] See attached Figure 4 In some embodiments, determining a target rock matrix according to the target turbidite sandstone reservoir information may include:

[0124] S401: Determine the mud content of the target turbidite sandstone reservoir using the well logging curve;

[0125] S402: determining rock matrix modulus interval information according to the rock matrix mud and sand component modulus information and the mud content; wherein the rock matrix mud and sand component modulus information includes the bulk modulus of the sandy mineral component constituting the rock matrix, the bulk modulus of the muddy mineral component, the shear modulus of the sandy mineral component, and the shear modulus of the muddy mineral component; the rock matrix modulus interval information includes the upper limit and lower limit of the rock matrix bulk modulus, and the upper limit and lower limit of the rock matrix shear modulus;

[0126] S403: Determine the bulk modulus and shear modulus of the target rock matrix according to the rock matrix modulus interval information.

[0127] It can be understood that, in some embodiments, the shale content is an important basic data for calculating and evaluating the parameters of turbidite sandstone reservoirs. It not only reflects the lithology of the formation, but is also closely related to the effective porosity, water saturation, irreducible water saturation, permeability and other parameters of the reservoir. The shale content of the formation can be determined by natural gamma logging, neutron and density logging, etc. The specific calculation method of the shale content is not limited in this article. Then, based on the shale content, the bulk modulus of the sandy mineral components constituting the rock matrix, the bulk modulus of the shale mineral components, the shear modulus of the sandy mineral components and the shear modulus of the shale mineral components, the rock matrix modulus interval information can be determined. The rock matrix modulus interval information includes the upper and lower limits of the rock matrix bulk modulus, and the upper and lower limits of the rock matrix shear modulus. Specifically, in some embodiments, the bulk modulus and shear modulus of the target rock matrix can be determined by the following formula, so as to accurately determine the target rock matrix:

[0128]

[0129] K V =V shale K shale +(1-V shale )K sand

[0130]

[0131]

[0132] G V =V shale G shale +(1-V shale )G sand

[0133]

[0134] Among them, K ma is the bulk modulus of the target rock matrix, K V is the upper limit of the rock matrix bulk modulus, V shale is the mud content, K shale is the bulk modulus of the argillaceous mineral components that constitute the rock matrix, K sand is the bulk modulus of the sandy mineral components that constitute the rock matrix, K R is the lower limit of the rock matrix bulk modulus, G ma is the shear modulus of the target rock matrix, G V is the upper limit of the rock matrix shear modulus, G R is the lower limit of rock matrix shear modulus, G shale is the shear modulus of the argillaceous mineral components that constitute the rock matrix, G sandis the shear modulus of the sandy mineral components that make up the rock matrix.

[0135] See attached Figure 5 In some embodiments, using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir to establish a dry rock skeleton model may include:

[0136] S501: Determine the bulk modulus and shear modulus of the cemented rock using the average coordination number and critical porosity of the sandy mineral particles;

[0137] S502: Establishing the dry rock skeleton model according to the bulk modulus and shear modulus of the cemented rock.

[0138] It can be understood that in some embodiments, the bulk modulus and shear modulus of cemented rock can be determined using the average coordination number and critical porosity of sandy mineral particles according to the following formula:

[0139] in,

[0140] in,

[0141] Among them, K b is the bulk modulus of cemented rock, C is the average coordination number of sandy mineral particles, φ c is the critical porosity of the rock in its initial sedimentary state, M shale is the longitudinal wave modulus of mud minerals in rocks, ρ shale is the density of argillaceous minerals, V Pshale is the longitudinal wave velocity of argillaceous minerals, S n is the normal flexibility of the cemented sandy mineral particles, G b is the shear modulus of cemented rock, G shale is the shear modulus of argillaceous minerals, V Sshale is the shear wave velocity of argillaceous minerals, S τ It is the longitudinal flexibility of the cemented sandy mineral particles. It can be understood that, in some embodiments, the average coordination number C of the sandy mineral particles means how many sandy mineral particles are in contact with each sandy mineral particle around it. In the process of seismic rock physics modeling of turbidite sandstone reservoirs, the coordination number of sandy mineral particles and the critical porosity are considered, which changes the traditional seismic rock physics modeling of turbidite sandstone reservoirs under a fixed critical porosity value, thereby improving the modeling accuracy of seismic rock physics of turbidite sandstone reservoirs.

[0142] Further, in some embodiments, after the bulk modulus and shear modulus of the cemented rock are determined, the dry rock skeleton model can be established using the following formula according to the bulk modulus and shear modulus of the cemented rock:

[0143]

[0144]

[0145]

[0146] Among them, K dry is the bulk modulus of the dry rock skeleton, φ is the porosity of the fluid-saturated rock. The porosity of the fluid-saturated rock φ is an important parameter when establishing the dry rock skeleton model. b is the cementation porosity, K ma is the bulk modulus of the target rock matrix, G ma is the shear modulus of the target rock matrix, G dry is the shear modulus of the dry rock skeleton, and z is an intermediate quantity used to simplify the formula expression. It can be understood that in some embodiments, the cementation porosity φ b The meaning is that the rock porosity gradually decreases with the cementation effect. When it decreases to φ b When the cementation is completed and compaction begins, the corresponding porosity is φ b In some typical embodiments, the dry rock skeleton model is mainly based on the partial cementation model, so that the dry rock skeleton model can reflect the mechanical characteristics of the target reservoir. Figure 6 The schematic diagram of the partial cementation model is shown in FIG. 1 , where the horizontal axis is Porosity and the vertical axis is P-wave velocity. Under the three different cementation modes of uncemented, contact cemented and partially cemented, the P-wave velocity on the vertical axis varies with the change of rock porosity. Specifically, in some embodiments, in order to improve the modeling accuracy of the dry rock skeleton model, the cementation porosity φ is calculated. b It is also necessary to consider the normal flexibility of the cemented sandy mineral particles, the bulk modulus of the cement, the shear modulus of the cement, the effective pressure of the reservoir, etc.

[0147] Furthermore, in some embodiments, the normal flexibility and longitudinal flexibility of the cementing material, i.e., the muddy mineral particles, and the cemented sandy mineral particles are related to the cementing mode of the sandy mineral particles, the shear modulus of the cementing material, and the shear modulus of the matrix mineral, and can be determined using the following formula:

[0148] S n =A n (Λ n )α 2 +B n (Λ n )α+C n (Λ n )

[0149] in,

[0150]

[0151]

[0152]

[0153] S τ =A τ (Λ τ ,v)α 2 +B τ (Λ τ ,v)α+C τ (Λ τ ,v)

[0154] in,

[0155]

[0156]

[0157]

[0158] Among them, S n is the normal flexibility of the cemented sandy mineral particles, A n ,Λ n , B n , C n , A τ ,Λ τ , B τ and C τ are all intermediate variables, α represents the cementation mode of sandy mineral particles, G shale is the shear modulus of argillaceous minerals, π is the circumference, G sand is the shear modulus of the matrix mineral, v sand represents the Poisson's ratio of sandy minerals, v shale represents the Poisson's ratio of mud minerals, S τ is the longitudinal flexibility of the cemented sandy mineral particles, and v represents the Poisson's ratio considering the bulk modulus and shear modulus of the cemented rock.

[0159] Furthermore, in some embodiments, the cementation mode α of the sandy mineral particles is related to the average coordination number of the sandy mineral particles and the effective pressure of the reservoir, and can be determined using the following formula:

[0160]

[0161] C=0.5*(777.1P c -0.001545 -770.7)+0.5*(-4.457Pc -0.2724 +9.401)

[0162] Among them, P c is the effective reservoir pressure, S is the percentage of pore space occupied by cement, φ is the porosity of fluid-saturated rock, and φ c The critical porosity is the porosity that takes into account the particle arrangement, mineral sorting degree and reservoir effective pressure, that is, the critical porosity in the initial deposition state. Specifically, in some embodiments, the average coordination number of sandy mineral particles is related to the effective pressure of the reservoir. Generally, the greater the effective pressure of the reservoir, the higher the average coordination number of sandy mineral particles. The average coordination number of sandy mineral particles in the reservoir is not a fixed value, and the particle arrangement, mineral sorting degree and reservoir effective pressure will affect the initial deposition porosity of the rock (that is, the critical porosity of the rock φ c ), it is traditionally believed that the critical porosity of rock is a fixed value, such as the critical porosity of sandstone is often taken as 0.4. In fact, the porosity of the rock at the initial deposition is affected by the particle arrangement, mineral sorting degree and effective reservoir pressure factors. Therefore, the critical porosity of sandstone is φ c It is a dynamically changing value. Compared with the traditional dry rock skeleton, the dry rock skeleton constructed by the dynamically changing critical porosity and the average coordination number of sandy mineral particles takes into account more comprehensive influencing factors and has higher modeling accuracy.

[0163] See attached Figure 7 In some embodiments, establishing a fluid-saturated rock model based on the dry rock skeleton model and the target rock matrix may include:

[0164] S701: Determine a mixed fluid corresponding to the fluid-saturated rock model according to the influence of pressure and temperature changes on the fluid;

[0165] S702: Determine the bulk modulus and shear modulus of the fluid-saturated rock based on the bulk modulus of the mixed fluid, the dry rock skeleton model, and the target rock matrix;

[0166] S703: According to the bulk modulus and shear modulus of the fluid-saturated rock and the density of the mixed fluid, the longitudinal wave velocity and the transverse wave velocity of the equivalent fluid-saturated rock after the rock physical fluid is replaced are determined to establish a fluid-saturated rock model.

[0167] It can be understood that in some embodiments, the pressure and temperature at different positions in the target turbidite sandstone reservoir are different. As the pressure and temperature change, the bulk modulus and density of the mixed fluid also change accordingly. According to the dry rock skeleton model and the target rock matrix, the mixed fluid is used for fluid replacement to determine the bulk modulus and shear modulus of the fluid-saturated rock model, as well as the longitudinal wave velocity and transverse wave velocity of the equivalent fluid-saturated rock. A fluid-saturated rock model that considers the average coordination number and critical porosity of sandy mineral particles is established, thereby improving the modeling accuracy of the fluid-saturated rock model corresponding to the target turbidite sandstone reservoir.

[0168] See attached Figure 8 In some embodiments, determining a mixed fluid corresponding to the fluid-saturated rock model based on the effect of pressure and temperature changes on the fluid may include:

[0169] S801: Determine the flow rate and density of oil, and the flow rate and density of salt water according to the effects of pressure and temperature changes on the fluid;

[0170] S802: Based on the flow rate and density of the oil and the flow rate and density of the brine, determine the bulk modulus and density of the mixed fluid consisting of the oil and the brine.

[0171] It can be understood that, in some embodiments, according to the influence of pressure and temperature changes on the fluid, the following formula can be used to accurately determine the flow rate and density of the oil and the flow rate and density of the brine as the pressure and temperature change:

[0172]

[0173] ρ Oil =[ρ0+(0.00277P-1.71×10 -7 P 3 )(ρ0-1.15) 2 +3.49×10 -4 P] / [0.972+3.81×10 -4 (T+17.78) 1.175 ]

[0174] Among them, V Oil is the oil flow rate considering the changes in temperature and pressure, ρ0 is the oil density at normal temperature and pressure, T is the temperature, P is the pressure, ρ Oil Oil density to account for temperature and pressure changes;

[0175] V brine =V water +S water 1.5 (780-10P+0.16P 2 )-820S water2 +S water (1170-9.6T+0.055T 2 -8.5×10 -5 T 3 +2.6P-0.0029TP-0.0476P 2 )

[0176] ρ water =1+1×10 -6 (-80T-3.3T 2 +0.00175T 3 +489P-2TP+0.016T 2 P-1.3×10 -5 T 3 P-0.333P 2 -0.002TP 2 )

[0177] ρ brine =ρ water +0.668S water +0.44S water 2 +10 -6 S water [300P-2400PS water +T(80+3T-3300S water -13P+47PS water )]

[0178] Among them, V brine V is the brine flow rate considering the changes in temperature and pressure. water is the pure water speed, S water is the salt concentration in the brine, ρ water is the density of pure water considering temperature and pressure changes, ρ brine is the density of salt water taking into account changes in temperature and pressure.

[0179] Further, in some embodiments, after obtaining the flow rate and density of oil and the flow rate and density of brine that vary with pressure and temperature, since the mixed fluid corresponding to the target turbidite sandstone reservoir is a mixture of oil and brine, the bulk modulus and density of the mixed fluid composed of oil and brine can be calculated using the following formula:

[0180]

[0181] ρ fl =S brine ρ brine +S Oil ρ Oil

[0182] Among them, K fl is the bulk modulus of the mixed fluid in the pore, S brine is the brine saturation, ρ brine V is the density of salt water considering the changes in temperature and pressure. brine is the brine flow rate considering the changes in temperature and pressure, S Oil is the oil saturation, ρ Oil V is the oil density considering the changes in temperature and pressure. Oil is the oil flow rate considering temperature and pressure changes, ρ fl is the density of the mixed fluid in the pores.

[0183] Further, in some embodiments, based on the bulk modulus of the mixed fluid, the dry rock skeleton model and the target rock matrix, the bulk modulus and shear modulus of the fluid-saturated rock can be determined using the following formula:

[0184]

[0185] G sat =G dry

[0186] Among them, K sat is the bulk modulus of fluid-saturated rock, K dry is the bulk modulus of the dry rock skeleton, K fl is the bulk modulus of the mixed fluid, K ma is the bulk modulus of the target rock matrix, φ is the porosity of the fluid-saturated rock, G dry is the shear modulus of the dry rock skeleton, G sat is the shear modulus of fluid-saturated rock.

[0187] Further, in some embodiments, based on the bulk modulus and shear modulus of the fluid-saturated rock and the density of the mixed fluid, the following formula can be used to determine the longitudinal wave velocity and transverse wave velocity of the equivalent fluid-saturated rock after the rock physical fluid is replaced:

[0188]

[0189]

[0190] Among them, V p and V s are the P-wave velocity and S-wave velocity of equivalent fluid saturated rock, K sat is the bulk modulus of fluid-saturated rock, G sat is the shear modulus of fluid-saturated rock, ρ sat is the density of fluid-saturated rock, thus completing the modeling of fluid-saturated rock model.

[0191] Further, see Attachment Figure 9a The conventional method is used to take empirical values ​​to carry out rock physics modeling and analyze the P- and S-wave velocities of a well A. Figure 9b The method described in the embodiment of this specification is used to perform rock physics modeling and analyze the P- and S-wave velocities of a well A, and the attached Figure 9a and attached Figure 9b In the figure, the dashed line corresponds to the rock physics modeling results, and the solid line corresponds to the actual logging data. Figure 9a And attached Figure 9b The horizontal axis units are all stratum depth, and the vertical axis units are located at the top of the figure. From left to right, they are density, mud content, porosity, saturation, P-wave velocity and S-wave velocity. Figure 9a and attached Figure 9b The error between the actual logging data and the rock physics modeling results can determine the Figure 9a There is a large error between the P-wave and S-wave velocities calculated from the empirical values ​​and the actual logging data, and the error is particularly obvious at the formation depth of 2830m-2900m. The rock physics modeling accuracy is low, and the Figure 9b It is shown that the accuracy of rock physics modeling obtained using the method described in the embodiments of this specification is significantly improved.

[0192] It should be noted that, although the operations of the method of the present invention are described in a specific order in the above embodiments and the accompanying drawings, this does not require or imply that the operations must be performed in the specific order, or that all the operations shown must be performed to achieve the desired results. Additionally or alternatively, some steps may be omitted, multiple steps may be combined into one step, and / or one step may be decomposed into multiple steps.

[0193] Corresponding to the above-mentioned turbidite sandstone reservoir seismic rock physical modeling method, some embodiments of this specification also provide a turbidite sandstone reservoir seismic rock physical modeling device, referring to Fig.10 As shown, in some embodiments, the apparatus may include:

[0194] Receiving module 1001, used to obtain target turbidite sandstone reservoir information;

[0195] A rock matrix determination module 1002 is used to determine a target rock matrix according to the target turbidite sandstone reservoir information;

[0196] A dry rock skeleton determination module 1003 is used to establish a dry rock skeleton model using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir;

[0197] The modeling module 1004 is used to establish a fluid-saturated rock model based on the dry rock skeleton model and the target rock matrix.

[0198] For the convenience of description, the above device is described in various units according to their functions. Of course, when implementing this specification, the functions of each unit can be implemented in the same or multiple software and / or hardware.

[0199] It should be noted that in the embodiments of this specification, the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved are all information and data authorized by the user and fully authorized by all parties.

[0200] The embodiments of this specification also provide a computer device. Fig.11 As shown, in some embodiments of the present specification, the computer device 1102 may include one or more processors 1104, such as one or more central processing units (CPUs) or graphics processing units (GPUs), and each processing unit may implement one or more hardware threads. The computer device 1102 may also include any memory 1106, which is used to store any kind of information such as code, settings, data, etc. In a specific embodiment, the computer program on the memory 1106 and can be run on the processor 1104, when the computer program is run by the processor 1104, the instructions of the method described in any of the above embodiments may be executed. Non-limitingly, for example, the memory 1106 may include any one or more combinations of the following: any type of RAM, any type of ROM, flash memory device, hard disk, optical disk, etc. More generally, any memory may use any technology to store information. Further, any memory may provide volatile or non-volatile retention of information. Further, any memory may represent a fixed or removable component of the computer device 1102. In one embodiment, when the processor 1104 executes the associated instructions stored in any memory or combination of memories, the computer device 1102 can perform any operation of the associated instructions. The computer device 1102 also includes one or more drive mechanisms 1108 for interacting with any memory, such as a hard disk drive mechanism, an optical disk drive mechanism, etc.

[0201] The computer device 1102 may also include an input / output interface 1110 (I / O) for receiving various inputs (via input devices 1112) and for providing various outputs (via output devices 1114). A specific output mechanism may include a presentation device 1116 and an associated graphical user interface 1118 (GUI). In other embodiments, the input / output interface 1110 (I / O), input device 1112, and output device 1114 may not be included, and the computer device 1102 may be used as a computer device in a network. The computer device 1102 may also include one or more network interfaces 1120 for exchanging data with other devices via one or more communication links 1122. One or more communication buses 1124 couple the components described above together.

[0202] The communication link 1122 may be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication link 1122 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc. governed by any protocol or combination of protocols.

[0203] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), computer-readable storage media, and computer program products of some embodiments of the present specification. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processor to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processor generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0204] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processor to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0205] These computer program instructions can also be loaded onto a computer or other programmable data processor so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for executing on the computer or other programmable device to implement the process. Figure 1A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0206] In a typical configuration, a computer device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0207] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0208] Computer-readable media include permanent and non-permanent, removable and non-removable media that can be used to store information by any method or technology. Information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computer device. As defined in this specification, computer-readable media does not include temporary computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0209] Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems or computer program products. Therefore, the embodiments of this specification may take the form of complete hardware embodiments, complete software embodiments or embodiments combining software and hardware. Moreover, the embodiments of this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0210] The present specification embodiments may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The present specification embodiments may also be practiced in distributed computing environments where tasks are performed by remote processors connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

[0211] It should also be understood that in the embodiments of this specification, the term "and / or" is only a description of the association relationship of the associated objects, indicating that three relationships can exist. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship.

[0212] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0213] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the embodiments of this specification. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0214] The above is only an embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims

1. A method for seismic rock physics modeling of turbidite sandstone reservoirs, characterized in that: The method comprises: Obtain target turbidite sandstone reservoir information; Determining a target rock matrix according to the target turbidite sandstone reservoir information; A dry rock skeleton model is established using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir; A fluid-saturated rock model is established based on the dry rock skeleton model and the target rock matrix.

2. The method according to claim 1, characterized in that: The target turbidite sandstone reservoir information at least includes well logging curves and rock matrix mud and sand component modulus information.

3. The method according to claim 2, characterized in that Determining a target rock matrix according to the target turbidite sandstone reservoir information includes: Determining the shale content of the target turbidite sandstone reservoir by using the well logging curve; According to the rock matrix mud and sand component modulus information and the mud content, the rock matrix modulus interval information is determined; wherein the rock matrix mud and sand component modulus information includes the bulk modulus of the sandy mineral component constituting the rock matrix, the bulk modulus of the muddy mineral component, the shear modulus of the sandy mineral component and the shear modulus of the muddy mineral component; the rock matrix modulus interval information includes the upper limit and lower limit of the rock matrix bulk modulus, and the upper limit and lower limit of the rock matrix shear modulus; The bulk modulus and shear modulus of the target rock matrix are determined according to the rock matrix modulus interval information.

4. The method according to claim 1, characterized in that The dry rock skeleton model is established by using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir, including: Determine the bulk modulus and shear modulus of the cemented rock using the average coordination number and critical porosity of the sandy mineral particles; The dry rock skeleton model is established according to the bulk modulus and shear modulus of the cemented rock.

5. The method according to claim 4, characterized in that The bulk modulus and shear modulus of the cemented rock are determined using the average coordination number and critical porosity of the sandy mineral particles according to the following formula: in, in, Among them, K b is the bulk modulus of cemented rock, C is the average coordination number of sandy mineral particles, φ c is the critical porosity of the rock in its initial sedimentary state, M shale is the longitudinal wave modulus of mud minerals in rocks, ρ shale is the density of argillaceous minerals, V Pshale is the longitudinal wave velocity of argillaceous minerals, S n is the normal flexibility of the cemented sandy mineral particles, G b is the shear modulus of cemented rock, G shale is the shear modulus of argillaceous minerals, V Sshale is the shear wave velocity of argillaceous minerals, S τ It is the longitudinal flexibility of cemented sandy mineral particles.

6. The method according to claim 5, characterized in that The normal and longitudinal flexibility of cemented sandy mineral particles are determined using the following formula: S n =A n (L n )a 2 +B n (L n )a+C n (L n ) in, S τ =A τ (L τ ,v)a 2 +B τ (L τ ,v)α+C τ (L τ ,v) in, Among them, A n ,Λ n , B n , C n , A τ ,Λ τ , B τ , C τ are all intermediate variables, α represents the cementation mode of sandy mineral particles, π is the circumference, G sand is the shear modulus of the matrix mineral, v sand represents the Poisson's ratio of sandy minerals, v shale represents the Poisson's ratio of mud minerals, and v represents the Poisson's ratio taking into account the bulk modulus and shear modulus of cemented rock.

7. The method according to claim 6, characterized in that The cementation mode of sandy mineral particles is determined according to the average coordination number of the sandy mineral particles and the effective reservoir pressure using the following formula: C=0.5*(777.1P c -0.001545 -770.7)+0.5*(-4.457P c -0.2724 +9.401) Among them, P c is the effective reservoir pressure, S is the percentage of pore space occupied by cement, and φ is the porosity of fluid-saturated rock.

8. The method according to claim 7, characterized in that According to the bulk modulus and shear modulus of the cemented rock, the dry rock skeleton model is established using the following formula: Among them, K dry is the bulk modulus of the dry rock skeleton, φ b is the cementation porosity, K ma is the bulk modulus of the target rock matrix, G ma is the shear modulus of the target rock matrix, G dry is the shear modulus of the dry rock skeleton, and z is an intermediate quantity.

9. The method according to claim 1, characterized in that: A fluid saturated rock model is established based on the dry rock skeleton model and the target rock matrix, including: Determining a mixed fluid corresponding to the fluid-saturated rock model according to the effects of pressure and temperature changes on the fluid; Determine the bulk modulus and shear modulus of the fluid-saturated rock based on the bulk modulus of the mixed fluid, the dry rock skeleton model, and the target rock matrix; According to the bulk modulus and shear modulus of fluid-saturated rock and the density of the mixed fluid, the longitudinal wave velocity and shear wave velocity of the equivalent fluid-saturated rock after the rock physical fluid replacement are determined to establish a fluid-saturated rock model.

10. The method according to claim 9, characterized in that Determining a mixed fluid corresponding to the fluid saturated rock model according to the influence of pressure and temperature changes on the fluid includes: Determine the flow rate and density of oil and the flow rate and density of salt water based on the effects of pressure and temperature changes on the fluid; Based on the flow rate and density of the oil, and the flow rate and density of the brine, the bulk modulus and density of the mixed fluid consisting of the oil and the brine are determined.

11. The method according to claim 10, characterized in that Based on the effects of pressure and temperature changes on the fluid, the flow rate and density of oil, as well as the flow rate and density of brine, are determined using the following formula: r Oil =[ρ0+(0.00277P-1.71×10 -7 P 3 (ρ0-1.15) 2 +3.49×10 -4 P] / [0.972+3.81×10 -4 (T+17.78) 1.175 ] Among them, V Oil is the oil flow rate considering the changes in temperature and pressure, ρ0 is the oil density at normal temperature and pressure, T is the temperature, P is the pressure, ρ Oil Oil density to account for temperature and pressure changes; V brine =V water +S water 1.5 (780-10P+0.16P2)-820S water 2 +S water (1170-9.6T+0.055T2-8.5×10-5T3+2.6P-0.0029TP-0.0476P2) ρ water =1+1×10 -6 (-80T-3.3T 2 +0.00175T 3 +489P-2TP+0.016T 2 P-1.3×10 -5 T 3 P-0.333P 2 -0.002TP 2 ) r brine =ρ water +0.668S water +0.44S water 2 +10 -6 S water [300P-2400PS water +T(80+3T-3300S water -13P+47PS water )] Among them, V brine V is the brine flow rate considering the changes in temperature and pressure. water is the pure water speed, S water is the salt concentration in the brine, ρ water is the density of pure water considering temperature and pressure changes, ρ brine is the density of salt water taking into account changes in temperature and pressure.

12. The method according to claim 10, characterized in that Based on the flow rate and density of the oil, and the flow rate and density of the brine, the bulk modulus and density of the mixed fluid consisting of oil and brine are calculated using the following formula: r fl =S brine r brine +S Oil r Oil Among them, K fl is the bulk modulus of the mixed fluid in the pore, S brine is the brine saturation, ρ brine V is the density of salt water considering the changes in temperature and pressure. brine S is the brine flow rate considering the changes in temperature and pressure. Oil is the oil saturation, ρ Oil V is the oil density considering the changes in temperature and pressure. Oil is the oil flow rate considering temperature and pressure changes, ρ fl is the density of the mixed fluid in the pores.

13. The method according to claim 9, characterized in that Based on the bulk modulus of the mixed fluid, the dry rock skeleton model and the target rock matrix, the bulk modulus and shear modulus of the fluid-saturated rock are determined using the following formula: G sat =G dry Among them, K sat is the bulk modulus of fluid-saturated rock, K dry is the bulk modulus of the dry rock skeleton, K ma is the bulk modulus of the target rock matrix, φ is the porosity of the fluid-saturated rock, K fl is the bulk modulus of the mixed fluid in the pore, G dry is the shear modulus of the dry rock skeleton, G sat is the shear modulus of fluid-saturated rock.

14. A seismic rock physics modeling device for turbidite sandstone reservoirs, characterized in that: The device comprises: A receiving module, used for acquiring target turbidite sandstone reservoir information; A rock matrix determination module, used to determine a target rock matrix according to the target turbidite sandstone reservoir information; A dry rock skeleton determination module is used to establish a dry rock skeleton model using the average coordination number and critical porosity of sandy mineral particles corresponding to the target turbidite sandstone reservoir; A modeling module is used to establish a fluid-saturated rock model based on the dry rock skeleton model and the target rock matrix.

15. A computer device comprising a memory, a processor, and a computer program stored in the memory, characterized in that: When the computer program is executed by the processor, the computer program executes the instructions of the method according to any one of claims 1 to 13.

16. A computer storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor of a computer device, the computer program executes instructions of the method according to any one of claims 1 to 13.

17. A computer program product, characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 13 is implemented.

Images