CN117854605B - Method, system, equipment and storage medium for simulating viscous finger-in phenomenon - Google Patents
Method, system, equipment and storage medium for simulating viscous finger-in phenomenon Download PDFInfo
- Publication number
- CN117854605B CN117854605B CN202410259203.6A CN202410259203A CN117854605B CN 117854605 B CN117854605 B CN 117854605B CN 202410259203 A CN202410259203 A CN 202410259203A CN 117854605 B CN117854605 B CN 117854605B
- Authority
- CN
- China
- Prior art keywords
- fracturing fluid
- simulation time
- calculation
- space
- fluid
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000003860 storage Methods 0.000 title claims abstract description 18
- 239000012530 fluid Substances 0.000 claims abstract description 407
- 238000004088 simulation Methods 0.000 claims abstract description 250
- 239000002245 particle Substances 0.000 claims abstract description 219
- 238000004364 calculation method Methods 0.000 claims abstract description 193
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims description 24
- 238000006467 substitution reaction Methods 0.000 claims description 12
- 238000005086 pumping Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- 238000009877 rendering Methods 0.000 claims description 6
- 238000005325 percolation Methods 0.000 claims description 2
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 5
- 230000001413 cellular effect Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 5
- 238000012935 Averaging Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention provides a method, a system, equipment and a storage medium for simulating a viscous finger-in phenomenon, which are applied to the technical field of computer simulation, wherein the method comprises the following steps: grid division is carried out on a model of a fracturing fluid flowing region in a target stratum, a plurality of calculation spaces are obtained, and under the condition that the current simulation time is not the termination simulation time, each calculation space is divided into: according to the velocity field and the fracturing fluid parameters of each type of fracturing fluid particles of the computing space and the adjacent computing space under the adjacent historical simulation time, calculating the position parameters and the velocity field of each type of fracturing fluid particles in the computing space under the current simulation time, and under the condition that the current simulation time is the termination simulation time, generating a simulation result of the viscous finger-in phenomenon based on the position parameters of each type of fracturing fluid particles of each computing space under the adjacent historical simulation time of the termination simulation time. The invention realizes the accurate simulation of the viscous finger-in phenomenon of the fracturing fluid in the acidizing and fracturing process.
Description
Technical Field
The present invention relates to the field of computer simulation technologies, and in particular, to a method, a system, an apparatus, and a storage medium for simulating a viscous finger-in phenomenon.
Background
The viscous fingering (viscous fingering) phenomenon refers to a phenomenon that during acidizing fracturing (Fracturing Acidizing), the physical properties of various fracturing fluids injected into the fracture are different, so that the various fracturing fluids generate finger-like interfacial instability during the flowing process in the fracture. The phenomenon of viscous fingering can cause uneven distribution of various fracturing fluids injected into cracks, affect the corrosion effect on the wall surfaces of the cracks, and form the cracks with severely uneven width, thereby affecting the subsequent production. Therefore, it is extremely important to simulate the viscous fingering of fracturing fluids during acidizing fracturing.
Disclosure of Invention
The embodiment of the invention aims to provide a method, a system, equipment and a storage medium for simulating the viscous fingering phenomenon, so as to realize accurate simulation of the viscous fingering phenomenon of fracturing fluid in the acidizing and fracturing process. The specific technical scheme is as follows:
a method of simulating a viscous finger-in phenomenon, the method comprising:
Performing grid division on a model of a fracturing fluid flowing region in a target stratum to obtain a plurality of calculation spaces;
In the case where the current simulation time is not the termination simulation time, for each calculation space: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation time according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the fracturing fluid parameters of the adjacent historical simulation time of the calculation space, wherein the adjacent historical simulation time is adjacent to the current simulation time, and the time value is earlier than one simulation time of the current simulation time;
and under the condition that the current simulation time is the termination simulation time, generating a simulation result of the viscous finger-in phenomenon of each type of fracturing fluid particle in the fracturing fluid flowing area based on the position parameters of each type of fracturing fluid particle of each calculation space at the adjacent historical simulation time of the termination simulation time.
Optionally, in the case where the current simulation time is not the termination simulation time, for each computation space: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the adjacent historical simulation moment of the calculation space and the fracturing fluid parameters, wherein the method comprises the following steps:
Performing particle calibration processing on each computing space, wherein the particle calibration processing comprises: calculating the position of each fracturing fluid particle at the current simulation moment according to the velocity field of each type of fracturing fluid particle of the calculation space at the adjacent historical simulation moment, deleting the fracturing fluid particle of which the position is not in the calculation space from the calculation space, and sending the type and the velocity field of the deleted fracturing fluid particle to the calculation space adjacent to the calculation space where the position is located;
For each of the computation spaces after the particle calibration process: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid.
Optionally, the calculating, according to the velocity field and the fracturing fluid parameter of each type of fracturing fluid particles at the adjacent historical simulation time by using a continuity equation of a preset fluid and a seepage darcy equation of the preset fluid, the position parameter and the velocity field of each type of fracturing fluid particles at the current simulation time in the calculation space includes:
For each of the fracturing fluid particles of the target type:
Substituting the fluid density rho in the fluid parameters of the fracturing fluid of the target type and the velocity field of the fracturing fluid particles at the adjacent historical simulation time t-1 into the continuity equation of the preset fluid to obtain the continuity equation of the preset fluid subjected to parameter substitution:
，
wherein t is the current simulation time, the Being the partial derivative sign, b being the crack width of the computation space, said/>For the velocity of the fracturing fluid particles in the x-axis direction of the computation space in the velocity field, the/>For the velocity of the fracturing fluid particles in the z-axis direction of the computation space in the velocity field, the/>For fluid loss in the fracturing fluid parameters of the target type, the/>Preset pumping speed of the fracturing fluid particles of the target type;
substituting the fluid density ρ, the fluid pressure p and the fluid viscosity μ in the fluid parameters of the fracturing fluid of the target type into a seepage darcy equation of the preset fluid to obtain a seepage darcy equation of the preset fluid subjected to parameter substitution:
，
Wherein the said For the velocity of the target fracturing fluid particle in the x-axis direction of the calculation space at the current simulation time t to be solved, the/>For the velocity of the target fracturing fluid particle in the z-axis direction of the calculation space at the current simulation time t to be solved, the/>The K is the crack permeability of the calculation space, and the g is the gravity acceleration;
simultaneously solving a continuity equation of the preset fluid substituted by the parameters and a seepage darcy equation of the preset fluid to obtain a model comprising the parameters And said/>Is based on the velocity field of the fracturing fluid particles at the current simulation time and based on the/>And said/>And obtaining the position parameters of the fracturing fluid particles in the calculation space.
Optionally, the computing space includes a plurality of grids, and when the current simulation time is a termination simulation time, generating a simulation result of a viscous fingering phenomenon of each type of fracturing fluid particle in the fracturing fluid flowing area based on the position parameter of each type of fracturing fluid particle of each computing space at an adjacent historical simulation time of the termination simulation time includes:
For each computation space: determining the distribution condition of the fracturing fluid points in each grid of the calculation space according to the position parameters of the fracturing fluid points of each type; rendering the grid into a display state corresponding to the type under the condition that the distribution condition characterizes that only one type of fracturing fluid particles exist in the grid; under the condition that the distribution condition represents that at least two types of fracturing fluid particles exist in the grid, interpolation average is carried out based on each type, and the grid is rendered into a display state corresponding to the interpolation average result; splicing the display states of the grids of the calculation space to obtain the display state of the calculation space;
And splicing the display states of the calculation spaces to generate simulation results of the viscous fingering phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing area.
A system for simulating the phenomenon of viscous finger-in, the system comprising:
the grid division module is used for carrying out grid division on a model of a fracturing fluid flowing region in the target stratum to obtain a plurality of calculation spaces;
The parameter calculation module is used for calculating the space for each calculation when the current simulation time is not the termination simulation time: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation time according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the fracturing fluid parameters of the adjacent historical simulation time of the calculation space, wherein the adjacent historical simulation time is adjacent to the current simulation time, and the time value is earlier than one simulation time of the current simulation time;
And the result generation module is used for generating simulation results of the viscous finger-in phenomenon of each type of fracturing fluid particles in the fracturing fluid flowing area based on the position parameters of each type of fracturing fluid particles of each calculation space at the adjacent historical simulation time of the termination simulation time under the condition that the current simulation time is the termination simulation time.
Optionally, the parameter calculation module is configured to:
Performing particle calibration processing on each computing space, wherein the particle calibration processing comprises: calculating the position of each fracturing fluid particle at the current simulation moment according to the velocity field of each type of fracturing fluid particle of the calculation space at the adjacent historical simulation moment, deleting the fracturing fluid particle of which the position is not in the calculation space from the calculation space, and sending the type and the velocity field of the deleted fracturing fluid particle to the calculation space adjacent to the calculation space where the position is located;
For each of the computation spaces after the particle calibration process: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid.
Optionally, the parameter calculation module calculates, according to the velocity field and the fracturing fluid parameters of each type of fracturing fluid particles at the adjacent historical simulation time, the position parameters and the velocity fields of each type of fracturing fluid particles at the current simulation time in the calculation space by using a continuity equation of a preset fluid and a seepage darcy equation of the preset fluid, where the parameters are set as follows:
For each of the fracturing fluid particles of the target type:
Substituting the fluid density rho in the fluid parameters of the fracturing fluid of the target type and the velocity field of the fracturing fluid particles at the adjacent historical simulation time t-1 into the continuity equation of the preset fluid to obtain the continuity equation of the preset fluid subjected to parameter substitution:
，
wherein t is the current simulation time, the Being the partial derivative sign, b being the crack width of the computation space, said/>For the velocity of the fracturing fluid particles in the x-axis direction of the computation space in the velocity field, the/>For the velocity of the fracturing fluid particles in the z-axis direction of the computation space in the velocity field, the/>For fluid loss in the fracturing fluid parameters of the target type, the/>Preset pumping speed of the fracturing fluid particles of the target type;
substituting the fluid density ρ, the fluid pressure p and the fluid viscosity μ in the fluid parameters of the fracturing fluid of the target type into a seepage darcy equation of the preset fluid to obtain a seepage darcy equation of the preset fluid subjected to parameter substitution:
，
Wherein the said For the velocity of the target fracturing fluid particle in the x-axis direction of the calculation space at the current simulation time t to be solved, the/>For the velocity of the target fracturing fluid particle in the z-axis direction of the calculation space at the current simulation time t to be solved, the/>The K is the crack permeability of the calculation space, and the g is the gravity acceleration;
simultaneously solving a continuity equation of the preset fluid substituted by the parameters and a seepage darcy equation of the preset fluid to obtain a model comprising the parameters And said/>Is based on the velocity field of the fracturing fluid particles at the current simulation time and based on the/>And said/>And obtaining the position parameters of the fracturing fluid particles in the calculation space.
Optionally, the computation space comprises a plurality of grids, and the result generation module is configured to:
For each computation space: determining the distribution condition of the fracturing fluid points in each grid of the calculation space according to the position parameters of the fracturing fluid points of each type; rendering the grid into a display state corresponding to the type under the condition that the distribution condition characterizes that only one type of fracturing fluid particles exist in the grid; under the condition that the distribution condition represents that at least two types of fracturing fluid particles exist in the grid, interpolation average is carried out based on each type, and the grid is rendered into a display state corresponding to the interpolation average result; splicing the display states of the grids of the calculation space to obtain the display state of the calculation space;
And splicing the display states of the calculation spaces to generate simulation results of the viscous fingering phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing area.
An apparatus for simulating the viscous finger phenomenon, the apparatus comprising:
A processor;
A memory for storing the processor-executable instructions;
Wherein the processor is configured to execute the instructions to implement a method of simulating a viscous fingering phenomenon as described in any of the above.
A computer readable storage medium, which when executed by a processor of a simulation device of a viscous fingering phenomenon, causes the simulation device of a viscous fingering phenomenon to perform a method of simulating a viscous fingering phenomenon as described in any one of the above.
According to the simulation method, system, equipment and storage medium for the viscous finger-in phenomenon, which are provided by the embodiment of the invention, the model of the fracturing fluid flowing region in the target stratum can be configured to be subjected to grid division, and the obtained position parameters and speed fields of various types of fracturing fluid particles in each calculation space are calculated, so that the position parameters representing the positions of various types of fracturing particles in each calculation space at each simulation moment are obtained. And then, under the condition that the current simulation time is the termination simulation time, determining the specific positions of the various types of fracturing fluid particles in the fracturing fluid flowing region in the target stratum according to the position parameters of the various types of fracturing fluid particles in the adjacent historical simulation time of each calculation space, so as to generate the simulation result of the viscous finger-in phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing region. Meanwhile, the method and the device determine the position parameters of the fracturing fluid particles in the calculation space by configuring the velocity field based on the fracturing fluid particles of each type and the parameters with strictly defined physical meanings such as the fracturing fluid parameters, and the calculation space can accurately represent the physical characteristics of the fracturing fluid flowing region in the target stratum, so that the accuracy of the finally obtained simulation result of the viscous fingering phenomenon is improved. Therefore, the invention realizes the accurate simulation of the viscous finger-in phenomenon of the fracturing fluid in the acidizing and fracturing process.
Of course, it is not necessary for any one product or method of practicing the invention to achieve all of the advantages set forth above at the same time.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for simulating the viscous finger-in phenomenon according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of distribution of particles of fracturing fluid at adjacent historical simulation moments provided by an alternative embodiment of the present invention;
FIG. 3 is a schematic diagram of the distribution of particles of a fracturing fluid at a current simulated time provided by an alternative embodiment of the present invention;
FIG. 4 is a flow chart of a method for simulating the viscous finger-in phenomenon according to another alternative embodiment of the present invention;
FIG. 5 is a block diagram of a simulation system for viscous finger-in according to an embodiment of the present invention;
fig. 6 is a block diagram of a simulation apparatus for viscous finger phenomenon according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. 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.
The embodiment of the invention provides a method for simulating a viscous finger-in phenomenon, as shown in fig. 1, which comprises the following steps:
S101, performing grid division on a model of a fracturing fluid flowing region in a target stratum to obtain a plurality of calculation spaces.
It should be noted that, in an actual application scenario, the model of the flowing area of the fracturing fluid in the target stratum may be a model constructed based on the combination of construction data and geological parameters after the target stratum is explored, where the types of the geological parameters include, but are not limited to: formation burial depth, cap layer location, young's modulus of rock, poisson's ratio of rock, gap size, perforation location, perforation depth, etc.
In the actual application scenario, the calculation space may be a partial model obtained by dividing the model of the flow passage area of the fracturing fluid subjected to grid division according to the actual application scenario requirement. As the existing simulation means mostly adopts full-model simulation calculation, the simulation efficiency is reduced due to the excessive calculated amount. The invention calculates and processes the obtained multiple calculation spaces by configuration, and integrates the calculation and processing results of each calculation space, thereby reducing the calculation amount and improving the simulation efficiency.
It will be appreciated by those skilled in the art that in a practical application scenario, the computation space may be implemented by an adaptive mesh subdivision algorithm, and the types of the adaptive mesh subdivision algorithm may be various, for example: an adaptive mesh subdivision algorithm based on error estimation, an adaptive mesh subdivision algorithm based on gradient information, an adaptive mesh subdivision algorithm based on characteristic information, and the like. The invention does not define and describe the specific process of obtaining the calculation space by using the self-adaptive mesh subdivision algorithm too much.
S102, when the current simulation time is not the termination simulation time, performing calculation on each calculation space: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation time according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the fracturing fluid parameters at the adjacent historical simulation time of the calculation space, wherein the adjacent historical simulation time is adjacent to the current simulation time, and the time value is earlier than the simulation time of the current simulation time.
In the actual application scenario, the termination simulation time may be a termination time in the fracturing fluid pumping duration set in the configuration file based on the model.
It should be noted that, in a practical application scenario, the parameters of the fracturing fluid may be parameters of various types of fracturing fluids adopted in the construction process, including but not limited to: mass, density, viscosity, etc.
In the practical application scenario, the velocity field refers to a physical field formed by velocity vectors of the current position in the calculation space in all directions at one simulation time of a fracturing fluid particle. For example, if the computation space is a two-dimensional spatial coordinate system consisting of an x-axis and a y-axis, the velocity field of the fracturing fluid particles in the computation space should include a velocity vector in the x-axis direction and a velocity vector in the y-axis direction.
In the practical application scenario, the simulation of the viscous finger advance is implemented by using a cellular automaton (Cellular Automata, CA). Cellular automata is a grid dynamics model with discrete time, space and state and local space interaction and time causal relationship. However, since cellular automata is not simulated by a strictly defined physical equation or parameter, the accuracy of the simulation result output is low when the viscous finger phenomenon is simulated. The method and the device determine the position parameters of the fracturing fluid particles in the calculation space by configuring the velocity field based on the fracturing fluid particles of each type and the parameters with strictly defined physical meanings such as the fracturing fluid parameters, and the calculation space can accurately represent the physical characteristics of the fracturing fluid flowing region in the target stratum, so that the accuracy of the finally obtained simulation result of the viscous fingering phenomenon is improved.
And S103, under the condition that the current simulation time is the termination simulation time, generating a simulation result of the viscous finger-in phenomenon of each type of fracturing fluid particle in the fracturing fluid flowing area based on the position parameters of each type of fracturing fluid particle of each calculation space at the adjacent historical simulation time of the termination simulation time.
It should be noted that, in an actual application scenario, the above-mentioned position parameter may be a parameter indicating a specific position of the fracturing fluid particles in the corresponding calculation space, for example, a coordinate in a coordinate system where the calculation space is located.
According to the method, the model of the fracturing fluid flowing region in the target stratum is subjected to grid division, and the obtained position parameters and the speed fields of the fracturing fluid particles of all types in each calculation space are calculated, so that the position parameters representing the positions of the fracturing particles of all types in each calculation space at each simulation moment are obtained. And then, under the condition that the current simulation time is the termination simulation time, determining the specific positions of the various types of fracturing fluid particles in the fracturing fluid flowing region in the target stratum according to the position parameters of the various types of fracturing fluid particles in the adjacent historical simulation time of each calculation space, so as to generate the simulation result of the viscous finger-in phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing region. Meanwhile, the method and the device determine the position parameters of the fracturing fluid particles in the calculation space by configuring the velocity field based on the fracturing fluid particles of each type and the parameters with strictly defined physical meanings such as the fracturing fluid parameters, and the calculation space can accurately represent the physical characteristics of the fracturing fluid flowing region in the target stratum, so that the accuracy of the finally obtained simulation result of the viscous fingering phenomenon is improved. Therefore, the invention realizes the accurate simulation of the viscous finger-in phenomenon of the fracturing fluid in the acidizing and fracturing process.
Optionally, in the case where the current simulation time is not the termination simulation time, for each computation space: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the adjacent calculation space of the calculation space at the adjacent historical simulation moment and the fracturing fluid parameters, wherein the method comprises the following steps:
Carrying out particle calibration processing on each calculation space, wherein the particle calibration processing comprises the following steps: calculating the position of each fracturing fluid particle under the current simulation moment according to the velocity field of each type of fracturing fluid particle of the calculation space under the adjacent historical simulation moment, deleting the fracturing fluid particle of which the position is not in the calculation space from the calculation space, and transmitting the type and the velocity field of the deleted fracturing fluid particle to the calculation space adjacent to the calculation space where the position is located;
For each calculation space after particle calibration treatment: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid.
It should be noted that, in a practical application scenario, there are various embodiments of the particle calibration process, and one of the embodiments is provided herein by way of example:
fig. 2 is a schematic diagram of distribution of fracturing fluid particles at adjacent historical simulation moments, and fig. 3 is a schematic diagram of distribution of fracturing fluid particles at a current simulation moment. The four calculation spaces divided by the thickened line segments which are perpendicular to each other are respectively a calculation space 1, a calculation space 2, a calculation space 3 and a calculation space 4, wherein white circles in the figures are fracturing fluid particles of the first type, and black circles are fracturing fluid particles of the second type.
Assume that the location of a fracturing fluid particle of the type a, numbered 201 in fig. 2, at adjacent historic simulation moments is within computation space 3. Assume that the location at the current simulation time calculated based on the velocity field of the fracturing fluid particles 201 at the adjacent historical simulation time is no longer within the computation space 3, but rather within the computation space 4. Then for computation space 3 it is necessary to delete the first type of fracturing fluid particles 201 from computation space 3 and send the type of first type of fracturing fluid particles 201 and the calculated velocity field to computation space 2 so that computation space 2 calculates the position of the first type of fracturing fluid particles 201 in computation space 2 at the current simulation time as shown in fig. 3 from the velocity field.
It will be appreciated by those skilled in the art that in an actual application scenario, the foregoing calculation space may be implemented based on a multi-point interface (Multi Point Interface, MPI) protocol by configuring a plurality of data processing units, setting one data processing unit for processing all calculations of one calculation space, and establishing data transmission for quality inspection of each data processing unit through the MPI protocol. The specific construction process of the MPI protocol and the data processing unit is not excessively limited and repeated.
It should be noted that, in a practical application scenario, the continuity equation of the preset fluid is an expression for representing conservation of mass of the fluid in the flowing process. The above-described percolation darcy equation for a preset fluid is an expression that characterizes the relationship between pressure and velocity of the fluid during flow. Because the key influencing factors influencing the flowing direction of the fracturing fluid are the moving speed of the fracturing fluid in all directions, the method has the advantages that two strictly defined physical equations, namely the continuity equation of the preset fluid and the seepage darcy equation of the preset fluid, are introduced, the position parameters and the speed fields of all types of fracturing fluid particles in the calculation space at the current simulation moment are calculated according to the speed fields and the fracturing fluid parameters of all types of fracturing fluid particles at the adjacent historical simulation moment, and the accurate positioning of the positions of the fracturing fluid particles at the current simulation moment is realized, so that compared with the simulation mode of the prior cellular automaton, the simulation precision is improved.
Optionally, calculating the position parameter and the velocity field of each type of fracturing fluid particle in the calculation space at the current simulation moment according to the velocity field and the fracturing fluid parameter of each type of fracturing fluid particle at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage darcy equation of the preset fluid, including:
For each fracturing fluid particle of the target type:
substituting the fluid density rho in the fluid parameters of the target type fracturing fluid and the velocity field of the fracturing fluid particles at the adjacent historical simulation time t-1 into a continuity equation of the preset fluid to obtain the continuity equation of the preset fluid substituted by the parameters:
，
Where t is the current simulation time instant, Is the partial derivative sign, b is the crack width of the computation space,/>Is the velocity of the fracturing fluid particles in the x-axis direction of the computation space in the velocity field,/>Is the velocity of the fracturing fluid particles in the z-axis direction of the computation space in the velocity field,/>For fluid loss in a target type of fracturing fluid parameter,/>The method comprises the steps of presetting pumping speed for fracturing fluid particles of a target type;
substituting the fluid density rho, the fluid pressure p and the fluid viscosity mu in the fluid parameters of the target type fracturing fluid into a seepage darcy equation of the preset fluid to obtain the seepage darcy equation of the preset fluid subjected to parameter substitution:
，
wherein, For the velocity of the target fracturing fluid particle in the x-axis direction of the calculation space at the current simulation time t to be solved,/>For the velocity of the target fracturing fluid particle in the z-axis direction of the calculation space at the current simulation time t to be solved,/>K is the crack permeability of the calculation space, g is the gravity acceleration;
simultaneously solving a continuity equation of the preset fluid subjected to parameter substitution and a seepage darcy equation of the preset fluid to obtain a model comprising And/>The velocity field of the fracturing fluid particles at the current simulation time is based on/>And/>And obtaining the position parameters of the fracturing fluid particles in the calculation space.
Optionally, the computing space includes a plurality of grids, and when the current simulation time is the termination simulation time, generating a simulation result of a viscous fingering phenomenon of each type of fracturing fluid particle in a fracturing fluid flowing area based on a position parameter of each type of fracturing fluid particle of each computing space at an adjacent historical simulation time of the termination simulation time, where the simulation result includes:
For each computation space: determining the distribution condition of fracturing fluid points in each grid of the calculation space according to the position parameters of the fracturing fluid points of each type; under the condition that the distribution condition characterizes that only one type of fracturing fluid particles exist in the grid, rendering the grid into a display state corresponding to the type; under the condition that the distribution condition represents that at least two types of fracturing fluid particles exist in the grid, interpolation averaging is carried out based on each type, and the grid is rendered into a display state corresponding to the interpolation averaging result; splicing the display states of the grids of the calculation space to obtain the display state of the calculation space;
And splicing the display states of the calculation spaces to generate simulation results of the viscous fingering phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing region.
It should be noted that, in a practical application scenario, there are various embodiments of the above-mentioned simulation method of the viscous finger-in phenomenon as shown in fig. 1, and an exemplary embodiment is provided herein:
As shown in fig. 4, a flowchart of a method for simulating the viscous finger phenomenon is shown, and the specific operation steps are as follows:
And S401, constructing a model of a fracturing fluid flowing region of the target stratum according to the imported configuration file, and performing grid division and calculation space distribution on the model to obtain a plurality of calculation spaces. Step S402 is triggered.
Step S402, judging whether the current simulation time is the termination simulation time, if not, triggering step S403, and if so, triggering step S405.
In step S403, according to the velocity field of each type of fracturing fluid particle at the adjacent historical simulation time, the particle calibration process is performed on each calculation space, and each calculation space after the particle calibration process is performed on each calculation space: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid. And triggers step S404.
Optionally, in an optional embodiment of the present invention, if the current simulation time is the time when the simulation is performed for the first time, the velocity field of each type of fracturing fluid particles at the adjacent simulation time may be selected according to the initial parameters in the configuration file of the model.
In step S404, the current simulation time is updated to the next simulation time of the current simulation time, and the current simulation time is determined to be the adjacent historical simulation time of the updated current simulation time. And returns to step S402.
In step S405, based on the position parameters of each type of fracturing fluid particles in each computation space at the adjacent historical simulation time of the termination simulation time, a simulation result of the viscous index phenomenon of each type of fracturing fluid particles in the fracturing fluid flowing region is generated.
It should be noted that, in the practical application scenario, the step S401 shown in fig. 4 is an alternative embodiment of the step S101 shown in fig. 1, the step S402 and the step S403 shown in fig. 4 are alternative embodiments of the step S102 shown in fig. 1, and the step S405 shown in fig. 4 is an alternative embodiment of the step S103 shown in fig. 1.
Corresponding to the above method embodiment, the present invention further provides a system for simulating the viscous fingering phenomenon, as shown in fig. 5, where the system for simulating the viscous fingering phenomenon includes:
the grid division module 501 is configured to grid-divide a model of a fracturing fluid flowing region in a target stratum to obtain a plurality of computation spaces;
The parameter calculation module 502 is configured to, in a case where the current simulation time is not the termination simulation time, calculate, for each calculation space: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation time according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the fracturing fluid parameters at the adjacent historical simulation time of the calculation space, wherein the adjacent historical simulation time is adjacent to the current simulation time, and the time value is earlier than one simulation time of the current simulation time;
The result generating module 503 is configured to generate, when the current simulation time is the termination simulation time, a simulation result of the viscous fingering phenomenon of each type of fracturing fluid particle in the fracturing fluid flowing region based on the position parameters of each type of fracturing fluid particle in each computation space at the adjacent historical simulation time of the termination simulation time.
Optionally, the parameter calculation module 502 is configured to:
Carrying out particle calibration processing on each calculation space, wherein the particle calibration processing comprises the following steps: calculating the position of each fracturing fluid particle under the current simulation moment according to the velocity field of each type of fracturing fluid particle of the calculation space under the adjacent historical simulation moment, deleting the fracturing fluid particle of which the position is not in the calculation space from the calculation space, and transmitting the type and the velocity field of the deleted fracturing fluid particle to the calculation space adjacent to the calculation space where the position is located;
For each calculation space after particle calibration treatment: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid.
Optionally, the parameter calculating module 502 calculates, according to the velocity field and the fracturing fluid parameters of each type of fracturing fluid particles at the adjacent historical simulation time, the position parameters and the velocity field of each type of fracturing fluid particles in the calculation space at the current simulation time by using the continuity equation of the preset fluid and the seepage darcy equation of the preset fluid, where the parameters and the velocity fields are set as follows:
For each fracturing fluid particle of the target type:
substituting the fluid density rho in the fluid parameters of the target type fracturing fluid and the velocity field of the fracturing fluid particles at the adjacent historical simulation time t-1 into a continuity equation of the preset fluid to obtain the continuity equation of the preset fluid substituted by the parameters:
，
Where t is the current simulation time instant, Is the partial derivative sign, b is the crack width of the computation space,/>Is the velocity of the fracturing fluid particles in the x-axis direction of the computation space in the velocity field,/>Is the velocity of the fracturing fluid particles in the z-axis direction of the computation space in the velocity field,/>For fluid loss in a target type of fracturing fluid parameter,/>The method comprises the steps of presetting pumping speed for fracturing fluid particles of a target type;
substituting the fluid density rho, the fluid pressure p and the fluid viscosity mu in the fluid parameters of the target type fracturing fluid into a seepage darcy equation of the preset fluid to obtain the seepage darcy equation of the preset fluid subjected to parameter substitution:
，
wherein, For the velocity of the target fracturing fluid particle in the x-axis direction of the calculation space at the current simulation time t to be solved,/>For the velocity of the target fracturing fluid particle in the z-axis direction of the calculation space at the current simulation time t to be solved,/>K is the crack permeability of the calculation space, g is the gravity acceleration;
simultaneously solving a continuity equation of the preset fluid subjected to parameter substitution and a seepage darcy equation of the preset fluid to obtain a model comprising And/>The velocity field of the fracturing fluid particles at the current simulation time is based on/>And/>And obtaining the position parameters of the fracturing fluid particles in the calculation space.
Optionally, the computation space includes a plurality of grids, and the result generating module 503 is configured to:
For each computation space: determining the distribution condition of fracturing fluid points in each grid of the calculation space according to the position parameters of the fracturing fluid points of each type; under the condition that the distribution condition characterizes that only one type of fracturing fluid particles exist in the grid, rendering the grid into a display state corresponding to the type; under the condition that the distribution condition represents that at least two types of fracturing fluid particles exist in the grid, interpolation averaging is carried out based on each type, and the grid is rendered into a display state corresponding to the interpolation averaging result; splicing the display states of the grids of the calculation space to obtain the display state of the calculation space;
And splicing the display states of the calculation spaces to generate simulation results of the viscous fingering phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing region.
The embodiment of the invention also provides a device for simulating the viscous finger-in phenomenon, as shown in fig. 6, which comprises:
A processor 601;
A memory 602 for storing instructions executable by the processor 601;
Wherein the processor 601 is configured to execute instructions to implement a method of simulating a sticky finger phenomenon as any one of the above.
The embodiment of the present invention also provides a computer-readable storage medium, which when executed by a processor of a device for simulating a viscous fingering phenomenon, causes the device for simulating a viscous fingering phenomenon to perform a method for simulating a viscous fingering phenomenon as any one of the above.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In one typical configuration, the device includes one or more processors (CPUs), memory, and a bus. The device may also include input/output interfaces, network interfaces, and the like.
The memory may include volatile memory, random Access Memory (RAM), and/or nonvolatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM), among other forms in computer readable media, the memory including at least one memory chip. Memory is an example of a computer-readable medium.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises an element.
In this specification, each embodiment is described in a related manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for system embodiments, since they are substantially similar to method embodiments, the description is relatively simple, as relevant to see a section of the description of method embodiments.
The foregoing is merely exemplary of the present application and is not intended to limit the present application. Various modifications and variations of the present application will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. which come within the spirit and principles of the application are to be included in the scope of the claims of the present application.
Claims (8)
1. A method of simulating the phenomenon of viscous finger-stick, the method comprising:
Performing grid division on a model of a fracturing fluid flowing region in a target stratum to obtain a plurality of calculation spaces;
In the case where the current simulation time is not the termination simulation time, for each calculation space: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation time according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the fracturing fluid parameters of the adjacent historical simulation time of the calculation space, wherein the adjacent historical simulation time is adjacent to the current simulation time, the time value is earlier than one simulation time of the current simulation time, and the parameters of the various types of fracturing fluid comprise quality, density and viscosity;
Under the condition that the current simulation time is the termination simulation time, generating a simulation result of the viscous finger-in phenomenon of each type of fracturing fluid particle in the fracturing fluid flowing area based on the position parameters of each type of fracturing fluid particle of each calculation space at the adjacent historical simulation time of the termination simulation time;
Wherein, in the case that the current simulation time is not the termination simulation time, for each calculation space: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the adjacent historical simulation moment of the calculation space and the fracturing fluid parameters, wherein the method comprises the following steps:
Performing particle calibration processing on each computing space, wherein the particle calibration processing comprises: calculating the position of each fracturing fluid particle at the current simulation moment according to the velocity field of each type of fracturing fluid particle of the calculation space at the adjacent historical simulation moment, deleting the fracturing fluid particle of which the position is not in the calculation space from the calculation space, and transmitting the type of the deleted fracturing fluid particle and the velocity field to the calculation space adjacent to the calculation space where the position is based on a multipoint interface protocol;
For each of the computation spaces after the particle calibration process: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid.
2. The method of claim 1, wherein calculating the location parameters and the velocity fields of each type of the fracturing fluid particles in the calculation space at the current simulation time from the velocity fields and the fracturing fluid parameters of each type of the fracturing fluid particles at the adjacent historical simulation time using a continuity equation of a preset fluid and a seepage darcy equation of the preset fluid comprises:
For each of the fracturing fluid particles of the target type:
Substituting the fluid density rho in the fluid parameters of the fracturing fluid of the target type and the velocity field of the fracturing fluid particles at the adjacent historical simulation time t-1 into the continuity equation of the preset fluid to obtain the continuity equation of the preset fluid subjected to parameter substitution:
wherein t is the current simulation time, the Taking the sign of partial derivative, b as the crack width of the calculation space, u x(t-1) as the speed of the fracturing fluid particles in the x-axis direction of the calculation space in the speed field, u z(t-1) as the speed of the fracturing fluid particles in the z-axis direction of the calculation space in the speed field, q l as the fluid loss in the fracturing fluid parameters of the target type, and q in as the preset pumping speed of the fracturing fluid particles of the target type;
substituting the fluid density ρ, the fluid pressure p and the fluid viscosity μ in the fluid parameters of the fracturing fluid of the target type into a seepage darcy equation of the preset fluid to obtain a seepage darcy equation of the preset fluid subjected to parameter substitution:
Wherein, v xt is the velocity of the target fracturing fluid particle in the x-axis direction of the calculation space at the current simulation time t to be calculated, v zt is the velocity of the target fracturing fluid particle in the z-axis direction of the calculation space at the current simulation time t to be calculated, and The K is the crack permeability of the calculation space, and the g is the gravity acceleration;
And simultaneously solving a continuity equation of the preset fluid and a seepage darcy equation of the preset fluid, which are substituted by the parameters, to obtain a velocity field of the fracturing fluid particles including v xt and v zt at the current simulation moment, and obtaining the position parameters of the fracturing fluid particles in the calculation space based on v xt and v zt.
3. The method of claim 1, wherein the computation space comprises a plurality of grids, and wherein the generating a simulation result of the viscous fingering of each type of fracturing fluid particle in the fracturing fluid flow-through region based on the location parameters of each type of fracturing fluid particle for each computation space at adjacent historical simulation times of the termination simulation time, if the current simulation time is a termination simulation time, comprises:
For each computation space: determining the distribution condition of the fracturing fluid points in each grid of the calculation space according to the position parameters of the fracturing fluid points of each type; rendering the grid into a display state corresponding to the type under the condition that the distribution condition characterizes that only one type of fracturing fluid particles exist in the grid; under the condition that the distribution condition represents that at least two types of fracturing fluid particles exist in the grid, interpolation average is carried out based on each type, and the grid is rendered into a display state corresponding to the interpolation average result; splicing the display states of the grids of the calculation space to obtain the display state of the calculation space;
And splicing the display states of the calculation spaces to generate simulation results of the viscous fingering phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing area.
4. A system for simulating the phenomenon of viscous finger-in, the system comprising:
the grid division module is used for carrying out grid division on a model of a fracturing fluid flowing region in the target stratum to obtain a plurality of calculation spaces;
The parameter calculation module is used for calculating the space for each calculation when the current simulation time is not the termination simulation time: calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation time according to the velocity fields of the various types of fracturing fluid particles in the calculation space and the fracturing fluid parameters of the adjacent historical simulation time of the calculation space, wherein the adjacent historical simulation time is adjacent to the current simulation time, the time value is earlier than one simulation time of the current simulation time, and the parameters of the various types of fracturing fluid comprise quality, density and viscosity;
The result generation module is used for generating simulation results of the viscous finger-in phenomenon of each type of fracturing fluid particles in the fracturing fluid flowing area based on the position parameters of each type of fracturing fluid particles of each calculation space at the adjacent historical simulation time of the termination simulation time under the condition that the current simulation time is the termination simulation time;
Wherein the parameter calculation module is configured to:
Performing particle calibration processing on each computing space, wherein the particle calibration processing comprises: calculating the position of each fracturing fluid particle at the current simulation moment according to the velocity field of each type of fracturing fluid particle of the calculation space at the adjacent historical simulation moment, deleting the fracturing fluid particle of which the position is not in the calculation space from the calculation space, and transmitting the type of the deleted fracturing fluid particle and the velocity field to the calculation space adjacent to the calculation space where the position is based on a multipoint interface protocol;
For each of the computation spaces after the particle calibration process: and calculating the position parameters and the velocity fields of the various types of fracturing fluid particles in the calculation space at the current simulation moment according to the velocity fields and the fracturing fluid parameters of the various types of fracturing fluid particles at the adjacent historical simulation moment by using a continuity equation of the preset fluid and a seepage Darcy equation of the preset fluid.
5. The system of claim 4, wherein the parameter calculation module, when using the continuity equation of the preset fluid and the percolation darcy equation of the preset fluid, calculates the position parameters and velocity fields of each type of fracturing fluid particle in the calculation space at the current simulation time from the velocity fields and the fracturing fluid parameters of each type of fracturing fluid particle at the adjacent historical simulation time, is configured to:
For each of the fracturing fluid particles of the target type:
Substituting the fluid density rho in the fluid parameters of the fracturing fluid of the target type and the velocity field of the fracturing fluid particles at the adjacent historical simulation time t-1 into the continuity equation of the preset fluid to obtain the continuity equation of the preset fluid subjected to parameter substitution:
wherein t is the current simulation time, the Taking the sign of partial derivative, b as the crack width of the calculation space, u x(t-1) as the speed of the fracturing fluid particles in the x-axis direction of the calculation space in the speed field, u z(t-1) as the speed of the fracturing fluid particles in the z-axis direction of the calculation space in the speed field, q l as the fluid loss in the fracturing fluid parameters of the target type, and q in as the preset pumping speed of the fracturing fluid particles of the target type;
substituting the fluid density ρ, the fluid pressure p and the fluid viscosity μ in the fluid parameters of the fracturing fluid of the target type into a seepage darcy equation of the preset fluid to obtain a seepage darcy equation of the preset fluid subjected to parameter substitution:
Wherein, v xt is the velocity of the target fracturing fluid particle in the x-axis direction of the calculation space at the current simulation time t to be calculated, v zt is the velocity of the target fracturing fluid particle in the z-axis direction of the calculation space at the current simulation time t to be calculated, and The K is the crack permeability of the calculation space, and the g is the gravity acceleration;
And simultaneously solving a continuity equation of the preset fluid and a seepage darcy equation of the preset fluid, which are substituted by the parameters, to obtain a velocity field of the fracturing fluid particles including v xt and v zt at the current simulation moment, and obtaining the position parameters of the fracturing fluid particles in the calculation space based on v xt and v zt.
6. The system of claim 4, wherein the computation space comprises a plurality of grids, the result generation module configured to:
For each computation space: determining the distribution condition of the fracturing fluid points in each grid of the calculation space according to the position parameters of the fracturing fluid points of each type; rendering the grid into a display state corresponding to the type under the condition that the distribution condition characterizes that only one type of fracturing fluid particles exist in the grid; under the condition that the distribution condition represents that at least two types of fracturing fluid particles exist in the grid, interpolation average is carried out based on each type, and the grid is rendered into a display state corresponding to the interpolation average result; splicing the display states of the grids of the calculation space to obtain the display state of the calculation space;
And splicing the display states of the calculation spaces to generate simulation results of the viscous fingering phenomenon of the various types of fracturing fluid particles in the fracturing fluid flowing area.
7. A device for simulating the phenomenon of viscous finger-in, said device comprising:
A processor;
A memory for storing the processor-executable instructions;
wherein the processor is configured to execute the instructions to implement a method of simulating the viscous finger phenomenon of any one of claims 1 to 3.
8. A computer readable storage medium, characterized in that instructions in the computer readable storage medium, when executed by a processor of a simulation device of the viscous fingering phenomenon, enable the simulation device of the viscous fingering phenomenon to perform the simulation method of the viscous fingering phenomenon according to any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410259203.6A CN117854605B (en) | 2024-03-07 | 2024-03-07 | Method, system, equipment and storage medium for simulating viscous finger-in phenomenon |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410259203.6A CN117854605B (en) | 2024-03-07 | 2024-03-07 | Method, system, equipment and storage medium for simulating viscous finger-in phenomenon |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117854605A CN117854605A (en) | 2024-04-09 |
| CN117854605B true CN117854605B (en) | 2024-05-24 |
Family
ID=90540467
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410259203.6A Active CN117854605B (en) | 2024-03-07 | 2024-03-07 | Method, system, equipment and storage medium for simulating viscous finger-in phenomenon |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117854605B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108427605A (en) * | 2018-02-09 | 2018-08-21 | 中国地质大学（北京） | The accelerated method of streamline simulation is realized based on particles trace algorithm |
| CN117113488A (en) * | 2023-08-23 | 2023-11-24 | 北京源烃泰克科技有限公司 | Fracturing simulator implementation method based on object point method |
| CN117307122A (en) * | 2022-06-20 | 2023-12-29 | 中国石油天然气股份有限公司 | Sequential mixed injection acid pressure injection process parameter optimization method based on acid liquor finger-feed experiment |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8788250B2 (en) * | 2007-05-24 | 2014-07-22 | Exxonmobil Upstream Research Company | Method of improved reservoir simulation of fingering systems |
| US9576084B2 (en) * | 2013-08-27 | 2017-02-21 | Halliburton Energy Services, Inc. | Generating a smooth grid for simulating fluid flow in a well system environment |
-
2024
- 2024-03-07 CN CN202410259203.6A patent/CN117854605B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108427605A (en) * | 2018-02-09 | 2018-08-21 | 中国地质大学（北京） | The accelerated method of streamline simulation is realized based on particles trace algorithm |
| CN117307122A (en) * | 2022-06-20 | 2023-12-29 | 中国石油天然气股份有限公司 | Sequential mixed injection acid pressure injection process parameter optimization method based on acid liquor finger-feed experiment |
| CN117113488A (en) * | 2023-08-23 | 2023-11-24 | 北京源烃泰克科技有限公司 | Fracturing simulator implementation method based on object point method |
Non-Patent Citations (1)
| Title |
|---|
| 周维垣等.《岩石力学数值计算方法》.中国电力出版社,2005,第451-456页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117854605A (en) | 2024-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pollock | User guide for MODPATH Version 7—A particle-tracking model for MODFLOW | |
| Tielen et al. | A high order material point method | |
| George | Adaptive finite volume methods with well‐balanced Riemann solvers for modeling floods in rugged terrain: Application to the Malpasset dam‐break flood (France, 1959) | |
| US11941331B2 (en) | Multi-phase flow visualizations based on fluid occupation time | |
| MX2014008897A (en) | Systems and methods for estimating fluid breakthrough times at producing well locations. | |
| CN107817525B (en) | A kind of near surface constant speed fill method and device | |
| CN111175818A (en) | CO2Method for judging position of gas drive front edge and method and device for training model of gas drive front edge | |
| Settgast et al. | Simulation of hydraulic fracture networks in three dimensions utilizing massively parallel computing resources | |
| CA2919871C (en) | A geostatistical procedure for simulation of the 3d geometry of a natural fracture network conditioned by well bore observations | |
| Zago et al. | Preliminary validation of lava benchmark tests on the GPUSPH particle engine | |
| CN104424294A (en) | Information processing method and information processing device | |
| CN117854605B (en) | Method, system, equipment and storage medium for simulating viscous finger-in phenomenon | |
| AU2013398344B2 (en) | Local updating of 3D geocellular model | |
| CN108595782B (en) | Calculation method for mass transfer between matrix and cracks in discrete cracks | |
| CN116882221B (en) | Geothermal exploitation numerical simulation method and system based on three-dimensional fracture type thermal storage model | |
| WO2024051131A1 (en) | Water environmental capacity calculation method and apparatus, electronic device, and storage medium | |
| CN113160399B (en) | Method, system and readable storage medium for constructing three-dimensional geological model | |
| CN116016225A (en) | Information processing method, signal positioning method, device, equipment and medium | |
| CN112419493B (en) | Shale reservoir three-dimensional attribute model building method and device | |
| Stosic et al. | Computational approach to improving the efficiency of river discharge measurement | |
| CN107256004A (en) | Product processing method and device | |
| CN113378445B (en) | Gas-liquid multiphase system calculation method and system based on discrete simulation | |
| CN109215134A (en) | Threedimensional model blocks determining method, apparatus, computer equipment and storage medium | |
| CN115168975B (en) | BIM technology-based pavement evenness quality control method and device | |
| CN117973236B (en) | Correction method, system, equipment and storage medium for simulation result |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |