CN109033672B - Dynamic crack determination method and device for displacement simulation mesoporous throat network model - Google Patents
Dynamic crack determination method and device for displacement simulation mesoporous throat network model Download PDFInfo
- Publication number
- CN109033672B CN109033672B CN201810901576.3A CN201810901576A CN109033672B CN 109033672 B CN109033672 B CN 109033672B CN 201810901576 A CN201810901576 A CN 201810901576A CN 109033672 B CN109033672 B CN 109033672B
- Authority
- CN
- China
- Prior art keywords
- throat
- network model
- dimensional pore
- pore
- unit
- 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
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 117
- 238000004088 simulation Methods 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 77
- 239000011148 porous material Substances 0.000 claims abstract description 170
- 230000008569 process Effects 0.000 claims abstract description 38
- 238000005315 distribution function Methods 0.000 claims description 22
- 238000009736 wetting Methods 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 239000012530 fluid Substances 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 3
- 238000005325 percolation Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 15
- 238000011161 development Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003938 response to stress Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the application provides a dynamic crack judgment method and a dynamic crack judgment device for a displacement simulation mesoporous and laryngeal network model, wherein the method comprises the following steps: acquiring characteristic parameters of the porous medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters; displacement simulation of the seepage rule of the two-dimensional pore throat network model; in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model. The method and the device can realize dynamic crack judgment of the displacement simulation mesoporous throat network model and improve the accuracy of displacement simulation.
Description
Technical Field
The application relates to the technical field of development and research of unconventional oil and gas reservoirs, in particular to a dynamic crack judgment method and device for a displacement simulation mesoporous and laryngeal network model.
Background
The development and research of unconventional oil and gas reservoirs have strategic significance on energy in China, the reservoirs have the characteristics of dual porous media, and the research on the influence of micropore distribution and dynamic micropore opening and expansion on seepage is the key for solving the development of the resources. For the development problem of unconventional oil and gas reservoirs with complex pore structure characteristics, the problem is difficult to solve by a macroscopic seepage theory, so that research on microstructure characteristics and a microscopic seepage mechanism is a new breakthrough.
At present, in unconventional oil and gas reservoir development simulation based on microstructure characteristics and a micro seepage mechanism, the used pore throat network models mainly comprise an equivalent continuum model, a discrete network model, a mixed model, a percolation model and a pore throat network model. The pore throat network model based on the percolation theory has shown certain advantages on the understanding of the conventional oil and gas reservoir development mechanism. However, since reservoir development is a dynamic process, and the pore throat network models do not consider the process of dynamic throat cracking, the displacement simulation process has certain errors from the actual displacement process. Therefore, a dynamic pore throat network model needs to be established, and how to dynamically judge the fracture in the displacement simulation in the establishment of the dynamic pore throat network model is a technical problem to be solved at present.
Disclosure of Invention
The embodiment of the application aims to provide a dynamic fracture judgment method and a dynamic fracture judgment device for a displacement simulation mesoporous throat network model, so as to realize dynamic fracture judgment on the displacement simulation mesoporous throat network model.
In order to achieve the above object, in one aspect, an embodiment of the present application provides a dynamic fracture determination method for a displacement simulation mesopore throat network model, including:
acquiring characteristic parameters of the porous medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters;
displacement simulation of the seepage rule of the two-dimensional pore throat network model;
in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model.
The dynamic fracture determination method for the pore throat network model in the displacement simulation of the embodiment of the application further comprises the following steps before the displacement simulation of the seepage rule of the two-dimensional pore throat network model:
determining that the wetting phase of the displacement simulation is water and the non-wetting phase is oil;
assigning the upper and lower boundary states of each grid in the two-dimensional pore throat network model as fixed values, and assigning the left and right boundary states as variable values; the fixed value represents that the state of the boundary is fixed during the displacement simulation, and the variable value represents that the state of the boundary can be changed during the displacement simulation;
and assigning the states of all pore units and throat units of the two-dimensional pore throat network model as 0 values, wherein the 0 values represent the closed state and are full of water.
The dynamic fracture determination method for the displacement simulation mesoporous aperture throat network model in the embodiment of the application, before confirming that the peak pressure of the throat unit in the two-dimensional aperture throat network model reaches the critical pressure value, further includes:
when the seepage rule is simulated, assigning the states of the pore units and the throat units of the displacement inlet boundary in the two-dimensional pore throat network model to be 1 value, wherein the 1 value represents an open state and is full of oil;
confirming whether a throat unit with a state of 0 exists in throat units connected with the pore unit with a state of 1;
and if so, judging whether the peak pressure of the throat unit reaches a critical pressure value.
The dynamic fracture determination method for the displacement simulation mesopore throat network model according to the embodiment of the application, when it is determined that the peak pressure of the throat unit in the two-dimensional pore throat network model reaches the critical pressure value, the corresponding throat unit is opened to form a fracture, and the method includes the following steps:
when it is confirmed that the peak pressure of the throat unit reaches the critical pressure value, the state of the throat unit is changed to 1 value, and the state of the pore unit connected to the throat unit is changed to 1 value.
The dynamic crack determination method for the displacement simulation mesopore throat network model in the embodiment of the application further comprises the following steps:
and when the vertex pressure of the throat unit is confirmed to not reach the critical pressure value, maintaining the current state of the throat unit.
The method for judging the dynamic cracks of the pore throat network model in the displacement simulation, which is provided by the embodiment of the application, comprises the following steps of:
determining a statistical distribution function for representing the distribution rule of the characteristic parameters;
and generating a pore throat network model according to the statistical distribution function, and assigning the characteristic parameters to throat units and pore units in the pore throat network model so as to form a two-dimensional pore throat network model.
The method for judging the dynamic fracture of the displacement simulation mesopore throat network model in the embodiment of the application, wherein the step of determining the statistical distribution function for representing the distribution rule of the characteristic parameters comprises the following steps:
based on the formula
Representing the distribution rule of pore radius;
based on the formula
Representing the distribution rule of throat radius;
wherein p (R) is a distribution function of pore radii, f (x) is a distribution function of throat radii, R is a pore radiusmaxMaximum value of pore radius, RminIs the maximum value of the pore radius, σ is the standard deviation of the distribution function, e is the natural constant, μ is the expected value of the distribution function, and x is the throat radius.
The dynamic fracture judgment method for the pore throat network model in the displacement simulation of the embodiment of the application, wherein the displacement simulation of the seepage rule of the two-dimensional pore throat network model, comprises the following steps:
according to the formula
Simulating the seepage rule of a throat unit in the two-dimensional pore throat network model;
according to the formula
Simulating the seepage rule of the crack in the two-dimensional pore throat network model;
wherein q isijIs the flow at the ith row and jth column position in the two-dimensional pore throat network model, rijIs the throat radius at the ith row and jth column position in the two-dimensional pore throat network model,
is the dynamic viscosity coefficient of the fluid, LijIs the fracture length from the ith row and jth column position in the two-dimensional pore throat network model, PjIs the pressure at the jth node, P, in a two-dimensional pore throat network modeliIs the pressure at the ith node in the two-dimensional pore throat network model,
theta is the contact angle, sigmawnIs the interfacial tension between a wetting phase and a non-wetting phase, r is the capillary radius, q is the flow in a pore unit, e is a natural constant, mu is the dynamic viscosity coefficient of fluid in the pore, and delta P is Pj-Pi-Pc。
On the other hand, the embodiment of the present application provides a dynamic fracture determination device for a displacement simulation mesopore throat network model, including:
the model construction module is used for acquiring the characteristic parameters of the porous medium and constructing a two-dimensional pore throat network model according to the characteristic parameters;
the seepage simulation module is used for displacement simulation of the seepage rule of the two-dimensional pore throat network model;
and the crack judging module is used for opening the corresponding throat unit to form a crack when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach the critical pressure value in the process of simulating the seepage rule so as to update the two-dimensional pore throat network model.
In another aspect, an embodiment of the present application provides another dynamic fracture determination apparatus for a displacement simulation mesoporous and laryngeal network model, including a memory, a processor, and a computer program stored on the memory, where the computer program is executed by the processor to perform the following steps:
acquiring characteristic parameters of the porous medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters;
displacement simulation of the seepage rule of the two-dimensional pore throat network model;
in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model.
The method comprises the steps of firstly obtaining characteristic parameters of a pore throat medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters; secondly, displacement simulation is carried out on the seepage rule of the two-dimensional pore throat network model; and then in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach the critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model, thereby realizing the dynamic crack judgment of the displacement simulation middle pore throat network model. And because the pore throat network model of the embodiment of the application considers that pores of the unconventional oil and gas reservoir are dynamically cracked in the displacement process, the displacement process simulation of the embodiment of the application is closer to the actual displacement process of the unconventional oil and gas reservoir, so that the errors of the displacement simulation process and the actual displacement process of the unconventional oil and gas reservoir are reduced, and the displacement process simulation closer to the actual displacement process of the unconventional oil and gas reservoir can provide a more objective and accurate reference basis for the subsequent displacement development scheme of the unconventional oil and gas reservoir.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a flow chart of a dynamic fracture determination method for a displacement simulation mesoporous and laryngeal network model according to an embodiment of the present application;
FIG. 2 is a plot of pore throat radius distribution obtained in an example of the present application;
FIG. 3 is a schematic diagram of fitting data based on matlab in one embodiment of the present application;
FIGS. 4a to 4d are schematic views of a displacement path with the number of cracks as a main control parameter in an embodiment of the present application;
FIGS. 5 a-5 d are schematic views of a displacement path with a fracture length as a master parameter in an embodiment of the present application;
fig. 6a to 6d are schematic views of a displacement path with a displacement pressure difference as a main control parameter in an embodiment of the present application;
FIG. 7 is a block diagram of a dynamic fracture determination apparatus for a displacement simulation mesoporous and laryngeal network model according to an embodiment of the present disclosure;
fig. 8 is a block diagram of a dynamic fracture determination device of a displacement simulation mesopore-throat network model according to another embodiment of the present application.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application. For example, in the following description, forming the second component over the first component may include embodiments in which the first and second components are formed in direct contact, embodiments in which the first and second components are formed in non-direct contact (i.e., additional components may be included between the first and second components), and so on.
Also, for ease of description, some embodiments of the present application may use spatially relative terms such as "above …," "below …," "top," "below," etc., to describe the relationship of one element or component to another (or other) element or component as illustrated in the various figures of the embodiments. It will be understood that the spatially relative terms are intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements or components described as "below" or "beneath" other elements or components would then be oriented "above" or "over" the other elements or components.
Referring to fig. 1, a dynamic fracture determination method for a displacement simulation mesopore throat network model according to an embodiment of the present application may include the following steps:
s101, obtaining characteristic parameters of a pore throat medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters.
In an embodiment of the application, the characteristic parameter of the pore throat medium is obtained by obtaining an actual core of an unconventional oil and gas reservoir. Wherein the characteristic parameters may include, for example, the distribution of pores and throats, coordination number and spatial correlation, etc. (e.g., as shown in fig. 2). In an exemplary embodiment of the present application, the characteristic parameters of the real core may be obtained by using nuclear magnetic equipment or the like.
In an embodiment of the present application, the constructing a two-dimensional pore throat network model according to the characteristic parameters may include the following steps:
1) determining a statistical distribution function for representing the distribution rule of the characteristic parameters; generally, any kind of characteristic parameter can select a distribution function matching with the distribution rule according to the distribution rule, so that the distribution rule of the characteristic parameter can be represented by the distribution function. In another embodiment of the present application, software such as matlab and the like may also be used to fit the characteristic parameters (for example, as shown in fig. 3), so that the distribution rule of the characteristic parameters may be represented by the fitted curve.
In an exemplary embodiment of the present application, the data may be based on a formula, for example
Representing the distribution rule of pore radius; in an exemplary embodiment of the present application, the data may be based on a formula, for example
Representing the distribution rule of the throat radius.
Wherein p (R) is a distribution function of pore radii, f (x) is a distribution function of throat radii, R is a pore radiusmaxMaximum value of pore radius, RminIs the maximum value of the pore radius, σ is the standard deviation of the distribution function, e is the natural constant, μ is the expected value of the distribution function, and x is the throat radius. 2) And generating a pore throat network model according to the statistical distribution function, and assigning the characteristic parameters to throat units and pore units in the pore throat network model so as to form a two-dimensional pore throat network model. In another embodiment of the present application, a pore throat network model may also be generated for a fitted curve used to characterize the distribution rule of the characteristic parameters, and the characteristic parameters are assigned to throat units and pore units in the pore throat network model, thereby forming a two-dimensional pore throat network model.
S102, displacement simulation of the seepage rule of the two-dimensional pore throat network model.
In an embodiment of the present application, the seepage law may also be referred to as a fluid seepage mechanism, which is a seepage mechanism including dynamic fractures. The dynamic fracture is characterized in that a closed throat is continuously fractured under the action of displacement pressure in the displacement simulation process, so that fluid flows in a solid fracture with obvious stress response.
In an exemplary embodiment of the present application, the pore unit may be represented by a small sphere, and the throat unit may be represented by a cylindrical pipe, thereby constituting a pore-throat basic unit. Wherein the throat unit can adopt a circular section unit, and the flow rule can be according to the formula
Characterizing; the flow rule of the crack unit can be simplified into flat plate flow and can be according to a formula
And (5) characterizing.
Wherein q isijIs the flow at the ith row and jth column position in the two-dimensional pore throat network model, rijIs the throat radius at the ith row and jth column position in the two-dimensional pore throat network model,
is the dynamic viscosity coefficient of the fluid, LijIs the fracture length from the ith row and jth column position in the two-dimensional pore throat network model, PjIs the pressure at the jth node, P, in a two-dimensional pore throat network modeliIs the pressure at the ith node, P, in a two-dimensional pore throat network modelcIs capillary force, and
theta is the contact angle, sigmawnIs the interfacial tension between a wetting phase and a non-wetting phase, r is the capillary radius, q is the flow in a pore unit, e is a natural constant, mu is the dynamic viscosity coefficient of fluid in the pore, and delta P is Pj-Pi-PcWhen the wetting phase displaces the non-wetting phase, Δ P takes a positive sign, and when the non-wetting phase displaces the wetting phase, Δ P takes a negative sign.
In another exemplary embodiment of the present application, the above formula
For two-phase multi-interface flow. When it is desired to characterize two-phase single-interface flow, the formula can be simplified as follows:
further, when it is desired to characterize single-phase flow, the formula can be further simplified as follows:
in an embodiment of the present application, the displacement master control parameters may be selected as needed during the displacement simulation. In an exemplary embodiment of the present application, fig. 4a to 4d show schematic diagrams of a displacement path with the number of cracks as a main control parameter. Wherein fig. 4a is a schematic view of a displacement path without a fracture (i.e., a crack); FIG. 4b is a schematic view of a displacement path comprising 5 fractures; FIG. 4c is a schematic of a displacement path comprising 10 fractures; fig. 4d is a schematic diagram of a displacement path comprising 40 fractures. In another exemplary embodiment of the present application, fig. 5 a-5 d show schematic views of a displacement path with a fracture length as a master parameter. Wherein, FIG. 5a is a schematic view of a crack-free displacement path; FIG. 5b is a schematic of a displacement path comprising 15 fractures of 600 microns in length; FIG. 5c is a schematic of a displacement path comprising 15 fractures of 1500 microns in length; figure 5d is a schematic of a displacement path comprising 15 fractures of 3000 microns in length. In another exemplary embodiment of the present application, fig. 6 a-6 d show displacement path schematics with displacement differential pressure as a master parameter. Wherein, FIG. 6a is a schematic diagram of a displacement path with a displacement differential pressure of 1 MPa; FIG. 6b is a schematic of the displacement path at a displacement differential pressure of 3 MPa; FIG. 6c is a schematic of the displacement path at a displacement differential pressure of 6 MPa; fig. 6d is a schematic diagram of the displacement path at a displacement differential pressure of 9 MPa.
S103, in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model. A split is herein understood to be a split throat unit.
In one embodiment of the application, in order to dynamically simulate the process that the throat of the reservoir is fractured under the action of displacement pressure, the critical pressure value of the top pressure of the throat unit can be preset; and for each throat unit in the model, when the peak pressure of the throat unit reaches a critical pressure value, opening the corresponding throat unit to form a crack, and updating the two-dimensional pore throat network model.
In an embodiment of the present application, before the displacement simulation of the seepage law of the two-dimensional pore throat network model, the method may further include the following steps:
determining that the wetting phase of the displacement simulation is water and the non-wetting phase is oil;
assigning the upper and lower boundary states of each grid in the two-dimensional pore throat network model as fixed values, and assigning the left and right boundary states as variable values; the fixed value represents that the state of the boundary is fixed during the displacement simulation, and the variable value represents that the state of the boundary can be changed during the displacement simulation;
and assigning the states of all pore units and throat units of the two-dimensional pore throat network model as 0 values, wherein the 0 values represent the closed state and are full of water.
Correspondingly, before the step of confirming that the peak pressure of the throat unit in the two-dimensional pore throat network model reaches the critical pressure value, the method also comprises the following steps:
when the seepage rule is simulated, assigning the states of the pore units and the throat units of the displacement inlet boundary in the two-dimensional pore throat network model to be 1 value, wherein the 1 value represents an open state and is full of oil;
confirming whether a throat unit with a state of 0 exists in throat units connected with the pore unit with a state of 1;
and if so, judging whether the peak pressure of the throat unit reaches a critical pressure value. And when the vertex pressure of the throat unit is confirmed to not reach the critical pressure value, maintaining the current state of the throat unit.
Correspondingly, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach the critical pressure value, the corresponding throat unit is opened to form a crack, and the method specifically comprises the following steps:
when it is confirmed that the peak pressure of the throat unit reaches the critical pressure value, the state of the throat unit is changed to 1 value, and the state of the pore unit connected to the throat unit is changed to 1 value.
Therefore, the characteristic parameters of the pore throat medium are obtained firstly, and a two-dimensional pore throat network model is constructed according to the characteristic parameters; secondly, displacement simulation is carried out on the seepage rule of the two-dimensional pore throat network model; and then in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach the critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model, thereby realizing the dynamic crack judgment of the displacement simulation middle pore throat network model. And because the pore throat network model of the embodiment of the application considers that pores of the unconventional oil and gas reservoir are dynamically cracked in the displacement process, the displacement process simulation of the embodiment of the application is closer to the actual displacement process of the unconventional oil and gas reservoir, so that the errors of the displacement simulation process and the actual displacement process of the unconventional oil and gas reservoir are reduced, and the displacement process simulation closer to the actual displacement process of the unconventional oil and gas reservoir can provide a more objective and accurate reference basis for the subsequent displacement development scheme of the unconventional oil and gas reservoir.
Referring to fig. 7, a dynamic fracture determination device for a displacement simulation mesopore-throat network model according to an embodiment of the present application may include:
the model construction module 71 may be configured to obtain characteristic parameters of a pore throat medium, and construct a two-dimensional pore throat network model according to the characteristic parameters;
the seepage simulation module 72 can be used for displacement simulation of the seepage rule of the two-dimensional pore throat network model;
the crack determination module 73 may be configured to, in the process of simulating the seepage rule, open the corresponding throat unit to form a crack when it is determined that the peak pressure of the throat unit in the two-dimensional pore throat network model reaches a critical pressure value, so as to update the two-dimensional pore throat network model;
referring to fig. 8, another dynamic fracture determination apparatus for a displacement simulation mesopore throat network model according to an embodiment of the present application includes a memory, a processor, and a computer program stored in the memory, where the computer program is executed by the processor to perform the following steps:
acquiring characteristic parameters of a pore throat medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters;
displacement simulation of the seepage rule of the two-dimensional pore throat network model;
in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams 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.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). 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 computer storage media 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 that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
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, 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, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method or apparatus that comprises the element.
As will be appreciated by one skilled in the art, 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.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (9)
1. A dynamic fracture judgment method for a displacement simulation mesopore throat network model is characterized by comprising the following steps:
acquiring characteristic parameters of the porous medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters;
displacement simulation of the seepage rule of the two-dimensional pore throat network model;
determining that the wetting phase of the displacement simulation is water and the non-wetting phase is oil;
assigning the upper and lower boundary states of each grid in the two-dimensional pore throat network model as fixed values, and assigning the left and right boundary states as variable values; the fixed value represents that the state of the boundary is fixed during the displacement simulation, and the variable value represents that the state of the boundary can be changed during the displacement simulation;
assigning the states of all pore units and throat units of the two-dimensional pore throat network model to be 0 values, wherein the 0 values represent closed states and are full of water;
in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model.
2. The method for determining dynamic fractures of a displacement simulation mesoporous network model according to claim 1, wherein before confirming that the peak pressure of the throat unit in the two-dimensional pore throat network model reaches the critical pressure value, the method further comprises:
when the seepage rule is simulated, assigning the states of the pore units and the throat units of the displacement inlet boundary in the two-dimensional pore throat network model to be 1 value, wherein the 1 value represents an open state and is full of oil;
confirming whether a throat unit with a state of 0 exists in throat units connected with the pore unit with a state of 1;
and if so, judging whether the peak pressure of the throat unit reaches a critical pressure value.
3. The method for determining dynamic fractures of a displacement simulation mesopore throat network model according to claim 2, wherein when it is confirmed that the peak pressure of the throat unit in the two-dimensional pore throat network model reaches the critical pressure value, the corresponding throat unit is opened to form the fractures, the method comprises the following steps:
when it is confirmed that the peak pressure of the throat unit reaches the critical pressure value, the state of the throat unit is changed to 1 value, and the state of the pore unit connected to the throat unit is changed to 1 value.
4. The method for determining dynamic fractures of a displacement simulation mesoporous throat network model according to claim 2, further comprising:
and when the vertex pressure of the throat unit is confirmed to not reach the critical pressure value, maintaining the current state of the throat unit.
5. The method for determining dynamic fractures of a pore throat network model in a displacement simulation according to claim 1, wherein the constructing a two-dimensional pore throat network model according to the characteristic parameters comprises:
determining a statistical distribution function for representing the distribution rule of the characteristic parameters;
and generating a pore throat network model according to the statistical distribution function, and assigning the characteristic parameters to throat units and pore units in the pore throat network model so as to form a two-dimensional pore throat network model.
6. The method for determining dynamic fractures of a displacement simulation mesopore throat network model according to claim 5, wherein the determining a statistical distribution function for characterizing the distribution rule of the characteristic parameters comprises:
based on the formula
Representing the distribution rule of pore radius;
based on the formula
Representing the distribution rule of throat radius;
wherein p (R) is a distribution function of pore radii, f (x) is a distribution function of throat radii, R is a pore radiusmaxMaximum value of pore radius, RminIs the minimum value of the pore radius, σ is the standard deviation of the distribution function, e is the natural constant, μ is the expected value of the distribution function, and x is the throat radius.
7. The method for determining dynamic fractures of a pore throat network model in a displacement simulation of claim 1, wherein the displacement simulation of the percolation law of the two-dimensional pore throat network model comprises:
according to the formula
Simulating the seepage rule of a throat unit in the two-dimensional pore throat network model;
according to the formula
Simulating the seepage rule of the crack in the two-dimensional pore throat network model;
wherein q isijIs the flow at the ith row and jth column position in the two-dimensional pore throat network model, rijIs the throat radius at the ith row and jth column position in the two-dimensional pore throat network model,
is the dynamic viscosity coefficient of the fluid, LijIs the fracture length from the ith row and jth column position in the two-dimensional pore throat network model, PjIs the pressure at the jth node, P, in a two-dimensional pore throat network modeliIs the pressure at the ith node in the two-dimensional pore throat network model,
theta is the contact angle, sigmawnIs the interfacial tension between a wetting phase and a non-wetting phase, r is the capillary radius, q is the flow in a pore unit, e is a natural constant, mu is the dynamic viscosity coefficient of fluid in the pore, and delta P is Pj-Pi-Pc。
8. A dynamic fracture determination device for a displacement simulation mesopore throat network model is characterized by comprising:
the model construction module is used for acquiring the characteristic parameters of the porous medium and constructing a two-dimensional pore throat network model according to the characteristic parameters;
a determination module for determining that a wetting phase of the displacement simulation is water and a non-wetting phase is oil;
the first assignment module is used for assigning the upper and lower boundary states of each grid in the two-dimensional pore throat network model as fixed values, and assigning the left and right boundary states as variable values; the fixed value represents that the state of the boundary is fixed during the displacement simulation, and the variable value represents that the state of the boundary can be changed during the displacement simulation;
the first assignment module is used for assigning the states of all pore units and throat units of the two-dimensional pore throat network model to be 0 values, wherein the 0 values represent closed states and are full of water;
the seepage simulation module is used for displacement simulation of the seepage rule of the two-dimensional pore throat network model;
and the crack judging module is used for opening the corresponding throat unit to form a crack when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach the critical pressure value in the process of simulating the seepage rule so as to update the two-dimensional pore throat network model.
9. A dynamic fracture determination apparatus for a displacement simulation of a mesoporous throat network model, comprising a memory, a processor, and a computer program stored on the memory, wherein the computer program when executed by the processor performs the steps of:
acquiring characteristic parameters of the porous medium, and constructing a two-dimensional pore throat network model according to the characteristic parameters;
determining that the wetting phase of the displacement simulation is water and the non-wetting phase is oil;
assigning the upper and lower boundary states of each grid in the two-dimensional pore throat network model as fixed values, and assigning the left and right boundary states as variable values; the fixed value represents that the state of the boundary is fixed during the displacement simulation, and the variable value represents that the state of the boundary can be changed during the displacement simulation;
assigning the states of all pore units and throat units of the two-dimensional pore throat network model to be 0 values, wherein the 0 values represent closed states and are full of water;
displacement simulation of the seepage rule of the two-dimensional pore throat network model;
in the process of simulating the seepage rule, when the vertex pressure of the throat unit in the two-dimensional pore throat network model is confirmed to reach a critical pressure value, the corresponding throat unit is opened to form a crack so as to update the two-dimensional pore throat network model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810901576.3A CN109033672B (en) | 2018-08-09 | 2018-08-09 | Dynamic crack determination method and device for displacement simulation mesoporous throat network model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810901576.3A CN109033672B (en) | 2018-08-09 | 2018-08-09 | Dynamic crack determination method and device for displacement simulation mesoporous throat network model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109033672A CN109033672A (en) | 2018-12-18 |
| CN109033672B true CN109033672B (en) | 2022-01-04 |
Family
ID=64633316
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810901576.3A Active CN109033672B (en) | 2018-08-09 | 2018-08-09 | Dynamic crack determination method and device for displacement simulation mesoporous throat network model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109033672B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110853138B (en) * | 2019-11-21 | 2023-08-18 | 科吉思石油技术咨询(北京)有限公司 | Construction method of dual-medium carbonate pore-crack dual-network model |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7966164B2 (en) * | 2005-12-05 | 2011-06-21 | Shell Oil Company | Method for selecting enhanced oil recovery candidate |
| CN103257081A (en) * | 2013-04-28 | 2013-08-21 | 北京大学 | Method and device for recovering oil and gas reservoir rock mechanics underground in-situ model |
| CN103822865A (en) * | 2014-03-20 | 2014-05-28 | 中国石油大学(华东) | High-resolution three-dimensional digital rock core modeling method |
| CN108153933A (en) * | 2017-11-30 | 2018-06-12 | 中国石油天然气股份有限公司 | Production prediction method and device for tight reservoir |
-
2018
- 2018-08-09 CN CN201810901576.3A patent/CN109033672B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7966164B2 (en) * | 2005-12-05 | 2011-06-21 | Shell Oil Company | Method for selecting enhanced oil recovery candidate |
| CN103257081A (en) * | 2013-04-28 | 2013-08-21 | 北京大学 | Method and device for recovering oil and gas reservoir rock mechanics underground in-situ model |
| CN103822865A (en) * | 2014-03-20 | 2014-05-28 | 中国石油大学(华东) | High-resolution three-dimensional digital rock core modeling method |
| CN108153933A (en) * | 2017-11-30 | 2018-06-12 | 中国石油天然气股份有限公司 | Production prediction method and device for tight reservoir |
Non-Patent Citations (1)
| Title |
|---|
| 基于孔隙网络模型的砂砾岩微观水驱油驱替特征;鲁国明;《特种油气藏》;20100421;第17卷(第1期);第64-67页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109033672A (en) | 2018-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102211374B1 (en) | Method and device for outputting risk information and constructing risk information | |
| CN111144055B (en) | Method, device and medium for determining concentration distribution of toxic heavy gas leakage in urban environment | |
| CN105704188B (en) | Using the dispositions method and device with service | |
| CN109190183B (en) | Method and device for determining macroscopic parameters in displacement simulation of pore throat network model | |
| CN109063346B (en) | Displacement simulation method and device of pore throat network model considering dynamic cracking | |
| CN106940553A (en) | Industrial flow control management method and device based on industry internet operating system | |
| CN109033672B (en) | Dynamic crack determination method and device for displacement simulation mesoporous throat network model | |
| CN109063347B (en) | Displacement simulation method and device based on pore throat network model | |
| CN112927769B (en) | MOFs material defect structure prediction method based on pore size distribution curve | |
| CN109063348B (en) | Displacement simulation method and device based on pore throat network model | |
| Khalid et al. | Development of a building information modelling (BIM)-based real-time data integration system using a building management system (BMS) | |
| CN108595782B (en) | Calculation method for mass transfer between matrix and cracks in discrete cracks | |
| CN109635405B (en) | Multiphase flow transient calculation method and system based on space conservation | |
| CN109542754A (en) | The cut-in method of mobile device, the adjustment method of application program, device and system | |
| CN109117551B (en) | Method and device for determining oil displacement efficiency in pore throat network model displacement simulation | |
| CN110969249A (en) | Production well yield prediction model establishing method, production well yield prediction method and related device | |
| CN114239216A (en) | Modelica-based hot fluid medium simulation real-time method | |
| Uekermann et al. | Towards partitioned fluid-structure interaction on massively parallel systems | |
| CN109992858B (en) | SPH-based ingress and egress flow boundary calculation method and device and storage medium | |
| CN111677486A (en) | Compact oil carbon dioxide huff and puff simulation method and device and storage medium | |
| CN116362031A (en) | Method and device for analyzing permeability evolution law of hydrate reservoir | |
| CN115522909A (en) | Single well model size determination method and device based on numerical simulation | |
| CN114818514A (en) | Method and processor for determining the reservoir capacity of a gas reservoir | |
| US11366943B2 (en) | Platform for design and prototyping of micro paper based devices | |
| Baiamonte | Discussion of “Hydraulic Model of Transition of Transient to Steady Flows in the Vadose Zone” by Yaguo Zhang, Tonglu Li, Wei Shen, and Yu Wang |
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 |