Visual experiment device and method for microscopic imbibition rate of shale reservoir

Abstract

The invention discloses a visual experiment device and method for the microscopic imbibition rate of a shale reservoir, and relates to the technical field of imbibition numerical simulation. In the visual experiment device, firstly, the pore radius of each pore of a target core column is scanned through a core pore radius scanning module, and a pore scanning result is obtained; calculating to obtain a pore radius distribution frequency curve so as to generate a core column finite element model; then the imbibition physical simulation experiment module is controlled, after simulation processing is conducted on the target core column, a static imbibition physical simulation experiment is conducted, an imbibition rate curve is obtained through calculation, and finally in the core simulation imaging module, a static imbibition numerical simulation experiment is conducted on the core column finite element model based on the imbibition rate curve. And the imbibition process in each pore is visually displayed. The scheme provided by the invention can be used for observing the recovery ratio and the imbibition rate of each pore of the core column in an imbibition physical simulation experiment, and is convenient for analyzing the imbibition influence factors and rules.

Description

Technical Field

[0001] The present application relates to the technical field of imbibition numerical simulation, and in particular to a visualization experimental device and method for the microscopic imbibition rate of a shale reservoir. Background Art

[0002] Many scholars at home and abroad have derived static imbibition physical simulation methods through the study of the imbibition laws of fractured reservoirs or low-permeability and ultra-low-permeability reservoirs. However, there is currently no visual research on the imbibition rate of the imbibition experiment. The current understanding of the imbibition rate only exists in the xy-axis imbibition rate curve drawn through the physical simulation experiment. It is impossible to conduct a stage-by-stage analysis of the imbibition state of the core column, and it is impossible to fully understand the contribution of the core pores of each pore radius to the imbibition rate at any imbibition rate, and it is difficult to fully understand the participation of the core pores of each pore radius in the entire imbibition process. To this end, the present application designs a visualization device and method for the microscopic imbibition rate of shale reservoirs, which can be used to study the relationship between the imbibition recovery degree and the imbibition rate of pores of each pore radius at different times. Summary of the invention

[0003] The purpose of this application is to provide a visualization experimental device and method for the microscopic imbibition rate of shale reservoirs, which can intuitively display the relationship between the pore imbibition recovery degree and the imbibition rate of each pore radius at different time points in the imbibition physical simulation experiment.

[0004] To achieve the above objectives, this application provides the following solutions: In the first aspect, the present application provides a visualization experimental device for the microscopic imbibition rate of a shale reservoir, comprising: a core pore radius scanning module, an imbibition physical simulation experiment module, a core simulation imaging module, and a data solution and process control module; the data solution and process control module is used to: control the core pore radius scanning module to scan the pore radius of each pore of a target core column to obtain a pore scanning result of the target core column, and solve the pore radius distribution frequency curve of the target core column according to the pore scanning result of the target core column; the pore radius distribution frequency curve is used to characterize the pores with different pore radii in the target core column. The frequency of occurrence in the simulation; control the imbibition physical simulation experiment module, simulate the target core column, conduct static imbibition physical simulation experiment on the simulated target core column, and solve the imbibition rate curve according to the experimental results of the static imbibition physical simulation experiment; the imbibition rate curve is used to characterize the relationship between the oil production rate per unit area and time; control the core simulation imaging module, generate the core column finite element model according to the pore radius distribution frequency curve, and conduct static imbibition numerical simulation experiment on the core column finite element model based on the imbibition rate curve, and visualize the imbibition process in each pore of the core column finite element model.

[0005] Optionally, in the core pore radius scanning module, the target core column is scanned by a CT scanner to obtain pore radius data within the macropore range of 1000nm~100000nm, and a high-speed mercury injection experiment is performed on the target core column to obtain pore radius data within the range of mesopores and macropores from 2nm~50000nm; the pore scanning result of the target core column includes pore radius data within the macropore range of 1000nm~100000nm and pore radius data within the range of mesopores and macropores from 2nm~50000nm; the data solution and process control module solves the pore radius distribution frequency curve of the target core column according to the pore radius data within the range of mesopores and macropores from 2nm~50000nm and the pore radius data within the range of macropores from 50000nm~100000nm.

[0006] Optionally, the imbibition physical simulation experiment module conducts a static imbibition physical simulation experiment on the treated target core column using a volumetric method based on an imbibition bottle; the treatment of the target core column includes introducing simulated oil into the pores of the target core column after vacuum drying and simulating reservoir temperature aging for a preset time to simulate the core column in a real reservoir environment.

[0007] Optionally, the data solution and process control module reads the oil output of the static imbibition physical simulation experiment at preset time intervals to calculate the imbibition rate at the current moment, and plots the imbibition rate at each moment of the static imbibition physical simulation experiment to obtain the imbibition rate curve of the target core column.

[0008] Optionally, the core simulation imaging module is used to: determine the proportion of pores with different pore radii in the target core column according to the pore radius distribution frequency curve; generate a core column finite element model according to the proportion of pores with different pore radii in the target core column; the core column finite element model and the target core column have exactly the same pore radius distribution frequency curve; fill the pores of the core column finite element model with simulated oil and set water at both ends of the pores; perform a static imbibition numerical simulation experiment on the core column finite element model according to the imbibition rate curve, and visualize the imbibition process in each pore of the core column finite element model.

[0009] Optionally, the core simulation imaging module is also used to: model pores with different pore radii respectively to obtain several pore models corresponding to different pore radii; place several pore models vertically, and draw a critical oil yield curve between pore radius and wetting angle through static imbibition numerical simulation experiments at different pore radii under a predetermined wetting angle and static imbibition numerical simulation experiments at different wetting angles under a predetermined pore radius; the abscissa and ordinate of the critical oil yield curve are the pore radius and the wetting angle, respectively, and any point on the critical oil yield curve indicates that under the condition of this pore radius and this wetting angle, the simulated oil in the pore model begins to displace in the pore model due to the action of gravity.

[0010] In a second aspect, the present application also provides a visualization experimental method for the microscopic imbibition rate of a shale reservoir, comprising the following steps: The pore radius of each pore of the acquired target core column is scanned to obtain the pore scanning result of the target core column.

[0011] According to the pore scanning results of the target core column, the pore radius distribution frequency curve of the target core column is calculated; the pore radius distribution frequency curve is used to characterize the occurrence frequency of pores with different pore radii in the target core column.

[0012] The target core column is simulated, and a static imbibition physical simulation experiment is carried out on the simulated target core column.

[0013] According to the experimental results of static imbibition physical simulation experiment, the imbibition rate curve is calculated; the imbibition rate curve is used to characterize the relationship between oil production rate per unit area and time.

[0014] According to the pore radius distribution frequency curve, a core column finite element model is generated.

[0015] Based on the imbibition rate curve, a static imbibition numerical simulation experiment was carried out on the core column finite element model to visualize the imbibition process in each pore of the core column finite element model.

[0016] Optionally, the pore scanning result of the target core column includes: pore radius data within the range of macropores from 1000nm to 100000nm and pore radius data within the range of mesopores and macropores from 2nm to 50000nm; scanning the pore radius of each pore of the acquired target core column specifically includes the following steps: The target core column was scanned by a CT scanner to obtain the pore radius data within the macropore range of 1000nm~100000nm.

[0017] A mercury injection experiment was carried out on the target core column using a high-speed mercury injection scanner to obtain pore radius data within the range of 2nm~50000nm for mesopores and macropores.

[0018] Optionally, generating a core column finite element model according to the pore radius distribution frequency curve specifically includes the following steps: According to the pore radius distribution frequency curve, the proportion of pores with different pore radii in the target core column is determined.

[0019] According to the proportion of pores with different pore radii in the target core column, a finite element model of the core column is generated; the finite element model of the core column and the target core column have exactly the same pore radius distribution frequency curve.

[0020] Optionally, simulated oil is filled in the pores of the core column finite element model and water is set at both ends of the pores, and a static imbibition numerical simulation experiment is performed on the core column finite element model according to the imbibition rate curve to visualize the imbibition process in each pore of the core column finite element model.

[0021] According to the specific embodiments provided in this application, this application discloses the following technical effects: The present application provides a visualization experimental device and method for the microscopic imbibition rate of a shale reservoir. In the visualization experimental device, the pore radius of each pore of a target core column is scanned by controlling a core pore radius scanning module to obtain a pore scanning result of the target core column, and a pore radius distribution frequency curve of the target core column is calculated based on the result to generate a finite element model of the core column; then, the imbibition physical simulation experimental module is controlled to simulate the target core column, and a static imbibition physical simulation experiment is performed on the simulated target core column to obtain an imbibition rate curve; finally, in a core simulation imaging module, a static imbibition numerical simulation experiment is performed on the finite element model of the core column based on the imbibition rate curve to visualize the imbibition process in each pore of the finite element model of the core column. The above-mentioned scheme provided in the present application can be used to visualize the imbibition process of the imbibition physical simulation experiment at different times, can be used to observe the recovery rate and imbibition rate of each pore in the core column, is convenient for analyzing the factors and laws affecting imbibition, is convenient for conducting a phased analysis, understands the contribution of pores of each pore radius to the imbibition rate at any imbibition rate, and fully understands the participation of pores of each pore radius in the entire imbibition process. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0023] Figure 1A functional module diagram of a visualization experimental device for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application.

[0024] Figure 2 A schematic diagram of a pore radius distribution frequency curve generated by a visualization experimental device for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application.

[0025] Figure 3 A schematic diagram of an imbibition rate curve generated by a visualization experimental device for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application.

[0026] Figure 4 A schematic diagram of a single pore model constructed by a visualization experimental device for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application.

[0027] Figure 5 A visualization experimental device for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application simulates a pore with a pore radius of 0.5 nm, and is a schematic diagram showing the change in area of ​​crude oil displaced by water from 0h to 1h.

[0028] Figure 6 A schematic diagram of a critical oil yield curve generated by a visualization experimental device for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application.

[0029] Figure 7 A flowchart of a visualization experimental method for the microscopic imbibition rate of a shale reservoir provided in one embodiment of the present application.

[0030] Figure 8 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0031] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0032] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0033] In an exemplary embodiment, Figure 1As shown, a visualization experimental device for the microscopic imbibition rate of a shale reservoir is provided, comprising: a core pore radius scanning module, an imbibition physical simulation experiment module, a core simulation imaging module and a data solution and process control module; the data solution and process control module is used to: control the core pore radius scanning module to scan the pore radius of each pore of a target core column to obtain a pore scanning result of the target core column, and solve the pore radius distribution frequency curve of the target core column according to the pore scanning result of the target core column; the pore radius distribution frequency curve is used to characterize the occurrence of pores with different pore radii in the target core column. frequency; control the imbibition physical simulation experiment module, simulate the target core column, conduct static imbibition physical simulation experiment on the simulated target core column, and solve the imbibition rate curve according to the experimental results of the static imbibition physical simulation experiment; the imbibition rate curve is used to characterize the relationship between the oil production rate per unit area and time; control the core simulation imaging module, generate the core column finite element model according to the pore radius distribution frequency curve, and conduct static imbibition numerical simulation experiment on the core column finite element model based on the imbibition rate curve, and visualize the imbibition process in each pore of the core column finite element model.

[0034] certainly, Figure 1 The device architecture shown is only exemplary and can be omitted according to actual needs when implementing different functions. Figure 1 One or at least two components of the device shown.

[0035] Specifically in this embodiment, in the core pore radius scanning module, the target core column is scanned by a CT scanner to obtain pore radius data within the macropore range of 1000nm~100000nm, and a high-speed mercury injection test is performed on the target core column to obtain pore radius data within the range of mesopores and macropores from 2nm to 50000nm; the pore scanning result of the target core column includes pore radius data within the macropore range of 1000nm~100000nm and pore radius data within the range of mesopores and macropores from 2nm~50000nm; the data solution and process control module solves the pore radius distribution frequency curve of the target core column according to the pore radius data within the range of mesopores and macropores from 2nm~50000nm and the pore radius data within the range of macropores from 50000nm~100000nm, as shown in FIG. Figure 2 As shown in the figure. In the macropore range of 1000nm to 100000nm, CT scanning using the X-ray method has better accuracy, while the high-speed mercury injection method has better accuracy in the mesopore and macropore range of 2nm to 50000nm. Therefore, the combination of CT scanner and high-speed mercury injection scanner can better describe the proportion of each pore radius in the core column.

[0036] In an exemplary embodiment, the imbibition physical simulation experiment module conducts a static imbibition physical simulation experiment on a treated target core column using a volumetric method based on an imbibition bottle; the treatment of the target core column includes introducing simulated oil into the pores of the target core column after vacuum drying and simulating reservoir temperature aging for a preset time to simulate the core column in a real reservoir environment.

[0037] Due to the small pore volume of ultra-low permeability cores, it was found during the indoor experimental research that the volume expansion of the oil phase and the adsorption of oil on the surface of the core had a greater impact on the experimental results. To avoid interference with the experimental results, the experimental core was soaked in simulated oil at reservoir temperature for 24 hours before the experiment began. Specifically in this embodiment, before conducting a static imbibition physical simulation experiment, a core column sample was prepared by the following process, and the process for preparing the core column sample includes the following steps: 1) Use a vernier caliper to measure the length and diameter of the core.

[0038] 2) Dry the core under vacuum at 110°C for 48 h.

[0039] 3) Weigh the mass of the dried core using an analytical balance.

[0040] 4) Place the core into a vacuum pressurized saturation device and evacuate the core for 4 hours.

[0041] 5) Turn off the vacuum pump and use the pressure difference to introduce the simulated oil into the vacuum chamber containing the core so that the simulated oil submerges the core.

[0042] 6) Use a hand pump to pressurize the simulated oil in the vacuum chamber to 18 MPa and maintain the pressurized state for 24 hours.

[0043] 7) Open the pressure relief valve, take out the core, transfer the core to a container and immerse it in simulated oil. Age it at simulated reservoir temperature (60°C) for 14 days, wipe off the simulated oil on the surface of the core, weigh the core using an analytical balance, calculate the mass of the simulated oil saturated into the core, and set it aside.

[0044] Before the core column sample is taken out from the simulated oil and transferred to the imbibition liquid phase, the oil phase adsorbed on the core surface is wiped clean and quickly transferred into the imbibition liquid at the reservoir temperature.

[0045] The data solution and process control module reads the oil output of the static imbibition physical simulation experiment at preset intervals to calculate the imbibition rate at the current moment, and plots the imbibition rate at each moment of the static imbibition physical simulation experiment to obtain the imbibition rate curve of the target core column, such as Figure 3 shown.

[0046] In this embodiment, a simple and easy-to-operate volume method is used to carry out the imbibition experiment, and the simulated reservoir temperature is 60°C. When the static imbibition physical simulation experiment is carried out in the imbibition bottle, the oil output is read every 1 hour, and the imbibition rate is plotted to obtain the imbibition rate under the physical simulation experiment of the core. The physical simulation experiment cannot realize the visualization of the imbibition rate. Specifically, the produced oil volume and its corresponding time are recorded, and the recovery rate-time curve is plotted; the oil production-time curve is differentially calculated and processed, and the oil production rate-time curve and the unit area oil production rate-time curve are plotted; according to the recovery rate-time curve, the oil production rate-time curve and the unit area oil production rate-time curve, the nanofluid imbibition and oil drainage performance can be evaluated.

[0047] In this embodiment, the core simulation imaging module uses Comsol simulation software. The core simulation imaging module is used to: determine the proportion of pores with different pore radii in the target core column according to the pore radius distribution frequency curve; generate a core column finite element model according to the proportion of pores with different pore radii in the target core column; the core column finite element model and the target core column have exactly the same pore radius distribution frequency curve; fill the pores of the core column finite element model with simulated oil and set water at both ends of the pores; perform a static imbibition numerical simulation experiment on the core column finite element model according to the imbibition rate curve, and visualize the imbibition process in each pore of the core column finite element model.

[0048] Specifically, Comsol simulation software is used to implement the following steps: 1) Divide the x-axis of the pore radius distribution frequency curve into several equal parts, for example, in the range of 0 to 100,000 nm, divide it into 100 equal parts.

[0049] 2) Normalize the corresponding areas of these parts so that we can know the proportion of each interval in the future, which is convenient for the subsequent modeling of the Comsol finite element model.

[0050] 3) Select a representative radius in the interval for modeling, for example: 0-1000nm select 500nm as the representative radius of the interval.

[0051] 4) Percentage the representative radius of each interval, for example: 0-1000nm accounts for 5%, then create 5 pores with a radius of 5nm in the finite element model of Comsol software.

[0052] 5) Repeat step 4) to create a total of 100 parallel pores with representative radii in their respective intervals.

[0053] 6) When performing numerical simulations, fill the pores with oil and then set water on both sides of the pores. On this basis, no pressure is added and the static imbibition rate is observed.

[0054] As an example, when Comsol software is modeling, each pore can be modeled, such as Figure 4 As shown in the figure, by filling the drawn pores with oil, setting the boundary pressure to 0 and the wetting angle to 30°, a visual simulation of the static spontaneous imbibition experiment is performed. The animation of Comsol software can observe the static spontaneous imbibition recovery degree in the pore model at various time points.

[0055] By simulating and observing the simulated images at different times, the imbibition rate and recovery degree of the core column at each time point can be intuitively observed. The simulation process is as follows: Figure 5 As shown, it characterizes the change in the area of ​​crude oil displaced by water from 0h to 1h in a pore with a pore radius of 0.5nm; the imbibition at 0.2h is also shown in this interval. The recovery degree is the area of ​​crude oil (dark gray area on the right) displaced by water (black area on the left) in the figure.

[0056] In another exemplary embodiment, in the discussion of the factors affecting gravity and capillary force, gravity is the attraction of the earth on a substance, and its direction always points downward toward the center of the earth. Oil and water in porous media are both acted downward by gravity, but the density of water is greater than that of oil. The difference in gravity between oil and water will generate a buoyancy, and the direction of the buoyancy is vertically upward. When the direction of the buoyancy is consistent with the direction of oil infiltration and drainage, the gravity effect is conducive to infiltration; if the opposite is true, it is not conducive to infiltration.

[0057] The Bond inverse is often used to describe the relative contributions of capillary force and gravity in the imbibition process. It is also applicable to the case of a single microtube in the Carman-Kozeny equation. NB -1 When ≥5, the dominant force in the oil-imbibition process is capillary force; NB -1 ≤1, the dominant force in the oil drainage process is gravity. NB -1 When <5, gravity and capillary force act together on oil absorption and drainage.

[0058] .

[0059] in, NB -1 is the reciprocal of the Bond number, dimensionless, C is a constant. For the single-pore model C =0.4, ϕ is the porosity, dimensionless, k is the permeability of rock, μm 2 , r is the pore radius, Δpis the density difference between oil and water, g is the acceleration due to gravity, h is the height difference between the oil and water layers.

[0060] For 1< NB -1 When the value is less than 5, we cannot use the overall concept to view the contribution of gravity and capillary force to the oil infiltration and drainage of the core column. At this time, we need to refine the state of pores with a certain radius affected by gravity and capillary force. For this reason, this application uses Comsol to simulate different pore radii at a specific wetting angle (different wetting angles characterize the wettability of the porous rock surface, that is, whether the rock surface is hydrophilic or oleophilic. It is used to characterize the degree of hydrophilicity of the rock wall) and the relationship equation between capillary force and wetting angle is as follows: In the formula, P c is the capillary pressure, Pa; σ is the interfacial tension between the non-wetting phase and the wetting phase, mN / m; r is the pore radius, μm; θ is the three-phase contact angle. The capillary pressure is the pressure difference between the wetting phase and the non-wetting phase on both sides of the fluid curved liquid surface. The capillary pressure always points to the concave side of the curved liquid surface. The capillary pressure is usually described by the Washburn formula. The capillary pressure is proportional to the interfacial tension and the cosine of the contact angle, and inversely proportional to the pore radius.

[0061] In this embodiment, the core simulation imaging module is also used to: model pores with different pore radii respectively to obtain several pore models corresponding to different pore radii; place several pore models vertically, and draw a critical oil yield curve between pore radius and wetting angle through static imbibition numerical simulation experiments of different pore radii under a predetermined wetting angle and static imbibition numerical simulation experiments of different wetting angles under a predetermined pore radius; the horizontal and vertical coordinates of the critical oil yield curve are the pore radius and the wetting angle, respectively, and any point on the critical oil yield curve represents that under the condition of this pore radius and this wetting angle, the simulated oil in the pore model begins to displace in the pore model due to the action of gravity.

[0062] In the Comsol simulation software, a pore model of a certain pore radius is placed vertically. By simulating different pore radii with a specific wetting angle and different wetting angles with a specific pore radius, the critical oil yield curve between pore radius and wetting angle is drawn, such as Figure 6 As shown, a point with a pore radius of 30 microns is used as an example for the experiment. The wetting angle is continuously changed during the simulation process until the crude oil is displaced in the pore due to gravity at a specific wetting angle. Because the capillary force at this time is a resistance, it prevents the crude oil from displacing in the pore due to gravity. Based on this, it is determined that the crude oil in the pore is affected by gravity as much as by capillary force.

[0063] Based on the same inventive concept, the present application also provides a method for applying the shale reservoir microscopic imbibition rate visualization experimental device provided in the above embodiment. The solution to the problem provided by the method is similar to the solution described in the above device. In an exemplary embodiment, Figure 7 As shown, a visualization experimental method for the microscopic imbibition rate of a shale reservoir is provided, comprising the following steps: S1. Scan the pore radius of each pore of the acquired target core column to obtain the pore scanning result of the target core column. Specifically in this embodiment, the pore scanning result of the target core column includes: pore radius data within the macropore range of 1000nm~100000nm and pore radius data within the range of mesopores and macropores of 2nm~50000nm; step S1 specifically includes the following steps: S11. Scan the target core column with a CT scanner to obtain pore radius data within the macropore range of 1000 nm to 100000 nm.

[0064] S12. Perform a mercury injection experiment on the target core column using a high-speed mercury injection scanner to obtain pore radius data within the range of 2nm to 50000nm for mesopores and macropores.

[0065] S2. According to the pore scanning result of the target core column, the pore radius distribution frequency curve of the target core column is calculated; the pore radius distribution frequency curve is used to characterize the occurrence frequency of pores with different pore radii in the target core column.

[0066] S3. Simulate the target core column and conduct a static imbibition physical simulation experiment on the simulated target core column.

[0067] S4. According to the experimental results of the static imbibition physical simulation experiment, the imbibition rate curve is calculated; the imbibition rate curve is used to characterize the relationship between the oil production rate per unit area and time.

[0068] S5. Generate a core column finite element model according to the pore radius distribution frequency curve. Specifically in this embodiment, step S5 includes the following steps: S51. Determine the proportion of pores with different pore radii in the target core column according to the pore radius distribution frequency curve.

[0069] S52. Generate a core column finite element model according to the proportion of pores with different pore radii in the target core column; the core column finite element model and the target core column have exactly the same pore radius distribution frequency curve.

[0070] S6. Based on the imbibition rate curve, a static imbibition numerical simulation experiment is performed on the core column finite element model to visualize the imbibition process in each pore of the core column finite element model. Specifically in this embodiment, simulated oil is filled in the pores of the core column finite element model and water is set at both ends of the pores. According to the imbibition rate curve, a static imbibition numerical simulation experiment is performed on the core column finite element model to visualize the imbibition process in each pore of the core column finite element model.

[0071] The above-mentioned scheme provided in the present application can be used to visualize the imbibition process of the imbibition physical simulation experiment at different times, can be used to observe the recovery rate and imbibition rate of each pore in the core column, is convenient for analyzing the factors and laws affecting imbibition, is convenient for conducting a phased analysis, understands the contribution of pores of each pore radius to the imbibition rate at any imbibition rate, and fully understands the participation of pores of each pore radius in the entire imbibition process.

[0072] In another exemplary embodiment, a computer device is also provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 8 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, the visualization experimental method of the microscopic imbibition rate of the shale reservoir provided in the above embodiment can be implemented.

[0073] Those skilled in the art will understand that Figure 8 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0074] In an exemplary embodiment, a computer device is further provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0075] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0076] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0077] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0078] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0079] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

[0080] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0081] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A visualization experimental device for the microscopic imbibition rate of shale reservoirs, characterized in that: include: A core pore radius scanning module, an imbibition physical simulation experiment module, a core simulation imaging module and a data solution and process control module; the data solution and process control module is used to: control the core pore radius scanning module to scan the pore radius of each pore of the target core column to obtain the pore scanning result of the target core column, and solve the pore radius distribution frequency curve of the target core column according to the pore scanning result of the target core column; the pore radius distribution frequency curve is used to characterize the occurrence frequency of pores with different pore radii in the target core column; control the imbibition physical simulation experiment module , simulate the target core column, conduct a static imbibition physical simulation experiment on the simulated target core column, and calculate the imbibition rate curve according to the experimental results of the static imbibition physical simulation experiment; the imbibition rate curve is used to characterize the relationship between the oil production rate per unit area and time; control the core simulation imaging module, generate a core column finite element model according to the pore radius distribution frequency curve, and conduct a static imbibition numerical simulation experiment on the core column finite element model based on the imbibition rate curve, and visualize the imbibition process in each pore of the core column finite element model.

2. The visualization experimental device for the microscopic imbibition rate of shale reservoir according to claim 1, characterized in that: In the core pore radius scanning module, the target core column is scanned by a CT scanner to obtain pore radius data within the macropore range of 1000nm~100000nm, and a high-speed mercury injection test is performed on the target core column to obtain pore radius data within the range of mesopores and macropores of 2nm~50000nm; the pore scanning result of the target core column includes pore radius data within the macropore range of 1000nm~100000nm and pore radius data within the range of mesopores and macropores of 2nm~50000nm; the data solution and process control module solves and obtains the pore radius distribution frequency curve of the target core column according to the pore radius data within the range of mesopores and macropores of 2nm~50000nm and the pore radius data within the range of macropores of 50000nm~100000nm.

3. The visualization experimental device for the microscopic imbibition rate of shale reservoir according to claim 1, characterized in that: The imbibition physical simulation experiment module uses the volume method based on the imbibition bottle to perform a static imbibition physical simulation experiment on the treated target core column; the treatment of the target core column includes introducing simulated oil into the pores of the target core column after vacuum drying and simulating the reservoir temperature aging for a preset time to simulate the core column in a real reservoir environment.

4. The visualization experimental device for the microscopic imbibition rate of shale reservoir according to claim 1, characterized in that: The data solution and process control module reads the oil output of the static imbibition physical simulation experiment at preset intervals to calculate the imbibition rate at the current moment, and plots the imbibition rate at each moment of the static imbibition physical simulation experiment to obtain the imbibition rate curve of the target core column.

5. The visualization experimental device for the microscopic imbibition rate of shale reservoir according to claim 1, characterized in that: The core simulation imaging module is used to: determine the proportion of pores with different pore radii in the target core column according to the pore radius distribution frequency curve; generate a core column finite element model according to the proportion of pores with different pore radii in the target core column; the core column finite element model and the target core column have exactly the same pore radius distribution frequency curve; fill the pores of the core column finite element model with simulated oil and set water at both ends of the pores; perform a static imbibition numerical simulation experiment on the core column finite element model according to the imbibition rate curve, and visually display the imbibition process in each pore of the core column finite element model.

6. The visualization experimental device for the microscopic imbibition rate of shale reservoir according to claim 1, characterized in that: The core simulation imaging module is also used to: model pores with different pore radii respectively to obtain several pore models corresponding to different pore radii; place several of the pore models vertically, and draw a critical oil yield curve between pore radius and wetting angle through static imbibition numerical simulation experiments at different pore radii under a predetermined wetting angle and static imbibition numerical simulation experiments at different wetting angles under a predetermined pore radius; the abscissa and ordinate of the critical oil yield curve are the pore radius and the wetting angle, respectively, and any point on the critical oil yield curve represents that under the condition of the pore radius and the wetting angle, the simulated oil in the pore model begins to displace in the pore model due to the action of gravity.

7. A visualization experimental method for microscopic imbibition rate of shale reservoirs, characterized in that: include: Scanning the pore radius of each pore of the acquired target core column to obtain a pore scanning result of the target core column; According to the pore scanning result of the target core column, a pore radius distribution frequency curve of the target core column is obtained by solving the pore radius distribution frequency curve; the pore radius distribution frequency curve is used to characterize the occurrence frequency of pores with different pore radii in the target core column; Perform simulation processing on the target core column, and conduct static imbibition physical simulation experiment on the simulated target core column; According to the experimental results of the static imbibition physical simulation experiment, an imbibition rate curve is calculated; the imbibition rate curve is used to characterize the relationship between the oil production rate per unit area and time; Generating a core column finite element model according to the pore radius distribution frequency curve; Based on the imbibition rate curve, a static imbibition numerical simulation experiment is performed on the core column finite element model to visualize the imbibition process in each pore of the core column finite element model.

8. The visualization experimental method of shale reservoir microscopic imbibition rate according to claim 7, characterized in that: The pore scanning results of the target core column include: pore radius data within the macropore range of 1000nm to 100000nm and pore radius data within the range of mesopores and macropores of 2nm to 50000nm; scanning the pore radius of each pore of the target core column obtained specifically includes: The target core column is scanned by a CT scanner to obtain pore radius data within the macropore range of 1000nm to 100000nm; A mercury injection experiment was carried out on the target core column using a high-speed mercury injection scanner to obtain pore radius data within the range of 2nm~50000nm for mesopores and macropores.

9. The visualization experimental method of shale reservoir microscopic imbibition rate according to claim 7, characterized in that: According to the pore radius distribution frequency curve, a core column finite element model is generated, which specifically includes: According to the pore radius distribution frequency curve, determining the proportion of pores with different pore radii in the target core column; A core column finite element model is generated according to the proportion of pores with different pore radii in the target core column; the core column finite element model and the target core column have exactly the same pore radius distribution frequency curve.

10. The visualization experimental method of shale reservoir microscopic imbibition rate according to claim 7, characterized in that: The pores of the core column finite element model are filled with simulated oil and water is set at both ends of the pores. A static imbibition numerical simulation experiment is performed on the core column finite element model according to the imbibition rate curve to visualize the imbibition process in each pore of the core column finite element model.

Images