CN106934185B - A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid - Google Patents
A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid Download PDFInfo
- Publication number
- CN106934185B CN106934185B CN201710285521.XA CN201710285521A CN106934185B CN 106934185 B CN106934185 B CN 106934185B CN 201710285521 A CN201710285521 A CN 201710285521A CN 106934185 B CN106934185 B CN 106934185B
- Authority
- CN
- China
- Prior art keywords
- multiple dimensioned
- matrix
- basic function
- scale
- large scale
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
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 invention discloses a kind of multiple dimensioned flow simulating methods of the fluid structurecoupling of elastic fluid.This method comprises: obtaining reservoir geology parameter and mechanics parameter, oil reservoir geometrical model is established, and multiple dimensioned grid dividing is carried out to oil reservoir geometrical model, obtain multiple dimensioned grid system;Based on multiple dimensioned grid system, establish mass-conservation equation and momentum conservation equation, realize the real simulation of changeability medium seepage flow situation, local momentum conservation equation is solved using multi-scale finite element method, obtain multiple dimensioned displacement basic functions, using multi-scale Simulation finite difference method Part mass conservation equation, multiple dimensioned speed basic function and multiple dimensioned pressure basic function are obtained, to guarantee the conservativeness calculated and accuracy;The flow simulating on large scale coarse grid is realized by multiple dimensioned basic function, and according to the mapping relations between large scale solution and small scale solution, the small scale solution of refined net unit is obtained, to reduce calculation amount.
Description
Technical field
The present invention relates to the field that porous media flow consolidates Coupled Flow simulation, in particular to a kind of stream of elastic fluid consolidates coupling
Close multiple dimensioned flow simulating method.
Background technique
Underground medium has strong stress sensitivity, and the stress in oil reservoir is to capability forecasting, uncertainty analysis and wind
Danger assessment etc. has great influence, needs to fully consider the influence of stress in reservoir numerical simulation.But it is answered in practical oil reservoir
In, since the multiple dimensioned property and oil reservoir space scale of reservoir geology and geometric properties are big, the time scale of simulation is long, makes
Conventional numeric method first to region carry out fine grid blocks subdivision, then carry out on a fine scale simulation calculation amount it is huge, surpass
Existing computing capability is gone out.In recent years, with the development of multi-scale method, domestic and international experts and scholars start multi-scale method
Applied in the non-linear flow simulating of heterogeneous reservoir.Multi-scale method constructs large scale by the multiple dimensioned basic function of local calculation
Equation group, after obtaining large scale solution, multi-scale method can map to obtain accuracy higher based on multiple dimensioned basic function
Small scale finely solves, and multi-scale method is the calculating speed for possessing scale upgrade method, and possesses very high computational accuracy, but
It is that now multiple dimensioned development focuses primarily upon Multiphase Flow simulation, especially ellipse pressure equation solves, and is seldom related to being situated between
The factor of the mechanical deformation of matter, and only the multiple dimensioned golden finite element method of gal the Liao Dynasty is used in asking for fluid structurecoupling problem at present
Solution, but multi-level finite element modeling method is the frame based on the golden finite element of gal the Liao Dynasty, therefore in the building of multiple dimensioned speed basic function
It cannot be guaranteed speed part conservativeness, be unsuitable for the calculating of multiphase flow.
Therefore, find a kind of algorithm can guarantee that the New Type of Numerical method of computational accuracy is imperative again.
Summary of the invention
The object of the present invention is to which the technological deficiency in order to overcome existing fluid structurecoupling oil deposit flow simulation method, obtains
One kind not only can be reduced calculation amount, but also possesses the fluid structurecoupling flow simulating method of higher computational accuracy, provides a kind of elasticity Jie
The multiple dimensioned flow simulating method of the fluid structurecoupling of matter.
To achieve the above object, the present invention provides following schemes:
A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid, includes the following steps:
Reservoir geology parameter and mechanics parameter are obtained, oil reservoir geometrical model is established;
Multiple dimensioned grid dividing is carried out to oil reservoir geometrical model, obtains including coarse net trellis system and fine-structure mesh trellis system
Multiple dimensioned grid system;
Based on multiple dimensioned grid system, local flow region is chosen, the Part mass conservation side of the flow region is established
Journey, and local momentum conservation equation is solved using multi-scale finite element method, obtain multiple dimensioned displacement basic functions;
Based on multiple dimensioned grid system, flow region is chosen, the local momentum conservation equation of the flow region is established, adopts
With multi-scale Simulation finite difference method Part mass conservation equation, multiple dimensioned speed basic function and multiple dimensioned pressure base are obtained
Function;
Based on multiple dimensioned principle and multiple dimensioned displacement basic functions, multiple dimensioned speed basic function and multiple dimensioned pressure basic function,
Macroscopical large scale solution of coarse grid cell is obtained, and obtains the mapping relations between macroscopical large scale solution and small scale solution;
According to the mapping relations between macroscopical large scale solution and large scale solution and small scale solution, the small of refined net unit is obtained
Scale solution.
Optionally, the specific steps of the multiple dimensioned grid dividing include:
According to survey region size, the large scale coarse grid step-length and quantity of each direction in space are determined, using orthogonal net
Lattice carry out small scale fine grid blocks to oil reservoir geometrical model and divide acquisition fine-structure mesh trellis system;
Coarse net trellis system, the coarse grid subsystem are constructed using load balancing algorithm on the basis of small scale refined net
Coarse grid cell in system is connected with each other by the refined net unit in the fine-structure mesh trellis system;
Multiple dimensioned grid system is formed by the coarse net trellis system and the fine-structure mesh trellis system.
Optionally, the fine-structure mesh trellis system includes reservoir rocks, the essential characteristic parameter of fluid and mechanics parameter.
Optionally, the specific steps for obtaining multiple dimensioned displacement basic functions include:
Each coarse grid cell is considered as a local flow region Ω, the momentum for establishing the local flow region is kept
Permanent equation:
Wherein, CdrFor elasticity tensor, l is time interval l=[0, t], and u and p respectively represent motion vector and pressure, and I is
Unit matrix, b are to finish wet coefficient,For symmetric gradient operator;
Choose boundary condition φi(xj)=δij, δijFor kronecker delta, it is dynamic that part is solved using multi-scale finite element method
Conservation equation is measured, the multiple dimensioned displacement basic functions of each coarse grid cell are obtainedAnd construct multiple dimensioned displacement base
Jacobian matrixWhereinFor small scale displacement basic functions, the small scale displacement basic functions
For carry out small scale fine grid blocks calculate when offset variable selected by basic function,For i-th of fine-structure mesh in coarse grid cell
Displacement discrete value on lattice node,For the quantity of refined net cell node in coarse grid cell.
Optionally, the multiple dimensioned speed basic function and the specific steps of multiple dimensioned pressure basic function of obtaining includes: by phase
Two adjacent coarse grid cells are considered as a flow region Ω ', to guarantee local conservativeness, establish the part of the flow region
Mass-conservation equation:
Wherein, b is to finish wet coefficient, and l is time interval l=[0, t], MbTo finish wet modulus,For leading for displacement versus time
Number,For the derivative of Pressure versus Time, λ is fluid flow coefficient.
Closed boundary condition is chosen, using multi-scale Simulation finite difference method Part mass conservation equation, is obtained every
The corresponding multiple dimensioned pressure basic function of one coarse grid cellIt is corresponding with each coarse grid boundary multiple dimensioned
Speed basic functionAnd construct multiple dimensioned pressure basic function matrixWith more rulers
Spend speed basic function matrixWherein,WithRespectively small scale pressure basic function and small
Scale speed basic function, the small scale pressure basic function are the base of pressure variations when carrying out the calculating of small scale fine grid blocks
Function, the small scale speed basic function are the basic function of speed variables when carrying out the calculating of small scale fine grid blocks,It is thick
The borderline pressure discrete value of i-th of refined net in grid cell,For the borderline speed of i-th of refined net in coarse grid cell
Discrete value is spent,For refined net element number in coarse grid cell;For refined net elementary boundary quantity in coarse net unit.
Optionally, the specific steps of the macroscopical large scale solution for obtaining coarse grid cell include:
By multiple dimensioned displacement basic functions matrix, multiple dimensioned pressure basic function matrix and multiple dimensioned speed basic function matrix assembly
As mapping operator Φ,
I.e.
Construct large scale stiffness matrix A=Φ AfΦT, wherein AfFor small scale stiffness matrix, the small scale stiffness matrix
To be formed by stiffness matrix, specific element after discrete momentum conservation equation and momentum conservation equation on small scale refined net
Including mechanics stiffness matrix, finite difference coefficient matrix, and compression term system matrix number are simulated, forms large scale side on this basis
Journey group Ax=b, wherein b is equation group right-hand vector, including source sink term, displacement boundary conditions and stress boundary condition;
Coupling large scale equation group is carried out to solve the macroscopical large scale solution for obtaining coarse grid cell, including large scale displacement
uc, large scale pressure pc, large scale speed vc。
Optionally, the specific steps of the small scale solution for obtaining refined net unit include:
Construct mapping matrix Nu, Np, Nv, wherein NuIt is the matrix using all multiple dimensioned displacement basic functions as column vector, Np
It is the matrix using all multiple dimensioned pressure basic functions as column vector, NvIt is using all multiple dimensioned speed basic functions as column vector
Matrix;
The corresponding relationship between macroscopical large scale solution and the small scale solution, which is obtained, according to mapping matrix is respectively as follows: uf
=Nuuc, pf≈Ipc+NpDλvc, vf=Nvvc, and macroscopical large scale solution of coarse grid cell is combined, obtain the small of refined net unit
Scale solution, wherein uf, pf, vfSmall scale displacement, small scale pressure and small scale speed are respectively indicated, I is unit matrix, DλFor
Mobility coefficient matrix.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid of disclosure of the invention, using mass-conservation equation and
Momentum conservation equation realizes the real simulation of changeability medium seepage flow situation, solves local momentum using multi-scale finite element method and keeps
Permanent equation obtains multiple dimensioned displacement basic functions, is obtained using multi-scale Simulation finite difference method Part mass conservation equation more
Scale speed basic function and multiple dimensioned pressure basic function are consolidated coupled wave equation by multiple dimensioned basic function building large-scale currents and are obtained
Large scale coarse grid solution, including large scale displacement solution, pressure solution and velocity solution;According to reflecting between large scale solution and small scale solution
Relationship is penetrated, the small scale solution of refined net unit is obtained, greatly reduces meter while guaranteeing fluid structurecoupling flow simulating precision
Calculation amount, can more fully simulation oil field development behavior, provide technical support for the Efficient Development of oil reservoir.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is a kind of one embodiment of the multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid provided by the invention
Flow chart.
Fig. 2 is a kind of multiple dimensioned grid of the multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid provided by the invention
System schematic.
Fig. 3 is a kind of multiple dimensioned displacement of the multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid provided by the invention
The acquisition schematic diagram of basic function.
Fig. 4 is a kind of another implementation of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid provided by the invention
The flow chart of example.
Specific embodiment
The object of the present invention is to provide a kind of multiple dimensioned flow simulating methods of the fluid structurecoupling of elastic fluid.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Fig. 1 is a kind of one embodiment of the multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid provided by the invention
Flow chart, as shown in Figure 1, as an embodiment, a kind of multiple dimensioned flow simulating side of the fluid structurecoupling of elastic fluid
Method includes the following steps:
S1 obtains reservoir geology parameter and mechanics parameter, establishes oil reservoir geometrical model;
S2 carries out multiple dimensioned grid dividing to oil reservoir geometrical model, obtains including coarse net trellis system and fine-structure mesh trellis system
Multiple dimensioned grid system;
S3 is based on multiple dimensioned grid system, chooses local flow region, establishes the Part mass conservation of the flow region
Equation, and local momentum conservation equation is solved using multi-scale finite element method, obtain multiple dimensioned displacement basic functions;
S4 is based on multiple dimensioned grid system, chooses flow region, establishes the local momentum conservation equation of the flow region,
Using multi-scale Simulation finite difference method Part mass conservation equation, multiple dimensioned speed basic function and multiple dimensioned pressure are obtained
Basic function;
S5 is based on multiple dimensioned principle and multiple dimensioned displacement basic functions, multiple dimensioned speed basic function and multiple dimensioned pressure base letter
Number, obtains macroscopical large scale solution of coarse grid cell, and obtain the mapping relations between macroscopical large scale solution and small scale solution;
S6 obtains refined net unit according to the mapping relations between macroscopical large scale solution and large scale solution and small scale solution
Small scale solution.
The specific steps of multiple dimensioned grid dividing described in step S2 include:
According to survey region size, the large scale coarse grid step-length and quantity of each direction in space are determined, using orthogonal net
Lattice carry out small scale fine grid blocks to oil reservoir geometrical model and divide acquisition fine-structure mesh trellis system;
Coarse net trellis system, the coarse grid subsystem are constructed using load balancing algorithm on the basis of small scale refined net
Coarse grid cell in system is connected with each other by the refined net unit in the fine-structure mesh trellis system;
Multiple dimensioned grid system, the multiple dimensioned grid are formed by the coarse net trellis system and the fine-structure mesh trellis system
System is as shown in Figure 2.
The coarse net trellis system includes the essential attribute of rock and fluid, and the fine-structure mesh trellis system includes oil deposit parameter
And mechanics parameter.
The specific steps of the multiple dimensioned displacement basic functions of acquisition described in step S3 include:
As shown in figure 3, each coarse grid cell is considered as a local flow region Ω, the local flow area is established
The momentum conservation equation in domain:
Wherein, CdrFor elasticity tensor, l is time interval l=[0, t], and u and p respectively represent motion vector and pressure, and I is
Unit matrix, b are to finish wet coefficient,For symmetric gradient operator, it is preferable that Ω=Ωi, wherein ΩiIndicate multiple dimensioned grid system
The region of i-th of coarse grid cell in system;
Choose boundary condition φi(xj)=δij, δijFor kronecker delta, it is dynamic that part is solved using multi-scale finite element method
Conservation equation is measured, the multiple dimensioned displacement basic functions of each coarse grid cell are obtainedAnd construct multiple dimensioned displacement base
Jacobian matrixWhereinFor small scale displacement basic functions, the small scale displacement basic functions
For carry out small scale fine grid blocks calculate when offset variable selected by basic function,For i-th of fine-structure mesh in coarse grid cell
Displacement discrete value on lattice node,For the quantity of refined net cell node in coarse grid cell.
Multiple dimensioned speed basic function is obtained described in step S4 and the specific steps of multiple dimensioned pressure basic function include:
Two adjacent coarse grid cells are considered as a flow region Ω ', to guarantee local conservativeness, establish the stream
The Part mass conservation equation in dynamic region:
Wherein, b is to finish wet coefficient, and l is time interval l=[0, t], MbTo finish wet modulus,For leading for displacement versus time
Number,For the derivative of Pressure versus Time, λ is fluid flow coefficient, it is preferable that Ω '=Ωi+Ωj, ΩjIndicate multiple dimensioned grid
The region of j-th of coarse grid cell adjacent with i-th of coarse grid cell in system;
Closed boundary condition is chosen, using multi-scale Simulation finite difference method Part mass conservation equation, is obtained every
The corresponding multiple dimensioned pressure basic function of one coarse grid cellIt is corresponding with each coarse grid boundary multiple dimensioned
Speed basic functionAnd construct multiple dimensioned pressure basic function matrixWith more rulers
Spend speed basic function matrixWherein,WithRespectively small scale pressure basic function and small
Scale speed basic function, the small scale pressure basic function are the base of pressure variations when carrying out the calculating of small scale fine grid blocks
Function, the small scale speed basic function are the basic function of speed variables when carrying out the calculating of small scale fine grid blocks,It is thick
The borderline pressure discrete value of i-th of refined net in grid cell,For the borderline speed of i-th of refined net in coarse grid cell
Discrete value is spent,For refined net element number in coarse grid cell;For refined net elementary boundary quantity in coarse net unit.
The specific steps of macroscopical large scale solution that coarse grid cell is obtained described in step S5 include:
By multiple dimensioned displacement basic functions matrix, multiple dimensioned pressure basic function matrix and multiple dimensioned speed basic function matrix assembly
As mapping operator Φ,
I.e.
Construct large scale stiffness matrix A=Φ AfΦT, wherein AfFor small scale stiffness matrix, the small scale stiffness matrix
On refined net, to be formed by stiffness matrix after discrete momentum conservation equation and momentum conservation equation, specific element includes power
Stiffness matrix is learned, finite difference coefficient matrix, and compression term system matrix number is simulated, forms large scale equation group Ax on this basis
=b, wherein b is equation group right-hand vector, including source sink term, displacement boundary conditions and stress boundary condition;
Coupling large scale equation group is carried out to solve the macroscopical large scale solution for obtaining coarse grid cell, including large scale displacement
uc, large scale pressure pc, large scale speed vc。
The specific steps of small scale solution that refined net unit is obtained described in step S6 include:
Construct mapping matrix Nu, Np, Nv, wherein NuIt is the matrix using all multiple dimensioned displacement basic functions as column vector,
The multiple dimensioned displacement basic functions are the corresponding multiple dimensioned displacement basic functions of coarse grid cell each in more grid systems, NpBe with
Matrix of all multiple dimensioned pressure basic functions as column vector, the multiple dimensioned pressure basic function are each thick in more grid systems
The corresponding multiple dimensioned pressure basic function of grid cell, NvIt is the matrix using all multiple dimensioned speed basic functions as column vector, institute
Multiple dimensioned speed basic function is stated as the corresponding multiple dimensioned pressure basic function in coarse grid boundary each in more grid systems;
The corresponding relationship between macroscopical large scale solution and the small scale solution, which is obtained, according to mapping matrix is respectively as follows: uf
=Nuuc, pf≈Ipc+NpDλvc, vf=Nvvc, and macroscopical large scale solution of coarse grid cell is combined, obtain the small of refined net unit
Scale solution, wherein uf, pf, vfSmall scale displacement, small scale pressure and small scale speed are respectively indicated, I is unit matrix, DλFor
Mobility coefficient matrix, the macroscopic view large scale solution are the corresponding solution of coarse grid, and the small scale solution is the corresponding solution of refined net.
Fig. 4 is a kind of another implementation of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid provided by the invention
The flow chart of example, as shown in figure 4, according to oil reservoir actual conditions, obtaining elastic variable in oil reservoir as another embodiment
Property medium seepage flow and stress parameters, establish seepage-pipe coupling model model, and carry out to seepage-pipe coupling model model multiple dimensioned
Mesh generation establishes multiple dimensioned grid system, and then is based on multiple dimensioned grid system, establishes momentum conservation equation and the conservation of mass
Equation, and momentum conservation equation is solved using multi-level finite element modeling hair, it is kept using multi-scale Simulation finite difference method quality
Permanent equation obtains multiple dimensioned basic function, including multiple dimensioned displacement basic functions, multiple dimensioned pressure basic function and multiple dimensioned displacement base letter
Then number solves coarse grid based on multiple dimensioned basic function building stiffness matrix, including coarse grid displacement solves and coarse net
Lattice Pressure solution, according between coarse grid and refined net be displaced and pressure solution mapping relations, carry out refined net displacement solve and
Refined net Pressure solution.
Specific examples are used herein to describe the principles and implementation manners of the present invention, the explanation of above embodiments
Method and its core concept of the invention are merely used to help understand, described embodiment is only that a part of the invention is real
Example is applied, instead of all the embodiments, based on the embodiments of the present invention, those of ordinary skill in the art are not making creation
Property labour under the premise of every other embodiment obtained, shall fall within the protection scope of the present invention.
Claims (4)
1. a kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid, which comprises the steps of:
Reservoir geology parameter and mechanics parameter are obtained, oil reservoir geometrical model is established;
Multiple dimensioned grid dividing is carried out to oil reservoir geometrical model, obtain include coarse net trellis system and fine-structure mesh trellis system more rulers
Spend grid system;
Based on multiple dimensioned grid system, local flow region is chosen, the local momentum conservation equation of the flow region is established, and
Local momentum conservation equation is solved using multi-scale finite element method, obtains multiple dimensioned displacement basic functions, specific steps include:
Each coarse grid cell is considered as a local flow region Ω, establishes the conservation of momentum side in the local flow region
Journey:
▽·(Cdr▽sU-bpI)=0 (1) on Ω × l
Wherein, CdrFor elasticity tensor, l is time interval l=[0, t], and u and p respectively represent motion vector and pressure, and I is unit
Matrix, b are to finish wet coefficient, ▽sFor symmetric gradient operator;
Choose boundary condition φi(xj)=δij, δijFor kronecker delta, local momentum is solved using multi-scale finite element method and is kept
Permanent equation obtains the multiple dimensioned displacement basic functions of each coarse grid cellAnd construct multiple dimensioned displacement basic functions
MatrixWhereinFor small scale displacement basic functions, the small scale displacement basic functions be
Carry out basic function selected by offset variable, φ when refined net calculatingi uFor the position in coarse grid cell on i-th of refined net node
Discrete value is moved,For the quantity of refined net cell node in coarse grid cell;
Based on multiple dimensioned grid system, flow region is chosen, the Part mass conservation equation of the flow region is established, using more
Scale simulation finite difference method Part mass conservation equation, obtains multiple dimensioned speed basic function and multiple dimensioned pressure base letter
Number, specific steps include:
Two adjacent coarse grid cells are considered as a flow region Ω ', to guarantee local conservativeness, establish the flow region
The Part mass conservation equation in domain:
Wherein, b is to finish wet coefficient, and l is time interval l=[0, t], MbTo finish wet modulus,For the derivative of displacement versus time,
For the derivative of Pressure versus Time, λ is fluid flow coefficient, and f is source sink term;
It chooses closed boundary condition and each is obtained using multi-scale Simulation finite difference method Part mass conservation equation
The corresponding multiple dimensioned pressure basic function of coarse grid cellMultiple dimensioned speed corresponding with each coarse grid boundary
Basic functionAnd construct multiple dimensioned pressure basic function matrixWith multiple dimensioned speed
Basic function matrixWherein,WithRespectively small scale pressure basic function and small scale speed
Basic function is spent, the small scale pressure basic function is the basic function of pressure variations when carrying out refined net calculating, the small scale
Speed basic function is the basic function of speed variables when carrying out refined net calculating,For i-th of refined net side in coarse grid cell
Pressure discrete value in boundary, φi vFor the borderline speed dispersion value of i-th of refined net in coarse grid cell,For coarse grid list
Refined net element number in member;For refined net elementary boundary quantity in coarse net unit;
Based on multiple dimensioned principle and multiple dimensioned displacement basic functions, multiple dimensioned speed basic function and multiple dimensioned pressure basic function, obtain
Macroscopical large scale solution of coarse grid cell, and obtain the mapping relations between macroscopical large scale solution and small scale solution;
According to the mapping relations between macroscopical large scale solution and large scale solution and small scale solution, the small scale of refined net unit is obtained
Solution, specific steps include: building mapping matrix Nu, Np, Nv, wherein NuIt is using all multiple dimensioned displacement basic functions as column vector
Matrix, NpIt is the matrix using all multiple dimensioned pressure basic functions as column vector, NvIt is to be made with all multiple dimensioned speed basic functions
For the matrix of column vector;
The corresponding relationship between macroscopical large scale solution and the small scale solution, which is obtained, according to mapping matrix is respectively as follows: uf=
Nuuc, pf≈Ipc+NpDλvc, vf=Nvvc, and macroscopical large scale solution of coarse grid cell is combined, obtain the small ruler of refined net unit
Degree solution, wherein uf, pf, vfRespectively indicate small scale displacement, small scale pressure and small scale speed, ucFor large scale displacement, pcFor
Large scale pressure, vcFor large scale speed, I is unit matrix, DλFor mobility coefficient matrix.
2. a kind of multiple dimensioned flow simulating method of fluid structurecoupling of elastic fluid according to claim 1, which is characterized in that
The specific steps of the multiple dimensioned grid dividing include:
According to survey region size, the large scale coarse grid step-length and quantity of each direction in space are determined, using orthogonal grid pair
Oil reservoir geometrical model carries out small scale fine grid blocks and divides acquisition fine-structure mesh trellis system;
Coarse net trellis system is constructed using load balancing algorithm on the basis of refined net, the coarse net in the coarse net trellis system
Lattice unit is connected with each other by the refined net unit in the fine-structure mesh trellis system;
Multiple dimensioned grid system is formed by the coarse net trellis system and the fine-structure mesh trellis system.
3. a kind of multiple dimensioned flow simulating method of fluid structurecoupling of elastic fluid according to claim 2, which is characterized in that
The fine-structure mesh trellis system includes reservoir rocks, the essential characteristic parameter of fluid and mechanics parameter.
4. a kind of multiple dimensioned flow simulating method of fluid structurecoupling of elastic fluid according to claim 1, which is characterized in that
The specific steps of macroscopical large scale solution for obtaining coarse grid cell include:
By multiple dimensioned displacement basic functions matrix, multiple dimensioned pressure basic function matrix and multiple dimensioned speed basic function matrix assembly become
Mapping operator φ,
I.e.
Construct large scale stiffness matrix A=φ AfφT, wherein AfFor small scale stiffness matrix, the small scale stiffness matrix be
On small scale refined net, stiffness matrix is formed by after discrete momentum conservation equation and momentum conservation equation, specific element includes
Mechanics stiffness matrix simulates finite difference coefficient matrix and compression term system matrix number, forms large scale equation group on this basis
Ax=B, wherein B is equation group right-hand vector, including source sink term, displacement boundary conditions and stress boundary condition;
Coupling large scale equation group is carried out to solve the macroscopical large scale solution for obtaining coarse grid cell, including large scale displacement components uc, greatly
Scale pressure pc, large scale speed vc。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710285521.XA CN106934185B (en) | 2017-04-27 | 2017-04-27 | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710285521.XA CN106934185B (en) | 2017-04-27 | 2017-04-27 | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106934185A CN106934185A (en) | 2017-07-07 |
CN106934185B true CN106934185B (en) | 2019-08-06 |
Family
ID=59437247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710285521.XA Active CN106934185B (en) | 2017-04-27 | 2017-04-27 | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106934185B (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107958124B (en) * | 2017-12-06 | 2021-04-02 | 哈尔滨工业大学深圳研究生院 | Mixed scale coupling method based on stress distribution mode optimization |
CN108829919B (en) * | 2018-05-02 | 2022-03-25 | 嘉兴南湖学院 | Method for multi-scale simulation of multiphase electrohydrodynamic of core-shell structure |
CN110728079B (en) * | 2018-06-26 | 2022-11-04 | 中国石油天然气股份有限公司 | Numerical simulation method and numerical simulation device for flow state of oil-water two-phase flow |
CN109241562B (en) * | 2018-08-02 | 2022-12-16 | 上海交通大学 | Method for measuring elastic property of microstructure material based on multi-scale finite element method |
CN109489745B (en) * | 2018-11-23 | 2020-06-02 | 宁波水表股份有限公司 | Flow measuring method based on data iteration |
US11488357B2 (en) * | 2018-12-20 | 2022-11-01 | Medicim Nv | Automated trimming of a surface mesh |
CN110263434A (en) * | 2019-06-20 | 2019-09-20 | 中国石油大学(华东) | A kind of flow unit method for numerical simulation based on multiple dimensioned mixed finite element |
CN113338218B (en) * | 2021-08-06 | 2021-10-26 | 西南交通大学 | Multi-scale multi-medium comprehensive inversion method for debris flow flexible protection |
CN114564775B (en) * | 2022-02-14 | 2022-10-11 | 中国地质大学(北京) | Double-scale coupling numerical simulation method suitable for building tunnel above goaf |
CN114950888B (en) * | 2022-07-08 | 2023-04-07 | 昆山丘钛光电科技有限公司 | Method, device, equipment and medium for determining pressure of needle cylinder of dispenser |
CN115995279B (en) * | 2023-03-22 | 2023-06-02 | 中国空气动力研究与发展中心计算空气动力研究所 | Material mechanical property evaluation method, device, equipment and readable storage medium |
CN116052820B (en) * | 2023-03-22 | 2023-06-02 | 中国空气动力研究与发展中心计算空气动力研究所 | Material thermal performance evaluation method, device, equipment and readable storage medium |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1673626B1 (en) * | 2003-10-14 | 2011-10-05 | Verseon | Method and apparatus for analysis of molecular combination based on computations of shape complementarity using basis expansions |
CN102224502A (en) * | 2008-10-09 | 2011-10-19 | 雪佛龙美国公司 | Iterative multi-scale method for flow in porous media |
EP1606760B1 (en) * | 2003-03-06 | 2013-01-09 | Chevron U.S.A. Inc. | Multi-scale finite-volume method for use in subsurface flow simulation |
CN104298797A (en) * | 2013-07-16 | 2015-01-21 | 中国石油化工股份有限公司 | Method for determining remaining oil of fracture and hole type oil reservoir under encirclement of high flow guide channel |
CN105631078A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Numerical value simulation method for adaptive medium of natural fractured reservoir |
CN105913494A (en) * | 2016-03-30 | 2016-08-31 | 北京大学 | Multi-scale fracture fine geological modeling and value simulation method and device |
-
2017
- 2017-04-27 CN CN201710285521.XA patent/CN106934185B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1606760B1 (en) * | 2003-03-06 | 2013-01-09 | Chevron U.S.A. Inc. | Multi-scale finite-volume method for use in subsurface flow simulation |
EP1673626B1 (en) * | 2003-10-14 | 2011-10-05 | Verseon | Method and apparatus for analysis of molecular combination based on computations of shape complementarity using basis expansions |
CN102224502A (en) * | 2008-10-09 | 2011-10-19 | 雪佛龙美国公司 | Iterative multi-scale method for flow in porous media |
CN104298797A (en) * | 2013-07-16 | 2015-01-21 | 中国石油化工股份有限公司 | Method for determining remaining oil of fracture and hole type oil reservoir under encirclement of high flow guide channel |
CN105631078A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Numerical value simulation method for adaptive medium of natural fractured reservoir |
CN105913494A (en) * | 2016-03-30 | 2016-08-31 | 北京大学 | Multi-scale fracture fine geological modeling and value simulation method and device |
Non-Patent Citations (4)
Title |
---|
An Efficient Multiscale Mixed Finite Element Method for Modeling Flow in Discrete Fractured Reserviors;Q.Zhang 等;《15th European Conference on the Mathematics of Oil Recovery》;20161231;第1132-1147页 |
Integrated simulation of multi-stage hydraulic fracturing in unconventional reservoirs;Huiying Tang 等;《Journal of Natural Gas Science and Engineering》;20161105;第875-892页 |
基于多尺度混合有限元的离散裂缝两相渗流数值模拟;张庆福 等;《科学通报》;20161104;第62卷(第13期);摘要、正文第1-2节 |
基于离散缝洞网络模型的缝洞型油藏混合多尺度有限元数值模拟;张娜 等;《计算机力学学报》;20150831;第32卷(第4期);摘要、正文第1-4节 |
Also Published As
Publication number | Publication date |
---|---|
CN106934185A (en) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106934185B (en) | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid | |
Chow et al. | Evaluation of turbulence closure models for large-eddy simulation over complex terrain: flow over Askervein Hill | |
Wang et al. | Large deformation finite element analyses in geotechnical engineering | |
CN107060746B (en) | A kind of method of complex fracture oil deposit flow simulation | |
CN106202746B (en) | Simulate the Yeh- multi-level finite element modeling method of Water in Porous Medium stream Darcy velocity | |
Piggott et al. | Anisotropic mesh adaptivity for multi-scale ocean modelling | |
CN108984874A (en) | Obtain the method for numerical simulation in the flow field of potential flows | |
Tootoonchi et al. | A cell-based smoothed point interpolation method for flow-deformation analysis of saturated porous media | |
Van den Berg et al. | Non-linear process based modelling of offshore sand waves | |
Hu et al. | Three-pattern decomposition of global atmospheric circulation: part I—decomposition model and theorems | |
CN107145671B (en) | A kind of numerical reservoir simulation method and system | |
Li et al. | Reduced-order thrust modeling for an efficiently flapping airfoil using system identification method | |
CN110263434A (en) | A kind of flow unit method for numerical simulation based on multiple dimensioned mixed finite element | |
CN112347678B (en) | Novel multi-scale finite element method for simultaneously simulating underground water flow and Darcy speed | |
CN111507026A (en) | Dual-grid multi-scale finite element method for simulating node Darcy permeation flow rate | |
CN106934093A (en) | The Trple grid multi-level finite element modeling method of simulation THREE DIMENSIONAL GROUNDWATER FLOW motion | |
CN109726465A (en) | The three-dimensional method for numerical simulation streamed without viscous low speed based on non-structural curl grid | |
CN103778298A (en) | Improved multi-scale finite element method for stimulating two-dimensional water flow movement in porous media | |
Skamarock et al. | A fully compressible nonhydrostatic deep-atmosphere equations solver for MPAS | |
Berntsen | A perfectly balanced method for estimating the internal pressure gradients in σ-coordinate ocean models | |
CN110083853A (en) | The limited bulk Yeh multi-scale finite element method of simulated groundwater stream movement | |
CN106445882A (en) | Improved CLSVOF method for quickly building signed distance function by VOF function | |
Zhang et al. | A moving mesh finite difference method for non-monotone solutions of non-equilibrium equations in porous media | |
Lamine et al. | Higher order multidimensional upwind convection schemes for flow in porous media on structured and unstructured quadrilateral grids | |
Higdon | Numerical modelling of ocean circulation |
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 |