Simulation method and device for multi-physical-field oil displacement process and medium

Abstract

The invention discloses a simulation method and device for a multi-physical-field oil displacement process and a medium. The method comprises the steps that a simulation system is determined according to a target area rock core, calculation grids are divided, and the simulation system is pores in a digital rock core; describing the acting force of the electric field on the water phase, the generation rate of the gas phase, the disappearance rate of the gas phase, the inlet condition of the water phase and the physical characteristics of the core in the target area through computational fluid mechanics software according to the computational grid; and after numerical simulation is completed, multiple flow states and recovery rates in the simulation system are determined through data post-processing. The method solves the problems that in the prior art, the effect of electric fields, ultrasonic waves and pulse liquid phase discharge on multiphase flow cannot be described, and the method cannot be used for predicting the oil displacement process under the effect of multiple physical fields, and the accuracy of determining the recovery ratio under the multiple physical fields is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of oilfield development, and in particular to a simulation method, equipment and medium for a multi-physical field oil displacement process. Background Art

[0002] Multiphase flow in porous formation media is crucial for understanding oil displacement mechanisms, optimizing displacement conditions, and developing novel enhanced oil recovery technologies. Numerical simulation is a key tool for revealing multiphase flow in porous media. However, existing numerical simulation methods for oil displacement processes are unable to describe the effects of electric fields, ultrasound, and pulsed liquid discharge on multiphase flow. Consequently, they are unable to predict oil displacement processes under multi-physics fields, making it difficult to provide theoretical support for the development and application of multi-physics enhanced oil recovery technologies. Summary of the Invention

[0003] The present invention provides a simulation method, equipment and medium for a multi-physical field oil recovery process to solve the problem in the prior art that the effects of electric fields, ultrasonic waves and pulsed liquid discharge on multiphase flow cannot be described and cannot be used to predict the oil recovery process under the action of multi-physical fields.

[0004] According to one aspect of the present invention, a method for simulating a multi-physics field oil recovery process is provided, comprising:

[0005] Determine a simulation system based on the target area core and divide the computational grid. The simulation system is the pores in the digital core. The digital core is a model obtained by 3D modeling the target area core.

[0006] Describing the force of the electric field on the water phase, the generation rate of the gas phase, the disappearance rate of the gas phase, the water phase inlet conditions, and the physical properties of the target area core using computational fluid dynamics software based on the computational grid, wherein the water phase inlet conditions are used to characterize the initial velocity of the pulsating water flow, and the physical properties of the target area core at least include: initial oil-water-gas volume fraction distribution, initial pressure, and initial oil mass within the target area core;

[0007] After the numerical simulation is completed, multiple flow states and recovery factors in the simulation system are determined through data post-processing. The multiphase flow state includes at least: the oil-water-gas volume fraction, velocity, pressure and cumulative outflow oil mass inside and at the outlet of the porous medium.

[0008] According to another aspect of the present invention, an electronic device is provided, comprising:

[0009] at least one processor; and

[0010] a memory communicatively connected to the at least one processor; wherein,

[0011] The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the simulation method of the multi-physics field oil recovery process described in any embodiment of the present invention.

[0012] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the simulation method of the multi-physics field oil recovery process described in any embodiment of the present invention when executed.

[0013] The technical solution of the embodiment of the present invention determines a simulation system based on the target area core and divides the computational grid into pores within the digital core. Based on the computational grid, computational fluid dynamics software is used to describe the electric field's effect on the water phase, the gas phase generation rate, the gas phase disappearance rate, the water phase inlet conditions, and the physical properties of the target area core. After the numerical simulation is completed, multiple flow states and recovery factors within the simulation system are determined through data post-processing. This method addresses the problem of existing technologies that cannot describe the effects of electric fields, ultrasonic waves, and pulsed liquid discharges on multiphase flow, making it unsuitable for predicting oil displacement processes under multi-physical fields. This method improves the accuracy of recovery factor determination under multi-physical fields.

[0014] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it intended to limit the scope of the present invention. Other features of the present invention will become readily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0016] Figure 1 A flowchart of a method for simulating a multi-physics field oil recovery process provided by an embodiment of the present invention;

[0017] Figure 2 A schematic structural diagram of an electronic device for implementing the multi-physics field oil recovery process simulation method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0018] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of the present invention.

[0019] It should be noted that the terms "first", "second", etc. in the description and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that the numbers used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0020] Figure 1 This is a flow chart of a method for simulating a multi-physics field oil recovery process provided by an embodiment of the present invention. This embodiment is applicable to the case of simulating a multi-physics field oil recovery process. The method can be executed by a simulation device for a multi-physics field oil recovery process. The simulation device for a multi-physics field oil recovery process can be implemented in the form of hardware and / or software. The simulation device for a multi-physics field oil recovery process can be configured in any electronic device with network communication function. Figure 1 As shown, the method includes:

[0021] S110. Determine the simulation system based on the core of the target area and divide the calculation grid.

[0022] The simulation system is the pores in the digital core, and the digital core is a model obtained by 3D modeling of the core in the target area.

[0023] Furthermore, the creation of digital cores can be divided into two main categories: physical testing methods and numerical reconstruction methods. Physical testing methods involve using testing instruments, such as high-magnification optical microscopes or X-ray CT scanners, to photograph or scan core samples to obtain a large number of two-dimensional core images. These images are then overlaid and reconstructed into a three-dimensional digital core using modeling programs or software. These methods primarily include serial slice imaging, laser scanning confocal microscopy, and X-ray CT scanning. Numerical reconstruction methods, based on a small number of two-dimensional thin-section images, utilize the information contained in these images through stochastic simulation or sedimentary rock process simulation to reconstruct a three-dimensional digital core.

[0024] The pores of the target area core are determined according to the rock structure of the target area core and the digital core modeling method, and the obtained pore model is meshed.

[0025] Optionally, a simulation system is determined based on the core of the target area and a computational grid is divided, including steps A1-A3:

[0026] Step A1: Create a digital core based on the target area core.

[0027] A physical test method or a data reconstruction method is used to create a three-dimensional digital core of the target area core based on the target area core.

[0028] Step A2: extract the pores in the digital core to obtain a simulation system.

[0029] Pore ​​structure refers to the type, size, distribution, and interconnectedness of pores and throats within a rock. A rock's pore system consists of two components: pores and throats. Pores are the bulky part of the system, while the smaller parts connecting the pores are called throats.

[0030] The pores in the digital core are extracted based on the established digital core to obtain a simulation system.

[0031] Furthermore, pores can be extracted by scanning the image of the digital core model through CT scanning, and the pores can be extracted based on the color difference between the pores and the rock entity.

[0032] Step A3: divide the simulation system into computational grids.

[0033] Meshing can divide the model into many small units. As the top priority of numerical simulation pre-processing, the degree of matching between meshing and calculation objectives and the quality of the mesh determine the quality of subsequent simulation calculations.

[0034] The obtained simulation system is divided into grids.

[0035] Furthermore, ICEM or OPENFOAM can be used to divide the computational grid.

[0036] S120. Describe the force of the electric field on the water phase, the generation rate of the gas phase, the disappearance rate of the gas phase, the water phase inlet conditions, and the physical properties of the core in the target area using computational fluid dynamics software based on the computational grid.

[0037] The inlet conditions of the water phase are used to characterize the initial velocity of the pulsating water flow. The physical properties of the target area core include at least: the initial oil, water and gas volume fraction distribution, initial pressure and initial oil quality in the target area core.

[0038] The computational grid is imported into the computational fluid dynamics software. The simulation algorithm configured in the computational fluid dynamics software is used to simulate the force of the electric field on the water phase, the generation rate of the gas phase, the disappearance rate of the gas phase, the water phase inlet conditions and the physical properties of the core in the target area based on the computational grid.

[0039] Furthermore, computational fluid dynamics software may include ANSYS Fluent, cfx, STAR-CCM, comsol, OpenFOAM, Phoenics, and the like.

[0040] Optionally, the force exerted by the electric field on the water phase and the rate of generation and disappearance of the gas phase are described using computational fluid dynamics software based on the computational grid, including steps B1-B4:

[0041] Step B1: Import the computational grid into computational fluid dynamics software and use the fluid volume method to establish a basic mathematical model.

[0042] Among them, the basic mathematical model is used to represent the data information in the computing grid.

[0043] The volume of fluid method (VoF) is an interface tracking method based on an Eulerian grid. In this method, mutually incompatible fluid components share a set of momentum equations, and the phase volume fraction is introduced to track the interfaces between phases within the computational domain.

[0044] The computational grid is imported into computational fluid dynamics software, and the fluid volume method is used to simulate the data information in the computational grid to obtain the basic mathematical model of the core of the target area represented by the computational grid.

[0045] Step B2: Based on the basic mathematical model, the Helmholtz-Smoluchowski theory is used to describe the force acting on the water phase by changing the potential gradient and setting the slip velocity at the wall.

[0046] The Helmholtz-Smoluchowski theory is used to describe the force of the electric field on the water phase by setting the slip velocity at the wall, that is, to drive the flow of water for conventional oil displacement.

[0047] Step B3: Use the Schnerr-Sauer cavitation model to describe the generation and disappearance rates of the gas phase under the action of ultrasound.

[0048] Cavitation is a process in which sound waves in a fluid are actually pressure waves. Mechanical vibrations in a liquid cause the internal pressure to change. When the pressure decreases, structural fractures occur inside the fluid or at some fluid-solid interfaces, forming cavities.

[0049] The Schnerr-Sauer cavitation model is mainly based on the dynamic behavior of liquid in supersonic airflow, and its basic principles mainly include: cavitation phenomenon, interface dynamics, and control equations. Among them, the cavitation phenomenon refers to the liquid cavitation phenomenon caused by the acceleration and shearing of the liquid in the high-speed airflow in the supersonic jet. The interface dynamics can be The Schnerr-Sauer cavitation model describes the behavior of droplets in supersonic airflow by considering the movement and deformation of the gas-liquid interface. The control equation can be The control equation of the Schnerr-Sauer cavitation model is a mathematical equation that describes the cavitation phenomenon and droplet motion. This equation can quantitatively describe the cavitation behavior and dynamic characteristics of droplets in supersonic jets.

[0050] Step B4: Determine the inlet conditions of the water phase through an indoor test of liquid-phase discharge plasma rock breaking.

[0051] The indoor test of liquid phase discharge plasma rock breaking is a new and efficient rock breaking technology that converts electrical energy into shock wave mechanical energy.

[0052] In the indoor test environment of liquid phase discharge plasma rock breaking, the pulse discharge voltage, frequency and pulse width are changed. The initial velocity of the pulsating water flow caused by the pulse discharge is determined through experiments. The relationship between the obtained initial velocity and the pulse discharge voltage, frequency and pulse width is fitted into a mathematical expression as the inlet condition of the water phase.

[0053] Optionally, the inlet conditions of the water phase are determined by a liquid phase discharge plasma rock breaking indoor test, including steps C1-C2:

[0054] Step C1: Determine the effect of pulse discharge voltage, pulse frequency and pulse width on the initial velocity of the pulsating water flow induced by the pulse discharge through a liquid phase discharge plasma rock breaking indoor test.

[0055] In the liquid phase discharge plasma rock breaking indoor test environment, the initial velocity of the pulsating water flow caused by the pulse discharge when breaking the target area core is measured by changing the pulse discharge voltage, pulse frequency and pulse width.

[0056] Step C2: Fitting the initial velocity of the pulsating water flow induced by the pulse discharge with the pulse discharge voltage, pulse frequency and pulse width to obtain the inlet conditions of the water phase.

[0057] The pulse discharge voltage, pulse frequency and pulse width are fitted to the initial velocity of the pulsating water flow caused by the pulse discharge, and the fitted mathematical expression is used as the inlet condition of the water phase.

[0058] S130. After the numerical simulation is completed, multiple flow states and recovery factors in the simulation system are determined through data post-processing.

[0059] The multiphase flow state includes at least: the oil-water-gas volume fraction, velocity, pressure and cumulative outflow oil mass inside the porous medium and at the outlet.

[0060] Numerical simulation is performed on a computer based on the force of the water phase, the generation rate of the gas phase, the disappearance rate of the gas phase, the water phase inlet conditions and the initial physical properties of the core in the target area. The simulation results are analyzed to obtain multiple flow states in the simulation system and calculate the recovery rate.

[0061] Optionally, after the numerical simulation is completed, multiple flow states and recovery factors in the simulation system are determined through data post-processing, including:

[0062] The recovery factor is determined based on the cumulative outflow oil mass and the initial oil mass of the core in the target area.

[0063] The recovery factor is the ratio of the amount of crude oil produced to the original geological reserves of the reservoir, usually expressed as a percentage.

[0064] The recovery factor is determined based on the cumulative outflow oil mass obtained from the simulation and the initial oil mass of the core in the target area.

[0065] Optionally, the recovery factor is determined based on the cumulative outflow oil mass and the initial oil mass of the core in the target area, including:

[0066] The recovery factor is equal to the cumulative outflow oil mass divided by the initial oil mass of the core in the target area.

[0067] The initial oil mass of the target area core may be the mass of the oil remaining in the pores of the target area core.

[0068] Optionally, after the numerical simulation is completed, after multiple flow states and recovery factors in the simulation system are determined through data post-processing, steps D1-D2 are included:

[0069] Step D1: By changing the potential gradient, ultrasonic pulse frequency, pulse discharge voltage, pulse frequency and pulse width, multiple flow states and recovery factors in the simulation system under different conditions are obtained.

[0070] By continuously changing the potential gradient, ultrasonic pulse frequency, pulse discharge voltage, pulse frequency and pulse width, multiple flow states and recovery rates in the simulation system under different conditions can be obtained.

[0071] Step D2: Optimize the displacement parameters using the optimization model to obtain displacement parameters that enable the recovery factor to meet preset requirements.

[0072] Among them, oil displacement parameters include: potential gradient, ultrasonic pulse frequency, pulse discharge voltage, pulse frequency and pulse width.

[0073] The optimization model is used to optimize the potential gradient, ultrasonic pulse frequency, pulse discharge voltage, pulse frequency and pulse width to obtain the oil displacement parameters that meet the preset requirements of the recovery factor.

[0074] Furthermore, the optimization model may adopt response surface or genetic algorithm.

[0075] The technical solution of this embodiment determines a simulation system based on the target area core and divides the computational grid into pores within the digital core. Based on the computational grid, computational fluid dynamics software is used to describe the electric field's effect on the water phase, the gas generation rate, the gas phase disappearance rate, the water phase inlet conditions, and the physical properties of the target area core. After the numerical simulation is complete, data post-processing is used to determine multiple flow states and recovery factors within the simulation system. This method addresses the existing problem of being unable to describe the effects of electric fields, ultrasonic waves, and pulsed liquid discharges on multiphase flow, making it incapable of predicting oil displacement processes under multi-physical fields. This method improves the accuracy of recovery factor determination under multi-physical fields.

[0076] Figure 2 The structural representation of the electronic device of the simulation method for the multi-physics field flooding process for realizing the embodiment of the present invention.The electronic device is intended to represent various forms of digital computers, such as, laptop computers, desktop computers, workbenches, personal digital assistants, servers, blade servers, mainframe computers and other applicable computers.The electronic device can also represent various forms of mobile devices, such as, personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches etc.) and other similar computing devices.The components shown herein, their connection and relationship and their functions are merely examples, and are not intended to limit the realization of the present invention described herein and / or required.

[0077] like Figure 2 As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., which is communicatively connected to the at least one processor 11. The memory stores a computer program that can be executed by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. Various programs and data required for the operation of the electronic device 10 can also be stored in the RAM 13. The processor 11, ROM 12, and RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0078] Multiple components in the electronic device 10 are connected to the I / O interface 15, including an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0079] The processor 11 can be a variety of general-purpose and / or specialized processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the simulation method of the multi-physics field flooding process.

[0080] In some embodiments, the simulation method of the multi-physics flooding process can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the simulation method of the multi-physics flooding process described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to perform the simulation method of the multi-physics flooding process by any other appropriate means (for example, by means of firmware).

[0081] Various embodiments of the systems and techniques described herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system comprising at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0082] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program may be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0083] In the context of the present invention, computer-readable storage media can be tangible media that can contain or store a computer program for use with an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. Computer-readable storage media can include but are not limited to electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, computer-readable storage media can be machine-readable signal media. More specific examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0084] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0085] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0086] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0087] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in the present invention can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved. This is not limited herein.

[0088] The above specific embodiments do not limit the scope of protection of the present invention. Those skilled in the art will appreciate that various modifications, combinations, sub-combinations, and substitutions may be made based on design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention are intended to be included within the scope of protection of the present invention.

Claims

1. A method for simulating a multi-physics field flooding process, characterized in that: include: Determine a simulation system based on the target area core and divide the computational grid. The simulation system is the pores in the digital core. The digital core is a model obtained by 3D modeling the target area core. Describing the force of the electric field on the water phase, the generation rate of the gas phase, the disappearance rate of the gas phase, the water phase inlet conditions, and the physical properties of the target area core using computational fluid dynamics software based on the computational grid, wherein the water phase inlet conditions are used to characterize the initial velocity of the pulsating water flow, and the physical properties of the target area core at least include: initial oil-water-gas volume fraction distribution, initial pressure, and initial oil mass within the target area core; After the numerical simulation is completed, multiple flow states and recovery factors in the simulation system are determined through data post-processing. The multiphase flow state includes at least: the oil-water-gas volume fraction, velocity, pressure and cumulative outflow oil mass inside and at the outlet of the porous medium.

2. The method according to claim 1, characterized in that Determine the simulation system based on the core of the target area and divide the computational grid, including: Create digital cores based on target area cores; Extracting pores in the digital core to obtain the simulation system; The simulation system is divided into computational grids.

3. The method according to claim 1, characterized in that The force of the electric field on the water phase and the rate of generation and disappearance of the gas phase are described using computational fluid dynamics software based on the computational grid, including: Importing the computational grid into computational fluid dynamics software, and using the fluid volume method to establish a basic mathematical model, wherein the basic mathematical model is used to characterize the data information in the computational grid; According to the basic mathematical model, the Helmholtz-Smoluchowski theory is used to describe the force acting on the water phase by changing the potential gradient and setting the slip velocity at the wall; The Schnerr-Sauer cavitation model is used to describe the generation and disappearance rates of the gas phase under ultrasonic action; The inlet conditions of the water phase were determined through indoor experiments on liquid phase discharge plasma rock breaking.

4. The method according to claim 3, characterized in that The inlet conditions of the water phase are determined through indoor tests of liquid phase discharge plasma rock breaking, including: The effects of pulse discharge voltage, pulse frequency and pulse width on the initial velocity of pulsating water flow caused by pulse discharge are determined by the liquid phase discharge plasma rock breaking indoor test; The pulse discharge voltage, pulse frequency and pulse width are fitted to the initial velocity of the pulsating water flow caused by the pulse discharge to obtain the inlet conditions of the water phase.

5. The method according to claim 1, wherein After the numerical simulation is completed, data post-processing is used to determine multiple flow states and recovery factors in the simulation system, including: The recovery factor is determined based on the cumulative outflow oil mass and the initial oil mass of the core in the target area.

6. The method according to claim 5, characterized in that The recovery factor is determined based on the cumulative outflow oil quality and the initial oil quality of the core in the target area, including: The recovery factor is equal to the cumulative outflow oil mass divided by the initial oil mass of the target area core.

7. The method according to claim 1, characterized in that After the numerical simulation is completed, multiple flow states and recovery factors in the simulation system are determined through data post-processing, including: By changing the potential gradient, the ultrasonic pulse frequency, the pulse discharge voltage, the pulse frequency and the pulse width, multiple flow states and recovery factors in the simulation system under different conditions are obtained; The oil displacement parameters are optimized using an optimization model to obtain oil displacement parameters that enable the recovery factor to meet preset requirements. The oil displacement parameters include: the potential gradient, the ultrasonic pulse frequency, the pulse discharge voltage, the pulse frequency, and the pulse width.

8. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor to enable the at least one processor to perform the multi-physics field oil recovery process simulation method according to any one of claims 1 to 7.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the simulation method of the multi-physical field oil recovery process according to any one of claims 1 to 7 when executed.

Images