CN109948171A - The processing method and its corresponding system of reservoir numerical simulation well model - Google Patents
The processing method and its corresponding system of reservoir numerical simulation well model Download PDFInfo
- Publication number
- CN109948171A CN109948171A CN201810711096.0A CN201810711096A CN109948171A CN 109948171 A CN109948171 A CN 109948171A CN 201810711096 A CN201810711096 A CN 201810711096A CN 109948171 A CN109948171 A CN 109948171A
- Authority
- CN
- China
- Prior art keywords
- value
- matrix
- equation
- formula
- module
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 33
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 24
- 230000000452 restraining effect Effects 0.000 claims description 5
- 230000010429 evolutionary process Effects 0.000 abstract description 6
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides the processing method and its corresponding system of a kind of reservoir numerical simulation well model.The processing method, including obtain the original reservoir simulation arrow-shaped sparse matrix；It divides original reservoir modeling arrow-shaped sparse matrix and forms mesh equation matrix, mesh equation correction matrix, well equation matrix and well equation correction matrix, and the matrix that aforementioned division is formed is coupled；Assign the y initial value of matrix equation, and to the mesh equation correction matrix in matrix equation using iterative method solve x value is the first value；By first value substitute into using iterative method solve y value is second value；Judge first value, second value whether and meanwhile the step of restrain.Using the processing method, processing result fast convergence can be made, the evolutionary process of oil field data at any time can be more quickly obtained in conjunction with other relevant parameters.
Description
Technical field
The invention belongs to oilfield exploitation technical fields, and in particular to a kind of processing method of reservoir numerical simulation well model and
Its corresponding system.
Background technique
The data that Need Hierarchy Theory simulation calculates before engineering construction are guided, and it is various that oil extraction process is calculated by simulation
Yield and corresponding intervening measure obtain sequence of operations and data that the maximization of economic benefit of oilfield exploitation needs.
At present common practice be using finite difference calculus by oil reservoir equation time and space it is discrete after to obtain seven diagonal
Linear matrix, well equation then passes through perforation and is coupled with grid, and data are placed on the left side of grid matrix and following, constitute one
A bigger arrowhead-shaped sparse matrix can obtain the evolutionary process of oil field data at any time by solving this sparse matrix.
And well equation data destroys seven diagonal arrangements of mesh equation, may cause master is not dominant diagonally, increases grid side
The conditional number of journey, so that direct solution convergence is relatively slow to be even not achieved convergence precision requirement in given iterative steps, especially
It is that then convergence is less susceptible in the biggish situation of well quantity, to make the acquisition process of the evolutionary process of oil field data at any time
It is especially sluggish.
Summary of the invention
Therefore, the technical problem to be solved in the present invention is that provide a kind of reservoir numerical simulation well model processing method and
Its corresponding system enables the original reservoir simulation arrow-shaped sparse matrix to restrain faster, can in the convergent situation of result
The evolutionary process of oil field data at any time is more quickly obtained in conjunction with other relevant parameters, and then can more quickly be obtained
The sequence of operations and data needed to the maximization of economic benefit of oilfield exploitation.
To solve the above-mentioned problems, the present invention provides a kind of processing method of reservoir numerical simulation well model, including as follows
Step:
A. the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) is obtained；
B. the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) is divided, mesh equation matrix is formed
A (n × n), mesh equation correction matrix R (n × m), well equation correction matrix D (m × n) and well equation matrix H (m × m), and will
The mesh equation matrix A (n × n), mesh equation correction matrix R (n × m) are coupled by formula Ax=b-Ry (1), will
The well equation correction matrix D (m × n) and well equation matrix H (m × m) are coupled by formula Hy=c-Dx (3), wherein
B, c is constant vector；
C. assign y initial value, and substitute into the formula (1), using iterative method solve x value is the first value；
D. first value is substituted into the formula (2), it is second value that direct solution, which obtains y value,；
E. judge first value, second value whether and meanwhile restrain, if convergence, terminate to calculate；If not restraining, by
Initial value of the two-value as y repeats step c, step d, until passing through formula (1), the x value that formula (2) obtains, y value while receiving
It holds back.
Preferably, the y initial value in step c is the solution of 0 or upper one time step.
Preferably, the iterative method is BICGstab.
Preferably, described the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) in step a is by having
The discrete three-dimensional oil reservoir simulation equilibrium equation of calculus of finite differences is limited to obtain.
The present invention also provides a kind of for implementing the processing system of above-mentioned method, including input module, computing module, sentences
Disconnected module, feedback module, the input module be used for acquisition the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) ×
(m+n)), the initial value of y is inputted, and is transmitted in the computing module；The computing module is for receiving the input
The original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) of module input is simultaneously divided, and will be at the beginning of the y
Initial value substitutes into formula (1) the method operation of being iterated and obtains the first value, and the first value substituted into the formula (2) and obtains the
Two-value, the judgment module judge whether to restrain with after second value for obtaining first value to first value, second value,
The feedback module is used to receive receiving simultaneously as a result, working as the judgment module and obtaining the first value with second value for the judgment module
When holding back, then operation is terminated, if the first value and second value do not restrain simultaneously, second value is sent to input module.
The processing method and its corresponding system of reservoir numerical simulation well model provided by the invention make the institute that data volume is huge
It states the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) and splits into four relatively small matrixes of data volume, and
Each matrix is the simple canonical matrix of form, can greatly simplify the calculating process of matrix equation, using iterative method to institute
The formula (1) stated calculate enable unknown number vector x in matrix equation, y restrain faster, which greatly improves meters
Treatment effeciency is calculated, that is, enabling the original reservoir simulation arrow-shaped sparse matrix to restrain faster, the case where result restrains
Under, the evolutionary process of oil field data at any time can be more quickly obtained in conjunction with other relevant parameters, and then can be more
Quickly obtain sequence of operations and data that the maximization of economic benefit of oilfield exploitation needs.
Detailed description of the invention
Fig. 1 is the mesh equation that the middle division of the processing method of the reservoir numerical simulation well model of the embodiment of the present invention is formed
Matrix, mesh equation correction matrix, well equation correction matrix and the signal of well equation matrix-vector product；
Fig. 2 is the step schematic diagram of the processing method of the reservoir numerical simulation well model of the embodiment of the present invention；
Fig. 3 is the module diagram of the processing system of the embodiment of the present invention.
Specific embodiment
In conjunction with referring to Fig. 1 to Fig. 3, according to an embodiment of the invention, providing a kind of processing of reservoir numerical simulation well model
Method includes the following steps:
A. the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) is obtained；
B. the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) is divided, mesh equation matrix is formed
A (n × n), mesh equation correction matrix R (n × m), well equation correction matrix D (m × n) and well equation matrix H (m × m), and will
The mesh equation matrix A (n × n), mesh equation correction matrix R (n × m) are coupled by formula Ax=b-Ry (1), will
The well equation correction matrix D (m × n) and well equation matrix H (m × m) are coupled by formula Hy=c-Dx (2), wherein
B, c is constant vector；
C. assign y initial value, and substitute into the formula (1), using iterative method solve x value is the first value；
D. first value is substituted into the formula (2), it is second value that direct solution, which obtains y value,；
E. judge first value, second value whether and meanwhile restrain, if convergence, terminate to calculate；If not restraining, by
Initial value of the two-value as y repeats step c, step d, until passing through formula (1), the x value that formula (2) obtains, y value while receiving
It holds back.
It can be understood that described the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) be usually it is refined can
Compare matrix.By above-mentioned processing step, make described the original reservoir simulation arrow-shaped sparse matrix A ' ((m+ that data volume is huge
N) × (m+n)) four relatively small matrixes of data volume are split into, and each matrix is the simple canonical matrix of form, energy
The calculating process for enough greatly simplifying matrix equation, calculate to the formula (1) using iterative method makes in matrix equation
Unknown number vector x, y can with restraining faster, and which greatly improves calculation processing efficiency, namely the original reservoir is made to simulate arrow
Capitiform sparse matrix can restrain faster, can be more quickly in conjunction with other relevant parameters in the convergent situation of result
The evolutionary process of oil field data at any time is obtained, and then the maximization of economic benefit that can more quickly obtain oilfield exploitation needs
The sequence of operations and data wanted.
Certainly, the y initial value in step c can be any given, it is preferable that y initial value is the solution of 0 or upper one time step.
The iterative method can be for example BICGstab, this can make x, y fast convergence.
For the acquisition modes of described the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) in step a
For example equilibrium equation can be simulated by the discrete three-dimensional oil reservoir of finite difference calculus and obtained.
The present invention also provides a kind of for implementing the processing system of above-mentioned method, including input module, computing module, sentences
Disconnected module, feedback module, the input module be used for acquisition the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) ×
(m+n)), the initial value of y is inputted, and is transmitted in the computing module；The computing module is for receiving the input
The original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) of module input is simultaneously divided, and will be at the beginning of the y
Initial value substitutes into formula (1) the method operation of being iterated and obtains the first value, and the first value substituted into the formula (2) and obtains the
Two-value, the judgment module judge whether to restrain with after second value for obtaining first value to first value, second value,
The feedback module is used to receive receiving simultaneously as a result, working as the judgment module and obtaining the first value with second value for the judgment module
When holding back, then operation is terminated, if the first value and second value do not restrain simultaneously, second value is sent to input module.
Those skilled in the art will readily recognize that above-mentioned each advantageous manner can be free under the premise of not conflicting
Ground combination, superposition.
The above is merely preferred embodiments of the present invention, be not intended to limit the invention, it is all in spirit of the invention and
Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within principle.Above only
It is the preferred embodiment of the present invention, it is noted that for those skilled in the art, do not departing from this hair
Under the premise of bright technical principle, several improvements and modifications can also be made, these improvements and modifications also should be regarded as guarantor of the invention
Protect range.
Claims (5)
1. a kind of processing method of reservoir numerical simulation well model, which comprises the steps of:
A. the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) is obtained；
B. the original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) is divided, mesh equation matrix A (n is formed
× n), mesh equation correction matrix R (n × m), well equation correction matrix D (m × n) and well equation matrix H (m × m), and by institute
It states mesh equation matrix A (n × n), mesh equation correction matrix R (n × m) to be coupled by formula Ax=b-Ry (1), by institute
Well equation correction matrix D (m × n) and well equation matrix H (m × m) is stated to be coupled by formula Hy=c-Dx (2), wherein b,
C is constant vector；
C. assign y initial value, and substitute into the formula (1), using iterative method solve x value is the first value；
D. first value is substituted into the formula (2), it is second value that direct solution, which obtains y value,；
E. judge first value, second value whether and meanwhile restrain, if convergence, terminate to calculate；If not restraining, by second value
As the initial value of y, step c, step d are repeated, until passing through formula (1), the x value that formula (2) obtains, y value while restraining.
2. the method according to claim 1, wherein the y initial value in step c is the solution of 0 or upper one time step.
3. the method according to claim 1, wherein the iterative method is BICGstab.
4. method according to claim 1, which is characterized in that described the original reservoir in step a simulates the sparse square of arrow-shaped
Battle array A ' ((m+n) × (m+n)) simulates equilibrium equation by the discrete three-dimensional oil reservoir of finite difference calculus and obtains.
5. a kind of for implementing the processing system of method described in any one of any one of claims 1 to 44, including input module, operation
Module, judgment module, feedback module, the input module are used to simulate arrow-shaped sparse matrix A ' to the original reservoir of acquisition
The initial value of ((m+n) × (m+n)), y are inputted, and are transmitted in the computing module；The computing module is for receiving
The original reservoir simulation arrow-shaped sparse matrix A ' ((m+n) × (m+n)) of input module input is simultaneously divided, and by institute
The initial value for stating y substitutes into formula (1) the method operation of being iterated and obtains the first value, and the first value is substituted into the formula (2)
Middle acquisition second value, the judgment module is for judging first value, second value after obtaining first value and second value
Whether restrain, the feedback module be used to receive the judgment module as a result, when the judgment module obtains the first value and the
When two-value restrains simultaneously, then operation is terminated, if the first value and second value do not restrain simultaneously, second value is sent to input mould
Block.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810711096.0A CN109948171A (en) | 2018-07-03 | 2018-07-03 | The processing method and its corresponding system of reservoir numerical simulation well model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810711096.0A CN109948171A (en) | 2018-07-03 | 2018-07-03 | The processing method and its corresponding system of reservoir numerical simulation well model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109948171A true CN109948171A (en) | 2019-06-28 |
Family
ID=67006280
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810711096.0A Pending CN109948171A (en) | 2018-07-03 | 2018-07-03 | The processing method and its corresponding system of reservoir numerical simulation well model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109948171A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210181375A1 (en) * | 2019-12-12 | 2021-06-17 | China Institute Of Water Resources And Hydropower Research | Numerical method for simulating a karez well in association with a groundwater model |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105844009A (en) * | 2016-03-22 | 2016-08-10 | 北京大学 | Efficient sparse matrix storage and numerical reservoir simulation method and apparatus |
| CN105956707A (en) * | 2016-05-09 | 2016-09-21 | 王君岳 | Method of optimizing oil reservoir monitoring system |
| CN107153755A (en) * | 2016-03-03 | 2017-09-12 | 中国石油化工股份有限公司 | A kind of method for solving of shale gas well numerical simulation |
-
2018
- 2018-07-03 CN CN201810711096.0A patent/CN109948171A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107153755A (en) * | 2016-03-03 | 2017-09-12 | 中国石油化工股份有限公司 | A kind of method for solving of shale gas well numerical simulation |
| CN105844009A (en) * | 2016-03-22 | 2016-08-10 | 北京大学 | Efficient sparse matrix storage and numerical reservoir simulation method and apparatus |
| CN105956707A (en) * | 2016-05-09 | 2016-09-21 | 王君岳 | Method of optimizing oil reservoir monitoring system |
Non-Patent Citations (3)
| Title |
|---|
| AYHAM ZAZA 等: "A CUDA based parallel multi-phase oil reservoir simulator", 《COMPUTER PHYSICS COMMUNICATIONS》 * |
| 李中锋等: "黑油模型的预处理正交极小化解法", 《新疆石油地质》 * |
| 李政: "精细油藏数值模拟中的高效求解器研究", 《中国博士学位论文全文数据库 (工程科技Ⅰ辑)》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210181375A1 (en) * | 2019-12-12 | 2021-06-17 | China Institute Of Water Resources And Hydropower Research | Numerical method for simulating a karez well in association with a groundwater model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chahlaoui et al. | Second-order balanced truncation | |
| CN105224984A (en) | A kind of data category recognition methods based on deep neural network and device | |
| Li et al. | Chaotifying linear Elman networks | |
| CN104182209B (en) | A kind of GCRO DR algorithm method for parallel processing based on PETSc | |
| Sezer et al. | Nested epsilon decompositions of linear systems: Weakly coupled and overlapping blocks | |
| CN107069696B (en) | A kind of parallel calculating method of Power system state estimation | |
| CN108897716B (en) | Data processing device and method for reducing calculation amount through memory read-write operation | |
| CN103593703A (en) | Neural network optimizing system and method based on genetic algorithms | |
| Interian et al. | Tumor growth modelling by cellular automata | |
| Düşünceli | Solutions for the Drinfeld-Sokolov equation using an IBSEFM method | |
| CN105629733A (en) | Fractional order cell neural network adaptive synchronization control and circuit design method | |
| Andreassen et al. | Parameter estimation using neural networks in the presence of detector effects | |
| CN109948171A (en) | The processing method and its corresponding system of reservoir numerical simulation well model | |
| RU2016147181A (en) | METHOD AND SYSTEM FOR DETERMINING THE CONFIGURATION OF THE MODEL CONTAINING A COMPLETE OF ELEMENTS AND MEETING A SET OF RESTRICTIONS | |
| Chen et al. | Inverse optimal control of hyperchaotic finance system | |
| CN104809258B (en) | The adaptive amending method of dough sheet normal vector in electromagnetic scattering simulation modeling | |
| Catabriga et al. | Compressible flow SUPG stabilization parameters computed from degree-of-freedom submatrices | |
| Van Gorder et al. | Classification of chaotic regimes in the T system by use of competitive modes | |
| Zhang et al. | Backward Euler-Maruyama method applied to nonlinear hybrid stochastic differential equations with time-variable delay | |
| CN104572588A (en) | Matrix inversion processing method and device | |
| CN109460587B (en) | Finite element calculation method for volcano and seismic visco-elastic deformation automatic modeling | |
| Simes et al. | Introduction to electrical engineering simulation | |
| Dieterich et al. | Observable-targeting global cluster structure optimization | |
| CN104320659A (en) | Background modeling method, device and apparatus | |
| Harris | A simulation method for macro-meteorological wind speeds with a Forward Weibull parent distribution of general index |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190628 |