CN110717295A - Method for tracking streamline distribution of tight sandstone reservoir by using finite element method - Google Patents
Method for tracking streamline distribution of tight sandstone reservoir by using finite element method Download PDFInfo
- Publication number
- CN110717295A CN110717295A CN201910954959.1A CN201910954959A CN110717295A CN 110717295 A CN110717295 A CN 110717295A CN 201910954959 A CN201910954959 A CN 201910954959A CN 110717295 A CN110717295 A CN 110717295A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- coordinates
- streamline
- sandstone reservoir
- real space
- 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.)
- Granted
Links
Images
Abstract
The invention relates to a method for tracking streamline distribution of a compact sandstone reservoir in a two-dimensional unstructured grid by using a finite element method, which comprises the following steps of: establishing a numerical simulation model of the tight sandstone reservoir; given a real space initial point P0Coordinates, obtaining P0The mesh number, the three vertex coordinates and the three vertex speeds of the triangular mesh in which the triangular mesh is positioned; converting the speed in the Master Element by using the three vertex speeds of the Jacobian inverse matrix; calculating an initial point P by Shape Function0The speed of (d); calculating an exit coordinate; converting exit coordinates to real space coordinates Pe(ii) a Judging whether to leave the system; if leaving the grid system, starting to track the next streamline, otherwise PeAssign to P0And repeating the steps until the flow line is traced. The method can track the streamline in the triangular mesh system by using a finite element method, and accurately track the streamline in the compact reservoirAnd (4) distribution situation.
Description
Technical Field
The invention relates to a method for tracking the streamline distribution of a compact sandstone reservoir by using a finite element method, belonging to the technical field of unconventional oil and gas exploration and development.
Background
In the 50 s of the 20 th century, with the development of numerical simulation methods such as a flow pipe, an overlay theory, a well mirror image method and the like, a streamline simulation method starts to be developed in the petroleum field, and the method can effectively track a fluid path and visually research fluid velocity distribution in the calculation field. In the 90's of the 20 th century, a number of scholars began to propose modern streamline simulation concepts and began to study modern streamline simulation methods which put forward the term "time of flight" (TOF, denoted by the symbol τ) for describing the time of motion of particles in the motion domain. According to this definition, a coordinate transfer function is proposed
The conversion function can convert a multidimensional operator ▽ into a one-dimensional TOF coordinate system, which generates a key principle, and the modern streamline simulation method can simplify a two-dimensional or three-dimensional problem into a series of one-dimensional subproblems along a streamline, wherein the multidimensional problem is a problem under a space (x, y) coordinate system, and the one-dimensional subproblems are problems under a time coordinate (TOF, tau) systemMigration and reaction migration.
According to most of the streamline simulation methods, the classical principle of streamline simulation mainly comprises the following steps:
(1) tracking a streamline by using the concept of 'flight time', and constructing a one-dimensional coordinate system related to the streamline;
(2) describing transport equations (e.g., mass or energy transfer equations) with TOF;
(3) periodically updating the streamlines to account for the evolution of the velocity field;
(4) numerically solving or analytically solving a one-dimensional expression of the transport equation along each streamline, and expressing the one-dimensional expression by TOF;
(5) computing the transverse flux (or other cross-streamline behavior) uses an operator splitting technique.
According to the above principle, the most important step is to trace the streamline based on the velocity field obtained by solving the continuity equation and the material equation in the porous medium by the finite element method (Darcy's law in the porous medium). Once the TOF grid system along each streamline is established, one-dimensional mass transfer equations for TOF can be solved. A multi-dimensional grid system formed by finite element methods can be used to solve the flow problem and trace streamlines, and each streamline can be discretized into a series of TOF-based computing nodes.
In the aspect of calculating the streamline, a Runge-Kutta type solver is generally used for numerical tracking, but because the difficulty of generating the TOF grid by the Runge-Kutta type solver is high, and the numerical tracking cannot be performed by a longge-coutta solver in streamline simulation, a semi-analytic method such as a Pollock algorithm is widely used. These methods are referred to as semi-analytic methods because the analytic solution is derived from the control body, not the entire domain. Pollock proposes a piecewise analytical solution, derived by TOF, when assuming that the velocity of each direction in the velocity field varies linearly, and is independent of the velocities of other directions. Given the start point coordinates and velocity for any hypothetical fluid particle, the algorithm can determine the exit coordinates from the grid and the time it takes to pass through the grid (this is TOF and is also the residence time of the particle within a grid). Although Pollock can efficiently generate TOF grids while tracing streamlines and can reasonably assume patterns for streamlines in each computational domain, it is only applicable to structured grids, and it cannot trace streamlines when there are sources or sinks (represented by the flow in the grid) within the grid. Although streamline simulation in unstructured grids has been studied by many scholars, these algorithms suffer from at least one of the following disadvantages:
(1) these formulas rely on the concept of a control volume, requiring flow to be applied to each face of the mesh, and a certain flow in each mesh;
(2) in this process, ordinary differential equation numerical solution is involved, and TOF grid cannot be calculated and obtained efficiently.
Therefore, the streamline obtained by the current method can only be used for visualizing the velocity field, and the one-dimensional mass transfer equation cannot be mapped on the streamline.
Disclosure of Invention
Aiming at the problems, the invention mainly overcomes the defects in the prior art and provides a method for tracking the streamline distribution of a compact sandstone reservoir by using a finite element method.
The technical scheme provided by the invention for solving the technical problems is as follows: a method for tracking the streamline distribution of a tight sandstone reservoir by using a finite element method comprises the following steps:
s1, establishing a numerical simulation model of the compact sandstone reservoir of the two-dimensional unstructured grid by using a finite element method, and acquiring a pressure field and a velocity field of the compact sandstone reservoir according to parameters of the compact sandstone reservoir;
step S2, setting the injection as the real space initial point P0(x0,y0) (ii) a Obtaining P0The number of the located grid is CellID, and the coordinates (x) of three vertexes of the grid are obtained1,y1)，(x2,y2)，(x3,y3) And vertex velocity (v)x1,vy1)，(vx2,vy2)，(vx3,vy3)；
Step S3, converting the real data into real data through a coordinate conversion formulaInitial point in space P0(x0,y0) Conversion to Master Element as (eta)0,ξ0)；
Step S4, converting the three vertex speeds of the grid into (v) by using a Jacobian inverse matrixξ1,vη1)，(vξ2,vη2)，(vξ3,vη3) (ii) a The velocity (v) of the initial point is calculated by Shape Functionξ0,vη0)；
Step S5, solving a track equation in the MasterElement, obtaining the time of the particles to reach three sides of the triangular unit through a Newton-Raphson iterative algorithm, taking the minimum positive value as the flight time delta tau in the unit, and calculating an exit coordinate Pe *(ηe,ξe)；
Step S6, utilizing Shape Function to convert Pe *(ηe,ξe) Conversion to real space coordinates Pe(xe,ye) And accumulating the flight time delta tau;
step S7, judgment Pe(xe,ye) Whether to leave the grid system;
step S8, if leaving the grid system, starting to track the next streamline; otherwise P will bee(xe,ye) Assign to P0(x0,y0) And repeating the steps S2 to S7 until the flow line is traced.
Further technical solution is that, in the step S2, the real space initial point P is obtained under the FEniCS system through programming0(x0,y0) Three vertex coordinates (x) of the located grid1,y1)，(x2,y2)，(x3,y3) And vertex velocity (v)x1,vy1)，(vx2,vy2)，(vx3,vy3)。
The further technical solution is that in step S3, the expression of η and ξ is solved by the following formula:
x＝x1(1-ξ-η)+x2ξ+x3η
y＝y1(1-ξ-η)+y2ξ+y3η
in the formula: x is the x coordinate of the real space, y is the y coordinate of the real space, and xi is the transverse coordinate in the Master Element (as shown in FIG. 3), corresponding to the x coordinate; η is the longitudinal coordinate in the Master Element (as shown in FIG. 3), corresponding to the y coordinate; (x)1，y1)，(x2，y2)，(x3，y3) The coordinates of three vertexes of the triangle in the real space are respectively converted into (0, 0), (1, 0), (0,1) in the MasterElement.
The further technical scheme is that the specific process of the step S4 is as follows:
converting the three vertex velocities of the grid into (v) in the MasterElement by using a coordinate conversion formulaξ1,vη1)，(vξ2,vη2)，(vξ3,vη3) (ii) a Calculating an initial point P by Shape Function0Velocity (v) ofξ0,vη0)。
Wherein the Jacobian inverse matrix is:
wherein
In the formula: (x)1，y1)，(x2，y2)，(x3，y3) Coordinates of three vertexes of a triangle in a real space are respectively; upsilon isxAnd upsilonyVelocity in real space x and y, respectively, upsilonξAnd upsilonηThe velocities in the xi and eta directions in the MasterElement.
The further technical scheme is that the specific process of the step S5 is as follows:
step S51, setting particle slave P0Starting from the triangle, P can be calculated separately from any side of the triangle0The time from the side 1, the side 2, and the side 3 is Δ τ1、Δτ2、Δτ3；
Step S52, three edge time delta tau is taken1、Δτ2、Δτ3Is the time of flight Δ τ within the cell;
step S53, determining departure coordinate P in MasterElemente(ξe,ηe) Wherein the flight time delta tau obtained in the second step is substituted into the following equation to calculate xie,ηe；
Wherein
In the formula: (v)ξ1,υη1)，(υξ2,υη2)，(υξ3,υη3) Respectively the speeds in the xi direction and the eta direction of three vertexes of a triangle in the MasterElement; lambda [ alpha ]1，2Is a characteristic equation
Two characteristic values of
e is the Euler number, and has a value of about 2.7183; Δ τ is the time of flight; xie,ηeThe exit coordinates of the particles in the Master Element from the triangle.
Step S54, determining departure coordinate P under real coordinate systeme *(ηe,ξe) Xi in the third stepe,ηeConversion to real coordinates P by Shape Functione *(ηe,ξe)。
The invention has the following beneficial effects:
(1) similar to Pollock algorithm, the exit coordinates can be determined from the entry coordinates in the triangular unit mesh;
(2) the calculation only needs the data (coordinates and flow) of the vertex of the triangular mesh, and does not need to consider the flow on the boundary of the triangular unit;
(3) the invention provides an analytic solution for describing the streamline locus in the triangular unit, so that the TOF can be accurately and efficiently determined.
Drawings
Obtaining a grid vertex coordinate code of an arbitrary point in the FENICS system in FIG. 1;
FIG. 2 is a graph showing a graph vertex velocity code of a mesh vertex where an arbitrary point is located in the FENICS system;
FIG. 3 is a triangular solid graph and coordinate distribution in Master element;
FIG. 4 is a solution process for streamline tracing in unstructured mesh triangle cells generated using finite elements;
FIG. 5 validation of results of streamline tracing in triangle unit in MasterElement;
FIG. 6 is a technical route diagram of a tight sandstone reservoir utilizing a finite element tracing streamline method according to an embodiment of the invention;
FIG. 7 is a computational grid diagram of a tight sandstone reservoir in an embodiment of the present invention;
figure 8 is a pressure field diagram of a tight sandstone reservoir in an embodiment of the present invention;
FIG. 9 is a velocity field plot for a tight sandstone reservoir in an embodiment of the present invention;
figure 10 is a tight sandstone reservoir flowline profile according to an embodiment of the present invention.
Detailed Description
The present invention will be further described with reference to the following examples and the accompanying drawings.
Example 1
The invention discloses a method for tracking streamline distribution of a compact sandstone reservoir by using a finite element method, which comprises the following steps of:
step S1, firstly, establishing a compact sandstone reservoir numerical simulation model of a two-dimensional unstructured grid (triangular grid) by using a finite element method;
wherein the size of the compact sandstone reservoir is 500m multiplied by 500m, the reservoir permeability is 5mD, 1 injection well and 4 production wells are arranged, the bottom pressure of the injection well is 14MPa, the bottom pressure of the production well is 12.5MPa, and a pressure field (shown in figure 8) and a velocity field (shown in figure 9) of the compact sandstone reservoir are obtained;
step S2, injection well is given as real space initial point P0(x0,y0) Obtaining P0The grid number CellID is located, and the coordinates (x) of three vertexes of the grid are read1,y1)，(x2,y2)，(x3,y3) And vertex velocity (v)x1,vy1)，(vx2,vy2)，(vx3,vy3)；
The vertex coordinates of the mesh where the point is located are obtained under the FEniCS system through programming, and the code is as follows (as shown in fig. 1, programmed by using Python language):
the vertex speed of the mesh where the point is located is obtained under the FEniCS system through programming, and the code is as follows (as shown in fig. 2, using Python language for programming):
s3, converting the real space initial point through a coordinate conversion formulaP0(x0,y0) Conversion to MasterElement of (η)0,ξ0)；
Wherein by solving for the expression of η and ξ:
x＝x1(1-ξ-η)+x2ξ+x3η
y＝y1(1-ξ-η)+y2ξ+y3η
that is, the initial point P in real space can be obtained0(x0,y0) Conversion to MasterElement of (η)0,ξ0)。
S4, converting the three vertex velocities (v) of the grid by using a Jacobian inverse matrixξ1,vη1)，(vξ2,vη2)，(vξ3,vη3) (ii) a The velocity (v) of the initial point is calculated by Shape Functionξ0,vη0)；
The MasterElement rates were:
wherein
S5, solving a track equation in Master Element, obtaining the time of the particles to reach three edges of a triangular unit through a Newton-Raphson iterative algorithm, taking the minimum positive value as the flight time delta tau in the unit, and calculating an exit coordinate Pe *(ηe,ξe)；
Master Element is a right-angled isosceles triangle with three vertexes as P1(0,0)，P2(1,0)，P3(0, 1); the three sides are respectively: edge 1 (P)1-P2) Edge 2 (P)2-P3) Edge 3 (P)3-P1) As shown in fig. 3.
P in one particle from Master Element Unit0From a point P on one side of the triangleeThe time taken for departure, namely the time of flight (TOF), is denoted by Δ τ.
The trajectory equation is solved as follows:
and solving a trajectory equation of any point P in a standard triangle in the Master element according to the first-order Lagrange Shape Function.
Three vertex coordinates P of triangle1(0,0),P2(1,0),P2(0,1), where any point P is (xi, eta), the initial condition is t ═ 0, xi ═ xi0,ξ'＝υξ0,η＝η0,η'＝υη0。
Velocity upsilon of arbitrary point P (xi, eta)ξIs composed of
Wherein
Since the velocity distribution in a velocity field is arbitrary, the coefficients of the equation have uncertainty, and it is necessary to discuss the equation coefficients to solve the equation:
(1) when b is2Is equal to 0, and c1When 0, the differential equation can be converted to a directly separable set of variable differential equations whose solution is:
(2) when b is2When not equal to 0, the differential equation is a non-homogeneous linear equation set,
when Δ ═ b1-c2)2+4b2c1＞0
Specific solutions of two distinct solid roots
Wherein
(3) When Δ ═ b1-c2)2+4b2c1＝0，
Special solution of two equal root
Wherein:
(4) when Δ ═ b1-c2)2+4b2c1＜0
Special solution of two conjugate complex roots
Wherein
Specifically, according to the master element trajectory equation, the coordinate PeAnd time of flight Δ τ can be calculated by the following four steps:
the first step is as follows: suppose that this particle is from P0Starting from the triangle, P can be calculated separately from any side of the triangle0The time from the side 1, the side 2, and the side 3 is Δ τ1、Δτ2、Δτ3。
Given a position, given equation (4) as an example, the time for the particle to reach this position can be determined, and thus given η equal to 0, Δ τ can be determined1(ii) a ξ + η ═ 1 can be determinedΔτ2(ii) a ξ ═ 0, Δ τ can be determined3。
Based on these principles, the following equation (7) can be summarized, where Δ τ1、Δτ2、Δτ3Obtained by the Newton-Raphson algorithm.
The second step is that: from Δ τ1、Δτ2、Δτ3Finding the true time of flight Δ τ. The calculated TOF value may be positive or negative and the true time of flight Δ τ is the smallest positive value, denoted Δ τ -min (Δ τ)1,Δτ2,Δτ3)。
The third step: determining departure coordinates P in a master elemente(ηe,ξe) Substituting the time of flight delta tau obtained in the second step into equation (8) to calculate xie,ηe
The fourth step: determining departure coordinates P in a real coordinate systeme *(ηe,ξe) I.e. xi in the third stepe,ηeConversion to real coordinates P by Shape Functione *(ηe,ξe)；
S6, using Shape Function to convert Pe *(ηe,ξe) Conversion to real space coordinates Pe(xe,ye) And accumulating the flight time delta tau;
converting any point P (eta, xi) in the MasterElement into a formula in the Actual Element:
x＝x1(1-ξ-η)+x2ξ+x3η
y＝y1(1-ξ-η)+y2ξ+y3η
s7, judgment Pe(xe,ye) Is thatIf not, leaving the system;
s8, if leaving the grid system, starting to track the next streamline, otherwise, Pe(xe,ye) Assign to P0(x0,y0) Steps S2 to S7 are repeated, with 200 streamlines traced starting from the injection well, with the results shown in fig. 10.
When initial conditions are given, the motion trajectory in a triangle unit is tracked by using the above formula and theory, as shown in fig. 4, fig. 4 shows the solving process of streamline tracing in an unstructured mesh triangle unit generated by using finite elements when Δ >0, fig. 4(a) shows the coordinates and the velocity of the vertex and the starting point in Actual Element, then the information in fig. 4(a) is converted from Actual Element to Master Element, and the motion trajectory of the particle is tracked, and as a result, as shown in fig. 4(b), the result needs to be converted into Actual Element after the Master Element tracing is finished, as shown in fig. 4 (c). The same applies when Δ is 0 and Δ < 0.
The results of the patent of the invention are verified through 4-order Runge-Kutta numerical algorithm results, as shown in FIG. 5, the lines in the graph represent the patent results of the invention, the black dots represent the 4-order Runge-Kutta numerical algorithm results, the errors of the two results are very small, and the accuracy of the patent method of the invention is verified.
Although the present invention has been described with reference to the above embodiments, it should be understood that the present invention is not limited to the above embodiments, and those skilled in the art can make various changes and modifications without departing from the scope of the present invention.
Claims (5)
1. A method for tracking the streamline distribution of a tight sandstone reservoir by using a finite element method is characterized by comprising the following steps:
s1, establishing a numerical simulation model of the compact sandstone reservoir of the two-dimensional unstructured grid by using a finite element method, and acquiring a pressure field and a velocity field of the compact sandstone reservoir according to parameters of the compact sandstone reservoir;
step S2, setting the position of the injection well as the real space initial point P0(x0,y0) (ii) a Obtaining P0The number Cell ID of the mesh where the grid is located, and the coordinates (x) of three vertexes of the mesh are obtained1,y1)，(x2,y2)，(x3,y3) And vertex velocity (v)x1,vy1)，(vx2,vy2)，(vx3,vy3)；
Step S3, converting the real space initial point P by the coordinate conversion formula0(x0,y0) Converting into Master Element space coordinate (eta)0,ξ0)；
Step S4, converting the three vertex speeds of the grid into (v) in Master Element by using Jacobian inverse matrixξ1,vη1)，(vξ2,vη2)，(vξ3,vη3) (ii) a Calculating an initial point P by Shape Function0Velocity (v) ofξ0,vη0)；
Step S5, solving a track equation in the Master Element, obtaining the time of the particles to reach three edges of the triangular unit in the Master Element through a Newton-Raphson iterative algorithm, taking the minimum positive value as the flight time delta tau in the unit, and calculating an outlet coordinate Pe *(ηe,ξe)；
Step S6, utilizing Shape Function to convert Pe *(ηe,ξe) Conversion to real space coordinates Pe(xe,ye) And accumulating the flight time delta tau;
step S7, judgment Pe(xe,ye) Whether to leave the grid system;
step S8, if leaving the grid system, starting to track the next streamline; otherwise P will bee(xe,ye) Assign to P0(x0,y0) Repeating steps S2 to S7 until theAnd finishing tracing the streamline.
2. The method for tracking streamline distribution of tight sandstone reservoir by using finite element method as claimed in claim 1, wherein said step S2 is implemented by obtaining real space initial point P under FENICS system0(x0,y0) Three vertex coordinates (x) of the located grid1,y1)，(x2,y2)，(x3,y3) And vertex velocity (v)x1,vy1)，(vx2,vy2)，(vx3,vy3)。
3. The method for tracking streamline distribution of tight sandstone reservoir according to claim 1, wherein the expression of η and ξ is solved in step S3 by the following formula:
x＝x1(1-ξ-η)+x2ξ+x3η
y＝y1(1-ξ-η)+y2ξ+y3η
in the formula: x is the x coordinate of the real space, y is the y coordinate of the real space, xi is the transverse coordinate in the Master Element, and corresponds to the x coordinate; eta is a longitudinal coordinate in the Master Element and corresponds to a y coordinate; (x)1，y1)，(x2，y2)，(x3，y3) The coordinates of three vertexes of the triangle in the real space are respectively converted into (0, 0), (1, 0), (0,1) in the Master Element.
4. The method for tracking streamline distribution of tight sandstone reservoir by finite element method according to claim 1, wherein the method is characterized in thatCharacterized in that, in the step S4, the three vertex speeds of the grid are converted into (v) in the Master Element by using a Jacobian inverse matrixξ1,vη1)，(vξ2,vη2)，(vξ3,vη3) (ii) a Calculating an initial point P by Shape Function0Velocity (v) ofξ0,vη0)。
The Jacobian inverse matrix is:
wherein
In the formula: (x)1，y1)，(x2，y2)，(x3，y3) Coordinates of three vertexes of a triangle in a real space are respectively; upsilon isxAnd upsilonyVelocity in real space x and y, respectively, upsilonξAnd upsilonηThe velocities in the xi and eta directions in the Master Element.
5. The method for tracking the streamline distribution of the tight sandstone reservoir by using a finite element method as claimed in claim 1, wherein the specific process of the step S5 is as follows:
step S51, setting particle slave P0Starting from the triangle, P can be calculated separately from any side of the triangle0The time from the side 1, the side 2, and the side 3 is Δ τ1、Δτ2、Δτ3；
Step S52, taking the time delta tau of three edges of the particle1、Δτ2、Δτ3Is the time of flight Δ τ within the cell;
step S53, determining departure coordinate P in Master Elemente(ξe,ηe) Wherein the flight time delta tau obtained in the second step is substituted into the following equation to calculate xie,ηe；
Wherein
In the formula: (v)ξ1,υη1)，(υξ2,υη2)，(υξ3,υη3) Respectively the speeds in the xi direction and the eta direction of three vertexes of a triangle in the Master Element; lambda [ alpha ]1，2Is a characteristic equation
Two characteristic values of
e is the Euler number, and has a value of about 2.7183; Δ τ is the time of flight; xie,ηeThe exit coordinates of the particles in the Master Element from the triangle.
Step S54, determining departure coordinate P under real coordinate systeme *(ηe,ξe) Xi in the third stepe,ηeConversion to real coordinates P by Shapefunctione *(ηe,ξe)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910954959.1A CN110717295B (en) | 2019-10-09 | 2019-10-09 | Method for tracking streamline distribution of tight sandstone reservoir by using finite element method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910954959.1A CN110717295B (en) | 2019-10-09 | 2019-10-09 | Method for tracking streamline distribution of tight sandstone reservoir by using finite element method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110717295A true CN110717295A (en) | 2020-01-21 |
| CN110717295B CN110717295B (en) | 2020-09-08 |
Family
ID=69212319
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910954959.1A Active CN110717295B (en) | 2019-10-09 | 2019-10-09 | Method for tracking streamline distribution of tight sandstone reservoir by using finite element method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110717295B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020199156A1 (en) * | 2001-05-23 | 2002-12-26 | Fuel Tech, Inc. | Hardware-adaptable data visualization tool for use in complex data analysis and engineering design |
| CN105067795A (en) * | 2015-07-29 | 2015-11-18 | 中国石油大学(北京) | Method for evaluating microstructure of compact reservoir rock |
| CN108828190A (en) * | 2018-06-22 | 2018-11-16 | 中国石油大学(北京) | A kind of results of fracture simulation method of Fractured tight sand oil-gas reservoir |
| CN108846224A (en) * | 2018-06-27 | 2018-11-20 | 中国人民解放军国防科技大学 | Supersonic flow channel design method and device |
| CN109102558A (en) * | 2018-07-09 | 2018-12-28 | 珠江水利委员会珠江水利科学研究院 | A kind of water flow streamline rapid generation based on the non-structural triangulation network |
| JP2019067276A (en) * | 2017-10-04 | 2019-04-25 | 日本電信電話株式会社 | Apparatus and method for fluid mass extraction |
-
2019
- 2019-10-09 CN CN201910954959.1A patent/CN110717295B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020199156A1 (en) * | 2001-05-23 | 2002-12-26 | Fuel Tech, Inc. | Hardware-adaptable data visualization tool for use in complex data analysis and engineering design |
| CN105067795A (en) * | 2015-07-29 | 2015-11-18 | 中国石油大学(北京) | Method for evaluating microstructure of compact reservoir rock |
| JP2019067276A (en) * | 2017-10-04 | 2019-04-25 | 日本電信電話株式会社 | Apparatus and method for fluid mass extraction |
| CN108828190A (en) * | 2018-06-22 | 2018-11-16 | 中国石油大学(北京) | A kind of results of fracture simulation method of Fractured tight sand oil-gas reservoir |
| CN108846224A (en) * | 2018-06-27 | 2018-11-20 | 中国人民解放军国防科技大学 | Supersonic flow channel design method and device |
| CN109102558A (en) * | 2018-07-09 | 2018-12-28 | 珠江水利委员会珠江水利科学研究院 | A kind of water flow streamline rapid generation based on the non-structural triangulation network |
Non-Patent Citations (1)
| Title |
|---|
| 唐洪明，等: "致密砂岩气层水锁机理及控制因素研究", 《中国科学:技术科学》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110717295B (en) | 2020-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kallinderis et al. | Adaptation methods for a new Navier-Stokes algorithm | |
| Bernard et al. | High-order h-adaptive discontinuous Galerkin methods for ocean modelling | |
| Bahar et al. | Numerical solution of advection-diffusion equation using operator splitting method | |
| Nishikawa et al. | A third-order fluctuation splitting scheme that preserves potential flow | |
| Pandya et al. | Accuracy, Scalability, and Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows | |
| CN110717295B (en) | Method for tracking streamline distribution of tight sandstone reservoir by using finite element method | |
| Cummings et al. | Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies | |
| Zhang et al. | Machine learning methods for turbulence modeling in subsonic flows over airfoils | |
| CN111210522B (en) | Method for tracking streamline distribution in three-dimensional unstructured grid flow field by using FEM (finite element modeling) | |
| Bisson et al. | Adjoint-based aerodynamic optimization framework | |
| Coirier et al. | A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations | |
| Paulete-Periáñez et al. | Surrogate modelling for aerodynamic coefficients prediction in aeronautical configurations | |
| Ito et al. | Solution adaptive mesh generation using feature-aligned embedded surface meshes | |
| Iaccarino et al. | Local grid refinement for an immersed boundary RANS solver | |
| Chandar | Development of a Parallel Overset Grid Framework for Moving Body Simulations in Open FOAM. | |
| Harris | Flow feature aligned mesh generation and adaptation | |
| CN106066912A (en) | A kind of generation method of multi partition structured grid | |
| Sadrehaghighi | Dynamic & Adaptive Meshing | |
| Villamizar et al. | Generation of smooth grids with line control for scattering from multiple obstacles | |
| Katz et al. | Edge-based meshless methods for compressible viscous flow with applications to overset grids | |
| Yamakawa et al. | Unstructured moving-grid finite-volume method for unsteady shocked flows | |
| McCaslin | Development And Application Of Numerical Methods For Interfacial Dynamics In Turbulent Liquid-Gas Flows | |
| Kandil et al. | Sonic boom computations for double-cone configuration using cfl3d, fun3d and full-potential codes | |
| Maruyama | Robust Multi-Objective Optimization in Aerodynamics using MGDA | |
| Pattinson | A cut-cell, agglomerated-multigrid accelerated, Cartesian mesh method for compressible and incompressible flow |
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 |