CN109117551A - The determination method and device of oil displacement efficiency in pore constriction network model displacement simulation - Google Patents
The determination method and device of oil displacement efficiency in pore constriction network model displacement simulation Download PDFInfo
- Publication number
- CN109117551A CN109117551A CN201810901833.3A CN201810901833A CN109117551A CN 109117551 A CN109117551 A CN 109117551A CN 201810901833 A CN201810901833 A CN 201810901833A CN 109117551 A CN109117551 A CN 109117551A
- Authority
- CN
- China
- Prior art keywords
- network model
- pore constriction
- dimentional
- pore
- displacement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The embodiment of the present application provides a kind of determination method and device of oil displacement efficiency in pore constriction network model displacement simulation, this method comprises: obtaining the characteristic parameter of pore constriction medium, and constructs two-dimentional pore constriction network model according to the characteristic parameter；The percolation law of two dimension pore constriction network model described in displacement simulation；During simulating the percolation law, when the vertex pressure for confirming venturi unit in the two-dimentional pore constriction network model reaches critical pressure value, corresponding venturi unit is made to open and form crack, to update the two-dimentional pore constriction network model；During simulating the percolation law, the oil displacement efficiency of the two-dimentional pore constriction network model of different phase is obtained in real time.The embodiment of the present application can determine oil displacement efficiency in pore constriction network model displacement simulation and the accuracy of displacement simulation can be improved.
Description
Technical field
This application involves untraditional reservoir developmental research technical fields, more particularly, to a kind of pore constriction network model
The determination method and device of oil displacement efficiency in displacement simulation.
Background technique
The developmental research of untraditional reservoir has strategic importance for China's energy, their reservoir all has dual more
The feature of hole medium, micropore distribution is opened with dynamic micropore, influence research of the extension for seepage flow is to solve these resources
The key of exploitation.For the development problem of the untraditional reservoir of these complicated pore structure characteristics, it is difficult merely from macroscopic view
Seepage theory go to solve, therefore, conduct a research from microstructure characteristic and microscopic percolation mechanism, will be one new prominent
Cut.
Currently, in the untraditional reservoir exploitation simulation based on microstructure characteristic and microscopic percolation mechanism, it is used
Pore constriction network model mainly have equivalent continuum model, discrete networks model, mixed model and percolation model and hole
Gap venturi network model.Wherein, understanding of the pore constriction network model of percolation-based theory in conventional oil gas reservoir development mecha-nism
On embodied certain advantage.However, since oil reservoir development is a dynamic process, and these pore constriction network models
There is no the process for considering dynamic venturi cracking, displacement simulation process and practical displacement process have certain error.For this reason, it may be necessary to
Dynamic pore constriction network model is established, and in the displacement simulation based on dynamic pore constriction network model, how really
Determining the oil displacement efficiency in displacement simulation is a technical problem to be solved urgently.
Summary of the invention
The embodiment of the present application is designed to provide in a kind of pore constriction network model displacement simulation oil displacement efficiency really
Determine method and device, determines oil displacement efficiency in pore constriction network model displacement simulation to realize.
In order to achieve the above objectives, on the one hand, the embodiment of the present application provides a kind of pore constriction network model displacement simulation
The determination method of middle oil displacement efficiency, comprising:
The characteristic parameter of pore constriction medium is obtained, and two-dimentional pore constriction network mould is constructed according to the characteristic parameter
Type；
The percolation law of two dimension pore constriction network model described in displacement simulation；
During simulating the percolation law, when venturi unit in the confirmation two-dimentional pore constriction network model
When vertex pressure reaches critical pressure value, corresponding venturi unit is made to open and form crack, to update the two-dimentional hole larynx
Road network model；
During simulating the percolation law, the drive of the two-dimentional pore constriction network model of different phase is obtained in real time
Oily efficiency.
The determination method of oil displacement efficiency, described to obtain in real time in the pore constriction network model displacement simulation of the embodiment of the present application
Take the oil displacement efficiency of the two-dimentional pore constriction network model of different phase, comprising:
According to formulaCalculate the oil displacement efficiency of the two-dimentional pore constriction network model of different phase；
Wherein, EdFor the oil displacement efficiency of two-dimentional pore constriction network model, SwiFor the beam of two-dimentional pore constriction network model
Tie up water saturation, SorFor the residual oil saturation of two-dimentional pore constriction network model.
The determination method of oil displacement efficiency in the pore constriction network model displacement simulation of the embodiment of the present application, further includes:
During simulating the percolation law, the injection PV of the two-dimentional pore constriction network model is obtained in real time
Number.
The determination method of oil displacement efficiency in the pore constriction network model displacement simulation of the embodiment of the present application, further includes:
Relationship fitting is carried out to the injection PV number and oil displacement efficiency of the two-dimentional pore constriction network model, obtains injection PV
Number is in the relation curve of oil displacement efficiency.
The determination method of oil displacement efficiency, described to obtain in real time in the pore constriction network model displacement simulation of the embodiment of the present application
Take the injection PV number of the two-dimentional pore constriction network model, comprising:
According to formulaObtain the two-dimentional pore constriction network model walked in the n-th displacement it is total
Injection rate；
By the total injection rate of the n-th displacement step divided by the total pore size volume of the two-dimentional pore constriction network model, institute is obtained
State the injection PV number that two-dimentional pore constriction network model is walked in the n-th displacement；
Wherein, VnFor the total injection rate that two-dimentional pore constriction network model is walked in the n-th displacement, VoIt is two-dimentional under original state
The sum of the pore cell being broken into pore constriction network model and the volume of venturi unit, Vi1For two-dimentional pore constriction network mould
The volume that venturi unit in type has been broken into after the i-th displacement step, Vi2For the hole list in two-dimentional pore constriction network model
The volume that member has been broken into after the i-th displacement step.
The determination method of oil displacement efficiency, described according to institute in the pore constriction network model displacement simulation of the embodiment of the present application
It states characteristic parameter and constructs two-dimentional pore constriction network model, comprising:
Determine the statistical distribution functions for characterizing the regularity of distribution of the characteristic parameter；
Pore constriction network model is generated according to the statistical distribution functions, and the characteristic parameter is assigned to the hole
Venturi unit and pore cell in gap venturi network model, to form two-dimentional pore constriction network model.
The determination method of oil displacement efficiency in the pore constriction network model displacement simulation of the embodiment of the present application, the determining use
In the statistical distribution functions for the regularity of distribution for characterizing the characteristic parameter, comprising:
Based on formulaCharacterize the regularity of distribution of pore radius；
Based on formulaCharacterize the regularity of distribution of throat radius；
Wherein, p (R) is the distribution function of pore radius, and f (x) is the distribution function of throat radius, and R is pore radius,
RmaxFor the maximum value of pore radius, RminFor the maximum value of pore radius, σ is the standard deviation of distribution function, and e is natural constant, μ
For the desired value of distribution function, x is throat radius.
The determination method of oil displacement efficiency, the displacement mould in the pore constriction network model displacement simulation of the embodiment of the present application
Intend the percolation law of the two-dimentional pore constriction network model, comprising:
According to formulaSimulate the venturi in the two-dimentional pore constriction network model
The percolation law of unit；
According to formulaSimulate the percolation law in the crack in the two-dimentional pore constriction network model；
Wherein, qijFor the flow in two-dimentional pore constriction network model at the i-th row jth column position, rijFor two-dimentional hole larynx
Throat radius in road network model at the i-th row jth column position,For the coefficient of dynamic viscosity of fluid, LijFor two-dimentional hole larynx
To the fracture length at the i-th row jth column position, P in road network modeljThe pressure arranged for jth in two-dimentional pore constriction network model
Power, PiFor the pressure of the i-th row in two-dimentional pore constriction network model,θ is contact angle, σwnFor wetting phase
With non-wetted alternate interfacial tension, r is capillary radius, and q is the flow in pore cell, and e is natural constant, and μ is fluid
Coefficient of dynamic viscosity, Δ P=Pj-Pi-Pc。
On the other hand, the embodiment of the present application also provides oil displacement efficiencies in a kind of pore constriction network model displacement simulation
Determining device, comprising:
Network model constructs module, for obtaining the characteristic parameter of pore constriction medium, and according to the characteristic parameter structure
Build two-dimentional pore constriction network model；
Percolation law analog module, the percolation law for two dimension pore constriction network model described in displacement simulation；
Network model update module is used for during simulating the percolation law, when the confirmation two-dimentional hole larynx
When the vertex pressure of venturi unit reaches critical pressure value in road network model, corresponding venturi unit is made to open and be formed and split
Seam, to update the two-dimentional pore constriction network model；
Oil displacement efficiency obtains module, for obtaining the two of different phase in real time during simulating the percolation law
Tie up the oil displacement efficiency of pore constriction network model.
On the other hand, the embodiment of the present application also provides oil displacement efficiencies in another pore constriction network model displacement simulation
Determining device, including memory, processor and the computer program being stored on the memory, the computer journey
Sequence executes following steps when being run by the processor:
The characteristic parameter of pore constriction medium is obtained, and two-dimentional pore constriction network mould is constructed according to the characteristic parameter
Type；
The percolation law of two dimension pore constriction network model described in displacement simulation；
During simulating the percolation law, when venturi unit in the confirmation two-dimentional pore constriction network model
When vertex pressure reaches critical pressure value, corresponding venturi unit is made to open and form crack, to update the two-dimentional hole larynx
Road network model；
During simulating the percolation law, the drive of the two-dimentional pore constriction network model of different phase is obtained in real time
Oily efficiency.
It can be seen that the embodiment of the present application obtains the characteristic parameter of pore constriction medium first, and according to characteristic parameter structure
Build two-dimentional pore constriction network model；Secondly the percolation law of displacement simulation two dimension pore constriction network model；Then it is simulating
During percolation law, when the vertex pressure for confirming venturi unit in two-dimentional pore constriction network model reaches critical pressure value
When, so that corresponding venturi unit is opened and is formed crack, to update two-dimentional pore constriction network model, then in simulation seepage flow rule
During rule, the oil displacement efficiency of the two-dimentional pore constriction network model of different phase can be obtained in real time, since the application is implemented
The pore constriction network model of example considers the hole of untraditional reservoir in displacement process Dynamic cracking, and therefore, the application is real
The displacement process simulation of example and the practical displacement process of untraditional reservoir are applied more closely, to reduce untraditional reservoir
Displacement simulation process and practical displacement process error.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings:
Fig. 1 is the stream of the determination method of oil displacement efficiency in the pore constriction network model displacement simulation of one embodiment of the application
Cheng Tu；
Fig. 2 is the pore throat radius distribution map obtained in one embodiment of the application；
Fig. 3 is that matlab fitting data schematic diagram is based in one embodiment of the application；
Fig. 4 a~Fig. 4 d is the displacement path schematic diagram in one embodiment of the application using Number of Fractures as main contral parameter；
Fig. 5 a~Fig. 5 d is the displacement path schematic diagram in one embodiment of the application using fracture length as main contral parameter；
Fig. 6 a~Fig. 6 d is the displacement path schematic diagram in one embodiment of the application using displacement pressure difference as main contral parameter；
Fig. 7 a is the influence schematic diagram in one embodiment of the application under different Number of Fractures to oil displacement efficiency；
Fig. 7 b is the influence schematic diagram in one embodiment of the application under different fracture lengths to oil displacement efficiency；
Fig. 7 c is the influence schematic diagram in one embodiment of the application under different displacement pressure differences to oil displacement efficiency；
Fig. 8 is the relation curve schematic diagram of the injection PV number and oil displacement efficiency that obtain in one embodiment of the application；
Fig. 9 is the knot of the determining device of oil displacement efficiency in the pore constriction network model displacement simulation of one embodiment of the application
Structure block diagram；
Figure 10 is the determining device of oil displacement efficiency in the pore constriction network model displacement simulation of another embodiment of the application
Structural block diagram.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.Such as in being described below, second component is formed above the first component, may include the first component and second component
It can also include the first component and second component in a manner of non-direct contact (i.e. the with the embodiment that the direct way of contact is formed
Can also include additional component between one component and second component) embodiment etc. that is formed.
Moreover, for ease of description, some embodiments of the application can be used such as " in ... top ", " ... under ",
The spatially relative terms such as " top ", " lower section ", with description such as each element shown in the drawings of embodiment or component and another
Relationship between (or other) element or component.It should be understood that space is opposite other than the orientation described in the attached drawing
Term also aims to the different direction including device in use or operation.If such as the device in attached drawing is reversed, and is described
For the element or component of " " other elements or component " below " or " under ", it will then be positioned as " " other elements or portion
Part " top " or " on ".
Refering to what is shown in Fig. 1, in the pore constriction network model displacement simulation of the embodiment of the present application oil displacement efficiency determination side
Method may comprise steps of:
S101, the characteristic parameter for obtaining pore constriction medium, and two-dimentional pore constriction net is constructed according to the characteristic parameter
Network model.
In one embodiment of the application, the characteristic parameter for obtaining pore constriction medium is to obtain untraditional reservoir
The characteristic parameter of practical rock core.Wherein the characteristic parameter for example may include distribution, ligancy and the space of hole and venturi
Correlation etc. (such as shown in Fig. 2).In one exemplary embodiment of the application, true core can be obtained with nuclear-magnetism equipment etc.
Characteristic parameter.
It is described to may include according to the two-dimentional pore constriction network model of characteristic parameter building in one embodiment of the application
Following steps:
1) statistical distribution functions for characterizing the regularity of distribution of the characteristic parameter, are determined；In general, any special
Parameter is levied according to a distribution function to match with the regularity of distribution can be selected according to its regularity of distribution, so as to use this
Distribution function characterizes the regularity of distribution of this feature parameter.It, can also be with softwares such as matlab in another embodiment of the application
(such as shown in Fig. 3) is fitted to characteristic parameter, so as to characterize the distribution rule of characteristic parameter by the curve fitted
Rule.
In one exemplary embodiment of the application, such as formula can be based on
Characterize the regularity of distribution of pore radius；In one exemplary embodiment of the application, such as formula can be based onCharacterize the regularity of distribution of throat radius.
Wherein, p (R) is the distribution function of pore radius, and f (x) is the distribution function of throat radius, and R is pore radius,
RmaxFor the maximum value of pore radius, RminFor the maximum value of pore radius, σ is the standard deviation of distribution function, and e is natural constant, μ
For the desired value of distribution function, x is throat radius.
2) pore constriction network model, is generated according to the statistical distribution functions, and the characteristic parameter is assigned to institute
Venturi unit and the pore cell in pore constriction network model are stated, to form two-dimentional pore constriction network model.In this Shen
For the curve of the regularity of distribution to characteristic feature parameter fitted, hole please equally can be generated in another embodiment
Venturi network model, and the characteristic parameter is assigned to venturi unit and pore cell in pore constriction network model, from
And form two-dimentional pore constriction network model.
The percolation law of two dimension pore constriction network model described in S102, displacement simulation.
In one embodiment of the application, percolation law is referred to as fluid neuron network mechanism, is a kind of containing dynamic crack
Seepage flow mechanism.Wherein, the dynamic crack refers to during displacement simulation, under displacement pressure effect, can constantly have
The venturi of closure is opened by pressure break, so that fluid flows in the solid crack with obvious stress response.
In one exemplary embodiment of the application, pore cell can be represented with spherula, larynx can be represented with cylindrical pipe
Road unit, to form hole-venturi basic unit.Wherein, circular section unit can be used in venturi unit, and flowing law can root
According to formulaCharacterization；The flowing law in crack can be reduced to sheet flow, and can basis
FormulaCharacterization.
Wherein, qijFor the flow in two-dimentional pore constriction network model at the i-th row jth column position, rijFor two-dimentional hole larynx
Throat radius in road network model at the i-th row jth column position,For the coefficient of dynamic viscosity of fluid, LijFor two-dimentional hole larynx
To the fracture length at the i-th row jth column position, P in road network modeljThe pressure arranged for jth in two-dimentional pore constriction network model
Power, PiFor the pressure of the i-th row in two-dimentional pore constriction network model,θ is contact angle, σwnFor wetting phase
With non-wetted alternate interfacial tension, r is capillary radius, and q is the flow in pore cell, and e is natural constant, and μ is fluid
Coefficient of dynamic viscosity, Δ P=Pj-Pi-Pc, when the non-wetted phase of wetting phase displacement, Δ P takes positive sign, when non-wetted phase displacement soaks
Xiang Shi, Δ P take negative sign.
In the application another exemplary embodiment, above-mentioned formulaFor two-phase
Multiple solutions flowing.As needed, when needing to characterize two-phase list interface motion, which be can simplify are as follows:Further, as needed, when needing to characterize single-phase flow, which can be further
Simplify are as follows:
In one embodiment of the application, displacement main contral parameter can be selected as needed in displacement simulation.Show in the application one
In example property embodiment, Fig. 4 a~Fig. 4 d shows the displacement path schematic diagram using Number of Fractures as main contral parameter.Wherein, Fig. 4 a
For fissureless displacement path schematic diagram；Fig. 4 b is the displacement path schematic diagram comprising 5 cracks；Fig. 4 c is to include 10 cracks
Displacement path schematic diagram；Fig. 4 d is the displacement path schematic diagram comprising 40 cracks.In the application another exemplary embodiment
In, Fig. 5 a~Fig. 5 d shows the displacement path schematic diagram using fracture length as main contral parameter.Wherein, Fig. 5 a is fissureless
Displacement path schematic diagram；Fig. 5 b is the displacement path schematic diagram in the crack for being 600 microns comprising 15 length；Fig. 5 c is to include 15
The displacement path schematic diagram in the crack that length is 1500 microns；Fig. 5 d is the drive in the crack for being 3000 microns comprising 15 length
For path schematic diagram.In the application another exemplary embodiment, Fig. 6 a~FIG. 6d shows that joined using displacement pressure difference as master control
Several displacement path schematic diagrames.Wherein, Fig. 6 a be displacement pressure difference be 1MPa under displacement path schematic diagram；Fig. 6 b is displacement pressure difference
For the displacement path schematic diagram under 3MPa；Fig. 6 c be displacement pressure difference be 6MPa under displacement path schematic diagram；Fig. 6 d is displacement pressure
Difference is the displacement path schematic diagram under 9MPa.
S103, during simulating the percolation law, when confirming venturi in the two-dimentional pore constriction network model
When the vertex pressure of unit reaches critical pressure value, corresponding venturi unit is set to open and form crack, to update the two dimension
Pore constriction network model.
It, can for the dynamic analog process that reservoir venturi is opened by pressure break under displacement pressure effect in one embodiment of the application
Pre-set the critical pressure value of the vertex pressure of venturi unit；For each of model venturi unit, when its vertex
When pressure reaches critical pressure value, corresponding venturi unit is set to open and form crack, with the two-dimentional pore constriction net that timely updates
Network model.The venturi unit that can be understood as cracking in crack herein.
In one embodiment of the application, before the percolation law of the two-dimentional holes gap fracture network model of analogous outcrop described in displacement simulation, also
It may include steps of:
The wetting phase for determining displacement simulation is water and non-wetted phase is oil；
It is fixed value by the up-and-down boundary state assignment of each grid in the two-dimentional holes gap fracture network model of analogous outcrop, and or so
Boundary condition is assigned a value of variable value；The fixed value indicates that the state on boundary is fixed during displacement simulation, described variable
Value indicates that the state on boundary during displacement simulation is changeable；
It is 0 value, institute by the state assignment of all pore cells of the two-dimentional holes gap fracture network model of analogous outcrop and venturi unit
Stating 0 value indicates closed state and full of water.
Correspondingly, the vertex pressure of venturi unit reaches critical in the confirmation the two-dimentional holes gap fracture network model of analogous outcrop
Before pressure value, it can also include the following steps:
When starting to simulate the percolation law, by the displacement entrance boundary in the two-dimentional holes gap fracture network model of analogous outcrop
The state assignment of pore cell and venturi unit is 1 value, and 1 value indicates opening state and full of oil；
Confirm in the venturi unit being connected with the pore cell that state is 1 value, if existence is the venturi list of 0 value
Member；
If it is present judging whether the vertex pressure of the venturi unit reaches critical pressure value.When the confirmation venturi list
When the vertex pressure of member is not up to critical pressure value, the current state of the venturi unit is kept.
Correspondingly, facing when described when the vertex pressure for confirming venturi unit in the two-dimentional holes gap fracture network model of analogous outcrop reaches
When boundary's pressure value, so that corresponding venturi unit is opened and is formed crack, specifically includes:
It is 1 by the Status Change of the venturi unit when the vertex pressure for confirming the venturi unit reaches critical pressure value
Value, and be 1 value by the Status Change for the pore cell being connected with the venturi unit.
S104, during simulating the percolation law, in real time obtain different phase two-dimentional pore constriction network mould
The oil displacement efficiency of type.
In one embodiment of the application, the displacement of reservoir oil effect of the real-time two-dimentional pore constriction network model for obtaining different phase
Rate, can be according to formulaIt is calculated.In some exemplary embodiments of the application, oil displacement efficiency is obtained
Such as it can be shown in Fig. 7 a~Fig. 7 c.
Wherein, EdFor the oil displacement efficiency of two-dimentional pore constriction network model, SwiFor the beam of two-dimentional pore constriction network model
Tie up water saturation, SorFor the residual oil saturation of two-dimentional pore constriction network model.
In one embodiment of the application, during simulating the percolation law, the two-dimentional holes can also be obtained in real time
The injection PV number of gap venturi network model.Correspondingly, can also injection PV number to the two-dimentional pore constriction network model and
Oil displacement efficiency carries out relationship fitting, obtains injection PV number in the relation curve of oil displacement efficiency, as shown in Figure 8.
In one embodiment of the application, the injection PV number for obtaining the two-dimentional pore constriction network model in real time, such as
It can be according to formulaObtain total note that the two-dimentional pore constriction network model is walked in the n-th displacement
Enter amount；Then the total injection rate of the n-th displacement step is obtained into institute divided by the total pore size volume of the two-dimentional pore constriction network model
State the injection PV number that two-dimentional pore constriction network model is walked in the n-th displacement；
Wherein, VnFor the total injection rate that two-dimentional pore constriction network model is walked in the n-th displacement, VoIt is two-dimentional under original state
The sum of the pore cell being broken into pore constriction network model and the volume of venturi unit, Vi1For two-dimentional pore constriction network mould
The volume that venturi unit in type has been broken into after the i-th displacement step, Vi2For the hole list in two-dimentional pore constriction network model
The volume that member has been broken into after the i-th displacement step.
It can be seen that the embodiment of the present application obtains the characteristic parameter of pore constriction medium first, and according to characteristic parameter structure
Build two-dimentional pore constriction network model；Secondly the percolation law of displacement simulation two dimension pore constriction network model；Then it is simulating
During percolation law, when the vertex pressure for confirming venturi unit in two-dimentional pore constriction network model reaches critical pressure value
When, so that corresponding venturi unit is opened and is formed crack, to update two-dimentional pore constriction network model, then in simulation seepage flow rule
During rule, the oil displacement efficiency of the two-dimentional pore constriction network model of different phase can be obtained in real time, since the application is implemented
The pore constriction network model of example considers the hole of untraditional reservoir in displacement process Dynamic cracking, and therefore, the application is real
The displacement process simulation of example and the practical displacement process of untraditional reservoir are applied more closely, to reduce untraditional reservoir
Displacement simulation process and practical displacement process error, this and untraditional reservoir practical displacement process is more close
Displacement process simulation, can for it is subsequent formulation untraditional reservoir displacement development plan provide more objective and accurate reference according to
According to.
Refering to what is shown in Fig. 9, oil displacement efficiency is really in a kind of pore constriction network model displacement simulation of the embodiment of the present application
Determining device may include:
Network model constructs module 91, can be used for obtaining the characteristic parameter of pore constriction medium, and according to the feature
Parameter constructs two-dimentional pore constriction network model；
Percolation law analog module 92 can be used for the seepage flow rule of two dimension pore constriction network model described in displacement simulation
Rule；
Network model update module 93 can be used for during simulating the percolation law, when the confirmation two dimension
When the vertex pressure of venturi unit reaches critical pressure value in pore constriction network model, open corresponding venturi unit and shape
At crack, to update the two-dimentional pore constriction network model；
Oil displacement efficiency obtains module 94, can be used for during simulating the percolation law, obtains not same order in real time
The oil displacement efficiency of the two-dimentional pore constriction network model of section.
Refering to what is shown in Fig. 10, oil displacement efficiency in another pore constriction network model displacement simulation of the embodiment of the present application
Determining device may include memory, processor and the computer program being stored on the memory, the computer journey
Sequence executes following steps when being run by the processor:
The characteristic parameter of pore constriction medium is obtained, and two-dimentional pore constriction network mould is constructed according to the characteristic parameter
Type；
The percolation law of two dimension pore constriction network model described in displacement simulation；
During simulating the percolation law, when venturi unit in the confirmation two-dimentional pore constriction network model
When vertex pressure reaches critical pressure value, corresponding venturi unit is made to open and form crack, to update the two-dimentional hole larynx
Road network model；
During simulating the percolation law, the drive of the two-dimentional pore constriction network model of different phase is obtained in real time
Oily efficiency.
Although procedures described above process includes the multiple operations occurred with particular order, it should however be appreciated that understand,
These processes may include more or fewer operations, these operations can be executed sequentially or be executed parallel (such as using parallel
Processor or multi-thread environment).
For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each unit can be realized in the same or multiple software and or hardware when application.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/or
The forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable medium
Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves
State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable
Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM),
Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices
Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates
Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
Include, so that process, method or equipment including a series of elements not only include those elements, but also including not having
There is the other element being expressly recited, or further includes for this process, method or the intrinsic element of equipment.Do not having more
In the case where more limitations, the element that is limited by sentence "including a ...", it is not excluded that including process, the side of the element
There is also other identical elements in method or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product.
Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application
Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code
The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art
For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal
Replacement, improvement etc., should be included within the scope of the claims of this application.
Claims (10)
1. a kind of determination method of oil displacement efficiency in pore constriction network model displacement simulation characterized by comprising
The characteristic parameter of pore constriction medium is obtained, and two-dimentional pore constriction network model is constructed according to the characteristic parameter；
The percolation law of two dimension pore constriction network model described in displacement simulation；
During simulating the percolation law, when the vertex for confirming venturi unit in the two-dimentional pore constriction network model
When pressure reaches critical pressure value, corresponding venturi unit is made to open and form crack, to update the two-dimentional pore constriction net
Network model；
During simulating the percolation law, the displacement of reservoir oil effect of the two-dimentional pore constriction network model of different phase is obtained in real time
Rate.
2. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as described in claim 1
In the oil displacement efficiency of the real-time two-dimentional pore constriction network model for obtaining different phase, comprising:
According to formulaCalculate the oil displacement efficiency of the two-dimentional pore constriction network model of different phase；
Wherein, EdFor the oil displacement efficiency of two-dimentional pore constriction network model, SwiFor the irreducible water of two-dimentional pore constriction network model
Saturation degree, SorFor the residual oil saturation of two-dimentional pore constriction network model.
3. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as described in claim 1
In, further includes:
During simulating the percolation law, the injection PV number of the two-dimentional pore constriction network model is obtained in real time.
4. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as claimed in claim 3
In, further includes:
Relationship fitting is carried out to the injection PV number and oil displacement efficiency of the two-dimentional pore constriction network model, obtain injection PV number in
The relation curve of oil displacement efficiency.
5. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as claimed in claim 3
In the injection PV number for obtaining the two-dimentional pore constriction network model in real time, comprising:
According to formulaObtain total injection that the two-dimentional pore constriction network model is walked in the n-th displacement
Amount；
By the total injection rate of the n-th displacement step divided by the total pore size volume of the two-dimentional pore constriction network model, described two are obtained
The injection PV number that dimension pore constriction network model is walked in the n-th displacement；
Wherein, VnFor the total injection rate that two-dimentional pore constriction network model is walked in the n-th displacement, VoFor hole two-dimentional under original state
The sum of the pore cell being broken into venturi network model and the volume of venturi unit, Vi1For in two-dimentional pore constriction network model
Venturi unit after the i-th displacement step the volume that has been broken into, Vi2Exist for the pore cell in two-dimentional pore constriction network model
The volume being broken into after i-th displacement step.
6. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as described in claim 1
In described to construct two-dimentional pore constriction network model according to the characteristic parameter, comprising:
Determine the statistical distribution functions for characterizing the regularity of distribution of the characteristic parameter；
Pore constriction network model is generated according to the statistical distribution functions, and the characteristic parameter is assigned to the hole larynx
Venturi unit and pore cell in road network model, to form two-dimentional pore constriction network model.
7. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as claimed in claim 6
In the determination is used to characterize the statistical distribution functions of the regularity of distribution of the characteristic parameter, comprising:
Based on formulaCharacterize the regularity of distribution of pore radius；
Based on formulaCharacterize the regularity of distribution of throat radius；
Wherein, p (R) is the distribution function of pore radius, and f (x) is the distribution function of throat radius, and R is pore radius, RmaxFor
The maximum value of pore radius, RminFor the maximum value of pore radius, σ is the standard deviation of distribution function, and e is natural constant, and μ is point
The desired value of cloth function, x are throat radius.
8. the determination method of oil displacement efficiency, feature exist in pore constriction network model displacement simulation as described in claim 1
In the percolation law of two dimension pore constriction network model described in the displacement simulation, comprising:
According to formulaSimulate the venturi unit in the two-dimentional pore constriction network model
Percolation law；
According to formulaSimulate the percolation law in the crack in the two-dimentional pore constriction network model；
Wherein, qijFor the flow in two-dimentional pore constriction network model at the i-th row jth column position, rijFor two-dimentional pore constriction net
Throat radius in network model at the i-th row jth column position,For the coefficient of dynamic viscosity of fluid, LijFor two-dimentional pore constriction net
To the fracture length at the i-th row jth column position, P in network modeljFor the pressure that jth in two-dimentional pore constriction network model arranges, Pi
For the pressure of the i-th row in two-dimentional pore constriction network model,θ is contact angle, σwnFor wetting phase and non-profit
Interfacial tension between wetting phase, r are capillary radius, and q is the flow in pore cell, and e is natural constant, and μ is that the power of fluid is viscous
Stagnant coefficient, Δ P=Pj-Pi-Pc。
9. the determining device of oil displacement efficiency in a kind of pore constriction network model displacement simulation characterized by comprising
Network model constructs module, for obtaining the characteristic parameter of pore constriction medium, and according to characteristic parameter building two
Tie up pore constriction network model；
Percolation law analog module, the percolation law for two dimension pore constriction network model described in displacement simulation；
Network model update module is used for during simulating the percolation law, when the confirmation two-dimentional pore constriction net
When the vertex pressure of venturi unit reaches critical pressure value in network model, corresponding venturi unit is made to open and form crack, with
Update the two-dimentional pore constriction network model；
Oil displacement efficiency obtains module, for obtaining the two-dimentional holes of different phase in real time during simulating the percolation law
The oil displacement efficiency of gap venturi network model.
10. the determining device of oil displacement efficiency in a kind of pore constriction network model displacement simulation, including memory, processor, with
And it is stored in the computer program on the memory, which is characterized in that when the computer program is run by the processor
Execute following steps:
The characteristic parameter of pore constriction medium is obtained, and two-dimentional pore constriction network model is constructed according to the characteristic parameter；
The percolation law of two dimension pore constriction network model described in displacement simulation；
During simulating the percolation law, when the vertex for confirming venturi unit in the two-dimentional pore constriction network model
When pressure reaches critical pressure value, corresponding venturi unit is made to open and form crack, to update the two-dimentional pore constriction net
Network model；
During simulating the percolation law, the displacement of reservoir oil effect of the two-dimentional pore constriction network model of different phase is obtained in real time
Rate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810901833.3A CN109117551B (en) | 2018-08-09 | 2018-08-09 | Method and device for determining oil displacement efficiency in pore throat network model displacement simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810901833.3A CN109117551B (en) | 2018-08-09 | 2018-08-09 | Method and device for determining oil displacement efficiency in pore throat network model displacement simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109117551A true CN109117551A (en) | 2019-01-01 |
| CN109117551B CN109117551B (en) | 2022-01-04 |
Family
ID=64851645
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810901833.3A Active CN109117551B (en) | 2018-08-09 | 2018-08-09 | Method and device for determining oil displacement efficiency in pore throat network model displacement simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109117551B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
| CN102540265A (en) * | 2011-12-21 | 2012-07-04 | 西南石油大学 | Porous medium water saturation calculation method based on network simulation |
| CN103278436A (en) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure |
| CN104265255A (en) * | 2014-09-26 | 2015-01-07 | 中国石油天然气股份有限公司 | Two-dimensional microcosmic visual thickened oil replacement simulation experiment system and use method thereof |
-
2018
- 2018-08-09 CN CN201810901833.3A patent/CN109117551B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
| CN102540265A (en) * | 2011-12-21 | 2012-07-04 | 西南石油大学 | Porous medium water saturation calculation method based on network simulation |
| CN103278436A (en) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure |
| CN104265255A (en) * | 2014-09-26 | 2015-01-07 | 中国石油天然气股份有限公司 | Two-dimensional microcosmic visual thickened oil replacement simulation experiment system and use method thereof |
Non-Patent Citations (2)
| Title |
|---|
| 李胜彪: "王集复杂断块油藏开发地质综合评价及调整治理研究", 《中国博士学位论文全文数据库》 * |
| 鲁国明: "基于孔隙网络模型的砂砾岩微观水驱油驱替特征", 《特种油气藏》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109117551B (en) | 2022-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104657510B (en) | A kind of method and device that steel truss arched bridge rapid modeling is realized using BIM | |
| CN105239976A (en) | Numerical simulation method taking dynamic changes of seepage parameter of water-drive reservoir into consideration | |
| CN105696986A (en) | Novel combination flooding oil flooding experiment/test simulating method | |
| CN109711106A (en) | A kind of collection transmission pipe network optimization method and device | |
| CN108876019A (en) | A kind of electro-load forecast method and system based on big data | |
| CN109829217A (en) | Pressure break Fractured Reservoir productivity simulation method and device | |
| CN109190183A (en) | The determination method and device of macroparameter in pore constriction network model displacement simulation | |
| San Cristóbal | Time, cost, and quality in a road building project | |
| CN107939366A (en) | A kind of volume fracturing Dynamic Analysis of Horizontal Well Pressure method and device | |
| CN109063346A (en) | Consider the displacement simulation method and device of the pore constriction network model of Dynamic cracking | |
| CN105528648A (en) | Dynamic production prediction method and device of fracture-cavity unit | |
| CN108363114B (en) | Fine and close oil dessert area's evaluation method and device | |
| CN104809275A (en) | Carbonatite microflow simulating method and device | |
| CN109063348A (en) | Displacement simulation method and device based on pore constriction network model | |
| CN109033672A (en) | The dynamic crack determination method and device of pore constriction network model in displacement simulation | |
| CN109063347A (en) | Displacement simulation method and device based on pore constriction network model | |
| CN108952676A (en) | A kind of shale gas reservoir heterogeneity evaluation method and its device | |
| CN109117551A (en) | The determination method and device of oil displacement efficiency in pore constriction network model displacement simulation | |
| CN109635405A (en) | A kind of Multiphase Flow transient state calculation method and system based on space conservation | |
| Belsnes et al. | Unit commitment in hydro power operation scheduling | |
| CN108005644A (en) | A kind of inclined seam Dynamic Permeability Forecasting Methodology and device | |
| CN111677486A (en) | Compact oil carbon dioxide huff and puff simulation method and device and storage medium | |
| CN109946997A (en) | A kind of emulation dam regulation method, apparatus and system | |
| CN110969249A (en) | Production well yield prediction model establishing method, production well yield prediction method and related device | |
| CN102708226A (en) | Water conservancy and hydropower engineering simulation modeling method based on process flow |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |