CN106321082B - For seeking the method and system of salt constituent stratum maximum outside casing load - Google Patents
For seeking the method and system of salt constituent stratum maximum outside casing load Download PDFInfo
- Publication number
- CN106321082B CN106321082B CN201510379772.5A CN201510379772A CN106321082B CN 106321082 B CN106321082 B CN 106321082B CN 201510379772 A CN201510379772 A CN 201510379772A CN 106321082 B CN106321082 B CN 106321082B
- Authority
- CN
- China
- Prior art keywords
- stress
- enhancement region
- crustal
- region
- coefficient
- 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
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of for seeking the method and system of salt constituent stratum maximum outside casing load, this method comprises: determining the stress enhancement region and unaffected zone on salt constituent stratum according to distribution of faults feature and geometric shape；The tectonic stress coefficient of the crustal stress situation identified sign enhancement region of drilling well position based on unaffected zone；Tectonic stress coefficient based on stress enhancement region calculates salt constituent stratum maximum outside casing load.The present invention can Exact Design by the saline bed casing strength of the ground stress anomally of braking effect, can be smoothly in the smooth wellbore construction of saline bed of the ground stress anomally by braking effect.
Description
Technical field
The invention belongs to petroleum engineering well drilling industries, specifically, being related to a kind of for seeking salt constituent stratum maximum set
The method and system of pipe external applied load.
Background technique
Crustal stress is internal stress present in stratum, is an important factor for influencing drilling engineering.In drilling engineering and open
In hair, the size and orientation of the tectonic stress in oil gas reservoir region are grasped, drilling engineering design can be optimized, oil gas field is opened
Well pattern arrangement and control wellbore stability are sent out, serious financial consequences and the person caused by accidents such as leakage is reduced or avoided, sprays, collapse, blocking
Accident etc..
Often assume that stratum is under the effect of three axis crustal stress, the principal stress of three of them principal direction is maximum horizontal principal stress
σH, minimum level crustal stress σhWith vertical stress σz, wherein vertical stress is mainly as caused by the self gravity of stratum, horizontally
Stress is changed by rock mass self weight, tectonic movement, formation fluid pressure and formation temperature to be generated.But it is strong carrying out casing
When degree is checked, the maximum casing external load under common formation condition seldom considers the influence of horizontal crustal stress, in some offices
The saline bed stratum of portion's fault development, ground stress anomally will lead to sleeve design intensity using the common maximum outer calculation method that carries
It is relatively low, cause casing collapse, the huge accident of drilling well failure.
Currently, for salt constituent stratigraphic anormaly crustal stress is influenced by local fracture calculation method and casing outside equivalent load
Check method has carried out some researchs.Document " influence based on measured data and numerical simulation tomography to crustal stress " and " tomography
The Finite-Element Study influenced on stress field " is pointed out to be influenced by being broken, and different well locations not far from one another show that gap is larger
Crustal stress numerical value, the difference of crustal stress states be due to tomography influence caused by, be based on FEM-software ANSYS, utilize three-dimensional
Numerical simulation method carries out inverting to stress field, shows in conjunction with Optimization inversion technology: since tomography influences and
The superposition of compund fault disturbs, and research area's crustal stress states show complicated non-uniform Distribution feature.Document " fault region
Stress field Study on Forecasting Method " points out that the distribution research of crustal stress containing fault region is the challenge with multi-solution, very
Hardly possible is calculated using numerical method, needs to attempt using mechanical model come sunykatuib analysis to apply discontinuous trend surface analysis,
From the angle of geometry deformation, using Thin plate under small deflection theory, fault region is analyzed according to discontinuous trend surface principal curvatures
Stress field Study on Forecasting Method is carried out in the distribution of landform trend surface and tectonic stress field." rift structure is to stress field for document
Influence and its engineering significance " is pointed out to determine that rift structure specifically includes that 1. rift structures influence to the parameter that stress field influences
Range；2. the relationship of local fault construction and regional stress field, including boundary stress ratio KbInfluence and boundary stress side
Influence to the angle α between fracture；3. the physico-mechanical properties for being broken the rock of two sides are (deformation modulus, modulus of shearing, interior
Angle of friction, cohesive force) influence；4. the mechanical property (fracture normal stiffness, shear stiffness, internal friction angle, cohesive force) of fracture
Influence；5. the geometric shape and combining case of tomography.
Currently, saline bed is classified as serious collapsing, expansion, slid, the stratum of wriggling by drilling well handbook, outside saline bed casing
It carries to calculate and use: Pee=ZG0, PeeFor casing external load；Z is the depth of well；G0For overlying formation pressure.This calculation method is applicable in
Make horizontally stress numerical for being influenced to generate larger tectonic stress by local fracture in the situation that tectonic movement is not strong
In the case where substantially exceeding overlying formation pressure, if being designed using the standard can make casing strength obviously too small, cause
The accidents such as casing collapse.
Document above, which all shows to be broken, bigger influence to local crustal stress, and impacted factor also compares more.But
It is the crustal stress Numerical Predicting Method all without a set of suitable field application of proposition, it is more not this to saline bed by crustal stress shadow
The realistic problem for ringing biggish stratum, which carries out the maximum outer calculating for carrying exact value of exact value analysis, especially casing, not to be had yet
It refers to.
Summary of the invention
In order to solve the above problem, the present invention provides a kind of methods for seeking salt constituent stratum maximum outside casing load
And system, accurately to seek the salt constituent stratum maximum outside casing load by braking effect.
According to an aspect of the invention, there is provided a kind of for seeking the side of salt constituent stratum maximum outside casing load
Method, comprising:
The stress enhancement region and unaffected zone on salt constituent stratum are determined according to distribution of faults feature and geometric shape；
Based on the drilling well position of the unaffected zone crustal stress situation determine the unaffected zone maximum and most
Small horizontal crustal stress；
Minimum and maximum horizontal crustal stress based on the unaffected zone determines regional structure stress coefficient, region
The angle of stress ratio and fault strike and region maximum horizontal stress direction；
Minimum and maximum horizontal crustal stress, the region crustal stress ratio and the tomography based on the unaffected zone are walked
Numerical simulator is established to the angle with region maximum horizontal stress direction；
Described in minimum and maximum horizontal crustal stress based on the numerical simulator and the unaffected zone obtains
The minimum and maximum horizontally stress variation by braking effect of stress enhancement region；
It is answered described in the minimum and maximum horizontally stress variation determination by braking effect based on the stress enhancement region
The minimum and maximum horizontal crustal stress by braking effect of power enhancement region；
The stress enhancing is determined based on the minimum and maximum horizontal crustal stress by braking effect of the stress enhancement region
The tectonic stress coefficient in area；It is calculated based on the regional structure stress coefficient and the tectonic stress coefficient of the stress enhancement region disconnected
Split influence coefficient；
Salt constituent stratum maximum casing load is calculated based on the regional structure stress coefficient and the braking effect coefficient
Lotus.
The crustal stress situation drilling well position of the crustal stress situation drilling well position of the crustal stress situation drilling well position of drilling well position
Crustal stress situation according to one embodiment of present invention, the numerical simulator includes the most flood of the stress enhancement region
Level land stress increment percentage-region crustal stress is than submodel, maximum horizontal stress-rift direction and maximum horizontal
The angle submodel of stress direction, minimum level crustal stress increment percentage-region crustal stress are than submodel and minimum level
The angle submodel of stress difference-rift direction and minimum level stress direction.
According to one embodiment of present invention, the maximum based on the numerical simulator and the unaffected zone and
Minimum level crustal stress obtains the minimum and maximum horizontally stress variation by braking effect of the stress enhancement region, drilling well
The crustal stress situation of position further comprises:
Institute is obtained than submodel by maximum horizontal crustal stress increment percentage-region crustal stress of the stress enhancement region
State the maximum horizontal crustal stress increment percentage by braking effect of stress enhancement region；
It is pressed from both sides by the maximum horizontal stress-rift direction and maximum horizontal stress direction of the stress enhancement region
Silver coin model obtains the maximum horizontal stress by braking effect in the stress enhancement region；
Institute is obtained than submodel by minimum level crustal stress increment percentage-region crustal stress of the stress enhancement region
State the minimum level crustal stress increment percentage by braking effect of stress enhancement region；
It is pressed from both sides by the minimum level stress-rift direction and minimum level stress direction of the stress enhancement region
Silver coin model obtains the minimum level stress by braking effect of the stress enhancement region；
By the stress enhancement region by the maximum horizontal crustal stress increment percentage of braking effect, by braking effect
The maximum horizontal principal stress of maximum horizontal stress and the unaffected zone calculates being broken for the stress enhancement region
The maximum horizontal crustal stress variable quantity of influence passes through the minimum level crustal stress increment by braking effect of the stress enhancement region
The minimum level crustal stress calculating of percentage, the minimum level stress and the unaffected zone by braking effect
The minimum level crustal stress variable quantity by braking effect of the stress enhancement region.
According to one embodiment of present invention, the tectonic stress coefficient based on the stress enhancement region calculates salt constituent stratum
Maximum outside casing load is calculate by the following formula to obtain:
Wherein, PeeFor by the salt constituent stratum maximum casing external load of braking effect, σ 'zTo cover pressure on stress enhancement region
Power, T ξ1With T ξ2For the tectonic stress coefficient of stress enhancement region, ξ1And ξ2For regional structure stress coefficient, T is braking effect system
Number, E 'sFor the elasticity modulus of stress enhancement region, υ 'sFor the Poisson's ratio of stress enhancement region, α ' P '0For acting on for stress enhancement region
The strata pressure of skeleton, α ' are the effective stress coefficient of stress enhancement region, P '0For strata pressure, P 'qFor the work of stress enhancement region
Strata pressure for casing.
According to one embodiment of present invention, the braking effect coefficient is calculate by the following formula to obtain:
Wherein, ξ '1With ξ '2For the tectonic stress coefficient of stress enhancement region, ξ1And ξ2For regional structure stress coefficient.
According to one embodiment of present invention, the tectonic stress coefficient of the stress enhancement region is calculate by the following formula to obtain:
Wherein, σ 'HThe maximum horizontal by braking effect for the stress enhancement region being calculated by numerical simulator
Crustal stress, σ 'zFor the burden pressure of stress enhancement region, σ 'hFor the stress enhancement region that is calculated by numerical simulator by
The minimum level crustal stress of braking effect, υ 'sFor the Poisson's ratio of stress enhancement region, E 'sFor the elasticity modulus of stress enhancement region, α '
For the effective stress coefficient of stress enhancement region, P 'pFor the pore pressure of stress enhancement region.
According to one embodiment of present invention, the crustal stress situation based on the unaffected zone determines described uninfluenced
The minimum and maximum horizontal crustal stress in area further comprises:
The crustal stress situation of drilling well position based on the unaffected zone obtains formation rock fracture pressure, instantaneous termination of pumping
Pressure, strata pressure, Tensile Strength of Rock, breakpoint depth, drilling fluid density, effective stress coefficient, break area；
Based on the instantaneous characteristics, the drilling fluid density, described ground breakpoint depth and described break areal calculation
Obtain minimum level crustal stress；
Based on the minimum level crustal stress, the formation rock fracture pressure, the effective stress coefficient, the stratum
Pressure and the Tensile Strength of Rock calculate maximum horizontal crustal stress.
According to another aspect of the present invention, it additionally provides a kind of for seeking salt constituent stratum maximum outside casing load
System, comprising:
Stress division module determines the stress enhancement region and not on salt constituent stratum according to distribution of faults feature and geometric shape
Affected area；
Tectonic stress coefficients calculation block, based on the unaffected zone drilling well position crustal stress situation determine described in
The minimum and maximum horizontal crustal stress of unaffected zone；
Also, the minimum and maximum horizontal crustal stress based on the unaffected zone, determines regional structure stress coefficient, area
Domain crustal stress is than the angle with fault strike and region maximum horizontal stress direction；
Minimum and maximum horizontal crustal stress, the region crustal stress ratio and the tomography based on the unaffected zone are walked
Numerical simulator is established to the angle with region maximum horizontal stress direction；
Described in minimum and maximum horizontal crustal stress based on the numerical simulator and the unaffected zone obtains
The minimum and maximum horizontally stress variation by braking effect of stress enhancement region；
It is answered described in the minimum and maximum horizontally stress variation determination by braking effect based on the stress enhancement region
The minimum and maximum horizontal crustal stress by braking effect of power enhancement region；
The stress enhancing is determined based on the minimum and maximum horizontal crustal stress by braking effect of the stress enhancement region
The tectonic stress coefficient in area；
Maximum outside casing load computing module, the construction based on the regional structure stress coefficient and the stress enhancement region
Stress coefficient calculates braking effect coefficient；Also, it is calculated based on the regional structure stress coefficient and the braking effect coefficient
Salt constituent stratum maximum outside casing load.Beneficial effects of the present invention:
The present invention is directed to by maximum load outside the casing on braking effect ground stress anomally saline bed stratum, according to drilling well standard
It carries out calculating the characteristics of being easy to appear problem, targetedly proposes the set on the stratum by braking effect ground stress anomally saline bed
Manage outer maximum load and seek model, can Exact Design by the saline bed casing strength of the ground stress anomally of braking effect, can
Smoothly in the smooth wellbore construction of saline bed of the ground stress anomally by braking effect.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention can be by specification, right
Specifically noted structure is achieved and obtained in claim and attached drawing.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is required attached drawing in technical description to do simple introduction:
Fig. 1 is method flow diagram according to an embodiment of the invention；
Fig. 2 is a fracture according to an embodiment of the invention and well location schematic diagram；
Fig. 3 a is maximum horizontal crustal stress increment percentage-area of stress enhancement region according to an embodiment of the invention
Domain crustal stress is than submodel schematic diagram；
Fig. 3 b is maximum horizontal stress-fault strike of stress enhancement region according to an embodiment of the invention
With the angle submodel schematic diagram of maximum horizontal stress direction；
Fig. 4 a is minimum level crustal stress increment percentage-area of stress enhancement region according to an embodiment of the invention
Domain crustal stress is than submodel schematic diagram；
Fig. 4 b is minimum level stress-fault strike of stress enhancement region according to an embodiment of the invention
With the angle submodel schematic diagram of minimum level stress direction；
Fig. 5 a is Paleogene System saline bed casing edge load schematic diagram according to an embodiment of the invention；
Fig. 5 b is Paleogene System saline bed casing edge load schematic diagram according to another embodiment of the invention；And
Fig. 6 is Paleogene System saline bed casing maximum external applied load area distribution signal according to an embodiment of the invention
Figure.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to apply to the present invention whereby
Technological means solves technical problem, and the realization process for reaching technical effect can fully understand and implement.It needs to illustrate
As long as not constituting conflict, each feature in each embodiment and each embodiment in the present invention can be combined with each other,
It is within the scope of the present invention to be formed by technical solution.
As shown in Figure 1 for according to the method flow diagram of one embodiment of the present of invention, below with reference to Fig. 1 come to the present invention into
Row is described in detail.
Firstly, in step s 110, determining that the stress on salt constituent stratum enhances according to distribution of faults feature and geometric shape
Area and unaffected zone.
In this step, by geological exploration techniques, form, shape and the combination side in salt constituent stratum breaking region are obtained
Formula, and the stress enhancement region in stratum breaking region and uninfluenced is found out according to the distribution of faults feature and geometric shape on stratum
Area.Wherein, it is stress enhancement region that crustal stress, which is significantly greater than the region of the normal crustal stress in this area,.Apart from tomography farther out, by tomography
Influencing lesser region is unaffected zone.
Next, in the step s 120, the crustal stress situation identified sign of the drilling well position based on unaffected zone enhances
The tectonic stress coefficient in area.
In this step, firstly, the crustal stress situation of the drilling well position based on unaffected zone, determines regional structure stress
Coefficient, region crustal stress are than the angle with fault strike and region maximum horizontal stress direction.Wherein, the brill of unaffected zone
The crustal stress situation that well location is set includes that the data of experiment and core experiment acquisition are broken by ground.The crustal stress situation packet of drilling well position
Include according to well data, survey logging data and by break the crustal stress that formation data that experiment and core experiment obtain etc. obtains
Size and Orientation, specifically, well data includes, well depth opens secondary, strata division situation；Surveying logging data includes sound wave, gal
Horse, density, lithology, wellbore shape etc.；Ground break experimental data include formation rock fracture pressure, instantaneous characteristics, be laminated
Power, Tensile Strength of Rock, breakpoint depth, drilling fluid density, effective stress coefficient, break area；The number that core experiment obtains
According to including the crustal stress size for the specific location that coring point rock core Kaiser effect obtains.
Specifically, in this step, determining area in conjunction with experiment and laboratory core experiment is broken lively using the Seventh Five-Year Plan mode
Domain construction stress coefficient ξ1And ξ2, fault strike and region maximum horizontal stress direction angle, region crustal stress ratio Kb.Its
In, ξ1For region maximum horizontal crustal stress coefficient, i.e. the maximum horizontal crustal stress coefficient of unaffected zone；ξ2For region minimum water
Level land stress coefficient, i.e. the minimum level crustal stress coefficient of unaffected zone.
The Seventh Five-Year Plan mode is to assume that stratum is homogeneous isotropic line elastomer, and assume that phase geological structure is transported after deposit
During dynamic, relative displacement does not occur between stratum and stratum, the strain of all two horizontal directions of stratum is constant, then maximum
Horizontal crustal stress σH, minimum level crustal stress σhWith tectonic stress coefficient ξ1And ξ2Relationship is expressed as follows:
In formula: ξ1And ξ2For regional structure stress coefficient；σH, σh, σzFor region maximum horizontal, minimum level crustal stress and on
Pressure is covered, maximum burden pressure is calculated according to the value of density log or acoustic logging in well-log information；PpFor hole
Gap pressure, is calculated according to log data；υs、EsFor regional stratum static state Poisson's ratio and elasticity modulus, according to well-log information into
Row calculates；α is region effective stress coefficient, is calculated according to well-log information.
It, can be with by breaking experimental data lively by fracture pressure and maximum, minimum level crustal stress numerical relation
Inverse goes out tectonic stress coefficient, not constructed coefficient as regional stress by the tectonic stress coefficient of braking effect well location.
Specifically, the fracture pressure of unaffected zone and maximum, minimum level crustal stress are as follows:
Pf=3 σh-σH-αPp+St (2)
Wherein, fracture pressure PfExperimental calculation is broken by ground to obtain.
By instantaneous characteristics, drilling fluid density, breakpoint depth and ground break areal calculation and answer with obtaining minimum level
Minimum level crustal stress is calculated based on following formula in power:
σh=Ps+10-6×ρmgHw (3)
It is strong by minimum level crustal stress, formation rock fracture pressure, effective stress coefficient, strata pressure and rock tension
Degree calculates maximum horizontal crustal stress, and maximum horizontal crustal stress is calculated based on following formula:
σH=3 σh-Pf-α·Pp+|St| (4)
Wherein, StIt is acquired for Tensile Strength of Rock by rock mechanics parameters calculating, PsInstantaneous characteristics break experiment by ground and break
The reading of curve data figure is split, Pf is by the formation rock fracture pressure that obtains on rupture curve, strata pressure Pp, effective stress coefficient
α is obtained, ρ by formation pressure calculation related datam, HwRespectively drilling fluid density and ground breakpoint depth.Thus it obtains not by shadow
Ring the ground breakpoint in area: σhMinimum level crustal stress, σHThe numerical value of maximum horizontal crustal stress, brings following formula into:
Regional structure stress coefficient ξ can be obtained1And ξ2。
At the same time it can also which zone boundary stress axis is calculated:
The fault strike of unaffected zone and the angle of maximum horizontal stress direction are calculated by log data.
Next, based on regional structure stress coefficient, region crustal stress ratio (zone boundary stress axis), fault strike with
The angle of region maximum horizontal stress direction, the tectonic stress coefficient of identified sign enhancement region.
In this step, the well drilling positions firstly, according to Theory of Fracture Mechanics, using discrete element method based on unaffected zone
The crustal stress situation set, region crustal stress are than establishing Numerical-Mode with the angle of fault strike and region maximum horizontal stress direction
Analog model is found out under various influence conditions, by the stress situation of braking effect.The numerical simulator includes stress
Maximum horizontal crustal stress increment percentage-region crustal stress of enhancement region is than submodel, maximum horizontal stress-fracture side
Compare submodule to angle submodel, the minimum level crustal stress increment percentage-region crustal stress with maximum horizontal stress direction
The angle submodel of type and minimum level stress-rift direction and minimum level stress direction.
It is obtained than submodel certain by maximum horizontal crustal stress increment percentage-region crustal stress of stress enhancement region
Maximum horizontal of the region crustal stress than the angle and uninfluenced region of, certain rift direction and maximum horizontal stress direction
The maximum horizontal crustal stress increment by braking effect of stress enhancement region under the conditions of stress (i.e. region maximum horizontal crustal stress)
Percentage.
Pass through maximum horizontal stress-maximum horizontal stress direction and rift direction angle of stress enhancement region
Answering under the conditions of submodel is obtained in certain rift direction with maximum horizontal stress direction angle, in certain area crustal stress ratio
The maximum horizontal stress by braking effect of power enhancement region.
It is obtained than submodel certain by minimum level crustal stress increment percentage-region crustal stress of stress enhancement region
Region crustal stress is answered than, certain rift direction with the minimum level in maximum horizontal stress direction angle and uninfluenced region
The minimum level crustal stress increment percentage by braking effect of stress enhancement region under the conditions of power.
Pass through minimum level stress-maximum horizontal stress direction and rift direction angle of stress enhancement region
Answering under the conditions of submodel is obtained in certain rift direction with maximum horizontal stress direction angle, in certain area crustal stress ratio
The minimum level stress by braking effect of power enhancement region.
By stress enhancement region by the maximum horizontal crustal stress increment percentage of braking effect, by the maximum of braking effect
The most flood by braking effect of the maximum horizontal principal stress of horizontal stress and unaffected zone calculating stress enhancement region
Level land stress variation, by stress enhancement region by the minimum level crustal stress increment percentage of braking effect, by fracture shadow
The minimum level crustal stress of loud minimum level stress and unaffected zone calculate stress enhancement region by braking effect
Minimum level crustal stress variable quantity.
The maximum horizontal crustal stress variable quantity and minimum level crustal stress by braking effect based on stress enhancement region changes
The maximum horizontal principal stress and minimum level crustal stress of stress enhancement region is calculated in amount.By the stress enhancement region being calculated
Maximum horizontal principal stress and minimum level crustal stress substitute into the Seventh Five-Year Plan model inverse, i.e. substitution formula (5) obtains the structure of stress enhancement region
Make stress coefficient ξ '1With ξ '2, stress coefficient ξ ' through this construction1Or ξ '2, the tectonic stress coefficient ξ of unaffected zone1Or ξ2, lead to
It crosses formula (7) and braking effect coefficient is calculated:
Wherein, the tectonic stress coefficient ξ ' of stress enhancement region is calculated by formula (5)1With ξ '2When:
Wherein σ 'H, the maximum horizontal by braking effect for the stress enhancement region being calculated by numerical simulator
Stress；σ′hThe minimum level crustal stress by braking effect for the stress enhancement region being calculated by numerical simulator；υ′s
For the stratum static state Poisson's ratio of stress enhancement region；E′sFor the elasticity modulus of stress enhancement region；α ' has effect for stress enhancement region
Force coefficient, P 'pFor the pore pressure of stress enhancement region.When stress enhancement region is with drilling well, the brill of drilling well can use
The crustal stress situation data that well location is set seek the corresponding tectonic stress coefficient of the well location, then the tectonic stress obtained based on the well location
Coefficient is to the tectonic stress coefficient ξ ' being calculated1Or ξ '2It is modified.Specifically, the flat of these tectonic stress coefficients can be used
Mean value, or weight coefficient is designed to seek final tectonic stress coefficient ξ '1With ξ '2。
Finally, in step s 130, the tectonic stress coefficient based on stress enhancement region calculates salt constituent stratum maximum casing
External applied load.
Salt constituent stratum maximum outside casing load is calculate by the following formula to obtain:
Wherein, PeeFor by the salt constituent stratum maximum casing external load of braking effect, T ξ1With T ξ2For stress enhancement region
Tectonic stress coefficient, ξ1And ξ2For regional structure stress coefficient, T is braking effect coefficient, E 'sFor the springform of stress enhancement region
Amount, υ 'sFor the Poisson's ratio of stress enhancement region, α ' P '0For the strata pressure for acting on skeleton of stress enhancement region, α ' is stress increasing
The effective stress coefficient in strong area, P '0For strata pressure, P 'qFor the strata pressure for acting on casing of stress enhancement region.
Feasibility of the invention is verified below by way of a specific example.It is illustrated in figure 2 and is locally broken
It splits the specific well location of salt constituent stratum block of influence and is broken distribution mode and geometry schematic diagram.As shown in Fig. 2, influencing wheat
The fracture of Gai Tiyi block is mainly the fracture of Ma south, and No. 5 structural belts, main body is two parallel fractures, and outside is lower wall, is broken it
Interior is upper disk.
It as shown in Fig. 2, jade Bei1Jing, Yu Bei 1-1 well, Yu Bei 1-2 well are between two parallel fractures, and is fracture
Upper disk position.Beautiful 6 well of Bei5Jing, Yu Bei is in the two sides of parallel fracture, and belongs to the lower wall position of fracture.According to fracture
Influence to local crustal stress can judge beautiful Bei5Jing and beautiful northern 6 well stress situations and jade Bei1Jing, Yu Bei 1-1 well, Yu Bei
The stress situation of 1-2 well has a difference, and the size of difference is by the corner dimension of rift direction and region biggest principal stress direction, disconnected
Split the influence of the rock strength between two sides and fracture, the internal friction angle size for being broken itself.
According to the actual conditions of beautiful northern 1 block, there is different degrees of casing collapse in 5 well of Yu Bei and jade Bei6Jing
Situation illustrates that the crustal stress for the two side areas being broken under conditions of local structure influences close to Ma south will be significantly greater than the ground
The normal crustal stress in area illustrates that this two well location enhances region in stress.And locating for jade Bei1Jing, Yu Bei 1-1 well and beautiful north 1-2 well
Two fracture belts inside influenced by local structure stress, the crustal stress in the region should be lower than the normal crustal stress in this area
Situation, this three well location is in stress weakened region.
Beautiful 8 well of Bei4Jing, Yu Bei and jade Bei9Jing, farther out apart from tomography, crustal stress is mainly by overburden pressure and region
Construction determines, is influenced by local structure smaller, belongs to the normal crustal stress situation in this area, the well not influenced by local structure
Position is located at unaffected zone.
Beautiful Bei7Jing is broken north side in Ma south, is broken northwest side in a single line, should answer on the whole with beautiful northern 5 wells
The trend of power is consistent, but relatively weak, by local structure influenced the result is that normal crustal stress should be slightly higher than.Beautiful north 3
Well from Ma south is main due to being broken farther out, and nearby offset well is less, and concrete condition also needs to continue to judge.
Each well is as shown in table 1 with respect to the specific location of tomography geometric shape in Fig. 2.
Table 1
Next, being calculated according to the crustal stress situation of the drilling well position in uninfluenced region based on formula (2), (3), (4)
To the minimum level crustal stress and maximum horizontal crustal stress of unaffected zone.It is then based on formula (5) and (6) is calculated not by shadow
Ring the tectonic stress coefficient ξ in area1=0.778 and ξ2=0.368, and then pass through formula (6)
Zone boundary stress axis is calculatedIt can also be obtained by the crustal stress situation of drilling well position
It is nearly North and South direction (N4 °/184 °) to region maximum horizontal stress direction, fracture direction and principal stress angle are about 45 °.Into
And region minimum level stress direction is obtained, fracture direction and minimum level crustal stress angle.
Based on numerical simulator according to beautiful backlands region stress ratio KbValue about 1.26, fault strike and maximum horizontal
Crustal stress angle is about 45 °, according to the numerical simulation result of Fig. 3 a, 3b and 4a, 4b, according to the concrete condition of embodiment block,
Know that influence of the fracture to crustal stress is about 25% or so of region crustal stress, by the maximum principal stress difference under the conditions of the angle
The maximum horizontal crustal stress variable quantity that stress enhancement region can be obtained multiplied by 25%.The minimum level crustal stress of stress enhancement region changes
Amount obtains by the above process.Bring the maximum horizontal crustal stress for the stress enhancement region being calculated and minimum level crustal stress into formula
(5) it obtains: the tectonic stress coefficient of stress enhancement region: ξ1=0.844, ξ2=0.399, braking effect is calculated based on formula (7)
Coefficient T.
Braking effect coefficient T will be calculated and substitute into formula (8), load value outside a maximum of stress enhancement region can be obtained.
Calculate the outer load value of all maximum and the therefrom maximum value that selects of stress enhancement region by formula (8), as required by take
Maximum outside casing load.As shown in two mouthfuls of well edge load schematic diagrames of Fig. 5 a and Fig. 5 b, all casing boundaries can be chosen and carried
Foundation of the maximum value of lotus as design outside casing load.Finally obtained external applied load numerical value section in figure as shown in fig. 6, as put down
Shown in the bumps in face, wherein the part more protruded in plane indicates that external applied load is maximum.External applied load largest portion corresponding diagram 2 in Fig. 6
Middle stress reinforces the northern well 5 of jade and 6 region of beautiful northern well in area.
According to another aspect of the present invention, it additionally provides a kind of for seeking salt constituent stratum maximum outside casing load
System, the system include stress division module, tectonic stress coefficients calculation block and maximum outside casing load computing module.
Wherein, stress division module determines that the stress on salt constituent stratum enhances according to distribution of faults feature and geometric shape
Area and unaffected zone；Tectonic stress coefficients calculation block, the crustal stress situation of the drilling well position based on the unaffected zone
Determine the tectonic stress coefficient of the stress enhancement region；Maximum outside casing load computing module, based on the stress enhancement region
Tectonic stress coefficient calculates salt constituent stratum maximum outside casing load.
The mode that the present invention is combined by the way of qualitative and quantitative proposes a set of reasonable relatively accurate crustal stress
Numerical value acquiring method, can Exact Design by the ground stress anomally of braking effect saline bed casing strength, can smoothly by
The smooth wellbore construction of the saline bed of the ground stress anomally of braking effect.
While it is disclosed that embodiment content as above but described only to facilitate understanding the present invention and adopting
Embodiment is not intended to limit the invention.Any those skilled in the art to which this invention pertains are not departing from this
Under the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details,
But scope of patent protection of the invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (8)
1. a kind of method for seeking salt constituent stratum maximum outside casing load, comprising:
The stress enhancement region and unaffected zone on salt constituent stratum are determined according to distribution of faults feature and geometric shape；
The minimum and maximum water of the unaffected zone is determined based on the crustal stress situation of the drilling well position of the unaffected zone
Level land stress；
Minimum and maximum horizontal crustal stress based on the unaffected zone, determines regional structure stress coefficient, region crustal stress
Than the angle with fault strike and region maximum horizontal stress direction；
Minimum and maximum horizontal crustal stress, the region crustal stress ratio and the fault strike based on the unaffected zone with
The angle of region maximum horizontal stress direction establishes numerical simulator；
Minimum and maximum horizontal crustal stress based on the numerical simulator and the unaffected zone obtains the stress
The minimum and maximum horizontally stress variation by braking effect of enhancement region；
Determine that the stress increases based on the minimum and maximum horizontally stress variation by braking effect of the stress enhancement region
The minimum and maximum horizontal crustal stress by braking effect in strong area；
The stress enhancement region is determined based on the minimum and maximum horizontal crustal stress by braking effect of the stress enhancement region
Tectonic stress coefficient；
Braking effect coefficient is calculated based on the regional structure stress coefficient and the tectonic stress coefficient of the stress enhancement region；
Salt constituent stratum maximum outside casing load is calculated based on the regional structure stress coefficient and the braking effect coefficient.
2. the method according to claim 1, wherein the numerical simulator includes the stress enhancement region
Maximum horizontal crustal stress increment percentage-region crustal stress is than submodel, maximum horizontal stress-rift direction and maximum
Horizontal stress direction angle submodel, minimum level crustal stress increment percentage-region crustal stress are than submodel and minimum water
Level land stress difference-rift direction and minimum level stress direction angle submodel.
3. according to the method described in claim 2, it is characterized in that, based on the numerical simulator and described uninfluenced
The minimum and maximum horizontal crustal stress in area obtains the minimum and maximum horizontal crustal stress by braking effect of the stress enhancement region
Variable quantity further comprises:
By maximum horizontal crustal stress increment percentage-region crustal stress of the stress enhancement region than being answered described in submodel acquisition
The maximum horizontal crustal stress increment percentage by braking effect of power enhancement region；
Pass through maximum horizontal stress-rift direction and maximum horizontal stress direction angle of the stress enhancement region
Model obtains the maximum horizontal stress by braking effect in the stress enhancement region；
By minimum level crustal stress increment percentage-region crustal stress of the stress enhancement region than being answered described in submodel acquisition
The minimum level crustal stress increment percentage by braking effect of power enhancement region；
Pass through minimum level stress-rift direction and minimum level stress direction angle of the stress enhancement region
Model obtains the minimum level stress by braking effect of the stress enhancement region；
By the stress enhancement region by the maximum horizontal crustal stress increment percentage, described by braking effect of braking effect
The maximum horizontal crustal stress of maximum horizontal stress and the unaffected zone calculates being broken for the stress enhancement region
The maximum horizontal crustal stress variable quantity of influence passes through the minimum level crustal stress increment by braking effect of the stress enhancement region
The minimum level crustal stress calculating of percentage, the minimum level stress and the unaffected zone by braking effect
The minimum level crustal stress variable quantity by braking effect of the stress enhancement region.
4. method according to any one of claim 1-3, which is characterized in that the construction based on the stress enhancement region is answered
Force coefficient calculates salt constituent stratum maximum outside casing load and is calculate by the following formula to obtain:
Wherein, PeeFor by the salt constituent stratum maximum casing external load of braking effect, σ 'zFor the burden pressure of stress enhancement region, T
ξ1With T ξ2For the tectonic stress coefficient of stress enhancement region, ξ1And ξ2For regional structure stress coefficient, T is braking effect coefficient, E 's
For the elasticity modulus of stress enhancement region, υ 'sFor the Poisson's ratio of stress enhancement region, α ' P '0For the skeleton that acts on of stress enhancement region
Strata pressure, α ' are the effective stress coefficient of stress enhancement region, P '0For strata pressure, P 'qSet is acted on for stress enhancement region
The strata pressure of pipe.
5. according to the method described in claim 4, it is characterized in that, the braking effect coefficient is calculate by the following formula to obtain:
Wherein, ξ '1With ξ '2For the tectonic stress coefficient of stress enhancement region, ξ1And ξ2For regional structure stress coefficient.
6. according to the method described in claim 5, it is characterized in that, the tectonic stress coefficient of the stress enhancement region passes through following formula
It is calculated:
Wherein, σ 'HFor the stress enhancement region being calculated by numerical simulator the maximum horizontal by braking effect answer
Power, σ 'zFor the burden pressure of stress enhancement region, σ 'hFor being broken for the stress enhancement region that is calculated by numerical simulator
The minimum level crustal stress of influence, υ 'sFor the Poisson's ratio of stress enhancement region, E 'sFor the elasticity modulus of stress enhancement region, α ' is to answer
The effective stress coefficient of power enhancement region, P 'pFor the pore pressure of stress enhancement region.
7. the method according to claim 1, wherein the crustal stress situation based on the unaffected zone determines institute
The minimum and maximum horizontal crustal stress for stating unaffected zone further comprises:
The crustal stress situation of drilling well position based on the unaffected zone obtains formation rock fracture pressure, instantaneous termination of pumping pressure
Power, strata pressure, Tensile Strength of Rock, breakpoint depth, drilling fluid density, effective stress coefficient, break area；
Based on the instantaneous characteristics, the drilling fluid density, described ground breakpoint depth and described break areal calculation and obtain
Minimum level crustal stress；
Based on the minimum level crustal stress, the formation rock fracture pressure, the effective stress coefficient, the strata pressure
Maximum horizontal crustal stress is calculated with the Tensile Strength of Rock.
8. a kind of system for seeking salt constituent stratum maximum outside casing load, comprising:
Stress division module determines the stress enhancement region on salt constituent stratum according to distribution of faults feature and geometric shape and not by shadow
Ring area；
Tectonic stress coefficients calculation block, based on the drilling well position of the unaffected zone crustal stress situation determine it is described not by
The minimum and maximum horizontal crustal stress of the zone of influence；
Also, minimum and maximum horizontal crustal stress based on the unaffected zone determines regional structure stress coefficient, region
The angle of stress ratio and fault strike and region maximum horizontal stress direction；
Minimum and maximum horizontal crustal stress, the region crustal stress ratio and the fault strike based on the unaffected zone with
The angle of region maximum horizontal stress direction establishes numerical simulator；
Minimum and maximum horizontal crustal stress based on the numerical simulator and the unaffected zone obtains the stress
The minimum and maximum horizontally stress variation by braking effect of enhancement region；
Determine that the stress increases based on the minimum and maximum horizontally stress variation by braking effect of the stress enhancement region
The minimum and maximum horizontal crustal stress by braking effect in strong area；
The stress enhancement region is determined based on the minimum and maximum horizontal crustal stress by braking effect of the stress enhancement region
Tectonic stress coefficient；
Maximum outside casing load computing module, the tectonic stress based on the regional structure stress coefficient and the stress enhancement region
Coefficient calculates braking effect coefficient, also, calculates salt cream based on the regional structure stress coefficient and the braking effect coefficient
Rock stratum maximum outside casing load.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379772.5A CN106321082B (en) | 2015-07-01 | 2015-07-01 | For seeking the method and system of salt constituent stratum maximum outside casing load |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379772.5A CN106321082B (en) | 2015-07-01 | 2015-07-01 | For seeking the method and system of salt constituent stratum maximum outside casing load |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106321082A CN106321082A (en) | 2017-01-11 |
| CN106321082B true CN106321082B (en) | 2019-05-10 |
Family
ID=57726863
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510379772.5A Active CN106321082B (en) | 2015-07-01 | 2015-07-01 | For seeking the method and system of salt constituent stratum maximum outside casing load |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106321082B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108412478A (en) * | 2018-01-27 | 2018-08-17 | 东北石油大学 | A kind of device and analogy method of simulation mud stone water suction casing damage |
| CN110107276B (en) * | 2018-01-30 | 2022-06-03 | 中国石油天然气股份有限公司 | Casing design method and apparatus |
| CN110032777B (en) * | 2019-03-25 | 2021-09-21 | 西南石油大学 | Drilling fluid density safety window estimation method based on uncertainty analysis |
| CN113550727B (en) * | 2021-08-25 | 2022-10-21 | 中国石油大学(北京) | Method for calculating strength of volume fracturing sleeve in fracture development area |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1594829A (en) * | 2004-06-18 | 2005-03-16 | 大庆油田有限责任公司 | Method for preventing casing wear in oil field |
| WO2012106227A3 (en) * | 2011-01-31 | 2012-11-15 | Baker Hughes Incorporated | System and method for determining pressure transition zones |
| CN103195374A (en) * | 2013-03-27 | 2013-07-10 | 西南石油大学 | Method for designing salt rock stratum irregular-shaped sleeve |
| EP2856387A1 (en) * | 2012-05-30 | 2015-04-08 | Landmark Graphics Corporation | System and method for reservoir simulation optimization |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7178391B2 (en) * | 2002-10-31 | 2007-02-20 | Battelle Energy Alliance, Llc | Insertion tube methods and apparatus |
-
2015
- 2015-07-01 CN CN201510379772.5A patent/CN106321082B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1594829A (en) * | 2004-06-18 | 2005-03-16 | 大庆油田有限责任公司 | Method for preventing casing wear in oil field |
| WO2012106227A3 (en) * | 2011-01-31 | 2012-11-15 | Baker Hughes Incorporated | System and method for determining pressure transition zones |
| EP2856387A1 (en) * | 2012-05-30 | 2015-04-08 | Landmark Graphics Corporation | System and method for reservoir simulation optimization |
| CN103195374A (en) * | 2013-03-27 | 2013-07-10 | 西南石油大学 | Method for designing salt rock stratum irregular-shaped sleeve |
Non-Patent Citations (4)
| Title |
|---|
| 大庆油田葡北X断块套管损坏预测方法及矿场试验;李进旺 等;《石油学报》;20080731;625-628页 |
| 深井盐层套管非均匀载荷计算与套管设计方法研究;邹灵战 等;《石油钻探技术》;20080131;23-27页 |
| 玉北地区套管变形分析及措施;胡广强;《中国西部科技》;20140531;33-34页 |
| 蠕变地层套管载荷分析研究;殷有泉 等;《岩石力学与工程学报》;20040731;2381-2384页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106321082A (en) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104500050B (en) | Fractured reservoir permeability tensor and anisotropy quantitative prediction method | |
| CN105672971B (en) | A kind of Reservoir Fracture Opening pressure, open order and water injection pressure Forecasting Methodology | |
| Guo et al. | Simulation of three-dimensional tectonic stress fields and quantitative prediction of tectonic fracture within the Damintun Depression, Liaohe Basin, northeast China | |
| Ju et al. | The present-day in-situ stress field within coalbed methane reservoirs, Yuwang Block, Laochang Basin, south China | |
| CN106484958B (en) | A kind of system of the three-dimensional cracking permeability tensor computation model based on pit shaft gap observation | |
| CN106321082B (en) | For seeking the method and system of salt constituent stratum maximum outside casing load | |
| Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
| CN106291755B (en) | A kind of areas Long Sheng low-grade fault law of development quantitative forecasting technique | |
| CN103914620B (en) | Method for computing distribution of opening spaces of fractures of fault rupture zones | |
| CN106014399A (en) | Method for establishing high-precision three-dimensional ground stress model for heterogeneous formation | |
| CN110969699B (en) | High-precision three-dimensional ground stress modeling method and application thereof | |
| Feng et al. | Near-surface stress measurements in the Longmenshan fault belt after the 2008 Wenchuan Ms8. 0 earthquake | |
| Wang et al. | Adaptive finite element‐discrete element analysis for microseismic modelling of hydraulic fracture propagation of perforation in horizontal well considering pre‐existing fractures | |
| Ren et al. | Simulation of stress fields and quantitative prediction of fractures distribution in upper Ordovician biological limestone formation within Hetianhe field, Tarim Basin, NW China | |
| US20140067269A1 (en) | Undisturbed fracture orientation | |
| Yan et al. | In-situ stress partition and its implication on coalbed methane occurrence in the basin–mountain transition zone: a case study of the Pingdingshan coalfield, China | |
| Sandiford et al. | Hydrogeological implications of active tectonics in the Great Artesian Basin, Australia | |
| CN106815412B (en) | Simulation method and device for structural stress field | |
| Lin et al. | In situ stress variation and coal reservoir permeability response of the Jurassic Yan'an formation in the southwestern Ordos basin, China: Its impact on coalbed methane development | |
| Ienaga et al. | Borehole image analysis of the Nankai Accretionary Wedge, ODP Leg 196: Structural and stress studies | |
| Dembo et al. | Mechanical contrast and asymmetric distribution of crustal deformation across plate boundaries: Insights from the northern Dead Sea fault system | |
| Zhao et al. | Present in-situ stress distribution characteristics of strike-slip in SH Oilfield, Tarim Basin | |
| Feng et al. | An improved geomechanical model for the prediction of fracture generation and distribution in brittle reservoirs | |
| Fay et al. | Small‐scale upper mantle convection and crustal dynamics in southern California | |
| Roberts | Fault orientation variations along the strike of active normal fault systems in Italy and Greece: Implications for predicting the orientations of subseismic-resolution faults in hydrocarbon reservoirs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |