Reservoir in-situ key rock mechanical parameter space inversion prediction method and system

Abstract

The invention provides a reservoir in-situ key rock mechanical parameter space inversion prediction method and a reservoir in-situ key rock mechanical parameter space inversion prediction system, and belongs to the field of shale oil and gas geology, and the reservoir in-situ key rock mechanical parameter space inversion prediction method specifically comprises the following steps: carrying out a nanoindentation experiment to obtain an indentation experiment database, and simulating rock indentation experiments under different conditions to obtain a macroscopic indentation numerical model; using a Latin hypercube sampling method to obtain a simulation parameter combination database, carrying out batch numerical experiments, and constructing an inversion data set; according to the inversion data set, training the deep neural network model subjected to Bayesian optimization to obtain a single-point inversion prediction model; further establishing a spatial inversion prediction model according to a Kriging interpolation algorithm; and inputting real-time detection data of the deep rock stratum into the spatial inversion prediction model to obtain real-time underground rock mechanical parameter spatial distribution data of the deep rock mass. The deep reservoir space key rock mechanical parameters can be accurately inverted, and a basis is provided for well spacing, well trajectory design, fracturing scheme optimization and the like.

Description

Technical Field

[0001] The present invention belongs to the field of shale oil and gas geology, and particularly relates to a method and system for spatially inverting and predicting key in-situ rock mechanical parameters of a reservoir. Background Art

[0002] In recent years, the focus of oil and gas exploration and development in China has gradually shifted to the deep part. The accurate inversion of key geological parameters such as the elastic modulus, hardness, and fracture toughness of deep rock masses is crucial for the evaluation of sweet spots, the design of wellbore trajectories, and the optimization of fracturing schemes. Currently, macro rock mechanics experiments are generally used to measure rock mechanical parameters. However, this method requires relatively large rock samples, which are difficult to obtain. At the same time, the amount of data obtained is relatively limited, and insufficient data support cannot be obtained. In recent years, nano / micro-indentation experiments have been favored in the measurement of rock mechanical parameters. On the one hand, this method has lower requirements for sample specimens and sizes. On the other hand, continuous mechanical data can be obtained by combining well logging data, providing a sufficient data basis for the effective development of reservoirs. However, current research results mainly focus on two aspects: one is the establishment of mechanical parameter models for different lithological rock samples. For example, the patent application document with the title "A method for obtaining macro rock mechanical parameters based on nano-indentation technology" and the publication number CN 119023472 A, which was published on November 6, 2024, records a method for obtaining macro rock mechanical parameters based on nano-indentation technology. This method is based on the macro hardness of rocks obtained by nano-indentation, and successively establishes relationship models for compressive strength, Young's modulus, and Poisson's ratio to obtain key macro rock mechanical parameters. The other is the upgrade of mechanical parameters from the nano scale to the macro scale. For example, the paper "Study on the mechanical properties of laminated shale rocks using nano-indentation experiments" published by Shi Xian et al. in 2019 mainly measured the basic rock mechanical parameters of shale using nano-indentation technology, proposed a three-component rock micro-mechanical model based on mineral strength classification, and used the Mori-Tanaka method to realize the upgrade of mechanical parameters from the nano scale to the centimeter scale.

[0003] Therefore, there is little research on the applicability of nano-indentation experiments under extreme conditions of deep formations and mechanical parameter models at present, and the mechanical parameters of deep reservoirs cannot be effectively obtained, which restricts the effective development of deep reservoir energy. Summary of the Invention

[0004] In order to solve the problem of spatially inverting the key in-situ rock mechanical parameters of a reservoir under extreme conditions, the present invention provides a method and system for spatially inverting and predicting key in-situ rock mechanical parameters of a reservoir.

[0005] To achieve the above object, the present invention provides the following technical solutions: A method for spatially inverting and predicting key in-situ rock mechanical parameters of a reservoir, comprising the following steps: Obtain the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment, and the DP criterion parameters to construct a data set, and construct a macroscopic indentation numerical model according to the data set; Sample the elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample by the Latin hypercube sampling method to form parameter combinations; input the parameter combinations into the macroscopic indentation numerical model for simulation analysis to obtain load-displacement curves under different conditions; extract multiple characteristic values from the load-displacement curves under different conditions as the training set; wherein, the characteristic values characterize the rock mechanics response parameters; Optimize the hyperparameters of the deep neural network based on the Bayesian optimization algorithm, input the training set data into the optimized deep neural network to train the model, and obtain a single-point inversion prediction model; based on the Kriging interpolation algorithm, convert the single-point inversion prediction model into a spatial inversion prediction model; Input the load-displacement data detected in real time of the rock stratum into the spatial inversion prediction model to obtain the spatial distribution data of the rock mechanics parameters in the deep rock stratum.

[0006] Preferably, the sampling of the elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample by the Latin hypercube sampling method is specifically through the following formula: ; Wherein, represents the estimated value; represents the expectation of; represents the number of samplings participating in the computer operation; represents the linear coefficient; represents the standard deviation of; represents the standard deviation of; represents the sampling result value, and the sampling result is any one of the elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample; represents the index subscript.

[0007] Preferably, the optimization of the hyperparameters of the deep neural network based on the Bayesian optimization algorithm is specifically: Based on the Bayesian optimization algorithm, with the mean absolute error as the objective function, optimize the hyperparameters of the deep neural network model, and screen out the optimal deep neural network structure parameter combination.

[0008] Preferably, the obtaining of the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment, and the DP criterion parameters specifically involves performing a mineral component experiment, a scanning electron microscope experiment, and an indentation experiment on the target rock sample to obtain the mineral composition of the target rock sample and the microscopic information of the sample rock after the indentation experiment; at the same time, performing a triaxial stress experiment to obtain the DP criterion parameters.

[0009] Preferably, the multiple characteristic values specifically include: temperature, pressure, maximum applied load, indentation depth corresponding to the maximum applied load, and indentation depth after complete unloading.

[0010] The present invention also provides a system for spatial inversion prediction of key in-situ rock mechanics parameters of a reservoir, specifically including: A data processing module, configured to obtain the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment, and the DP criterion parameters to construct a data set, and construct a macroscopic indentation numerical model according to the data set.

[0011] A parameter calculation module, configured to sample the elastic modulus, Poisson's ratio, compressive strength, temperature, and pressure parameters of the target rock sample by the Latin hypercube sampling method to form a parameter combination; input the parameter combination into the macroscopic indentation numerical model for simulation analysis to obtain load-displacement curves under different conditions; extract multiple characteristic values from the load-displacement curves under different conditions as a training set; wherein, the characteristic values characterize the rock mechanics response parameters.

[0012] A training module, configured to optimize the hyperparameters of the deep neural network based on the Bayesian optimization algorithm, input the training set data into the optimized deep neural network to train the model, and obtain a single-point inversion prediction model; based on the Kriging interpolation algorithm, convert the single-point inversion prediction model into a spatial inversion prediction model.

[0013] A parameter inversion module, configured to input the load-displacement data detected in real time in the rock formation into the spatial inversion prediction model to obtain the spatial distribution data of the rock mechanics parameters in the deep rock formation.

[0014] The present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory, and the processor executes the computer program to implement the steps in the method for spatial inversion prediction of key in-situ rock mechanics parameters of a reservoir.

[0015] The present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is loaded by a processor, it can execute the steps in the method for spatial inversion prediction of key in-situ rock mechanics parameters of a reservoir.

[0016] The reservoir in-situ key rock mechanical parameter spatial inversion prediction method and system provided by the present invention have the following beneficial effects: The present invention establishes a macroscopic indentation numerical model through the experimental parameter information of target rock samples, simulates and analyzes the combined parameters of the target rock samples through the macroscopic indentation numerical model, constructs a training set, and obtains the mechanical parameters of the target rock samples that are more accurate under extreme conditions of deep formations, providing a data basis for subsequent model training. The deep neural network model is trained with the training set data and combined with the Kriging interpolation algorithm to obtain a spatial inversion prediction model. The real-time detection data of deep rock formations are input into the spatial inversion prediction model to obtain the spatial distribution data of rock mechanical parameters in deep rock formations. It can more accurately invert the key rock mechanical parameters in the deep reservoir space, providing a basis for well placement, well trajectory design, and fracturing scheme optimization of deep reservoir energy. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the embodiments of the present invention and their design schemes, the accompanying drawings required for this embodiment will be briefly introduced below. The accompanying drawings in the following description are only partial embodiments of the present invention. For those of ordinary skill in the art, other accompanying drawings can be obtained based on these drawings without creative efforts.

[0018] Figure 1 It is a flowchart of a method for spatial inversion prediction of key rock mechanical parameters in-situ in a reservoir according to an embodiment of the present invention.

[0019] Figure 2 It is a schematic diagram of the DP criterion and the principle of rock damage discrimination in an embodiment of the present invention, where Figure 2 (a) is the p-t relationship curve, Figure 2 (b) is the D-P yield surface, Figure 2 and (c) is the stress-strain relationship curve during the rock damage process.

[0020] Figure 3 It is the nano-indentation numerical model in an embodiment of the present invention.

[0021] Figure 4 It is the macroscopic indentation numerical model in an embodiment of the present invention.

[0022] Figure 5 It is a comparison chart of the predicted elastic modulus and the true elastic modulus in an embodiment of the present invention.

[0023] Figure 6 It is the spatial distribution characteristic of the elastic modulus in an embodiment of the present invention.

[0024] Figure 7 It is a flowchart of a method for spatial inversion prediction of key rock mechanical parameters in-situ in a reservoir according to the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0025] To enable those skilled in the art to better understand the technical solution of the present invention and be able to implement it, the present invention will be described in detail below in conjunction with the accompanying drawings and specific embodiments. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and cannot be used to limit the protection scope of the present invention.

[0026] Embodiment In Area A, the general burial depth of the favorable producing formation is greater than 3500m. The high-temperature and high-pressure geological environment has a huge impact on the mechanical properties of reservoir rocks. It is crucial to clarify the distribution characteristics of mechanical parameters in the in-situ environment of the reservoir for the efficient development of the reservoir.

[0027] The present invention provides a method for spatially inverting and predicting key rock mechanical parameters in-situ in a reservoir. Taking Area A as an example, as Figure 1 shown, it specifically includes the following steps: Step 1: Select typical deep wells in Area A, continuously collect cuttings samples of the target layer in the horizontal section at intervals of 10m (at least 3 samples are collected at each sampling point), prepare them into square specimens of 15mm×15mm×8mm, then polish the test surface with sandpapers of different mesh numbers in turn to reduce the surface roughness, and finally use argon ion polishing to further reduce the surface roughness to complete the specimen preparation work. According to the actual geothermal conditions of the target layer, set high-temperature conditions of 100 - 200°C, carry out nanoindentation experiments under high temperature, obtain load-displacement relationship curves at different temperatures and depths, and establish an indentation experiment database.

[0028] Step 2: For the samples in Step 1, carry out mineral composition experiments and scanning electron microscope experiments, obtain the mineral composition of the target rock samples and conduct microtopography analysis of the cross-section after the indentation experiment, and carry out triaxial stress experiments to obtain the key parameters required by the DP criterion. Based on the microscopic information of rock minerals and the relevant parameters of the DP criterion, use Abaqus finite element numerical simulation software to establish a nanoindentation numerical model, as Figure 3 shown. The nanoindentation numerical model includes two parts: the indenter and the specimen. Among them, referring to the indoor experimental process, during the nanoindentation test, the specimen is fixed on the stage and the indenter only moves in the loading direction. Therefore, a positive displacement constraint in the Y direction is applied at the bottom of the numerical model, and an axisymmetric constraint is applied on the axis of symmetry. At the same time, the horizontal movement and rotation of the indenter in the model are constrained so that the indenter can only move in the vertical direction. The numerical specimen considers the basic material parameters and thermodynamic parameters of the rock, and the specific parameters are shown in Table 1.

[0029] Table 1 Basic material parameters of the nanoindentation numerical model

[0030] The nanoindentation numerical model uses the Drucker-Prager criterion (DP criterion for short) to characterize the elastic-plastic deformation process of rocks during the indentation process. The DP criterion is an optimization and generalization of the Mohr-Coulomb criterion and the Mises criterion. The DP criterion and the rock damage discrimination principle are as follows: Figure 2 As shown, Figure 2 The PT curve in (a) refers to the yield surface of the linear DP model on the meridian plane, where is the friction angle, The cohesion of the material. Figure 2 (b) is the DP yield surface, which is the typical yield / flow surface of the DP linear model in the deviatoric plane, where K is the ratio of the triaxial tensile yield stress to the triaxial compressive yield stress. Figure 2 (c) is the stress-strain relationship curve during the rock damage process, where σ 0 is the elastic limit stress of rock; D is the damage degree of rock; y0 is the yield stress when the rock begins to be damaged, at which time the damage factor D=0; when the plastic strain increases to ε f When the damage factor D=1, the rock fails completely and falls off. Its expression is: ; ; ; ; ; ; in, is the internal friction angle of geotechnical materials; is the cohesion of geotechnical materials; and Both are cohesion with rock and soil materials and internal friction angle The relevant constants; is the first invariant of the stress tensor; is the second invariant of the stress deviator tensor; , , are the first principal stress, the second principal stress and the third principal stress respectively; D is the damage factor; , are Young’s modulus of the material before and after, respectively; For strain; is the equivalent plastic strain when the material begins to be damaged; for stress; is the stress after material damage.

[0031] Compare the load-displacement curves at different temperatures and different depths obtained from the nanoindentation numerical model with those from nanoindentation experiments to verify the accuracy of the nanoindentation numerical model.

[0032] Step 3: Since in actual engineering, the scale is larger, generally at the centimeter scale. Therefore, based on the verified accurate nanoindentation numerical model parameters, further establish a macroscopic indentation numerical model, as Figure 4 shown. Combine the Latin hypercube sampling method to sample the reservoir temperature and pressure conditions and basic rock mechanics parameters that affect the experimental results. The reservoir temperature and pressure conditions and basic rock mechanics parameters include elastic modulus, Poisson's ratio, compressive strength, temperature, and pressure. According to the actual formation conditions in Area A, the temperature setting range is 100 - 200 °C, the confining pressure setting range is 10 - 80 MPa, the elastic modulus setting range is 15 - 60 GPa, the Poisson's ratio setting range is 0.15 - 0.3, and the compressive strength setting range is 100 - 500 MPa. Randomly select 200 groups of parameters to form a parameter combination database as the conditions for carrying out numerical experiments.

[0033] Latin hypercube sampling method: ; where: represents the estimated value; represents the expectation of; represents the number of samples participating in the computer operation; represents the linear coefficient; represents the standard deviation of; represents the standard deviation of; represents the sampling result value, and the sampling result is any one of the elastic modulus, Poisson's ratio, compressive strength, temperature, and pressure parameters of the target rock sample; represents the index subscript.

[0034] Apply these combined parameters to the macroscopic indentation numerical model for simulation analysis to obtain the load-displacement curves under different conditions.

[0035] Step 4: Extract three parameters, namely the maximum applied load, the corresponding indentation depth during loading, and the indentation depth after complete unloading, from the inverted elastic modulus target parameter according to the load-displacement curves under different conditions. Set these three parameters together with the two parameters of temperature and pressure as characteristic values, set the elastic modulus data as the target value, and store the data in the database as the inversion training data set for the trained neural network.

[0036] Step 5: Based on the Bayesian optimization algorithm, with the mean absolute error as the objective function, optimize the hyperparameter problem of the deep neural network method, screen out the optimal combination of deep neural network parameters (as shown in Table 2), and build a deep neural network model accordingly. Input the training set data into the optimized deep neural network to train the model and obtain a single-point inversion prediction model. Taking the elastic modulus output by the model as a reference, analyze the accuracy of the model training results. After 1000 training times, the final result is obtained, and the MAE value of the model is 32.143. The comparison result between the predicted elastic modulus and the true elastic modulus is as Figure 5 shown, indicating that the model predicts the elastic modulus with a very small loss, and thus a single-point inversion prediction model is obtained.

[0037] The calculation method of the inversion error is as follows: ; where MAE is the mean absolute error, and are the target value and the predicted value respectively, is the number of data.

[0038] Table 2 Optimal combination parameter table of the deep neural network

[0039] Based on the Kriging interpolation algorithm, convert the single-point inversion prediction model into a spatial inversion prediction model. Specifically: obtain the elastic model parameters output by the single-point inversion prediction model; select a variogram model for fitting; construct and solve the Kriging equations according to the fitting results to obtain the weight coefficients of each known point; use the weight coefficients for spatial prediction, generate a spatial grid, and obtain the distribution law and variation characteristics of the elastic modulus in space, so as to realize the conversion of the single-point inversion prediction model into a spatial inversion prediction model.

[0040] Step 6: The on-site drilling technology can realize the downhole placement of the indentation detection device while drilling, and collect the load-displacement data of the wellbore rock detection in real time. Through the in-situ point method loading test, obtain the downhole rock mechanics parameters in real time according to the spatial inversion prediction model, realize the real-time establishment and update of the spatial distribution data of the deep rock mass real-time underground rock mechanics parameters, and a spatial rock mechanics parameter field can be established, and the result is as Figure 6 shown.

[0041] The present invention also provides a system for the spatial inversion prediction method of the key in-situ rock mechanics parameters of a reservoir, specifically including: A data processing module, used to obtain the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment, and the DP criterion parameters to construct a data set, and construct a macroscopic indentation numerical model according to the data set.

[0042] A parameter calculation module is used to sample the elastic modulus, Poisson's ratio, compressive strength, temperature, and pressure parameters of a target rock sample through the Latin hypercube sampling method to form parameter combinations; input the parameter combinations into a macroscopic indentation numerical model for simulation analysis to obtain load-displacement curves under different conditions; extract multiple characteristic values from the load-displacement curves under different conditions as a training set; where the characteristic values characterize the rock mechanics response parameters.

[0043] A training module is used to optimize the hyperparameters of a deep neural network based on the Bayesian optimization algorithm, input the training set data into the optimized deep neural network to train the model, and obtain a single-point inversion prediction model; based on the Kriging interpolation algorithm, convert the single-point inversion prediction model into a spatial inversion prediction model.

[0044] A parameter inversion module is used to input the load-displacement data detected in real time in the rock formation into the spatial inversion prediction model to obtain the spatial distribution data of the rock mechanics parameters in the deep rock formation.

[0045] Each module in the above system for spatially inverting and predicting key rock mechanics parameters in-situ in a reservoir can be implemented in whole or in part by software, hardware, and their combination. The above modules can be embedded in the processor of a computer device in hardware form or be independent of it, or can be stored in the memory of the computer device in software form to facilitate the processor to call and execute the operations corresponding to the above modules.

[0046] The present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory. The processor executes the computer program to implement the steps in an embodiment of a method for spatially inverting and predicting key rock mechanics parameters in-situ in a reservoir. The specific implementation method can refer to the method embodiment and will not be elaborated here.

[0047] Furthermore, the present invention also provides a non-transitory computer-readable storage medium containing instructions, and a computer program is stored on the storage medium. For example, a memory containing instructions, and the above instructions can be executed by the processor of a computer device to complete the above method. For example, the non-transitory computer-readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, and an optical data storage device, etc. When the computer program is executed by the processor, it can implement the steps in an embodiment of a method for spatially inverting and predicting key rock mechanics parameters in-situ in a reservoir. The specific implementation method can refer to the method embodiment and will not be elaborated here.

[0048] Those skilled in the art should understand that the embodiments of the present invention may provide a method, a system or a computer program product. Therefore, the present invention may take the form of an all-hardware embodiment, an all-software embodiment or an embodiment combining software and hardware aspects. Moreover, the present invention 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.) that contain computer-usable program code.

[0049] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems) and computer program products according to embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of flows and / or blocks in the flowchart and / or block diagram can also be implemented. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for realizing the functions specified in Figure 1 one or more flows and / or blocks Figure 1 one or more blocks.

[0050] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including an instruction device, and the instruction device realizes the functions specified in Figure 1 one or more flows and / or blocks Figure 1 one or more blocks.

[0051] These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for realizing the functions specified in Figure 1 one or more flows and / or blocks Figure 1 one or more blocks.

[0052] It should be pointed out that the specific implementation methods described above can enable those skilled in the art to understand the invention more comprehensively, but do not limit the invention in any way. Therefore, although the invention has been described in detail in this specification and embodiments, those skilled in the art should understand that the invention can still be modified or replaced by equivalents; and all technical solutions and improvements that do not deviate from the spirit and scope of the invention are included in the protection scope of the patent for the invention. Any figure mark in the claims should not be regarded as limiting the claims involved. Any simple change or equivalent replacement of the technical solution that can be obviously obtained by any technician familiar with the field within the technical scope disclosed in the present invention belongs to the protection scope of the present invention.

Claims

1. A spatial inversion prediction method for in-situ key rock mechanical parameters of a reservoir, characterized in that: The following steps are involved: Acquire the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment, and the DP criterion parameters to construct a data set, and construct a macroscopic indentation numerical model based on the data set; The elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample are sampled by the Latin hypercube sampling method to form a parameter combination; the parameter combination is input into the macro indentation numerical model for simulation analysis to obtain the load-displacement curve under different conditions; multiple eigenvalues ​​are extracted from the load-displacement curve under the different conditions as a training set; wherein the eigenvalues ​​represent the rock mechanical response parameters; The hyperparameters of the deep neural network are optimized based on the Bayesian optimization algorithm, and the training set data is input into the optimized deep neural network to train the model to obtain a single-point inversion prediction model; based on the Kriging interpolation algorithm, the single-point inversion prediction model is converted into a spatial inversion prediction model; The load-displacement data detected in real time in the rock formation is input into the spatial inversion prediction model to obtain the spatial distribution data of rock mechanical parameters in the deep rock formation.

2. A spatial inversion prediction method for in-situ key rock mechanical parameters of a reservoir according to claim 1, characterized in that: The elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample are sampled by the Latin hypercube sampling method, specifically by the following formula: ； in, Indicates an estimated value; express expectations; Indicates the number of samples participating in the computer operation; represents the linear coefficient; express The standard deviation of express The standard deviation of Indicates the sampling result value, where the sampling result is any one of the elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample; Indicates an index subscript.

3. The spatial inversion prediction method for in-situ key rock mechanical parameters of a reservoir according to claim 1, characterized in that: The hyperparameters of the deep neural network are optimized based on the Bayesian optimization algorithm, specifically: Based on the Bayesian optimization algorithm, the mean absolute error is used as the objective function to optimize the hyperparameters of the deep neural network model and screen out the optimal combination of deep neural network structure parameters.

4. The spatial inversion prediction method for in-situ key rock mechanical parameters of a reservoir according to claim 1, characterized in that: The method of obtaining the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment and the DP criterion parameters specifically comprises performing a mineral composition experiment, a scanning electron microscope experiment and an indentation experiment on the target rock sample to obtain the mineral composition of the target rock sample and the microscopic information of the sample rock after the indentation experiment; and performing a triaxial stress experiment at the same time to obtain the DP criterion parameters.

5. The spatial inversion prediction method for in-situ key rock mechanical parameters of a reservoir according to claim 1, characterized in that: The multiple characteristic values ​​specifically include: temperature, pressure, maximum applied load, the indentation depth corresponding to the maximum applied load, and the indentation depth after complete unloading.

6. A spatial inversion prediction system for in-situ key rock mechanical parameters of reservoirs, characterized in that: include: A data processing module is used to obtain the mineral composition of the target rock sample, the microscopic information of the sample rock after the indentation experiment, and the DP criterion parameters to construct a data set, and to construct a macroscopic indentation numerical model based on the data set; The parameter calculation module is used to sample the elastic modulus, Poisson's ratio, compressive strength, temperature and pressure parameters of the target rock sample by Latin hypercube sampling method to form a parameter combination; the parameter combination is input into the macro indentation numerical model for simulation analysis to obtain the load-displacement curve under different conditions; a plurality of characteristic values ​​are extracted from the load-displacement curve under the different conditions as a training set; wherein the characteristic value represents the rock mechanical response parameter; A training module is used to optimize the hyperparameters of the deep neural network based on the Bayesian optimization algorithm, input the training set data into the optimized deep neural network to train the model, and obtain a single-point inversion prediction model; based on the Kriging interpolation algorithm, the single-point inversion prediction model is converted into a spatial inversion prediction model; The parameter inversion module is used to input the load-displacement data detected in real time in the rock formation into the spatial inversion prediction model to obtain the spatial distribution data of rock mechanics parameters in the deep rock formation.

7. A computer device comprising a memory, a processor and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is loaded into a processor, it can execute the steps of the method according to any one of claims 1 to 5.

Images