CN109376481A - Calculation method, device and the computer equipment of digital cores phase percolation curve based on more GPU - Google Patents
Calculation method, device and the computer equipment of digital cores phase percolation curve based on more GPU Download PDFInfo
- Publication number
- CN109376481A CN109376481A CN201811455946.1A CN201811455946A CN109376481A CN 109376481 A CN109376481 A CN 109376481A CN 201811455946 A CN201811455946 A CN 201811455946A CN 109376481 A CN109376481 A CN 109376481A
- Authority
- CN
- China
- Prior art keywords
- oily
- water
- phase
- hole
- parameter
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The disclosure discloses the calculation method and device of a kind of digital cores phase percolation curve based on more GPU, comprising: to the digital cores according to constructed by micro- rock core in mesh space, the grid data of porosity points where obtaining and storing each hole；Oily displacement is carried out to digital cores and is saturated Fluid Dynamics, numerical value calculating is carried out according to grid data in oily displacement saturation Fluid Dynamics, obtains the oily phase parameter and water phase parameter at each moment, until oily displacement saturation Fluid Dynamics reach stable state；Oily phase parameter and water phase parameter when using stable state are as initial parameter, water drive is carried out for oil simulation to digital cores, in simulations, the phase percolation curve predicted by micro- rock core is calculated according to initial parameter and grid data, the phase percolation curve is to replace the water saturation at each moment, oleic permeability and water phase permeability in oily simulation process obtained according in the water drive.So as to be quickly obtained the phase percolation curve for supporting development decision-making.
Description
It was submitted this application claims on August 16th, 2018, entitled " a kind of digital cores based on more GPU mutually seep song
The priority of the Chinese patent application CN2018109334379 of line computation method " is herein merged entire contents by quoting
In this.
Technical field
This disclosure relates to computer application technology, in particular to a kind of digital cores phase percolation curve based on more GPU
Calculation method, device and computer equipment.
Background technique
Phase percolation curve includes oil relative permeability curve and water phase permeability saturation curve, wherein oil mutually opposite infiltration
Rate curve refers to that the relation curve between oleic permeability and water saturation, water phase permeability saturation curve refer to that water phase is permeated
Relation curve between rate and water saturation.Phase percolation curve is able to reflect oily phase, water phase fluid in porous media (such as micro- rock
The heart) in percolation law, be the significant data in oil reservoir development, it is most important for development decision-making.
In the prior art, it mainly tests, i.e., is obtained by experiment to from micro- rock core of crude oil reservoir collection in worksite
Water saturation, water phase permeability in micro- rock core, oleic permeability obtain phase percolation curve to draw.On the one hand, experiment
Time is long, supports the time of the phase percolation curve of development decision-making long so obtaining；On the other hand, experiment is carried out dependent on testing
Pressure environment, temperature environment are constructed in room.
It follows that there are also to be solved for the problem of how being quickly obtained the phase percolation curve for supporting oil extraction decision.
Summary of the invention
In order to solve the problems, such as present in the relevant technologies, present disclose provides a kind of digital cores based on more GPU mutually to seep
Calculation method, device and the computer equipment of curve.
In a first aspect, a kind of calculation method of the digital cores phase percolation curve based on more GPU, comprising:
To the digital cores according to constructed by micro- rock core in mesh space, porosity points where obtaining and storing each hole
Grid data；
Oily displacement is carried out to the digital cores and is saturated Fluid Dynamics, according to the net in the oily displacement saturation Fluid Dynamics
Lattice data carry out numerical value calculating, obtain the oily phase parameter and water phase parameter at each moment, until the oily displacement is saturated Fluid Dynamics
Reach stable state；
Oily phase parameter and water phase parameter when using the stable state carry out water to the digital cores as initial parameter
The simulation of displacement oil, in the simulation, is calculated according to the initial parameter and the grid data as micro- rock core
The phase percolation curve predicted, the phase percolation curve be according to the water drive replace each moment in oily simulation process containing water saturation
Degree, oleic permeability and water phase permeability are obtained.；
Second aspect, a kind of computing device of the digital cores phase percolation curve based on more GPU, comprising:
Module is obtained, is configured as: to the digital cores of constructed micro- rock core in mesh space, obtaining and stores each
The grid data of porosity points where hole；
Oily displacement saturated water analog module, is configured as: oily displacement carried out to the digital cores and is saturated Fluid Dynamics, according to
The grid data carries out the oily phase parameter and water phase parameter that each moment is calculated in numerical value, until the oily displacement saturated water
Simulation reaches stable state；
Water drive replaces oily analog module, is configured as: the oily phase parameter and water phase parameter when using the stable state are as at the beginning of
Beginning parameter carries out water drive to the digital cores for oil simulation, is calculated according to the initial parameter and the grid data
To the phase percolation curve predicted by micro- rock core, the phase percolation curve is each in oily simulation process according to replacing in the water drive
The water saturation at moment, oleic permeability and water phase permeability are obtained.；
In one embodiment, the computing device of the digital cores phase percolation curve based on more GPU, further includes:
Core three-dimensional gray level image obtains module, is configured as: acquisition is scanned obtained rock to micro- rock core
Heart three dimensional grey scale image；
Picture recognition module is configured as: being carried out image recognition to the core three-dimensional gray level image, is obtained micro- rock
The pore framework three-dimensional data of the heart；
Three-dimensional building module, is configured as: carrying out micro- rock according to the pore framework three-dimensional data in mesh space
The three-dimensional building of the heart, obtains the digital cores.
In one embodiment, the acquisition module includes:
Model unit acquiring unit, is configured as: according to structure class corresponding to model unit each in the digital cores
Type sign value obtains the model unit where hole；
Porosity points acquiring unit, is configured as: according to the model unit where the hole, obtaining the hole where hole
Point, the porosity points are the mesh points occupied in the mesh space of the model unit where the hole；
Grid data writing unit, is configured as: grid is written in the hole data group distributed by the porosity points
Data, the grid data include the distribution function and density function of the coordinate of the porosity points, initial configuration；
Hole array storage unit, is configured as: the hole array is stored in the corresponding GPU distributed, with
The oil displacement saturation Fluid Dynamics and the water drive are replaced in oily simulation process and are carried out pair by the hole array in the calling GPU
Answer the calculating of parameter.
In one embodiment, the computing device of the digital cores phase percolation curve based on more GPU, further includes:
Index construct module, is configured as: according to direct with the porosity points in the model unit where the hole
Adjacent several immediate neighbor porosity points, construct hole array corresponding to the porosity points and the immediate neighbor porosity points institute is right
Answer the index of hole data；
Immediate neighbor hole index information generation module, is configured as: generating immediate neighbor hole according to constructed index
Gap index information, and immediate neighbor hole index information generated is written in hole array corresponding to the hole.
In one embodiment, the computing device of the digital cores phase percolation curve based on more GPU, further includes:
Mercury analog module is pressed, is configured as: pressure mercury simulation being carried out to the digital cores, is obtained occupied by pore throat radius
Number of grid；
Judgment module is configured as: judging whether the occupied number of grid is lower than the threshold value of setting；
Execution module is configured as: if it is determined that module judges threshold value of the occupied number of grid lower than setting,
Then the grid in the mesh space is encrypted；It is set if it is determined that module judges that the occupied number of grid is not less than
Fixed threshold value then goes to oily displacement saturated water analog module.
In one embodiment, the oily displacement saturated water analog module includes:
First configuration unit, is configured as: for the digital cores configure the entrance oily phase pressure P1 and it is described go out
Mouthful water phase pressure P2, wherein P1 > P2；
Oily displacement saturated water analogue unit, is configured as: in oily phase pressure P1 and the exit of the inlet
Under the driving of the formed pressure difference of water phase low pressure P2, carries out the oily displacement and be saturated Fluid Dynamics；
First numerical calculation unit, is configured as: in the oily displacement saturated water simulation process, according to the grid number
According to numerical value calculating is carried out, the oily phase parameter and water phase parameter of each porosity points of each moment are obtained；
Stable state determination unit, is configured as: if calculating obtained oily phase parameter and water phase parameter in setting time
It inside remains unchanged, it is determined that the oil displacement saturation Fluid Dynamics reach stable state.
In one embodiment, the water drive includes: for oily analog module
Second configuration unit, is configured as: being the number after the oily displacement saturation Fluid Dynamics reach stable state
The oily phase pressure P2 of aqueous phase flow rate V1 and outlet of the rock core with posting port；
Water drive replaces oily analogue unit, is configured as: under the driving of the aqueous phase flow rate V1 of the entrance, with the stable shape
Oily phase parameter and water phase parameter when state carry out the water drive as initial parameter for oil simulation；
Second value computing unit, is configured as: replaced in oily simulation process in the water drive, according to the initial parameter and
The grid data carries out numerical value calculating, obtains water saturation, oleic permeability and the water phase permeability at each moment；
Phase percolation curve construction unit, is configured as: according to the water saturation, oleic permeability and water phase obtained
Permeability constructs to obtain the phase percolation curve predicted by micro- rock core.
The third aspect, a kind of computer equipment, comprising:
Processor；And
Memory is stored with computer-readable instruction on the memory, and the computer-readable instruction is by the processing
Device realizes method as described above when executing.
Fourth aspect, a kind of computer readable storage medium are stored thereon with computer program, the computer program quilt
Processor realizes method as described above when executing.
By the technical solution of the disclosure, using the digital cores according to constructed by micro- rock core carry out oily displacement saturated water,
Water drive obtains for oily simulation to the highly important phase percolation curve of oil extraction decision；Without being carried out to micro- rock core
Oily displacement saturated water test, water drive obtain phase percolation curve for oil experiment, according to experimental data.Thus on the one hand can basis
Digital cores are quickly obtained phase percolation curve, on the other hand the experiment condition independent of laboratory building, immediately in exploration and development
Scene can also carry out logical operation according to the disclosed method by computer equipment, obtain phase percolation curve, meet engineering
In real time, high-speed decision needs.
It should be understood that the above general description and the following detailed description are merely exemplary, this can not be limited
It is open.
Detailed description of the invention
The drawings herein are incorporated into the specification and forms part of this specification, and shows and meets implementation of the invention
Example, and in specification together principle for explaining the present invention.
Fig. 1 is the schematic diagram of implementation environment involved in disclosure shown according to an exemplary embodiment；
Fig. 2 is a kind of block diagram of server shown according to an exemplary embodiment；
Fig. 3 is a kind of calculating side of digital cores phase percolation curve based on more GPU shown according to an exemplary embodiment
The flow chart of method；
Fig. 4 be Fig. 3 corresponding embodiment step S110 before step flow chart；
Fig. 5 is the flow chart of the step S110 of Fig. 3 corresponding embodiment；
Fig. 6 be Fig. 5 corresponding embodiment step S113 after step flow chart；
Fig. 7 be Fig. 3 corresponding embodiment step S130 before step flow chart；
Fig. 8 is the flow chart of the step S130 of Fig. 3 corresponding embodiment；
Fig. 9 is the flow chart of the step S150 of Fig. 3 corresponding embodiment；
Figure 10 is a kind of calculating dress of digital cores phase percolation curve based on more GPU shown according to an exemplary embodiment
The block diagram set.
Specific embodiment
Here will the description is performed on the exemplary embodiment in detail, the example is illustrated in the accompanying drawings.Following description is related to
When attached drawing, unless otherwise indicated, the same numbers in different drawings indicate the same or similar elements.Following exemplary embodiment
Described in embodiment do not represent all embodiments consistented with the present invention.On the contrary, they be only with it is such as appended
The example of device and method being described in detail in claims, some aspects of the invention are consistent.
Fig. 1 is the schematic diagram of implementation environment involved in the disclosure shown according to an exemplary embodiment.Show at one
In example property embodiment, the executing subject portion of the calculation method for the digital cores phase percolation curve based on more GPU that the present invention is realized
Administration carries out oily displacement saturation Fluid Dynamics and water drive in computer equipment, to the digital cores of micro- rock core for oil simulation, and in mould
Numerical value calculating is carried out during intending, and replaces the water saturation that each moment is obtained in oily simulation process, water phase infiltration according in water drive
Saturating rate and oleic permeability, obtain the phase percolation curve of reflection grease two-phase percolation law in micro- rock core.
As shown in Figure 1, embodiment of the disclosure environment includes projection device 100 and computer equipment 300.Projection device
Wired or wireless network connection is established between 100 and computer equipment 300.
Wherein projection device 100 is for being scanned given micro- rock core, to obtain constructing the digital rock of micro- rock core
The significant data of the heart, i.e. core three-dimensional gray level image.
Computer equipment 300 carries out three-dimensional reconstruction according to core three-dimensional gray level image and obtains the hole firmware three of micro- rock core
Dimension module, and oily displacement saturation Fluid Dynamics, water drive are carried out based on digital cores for oil simulation, exist to obtain reflection grease two-phase
The phase percolation curve of percolation law in micro- rock core.
Fig. 2 is a kind of hardware block diagram of server shown according to an exemplary embodiment.Server 200 can be made
For computer equipment shown in FIG. 1, to execute disclosed method.
It should be noted that the server, which is one, adapts to example of the invention, it must not believe that there is provided to this
Any restrictions of the use scope of invention.The server can not be construed to need to rely on or must have in Fig. 2 to show
Illustrative server 200 in one or more component.
The hardware configuration of the server can generate biggish difference due to the difference of configuration or performance, as shown in Fig. 2, clothes
Business device 200 includes: power supply 210, interface 230, at least a memory 250 and at least central processing unit (CPU, a Central
Processing Units) 270, GPU (Graphics Processing Unit, graphics processor) 280.
Wherein, power supply 210 is used to provide operating voltage for each hardware device on server 200.
Interface 230 includes an at least wired or wireless network interface 231, at least a string and translation interface 233, at least one defeated
Enter output interface 235 and at least USB interface 237 etc., is used for and external device communication.
The carrier that memory 250 is stored as resource, can be read-only memory, random access memory, disk or CD
Deng the resource stored thereon includes operating system 251, application program 253 and data 255 etc., and storage mode can be of short duration
It stores or permanently stores.Wherein, operating system 251 is for managing and each hardware device in control server 200 and answering
It can be Windows with program 253 to realize calculating and processing of the central processing unit 270 to mass data 255
ServerTM, Mac OS XTM, UnixTM, LinuxTM, FreeBSDTM etc..Application program 253 be based on operating system 251 it
The upper computer program for completing at least one particular job, may include an at least module (being not shown in Fig. 2), each module
The series of computation machine readable instruction to server 200 can be separately included.Data 255 can be stored in disk
File etc..
Central processing unit 270 may include the processor of one or more or more, and be set as through bus and memory
250 communications, for the mass data 255 in operation and processing memory 250.
The cluster that GPU280 can be made of multiple GPU, and be set as through bus and memory 250 and central processing
Device 270 communicates, for being saturated Fluid Dynamics and water drive in oily displacement for the image procossing of progress magnanimity in oily simulation process, thus
The dependence to central processing unit 270 is reduced, guarantees the calculation process for carrying out high speed.
In addition, also can equally realize the present invention by hardware circuit or hardware circuit combination software, therefore, this hair is realized
The bright combination for being not limited to any specific hardware circuit, software and the two.
Fig. 3 is a kind of calculating side of digital cores phase percolation curve based on more GPU shown according to an exemplary embodiment
The flow chart of method.This implements numerical simulation to micro- rock core and obtains phase percolation curve method can be executed by computer equipment 300, can be with
The following steps are included:
Step S110 is obtained to the digital cores according to constructed by micro- rock core in mesh space and is stored each hole institute
In the grid data of porosity points.
Micro- rock core comes from crude oil exploration and development scene, carries out the obtained phase percolation curve of numerical simulation by micro- rock core
Reflect the phase percolation curve in micro- rock core source place (i.e. source exploration activity exploitation scene), to survey according to phase percolation curve to carry out source oil
Visit the development decision-making at exploitation scene.
Micro- rock core is a kind of porous media, is made of hole and skeleton, and skeleton is the solid phase of micro- rock core, and hole is to flow
Carrier of the body in micro- rock core.At exploration scene, oil, water, air etc. are present in the hole of micro- rock core.
Digital cores are the i.e. digital cores equal proportion three-dimensional reconstruction that is micro- rock core according to constructed by micro- rock core, number
Rock core accurately expresses the surface and internal structure of micro- rock core, the hole inside shape, micro- rock core including micro- rock core
Connected relation etc. between gap and skeletal distribution, hole.
In mesh space, grid dividing is carried out respectively along X, Y and Z three-dimensional, to be formed in mesh space multiple
Mesh point, mesh space are the three-dimensional space for being divided into multiple units.And it is located at the also phase of the digital cores in mesh space
When several units in the then space occupied by digital cores are constituted, the unit in space occupied by digital cores is referred to as
For model unit.
It necessarily include multiple mesh points in space occupied by each hole for the digital cores in mesh space,
Mesh point in space occupied by hole is porosity points.
Oily displacement saturation Fluid Dynamics and water drive hereafter are simulated for oil, are all the holes with oily and water in digital cores
In flow simulating based on, to need through fluid (oil, water) in hole at each point pressure, speed to reflect fluid
Flow regime.And in carrying out simulation process, using the fluid at porosity points as object, the corresponding oily phase carried out at porosity points
Or the calculating of water phase fluid, obtain water phase parameter, oily phase parameter and water saturation, water phase permeability, oleic permeability
Deng.
The description of flowing of the fluid in hole needs to carry out the dynamics and energy in flow process by kinematic parameter
Description.In the technical solution of the disclosure, grid data is the variable configured in fenestral porosity point, such as oil is mutually distributed letter
Number, water phase distribution function, density function etc., to be replaced in oily simulation process in oily displacement saturation Fluid Dynamics and water drive, hole
Variate-value is calculated in the correspondence that carries out of flow simulating of the grid data according to fluid in hole at point.
Step S130 carries out oily displacement to the digital cores and is saturated Fluid Dynamics, in the oily displacement saturation Fluid Dynamics
Numerical value calculating is carried out according to the grid data, obtains the oily phase parameter and water phase parameter at each moment, until the oily displacement
Saturation Fluid Dynamics reach stable state.
Oily displacement saturated water with oil with certain flow, utilizes that is, under the action of driving force (such as temperature, pressure)
Osmosis replaces the water in hole, and saturated water refers to numerically configures porosity points to water phase point entirely.To what is carried out
Oily displacement saturation Fluid Dynamics are oil and flow simulating of the water in hole during digital cores carry out this, until oil and water
Flowing in hole reaches stable state.
In the technical solution of the disclosure, the simulation of the oily displacement saturated water carried out is the simulation to oil generation process, from
And when oily displacement saturation Fluid Dynamics reach stable state, oil and water zonation when stable state in digital cores reflects exploration
Oil and water zonation in the micro- rock core in scene, simulates the initial configuration of oil-bearing rock.And in oil reservoir exploitation field, and according to first
The oil-bearing rock of beginning form assesses oil reservoir.
Oily phase parameter is parameter of the oil at each porosity points, and corresponding, water phase parameter is water at each porosity points
Parameter.Wherein oily phase parameter includes oil phase velocity, oily phase pressure, oleic permeability, water saturation etc., corresponding, water phase
Parameter includes water phase speed, water phase pressure, water phase permeability, water saturation etc..
In oily displacement saturated water simulation process, the mistake for the saturated water that oily displacement goes out in the hole of digital cores is simulated
Journey, in other words, due to simulating oil and flow process of the water in hole, so the oil in hole, at each porosity points
Phase parameter, water phase parameter are in continuous renewal.With the simulation that oil and water flow in hole, the grid data of each porosity points
It is updated, such as density function, water phase distribution function, oily phase distribution function, then according to by being configured in mesh space
Water phase parameter, the oily phase parameter of each porosity points is calculated in computation model.
It in a preferred embodiment, is numerical solution Na Wei-Stokes by the computation model configured in mesh space
The D3Q19 Lattice Boltzmann Shan-Chen model of equation (Navier-Stokes equation, abbreviation N-S equation).At it
In his embodiment, D3Q15 or D3Q27 Lattice Boltzmann Shan-Chen model can also be configured, herein without specifically limiting
It is fixed.
To which the grid data at each porosity points be substituted into computation model and be calculated at each porosity points of each moment
Water phase, oily phase parameter.
Step S150, the oily phase parameter and water phase parameter when using the stable state are as initial parameter, to the number
Rock core carry out water drive for oil simulation, in the simulation, according to the initial parameter and the grid data be calculated for
The phase percolation curve that micro- rock core is predicted, the phase percolation curve are to replace each moment in oily simulation process according in the water drive
Water saturation, oleic permeability and water phase permeability it is obtained.
When oily displacement saturation Fluid Dynamics reach stable state, oil and water in micro- rock core are simulated in the hole of digital cores
The initial configuration of distribution, in the stable state, the stable pressure field of stroke in hole, i.e., in each porosity points, oil pressure,
The oil phase such as water pressure parameter, water phase parameter reach stable.
The water drive carried out for the process that oil simulation is that simulation is recovered the oil, thus need with explore oil in the micro- rock core in scene and
The initial configuration of water distribution is starting point.And the simulation that step S130 is carried out reaches this effect.
It is similar with oily displacement saturated water simulation process, it is also based on and is counted by the computation model configured in mesh space
It calculates, i.e., initial parameter and mesh parameter is substituted into computation model, carry out oily phase parameter, the water phase parameter of each porosity points
Calculating, and then water saturation, the oily phase at each moment is calculated according to the oily phase parameter, the water phase parameter that are calculated
Permeability, water phase permeability, and according to the water saturation at each moment, oleic permeability, water phase permeability construct to obtain for
The phase percolation curve that micro- rock core is predicted.
Wherein, water saturation=water phase porosity points/porosity points；
A certain moment oleic permeability k=- μ v0/dP, wherein μ is oil viscosity, v0It is oil along the average speed in displacement direction,
dPFor oily phase pressure gradient；It can similarly be calculated according to water viscosity, water along the average speed in displacement direction, water phase barometric gradient
To the water phase permeability at each moment.
Water drive replaces the computation model configured in oily simulation process and configured in oily displacement saturated water simulation process
Computation model is identical.To which in a preferred embodiment, the computation model configured is numerical solution Na Wei-Stokes side
The D3Q19 Lattice Boltzmann Shan-Chen model of journey (Navier-Stokes equation, abbreviation N-S equation).
It is noted that since phase percolation curve is the relationship of water phase permeability, oleic permeability and water saturation, institute
To be replaced in oily simulation process in water drive, need to be calculated water drive for the water saturation of oil simulation original state, as initially
Water saturation.And the water saturation at the moment be also equal to be saturated in oily displacement it is aqueous full when Fluid Dynamics reach stable state
And degree, so that initial water saturation (also known as irreducible water can also be calculated when oily displacement saturated water reaches stable state
Saturation degree), which describes the initial oil-containing situation at the scene of recovering the oil.
As described above, phase percolation curve includes oil relative permeability curve and water phase permeability saturation curve, thus,
Constructed respectively according to water saturation obtained, water phase permeability and oleic permeability obtain oil relative permeability curve and
Water phase permeability saturation curve, and then obtain the phase percolation curve predicted by micro- rock core.So as to anti-according to phase percolation curve institute
The oil reservoir mirrored carries out oil recovery development decision-making.Such as it can determine that irreducible water saturation, calculating grease are being stored up using phase percolation curve
Vertical characteristics in layer determine table, calculate oil displacement efficiency and reservoir water drive recovery ratio, judge reservoir wettability etc., and
These data are particularly significant to oil field development.
It in one embodiment, can be further bent to mutually seeping to keep obtained phase percolation curve smooth and parsingization
Line is fitted, such as on the basis of Corey equation, i.e., the mutually infiltration using Corey equation as fitting template, after being fitted
Curve.
By the technical solution of the disclosure, using the digital cores according to constructed by micro- rock core carry out oily displacement saturated water,
Water drive obtains for oily simulation to the highly important phase percolation curve of oil extraction decision；Without being carried out to micro- rock core
Oily displacement saturated water test, water drive obtain phase percolation curve for oil experiment, according to experimental data.Thus on the one hand can basis
Digital cores are quickly obtained phase percolation curve, immediately can also at exploration and development scene on the other hand independent of experimental site requirement
To carry out logical operation according to the disclosed method by computer equipment, phase percolation curve is obtained, the real-time, quick of engineering is met
Decision needs.
In one embodiment, as shown in figure 4, before step S110, further includes:
Step S010, acquisition are scanned obtained core three-dimensional gray level image to micro- rock core.
The scanning of micro- rock core is wherein carried out using high-resolution computed tomography equipment, the resolution ratio of CT scan equipment directly affects
To constructed digital cores, thus, in order to guarantee to be saturated Fluid Dynamics, water drive in oily displacement for the calculating essence in oily simulation process
Degree, is scanned micro- rock core using the CT scan equipment for meeting resolution requirement.In one embodiment, CT scan equipment
The size of sweep object is required equipped with size, i.e., the sweep object (such as micro- rock core) that size requires is exceeded for size, no
The rock core unit gray level image of micro- rock core structure can completely be reacted.For example, if the size of setting requires to be length
No more than 5 millimeters, then carrying out the length of micro- rock core of numerical simulation object no more than 5 millimeters.
It should further illustrate, used CT scan equipment is micron CT scan equipment, during the scanning process, by
The X-ray of micron CT scan equipment transmitting penetrates micro- rock core, due to micro- rock core different parts to the absorptivity of X-ray not
Together, so that correspondence is imaged, the core three-dimensional gray level image of micro- rock core is obtained.
Step S030 carries out image recognition to the core three-dimensional gray level image, obtains the pore framework of micro- rock core
Three-dimensional data.
Core three-dimensional gray level image presents the three-dimensional of hole and skeleton point in the image information and micro- rock core of micro- rock core
Cloth, so as to carry out image recognition to core three-dimensional gray level image, the pore framework obtained for constructing digital cores is three-dimensional
Data.
Carrying out image recognition to core three-dimensional gray level image is to carry out the segmentation and identification of hole and skeleton, and carry out in turn
The mark of hole and skeleton obtains pore framework three-dimensional data.In other words, in pore framework three-dimensional data, to hole and
Skeleton has carried out corresponding structure type mark, such as with 1 mark skeleton, with 0 mark hole etc..Wherein image recognition is used
Algorithm can be watershed algorithm and decision Tree algorithms and other thresholding algorithms, interactive machine learning algorithm etc., herein
Without specifically limiting.
Step S050 carries out the three-dimensional of micro- rock core according to the pore framework three-dimensional data in mesh space and constructs,
Obtain the digital cores.
Micro- rock core surface, the distribution of pores, skeletal distribution reflected according to pore framework three-dimensional data, in grid sky
Between the middle three-dimensional building for carrying out micro- rock core, obtain digital cores, that is, realize three-dimensional space of micro- rock core constructed by computer
In description.To replace the simulation of oily process in oily displacement saturated water and water drive, can based on constructed digital cores into
Row, without carrying out fluid in the flow simulating of micro- rock core using micro- rock core in the lab.
In one embodiment, as shown in figure 5, step S110 includes:
Step S111 obtains hole according to structure type sign value corresponding to model unit each in the digital cores
The model unit at place.
As described above, having carried out the mesh space of grid dividing, digital cores also accordingly are carried out being divided into several
Model unit, i.e. digital cores are made of several model units.
Structure type designation is carried out to hole and skeleton in pore framework three-dimensional data, thus in constructed digital rock
Building hole and skeleton are carried out according to structure type mark in the heart, so as to directly tell hole and bone from digital cores
Frame, in other words, for constructed digital cores, hole and skeleton are corresponding with structure type sign value.To
According to structure type sign value, the model unit where hole is determined.
Step S112 obtains the porosity points where hole according to the model unit where the hole, and the porosity points are
Model unit where hole mesh point occupied in the mesh space.
As before, due to the grid dividing carried out to three-dimensional space, several mesh points are formd in three-dimensional space, and hole institute
Model unit in mesh point be porosity points, corresponding, the mesh point in model unit where skeleton is skeletal point.
Grid data, the grid data is written in step S113 in the hole data group distributed by the porosity points
The distribution function and density function of coordinate, initial configuration including the porosity points.
The hole array is stored in the corresponding GPU distributed by step S114, in the oily displacement saturated water mould
The water drive is fitted to replace in oily simulation process by calling the hole array in the GPU to carry out the calculating of corresponding parameter.
In order to carry out storing grid data corresponding to each porosity points, each porosity points configuration hole array is corresponded to, is led to
Hole array corresponding to porosity points is crossed to correspond to the grid data for getting porosity points.
In calculating in step S130, S150, it is related to millions of secondary loop iterations, each circulation needs are several hundred million
Complicated floating-point operation is done on to billions of grids, and grid data is stored in the corresponding GPU distributed, thus by more
A GPU Parallel Implementation numerical procedure, the computational efficiency of equivalent high-performance CPU cluster, improves arithmetic speed.Moreover, adopting
Multiple streams are opened according to the number of GPU with the distinctive stream concurrent technique of GPU, the asynchronous communication between GPU are realized, to guarantee numerical value
The progress of calculating.
Specifically, hole array corresponding to each porosity points is segmented according to GPU uniform amount, thus by hole array pair
It should be stored in distributed GPU, also ensure the load balancing of each GPU.
Due in the simulation process of step S130 and step S150, it is only necessary to the grid data of porosity points is calculated,
To only store the grid data of porosity points, for being not involved in data corresponding to the skeletal point of calculating without storage, thus
The storage burden of GPU can be reduced.
In one embodiment, due to storing hole array by multiple GPU, in the simulation of step S130 and step S150
In the calculating of journey, the data exchange between GPU is certainly existed, thus, it is right in hole array in order to improve data exchange efficiency
The address information that different GPU need to transmit and receive data should be stored, thus after the completion of GPU operation, according to address information into
The storage of row corresponding data realizes data and is carried out in GPU without finding calculated result storage address after the completion of operation
Floating-point operation before and after high-speed transfer.
It in one embodiment, further include immediate neighbor porosity points index information in hole array, as shown in fig. 6, step
Before S114, further includes:
Step S210, according to several direct with the porosity points direct neighbor in the model unit where the hole
Neighbours' porosity points construct the rope of hole data corresponding to hole array corresponding to the porosity points and the immediate neighbor porosity points
Draw.
Step S220 generates immediate neighbor hole index information according to constructed index, and by direct neighbour generated
Hole index information is occupied to be written in hole array corresponding to the hole.
In calculating process, by the grid data of the direct porosity points of hole index information direct index porosity points, from
And data read operation is carried out, further improve the speed of data processing.
Flow simulating of the fluid in digital cores certainly exists flow direction of the fluid in hole, in other words flows
Body moves to next porosity points according to a certain direction at porosity points.In one embodiment, due to constructed Bohr
Hereby graceful grid, such as D3Q19 grid model, at a specific porosity points, such as the hole positioned at D3Q19 grid model center
Point, fluid can past 18 porosity points with the porosity points direct neighbor move, to correspond to each direct neighbor porosity points,
Corresponding that migratory direction information is written in immediate neighbor hole index information, it is straight at this which indicates fluid
The adjoining direction of motion information for occupying porosity points movement.
In one embodiment, as shown in fig. 7, before step S130, further includes:
Step S310 carries out pressure mercury simulation to digital cores, obtains number of grid occupied by pore throat radius.
Step S320, judges whether occupied number of grid is lower than the threshold value of setting.
Step S330, if it has, then being encrypted to the grid in mesh space；If it has not, then executing to digital cores
The step of carrying out oily displacement saturation Fluid Dynamics.
In micro- rock core, it is connected between hole and hole by throat, the throat of pore throat radius, that is, interconnected pore and hole
Radius, wherein pore throat radius is can be measured by the maximum radius of sphericity of venturi.
In one embodiment, mercury process is pressed by hole waveforms method numerical simulation, to obtain in digital cores
Pore throat radius distribution.
By setting the threshold value of number of grid, the precision in the numerical value calculating of step S130, S150 thereby may be ensured that.
In one embodiment, set threshold value is 3 grids, that is, requires grid occupied by pore throat radius not low
In 3 grids, if the requirements are not met, then encrypts to grid.
Mercury simulation is pressed, that is, simulates mercury from the entrance of digital cores through the process of flow of pore to outlet, thus, pass through and presses mercury
Simulation may determine that in digital cores whether each hole is connected to by throat, thus in the micro- rock core of corresponding reaction each hole whether
It is connected to by throat.And the simulation that step S130, S150 is carried out is based on each hole in hole firmware threedimensional model by shouting
Road connection, thus, corresponding to the digital cores not being connected to by throat, it is impossible to be used in the simulation in step S130, S150 is carried out,
So need replacing for throat connection micro- rock core rebuild digital cores carry out disclosed technique scheme in correspondence step.
In one embodiment, as shown in figure 8, step S130 includes:
Step S131 matches the oily phase pressure P1 of posting port and the water phase pressure P2 of outlet for digital cores, wherein P1 > P2.
Step S132, in the oily phase pressure P1 of the inlet and the formed pressure difference of water phase low pressure P2 in the exit
Driving under, carry out the oily displacement and be saturated Fluid Dynamics.
Step S133 carries out numerical value calculating according to the grid data, obtains in the oily displacement saturated water simulation process
To the oily phase parameter and water phase parameter of each porosity points of each moment.
Step S134, if the obtained oily phase parameter of calculating and water phase parameter remain unchanged within the set time, it is determined that
The oil displacement saturation Fluid Dynamics reach stable state.
It is noted that the corresponding embodiment of Fig. 8 and Fig. 9, the computation model configured is that numerical value mentioned above is asked
Solve the D3Q19 Lattice Boltzmann of Navier Stokes equation (Navier-Stokes equation, abbreviation N-S equation)
Shan-Chen model.The oily phase pressure P1 of the inlet configured in step S131 the and water phase pressure P2 in exit is basis
What the requirement of the computation model was configured.It certainly, can be by adjusting the pressure in step S131 in numerical procedure
P1, P2 re-start the simulation process of oily displacement saturated water, so that the result being calculated meets the actual conditions of oil reservoir storage.
In one embodiment, P1 and P2 are set with P1:P2=10:6.
In one embodiment, as shown in figure 9, step S150 includes:
Step S151 matches posting port after the oily displacement saturation Fluid Dynamics reach stable state for the digital cores
Aqueous phase flow rate V1 and outlet oily phase pressure P2；
Step S152, under the driving of the aqueous phase flow rate V1 of the entrance, oily phase parameter when with the stable state and
Water phase parameter carries out the water drive as initial parameter for oil simulation；
Step S153 is replaced in oily simulation process in the water drive, is carried out according to the initial parameter and the grid data
Numerical value calculates, and obtains water saturation, oleic permeability and the water phase permeability at each moment；
Step S154, according to the water saturation, oleic permeability and water phase permeability obtained construct to obtain for
The phase percolation curve that micro- rock core is predicted.
The water drive that Fig. 9 corresponding embodiment is realized is infused with constant speed from the entrance of digital cores for oil simulation i.e. simulation water
Enter, thus with water drive for oil.
It is similar with the embodiment of Fig. 8, the aqueous phase flow rate V1 and exit oil phase of the inlet configured in step S151
Pressure P2 is also to be configured according to the requirement of computation model.Similarly, simulation can also be re-started by adjusting V1 and P2, makes mould
Quasi- result more meets the actual conditions in oil reservoir.
In one embodiment, it in order to reduce the computing load of computer equipment, extracts obtained at set time intervals
Water saturation, oleic permeability and water phase permeability, according to water saturation corresponding to extracted various time points,
Oleic permeability, i.e. water phase permeability construct to obtain the phase percolation curve predicted by micro- rock core.Certainly, which is only to show
Example property citing, is not construed as the limitation to disclosure use scope.
Following is embodiment of the present disclosure, can be used for executing the above-mentioned computer equipment 300 of the disclosure executes based on
The calculation method embodiment of the digital cores phase percolation curve of more GPU.For those undisclosed details in the apparatus embodiments, it asks
Referring to the calculation method embodiment of digital cores phase percolation curve of the disclosure based on more GPU.
Figure 10 is a kind of calculating dress of digital cores phase percolation curve based on more GPU shown according to an exemplary embodiment
The block diagram set, the device can be used in the computer equipment 300 of implementation environment shown in Fig. 1, execute in above method embodiment
All or part of step shown in any.As shown in Figure 10, which includes but is not limited to: obtaining module 110, oily displacement is satisfied
Oily analog module 150 is replaced with Fluid Dynamics module 130 and water drive, in which:
Module 110 is obtained, is configured as: to the digital cores according to constructed by micro- rock core in mesh space, being obtained simultaneously
The grid data of porosity points where storing each hole.
Oily displacement saturated water analog module 130, is configured as: oily displacement carried out to the digital cores and is saturated Fluid Dynamics,
Numerical value calculating is carried out according to the grid data in the oily displacement saturation Fluid Dynamics, obtain each moment oily phase parameter and
Water phase parameter, until the oily displacement saturation Fluid Dynamics reach stable state.
Water drive replaces oily analog module 150, is configured as: oily phase parameter and water phase parameter when using the stable state as
Initial parameter carries out water drive for oil simulation, in the simulation, according to the initial parameter and described to the digital cores
The phase percolation curve predicted by micro- rock core is calculated in grid data, and the phase percolation curve is to replace oil according in the water drive
The water saturation at each moment, oleic permeability and water phase permeability are obtained in simulation process.
The function of modules and the realization process of effect are specifically detailed in the above-mentioned digital rock based on more GPU in above-mentioned apparatus
The realization process of step is corresponded in the calculation method of heart phase percolation curve, details are not described herein.
It is appreciated that these modules can by hardware, software, or a combination of both realize.When realizing in hardware
When, these modules may be embodied as one or more hardware modules, such as one or more specific integrated circuits.When with software side
When formula is realized, these modules may be embodied as the one or more computer programs executed on the one or more processors, example
The program of storage in memory 204 as performed by the processor 218 of Fig. 2.
In one embodiment, the computing device of the digital cores phase percolation curve based on more GPU, further includes:
Core three-dimensional gray level image obtains module, is configured as: acquisition is scanned obtained rock core three to micro- rock core
Tie up gray level image.
Picture recognition module is configured as: carrying out image recognition to core three-dimensional gray level image, it is three-dimensional to obtain pore framework
Data.
Three-dimensional building module, is configured as: carrying out the three-dimensional of micro- rock core according to pore framework three-dimensional data in mesh space
Building, obtains the digital cores of micro- rock core.
In one embodiment, obtaining module 110 includes:
Model unit acquiring unit, is configured as: according to structure class corresponding to model unit each in the digital cores
Type sign value obtains the model unit where hole.
Porosity points acquiring unit, is configured as: according to the model unit where the hole, obtaining the hole where hole
Point, the porosity points are the mesh points occupied in the mesh space of the model unit where the hole.
Grid data writing unit, is configured as: grid data is written in the hole data group distributed by porosity points,
Grid data includes the distribution function and density function of the coordinate of mesh point, initial configuration.
Hole array storage unit, is configured as: hole array being stored in the corresponding GPU distributed, to drive in oil
It replaces in oily simulation process for saturation Fluid Dynamics and water drive by calling the hole array in GPU to carry out the calculating of corresponding parameter.
In one embodiment, the computing device of the digital cores phase percolation curve based on more GPU, further includes:
Index construct module, is configured as: if according in the model unit where hole with porosity points direct neighbor
Dry immediate neighbor porosity points construct the rope of hole data corresponding to hole array corresponding to porosity points and immediate neighbor porosity points
Draw.
Immediate neighbor hole index information generation module, is configured as: generating immediate neighbor hole according to constructed index
Gap index information, and immediate neighbor hole index information generated is written in hole array corresponding to hole.
In one embodiment, the computing device of the digital cores phase percolation curve based on more GPU, further includes:
Mercury analog module is pressed, is configured as: pressure mercury simulation being carried out to digital cores, obtains grid occupied by pore throat radius
Quantity.
Judgment module is configured as: judging whether occupied number of grid is lower than the threshold value of setting.
Execution module is configured as: if it is determined that threshold value of the number of grid occupied by module judgement lower than setting, then right
Grid in mesh space is encrypted；If it is determined that threshold value of the number of grid occupied by module judgement not less than setting, then
Go to oily displacement saturated water analog module.
In one embodiment, oily displacement saturated water analog module 130 includes:
First configuration unit, is configured as: the water phase pressure of oily phase pressure P1 and outlet for digital cores with posting port
P2, wherein P1 > P2.
Oily displacement saturated water analogue unit, is configured as: in the oily phase pressure P1 of inlet and the water phase low pressure in exit
Under the driving of the formed pressure difference of P2, carries out oily displacement and be saturated Fluid Dynamics.
First numerical calculation unit, is configured as: in oily displacement saturated water simulation process, being counted according to grid data
Value calculates, and obtains the oily phase parameter and water phase parameter of each porosity points of each moment.
Stable state determination unit, is configured as: if calculating obtained oily phase parameter and water phase parameter in setting time
It inside remains unchanged, it is determined that oily displacement saturation Fluid Dynamics reach stable state.
In one embodiment, water drive includes: for oily analog module 170
Second configuration unit, is configured as: after oily displacement saturation Fluid Dynamics reach stable state, configuring for digital cores
The aqueous phase flow rate V1 of the entrance and oily phase pressure P2 of outlet.
Water drive replaces oily analogue unit, is configured as: the oil under the driving of the aqueous phase flow rate V1 of entrance, when with stable state
Phase parameter and water phase parameter carry out water drive as initial parameter and simulate for oil.
Second value computing unit, is configured as: replacing in oily simulation process in water drive, according to initial parameter and grid data
Numerical value calculating is carried out, water saturation, oleic permeability and the water phase permeability at each moment are obtained；
Phase percolation curve construction unit, is configured as: according to the water saturation, oleic permeability and water phase obtained
Permeability constructs to obtain the phase percolation curve predicted by micro- rock core.
Modules/unit function and the realization process of effect are specifically detailed in above-mentioned based on more in above-mentioned apparatus embodiment
The realization process of step is corresponded in the calculation method of the digital cores phase percolation curve of GPU, details are not described herein.
A kind of computer equipment, comprising:
Processor；And
Memory is stored with computer-readable instruction on memory, realization when computer-readable instruction is executed by processor
The method of any one of above method embodiment.
The processor of device in the embodiment executes the concrete mode of operation in relation to being somebody's turn to do the number based on more GPU
Detailed description is performed in the embodiment of the calculation method of rock core phase percolation curve, no detailed explanation will be given here.
Optionally, the disclosure also provides a kind of computer readable storage medium, is stored thereon with computer program, computer
The method of any one of above method embodiment is realized when program is executed by processor.
It should be understood that the present invention is not limited to the precise structure already described above and shown in the accompanying drawings, and
And various modifications and change can executed without departing from the scope.The scope of the present invention is limited only by the attached claims.
Claims (10)
1. a kind of calculation method of the digital cores phase percolation curve based on more GPU characterized by comprising
To the digital cores according to constructed by micro- rock core in mesh space, the grid of porosity points where obtaining and storing each hole
Data；
Oily displacement is carried out to the digital cores and is saturated Fluid Dynamics, according to the grid number in the oily displacement saturation Fluid Dynamics
According to numerical value calculating is carried out, the oily phase parameter and water phase parameter at each moment are obtained, until the oily displacement saturation Fluid Dynamics reach
Stable state；
Oily phase parameter and water phase parameter when using the stable state carry out water drive to the digital cores and replace as initial parameter
Oil simulation is calculated in the simulation according to the initial parameter and the grid data as micro- rock core institute in advance
The phase percolation curve of survey, the phase percolation curve be according to the water drive replace the water saturation at each moment in oily simulation process,
Oleic permeability and water phase permeability are obtained.
2. the method according to claim 1, wherein it is described in mesh space according to constructed by micro- rock core
Digital cores, before the grid data of porosity points where obtaining and storing each hole, further includes:
Acquisition is scanned obtained core three-dimensional gray level image to micro- rock core；
Image recognition is carried out to the core three-dimensional gray level image, obtains the pore framework three-dimensional data of micro- rock core；
The three-dimensional of micro- rock core is carried out according to the pore framework three-dimensional data in mesh space to construct, and obtains the digital rock
The heart.
3. the method according to claim 1, wherein it is described in mesh space according to constructed by micro- rock core
Digital cores, the grid data of porosity points where obtaining and storing each hole, comprising:
According to structure type sign value corresponding to model unit each in the digital cores, the model list where hole is obtained
Member；
According to the model unit where the hole, the porosity points where hole are obtained, the porosity points are the hole places
Model unit mesh point occupied in the mesh space；
Grid data is written in the hole data group distributed by the porosity points, the grid data includes the porosity points
Coordinate, initial configuration distribution function and density function；
The hole array is stored in the corresponding GPU distributed, in the oily displacement saturation Fluid Dynamics and the water drive
For the calculating for carrying out corresponding parameter in oily simulation process by the hole array in the calling GPU.
4. according to the method described in claim 3, it is characterized in that, further including immediate neighbor hole index in the hole data
Information, it is described the hole array is stored in the corresponding GPU distributed before, further includes:
According to several immediate neighbor porosity points in the model unit where the hole with the porosity points direct neighbor, structure
Build the index of hole data corresponding to hole array corresponding to the porosity points and the immediate neighbor porosity points；
Immediate neighbor hole index information is generated according to constructed index, and by immediate neighbor hole index information generated
It is written in hole array corresponding to the hole.
5. the method according to claim 1, wherein described carry out oily displacement saturated water mould to the digital cores
It is quasi-, the oily phase parameter that each moment is calculated in numerical value is carried out according to the grid data in the oily displacement saturation Fluid Dynamics
With water phase parameter, until before the oily displacement saturation Fluid Dynamics reach stable state, further includes:
Pressure mercury simulation is carried out to the digital cores, obtains number of grid occupied by pore throat radius；
Judge whether the occupied number of grid is lower than the threshold value of setting；
If it has, then being encrypted to the grid in the mesh space；Oil is carried out to the digital cores if it has not, then executing
Displacement is saturated the step of Fluid Dynamics.
6. described to institute the method according to claim 1, wherein the digital cores include an inlet and an outlet
It states digital cores and carries out oily displacement saturation Fluid Dynamics, counted in the oily displacement saturation Fluid Dynamics according to the grid data
The oily phase parameter and water phase parameter at each moment is calculated in value, until the oily displacement saturation Fluid Dynamics reach stable state,
Include:
The oily phase pressure P1 of the entrance and the water phase pressure P2 of the outlet are configured for the digital cores, wherein P1 > P2；
Under the driving of the formed pressure difference of water phase low pressure P2 of the oily phase pressure P1 and exit of the inlet, carry out
The oil displacement is saturated Fluid Dynamics；
In the oily displacement saturated water simulation process, numerical value calculating is carried out according to the grid data, it is each to obtain each moment
The oily phase parameter and water phase parameter of porosity points；
If calculating obtained oily phase parameter and water phase parameter remaining unchanged within the set time, it is determined that the oil displacement saturation
Fluid Dynamics reach stable state.
7. the method according to claim 1, wherein it is described with the stable state when oily phase parameter and water phase
Parameter carries out water drive for oil simulation, according to the initial parameter and the grid as initial parameter, to the digital cores
The phase percolation curve predicted by micro- rock core is calculated in data, comprising:
After the oily displacement saturation Fluid Dynamics reach stable state, aqueous phase flow rate V1 for the digital cores with posting port and
The oily phase pressure P2 of outlet；
Under the driving of the aqueous phase flow rate V1 of the entrance, the oily phase parameter and water phase parameter when using the stable state are as just
Beginning parameter carries out the water drive for oil simulation；
It is replaced in oily simulation process in the water drive, numerical value calculating is carried out according to the initial parameter and the grid data, is obtained
Water saturation, oleic permeability and the water phase permeability at each moment；
It constructs to obtain as micro- rock core institute in advance according to the water saturation, oleic permeability and water phase permeability obtained
The phase percolation curve of survey.
8. a kind of computing device of the digital cores phase percolation curve based on more GPU characterized by comprising
Module is obtained, is configured as: to the digital cores according to constructed by micro- rock core in mesh space, being obtained and store each hole
The grid data of porosity points where gap；
Oily displacement saturated water analog module, is configured as: carrying out oily displacement to the digital cores and is saturated Fluid Dynamics, in the oil
Displacement is saturated in Fluid Dynamics and carries out numerical value calculating according to the grid data, obtains the oily phase parameter and water phase ginseng at each moment
Number, until the oily displacement saturation Fluid Dynamics reach stable state；
Water drive replaces oily analog module, is configured as: oily phase parameter and water phase parameter when using the stable state are as initially joining
Number carries out water drive for oil simulation, in the simulation, according to the initial parameter and the grid number to the digital cores
According to the phase percolation curve predicted by micro- rock core is calculated, the phase percolation curve is to simulate according in the water drive for oil
The water saturation at each moment, oleic permeability and water phase permeability are obtained in journey.
9. a kind of computer equipment characterized by comprising
Processor；And
Memory is stored with computer-readable instruction on the memory, and the computer-readable instruction is held by the processor
The method as described in any one of claims 1 to 7 is realized when row.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The method as described in any one of claims 1 to 7 is realized when being executed by processor.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2018109334379 | 2018-08-16 | ||
| CN201810933437 | 2018-08-16 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109376481A true CN109376481A (en) | 2019-02-22 |
| CN109376481B CN109376481B (en) | 2022-12-09 |
Family
ID=65376314
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811455946.1A Active CN109376481B (en) | 2018-08-16 | 2018-11-30 | Multi-GPU-based digital core-facies permeability curve calculation method and device and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109376481B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109916934A (en) * | 2019-04-03 | 2019-06-21 | 青岛科技大学 | A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data |
| CN111141768A (en) * | 2020-01-10 | 2020-05-12 | 清能艾科（深圳）能源技术有限公司 | Method and device for acquiring rock wettability change in crude oil displacement |
| CN111242009A (en) * | 2020-01-10 | 2020-06-05 | 清能艾科（深圳）能源技术有限公司 | Core fingerprint correlation method and correlation device |
| CN111624147A (en) * | 2020-04-16 | 2020-09-04 | 中国石油天然气股份有限公司 | Relative permeability measuring method and device for rock core |
| CN112666059A (en) * | 2020-12-14 | 2021-04-16 | 中国石油大学(华东) | Method for determining gas-water relative permeability of porous medium in gas hydrate decomposition process |
| CN113158594A (en) * | 2021-04-19 | 2021-07-23 | 中国海洋石油集团有限公司 | Oil displacement efficiency analysis method based on cast sheet microcosmic displacement simulation |
| CN113569407A (en) * | 2021-07-26 | 2021-10-29 | 中国石油大学(北京) | Method and device for calculating capillary force and relative permeability curve |
| CN115455794A (en) * | 2022-08-10 | 2022-12-09 | 武汉理工大学 | LBM parallel optimization method and device based on connected pore division calculation area and storage medium |
| CN117269000A (en) * | 2023-11-23 | 2023-12-22 | 中国石油大学(华东) | Method for measuring two-phase relative permeability of compact rock core |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010148628A1 (en) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | Quantitative calculation method for hydrocarbon saturation of fractured reservoir |
| CN102915406A (en) * | 2012-10-26 | 2013-02-06 | 中国石油大学(华东) | Calculation method for relative permeability curve of oil and water in radial flow condition |
| CN105606787A (en) * | 2015-12-31 | 2016-05-25 | 中国石油天然气股份有限公司 | Device and method for testing capillary pressure curve of rock core |
| CA2974829A1 (en) * | 2015-02-03 | 2016-08-11 | Schlumberger Canada Limited | Modeling of fluid introduction and/or fluid extraction elements in simulation of coreflood experiment |
| CN106093083A (en) * | 2016-07-01 | 2016-11-09 | 龙威 | A kind of method using digital cores simulation to set up the relationship between lithology and logging |
-
2018
- 2018-11-30 CN CN201811455946.1A patent/CN109376481B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010148628A1 (en) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | Quantitative calculation method for hydrocarbon saturation of fractured reservoir |
| CN102915406A (en) * | 2012-10-26 | 2013-02-06 | 中国石油大学(华东) | Calculation method for relative permeability curve of oil and water in radial flow condition |
| CA2974829A1 (en) * | 2015-02-03 | 2016-08-11 | Schlumberger Canada Limited | Modeling of fluid introduction and/or fluid extraction elements in simulation of coreflood experiment |
| CN105606787A (en) * | 2015-12-31 | 2016-05-25 | 中国石油天然气股份有限公司 | Device and method for testing capillary pressure curve of rock core |
| CN106093083A (en) * | 2016-07-01 | 2016-11-09 | 龙威 | A kind of method using digital cores simulation to set up the relationship between lithology and logging |
Non-Patent Citations (2)
| Title |
|---|
| 尹洪军等: "水湿储层水驱油微观驱替特征的网络模拟研究", 《数学的实践与认识》 * |
| 祝春生等: "特低渗透砂岩油藏渗流特性研究", 《油气地质与采收率》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109916934A (en) * | 2019-04-03 | 2019-06-21 | 青岛科技大学 | A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data |
| CN111242009B (en) * | 2020-01-10 | 2023-06-13 | 清能艾科（深圳）能源技术有限公司 | Core fingerprint correlation method and equipment |
| CN111141768A (en) * | 2020-01-10 | 2020-05-12 | 清能艾科（深圳）能源技术有限公司 | Method and device for acquiring rock wettability change in crude oil displacement |
| CN111242009A (en) * | 2020-01-10 | 2020-06-05 | 清能艾科（深圳）能源技术有限公司 | Core fingerprint correlation method and correlation device |
| CN111624147A (en) * | 2020-04-16 | 2020-09-04 | 中国石油天然气股份有限公司 | Relative permeability measuring method and device for rock core |
| CN112666059A (en) * | 2020-12-14 | 2021-04-16 | 中国石油大学(华东) | Method for determining gas-water relative permeability of porous medium in gas hydrate decomposition process |
| CN113158594A (en) * | 2021-04-19 | 2021-07-23 | 中国海洋石油集团有限公司 | Oil displacement efficiency analysis method based on cast sheet microcosmic displacement simulation |
| CN113569407A (en) * | 2021-07-26 | 2021-10-29 | 中国石油大学(北京) | Method and device for calculating capillary force and relative permeability curve |
| CN113569407B (en) * | 2021-07-26 | 2022-08-26 | 中国石油大学(北京) | Method and device for calculating capillary force and relative permeability curve |
| CN115455794A (en) * | 2022-08-10 | 2022-12-09 | 武汉理工大学 | LBM parallel optimization method and device based on connected pore division calculation area and storage medium |
| CN115455794B (en) * | 2022-08-10 | 2024-03-29 | 武汉理工大学 | LBM parallel optimization method, device and storage medium based on communication pore division calculation region |
| CN117269000A (en) * | 2023-11-23 | 2023-12-22 | 中国石油大学(华东) | Method for measuring two-phase relative permeability of compact rock core |
| CN117269000B (en) * | 2023-11-23 | 2024-02-09 | 中国石油大学(华东) | Method for measuring two-phase relative permeability of compact rock core |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109376481B (en) | 2022-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109376481A (en) | Calculation method, device and the computer equipment of digital cores phase percolation curve based on more GPU | |
| Groen et al. | Analysing and modelling the performance of the HemeLB lattice-Boltzmann simulation environment | |
| Butts et al. | Geographical variability and network structure | |
| An et al. | GPU-accelerated volumetric lattice Boltzmann method for porous media flow | |
| CN105354873B (en) | For the pattern density function analogy method of porous media three-dimensionalreconstruction | |
| CN105938628A (en) | Direct computation of image-derived biomarkers | |
| Chen et al. | Local curvature entropy-based 3D terrain representation using a comprehensive Quadtree | |
| AU2013297036B2 (en) | Methods and systems related to hydrocarbon recovery strategy development | |
| Hernandez et al. | Simulating and visualizing real-time crowds on GPU clusters | |
| Bernaschi et al. | Petaflop hydrokinetic simulations of complex flows on massive GPU clusters | |
| CN109102568A (en) | Vascular flow analogy method and relevant apparatus based on Region Decomposition | |
| Raimbault et al. | A hybrid network/grid model of urban morphogenesis and optimization | |
| Ahmadi-Badejani et al. | An image-based geometric model for numerical simulation of blood perfusion within the liver lobules | |
| Pothuaud et al. | A new computational efficient approach for trabecular bone analysis using beam models generated with skeletonized graph technique | |
| Seifollahi et al. | An enhanced stochastic optimization in fracture network modelling conditional on seismic events | |
| Elshakhs et al. | A comprehensive survey on Delaunay triangulation: applications, algorithms, and implementations over CPUs, GPUs, and FPGAs | |
| Varga et al. | An extensible, generic environmental process modelling framework with an example for a watershed of a shallow lake | |
| Zhang et al. | High-fidelity tetrahedral mesh generation from medical imaging data for fluid-structure interaction analysis of cerebral aneurysms | |
| König et al. | Cognitive computing for urban planning | |
| Jung et al. | VoroCrack3d: An annotated semi-synthetic 3d image data set of cracked concrete | |
| Kemgue et al. | From spheres to ellipsoids: Speeding up considerably the morphological modeling of pore space and water retention in soils | |
| Huang et al. | A novel method for multiscale digital core reconstruction based on regional superposition algorithm | |
| Veselý | Building massing generation using GAN trained on Dutch 3D city models | |
| McCaskill et al. | Pore-to-Core Upscaling of Two-Phase Flow in Mixed-Wet Porous Media: Part I—Seamless Pore-Network Extraction | |
| Yazici | A parametric landscape urbanism method: The search for an optimal solution |
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 |