Method and device for determining main control factors of mechanical properties of heterogeneous dolomite

Abstract

The application discloses a method and device for determining main control factors of mechanical properties of non-homogeneous dolomite, relates to the technical field of rock mechanics, and aims to improve the accuracy and efficiency of determining main control factors of mechanical properties of dolomite. The main technical scheme of the application is as follows: the discrete element method is used to obtain the particle contact force distribution according to the contact information between particles in the microstructure model of dolomite; the numerical method is used to simulate the flow of pore water in the pore network structure of the microstructure model of dolomite, and the pore pressure distribution is obtained; the thermodynamic equation is used to determine the temperature change effect according to temperature data; the multi-physical field coupling model is constructed according to the particle contact force distribution, the pore pressure distribution and the temperature change effect; if the calculation complexity of the multi-physical field coupling model exceeds a preset threshold, the multi-physical field coupling model is simplified to obtain a simplified model; and the multi-physical field coupling model and the simplified model are used to determine the main control factors of the mechanical properties of dolomite under different environmental conditions.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of rock mechanics, and in particular to a method and device for determining main control factors of mechanical properties of heterogeneous dolomite. BACKGROUND

[0002] Dolomite is a kind of heterogeneous sedimentary rock formed by different sizes of chert particles through different ways of nucleation or under the action of geological movement. Its heterogeneity and anisotropy make its mechanical behavior in engineering particularly complex. These characteristics directly affect the stability of the operation, especially in the application of dolomite formation drilling engineering. Since the mechanical properties of heterogeneous dolomite are affected by many factors such as inter-particle contact force, pore water pressure and temperature change, determining these main control factors is crucial to prevent potential risks in engineering.

[0003] At present, in the method for determining the main control factors of the mechanical properties of heterogeneous dolomite, the average processing method is used to process the contact force distribution between particles according to the macroscopic model. In addition, the pore water pressure is usually estimated by the fluid mechanics equation, and the temperature change is processed by the coupling model of thermodynamics and mechanics.

[0004] However, this method has several major problems: first, in the problem of inter-particle contact force, the macroscopic model may cause deviation in predicting the actual mechanical behavior because the average processing ignores the real complexity of the inter-particle contact force; second, it may be difficult to depict the complexity of the pore structure and the flow path of the fluid using the fluid mechanics equation, which makes the calculation of the pore water pressure not accurate enough; third, this method processes the particle contact force, pore water pressure and temperature change separately, which to some extent increases the processing flow of the scheme and reduces the processing efficiency; not to mention that for the temperature change, a multi-physical field coupling model of thermodynamics and mechanics is used, which often requires a large amount of computing resources, further affecting the processing efficiency and accuracy. In summary, the existing method has obvious deficiencies in accuracy and efficiency. In order to more accurately simulate the mechanical properties of dolomite and improve the processing efficiency, a new method is needed to overcome the above challenges. SUMMARY

[0005] In view of the above problems, the present application provides a method and device for determining main control factors of the mechanical properties of heterogeneous dolomite, the main purpose of which is to improve the accuracy and efficiency of determining the main control factors of the mechanical properties of heterogeneous dolomite.

[0006] To solve the above technical problems, the present application proposes the following solutions:

[0007] In a first aspect, the present application provides a method for determining main control factors of the mechanical properties of heterogeneous dolomite, the method comprising:

[0008] a microstructure model of the heterogeneous dolomite is reconstructed according to microstructure data of the heterogeneous dolomite sample by using a three-dimensional scanning technology;

[0009] a contact force distribution of the particles is obtained by calculating the contact force between the particles according to contact information between the particles in the microstructure model by using a discrete element method;

[0010] a pore pressure distribution is obtained by simulating the flow of pore water in the pore network structure of the microstructure model by using a numerical method;

[0011] a temperature variation effect is obtained according to the temperature variation quantity by using a thermodynamic equation in combination with collected temperature data;

[0012] a multi-physics field coupling model is constructed according to the contact force distribution of the particles, the pore pressure distribution and the temperature variation effect, so as to simulate the mechanical response of the heterogeneous dolomite;

[0013] if the calculation complexity of the multi-physics field coupling model exceeds a preset threshold, the multi-physics field coupling model is simplified by using a dimension reduction technology to obtain a simplified model;

[0014] stress distribution and strain distribution data of the heterogeneous dolomite under different environmental conditions are obtained by using the multi-physics field coupling model;

[0015] mechanical characteristic parameters of the heterogeneous dolomite are extracted from the stress distribution and strain distribution data based on the simplified model, and a main control factor of the mechanical characteristics of the heterogeneous dolomite under different environmental conditions is determined according to the mechanical characteristic parameters.

[0016] In a second aspect, the present application provides a device for determining a main control factor of mechanical characteristics of a heterogeneous dolomite, the device comprising:

[0017] a first model construction unit configured to reconstruct a microstructure model of the heterogeneous dolomite according to microstructure data of a heterogeneous dolomite sample by using a three-dimensional scanning technology;

[0018] a contact force calculation unit configured to calculate the contact force between the particles according to contact information between the particles in the microstructure model constructed by the first model construction unit by using a discrete element method, so as to obtain a contact force distribution of the particles;

[0019] a pore pressure calculation unit configured to simulate the flow of pore water in the pore network structure of the microstructure model constructed by the first model construction unit by using a numerical method, so as to obtain a pore pressure distribution;

[0020] a temperature effect calculation unit configured to determine a temperature change amount by using a thermodynamic equation in combination with the collected temperature data, and obtain a temperature change effect according to the temperature change amount;

[0021] a second model construction unit configured to construct a multi-physics field coupling model according to the particle contact force distribution calculated by the contact force calculation unit, the pore pressure distribution calculated by the pore pressure calculation unit, and the temperature change effect calculated by the temperature effect calculation unit, so as to simulate the mechanical response of the heterogeneous dolomite;

[0022] a third model construction unit configured to simplify the multi-physics field coupling model by using a dimension reduction technique if the calculation complexity of the multi-physics field coupling model constructed by the second model construction unit exceeds a preset threshold, so as to obtain a simplified model;

[0023] a coupling model application unit configured to obtain stress distribution and strain distribution data of the heterogeneous dolomite under different environmental conditions by using the multi-physics field coupling model obtained by the second model construction unit;

[0024] a main control factor determination unit configured to extract mechanical characteristic parameters of the heterogeneous dolomite from the stress distribution and strain distribution data obtained by the coupling model application unit based on the simplified model obtained by the third model construction unit, and determine the main control factor of the mechanical characteristics of the heterogeneous dolomite under different environmental conditions according to the mechanical characteristic parameters.

[0025] In order to achieve the above-mentioned purpose, according to a third aspect of the present application, a computer device is provided, comprising a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the non-heterogeneous dolomite mechanical characteristic main control factor determination method of the first aspect.

[0026] In order to achieve the above-mentioned purpose, according to a fourth aspect of the present application, a computer readable storage medium is provided, which stores a computer program, wherein the computer program is executed by a processor to implement the non-heterogeneous dolomite mechanical characteristic main control factor determination method of the first aspect.

[0027] In order to achieve the above-mentioned purpose, according to a fifth aspect of the present application, a computer program product is provided, comprising a computer program, wherein the computer program is executed by a processor to implement the non-heterogeneous dolomite mechanical characteristic main control factor determination method of the first aspect.

[0028] By the technical scheme, the method and device for determining the main control factor of the mechanical properties of the heterogeneous dolomite are provided, first, the microstructure data of the heterogeneous dolomite sample is obtained based on the three-dimensional scanning technology, and the microstructure model is reconstructed according to the microstructure data. Then, the discrete element method is used to calculate the contact force distribution between the particles in the model. Compared with the traditional macroscopic model, this method can more truly reflect the complexity of the contact force between the particles. Next, the numerical method is used to simulate the pore water flow in the pore network of the microstructure model. Compared with the traditional fluid mechanics equation, this method can more accurately capture the complexity of the pore structure and the flow path of the fluid, thereby obtaining a more accurate pore pressure distribution. In addition, combined with the collected temperature data, the amount of temperature change is calculated by the thermodynamic equation, and the effect is analyzed. Based on the particle contact force distribution, the pore pressure distribution and the temperature change effect, a multi-physical field coupling model is constructed to simulate the mechanical response of the heterogeneous dolomite. If the calculation complexity of the multi-physical field coupling model exceeds the preset threshold, the dimension reduction technique can be used to simplify the model to obtain a simplified model. The stress distribution and strain distribution data of the heterogeneous dolomite under different environmental conditions are obtained by the multi-physical field coupling model, and the mechanical characteristic parameters are extracted from the data by the simplified model, and then the main control factor of the mechanical characteristics of the heterogeneous dolomite under different environmental conditions is determined. The scheme adopts the multi-physical field coupling model in the key part to ensure the accuracy, and adopts the simplified model in other parts to improve the efficiency, so as to ensure the accuracy of the results while significantly improving the work efficiency. Through the improvement of the accuracy and the efficiency, it is also helpful to better predict and prevent potential risks in engineering, thereby improving the stability of the operation.

[0029] The above description is only a summary of the technical scheme of the present application. In order to more clearly understand the technical means of the present application, the content of the specification can be implemented, and in order to make the above and other purposes, characteristics and advantages of the present application more obvious and easy to understand, the specific embodiments of the present application are described below. BRIEF DESCRIPTION OF DRAWINGS

[0030] By reading the detailed description of the preferred embodiments below, various other advantages and benefits will become apparent to those of ordinary skill in the art. The drawings are only for the purpose of illustrating the preferred embodiments and are not considered to be limiting on the present application. Moreover, the same reference symbols are used throughout the drawings to represent the same parts. In the drawings:

[0031] Figure 1 A flow chart of a method for determining the main control factor of the mechanical properties of the heterogeneous dolomite is shown;

[0032] Figure 2 A flow chart of another method for determining the main control factor of the mechanical properties of the heterogeneous dolomite is shown;

[0033] Figure 3 A composition block diagram of a device for determining a main control factor of mechanical properties of heterogeneous dolomite according to an embodiment of the present application is shown;

[0034] Figure 4 A composition block diagram of another device for determining a main control factor of mechanical properties of heterogeneous dolomite according to an embodiment of the present application is shown. DETAILED DESCRIPTION

[0035] Exemplary embodiments of the present disclosure will be described more fully hereinafter with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it is understood that the present disclosure can be embodied in various forms and should not be limited by the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the present disclosure to those skilled in the art.

[0036] In view of the limitations in the current method for determining the main control factor of the mechanical properties of heterogeneous dolomite, the inventors have realized that by using the microstructure data of the heterogeneous dolomite sample to reconstruct its microstructure model and obtaining the particle contact force distribution and pore pressure distribution based on this fine model, the internal state of the heterogeneous dolomite can be more accurately reflected compared to traditional macroscopic models and fluid mechanics equations. Further, the inventors have proposed a multi-physical field coupling model that combines the contact force distribution between particles, pore pressure distribution, and temperature variation effect to comprehensively capture the complex behavior of heterogeneous dolomite under various physical fields.

[0037] However, multi-physical field coupling models often have high requirements for computing resources. In order to effectively address this challenge while maintaining accuracy and computational precision, the inventors have introduced the concept of a simplified model, which is a reasonable simplification of the multi-physical field coupling model that retains key information while reducing computational complexity. This strategy allows the use of the multi-physical field coupling model when high-precision analysis is required, while using the simplified model for rapid calculation in other scenarios, thereby achieving both accuracy and efficiency.

[0038] In summary, the method for determining the main control factor of the mechanical properties of heterogeneous dolomite provided by the present application combines microstructure modeling, multi-physical field coupling analysis, and model simplification techniques, which not only deepens the understanding of the main control factor of the mechanical properties of heterogeneous dolomite, but also optimizes the calculation process, achieving both precision and efficiency.

[0039] Next, a method for determining the main control factor of the mechanical properties of heterogeneous dolomite according to an embodiment of the present application will be described in conjunction with Figure 1 A method for determining the main control factor of the mechanical properties of heterogeneous dolomite according to an embodiment of the present application will be described, and the specific execution steps are as shown in Figure 1 ​

[0040] 101. Reconstruct the microstructure model of the heterogeneous dolomite using three-dimensional scanning technology based on the microstructure data of the heterogeneous dolomite sample.

[0041] In this process, first, a representative heterogeneous dolomite sample is selected to ensure that it matches the scanning range of the three-dimensional scanning device. After selecting the sample, it is cleaned to remove dust and other impurities to ensure the accuracy of the scanning data.

[0042] After cleaning, the parameters of the high-precision three-dimensional scanning device are configured, the scanning resolution is set to the micron level, and the scanning path is planned according to the specific shape of the heterogeneous dolomite sample to ensure full coverage of the sample.

[0043] Next, an advanced laser scanner or structured light scanner is used to follow the pre-set scanning resolution and path to perform fine scanning of the heterogeneous dolomite sample, aiming to capture the details of its microstructure in detail.

[0044] After scanning, the three-dimensional point cloud data of the sample surface is obtained and filtered to effectively remove noise points and abnormal points to improve data quality.

[0045] Subsequently, three-dimensional reconstruction algorithms such as Poisson reconstruction algorithm or triangular meshing algorithm are used to convert the processed point cloud data into a three-dimensional surface model, providing a basis for extracting the microstructure features of the heterogeneous dolomite sample.

[0046] Based on the constructed three-dimensional surface model, the microstructure features of the heterogeneous dolomite sample are further extracted, including particle shape, pore distribution, etc., and these features are quantitatively analyzed to generate microstructure parameters such as porosity and particle size distribution, thereby constructing a complete microstructure model.

[0047] Finally, the microstructure model is exported in a common file format (such as STL, OBJ, etc.) to facilitate subsequent analysis and application. This process not only ensures the accuracy of the operation, but also significantly improves work efficiency and data utility.

[0048] 102. Calculate the contact force between particles based on the contact information between particles in the microstructure model using the discrete element method to obtain the distribution of particle contact forces.

[0049] In this step, the heterogeneous dolomite microstructure model reconstructed by three-dimensional scanning technology can be imported into the discrete element analysis software. This model contains information about the shape, size, and position of the particles, ensuring that it accurately reflects the microstructure of the heterogeneous dolomite.

[0050] Next, physical properties such as elastic modulus, Poisson's ratio, and density can be set for each particle in the model. These properties need to be set according to the actual mineral composition of the heterogeneous dolomite to ensure the accuracy of the calculation. Furthermore, the calculation of contact forces between particles depends on the choice of contact model. For complex multi-mineral materials like heterogeneous dolomite, the Hertzian contact model can be used because it better describes the elastic deformation behavior between particles.

[0051] Based on experimental data or literature references, reasonable parameters can be set for the contact model. For example, the normal stiffness can be set to 1e8 N / m, the tangential stiffness to 5e7 N / m, and the friction coefficient to 0.5. These parameters directly affect the interaction between particles and therefore need to be adjusted according to the actual situation.

[0052] To simulate actual stress conditions, the bottom boundary of the model can be set as a fixed boundary, while the top boundary can be set as a movable boundary. This setup can simulate the mechanical behavior of heterogeneous dolomite under compressive loads. Subsequently, a uniform compressive load of 10 MPa is applied to the top of the model. This load can be adjusted according to actual research needs, such as simulating formation pressure at different depths.

[0053] After completing the above steps, the calculation can be initialized in the discrete element method (DEM) software to ensure that the contact relationships between particles are correctly established. This step is crucial because the initial conditions directly affect the accuracy of the calculation results. Once the calculation starts, the software iteratively solves the particle motion equations and contact force equations. During the calculation, the particle displacement and energy changes can be monitored in real time to ensure the convergence of the calculation results.

[0054] After calculation, the contact force data between each particle can be extracted, including normal and tangential forces. This data can be exported as a CSV file for subsequent analysis and visualization. Then, post-processing software (such as ParaView or MATLAB) can be used to visualize the contact force data. By plotting the distribution of contact forces between particles, the force transmission path and concentration areas can be visually observed, thus obtaining the particle contact force distribution.

[0055] The distribution characteristics of contact forces can also be statistically analyzed, such as the maximum contact force, average contact force, and the range of contact force distribution. This data can help understand the mechanical behavior of the microstructure of heterogeneous dolomite. For example, it may be found that the contact force in some areas is significantly higher than in others, which could be related to the shape and arrangement of the grains.

[0056] To ensure the accuracy of the calculation results, the contact force distribution calculated by the discrete element method can be compared with experimental results. For example, the rationality of the contact force distribution can be verified through indentation experiments or acoustic emission experiments. If a large error is found, the contact model parameters or particle properties can be adjusted, and the calculation can be repeated.

[0057] Through the above-described embodiments, this invention can efficiently and accurately calculate the distribution of contact forces between particles in the microstructure of heterogeneous dolomite. This method not only reveals the laws governing the micromechanical behavior of heterogeneous dolomite but also provides important evidence for the study of its macroscopic mechanical properties. Compared with traditional macroscopic models, this invention has higher accuracy and broader application prospects.

[0058] 103. Using numerical methods, the flow of pore water in the pore network structure of the microstructure model is simulated to obtain the pore pressure distribution.

[0059] In this step, the microstructure model of the heterogeneous dolomite reconstructed using 3D scanning technology can be imported into numerical analysis software. This model contains information on the shape, size, and connectivity of the pore network, ensuring that it accurately reflects the pore structure of the heterogeneous dolomite.

[0060] Next, fluid properties, such as water density, viscosity, and compressibility, can be set for the pore network in the model. These properties need to be set according to the physical characteristics of the actual fluid to ensure the accuracy of the calculation. Furthermore, to simulate the actual flow behavior of water, the permeability and porosity parameters of the pore network can be set. These parameters directly affect the water flow velocity and pressure distribution, and therefore need to be adjusted according to the actual situation.

[0061] To simulate real-world flow conditions, the model's inlet boundary can be set as a pressure inlet, and the outlet boundary as a pressure outlet. This setup simulates the flow behavior of pore water under pressure. Subsequently, a pressure of 1 MPa is applied at the model inlet, and the outlet is set to atmospheric pressure. This pressure can be adjusted according to actual research needs, such as simulating groundwater flow at different depths.

[0062] After completing the above steps, you can initialize the calculation in the numerical software to ensure the connectivity of the porous network and the correct setting of boundary conditions. This step is crucial because the initial conditions directly affect the accuracy of the calculation results. Once the calculation starts, the software will simulate the flow of water through the porous network. During the calculation, you can monitor the changes in water velocity and pressure in real time to ensure the convergence of the calculation results.

[0063] After calculation, pressure distribution data in the pore network can be extracted. This data can be exported as a CSV file for subsequent analysis and visualization. Then, post-processing software (such as ParaView or MATLAB) can be used to visualize the pressure distribution data. By plotting the pore pressure distribution map, the pressure transmission path and concentration areas can be visually observed, thus obtaining the pore pressure distribution.

[0064] We can also statistically analyze the characteristics of pressure distribution, such as maximum pressure, average pressure, and pressure gradient. This data can help us understand the patterns of pore water flow. For example, we might find that the pressure in some areas is significantly higher than in others, which could be related to the connectivity and shape of the pores.

[0065] To ensure the accuracy of the calculation results, the pressure distribution obtained from the numerical calculation can be compared with experimental results. For example, the rationality of the pressure distribution can be verified through permeation experiments or pressure sensor measurements. If a large error is found, the fluid property parameters or boundary conditions can be adjusted, and the calculation can be repeated.

[0066] Through the above-described embodiments, this invention can efficiently and accurately simulate the actual flow of pore water in the microstructure of heterogeneous dolomite and obtain the pore pressure distribution. This method not only reveals the laws governing pore water flow but also provides important evidence for studying the seepage characteristics of heterogeneous dolomite. Compared with traditional macroscopic models, this invention has higher accuracy and broader application prospects.

[0067] 104. Use thermodynamic equations in conjunction with collected temperature data to determine the amount of temperature change, and obtain the temperature change effect based on the amount of temperature change.

[0068] In this step, multiple temperature data points can be acquired using a pre-set temperature acquisition device. These temperature data points can come from experimental measurements or on-site monitoring to ensure that they accurately reflect the temperature changes of heterogeneous dolomite under different environments.

[0069] Next, we can use thermodynamic equations to calculate the temperature change based on multiple temperature data, and then obtain the temperature change effect based on the temperature change.

[0070] Among them, the temperature change effect can refer to the changes in the distribution of interparticle contact force and pore pressure caused by temperature changes.

[0071] 105. A multi-physics coupling model was constructed based on the particle contact force distribution, pore pressure distribution, and temperature change effect to simulate the mechanical response of heterogeneous dolomite.

[0072] 106. If the computational complexity of the multiphysics coupling model exceeds a preset threshold, the multiphysics coupling model is simplified using dimensionality reduction techniques to obtain a simplified model.

[0073] 107. Using a multiphysics coupling model, obtain stress and strain distribution data of heterogeneous dolomite under different environmental conditions.

[0074] 108. Based on a simplified model, mechanical characteristic parameters of heterogeneous dolomite are extracted from stress distribution and strain distribution data, and the main controlling factors of mechanical characteristics of heterogeneous dolomite under different environmental conditions are determined based on the mechanical characteristic parameters.

[0075] In steps 105 to 108, a multiphysics coupling model was constructed based on detailed data such as particle contact force distribution, pore pressure distribution, and temperature change effects. However, when the computational complexity of this model exceeds a preset threshold (e.g., excessive computation time or high resource consumption), dimensionality reduction techniques are needed to simplify the model. Specifically, principal component analysis (PCA) can be used to extract the core features from the multiphysics coupling model, thereby effectively reducing the data dimensionality and constructing a simplified model. Another approach is to decompose the data of the multiphysics coupling model, retaining only the main modes, thus reducing the computational burden and obtaining a simplified model as well. These simplified models, while ensuring computational accuracy, can significantly improve computational efficiency, making them suitable for large-scale simulation and real-time analysis scenarios.

[0076] After obtaining the simplified model, a multiphysics coupling model can be used to simulate the mechanical response of heterogeneous dolomite under various environmental conditions (such as different temperatures, pressures, or fluid flow conditions). Then, using the simplified model, mechanical characteristic parameters of the heterogeneous dolomite are extracted from stress and strain distribution data. Based on these mechanical characteristic parameters, the main controlling factors of the mechanical characteristics of heterogeneous dolomite under different environmental conditions can be identified.

[0077] Through the above process, this invention can efficiently and accurately simulate the mechanical response of heterogeneous dolomite under multiphysics coupling environments, and by utilizing simplified models and in-depth data analysis, reveals the main controlling factors of its mechanical characteristics. This method not only significantly improves computational efficiency but also provides a solid and scientific basis for engineering practice.

[0078] Based on the above Figure 1As can be seen from the implementation method, the present invention provides a method for determining the main controlling factors of the mechanical properties of heterogeneous dolomite. First, the microstructure data of the heterogeneous dolomite sample is acquired based on three-dimensional scanning technology, and its microstructure model is reconstructed accordingly. Then, the discrete element method is used to calculate the distribution of contact forces between particles in the model. Compared with traditional macroscopic models, this method can more realistically reflect the complexity of contact forces between particles. Next, numerical methods are used to simulate pore water flow in the pore network of the microstructure model. Compared with traditional fluid dynamics equations, this method can more accurately capture the complexity of the pore structure and the flow path of the fluid, thereby obtaining a more accurate pore pressure distribution. Furthermore, combined with the collected temperature data, the temperature change is calculated through thermodynamic equations, and its effect is analyzed. Based on the particle contact force distribution, pore pressure distribution, and temperature change effect, a multiphysics coupling model is constructed to simulate the mechanical response of heterogeneous dolomite. If the computational complexity of the multiphysics coupling model exceeds a preset threshold, dimensionality reduction techniques can be used to simplify the model, resulting in a simplified model. Stress and strain distribution data of heterogeneous dolomite under different environmental conditions were obtained by using a multiphysics coupling model. A simplified model was then used to extract mechanical characteristic parameters from the data, thereby determining the main controlling factors of the mechanical characteristics of heterogeneous dolomite under different environmental conditions. This approach employs a multiphysics coupling model in key parts to ensure accuracy, while simplifying the model in other parts to improve efficiency, thus significantly improving work efficiency while ensuring the accuracy of the results.

[0079] Furthermore, as a response to Figure 1 Further refinement and extension of the illustrated embodiment, this invention also provides another method for determining the main controlling factors of the mechanical properties of heterogeneous dolomite, such as... Figure 2 As shown, the specific steps are as follows:

[0080] 201. Using three-dimensional scanning technology, reconstruct the microstructure model of heterogeneous dolomite based on the microstructure data of heterogeneous dolomite samples.

[0081] 202. The contact force between particles is calculated using the discrete element method based on the contact information between particles in the microstructure model, and the particle contact force distribution is obtained.

[0082] 203. Numerical methods were used to simulate the flow of pore water in the pore network structure of the microstructure model to obtain the pore pressure distribution.

[0083] The implementation methods of steps 201 to 203 are the same as those of steps 101 to 103, and can achieve the same technical effect and solve the same technical problem, so they will not be repeated here.

[0084] 204. Use thermodynamic equations in conjunction with the collected temperature data to determine the amount of temperature change, and obtain the temperature change effect based on the amount of temperature change.

[0085] In this step, multiple temperature data points can be acquired using a preset temperature acquisition device. Then, using thermodynamic equations, the temperature change can be calculated based on the multiple temperature data points. Subsequently, the temperature change can be input into a microstructure model to simulate the effect of temperature change on interparticle contact force and pore structure, thereby obtaining data on particle contact force distribution and pore structure deformation under temperature change.

[0086] Subsequently, the particle contact force distribution and pore structure deformation data can be analyzed using a preset algorithm to calculate the changes in interparticle contact force and pore structure deformation caused by temperature changes. Based on the changes in interparticle contact force and pore structure deformation, the comprehensive effect of temperature changes on the mechanical properties of heterogeneous dolomite can be quantified, and a numerical table of temperature change effects can be generated.

[0087] In addition, after generating a numerical table of temperature change effects in the embodiments of the present invention, a machine learning algorithm can be used to fit the numerical table of temperature change effects to establish a predictive model between temperature change and microstructure change. Based on the predictive model, the microstructure change under different temperature conditions is predicted. Finally, a comprehensive report on the impact of temperature change on microstructure is generated based on the microstructure change under different temperature conditions.

[0088] 205. A multiphysics coupling model was constructed based on the particle contact force distribution, pore pressure distribution, and temperature change effect to simulate the mechanical response of heterogeneous dolomite.

[0089] In this step, the stress field distribution of heterogeneous dolomite can be calculated first using numerical methods based on the particle contact force distribution; at the same time, based on the pore pressure distribution and fluid seepage data, the pore pressure field distribution and seepage field distribution of heterogeneous dolomite can be calculated using fluid dynamics equations; and based on the temperature change effect, the temperature field distribution of heterogeneous dolomite can be calculated using heat conduction equations.

[0090] Then, the coupling field equations between the stress field, pore pressure field, seepage field and temperature field are solved by iterative method to obtain the multiphysics coupling model.

[0091] Furthermore, the mechanical parameters of the heterogeneous dolomite numerical model can be updated based on the solution of the multiphysics coupling model. Then, the support vector machine algorithm in machine learning is used to predict the mechanical response of the heterogeneous dolomite based on the numerical model of the heterogeneous dolomite with updated mechanical parameters. The mechanical response includes the stress-strain relationship. Finally, if the error between the prediction result and the experimental data exceeds a preset threshold, the mechanical parameters are adjusted and the multiphysics coupling field equations are solved again until the error meets the requirements, thus obtaining the final multiphysics coupling model.

[0092] 206. If the computational complexity of the multiphysics coupling model exceeds a preset threshold, the multiphysics coupling model is simplified using dimensionality reduction techniques to obtain a simplified model.

[0093] The implementation method of step 206 is the same as that of step 106, and can achieve the same technical effect and solve the same technical problem, so it will not be repeated here.

[0094] 207. Using a multiphysics coupling model, obtain stress and strain distribution data of heterogeneous dolomite under different environmental conditions.

[0095] In this step, different environmental conditions can be set according to research needs, such as temperature range (e.g., -20°C to 60°C), pressure range (e.g., 0.1 MPa to 50 MPa), and fluid flow conditions (e.g., hydrostatic pressure or dynamic water flow). These environmental conditions can be determined through experimental data or engineering requirements.

[0096] Subsequently, data on particle contact force distribution, pore pressure distribution, and temperature change effects can be input into the multiphysics coupling model to ensure that the initial conditions of the model are consistent with the actual environment. The multiphysics coupling model is then activated to simulate the mechanical response of heterogeneous dolomite under different environmental conditions.

[0097] The model iteratively solves for the coupling effects of interparticle contact forces, pore water flow, and temperature changes. During the calculation, the model's convergence is monitored in real time to ensure the accuracy of the results.

[0098] After the calculations are completed, stress and strain distribution data are extracted from the multiphysics coupling model. This data may include: stress distribution: maximum principal stress, minimum principal stress, stress concentration regions, etc.; strain distribution: volumetric strain, shear strain, strain concentration regions, etc.

[0099] Finally, the extracted stress and strain distribution data can be exported to a common format (such as CSV or VTK file) for easy subsequent analysis and visualization.

[0100] During visualization, post-processing software (such as ParaView or MATLAB) can be used to visualize the extracted stress and strain distribution data. This includes drawing stress cloud maps, strain cloud maps, and stress-strain curves to intuitively display the mechanical response of heterogeneous dolomite.

[0101] 208. Based on a simplified model, mechanical characteristic parameters of heterogeneous dolomite are extracted from stress distribution and strain distribution data, and the main controlling factors of mechanical characteristics of heterogeneous dolomite under different environmental conditions are determined based on the mechanical characteristic parameters.

[0102] In this step, stress and strain distribution data output from the multiphysics coupling model can be input into the simplified model. Then, using a preset algorithm in the simplified model, mechanical characteristic parameters of the heterogeneous dolomite are extracted from the stress and strain distribution data. For example, the elastic modulus is calculated using the stress-strain curve; Poisson's ratio is calculated using the ratio of transverse strain to longitudinal strain. The compressive strength is determined using the maximum principal stress value. The tensile strength is determined using the minimum principal stress value.

[0103] After obtaining the mechanical characteristic parameters, the variation patterns of these parameters can be analyzed based on temperature and humidity field data to determine the main controlling factors influencing the mechanical characteristics of heterogeneous dolomite. Then, principal component analysis is used to reduce the dimensionality of these main controlling factors, obtaining their weight values.

[0104] If the number of controlling factors exceeds a preset threshold, the random forest algorithm is used to rank the controlling factors by feature importance. Finally, based on the feature importance ranking results, target controlling factors with weight values ​​greater than the preset value are selected.

[0105] Finally, a quantitative relationship model between the main controlling factors and the mechanical characteristics of heterogeneous dolomite can be established through a multiple linear regression model.

[0106] 209. Based on the weight values ​​of each controlling factor, update the parameters of the heterogeneous dolomite numerical model to obtain the adjusted heterogeneous dolomite numerical model.

[0107] In this embodiment of the invention, after obtaining the weight values ​​of each main control factor, the parameters of the heterogeneous dolomite numerical model can be adjusted according to the weight values ​​of each main control factor to obtain the adjusted heterogeneous dolomite numerical model. In this way, more dolomite research can be carried out using the heterogeneous dolomite numerical model in the future.

[0108] Finally, a summary description of the embodiments of the present invention is provided:

[0109] This invention constructs a microstructural model of heterogeneous dolomite using three-dimensional scanning technology. Based on this model, the distribution of interparticle contact forces is calculated using the discrete element method, and the pore water flow is simulated numerically to obtain the pore pressure distribution. Then, based on the interparticle contact force distribution and pore pressure distribution, combined with the effect of temperature change, a multiphysics coupling model is constructed to simulate the mechanical response of heterogeneous dolomite under different environmental conditions, obtaining stress and strain distribution data.

[0110] When the computational complexity of a multiphysics coupling model is too high, dimensionality reduction techniques are used to simplify the model and improve computational efficiency. This allows for the extraction of mechanical characteristic parameters (such as elastic modulus and compressive strength) from the output data of the simplified multiphysics coupling model. These mechanical characteristic parameters then reveal the main controlling factors influencing the mechanical properties of heterogeneous dolomite. Subsequently, machine learning algorithms are used to establish a quantitative relationship model between these main controlling factors (such as temperature, pore pressure, and particle contact force) and the mechanical characteristic parameters.

[0111] The functions of each part and the relationships between them will be explained next:

[0112] Microstructure model: Constructed using 3D scanning technology, reflecting the microscopic characteristics of heterogeneous dolomite, such as grain shape and pore distribution. It provides fundamental data for multiphysics coupling models, supporting the calculation of interparticle contact forces and pore pressure distribution.

[0113] Multiphysics Coupled Model: Based on a microstructure model, this model couples particle contact forces, pore pressure, and temperature variation effects to simulate the mechanical response of heterogeneous dolomite. It outputs stress and strain distribution data, providing a basis for extracting mechanical characteristic parameters.

[0114] Simplified Model: When the computational complexity of a multiphysics coupling model is too high, dimensionality reduction techniques can be used to simplify the model and improve computational efficiency. The simplified model can still accurately reflect the mechanical behavior of heterogeneous dolomite and is suitable for large-scale simulations.

[0115] Numerical Model for Heterogeneous Dolomite: Based on stress-strain data output from a multiphysics coupled model, mechanical characteristic parameters (such as elastic modulus and compressive strength) are extracted. A quantitative relationship model between the main controlling factors and mechanical characteristic parameters is established using machine learning algorithms to predict the mechanical behavior of heterogeneous dolomite.

[0116] Mechanical characteristic parameters and controlling factors: Mechanical characteristic parameters (such as elastic modulus and Poisson's ratio) are quantitative indicators of the mechanical behavior of heterogeneous dolomite. Controlling factors (such as temperature and pore pressure) are key variables affecting mechanical characteristic parameters. By analyzing the weight values ​​of controlling factors, their degree of influence on the mechanical properties of heterogeneous dolomite can be determined.

[0117] More specifically:

[0118] Microstructure model: Representation: Usually presented as three-dimensional point cloud data or mesh model (such as STL, OBJ format), showing the microscopic characteristics of heterogeneous dolomite such as grain shape and pore distribution.

[0119] Function: It forms the basis of the research, providing microstructural data of heterogeneous dolomite and serving as input for subsequent multiphysics coupling models and numerical models of heterogeneous dolomite.

[0120] Multiphysics coupling model: Representation: Presented in the form of numerical simulation results, including stress distribution map, strain distribution map, pore pressure distribution map, etc.

[0121] Function: By simulating the mechanical behavior of heterogeneous dolomite, it outputs stress and strain distribution data, providing a basis for extracting mechanical characteristic parameters.

[0122] Numerical Model for Heterogeneous Dolomite: Presented as a mathematical or machine learning model, including quantitative relationships between mechanical characteristic parameters (such as elastic modulus and compressive strength) and key controlling factors (such as temperature and pore pressure). Purpose: To predict the mechanical behavior of heterogeneous dolomite under different environmental conditions, providing a scientific basis for engineering design and risk assessment.

[0123] It should be noted that the "microstructure model" or "numerical model of heterogeneous dolomite" mentioned in this invention are merely custom names defined in the context of this invention's embodiments.

[0124] Furthermore, as a response to the above Figure 1 In addition to the method shown, this embodiment of the invention also provides a device for determining the main controlling factors of the mechanical properties of heterogeneous dolomite, used for determining the main controlling factors of the above-mentioned mechanical properties. Figure 1 The method shown is implemented accordingly. This device embodiment corresponds to the foregoing method embodiment. For ease of reading, this device embodiment will not repeat the details of the foregoing method embodiment, but it should be clear that the device in this embodiment can implement all the contents of the foregoing method embodiment. Figure 3 As shown, the device includes:

[0125] The first model construction unit 301 is used to reconstruct the microstructure model of heterogeneous dolomite based on the microstructure data of heterogeneous dolomite samples using three-dimensional scanning technology.

[0126] The contact force calculation unit 302 is used to calculate the contact force between particles based on the contact information between particles in the microstructure model constructed by the first model construction unit 301 using the discrete element method, and to obtain the particle contact force distribution.

[0127] The pore pressure calculation unit 303 is used to simulate the flow of pore water in the pore network structure of the microstructure model constructed by the first model construction unit 301 using numerical methods, and to obtain the pore pressure distribution.

[0128] The temperature effect calculation unit 304 is used to determine the amount of temperature change by combining thermodynamic equations with the collected temperature data, and to obtain the temperature change effect based on the amount of temperature change.

[0129] The second model construction unit 305 is used to construct a multi-physics coupling model based on the particle contact force distribution calculated by the contact force calculation unit 302, the pore pressure distribution obtained by the pore pressure calculation unit 303, and the temperature change effect obtained by the temperature effect calculation unit 304, so as to simulate the mechanical response of heterogeneous dolomite.

[0130] The third model building unit 306 is used to simplify the multiphysics coupling model by using dimensionality reduction techniques to obtain a simplified model if the computational complexity of the multiphysics coupling model built by the second model building unit 305 exceeds a preset threshold.

[0131] The coupling model application unit 307 is used to obtain stress and strain distribution data of heterogeneous dolomite under different environmental conditions by using the multiphysics coupling model obtained by the second model construction unit 305.

[0132] The main control factor determination unit 308 is used to extract the mechanical characteristic parameters of heterogeneous dolomite from the stress distribution and strain distribution data obtained from the coupling model application unit 307 based on the simplified model obtained by the third model construction unit 306, and determine the main control factors of the mechanical characteristics of heterogeneous dolomite under different environmental conditions based on the mechanical characteristic parameters.

[0133] Furthermore, as a response to the above Figure 2 In addition to the method shown, this embodiment of the invention also provides another device for determining the main controlling factors of the mechanical properties of heterogeneous dolomite, used for determining the above-mentioned... Figure 2 The method shown is implemented accordingly. This device embodiment corresponds to the foregoing method embodiment. For ease of reading, this device embodiment will not repeat the details of the foregoing method embodiment, but it should be clear that the device in this embodiment can implement all the contents of the foregoing method embodiment. Figure 4 As shown, the device includes:

[0134] The first model construction unit 301 is used to reconstruct the microstructure model of heterogeneous dolomite based on the microstructure data of heterogeneous dolomite samples using three-dimensional scanning technology.

[0135] The contact force calculation unit 302 is used to calculate the contact force between particles based on the contact information between particles in the microstructure model constructed by the first model construction unit 301 using the discrete element method, and to obtain the particle contact force distribution.

[0136] The pore pressure calculation unit 303 is used to simulate the flow of pore water in the pore network structure of the microstructure model constructed by the first model construction unit 301 using numerical methods, and to obtain the pore pressure distribution.

[0137] The temperature effect calculation unit 304 is used to determine the amount of temperature change by combining thermodynamic equations with the collected temperature data, and to obtain the temperature change effect based on the amount of temperature change.

[0138] The second model construction unit 305 is used to construct a multi-physics coupling model based on the particle contact force distribution calculated by the contact force calculation unit 302, the pore pressure distribution obtained by the pore pressure calculation unit 303, and the temperature change effect obtained by the temperature effect calculation unit 304, so as to simulate the mechanical response of heterogeneous dolomite.

[0139] The third model building unit 306 is used to simplify the multiphysics coupling model by using dimensionality reduction techniques to obtain a simplified model if the computational complexity of the multiphysics coupling model built by the second model building unit 305 exceeds a preset threshold.

[0140] The coupling model application unit 307 is used to obtain stress and strain distribution data of heterogeneous dolomite under different environmental conditions by using the multiphysics coupling model obtained by the second model construction unit 305.

[0141] The main control factor determination unit 308 is used to extract the mechanical characteristic parameters of heterogeneous dolomite from the stress distribution and strain distribution data obtained from the coupling model application unit 307 based on the simplified model obtained by the third model construction unit 306, and determine the main control factors of the mechanical characteristics of heterogeneous dolomite under different environmental conditions based on the mechanical characteristic parameters.

[0142] In one optional implementation, the temperature effect calculation unit 304 is specifically used for:

[0143] Acquire multiple temperature data points collected using a preset temperature acquisition device;

[0144] Thermodynamic equations are used to calculate the temperature change based on multiple temperature data.

[0145] The temperature change is input into the microstructure model to simulate the effect of temperature change on interparticle contact force and pore structure, and to obtain the distribution of particle contact force and pore structure deformation data under temperature change.

[0146] The particle contact force distribution and pore structure deformation data are analyzed by a preset algorithm to calculate the change in particle contact force and pore structure deformation caused by temperature changes.

[0147] Based on the changes in interparticle contact force and pore structure deformation, the comprehensive effect of temperature change on the mechanical properties of heterogeneous dolomite is quantified, and a numerical table of temperature change effects is generated.

[0148] In one optional embodiment, after the temperature effect calculation unit 304 quantifies the comprehensive effect of temperature change on the mechanical properties of heterogeneous dolomite based on the change in interparticle contact force and the deformation of pore structure, and generates a temperature change effect numerical table, the device further includes a report generation unit 309, which is specifically used for:

[0149] Machine learning algorithms were used to fit the numerical table of temperature change effects in order to establish a predictive model between temperature change and microstructure change.

[0150] Based on the prediction model, the microstructure changes under different temperature conditions are predicted;

[0151] Based on the changes in microstructure under different temperature conditions, a comprehensive report on the impact of temperature changes on microstructure is generated.

[0152] In one optional implementation, the second model building unit 305 is specifically used for:

[0153] Based on the particle contact force distribution, the stress field distribution of heterogeneous dolomite was calculated using numerical methods.

[0154] Based on the pore pressure distribution and fluid seepage data, the pore pressure field distribution and seepage field distribution of heterogeneous dolomite are calculated using fluid dynamics equations.

[0155] Based on the aforementioned temperature change effect, the temperature field distribution of heterogeneous dolomite is calculated using the heat conduction equation.

[0156] The coupling field equations between the stress field, the pore pressure field, the seepage field, and the temperature field are solved by using an iterative method to obtain a multiphysics coupling model.

[0157] In one optional implementation, after the second model building unit 305 obtains the multiphysics coupling model, the device further includes a model tuning unit 310, which is specifically used for:

[0158] Based on the solution of the multiphysics coupling model, update the mechanical parameters of the numerical model of heterogeneous dolomite;

[0159] The mechanical response of heterogeneous dolomite, including stress-strain relationship, is predicted using a support vector machine (SVM) algorithm from the machine learning algorithm based on a numerical model of heterogeneous dolomite with updated mechanical parameters.

[0160] If the error between the predicted result and the experimental data exceeds the preset threshold, the mechanical parameters are adjusted and the multiphysics coupling field equations are solved again until the error meets the requirements, and the final multiphysics coupling model is obtained.

[0161] In one optional implementation, the main control factor determination unit 308 is specifically used for:

[0162] Based on temperature and humidity field data, the variation patterns of the mechanical characteristic parameters are analyzed to determine the main controlling factors affecting the mechanical characteristics of heterogeneous dolomite.

[0163] Principal component analysis was used to reduce the dimensionality of the main controlling factors and obtain the weight values ​​of each main controlling factor.

[0164] If the number of controlling factors exceeds a preset threshold, the random forest algorithm is used to rank the controlling factors by feature importance.

[0165] Based on the ranking of feature importance, target controlling factors with weight values ​​greater than preset values ​​are selected.

[0166] In one optional embodiment, the device further includes a numerical model building unit 311, which is used to update the parameters of the heterogeneous dolomite numerical model according to the weight values ​​of each main control factor, so as to obtain the adjusted heterogeneous dolomite numerical model.

[0167] Furthermore, embodiments of the present invention also provide a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the above-described... Figures 1-2 The method for determining the main controlling factors of the mechanical properties of heterogeneous dolomite described in the paper.

[0168] Furthermore, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the above-described... Figures 1-2 The method for determining the main controlling factors of the mechanical properties of heterogeneous dolomite described in the paper.

[0169] Furthermore, embodiments of the present invention also provide a computer program product, including a computer program, wherein the computer program, when executed by a processor, implements the above-described... Figures 1-2 The method for determining the main controlling factors of the mechanical properties of heterogeneous dolomite described in the paper.

[0170] In the above embodiments, the descriptions of each embodiment have different focuses. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions in other embodiments.

[0171] It is understood that the relevant features in the above methods and apparatus can be referenced interchangeably. Furthermore, the terms "first," "second," etc., in the above embodiments are used to distinguish between embodiments and do not represent the superiority or inferiority of any particular embodiment.

[0172] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0173] The algorithms and displays provided herein are not inherently related to any particular computer, virtual system, or other device. Various general-purpose systems can also be used in conjunction with the teachings herein. The required structure for constructing such systems is apparent from the above description. Furthermore, this invention is not directed to any particular programming language. It should be understood that the contents of the invention described herein can be implemented using various programming languages, and the above description of specific languages ​​is for the purpose of disclosing the best mode of implementation of the invention.

[0174] In addition, the memory may include non-permanent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0175] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0176] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxesFigure 1 A device that provides the functions specified in one or more boxes.

[0177] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0178] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

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

[0180] Memory may include non-persistent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0181] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information using any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.

[0182] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0183] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0184] The above are merely embodiments of this application and are not intended to limit the scope of this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of the claims of this application.

Claims

1. A method for determining the main controlling factors of the mechanical properties of heterogeneous dolomite, characterized in that, The method includes: Using 3D scanning technology, a microstructure model of heterogeneous dolomite was reconstructed based on the microstructure data of heterogeneous dolomite samples. The contact force between particles is calculated using the discrete element method based on the contact information between particles in the microstructure model, and the particle contact force distribution is obtained. The flow of pore water in the pore network structure of the microstructure model was simulated using numerical methods to obtain the pore pressure distribution. The temperature change is determined by combining thermodynamic equations with collected temperature data, and the temperature change effect is obtained based on the temperature change. A multiphysics coupling model is constructed based on the particle contact force distribution, the pore pressure distribution, and the temperature change effect to simulate the mechanical response of heterogeneous dolomite. If the computational complexity of the multiphysics coupling model exceeds a preset threshold, then the multiphysics coupling model is simplified using dimensionality reduction techniques to obtain a simplified model. Using the multiphysics coupling model, stress and strain distribution data of heterogeneous dolomite under different environmental conditions were obtained; Based on a simplified model, mechanical characteristic parameters of heterogeneous dolomite are extracted from stress and strain distribution data, and the main controlling factors of mechanical characteristics of heterogeneous dolomite under different environmental conditions are determined based on the mechanical characteristic parameters.

2. The method according to claim 1, characterized in that, The temperature change is determined using thermodynamic equations combined with collected temperature data, and the temperature change effect is obtained based on the temperature change, including: Acquire multiple temperature data points collected using a preset temperature acquisition device; Thermodynamic equations are used to calculate the temperature change based on multiple temperature data. The temperature change is input into the microstructure model to simulate the effect of temperature change on interparticle contact force and pore structure, and to obtain the distribution of particle contact force and pore structure deformation data under temperature change. The particle contact force distribution and pore structure deformation data are analyzed by a preset algorithm to calculate the change in particle contact force and pore structure deformation caused by temperature changes. Based on the changes in interparticle contact force and pore structure deformation, the comprehensive effect of temperature change on the mechanical properties of heterogeneous dolomite is quantified, and a numerical table of temperature change effects is generated.

3. The method according to claim 2, characterized in that, After quantifying the comprehensive effect of temperature change on the mechanical properties of heterogeneous dolomite based on the changes in interparticle contact force and pore structure deformation, and generating a numerical table of temperature change effects, the method further includes: Machine learning algorithms were used to fit the numerical table of temperature change effects in order to establish a predictive model between temperature change and microstructure change. Based on the prediction model, the microstructure changes under different temperature conditions are predicted; Based on the changes in microstructure under different temperature conditions, a comprehensive report on the impact of temperature changes on microstructure is generated.

4. The method according to claim 1, characterized in that, A multiphysics coupling model is constructed based on the particle contact force distribution, the pore pressure distribution, and the temperature change effect to simulate the mechanical response of heterogeneous dolomite, including: Based on the particle contact force distribution, the stress field distribution of heterogeneous dolomite was calculated using numerical methods. Based on the pore pressure distribution and fluid seepage data, the pore pressure field distribution and seepage field distribution of heterogeneous dolomite are calculated using fluid dynamics equations. Based on the aforementioned temperature change effect, the temperature field distribution of heterogeneous dolomite is calculated using the heat conduction equation. The coupling field equations between the stress field, the pore pressure field, the seepage field, and the temperature field are solved by using an iterative method to obtain a multiphysics coupling model.

5. The method according to claim 4, characterized in that, After obtaining the multiphysics coupling model, the method further includes: Based on the solution of the multiphysics coupling model, update the mechanical parameters of the numerical model of heterogeneous dolomite; The mechanical response of heterogeneous dolomite, including stress-strain relationship, is predicted using a support vector machine (SVM) algorithm from the machine learning algorithm based on a numerical model of heterogeneous dolomite with updated mechanical parameters. If the error between the predicted result and the experimental data exceeds the preset threshold, the mechanical parameters are adjusted and the multiphysics coupling field equations are solved again until the error meets the requirements, and the final multiphysics coupling model is obtained.

6. The method according to claim 1, characterized in that, Based on the aforementioned mechanical characteristic parameters, the main controlling factors of the mechanical characteristics of heterogeneous dolomite under different environmental conditions are determined, including: Based on temperature and humidity field data, the variation patterns of the mechanical characteristic parameters are analyzed to determine the main controlling factors affecting the mechanical characteristics of heterogeneous dolomite. Principal component analysis was used to reduce the dimensionality of the main controlling factors and obtain the weight values ​​of each main controlling factor. If the number of controlling factors exceeds a preset threshold, the random forest algorithm is used to rank the controlling factors by feature importance. Based on the ranking of feature importance, target controlling factors with weight values ​​greater than preset values ​​are selected.

7. The method according to claim 6, characterized in that, Based on the weight values ​​of each controlling factor, the parameters of the heterogeneous dolomite numerical model are updated to obtain the adjusted heterogeneous dolomite numerical model.

8. A device for determining the main controlling factors of the mechanical properties of heterogeneous dolomite, characterized in that, The device includes: The first model building unit is used to reconstruct the microstructure model of heterogeneous dolomite based on the microstructure data of heterogeneous dolomite samples using three-dimensional scanning technology. The contact force calculation unit is used to calculate the contact force between particles using the discrete element method based on the contact information between particles in the microstructure model constructed by the first model construction unit, and to obtain the particle contact force distribution. The pore pressure calculation unit is used to simulate the flow of pore water in the pore network structure of the microstructure model constructed by the first model construction unit using numerical methods, and to obtain the pore pressure distribution. The temperature effect calculation unit is used to determine the amount of temperature change by combining thermodynamic equations with the collected temperature data, and to obtain the temperature change effect based on the amount of temperature change. The second model construction unit is used to construct a multi-physics coupling model based on the particle contact force distribution calculated by the contact force calculation unit, the pore pressure distribution calculated by the pore pressure calculation unit, and the temperature change effect calculated by the temperature effect calculation unit, so as to simulate the mechanical response of heterogeneous dolomite. The third model building unit is used to simplify the multiphysics coupling model by using dimensionality reduction techniques to obtain a simplified model if the computational complexity of the multiphysics coupling model built by the second model building unit exceeds a preset threshold. The coupling model application unit is used to obtain stress and strain distribution data of heterogeneous dolomite under different environmental conditions by using the multiphysics coupling model obtained by the second model construction unit. The main control factor determination unit is used to extract the mechanical characteristic parameters of heterogeneous dolomite from the stress distribution and strain distribution data obtained from the coupling model application unit based on the simplified model obtained from the third model construction unit, and determine the main control factors of the mechanical characteristics of heterogeneous dolomite under different environmental conditions based on the mechanical characteristic parameters.

9. 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-7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the method described in any one of claims 1-7.

11. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the steps of the method described in any one of claims 1-7.

Images