CN109063383A - Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model - Google Patents
Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model Download PDFInfo
- Publication number
- CN109063383A CN109063383A CN201811092588.2A CN201811092588A CN109063383A CN 109063383 A CN109063383 A CN 109063383A CN 201811092588 A CN201811092588 A CN 201811092588A CN 109063383 A CN109063383 A CN 109063383A
- Authority
- CN
- China
- Prior art keywords
- heat
- model
- stream
- fluid
- rock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
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
Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model that this application discloses a kind of.This method utilizes the micro- CT image of rock, blowhole mesostructure grid model based on reconstruction, Gu constructing heat-stream-THM coupling mathematical model, and the prediction of stress and temperature change to rock pore structure evolution and flow behavior is realized by numerical simulation study.Different from heat-THEORY OF FLUID-SOLID COUPLING FLOW THROUGH research of traditional macro scale, reconstruction model of this method based on pore scale, Gu heat-stream-multi- scenarios method interaction mechanism is deeply disclosed from micro-scale.Not only facilitate the prediction of stress and temperature parameter for development effect influence in oilfield exploitation procedure, simultaneously, also can promote underground nuclear waste disposal, CO2 seal up for safekeeping and the geotechnical engineerings fields such as geothermal exploitation in be related to the research of porous media multi- scenarios method problem, preferably instruct the development of engineering practice.
Description
Technical field
The present invention relates to rock physics technical fields, micro- CT image are utilized more particularly to a kind of, in conjunction with minute yardstick model
Gu rock fever-stream-multi- scenarios method analogy method of reconstruction technique.
Background technique
Studies have shown that the rock in nature is as a kind of strong heterogeneity porous media, wherein often there is gas phase,
The fluid-wall interaction system of liquid phase and solid phase interaction.Macroscopic physical mechanical characteristic and fluid of the fluid-wall interaction system to rock
Transport property play important control action, while this kind of fluid-wall interaction system also tends to be controlled by temperature in nature and becomes
The influence of change.A large amount of engineering practices show: such as natural petroleum gas field, geothermal energy utilization, CO2 geological storage and water power and tunnel
In the construction of the engineering field relevant to rock such as road, the phase there is multi- scenarios methods such as stress field, seepage field and temperature fields
Interaction, interactional dynamic process, so that the influence factor of Process of Engineering Construction more complicates.Therefore, carry out
The further investigation of rock fever-fluid-wall interaction mechanism is more conducive to the actual demand of engineering practice.
All in all, it is based on existing research achievement, rock fever-THEORY OF FLUID-SOLID COUPLING FLOW THROUGH is mainly based upon at present
The basis of Terzaghi principle of effective stress and Biot three-dimensional soil solidifying theory is made by the amendment to original basic assumption
Equation is closer to true rock environmental effect.However, at this stage often for the research of rock fever-THEORY OF FLUID-SOLID COUPLING FLOW THROUGH
Based on the continuous media thought of macroscopic view, rock matrix and pore-fluid are characterized into porous media, office as a continuous volume elements
It is limited to macro-scale, to have ignored in rock a large amount of existing complicated unordered pore structures for material properties of rock, stream
The influence of body transport property and thermodynamic property.Gu and these exactly carry out rock fever-stream-multi- scenarios method mechanism study base
Plinth.Since in the rock of true environment, solid skeletal and pore-fluid respectively occupy certain space, and skeleton is in external environment
Under the action of deform, so as to cause the continuous change of pore shape；Simultaneously in the case where heterogeneous fluid preservation, fluid
Property also constantly changes with distributional pattern.Therefore, the application discloses a kind of based on rock CT scan image, utilizes pore scale weight
Gu established model carries out rock fever-stream-multi-field coupling theory research numerical analysis method.This method establishes skeleton and hole respectively
Pore scale structured grid model, preferably can disclose heat-stream-solid phase interaction machine in porous media from microcosmic angle
System provides the technological guidance of science for engineering practice relevant to rock.
Summary of the invention
Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model that the purpose of the present invention is to provide a kind of.
In order to achieve the above object, the application adopts the following technical solutions to realize:
(1) it is basic data with the micro- CT image of rock core, pore scale is constructed by the Model Reconstruction algorithm of independent development
Structured grid model constructs the grid model of solid skeletal and pore-fluid respectively.
(2) on the basis of permeation fluid mechanics and THEORY OF FLUID-SOLID COUPLING FLOW THROUGH, by introducing parameter pair relevant to temperature change
Rock deformation is modified with fluid item, establishes heat-fluid-wall interaction mathematical model, the governing equation as subsequent numerical simulation.
(3) this method uses the calculation method of bidirectional couple, on the rock matrix model and pore structure model of reconstruction
The solution of solid deformation and fluid field is carried out respectively, is then realized the transmitting of coupling variable by intermediate data platform and is solved and controls
System.
(4) it is solved in the form of non-coupled equation in this method, therefore the specific manifestation of coupling effect are as follows: temperature
Thermal strain caused by spending leads to the deformation of skeleton pattern and pore model grid respectively and repartitions；Framework deformation causes fluid
The change in channel, Fluid pressure lead to the deformation of skeleton again；The deformation of skeleton will affect the distribution situation in temperature field, pore-fluid
Flowing change the temperature at skeleton wall surface in the form of heat convection, to influence the heat transfer of skeleton.The mutual shadow of three
It rings and final balance is reached by iteration convergence continuous in numerical value calculates.
(5) due to the pore scale model used in the application, physical size is smaller, therefore the temperature on skeleton pattern
Degree gradient can be ignored, and the heat transfer between skeleton and fluid also reaches balance in the extremely short time, so this method is main
Study influence of the different temperatures to rock deformation and fluid behaviour.
(6) pore scale skeleton pattern and pore model based on reconstruction, in conjunction with heat-fluid-wall interaction mathematical modulo of building
The solution platform that ANASYS and CFX software is flowed as solid deformation and fluid is respectively adopted in type, and MPCCI is as coupling data
Transmitting platform carry out multi- scenarios method numerical simulation calculation.
Gu compared with existing heat-stream-multi-field coupling theory research, the beneficial effects of the present invention are:
This method is respectively formed two sets of independent grids of solid domain and fluid domain by building skeleton and hole grid model
System realizes the transmitting of coupled physic quantity by the coupling boundary of definition, discloses heat-stream-Gu Duochang by numerical simulation means
The microcosmic effect mechanism of coupling.
It is integrated in the advantage in ANASYS using CFX, the bidirectional couple calculating proposed in this method does not need third party's coupling
The auxiliary of software is closed, need to only can be realized by not having the MPCCI software for solving function as the transmitting platform of coupling data
The solution calculated is coupled, couple solution realizes that process is simple, and numerical simulation software is easy to operate.
Detailed description of the invention
In order to illustrate more clearly of the technical solution of the method for the present invention, come below in conjunction with the drawings and specific embodiments into one
Step illustrates the method for the present invention.
Fig. 1 is the theoretical basis concept map of the application method.
Fig. 2 is the application method implementation process diagram.
Fig. 3 is the original micro- CT image of rock core used in this method.
Fig. 4 is the pore scale grid model based on micro- CT image reconstruction, and (a) is solid skeletal model, (b) is fluid bore
Gap model.
Fig. 5 is that model boundary condition setting and coupling boundary define the setting of (a) skeleton pattern boundary；(b) pore model side
Boundary's setting.
Fig. 6 is the geometrical model of skeleton pattern and pore model.
Fig. 7 is changing rule of the model porosity with confining pressure.
Fig. 8 is changing rule of the model permeability with confining pressure.
Fig. 9 is model porosity variation with temperature rule.
Figure 10 is model permeability variation with temperature rule.
Figure 11 is changing rule of the model permeability with temperature and effective pressure, the direction (a) x；(b) direction y；(c) side z
To.
Figure 12 is changing rule of the model relative permeability with confining pressure.
Figure 13 is model relative permeability variation with temperature rule.
Specific embodiment
In order to realize the present invention basic theory, technological means and reach purpose ease of explanation, below in conjunction with this Shen
Please attached drawing in embodiment, and the specific case technical solutions in the embodiments of the present application that calculates carries out clearly and completely
Description.It needs to know, described embodiment is only some embodiments of the present application, rather than whole embodiments.It is based on
Embodiments herein, others skilled in the art obtain every other under the premise of other no innovative labors
The application protection scope all should belong in embodiment.
The method of the present invention utilizes the micro- CT image of rock, Gu rock fever-stream-multi- scenarios method is carried out based on minute yardstick reconstruction model
Numerical analysis method embodiment combination attached drawing details are as follows:
It is specific as described in Figure 1, Gu in rock fever-stream-multi-field coupling theory described herein, temperature field, stress field and
Interaction, interactional relationship between deformation field are as follows: influence of the fluid to temperature is mainly reflected on heat transfer, temperature pair
The change in flow field mainly passes through influence convective term；The change of temperature can cause the thermal strain of solid, and the deformation of solid will lead to again
The redistribution in temperature field；The deformation of solid can directly cause the change of seepage channel, Fluid pressure will affect again effective stress from
And change solid deformation amount.Therefore, the interaction between these three physical fields and influence control the Macroscopic physical spy of rock
Property, to decide engineering construction security development.
The specific implementation process of this method is as shown in Figure 2: mainly including pore scale Model Reconstruction, multi- scenarios method mathematical modulo
Three type building, numerical simulation calculation main steps.
S1: being illustrated in figure 3 the original micro- CT image data used in the present embodiment, and rock specimens as shown in the figure are artificial
Sandstone (MS1), remaining rock specimens also: carbonate rock (C1), Berea sandstone (B1) and loose sand (S6), the present invention are real
The core microcosmic CT experiment of example is applied mainly using Southwest Petrol University's hydrocarbon reservoirs and Development Engineering National Key Laboratory
Zeiss Xradia MicroXCT-400 microscopic imaging systems.
Using the micro- CT image data of the rock core of acquisition, pore scale structure is realized using the model reconstruction techniques of independent development
Change the building of grid model, in order to reduce the number of grid of reconstruction model, the generation and subsequent numerical value for improving Model Reconstruction are calculated
Rate of convergence, construct cube model as shown in Figure 4 (generally 4003~6003Pixel Dimensions).
S2: the method proposed in the application is different from traditional multi-field coupling theory research, is not by rock matrix and hole
Fluid is regarded as a continuous volume elements to characterize rock porous media, but framework skeleton pattern (Fig. 4 a) and pore model respectively
(Fig. 4 b) realizes the independent calculating of deformation field and flow field, then the analysis of multi- scenarios method is realized by the transmitting of coupling parameter, calculates
It is more accurate, also more meet reality.
Methods described herein be respectively adopted by the way of bidirectional couple ANASYS and CFX software solve solid deformation field and
Seepage field realizes the transmitting of intermediate coupling physical quantity using MPCCI software as data interchange platform.Each step of whole system
Rapid calculating is referred to as primary big iteration, and solid field and fluid field are iterated operation, single game operation respectively in big iteration each time
Convergence carries out coupled physic quantity transmitting again, until the calculating convergence of whole system.
S3: the numerical simulation calculation of fluid field uses CFX software in this method, therefore using based on N-S equation group
Governing equation describes the seepage flow of pore scale.By taking oil-water two-phase flow as an example, not miscible fluids continuity equation such as formula (1) institute
Show:
Wherein αiFor i in computing unitthThe volume fraction of item fluid, ρiWithIt is ithThe density and speed of item fluid, Sai
It can be the customized quality source item of user for source sink term.
It is controlled in entire fluid calculation domain using the single equation of momentum, as shown in formula (2).Velocity field shares in domain,
The form of expression is that each phase physical parameter (density and viscosity) takes volume average in control volume.
Wherein ρ and μ is respectively density and viscosity, p representative pressure gradient,Represent body force, pore scale model due to
Size is too small, can be ignored, itemLaplace pressure effect is described, value is only non-at fluid phase interface
0.It is associated with volume fraction every in domain, the volume of density and viscosity can averagely respectively indicate are as follows:
Wherein the viscosity variation with temperature of water can indicate are as follows:
The viscosity variation with temperature relational expression of oil is related to the component of oil, but generally referring to following relationship:
lgμ0=A+BT (6)
Wherein coefficient A and B can be fitted to obtain by the viscosity-temperature curve of oil.
Every shared energy equation in computational domain, as shown in formula (7):
WhereinFor i fluid effective thermal conductivities, ShFor heat radiation and the contribution of other volumetric sources, E
It is respectively the weighted average of each phase energy and temperature with T, is shown below respectively:
S4: rock deformation field in this method carries out numerical simulation using ANASYS software, rock using etc. moulded to ideal
Property strengthen model, the framework deformation field mathematical model of elastic stage is mainly by equilibrium equation.Geometric equation and physical equation are constituted,
Respectively such as formula (10), (11) and (12) are shown:
Wherein σijFor stress tensor, fiFor body force.
Wherein εx, εyAnd εzFor normal strain, γxy, γyzAnd γzxFor shearing strain, u and v and ω are displacement component
Wherein E and μ is respectively elasticity modulus and Poisson's ratio.
In deformation field, the heat transfer of skeleton follows Fourier heat equation:
Wherein qnIt " is heat flow density, k is thermal coefficient.
Thermal convection current phenomenon between rock matrix and pore-fluid then follows newton criterion:
qn"=h (Ts-Tf) (14)
Wherein h is convection transfer rate, TsFor skeleton surface temperature, TfFor surrounding fluid temperature.
S4: before carrying out numerical simulation calculation, need to carry out the setting of model boundary condition and the definition of coupling boundary.Gu
Body deformation solved by ANASYS, fluid flows through CFX and is solved, and based on Workbench platform into
The addition of row model boundary condition.As shown in Fig. 5 (a), four sides of skeleton pattern apply confining pressure respectively, and Fig. 5 (b) show fluid
Field upper and lower surface applies inlet and outlet pressure respectively, while the interface for defining solid and fluid is fluid-wall interaction interface.Fig. 6 institute
It is shown as the geometrical model file of the skeleton established in Workbench and hole.
Temperature gradient by moulded dimension very little in this present embodiment, therefore in skeleton pattern can be ignored；Simultaneously
Since the heat exchange between solid and fluid also reaches thermal balance in a very short period of time, so mainly being considered in embodiment different
At a temperature of the variation of pore structure caused by thermal strain and the influence to seepage flow, therefore steady temperature perimeter strip is arranged to skeleton pattern
Part.Coupling is calculated and is solved since seepage field, and the pore pressure being calculated is passed to solid field solver；Solid field meter
Obtained boundary node deflection passes to the wall surface of fluid field, so until system iteration convergence, completes to calculate.
S5: Gu heat-stream-multi- scenarios method mathematical model of model and building based on reconstruction, embodiment 1 pass through numerical simulation
Technique study stress and temperature develops on rock pore structure and absolute permeability influences.
The microcosmic rock mechanics parameters for wherein simulating use are measured by rock micron indentation test, the thermal coefficient of rock and
The coefficient of expansion uses empirical value, is as shown in table 1 the physical parameter of the different rock samples used in the present embodiment.
1. petrophysical parameter table of table
As shown in fig. 7, in the case where pore pressure is constant, model porosity increases with confining pressure and is reduced, when confining pressure increases
It is big to 30MPa when downward trend slow down because the plastic region inside rock sample increases.Fig. 8 show permeability with the change of confining pressure
Law becomes smaller as confining pressure increases permeability, and mainly seepage channel compressive deformation causes, furthermore with the increasing of plastic zone
Greatly, permeability decrease amplitude slows down, because rock sample is assumed to be isotropic hardening elastic-plastic material herein, plastic zone
Strain will not continue to increase.
Assuming that in the case where confining pressure and constant pore pressure, while the factors such as crack initiation caused by thermal strain are not considered
(20 DEG C~100 DEG C of set temperature range).It is respectively porosity and permeability variation with temperature relationship shown in Fig. 9 and Figure 10.
As can be seen that temperature increases, as skeleton expansion causes porosity to reduce, but the decline of this range of temperature internal porosity
Limited extent, only 2% or so.It can be seen that temperature to swollen by the change curve of different rock sample model permeabilities and temperature
Swollen coefficient is larger and the lower model influence of original permeability is bigger.Figure 11 (a), (b) and (c) can be seen that core permeability exists
Heterogeneity spatially, three directional permeabilitys illustrate with the situation of change of stress and temperature: pore-size is small, poor connectivity
Lower rock sample is higher to the sensibility of stress and temperature with benchmark permeability.
S6: embodiment 2 has studied the influence of stress and temperature for water drive oil effect.Select pore throat size larger, even
The general character is good, and the artificial sand rock sample MS1 that space transfusion rate heterogeneity is weaker.
Stress carries out the research that relative permeability influences under the hypothesis of steady temperature and pore pressure.In numerical simulation
The grease physical parameter of use is as shown in table 2.
2. grease physical parameter table of table
As shown in figure 12, in the case where pore pressure is certain, permeability saturation curve with confining pressure increase gradually to
Lower right is mobile, and isotonic point moves to right, and two-phase permeation area reduces, and water-oil phase permeability is on a declining curve, and residual oil is full
It is increased with degree.This is mainly due to increasing with model effective stress, pore structure compressive deformation, part narrowing throat is shunk very
To closure, model permeability decrease, the decline of grease fluid ability causes more residual oils can not be by displacement.
The research of Temperature Effect On Relative Permeability carries out under the hypothesis of constant confining pressure and pore pressure.In numerical simulation
The grease physical parameter of use is as shown in table 3.
3. grease physical parameter table of table
It simulates fluid media (medium) and uses water flooding and crude oil, the viscosity-temperature curve of water flooding uses the viscosity-temperature curve of water；Model is close
The viscosity coefficient that wall surface goes out is modified using miniflow boundary layer theory, i.e., there is given below the viscosity coefficient of water in model:
Wherein μwIt, can be with when ignoring temperature to influence for laminar boundary layer for dynamic viscosity of the water flooding in 273K
Using the function coefficient Φ not varied with temperature.
The rock sample buried depth selected in the present embodiment is shallower, and formation temperature is lower than 100 DEG C, therefore with common within 100 DEG C
Carry out analog study based on the viscosity-temperature curve of viscous crude sample, ignores influence of the temperature to miniflow boundary layer.Fluid temperature (F.T.) is by bone
Frame temperature determines, and corresponding grease viscosity coefficient is calculated.Foundation crude oil viscosity-temperature curve combination Bingham fluid fundamental formular,
Comprehensively consider miniflow effect of boundary layer, when Rheology Relationship in available microchannel:
Wherein τ0And μpLimiting shear stress and structural viscosity are respectively represented, can be fitted by crude oil rheological curve.
As shown in figure 13, as the temperature rises, curve right endpoint and isotonic point move to right, and oily phase fluid ability mentions
Height, relative permeability increase.Crude viscosity is increased with temperature and is reduced simultaneously, and grease mobile performance enhances, and fetters water saturation
Degree increases, and residual oil saturation reduces.Oil and water mobility ratio reduces, and weakens the fingering effect of water drive process, improves and finally adopt
Yield.Although the raising of temperature causes pore throat to shrink, the promotion of the reduction bring waterflooding effect of viscosity is apparently higher than heat
Strain bring negative effect.
The foregoing is merely presently preferred embodiments of the present invention, basic principle for describing the present invention, feature and main excellent
Point, is not intended to limit the invention, any modification made according to the technical essence of the invention to embodiment of above, etc.
With variation and modification, belong in the range of technical solution of the present invention.
Claims (9)
1. Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model, it is characterised in that:
(1) pore scale structured grid model is rebuild based on the micro- CT image of rock core；
(2) it is based on two-phase fluid mechanics and heat-THEORY OF FLUID-SOLID COUPLING FLOW THROUGH, Gu establish heat-stream-two phase fluid flow mathematical model；
(3) using ANASYS and CFX software as Coupled Numerical computing platform, realize that the numerical value of deformation domain and fluid domain is asked respectively
Solution；Data exchange and solver control are realized using MPCCI software as coupling data switching plane.
2. Gu according to heat-stream-multi- scenarios method analogy method described in claim 1 based on minute yardstick reconstruction model, feature exists
In: in the step (1), original CT image is handled by digital image processing techniques, respectively extract rock matrix and
Pore structure, and utilize the pore scale structured grid model of the Model Reconstruction algorithm of independent development building rock.
3. Gu according to heat-stream-multi- scenarios method analogy method described in claim 1 based on minute yardstick reconstruction model, feature exists
In: in the step (2), based on traditional permeation fluid mechanics and THEORY OF FLUID-SOLID COUPLING FLOW THROUGH, by introducing temperature term to fluid field and change
Shape field is modified, Gu building heat-stream-multi- scenarios method mathematical model.
4. Gu according to heat-stream-multi- scenarios method analogy method described in claim 1 based on minute yardstick reconstruction model, feature exists
In: in the step (3), based on Workbench platform, respectively to rock matrix grid model and hole grid model into
The addition of row bound condition and coupling boundary definition.
5. Gu according to heat-stream-multi- scenarios method analogy method described in claim 3 based on minute yardstick reconstruction model, feature exists
In: the mathematical model of building is realized in numerical simulation software using UDF (User Define Function) and is compiled.
6. Gu according to heat-stream-multi- scenarios method analogy method described in claim 1 based on minute yardstick reconstruction model, feature exists
In: in the step (3), rock matrix model is subjected to numerical simulation in ANASYS software, Gu while introducing heat-coupled mode
Solid deformation domain and solution and analysis that temperature influences are realized in the calculating of block.
7. Gu according to heat-stream-multi- scenarios method analogy method described in claim 1 based on minute yardstick reconstruction model, feature exists
In: in the step (3), pore model is subjected to numerical simulation in CFX software, while introducing the meter of heat-flow coupling module
It calculates, realizes fluid domain and solution and analysis that temperature influences.
8. Gu according to heat-stream-multi- scenarios method analogy method described in claim 1 based on minute yardstick reconstruction model, feature exists
In: in the step (3), using MPCCI software as data exchange and solver control platform, realize solid domain and fluid domain
The transmitting of calculating parameter.
9. Gu according to heat-stream-multi- scenarios method analogy method described in claim 8 based on minute yardstick reconstruction model, feature exists
In: by the post-processing of logarithm calculated result, analyzes stress condition and temperature parameter and rock pore structure is developed and transported
The influence of characteristic.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811092588.2A CN109063383A (en) | 2018-09-19 | 2018-09-19 | Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811092588.2A CN109063383A (en) | 2018-09-19 | 2018-09-19 | Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109063383A true CN109063383A (en) | 2018-12-21 |
Family
ID=64762244
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811092588.2A Pending CN109063383A (en) | 2018-09-19 | 2018-09-19 | Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109063383A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111179296A (en) * | 2020-01-08 | 2020-05-19 | 青岛科技大学 | Novel method for researching heat conduction characteristic of rock based on digital rock core technology |
| CN111948099A (en) * | 2020-08-13 | 2020-11-17 | 西南石油大学 | Testing device and method for evaluating backflow opportunity of fracturing fluid based on hydration |
| CN112049628A (en) * | 2020-10-20 | 2020-12-08 | 中国石油集团渤海钻探工程有限公司 | Method for judging and identifying oil-water interface of buried hill by using water content index |
| CN113109162A (en) * | 2021-04-13 | 2021-07-13 | 西南石油大学 | Rock fracture initiation pressure calculation method based on thermo-fluid-solid coupling |
| CN113189305A (en) * | 2021-04-29 | 2021-07-30 | 中国华能集团清洁能源技术研究院有限公司 | Simulation method and system for porous permeable rock |
| CN114861423A (en) * | 2022-04-22 | 2022-08-05 | 西北核技术研究所 | Method for simulating dynamic response numerical value of underground medium under action of thermal radiation |
| CN117408191A (en) * | 2023-12-15 | 2024-01-16 | 山东大学 | Fracture rock mass permeation-flow grouting simulation method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120221306A1 (en) * | 2009-04-08 | 2012-08-30 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
| CN104318032A (en) * | 2014-11-01 | 2015-01-28 | 西南石油大学 | Method for calculating oil field casing damage under fluid-solid coupling effect |
| CN104345133A (en) * | 2014-09-22 | 2015-02-11 | 西南石油大学 | Numerical analysis method applied to leaking and channeling mechanism and prevention and control method of injection fluid |
| CN105393110A (en) * | 2013-08-06 | 2016-03-09 | Bp北美公司 | Image-based direct numerical simulation of petrophysical properties under simulated stress and strain conditions |
-
2018
- 2018-09-19 CN CN201811092588.2A patent/CN109063383A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120221306A1 (en) * | 2009-04-08 | 2012-08-30 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
| CN105393110A (en) * | 2013-08-06 | 2016-03-09 | Bp北美公司 | Image-based direct numerical simulation of petrophysical properties under simulated stress and strain conditions |
| CN104345133A (en) * | 2014-09-22 | 2015-02-11 | 西南石油大学 | Numerical analysis method applied to leaking and channeling mechanism and prevention and control method of injection fluid |
| CN104318032A (en) * | 2014-11-01 | 2015-01-28 | 西南石油大学 | Method for calculating oil field casing damage under fluid-solid coupling effect |
Non-Patent Citations (3)
| Title |
|---|
| RUISONG等: "A new method to reconstruct structured mesh model from micro-computed tomography images of porous media and its application", 《INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER》 * |
| 宋睿: "基于微尺度重建模型的岩石热—流—固耦合细观机理研究", 《中国优秀博硕士学位论文全文数据库(博士)(工程科技Ⅰ辑)》 * |
| 邓英尔等: "低渗透微尺度孔隙气体渗流规律", 《力学与实践》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111179296A (en) * | 2020-01-08 | 2020-05-19 | 青岛科技大学 | Novel method for researching heat conduction characteristic of rock based on digital rock core technology |
| CN111948099A (en) * | 2020-08-13 | 2020-11-17 | 西南石油大学 | Testing device and method for evaluating backflow opportunity of fracturing fluid based on hydration |
| CN111948099B (en) * | 2020-08-13 | 2022-02-08 | 西南石油大学 | Testing device and method for evaluating backflow opportunity of fracturing fluid based on hydration |
| CN112049628A (en) * | 2020-10-20 | 2020-12-08 | 中国石油集团渤海钻探工程有限公司 | Method for judging and identifying oil-water interface of buried hill by using water content index |
| CN113109162A (en) * | 2021-04-13 | 2021-07-13 | 西南石油大学 | Rock fracture initiation pressure calculation method based on thermo-fluid-solid coupling |
| CN113189305A (en) * | 2021-04-29 | 2021-07-30 | 中国华能集团清洁能源技术研究院有限公司 | Simulation method and system for porous permeable rock |
| CN114861423A (en) * | 2022-04-22 | 2022-08-05 | 西北核技术研究所 | Method for simulating dynamic response numerical value of underground medium under action of thermal radiation |
| CN117408191A (en) * | 2023-12-15 | 2024-01-16 | 山东大学 | Fracture rock mass permeation-flow grouting simulation method and system |
| CN117408191B (en) * | 2023-12-15 | 2024-04-02 | 山东大学 | Fracture rock mass permeation-flow grouting simulation method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109063383A (en) | Gu heat-stream-multi- scenarios method analogy method based on minute yardstick reconstruction model | |
| CN113826099B (en) | Flow simulation and transient well analysis method based on generalized pipe flow seepage coupling | |
| Schaap et al. | Comparison of pressure‐saturation characteristics derived from computed tomography and lattice Boltzmann simulations | |
| Xiong et al. | Experimental and numerical investigation on nonlinear flow behaviour through three dimensional fracture intersections and fracture networks | |
| Lacey et al. | Effects of pore and grain size on water and polymer flooding in micromodels | |
| Deng et al. | Simulation of grouting process in rock masses under a dam foundation characterized by a 3D fracture network | |
| Yan et al. | An enhanced discrete fracture network model for multiphase flow in fractured reservoirs | |
| Ma et al. | Numerical investigation of heat transfer performance of water flowing through a reservoir with two intersecting fractures | |
| CN104298797A (en) | Method for determining remaining oil of fracture and hole type oil reservoir under encirclement of high flow guide channel | |
| Elfeel et al. | Multi-scale simulation of WAG flooding in naturally fractured reservoirs | |
| Hussain et al. | Multiphase fluid flow through fractured porous media supported by innovative laboratory and numerical methods for estimating relative permeability | |
| Shimizu | Microscopic numerical model of fluid flow in granular material | |
| Mei et al. | Fractal analysis of shape factor for matrix-fracture transfer function in fractured reservoirs | |
| Zhao et al. | A dynamic pore network model for imbibition simulation considering corner film flow | |
| Wang et al. | Theoretical analysis and semi-analytical formulation for capturing the coupled thermal-hydraulic-mechanical process using the stress formulation | |
| Wolff et al. | TREATMENT OF TENSORIAL RELATIVE PERMEABILITIES WITH MULTIPOINT FLUX APPROXIMATION. | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| Liu et al. | Impact of CO2 injection rate on heat extraction at the HDR geothermal field of Zhacanggou, China | |
| Heeremans et al. | Feasibility study of WAG injection in naturally fractured reservoirs | |
| Guo et al. | Vertically integrated dual-porosity and dual-permeability models for CO2 sequestration in fractured geological formation | |
| Wang et al. | Simulation of heat and fluid flow in porous medium and fractures by material point method | |
| Huang et al. | Water Invasion Into Multi-Layer and Multi-Pressure Carbonate Reservoir: A Pore-Scale Simulation | |
| Lei et al. | A dynamic discrete fracture model for fluid flow in fractured low-permeability reservoirs | |
| Al-Husseini et al. | Hydrodynamics Behaviour of Single and Multi Fracture with Different Orientations in Petroleum Reservoir | |
| Tachibana et al. | Characterization of transition from Darcy to non-Darcy flow with 3D pore-level simulations |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181221 |
|
| WD01 | Invention patent application deemed withdrawn after publication |