CN117390997A - Digital experimental method and device for fracture-cavity oil reservoir physical model - Google Patents

Digital experimental method and device for fracture-cavity oil reservoir physical model Download PDF

Info

Publication number
CN117390997A
CN117390997A CN202311532241.6A CN202311532241A CN117390997A CN 117390997 A CN117390997 A CN 117390997A CN 202311532241 A CN202311532241 A CN 202311532241A CN 117390997 A CN117390997 A CN 117390997A
Authority
CN
China
Prior art keywords
model
simulation
digital model
digital
grid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202311532241.6A
Other languages
Chinese (zh)
Inventor
李俊超
畅新鸽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Shiyou University
Original Assignee
Xian Shiyou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN202311532241.6A priority Critical patent/CN117390997A/en
Publication of CN117390997A publication Critical patent/CN117390997A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application discloses a method and a device for digitized experiment of a fracture-cave oil reservoir physical model, wherein the method comprises the following steps: constructing a digital model of the physical model; gridding the digital model to generate an unstructured grid file; introducing the unstructured grid file into fluid mechanics simulation software, and setting simulation parameters to generate a simulation model; initializing a simulation model, and setting solving parameters to perform simulation calculation. By implementing the method, the problems that the physical simulation method in the prior art is poor in visual effect, long in experimental period and low in precision of multi-phase fluid flow and multi-phase fluid interface depiction are effectively solved. Furthermore, the method for digitally testing the fracture-cavity oil reservoir physical model can realize model visualization, can accelerate the experimental process and saves the experimental time.

Description

Digital experimental method and device for fracture-cavity oil reservoir physical model
Technical Field
The application relates to the technical field of oil reservoir numerical simulation, in particular to a method and a device for digitized experiment of a fracture-cavity oil reservoir physical model.
Background
In the fracture-cavity type oil reservoir space, the porous medium seepage flow and the large-space pipe flow are complex multi-fluid coupling flow. The flow rule of the multi-phase fluid of the fracture-cavity type carbonate reservoir is researched, and the method has important data guiding significance for the development of oil and gas fields. In the prior art, a physical simulation method and a simulation method are generally used.
The physical simulation method is to scale down the dynamic parameters and the static parameters of the real model according to the similar principle so as to measure the required data on one observation scale in a short time. However, in the physical simulation process, the problems of complex manufacturing of a physical model, poor visual effect, long experimental period and the like are faced, and the temperature and pressure conditions of the physical model are difficult to meet. In addition, a physical model is generally only capable of carrying out research on a certain type of fracture-cavity unit under a specific displacement condition, and knowledge and conclusion of universality are difficult to obtain. If the mining rules or the flow rules of different types of fracture-cavity units are to be researched, a large amount of physical simulation researches are required to be carried out, and a large amount of manpower and material resources are consumed.
The simulation method is an oil reservoir numerical simulation method based on an equivalent medium model to simulate the displacement process of the fracture-cavity oil reservoir physical model. The equivalent medium model is visual, simple in modeling and high in simulation speed, and is well applied to the processing of the oil reservoir scale model. However, the method is based on the Darcy percolation equation, has lower precision for describing multiphase fluid flow and multiphase fluid interfaces, and cannot meet the precision requirement of a digital experiment.
Disclosure of Invention
The embodiment of the application solves the problems of complex model manufacture, poor visual effect, longer test period and lower precision of a simulation method of a physical simulation method in the prior art by providing the digital experiment method for the physical model of the hole oil reservoir. An experimental method for digitizing a physical force model is realized, which can solve the problems.
In a first aspect, an embodiment of the present application provides a method for digitized experiment of a physical model of a hole reservoir, including: constructing a digital model of the physical model; gridding the digital model to generate an unstructured grid file; introducing the unstructured grid file into fluid mechanics simulation software, and setting simulation parameters to generate a simulation model; initializing the simulation model, and setting solving parameters to perform simulation calculation.
With reference to the first aspect, in a first possible implementation manner, the constructing a digital model of a physical model includes: inserting the view of the physical model into three-dimensional software for copying modeling, and adjusting the model to be consistent with the size of the physical model; and adjusting the width of a crack channel of the model according to the effective volume of the physical model until the volume of the model is close to the physical model to obtain the digital model.
With reference to the first aspect, in a second possible implementation manner, the meshing the digital model to generate an unstructured mesh file includes: determining the gridding grid type and the target grid number of the digital model; calculating the grid resolution of the digital model; and generating and outputting an unstructured grid file of the digital model according to the grid type, the target grid number and the grid resolution.
With reference to the second possible implementation manner of the first aspect, in a third possible implementation manner, the calculation formula for calculating the grid resolution of the digital model is as follows:the method comprises the steps of carrying out a first treatment on the surface of the In the method, in the process of the invention,for the grid resolution, V is the volume of the digital model, w is the thickness of the digital model,and the target grid number.
With reference to the second possible implementation manner of the first aspect, in a fourth possible implementation manner, the meshing the digital model to generate an unstructured grid file further includes: checking the dividing quality of the unstructured grid file; and if the negative volume appears in the unstructured grid file, the digital model is gridded again until the generated unstructured grid file has no negative volume.
With reference to the first aspect, in a fifth possible implementation manner, the setting simulation parameters includes: setting gravity acceleration for the digital model, and selecting a model type; setting the number of fluid phases in the digital model, and setting the density and viscosity of each phase of fluid; and setting an interfacial tension coefficient of the fluid for the digital model, and selecting a viscous flow model type.
With reference to the first aspect, in a sixth possible implementation manner, the solution parameter includes a pressure velocity coupling scheme, a gradient space discrete type, a pressure space discrete type, a momentum space discrete type, a volume fraction space discrete type, a turbulence kinetic energy space discrete type, a specific dissipation ratio space discrete type, an inlet boundary condition, an outlet boundary condition, a time step, a number of time steps, and a number of iterations of the time step.
In a second aspect, an embodiment of the present application provides a hole reservoir physical model digitizing experimental apparatus, including: the construction module is used for constructing a digital model of the physical model; the generation module is used for meshing the digital model to generate an unstructured grid file; the simulation module is used for guiding the unstructured grid file into fluid mechanics simulation software and setting simulation parameters to generate a simulation model; and the calculation module is used for initializing the simulation model and setting solving parameters to perform simulation calculation.
With reference to the second aspect, in a first possible implementation manner, the constructing a digital model of a physical model includes: inserting the view of the physical model into three-dimensional software for copying modeling, and adjusting the model to be consistent with the size of the physical model; and adjusting the width of a crack channel of the model according to the effective volume of the physical model until the volume of the model is close to the physical model to obtain the digital model.
With reference to the second aspect, in a second possible implementation manner, the meshing the digital model to generate an unstructured grid file includes: determining the gridding grid type and the target grid number of the digital model; calculating the grid resolution of the digital model; and generating and outputting an unstructured grid file of the digital model according to the grid type, the target grid number and the grid resolution.
Second aspect combined with the second aspectIn a third possible implementation manner, the calculation formula for calculating the grid resolution of the digital model is as follows:the method comprises the steps of carrying out a first treatment on the surface of the In the method, in the process of the invention,for the grid resolution, V is the volume of the digital model, w is the thickness of the digital model,and the target grid number.
With reference to the second possible implementation manner of the second aspect, in a fourth possible implementation manner, the meshing the digital model to generate an unstructured grid file further includes: checking the dividing quality of the unstructured grid file; and if the negative volume appears in the unstructured grid file, the digital model is gridded again until the generated unstructured grid file has no negative volume.
With reference to the second aspect, in a fifth possible implementation manner, the setting simulation parameters includes: setting gravity acceleration for the digital model, and selecting a model type; setting the number of fluid phases in the digital model, and setting the density and viscosity of each phase of fluid; and setting an interfacial tension coefficient of the fluid for the digital model, and selecting a viscous flow model type.
With reference to the second aspect, in a sixth possible implementation manner, the solution parameter includes a pressure velocity coupling scheme, a gradient space discrete type, a pressure space discrete type, a momentum space discrete type, a volume fraction space discrete type, a turbulence kinetic energy space discrete type, a specific dissipation ratio space discrete type, an inlet boundary condition, an outlet boundary condition, a time step, a number of time steps, and a number of iterations of the time step.
In a third aspect, embodiments of the present application provide an apparatus, including: a processor; a memory for storing processor-executable instructions; the processor, when executing the executable instructions, implements a method as described in the first aspect or any one of the possible implementations of the first aspect.
In a fourth aspect, embodiments of the present application provide a non-transitory computer readable storage medium comprising instructions for storing a computer program or instructions which, when executed, cause a method as described in the first aspect or any one of the possible implementations of the first aspect to be implemented.
One or more technical solutions provided in the embodiments of the present application at least have the following technical effects or advantages:
according to the embodiment of the application, the digital model for constructing the physical model is adopted, and the fluid mechanics simulation software is utilized to construct the simulation model, so that the problems that in the prior art, the physical simulation method is poor in visualization effect and long in experimental period, and the simulation method is low in precision for describing multiphase fluid flow and multiphase fluid interfaces are effectively solved. The hole oil reservoir physical model digital experiment method which can realize visualization and has the experiment precision similar to that of a physical simulation method can be further realized, the experiment process can be quickened, the experiment time is saved, and the consumption of manpower and material resources in the experiment process is reduced.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings used in the embodiments of the present application or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for digitized experiments on a physical model of a hole reservoir according to an embodiment of the present application;
FIG. 2 is a flowchart of gridding a digital model to generate an unstructured grid file according to an embodiment of the present application;
FIG. 3 is a block diagram of a hole reservoir physical model digital experimental device provided in an embodiment of the present application;
FIG. 4a is a schematic diagram of the ratio of water to oil before water flooding for a simulation model according to an embodiment of the present disclosure;
FIG. 4b is a schematic diagram of the ratio of water to oil after water flooding of the simulation model according to the embodiment of the present application;
fig. 5 is an injection-production curve of a physical experiment provided in an embodiment of the present application and a simulated injection-production curve of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. It will be apparent that the described embodiments are some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Some of the techniques involved in the embodiments of the present application are described below to aid understanding, and they should be considered as merely exemplary. Accordingly, one of ordinary skill in the art will recognize that various changes and modifications of the embodiments described herein can be made without departing from the scope and spirit of the present application. Also, for the sake of clarity and conciseness, descriptions of well-known functions and constructions are omitted in the following description.
Fig. 1 is a flowchart of a method for digitized experiment of a physical model of a hole reservoir according to an embodiment of the present application, including steps 101 to 104. Wherein fig. 1 is only one execution sequence shown in the embodiment of the present application, and does not represent the only execution sequence of the hole reservoir physical model digital experiment method, and the steps shown in fig. 1 may be executed in parallel or in reverse in case that the final result can be achieved.
Step 101: a digital model of the physical model is constructed. Specifically, a view of the physical model is inserted into three-dimensional software for copy modeling, and the model is adjusted to conform to the size of the physical model. In the present embodiment, three-dimensional software is exemplified by CAD software. After the view of the physical model is inserted into CAD software, the contrast of the picture is firstly adjusted to make the lines obvious. And copying by using spline fitting curves, merging the lines, and deleting the two middle surfaces. Finally, the copy and the combined model are scaled to be equal to the physical model in size.
And adjusting the width of the crack channel of the model according to the effective volume of the physical model until the volume of the model is close to the physical model to obtain a digital model. Specifically, the relative error due to the slit channel width is large. Therefore, the volume of the model is calculated first, and is compared with the effective volume of the physical model, so that the width of a crack channel in the model is adjusted until the volume of the model is close to the volume of the real physical model. And exporting the adjusted model to obtain a digital model.
Step 102: the digital model is gridded to generate an unstructured grid file. The specific implementation method of step 102 is shown in fig. 2, and includes steps 201 to 203, which are specifically as follows.
Step 201: the gridded grid type and the target grid number of the digital model are determined. In the embodiment of the application, since the boundary of the visualized flat panel fracture-cavity model is very complex, the triangle mesh type is adopted to carry out mesh division on the digital model, and the target mesh number for dividing the digital model is determined. It should be understood by those skilled in the art that the triangular mesh is only one example shown in the present application, and is not intended to limit the scope of the present application, and those skilled in the art may make modifications or changes according to practical experimental needs, and still fall within the scope of the present application.
Step 202: the grid resolution of the digital model is calculated. In the embodiment of the present application, a calculation formula for calculating the grid resolution of the digital model is as follows:
. In the method, in the process of the invention,for the resolution of the grid, V is the volume of the digital model, w is the thickness of the digital model,is the target grid number.
Step 203: and generating and outputting an unstructured grid file of the digital model according to the grid type, the target grid number and the grid resolution. Specifically, an unstructured grid file of the digital model is generated according to the grid type, the target grid number and the grid resolution obtained through calculation in the steps and is output.
In the embodiment of the application, when the unstructured grid file is solved by using the fluid mechanics simulation software, the calculation result is set to be of a double-precision type. Before solving the unstructured grid file, the division quality of the unstructured grid file is checked. Judging whether negative volumes appear in the unstructured grid file, and if so, re-gridding the digital model until no negative volumes exist in the generated unstructured grid file. The occurrence of solving failure caused by negative volume can be avoided through the steps.
Step 103: the unstructured grid file is led into fluid mechanics simulation software, and simulation parameters are set to generate a simulation model. In the embodiment of the present application, the fluid mechanics simulation software used is Fluent. Unstructured grid files are imported in hydrodynamic simulation software. The gravity acceleration is set for the digital model and the model type is selected. Specifically, the gravitational acceleration is exemplarily set to-9.81 m/s 2 The direction is the Y direction. Meanwhile, a multiphase flow model-Euler model is selected as a model type of a digital model in fluid mechanics simulation software. And the interface type at this time is set to be Sharp type, so that a more accurate phase interface can be obtained. Further, in the embodiment of the present application, the display format of the volume count of the digital model is set to a discrete format.
The number of fluid phases in the digital model is set, and the density and viscosity of each phase of fluid are set. In the embodiment of the application, the fluid phase number is set according to the kind number of the actual fluid used in the displacement process. For example, in one displacement process, three fluids, oil (oil phase), water (liquid phase), and air (gas phase), respectively, are present, and the number of fluid phases is set to three, and the main phase therein is set to be the oil phase. The density and viscosity of each phase of fluid is set according to the experimental requirements, and is specific data during the experiment and will not be described in detail herein.
The interfacial tension coefficient of the fluid is set for the digital model and the viscous flow model type is selected. In the embodiment of the application, the interfacial tension coefficient of the fluid is set, and the oil-water interfacial tension is set to be 0.072N/m and the gas-liquid interfacial tension is set to be 0.02N/m by way of example. And selectThe model is a viscous flow model type.
Step 104: initializing a simulation model, and setting solving parameters to perform simulation calculation. In the embodiment of the application, the height proportion of the different-phase fluid in the view of the physical model is calculated respectively, and then the space height occupied by the different-phase fluid in the simulation model is calculated according to the total height of the simulation model. And calculating Y coordinates of phase interfaces of the different phase fluids according to the spatial distribution condition of the different phase fluids and the space height occupied by the different phase fluids in the simulation model. In the embodiment of the application, the spatial distribution condition of different phase fluids is that water is at the bottom layer, oil is at the middle layer and air is at the upper layer. The Y coordinate at the top of the space where the water is located is equal to the Y coordinate at the bottom of the space where the oil is located, specifically equal to the Y coordinate at the bottom of the space where the water is located plus the space height of the space where the water is located. Similarly, the Y coordinate of the top of the space where the oil is located is equal to the Y coordinate of the bottom of the space where the air is located, and specifically equal to the Y coordinate of the bottom of the space where the water is located plus the space height of the space where the oil is located. The Y coordinate of the top of the space where the air is located is the Y coordinate of the bottom of the space where the water is located plus the space height of the space where the oil is located plus the space height of the space where the air is located.
Global initialization is set, the initial phase of the simulation model is set as the water phase, i.e. the volume fraction of the initial water phase is set as 1. And then carrying out local initialization setting, setting the volume fraction of the region where the oil is located as 1 according to the spatial distribution condition of the different phase fluids, and setting the volume fraction of the region where the water is located as 1.
The solving parameters comprise a pressure velocity coupling scheme, a gradient space discrete type, a pressure space discrete type, a momentum space discrete type, a volume fraction space discrete type, a turbulence kinetic energy space discrete type, a specific dissipation ratio space discrete type, an inlet boundary condition, an outlet boundary condition, a time step number and an iteration number of the time step.
In the embodiment of the present application, the pressure-velocity coupling scheme is set to be the Coupled scheme. The gradient spatial discrete type is set to Least Squarries Cell Based. The pressure space discrete type is set to prest. The momentum space discrete type is set as Second Order Upwind. The volume fraction spatial discrete type is set to compression. The turbulent kinetic energy space discrete type is set as Second Order Upwind. The specific dissipation ratio spatial discrete type is set to Second Order Upwind. And inlet boundary conditions are set according to experimental injection speed and outlet boundary conditions are set according to experimental outlet conditions. Setting the time step and the time step number to enable the Kurthan numberThe following formula is obtained:
. In the method, in the process of the invention,for the time step size of the time step,for the resolution of the grid it is advantageous,for the purpose of experimental injection rates,for the number of kulange's,in order to make the number of time steps,is the total duration of the experiment. And the number of iterations per time step is illustratively set to 20. It should be noted that the above setting of the solution parameters is only an example of the embodiments of the present application, and is not intended to limit the scope of the present application, and those skilled in the art may make various modifications and changes to the solution parameters according to the actual usage scenario, which still falls within the scope of the present application.
And carrying out numerical simulation calculation according to the set solving parameters to obtain a pressure field, a speed field and a volume fraction distribution field of the simulation model, and a change curve of flow and pressure at an outlet and the volume fraction of the oil, gas and water phases along with time. The data obtained by the solution can be used for guiding oil and gas exploitation work. Fig. 4a to fig. 4b are schematic diagrams showing the ratio change of water to oil of the simulation model provided in the present application before water flooding and after water flooding, wherein the black area is oil, and the white area is water. The flow rule of the multiphase fluid studied by the application is known to have guiding significance on oil and gas exploitation. Fig. 5 shows an injection-production curve of a physical experiment provided in an embodiment of the present application and a simulated injection-production curve of the present application. It can be seen from fig. 5 that the accuracy of the method used in the present application is close to that of the physical simulation method.
Although the present application provides method operational steps as described in the examples or flowcharts, more or fewer operational steps may be included based on conventional or non-inventive labor. The order of steps recited in the present embodiment is only one way of performing the steps in a plurality of steps, and does not represent a unique order of execution. When implemented by an actual device or client product, the method of the present embodiment or the accompanying drawings may be performed sequentially or in parallel (e.g., in a parallel processor or a multithreaded environment).
As shown in fig. 3, the embodiment of the present application further provides a hole reservoir physical model digitizing experimental apparatus 300. The device comprises: the construction module 301, the generation module 302, the simulation module 303 and the calculation module 304 are specifically as follows.
The construction module 301 is used for constructing a digital model of the physical model. The building module 301 is specifically configured to insert a view of the physical model into three-dimensional software for copy modeling, and adjust the model to conform to the physical model. In the present embodiment, three-dimensional software is exemplified by CAD software. After the view of the physical model is inserted into CAD software, the contrast of the picture is firstly adjusted to make the lines obvious. And copying by using spline fitting curves, merging the lines, and deleting the two middle surfaces. Finally, the models after copying and merging are scaled down to be equal to the physical model in size.
And adjusting the width of the crack channel of the model according to the effective volume of the physical model until the volume of the model is close to the physical model to obtain a digital model. Specifically, the relative error due to the slit channel width is large. Therefore, the volume of the model is calculated first, and is compared with the effective volume of the physical model, so that the width of a crack channel in the model is adjusted until the volume of the model is close to the volume of the real physical model. And exporting the adjusted model to obtain a digital model.
The generation module 302 is configured to grid the digital model to generate an unstructured grid file. The generating module 302 is specifically configured to determine a gridded grid type and a target grid number of the digital model. In the embodiment of the application, since the boundary of the visualized flat panel fracture-cavity model is very complex, the triangle mesh type is adopted to carry out mesh division on the digital model, and the target mesh number for dividing the digital model is determined.
The grid resolution of the digital model is calculated. In the embodiment of the present application, a calculation formula for calculating the grid resolution of the digital model is as follows:
. In the method, in the process of the invention,for the resolution of the grid, V is the volume of the digital model, w is the thickness of the digital model,for a target gridA number.
And generating and outputting an unstructured grid file of the digital model according to the grid type, the target grid number and the grid resolution. Specifically, an unstructured grid file of the digital model is generated according to the grid type, the target grid number and the grid resolution obtained through calculation in the steps and is output.
In the embodiment of the application, when the unstructured grid file is solved by using the fluid mechanics simulation software, the calculation result is set to be of a double-precision type. Before solving the unstructured grid file, the division quality of the unstructured grid file is checked. Judging whether negative volumes appear in the unstructured grid file, and if so, re-gridding the digital model until no negative volumes exist in the generated unstructured grid file. The occurrence of solving failure caused by negative volume can be avoided through the steps.
The simulation module 303 is used for importing the unstructured grid file into the fluid mechanics simulation software and setting simulation parameters to generate a simulation model. The simulation module 303 is specifically configured to use Fluent simulation software in the embodiment of the present application. Unstructured grid files are imported in hydrodynamic simulation software. The gravity acceleration is set for the digital model and the model type is selected. Specifically, the gravitational acceleration is exemplarily set to-9.81 m/s 2 The direction is the Y direction. Meanwhile, a multiphase flow model-Euler model is selected as a model type of a digital model in fluid mechanics simulation software. And the interface type at this time is set to be Sharp type, so that a more accurate phase interface can be obtained. Further, in the embodiment of the present application, the display format of the volume count of the digital model is set to a discrete format.
The number of fluid phases in the digital model is set, and the density and viscosity of each phase of fluid are set. In the embodiment of the application, the fluid phase number is set according to the kind number of the actual fluid used in the displacement process. For example, in one displacement process, three fluids, oil (oil phase), water (liquid phase), and air (gas phase), respectively, are present, and the number of fluid phases is set to three, and the main phase therein is set to be the oil phase. The density and viscosity of each phase of fluid is set according to the experimental requirements, and is specific data during the experiment and will not be described in detail herein.
The interfacial tension coefficient of the fluid is set for the digital model and the viscous flow model type is selected. In the embodiment of the application, the interfacial tension coefficient of the fluid is set, and the oil-water interfacial tension is set to be 0.072N/m and the gas-liquid interfacial tension is set to be 0.02N/m by way of example. And selectThe model is a viscous flow model type.
The calculation module 304 is used for initializing a simulation model and setting solution parameters to perform simulation calculation. In the embodiment of the application, the height proportion of the different-phase fluid in the view of the physical model is calculated respectively, and then the space height occupied by the different-phase fluid in the simulation model is calculated according to the total height of the simulation model. And calculating Y coordinates of phase interfaces of the different phase fluids according to the spatial distribution condition of the different phase fluids and the space height occupied by the different phase fluids in the simulation model. In the embodiment of the application, the spatial distribution condition of different phase fluids is that water is at the bottom layer, oil is at the middle layer and air is at the upper layer. The Y coordinate at the top of the space where the water is located is equal to the Y coordinate at the bottom of the space where the oil is located, specifically equal to the Y coordinate at the bottom of the space where the water is located plus the space height of the space where the water is located. Similarly, the Y coordinate of the top of the space where the oil is located is equal to the Y coordinate of the bottom of the space where the air is located, and specifically equal to the Y coordinate of the bottom of the space where the water is located plus the space height of the space where the oil is located. The Y coordinate of the top of the space where the air is located is the Y coordinate of the bottom of the space where the water is located plus the space height of the space where the oil is located plus the space height of the space where the air is located.
Global initialization is set, the initial phase of the simulation model is set as the water phase, i.e. the volume fraction of the initial water phase is set as 1. And then carrying out local initialization setting, setting the volume fraction of the region where the oil is located as 1 according to the spatial distribution condition of the different phase fluids, and setting the volume fraction of the region where the water is located as 1.
The solving parameters comprise a pressure velocity coupling scheme, a gradient space discrete type, a pressure space discrete type, a momentum space discrete type, a volume fraction space discrete type, a turbulence kinetic energy space discrete type, a specific dissipation ratio space discrete type, an inlet boundary condition, an outlet boundary condition, a time step number and an iteration number of the time step.
In the embodiment of the present application, the pressure-velocity coupling scheme is set to be the Coupled scheme. The gradient spatial discrete type is set to Least Squarries Cell Based. The pressure space discrete type is set to prest. The momentum space discrete type is set as Second Order Upwind. The volume fraction spatial discrete type is set to compression. The turbulent kinetic energy space discrete type is set as Second Order Upwind. The specific dissipation ratio spatial discrete type is set to Second Order Upwind. And inlet boundary conditions are set according to experimental injection speed and outlet boundary conditions are set according to experimental outlet conditions. Setting the time step and the time step number to enable the Kurthan numberThe following formula is obtained:
. In the method, in the process of the invention,for the time step size of the time step,for the resolution of the grid it is advantageous,for the purpose of experimental injection rates,for the number of kulange's,in order to make the number of time steps,is the total duration of the experiment. And the number of iterations per time step is illustratively set to 20. It should be noted that the above setting of the solution parameters is only an example of the embodiments of the present application, and is not intended to limit the scope of the present application, and those skilled in the art may make various modifications and changes to the solution parameters according to the actual usage scenario, which still falls within the scope of the present application.
Some of the modules of the apparatus described herein may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The apparatus or module set forth in the embodiments of the application may be implemented in particular by a computer chip or entity, or by a product having a certain function. For convenience of description, the above devices are described as being functionally divided into various modules, respectively. The functions of the modules may be implemented in the same piece or pieces of software and/or hardware when implementing the embodiments of the present application. Of course, a module that implements a certain function may be implemented by a plurality of sub-modules or a combination of sub-units.
The methods, apparatus or modules described herein may be implemented in computer readable program code means and in any suitable manner, for example, the controller may take the form of, for example, a microprocessor or processor and a computer readable medium storing computer readable program code (e.g., software or firmware) executable by the (micro) processor, logic gates, switches, application specific integrated circuits (english: application Specific Integrated Circuit; abbreviated: ASIC), programmable logic controllers and embedded microcontrollers, examples of which include, but are not limited to, the following microcontrollers: ARC 625D, atmel AT91SAM, microchip PIC18F26K20, and Silicone Labs C8051F320, the memory controller may also be implemented as part of the control logic of the memory. Those skilled in the art will also appreciate that, in addition to implementing the controller in a pure computer readable program code, it is well possible to implement the same functionality by logically programming the method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers, etc. Such a controller can be regarded as a hardware component, and means for implementing various functions included therein can also be regarded as a structure within the hardware component. Or even means for achieving the various functions may be regarded as either software modules implementing the methods or structures within hardware components.
The embodiment of the application also provides equipment, which comprises: a processor; a memory for storing processor-executable instructions; the processor, when executing the executable instructions, implements a method as described in embodiments of the present application.
The embodiments also provide a non-transitory computer readable storage medium having stored thereon a computer program or instructions which, when executed, cause a method as described in the embodiments of the present application to be implemented.
In addition, each functional module in the embodiments of the present invention may be integrated into one processing module, each module may exist alone, or two or more modules may be integrated into one module.
The storage medium includes, but is not limited to, a random access Memory (English: random Access Memory; RAM), a Read-Only Memory (ROM), a Cache Memory (English: cache), a Hard Disk (English: hard Disk Drive; HDD), or a Memory Card (English: memory Card). The memory may be used to store computer program instructions.
From the description of the embodiments above, it will be apparent to those skilled in the art that the present application may be implemented in software plus necessary hardware. Based on such understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art in the form of a software product, or may be embodied in the implementation of data migration. The computer software product may be stored on a storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., comprising instructions for causing a computer device (which may be a personal computer, mobile terminal, server, or network device, etc.) to perform the methods described in various embodiments or portions of embodiments herein.
In this specification, each embodiment is described in a progressive manner, and the same or similar parts of each embodiment are referred to each other, and each embodiment is mainly described as a difference from other embodiments. All or portions of the present application can be used in a number of general purpose or special purpose computer system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet devices, mobile communication terminals, multiprocessor systems, microprocessor-based systems, programmable electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The above embodiments are only for illustrating the technical solution of the present application, and not for limiting the present application; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced with equivalents; such modifications and substitutions do not depart from the spirit of the corresponding technical solutions.

Claims (10)

1. A digital experimental method for a physical model of a hole oil reservoir is characterized by comprising the following steps:
constructing a digital model of the physical model;
gridding the digital model to generate an unstructured grid file;
introducing the unstructured grid file into fluid mechanics simulation software, and setting simulation parameters to generate a simulation model;
initializing the simulation model, and setting solving parameters to perform simulation calculation.
2. The method of claim 1, wherein said constructing a digital model of a physical model comprises:
inserting the view of the physical model into three-dimensional software for copying modeling, and adjusting the model to be consistent with the size of the physical model;
and adjusting the width of a crack channel of the model according to the effective volume of the physical model until the volume of the model is close to the physical model to obtain the digital model.
3. The method of claim 1, wherein the meshing the digital model to generate an unstructured grid file comprises:
determining the gridding grid type and the target grid number of the digital model;
calculating the grid resolution of the digital model;
and generating and outputting an unstructured grid file of the digital model according to the grid type, the target grid number and the grid resolution.
4. A method according to claim 3, wherein the calculation formula for calculating the grid resolution of the digital model is as follows:
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->For the grid resolution, V is the volume of the digital model, w is the thickness of the digital model, +.>And the target grid number.
5. The method of claim 3, wherein the meshing the digital model to generate an unstructured grid file further comprises:
checking the dividing quality of the unstructured grid file;
and if the negative volume appears in the unstructured grid file, the digital model is gridded again until the generated unstructured grid file has no negative volume.
6. The method of claim 1, wherein the setting simulation parameters comprises:
setting gravity acceleration for the digital model, and selecting a model type;
setting the number of fluid phases in the digital model, and setting the density and viscosity of each phase of fluid;
and setting an interfacial tension coefficient of the fluid for the digital model, and selecting a viscous flow model type.
7. The method of claim 1, wherein the solution parameters include a pressure velocity coupling scheme, a gradient spatial discrete type, a pressure spatial discrete type, a momentum spatial discrete type, a volume fraction spatial discrete type, a turbulence kinetic energy spatial discrete type, a specific dissipation ratio spatial discrete type, an inlet boundary condition, an outlet boundary condition, a time step, a number of time steps, and a number of iterations of the time step.
8. The utility model provides a hole oil reservoir physical model digital experiment device which is characterized in that includes:
the construction module is used for constructing a digital model of the physical model;
the generation module is used for meshing the digital model to generate an unstructured grid file;
the simulation module is used for guiding the unstructured grid file into fluid mechanics simulation software and setting simulation parameters to generate a simulation model;
and the calculation module is used for initializing the simulation model and setting solving parameters to perform simulation calculation.
9. An apparatus for performing a hole reservoir physical model digitization experimental method, comprising:
a processor;
a memory for storing processor-executable instructions;
the processor, when executing the executable instructions, implements the method of any one of claims 1 to 7.
10. A non-transitory computer readable storage medium comprising instructions for storing a computer program or instructions which, when executed, cause the method of any one of claims 1 to 7 to be implemented.
CN202311532241.6A 2023-11-17 2023-11-17 Digital experimental method and device for fracture-cavity oil reservoir physical model Pending CN117390997A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311532241.6A CN117390997A (en) 2023-11-17 2023-11-17 Digital experimental method and device for fracture-cavity oil reservoir physical model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311532241.6A CN117390997A (en) 2023-11-17 2023-11-17 Digital experimental method and device for fracture-cavity oil reservoir physical model

Publications (1)

Publication Number Publication Date
CN117390997A true CN117390997A (en) 2024-01-12

Family

ID=89468404

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311532241.6A Pending CN117390997A (en) 2023-11-17 2023-11-17 Digital experimental method and device for fracture-cavity oil reservoir physical model

Country Status (1)

Country Link
CN (1) CN117390997A (en)

Similar Documents

Publication Publication Date Title
Larsen et al. Simulation of wave-plus-current scour beneath submarine pipelines
Begnudelli et al. Unstructured grid finite-volume algorithm for shallow-water flow and scalar transport with wetting and drying
Zhao et al. Numerical modeling of local scour below a piggyback pipeline in currents
Shakibaeinia et al. MPS-based mesh-free particle method for modeling open-channel flows
Toombes et al. Numerical limitations of hydraulic models
Huang et al. Well-balanced finite volume scheme for shallow water flooding and drying over arbitrary topography
Gutiérrez et al. Numerical study of Taylor bubbles rising in a stagnant liquid using a level-set/moving-mesh method
GB2506271A (en) Modeling a fractured medium
Kurbanaliev et al. Using OpenFOAM multiphase solver interFoam for large scale modeling
Moreira et al. Dam spillways and the SPH method: Two case studies in Portugal
Qiu Two-dimensional SPH simulations of landslide-generated water waves
Lanetc et al. Coupling of pore network modelling and volume of fluid methods for multiphase flow in fractured media
CN114065593A (en) Method and device for underground water coupling simulation in multi-scale fracture matrix system
Shao et al. Rotation and orientation of irregular particles in viscous fluids using the gradient smoothed method (GSM)
Meng et al. Numerical simulation of gas-liquid two-phase flow impacting fixed structure by CLSVOF/IB method based on OpenFOAM
CN116702648A (en) Vortex-induced vibration calculation method and device for three-dimensional drill string
CN117390997A (en) Digital experimental method and device for fracture-cavity oil reservoir physical model
Lee et al. Development of a numerical method for multiphase flows using an electrostatic model in a wire-mesh sensor
Sabisch et al. 3D volume-of-fluid simulation of a wobbling bubble in a gas-liquid system of low Morton number
Kanfoudi et al. Numerical investigation for steady and unsteady cavitating flows
Abdi et al. Simulations of high-pressure fluid flow in a pre-cracked rock specimen composed of densely packed bonded spheres using a 3D CFD model and simplified 2D coupled CFD-DEM approach
Ma et al. Non-hydrostatic wave model NHWAVE: Documentation and user’s manual (version 2.0)
CN117077575B (en) Method and device for calculating retention time of acid fracturing acid liquor of carbonate reservoir
CN117077574B (en) Quantitative characterization method and device for displacement mechanism of fracture-cavity oil reservoir model
Pelaez A flexible approach to 2D-3D coupling of a shallow-water equation solver to OpenFOAM

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination