CN108804750A - The acquisition methods of trivector permeability suitable for reservoir numerical simulation - Google Patents
The acquisition methods of trivector permeability suitable for reservoir numerical simulation Download PDFInfo
- Publication number
- CN108804750A CN108804750A CN201810356830.6A CN201810356830A CN108804750A CN 108804750 A CN108804750 A CN 108804750A CN 201810356830 A CN201810356830 A CN 201810356830A CN 108804750 A CN108804750 A CN 108804750A
- Authority
- CN
- China
- Prior art keywords
- node unit
- permeability
- axis
- seepage flow
- trivector
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of acquisition methods of the trivector permeability suitable for reservoir numerical simulation, including：Step 1, the seepage area and seepage flow length for calculating adjacent node unit, determine the direction of displacement pressure gradient；Step 2, the included angle cosine value between the normal vector and horizontal X axis, Y-axis and vertical Z axis in adjacent node unit seepage flow face is calculated；Step 3, the included angle cosine value between seepage flow displacement direction and adjacent node unit seepage flow face normal vector is calculated；Step 4, according to the level of node unit, vertical permeability value and vector penetration rate model, the vector penetration value on displacement pressure gradient direction between adjacent node unit is calculated.The acquisition methods of the trivector permeability suitable for reservoir numerical simulation can be quantitatively characterizing reservoir numerical simulation on different displacement pressure gradients and percolation path permeability size, lay the foundation to improve Research Numerical Simulation Techique.
Description
Technical field
The present invention relates to oil reservoir development technical fields, especially relate to a kind of three-dimensional arrow suitable for reservoir numerical simulation
Measure the acquisition methods of permeability.
Background technology
Permeability is an important attribute parameter in Research Numerical Simulation Techique.Permeability has vectorial property, size
With apparent directionality.Development has the crack of different directions in Natural fractures system, and reservoir permeability is made to have centainly
Directionality shows as the directionality of the flowing of fluid；In addition to this, reservoir fracturing improvement is to improve fine and close low-permeability oil deposit to open
The effective means of effect is sent out, and improved reservoir permeability shows as direction characteristic.For such oil reservoir of more efficient exploitation,
To Research Numerical Simulation Techique, more stringent requirements are proposed, predominantly accurately characterizes each crack, and then can be accurate
Describing reservoir fluid flow performance.Research Numerical Simulation Techique based on unstrctured grid technology, which is widely applied to, to be split
In the reservoir numerical simulation of seam property, the grid system used is no longer limited to traditional rectangular mesh.Therefore, space vector permeability
Acquisition be particularly important.Currently, there are no a kind of effective methods for obtaining oil reservoir vector permeability on three dimensions, more
The accurately flowing in description reservoir numerical simulation between grid cell.For this purpose, we have invented one kind being suitable for numerical reservoir mould
The acquisition methods of quasi- trivector permeability, to solve above-mentioned problem.
Invention content
It can be quantified the object of the present invention is to provide one kind and calculate space three-dimensional vector permeability, to improve numerical reservoir mould
The acquisition methods of the quasi- trivector permeability suitable for reservoir numerical simulation that technical support is provided.
The purpose of the present invention can be achieved by the following technical measures：Trivector suitable for reservoir numerical simulation is permeated
The acquisition methods of rate, the acquisition methods for being somebody's turn to do the trivector permeability suitable for reservoir numerical simulation include：Step 1, phase is calculated
The seepage area and seepage flow length of neighbors unit, determine the direction of displacement pressure gradient；Step 2, adjacent node unit is calculated
The normal vector and horizontal X axis in seepage flow face, the included angle cosine value between Y-axis and vertical Z axis；Step 3, calculate seepage flow displacement direction with
Included angle cosine value between the normal vector of adjacent node unit seepage flow face；Step 4, according to the level of node unit, vertical permeability
Value and vector penetration rate model calculate the vector penetration value on displacement pressure gradient direction between adjacent node unit.
The purpose of the present invention can be also achieved by the following technical measures：
Step 1 includes：
A, node unit I and node unit J are spatially adjacent two node units, SABCDIt is node unit I and section
The shared seepage area of dot element J, the body centre coordinate I (xi, yi, zi) of calculate node unit I；
B calculates the space coordinate O (x of the seepage flow face faces ABCD center Oo,yo,zo)；
C calculates the flowing distance LIO, and according to the flowing distance LIOComputational methods calculate the flowing distance LJO。
In step a, spatial point A, B, C, D, E, F, H, G are the space coordinate point of configuration node unit I, node unit I
The calculation formula of body centre coordinate xi be
The value of computational methods coordinates computed point yi and zi according to coordinate points x.
In stepb, space coordinate O (xo,yo,zo) computational methods with node unit I body centre coordinate xi calculating
Method.
In step c, the flowing distance L is calculatedIOUsed calculation formula is：
In step 1, Δ PiFor the pressure of node unit I, Δ PjFor the pressure of node unit J, as Δ Pi>ΔPjWhen,
Displacement pressure gradient direction is from node unit I to node unit J, as Δ Pi<ΔPjWhen, displacement pressure gradient direction be from
Node unit J to node unit I.
In step 2, the normal vector in the seepage flow face of calculate node unit I and J For oozing for adjacent node unit I and J
The normal vector of stream interface △ ABC calculates normal vector using the cross multiplication of vectorCalculation formula isSeepage flow face △
The normal vector of ABCSpace coordinate be N (xn,yn,zn)。
In step 2, a point coordinates X1 (1,0,0) is taken in reference axis X-axis, calculates the normal vector of seepage flow face △ ABCWith
Included angle cosine between X-axis：
A point coordinates Y1 (0,1,0) is taken in reference axis Y-axis, calculates the normal vector of seepage flow face △ ABCFolder between Y-axis
Angle cosine：
A point coordinates X1 (0,0,1) is taken on reference axis Z axis, calculates the normal vector of seepage flow face △ ABCFolder between Z axis
Angle cosine：
In step 3, vectorialIt is node unit I along the vector of seepage flow displacement gradient direction, concrete form isThe normal vector of vector and seepage flow face △ ABC on the seepage flow displacement directionBetween
The calculation formula of included angle cosine value is
In step 4, permeability Ks of the acquisition node unit I in X-direction, Y-direction and Z-directionxi、Kyi、Kzi；Obtain node
Permeability Ks of the unit J in X-direction, Y-direction and Z-directionxj、Kyj、Kzj；Calculate node unit I and adjacent node unit J is in displacement
Trivector permeability on direction.
In step 4, as Δ Pi>ΔPjWhen, displacement pressure gradient direction is from node unit I to adjacent node unit J
Trivector permeability K on displacement directionpi, calculation formula is
Kpi=(Kxicos2α+K yicos2β+Kzicos2ψ)cosθ (8)
As Δ Pi<ΔPjWhen, the trivector from node unit J to adjacent node unit I along displacement pressure gradient direction is oozed
Saturating rate Kpj, calculation formula is
Kpj=(Kxjcos2α+Kyjcos2β+Kzjcos2ψ)cosθ (9)
In formula, Kpi, KpjFor permeabilities of the two neighboring node unit I and J on displacement pressure gradient direction, md；Kxi、
KxjTo be respectively the permeability of node unit I and J horizontal direction X-axis, md；Kyi、KyjRespectively node unit I and J are in level side
To the permeability of Y-axis, md；Kzj、KzjRespectively node unit I and J is in the permeability of vertical direction Z axis, md；α, β, ψ are respectively
Angle between the normal vector and horizontal X axis, Y-axis and vertical Z axis in the shared seepage flow face of node unit I and adjacent node unit J,
Radian；Angles of the θ between displacement pressure gradient direction and the seepage flow face normal vector.
The acquisition methods of the trivector permeability suitable for reservoir numerical simulation in the present invention are research oil, day
The important method of the seepage flow characteristics and Research Numerical Simulation Techique of the underground fluids in LOW PERMEABILITY POROUS MEDIA such as right gas.With multiple
Miscellaneous medium oil pool developmental research deepens continuously and deepens, and the research of multi-dielectric Research Numerical Simulation Techique is also grown with applying
The development of foot.Permeability is an important attribute parameter in Research Numerical Simulation Techique.Permeability has vectorial property, size
With apparent directionality.In reservoir numerical simulation, the directional difference of permeability directly affects conductivity between adjacent node
Size, and then influence the propagation of pressure and the distribution of saturation degree under Model Condition.This method is directed in reservoir numerical simulation
The problem of vector permeability acquisition methods, in the base for obtaining each node unit level of reservoir numerical simulation, vertical permeability numerical value
On plinth, using vector calculation model of permeability, the penetration value size on different seepage flow faces is calculated.This method establishes can
The method for quantitatively calculating vector permeability in reservoir numerical simulation, lays the foundation to improve Research Numerical Simulation Techique.
Description of the drawings
Fig. 1 is a specific implementation of the acquisition methods of the trivector permeability suitable for reservoir numerical simulation of the present invention
The flow chart of example；
Fig. 2 is neighboring mesh nodes unit flow path unicom figure in the specific embodiment of the present invention；
Fig. 3 be the present invention a specific embodiment in seepage flow face, seepage direction spatial relation graph.
Specific implementation mode
For enable the present invention above and other objects, features and advantages be clearer and more comprehensible, it is cited below particularly go out preferable implementation
Example, and coordinate shown in attached drawing, it is described in detail below.
As shown in FIG. 1, FIG. 1 is the acquisition methods of the trivector permeability suitable for reservoir numerical simulation of the present invention
Flow chart.
In a step 101, the seepage area and seepage flow length for calculating adjacent node unit, determine the side of displacement pressure gradient
To.The determination method in the direction of seepage area, seepage flow length and displacement pressure gradient is：Node unit I and node unit J is empty
Between upper two adjacent node units, as shown in Figure 2.
SABCDIt is the shared seepage area of node unit I and node unit J, which can be according to the specific side in seepage flow face
Number is divided into several triangles, calculates separately the area of each triangle, and all triangle areas are added up and are obtained
Seepage flow face, the ABCD of seepage flow face shown in Fig. 2 are a quadrangle；Two triangles can be divided into, calculate separately two three
Angular area is to obtain the area of seepage flow face ABCD；
LIO、LJOIt is the distance of node unit I and node unit J to the face center O of seepage flow face ABCD respectively.With seepage flow away from
From LIOIts computational methods is illustrated for example：The body centre coordinate I (xi, yi, zi) of calculate node unit I first, spatial point
A, B, C, D, E, F, H, G are the space coordinate point of configuration node unit I, the calculation formula of the body centre coordinate xi of node unit I
For
The value of computational methods coordinates computed point yi and zi according to coordinate points x；Secondly, the seepage flow face faces ABCD center O is calculated
Space coordinate O (xo,yo,zo), computational methods of the computational methods with the body centre coordinate xi of node unit I；Finally, it calculates
The flowing distance LIO, used calculation formula is：
According to the flowing distance LIOComputational methods calculate the flowing distance LJO；
ΔPiFor the pressure of node unit I, Δ PjFor the pressure of node unit J, as Δ Pi>ΔPjWhen, displacement pressure ladder
It is from node unit I to node unit J, as Δ P to spend directioni<ΔPjWhen, displacement pressure gradient direction be from node unit J to
Node unit I.
In a step 102, between the normal vector and horizontal X axis, Y-axis and vertical Z axis that calculate adjacent node unit seepage flow face
Included angle cosine value；
The normal vector in the seepage flow face of calculate node unit I and J firstAs shown in figure 3,For adjacent node unit I and J
Seepage flow face △ ABC normal vector, using vector cross multiplication calculate normal vectorCalculation formula isSeepage flow
The normal vector of face △ ABCSpace coordinate be N (xn,yn,zn)。
A point coordinates X1 (1,0,0) is taken in reference axis X-axis, calculates the normal vector of seepage flow face △ ABCFolder between X-axis
Angle cosine：
A point coordinates Y1 (0,1,0) is taken in reference axis Y-axis, calculates the normal vector of seepage flow face △ ABCFolder between Y-axis
Angle cosine：
A point coordinates X1 (0,0,1) is taken on reference axis Z axis, calculates the normal vector of seepage flow face △ ABCFolder between Z axis
Angle cosine：
In step 103, it calculates between seepage flow displacement pressure gradient direction and adjacent node unit seepage flow face normal vector
Included angle cosine value；
As shown in figure 3, vectorialIt is node unit I along the vector of seepage flow displacement gradient direction, concrete form isThe normal vector of vector and seepage flow face △ ABC on the seepage flow displacement directionBetween
The calculation formula of included angle cosine value is
At step 104, permeability Ks of the acquisition node unit I in X-direction, Y-direction and Z-direction firstxi、Kyi、Kzi；It obtains
Take node unit J in the permeability K of X-direction, Y-direction and Z-directionxj、Kyj、Kzj；
Calculate node unit I and trivector permeabilities of the adjacent node unit J on displacement direction.As Δ Pi>ΔPj
When, displacement pressure gradient direction is the trivector permeability from node unit I to adjacent node unit J on displacement direction
Kpi, calculation formula is
Kpi=(Kxicos2α+Kyicos2β+Kzicos2ψ)cosθ (8)
As Δ Pi<ΔPjWhen, the trivector from node unit J to adjacent node unit I along displacement pressure gradient direction is oozed
Saturating rate Kpj, calculation formula is
Kpj=(Kxjcos2α+Kyjcos2β+Kzjcos2ψ)cosθ (9)
In formula, Kpi, KpjFor permeabilities of the two neighboring node unit I and J on displacement pressure gradient direction, md；Kxi、
KxjTo be respectively the permeability of node unit I and J horizontal direction X-axis, md；Kyi、KyjRespectively node unit I and J are in level side
To the permeability of Y-axis, md；Kzj、KzjRespectively node unit I and J is in the permeability of vertical direction Z axis, md；α, β, ψ are respectively
Angle between the normal vector and horizontal X axis, Y-axis and vertical Z axis in the shared seepage flow face of node unit I and adjacent node unit J,
Radian；Angles of the θ between displacement pressure gradient direction and the seepage flow face normal vector.Flow terminates.
Claims (11)
1. the acquisition methods of the trivector permeability suitable for reservoir numerical simulation, which is characterized in that oil reservoir number should be suitable for
The acquisition methods of trivector permeability of value simulation include：
Step 1, the seepage area and seepage flow length for calculating adjacent node unit, determine the direction of displacement pressure gradient；
Step 2, the included angle cosine between the normal vector and horizontal X axis, Y-axis and vertical Z axis in adjacent node unit seepage flow face is calculated
Value；
Step 3, the included angle cosine value between seepage flow displacement direction and adjacent node unit seepage flow face normal vector is calculated；
Step 4, it according to the level of node unit, vertical permeability value and vector penetration rate model, calculates between adjacent node unit
Vector penetration value on displacement pressure gradient direction.
2. the acquisition methods of the trivector permeability according to claim 1 suitable for reservoir numerical simulation, feature
It is, step 1 includes：
A, node unit I and node unit J are spatially adjacent two node units, SABCDIt is node unit I and node unit
The shared seepage area of J, the body centre coordinate I (xi, yi, zi) of calculate node unit I；
B calculates the space coordinate O (x of the seepage flow face faces ABCD center Oo,yo,zo)；
C calculates the flowing distance LIO, and according to the flowing distance LIOComputational methods calculate the flowing distance LJO。
3. the acquisition methods of the trivector permeability according to claim 2 suitable for reservoir numerical simulation, feature
It is, in step a, spatial point A, B, C, D, E, F, H, G are the space coordinate point of configuration node unit I, the body of node unit I
The calculation formula of centre coordinate xi is
The value of computational methods coordinates computed point yi and zi according to coordinate points x.
4. the acquisition methods of the trivector permeability according to claim 3 suitable for reservoir numerical simulation, feature
It is, in stepb, space coordinate O (xo,yo,zo) computational methods with node unit I body centre coordinate xi calculating side
Method.
5. the acquisition methods of the trivector permeability according to claim 2 suitable for reservoir numerical simulation, feature
It is, in step c, calculates the flowing distance LIOUsed calculation formula is：
6. the acquisition methods of the trivector permeability according to claim 2 suitable for reservoir numerical simulation, feature
It is, in step 1, Δ PiFor the pressure of node unit I, Δ PjFor the pressure of node unit J, as Δ Pi>ΔPjWhen, it drives
It is from node unit I to node unit J, as Δ P for barometric gradient directioni<ΔPjWhen, displacement pressure gradient direction is from section
Dot element J to node unit I.
7. the acquisition methods of the trivector permeability according to claim 2 suitable for reservoir numerical simulation, feature
It is, in step 2, the normal vector in the seepage flow face of calculate node unit I and JFor the seepage flow of adjacent node unit I and J
The normal vector of face △ ABC calculates normal vector using the cross multiplication of vectorCalculation formula isSeepage flow face △ ABC
Normal vectorSpace coordinate be N (xn,yn,zn)。
8. the acquisition methods of the trivector permeability according to claim 7 suitable for reservoir numerical simulation, feature
It is, in step 2, a point coordinates X1 (1,0,0) is taken in reference axis X-axis, calculate the normal vector of seepage flow face △ ABCWith X-axis
Between included angle cosine：
A point coordinates Y1 (0,1,0) is taken in reference axis Y-axis, calculates the normal vector of seepage flow face △ ABCMore than angle between Y-axis
String：
A point coordinates X1 (0,0,1) is taken on reference axis Z axis, calculates the normal vector of seepage flow face △ ABCMore than angle between Z axis
String：
9. the acquisition methods of the trivector permeability according to claim 8 suitable for reservoir numerical simulation, feature
It is, in step 3, vectorIt is node unit I along the vector of seepage flow displacement gradient direction, concrete form isThe normal vector of vector and seepage flow face △ ABC on the seepage flow displacement directionBetween
The calculation formula of included angle cosine value is
10. the acquisition methods of the trivector permeability according to claim 9 suitable for reservoir numerical simulation, feature
It is, in step 4, permeability Ks of the acquisition node unit I in X-direction, Y-direction and Z-directionxi、Kyi、Kzi；Obtain node list
Permeability Ks of first J in X-direction, Y-direction and Z-directionxj、Kyj、Kzj；Calculate node unit I and adjacent node unit J is in displacement side
Upward trivector permeability.
11. the acquisition methods of the trivector permeability according to claim 10 suitable for reservoir numerical simulation, special
Sign is, in step 4, as Δ Pi>ΔPjWhen, displacement pressure gradient direction is from node unit I to adjacent node unit J
Trivector permeability K on displacement directionpi, calculation formula is
Kpi=(Kxicos2α+Kyicos2β+Kzicos2ψ)cosθ (8)
As Δ Pi<ΔPjWhen, from node unit J to adjacent node unit I along the trivector permeability in displacement pressure gradient direction
Kpj, calculation formula is
Kpj=(Kxjcos2α+Kyjcos2β+Kzjcos2ψ)cosθ (9)
In formula, Kpi, KpjFor permeabilities of the two neighboring node unit I and J on displacement pressure gradient direction, md；Kxi、KxjFor
The permeability of respectively node unit I and J horizontal direction X-axis, md；Kyi、KyjRespectively node unit I and J Y in the horizontal direction
The permeability of axis, md；Kzj、KzjRespectively node unit I and J is in the permeability of vertical direction Z axis, md；α, β, ψ are respectively to save
Angle between the normal vector and horizontal X axis, Y-axis and vertical Z axis in the shared seepage flow face of dot element I and adjacent node unit J, arc
Degree；Angles of the θ between displacement pressure gradient direction and the seepage flow face normal vector.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810356830.6A CN108804750B (en) | 2018-04-19 | 2018-04-19 | Three-dimensional vector permeability obtaining method suitable for numerical reservoir simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810356830.6A CN108804750B (en) | 2018-04-19 | 2018-04-19 | Three-dimensional vector permeability obtaining method suitable for numerical reservoir simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108804750A true CN108804750A (en) | 2018-11-13 |
| CN108804750B CN108804750B (en) | 2021-11-12 |
Family
ID=64092986
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810356830.6A Active CN108804750B (en) | 2018-04-19 | 2018-04-19 | Three-dimensional vector permeability obtaining method suitable for numerical reservoir simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108804750B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6023656A (en) * | 1996-12-30 | 2000-02-08 | Institut Francais Du Petrole | Method for determining the equivalent fracture permeability of a fracture network in a subsurface multi-layered medium |
| CN104102802A (en) * | 2013-04-03 | 2014-10-15 | 中国石油化工股份有限公司 | Method for representing oil-water relative permeability curve under condition of different driving displacement pressure gradients |
| CN105510203A (en) * | 2015-11-18 | 2016-04-20 | 中国石油新疆油田分公司勘探开发研究院 | Method for determination of sandstone oil reservoir oil-water relative permeability under different temperature gradients |
| CN106979918A (en) * | 2017-04-20 | 2017-07-25 | 中国石油大学(北京) | A kind of method and device for the Test Liquid Permeability of Core for obtaining compact oil reservoir rock core |
| CN107346518A (en) * | 2016-05-06 | 2017-11-14 | 中国石油化工股份有限公司 | The acquisition methods of fine and close low-permeability oil deposit oil-water two-phase flow maximum filtrational resistance gradient |
-
2018
- 2018-04-19 CN CN201810356830.6A patent/CN108804750B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6023656A (en) * | 1996-12-30 | 2000-02-08 | Institut Francais Du Petrole | Method for determining the equivalent fracture permeability of a fracture network in a subsurface multi-layered medium |
| CN104102802A (en) * | 2013-04-03 | 2014-10-15 | 中国石油化工股份有限公司 | Method for representing oil-water relative permeability curve under condition of different driving displacement pressure gradients |
| CN105510203A (en) * | 2015-11-18 | 2016-04-20 | 中国石油新疆油田分公司勘探开发研究院 | Method for determination of sandstone oil reservoir oil-water relative permeability under different temperature gradients |
| CN107346518A (en) * | 2016-05-06 | 2017-11-14 | 中国石油化工股份有限公司 | The acquisition methods of fine and close low-permeability oil deposit oil-water two-phase flow maximum filtrational resistance gradient |
| CN106979918A (en) * | 2017-04-20 | 2017-07-25 | 中国石油大学(北京) | A kind of method and device for the Test Liquid Permeability of Core for obtaining compact oil reservoir rock core |
Non-Patent Citations (3)
| Title |
|---|
| 周涌沂 等: "渗透率的矢量性研究", 《新疆石油地质》 * |
| 孙业恒: "史南油田史深100块裂缝性砂岩油藏建模及数值模拟研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
| 苏海波 等: "低渗透储层水驱油渗流阻力特征", 《油气地质与采收率》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108804750B (en) | 2021-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107060746B (en) | A kind of method of complex fracture oil deposit flow simulation | |
| CN104533370B (en) | Pressure break horizontal well oil reservoir, crack, pit shaft coupled model method | |
| Ren et al. | Pipe network model for unconfined seepage analysis in fractured rock masses | |
| CN105178939B (en) | A kind of prediction technique for channel pressure break flow conductivity | |
| CN105574320B (en) | The evaluation method of the effective percolation ability of low permeability sandstone reservoirs | |
| AU2012375233B2 (en) | System and method for automatic local grid refinement in reservoir simulation systems | |
| CN106649963B (en) | Volume fracturing complexity seam net average crack length and equivalent fissure item number determine method | |
| CN109977469A (en) | A kind of Two-dimensional Porous Medium model building method based on Voronoi diagram | |
| CN103907118A (en) | System and method for coarsening in reservoir simulation system | |
| CN110738001B (en) | Unconventional reservoir fracturing yield-increasing transformation area calculation method | |
| Aleroev et al. | Features of seepage of a liquid to a chink in the cracked deformable layer | |
| CN105808862A (en) | Displacement analysis method for determining slope critical sliding surface | |
| CN113836695B (en) | Oil reservoir numerical simulation method based on gridless connecting element | |
| CN102288986A (en) | Method for obtaining elastic modulus of carbonate rock reservoir under earthquake scale | |
| CN103455667A (en) | Numerical simulation method for controlling confined aquifer seawater invasion through inflation process | |
| CN115114834B (en) | Fracturing well test simulation method under complex condition | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| CN111125905B (en) | Two-dimensional fracture network expansion model for coupling oil reservoir fluid flow and simulation method thereof | |
| Wei et al. | Equivalent fracture network model for steady seepage problems with free surfaces | |
| CN107357759B (en) | Seepage solving method based on seepage boundary and motion differential equation condition | |
| CN108804750A (en) | The acquisition methods of trivector permeability suitable for reservoir numerical simulation | |
| Cao et al. | Unsaturated seepage analysis of cracked soil including development process of cracks | |
| CN109558614A (en) | The analogy method and system that gas flows in shale gas reservoir multi-scale facture | |
| CN105160088B (en) | The computational methods of the rate of discharge of groundwater | |
| Gioda et al. | Some numerical techniques for free-surface seepage analysis |
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 |