CN110362842A - Random pore network model modeling method based on various shapes pore throat - Google Patents
Random pore network model modeling method based on various shapes pore throat Download PDFInfo
- Publication number
- CN110362842A CN110362842A CN201810309123.1A CN201810309123A CN110362842A CN 110362842 A CN110362842 A CN 110362842A CN 201810309123 A CN201810309123 A CN 201810309123A CN 110362842 A CN110362842 A CN 110362842A
- Authority
- CN
- China
- Prior art keywords
- network model
- random
- pore
- node
- radius
- 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.)
- Pending
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
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 random pore network model modeling method based on various shapes pore throat that the invention discloses a kind of, this method comprises the following steps: 1) setting the size of random pore network model；2) node coordinate of random pore network model is sought；3) the pore throat shape of random pore network model is selected；4) random pore network model connected probability is set；5) pore throat radius is distributed；6) Ohm's law in random pore network model is established；7) each node voltage of random pore network model is solved；8) random pore network model is generated.Modeling method of the invention considers various shapes pore throat, further to the pore throat type of true core, can deepen the understanding to reservoir rock internal pore structure, further investigation reservoir rock fluid flowing.
Description
Technical field
The present invention relates to the technical fields of reservoir rock random pore network model, and in particular to one kind is based on various shapes
The random pore network model modeling method of pore throat.
Background technique
Pore network model is a kind of computer modeling technique, for stating complicated interstitial space internal structure；And base
The flowing of theoretical and microscopic percolation mechanism Study of Fluid is seeped in exceeding.It is divided into live network model and random pore according to its developing direction
Network model, it is contemplated that live network model previous experiments are at high cost, calculation method is complicated, modeling accuracy is again by experimental facilities
Influence, and random pore network model modeling it is time saving save trouble while having many advantages, such as that model structure is simple, calculation amount is small obtain
It is widely applied.
In the prior art, application No. is 201410805441.9 Chinese invention patents to disclose one kind based on random fractal
Theoretical digital cores and pore network model reconstructing method, digital cores and pore network model based on random fractal theory
Reconstructing method, the capillary pressure curve that conventional mercury injection method is obtained carry out fractal characterization, using random distribution theory deduction with
Machine distribution density function, mean value and variance divide shape expression formula, and random theory and multi-fractal Theory combine the digital rock of building
The heart easily establishes out three-dimensional microcosmic network model accordingly.However, only being considered using the random pore network model established in the past
Simplest round pore throat (only use round tube bank), and the pore throat for finding that true rock has various shapes is studied, such as
Circle, ellipse, cone, hyperbolic side triangle, star etc., it is therefore desirable to establish consider common various shapes pore throat with
Machine pore network model, the random pore network model by establishing various shapes pore throat can be deepened to reservoir rock internal holes
The understanding of gap structure, further investigation reservoir rock fluid flowing.
Summary of the invention
Present invention aims to overcome that the deficiency of above-mentioned background technique, and provide it is a kind of it is applied widely, computational accuracy is high
The random pore network model modeling method based on various shapes pore throat, can establish different shape pore throat, different aspect ratio
Ratio and the scheme of ratio combination.
To achieve the above object, a kind of random pore network model based on various shapes pore throat provided by the present invention is built
Mould method, includes the following steps:
1) size of random pore network model is set:
The length for setting the transverse and longitudinal direction of two-dimension square shape pore network model is L, and number of nodes is n；
2) node coordinate of random pore network model is sought:
For square pore network model, the length in transverse and longitudinal direction, number of nodes are equal, respectively save in pore network model
Point coordinate representation are as follows:
In formula (1), i is the direction x node ID；J is the direction y node ID, and L is the length in transverse and longitudinal direction, and n is node
Number；
3) the pore throat shape of random pore network model is selected:
According to simulation needs, a kind of be used as is chosen from circle, ellipse, cone, hyperbolic side triangle, star and is built
Pore throat shape in vertical random pore network model, and according to the ratio of modeling scheme selection random pore network model aspect ratio
And ratio combination；
4) random pore network model connected probability is set:
A probability function is set in random pore network model program, and adjacent node is determined by random generator
Between restrain and whether be connected to；When there is tube bank in connection, i.e. pore throat radius is randomly assigned by another random function；
5) pore throat radius is distributed:
Random function is used to be randomly assigned tube bank radius to realize that heterogeneous random pore network model, pore throat radius r are logical
Following random function is crossed to generate:
In formula (2), e is the nature truth of a matter；Rand () % is random number；R is pore throat radius, μm；rmaxFor maximum pore throat half
Diameter, μm；rminFor minimum pore throat radius, μm；
6) Ohm's law in random pore network model is established:
By setting the both end voltage of pore network model as V1、V2, network model all-in resistance can be sought by Ohm's law:
In formula (3), R is all-in resistance, Ω；Δ V is voltage, V；I is electric current, A；V1、V2For model both end voltage, V；
7) each node voltage of random pore network model is solved:
According to the conductibility of Kirchoff voltage law, node and line, using the side of following n meshed network model
Journey solves each node voltage of random pore network model:
(q)=AKAT(V) (4)
Wherein
(q)=[q1q2q3...qn]T (5)
(V)=[V1V2V3...Vn]T (6)
In formula, q is total flow, cm3/s；V is total voltage, V；qiFor the volume flow of i-node, (i=1,2,3 ... n), cm3/
s；ViFor the pressure of i-node, (i=1,2,3 ... n), V；K is diagonal matrix (n × n rank)；A is incidence matrix, ATFor turning for matrix A
Set matrix；
8) random pore network model is generated:
Connection situation is restrained according between node coordinate each in network model, adjacent node, random pore network can be constructed
Model framework distributes pore throat shape in random pore network model, in length and breadth in conjunction with the modeling scheme of random pore network model
Than ratio and ratio combination, pore throat radius, the random pore network model for considering various shapes pore throat is generated, and can root
Random pore network model carries out fluid-flow analogy accordingly.
In above-mentioned technical proposal, in the step 3), the ratio of random pore network model aspect ratio is less than 1, random hole
The sum of ratio combination of gap network model aspect ratio is 100%.
In above-mentioned technical proposal, in the step 4), when tube bank, which is in, not to be connected to, distribution pipe beam radius, does not repeat to walk
It is rapid 1)~step 3).
In above-mentioned technical proposal, in the step 5), in pore throat radius value, keeps hydraulic radius value constant, press
Two kinds of distribution modes are uniformly distributed and be uniformly distributed according to logarithm and choose radius, wherein the equally distributed normalization standard deviation of logarithm
σrRespectively 0.05,0.30,0.55,0.80 and 1.05, equally distributed normalization standard deviation are respectively 0.05,0.30 and
0.55。
In above-mentioned technical proposal, in the step 5), when choosing radius according to logarithm even distribution pattern, hydraulic radius
Meet following relationship with normalization standard deviation:
In formula (7-1), rHFor hydraulic radius, μm；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm；
In formula (7-2), σrTo normalize standard deviation；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm.
In above-mentioned technical proposal, in the step 5), when according to even distribution pattern choose radius when, hydraulic radius with return
One change standard deviation meets following relationship:
In formula (8-1), rHFor hydraulic radius, μm；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm；
In formula (8-2), σrTo normalize standard deviation；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm.
In above-mentioned technical proposal, in the step 7), each of random pore network model is solved using successive iteration method
Node voltage:
7.1) initial value is assigned to node:
Being set in flow direction in square mesh is from left to right Far Left boundary node assignment voltage V1, rightmost
Boundary node assignment voltage V2, remaining node valuation is 0；
7.2) equation is established:
For being assigned a value of 0 node in square mesh, following relationship is met according to Kirchoff current law:
In formula (9), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j), S/
m；
Then, relational expression all is listed according to Kirchoff current law to each of network node, obtains equation number
The system of linear equations equal with number of nodes:
In formula (10), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j),
S/m；
Formula (10) is deformed again, is obtained:
In formula (11), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j),
S/m；
7.3) it iteratively solves:
After every progress an iteration solution, V will be obtained1Value and V2Value substitutes into formula (3) and calculates input current in network
With output size of current, when equal or its difference meets the error precision of setting to input current with output electric current, iterative calculation
Terminate.
In above-mentioned technical proposal, in the step 7.3), at calculate node (i, j), left side point (i-1, j) and lower edge point
The voltage of (i, j-1) is replaced using the voltage value that previous step is sought, and meets following relationship:
In formula (12), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j),
S/m。
In above-mentioned technical proposal, in the step 7.3), at calculate node (i, j), left side point (i-1, j) and lower edge point
The voltage of (i, j-1) is replaced using the voltage value that previous step is sought, and is introduced relaxation factor e, is met following relationship:
In formula (13), e indicates relaxation factor, e=1~2；Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For node
The tube bank conductivity of (i, j~i, j), S/m.
Compared with prior art, there are following advantages by the present invention:
First, random pore network model of the invention considers various shapes pore throat, further to true core
Pore throat type can deepen the understanding to reservoir rock internal pore structure, further investigation reservoir rock fluid flowing.
It saves trouble second, the random pore network model modeling established of the present invention is time saving while having that model structure is simple, counts
The advantages that calculation amount is small, overcomes in the prior art that live network model previous experiments are at high cost, calculation method is complicated, modeling accuracy
It is influenced by experimental facilities.
Third, the pore throat character for the more close true core of model that the present invention establishes, therefore on the basis of the model
The various physical characteristics (porosity, permeability, conductive characteristic, thermal property) and fluid neuron network mechanism (hollow billet pressure of study of rocks
Power, relative permeability) it is more accurate, the mistakes and omissions and missing of actual physical experimental data can be corrected and be made up to result, be oil reservoir
Engineer provides reliable reference.
Detailed description of the invention
Fig. 1 is two-dimension square shape pore network model schematic diagram；
Fig. 2 is a node unit of two-dimension square shape pore network model.
Specific embodiment
Below with reference to the embodiment performance that the present invention will be described in detail, but they and do not constitute a limitation of the invention,
It is only for example.Simultaneously by illustrating that advantages of the present invention will become clearer and be readily appreciated that.
As depicted in figs. 1 and 2, the present invention illustrates random pore network model by taking two-dimension square shape pore network model as an example
Modeling, for other types of pore network model, modeling approach is almost the same.
1) size of random pore network model is set
By taking two-dimension square shape pore network model as an example, the length for setting its transverse and longitudinal direction is L, and number of nodes is n；
2) node coordinate of random pore network model is sought
For square pore network model, the length in transverse and longitudinal direction, number of nodes are equal, respectively save in pore network model
Point coordinate may be expressed as:
In formula (1), i is the direction x node ID；J is the direction y node ID, and L is the length in transverse and longitudinal direction, and n is node
Number；
3) the pore throat shape of random pore network model is selected
According to simulation needs, a kind of be used as is chosen from circle, ellipse, cone, hyperbolic side triangle, star and is built
Pore throat shape in vertical random pore network model, and combined according to the ratio and ratio of modeling scheme selection aspect ratio, at random
For the ratio of pore network model aspect ratio less than 1, the sum of ratio combination of random pore network model aspect ratio is 100%.
4) random pore network model connected probability is set.
A probability function is set in random pore network model program, and adjacent node is determined by random generator
Between restrain and whether be connected to；When there is tube bank in connection, i.e. pore throat radius is randomly assigned by another random function；Work as tube bank
In when not being connected to, distribution pipe beam radius, does not repeat step 1)~step 3).
5) pore throat radius is distributed
In view of reservoir rock pore throat radius random distribution, and there is strong heterogeneity, therefore random pore network model
It is middle to use random function to be randomly assigned tube bank radius to realize that heterogeneous random pore network model, pore throat radius r pass through following
Random function generates:
In formula (2), e is the nature truth of a matter；Rand () % is random number；R is pore throat radius, μm；rmaxFor maximum pore throat half
Diameter, μm；rminFor minimum pore throat radius, μm；Wherein pore throat radius assignment, maximum pore throat radius rmaxIt is 41 μm；Minimum pore throat half
Diameter rminIt is 39 μm.
Known by above formula, any pore throat radius is all between maximum pore throat radius to minimum pore throat radius in pore network model
Between (input of minimax pore throat radius program interface), need exist for it is specifically intended that in pore throat radius value, it is necessary to
Keep hydraulic radius value constant.This research is uniformly distributed and is uniformly distributed two kinds of distribution modes according to logarithm and chooses radius,
The wherein equally distributed normalization standard deviation of logarithmrRespectively 0.05,0.30,0.55,0.80 and 1.05, it is equally distributed
Normalizing standard deviation is respectively 0.05,0.30 and 0.55.
Logarithm is uniformly distributed and is uniformly distributed down, and hydraulic radius and normalization standard deviation meet table 1.
Table 1
In above-mentioned table 1, rHFor hydraulic radius, μm；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm；
In formula (7-2), σrTo normalize standard deviation；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm.
6) Ohm's law in random pore network model is established
Since the basis of random pore network model simulation is the water power principle of similitude, analysis fluid flows in porous media
Sunykatuib analysis can be carried out using the network structure being similar in circuit, therefore by setting pore network model both end voltage,
Network model all-in resistance can be sought by Ohm's law:
In formula (3), R is all-in resistance, Ω；Δ V is voltage, V；I is electric current, A；V1、V2For model both end voltage, V；
When carrying out network analog, it is known that network model both end voltage, total resistance value in network model to be asked need to solve
Size of current in network model；And before solving the size of current in network model, need to solve the voltage of each node first.
After obtaining the voltage of all nodes in network model, so that it may calculate the numerical value of network model input/output electric current.
7) each node voltage of random pore network model is solved
Building operation expression be network model simulation core, basic thought be Kirchoff voltage law,
The conductibility of node and line.According to the conductibility of Kirchoff voltage law, node and line, saved using following n
The equation of spot net model solves each node voltage of random pore network model:
(q)=AKAT(V) (4)
Wherein
(q)=[q1q2q3...qn]T (5)
(V)=[V1V2V3...Vn]T (6)
In formula, q is total flow, cm3/s；V is total voltage, V；qiFor the volume flow of i-node, (i=1,2,3 ... n), cm3/
s；ViFor the pressure of i-node, (i=1,2,3 ... n), V；K is diagonal matrix (n × n rank)；A is incidence matrix, ATFor turning for matrix A
Set matrix；
Network model is the analogy method based on Statistics, is known by Statistics, node, session number in network model
Need to reach the sufficiently large reliability that just can guarantee simulation.In this way, for sufficiently large node, line, it cannot be using simple
Algebraic solution.And the written in code for using Cholesky to decompose solution coefficient matrix is verbose, precision will receive rounding error
It influences and runs slowly, for this purpose, using solution by iterative method coefficient matrix.
Iterative method, also known as tossing method, it is corresponding with direct method, it is a kind of method of gradually approaching to reality solution, starting point
It is to carry out subsequent iteration using assumed value as solution to obtain non trivial solution until meeting the condition of convergence.For the convergence system of setting
System, the error solved after iteration each time will be reduced, and solution Xiang Zhenshi solution is also more approached.Iterative method can also automatically correct
The accidental calculating error occurred in iteration, using each node voltage of the solution random pore network model of successive iteration method:
7.1) initial value is assigned to node:
Since last solution is unrelated with given initial value, initial value can arbitrarily be set, but convergent speed in iteration
Degree depends on the order of accuarcy of initial value estimation.It is set in square mesh, flow direction (including electric current flowing and fluid stream
It is dynamic) from left to right, in this way, Far Left boundary node assignment voltage V1, rightmost boundary node assignment voltage V2, the tax of remaining node
Value is 0；
7.2) equation is established:
For being assigned a value of 0 node in square mesh, following relationship is met according to Kirchoff current law:
In formula (9), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j), S/
m；
Then, relational expression all is listed according to Kirchoff current law to each of network node, obtains equation number
The system of linear equations equal with number of nodes:
In formula (10), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j),
S/m；
It is solved to be iterated, formula (10) is deformed, is obtained:
In formula (11), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j),
S/m；
7.3) it iteratively solves:
After every progress an iteration solution, V will be obtained1Value and V2Value substitutes into formula (3) and calculates input current in network
With output size of current, when equal or its difference meets the error precision of setting to input current with output electric current, iterative calculation
Terminate.
In order to accelerate to restrain, the present invention understands the adjacent node voltage value band that previous step is sought in calculate node
Enter to calculate, in two-dimension square shape random pore network model, when calculating point (i, j), left side point (i-1, j) and lower edge point
The voltage value that the voltage of (i, j-1) can be sought with previous step replaces, and meets following relationship:
In formula (12), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j),
S/m。
In order to further speed up convergence, introducing relaxation factor e (between value 1~2) meets following relationship:
In formula (13), e is relaxation factor, e=1~2；Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For node (i, j
~i, j) tube bank conductivity, S/m.
Relaxation factor principle is that give the increment of each point be more than required value when equation being made to reach local equilibrium, thus plus
The convergence rate of speed solution.Relaxation factor e usually has an optimal solution, if selection is appropriate, convergence rate will also be further speeded up.
It calculates repeatedly in this way, until meeting the convergence of given required precision.Network analog of the invention is using super loose
When relaxation solution by iterative method, it is capable of the cube network model of solution node number 100 × 100 × 100.After iteration, find out
Input/output electric current calculates resistance sizes further according to Ohm's law, and then can seek its in random pore network model
His parameter, such as resistivity.
8) random pore network model is generated:
Connection situation is restrained according between node coordinate each in network model, adjacent node, random pore network can be constructed
Model framework distributes pore throat shape in random pore network model, in length and breadth in conjunction with the modeling scheme of random pore network model
Than ratio and ratio combination, pore throat radius, the random pore network model for considering various shapes pore throat is generated, and can root
Random pore network model carries out fluid-flow analogy accordingly.
The content that this specification is not described in detail belongs to the prior art well known to professional and technical personnel in the field.
Claims (9)
1. a kind of random pore network model modeling method based on various shapes pore throat, which comprises the steps of:
1) size of random pore network model is set:
The length for setting the transverse and longitudinal direction of two-dimension square shape pore network model is L, and number of nodes is n；
2) node coordinate of random pore network model is sought:
For square pore network model, the length in transverse and longitudinal direction, number of nodes are equal, and each node is sat in pore network model
Mark indicates are as follows:
In formula (1), i is the direction x node ID；J is the direction y node ID, and L is the length in transverse and longitudinal direction, and n is number of nodes；
3) the pore throat shape of random pore network model is selected:
According to simulation needs, choose from circle, ellipse, cone, hyperbolic side triangle, star a kind of as being established
Pore throat shape in random pore network model, and according to the ratio and ratio of modeling scheme selection random pore network model aspect ratio
Example combination；
4) random pore network model connected probability is set:
A probability function is set in random pore network model program, determines to manage between adjacent node by random generator
Whether beam is connected to；When there is tube bank in connection, i.e. pore throat radius is randomly assigned by another random function；
5) pore throat radius is distributed:
Random function is used to be randomly assigned tube bank radius to realize heterogeneous random pore network model, pore throat radius r is by such as
Lower random function generates:
In formula (2), e is the nature truth of a matter；Rand () % is random number；R is pore throat radius, μm；rmaxFor maximum pore throat radius, μm；
rminFor minimum pore throat radius, μm；
6) Ohm's law in random pore network model is established:
By setting the both end voltage of pore network model as V1、V2, network model all-in resistance can be sought by Ohm's law:
In formula (3), R is all-in resistance, Ω；Δ V is voltage, V；I is electric current, A；V1、V2For model both end voltage, V；
7) each node voltage of random pore network model is solved:
According to the conductibility of Kirchoff voltage law, node and line, using the equation of following n meshed network model,
Solve each node voltage of random pore network model:
(q)=AKAT(V) (4)
Wherein
(q)=[q1q2q3...qn]T (5)
(V)=[V1V2V3...Vn]T (6)
In formula, q is total flow, cm3/s；V is total voltage, V；qiFor the volume flow of i-node, (i=1,2,3 ... n), cm3/s；Vi
For the pressure of i-node, (i=1,2,3 ... n), V；K is diagonal matrix (n × n rank)；A is incidence matrix, ATFor the transposition of matrix A
Matrix；
8) random pore network model is generated:
Connection situation is restrained according between node coordinate each in network model, adjacent node, random pore network model can be constructed
Frame distributes pore throat shape, aspect ratio ratio in random pore network model in conjunction with the modeling scheme of random pore network model
Value and ratio combination, pore throat radius, generate the random pore network model for considering various shapes pore throat, and can be according to this
Random pore network model carries out fluid-flow analogy.
2. the random pore network model modeling method according to claim 1 based on various shapes pore throat, feature exist
In: in the step 3), the ratio of random pore network model aspect ratio is less than 1, the ratio of random pore network model aspect ratio
The sum of example combination is 100%.
3. the random pore network model modeling method according to claim 1 based on various shapes pore throat, feature exist
In: in the step 4), when tube bank, which is in, not to be connected to, distribution pipe beam radius, does not repeat step 1)~step 3).
4. the random pore network model modeling method according to claim 1 based on various shapes pore throat, feature exist
In: in the step 5), in pore throat radius value, keeps hydraulic radius value constant, be uniformly distributed and uniformly according to logarithm
It is distributed two kinds of distribution modes and chooses radius, wherein the equally distributed normalization standard deviation of logarithmrRespectively 0.05,0.30,
0.55,0.80 and 1.05, equally distributed normalization standard deviation is respectively 0.05,0.30 and 0.55.
5. the random pore network model modeling method according to claim 4 based on various shapes pore throat, feature exist
In: in the step 5), when choosing radius according to logarithm even distribution pattern, hydraulic radius and normalization standard deviation meet
Following relationship:
In formula (7-1), rHFor hydraulic radius, μm；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm；Formula (7-
2) in, σrTo normalize standard deviation；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm.
6. the random pore network model modeling method according to claim 4 based on various shapes pore throat, feature exist
In: in the step 5), when choosing radius according to even distribution pattern, hydraulic radius and normalization standard deviation satisfaction are as follows
Relationship:
In formula (8-1), rHFor hydraulic radius, μm；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm；Formula (8-
2) in, σrTo normalize standard deviation；rmaxFor maximum pore throat radius, μm；rminFor minimum pore throat radius, μm.
7. the random pore network model modeling method according to claim 1 based on various shapes pore throat, feature exist
In: in the step 7), using each node voltage of the solution random pore network model of successive iteration method:
7.1) initial value is assigned to node:
Being set in flow direction in square mesh is from left to right Far Left boundary node assignment voltage V1, rightmost boundary section
Point assignment voltage V2, remaining node valuation is 0；
7.2) equation is established:
For being assigned a value of 0 node in square mesh, following relationship is met according to Kirchoff current law:
In formula (9), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j), S/m；
Then, relational expression all is listed according to Kirchoff current law to each of network node, obtains equation number and section
It counts equal system of linear equations:
In formula (10), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j), S/m；
Formula (10) is deformed again, is obtained:
In formula (11), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j), S/m；
7.3) it iteratively solves:
After every progress an iteration solution, V will be obtained1Value and V2Value, which substitutes into formula (3), calculates in network input current and defeated
Size of current out, when equal or its difference meets the error precision of setting to input current with output electric current, iterative calculation terminates.
8. the random pore network model modeling method according to claim 7 based on various shapes pore throat, feature exist
In: in the step 7.3), at calculate node (i, j), the voltage of left side point (i-1, j) and lower edge point (i, j-1) is using upper
The voltage value that one step is sought replaces, and meets following relationship:
In formula (12), Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For the tube bank conductivity of node (i, j~i, j), S/m.
9. the random pore network model modeling method according to claim 7 based on various shapes pore throat, feature exist
In: in the step 7.3), at calculate node (i, j), the voltage of left side point (i-1, j) and lower edge point (i, j-1) is using upper
The voltage value that one step is sought replaces, and introduces relaxation factor e, meets following relationship:
In formula (13), e is relaxation factor, e=1~2；Vi,jFor the voltage of node (i, j), V；g(i, j~i, j)For node (i, j~i,
J) tube bank conductivity, S/m.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810309123.1A CN110362842A (en) | 2018-04-09 | 2018-04-09 | Random pore network model modeling method based on various shapes pore throat |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810309123.1A CN110362842A (en) | 2018-04-09 | 2018-04-09 | Random pore network model modeling method based on various shapes pore throat |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110362842A true CN110362842A (en) | 2019-10-22 |
Family
ID=68213541
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810309123.1A Pending CN110362842A (en) | 2018-04-09 | 2018-04-09 | Random pore network model modeling method based on various shapes pore throat |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110362842A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111950193A (en) * | 2020-07-15 | 2020-11-17 | 中海油田服务股份有限公司 | Modeling method and device of pore network model based on reservoir |
| CN112082917A (en) * | 2020-08-03 | 2020-12-15 | 西南石油大学 | Gas-water unsteady two-phase seepage simulation method based on dynamic network simulation |
| CN112098293A (en) * | 2020-08-03 | 2020-12-18 | 西南石油大学 | Unsteady gas-water two-phase seepage simulation method based on pore fracture dual-medium gas reservoir |
| CN112329358A (en) * | 2020-11-09 | 2021-02-05 | 王立佳 | Method for researching sulfur deposition pore network model of high-sulfur-content gas reservoir |
| CN112796743A (en) * | 2021-01-06 | 2021-05-14 | 中国石油大学(华东) | Core oil accumulation structure generation method and system, computer equipment, terminal and application |
| CN115114787A (en) * | 2022-06-30 | 2022-09-27 | 河南理工大学 | Reservoir fractal pore structure complex texture mode characterization method |
| CN115254213A (en) * | 2022-06-24 | 2022-11-01 | 中国计量大学 | Micro-fluidic chip device based on real soil pore network |
-
2018
- 2018-04-09 CN CN201810309123.1A patent/CN110362842A/en active Pending
Non-Patent Citations (1)
| Title |
|---|
| 陶正武: "低滲砂岩孔隙结构特征值对有效应力变化的响应" * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111950193A (en) * | 2020-07-15 | 2020-11-17 | 中海油田服务股份有限公司 | Modeling method and device of pore network model based on reservoir |
| WO2022011893A1 (en) * | 2020-07-15 | 2022-01-20 | 中海油田服务股份有限公司 | Reservoir-based modeling method and device for pore network model |
| CN111950193B (en) * | 2020-07-15 | 2023-03-14 | 中海油田服务股份有限公司 | Modeling method and device of pore network model based on reservoir |
| CN112082917A (en) * | 2020-08-03 | 2020-12-15 | 西南石油大学 | Gas-water unsteady two-phase seepage simulation method based on dynamic network simulation |
| CN112098293A (en) * | 2020-08-03 | 2020-12-18 | 西南石油大学 | Unsteady gas-water two-phase seepage simulation method based on pore fracture dual-medium gas reservoir |
| CN112329358A (en) * | 2020-11-09 | 2021-02-05 | 王立佳 | Method for researching sulfur deposition pore network model of high-sulfur-content gas reservoir |
| CN112796743A (en) * | 2021-01-06 | 2021-05-14 | 中国石油大学(华东) | Core oil accumulation structure generation method and system, computer equipment, terminal and application |
| CN115254213A (en) * | 2022-06-24 | 2022-11-01 | 中国计量大学 | Micro-fluidic chip device based on real soil pore network |
| CN115254213B (en) * | 2022-06-24 | 2024-05-03 | 中国计量大学 | Microfluidic chip device based on true soil pore network |
| CN115114787A (en) * | 2022-06-30 | 2022-09-27 | 河南理工大学 | Reservoir fractal pore structure complex texture mode characterization method |
| CN115114787B (en) * | 2022-06-30 | 2023-12-01 | 河南理工大学 | Complex organization mode characterization method for reservoir fractal pore structure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110362842A (en) | Random pore network model modeling method based on various shapes pore throat | |
| Gugelmann et al. | Random hyperbolic graphs: degree sequence and clustering | |
| Csányi et al. | Fractal–small-world dichotomy in real-world networks | |
| CN106991244B (en) | A method of Fracture Networks connectivity and seepage calculation based on graph theory | |
| CN110689283B (en) | River system and canal system staggered water supply simulation method and device based on underground water model | |
| Tranter et al. | Pore network modeling of compressed fuel cell components with OpenPNM | |
| CN107680131A (en) | A kind of quick method for determining porous media characterization unit volume size | |
| US20210181375A1 (en) | Numerical method for simulating a karez well in association with a groundwater model | |
| CN105868447B (en) | User communication behavioural analysis and model emulation system based on double-layer network | |
| CN109087301B (en) | Core throat segmentation method based on centering axis and surface model | |
| CN109632604B (en) | Method for coarsening relative permeability of polymer flooding from pore size to core size | |
| CN108880872B (en) | Method and device for decomposing topological structure of Internet test bed | |
| CN104063359A (en) | Implementation method for personalized Chinese character word library | |
| CN104578055A (en) | Power system transient stability simulation method based on network node numbering optimization | |
| Fan et al. | Establishment and verification of the prediction model of soil wetting pattern size in vertical moistube irrigation | |
| CN108510010A (en) | A kind of density peaks clustering method and system based on prescreening | |
| CN104850671B (en) | The grid cutting algorithm in rock core complexity duct | |
| Carlini et al. | Semi-Lagrangian schemes for mean field game models | |
| CN110032818B (en) | Borderless swimming pool water drop flow calculation method | |
| CN111431694B (en) | Multi-scroll chaotic circuit based on sawtooth wave control | |
| CN112613248B (en) | Method and system for improving embedded micro-fluid structure design | |
| CN109472072B (en) | Seasonal river and underground water interaction prediction method based on river simulation | |
| CN113128009A (en) | Sub-basin unit division method considering mountain area plain landform difference | |
| CN112287629A (en) | Impedance network model extraction method and system for very large scale integrated circuit | |
| CN112131803A (en) | Multi-scale pore network multiphase seepage simulation method based on GPU distributed parallel computing |
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 |