CN108979612A - A kind of densification oil-gas reservoir fracture acidizing complex fracture fluid ability optimization method - Google Patents
A kind of densification oil-gas reservoir fracture acidizing complex fracture fluid ability optimization method Download PDFInfo
- Publication number
- CN108979612A CN108979612A CN201810906243.XA CN201810906243A CN108979612A CN 108979612 A CN108979612 A CN 108979612A CN 201810906243 A CN201810906243 A CN 201810906243A CN 108979612 A CN108979612 A CN 108979612A
- Authority
- CN
- China
- Prior art keywords
- fracture
- reservoir
- crack
- flow
- fluid
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000005457 optimization Methods 0.000 title claims abstract description 25
- 238000000280 densification Methods 0.000 title description 2
- 230000000694 effects Effects 0.000 claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 230000008878 coupling Effects 0.000 claims abstract description 8
- 238000010168 coupling process Methods 0.000 claims abstract description 8
- 238000005859 coupling reaction Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 7
- 206010017076 Fracture Diseases 0.000 claims description 242
- 208000010392 Bone Fractures Diseases 0.000 claims description 240
- 239000007789 gas Substances 0.000 claims description 96
- 230000009466 transformation Effects 0.000 claims description 54
- 238000004519 manufacturing process Methods 0.000 claims description 48
- 230000035699 permeability Effects 0.000 claims description 43
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 36
- 239000003345 natural gas Substances 0.000 claims description 18
- 238000009826 distribution Methods 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 9
- 238000012913 prioritisation Methods 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 230000001052 transient effect Effects 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 210000003491 skin Anatomy 0.000 claims description 6
- 230000007480 spreading Effects 0.000 claims description 5
- 238000003892 spreading Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 claims description 4
- 230000001788 irregular Effects 0.000 claims description 4
- 230000000052 comparative effect Effects 0.000 claims description 3
- 210000002615 epidermis Anatomy 0.000 claims description 3
- 238000004613 tight binding model Methods 0.000 claims description 3
- 230000020477 pH reduction Effects 0.000 claims description 2
- 230000007812 deficiency Effects 0.000 abstract description 3
- 239000002253 acid Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 206010016256 fatigue Diseases 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005530 etching Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Business, Economics & Management (AREA)
- Primary Health Care (AREA)
- Human Resources & Organizations (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Marketing (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses the purpose of the present invention is to provide a kind of fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability optimization methods, it can be with overcome the deficiencies in the prior art using this method, the reservoir boundaries effect of fracture acidizing well is fully considered, stress sensitive, fluid non-Darcy flow influences in the heterogeneous variation of fracture condudtiviy and crack, the fluid flowing of fracture acidizing well is imported into secondary crack by matrix first respectively, major fracture, the fluid for importing secondary crack is flowed into the physical process in shaft bottom using major fracture, by coupling the flowing relation between them, to provide rational basis for the optimization of the flow in fracture ability of fracture acidizing well, improve reservoir reconstruction effect.
Description
Technical field
The present invention relates to oil and gas exploration and development fields, and in particular to multiple to a kind of fine and close oil-gas reservoir fracture acidizing
Miscellaneous flow in fracture ability optimization method.
Background technique
As oil and gas industry develops, the exploitation of tight sand and carbonate rock hydrocarbon reservoir is increasingly becoming emphasis.Preliminary assessment
The results show that the fine and close oily mining resources amount of China is 13 × 108~14 × 108T, tight gas mining resources amount are 8.8 × 1012~
12.1×1012m3, shale gas mining resources amount is 12.85 × 1012m3.Fracture acidizing transformation is fine and close reservoir exploration and development
Core indispensability technological means.Proppant is laid by hydraulic fracturing or acid fracturing etching generates complex fracture in reservoir and forms oil
The flow at high speed channel of gas is the key technology for improving reservoir reconstruction.During reservoir reconstruction, in addition to forming proppant support
Outside acid etching major fracture, since stress disturbance can also form a series of secondary crack around major fracture.In order to obtain pressure break acid
Change the maximum production of well, the man-made fracture fluid ability to be formed optimization is transformed in fracture acidizing should meet fluid in reservoir base simultaneously
Matter, secondary crack and major fracture are mutually matched.
At present for the optimization of fine and close oil-gas reservoir fracture acidizing fracture condudtiviy mainly consider matrix trandfer fluid with
Crack conducting power matches the corresponding research of development, yet there are no for reservoir base after fine and close oil-gas reservoir fracture acidizing transformation
Matter, secondary crack and complicated major fracture flowing are coupled, and are ground with obtaining the fracture condudtiviy optimization method of maximum oil and gas production
Study carefully.
In order to optimize fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability, the present invention considers fracture acidizing major fracture week
It encloses in the presence of secondary crack, fluid passes through the secondary crack of the non-homogeneous remittance of matrix first, further flowed into and led by secondary crack
Crack, the fluid for importing major fracture are flowed into shaft bottom further along major fracture；Optimizing fine and close oil-gas reservoir fracture acidizing complex fracture
During fluid ability, the non-homogeneous inflow time of complex fracture form, reservoir boundaries effect, stress sensitive, fluid has been fully considered
Grade crack and major fracture, the influence of the factors such as heterogeneous flow conductivity variation and high speed non-Darcy flow in crack, so that excellent
Change flow in fracture ability and is more in line with truth.
Summary of the invention
The object of the present invention is to provide a kind of fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability optimization methods.
In order to achieve the above object, providing a kind of flow in fracture ability optimization method of fracture acidizing well in the present invention, specifically
The following steps are included:
(S100) reservoir geology parameter, fluid basic parameter and wellbore parameters are collected；
(S200) fracture acidizing well secondary crack and major fracture parameter are collected；
(S300) fracture acidizing well major fracture is successively numbered, major fracture single-blade is divided equally by spatial spreading method
At the ns congruence of equal length, the Reservoir Seepage model of coupling reservoir boundaries effect, stress sensitive is established for each congruence；
(S400) consider that the secondary crack formed around fracture acidizing well major fracture improves reservoir permeability energy, work as fluid
Some pressure can be consumed less when flowing into major fracture by secondary crack, characterize fluid well by introducing a negative skin coefficient
The fluid ability in secondary crack establishes the Reservoir Seepage model of coupled secondary flow in fracture；
(S500) consider that fluid non-uniform flow along irregular cracks face behind secondary crack enters major fracture, consider major fracture
Pressure Drop Model of the fluid in hydraulic fracture is established in interior heterogeneous flow conductivity distribution, the influence of high speed non-Darcy effects；
(S600) the pressure break acid of coupled fluid high speed non-Darcy flow in reservoir matrix, secondary crack and major fracture is established
Change well production computation model；
(S700) under the conditions of given reservoir conditions, flowing bottomhole pressure (FBHP), to produce 1 year cumulative production comparative situation, optimization is caused
Close oil-gas reservoir fracture acidizing complex fracture fluid ability.
Further, reservoir geology parameter, fluid basic parameter, wellbore parameters are collected in the step (S100), wherein
The reservoir geology parameter specifically includes gas reservoir length, gas reservoir width, gas reservoir thickness, gas reservoir stress sensitive coefficients, gas reservoir constraint
Water saturation, reservoir temperature, reservoir permeability, reservoir porosity, original formation pressure；The fluid basic parameter specifically includes
The quasi- critical pressure of natural gas pseudo-critical temperature, natural gas, natural gas pseudoreduced temperature, gas compressive coefficient, natural gas are opposite
Density, natural gas density, Natural Gas Viscosity；The wellbore parameters specifically include flowing bottomhole pressure (FBHP) after wellbore radius, pressure.
Further, fracture acidizing well secondary crack is collected in the step (S200) and major fracture parameter specifically includes pressure
Split acidification well major fracture length, major fracture permeability, secondary crack transformation depth and secondary fracture permeabgility.
Further, the step (S300) specifically includes:
(S310) fracture acidizing well crack single-blade crack is divided into the ns congruence, constructs and closes for each congruence
The Green function expression of boundary box-shaped gas reservoir Point Source Function solution；
(S320) consider that there are stress sensitive effects for fine and close fracture acidizing well, establish fracture acidizing well reservoir matrix seepage flow mould
Type.
Further, the step (S400) comprising steps of
(S410) during fracture acidizing, since stress disturbance etc. causes the secondary crack of formation around major fracture to make reservoir
Effective permeability increases, and considers that secondary crack transformation depth perpendicular to fracture surface, is transformed depth using secondary crack and secondary is split
Seam area transformation permeability adds epidermis pressure drop to characterize to calculate as caused by secondary crack.
Further, the step (S500) comprising steps of
(S510) it establishes non-homogeneous water conservancy diversion crack and stitches interior high speed non-Darcy flow equation.
Further, the step (S600) comprising steps of
(S610) fracture acidizing well reservoir matrix-secondary crack-major fracture Coupled Flow model is established, fluid is seeped from reservoir
The physical process for flowing to wellbore is divided into Reservoir Seepage, three ranks of high speed non-Darcy flow in flowing and crack in secondary crack
Section, and gas from reservoir along the non-homogeneous inflow crack of fracture surface, according to pressure is continuous at the wall surface of crack and equal principle, establish
Pressure continuity equation；
(S620) for bottom pressure as definite value, the pressure at crack and wellbores is shaft bottom stream when setting flowing bottomhole pressure (FBHP) production
Pressure:
(S630) fracture acidizing well matrix-crack Coupled Flow Transient Flow model is established.
Further, the step (S700) comprising steps of
(S710) different fracture acidizing flow in fracture ability prioritization schemes are designed；
(S720) cumulative production for calculating different fracture acidizing flow in fracture ability prioritization scheme production a period of times, in conjunction with
Cumulative production result simultaneously considers engineering technology condition, selects best prioritization scheme, realizes that Optimum Fracturing is acidified complex fracture stream
Kinetic force.
In conclusion the invention has the following advantages that
The present invention provides a kind of fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability optimization methods, utilize this method
It can be with overcome the deficiencies in the prior art.By considering fine and close oil-gas reservoir fracture acidizing complex fracture form, and respectively, consideration master is split
Seam and secondary fracture pattern, while having also contemplated the secondary crack of the non-homogeneous inflow of reservoir boundaries effect, stress sensitive, fluid and master
Crack, the influence of the factors such as heterogeneous flow conductivity variation and high speed non-Darcy flow in crack, so that the model established is examined
Worry factor is more perfect, to be more in line with actual conditions, optimizes fine and close oil-gas reservoir fracture acidizing complex fracture flowing to reach
Ability.
Detailed description of the invention
Fig. 1 is fracture acidizing well complex fracture form schematic diagram；
Fig. 2 is that area's secondary permeability is transformed to the influence diagram of tired gas production；
Fig. 3 is influence diagram of the transformation with permeability size to discrete unit flow distribution；
Fig. 4 is influence diagram of the transformation with depth to tired gas production；
Fig. 5 is influence diagram of the transformation with depth to discrete unit flow distribution；
Specific embodiment
The purpose of the present invention is to provide a kind of fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability optimization method, benefits
The reservoir boundaries effect of fracture acidizing well with overcome the deficiencies in the prior art, can have been fully considered with this method, stress sensitive, split
Stitching fluid non-Darcy flow in the heterogeneous variation of flow conductivity and crack influences, and the fluid flowing of fracture acidizing well is logical first
It crosses matrix and imports branch fractures, major fracture respectively, the fluid for importing branch fractures is flowed into the physics in shaft bottom using major fracture
Process, by coupling the flowing relation between them, thus for fracture acidizing well flow in fracture ability optimization provide rationally according to
According to raising reservoir reconstruction effect.
A kind of flow in fracture ability optimization method of fracture acidizing well, mainly comprises the steps that
(S100) reservoir geology parameter, fluid basic parameter and wellbore parameters are collected；
(S200) fracture acidizing well secondary crack and major fracture parameter are collected；
(S300) fracture acidizing well major fracture is successively numbered, major fracture single-blade is divided equally by spatial spreading method
At the ns congruence of equal length, the Reservoir Seepage model of coupling reservoir boundaries effect, stress sensitive is established for each congruence；
(S400) consider that the secondary crack formed around fracture acidizing well major fracture improves reservoir permeability energy, work as fluid
Some pressure can be consumed less when flowing into major fracture by secondary crack, characterize fluid well by introducing a negative skin coefficient
The fluid ability in secondary crack establishes the Reservoir Seepage model of coupled secondary flow in fracture；
(S500) consider that fluid non-uniform flow along irregular cracks face behind secondary crack enters major fracture, consider major fracture
Pressure Drop Model of the fluid in hydraulic fracture is established in interior heterogeneous flow conductivity distribution, the influence of high speed non-Darcy effects；
(S600) the pressure break acid of coupled fluid high speed non-Darcy flow in reservoir matrix, secondary crack and major fracture is established
Change well production computation model；
(S700) under the conditions of given reservoir conditions, flowing bottomhole pressure (FBHP), to produce 1 year cumulative production comparative situation, optimization is caused
Close oil-gas reservoir fracture acidizing complex fracture fluid ability.
In the present invention, reservoir geology parameter, fluid basic parameter, wellbore parameters are collected in the step (S100),
Described in reservoir geology parameter specifically include gas reservoir length, gas reservoir width, gas reservoir thickness, gas reservoir stress sensitive coefficients, gas reservoir beam
Tie up water saturation, reservoir temperature, reservoir permeability, reservoir porosity, original formation pressure；The fluid basic parameter specifically wraps
Include natural gas pseudo-critical temperature, natural gas intends critical pressure, natural gas pseudoreduced temperature, gas compressive coefficient, natural gas phase
To density, natural gas density, Natural Gas Viscosity；The wellbore parameters specifically include flowing bottomhole pressure (FBHP) after wellbore radius, pressure.
In the present invention, fracture acidizing well secondary crack is collected in the step (S200) and major fracture parameter specifically includes
Fracture acidizing well major fracture length, major fracture permeability, secondary crack transformation depth and secondary fracture permeabgility.
In the present invention, fracture acidizing well major fracture is successively numbered the step (S300), by major fracture single-blade
It is divided into the ns congruence of equal length by spatial spreading method, establishes coupling reservoir boundaries effect for each congruence, answer
The Reservoir Seepage model of power sensitivity, specifically includes:
(S310) fracture acidizing well crack single-blade crack is divided into the ns congruence, constructs and closes for each congruence
The Green function expression of boundary box-shaped gas reservoir Point Source Function solution:
Wherein:
In formula:Indicate coordinate points (x in infinity plane0,y0) withDetermine mass flow production (t-t0) time
Afterwards in the instantaneous pseudopressure of coordinate points (x, y), MPa2/(Pa·s)；Indicate prime stratum pseudopressure, MPa2/(Pa·s)；Indicates coordinate point (x0,y0) constant flow production when yield, kg/ks；φ indicate reservoir matrix porosity, it is no because
It is secondary；CtIndicate fluid compressibility, MPa-1；N indicates digit, dimensionless；τ indicates quantity-produced duration, ks；χ
Indicate piezometric conductivity, m2MPa/ (Pas), χ=K/ (μ ctφ), K indicates reservoir original permeability, m2；μ indicates that fluid is viscous
Degree, Pas；T is indicated from the production time for starting to measure when production, ks；xeIndicate closed boundary box-shaped gas reservoir region in the direction x
On two boundaries be located at x=0 and x=xe；xwIndicate the coordinate of the congruence in the x direction, m；yeIndicate closed boundary box-shaped
Two boundaries of gas reservoir region in y-direction are located at y=0 and y=ye；ywIndicate the coordinate of the congruence in y-direction, m.
According to actual gas state equation, yield under the mark condition of ground can be calculated:
In formula:
ρscGas density under-mark condition, kg/m3；
The current strata pressure of p-, MPa；
Convolution (1) and (3), can be obtained the Point Source Function of closed boundary box-shaped gas reservoir:
In formula:
P --- strata pressure at present, MPa；
Mark the volume flow under condition, m in q --- ground3/ks；
psc--- mark condition pressure, MPa；
ρsc--- gas density under mark condition, kg/m3；
Tsc--- mark condition temperature, K；
T --- reservoir temperature, K；
Z --- the gas deviation factor under current reservoir pressure, dimensionless；
Zsc--- the gas deviation factor under the status of criterion, dimensionless.
In fracture acidizing well production process, since reservoir pressure decline leads to permeability reduction, any time storage
Layer permeability is the function of strata pressure, i.e. KpIndicate the reservoir permeability under stress sensitive effect.
Kp=Kexp [- α (pi-p)] (5)
In formula: K --- reservoir original permeability, mD；
Kp--- the current permeability of reservoir, mD；
α --- reservoir stress sensitive coefficients, MPa-1；
pi--- original formation pressure, MPa；
P --- strata pressure at present, MPa；
The present invention consider tight gas reservoir be closing box-shaped gas reservoir, at present strata pressure p by closing box-shaped constant volume gas reservoir object
Matter equilibrium equation formula (6) is calculated:
Deviation factor for gas under the current strata pressure of Z-, zero dimension；
ZiDeviation factor for gas under-original formation pressure, zero dimension；
Gp- fracture acidizing gas well cumulative production, m3；
G-original oil in place, m3；G=xe·ye·h·(1-sw)/Bg；
H --- closing box-shaped gas reservoir height, m；
sw- water saturation, %.
(S320) consider that there are stress sensitive effects for fine and close fracture acidizing well, establish fracture acidizing well reservoir matrix seepage flow mould
Type；
Using spatial spreading crack method, can by major fracture single-blade it is discrete be ns congruence source, for each congruence source
Pressure response when pressure response in process of production can be produced by each congruence source is superimposed to obtain.Any bit on crack
Set M (xf,i,yf,i) at by yield be qf,jCongruence source N (xf,j,yf,j) generate pressure response are as follows:
In formula: pf,i--- pressure, MPa at jth infinitesimal section center on crack；
qf,j--- the i-th infinitesimal section marks the volume flow under condition, m on ground on crack3/ks；
(xf,i, yf,i) --- the coordinate of i-th of discrete unit, m on crack；
(xf,j, yf,j) --- the coordinate of j-th of discrete unit, m on crack；
I, j --- crack discrete unit number, zero dimension.
It is ns infinitesimal section since fracture acidizing single-blade crack is discrete, using above-mentioned thought, then it is total that whole crack can be obtained
2ns discrete unit is when t moment is produced in stratum point M (xf,i,yf,i) at generate pressure response equation:
In formula: ns --- fracture acidizing single-blade crack discrete unit number, it is a；
In the present invention, the step (S400) considers that the secondary crack formed around fracture acidizing well major fracture improves
Reservoir permeability energy characterizes the fluid ability in fluid well secondary crack by introducing a negative skin coefficient, establishes coupled secondary
The Reservoir Seepage model of flow in fracture, specifically comprises the following steps；
(S410) during fracture acidizing, since stress disturbance etc. causes the secondary crack of formation around major fracture to make reservoir
Effective permeability increases, and considers that secondary crack transformation depth perpendicular to fracture surface, is transformed depth using secondary crack and secondary is split
Seam area transformation permeability adds epidermis pressure drop to characterize to calculate as caused by secondary crack；
In formula: sf--- skin factor, zero dimension is transformed in fracture surface secondary crack area；
D --- regional depth, m is transformed in secondary crack；
xf--- man-made fracture length, m；
Kc--- region permeability, mD is transformed in secondary crack.
According to the different actual conditions of secondary crack transformation depth at fracture acidizing major fracture different location, consider along crack
Depth is transformed on length direction to be linearly gradually reduced, in crack, heel transformation depth is up to dfmax, fracture tip transformation is deeply
Spend minimum dfmin, then the transformation depth of i-th of fracture acidizing crack discrete unit may be expressed as: on crack
In formula: di--- the transformation depth of discrete unit, m on i-th of fracture acidizing crack；
dfmin--- the transformation depth of fracture tip discrete unit, m on fracture acidizing crack；
dfmax--- the transformation depth of crack heel discrete unit, m on fracture acidizing crack；
xf,i--- distance of i-th of fracture acidizing crack discrete unit to crack toe-end, m.
Formula (10) are substituted into formula (9), and consider that the length of i-th of discrete unit is Δ xf,i, obtain i-th discrete list
The transformation band skin factor s of memberf,iAre as follows:
Compare the pseudopressure difference variation of secondary crack the reconstruction scope intrinsic permeability variation front and back, it can be deduced that i-th crack
The parasitic pressure drop that discrete unit is generated by secondary transformation crackAre as follows:
Formula (8), formula (12) are combined to the reservoir matrix percolation equationk for obtaining considering secondary flow in fracture ability, formula (13):
In the present invention, the step (S500) considers fluid behind secondary crack along irregular cracks face non-uniform flow
Enter major fracture, considers that heterogeneous flow conductivity distribution, the influence of high speed non-Darcy effects establish fluid in hydraulic fracture in major fracture
Interior Pressure Drop Model, specifically includes step:
(S510) it is based on Forchheimer equation, non-homogeneous water conservancy diversion crack is established and stitches interior high speed non-Darcy flow equation are as follows:
In formula: pf,i--- crack fluid pressure, Pa in the middle part of i-th of discrete unit；
vf,i--- crack fluid velocity inside, m/s in the middle part of i-th of discrete unit；
βg,i--- i-th of discrete unit crack fluid velocity inside coefficient, m-1；
ρg,i--- fluid density in i-th of discrete unit crack, kg/m3；
Kf,i--- i-th of discrete unit fracture permeabgility, m2；
Wherein βi,jNumerical values recited can be calculated with formula (15), consider acid fracturing during major fracture permeability
It is constant, therefore velocity coeffficient can be further written as:
βg,i=β=7.644 × 1010/Kf,i 1.5=7.644 × 1010/K1.5 (15)
In formula (15), gross pressure gradient delta pf,i/Δxf,iIt is made of two parts, first part is equation right end first item
Seam in Darcy Flow pressure drop, second part be equation right end Section 2 seam in high speed non-Darcy effects generate flowing pressure
Drop, by the Section 2 in formula (14), that is, non-Darcy flow pressure drop symbol pDfi,jIt indicates, then has:
In formula
γ in formulag--- gas relative density, dimensionless；
Mair--- air molecule quality, g/mol；
R --- gas constant, dimensionless；
Z --- Gas Compression Factor, dimensionless；
T --- reservoir temperature, K；
wf,i--- the slit width on crack at the i-th discrete unit, m；
Bg--- gas volume factor, dimensionless；
qsc--- gas yield under standard state, m3/d；
psc--- gas reservoir pressure under standard state, m3/d；
Tsc--- reservoir temperature under standard state, K.
Actual conditions, the discrete side in application space are become narrow gradually by heel to toe according to fracture acidizing crack single-blade slit width
Method handles each crack infinitesimal for isosceles trapezoid, i.e. every man-made fracture single-blade is made of ns isosceles trapezoid, thus real
The wedge shape variation long along seam of existing slit width.The slit width w of i-th of crack discrete unitf,iIt may be expressed as:
W in formulai,j--- the i-th crack numbers the middle part width of upper j-th of discrete unit, mm；
wmin,i--- the i-th crack toe-end width, mm；
wmax--- the i-th crack heel end width, mm.
Gas can generate non-linear after injuring the non-homogeneous inflow crack of fracture surface because of fluid high-speed non-Darcy effects
It flows, then any position O (x on crackf,i, yf,i) arrive wellbore O (xf,0, yf,0) generate total droop lossAnd it converts
Yield under marking condition to ground are as follows:
In the present invention, the step (S600) establishes coupling stream high speed in reservoir matrix, secondary crack and major fracture
The fracture acidizing well production computation model of non-Darcy flow, specifically includes:
(S610) fracture acidizing well reservoir matrix-secondary crack-major fracture Coupled Flow model is established, fluid is seeped from reservoir
The physical process for flowing to wellbore is divided into Reservoir Seepage, three ranks of high speed non-Darcy flow in flowing and crack in secondary crack
Section, and gas from reservoir along the non-homogeneous inflow crack of fracture surface, according to pressure is continuous at the wall surface of crack and equal principle
Point O (x according to the observationf,j,yf,j) at pressure it is continuous, and pressure continuity equation is established by formula (13) and (21)；
(S620) bottom pressure is definite value when setting stable bottom hole pressure production, O (x at crack and wellboresf,0, yf,0)
Pressure be flowing bottomhole pressure (FBHP):
pf,0=pwf (22)
P in formula0--- the pressure at man-made fracture and wellbores, MPa；
pwf--- flowing bottomhole pressure (FBHP), MPa.
(S630) fracture acidizing well matrix-crack Coupled Flow Transient Flow model is established；
Matrix-crack coupling of the discrete infinitesimal section in the i-th crack of t moment can be obtained in simultaneous equations (13), (21) and (22)
The Transient Flow continuity equation of flowing:
It for formula (23), is made of 2ns equation altogether, wherein each discrete segments flow is unknown number, that is, has 2ns
Unknown number, due to equation number with unknown number be it is equal, mathematical model can solve.It can solve to obtain every under transient state
The flow of a discrete unit, so that superposition obtains fracture acidizing well yield:
Q --- the daily output of fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability, m3/d。
Formula (24) just establishes the transient state Production rate method of fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability.
In order to calculate the transient state yield formula under any production time t=Δ t.According to time principle of stacking, it can solve and take the post as
The unstable state deliverability equation anticipated under moment t=n Δ t.The yield that daily crack is calculated is overlapped, so that it may be tired out
Count yield.
In formula:
Gp--- cumulative production, m3；
T --- production time calculator, day；
D --- accumulative production time, day；
In the present invention, the step (S700) is accumulative with production 1 year under the conditions of given reservoir conditions, flowing bottomhole pressure (FBHP)
Yield comparison situation optimizes fine and close oil-gas reservoir fracture acidizing complex fracture fluid ability, wherein specifically including:
(S710) different fracture acidizing flow in fracture ability prioritization schemes are designed；
(S720) cumulative production that different fracture acidizing flow in fracture ability prioritization schemes produce 1 year is calculated, in conjunction with accumulative
Yield result simultaneously considers engineering technology condition, selects best prioritization scheme, realizes that Optimum Fracturing is acidified the complex fracture energy of flow
Power.
It is of the invention in order to facilitate understanding by those skilled in the art, it has carried out calculating analysis using actual optimization case.
Table 4-1 gas reservoir and fracture acidizing crack basic parameter table
(1) transformation band permeability
Consider that tight gas reservoir fracture acidizing forms complex fracture, is to form a series of transformation areas secondary around major fracture
Crack.Shown in its basic physical model Fig. 1, area's permeability is transformed in analysis first, i.e., secondary fracture permeabgility is to fracture acidizing
Well tires out the influence of gas production.
(a) influence of the transformation area's secondary fracture permeabgility to tired gas production
The permeability that transformation area is arranged is set as tetra- kinds of situations of 1mD, 3mD, 5mD and 7mD, analyzes transformation area's secondary crack
Permeability tires out the influence of gas production to fractured well.
The influence of gas production is tired out to fine and close oil-gas reservoir fracture acidizing complex fracture well from Fig. 2 difference transformation area's fracture permeabgility
It can be seen that under the conditions of other parameters identical parameters, with the increase of transformation area's fracture permeabgility, the tired gas production of fractured well
It is higher.This is because transformation area's permeability is bigger, gas filtrational resistance when flowing through transformation area is smaller, so as to cause tired gas production
It is higher.Influence from different transformation area's fracture permeabgilitys to tired gas production can be seen that when transformation area's permeability increases to centainly
It after value, further increases transformation area's permeability and becomes smaller on the influence of the amplification of tired gas production, preferred transformation area's permeability is accordingly
5mD。
(b) influence of the transformation area's permeability to discrete segments flow distribution
It is discrete to fine and close oil-gas reservoir fracture acidizing complex fracture well when producing 1d from Fig. 3 difference transformation area's fracture permeabgility
The influence of unit of flow distribution is it can be seen that under the conditions of other parameters identical parameters, and transformation area's permeability is bigger, then crack is each
The yield of discrete unit is higher；The yield distribution of each discrete unit is showed close to yield highest at pit shaft,
Far from flow at pit shaft in first reducing the trend increased afterwards, and between each crack discrete unit yield difference in size with transformation
The increase of area's permeability and increase, i.e., each crack discrete unit yield difference is smaller under lower transformation area permeability；
Fine and close oil-gas reservoir fracture acidizing complex fracture well yield can be significantly improved with fracture permeabgility by increasing to be transformed.
(2) influence of the reconstruction scope
(a) influence of the transformation depth to tired gas production
The influence of gas production, setting are tired out further to study transformation depth to fine and close oil-gas reservoir fracture acidizing complex fracture well
Be transformed area permeability be 5mD, the reconstruction scope that area is transformed be three kinds of situations: comparison setting be transformed the penetration of fracture be 10m, 20m,
The influence of 30m, 40m to tired gas production.
The influence for tiring out gas production to fine and close oil-gas reservoir fracture acidizing complex fracture well from Fig. 4 difference transformation area's depth can
Out: under the conditions of other parameters identical parameters, with the increase of transformation depth (range), the tired gas production of fractured well is higher.It is right
Than transformation area's depth and area's permeability is transformed to the tired influence for producing gas as can be seen that influence of the transformation area's permeability to tired gas production
It is significantly higher than influence of the transformation area's depth to tired gas production, therefore pays the utmost attention to improve transformation area's permeability in the construction process
To further increase yield.
(b) influence of the transformation area's depth to discrete segments flow distribution
From Fig. 5 difference transformation area's depth when producing 1d to fine and close oil-gas reservoir fracture acidizing complex fracture well discrete unit stream
The influence of distribution is measured it can be seen that transformation area's depth is bigger, then each discrete unit in crack under the conditions of other parameters identical parameters
Yield it is higher；The yield distribution of each discrete unit is showed close to yield highest at pit shaft, far from pit shaft
Place's flow in first reducing the trend increased afterwards, and between each crack discrete unit yield difference in size with being transformed with depth
Increase and increase, i.e., each crack discrete unit yield difference is smaller under smaller transformation band depth.
Although being described in detail in conjunction with attached drawing to a specific embodiment of the invention, should not be construed as special to this
The restriction of the protection scope of benefit.In range described by claims, those skilled in the art are without creative work
The various modifications and deformation that can make still belong to the protection scope of this patent.
Claims (8)
1. a kind of flow in fracture ability optimization method of fracture acidizing well, comprising the following steps:
(S100) reservoir geology parameter, fluid basic parameter and wellbore parameters are collected；
(S200) fracture acidizing well secondary crack and major fracture parameter are collected；
(S300) fracture acidizing well major fracture is successively numbered, major fracture single-blade is divided equally by spatial spreading method and is grown up
The equal ns congruence is spent, the Reservoir Seepage model of coupling reservoir boundaries effect, stress sensitive is established for each congruence；
(S400) consider that the secondary crack formed around fracture acidizing well major fracture improves reservoir permeability energy, when fluid passes through
Secondary crack can consume some pressure less when flowing into major fracture, characterize fluid well secondary by introducing a negative skin coefficient
The fluid ability in crack establishes the Reservoir Seepage model of coupled secondary flow in fracture；
(S500) consider that fluid non-uniform flow along irregular cracks face behind secondary crack enters major fracture, consider non-in major fracture
Pressure Drop Model of the fluid in hydraulic fracture is established in the distribution of homogeneous flow conductivity, the influence of high speed non-Darcy effects；
(S600) the fracture acidizing well of coupled fluid high speed non-Darcy flow in reservoir matrix, secondary crack and major fracture is established
Production rate model；
(S700) under the conditions of given reservoir conditions, flowing bottomhole pressure (FBHP), to produce 1 year cumulative production comparative situation, optimize fine and close oil
Gas reservoir fracture acidizing complex fracture fluid ability.
2. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, in the step (S100)
Collect reservoir geology parameter, fluid basic parameter, wellbore parameters, wherein the reservoir geology parameter specifically include gas reservoir length,
Gas reservoir width, gas reservoir thickness, gas reservoir stress sensitive coefficients, gas reservoir irreducible water saturation, reservoir temperature, reservoir permeability, reservoir
Porosity, original formation pressure；The fluid basic parameter specifically includes natural gas pseudo-critical temperature, natural gas intends critical pressure
Power, natural gas pseudoreduced temperature, gas compressive coefficient, natural gas relative density, natural gas density, Natural Gas Viscosity；It is described
Wellbore parameters specifically include flowing bottomhole pressure (FBHP) after wellbore radius, pressure.
3. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, in the step (S200)
It collects fracture acidizing well secondary crack and major fracture parameter specifically includes fracture acidizing well major fracture length, major fracture permeability,
Secondary crack transformation depth and secondary fracture permeabgility.
4. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, step (S300) tool
Body includes:
(S310) fracture acidizing well crack single-blade crack is divided into the ns congruence, constructs closed boundary for each congruence
The Green function expression of box-shaped gas reservoir Point Source Function solution；
(S320) consider that there are stress sensitive effects for fine and close fracture acidizing well, establish fracture acidizing well reservoir matrix flow model in porous media.
5. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, step (S400) packet
Include step:
(S410) during fracture acidizing, since stress disturbance etc. causes the secondary crack of formation around major fracture to make reservoir effective
Permeability increases, and considers secondary crack transformation depth perpendicular to fracture surface, using secondary crack transformation depth and secondary crack area
Transformation permeability adds epidermis pressure drop to characterize to calculate as caused by secondary crack.
6. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, step (S500) packet
Include step:
(S510) it establishes non-homogeneous water conservancy diversion crack and stitches interior high speed non-Darcy flow equation.
7. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, step (S600) packet
Include step:
(S610) establish fracture acidizing well reservoir matrix-secondary crack-major fracture Coupled Flow model, fluid from Reservoir Seepage to
The physical process of wellbore is divided into Reservoir Seepage, high speed non-Darcy flow three phases in flowing and crack in secondary crack, and
Gas along the non-homogeneous inflow crack of fracture surface, according to pressure is continuous at the wall surface of crack and equal principle, establishes pressure from reservoir
Continuity equation；
(S620) for bottom pressure as definite value, the pressure at crack and wellbores is flowing bottomhole pressure (FBHP) when setting flowing bottomhole pressure (FBHP) production:
(S630) fracture acidizing well matrix-crack Coupled Flow Transient Flow model is established.
8. a kind of flow in fracture ability optimization method of fracture acidizing well as described in claim 1, step (S700) packet
Include step:
(S710) different fracture acidizing flow in fracture ability prioritization schemes are designed；
(S720) cumulative production for calculating different fracture acidizing flow in fracture ability prioritization scheme production a period of times, in conjunction with accumulative
Yield result simultaneously considers engineering technology condition etc., selects best prioritization scheme, realizes Optimum Fracturing acidification complex fracture flowing
Ability.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810906243.XA CN108979612B (en) | 2018-08-09 | 2018-08-09 | Method for optimizing flow capacity of complex fracture of fracturing and acidizing of compact oil and gas reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810906243.XA CN108979612B (en) | 2018-08-09 | 2018-08-09 | Method for optimizing flow capacity of complex fracture of fracturing and acidizing of compact oil and gas reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108979612A true CN108979612A (en) | 2018-12-11 |
| CN108979612B CN108979612B (en) | 2020-07-07 |
Family
ID=64555627
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810906243.XA Active CN108979612B (en) | 2018-08-09 | 2018-08-09 | Method for optimizing flow capacity of complex fracture of fracturing and acidizing of compact oil and gas reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108979612B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109653740A (en) * | 2018-12-26 | 2019-04-19 | 西南石油大学 | A kind of Carbonate Acid Fracturing self-supporting fracture permeabgility prediction technique |
| CN109882164A (en) * | 2019-03-28 | 2019-06-14 | 中国石油大学(华东) | A kind of large scale acidification analogy method of fractured carbonate rock oil reservoir |
| CN110397430A (en) * | 2019-07-08 | 2019-11-01 | 西南石油大学 | A method of prediction carbonate rock fore negative bed fracture condudtiviy distribution |
| CN110410054A (en) * | 2019-06-25 | 2019-11-05 | 中联煤层气有限责任公司 | A kind of coal bed gas well fracturing fracture becomes the prediction technique of space-time flow conductivity |
| CN110454152A (en) * | 2019-08-12 | 2019-11-15 | 成都北方石油勘探开发技术有限公司 | The non-homogeneous injury skin factor calculation method of fracture-type reservoir |
| CN112878982A (en) * | 2020-12-31 | 2021-06-01 | 西南石油大学 | Deep shale gas productivity prediction method considering long-term fracture conductivity |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104594872A (en) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well |
| RU2578134C1 (en) * | 2015-03-11 | 2016-03-20 | Открытое акционерное общество "Татнефть" имени В.Д. Шашина | Method of developing oil deposits in fractured reservoirs with water oil zones |
| CN105422071A (en) * | 2015-12-07 | 2016-03-23 | 西南石油大学 | Method for evaluating rationality of low-permeability non-homogeneous gas reservoir fracturing horizontal well fracture parameters |
| CN107220493A (en) * | 2017-05-24 | 2017-09-29 | 王欣 | Shale gas horizontal well network fracture modeling method based on micro-seismic event |
| CN107480411A (en) * | 2017-10-20 | 2017-12-15 | 合肥工业大学 | Reservoir fracturing effect evaluation method and evaluation system |
-
2018
- 2018-08-09 CN CN201810906243.XA patent/CN108979612B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104594872A (en) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well |
| RU2578134C1 (en) * | 2015-03-11 | 2016-03-20 | Открытое акционерное общество "Татнефть" имени В.Д. Шашина | Method of developing oil deposits in fractured reservoirs with water oil zones |
| CN105422071A (en) * | 2015-12-07 | 2016-03-23 | 西南石油大学 | Method for evaluating rationality of low-permeability non-homogeneous gas reservoir fracturing horizontal well fracture parameters |
| CN107220493A (en) * | 2017-05-24 | 2017-09-29 | 王欣 | Shale gas horizontal well network fracture modeling method based on micro-seismic event |
| CN107480411A (en) * | 2017-10-20 | 2017-12-15 | 合肥工业大学 | Reservoir fracturing effect evaluation method and evaluation system |
Non-Patent Citations (1)
| Title |
|---|
| 柯玉彪: ""多重运移机制页岩气藏压裂水平井产能研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109653740A (en) * | 2018-12-26 | 2019-04-19 | 西南石油大学 | A kind of Carbonate Acid Fracturing self-supporting fracture permeabgility prediction technique |
| CN109882164A (en) * | 2019-03-28 | 2019-06-14 | 中国石油大学(华东) | A kind of large scale acidification analogy method of fractured carbonate rock oil reservoir |
| CN110410054A (en) * | 2019-06-25 | 2019-11-05 | 中联煤层气有限责任公司 | A kind of coal bed gas well fracturing fracture becomes the prediction technique of space-time flow conductivity |
| CN110410054B (en) * | 2019-06-25 | 2022-04-05 | 中联煤层气有限责任公司 | Prediction method for time-space flow conductivity of fracturing fracture of coal-bed gas well |
| CN110397430A (en) * | 2019-07-08 | 2019-11-01 | 西南石油大学 | A method of prediction carbonate rock fore negative bed fracture condudtiviy distribution |
| CN110454152A (en) * | 2019-08-12 | 2019-11-15 | 成都北方石油勘探开发技术有限公司 | The non-homogeneous injury skin factor calculation method of fracture-type reservoir |
| CN110454152B (en) * | 2019-08-12 | 2020-12-01 | 成都北方石油勘探开发技术有限公司 | Method for calculating heterogeneous damage skin coefficient of fractured reservoir |
| CN112878982A (en) * | 2020-12-31 | 2021-06-01 | 西南石油大学 | Deep shale gas productivity prediction method considering long-term fracture conductivity |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108979612B (en) | 2020-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108979612A (en) | A kind of densification oil-gas reservoir fracture acidizing complex fracture fluid ability optimization method | |
| CN102606129B (en) | Method and system for thin interbed oilfield development | |
| Qun et al. | The “fracture-controlled reserves” based stimulation technology for unconventional oil and gas reservoirs | |
| CN109025942A (en) | A kind of irregular cracky Production rate method of tight gas reservoir Deviated Well Fracturing | |
| CN110334431A (en) | A kind of low permeability tight gas reservoir single well controlled reserves calculating and remaining gas analysis method | |
| CN111236908A (en) | Multi-stage fractured horizontal well productivity prediction model and productivity sensitivity analysis method suitable for low-permeability tight gas reservoir | |
| CN109033677A (en) | A kind of fracture acidizing well fracture condudtiviy optimization method | |
| CN109033674B (en) | Method for optimizing fracture parameters of target fracturing acidizing well | |
| CN104594872A (en) | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well | |
| CN104453834A (en) | Injection-production relation optimizing and adjusting method for well group | |
| CN103577886A (en) | Staged fracturing yield prediction method of low-permeability gas reservoir horizontal well | |
| CN103590824A (en) | Productivity calculation method for tight gas horizontal wells modified by multi-stage fracturing | |
| CN106437674A (en) | Imitation horizontal well waterflooding well pattern adapting method | |
| CN105089612A (en) | Determining method for distance of well-drain and length of pressure break of low penetration oil reservoir artificial fracture | |
| CN104879104B (en) | Oil reservoir water injection method | |
| CN103400020A (en) | Numerical reservoir simulation method for measuring flowing situation of a plurality of intersected discrete cracks | |
| CN107630686A (en) | Interlock between net horizontal section the fine and close oily energy supplement method of displacement and the imbibition displacement of reservoir oil | |
| CN107066769A (en) | Suitable for the efficient acidifying design method of ultra deep slit formation carbonate reservoir | |
| CN113177322B (en) | Fracturing single well control reserve calculation method | |
| CN108614902A (en) | A kind of shale gas well creation data analysis method producing system variation | |
| CN103399970A (en) | Method for oil deposit flow state digifax calculating with discrete fracture line processing | |
| CN106321051A (en) | Method for optimizing multi-section fractured horizontal well network crack parameter | |
| CN206309375U (en) | The improvement well pattern structure of rhombic inverted nini-spot well pattern | |
| CN107503739A (en) | A kind of pressure monitoring method that inflow direction is differentiated for horizontal well | |
| CN107832540A (en) | A kind of compact oil reservoir technical limit well space determines method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |