A hydrate saturation calculation method and device, electronic equipment and storage medium

Abstract

The application discloses a hydrate saturation calculation method and device, electronic equipment and a storage medium. The method comprises the following steps: determining the friction coefficient of the skeleton particle roughness according to the logging observation data; calculating the longitudinal wave velocity of the hydrate-bearing sediment and the transverse wave velocity of the hydrate-bearing sediment according to the friction coefficient of the skeleton particle roughness and the rock physical model of the rough particle model; updating the hydrate saturation based on the longitudinal wave velocity of the hydrate-bearing sediment and the transverse wave velocity of the hydrate-bearing sediment to obtain optimal hydrate saturation data. The technical scheme can improve the accuracy of hydrate saturation calculation.

Description

Technical Field

[0001] This invention relates to the field of hydrate saturation calculation technology, and in particular to a method, apparatus, electronic device, and storage medium for calculating hydrate saturation. Background Technology

[0002] Natural gas hydrates are cage-like crystalline compounds formed from natural gas and water under high pressure and low temperature conditions. They are widely distributed in sediments of terrestrial permafrost zones and deep-water basins on continental margins. Due to their high energy density and large reserves, they are considered a potential alternative energy source. However, their decomposition can induce marine geological disasters. Therefore, predicting hydrate content has significant scientific and social value.

[0003] The presence of hydrates in sediments typically increases the elastic wave velocity of the formation. While the linear intersection of velocity and saturation can be used to predict hydrate saturation, this method is simple and easy to use, but its accuracy is low and its error is large. Furthermore, the lithology, porosity, and effective pressure of sediments vary under different geological backgrounds, leading to changes in the linear relationship between velocity and saturation, making it difficult to apply uniformly. Since hydrate saturation is a key parameter for calculating resource quantity, a more applicable and accurate saturation prediction method is needed. Rock physics is a relatively applicable and accurate method for predicting hydrate saturation. The principle of rock physics for estimating hydrate saturation is as follows: Based on initial assumptions about mineral composition, porosity, and saturation parameters, the elastic modulus of dry rock and pore fluid is calculated. Then, the elastic modulus of the saturated fluid is calculated using the Gassmann equation, and velocity data is predicted. By comparing this with measured velocities and continuously adjusting the parameters to achieve the best fit with the measured velocities, parameters such as saturation can be predicted. This method, by varying lithological combinations and porosity parameters, is applicable to different geological backgrounds, and the best fit with observational data also improves the accuracy of the prediction.

[0004] Current rock physics modeling of hydrates relies on particle contact models that assume the framework particles are completely rough and without slippage. In reality, slippage does occur at the edges of particle contact regions. Therefore, the shear modulus predicted by this model is too high, and the predicted velocity is much higher than the actual observed velocity, leading to inaccurate predictions of hydrate saturation. Summary of the Invention

[0005] This invention provides a method, apparatus, electronic device, and storage medium for calculating hydrate saturation, which can improve the accuracy of hydrate saturation calculation.

[0006] According to one aspect of the present invention, a method for calculating hydrate saturation is provided, the method comprising:

[0007] The friction coefficient of the skeleton particle roughness was determined based on well logging data;

[0008] Based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model, the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediments are calculated.

[0009] Based on the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediments, the hydrate saturation is updated to obtain the optimal hydrate saturation data.

[0010] According to another aspect of the present invention, a hydrate saturation calculation apparatus is provided, the apparatus comprising:

[0011] The friction coefficient determination module is used to determine the friction coefficient of the skeleton particle roughness based on well logging observation data;

[0012] The velocity calculation module is used to calculate the longitudinal wave velocity and transverse wave velocity of hydrate-bearing sediments based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model.

[0013] The optimal hydrate saturation data acquisition module is used to update the hydrate saturation based on the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment to obtain the optimal hydrate saturation data.

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

[0015] At least one processor; and

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

[0017] The memory stores a computer program that can be executed by the at least one processor, which enables the at least one processor to perform a hydrate saturation calculation method according to any embodiment of the present invention.

[0018] According to another aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing computer instructions for causing a processor to execute and implement a hydrate saturation calculation method according to any embodiment of the present invention.

[0019] The technical solution of this invention determines the friction coefficient of the skeleton particle roughness based on well logging data. Then, based on the friction coefficient of the skeleton particle roughness and the rock physics model of the rough particle model, it calculates the P-wave velocity and S-wave velocity of the hydrate-bearing sediment. Based on these P-wave and S-wave velocities, the hydrate saturation is updated to obtain optimal hydrate saturation data. This technical solution improves the accuracy of hydrate saturation calculation.

[0020] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0022] Figure 1 This is a flowchart of a hydrate saturation calculation method provided in Embodiment 1 of the present invention;

[0023] Figure 2 This is a schematic diagram of a hydrate saturation calculation process provided in Embodiment 2 of the present invention;

[0024] Figure 3 This is a flowchart of the hydrate saturation treatment process provided in Embodiment 2 of this application;

[0025] Figure 4 This is a schematic diagram of another hydrate saturation calculation process provided in Embodiment 3 of the present invention;

[0026] Figure 5 (a) is a schematic diagram of the friction coefficient provided in Embodiment 3 of this application;

[0027] Figure 5 (b) is a comparison chart of the longitudinal wave velocity calculation results provided in Embodiment 3 of this application;

[0028] Figure 5 (c) is a comparison chart of the shear wave velocity calculation results provided in Embodiment 3 of this application;

[0029] Figure 6 (a) is a schematic diagram of the hydrate saturation calculated by the method provided in Embodiment 3 of this application;

[0030] Figure 6(b) is a schematic diagram of hydrate saturation calculated by the conventional method provided in Embodiment 3 of this application;

[0031] Figure 7 This is a schematic diagram of the structure of a hydrate saturation calculation device provided in Embodiment 4 of the present invention;

[0032] Figure 8 This is a schematic diagram of the structure of an electronic device that implements a hydrate saturation calculation method according to an embodiment of the present invention. Detailed Implementation

[0033] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0034] It should be noted that the terms "initial," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0035] Example 1

[0036] Figure 1 This is a flowchart illustrating a method for calculating hydrate saturation according to Embodiment 1 of the present invention. This embodiment is applicable to situations where hydrate saturation is calculated based on a coarse particle model. The method can be executed by a hydrate saturation calculation device, which can be implemented in hardware and / or software and can be configured within a device. For example, the device can be a backend server or other equipment with communication and computing capabilities. Figure 1 As shown, the method includes:

[0037] S110. Determine the friction coefficient of the skeleton particle roughness based on well logging observation data.

[0038] In this scheme, well logging data includes P-wave velocity, S-wave velocity, density data, initial hydrate saturation, structural composition, and mineral content. Well logging data can be acquired through well logging.

[0039] In this embodiment, the friction coefficient is used to characterize the roughness of the skeleton particles. By determining the friction coefficient, the accuracy of hydrate saturation calculation can be improved.

[0040] Specifically, the friction coefficient of the skeleton particles roughness can be calculated based on the P-wave velocity, S-wave velocity, density data, initial hydrate saturation, structural components and mineral content in the well logging observation data, according to a predetermined calculation formula.

[0041] S120. Calculate the longitudinal wave velocity and transverse wave velocity of hydrate-bearing sediments based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model.

[0042] In this scheme, the friction coefficient of the skeleton particle roughness can be calculated through well logging observation data, and the friction coefficient considering the particle roughness is introduced into the rock physics model of the rough particle model to calculate the P-wave velocity and S-wave velocity of hydrate-bearing sediments, which is more applicable than the traditional rock physics model.

[0043] Furthermore, by varying the friction coefficient, different roughnesses of the skeletal particles can be represented, making it applicable to various geological conditions. The friction coefficient ranges from 0 to 1. When the friction coefficient is 1, it represents complete roughness, and the rock physics model of the rough particle model is equivalent to the traditional model. When the friction coefficient is 0, it represents complete smoothness, which is applicable to loose and unconsolidated strata.

[0044] S130. Based on the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment, the hydrate saturation is updated to obtain the optimal hydrate saturation data.

[0045] In this scheme, the hydrate saturation can be updated using the least squares method to minimize the error between the observed values ​​and the P-wave velocity and S-wave velocity of the hydrate-bearing sediments, thereby obtaining the optimal hydrate saturation data.

[0046] The technical solution of this invention determines the friction coefficient of the skeleton particle roughness based on well logging observation data. Then, based on the friction coefficient of the skeleton particle roughness and the rock physics model of the rough particle model, it calculates the P-wave velocity and S-wave velocity of the hydrate-bearing sediment. Based on these P-wave and S-wave velocities, the hydrate saturation is updated to obtain optimal hydrate saturation data. This technical solution is highly practical, the established rock physics model has wider applicability, can be used to characterize hydrate rock physics models with different skeleton particle roughness, and can improve the accuracy of hydrate saturation calculation.

[0047] Example 2

[0048] Figure 2 This is a schematic diagram illustrating the hydrate saturation calculation process provided in Embodiment 2 of the present invention. The relationship between this embodiment and the above embodiments is illustrated in the detailed description of the process for determining the friction coefficient of the skeleton particle roughness. For example... Figure 2 As shown, the method includes:

[0049] S210. Calculate the bulk modulus and shear modulus of the saturated fluid based on the P-wave velocity, S-wave velocity, and density data from the well logging observation data.

[0050] In this embodiment, Figure 3 This is a flowchart of the hydrate saturation treatment process provided in Embodiment 2 of this application, as follows: Figure 3 As shown, it can be based on the longitudinal wave velocity V p Shear wave velocity V s Calculate the bulk modulus K of the saturated fluid using the density data ρ. sat and shear modulus G sat .

[0051] In this scheme, the bulk modulus and shear modulus of the saturated fluid are calculated using the following formulas;

[0052]

[0053] Among them, K sat G is the bulk modulus of a saturated fluid. sat Here, ρ represents the shear modulus of the saturated fluid, ρ represents the density data, and V represents the shear modulus of the fluid. p V is the longitudinal wave velocity. s This represents the transverse wave velocity.

[0054] S220. Based on the initial hydrate saturation, structural components, and mineral content in the well logging observation data, calculate the matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus.

[0055] In this plan, such as Figure 3As shown, the initial hydrate saturation S in the well logging observation data can be obtained. h Structural components (f) rm f pf ) and parameters such as mineral content. Among them, f rm It is the volume fraction of hydrates supported by the skeleton, f pf It is the volume fraction of hydrates filling the pores, f rm +f pf =1,

[0056] Furthermore, such as Figure 3 As shown, it can be based on the initial hydrate saturation S h Structural components (f) rm f pf ) and mineral content and other parameters to calculate matrix bulk modulus K, matrix shear modulus G and pore fluid bulk modulus K f .

[0057] In this embodiment, the matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus are calculated using the following formulas;

[0058]

[0059] Where φ is porosity, S h It is the initial hydrate saturation, f rm It is the volume fraction of hydrates supported by the skeleton, f pf It is the volume fraction of hydrates filling the pores, f rm +f pf =1,K h It is the bulk modulus of hydrates, K i f is the bulk modulus of the i-th mineral. i G is the volume fraction of the i-th mineral. h It is the shear modulus of hydrate, G i K is the shear modulus of the i-th mineral, K is the matrix bulk modulus, and G is the matrix shear modulus. f This represents the bulk modulus of the pore fluid.

[0060] S230. Calculate the bulk modulus and shear modulus of dry rock using the bulk modulus and shear modulus of the saturated fluid, the bulk modulus of the matrix, the shear modulus of the matrix, and the bulk modulus of the pore fluid.

[0061] In this embodiment, as Figure 3 As shown, it can be based on the bulk modulus K of the saturated fluid. sat and shear modulus G sat Matrix bulk modulus K, matrix shear modulus G, and pore fluid bulk modulus K f Calculate the bulk modulus of dry rock dry rock shear modulus

[0062] In this scheme, the bulk modulus and shear modulus of dry rock are calculated using the following formulas;

[0063]

[0064] in, The bulk modulus of dry rock. This is the dry rock shear modulus.

[0065] S240. Based on the bulk modulus and shear modulus of the dry rock, calculate the friction coefficient of the skeleton particles roughness.

[0066] In this plan, such as Figure 3 As shown, the calculated bulk modulus of dry rock can be used as a basis. dry rock shear modulus Calculate the friction coefficient α of the skeleton particle roughness. Specifically, this needs to be based on the dry rock bulk modulus. dry rock shear modulus The matrix bulk modulus K and matrix shear modulus G are used to calculate the friction coefficient α of the skeleton particle roughness.

[0067] Furthermore, the friction coefficient of the skeleton particle roughness is calculated using the following formula;

[0068]

[0069] Where α is the friction coefficient of the skeleton particles, and v is the Poisson's ratio of the skeleton particles.

[0070] S250. Based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model, calculate the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment.

[0071] S260. Based on the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment, the hydrate saturation is updated to obtain the optimal hydrate saturation data.

[0072] The technical solution of this invention calculates the bulk modulus and shear modulus of saturated fluid based on P-wave velocity, S-wave velocity, and density data from well logging observations. It also calculates the matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus based on the initial hydrate saturation, structural composition, and mineral content from the well logging observations. Then, using the bulk modulus and shear modulus of saturated fluid, matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus, the dry rock bulk modulus and dry rock shear modulus are calculated. Based on these dry rock bulk modulus and dry rock shear modulus, the friction coefficient of the framework particle roughness is calculated. Then, based on the friction coefficient of the framework particle roughness and the rock physics model of the rough particle model, the P-wave velocity and S-wave velocity of the hydrate-bearing sediment are calculated. Based on these P-wave velocity and S-wave velocity, the hydrate saturation is updated to obtain optimal hydrate saturation data. By implementing this technical solution, a friction coefficient that takes into account particle roughness is introduced into the rock physics modeling of hydrates. This makes it more applicable than traditional rock physics models and can improve the accuracy of hydrate saturation calculation.

[0073] Example 3

[0074] Figure 4 This is a schematic diagram illustrating another hydrate saturation calculation process provided in Embodiment 3 of the present invention. The relationship between this embodiment and the above embodiments is described in detail regarding the hydrate saturation calculation process. Figure 4 As shown, the method includes:

[0075] S410. Determine the friction coefficient of the skeleton particle roughness based on well logging observation data.

[0076] S420, Constructing a rock physics model with coarse grains.

[0077] In this scheme, the rock physics model for constructing the coarse particle model is constructed using the following formula;

[0078]

[0079] Where, φ c For critical porosity, K HM G is the bulk modulus of dry rock at critical porosity. HM Where is the shear modulus of critical porosity, C is the coordination number of the framework particles (valued at 8.5), and P is the effective pressure, P = (ρ - ρ w )·g·h, where ρ is the formation density, ρ w ρ is the density of seawater, g is the acceleration due to gravity, and h is the depth below the seabed.

[0080] Furthermore, when α = 1, it represents a completely rough model; when α = 0, it represents a completely smooth model. The introduction of a friction coefficient that considers particle roughness makes the rock physics model with rough particles more applicable than the traditional rock physics model.

[0081] S430. The friction coefficient of the roughness of the skeleton particles is substituted into the rock physics model of the rough particle model to calculate the dry rock bulk modulus and dry rock shear modulus of the critical porosity.

[0082] In this scheme, the friction coefficient α of the calculated skeleton particle roughness can be substituted into the rock physics model of the rough particle model to calculate the critical porosity φ. c Dry rock bulk modulus K HM Dry rock shear modulus G HM Among them, the critical porosity φ c For unconsolidated sandstone, the value is 0.4.

[0083] S440. Based on the dry rock bulk modulus and dry rock shear modulus of the critical porosity, calculate the dry rock bulk modulus and dry rock shear modulus of any porosity.

[0084] In this scheme, the modified Hashin-Shtrikman limit is applied to calculate the dry rock bulk modulus K for arbitrary porosity. dry and shear modulus G dry .

[0085] Specifically, the bulk modulus and shear modulus of dry rock with arbitrary porosity are calculated using the following formulas;

[0086]

[0087]

[0088] Among them, K dry G is the bulk modulus of dry rock with arbitrary porosity. dry For dry rock with arbitrary porosity, the shear modulus is given.

[0089] S450. Based on the dry rock bulk modulus and dry rock shear modulus of the arbitrary porosity, calculate the bulk modulus and shear modulus of the initial parameter saturated fluid.

[0090] In this embodiment, the bulk modulus of the saturated fluid with initial parameters can be obtained by applying the Gassmann equation. and shear modulus

[0091] Specifically, the bulk modulus and shear modulus of the saturated fluid with initial parameters are calculated using the following formulas;

[0092]

[0093] in, The initial parameter is the bulk modulus of the saturated fluid. The initial parameter is the shear modulus of the saturated fluid.

[0094] S460. Calculate the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment based on the initial parameters, the bulk modulus and shear modulus of the saturated fluid.

[0095] Furthermore, the P-wave velocity and S-wave velocity of hydrate-bearing sediments are calculated using the following formulas;

[0096]

[0097] in, For the longitudinal wave velocity and V of hydrate-bearing sediments est s ρ represents the shear wave velocity of hydrate-bearing sediments. b est For density, ρ i Let ρ be the density of the i-th mineral in the skeleton. h The density of the pure hydrate is taken as 0.9 g / cm³. 3 , ρ w The density of seawater is taken as 1.032 g / cm³. 3 .

[0098] S470. Construct an optimal function based on the longitudinal wave velocity of the hydrate-bearing sediment, the transverse wave velocity of the hydrate-bearing sediment, the longitudinal wave velocity, the transverse wave velocity, the friction coefficient of the skeleton particle roughness, the hydrate saturation, the volume fraction of hydrate supporting the skeleton, and the volume fraction of hydrate filling the pores.

[0099] In this scheme, it can be seen from the calculation formulas for the P-wave and S-wave velocities of hydrate-bearing sediments that the P-wave velocity... and transverse wave velocity V est s It is the hydrate saturation S h And a function of the friction coefficient α. Among these, the hydrate saturation can be divided into the skeletal support portion f. rm and pore filling part f pf Therefore, given a set (S) h f rm f pfThe parameters α) can be used to obtain estimates of the P-wave and S-wave velocities. To determine hydrate saturation from actual observed velocity data, the optimal hydrate saturation can be calculated by minimizing the residual between the estimated and observed values, achieving the best fit between the estimated and observed velocity data. Therefore, an optimization function is established based on the P-wave velocity of hydrate-bearing sediments, the S-wave velocity of hydrate-bearing sediments, the P-wave velocity, the S-wave velocity, the friction coefficient of the framework particle roughness, the hydrate saturation, the volume fraction of hydrates supported by the framework, and the volume fraction of hydrates filling the pores.

[0100] minfun(S h ,f rm ,f pf ,α)=||V p est -V p V s est -V s ||.

[0101] S480. Based on the optimal function, update the hydrate saturation to obtain the optimal hydrate saturation data.

[0102] Furthermore, the initial value can be updated using an iterative least squares method, while limiting (S) h f rm f pf , α) parameters satisfy f rm +f pf =1 and S h Given the constraint α∈[0,1], the optimal saturation and the content of hydrates in different component structures can be obtained when the optimal function value is minimized, thus obtaining the optimal hydrate saturation data.

[0103] The technical solution of this invention determines the friction coefficient of the skeleton particle roughness based on well logging observation data. A rock physics model of the rough particle model is constructed; the friction coefficient of the skeleton particle roughness is substituted into the rock physics model of the rough particle model to calculate the dry rock bulk modulus and dry rock shear modulus for critical porosity; based on the dry rock bulk modulus and dry rock shear modulus for critical porosity, the dry rock bulk modulus and dry rock shear modulus for arbitrary porosity are calculated; based on the dry rock bulk modulus and dry rock shear modulus for arbitrary porosity, the bulk modulus and shear modulus of the initial parameter saturated fluid are calculated; based on the bulk modulus and shear modulus of the initial parameter saturated fluid, the P-wave velocity and S-wave velocity of the hydrate-bearing sediment are calculated. An optimal function is constructed based on the P-wave velocity, S-wave velocity, P-wave velocity, S-wave velocity, friction coefficient of framework particle roughness, hydrate saturation, volume fraction of framework-supported hydrate, and volume fraction of pore-filled hydrate. Based on this optimal function, the hydrate saturation is updated to obtain optimal hydrate saturation data. This technical solution demonstrates strong practicality, and the established rock physics model has wider applicability. It can be used to characterize rock physics models with different roughness levels of framework particles. The calculated P-wave velocity and S-wave velocity are in better agreement with measured data, and the inverted saturation parameters are more accurate, thus improving the accuracy of hydrate saturation calculation.

[0104] In this plan, Figure 5 (a) is a schematic diagram of the friction coefficient provided in Embodiment 3 of this application. Figure 5 (b) is a comparison chart of the P-wave velocity calculation results provided in Embodiment 3 of this application. Figure 5 (c) is a comparison chart of the shear wave velocity calculation results provided in Embodiment 3 of this application, such as Figure 5 As shown in (a), (b), and (c), the P-wave and S-wave velocities predicted by conventional methods are higher than the actual measured values, especially the S-wave velocity, which is significantly higher than the measured value. The main reason is that the slippage between skeletal particles is not considered, resulting in an overestimation of the calculated shear modulus. Based on the calculation method of this scheme, the friction coefficient is calculated first, and then the friction coefficient is introduced into the rock physics model, taking into account the slippage between particles. The calculated P-wave and S-wave velocities are in better agreement with the actual measured values.

[0105] In this embodiment, Figure 6 (a) is a schematic diagram of the hydrate saturation calculated by the method provided in Embodiment 3 of this application. Figure 6 (b) is a schematic diagram of the hydrate saturation calculated by the conventional method provided in Embodiment 3 of this application, as shown in the figure. Figure 6As shown in (a) and 6(b), the hydrate saturation predicted by conventional rock physics methods is significantly lower than that predicted by resistivity, while the hydrate saturation predicted by this method is closer to that predicted by resistivity, indicating that the method is practical and effective.

[0106] Example 4

[0107] Figure 7 This is a schematic diagram of a hydrate saturation calculation device provided in Embodiment 4 of the present invention. Figure 7 As shown, the device includes:

[0108] Friction coefficient determination module 710 is used to determine the friction coefficient of the skeleton particle roughness based on well logging observation data;

[0109] The velocity calculation module 720 is used to calculate the longitudinal wave velocity and transverse wave velocity of hydrate-bearing sediments based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model.

[0110] The optimal hydrate saturation data acquisition module 730 is used to update the hydrate saturation based on the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment to obtain the optimal hydrate saturation data.

[0111] Optional, the friction coefficient determination module 710 is specifically used for:

[0112] Based on the P-wave velocity, S-wave velocity, and density data from the well logging observation data, calculate the bulk modulus and shear modulus of the saturated fluid;

[0113] Based on the initial hydrate saturation, structural components, and mineral content in the well logging data, the matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus are calculated.

[0114] Using the bulk modulus and shear modulus of the saturated fluid, the bulk modulus of the matrix, the shear modulus of the matrix, and the bulk modulus of the pore fluid, calculate the bulk modulus of dry rock and the shear modulus of dry rock;

[0115] The friction coefficient of the skeleton particles roughness is calculated based on the dry rock bulk modulus and dry rock shear modulus.

[0116] Optionally, the friction coefficient determination module 710 is also used for:

[0117] The bulk modulus and shear modulus of a saturated fluid are calculated using the following formulas;

[0118]

[0119] Among them, K sat G is the bulk modulus of a saturated fluid.sat Here, ρ represents the shear modulus of the saturated fluid, ρ represents the density data, and V represents the shear modulus of the fluid. p V is the longitudinal wave velocity. s The transverse wave velocity;

[0120] The matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus are calculated using the following formulas;

[0121]

[0122] Where φ is porosity, S h It is the initial hydrate saturation, f rm It is the volume fraction of hydrates supported by the skeleton, f pf It is the volume fraction of hydrates filling the pores, f rm +f pf =1,K h It is the bulk modulus of hydrate, K i f is the bulk modulus of the i-th mineral. i G is the volume fraction of the i-th mineral. h It is the shear modulus of hydrate, G i K is the shear modulus of the i-th mineral, K is the matrix bulk modulus, and G is the matrix shear modulus. f The bulk modulus of pore fluid;

[0123] The bulk modulus and shear modulus of dry rock are calculated using the following formulas;

[0124]

[0125] in, The bulk modulus of dry rock. Shear modulus of dry rock;

[0126] The friction coefficient of the skeleton particle roughness is calculated using the following formula;

[0127]

[0128] Where α is the friction coefficient of the skeleton particles, and v is the Poisson's ratio of the skeleton particles.

[0129] Optional, speed calculation module 720, specifically used for:

[0130] Constructing a rock physics model based on coarse grains;

[0131] The friction coefficient of the roughness of the skeleton particles is substituted into the rock physics model of the rough particle model to calculate the dry rock bulk modulus and dry rock shear modulus of critical porosity.

[0132] Based on the dry rock bulk modulus and dry rock shear modulus of the critical porosity, calculate the dry rock bulk modulus and dry rock shear modulus of any porosity.

[0133] Based on the bulk modulus and shear modulus of dry rock with arbitrary porosity, calculate the bulk modulus and shear modulus of the initial parameter saturated fluid.

[0134] Based on the initial parameters of the bulk modulus and shear modulus of the saturated fluid, the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment are calculated.

[0135] Optionally, the speed calculation module 720 is also used for:

[0136] The following formula is used to construct the rock physics model for the coarse-grained model;

[0137]

[0138] Where, φ c For critical porosity, K HM G is the dry rock bulk modulus at critical porosity. HM Where C is the shear modulus of critical porosity, C is the coordination number of the framework particles, and P is the effective pressure, P = (ρ - ρ w )·g·h, where ρ is the formation density, ρ w ρ is the density of seawater, g is the acceleration due to gravity, and h is the depth below the seabed.

[0139] Optionally, the speed calculation module 720 is also used for:

[0140] The bulk modulus and shear modulus of dry rock with arbitrary porosity are calculated using the following formulas;

[0141]

[0142]

[0143] Among them, K dry G is the bulk modulus of dry rock with arbitrary porosity. dry For dry rock with arbitrary porosity;

[0144] The bulk modulus and shear modulus of the saturated fluid with initial parameters are calculated using the following formulas;

[0145]

[0146] in, The initial parameter is the bulk modulus of the saturated fluid. The initial parameter is the shear modulus of the saturated fluid.

[0147] The P-wave velocity and S-wave velocity of hydrate-bearing sediments are calculated using the following formulas;

[0148]

[0149] in, For the longitudinal wave velocity and V of hydrate-bearing sediments est s ρ represents the shear wave velocity of hydrate-bearing sediments. b est For density, ρ i Let ρ be the density of the i-th mineral in the skeleton. h The density of the pure hydrate is ρ. w This refers to the density of seawater.

[0150] Optionally, the optimal hydrate saturation data acquisition module 730 is specifically used for:

[0151] An optimal function is constructed based on the longitudinal wave velocity of the hydrate-bearing sediment, the transverse wave velocity of the hydrate-bearing sediment, the longitudinal wave velocity, the transverse wave velocity, the friction coefficient of the skeleton particle roughness, the hydrate saturation, the volume fraction of hydrate supported by the skeleton, and the volume fraction of hydrate filling the pores.

[0152] Based on the optimal function, the hydrate saturation is updated to obtain the optimal hydrate saturation data.

[0153] The hydrate saturation calculation device provided in this embodiment of the invention can execute the hydrate saturation calculation method provided in any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the method.

[0154] Example 5

[0155] Figure 8 A schematic diagram of an electronic device 10 that can be used to implement embodiments 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 processors, cellular phones, smartphones, wearable devices (e.g., helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.

[0156] like Figure 8As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 may also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

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

[0158] Processor 11 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 11 performs the various methods and processes described above, such as a hydrate saturation calculation method.

[0159] In some embodiments, a hydrate saturation calculation method may be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the hydrate saturation calculation method described above may be performed. Alternatively, in other embodiments, processor 11 may be configured to perform a hydrate saturation calculation method by any other suitable means (e.g., by means of firmware).

[0160] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0161] Computer programs used to implement the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

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

[0163] 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 provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, 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 sound input, voice input, or tactile input).

[0164] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or computing systems that include middleware components (e.g., application servers), or computing systems that include frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0165] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.

[0166] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0167] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.

Claims

1. A method for calculating hydrate saturation, characterized in that, include: The friction coefficient of the skeleton particle roughness was determined based on well logging data; Based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model, the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediments are calculated. Based on the P-wave velocity and S-wave velocity of the hydrate-bearing sediments, the hydrate saturation is updated to obtain the optimal hydrate saturation data. The friction coefficient of the skeleton particle roughness is determined based on well logging data, including: Based on the P-wave velocity, S-wave velocity, and density data from the well logging observation data, calculate the bulk modulus and shear modulus of the saturated fluid; Based on the initial hydrate saturation, structural components, and mineral content in the well logging data, the matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus are calculated. Using the bulk modulus and shear modulus of the saturated fluid, the bulk modulus of the matrix, the shear modulus of the matrix, and the bulk modulus of the pore fluid, calculate the bulk modulus of dry rock and the shear modulus of dry rock; Based on the dry rock bulk modulus and dry rock shear modulus, the friction coefficient of the skeleton particles roughness is calculated. Specifically, based on the friction coefficient of the skeleton particle roughness and the rock physics model of the rough particle model, the longitudinal wave velocity and transverse wave velocity of the hydrate are calculated, including: Constructing a rock physics model based on coarse grains; The friction coefficient of the roughness of the skeleton particles is substituted into the rock physics model of the rough particle model to calculate the dry rock bulk modulus and dry rock shear modulus of critical porosity. Based on the dry rock bulk modulus and dry rock shear modulus of the critical porosity, calculate the dry rock bulk modulus and dry rock shear modulus of any porosity. Based on the bulk modulus and shear modulus of dry rock with arbitrary porosity, calculate the bulk modulus and shear modulus of the initial parameter saturated fluid. Based on the initial parameters of the bulk modulus and shear modulus of the saturated fluid, the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment are calculated.

2. The method according to claim 1, characterized in that, include: The bulk modulus and shear modulus of a saturated fluid are calculated using the following formulas; ； in, The bulk modulus of a saturated fluid. The shear modulus of the saturated fluid. For density data, For the longitudinal wave velocity, The transverse wave velocity; The matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus are calculated using the following formulas; ； in, Porosity It is the initial hydrate saturation. It is the skeletal support for the hydrate volume fraction. It is the volume fraction of hydrates filling the pores. , It is the bulk modulus of hydrates. It is the first The bulk modulus of a mineral. It is the first Volume fraction of the minerals It is the shear modulus of hydrates. It is the first The shear modulus of a mineral, The bulk modulus of the matrix. For matrix shear modulus, The bulk modulus of pore fluid; The bulk modulus and shear modulus of dry rock are calculated using the following formulas; ； in, The bulk modulus of dry rock. Shear modulus of dry rock; The friction coefficient of the skeleton particle roughness is calculated using the following formula; ； in, The coefficient of friction represents the roughness of the skeletal particles. It is the Poisson's ratio of the skeletal particles.

3. The method according to claim 2, characterized in that, The rock physics model for constructing a coarse-grained model includes: The following formula is used to construct the rock physics model for the coarse-grained model; ； in, Critical porosity The bulk modulus of dry rock at critical porosity. The shear modulus of critical porosity. For the coordination number of skeletal particles, For effective pressure, , It is the density of the formation. It is the density of seawater. It is gravitational acceleration. It refers to the depth below the seabed.

4. The method according to claim 3, characterized in that, include: The bulk modulus and shear modulus of dry rock with arbitrary porosity are calculated using the following formulas; ； ； in, For dry rock with arbitrary porosity, For dry rock with arbitrary porosity; The bulk modulus and shear modulus of the saturated fluid with initial parameters are calculated using the following formulas; ； in, The initial parameter is the bulk modulus of the saturated fluid. The initial parameter is the shear modulus of the saturated fluid. The P-wave velocity and S-wave velocity of hydrate-bearing sediments are calculated using the following formulas; ； in, For the longitudinal wave velocity of hydrate-bearing sediments, The transverse wave velocity of hydrate-bearing sediments. For density, For the first in the skeleton The density of the mineral, The density of pure hydrate, This refers to the density of seawater.

5. The method according to claim 1, characterized in that, Based on the P-wave velocity and S-wave velocity of the hydrate-bearing sediments, the hydrate saturation is updated to obtain optimal hydrate saturation data, including: An optimal function is constructed based on the longitudinal wave velocity of the hydrate-bearing sediment, the transverse wave velocity of the hydrate-bearing sediment, the longitudinal wave velocity, the transverse wave velocity, the friction coefficient of the skeleton particle roughness, the hydrate saturation, the volume fraction of hydrate supported by the skeleton, and the volume fraction of hydrate filling the pores. Based on the optimal function, the hydrate saturation is updated to obtain the optimal hydrate saturation data.

6. A device for calculating hydrate saturation, characterized in that, include: The friction coefficient determination module is used to determine the friction coefficient of the skeleton particle roughness based on well logging observation data; The velocity calculation module is used to calculate the longitudinal wave velocity and transverse wave velocity of hydrate-bearing sediments based on the friction coefficient of the roughness of the skeleton particles and the rock physics model of the rough particle model. The optimal hydrate saturation data acquisition module is used to update the hydrate saturation based on the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment to obtain the optimal hydrate saturation data. The friction coefficient determination module is specifically used for: Based on the P-wave velocity, S-wave velocity, and density data from the well logging observation data, calculate the bulk modulus and shear modulus of the saturated fluid; Based on the initial hydrate saturation, structural components, and mineral content in the well logging data, the matrix bulk modulus, matrix shear modulus, and pore fluid bulk modulus are calculated. Using the bulk modulus and shear modulus of the saturated fluid, the bulk modulus of the matrix, the shear modulus of the matrix, and the bulk modulus of the pore fluid, calculate the bulk modulus of dry rock and the shear modulus of dry rock; Based on the dry rock bulk modulus and dry rock shear modulus, the friction coefficient of the skeleton particles roughness is calculated. The speed calculation module is specifically used for: Constructing a rock physics model based on coarse grains; The friction coefficient of the roughness of the skeleton particles is substituted into the rock physics model of the rough particle model to calculate the dry rock bulk modulus and dry rock shear modulus of critical porosity. Based on the dry rock bulk modulus and dry rock shear modulus of the critical porosity, calculate the dry rock bulk modulus and dry rock shear modulus of any porosity. Based on the bulk modulus and shear modulus of dry rock with arbitrary porosity, calculate the bulk modulus and shear modulus of the initial parameter saturated fluid. Based on the initial parameters of the bulk modulus and shear modulus of the saturated fluid, the longitudinal wave velocity and transverse wave velocity of the hydrate-bearing sediment are calculated.

7. An electronic device, characterized in that, The electronic device includes: At least one processor; and a memory communicatively connected to said at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, the computer program being executed by the at least one processor to enable the at least one processor to perform a hydrate saturation calculation method according to any one of claims 1-5.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions that cause a processor to execute a method for calculating hydrate saturation according to any one of claims 1-5.

Images