Porous rock two-dimensional numerical model processing method and device, equipment and medium

Abstract

The invention discloses a porous rock two-dimensional numerical model processing method and device, equipment and a medium, and can be applied to the technical field of rock mechanics. The method comprises the following steps: importing a wafer radius and coordinate data of an initial rock sample generated by first software into second software, and obtaining user demand pore feature information including target porosity, target pore size distribution and target pore direction distribution probability; controlling the second software to generate a total wafer set according to the pore feature information required by the user and the wafer radius and coordinate data of the initial rock sample, importing the total pore wafer set into the first software, and then controlling the first software to process the initial rock sample according to the total wafer set, the discrete element model meeting the pore characteristic information required by the user is obtained, and then the mechanical behavior simulation of the porous rock microstructure characteristics is performed according to the discrete element model, so that the accuracy of porous rock simulation can be effectively improved, and the technical cost is reduced.

Description

Technical Field

[0001] The present application relates to the technical field of rock mechanics, and in particular to a method and device, equipment and medium for processing a two-dimensional numerical model of porous rocks. Background Art

[0002] Among the related technologies, in the field of geotechnical engineering, numerical simulation technology has become a key tool for exploring the mechanical behavior of porous media materials. The current mainstream numerical simulation methods include discrete element method (DEM), finite element method (FEM) and the coupling method of the two. Among them, the finite element method is restricted in its computational efficiency by factors such as the need for finite element meshing and the complex algorithm required to simulate the cracking of unit cells; discrete element numerical simulation usually relies on software to randomly generate pores when simulating pores, resulting in low simulation accuracy; and the application cost of CT scanning simulation technology is high and requires high technical mastery, which limits its application process.

[0003] In summary, the technical problems existing in the relevant technologies need to be improved. Summary of the invention

[0004] The main purpose of the embodiments of the present application is to propose a processing method, device, equipment and medium for a two-dimensional numerical model of porous rock, which can improve the accuracy of porous rock simulation and can be widely used.

[0005] To achieve the above object, one aspect of an embodiment of the present application provides a method for processing a two-dimensional numerical model of porous rock, the method comprising the following steps:

[0006] Controlling the first software to generate an initial rock sample, wherein the initial rock sample is formed by randomly and closely arranging discs;

[0007] Importing the disc radius and coordinate data of the initial rock sample into the second software;

[0008] Acquiring pore characteristic information required by a user, wherein the pore characteristic information required by the user includes a target porosity, a target pore size distribution, and a target pore direction distribution probability;

[0009] According to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample, the second software is controlled to generate a total set of discs corresponding to the pore characteristic information required by the user;

[0010] importing the total set of aperture discs into the first software;

[0011] Controlling the first software to process the initial rock sample according to the total set of discs to obtain a discrete element model;

[0012] The mechanical behavior of porous rock is simulated according to the discrete element model to obtain a stress-strain curve.

[0013] In some embodiments, the step of controlling the second software to generate a total set of discs corresponding to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample according to the pore characteristic information required by the user includes:

[0014] Determine the total number of pore discs according to the target porosity and the disc radius and coordinate data of the initial rock sample;

[0015] Calculating a target number of pores according to the target pore size distribution and the total number of pore discs;

[0016] The iteration step is executed, and when the number of iterations meets the iteration threshold, a total set of wafers corresponding to the pore characteristic information required by the user is generated. The iteration step includes:

[0017] Disc screening is performed according to the target pore number, the target pore direction distribution probability, and the disc radius and coordinate data of the initial rock sample.

[0018] In some embodiments, the disc screening according to the number of target pores, the probability of distribution of the target pore directions, and the disc radius and coordinate data of the initial rock sample includes:

[0019] Constructing a first matrix according to the disc radius and coordinate data of the initial rock sample;

[0020] generating a pore set according to the target pore number and the first matrix;

[0021] Selecting a number of discs from the pore set corresponding to the previous iteration step to form a disc set for the next iteration step;

[0022] Selecting a number of discs from the disc set of the next iteration step to form a central disc set;

[0023] According to the target pore direction distribution probability, randomly assign the pore growth direction to the central disc point in the central disc set, calculate the direction matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set, and calculate the distance matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set;

[0024] Screening the discs according to the direction matrix and the distance matrix to obtain a set of discs in each direction;

[0025] According to the disc sets in each direction, a disc set with pores in different directions is formed.

[0026] In some embodiments, the step of calculating the direction matrix of the center of the central disk and the center of each disk in the remaining non-porous disk set includes:

[0027] Calculating the azimuth cosine value of the center of the central disc and the center of each disc in the remaining non-porous disc set;

[0028] The direction matrix is ​​constructed according to all the azimuth cosine values; each row in the direction matrix corresponds to a central disk, and each column corresponds to a disk in the remaining non-porous disk set.

[0029] In some embodiments, the step of calculating the distance matrix between the center of the central disk and the centers of each disk in the remaining non-porous disk set includes:

[0030] Calculating the distance between the center of the central disc and the center of each disc in the remaining non-porous disc set;

[0031] The distance matrix is ​​constructed according to all the connection distances; each row in the distance matrix corresponds to a central disk, and each column corresponds to a disk in the remaining non-porous disk set.

[0032] In some embodiments, the processing of the initial rock sample according to the total set of discs to obtain a discrete element model includes:

[0033] The corresponding sphere in the initial rock sample is deleted according to the coordinate data in the total set of discs and the radius of the discs to obtain a discrete element model.

[0034] In some embodiments, simulating the mechanical behavior of porous rock according to the discrete element model to obtain a stress-strain curve includes:

[0035] Loading plates are arranged on the upper and lower surfaces of the discrete element model and loaded at a preset speed to obtain deformation and failure state data of the discrete element model under different stress conditions;

[0036] The stress-strain curve is generated based on the deformation and failure state data.

[0037] To achieve the above-mentioned purpose, another aspect of the embodiment of the present application provides a processing device for a two-dimensional numerical model of porous rock, the device comprising:

[0038] The first module is used to control the first software to generate an initial rock sample, wherein the initial rock sample is formed by randomly and closely arranging discs;

[0039] The second module is used to import the disc radius and coordinate data of the initial rock sample into the second software;

[0040] The third module is used to obtain the pore characteristic information required by the user, wherein the pore characteristic information required by the user includes the target porosity, the target pore size distribution and the target pore direction distribution probability;

[0041] The fourth module is used to control the second software to generate a total set of discs corresponding to the pore characteristic information required by the user according to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample;

[0042] A fifth module, for importing the total set of aperture discs into the first software;

[0043] A sixth module is used to control the first software to process the initial rock sample according to the total set of discs to obtain a discrete element model;

[0044] The seventh module is used to simulate the mechanical behavior of porous rock according to the discrete element model to obtain a stress-strain curve.

[0045] To achieve the above object, another aspect of an embodiment of the present application provides an electronic device, including:

[0046] at least one processor;

[0047] at least one memory for storing at least one program;

[0048] When the at least one program is executed by the at least one processor, the at least one processor implements the above method.

[0049] To achieve the above objective, another aspect of an embodiment of the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program implements the above method when executed by a processor.

[0050] The embodiments of the present application include at least the following beneficial effects: The present application provides a method and device, equipment and medium for processing a two-dimensional numerical model of porous rock. The scheme controls the first software to generate an initial rock sample formed by randomly and closely arranged discs, imports the disc radius and coordinate data of the initial rock sample into the second software, and then obtains the user-required pore characteristic information including the target porosity, target pore size distribution and target pore direction distribution probability. According to the user-required pore characteristic information and the disc radius and coordinate data of the initial rock sample, the second software is controlled to generate a total set of discs corresponding to the user-required pore characteristic information, and then after importing the total set of pore discs into the first software, the first software is controlled to process the initial rock sample according to the total set of discs to obtain a discrete element model that meets the pore characteristic information required by the user, and then simulates the mechanical behavior of the microstructure characteristics of porous rock according to the discrete element model, thereby effectively improving the accuracy of porous rock simulation, without relying on CT scanning technology, effectively reducing technical costs, and thus can be widely used. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a flow chart of a method for processing a two-dimensional numerical model of porous rock provided in an embodiment of the present application;

[0052] Figure 2 It is an application flow chart of the method for processing a two-dimensional numerical model of porous rocks provided in an embodiment of the present application;

[0053] Figure 3 is a pore size distribution histogram of a simulation provided in an embodiment of the present application;

[0054] Figure 4 It is a directional probability rose diagram for a simulation provided by an embodiment of the present application;

[0055] Figure 5 This is a schematic diagram showing the color separation of pore discs formed at different iteration steps during the simulation provided by the embodiment of the present application;

[0056] Figure 6 is a schematic diagram of a discrete element model of a porous medium material provided in an embodiment of the present application;

[0057] Figure 7 It is a failure diagram of the sample after loading provided in the embodiment of the present application;

[0058] Figure 8 is a crack distribution diagram of a sample after loading provided in an embodiment of the present application;

[0059] Fig. 9 is a stress-strain curve diagram output after loading provided in an embodiment of the present application;

[0060] Fig.10It is a structural schematic diagram of a processing device for a two-dimensional numerical model of porous rock provided in an embodiment of the present application;

[0061] Fig.11 It is a schematic diagram of the hardware structure of the electronic device provided in the embodiment of the present application. DETAILED DESCRIPTION

[0062] In order to make the purpose, technical solutions and advantages of the present application clearer, the present application is further described in detail below in conjunction with the accompanying drawings and examples. It should be understood that the specific embodiments described herein are only used to explain the present application and are not intended to limit the present application. When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the embodiments of the present application, they are only examples of devices and methods consistent with some aspects of the embodiments of the present application.

[0063] It is understood that the terms "first", "second", etc. used in this application can be used to describe various concepts in this article, but unless otherwise specified, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another concept. For example, without departing from the scope of the embodiment of the present application, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the words "if" and "if" as used herein can be interpreted as "at the time of" or "when" or "in response to determination".

[0064] The terms "at least one", "multiple", "each", "any", etc. used in this application, at least one includes one, two or more, multiple includes two or more, each refers to each of the corresponding multiple, and any refers to any one of the multiple.

[0065] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which this application belongs. The terms used herein are only for the purpose of describing the embodiments of this application and are not intended to limit this application.

[0066] Before describing the embodiments of the present application in detail, some nouns and terms involved in the embodiments of the present application are first described. The nouns and terms involved in the embodiments of the present application are subject to the following interpretations:

[0067] Rock: It is a solid aggregate with a stable shape, composed of one or more minerals and natural glass. It is the main material that constitutes the upper layer of the earth (the crust and at least the upper mantle). It is naturally produced and combined in a certain way. It is the product of various geological processes at a certain stage of the earth's development.

[0068] Porosity: refers to the spaces between particles or aggregates of particles in loose rocks. These spaces are an important condition for rocks to store groundwater and other fluids. The volume and distribution of pores have an important influence on the water storage capacity, permeability and other physical properties of rocks. For example, sedimentary rocks usually have the highest porosity, which makes them an important medium for oil and gas storage.

[0069] Disc: It refers to the name of the circular scattered particle unit when simulating two-dimensional rock samples using the discrete element method. Through these discs, the microstructure and particle composition of the rock can be represented in the discrete element software. By deleting the discs at specific positions in the sample, the pore shape, size and distribution of the rock can be simulated in the discrete element software.

[0070] The present application provides a method, device, equipment and medium for processing a two-dimensional numerical model of porous rock. The present application controls a first software to generate an initial rock sample formed by randomly and closely arranged discs, imports the disc radius and coordinate data of the initial rock sample into a second software, and then obtains user-required pore characteristic information including target porosity, target pore size distribution and target pore direction distribution probability. According to the user-required pore characteristic information and the disc radius and coordinate data of the initial rock sample, the second software is controlled to generate a total set of discs corresponding to the user-required pore characteristic information. After the total set of pore discs is imported into the first software, the first software is controlled to process the initial rock sample according to the total set of discs to obtain a discrete element model that meets the pore characteristic information required by the user, and then the mechanical behavior of the microstructure characteristics of the porous rock is simulated according to the discrete element model, thereby effectively improving the accuracy of the porous rock simulation, without relying on CT scanning technology, effectively reducing the technical cost, and thus being widely used.

[0071] The processing method of the two-dimensional numerical model of porous rock provided in the embodiment of the present application relates to the field of rock mechanics technology. The processing method of the two-dimensional numerical model of porous rock provided in the embodiment of the present application can be applied to a terminal, can also be applied to a server, and can also be software running in a terminal or a server. In some embodiments, the terminal can be a smart phone, a tablet computer, a laptop computer, a desktop computer, a smart speaker, a smart watch, and a car terminal, etc., but is not limited to this; the server side can be configured as an independent physical server, or it can be configured as a server cluster or distributed system composed of multiple physical servers, and can also be configured to provide cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms and other basic cloud computing services. The server can also be a node server in a blockchain network; the software can be an application of the processing method of the two-dimensional numerical model of porous rocks, etc., but is not limited to the above forms.

[0072] The present application can be used in many general or special computer system environments or configurations. For example: personal computers, server computers, handheld or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, etc. The present application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The present application can also be practiced in distributed computing environments, in which tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media including storage devices.

[0073] Figure 1 is an optional flow chart of a method for processing a two-dimensional numerical model of porous rock provided in an embodiment of the present application. Figure 1 The method may include but is not limited to steps S110 to S170:

[0074] Step S110, controlling the first software to generate an initial rock sample, wherein the initial rock sample is formed by randomly and closely arranging discs;

[0075] Step S120, importing the disc radius and coordinate data of the initial rock sample into the second software;

[0076] Step S130, obtaining pore characteristic information required by the user, wherein the pore characteristic information required by the user includes target porosity, target pore size distribution and target pore direction distribution probability;

[0077] Step S140, according to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample, controlling the second software to generate a total set of discs corresponding to the pore characteristic information required by the user;

[0078] Step S150, importing the total set of pore discs into the first software;

[0079] Step S160, controlling the first software to process the initial rock sample according to the total set of discs to obtain a discrete element model;

[0080] Step S170: simulate the mechanical behavior of porous rock according to the discrete element model to obtain a stress-strain curve.

[0081] It is understandable that if Figure 2As shown, in this embodiment, the corresponding data processing process can be performed by the first software and the second software respectively, wherein the first software can be PFC software and the second software can be MATLAB software. After the PFC software generates the initial rock sample composed of discs randomly and closely arranged, the radius and coordinate data of all discs in the initial rock sample are imported into the MATLAB software.

[0082] Specifically, when the user needs to construct a simulation model corresponding to specific pore characteristic information, the user's required pore characteristic information including the target porosity, target pore size distribution and target pore direction distribution probability can be input in the input interface of the MATLAB software. After obtaining the user's required pore characteristic information, the MATLAB software can generate a total set of discs corresponding to the user's required pore characteristic information based on the disc radius and coordinate data of the initial rock sample. It can be understood that the generation process of the total set of discs includes but is not limited to the following steps:

[0083] Step S210, determining the total number of pore discs according to the target porosity and the disc radius and coordinate data of the initial rock sample;

[0084] Step S220, calculating the target number of pores according to the target pore size distribution and the total number of pore discs;

[0085] Step S230, executing the iteration step, when the number of iterations meets the iteration threshold, generating a total set of discs corresponding to the pore characteristic information required by the user, wherein the iteration step includes: screening the discs according to the number of target pores, the probability of distribution of the target pore direction, and the disc radius and coordinate data of the initial rock sample.

[0086] In the embodiment of the present application, the process of disc screening according to the target pore number, the target pore direction distribution probability and the disc radius and coordinate data of the initial rock sample includes but is not limited to the following steps:

[0087] Step S310, constructing a first matrix according to the disc radius and coordinate data of the initial rock sample;

[0088] Step S320, generating a pore set according to the target pore number and the first matrix;

[0089] Step S330, selecting a number of discs from the pore set corresponding to the previous iteration step to form a disc set for the next iteration step;

[0090] Step S340, selecting a number of discs from the disc set of the next iteration step to form a central disc set;

[0091] Step S350: According to the target pore direction distribution probability, randomly assign pore growth directions to the central wafer points in the central wafer set, calculate the direction matrix of the central wafer center and the centers of each wafer in the remaining non-pore wafer set, and calculate the distance matrix of the central wafer center and the centers of each wafer in the remaining non-pore wafer set;

[0092] Step S360: Screen the wafers according to the direction matrix and the distance matrix to obtain wafer sets in each direction;

[0093] Step S370: Form wafer sets of different direction pores according to the wafer sets in each direction.

[0094] It can be understood that, as Figure 2 shown, after the MATLAB software obtains the wafer radius and coordinate data of the initial rock specimen, it generates the first matrix H corresponding to the wafer coordinates and wafer radii to provide corresponding data for subsequent steps. Specifically, in this embodiment, according to the first matrix H and the target porosity δ, the total number of pore wafers P required is determined, and then the target pore number N generated by different numbers of wafers is calculated i . Among them, i is up to the maximum pore number M specified by the user. After obtaining the target pore number N i , when i = 1, randomly select wafers to form the initial pore set A1 as the starting point for subsequent pore development. Then iteratively generate pore wafers.

[0095] Specifically, for the case where 1 < i ≤ M, randomly select i-1 wafers from the pore set A corresponding to the previous iteration step as the wafer set A i of the next iteration step. Subsequently, according to the direction distribution probability θ j , randomly assign the growth direction of the pore to S i *θ j central wafer points. Subsequently, form a direction matrix composed of the cosine values of the azimuth angles of the central wafer centers and the centers of the remaining non-pore wafer sets and a distance matrix composed of the connection distances. Then, through two screening processes of first selecting the wafer points that meet the specified direction in the direction matrix and then selecting the wafer with the smallest element value in the distance matrix, wafer sets B iθj in each direction are obtained, and they are combined into the wafer set B i for generating the i-th pore. This process is iteratively carried out until the specified maximum pore number M is reached, and the total wafer set Q of pores is formed, Q = A1 ∪ B1... ∪ B M. In this embodiment, the calculation process of the direction matrix can be constructed according to all the azimuth cosine values ​​after calculating the azimuth cosine values ​​between the center of the central disc and the center of each disc in the remaining non-porous disc set. Wherein, each row in the direction matrix corresponds to the central disc, and each column corresponds to the disc in the remaining non-porous disc set. The calculation process of the distance matrix can be constructed according to all the line distances after calculating the line distance between the center of the central disc and the center of each disc in the remaining non-porous disc set. Wherein, each row in the distance matrix corresponds to the central disc, and each column corresponds to the disc in the remaining non-porous disc set.

[0096] In an embodiment of the present application, after the total set of obtained discs is imported into the PFC software, the processor in the PFC software can delete the corresponding specified spheres in the initial rock sample according to the coordinates and diameter data of each point, thereby constructing a discrete element model containing an accurate pore structure.

[0097] In the embodiment of the present application, in order to simulate the mechanical behavior of porous rock in the actual environment, the present embodiment further sets and loads the model. Specifically, the present embodiment uses the parallel adhesive contact model to simulate the internal bonding strength of the porous medium material. The simulation results are as follows: Figure 3 , Figure 4 and Figure 5 In addition, the present embodiment can also set loading plates on the upper and lower surfaces of the discrete element model and load them at a preset speed to simulate the deformation and failure process of porous rock under different stress conditions. After the deformation and failure state data of the discrete element model under different stress conditions are obtained, a stress-strain curve is generated according to the deformation and failure state data for subsequent analysis and research. Figure 6 Taking the porous medium material with pores as an example for simulation, we can get Figure 7 The specimen failure diagram after loading is shown. Figure 8 The crack distribution diagram and Fig. 9 Schematic diagram of the stress-strain curve shown.

[0098] From the above content, it can be known that when the embodiment of the present application is used, the user only needs to input key parameters such as pore size distribution and probability in each direction, and then the simulation can be easily performed without relying on CT scanning technology, which effectively reduces the technical cost and can be widely used; and by using PFC software to generate an initial set of discs, and then combining MATLAB for complex data processing and analysis, and finally returning to PFC for model construction, the accuracy of porous rock simulation can be effectively improved. In addition, the present application simulates the internal bonding strength of porous medium materials by setting a parallel bonding contact model, which further enhances the physical authenticity of the numerical model. In addition, the method of the present application can also be applied to mining, geological engineering and other fields.

[0099] Reference Fig.10 The embodiment of the present application provides a device for processing a two-dimensional numerical model of porous rock, the device comprising:

[0100] The first module 1010 is used to control the first software to generate an initial rock sample, wherein the initial rock sample is formed by randomly and closely arranging discs;

[0101] The second module 1020 is used to import the disc radius and coordinate data of the initial rock sample into the second software;

[0102] The third module 1030 is used to obtain the pore characteristic information required by the user, wherein the pore characteristic information required by the user includes the target porosity, the target pore size distribution and the target pore direction distribution probability;

[0103] The fourth module 1040 is used to control the second software to generate a total set of discs corresponding to the pore characteristic information required by the user according to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample;

[0104] A fifth module 1050 is used to import the total set of aperture discs into the first software;

[0105] The sixth module 1060 is used to control the first software to process the initial rock sample according to the total set of discs to obtain a discrete element model;

[0106] The seventh module 1070 is used to simulate the mechanical behavior of porous rock according to the discrete element model to obtain the stress-strain curve.

[0107] It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0108] The embodiment of the present application also provides an electronic device, the electronic device includes a memory and a processor, the memory stores a computer program, and the processor implements the above-mentioned method for processing the two-dimensional numerical model of the hole rock when executing the computer program. The electronic device can be any intelligent terminal including a tablet computer, a car computer, etc.

[0109] It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0110] See also Fig.11 , Fig.11The hardware structure of an electronic device of another embodiment is illustrated, and the electronic device includes:

[0111] The processor 1110 may be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of the present application.

[0112] The memory 1120 can be implemented in the form of a read-only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM). The memory 1120 can store an operating system and other application programs. When the technical solution provided in the embodiment of this specification is implemented by software or firmware, the relevant program code is stored in the memory 1120, and the processor 1110 calls and executes the processing method of the two-dimensional numerical model of the hole rock in the embodiment of this application;

[0113] Input / output interface 1130, used to implement information input and output;

[0114] The communication interface 1140 is used to realize the communication interaction between the device and other devices. The communication can be realized through a wired manner (such as USB, network cable, etc.) or a wireless manner (such as mobile network, WI FI, Bluetooth, etc.);

[0115] bus 1150 , which transmits information between the various components of the device (e.g., processor 1110 , memory 1120 , input / output interface 1130 , and communication interface 1140 );

[0116] The processor 1110 , the memory 1120 , the input / output interface 1130 , and the communication interface 1140 are connected to each other in communication within the device via the bus 1150 .

[0117] An embodiment of the present application also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the processing method of the above-mentioned two-dimensional numerical model of the porous rock is implemented.

[0118] It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiments, the functions specifically implemented by the present storage medium embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0119] The memory, as a non-transient computer-readable storage medium, can be used to store non-transient software programs and non-transient computer executable programs. In addition, the memory may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some embodiments, the memory may optionally include a memory remotely disposed relative to the processor, and these remote memories may be connected to the processor via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0120] The embodiments described in the embodiments of the present application are intended to more clearly illustrate the technical solutions of the embodiments of the present application and do not constitute a limitation on the technical solutions provided in the embodiments of the present application. Those skilled in the art will appreciate that with the evolution of technology and the emergence of new application scenarios, the technical solutions provided in the embodiments of the present application are also applicable to similar technical problems.

[0121] Those skilled in the art will appreciate that the technical solutions shown in the figures do not constitute a limitation on the embodiments of the present application, and may include more or fewer steps than shown in the figures, or a combination of certain steps, or different steps.

[0122] The device embodiments described above are merely illustrative, and the units described as separate components may or may not be physically separated, that is, they may be located in one place or distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0123] Those skilled in the art will appreciate that all or some of the steps in the methods disclosed above, and the functional modules / units in the systems and devices may be implemented as software, firmware, hardware, or a suitable combination thereof.

[0124] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0125] It should be understood that in the present application, "at least one (item)" means one or more, and "plurality" means two or more. "And / or" is used to describe the association relationship of associated objects, indicating that three relationships may exist. For example, "A and / or B" can mean: only A exists, only B exists, and A and B exist at the same time, where A and B can be singular or plural. The character " / " generally indicates that the objects associated before and after are in an "or" relationship. "At least one of the following" or similar expressions refers to any combination of these items, including any combination of single or plural items. For example, at least one of a, b or c can mean: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, c can be single or multiple.

[0126] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the above units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0127] The units described above as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0128] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0129] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including multiple instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of various embodiments of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (Read-Only Memory, referred to as ROM), random access memory (Random Access Memory, referred to as RAM), disk or optical disk and other media that can store programs.

[0130] The preferred embodiments of the present invention are described above with reference to the accompanying drawings, but the scope of the rights of the present invention is not limited thereto. Any modification, equivalent substitution and improvement made by a person skilled in the art without departing from the scope and essence of the present invention should be within the scope of the rights of the present invention.

Claims

1. A method for processing a two-dimensional numerical model of porous rock, characterized in that: The method comprises the following steps: Controlling the first software to generate an initial rock sample, wherein the initial rock sample is formed by randomly and closely arranging discs; Importing the disc radius and coordinate data of the initial rock sample into the second software; Acquiring pore characteristic information required by a user, wherein the pore characteristic information required by the user includes a target porosity, a target pore size distribution, and a target pore direction distribution probability; According to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample, the second software is controlled to generate a total set of discs corresponding to the pore characteristic information required by the user; importing the total set of aperture discs into the first software; Controlling the first software to process the initial rock sample according to the total set of discs to obtain a discrete element model; The mechanical behavior of porous rock is simulated according to the discrete element model to obtain a stress-strain curve.

2. The method according to claim 1, characterized in that The step of controlling the second software to generate a total set of discs corresponding to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample according to the pore characteristic information required by the user includes: Determine the total number of pore discs according to the target porosity and the disc radius and coordinate data of the initial rock sample; Calculating a target number of pores according to the target pore size distribution and the total number of pore discs; The iteration step is executed, and when the number of iterations meets the iteration threshold, a total set of wafers corresponding to the pore characteristic information required by the user is generated. The iteration step includes: Disc screening is performed according to the target pore number, the target pore direction distribution probability, and the disc radius and coordinate data of the initial rock sample.

3. The method according to claim 2, characterized in that The disc screening is performed according to the number of target pores, the probability of distribution of the target pore directions, and the disc radius and coordinate data of the initial rock sample, including: Constructing a first matrix according to the disc radius and coordinate data of the initial rock sample; generating a pore set according to the target pore number and the first matrix; Selecting a number of discs from the pore set corresponding to the previous iteration step to form a disc set for the next iteration step; Selecting a number of discs from the disc set of the next iteration step to form a central disc set; According to the target pore direction distribution probability, randomly assign the pore growth direction to the central disc point in the central disc set, calculate the direction matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set, and calculate the distance matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set; Screening the discs according to the direction matrix and the distance matrix to obtain a set of discs in each direction; According to the disc sets in each direction, a disc set with pores in different directions is formed.

4. The method according to claim 3, characterized in that The calculation of the direction matrix of the center of the central disk and the center of each disk in the remaining non-porous disk set includes: Calculating the azimuth cosine value of the center of the central disc and the center of each disc in the remaining non-porous disc set; The direction matrix is ​​constructed according to all the azimuth cosine values; each row in the direction matrix corresponds to a central disk, and each column corresponds to a disk in the remaining non-porous disk set.

5. The method according to claim 3, characterized in that: The calculation of the distance matrix between the center of the central disk and the centers of each disk in the remaining non-porous disk set includes: Calculating the distance between the center of the central disc and the center of each disc in the remaining non-porous disc set; The distance matrix is ​​constructed according to all the connection distances; each row in the distance matrix corresponds to a central disk, and each column corresponds to a disk in the remaining non-porous disk set.

6. The method according to claim 1, characterized in that The initial rock sample is processed according to the total set of discs to obtain a discrete element model, including: The corresponding sphere in the initial rock sample is deleted according to the coordinate data in the total set of discs and the radius of the discs to obtain a discrete element model.

7. The method according to claim 1, characterized in that The mechanical behavior simulation of porous rock is performed according to the discrete element model to obtain a stress-strain curve, including: Loading plates are arranged on the upper and lower surfaces of the discrete element model and loaded at a preset speed to obtain deformation and failure state data of the discrete element model under different stress conditions; The stress-strain curve is generated based on the deformation and failure state data.

8. A processing device for a two-dimensional numerical model of porous rocks, characterized in that: The device comprises: The first module is used to control the first software to generate an initial rock sample, wherein the initial rock sample is formed by randomly and closely arranging discs; The second module is used to import the disc radius and coordinate data of the initial rock sample into the second software; The third module is used to obtain the pore characteristic information required by the user, wherein the pore characteristic information required by the user includes the target porosity, the target pore size distribution and the target pore direction distribution probability; The fourth module is used to control the second software to generate a total set of discs corresponding to the pore characteristic information required by the user according to the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample; A fifth module, for importing the total set of aperture discs into the first software; A sixth module is used to control the first software to process the initial rock sample according to the total set of discs to obtain a discrete element model; The seventh module is used to simulate the mechanical behavior of porous rock according to the discrete element model to obtain a stress-strain curve.

9. An electronic device, characterized in that: include: at least one processor; at least one memory for storing at least one program; When the at least one program is executed by the at least one processor, the at least one processor implements the method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images