CN107169248B - Special stratum safe mud density window determination method - Google Patents

Special stratum safe mud density window determination method Download PDF

Info

Publication number
CN107169248B
CN107169248B CN201710542811.8A CN201710542811A CN107169248B CN 107169248 B CN107169248 B CN 107169248B CN 201710542811 A CN201710542811 A CN 201710542811A CN 107169248 B CN107169248 B CN 107169248B
Authority
CN
China
Prior art keywords
stress
well
pressure
well wall
mud
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.)
Active
Application number
CN201710542811.8A
Other languages
Chinese (zh)
Other versions
CN107169248A (en
Inventor
袁俊亮
周建良
罗洪斌
许亮斌
殷志明
郝希宁
张玉亭
高飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
Beijing Research Center of CNOOC China Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
Beijing Research Center of CNOOC China Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, Beijing Research Center of CNOOC China Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201710542811.8A priority Critical patent/CN107169248B/en
Publication of CN107169248A publication Critical patent/CN107169248A/en
Application granted granted Critical
Publication of CN107169248B publication Critical patent/CN107169248B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Structural Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Civil Engineering (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention relates to a method for determining a safe mud density window of a special stratum, which comprises the following steps: 1) determining the mechanical parameters of sedimentary rock stratum rocks according to the acoustic time difference, density and natural gamma data obtained by logging of a target well and an indoor rock mechanical experiment; 2) determining geomechanical parameters according to natural gamma logging information and a measured value of the stratum pressure-bearing capacity of a drilling site; 3) calculating the stress state around the well by using geomechanical parameters; 4) judging whether the stratum has an upper collapse pressure limit; 5) according to the principle of combination of numbers and shapes, aiming at three deformation damage forms, the lower collapse pressure limit, the fracture pressure and the upper collapse pressure limit are determined, and a safe mud density window of a special stratum is obtained. By the method, the special stratum with the upper collapse pressure limit can be effectively and accurately identified, and the safe mud density window of the special stratum is determined.

Description

Special stratum safe mud density window determination method
Technical Field
The invention relates to the technical field of petroleum drilling, in particular to a method for determining a special stratum safe mud density window.
Background
In the field of petroleum industry, the problem of borehole wall instability in the process of oil and gas field exploration and development severely restricts the drilling speed and causes huge economic loss. According to the drilling industry, the economic loss caused by borehole wall instability worldwide is over $ 1 billion per year. During actual drilling, formations in which borehole wall instability problems typically occur include: hard brittle shale formations, bedding formations, fractured formations, highly steep complex formations, and high temperature high pressure formations. There are two main forms of damage to borehole wall instability: borehole wall collapse and borehole wall rupture. The former is mainly caused by that the mud density is too low, so that the structural compression stress borne by the well wall exceeds the compressive strength of rock, and the well wall is subjected to shear failure, and the mud density at the moment is called as collapse pressure Pt|lower(ii) a The latter is mainly caused by the fact that the mud density is too high, so that the tensile stress brought by the liquid column born by the well wall exceeds the tensile strength of the stratum, and the mud density at this moment is called as the fracture pressure Pf
The method for analyzing the stability of the well wall generally comprises the step of calculating the density of slurry for maintaining the stability of the well wall by simulating the stress state bonding strength criterion around the well, wherein the minimum density of the slurry required for ensuring that the well wall is not collapsed and damaged is the collapse pressure. It is therefore common to maintain the borehole wall stable during drilling with an appropriate mud density and it is believed that the borehole wall will not experience shear nor tension damage as long as the mud density is greater than the collapse pressure and does not exceed the fracture pressure. However, the actual drilling phenomenon of the part of the formation ahead of the Tarim mountain with the field reflects: there is a special formation and there is an upper limit P of the collapse pressure during drillingt|upOnce the mud density is too high, the upper limit of the collapse pressure is exceededPt|upSevere borehole wall collapse can also occur at this time. The safe mud density window for this particular type of formation should not be limited by the fracture pressure curve alone. There is no current method to identify such specific formations and calculate a safe mud density window for such specific formations.
Disclosure of Invention
In view of the above problems, the present invention provides a method for determining a safe mud density window of a special formation, by which the special formation with an upper collapse pressure limit can be effectively and accurately identified, and the safe mud density window of the special formation can be determined.
In order to achieve the purpose, the invention adopts the following technical scheme: a method for determining a safe mud density window of a special stratum is characterized by comprising the following steps:
1) determining the mechanical parameters of sedimentary rock stratum rocks according to the acoustic time difference, density and natural gamma data obtained by logging of a target well and an indoor rock mechanical experiment;
2) determining geomechanical parameters according to the natural gamma logging data obtained in the step 1) and the measured value of the stratum pressure-bearing capacity of the drilling site;
3) calculating the stress state around the well by using the geomechanical parameters obtained in the step 2);
4) determining three deformation failure modes generated in the process of bearing different mud liquid column pressures by the well wall according to Mohr-Coulomb yield criterion and combining the well peripheral stress state and geomechanical parameters, and judging whether the stratum has an upper limit of collapse pressure;
5) determining the lower limit P of the collapse pressure aiming at three deformation damage forms according to the principle of combination of numbers and shapest|lowerRupture pressure PfAnd upper limit of collapse pressure Pt|upAnd obtaining a safe mud density window of the special stratum as follows:
Pt|lower<ρmud<min(Pt|up,Pf) (1)
in the above formula, ρmudIs the mud density.
In the step 1), the national standard engineering rock mass test method standard (GB/T50266-99) is adopted in the indoor rock mechanics experiment for determining the sedimentary rock formation rock mechanics parameters; the sedimentary rock stratum lithology is sandstone-mudstone, the mechanical property is an ideal linear elastomer with homogeneous isotropy, and a borehole is drilled to meet the Mohr-Coulomb yield criterion:
τ=C+σ·tgφ (2)
wherein tau is the shear stress of sedimentary rock; c is cohesive force of sedimentary rock; phi is the internal friction angle of the sedimentary rock; σ is the positive stress of the sedimentary rock.
In step 2) above, the geomechanical parameter comprises overburden pressure σVHorizontal maximum stress σHAnd horizontal minimum ground stress σh
Figure BDA0001342262500000021
Where rho is the density of sedimentary rock, α is the effective stress coefficient, omega12To construct the stress coefficient; 0-H is a depth integration interval; h is the well depth of the target well; g is the acceleration of gravity; ppIs the pore pressure; v is the Poisson's ratio.
In the step 3), according to the effective stress theory of the porous medium, the thick-wall cylinder stressed unevenly has the well-periphery stress state expression as follows:
Figure BDA0001342262500000022
in the formula, theta is an included angle between a well circumferential radial vector and a horizontal maximum principal stress direction; sigma'ZAxial effective stress of the well wall; sigma'θIs the tangential effective stress of the well wall; sigma'rRadial effective stress of the well wall; piThe pressure of the mud-liquid column, R the well-to-well distance, R the radius of the well bore, f the porosity, η the nonlinear correction coefficient, 0.95, delta the seepage coefficient, ξ the intermediate variable,
Figure BDA0001342262500000031
in the step 4), for the hard and brittle shale stratum, the well wall has the following three forms of deformation damage:
(1) when the density of the slurry is lower than the conventional lower limit P of the collapse pressuret|lowerWhen the most dangerous azimuth around the well is the azimuth theta of 90 degrees or 270 degrees, the difference value between the maximum principal stress and the minimum principal stress around the well is sigma'θ-σ′rThe well wall will be damaged by lower limit shearing;
(2) when the mud density is higher than the fracture pressure PfIn this case, the most dangerous azimuth around the well is the azimuth of θ ═ 0 ° or 180 °, and the stress state around the well is σ'θIf the well wall is less than 0, the well wall will be subjected to tensile damage;
(3) for a particular type of formation, at mud column pressure PiReaches a rupture pressure PfPreviously, the difference between the maximum principal stress and the minimum principal stress of the well wall at the positions of 0 DEG and 180 DEG is sigma DEG'Z-σ′θIf the difference is too high, the upper limit shearing damage will occur to the well wall, and the pressure P of the mud-liquid column at the momentiIs the upper limit P of the collapse pressuret|up
In the step 5), the lower limit of collapse pressure P in the three deformation damage forms existing on the well wallt|lowerRupture pressure PfUpper limit of collapse pressure Pt|upThe relative size relationship of the three is as follows:
(1) the tangential to radial stress difference curve for the borehole wall represents the maximum principal stress σ 'on the borehole wall at the orientation of θ -90 ° or 270 ° as the mud column pressure increases'θAnd minimum principal stress σ'rThe difference between the two values and the equation expression of the tangential and radial stress difference curve of the well wall are as follows:
y=3σHh-2x (5)
in the formula, x is the pressure of a mud slurry column; y is the stress difference;
the abscissa of the intersection point of the tangential and radial stress difference curve of the well wall and the stress difference curve required by shear failure is collapse pressureLower limit of force Pt|lowerThe calculation expression is as follows:
Figure BDA0001342262500000032
wherein UCS is the uniaxial compressive strength of sedimentary rock, and UCS is 2. C.K; k is an intermediate variable, K ═ tan (pi/4 + phi/2);
(2) the axial and tangential stress difference curves of the borehole wall represent the maximum principal stress σ 'on the borehole wall at an orientation of 0 ° or 180 ° with increasing mud column pressure'ZAnd minimum principal stress σ'θThe difference between the two, axial and tangential stress difference curve equation expression of the well wall is as follows:
y=σV-(2v-1)σH+(2v-3)σh+x (7)
the abscissa of the intersection point of the axial and tangential stress difference curves of the well wall and the stress difference curve required by shear failure is the upper limit P of collapse pressuret|upSubstituting the formula (4) into the Mohr-Coulomb yield criterion, and finishing to obtain the upper limit P of the collapse pressuret|upThe calculation expression of (a) is as follows:
Figure BDA0001342262500000041
the expression of the transition point (Xo, Yo) between the lower limit shear failure and the upper limit shear failure is as follows:
Figure BDA0001342262500000042
(3) the tangential stress curve of the borehole wall represents the tangential effective stress σ 'at an orientation of 0 or 180 as the mud slurry column pressure increases'θWhen the pressure is reduced to be less than 0, the well wall is subjected to tensile failure, and the abscissa of the intersection point of the tangential stress curve of the well wall and the x axis is fracture pressure PfThe calculation expression is as follows:
Pf=3σhH-αPp+St(10)
wherein St is the tensile strength of the well wall rock.
Less than critical uniaxial compressive strength UCS for uniaxial compressive strength UCSpThe formation of (1), which is a special formation that must consider the upper collapse pressure limit, the critical uniaxial compressive strength UCSpCalculated by equation (11):
UCSP=σV-2v(σHh)-αPp(11)。
due to the adoption of the technical scheme, the invention has the following advantages: 1. the invention provides a concept and a calculation method of an upper collapse pressure limit, and considers that a safety density window is limited by a lower collapse pressure limit and a fracture pressure, and the upper collapse pressure limit also exists for a special stratum, and a real safety density window is determined by the three. Compared with the prior art, the invention expands the knowledge of the safe density window. 2. The invention provides a method for judging a special stratum with a collapse pressure upper limit, which can calculate the critical Uniaxial Compressive Strength (UCS) through overburden pressure, horizontal ground stress and stratum pore pressurep) The stratum with the uniaxial strength lower than the critical value belongs to a special stratum with the collapse pressure upper limit needing to be considered, so that the invention also provides the applicable conditions while obtaining the determination method.
Drawings
FIG. 1 is a schematic diagram of a borehole wall shear failure at a lower collapse pressure limit;
FIG. 2 is a schematic view of a borehole wall in tension failure;
FIG. 3 is a schematic diagram of borehole wall shear failure at an upper collapse pressure limit;
FIG. 4 is a graph showing the pressure profile of a slurry column in which lower limit shear failure, upper limit shear failure and tension failure occur, respectively;
FIG. 5 is a flow chart of the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and examples. It is to be understood, however, that the drawings are provided solely for the purposes of promoting an understanding of the invention and that they are not to be construed as limiting the invention.
As shown in fig. 5, the method for determining the safe mud density window of the special formation provided by the invention comprises the following steps:
1) determining mechanical parameters of sedimentary rock stratum rocks including uniaxial compressive strength, cohesive force and internal friction angle of the sedimentary rock stratum rocks according to acoustic time difference, density and natural gamma data obtained by logging of a target well and an indoor rock mechanical experiment;
in the embodiment, the indoor rock mechanics experiment for determining the sedimentary rock formation rock mechanics parameters adopts the national standard engineering rock mass test method standard (GB/T50266-99). Wherein, the sedimentary rock stratum lithology is sand shale, the mechanical property is an ideal line elastomer with homogeneous isotropy, and the borehole is drilled to meet the Mohr-Coulomb (Mohr-Coulomb) yield criterion:
τ=C+σ·tgφ (1)
wherein tau is the shear stress of sedimentary rock; c is cohesive force of sedimentary rock; phi is the internal friction angle of the sedimentary rock; σ is the positive stress of the sedimentary rock.
2) Determining geomechanical parameters according to the natural gamma logging data obtained in the step 1) and the measured value of the stratum pressure-bearing capacity of the drilling site; wherein the geomechanical parameter comprises overburden pressure sigmaVHorizontal maximum stress σHAnd horizontal minimum ground stress σh
Figure BDA0001342262500000051
Where rho is the density of sedimentary rock, α is the effective stress coefficient, PpIs the pore pressure; omega12The formation stress coefficient is determined by a drilling site stratum bearing capacity test; 0-H is a depth integration interval; h is the well depth of the target well; g is the acceleration of gravity; v is the Poisson's ratio.
3) Calculating the stress state (sigma ') around the well by using the geomechanical parameters obtained in the step 2)'Z、σ′r、σ′θ):
According to the effective stress theory of porous media, the thick-wall cylinder under the action of non-uniform ground stress has the following expression of the stress state around the well:
Figure BDA0001342262500000052
in the formula, theta is an included angle between a well circumferential radial vector and a horizontal maximum principal stress direction; sigma'ZAxial effective stress of the well wall; sigma'θIs the tangential effective stress of the well wall; sigma'rRadial effective stress of the well wall; piThe pressure of the mud-liquid column, R the well-to-well distance, R the radius of the well bore, f the porosity, η the nonlinear correction coefficient, 0.95, delta the seepage coefficient, ξ the intermediate variable,
Figure BDA0001342262500000053
4) according to Mohr-Coulomb yield criterion, three deformation failure modes generated in the process of bearing different mud liquid column pressures on the well wall are determined by combining the stress state around the well and geomechanical parameters:
for a hard and brittle shale stratum, the seepage effect and the nonlinear characteristics can be simplified, and the well wall has the following three deformation damages:
(1) when the density of the slurry is lower than the conventional lower limit P of the collapse pressuret|lowerWhen the most dangerous azimuth around the well is the azimuth theta of 90 degrees or 270 degrees, the difference value between the maximum principal stress and the minimum principal stress around the well is sigma'θ-σ′rThe well wall will have lower limit shear failure, the failure form is shown in figure 1;
(2) when the mud density is higher than the fracture pressure PfIn this case, the most dangerous azimuth around the well is the azimuth of θ ═ 0 ° or 180 °, and the stress state around the well is σ'θIf the pressure is less than 0, the well wall will be subjected to tensile damage, and the damage form is shown in figure 2;
(3) for a particular type of formation, at mud column pressure PiReaches a rupture pressure PfPreviously, the difference between the maximum principal stress and the minimum principal stress of the well wall at the positions of 0 DEG and 180 DEG is sigma DEG'Z-σ′θIf the difference is too high, the upper limit shearing damage will occur to the well wall,the failure is shown in FIG. 3, where the mud column pressure P isiIs the upper limit P of the collapse pressuret|up
5) According to the principle of combination of numbers and shapes, the lower limit P of collapse pressure is given for three deformation failure formst|lowerRupture pressure PfUpper limit of collapse pressure Pt|upAnd giving a safe mud density window:
as shown in FIG. 4, the lower limit of collapse pressure P in the three types of deformation damage existing on the well wallt|lowerRupture pressure PfUpper limit of collapse pressure Pt|upThe relative size relationship of the three is as follows:
(1) the tangential to radial stress difference curve for the borehole wall represents the maximum principal stress σ 'on the borehole wall at the orientation of θ -90 ° or 270 ° as the mud column pressure increases'θAnd minimum principal stress σ'rThe difference between the two values and the equation expression of the tangential and radial stress difference curve of the well wall are as follows:
y=3σHh-2x (4)
in the formula, x is the pressure of a mud slurry column; and y is the stress difference.
The abscissa of the intersection point of the tangential and radial stress difference curve of the well wall and the stress difference curve required by shear failure is the lower limit P of collapse pressuret|lowerThe calculation expression is as follows:
Figure BDA0001342262500000061
wherein UCS is the uniaxial compressive strength of sedimentary rock, and UCS is 2. C.K; k is an intermediate variable, K ═ tan (pi/4 + phi/2);
(2) the axial and tangential stress difference curves of the borehole wall represent the maximum principal stress σ 'on the borehole wall at an orientation of 0 ° or 180 ° with increasing mud column pressure'ZAnd minimum principal stress σ'θThe difference between the two, axial and tangential stress difference curve equation expression of the well wall is as follows:
y=σV-(2v-1)σH+(2v-3)σh+x (6)
the abscissa of the intersection point of the axial and tangential stress difference curves of the well wall and the stress difference curve required by shear failure is the upper limit P of collapse pressuret|upSubstituting the formula (3) into the Mohr-Coulomb yield criterion, and finishing to obtain the upper limit P of the collapse pressuret|upThe calculation expression of (a) is as follows:
Figure BDA0001342262500000071
the expression of the transition point (Xo, Yo) between the lower limit shear failure and the upper limit shear failure is as follows:
Figure BDA0001342262500000072
(3) the tangential stress curve of the borehole wall represents the tangential effective stress σ 'at an orientation of 0 or 180 as the mud slurry column pressure increases'θWhen the pressure is reduced to be less than 0, the well wall is subjected to tensile failure, and the abscissa of the intersection point of the tangential stress curve of the well wall and the x axis is fracture pressure PfThe calculation expression is as follows:
Pf=3σhH-αPp+St(9)
wherein St is the tensile strength of the well wall rock.
If the upper limit of the collapse pressure Pt|upGreater than the rupture pressure PfThen there is no upper collapse pressure limit for that type of formation (because the formation will first experience tensile failure); if the upper limit of the collapse pressure Pt|upLess than the cracking pressure PfAnd judging that the stratum has the upper limit of collapse pressure, namely judging the stratum to be the special stratum.
Comparing the collapse pressure upper limit with the rupture pressure, and taking the smaller value of the two, i.e. min (P)t|up,Pf) As an upper limit on the safe mud density window for a particular formation. In summary, the safe mud density window for a particular formation is as follows:
Pt|lower<ρmud<min(Pt|up,Pf) (10)
in the formula, ρmudIs the mud density.
6) To clarify the application range of the present invention, the critical uniaxial compressive strength UCS can be calculated by analyzing the mechanical properties of the special stratum, i.e. making formula (7) and formula (9) equalp
UCSP=σV-2v(σHh)-αPp(11)
Less than critical uniaxial compressive strength UCS for uniaxial compressive strength UCSpThe formation of (2) is a specific formation for which the upper collapse pressure limit must be considered.
The above embodiments are only used for illustrating the present invention, and the structure, connection mode, manufacturing process, etc. of the components may be changed, and all equivalent changes and modifications performed on the basis of the technical solution of the present invention should not be excluded from the protection scope of the present invention.

Claims (2)

1. A method for determining a safe mud density window of a special stratum is characterized by comprising the following steps:
1) determining the mechanical parameters of sedimentary rock stratum rocks according to the acoustic time difference, density and natural gamma data obtained by logging of a target well and an indoor rock mechanical experiment; wherein, the indoor rock mechanics experiment for determining the sedimentary rock stratum rock mechanics parameters adopts the national standard of engineering rock mass test method standard GB/T50266-99; the sedimentary rock stratum lithology is sandstone-mudstone, the mechanical property is an ideal linear elastomer with homogeneous isotropy, and a borehole is drilled to meet the Mohr-Coulomb yield criterion:
τ=C+σ·tgφ (1)
wherein tau is the shear stress of sedimentary rock; c is cohesive force of sedimentary rock; phi is the internal friction angle of the sedimentary rock; σ is the normal stress of the sedimentary rock;
2) determining geomechanical parameters according to the natural gamma logging data obtained in the step 1) and the measured value of the stratum pressure-bearing capacity of the drilling site; wherein the geomechanical parameter comprises overburden pressure sigmaVHorizontal maximum stress σHAnd horizontal minimum ground stress σh
Figure FDA0002391239950000011
Where rho is the density of sedimentary rock, α is the effective stress coefficient, omega12To construct the stress coefficient; 0-H is a depth integration interval; h is the well depth of the target well; g is the acceleration of gravity; ppIs the pore pressure;vis the poisson ratio;
3) calculating the stress state around the well by using the geomechanical parameters obtained in the step 2), wherein the specific process is as follows:
according to the effective stress theory of porous media, the thick-wall cylinder under the action of non-uniform ground stress has the following expression of the stress state around the well:
Figure FDA0002391239950000012
in the formula, theta is an included angle between a well circumferential radial vector and a horizontal maximum principal stress direction; sigma'ZAxial effective stress of the well wall; sigma'θIs the tangential effective stress of the well wall; sigma'rRadial effective stress of the well wall; piThe pressure of the mud-liquid column, R the well-to-well distance, R the radius of the well bore, f the porosity, η the nonlinear correction coefficient, 0.95, delta the seepage coefficient, ξ the intermediate variable,
Figure FDA0002391239950000021
4) determining three deformation failure modes generated in the process of bearing different mud liquid column pressures by the well wall according to Mohr-Coulomb yield criterion and combining the well peripheral stress state and geomechanical parameters, and judging whether the stratum has an upper collapse pressure limit:
for hard and brittle shale formations, the borehole wall has three forms of deformation damage:
(1) when the density of the slurry is lower than the conventional lower limit P of the collapse pressuret|lowerWhen the most dangerous orientation around the well is θ ═90 degrees or 270 degrees, and the difference value of the maximum main stress of the well wall and the minimum main stress of the well wall is sigma'θ-σ'rThe well wall will be damaged by lower limit shearing;
(2) when the mud density is higher than the fracture pressure PfIn the meantime, the most dangerous azimuth around the well is the azimuth of 0 ° or 180 °, and the well wall stress state is σ'θ<0, the well wall will be subjected to tensile damage;
(3) for a particular type of formation, at mud column pressure PiReaches a rupture pressure PfPreviously, the difference between the borehole wall maximum principal stress and the borehole wall minimum principal stress was σ 'at the positions of 0 ° and 180 °'Z-σ'θIf the difference is too high, the upper limit shearing damage will occur to the well wall, and the pressure P of the mud-liquid column at the momentiIs the upper limit P of the collapse pressuret|up
5) Determining the lower limit P of the collapse pressure aiming at three deformation damage forms according to the principle of combination of numbers and shapest|lowerRupture pressure PfAnd upper limit of collapse pressure Pt|upAnd obtaining a safe mud density window of the special stratum as follows:
Pt|lowermud<min(Pt|up,Pf) (4)
in the above formula, ρmudIs the slurry density;
wherein, the lower limit P of the collapse pressure in three deformation damage modes of the well wallt|lowerRupture pressure PfUpper limit of collapse pressure Pt|upThe relative size relationship of the three is as follows:
(1) the tangential and radial stress difference curve of the well wall represents the difference between the maximum main stress of the well wall and the minimum main stress of the well wall at the position of theta being 90 degrees or 270 degrees along with the increase of the pressure of the mud liquid column, and the equation expression of the tangential and radial stress difference curve of the well wall is as follows:
y=3σHh-2x (5)
in the formula, x is a mud liquid column pressure value; y is the stress difference;
tangential and radial stress difference curves and shear failure of well wallsThe abscissa of the intersection point of the required stress difference curve is the lower limit P of the collapse pressuret|lowerThe calculation expression is as follows:
Figure FDA0002391239950000022
wherein UCS is the uniaxial compressive strength of sedimentary rock, UCS is 2. CKKIs the intermediate variable(s) of the variable,K=tan(π/4+φ/2);
(2) the axial and tangential stress difference curve of the well wall represents the difference between the maximum main stress of the well wall and the minimum main stress of the well wall at the position of theta equal to 0 DEG or 180 DEG along with the increase of the pressure of the mud liquid column, and the equation expression of the axial and tangential stress difference curve of the well wall is as follows:
y=σV-(2v-1)σH+(2v-3)σh+x (7)
the abscissa of the intersection point of the axial and tangential stress difference curves of the well wall and the stress difference curve required by shear failure is the upper limit P of collapse pressuret|upSubstituting the formula (4) into the Mohr-Coulomb yield criterion, and finishing to obtain the upper limit P of the collapse pressuret|upThe calculation expression of (a) is as follows:
Figure FDA0002391239950000031
the expression of the transition point (Xo, Yo) between the lower limit shear failure and the upper limit shear failure is as follows:
Figure FDA0002391239950000032
(3) the tangential stress curve of the borehole wall represents the tangential effective stress σ 'at an orientation of 0 or 180 as the mud slurry column pressure increases'θWhen the pressure is reduced to be less than 0, the well wall is subjected to tensile failure, and the abscissa of the intersection point of the tangential stress curve of the well wall and the x axis is fracture pressure PfThe calculation expression is as follows:
Pf=3σhH-αPp+St(10)
wherein St is the tensile strength of the well wall rock.
2. The method of claim 1, wherein the UCS is less than the UCS for a compressive uniaxial strengthpThe formation of (1), which is a special formation that must consider the upper collapse pressure limit, the critical uniaxial compressive strength UCSpCalculated by equation (11):
UCSP=σV-2ν(σHh)-αPp(11)。
CN201710542811.8A 2017-07-05 2017-07-05 Special stratum safe mud density window determination method Active CN107169248B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710542811.8A CN107169248B (en) 2017-07-05 2017-07-05 Special stratum safe mud density window determination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710542811.8A CN107169248B (en) 2017-07-05 2017-07-05 Special stratum safe mud density window determination method

Publications (2)

Publication Number Publication Date
CN107169248A CN107169248A (en) 2017-09-15
CN107169248B true CN107169248B (en) 2020-06-30

Family

ID=59822822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710542811.8A Active CN107169248B (en) 2017-07-05 2017-07-05 Special stratum safe mud density window determination method

Country Status (1)

Country Link
CN (1) CN107169248B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108104804A (en) * 2017-12-11 2018-06-01 西安石油大学 A kind of hard brittle shale Fracturing Pressure Prediction method
CN108843313B (en) * 2018-06-23 2019-07-12 东北石油大学 Shale formation drilling safety drilling fluid density window design method
CN110761784B (en) * 2018-07-26 2023-01-10 中国石油天然气股份有限公司 Method for determining chemical window of well bore working fluid
CN109858147B (en) * 2019-01-30 2022-07-22 西南石油大学 Reliability theory-based quantitative evaluation method for borehole wall instability risk
CN112651113A (en) * 2020-12-17 2021-04-13 中石化石油工程技术服务有限公司 Safe mud density window judgment method suitable for east China sea pressure failure stratum
CN113175322B (en) * 2021-04-28 2023-07-21 中海石油(中国)有限公司 Formation loss pressure profile establishing method, computer device and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102230373A (en) * 2011-05-19 2011-11-02 中国石油大学(北京) Method for determining safe drilling fluid density of coal bed based on structural element analysis model
CN105574609A (en) * 2015-12-11 2016-05-11 中国石油天然气集团公司 Method and device for predicting coal bed safety drilling fluid density window
CN106351650A (en) * 2015-07-16 2017-01-25 中国石油化工股份有限公司 Calculation method of borehole sloughing pressure applicable to the bedding fractured stratum

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070093996A1 (en) * 2005-10-25 2007-04-26 Smith International, Inc. Formation prioritization optimization
CN101392647B (en) * 2008-11-14 2012-11-21 北京石大联创石油新技术有限公司 Borehole wall stability prediction method suitable for gas drilling
BR112016011902A2 (en) * 2013-12-06 2017-08-08 Halliburton Energy Services Inc METHOD AND CONTROL SYSTEM FOR DOWN WELL DRILLING EQUIPMENT, AND COMPUTER READABLE STORAGE MEDIA

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102230373A (en) * 2011-05-19 2011-11-02 中国石油大学(北京) Method for determining safe drilling fluid density of coal bed based on structural element analysis model
CN106351650A (en) * 2015-07-16 2017-01-25 中国石油化工股份有限公司 Calculation method of borehole sloughing pressure applicable to the bedding fractured stratum
CN105574609A (en) * 2015-12-11 2016-05-11 中国石油天然气集团公司 Method and device for predicting coal bed safety drilling fluid density window

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Probabilistic Approach in Wellbore Stability Analysis during Drilling;Mahmood R. Al-Khayari 等;《Journal of Petroleum Engineering》;20161231;第1-14页 *
伊通地区易坍塌地层钻井液安全密度窗口的确定;吴磊 等;《内蒙古石油化工》;20121231;第2012年卷(第12期);摘要、正文第1-3节 *
砂泥岩地层井壁力学稳定性研究及软件编制;雷正义;《中国优秀硕士学位论文全文数据库.工程科技Ⅰ辑》;20040315;第2005年卷(第1期);摘要、第2-6章 *

Also Published As

Publication number Publication date
CN107169248A (en) 2017-09-15

Similar Documents

Publication Publication Date Title
CN107169248B (en) Special stratum safe mud density window determination method
CN108952700B (en) Method for determining anisotropic stratum well wall fracture pressure
CN108468538B (en) Shale hydraulic fracture propagation prediction method
CN104806233B (en) A kind of method for predicting plane of weakness formation collapsed pressure equal yield density window
CN109356567B (en) Method for predicting stability of deep water shallow stratum well wall
CN106126866A (en) Oil-gas well sand prevention method and device based on geomechanical model
CN110821481A (en) Method for evaluating stability of air drilling well wall
CN110580401A (en) method for judging temporary plugging times of segmented multi-cluster fractured well shafts of directional well and horizontal well
Klimov et al. Mechanical-mathematical and experimental modeling of well stability in anisotropic media
CN113109162B (en) Rock fracture initiation pressure calculation method based on thermo-fluid-solid coupling
CN108254262B (en) Rock bedding crack shearing parameter prediction method and device
US10767473B2 (en) Systems and methods for detection of induced micro fractures
CN109726450B (en) Method and equipment for determining repeated fracturing time of shale gas reservoir horizontal well
张杜杰 et al. Mechanisms of sand production from fracture wall and its effect on stress sensitivity in ultra-deep tight sandstone reservoirs
Halafawi et al. Borehole insitu stress stability analysis of RBS-9 field utilizing the inversion technique
CN109033698B (en) Method for calculating fracture pressure of horizontal well in layered stratum
CN116341174A (en) Method and system for predicting critical sand production pressure difference of fractured sandstone reservoir
CN114647913A (en) Maximum level principal ground stress acquisition method based on Mohr-Coulomb criterion
CN114169204A (en) Sand prevention opportunity determination method for offshore oil and gas field development and production
CN113449353A (en) Horizontal well wall stability analysis method considering stratum fluid seepage effect
CN114036703A (en) Method for analyzing strength of salt gypsum rock stratum deep well petroleum casing
CN112651113A (en) Safe mud density window judgment method suitable for east China sea pressure failure stratum
Wei et al. Stability of horizontal coalbed well
CN107916917B (en) Method for improving hydraulic fracturing double-wing crack length
Zhang et al. Effect of stratum on the interfacial stress and structural integrity of cement sheath

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