CN108829945A - One kind assessing storey increase design volume method based on Fracturing Pressure in real time - Google Patents
One kind assessing storey increase design volume method based on Fracturing Pressure in real time Download PDFInfo
- Publication number
- CN108829945A CN108829945A CN201810530032.0A CN201810530032A CN108829945A CN 108829945 A CN108829945 A CN 108829945A CN 201810530032 A CN201810530032 A CN 201810530032A CN 108829945 A CN108829945 A CN 108829945A
- Authority
- CN
- China
- Prior art keywords
- fracturing
- pressure
- mpa
- frictional resistance
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The invention belongs to oilfield prospecting developing hydraulic fracturing technology fields, are related to one kind based on Fracturing Pressure and assess storey increase design volume method in real time, include the following steps:S1, target well reservoir geology and engineering parameter, pressing crack construction basic parameter are collected;S2, fracturing string flowing frictional resistance accounting equation is established, flowing frictional resistance is calculated using hydrodynamic methods;S3, fracturing parameter, equivalent hole frictional resistance and nearly well bending friction are based on;S4, it is based on pressure balance in fracturing string, calculates the real-time bottom pressure of pressure break;S5, extension fracture and shear fracture calculation of pressure critical pressure are based on;S6, it is based on bottom pressure Numerical Discrimination storey increase design volume expansion validity;S7, range is involved based on injection process linear diffusion equation calculating storey increase design.Beneficial effects of the present invention:Volume transformation is calculated through the invention and involves range and SRV, economic and reliable has filled up the blank that pressing crack construction assesses the storey increase design volume prior art in real time.
Description
Technical field
The present invention relates to one kind to assess storey increase design volume method based on Fracturing Pressure in real time, belongs to exploration of oil and gas field
Develop hydraulic fracturing technology field.
Background technique
In recent years, huge success is achieved in North America with shale gas revolution, China is also accelerating shale gas exploitation
Paces.In October, 2015, Chinese World Mining Industry conference are announced, and after the U.S., Canada, China becomes third and realizes shale
The country of gas commercial development.It forms Fuling, Changning, Weiyuan, extend four big shale gas producing regions.China's shale gas produces within 2015
Amount is 44.71 billion cubic meters, embodies wide development prospect.Since the gas supply capacity in shale reservoir matrix to crack is poor,
Only it is difficult to obtain expected effect of increasing production by the main seam of traditional single pressure break, volume, which is transformed, to be formed fracture network and can just make fluid
" shortest distance " seepage flow is realized from matrix to crack.Therefore it there has been proposed low hole, hypotonic oil and gas reservoir is suitble to, improves reservoir and changes
Volume is made to increase " volume fracturing " technology of oil and gas production.The technology can not only increase substantially well yield, additionally it is possible to
The effective Lower Limit of Producing of Ultra-low of reservoir is reduced, improves reservoir use rate and recovery ratio to greatest extent.
Reservoir reconstruction volume (SRV) is used to describe the fracture network form of volume fracturing generation and involves range.Currently, meter
The method for calculating SRV mainly has micro-seismic monitoring method, dipmeter survey method and discrete fracture network simulation.Wherein micro seismic method compared with
To be reliable, but microsei-smic testing and dipmeter testing cost are high, less economical, and large-scale application is limited.Discrete fracture network
Simulation does not have generalization as it is assumed that condition causes practical calculating to differ greatly with simulation calculated result.
Summary of the invention
For the deficiencies in the prior art, the present invention provides a kind of independent of specially test, directly application
What pressing crack construction process, that is, retrievable pressure, discharge capacity parameter were directly explained assesses storey increase design based on Fracturing Pressure in real time
Volume method.
Technical scheme is as follows:
One kind assessing storey increase design volume method based on Fracturing Pressure in real time, which is characterized in that includes the following steps:
S1, target well reservoir geology and engineering parameter, casing programme data, fracturing material data, pump note program base are collected
This parameter;
S2, fracturing string flowing frictional resistance accounting equation is established, using hydrodynamic methods combination laboratory experiment test result
It calculates fluid in fracturing string and flows frictional resistance;
S3, fluid properties and construction parameter are injected based on real-time, calculates nearly well frictional resistance, acquires always rubbing for ground to crack body
Resistance;
S4, it is based on pressure balance in fracturing string, calculates the real-time bottom pressure of hydraulic fracturing;
S5, the extension fracture of structure plane of weakness and shear fracture calculation of pressure critical pressure are based on;
S6, it is based on bottom pressure Numerical Discrimination storey increase design volume expansion validity;
S7, range is involved based on injection process linear diffusion equation calculating storey increase design.
Further, target well reservoir geology and engineering parameter include well depth, vertical depth and horizontal segment length, target spot in the S1
Position, thickness, cluster spacing, Original strata stress position, the intrinsic fracture angle of intersection, cohesive force and internal friction angle, perforation hole count, aperture,
Eyelet discharge coefficient, preforation tunnel diameter, effective perforation hole count, the casing programme data include hole deviation, orientation, fracturing string
Diameter, the fracturing material include that fracturing fluid net liquid density, consistency coefficient, flow stance index, proppant partial size, proppant particles are close
Degree, proppant apparent density, drag reduction ratio, the peaceful coefficient of friction resistance of model, clear water density.
Further, fracturing string flowing frictional resistance accounting equation uses common resistance-reducing yield method meter in engineering in the S2
Calculate flowing frictional resistance of the fracturing working fluid in round tube;Based on flowing frictional resistance i.e. by clear water in round tube, pass through experiment test
The ratio of practical fracturing working fluid and clear water frictional resistance is modified, and flowing frictional resistance expression formula of the fracturing working fluid in tubing string is such as
Under:
Wherein:When fracturing working fluid is net liquid body:Δpfr=(1- λ) Δ pwater
When fracturing working fluid is mulling liquid:Δpfr=fCorrection(1-λ)Δpwater
In formula:ΔpfrThe flowing frictional resistance of-fracturing working fluid in fracturing string, MPa;Mono- drag reducing efficiency of λ, zero dimension;
fCorrection- mulling liquid frictional resistance correction coefficient, zero dimension;ΔpwaterThe flowing frictional resistance of-clear water in fracturing string, MPa;
Wherein clear water is as follows in the calculation expression of tubing string friction drag:
In formula:fwaterThe peaceful coefficient of friction resistance of-model, zero dimension;Flow velocity of the ν-clear water in tubing string, m/s;D-fracturing string
Internal diameter, m;ρwaterThe density of-clear water, kg/m3;L-fracturing string total length, m.
Wherein the calculation expression of the peaceful coefficient of friction resistance of model is as follows:
Nrer=106Dv
In formula:fwaterThe peaceful coefficient of friction resistance of-model, zero dimension;Flow velocity of the ν-clear water in tubing string, m/s;D-fracturing string
Internal diameter, m;Nre-Reynolds number, zero dimension.
When fracturing working fluid is mulling liquid in the pressure break pipe tubing string, then introduces mulling liquid zero dimension density and rub to mulling liquid
Resistance is modified.Wherein mulling liquid density, mulling liquid zero dimension density distinguish expression formula and are:
In formula:ρm- mulling liquid density, kg/m3;ρf- net liquid volume density, kg/m3;C-proppant bulk concentration, zero dimension;
ρpr- proppant particles density, kg/m3;ρps- proppant apparent density, kg/m3;ρD- mulling liquid zero dimension density, kg/m3。
In conjunction with the great mass of data that pressing crack construction scene enrolls, data processing is carried out using the method for fitting of a polynomial, is obtained
Mulling liquid frictional resistance correction coefficient:
fCorrection=-0.1448 ρD 2+0.1094ρD+1.1354。
Further, nearly well frictional resistance includes Kong Zhaoxia and nearly well bending friction, the meter of the Kong Zhaoxia in the S3
Operator expression formula is as follows:
In formula:Δpper- Kong Zhaoxia, MPa;ρ-fracturing working fluid density, kg/m3;Q-operational discharge capacity, m3/min;
Np- effective perforation hole count, zero dimension;Dp- preforation tunnel diameter, m;Cd- eyelet discharge coefficient;
After injecting load fluid, eyelet abrasion causes discharge coefficient to increase, the empirical equation expression of equivalent hole discharge coefficient
It is as follows:
In formula:M-passes through the proppant quality of eyelet, kg;
The nearly well bending friction drops discharge capacity by ladder and acquires, and expression formula is as follows:
In formula:Knear- nearly well crack tortuosity the coefficient of friction resistance, MPa/ (m3/min)1/2;Q-operational discharge capacity, m3/min;
Nearly well frictional resistance is further acquired, expression formula is as follows:
Δptotal=Δ pper+Δpnear
In formula:Δptotal- nearly well frictional resistance, MPa;Δpper- Kong Zhaoxia, MPa;Δpnear- nearly well bending friction,
MPa。
Further, the expression formula of the real-time bottom pressure of S4 hydraulic fracturing is as follows:
pinj=pt+ph-Δpfr-Δptotal
In formula:pinjThe real-time bottom pressure of-pressure break, MPa;ptThe real-time well head pressure of-pressing crack construction, MPa;ph- pressure break is applied
Work tubing string liquid column hydrostatic pressure, MPa;ΔpfrThe flowing frictional resistance of-fracturing working fluid in fracturing string, MPa;Δptotal- nearly well
Frictional resistance, MPa.
Further, the expression formula of the critical pressure based on structure plane of weakness extension fracture in the S5 is as follows:
In formula:ptrg- extensional destroys critical pressure, MPa;σmax- horizontal direction maximum principal stress, MPa;σmin- horizontal
Direction minimum principal stress, MPa;
Expression of critical pressure based on shear fracture is as follows:
In formula:ptrg- failure by shear critical pressure, MPa;σmax- horizontal direction maximum principal stress, MPa;σmin- horizontal
Direction minimum principal stress, MPa;KfThe tangent value of-angle of friction, dimensionless;τ0- cohesive force, MPa.
Further, Numerical Discrimination storey increase design volume expansion validity is based in the S6, when real-time bottom pressure is small
In previous moment bottom pressure, illustrate after the extension of crack that storey increase design volume maintains not because fluid filling causes pressure not increase
Become;Otherwise new bottom-hole pressure reservoir reconstruction volume is pressed.
Further, linear diffusion equation is injected in the S7 using single hop perforation source as coordinate origin, establishes three
Cartesian coordinate system is tieed up, wherein X-direction is main fracture extension direction, and Y-direction is parallel to pit shaft direction, and Z-direction is flat perpendicular to XY
Upwardly.Hydraulic fracture is obtained according to diffusion equation and extends maximum distance on X, Y, Z different directions, and is calculated and increased by spheroid
Produce transformation volume:
Wherein Δ ptrg=min (ptrg1,,ptrg2,)-piΔpinj=pinj-pi
In formula:ΔptrgThe minimum pressure increment that-intrinsic fracture destroys, MPa;Δpinj- bottom pressure increment, MPa;
pi- reservoir initial pressure, MPa;ηx, ηy, ηz- diffusion coefficient, m2/s。kx, ky, kz- permeability tensor principal direction component,
mD;- effecive porosity, dimensionless;μ-liquid viscosity, mPas;Ctol- stratum total compression coefficient, MPa-1;T-time,
s;X, y, z-coordinate components, m.
Beneficial effects of the present invention are:
The present invention is based on diffusion effect of the crack fluid pressure in reservoir in shale volume fracturing process, in conjunction with rock
Failure criteria, proposes the computation model and method of a kind of Fast Evaluation calculating SRV, and the present invention considers that volume fracturing injection is more
The complexity of seed type fluid and the switching of material real-time, interactive is rubbed in conjunction with fracturing string fluid flowing frictional resistance computation model, eyelet
Resistance calculates with bending friction model and obtains real-time net pressure, compages plane of weakness extensional and shear fracture calculation of pressure model and stream
The linear diffusion equation of body calculates assessment storey increase design volume in real time;The present invention can be implemented optimization volume transformation scale and be based on certain page
Rock example well parameter, the model calculation and microseism cloud atlas are compared, the reliability of model in text is demonstrated.
Detailed description of the invention
It, below will be to use required in embodiment in order to illustrate more clearly of the technical solution of embodiment of the present invention
Attached drawing be briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not to be seen as
It is the restriction to range, it for those of ordinary skill in the art, without creative efforts, can be with root
Other relevant attached drawings are obtained according to these attached drawings.
Discharge capacity schematic diagram drops in Fig. 1 ladder provided by the invention;
Fig. 2 fracturing fluid provided by the invention flows frictional resistance calculated result figure;
Fig. 3 preforation tunnel frictional resistance calculated result figure provided by the invention;
Fig. 4 eyelet discharge coefficient change curve provided by the invention;
Fig. 5 structure plane of weakness extension fracture provided by the invention and shear fracture calculation of pressure critical pressure are tried hard to;
Fig. 6 X well section microseism cloud atlas provided by the invention;
Fig. 7 model provided by the invention calculates SRV schematic three dimensional views.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples.
As shown in Figure 1, a kind of assess storey increase design volume method based on Fracturing Pressure in real time, include the following steps:
S1, target well reservoir geology and engineering parameter, casing programme data, fracturing material data, pump note program base are collected
This parameter;
The target well reservoir geology and engineering parameter include well depth, vertical depth and horizontal segment length, target position, thickness, cluster
Spacing, Original strata stress position, the intrinsic fracture angle of intersection, cohesive force and internal friction angle, perforation hole count, aperture, eyelet discharge coefficient,
Preforation tunnel diameter, effective perforation hole count, the casing programme data include hole deviation, orientation, fracturing string diameter, the pressure break
Material include fracturing fluid net liquid density, consistency coefficient, flow stance index, proppant partial size, proppant particles density, proppant view it is close
The peaceful coefficient of friction resistance of degree, drag reducing efficiency, model, clear water density.
S2, fracturing string flowing frictional resistance accounting equation is established, using hydrodynamic methods combination laboratory experiment test result
It calculates fluid in fracturing string and flows frictional resistance;
The fracturing string flowing frictional resistance accounting equation calculates fracturing working fluid using drag reducing efficiency method common in engineering
Flowing frictional resistance in round tube;Based on flowing frictional resistance i.e. by clear water in round tube, passes through experiment and test practical fracture job
The ratio of liquid and clear water frictional resistance is modified, and flowing frictional resistance expression formula of the fracturing working fluid in tubing string is as follows:
Wherein:When fracturing working fluid is net liquid body:
Δpfr=(1- λ) Δ pwater (1)
When fracturing working fluid is mulling liquid:
Δpfr=fCorrection(1-λ)Δpwater (2)
In formula:ΔpfrThe flowing frictional resistance of-fracturing working fluid in fracturing string, MPa;Mono- drag reducing efficiency of λ, zero dimension;
fCorrection- mulling liquid frictional resistance correction coefficient, zero dimension;ΔpwaterThe flowing frictional resistance of-clear water in fracturing string, MPa;
Wherein clear water is as follows in the calculation expression of tubing string friction drag:
In formula:fwaterThe peaceful coefficient of friction resistance of-model, zero dimension;Flow velocity of the ν-clear water in tubing string, m/s;D-fracturing string
Internal diameter, m;ρwaterThe density of-clear water, kg/m3;L-fracturing string total length, m.
Wherein the calculation expression of the peaceful coefficient of friction resistance of model is as follows:
Nrer=106Dv
In formula:fwaterThe peaceful coefficient of friction resistance of-model, zero dimension;Flow velocity of the ν-clear water in tubing string, m/s;D-fracturing string
Internal diameter, m;Nre-Reynolds number, zero dimension.
When fracturing working fluid is mulling liquid in the pressure break pipe tubing string, then introduces mulling liquid zero dimension density and rub to mulling liquid
Resistance is modified.Wherein mulling liquid density, mulling liquid zero dimension density distinguish expression formula and are:
In formula:ρm- mulling liquid density, kg/m3;ρf- net liquid volume density, kg/m3;C-proppant bulk concentration, zero dimension;
ρpr- proppant particles density, kg/m3;ρps- proppant apparent density, kg/m3;ρD- mulling liquid zero dimension density, kg/m3。
In conjunction with the great mass of data that pressing crack construction scene enrolls, data processing is carried out using the method for fitting of a polynomial, is obtained
Mulling liquid frictional resistance correction coefficient:
fCorrection=-0.1448 ρD 2+0.1094ρD+1.1354 (6)
S3, fluid properties and construction parameter are injected based on real-time, calculates nearly well frictional resistance, acquires always rubbing for ground to crack body
Resistance;
The nearly well frictional resistance includes Kong Zhaoxia and nearly well bending friction, and the calculation expression of the Kong Zhaoxia is as follows:
In formula:Δpper- Kong Zhaoxia, MPa;ρ-fracturing working fluid density, kg/m3;Q-operational discharge capacity, m3/min;
Np- effective perforation hole count, zero dimension;Dp- preforation tunnel diameter, m;Cd- eyelet discharge coefficient;
After injecting load fluid, eyelet abrasion causes discharge coefficient to increase, the empirical equation expression of equivalent hole discharge coefficient
It is as follows:
In formula:M-passes through the proppant quality of eyelet, kg;
It is acquired as shown in Fig. 2, the nearly well bending friction drops discharge capacity by ladder, expression formula is as follows:
In formula:Knear- nearly well crack tortuosity the coefficient of friction resistance, MPa/ (m3/min)1/2;Q-operational discharge capacity, m3/min;
The nearly total frictional resistance of well is further acquired, expression formula is as follows:
Δptotal=Δ pper+Δpnear (10)
In formula:Δptotal- nearly well frictional resistance, MPa;Δpper- Kong Zhaoxia, MPa;Δpnear- nearly well bending friction,
MPa。
S4, it is based on the conservation of energy, calculates the real-time bottom pressure of hydraulic fracturing;
The expression formula of the real-time bottom pressure of hydraulic fracturing is as follows:
pinj=pt+ph-Δpfr-Δptotal (11)
In formula:pinjThe real-time bottom pressure of-pressure break, MPa;ptThe real-time well head pressure of-pressing crack construction, MPa;ph- pressure break is applied
Work tubing string liquid column hydrostatic pressure, MPa;ΔpfrThe flowing frictional resistance of-fracturing working fluid in fracturing string, MPa;Δptotal- nearly well
Frictional resistance, MPa
S5, the extension fracture of structure plane of weakness and shear fracture calculation of pressure critical pressure are based on;
The expression formula of the critical pressure based on structure plane of weakness extension fracture is as follows:
In formula:ptrg- extensional destroys critical pressure, MPa;σmax- horizontal direction maximum principal stress, MPa;σmin- horizontal
Direction minimum principal stress, MPa;
Expression of critical pressure based on shear fracture is as follows:
In formula:ptrg- failure by shear critical pressure, MPa;σmax- horizontal direction maximum principal stress, MPa;σmin- horizontal
Direction minimum principal stress, MPa;Kf- coefficient of internal friction, dimensionless;τ0- cohesive force, MPa.
S6, it is based on bottom pressure Numerical Discrimination storey increase design volume expansion validity;
It is described to be based on Numerical Discrimination storey increase design volume expansion validity, when real-time bottom pressure is less than previous moment shaft bottom
Pressure, illustrate crack extension after because fluid filling cause pressure not increase, storey increase design volume remains unchanged;Otherwise new well is pressed
Bottom pressure calculates reservoir reconstruction volume.
S7, range is involved based on injection process linear diffusion equation calculating storey increase design;
The injection linear diffusion equation establishes three-dimensional cartesian coordinate using single hop perforation source as coordinate origin
System, wherein X-direction is main fracture extension direction, and Y-direction is parallel to pit shaft direction, and Z-direction is upward perpendicular to X/Y plane.Assuming that storage
Layer is homogeneous anisotropy poroelasticity medium, and formation fluid is gas phase, and nearly well bending friction is ignored, and injects fluid not shadow
Gas reservoir pore volume is rung, expression formula is as follows:
In formula:- formation pore pressure, MPa;kx, ky, kzComponent of-the permeability tensor in principal direction, mD;- have
Imitate porosity, dimensionless;μ-liquid viscosity, mPas;Ctol- stratum total compression coefficient, MPa-1;T-time, s;X, y, z-
Coordinate components, m;
Primary condition is with boundary condition:
P=pi, 0≤x≤∞, 0≤y≤∞, 0≤z≤∞, t=0 (15)
P=pinj, x=0, y=0, z=0, t > 0 (16)
In formula:pi- reservoir initial pressure, MPa;pinj- shaft bottom pumps injection pressure power, MPa;
Further acquire:
In formula:ΔpresCertain point pressure increases numerical value, MPa in-reservoir;Δpinj- bottom pressure increases numerical value, MPa;
ηx,ηy,ηz- diffusion coefficient, m2/s;
Wherein:
As reservoir point pressure amplification Δ presThe critical pressure amplification Δ p opened greater than intrinsic fracturetrgWhen, day at the point
Right crack will be destroyed;Assuming that the intrinsic fracture mechanical property of reservoir everywhere is identical, then intrinsic fracture destroys in reservoir
Condition it is same, then the Δ p inside reservoirtrgIt is equal everywhere, then:
In formula:Δptrg- it is the minimum pressure increment that intrinsic fracture destroys;
The SRV extension maximum distance on different directions is finally calculated:
Wherein Δ ptrg=min (ptrg1,,ptrg2,)-piΔpinj=pinj-pi (23)
In formula:ΔptrgThe minimum pressure increment that-intrinsic fracture destroys, MPa;Δpinj- bottom pressure increment, MPa;
pi- reservoir initial pressure, MPa;ηx, ηy, ηz- diffusion coefficient, m2/s。kx, ky, kz- permeability tensor principal direction component,
mD;- effecive porosity, dimensionless;μ-liquid viscosity, mPas;Ctol- stratum total compression coefficient, MPa-1;T-time,
s;X, y, z-coordinate components, m.
Embodiment:
The X well section of selection Southwestern China shale exploitation demonstration area, horizontal well completion vertical depth 4276m, horizontal segment length 800m,
Effective thickness 32m, formation fluid pressure 42MPa.It is respectively 2.1 × 10 according to core experiment analysis level and vertical permeability-3μ
m2、4.9×10-3μm2、2.8×10-3μm2, rock system compressibility 1 × 10-4MPa-1;It is vertical to be respectively with horizontal principal stress
145MPa,121.6MPa,105.0MPa;25 ° of intrinsic fracture cohesive force 17.5MPa, internal friction angle, 20 ° of the angle of intersection;Average tensile
Intensity 6.56MPa.Now use the smooth casing of 177.8 δ of φ 12.65, slippery water fracturing fluid (flow pattern index is 1.0, drag reducing efficiency 70%)
Divide 15 sections, 3 clusters, crack cluster spacing 15m;Effective 20 hole of eyelet of perforation, aperture φ 12;Preamble fracturing fracture height 50m.Pressure break is applied
Work injection 2hr amounts to total net liquid amount 1700m3, total to support dosage 60m3(40/70 mesh, grain density 2800Kg/m3, apparent density
1800Kg/m3), discharge capacity 15m3/min.Well head operation pressure is 90MPa when fluid injection 80min, sand are than 5%, determines that volume increase changes at this time
Make volume.
As shown in Fig. 2, being calculated by formula (4), mulling liquid density of the sand than 5% is calculated as 1058kg/m3, pass through
Net liquid frictional resistance (this well is sliding to slip water, friction drag 1.5MPa/km using drag reducing efficiency 70%) is calculated in formula (1), passes through public affairs
It is 8.29MPa that mulling liquid frictional resistance calculated result mulling liquid frictional resistance, which is calculated, in formula (2), wherein nothing is calculated by formula (5)
Dimension density is 1.058, and f is calculated by formula (6)Correction=1.089.
As shown in figs. 34, preforation tunnel frictional resistance 1.3MPa is calculated by formula (8), is calculated by formula 10
Nearly well bending friction is 2.68MPa.
It is 120.48MPa by the real-time bottom pressure that hydraulic fracturing nearly well crack end is calculated in formula (12).
As shown in figure 5, structure plane of weakness, which is calculated, by formula (13), formula (14) generates extensional destruction and failure by shear
The critical pressure needed.
As shown in fig. 6, according to the micro-seismic monitoring in work progress, as shown in Figure 1.(t=2h), reservoir at the end of construction
The a length of 125m of volume average crack half is transformed, laterally involving range is 120m, and average crack height is 100m, and transformation volume is
157.08×104m3。
As shown in fig. 6, three shaft lengths according to spheroid are drawn in conjunction with the microseism data monitored in work progress
Three-dimensional S RV6;As can be seen from Figure, the calculated three-dimensional S RV of model and the fitting of microseism data are preferable, and accuracy is higher.
The above is not intended to limit the present invention in any form, although the present invention is disclosed such as by embodiment
On, however, it is not intended to limit the invention, any person skilled in the art, is not departing from technical solution of the present invention range
It is interior, when the technology contents using the disclosure above make a little change or are modified to the equivalent embodiment of equivalent variations, but it is all
It is to the above embodiments according to the technical essence of the invention any simply to repair without departing from the content of technical solution of the present invention
Change, equivalent variations and modification, all of which are still within the scope of the technical scheme of the invention.
Claims (8)
1. one kind assesses storey increase design volume method based on Fracturing Pressure in real time, which is characterized in that include the following steps:
S1, collection target well reservoir geology and engineering parameter, casing programme data, fracturing material data, pump note program are joined substantially
Number;
S2, fracturing string flowing frictional resistance accounting equation is established, using hydrodynamic methods combination laboratory experiment test result calculations
Fluid flows frictional resistance in fracturing string;
S3, fluid properties and construction parameter are injected based on real-time, calculates nearly well frictional resistance, acquires total frictional resistance of ground to crack body;
S4, it is based on pressure balance in fracturing string, calculates the real-time bottom pressure of hydraulic fracturing;
S5, the extension fracture of structure plane of weakness and shear fracture calculation of pressure critical pressure are based on;
S6, it is based on bottom pressure Numerical Discrimination storey increase design volume expansion validity;
S7, range is involved based on injection process linear diffusion equation calculating storey increase design.
2. a kind of Fracturing Pressure that is based on according to claim 1 assesses storey increase design volume method, feature in real time
It is, target well reservoir geology and engineering parameter include well depth, vertical depth and horizontal segment length, target position, thickness, cluster in the S1
Spacing, Original strata stress position, the intrinsic fracture angle of intersection, cohesive force and internal friction angle, perforation hole count, aperture, eyelet discharge coefficient,
Preforation tunnel diameter, effective perforation hole count, the casing programme data include hole deviation, orientation, fracturing string diameter, the pressure break
Material include fracturing fluid net liquid density, consistency coefficient, flow stance index, proppant partial size, proppant particles density, proppant view it is close
Degree, drag reducing efficiency, clear water density.
3. a kind of Fracturing Pressure that is based on according to claim 3 assesses storey increase design volume method, feature in real time
It is, fracturing string flowing frictional resistance accounting equation calculates fracturing working fluid using drag reducing efficiency method common in engineering in the S2
Flowing frictional resistance in round tube;Based on flowing frictional resistance i.e. by clear water in round tube, passes through experiment and test practical fracture job
The ratio of liquid and clear water frictional resistance is modified, and flowing frictional resistance expression formula of the fracturing working fluid in tubing string is as follows:
Wherein:When fracturing working fluid is net liquid body:Δpfr=(1- λ) Δ pwater
When fracturing working fluid is mulling liquid:Δpfr=fCorrection(1-λ)Δpwater
In formula:ΔpfrThe flowing frictional resistance of-fracturing working fluid in fracturing string, MPa;Mono- drag reducing efficiency of λ, zero dimension;fCorrection- mixed
Sand liquid frictional resistance correction coefficient, zero dimension;ΔpwaterThe flowing frictional resistance of-clear water in fracturing string, MPa;
Wherein clear water is as follows in the calculation expression of tubing string friction drag:
In formula:fwaterThe peaceful coefficient of friction resistance of-model, zero dimension;Flow velocity of the ν-clear water in tubing string, m/s;D-fracturing string internal diameter,
m;ρwaterThe density of-clear water, kg/m3;L—.Fracturing string total length, m.
Wherein the calculation expression of the peaceful coefficient of friction resistance of model is as follows:
Nrer=106Dv
In formula:fwaterThe peaceful coefficient of friction resistance of-model, zero dimension;Flow velocity of the ν-clear water in tubing string, m/s;D-fracturing string internal diameter,
m;Nre-Reynolds number, zero dimension.
In the pressure break pipe tubing string fracturing working fluid be mulling liquid when, then introduce mulling liquid zero dimension density to mulling liquid frictional resistance into
Row amendment.Wherein mulling liquid density, mulling liquid zero dimension density distinguish expression formula and are:
In formula:ρm- mulling liquid density, kg/m3;ρf- net liquid volume density, kg/m3;C-proppant bulk concentration, zero dimension;ρpr—
Proppant particles density, kg/m3;ρps- proppant apparent density, kg/m3;ρD- mulling liquid zero dimension density, kg/m3。
In conjunction with the great mass of data that pressing crack construction scene enrolls, data processing is carried out using the method for fitting of a polynomial, obtains mulling
Liquid frictional resistance correction coefficient:
fCorrection=-0.1448 ρD 2+0.1094ρD+1.1354
In formula:fCorrection- mulling liquid frictional resistance correction coefficient, zero dimension;ρD- mulling liquid zero dimension density, kg/m3。
4. a kind of Fracturing Pressure that is based on according to claim 1 assesses storey increase design volume method, feature in real time
It is, nearly well frictional resistance includes Kong Zhaoxia and nearly well bending friction in the S3, and the calculation expression of the Kong Zhaoxia is as follows:
In formula:Δpper- Kong Zhaoxia, MPa;ρ-fracturing working fluid density, kg/m3;Q-operational discharge capacity, m3/min;Np- have
Imitate perforation hole count, zero dimension;Dp- preforation tunnel diameter, m;Cd- eyelet discharge coefficient;
After injecting mulling liquid, eyelet abrasion causes discharge coefficient to increase, and the empirical equation of equivalent hole discharge coefficient is expressed as follows:
In formula:M-passes through the proppant quality of eyelet, kg;
The nearly well bending friction drops discharge capacity by ladder and acquires, and expression formula is as follows:
In formula:Knear- nearly well crack tortuosity the coefficient of friction resistance, MPa/ (m3/min)1/2;Q-operational discharge capacity, m3/min;
The nearly total frictional resistance of well is further acquired, expression formula is as follows:
Δptotal=Δ pper+Δpnear
In formula:Δptotal- nearly well frictional resistance, MPa;Δpper- Kong Zhaoxia, MPa;Δpnear- nearly well bending friction, MPa.
5. a kind of Fracturing Pressure that is based on according to claim 1 assesses storey increase design volume method, feature in real time
It is, the expression formula of the real-time bottom pressure of S4 hydraulic fracturing is as follows:
pinj=pt+ph-Δpfr-Δptotal
In formula:pinjThe real-time bottom pressure of-pressure break, MPa;ptThe real-time well head pressure of-pressing crack construction, MPa;ph- pressing crack construction pipe
Column liquid column hydrostatic pressure, MPa;ΔpfrThe flowing frictional resistance of-fracturing working fluid in fracturing string, MPa;Δptotal- closely well rubs
Resistance, MPa.
6. a kind of Fracturing Pressure that is based on according to claim 1 assesses storey increase design volume method, feature in real time
It is, the expression formula of the critical pressure based on structure plane of weakness extension fracture in the S5 is as follows:
In formula:ptrg-1- extensional destroys critical pressure, MPa;σmax- horizontal direction maximum principal stress, MPa;σmin- horizontal direction
Minimum principal stress, MPa.
Expression of critical pressure based on shear fracture is as follows:
In formula:ptrg-2- failure by shear critical pressure, MPa;σmax- horizontal direction maximum principal stress, MPa;σmin- horizontal direction
Minimum principal stress, MPa;Kf- coefficient of internal friction, dimensionless;τ0- cohesive force, MPa.
7. a kind of Fracturing Pressure that is based on according to claim 1 assesses storey increase design volume method, feature in real time
It is, Numerical Discrimination storey increase design volume expansion validity is based in the S6, when real-time bottom pressure is less than previous moment well
Bottom pressure, illustrate crack extension after because fluid filling cause pressure not increase, storey increase design volume remains unchanged;Otherwise it presses newly
Bottom-hole pressure reservoir reconstruction volume.
8. a kind of Fracturing Pressure that is based on according to claim 1 assesses storey increase design volume method, feature in real time
It is, linear diffusion equation is injected in the S7 using single hop perforation source as coordinate origin, establishes three-dimensional cartesian coordinate
System, wherein X-direction is main fracture extension direction, and Y-direction is parallel to pit shaft direction, and Z-direction is upward perpendicular to X/Y plane.According to expansion
Scattered equation obtains hydraulic fracture and extends maximum distance on X, Y, Z different directions, and calculates storey increase design volume by spheroid:
Wherein Δ ptrg=min (ptrg1,ptrg2)-pi Δpinj=pinj-pi
In formula:ΔptrgThe minimum pressure increment that-intrinsic fracture destroys, MPa;Δpinj- bottom pressure increment, MPa;pi- storage
Layer initial pressure, MPa;ηx, ηy, ηz- diffusion coefficient, m2/s。kx, ky, kzComponent of-the permeability tensor in principal direction, mD;
- effecive porosity, dimensionless;μ-liquid viscosity, mPas;Ctol- stratum total compression coefficient, MPa-1;T-time, s;X,
Y, z-coordinate components, m.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810530032.0A CN108829945A (en) | 2018-05-29 | 2018-05-29 | One kind assessing storey increase design volume method based on Fracturing Pressure in real time |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810530032.0A CN108829945A (en) | 2018-05-29 | 2018-05-29 | One kind assessing storey increase design volume method based on Fracturing Pressure in real time |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108829945A true CN108829945A (en) | 2018-11-16 |
Family
ID=64146157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810530032.0A Pending CN108829945A (en) | 2018-05-29 | 2018-05-29 | One kind assessing storey increase design volume method based on Fracturing Pressure in real time |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108829945A (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110017943A (en) * | 2019-05-05 | 2019-07-16 | 洛阳理工学院 | A kind of error correcting method of discrete media pressure measurement |
CN110056336A (en) * | 2019-05-31 | 2019-07-26 | 西南石油大学 | A kind of shale air cleft network pressure splits operation pressure curve automatic diagnosis method |
CN110080739A (en) * | 2019-05-08 | 2019-08-02 | 中国石油天然气股份有限公司 | A kind of straight well multiple-limb seam pressure break site operation judgment method |
CN110469303A (en) * | 2019-07-04 | 2019-11-19 | 西南石油大学 | A kind of volume fracturing method for optimally designing parameters based on four classes transformation volume |
CN110905472A (en) * | 2019-10-29 | 2020-03-24 | 中国石油集团川庆钻探工程有限公司 | Method for determining real-time steering fracturing parameters based on composite temporary plugging system |
CN111734382A (en) * | 2020-08-05 | 2020-10-02 | 西南石油大学 | Method for explaining multiple parameters by testing fracturing through stepped displacement reduction |
CN111980698A (en) * | 2020-09-30 | 2020-11-24 | 西南石油大学 | Deep high-temperature shale formation fracture pressure calculation method considering thermal damage |
CN112145167A (en) * | 2020-11-10 | 2020-12-29 | 西南石油大学 | Fracturing dynamic bottom hole pressure calculation method based on construction data |
CN112347599A (en) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
CN116733432A (en) * | 2023-05-09 | 2023-09-12 | 陕西省煤层气开发利用有限公司 | Transformation process method for deep coal bed gas fracturing |
CN116877067A (en) * | 2023-07-18 | 2023-10-13 | 重庆地质矿产研究院 | Method for predicting hydraulic fracturing generated cracks and swept area fluid pressure |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107742020A (en) * | 2017-10-09 | 2018-02-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Shale gas reservoir fracturing improvement method for predicting volume |
CN108021728A (en) * | 2016-10-28 | 2018-05-11 | 中国石油化工股份有限公司 | A kind of computational methods of fracturing fluid recovery (backflow) system |
-
2018
- 2018-05-29 CN CN201810530032.0A patent/CN108829945A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021728A (en) * | 2016-10-28 | 2018-05-11 | 中国石油化工股份有限公司 | A kind of computational methods of fracturing fluid recovery (backflow) system |
CN107742020A (en) * | 2017-10-09 | 2018-02-27 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Shale gas reservoir fracturing improvement method for predicting volume |
Non-Patent Citations (4)
Title |
---|
任岚等: "《基于各向异性扩散方程的页岩气井改造体积计算模型》", 《大庆石油地质与开发》 * |
刘登峰: "《煤层气井压裂施工资料反演岩石力学参数及压后产能预测研究》", 《中国优秀硕士学位论文全文数据库(工程科技I辑)》 * |
张彦杰等: "《限流压裂过程中射孔孔眼摩阻变化研究》", 《勘探开发》 * |
路智勇: "《压裂施工动态摩阻模型建立及其敏感性分析》", 《断块油气田》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110017943A (en) * | 2019-05-05 | 2019-07-16 | 洛阳理工学院 | A kind of error correcting method of discrete media pressure measurement |
CN110017943B (en) * | 2019-05-05 | 2021-02-26 | 洛阳理工学院 | Error correction method for discrete medium pressure measurement |
CN110080739A (en) * | 2019-05-08 | 2019-08-02 | 中国石油天然气股份有限公司 | A kind of straight well multiple-limb seam pressure break site operation judgment method |
CN110080739B (en) * | 2019-05-08 | 2021-03-26 | 中国石油天然气股份有限公司 | Method for judging construction site of vertical well multi-branch fracture |
CN110056336A (en) * | 2019-05-31 | 2019-07-26 | 西南石油大学 | A kind of shale air cleft network pressure splits operation pressure curve automatic diagnosis method |
CN110469303A (en) * | 2019-07-04 | 2019-11-19 | 西南石油大学 | A kind of volume fracturing method for optimally designing parameters based on four classes transformation volume |
CN112347599A (en) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
CN110905472A (en) * | 2019-10-29 | 2020-03-24 | 中国石油集团川庆钻探工程有限公司 | Method for determining real-time steering fracturing parameters based on composite temporary plugging system |
CN110905472B (en) * | 2019-10-29 | 2021-10-22 | 中国石油集团川庆钻探工程有限公司 | Method for determining real-time steering fracturing parameters based on composite temporary plugging system |
CN111734382A (en) * | 2020-08-05 | 2020-10-02 | 西南石油大学 | Method for explaining multiple parameters by testing fracturing through stepped displacement reduction |
CN111734382B (en) * | 2020-08-05 | 2021-06-01 | 西南石油大学 | Method for explaining multiple parameters by testing fracturing through stepped displacement reduction |
CN111980698A (en) * | 2020-09-30 | 2020-11-24 | 西南石油大学 | Deep high-temperature shale formation fracture pressure calculation method considering thermal damage |
CN111980698B (en) * | 2020-09-30 | 2022-03-18 | 西南石油大学 | Deep high-temperature shale formation fracture pressure calculation method considering thermal damage |
CN112145167A (en) * | 2020-11-10 | 2020-12-29 | 西南石油大学 | Fracturing dynamic bottom hole pressure calculation method based on construction data |
CN116733432A (en) * | 2023-05-09 | 2023-09-12 | 陕西省煤层气开发利用有限公司 | Transformation process method for deep coal bed gas fracturing |
CN116877067A (en) * | 2023-07-18 | 2023-10-13 | 重庆地质矿产研究院 | Method for predicting hydraulic fracturing generated cracks and swept area fluid pressure |
CN116877067B (en) * | 2023-07-18 | 2024-03-12 | 重庆地质矿产研究院 | Method for predicting hydraulic fracturing generated cracks and swept area fluid pressure |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108829945A (en) | One kind assessing storey increase design volume method based on Fracturing Pressure in real time | |
CN110905472B (en) | Method for determining real-time steering fracturing parameters based on composite temporary plugging system | |
CN108868748B (en) | Method for calculating repeated fracturing fracture opening pressure of shale gas horizontal well | |
Lian et al. | A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells | |
CN102925133B (en) | Fracturing fluid for controlling extension of fracture and fracturing method | |
CN107705215B (en) | A kind of shale reservoir refracturing selects well selections method | |
CN109958416B (en) | Multi-cluster perforation fracturing method for uniformly feeding liquid and sand with variable aperture and variable pore density | |
CN103867184B (en) | A kind of gas well critical liquid carrying flow rate determines method and device | |
CN109931045B (en) | Self-supporting acid fracturing method of double-seam system | |
CN107605451B (en) | A kind of ladder discharge capacity fracturing pump injecting method based on combined perforation | |
CN104234677B (en) | A kind of vertical displacement of gas injection improves gas condensate reservoir condensate recovery ratio method | |
CN103437746A (en) | Volume fracturing method for multiple cracks in multiple sections of horizontal well | |
CN109958411B (en) | Horizontal well cluster perforation staged fracturing method | |
CN103089224A (en) | Fracturing method for comprehensively controlling fracture height | |
CN107545088A (en) | A kind of normal pressure shale gas horizontal well volume fracturing method | |
CN112253066B (en) | Method for improving fracture complexity and transformation volume of deep shale gas | |
CN106545336A (en) | Consider the Productivity of tight gas reservoir seepage flow mechanism | |
CN105604534A (en) | Hydraulically affected fracturing process method for increasing production of coal-bed gas reservoir | |
CN110984949B (en) | Shale continuous sand-adding fracturing process | |
CN109323971A (en) | A kind of rock stratum bottom plate seepage flow gushing water experimental rig and test method | |
Zhao et al. | Experimental investigation on the fracture propagation of three-stage acid fracturing of tight sandstone gas reservoir | |
CN104234682B (en) | A kind ofly be applicable to separate stratum fracturing of continuous oil pipe method that is many, girdle | |
CN109187321B (en) | Large-size salt level communication well group water-soluble cavity-making simulation experiment method | |
CN103048184A (en) | Testing method of breakthrough pressure of repeated fracturing plugging agent | |
CN114592840B (en) | Temporary plugging fracturing method and application thereof |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181116 |