Method, device, equipment and medium for processing two-dimensional numerical model of porous rock

Abstract

The present application discloses a method, device, equipment and medium for processing a two-dimensional numerical model of porous rock, which can be applied to the field of rock mechanics technology. The present application imports the disc radius and coordinate data of the initial rock sample generated by the first software into the second software, and 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 disc set. After the total pore disc set is imported into the first software, the first software is controlled to process the initial rock sample according to the total disc set to obtain a discrete element model that meets the user's pore characteristic information. The mechanical behavior simulation of the microstructural characteristics of the porous rock is then performed based on the discrete element model, thereby effectively improving the accuracy of the porous rock simulation and reducing the technical cost.

Description

Technical Field

[0001] The present application relates to the field of rock mechanics technology, and in particular to a method and apparatus, 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 mesh division and the complex algorithm required to simulate unit body cracking; 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 method, device, equipment and medium for processing 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-mentioned purpose, 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 user-required pore characteristic information, wherein the user-required pore characteristic information includes target porosity, target pore size distribution, and 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, 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] Determining the total number of pore discs based on the target porosity and the disc radius and coordinate data of the initial rock sample;

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

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

[0017] 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.

[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 based on the disc radius and coordinate data of the initial rock sample;

[0020] generating a pore set according to the target number of pores 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, the pore growth direction is randomly assigned to the central disc point in the central disc set, and the direction matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set is calculated, as well as the distance matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set.

[0024] Performing disc screening 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 set of non-aperture disks includes:

[0027] Calculating the cosine value of the azimuth angle between 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, calculating the distance matrix between the center of the central disc and the centers of the discs in the remaining set of non-aperture discs 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 line distances; each row in the distance matrix corresponds to a central disk, and each column corresponds to a disk in the remaining non-pore disk set.

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

[0033] The corresponding spheres in the initial rock sample are deleted according to the coordinate data and the disc radius in the total disc set 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 present application provides a device for processing 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 arranged 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 target porosity, target pore size distribution and target pore direction distribution probability;

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

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

[0043] A sixth module is configured 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 objectives, another aspect 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-mentioned purpose, 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-mentioned 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, 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, and then imports the disc radius and coordinate data of the initial rock sample into the second software. Then, after obtaining the user-required pore characteristic information including the target porosity, target pore size distribution and target pore direction distribution probability, the second software is controlled to generate a total set of discs corresponding to the user-required pore characteristic information according to the user-required pore characteristic information and the disc radius and coordinate data of the initial rock sample. 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 user-required pore characteristic information, and then simulates the mechanical behavior of the microstructural characteristics of the porous rock based on the discrete element model, thereby effectively improving the accuracy of the 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 This 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 This is an application flow chart of the method for processing a two-dimensional numerical model of porous rock 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 of a simulation provided by an embodiment of the present application;

[0055] Figure 5 This is a schematic diagram showing the color separation of the 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 This is a diagram of the sample failure after loading provided in the embodiment of the present application;

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

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

[0060] Figure 10Schematic diagram of the structure of a processing device for a two-dimensional numerical model of porous rock provided in an embodiment of the present application;

[0061] Figure 11 This is a schematic diagram of the hardware structure of the electronic device provided in an 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 with reference to 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 merely examples of devices and methods consistent with some aspects of the embodiments of the present application.

[0063] It will be understood that the terms "first", "second", etc. used in this application may be used herein to describe various concepts, but unless otherwise specified, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another. For example, without departing from the scope of the embodiments 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 may be interpreted as "at the time of" or "when" or "in response to determining".

[0064] The terms "at least one", "plurality", "each", "any", etc. used in this application include "at least one", "two" or more, "plurality" or "each", "any" or "any one", "each" or "any one" as used herein.

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

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

[0067] Rock: A stable, solid aggregate composed of one or more minerals and natural glass. It is the primary material that makes up the Earth's upper layers (the crust and at least the upper mantle). It is naturally occurring, formed in a specific way, and is the product of various geological processes at a certain stage in the Earth's development.

[0068] Porosity refers to the spaces between particles or aggregates of particles in loose rock. These spaces are crucial for rock's ability to store groundwater and other fluids. The volume and distribution of pores significantly influence a rock's water storage capacity, permeability, and other physical properties. For example, sedimentary rocks typically have the highest porosity, making them ideal media for oil and gas storage.

[0069] Discs: This is a name for the circular, scattered particles used in the discrete element method (DEM) simulation of two-dimensional rock specimens. These discs can be used to represent the rock's microstructure and particle composition within discrete element software. By deleting discs at specific locations within the specimen, the pore shape, size, and distribution of the rock can be simulated within the DEM software.

[0070] The embodiment of the present application provides a processing method, device, equipment and medium for a two-dimensional numerical model of porous rock. The present application controls the first software to generate an initial rock sample formed by randomly and closely arranged discs, and then imports the disc radius and coordinate data of the initial rock sample into the second software. Then, after obtaining the user-required pore characteristic information including the target porosity, target pore size distribution and target pore direction distribution probability, the second software is controlled to generate a total set of discs corresponding to the user-required pore characteristic information according to the user-required pore characteristic information and the disc radius and coordinate data of the initial rock sample. 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 user-required pore characteristic information, and then simulates the mechanical behavior of the microstructural characteristics of porous rock based on 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.

[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 be applied to a server, or can 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 can be configured as a server cluster or distributed system composed of multiple physical servers, and can also be configured as a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. 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 rock, 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, and the like. 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, and the like 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 via 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 This is an optional flow chart of the 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: Control 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: Acquire user-required pore characteristic information, wherein the user-required pore characteristic information includes target porosity, target pore size distribution, and target pore direction distribution probability;

[0077] Step S140: Based on the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample, control 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: Control 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, corresponding data processing processes can be performed by first software and second software, respectively. The first software can be PFC software, and the second software can be MATLAB software. In this embodiment, after generating an initial rock sample composed of randomly and closely arranged discs using PFC software, the radius and coordinate data of all discs in the initial rock sample are imported into MATLAB software.

[0082] Specifically, when a 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 into the MATLAB software input interface. 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 is understood that the process of generating 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, a total set of discs corresponding to the pore characteristic information required by the user is generated. 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 based on the target pore number, target pore direction distribution probability, and 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 based on 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: Select several discs from the pore set corresponding to the previous iteration step to form a disc set for the next iteration step;

[0090] Step S340: Select several 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 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 P of required pore wafers 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 of 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 each wafer in the remaining non-pore wafer set 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 performed until the specified maximum pore number M is reached, and the total wafer set Q of pores is formed as Q = A1 ∪ B1... ∪ B MIn this embodiment, the direction matrix can be calculated by calculating the azimuth cosine values of the center of the central disc and the center of each disc in the remaining non-porous disc set, and then constructing the direction matrix based on all the azimuth cosine values. In which, each row in the direction matrix corresponds to the central disc, and each column corresponds to the discs in the remaining non-porous disc set. The distance matrix can be calculated by calculating the line distance between the center of the central disc and the center of each disc in the remaining non-porous disc set, and then constructing the distance matrix based on all the line distances. In which, each row in the distance matrix corresponds to the central disc, and each column corresponds to the discs in the remaining non-porous disc set.

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

[0097] In this embodiment of the present application, in order to simulate the mechanical behavior of porous rocks in a real environment, this embodiment further sets and loads the model. Specifically, this embodiment uses a parallel adhesive contact model to simulate the internal bonding strength of porous media materials. The simulation results are as follows: Figure 3 、 Figure 4 and Figure 5 In addition, this 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 obtaining the deformation and failure state data of the discrete element model under different stress conditions, a stress-strain curve is generated based on the deformation and failure state data for subsequent analysis and research. For example, Figure 6 Taking the discrete element model of porous medium material with pores as an example, we can get Figure 7 The specimen failure diagram after loading is shown. Figure 8 The crack distribution map shown and Figure 9 Schematic diagram of the stress-strain curve shown.

[0098] From the above content, it can be seen 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 to easily perform simulation, without relying on CT scanning technology, effectively reducing technical costs, and thus can be widely used; and by using PFC software to generate an initial set of discs, 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 media materials by setting a parallel bonding contact model, further enhancing the physical reality 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 Figure 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 arranged 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 target porosity, target pore size distribution and target pore direction distribution probability;

[0103] The fourth module 1040 is configured to control 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;

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

[0105] The sixth module 1060 is configured 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 a 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 present application also provides an electronic device comprising a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the aforementioned method for processing a two-dimensional numerical model of a borehole rock. The electronic device can be any smart terminal, including a tablet computer and an in-vehicle computer.

[0109] It can be understood that the contents of the above method embodiments are 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 Figure 11 , Figure 11The hardware structure of an electronic device according to another embodiment is shown. The electronic device includes:

[0111] The processor 1110 may be implemented using a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is configured 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 solutions provided in the embodiments of this specification are implemented through software or firmware, the relevant program codes are stored in the memory 1120, and the processor 1110 calls and executes the method for processing the two-dimensional numerical model of the pore rock in the embodiments of this application.

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

[0114] Communication interface 1140, used to implement communication interaction between this device and other devices, which can be achieved through wired means (such as USB, network cable, etc.) or wireless means (such as mobile network, Wi-Fi, Bluetooth, etc.);

[0115] bus 1150 , which transmits information between 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 communicatively connected to each other within the device via a bus 1150 .

[0117] An embodiment of the present application further provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the above-mentioned method for processing the two-dimensional numerical model of the hole rock.

[0118] It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiment, the functions specifically implemented by the present storage medium embodiment 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 arranged 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 this application are intended to more clearly illustrate the technical solutions of the embodiments of this application and do not constitute a limitation on the technical solutions provided by the embodiments of this 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 this application are also applicable to similar technical problems.

[0121] Those skilled in the art will understand 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. The units described as separate components may or may not be physically separate, i.e., they may be located in one place or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of this embodiment.

[0123] Those skilled in the art will appreciate that all or some of the steps in the methods, systems, and functional modules / units in the devices disclosed above may be implemented as software, firmware, hardware, or appropriate combinations 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 sequential order. 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 variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0125] It should be understood that in this 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 previous and next associated objects are in an "or" relationship. "At least one of the following items" or similar expressions refers to any combination of these items, including any combination of single items 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 this 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 merely schematic. For example, the division of the above-mentioned units is only a logical function division. In actual implementation, there may be other division methods, 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 separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0128] In addition, the functional units in the various embodiments of the present application may be integrated into a single processing unit, or each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or 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 for enabling 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 (ROM), random access memory (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 are not intended to limit the scope of the present invention. Any modifications, equivalent substitutions, and improvements made by those skilled in the art without departing from the scope and essence of the present invention should be within the scope 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 user-required pore characteristic information, wherein the user-required pore characteristic information includes target porosity, target pore size distribution, and 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; Simulating the mechanical behavior of porous rock according to the discrete element model to obtain a stress-strain curve; 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: Determining the total number of pore discs based on the target porosity and the disc radius and coordinate data of the initial rock sample; Calculating a target number of pores based on the target pore size distribution and the total number of pore discs; The iterative step is executed. When the number of iterations meets the iteration threshold, a total set of wafers corresponding to the pore feature information required by the user is generated. The iterative step includes: Performing 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; 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: Constructing a first matrix based on the disc radius and coordinate data of the initial rock sample; generating a pore set according to the target number of pores 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, the pore growth direction is randomly assigned to the central disc point in the central disc set, and the direction matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set is calculated, as well as the distance matrix between the center of the central disc and the center of each disc in the remaining non-porous disc set. Performing disc screening 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.

2. The method according to claim 1, characterized in that The calculation of the direction matrix of the center of the central disk and the centers of each disk in the remaining non-porous disk set includes: Calculating the cosine value of the azimuth angle between 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.

3. The method according to claim 2, characterized in that The calculation of the distance matrix between the center of the central disc and the centers of the discs in the remaining non-porous disc 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 line distances; each row in the distance matrix corresponds to a central disk, and each column corresponds to a disk in the remaining non-pore disk set.

4. The method according to claim 1, wherein The processing of the initial rock sample according to the total set of discs to obtain a discrete element model includes: The corresponding spheres in the initial rock sample are deleted according to the coordinate data and the disc radius in the total disc set to obtain a discrete element model.

5. The method according to claim 1, wherein The mechanical behavior simulation of porous rock according to the discrete element model to obtain a stress-strain curve includes: 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.

6. A device for processing a two-dimensional numerical model of porous rock applied to the method according to any one of claims 1 to 5, 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 arranged 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 target porosity, target pore size distribution and target pore direction distribution probability; A fourth module is configured to control the second software to generate a total set of discs corresponding to the pore characteristic information required by the user, based on the pore characteristic information required by the user and the disc radius and coordinate data of the initial rock sample; A fifth module is used to import the total set of aperture discs into the first software; A sixth module is configured 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.

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

8. 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 5 is implemented.

Images