CN107169192A - Method and device for identifying size of karst cave while drilling - Google Patents
Method and device for identifying size of karst cave while drilling Download PDFInfo
- Publication number
- CN107169192A CN107169192A CN201710326895.1A CN201710326895A CN107169192A CN 107169192 A CN107169192 A CN 107169192A CN 201710326895 A CN201710326895 A CN 201710326895A CN 107169192 A CN107169192 A CN 107169192A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- solution cavity
- depth
- size
- energy density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005553 drilling Methods 0.000 title abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 83
- 239000011435 rock Substances 0.000 claims description 56
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 25
- 238000002474 experimental method Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000004323 axial length Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 206010019233 Headaches Diseases 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000004209 hair Anatomy 0.000 description 1
- 231100000869 headache Toxicity 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a method and a device for identifying the size of a karst cave while drilling, wherein a mechanical model of a reservoir containing the karst cave with a preset size is established; according to the mechanical model of the reservoir, with the karst cave as a center, calculating the ground stress field of the reservoir distributed around the karst cave; calculating the relation between the depth of the reservoir and the distortion energy density of the reservoir according to the ground stress field; determining a theoretical value between the bit reaction torque and the reservoir depth according to the relation between the reservoir depth and the reservoir distortion energy density and the relation between experimentally determined distortion energy density and the bit reaction torque; obtaining an actual measurement value of the bit reaction torque of the target karst cave reservoir stratum, comparing the actual measurement value with the theoretical value, determining a mechanical model corresponding to the target karst cave reservoir stratum, and determining the size of the karst cave. The method can accurately predict the size of the karst cave in the reservoir.
Description
Technical field
The present invention relates to oil-gas exploration and development field, more particularly to a kind of solution cavity size with boring recognition methods and device.
Background technology
Reservoir is the place of Oil And Gas Occurrence, is also the direct target zone of oil-gas exploration and development.For example, in carbonate reservoir
In the development of common solution cavity, the size of solution cavity is from several meters to tens of rice, due to that whether can not predict front in drilling process
The size for meeting and being bored in front of solution cavity or None- identified and meeting solution cavity will be bored, serious production accident and economic loss is frequently resulted in.
Recognize that solution cavity size is always international headache in reservoir in drilling process, therefore, need a kind of accurate predicting reservoir badly
The method of middle solution cavity size.
The content of the invention
The present invention provide a kind of solution cavity size with recognition methods and device is bored, to solve to recognize reservoir in drilling process
Middle this problem of solution cavity size, realizes solution cavity size in accurate predicting reservoir.
One side of the invention provide a kind of solution cavity size with boring recognition methods, including：
According to multiple default solution cavity sizes, the mechanical model of its reservoir is set up to the solution cavity of every kind of size；
According to the mechanical model of the reservoir, centered on the solution cavity, calculate and surround the reservoir that the solution cavity is distributed
Stress field；
Relation between depth of reservoirs and reservoir distortion energy density is obtained according to stress field calculating；
Relation that can be between density according to the depth of reservoirs and reservoir distortion, and measuring distortion energy density with
Relation between drill bit reaction torque, determines the theoretical values between the drill bit reaction torque and the depth of reservoirs；
Obtain target carbonate rock drill bit reaction torque and depth of reservoirs between measured value, by the measured value with
The theoretical values are compared, and determine the corresponding mechanical model of the target carbonate rock, obtain the mechanical model correspondence
Solution cavity size.
Another aspect of the invention provide a kind of solution cavity size with boring identifying device, including build module, first calculate
Module, the second computing module, the 3rd computing module, identification module；
Wherein, the structure module, for according to multiple default solution cavity sizes, being set up to the solution cavity of every kind of size
The mechanical model of its reservoir；
First computing module, for the mechanical model according to the reservoir, centered on the solution cavity, calculating is surrounded
The stress field of the reservoir of the solution cavity distribution；
Second computing module, depth of reservoirs and reservoir distortion energy density are obtained for being calculated according to the stress field
Between relation；
3rd computing module, for the relation between the energy density that distorted according to the depth of reservoirs and reservoir, and
The distortion of measuring can be between density and drill bit reaction torque relation, determine the drill bit reaction torque and the depth of reservoirs it
Between theoretical values；
The identification module, for obtaining the actual measurement number between the drill bit reaction torque of target carbonate rock and depth of reservoirs
Value, the measured value is compared with the theoretical values, the corresponding mechanical model of the target carbonate rock is determined, obtains
To the corresponding solution cavity size of the mechanical model.
As shown from the above technical solution, the mechanical model that the present invention passes through reservoir of the foundation comprising pre-set dimension solution cavity；Root
According to the mechanical model of the reservoir, centered on the solution cavity, the stress field for the reservoir being distributed around the solution cavity is calculated；Root
The relation between depth of reservoirs and reservoir distortion energy density is calculated according to the stress field；It is abnormal according to the depth of reservoirs and reservoir
Relation between the relation become between energy density, and the distortion energy density and drill bit reaction torque of measuring, determines the brill
Theoretical values between head reaction torque and the depth of reservoirs；The measured value of the drill bit reaction torque of target carbonate rock is obtained,
The measured value is compared with the theoretical values, the corresponding mechanical model of the target carbonate rock is determined, it is determined that
Solution cavity size.The present invention can accurate solution cavity size in predicting reservoir.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are this hairs
Some bright embodiments, for those of ordinary skill in the art, without having to pay creative labor, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is a kind of solution cavity size that the embodiment of the present invention one is provided with the schematic flow sheet for boring recognition methods；
Fig. 2 is a kind of solution cavity size that the embodiment of the present invention two is provided with the schematic flow sheet for boring recognition methods；
Fig. 3 is a kind of solution cavity size that the embodiment of the present invention three is provided with the structural representation for boring identifying device；
Fig. 4 is a kind of solution cavity size that the embodiment of the present invention four is provided with the structural representation for boring identifying device.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art
Acquired every other embodiment, belongs to the scope of protection of the invention under the premise of creative work is not made.
First, to reservoir of the present invention stress field, distortion can the relation of density and drill bit reaction torque say
It is bright.
Rock by factors such as underground heat, rotation, gravitation, extruding because being influenceed, and its original shape to be recovered outwardly produced
The active force of state, referred to as crustal stress.Reservoir rock is acted on by crustal stress, in the certain area centered on reservoir solution cavity hole
Form stress field.Stress field is relevant with geological form, for example the weight of overlying material, depth where from earth's surface to reservoir,
Variety classes rock is distributed with different depth, and its rock density and thickness are different, therefore affects the crustal stress of different depth
.Meanwhile, solution cavity is of different sizes, and the stress field on its periphery is also different.
When rock is acted on by external force, its inner product accumulation of energy amount has stored certain deformation energy, the distortion of meaning can density.
The difference of the stress field on solution cavity periphery, the distortion energy density of correspondence solution cavity periphery each point is different.It is real by using the principle of similitude
Test the affecting laws for determining reservoir distortion energy density to drill bit reaction torque, it is found that reservoir distortion energy density is big, and reservoir is fragile
Broken, the reaction torque that drill bit is produced is larger；Reservoir distortion energy density is small, and reservoir is non-breakable, and the reaction torque that drill bit is produced is smaller.By
This sets up the contact of the drill bit reaction torque of solution cavity size corresponding thereto.
Based on this feature, the drill bit reaction torque of experimental analysis different size solution cavity periphery each point, and in actual well drilled mistake
Drill bit reaction torque changing rule, the drill bit reaction torque pair with the different size solution cavity periphery each point obtained by experimental analysis are monitored in journey
Than so that solution cavity size is accurately recognized, while proposing targetedly drilling well and development plan.
Embodiment one
Fig. 1 is a kind of solution cavity size that the embodiment of the present invention one is provided with the schematic flow sheet for boring recognition methods, such as Fig. 1
Shown, the method for the present embodiment can include：
Step 101：According to multiple default solution cavity sizes, the mechanics mould of its reservoir is set up to the solution cavity of every kind of size
Type；
The mechanical model of reservoir can be the geometrical model of various shapes, for example, can include different size solution cavity to be internal
Square body Model.The solution cavity included in square body Model is elliposoidal, and different size of solution cavity has different long axis lengths
And minor axis length, according to the combination of multiple default long axis lengths and minor axis length, solution cavity size is determined, and then to every kind of
The solution cavity of the combination of long axis length and minor axis length sets up the mechanical model of its reservoir.Meanwhile, comprising different size solution cavity just
The length of side of cube model is identical.
X, y-axis are horizontal direction axle in square body Model, and z-axis is vertical axes, the geometrical model and border of solution cavity and country rock
Condition, square body Model side is along x, y-axis direction by the squeezing action of surrounding rock, and square bottom surface is along z-axis negative direction by weight
Power is acted on；Square body Model top surface is by overburden pressure.Therefore five faces constraint of the square body Model in addition to top surface is square
Body Model is along the displacement on three directions of x, y, z.
Meanwhile, using finite element method on square body Model grid division, by square model partition to be limited multiple
The junior unit of interconnection, by being defined condition to each junior unit, for example, limits the physical and mechanical parameter of each junior unit, can
With the accurate physico mechanical characteristic for obtaining square each region of body Model.
Step 102：According to the mechanical model of reservoir, centered on solution cavity, the ground for calculating the reservoir being distributed around solution cavity should
The field of force；
Rock by factors such as underground heat, rotation, gravitation, extruding because being influenceed, and its original shape to be recovered outwardly produced
The active force of state, referred to as crustal stress.Reservoir rock is acted on by crustal stress, in the certain area centered on reservoir solution cavity hole
Form stress field.Stress field is relevant with geological form, for example the weight of overlying material, depth where from earth's surface to reservoir,
Variety classes rock is distributed with different depth, and its rock density and thickness are different, therefore affects the crustal stress of different depth
.Meanwhile, solution cavity is of different sizes, and the stress field on its periphery is also different.
The stress field of solution cavity peripheral rock changes with the distance apart from solution cavity hole, and solution cavity is chosen on solution cavity periphery
Country rock in the range of 5 times of axial length, the physical and mechanical parameter of each junior unit is limited using finite element method, is calculated and is surrounded solution cavity
The stress field of the reservoir of distribution.
Square body Model side is acted on along x, y-axis direction stress field, is distributed by the stress field for calculating solution cavity periphery,
It can be seen that there is significant difference in the stress field for the reservoir being distributed around different size solution cavity.
Step 103：Base area stress field calculation obtains the relation between depth of reservoirs and reservoir distortion energy density；
When rock is acted on by external force, its inner product accumulation of energy amount has stored certain deformation energy, the distortion of meaning can density.
Due to the difference of the stress field of reservoir being distributed around different size solution cavity, the distortion energy density of correspondence solution cavity periphery each point is not
Together.Calculating its corresponding reservoir distortion according to the stress field for the reservoir being distributed around solution cavity can density.Obtain around solution cavity point
The corresponding depth of reservoirs of stress field of the reservoir of cloth, can density while being distorted according to reservoir corresponding with the stress field of reservoir
Numerical value, determines the relation between depth of reservoirs and reservoir distortion energy density.
Step 104：According to the relation between depth of reservoirs and reservoir distortion energy density, and the distortion of measuring can be close
Relation between degree and drill bit reaction torque, determines the theoretical values between drill bit reaction torque and depth of reservoirs；
Drill bit reaction torque be drill bit in reservoir in drilling process, bit contact rock and the active force that is subject to.Storage
Layer distortion energy density is big, and reservoir is more broken, then the reaction torque that drill bit is subject to is bigger.On the contrary, reservoir distortion energy density is small, reservoir
More non-breakable, then the reaction torque that drill bit is subject to is smaller.According to reservoir distort can be between density and drill bit reaction torque relation, with
And the relation between depth of reservoirs and reservoir distortion energy density, it can calculate between depth of reservoirs and reservoir distortion energy density
Relation.Different size solution cavities are in identical depth of reservoirs, and its corresponding drill bit reaction torque numerical value is different, obtains drill bit reaction torque
Theoretical values.
Step 105：The measured value between the drill bit reaction torque and depth of reservoirs of target carbonate rock is obtained, number will be surveyed
Value is compared with theoretical values, determines the corresponding mechanical model of target carbonate rock, obtains the corresponding solution cavity chi of mechanical model
Very little size.
In the reservoir actual well drilled containing target solution cavity, strain monitoring element is installed on drill bit.Strain monitoring element
Torque signal is changed into electric signal, and reaches ground.Monitor and record in real time reality of the drill bit reaction torque in different reservoir depth
Survey numerical value.
Drill bit reaction torque is deep in the measured value reservoir identical with the mechanical model pre-established of a certain depth of reservoirs
The theoretical values contrast of corresponding reaction torque is spent, the measured value and the reaction torque in any mechanical model of reaction torque is determined
Theoretical values are identical, you can determine the size of target solution cavity.
The mechanical model that the present embodiment passes through reservoir of the foundation comprising pre-set dimension solution cavity；According to the mechanics mould of the reservoir
Type, centered on the solution cavity, calculates the stress field for the reservoir being distributed around the solution cavity；Calculated according to the stress field
Relation between depth of reservoirs and reservoir distortion energy density；According to the pass between the depth of reservoirs and reservoir distortion energy density
Relation between system, and the distortion energy density and drill bit reaction torque of measuring, determines the drill bit reaction torque and the storage
Theoretical values between layer depth；Obtain target carbonate rock drill bit reaction torque measured value, by the measured value with
The theoretical values are compared, and determine the corresponding mechanical model of the target carbonate rock, determine solution cavity size.This reality
Applying example can accurate solution cavity size in predicting reservoir.
Embodiment two
Fig. 2 is a kind of solution cavity size that the embodiment of the present invention two is provided with the schematic flow sheet for boring recognition methods, such as Fig. 2
Shown, the method for the present embodiment can include：
Step 201：According to multiple default solution cavity sizes, the physics of the corresponding reservoir of various sizes of solution cavity is determined
Mechanics parameter；
The mechanical model of reservoir can be the geometrical model of various shapes, for example, can include different size solution cavity to be internal
Square body Model.The solution cavity included in square body Model is elliposoidal, and different size of solution cavity has different long axis lengths
And minor axis length, according to the combination of multiple default long axis lengths and minor axis length, solution cavity size is determined, and then to every kind of
The solution cavity of the combination of long axis length and minor axis length sets up the mechanical model of its reservoir.Meanwhile, comprising different size solution cavity just
The length of side of cube model is identical.
X, y-axis are horizontal direction axle in square body Model, and z-axis is vertical axes, the geometrical model and border of solution cavity and country rock
Condition, square body Model side is along x, y-axis direction by the squeezing action of surrounding rock, and square bottom surface is along z-axis negative direction by weight
Power is acted on；Square body Model top surface is by overburden pressure.Therefore five faces constraint of the square body Model in addition to top surface is square
Body Model is along the displacement on three directions of x, y, z.
Optionally, according to the long axis length and minor axis length of solution cavity, the thing of the corresponding reservoir of various sizes of solution cavity is determined
Manage mechanics parameter.For example determine the modulus of elasticity of the corresponding reservoir of various sizes of solution cavity, Poisson's ratio, internal friction angle, cohesion,
Tensile strength.
Step 202：According to the physical and mechanical parameter, the mechanical model of the corresponding reservoir of solution cavity of the size is set up；
Using finite element method on square body Model grid division, by square model partition for limited multiple interconnections
Junior unit, by being defined condition to each junior unit, for example, limits the physical and mechanical parameter of each junior unit.Each small list
The physical and mechanical parameter of member includes：Modulus of elasticity, Poisson's ratio, internal friction angle, cohesion, tensile strength.
By using finite element method on square body Model grid division, can accurately set up each region of square body Model
Physical and mechanical parameter.
Step 203：According to the mechanical model of reservoir, centered on solution cavity, calculate and surround the vertical of the reservoir that solution cavity is distributed
The horizontal maximum principal stress of crustal stress, the horizontal minimum principal stress of reservoir and reservoir；
Rock by factors such as underground heat, rotation, gravitation, extruding because being influenceed, and its original shape to be recovered outwardly produced
The active force of state, referred to as crustal stress.Reservoir rock is acted on by crustal stress, in the certain area centered on reservoir solution cavity hole
Form stress field.Stress field is relevant with geological form, for example the weight of overlying material, depth where from earth's surface to reservoir,
Variety classes rock is distributed with different depth, and its rock density and thickness are different, therefore affects the crustal stress of different depth
.Meanwhile, solution cavity is of different sizes, and the stress field on its periphery is also different.
The stress field of solution cavity peripheral rock changes with the distance apart from solution cavity hole, and solution cavity is chosen on solution cavity periphery
Country rock in the range of 5 times of axial length, the physical and mechanical parameter of each junior unit is limited using finite element method, is calculated and is surrounded solution cavity
The stress field of the reservoir of distribution.
Square body Model side is acted on along x, y-axis direction stress field, is distributed by the stress field for calculating solution cavity periphery,
It can be seen that there is significant difference in the stress field for the reservoir being distributed around different size solution cavity.
Optionally, the stress field for the reservoir being distributed around different size solution cavity includes vertical crustal stress, the reservoir of reservoir
Horizontal minimum principal stress and reservoir horizontal maximum principal stress.According to the physical and mechanical parameter of the mechanical model of reservoir, and
Formula one, formula two and formula three, can calculate the vertical crustal stress of reservoir, the horizontal minimum principal stress of reservoir and reservoir respectively
Horizontal maximum principal stress.
Formula one：
Wherein, σvFor the vertical crustal stress of reservoir；ρiFor stratigraphic rock density；Represent to pass through a variety of rocks from earth's surface to reservoir
Layer, the density of every kind of rock stratum is different；G is acceleration of gravity, takes 9.8m/s2；hiFor formation rock thickness, represent from earth's surface to storage
Layer passes through a variety of rock stratum, and the thickness of every kind of rock stratum is different.
Formula two：
Wherein, σHFor the horizontal maximum principal stress of reservoir；σvFor the vertical crustal stress of reservoir；E is the modulus of elasticity of reservoir；μ
For the Poisson's ratio of reservoir；PpFor overlying formation pressure；α is Biot coefficients；εH、εhFor tectonic strain coefficient.
Formula three：
Wherein, σhFor the horizontal minimum principal stress of reservoir；σvFor the vertical crustal stress of reservoir；E is the modulus of elasticity of reservoir；μ
For the Poisson's ratio of reservoir；PpFor overlying formation pressure；α is Biot coefficients；εH、εhFor tectonic strain coefficient.
The vertical crustal stress for the reservoir being distributed by the calculating of formula one, formula two and formula three around different size solution cavity,
The horizontal minimum principal stress of reservoir and the horizontal maximum principal stress size of reservoir are accurate to obtain what is be distributed around different size solution cavity
Stress field cloud atlas, intuitively obtains the stress field size of each region reservoir.Different size solution cavities in identical depth of reservoirs, its
The vertical crustal stress of corresponding reservoir, the horizontal minimum principal stress of reservoir are different with the horizontal maximum principal stress of reservoir.
Step 204：Preset position on well track is obtained, depth of the preset position in reservoir is determined；
According to well track, the calculating track of reservoir distortion energy density is determined.Well is drill bit from earth drilling reservoir institute
The track of process.On the direction that well track is moved closer to solution cavity, choose multiple distortion energy density and calculate point as default
Point position.
Optionally, from 10m positions are pushed up away from solution cavity hole to calculating points away from choosing multiple distortion energy density the 2m positions of hole top
It is used as preset position.
By flexibly designing and choosing preset position, the distortion energy density of different reservoir depth can be obtained, makes meter
The well track covered during distortion energy density is calculated more, so that result of calculation is more accurate.
Step 205：Answered according to the level of the vertical crustal stress of reservoir, the horizontal minimum principal stress of reservoir and reservoir maximum is main
Power, calculates the distortion energy density at the position of preset；
When rock is acted on by external force, its inner product accumulation of energy amount has stored certain deformation energy, the distortion of meaning can density.
Due to the difference of the stress field of reservoir being distributed around different size solution cavity, the distortion energy density of correspondence solution cavity periphery each point is not
Together.Calculating its corresponding reservoir distortion according to the stress field for the reservoir being distributed around solution cavity can density.
Distortion energy density at the position of preset depends on the vertical crustal stress at the position of preset, at the position of preset
Corresponding physical and mechanical parameter at horizontal maximum principal stress and preset position at horizontal minimum principal stress, preset position.
Distortion energy density at the position of preset is calculated according to formula four.
Formula four：
Wherein, UdCan density for distortion；σHFor the horizontal maximum principal stress of reservoir；σhFor the horizontal minimum principal stress of reservoir；
σvFor the vertical crustal stress of reservoir；E is the modulus of elasticity of reservoir；μ is the Poisson's ratio of reservoir.
Step 206：By corresponding depth of reservoirs at the position of preset with distortion can density be combined, obtain depth of reservoirs with
Relation between reservoir distortion energy density.
According to depth of the preset position in reservoir, and reservoir corresponding with the stress field of reservoir distortion can density
Numerical value, determines the relation between depth of reservoirs and reservoir distortion energy density.Different size solution cavities are in identical depth of reservoirs, and its is right
The reservoir distortion answered can density difference.
Step 207：According to the relation between depth of reservoirs and reservoir distortion energy density, and the distortion of measuring can be close
Relation between degree and drill bit reaction torque, determines the theoretical values between drill bit reaction torque and depth of reservoirs；
Drill bit reaction torque be drill bit in reservoir in drilling process, bit contact rock and the active force that is subject to.Storage
Layer distortion energy density is big, and reservoir is more broken, then the reaction torque that drill bit is subject to is bigger.On the contrary, reservoir distortion energy density is small, reservoir
More non-breakable, then the reaction torque that drill bit is subject to is smaller.
The quantitative relationship of drill bit reaction torque and reservoir distortion energy density can be determined in several ways.For example, may be used
With according to the principle of similitude, by scaled down mechanical model, experiment is determined under the conditions of different stress fields, the reaction torque of drill bit
The relation for the energy density that distorted with reservoir.
Experiment determines the quantitative relationship of drill bit reaction torque and reservoir distortion energy density, specifically includes：
According to the mechanical model of reservoir, the rock block of selection scaled down, which is set up in sillar model, sillar model, to be included
The solution cavity of scaled down.The multiple drill bit reaction torques of selection on the solution cavity periphery of scaled down calculate point, are bored using actual
Multiple drill bit reaction torques on well device drilling sillar model calculate point, obtain the corresponding drill bit of different depth in sillar model anti-
Moment of torsion.Meanwhile, the multiple drill bit reaction torques chosen in the mechanical model of reservoir in the sillar model with scaled down are calculated
The corresponding multiple reservoirs distortion of point can density calculate points, according to multiple reservoirs distort can density calculate the corresponding physical mechanics of point and join
Number calculate the distortion of multiple reservoirs can density calculate the stress field of point, according to multiple reservoirs distort can density calculate the crustal stress of point
The distortion of field computation multiple reservoirs can density calculate point reservoir distortion can density, distorted energy with obtaining the reservoir of different reservoir depth
Density.Finally, it is anti-according to the drill bit of correspondence depth in the reservoir distortion energy density and sillar model of different depth in mechanical model
Moment of torsion, fitting obtains the quantitative relationship of drill bit reaction torque and reservoir distortion energy density.
The quantitative relationship of the drill bit reaction torque being determined by experiment and distortion energy density, as shown in formula five：
Formula five：Mbt=kUd+c；
Wherein, MbtFor drill bit reaction torque；UdBeing distorted for reservoir can density；K, c are experiment parameter；Specifically, k value exists
0.7 or so, c value are between -1~1.
According to the quantitative relationship between reservoir distortion energy density and drill bit reaction torque, and depth of reservoirs and reservoir distortion energy
Relation between density, can calculate the relation between depth of reservoirs and reservoir distortion energy density.Different size solution cavities are in phase
Same depth of reservoirs, its corresponding drill bit reaction torque numerical value is different, thus obtains the reason of the drill bit reaction torque of different reservoir depth
By numerical value.
The quantitative relationship of drill bit reaction torque and reservoir distortion energy density is determined by experiment, makes the relationship of the two more accurate,
Make the theoretical values of drill bit reaction torque corresponding with depth of reservoirs calculated more accurate.
Step 208：The measured value between the drill bit reaction torque and depth of reservoirs of target carbonate rock is obtained, number will be surveyed
Value is compared with theoretical values, determines the corresponding mechanical model of target carbonate rock, obtains the corresponding solution cavity chi of mechanical model
Very little size.
In the reservoir actual well drilled containing target solution cavity, strain monitoring element is installed on drill bit.Strain monitoring element
Torque signal is changed into electric signal, and reaches ground.Monitor and record in real time reality of the drill bit reaction torque in different reservoir depth
Survey numerical value.
Drill bit reaction torque is deep in the measured value reservoir identical with the mechanical model pre-established of a certain depth of reservoirs
The theoretical values contrast of corresponding reaction torque is spent, the measured value and the reaction torque in any mechanical model of reaction torque is determined
Theoretical values are identical, you can determine the size of target solution cavity.
The mechanical model that the present embodiment passes through reservoir of the foundation comprising pre-set dimension solution cavity；According to the mechanics mould of the reservoir
Type, centered on the solution cavity, calculates the stress field for the reservoir being distributed around the solution cavity；Calculated according to the stress field
Relation between depth of reservoirs and reservoir distortion energy density；According to the pass between the depth of reservoirs and reservoir distortion energy density
Relation between system, and the distortion energy density and drill bit reaction torque of measuring, determines the drill bit reaction torque and the storage
Theoretical values between layer depth；Obtain target carbonate rock drill bit reaction torque measured value, by the measured value with
The theoretical values are compared, and determine the corresponding mechanical model of the target carbonate rock, determine solution cavity size.This reality
Applying example can accurate solution cavity size in predicting reservoir.
Embodiment three
Fig. 3 be the embodiment of the present invention three provide a kind of solution cavity size with bore identifying device schematic diagram, as shown in figure 3,
The device of the present embodiment can include：Build module 31, the first computing module 32, the second computing module 33, the 3rd computing module
34th, identification module 35；
Wherein, module 31 is built, for according to multiple default solution cavity sizes, it to be set up to the solution cavity of every kind of size
The mechanical model of reservoir；
First computing module 32, for the mechanical model according to reservoir, centered on solution cavity, calculates what is be distributed around solution cavity
The stress field of reservoir；
Second computing module 33, is obtained between depth of reservoirs and reservoir distortion energy density for base area stress field calculation
Relation；
3rd computing module 34, for being surveyed according to the relation between depth of reservoirs and reservoir distortion energy density, and experiment
Relation between fixed distortion energy density and drill bit reaction torque, determines the theoretical values between drill bit reaction torque and depth of reservoirs；
Identification module 35, for obtaining the measured value between the drill bit reaction torque of target carbonate rock and depth of reservoirs,
Measured value is compared with theoretical values, the corresponding mechanical model of target carbonate rock is determined, mechanical model correspondence is obtained
Solution cavity size.
On the device in the present embodiment, wherein modules perform the concrete mode of operation in relevant this method
It is described in detail in embodiment, explanation will be not set forth in detail herein.
The mechanical model that the present embodiment passes through reservoir of the foundation comprising pre-set dimension solution cavity；According to the mechanical model of reservoir,
Centered on solution cavity, the stress field for the reservoir being distributed around solution cavity is calculated；Base area stress field calculation depth of reservoirs and reservoir
Relation between distortion energy density；Relation between the energy density that distorted according to depth of reservoirs and reservoir, and measuring are abnormal
The relation become between energy density and drill bit reaction torque, determines the theoretical values between drill bit reaction torque and depth of reservoirs；Obtain mesh
The measured value of the drill bit reaction torque of carbonate rock is marked, measured value is compared with theoretical values, determines that target solution cavity is stored up
The corresponding mechanical model of layer, determines solution cavity size.The present embodiment can accurate solution cavity size in predicting reservoir.
Example IV：
Fig. 4 be the embodiment of the present invention four provide a kind of solution cavity size with bore identifying device schematic diagram, as shown in figure 4,
The device of the present embodiment can include：Build module 41, the first computing module 42, the second computing module 43, the 3rd computing module
44th, identification module 45；
Wherein, module 41 is built, for according to multiple default solution cavity sizes, it to be set up to the solution cavity of every kind of size
The mechanical model of reservoir；
First computing module 42, for the mechanical model according to reservoir, centered on solution cavity, calculates what is be distributed around solution cavity
The stress field of reservoir；
Second computing module 43, is obtained between depth of reservoirs and reservoir distortion energy density for base area stress field calculation
Relation；
3rd computing module 44, for being surveyed according to the relation between depth of reservoirs and reservoir distortion energy density, and experiment
Relation between fixed distortion energy density and drill bit reaction torque, determines the theoretical values between drill bit reaction torque and depth of reservoirs；
Identification module 45, for obtaining the measured value between the drill bit reaction torque of target carbonate rock and depth of reservoirs,
Measured value is compared with theoretical values, the corresponding mechanical model of target carbonate rock is determined, mechanical model correspondence is obtained
Solution cavity size.
Optionally, module 41 is built to specifically include：Parameter determination unit 411, processing unit 412；
Wherein, parameter determination unit 411, for the size according to solution cavity, determine the thing of the corresponding reservoir of solution cavity of size
Manage mechanics parameter.
Specifically, physical and mechanical parameter is specifically included：It is the modulus of elasticity of reservoir, Poisson's ratio, internal friction angle, cohesion, anti-
Tensile strength.
Processing unit 412, for according to physical and mechanical parameter, setting up the mechanical model of the corresponding reservoir of solution cavity of size.
First computing module 42 is specifically included：First computing unit 421, the second computing unit 422, the 3rd computing unit
423；
Wherein, the first computing unit 421, the vertical crustal stress for calculating reservoir；
Second computing unit 422, the horizontal minimum principal stress for calculating reservoir；
3rd computing unit 423, the horizontal maximum principal stress for calculating reservoir.
Second computing module 43 is specifically included：First determining unit 431, the second determining unit 432, the 3rd determining unit
433；
Wherein, the first determining unit 431, for obtaining preset position on well track, and determines that preset position exists
Depth in reservoir；
Second determining unit 432, for according to stress field, calculating the distortion energy density at the position of preset；
3rd determining unit 433, for corresponding depth of reservoirs at the position of preset to be combined with distortion energy density, is obtained
Relation between depth of reservoirs and reservoir distortion energy density.
On the device in the present embodiment, wherein modules perform the concrete mode of operation in relevant this method
It is described in detail in embodiment, explanation will be not set forth in detail herein.
The mechanical model that the present embodiment passes through reservoir of the foundation comprising pre-set dimension solution cavity；According to the mechanical model of reservoir,
Centered on solution cavity, the stress field for the reservoir being distributed around solution cavity is calculated；Base area stress field calculation depth of reservoirs and reservoir
Relation between distortion energy density；Relation between the energy density that distorted according to depth of reservoirs and reservoir, and measuring are abnormal
The relation become between energy density and drill bit reaction torque, determines the theoretical values between drill bit reaction torque and depth of reservoirs；Obtain mesh
The measured value of the drill bit reaction torque of carbonate rock is marked, measured value is compared with theoretical values, determines that target solution cavity is stored up
The corresponding mechanical model of layer, determines solution cavity size.The present embodiment can accurate solution cavity size in predicting reservoir.
One of ordinary skill in the art will appreciate that：Realizing all or part of step of above method embodiment can pass through
Programmed instruction related hardware is completed, and foregoing program can be stored in a computer read/write memory medium, the program
Upon execution, the step of including above method embodiment is performed；And foregoing storage medium includes：ROM, RAM, magnetic disc or light
Disk etc. is various can be with the medium of store program codes.
Finally it should be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
The present invention is described in detail with reference to the foregoing embodiments, it will be understood by those within the art that：It still may be used
To be modified to the technical scheme described in foregoing embodiments, or equivalent substitution is carried out to which part technical characteristic；
And these modifications or replacement, the essence of appropriate technical solution is departed from the scope of various embodiments of the present invention technical scheme.
Claims (10)
1. a kind of solution cavity size with bore recognition methods, it is characterised in that including：
According to multiple default solution cavity sizes, the mechanical model of its reservoir is set up to the solution cavity of every kind of size；
According to the mechanical model of the reservoir, centered on the solution cavity, the ground for calculating the reservoir being distributed around the solution cavity should
The field of force；
Relation between depth of reservoirs and reservoir distortion energy density is obtained according to stress field calculating；
According to the distortion energy density and drill bit of the relation between the depth of reservoirs and reservoir distortion energy density, and measuring
Relation between reaction torque, determines the theoretical values between the drill bit reaction torque and the depth of reservoirs；
Obtain target carbonate rock drill bit reaction torque and depth of reservoirs between measured value, by the measured value with it is described
Theoretical values are compared, and determine the corresponding mechanical model of the target carbonate rock, obtain the mechanical model corresponding molten
Hole size.
2. according to the method described in claim 1, it is characterised in that the solution cavity to every kind of size sets up the mechanics of its reservoir
Model, including：
According to the size of solution cavity, the physical and mechanical parameter of the corresponding reservoir of solution cavity of the size is determined；
According to the physical and mechanical parameter, the mechanical model of the corresponding reservoir of solution cavity of the size is set up.
3. according to the method described in claim 1, it is characterised in that the stress field includes：
The horizontal maximum principal stress of the vertical crustal stress of reservoir, the horizontal minimum principal stress of reservoir and reservoir.
4. according to the method described in claim 1, it is characterised in that described calculated according to the stress field obtains depth of reservoirs
Relation between the energy density that distorted with reservoir, including：
Preset position on well track is obtained, depth of the preset position in the reservoir is determined；
According to the stress field, the distortion energy density at the preset position is calculated；
Corresponding depth of reservoirs at the preset position is combined with distortion energy density, depth of reservoirs is obtained and is distorted with reservoir
Relation between energy density.
5. method according to claim 2, it is characterised in that the size according to solution cavity, determines the molten of the size
The physical and mechanical parameter of the corresponding reservoir in hole, including：
Determine the modulus of elasticity of the corresponding reservoir of solution cavity of the size, Poisson's ratio, internal friction angle, cohesion, tensile strength.
6. a kind of solution cavity size with bore identifying device, it is characterised in that including：Build module, the first computing module, the second meter
Calculate module, the 3rd computing module, identification module；
Wherein, the structure module, for according to multiple default solution cavity sizes, its storage to be set up to the solution cavity of every kind of size
The mechanical model of layer；
First computing module, for the mechanical model according to the reservoir, centered on the solution cavity, is calculated around described
The stress field of the reservoir of solution cavity distribution；
Second computing module, is obtained between depth of reservoirs and reservoir distortion energy density for being calculated according to the stress field
Relation；
3rd computing module, for according to the relation between the depth of reservoirs and reservoir distortion energy density, and experiment
Relation between the distortion energy density and drill bit reaction torque of measure, is determined between the drill bit reaction torque and the depth of reservoirs
Theoretical values；
The identification module, will for obtaining the measured value between the drill bit reaction torque of target carbonate rock and depth of reservoirs
The measured value is compared with the theoretical values, determines the corresponding mechanical model of the target carbonate rock, obtains institute
State the corresponding solution cavity size of mechanical model.
7. device according to claim 6, it is characterised in that the structure module includes：Parameter determination unit, processing are single
Member；
Wherein, the parameter determination unit, for the size according to solution cavity, determines the thing of the corresponding reservoir of solution cavity of the size
Manage mechanics parameter；
The processing unit, for according to the physical and mechanical parameter, setting up the mechanics of the corresponding reservoir of solution cavity of the size
Model.
8. device according to claim 6, it is characterised in that first computing module includes：First computing unit,
Two computing units, the 3rd computing unit；
Wherein, first computing unit, the vertical crustal stress for calculating reservoir；
Second computing unit, the horizontal minimum principal stress for calculating reservoir；
3rd computing unit, the horizontal maximum principal stress for calculating reservoir.
9. device according to claim 6, it is characterised in that second computing module includes：First determining unit,
Two determining units, the 3rd determining unit；
Wherein, first determining unit, for obtaining preset position on well track, and determines that the preset position exists
Depth in the reservoir；
Second determining unit, for according to the stress field, calculating the distortion energy density at the preset position；
3rd determining unit, for corresponding depth of reservoirs at the preset position to be combined with distortion energy density,
Obtain the relation between depth of reservoirs and reservoir distortion energy density.
10. device according to claim 7, it is characterised in that the parameter determination unit, specifically for：According to solution cavity
Size, determine that the modulus of elasticity of the corresponding reservoir of solution cavity of the size, Poisson's ratio, internal friction angle, cohesion, tension are strong
Degree.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710326895.1A CN107169192B (en) | 2017-05-10 | 2017-05-10 | Method and device for identifying size of karst cave while drilling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710326895.1A CN107169192B (en) | 2017-05-10 | 2017-05-10 | Method and device for identifying size of karst cave while drilling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107169192A true CN107169192A (en) | 2017-09-15 |
| CN107169192B CN107169192B (en) | 2020-06-09 |
Family
ID=59812823
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710326895.1A Active CN107169192B (en) | 2017-05-10 | 2017-05-10 | Method and device for identifying size of karst cave while drilling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107169192B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009127885A2 (en) * | 2008-04-17 | 2009-10-22 | Darko Horvat | Flexible crankshaft for software controllable aircraft wheel drive |
| CN102373923A (en) * | 2010-08-20 | 2012-03-14 | 中国石油天然气股份有限公司 | Reservoir stratum identification method |
| CN103278614A (en) * | 2013-04-28 | 2013-09-04 | 北京大学 | Method and device for correcting dynamic and static rock mechanical parameters |
| CN205858201U (en) * | 2016-07-27 | 2017-01-04 | 中国石油化工股份有限公司 | The reaction torque of a kind of mud motor controls joint |
| CN106446396A (en) * | 2016-09-20 | 2017-02-22 | 中国石油大学(华东) | Method and device for determining influences of distributions of fractures and karst caves on reservoir permeability |
-
2017
- 2017-05-10 CN CN201710326895.1A patent/CN107169192B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009127885A2 (en) * | 2008-04-17 | 2009-10-22 | Darko Horvat | Flexible crankshaft for software controllable aircraft wheel drive |
| CN102373923A (en) * | 2010-08-20 | 2012-03-14 | 中国石油天然气股份有限公司 | Reservoir stratum identification method |
| CN103278614A (en) * | 2013-04-28 | 2013-09-04 | 北京大学 | Method and device for correcting dynamic and static rock mechanical parameters |
| CN205858201U (en) * | 2016-07-27 | 2017-01-04 | 中国石油化工股份有限公司 | The reaction torque of a kind of mud motor controls joint |
| CN106446396A (en) * | 2016-09-20 | 2017-02-22 | 中国石油大学(华东) | Method and device for determining influences of distributions of fractures and karst caves on reservoir permeability |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107169192B (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106814407B (en) | The method for determining Three-dimensional Rock crustal stress based on single borehole deformation measurement | |
| CN105672971B (en) | A kind of Reservoir Fracture Opening pressure, open order and water injection pressure Forecasting Methodology | |
| CN108680952A (en) | A kind of strike-slip fault structural evolution analytic method | |
| US20170212276A1 (en) | Method for characterizing and exploiting a subterranean formation comprising a network of fractures | |
| MX2008015872A (en) | Method and system for designing and optimizing drilling and completion operations in hydrocarbon reservoirs. | |
| CN105974056A (en) | Tunnel water bursting disaster precursor information monitoring model test system and test method thereof | |
| CN101421640A (en) | Methods and apparatus for predicting the hydrocarbon production of a well location | |
| CN102253415A (en) | Method for establishing earthquake response mode based on fracture equivalent medium model | |
| CN108572401B (en) | Construction method of fracture-cavity combined model and method for detecting deformation of reservoir fracture-cavity | |
| MX2013014709A (en) | Systems and methods for measuring parameters of a formation. | |
| CN108842821B (en) | Calculation method for reasonable buried depth of submarine tunnel constructed by drilling and blasting method | |
| WO2017031615A1 (en) | Method of constructing ground subsidence prediction model of two-cavity salt cavern reservoir | |
| CN107742020A (en) | Shale gas reservoir fracturing improvement method for predicting volume | |
| CN103389489A (en) | Micro earthquake monitoring and positioning method based on highly-deviated well | |
| CN116703244B (en) | Mining subsidence area treatment effect and comprehensive evaluation method | |
| CN113779811A (en) | A Quantitative Analysis Method of Tectonic Deformation to Present Orientation Disturbance of In-situ Stress | |
| Hekmatnejad et al. | A hybrid predictive model of unstable rock blocks around a tunnel based on estimated volumetric fracture intensity and circular variance from borehole data sets | |
| CN103390108B (en) | A kind of well type method for designing based on fracturing geology body compressibility and device | |
| CN106199712B (en) | A kind of method and device in definite pressure break casing deformation region | |
| CN114547927A (en) | Hydrothermal deposit three-dimensional quantitative prediction evaluation method based on numerical simulation | |
| CN107169192A (en) | Method and device for identifying size of karst cave while drilling | |
| CN109655400B (en) | Roadway plastic region analysis method | |
| CN114200538B (en) | Ground stress direction determining method, device, storage medium and computer equipment | |
| Leary et al. | Coupling crustal seismicity to crustal permeability-Power-law spatial correlation for EGS-induced and hydrothermal seismicity | |
| Grasselli et al. | Rock fabric not principal stress dictates SRV: The Story of how a~ 70 Year-old discounted data point still plagues our industry and how true triaxial testing finally confirms it |
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 |