CN106150484A - A kind of drilling well leakage Forecasting Methodology based on numeral rock mass - Google Patents
A kind of drilling well leakage Forecasting Methodology based on numeral rock mass Download PDFInfo
- Publication number
- CN106150484A CN106150484A CN201510163437.1A CN201510163437A CN106150484A CN 106150484 A CN106150484 A CN 106150484A CN 201510163437 A CN201510163437 A CN 201510163437A CN 106150484 A CN106150484 A CN 106150484A
- Authority
- CN
- China
- Prior art keywords
- density
- parameter
- credibility
- pressure
- drilling
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a kind of drilling well leakage Forecasting Methodology based on numeral rock mass, the method includes: gather the borehole track data of current well, to obtain the coordinate at borehole track each position;According to the coordinate at each position in borehole track, from the digital rock mass of current well, obtain the geologic parameter information at each position;According to the coordinate at each position in borehole track, obtain the engineering parameter information at each position;Drilling well leakage value-at-risk according to the geologic parameter information at each position and this position of engineering parameter acquisition of information.The inventive method can be based on relevant geologic parameter and drilling engineering parameter, well section complete to drilling well carries out leakage risk profile, both can be used to carry out evaluation and the optimization of Drilling Design scheme, the leakage symptom of a trend can be identified again during wellbore construction in time, thus avoid the generation of leakage to greatest extent.
Description
Technical field
The invention belongs to oil gas well drilling risk profile technical field, particularly relate to a kind of brill based on numeral rock mass
Well leakage Forecasting Methodology.
Background technology
Leakage is one of most common complexity event in deep & ultra-deep well drilling engineering.In order to avoid there is leakage
And reservoir protec-tion, it is necessary to the probability of the full well of look-ahead each interval generation leakage risk, carry out drilling well in advance
The optimization of scheme or take preventive measures, pays close attention to excessive risk well section in work progress, or identifies in time i.e.
The leakage risk that will occur, shoots the arrow at the target, it is ensured that drilling safety.
For solving an above-mentioned difficult problem, domestic and international experts and scholars have carried out some research, the Forecasting recognition of existing leakage
Method, or observe on-the-spot phenomenon, or the various characteristics according to leakage after leakage occurs by sensor
Determine, or the geologic feature causing leakage to occur carries out fuzzy evaluation macroscopic view obtains the potential well in somewhere
Leak out danger, or utilize the side that the partial parameters of well to be predicted matches with some " eigenvalues " of offset well leakage
Method predicts leakage.
But, most prior art is the method afterwards judged, although can control risk, and retrieves part and damages
Lose, but so inevitably lead to drilling fluid leakage and hydrocarbon contamination.And the method predicted in advance does not has
Consider the interaction relationship between concrete drilling engineering design information and geological condition, therefore cannot be for tool
The prevention that specific aim measure is optimized taked by the well of body.
Summary of the invention
One of the technical problem to be solved is to need to provide a kind of drilling well leakage based on numeral rock mass
Forecasting Methodology, the method can be in conjunction with geology parameter information and engineering parameter information well complete to the drilling well potential leakage of section
Risk is predicted.
In order to solve above-mentioned technical problem, embodiments herein provides a kind of drilling well well based on numeral rock mass
Leakage Forecasting Methodology, the method includes: gather the borehole track data of current well, with obtain borehole track each
Put the coordinate at place;According to the coordinate at each position in borehole track, obtain every from the digital rock mass of current well
The geologic parameter information of one position;According to the coordinate at each position in borehole track, obtain at each position
Engineering parameter information;According to the geologic parameter information at each position and this position of engineering parameter acquisition of information
Drilling well leakage value-at-risk.
Preferably, in the step of the geologic parameter information obtained from the digital rock mass of current well at each position
In, including: from the digital rock mass of current well, extract the numerical value of geologic parameter at each position;According to numeral
The Data Source of rock mass and building process, the confidence value of the geologic parameter at acquisition each position, wherein, institute
State geologic parameter include the following: solution cavity information, fault information, crack information, rock type, permeability,
Layer pore pressure equal yield density and formation fracture pressure equal yield density.
Preferably, in the step of the engineering parameter information obtained at each position, including: set from drilling engineering
Meter scheme or the numerical value with the engineering parameter at extracting data each position of drilling collection;True according to Data Source
The credibility of fixed every kind of engineering parameter, wherein, described engineering parameter includes the following: drilling fluid density, drilling fluid
Plastic viscosity, discharge capacity, drilling speed, landwaste size, landwaste, borehole diameter, drilling tool internal diameter, drilling tool external diameter, water
Eye diameter and surge pressure coefficient or surge pressure equal yield density;Based on the engineering ginseng at each position extracted
Number information is calculated the equivalent circulating density information at each position.
Preferably, in being calculated the step of the pressure equivalent circulating density information at each position, including:
According to the engineering parameter information at each position, the pressure equivalent selecting respective algorithms to calculate at each position follows
Ring density;The algorithm used according to calculating pressure equivalent circulating density and the parameter participating in calculating obtain pressure and work as
The credibility of amount circulating density.
Preferably, obtain in algorithm and the parameter participating in calculating used according to calculating pressure equivalent circulating density
In the step of the credibility of pressure equivalent circulating density ECD, farther include: that assumes participation calculating has n
Individual parameter E1 ... En, calculate independent credibility CF (Ei) that each parameter calculates for ECD=
CF (E) [i] * CF (ECD algorithm), CF (ECD algorithm) represent the credibility of the algorithm calculating ECD, CF (E) [i]
Represent the credibility of parameter Ei;Synthetic reliability CF (E1E2)=CF (the E1)+CF (E2) of calculating E1 and E2-
CF (E1) * CF (E2), then using CF (E1E2) as independent credibility, is calculated it with CF (E3)
CF (E1E2E3), iterative computation is until last parameter En successively, and then obtains final synthetic reliability
CF (ECD), as the credibility of pressure equivalent circulating density ECD.
Preferably, according to the geologic parameter information at each position and this position of engineering parameter acquisition of information
In the step of drilling well leakage value-at-risk, including: believe according to the geologic parameter information at each position and engineering parameter
Breath judges whether to meet pre-conditioned: if meeting, then determine initial risk values according to pre-conditioned, and calculate should
The credibility of initial risk values, obtains final risk value based on initial risk values and credibility thereof;If being unsatisfactory for pre-
If condition, it is determined that final risk value is 0.
Preferably, the described pre-conditioned at least one of that includes: there is solution cavity or there is tomography;Pressure equivalent
Circulating density is close more than formation fracture pressure equivalent with surge pressure coefficient or surge pressure equal yield density sum
Degree;Crack for grow very much, and pressure equivalent circulating density is close with surge pressure coefficient or surge pressure equivalent
Degree sum is more than formation pore pressure equal yield density;Crack is general growth, and pressure equivalent circulating density is with sharp
Dynamic pressure coefficient or surge pressure equal yield density sum are more than formation pore pressure equal yield density;Rock type is
Sandstone, conglomerate or glutenite, and permeability is more than certain value, pressure equivalent circulating density and surge pressure coefficient
Or surge pressure equal yield density sum is more than formation pore pressure equal yield density;Rock type is sandstone, conglomerate
Or glutenite, and pressure equivalent circulating density is more than with surge pressure coefficient or surge pressure equal yield density sum
Formation pore pressure equal yield density;Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density
More than or equal to formation pore pressure equal yield density;Pressure equivalent circulating density and surge pressure coefficient or exciting pressure
Power equal yield density sum is more than formation pore pressure equal yield density;Pressure equivalent circulating density is more than formation pore pressure
Power equal yield density.
Preferably, according to met pre-conditioned in the credibility of each parameter calculate the credible of initial risk values
Degree.
Preferably, also include: based on geologic parameter information, engineering parameter information and drilling well leakage value-at-risk, with
Leakage risk profile section is shown with curve form on the basis of same well depth.
Preferably, also include: the drilling well leakage value-at-risk in a position is more than or equal to when setting thresholding, to this
The engineering parameter information of position is adjusted, so that the drilling well leakage value-at-risk of this position is less than setting door
Limit.
Compared with prior art, the one or more embodiments in such scheme can have the advantage that or have
Benefit effect.
The drilling well leakage Forecasting Methodology based on numeral rock mass of the present invention, it is possible to based on relevant geologic parameter and brill
Well engineering parameter, well section complete to drilling well carries out leakage risk profile, has not only been embodied leakage severity but also has embodied and has sent out
The quantization risk class of raw probability, both can be used to carry out evaluation and the optimization of Drilling Design scheme, can bore again
Identify the leakage symptom of a trend during well construction in time, thus avoid the generation of leakage to greatest extent.
Other features and advantages of the present invention will illustrate in the following description, and, partly from explanation
Book becomes apparent, or understands by implementing technical scheme.The purpose of the present invention and its
He can come real by structure specifically noted in description, claims and accompanying drawing and/or flow process at advantage
Now and obtain.
Accompanying drawing explanation
Accompanying drawing is used for providing being further appreciated by of the technical scheme to the application or prior art, and constitutes
A part for bright book.Wherein, the accompanying drawing expressing the embodiment of the present application is used for explaining together with embodiments herein
The technical scheme of the application, but it is not intended that the restriction to technical scheme.
Fig. 1 is the schematic flow sheet of the drilling well leakage Forecasting Methodology based on numeral rock mass of the embodiment of the present invention one.
Fig. 2 is the schematic flow sheet of the drilling well leakage Forecasting Methodology based on numeral rock mass of one example of the present invention.
Detailed description of the invention
Embodiments of the present invention are described in detail, whereby to the present invention how below with reference to drawings and Examples
Application technology means solve technical problem, and the process that realizes reaching relevant art effect can fully understand and evidence
To implement.Each feature in the embodiment of the present application and embodiment, can mutually tie under not colliding premise
Closing, the technical scheme formed is all within protection scope of the present invention.
It addition, the step shown in the flow chart of accompanying drawing can be at the computer of such as one group of computer executable instructions
System performs.And, although show logical order in flow charts, but in some cases, permissible
To be different from the step shown or described by order execution herein.
(embodiment one)
Fig. 1 is the schematic flow sheet of the drilling well leakage Forecasting Methodology based on numeral rock mass of the embodiment of the present invention, under
Each step of process in detail is carried out with reference to Fig. 1 in face.
In step s 110, gather the borehole track data of current well, to obtain at borehole track each position
Coordinate.
Borehole track data include the information such as well depth, hole angle, azimuth.Typically, preferably every relatively closely spaced
Take a bit, take 1 point every 1 meter in this example, gather the data of each point successively, then by each point
Data are converted to cartesian coordinate value X, Y, Z by calculating, finally give a well rail about current well
The coordinate sequence in road.
In the step s 120, according to the coordinate at each position in borehole track, from the digital rock mass of current well
Geologic parameter information (including geologic parameter value and confidence value thereof) at middle acquisition each position.
It should be noted that numeral rock mass refers to series three-dimensional mesh data body, each grid node has coordinate
Value and the property value (such as solution cavity, crack, porosity, pore pressure etc.) with rock mass.Typically, based on
Existing geological structure model, it becomes possible to construct the digital rock mass for drilling well.And, from numeral rock mass
Obtain the various engineering geology parameters along well track, it is possible to fine optimization wellbore construction scheme, improve prediction and bore
The precision of well leakage.
Specifically, in the step of the geologic parameter information obtained from the digital rock mass of current well at each position
In, including: from the digital rock mass of current well, extract the numerical value of geologic parameter at each position, according to numeral
The Data Source of rock mass and building process, the confidence value of the geologic parameter at acquisition each position.Wherein, ground
Matter parameter includes the following: solution cavity information, fault information, crack information, rock type, permeability, earth bore
Gap pressure equivalent density and formation fracture pressure equal yield density etc..
In step s 130, according to the coordinate at each position in borehole track, the work at each position is obtained
Journey parameter information (includes engineering parameter values and confidence value thereof).
In particular it is required that execution following steps: from drilling engineering design scheme or with the extracting data of drilling collection
Engineering parameter at each position, determines the credibility of every kind of engineering parameter according to Data Source.Wherein, engineering
Parameter include the following: drilling fluid density, drilling fluid plastic viscosity, discharge capacity, drilling speed, landwaste size, landwaste,
Surge pressure coefficient or the calculating of borehole diameter, drilling tool internal diameter, drilling tool external diameter, hydrophthalmia diameter and setting obtain
Surge pressure equal yield density etc..Then, according to the engineering parameter information at each position extracted, calculate
Pressure equivalent circulating density (being also an engineering parameter) information at each position.
In being calculated the step of the pressure equivalent circulating density information at each position, specifically, according to often
The engineering parameter information of one position, the pressure equivalent circulation selecting respective algorithms to calculate at each position is close
Degree, the algorithm then used according to calculating pressure equivalent circulating density and the parameter participating in calculating obtain pressure and work as
The credibility of amount circulating density.
Obtain pressure equivalent follow at the algorithm used according to calculating pressure equivalent circulating density and the parameter used
In the step of the credibility of ring density ECD, farther include: that assumes participation calculating has n parameter E1 ...
En, the most first calculates independent credibility CF (Ei)=CF (E) [i] * CF (ECD calculation that each parameter calculates for ECD
Method), CF (ECD algorithm) represents the credibility of the algorithm calculating ECD, and CF (E) [i] represents the credible of parameter Ei
Degree.Then synthetic reliability CF (E1E2)=CF (E1)+CF (E2)-CF (E1) the * CF (E2) of E1 and E2 is calculated,
Then CF (E1E2) is calculated CF (E1E2E3) by similar method with CF (E3) as independent credibility,
I.e. CF (E1E2E3)=CF (E1E2)+CF (E3)-CF (E1E2) * CF (E3) ... successively iterative computation until
The position of last parameter En, obtains final synthetic reliability CF (ECD), as pressure equivalent
The credibility of circulating density ECD.
It should be noted that existing method all assumes that the credibility of used data is 100% accurately, institute
The calculating used and determination methods are also 100% reliable, but these data and method unreasonably in real work
Wanting, compared to existing technology, the present invention considers the reliability of data and computational methods, further increases pre-
Survey the degree of accuracy of drilling well leakage.
In step S140, believe according to the geologic parameter information at each position of borehole track and engineering parameter
Breath obtains the drilling well leakage value-at-risk of this position.
Specifically, satisfied presetting is judged whether according to the geologic parameter information at each position and engineering parameter information
Condition, if meeting, then determines initial risk values according to pre-conditioned, calculates the credibility of this initial risk values,
Final risk value is obtained based on initial risk values and credibility thereof.If being unsatisfactory for pre-conditioned, it is determined that final wind
Danger value is 0.In this example, according to met pre-conditioned in the credibility of each parameter calculate initial risks
The credibility of value.
Further, the above-mentioned pre-conditioned at least one of that can include: (1) exists solution cavity or there is tomography;
(2) pressure equivalent circulating density breaks more than stratum with surge pressure coefficient or surge pressure equal yield density sum
Split pressure equal yield density;(3) crack for grow very much, and pressure equivalent circulating density and surge pressure coefficient
Or surge pressure equal yield density sum is more than formation pore pressure equal yield density;(4) crack is general growth,
And pressure equivalent circulating density is more than formation pore with surge pressure coefficient or surge pressure equal yield density sum
Pressure equivalent density;(5) rock type is sandstone, conglomerate or glutenite, and permeability is more than certain value,
Pressure equivalent circulating density is more than formation pore pressure with surge pressure coefficient or surge pressure equal yield density sum
Power equal yield density;(6) rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density is with sharp
Dynamic pressure coefficient or surge pressure equal yield density sum are more than formation pore pressure equal yield density;(7) rock
Type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density is close more than or equal to formation pore pressure equivalent
Degree;(8) pressure equivalent circulating density is more than ground with surge pressure coefficient or surge pressure equal yield density sum
Layer pore pressure equal yield density;(9) pressure equivalent circulating density is more than formation pore pressure equal yield density.
Owing to the present invention not only allows for initial risk values when determining value-at-risk, it is also contemplated that initial risk values can
Reliability, the degree of accuracy of therefore obtained drilling well leakage value-at-risk is higher.Although prior art judges afterwards
Method can be controlled risk as far as possible, retrieves partial loss, but so inevitably leads to drilling fluid leakage and oil
Gas-bearing formation pollute, and if take measures not in time, easily make the state of affairs deteriorate, lead to a disaster.And the present invention exists
The design initial stage just can accurately predict the risk that leakage occurs, it is possible to early avoids the generation of risk.
It addition, as it is shown in figure 1, the method can also include step S150.
In step S150, based on geologic parameter information, engineering parameter information and drilling well leakage value-at-risk, with
Leakage risk profile section is shown with curve form on the basis of same well depth axle.Cut open at this leakage risk profile
Face figure can demonstrating, the geologic parameter curve of the borehole track of this well, engineering parameter curve and risk class are bent
Line.If being provided with alarm threshold, then this leakage risk profile section can also demonstrate that risk alarm is pointed out
Information.
When drilling well leakage value-at-risk in a position is more than or equal to setting thresholding, the engineering parameter to this position
Information is adjusted, so that the drilling well leakage value-at-risk of this position is less than setting thresholding.
Support technician can check geology, project situation and the correspondence of full well section by this step comprehensively
Risk class, the well section that especially risk class is high, by analyzing reason further, formulate improve or reply
Measure, it is possible to the generation controlled risk.
Fig. 2 is the schematic flow sheet of the drilling well leakage Forecasting Methodology based on numeral rock mass of one example of the present invention.
This example carries out leakage risk profile to Drilling Design scheme, describes acquisition drilling well in detail below with reference to Fig. 2
The method of leakage value-at-risk.
In step 1, from the beginning of 0m, with well depth 1m for interval, from well designed path extracting data
The coordinate figure (well depth, north coordinate, east coordinate) of each point, obtains coordinate array T [N] [3] and (assumes altogether
There is N number of coordinate points).
In step 2, the coordinate array obtained according to step 1, extracts each coordinate points from numeral rock mass
Geologic parameter data, form geologic parameter array G [N] [I] (I=7), and wherein i represents solution cavity respectively from 0 to 6 (has
Or nothing), tomography (with or without), crack (grow very much, general grow, nothing), rock type, infiltration
Rate K, formation pore pressure equal yield density rp, formation fracture pressure equal yield density rf.Meanwhile, according to numeral rock
The primary data source of body and building process, credibility CF (G) [n] [i] of the given each attribute of each degree of depth,
To credibility array CF (G) [N] [I] (I=7).
It should be noted that what the historical data that digital rock mass is mainly based upon offset well built, therefore offset well number
Amount, offset well and the average distance of this well, offset well the data precision are all by the credibility of impact numeral rock mass.
Typically, the independent credibility of three can be set according to expertise and certain rule: 1. build numeral rock
The credibility of the offset well the data precision of body (by expert consider each item data of these offset wells original accurately
Degree, the Maturity of computational methods and give), the such as credibility of this block permeability K is defined as CF (K0).
2. offset well quantity credibility CF (W) can rule of thumb delimit degree of membership, as W < 3 sets CF (W)=0.4,3≤
W < 10 then sets CF (W)=0.5, and 10≤W < 20 then sets CF (W)=0.6, and 20≤W < 50 then sets
CF (W)=0.7,50≤W < 200 then set CF (W)=0.8, and 200≤W < 800 then sets CF (W)=0.9, W > 800
Then set CF (W)=1.3. average distance credibility CF (D) can rule of thumb delimit degree of membership, such as D < 100m
Then setting CF (D)=1,100≤D < 2000 then sets CF (D)=0.9, and 2000≤D < 5000 then sets CF (D)=0.8,
5000≤D < 10000 then sets CF (D)=0.7, and 10000≤D < 20000 then sets CF (D)=0.6,20000≤D
< 50000 then sets CF (D)=0.5, and 50000≤D < 100000 then sets CF (D)=0.4 ..., finally calculate three
Geometrical mean (or weighted geometric mean) obtain credibility CF (K) of this point.Different is subordinate to interval
Division can affect the size of CF (K) value, but does not interferes with its relative Changing Pattern, thus without affecting its engineering
Using value.
In step 3, the coordinate array obtained according to step 1, from drilling engineering design scheme or with brill
The engineering parameter of each coordinate points of extracting data gathered, forms engineering parameter array E [N] [I] (I=11), its
Middle i represents drilling fluid density, drilling fluid plastic viscosity, design discharge capacity, design drilling speed, rock respectively from 0 to 10
Bits size, landwaste, borehole diameter, drilling tool internal diameter, drilling tool external diameter, hydrophthalmia diameter, the surge pressure system of setting
Number Sg.Confidence value CF (E) [i] of every kind of engineering parameter is given according to Data Source situation.
If these data come from engineering design, it is taken as that all parameters are the most credible, set each parameter
Confidence value CF (E) [i]=1;If from brill, again can set according to the reliability of collection in worksite data
The value of the CF (E) [i] of some parameter fixed.
It should be noted that for the design or real brill of a bite well, Data Source mode is identical, the most right
There is no difference in same parameter in the credibility of different depth, ignore the most here [N] in array, the most any
CF (E0) [I] this group credibility is all used in depth point.
In step 4, based on the data in step 3, the algorithm that SY/T 6613-2005 standard provides is utilized (also
Other similar algorithms can be used), calculate the pressure equivalent circulating density (i.e. ECD) of each coordinate points, use
rcRepresent, form ECD array rc[N]。
If using industry standard algorithm, then it is assumed that this algorithm is the most credible, the credibility of this algorithm can be set
(i.e. rule intensity) CF (ECD algorithm)=0.98;If, with other algorithms, can be according to the accreditation to this algorithm
The given corresponding CF value of degree, CF ∈ [0,1].
The each engineering parameter used for this algorithm, calculates the credibility of final ECD as follows:
Assume participate in calculate have n parameter E1 ... En, the most first calculate each parameter for ECD calculate only
Vertical credibility CF (Ei)=CF (E) [i] * CF (ECD algorithm), CF (ECD algorithm) represent the algorithm calculating ECD
Credibility, CF (E) [i] represent parameter Ei credibility.Then the synthetic reliability of E1 and E2 is calculated
CF (E1E2)=CF (E1)+CF (E2)-CF (E1) * CF (E2), then uses CF (E1E2) as independent credibility
Similar method and CF (E3) are calculated CF (E1E2E3), and iterative computation is until last parameter bit successively
Put, obtain final summation credibility CF (ECD).Comparing additive method, utilize that said method calculates can
Reliability precision is higher.
In steps of 5, start the prediction of leakage risk, make n=0.
In step 6, according to n value, extract data G [n] from geologic parameter array G [N] [I] respectively, from ECD
Array rc[N] extracts rc[n], extracts array CF (G) [n] from geologic parameter credibility array CF (G) [N] [I],
Array E [n] is extracted from engineering parameter array E [N] [I].
In step 7, if (solution cavity=have) or (tomography=have), then initial risk values R0=5, then base
Credibility (i.e. CF (solution cavity) and CF (tomography)) in solution cavity and tomography calculates credibility CF (R0) of R0,
Final risk value R=R0*CF (R0), then performs step 17.Wherein, (CF (solution cavity), CF are (disconnected for CF (R0)=max
Layer)).
In step 8, if rc+Sg>rf, then R0=5, it is then based on rc、Sg、rfCredibility calculate R0
Credibility CF (R0), final risk value R=R0*CF (R0), then perform step 17.Due to S hereing
For design load, the problem being therefore not related to its credibility, then CF (R0)=CF (rc)+CF(rf)-CF(rc)*
CF(rf).In other examples, if the real going down process bored, then SgSurge pressure equivalent for Practical Calculation
Density, now need to iterative computation r according to step 4c、Sg、rfThe synthetic reliability of these three factor.
In step 9, if (crack=very grow) and (rc+Sg>rp), then R0=5, be then based on crack,
rc、Sg、rfCredibility calculate credibility CF (R0) of R0, final risk value R=R0*CF (R0), so
Rear execution step 17.Here the credibility of R0 can come iterative computation crack, r according to the method for step 4c、rf
The synthetic reliability of these three factor.Under the lower brill state of real drill-through journey, substitute with surge pressure equal yield density
Sg, now considering its credibility, the credibility of R0 can come iterative computation crack, r according to the method for step 4c、
Sg、rfThe synthetic reliability of these four factors.
In step 10, if (crack=typically grow) and (rc+Sg>rp), then R0=4, it is then based on splitting
Seam, rc、Sg、rfCredibility calculate credibility CF (R0) of R0, final risk value R=R0*CF (R0),
Then step 17 is performed.The circular of credibility CF (R0) of R0 can be identical with step 9.
In a step 11, if (rock type=sandstone or conglomerate or glutenite) and (K > 14 μm2) and
(rc+Sg>rp), then R0=4, it is then based on rock type, K, rc、Sg、rfCredibility calculate R0
Credibility CF (R0), final risk value R=R0*CF (R0), then perform step 17.The credibility of R0
The circular of CF (R0) can be similar with step 10.
By this step it can be seen that this example not only allows for stitching the leakage of hole property, but also permeability can be leaked
Mistake is identified.
In step 12, if (rock type=sandstone or conglomerate or glutenite) and (rc+Sg>rp), then R0=3,
It is then based on rock type, rc、Sg、rfCredibility calculate credibility CF (R0) of R0, final risk value
R=R0*CF (R0), then performs step 17.The circular of credibility CF (R0) of R0 can be with step
10 are similar to.
In step 13, if (rock type=sandstone or conglomerate or glutenite) and rc≥rp, then R0=2, so
After based on rock type, rc、rfCredibility calculate credibility CF (R0) of R0, final risk value R=R0*
CF (R0), then performs step 17.Here the credibility of R0 can carry out iterative computation rock according to the method for step 4
Stone type, rc、rfThe synthetic reliability of these three factor.
At step 14, if rc+Sg>rp, then R0=2, it is then based on rc、Sg、rfCredibility calculate
Credibility CF (R0) of R0, final risk value R=R0*CF (R0);Perform step 17.Due to S hereingFor
Design load, the problem being therefore not related to its credibility, then CF (R0)=CF (rc)+CF(rf)-CF(rc)*CF(rf)。
If the real going down process bored, then SgFor the surge pressure equal yield density of Practical Calculation, now need to be according to step 4
Iterative computation r comec、Sg、rfThe synthetic reliability of these three factor.
In step 15, if rc>rp, then R0=1, it is then based on rc、rfCredibility calculate the credible of R0
Degree CF (R0), final risk value R=R0*CF (R0), then perform step 17.Wherein,
CF (R0)=max (CF (rc),CF(rf))。
In step 16, if above-mentioned condition is all unsatisfactory for, then R=0.
In step 17, R [n]=R is made.
In step 18, if n < N, then n=n+1, return and perform step 6~17, otherwise terminate.
Except utilizing above-mentioned alternative manner to calculate credibility CF (R0) of R0, it is also possible to utilize each factor can
The geometrical mean of reliability obtains credibility CF (R0) of R0, in this no limit.
If it should be noted that for brill, engineering parameter is real brill data rather than design data, partly
Matter parameter information also can be updated by measurement and calculating.During real boring, when being in " lower brill " state,
Above-mentioned SgThe surge pressure equal yield density of unification Practical Calculation substitutes.
The present invention is applicable to drilling engineering design personnel, field technician, rear decision-making expert and production management
Personnel carry out the potential leakage risk profile of full well section, and the monitoring risk symptom of a trend in real time during real boring, thus
Reduce construction risk.
Example 1
In the YB*** well Drilling Design stage, according to conventional design cycle, by borehole track, casing programme,
BHA, drill bit, Hydraulics design complete (well depth 0m~6756m).Carry from this design
Take the engineering parameter collection needed for this method, simultaneously from the digital rock mass of this well place YB block, extract this well
The rock mass property value that eye track is passed, i.e. geologic parameter collection.Above-mentioned engineering parameter collection and geologic parameter collection are respectively
With identical depth interval extraction, such as 1m.Then utilize method of the present invention to differentiate and obtain well rail
The risk class of each point on road, finds at 4801m to 5112m well section leakage risk class R > 4, design
Personnel think that this potential risk is the highest, then pass through and carefully study, and repeatedly adjust relevant design parameter, every time
Re-start risk profile, until value-at-risk R < 3, reach expected design, the most finally finalize a text.
This example considers the interaction relationship between concrete drilling engineering design information and geological condition, therefore
Specific aim measure can be taked to be optimized and prevent for concrete well.
Example 2
During SB*** well wellbore construction, utilize the well track information and mud logging technique measured in real time, real
Time obtain engineering parameter collection, and calculate ECD and surge pressure in real time;Digital rock from this well place block
In body, extract the rock mass property value that this drilled wellbore trajectories is passed or passing, utilize simultaneously and measure in real time
Engineering parameter calculate formation pore pressure with revise former numeral rock mass extraction value formation pressure.Then this is utilized
Method described in invention differentiates and obtains well track foremost or drill bit position or the leakage wind of whole well section
Danger rank.One day, in the going down process that makes up a joint, Field Force sets running speed, starts lower brill, logical
Cross this method to judge, show leakage value-at-risk R=4.6 of current bottom hole location (5263.2m), represent risk very
Height, technical staff think may because of under drill-through cause soon caused by pressure oscillation, slow down running speed immediately, renewal
Data differentiate again, and risk class is down to less than 3, effectively prevent malignant event and occur.
The present invention devises a kind of drilling well leakage Forecasting Methodology based on numeral rock mass, extracts from numeral rock mass
Geologic parameter, and drilling engineering parameter coupling analysis, before brill, carry out full well section or the potential well of local well section
Leakage risk profile, in drilling course, the symptom of a trend to leakage risk carries out Real time identification early warning, can aided drilling design
Personnel and construction technical staff control drilling well risk.
The present invention can be used for the Drilling Design stage and full well carry out leakage prediction, and then optimizes design, and proposes
The prediction scheme of reply risk;In work progress, carry out revising the most in real time and early warning to the risk of prediction before drilling,
Can at utmost control the generation of leakage or find risk in the very first time and take effective treatment measures, will loss
It is down to minimum.
Although the embodiment that disclosed herein is as above, but described content only readily appreciates the technology of the present invention
Scheme and the embodiment that uses, be not limited to the present invention.Technology people in any art of the present invention
Member, on the premise of without departing from the spirit and scope that disclosed herein, can be in the form implemented and details
Carry out any amendment and change, but the scope of patent protection of the present invention, still must be with appending claims institute
Define in the range of standard.
One of ordinary skill in the art will appreciate that all or part of step realizing in above-described embodiment method is
Can instruct relevant hardware by program to complete, described program can be stored in an embodied on computer readable
In storage medium, this program upon execution, including the step of method in this embodiment, described storage medium,
As: ROM/RAM, magnetic disc, CD etc..
Claims (10)
1. a drilling well leakage Forecasting Methodology based on numeral rock mass, the method includes:
Gather the borehole track data of current well, to obtain the coordinate at borehole track each position;
According to the coordinate at each position in borehole track, obtain at each position from the digital rock mass of current well
Geologic parameter information;
According to the coordinate at each position in borehole track, obtain the engineering parameter information at each position;
According to the geologic parameter information at each position and the drilling well leakage of this position of engineering parameter acquisition of information
Value-at-risk.
Method the most according to claim 1, it is characterised in that obtaining from the digital rock mass of current well
Take in the step of the geologic parameter information at each position, including:
The numerical value of geologic parameter at each position is extracted from the digital rock mass of current well;
Data Source according to numeral rock mass and building process, the credibility of the geologic parameter at acquisition each position
Value,
Wherein, described geologic parameter include the following: solution cavity information, fault information, crack information, rock type,
Permeability, formation pore pressure equal yield density and formation fracture pressure equal yield density.
Method the most according to claim 2, it is characterised in that obtaining the ginseng of the engineering at each position
In the step of number information, including:
From drilling engineering design scheme or with the number of the engineering parameter at extracting data each position of drilling collection
Value;
Determine the credibility of every kind of engineering parameter according to Data Source, wherein, described engineering parameter includes the following:
Drilling fluid density, drilling fluid plastic viscosity, discharge capacity, drilling speed, landwaste size, landwaste, borehole diameter, drilling tool
Internal diameter, drilling tool external diameter, hydrophthalmia diameter and surge pressure coefficient or surge pressure equal yield density;
The equivalent circulation being calculated at each position based on the engineering parameter information at each position extracted is close
Degree information.
Method the most according to claim 3, it is characterised in that be calculated the pressure at each position
In the step of power equivalent circulating density information, including:
According to the engineering parameter information at each position, select respective algorithms to calculate the pressure at each position and work as
Amount circulating density;
Pressure equivalent is obtained according to calculating pressure equivalent the circulating density algorithm used and the parameter participating in calculating
The credibility of circulating density.
Method the most according to claim 4, it is characterised in that close according to calculating pressure equivalent circulation
The algorithm that degree is used obtains the step of the credibility of pressure equivalent circulating density ECD with the parameter participating in calculating
In, farther include:
That assumes participation calculating has n parameter E1 ... En, calculating each parameter can for the independence that ECD calculates
Reliability CF (Ei)=CF (E) [i] * CF (ECD algorithm), CF (ECD algorithm) represent calculate ECD algorithm can
Reliability, CF (E) [i] represents the credibility of parameter Ei;
Calculate synthetic reliability CF (E1E2)=CF (E1)+CF (E2)-CF (E1) the * CF (E2) of E1 and E2,
Then using CF (E1E2) as independent credibility, carry out being calculated CF (E1E2E3) with CF (E3) by it, depend on
Secondary iterative computation is until last parameter En, and then obtains final synthetic reliability CF (ECD), by it
Credibility as pressure equivalent circulating density ECD.
6. according to the method described in any one of claim 1-5, it is characterised in that at according to each position
In the step of the drilling well leakage value-at-risk of geologic parameter information and this position of engineering parameter acquisition of information, including:
Judge whether to meet pre-conditioned according to the geologic parameter information at each position and engineering parameter information:
If meeting, then determine initial risk values according to pre-conditioned, and calculate the credibility of this initial risk values,
Final risk value is obtained based on initial risk values and credibility thereof;If being unsatisfactory for pre-conditioned, it is determined that final wind
Danger value is 0.
Method the most according to claim 6, it is characterised in that described pre-conditioned include following at least
One of:
There is solution cavity or there is tomography;
Pressure equivalent circulating density breaks more than stratum with surge pressure coefficient or surge pressure equal yield density sum
Split pressure equal yield density;
Crack for grow very much, and pressure equivalent circulating density is close with surge pressure coefficient or surge pressure equivalent
Degree sum is more than formation pore pressure equal yield density;
Crack is general growth, and pressure equivalent circulating density is close with surge pressure coefficient or surge pressure equivalent
Degree sum is more than formation pore pressure equal yield density;
Rock type is sandstone, conglomerate or glutenite, and permeability is more than certain value, pressure equivalent circulating density
It is more than formation pore pressure equal yield density with surge pressure coefficient or surge pressure equal yield density sum;
Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density and surge pressure coefficient or
Surge pressure equal yield density sum is more than formation pore pressure equal yield density;
Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density is more than or equal to formation pore pressure
Power equal yield density;
Pressure equivalent circulating density is more than earth bore with surge pressure coefficient or surge pressure equal yield density sum
Gap pressure equivalent density;
Pressure equivalent circulating density is more than formation pore pressure equal yield density.
Method the most according to claim 6, it is characterised in that
According to met pre-conditioned in the credibility of each parameter calculate the credibility of initial risk values.
Method the most according to claim 1, it is characterised in that also include:
Based on geologic parameter information, engineering parameter information and drilling well leakage value-at-risk, on the basis of same well depth
Leakage risk profile section is shown with curve form.
10. according to the method described in claim 1 or 9, it is characterised in that also include:
When drilling well leakage value-at-risk in a position is more than or equal to setting thresholding,
The engineering parameter information of this position is adjusted, so that the drilling well leakage value-at-risk of this position is less than
Set thresholding.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510163437.1A CN106150484B (en) | 2015-04-08 | 2015-04-08 | A kind of drilling well leakage prediction technique based on digital rock mass |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510163437.1A CN106150484B (en) | 2015-04-08 | 2015-04-08 | A kind of drilling well leakage prediction technique based on digital rock mass |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106150484A true CN106150484A (en) | 2016-11-23 |
CN106150484B CN106150484B (en) | 2019-03-15 |
Family
ID=57336405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510163437.1A Active CN106150484B (en) | 2015-04-08 | 2015-04-08 | A kind of drilling well leakage prediction technique based on digital rock mass |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106150484B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108694459A (en) * | 2017-04-10 | 2018-10-23 | 中国石油化工股份有限公司 | A kind of drilling plan preferred method based on virtual emulation |
CN108804722A (en) * | 2017-04-26 | 2018-11-13 | 中国石油化工股份有限公司 | A kind of calculation method of parameters and device for drilling well emulation |
CN110656932A (en) * | 2019-10-24 | 2020-01-07 | 西南石油大学 | Drilling fluid loss position estimation method in drilling process |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101139925A (en) * | 2006-09-08 | 2008-03-12 | 西南石油大学 | Method for while-drilling testing reservoir parameter property and adjusting well drilling action in real time |
CN101726761A (en) * | 2008-10-15 | 2010-06-09 | 中国石油天然气股份有限公司 | Risk-constrained oil-gas resource spatial distribution prediction method |
CN102128022A (en) * | 2010-12-30 | 2011-07-20 | 中国电子科技集团公司第二十二研究所 | Drilling engineering early warning method and system thereof |
CN202975361U (en) * | 2012-12-10 | 2013-06-05 | 山东大学 | Advance detection apparatus for unfavorable geology and rock mass mechanical property of deep and long tunnel |
CN104271882A (en) * | 2012-03-19 | 2015-01-07 | 哈里伯顿能源服务公司 | Drilling system failure risk analysis method |
-
2015
- 2015-04-08 CN CN201510163437.1A patent/CN106150484B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101139925A (en) * | 2006-09-08 | 2008-03-12 | 西南石油大学 | Method for while-drilling testing reservoir parameter property and adjusting well drilling action in real time |
CN101726761A (en) * | 2008-10-15 | 2010-06-09 | 中国石油天然气股份有限公司 | Risk-constrained oil-gas resource spatial distribution prediction method |
CN102128022A (en) * | 2010-12-30 | 2011-07-20 | 中国电子科技集团公司第二十二研究所 | Drilling engineering early warning method and system thereof |
CN104271882A (en) * | 2012-03-19 | 2015-01-07 | 哈里伯顿能源服务公司 | Drilling system failure risk analysis method |
CN202975361U (en) * | 2012-12-10 | 2013-06-05 | 山东大学 | Advance detection apparatus for unfavorable geology and rock mass mechanical property of deep and long tunnel |
Non-Patent Citations (1)
Title |
---|
张庆华: "油气探井安全评价技术应用研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108694459A (en) * | 2017-04-10 | 2018-10-23 | 中国石油化工股份有限公司 | A kind of drilling plan preferred method based on virtual emulation |
CN108694459B (en) * | 2017-04-10 | 2021-11-02 | 中国石油化工股份有限公司 | Drilling scheme optimization method based on virtual simulation |
CN108804722A (en) * | 2017-04-26 | 2018-11-13 | 中国石油化工股份有限公司 | A kind of calculation method of parameters and device for drilling well emulation |
CN108804722B (en) * | 2017-04-26 | 2021-07-27 | 中国石油化工股份有限公司 | Parameter calculation method and device for drilling simulation |
CN110656932A (en) * | 2019-10-24 | 2020-01-07 | 西南石油大学 | Drilling fluid loss position estimation method in drilling process |
CN110656932B (en) * | 2019-10-24 | 2020-04-24 | 西南石油大学 | Drilling fluid loss position estimation method in drilling process |
Also Published As
Publication number | Publication date |
---|---|
CN106150484B (en) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106150485B (en) | A kind of drilling well leakage forecasting system based on digital rock mass | |
CN104730596B (en) | A kind of discrete fractures modeling method based on multiple dimensioned factor constraints | |
CN104790926B (en) | A kind of fracture-pore reservoir waterflooding development effect evaluation method | |
US10495782B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
EP2948884B1 (en) | Hazard avoidance analysis | |
US10168447B2 (en) | Automatic geosteering and evolutionary algorithm for use with same | |
CN105484724A (en) | Drilling downhole anomaly monitoring method | |
CN106150484A (en) | A kind of drilling well leakage Forecasting Methodology based on numeral rock mass | |
CN110646847B (en) | Method for identifying low-order fault breakpoint position based on well logging data of close well pattern area | |
CN105484725A (en) | Drilling downhole anomaly monitoring device | |
CN108952657A (en) | Horizontal well platform fracturing fracture length determination method and device | |
CN107153881B (en) | Drilling risk prediction method | |
CN106150476A (en) | A kind of system of the viscous suction bit freezing risk predicting drill string | |
Singh et al. | A comprehensive review of fracture-driven interaction in unconventional oil and gas plays: Characterization, real-time diagnosis, and impact on production | |
Cai et al. | Developing a geomechanics-modeling based method for lost circulation risk assessment: A case study in Bohai Bay, China | |
Jelušič et al. | An adaptive network fuzzy inference system approach for site investigation | |
CN106156385A (en) | A kind of method of the viscous suction bit freezing risk predicting drill string | |
Yi et al. | Application of Formation Pressure While Drilling Prediction Technology in Offshore HTHP Wells | |
Tyree | Predicting Frac Hits Using Artificial Intelligence; Application in Marcellus Shale | |
CN114033352B (en) | Method and equipment for estimating density of well Zhou Liefeng | |
Artun et al. | A pattern-based approach to waterflood performance prediction using knowledge management tools and classical reservoir engineering forecasting methods | |
CN114169558A (en) | Method and system for predicting well kick risk | |
Sitinjak | SURFACE PRODUCTION NETWORK OPTIMIZATION BY MAXIMIZING WELL PERFORMANCE TO IMPROVE GAS RECOVERY IN MULTILAYER MARGINAL GAS FIELD | |
Carpenter | Automated Reservoir Model Calibration Applied to a Complex Multizone Heavy Oil Field | |
CN117189042A (en) | Marine thickened oil thermochemical oil extraction method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |