CN105089595A - Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion - Google Patents
Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion Download PDFInfo
- Publication number
- CN105089595A CN105089595A CN201510280351.7A CN201510280351A CN105089595A CN 105089595 A CN105089595 A CN 105089595A CN 201510280351 A CN201510280351 A CN 201510280351A CN 105089595 A CN105089595 A CN 105089595A
- Authority
- CN
- China
- Prior art keywords
- fracture
- well
- fluid
- model
- wellbore
- 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
Abstract
The invention provides an oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion. The reservoir numerical simulation method comprises the following steps: obtaining dynamic and static data of a well and an oil reservoir, and determining fracture parameters; establishing a flow equation of fluid in a horizontal fracture according to fluid characteristics and diversion capacity in the fracture parameters; establishing a wellbore flow equation according to wellbore pipe-flow characteristics; establishing an extension well model according to a well track, a well completion position, the fluid flow equation, the wellbore flow equation and the fracture parameters; establishing a geologic model according to static parameters of the oil reservoir, and initializing the geologic model according to initial formation pressure distribution or the relationship between pressure and depth and attributes of rocks and the fluid; and coupling the extension well model and the initialized geological model, and solving a coupling result. The reservoir numerical simulation method realizes the high-precision numerical simulation of the hydraulic fracturing well and the oil reservoir, provides the basis for scientific use of downhole flow control equipment and achieves the purpose of increasing a recovery ratio.
Description
Technical field
The present invention relates to Research Numerical Simulation Techique, particularly relate to the numerical reservoir simulation method under the effect of a kind of horizontal fractures fracture guide and device.
Background technology
Middle-Low Permeability Reservoirs has become the main potential of ultra-high water cut stage multilayer heterogeneous reservoir, as grand celebration Xingshugang oilfield, Dongying Tuo21 fault block and Pucheng oil field in Zhongyuan etc.When adopting the Non-primary Trunk Die Cavity of straight well exploitation poor properties, usually need fracturing to increase production to obtain economic productivity, the crack of formation mainly comprise perpendicular to pit shaft horizontal joint and be parallel to the vertical lap seam of pit shaft.Easily horizontal fracture is formed during shallow at buried depth, oil reservoir fracture pressure gradient is little Middle-Low Permeability Reservoirs pressure break, as Daqing oil field, Yumen Oilfield, Karamay oilfield etc., this is very large for water flooding effectiveness impact, thus needs a kind of accurate numerical reservoir simulation method badly.
For above problem, domestic and international many scholars have carried out large quantity research for the production capacity of vertical lap seam, as M.PRATS (M.Prats, MemberAime, P.hazebroek, W.R.Strickler.EffectofVerticalFracturesonReservoirBehavi or-Compressible-FluidCase [C] .SPE98, 1962.), KyleE.Friehasuf, AjaySuri (KyleE.Friehauf, AjaySuri, MukulM.Sharma.ASimpleandAccurateModelforWellProductivity forHydraulicallyFracturedWells [C] .SPEProduction & Operations, 2010, 25 (4): 453-460.), Deng Yinger and Liu Ciqun (Deng Yinger, Liu Ciqun. vertically fractured well low-permeability reservoir pressure Analysis for Non-Linear Percolation [J]. Petroleum finance, 2003, 30 (1): 81-83.) etc., but it is little for the Study on Productivity of horizontal joint, mainly concentrate on Study on Productivity aspect, method as W.Sung numerical simulation in 1987 compared for the effect of increasing production (W.Sung of horizontal fracture and vertical fracture, T.Ertekin.PerformanceComparisonofVerticalandHorizontalHy draulicFracturesandHorizontalBoreholesinLowPermeabilityR eservoirs:ANumericalStudy [C] .SPE16407, 1987.), the subject matter existed comprises: 1. pressure break reservoir numerical simulation research is main concentrates vertical lap seam, and business-like numerical simulation software can only consider vertical fracture, 2. the research for horizontal fracture is all simple study mechanism, can not hide numerical simulation to carrying out after the whole district's pressure break Large Oil, 3. very low for flow simulating precision in horizontal fracture, normally direct is high permeability band by crack treatment, and the flowing in crack is still considered as Darcy Flow, 4. the multiphase pipe flow in pit shaft is not considered.
Summary of the invention
Embodiments provide the numerical reservoir simulation method under the effect of a kind of horizontal fractures fracture guide, to realize carrying out high resolution numerical simulation to hydraulically fractured wells and oil reservoir, for the use downhole flow control appliance of science provides foundation, reach the object improving recovery ratio.
To achieve these goals, the embodiment of the present invention provides the numerical reservoir simulation method under the effect of a kind of horizontal fractures fracture guide, and described numerical reservoir simulation method comprises:
Obtain the dynamic static data of well and oil reservoir, determine fracture parameters;
The liquid flow equation in horizontal fracture is set up according to the flow conductivity in characteristic of fluid and described fracture parameters;
Wellbore Flow equation is set up according to wellbore tubular stream feature;
Expansion well model is set up according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters;
Static parameter according to oil reservoir sets up geological model, and according to original formation pressure distribution or the relation of pressure and the degree of depth, and rock, fluid properties initialize described geological model;
Described expansion well model and the described geological model after initializing are coupled and solve coupling result.
In one embodiment, describedly determine fracture parameters, comprise: determine described fracture parameters according to the well test analysis data in described dynamic static data and microseism decryption, wherein, described fracture parameters comprises: fracture support width, fracture spacing, fracture half-length, fracture condudtiviy, the extension feature in crack and distribution layer position.
In one embodiment, set up the liquid flow equation in horizontal fracture according to the flow conductivity in characteristic of fluid and described fracture parameters, comprising:
According to the characteristic of fluid in described dynamic static data and described fracture condudtiviy, Reynolds number flow pattern judgment criterion is utilized to determine the crack inner fluid type of flow;
The liquid flow equation in horizontal fracture is set up according to the described crack inner fluid type of flow;
Calculate according to described liquid flow equation and split pressure distribution in fracture and fluid flowing.
In one embodiment, calculate according to described liquid flow equation and split pressure distribution in fracture and fluid flowing, comprising:
Utilize binomial high speed non linear fluid flow through porous medium model to calculate and split pressure distribution in fracture and fluid flowing, ignore the flowing in z direction, described binomial high speed non linear fluid flow through porous medium model is as follows:
Wherein, x, y represent x direction and y direction respectively, and μ is fluid viscosity, k
ffor fracture permeabgility, ρ is fluid density, and v is seepage velocity, and ζ is high speed non linear fluid flow through porous medium coefficient.
In one embodiment, set up Wellbore Flow equation according to wellbore tubular stream feature, comprising:
M phase Wellbore Flow equation is set up according to wellbore tubular stream feature:
According to described Wellbore Flow equation calculating pressure distribution;
Wherein, ρ
mgsin (θ) is hydrostatic power,
for friction pressure loss,
for the acceleration pressure loss, g is acceleration of gravity; θ is the angle between pit shaft and stratum, and A is that wellbore section amasss, and D is mineshaft diameter, ρ
mfor m phase fluid density, f
tpfor friction pressure loss coefficient, V
mfor m phase unit volume.
In one embodiment, set up expansion well model according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters, comprising:
According to flaw size, crack extension feature and divide crack grid with the way of contact of pit shaft, and subdivision is carried out to described crack grid;
Liquid flow equation discretization in the grid of described crack is utilized to form numerical model;
By fracture number value model be described expansion well model based on the well Model coupling of described Wellbore Flow equation.
In one embodiment, described expansion well model and the described geological model after initializing are coupled and solve coupling result, comprising:
The Nonlinear System of Equations with sparse coefficient matrix is set up based on described expansion well model and the described geological model after initializing;
Adopt Nonlinear System of Equations described in decoupling zero fully implicit solution Algorithm for Solving.
In one embodiment, described dynamic static packet is drawn together: the well location coordinate of every mouthful of well, completion mode, perforating depth, well head, well depth, hole deviation and well track data, the structural configuration of oil reservoir target zone position, log analysis data, seismic interpretation data, layer data, the degree of porosity of target zone position, permeability, oil saturation and original formation pressure parameter, high pressure property data, characteristic of fluid data, the wellbore tubular stream feature of rock and fluid in stratum, individual well Production development data, well test analysis data and microseism decryption.
To achieve these goals, the embodiment of the present invention additionally provides the reservoir numerical simulation device under the effect of a kind of horizontal fractures fracture guide, and described reservoir numerical simulation device comprises:
Fracture parameters determining unit, for obtaining the dynamic static data of well and oil reservoir, determines fracture parameters;
Liquid flow equation sets up unit, for setting up the liquid flow equation in horizontal fracture according to the flow conductivity in characteristic of fluid and described fracture parameters;
Wellbore Flow establishing equation unit, for setting up Wellbore Flow equation according to wellbore tubular stream feature;
Unit set up by expansion well model, for setting up expansion well model according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters;
Geological model generation unit, sets up geological model for the static parameter according to oil reservoir, and according to original formation pressure distribution or the relation of pressure and the degree of depth, and rock, fluid properties initialize described geological model;
Coupling unit, for being coupled described expansion well model and the described geological model after initializing and solving coupling result.
In one embodiment, described fracture parameters determining unit specifically for: determine described fracture parameters according to the well test analysis data in described dynamic static data and microseism decryption, wherein, described fracture parameters comprises: fracture support width, fracture spacing, fracture half-length, fracture condudtiviy, the extension feature in crack and distribution layer position.
In one embodiment, described liquid flow equation is set up unit and is comprised:
Type of flow determination module, for according to the characteristic of fluid in described dynamic static data and described fracture condudtiviy, utilizes Reynolds number flow pattern judgment criterion to determine the crack inner fluid type of flow;
Liquid flow equation sets up module, for setting up the liquid flow equation in horizontal fracture according to the described crack inner fluid type of flow;
First computing module, splits pressure distribution in fracture and fluid flowing for calculating according to described liquid flow equation.
In one embodiment, described first computing module specifically for:
Utilize binomial high speed non linear fluid flow through porous medium model to calculate and split pressure distribution in fracture and fluid flowing, ignore the flowing in z direction, described binomial high speed non linear fluid flow through porous medium model is as follows:
Wherein, x, y represent x direction and y direction respectively, and μ is fluid viscosity, k
ffor fracture permeabgility, ρ is fluid density, and v is seepage velocity, and ζ is high speed non linear fluid flow through porous medium coefficient.
In one embodiment, described Wellbore Flow establishing equation unit comprises:
Establishing equation module, for setting up m phase Wellbore Flow equation according to wellbore tubular stream feature:
Second computing module, for distributing according to described Wellbore Flow equation calculating pressure;
Wherein, ρ
mgsin (θ) is hydrostatic power,
for friction pressure loss,
for the acceleration pressure loss, g is acceleration of gravity; θ is the angle between pit shaft and stratum, and A is that wellbore section amasss, and D is mineshaft diameter, ρ
mfor m phase fluid density, f
tpfor friction pressure loss coefficient, V
mfor m phase unit volume.
In one embodiment, described expansion well model is set up unit and is comprised:
Stress and strain model module, for according to flaw size, crack extension feature and divide crack grid with the way of contact of pit shaft, and carries out subdivision to described crack grid;
Descretization module, forms numerical model for utilizing liquid flow equation discretization in the grid of described crack;
Coupling module, for by fracture number value model be described expansion well model based on the well Model coupling of described Wellbore Flow equation.
In one embodiment, described coupling unit comprises:
Equation group sets up module, for setting up the Nonlinear System of Equations with sparse coefficient matrix based on described expansion well model and the described geological model after initializing;
Equation solution module, for adopting Nonlinear System of Equations described in decoupling zero fully implicit solution Algorithm for Solving.
In one embodiment, described dynamic static packet is drawn together: the well location coordinate of every mouthful of well, completion mode, perforating depth, well head, well depth, hole deviation and well track data, the structural configuration of oil reservoir target zone position, log analysis data, seismic interpretation data, layer data, the degree of porosity of target zone position, permeability, oil saturation and original formation pressure parameter, high pressure property data, characteristic of fluid data, the wellbore tubular stream feature of rock and fluid in stratum, individual well Production development data, well test analysis data and microseism decryption.
Present invention achieves and high resolution numerical simulation is carried out to hydraulically fractured wells and oil reservoir, for the use downhole flow control appliance of science provides foundation, reach the object improving recovery ratio.Particularly, the horizontal fracture formed after the present invention is directed to low permeability formation fracturing, propose the method for numerical simulation based on expansion well model and reservoir model coupling, can non-darcy Multiphase Flow in dummy level crack, consider that the pitometer in pit shaft calculates pressure distribution, improve the simulation precision of well in Research Numerical Simulation Techique.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the numerical reservoir simulation method flow chart of one embodiment of the invention;
Fig. 2 is the liquid flow equation method for building up flow chart of the embodiment of the present invention;
Fig. 3 is the expansion well method for establishing model flow chart of the embodiment of the present invention;
Fig. 4 is the expansion well model schematic of the embodiment of the present invention;
Fig. 5 is the coupling result method for solving flow chart of the embodiment of the present invention;
Fig. 6 is the numerical reservoir simulation method flow chart of another embodiment of the present invention;
Fig. 7 is the structured flowchart of the reservoir numerical simulation device of the embodiment of the present invention;
Fig. 8 is the structured flowchart that the liquid flow equation of the embodiment of the present invention sets up unit 702;
Fig. 9 is the structured flowchart of the Wellbore Flow establishing equation unit 703 of the embodiment of the present invention;
Figure 10 is the structured flowchart that unit 704 set up by the expansion well model of the embodiment of the present invention;
Figure 11 is the structured flowchart of the coupling unit 706 of the embodiment of the present invention.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
In order to solve problems of the prior art, embodiments provide the numerical reservoir simulation method under the effect of a kind of horizontal fractures fracture guide, as shown in Figure 1, described numerical reservoir simulation method comprises:
Step 101: the dynamic static data obtaining well and oil reservoir, determines fracture parameters;
Step 102: set up the liquid flow equation in horizontal fracture according to the flow conductivity in characteristic of fluid and described fracture parameters;
Step 103: set up Wellbore Flow equation according to wellbore tubular stream feature;
Step 104: set up expansion well model according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters;
Step 105: the static parameter according to oil reservoir sets up geological model, and according to original formation pressure distribution or the relation of pressure and the degree of depth, and rock, fluid properties initialize described geological model;
Step 106: described expansion well model and the described geological model after initializing are coupled and solve coupling result.
Flow process is as shown in Figure 1 known, in the numerical reservoir simulation method of the embodiment of the present invention, sets up the liquid flow equation in horizontal fracture, Wellbore Flow equation and expansion well model first respectively, then sets up geological model of oil accumulation and numerical simulator
Described dynamic static packet in step 101 is drawn together: the well location coordinate of every mouthful of well, completion mode, perforating depth, well head, well depth, hole deviation and well track data, structural configuration (the end face of oil reservoir target zone position, bottom surface constructional drawing), log analysis data, seismic interpretation data, layer data, the degree of porosity of target zone position, permeability, oil saturation and original formation pressure parameter, high-voltage physics data (the phase percolation curve of rock and fluid in stratum, capillary pressure curve, PTV data), characteristic of fluid data, wellbore tubular stream feature, individual well Production development data (oil production, aquifer yield, gas production, gas-oil ratio, moisture content, oil pressure, casing pressure, flowing bottomhole pressure (FBHP), injection allocation amount), well test analysis data and microseism decryption.
In step 101, determine that fracture parameters specifically comprises: determine described fracture parameters according to the well test analysis data in described dynamic static data and microseism decryption, wherein, described fracture parameters comprises: fracture support width, fracture spacing, fracture half-length, fracture condudtiviy, the extension feature in crack and distribution layer position.
In one embodiment, as shown in Figure 2, set up the liquid flow equation in horizontal fracture in step 102 according to the flow conductivity in characteristic of fluid and described fracture parameters, comprise the steps:
Step 201: according to the characteristic of fluid in described dynamic static data and described fracture condudtiviy, utilizes Reynolds number flow pattern judgment criterion to determine the crack inner fluid type of flow.
Seepage flow results from capillary tube, and in reservoir, fluidised form is generally seepage flow, and due to physical properties of fluids (as viscosity), reservoir properties difference, rate of flow of fluid is different with pressure dependence, and darcy is linear, and non-darcy is non-linear.
Step 202: set up the liquid flow equation in horizontal fracture according to the described crack inner fluid type of flow;
Step 203: calculate according to described liquid flow equation and split pressure distribution in fracture and fluid flowing.
In step 202 and step 203; non-Darcy flow can be caused usually in high flow conductivity crack; binomial high speed non linear fluid flow through porous medium model can be utilized to calculate and split pressure distribution in fracture and fluid flowing, ignore the flowing in z direction, described binomial high speed non linear fluid flow through porous medium model is as follows:
Wherein, x, y represent x direction and y direction respectively, and μ is fluid viscosity, k
ffor fracture permeabgility, ρ is fluid density, and v is seepage velocity, and ζ is high speed non linear fluid flow through porous medium coefficient.
In one embodiment, set up Wellbore Flow equation according to wellbore tubular stream feature in step 103, specifically comprise:
M phase Wellbore Flow equation is set up according to wellbore tubular stream feature:
According to described Wellbore Flow equation calculating pressure distribution;
Wherein, ρ
mgsin (θ) is hydrostatic power,
for friction pressure loss,
for the acceleration pressure loss, g is acceleration of gravity; θ is the angle between pit shaft and stratum, and A is that wellbore section amasss, and D is mineshaft diameter, ρ
mfor m phase fluid density, f
tpfor friction pressure loss coefficient, V
mfor m phase unit volume.
As shown in Figure 3, in an embodiment, set up expansion well model according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters in step 104, specifically comprise:
Step 301: according to flaw size, crack extension feature and divide crack grid with the way of contact of pit shaft, and subdivision is carried out to described crack grid, crack grid much smaller than oil reservoir grid, as shown in Figure 4.In Fig. 4,1 is pit shaft, and 2 is expansion well area coarsening grids, and 3 is oil reservoir area grid.
Step 302: utilize liquid flow equation discretization in the grid of described crack to form numerical model.
Step 303: by fracture number value model be described expansion well model based on the well Model coupling of described Wellbore Flow equation.
In step 105, need to arrange nonlinear analysis controling parameters: arrange computational accuracy, iterative number of times initializes and calculate data (initializing the three-dimensional pressure of geological model and saturation field data and time data).
In one embodiment, when specifically implementing, step 106 comprises:
Step 501: set up the Nonlinear System of Equations with sparse coefficient matrix based on described expansion well model and the described geological model after initializing;
Step 502: adopt Nonlinear System of Equations described in decoupling zero fully implicit solution Algorithm for Solving.
During concrete enforcement, to time parameter upgrade and initialize assignment as follows: calculated sometime step after, each grid that the pressure that this time step is calculated and saturation data assignment walk to future time, as future time walk initial value.
Fig. 1, Fig. 2, Fig. 3 and Fig. 5 can be summarized by the numerical reservoir simulation method flow chart shown in Fig. 6, and the method flow in Fig. 6 is divided into two, left and right, has been described in detail, does not repeat them here in figure.
Present invention achieves and high resolution numerical simulation is carried out to hydraulically fractured wells and oil reservoir, for the use downhole flow control appliance of science provides foundation, reach the object improving recovery ratio.Particularly, the horizontal fracture formed after the present invention is directed to low permeability formation fracturing, propose the method for numerical simulation based on expansion well model and reservoir model coupling, can non-darcy Multiphase Flow in dummy level crack, consider that the pitometer in pit shaft calculates pressure distribution, improve the simulation precision of well in Research Numerical Simulation Techique.
Embodiments provide the reservoir numerical simulation device under the effect of a kind of horizontal fractures fracture guide, as shown in Figure 7, described reservoir numerical simulation device comprises: fracture parameters determining unit 701, liquid flow equation sets up unit 702, Wellbore Flow establishing equation unit 703, expansion well model sets up unit 704, geological model generation unit 705 and coupling unit 706.
Fracture parameters determining unit 701, for obtaining the dynamic static data of well and oil reservoir, determines fracture parameters;
Liquid flow equation sets up unit 702 for setting up the liquid flow equation in horizontal fracture according to the flow conductivity in characteristic of fluid and described fracture parameters;
Wellbore Flow establishing equation unit 703 is for setting up Wellbore Flow equation according to wellbore tubular stream feature;
Expansion well model sets up unit 704 for setting up expansion well model according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters;
Geological model generation unit 705 sets up geological model for the static parameter according to oil reservoir, and according to original formation pressure distribution or the relation of pressure and the degree of depth, and rock, fluid properties initialize described geological model;
Coupling unit 706 is for being coupled described expansion well model and the described geological model after initializing and solving coupling result.
In one embodiment, described dynamic static packet is drawn together: the well location coordinate of every mouthful of well, completion mode, perforating depth, well head, well depth, hole deviation and well track data, structural configuration (the end face of oil reservoir target zone position, bottom surface constructional drawing), log analysis data, seismic interpretation data, layer data, the degree of porosity of target zone position, permeability, oil saturation and original formation pressure parameter, high-voltage physics data (the phase percolation curve of rock and fluid in stratum, capillary pressure curve, PTV data), characteristic of fluid data, wellbore tubular stream feature, individual well Production development data (oil production, aquifer yield, gas production, gas-oil ratio, moisture content, oil pressure, casing pressure, flowing bottomhole pressure (FBHP), injection allocation amount), well test analysis data and microseism decryption.
In one embodiment, fracture parameters determining unit 701 specifically for: determine described fracture parameters according to the well test analysis data in described dynamic static data and microseism decryption, wherein, described fracture parameters comprises: fracture support width, fracture spacing, fracture half-length, fracture condudtiviy, the extension feature in crack and distribution layer position.
In one embodiment, as shown in Figure 8, liquid flow equation is set up unit 702 and is comprised: type of flow determination module 801, and liquid flow equation sets up module 802 and the first computing module 803.
Type of flow determination module 801, for according to the characteristic of fluid in described dynamic static data and described fracture condudtiviy, utilizes Reynolds number flow pattern judgment criterion to determine the crack inner fluid type of flow;
Liquid flow equation sets up module 802 for setting up the liquid flow equation in horizontal fracture according to the described crack inner fluid type of flow;
First computing module 803 splits pressure distribution in fracture and fluid flowing for calculating according to described liquid flow equation.
In one embodiment, the first computing module 803 specifically for:
Utilize binomial high speed non linear fluid flow through porous medium model to calculate and split pressure distribution in fracture and fluid flowing, ignore the flowing in z direction, described binomial high speed non linear fluid flow through porous medium model is as follows:
Wherein, x, y represent x direction and y direction respectively, and μ is fluid viscosity, k
ffor fracture permeabgility, ρ is fluid density, and v is seepage velocity, and ζ is high speed non linear fluid flow through porous medium coefficient.
In one embodiment, as shown in Figure 9, Wellbore Flow establishing equation unit 703 comprises: establishing equation module 901 and the second computing module 902.
Establishing equation module 901 is for setting up m phase Wellbore Flow equation according to wellbore tubular stream feature:
Second computing module 902 is for distributing according to described Wellbore Flow equation calculating pressure;
Wherein, ρ
mgsin (θ) is hydrostatic power,
for friction pressure loss,
for the acceleration pressure loss, g is acceleration of gravity; θ is the angle between pit shaft and stratum, and A is that wellbore section amasss, and D is mineshaft diameter, ρ
mfor m phase fluid density, f
tpfor friction pressure loss coefficient, V
mfor m phase unit volume.
In one embodiment, as shown in Figure 10, described expansion well model is set up unit 704 and is comprised: stress and strain model module 1001, descretization module 1002 and coupling module 1003.
Stress and strain model module 1001 for according to flaw size, crack extension feature and divide crack grid with the way of contact of pit shaft, and carries out subdivision to described crack grid;
Descretization module 1002 forms numerical model for utilizing liquid flow equation discretization in the grid of described crack;
Coupling module 1003 for by fracture number value model be described expansion well model based on the well Model coupling of described Wellbore Flow equation.
In one embodiment, as shown in figure 11, described coupling unit 706 comprises: equation group sets up module 1101 and equation solution module 1102.
Establishing equation module 1101 is for setting up the Nonlinear System of Equations with sparse coefficient matrix based on described expansion well model and the described geological model after initializing;
Equation solution module 1102 is for adopting Nonlinear System of Equations described in decoupling zero fully implicit solution Algorithm for Solving.
Present invention achieves and high resolution numerical simulation is carried out to hydraulically fractured wells and oil reservoir, for the use downhole flow control appliance of science provides foundation, reach the object improving recovery ratio.Particularly, the horizontal fracture formed after the present invention is directed to low permeability formation fracturing, propose the method for numerical simulation based on expansion well model and reservoir model coupling, can non-darcy Multiphase Flow in dummy level crack, consider that the pitometer in pit shaft calculates pressure distribution, improve the simulation precision of well in Research Numerical Simulation Techique.
Those skilled in the art should understand, embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt the form of complete hardware embodiment, completely software implementation or the embodiment in conjunction with software and hardware aspect.And the present invention can adopt in one or more form wherein including the upper computer program implemented of computer-usable storage medium (including but not limited to magnetic disc store, CD-ROM, optical memory etc.) of computer usable program code.
The present invention describes with reference to according to the flow chart of the method for the embodiment of the present invention, equipment (system) and computer program and/or block diagram.Should understand can by the combination of the flow process in each flow process in computer program instructions realization flow figure and/or block diagram and/or square frame and flow chart and/or block diagram and/or square frame.These computer program instructions can being provided to the processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing device to produce a machine, making the instruction performed by the processor of computer or other programmable data processing device produce device for realizing the function of specifying in flow chart flow process or multiple flow process and/or block diagram square frame or multiple square frame.
These computer program instructions also can be stored in can in the computer-readable memory that works in a specific way of vectoring computer or other programmable data processing device, the instruction making to be stored in this computer-readable memory produces the manufacture comprising command device, and this command device realizes the function of specifying in flow chart flow process or multiple flow process and/or block diagram square frame or multiple square frame.
These computer program instructions also can be loaded in computer or other programmable data processing device, make on computer or other programmable devices, to perform sequence of operations step to produce computer implemented process, thus the instruction performed on computer or other programmable devices is provided for the step realizing the function of specifying in flow chart flow process or multiple flow process and/or block diagram square frame or multiple square frame.
Apply specific embodiment in the present invention to set forth principle of the present invention and embodiment, the explanation of above embodiment just understands method of the present invention and core concept thereof for helping; Meanwhile, for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention.
Claims (16)
1. the numerical reservoir simulation method under the effect of horizontal fractures fracture guide, is characterized in that, described numerical reservoir simulation method comprises:
Obtain the dynamic static data of well and oil reservoir, determine fracture parameters;
The liquid flow equation in horizontal fracture is set up according to the flow conductivity in characteristic of fluid and described fracture parameters;
Wellbore Flow equation is set up according to wellbore tubular stream feature;
Expansion well model is set up according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters;
Static parameter according to oil reservoir sets up geological model, and according to original formation pressure distribution or the relation of pressure and the degree of depth, and rock, fluid properties initialize described geological model;
Described expansion well model and the described geological model after initializing are coupled and solve coupling result.
2. numerical reservoir simulation method according to claim 1, it is characterized in that, describedly determine fracture parameters, comprise: determine described fracture parameters according to the well test analysis data in described dynamic static data and microseism decryption, wherein, described fracture parameters comprises: fracture support width, fracture spacing, fracture half-length, fracture condudtiviy, the extension feature in crack and distribution layer position.
3. numerical reservoir simulation method according to claim 2, is characterized in that, sets up the liquid flow equation in horizontal fracture, comprising according to the flow conductivity in characteristic of fluid and described fracture parameters:
According to the characteristic of fluid in described dynamic static data and described fracture condudtiviy, Reynolds number flow pattern judgment criterion is utilized to determine the crack inner fluid type of flow;
The liquid flow equation in horizontal fracture is set up according to the described crack inner fluid type of flow;
Calculate according to described liquid flow equation and split pressure distribution in fracture and fluid flowing.
4. numerical reservoir simulation method according to claim 3, is characterized in that, calculates and splits pressure distribution in fracture and fluid flowing, comprising according to described liquid flow equation:
Utilize binomial high speed non linear fluid flow through porous medium model to calculate and split pressure distribution in fracture and fluid flowing, ignore the flowing in z direction, described binomial high speed non linear fluid flow through porous medium model is as follows:
Wherein, x, y represent x direction and y direction respectively, and μ is fluid viscosity, k
ffor fracture permeabgility, ρ is fluid density, and v is seepage velocity, and ζ is high speed non linear fluid flow through porous medium coefficient.
5. numerical reservoir simulation method according to claim 4, is characterized in that, sets up Wellbore Flow equation, comprising according to wellbore tubular stream feature:
M phase Wellbore Flow equation is set up according to wellbore tubular stream feature:
According to described Wellbore Flow equation calculating pressure distribution;
Wherein, ρ
mgsin (θ) is hydrostatic power,
for friction pressure loss,
for the acceleration pressure loss, g is acceleration of gravity; θ is the angle between pit shaft and stratum, and A is that wellbore section amasss, and D is mineshaft diameter, ρ
mfor m phase fluid density, f
tpfor friction pressure loss coefficient, V
mfor m phase unit volume.
6. numerical reservoir simulation method according to claim 5, is characterized in that, sets up expansion well model, comprising according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters:
According to flaw size, crack extension feature and divide crack grid with the way of contact of pit shaft, and subdivision is carried out to described crack grid;
Liquid flow equation discretization in the grid of described crack is utilized to form numerical model;
By fracture number value model be described expansion well model based on the well Model coupling of described Wellbore Flow equation.
7. numerical reservoir simulation method according to claim 6, is characterized in that, described expansion well model and the described geological model after initializing is coupled and solves coupling result, comprising:
The Nonlinear System of Equations with sparse coefficient matrix is set up based on described expansion well model and the described geological model after initializing;
Adopt Nonlinear System of Equations described in decoupling zero fully implicit solution Algorithm for Solving.
8. the numerical reservoir simulation method according to any one of claim 1-7, it is characterized in that, described dynamic static packet is drawn together: the well location coordinate of every mouthful of well, completion mode, perforating depth, well head, well depth, hole deviation and well track data, the structural configuration of oil reservoir target zone position, log analysis data, seismic interpretation data, layer data, the degree of porosity of target zone position, permeability, oil saturation and original formation pressure parameter, the high pressure property data of rock and fluid in stratum, characteristic of fluid data, wellbore tubular stream feature, individual well Production development data, well test analysis data and microseism decryption.
9. the reservoir numerical simulation device under the effect of horizontal fractures fracture guide, is characterized in that, described reservoir numerical simulation device comprises:
Fracture parameters determining unit, for obtaining the dynamic static data of well and oil reservoir, determines fracture parameters;
Liquid flow equation sets up unit, for setting up the liquid flow equation in horizontal fracture according to the flow conductivity in characteristic of fluid and described fracture parameters;
Wellbore Flow establishing equation unit, for setting up Wellbore Flow equation according to wellbore tubular stream feature;
Unit set up by expansion well model, for setting up expansion well model according to well track, completion position, described liquid flow equation, Wellbore Flow equation and fracture parameters;
Geological model generation unit, sets up geological model for the static parameter according to oil reservoir, and according to original formation pressure distribution or the relation of pressure and the degree of depth, and rock, fluid properties initialize described geological model;
Coupling unit, for being coupled described expansion well model and the described geological model after initializing and solving coupling result.
10. reservoir numerical simulation device according to claim 9, it is characterized in that, described fracture parameters determining unit specifically for: determine described fracture parameters according to the well test analysis data in described dynamic static data and microseism decryption, wherein, described fracture parameters comprises: fracture support width, fracture spacing, fracture half-length, fracture condudtiviy, the extension feature in crack and distribution layer position.
11. reservoir numerical simulation devices according to claim 10, is characterized in that, described liquid flow equation is set up unit and comprised:
Type of flow determination module, for according to the characteristic of fluid in described dynamic static data and described fracture condudtiviy, utilizes Reynolds number flow pattern judgment criterion to determine the crack inner fluid type of flow;
Liquid flow equation sets up module, for setting up the liquid flow equation in horizontal fracture according to the described crack inner fluid type of flow;
First computing module, splits pressure distribution in fracture and fluid flowing for calculating according to described liquid flow equation.
12. reservoir numerical simulation devices according to claim 11, is characterized in that, described first computing module specifically for:
Utilize binomial high speed non linear fluid flow through porous medium model to calculate and split pressure distribution in fracture and fluid flowing, ignore the flowing in z direction, described binomial high speed non linear fluid flow through porous medium model is as follows:
Wherein, x, y represent x direction and y direction respectively, and μ is fluid viscosity, k
ffor fracture permeabgility, ρ is fluid density, and v is seepage velocity, and ζ is high speed non linear fluid flow through porous medium coefficient.
13. reservoir numerical simulation devices according to claim 12, is characterized in that, described Wellbore Flow establishing equation unit comprises:
Establishing equation module, for setting up m phase Wellbore Flow equation according to wellbore tubular stream feature:
Second computing module, for distributing according to described Wellbore Flow equation calculating pressure;
Wherein, ρ
mgsin (θ) is hydrostatic power,
for friction pressure loss,
for the acceleration pressure loss, g is acceleration of gravity; θ is the angle between pit shaft and stratum, and A is that wellbore section amasss, and D is mineshaft diameter, ρ
mfor m phase fluid density, f
tpfor friction pressure loss coefficient, V
mfor m phase unit volume.
14. reservoir numerical simulation devices according to claim 13, is characterized in that, described expansion well model is set up unit and comprised:
Stress and strain model module, for according to flaw size, crack extension feature and divide crack grid with the way of contact of pit shaft, and carries out subdivision to described crack grid;
Descretization module, forms numerical model for utilizing liquid flow equation discretization in the grid of described crack;
Coupling module, for by fracture number value model be described expansion well model based on the well Model coupling of described Wellbore Flow equation.
15. reservoir numerical simulation devices according to claim 14, it is characterized in that, described coupling unit comprises:
Equation group sets up module, for setting up the Nonlinear System of Equations with sparse coefficient matrix based on described expansion well model and the described geological model after initializing;
Equation solution module, for adopting Nonlinear System of Equations described in decoupling zero fully implicit solution Algorithm for Solving.
16. reservoir numerical simulation devices according to any one of claim 9-15, it is characterized in that, described dynamic static packet is drawn together: the well location coordinate of every mouthful of well, completion mode, perforating depth, well head, well depth, hole deviation and well track data, the structural configuration of oil reservoir target zone position, log analysis data, seismic interpretation data, layer data, the degree of porosity of target zone position, permeability, oil saturation and original formation pressure parameter, the high pressure property data of rock and fluid in stratum, characteristic of fluid data, wellbore tubular stream feature, individual well Production development data, well test analysis data and microseism decryption.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510280351.7A CN105089595B (en) | 2015-05-27 | 2015-05-27 | Numerical reservoir simulation method and device under the effect of horizontal fractures fracture guide |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510280351.7A CN105089595B (en) | 2015-05-27 | 2015-05-27 | Numerical reservoir simulation method and device under the effect of horizontal fractures fracture guide |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105089595A true CN105089595A (en) | 2015-11-25 |
| CN105089595B CN105089595B (en) | 2019-04-12 |
Family
ID=54570828
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510280351.7A Active CN105089595B (en) | 2015-05-27 | 2015-05-27 | Numerical reservoir simulation method and device under the effect of horizontal fractures fracture guide |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105089595B (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106223939A (en) * | 2016-07-19 | 2016-12-14 | 中国石油天然气股份有限公司 | A kind of numerical reservoir determines method and apparatus |
| CN107844637A (en) * | 2017-10-18 | 2018-03-27 | 中国石油天然气股份有限公司 | A kind of numerical reservoir determines method and device |
| CN108019196A (en) * | 2017-11-22 | 2018-05-11 | 中国石油大学(北京) | A kind of volume fracturing horizontal well well test analysis processing method and processing device |
| CN108829902A (en) * | 2017-05-02 | 2018-11-16 | 中国石油化工股份有限公司 | A kind of appraisal procedure for tubing string hydraulic fracturing erosion damage |
| CN110080745A (en) * | 2019-05-16 | 2019-08-02 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | Separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and device |
| CN110630243A (en) * | 2018-06-21 | 2019-12-31 | 中国石油化工股份有限公司 | Method for determining fracturing fracture parameters based on fracturing well production data |
| CN110954419A (en) * | 2019-11-15 | 2020-04-03 | 山东大学 | Hydraulic fatigue fracture design method for preset fracture |
| CN112101710A (en) * | 2020-07-31 | 2020-12-18 | 中海油能源发展股份有限公司 | Quantitative injection-production balance adjustment method based on water drive front edge control |
| CN112282744A (en) * | 2020-10-23 | 2021-01-29 | 中国石油天然气股份有限公司 | Unconventional oil and gas reservoir well pattern deployment optimization method and device |
| CN112541287A (en) * | 2020-12-04 | 2021-03-23 | 中海石油(中国)有限公司 | Loose sandstone fracturing filling sand control production increase and profile control integrated design method |
| CN112580271A (en) * | 2020-12-09 | 2021-03-30 | 中国石油大学(华东) | Modeling method for integrating staged multi-cluster fracturing fracture expansion and proppant migration of horizontal well |
| CN112818503A (en) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | Method for determining moving speed of temporary plugging ball |
| CN113505472A (en) * | 2021-02-06 | 2021-10-15 | 中国石油大学(华东) | Numerical simulation method for repairing reservoir guanidine gum fracturing fluid damage by biological enzyme |
| CN112541287B (en) * | 2020-12-04 | 2024-04-23 | 中海石油(中国)有限公司 | Loose sandstone fracturing filling sand-prevention yield-increase profile control design method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5515923A (en) * | 1994-08-26 | 1996-05-14 | Loree; Dwight N. | Oil and gas well productivity |
| CN102203638A (en) * | 2008-09-19 | 2011-09-28 | 雪佛龙美国公司 | Computer-implemented systems and methods for use in modeling a geomechanical reservoir system |
| CN102418511A (en) * | 2011-07-25 | 2012-04-18 | 中国石油天然气股份有限公司 | Draw-down test analysis method for underground shut-in water well with low-permeability storage layer |
| CN104533370A (en) * | 2014-11-06 | 2015-04-22 | 中国石油大学(北京) | Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well |
| CN104594872A (en) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well |
-
2015
- 2015-05-27 CN CN201510280351.7A patent/CN105089595B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5515923A (en) * | 1994-08-26 | 1996-05-14 | Loree; Dwight N. | Oil and gas well productivity |
| CN102203638A (en) * | 2008-09-19 | 2011-09-28 | 雪佛龙美国公司 | Computer-implemented systems and methods for use in modeling a geomechanical reservoir system |
| CN102418511A (en) * | 2011-07-25 | 2012-04-18 | 中国石油天然气股份有限公司 | Draw-down test analysis method for underground shut-in water well with low-permeability storage layer |
| CN104533370A (en) * | 2014-11-06 | 2015-04-22 | 中国石油大学(北京) | Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well |
| CN104594872A (en) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106223939A (en) * | 2016-07-19 | 2016-12-14 | 中国石油天然气股份有限公司 | A kind of numerical reservoir determines method and apparatus |
| CN108829902A (en) * | 2017-05-02 | 2018-11-16 | 中国石油化工股份有限公司 | A kind of appraisal procedure for tubing string hydraulic fracturing erosion damage |
| CN108829902B (en) * | 2017-05-02 | 2021-06-25 | 中国石油化工股份有限公司 | Evaluation method for hydraulic fracturing erosion damage of tubular column |
| CN107844637A (en) * | 2017-10-18 | 2018-03-27 | 中国石油天然气股份有限公司 | A kind of numerical reservoir determines method and device |
| CN108019196A (en) * | 2017-11-22 | 2018-05-11 | 中国石油大学(北京) | A kind of volume fracturing horizontal well well test analysis processing method and processing device |
| CN108019196B (en) * | 2017-11-22 | 2019-11-05 | 中国石油大学(北京) | A kind of volume fracturing horizontal well well test analysis processing method and processing device |
| CN110630243A (en) * | 2018-06-21 | 2019-12-31 | 中国石油化工股份有限公司 | Method for determining fracturing fracture parameters based on fracturing well production data |
| CN110630243B (en) * | 2018-06-21 | 2023-04-28 | 中国石油化工股份有限公司 | Method for determining fracturing fracture parameters based on fracturing well production data |
| CN110080745B (en) * | 2019-05-16 | 2023-04-07 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | Method and device for predicting productivity of vertical well through separate-layer fracturing |
| CN110080745A (en) * | 2019-05-16 | 2019-08-02 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | Separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and device |
| CN110954419A (en) * | 2019-11-15 | 2020-04-03 | 山东大学 | Hydraulic fatigue fracture design method for preset fracture |
| CN110954419B (en) * | 2019-11-15 | 2020-11-10 | 山东大学 | Hydraulic fatigue fracture design method for preset fracture |
| CN112818503B (en) * | 2019-11-15 | 2022-11-01 | 中国石油天然气股份有限公司 | Method for determining moving speed of temporary plugging ball |
| CN112818503A (en) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | Method for determining moving speed of temporary plugging ball |
| CN112101710A (en) * | 2020-07-31 | 2020-12-18 | 中海油能源发展股份有限公司 | Quantitative injection-production balance adjustment method based on water drive front edge control |
| CN112282744A (en) * | 2020-10-23 | 2021-01-29 | 中国石油天然气股份有限公司 | Unconventional oil and gas reservoir well pattern deployment optimization method and device |
| CN112282744B (en) * | 2020-10-23 | 2023-11-28 | 中国石油天然气股份有限公司 | Unconventional oil and gas reservoir well pattern deployment optimization method and device |
| CN112541287A (en) * | 2020-12-04 | 2021-03-23 | 中海石油(中国)有限公司 | Loose sandstone fracturing filling sand control production increase and profile control integrated design method |
| CN112541287B (en) * | 2020-12-04 | 2024-04-23 | 中海石油(中国)有限公司 | Loose sandstone fracturing filling sand-prevention yield-increase profile control design method |
| CN112580271B (en) * | 2020-12-09 | 2022-11-04 | 中国石油大学(华东) | Modeling method for integrating staged multi-cluster fracturing fracture expansion and proppant migration of horizontal well |
| CN112580271A (en) * | 2020-12-09 | 2021-03-30 | 中国石油大学(华东) | Modeling method for integrating staged multi-cluster fracturing fracture expansion and proppant migration of horizontal well |
| CN113505472B (en) * | 2021-02-06 | 2022-09-02 | 中国石油大学(华东) | Numerical simulation method for repairing reservoir guanidine gum fracturing fluid damage by biological enzyme |
| CN113505472A (en) * | 2021-02-06 | 2021-10-15 | 中国石油大学(华东) | Numerical simulation method for repairing reservoir guanidine gum fracturing fluid damage by biological enzyme |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105089595B (en) | 2019-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105089595A (en) | Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion | |
| CN102022107B (en) | Method for establishing physical model capable of predicting waterflooding of fractured anisotropic oil reservoirs | |
| CN105808793B (en) | A kind of horizontal well fracturing method for numerical simulation based on unstrctured grid | |
| US9836561B2 (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics | |
| Fiallos Torres et al. | Modeling interwell fracture interference and Huff-N-Puff pressure containment in eagle ford using EDFM | |
| Bakhshi et al. | Numerical modeling and lattice method for characterizing hydraulic fracture propagation: a review of the numerical, experimental, and field studies | |
| CN104879103A (en) | Layered water injection effect analysis method | |
| CN104963657A (en) | Oil reservoir numerical simulation method and apparatus based on high water-cut period monitoring data constraints | |
| CN104360389B (en) | Tight sandstone reservoir rock elasticity modulus calculation method | |
| CN104453834A (en) | Injection-production relation optimizing and adjusting method for well group | |
| Pan et al. | TOUGH–RDCA modeling of multiple fracture interactions in caprock during CO2 injection into a deep brine aquifer | |
| US20190218895A1 (en) | Selective diffusion inclusion for a reservoir simulation for hydrocarbon recovery | |
| CN109577945A (en) | A kind of experimental provision and method of hypotonic-Oil in Super-low Permeability oil reservoir fluid-channeling channel differentiation | |
| Yu et al. | Embedded discrete fracture model assisted study of gas transport mechanisms and drainage area for fractured shale gas reservoirs | |
| Zhang et al. | Permeability enhancement and methane drainage capacity of tree-type boreholes to stimulate low-permeability coal seams | |
| US10846445B2 (en) | Method of waste distribution and formation pressure profile determination | |
| Cai et al. | Characteristics of hydraulic fracture penetration behavior in tight oil with multi-layer reservoirs | |
| Lei et al. | A dynamic discrete fracture model for fluid flow in fractured low-permeability reservoirs | |
| CN104392131A (en) | Calculating method of broken rock seepage field in process of driving sand by water | |
| Jin et al. | Investigation of produced gas injection in the Bakken for enhanced oil recovery considering well interference | |
| CN107766639A (en) | The computational methods of the natural gas lateral migration ultimate range of coefficient are reduced based on pressure | |
| CN104712299B (en) | It is adapted to gas well control water and increases the design method that air pressure splits | |
| Lisjak et al. | Development of a fully-coupled, hydro-mechanical model for finite-discrete element simulations of fluid-driven fracturing | |
| Jin et al. | Metric space mapping of oil sands reservoirs using streamline simulation | |
| Yao et al. | A numerical investigation on the hydraulic fracturing efficiency in radial well |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |