CN104392109A - Method and device for simulating dynamic discrete cracks of oil deposit - Google Patents
Method and device for simulating dynamic discrete cracks of oil deposit Download PDFInfo
- Publication number
- CN104392109A CN104392109A CN201410612927.0A CN201410612927A CN104392109A CN 104392109 A CN104392109 A CN 104392109A CN 201410612927 A CN201410612927 A CN 201410612927A CN 104392109 A CN104392109 A CN 104392109A
- 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.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
An embodiment of the invention discloses a method and a device for simulating dynamic discrete cracks of an oil deposit. The method comprises the following steps of performing the mesh generation on the oil deposit so as to obtain meshes of the oil deposit; further dividing the meshes of the oil deposit into multiple secondary meshes; calculating conductivities between matrices, between each matrix and cracks as well as between the cracks, calculating a fluid exchange coefficient between a storage layer and a pitshaft, and further obtaining pressure and saturability under every time step; building discriminant criteria of crack formation and dynamic extension in the water injection process, building coupling relationships between natural cracks and water injection dynamic cracks as well as between the natural cracks and pressure cracks, and judging crack extension situations at different time points; dynamically dividing the meshes by means of taking an extended new crack as a boundary condition, and meanwhile assigning the attribute information of the new crack to a mesh boundary so as to update the mesh boundary; correcting mechanical parameters of the oil deposit and rocks, performing the historical data fitting, and calculating the saturability and pressure situation of the present time step.
Description
Technical field
The application relates to Research Numerical Simulation Techique field, particularly relates to a kind of dynamic discrete results of fracture simulation method and apparatus of oil reservoir.
Background technology
In flow through oil reservoir situation calculating field, results of fracture simulation method is widely used.Results of fracture simulation method of the prior art mainly contains two kinds, and namely diplopore is two blends discrete fractures analogy method.
The two infiltration method of diplopore utilizes the continuum Model of equivalence to carry out the coupling of simulation fracture and matrix, crack is done the process of homogenization.In the model, Rock Matrix space provides the main storage space of fluid, and main fluid flowing then occurs between crack and crack.Fluid communication between crack and matrix uses " transmission equation " or " form factor " to describe.Although this model counting yield improves, model is too idealized, and it has many restrictions, can not accurately portray reservoir fractures state and flow condition, when fissured structure is complicated, then needs large-scale grid.In addition due to negligible amounts, comparatively wide across region, and very inadaptable to the determinative results of fracture simulation of flowing, cause the transmission equation in matrix and crack to be difficult to accurate Calculation.
Complete and dominant description is carried out according to its physical size and distributional pattern in all cracks by discrete fractures analogy method, and the flowing of crack and matrix is all controlled by Darcy's law.In this analogy method, Fractured Reservoir system is presented by high-resolution grid, and therefore discrete fractures analogy method has the accuracy of certain applicability and height.But the shortcoming of discrete fractures analogy method is that calculated amount is larger.
Summary of the invention
The object 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 accurately but also simulation fracture fast and effectively.
For solving the problems of the technologies described above, the application provides a kind of dynamic discrete results of fracture simulation method of oil reservoir on the one hand, and the method comprises:
Mesh generation is carried out to oil reservoir, obtains the grid of oil reservoir;
Be multiple secondary grid to the grid Further Division of described oil reservoir;
Based on grid and the multiple secondary grid of described oil reservoir, calculate matrix and matrix, matrix and crack, conductivity between crack and crack, and calculate the fluid communication coefficient of reservoir and pit shaft, and the pressure obtained further under each time step and saturation degree;
The criterion set up cracking initiation in the injecting process and dynamically extend, and set up the coupling relation between natural fracture and water filling dynamic crack, pressure-break, judge the fracture propagation situation of different time points;
With the new crack extended for the dynamic grid division of boundary condition, the attribute information in described new crack is assigned to net boundary, to upgrade described net boundary simultaneously;
Pressure field and saturation field is recalculated with the net boundary after described renewal;
Revise oil reservoir and rock mechanics parameters, carry out historical data matching, calculate saturation degree and the pressure condition of current time step.
In a preferred embodiment, describedly mesh generation carried out to oil reservoir comprise:
Employing is that the Delaunary subdivision method of geometrical constraint carries out mesh generation to described oil reservoir with fracture plane, and being curved surface by crack treatment, is tetrahedral grid by matrix treatments.
In a preferred embodiment, describedly mesh generation carried out to oil reservoir also comprise:
Matrix mesh volume is deducted the half of the volume of adjacent crack grid.
In a preferred embodiment, described method also comprises:
Judge whether to occur new crack.
In a preferred embodiment, judge whether described in occur that new crack comprises:
Using the most advanced and sophisticated coordinate of pre-existing fractures as the starting point coordinate in new crack, judge at each time step whether the pressure of the most advanced and sophisticated grid of pre-existing fractures be connected with new crack is greater than its emergent pressure, if be greater than its emergent pressure, is then judged as occurring new crack.
In a preferred embodiment, the attribute information in described new crack is assigned to net boundary and comprises:
If the crack grid element center of system after new crack produces is with the crack grid center superposition of original system or on the crack grid of original system, then the grid pressure of original system and saturation degree indirect assignment are given on the crack grid of the system after the generation of new crack;
If the crack grid element center of the system after new crack produces is in the matrix grid of original system, then on the crack grid of the system after the crack grid pressure of original system and saturation degree assignment being produced to new crack;
If the matrix grid element center of system after new crack produces overlaps with the matrix grid element center of original system or in the matrix grid of original system, then by the matrix grid of the system after the pressure of the matrix grid of original system and saturation degree indirect assignment to new crack generation.
The application also provides a kind of dynamic discrete board slit simulation apparatus of oil reservoir on the other hand, comprising:
Oil reservoir stress and strain model unit, for carrying out stress and strain model to oil reservoir, obtains the grid of oil reservoir;
Secondary stress and strain model unit, for being multiple secondary grid to the grid Further Division of described oil reservoir;
Fluid exchange systems computing unit, for based on the grid of described oil reservoir and multiple secondary grid, calculate matrix and matrix, matrix and crack, conductivity between crack and crack, and calculate the fluid communication coefficient of reservoir and pit shaft, and the pressure obtained further under each time step and saturation degree;
Fracture propagation situation judging unit, for the criterion setting up cracking initiation in the injecting process and dynamically extend, and sets up the coupling relation between natural fracture and water filling dynamic crack, pressure-break, judges the fracture propagation situation of different time points;
Attribute information assignment unit, for extend new crack for the dynamic grid division of boundary condition, the attribute information in described new crack is assigned to net boundary, to upgrade described net boundary simultaneously;
Pressure saturation degree amending unit, for recalculating pressure field and saturation field with the net boundary after described renewal;
Current time step computing unit, for revising oil reservoir and rock mechanics parameters, carries out historical data matching, calculates saturation degree and the pressure condition of current time step.
In a preferred embodiment, described oil reservoir stress and strain model unit is that the Delaunary subdivision method of geometrical constraint carries out mesh generation to described oil reservoir with fracture plane specifically for adopting, and being curved surface by crack treatment, is tetrahedral grid by matrix treatments.
In a preferred embodiment, described oil reservoir stress and strain model unit is also for deducting the half of the volume of adjacent crack grid by matrix mesh volume.
In a preferred embodiment, described device also comprises: judging unit appears in new crack, occurs new crack for judging whether.
According to the scheme of the embodiment of the present application on the basis that conventional discrete results of fracture simulation is right, secondary grid is introduced in local, describe complete information according to discrete fractures and calculate correlation parameter, thus set up not only accurate but also effective flow model, the basis of not losing accuracy solves the problem that the large calculating of discrete fractures number of grid is slow, improves the efficiency of the numerical simulation of fractured reservoir.In addition, the scheme provided in the embodiment of the present application introduces crack occurrence and expands accurately and grid property mapping, achieve the automatic division of dynamic grid, in conjunction with efficient discrete technology and efficient solution technique, can to the simulation carrying out overall process in extra-low permeability oil reservoirs the injecting process at crack occurrence and dynamic extension.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present application or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, the accompanying drawing that the following describes is only some embodiments recorded in the application, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the schematic diagram of the dynamic discrete results of fracture simulation method of the oil reservoir that the application provides;
Fig. 2 is the specific implementation schematic flow sheet of the dynamic discrete results of fracture simulation method of a kind of oil reservoir that the embodiment of the present application provides;
Fig. 3 is the schematic diagram being represented the crack grid in grid contents portion and the connected relation of the secondary grid of matrix in the embodiment of the present application by connection;
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 of the oil reservoir that the application provides.
Embodiment
Technical scheme in the application is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present application, technical scheme in the embodiment of the present application is clearly and completely described, obviously, described embodiment is only some embodiments of the present application, instead of whole embodiments.Based on the embodiment in the application, those of ordinary skill in the art are not making the every other embodiment obtained under creative work prerequisite, all should belong to the scope of the application's protection.
Fig. 1 shows the dynamic discrete results of fracture simulation method of a kind of oil reservoir that the embodiment of the present application provides, and as shown in Figure 1, the method comprises the steps:
Step S101: carry out mesh generation to oil reservoir, obtains the grid of oil reservoir;
Step S102: be multiple secondary grid to the grid Further Division of oil reservoir;
Step S103: based on grid and the multiple secondary grid of oil reservoir, calculates matrix and matrix, matrix and crack, coupling relation between crack and crack, judges the fracture propagation situation of different time points;
Step S104: with the new crack extended for boundary condition, dynamic grid division, is assigned to net boundary by the attribute information in new crack simultaneously, to upgrade net boundary.
Step S105: revise oil reservoir and rock mechanics parameters, carry out historical data matching, calculates saturation degree and the pressure condition of current time step.
The method provided in the embodiment of the present application is on the basis that conventional discrete results of fracture simulation is right, secondary grid is introduced in local, describe complete information according to discrete fractures and calculate correlation parameter, thus set up not only accurate but also effective flow model, the basis of not losing accuracy solves the problem that the large calculating of discrete fractures number of grid is slow, improves the efficiency of the numerical simulation of fractured reservoir.In addition, the method provided in the embodiment of the present application introduces crack occurrence and expands accurately and grid property mapping, achieve the automatic division of dynamic grid, in conjunction with efficient discrete technology and efficient solution technique, can to the simulation carrying out overall process in extra-low permeability oil reservoirs the injecting process at crack occurrence and dynamic extension.
The specific implementation schematic flow sheet of the dynamic discrete results of fracture simulation method of a kind of oil reservoir that Fig. 2 provides for the embodiment of the present application, as shown in Figure 2, the method comprises the steps:
Step S201: initialization model, adopts Delaunary subdivision method to carry out Tetrahedron subdivision to oil reservoir, sets up three-dimensional oil reservoir grid.
Wherein, Delaunary subdivision process take fracture plane as geometry constraint conditions.Carry out Tetrahedron subdivision by fracture, crack by discrete be several little, matrix by discrete be tetrahedral grid.
When calculating matrix mesh volume, the half of the volume of adjacent crack grid can be deducted from matrix mesh volume, once realize by doing volume correction to the matrix grid adjacent with crack.
Step S202: according to the characteristic dimension of matrix and discrete fractures, be multiple secondary grid by three-dimensional oil reservoir grid Further Division in step S201.
The explicit contact of crack discrete model and continuous model Municipal Quality Supervision Bureau can be set up by multiple secondary grid.Like this, introducing secondary grid by local, can calculate according to the complete information of the description of discrete fractures and meet the two flow parameter oozing continuous model of diplopore, then be input in numerical simulator according to the form through connection table and calculate.
When calculating flow parameter, first can define a closed boundary condition, injecting monophasic fluid in a fluid simultaneously, solve the pressure distribution under plan steady state (SS) and make the division of corresponding pressure choropleth.
Above-mentioned connection table is used for representing the crack grid in grid contents portion and the connected relation of the secondary grid of matrix, see the schematic diagram shown in Fig. 3.
Step S203: the fluid communication coefficient calculating conductivity between matrix and matrix, between matrix and crack, between crack and crack and reservoir and pit shaft according to Karimi-Fard method, by solving based on the numerical simulation of connection table, obtain the pressure under each time step and intensity value.
Step 204: the criterion setting up cracking initiation and dynamically extension in long-term injecting water process, and natural fracture, water filling produce dynamic crack, pressure produce crack (hereinafter referred to as pressure-break) between coupling relational expression, judge the fracture propagation situation of different time points.
Following formula can be adopted:
Wherein:
Saturation degree field pattern during the dynamic crack that water filling as shown in Figure 4 produces, which show the result of cracking initiation in the injecting process.
When running into natural fracture and pressure-break for dynamic crack in the injecting process, compare judgement according to the pressure of crack tip and the hydrodynamic pressure of natural fracture and pressure-break, if tip pressure is greater than hydrodynamic pressure, then continue along water filling dynamic crack Directional Extension; If tip pressure is less than hydrodynamic pressure, then along natural fracture or pressure-break Directional Extension.
Step S205: according to the judgment criterion of cracking initiation expansion in step S204, judge in different time points, formative dynamics unstrctured grid.
In the process, if rock burst, then with new crack for tessellated mesh subdivision is continued on border, simultaneously by the attribute information assignment in new crack to net boundary, simulate as net boundary, and complete best property of attribute mapping.
Particularly, best property of attribute mapping is divided into three kinds of situations:
One, the crack grid element center of new system (system namely after the generation of new crack) and the crack grid center superposition of original system or on the crack grid of original system, then by the crack grid pressure of original system and saturation degree indirect assignment on the crack grid of new system;
Two, 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 and saturation degree indirect assignment on the crack grid of new system;
Three, the new matrix grid element center of system overlaps with the matrix grid element center of original system or in the matrix grid of original system, by the pressure of the matrix grid of original system and saturation degree indirect assignment on the matrix grid of new system.
Step S206: build according to the grid system of new system and original system and new restart calculation element, calculating pressure field and saturation field, and store.
Step S207: revise oil reservoir and rock mechanics parameters, carry out matching to historical data, calculates the saturation degree of current time step and the situation of pressure field according to the result of matching.
It should be noted that, the process judging whether to occur new crack in such scheme, can also be comprised, and recalculate fluid flowing based on the new crack produced.
Wherein the starting point coordinate in new crack is as the most advanced and sophisticated coordinate of pre-existing fractures, and each time point step needs the pressure of the pre-existing fractures tip grid judging to be connected with new crack whether to be greater than emergent pressure, if be greater than emergent pressure, then judges that new crack produces.
Bearing of trend for new crack can be determined in the following way: bearing of trend perpendicular to least principal stress direction along biggest principal stress direction.If survey region only records the stress of the three unities, then the bearing of trend in all new cracks is all parallel to biggest principal stress direction.If know the stress field in whole survey region, then by measuring the size and Orientation of several places stress of appointed area, using these place stress as boundary condition, adopt finite element method can simulate the size and Orientation of the stress field in whole region.
Development length for new crack can be determined in the following way: the length of new fracture propagation depends on the difference of crack fluid pressure and least principal stress.Generation demand fulfillment two conditions in every a bit of new crack, condition is whether the hydrodynamic pressure on crack is played is greater than its emergent pressure, and whether closer to the crack of pre-existing fractures breakpoint all produce another condition if being this little new crack.
The corresponding dynamic discrete board slit simulation apparatus that a kind of oil reservoir is also provided of the embodiment of the present application, as shown in Figure 5, this device comprises: oil reservoir stress and strain model unit 501, secondary stress and strain model unit 502, fluid exchange systems computing unit 503, fracture propagation 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 stress and strain model unit 501, for carrying out stress and strain model to oil reservoir, obtains the grid of oil reservoir.Secondary stress and strain model unit 502 is for being multiple secondary grid to the grid Further Division of described oil reservoir.Fluid exchange systems computing unit 503 is for based on the grid of described oil reservoir and multiple secondary grid, calculate matrix and matrix, matrix and crack, conductivity between crack and crack, and calculate the fluid communication coefficient of reservoir and pit shaft, and the pressure obtained further under each time step and saturation degree.Fracture propagation situation judging unit 504 for the criterion setting up cracking initiation in the injecting process and dynamically extend, and sets up the coupling relation between natural fracture and water filling dynamic crack, pressure-break, judges the fracture propagation situation of different time points.Attribute information assignment unit 505 for extend new crack for the dynamic grid division of boundary condition, the attribute information in described new crack is assigned to net boundary, to upgrade described net boundary simultaneously.Pressure saturation degree amending unit 506 is for recalculating pressure field and saturation field with the net boundary after described renewal.Current time step computing unit 507, for revising oil reservoir and rock mechanics parameters, carries out historical data matching, calculates saturation degree and the pressure condition of current time step.
Particularly, it is that the Delaunary subdivision method of geometrical constraint carries out mesh generation to described oil reservoir that above-mentioned oil reservoir stress and strain model unit 501 specifically can adopt with fracture plane, be curved surface by crack treatment, be tetrahedral grid by matrix treatments, and matrix mesh volume can also be deducted the half of the volume of adjacent crack grid.
In addition, said apparatus can also comprise for judging whether that judging unit appears in the new crack occurring new crack.
In the nineties in 20th century, it is improvement (such as, to the improvement of the circuit structures such as diode, transistor, switch) on hardware or the improvement (improvement for method flow) on software that improvement for a technology can clearly be distinguished.But along with the development of technology, the improvement of current a lot of method flows can be considered as the direct improvement of hardware circuit.Designer nearly all obtains corresponding hardware circuit by being programmed in hardware circuit by the method flow of improvement.Therefore, the improvement of not talkative method flow just can not realize by hardware entities module.Such as, programmable logic device (PLD) (Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable GateArray, FPGA) be exactly) so a kind of integrated circuit, its logic function is determined device programming by user.Programmed voluntarily a digital display circuit " integrated " on a slice PLD by designer, and do not need chip maker to carry out the special integrated circuit (IC) chip of designing and making 2.And, nowadays, replace and manually make integrated circuit (IC) chip, this programming is also used instead mostly " logic compiler (logic compiler) ", and software realizes, it and program development is write time software compiler used similar, and source code before will compiling also handy specific programming language write, this is referred to as hardware description language (Hardware Description Language, HDL), and HDL also not only has one, but have many kinds, as ABEL (Advanced Boolean Expression Language), AHDL (Altera Hardware Description Language), Confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java HardwareDescription Language), Lava, Lola, MyHDL, PALASM, RHDL (Ruby Hardware DescriptionLanguage) etc., the most generally use VHDL (Very-High-Speed Integrated Circuit HardwareDescription Language) and Verilog2 at present.Those skilled in the art also should be clear, only needs method flow is slightly made programming in logic with above-mentioned several hardware description language and is programmed in integrated circuit, just can be easy to the hardware circuit of this logical method flow process accomplished.
Controller can realize by any suitable mode, such as, controller can be taked such as microprocessor or processor and store the computer-readable medium of the computer readable program code (such as software or firmware) that can be performed by this (micro-) processor, logic gate, switch, special IC (Application Specific Integrated Circuit, ASIC), the form of programmable logic controller (PLC) and embedding microcontroller, the example of controller includes but not limited to following microcontroller: ARC 625D, AtmelAT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320, Memory Controller can also be implemented as a part for the steering logic of storer.
Those skilled in the art also know, except realizing except controller in pure computer readable program code mode, controller can be made to realize identical function with the form of logic gate, switch, special IC, programmable logic controller (PLC) and embedding microcontroller etc. by method step being carried out programming in logic completely.Therefore this controller can be considered to a kind of hardware component, and to the structure that also can be considered as the device realizing various function in hardware component comprised in it.Or even, the device being used for realizing various function can be considered as not only can be implementation method software module but also can be structure in hardware component.
System, device, module or unit that above-described embodiment is illustrated, specifically can be realized by computer chip or entity, or be realized by the product with certain function.
For convenience of description, various unit is divided into describe respectively with function when describing above device.Certainly, the function of each unit can be realized in same or multiple software and/or hardware when implementing the application.
As seen through the above description of the embodiments, those skilled in the art can be well understood to the mode that the application can add required general hardware platform by software and realizes.Based on such understanding, the technical scheme of the application can embody with the form of software product the part that prior art contributes in essence in other words, in one typically configuration, computing equipment comprises one or more processor (CPU), input/output interface, network interface and internal memory.This computer software product can comprise the method some part described in of some instructions in order to make a computer equipment (can be personal computer, server, or the network equipment etc.) perform each embodiment of the application or embodiment.This computer software product can be stored in internal memory, internal memory may comprise the volatile memory in computer-readable medium, the forms such as random access memory (RAM) and/or Nonvolatile memory, as ROM (read-only memory) (ROM) or flash memory (flash RAM).Internal memory is the example of computer-readable medium.Computer-readable medium comprises permanent and impermanency, removable and non-removable media can be stored to realize information by any method or technology.Information can be computer-readable instruction, data structure, the module of program or other data.The example of the storage medium of computing machine comprises, but be not limited to phase transition internal memory (PRAM), static RAM (SRAM), dynamic RAM (DRAM), the random access memory (RAM) of other types, ROM (read-only memory) (ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc ROM (read-only memory) (CD-ROM), digital versatile disc (DVD) or other optical memory, magnetic magnetic tape cassette, tape magnetic rigid disk stores or other magnetic storage apparatus or any other non-transmitting medium, can be used for storing the information can accessed by computing equipment.According to defining herein, computer-readable medium does not comprise of short duration computer readable media (transitory media), as data-signal and the carrier wave of modulation.
Each embodiment in this instructions all adopts the mode of going forward one by one to describe, between each embodiment identical similar part mutually see, what each embodiment stressed is the difference with other embodiments.Especially, for system embodiment, because it is substantially similar to embodiment of the method, so description is fairly simple, relevant part illustrates see the part of embodiment of the method.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, server computer, handheld device or portable set, laptop device, multicomputer system, system, set top box, programmable consumer-elcetronics devices, network PC, small-size computer, mainframe computer, the distributed computing environment comprising above any system or equipment etc. based on microprocessor.
The application can describe in the general context of computer executable instructions, such as program module.Usually, program module comprises the routine, program, object, assembly, data structure etc. that perform particular task or realize particular abstract data type.Also can put into practice the application in a distributed computing environment, in these distributed computing environment, be executed the task by the remote processing devices be connected by communication network.In a distributed computing environment, program module can be arranged in the local and remote computer-readable storage medium comprising memory device.
Although depict the application by embodiment, those of ordinary skill in the art know, the application has many distortion and change and do not depart from the spirit of the application, and the claim appended by wishing comprises these distortion and change and do not depart from the spirit of the application.
Claims (10)
1. a dynamic discrete results of fracture simulation method for oil reservoir, is characterized in that, comprising:
Mesh generation is carried out to oil reservoir, obtains the grid of oil reservoir;
Be multiple secondary grid to the grid Further Division of described oil reservoir;
Based on grid and the multiple secondary grid of described oil reservoir, calculate matrix and matrix, matrix and crack, conductivity between crack and crack, and calculate the fluid communication coefficient of reservoir and pit shaft, and the pressure obtained further under each time step and saturation degree;
The criterion set up cracking initiation in the injecting process and dynamically extend, and set up the coupling relation between natural fracture and water filling dynamic crack, pressure-break, judge the fracture propagation situation of different time points;
With the new crack extended for the dynamic grid division of boundary condition, the attribute information in described new crack is assigned to net boundary, to upgrade described net boundary simultaneously;
Pressure field and saturation field is recalculated with the net boundary after described renewal;
Revise oil reservoir and rock mechanics parameters, carry out historical data matching, calculate saturation degree and the pressure condition of current time step.
2. method according to claim 1, is characterized in that, describedly carries out mesh generation to oil reservoir and comprises:
Employing is that the Delaunary subdivision method of geometrical constraint carries out mesh generation to described oil reservoir with fracture plane, and being curved surface by crack treatment, is tetrahedral grid by matrix treatments.
3. method according to claim 2, is characterized in that, describedly carries out mesh generation to oil reservoir and also comprises:
Matrix mesh volume is deducted the half of the volume of adjacent crack grid.
4. method according to claim 1, is characterized in that, described method also comprises:
Judge whether to occur new crack.
5. method according to claim 4, is characterized in that, described in judge whether to occur that new crack comprises:
Using the most advanced and sophisticated coordinate of pre-existing fractures as the starting point coordinate in new crack, judge at each time step whether the pressure of the most advanced and sophisticated grid of pre-existing fractures be connected with new crack is greater than its emergent pressure, if be greater than its emergent pressure, is then judged as occurring new crack.
6. method according to claim 1, is characterized in that, the attribute information in described new crack is assigned to net boundary and comprises:
If the crack grid element center of system after new crack produces is with the crack grid center superposition of original system or on the crack grid of original system, then the grid pressure of original system and saturation degree indirect assignment are given on the crack grid of the system after the generation of new crack;
If the crack grid element center of the system after new crack produces is in the matrix grid of original system, then on the crack grid of the system after the crack grid pressure of original system and saturation degree assignment being produced to new crack;
If the matrix grid element center of system after new crack produces overlaps with the matrix grid element center of original system or in the matrix grid of original system, then by the matrix grid of the system after the pressure of the matrix grid of original system and saturation degree indirect assignment to new crack generation.
7. a dynamic discrete board slit simulation apparatus for oil reservoir, is characterized in that, comprising:
Oil reservoir stress and strain model unit, for carrying out stress and strain model to oil reservoir, obtains the grid of oil reservoir;
Secondary stress and strain model unit, for being multiple secondary grid to the grid Further Division of described oil reservoir;
Fluid exchange systems computing unit, for based on the grid of described oil reservoir and multiple secondary grid, calculate matrix and matrix, matrix and crack, conductivity between crack and crack, and calculate the fluid communication coefficient of reservoir and pit shaft, and the pressure obtained further under each time step and saturation degree;
Fracture propagation situation judging unit, for the criterion setting up cracking initiation in the injecting process and dynamically extend, and sets up the coupling relation between natural fracture and water filling dynamic crack, pressure-break, judges the fracture propagation situation of different time points;
Attribute information assignment unit, for extend new crack for the dynamic grid division of boundary condition, the attribute information in described new crack is assigned to net boundary, to upgrade described net boundary simultaneously;
Pressure saturation degree amending unit, for recalculating pressure field and saturation field with the net boundary after described renewal;
Current time step computing unit, for revising oil reservoir and rock mechanics parameters, carries out historical data matching, calculates saturation degree and the pressure condition of current time step.
8. device according to claim 7, it is characterized in that, described oil reservoir stress and strain model unit is that the Delaunary subdivision method of geometrical constraint carries out mesh generation to described oil reservoir with fracture plane specifically for adopting, and being curved surface by crack treatment, is tetrahedral grid by matrix treatments.
9. device according to claim 8, is characterized in that, described oil reservoir stress and strain model unit is also for deducting the half of the volume of adjacent crack grid by matrix mesh volume.
10. device according to claim 7, is characterized in that, described device also comprises: judging unit appears in new crack, occurs new crack for judging whether.
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 true CN104392109A (en) | 2015-03-04 |
| CN104392109B 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 (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104989360A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Characterization method for dynamic fractures of low-permeability oil reservoir |
| CN106326591A (en) * | 2016-08-31 | 2017-01-11 | 西南石油大学 | Method and device for obtaining pressure field of fracturing fluid in crack during waterpower cracking process |
| CN106326517A (en) * | 2015-07-03 | 2017-01-11 | 中国石油化工股份有限公司 | Layered fracture-matrix hybrid grid modeling method and device |
| CN106407503A (en) * | 2016-08-22 | 2017-02-15 | 南京特雷西能源科技有限公司 | Reservoir fracture prediction method and apparatus |
| CN107169227A (en) * | 2017-06-09 | 2017-09-15 | 中国石油大学(华东) | The coarse grid analogy method and system of a kind of staged fracturing horizontal well |
| CN107193042A (en) * | 2017-05-11 | 2017-09-22 | 中国石油大学(北京) | Man-made fracture modeling method and device, computer-readable storage medium based on pressure break phase |
| CN107704685A (en) * | 2017-10-10 | 2018-02-16 | 中国石油大学(北京) | A kind of Meshing Method and device |
| CN107844637A (en) * | 2017-10-18 | 2018-03-27 | 中国石油天然气股份有限公司 | A kind of numerical reservoir determines method and device |
| CN108171420A (en) * | 2017-12-28 | 2018-06-15 | 美国德州模拟技术公司 | Non-intrusion type simulates the EDFM method and devices of complex fracture |
| CN108710734A (en) * | 2018-05-04 | 2018-10-26 | 南京特雷西能源科技有限公司 | Method for numerical simulation and device based on mesh adaption encryption and coarsening technique |
| CN108999608A (en) * | 2018-06-14 | 2018-12-14 | 中国石油天然气股份有限公司 | A kind of low permeable glutenite rock reservoir predominant pathway recognition methods and system |
| CN109670220A (en) * | 2018-12-05 | 2019-04-23 | 西南石油大学 | A kind of horizontal well air water two-phase method for numerical simulation based on unstrctured grid |
| CN111079335A (en) * | 2019-12-20 | 2020-04-28 | 西南石油大学 | Mathematical derivation and numerical calculation method of embedded discrete fracture model |
| CN112507551A (en) * | 2020-12-09 | 2021-03-16 | 中国石油天然气股份有限公司 | Unstructured dynamic mesh generation method and device |
| CN112949119A (en) * | 2021-02-03 | 2021-06-11 | 中国石油天然气股份有限公司 | Crack form simulation method, device, equipment and medium |
| CN113901661A (en) * | 2021-10-11 | 2022-01-07 | 中国石油大学(华东) | Load balancing method, device, equipment and storage 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 |
|---|
| M.KARIMI-FARD等: "an efficient discrete-fracture model applicable for general-purpose reservoir simulators", 《SOCIETY OF PETROLEUM ENGINEERS》 * |
| 程时清等: "裂缝识别的动态模型与应用", 《石油钻采工艺》 * |
| 肖玲玲等: "一种矿井岩体动态裂缝快速识别算法", 《计算机技术应用》 * |
| 魏浩: "基于岩体动态裂缝识别的矿井预警系统", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
Cited By (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106326517A (en) * | 2015-07-03 | 2017-01-11 | 中国石油化工股份有限公司 | Layered fracture-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 |
| CN106407503A (en) * | 2016-08-22 | 2017-02-15 | 南京特雷西能源科技有限公司 | Reservoir fracture prediction method and apparatus |
| CN106407503B (en) * | 2016-08-22 | 2019-10-29 | 南京特雷西能源科技有限公司 | Forecast Means of Reservoir Fractures and device |
| CN106326591A (en) * | 2016-08-31 | 2017-01-11 | 西南石油大学 | Method and device for obtaining pressure field of fracturing fluid in crack during waterpower cracking process |
| CN106326591B (en) * | 2016-08-31 | 2019-09-06 | 西南石油大学 | In hydraulic fracturing process in crack fracturing fluid pressure field acquisition methods and device |
| CN107193042A (en) * | 2017-05-11 | 2017-09-22 | 中国石油大学(北京) | Man-made fracture modeling method and device, computer-readable storage medium based on pressure break phase |
| CN107193042B (en) * | 2017-05-11 | 2019-07-26 | 中国石油大学(北京) | Man-made fracture modeling method and device, computer storage medium based on pressure break phase |
| CN107169227A (en) * | 2017-06-09 | 2017-09-15 | 中国石油大学(华东) | The coarse grid analogy method and system of a kind of staged fracturing horizontal well |
| CN107169227B (en) * | 2017-06-09 | 2019-04-05 | 中国石油大学(华东) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well |
| CN107704685A (en) * | 2017-10-10 | 2018-02-16 | 中国石油大学(北京) | A kind of Meshing Method and device |
| CN107704685B (en) * | 2017-10-10 | 2020-04-17 | 中国石油大学(北京) | Mesh division method and device |
| CN107844637A (en) * | 2017-10-18 | 2018-03-27 | 中国石油天然气股份有限公司 | A kind of numerical reservoir determines method and device |
| CN108171420A (en) * | 2017-12-28 | 2018-06-15 | 美国德州模拟技术公司 | Non-intrusion type simulates the EDFM method and devices of complex fracture |
| CN108710734A (en) * | 2018-05-04 | 2018-10-26 | 南京特雷西能源科技有限公司 | Method for numerical simulation and device based on mesh adaption encryption and coarsening technique |
| 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 |
| CN109670220A (en) * | 2018-12-05 | 2019-04-23 | 西南石油大学 | A kind of horizontal well air water two-phase method for numerical simulation based on unstrctured grid |
| CN111079335A (en) * | 2019-12-20 | 2020-04-28 | 西南石油大学 | Mathematical derivation and numerical calculation method of embedded discrete fracture model |
| CN112507551A (en) * | 2020-12-09 | 2021-03-16 | 中国石油天然气股份有限公司 | Unstructured dynamic mesh generation method and device |
| CN112507551B (en) * | 2020-12-09 | 2022-11-04 | 中国石油天然气股份有限公司 | Unstructured dynamic mesh generation method and device |
| CN112949119A (en) * | 2021-02-03 | 2021-06-11 | 中国石油天然气股份有限公司 | Crack form simulation method, device, equipment and medium |
| CN112949119B (en) * | 2021-02-03 | 2022-11-04 | 中国石油天然气股份有限公司 | Crack form simulation method, device, equipment and medium |
| CN113901661A (en) * | 2021-10-11 | 2022-01-07 | 中国石油大学(华东) | Load balancing method, device, equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104392109B (en) | 2017-10-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104392109A (en) | Method and device for simulating dynamic discrete cracks of oil deposit | |
| Simons et al. | A model for overland flow and associated processes within the Hydroinformatics Modelling System | |
| CA2783977C (en) | Method and system for partitioning parallel simulation models | |
| GB2506271A (en) | Modeling a fractured medium | |
| Bazyar et al. | Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method | |
| US10445446B2 (en) | Discretization scheme for simulating proppant transport in dynamic fracture networks | |
| US10626706B2 (en) | Junction models for simulating proppant transport in dynamic fracture networks | |
| CN108153933A (en) | A kind of compact reservoir production forecast method and device | |
| CN104462755A (en) | Electronic equipment spare part configuration and calculation method based on reliability model | |
| CN106775474A (en) | A kind of Nand Flash abrasion equilibrium methods, device and memory | |
| US9129075B2 (en) | Mesh generation system | |
| Bernard et al. | High-order h-adaptive discontinuous Galerkin methods for ocean modelling | |
| CN107657075A (en) | The Region Decomposition FInite Element of Darcy velocity at simulated groundwater dielectric interface | |
| Sorek et al. | Model order reduction and control polynomial approximation for well-control production optimization | |
| CN104504472A (en) | Method and device for predicting productivity of gas storage gas water interaction area | |
| Yu et al. | Efficient finite-volume model for shallow-water flows using an implicit dual time-stepping method | |
| US20140365186A1 (en) | System and method for load balancing for parallel computations on structured multi-block meshes in cfd | |
| CN107246900A (en) | A kind of method and device for determining oil-water interface position | |
| CN108595782B (en) | Calculation method for mass transfer between matrix and cracks in discrete cracks | |
| CN107146267A (en) | A kind of method and system of drawing isoline figure | |
| US20160196369A1 (en) | Relative permeability inversion from historical production data using viscosity ratio invariant step-function relative permeability approximations | |
| CN109829078A (en) | A kind of data processing method and device of raster data | |
| CN107939371B (en) | A kind of method and device of determining well pattern thickening feasibility | |
| Arbogast et al. | Von Neumann stable, implicit, high order, finite volume WENO schemes | |
| CN104519082A (en) | Expansion method and device of cloud computation |
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 |