CN104392109B - A kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir - Google Patents
A kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir Download PDFInfo
- Publication number
- CN104392109B CN104392109B CN201410612927.0A CN201410612927A CN104392109B CN 104392109 B CN104392109 B CN 104392109B CN 201410612927 A CN201410612927 A CN 201410612927A CN 104392109 B CN104392109 B CN 104392109B
- Authority
- CN
- China
- Prior art keywords
- crack
- grid
- oil reservoir
- matrix
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the present application discloses a kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir, and this method includes：Mesh generation is carried out to oil reservoir, the grid of oil reservoir is obtained；Multiple secondary grid is further divided into the grid of oil reservoir；The conductivity between matrix and matrix, matrix and crack, crack and crack is calculated, and calculates the fluid communication coefficient of reservoir and pit shaft, and further obtains pressure and saturation degree under each time step；Cracking initiation and the criterion of dynamic extension in the injecting process, and the coupling relation set up between intrinsic fracture and water filling dynamic crack, pressure-break are set up, the fracture extension situation of different time points is judged；The dynamic grid division by boundary condition of the new crack of extension, while the attribute information in the new crack is assigned into net boundary, to update the net boundary；Oil reservoir and rock mechanics parameters are corrected, historical data fitting is carried out, the saturation degree and pressure condition of current time step is calculated.
Description
Technical field
The application is related to Research Numerical Simulation Techique field, more particularly to a kind of oil reservoir dynamic discrete results of fracture simulation method
And device.
Background technology
In flow through oil reservoir situation calculating field, results of fracture simulation method is widely used.Results of fracture simulation of the prior art
Method mainly has two kinds, i.e. diplopore are double to blend discrete fractures analogy method.
The double infiltration methods of diplopore, come the coupling of simulation fracture and matrix, crack are done using equivalent continuum Model
The processing of homogenization.In the model, Rock Matrix space provides the main memory space of fluid, and main flow of fluid is then
Occur between crack and crack.Fluid communication between crack and matrix is retouched using " transmission equation " or " form factor "
State.Although this model computational efficiency is improved, model is excessively idealized, and it has many limitations, it is impossible to accurately portray oil reservoir
Crack state and flow condition, when fissured structure is complicated, then need large-scale grid.Additionally, due to negligible amounts, across
Region is wider, and is not suitable with very much to flowing determinative results of fracture simulation, causes the transmission equation in matrix and crack to be difficult to
It is accurate to calculate.
Discrete fractures analogy method carries out in all cracks according to its actual size and distributional pattern complete and dominant
Description, the flowing of crack and matrix is all by Darcy's law control.In the analogy method, Fractured Reservoir system passes through high-resolution
Grid present, therefore discrete fractures analogy methods has the accuracy of certain applicability and height.But discrete fractures
The shortcoming of analogy method is amount of calculation than larger.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir, can realize
Not only accurate but also fast and effectively simulation fracture.
In order to solve the above technical problems, on the one hand the application provides a kind of dynamic discrete results of fracture simulation method of oil reservoir, should
Method includes：
Mesh generation is carried out to oil reservoir, the grid of oil reservoir is obtained；
Multiple secondary grid is further divided into the grid of the oil reservoir；
Grid and multiple secondary grid based on the oil reservoir, calculate matrix and matrix, matrix and crack, crack and crack
Between conductivity, and calculate the fluid communication coefficient of reservoir and pit shaft, and further obtain the pressure under each time step
Power and saturation degree；
Cracking initiation and the criterion of dynamic extension in the injecting process are set up, and sets up intrinsic fracture and is dynamically split with water filling
Coupling relation between seam, pressure-break, judges the fracture extension situation of different time points；
The dynamic grid division by boundary condition of the new crack of extension, while the attribute information in the new crack is assigned into net
Lattice border, to update the net boundary；
Pressure field and saturation field are recalculated with the net boundary after the renewal；
Oil reservoir and rock mechanics parameters are corrected, historical data fitting is carried out, the saturation degree and pressure of current time step is calculated
Situation.
In a preferred embodiment, it is described that oil reservoir progress mesh generation is included：
Use the Delaunary subdivision methods using fracture surface as geometrical constraint to carry out mesh generation to the oil reservoir, will split
Seam is processed as curved surface, is tetrahedral grid by matrix treatments.
In a preferred embodiment, it is described that oil reservoir progress mesh generation is also included：
Matrix mesh volume is subtracted to the half of the volume of adjacent crack grid.
In a preferred embodiment, methods described also includes：
Judge whether new crack occur.
In a preferred embodiment, it is described to judge whether that new crack occur includes：
Using the sophisticated coordinate of pre-existing fracture as the starting point coordinate in new crack, judge to be connected with new crack in each time step
The pressure of the sophisticated grid of pre-existing fracture whether be more than its critical pressure, if more than its critical pressure, being judged as occurring newly splitting
Seam.
In a preferred embodiment, the attribute information in the new crack, which is assigned to net boundary, includes：
If new crack produce after the crack grid element center of system and the crack grid center superposition of original system or
On the crack grid of original system, then by the grid pressure and saturation degree indirect assignment of original system to the system after the generation of new crack
On the grid of crack；
If the crack grid element center of the system after new crack generation is in the matrix grid of original system, by original system
On the crack grid of system after crack grid pressure and saturation degree assignment to new crack generation；
If the matrix grid element center of system after new crack is produced is overlapped or in original with the matrix grid element center of original system
In the matrix grid of system, then after producing the pressure and saturation degree indirect assignment of the matrix grid of original system to new crack is
On the matrix grid of system.
On the other hand the application also provides a kind of dynamic discrete board slit simulation apparatus of oil reservoir, including：
Oil reservoir mesh generation unit, for carrying out mesh generation to oil reservoir, obtains the grid of oil reservoir；
Secondary mesh generation unit, multiple secondary grid is further divided into for the grid to the oil reservoir；
Fluid exchange systems computing unit, for the grid based on the oil reservoir and multiple secondary grid, calculate matrix with
Conductivity between matrix, matrix and crack, crack and crack, and the fluid communication coefficient of reservoir and pit shaft is calculated, go forward side by side
One step obtains pressure and saturation degree under each time step；
Fracture extension situation judging unit, for setting up cracking initiation and the criterion of dynamic extension in the injecting process,
And the coupling relation set up between intrinsic fracture and water filling dynamic crack, pressure-break, judge that the crack of different time points is prolonged
Stretch situation；
Attribute information assignment unit, for the dynamic grid division by boundary condition of the new crack of extension, while will be described
The attribute information in new crack is assigned to net boundary, to update the net boundary；
Pressure saturation degree amending unit, for recalculating pressure field and saturation degree with the net boundary after the renewal
；
Current time walks computing unit, for correcting oil reservoir and rock mechanics parameters, carries out historical data fitting, calculates and work as
The saturation degree and pressure condition of preceding time step.
In a preferred embodiment, the oil reservoir mesh generation unit specifically for use using fracture surface as geometry about
The Delaunary subdivision methods of beam carry out mesh generation to the oil reservoir, are curved surface by crack treatment, are four sides by matrix treatments
Volume mesh.
In a preferred embodiment, the oil reservoir mesh generation unit is additionally operable to subtract matrix mesh volume adjacent
Crack grid volume half.
In a preferred embodiment, described device also includes：There is judging unit in new crack, for judging whether
Existing new crack.
It is local to introduce secondary grid according to the scheme of the embodiment of the present application on the basis of conventional discrete results of fracture simulation pair,
Complete information is described according to discrete fractures and calculates relevant parameter, so that not only accurate but also effective flow model is set up, not
The problem of discrete fractures number of grid is big to calculate slow is solved on the basis of mistake accuracy, the Numerical-Mode of fractured reservoir is improved
The efficiency of plan.In addition, the scheme provided in the embodiment of the present application is introduced, crack occurrence extension is accurate and grid property maps, real
Show the automatic division of dynamic grid, can be to extra-low permeability oil reservoirs with reference to efficient discrete technology and efficient solution technique
The simulation of overall process is carried out in the injecting process in crack occurrence and dynamic extension.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is the schematic diagram of the dynamic discrete results of fracture simulation method for the oil reservoir that the application is provided；
Fig. 2 is that a kind of flow that implements of the dynamic discrete results of fracture simulation method for oil reservoir that the embodiment of the present application is provided is shown
It is intended to；
Fig. 3 is to represent the crack grid in grid contents portion and the secondary grid of matrix in the embodiment of the present application by connection table
Connected relation schematic diagram；
The schematic diagram of saturation field when Fig. 4 is the dynamic crack of water filling generation in the embodiment of the present application；
Fig. 5 is the schematic diagram of the dynamic discrete board slit simulation apparatus for the oil reservoir that the application is provided.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
Fig. 1 shows a kind of dynamic discrete results of fracture simulation method for oil reservoir that the embodiment of the present application is provided, as shown in figure 1,
This method comprises the following steps：
Step S101：Mesh generation is carried out to oil reservoir, the grid of oil reservoir is obtained；
Step S102：Multiple secondary grid is further divided into the grid of oil reservoir；
Step S103：Grid and multiple secondary grid based on oil reservoir, calculate matrix and matrix, matrix and crack, crack
Coupling relation between crack, judges the fracture extension situation of different time points；
Step S104：Boundary condition is in new crack using extension, dynamic grid division, while by the attribute information in new crack
Net boundary is assigned to, to update net boundary.
Step S105：Oil reservoir and rock mechanics parameters are corrected, historical data fitting is carried out, the saturation of current time step is calculated
Degree and pressure condition.
The method provided in the embodiment of the present application is local to introduce secondary net on the basis of conventional discrete results of fracture simulation pair
Lattice, describe complete information according to discrete fractures and calculate relevant parameter, so that not only accurate but also effective flow model is set up,
The problem of discrete fractures number of grid is big to calculate slow is solved on the basis of true property without misalignment, the numerical value of fractured reservoir is improved
The efficiency of simulation.In addition, the method provided in the embodiment of the present application is introduced, crack occurrence extension is accurate and grid property maps,
The automatic division of dynamic grid is realized, can be to extra-low permeability oil with reference to efficient discrete technology and efficient solution technique
Hide the simulation in crack occurrence and dynamic extension progress overall process in the injecting process.
A kind of flow that implements of the dynamic discrete results of fracture simulation method for oil reservoir that Fig. 2 provides for the embodiment of the present application is shown
It is intended to, as shown in Fig. 2 this method comprises the following steps：
Step S201：Initialization model, carries out Tetrahedron subdivision to oil reservoir using Delaunary subdivision methods, sets up three
Tie up oil reservoir grid.
Wherein, Delaunary subdivisions process is using fracture surface as geometry constraint conditions.Tetrahedron is carried out by fracture to cut open
Point, crack by it is discrete be several facets, matrix by it is discrete be tetrahedral grid.
When calculating matrix mesh volume, the one of the volume of adjacent crack grid can be subtracted from matrix mesh volume
Half, once realize and volume amendment is done by a pair matrix grid adjacent with crack.
Step S202：It is according to the characteristic dimension of matrix and discrete fractures, three-dimensional oil reservoir grid in step S201 is further
It is divided into multiple secondary grid.
The explicit contact of crack discrete model and continuous model Municipal Quality Supervision Bureau can be set up by multiple secondary grid.This
Sample, secondary grid is being introduced by local, can be calculated according to the complete information of the description of discrete fractures and be met that diplopore is double to be oozed
The flow parameter of continuous model, is then input in numerical simulator according to the form through connection table and is calculated.
When calculating flow parameter, a closed boundary condition can be defined first, while injecting single-phase flow in a fluid
Body, solves the pressure distribution intended under stable state and makes corresponding pressure choropleth division.
Above-mentioned connection table is used for representing the connected relation of the crack grid in grid contents portion and the secondary grid of matrix, referring to
Schematic diagram shown in Fig. 3.
Step S203：Calculated according to Karimi-Fard methods between matrix and matrix, between matrix and crack, crack with
The fluid communication coefficient of conductivity and reservoir and pit shaft between crack, is solved by the numerical simulation based on connection table,
Obtain the pressure and intensity value under each time step.
Step 204：Set up cracking initiation and the criterion of dynamic extension during long-term injecting water, and intrinsic fracture,
Coupling relational expression between the crack (hereinafter referred to as pressure-break) that the dynamic crack of water filling generation, pressure are produced, judges not
With the fracture extension situation at time point.
Equation below can be used：
Wherein：
Saturation field figure during the dynamic crack that water filling as shown in Figure 4 is produced, which show crack shape in the injecting process
Into result.
When running into intrinsic fracture and pressure-break for dynamic crack in the injecting process, according to the pressure of crack tip and naturally
The Fluid pressure of crack and pressure-break is compared judgement, if tip pressure is more than Fluid pressure, continues along water filling dynamic
Fractuer direction extends；If tip pressure is less than Fluid pressure, along intrinsic fracture or pressure-break Directional Extension.
Step S205：The judgment criterion extended according to cracking initiation in step S204, is judged, shape in different time points
Into dynamic unstrctured grid.
In the process, if rock rupture, tessellated mesh subdivision is continued by border of new crack, while by new crack
Attribute information be assigned to net boundary, simulated as net boundary, and complete attribute mapping.
Specifically, attribute mapping is divided into three kinds of situations：
First, the crack grid element center weight of the crack grid element center of new system (system after i.e. new crack generation) and original system
Close or on the crack grid of original system, then by crack grid pressure and saturation degree indirect assignment the splitting to new system of original system
Stitch on grid；
2nd, the crack grid element center of new system is in the matrix grid of original system, by the pressure of the matrix grid of original system
On the crack grid of saturation degree indirect assignment to new system；
3rd, the matrix grid element center of new system is overlapped with the matrix grid element center of original system or in the matrix grid of original system
In, by the pressure and saturation degree indirect assignment of the matrix grid of original system to the matrix grid of new system.
Step S206：New computing device of restarting is built according to the grid system of new system and original system, pressure field is calculated
And saturation field, and stored.
Step S207：Oil reservoir and rock mechanics parameters are corrected, historical data is fitted, are calculated according to the result of fitting
The saturation degree of current time step and the situation of pressure field.
It should be noted that the process for new crack occur can also be comprised determining whether in such scheme, and based on new production
Recalculate flow of fluid in raw crack.
The starting point coordinate in wherein new crack is as the sophisticated coordinate of pre-existing fracture, and each time point step needs judgement is with newly splitting
Whether the pressure of the sophisticated grid of the connected pre-existing fracture of seam is more than critical pressure, if greater than critical pressure, then judges new crack
Produce.
It can be determined in the following way for the bearing of trend in new crack：Extend perpendicularly to minimum principal stress direction
And along biggest principal stress direction.If survey region only measures the stress of one place, the bearing of trend in all new cracks is all
Parallel to biggest principal stress direction.If it is known that the stress field in whole survey region, then by measuring the several of designated area
The size and Orientation of point stress, using these place stress as boundary condition, whole region can be simulated using FInite Element
Stress field size and Orientation.
It can be determined in the following way for the development length in new crack：The length of new fracture extension is to rely on crack
The difference of fluid pressure and minimum principal stress.Generation per a bit of new crack needs to meet two conditions, and a condition is
Whether the Fluid pressure that crack is played is more than its critical pressure, and another condition is this small new crack closer to pre-existing fracture
Whether the crack of breakpoint has all produced.
The embodiment of the present application accordingly also provides a kind of dynamic discrete board slit simulation apparatus of oil reservoir, as shown in figure 5, the device
Including：Oil reservoir mesh generation unit 501, secondary mesh generation unit 502, fluid exchange systems computing unit 503, fracture extension
Situation judging unit 304, attribute information assignment unit 505, pressure saturation degree amending unit 506 and current time step computing unit
507。
Wherein, oil reservoir mesh generation unit 501 is used to carry out mesh generation to oil reservoir, obtains the grid of oil reservoir.Secondary net
Lattice division unit 502 is used to be further divided into multiple secondary grid to the grid of the oil reservoir.Fluid exchange systems calculate single
Member 503 is used for grid and multiple secondary grid based on the oil reservoir, calculates matrix with matrix, matrix and crack, crack with splitting
Conductivity between seam, and the fluid communication coefficient of reservoir and pit shaft is calculated, and further obtain under each time step
Pressure and saturation degree.Fracture extension situation judging unit 504 is used to set up cracking initiation and sentencing that dynamic extends in the injecting process
Other criterion, and the coupling relation set up between intrinsic fracture and water filling dynamic crack, pressure-break, judge different time points
Fracture extension situation.Attribute information assignment unit 505 is used for the dynamic grid division by boundary condition of the new crack of extension, simultaneously
The attribute information in the new crack is assigned to net boundary, to update the net boundary.Pressure saturation degree amending unit 506 is used
Net boundary after with the renewal recalculates pressure field and saturation field.Current time step computing unit 507 is used to repair
Positive oil reservoir and rock mechanics parameters, carry out historical data fitting, calculate the saturation degree and pressure condition of current time step.
Specifically, above-mentioned oil reservoir mesh generation unit 501 can be specifically used using fracture surface as geometrical constraint
Delaunary subdivision methods carry out mesh generation to the oil reservoir, are curved surface by crack treatment, are tetrahedron net by matrix treatments
Lattice, and matrix mesh volume can also be subtracted to the half of the volume of adjacent crack grid.
In addition, said apparatus can also include being used to judge whether that judging unit occurs in the new crack for new crack occur.
In the 1990s, for a technology improvement can clearly distinguish be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (for the improvement of method flow).So
And, with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow is programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make
Special IC chip 2.Moreover, nowadays, substitution manually makes IC chip, and this programming is also used instead mostly
" logic compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development
Seemingly, and the source code before compiling also handy specific programming language is write, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only needing to method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
Controller can be implemented in any suitable manner, for example, controller can take such as microprocessor or processing
Device and storage can by the computer of the computer readable program code (such as software or firmware) of (micro-) computing device
Read medium, gate, switch, application specific integrated circuit (Application Specific Integrated Circuit,
ASIC), the form of programmable logic controller (PLC) and embedded microcontroller, the example of controller includes but is not limited to following microcontroller
Device：ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320, are deposited
Memory controller is also implemented as a part for the control logic of memory.
It is also known in the art that in addition to realizing controller in pure computer readable program code mode, it is complete
Controller can be caused with gate, switch, application specific integrated circuit, programmable by the way that method and step is carried out into programming in logic entirely
Logic controller realizes identical function with the form of embedded microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure in hardware component can also be considered as to the device for realizing various functions included in it.Or
Even, it not only can be able to will be the software module of implementation method for realizing that the device of various functions is considered as but also can be Hardware Subdivision
Structure in part.
System, device, module or unit that above-described embodiment is illustrated, can specifically be realized by computer chip or entity,
Or realized by the product with certain function.
For convenience of description, it is divided into various units during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Understood based on such, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft
Part product can include some instructions to cause a computer equipment (can be personal computer, server, or network
Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can
To be stored in internal memory, internal memory potentially includes the volatile memory in computer-readable medium, random access memory
And/or the form, such as read-only storage (ROM) or flash memory (flash RAM) such as Nonvolatile memory (RAM).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer-readable instruction, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only storage
(ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, the information that can be accessed by a computing device available for storage.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.It is real especially for system
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.For example：Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
Positioned at including in the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (9)
1. a kind of dynamic discrete results of fracture simulation method of oil reservoir, it is characterised in that including：
Mesh generation is carried out to oil reservoir, the grid of oil reservoir is obtained；
Multiple secondary grid is further divided into the grid of the oil reservoir；
Grid and multiple secondary grid based on the oil reservoir, calculate matrix and matrix, matrix and crack and crack and crack it
Between conductivity, and calculate the fluid communication coefficient of reservoir and pit shaft, and further obtain the pressure under each time step
And saturation degree；
Set up the criterion of cracking initiation in the injecting process and dynamic extension, and set up intrinsic fracture and water filling dynamic crack it
Between coupling relation, and the coupling relation between intrinsic fracture and pressure-break judges the crack of different time points
Extension situation；
The dynamic grid division by boundary condition of the new crack of extension, while the attribute information in the new crack is assigned into Grid Edge
Boundary, to update the net boundary；
Pressure field and saturation field are recalculated with the net boundary after the renewal；
Oil reservoir and rock mechanics parameters are corrected, historical data fitting is carried out, the saturation degree and pressure condition of current time step is calculated,
Wherein, the attribute information in the new crack, which is assigned to net boundary, includes：
If new crack produce after the crack grid element center of system and the crack grid center superposition of original system or in original system
On the crack grid of system, then the crack of the system after producing the grid pressure and saturation degree indirect assignment of original system to new crack
On grid；
If the crack grid element center of the system after new crack generation is in the matrix grid of original system, by the crack of original system
On the crack grid of system after grid pressure and saturation degree assignment to new crack generation；
If the matrix grid element center of system after new crack is produced is overlapped or in original system with the matrix grid element center of original system
Matrix grid in, then the system after producing the pressure and saturation degree indirect assignment of the matrix grid of original system to new crack
On matrix grid.
2. according to the method described in claim 1, it is characterised in that described that oil reservoir progress mesh generation is included：
The Delaunary subdivision methods using fracture surface as geometrical constraint are used to carry out mesh generation to the oil reservoir, at crack
Manage as curved surface, be tetrahedral grid by matrix treatments.
3. method according to claim 2, it is characterised in that described also to include to oil reservoir progress mesh generation：
Matrix mesh volume is subtracted to the half of the volume of adjacent crack grid.
4. according to the method described in claim 1, it is characterised in that methods described also includes：
Judge whether new crack occur.
5. method according to claim 4, it is characterised in that described to judge whether that new crack occur includes：
Using the sophisticated coordinate of pre-existing fracture as the starting point coordinate in new crack, the original for judging to be connected with new crack in each time step
Whether the pressure of raw crack tip grid is more than its critical pressure, if more than its critical pressure, being judged as new crack occur.
6. a kind of dynamic discrete board slit simulation apparatus of oil reservoir, it is characterised in that including：
Oil reservoir mesh generation unit, for carrying out mesh generation to oil reservoir, obtains the grid of oil reservoir；
Secondary mesh generation unit, multiple secondary grid is further divided into for the grid to the oil reservoir；
Fluid exchange systems computing unit, for the grid based on the oil reservoir and multiple secondary grid, calculate matrix and matrix,
Conductivity between matrix and crack and crack and crack, and the fluid communication coefficient of reservoir and pit shaft is calculated, and further
Obtain the pressure and saturation degree under each time step；
Fracture extension situation judging unit, for setting up cracking initiation and the criterion of dynamic extension in the injecting process, and builds
Vertical coupling relation between intrinsic fracture and water filling dynamic crack, and the coupling between intrinsic fracture and pressure-break
Relation, judges the fracture extension situation of different time points；
Attribute information assignment unit, for the dynamic grid division by boundary condition of the new crack of extension, while newly being split described
The attribute information of seam is assigned to net boundary, to update the net boundary；
Pressure saturation degree amending unit, for recalculating pressure field and saturation field with the net boundary after the renewal；
Current time walks computing unit, for correcting oil reservoir and rock mechanics parameters, historical data fitting is carried out, when calculating current
The saturation degree and pressure condition of spacer step,
Wherein, the attribute information in the new crack, which is assigned to net boundary, includes：
If new crack produce after the crack grid element center of system and the crack grid center superposition of original system or in original system
On the crack grid of system, then the crack of the system after producing the grid pressure and saturation degree indirect assignment of original system to new crack
On grid；
If the crack grid element center of the system after new crack generation is in the matrix grid of original system, by the crack of original system
On the crack grid of system after grid pressure and saturation degree assignment to new crack generation；
If the matrix grid element center of system after new crack is produced is overlapped or in original system with the matrix grid element center of original system
Matrix grid in, then the system after producing the pressure and saturation degree indirect assignment of the matrix grid of original system to new crack
On matrix grid.
7. device according to claim 6, it is characterised in that the oil reservoir mesh generation unit is specifically for using to split
Seam face carries out mesh generation for the Delaunary subdivision methods of geometrical constraint to the oil reservoir, is curved surface by crack treatment, by base
Matter is processed as tetrahedral grid.
8. device according to claim 7, it is characterised in that the oil reservoir mesh generation unit is additionally operable to matrix grid
Volume subtracts the half of the volume of adjacent crack grid.
9. device according to claim 6, it is characterised in that described device also includes：There is judging unit in new crack, uses
In judging whether new crack occur.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410612927.0A CN104392109B (en) | 2014-11-04 | 2014-11-04 | A kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410612927.0A CN104392109B (en) | 2014-11-04 | 2014-11-04 | A kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104392109A CN104392109A (en) | 2015-03-04 |
| CN104392109B true CN104392109B (en) | 2017-10-17 |
Family
ID=52610011
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410612927.0A Active CN104392109B (en) | 2014-11-04 | 2014-11-04 | A kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104392109B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108171420A (en) * | 2017-12-28 | 2018-06-15 | 美国德州模拟技术公司 | Non-intrusion type simulates the EDFM method and devices of complex fracture |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106326517B (en) * | 2015-07-03 | 2019-07-30 | 中国石油化工股份有限公司 | Stratiform crack-matrix hybrid grid modeling method and device |
| CN104989360A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Characterization method for dynamic fractures of low-permeability oil reservoir |
| CN106407503B (en) * | 2016-08-22 | 2019-10-29 | 南京特雷西能源科技有限公司 | Forecast Means of Reservoir Fractures and device |
| CN106326591B (en) * | 2016-08-31 | 2019-09-06 | 西南石油大学 | In hydraulic fracturing process in crack fracturing fluid pressure field acquisition methods and device |
| CN107193042B (en) * | 2017-05-11 | 2019-07-26 | 中国石油大学(北京) | Man-made fracture modeling method and device, computer storage medium based on pressure break phase |
| CN107169227B (en) * | 2017-06-09 | 2019-04-05 | 中国石油大学(华东) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well |
| CN107704685B (en) * | 2017-10-10 | 2020-04-17 | 中国石油大学(北京) | Mesh division method and device |
| CN107844637B (en) * | 2017-10-18 | 2021-03-30 | 中国石油天然气股份有限公司 | Oil deposit numerical value determination method and device |
| CN108710734B (en) * | 2018-05-04 | 2022-12-20 | 特雷西能源科技(杭州)有限公司 | Numerical simulation method and device based on grid adaptive encryption and coarsening technology |
| CN108999608A (en) * | 2018-06-14 | 2018-12-14 | 中国石油天然气股份有限公司 | A kind of low permeable glutenite rock reservoir predominant pathway recognition methods and system |
| CN109670220B (en) * | 2018-12-05 | 2019-08-27 | 西南石油大学 | A kind of horizontal well air water two-phase method for numerical simulation based on unstrctured grid |
| CN111079335B (en) * | 2019-12-20 | 2020-10-30 | 西南石油大学 | Mathematical derivation and numerical calculation method of embedded discrete fracture model |
| CN112507551B (en) * | 2020-12-09 | 2022-11-04 | 中国石油天然气股份有限公司 | Unstructured dynamic mesh generation method and device |
| CN112949119B (en) * | 2021-02-03 | 2022-11-04 | 中国石油天然气股份有限公司 | Crack form simulation method, device, equipment and medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102243680A (en) * | 2011-07-21 | 2011-11-16 | 中国科学技术大学 | Grid partitioning method and system |
| CN102507412A (en) * | 2011-11-07 | 2012-06-20 | 中国石油集团川庆钻探工程有限公司 | Fracture-matrix permeability differential judgment method for carbonate reservoir equivalent model |
| CN103279991A (en) * | 2013-04-16 | 2013-09-04 | 西南石油大学 | Method for improving oil reservoir exploitation effect by utilizing discrete fracture endpoint deformation numerical simulation |
| CN103400020A (en) * | 2013-04-16 | 2013-11-20 | 西南石油大学 | Numerical reservoir simulation method for measuring flowing situation of a plurality of intersected discrete cracks |
-
2014
- 2014-11-04 CN CN201410612927.0A patent/CN104392109B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102243680A (en) * | 2011-07-21 | 2011-11-16 | 中国科学技术大学 | Grid partitioning method and system |
| CN102507412A (en) * | 2011-11-07 | 2012-06-20 | 中国石油集团川庆钻探工程有限公司 | Fracture-matrix permeability differential judgment method for carbonate reservoir equivalent model |
| CN103279991A (en) * | 2013-04-16 | 2013-09-04 | 西南石油大学 | Method for improving oil reservoir exploitation effect by utilizing discrete fracture endpoint deformation numerical simulation |
| CN103400020A (en) * | 2013-04-16 | 2013-11-20 | 西南石油大学 | Numerical reservoir simulation method for measuring flowing situation of a plurality of intersected discrete cracks |
Non-Patent Citations (4)
| Title |
|---|
| an efficient discrete-fracture model applicable for general-purpose reservoir simulators;M.Karimi-Fard等;《society of petroleum engineers》;20040630;第9卷(第2期);摘要部分、几何离散部分、传导能力评估部分、多重交叉裂缝部分、流动模拟部分、毛细管压力效应部分，图1、图9-11 * |
| 一种矿井岩体动态裂缝快速识别算法;肖玲玲等;《计算机技术应用》;20120406;第38卷(第4期);第106-109页 * |
| 基于岩体动态裂缝识别的矿井预警系统;魏浩;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20140415;B021-134 * |
| 裂缝识别的动态模型与应用;程时清等;《石油钻采工艺》;20070415;第29卷(第2期);摘要 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108171420A (en) * | 2017-12-28 | 2018-06-15 | 美国德州模拟技术公司 | Non-intrusion type simulates the EDFM method and devices of complex fracture |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104392109A (en) | 2015-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104392109B (en) | A kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir | |
| Marco et al. | Exact 3D boundary representation in finite element analysis based on Cartesian grids independent of the geometry | |
| Kesserwani et al. | Well-balanced RKDG2 solutions to the shallow water equations over irregular domains with wetting and drying | |
| US9361414B2 (en) | Hybrid local nonmatching method for multiphase flow simulations in heterogeneous fractured media | |
| Bazyar et al. | Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method | |
| Murillo et al. | A conservative 2D model of inundation flow with solute transport over dry bed | |
| CN104346753A (en) | Cutting optimization processing system and cutting optimization processing method | |
| CN108153933A (en) | A kind of compact reservoir production forecast method and device | |
| WO2016080983A1 (en) | Discretization scheme for simulating proppant transport in dynamic fracture networks | |
| Liang | A simplified adaptive Cartesian grid system for solving the 2D shallow water equations | |
| CA2964250A1 (en) | Junction models for simulating proppant transport in dynamic fracture networks | |
| CA2961562A1 (en) | Formation fracture flow monitoring | |
| JP7435951B2 (en) | Floating point number generation method, apparatus, electronic device, storage medium and computer program for integrated circuit chip verification | |
| CN107657075A (en) | The Region Decomposition FInite Element of Darcy velocity at simulated groundwater dielectric interface | |
| CN102646167B (en) | The method for computing work amount of a kind of engineering flowing water section and device | |
| CN105513006B (en) | A kind of TrueType font profile thickness method of adjustment and device | |
| RU2017140845A (en) | SEGMENTATION OF ANATOMIC STRUCTURE BASED ON THE MODEL | |
| CA2921492C (en) | Relative permeability inversion from historical production data using viscosity ratio invariant step-function relative permeability approximations | |
| CN108595782B (en) | Calculation method for mass transfer between matrix and cracks in discrete cracks | |
| CN108090656A (en) | A kind of method and device of definite sand body connectedness | |
| CN106558346A (en) | The method and device that a kind of solid-state disk service life based on RAIM frameworks is calculated | |
| CN107939371B (en) | A kind of method and device of determining well pattern thickening feasibility | |
| CN103150445B (en) | The method for decomposing of a kind of MATLAB model variable and device | |
| CN109558614A (en) | The analogy method and system that gas flows in shale gas reservoir multi-scale facture | |
| CN107704685A (en) | A kind of Meshing Method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |