CN109138949A - A kind of cloth acid method based on horizontal well skin factor distribution characteristics - Google Patents
A kind of cloth acid method based on horizontal well skin factor distribution characteristics Download PDFInfo
- Publication number
- CN109138949A CN109138949A CN201811026178.8A CN201811026178A CN109138949A CN 109138949 A CN109138949 A CN 109138949A CN 201811026178 A CN201811026178 A CN 201811026178A CN 109138949 A CN109138949 A CN 109138949A
- Authority
- CN
- China
- Prior art keywords
- horizontal well
- well
- unit
- horizontal
- 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
- 239000002253 acid Substances 0.000 title claims abstract description 50
- 239000004744 fabric Substances 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005553 drilling Methods 0.000 claims abstract description 84
- 230000006378 damage Effects 0.000 claims abstract description 79
- 238000004519 manufacturing process Methods 0.000 claims abstract description 38
- 230000035699 permeability Effects 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 15
- 239000012224 working solution Substances 0.000 claims description 9
- 239000000243 solution Substances 0.000 claims description 8
- 208000027418 Wounds and injury Diseases 0.000 claims description 7
- 239000012530 fluid Substances 0.000 claims description 6
- 208000014674 injury Diseases 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 239000000443 aerosol Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 6
- 230000020477 pH reduction Effects 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002699 waste material Substances 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/16—Enhanced recovery methods for obtaining hydrocarbons
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a kind of cloth acid methods based on horizontal well skin factor distribution characteristics, comprising the following steps: 1), establishes individual well geological model under primary condition；2) the individual well geological model, based on step 1) building, the individual well geological model after calculating drilling damage；3) the individual well geological model, based on step 1) building, the horizontal well production profile before calculating drilling damage；4), individual well geological model after the drilling damage being calculated based on step 2), the horizontal well production profile after calculating drilling damage；5), the horizontal well production profile after the drilling damage that the horizontal well production profile before the drilling damage obtained according to step 3) and step 4) obtain calculates the skin factor on entire net horizontal section different location；6) the cloth acid amount on horizontal well different location, is calculated along the distribution characteristics of net horizontal section according to skin factor.The present invention solves the problems, such as that acidification low efficiency caused by existing blind cloth acid, acidizing effect are poor.
Description
Technical field
The present invention relates to oil-gas field development technical fields, and in particular to a kind of based on horizontal well skin factor distribution characteristics
Cloth acid method.
Background technique
Oil reservoir drainage area has been significantly greatly increased in the casing programme feature of horizontal well, with the incomparable advantage of straight well.But
It faces that reservoir heterogeneity is strong, different well section pollution differences are big, production capacity is unable to get the utmostly engineerings skill such as release simultaneously
Art problem.
Acidification measure is to solve depollution, give full play to horizontal well potential, realize the effective means that fast run-up produces.But due to water
Horizontal well section is generally longer, and the means such as current existing well testing, dynamic analysis and well log interpretation can not be obtained along entire horizontal well
The damage character parameter distribution situation (including: the parameters such as porosity, permeability, skin factor) of section, causes conventional stage to be acidified
Or large-scale general acidification does not have specific aim, acidizing effect is undesirable.
The purpose of acidification is to release reservoir damage, restores reservoir production capacity.Do not have for reservoir local on net horizontal section
The position to come to harm, blind cloth acid not only cannot effectively restore reservoir productivity, while can also aggravate the waste of acid solution, to lead
It causes under limited acid solution dosage, horizontal well acidizing effect is had a greatly reduced quality.
Summary of the invention
The purpose of the present invention is to provide a kind of cloth acid methods based on horizontal well skin factor distribution characteristics, solve existing
Low efficiency, the problem of acidizing effect difference are acidified caused by blind cloth acid.
The present invention is achieved through the following technical solutions:
A kind of cloth acid method based on horizontal well skin factor distribution characteristics, comprising the following steps:
1) the individual well geological model under primary condition, is established；
2) the individual well geological model, based on step 1) building, the individual well geological model after calculating drilling damage；
3) the individual well geological model, based on step 1) building, the horizontal well production profile before calculating drilling damage；
4), individual well geological model after the drilling damage being calculated based on step 2), the horizontal well after calculating drilling damage
Production profile；
5), after the drilling damage that the horizontal well production profile before the drilling damage obtained according to step 3) and step 4) obtain
Horizontal well production profile, calculate skin factor on entire net horizontal section different location；
6), the skin factor obtained according to step 5) calculates on horizontal well different location along the distribution characteristics of net horizontal section
Cloth acid amount.
The present invention carries out quantitative assessment by injuring to horizontal well drilling, obtains the skin factor on entire net horizontal section
Distribution characteristics, and using skin factor as design considerations, optimize the cloth acid amount on entire horizontal well different location, to improve acid
Liquid utilization rate, it is whole to improve horizontal well acidizing effect.
Further, individual well geological model described in step 1) is the grid system established under polar coordinates, and well logging is solved
The reservoir permeability and porosity attribute for releasing acquisition are assigned to grid system by the following method:
The assignment model of grid system porosity:
φ(r,θ,z,t0)=φ0(z)+φ0(z)γG(r,θ,z,t0) (1)
The assignment model of permeability in grid system:
K(r,θ,z,t0)=K0(z)+K0(z)γG(r,θ,z,t0) (2)
In formula,It is the porosity in initial time polar grid system on different location, wherein r is water
The wellbore radial direction of horizontal well, θ are the wellbore circumferential directions of horizontal well, and z is the wellbore axial direction of horizontal well, t0When being initial
It carves, unit min；It is the reservoir pore space on being parallel to horizontal well length direction on different location that well log interpretation obtains
Degree；K(r,θ,z,t0) it is permeability in initial time polar grid system on different location, unit mD；K0It (z) is to survey
Well explains the reservoir permeability on being parallel to horizontal well length direction on different location obtained, unit mD；G(r,θ,z,
t0) it is initial time random normal distribution array, range is -1~1；γ is standard deviation factor, and range is 0~1.
Further, the calculation of the individual well geological model after drilling damage described in step 2) is as follows:
Firstly, horizontal well drilling injury analog synthesis governing equation is established, such as following formula:
In formula, cpwFor the aerosol concentration in working solution；For reservoir porosity；T is time step, unit s；u
It is the flow velocity in the direction r, unit m/s；V is the flow velocity in the direction θ, unit m/s；ξ is the flow velocity in the direction z, unit m/s；DL
The diffusion coefficient of working solution, unit m2/s；σ is trapped particle amount, unit m；U is mean flow rate, unit m/s；α is wound
Evil intensive parameter, unit 1/m；
Then, primary condition and boundary condition are brought into formula (3) and solve the individual well geological model after obtaining drilling damage,
Wherein, the primary condition and boundary condition are as follows:
Primary condition:
Boundary condition:
Assuming that bit speed is constant, then the time that different well sections are impregnated by working solution are as follows:
In formula, tdIt is mud soak time, unit s；zdIt is net horizontal section position, wherein assuming horizontal well total length
For L, unit m；vdIt is bit speed, unit m/s；
t>td(zd) when pressure boundary condition:
P (r, θ, z)=P (r, θ+2 π, z) (8)
t>td(zd) when concentration boundary conditions:
cpw(r, θ, z)=cpw(r,θ+2π,z) (12)
In formula, μ is fluid dynamic viscosity, unit mPas；P is pressure, unit MPa；rwIt is wellbore radius, unit
For m；reIt is oil reservoir radius, unit m；PeIt is stratum reset pressure, unit MPa；PwWellbore pressure when being drilling well, unit are
MPa。
Further, the calculation of the horizontal well production profile in step 3) before drilling damage is as follows:
Firstly, establishing the horizontal well Reservoir Seepage continuity equation before drilling damage are as follows:
Then, the horizontal well production profile computation model before establishing drilling damage:
In formula, Q (z, t0) it is flow before drilling damage on horizontal well different location, unit m3/s；
Then, primary condition and boundary condition are brought into formula (14) and formula (15) solves the water before obtaining drilling damage
Horizontal well production profile, wherein the primary condition and boundary condition are as follows:
Primary condition:
P|r,θ,z=Pe (16)
Boundary condition:
In formula, PwfIt is the stream pressure in horizontal wellbore, unit MPa.
Further, the calculation of the horizontal well production profile in step 4) after drilling damage is as follows:
Firstly, establishing the horizontal well Reservoir Seepage continuity equation after drilling damage are as follows:
(42) the horizontal well production profile computation model after establishing drilling damage:
In formula, Q (z, t) is the flow before drilling damage on horizontal well different location, unit m3/s；
Then, primary condition and boundary condition are brought into formula (18) and (19) solves the horizontal well after obtaining drilling damage
Production profile, wherein the primary condition and boundary condition are as follows:
Primary condition:
P|r,θ,z=Pe (20)
Boundary condition:
Further, the calculation of the skin factor in step 5) on entire net horizontal section different location is as follows:
The skin factor of drilling well pollution:
In formula, Sd(z) it is skin factor on entire net horizontal section different location.
Further, the calculation of the cloth acid amount in step 6) on horizontal well different location is as follows:
In formula, VacidIt (z) is cloth acid amount on entire net horizontal section different location, unit m3；VtotalIt is individual well acidification
Total acid solution dosage, unit m3。
Compared with prior art, the present invention having the following advantages and benefits:
It is comprehensive that cloth acid method of the present invention and current existing well testing, dynamic analysis and well log interpretation can only obtain individual well
Skin factor is compared, and cloth acid method of the present invention sufficiently combines well log interpretation achievement, anti-using method for numerical simulation dynamic
Drilling fluid is drilled to the whole process of reservoir damage, the non-homogeneous skin factor on entire net horizontal section different location can be carried out
Fine quantitative evaluation；(uniform acid distribution, oval taper type cloth acid are divided based on permeability with existing level well cloth acid design method
Cloth rule cloth acid) it compares, this technology method can utmostly reduce horizontal well cloth using non-homogeneous skin factor as design considerations
The blindness of acid can effectively improve acid solution utilization rate, whole to improve horizontal well acidizing effect.
Detailed description of the invention
Attached drawing described herein is used to provide to further understand the embodiment of the present invention, constitutes one of the application
Point, do not constitute the restriction to the embodiment of the present invention.In the accompanying drawings:
Fig. 1 is Horizontal Well Log Interpretation porosity curve；
Fig. 2 is horizontal well individual well porosity geological model；
Fig. 3 is the permeability decrease amplitude profiles after drilling damage；
Fig. 4 is the horizontal well production profile before drilling damage；
Fig. 5 is the horizontal well production profile after drilling damage；
Fig. 6 is the skin factor profile on net horizontal section different location；
Fig. 7 is the cloth acid strength section on net horizontal section different location.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below with reference to embodiment and attached drawing, to this
Invention is described in further detail, and exemplary embodiment of the invention and its explanation for explaining only the invention, are not made
For limitation of the invention.
Embodiment 1:
As shown in Figures 1 to 7, a kind of cloth acid method based on horizontal well skin factor distribution characteristics, comprising the following steps:
1) the individual well geological model under primary condition, is established, the individual well geological model is the net established under polar coordinates
Reservoir permeability and porosity attribute that well log interpretation obtains are assigned to grid system by case system by the following method:
The assignment model of grid system porosity:
φ(r,θ,z,t0)=φ0(z)+φ0(z)γG(r,θ,z,t0) (1)
The assignment model of permeability in grid system:
K(r,θ,z,t0)=K0(z)+K0(z)γG(r,θ,z,t0) (2)
In formula,It is the porosity (zero dimension) in initial time polar grid system on different location,
In, r is the wellbore radial direction of horizontal well, and θ is the wellbore circumferential direction of horizontal well, and z is the wellbore axial direction of horizontal well, t0
It is initial time, unit min；Be well log interpretation obtain on being parallel to horizontal well length direction on different location
Reservoir porosity (zero dimension)；K(r,θ,z,t0) it is permeability in initial time polar grid system on different location, it is single
Position is mD；K0It (z) is the reservoir permeability on being parallel to horizontal well length direction on different location that well log interpretation obtains, it is single
Position is mD；G(r,θ,z,t0) it is initial time random normal distribution array, range is -1~1 (zero dimension)；γ is standard deviation
Coefficient, range are 0~1 (zero dimension)；
By taking Iraq's oil field a bite horizontal well as an example, the horizontal segment well log interpretation porosity curve of the well is as shown in Figure 1:
It brings porosity value into formula (1), obtains the individual well porosity geological model of the well, as shown in Fig. 2, using same
Quadrat method brings well log interpretation penetration value into formula (2), obtains the individual well permeability geological model of the well, no longer superfluous herein
It states.
2) the individual well geological model, based on step 1) building, the individual well geological model after calculating drilling damage, drilling damage
The calculation of individual well geological model afterwards is as follows:
Firstly, horizontal well drilling injury analog synthesis governing equation is established, such as following formula:
In formula, cpwFor the aerosol concentration (zero dimension) in working solution；For reservoir porosity (zero dimension)；When t is
Between step-length, unit s；U is the flow velocity in the direction r, unit m/s；V is the flow velocity in the direction θ, unit m/s；ξ is the stream in the direction z
Speed, unit m/s；DLThe diffusion coefficient of working solution, unit m2/s；σ is trapped particle amount, unit m；U is mean flow rate,
Unit is m/s；α is injury intensive parameter, unit 1/m；
Then, primary condition and boundary condition are brought into formula (3) and solve the individual well geological model after obtaining drilling damage,
Wherein, the primary condition and boundary condition are as follows:
Primary condition:
Boundary condition:
Assuming that bit speed is constant, then the time that different well sections are impregnated by working solution are as follows:
In formula, tdIt is mud soak time, unit s；zdIt is net horizontal section position, wherein assuming horizontal well total length
For L, unit m；vdIt is bit speed, unit m/s；
t>td(zd) when pressure boundary condition:
P (r, θ, z)=P (r, θ+2 π, z) (8)
t>td(zd) when concentration boundary conditions:
cpw(r, θ, z)=cpw(r,θ+2π,z) (12)
In formula, μ is fluid dynamic viscosity, unit mPas；P is pressure, unit MPa；rwIt is wellbore radius, unit
For m；reIt is oil reservoir radius, unit m；PeIt is stratum reset pressure, unit MPa；PwWellbore pressure when being drilling well, unit are
MPa。
Or by taking the well as an example, based on the well individual well geological model that step 1) obtains, drilling damage geological model is calculated,
Fig. 3 is the permeability decrease amplitude profiles after drilling damage:
3) the individual well geological model, based on step 1) building, the horizontal well production profile before calculating drilling damage, drilling well wound
The calculation of horizontal well production profile before evil is as follows:
Firstly, establishing the horizontal well Reservoir Seepage continuity equation before drilling damage are as follows:
Then, the horizontal well production profile computation model before establishing drilling damage:
In formula, Q (z, t0) it is flow before drilling damage on horizontal well different location, unit m3/s；
Then, primary condition and boundary condition are brought into formula (14) and formula (15) solves the water before obtaining drilling damage
Horizontal well production profile, wherein the primary condition and boundary condition are as follows:
Primary condition:
P|r,θ,z=Pe (16)
Boundary condition:
In formula, PwfIt is the stream pressure in horizontal wellbore, unit MPa；
Based on above formula, the well is calculated after drilling damage, and the horizontal well under same constant producing pressure differential 3MPa produces liquid
Section, as shown in Figure 5；
4), individual well geological model after the drilling damage being calculated based on step 2), the horizontal well after calculating drilling damage
The calculation of production profile, the horizontal well production profile after drilling damage is as follows:
Firstly, establishing the horizontal well Reservoir Seepage continuity equation after drilling damage are as follows:
(42) the horizontal well production profile computation model after establishing drilling damage:
In formula, Q (z, t) is the flow before drilling damage on horizontal well different location, unit m3/s；
Then, primary condition and boundary condition are brought into formula (18) and (19) solves the horizontal well after obtaining drilling damage
Production profile, wherein the primary condition and boundary condition are as follows:
Primary condition:
P|r,θ,z=Pe (20)
Boundary condition:
Based on above formula, the well is calculated after drilling damage, and the horizontal well under same constant producing pressure differential 3MPa produces liquid
Section, as shown in Figure 5；
5), after the drilling damage that the horizontal well production profile before the drilling damage obtained according to step 3) and step 4) obtain
Horizontal well production profile, skin factor on entire net horizontal section different location is calculated, along entire net horizontal section difference position
The calculation for the skin factor set is as follows:
The skin factor of drilling well pollution:
In formula, Sd(z) it is skin factor (zero dimension) on entire net horizontal section different location；
In conjunction with the forward and backward horizontal well production profile of drilling damage that step 3) and step 4) obtain, calculated using formula (22)
The skin factor distribution obtained on entire net horizontal section different location is as shown in Figure 6；
6), the skin factor obtained according to step 5) calculates on horizontal well different location along the distribution characteristics of net horizontal section
Cloth acid amount, the calculation of the cloth acid amount on horizontal well different location is as follows:
In formula, VacidIt (z) is cloth acid amount on entire net horizontal section different location, unit m3；VtotalIt is individual well acidification
Total acid solution dosage, unit m3。
Divide 10 sections of average out to for horizontal well, cloth acid strength such as Fig. 7 of every section of net horizontal section is calculated according to formula (23)
It is shown.
In the present embodiment, cloth acid method of the present invention is based on individual well geological model, and simulation calculates practical drilling speed and brill
Filter loss of the drilling fluid of pressure at different level well section position, the reservoir after drilling damage is determined by drilling fluid filtration
Permeability.It brings the forward and backward permeability of drilling damage into darcy flow equation simulation and calculates horizontal well under constant producing pressure differential
Production profile, so that it is determined that taking the serious position of injury, emphasis along the non-homogeneous skin factor distribution situation of entire net horizontal section
The principle of cloth acid rationally designs the cloth acid amount on net horizontal section different location.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects
It is described in detail, it should be understood that being not intended to limit the present invention the foregoing is merely a specific embodiment of the invention
Protection scope, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should all include
Within protection scope of the present invention.
Claims (7)
1. a kind of cloth acid method based on horizontal well skin factor distribution characteristics, which comprises the following steps:
1) the individual well geological model under primary condition, is established；
2) the individual well geological model, based on step 1) building, the individual well geological model after calculating drilling damage；
3) the individual well geological model, based on step 1) building, the horizontal well production profile before calculating drilling damage；
4), individual well geological model after the drilling damage being calculated based on step 2), the horizontal well after calculating drilling damage produce liquid
Section；
5), the water after the drilling damage that the horizontal well production profile before the drilling damage obtained according to step 3) and step 4) obtain
Horizontal well production profile calculates the skin factor on entire net horizontal section different location；
6), the skin factor obtained according to step 5) calculates the cloth on horizontal well different location along the distribution characteristics of net horizontal section
Acid amount.
2. a kind of cloth acid method based on horizontal well skin factor distribution characteristics according to claim 1, which is characterized in that
Individual well geological model described in step 1) is the grid system established under polar coordinates, the reservoir permeability that well log interpretation is obtained
And porosity attribute is assigned to grid system by the following method:
The assignment model of grid system porosity:
φ(r,θ,z,t0)=φ0(z)+φ0(z)γG(r,θ,z,t0) (1)
The assignment model of permeability in grid system:
K(r,θ,z,t0)=K0(z)+K0(z)γG(r,θ,z,t0) (2)
In formula,It is the porosity in initial time polar grid system on different location, wherein r is horizontal well
Wellbore radial direction, θ are the wellbore circumferential directions of horizontal well, and z is the wellbore axial direction of horizontal well, t0It is initial time, it is single
Position is min；It is the reservoir porosity on being parallel to horizontal well length direction on different location that well log interpretation obtains；K
(r,θ,z,t0) it is permeability in initial time polar grid system on different location, unit mD；K0It (z) is well logging solution
Release the reservoir permeability on being parallel to horizontal well length direction on different location of acquisition, unit mD；G(r,θ,z,t0) be
Initial time random normal is distributed array, and range is -1~1；γ is standard deviation factor, and range is 0~1.
3. a kind of cloth acid method based on horizontal well skin factor distribution characteristics according to claim 1, which is characterized in that
The calculation of individual well geological model after drilling damage described in step 2) is as follows:
Firstly, horizontal well drilling injury analog synthesis governing equation is established, such as following formula:
In formula, cpwFor the aerosol concentration in working solution；For reservoir porosity；T is time step, unit s；U is the side r
To flow velocity, unit m/s；V is the flow velocity in the direction θ, unit m/s；ξ is the flow velocity in the direction z, unit m/s；DLWorking solution
Diffusion coefficient, unit m2/s；σ is trapped particle amount, unit m；U is mean flow rate, unit m/s；α is injury intensity
Parameter, unit 1/m；
Then, primary condition and boundary condition are brought into formula (3) and solve the individual well geological model after obtaining drilling damage,
In, the primary condition and boundary condition are as follows:
Primary condition:
Boundary condition:
Assuming that bit speed is constant, then the time that different well sections are impregnated by working solution are as follows:
In formula, tdIt is mud soak time, unit s；zdIt is net horizontal section position, wherein assume that horizontal well total length is L,
Unit is m；vdIt is bit speed, unit m/s；
t>td(zd) when pressure boundary condition:
P (r, θ, z)=P (r, θ+2 π, z) (8)
t>td(zd) when concentration boundary conditions:
cpw(r, θ, z)=cpw(r,θ+2π,z) (12)
In formula, μ is fluid dynamic viscosity, unit mPas；P is pressure, unit MPa；rwIt is wellbore radius, unit m；
reIt is oil reservoir radius, unit m；PeIt is stratum reset pressure, unit MPa；PwWellbore pressure when being drilling well, unit MPa.
4. a kind of cloth acid method based on horizontal well skin factor distribution characteristics according to claim 1, which is characterized in that
The calculation of horizontal well production profile in step 3) before drilling damage is as follows:
Firstly, establishing the horizontal well Reservoir Seepage continuity equation before drilling damage are as follows:
Then, the horizontal well production profile computation model before establishing drilling damage:
In formula, Q (z, t0) it is flow before drilling damage on horizontal well different location, unit m3/s；
Then, primary condition and boundary condition are brought into formula (14) and formula (15) solves the horizontal well before obtaining drilling damage
Production profile, wherein the primary condition and boundary condition are as follows:
Primary condition:
P|r,θ,z=Pe (16)
Boundary condition:
In formula, PwfIt is the stream pressure in horizontal wellbore, unit MPa.
5. a kind of cloth acid method based on horizontal well skin factor distribution characteristics according to claim 1, which is characterized in that
The calculation of horizontal well production profile in step 4) after drilling damage is as follows:
Firstly, establishing the horizontal well Reservoir Seepage continuity equation after drilling damage are as follows:
(42) the horizontal well production profile computation model after establishing drilling damage:
In formula, Q (z, t) is the flow before drilling damage on horizontal well different location, unit m3/s；
Then, primary condition and boundary condition are brought into formula (18) and (19) solves the horizontal well after obtaining drilling damage and produce liquid
Section, wherein the primary condition and boundary condition are as follows:
Primary condition:
P|r,θ,z=Pe (20)
Boundary condition:
6. a kind of cloth acid method based on horizontal well skin factor distribution characteristics according to claim 1, which is characterized in that
The calculation of skin factor in step 5) on entire net horizontal section different location is as follows:
The skin factor of drilling well pollution:
In formula, Sd(z) it is skin factor on entire net horizontal section different location.
7. a kind of cloth acid method based on horizontal well skin factor distribution characteristics according to claim 1, which is characterized in that
The calculation of cloth acid amount in step 6) on horizontal well different location is as follows:
In formula, VacidIt (z) is cloth acid amount on entire net horizontal section different location, unit m3；VtotalIt is that individual well is acidified total acid
Liquid dosage, unit m3。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811026178.8A CN109138949A (en) | 2018-09-28 | 2018-09-28 | A kind of cloth acid method based on horizontal well skin factor distribution characteristics |
| PCT/CN2019/107628 WO2020063604A1 (en) | 2018-09-28 | 2019-09-25 | Method for distributing acid on basis of distribution characteristics of skin coefficient of horizontal well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811026178.8A CN109138949A (en) | 2018-09-28 | 2018-09-28 | A kind of cloth acid method based on horizontal well skin factor distribution characteristics |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109138949A true CN109138949A (en) | 2019-01-04 |
Family
ID=64826726
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811026178.8A Pending CN109138949A (en) | 2018-09-28 | 2018-09-28 | A kind of cloth acid method based on horizontal well skin factor distribution characteristics |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109138949A (en) |
| WO (1) | WO2020063604A1 (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109854215A (en) * | 2019-01-28 | 2019-06-07 | 西南石油大学 | A method of acidification water injection well sophisticated category |
| CN110454152A (en) * | 2019-08-12 | 2019-11-15 | 成都北方石油勘探开发技术有限公司 | The non-homogeneous injury skin factor calculation method of fracture-type reservoir |
| CN110608036A (en) * | 2019-07-24 | 2019-12-24 | 王新海 | Method for calculating total skin coefficient of multilayer oil reservoir |
| CN110826250A (en) * | 2019-11-19 | 2020-02-21 | 成都北方石油勘探开发技术有限公司 | Drilling fluid damage calculation method for horizontal well |
| CN110910966A (en) * | 2019-11-19 | 2020-03-24 | 成都北方石油勘探开发技术有限公司 | Constant-displacement acid-squeezing calculation method for three-dimensional porous continuous medium discrete model |
| WO2020063604A1 (en) * | 2018-09-28 | 2020-04-02 | 成都北方石油勘探开发技术有限公司 | Method for distributing acid on basis of distribution characteristics of skin coefficient of horizontal well |
| CN112377149A (en) * | 2020-12-09 | 2021-02-19 | 中国石油天然气集团有限公司 | Method for determining acid strength for washing horizontal water injection well |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101298835A (en) * | 2007-04-30 | 2008-11-05 | 陈菲 | Controllable acid distributing tube and controllable acid distributing method |
| US20090166024A1 (en) * | 2004-06-18 | 2009-07-02 | Schlumberger Technology Corporation | Completion apparatus for measuring streaming potentials and determining earth formation characteristics |
| CN104314547A (en) * | 2014-08-18 | 2015-01-28 | 北京运嘉普泰石油天然气技术发展有限公司 | HAFD segmented optimization method and system |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108194053B (en) * | 2017-12-07 | 2019-12-06 | 中国石油化工股份有限公司华北油气分公司石油工程技术研究院 | method and device for determining matrix acidizing acid liquid dosage of open-hole horizontal well |
| CN109138949A (en) * | 2018-09-28 | 2019-01-04 | 成都北方石油勘探开发技术有限公司 | A kind of cloth acid method based on horizontal well skin factor distribution characteristics |
-
2018
- 2018-09-28 CN CN201811026178.8A patent/CN109138949A/en active Pending
-
2019
- 2019-09-25 WO PCT/CN2019/107628 patent/WO2020063604A1/en active Application Filing
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090166024A1 (en) * | 2004-06-18 | 2009-07-02 | Schlumberger Technology Corporation | Completion apparatus for measuring streaming potentials and determining earth formation characteristics |
| CN101298835A (en) * | 2007-04-30 | 2008-11-05 | 陈菲 | Controllable acid distributing tube and controllable acid distributing method |
| CN104314547A (en) * | 2014-08-18 | 2015-01-28 | 北京运嘉普泰石油天然气技术发展有限公司 | HAFD segmented optimization method and system |
Non-Patent Citations (3)
| Title |
|---|
| 付宣等: "非均质碳酸盐岩油藏水平井分段", 《油气地质与采收率》 * |
| 柳明等: "不同注入条件下地层酸蚀蚓孔扩展规律", 《油气地质与采收率》 * |
| 薛衡等: "碳酸盐岩水平井酸化模式及工艺参数优化研究", 《油气藏评价与开发》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020063604A1 (en) * | 2018-09-28 | 2020-04-02 | 成都北方石油勘探开发技术有限公司 | Method for distributing acid on basis of distribution characteristics of skin coefficient of horizontal well |
| CN109854215A (en) * | 2019-01-28 | 2019-06-07 | 西南石油大学 | A method of acidification water injection well sophisticated category |
| CN110608036A (en) * | 2019-07-24 | 2019-12-24 | 王新海 | Method for calculating total skin coefficient of multilayer oil reservoir |
| 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 |
| CN110826250A (en) * | 2019-11-19 | 2020-02-21 | 成都北方石油勘探开发技术有限公司 | Drilling fluid damage calculation method for horizontal well |
| CN110910966A (en) * | 2019-11-19 | 2020-03-24 | 成都北方石油勘探开发技术有限公司 | Constant-displacement acid-squeezing calculation method for three-dimensional porous continuous medium discrete model |
| CN110910966B (en) * | 2019-11-19 | 2022-04-22 | 成都北方石油勘探开发技术有限公司 | Constant-displacement acid-squeezing calculation method for three-dimensional porous continuous medium discrete model |
| CN112377149A (en) * | 2020-12-09 | 2021-02-19 | 中国石油天然气集团有限公司 | Method for determining acid strength for washing horizontal water injection well |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020063604A1 (en) | 2020-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109138949A (en) | A kind of cloth acid method based on horizontal well skin factor distribution characteristics | |
| CN109209333B (en) | Shale gas multi-well group efficient mining interval optimization method | |
| CN112343587A (en) | Ultra-low permeability reservoir dominant seepage channel identification and characterization method | |
| CN109342093A (en) | A kind of underground water-seal oil depot country rock dominant flowing path discrimination method | |
| CN104834807A (en) | Stress sensitive reservoir stratum relative permeability calculation method based on fractal theory | |
| CN107526891A (en) | A kind of polymer flooding macropore oil reservoir well test analysis method | |
| CN109577929A (en) | Quantitative evaluation method for establishing effective displacement of ultra-low permeability tight reservoir horizontal well | |
| CN109958404A (en) | Bottom water reservoir bottom water coning analogy method and water plugging profile control agent calculation method of parameters | |
| CN104533519B (en) | Management method for water burst water disasters in event of vertical shaft wellbore passing through strong water-containing thick rock layer | |
| CN103643928A (en) | Step-by-step deep profile control method based on pressure field and flow velocity field distribution | |
| CN110984974B (en) | Method for determining reasonable gas production speed of water-gas reservoir based on water invasion speed | |
| CN112377182B (en) | Method and device for determining parameters of large-scale hole crack type carbonate reservoir | |
| CN110991933A (en) | Mountain shale gas resource evaluation method and system | |
| CN114186440A (en) | Geological-engineering double-track shale compressibility comprehensive evaluation method | |
| CN104594854A (en) | Bottom water reservoir horizontal well variable-parameter perforation water-controlled completion method and device | |
| CN106022626A (en) | Water flooding reservoir positioning plugging water flooding virtual front edge calculation method | |
| CN110390154A (en) | A method of improving Complex reservoir reservoir numerical simulation efficiency | |
| CN106481315A (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
| CN116291406A (en) | Method and system for identifying large pore canal of offshore loose sandstone reservoir | |
| CN116629160A (en) | Dynamic prediction method for water inflow of mine roof | |
| CN114252381B (en) | Method for evaluating pollution degree of fractured reservoir horizontal well drilling fluid after solid-phase pollution | |
| CN112685970B (en) | Quantitative characterization method and system for seepage interface of flow unit of water-drive reservoir | |
| CN110714742A (en) | Method for improving recovery ratio of bottom water condensate gas reservoir | |
| CN112487594B (en) | Oil reservoir water body multiple calculation method and device | |
| CN113832999A (en) | Intelligent control system and control method for sand-gravel stratum deep foundation pit dewatering |
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: 20190104 |